Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic 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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Resolving Lonsdaleite's decade-long controversy: Atomistic insights into a metastable diamond polymorph

Lonsdaleite, a theoretically proposed hexagonal diamond polymorph, has remained at the center of a five-decade scientific controversy since its 1967 identification. While some studies claim it exhibits superior hardness through compression-induced structural changes, others contend it is merely a stacking-faulted cubic diamond. Meteoritic samples and synthetic preparations have yielded conflicting evidence, with even advanced characterisation techniques like XRD and TEM failing to provide definitive proof. In this work, we employ first-principles density functional theory (DFT) and molecular dynamics (MD) simulations to generate unambiguous theoretical fingerprints through XRD, Raman, and SAED patterns that distinguish true Lonsdaleite from cubic diamond and its defective variants. Our atomistic approach quantifies the thermodynamic metastability of Lonsdaleite under realistic pressure-temperature conditions, reveals distinct spectral signatures through simulated Raman and resolves the structural ambiguity through generalised stacking fault energy analysis. By establishing clear criteria for definitive identification, this study provides long-awaited clarity to the Lonsdaleite debate while offering a robust computational framework for characterising metastable carbon phases in meteoritic, synthetic and industrial materials.

DFT↗
DOE OSTI2024

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Neural Thermal Scattering (NeTS) Modules for Graphite & Beryllium [Slides]

This presentation discusses the motivations behind the project, which include facilitating compact formulation for TSL data, extending AI knowledge, and providing advanced reactor simulation framework. It also discusses how the project accounts for beryllium and graphite complexity, including atomistic, dynamical, and neural. Additionally, new material-informed neural thermal scattering (NeTS) Modules are examined. The implications, conclusions, and plan future work for the project are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Crystallization of the P 3 Sn 4 Phase upon Cooling P 2 Sn 5 Liquid by Molecular Dynamics Simulation Using a Machine Learning Interatomic Potential

We performed molecular dynamics simulations to study the crystallization of the P3Sn4 phase from P 2 Sn 5 liquid using a machine learning (ML) interatomic potential with desirable efficiency and accuracy. Our results capture the liquid properties of P 2 Sn 5 at 1300 K, which is well above the melting temperature. The phase separation and crystallization are observed when P 2 Sn 5 liquid is cooled down below 832 and 505 K, respectively. The simulation results are in good agreement with the experimentally observed phase transformation behaviors and provide useful insights into the complex nucleation and crystallization process at the details of atomistic scale. Our work also demonstrated that ML interatomic potentials based on neural network deep learning are robust and capable of accurately describing the energetics and kinetics of complex materials through molecular dynamics simulations.

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

Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel

Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hole formation effect on surface morphological response of plasma-facing tungsten

We report simulation results on the effect of helium (He) bubble bursting-mediated surface hole formation on the surface morphological response of tungsten plasma-facing components (PFCs) in nuclear fusion devices. Our analysis is based on an atomistically informed, continuum-scale model, which is capable of accessing the spatiotemporal scales relevant to the fuzz nanostructure formation process on the surface of PFC tungsten. Our simulations account, in an empirical fashion, for two types of subsurface bubble dynamical phenomena in the nanobubble region of PFC tungsten during He plasma irradiation, involving bubble bursting and surface crater formation. We demonstrate that the hole formation effect on the PFC tungsten surface accelerates the growth rate of nanotendrils and the onset of fuzz formation. As a result, the predicted incubation time for surface nanotendril growth is shortened, in agreement with experimental data of incubation fluence at comparable plasma exposure conditions. We also explore systematically the dependence of the PFC surface morphological response on the areal density of holes introduced at regular time intervals onto the He-implanted tungsten surface, a parameter in our analysis that serves as a proxy for the rate of He bubble bursting. More importantly, our simulations capture fine surface features in the PFC tungsten surface morphology and predict that the average spacing between nanotendrils is on the order of 100 nm, consistent with the experimental findings.

Chen, Chao-Shou (ORCID:0000000307910042)↗
DOE OSTI2026

Tracking the protein conformational motions driving HIV-1 membrane fusion

