Pith. sign in

REVIEW 7 minor 102 references

Analytical Forces from the Bethe-Salpeter Equation for Large-Scale Excited-State Relaxation

T0 review · 0 major / 7 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read An analytical-force implementation of the Bethe-Salpeter equation replaces one response calculation per atomic motion with a single Z-vector equation, bringing excited-state relaxation to defect supercells of hundreds of atoms.

desk verdict Solid methods letter: plane-wave BSE forces with an honest, well-tested scissor/frozen-W scope; the fixed-scissor caveat is real but the paper does the work to show when it holds. read the letter →

arxiv 2607.20728 v1 pith:TUSVTGQO submitted 2026-07-22 cond-mat.mtrl-sci

Yu Jin , Victor Wen-zhe Yu , Marco Govoni , Giulia Galli This is my paper · ORCID
classification cond-mat.mtrl-sci PACS 71.15.-m71.35.-y
keywords Bethe-Salpeterequationexcited-stateforcesanalyticalnucleargradientsZ-vectordensity-matrixperturbationtheoryprojectivedielectriceigenpotentialpointdefectsinsemiconductorsphotoluminescencelineshapes
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper's goal is to make structural relaxation of electronic excited states affordable at the Bethe-Salpeter equation (BSE) level for solids, by computing analytical forces rather than differentiating energies numerically. The key algorithmic claim is that a Z-vector Lagrangian combined with a low-rank representation of the screened Coulomb interaction reduces the cost from one response calculation per nuclear displacement to a single equation, with no explicit empty-state sums, so that GPU-accelerated supercells of several hundred atoms become feasible. On two point defects, the paper argues when that machinery changes the answer: for the nitrogen-vacancy center in diamond, BSE and hybrid-functional TDDFT agree closely, while for the carbon-dimer defect in two-dimensional h-BN, only the screened electron-hole attraction in BSE stabilizes the localized defect excitation and corrects a relaxation pattern that semilocal TDDFT gets qualitatively wrong. A reader would care because these forces control photoluminescence line shapes, self-trapped excitons, and the behavior of defect-based quantum emitters — quantities previously out of reach of BSE-level theory for realistic supercells.

What carries the argument

Central object: the Z-vector (Handy-Schaefer) response equation — one Sternheimer-like solve for orbitals Z_v whose source U is the derivative of the BSE excitation energy with respect to occupied KS orbitals. It yields the relaxed density Δρ(z); with the unrelaxed Δρ(x) from BSE amplitudes, dω/dR = ∫ (∂V_ext/∂R)(Δρ(x)+Δρ(z)), replacing 3N per-atom DFPT solves. Companion mechanism: the PDEP (projective dielectric eigenpotential) low-rank form of the screened interaction W, avoiding full dielectric matrices and empty-state sums. Two explicit approximations carry it: a fixed scissor shift in place of the QP-Hamiltonian derivative, and δW/δ⟨φ_v|≈0. Wannier localization, an inexact Krylov solver

What would settle it

Relax the carbon-dimer defect in 2D h-BN from a PBE starting point using fully analytical BSE-G0W0 forces, available in Gaussian-basis implementations, and compare with the fixed-scissor BSE relaxation: the paper's Fig. 3 already shows the two potential-energy curves diverging, so the G0W0-relaxed geometry would differ measurably in bond lengths and photoluminescence sideband. Applying the same comparison to any system with a strongly geometry-dependent self-energy correction at a hybrid starting point would delimit where the approximations bind.

Watch

Extended reading notes

Core claim

Central claim: the derivative of a BSE excitation energy is the integral of the external-potential derivative against a differential density built from BSE response orbitals plus Z-vector orbitals from one Handy-Schaefer response equation, at a cost only a few times that of TDDFT forces. After validation, the method relaxes the triplet 3E state of NV- in a 511-atom diamond supercell and the first singlet state of the CBCN defect in a 288-atom 2D-hBN supercell. For NV-, BSE and TDDFT give nearly identical displacements and photoluminescence spectra; for CBCN, only BSE with a dielectric-dependent hybrid starting point reproduces experiment, the screened electron-hole interaction keeping the re

Load-bearing premise

The load-bearing premise is that the quasiparticle correction behaves as a fixed, geometry-independent scissor shift and that screening barely changes when occupied orbitals move; where this fails, the forces are not derivatives of the actual BSE energy surface — a breakdown the paper's own Fig. 3 shows for the PBE starting point in 2D h-BN.

Editorial extensions

