Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles calculations”

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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Lattice vibrations and energy landscape of the isoelectronic semiconductor series CuBr, ZnSe, GaAs, and Ge: The special case of CuBr and its $\mathcal{d}$-level chemistry

Here we have examined the lattice vibrations and the energy landscape of the isoelectronic diamond and zincblende semiconductor series CuBr, ZnSe, GaAs, and Ge. Vibrations are found to be an increasing function of ionicity, with the cation sublattice always vibrating more strongly than the anion sublattice. These findings are consistent with density functional theory (DFT) calculations of the energy landscape and temperature-dependent molecular dynamics simulations of the atomic-position fluctuations. For CuBr, inclusion of the Cu 3d Hubbard U term is found necessary to stabilize the zincblende structure and to bring its vibrational amplitudes into agreement with experiment. In addition, vibrations are found to strongly affect the CuBr near-edge x-ray absorption fine structure that we have successfully modeled by including displacements in our theoretical Bethe-Salpeter equation calculations. Reverse Monte Carlo structural refinements using large atomic configurations to simultaneously fit x-ray absorption and x-ray total-scattering data support these conclusions, and they reveal strong Cu-Br first-neighbor correlations and asymmetric distributions of interatomic distances in the temperature ranges of both negative and positive thermal expansion. Delineation of the CuBr valence band photoelectron spectrum into its Cu 3d and Br 4p states uniquely reveals their covalent mixing and further supports the DFT results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of ionic disorder on the principal shock Hugoniot

Here, the effect of ionic disorder on the principal Hugoniot is investigated using multiple scattering theory to very high pressure (Gbar). Calculations using molecular dynamics to simulate ionic disorder are compared to those with a fixed crystal lattice, for both carbon and aluminum. For the range of conditions considered here we find that ionic disorder has a relatively minor influence. It is most important at the onset of shell ionization and we find that, at higher pressures, the subtle effect of the ionic environment is overwhelmed by the larger number of ionized electrons with higher thermal energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Two-Dimensional Moiré Polaronic Electron Crystals

Two-dimensional moiré materials have emerged as the most versatile platform for realizing quantum phases of electrons. Here, we explore the stability origins of correlated states in WSe 2 /WS 2 moiré superlattices. We find that ultrafast electronic excitation leads to partial melting of the Mott states on timescales 5 times longer than predictions from the charge hopping integrals and that the melting rates are thermally activated, with activation energies of 1⁢8±3 and 1⁢3±2 meV for the one- and two-hole Mott states, respectively, suggesting significant electron-phonon coupling. A density functional theory calculation of the one-hole Mott state confirms polaron formation and yields a hole-polaron binding energy of 16 meV. These findings reveal a close interplay of electron-electron and electron-phonon interactions in stabilizing the polaronic Mott insulators at transition metal dichalcogenide moiré interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

Fitting Vibrational Energy Parameters to Falloff Experiments: CH3 + CH3 + He (0.6-2 Torr 200-298 K)

There are many reported experiments and theoretical calculations for the combination of methyl free radicals diluted in a non-reactive 'deactivator. At high pressure the ethane (C2H6) combination product is the sole product; with decreasing pressure the chemically activated C2H6 will decompose. The falloff in the observed rate coefficient is the result of the competition between collisional stabilization of the chemically activated CA and unimolecular decomposition. The dependence of the rate coefficient on pressure, temperature and collision properties is complex and can not be calculated from first principles. The understanding of this system is not only of fundamental importance but is relevant to the recent detection of methyl free radicals in the atmospheres of Saturn - and Neptune. The temperatures of these outer planet atmospheres are in the 140-200 K region with total pressures (predominately H2 and He) less than 0.2 Torr. Experimentally determined rate coefficients have been reported for this reaction at T = 296-906 K and T = 200-408 K mostly with argon as the deactivator. At T = 200 K only the high pressure rate coefficient has been determined. Complete falloff curves over a wide temperature range (200-1600 K) with a variety of weak collider models used to simulate argon as the deactivator have also been reported by Klippenstein and Harding (KH). More recently we have reported the experimental rate coefficients in the falloff region with helium as the deactivator at 200 and 298 K. In this paper we have used the calculated falloff curves reported by KH for argon to determine the average energy transferred per collision for helium in our recently reported experiments. Collision rates were converted using Lennard Jones parameters; the temperature dependence of this conversion factor is noted. The helium experiments were consistent with a down of approximately 100 cm (exp-1); the temperature dependence was slight. The magnitude of down and its temperature dependence will be compared to other systems with similar substrate complexity and reaction energetics.

Tardy, Dwight C.↗
DOE OSTI2020

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Stability, electronic quantum states, and magnetic interactions of Er 3+ ions in Ga 2 ⁢O 3