HIV-1 Env (trimeric gp120/gp41) is the surface protein responsible for membrane fusion. The Env binds to the receptor proteins, which induces gp120 shedding leading to conformational changes of gp41 from the pre-fusion to post-fusion state, allowing its fusion peptide to embed in the host cell membrane and bringing the viral and host cell membranes together. The gp41 refolding is a target of several peptide inhibitors. Yet, the molecular mechanism of this dynamic process is still not well understood. In this study, we successfully simulate the conformational change of gp41 from pre-fusion to post-fusion state in atomistic resolution using all-atom structure-based models. We reveal that maintaining the directionality of protomer interactions in both pre-fusion and post-fusion states is crucial for gp41 refolding. Additionally, we find that HR1 inherently extends as a three-helical bundle toward the host-cell membrane without any bias. Importantly, we identify native contacts in the pre-fusion state that are critical for the proper refolding of gp41 towards the post-fusion state. Lastly, by incorporating the membrane-fusion inhibitors, T20 and SFT, we identify the most vulnerable stage in the fusion pathway that exhibits the greatest sensitivity to these drugs, which could aid in a better understanding of drug resistance mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Localization of Zn 2+ ions affects the structural folding and mechanics of Nereis virens Nvjp-1

Several biological organisms utilize metal-coordination bonds to produce remarkable materials, such as the jaw of the marine worm Nereis virens, where metal-coordination bonds yield remarkable hardness without mineralization. Though the structure of a major component of the jaw, the Nvjp-1 protein, has recently been resolved, a detailed nanostructural understanding of the role of metal ions on the structural and mechanical properties of the protein is missing, especially with respect to the localization of metal ions. In this work, atomistic replica exchange molecular dynamics with explicit water and Zn 2+ ions and steered molecular dynamics simulations were used to explore how the initial localization of the Zn 2+ ions impacts the structural folding and mechanical properties of Nvjp-1. We found that the initial distribution of metal ions for Nvjp-1, and likely for other proteins with high amounts of metal-coordination, has important effects on the resulting structure, with larger metal ion quantity resulting in a more compact structure. These structural compactness trends, however, are independent from the mechanical tensile strength of the protein, which increases with greater hydrogen bond content and uniform distribution of metal ions. Our results indicate that different physical principles underlie the structure or mechanics of Nvjp-1, with broader implications in the development optimized hardened bioinspired materials and the modeling of proteins with significant metal ion content.

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

Large-scale atomistic studies of sliding friction in polycrystalline aluminum interfaces

In this work, we discuss large-scale non-equilibrium molecular dynamics (NEMD) simulations of ductile metal sliding comprising up to 1.8 × 10 9 atoms over time scales of 100 ns. The results of these simulations have identified a variety of physical mechanisms that are important in determining the steady-state frictional force for a wide range of velocities at compressed metal–metal interfaces. These include grain growth and refinement, the evolution of large plastic strains and strain rates, material mixing, and melting. These phenomena can be included in a strain, strain rate, and grain size model that gives good agreement with the NEMD simulations and can be applied to macroscopic continua.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Defect engineering of porous aromatic frameworks via end capping improves dioxane removal from water

Amorphous porous organic polymers show promise for energy-efficient adsorptive separations, but it is difficult to understand or improve their performance through intentional structural modification. Herein, we report the synthesis of porous aromatic frameworks (PAFs) with pore structures tailored by the incorporation of a monofunctionalized end-capping monomer that disrupts framework topology. Combining experimental characterization with molecular simulations, we show that this defect engineering strategy yields less densely crosslinked networks, which leads to pore collapse and the presence of unique adsorption sites. These defect-engineered PAFs exhibit enhanced removal of 1,4-dioxane, an important environmental pollutant, from water. Increasing the concentration of end-capping monomers produces PAFs with narrower pore size distributions and improved 1,4-dioxane uptake. These results illustrate that defect engineering can effectively modulate polymer connectivity and porosity for applications in selective adsorptive separations. In conclusion, this technique avoids post-synthetic treatments and presents another approach to tailor amorphous polymeric adsorbents.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nanotwin assisted reversible formation of low angle grain boundary upon reciprocating shear load

Severe plastic deformation of metals is known to lead to superior properties that cannot be achieved by any traditional metallurgic process. Origin of the superior properties is perceived to be closely associated with grain refinement, a fundamental process during the severe plastic deformation, which is essentially the formation of new grain boundaries. However, the atomistic mechanism of grain boundary formation remains largely obscure. Here, by using in-situ transmission electron microscopy and molecular dynamic simulation, we reveal, for the first time at atomic level, the shear-induced low-angle grain boundary (LAGB) formation processes in Au nanocrystal. In this work, we discover the LAGB formation is accomplished through inward propagation of nanotwins accompanied by dislocations gliding on the twin boundaries, a nanotwin-mediated dislocation slip mechanism, which shows reversible characteristic under reciprocating shear load and is affected by the nanocrystal microstructure and orientation. Our result unveils unprecedented atomistic insights on shear driven grain refinement towards nanostructure of superior properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Chemistry↗
DOE OSTI2020

