{"id":"b3da340e-65e4-417d-8035-c45ef1825952","arxiv_id":"2507.12390","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The rescaled Schrödinger dynamics of N dense fermions with C^2 pair interactions converges to the time-dependent Hartree dynamics with explicit rate N^{-1/24}.","lead":"A rigorous proof shows that for N fermions packed into a fixed volume with strong pair interactions, the full many-body Schrödinger dynamics converges to the simpler Hartree mean-field equations on the natural rescaled time scale. The new ingredient is a time-dependent gauge transformation that removes the large interaction term, a technique that may transfer to other strongly interacting many-body models.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The gauge transformation (2.11) with f=-∇v produces H_g=∑(i∇-tε∑f)^2, not (2.12); likewise (2.14) does not give (2.16). The proof estimates the wrong gauged dynamics, so the final comparison to Hartree evolution is not established.","rationale":"The reader's verdict is CONDITIONAL, and I agree that the paper should not be rejected outright: the proof strategy is coherent and the estimates are mostly sign-insensitive, so a systematic correction of the gauge-phase sign could repair the argument. However, I do not share the reader's assessment that the sign inconsistency in Eqs. (2.11)–(2.12) and (2.14)–(2.16) is merely a harmless typo. Because the entire proof is built on the gauged Hamiltonians, a mismatch between the gauge transformation and the generator means the estimates control a different dynamics, not the one linked to the physical Schrödinger and Hartree evolutions. The concrete differentiation test would settle whether this is a purely typographical issue or a genuine obstruction. The C^2 regularity restriction and the exclusion of Coulomb are properly acknowledged limitations, not internal errors. The final verdict remains CONDITIONAL pending correction of the gauge identities.","tokens_in":52348,"tokens_out":43234,"duration_ms":441587,"concrete_test":"Perform the explicit differentiation of (2.11): for H=-Δ+V, compute i∂t(e^{itεV}e^{-itεH}Φ_0) and compare the coefficient of the cross term tε∑_{j≠i}f_{ij}·∇_i with the corresponding term in ε∑(i∇_i+tε∑_{j≠i}f_{ij})^2. The signs will differ, proving that (2.12) corresponds to the opposite gauge phase. Also substitute the gauge relation (2.14) into (2.15)–(2.16) for a simple smooth v and Gaussian initial orbitals; if the equation fails, the gauged Hartree orbitals used in the counting estimates are not the gauge transforms of the true Hartree orbitals.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Direct differentiation of (2.11) with V=∑_{i<j}v(x_i-x_j) gives i∂tΨ_t = e^{itεV}ε(-Δ)e^{-itεV}Ψ_t = ε∑_i(i∇_i + tε∇_i V)^2Ψ_t = ε∑_i(i∇_i - tε∑_{j≠i}f_{ij})^2Ψ_t, because f=-∇v. Eq. (2.12) instead has (i∇_i + tε∑_{j≠i}f_{ij})^2, which is the generator for the opposite phase e^{-itεV}Φ_t. Similarly, differentiating ψ_t=e^{itε(v*ρ_t)}φ_t using (1.6) yields a gauged Hartree generator with magnetic term (i∇-tε\\bar f)^2 and with the scalar f·f term entering with the opposite sign to (2.16); the claimed identity (2.19) is consistent with neither sign. All subsequent constructions—H_g, ~H_g, R(t), W(t), and the cancellation estimates in §§3–4—are built on the wrong generators. Hence the bound (4.30) concerns a wave function that is not the gauge transform of Φ_t, and the equality Tr(Mγ^{Φ_t}) - 1/N Tr(Mp^{φ_t}) = Tr(Mγ^{Ψ_t}) - 1/N Tr(Mp^{ψ_t}) in (4.33) is not justified. The theorem as written is therefore not proved. The error is systematic and likely repairable by flipping the gauge-phase sign consistently, but the displayed identities must be corrected before acceptance.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims a rigorous derivation of the time-dependent Hartree equations for N fermions in a volume of order one with an N-independent, strongly interacting C^2 radial pair potential, on the rescaled time scale t = tau N^{-2/3}. The main result, Theorem 1.1, asserts that expectation values of bounded multiplication operators in the microscopic Schr\\\"odinger state are approximated by the corresponding expectation values in the Slater determinant of Hartree orbitals, up to an error of order exp(C e^{(1+t)^2}) times a vanishing power of N, provided the initial data satisfy Assumptions 1.2 and (1.10). The proof proceeds through time-dependent gauge transformations (2.11) and (2.14), an auxiliary quadratic Hamiltonian \\tilde H_g, counting-functional and bad-particle kinetic estimates for the auxiliary dynamics, a norm approximation between the gauged and auxiliary evolutions, and a final reduction from counting estimates to one-particle reduced density estimates. The logical structure is clear and the technical machinery is substantial, but a sign inconsistency in the gauge transformation affects the central comparison and must be corrected before the theorem can be regarded as proved.","tokens_in":52667,"tokens_out":6903,"duration_ms":82406,"significance":"If the proof is repaired, the result would be a significant advance: it would provide the first derivation of fermionic mean-field dynamics in a dense, strongly interacting, non-semiclassical regime in which kinetic and interaction terms both contribute at leading order. The gauge-removal strategy and the auxiliary Bogoliubov-type quadratic approximation are original and appear well suited to the strong-coupling problem. The paper also gives explicit convergence rates and clearly identifies the physical restrictions: the pair potential must be C^2 and radial, the Coulomb potential is excluded, and the approximation is stated only for multiplication observables rather than in trace norm. The authors are also transparent about the fact that singular potentials |x|^{-s} with s<5/8 would require additional technical assumptions. These strengths are conditional, however, because the gauge sign error invalidates the comparison of the gauged and ungauged dynamics as written.","major_comments":[{"comment":"The gauge transformation is sign-inconsistent, and the subsequent proof estimates the wrong dynamics. Direct differentiation of Ψ_t = exp(it ε Σ_{i<j} v(x_i-x_j)) Φ_t using i∂_t Φ_t = ε H Φ_t gives i∂_t Ψ_t = ε Σ_i (i∇_i - tε Σ_{j≠i} f_{ij})^2 Ψ_t, because ∇_i Σ_{j≠i} v(x_i-x_j) = -Σ_{j≠i} f_{ij} with f = -∇v. Equation (2.12) instead has (i∇_i + tε Σ_{j≠i} f_{ij})^2. Likewise, differentiating ψ_t^k = exp(it ε (v*ρ_t)) φ_t^k using (1.6) yields a gauged Hartree generator with (i∇ - tε \\bar f)^2, not the plus sign in (2.16), and the identity (2.19) has the opposite sign as well. Since the auxiliary Hamiltonian (3.10), the counting-functional estimates in §3, and the norm approximation in §4 are all built on the generators (2.12) and (2.16), the bound (4.4) and the equality (4.33) are not established for the wave functions defined in (2.11) and (2.14). This is a load-bearing error rather than a typographical slip: the sign propagates through (2.19), (3.10), (3.70), and the estimates of Lemmas 3.3–3.7. The authors should correct the gauge convention consistently—for example by using the phase exp(-it ε V) in (2.11) and (2.14) or by changing the sign of f throughout—and re-verify all affected displayed identities and estimates.","section":"§2.2, Eqs. (2.11)–(2.16)"}],"minor_comments":[{"comment":"In the second displayed line of the computation, the term t²ε²(~w∇f)_ij should presumably be t²ε²(~wf)_ij; as printed it duplicates the preceding term and omits the f·f interaction term that appears in the definition of \\tilde H_g in (3.10).","section":"§3.1, Eq. (3.20)"},{"comment":"The notation ⟨w^{(γ)}_{-1}Ψ_t, ...⟩ is missing the hat on the weight operator; the context indicates it should read ⟨\\hat w^{(γ)}_{-1}Ψ_t, ...⟩.","section":"§4, proof of Lemma 4.1, Eq. (4.25)"},{"comment":"The labels (IIIc1) and (IIIc1) are duplicated in the two successive estimates; the second should be labeled differently, for instance (IIIc2), to keep the enumeration consistent.","section":"§4, proof of Lemma 4.1, Eqs. (4.26)–(4.27)"},{"comment":"The restriction to C^2 radial pair potentials and the explicit exclusion of the Coulomb potential are important physical limitations; since the introduction emphasizes applications such as electrons in molecules and dense matter, the abstract or introduction should state clearly that the Coulomb case is not covered and that the paper only treats bounded forces.","section":"§1.2, Remark after Theorem 1.1"}],"recommendation":"major_revision","confidential_remarks":"The sign error in §2.2 is systematic and affects the central claim, so the paper cannot be accepted in its current form. The error appears repairable by a consistent change of gauge-phase sign, and the rest of the proof structure is coherent enough that a corrected version could be publishable. I would ask the authors to re-derive the gauged generators, update all dependent identities, and confirm that the bounds in §§3–4 go through unchanged or with the necessary modifications."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read the paper by Hoang, Mitrouskas, and Pickl on deriving time-dependent Hartree equations for dense strongly interacting fermions. The headline: this is the first rigorous mean-field derivation in a non-semiclassical regime where interactions are leading order, and the time-dependent gauge transformation is genuinely new. The proof is long but carefully structured: counting-functional bounds, an auxiliary quadratic dynamics, norm approximation, and a final trace-norm comparison. I found no issue with the overall architecture, and the cited tools from Petrat–Pickl are appropriate.\n\nWhat the paper does well: the heuristics in the introduction (Ehrenfest scaling, localized vs. delocalized orbitals) are physically illuminating. The gauge idea—removing the O(N^{4/3}) potential by conjugating the kinetic operator rather than the interaction—is clever and seems transferable. The technical lemmas in Section 5 are detailed and appear internally consistent.\n\nNow the soft spots. The most serious is a systematic sign inconsistency. The gauge in (2.11) is e^{+itεV}; differentiating it with V = ∑_{i<j} v(x_i-x_j) and f = -∇v gives H_g = ∑(i∇_i - tε∑_{j≠i} f_{ij})^2, not the plus sign in (2.12). Similarly, (2.14) produces a Hartree generator with (i∇ - tε \\bar f)^2, not the plus sign in (2.16). This matters: all subsequent estimates—the auxiliary Hamiltonian, the counting-functional derivatives, the norm approximation—are built on the displayed (wrong) generator. As written, the proof compares the true Φ_t to a Slater state via a wave function Ψ_t that is not its gauge transform. The good news: the error looks global and sign-insensitive. The bounds in the paper use norms and absolute values of f, so flipping the sign of f everywhere (or equivalently using e^{-itεV}) should leave the estimates unchanged. So I expect the argument is repairable, but the displayed identities must be corrected before the theorem can be accepted.\n\nThe restriction to C^2 radial potentials is real—Coulomb is excluded—but the authors are transparent about it and mention an extension to |x|^{-s}, s<5/8. The final bound has a double exponential in t, so the effective time window is tiny, but that is par for the course in this kind of work.\n\nThis paper deserves a serious referee. It is a substantial advance, and the flaw, while load-bearing at the level of displayed equations, appears to be a fixable sign convention. I would send it back for careful revision and re-checking of Sections 2–4, not desk-reject it.","headline":"Dense strongly-interacting fermionic mean-field derivation with a novel gauge method, but a systematic sign error in the displayed gauge dynamics must be corrected before the theorem as stated is proved.","tokens_in":53250,"tokens_out":5375,"would_cite":true,"duration_ms":58711,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V70","35Q41","81Q05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that the time-dependent Hartree equations are the large-$N$ limit of the $N$-fermion Schr\\\"odinger dynamics in a dense, strongly interacting regime, and does so through a time-dependent gauge transformation that removes…","keywords":["time-dependent Hartree equations","fermionic mean-field dynamics","dense strongly interacting fermions","gauge transformation","counting functional","Schrodinger dynamics","Slater determinant approximation","many-body quantum dynamics"],"falsifier":"Run the microscopic Schr\\\"odinger evolution numerically for a dense system of $N$ fermions with a compactly supported, smooth radial pair potential, starting from an exact Slater determinant built from orbitals satisfying Assumption 1.2, and compare the one-particle density with the Hartree prediction at fixed rescaled time: the theorem predicts the difference of expectation values decays roughly like $N^{-1/24}$, so a clearly different decay rate or non-convergence would refute the claimed bound. A second, more targeted test uses a singular potential such as $v(x)=|x|^{-1/2}$, which satisfies $s=1/2<5/8$ and which the paper says the proof can handle with modifications; failure of convergence there would pinpoint where the $L^\\infty$ force estimates are the load-bearing step.","tokens_in":52103,"feed_emoji":"⚛️","tokens_out":13624,"duration_ms":129766,"temperature":0.7,"pith_summary":"The paper aims to prove that the time-dependent Hartree equations\\u2014the standard single-particle mean-field description of fermions in atoms, molecules, and nuclei\\u2014follow rigorously from the microscopic Schr\\\"odinger dynamics of $N$ fermions in a regime where the interaction is strong, not weak. The setting is a dense gas: $N$ fermions in a volume of order one, governed by the rescaled Hamiltonian $\\varepsilon H$ with $\\varepsilon = N^{-2/3}$, on time scales on which each particle feels a force of order $N$ and moves a distance of order $N^{-1/3}$. Earlier derivations of fermionic mean-field dynamics required either weak coupling, where the interaction is subleading, or a semiclassical limit; here the pair potential is $N$-independent and contributes at leading order together with the kinetic energy. The central device is a time-dependent gauge transformation that removes the large interaction potential from both the Schr\\\"odinger and Hartree generators, yielding magnetic-type terms that are effectively of order one per particle. If the theorem is right, Hartree dynamics\\u2014not free or Vlasov evolution\\u2014is the correct effective theory on the natural short time scale of dense fermionic systems, including localized orbitals whose motion amounts to genuine macroscopic quantum transport.","feed_headline":"Gauge trick yields Hartree equations for dense fermions","feed_subtitle":"No weak-coupling assumption needed: the N-body Schrödinger state tracks the Hartree mean-field state as N grows.","key_machinery":"The load-bearing mechanism is a time-dependent gauge transformation: multiplying the Schr\\\"odinger wave function by the phase $\\exp\\bigl(i t \\varepsilon \\sum_{i<j} v(x_i-x_j)\\bigr)$ removes the large potential from the microscopic Hamiltonian and replaces it with magnetic-type kinetic terms $(i\\nabla_i + t\\varepsilon \\sum_{j\\ne i} f_{ij})^2$, with $f = -\\nabla v$, which carry extra factors of $\\varepsilon$ and are effectively of order one per particle; the identical transformation is applied to the Hartree orbitals. Since the old counting-functional method would produce an $O(N^{1/3})$ growth rate, the proof introduces an auxiliary Hamiltonian $\\widetilde{H}_g(t)$\\u2014the quadratic approximation to the gauged generator obtained by discarding all terms with three or more $q$-projections\\u2014and controls the number of 'bad' particles outside the Hartree orbitals via counting functionals with weight functions $m^{(\\gamma)}(k)=\\min\\{1, k/N^\\gamma\\}$ and $w^{(\\gamma)}=1-m^{(\\gamma)}$. Diagonalization estimates, which express operators of the form $p_2 h_{12} p_2$ in a basis where they are diagonal and subtract the mean-field contribution, provide the cancellations needed for the zero- and one-excitation terms.","core_discovery":"The central claim (Theorem 1.1) is quantitative: for a real-valued, radial pair potential $v \\in C^2(\\mathbb{R}^3)$ and initial data satisfying Assumptions 1.2 and (1.10), for all bounded multiplication operators $M$, $$\\sup_{\\|M\\|\\le 1} \\left| \\operatorname{Tr}(M\\$gamma^{{\\Phi_t}}$) - \\frac{1}{N}\\operatorname{Tr}(M $p^{{\\varphi_t}}$) \\right| \\le \\exp\\bigl(C $e^{{(1+t)^2}}$\\bigr) \\max\\bigl\\{ $N^{{5/24-\\delta_1/4}}$,\\, $N^{{1/12-\\delta_2/4}}$,\\, $N^{{1/12-\\delta_1/8}}$,\\, $N^{{-1/24}}$ \\bigr\\},$$ which vanishes as $N\\to\\infty$ for fixed $t$, with $\\gamma^{\\Phi_t}$ the one-particle reduced density of the Schr\\\"odinger evolution and $p^{\\varphi_t}$ the projector onto the Hartree orbitals. The orbitals solve the rescaled Hartree equations $i\\partial_t \\varphi_t^k = \\varepsilon(-\\Delta + v\\ast \\rho_t)\\varphi_t^k$ with density $\\rho_t = \\sum_{k=1}^N |\\varphi_t^k|^2$. The paper presents this as the first derivation of fermionic mean-field dynamics in which both quantum effects and inter-particle forces are leading order, with no $N$-dependent coupling constant. For exact Slater initial data the convergence rate becomes $N^{-1/24}$.","pith_inferences":["Because the error bound carries a factor $\\exp(C e^{(1+t)^2})$, the theorem guarantees convergence only for fixed $t$ as $N\\to\\infty$; an open question the paper does not address is whether the valid time interval can grow with $N$ while keeping the error small.","The theorem is deliberately restricted to multiplication operators, since these commute with the gauge transformation; an extension to trace-norm closeness of full density matrices would require controlling the gauged dynamics against momentum-sensitive observables, which the present estimates do not yet reach.","The regularity threshold is the most exposed boundary: the proof needs $L^\\infty$ control of $f$ and $\\nabla f$, so potentials with singularities in the range the paper says are reachable would cleanly test where the diagonalization machinery breaks down.","The gauge-elimination strategy is in principle model-agnostic: any strongly interacting Hamiltonian whose dominant potential can be absorbed by a phase is a candidate for the same treatment, provided the generated magnetic terms can be controlled."],"forward_implications":["For dense fermionic systems with $C^2$ pair potentials, the true $N$-body state remains close to a Slater determinant of Hartree orbitals for times $t=O(1)$ in the rescaled variable, with error vanishing like $N^{-1/24}$ for exact Slater initial data.","The approximation holds for both delocalized orbitals, whose density varies on macroscopic scales, and localized orbitals, whose density varies on microscopic scales; in the localized case the Hartree dynamics produces observable macroscopic transport.","The strong-interaction regime is distinct from the semiclassical one: no small Planck constant and no weak coupling are needed, so the Hartree equations, rather than the Vlasov equation, are the correct leading-order effective dynamics.","The gauge transformation plus quadratic Hamiltonian provides a concrete analytic template that the authors indicate extends to singular potentials $|x|^{-s}$ with $s<5/8$, with the Coulomb potential flagged as future work."],"supporting_citations":[{"why":"Supplies the counting-functional method, the weight-function estimates, and the old scaling that the paper shows fails for strong interactions.","marker":"[35]"},{"why":"Introduces the notion of dense, strongly interacting fermionic systems that defines the regime studied here.","marker":"[17]"},{"why":"Provides global existence and uniqueness of the Hartree orbital evolution used in the theorem.","marker":"[7]"},{"why":"Supplies kinetic-energy estimates and well-posedness for Hartree solutions relied on throughout the analysis.","marker":"[1]"},{"why":"The well-posedness theorem for non-autonomous evolution equations invoked to define the auxiliary gauged dynamics.","marker":"[39]"},{"why":"The diagonalization lemmas for mean-field fermion derivations that Section 5.2 generalizes for the cancellation estimates.","marker":"[33]"},{"why":"The weakly coupled fermion Hartree derivation whose regime this work explicitly contrasts with the strong-interaction setting.","marker":"[14]"}],"fun_headline_variants":["Gauge trick proves Hartree for strongly interacting fermions","No weak coupling needed: fermionic Hartree derived","Strong-interaction fermions: Hartree equations in large-N","Gauge transform yields fermionic mean-field limit","Dense fermions: Hartree equations derived rigorously"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the regularity of the pair potential: $v$ must be real-valued, radial, and twice continuously differentiable, so that the force $f=-\\nabla v$ and its derivative are bounded; the diagonalization and cancellation estimates rely on that boundedness, which is why the Coulomb potential is explicitly excluded, with $|x|^{-s}$, $s<5/8$, noted as reachable only through technical modifications.","fun_headline_variants_meta":{"raw":{"variants":["Gauge trick proves Hartree for strongly interacting fermions","No weak coupling needed: fermionic Hartree derived","Strong-interaction fermions: Hartree equations in large-N","Gauge transform yields fermionic mean-field limit","Dense fermions: Hartree equations derived rigorously"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000881,"raw_usage":{"total_tokens":3820,"prompt_tokens":972,"completion_tokens":2848,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":588,"completion_tokens_details":{"reasoning_tokens":2770}},"tokens_in":588,"tokens_out":2848,"duration_ms":26213,"temperature":1.0,"reasoning_tokens":2770,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:47:35.968995+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the microscopic Schr\\\"odinger evolution numerically for a dense system of $N$ fermions with a compactly supported, smooth radial pair potential, starting from an exact Slater determinant built from orbitals satisfying Assumption 1.2, and compare the one-particle density with the Hartree prediction at fixed rescaled time: the theorem predicts the difference of expectation values decays roughly like $N^{-1/24}$, so a clearly different decay rate or non-convergence would refute the claimed bound. A second, more targeted test uses a singular potential such as $v(x)=|x|^{-1/2}$, which satisfies $s=1/2<5/8$ and which the paper says the proof can handle with modifications; failure of convergence there would pinpoint where the $L^\\infty$ force estimates are the load-bearing step.","supporting_citations":[{"cited_title":"Petrat and P","cited_arxiv_id":null,"evidence_quote":"Supplies the counting-functional method, the weight-function estimates, and the old scaling that the paper shows fails for strong interactions."},{"cited_title":"Fournais, B","cited_arxiv_id":null,"evidence_quote":"Introduces the notion of dense, strongly interacting fermionic systems that defines the regime studied here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides global existence and uniqueness of the Hartree orbital evolution used in the theorem."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies kinetic-energy estimates and well-posedness for Hartree solutions relied on throughout the analysis."},{"cited_title":"Schmid and M","cited_arxiv_id":null,"evidence_quote":"The well-posedness theorem for non-autonomous evolution equations invoked to define the auxiliary gauged dynamics."},{"cited_title":"Petrat,Derivation of mean-field dynamics for fermions, Ph.D","cited_arxiv_id":null,"evidence_quote":"The diagonalization lemmas for mean-field fermion derivations that Section 5.2 generalizes for the cancellation estimates."},{"cited_title":"Elgart, L","cited_arxiv_id":null,"evidence_quote":"The weakly coupled fermion Hartree derivation whose regime this work explicitly contrasts with the strong-interaction setting."}],"review_version":1}