Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical molecular dynamic simulations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Molecular Simulations of Ultrafast Radiation Induced Melting at Metal-Semiconductor Interfaces

Understanding radiation induced ultrafast melting at material interfaces is essential in designing robust electronic devices for aviation/space applications and in laser machining. While it is difficult to achieve the spatial and temporal resolution required to quantify the phenomenon experimentally, simulations can provide the detailed mechanisms of the structural changes that happen during phase transition. In this work, we use molecular simulations to study the effect of radiation damage on silicon carbide (SiC) - tungsten (W) interfaces which is of interest in high power electronics. A multi-scale approach is involved wherein the reactions at the interfaces are quantified using ab-initio molecular dynamics (MD) simulations and classical MD simulations are employed to understand the structural and diffusional changes across the material interface. Finally, coarse-grained Lennard-Jones type models are used to study the larger scale mechanisms and structures obtained due to the induced damages. We show that the response of the material to radiation damage depends on factors such as energy of the incident radiation, thermal properties, and molecular structure of the material.

Ravichandran, Ashwin↗
DOE OSTI2024

Coordinate versus momentum cuts and effects of collective flow on critical fluctuations

We analyze particle number fluctuations in the crossover region near the critical endpoint of a first-order phase transition by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our previous study [V. A. Kuznietsov , ] by incorporating longitudinal collective flow. The scaled variance of particle number distribution inside different coordinate and momentum space acceptances is computed through ensemble averaging and found to agree with earlier results obtained using time averaging, validating the ergodic hypothesis for fluctuation observables. Presence of a sizable collective flow is found to be essential for observing large fluctuations from the critical point in momentum space acceptances. We discuss our findings in the context of heavy-ion collisions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Theoretical upper limits of the thermal conductivity of Si3N4

Silicon nitride (Si3N4) is a promising substrate for high-power electronics due to its superior mechanical properties and potential outstanding thermal conductivity (κ). As experiments keep pushing the upper limit of κ of Si3N4, it is believed that it can reach 450 W/mK, similar to SiC, based on classical models and molecular dynamics simulations. In this work, we reveal from first principles that the theoretical κ upper limits of β-Si3N4 are only 169 and 57 W/mK along the c and a axes at room temperature, respectively. Those of α-Si3N4 are about 116 and 87 W/mK, respectively. The predicted temperature-dependent κ matches well with the highest available experimental data, which supports the accuracy of our calculations, and suggests that the κ upper limit of Si3N4 has already been reached in the experiment. Compared to other promising semiconductors (e.g., SiC, AlN, and GaN), Si3N4 has a much lower κ than expected even though the chemical bonding and mechanical strengths are close or even stronger. We find the underlying reason is that Si3N4 has much lower phonon lifetimes and mean free paths (<0.5 μm) due to the larger three-phonon scattering phase space and stronger anharmonicity. Interestingly, we find that the larger unit cell (with more basis atoms) that leads to a smaller fraction of acoustic phonons is not the reason for lower κ. Grain size-dependent κ indicates that the grain boundary scattering plays a negligible role in most experimental samples. This work clarifies the theoretical κ upper limits of Si3N4 and can guide experimental research.

Physics↗
DOE OSTI2025

Water under hydrophobic confinement: entropy and diffusion

The properties of liquid water are known to change drastically in confined geometries. A most interesting and intriguing phenomenon is that the diffusion of water is found to be strongly enhanced by the proximity of a hydrophobic confining wall relative to the bulk diffusion. We report a molecular dynamics simulation using a classical water model investigating the water diffusion near a non-interacting smooth confining wall, which is assumed to imitate a hydrophobic surface, revealing a pronounced diffusion enhancement within several water layers adjacent to the wall. We present evidence that the observed diffusion enhancement can be accounted for, with a quantitative accuracy, using the universal scaling law for liquid diffusion that relates the diffusion rate to the excess entropy. These results show that the scaling law, which has so far only been used for the description of the diffusion in simple liquids, can successfully describe the diffusion in water. It is shown that the law can be used for the analysis of water dynamics under nanoscale hydrophobic confinement, which is currently a subject of intense research activity.

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

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2008

State of the Art Assessment of Simulation in Advanced Materials Development

Advances in both the underlying theory and in the practical implementation of molecular modeling techniques have increased their value in the advanced materials development process. The objective is to accelerate the maturation of emerging materials by tightly integrating modeling with the other critical processes: synthesis, processing, and characterization. The aims of this report are to summarize the state of the art of existing modeling tools and to highlight a number of areas in which additional development is required. In an effort to maintain focus and limit length, this survey is restricted to classical simulation techniques including molecular dynamics and Monte Carlo simulations.

Wise, Kristopher E.↗
DOE OSTI2021

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

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

Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature

It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. Here, we find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.

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

Modeling of Peptides with Classical and Novel Machine Learning Force Fields: A Comparison

