{"id":"003c1b00-79f6-409a-bf12-768a03ddfa85","arxiv_id":"2411.13273","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Charge current in particle-hole symmetric graphene at the neutrality point obeys Navier-Stokes-like diffusion equations controlled by a current diffusion coefficient, not by momentum viscosity.","lead":"Using exact quantum Monte Carlo simulations and Boltzmann kinetic theory, this paper shows that electric current in clean, particle-hole symmetric graphene at half-filling can flow hydrodynamically even when the total momentum stays zero. It introduces a current diffusion coefficient that replaces viscosity and may change how graphene hydrodynamic flow experiments at the charge neutrality point are interpreted.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The QMC validation of the central claim is circular: the current-flux closure (Eq. 5) is an input to Eq. (6), so profile fits cannot verify it, and the only independent derivation is a leading-log variational calculation that overestimates C by roughly a factor of 4.5.","rationale":"The central claim requires that the exact first-moment equation (3) closes at first gradient order with a well-defined current diffusion coefficient ζ, and that the QMC simulation independently verifies the kinetic-theory prediction. The weakest point is precisely the closure (5): it is an input to the equation whose solution is fitted to the QMC profiles, so the profile agreement cannot validate the closure. The SM derivation of ζ uses the variational ansatz (39), which already enforces J_xy ∝ ∇_x j_y; it is therefore not an independent confirmation. The factor-of-4.5 discrepancy between the QMC and leading-log kinetic values of C further shows that the kinetic-theory side of the validation is only qualitative. These observations support the reader's conditional verdict: the conceptual proposal is novel and plausible, and the QMC data do show smooth, temperature-dependent profiles with an approximately 1/T scaling of ζ, but the central constitutive relation remains unverified. A direct QMC computation of J_xy(x), or a width/boundary-dependence study of the extracted ζ, would settle whether ζ is a true transport coefficient. The sign difference between Eq. (5) and SM Eq. (41) is noted as a secondary internal inconsistency that should be corrected in revision. None of these issues amount to proof that the claim is false, so the appropriate verdict remains CONDITIONAL rather than ACCEPT or REJECT.","tokens_in":20708,"tokens_out":16343,"duration_ms":184744,"concrete_test":"Compute in QMC the Kubo response of the local current-flux operator J_xy(x) to the applied y-current (i.e., the correlator between the bond-flux operator and the y-current), and test pointwise whether J_xy(x) = -ζ ∇_x j_y(x) with the same ζ obtained from the profile fit. Alternatively, or additionally, run the same temperature for at least two different sample widths or adatom densities: if the extracted ζ changes with width or boundary conditions, the first-order closure (5) is not a bulk constitutive law.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (5) asserts a Fick-type constitutive law J_xy = -ζ ∇_x j_y. This is the step that turns the exact first-moment equation (3) into the closed diffusion equation (6); without it, ζ is not defined. The paper's stated confirmation of Eq. (5) in the SM consists of evaluating both J_xy and ∇_x j_y from the variational solution of the linearized Boltzmann equation with the shear-force ansatz (39). That ansatz builds the proportionality into the trial function, so it cannot independently establish the closure. The QMC profiles are then fitted with the catenary solution (7) of Eq. (6), i.e., with the same closure assumed. Thus the agreement between QMC profiles and Eq. (7) is a consistency check, not a test of Eq. (5). If the true current flux contains higher-gradient or nonlinear terms, the ζ extracted from a catenary fit would be an effective, geometry-dependent parameter rather than a genuine transport coefficient. The large prefactor mismatch (C_QMC = 5.67 vs C_kin = 25.76) reinforces that the kinetic-theory derivation of the closure is quantitatively unreliable. As a minor internal red flag, the SM constitutive relation (Eq. 41) is written with the opposite sign (J_xy = +ζ ∇_x j_y) from Eq. (5); even if typographical, it indicates that the sign and structure of the closure are not robustly pinned down.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a hydrodynamic description of charge current in particle-hole symmetric systems at half-filling. Starting from the Boltzmann equation with a relaxation-time collision integral, the authors derive an exact first-moment equation (Eq. (3)), propose a Fick-type constitutive closure J_xy = -ζ ∇_x j_y (Eq. (5)), and obtain a Navier-Stokes-like diffusion equation for the current (Eq. (6)) whose stationary solution is the catenary profile (Eq. (7)). The paper then uses sign-problem-free auxiliary-field quantum Monte Carlo simulations of graphene strips with disordered edges to compute local DC conductivities via stochastic analytic continuation and middle-point estimators, fits the profiles to Eq. (7), and extracts a current diffusion coefficient ζ(T) ~ 1/T^{1.09±0.14} with prefactor C = 5.67±0.46. This is compared with a kinetic-theory leading-log calculation giving C = 25.76, and the qualitative agreement is presented as evidence for the Boltzmann description.","tokens_in":21076,"tokens_out":3831,"duration_ms":40975,"significance":"The central idea is interesting and potentially important. If Eq. (6) is correct, it resolves a conceptual puzzle: at charge neutrality, hydrodynamic-looking current profiles can arise even though the electric field does not generate a momentum response (SM Eq. (22)). The first-moment equation (3) is derived cleanly, and the QMC simulations are sign-problem-free, with code and data availability statements, and with useful cross-checks between SAC and middle-point estimators as well as free-case analytical tests. The new transport coefficient ζ is a plausible analogue of viscosity for current flow. However, the load-bearing closure (5) is not independently tested by the profile fits, because Eq. (7) is the solution of Eq. (6) with that closure imposed. The large prefactor mismatch between QMC and kinetic theory also raises quantitative concerns. A direct numerical test of Eq. (5) or a demonstration of width independence of the extracted ζ would substantially strengthen the claim.","major_comments":[{"comment":"The sign of the constitutive relation is inconsistent between the main text and the supplementary material: Eq. (5) states J_xy = -ζ ∇_x j_y, while SM Eq. (41) states J_xy = +ζ ∇_x j_y. Since this sign converts the exact first-moment equation (3) into Eq. (6) and determines the catenary shape (7), the correct sign must be established and the derivation made internally consistent.","section":"SM §II.3, Eq. (41)"},{"comment":"The QMC profile fits do not independently test the closure (5), because Eq. (7) is the solution of Eq. (6) with that closure imposed; the agreement between QMC profiles and Eq. (7) is therefore a consistency check rather than a confirmation of the closure. The only independent derivation in the paper is the variational Boltzmann calculation in SM §II.3, but the trial function (39) builds proportionality between J_xy and ∇_x j_y into the ansatz. The authors should test Eq. (5) directly from the QMC data (for example, compute the current-flux tensor J_xy and compare it with ∇_x j_y stripe by stripe) and/or demonstrate that the ζ extracted from catenary fits does not depend on sample width at fixed temperature.","section":"§II.1 and Eq. (6)"},{"comment":"The factor-4.5 discrepancy between C_QMC = 5.67±0.46 and C_kin = 25.76 is large for the quantitative claim that the fitted ζ is the kinetic-theory coefficient. The paper attributes this to Dirac-cone and tree-level approximations, but if the closure (5) is only approximate, the fitted value could be an effective, geometry-dependent parameter rather than a universal transport coefficient. The authors should provide a quantitative estimate of systematic uncertainties, for instance by varying the fit window, by testing width independence at fixed temperature, and by discussing the range of validity of the leading-log value.","section":"Fig. 4 and Eq. (9)"}],"minor_comments":[{"comment":"The sentence 'Later, the linear connection (7) is confirmed' should refer to Eq. (5), not Eq. (7), which is the catenary solution of Eq. (6).","section":"Main text after Eq. (5)"},{"comment":"There is a typo: 'the collision integral is treated in treated in the relaxation time approximation' should read 'the collision integral is treated in the relaxation time approximation'.","section":"SM §II.2"},{"comment":"The caption contains the typo 'errobars' instead of 'error bars'.","section":"Fig. 6 caption"},{"comment":"References [4] and [48] are the same work, as are [5] and [22]; these duplicate entries should be consolidated.","section":"References"},{"comment":"The notation σ0 = 1/4 in units of e^2/ℏ appears only in the figure caption; it should be defined in the main text where the conductivity profiles are introduced.","section":"Fig. 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal's scope and the QMC effort is substantial. My main reservation is the circularity of the central test: the profile fits use Eq. (6) whose closure (5) is not independently verified by the QMC data. That said, the kinetic-theory derivation is independent, and the data could still support the claim if the authors add a direct test of the closure or a width-scaling study. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper makes a real conceptual step: instead of trying to link momentum flow to current at the charge neutrality point, it derives and tests a Navier-Stokes-like equation directly for the charge current, with a new transport coefficient, the current diffusion coefficient ζ, replacing viscosity. The first-moment equation from the relaxation-time Boltzmann equation is clean, and the particle-hole symmetry argument for why the electric field does not pump momentum is nicely presented. The QMC work is careful: sign-problem-free simulations, SAC cross-checked against the correlator midpoint, multiple lattice sizes and adatom configurations, and openly available data. Those are serious assets.\n\nThe soft spots are real but not fatal. The constitutive relation (5), J_xy = -ζ ∇ j_y, is assumed and then the QMC profiles are fitted with the solution of the resulting equation (6). That makes the profile agreement a consistency check, not an independent confirmation of the closure. The kinetic-theory 'confirmation' in the SM uses a variational shear-force ansatz that essentially builds the proportionality into the trial function, so it does not settle the question either. The prefactor mismatch, C = 5.67 from QMC versus 25.76 from kinetic theory, is a quantitative gap of a factor 4.5; the paper attributes it to Dirac dispersion and tree-level diagrams, which is plausible, but it means the coefficient is not pinned down. There is also a small sign inconsistency between Eq. (5) and SM Eq. (41) for J_xy; likely a typo, but it should be fixed because the sign carries physical content.\n\nHaving said that, the central claim deserves a serious referee. The observation that current profiles at half-filling can be described without momentum flow is novel and addresses a genuine puzzle in graphene hydrodynamics. The temperature dependence ζ ~ 1/T is robust across estimators, and the framework is generalizable beyond graphene. The natural next step is a direct QMC measurement of J_xy to test the closure explicitly; the authors should be pushed to do that or to state plainly that the closure is assumed.\n\nThe paper is written for condensed-matter theorists and computational transport people. It is not a finished quantitative theory, but it is a solid, interesting contribution that should go to peer review: the math is coherent, the numerics are reproducible, and the central idea is worth engaging. I would recommend accepting the paper for review with a request for revision: fix the sign typo, add a direct discussion of the closure's testability, and either improve the prefactor comparison or explain why the 4.5 mismatch is expected.","headline":"A genuinely new current-diffusion framework for charge-neutral graphene hydrodynamics, with clean QMC profiles, but the fit-based validation cannot independently prove the central closure and the kinetic-theory prefactor is off by a factor of 4.5.","tokens_in":21573,"tokens_out":5021,"would_cite":true,"duration_ms":62301,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Charge current in particle-hole symmetric systems obeys its own Navier-Stokes equation, with a current diffusion coefficient in place of viscosity, even though the electric field creates no momentum.","keywords":["charge neutrality","particle-hole symmetry","graphene","hydrodynamic transport","current diffusion coefficient","Boltzmann transport theory","quantum Monte Carlo","Navier-Stokes equation"],"falsifier":"Compute the off-diagonal current flux $J_{xy}$ and the current gradient $\\nabla_x j_y$ independently from the QMC or kinetic-theory distribution functions and check whether their ratio is a position-independent constant at fixed temperature and field; if $J_{xy}/\\nabla_x j_y$ varies across the strip or with the driving field strength, the linear closure fails and $\\zeta$ is not a well-defined transport coefficient.","tokens_in":20492,"feed_emoji":"🌊","tokens_out":9539,"duration_ms":93143,"temperature":0.7,"pith_summary":"At the charge-neutrality point of a particle-hole symmetric system, an electric field produces charge current but no net momentum, so the usual Navier-Stokes picture in terms of momentum flow cannot directly explain measured current profiles. This paper claims that hydrodynamics still emerges for the charge current itself: taking the first moment of the two-fluid electron-hole Boltzmann equation and closing it with a linear relation between current flux and current gradient yields a Navier-Stokes-type diffusion equation for the current, with a new coefficient, the current diffusion coefficient $\\zeta$, replacing viscosity. The predicted catenary current profiles are reproduced by numerically exact quantum Monte Carlo simulations of clean graphene strips with disordered edges, and the extracted $\\zeta(T)$ follows roughly $1/T$, in qualitative agreement with a leading-log Boltzmann calculation. If correct, experimental current-imaging data at half-filling can be interpreted without any momentum-current coupling, and the same description should apply to any particle-hole symmetric system.","feed_headline":"Hydrodynamic flow emerges even with zero momentum","feed_subtitle":"Quantum Monte Carlo shows charge current obeys its own diffusion equation, with a new coefficient replacing viscosity.","key_machinery":"The central object is the first-moment equation for the charge current obtained from the Boltzmann equation in the electron-hole basis: $\\partial_t j_n + \\partial_m J^{mn} + eE_m E^{mn} = -j_n/(\\tau_{ee}/2)$, where $J^{mn}$ is the current flux tensor (the current analogue of the stress tensor) and $E^{mn}$ couples to the external field and stays nonzero at charge neutrality, $E^{mn}=-(2Te\\ln 2/\\pi)\\delta^{mn}$. The argument is carried by the phenomenological closure $J^{xy}=-\\zeta\\nabla_x j_y$, which converts the exact first-moment equation into the closed diffusion equation whose solution is the catenary profile used to fit the QMC data, and by a leading-log variational solution of the full linearized Boltzmann equation that independently yields $\\zeta = C(\\epsilon a t)^2 t/T$.","core_discovery":"The paper claims that at charge neutrality in a particle-hole symmetric system, the electric field induces a charge current while leaving momentum strictly zero, and that this current obeys its own hydrodynamic equation. Combining the two-fluid (electrons and holes) Boltzmann equation with the linear constitutive relation $J^{xy}=-\\zeta\\nabla_x j_y$, the stationary current satisfies $-\\zeta\\nabla_x^2 j_y(x)+eE_y E^{yy}=-j_y/(\\tau_{ee}/2)$, a Navier-Stokes-type diffusion equation in which $\\zeta$ plays the role viscosity plays for momentum. Its solution for a strip is a catenary profile, and the paper shows that numerically exact QMC conductivity profiles in clean graphene strips follow this shape. From the fits it extracts $\\zeta(T)$, finding the same power law as the leading-log variational Boltzmann prediction, $\\zeta\\sim T^{-1.09\\pm0.14}$ with $C=5.67\\pm0.46$ versus the theoretical $C=25.76$, a difference attributed to the Dirac-cone approximation and to restricting the collision integral to tree-level diagrams.","pith_inferences":["Editorial inference: if the current-diffusion picture is correct, previously reported \"viscosity\" values extracted from current-imaging experiments at the charge-neutrality point should be reinterpreted as the current diffusion coefficient, since the standard linear-response connection between electric field and momentum vanishes by particle-hole symmetry.","Editorial inference: the same two-fluid closure could be applied to other conserved or approximately conserved channels in particle-hole symmetric systems, such as spin or valley currents, each with its own analogue of $\\zeta$ replacing the corresponding viscosity.","Editorial inference: a direct numerical test of the closure $J^{xy}=-\\zeta\\nabla_x j_y$, computing both sides from the QMC distribution functions rather than from a fit to the catenary profile, would independently verify that $\\zeta$ is a well-defined transport coefficient and reveal whether gradient corrections become important.","Editorial inference: the factor of roughly 4.5 between the Boltzmann prefactor ($C=25.76$) and the QMC prefactor ($C=5.67$) gives a quantitative target for improved collision integrals; a next-order calculation that closes this gap would strongly validate the two-fluid current-diffusion description."],"forward_implications":["In any particle-hole symmetric system at half-filling, hydrodynamic current profiles can be described directly by the current-diffusion equation; the paper explicitly names AB-stacked bilayer graphene with quadratic band touching and narrow-gap semiconductors with symmetric bands as candidates.","Experimental current-imaging data at the graphene charge-neutrality point can be analyzed with the current diffusion coefficient instead of shear viscosity, removing the need for an artificial chemical potential $\\mu=k_B T$ to couple momentum to current.","Within the hydrodynamic window, $\\zeta$ grows roughly as $1/T$ toward lower temperatures, with a prefactor set by the Coulomb interaction strength, so the temperature dependence of the extracted profiles is a testable prediction of the theory.","Quantum Monte Carlo can serve as a controlled substitute for experiments in transport: because all scattering mechanisms are fixed by a microscopic Hamiltonian, it can test and refine Boltzmann collision integrals, including corrections beyond the tree-level approximation."],"supporting_citations":[{"why":"Supplies the leading-log variational calculation of shear viscosity at charge neutrality on which the current-diffusion calculation is directly patterned.","marker":"[4]"},{"why":"Provides the current-imaging experiment at the charge-neutrality point whose interpretation motivated replacing viscosity with the current diffusion coefficient.","marker":"[12]"},{"why":"Provides the finite-temperature auxiliary-field quantum Monte Carlo method used to obtain exact current-current correlators for the microscopic model.","marker":"[18]"},{"why":"Derives the relaxation-time Boltzmann treatment of the homogeneous current equation whose stationary limit gives the critical conductivity.","marker":"[25]"},{"why":"Supplies the Coulomb collision integral and the frequency-temperature criteria defining the hydrodynamic regime used here.","marker":"[26]"},{"why":"Provides the variational method for solving the linearized Boltzmann equation used to extract the current diffusion coefficient.","marker":"[27]"},{"why":"Supplies the long-range Coulomb interaction matrix for free-standing graphene entering the microscopic Hamiltonian.","marker":"[30]"},{"why":"Provides the lattice-fermion quantum Monte Carlo implementation used to compute the current-current correlator and conductivity profiles.","marker":"[36]"}],"fun_headline_variants":["Charge current flows hydrodynamically even with zero momentum","New diffusion coefficient replaces viscosity at charge neutrality","Hydrodynamic current emerges in particle-hole symmetric systems","Zero momentum, but charge flow behaves like a fluid","QMC reveals current's own Navier-Stokes equation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole extraction of $\\zeta$ rests on the linear closure $J^{xy}=-\\zeta\\nabla_x j_y$, which turns the exact first-moment equation into the closed diffusion equation; the paper itself notes that this relation is confirmed only within the leading-log variational solution of the Boltzmann equation, so the QMC profile fits, which already assume the solution of that closed equation, cannot independently test the closure.","fun_headline_variants_meta":{"raw":{"variants":["Charge current flows hydrodynamically even with zero momentum","New diffusion coefficient replaces viscosity at charge neutrality","Hydrodynamic current emerges in particle-hole symmetric systems","Zero momentum, but charge flow behaves like a fluid","QMC reveals current's own Navier-Stokes equation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000417,"raw_usage":{"total_tokens":2167,"prompt_tokens":980,"completion_tokens":1187,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":1112}},"tokens_in":596,"tokens_out":1187,"duration_ms":10258,"temperature":1.0,"reasoning_tokens":1112,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:38:19.517653+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the off-diagonal current flux $J_{xy}$ and the current gradient $\\nabla_x j_y$ independently from the QMC or kinetic-theory distribution functions and check whether their ratio is a position-independent constant at fixed temperature and field; if $J_{xy}/\\nabla_x j_y$ varies across the strip or with the driving field strength, the linear closure fails and $\\zeta$ is not a well-defined transport coefficient.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the current-imaging experiment at the charge-neutrality point whose interpretation motivated replacing viscosity with the current diffusion coefficient."},{"cited_title":"Sorella and E","cited_arxiv_id":null,"evidence_quote":"Provides the finite-temperature auxiliary-field quantum Monte Carlo method used to obtain exact current-current correlators for the microscopic model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the relaxation-time Boltzmann treatment of the homogeneous current equation whose stationary limit gives the critical conductivity."},{"cited_title":"Pongsangangan, T","cited_arxiv_id":null,"evidence_quote":"Provides the variational method for solving the linearized Boltzmann equation used to extract the current diffusion coefficient."},{"cited_title":"Pongsangangan, P","cited_arxiv_id":null,"evidence_quote":"Supplies the long-range Coulomb interaction matrix for free-standing graphene entering the microscopic Hamiltonian."}],"review_version":1}