Search NASA⌕ Search

SEARCH · Search NASA

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.

663 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Neutralizing Antibodies to SARS-CoV-2 Selected from a Human Antibody Library Constructed Decades Ago

Combinatorial antibody libraries not only effectively reduce antibody discovery to a numbers game, but enable documentation of the history of antibody responses in an individual. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has prompted a wider application of this technology to meet the public health challenge of pandemic threats in the modern era. Herein, a combinatorial human antibody library constructed 20 years before the coronavirus disease 2019 (COVID-19) pandemic is used to discover three highly potent antibodies that selectively bind SARS-CoV-2 spike protein and neutralize authentic SARS-CoV-2 virus. Compared to neutralizing antibodies from COVID-19 patients with generally low somatic hypermutation (SHM), these three antibodies contain over 13–22 SHMs, many of which are involved in specific interactions in their crystal structures with SARS-CoV-2 spike receptor binding domain. The identification of these somatically mutated antibodies in a pre-pandemic library raises intriguing questions about the origin and evolution of these antibodies with respect to their reactivity with SARS-CoV-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Overview of Diverse Methyl/Alkyl-Coenzyme M Reductases and Considerations for Their Potential Heterologous Expression

Methyl-coenzyme M reductase (MCR) is an archaeal enzyme that catalyzes the final step of methanogenesis and the first step in the anaerobic oxidation of methane, the energy metabolisms of methanogens and anaerobic methanotrophs (ANME), respectively. Variants of MCR, known as alkyl-coenzyme M reductases, are involved in the anaerobic oxidation of short-chain alkanes including ethane, propane, and butane as well as the catabolism of long-chain alkanes from oil reservoirs. MCR is a dimer of heterotrimers (encoded by mcrABG) and requires the nickel-containing tetrapyrrole prosthetic group known as coenzyme F430. MCR houses a series of unusual post-translational modifications within its active site whose identities vary depending on the organism and whose functions remain unclear. Methanogenic MCRs are encoded in a highly conserved mcrBDCGA gene cluster, which encodes two accessory proteins, McrD and McrC, that are believed to be involved in the assembly and activation of MCR, respectively. The requirement of a unique and complex coenzyme, various unusual post-translational modifications, and many remaining questions surrounding assembly and activation of MCR largely limit in vitro experiments to native enzymes with recombinant methods only recently appearing. Production of MCRs in a heterologous host is an important step toward developing optimized biocatalytic systems for methane production as well as for bioconversion of methane and other alkanes into value-added compounds. This review will first summarize MCR catalysis and structure, followed by a discussion of advances and challenges related to the production of diverse MCRs in a heterologous host.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Community Structure of Arbuscular Mycorrhizal Fungi in Soils of Switchgrass Harvested for Bioenergy

We assessed the different species of beneficial fungi living in agricultural fields of switchgrass, a large grass grown for biofuels, using high-resolution DNA sequencing. Contrary to our expectations, the fungi were not greatly affected by fertilization. However, we found a positive relationship between plant productivity and the number of families of beneficial fungi at one site. Furthermore, we sequenced many species that could not be identified with existing reference databases. One group of fungi was highlighted in an earlier study for being widely distributed but of unknown taxonomy. We discovered that this group belonged to a family called Pervetustaceae , which may benefit switchgrass in stressful environments. To produce higher-yielding switchgrass in a more sustainable manner, it could help to study these undescribed fungi and the ways in which they may contribute to greater switchgrass yield in the absence of fertilization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Fock-Space Schrieffer–Wolff Transformation: Classically-Assisted Rank-Reduced Quantum Phase Estimation Algorithm

We present an extension of many-body downfolding methods to reduce the resources required in the quantum phase estimation (QPE) algorithm. In this paper, we focus on the Schrieffer–Wolff (SW) transformation of the electronic Hamiltonians for molecular systems that provides significant simplifications of quantum circuits for simulations of quantum dynamics. We demonstrate that by employing Fock-space variants of the SW transformation (or rank-reducing similarity transformations (RRST)) one can significantly increase the locality of the qubit-mapped similarity-transformed Hamiltonians. The practical utilization of the SW-RRST formalism is associated with a series of approximations discussed in the manuscript. In particular, amplitudes that define RRST can be evaluated using conventional computers and then encoded on quantum computers. The SW-RRST QPE quantum algorithms can also be viewed as an extension of the standard state-specific coupled-cluster downfolding methods to provide a robust alternative to the traditional QPE algorithms to identify the ground and excited states for systems with various numbers of electrons using the same Fock-space representations of the downfolded Hamiltonian. The RRST formalism serves as a design principle for developing new classes of approximate schemes that reduce the complexity of quantum circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

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

