Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum

The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl - , Na + , and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl - .

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

Molecular dynamics on quantum annealers

Abstract In this work we demonstrate a practical prospect of using quantum annealers for simulation of molecular dynamics. A methodology developed for this goal, dubbed Quantum Differential Equations (QDE), is applied to propagate classical trajectories for the vibration of the hydrogen molecule in several regimes: nearly harmonic, highly anharmonic, and dissociative motion. The results obtained using the D-Wave 2000Q quantum annealer are all consistent and quickly converge to the analytical reference solution. Several alternative strategies for such calculations are explored and it was found that the most accurate results and the best efficiency are obtained by combining the quantum annealer with classical post-processing (greedy algorithm). Importantly, the QDE framework developed here is entirely general and can be applied to solve any system of first-order ordinary nonlinear differential equations using a quantum annealer.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field

We have implemented the Martini force field within Lawrence Livermore National Laboratory’s molecular dynamics program, ddcMD. The program is extended to a heterogeneous programming model so that it can exploit graphics processing unit (GPU) accelerators. In addition to the Martini force field being ported to the GPU, the entire integration step, including thermostat, barostat, and constraint solver, is ported as well, which speeds up the simulations to 278-fold using one GPU vs one central processing unit (CPU) core. A benchmark study is performed with several test cases, comparing ddcMD and GROMACS Martini simulations. The average performance of ddcMD for a protein–lipid simulation system of 136k particles achieves 1.04 µs/day on one NVIDIA V100 GPU and aggregates 6.19 µs/day on one Summit node with six GPUs. The GPU implementation in ddcMD offloads all computations to the GPU and only requires one CPU core per simulation to manage the inputs and outputs, freeing up remaining CPU resources on the compute node for alternative tasks often required in complex simulation campaigns.

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

Spin-orbit coupling induced demagnetization in Ni: Ab initio nonadiabatic molecular dynamics perspective

Spin-orbit coupling (SOC), which can induce spin flip during the relaxation of photoexcited charge carrier, plays a crucial role in spin dynamics. Here in this work, we have used time-domain ab initio nonadiabatic molecular dynamics (NAMD) method to study the SOC induced ultrafast demagnetization in Ni at 300K. The spin-diabatic representation using spin-polarized Kohn-Sham (KS) basis sets and spin-adiabatic representation using spinor basis sets have been applied, and both of them achieve demagnetization in Ni with a timescale around 100fs. The spin-diabatic representation suggests a picture that the electron-phonon coupling (EPC) provides direct energy relaxation channel among the same-spin states, while the SOC can induce spin flip. After photoexcitation, it is found the spin-minority electrons relax to the same-spin states rather than the opposite-spin states, since EPC is larger than SOC by one order of magnitude. By contrast, for the spin-majority electrons, spin flip occurs since there are no empty same-spin states as electron acceptor above the Fermi level. The different relaxation pathways for spin-majority and spin-minority electrons induce the demagnetization. The spin-adiabatic representation provides an Elliott-Yafet spin-phonon scattering picture. The SOC induced reduction of magnetic moment in Ni may induce magnon to drive further demagnetization. The ab initio NAMD simulation provides a critical angle to understand how the SOC and EPC affect demagnetization process in Ni.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
NASA NTRS1997

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗
DOE OSTI2024

Data and Code for Atomic Scale Etching of Diamond: Insights from Molecular Dynamics Simulations

This work investigates the effects of argon ions, hydrogen atoms, and hydrogen ions on the diamond (100) surface using classical molecular dynamics simulations. The purpose of this investigation was to asses plasma processing techniques for applications in quantum device manufacturing. The simulations suggest that combining argon ion smoothing with selective, near threshold energy H removal of amorphous C could be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices. Results were found to differ significantly with interatomic potential, and an analysis of these differences was also carried out. Included in this repository are LAMMPS source files, input scripts, and plotting scripts required to reproduce the data. Also included are the output data required to make all the plots included in the associated publication.

Brenner↗
DOE OSTI2022

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

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

TeraChem: A graphical processing unit-accelerated electronic structure package for large-scale ab initio molecular dynamics

TeraChem was born in 2008 with the goal of providing fast on-the-fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low-cost parallel computer architecture that became more accessible for general-purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Here, electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large-scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large-scale nonadiabatic dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Molecular dynamics study of interstitial He clusters in nickel

This study presents a molecular dynamics analysis focusing on the behavior of interstitial helium (He) clusters in nickel (Ni), examining their formation, stability, and migration energetics. Consistent with previous research, we found that the binding energies of interstitial helium within a helium cluster are positive and increase with the cluster size, indicating a preference for helium atoms to cluster together. However, our findings also reveal that while the formation energy increases monotonically with cluster size, the increase in binding energy is non-monotonic. Importantly, small He clusters were observed to be thermally unstable at reactor operational temperatures (approximately 600 K), with the He 2 cluster exhibiting instability even at room temperature. With a binding energy of 0.49 eV for a He 4 cluster, we hypothesize that for helium bubbles to form via homogeneous nucleation (i.e., through trap mutation) at reactor operating temperatures, the helium concentration must be high enough to facilitate the formation of helium clusters of at least size 4 or larger. As expected, interstitial helium and small helium clusters are highly mobile. This mobility was observed not only at room temperature but also at temperatures as low as approximately 200 K. Furthermore, the mean squared displacement method has been utilized to determine the migration barriers and the corresponding prefactors for clusters ranging from He 1 to He 6

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Exploration of a Novel Technique for Waste Heat Recovery Through Molecular Dynamics: Influence of Wettability and Electric Field on Water and Water-Based Nanofluids

