Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamical 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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Solution-state structure and affinities of cyclodextrin: Fentanyl complexes by nuclear magnetic resonance spectroscopy and molecular dynamics simulation

Cyclodextrins (CDs) are investigated for their ability to form inclusion complexes with the analgesic fentanyl and three similar molecules: acetylfentanyl, thiofentanyl, and acetylthiofentanyl. Stoichiometry, binding strength, and complex structure are revealed through nuclear magnetic resonance (NMR) techniques and discussed in terms of molecular dynamics (MD) simulations. It was found that β-cyclodextrin is generally capable of forming the strongest complexes with the fentanyl panel. Two-dimensional NMR data and computational chemical calculations are used to derive solution-state structures of the complexes. Binding of the fentanyls to the CDs occurs at the amide phenyl ring, leaving the majority of the molecule solvated by water, an observation common to all four fentanyls. This finding suggests a universal binding behavior, as the vast majority of previously synthesized fentanyl analogues contain this structural moiety. Furthermore, this baseline study serves as the most complete work on CD:fentanyl complexes to date and provides the insights into strategies for producing future generations of designer cyclodextrins capable of stronger and more selective complexation of fentanyl and its analogues.

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

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi’s golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO 2 under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi’s golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. Finally, we also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

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

Iron surface corrosion in supercritical CO2 at atomic scale investigated by molecular dynamics simulations

Understanding the corrosion behavior of steels in supercritical carbon dioxide (S-CO2) is essential for ensuring the safe application of S-CO2 as a heat-transfer fluid in high-temperature energy systems, including advanced nuclear reactors. In this work, molecular dynamics (MD) simulations using ReaxFF potential are performed to explore the atomic-scale corrosion mechanisms of body-centered cubic iron (BCC-Fe) in S-CO2. The results show that CO2 molecules in S-CO2 decompose at the Fe surface, generating free C and O atoms that form Fe-C and Fe-O bonds and subsequently produce oxides and carbides. Concurrently, Fe atoms dissolve from the surface and diffuse into the S-CO2 region, resulting in interdiffusion of Fe, C and O atoms at the interface. The corrosion-layer thickness calculations show that high pressure and temperature induced by S-CO2 have stronger effects than surface orientation on the corrosion process. In addition, surface Fe atoms undergo substantial displacement under S-CO2 exposure, further accelerating corrosion. When a radiation-induced void is introduced near the Fe surface, the corrosion is enhanced. The void-matrix interface expands the reaction surface area and simultaneously induces corrosion reactions inside the bulk, resulting in a deeper penetration of C and O and thicker corrosion layers. All these results indicate that high-temperature, high-pressure and radiation-induced voids can seriously affect the corrosion of Fe in S-CO2, and must be considered to better use S-CO2 in nuclear facilities.

Li, Wenhua↗
DOE OSTI2018

EXAFS and molecular dynamics simulation studies of Cu-Zr metallic glass: Short-to-medium range order and glass forming ability

Highlights: • Atomic structures of Cu{sub 100−x}Zr{sub x} MGs were studied by EXAFS and MD simulations. • The change of GFA with composition is not dominated by the SRO and MROs. • An apparent reduction of “soft region” at the local high-GFA compositions. - Abstract: Atomic structures of Cu{sub 100−x}Zr{sub x} (x = 41.18–66.67) metallic glasses were investigated by extended X-ray absorption fine structure (EXAFS) spectroscopy and molecular dynamics (MD) simulations. It is demonstrated that both the fraction and average volume of the typical atomic clusters, characterized by Voronoi cells, are of compositional dependence. With the increasing of Zr content, the fraction of Cu and Zr centered Voronoi clusters with coordination number lower than 11 and 14 respectively increases, while that of clusters with larger coordination numbers reduces. Among the major Voronoi clusters, the Cu-centered icosahedra are distinctive. Their average volume is far smaller than that of the other Cu-centered Voronoi cells with the same coordination number. More importantly, they prefer to interlink into pentagon-rich icosahedral super-clusters. However, compared with their neighbors, the glasses with x = 44 and x = 50 exhibit no other particularity in Cu-centered icosahedra than a slightly enhanced interlink at x = 50 although glass forming ability (GFA) reaches a local maximum at these two compositions, indicating that there are other factors dominating the GFA change with composition. A significant reduction in free volume for the loosely packed regions takes place at x = 44 and x = 50, due to which the corresponding Cu-Zr metallic glasses have the maximum microhardness.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Interactions between Ether Phospholipids and Cholesterol as Determined by Scattering and Molecular Dynamics Simulations

