Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical molecular 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.

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dipolar Aleppo lattice: Ground state ordering and ergodic dynamics in the absence of vertex frustration

We introduce the Aleppo spin ice geometry, another variation of decimated square ice patterns, which in contrast to similar systems previously studied, does not exhibit vertex frustration. Using synchrotron-based photoemission electron microscopy, we directly visualize low-energy states achieved after thermal annealing, in addition to temperature-dependent moment fluctuations. The results reveal the observation of ground state patterns and the absence of ergodicity-breaking dynamics. Our observations further confirm vertex frustration to be an important criterion for the emergence of ergodicity transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quantum simulation: From spin models to gauge-gravity correspondence

Quantum many-body systems with coherent and controllable interactions enable the realization of novel quantum materials, the quantum simulation of problems that cannot be simulated on classical computers, computational systems that are exponentially faster than existing classical algorithms, and the direct investigation of the conjectured duality between gravitational theories and quantum field theories. This project focusses on quantum phase transitions in strongly interacting many-body spin systems, and in particular on the dynamical evolution near the critical point and across the phase transition. The dynamics of a large quantum system containing more than 50 qubits is investigated in various critical regimes. In particular, the Kibble-Zurek transition has been characterized and critical exponents have been measured. Furthermore, quantum many-body scars have been observed and characterized for the first time. We are investigating the Kibble-Zurek mechanism in the quantum phase transition in one dimension from a disordered system to a system with ordering induced by Rydberg interactions between neighboring atoms. The Kibble-Zurek mechanism describes the formation of topological defects as a system evolves across a phase transitions, and is expected to apply to a large number of problems, including cosmology, where it should describe the structure of domain formation in the early universe.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

A new interatomic potential for mixed Mg-Al-Ga-In spinels

While density functional theory (DFT) has become the de facto approach for accurate simulation of materials at the atomic scale, there are many aspects of materials that are simply out of reach of DFT methods. In particular, finite temperature properties such as diffusivities, the structure and properties of grain boundaries and interfaces, and the study of defect properties in complex alloys are computationally challenging for DFT methods. Recently, a new class of spinels in which three cations order over two sublattices was discovered. In order to predict the properties of these types of structures, classical potentials are a must. Here, in this work, we derive a new classical potential for Mg-bearing spinels in which the B cations are Al, Ga, and/or In. The potential does well in describing the DFT energetics of various spinel structures as a function of chemistry and inversion. In particular, it reproduces the thermodynamically favorable MgAlGaO 4 structure while correctly predicting that neither MgAlInO 4 nor MgGaInO 4 are stable. Further, it reproduces physical trends in elastic properties as compared against experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coordination and thermophysical properties of select trivalent lanthanides in LiCl–KCl

The coordination chemistry of various fission and decay products, such as actinides and lanthanides, are crucial to the commercial deployment of molten salt reactors as they can affect the thermophysical properties. In this report we examined the structure, coordination environment, and physical properties such as the density and the vibrational density of states for three lanthanide species, namely Ce, Eu, and Sm in the LiCl–KCl eutectic system using a combination of quantum mechanics simulations and spectroscopic experiments. Quantum mechanics molecular dynamics (QM-MD) modelling was employed to determine the physical properties of each system resulting in accurate local coordination of each species. Then, the vibrational density of states (DOS) was determined using a two-phase thermodynamic modelling which was then compared to the experimentally obtained Raman spectra of the species in molten LiCl–KCl having the eutectic composition. We find that Ce 3+ , Eu 3+ and Sm 3+ all adopt octahedral local coordination environments in the eutectic salt composition in good agreement with experimental results. Ce3+ is found to fluctuate between an octahedral six-coordinated and a seven-coordinated structure due to the increased local proximity of Cl in the eutectic salt, resulting in a lower fluidicity/diffusivity than the other trivalent lanthanides studied. The thermophysical properties of the eutectic composition with trivalent lanthanides were not significantly different from the pure eutectic salt composition, but several changes were noted.

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

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗
NASA NTRS2005

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗
DOE OSTI2021

