Search NASA⌕ Search

SEARCH · Search NASA

Results for “root chemistry”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Data for "Depth of nutrient uptake by deep-rooted plants is regulated by water availability"

The data set consists of strontium (Sr) isotope ratios (87Sr/86Sr), water isotopes, soil cation concentrations, soil water potential sensor data, and results of 87Sr/86Sr mixing model. The plant canopy size files include the dataset of canopy dimension of sagebrush, lupine, and sunflower. The soil and plant ICPMS (Inductively Coupled Plasma Mass Spectrometry) data file includes both of 87Sr/86Sr, and cation concentration dataset from soil exchangeable pool, apatite pool, silicate extract, atmospheric rain deposition, and plant leaf and stem tissues. The plant dendrochronology file includes the dendrochronogical ring width of several sagebrush, and dendrochemical sample data includes the 87Sr/86Sr for each separated growth ring. The modeling result gives the proportion of nutrient sources of each plants (based on their 87Sr/86Sr in leaf tissues and growth rings) from atmospheric deposition and mineral weathering. Soil water potential data includes continuous collection of soil water potential dataset at 2 depths (30 cm and 60 cm, from Nov 24 - Jun 25) of the sampling site. All the samples were collected from 2 sampling campaign June and July 2023, and rain water is a separate sampling from Aug - Sept 2023, at north-facing hillslope near pumphouse site. The data showed that the depth of cation nutrient acquisition is thus tightly coupled with, and likely determined by, water availability in soil, saprolite and bedrock. The enhanced uptake of cations and water from regions of mineral weathering could confer plant and ecosystem resilience during low water years and may impact the rate of bedrock weathering and watershed chemistry during drought. This dataset includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type; a location metadata file (locations.csv); and a samples metadata file (samples.csv). All files are provided as comma-separated values (CSV) files (.csv). 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.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Leveraging design of experiments to build chemometric models for the quantification of uranium (VI) and HNO3 by Raman spectroscopy

Partial least squares regression (PLSR) and support vector regression (SVR) models were optimized for the quantification of U(VI) (10–320 g L −1 ) and HNO 3 (0.6–6 M) by Raman spectroscopy with optimized calibration sets chosen by optimal design of experiments. The designed approach effectively minimized the number of samples in the calibration set for PLSR and SVR by selecting sample concentrations with a quadratic process model, despite complex confounding and covarying spectral features in the spectra. The top PLS2 model resulted in percent root mean square errors of prediction for U(VI), HNO 3 , and NO 3 − of 3.7%, 3.6%, and 2.9%, respectively. PLS1 models performed similarly despite modeling an analyte with a majority linear response (i.e., uranyl symmetric stretch) and another with more covarying vibrational modes (i.e., HNO 3 ). Partial least squares (PLS) model loadings and regression coefficients were evaluated to better understand the relationship between weaker Raman bands and covarying spectral features. Support vector machine models outperformed PLS1 models, resulting in percent root mean square error of prediction values for U(VI) and HNO 3 of 1.5% and 3.1%, respectively. The optimal nonlinear SVR model was trained using a similar number of samples (11) compared with the PLSR model, even though PLS is a linear modeling approach. The generic D-optimal design presented in this work provides a robust statistical framework for selecting training set samples in disparate two-factor systems. This approach reinforces Raman spectroscopy for the quantification of species relevant to the nuclear fuel cycle and provides a robust chemometric modeling approach to bolster online monitoring in challenging process environments.

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

Plants and mycorrhizal symbionts acquire substantial soil nitrogen from gaseous ammonia transport

Nitrogen (N) is an essential nutrient that limits plant growth in many ecosystems. Here we investigate an overlooked component of the terrestrial N cycle – subsurface ammonia (NH 3 ) gas transport and its contribution to plant and mycorrhizal N acquisition. In this work, we used controlled mesocosms, soil incubations, stable isotopes, and imaging to investigate edaphic drivers of NH 3 gas efflux, track lateral subsurface N transport originating from 15 NH 3 gas or 15 N-enriched organic matter, and assess plant and mycorrhizal N assimilation from this gaseous transport pathway. NH3 is released from soil organic matter, travels belowground, and contributes to root and fungal N content. Abiotic soil properties (pH and texture) influence the quantity of NH 3 available for subsurface transport. Mutualisms with arbuscular mycorrhizal (AM) fungi can substantially increase plant NH 3 -N uptake. The grass Brachypodium distachyon acquired 6–9% of total plant N from organic matter-N that traveled as a gas belowground. Colonization by the AM fungus Rhizophagus irregularis was associated with a two-fold increase in total plant N acquisition from subsurface NH 3 gas. NH 3 gas transport and uptake pathways may be fundamentally different from those of more commonly studied soil N species and warrant further research.

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