If this is right

  • Excited-state geometry relaxation at the BSE level becomes routine for defect supercells of hundreds of atoms; the reported cost is only a few times that of hybrid-functional TDDFT forces (hundreds of GPU-hours per force evaluation).
  • For localized deep defects in quasi-homogeneous hosts such as NV- in diamond, hybrid-functional TDDFT is a reliable and much cheaper surrogate: BSE and TDDFT agree on displacements and photoluminescence line shapes within the starting-point dependence.
  • For defects in inhomogeneous dielectric environments such as the carbon dimer in 2D h-BN, semilocal TDDFT relaxes toward an excited state of the wrong character; the screened electron-hole attraction of BSE is required to keep the excitation localized.
  • The fixed-scissor approximation used in the forces is reliable only when the starting point leaves a small, geometry-independent quasiparticle correction; a dielectric-dependent hybrid starting point restores that condition, whereas a semilocal starting point does not (the paper's Fig. 3).

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same Lagrangian machinery should extend to analytical non-adiabatic couplings and excited-state vibrational line widths, since the Z-vector source U already contains the orbital-response information those derivatives require; the paper lists both as future work.
  • A practical screening rule follows from the two defect case studies: start with hybrid-functional TDDFT forces and escalate to BSE forces only when the dielectric environment is strongly inhomogeneous or the excited state delocalizes — a policy directly applicable to computational searches for quantum emitters.
  • The paper's Fig. 3 pattern implies a precise prediction: the scissor approximation will fail exactly where the self-energy correction is large and geometry-dependent, for example a defect near an interface or surface even at a hybrid starting point; testing such a system would sharpen the method's domain of validity.
  • Because the frozen-W approximation fixes screening at the ground-state geometry, the forces inherit a residual starting-point dependence that only a self-consistent GW workflow would remove; until then, differences in relaxed geometry between starting points may rival the differences between TDDFT and BSE themselves.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

0 major / 7 minor

Summary. This Letter reports a plane-wave implementation of analytical nuclear forces for BSE excited states, built on density-matrix perturbation theory, a Z-vector Lagrangian, the projective dielectric eigenpotential (PDEP) representation of screening, and GPU acceleration. The method requires a single response equation rather than 3N_atom Sternheimer solves, and it avoids explicit empty-state sums. The authors validate the implementation against finite-difference forces for the NV- center in diamond with PBE, DDH, and HSE starting points, against small-molecule BSE geometries from a fully analytical molecular implementation, and against finite-difference derivatives for CO under the same approximations. They then apply the method to NV- in diamond and the carbon-dimer defect C_BC_N in 2D h-BN, comparing TDDFT and BSE relaxed geometries and photoluminescence spectra. The two central approximations—a fixed scissor for the quasiparticle Hamiltonian derivative and neglect of the orbital derivative of W—are stated explicitly and are the subject of the CO test and the PEC comparison in Fig. 3.

Significance. The paper makes a substantial methodological contribution: it demonstrates a scalable path to BSE-level excited-state geometry relaxation in condensed systems, a capability that has so far been limited to small molecules. The implementation is carefully validated, with analytical finite-difference force agreement at the ~10^-5 Ry/bohr level, molecular benchmarks against fully analytical BSE-G0W0 forces, and an explicit test of the fixed-scissor approximation in the end matter and Fig. 3. The two applications are meaningful: for NV- in diamond the comparison with hybrid-TDDFT delineates where simpler methods suffice, while for the 2D h-BN defect the BSE result identifies the qualitative failure of semilocal TDDFT and the importance of screened electron-hole attraction. The authors are transparent about the surrogate nature of the forces—they are derivatives of a BSE Hamiltonian with a fixed scissor and frozen W, not of the full BSE-G0W0 functional—and they provide evidence on the conditions under which this surrogate is reliable.

minor comments (7)
  1. [Title/Abstract] The title and abstract describe the method as 'analytical forces from the Bethe-Salpeter equation' without qualification. Since the forces are derivatives of a BSE Hamiltonian with a fixed-scissor approximation for QP derivatives and with δW/δ<φ_v| set to zero, I recommend adding a qualifier such as 'with fixed-scissor/frozen-W approximations' in the abstract (and at least a clarifying sentence in the introduction) to avoid readers interpreting the forces as full BSE-G0W0 gradients. The paper itself already states the scope in the theoretical framework, but the abstract currently overstates it.
  2. [Theoretical framework, Eq. (4)-(7)] The notation dω/dR in Eq. (4) should be explicitly defined as the derivative at fixed KS orbitals; otherwise the reader may not immediately see that the Z-vector formalism is needed because the orbitals are geometry-dependent. A one-sentence clarification would help.
  3. [Tables I and II] The BSE VEEs are fitted by construction: the scissor is chosen to reproduce the BSE-G0W0 VEE at the ground-state geometry. The AEEs are therefore not fully independent predictions. The text notes this for VEEs, but it would be useful to add a footnote to both tables stating that the BSE AEEs inherit the geometry-independent scissor fitted at the GS geometry, so the apparent agreement with experiment is a combined test of the BSE model and the scissor approximation.
  4. [Fig. 3 and discussion] The potential energy curve test of the fixed-scissor approximation is performed for only one interpolation coordinate in a 128-atom supercell, while the production CBCN calculations use a 288-atom supercell. The authors should state whether the smaller cell changes the conclusions, and ideally report the maximum deviation between fixed-scissor and BSE-G0W0 forces or AEEs along the tested path. This would make the 'reliable for DDHα' conclusion more quantitative.
  5. [SM S4] There is a typo: 'for the save of computational cost' should read 'for the sake of computational cost.'
  6. [Table SIII and main text] The sentence 'at a cost only a few times larger than the corresponding TDDFT forces calculations for hybrid-functional starting points' is supported for DDH (~9x for NV-, ~4x for hBN) but not for PBE (~73x for NV-). Since the sentence already limits the claim to hybrid starting points, it may be worth adding the PBE comparison explicitly to avoid an unintended implication that BSE-PBE is similarly close to TDDFT-PBE.
  7. [SM S5, Fig. S2/S3] The statement 'consistent with the absence of an explicit screened electron-hole attraction in TDDFT with semilocal functionals' is useful. Consider also stating that the DDHα TDDFT result includes an analogous effective attraction through the exact-exchange term, which is why TDDFT-DDHα and BSE-DDHα differential densities are nearly identical. This helps the reader understand why the starting point, not just the BSE kernel, is the controlling factor.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the BSE-force derivation is self-contained, validated against finite differences and external benchmarks; the fixed-scissor calibration is a disclosed approximation, not a hidden fit of the predicted relaxation.

full rationale

The paper's central derivation is not circular. The analytical BSE forces follow from a Z-vector Lagrangian (Eqs. 4-7), and the implementation is checked against finite-difference forces at the same level of approximation (SM Table SII) and against fully analytical BSE-G0W0 and CC3/CCSDR(3) geometries for small molecules (SM S3). The fixed-scissor operator is fitted only to the vertical excitation energy at the ground-state geometry, and the paper explicitly states this calibration: 'The BSE excitation energies reported here are obtained using a fixed scissor operator chosen to reproduce the VEE from the BSE-G0W0 calculation at the GS geometry.' The relaxed geometries, AEEs, and PL line shapes are not fit to this quantity; they are computed from forces and are independently tested against BSE-G0W0 potential-energy curves for CO (End Matter) and for the CBCN defect in hBN (Fig. 3), which also reveals where the approximation fails. The two approximations used in the force expression (delta H^QP ~ delta H^KS, delta W/delta <phi_v| = 0) are stated, not hidden, and their geometry-dependence is assessed. Self-citations to WEST, PDEP, and DDH functionals refer to established computational methods and nonempirical functionals; they are not invoked as a uniqueness theorem or to forbid alternatives. Therefore no load-bearing step in the derivation reduces to its own input by construction.

Assumptions & free parameters 1 free parameters · 6 assumptions · 0 invented entities

The method relies on established DFT/BSE formalism and two explicit approximations: fixed-scissor QP derivative and frozen W. No new physical entities are introduced. The scissor shift is the only free parameter fitted by the authors; the DDH EXX fractions are inherited from prior parameterization conditions.

free parameters (1)
  • Fixed scissor shift (Δsci) = NV− diamond: 1.17 eV (PBE), 0.45 eV (DDH); CBCN 2D h-BN: 2.77 eV (PBE), 1.01 eV (DDHα)
    Fitted per system and starting point so that the BSE vertical excitation energy at the ground-state geometry matches the BSE-G0W0 reference; this parameter enters the force calculation via the fixed-scissor approximation.
assumptions (6)
  • domain assumption BSE in the Tamm-Dancoff approximation (TDA), B_s = 0
    Used in Eq. (1) of the main text; standard approximation for neutral excitations in solids.
  • domain assumption Statically screened Coulomb interaction W evaluated in RPA using PDEP
    Eq. (3); W is static and constructed from the leading eigenpairs of the symmetrized density-density response function.
  • ad hoc to paper Fixed-scissor approximation for the QP Hamiltonian derivative: δH_QP/δ⟨φv| ≈ δH_KS/δ⟨φv|
    Stated in the Theoretical Framework: the derivative of the QP Hamiltonian is replaced by that of the KS Hamiltonian plus a constant shift. The paper validates this for CO and via potential-energy-curve comparisons.
  • ad hoc to paper Neglect of δW/δ⟨φv| ≈ 0
    Explicitly set to zero in the derivation of Eq. (S12). This is the main approximation; the paper assesses its effect on CO and in Fig. 3 for 2D h-BN.
  • domain assumption KS orbitals from semi-local or hybrid DFT starting points are not further optimized at the BSE level
    The BSE eigenvalue problem uses KS orbitals; the relaxation results depend on the chosen starting functional, which the paper discusses as a source of starting-point dependence.
  • domain assumption PBE frozen-phonon vibrational modes are used for all PL spectra regardless of functional
    Supplemental Material S4: vibrational modes are computed at the PBE level for computational cost, introducing an approximation in the Huang-Rhys line shapes.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Analytical Forces from the Bethe-Salpeter Equation for Large-Scale Excited-State Relaxation." pith.science (2026). https://pith.science/paper/TUSVTGQO

@misc{pith2026260720728,
  author       = {Pith},
  title        = {Pith review of: Analytical Forces from the Bethe-Salpeter Equation for Large-Scale Excited-State Relaxation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TUSVTGQO}},
  note         = {Machine review of arXiv:2607.20728}
}
read the original abstract

