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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamical correlations leading to site and orbital selective Mott insulator transition in hydrogen doped SmNiO 3

Electron doping induces metal -to -insulator transition (MIT) in SmNiO 3 as realized by experiments. While earlier density functional theory (DFT) studies with static correlations fell short of explaining the recent MIT observations at lower hydrogen concentrations, we present a comprehensive computational investigation employing an advanced approach. We combine DFT with dynamical mean field theory (DFT + DMFT) to efficiently analyze the insulating behavior of hydrogen -doped SmNiO 3 . In contrast to previous theoretical works, our calculations predict an insulator transition occurring at a reduced doping level of H:Ni = 0.5:1. Specifically, while the DFT + U method reveals a gap opening between p-to-d orbitals, the DMFT approach highlights a gap opening between d-to-d orbitals. Our findings uncover a selective Mott transition in site and orbital characteristics, with the Ni ions proximate to the doped hydrogen exhibiting Mott -like traits. Notably, DMFT calculations highlight a pronounced dependence on Hund's parameter J, implying the presence of Hundness in the Mott insulator. Finally, this study underscores the necessity of accounting for dynamical correlations to accurately describe the electronic structure of strongly correlated electron-doped rare-earth nickelates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Carbon complexes in highly C-doped GaN

We investigate the properties of heavily C-doped GaN grown by hydride vapor phase epitaxy using both optical experiments and hybrid density functional theory calculations. Previous work has established that carbon acceptors (C N ) give rise to a yellow luminescence band near 2.2 eV along with a blue luminescence band near 2.9 eV. Photoluminescence measurements show the yellow band shifting as a function of carbon concentration, suggesting a change in the behavior of carbon species as carbon content increases. With hybrid density functional theory we calculate the electrical and optical behavior of carbon centers containing multiple carbon impurities, which may arise in heavily doped material. We compare the behavior of these complexes to the isolated centers, and find that the dicarbon donor-acceptor (C Ga –C N ) complex is a candidate to explain the shift in the yellow luminescence peak. Tricarbon complexes have high formation energies and modest binding energies, and also give rise to optical transitions that are inconsistent with the observed spectra. We also identify the split dicarbon interstitial on the gallium site as a low-energy species with a large binding energy that may act to compensate carbon acceptors. Here, local vibrational modes are calculated for carbon impurity centers, and we compare these results to recent experiments. Dicarbon and tricarbon complexes involving C Ga and C N exhibit modes that are only slightly higher than the isolated species, while carbon interstitials and related complexes give rise to vibrational modes much higher than C Ga and C N .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Theoretical Calculation Study of Nanomaterials

In recent years, theoretical simulations using computational physics and chemistry models have been used widely to predict the physical, chemical, and structural properties of nanomaterials. To reflect such rapid progress, two Special Issues on the theoretical study of nanomaterials were published. The 36 articles collected cover broad re-search fields and applications. In this Editorial, we provide an overview of these pub-lications and summarize the main findings. Based on these research outcomes, we lay out some future research directions.

First-principles calculations↗
DOE OSTI2022

$O(N)$ ab initio calculation scheme for large-scale moiré structures

Here we present a two-step method specifically tailored for band structure calculation of the small-angle moiré-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for the ground state is performed with the O(N) Krylov subspace method implemented in openmx. Second, the crystal momentum-dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of the Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures down to the first magic angle (θ = 1.08°) with high enough accuracy at affordable costs. The band structures are in good agreement with those from tight-binding models, continuum models, plane-wave pseudopotential based ab initio calculations, and experimental observations. This method is also shown to be efficient in twisted double-bilayer graphene and bilayer WSe 2 . We think this two-step method can play a crucial role in other twisted two-dimensional materials, especially those with much more complex band structure and where the effective model is hard to construct.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Revealing Hidden Orbital Pseudospin Texture with Time-Reversal Dichroism in Photoelectron Angular Distributions

