Pith. sign in

REVIEW 3 major objections 3 minor 27 references

Spatial-Order Hierarchy of Time-Dependent Exchange-Correlation Potential

T0 review · 3 major / 3 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read In exact TDDFT, a linear off-diagonal deviation of the density matrix from a Hartree-Fock reference already produces a non-HF current, while kinetic correlation activates only at second order in the relative coordinate.

desk verdict The hierarchy claim is a good idea, but the proof swaps the KS 1RDM for the TDHF one, so the main theorem doesn't go through. read the letter →

arxiv 2608.01483 v1 pith:PQLMFK2Q submitted 2026-08-02 quant-ph physics.comp-ph

Naoki Negishi This is my paper · ORCID
classification quant-phphysics.comp-ph
keywords time-dependentdensityfunctionaltheoryexchange-correlationpotentialkinetic-correlationreducedmatrixoff-diagonalexpansionTDHFreferencenon-adiabaticspatial-orderhierarchy
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 claims that the two pieces of the time-dependent exchange-correlation potential in TDDFT—the electron-interaction term $v^W_{xc}$ and the kinetic-correlation term $v^T_c$—are locked together by the off-diagonal structure of the one-electron reduced density matrix. Expanding the difference between the exact 1RDM and a time-dependent Hartree-Fock reference in the relative coordinate $r-r'$, the author shows that a non-zero first-order coefficient already produces a current beyond TDHF while $v^T_c$ retains its Hartree-Fock form; only quadratic and higher off-diagonal terms force kinetic correlation to change. If correct, this gives an exact, spatially local constraint for building non-adiabatic functionals and explains when the adiabatic picture is acceptable. A numerical two-electron quench backs the claim by showing that dropping $v^T_c$ freezes the density broadening seen in the exact evolution.

What carries the argument

The central object is the relative-coordinate expansion of the one-body reduced density matrix, $\gamma(r,r',t) = \gamma_{\rm HF}(r,r',t) + n[\Delta\hat\gamma] + i(r-r')\cdot j[\Delta\hat\gamma] + \cdots$, with the kinetic-correlation potential $v^T_c$ defined through the second-order $(\nabla-\nabla')^2$ coefficient of $\gamma-\gamma_s$. The off-diagonal distance $|r-r'|$ is the bookkeeping parameter: linear order controls the current, quadratic order controls kinetic correlation, so the hierarchy is read off from a single expansion.

What would settle it