Evaluation of CMIP6 models in simulating the statistics of extreme precipitation over Eastern Africa

We report the Eastern Africa region experiences frequent extreme precipitation events that can cause destruction of property and environment, and loss of lives. Thus, there is a need to understand how these events may change in the future and how well the global climate models that are used to make projections can simulate precipitation extremes in this region before they can be used in downscaling or flood and drought impact assessment studies. In this work, we evaluated the ability of sixteen Coupled Model Intercomparison Project Phase 6 (CMIP6) models to simulate present-day precipitation extremes over the Eastern Africa region during the two rainy seasons (March–May and September–November). We used nine extreme precipitation indices (including seven (one) indices of wet (dry) extremes) defined by the Expert Team on Climate Change Detection and Indices. The CMIP6 models were evaluated against two gridded observation datasets: Global Precipitation Climatology Project One-Degree Daily Dataset and Tropical Rainfall Measuring Mission Multi-satellite Precipitation Analysis 3B42. Three model performance metrics (percentage bias, normalized root-mean-square error, and pattern correlation coefficient) were employed to further assess the strengths and weakness of the models. Our results show that the multi-model ensemble mean generally provides a better representation of observed precipitation and related extremes compared to individual models when considering all metrics and seasons. Several consistent biases are evident across CMIP6 models, which tend to overestimate the total-wet day precipitation and consecutive wet days, and underestimate very wet days and maximum 5-day precipitation in both seasons. Furthermore, no single model consistently performs best, model performance varies with the season and index under consideration and is generally independent of horizontal resolution.

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

Programmable bionanocomposite coated fertilizers for prolonged controlled release of nitrogen

Controlled-release fertilizers (CRFs) present a promising solution for alleviating food and nutrient scarcity. However, their development has been mainly hindered by both rapid and unsynchronized nutrient release and unsustainable coating materials. In this study, we address this issue by developing a new coating layer for CRFs using a programmable biopolyurethane nanocomposite. This nanocomposite is prepared from diphenylmethane diisocyanate (MDI)-functionalized bentonite nanoclay (BNT-MDI) and a biopolyol from biomass waste. The results show that the BNT-MDI-doped CRFs (BCRFs) exhibit an impressive nitrogen (N)-release longevity of approximately 120 days at a 4 wt% coating ratio, surpassing previous CRF formulations. In a 30-day snap bean cultivation study, BCRF significantly improved root length (1000%), leaf length (257%), leaf width (400%), and plant height (1400%) compared to the control. The superior performance of BCRF is attributed to the PU-nanoclay biocomposite film, with full nano-exfoliation, controllable porosity, and high crosslinking density. Furthermore, we introduce a new dynamic release mechanism and establish a quantitative relationship between the nanostructure, property, and release performance of BCRF by combining the multiplicative and diffusion models for porous materials. Finally, this study provides a theoretical framework and a straightforward methodology for designing programmable nanocomposite structures for future biobased CRFs.

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

Nutrient and moisture limitations reveal keystone metabolites linking rhizosphere metabolomes and microbiomes

