Results for “first principles calculations” — Page 70

DOE OSTI2023

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides (Final Report)

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides. Final Report

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Advances in machine-learning-based sampling motivated by lattice quantum chromodynamics

Sampling from known probability distributions is a ubiquitous task in computational science, underlying calculations in domains from linguistics to biology and physics. Generative machine-learning (ML) models have emerged as a promising tool in this space, building on the success of this approach in applications such as image, text, and audio generation. Often, however, generative tasks in scientific domains have unique structures and features—such as complex symmetries and the requirement of exactness guarantees—that present both challenges and opportunities for ML. This Perspective outlines the advances in ML-based sampling motivated by lattice quantum field theory, in particular for the theory of quantum chromodynamics. Enabling calculations of the structure and interactions of matter from our most fundamental understanding of particle physics, lattice quantum chromodynamics is one of the main consumers of open-science supercomputing worldwide. Here, the design of ML algorithms for this application faces profound challenges, including the necessity of scaling custom ML architectures to the largest supercomputers, but also promises immense benefits, and is spurring a wave of development in ML-based sampling more broadly. In lattice field theory, if this approach can realize its early promise it will be a transformative step towards first-principles physics calculations in particle, nuclear and condensed matter physics that are intractable with traditional approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Spin-liquid behavior of the three-dimensional magnetic system Ba 3 NiIr 2 O 9 with S = 1

The quantum spin liquid (QSL) is an exotic phase of magnetic materials where the spins continue to fluctuate without any symmetry breaking down to zero temperature. Among the handful reports of QSL with spin S ≥ 1, examples with magnetic ions on a three-dimensional (3D) magnetic lattice are extremely rare since both larger spin and higher dimension tend to suppress quantum fluctuations. In this work, we offer a new strategy to achieve 3D QSL with high spin by utilizing two types of transition metal ions; both are magnetically active but located at crystallographically inequivalent positions. Furthermore, we design a 3D magnetic system Ba 3 NiIr 2 O 9 consisting of interconnected corner-shared NiO 6 octahedra and face-shared Ir 2 O 9 dimer, both having triangular arrangements in a-b plane. X-ray absorption spectroscopy measurements confirm the presence of Ni 2+ (S = 1). Furthermore, our detailed thermodynamic and magnetic measurements reveal that this compound is a realization of gapless QSL state down to at least 100 mK. Ab initio calculations find a strong magnetic exchange between Ir and Ni sublattices and in-plane antiferromagnetic coupling between the dimers, resulting in dynamically fluctuating magnetic moments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Quantum oscillations in the field-induced ferromagnetic state of MnBi 2-x Sb x Te 4

The intrinsic antiferromagnetic topological insulator MnBi2Te4 undergoes a metamagnetic transition in a c-axis magnetic field. It has been predicted that ferromagnetic MnBi2Te4 is an ideal Weyl semimetal with a single pair of Weyl nodes. Here we report measurements of quantum oscillations detected in the field-induced ferromagnetic phase of MnBi 2-x Sb x Te 4 , where Sb substitution tunes the majority carriers from electrons to holes. Single-frequency Shubnikov–de Haas oscillations were observed in a wide range of Sb concentrations (0.54≤x≤1.21). The evolution of the oscillation frequency and the effective mass shows reasonable agreement with the Weyl semimetal band structure of ferromagnetic MnBi 2 Te 4 predicted by density functional calculations. Intriguingly, the quantum oscillation frequency shows a strong temperature dependence, indicating that the electronic structure depends sensitively on magnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

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↗

DOE OSTI2021

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↗

DOE OSTI2021

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↗

DOE OSTI2022

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↗

DOE OSTI2022

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↗

DOE OSTI2025

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↗

DOE OSTI2021

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↗

DOE OSTI2023

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↗

DOE OSTI2023

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↗

DOE OSTI2022

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↗

DOE OSTI2023

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↗

DOE OSTI2022

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↗

DOE OSTI2022

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↗