Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structural variant”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Development and Evaluation of Embrittlement Resistant Alloys for Advanced LWR Cladding

The primary challenge preventing the deployment of reactor designs that leverage wrought FeCrAl as an advanced light-water reactor (LWR) cladding is irradiation hardening and embrittlement. Wrought FeCrAl alloys experience a loss of post-irradiation ductility and fracture toughness under low-temperature neutron irradiation (< 350°C) resulting from the combined effects of dislocation loop formation and the irradiation-enhanced precipitation of Cr-rich alpha-prime precipitates throughout the microstructure. Over the past decade, significant improvements in our understanding regarding the effect of Cr and Al content have been enabled through neutron irradiations, but even optimized wrought alloys such as C26M remain vulnerable to brittle failure during storage, transportation, and handling following irradiation. This report summarizes multiple irradiation campaigns initiated over the past 7 years and provides recent insights into the effect of (1) minor alloying composition, (2) alloy processing strategy, (3) crystal structure, and (4) the use of advanced oxide dispersion strengthened (ODS) alloys. The first irradiation campaign irradiated body-centered cubic FeCrAl alloy C26M with and without elements such as Mo and Y. It also included face-centered cubic alumina-forming austenitic (AFA) alloys that have a steeper neutronic penalty but that have far superior high-temperature strength than wrought FeCrAl. Finally, this campaign compared conventional wrought C26M with the same material produced using modern powder metallurgical hotisostatic pressing (PM-HIP). The second irradiation campaign, initiated in 2018, irradiated wrought and FeCrAl-ODS materials up to 50 dpa to assess their irradiation resistance to high-doses. The results of this work indicate that wrought C26M suffers significant degradation following neutron irradiation at LWR-relevant temperatures regardless of Mo and Y content. The PM-HIP variant appears to be more resilient to irradiation-induced ductility loss in comparison with wrought variants but still loses almost all ductility by 8 dpa. Notably, the AFAs retained at least 5% total elongation after irradiation at the same dose and temperature condition, although additional deterioration may be expected at higher fluence levels. FeCrAl-ODS materials show the greatest promise as a transformative longer-term accident-tolerant fuel (LT-ATF) cladding material. The extruded 106ZY10C alloy retained 10% total elongation after 16 dpa irradiation and retained 8% total elongation after 50 dpa. Additional challenges with FeCrAl-ODS alloys also remain, including optimization of end cap joining methodologies, enhancement of fracture toughness, scaling production to prove economic viability, and optimizing post-pilger heat treatments to maximize ductility and irradiation resistance. This series of irradiations demonstrates the efficacy of the LWR testbed available within the United States for the rapid irradiation and down-selection of LT-ATF candidates. As the country rapidly accelerates its timeline for the deployment of advanced reactor concepts, the effective utilization of the High Flux Isotope Reactor for separate-effects style irradiations should continue to be prioritized to answer the final questions pertaining to LT-ATF candidates necessary for the deployment of advanced boiling water reactors and small modular reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Prospects for additive manufacturing of nuclear fuel forms

Additive manufacturing (AM) is viewed as a critical enabling technology for achieving superior performance and improved economics across many sectors. Although the past decade has seen increased application of AM methods to nuclear reactor cladding and structural materials, exploration of AM as applied to the uranium-bearing nuclear fuel forms has been limited. The two major families of nuclear fuel forms (monolithic, particle/dispersion) are utilized with their own set of core objectives in mind, and their differences require distinct fabrication infrastructures. These reference fabrication methods impose many limitations on nuclear fuels. Further, both currently operating reactors and future concepts have the potential for improved performance if these accepted limitations are relaxed or removed entirely. The primary limitations of reference fabrication processes for the common nuclear fuel forms are outlined in this paper. This groundwork is then used to identify avenues of fuel performance specific to each of these fuel architectures that could be exploited if the restrictions of conventional fuel fabrication are removed. Moreover, multiple targets for AM studies are laid out for each of the major nuclear fuel variants. Finally, key strategic components to guide research activities in AM of nuclear fuels are outlined, with an emphasis on use of modeling and simulation to motivate research aims and embrace of an accelerated testing methodology to screen and quality new fuel forms.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Timely Reporting of Heavy Hitters Using External Memory

