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 1,261 records · Page 70

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗

Relationship between A-site cation and magnetic structure in 3d–5d–4f double perovskite iridates Ln 2 NiIrO 6 ( Ln = La, Pr, Nd)

We report a comprehensive investigation of Ln 2 NiIrO 6 (Ln=La, Pr, Nd) using thermodynamic and transport properties, neutron powder diffraction, resonant inelastic x-ray scattering, and density-functional theory (DFT) calculations to investigate the role of A-site cations on the magnetic interactions in this family of hybrid 3d–5d–4f compositions. Magnetic structure determination using neutron diffraction reveals antiferromagnetism for La 2 NiIrO 6 , a collinear ferrimagnetic Ni and Ir state that is driven to long-range antiferromagnetism upon the onset of Nd ordering in Nd 2 NiIrO 6 , and a noncollinear ferrimagnetic Ni and Ir sublattice interpenetrated by a ferromagnetic Pr lattice for Pr 2 NiIrO 6 . For Pr 2 NiIrO 6 , heat-capacity results reveal the presence of two independent magnetic sublattices, and transport resistivity indicates insulating behavior and a conduction pathway that is thermally mediated. Here, a first principles DFT calculation elucidates the existence of the two independent magnetic sublattices within Pr 2 NiIrO 6 and offers insight into the behavior in La 2 NiIrO 6 and Nd 2 NiIrO 6 . Resonant inelastic x-ray scattering is consistent with spin-orbit coupling splitting the t 2 g manifold of octahedral Ir 4+ into a J eff = 1/2 and J eff = 3/2 state for all members of the series considered.

36 MATERIALS SCIENCE↗

Fully self-consistent finite-temperature $GW$ in Gaussian Bloch orbitals for solids

In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc ⁢$GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc⁢ $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dense plasma opacity via the multiple-scattering method

The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. Here we calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dark matter absorption via electronic excitations

We revisit the calculation of bosonic dark matter absorption via electronic excitations. Working in an effective field theory framework and consistently taking into account in-medium effects, we clarify the relation between dark matter and photon absorption. As is well-known, for vector (dark photon) and pseudoscalar (axion-like particle) dark matter, the absorption rates can be simply related to the target material’s optical properties. However, this is not the case for scalar dark matter, where the dominant contribution comes from a different operator than the one contributing to photon absorption, which is formally next-to-leading-order and does not suffer from in-medium screening. It is therefore imperative to have reliable first-principles numerical calculations and/or semi-analytic modeling in order to predict the detection rate. We present updated sensitivity projections for semiconductor crystal and superconductor targets for ongoing and proposed direct detection experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Occupational disorder as the origin of flattening of the acoustic phonon branches in the clathrate Ba 8 ⁢Ga 16 ⁢Ge 30

In the search for high-performance thermoelectrics, materials such as clathrates have drawn attention due to having both glasslike low phonon thermal conductivity and crystal-like high electrical conductivity. Ba 8 ⁢Ga 16 ⁢Ge 30 (BGG) has a loosely bound guest Ba atom trapped inside rigid Ga-Ge cage structures. Avoided crossings between acoustic phonons and the flat guest atom branches have been proposed to be the source of the low lattice thermal conductivity of BGG. Ga-Ge site disorder with Ga and Ge exchanging places in different unit cells has also been reported. We used time-of-flight neutron scattering to measure the complete phonon spectrum in a large single crystal of BGG and compared these results with predictions of density functional theory to elucidate the effect of the disorder on heat-carrying phonons. Experimental results agreed much better with the calculation assuming the disorder than with the calculation assuming the ordered configuration. Although the atomic masses of Ga and Ge are nearly identical, we found that disorder strongly reduces phonon group velocities, which significantly reduces thermal conductivity. Finally, our work points to a path towards optimizing thermoelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental Characterization of Photoemission from Plasmonic Nanogroove Arrays

Metal photocathodes are an important source of high-brightness electron beams, ubiquitous in the operation of both large-scale accelerators and table-top microscopes. When the surface of a metal is nanoengineered with patterns on the order of the optical wavelength, it can lead to the excitation and confinement of surface-plasmon-polariton waves that drive nonlinear photoemission. In this work, we aim to evaluate gold plasmonic nanogrooves as a concept for producing bright electron beams for accelerators via nonlinear photoemission. Here, we do this by first comparing their optical properties to numerical calculations from first principles to confirm our ability to fabricate these nanoscale structures. Their nonlinear photoemission yield is found by measuring emitted photocurrent as the intensity of their driving laser is varied. Finally, the mean transverse energy of this electron source is found using the solenoid-scan technique. Our data demonstrate the ability of these cathodes to provide a tenfold enhancement in the efficiency of photoemission over flat metals driven with a linear process. We find that these cathodes are robust and capable of reaching sustained average currents over 100 nA at optical intensities larger than 2 GW cm –2 with no degradation of performance. The emittance of the generated beam is found to be highly asymmetric, a fact we can explain with calculations involving the also asymmetric roughness of the patterned surface. These results demonstrate the use of nanoengineered surfaces as enhanced photocathodes, providing a robust air-stable source of high-average-current electron beams hopefully with great potential for industrial and scientific applications.

