Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quasiparticle approximation”

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 73 records · Page 4

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

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

36 MATERIALS SCIENCE↗

Mixed Stochastic-Deterministic Approach for Many-Body Perturbation Theory Calculations

Here, we present an approach for GW calculations of quasiparticle energies with quasi-quadratic scaling by approximating high-energy contributions to the Green's function in its Lehmann representation with effective stochastic vectors. The method is easy to implement without altering the GW code, converges rapidly with stochastic parameters, and treats systems of various dimensionality and screening response. Our calculations on a 5.75° twisted MoS 2 bilayer show how large-scale GW methods include geometry relaxations and electronic correlations on an equal basis in structurally nontrivial materials.

36 MATERIALS SCIENCE↗

Reshaped Weyl fermionic dispersions driven by Coulomb interactions in MoTe 2

In this work, we report the direct evidence of impacts of the Coulomb interaction in a prototypical Weyl semimetal, MoTe 2 , that alter its bare bands in a wide range of energy and momentum. Our quasiparticle interference patterns measured using scanning tunneling microscopy are shown to match the joint density of states of quasiparticle energy bands including momentum-dependent self-energy corrections, while electronic energy bands based on the other simpler local approximations of the Coulomb interaction fail to explain neither the correct number of quasiparticle pockets nor the shape of their dispersions observed in our spectrum. With this, we predict a transition between type-I and type-II Weyl fermions with doping and resolve its disparate quantum oscillation experiments, thus highlighting the critical roles of Coulomb interactions in layered Weyl semimetals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Kinetic Inductance Photodetectors Based on Nonequilibrium Response in Superconducting Thin-Film Structures

While experimental studies of kinetic-inductance sensors have been limited so far by the temperature range near the superconducting transition, these detectors can be very sensitivity at temperatures well below the transition, where the number of equilibrium quasiparticles is exponentially small. In this regime, a shift of the quasiparticle chemical potential under radiation results in the change of the kinetic inductance, which can be measured by a sensitive SQUID readout. We modeled the kinetic inductance response of detectors made from disordered superconducting Nb, NbC, and MoRe films. Low phonon transparency of the interface between the superconductor and the substrate causes substantial re-trapping of phonons providing high quantum efficiency and the operating time of approximately 1 ms at 1 K. Due to the small number of quasiparticles, the noise equivalent power of the detector determined by the quasiparticle generation-recombination noise can be as small as approximately 10(exp -19) W/Hz(exp 1/2) at He4 temperatures.

Sergeev, A. V.↗

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗

GW band structure of monolayer MoS 2 using the SternheimerGW method and effect of dielectric environment

Monolayers of transition-metal dichalcogenides (TMD) hold great promise as future nanoelectronic and optoelectronic devices. An essential feature for achieving high device performance is the use of suitable supporting substrates, which can affect the electronic and optical properties of these two-dimensional (2D) materials. Here, we perform many-body GW calculations using the SternheimerGW method to investigate the quasiparticle band structure of monolayer MoS 2 subject to an effective dielectric screening model, which is meant to approximately describe substrate polarization in real device applications. We show that, within this model, the dielectric screening has a sizeable effect on the quasiparticle band gap, for example the gap renormalization is as large as 250 meV for MoS 2 with model screening corresponding to SiO 2 . Within the G 0 W 0 approximation, we also find that the inclusion of the effective screening induces a direct band gap, in contrast to the unscreened monolayer. We also find that the dielectric screening induces an enhancement of the carrier effective masses by as much as 27% for holes, shifts plasmon satellites, and redistributes quasiparticle weight. Our results highlight the importance of the dielectric environment in the design of 2D TMD-based devices.

36 MATERIALS SCIENCE↗

Fine structure of the isovector giant dipole resonance in Nd 142 – 150 and Sm 152

