Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Design of high polarization low switching barrier hybrid improper ferroelectric perovskite oxide superlattices

Hybrid improper ferroelectricity is a useful tool to design ABO 3 /A'BO 3 polar superlattices from non polar building blocks. In this study, we have designed high polarization-low switching barrier hybrid improper ferroelectric superlattices with efficient polarization, and polarization-magnetization switching properties above room temperature, using density functional theory and ab initio molecular dynamics simulations. Superlattices with a chemical formula of (AAlO 3 ) m /(A'AlO 3 ) n , where m/n = 1/1, 1/3, 3/1, 1/5 and 5/1, A, A' = Lanthanide and Y cations are considered to outline the design principles behind polarization switching and (LaFeO 3 ) 3 /(CeFeO 3 ) 1 is investigated for polarization-magnetization switching. We find that the unconventional switching paths via out-of-phase rotation Q R- (a 0 a 0 c - ) and tilt precession QTP always yield lower switching barrier compared to those via in-phase rotation Q R+ (a 0 a 0 c + ) and tilt QT (a - a - c 0 ) of BO 6 octahedra. Results from ab initio molecular dynamics simulations estimate the temperature at which the lowest energy barrier can be overcome. Furthermore, it is possible to tune the polarization switching barrier by tuning the tolerance factor, A,A' cation radius mismatch and super lattice periodicity. For switching via Q R- , the switching barrier varies exponentially with rotation angle, indicating how high switching barrier is expected for systems, away from cubic symmetry. We provide a recipe to overcome such a bottleneck by tuning superlattice periodicity. Finally, we have proposed the multiferroic device application concept through a proposed polarization-temperature hysteresis loop and magnetization switching.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anomalous entropy-driven kinetics of dislocation nucleation

The kinetics of dislocation reactions, such as dislocation multiplication, controls the plastic deformation in crystals beyond their elastic limit, therefore critical mechanisms in a number of applications in materials science. We present a series of large-scale molecular dynamics simulations that shows that one such type of reactions, the nucleation of dislocation at free surfaces, exhibit unconventional kinetics, including unexpectedly large nucleation rates under compression, very strong entropic stabilization under tension, as well as strong non-Arrhenius behavior. These unusual kinetics are quantitatively rationalized using a variational transition state theory approach coupled with an efficient numerical scheme for the estimation of vibrational entropy changes. These results highlight the need for a variational treatment of the kinetics to quantitatively capture dislocation reaction kinetics, especially at low-to-moderate strains where large deformations are required to activate reactions. These observations suggest possible explanations to previously observed unconventional deformation kinetics in both molecular dynamics simulations and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Dataset_for_Conserved_macromolecular_architecture_of_Poplar_secondary_cell_walls_revealed_by_ssNMR_and_atomistic_modeling

This dataset contains solid-state 13C NMR data and atomistic molecular dynamics simulation files supporting the study of nanoscale secondary cell wall architecture across 13 genetically diverse Populus trichocarpa genotypes grown under uniform greenhouse conditions in 13C-enriched CO2 atmospheres (~89% 13C enrichment).The dataset contains two collections of solid-state 13C NMR data. (1) 200 MHz data (Bruker Avance III HD, 4 mm HX probe, 10 kHz MAS): raw Bruker TopSpin experiment folders and DMFIT-exported ascii spectra for selective and non-selective 1D 13C-13C spin diffusion experiments (3000 ms mixing) used to quantify inter-polymer spatial proximities, and short-mixing (1 ms) reference spectra used for polymeric abundance quantification by spectral deconvolution. (2) 600 MHz data (Bruker Avance III, 1.6 mm PhoenixNMR HXY probe, 30 kHz MAS): raw Bruker TopSpin experiment folders containing 2D CORD, 2D CP-INADEQUATE, and 13C/1H relaxation (T1, T1rho) experiments for all 13 genotypes, with processed Excel workbooks per experiment type. Molecular dynamics simulation code, coordinate files, and analysis scripts (NAMD/CHARMM/Python) for six atomistic cell wall models are included. Summarized ssNMR data are compiled into a single excel file and subjected to statistical analysis. Multivariate analysis code (PCA, Pearson correlation) and summary data are provided as excel worksheets and Jupyter notebooks (Python 3).