Take a two-electron model, propagate the exact TDSE and TDHF from a common HF starting point, and expand $\Delta\gamma(x,x',t)$ in $x-x'$. If there exists a time interval where the first-order coefficient is non-zero and the second-order coefficient vanishes while $n\nabla v^T_c$ computed from Eq. (6) differs from its HF reference value, Theorem 1's characterization fails. Alternatively, construct a KS system with the exact density but $\gamma_s \neq \gamma_{\rm HF}$; if Eq. (10) then fails, the identification premise is the weak point.

Watch

Extended reading notes

Core claim

On the paper's own terms: Theorem 1 sets a sufficient condition on the exact 1RDM—non-idempotent, with a non-trivial first-order and trivial second-order difference from the TDHF 1RDM in $r-r'$—under which the exact current deviates from the TDHF current, the kinetic-correlation potential stays representable in the HF sense (Eq. 10), and the interaction exchange-correlation force cannot equal the TDHF exchange form (Eq. 11). The consequence is a spatial-order hierarchy: interaction-driven correlation enters first, kinetic correlation at second order, giving a representability constraint on exact functional decompositions.

Load-bearing premise

The proof relies on identifying the Kohn-Sham one-body density matrix with the TDHF one; if that identification is not valid, the claim that kinetic correlation stays 'HF-representable' does not follow from the stated premises.

Editorial extensions

If this is right

  • If the hierarchy is exact, functionals that capture interaction correlation to infinite order but set $v^T_c=0$ are incomplete for non-stationary dynamics; kinetic correlation is required whenever the 1RDM develops quadratic off-diagonal structure.
  • The result gives a classification of density equations of motion into three classes (mean-field, high-accuracy $v^W_{xc}$ with trivial $v^T_c$, and exact), making the approximation trade-off explicit.
  • Adiabatic approximations are on firmer ground for small off-diagonal deviations, since first-order truncation keeps kinetic correlation trivial.
  • The theorem supplies a concrete necessary condition: any non-adiabatic functional whose $v^W_{xc}$ and $v^T_c$ violate the spatial-order relation cannot reproduce exact density dynamics.

Reading between the lines

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

  • A natural extension is to use the order of the first non-vanishing off-diagonal coefficient as a diagnostic of non-adiabaticity in real-time simulations, something the paper does not explicitly propose.
  • The same expansion could be applied in current-density functional theory, where the current (first-order term) is a basic variable; the hierarchy would then constrain the kinetic-correlation functional directly.
  • One could test the theorem's boundary by weakening 'trivial second-order coefficient' to 'small compared to first order' and checking numerically whether $v^T_c$ remains approximately HF-representable, giving a practical error bound.
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

3 major / 3 minor

Summary. The paper proposes a spatial-order hierarchy relating the two components of the time-dependent exchange-correlation potential in TDDFT: the electron-interaction component v_W_xc and the kinetic-correlation component v_T_c. Theorem 1 claims that if the exact 1RDM differs from the TDHF 1RDM by a non-zero linear term in the relative coordinate but a zero quadratic term, then (i) the exact current differs from the TDHF current (Eq. 9), (ii) the kinetic-correlation component remains 'HF representable' (Eq. 10), and (iii) the interaction xc potential cannot take the TDHF exchange form (Eq. 11). The proof is given in an Appendix using an expansion of the 1RDM in (r-r') and (t-t0). A numerical experiment on a one-dimensional two-electron soft-Coulomb quench compares exact TDSE propagation with a truncated propagation that sets v_T_c=0, showing that the truncation suppresses the density broadening seen in the exact evolution.

Significance. If correct, the claimed constraint would be a useful exact structural relation for constructing non-adiabatic TDDFT functionals: it would specify the off-diagonal spatial order at which kinetic correlation becomes unavoidable. The paper is clearly written, and the numerical illustration is concrete and reproducible in principle. However, the central theorem is not established: Eq. (10) is derived by conflating the Kohn-Sham 1RDM with the TDHF 1RDM, and Eq. (11) rests on an unproved analyticity/induction argument. The numerical result is a single two-electron trajectory and cannot compensate for the proof gap. The conceptual decomposition and the question addressed are important, but the main claim is unsupported as stated.

major comments (3)
  1. The proof conflates the KS 1RDM γ_s with the TDHF 1RDM γ_HF. From Eq. (6), n∇v_T_c[γ] = -(1/4)L(γ - γ_s[γ]) and similarly for γ_HF, with L=(∇-∇')(∇²-∇'²)|_{r=r'}. The assumption that Δγ has a trivial quadratic coefficient gives L(γ-γ_HF)=0. Eq. (10) then reduces to L(γ_s[γ])=L(γ_s[γ_HF]). No stated premise implies this. The Appendix's proof of Eq. (10) never mentions γ_s and asserts that Eq. (6) 'equivalently leads to a trivial solution', which is true only after setting γ_s=γ_HF. In TDDFT, γ_s is determined by the density and is generally not the TDHF 1RDM; for the two-electron singlet it is rank-one, but with a different orbital than γ_HF. Thus Eq. (10) does not follow from the hypotheses, and the central 'kinetic component remains HF representable' claim fails.
  2. The proof of Eq. (11) depends on an expansion in (t-t0) and on an induction that ∂_t^m(n[γ]-n[γ_HF])=0 and ∂_t^m(j[γ]-j[γ_HF])=0 for all m below the first non-zero order, together with an implicit time-analyticity assumption. Neither the analyticity of n[Δγ] and j[Δγ] near t0 nor the induction is proved; the text only asserts it ('This is because it can be proved inductively'). For general many-body dynamics with soft-Coulomb interactions, analyticity in t is not automatic. Without these steps, the existence of a time region where ∇v_W_xc departs from the TDHF exchange form is not established.
  3. The theorem is stated as 'necessary and sufficient', but the appendix provides only a conditional/sufficient argument for (9) and (11) and gives no necessity proof. Moreover, Eq. (9) is essentially a restatement of the assumption: in the expansion (A1), the first-order coefficient is by definition j[Δγ], so a non-trivial first-order coefficient is exactly j[γ]≠j[γ_HF]. The continuity-equation argument adds nothing. The 'theorem' therefore overstates its logical content.
minor comments (3)
  1. The TDHF 2RDM subtraction appears to be written as -γ_HF(r',r'')n_HF(r''); the standard factorization is Γ_HF(rr''|r'r'')=γ_HF(r,r')n_HF(r'')-γ_HF(r,r'')γ_HF(r'',r'), so the first subtraction should involve γ_HF(r,r'), not γ_HF(r',r''). Please check signs and variables.
  2. The identification 'γ_HF(t) with γ_s(t)' that yields v_T_c=0 is presented as a consequence of Eq. (10). In the context of the numerical model this is an extra assumption and should be flagged as such; without it, the numerical comparison does not verify Eq. (10).
  3. Minor wording: 'a one-dimensional helium model modeled after' is redundant; suggest 'a one-dimensional helium-like model'.

Circularity Check

2 steps flagged · score 6.0 of 10

Theorem 1's Eqs. (9) and (10) restate its own premises: the 'non-trivial first-order coefficient' is defined as j[Δγ̂], and the 'trivial second-order coefficient' becomes n∇v_T_c equality only after silently identifying γ_s with γ_HF, which the proof never establishes.

  1. self definitional [Theorem 1, Eq. (9) and Appendix Eq. (A1)]
    "Given that the exact 1RDM is non-idempotent and Δγ(r,r′,t) shows a non-trivial first-order coefficient and a trivial second-order coefficient in the relative spatial coordinate r−r′ at any time t>t0, the assumptions are necessary and sufficient for the following relations: j[γ̂(t),r] = 1/2i (∇−∇′)γ(r,r′,t)|_{r=r′} ≠ j[γ̂_HF(t),r] (9) … γ(r,r′,t) = γHF(r,r′,t) + n[Δγ̂(t),rG] + i(r−r′)·j[Δγ̂(t),rG] + ((r−r′)·(∇−∇′)/2)^2 Δγ(r,r′,t)|_{r=r′} + O(|r−r′|^3). (A1)"

    The appendixed expansion (A1) defines the first-order coefficient of Δγ = γ − γ_HF in the relative coordinate to be i(r−r′)·j[Δγ̂(t),rG]. Therefore the assumption of a 'non-trivial first-order coefficient' is exactly the statement j[Δγ̂] ≠ 0, i.e., j[γ̂] ≠ j[γ̂_HF]. Equation (9) asserts this same inequality as a conclusion. The theorem does not derive the non-HF current from the assumption; it merely restates the definition of the assumed non-trivial first-order coefficient as the conclusion.

  2. self definitional [Theorem 1, Eqs. (6), (10), and Appendix after Eq. (A1)]
    "On the other hand, from the definition in Eq. (6), n∇v_T_c corresponds to the second-order coefficient in r−r′ and equivalently leads to a trivial solution n[γ̂(t),r]∇v_T_c[γ̂(t),r] = n[γ̂_HF(t),r]∇v_T_c[γ̂_HF(t),r]. … ∇v_T_c(r,t) = −1/(4n(r,t)) (∇−∇′)(∇²−∇′²) { γ(r,r′,t) − γ_s(r,r′,t) }|_{r=r′}. (6)"

    Equation (6) defines n∇v_T_c through γ − γ_s, not through γ − γ_HF. The appendix's 'second-order coefficient' argument applies to the expansion of Δγ = γ − γ_HF in (A1), whose quadratic term is assumed trivial. For Eq. (10) to follow, one needs L(γ_s[γ]) = L(γ_s[γ_HF]) (with L the fourth-order operator in (6)); no stated premise supplies this. The equality (10) becomes equivalent to the assumed trivial quadratic coefficient only if γ_s is silently identified with γ_HF, and then it is the premise restated. The proof never mentions γ_s; the numerical section later explicitly says 'identifying γ̂_HF(t) with γ̂_s(t)', but that identification is an extra assumption not justified in TDDFT.

full rationale

The paper is not globally circular: Eq. (11) about ∇v_W_xc departing from the TDHF exchange form is a substantive claim, and the one-dimensional quench simulation is an external numerical experiment. However, the two central relations of Theorem 1 that establish the claimed hierarchy — non-HF current at first order and HF-representable kinetic component — reduce to the theorem's hypotheses by construction. Eq. (A1) defines the first-order off-diagonal coefficient as i(r−r′)·j[Δγ̂], so the premise 'non-trivial first-order coefficient' is literally Eq. (9). Similarly, Eq. (6) defines n∇v_T_c via γ−γ_s, while the proof's 'second-order coefficient' is read off the γ−γ_HF expansion; Eq. (10) follows only if γ_s=γ_HF, and in that case it restates the assumed trivial quadratic coefficient. This identification is not proven and is only admitted later in the numerical section. The independent content, Eq. (11), has a separate proof gap (induction/analyticity) but is not circular per se. Because the headline 'hierarchy' for the kinetic component is forced by definition rather than derived, the circularity score is 6.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The theorem rests on the prior decomposition of Luo et al. [20], a silent identification of the KS and TDHF 1RDMs, an asserted induction with unstated time-analyticity, and the rank-one KS ansatz of the numerics. No new physical entities (particles, forces, dimensions) are postulated. Hand-chosen model and grid parameters are listed.

free parameters (2)
  • soft-Coulomb model parameters = Z = 2 in v_ext = -2/sqrt(1+x^2); w_ee = 1/sqrt(1+(x-y)^2)
    hand-chosen model potentials for the 1D helium analog; no sensitivity study is reported.
  • numerical grid and time-step = dx = 0.1 a.u.; box x in [-15,15]; dt = 0.01 a.u.; T = 20 a.u.
    hand-chosen propagation parameters; no convergence tests or error estimates are given.
assumptions (5)
  • domain assumption The exact decomposition v_xc = v_W_xc + v_T_c with the explicit forms of Eqs. (5)-(6)
    adopted from Luo, Fuks, Sandoval, Elliott, Maitra (J. Chem. Phys. 140, 18A515 (2014), ref. [20]); the paper builds its theorem on it without re-deriving it.
  • domain assumption Kohn-Sham 1RDM gamma_s identified with the TDHF 1RDM gamma_HF in Theorem 1 and its proof
    v_T_c is defined through gamma_s in Eq. (6), but the theorem's premise and Eq. (10) are phrased through gamma_HF; the appendix switches between the two without justification. In TDDFT the KS system reproduces only the density, so this identification is non-trivial.
  • ad hoc to paper Time-analyticity of n[Δγ] and j[Δγ] in (t-t0) around t0, and the induction that all time derivatives of (n-n_HF) and (j-j_HF) vanish up to order m
    asserted in the appendix ('it can be proved inductively from m=0') without a proof; it is the load-bearing step connecting the LHS and RHS expansions that yields Eq. (11).
  • domain assumption Rank-one KS orbital ansatz phi(x,t) = sqrt(n(x,t)/2) exp(i ∫^x u(y,t) dy) for the two-electron singlet
    Eq. (12), from Wilken and Bauer (ref. [25]); the ansatz forces v_T_c = 0 by construction, so the numerical 'truncation' is a restriction on the KS reference space, not a general functional truncation.
  • standard math Standard one-body density matrix equation of motion and continuity equation
    used in the appendix (Eqs. (A2)-(A4)); standard quantum-mechanical results invoked without proof.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Spatial-Order Hierarchy of Time-Dependent Exchange-Correlation Potential." pith.science (2026). https://pith.science/paper/PQLMFK2Q

@misc{pith2026260801483,
  author       = {Pith},
  title        = {Pith review of: Spatial-Order Hierarchy of Time-Dependent Exchange-Correlation Potential},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PQLMFK2Q}},
  note         = {Machine review of arXiv:2608.01483}
}
read the original abstract

Exact time-dependent density-functional theory separates the exchange-correlation potential into interaction and kinetic-correlation components, but the structural relation between them remains unknown. We establish a representability constraint based on the off-diagonal expansion of the one-electron reduced density matrix relative to a time-dependent Hartree-Fock (TDHF) reference. A non-zero linear term generates a non-HF current density while the kinetic-correlation component remains HF representable; higher-order off-diagonal structure activates the kinetic component. In a one-dimensional two-electron correlation quench, retaining the exact interaction component while setting the kinetic component to zero suppresses the density broadening of the exact evolution. These results establish a hierarchy of density equations of motion and provide an exact constraint for non-adiabatic functional construction.

Figures

Figures reproduced from arXiv: 2608.01483 by the authors.

Figure 1
Figure 1. The electron density is denser in regions where the [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Snapshots of kinetic-correlation force [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 22 canonical work pages

  1. [1]

    Hohenberg and W

    P. Hohenberg and W. Kohn, Phys. Rev.136, B864 (1964)

    1964

  2. [2]

    Kohn and L

    W. Kohn and L. J. Sham, Phys. Rev.140, A1133 (1965)

    1965

  3. [3]

    R. G. Parr, inHorizons of Quantum Chemistry: Pro- ceedings of the Third International Congress of Quantum Chemistry Held at Kyoto, Japan, October 29-November 3, 1979(Springer, 1989) pp. 5–15

    work page 1979

  4. [4]

    Burke, J

    K. Burke, J. Chem. Phys.136, 150901 (2012)

    work page 2012

  5. [5]

    Runge and E

    E. Runge and E. K. Gross, Phys. Rev. Lett.52, 997 (1984)

    work page 1984

  6. [6]

    T. Kato, K. Nobusada, and S. Saito, J. Phys. Soc. Jpn. 89, 024301 (2020)

    work page 2020

  7. [7]

    van Leeuwen and E

    R. van Leeuwen and E. J. Baerends, Phys. Rev. A49, 2421 (1994)

    1994

  8. [8]

    E. K. U. Gross and W. Kohn, Phys. Rev. Lett.55, 2850 (1985)

    1985

Show all 27 references
  1. [9]

    M. E. Casida,Recent Developments and Applications of Modern Density Functional Theory, edited by J. Semi- nario, Theo. Comput. Chem., Vol. 4 (Elsevier, 1996) pp. 391–439

    1996

  2. [10]

    Petersilka, U

    M. Petersilka, U. J. Gossmann, and E. K. U. Gross, Phys. Rev. Lett.76, 1212 (1996)

    1996

  3. [11]

    Bauernschmitt and R

    R. Bauernschmitt and R. Ahlrichs, Chem. Phys. Lett. 256, 454 (1996)

    1996

  4. [12]

    Lein and S

    M. Lein and S. K¨ ummel, Phys. Rev. Lett.94, 143003 (2005)

    2005

  5. [13]

    Brics, J

    M. Brics, J. Rapp, and D. Bauer, Phys. Rev. A90, 053418 (2014)

    2014

  6. [14]

    Abedi, N

    A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett.105, 123002 (2010)

    2010

  7. [15]

    Schild and E

    A. Schild and E. K. U. Gross, Phys. Rev. Lett.118, 163202 (2017)

    2017

  8. [16]

    Yabana and G

    K. Yabana and G. F. Bertsch, Phys. Rev. B54, 4484 (1996)

    1996

  9. [17]

    Selstø and S

    S. Selstø and S. Kvaal, J. Phys. B: At. Mol. Opt. Phys. 43, 065004 (2010)

    2010

  10. [18]

    Kato and H

    T. Kato and H. Kono, Chem. Phys. Lett.392, 533 (2004)

    2004

  11. [19]

    Sato and K

    T. Sato and K. L. Ishikawa, J. Phys. B: At. Mol. Opt. Phys.47, 204031 (2014)

    2014

  12. [20]

    K. Luo, J. I. Fuks, E. D. Sandoval, P. Elliott, and N. T. Maitra, J. Chem. Phys.140, 18A515 (2014)

    2014

  13. [21]

    Suzuki, L

    Y. Suzuki, L. Lacombe, K. Watanabe, and N. T. Maitra, Phys. Rev. Lett.119, 263401 (2017)

    2017

  14. [22]

    T. L. Gilbert, Phys. Rev. B12, 2111 (1975)

    1975

  15. [23]

    Schilling, J

    C. Schilling, J. Chem. Phys.149, 231102 (2018)

    2018

  16. [24]

    D. A. Mazziotti,Reduced-density-matrix mechanics: with applications to many-electron atoms and molecules, Adv. Chem. Phys., Vol. 134 (John Wiley and Sons, New York, 2007)

    2007

  17. [25]

    Wilken and D

    F. Wilken and D. Bauer, Phys. Rev. Lett.97, 203001 (2006)

    2006

  18. [26]

    H. F. Trotter, Proc. Am. Math. Soc.10, 545–551 (1959)

    1959

  19. [27]

    Suzuki, Commun

    M. Suzuki, Commun. Math. Phys.51, 183–190 (1976). Appendix: Proof of the Theorem For timet > t0, the 1RDM difference ∆ˆγ(t) = ˆγ(t)−ˆγHF(t) can be expanded in the relative coordinater−r ′ as follows: γ(r,r ′, t) =γHF(r,r ′, t) +n ∆ˆγ(t),rG +i(r−r ′)·j ∆ˆγ(t),rG + (r−r ′)·(∇ − ...

    1976

Pith tools

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