Here, we report an ab initio study of phase stability, defect formation, electronic structure, and multiple magnetic, Dzyaloshinskii-Moriya, optical, hyperfine, and crystal field interactions in erbium (Er)-doped wide band gap 𝛼- and 𝛽-gallium oxides (Ga 2 ⁢O 3 ), critically important to make a foundation for both optoelectronic and quantum information applications. The chemical, structural, mechanical, and dynamical stabilities of the pristine phases are confirmed from respective negative formation energies, negative cohesive energies, favorable elastic constants, and positive phonon frequencies. The phonon dispersions indicate that the Ga-O bonds are uniform in the 𝛼-phase, while they vary in the 𝛽-phase due to the anisotropic polyhedral movement. The defect formation energy analysis confirms that both Er-doped 𝛼- and 𝛽−Ga 2 ⁢O 3 prefer Er 3+ (neutral) state. The underestimated band gaps of the pristine phases from standard density functional theory (DFT) calculations as compared to experimental values are corrected by employing the hybrid functional calculations, resulting in the indirect band gaps of 5.21 eV in 𝛼−Ga 2 ⁢O 3 and 4.94 eV in 𝛽−Ga 2 ⁢O 3 . The site preference energy analysis indicates partial occupation of Er in the octahedral site of Ga. The anisotropic nature of hyperfine tensor coefficients of Er is similar in both phases, which may be due to the occupation of Er in the same octahedral Ga site. On the other hand, the calculated magnetic exchange interaction between two Er dopants is negative for 𝛼 and positive for 𝛽, indicating an antiferromagnetic (AFM) ground state in the former and a ferromagnetic (FM) ground state in the latter. Large values of Dzyaloshinskii-Moriya interactions (DMIs) are obtained along the 𝑥 direction in the 𝛼 and along the 𝑦 direction in the 𝛽. The large DMI may support exotic magnetic textures, a promising direction for spintronic applications. The analysis of dielectric constants and refractive indices of both pristine and Er-doped phases shows a good agreement with available experimental values. The calculated optical anisotropy is slightly higher in 𝛽 than those in 𝛼, which is due to the involvement of lower symmetry in 𝛽. The crystal field coefficients (CFCs) calculated from DFT are used to analyze 4⁢𝑓 multiplets and 4⁢𝑓 −4⁢𝑓 transitions. Thus calculated lowest energy level of the first excited state to the lowest energy level of the ground state is about 1.53 µ⁢m, which is in a good agreement with available experiments, and it falls within the quantum telecommunication wavelength range.

3-dimensional systems↗
DOE OSTI2021

Mirror symmetry breaking and lateral stacking shifts in twisted trilayer graphene

Here, we construct a continuum model of twisted trilayer graphene using ab initio density-functional-theory calculations, and we apply it to address twisted trilayer electronic structure. Our model accounts for moiré variation in site energies, hopping between outside layers, and hopping within layers. We focus on the role of a mirror symmetry present in ABA graphene trilayers with a middle layer twist. The mirror symmetry is lost intentionally when a displacement field is applied between layers, and unintentionally when the top layer is shifted laterally relative to the bottom layer. We calculate two band-structure characteristics that are directly relevant to transport measurements, namely the Drude weight and the weak-field Hall conductivity. We relate them via the Hall density to assess the influence of the accidental lateral stacking shifts currently present in all experimental devices on electronic properties. Finally, we comment on the role of the possible importance of accidental lateral stacking shifts for superconductivity in twisted trilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Efficient time-dependent orbital-free density functional theory: Semilocal adiabatic response

Orbital-free density functional theory and its time-dependent extension are efficient ab initio methods for calculating the electronic structure and dynamics of large systems. Through the calculation of the optical spectra of selected clusters, we reach three important conclusions: (1) The quality of the spectra is strongly affected by the quality of the corresponding ground-state electron density; (2) the adiabatic part of the electronic response to external perturbations can be safely evaluated at the semilocal level; and (3) the nonadiabatic, current-dependent part of the time-dependent Pauli potential is key to recover correct spectral envelopes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shell-model study of calcium isotopes toward their drip line

Here, we report in this paper a study in terms of the nuclear shell model about the location of the calcium isotopes drip line. The starting point is considering the realistic two-body potential derived by Entem and Machleidt within chiral perturbation theory at next-to-next-to-next-to-leading order (N 3 LO), as well as a chiral three-body force at next-to-next-to-leading order (N 2 LO) whose structure and low-energy constants are consistent with the two-body potential. Then we construct the effective single-particle energies and residual interaction needed to diagonalize the shell-model Hamiltonian. The calculated two-neutron separation energies agree nicely with experiment until 56 Ca, which is the heaviest isotope whose mass has been measured, and do not show any sign of two-neutron emission until 70 Ca . We discuss the role of the choice of the model space in determining the neutron drip line, and also the dependence of the results on the parameters of the shell-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Origin of light-induced metastability in ZrTe 5