Two-dimensional phase cartography for high-harmonic spectroscopy

Extremely nonlinear spectroscopy based on high-order-harmonic generation has become a powerful investigation method for attosecond dynamics in gas and solid targets. In particular, the phase of harmonic emission was shown to carry profound insight into atomic and molecular structure and dynamics. However, current techniques offer phase measurements only along specific directions, thus providing partial characterization. Here we report on a new approach combining optical and quantum interferometers measuring along two dimensions the intensity and phase of harmonic emission from aligned molecules in the exact same experimental conditions. This two-dimensional cartography technique measures the phase with no arbitrary offset and no uncertainty on its sign. Measurements along different dimensions can be combined in two ways: either a single mapping or a redundant mapping allowing high-precision phase recovery using a Shack–Hartmann-like algorithm. We demonstrate both methods in a nitrogen test case, which allows disentangling structural and dynamical effects. Two-dimensional phase cartography paves the way to high-resolution high-harmonic spectroscopy for applications such as quantum orbital tomography and attosecond charge migration in molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2012

Ground-Based Centimeter, Millimeter, and Submillimeter Observations of Recent Comets

Comets provide important clues to the physical and chemical processes that occurred during the formation and early evolution of the Solar System, and could also have been important for initiating prebiotic chemistry on the early Earth [I]. Comets are comprised of molecular ices, that may be pristine interstellar remnants of Solar System formation, along with high-temperature crystalline silicate dust that is indicative of a more thermally varied history in the protosolar nebula [2]. Comparing abundances of cometary parent volatiles, and isotopic fractionation ratios, to those found in the interstellar medium, in disks around young stars, and between cometary families, is vital to understanding planetary system formation and the processing history experienced by organic matter in the so-called interstellar-comet connection [3]. In the classical picture, the long-period comets probably formed in the nebular disk across the giant planet formation region (5-40 AU) with the majority of them originating from the Uranus-Neptune region. They were subsequently scattered out to the Oort Cloud (OC) by Jupiter. The short-period comets (also known as ecliptic or Jupiter Family Comets - JFC) reside mainly in the Edgeworth-Kuiper belt where they were formed. Given the gradient in physical conditions expected across this region of the nebula, chemical diversity in this comet population is to be expected [4,5]. We have conducted observations of comets I 03P/Hartley 2 (JFC) and C/2009 PI (Garradd) (OC), at primarily millimeter and submillimeter wavelengths, to determine important cosmogonic quantities, such as the ortho:para ratio and isotope ratios, as well as probe the origin of cometary organics and if they vary between the two dynamic reservoirs.

Milam, S. N.↗
DOE OSTI2021

On the molecular correlations that result in field-dependent conductivities in electrolyte solutions

Employing recent advances in response theory and nonequilibrium ensemble reweighting, we study the dynamic and static correlations that give rise to an electric field-dependent ionic conductivity in electrolyte solutions. We consider solutions modeled with both implicit and explicit solvents, with different dielectric properties, and at multiple concentrations. Implicit solvent models at low concentrations and small dielectric constants exhibit strongly field-dependent conductivities. We compare these results to Onsager-Wilson theory of the Wien effect, which provides a qualitatively consistent prediction at low concentrations and high static dielectric constants but is inconsistent away from these regimes. The origin of the discrepancy is found to be increased ion correlations under these conditions. Explicit solvent effects act to suppress nonlinear responses, yielding a weakly field-dependent conductivity over the range of physically realizable field strengths. By decomposing the relevant time correlation functions, we find that the insensitivity of the conductivity to the field results from the persistent frictional forces on the ions from the solvent. Our findings illustrate the utility of nonequilibrium response theory in rationalizing nonlinear transport behavior.

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

Automated Development of Molten Salt Machine Learning Potentials: Application to LiCl