Given an input stream S of size N, a Φ-heavy hitter is an item that occurs at least ΦN times in S. The problem of finding heavy-hitters is extensively studied in the database literature. In this work, we study a real-time heavy-hitters variant in which an element must be reported shortly after we see its T = Φ N-th occurrence (and hence it becomes a heavy hitter). We call this the Timely Event Detection (TED) Problem. The TED problem models the needs of many real-world monitoring systems, which demand accurate (i.e., no false negatives) and timely reporting of all events from large, high-speed streams with a low reporting threshold (high sensitivity). Like the classic heavy-hitters problem, solving the TED problem without false-positives requires large space (Ω (N) words). Thus in-RAM heavy-hitters algorithms typically sacrifice accuracy (i.e., allow false positives), sensitivity, or timeliness (i.e., use multiple passes). We show how to adapt heavy-hitters algorithms to external memory to solve the TED problem on large high-speed streams while guaranteeing accuracy, sensitivity, and timeliness. Our data structures are limited only by I/O-bandwidth (not latency) and support a tunable tradeoff between reporting delay and I/O overhead. With a small bounded reporting delay, our algorithms incur only a logarithmic I/O overhead. We implement and validate our data structures empirically using the Firehose streaming benchmark. Multi-threaded versions of our structures can scale to process 11M observations per second before becoming CPU bound. In comparison, a naive adaptation of the standard heavy-hitters algorithm to external memory would be limited by the storage device’s random I/O throughput, i.e., ≈100K observations per second.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Mapping the influence of ligand electronics on the spectroscopic and 1 O 2 sensitization characteristics of Pd( ii ) biladiene complexes bearing phenyl–alkynyl groups at the 2- and 18-positions

Photodynamic therapy (PDT) is a promising treatment for certain cancers that proceeds via sensitization of ground state 3 O 2 to generate reactive 1 O 2 . Classic macrocyclic tetrapyrrole ligand scaffolds, such as porphyrins and phthalocyanines, have been studied in detail for their 1 O 2 photosensitization capabilities. Despite their compelling photophysics, these systems have been limited in PDT applications because of adverse biological side effects. Conversely, the development of non-traditional oligotetrapyrrole ligands metalated with palladium (Pd[DMBil1]) have established new candidates for PDT that display excellent biocompatibility. Herein, the synthesis, electrochemical, and photophysical characterization of a new family of 2,18-bis(phenylalkynyl)-substituted Pd II 10,10-dimethyl-5,15-bis(pentafluorophenyl)-biladiene (Pd[DMBil2-R]) complexes is presented. These second generation biladienes feature extended conjugation relative to previously characterized Pd II biladiene scaffolds (Pd[DMBil1]). We show that these new derivatives can be prepared in good yield and, that the electronic nature of the phenylalkynyl appendages dramatically influence the Pd II biladiene photophysics. Extending the conjugation of the Pd[DMBil1] core through installation of phenylacetylene resulted in a ~75 nm red-shift of the biladiene absorption spectrum into the phototherapeutic window (600–900 nm), while maintaining the Pd II biladiene's steady-state spectroscopic 1 O 2 sensitization characteristics. Varying the electronics of the phenylalkyne groups via installation of electron donating or withdrawing groups dramatically influences the steady-state spectroscopic and photophysical properties of the resulting Pd[DMBil2-R] family of complexes. The most electron rich variants (Pd[DMBil2-N(CH3)2]) can absorb light as far red as ~700 nm but suffer from significantly reduced ability to sensitize formation of 1 O 2 . By contrast, Pd[DMBil2-R] derivatives bearing electron withdrawing functionalities (Pd[DMBil2-CN] and Pd[DMBil2-CF3]) display 1 O 2 quantum yields above 90%. The collection of results we report suggest that excited state charge transfer from more electron-rich phenyl-alkyne appendages to the electron deficient biladiene core circumvents triplet sensitization. The spectral and redox properties, as well as the triplet sensitization efficiency of each Pd[DMBil2-R] derivative is considered in relation to the Hammett value (σp) for each biladiene's R-group. More broadly, the results reported in this study clearly demonstrate that biladiene redox properties, spectral properties, and photophysics can be perturbed greatly by relatively minor alterations to biladiene structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Dissecting neurofilament tail sequence-phosphorylation-structure relationships with multicomponent reconstituted protein brushes