Two-dimensional heteronuclear single quantum coherence (HSQC) NMR spectra of lignin isolated from Populus trichocarpa residues after CELF pretreatment and CBP fermentation

Here we present a curated dataset of a series of two-dimensional heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectra of lignin isolated from a woody energy crop (Populus trichocarpa) residues after co-solvent enhanced lignocellulosic fractionation (CELF) pretreatment and consolidated bioprocessing (CBP) process. The natural poplar variant GW-9947 from the Center for Bioenergy Innovation (CBI) was used. The poplar was knife milled and passed through a 1 mm sieve. The CELF pretreatment was performed in a Parr autoclave reactor with 7.5 wt % solids loading, 0.5 wt% H2SO4 as catalyst at 150°C with 15, 25 and 30 minutes, respectively. Tetrahydrofuran was added in a 1:1 mass ratio with water as the pretreatment solvent. The residues from CELF pretreatment were then subjected to CBP using the bacterium C. thermocellum DSM 1313. CBP fermentations were performed at 60 °C in a shaker at 50 grams/L solids loadings. Lignin was isolated from the pretreated samples after ball-milling in a porcelain jar with ceramic balls via Retsch PM 200 at 580 rpm for 2.5 h followed by enzymatic hydrolysis in acetate buffer (pH 4.8, 50 °C) for 48 h. The lignin samples were characterized using 13C–1H HSQC experiments which were performed in a Bruker Avance III HD 500 MHz NMR spectrometer operating at a frequency of 125.12 MHz for the 13C nucleus. A standard Bruker pulse sequence was used on a Prodigy platform cryoprobe. The dry lignin samples were dissolved in deuterated dimethylsulfoxide for HSQC experiments. The spectra were acquired under the following acquisition conditions: 210 ppm spectral width in F1 (13C) dimension with 256 data points and 11 ppm spectral width in F2 (1H) dimension with 1024 data points, a 90° pulse, a one bond C–H coupling constant of 145 Hz, a 1.0 s pulse delay, and 64 scans. All the data was processed using the TopSpin 3.6 software (Bruker BioSpin). The NMR spectra provides structural characteristics information about lignin remaining in solids after CELF (150 °C with 15, 25 and 30 minutes) process and C. thermocellum CBP.

Lignin structure, HSQC, poplar, CELF, CBP, CBI↗
DOE OSTI2023

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low-rank approximation of the dielectric matrix

We present an efficient preconditioning technique for accelerating the fixed-point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low-rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed-point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous self-consistent field iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from ~100 to 1100 atoms (~500–20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized d orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20–30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed-point iteration in some cases. Lastly, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4-fold reduction in computational cost for the total ground-state calculation compared to other preconditioners.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Hepatitis C virus NS3/4A inhibitors and other drug-like compounds as covalent binders of SARS-CoV-2 main protease

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), threatens global public health. The world needs rapid development of new antivirals and vaccines to control the current pandemic and to control the spread of the variants. Among the proteins synthesized by the SARS-CoV-2 genome, main protease (M pro also known as 3CL pro ) is a primary drug target, due to its essential role in maturation of the viral polyproteins. In this study, we provide crystallographic evidence, along with some binding assay data, that three clinically approved anti hepatitis C virus drugs and two other drug-like compounds covalently bind to the M pro Cys145 catalytic residue in the active site. Also, molecular docking studies can provide additional insight for the design of new antiviral inhibitors for SARS-CoV-2 using these drugs as lead compounds. One might consider derivatives of these lead compounds with higher affinity to the M pro as potential COVID-19 therapeutics for further testing and possibly clinical trials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ex vivo to in vivo model of malignant peripheral nerve sheath tumors for precision oncology

Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas that often develop in patients with neurofibromatosis type 1 (NF1). To address the critical need for novel therapeutics in MPNST, we aimed to establish an ex vivo 3D platform that accurately captured the genomic diversity of MPNST and could be utilized in a medium-throughput manner for drug screening studies to be validated in vivo using patient-derived xenografts (PDX). Genomic analysis was performed on all PDX-tumor pairs. Selected PDX were harvested for assembly into 3D microtissues. Based on prior work in our labs, we evaluated drugs (trabectedin, olaparib, and mirdametinib) ex vivo and in vivo. For 3D microtissue studies, cell viability was the endpoint as assessed by Zeiss Axio Observer. For PDX drug studies, tumor volume was measured twice weekly. Bulk RNA sequencing was performed to identify pathways enriched in cells. We developed 13 NF1-associated MPNST-PDX and identified mutations or structural abnormalities in NF1 (100%), SUZ12 (85%), EED (15%), TP53 (15%), CDKN2A (85%), and chromosome 8 gain (77%). We successfully assembled PDX into 3D microtissues, categorized as robust (>90% viability at 48 h), good (>50%), or unusable (<50%). We evaluated drug response to “robust” or “good” microtissues, namely MN-2, JH-2-002, JH-2-079-c, and WU-225. Drug response ex vivo predicted drug response in vivo, and enhanced drug effects were observed in select models. These data support the successful establishment of a novel 3D platform for drug discovery and MPNST biology exploration in a system representative of the human condition.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

