REVIEW 3 major objections 6 minor 58 references
Real-Time Approach to the Dynamical Bethe-Salpeter Equation for Finite Systems
T0 review · 3 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper claims that the frequency-dependent kernel of the dynamical Bethe-Salpeter equation can be built by real-time propagation and a time-domain product, avoiding frequency-grid integration and storage of the screened Coulomb…
desk verdict A real algorithmic step for dynamical BSE, but the reported shifts test a Shindo-approximated kernel, not the full-frequency BSE. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is Eq. (17), the real-time expression for the polarization contribution to the dynamical kernel: it converts the convolution over frequency in Eq. (3) into multiplication in time by the phase factors $e^{-i(\varepsilon_b-\varepsilon_i+\Delta)t}$ and $e^{i(\varepsilon_a-\varepsilon_j+\Delta)t}$, followed by a Fourier transform. The retarded response $W^{\mathrm{pol},ij}_{ab}(t)$ is generated by time-dependent Hartree propagation: an impulsive perturbation from each occupied pair $kl$ produces an induced density matrix, which is contracted with the Coulomb tensor to give the polarization part of the screened interaction. A causal-to-time-ordered transformation puts this response on the full time axis. The kernel itself rests on the Shindo-type factorization of the four-point correlation function, Eq. (A15), which allows the two internal fermionic frequency integrals to be evaluated by contour integration so that the BSE closes.
What would settle it
Compute the dynamical excitation energies of a small molecule whose low-lying polarization mode lies near the optical excitation energy, using an unfactorized frequency-domain solution of the Bethe-Salpeter equation that retains full frequency dependence in W, and compare with the real-time time-product result. Agreement within the roughly 50 meV numerical broadening would support the Shindo factorization; a larger mismatch would show that Eq. (17) is not the full-frequency BSE kernel.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that Eq. (17) is an equivalent time-domain rendering of the frequency convolution defining the dynamical BSE kernel. The retarded screened interaction is obtained by propagating the time-dependent Hartree density matrix after impulsive perturbations built from occupied-orbital pair densities; the causal response is converted to a time-ordered interaction, multiplied by phase factors carrying the bare electron-hole transition energies, and Fourier transformed to give the full-frequency-dependent kernel without explicit frequency integration or storage of the full real-space screened Coulomb operator. The construction is validated on silane in two ways: static BSE from real-time response matches a standard dielectric-matrix-inversion implementation with mean absolute deviations below 1 meV, and the time-product evaluation of the dynamical interaction reproduces the direct frequency convolution. Dynamical screening shifts the lowest singlet excitation down by 0.12 eV and the triplet by 0.35 eV, a redshift in the same direction as earlier dynamical-BSE results.
Load-bearing premise
The load-bearing premise is the Shindo-type factorization of the four-point correlation function in Eq. (A15): if that factorization is not quantitatively accurate for molecular screening, the computed dynamical shifts are not the true full-frequency BSE shifts, and the paper gives no benchmark against a dynamical BSE solver that avoids it.
Editorial extensions
If this is right
- For systems where polarization modes lie near optical excitation energies, the static approximation can misestimate singlet and triplet energies by roughly 0.1-0.3 eV; the real-time dynamical kernel corrects this at the cost of propagating TDH equations rather than storing a frequency grid.
- The nonlinear BSE can be solved by diagonalizing the dynamical matrix on a sparse grid of trial frequencies and locating roots by secant interpolation; for silane the matrix varies smoothly, so few sampling frequencies suffice.
- The same real-time response supplies both the static and the dynamical kernel, so a single code path covers both approximations and can be checked against independent static implementations.
- The time-domain form is compatible with stochastic sampling, in which random occupied and unoccupied orbital combinations carry the orbital-energy phase factors, offering a route to dynamical BSE for systems with hundreds to thousands of valence electrons.
- The time-ordered interaction can also drive a memory-based density-matrix evolution, outlined in Appendix B, potentially yielding spectra or steady states without a frequency-domain diagonalization.
Reading between the lines
- A natural test of the Shindo factorization would be to compare the real-time kernel with an unfactorized dynamical BSE calculation on a small molecule whose low-lying polarization mode overlaps the optical excitation; the paper does not provide such a benchmark.
- The real-time construction may be most beneficial where memory rather than floating-point operations is the bottleneck, since no screened-Coulomb operator is stored on a frequency grid; profiling memory against system size would identify the crossover.
- If the stochastic extension succeeds, the dynamical BSE could become available precisely for the molecules where the static approximation is most questionable, such as weakly screened systems and triplet excitations, because the phase factors factorize onto occupied and unoccupied subspaces.
- The Appendix B memory-based evolution suggests that the same time-dependent interaction could bypass the nonlinear root search entirely by propagating a density matrix; whether that evolution is numerically stable is left open.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a real-time linear-response formulation of the dynamical Bethe-Salpeter equation for finite systems. The screened interaction entering the BSE kernel is constructed by propagating time-dependent Hartree equations under impulsive pair-density perturbations, and the frequency convolution that defines the dynamical kernel is evaluated as a product in the time domain followed by Fourier transformation. This avoids explicit frequency grids and storage of the full screened Coulomb operator. The method is validated on SiH4 in two ways: the static BSE from real-time TDH response agrees with a frequency-domain dielectric-matrix inversion to below 1 meV, and the time-product evaluation of the dynamical interaction reproduces direct numerical frequency convolution. The dynamical BSE redshifts the lowest singlet by 0.12 eV and triplet by 0.35 eV, in qualitative agreement with earlier work by Rohlfing and Louie. The paper also sketches a stochastic extension and a memory-based time-evolution formulation.
Significance. If the central claims hold, the method gives a practical route to dynamical BSE calculations that avoids explicit frequency integration and storage of the screened interaction, and its compatibility with stochastic sampling is a genuine algorithmic advantage. The algebraic transformation from the frequency convolution to a time-domain product is exact and internally consistent, and the reproduction of an independent static implementation is a useful validation. However, the physical significance of the reported dynamical shifts is conditional on the Shindo-type factorization adopted in Appendix A, which is not benchmarked against a reference that avoids that approximation. The validation is limited to a single small molecule with no convergence tests, so the quantitative redshifts should be treated as preliminary.
major comments (3)
- [II.A and Appendix A, Eq. (A15)] The central physical claim—that the method goes beyond the static screening approximation with full-frequency dependence in the screened interaction—rests on the Shindo-type factorization in Eq. (A15), but this factorization is not validated. The agreement between the time-product and frequency-convolution routes reported in Section III only confirms that Eq. (17) correctly implements Eq. (3); it does not test the approximation in Eq. (3) itself. Since the dynamical shifts in Table I are presented as physical predictions, a comparison against a dynamical BSE solver that avoids the Shindo factorization, using the same LDA+scissor input, is needed to support the attribution of these shifts to full-frequency screening.
- [III, Table I] The comparison with Ref. [13] cannot serve as a quantitative validation because the reference uses GW quasiparticle states and a plasmon-pole model, as the authors acknowledge in the text. The statement that dynamical screening redshifts the excitations is supported, but the magnitudes (0.12 eV singlet, 0.35 eV triplet) are not benchmarked. Please either provide a same-input reference calculation or explicitly label these magnitudes as model-dependent within the Shindo approximation.
- [III, Eqs. (13)-(19)] The reported numerical results have no convergence analysis. The choices N_c=24, dt=0.05 au, T_max=20 fs, gamma~100 meV, and eta~50 meV are stated, but their effect on the dynamical shifts is not quantified. Because the shifts are of order 0.1-0.35 eV, the numerical error from Gaussian damping, finite propagation time, and the frequency grid used for secant interpolation could be a substantial fraction of the effect. Please report convergence with respect to T_max (or gamma), eta, basis size, and number of frequency samples.
minor comments (6)
- [Eqs. (14) and (17)] The retarded response in Eq. (14) uses Gaussian damping e^{-gamma^2 t^2/2}, while the time-domain convolution in Eq. (17) introduces e^{-eta|t|}. The relationship between gamma and eta, and how both appear in the final time-ordered W(t) entering Eq. (17), should be stated explicitly.
- [Eq. (20)] Eq. (20) writes a relation for the full W(t), whereas Eq. (17) uses the polarization part W_pol(t), with the bare Coulomb term added separately in Eq. (18). Please clarify whether Eq. (20) includes the instantaneous delta-function term and how the delta function is handled in the stochastic sampling context.
- [Section III] The secant root-finding procedure is described in words but the frequency grid spacing and the number of sampled frequencies are not given. This information is needed to assess the numerical stability of the reported roots.
- [References] References [52] and [53] are identical (Romaniello, Guyot, and Reining, J. Chem. Phys. 131, 154111 (2009)); the duplicate should be removed or replaced.
- [Table I caption] The column layout of Table I is confusing: the columns labeled 'Static', 'Time Product', and 'Frequency Convolution' under the present work, followed by 'Static' and 'Dynamic' under Ref. [13], are not clearly separated. The caption should explicitly indicate which columns refer to which implementation.
- [Eq. (21)] The stochastic notation, particularly the distinction between overlined and double-overlined random states, is introduced abruptly. A brief definition of the sampling distribution and the normalization of the random states would improve reproducibility.
Circularity Check
No significant circularity: Eq. (17) is an exact time-domain reformulation of Eq. (3), and no fitted parameter is renamed as a prediction.
full rationale
The derivation chain is self-contained. Eq. (17) is an exact Fourier/convolution identity relative to Eq. (3): the paper explicitly states that the orbital-energy phase factors originate from the energy denominators in Eq. (3), so the time-product evaluation is mathematically equivalent to the frequency convolution. This is an algebraic reformulation, not an output-dependent fit or a self-referential derivation. The paper also explicitly labels Eq. (3) as a Shindo-like approximation and derives it in Appendix A through the stated factorization (A15), rather than importing the kernel as a fitted or externally imposed result. The screened interaction W(omega) is obtained from time-dependent Hartree/RPA response, an independent model, and no parameter is fitted to the reported excitation energies. The only external fit, the scissor shift Delta, is fixed to the GW gap of Ref. [13], not to the BSE excitation energies that are the outputs. The validation that the time product reproduces the frequency convolution is tautological in the sense that both evaluate the same kernel, but it is a numerical consistency check, not a circular claim. Self-citations appear only in describing the causal-to-time-ordered transformation and in speculative future stochastic extensions; they are not load-bearing for the central derivation. The absence of a non-Shindo dynamical BSE benchmark is a validation gap and a correctness risk, but it is not circularity under the criteria used here.
Assumptions & free parameters
free parameters (3)
- Scissor correction Delta =
7.18 eV
- Gaussian damping gamma =
3/T_max ~ 100 meV
- Numerical broadening eta =
~50 meV
assumptions (6)
- domain assumption The screened interaction W(omega) is obtained at the time-dependent Hartree/RPA level via Eq. (11).
- domain assumption The Tamm-Dancoff approximation (TDA), neglecting the resonant-antiresonant coupling block, is used for the BSE.
- domain assumption The Shindo-like factorization, Eq. (A15), is accurate enough for the dynamical BSE.
- domain assumption LDA-DFT orbitals with a rigid scissor sufficiently approximate GW quasiparticle states for the BSE.
- domain assumption The finite-grid plane-wave LDA ground state is converged sufficiently for BSE excitations.
- standard math Fourier and convolution identities, plus contour integration, used in Appendix A are valid.
Cite this review
Pith. "Pith review of Real-Time Approach to the Dynamical Bethe-Salpeter Equation for Finite Systems." pith.science (2026). https://pith.science/paper/R72QPSWY
@misc{pith2026260808948,
author = {Pith},
title = {Pith review of: Real-Time Approach to the Dynamical Bethe-Salpeter Equation for Finite Systems},
year = {2026},
howpublished = {\url{https://pith.science/paper/R72QPSWY}},
note = {Machine review of arXiv:2608.08948}
}
abstract
We present a real-time linear-response approach to solving the dynamical Bethe-Salpeter equation (BSE). The polarization part of the screened interaction is obtained from time-dependent Hartree propagation of the one-particle density matrix within an orbital basis-set representation. The frequency convolution defining the dynamical kernel is evaluated as a product in time, avoiding numerical integration and storage of the full screened Coulomb operator. The resulting nonlinear eigenvalue problem is then solved directly, going beyond the static screening approximation with full-frequency dependence in the screened interaction. While the deterministic scaling remains $\mathcal{O}(N^6)$, the time propagation formulation is readily compatible with grid-based stochastic sampling methods, which will open the possibility for dynamical BSE calculations of very large systems.
Figures
Reference graph
Works this paper leans on
-
[13]
Rohlfing and S
M. Rohlfing and S. G. Louie, Phys. Rev. B62, 4927 (2000)
2000
-
[1]
C. A. Ullrich, APL Computational Physics1, 020901 (2025)
2025
-
[2]
N. T. Maitra, F. Zhang, R. J. Cave, and K. Burke, The Journal of Chemical Physics120, 5932 (2004)
work page 2004
-
[3]
B. G. Levine, C. Ko, J. Quenneville, and T. J. Mart ´Inez, Molecular Physics104, 1039 (2006)
work page 2006
-
[4]
P. Romaniello, D. Sangalli, J. A. Berger, F. Sottile, L. G. Molinari, L. Reining, and G. Onida, The Journal of Chemical Physics130(2009)
work page 2009
-
[5]
D. Sangalli, P. Romaniello, G. Onida, and A. Marini, The Journal of Chemical Physics134 (2011)
work page 2011
- [6]
-
[7]
E. Rebolini and J. Toulouse, The Journal of Chemical Physics144, 094107 (2016)
work page 2016
Show all 58 references
-
[8]
Lacombe and N
L. Lacombe and N. T. Maitra, npj Computational Materials9(2023)
2023
-
[9]
Hedin, Phys
L. Hedin, Phys. Rev.139, A796 (1965)
1965
-
[10]
M. S. Hybertsen and S. G. Louie, Physical Review Letters55, 1418 (1985)
1985
-
[11]
Strinati, La Rivista del Nuovo Cimento11, 1 (1988)
G. Strinati, La Rivista del Nuovo Cimento11, 1 (1988)
1988
-
[12]
Blase, I
X. Blase, I. Duchemin, D. Jacquemin, and P.-F. Loos, The Journal of Physical Chemistry Letters11, 7371 (2020)
2020
-
[14]
Onida, L
G. Onida, L. Reining, and A. Rubio, Rev. Mod. Phys.74, 601 (2002)
2002
-
[15]
Jacquemin, I
D. Jacquemin, I. Duchemin, and X. Blase, The Journal of Physical Chemistry Letters8, 1524 (2017)
2017
-
[16]
Li and V
J. Li and V. Olevano, Journal of Photochemistry and Photobiology B: Biology232, 112475 (2022)
2022
-
[17]
Rocca, Y
D. Rocca, Y. Ping, R. Gebauer, and G. Galli, Physical Review B85, 045116 (2012)
2012
-
[18]
Jacquemin, I
D. Jacquemin, I. Duchemin, and X. Blase, Journal of Chemical Theory and Computation11, 3290 (2015)
2015
-
[19]
X. Gui, C. Holzer, and W. Klopper, Journal of Chemical Theory and Computation14, 2127 19 (2018)
2018
-
[20]
F¨ orster and L
A. F¨ orster and L. Visscher, Journal of Chemical Theory and Computation18, 6779 (2022)
2022
-
[21]
Marini and R
A. Marini and R. Del Sole, Physical Review Letters91, 176402 (2003)
2003
-
[22]
Bechstedt,Many-Body Approach to Electronic Excitations: Concepts and Applications (Springer Berlin Heidelberg, 2015)
F. Bechstedt,Many-Body Approach to Electronic Excitations: Concepts and Applications (Springer Berlin Heidelberg, 2015)
2015
-
[23]
Y. Ma, M. Rohlfing, and C. Molteni, Physical Review B80, 241405 (2009)
2009
-
[24]
Loos and X
P.-F. Loos and X. Blase, The Journal of Chemical Physics153, 114120 (2020)
2020
-
[25]
Zhang, J
X. Zhang, J. A. Leveillee, and A. Schleife, Physical Review B107, 235205 (2023)
2023
-
[26]
M. Wen, G. Harsha, and D. Zgid, Journal of Chemical Theory and Computation (2026)
2026
-
[27]
All-electron dynamical Bethe-Salpeter Equa- tion for extended systems with atom-centered orbital basis set,
R. Zhou, S. Liu, J. Xu, Y. Yao, and Y. Kanai, “All-electron dynamical Bethe-Salpeter Equa- tion for extended systems with atom-centered orbital basis set,” (2026), arXiv:2606.08350 [physics.chem-ph]
2026 arXiv
-
[28]
S. J. Bintrim and T. C. Berkelbach, The Journal of Chemical Physics156(2022)
2022
-
[29]
Kurzweil and R
Y. Kurzweil and R. Baer, Physical Review B77(2008)
2008
-
[30]
K. J. Daas, S. Crisostomo, and K. Burke, Phys. Rev. Lett.137, 028002 (2026)
2026
-
[31]
Casanova and A
D. Casanova and A. I. Krylov, Physical Chemistry Chemical Physics22, 4326 (2020)
2020
-
[32]
W. Park, J. Shen, S. Lee, P. Piecuch, M. Filatov, and C. H. Choi, The Journal of Physical Chemistry Letters12, 9720 (2021)
2021
-
[33]
D. B. Dar and N. T. Maitra, The Journal of Physical Chemistry Letters16, 703 (2025)
2025
-
[34]
Authier and P.-F
J. Authier and P.-F. Loos, The Journal of Chemical Physics153, 184105 (2020)
2020
-
[35]
Loos and P
P.-F. Loos and P. Romaniello, The Journal of Chemical Physics156, 164101 (2022)
2022
-
[36]
H. F. Wilson, F. m. c. Gygi, and G. Galli, Phys. Rev. B78, 113303 (2008)
2008
-
[37]
X. Ren, P. Rinke, V. Blum, J. Wieferink, A. Tkatchenko, A. Sanfilippo, K. Reuter, and M. Scheffler, New Journal of Physics14, 053020 (2012)
2012
-
[38]
Neuhauser, Y
D. Neuhauser, Y. Gao, C. Arntsen, C. Karshenas, E. Rabani, and R. Baer, Phys. Rev. Lett. 113, 076402 (2014)
2014
-
[39]
Duchemin and X
I. Duchemin and X. Blase, Journal of Chemical Theory and Computation17, 2383 (2021)
2021
-
[40]
C.-N. Yeh, S. Iskakov, D. Zgid, and E. Gull, Physical Review B106, 235104 (2022)
2022
-
[41]
Rabani, R
E. Rabani, R. Baer, and D. Neuhauser, Physical Review B91, 235302 (2015)
2015
-
[42]
N. C. Bradbury, T. Allen, M. Nguyen, K. Z. Ibrahim, and D. Neuhauser, The Journal of 20 Chemical Physics158, 154104 (2023)
2023
-
[43]
Hillenbrand, J
C. Hillenbrand, J. Li, and T. Zhu, The Journal of Chemical Physics162, 174117 (2025)
2025
-
[44]
Vlˇ cek, W
V. Vlˇ cek, W. Li, R. Baer, E. Rabani, and D. Neuhauser, Physical Review B98, 075107 (2018)
2018
-
[45]
N. C. Bradbury, M. Nguyen, J. R. Caram, and D. Neuhauser, The Journal of Chemical Physics 157, 031104 (2022)
2022
-
[46]
N. C. Bradbury, B. Y. Li, T. Allen, J. R. Caram, and D. Neuhauser, The Journal of Chemical Physics161, 141101 (2024)
2024
-
[47]
Allen, B
T. Allen, B. Y. Li, N. C. Bradbury, and D. Neuhauser, Journal of Chemical Theory and Computation , acs.jctc.6c00440 (2026)
2026
-
[48]
J. P. Perdew and Y. Wang, Phys. Rev. B45, 13244 (1992)
1992
-
[49]
A. L. Fetter and J. D. Walecka,Quantum Theory of Many-Particle Systems, Dover Books on Physics (Dover Publications, Newburyport, 2012)
2012
-
[50]
Troullier and J
N. Troullier and J. L. Martins, Phys. Rev. B43, 1993 (1991)
1991
-
[51]
G. J. Martyna and M. E. Tuckerman, The Journal of Chemical Physics110, 2810 (1999)
1999
-
[53]
Romaniello, S
P. Romaniello, S. Guyot, and L. Reining, The Journal of Chemical Physics131, 154111 (2009)
2009
-
[54]
Neuhauser, E
D. Neuhauser, E. Rabani, and R. Baer, Journal of Chemical Theory and Computation9, 24 (2012)
2012
-
[55]
R. Baer, D. Neuhauser, and E. Rabani, Phys. Rev. Lett.111, 106402 (2013)
2013
-
[56]
N. C. Bradbury, T. Allen, M. Nguyen, and D. Neuhauser, Journal of Chemical Theory and Computation19, 9239 (2023)
2023
-
[57]
Y. Gao, D. Neuhauser, R. Baer, and E. Rabani, The Journal of Chemical Physics142, 034106 (2015)
2015
-
[58]
Allen, M
T. Allen, M. Nguyen, and D. Neuhauser, The Journal of Chemical Physics161, 114116 (2024)
2024
-
[59]
P. S. Thomas, M. Nguyen, D. Bazile, T. Allen, B. Y. Li, W. Li, M. Del Ben, J. Deslippe, and D. Neuhauser, Journal of Chemical Theory and Computation22, 3960 (2026). 21
2026
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.