Here we study the nonequilibrium electronic structure of a model Dirac semimetal ZrTe 5 by using time-and-angle resolved photoemission spectroscopy and density functional theory–based electron and phonon calculations. By measuring the electronic dispersion near the Γ point at time delays up to 10 picoseconds, we discovered that the band spectral weight does not recover during the measured temporal window, revealing the existence of a light-induced metastable state in the electronic structure of this material. Our calculations find that the photoexcited A 1⁢g phonon mode leads to a band renormalization that both supports our experimental observations at the zone center and predicts changes to the band structure outside of our experimental window, ultimately showing the evolution from a direct to an indirect gap semimetal; such band renormalization dramatically reduces the electron-hole recombination rate giving rise to the metastability in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Elucidation of Structure and Physical Properties of Pyrolytic Sugar Oligomers Derived from Cellulose Depolymerization/Dehydration Reactions: A Density Functional Theory Study

Fast pyrolysis of lignocellulosic materials is a promising research area to produce renewable fuels and chemicals. Dehydration is known to be among the most important reaction families during cellulose pyrolysis; water is the most important product. Together with water, dehydration reactions also form a range of poorly known oligomer species of varying molecular sizes, often collected as part of bio-oil water-soluble (WS) fraction. In this work, we used electronic structure calculations to evaluate the relative thermodynamic stabilities of several oligomer species resulting from up to three consecutive dehydration events from cellulose depolymerization intermediates. A library of the thermodynamically favored candidate molecular structures was compiled. Results revealed that most of the water molecules are eliminated from the non-reducing end, forming thermodynamically more stable conjugated compounds. This is consistent with results reported by other researchers in literature where dehydration reactions occur preferably at the non-reducing ends of oligomers. The physical-chemical properties of the proposed structures were estimated using quantitative structure-property relationships (QSPRs) and quantitative property-property relationships (QPPRs). The anhydro-sugars derived from cellulose are often blamed for coke formation during bio-oil hydrotreatment. Understanding their chemical structure could help to develop rational strategies to mitigate coke formation. Furthermore, the thermo-physical properties reported (boiling point, melting point, Gibb’s free energy of formation, enthalpy of formation, and solubility parameters among others) are also fundamental to conducting first principle engineering calculations to design and analyze new pyrolysis reactors and bio-oil up-grading units.

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

Exchange-Driven Intermixing of Bulk and Topological Surface States by Chiral Excitons in Bi 2 Se 3

Topological surface states (TSS) in the prototypical topological insulator (TI) Bi 2 Se 3 are frequently characterized using optical probes, but electron-hole interactions and their effect on surface localization and optical response of the TSS remain unexplored. Here, we use ab initio calculations to understand excitonic effects in the bulk and surface of Bi 2 Se 3 . We identify multiple series of chiral excitons that exhibit both bulk and TSS character, due to exchange-driven mixing. Further, our results address fundamental questions about the degree to which electron-hole interactions can relax the topological protection of surface states and dipole selection rules for circularly polarized light in TIs by elucidating the complex intermixture of bulk and surface states excited in optical measurements and their coupling to light.

42 ENGINEERING↗
DOE OSTI2023

Band structure, superconductivity, and polytypism in AuSn 4

The orthorhombic compound AuSn 4 is compositionally similar to the Dirac node arc semimetal PtSn 4 . AuSn 4 is, contrary to PtSn 4 , superconducting with a critical temperature of T c = 2.35 K. Recent measurements present indications for quasi-two-dimensional superconducting behavior in AuSn 4 . Here we present measurements of the superconducting density of states and the band structure of AuSn 4 through scanning tunneling microscopy and angular resolved photoemission spectroscopy (ARPES). The superconducting gap values in different portions of the Fermi surface are spread around Δ 0 = 0.4 meV, which is close to but somewhat larger than Δ = 1.76k B T c expected from BCS theory. We observe superconducting features in the tunneling conductance at the surface up to temperatures about 20% larger than bulk T c . The band structure calculated with density functional theory follows well the results of ARPES. The crystal structure presents two possible stackings of Sn layers, giving two nearly degenerate polytypes. This makes AuSn 4 a rather unique case with a three-dimensional electronic band structure but properties ressembling those of low-dimensional layered compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Accuracy of DFT computed oxygen-vacancy formation energies and high-throughput search of solar thermochemical water-splitting compounds