Cholesterol and ether lipids are ubiquitous in mammalian cell membranes, and their interactions are crucial in ether lipid mediated cholesterol trafficking. We report on cholesterol s molecular interactions with ether lipids as determined using a combination of small-angle neutron and Xray scattering, and all-atom molecular dynamics (MD) simulations. A scattering density profile model for an ether lipid bilayer was developed using MD simulations, which was then used to simultaneously fit the different experimental scattering data. From analysis of the data the various bilayer structural parameters were obtained. Surface area constrained MD simulations were also performed to reproduce the experimental data. This iterative analysis approach resulted in good agreement between the experimental and simulated form factors. The molecular interactions taking place between cholesterol and ether lipids were then determined from the validated MD simulations. We found that in ether membranes cholesterol primarily hydrogen bonds with the lipid headgroup phosphate oxygen, while in their ester membrane counterparts cholesterol hydrogen bonds with the backbone ester carbonyls. This different mode of interaction between ether lipids and cholesterol induces cholesterol to reside closer to the bilayer surface, dehydrating the headgroup s phosphate moiety. Moreover, the three-dimensional lipid chain spatial density distribution around cholesterol indicates anisotropic chain packing, causing cholesterol to tilt. These insights lend a better understanding of ether lipid-mediated cholesterol trafficking and the roles that the different lipid species have in determining the structural and dynamical properties of membrane associated biomolecules.

EQ-SANS↗
DOE OSTI2022

Predicting initial dissolution rates using structural features from molecular dynamics simulations

Predicting chemical durability of glass materials is important for various applications from daily life such as drink glass and kitchen ware to advanced technologies such as nuclear waste disposal and biomedicine. In this work, we explored prediction of initial dissolution rate through structural features from molecular dynamics (MD) simulations for a wide range of glass compositions (total 28) including borosilicate and aluminosilicate glasses, ZrO 2 -containing and V 2 O 5 -containing boroaluminosilicate glasses. The initial dissolution rates (r 0 ) measured experimentally at 90 °C with varying solution conditions were correlated with structural features (e.g., polyhedral linkages and non-bridging oxygen species) obtained from MD simulations, either from this study or from literature. Since hydrolysis of the glass network through breaking of the network former linkages (e.g., Si-O-Si, Si-O-Al, etc.) is a critical step of network glass dissolution, the statistics of these linkages obtained from MD were correlated to r0 through linear regression, where the coefficient of determination (R 2 ) and root mean square error are found to be 0.949 and 0.681, respectively. This model was compared and discussed with existing models developed by various approaches including machine learning, the kinetic rate equation, topological constraint theory, and other descriptors from MD simulations. The discussion provides insights on future model improvements to predict glass dissolution. In addition, the impact of V 2 O 5 on the glass dissolution was examined in detail, implicating that the impact is not the same across all glass compositions and test conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Crystallographic and Molecular Dynamics Simulation Analysis of Escherichia Coli Dihydroneopterin Aldolase