Inelastic proton scattering at energies of a few hundred MeV and very-forward scattering angles including 0° has been established as a tool for the study of electric-dipole strength distributions in nuclei. The present work reports a systematic investigation of the chain of stable even-mass Nd isotopes representing a transition from spherical to quadrupole-deformed nuclei. Extraction of the equivalent photo-absorption cross sections and analysis of their fine structure in the energy region of the isovector giant dipole resonance (IVGDR). Method: Proton inelastic scattering reactions of 200 MeV protons were measured at the iThemba Laboratory for Accelerator Based Sciences in Cape Town, South Africa. The scattering products were momentum-analyzed by the K600 magnetic spectrometer positioned at θ Lab = 0°. Using dispersion-matching techniques, energy resolutions of ΔE ≈ 40–50 keV (full width at half maximum) were obtained. After subtraction of background and contributions from other multipoles, the spectra were converted to photoabsorption cross sections using the equivalent virtual-photon method. Wavelet-analysis techniques are used to extract characteristic energy scales of the fine structure of the IVGDR from the experimental data. Fine structure of the IVGDR is observed even for the most deformed nuclei studied. Comparisons between the extracted experimental energy scales and those energy scales obtained from the quasiparticle-phonon model (QPM) and Skyrme separable random phase approximation (SSRPA) predictions provide insight into the role of different giant-resonance damping mechanisms. It can be seen that the scales in the spherical and most likely also in the deformed nuclei mainly result from the fragmentation of the one-particle-one-hole (1p1h) strength into several dominant transitions serving as doorway states. In cases where calculations beyond the 1p1h level are available, some impact of the spreading due to coupling of the two-particle-two-hole (2p2h) states to the 1p1h doorway states is observed. New virtual-photon absorption data for the chain of stable Nd isotopes and 152 Sm are presented, with a focus on the phenomenon of nonstatistical cross-section fluctuations, referred to as fine structure, in the energy region of the IVGDR. The wavelet-analysis techniques used allowed for the features of the fine structure to be quantified in the form of characteristic scales. Finally, comparisons between experimental results and model predictions indicate that Landau damping seems to be the main source of the fine structure in both the spherical and deformed nuclei, but calculations including 2p2h degrees of freedom would be beneficial to confirm this for the deformed cases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonequilibrium transport in superconducting filaments

The step-like current-voltage characteristics of highly homogeneous single-crystalline tin and indium thin filaments has been measured. The length of the samples L approximately 1 cm was much greater than the nonequilibrium quasiparticle relaxation length Lambda. It was found that the activation of a successive i-th voltage step occurs at current significantly greater than the one derived with the assumption that the phase slip centers are weakly interacting on a scale L much greater than Lambda. The observation of 'subharmonic' fine structure on the voltage-current characteristics of tin filaments confirms the hypothesis of the long-range phase slip centers interaction.

Arutyunov, K. YU.↗

Lattice collective modes from a continuum model of magic-angle twisted bilayer graphene

Here, we show that the insulating states of magic-angle twisted bilayer graphene support a series of collective modes corresponding to local particle-hole excitations on triangular lattice sites. Our theory is based on a continuum model of the magic angle flat bands. When the system is insulating at moiré band filling ν = -3, our calculations show that the ground state supports seven low-energy modes that lie well below the charge gap throughout the moiré Brillouin zone, one of which couples strongly to THz photons. The low-energy collective modes are faithfully described by a model with a local SU(8) degree of freedom in each moiré unit cell that we identify as the direct product of spin, valley, and an orbital pseudospin. Apart from spin and valley-wave modes, the collective mode spectrum includes a low-energy intraflavor exciton mode associated with transitions between flat valence and conduction band orbitals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Phonon-induced collective modes in spin-orbit coupled polar metals

We study the interplay between collective electronic and lattice modes in polar metals in an applied magnetic field aligned with the polar axis. Static spin-orbit coupling leads to the appearance of a particle-hole spin-flip continuum that is gapped at low energies in a finite field. We find that a weak spin-orbit assisted coupling between electrons and polar phonons induces the emergence of electronic collective modes. The strength of the applied magnetic field tunes the number of modes and their energies, which can lie both above and below the particle-hole continuum. For a range of field values, we identify Fano-type interference between the electronic continuum and phonons. Here, we show that signatures of these collective modes can be observed in optical spectroscopy experiments, and we provide the corresponding theoretical predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anderson Localization in the Fractional Quantum Hall Effect

The interplay between interaction and disorder-induced localization is of fundamental interest. This article addresses localization physics in the fractional quantum Hall state, where both interaction and disorder have nonperturbative consequences. We provide compelling theoretical evidence that the localization of a single quasiparticle of the fractional quantum Hall state at filling factor ν=n/(2n+1) has a striking quantitative correspondence to the localization of a single electron in the (n+1)th Landau level. By analogy to the dramatic experimental manifestations of Anderson localization in integer quantum Hall effect, this leads to predictions in the fractional quantum Hall regime regarding the existence of extended states at a critical energy, and the nature of the divergence of the localization length as this energy is approached. Within a mean field approximation, these results can be extended to situations where a finite density of quasiparticles is present.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

RT-EOM-CCSD Calculations of Inner and Outer Valence Ionization Energies and Spectral Functions

Photoelectron spectroscopy (PES) is a standard experimental method for material characterization, but its interpretation can be hampered by its reliance on standard materials. To facilitate the study of unknown systems, theoretical methods are desirable. Here we present a real-time equation-of-motion coupled cluster (RT-EOM-CC) approach for valence PES, extending our core-level development. Here, we demonstrate that RT-EOM-CC yields ionization energies and spectral functions in good agreement with experiment and CI-based methods, even for some more correlated cases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio calculations of low-energy quasiparticle lifetimes in bilayer graphene