We performed angle-resolved photoemission spectroscopy (ARPES) of bulk 2H-WSe 2 for different crystal orientations linked to each other by time-reversal symmetry. We introduce a new observable called time-reversal dichroism in photoelectron angular distributions (TRDAD), which quantifies the modulation of the photoemission intensity upon effective time-reversal operation. We demonstrate that the hidden orbital pseudospin texture leaves its imprint on TRDAD, due to multiple orbital interference effects in photoemission. Our experimental results are in quantitative agreement with both the tight-binding model and state-of-the-art fully relativistic calculations performed using the one-step model of photoemission. While spin-resolved ARPES probes the spin component of entangled spin-orbital texture in multiorbital systems, we unambiguously demonstrate that TRDAD reveals its orbital pseudospin texture counterpart.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2026

Impact of absorption due to zero-field splitting on loss in dielectrics: A case study in sapphire

The coherence times of superconducting qubits are limited by loss mechanisms, whose microscopic origins have remained elusive. We propose a mechanism caused by transitions between zero-field-split states of paramagnetic impurities or defects. We derive the absorption cross section for a magnetic dipole transition and apply it to calculate the bulk loss tangent. For Cr, Fe, and V impurities in sapphire, we find loss tangents at 4.5 GHz in the range of 10 −9 –10 −8 , comparable to the loss measured in experiments. This value suggests that magnetic loss may be a limiting factor in the coherence times of superconducting qubits.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Size-dependent attraction of Cu solutes to clusters formed at Ag grain boundaries

We report a size-dependent solute clustering mechanism at grain boundaries in a sputtered ultrafine-grained Ag-Cu alloy, where large Cu clusters form despite weak individual solute-solute interactions. X-ray diffraction confirms limited Cu solubility in the Ag matrix, while scanning transmission electron microscopy reveals Cu clustering at both ordinary GBs and GB junctions. Density functional theory calculations show that 12-atom Cu clusters are energetically preferred, while smaller three-atom clusters are significantly less stable. Additional calculations demonstrate a marked increase in solute-cluster attraction energy with cluster size. As a result, these findings point to a previously unrecognized pathway for grain-boundary solute clustering in immiscible systems, driven by collective solute-cluster interactions, with implications for segregation behavior and stability in nanocrystalline and ultrafine-grained alloys.

Alloys↗
DOE OSTI2021

Role of Coulomb interaction in the phase formation of fcc Ce: Correlation matrix renormalization theory

The effect of electronic Coulomb interaction on the phase formation of fcc Ce lattice is investigated by full ab initio calculations without adjustable Coulomb U and J parameters using the Gutzwiller wavefunction-based correlation matrix renormalization theory (CMRT). Its total energy and pressure as a function of volume agree reasonably well with existing DFT+Gutzwiller calculations and experiments, indicating correct capture of electronic correlation and screening effects within the CMRT formalism. Here, a stable phase is found in line with the experimental α-Ce phase, and a lurking phase is identified supposedly linked with the experimental γ-Ce phase. A criterion based on the local 4f electron charge fluctuation is introduced to confirm the distinct electronic correlation natures of both phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Stability and distortion of fcc LaH 10 with path-integral molecular dynamics

The synthesis of the high-temperature superconductor LaH 10 requires pressures in excess of 100 GPa, wherein it adopts a face-centered cubic structure. Upon decompression, this structure undergoes a distortion, which still supports superconductivity, but with a lower critical temperature. Previous calculations have shown that quantum and anharmonic effects are necessary to stabilize the cubic structure, but have not resolved the low pressure distortion. Using large scale path-integral molecular dynamics enabled by a machine learned potential, we show that a rhombohedral distortion appears at sufficiently low pressures. Here, we also highlight the importance of quantum zero-point motion in stabilizing the cubic structure.

08 HYDROGEN↗
DOE OSTI2021

Nonlocal and nonadiabatic Pauli potential for time-dependent orbital-free density functional theory

Time-dependent orbital-free density functional theory is an efficient ab initio method for calculating the electronic dynamics of large systems. In comparison to standard time-dependent density functional theory, it computes only a single electronic state regardless of system size, but it requires an additional time-dependent Pauli potential term. Herein we propose a nonadiabatic and nonlocal Pauli potential whose main ingredients are the time-dependent particle and current densities. Our calculations of the optical spectra of metallic and semiconductor clusters indicate that nonlocal and nonadiabatic time-dependent orbital-free density functional theory performs accurately for metallic systems and semiquantitatively for semiconductors. This paper opens the door to wide applicability of time-dependent orbital-free density functional theory for nonequilibrium electron and electron-nuclear dynamics of complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