To understand the structural basis for the biochemical differences and further investigate the catalytic mechanism of DHNA, we have determined the structure of EcDHNA complexed with NP at 1.07-Å resolution [PDB:2O90], built an atomic model of EcDHNA complexed with the substrate DHNP, and performed molecular dynamics (MD) simulation analysis of the substrate complex. EcDHNA has the same fold as SaDHNA and also forms an octamer that consists of two tetramers, but the packing of one tetramer with the other is significantly different between the two enzymes. Furthermore, the structures reveal significant differences in the vicinity of the active site, particularly in the loop that connects strands β3 and β4, mainly due to the substitution of nearby residues. The building of an atomic model of the complex of EcDHNA and the substrate DHNP and the MD simulation of the complex show that some of the hydrogen bonds between the substrate and the enzyme are persistent, whereas others are transient. The substrate binding model and MD simulation provide the molecular basis for the biochemical behaviors of the enzyme, including noncooperative substrate binding, indiscrimination of a pair of epimers as the substrates, proton wire switching during catalysis, and formation of epimerization product.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2600+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6,000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with great emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with a suite of experimental data obtained from x-ray diffraction, infrared spectroscopy, spectroscopic ellipsometry and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stability of Calcium Ion Battery Electrolytes: Predictions from Ab Initio Molecular Dynamics Simulations

Multivalent batteries, such as magnesium-ion, calcium-ion, and zinc-ion batteries, have attracted significant attention as next-generation electrochemical energy storage devices to complement conventional lithium-ion batteries (LIBs). Among them, calcium-ion batteries (CIBs) are the least explored due to the difficult reversible Ca deposition-dissolution. In this work, we examined the stability of four different Ca salts with weakly coordinating anions and three different solvents commonly employed in existing battery technologies to identify suitable candidates for CIBs. By employing Born-Oppenheimer molecular dynamics (BOMD) simulations on salt-Ca and solvent-Ca interfaces, we find that the tetraglyme solvent and carborane salt are promising candidates for CIBs. Due to the strong reducing nature of the calcium surface, the other salts and solvents readily decompose. Further, we explain the microscopic mechanisms of salt/solvent decomposition on the Ca surface using time-dependent projected density of states, time-dependent charge-transfer plots, and climbing-image nudged elastic band calculations. Collectively, this work presents the first mechanistic assessment of the dynamical stability of candidate salts and solvents on a Ca surface using BOMD simulations, and provides a predictive path toward designing stable electrolytes for CIBs.

Born-Oppenheimer molecular dynamics↗
DOE OSTI2020

Nonadiabatic molecular dynamics simulation of C2H22+ in a strong laser field*

We investigate the alignment dependence of the strong laser dissociation dynamics of molecule C 2 H 2 2 + in the frame of real-time and real-space time-dependent density function theory coupled with nonadiabatic quantum molecular dynamics (TDDFT-MD) simulation. This work is based on a recent experiment study “ultrafast electron diffraction imaging of bond breaking in di-ionized acetylene” [Wolter et al , Science 354 , 308–312 (2016)]. Our simulations are in excellent agreement with the experimental data and the analysis confirms that the alignment dependence of the proton dissociation dynamics comes from the electron response of the driving laser pulse. Our results validate the ability of the TDDFT-MD method to reveal the underlying mechanism of experimentally observed and control molecular dissociation dynamics.

Physics↗
DOE OSTI2014

Molecular dynamics simulation of shock induced ejection on fused silica surface

Shock response and surface ejection behaviors of fused silica are studied by using non-equilibrium molecular dynamics combining with the Tersoff potential. First, bulk structure and Hugoniot curves of fused silica are calculated and compared with experimental results. Then, the dynamical process of surface ejection behavior is simulated under different loading velocities ranging from 3.5 to 5.0 km∕s, corresponding to shock wave velocities from 7.1 to 8.8 km∕s. The local atomistic shear strain parameter is used to describe the local plastic deformation under conditions of shock compression or releasing. Our result shows that the shear strain is localized in the bottom area of groove under the shock compression. Surface ejection is observed when the loading velocity exceeds 4.0 km∕s. Meanwhile, the temperature of the micro-jet is ∼5574.7 K, which is close to experiment measurement. Several kinds of structural defects including non-bridging oxygen are found in the bulk area of the sample after ejection.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Adsorption mechanisms of cesium at calcium-silicate-hydrate surfaces using molecular dynamics simulations