Neurofilaments (NFs) are multisubunit, bottlebrush-shaped intermediate filaments abundant in the axonal cytoskeleton. Each NF subunit contains a long intrinsically disordered tail domain, which protrudes from the NF core to form a “brush” surrounding each NF. Precisely how the tails’ variable charge patterns and repetitive phosphorylation sites mediate their conformation within the brush remains an open question in axonal biology. We address this problem by grafting recombinant NF tail protein constructs NF-Light, -Medium, and -Heavy (NFL, NFM, and NFH) to surfaces, yielding protein brushes of defined stoichiometry that can be phosphorylated in vitro. Atomic force microscopy measurements reveal that brush height depends on composition monotonically but not always linearly for binary NFL:NFM or NFL:NFH systems, and that NFM-based brushes are highly extended, while brushes incorporating the much larger NFH are surprisingly compact even after multisite phosphorylation. Complementary self-consistent field theory (SCFT) predicts multilayer brush morphologies for NFM and phosphorylated NFH brushes. Further experiments and SCFT analysis with designed mutants reveal that N-terminal negative charges in the NFH tail repel phosphorylated residues to generate the multilayer morphology, while the C-terminal charge-neutral region contributes to multilayer brush morphology but not total brush height. Charge-shuffled NFM variants show that charge segregation promotes brush collapse near physiological ionic strengths. Collectively, this study supports a role for NFM in establishing a dynamic range for NF brush conformation, lending insight into previous in vitro and in vivo findings. More broadly, this work establishes a platform for dissecting contributions of disordered protein sequence to conformation at interfaces.

Science & Technology - Other Topics↗
DOE OSTI2021

Peptidomimetic nitrile warheads as SARS-CoV-2 3CL protease inhibitors

Tragically, the death toll from the COVID-19 pandemic continues to rise, and with variants being observed around the globe new therapeutics, particularly direct-acting antivirals that are easily administered, are desperately needed. Studies targeting the SARS-CoV-2 3CL protease, which is critical for viral replication, with different peptidomimetics and warheads is an active area of research for development of potential drugs. To date, however, only a few publications have evaluated the nitrile warhead as a viral 3CL protease inhibitor, with only modest activity reported. This article describes our investigation of P3 4-methoxyindole peptidomimetic analogs with select P1 and P2 groups with a nitrile warhead that are potent inhibitors of SARS-CoV-2 3CL protease and demonstrate in vitro SARS-CoV-2 antiviral activity. A selectivity for SARS-CoV-2 3CL protease over human cathepsins B, S and L was also observed with the nitrile warhead, which was superior to that with the aldehyde warhead. A co-crystal structure with SARS-CoV-2 3CL protease and a reversibility study indicate that a reversible, thioimidate adduct is formed when the catalytic sulfur forms a covalent bond with the carbon of the nitrile. Furthermore, this effort also identified efflux as a property limiting antiviral activity of these compounds, and together with the positive attributes described these results provide insight for further drug development of novel nitrile peptidomimetics targeting SARS-CoV-2 3CL protease.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Multiple Critical Unresolved Region Integral Experiment (MCURIE): A Proposed Integral Critical Experimental Framework for Unresolved Region and Intermediate Energies

