Official Resources
- Homepage: https://www.vasp.at/
- Documentation: https://www.vasp.at/wiki/Bethe-Salpeter-equations_calculations
- License: Commercial (VASP license required)
Overview
The Bethe-Salpeter Equation (BSE) module in VASP provides a powerful tool for computing the frequency-dependent dielectric function with excitonic effects. The BSE can be solved on top of ground-state electronic structure from DFT, hybrid functionals, or the GW approximation, making it a state-of-the-art method for ab initio simulation of optical absorption spectra.
The BSE in VASP accounts for electron-hole interactions (excitonic effects) beyond the independent-particle and random-phase approximations. The calculation requires a preceding GW step to determine the screened Coulomb kernel W(ω→0), which is stored in temporary files (Wxxxx.tmp). VASP supports both the full BSE and the Tamm-Dancoff approximation (TDA), which neglects coupling between excitations and de-excitations for computational efficiency. The module also supports TDHF (time-dependent Hartree-Fock) calculations that do not require a preceding GW step.
Scientific domain: Optical spectroscopy, excitonic effects, many-body perturbation theory
Target user community: Researchers studying optical properties and excitons in materials
Theoretical Methods
- Bethe-Salpeter equation for electron-hole interactions
- Tamm-Dancoff approximation (TDA)
- Time-dependent Hartree-Fock (TDHF) as alternative
- Screened Coulomb interaction W from GW (static approximation)
- Bare Coulomb interaction V (exchange term)
- RPA dielectric function for screening
- Independent-particle (IP) and RPA approximations as references
- k·p perturbation theory for optical matrix elements
Capabilities (CRITICAL)
- Frequency-dependent dielectric function with excitonic effects
- Optical absorption spectra
- Exciton energies and wavefunctions
- EELS (electron energy loss spectroscopy) spectra
- BSE on top of DFT, hybrid functionals, or GW
- TDA for reduced computational cost
- TDHF without requiring GW step
- Configurable: LHARTREE, LADDER flags for different approximations
- Exciton analysis and characterization
Inputs & Outputs
Input formats:
- VASP INCAR with ALGO=BSE
- WAVECAR from ground-state calculation
- WAVEDER for optical matrix elements (LOPTICS=.TRUE.)
- Wxxxx.tmp files from preceding GW calculation (for BSE)
Output data types:
- Frequency-dependent dielectric function (real and imaginary parts)
- Optical absorption spectra
- Exciton energies and oscillator strengths
- Exciton wavefunctions and analysis
- EELS spectra
Interfaces & Ecosystem
- Programming language: Fortran (VASP core)
- Prerequisite: DFT ground-state + GW calculation (for BSE with screened W)
- Part of VASP suite: Integrated with DFT and GW modules
- Parallel computing: MPI parallelization
- Analysis: Built-in exciton analysis tools
Limitations & Known Constraints
- Commercial VASP license required
- BSE requires preceding GW step (expensive)
- Hamiltonian size scales quadratically with number of transitions
- Limited to relatively small unit cells
- Static W approximation (dynamical effects neglected)
Performance Characteristics
- BSE diagonalization scales with number of electron-hole transitions
- TDA reduces problem size by factor of ~2
- Memory-intensive for large transition spaces
- NBANDSO and NBANDSV control valence/conduction space size
- Parallel over k-points and transitions
Comparison with Other Codes
- vs EXC: EXC is a standalone BSE code; VASP-BSE is integrated
- vs Yambo: Yambo works with QE/ABINIT; VASP-BSE is self-contained
- vs BerkeleyGW: BerkeleyGW has separate BSE solver; VASP integrates everything
- vs ABINIT-BSE: ABINIT is open-source; VASP is commercial
Best Practices
- Use ALGO=BSE with preceding GW calculation for accurate screening
- Converge NBANDSO (occupied) and NBANDSV (empty) for transition space
- Use TDA for large systems to reduce computational cost
- Set LADDER=.TRUE. for excitonic effects; test with LHARTREE for exchange
- Use sufficient k-point sampling for exciton convergence
- Compare IP, RPA, and BSE spectra to assess excitonic contributions
Verification & Sources
Primary sources:
- VASP BSE documentation: https://www.vasp.at/wiki/Bethe-Salpeter-equations_calculations
- VASP BSE tutorial: https://vasp.at/tutorials/latest/bse/
- M. Gatti et al., Phys. Rev. B (various papers on VASP BSE implementation)
Confidence: VERIFIED - Official VASP documentation and tutorials confirmed