Plants release a wealth of metabolites into the rhizosphere that can shape the composition and activity of microbial communities in response to environmental stress. The connection between rhizodeposition and rhizosphere microbiome succession has been suggested, particularly under environmental stress conditions, yet definitive evidence is scarce. In this study, we investigated the relationship between rhizosphere chemistry, microbiome dynamics, and abiotic stress in the bioenergy crop switchgrass grown in a marginal soil under nutrient-limited, moisture-limited, and nitrogen (N)-replete, phosphorus (P)-replete, and NP-replete conditions. We combined 16S rRNA amplicon sequencing and LC-MS/MS-based metabolomics to link rhizosphere microbial communities and metabolites. We identified significant changes in rhizosphere metabolite profiles in response to abiotic stress and linked them to changes in microbial communities using network analysis. N-limitation amplified the abundance of aromatic acids, pentoses, and their derivatives in the rhizosphere, and their enhanced availability was linked to the abundance of bacterial lineages from Acidobacteria, Verrucomicrobia, Planctomycetes, and Alphaproteobacteria. Conversely, N-amended conditions increased the availability of N-rich rhizosphere compounds, which coincided with proliferation of Actinobacteria. Treatments with contrasting N availability differed greatly in the abundance of potential keystone metabolites; serotonin and ectoine were particularly abundant in N-replete soils, while chlorogenic, cinnamic, and glucuronic acids were enriched in N-limited soils. Serotonin, the keystone metabolite we identified with the largest number of links to microbial taxa, significantly affected root architecture and growth of rhizosphere microorganisms, highlighting its potential to shape microbial community and mediate rhizosphere plant–microbe interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Fermi-Based Kinetic Model for the Sabatier Reaction: Sabatier Principle and Beyond It

The CO 2 methanation reaction (aka the Sabatier reaction) is one of the important carbon utilization technologies to reduce atmospheric greenhouse gases. In this study, kinetic modeling of the Sabatier reaction proceeding via the CO route has been presented by considering weak and strong chemisorption concepts rooted in electronic theory. Based on the single-route mechanism of the Sabatier reaction, the dependence of the reaction rate on the position of the catalyst Fermi level with respect to the antibonding states of the reagent (CO 2 ) and the product (CH 4 ), is demonstrated with the classic Sabatier volcano-shape. Furthermore, the model was applied to explain experimental results of reaction rates for Ni- and Ru-based catalytic systems and to understand origin of enhancement of turnover frequencies established for dynamic catalysis.

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

Molecular Dynamics Simulation of Atomic Interactions in the Vancomycin Binding Site

Vancomycin is a glycopeptide antibiotic produced by Amycolaptopsis orientalis used to treat serious infections by Grampositive pathogens including methicillin-resistant Staphylococcus aureus. Vancomycin inhibits cell wall biosynthesis by targeting lipid II, which is the membrane-bound peptidoglycan precursor. The heptapeptide aglycon structure of vancomycin binds to the D-Ala-D-Ala of the pentapeptide stem structure in lipid II. The third residue of vancomycin aglycon is asparagine, which is not directly involved in the dipeptide binding. Nonetheless, asparagine plays a crucial role in substrate recognition, as the vancomycin analogue with asparagine substituted by aspartic acid (V D ) shows a reduction in antibacterial activities. To characterize the function of asparagine, binding of vancomycin and its aspartic-acid-substituted analogue V D to L-Lys-D-Ala-D-Ala and L-Lys-D-Ala-D-Lac was investigated using molecular dynamic simulations. Binding interactions were analyzed using root-mean-square deviation (RMSD), two-dimensional (2D) contour plots, hydrogen bond analysis, and free energy calculations of the complexes. The analysis shows that the aspartate substitution introduced a negative charge to the binding cleft of V D , which altered the aglycon conformation that minimized the repulsive lone pair interaction in the binding of a depsipeptide. Our findings provide new insight for the development of novel glycopeptide antibiotics against the emerging vancomycin-resistant pathogens by chemical modification at the third residue in vancomycin to improve its binding affinity to the D-Ala-D-Lac-terminated peptidoglycan in lipid II found in vancomycin-resistant enterococci and vancomycinresistant S. aureus.

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

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

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

Comparing Designed Training Sets to Optimize Multivariate Regression Models for Pr, Nd, and Nitric Acid Using Spectrophotometry

Chemometric regression models were developed for the quantification of praseodymium (Pr, 0–1000 µg/mL), neodymium (Nd, 0–1000 µg/mL), and nitric acid (HNO 3 , 0.1–5 M) using spectrophotometry. Designed calibration sets were composed of 20 samples each: 10 model points and 10 lack-of-fit (LOF) points. The D-optimal designs effectively minimized the number of samples required to build models, and each design resulted in similar prediction performance, suggesting that statistical design of experiments can provide a reliable framework for selecting training set samples in three-variable systems. Partial least squares regression (PLSR) models were validated against a one-factor-at-a-time validation set composed of 125 samples (three variables, five levels). The top PLS-1 models resulted in average percent root mean square error of prediction error values of 3.5%, 1.7%, and 1.2% for Pr(III), Nd(III), and HNO 3 , respectively. Power set augmentations of the model and LOF samples were investigated to optimize the number of training set samples. PLSR models built using just required model points (10) had similar predictive capabilities as models including the LOF points (20) but with fewer samples. The number of validation samples was also varied systematically to learn how many samples are needed to validate regression models. This work addresses long-standing questions in the field of chemometrics to help make this approach amenable to the near-real-time quantification of hazardous species in remote settings.

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