We present an efficient plane-wave implementation of analytical nuclear forces for electronic excited states described by the Bethe-Salpeter equation (BSE). The formulation combines density-matrix perturbation theory with a Lagrangian approach, and avoids both explicit empty-state summations and the response calculations for each atomic displacement, required by conventional approaches based on density functional perturbation theory. Together with GPU acceleration, these advances make BSE forces calculations tractable for solid-state systems containing hundreds of atoms. We demonstrate the method on two point defects with distinct dielectric environments: the nitrogen-vacancy center in diamond, where BSE and time-dependent density functional theory (TDDFT) yield consistent excited-state relaxations, and the carbon-dimer defect in two-dimensional hexagonal boron nitride, where the screened electron-hole interaction included in the BSE stabilizes the localized defect excitation and corrects the relaxation pattern predicted by semilocal TDDFT. These results establish a scalable framework for BSE-level studies of excited-state relaxation and vibronic coupling in heterogeneous condensed systems.

Figures

Figures reproduced from arXiv: 2607.20728 by the authors.

Figure 1
Figure 1. FIG. 1. Excited-state relaxation of NV [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Excited-state relaxation of C [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Potential energy curves computed with different [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Validation of BSE forces for the lowest singlet excited [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

102 extracted references · 3 linked inside Pith

  1. [1]

    Huang, A

    K. Huang, A. Rhys, and N. F. Mott, Theory of light absorption and non-radiative transitions inF-centres, Proc. R. Soc. Lond. A204, 406 (1950)

    1950

  2. [2]

    Alkauskas, B

    A. Alkauskas, B. B. Buckley, D. D. Awschalom, and C. G. Van de Walle, First-principles theory of the lumi- nescence lineshape for the triplet transition in diamond nv centres, New J. Phys.16, 073026 (2014)

    2014

  3. [3]

    Alkauskas, M

    A. Alkauskas, M. D. McCluskey, and C. G. Van de Walle, Tutorial: Defects in semiconductors—combining experiment and theory, J. Appl. Phys.119(2016)

    2016

  4. [4]

    S. Li, J. Luo, J. Liu, and J. Tang, Self-trapped exci- tons in all-inorganic halide perovskites: fundamentals, status, and potential applications, J. Phys. Chem. Lett. 10, 1999 (2019)

    1999

  5. [5]

    Z. Li, J. Xie, M. Faizan, S. Wang, X. Chen, and L. Zhang, Luminescence from localized states in solids: a first-principles perspective, Adv. Funct. Mater.35, 2422716 (2025)

    2025

  6. [6]

    D. K. Lewis and S. Sharifzadeh, Modeling excited states of point defects in materials from many-body perturba- tion theory, ACS Materials Lett.3, 862 (2021)

    2021

  7. [7]

    H.-Y. Sun, L. Xiong, and H. Jiang, Toward first- principles approaches for mechanistic study of self- trapped exciton luminescence, Chem. Phys. Rev.4, 031302 (2023)

    2023

  8. [8]

    H. Seo, V. Iv´ ady, and Y. Ping, First-principles computa- tional methods for quantum defects in two-dimensional materials: A perspective, Appl. Phys. Lett.125, 140501 (2024)

    2024

Show all 102 references
  1. [9]

    Gali, Recent advances in theab initiotheory of solid- state defect qubits, Nanophotonics12, 359 (2023)

    ´A. Gali, Recent advances in theab initiotheory of solid- state defect qubits, Nanophotonics12, 359 (2023)

    2023

  2. [10]

    Ping and T

    Y. Ping and T. J. Smart, Computational design of quan- tum defects in two-dimensional materials, Nat. Comput. Sci.1, 646 (2021)

    2021

  3. [11]

    Xiong, Y

    Y. Xiong, Y. Zhu, S. M. Griffin, and G. Hautier, Finding the perfect imperfection: Accelerated, computationally driven discovery and design of quantum defects, MRS Bull.51, 276– (2026)

    2026

  4. [12]

    A. Gali, E. Janz´ en, P. De´ ak, G. Kresse, and E. Kaxi- ras, Theory of spin-conserving excitation of then−V − center in diamond, Phys. Rev. Lett.103, 186404 (2009)

    2009

  5. [13]

    X. Wang, W. Meng, W. Liao, J. Wang, R.-G. Xiong, and Y. Yan, Atomistic mechanism of broadband emission in metal halide perovskites, J. Phys. Chem. Lett.10, 501 (2019)

    2019

  6. [14]

    Y. Jin, M. Govoni, and G. Galli, Vibrationally resolved optical excitations of the nitrogen-vacancy center in di- amond, npj Comput. Mater.8, 238 (2022)

    2022

  7. [15]

    Y. Jin, V. W.-z. Yu, M. Govoni, A. C. Xu, and G. Galli, Excited state properties of point defects in semiconduc- tors and insulators investigated with time-dependent density functional theory, J. Chem. Theory Comput.19, 8689 (2023)

    2023

  8. [16]

    Botti, A

    S. Botti, A. Schindlmayr, R. Del Sole, and L. Reining, Time-dependent density-functional theory for extended systems, Rep. Prog. Phys.70, 357 (2007)

    2007

  9. [17]

    Sun, S.-X

    H.-Y. Sun, S.-X. Li, and H. Jiang, Pros and cons of the time-dependent hybrid density functional approach for calculating the optical spectra of solids: a case study of ceo 2, Phys. Chem. Chem. Phys.23, 16296 (2021)

    2021

  10. [18]

    Xiong and G

    Y. Xiong and G. Hautier, ∆SCF inV ASPfor excited-state defect computations: tips and pitfalls, arXiv:2505.04748 10.48550/arXiv.2505.04748 (2025)

    2025 doi arXiv

  11. [19]

    Ismail-Beigi and S

    S. Ismail-Beigi and S. G. Louie, Excited-state forces within a first-principles green’s function formalism, Phys. Rev. Lett.90, 076401 (2003)

    2003

  12. [20]

    Ismail-Beigi and S

    S. Ismail-Beigi and S. G. Louie, Self-trapped excitons in silicon dioxide: Mechanism and properties, Phys. Rev. Lett.95, 156401 (2005)

    2005

  13. [21]

    R. R. Del Grande and D. A. Strubbe, Revisiting ab- initio excited state forces from many-body green’s func- tion formalism: approximations and benchmark, arXiv preprint arXiv:2502.05144 (2025)

    2025

  14. [22]

    Villalobos-Castro, I

    J. Villalobos-Castro, I. Knysh, D. Jacquemin, I. Duchemin, and X. Blase, Lagrangian z-vector approach to bethe–salpeter analytic gradients: As- sessing approximations, J. Chem. Phys.159, 024116 (2023)

    2023

  15. [23]

    Tolle, M.-P

    J. Tolle, M.-P. Kitsaras, and P.-F. Loos, Fully analytic nuclear gradients for the bethe–salpeter equation, J. Phys. Chem. Lett.16, 11134 (2025)

    2025

  16. [24]

    Z. Dai, C. Lian, J. Lafuente-Bartolome, and F. Giustino, Excitonic polarons and self-trapped excitons from first- principles exciton-phonon couplings, Phys. Rev. Lett. 132, 036902 (2024)

    2024

  17. [25]

    Z. Dai, C. Lian, J. Lafuente-Bartolome, and F. Giustino, Theory of excitonic polarons: From models to first- principles calculations, Phys. Rev. B109, 045202 (2024)

    2024

  18. [26]

    Z. Dai, J. Lafuente-Bartolome, and F. Giustino, Po- larons from first principles (2025), arXiv:2512.06176 [cond-mat.mtrl-sci]

    2025

  19. [27]

    Y. Bai, Y. Wang, and S. Meng, Ab initio self-trapped excitons, Phys. Rev. Lett.133, 046903 (2024)

    2024

  20. [28]

    Yang and C

    M. Yang and C. Draxl, Novel approach to structural relaxation of materials in optically excited states, arXiv preprint arXiv:2212.13645 (2022)

    2022 arXiv

  21. [29]

    Govoni and G

    M. Govoni and G. Galli, Large scale GW calculations, J. Chem. Theory Comput.11, 2680 (2015)

    2015

  22. [30]

    V. W.-z. Yu and M. Govoni, GPU acceleration of large- scale full-frequency GW calculations, J. Chem. Theory 7 Comput.18, 4690 (2022)

    2022

  23. [31]

    V. W. zhe Yu, S. Chen, Y. Jin, V. Somjit, S. P. Villani, J. Zhan, M. Govoni, and G. Galli, The west code for large-scale excited-state materials simulations (2026), arXiv:2607.14025 [cond-mat.mtrl-sci]

    2026 arXiv

  24. [32]

    Rocca, D

    D. Rocca, D. Lu, and G. Galli, Ab initio calculations of optical absorption spectra: Solution of the bethe– salpeter equation within density matrix perturbation theory, J. Chem. Phys.133, 164109 (2010)

    2010

  25. [33]

    Rocca, Y

    D. Rocca, Y. Ping, R. Gebauer, and G. Galli, Solution of the bethe-salpeter equation without empty electronic states: Application to the absorption spectra of bulk systems, Phys. Rev. B85, 045116 (2012)

    2012

  26. [34]

    N. L. Nguyen, H. Ma, M. Govoni, F. Gygi, and G. Galli, Finite-field approach to solving the bethe-salpeter equa- tion, Phys. Rev. Lett.122, 237402 (2019)

    2019

  27. [35]

    V. W.-z. Yu, Y. Jin, G. Galli, and M. Govoni, Gpu- accelerated solution of the bethe–salpeter equation for large and heterogeneous systems, J. Chem. Theory Comput.20, 10899 (2024)

    2024

  28. [36]

    H. F. Wilson, F. m. c. Gygi, and G. Galli, Efficient iterative method for calculations of dielectric matrices, Phys. Rev. B78, 113303 (2008)

    2008

  29. [37]

    H. F. Wilson, D. Lu, F. m. c. Gygi, and G. Galli, Iter- ative calculations of dielectric eigenvalue spectra, Phys. Rev. B79, 245106 (2009)

    2009

  30. [38]

    Nguyen, T

    H.-V. Nguyen, T. A. Pham, D. Rocca, and G. Galli, Im- proving accuracy and efficiency of calculations of pho- toemission spectra within the many-body perturbation theory, Phys. Rev. B85, 081101(R) (2012)

    2012

  31. [39]

    T. A. Pham, H.-V. Nguyen, D. Rocca, and G. Galli,gw calculations using the spectral decomposition of the di- electric matrix: Verification, validation, and comparison of methods, Phys. Rev. B87, 155148 (2013)

    2013

  32. [40]

    Rohlfing and S

    M. Rohlfing and S. G. Louie, Electron-hole excitations and optical spectra from first principles, Phys. Rev. B 62, 4927 (2000)

    2000

  33. [41]

    Hutter, Excited state nuclear forces from the tamm– dancoff approximation to time-dependent density func- tional theory within the plane wave basis set framework, J

    J. Hutter, Excited state nuclear forces from the tamm– dancoff approximation to time-dependent density func- tional theory within the plane wave basis set framework, J. Chem. Phys.118, 3928 (2003)

    2003

  34. [42]

    Zhang and G

    X. Zhang and G. Lu, Subspace formulation of time- dependent density functional theory for large-scale cal- culations, J. Chem. Phys.143, 064110 (2015)

    2015

  35. [43]

    N. C. Handy and H. F. Schaefer III, On the evaluation of analytic energy derivatives for correlated wave func- tions, J. Chem. Phys.81, 5031 (1984)

    1984

  36. [44]

    [86–103]

    See the Supplemental Material for additional details on the analytical BSE forces formalism, validation of the implementation, applications to small molecules, com- putational parameters for the point-defect systems, fur- ther analysis of the carbon dimer defect in 2D-hBN, and...

  37. [45]

    M. J. Alrahamneh, I. Tolbatov, and P. Umari, Excited- state forces with gw-bse through the hellmann–feynman theorem, Int. J. Mol. Sci.26, 2306 (2025)

    2025

  38. [46]

    Gygi, J.-L

    F. Gygi, J.-L. Fattebert, and E. Schwegler, Computa- tion of maximally localized wannier functions using a simultaneous diagonalization algorithm, Comput. Phys. Commun.155, 1 (2003)

    2003

  39. [47]

    Simoncini and D

    V. Simoncini and D. B. Szyld, Theory of inexact Krylov subspace methods and applications to scientific comput- ing, SIAM J. Sci. Comput.25, 454 (2003)

    2003

  40. [48]

    Van Den Eshof and G

    J. Van Den Eshof and G. L. G. Sleijpen, Inexact Krylov subspace methods for linear systems, SIAM J. Matrix Anal. Appl.26, 125 (2004)

    2004

  41. [49]

    Lin, Adaptively compressed exchange operator, J

    L. Lin, Adaptively compressed exchange operator, J. Chem. Theory Comput.12, 2242 (2016)

    2016

  42. [50]

    V. W.-z. Yu and M. Govoni, Many-body perturbation theory with hybrid density functional theory starting points accelerated by adaptively compressed exchange, J. Chem. Phys.163, 024111 (2025)

    2025

  43. [51]

    M. W. Doherty, N. B. Manson, P. Delaney, F. Jelezko, J. Wrachtrup, and L. C. Hollenberg, The nitrogen- vacancy colour centre in diamond, Phys. Rep.528, 1 (2013)

    2013

  44. [52]

    Gali,Ab initiotheory of the nitrogen-vacancy center in diamond, Nanophotonics8, 1907 (2019)

    ´A. Gali,Ab initiotheory of the nitrogen-vacancy center in diamond, Nanophotonics8, 1907 (2019)

    1907

  45. [53]

    Schirhagl, K

    R. Schirhagl, K. Chang, M. Loretz, and C. L. Degen, Nitrogen-vacancy centers in diamond: nanoscale sensors for physics and biology, Annu. Rev. Phys. Chem.65, 83 (2014)

    2014

  46. [54]

    J. F. Barry, J. M. Schloss, E. Bauch, M. J. Turner, C. A. Hart, L. M. Pham, and R. L. Walsworth, Sensitivity optimization for nv-diamond magnetometry, Rev. Mod. Phys.92, 015004 (2020)

    2020

  47. [55]

    Childress and R

    L. Childress and R. Hanson, Diamond nv centers for quantum computing and quantum networks, MRS Bull. 38, 134 (2013)

    2013

  48. [56]

    Weber, W

    J. Weber, W. Koehl, J. Varley, A. Janotti, B. Buckley, C. Van de Walle, and D. D. Awschalom, Quantum com- puting with defects, Proc. Natl. Acad. Sci. U.S.A.107, 8513 (2010)

    2010

  49. [57]

    Waldherr, Y

    G. Waldherr, Y. Wang, S. Zaiser, M. Jamali, T. Schulte- Herbr¨ uggen, H. Abe, T. Ohshima, J. Isoya, J. F. Du, P. Neumann, and J. Wrachtrup, Quantum error correc- tion in a solid-state hybrid spin register, Nature506, 204 (2014)

    2014

  50. [58]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

    1996

  51. [59]

    J. H. Skone, M. Govoni, and G. Galli, Self-consistent hybrid functional for condensed systems, Phys. Rev. B 89, 195112 (2014)

    2014

  52. [60]

    J. H. Skone, M. Govoni, and G. Galli, Nonempiri- cal range-separated hybrid functionals for solids and molecules, Phys. Rev. B93, 235106 (2016)

    2016

  53. [61]

    S. Chen, V. W.-z. Yu, Y. Jin, M. Govoni, and G. Galli, Advances in quantum defect embedding theory, J. Chem. Theory Comput.21, 7797 (2025)

    2025

  54. [62]

    Davies and M

    G. Davies and M. Hamer, Optical studies of the 1.945 ev vibronic band in diamond, Proc. Math. Phys. Eng. Sci. P ROY SOC A-MATH PHY348, 285 (1976)

    1976

  55. [63]

    Museur, E

    L. Museur, E. Feldbach, and A. Kanaev, Defect-related photoluminescence of hexagonal boron nitride, Phys. Rev. B78, 155204 (2008)

    2008

  56. [64]

    Bourrellier, S

    R. Bourrellier, S. Meuret, A. Tararan, O. St´ ephan, M. Kociak, L. H. Tizei, and A. Zobelli, Bright uv single photon emission at point defects in h-bn, Nano lett.16, 4317 (2016)

    2016

  57. [65]

    Mackoit-Sinkeviˇ cien˙ e, M

    M. Mackoit-Sinkeviˇ cien˙ e, M. Maciaszek, C. G. Van de Walle, and A. Alkauskas, Carbon dimer defect as a source of the 4.1 ev luminescence in hexagonal boron nitride, Appl. Phys. Lett.115, 212101 (2019)

    2019

  58. [66]

    T. J. Smart, F. Wu, M. Govoni, and Y. Ping, Fun- damental principles for calculating charged defect ion- 8 ization energies in ultrathin two-dimensional materials, Phys. Rev. Mater.2, 124002 (2018)

    2018

  59. [67]

    Winter, M

    M. Winter, M. H. E. Bousquet, D. Jacquemin, I. Duchemin, and X. Blase, Photoluminescent properties of the carbon-dimer defect in hexagonal boron-nitride: A many-body finite-size cluster approach, Phys. Rev. Mater.5, 095201 (2021)

    2021

  60. [68]

    J. Zhan, M. Govoni, and G. Galli, Nonempirical range- separated hybrid functional with spatially dependent screened exchange, J. Chem. Theory Comput.19, 5851 (2023)

    2023

  61. [69]

    J. Zhan, M. Govoni, and G. Galli, Dielectric-dependent range-separated hybrid functional calculations for metal oxides, Phys. Rev. Mater.9, 053808 (2025)

    2025

  62. [70]

    Zhan and G

    J. Zhan and G. Galli, Hybrid functionals with spatially resolved screening for excited-state calculations (2026), in preparation

    2026

  63. [71]

    Iv´ ady, G

    V. Iv´ ady, G. Barcza, G. Thiering, S. Li, H. Hamdi, J.-P. Chou, ¨O. Legeza, and A. Gali, Ab initio theory of the negatively charged boron vacancy qubit in hexagonal boron nitride, npj Comput. Mater.6, 41 (2020)

    2020

  64. [72]

    Y. Lee, Y. Hu, X. Lang, D. Kim, K. Li, Y. Ping, K.-M. C. Fu, and K. Cho, Spin-defect qubits in two- dimensional transition metal dichalcogenides operating at telecom wavelengths, Nat. Commun.13, 7501 (2022)

    2022

  65. [73]

    Bertoldo, S

    F. Bertoldo, S. Ali, S. Manti, and K. S. Thygesen, Quan- tum point defects in 2d materials-the qpod database, npj Comput. Mater.8, 56 (2022)

    2022

  66. [74]

    Y. Jin, M. Rusishvili, M. Govoni, and G. Galli, Self- trapped excitons in metal-halide perovskites investi- gated by time-dependent density functional theory, J. Phys. Chem. Lett.15, 3229 (2024)

    2024

  67. [75]

    Davidsson, M

    J. Davidsson, M. Onizhuk, C. Vorwerk, and G. Galli, Discovery of atomic clock-like spin defects in simple ox- ides from first principles, Nat. Comm.15, 4812 (2024)

    2024

  68. [76]

    Somjit, J

    V. Somjit, J. Davidsson, Y. Jin, and G. Galli, An nv- center in magnesium oxide as a spin qubit for hy- brid quantum technologies, npj Comput. Mater.11, 74 (2025)

    2025

  69. [77]

    Zhang, E

    S. Zhang, E. Perez, T. Park, K. Li, X. Wang, M. Man- souri, Y. Wang, J. D. V. Bazantes, R. Zhang, J. Sun, K.-M. C. Fu, H. Seo, and Y. Ping, Deep spin defects in zinc oxide for high-fidelity single-shot readout, PRX Quantum7, 020348 (2026)

    2026

  70. [78]

    S. P. Villani, Y. Jin, and G. Galli, First-principles cal- culations of internal conversion processes in spin defects (2026), arXiv:2606.14545 [cond-mat.mtrl-sci]

    2026

  71. [79]

    Malosso, W

    C. Malosso, W. B. How, G. D. Mir´ on, A. Hassanali, and M. Ceriotti, Transferable machine learning of excited-state dynamics with extremal pooling (2026), arXiv:2606.16859 [physics.chem-ph]

    2026

  72. [80]

    S. S. Dong, M. Govoni, and G. Galli, Machine learning dielectric screening for the simulation of excited state properties of molecules and materials, Chem. Sci.12, 4970 (2021)

    2021

  73. [81]

    M. G. Zauchner, A. Horsfield, and J. Lischner, Accel- erating gw calculations through machine-learned dielec- tric matrices, npj Comput. Mater.9, 184 (2023)

    2023

  74. [82]

    Giannozzi, O

    P. Giannozzi, O. Baseggio, P. Bonf` a, D. Brunato, R. Car, I. Carnimeo, C. Cavazzoni, S. de Gironcoli, P. Delugas, F. F. Ruffino, A. Ferretti, N. Marzari, I. Timrov, A. Urru, and S. Baroni, Quantum ESPRESSO toward the exascale, J. Chem. Phys.152, 154105 (2020)

    2020

  75. [83]

    Carnimeo, F

    I. Carnimeo, F. Affinito, S. Baroni, O. Baseggio, L. Bel- lentani, R. Bertossa, P. D. Delugas, F. F. Ruffino, S. Or- landini, F. Spiga, and P. Giannozzi, Quantum espresso: One further step toward the exascale, J. Chem. Theory Comput.19, 6992 (2023)

    2023

  76. [84]

    Y. Jin, V. Yu, M. Govoni, and G. Galli, PyPL: Python code to compute vibrationally resolved optical spectra (2026), a Python package for computing vibrationally resolved optical spectra of point defects. [85]https://github.com/jinyuchem/supporting_data/ tree/main/2026/arXiv

    2026

  77. [86]

    Budzak, G

    S. Budzak, G. Scalmani, and D. Jacquemin, Accurate excited-state geometries: A caspt2 and coupled-cluster reference database for small molecules, J. Chem. Theory Comput.13, 6237 (2017)

    2017

  78. [87]

    Adamo and V

    C. Adamo and V. Barone, Toward reliable density func- tional methods without adjustable parameters: The pbe0 model, J. Chem. Phys.110, 6158 (1999)

    1999

  79. [88]

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid func- tionals based on a screened coulomb potential, J. Chem. Phys.118, 8207 (2003)

    2003

  80. [89]

    A. V. Krukau, O. A. Vydrov, A. F. Izmaylov, and G. E. Scuseria, Influence of the exchange screening parameter on the performance of screened hybrid functionals, J. Chem. Phys.125, 224106 (2006)

    2006

  81. [90]

    D. R. Hamann, Optimized norm-conserving Vanderbilt pseudopotentials, Phys. Rev. B88, 085117 (2013)

    2013

  82. [91]

    Schlipf and F

    M. Schlipf and F. Gygi, Optimization algorithm for the generation of ONCV pseudopotentials, Comput. Phys. Commun.196, 36 (2015)

    2015

  83. [92]

    A. Togo, L. Chaput, T. Tadano, and I. Tanaka, Imple- mentation strategies in phonopy and phono3py, J. Phys. Condens. Matter35, 353001 (2023)

    2023

  84. [93]

    Togo, First-principles phonon calculations with phonopy and phono3py, J

    A. Togo, First-principles phonon calculations with phonopy and phono3py, J. Phys. Soc. Jpn.92, 012001 (2023)

    2023

  85. [94]

    Razinkovas, M

    L. Razinkovas, M. W. Doherty, N. B. Manson, C. G. Van de Walle, and A. Alkauskas, Vibrational and vi- bronic structure of isolated point defects: The nitrogen- vacancy center in diamond, Phys. Rev. B104, 045303 (2021)

    2021

  86. [95]

    Y. Jin, M. Govoni, G. Wolfowicz, S. E. Sullivan, F. J. Heremans, D. D. Awschalom, and G. Galli, Photolu- minescence spectra of point defects in semiconductors: Validation of first-principles calculations, Phys. Rev. Mater.5, 084603 (2021)

    2021

  87. [96]

    Muechler, D

    L. Muechler, D. I. Badrtdinov, A. Hampel, J. Cano, M. R¨ osner, and C. E. Dreyer, Quantum embedding methods for correlated excited states of point defects: Case studies and challenges, Phys. Rev. B105, 235104 (2022)

    2022

  88. [97]

    Haldar, A

    S. Haldar, A. Mitra, M. R. Hermes, and L. Gagliardi, Local excitations of a charged nitrogen vacancy in di- amond with multireference density matrix embedding theory, J. Phys. Chem. Lett.14, 4273 (2023)

    2023

  89. [98]

    Benedek, ´A

    Z. Benedek, ´A. Ganyecz, A. Pershin, V. Iv´ ady, and G. Barcza, Accurate and convergent energetics of color centers by wavefunction theory, npj Comput. Mater.11, 346 (2025)

    2025

  90. [99]

    J. M. P. Martirez, Optical properties of a diamond nv color center from capped embedded multiconfigu- rational correlated wavefunction theory, npj Comput. 9 Mater.12, 113 (2026)

    2026

  91. [100]

    B. T. Lau, B. Busemeyer, and T. C. Berkelbach, Optical properties of defects in solids via quantum embedding with good active space orbitals, J. Phys. Chem. C128, 2959 (2024)

    2024

  92. [101]

    Korona and M

    T. Korona and M. Chojecki, Exploring point defects in hexagonal boron-nitrogen monolayers, Int. J. Quantum Chem.119, e25925 (2019)

    2019

  93. [102]

    J. Li, Y. Jin, J. Yu, W. Yang, and T. Zhu, Accu- rate excitation energies of point defects from fast parti- cle–particle random phase approximation calculations, J. Phys. Chem. Lett.15, 2757 (2024)

    2024

  94. [103]

    Analytical Forces from the Bethe-Salpeter Equation for Large-Scale Excited-State Relaxation

    J. Li, Y. Jin, J. Yu, W. Yang, and T. Zhu, Particle– particle random phase approximation for predicting cor- related excited states of point defects, J. Chem. Theory Comput.20, 7979 (2024). END MA TTER 6 7 8 (eV ) R0 = 1.26 /Aring.low (a) BSE + scissor B S E + G0W0 Q P 1.10 1....

    2024

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.