Bethe-Salpeter Research Articles

Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory.

Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe-Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.

Unveiling the Electronic Structure of SnS for Photocatalytic Applications: A DFT Approach†

AbstractPhotocatalytic materials are of great scientific interest due to their potential applications in sustainable energy sources. Band gap and optical absorption of the materials are two core characteristics for photocatalysis that eventually depend on electronic structure. In order to accurately compute (predict) the electronic‐scale properties of a crystalline solid using first principles calculations, which might be costly to measure experimentally, large‐scale atomic level calculation frameworks have been made possible through the use of density functional theory (DFT). We have used DFT with many body perturbation theory (MBPT) by solving the bethe‐salpeter equation (BSE) to priorly predict the properties of bulk tin sulfide (SnS) as it meets the necessary requirements for photocatalytic applications. We fixed the DFT band gap underestimation by including the Hubbard parameter (U) and used projector augmented wave pseudo‐potentials with generalized‐gradient approximation (GGA) to calculate the electronic and optical properties of SnS. Optical properties were computed using independent particle approximation (IPA) and BSE with DFT+U wave functions. Our theoretical results align with experimental data for this material, demonstrating the potential for further research in validating experimental results with theoretical approximations and empirical relations.

Significant improvement in structural, electronic, optical and thermoelectric properties of PdTe2 in bulk and monolayer phase: A G0W0+BSE approach

Using first principles calculations structural, electronic, thermoelectric and optical properties of PdTe2 in bulk and monolayer were investigated. The calculated lattice parameters with different exchange correlation functionals show that the results obtained with addition of van der Waals corrections (vdW-DFC09x) are very well matched with the results obtained experimental. For correct prediction of electronic band gap, G0W0 approximation is used and the calculated band structure reveal that PdTe2 in bulk form has metallic nature as attained in experimental studies. In addition, the electronic band gap results with G0W0 approximation depicted that PdTe2 in monolayer phase is a semiconductor material with a direct band gap of 1.12 eV. The optical spectra have been calculated via many-body perturbation theory (MBPT) by solving Bethe-Salpeter equation (BSE) (G0W0+BSE). Ultimately, the results of optical spectra demonstrated that PdTe2 in monolayer phase has strong optical absorption within visible light wavelengths than the bulk phase which is good for solar cell applications. Thermoelectric (TE) properties of PdTe2 in bulk and monolayer structures have been studied using BoltzTraP code. The TE properties calculations indicated that better performance can be obtained by reducing the dimensionality of materials. The figure of merit was found to be above unity for monolayer phase. This paper successfully exploited the enhancement of the figure of merit values with the monolayer structure of PdTe2 material.

Tailoring the optoelectronic properties of MoS2 for broadband photodetection: Showcasing an Ab-into study involving the quasi-particle correction within the Green’s function-based approximation

In this study, the electronic and structural characteristics of both bulk and monolayer hexagonal MoS2 material are studied within and beyond density functional theory (DFT) by considering one-shot GW corrections for an accurate band gap prediction. Interestingly, our computed quasiparticle (QP) bandgap of 1.36 eV and 3.33 eV for bulk and monolayer, respectively, is in excellent agreement with the experimental finding. The elucidation of the optical properties was estimated with the aid of the Bethe–Salpeter equation within the many-body perturbation theory (MPB). The computed optical transitions such as absorption coefficient, extinction coefficient, refractive index, reflectivity, and electron energy loss function are in better agreement with the experimental data. Moreover, the new approach allows for the efficient inclusion of all conduction bands in GW calculations, significantly reducing computational costs compared to traditional methods. A good agreement of the lattice parameters, as well as the interlayer distance with the experimental findings, is observed by the inclusion of van der Waals (vdW) corrections. Besides the agreement with previously reported values, the evaluated band structure determined with the quasiparticle corrections shows that both bulk and monolayer MoS2 material exhibit semiconducting properties with direct and indirect band gaps, respectively. The calculation of the density of states reveals that the s-orbitals of S atoms and the d-orbitals of Mo are predominantly at the origin of the essential material properties secured at the Fermi level. The results of the study also point to the possibility of a foundation for creating physical structures with swift, decisive, directed, and long-range broadband photodetection capabilities.

Multiple scattering, radiative transfer, and parameter regimes

