{"id":"10167f89-c163-428e-8e6a-e25646d85c0b","arxiv_id":"2608.11184","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Starting from a kinetic theory of hard-core active Brownian particles, the authors derive the deterministic part of Active Model B+ and give explicit coefficients, using a fitted pair-correlation function to obtain phase separation.","lead":"This paper derives Active Model B+, a standard equation for how self-propelled particle densities change, from a microscopic kinetic theory of hard active Brownian particles. It gives explicit formulas for the model's coefficients, though one ingredient, the pair correlation function, is chosen by hand to make the theory work.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unproved parity cancellation in the fourth-order Chapman–Enskog step (Supplement Eqs. A60–A63) is the most load-bearing gap: if false, the density equation may not reduce to AMB+.","rationale":"The reader's weakest assumption focused on the ad hoc parametrization chi_eff(ρ), which indeed controls the quantitative coefficients and the signs of λ and ζ. That is a legitimate concern, and it supports a CONDITIONAL verdict. However, the single most load-bearing issue is the unproved parity cancellation in the Supplement, because it threatens the formal derivation of the AMB+ structure itself, not merely the numerical values of the coefficients. If the cancellation is false, the derived equation may contain extra gradient terms that cannot be absorbed into AMB+, undermining the paper's central claim. Since the authors explicitly acknowledge the missing proof, this is not a manufactured objection; it is an identified gap in the manuscript. The verdict remains CONDITIONAL: the paper should be accepted only if the cancellations are proven or the derivation adjusted. The reader's verdict is unchanged because it already reflects a significant degree of conditionality, but the specific load-bearing location differs from the chi_eff closure. The concrete test is feasible and would settle the matter directly, without requiring simulations or external input.","tokens_in":18158,"tokens_out":4675,"duration_ms":41543,"concrete_test":"Independently verify the two cancellations by evaluating the uncontracted integrals with the explicit collision displacement Δ_1 from Eq. (6) and Fig. 1. Specifically, compute I(1) = ∫ Δ_α n_{2β} d n1 dH and I(2) = ∫ Δ_α n_{1β} d n1 dH (with dH defined in Eq. A7), and check whether I(1)+I(2) = 0 identically for all α,β. Similarly verify the pair in Eqs. (A62)–(A63). This can be done analytically by symbolic integration over n1, n2, σ with the Heaviside step, or numerically to high precision. If the sums are nonzero, the fourth-order equation acquires extra terms and the AMB+ mapping fails; if they vanish exactly for all α,β, the structural derivation is restored.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of the fourth-order density equation, which is the step that produces the AMB+ structure, relies on two cancellations: I^1_03 + I^1_30 = 0 and I^1_12 + I^1_21 = 0. These are stated in the Supplemental Material immediately after Eqs. (A60)–(A63) and explicitly left unproved: the authors write that the cancellation 'comes from a parity symmetry... which can probably be proven using the explicit expression for the displacement Δ. The proof falls out of the scope of this article so we leave it for any interested reader.' The final AMB+ form (Eqs. 1–2, with the specific gradient structure of the non-equilibrium fluxes) depends on these cancellations; if they fail, additional terms of order ∇^2(χρ(...)) and similar would survive in Eq. (13), producing terms that are not representable as −∇·(λ(∇φ)^2∇φ + ζ(∇^2φ)∇φ) plus the free-energy flux. This is a correctness risk at the heart of the claim that the kinetic theory 'derives AMB+ from first principles.' Unlike the chi_eff parametrization, which the authors openly describe as a proposed closure, this cancellation is presented as a mathematical fact without proof. The paper itself flags it as missing support, so the derivation is incomplete as written. No amount of fitting the pair correlation function can repair this, because the cancellation is needed for the form of the equation regardless of chi.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a Chapman-Enskog expansion of a kinetic equation for hard-core active Brownian particles (ABP) in the high-persistence regime, aiming to derive a fourth-order density field equation with the structure of deterministic Active Model B+ (AMB+). The authors compute explicit expressions for the transport coefficients kappa, xi, lambda, and zeta in terms of the microscopic parameters (V, sigma, Peclet number) to leading order in Pe, using an effective parametrization chi_eff(rho) of the pair correlation function to ensure a double-well free energy. They also provide the resulting free energy and Ginzburg coefficient. The derivation is presented in the main text with extensive algebraic details relegated to the Supplemental Material.","tokens_in":18548,"tokens_out":5095,"duration_ms":44952,"significance":"If the central derivation is correct, this would be an important bottom-up connection between ABP kinetic theory and the phenomenological AMB+ model, giving explicit microscopic expressions for the nonequilibrium transport coefficients and a falsifiable prediction (opposite signs of lambda and zeta, implying coarsening and complete phase separation). The paper is unusually transparent about its limitations: it explicitly flags the ad hoc nature of chi_eff, the missing proof of the parity cancellations, and the discrepancy between the negative kappa and the known short-wavelength stability of the kinetic theory. The detailed Supplemental Material provides many of the intermediate integrals, which is a strength. However, the load-bearing gaps (unproved cancellations and the fitted closure) mean that the 'first-principles derivation' claim is not fully supported as written.","major_comments":[{"comment":"The derivation of the fourth-order density equation relies on the equalities I^1_03 + I^1_30 = 0 and I^1_12 + I^1_21 = 0, stated immediately after these equations. The authors explicitly defer the proof, writing that the cancellation 'comes from a parity symmetry... which can probably be proven using the explicit expression for the displacement Δ' and 'falls out of the scope of this article.' This is a load-bearing step: if either cancellation fails, additional terms of order ∇·(χρ ...) survive in Eq. (A59), and the density equation would not reduce to the AMB+ form of Eqs. (1)-(2). Since the central claim of the paper is the derivation of AMB+ from first principles, this omitted proof leaves the conclusion incomplete. I request a rigorous proof of these identities (or a numerical verification using the explicit Δ for the model parameters used), or a reformulation of the derivation that does not require them.","section":"Supplemental Material, after Eqs. (A60)-(A63)"},{"comment":"The claim of a first-principles derivation is weakened by the effective pair correlation function chi_eff(rho) = 1 + A rho - B rho^2 + C rho^3 with A = 3.14 sigma^2, B = 3.7 sigma^4, C = 3.1 sigma^6. As the authors state, these coefficients are chosen so that the spinodal density is close to simulation values and so that D''(rho_s) > 0; the explicit expressions for kappa, xi, lambda, and zeta in Eqs. (17)-(20) are therefore outputs of this fitted closure, not of the kinetic theory alone. While the paper is transparent about this, the abstract's phrase 'derive AMB+ from first principles' overstates the result. The authors should either derive chi_eff from a controlled microscopic approximation or clearly qualify the result as a closure-based derivation whose quantitative predictions depend on the choice of chi.","section":"Main text, paragraph after Eq. (16)"},{"comment":"The paper acknowledges that the negative kappa in Eq. (17) would imply infinitely rough states beyond the spinodal, in contrast to the numerical stability of small-wavelength modes found in Ref. [37], attributing this to truncation at fourth order. This is not a minor caveat: it indicates that the derived AMB+ equation does not reproduce the linear stability of the underlying kinetic theory in the regime where it is intended to apply. The authors should provide evidence that higher-order terms indeed restore stability, for example by estimating the next-order contribution to the linearized dispersion relation or by identifying the range of validity of the fourth-order gradient expansion. Without such a check, the predictive content of the derived coefficients, including the predicted sign structure of lambda and zeta, remains to be established.","section":"Main text, paragraph after Eq. (20)"}],"minor_comments":[{"comment":"The text reads 'Ginzbug coefficient'; this should be 'Ginzburg coefficient.'","section":"Eq. (15)"},{"comment":"The sentence 'needing for one colliding particle to to be at r1 - Delta1 and r1, respectively' contains a duplicated 'to' and is ungrammatical; please rephrase.","section":"Main text, paragraph after Eq. (6)"},{"comment":"The accent in 'P\\'eclet' appears as an escaped apostrophe; please use the proper Unicode character (Péclet) to avoid typesetting issues.","section":"Abstract"},{"comment":"The notation for the integration measure dH = V sigma |(n2 - n1)·sigma| Theta(...) dn2 dsigma is introduced but the hat notation for unit vectors is not applied consistently in the surrounding text; please standardize.","section":"Supplemental Material, Eq. (A7)"},{"comment":"The coefficient values for the three closures are quoted with differing numbers of significant digits; a short statement about numerical precision would improve reproducibility.","section":"End Matter, Eqs. (21)-(28)"}],"recommendation":"major_revision","confidential_remarks":"The paper is technically ambitious and unusually candid about its limitations, which I appreciate. The most serious issue is the unproved parity cancellation in the fourth-order Chapman-Enskog step; without a proof or numerical check, the derivation of the AMB+ structure is incomplete. The fitted chi_eff is a significant caveat to the 'first-principles' claim, but it is openly disclosed and could be addressed by reframing the claims. The negative-kappa discrepancy with the kinetic-theory stability is also concerning, though the authors' truncation-based explanation is plausible. I do not recommend rejection because these issues are fixable within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news: this is the first derivation I know of that takes a kinetic theory for hard-core ABP all the way to fourth order in gradients and lands on the deterministic AMB+ structure with explicit formulas for kappa, xi, lambda and zeta in terms of V, sigma and Pe. The diffusive-order work stopped short, the quorum-sensing derivations miss zeta, and the DDFT route needs soft potentials. So the structural achievement is genuine and the coefficient formulas are a new result, with the caveat that the numbers in Eqs. (17)-(20) come from a hand-fitted closure.\n\nThe derivation is mostly transparent. The Chapman-Enskog hierarchy is laid out carefully, the integral reduction is systematic, and the authors are honest about what they did and did not do. They say plainly that chi_eff is an effective parametrization, not the real correlation function, and they acknowledge that the other closures fail to give two minima. I also believe them when they say the negative kappa is a truncation artifact; that kind of admission usually signals the authors know where the error bars are.\n\nThe soft spots are real but not equally soft. The hand-fitted chi_eff is the weaker of the two in my view, because it is doing load-bearing work: without it there is no two-minimum free energy and hence no phase separation. The authors are candid about this, which does not remove the arbitrariness but does mean the paper is not hiding anything. The bigger formal gap is the claimed parity cancellation in the Supplemental Material (Eqs. A60-A63). They state that I^1_03 + I^1_30 = 0 and I^1_12 + I^1_21 = 0, and they explicitly leave the proof to the reader. Those cancellations are needed to get the AMB+ form, so as written the derivation is incomplete at the exact step that justifies the headline claim. That is a genuine referee-level concern, not a nitpick.\n\nI am less worried than the stress-test note about whether this kills the paper. The cancellation is plausible from the symmetry they describe, and the structure of the remaining terms is consistent with what a fourth-order gradient expansion should produce. But 'plausible' is not 'proven', and the paper itself flags the gap. A serious referee should ask for the proof or an independent check of those two equalities. The chi_eff fitting is less fixable, since any explicit coefficient set will inherit its choice, but the authors have already framed the paper as a route and a framework rather than a final quantitative prediction, and I think that framing is honest.\n\nBottom line: this deserves a serious referee. The derivation is a substantive step, the limitations are stated, and the one mathematical gap is localized and identifiable. I would send it out, with a recommendation that the parity cancellations be proven or verified numerically before publication.","headline":"A genuine fourth-order Chapman-Enskog derivation of deterministic AMB+ for hard-core ABP, worth refereeing despite an unproved parity cancellation and a hand-fitted pair correlation function.","tokens_in":19093,"tokens_out":1397,"would_cite":false,"duration_ms":12036,"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":"This paper derives Active Model B+ for hard-core active Brownian particles from a kinetic equation, with all transport coefficients explicit in the microscopic parameters.","keywords":["active Brownian particles","Active Model B+","motility-induced phase separation","kinetic theory","Chapman-Enskog expansion","pair correlation function","Péclet number","hard disks"],"falsifier":"Simulate hard-disk active Brownian particles at $\\mathrm{Pe}\\gg1$, measure the pair correlation at contact $g(\\sigma;\\rho)$ across densities, and check whether it makes $D''(\\rho_s)>0$ with spinodal near $\\rho_s\\sigma^2\\approx0.32$; if it does not, the two-minimum free energy and the coefficient values in Eqs. (17)–(20) do not follow. A complementary check is to measure coarsening or interfacial tension and test the predicted $\\lambda>0$, $\\zeta<0$ sign pattern.","tokens_in":17883,"feed_emoji":"🔬","tokens_out":13744,"duration_ms":103276,"temperature":0.7,"pith_summary":"This paper seeks a first-principles route from the microscopic dynamics of hard-core active Brownian particles (ABP) to Active Model B+ (AMB+), the standard field theory used to describe active phase separation. Starting from an Enskog-like kinetic equation valid in the high-persistence regime, the authors carry out a Chapman–Enskog expansion to fourth order in gradients and obtain a deterministic equation of the AMB+ form for the density. All four transport coefficients—the stiffness $\\kappa$, the density-dependent stiffness $\\xi$, and the two time-reversal-breaking fluxes $\\lambda$ and $\\zeta$—are given explicitly in terms of particle speed $V$, diameter $\\sigma$, Péclet number $\\mathrm{Pe}$, and a proposed effective pair-correlation function. Because $\\lambda$ and $\\zeta$ come out with opposite signs, the AMB+ phase diagram predicts coarsening and complete phase separation for ABP. The derivation matters because it connects qualitative top-down phenomenology to tunable microscopic parameters.","feed_headline":"Kinetic theory yields AMB+ equations for active Brownian particles","feed_subtitle":"All four transport coefficients follow from particle speed, size, and Péclet number; the signs predict full phase separation.","key_machinery":"The load-bearing mechanism is the Enskog-like collision operator (Eq. 6) in which a collision between two high-persistence ABP displaces each particle by an effective vector $\\boldsymbol{\\Delta}_i$ without changing its director; this operator is expanded in gradients to fourth order. The diffusion-order mobility $D(\\rho)=D_0[1-(\\rho^2\\chi)'/(2\\rho^*)]$ sets the spinodal through $D(\\rho_s)=0$, and the Chapman–Enskog hierarchy is solved with the normalization $f^{(0)}=\\rho/2\\pi$ and isotropic-tensor reduction of the angular integrals at each order. To close the theory at the nonlinear level the paper proposes the effective pair-correlation parametrization $\\chi_{\\mathrm{eff}}(\\rho)=1+A\\rho-B\\rho^2+C\\rho^3$ with $A=3.14\\sigma^2$, $B=3.7\\sigma^4$, $C=3.1\\sigma^6$; this object is what makes $D''(\\rho_s)>0$ and gives the two-minimum free energy, and the numerical values of the AMB+ coefficients in Eqs. (17)–(20) are outputs of that choice.","core_discovery":"The central claim is that the deterministic part of AMB+ is not merely a phenomenological fit but follows from the kinetic theory of hard-core ABP when interactions are treated as instantaneous effective displacements at infinite persistence. Working at leading order in $\\mathrm{Pe}\\gg 1$, the authors solve the kinetic equation by Chapman–Enskog, close it at fourth order in gradients, and linearize around the spinodal density $\\rho_s$ fixed by $D(\\rho_s)=0$. The resulting equation has the AMB+ structure $\\partial_t\\phi=-\\nabla\\cdot\\mathbf{J}$ with $\\mathbf{J}=-\\nabla(\\mu_{\\mathrm{eq}}+\\lambda(\\nabla\\phi)^2)+\\zeta(\\nabla^2\\phi)\\nabla\\phi$, and the coefficients take the explicit forms $\\kappa$, $\\xi$, $\\lambda$, $\\zeta$ in terms of $V$, $\\sigma$, and $\\mathrm{Pe}$ (Eqs. 17–20). To obtain a free energy with two minima the paper introduces the effective pair correlation $\\chi_{\\mathrm{eff}}(\\rho)=1+3.14\\sigma^2\\rho-3.7\\sigma^4\\rho^2+3.1\\sigma^6\\rho^3$, which satisfies $\\chi(0)=1$, is monotone, puts the spinodal near $\\rho_s\\sigma^2\\approx0.32$, and makes $D''(\\rho_s)>0$; with this choice $\\lambda>0$ and $\\zeta<0$, which the AMB+ phase diagram reads as coarsening and complete phase separation. The paper also checks that the hard-disk closure gives coefficients of similar magnitude near the spinodal.","pith_inferences":["A direct simulation test of the sign pattern exists: measure the coarsening exponent or interfacial tension of hard-disk ABP at $\\mathrm{Pe}\\gg1$ and check that it follows the complete-separation branch of the AMB+ phase diagram rather than the bubbly branch.","Measuring $g(\\sigma;\\rho)$ in ABP simulations and substituting the measured pair correlation would make the prediction parameter-free; agreement or disagreement with the proposed cubic $\\chi_{\\mathrm{eff}}$ would settle whether the chosen parametrization is the right closure.","Because the kinetic theory contains no noise, the stochastic part of AMB+ is not derived; introducing a fluctuating collision operator is the natural next step and could determine whether microphase-separated states in AMB+ require the noise terms.","The noncommutation of adiabatic elimination and gradient expansion suggests that lower-order or dynamical-density-functional closures that eliminate the polarization first will systematically miss terms of the $\\zeta$ type, which appear only at fourth order in gradients."],"forward_implications":["The four AMB+ coefficients are no longer free: $\\kappa$, $\\xi$, $\\lambda$, and $\\zeta$ are fixed by $V$, $\\sigma$, $\\mathrm{Pe}$, and the chosen $\\chi$, so the coarse-grained model is fully determined from the microscopic parameters.","With $\\lambda>0$ and $\\zeta<0$, the AMB+ parameter map places hard-core ABP in the regime of coarsening and complete phase separation, not microphase separation or reversed Ostwald ripening.","Because the leading terms in $\\kappa$ and $\\zeta$ change sign across the spinodal through their $\\epsilon\\mathrm{Pe}^3$ contributions, the character of the gradient terms reorganizes exactly at the spinodal density.","The negative $\\kappa$ would produce infinitely rough short-wavelength states, and the paper attributes this to truncating the gradient expansion at fourth order, with higher-order terms expected to restore stability.","Using the hard-disk $\\chi_{\\mathrm{hd}}$ instead of the effective $\\chi$ gives nearly the same coefficient magnitudes near the spinodal, indicating the leading coefficients depend mainly on the pair correlation in the spinodal region."],"supporting_citations":[{"why":"Supplies the effective-displacement kinetic theory and the high-persistence spinodal density that the derivation extends to fourth order.","marker":"[37]"},{"why":"Defines AMB+ and its $\\lambda$–$\\zeta$ phase diagram, used to translate the computed signs into coarsening and complete phase separation.","marker":"[22]"},{"why":"The authors' earlier diffusive-order hydrodynamic equations from the same kinetic theory, which the fourth-order expansion here corrects by keeping gradient and adiabatic-elimination orders together.","marker":"[40]"},{"why":"Provides the equilibrium hard-disk pair correlation $\\chi_{\\mathrm{hd}}$ used as an alternative closure to test the stability of the coefficient values.","marker":"[42]"},{"why":"Simulation phase diagram of active Brownian disks that sets the target spinodal density $\\rho_s\\sigma^2\\approx0.32$ used to choose the effective $\\chi$.","marker":"[44]"},{"why":"The Chapman-Enskog method, the gradient-expansion machinery on which the entire derivation relies.","marker":"[38]"}],"fun_headline_variants":["Kinetic theory pins down all AMB+ coefficients from first principles","ABP kinetic theory yields full AMB+ with explicit transport laws","First-principles AMB+ from active Brownian kinetics","All AMB+ coefficients derived from particle parameters","Kinetic theory gives the full AMB+ transport coefficients"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative content rests on the proposed polynomial pair-correlation function $\\chi_{\\mathrm{eff}}(\\rho)=1+3.14\\sigma^2\\rho-3.7\\sigma^4\\rho^2+3.1\\sigma^6\\rho^3$, whose coefficients are chosen by hand to satisfy four reasonable conditions; the standard Boltzmann and hard-disk closures fail to give the required second derivative, and the paper does not claim $\\chi_{\\mathrm{eff}}$ is the true ABP correlation function.","fun_headline_variants_meta":{"raw":{"variants":["Kinetic theory pins down all AMB+ coefficients from first principles","ABP kinetic theory yields full AMB+ with explicit transport laws","First-principles AMB+ from active Brownian kinetics","All AMB+ coefficients derived from particle parameters","Kinetic theory gives the full AMB+ transport coefficients"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000502,"raw_usage":{"total_tokens":2461,"prompt_tokens":962,"completion_tokens":1499,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":578,"completion_tokens_details":{"reasoning_tokens":1416}},"tokens_in":578,"tokens_out":1499,"duration_ms":10086,"temperature":1.0,"reasoning_tokens":1416,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:29:16.217605+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate hard-disk active Brownian particles at $\\mathrm{Pe}\\gg1$, measure the pair correlation at contact $g(\\sigma;\\rho)$ across densities, and check whether it makes $D''(\\rho_s)>0$ with spinodal near $\\rho_s\\sigma^2\\approx0.32$; if it does not, the two-minimum free energy and the coefficient values in Eqs. (17)–(20) do not follow. A complementary check is to measure coarsening or interfacial tension and test the predicted $\\lambda>0$, $\\zeta<0$ sign pattern.","supporting_citations":[{"cited_title":"Pinto-Goldberg and R","cited_arxiv_id":null,"evidence_quote":"The authors' earlier diffusive-order hydrodynamic equations from the same kinetic theory, which the fourth-order expansion here corrects by keeping gradient and adiabatic-elimination orders together."},{"cited_title":"Henderson, A simple equation of state for hard discs, Molecular Physics30, 971 (1975)","cited_arxiv_id":null,"evidence_quote":"Provides the equilibrium hard-disk pair correlation $\\chi_{\\mathrm{hd}}$ used as an alternative closure to test the stability of the coefficient values."},{"cited_title":"Digregorio, D","cited_arxiv_id":null,"evidence_quote":"Simulation phase diagram of active Brownian disks that sets the target spinodal density $\\rho_s\\sigma^2\\approx0.32$ used to choose the effective $\\chi$."},{"cited_title":"Chapman and T","cited_arxiv_id":null,"evidence_quote":"The Chapman-Enskog method, the gradient-expansion machinery on which the entire derivation relies."}],"review_version":1}