Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics computer simulation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Variable Precision Computing

This report summarizes the activities and major accomplishments of the Variable Precision Computing Strategic Initiative project. The overarching goal of this project was to initiate and promote a new paradigm in High Performance Computing (HPC) that would fundamentally change how we represent and make use of representations of real numbers in finite precision and that would establish Lawrence Livermore National Laboratory (LLNL) as a leader in precision-related research for HPC. We pursued three integrated and concurrent research thrusts: new, more dynamic data representations to address data motion and capacity limitations; improved mixed precision algorithms to accelerate computation and to leverage new hardware; and new tools to help developers reason about precision and to automate code transformations. Within these thrust areas, we made significant advances in the use of compressed array data types in numerical calculations, in the theoretical justification, in improved performance, and in the prospects of hardware implementation; in the development of metrics to significantly reduce the storage needs of molecular dynamics simulations; in the development and analysis of efficient mixed-precision algorithms for time-dependent partial differential equations and for orthogonal factorization; and in the capability of tools that can analyze sensitivity of computed results to choices in precision and that can make automated code transformations based on such analysis. We established that there are definite opportunities to realize significant improvements in HPC above and beyond the 2x limits of traditional mixed precision using the technologies we have developed. The project produced numerous papers, presentations, and posters in important venues as well as open-source software contributions to the community, and many of the new capabilities have initiated or are being leveraged in ongoing projects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Uncertainty-aware molecular dynamics from Bayesian active learning for phase transformations and thermal transport in SiC

Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.

36 MATERIALS SCIENCE↗
DOE OSTI2022

E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials

Abstract This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.

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

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

ISiCLE: A Quantum Chemistry Pipeline for Establishing in Silico Collision Cross Section Libraries

High throughput, comprehensive, and confident identifications of metabolites and other chemicals in biological and envi-ronmental samples will revolutionize our understanding of the role of these intriguing and chemically diverse molecules in biological systems. Despite recent technological advances, a metabolomics study still results in the detection of a dispropor-tionate number of features than cannot be confidently assigned to a chemical structure. This inadequacy is driven by the sin-gle most significant limitation in metabolomics: the reliance on reference libraries constructed through the analysis of au-thentic reference chemicals and the limited commercial availability of these compounds. To this end, we have developed the in silico chemical library engine (ISiCLE), a high-performance computing-friendly cheminformatics workflow for generating chemical properties. In the instantiation described here, we predict probable three-dimensional molecular conformers (i.e. conformational isomers) using chemical identifiers as input, from which collision cross sections (CCS) are derived. The ap-proach employs state-of-the-art first-principles simulation, distinguished by use of molecular dynamics, quantum chemistry, and ion mobility calculations to generate predictions, all without training data. Importantly, optimization of ISiCLE includ-ed a refactoring of the popular MOBCAL code for trajectory-based mobility calculations, improving its computational effi-ciency by over two orders of magnitude. Calculated CCS values were validated against 1,983 experimental CCS values and compared to previously reported CCS calculation approaches. Average CCS calculation absolute error is 3.2% with standard deviation 3.1%. An online database is introduced for sharing both calculated and experimental CCS values (metabolom-ics.pnnl.gov), initially including a CCS library with over 0.8M entries. Finally, three successful applications of molecule characterization using calculated CCS are described, including the identification of an environmental degradation product, the separation of molecular isomers, and the decoding of complex blinded mixtures of exposure chemicals. This work represents a promising method to address the limitations of small molecule identification, and offers alternatives to standards-based chemical identification amenable to high-performance computation.

metabolomics, standards-free, collision cross sect↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

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

Insights into the Properties of MXenes and MXene Analogs from Atomistic Simulation

We review the properties that have been predicted for MXenes and MXene analogs from computational simulation with a focus on structural and electronic properties, energy storage, and ion transport. Methods considered range from quantum mechanical approaches to classical molecular dynamics. We conclude by reviewing current limitations and outstanding questions for simulation and MXene properties that have been little explored to-date.

Muraleedharan, Murali Gopal↗
DOE OSTI2026

Elastic strain engineering of lattice thermal conductivity of silicon: An ab-initio study

Silicon (Si) is the most essential material in the semiconductor industry. It is important to manage the thermal properties of crystalline Si. Elastic strain engineering (ESE) has proven to be an effective tool in controlling the electrical conductivity of Si in strained-silicon technology; its effects on the thermal conductivity of silicon, therefore, warrants careful investigation. The ESE effect is much more pronounced for nanostructured materials due to the ultralarge elastic strains (on the order of 10%) achievable at the nanoscale. In this work, the lattice thermal conductivity (κ L ) of Si under hydrostatic, biaxial, and uniaxial strain states is studied with ab-initio simulations, and the values of strain-dependent κ L compare well with experimental results and existing molecular dynamics simulations. To understand the mechanisms of strain-modulated κ L , the phonon bands, scattering rate, and Grüneisen parameters of phonon modes are computed. It is shown that strain can significantly change the anharmonicity of the crystal system, thus changing phonon scattering rates and κ L . Our results demonstrate that ESE can reduce silicon κ L by up to approximately 90%. Furthermore, uniaxial and biaxial strains can induce highly anisotropic thermal conductivity in Si, with relative variations up to 58.5% and 14.5%, respectively.

