Computation of the expectation value of the spin operator ${\hat{S}}^2$ for the Spin-Flip Bethe-Salpeter Equation

Abstract

Spin-flip methods applied to excited-state approaches like the Bethe-Salpeter Equation allow access to the excitation energies of open-shell systems, such as molecules and defects in solids. The eigenstates of these solutions, however, are generally not eigenstates of the spin operator ${\hat{S}}^2$. Even for simple cases where the excitation vector is expected to be, for example, a triplet state, the value of $⟨{\hat{S}}^2⟩$ may be found to differ from 2.00; this difference is called ``spin contamination.'' The expectation values $⟨{\hat{S}}^2⟩$ must be computed for each excitation vector, to assist with the characterization of the particular excitation and to determine the amount of spin contamination of the state. Our aim is to provide for the first time in the spin-flip methods literature a comprehensive resource on the derivation of the formulas for $⟨{\hat{S}}^2⟩$ as well as its computational implementation. After a brief discussion of the theory of the Spin-Flip Bethe-Salpeter Equation and some examples further illustrating the need for calculating $⟨{\hat{S}}^2⟩$, we present the derivation for the general equation for computing $⟨{\hat{S}}^2⟩$ with the eigenvectors from an SF-BSE calculation, how it is implemented in a Python script, and timing information on how this calculation scales with the size of the SF-BSE Hamiltonian.

arXiv comment21 pages, 6 figures, 1 table

References in corpus (11)

• Optimized norm-conserving Vanderbilt pseudopotentials
• The PseudoDojo: Training and grading a 85 element optimized norm-conserving pseudopotential table
• BerkeleyGW: A Massively Parallel Computer Package for the Calculation of the Quasiparticle and Optical Properties of Materials and Nanostructures
• Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems
• Properties of nitrogen-vacancy centers in diamond: group theoretic approach
• Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems
• SNEG - Mathematica package for symbolic calculations with second-quantization-operator expressions
• Excited state properties of point defects in semiconductors and insulators investigated with time-dependent density functional theory
• A substitutional quantum defect in WS${}_2$ discovered by high-throughput computational screening and fabricated by site-selective STM manipulation
• First-principles calculations of the electronic structure of open-shell condensed matter systems
• Spin-Conserved and Spin-Flip Optical Excitations From the Bethe-Salpeter Equation Formalism