The Critical Unresolved Region Integral Experiment (CURIE) critical experiment was performed at the National Criticality Experiments Research Center (NCERC) at the Device Assembly Facility (DAF) at the Nevada Nuclear Security Sites (NNSS) in 2020. The objective of CURIE was to improve the quality of integral nuclear data in the uranium-235 ( 235 U) unresolved resonance region (URR) by performing benchmark integral experiments that were sensitive to the URR energy ranges. The CURIE experiment was evaluated for the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook and the benchmark evaluation was accepted in 2022. The URR is a region within the intermediate neutron energy range (the intermediate energy ranges from 0.7 eV to 100 keV). The observed resonance structure in neutron cross sections is due to discrete energy levels in the nucleus and are characterized by resonance parameters. In the URR region, the resonance parameters are only partially resolved as the resolution of the experiment techniques becomes comparable to the average width of the resonances themselves, and the resonances are so close to one another that the structure cannot be determined empirically. The ENDF/B-VIII.0 and JEFF-3.3 nuclear data libraries define the URR as beginning at 2.25 keV and continuing until 25 keV. There are minimal intermediate neutron energy benchmarks available in the ICSBEP benchmark handbook, and aside from CURIE there are none that are highly sensitive in the URR energy region. There is a current need for additional experiments sensitive to the URR energy region. There is a new proposed subgroup for the Organization for Economic Co-operation and Development and Nuclear Energy Agency (OECD NEA) Working Party on International Nuclear Data Evaluation Cooperation (WPEC), so any new models or information will need experimental validation and testing. The NEA Working Party on Nuclear Criticality Safety (WPNCS) recent experimental needs and priority list includes intermediate energy 235 U and 238 U experiments, as does the recent Integral Experiments to Address Nuclear Criticality Safety Needs Meeting (May 2023). A variant on the original CURIE experiment called the Multiple Critical Unresolved Region Integral Experiment (MCURIE) is proposed to provide further investigation of the intermediate and URR energy region for uranium. MCURIE will utilize existing fuels at NCERC but use alternative moderators and reflectors to modify the neutron absorption, scattering, and fission spectra of the experiments, allowing for precise targeting of nuclear data sensitivities. The overall goal of MCURIE is to develop a framework for designing and performing integral benchmark experiments with high sensitivities in the intermediate and URR energy regions for nuclear data validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Variant Nelson-Barr mechanism with minimal flavor violation

Within the general framework of using spontaneous CP violation to solve the strong CP problem, we construct a variant Nelson-Barr model in which the Standard Model (SM) quark contribution to the strong CP phase is cancelled by new heavy QCD-charged fermions. This cancellation is ensured by choosing conjugate representations for the new colored states under the same global flavor symmetry of SM quarks. Choosing the global flavor symmetry to be that of minimal flavor violation, we suppress higher-order corrections to the strong CP phase to well below current experimental constraints. More than two dozen massless Goldstone bosons emerge from spontaneous flavor symmetry breaking, which yield strong astrophysical constraints on the symmetry breaking scale. In the early universe, the Goldstone bosons can be thermally produced from their interactions with the heavy colored fermions and contribute to ΔN eff at a measurable level. As a function of reheating temperature, the predicted ΔN eff shows an interesting plateau behavior we dub the “flavor stairway”, which encodes information about the SM quark flavor structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Noonan syndrome patient-specific induced cardiomyocyte model carrying SOS1 gene variant c.1654A>G

Noonan syndrome (NS) is a dominant autosomal genetic disorder, associated with mutations in several genes that exhibit multisystem abnormal development including cardiac defects. NS associated with the Son of Sevenless homolog 1 (SOS1) gene mutation attributes to the development of cardiomyopathy and congenital heart defects. Since the treatment option for NS is very limited, an in vitro disease model with SOS1 gene mutation would be beneficial for exploring therapeutic possibilities for NS. We reprogrammed cardiac fibroblasts obtained from a NS patient and normal control skin fibroblasts (C-SF) into induced pluripotent stem cells (iPSCs). We identified NS-iPSCs carry a heterozygous single nucleotide variation in the SOS1 gene at the c.1654A > G. Furthermore, the control and NS-iPSCs were differentiated into induced cardiomyocytes (iCMCs), and the electron microscopic analysis showed that the sarcomeres of the NS-iCMCs were highly disorganized. FACS analysis showed that 47.5% of the NS-iCMCs co-expressed GATA4 and cardiac troponin T proteins, and the mRNA expression levels of many cardiac related genes, studied by qRT-PCR array, were significantly reduced when compared to the control C-iCMCs. We report for the first time that NS-iPSCs carry a single nucleotide variation in the SOS1 gene at the c.1654A>G were showing significantly reduced cardiac genes and proteins expression as well as structurally and functionally compromised when compared to C-iCMCs. These iPSCs and iCMCs can be used as a modeling platform to unravel the pathologic mechanisms and also the development of novel drug for the cardiomyopathy in patients with NS.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2024

Polymorphism in the Ruddlesden–Popper Nickelate La 3 Ni 2 O 7 : Discovery of a Hidden Phase with Distinctive Layer Stacking

