Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Alpha-relaxation by scattering and medium-range atomic correlation in simple liquids

The relaxation dynamics of liquid and glass can be studied by inelastic x-ray or neutron scattering through the intermediate scattering function F(Q, t), where Q is the momentum transfer of scattering. Because of the time-consuming nature of these measurements, F(Q, t) is usually measured only at the first peak of the structure function S(Q), and its principal decay time is referred to as the α-relaxation time τ α . τ α is generally considered to describe the relaxation time of the bulk, which is related to viscosity and is controlled by the atomic cage around an atom. Here, through simulations on metallic liquids, we show that the α-relaxation time determined by scattering experiments does not purely reflect viscosity but is influenced by changes in spatial cooperativity. We also demonstrate that atomic caging is not exerted by the nearest neighbors but involves more cooperative atomic dynamics of the atomic medium-range order.

Glass transitions↗
DOE OSTI2020

Formula for the ionisation rate of an atom or ion in a strong electromagnetic field for numerical simulation

A formula is proposed for calculating the rate of field ionisation of an atom or ion, taking into account both tunnel ionisation and barrier suppression ionisation. Compared with the previous formula proposed in 2018, it more accurately describes the transition region between both ionisations and is determined mainly by two parameters: the ionisation potential of an atom or ion and the amplitude of an external electric field. This makes the presented formula suitable for use in numerical packages simulating the interaction of high-power laser radiation with matter by the particle-in-cell (PIC) method. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Elucidating ion capture and transport mechanisms of Preyssler anions in aqueous solutions using biased MACE-accelerated MD simulations

Equilibrium and biased multi-atomic cluster expansion (MACE) accelerated molecular dynamics (MD) simulations in aqueous solutions are performed to investigate the ion capture and transport mechanisms of the {P 5 W 30 } Preyssler anion (PA) as the smallest representative member of the extended polyoxometalate (POM) family with an internal cavity. The unique interatomic interactions present in the internal cavity vs the exterior of PA are carefully investigated using equilibrium MACE MD simulations for two representative Na(H 2 O)@PA and Na@PA complexes in aqueous solutions. Our careful analyses of radial distribution functions and coordination numbers show that the presence of confined water in Na(H 2 O)@PA has profound modulating effects on the nature of the interactions of the encapsulated ion with the oxygens of the PA cavity. Using well-converged MACE-accelerated multiple walker well-tempered metadynamics simulations with nanosecond timescales, two different associative (ion exchange) and dissociative (ion ejection) ion transport mechanisms were carefully investigated for Na + as one of the most abundant and representative ions present in seawater and saline solutions. By comparing systems with and without confined water, it was found that the presence of only one pre-encapsulated confined water in Na(H 2 O)@PA dramatically changes the free energy landscape of ion transport processes. It was also found that the contraction and dilation of the two windows present in PA directly influence the Na + and H 2 O transport. Furthermore, the results from this work are helpful, as they show a viable path toward tuning the ion exchange and transport phenomena in aqueous solutions of POM molecular clusters and frameworks.

25 ENERGY STORAGE↗
DOE OSTI2021

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA’s diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Atomic Edge-Guided Polyethylene Crystallization on Monolayer Two-Dimensional Materials

In this report we combine an advanced synthesis of two-dimensional (2D) materials (MoSe 2 ) having well-defined atomic edge configurations with ab initio and atomistic molecular dynamics (MD) simulations to study how atomic edges interact with polyethylene (HDPE) chains in a dilute solution assembly process. Our results reveal that Mo-terminated zigzag (Mo-ZZ) edges act as preferred nucleation sites and strongly interact with HDPE chains. The HDPE chains align in parallel with the Mo-ZZ edges and form arrays of lamellae that are perpendicular to the edges. Interestingly, atomic edge configurations are observed to dramatically change such interactions. The crystallization discrepancy at different edges was demonstrated on the same piece of MoSe 2 with different types of edges. The ab initio and MD simulations between n-alkane (n = 5 and 25), a segment of HDPE, and MoSe 2 suggest that the atomic structures of MoSe 2 can affect their interactions with n-alkane chains. Following the Mo-ZZ edge preferred nucleation principle, controlled long-range alignment of HDPE lamellae can be realized by creating multilayer MoSe 2 with parallel atomic steps. This research opens a pathway toward an atomic level understanding of polymer–2D nanomaterial interactions. It also bridges the gap between atomic-level and long-range mesoscopic structures and introduces a novel strategy for long-range structural control.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Towards Ab-Initio Simulations of Crystalline Defects at the Exascale Using Spectral Quadrature Density Functional Theory

