Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic”

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 451 records · Page 25

Probing Accuracy-Speedup Tradeoff in Machine Learning Surrogates for Molecular Dynamics Simulations

The performance promise of machine learning surrogates of molecular dynamics simulations of soft materials is significant but generally comes at the cost of acquiring large training datasets to learn the complex relationships between input soft material attributes and output properties. Under the constraint of limited high-performance computing resources, optimizing the size of the training datasets becomes paramount. Using an artificial neural network based surrogate for molecular dynamics simulations of confined electrolytes, we explore the tradeoff between surrogate accuracy and computational gains. Accuracy is assessed by computing the root-mean-square errors between the surrogate predictions and the ground truth results obtained via molecular dynamics simulations. The computational performance is judged by evaluating the speedup which incorporates the training dataset creation time. Improvement in accuracy occurs with a loss of speedup, which scales as the inverse of the training dataset size. Furthermore, the link between surrogate generalizability and the accuracy-speedup tradeoff is assessed by examining the errors incurred in surrogate predictions on unseen, interpolated input variables and developing a net speedup metric to capture the associated gains.

Anions↗

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonadiabatic Ab Initio Molecular Dynamics with the Floating Occupation Molecular Orbital-Complete Active Space Configuration Interaction Method [Non-Adiabatic Ab Initio Molecular Dynamics with Floating Occupation Molecular Orbitals CASCI Method]

Here, we show that the floating occupation molecular orbital complete active space configuration interaction (FOMO-CASCI) method is a promising alternative to the widely used complete active space self-consistent field (CASSCF) method in direct nonadiabatic dynamics simulations. We have simulated photodynamics of three archetypal molecules in photodynamics: ethylene, methaniminium cation, and malonaldehyde. We compared the time evolution of electronic populations and reaction mechanisms as revealed by the FOMO-CASCI and CASSCF approaches. Generally, the two approaches provide similar results. Some dynamical differences are observed, but these can be traced back to energetically minor differences in the potential energy surfaces. We suggest that the FOMO-CASCI method represents, due to its efficiency and stability, a promising approach for direct ab initio dynamics in the excited state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of Water Molecules and Ions in Concentrated Lithium Chloride Solutions Probed with Ultrafast 2D IR Spectroscopy

Water and ion dynamics in concentrated LiCl solutions were studied using ultrafast 2D IR spectroscopy with the methyl thiocyanate (MeSCN) CN stretch as the vibrational probe. The IR absorption spectrum of MeSCN has two peaks, one peak for water associated with the nitrogen lone pair of MeSCN (W) and the other peak corresponding to Li + associated with the lone pair (L). To extract the spectral diffusion (structural dynamics) of W and L species, we developed a method that isolates the peak of interest by subtracting the 2D Gaussian proxies of multiple interfering peaks. Center line slope data (normalized frequency–frequency correlation function) for 2D bands from the W and L are fit with triexponential functions. The fastest component (1.1–1.6 ps) is assigned to local hydrogen bond length fluctuations. The intermediate timescale (~4.0 ps) corresponds to the hydrogen bond network rearrangement. The slowest component decays in ~40 ps and corresponds to ion pair and ion cluster dynamics. The very similar W and L spectral diffusion indicates that the motions of the water and ions are strongly coupled. Orientational relaxations of the W and L species were extracted using a new method to eliminate the effects of overlapping peaks. The results show that MeSCN bound to water undergoes orientational relaxation significantly faster than MeSCN bound to Li + . The orientational and spectral diffusion results are compared. In conclusion, a Stark coupling model is used to extract the root mean square average electric field caused by the ion clouds along the CN moiety as a function of concentration.

14 SOLAR ENERGY↗

Molecular Dynamics Simulations of Water Structure and Diffusion in a 1 nm Diameter Silica Nanopore as a Function of Surface Charge and Alkali Metal Counterion Identity

Water confined in nanopores-particularly in pores narrower than 2 nm-displays distinct physicochemical properties that remain incompletely examined despite their importance in nanofluidics, molecular biology, geology, and materials sciences. In this study, we use molecular dynamics simulations to investigate the coordination structure and mobility of water and alkali metals (Li, Na, K, Cs) inside a 1 nm diameter cylindrical silica nanopore as a function of surface charge density, a model system particularly relevant to the alteration kinetics of silicate glasses and minerals in geologic formations. We find that the presence of a negative surface charge and adsorbed counterions within the pore strongly impacts water structure and dynamics. In particular, it significantly orients water O-H bonds toward the surface and slows water diffusion by almost 1 order of magnitude. Ion crowding in the charged nanopore enhances the tendency of counterions to coordinate closely with the silica surface, which moderates the impact of ions on water dynamics. Co-ions are strongly excluded from the nanopore at all surface charges, suggesting that the 1 nm diameter cylindrical silica nanopores likely exhibit nearly ideal semipermeable membrane transport properties.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Exact-Factorization-Based Surface Hopping for Multistate Dynamics