The Uranium-Containing and Thorium-Containing Anions Studied by Photoelectron Spectroscopy

An in-depth knowledge of actinide chemistry is fundamental to many aspects of nuclear science and technology, including the synthesis and processing of materials and the remediation of waste disposal sites. Among the actinides, the chemical bonding behaviors of actinium and thorium resemble those of the transition metals; the 5f-electrons of protactinium, uranium, neptunium, and plutonium often play important roles in their bonding; and among the still heavier elements, their bonding tends to mimic the lanthanide elements in terms of electron shielding and their f-electron contributions. Bonding that involves 5f-electrons, however, is especially important, in part because of the significance of uranium and plutonium, but also because these elements are among the few where f-electron participation in bonding is relatively common. This work focused on studying uranium-containing and thorium-containing anions in the gas phase using negative ion photoelectron spectroscopy. Since this technique directly probed valence electrons, it was uniquely positioned to address open questions regarding molecular bonding and electron configurations. A particularly important issue concerned how bonding in actinide-containing molecules was affected by modifications to their actinide atoms’ environment, i.e., due to their interaction with ligands. A closely related question was how actinide atoms’ suborbitals were qualitatively reordered and their energies quantitatively shifted as a result of their ligated environments. These were especially relevant issues in regard to uranium due to it having multiple possible oxidation states (OS) and the potential for 5f electron participation in bonding. The effects of ligands on oxidation states and 5f-orbital energies in uranium bonding was expected to be pronounced. Both ligands and excess electrons were seen as probes of actinide atoms within actinide-containing molecules. Our strategy for advancing knowledge of chemical bonding in the actinide-containing species utilized the synergy between experiments and theory, where in some cases experimental results validated theory and where in others computational results assisted in interpreting experiments. Calculations on actinide systems are terrifically challenging due to large spin-orbit interactions, relativistic effects, and just the sheer number of electrons involved. Even in the simplest species, e.g., U and U2, the most sophisticated, modern calculations carried out by the most experienced theorists often only approximate experimentally-measured values, such as electron affinities. For theory to provide confident predictions that can be used to solve real problems it needed an iterative and ultimately corrective mechanism by which its methods can develop further. Experiments can be used to identify when theory has failed; whereupon the subsequent process of using the experiment-theory interplay can be used to find the cause of the failure. Upon fixing it in one case, different test species can be proposed and studied by the experiment-theory combination to determine whether the problem has been corrected. Thus, experiments not only measure the values of molecular properties, they also provide navigational 3 beacons that keep computations off the reefs in an otherwise dark sea with few reference points. Experimental measurements in the actinide field are not only important, they are in actuality essential to computational progress. While it was not always possible to compare the theoreticallydetermined quantity of interest directly with the same experimentally-measured observable, it was usually possible to compare consequential properties that are both calculable and measurable. In the work completed here electron affinities and electronic state spacings were often sensitive consequential parameters. Reasonable agreement between measured and computational values signaled that a calculation that was very likely to be on-track. We had established collaborative relationships with five computational groups, all of which have expertise in computational actinide chemistry. Their PI’s are L. Cheng, D. Dixon, L. Gagliardi, K. Peterson, and B. Vlaisavljevich. Our close interaction with our theory partners led to us suggesting systems to them and them to us. This reciprocal interaction between our experimental and their computational results was among the most important strengths of this work and was a thread woven throughout. Even though anion photoelectron spectroscopic studies are conducted on anions, much of the information that they provide, pertains to the electronic structure of the neutral counterparts of those anions; among these are electron affinities and electronically excited state spacings. Our experimental tools included several specialized ion sources for forming the anionic species of interest, a mass spectrometer for identifying and mass-selecting them, and an anion photoelectron spectrometer for determining their electron affinities (EA) and characterizing the electronic states of the selected anions’ neutral counterparts. Anion photoelectron spectroscopy is conducted by crossing a mass-selected beam of anions with a fixed-frequency laser beam and energy-analyzing the resultant photodetached electrons. The photodetachment process is governed by the energyconserving relationship: hν = EBE + EKE, where hν is the photon’s energy, EBE is the electron binding (photodetachment transition) energy, and EKE is the electron’s kinetic energy. In our apparatus mass-selection is accomplished via time-of-flight mass spectrometry (TOF-MS), electron energy analysis is achieved with either a magnetic bottle or by velocity mapped imaging. Photodetachment of electrons from anions is implemented via either Nd:YAG or excimer lasers. The photodetachment transition energy, i.e., the EBE, between the ground vibrational and electronic state of an anion and the ground vibrational and electronic state of that anion’s neutral counterpart is the adiabatic electron affinity (EA) of that neutral molecule. Likewise, photodetachment transitions between the ground vibrational and electronic state of an anion and the various electronically-excited states of that anion’s corresponding neutral map the electronic spectrum of that neutral species, i.e., the spectral spacings in the photoelectron spectrum are a mirror image of the neutral’s electronic spectrum. It was, of course, crucial to be able to form the anionic species of interest. There, we had a particularly broad field of anion sources from which to choose. These included several variants of pulsed laser vaporization (LV), laser photoemission, infrared desorption plus photoemission, pulsed arc discharge (PACIS), electrospray ionization (ESI), and Rydberg electron transfer (RET). Each of these anion sources were readily combined with, i.e., connected to, the anion photoelectron spectroscopic portion of our apparatus as described above. Among the sources that utilize lasers, visible light for LV sources as well as IR for desorption sources are provided by Nd:YAG lasers. Ultraviolet photons are provided by both Nd:YAG and excimer lasers, whereas the excitation wavelengths for RET experiments come from two Nd:YAG-pumped dye lasers.

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

