{"id":"6b51f9bb-7f1f-41d5-a386-c521575c7b5f","arxiv_id":"2411.14028","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The Cauchy problem for the time-dependent Bogoliubov-Dirac-Fock evolution of graphene has a unique global solution for Fermi velocities above a critical threshold.","lead":"This mathematics paper proves global well-posedness for the time-dependent Hartree-Fock equations describing electrons in graphene, allowing for time-varying external electric potentials from local charge defects. The result extends earlier well-posedness theorems from three-dimensional quantum electrodynamics to the two-dimensional massless Dirac setting of graphene, under a condition on the Fermi velocity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Global well-posedness for graphene hinges on the numerical estimate vc≈0.36; the rigorous bound only gives vc<2.056, so the theorem's hypothesis is not proven for vF≈1.1.","rationale":"I focused on the numerical threshold vc rather than the missing projection hypothesis because the latter is a statement-level omission: physical initial data are always projections, Lemma 3.3 already proves preservation, and adding 'γI² = γI' fixes it. The vc issue instead determines whether the theorem's hypothesis is satisfied at the physical Fermi velocity vF ≈ 1.1. The paper is transparent about the numerical nature of vc, but the abstract and title do not carry the caveat, so the central claim about graphene is only as strong as the numerical estimate. This agrees with the reader's weakest_assumption. If a high-precision computation confirmed h(1.1) > 2, the main theorem would remain true but would not establish well-posedness for graphene, requiring a substantial downgrade of the paper's physical claims; if it confirms h(1.1) ≤ 2, the concern is resolved and the current conditional verdict is appropriate.","tokens_in":22427,"tokens_out":22724,"duration_ms":212061,"concrete_test":"Evaluate h(1.1) for the physical cutoff Λ ≈ 0.1 Å^{-1} (or equivalently verify the value vc ≈ 0.36) using high-precision/rigorous numerics for the variational problem (15), e.g., interval-arithmetic computation of the sup over φ. If h(1.1) ≤ 2, the theorem's hypothesis is satisfied for graphene and the concern is resolved; if h(1.1) > 2, Theorem 2.7 does not apply to graphene and the abstract must be revised to state the condition vF ≥ vc explicitly. A cheaper cross-check is to compute αc = 1/vc and test whether 1/vF ≈ 0.909 ≤ αc.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The theorem's hypothesis vF ≥ vc = h^{-1}(2) (Thm 2.7) is the only point where graphene's physical parameters enter. Global existence relies on coercivity inequality (48) in Lemma 3.4, which requires h(vF) ≤ 2 so that the exchange term is dominated by the kinetic term. The paper's own rigorous estimate is only vc < 2.0560 (αc > 0.48637, stated in the introduction), while the value vc ≈ 0.36 that would include graphene with vF ≈ 1.1 is a numerical computation reported in [34] and not proved here. Since the theorem is conditional on this numerical value, the abstract's 'global well-posedness result for graphene' is not a rigorously established statement about graphene unless h(1.1) ≤ 2; the paper's rigorous bound does not imply it. This makes the numerical estimate load-bearing for the central physical claim, even though the theorem itself is a valid conditional result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a time-dependent Hartree-Fock (Bogoliubov-Dirac-Fock) model for graphene with a sharp ultraviolet cutoff Λ, including a time-dependent external Coulomb potential generated by a charge density ν(t). The main result, Theorem 2.7, states that if the Fermi velocity satisfies vF ≥ vc = h^{-1}(2), then the Cauchy problem (28) for the evolution of the one-body density matrix has a unique global solution γ(t) ∈ C^1([0,∞), P_-^0 + Y). The proof proceeds through four steps: a local Lipschitz analysis of the mean-field vector field in the space Y (Lemma 3.2), preservation of the projection property (Lemma 3.3), a Grönwall-type bound for the BDF energy with a correcting term (Lemma 3.4), and the resulting global bound on the Y-norm (Lemma 3.5). The argument is a 2D adaptation of the 3D QED framework of Hainzl–Lewin–Sparber.","tokens_in":22658,"tokens_out":15334,"duration_ms":123859,"significance":"The result is a rigorous global well-posedness statement for a non-perturbative mean-field evolution relevant to graphene, extending the BDF framework to two dimensions with massless Dirac dispersion and time-dependent external fields. The proof is detailed and uses the coercivity of the BDF energy, avoiding smallness assumptions on the initial data. The reliance on established results from [34] and [33] is transparent, and the paper explicitly identifies the advantage over nonlinear Dirac models. The main caveat is that the theorem's hypothesis on vF is not rigorously verified for the physical value vF≈1.1; the rigorous bound vc<2.056 is insufficient, and the condition relies on a numerical estimate from [34]. Additionally, the critical case vF=vc is not handled by the proof.","major_comments":[{"comment":"The theorem is stated under the condition vF ≥ vc = h^{-1}(2), but the paper only proves the rough bound vc < 2.056 (Section 1, p.3) and cites a numerical estimate vc ≈ 0.36 from [34]. Since the physical Fermi velocity of graphene is vF ≈ 1.1, it is not rigorously established that the theorem applies to graphene. The abstract and title claim a 'global well-posedness result for graphene' without this caveat. Please state explicitly that the physical regime is covered only under the numerically supported assumption h(1.1) ≤ 2, or strengthen the rigorous bound.","section":"Section 1 and Theorem 2.7"},{"comment":"The inequality (48) contains the factor (1 - h(vF)/2). At the endpoint vF = vc, one has h(vF)=2 and this factor vanishes, so the energy bound from Lemma 3.4 does not control tr(D0 Q). Lemma 3.5 then cannot bound the Schatten norms of Q(t) as claimed. The theorem's hypothesis vF ≥ vc therefore needs to be strengthened to vF > vc, or a separate argument must be supplied for the critical case.","section":"Theorem 2.7 and Lemmas 3.4–3.5"},{"comment":"The estimate for the S2 component of [V'_Q, P_-^0] is displayed as || [V'_Q,P_-^0] |D0| ||_{S2} ≤ 2 ||D0||_∞ || Q(x,y)/|x-y| ||_{S2}. The norm required by (30) is || |D0|^{1/2} [V'_Q,P_-^0] ||_{S2}; the displayed bound is not sufficient because |D0| may vanish near p=0. Replacing |D0| by |D0|^{1/2} and ||D0||_∞ by || |D0|^{1/2} ||_∞ gives a valid estimate by the same Hardy argument. Please correct and confirm that the subsequent bounds go through.","section":"Lemma 3.2, Step II"}],"minor_comments":[{"comment":"The hypothesis on γI should explicitly require that γI is an orthogonal projection, i.e., γI^2 = γI (and 0 ≤ γI ≤ 1), since Lemma 3.3 uses A(0)=γI^2-γI=0 and the system (28) imposes γ(t)^2 = γ(t).","section":"Theorem 2.7 and Lemma 3.3"},{"comment":"The second term in (2) is written -iσ1∂x2; it should be -iσ2∂x2.","section":"Equation (2)"},{"comment":"The text says 'where PI is the initial datum'; this should be γI.","section":"After Eq. (28)"},{"comment":"The notation P(t) is used without definition; it appears to mean γ(t) = Q(t) + P_-^0. Please define it.","section":"Section 3, proof of Lemma 3.4"},{"comment":"The paper uses D0 for both the free Dirac operator (2) and the effective operator (10); the distinction is occasionally blurred. Consider using a different symbol, e.g., \\mathcal D_0, for the effective operator.","section":"Notation"},{"comment":"The title contains 'HAR TREE-FOCK' due to a line break; ensure the final version uses the correct spacing.","section":"Typesetting"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid adaptation of the BDF framework to graphene; the main theorem is plausible and the proof is detailed. The two substantive issues are the critical-case gap and the unverified numerical hypothesis for the physical regime. Both are fixable without changing the core argument. The dependence on [34] is heavy but acknowledged. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this is a competent extension of the 3D massive BDF evolution results to the 2D massless graphene model, and the main theorem is valid as stated with its hypothesis. The one thing to keep in mind is that the \"graphene\" part of the abstract is not fully covered by the proof: the theorem requires vF ≥ vc = h^{-1}(2), and for graphene's vF ≈ 1.1 there is only the numerical estimate vc ≈ 0.36 from [34]. The rigorous bound vc < 2.056 is too weak to imply it. The authors say exactly this in the introduction, so the mathematical statement is honest; the abstract is not.\n\nWhat is actually new: this is the first global well-posedness result for the time-dependent Hartree-Fock/BDF equation in 2D massless Dirac graphene with time-dependent external potentials. It adapts the strategy of [33,41] from massive 3D to massless 2D, and the main nontrivial new ingredient is handling the momentum-dependent Fermi velocity veff(p) and the associated Kohn anomaly. The proof is standard but careful: local Lipschitz flow in the Y norm, preservation of the projector structure, and an energy bound with Grönwall's lemma. I see no circularity: the coercivity of the energy comes from [34], and the evolution framework from [33], both independent of the present claim. The numerics are not fitted here, just imported from [34].\n\nSoft spots: Theorem 2.7 omits a hypothesis on the initial datum—it should say explicitly that γI is an orthogonal projection (equivalently −P0− ≤ QI ≤ P0+), not just QI ∈ Y. The proof uses it and the statement needs it. There are a few typographical slips; Step II of Lemma 3.2 has a mismatched heading. These are minor. The main limitation is the physical interpretation: the rigorous theorem applies to a range of parameters that is not shown to include graphene. The numerical estimate vc≈0.36 is load-bearing for the abstract's claim. This is not a defect in the theorem but it is a real gap between title/abstract and result, and it should be fixed before publication.\n\nThis paper is for mathematical physicists working in the BDF programme; it deserves a serious referee. I would accept it for review, and ask for a revised abstract plus the missing initial-data hypothesis.","headline":"A valid, well-executed conditional global well-posedness theorem for 2D massless BDF dynamics; the abstract overstates the case for graphene because the physical application leans on a numerical Fermi-velocity threshold that the paper does not prove.","tokens_in":23157,"tokens_out":4330,"would_cite":true,"duration_ms":41798,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q41","81V70"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that the time-dependent Hartree-Fock dynamics of graphene has a unique solution for all times whenever the Fermi velocity exceeds a critical value.","keywords":["graphene","Hartree-Fock","Bogoliubov-Dirac-Fock","global well-posedness","massless Dirac operator","mean-field approximation","Coulomb interaction","Fermi velocity"],"falsifier":"A numerical evaluation of the function $h(v_F)$ defined in the paper, using the explicit expression for the function $g$, would determine the critical velocity $v_c=h^{-1}(2)$ to enough precision to know whether $v_F\\simeq 1.1$ satisfies the theorem's hypothesis; if $v_c$ turns out to be above $1.1$, the theorem's condition fails for graphene.","tokens_in":22256,"feed_emoji":"⚛️","tokens_out":10787,"duration_ms":99383,"temperature":0.7,"pith_summary":"This paper proves a global well-posedness theorem for the time-dependent Hartree-Fock equations that describe electrons in a graphene sheet interacting through instantaneous Coulomb forces and an external, time-dependent charge density. Such a density models local defects placed on the carbon monolayer. Because graphene's electrons behave as massless Dirac fermions with a small Fermi velocity, the effective coupling is strong and perturbative methods fail, so the authors use the non-perturbative Bogoliubov-Dirac-Fock framework, adapted from three-dimensional quantum electrodynamics to two dimensions with the Fermi sea in place of the Dirac sea. The main result says that, whenever the Fermi velocity is at least a critical value depending only on the ultraviolet cutoff, the Cauchy problem has a unique solution for all times. If correct, this means the electron state does not blow up in finite time, and questions such as scattering can be posed on a rigorous footing.","feed_headline":"Global solution proven for graphene's Hartree-Fock dynamics","feed_subtitle":"Above the critical Fermi velocity, the electron state cannot blow up, making scattering theory tractable.","key_machinery":"The carrying mechanism is the Bogoliubov-Dirac-Fock (BDF) energy $E^\\phi_{\\mathrm{BDF}}(Q)$, the Hartree-Fock energy of a state relative to the free Fermi sea, written in terms of the renormalized density matrix $Q=\\gamma-P^0_-$, together with the associated mean-field operator $D_Q=v_F\\sigma\\cdot(-i\\nabla)-\\nu*|\\cdot|^{-1}+\\rho_Q*|\\cdot|^{-1}-Q(x,y)/|x-y|$. The energy is bounded below by a coercive expression whose leading term is the kinetic energy, provided the Fermi velocity satisfies $v_F\\ge v_c=h^{-1}(2)$; that same condition, through a Hardy-type inequality, controls the exchange term by the kinetic term. Along the evolution the energy is not conserved, but its time derivative is controlled by the external charge $\\nu$, and an exponential-growth bound keeps it finite on finite time intervals; combined with coercivity this keeps the norm on the space $Y$ finite for all times. The space $Y$ encodes finite kinetic energy, finite Hilbert-Schmidt regularity of $Q|D_0|^{1/2}$, and finite Coulomb energy of the density.","core_discovery":"On the paper's own terms, the central discovery is Theorem 2.7: for a fixed ultraviolet cutoff $\\Lambda>0$, a time-dependent charge density $\\nu\\in C^1([0,+\\infty),\\mathcal{C})$, and an initial state whose renormalized density matrix $Q_I=\\gamma_I-P^0_-$ lies in the space $Y$, the Cauchy problem for $\\gamma(t)=Q(t)+P^0_-$, namely $i\\frac{d}{dt}\\gamma(t)=[D_{Q(t)},\\gamma(t)]$ with $\\gamma(0)=\\gamma_I$ and $\\gamma(t)^2=\\gamma(t)$, admits a unique global solution $\\gamma\\in C^1([0,+\\infty),P^0_-+Y)$ provided $v_F\\ge v_c=h^{-1}(2)$. Here $P^0_-$ is the projector onto the negative spectral subspace of the massless Dirac operator $D_0=-i\\sigma\\cdot\\nabla$, and $D_Q$ is the mean-field operator containing kinetic, external, direct, and exchange terms. The proof works because the Bogoliubov-Dirac-Fock energy, although not conserved in the time-dependent external field, stays bounded on finite time intervals and controls the natural norm of the state; this is exactly where the condition on the Fermi velocity enters, through an inequality that lets kinetic energy dominate the exchange term.","pith_inferences":["Editorial inference: the same scheme should transfer to other two-dimensional Dirac materials with a conical band dispersion, provided the exchange-control condition can be verified with the material's own Fermi velocity; this is a testable modification of the function $h$.","Editorial inference: the exponential energy bound is likely not sharp; under extra decay assumptions on $\\dot{\\nu}$, the method may yield polynomial or uniform bounds, which would be relevant for scattering theory.","Editorial inference: the theorem concerns instantaneous Coulomb interactions and neglects photons, spin, and lattice-scale effects; adding those would require extra terms in $D_Q$ and is a natural stress test of the global well-posedness result.","Editorial inference: the rigorous bound $v_c<2.056$ does not by itself cover graphene, so closing the gap between that bound and the numerical value $v_c\\simeq 0.36$ is a concrete numerical-analytical challenge."],"forward_implications":["For any initial state satisfying the hypotheses, the Hartree-Fock evolution is defined for all $t\\ge 0$ and is unique, so no finite-time blow-up can occur.","Time-dependent external charge densities, modeling local defects that move or change in time, are admitted as long as they are $C^1$ with values in the Coulomb space.","The ground state of the stationary problem is a stationary solution of the evolution equation, so the theorem covers perturbations of it.","Because the solution's energy and norm remain controlled, long-time questions such as scattering are well-posed and can be investigated rigorously.","If the numerical value $v_c\\simeq 0.36$ is correct, the theorem covers graphene at its physical Fermi velocity $v_F\\simeq 1.1$, as the paper argues."],"supporting_citations":[{"why":"Supplies the graphene Hartree-Fock model, the effective mean-field operator, the function $h$, and the stationary ground-state results that the dynamical theorem builds on.","marker":"[34]"},{"why":"Establishes the no-photon QED/BDF framework, the energy subtraction via thermodynamic limit, and the stability arguments adapted here.","marker":"[32]"},{"why":"Provides the trace identity used in the energy-derivative computation and the existence theory for the polarized vacuum in the BDF approximation.","marker":"[29]"},{"why":"Proves the analogous global-in-time result for the massive three-dimensional Dirac-Fock evolution whose strategy this paper follows.","marker":"[33]"},{"why":"Supplies the inequality $Q_{++}-Q_{--}\\ge Q^2$ used to relate the kinetic energy to the Schatten norms in $Y$.","marker":"[2]"},{"why":"Introduces the Bogoliubov-Dirac-Fock formalism and identifies the evolution equation as the derivative of the energy.","marker":"[13]"}],"fun_headline_variants":["Graphene's Hartree-Fock equations have global solutions","Global well-posedness for graphene's Hartree-Fock model","Fermi velocity critical for global graphene Hartree-Fock","No blow-up in graphene's Hartree-Fock evolution","Global existence proven for graphene's Hartree-Fock dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole result rests on the condition that the Fermi velocity is at least the critical value $v_F\\ge v_c=h^{-1}(2)$, which guarantees the exchange energy is dominated by the kinetic energy; the paper's rigorous estimate leaves $v_c<2.056$, so whether real graphene with $v_F\\simeq 1.1$ is covered still depends on the numerical value $v_c\\simeq 0.36$ being correct.","fun_headline_variants_meta":{"raw":{"variants":["Graphene's Hartree-Fock equations have global solutions","Global well-posedness for graphene's Hartree-Fock model","Fermi velocity critical for global graphene Hartree-Fock","No blow-up in graphene's Hartree-Fock evolution","Global existence proven for graphene's Hartree-Fock dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000557,"raw_usage":{"total_tokens":2656,"prompt_tokens":957,"completion_tokens":1699,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":573,"completion_tokens_details":{"reasoning_tokens":1614}},"tokens_in":573,"tokens_out":1699,"duration_ms":13496,"temperature":1.0,"reasoning_tokens":1614,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:36:52.153733+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A numerical evaluation of the function $h(v_F)$ defined in the paper, using the explicit expression for the function $g$, would determine the critical velocity $v_c=h^{-1}(2)$ to enough precision to know whether $v_F\\simeq 1.1$ satisfies the theorem's hypothesis; if $v_c$ turns out to be above $1.1$, the theorem's condition fails for graphene.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the graphene Hartree-Fock model, the effective mean-field operator, the function $h$, and the stationary ground-state results that the dynamical theorem builds on."},{"cited_title":"Hainzl, M","cited_arxiv_id":null,"evidence_quote":"Establishes the no-photon QED/BDF framework, the energy subtraction via thermodynamic limit, and the stability arguments adapted here."},{"cited_title":"Hainzl, M","cited_arxiv_id":null,"evidence_quote":"Provides the trace identity used in the energy-derivative computation and the existence theory for the polarized vacuum in the BDF approximation."},{"cited_title":"Hainzl, M","cited_arxiv_id":null,"evidence_quote":"Proves the analogous global-in-time result for the massive three-dimensional Dirac-Fock evolution whose strategy this paper follows."},{"cited_title":"Bach, J.-M","cited_arxiv_id":null,"evidence_quote":"Supplies the inequality $Q_{++}-Q_{--}\\ge Q^2$ used to relate the kinetic energy to the Schatten norms in $Y$."},{"cited_title":"Chaix and D","cited_arxiv_id":null,"evidence_quote":"Introduces the Bogoliubov-Dirac-Fock formalism and identifies the evolution equation as the derivative of the energy."}],"review_version":1}