A surface-hopping algorithm recently derived from the exact factorization approach, SHXF, introduces an additional term in the electronic equation of surface hopping that couples electronic states through the quantum momentum. Furthermore, this term not only provides a first-principles description of decoherence, but here we show it is crucial to accurately capture nonadiabatic dynamics when more than two states are occupied at any given time. Using a vibronic coupling model of the uracil cation, we show that the lack of this term in traditional surface-hopping methods, including those with decoherence corrections, leads to failure to predict the dynamics through a three-state intersection, while SHXF performs similarly to the multiconfiguration time-dependent Hartree quantum dynamics benchmark.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Role of Excited State Proton Relays on the Photochemical Dynamics of Water Nanodroplets

In this work, we applied non-adiabatic excited state molecular dynamics in tandem with ab initio electronic structure theory to illustrate a complete mechanistic landscape underpinning the UV absorption initiated photochemical dynamics in water nanodroplets. The goal is to understand the nonequilibrium excited state molecular dynamics initiated by the relaxation of solvated photoelectron, and consequential photochemical processes. The lowest-lying excited state shows the proton dissociation for a single water molecule forming intermediate hydronium complexes through a proton relay. At approximately 100 femtoseconds, the proton relay process gives rise to the relaxation of the excited state accompanied with a rapid increase in the non-adiabatic coupling strength with the ground state, and the nanodroplet non-radiatively decays. The non-adiabatic transition to the ground state produces excited vibrational states that facilitate the recombination of the dissociated proton and hydroxyl group, eventually leading to the desorption of water molecules from the nanodroplet. Additionally, lifetimes of transient photochemical events are also resolved for the relaxation of solvated electron, excited state proton relay, and non-radiative transition.

Stetina, Torin F.↗

Role of Intradomain Heterogeneity on Ion and Polymer Dynamics in Block Polymer Electrolytes: Investigating Interfacial Mobility and Ion-Specific Dynamics and Transport

Block polymers show promise as solid-state battery electrolytes due to the optimization of conductive and mechanical properties enabled via tuning of block chemistry and length. We investigate a polystyrene-block-poly(oligo-oxyethylene methacrylate) (PS-b-POEM) electrolyte doped with various lithium salts to investigate the role of molecular structure on ion transport properties and on local ion dynamics and associations. Anion charge becomes more delocalized with increasing size, reducing the coupling between salt ions while increasing coupling between ion and polymer chain motions and creating a more mobile overall environment. We observe support for this ion-polymer coupling via 1 H, 7 Li and 19 F NMR spectroscopy, from which we obtain ion-specific mobility transition temperatures that differ from the polymer glass transition temperature. We also note faster transport and weaker local energetic interactions with anion size using temperature-dependent NMR diffusometry. 1 H NMR spectroscopy further elucidates polymer chain dynamics and enables quantification of the temperature-dependent fraction of the conducting block that is immobile near the PS-POEM domain interface. As a result, NMR thus represents a species-specific and timescale-specific platform to quantify phase and interface behavior, and to correlate ion-specific transport with polymer chain dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Intra-Domain Heterogeneity on Ion and Polymer Dynamics in Block Polymer Electrolytes: An Approach for Spatially Resolving Dynamics and Ion Transport

The design of safe and high-performance, nanostructured, block polymer (BP) electrolytes for lithium-ion batteries requires a thorough understanding of the key parameters that govern local structure and dynamics. Yet, the interfaces between microphase separated domains can introduce complexities in this local behavior that can be challenging to quantify. Herein, the local polymer, cation (Li + ), and anion dynamics were described in salt-doped polystyrene block poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) through a quantitative framework that considered the effects of polymer architecture, segmental mixing, chain stretching, and confinement on polymer mobility and ion transport. This framework was validated through nuclear magnetic resonance (NMR) spectroscopy measurements on solid (dry) polymer electrolyte samples. Notably, a mobility transition temperature (T mobility ) was identified through NMR spectroscopy that captured the local dynamics more accurately than the thermal glass transition temperature. Additionally, the approach quantitatively described the mobility gradient across a domain when segmental mixing effects were combined with chain stretching and confinement information, especially at higher segregation strengths – facilitating the assessment of local ion diffusion and conductivity. As a result, spatially averaged local ion diffusion predictions quantitatively matched NMR-measured ion diffusivities in the BP samples, while spatially summed ionic conductivity predictions across a domain qualitatively captured trends in the measured ionic conductivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Universal Dynamics of Molecular Reorientation in Hybrid Lead Iodide Perovskites