EVAP: A two-photon imaging tool to study conformational changes in endogenous Kv2 channels in live tissues

A primary goal of molecular physiology is to understand how conformational changes of proteins affect the function of cells, tissues, and organisms. Here, we describe an imaging method for measuring the conformational changes of the voltage sensors of endogenous ion channel proteins within live tissue, without genetic modification. We synthesized GxTX-594, a variant of the peptidyl tarantula toxin guangxitoxin-1E, conjugated to a fluorophore optimal for two-photon excitation imaging through light-scattering tissue. We term this tool EVAP (Endogenous Voltage-sensor Activity Probe). GxTX-594 targets the voltage sensors of Kv2 proteins, which form potassium channels and plasma membrane–endoplasmic reticulum junctions. GxTX-594 dynamically labels Kv2 proteins on cell surfaces in response to voltage stimulation. To interpret dynamic changes in fluorescence intensity, we developed a statistical thermodynamic model that relates the conformational changes of Kv2 voltage sensors to degree of labeling. We used two-photon excitation imaging of rat brain slices to image Kv2 proteins in neurons. We found puncta of GxTX-594 on hippocampal CA1 neurons that responded to voltage stimulation and retain a voltage response roughly similar to heterologously expressed Kv2.1 protein. Our findings show that EVAP imaging methods enable the identification of conformational changes of endogenous Kv2 voltage sensors in tissue.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Mitigate Stress Corrosion Cracking (SCC) in High-strength Al Castings (CRADA 508)