The replacement of classical force fields (FFs) with novel neural-network-based frameworks is an emergent topic in molecular dynamics (MD) simulations. In contrast to classical FFs, which have proven their capability to provide insights into complex soft matter systems at an atomistic resolution, the machine learning (ML) potentials have yet to demonstrate their applicability for soft materials. However, the underlying philosophy, which is learning the energy of an atom in its surrounding chemical environment, makes this approach a promising tool. In particular for the exploration of novel chemical compounds, which have not been considered in the original parametrization of classical FFs. In this article, we study the performance of the ANI-2x ML model and compare the results with those of two classical FFs, namely, CHARMM27 and the GROMOS96 43a1 FF. We explore the performance of these FFs for bulk water and two model peptides, trialanine and a 9-mer of the α-aminoisobutyric acid, in vacuum and water. The results for water describe a highly ordered water structure, with a structure similar to those using ab initio molecular dynamics simulations. The energy landscape of the peptides described by Ramachandran maps show secondary structure basins similar to those of the classical FFs but differ in the position and relative stability of the basins. Details of the sampled structures show a divergent performance of the different models, which can be related either to the short-ranged nature of the ML potentials or to shortcomings of the underlying data set used for training. These findings highlight the current state of the applicability of ANI-2x ML potential for MD simulations of soft matter systems. Simultaneously, they provide insights for future improvements of current ML potentials.

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

Structural and Dynamic Heterogeneity of Deep Eutectic Solvents Composed of Choline Chloride and Ortho-Phenol Derivatives

Structural, thermal, and dynamic properties of four deep eutectic solvents comprising choline chloride paired with ortho-phenolic derivative hydrogen-bond donors were probed using experiments and molecular simulations. The hydrogen-bond donors include phenol, catechol, o-chlorophenol, and o-cresol, in a 3:1 mixture with the hydrogen-bond acceptor choline chloride. Density, viscosity, and pulsed-field gradient NMR diffusivity measurements were conducted over a range of temperatures. Classical and ab initio molecular dynamics simulation results match experimental data reasonably well. Furthermore, the simulation results were then used to perform a more detailed analysis of the local structure and dynamics of these systems.

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

Variable Precision Computing (VPC) (Final Report)

The overall goal of the Variable Precision Computing project is exploring and identifying ways to reduce internodal communication cost in production physics simulation experiments running on a distributed-memory systems by compressing data, either in a lossy or lossless way, thus resulting in a significant computational performance increase. Among the several approaches proposed, the one discussed here has been proposed, designed and implemented by me. It concerns the application of information theory-derived metrics to classical molecular dynamics and computational fluid dynamics simulation experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Visualizing the Three-Dimensional Arrangement of Hydrogen Atoms in Organic Molecules by Coulomb Explosion Imaging

Structure-sensitive methods based on femtosecond light or electron pulses are now making it possible to measure how molecular structures change during light-induced processes. Despite significant progress, high-fidelity imaging of nuclear positions remains a challenge even for relatively small molecular systems and, notably, regarding the positions of hydrogen atoms. As demonstrated in recent work, X-ray-induced Coulomb explosion imaging (CEI) may overcome this obstacle, as its sensitivity does not depend on the mass of the imaged atoms. The photoinduced ring opening of the heterocyclic molecule 2(5 H )-thiophenone has attracted recent interest. Here, in this work, we show that CEI offers a powerful route to imaging the peripheral H atoms in this molecule and thus, more generally, to tracking detailed nuclear motions (e.g., isomerizations) in organic molecules on ultrafast time scales. Specifically, we record momentum-space Coulomb explosion images that report on the three-dimensional positioning of all nuclei within the molecule, for instance, distinguishing H atoms in C–H bonds that lie within or are directed out of the plane defined by the heavy atoms. The prospect of imaging peripheral H atoms to probe photochemical dynamics is explored by coupling ab initio molecular dynamics with classical Coulomb explosion simulations, thereby differentiating potential photoproduct isomers, including those whose structures primarily differ in the position of the hydrogens.

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

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

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

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

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

Long-Chain Lipids Facilitate Insertion of Large Nanoparticles into Membranes of Small Unilamellar Vesicles

Insertion of hydrophobic nanoparticles into phospholipid bilayers is limited to small particles that can incorporate into a hydrophobic membrane core between two lipid leaflets. Incorporation of nanoparticles above this size limit requires the development of challenging surface engineering methodologies. In principle, increasing the long-chain lipid component in the lipid mixture should facilitate incorporation of larger nanoparticles. Here, we explore the effect of incorporating very long phospholipids (C24:1) into small unilamellar vesicles on the membrane insertion efficiency of hydrophobic nanoparticles that are 5–11 nm in diameter. To this end, we improve an existing vesicle preparation protocol and utilized cryogenic electron microscopy imaging to examine the mode of interaction and evaluate the insertion efficiency of membrane-inserted nanoparticles. We also perform classical coarse-grained molecular dynamics simulations to identify changes in lipid membrane structural properties that may increase insertion efficiency. Furthermore, our results indicate that long-chain lipids increase the insertion efficiency by preferentially accumulating near membrane-inserted nanoparticles to reduce the thermodynamically unfavorable disruption of the membrane.

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

Examining normal modes as fundamental heat carriers in amorphous solids: The case of amorphous silicon

Normal mode decomposition of atomic vibrations has been used to provide microscopic understanding of thermal transport in amorphous solids for decades. In normal mode methods, it is naturally assumed that atoms vibrate around their equilibrium positions, and that individual normal modes are the fundamental vibrational excitations transporting heat. With the abundance of predictions from normal mode methods and experimental measurements now available, we carefully analyze these calculations in amorphous silicon, a model amorphous solid. We find a number of discrepancies, suggesting that treating individual normal modes as fundamental heat carriers may not be accurate in amorphous solids. Furthermore, our classical and ab initio molecular dynamics simulations of amorphous silicon demonstrate a large degree of atomic diffusion, especially at high temperatures, leading to the conclusion that thermal transport in amorphous solids could be better described starting from the perspectives of liquid physics rather than from crystalline solids.

36 MATERIALS SCIENCE↗