Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bethe-Salpeter equation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗
NASA NTRS1992

Solution of two-body relativistic bound state equations with confining plus Coulomb interactions

Studies of meson spectroscopy have often employed a nonrelativistic Coulomb plus Linear Confining potential in position space. However, because the quarks in mesons move at an appreciable fraction of the speed of light, it is necessary to use a relativistic treatment of the bound state problem. Such a treatment is most easily carried out in momentum space. However, the position space Linear and Coulomb potentials lead to singular kernels in momentum space. Using a subtraction procedure we show how to remove these singularities exactly and thereby solve the Schroedinger equation in momentum space for all partial waves. Furthermore, we generalize the Linear and Coulomb potentials to relativistic kernels in four dimensional momentum space. Again we use a subtraction procedure to remove the relativistic singularities exactly for all partial waves. This enables us to solve three dimensional reductions of the Bethe-Salpeter equation. We solve six such equations for Coulomb plus Confining interactions for all partial waves.

Maung, Khin Maung↗
DOE OSTI2025

Renormalization-group approach to Kohn-Luttinger superconductivity: Amplification of the pairing gap from ℓ 4 to ℓ

Here, we revisit the renormalization group (RG) analysis of the Kohn-Luttinger (KL) mechanism for superconductivity. The KL mechanism leads to superconductivity in a system with a repulsive bare interaction. The key ingredient is the screening effect that renders the induced interaction attractive in channels with nonzero angular momentum ℓ ≠ 0, thereby triggering the Bardeen-Cooper-Schrieffer (BCS) instability. According to the original argument, the resulting gap is exponentially small, with its exponent scaling as −ℓ 4 . However, the KL mechanism was originally formulated within perturbation theory, where the series is known to converge poorly in certain cases—most notably, for the 𝑝-wave paring gap induced by a repulsive 𝑠-wave contact interaction. This poor convergence may be attributed to a divergent integrand in a specific class of diagrams containing both the BCS logarithm and the Kohn anomaly, suggesting that one must resum the Kohn anomaly contributions separately from the BCS logarithm. In this work, we incorporate the Kohn anomaly contribution into the 𝛽 function of the RG equation governing the BCS instability near the Fermi surface. Our solution shows that the KL gap exponent is then proportional to −ℓ, indicating a significant enhancement of the KL mechanism beyond the previously known result. To illustrate this, we study the spin-triplet 𝑝-wave pairing gap arising from a repulsive 𝑠-wave contact interaction and compare our RG-based results with those obtained from the Bethe-Salpeter equation in perturbation theory.

nuclear matter in neutron stars↗
DOE OSTI2025

Real-Space Quantification of Exciton Localization in Acene Crystals Using Wannier Function Decomposition

Here, we introduce the Wannier function decomposition of excitons (WFDX) method to quantify exciton localization in solids within the ab initio Bethe-Salpeter equation framework. By decomposing each Bloch exciton wave function into products of single-particle electron and hole maximally localized Wannier functions, this real-space approach provides well-defined orbital- and spatial-resolved measures of both Frenkel and charge-transfer excitons at low computational cost. We apply WFDX to excitons in acene crystals, quantifying how the number of rings, the exciton spin state, and the center-of-mass momentum affect spatial localization. Additionally, we show how this real-space representation reflects structural nonsymmorphic symmetries that are hidden in standard reciprocal-space descriptions. We demonstrate how the WFDX framework can be used to efficiently interpolate exciton expansion coefficients in reciprocal-space and outline how it may facilitate evaluation of observables involving position operators, highlighting its potential as a general tool for both analyzing and computing excitonic properties in solids.

