BSE
Parameters for a calculation solving the Bethe-Salpeter equation (BSE) for electronic excitations. The full BSE \(\left( \begin{array}{cc}A & B\\B & A\end{array} \right)\) \(\left( \begin{array}{cc}\mathbf{X}^{(n)}\\\mathbf{Y}^{(n)}\end{array} \right) = \Omega^{(n)}\left(\begin{array}{cc}1&0\\0&-1\end{array}\right)\) \(\left(\begin{array}{cc}\mathbf{X}^{(n)}\\\mathbf{Y}^{(n)}\end{array}\right)\) enables, for example, the computation of electronic excitation energies \(\Omega^{(n)}\) as well as optical properties. The BSE can be solved by diagonalizing the full ABBA-matrix or by setting B=0, i.e. within the Tamm-Dancoff approximation (TDA). Preliminary reference: Eq. (35) in PRB 92, 045209 (2015); http://dx.doi.org/10.1103/PhysRevB.92.045209 [Edit on GitHub]
Keywords
Keyword descriptions
- SECTION_PARAMETERS: logical = F
Lone keyword:
T
Usage: &BSE .TRUE.
Activates BSE calculations. [Edit on GitHub]
- BSE_DEBUG_PRINT: logical = F
Lone keyword:
T
Usage: &BSE_DEBUG_PRINT .TRUE.
Activates debug print statements in the BSE calculation. [Edit on GitHub]
- BSE_DIAG_METHOD: enum = FULLDIAG
Usage: &BSE_DIAG_METHOD FULLDIAG
Valid values:
FULLDIAG
Fully diagonalizes the BSE matrices within the chosen level of approximation.ITERDIAG
Iterative diagonalization has not been implemented yet.
Method for BSE calculations. Choose between full or iterative diagonalization. [Edit on GitHub]
- ENERGY_CUTOFF_EMPTY: real = -2.72113839E+001 [eV]
Usage: ENERGY_CUTOFF_EMPTY 10.0
Remove all orbitals with indices a,b from A_ia,jb and B_ia,jb with energy difference to LUMO level larger than the given energy cutoff, i.e. \(\varepsilon_a\in[\varepsilon_{a=\text{LUMO}}^{GW},\varepsilon_{a=\text{LUMO}}^{GW}+E_\text{cut}^\text{empty}]\). Can be used to accelerate runtime and reduce memory consumption. [Edit on GitHub]
- ENERGY_CUTOFF_OCC: real = -2.72113839E+001 [eV]
Usage: ENERGY_CUTOFF_OCC 10.0
Remove all orbitals with indices i,j from A_ia,jb and B_ia,jb with energy difference to HOMO level larger than the given energy cutoff, i.e. \(\varepsilon_i\in[\varepsilon_{i=\text{HOMO}}^{GW}-E_\text{cut}^\text{occ},\varepsilon_{i=\text{HOMO}}^{GW}]\). Can be used to accelerate runtime and reduce memory consumption. [Edit on GitHub]
- EPS_X: real = 1.00000000E-001
Usage: EPS_X 0.1
Threshold for printing contributions of singleparticle transitions, i.e. elements of the eigenvectors \(X_{ia}^{(n)}\) and \(Y_{ia}^{(n)}\). [Edit on GitHub]
- NUM_PRINT_EXC: integer = 25
Usage: NUM_PRINT_EXC 25
Number of printed excitation levels with respective energies and oscillator strengths. Does not affect computation time. [Edit on GitHub]
- NUM_PRINT_EXC_DESCR: integer = 0
Usage: NUM_PRINT_EXC_DESCR 5
Number of excitation levels for which the exciton descriptors are computed. Negative or too large NUM_PRINT_EXC_DESCR defaults to NUM_PRINT_EXC. [Edit on GitHub]
- PRINT_DIRECTIONAL_EXC_DESCR: logical = F
Lone keyword:
T
Usage: &PRINT_DIRECTIONAL_EXC_DESCR .TRUE.
Activates printing of exciton descriptors per direction. [Edit on GitHub]
- SPIN_CONFIG: enum = SINGLET
Usage: SPIN_CONFIG TRIPLET
Valid values:
SINGLET
Computes singlet excitations.TRIPLET
Computes triplet excitations.
Choose between calculation of singlet or triplet excitation (cf. given Reference above). [Edit on GitHub]
- TDA: enum = ON
Usage: &TDA ON
Valid values:
ON
The TDA is applied, i.e. B=0.OFF
The ABBA-matrix is diagonalized, i.e. the TDA is not applied.TDA+ABBA
The BSE is solved within the TDA (B=0) as well as for the full ABBA-matrix.
Level of approximation applied to BSE calculations. Choose between Tamm Dancoff approximation (TDA) and/or diagonalization of the full ABBA-matrix. [Edit on GitHub]