Defects in crystalline solids play a crucial role in determining properties of materials at the nano, meso- and macroscales, such as the coalescence of vacancies at the nanoscale to form voids and prismatic dislocation loops or diffusion and segregation of solutes to nucleate precipitates, phase transitions in magnetic materials via disorder and doping. First principles Density Functional Theory (DFT) simulations can provide a detailed understanding of these phenomena. However, the number of atoms needed to correctly simulate these systems is often beyond the reach of many widely used DFT codes. The aim of this article is to discuss recent advances in first principles modeling of crystal defects using the spectral quadrature method. The spectral quadrature method is linear scaling with respect to the number of atoms, permits spatial coarse-graining, and is capable of simulating non-periodic systems embedded in a bulk environment, which allows the application of appropriate boundary conditions for simulations of crystalline defects. In this article, we discuss the state-of-the-art in ab-initio modeling of large metallic systems of the order of several thousand atoms that are suitable for utilizing exascale computing resourses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Local ordering in Ge/Ge–Sn semiconductor alloy core/shell nanowires revealed by extended x-ray absorption fine structure (EXAFS)

Short-range atomic order in semiconductor alloys is a relatively unexplored topic that may promote design of new materials with unexpected properties. Here, local atomic ordering is investigated in Ge–Sn alloys, a group-IV system that is attractive for its enhanced optoelectronic properties achievable via a direct gap for Sn concentrations exceeding ≈10 at. %. The substantial misfit strain imposed on Ge–Sn thin films during growth on bulk Si or Ge substrates can induce defect formation; however, misfit strain can be accommodated by growing Ge–Sn alloy films on Ge nanowires, which effectively act as elastically compliant substrates. In this work, Ge core/Ge 1–x Sn x (x ≈ 0.1) shell nanowires were characterized with extended x-ray absorption fine structure (EXAFS) to elucidate their local atomic environment. Simultaneous fitting of high-quality EXAFS data collected at both the Ge K-edge and the Sn K-edge reveals a large (≈ 40%) deficiency of Sn in the first coordination shell around a Sn atom relative to a random alloy, thereby providing the first direct experimental evidence of significant short-range order in this semiconductor alloy system. Furthermore, comparison of path length data from the EXAFS measurements with density functional theory simulations provides alloy atomic structures consistent with this conclusion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Recent advances in coupled cluster computations of open-shell atomic nuclei

In this contribution, we report on recent progress in coupled-cluster simulations of open-shell atomic nuclei using interactions consistently derived from chiral effective field theory. In particular, we compare different coupled-cluster approaches by computing binding energies and electric dipole polarizabilities in medium-mass calcium isotopes.