The role of organic molecular cations in the high-performance perovskite photovoltaic absorbers, methylammonium lead iodide (MAPbI 3 ) and formamidinium lead iodide (FAPbI 3 ), has been an enigmatic subject of great interest. Beyond aiding in the ease of processing of thin films for photovoltaic devices, there have been suggestions that many of the remarkable properties of the halide perovskites can be attributed to the dipolar nature and the dynamic behavior of these cations. Herein, we determine the dynamics of the molecular cations in FAPbI 3 between 4 K and 340 K and the nature of their interaction with the surrounding inorganic cage using a combination of solid state nuclear magnetic resonance and dielectric spectroscopies, neutron scattering, calorimetry, and ab initio calculations. Detailed comparisons with the reported temperature dependence of the dynamics of MAPbI 3 are then carried out which reveal the molecular ions in the two different compounds to exhibit very similar rotation rates (≈8 ps) at room temperature, despite differences in other temperature regimes. For FA, rotation about the N···N axis, which reorients the molecular dipole, is the dominant motion in all phases, with an activation barrier of ≈21 meV in the ambient phase, compared to ≈110 meV for the analogous dipole reorientation of MA. Geometrical frustration of the molecule–cage interaction in FAPbI 3 produces a disordered γ-phase and subsequent glassy freezing at yet lower temperatures. Hydrogen bonds suggested by atom–atom distances from neutron total scattering experiments imply a substantial role for the molecules in directing structure and dictating properties. The temperature dependence of reorientation of the dipolar molecular cations systematically described here can clarify various hypotheses including those of large-polaron charge transport and fugitive electron spin polarization that have been invoked in the context of these unusual materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lattice dynamics and elasticity for ε-plutonium [First-principles lattice dynamics for ε-plutonium]

Here, lattice dynamics and elasticity for the high-temperature ε phase (body-centered cubic; bcc) of plutonium is predicted utilizing first-principles electronic structure coupled with a self-consistent phonon method that takes phonon-phonon interaction and strong anharmonicity into account. These predictions establish the first sensible lattice-dynamics and elasticity data on ε-Pu. The atomic forces required for the phonon scheme are highly accurate and derived from the total energies obtained from relativistic and parameter-free density-functional theory. The results appear reasonable but no data exist to compare with except those from dynamical mean-field theory that suggest ε-plutonium is mechanically unstable. Fundamental knowledge and understanding of the high-temperature bcc phase, that is generally present in all actinide metals before melting, is critically important for a proper interpretation of the phase diagram as well as practical modeling of high-temperature properties.

36 MATERIALS SCIENCE↗

Modulating the dynamics of Brønsted acid sites on PtWOx inverse catalyst

Here metal-metal oxide (M-MO) inverse catalysts are broadly applied. Brønsted acid sites on the oxide overlayers are often hypothesized to drive selective C-O bond activation. However, the Brønsted acid site nature and dynamics under working conditions remain poorly understood due to multiple materials functionalities. Here, we investigate the formation and the dynamics of Brønsted acid and redox sites on PtWO x /C under working conditions. DFT-based thermodynamic calculations and microkinetic modeling reveal a complex interplay between Brønsted acid and redox sites and potentially fast catalyst dynamics at comparable time scales to the chemistry. Combining in situ characterization and probe chemistry, we demonstrate that the density of Brønsted acid sites on the PtWO x /C inverse catalyst could be modulated by up to two orders of magnitude by altering the reaction parameters and by the chemistry itself. We elicit an order of magnitude increase in the acid-catalyzed dehydration average reaction rate by periodic hydrogen pulsing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure, dynamics and stability of water/scCO 2 /mineral interfaces from ab initio molecular dynamics simulations