09 BIOMASS FUELS↗
DOE OSTI2022

Recovery of MOF-5 from Extreme High-Pressure Conditions Facilitated by a Modern Pressure Transmitting Medium

Here mechanisms underlying the mechanically induced amorphization of metal-organic frameworks (MOFs) are of current interest, and both high-pressure experimentation and molecular dynamics simulations have been used to reveal the fundamentals of load bearing, deformation, and pressure-induced amorphization (PIA) in these highly porous materials. Unfortunately, MOFs are typically highly susceptibility to amorphization, which limits the conditions under which they can be processed and used. However, their flexible structures can be stabilized to high pressures by incorporating guest species into the framework matrix. In this study, a large molecule pressure transmitting medium (DAPHNE 7575 ) is used as a structure-fortifying guest species to stabilize the prototypical MOF -5 to ultra-high pressures ( >9 GPa) and enable the recovery of crystalline material upon decompression. Structural changes associated with the penetration of the pressure transmitting medium on compression are examined using a combination of high pressure synchrotron powder diffraction and molecular dynamics simulations. This work enhances the understanding of PIA in MOFs while showcasing a potential route for the stabilization of MOFs to surprisingly high-pressures.

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

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dueling Backbones: Comparing Peptoid and Peptide Analogues of a Mussel Adhesive Protein

Ensembles of amino acid side chains often dominate the interfacial interactions of intrinsically disordered proteins; however, backbone contributions are far from negligible. Using a combination of nanoscale force measurements and molecular dynamics simulations, we demonstrated with analogous mussel-mimetic adhesive peptides and peptoids both 34 residues long that highly divergent adhesive/cohesive outcomes can be achieved on mica surfaces by altering backbone chemistry only. The Phe, Tyr, and Dopa containing peptoid variants used in this study deposited as dehydrated and incompressible films that facilitated analysis of peptoid side chain contributions to adhesion and cohesion. For example, whereas Phe and Dopa peptoids exhibited similar cohesion, Dopa peptoids were ~3 times more adhesive than Phe peptoids on mica. Compared with the peptides, Phe peptoid achieved only ~20% of Phe containing peptide adhesion, but the Dopa peptoids were >2-fold more adhesive than the Dopa peptides. Cation-π interactions accounted for some but not all of the cohesive interactions. Finally, our results were corroborated by molecular dynamics simulations and highlight the importance of backbone chemistry and the potential of peptoids or peptoid/peptide hybrids as wet adhesives and primers.

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

Mechanical behavior of ultralight nickel metamaterial

The mechanical response of ultralight kagomé structures consisting of hollow nickel (Ni) nanotubes and solid Ni nanorods to compression is studied using molecular dynamics simulations. In both kagomé architectures, 16[112] Shockley partial dislocations and twin formation are observed under compression. The structure made from solid nanorods shows deformation near both the nodes and beams of the kagomé lattice. The hollow kagomé architecture has a higher yield point than the solid kagomé structure. The deformation in the hollow nanotube structure is mostly localized in the nodal region for strains less than 11%. At higher strains, the deformation sets in all the struts and nodes of the hollow kagomé lattice. Owing to this two-stage deformation mechanism, the hollow Ni nanotube kagomé structure shows less bending and greater toughness than the solid Ni nanorod kagomé architecture.

Rajak, Pankaj (ORCID:0000000263446056)↗
DOE OSTI2020

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Raman studies of hydrogen trapped in As 4 O 6 ·2H 2 at high pressure and low temperature