We have established the radiative transfer regime in a four-parameter space spanned by amplitude and scale of fluctuations, convergence rate, and spatial range. The derivation of radiative transfer equations from wave equations offers a better understanding of the underlying physics in radiative transfer theory. One approach is based on tools borrowed from quantum field theory (QFT) and employs the Dyson and Bethe–Salpeter equations as the foundation. This has allowed us to investigate the small scale fluctuation regime, kl = O ( ε ) . We also investigated this regime using the scale-separated asymptotics approach and found that the conditions on the parameters are the same as in the QFT approach if the convergence rate of the Dyson series is O ( ε 2 ) . Furthermore, we found that there is a multitude of regimes where radiative transfer manifests. We hence infer that for a fixed spatial range, increase in the scale size of fluctuations leads to decrease in the upper bound of variance of permittivity fluctuations along with diminished convergence rate of the Wigner series. Also, for a fixed scale size and convergence rate, the spatial range where radiative transport manifests decreases with increase in strength of fluctuations.

Dual Bethe-Salpeter equation for the multiorbital lattice susceptibility within dynamical mean-field theory

Dual Bethe-Salpeter equation for the multiorbital lattice susceptibility within dynamical mean-field theory

Elucidating Solvatochromic Shifts in Two-Dimensional Photocatalysts by Solving the Bethe-Salpeter Equation Coupled with Implicit Solvation Method.

Many studies have focused on tailoring the photophysical properties of two-dimensional (2D) materials for photocatalytic (PC) or photoelectrochemical (PEC) applications. To understand the optical properties of 2D materials in solution, we established a computational method that combined the Bethe-Salpeter equation (BSE) calculations with our GW-GPE method, allowing for GW/BSE-level calculations with implicit solvation described using the generalized Poisson equation (GPE). We applied this method to MoS2, phosphorene (PP), and g-C3N4 and found that when the solvent dielectric increased, it reduced the exciton binding energy and quasiparticle bandgap, resulting in almost no solvatochromic shift in the excitonic peaks of MoS2 and PP, which is consistent with previous experiments. However, our calculations predicted that the solvent dielectric had a significant impact on the excitonic properties of g-C3N4, exhibiting a large solvatochromic shift. We expect that our GW/BSE-GPE method will offer insights into the design of 2D materials for PC and PEC applications.

Weakly Confined Organic-Inorganic Halide Perovskite Quantum Dots as High-Purity Room-Temperature Single Photon Sources.

Colloidal perovskite quantum dots (PQDs) have emerged as highly promising single photon emitters for quantum information applications. Presently, most strategies have focused on leveraging quantum confinement to increase the nonradiative Auger recombination (AR) rate to enhance single-photon (SP) purity in all-inorganic CsPbBr3 QDs. However, this also increases the fluorescence intermittency. Achieving high SP purity and blinking mitigation simultaneously remains a significant challenge. Here, we transcend this limitation with room-temperature synthesized weakly confined hybrid organic-inorganic perovskite (HOIP) QDs. Superior single photon purity with a low g(2)(0) < 0.07 ± 0.03 and a nearly blinking-free behavior (ON-state fraction >95%) in 11 nm FAPbBr3 QDs are achieved at room temperature, attributed to their long exciton lifetimes (τX) and short biexciton lifetimes (τXX). The significance of the organic A-cation is further validated using the mixed-cation FAxCs1-xPbBr3. Theoretical calculations utilizing a combination of the Bethe-Salpeter (BSE) and k·p approaches point toward the modulation of the dielectric constants by the organic cations. Importantly, our findings provide valuable insights into an additional lever for engineering facile-synthesized room-temperature PQD single photon sources.

Many-body theory calculations of positron binding to hydrogen cyanide

Positron bound state properties in hydrogen cyanide are studied via many-body theory calculations that account for strong positron-electron correlations including positron-induced polarization, screening of the electron–positron Coulomb interaction, virtual-positronium formation and positron–hole repulsion. Specifically, the Dyson equation is solved using a Gaussian basis, with the positron self-energy in the field of the molecule calculated using the Bethe–Salpeter equations for the two-particle and particle–hole propagators. The present results suggest near cancellation of screening corrections to the bare polarization, and the non-negligible role of the positron–hole interaction. There are no existing measurements to compare to for HCN. Previous configuration interaction (CI) and fixed-node diffusion Monte Carlo (FN-DMC) calculations give positron binding energies in the range 35–44 meV, most of which used a single even-tempered basis centred near the nitrogen atom. Using a similar single-centre positron basis we calculate a positron binding energy of 41 meV, in good agreement. However, we find that including additional basis centres gives an improved description of the positron wave function near the nuclei, and results in a converged binding energy in the range 63–73 meV (depending on geometry and approximation to the positron–molecule correlation potential used).Graphical