Marino, Francesco [Johannes Gutenberg Univ., Mainz↗
DOE OSTI2020

ReaxFF molecular dynamics simulations of electrolyte–water systems at supercritical temperature

We have performed ReaxFF molecular dynamics simulations of alkali metal–chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li+ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na+ and K+ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

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

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Ab Initio Bulk Free Energy Surface of Proper Ferroelectrics

We report a systematic and accurate approach for deriving the bulk free energy surface (FES), a function of temperature, polarization, and strain, from the first-principles density functional theory (DFT) of proper ferroelectrics. The core of our approach is the metadynamics algorithm that extracts the polarization dependence of the FES from all-atom molecular dynamics simulations without an a priori ansatz. The rest of the FES is derived from the metadynamics trajectories that span the relevant phase space. We demonstrate our approach in the case of lead titanate. The errors across the phase transition, due to DFT numerics, all-atom molecular dynamics, and free energy evaluation by enhanced sampling, can be systematically controlled and are of the order of 1 meV/atom. The accuracy of the resulting ab initio FES is only limited by the adopted functional approximation of DFT.

Xie, Pinchen [Lawrence Berkeley National Laborator↗
DOE OSTI2021

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Anomalous two-photon Compton scattering

Abstract X-ray free-electron lasers can generate radiation pulses with extreme peak intensities at short wavelengths. This enables the investigation of laser–matter interactions in a regime of high fields, yet at a non-relativistic ponderomotive potential, where ordinary rules of light–matter interaction may no longer apply and nonlinear processes are starting to become observable. Despite small cross-sections, first nonlinear effects in the hard x-ray regime have recently been observed in solid targets, including x-ray-optical sum-frequency generation (XSFG), x-ray second harmonic generation (XSHG) and two-photon Compton scattering (2PCS). Nonlinear interactions of bound electrons in the x-ray range are fundamentally different from those dominating at optical frequencies. Whereas in the optical regime nonlinearities are predominantly caused by anharmonicities of the atomic potential in the chemical bonds, x-ray nonlinearities far above atomic resonances are expected to be due to nonlinear oscillations of quasi-free electrons, including inner-shell atomic electrons. While the quasi-free-electron model agrees reasonably well with the experimental data for XSFG and XSHG, 2PCS measurements have led to unexpected results: the energy of the nonlinearly scattered photons from non-relativistic electrons shows a substantial unexpected red shift in addition to the Compton shift that is well beyond that predicted by a nonlinear quantum electrodynamics model for free electrons. A potential explanation for the spectral broadening is based on a previously unexplored scattering process that involves the whole atom rather than just quasi-free electrons. A first simulation that includes the atomic binding potential was successful in describing a broadening of the spectrum of the nonlinearly scattered photons to longer wavelengths for soft x-rays. However, the same model does not show any broadening at hard x-ray wavelengths, which is in agreement with other simulation approaches. To this point no calculation has been able to reproduce the experimentally observed broadening. Here we present further experimental data of 2PCS for an extended parameter range using additional diagnostics. In particular, we present measurements of the electron momentum distribution during the interaction that strongly suggest that the spectral broadening is not caused by an increased plasma temperature. We extend our measurement of the magnitude of the red shift in beryllium to > 1.9 k e V in addition to the Compton shift expected for free electrons and expand the measurement of the angular distribution to include forward scattering angles. We also present first measurements of 2PCS from diamond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Investigation of the Effect of Temperature on the Structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations

Statistical and epidemiological data imply temperature sensitivity of the SARS-CoV2 coronavirus. However, the molecular level understanding of the virus structure at different temperature is still not clear. Spike protein is the outermost structural protein of the SARS-CoV-2 virus which interacts with the Angiotensin Converting Enzyme 2 (ACE2), a human receptor, and enters the respiratory system. In this study, we performed an all atom molecular dynamics simulation to study the effect of temperature on the structure of the Spike protein. After 200 ns of simulation at different temperatures, we came across some interesting phenomena exhibited by the protein. We found that the solvent exposed domain of Spike protein, namely S1, is more mobile than the transmembrane domain, S2. Structural studies implied the presence of several charged residues on the surface of N-terminal Domain of S1 which are optimally oriented at 10– 30°C. Bioinformatics analyses indicated that it is capable of binding to other human receptors and should not be disregarded. Additionally, we found that receptor binding motif (RBM), present on the receptor binding domain (RBD) of S1, begins to close around temperature of 40°C and attains a completely closed conformation at 50°C. We also found that the presence of glycan moieties did not influence the observed protein dynamics. Nevertheless, the closed conformation disables its ability to bind to ACE2, due to the burying of its receptor binding residues. Our results clearly show that there are active and inactive states of the protein at different temperatures. This would not only prove beneficial for understanding the fundamental nature of the virus, but would be also useful in the development of vaccines and therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Sputter-deposited low-stress boron carbide films

Significant challenges have been faced in the manufacturing of low-stress B 4 C films. Here, we demonstrate a set of process parameters for direct-current magnetron sputter deposition of smooth, high-purity, amorphous B 4 C films with near-zero total residual stress and with thicknesses up to 10 μm. Films are characterized by a combination of high-energy ion scattering, x-ray diffraction, electron microscopy, and stress measurements based on substrate curvature. In order to facilitate the process transfer to other sputtering geometries, the favorable deposition conditions are correlated with distributions of landing energies and incident angles of depositing species estimated by Monte Carlo simulations of ballistic collisions and gas phase atomic transport. Based on such simulations, a decrease in compressive stress with increasing Ar working gas pressure is attributed to the corresponding broadening of the angular distribution of depositing species and associated atomic shadowing effects.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Reinforcement Learning-Guided Long-Timescale Simulation of Hydrogen Transport in Metals

Diffusion in alloys is an important class of atomic processes. However, atomistic simulations of diffusion in chemically complex solids are confronted with the timescale problem: the accessible simulation time is usually far shorter than that of experimental interest. In this work, long-timescale simulation methods are developed using reinforcement learning (RL) that extends simulation capability to match the duration of experimental interest. Two special limits, RL transition kinetics simulator (TKS) and RL low-energy states sampler (LSS), are implemented and explained in detail, while the meaning of general RL are also discussed. As a testbed, hydrogen diffusivity is computed using RL TKS in pure metals and a medium entropy alloy, CrCoNi, and compared with experiments. The algorithm can produce counter-intuitive hydrogen-vacancy cooperative motion. We also demonstrate that RL LSS can accelerate the sampling of low-energy configurations compared to the Metropolis–Hastings algorithm, using hydrogen migration to copper (111) surface as an example.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Insights into radiation resistance of titanium alloys from displacement cascade simulations

Radiation damage in beam window materials limits the use of high power proton beams in high energy physics research. The alloy Ti-6Al-4V is presently used as a beam window material but a prospective alternative, Ti-15V-3Cr-3Sn-3Al has been proposed. Since both these alloys contain dual phases at room temperature, we compare the radiation damage in the α and ß-phases of these two materials via primary knock-on atom (PKA) cascade simulations in the 10-40keV energy range. At PKA energies 30 and 40keV, the number of Frenkel pairs in the ballistic stage is higher in the ß-phase of Ti-6Al-4V than that in the ß-phase of Ti-15V-3Cr-3Sn-3Al almost by a factor of 2. The α-phase, of both these alloys, by far outperforms the ß-phases of the two alloys, both in terms of damage during the ballistic stage and in terms of the surviving defects. The average displacement threshold energy (E d ) in the α-phase of both alloys was found to be 66eV while that in the ß-phase of Ti-15-3 was 55 eV and in the ß-phase of Ti-6-4 was 46. Further, while the number of surviving defects is almost equal in both alloys, the vacancy and interstitial clustering mechanisms differ notably, which can impact the degree of radiation hardening and loss of ductility. Our simulations show larger vacancy and interstitial clusters form in Ti-6Al-4V as compared to that in Ti-15-3-3-3 alloy. These results indicate that Ti-15-3-3-3 alloys may be a promising candidate for next generation beam window material with a higher radiation tolerance than the existing Ti-6-4 alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