In this work, Raman spectroscopic measurements of the arsenolite–hydrogen inclusion compound As 4 O 6 ·2H 2 were performed in diamond anvil cells at high pressure and variable temperature down to 80 K. The experimental results were complemented by ab initio molecular dynamics simulations and phonon calculations. Observation of three hydrogen vibrons in As 4 O 6 ·2H 2 is reported in the entire temperature and pressure range studied (up to 24 GPa). While the experiments performed with protium and deuterium at variable temperatures allowed for the assignment of two vibrons as Q 1 (1) and Q 1 (0) transitions of ortho and para spin isomers of hydrogen trapped in the inclusion compound, the origin of the third vibron could not be unequivocally established. Low-temperature spectra revealed that the lowest-frequency vibron is actually composed of two overlapping bands of A g and T 2g symmetries dominated by H 2 stretching modes as predicted by our previous density functional theory calculations. We observed low-frequency modes of As 4 O 6 ·2H 2 vibrations dominated by H 2 “librations,” which were missed in a previous study. A low-temperature fine structure was observed for the J = 0 → 2 and J = 1 → 3 manifolds of hydrogen trapped in As 4 O 6 ·2H 2 , indicating the lifting of degeneracy due to an anisotropic environment. A non-spherical distribution was captured by molecular dynamics simulations, which revealed that the trajectory of H 2 molecules is skewed along the crystallographic $\langle$111$\rangle$ direction. Last but not least, low-temperature synchrotron powder x-ray diffraction measurements on As 4 O 6 ·2H 2 revealed that the bulk structure of the compound is preserved down to 5 K at 1.6 GPa.

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

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ~5-50 cm -1 . This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effect of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal/mol -1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

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

Three-dimensional atomic packing in amorphous solids with liquid-like structure

Liquids and solids are two fundamental states of matter. However, our understanding of their three-dimensional atomic structure is mostly based on physical models. Here we use atomic electron tomography to experimentally determine the three-dimensional atomic positions of monatomic amorphous solids, namely a Ta thin film and two Pd nanoparticles. We observe that pentagonal bipyramids are the most abundant atomic motifs in these amorphous materials. Instead of forming icosahedra, the majority of pentagonal bipyramids arrange into pentagonal bipyramid networks with medium-range order. Molecular dynamics simulations further reveal that pentagonal bipyramid networks are prevalent in monatomic metallic liquids, which rapidly grow in size and form more icosahedra during the quench from the liquid to the glass state. Furthermore, these results expand our understanding of the atomic structures of amorphous solids and will encourage future studies on amorphous–crystalline phase and glass transitions in non-crystalline materials with three-dimensional atomic resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shape-preserving amorphous-to-crystalline transformation of CaCO3 revealed by in situ TEM

Organisms use inorganic ions and macromolecules to regulate crystallization from amorphous- precursors, endowing natural biominerals with complex morphologies and enhanced properties. The mechanisms by which modifiers enable these shape preserving transformations are poorly understood. We used in situ liquid-phase transmission electron microscopy to follow the evolution from amorphous calcium carbonate to calcite in the presence of additives. A combination of contrast analysis and infrared spectroscopy shows that Mg ions, which are widely present in seawater and biological fluids, alter the transformation pathway in a concentration-dependent manner. The ions bring excess (structural) water into the amorphous bulk so that a direct transformation is triggered by dehydration in the absence of morphological changes. Molecular dynamics simulations suggest Mg-incorporated water induces structural fluctuations, allowing transformation without the need to nucleate a separate crystal. Thus, the obtained calcite retains the original morphology of the amorphous state, biomimetically achieving the morphological control of crystals seen in biominerals

crystallization, biomineralization, calcium carbon↗
DOE OSTI2020

Facile diamond synthesis from lower diamondoids

Carbon-based nanomaterials have exceptional properties that make them attractive for a variety of technological applications. Here, we report on the use of diamondoids (diamond-like, saturated hydrocarbons) as promising precursors for laser-induced high-pressure, high-temperature diamond synthesis. The lowest pressure and temperature $(P-T)$ conditions that yielded diamond were 12 GPa (at ~2000 K) and 900 K (at ~20 GPa), respectively. This represents a substantially reduced transformation barrier compared with diamond synthesis from conventional (hydro)carbon allotropes, owing to the similarities in the structure and full sp 3 hybridization of diamondoids and bulk diamond. At 20 GPa, diamondoid-to-diamond conversion occurs rapidly within <19 μs. Molecular dynamics simulations indicate that once dehydrogenated, the remaining diamondoid carbon cages reconstruct themselves into diamond-like structures at high $P-T$. This study is the first successful mapping of the $P-T$ conditions and onset timing of the diamondoid-to-diamond conversion and elucidates the physical and chemical factors that facilitate diamond synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗
DOE OSTI2020