The in silico modeling of molten salts is critical for emerging "carbon-free" energy applications but is inhibited by the cost of quantum mechanically treating the high polarizabilities of molten salts. Here, we integrate configurational sampling using classical force fields with active learning to automate and accelerate the generation of Gaussian approximation potentials (GAP) for molten salts. This methodology reduces the number of expensive ab initio evaluations required for training set generation to O(100), enabling the facile parametrization of a molten LiCl GAP model that exhibits a 19 000-fold speedup relative to AIMD. The developed molten LiCl GAP model is applied to sample extended spatiotemporal scales, permitting new physical insights into molten LiCl's coordination structure as well as experimentally validated predictions of structures, densities, self-diffusion constants, and ionic conductivities. The developed methodology significantly lowers the barrier to the in silico understanding and design of molten salts across the periodic table.

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

Interplay of coulomb and exciton–phonon coupling controls singlet fission dynamics in two pentacene polymorphs

Pentacene is an important model organic semiconductor in both the singlet exciton fission (SF) and organic electronics communities. We have investigated the effect of changing crystal structure on the SF process, generating multiple triplet excitons from an initial singlet exciton, and subsequent triplet recombination. Unlike for similar organic semiconductors that have strong SF sensitive to polymorphism, we find almost no quantitative difference between the kinetics of triplet pair (TT) formation in the two dominant polymorphs of pentacene. Both pairwise dimer coupling and momentum-space crystal models predict much faster TT formation from the bright singlet excited state of the Bulk vs ThinFilm polymorph, contrasting with the experiment. GW and Bethe–Salpeter equation calculations, including exciton–phonon coupling, reveal that ultrafast phonon-driven transitions in the ThinFilm polymorph compensate the intrinsically slower purely Coulomb-mediated TT formation channel, rationalizing the similarity in observed rates. Taking into account the influence of subtle structural distinctions on both the direct and phonon-mediated SF pathways reveals a predictive capability to these methods, expected to be applicable to a wide variety of molecular crystals.

14 SOLAR ENERGY↗
DOE OSTI2021

Motional Squeezing for Trapped Ion Transport and Separation

Transport, separation, and merging of trapped ion crystals are essential operations for most large-scale quantum computing architectures. In this Letter, we develop a theoretical framework that describes the dynamics of ions in time-varying potentials with a motional squeeze operator, followed by a motional displacement operator. Using this framework, we develop a new, general protocol for trapped ion transport, separation, and merging. We show that motional squeezing can prepare an ion wave packet to enable transfer from the ground state of one trapping potential to another. Here, the framework and protocol are applicable if the potential is harmonic over the extent of the ion wave packets at all times. As illustrations, we discuss two specific operations: changing the strength of the confining potential for a single ion and separating same-species ions with their mutual Coulomb force. Both of these operations are, ideally, free of residual motional excitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Relevance of hydrogen bonded associates to the transport properties and nanoscale dynamics of liquid and supercooled 2-propanol

2-Propanol was investigated, in both the liquid and supercooled states, as a model system to study how hydrogen bonds affect the structural relaxation and the dynamics of mesoscale structures, of approximately several Ångstroms, employing static and quasi-elastic neutron scattering and molecular dynamics simulation. Dynamic neutron scattering measurements were performed over an exchanged wave-vector range encompassing the pre-peak, indicative of the presence of H-bonding associates, and the main peak. The dynamics observed at the pre-peak is associated with the formation and disaggregation of the H-bonded associates and is measured to be at least one order of magnitude slower than the dynamics at the main peak, which is identified as the structural relaxation. Additionally, the measurements indicate that the macroscopic shear viscosity has a similar temperature dependence as the dynamics of the H-bonded associates, which highlights the important role played by these structures, together with the structural relaxation, in defining the macroscopic rheological properties of the system. Importantly, the characteristic relaxation time at the pre-peak follows an Arrhenius temperature dependence whereas at the main peak it exhibits a non-Arrhenius behavior on approaching the supercooled state. The origin of this differing behavior is attributed to an increased structuring of the hydrophobic domains of 2-propanol accommodating a more and more encompassing H-bond network, and a consequent set in of dynamic cooperativity.

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

Predicting 1 H NMR relaxation in Gd 3+ -aqua using molecular dynamics simulations