The enthalpy change involved in metal oxide reduction is a key quantity in various processes related to energy conversion and storage, and is of particular interest for computational prediction. Often this prediction involves the simulation of a high temperature reduction process with a 0K methodology like density functional theory (DFT), and it is not infrequent for the high temperature and 0K stable crystal structures to differ. This introduces a conundrum with regards to the choice of crystal structure to utilize in the computation, with approaches in the literature varying and experimental validation remaining scarce. In this work we address both the crystal structure conundrum and the experimental validation, and then apply the insights we gain to guide a high-throughput search for new materials for solar thermochemical water-splitting applications. By computing the DFT+U oxygen vacancy formation energy (ΔE vf ) of a selection of ABO 3 compounds and comparing different crystal structures for each composition, we highlight the issues that arise when the structure utilized in the computation is dynamically unstable at 0K, namely the presence of an artificial lowering of ΔE vf , and the lack of convergence of ΔE vf with cell size. We solve these limitations by identifying and employing a suitable surrogate dynamically stable structure. We then validate the predictive power of our calculations against appositely generated experimental measurements of reduction enthalpy for a series of Hubbard U values, finding an accuracy ranging between 0.2-0.6 eV/O. In light of such conclusions, we revise and expand a previous a high-throughput DFT study on ABO 3 perovskite oxides. As a result, we provide a list of candidate STCH materials, highlight trends with redox-active cation and structural distortion, and identify Mn 4+ , Mn 3+ and Co 3+ as the most promising redox-active cations.

08 HYDROGEN↗
DOE OSTI2023

Parameter-free treatment of a layered correlated van der Waals magnet: CrPS 4

The electronic and magnetic structure of CrPS 4 , a two-dimensional (2D) magnetic semiconductor is examined by employing the SCAN meta-GGA density functional. We find the resulting magnetic moment and band gap are in excellent agreement with experiment. From the bulk magnetic configurations, we confirm the experimentally observed A-type antiferromagnetic (A-AFM) ordered ground state with a magnetic moment of 2.78 µ B per chromium atom and band gap of 1.34 eV. To gain insight into the evolution of the ground state with layers, the total energy of each magnetic configuration is calculated for a variety of thicknesses. Monolayer CrPS 4 is predicted to be a ferromagnetic insulator with a band gap of 1.37 eV, and A-AFM for bilayer and trilayer, with band gaps of 1.35 and 1.30 eV, respectively. The electronic structure is reported for the single-, two-, and three-layer and bulk CrPS 4 . Finally, we explore the optical properties of the 2D structure and report the dielectric tensor components and Kerr parameters for the monolayer.

2-dimensional systems↗
DOE OSTI2022

Convergent close-coupling calculations of electron scattering on HeH +

Here, we use the molecular convergent close-coupling (MCCC) method to perform calculations of 10–1000 eV electron scattering on the ground state of HeH + . Cross sections are presented for excitation of the n = 2 – 3 singlet and triplet states (where n is the united-atoms-limit principle quantum number), as well as ionization. We also present cross sections for He + and H + ion production following dissociative excitation and ionization. The He + production cross section is compared with the measurements of Lecointre, Jureta, Urbain, and Defrance [J. Lecointre, J. J. Jureta, X. Urbain, and P. Defrance, J. Phys. B At. Mol. Opt. Phys. 47, 015203 (2014)]. We find that the MCCC results are up to 30% higher than experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Accurate tight-binding model for twisted bilayer graphene describes topological flat bands without geometric relaxation

A major hurdle in understanding the phase diagram of twisted bilayer graphene is the roles of lattice relaxation and electronic structure on isolated band flattening near magic twist angles. Here in this work, the authors develop an accurate local environment tight-binding model fit to tight-binding parameters computed from ab initio density-functional theory calculations across many atomic configurations. With the accurate parametrization, it is found that the magic angle shifts to slightly lower angles than often quoted, from around 1.05° to around 0.99°, and that isolated flat bands appear for rigidly rotated graphene layers, with enhancement of the flat bands when the layers are allowed to distort. Study of the orbital localization supports the emergence of fragile topology in the isolated flat bands without the need for lattice relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomically Precise PdSe 2 Pentagonal Nanoribbons

We report atomically precise pentagonal PdSe 2 nanoribbons (PNRs) fabricated on a pristine PdSe 2 substrate with a hybrid method of top-down and bottom-up processes. The PNRs form a uniform array of dimer structure with a width of 2.4 nm and length of more than 200 nm. In situ four-probe scanning tunneling microscopy (STM) reveals metallic behavior of PNRs with ballistic transport for at least 20 nm in length. Density functional theory calculations produce a semiconducting density of states of isolated PNRs and find that the band gap narrows and disappears quickly once considering coupling between PNR stacking layers or interaction with the PdSe 2 substrate. The coupling of PNRs is further corroborated by Raman spectroscopy and field-effect transistor measurements. The facile method of fabricating atomically precise PNRs offers an air-stable functional material for dimensional control.

36 MATERIALS SCIENCE↗