Innovative Method for Welding in Generation 3 CSP to Enable Reliable Manufacturing of Solar Receivers to withstand Daily Cycling at Temperatures Above 700°C (Final Technical Report)

Inconel® Alloy 740H® (alloy 740H) was the first age-hardenable nickel-based alloy approved by the ASME Boiler & Pressure Vessel Code for use in pressure-boundary applications. Over the past ~20 years the alloy has been optimized for weldability and high-temperature stability, approved for use in different applications. Development of a supply chain combined with the advantageous properties of the alloy (high-temperature creep strength, oxidation and corrosion resistance, etc.) have resulted in the alloy being applied to new high-temperature power cycle demonstration projects, and of particular interest are applications to concentrating solar power (CSP) to enable higher-efficiency Generation 3 CSP systems and the corresponding supercritical CO 2 (sCO 2 ) power cycle components (heat exchangers, piping, etc.). The high allowable stresses of alloy 740H also make it a desirable material for current Generation 2 CSP solar power receivers to improve cyclic capability and/or reduce receiver height. Recent experiences in demonstration projects utilizing alloy 740H identified cracking issues during welding and fabrication. In this project, a detailed study was done to confirm and clarify the Stress Relaxation Cracking (SRxC) mechanism, also known as stress relief cracking or strain-age cracking (SAC), during post-weld heat-treatment (PWHT). This involved detailed microscopy and advanced characterization to understand the root cause(s) of three failures obtained from industry. Based in-part on these findings, a targeted laboratory based SRxC test method was utilized to evaluate variables such as heat-to-heat variations, strain level, PWHT temperature, and starting material condition on three heats of alloy 740H. Industrial shop welding of cold-worked plates was also conducted. The research showed the following: SRxC was confirmed as the cracking mechanism for all field failures; Stress state (from residual stresses, constraint, deformation, and local stress concentrations) was playing a significant role in field failures and laboratory testing confirmed increasing susceptibility for all heats with increasing strain levels. High levels of microstructural strain were identified at crack initiation locations, in some cases leading to local recrystallization; Precipitate free zones (PFZs) at grain boundaries were found at relaxation cracks and crack initiation locations uniquely associated with SRxC in alloy 740H. Laboratory testing reproduced this microstructural feature which had only previously been reported in long-term creep testing of weldments. Advanced nano-scale characterization confirmed the presence of a moving boundary leading to coarsening of precipitates and PFZs where damage accumulated; The research suggested heat-to-heat variations due to local chemistry and processing may influence SRxC susceptibility, but more work is needed to fully clarify these effects. To disseminate the key learnings from this research to the scientific and engineering communities and alloy 740H end users, multiple technical publications and presentations were made, an industrial alloy 740H users meeting was held, and a new industry guideline specification document which can be directly implemented by end-users of alloy 740H was produced.

14 SOLAR ENERGY↗
DOE OSTI2020

General Protocol for the Accurate Prediction of Molecular 13 C/ 1 H NMR Chemical Shifts via Machine Learning Augmented DFT