47 OTHER INSTRUMENTATION↗

Quaternary MgSiN 2 -GaN alloy semiconductors for deep UV applications

Ultrawide direct band gap semiconductors hold great promise for deep ultraviolet optoelectronic applications. Here we evaluate the potential of MgSiN 2 -GaN alloys for this purpose. Although MgSiN 2 itself has an indirect gap ~0.4 eV below its direct gap of ~6.5 eV, its different sign lattice mismatch from GaN in two different basal plane directions could avoid the tensile strain which limits Al x Ga 1–x N on GaN for high x. Two octet-rule-preserving structures (with space groups Pmn2 1 and P1n1) of a 50% alloy of MgSiN 2 and GaN are investigated and are both found to have gaps larger than 4.75 eV using quasiparticle self-consistent GW calculations. Both are nearly direct gap in the sense that the indirect gap is less than 0.1 eV lower than the direct gap. Furthermore, their mixing energies are positive yet small, with values of 8 (31) meV/atom for Pmn2 1 (P1n1), indicating only a small driving force toward phase separation.

36 MATERIALS SCIENCE↗

Efficiency of charge transfer in changing the dissociation dynamics of OD + transients formed after the photo-fragmentation of D 2 O

Here, we present an investigation of the relaxation dynamics of deuterated water molecules after direct photo-double ionization at 61 eV. We focus on the very rare D + + O + + D reaction channel in which the sequential fragmentation mechanisms were found to dominate the dynamics. Aided by theory, the state-selective formation and breakup of the transient OD + (a 1 Δ, b 1 Σ + ) is traced, and the most likely dissociation path—OD + : a 1 Δ or b 1 Σ + → A 3 Π → X 3 Σ – → B 3 Σ – —involving a combination of spin–orbit and non-adiabatic charge transfer transitions is determined. The multi-step transition probability of this complex transition sequence in the intermediate fragment ion is directly evaluated as a function of the energy of the transient OD + above its lowest dissociation limit from the measured ratio of the D + + O + + D and competing D + + D + + O sequential fragmentation channels, which are measured simultaneously. Our coupled-channel time-dependent dynamics calculations reproduce the general trends of these multi-state relative transition rates toward the three-body fragmentation channels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Single-crystal growth and superconductivity in RbNi 2 Se 2

Herein we report the synthesis and characterization of RbNi 2 Se 2 , an analog of the iron chalcogenide superconductor Rb x Fe 2 Se 2 , via transport, angle-resolved photoemission spectroscopy, and density functional theory calculations. A superconducting transition at T c =1.20 K is identified. In the normal state, RbNi 2 Se 2 shows paramagnetic and Fermi-liquid behaviors. A large Sommerfeld coefficient yields an effective electron mass of m* ≈ 6m e . In the superconducting state, zero-field electronic specific-heat data C es can be described by a two-gap BCS model, indicating that RbNi 2 Se 2 is a possible multigap superconductor. Our density functional theory calculations and angle-resolved photoemission spectroscopy measurements demonstrate that RbNi 2 Se 2 exhibits relatively weak correlations and multiband characteristics, consistent with the multigap superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of realistic out-of-plane dopant potentials on the superfluid density of overdoped cuprates

Recent experimental papers on hole-doped overdoped cuprates have argued that a series of observations showing unexpected behavior in the superconducting state imply the breakdown of the quasiparticle-based Landau–BCS paradigm in that doping range. In contrast, some of the present authors have argued that a phenomenological “dirty d-wave” theoretical analysis explains essentially all aspects of thermodynamic and transport properties in the superconducting state, provided the unusual effects of weak, out-of-plane dopant impurities are properly accounted for. Here we attempt to place this theory on a more quantitative basis by performing ab initio calculations of dopant impurity potentials for LSCO and Tl-2201. These potentials are more complex than the pointlike impurity models considered previously, and require calculation of forward scattering corrections to transport properties. Including realistic, ARPES-derived band structures, Fermi liquid renormalizations, and vertex corrections, we show that the theory can explain semiquantitatively the unusual superfluid density measurements of the two most studied overdoped materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

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↗

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↗

$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↗

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↗