{"id":"1a4aa824-bac2-4bb7-aa08-7f885edd7be7","arxiv_id":"2608.01483","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"The paper claims the kinetic-correlation part of the TDDFT exchange-correlation potential activates only beyond first order in the off-diagonal expansion of the one-electron density matrix, with a 1D two-electron test where omitting it suppresses density broadening.","lead":"Exact time-dependent density-functional theory splits the exchange-correlation potential into interaction and kinetic-correlation parts; this paper claims a spatial-order rule linking that split to the off-diagonal structure of the one-electron density matrix, and tests it on a one-dimensional two-electron quench. The rule would help non-adiabatic functional builders decide when kinetic correlation must be included, but as presented the proof and numerics do not fully support","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 1's Eq. (10) silently identifies the KS 1RDM γ_s with γ_HF; the proof never uses γ_s, so the key 'kinetic component remains HF-representable' claim does not follow from the stated premises.","rationale":"I read the paper as aiming to establish an exact structural constraint between v_W_xc and v_T_c in TDDFT. The proof of Theorem 1 is the load-bearing element. Eq. (9) follows from the premises, and the numerical demonstration has some qualitative value (though class (ii) knowingly sets γ_s=γ_HF). The central problem is Eq. (10). The proof's one-line justification ignores the γ_s term in Eq. (6). If γ_s=γ_HF is assumed, the theorem's statement should say so and the result becomes nearly tautological; if not, the equality is unjustified. This is not a consensus disagreement; it is an internal gap between Eq. (6) and the theorem. The appendix's proof of Eq. (11) also relies on an unproven induction and time-analyticity, but the γ_s issue is sufficient to invalidate the 'hierarchy' as stated. The paper's limitation paragraph is honest and the numerical setup is concrete, but these do not repair the proof. I therefore agree with the reader's REJECT. The recommended check — comparing both sides of Eq. (10) in the paper's own model with the correct γ_s — would settle the matter quickly. Since my conclusion matches the reader's verdict, I recommend no change.","tokens_in":9062,"tokens_out":11727,"duration_ms":103901,"concrete_test":"Use the paper's 1D two-electron model. Propagate exact TDSE and TDHF from the same HF initial state; construct the KS 1RDM γ_s(x,x',t) from Eq. (12) with n and u=j/n from the exact propagation. At times where j[γ]≠j[γ_HF] and the second-order coefficient of γ-γ_HF is zero, compute both sides of Eq. (10) (with γ_s[γ] and γ_s[γ_HF] defined by the corresponding densities) from polynomial fits in d=x-x' to second order. If they differ, Theorem 1 is false as stated. Independently, re-derive Eq. (10) with the γ_s terms of Eq. (6) retained, and verify whether the equality imposes an unstated condition Q_{γ_s[γ]}=Q_{γ_s[γ_HF]}.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing weakness is in Theorem 1, Eq. (10). From Eq. (6), n∇v_T_c = -1/4 L(γ - γ_s), where L=(∇-∇')(∇²-∇'²)|_{r=r'}. Eq. (10) therefore requires L(γ - γ_s[γ]) = L(γ_HF - γ_s[γ_HF]). The premise that Δγ has a trivial second-order coefficient gives L(γ-γ_HF)=0, so Eq. (10) reduces to L(γ_s[γ]) = L(γ_s[γ_HF]). No stated premise ensures this. In the appendix proof, Eq. (10) is dismissed as 'from the definition ... equivalently leads to a trivial solution'; γ_s is never mentioned. That is only true if one has already set γ_s=γ_HF. But in TDDFT the KS 1RDM is fixed by the density, and γ_s is generally not the TDHF 1RDM (for the two-electron singlet, γ_s is rank-one via Eq. (12) but with a different orbital than γ_HF). Absent this identification, Eq. (10) — the claim that kinetic correlation remains HF-representable at linear off-diagonal order — is not a consequence of the theorem's hypotheses. The proof of Eq. (11) has a separate gap (the asserted induction and analyticity in the appendix), but Eq. (10) already blocks the central hierarchy claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":9296,"tokens_out":12652,"duration_ms":111760,"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":[{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"minor_comments":[{"comment":"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.","section":null},{"comment":"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).","section":null},{"comment":"Minor wording: 'a one-dimensional helium model modeled after' is redundant; suggest 'a one-dimensional helium-like model'.","section":null}],"recommendation":"reject","confidential_remarks":"The main theorem's proof gap is severe and not a local fix: the γ_s=γ_HF identification is invalid in general TDDFT, so Eq. (10) would require a different theorem. The numerical demonstration could potentially be developed into a separate computational study, but as submitted the central claim is not established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read the Negishi paper. The short take: the spatial-order hierarchy is a nice idea, but the proof of the main claim doesn't hold. The new ingredient is the explicit claim that a non-trivial linear off-diagonal term in the exact 1RDM relative to TDHF forces the current to differ while leaving the kinetic-correlation potential 'HF representable,' and that only quadratic terms activate v_T_c. That's a clean statement, and the numerical example—one-dimensional two-electron quench—does show that dropping v_T_c suppresses density broadening. Credit where due: the decomposition follows Luo et al., the setup is careful, and the conclusion includes an honest limitation paragraph.\n\nThe soft spot is load-bearing. Eq. (10) claims n∇v_T_c[γ] = n∇v_T_c[γ_HF]. By Eq. (6), n∇v_T_c involves the difference γ - γ_s, where γ_s is the KS 1RDM. The theorem's premises only constrain γ - γ_HF. To get Eq. (10) you need L(γ_s[γ]) = L(γ_s[γ_HF]), and nothing in the hypotheses gives that. The appendix never mentions γ_s; it just says 'from the definition' the term is trivial. That's true only if you have already set γ_s = γ_HF, which is not a valid TDDFT identification. The KS system reproduces the density, not the off-diagonal 1RDM, and γ_s is generally different from the TDHF density matrix. Without this identification, the hierarchy claim doesn't follow.\n\nThe proof of Eq. (11) has its own gap: an induction is asserted without proof, time-analyticity is assumed, and the same γ_s/γ_HF swap appears in rewriting the right-hand side. The numerics are one trajectory with no convergence tests, and the claim that class (ii) matches TDHF is unquantified.\n\nSo the paper is not ready as is. I don't think it's crank material—the question is real, and classifying by off-diagonal order is worth pursuing. But the current version doesn't establish the theorem. A serious referee could be useful if the editor expects heavy revision and a real proof. I'd send it to a referee who can check the appendix, not desk-reject it outright.","headline":"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.","tokens_in":9924,"tokens_out":4521,"would_cite":false,"duration_ms":39523,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["time-dependent density functional theory","exchange-correlation potential","kinetic-correlation potential","reduced density matrix","off-diagonal expansion","TDHF reference","non-adiabatic functional","spatial-order hierarchy"],"falsifier":"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.","tokens_in":8741,"feed_emoji":"⚛️","tokens_out":4785,"duration_ms":41478,"temperature":0.7,"pith_summary":"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.","feed_headline":"Off-diagonal order sets when kinetic correlation activates","feed_subtitle":"Linear deviations drive non-HF currents; quadratic deviations switch on kinetic correlation.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduces the separation of the exchange-correlation potential into interaction and kinetic-correlation components and demonstrates their distinct dynamical roles numerically; the theorem extends this to a proven constraint.","marker":"[20]"},{"why":"Provides exact structural decompositions of the time-dependent exchange-correlation potential that the hierarchy builds on.","marker":"[21]"},{"why":"Runge-Gross theorem, the foundation that the TDKS density evolves exactly and the framework for the density equation of motion.","marker":"[5]"},{"why":"Kohn-Sham theory, defining the KS 1RDM $\\gamma_s$ used in the definition of $v^T_c$.","marker":"[2]"},{"why":"Supplies the rank-one TDKS orbital form for a one-dimensional two-electron system used in the numerical test.","marker":"[25]"},{"why":"Establishes the 1RDM as the central object whose off-diagonal structure encodes quantum non-local correlation.","marker":"[22]"}],"fun_headline_variants":["Kinetic correlation turns on at second order in 1RDM","Spatial-order hierarchy links current and kinetic terms","Exact constraint: kinetic part activates at quadratic order","Off-diagonal order gates kinetic correlation","Hierarchy: interaction first, kinetic second"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Kinetic correlation turns on at second order in 1RDM","Spatial-order hierarchy links current and kinetic terms","Exact constraint: kinetic part activates at quadratic order","Off-diagonal order gates kinetic correlation","Hierarchy: interaction first, kinetic second"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000149,"raw_usage":{"total_tokens":977,"prompt_tokens":642,"completion_tokens":335,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":386,"completion_tokens_details":{"reasoning_tokens":261}},"tokens_in":386,"tokens_out":335,"duration_ms":3939,"temperature":1.0,"reasoning_tokens":261,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T00:08:47.088408+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the separation of the exchange-correlation potential into interaction and kinetic-correlation components and demonstrates their distinct dynamical roles numerically; the theorem extends this to a proven constraint."},{"cited_title":"Suzuki, L","cited_arxiv_id":null,"evidence_quote":"Provides exact structural decompositions of the time-dependent exchange-correlation potential that the hierarchy builds on."},{"cited_title":"Runge and E","cited_arxiv_id":null,"evidence_quote":"Runge-Gross theorem, the foundation that the TDKS density evolves exactly and the framework for the density equation of motion."},{"cited_title":"Wilken and D","cited_arxiv_id":null,"evidence_quote":"Supplies the rank-one TDKS orbital form for a one-dimensional two-electron system used in the numerical test."}],"review_version":1}