Here, the interface between a solid and a complex multi-component liquid forms a unique reaction environment whose structure and composition can significantly deviate from either bulk or liquid phase and is poorly understood due the innate difficulty to obtain molecular level information. Feldspar minerals, as typified by the Ca-end member Anorthite, serve as prototypical model systems to assess the reactivity and ion mobility at solid/water-bearing supercritical fluid (WBSF) interfaces due to recent X-ray based measurements that provide information on water-film formation, and cation vacancies at these surfaces. Using density functional theory based molecular dynamics, which allows the evaluation of reactivity and condensed phase dynamics on equal footing, we report on the structure and dynamics of water nucleation and surface aggregation, carbonation and Ca mobilization under geologic carbon sequestration scenarios (T = 323 K and P = 90 bar). We find that water has a strong enthalpic preference for aggregation on a Ca-rich, O-terminated anorthite (001) surface, but entropy strongly hinders the film formation at very low water concentrations. Carbonation reactions readily occur at electron-rich terminal Oxygen sites adjacent to cation vacancies, when in contact with supercritical CO 2 . Cation vacancies of this type can form readily in the presence of a water layer that allows for facile and enthalpicly favorable Ca 2+ extraction and solvation. Apart from providing unprecedented molecular level detail of a complex three component (mineral, water and scCO 2 ) system), this work highlights the ability of modern capabilities of AIMD methods to begin to qualitatively and quantitatively address structure and reactivity at solid-liquid interfaces of high chemical complexity. This work was supported by the US Department of Energy, Office of Fossil Energy (M.-S. L., B. P. M. and V.-A. G.) and the Office of Basic Energy Science, Division of Chemical Sciences, Geosciences and Biosciences (R.R.), and performed at the Pacific Northwest National Laboratory (PNNL). PNNL is a multi-program national laboratory operated for DOE by Battelle. Computational resources were provided by PNNL’s Platform for Institutional Computing (PIC), the W. R. Wiley Environmental Molecular Science Laboratory (EMSL), a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research located at PNNL and the National Energy Research Scientific Computing Center (NERSC) at Lawrence Berkeley National Laboratory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupled wave-packets for non-adiabatic molecular dynamics: a generalization of Gaussian wave-packet dynamics to multiple potential energy surfaces

Accurate simulation of the non-adiabatic dynamics of molecules in excited electronic states is key to understanding molecular photo-physical processes. Here we present a novel method, based on a semiclassical approximation, that is as efficient as the commonly used mean field Ehrenfest or ad hoc surface hopping methods and properly accounts for interference and decoherence effects. This novel method is an extension of Heller's thawed Gaussian wave-packet dynamics that includes coupling between potential energy surfaces. By studying several standard test problems we demonstrate that the accuracy of the method can be systematically improved while maintaining high efficiency. The method is suitable for investigating the role of quantum coherence in the non-adiabatic dynamics of many-atom molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Water dynamics in hydrated amorphous materials: a molecular dynamics study of the effects of dehydration in amorphous calcium carbonate

Amorphous calcium carbonate (ACC) is often a critical transient phase in the formation of crystalline phases of CaCO 3 via dehydration of hydrated ACC. The behavior of the water molecules plays a pivotal role in this transformation. We report here the dynamics of water molecules in ACC at hydration levels from CaCO 3 ·1H 2 O to CaCO 3 ·0.25H 2 O using molecular dynamics (MD) simulations. Due to the presence of highly hydrophilic Ca 2+ and the strong H-bond acceptor CO 3 2- , most of the water molecules in our simulations have restricted translational and orientational dynamics. These are referred here as ‘slow water’. However, a small fraction of them show high diffusivity at all hydration states and are referred here as ‘fast water’. The computed diffusion coefficients of the slow waters, as extrapolated from simulation results, yield diffusion distances of the order of μm on the 1 hour time scale, consistent with rapid dehydration of ACC nano-particles in the laboratory. We correlate our fast waters with the water molecules in the percolating water-rich, Ca 2+ -deficient nano-pores in the structure of Goodwin et al. [Chem. Mater., 2010, 22, 3197], obtained by analysis of X-ray scattering data, and our slow waters with those in the Ca 2+ -rich volumes with less water in their model. The fast waters can be considered to be free rotors on the ns time scale and have orders of magnitude shorter rotational relaxation times than the slow waters.

Saharay, Moumita↗

Quantized Hamiltonian dynamics captures the low-temperature regime of charge transport in molecular crystals

The quantized Hamiltonian dynamics (QHD) theory provides a hierarchy of approximations to quantum dynamics in the Heisenberg representation. We apply the first-order QHD to study charge transport in molecular crystals and find that the obtained equations of motion coincide with the Ehrenfest theory, which is the most widely used mixed quantum-classical approach. Quantum initial conditions required for the QHD variables make the dynamics surpass Ehrenfest. Most importantly, the first-order QHD already captures the low-temperature regime of charge transport, as observed experimentally. We expect that simple extensions to higher-order QHDs can efficiently represent other quantum effects, such as phonon zero-point energy and loss of coherence in the electronic subsystem caused by phonons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamics of skyrmions in chiral magnets: Dynamic phase transitions and equation of motion

We study the dynamics of skyrmions in a metallic chiral magnet. First, we show that skyrmions can be created dynamically by destabilizing the ferromagnetic background state through a spin polarized current. We then treat skyrmions as rigid particles and derive the corresponding equation of motion. The dynamics of skyrmions is dominated by the Magnus force, which accounts for the weak pinning of skyrmions observed in experiments. Finally, we discuss the quantum motion of skyrmions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