Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculation”

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

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Machine-learning-aided density functional theory calculations of stacking fault energies in steel

A combined large-scale first principles approach with machine learning and materials informatics is proposed to quickly sweep the chemistry-composition space of advanced high strength steels (AHSS). AHSS are composed of iron and key alloying elements such as aluminum and manganese. A systematic exploration of the distribution of aluminum and manganese atoms in iron is used to investigate low stacking fault energies configurations using first principles calculations. To overcome the computational cost of exploring the composition space, this process is sped up using an automated machine learning tool: DeepHyper. Here our results predict that it is energetically favorable for Al to stay away from a stacking fault, but Mn atoms do not affect the stacking fault energy and can stay in the vicinity of the fault. The distribution of Al and Mn atoms in systems containing stacking faults and the effects of their interactions on the equilibrium distribution are systematically analyzed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Symmetry progression and band vs Mott character of CdPS 3 under pressure

Complex chalcogenides are renowned for their tunable electronic, magnetic, and optical properties under external stimuli. The MPX 3 family (M = Mn, Ni, Co, V; X = S, Se) is a platform for many exciting discoveries—especially under compression—although CdPS3 is thought to be different because the Cd center possesses a filled 4d shell, which precludes Mottness. Here, we combine synchrotron-based infrared absorbance and Raman scattering spectroscopies with diamond anvil cell techniques, complementary lattice dynamics calculations, and an analysis of the energy landscape to reveal a series of structural phase transitions in CdPS 3 . We find four distinct pressure-driven transitions, with low frequency modes detectable over the full 35 GPa range of our investigation. A group–subgroup analysis along with our first-principles calculations allows us to partially unravel the space group sequence. For instance, the first critical pressure is a monoclinic C2/m to trigonal $\overline{R}$3 transition at 10 GPa. Despite the softness and overall sensitivity to pressure, we do not locate an insulator-to-metal transition in this pressure range, indicating that the energy scale for gap closure is significantly higher than expected. We discuss these findings in terms of force-induced color change and Mott vs band character in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Data-Driven Compression of Electron-Phonon Interactions

First-principles calculations of electron interactions in materials have seen rapid progress in recent years, with electron-phonon ( e − ph ) interactions being a prime example. However, these techniques use large matrices encoding the interactions on dense momentum grids, which reduces computational efficiency and obscures interpretability. For e − ph interactions, existing interpolation techniques leverage locality in real space, but the high dimensionality of the data remains a bottleneck to balance cost and accuracy. Here we show an efficient way to compress e − ph interactions based on singular value decomposition (SVD), a widely used matrix and image compression technique. Leveraging (un)constrained SVD methods, we accurately predict material properties related to e − ph interactions—including charge mobility, spin relaxation times, band renormalization, and superconducting critical temperature—while using only a small fraction (1%–2%) of the interaction data. These findings unveil the hidden low-dimensional nature of e − ph interactions. Furthermore, they accelerate state-of-the-art first-principles e − ph calculations by about 2 orders of magnitude without sacrificing accuracy. Our Pareto-optimal parametrization of e − ph interactions can be readily generalized to electron-electron and electron-defect interactions, as well as to other couplings, advancing quantitative studies of condensed matter. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

From closed shells to open shells: Coupled-cluster calculations of atomic nuclei

Coupled-cluster theory is a powerful tool for first-principles calculations of atomic nuclei, enabling accurate predictions of nuclear observables across the Segrè chart. While coupled-cluster computations are especially efficient at shell closures, extensions have been developed to tackle open-shell nuclei, by exploiting the equation-of-motion method or by expanding the coupled-cluster wave function on top of a symmetry-breaking (either deformed or superfluid) reference state. In this study, we provide a comprehensive comparison of these different formulations applied to the calcium and nickel isotopes using nuclear two-and three-body interactions from chiral effective field theory. Here, based on ground-state energies, two-neutron separation energies, and two-neutron shell gaps, different coupled-cluster computations—based on symmetry-broken reference states and equationof-motion techniques— offer consistent descriptions of bulk properties across medium-mass isotopic chains.

