Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initio molecular dynamics”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Resonant x-ray absorption of strong-field-ionized CF 3 Br

We report on an experimental and theoretical study of strong-field laser ionization of CF 3 Br followed by resonant x-ray absorption at the Br K-edge. Distinct 1s -> 4p, 5p Rydberg transitions of Br q+ (q = 1-4) atomic ions are observed and identified with Hartree-Fock-Slater and relativistic configuration interaction calculations. Time-dependent density functional theory and ab initio molecular dynamics calculations were performed to simulate the dissociative ionization process and the molecular orbitals for the q = 1-4 charge states. Measurements were made with both parallel and perpendicular linear polarizations of the laser and x-rays, but dichroism was not observed, indicating negligible alignment by the laser ionization process. This result is explained by calculations on atomic Br and the molecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Impacts of vacancy-induced polarization and distortion on diffusion in Solid Electrolyte Li 3 OCl

We report that Lithium-rich oxychloride antiperovskites are promising solid electrolytes for enabling next-generation batteries. Here, we report a comprehensive study varying Li + concentrations in Li 3 OCl using ab initio molecular dynamics simulations. The simulations accurately capture the complex interactions between Li + vacancies (V$^'_{Li}$), the dominant mobile species in Li 3 OCl. The V$^'_{Li}$ polarize and distort the host lattice, inducing additional non-vacancy-mediated diffusion mechanisms and correlated diffusion events that reduce the activation energy barrier at concentrations as low as 1.5% V$^'_{Li}$. Our analyses of discretized diffusion events in both space and time illustrate the critical interplay between correlated dynamics, polarization and local distortion in promoting ionic conductivity in Li 3 OCl.

25 ENERGY STORAGE↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗
DOE OSTI2020

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Direct time-domain determination of electron-phonon coupling strengths in chromium

In this work, we report the results of an ultrafast, direct structural measurement of optically pumped phonons in a Cr thin film using ultrashort x-ray pulses from a free-electron laser. In addition to measuring and confirming the known long-wavelength dispersion relation of Cr along a particular acoustic branch, we are able to determine the relative phase of the phonons as they are generated. The Cr sample exhibits two generation mechanisms for the phonons: the releasing of a preexisting charge density wave at higher frequencies, and the creation of an acoustic strain pulse via laser heating that dominates at lower frequencies. Furthermore, for the latter mechanism, we are able to measure the frequency dependence of the time required to generate the phonons. To explain the observed magnitude and slope of the delays, we perform first-principles simulations in the framework of density functional perturbation theory and ab initio molecular dynamics to fit anharmonic phonon models. In conclusion, these results show that the wave-vector dependence of the electron-phonon coupling is the driving mechanism behind the delay times: Phase-space limitation leads to higher times near the zone center. The absolute magnitudes of the delay times measured are found to be much shorter than the equilibrium electron-phonon coupling times we compute, indicating that the coupling strength is greatly enhanced when the electronic system is out of equilibrium with the lattice, as has been seen in bismuth and other systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Shock compression of vanadium at extremes: Theory and experiment

Here, the equation of state (EOS) and shock compression of bulk vanadium were investigated using canonical ab initio molecular dynamic simulations, with experimental validation to 865 GPa from shock data collected at Sandia's Z Pulsed Power Facility. In simulations the phase space was sampled along isotherms ranging from 3000 K to 50 000 K, for densities between ρ = 3 and 15 g / cm 3 , with a focus on the liquid regime and the body-centered-cubic phase in the vicinity of the melting limit. The principal Hugoniot predicted from first principles is overall consistent with shock data, while it showed that current multiphase SESAME-type EOS for vanadium needed revision in the liquid regime. A more accurate SESAME EOS was developed using constraints from experiments and simulations. This work emphasizes the need to use a combined theoretical and experimental approach to develop high-fidelity EOS models for extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Onset of anharmonicity and thermal conductivity in SnSe

The anharmonicity in SnSe is investigated through the analysis of moments of 119Sn nuclear resonant inelastic x-ray scattering and ab initio molecular dynamics calculations. Experimental evidences show that the anharmonic behavior started around 300 K, substantially lower than the usually suggested structural transition at 800 K. Both theory and experiments reveal substantial lifetime broadening and frequency renormalization of the optical phonons. Additionally, thermal conductivities calculated from the temporal energy moment using the Einstein diffusion equation are in good agreement with previous experiments. The abrupt increase of the thermal power near 800 K is driven by an electronic factor and not by the enhanced anharmonicity due to structural change.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Elasticity of two-dimensional ferroelectrics across their paraelectric phase transformation

The mechanical behavior of two-dimensional (2D) materials across 2D phase changes is unknown, and the finite temperature (T) elasticity of paradigmatic SnSe monolayers—ferroelectric 2D materials turning paraelectric as their unit cell turns from a rectangle into a square—is described here in a progressive manner. To begin with, their zero–T total energy landscape gives way to (Boltzmann-like) averages from which the elastic behavior is determined. Furthermore, these estimates are complemented with results from the strain-fluctuation method, which employs the energy landscape or ab initio molecular dynamics data. All approaches capture the coalescence of elastic moduli < C 11 (T) > = < C 22 (T) > due to the structural transformation. The broad evolution and sudden changes of elastic parameters < C 11 (T) >, < C 22 (T) >, and < C 12 (T) > of these atomically thin phase-change membranes establishes a heretofore overlooked connection among 2D materials and soft matter.

2-dimensional systems↗
DOE OSTI2023

Ti-6Al-4V to over 1.2 TPa: Shock Hugoniot experiments, ab initio calculations, and a broad-range multiphase equation of state

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia’s STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 10 5 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Shock-induced metallization of polystyrene along the principal Hugoniot investigated by advanced thermal density functionals

To date, none of the ab-initio molecular dynamics simulations of polystyrene, often used as an ablator material in inertial confinement fusion targets, with the standard ground-state exchange-correlation (XC) functional in density-functional theory can satisfactorily agree with experiments in terms of reflectivity measurements. Here, we use recently developed thermal strongly constrained and appropriately normed Laplacian dependent meta-generalized gradient approximation XC density functional (T-SCAN-L) and thermal hybrid XC density functional (KDT0) to show that the inclusion of thermal and inhomogeneity effects is crucial for accurate prediction of structural evolution and corresponding insulator-metal transition (IMT) during shock compression. Optical reflectivity calculated as an indicator of IMT is in perfect accord with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Temperature-induced gradual polymorphic structural crossover in liquid indium, tin, and antimony

Temperature-induced polymorphic structural crossover in metallic liquid indium (In), tin (Sn), and antimony (Sb) is investigated by in situ high-energy x-ray diffraction and ab initio molecular dynamics simulations. The results demonstrate the existence of a temperature-induced reversible and gradual structural crossover in these three liquids, reflecting from both “static” structures, e.g., peak positions of the structure factor and pair distribution function, bond angle distribution, number of tetrahedra, cluster connection, free volume, and dynamical behaviors, including the relaxation time, diffusion coefficient, and non-Gaussian parameter. Further, the main difference from liquid In, Sn, to Sb is that not only the large-sized free volume but also the anharmonicity and the dynamical heterogeneity increase significantly with the appearance of strong covalent bonds. These findings will deepen the understanding of liquid structure and the liquid-liquid transition in metallic liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Localized basis set for plutonium

The implementation of optimal strictly localized atomic orbitals basis for plutonium (Pu) using norm-conserving pseudopotential density-functional theory (DFT) is presented. The basis set was applied to the α, β, γ, δ, δ', and ε phases of Pu, δ-Pu surface, and δ-PuGa alloys. The computed properties of the Pu phases and δ-Pu surface were in good agreement with both available experimental data and prior DFT calculations based on plane-wave methodologies. Results for the δ-PuGa alloys were also in good agreement with experimental data. Finally, the reliability of the basis set was further demonstrated by using ab initio molecular dynamics to model the diffusion coefficient and activation barrier for atomic diffusion in a δ-PuGa alloy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

First-Principles Treatment of Vibrational Broadening in X-Ray Excited Valence Band Spectra for n-SrTiO 3 (001)

The valence band maximum (VBM) is an important quantity for semiconductors as it locates the Fermi level relative to the band edge. Accurate measurement of this quantity in near-surface regions of semiconductors by photoemission is a first step toward determining the electronic properties of heterostructures involving these materials. While extrapolating the leading edge of the valence band to the energy axis in photoemission spectra is a widely used way to find the VBM, this method can be ambiguous if the leading edge exhibits multiple slopes. Another way to determine the VBM is to fit the leading edge to an appropriately broadened, cross-section modulated theoretical density of states (DOS). Three kinds of broadening that should be included for maximum accuracy are those due to: (1) finite instrumental resolution, (2) valence hole lifetime, and (3) vibrational excitations. While steps (1) and (2) are straightforward to implement, (3) is more difficult because the appropriate amount of broadening is not known a priori. Here, we demonstrate that explicit inclusion of vibrational broadening using ab initio molecular dynamics facilitates accurate VBM determination for n-SrTiO 3 (001). The total DOS is constructed by summing time-averaged projections at elevated temperature onto s - , p - , and d orbitals for the constituent atoms and modulating with the associated photoemission cross sections. Subsequent convolutions of the total DOS, first with a Gaussian of width equal to the experimental energy resolution and second with a Lorentzian to simulate valence hole lifetime effects, yield line shapes that reproduce the experimental leading edges rather well. The VBM is then given by the energy at which the vibrationally broadened total DOS (prior to the convolutions) goes to zero. The VBMs generated by this method quantitatively agree with those resulting from extrapolating from the middle of the measured leading edge for SrTiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unraveling the intrinsic atomic physics behind x-ray absorption line shifts in warm dense silicon plasmas

Here, we present a novel free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter to cover a wide range of thermodynamic conditions and photon energies including the entire x-ray range. It uses Mermin–Kohn–Sham density functional theory with exchange-correlation (XC) thermal effects taken into account via a fully temperature dependent generalized gradient approximation XC functional. The methodology incorporates a combination of the ab initio molecular dynamics (AIMD) snapshotted Kubo–Greenwood optic data with a single-atom in simulation cell calculations to close the photon energy gap between the L and K edges and extend the K-edge tail toward many-keV photon energies. This gap arises in the standard scheme due to a prohibitively large number of bands required for the Kubo–Greenwood calculations with AIMD snapshots. Kubo–Greenwood data on snapshots provide an accurate description of optic properties at low photon frequencies slightly beyond the L-edge and x-ray absorption near edges structure (XANES) spectra, while data from periodic calculations with single atoms cover the tail regions beyond the edges. To demonstrate its applicability to mid-Z materials where the standard DFT-based approach is not computationally feasible, we have applied it to opacity calculations of warm dense silicon plasmas. These first-principles calculations revealed a very interesting phenomon of red-to-blue shift in K–L (1s → 2p) and K-edge absorptions along both isotherm and isochore, which are absent in most continuum-lowering models of traditional plasma physics. This new physics phenomenon can be attributed to the underlying competition between the screening of deeply bound core electrons and the screening of outer-shell electrons caused by warm-dense-plasma conditions. We further demonstrate that the ratio of 1s → 2p to the K-edge x-ray absorptions can be used to characterize warm-dense-plasma conditions. Eventually, based on our absorption calculations, we have established a first-principles opacity table (FPOT) for silicon in a wide range of material densities and temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Shear Softening in a Metallic Glass: First-Principles Local-Stress Analysis

Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this Letter using ab initio molecular dynamics simulations of the Cu 5 0Zr 50 metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the nonaffine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stabilization of Hydroxide Ions at the Interface of a Hydrophobic Monolayer on Water via Reduced Proton Transfer

We report a joint study using surface-specific sum-frequency vibrational spectroscopy and ab initio molecular dynamics simulations, respectively, on a pristine hydrophobic (sub)monolayer hexane-water interface, namely, the hexane/water interface with varied vapor pressures of hexane and different pHs in water. We show clear evidence that hexane on water revises the interfacial water structure in a way that stabilizes the hypercoordinated solvation structure and slows down the migration of hydroxide ion (OH – ) relative to that in bulk water. Furthermore, this mechanism effectively attracts the OH – to the water-hydrophobic interface with respect to its counterion. The result illustrates the striking difference of proton transfer of hydrated OH – at the interface and in the bulk, which is responsible for the intrinsic charging effect at the hydrophobic interface.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Understanding the lattice thermal conductivity of SrTiO 3 from an ab initio perspective

Here, we present a detailed analysis of the structure dependence of the lattice thermal conductivity of SrTiO 3 . We have used both ab initio molecular dynamic simulations and density functional theory calculations to decouple the effect of different structural distortions on the thermal conductivity. We have identified two main mechanisms for tuning the thermal conductivity when a distortion is applied. First, the modification of the acoustic-modes energy dispersion when a change in the lattice parameters is imposed and second, the low energy polar modes. In particular and counterintuitively, we have found that an increase in the angle of the oxygen octahedral rotations increases the thermal conductivity due to its coupling to these polar modes.

36 MATERIALS SCIENCE↗