Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBethe-Salpeter Equation Calculations of Core Excitation Spectra
127
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a hybrid approach for GW/Bethe-Salpeter Equation (BSE) calculations of core excitation spectra, including x-ray absorption (XAS), electron energy loss spectra (EELS), and non-resonant inelastic x-ray scattering (NRIXS). The method is based on {it ab initio} wavefunctions from the plane-wave pseudopotential code ABINIT; atomic core-level states and projector augmented wave (PAW) transition matrix elements; the NIST core-level BSE solver; and a many-pole GW self-energy model to account for final-state broadening and self-energy shifts. Multiplet effects are also accounted for. The approach is implemented using an interface dubbed OCEAN (Obtaining Core Excitations using ABINIT and NBSE). To demonstrate the utility of the code we present results for the K-edges in LiF as probed by XAS and NRIXS, the K-edges of KCl as probed by XAS, the Ti L_2,3-edge in SrTiO_3 as probed by XAS, and the Mg L_2,3-edge in MgO as probed by XAS. We compare the results to experiments and results obtained using other theoretical approaches.
rate research
Read More
We present an efficient implementation of the Bethe-Salpeter equation (BSE) method for obtaining core-level spectra including x-ray absorption (XAS), x-ray emission (XES), and both resonant and non-resonant inelastic x-ray scattering spectra (N/RIXS). Calculations are based on density functional theory (DFT) electronic structures generated either by abinit or Quantumespresso, both plane-wave basis, pseudopotential codes. This electronic structure is improved through the inclusion of a GW self energy. The projector augmented wave technique is used to evaluate transition matrix elements between core-level and band states. Final two-particle scattering states are obtained with the NIST core-level BSE solver (NBSE). We have previously reported this implementation, which we refer to as ocean (Obtaining Core Excitations from Ab initio electronic structure and NBSE) [Phys. Rev. B 83, 115106 (2011)]. Here, we present additional efficiencies that enable us to evaluate spectra for systems ten times larger than previously possible; containing up to a few thousand electrons. These improvements include the implementation of optimal basis functions that reduce the cost of the initial DFT calculations, more complete parallelization of the screening calculation and of the action of the BSE Hamiltonian, and various memory reductions. Scaling is demonstrated on supercells of SrTiO_3 and example spectra for the organic light emitting molecule Tris-(8-hydroxyquinoline)aluminum (Alq_3 ) are presented. The ability to perform large-scale spectral calculations is particularly advantageous for investigating dilute or non-periodic systems such as doped materials, amorphous systems, or complex nano-structures.
We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation (BSE) and the calculation of the screened Coulomb interaction in finite field. The method does not require the explicit evaluation of dielectric matrices nor of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wavefunctions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.
The Bethe-Salpeter equation (BSE) is currently the state of the art in the description of neutral electron excitations in both solids and large finite systems. It is capable of accurately treating charge-transfer excitations that present difficulties for simpler approaches. We present a local basis set formulation of the BSE for molecules where the optical spectrum is computed with the iterative Haydock recursion scheme, leading to a low computational complexity and memory footprint. Using a variant of the algorithm we can go beyond the Tamm-Dancoff approximation (TDA). We rederive the recursion relations for general matrix elements of a resolvent, show how they translate into continued fractions, and study the convergence of the method with the number of recursion coefficients and the role of different terminators. Due to the locality of the basis functions the computational cost of each iteration scales asymptotically as $O(N^3)$ with the number of atoms, while the number of iterations is typically much lower than the size of the underlying electron-hole basis. In practice we see that , even for systems with thousands of orbitals, the runtime will be dominated by the $O(N^2)$ operation of applying the Coulomb kernel in the atomic orbital representation
We present first-principles many-body perturbation theory calculations of the quasiparticle electronic structure and of the optical response of HfO$_2$ polymorphs. We use the $GW$ approximation including core electrons by the projector augmented wave (PAW) method and performing a quasiparticle self-consistency also on wavefunctions (QS$GW$). In addition, we solve the Bethe-Salpeter equation on top of $GW$ to calculate optical properties including excitonic effects. For monoclinic HfO$_2$ we find a fundamental band gap of $E_g = 6.33$ eV (with the direct band gap at $E_g^d = 6.41$ eV), and an exciton binding energy of 0.57 eV, which situates the optical gap at $E^o_g = 5.85$ eV. The latter is in the range of spectroscopic ellipsometry (SE) experimental estimates (5.5-6 eV), whereas our electronic band gap is well beyond experimental photoemission (PE) estimates ($< 6$ eV) and previous $GW$ works. Our calculated density of states and optical absorption spectra compare well to raw PE and SE spectra. This suggests that our predictions of both optical and electronic gaps are close to, or at least lower bounds of, the real values.
We check the ab initio GW approximation and Bethe-Salpeter equation (BSE) many-body methodology against the exact solution benchmark of the hydrogen molecule H$_2$ ground state and excitation spectrum, and in comparison with the configuration interaction (CI) and time-dependent Hartree-Fock methods. The comparison is made on all the states we could unambiguously identify from the excitonic wave functions symmetry. At the equilibrium distance $R = 1.4 , a_0$, the GW+BSE energy levels are in good agreement with the exact results, with an accuracy of 0.1~0.2 eV. GW+BSE potential-energy curves are also in good agreement with the CI and the exact result up to $2.3 , a_0$. The solution no longer exists beyond $3.0 , a_0$ for triplets ($4.3 , a_0$ for singlets) due to instability of the ground state. We tried to improve the GW reference ground state by a renormalized random-phase approximation (r-RPA), but this did not solve the problem.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