Anisotropy in Thermal Conductivity↗
DOE OSTI2025

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗
DOE OSTI2019

Graeme A. Bird

Few of those who read the 1963 research note by Graeme A. Bird in Physics of Fluids could have imagined that 50 years later the proposed new numerical technique would have become the dominant numerical technique in molecular gas dynamics. The introduction of the Direct Simulation Monte Carlo (DSMC) method not only has altered the field of molecular gas dynamics but also has influenced fields such as physical chemistry, mathematics, computer science, and aerothermodynamics. Further, the DSMC method has been used to probe into previously uninvestigated theoretical aspects of the Boltzmann equation and has also served as a platform for the development of nonequilibrium chemistry models. The DSMC method’s most noteworthy achievement is that molecular gas dynamics became a practical tool in the hands of aerospace engineers in many situations of spacecraft design and mission analysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Horizons of modern molecular dynamics simulation in digitalized solid freeform fabrication with advanced materials

Our ability to shape and finish a component by combined methods of fabrication including (but not limited to) subtractive, additive and/or no theoretical mass-loss/addition during the fabrication is now popularly known as solid free form fabrication. Fabrication of a telescope mirror is a typical example where grinding and polishing processes are first applied to shape the mirror and thereafter an optical coating is usually applied to enhance its optical performance. The area of nanomanufacturing cannot grow without a deep knowledge of the fundamentals of materials and consequently, the use of computer simulations is becoming ubiquitous. This article is intended to introduce the most recent advances in the computation benefit specific to the area of solid free form fabrication as these systems are traversing through the journey of digitalisation and Industry-4.0. Specifically, this article demonstrates that the application of the latest materials modelling approaches, based on techniques such as molecular dynamics, are enabling breakthroughs in applied precision manufacturing techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

N-Doping improves charge transport and morphology in the organic non-fullerene acceptor O-IDTBR

Molecular doping has been shown to improve the performance of various organic (opto)electronic devices. When compared to p-doped systems, research into n-doped organic small-molecules is relatively limited, primarily due to the lack of suitable dopants and the often encountered unfavourable microstructural effects. These factors have prevented the use of n-doping in a wider range of existing materials, such as non-fullerene acceptors (NFAs), that have already shown great promise for a range of (opto)electronic applications. Here, we show that several different molecular n-dopants, namely [1,2-b:2',1'-d]benzo[i][2.5]benzodiazocine potassium triflate adduct (DMBI-BDZC), tetra-n-butylammonium fluoride (TBAF) and 4-(2,3-dihydro-1,3-dimethyl-1H-benzimidazol-2-yl)-N,N-dimethylbenzenamine (N-DMBI), can be used to n-dope the molecular semiconductor O-IDTBR, a promising NFA, and increase the electron field-effect mobility to >1 cm 2 V -1 s -1 . By combining complementary experimental techniques with computer simulations of doping and charge carrier dynamics, we show that improved charge transport arises from synergistic effects of n-type doping and morphological changes. Specifically, a new, previously unreported dopant-induced packing orientation results in one of the highest electron mobility values reported to-date for an NFA molecule. Overall, this work highlights the importance of dopant–semiconductor interactions and their impact on morphology, showing that dopant-induced molecular packing motifs may be generic and a key element of the charge transport enhancement observed in doped organics.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Anisotropic shock response of single-crystalline β-phase tin shock physics, molecular dynamics, plasticity

Mesoscale simulations of the dynamic response of polycrystalline metals to shockwave compression can provide unique insight in to the nature of the various physical mechanisms responsible for material failure. This approach requires a constitutive description for individual grains and boundaries, including defects such as dislocations, within an explicit representation of the microstructure geometry and evolving deformation fields. Computational models of the single-crystal constituents cannot be unambiguously constrained by traditional measurements of the shock or stress-strain response of polycrystalline metals. Instead, these models require comprehensive measurements of the anisotropic shock response of single crystals for their calibration and validation. We present a coordinated experimental and simulation campaign on the shock response of single-crystalline β-phase tin demonstrating a remarkable anisotropic elastic-plastic response of the metal. This anisotropy will be explained with the help of molecular dynamics simulations that show preferred twinning in one orientation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effects of Fe atoms on hardening of a nickel matrix: Nanoindentation experiments and atom-scale numerical modeling

Significant hardening effect due to Fe concentrations in Ni-based alloys with face-centered-cubic structure has been studied by using a combined experimental and atomistic-based computational approach via nanoindentation tests. The obtained experimental load–displacement data for the [0 0 1] crystal orientation reached a qualitative good agreement with molecular dynamics simulations results, leading to strong evidence that the main strengthening factors are associated to sluggish dislocation diffusion, reduced defect sizes and the nucleation of tetrahedral stacking faults. Here, interstitial type prismatic dislocation loops mainly formed by $\frac{1}{6}$$\langle$1 1 2$\rangle$ Shockley dislocations are nucleated during the loading process, where their interaction leads to the formation of pyramidal shaped stacking fault which are mainly created by $\frac{1}{3}$$\langle$1 0 0$\rangle$ Hirth dislocations lines. Observing both types of defects coexisting in the same plastic deformation zone by both approaches. Reported mechanical data, measured experimentally and interpreted numerically, are also in accordance with microstructural SEM and TEM investigations.

36 MATERIALS SCIENCE↗