We report understanding the mechanisms of interactions and immobilization of radionuclides in cement systems at the molecular level is key to developing safe and novel solutions for the long-term storage of nuclear wastes. The molecular level adsorption mechanisms of cesium (Cs + ) ions onto tobermorite 9 Å, tobermorite 14 Å, and jennite – three crystalline structural analogues for calcium-silicate-hydrates, the main components in cement-based materials, were investigated using molecular dynamics (MD) simulations. Convergence monitoring of the simulations, using the root mean square displacement with average correlation analysis, indicated that MD trajectories of 10–15 ns were needed given the size of the selected computational cell to properly capture the dynamic process of adsorption. Cs + ions adsorbed as inner-sphere complexes at the tobermorite surfaces, while they showed lower affinity for jennite. The charge and structure of the surface layer strongly influenced the adsorption mechanisms of Cs + ions. Strong association with the hexagonal cavities of the silica sites were seen at the tetrahedral SiO 4 surface of the tobermorite structures. Cs + ion adsorption was attributed to co-ion adsorption at the octahedral CaO 6 surface of tobermorite 9 Å, cation exchange at that of tobermorite 14 Å, and electrostatic interaction with the surface hydroxyl groups for jennite. Distribution coefficients of Cs + ions onto cement pastes derived from molecular dynamics by taking the surface structure and dynamics of Cs + ions directly into account showed a range of overlap with experimental results.

42 ENGINEERING↗
DOE OSTI2015

Combining molecular dynamics with mesoscopic Green’s function reaction dynamics simulations

In many reaction-diffusion processes, ranging from biochemical networks, catalysis, to complex self-assembly, the spatial distribution of the reactants and the stochastic character of their interactions are crucial for the macroscopic behavior. The recently developed mesoscopic Green’s Function Reaction Dynamics (GFRD) method enables efficient simulation at the particle level provided the microscopic dynamics can be integrated out. Yet, many processes exhibit non-trivial microscopic dynamics that can qualitatively change the macroscopic behavior, calling for an atomistic, microscopic description. We propose a novel approach that combines GFRD for simulating the system at the mesoscopic scale where particles are far apart, with a microscopic technique such as Langevin dynamics or Molecular Dynamics (MD), for simulating the system at the microscopic scale where reactants are in close proximity. This scheme defines the regions where the particles are close together and simulated with high microscopic resolution and those where they are far apart and simulated with lower mesoscopic resolution, adaptively on the fly. The new multi-scale scheme, called MD-GFRD, is generic and can be used to efficiently simulate reaction-diffusion systems at the particle level.

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

Ballistic-diffusive phonon transport in cellulose nanocrystals by ReaxFF molecular dynamics simulations

Cellulose nanocrystals (CNCs), a ubiquitous nano-sized natural biopolymer, have drawn considerable attention in recent years due to their excellent mechanical and chemical properties. However, thermal properties of CNCs, which are critical for their potential applications as structural or functional materials, are less investigated. In this work, their thermal properties were systematically studied by molecular dynamics (MD) simulations in terms of the effects of polymorphs, size, temperature, and strain. It was found that the thermal conductivities depend on the polymorphs of CNCs, but the thermal conductance of unit chain of each polymorph is similar. Moreover, the strong dependence of the thermal conductivities of CNCs on their length and cross-sectional area indicates the existence of remarkable ballistic-diffusive phonon transport. Single chain CNCs can achieve a thermal conductivity of ~6 W/(m∙K), implying their potential as high thermal conductive polymer. The mechanical strain can enhance the alignment of cellulose skeleton, which reduces phonon scattering and thus increases thermal conductivity. Possible underlying mechanism of the varying thermal conductivity was further discussed and attributed to the shift of phonon frequencies and varied degree of orientation. This study reveals the sophisticated ballistic-diffusive phonon transport in CNCs, and pave path to future design of CNC-based materials with desirable thermal properties.