We report the discovery of a novel form of Ruddlesden–Popper (RP) nickelate that stands as the first example of long-range, coherent polymorphism in this class of inorganic solids. Rather than the well-known, uniform stacking of perovskite blocks ubiquitously found in RP phases, this newly discovered polymorph of the bilayer RP phase La 3 Ni 2 O 7 adopts a novel stacking sequence in which single-layer and trilayer blocks of NiO 6 octahedra alternate in a “1313” sequence. Crystals of this new polymorph are described in space group Cmmm, although we note evidence for a competing Imam variant. Transport measurements at ambient pressure reveal metallic character with evidence of a charge density wave transition with an onset at T ≈ 134 K. In conclusion, the discovery of such polymorphism could reverberate to the expansive range of science and applications that rely on RP materials, particularly the recently reported signatures of superconductivity in bilayer La 3 Ni 2 O 7 with T c as high as 80 K above 14 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Open-Shell Tensor Hypercontraction

The extension of least-squares tensor hypercontracted second- and third-order Møller–Plesset perturbation theory (LS-THC-MP2 and LS-THC-MP3) to open-shell systems is an important development due to the scaling reduction afforded by THC and the ubiquity of molecular ions, radicals, and other open-shell reactive species. The complexity of wavefunction-based quantum-chemical methods such as Møller–Plesset and coupled cluster theory is reflected in the steep scaling of the computational costs with the molecular size. The least-squares tensor hypercontraction (LS-THC) method is an efficient, single-step factorization for the two-electron integral tensor but can also be used to factorize the double excitation amplitudes, leading to significant scaling reduction. Herein we extend this promising method to open-shell variants of LS-THC-MP2 and -MP3 by using diagrammatic techniques and explicit spin summation. The accuracy of the resulting methods for open-shell species is benchmarked on standard test systems such as regular alkanes as well as realistic systems involving bond breaking, radical stabilization, and other effects. We find that open-shell LS-THC-MPn methods exhibit errors highly comparable to those produced by closed-shell LS-THC-MPn and are highly insensitive to particular chemical interactions, geometries, or even moderate spin contamination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Self-assembling Shell Proteins PduA and PduJ have Essential and Redundant Roles in Bacterial Microcompartment Assembly

Protein self-assembly is a common and essential biological phenomenon, and bacterial microcompartments present a promising model system to study this process. Bacterial microcompartments are large, protein-based organelles which natively carry out processes important for carbon fixation in cyanobacteria and the survival of enteric bacteria. These structures are increasingly popular with biological engineers due to their potential utility as nanobioreactors or drug delivery vehicles. However, the limited understanding of the assembly mechanism of these bacterial microcompartments hinders efforts to repurpose them for non-native functions. Here, we comprehensively investigate proteins involved in the assembly of the 1,2-propanediol utilization bacterial microcompartment from Salmonella enterica serovar Typhimurium LT2, one of the most widely studied microcompartment systems. We first demonstrate that two shell proteins, PduA and PduJ, have a high propensity for self-assembly upon overexpression, and we provide a novel method for self-assembly quantification. Using genomic knock-outs and knock-ins, we systematically show that these two proteins play an essential and redundant role in bacterial microcompartment assembly that cannot be compensated by other shell proteins. At least one of the two proteins PduA and PduJ must be present for the bacterial microcompartment shell to assemble. We also demonstrate that assembly-deficient variants of these proteins are unable to rescue microcompartment formation, highlighting the importance of this assembly property. Our work provides insight into the assembly mechanism of these bacterial organelles and will aid downstream engineering efforts.

1,2-propanediol utilization MCP↗
DOE OSTI2024

Decomposition and Algorithmic Approaches for Solving Large-Scale Process Family Design Problems

Our most recent work expands the water desalination case study from 76 variants to 10,897 variants using the equation-oriented model built in Pyomo as part of the PARETO project. Using the discretization formulation presented in Stinchfield (2024a), rather than solving for all 10,897 variants simultaneously, we decompose the formulation into subproblems containing subsets of variants from the process family. We solve the overall problem with Progressive Hedging (PH) deployed in parallel on a distributed HPC cluster using the open-source Python package mpi-sppy (Knueven et al., 2023). This approach allowed us to solve this process family design problem to ~1.5% relative optimality gap in about 5 hours; in comparison, Gurobi reached ~50% relative optimality gap in about 6 hours (Stinchfield et al., 2024b). However, this approach still requires discretization of the common unit module design ranges; additionally, PH acts as a heuristic for MILP’s with gap-closing capabilities. Ideally, we would not have to use ML surrogates or discretization to solve this problem, instead solving the process family design problem with the equation-oriented model directly to achieve the most accurate results. However, recall that we did not consider solving the MINLP directly due to complexity and size. In this work, we aim to decompose and solve this large-scale MINLP using a Structured Nonlinear Global Optimization algorithm presented by Cao and Zavala (2019).