Tao, Zui [University of California, Berkeley, CA (↗
DOE OSTI2026

Hybridization of Excited Interlayer Excitons with Intralayer Excitons in Transition-Metal Dichalcogenide Heterostructures: Interplay of Orbital and Structural Symmetry

Hybridization between interlayer exciton (ILX) and intralayer excitons offers a powerful route to engineer light–matter interactions in transition-metal dichalcogenide (TMDC) heterostructures, yet the underlying mechanism and selection rules remain elusive. Here we probe exciton hybridization in MoTe2/MoSe2 heterobilayers under a tunable out-of-plane electric field and observe distinct avoided crossings in the electric-field dependent reflection spectra. We associate them with coupling between excited ILX states, including the 2p and 2s Rydberg states, and MoTe2 intralayer exciton. Supported by ab initio GW plus Bethe-Salpeter equation (GW-BSE) calculations, we identify the hybridizing states and establish optical selection rules based on the effective angular momentum of the exciton states. These findings highlight the role of excitonic binding in the observed hybridization, rather than that of simple band-mixing models, and provide a new pathway for brightening high-energy dark states.

Yao, Helen [Department of Materials Science and En↗
DOE OSTI2025

Breakdown of the Static Dielectric Screening Approximation of Coulomb Interactions in Atomically Thin Semiconductors

Coulomb interactions in atomically thin materials are remarkably sensitive to variations in the dielectric screening of the environment, which can be used to control exotic quantum many-body phases and engineer exciton potential landscapes. For decades, static or frequency-independent approximations of the dielectric response, where increased dielectric screening is predicted to cause an energy redshift of the exciton resonance, have been sufficient. These approximations were first applied to quantum wells and were more recently extended with initial success to layered transition metal dichalcogenides (TMDs). Here, we use charge-tunable exciton resonances to investigate screening effects in TMD monolayers embedded in materials with low-frequency dielectric constants ranging from 4 to more than 1000, a range of 2 orders of magnitude larger than in previous studies. In contrast to the redshift predicted by static models, we observe a blueshift of the exciton resonance exceeding 30 meV in higher dielectric constant environments. We explain our observations by introducing a dynamical screening model based on a solution to the Bethe-Salpeter equation (BSE). When dynamical effects are strong, we find that the exciton binding energy remains mostly controlled by the low-frequency dielectric response, while the exciton self-energy is dominated by the high-frequency one. Our results supplant the understanding of screening in layered materials and their heterostructures, introduce a knob to tune selected many-body effects, and reshape the framework for detecting and controlling correlated quantum many-body states and designing optoelectronic and quantum devices.

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

Tunable Electronic Energy Level Alignment and Exciton Diversity in Organic–Inorganic van der Waals Heterostructures

van der Waals stacking of two-dimensional (2D) materials offers a powerful platform for engineering material interfaces with tailored electronic and optical properties. While most van der Waals multilayers have featured inorganic monolayers, incorporating molecular monolayers introduces additional degrees of tunability and functionality. Here, in this study, we investigate hybrid bilayers composed of atomically thin perylene-based molecular crystals interfaced with monolayer transition metal dichalcogenides (TMDs), specifically MoS 2 and WS 2 . Using the ab initio many-body perturbation theory within the GW approximation and the Bethe-Salpeter equation approach, we predict emergent properties beyond those of the isolated constituent systems. Notably, we find substantial renormalization of monolayer molecular crystal band gap due to TMD-induced polarization. Furthermore, by varying the TMD monolayer, we demonstrate tuning of the energy level alignment of the bilayer and subsequent control over a diversity of lowest-energy excitons, which include strongly bound hybrid excitons and long-lived charge-transfer excitons. These findings establish organic-inorganic van der Waals heterostructures as a promising class of materials for tunable optoelectronic devices and quantum excitonic phenomena, expanding the design space for low-dimensional systems.

GW-BSE calculations↗
DOE OSTI2025

PAH101: A GW+BSE Dataset of 101 Polycyclic Aromatic Hydrocarbon (PAH) Molecular Crystals

Abstract The excited-state properties of molecular crystals are important for applications in organic electronic devices. TheGWapproximation and Bethe-Salpeter equation (GW+BSE) is the state-of-the-art method for calculating the excited-state properties of crystalline solids with periodic boundary conditions. We present the PAH101 dataset ofGW+BSE calculations for 101 molecular crystals of polycyclic aromatic hydrocarbons (PAHs) with up to ~500 atoms in the unit cell. To the best of our knowledge, this is the firstGW+BSE dataset for molecular crystals. The data records include theGWquasiparticle band structure, the fundamental band gap, the static dielectric constant, the first singlet exciton energy (optical gap), the first triplet exciton energy, the dielectric function, and optical absorption spectra for light polarized along the three lattice vectors. The dataset can be used to (i) discover materials with desired electronic/optical properties, (ii) identify correlations between DFT andGW+BSE quantities, and (iii) train machine learned models to help in materials discovery efforts.

Science & Technology - Other Topics↗
DOE OSTI2024

High-Temperature Observation of Intralayer, Interlayer, and Rydberg Excitons in Bulk Van Der Waals Alloy Single Crystals

Transition metal dichalcogenides exhibit remarkable optical properties due to the diverse number of strongly bound excitons, which can be fine-tuned by alloying. Despite a flurry of research activity in characterizing these excitons, a comprehensive and profound understanding of their behavior with temperature is lacking. Here, we report the rich spectrum of excitonic features within bulk van der Waals alloy Mo0.5 W0.5 S2 and Mo0.5 W0.5 Se2 single crystals through temperature-dependent reflectance spectroscopy and first-principles calculations. We observed Rydberg excitons and interlayer excitons in both the single crystals. Notably, we provide the first experimental evidence of highly energetic A' and B' excitons in Mo0.5 W0.5 S2 at room temperature. The strong carrier-phonon scattering significantly broadens the A', B', and interlayer excitons at room temperature in bulk Mo0.5 W0.5 S2 single crystal compared to its selenide. Our findings, supported by density functional theory and Bethe-Salpeter equation calculations, signify the crucial role of carrier-phonon interactions. These results open pathways for next-generation optoelectronic devices and quantum technologies operating at high temperature.

excitons↗
DOE OSTI2024

Electronic topological transitions in cadmium under pressure studied via theoretical and experimental x-ray absorption spectroscopy

Here, an electronic topological transition (ETT) in cadmium below 1 GPa is investigated in situ with experimental x-ray absorption spectroscopy and projecting calculated core-valence excitons onto the band structure. These projections are a useful application of the Bethe-Salpeter equation approach that considers many-body effects. The method described herein can be used for systems that are otherwise difficult to probe in situ; therefore, it provides a generalizable approach to identifying and understanding ETTs under high pressure. Although pressure-induced ETTs are often probed using indirect structural responses, our own x-ray diffraction and Raman studies suggest a second-order structural transition around 3 GPa but are largely insensitive to or inconclusive for the previously studied ETT in this region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Lyapunov exponent as a signature of dissipative many-body quantum chaos

A distinct feature of Hermitian quantum chaotic dynamics is the exponential increase of certain out-of-time-order correlation (OTOC) functions around the Ehrenfest time with a rate given by a Lyapunov exponent. Physically, the OTOCs describe the growth of quantum uncertainty that crucially depends on the nature of the quantum motion. Here, we employ the OTOC in order to provide a precise definition of dissipative quantum chaos. For this purpose, we compute analytically the Lyapunov exponent for the vectorized formulation of the large- q limit of a q -body Sachdev-Ye-Kitaev model coupled to a Markovian bath. These analytic results are confirmed by an explicit numerical calculation of the Lyapunov exponent for several values of q ≥ 4 based on the solutions of the Schwinger-Dyson and Bethe-Salpeter equations. We show that the Lyapunov exponent decreases monotonically as the coupling to the bath increases and eventually becomes negative at a critical value of the coupling signaling a transition to a dynamics which is no longer quantum chaotic. Therefore, a positive Lyapunov exponent is a defining feature of dissipative many-body quantum chaos. The observation of the breaking of the exponential growth for sufficiently strong coupling suggests that dissipative quantum chaos may require in certain cases a sufficiently weak coupling to the environment. Published by the American Physical Society 2024

García-García, Antonio M.↗
DOE OSTI2025

Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange

We report on the use of the adaptively compressed exchange (ACE) operator to accelerate many-body perturbation theory (MBPT) calculations, including G 0 W 0 and the Bethe–Salpeter equation (BSE), for hybrid density functional theory starting points. We show that by approximating the exact exchange operator with the low-rank ACE operator, substantial computational savings can be achieved with systematically controllable errors in the quasiparticle energies computed with full-frequency G 0 W 0 and the optical absorption spectra and vertical excitation energies computed by solving the BSE within density matrix perturbation theory. Our implementation makes use of the ACE-accelerated electronic Hamiltonian to carry out both G 0 W 0 and BSE without explicitly computing empty states. We show the robustness of the approach and present the computational gains obtained on both the central processing unit and graphics processing unit nodes. In conclusion, our work will facilitate the exploration and evaluation of fine-tuned hybrid starting points aimed at enhancing the accuracy of MBPT calculations without involving computationally demanding self-consistency in Hedin’s equations.

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

Phonon Screening of Excitons in Atomically Thin Semiconductors

Atomically thin semiconductors, encompassing both 2D materials and quantum wells, exhibit a pronounced enhancement of excitonic effects due to geometric confinement. Consequently, these materials have become foundational platforms for the exploration and utilization of excitons. Recent ab initio studies have demonstrated that phonons can substantially screen electron-hole interactions in bulk semiconductors and strongly modify the properties of excitons. While excitonic properties of atomically thin semiconductors have been the subject of extensive theoretical investigations, the role of phonon screening on excitons in atomically thin structures remains unexplored. In this Letter, we demonstrate via ab initio GW-Bethe-Salpeter equation calculations that phonon screening can have a significant impact on optical excitations in atomically thin semiconductors. We further show that the degree of phonon screening can be tuned by structural engineering. We focus on atomically thin GaN quantum wells embedded in AlN and identify specific phonons in the surrounding material, AlN, that dramatically alter the lowest-lying exciton in monolayer GaN via screening. Our studies provide new intuition beyond standard models into the interplay among structural properties, phonon characteristics, and exciton properties in atomically thin semiconductors, and have implications for future experiments.

2-dimensional systems↗
DOE OSTI2025

Electronic structure and exchange interactions in altermagnetic MnGeP 2 in the quasiparticle self-consistent G W approach

The quasiparticle self-consistent 𝐺⁡𝑊 method is used to study the electronic band structure, optical dielectric function, and exchange interactions in chalcopyrite, 𝐼⁢ ‾ 4 ⁢2⁢𝑑, structure MnGeP 2 . The material is found to be an antiferromagnetic semiconductor with a lowest direct gap of 2.44 eV at the Γ point and a lower indirect gap of 1.87 eV from Γ to 𝑀. The material is an altermagnet because the two magnetic atoms of opposite spin are related by a twofold rotation operation perpendicular to the main fourfold rotation inversion axis. The spin splittings along a low symmetry line like 𝑃⁢𝑁 is sizable while at k points on the diagonal mirror planes or on the twofold symmetry axes the spin splitting is zero. The exchange interactions are calculated using a linear response approach. The antiferromagnetic exchange interaction between nearest neighbors in the primitive unit cell is dominating and found to be slightly decreasing upon carrier doping but not sufficiently to change the interaction to become ferromagnetic. The bare (noninteracting) and interacting transverse spin susceptibilities, which provide interatomic site exchange interactions after averaging over the muffin-tin spheres, are calculated from the 𝐺⁡𝑊 band structure and wave functions. From these exchange interactions, the spin wave spectra are obtained along the high symmetry lines and the Néel temperature is calculated using the mean-field and Tyablikov estimations. The dielectric function and the optical absorption spectra are calculated including excitonic effects using the Bethe Salpeter equation. The exchange interactions around Mn Ge defect sites is also studied. While we find it can generate ferromagnetic interactions with neighboring spins, we did not find direct evidence of producing an overall ferromagnetic phase. First, if Mn antisites are introduced by exchanging Mn with a nearby Ge, the interactions stay largely antiferromagnetic. Second, when we add additional Mn, in other words in Mn-rich stoichiometry, the Mn Ge antisites produce a strong ferromagnetic interaction primarily with the Mn in the same basal plane but weaker ferromagnetic interaction with adjacent plane Mn. The interactions between regular lattice Mn stay antiferromagnetic as before and thus favor keeping the antiferromagnetic order along the [001] direction. Adding Mn antisites, however, does lead to a metallic band structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗
DOE OSTI2025

𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination

Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.

Band structure methods↗
DOE OSTI2025

Five-point functions and the permutation group 𝑆 5

Five-point functions and five-body wave functions play an important role in many areas of nuclear and particle physics, e.g., in 2 →3 scattering processes, in the five-gluon vertex, or in the study of pentaquarks. In this work we consider the permutation group 𝑆 5 to facilitate the description of such objects. We work out the multiplets transforming under irreducible representations of 𝑆 5 and provide compact formulas allowing one to cast the permutations of an object 𝑓 12345 into combinations with definite permutation symmetry. We also give the explicit expressions for the irreducible multiplet products. We consider several practical applications as examples: We arrange the four-momenta and Lorentz invariants of a five-point function into the multiplet structure, we work out the color tensors of the five-gluon vertex in the multiplet notation, and we discuss applications for five-body wave functions like those of pentaquarks.

Bethe-Salpeter equation↗