Nanometer-Scale Correlations in Aqueous Salt Solutions

The intermediate-range (1–4 nm) correlation of cations, an- ions, and water in aqueous alkaline earth salt solutions is measured using synchrotron X-ray diffraction. We differentiate from the entire solution structure factor, S(Q), a separate region at low Q (<1 °A-1) containing local diffraction maxima (pre-peaks) that indicate nanometer-scale oscillatory behavior. These features are quantitatively reproduced with classical molecular dynamics simulations. At high concentrations, the pre-peaks emerge from correlations arising from existence of small quasi close-packed lattice- like structures comprised of cation hydration spheres. We also analyze the concentration dependence of the pre-peak and discuss the overall results in light of the rich literature dealing with intermediate-range correlations underlying universal phenomena in concentrated electrolytes. The authors were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is operated by Battelle for the U.S. Department of Energy under Contract No. DE-AC05-76RL01830. Computing resources were generously allocated by PNNL's Institutional Computing program. DOE Contract No. DEAC02- 06CH11357 supports operation of the Advanced Photon Source at Argonne National Laboratory.

Chemical physics, nucleation, Madelung potential↗
DOE OSTI2020

Polymer Conformations and Diffusion Through a Monolayer of Confining Nanoparticles

We present coarse-grained molecular dynamics simulations to probe chain-scale polymer conformations and diffusion between confining nanoparticles (NPs). By constructing a monolayer of hexagonally-packed NPs in a polymer melt with athermal interactions, we observe the magnitude and length-scale over which homogeneously confining NPs impact the polymer behavior, which provides fundamental insights into more complex polymer nanocomposites. We show that polymer conformations are more perturbed under strong confinement (i.e. when the interparticle distance, ID, is less than twice the polymer radius of gyration, 2R g ) as compared to around an isolated NP, and the effect depends on the ratio of R NP /R g rather than either independently. In fact, these conformations can be quantitatively replicated by executing a simple random walk in a similarly confining environment. We also show that polymer diffusion is slowed by the presence of NPs and the slowing persists far beyond the length-scale over which polymer conformations are perturbed. Although the slowing is strongest ~R g from the NPs, the diffusion coefficient is slower even beyond ~5R g from the NPs. Furthermore, by analyzing the directional van Hove distributions, we show polymer diffusion away from the NP monolayer is bulk-like while diffusion through the monolayer is slower with increasing NP confinement. Furthermore, these molecular dynamics simulations provide fundamental insights into the temporal and spatial effect of confining, athermal NPs on chain-scale polymer conformations and diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Band Structure-Driven Design of a α-CsPbI 3 Ammonia Sensor for Industrial Applications

We investigate the defect-dependent electronic structure and gas-sensing potential of cubic α-CsPbI 3 using first-principles density functional theory and nonadiabatic molecular dynamics. Among the intrinsic defects, interstitials, vacancies, antisites, and switches studied, the I Pb and Pb I antisite defects exhibit transition energy levels near the middle of the band gap, thus functioning as deep traps. Short-term adsorption of ammonia selectively modifies the electronic structure, coordinating with Pb at Pb I sites and Cs at I Pb sites, significantly altering recombination pathways. Detailed analysis reveals that NH 3 reduces anharmonicity at I Pb defects, enabling enhanced recombination at elevated temperatures, while trap-assisted recombination dominates at room temperature. Other analytes, including CH 3 NH 2 and NO 2 , show negligible impact on the band gap or recombination dynamics, highlighting the potential selectivity of NH 3 interactions. Ab initio nonadiabatic molecular dynamics simulations at 300 K and 600 K further demonstrate temperature-dependent modulation of carrier lifetimes, with NH 3 accelerating recombination at ambient conditions and suppressing certain pathways at higher temperatures. These findings suggest that α-CsPbI 3 can serve as a selective and sensitive ammonia sensor over a broad temperature range and offer insights for ammonia detection under industrially relevant conditions.

ammonia sensor↗