Marino, Francesco [Johannes Gutenberg-Universität ↗
DOE OSTI2025

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory↗
DOE OSTI2024

How spin relaxes and dephases in bulk halide perovskites

Abstract Spintronics in halide perovskites has drawn significant attention in recent years, due to their highly tunable spin-orbit fields and intriguing interplay with lattice symmetry. Here, we perform first-principles calculations to determine the spin relaxation time ( T 1 ) and ensemble spin dephasing time ( $${T}_{2}^{*}$$ T 2 * ) in a prototype halide perovskite, CsPbBr 3 . To accurately capture spin dephasing in external magnetic fields we determine the Landé g -factor from first principles and take it into account in our calculations. These allow us to predict intrinsic spin lifetimes as an upper bound for experiments, identify the dominant spin relaxation pathways, and evaluate the dependence on temperature, external fields, carrier density, and impurities. We find that the Fröhlich interaction that dominates carrier relaxation contributes negligibly to spin relaxation, consistent with the spin-conserving nature of this interaction. Our theoretical approach may lead to new strategies to optimize spin and carrier transport properties.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Exploring the interfacial structure and crystallinity for direct growth of Mn 3 Sn(0001) on sapphire (0001) by molecular beam epitaxy

The Kagome antiferromagnet Mn 3 Sn has garnered a lot of attention due to the presence of interesting properties such as anomalous Hall effects and Nernst effects. Until now, few papers have been reported to grow using molecular beam epitaxy with a buffer layer or on other substrates. In this paper, we discuss the synthesis of crystalline Mn 3 Sn layers, prepared on Al 2 O 3 (0001) without a buffer layer using molecular beam epitaxy. The growth is monitored in-situ using reflection high energy electron diffraction and measured ex-situ using X-ray diffraction, Rutherford back-scattering, and cross-sectional scanning transmission electron microscopy. The samples were deposited at 524 ± 5 °C, with an Mn: Sn atomic flux ratio of 3.2:1 for 90 minutes. Orientation relationships between the Mn 3 Sn films and the sapphire substrates are determined from in-plane and out-of-plane measurements. Our analysis indicates that the resulting film is predominantly c−plane oriented. Lastly, the samples prepared in this way were found to be discontiguous, showing a 3-dimensional morphology. According to first-principles calculations, the Mn 3 Sn exhibits a displaced Kagome structure in the very first stages of growth, for 2ML and 4ML growth on Al 2 O 3 (0001). This result is corroborated by calculating the surface formation energies and explains the observed RHEED patterns.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Plasmonically assisted channels of photoemission from metals

We analyze recently measured nonlinear photoemission spectra from Ag surfaces that reveal resonances whose energies do not scale with the applied photon energy but stay pinned to multiples of bulk plasmon energy $\hbar$ω p above the Fermi level. To elucidate these unexpected and peculiar features we investigate the spectra of plasmons generated in a solid by the optically pumped electronic polarization and their effect on photoemission. By combining quadratic response formalism for calculations of photoemission yield, a nonperturbative approach to inelastic electron scattering, and first-principles calculations for the electronic structure, we demonstrate the dependence of probability amplitude for single- and multiplasmon excitations on the basic parameters characterizing the photon pulse and the system. The resulting multiexcitation spectrum evolves towards a truncated plasmonic coherent state. Analogous concept is extrapolated to interpret plasmon generation by multiphoton excited electronic polarization. Based on this we elaborate a scenario that the thus created real plasmons act as supplementary frequency-locked pump field for non-Einsteinian plasmonically assisted channels of photoemission from metals. The established paradigm enables assignment and assessment of the observed linear $\hbar$ω p and nonlinear 2$\hbar$ω p electron yields from Ag. Such effects may be exploited for selective filtering of optical energy conversion in electronic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

pyGWBSE: a high throughput workflow package for GW-BSE calculations

Abstract We develop an open-source python workflow package, py GWBSE to perform automated first-principles calculations within the GW-BSE (Bethe-Salpeter) framework. GW-BSE is a many body perturbation theory based approach to explore the quasiparticle (QP) and excitonic properties of materials. GW approximation accurately predicts bandgaps of materials by overcoming the bandgap underestimation issue of the more widely used density functional theory (DFT). BSE formalism produces absorption spectra directly comparable with experimental observations. py GWBSE package achieves complete automation of the entire multi-step GW-BSE computation, including the convergence tests of several parameters that are crucial for the accuracy of these calculations. py GWBSE is integrated with Wannier90 , to generate QP bandstructures, interpolated using the maximally-localized wannier functions. py GWBSE also enables the automated creation of databases of metadata and data, including QP and excitonic properties, which can be extremely useful for future material discovery studies in the field of ultra-wide bandgap semiconductors, electronics, photovoltaics, and photocatalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Vibrational response and temperature of shock-compressed Pt: In situ extended x-ray absorption fine structure measurements to 325 GPa

Single pulse (~100 ps duration) extended x-ray absorption fine structure (EXAFS) measurements were obtained in laser shocked Pt to examine atomic vibrational properties and to determine temperatures for shock pressures from 72-325 GPa. Pt x-ray absorption spectra were recorded while a planar shock wave was propagating through the Pt sample. The absorption spectrum corresponding to the shocked Pt was obtained by subtracting the ambient Pt contribution from the measured absorption spectrum. Fits to the shocked state EXAFS data provided Pt lattice parameters consistent with the known Pt Hugoniot and the mean squared relative displacements (MSRDs) for nearest neighbor Pt atoms. Pt temperatures in the shocked state, estimated from the MSRDs using the correlated Debye model and correlated Einstein models, were consistent with both the Hugoniot temperatures calculated by integration along the Hugoniot and with the Hugoniot temperatures from published first-principles calculations. However, the agreement between the Hugoniot temperatures from EXAFS measurements and the calculated Hugoniot temperatures is not as good above 200 GPa, likely due to anharmonic effects. Furthermore, the present results demonstrate that single-pulse synchrotron EXAFS measurements in laser-shocked solids are useful for quantitative temperature determination and for examination of vibrational properties, including anharmonicity, at extreme pressure-temperature conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data

Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles investigation of high capacity, rechargeable CF x cathode batteries based on graphdiyne and “holey” graphene carbon allotropes

Batteries composed of CF x cathodes have high theoretical specific capacities (>860 mA h g -1 ). Attempts at realizing such batteries coupled with Li anodes have failed to deliver on this promise, however, due to a discharge voltage plateau below the theoretical maximum lowering the realized energy density and difficulties with recharging the system. Here, in this study, we use first-principles calculations to investigate novel carbon allotropes for these battery systems: graphdiyne and “holey” graphene. We first identify stable flourination structures and calculate their band gaps. We demonstrate that the holes in these carbon allotropes can induce the formation of an amorphous LiF network within the carbon and that this formation may, in fact, be kinetically favored. For structures where amorphous LiF forms within the carbon, we predict it is easier to recharge and higher discharge voltages can be achieved. If the LiF forms outside the carbon product, however, it will be crystalline in form and lead to lower discharge voltages and more difficulty in recharging the systems. Finally, we simulate XPS spectra of representative cases, demonstrating an experimental pathway for determining the reaction pathway of these systems. Our work suggests CF x allotropes with holes in them as potential targets for high capacity, rechargeable cathodes for Li batteries, provided they lead to the formation of amorphous LiF within the C structure.

25 ENERGY STORAGE↗
DOE OSTI2021

Intersystem crossing and exciton–defect coupling of spin defects in hexagonal boron nitride

Despite the recognition of two-dimensional (2D) systems as emerging and scalable host materials of single-photon emitters or spin qubits, the uncontrolled, and undetermined chemical nature of these quantum defects has been a roadblock to further development. Leveraging the design of extrinsic defects can circumvent these persistent issues and provide an ultimate solution. Here, we established a complete theoretical framework to accurately and systematically design quantum defects in wide-bandgap 2D systems. With this approach, essential static and dynamical properties are equally considered for spin qubit discovery. In particular, many-body interactions such as defect–exciton couplings are vital for describing excited state properties of defects in ultrathin 2D systems. Meanwhile, nonradiative processes such as phonon-assisted decay and intersystem crossing rates require careful evaluation, which competes together with radiative processes. From a thorough screening of defects based on first-principles calculations, we identify promising single-photon emitters such as Si VV and spin qubits such as Ti VV and Mo VV in hexagonal boron nitride. This work provided a complete first-principles theoretical framework for defect design in 2D materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Searching for ternary magnetic Zr–Fe–B compounds through deep machine learning

We use deep machine learning (ML) combined with first-principles calculations to search for energetically favorable ternary magnetic zirconium–iron borides. We show that an iterative ML approach enables efficient screening of vast structural libraries, effectively selecting promising candidates for subsequent first-principles investigations. Twenty-two new Fe-rich ternary compounds with formation energies within 60 meV atom −1 above the known ternary convex hull and with magnetic polarization larger than 0.6 T are identified, among which ten structures exhibit significant uniaxial anisotropy with magnetocrystalline anisotropy constant K 1 ⩾ 0.8 MJ m −3 , including a Zr 2 Fe 14 B phase. Such an ML-guided approach dramatically accelerates the discovery of rare-earth-free permanent magnetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