Motivated by recent experimental results we calculate from first-principles the lifetime of low-energy quasiparticles in bilayer graphene (BLG). Here, we take into account the scattering rate arising from electron-electron interactions within the GW approximation for the electron self-energy and consider several p-type doping levels ranging from 0 to ρ ≈ 2.4 × 10 12 holes/cm 2 . In the undoped case we find that the average inverse lifetime scales linearly with energy away from the charge neutrality point, with values in good agreement with experiments. The decay rate is approximately three times larger than in monolayer graphene, a consequence of the enhanced screening in BLG. In the doped case, the dependence of the inverse lifetime on quasiparticle energy acquires a non-linear component due to the opening of an additional decay channel mediated by acoustic plasmons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Solving the Bethe-Salpeter equation on a subspace: Approximations and consequences for low-dimensional materials

It is well known that the ambient environment can dramatically renormalize the quasiparticle gap and exciton binding energies in low-dimensional materials, but the effect of the environment on the energy splitting of the spin-singlet and spin-triplet exciton states is less understood. A prominent effect is the renormalization of the exciton binding energy and optical strength (and hence the optical spectrum) through additional screening of the direct Coulomb term describing the attractive electron-hole interaction in the kernel of the Bethe-Salpeter equation. The repulsive exchange interaction responsible for the singlet-triplet splitting, on the other hand, is unscreened within formal many-body perturbation theory. However, Benedict argued that in practical calculations restricted to a subspace of the full Hilbert space, the exchange interaction should be appropriately screened by states outside of the subspace, the so-called S approximation [L. X. Benedict, Phys. Rev. B 66, 193105 (2002)PRBMDO0163-182910.1103/PhysRevB.66.193105]. Here, we systematically explore the accuracy of the S approximation for different confined systems, including a molecule and heterostructures of semiconducting and metallic layered materials. We show that the S approximation is actually exact in the limit of small exciton binding energies (i.e., small direct term) and can be used to significantly accelerate convergence of the exciton energies with respect to the number of empty states, provided that a particular effective screening consistent with the conventional Tamm-Dancoff approximation is employed. We further find that the singlet-triplet splitting in the energy of the excitons is largely unaffected by the external dielectric environment for most quasi-two-dimensional materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High-energy magnetic excitations from heavy quasiparticles in CeCu 2 Si 2

Magnetic fluctuations is the leading candidate for pairing in cuprate, iron-based, and heavy fermion superconductors. This view is challenged by the recent discovery of nodeless superconductivity in CeCu 2 Si 2 , and calls for a detailed understanding of the corresponding magnetic fluctuations. Here, we mapped out the magnetic excitations in superconducting (S-type) CeCu 2 Si 2 using inelastic neutron scattering, finding a strongly asymmetric dispersion for E ≲ 1.5 meV, which at higher energies evolves into broad columnar magnetic excitations that extend to E ≳ 5 meV. While low-energy magnetic excitations exhibit marked three-dimensional characteristics, the high-energy magnetic excitations in CeCu 2 Si 2 are almost two-dimensional, reminiscent of paramagnons found in cuprate and iron-based superconductors. By comparing our experimental findings with calculations in the random-phase approximation, we find that the magnetic excitations in CeCu 2 Si 2 arise from quasiparticles associated with its heavy electron band, which are also responsible for superconductivity. Our results provide a basis for understanding magnetism and superconductivity in CeCu 2 Si 2 , and demonstrate the utility of neutron scattering in probing band renormalization in heavy fermion metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reproducibility in G 0 W 0 calculations for solids

Ab initio many-body perturbation theory within the GW approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn–Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange–correlation functional “starting point”, are used to construct G and W, and then perturbatively corrected by the resultant GW self-energy. In practice, there are several ways to construct the GW self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO 2 , the GW fundamental gaps reported in the literature can vary by more than 1 eV depending on the GW code used. In this work, we calculate and analyze GW quasiparticle (QP) energies of these and other systems with three different GW codes: BERKELEYGW, ABINIT and YAMBO. Through a systematic analysis of the GW implementation of these three codes, we identify the primary origin of major discrepancies between codes reported in prior literature to be the different implementations the Coulomb divergence in the Fock exchange term and the frequency integration scheme of the GW self-energy. We then eliminate these discrepancies by using common numerical methods and algorithms, demonstrating that the same quasiparticle energies for a given material can be obtained with different codes, within numerical differences ascribable to the technical details of the underling implementations. This work will be important for users and developers in assessing the precision of future GW applications and methods.

36 MATERIALS SCIENCE↗