Excitons, optical spectra, and electronic properties of semiconducting Hf-based MXenes.

Semiconducting MXenes are an intriguing two-dimensional (2D) material class with promising electronic and optoelectronic properties. Here, we focused on recently prepared Hf-based MXenes, namely, Hf3C2O2 and Hf2CO2. Using the first-principles calculation and excited state corrections, we proved their dynamical stability, reconciled their semiconducting behavior, and obtained fundamental gaps by using the many-body GW method (indirect 1.1 and 2.2eV; direct 1.4 and 3.5eV). Using the Bethe-Salpeter equation, we subsequently provided optical gaps (0.9 and 2.7eV, respectively), exciton binding energies, absorption spectra, and other properties of excitons in both Hf-based MXenes. The indirect character of both 2D materials further allowed for a significant decrease of excitation energies by considering indirect excitons with exciton momentum along the Γ-M path in the Brillouin zone. The first bright excitons are strongly delocalized in real space while contributed by only a limited number of electron-hole pairs around the M point in the k-space from the valence and conduction band. A diverse range of excitonic states in Hf3C2O2 MXene lead to a 4% and 13% absorptance for the first and second peaks in the infrared region of absorption spectra, respectively. In contrast, a prominent 28% absorptance peak in the visible region appears in Hf2CO2 MXene. Results from radiative lifetime calculations indicate the promising potential of these materials in optoelectric devices requiring sustained and efficient exciton behavior.

Efficient Full-Frequency GW Calculations Using a Lanczos Method.

The GW approximation is widely used for reliable and accurate modeling of single-particle excitations. It also serves as a starting point for many theoretical methods, such as its use in the Bethe-Salpeter equation (BSE) and dynamical mean-field theory. However, full-frequency GW calculations for large systems with hundreds of atoms remain computationally challenging, even after years of efforts to reduce the prefactor and improve scaling. We propose a method that reformulates the correlation part of the GW self-energy as a resolvent of a Hermitian matrix, which can be efficiently and accurately computed using the standard Lanczos method. This method enables full-frequency GW calculations of material systems with a few hundred atoms on a single computing workstation. We further demonstrate the efficiency of the method by calculating the defect-state energies of silicon quantum dots with diameters up to 4nm and nearly 2,000 silicon atoms using only 20 computational nodes.

GPAW: An open Python package for electronic structure calculations.

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Accelerating Analytic-Continuation GW Calculations with a Laplace Transform and Natural Auxiliary Functions.

We present a simple and accurate GW implementation based on a combination of a Laplace transform (LT) and other acceleration techniques used in post-self-consistent field quantum chemistry, namely, natural auxiliary functions and the frozen-core approximation. The LT-GW approach combines three major benefits: (a) a small prefactor for computational scaling, (b) easy integration into existing molecular GW implementations, and (c) significant performance improvements for a wide range of possible applications. Illustrating these advantages for systems consisting of up to 352 atoms and 7412 basis functions, we further demonstrate the benefits of this approach combined with an efficient implementation of the Bethe-Salpeter equation.

Adsorption of cadmium selenide clusters: A novel approach to enhance solar energy conversion using armchair graphene nanoribbons

Tailoring the electronic, optical, and transport properties of low-dimensional semiconductor materials is essential to improve the light-conversion efficiency of thin-film solar cell materials. Here, using first-principles calculations and non-equilibrium Green functions, we investigate the enhancement of optoelectronic and transport properties of armchair graphene nanoribbons (AGNRs) upon adsorption of cadmium selenide clusters. Upon adsorption of a CdSe diatomic molecule on an AGNR, the most energetically favorable configuration is the cadmium end sitting on top of a carbon atom. The corresponding electronic bandgap reduces ∼5 times with respect to that of the pristine system, thanks to the formation of a polaron state formed by the p-orbital of the selenide atom. Upon adsorption of CdSe cyclohexane molecules, the bandgap of this system slightly shrinks by 0.121 eV with respect to the pristine system. The charge accumulation induced by these clusters significantly enhances the absorption coefficient of the adsorbed systems, resulting in a red shift of the optical spectra toward the infrared region. More interestingly, by solving the Bethe–Salpeter equations with the Tamm–Dancoff approximation, we provide a direct link between the first-principles optical prediction and experimental observations. In addition, the electron transfer from these molecules to the hosted systems increases the transmission spectra in the vicinity of the Fermi level, leading to a remarkable electronic current passing through these scattering regions. These results highlight the role of cadmium selenide clusters in enhancing the light-to-energy conversion efficiency of next-generation solar cell devices.