42 ENGINEERING↗
DOE OSTI2014

Communication: On the consistency of approximate quantum dynamics simulation methods for vibrational spectra in the condensed phase

Including quantum mechanical effects on the dynamics of nuclei in the condensed phase is challenging, because the complexity of exact methods grows exponentially with the number of quantum degrees of freedom. Efforts to circumvent these limitations can be traced down to two approaches: methods that treat a small subset of the degrees of freedom with rigorous quantum mechanics, considering the rest of the system as a static or classical environment, and methods that treat the whole system quantum mechanically, but using approximate dynamics. Here, we perform a systematic comparison between these two philosophies for the description of quantum effects in vibrational spectroscopy, taking the Embedded Local Monomer model and a mixed quantum-classical model as representatives of the first family of methods, and centroid molecular dynamics and thermostatted ring polymer molecular dynamics as examples of the latter. We use as benchmarks D{sub 2}O doped with HOD and pure H{sub 2}O at three distinct thermodynamic state points (ice Ih at 150 K, and the liquid at 300 K and 600 K), modeled with the simple q-TIP4P/F potential energy and dipole moment surfaces. With few exceptions the different techniques yield IR absorption frequencies that are consistent with one another within a few tens of cm{sup −1}. Comparison with classical molecular dynamics demonstrates the importance of nuclear quantum effects up to the highest temperature, and a detailed discussion of the discrepancies between the various methods let us draw some (circumstantial) conclusions about the impact of the very different approximations that underlie them. Such cross validation between radically different approaches could indicate a way forward to further improve the state of the art in simulations of condensed-phase quantum dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

All-Atom Molecular Dynamics Simulations of Cationic Polyelectrolyte Brushes in the Presence of Halide Counterions

Understanding the response of the charged polyelectrolyte (PE) brushes and brush-supported water and ions to the changes in the nature of screening counterions is significant in developing strategies for utilizing PE brushes in various applications. In this paper, we employ all-atom molecular dynamics (MD) simulations for studying the behavior of cationic [poly(2-(methacryloyloxy)ethyl) trimethylammonium] (PMETA) brushes and brush-supported water and ions in the presence of different halide (X: I – , Br – , Cl – , and F – ) screening counterions. We find that despite the F– ion having the largest charge density, the extent of binding of the counterions on the PMETAX brushes varies as I – > Br – > Cl – > F – , leading to PMETAX brush height being least with I– counterions and greatest with F – counterions. This trend in the binding of the halide ions matches the previous experimental result and can be explained by identifying the chaotropic nature of the I – and Br – ions that promote a disruption of water structure and a more favorable binding of the ions to the polymer chains. Such a binding trend also ensures that the order of the water molecules around the PMETAX chains or the counterions as well as the number of water–water hydrogen bonds inside the brush layer increases in the following order of the counterion-specific PMETAX brushes: F – > Cl – > Br – > I – . Furthermore, halide-ion-PMETAX-chain binding takes place via both interchain and intrachain bridging: intrachain bridging dominates for the case of counterions that show enhanced binding (I – and Br – ), while interchain bridging is more favored for the case of counterions that show weakened binding (F – and Cl – ). Also, a greater degree of intrachain bridging leads to greater compressibility and flexibility of the brush layer. Lastly, we show that the mobility of the halide ions follows a nonmonotonic trend with the charge density: the mobility decreases as I – < Br – < Cl – (as their binding to the PMETAX chains varies as I – > Br – > Cl – ), but the mobility of F – ions is in between that of I – and Br – ions. We argue that the strongly attached hydration layer and the ensuing friction inside the brush layer lead to such a reduced mobility of the F – ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