Atomistic insights into the synergistic effects of tensile strain and hydroxyl group on increasing the thermal conductivity of monohydric alcohols as latent heat storage materials

The low thermal conductivity of linear monohydric alcohols significantly limits their performance when used as latent heat storage materials. It was found that imposing certain mechanical strains on nanomaterials with highly ordered structures can lead to appreciable thermal conductivity enhancement. To explore the potential of enhancing the intrinsic thermal conductivity of monohydric alcohols, a simplified ideal crystal model was employed for molecular arrangement in this work. The variation of atomistic heat transfer as a function of tensile strains was exploited by non-equilibrium molecular dynamics simulations. Furthermore, the influence of stretching on the vibration density of states and the corresponding molecular morphology were revealed. The results suggested that the thermal conductivity of monohydric alcohols with such simplified ideal crystalline structures can increase by about 170 % at the strain of 0.1. The distance between neighboring interfaces was found to determine the energy of inter-molecular interactions at the crystal grain boundaries. The hydrogen bonds energy has not reached the critical value when increasing tensile strains from 0 to 0.1. In contrast, the strongest van der Walls interactions for alkanes were obtained at the strain of 0.06 where the interfacial distances of alkanes have achieved their equilibrium distances. In addition, the interfacial heat transfer coefficient of ideal crystal monohydric alcohols at hydroxyl interfaces was found to be threefold than that between methyl interfaces, in consistent with the strength of inter-molecular interactions between functional groups at the interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Inter-Subunit Dynamics Controls Tunnel Formation During the Oxygenation Process in Hemocyanin Hexamers

Hemocyanin from horseshoe crab in its active form is a homo-hexameric protein. It exists in open and closed conformations when transitioning between deoxygenated and oxygenated states. Here, we present a detailed dynamic atomistic investigation of the oxygenated and deoxygenated states of the hexameric hemocyanin using explicit solvent molecular dynamics simulations. We focus on the variation in solvent cavities and the formation of tunnels in the two conformational states. By employing principal component analysis and CVAE-based deep learning, we are able to differentiate between the dynamics of the deoxy- and oxygenated states of hemocyanin. Finally, our results identify the deoxygenated open conformation, which adopts a stable, closed conformation after the oxygenation process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Atomistic Characterization of Gramicidin Channel Formation

We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA–gA stacking orientations that were stabilized by four or two hydrogen bonds (4HB or 2HB). When examining the temporal evolution of the dimerization, we found that two bilayer-inserted gA subunits can form the 6HB dimer directly, with no discernible intermediate states, as well as through paths that involve the 2HB and 4HB dimers.

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

Multiscale thermal properties prediction in the Multiphysics Object Oriented Simulation Environment (MOOSE) via a general Boltzmann solver [Poster]

The project objectives were: (1) Engineering-scale fuel performance modeling relies on accurate thermal properties; (2) Thermal properties are inherently multi-scale, arising from atomistic processes and interactions with a material’s microstructure; (3) Heat transport in solids via conduction occurs through transport and scattering of electrons and phonons; (4) Use the Boltzmann transport equation (BTE) to predict the macroscopic behavior of a materials system in terms of the microscopic dynamics of its heat carriers; and (5) This project establishes a new MOOSE (Multiphysics Object Oriented Simulation Environment) module, Boltzmann, dedicated to phonon and thermal electron transport.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Revealing Au13 as Elementary Clusters During the Early Formation of Au Nanocrystals

Understanding the formation mechanism of nanocrystals in solution is very fundamental to the development of materials science. For a metal nanocrystal, the cluster mediated formation mechanism is still poorly understood. Particularly, identifying what types of clusters are dominant and how they evolve into a nanocrystal in the early nucleation stage remains a great challenge. Here, using liquid-cell transmission electron microscopy, we directly observe the formation of ultra-small Au clusters (~0.84 nm) in the presence of PAA-Na. These clusters, which correspond to the size of Au13 cluster, coalesce to form nanocrystals. Our molecular dynamics simulations suggest Au13 in an aqueous environment have higher stability when compared to other cluster sizes and provide atomistic details of growth by cluster coalescence. Collectively, our demonstration of Au13 as the dominant species with elaboration of their coalescence kinetics sheds light on nonclassical nanocrystal formation mechanisms and offers useful guidelines for designing innovative pathways for synthesis of nanomaterials.

Au13 clusters, cluster coalescence, liquid-cell TE↗