Absolute energy levels of liquid water from many-body perturbation theory with effective vertex corrections

We demonstrate the importance of addressing the Γ vertex and thus going beyond the GW approximation for achieving the energy levels of liquid water in many-body perturbation theory. In particular, we consider an effective vertex function in both the polarizability and the self-energy, which does not produce any computational overhead compared with the GW approximation. We yield the band gap, the ionization potential, and the electron affinity in good agreement with experiment and with a hybrid functional description. The achieved electronic structure and dielectric screening further lead to a good description of the optical absorption spectrum, as obtained through the solution of the Bethe-Salpeter equation. In particular, the experimental peak position of the exciton is accurately reproduced.

Two‐Time Quantum Fluctuations Approach and Its Relation to the Bethe–Salpeter Equation

Correlated quantum many‐particle systems out of equilibrium are of high interest in many fields, including correlated solids, ultracold atoms, or dense plasmas. Accurate theoretical description of these systems is challenging both, conceptionally and with respect to computational resources. A quantum fluctuations approach is recently presented, which is equivalent to the nonequilibrium GW approximation that promises high accuracy at low computational cost. The method exhibits process time scaling that is linear in the number of time steps, like the G1–G2 scheme, however, at a much reduced computer memory cost. In a second publication, this approach is extended to the two‐time exchange–correlation functions and the dynamic density response properties. Herein, the properties of this approach are analyzed in more detail. The physical meaning of the central approximation, the quantum polarization approximation, is established. It is demonstrated that the method is equivalent to the Bethe–Salpeter equation for the two‐time exchange–correlation function when the generalized Kadanoff–Baym ansatz with Hartree–Fock propagators is applied.

Nonequilibrium Dynamics of the Hubbard Dimer

Electron dynamics in a two‐sites Hubbard model is studied using the nonequilibrium Green's function approach. The study is motivated by the empirical observation that a full solution of the integro‐differential Kadanoff–Baym equation (KBE) is more stable and often accompanied by artificial damping [Marc Puig von Friesen, C. Verdozzi, and C.‐O. Almbladh (2009)] than its time‐linear reformulations relying on the generalized Kadanoff–Baym ansatz (GKBA). Additionally, for conserving theories, numerical simulations suggest that KBE produces natural occupations bounded by one and zero in agreement with the Pauli exclusion principle, whereas, in some regimes, GKBA‐based theories violate this principle. As the first step for understanding these issues, the electron dynamics arising in the adiabatic switching scenario is studied. Many‐body approximations are classified according to the channel of the Bethe–Salpeter equation in which electronic correlations are explicitly treated. They give rise to the so‐called second Born, T‐matrix, and GW approximations. In each of these cases, the model is reduced to a system of ordinary differential equations, which resemble equations of motion for a driven harmonic oscillator with time‐dependent frequencies. A more complete treatment of electronic correlations is achieved by combining different correlation channels, with parquet theory serving as a starting point.

Spinel LiGa5O8 prospects as ultra-wideband-gap semiconductor: Band structure, optical properties, and doping

LiGa5O8 in the spinel type structure is investigated as a potential ultra-wideband-gap semiconductor. The band structure is determined using the quasiparticle self-consistent GW method, and the optical properties are calculated at the Bethe Salpeter Equation level including electron-hole interaction effects. The optical gap including exciton effects and an estimate of the zero-point motion electron phonon coupling renormalizations is estimated to be about 5.2±0.1 eV with an exciton binding energy of about 0.4 eV. Si doping as potential n-type dopant is investigated and found to be a promising shallow donor.

Non-linear light–matter interactions from the Bethe–Salpeter equation

A route to assess non-linear light–matter interactions from the increasingly popular GW-Bethe–Salpeter equation (GW-BSE) method is outlined. In the present work, the necessary analytic expressions within the static-screened exchange approximation of the BSE are derived. This enables a straightforward implementation of the computation of the first hyperpolarizability as well as two-photon absorption processes for molecular systems. Benchmark calculations on small molecular systems reveal that the GW-BSE method is intriguingly accurate for predicting both first hyperpolarizabilities and two-photon absorption strengths. Using state-of-the-art Kohn–Sham references as a starting point, the accuracy of the GW-BSE method rivals that of the coupled-cluster singles-and-doubles method, outperforming both second-order coupled-cluster and time-dependent density-functional theory.

Solving the Bethe-Salpeter equation in real frequencies at finite temperature

Solving the Bethe-Salpeter equation in real frequencies at finite temperature