Most of the energy produced globally comes by way of a heat engine. The Carnot principle places a limit as to how thermodynamically efficient a heat engine can be. There is no heat engine that can be 100% thermodynamically efficient and as such a substantial proportion of all heat supplied to a heat engine is lost as waste heat. Waste heat therefore is a large energy source ready to be properly utilized. Herein, a novel approach for converting waste heat to electricity is discussed. It involves the use of the liquid to vapor phase change of a material dielectric (water) or electrolyte (nanofluid) in the embodiment of a capacitor for direct thermal to electrostatic energy conversion. While this method of waste heat recovery could potentially be added to the ever expanding portfolio of energy conversion techniques, a number of aspects must be addressed before it can be brought into practice. Water was seen as an ideal dielectric phase change material given its high relative permittivity ratio when in the liquid form as compared to its vapor form. However, given its short voltage holdoff time the phase change of water would need to occur rapidly. This brings up concerns of explosive boiling. Herein, molecular dynamics analysis into the explosive boiling behavior of thin water films gave more insight into how the interaction between the surface and liquid affected explosive boiling onset time. A Lennard-Jones potential with one interaction site and a Morse potential with three interaction sites between water and solid substrate were used. It was found generally that a stronger interaction between water film and substrate led to faster explosive boiling onset times but an increase in the number of interaction sites delayed explosive boiling, even at the same wettability (contact angle). Understanding changes in the density and enthalpy of vaporization of a liquid dielectric such as water in the presence of an electric field is of importance due to the electrostatic nature of the waste heat conversion method under consideration. Specifically, if both density and enthalpy of vaporization are increased, the thermodynamic efficiency of the waste heat conversion method under consideration is decreased. Electric field effects are explored herein via molecular dynamics using two water models, the TIP4P-Ew and SWM4-NDP. The SWM4-NDP model is polarizable while the TIP4P-Ew model is not, which allows for a determination of the importance of model polarizability (i.e. variation in water model dipole moment) on these two properties of water when subjected to an electric field. Herein it was found that both water models respond similarly in terms of density and vaporization enthalpy variance upon the introduction of an electric field. Comparison was also made to the pressure induced by the electric field (electrostriction pressure) by way of a density comparison and it was found that the predicted electrostriction pressure overestimates the pressure experienced by water. Water by itself has a high enthalpy of vaporization, which limits the efficiency of the newly proposed conversion method. Research both experimental and through simulation has shown that the vaporization enthalpy of nanofluids can be engineered via nanoparticle size and material selection. An avenue less explored is manipulating the enthalpy of vaporization by altering the interaction strength between the nanoparticles and the base fluid. In practice this could be achieved through the addition of coatings to the nanoparticles to alter their wettability to the base fluid. This was explored by using a Lennard-Jones potential and Morse potential to model the interaction between base fluid (water) and the nanoparticle. For nanoparticles 2nm in diameter and at weight percentages up to 6%, the change in vaporization enthalpy due to alterations of the interaction strength between nanoparticle and base fluid was not significant (less than a 1% difference) when compared to the effect of altering the weight percentage of nanoparticles in the nanofluid or introducing an electric field. However, the effect of wettability may still become important at other nanoparticle concentrations and sizes. In all, the studies presented here further the understanding of phase change and thermodynamic properties of water and water based nanofluids under an electrostatic field which will help inform the development of a novel approach to waste heat conversion. The reduction of waste heat will improve energy sustainability outlooks.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Discussion on molecular dynamics (MD) simulations of the asphalt materials

The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effects of Surface Orientation and Termination Plane on Glass‐to‐Crystal Transformation of Lithium Disilicate by Molecular Dynamics Simulations

Glass‐to‐crystal transformation of lithium disilicate is studied using molecular dynamics simulations using an effective partial charge potential. The structural evolution of the interface between glassy and crystalline lithium disilicate is analyzed to simulate crystallization of glass on pre‐existing crystal seeds. Besides previously used atomic number density, the distribution of Q n species (Si tetrahedra with n bridging oxygen) is shown to be an effective parameter for following this transformation quantitatively. The early stages of crystal growth are significantly affected by the orientation and termination of the surface of adjacent crystal, as indicated by calculated atomic density, partial ordering, atomic segregation, and an increase in Q 3 concentration. In particular, under‐coordinated Si within the outer crystal layer is found to be most effective in transforming the amorphous structure toward crystallinity. The increase in Q 3 in the glass close to interface region most clearly shows the initial stage of lithium disilicate crystal growth.

Sun, Wei↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1981

Cluster adsorption on amorphous and crystalline surfaces - A molecular dynamics study of model Pt on Cu and model Pd on Pt

Molecular dynamics was used to study the structure, dispersion and short-time behavior of ten-atom clusters adsorbed onto amorphous and crystalline substrates, in which the cluster atoms differed from the substrate atoms. Two adatom-substrate model systems were chosen; one, in which the interaction energy between adatom pairs was greater than that between substrate pairs, and the other, in which the reverse was true. At relatively low temperature ranges, increased dispersion of cluster atoms occurred: (a) on the amorphous substrate as compared to the FCC(100) surface, (b) with increasing reduced temperature, and (c) with adatom-substrate interaction energy stronger than adatom-adatom interaction. Two-dimensional clusters (rafts) on the FCC(100) surface displayed migration of edge atoms only, indicating a mechanism for the cluster rotation and shape changes found in experimental studies.

Garofalini, S. H.↗