Soil carbon dynamics during drying vs. rewetting: Importance of antecedent moisture conditions
Not Available
Patel, Kaizad F. (ORCID:0000000164377809)↗
SEARCH · Engineering Papers
Results for “DYNAMICS”
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.
Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data
The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.
Prange, Micah P.↗
Molecular Dynamics Simulation of the Structures, Dynamics, and Aggregation of Dissolved Organic Matter
Dissolved organic matter (DOM) plays a significant role in the transport and transformation of pollutants in the aquatic environment. However, the experimental characterization of DOM has been limited mainly to bulk properties, and the molecular-level interactions among various components of DOM remain to be fully characterized. In this study, we use molecular dynamics (MD) simulations to probe the structural properties of model DOM systems at atomic detail. The 200 ns simulations, validated by available experimental data, reveal processes and mechanisms by which chemical species (cations, peptides, lipids, lignin, carbohydrates, and some low-molecular-weight aliphatic and aromatic compounds) aggregate to form complex DOM. The DOM aggregates are dynamic, consisting of a hydrophobic core and amphiphilic exterior. The lipid tails and other hydrophobic fragments form the core, with hydrophilic and amphiphilic groups exposed to water, making DOM accessible to both polar and nonpolar species. Thus, the lipid component acts as a nucleator, whereas cations (especially Ca 2+ ) connect the molecular fragments on the surface by coordinating with the O-containing functional groups of DOM. The structural details revealed here provide new insights including surface accessible atoms, overall assemblage, and interactions among the molecules of DOM for understanding the kinetics and mechanisms through which DOM interacts with metal and other contaminants.
54 ENVIRONMENTAL SCIENCES↗
Growth Kinetics and Size Distribution Dynamics of Viscous Secondary Organic Aerosol
Low bulk diffusivity inside viscous semisolid atmospheric secondary organic aerosol (SOA) can prolong equilibration time scale, but its broader impacts on aerosol growth and size distribution dynamics are poorly understood. In this article, we present quantitative insights into the effects of bulk diffusivity on the growth and evaporation kinetics of SOA formed under dry conditions from photooxidation of isoprene in the presence of a bimodal aerosol consisting of Aitken (ammonium sulfate) and accumulation (isoprene or α-pinene SOA) mode particles. Aerosol composition measurements and evaporation kinetics indicate that isoprene SOA is composed of several semivolatile organic compounds (SVOCs), with some reversibly reacting to form oligomers. Model analysis shows that liquid-like bulk diffusivities can be used to fit the observed evaporation kinetics of accumulation mode particles but fail to explain the growth kinetics of bimodal aerosol by significantly under-predicting the evolution of the Aitken mode. In contrast, the semisolid scenario successfully reproduces both evaporation and growth kinetics, with the interpretation that hindered partitioning of SVOCs into large viscous particles effectively promotes the growth of smaller particles that have shorter diffusion time scales. This effect has important implications for the growth of atmospheric ultrafine particles to climatically active sizes.
54 ENVIRONMENTAL SCIENCES↗
Modeling Nonperturbative Field-Driven Vibronic Dynamics: Selective State Preparation and Nonlinear Spectroscopy
The partially linearized density matrix formalism for nonadiabatic dynamics is adapted to incorporate a classical external electromagentic field into the system Hamiltonian. This advancement encompasses the possibility of describing field-driven dynamics and computing a variety of linear and nonlinear spectroscopic signals beyond the perturbative limit. Here, the capabilities of the developed approach are demonstrated on a simple two-state vibronic model coupled to a bath, for which we (a) perform an exhaustive search in the field parameter space for optimal state preparation and (b) compute time-resolved transient absorption spectroscopy to monitor the effect of different pulse shapes on measurable experimental signals. While no restrictions on the form of the field have to be assumed, we focus here on Gaussian shaped (linearly) chirped pulses.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Modeling the Ultrafast Electron Attachment Dynamics of Solvated Uracil
Electron attachment to DNA by low energy electrons can lead to DNA damage, so a fundamental understanding of how electrons interact with the components of nucleic acids in solution is an open challenge. In solution, low energy electrons can generate presolvated electrons, e pre – , which are efficiently scavanged by pyrimidine nucleobases to form transient negative ions, able to relax to either stable valence bound anions or undergo dissociative electron detachment or transfer to other parts of DNA/RNA leading to strand breakages. In order to understand the initial electron attachment dynamics, this paper presents a joint molecular dynamics and high-level electronic structure study into the behavior of the electronic states of the solvated uracil anion. Both the valence π* and nonvalence e pre – states of the solvated uracil system are studied, and the effect of the solvent environment and the geometric structure of the uracil core are uncoupled to gain insight into the physical origin of the stabilization of the solvated uracil anion. As a result, solvent reorganization is found to play a dominant role followed by relaxation of the uracil core.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Solvation Dynamics of Wet Ethaline: Water is the Magic Component
The past two decades witnessed the development of a new type of solvent system, named deep eutectic solvents, which have become increasingly investigated because they offer new and potentially favorable properties, such as wide tunability in electrochemical, mechanical, and transport properties. Deep eutectic solvent (DES) systems are composed of at least one main solvent and an additional component that is meant to interrupt the original solvent/solvent interactions, thereby introducing lower melting points relative to each individual component. Ethaline (a 1:2 mole% mixture of choline chloride and ethylene glycol) is one of the most promising DES systems. However, it is also known to be very hygroscopic, which is a constant concern because water absorption during the use of ethaline alters its properties. Within this work, we demonstrate that modest amounts of water addition (1-10%) to ethaline are of little concern for practical use and can even lead to performance improvements, such as accelerated relaxation and solvation. In contrast, very small amounts of <1% of water lead to additional slowing of the solvent response. Thus, we suggest that the attempt to dry ethaline below 1% moisture is rather counter-productive if one attempts to achieve effective solvation and charge transport properties from DESs. This study investigates the effect of water content on the diffusional relaxation dynamics of ethaline. A set of independent spectroscopic experiments and computational simulations are aimed to provide insight into the solvent response of the DES system using femtosecond time-resolved absorption spectroscopy (fs-TA), broadband dielectric spectroscopy (BDS), nuclear magnetic resonance (NMR) diffusometry and broadband relaxometry, and molecular dynamics simulations (MDS) on ethaline with 0, 0.1, 1, 10, and 28.5wt% added water. For dry ethaline, we identify choline chloride as the rate-limiting solvation component in ethaline. However, the role of the solvent components changes gradually as water is added. Further, we provide quantitative solvent relaxation rates using the different presented time-resolved spectroscopic techniques and find remarkable agreement between them. Based on the solvent relaxation rates and combined with MDS, we develop a molecular understanding of the individual solvent components and their interactions in dry and wet ethaline with varying amounts of water content.
25 ENERGY STORAGE↗
Molecular Dynamics Simulations of the Interfacial Region between Boehmite and Gibbsite Basal Surfaces and High Ionic Strength Aqueous Solutions
Classical molecular dynamics (MD) simulations were used to study the interactions of up to 2 M NaCl and NaNO3 aqueous solutions with the presumed inert boehmite (010) and gibbsite (001) surfaces. The force field parameters used in these simulations were validated against density functional theory calculations of Na+ and Cl- hydrated complexes adsorbed at the boehmite (010) surface. In all the classical MD simulations and regardless of the ionic strength or the nature of the anion, Na+ ions were found to preferably form inner-sphere complexes over outer-sphere complexes at the aluminum (oxy)hydroxide surfaces, adsorbing closer to the surface than both water molecules and anions. In contrast, Cl- ions were distributed almost equally between inner- and outer-sphere positions. The resulting asymmetry in adsorption strengths offers molecular-scale evidence for the observed isoelectric point (IEP) shift to higher pH at high ionic strength for aluminum (oxy)hydroxides. As such, the MD simulations also provided clear evidence against the assumption that the basal surfaces of boehmite and gibbsite are inert to background electrolytes. Finally, the MD simulations indicated that, although the adsorption behavior of Na+ in NaNO3 and NaCl solutions was similar, the different affinities of NO3- and Cl- for the aluminum (oxy)hydroxide surfaces might have macroscopic consequences, such as difference in the sensitivity of the IEP to the electrolyte concentration.
Shen, Zhizhang↗
Machine Learning Classical Interatomic Potentials for Molecular Dynamics from First-Principles Training Data
The ever-increasing power of modern super-computers, along with the availability of highly scalable atomistic simulation codes, has begun to revolutionize predictive modeling of materials. In particular, molecular dynamics (MD) has led to breakthrough advances in diverse fields, including tribology, catalysis, sensing, and nanoparticle self-assembly. Furthermore, recent integration of MD simulations with X-ray characterization has demonstrated promise in real-time 3-D characterization of materials on the atomic scale. The popularity of MD is driven by its applicability at disparate length/time scales, ranging from ab initio MD (hundreds of atoms and tens of picoseconds) to all-atom classical MD (millions of atoms over microseconds), and coarse-grained (CG) models (micrometers and tens of microseconds). Nevertheless, a substantial gap persists between AIMD, which is highly accurate but restricted to extremely small sizes, and those based on classical force fields (atomistic and CG) with limited accuracy but access to larger length/time scales. The accuracy and predictive power of classical MD simulations is dictated by the empirical force fields, and their capability to capture the relevant physics. Here, we discuss some of our recent work on the use of machine learning (ML) to combine the accuracy and flexibility of electronic structure calculations with the speed of classical potentials. Our ML framework attempts to bridge the significant gulf that exists between the handful of research groups that develop new interatomic potential models (often requiring several years of effort), and the increasingly large user community from academia and industry that applies these models. Our data-driven approach represents significant departure from the status quo and involves several steps including generation and manipulation of extensive training data sets through electronic structure calculations, defining novel potential functional forms, employing state-of-the-art ML algorithms to formulate highly optimized training procedures, and subsequently developing user-friendly workflow tools integrating these algorithms on high-performance computers (HPCs). Our ML approach shows marked success in developing force fields for a wide range of materials from metals, oxides, nitrides, and heterointerfaces to two-dimensional (2D) materials.
Chan, Henry↗
Molecular Rotational Dynamics in Mixed CH 4 –CO 2 Hydrates: Insights from Molecular Dynamics Simulations
The impact of guest molecule composition on the rotational dynamics in CH 4 , CO 2 , and mixed CH 4 –CO 2 gas hydrates is investigated with classical molecular dynamics simulations. Rotational autocorrelation functions are calculated for the guest and host molecules in each hydrate composition from simulation trajectories at 10, 40, 190, and 270 K. Analysis of these functions for each molecule is further decomposed into cage type for each CH 4 and CO 2 guest and cage face for the H 2 O host. CH 4 becomes more constrained, and CO 2 gains freedom in the mixed guest systems. Mixing guest species in gas hydrates alters the intermolecular interaction environment, impacting the rotational motion of the guest molecules. This effect is also seen in the host lattice H 2 O molecules, as evidence of molecular rotations is seen in the CO 2 hydrate at 270 K during longer simulations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Exciton Coherence Length and Dynamics in Graphene Quantum Dot Assemblies
Abstract not provided
Singh, Varun↗
Protein–Mineral Interactions: Molecular Dynamics Simulations Capture Importance of Variations in Mineral Surface Composition and Structure
Molecular dynamics simulations, conventional and metadynamics, were performed to determine the interaction of model protein Gb1 over kaolinite (001), Na+-montmorillonite (001), Ca2+-montmorillonite (001), goethite (100), and Na+-birnessite (001) mineral surfaces. Gb1, a small (56 residue) protein with a well-characterized solution-state nuclear magnetic resonance (NMR) structure and having α-helix, four-fold β-sheet, and hydrophobic core features, is used as a model protein to study protein soil mineral interactions and gain insights on structural changes and potential degradation of protein. From our simulations, we observe little change to the hydrated Gb1 structure over the kaolinite, montmorillonite, and goethite surfaces relative to its solvated structure without these mineral surfaces present. Over the Na+-birnessite basal surface, however, the Gb1 structure is highly disturbed as a result of interaction with this birnessite surface. Unraveling of the Gb1 β-sheet at specific turns and a partial unraveling of the α-helix is observed over birnessite, which suggests specific vulnerable residue sites for oxidation or hydrolysis possibly leading to fragmentation.
Environmental Molecular Sciences Laboratory↗
Origin, nature, and the dynamic behavior of nanoscale vacancy clusters in Ni-rich layered oxide cathodes
Oxygen and transition metal (TM) vacancies in transition metal oxides (TMOs) play an important role in materials functionalities. For TMO cathodes of Li-ion batteries, vacancies can be introduced, either deliberately or inevitably, during materials synthesis or upon electrochemical cycling. Effects of vacancies on the electrochemical properties of cathodes critically depend on the dynamic characteristics of vacancies during the battery cycling. However, a fundamental understanding of such characteristics in the layer structured cathode remains elusive. Here, using scanning transmission electron microscopy (STEM), we reveal a cycling-induced aggregation behavior of oxygen and TM vacancies in a Ni-rich layer structured cathode of secondary particles. We discover that during the initial charging, vacancies aggregate to form visible vacancy-rich nanoregions (vacancy clusters) firstly at the outer layer of the secondary particle, and then the vacancy-cluster populated region extends to the inner part of particle upon fully charged. The spatial distinction of preferential nucleation of vacancy clusters at the outer layer indicates that the vacancy migration and clustering is associated with Li migration flux, which is larger at the outer layer, as all the Li from the inner part has to go through the outer layer. With extended cycling (> 50), these vacancy clusters become immobilized. We further reveal that the generation of these vacancy clusters is correlated to the material-synthesis conditions. Our findings solve a long-standing puzzle on the origin, nature, and behavior of the commonly visible vacancy clusters in the NMC cathode, providing insights on correlation between properties and dynamic behaviors of atomic-scale defects in layered oxide cathodes.
Li-ion batteries, layered oxide cathodes, NMC, cat↗
Interfacial Reactivity and Speciation Emerging from Na-Montmorillonite Interactions with Water and Formic Acid at 200 °C: Insights from Reactive Molecular Dynamics Simulations, Infrared Spectroscopy, and X-ray Scattering Measurements
Reactive organic fluid - mineral interactions at elevated temperatures contribute to the evolution of planetary matter. One of the less studied but important transformations in this regard involves the reactions of formic acid with naturally occurring clays such as sodium montmorillonite. To advance a mechanistic understanding of these interactions, we use ReaxFF reactive molecular dynamics simulations in conjunction with infrared (IR) spectroscopy and X-ray scattering experiments to investigate the speciation behavior of water-formic acid mixtures on sodium montmorillonite interfaces at 473 K and 1 atm. Using a newly developed reactive forcefield, we show that the experimental IR spectra of unreacted and reacted mixture can be accurately reproduced by ReaxFF/MD. We further benchmark the simulation predictions of sodium carbonate and bicarbonate formation in the clay interlayers using Small and Wide-Angle X-ray Scattering measurements. Subsequently, leveraging the benchmarked forcefield, we interrogate the pathway of speciation reactions with emphasis on carbonate, formate, and hydroxide groups elucidating the energetics, transition states, intermediates, and preferred products. Further, we also delineate the differences in reactivities and catalytic effects of clay edges, facets, and interlayers owing to their local chemical environments, which have far reaching consequences in their speciation behavior. The experimental and simulation approaches described in this study and the transferable forcefields can be applied translationally to advance the science of clay-fluid interactions for several applications including subsurface fluid storage and recovery and clay-pollutant dynamics.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Dynamic Interplay between Copper Tetramers and Iron Oxide Boosting CO2 Conversion to Methanol and Hydrocarbons under Mild Conditions.
Atomically precise subnanometer catalysts are of significant interest because of their remarkable efficiency in a variety of catalytic reactions. However, the dynamic changes of active sites under reaction conditions, in particular, the transition of cluster-oxide interface structure have not yet been well-elucidated, lacking in situ measurements. By using multiple state-of-the-art in situ characterizations, here we show a dynamic interplay between copper tetramers and iron oxides in a single-size Cu-4/Fe2O3 catalyst, yielding an enrichment of surface Cu-4-Fe2+ species under reaction conditions that boosts CO2 hydrogenation at near-atmospheric pressures. During reaction, Cu-4 clusters facilitate the reduction of Fe2O3 producing surface-rich Fe2+ species in the proximate sites. The as-formed Fe2+ species in return promotes CO2 activation and transformation over Cu4 cluster, resulting in strikingly high methanol synthesis at low temperatures and C-1/C-3 hydrocarbon production in a high-temperature regime. The discovery of highly active Cu-4-Fe2+ sites thus provides new insights for the atomic-level design of copper catalyst toward high-efficiency CO2 conversion under mild conditions.
CO2 conversion↗