Development of high strength, corrosion resistant, aluminum alloys is important to a number of industries, including automotive, aerospace, power, etc. Understanding how processing methods such as casting, chemistry, additives and recycle content affect an alloy’s microstructure and performance in real world environments is crucial to their implementation. One major concern is how the repeated exposure to saline solutions can be detrimental to the structural integrity of aluminum alloy components. This study was undertaken as a collaboration between Eck Industries and PNNL to investigate the stress corrosion cracking (SCC) and corrosion behavior of aluminum alloys. This work examines the SCC response of three different materials set. First is cast A206 Al alloy produced by three different casting techniques (i.e., permanent mold, and sand casting with and without chill) and with or without nano forming additives as potential grain refiners. The second materials set comprises Al alloy tubes extruded by PNNL’s ShAPE process and produced using raw material with three different ratios of twitch and pre-consumer Al 6061 scrap. Finally, the third materials set comprises two Al-Si-Mg-Fe alloys with controlled Si% and Fe% to mimic recycle grade Al alloys. The corrosion behavior of these materials was characterized using SCC testing and scanning electrochemical cell microscopy (SECCM) technique; mechanical behavior was characterized using tensile testing, and microstructure was characterized using optical and electron microscopy techniques. Of the A206 Al alloys, the samples with nano forming additives performed poorly compared to those without. The best SCC performance (i.e. least degradation in mechanical properties) was demonstrated by A206 alloy sand cast with chill while the permanent mold sample showed the worst SCC performance. The SECCM was performed on A206 alloy sand cast with chill to understand the effect of microstructure and determine the location-dependent corrosion properties on grain boundaries (GB) and inside grains (IG) on the sample. The corrosion potential and current measured on GB and IG at various points on the sample established that the GB locations are more cathodic and corrosive than the IG locations. Further work needs to be conducted to investigate the mechanisms behind the observed dependence of SCC on casting technique and nano forming additives. Among the ShAPE extruded tubes, SECCM data suggests the decreasing order of tendency to corrode as follows: 100% twitch > 75% twitch > 50% twitch. Finally, both variants of the Al-Si-Mg-Fe alloys showed a large volume fraction of Fe- and Si-containing intermetallics and additional work is needed to discern differences in their respective corrosion responses. In summary, correlating local electrochemical behavior (e.g. via SECCM) with micro/nanostructure, in conjunction with bulk mechanical behavior, can help identify the right fabrication method and alloy chemistry to minimize SCC issues in Al alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Joint X-Ray, Kinetic Sunyaev–Zeldovich, and Weak Lensing Measurements: Toward a Consensus Picture of Efficient Gas Expulsion from Groups and Clusters

There is no consensus on how baryon feedback shapes the underlying matter distribution from either simulations or observations. We confront the uncertain landscape by jointly analyzing new measurements of the gas distribution around groups and clusters—DESI+ACT kinetic Sunyaev–Zel’dovich (kSZ) effect profiles and eROSITA X-ray gas masses—with mean halo masses characterized by galaxy–galaxy lensing. Across a wide range of halo masses ( M 500 = 10 13−14 M ⊙ ) and redshifts (0 < z < 1), we find evidence of more efficient gas expulsion beyond several R 500 than predicted by most state-of-the-art simulations. A like-with-like comparison reveals all kSZ and X-ray observations are inconsistent with the fiducial 1 Gpc 3 hydrodynamical FLAMINGO simulation, which was calibrated to reproduce pre-eROSITA X-ray gas fractions: eROSITA X-ray gas fractions are 2 × lower than the simulation, and the kSZ measurements are combined >8σ discrepant. The FLAMINGO simulation variant with the most gas expulsion, and therefore the most suppression of the matter power spectrum relative to a dark-matter-only simulation, provides a good description of how much gas is expelled and how far it extends; the enhanced gas depletion is achieved by more powerful but less frequent AGN outbursts. Joint kSZ, X-ray, and lensing measurements form a consistent picture of gas expulsion beyond several R 500 , implying a more suppressed matter power spectrum than predicted by most recent simulations. Complementary observables (e.g., thermal Sunyaev–Zel’dovich effect and fast radio bursts) and next-generation simulations are critical to understanding the physical mechanism behind this extreme gas expulsion and mapping its impact on the large-scale matter distribution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-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 near Trail Creek. 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 TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores 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 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 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

Integrative analysis of the 3D genome and epigenome in mouse embryonic tissues

While a rich set of putative cis-regulatory sequences involved in mouse fetal development have been annotated recently on the basis of chromatin accessibility and histone modification patterns, delineating their role in developmentally regulated gene expression continues to be challenging. To fill this gap, here we mapped chromatin contacts between gene promoters and distal sequences across the genome in seven mouse fetal tissues and across six developmental stages of the forebrain. We identified 248,620 long-range chromatin interactions centered at 14,138 protein-coding genes and characterized their tissue-to-tissue variations and developmental dynamics. Integrative analysis of the interactome with previous epigenome and transcriptome datasets from the same tissues revealed a strong correlation between the chromatin contacts and chromatin state at distal enhancers, as well as gene expression patterns at predicted target genes. We predicted target genes of 15,098 candidate enhancers and used them to annotate target genes of homologous candidate enhancers in the human genome that harbor risk variants of human diseases. We present evidence that schizophrenia and other adult disease risk variants are frequently found in fetal enhancers, providing support for the hypothesis of fetal origins of adult diseases.

59 BASIC BIOLOGICAL SCIENCES↗