Atomistic molecular dynamics simulations are used to predict 1 H NMR T 1 relaxation of water from paramagnetic Gd 3+ ions in solution at 25 °C. Simulations of the T 1 relaxivity dispersion function r 1 computed from the Gd 3+ – 1 H dipole–dipole autocorrelation function agree within ≃8% of measurements in the range f 0 ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. Here the simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f 0 ≲ 5 MHz, the simulation overestimates r 1 compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r 1 at high frequencies in chelated Gd 3+ contrast-agents used for clinical MRI.

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

Mechanism and Dynamics of Photodecarboxylation Catalyzed by Lactate Monooxygenase

Photoenzymes are a rare class of biocatalysts that use light to facilitate chemical reactions. Many of these catalysts utilize a flavin cofactor to absorb light, suggesting that other flavoproteins might have latent photochemical functions. Lactate monooxygenase is a flavin-dependent oxidoreductase previously reported to mediate the photodecarboxylation of carboxylates to afford alkylated flavin adducts. While this reaction holds a potential synthetic value, the mechanism and synthetic utility of this process are unknown. Here, we combine femtosecond spectroscopy, site-directed mutagenesis, and a hybrid quantum-classical computational approach to reveal the active site photochemistry and the role the active site amino acid residues play in facilitating this decarboxylation. Light-induced electron transfer from histidine to flavin was revealed, which has not been reported in other proteins. These mechanistic insights enable the development of catalytic oxidative photodecarboxylation of mandelic acid to produce benzaldehyde, a previously unknown reaction for photoenzymes. Furthermore, our findings suggest that a much wider range of enzymes have the potential for photoenzymatic catalysis than has been realized to date.

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

Structure of Novel Phosphonium-Based Ionic Liquids with S and O Substitutions from Experiments and a Mixed Quantum-Classical Approach

This article presents experimental characterization information and synchrotron X-ray scattering measurements on a set of novel O- and S-substituted phosphonium-based ionic liquids (ILs) all coupled with the bis(fluorosulfonyl)imide (FSI-) anion. The ILs include the ethoxyethyltriethylphosphonium (P 222(2O2) + ) and triethyl[2-(ethylthio)ethyl]phosphonium (P 222(2S2) + ) cations, and we contrast results on these with those for unsubstituted triethylpentylphosphonium (P 2225 + ). The article also introduces a physics-based protocol that combines classical force field studies on larger simulation boxes with classical and first-principles studies on smaller boxes. The method produces significantly improved S(q) functions in the regime which in prior publications we have associated with inter- and intraionic adjacency correlations. By understanding which shorter-range structural changes improve S(q) in the q-regime of interest, we are also able to pinpoint specific deficiencies in the classical force field model. The approach we take should be quite general and could help study other complex liquids on different length scales.

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

Molecular Dynamics Simulations of Supercritical Carbon Dioxide and Water using TraPPE and SWM4-NDP Force Fields

The increased levels of carbon dioxide (CO 2 ) emissions due to the combustion of fossil fuels and the consequential impact on global climate change have made CO 2 capture, storage, and utilization a significant area of focus for current research. In most electrochemical CO 2 applications, water is used as a proton donor due to its high availability and mobility and use as a polar solvent. Additionally, supercritical CO 2 is a promising avenue for electrochemical applications due to its unique chemical and physical properties. Consequently, understanding the interactions between water and supercritical CO 2 is of great importance for future electrochemical applications. Molecular dynamics (MD) simulation is a powerful tool that enables atomistic-resolution dynamics of molecular systems, which can complement and guide future experimental investigations. This study employed atomistic MD to study the cosolubilities, codiffusivities, and structure of supercritical CO 2 and water systems, with a polarizable water model (SWM4-NDP) and a nonpolarizable CO 2 model (TraPPE). Additionally, ab initio MD simulations were used to better understand how atomistic polarizable/nonpolarizable models compare to explicit modeling of electron densities. The polarizable water model exhibited substantial improvement in water-associated properties. In conclusion, we anticipate the development of a compatible polarizable CO 2 model to yield similar improvement, providing a pathway for realizing novel high-pressure electrochemical systems.

25 ENERGY STORAGE↗