Stinchfield, Georgia↗
DOE OSTI2022

Single- and two-particle finite size effects in interacting lattice systems

Simulations of extended quantum systems are typically performed by extrapolating results of a sequence of finite-system-size simulations to the thermodynamic limit. In the quantum Monte Carlo community, twist-averaging was pioneered as an efficient strategy to eliminate one-body finite size effects. In the dynamical mean field community, cluster generalizations of the dynamical mean field theory were formulated to study systems with nonlocal correlations. In this work, we put the twist-averaging and the dynamical cluster approximation variant of the dynamical mean field theory onto equal footing, discuss commonalities and differences, and compare results from both techniques to the standard periodic boundary technique. At the example of Hubbard-type models with local, short-range and Yukawa-like longer range interactions we show that all methods converge to the same limit, but that the convergence speed differs in practice. We show that embedding theories are an effective tool for managing both one-body and two-body finite size effects, in particular if interactions are averaged over twist angles.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Analysis of turbulent underexpanded jets. II - Shock noise features using SCIPVIS

SCIPVIS, the computational model discussed by Dash et al. (1985), is assessed in predicting the complicated flow structure associated with shock-containing plumes. In addition, the analysis in this study examines this code's applicability as a basic part of a program for estimating broadband shock noise radiation. The results of this study show that excellent agreement exists between predicted and measured static pressure distributions for both underexpanded and overexpanded flow cases considered. Of the three turbulence closure models incorporated in the SCIPVIS code, the kW model of Spalding produces the most uniform agreement with measurement. The k-epsilon-2 model of Launder consistently overestimates plume spreading for supersonic jets with exit Mach numbers in the 1-2 range. Dash's (1983) k-epsilon-2-cc, compressibility-corrected version of Launder's model underestimates plume spreading. Good qualitative agreement was also obtained between the measured longitudinal turbulence intensity and that predicted by the code for the same trial case. Comparison of measured and predicted broadband shock noise spectrum peak values were found to be in excellent agreement. This utilized a variant of the Harper-Bourne and Fisher (1973) phase-array model: the effective shock spacing was reinterpreted as the value of the end of the plume potential core, determined herein by the SCIPVIS code.

Seiner, J. M.↗
DOE OSTI2024

Dual engagement of the nucleosomal acidic patches is essential for deposition of histone H2A.Z by SWR1C

The yeast SWR1C chromatin remodeling enzyme catalyzes the ATP-dependent exchange of nucleosomal histone H2A for the histone variant H2A.Z, a key variant involved in a multitude of nuclear functions. How the 14-subunit SWR1C engages the nucleosomal substrate remains largely unknown. Studies on the ISWI, CHD1, and SWI/SNF families of chromatin remodeling enzymes have demonstrated key roles for the nucleosomal acidic patch for remodeling activity, however a role for this nucleosomal epitope in nucleosome editing by SWR1C has not been tested. Here, we employ a variety of biochemical assays to demonstrate an essential role for the acidic patch in the H2A.Z exchange reaction. Utilizing asymmetrically assembled nucleosomes, we demonstrate that the acidic patches on each face of the nucleosome are required for SWR1C-mediated dimer exchange, suggesting SWR1C engages the nucleosome in a ‘pincer-like’ conformation, engaging both patches simultaneously. Loss of a single acidic patch results in loss of high affinity nucleosome binding and nucleosomal stimulation of ATPase activity. We identify a conserved arginine-rich motif within the Swc5 subunit that binds the acidic patch and is key for dimer exchange activity. In addition, our cryoEM structure of a Swc5–nucleosome complex suggests that promoter proximal, histone H2B ubiquitylation may regulate H2A.Z deposition. Together these findings provide new insights into how SWR1C engages its nucleosomal substrate to promote efficient H2A.Z deposition.

59 BASIC BIOLOGICAL SCIENCES↗