An accurate prediction of NMR chemical shifts at affordable computational cost is very important for different types of structural assignments in experimental studies. Density functional theory (DFT) and gauge-including atomic orbital (GIAO) are two of the most popular computational methods for NMR calculation, yet they often fail to resolve ambiguities in structural assignments. In this work, we present a new method that uses machine learning (ML) techniques (DFT + ML) that significantly increases the accuracy of 13 C/ 1 H NMR chemical shift prediction for a variety of organic molecules. The input of the generalizable DFT + ML model contains two critical parts: one is a vector providing insights into chemical environments, which can be evaluated without knowing the exact geometry of the molecule; the other one is the DFT-calculated isotropic shielding constant. The DFT + ML model was trained with a data set containing 476 13 C and 270 1 H experimental chemical shifts. For the DFT methods used here, the root mean square deviations (RMSDs) for the errors between predicted and experimental 13 C/ 1 H chemical shifts can be as small as 2.10/0.18 ppm, which is much lower than those from simple DFT (5.54/0.25 ppm), or DFT + linear regression (LR) (4.77/0.23 ppm) approaches. It also has a smaller maximum absolute error than two previously proposed NMR-predicting ML models. The robustness of the DFT + ML model is tested on two classes of organic molecules (TIC10 and hyacinthacines), where the correct isomers were unambiguously assigned to the experimental ones. Overall, the DFT + ML model shows promise for structural assignments in a variety of systems, including stereoisomers, that are often challenging to determine experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Metallographic Preparation of Space Shuttle Reaction Control System Thruster Electron Beam Welds for Electron Backscatter Diffraction

A Space Shuttle Reaction Control System (RCS) thruster failed during a firing test at the NASA White Sands Test Facility (WSTF), Las Cruces, New Mexico. The firing test was being conducted to investigate a previous electrical malfunction. A number of cracks were found associated with the fuel closure plate/injector assembly (Fig 1). The firing test failure generated a flight constraint to the launch of STS-133. A team comprised of several NASA centers and other research institutes was assembled to investigate and determine the root cause of the failure. The JSC Materials Evaluation Laboratory was asked to compare and characterize the outboard circumferential electron beam (EB) weld between the fuel closure plate (Titanium 6Al-4V) and the injector (Niobium C-103 alloy) of four different RCS thrusters, including the failed RCS thruster. Several metallographic challenges in grinding/polishing, and particularly in etching were encountered because of the differences in hardness, ductility, and chemical resistance between the two alloys and the bimetallic weld. Segments from each thruster were sectioned from the outboard weld. The segments were hot-compression mounted using a conductive, carbon-filled epoxy. A grinding/polishing procedure for titanium alloys was used [1]. This procedure worked well on the titanium; but a thin, disturbed layer was visible on the niobium surface by means of polarized light. Once polished, each sample was micrographed using bright field, differential interference contrast optical microscopy, and scanning electron microscopy (SEM) using a backscatter electron (BSE) detector. No typical weld anomalies were observed in any of the cross sections. However, areas of large atomic contrast were clearly visible in the weld nugget, particularly along fusion line interfaces between the titanium and the niobium. This prompted the need to better understand the chemistry and microstructure of the weld (Fig 2). Energy Dispersive X-Ray Spectroscopy (EDS) was used to confirm the chemical composition of the variations in contrast in these areas. Niobium alloys generally require exposure to more aggressive chemical reagents than titanium alloys for etching because of niobium s chemical resistance; therefore, the titanium portion of the sample was etched first. A five second immersion in Kroll s reagent revealed a general microstructure on the titanium portion of the sample; however, the titanium heat affected zone closest to the weld, was over-etched due to higher concentrations of refined grains and an increase in eta-phase. The Kroll s etchant also revealed some microstructure in the weld nugget itself; the niobium portion of the sample remained unetched.

Martinez, James↗
DOE OSTI2020

Deep Eutectic Solvents: A Review of Fundamentals and Applications

Deep eutectic solvents (DESs) are an emerging class of mixtures characterized by significant depressions in melting points compared to those of the neat constituent components. These green materials are promising for applications as inexpensive ’designer’ solvents exhibiting a host of tunable physicochemical properties. Here, a detailed review of the current literature reveals the lack of predictive understanding of the microscopic mechanisms that govern the structure-property relationships in this class of solvents. Complex hydrogen bonding behavior is postulated as the root cause of their melt-ing point depressions and physicochemical properties; to understand these hydrogen bonded networks, it is imperative to study these systems as dynamic entities using both simulations and experiments. This review emphasizes recent research efforts in order to elucidate the next steps needed to develop a fundamental framework for understanding DESs. It covers recent developments in DES research, frames outstanding scientific questions, and identifies gainful research thrusts that would advance the field towards predictive models and fundamental understanding of these solvents.

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

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Analysis of Dynamical Peculiarities in Nanoalloys at Subsystems Level: Dynamical Degrees of Freedom, Temperature Differences, and the Chameleon Effect

Abstract A novel analysis of the dynamical behavior of nanoalloy systems, as represented by model Ni/Al 13‐atom clusters, over a broad range of energies that cover the stage‐wise transition of the systems from their solid‐like to liquid‐like state is presented. Conceptually, the analysis is rooted in partitioning the systems into judiciously chosen subsystems and characterizing the latter in terms of subsystem‐specific dynamical descriptors that include dynamical degrees of freedom , root‐mean‐square bond‐length fluctuation , and element‐specific subsystem temperature . The analysis reveals a host of intriguing new peculiarities in the dynamical behavior of the Ni/Al 13‐mers, among which are what we call the chameleon effect and the difference in the temperatures of the Ni and Al subsystems at high energies, a difference that strongly depends on the cluster composition and also changes with energy. These do not have an analog in pure Ni 13 and Al 13 and are explained in terms of the coupled effects of the difference between the masses of the Ni and Al atoms (the mass effect) and of the difference in the anharmonicity of the overall interaction potential as experienced by the Ni and Al subsystems of the clusters (the potential effect).

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

A classical model for three-body interactions in aqueous ionic systems

We present a classical induction model to evaluate the three-body ion–water–water (I–W–W) and water–water–water (W–W–W) interactions in aqueous ionic systems. The classical description of the induction energy is based on electrostatic distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole on the O and H atoms of water. The monatomic ions were described by a point charge and a dipole–dipole polarizability, while for the polyatomic ions, distributed multipoles up to hexadecapole and distributed polarizabilities up to quadrupole–quadrupole were used. The accuracy of the classical model is benchmarked against an accurate dataset of 936 (I–W–W) and 2184 (W–W–W) three-body terms for 13 different monatomic and polyatomic cation and anion systems. The classical model shows excellent agreement with the reference second order Moller–Plesset and coupled-cluster single double and perturbative triple [CCSD(T)] three-body energies. The Root-Mean-Square-Errors (RMSEs) for monatomic cations, monatomic anions, and polyatomic ions were 0.29, 0.25, and 0.12 kcal/mol, respectively. The corresponding RMSE for 1744 CCSD(T)/aVTZ three-body (W–W–W) energies, used to train MB-pol, was 0.12 kcal/mol. The accuracy of the proposed classical model demonstrates that the three-body term for aqueous ionic systems can be accurately modeled classically. This approach provides a fast, efficient, and as-accurate path toward modeling the three-body term in aqueous ionic systems that is fully transferable across systems with different ions without the need to fit to tens of thousands of ab initio calculations for each ion to extend existing many-body force fields to interactions between water and ions.

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

Microbially mediated iron redox cycling of subsurface sediments from Hanford Site, Washington State, USA

Iron-bearing clay minerals are predominant in soils and sediments, and they undergo oxidation-reduction cycles as a result of natural processes such as wetting/drying cycles and plant root respiration. However, the kinetics and mechanisms of multiple redox cycles of iron in clay-rich sediments and consequences of such cycling on sediment properties are poorly understood. The objective of this study was to understand how multiple redox cycles of Fe in clay-rich sediment affect the rate and extent of Fe bioreduction and the physicochemical properties of the sediment. A natural sediment sample containing Fe-bearing montmorillonite from Hanford, Washington, USA was size-fractionated [2.0 to 0.5 μm (Hanford-C) and 0.5 to 0.02 μm (Hanford-F)] and redox-cycled for four times. Bioreduction was achieved utilizing Geobacter sulfurreducens and re-oxidization was performed with sparged air. Time-course change of total Fe (II) was monitored to measure the rate and extent of Fe (III) bioreduction. Redox-cycled sediments were characterized to determine the physicochemical changes. Both the initial rate and extent of bioreduction fluctuated across the four redox cycles, but they ultimately decreased from 5.3 μmol g -1 h -1 and 22.9% to nearly zero by the fourth cycle. These fluctuation patterns were likely due to a combined effect of reductive dissolution of small/poorly crystalline clay particles (by 3–5%) and clay mineral structural changes, as evidenced by redox induced changes of aqueous chemistry, surface area, cation exchange capacity, mineralogy, and Mössbauer parameters. Once these small/poorly crystalline clay particles were dissolved, structural Fe in residual larger and more crystalline clay particles was largely reversible across additional redox cycles, as revealed by Mössbauer spectroscopy through the first three cycles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