{"id":"ce7c9db6-6187-41ac-9b4c-afcd70630c21","arxiv_id":"2411.14536","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":25,"one_line_summary":"A nonequilibrium coexistence theory with a crystallinity order parameter produces a phase diagram for active Brownian spheres that matches simulated solid-fluid and liquid-gas binodals.","lead":"This paper builds a statistical mechanics theory of crystallization for active Brownian spheres, combining derived field dynamics with simulation-guided equations of state. It predicts a phase diagram with solid-fluid and liquid-gas coexistence curves and claims quantitative agreement with simulations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The crystallinity-flux neglect in SM Eqs. (S11)-(S12) is the load-bearing assumption; it is explicitly admitted but never tested, and if it fails the generalized Maxwell construction loses its foundation.","rationale":"The reader's weakest_assumption identifies exactly the condition I find most load-bearing: the neglect of crystallinity flux, |∂z J_ψ| << |s_ψ|, which converts the crystallinity dynamics into Model A form and licenses the coexistence criteria of Ref. [31]. The paper itself flags this approximation in SM Eqs. (S11)-(S12), and the main text only expresses an expectation that it holds. Because the entire nonequilibrium Maxwell construction, and hence the central phase diagram, depends on this inequality, it is the correct target for a decisive test. I do not think the concern should overturn the verdict at this stage: the comparison with simulation is encouraging, and the authors are transparent about the approximate nature of their criteria. But a direct numerical measurement of Jψ in a coexisting solid-fluid interface is a concrete, feasible check that would either validate the Model A reduction or show that the theory is not yet grounded. Other issues, such as the many fitted EOS constants and the failure to recover the passive hard-sphere limit at low activity, are real but secondary; they would matter less if the interfacial flux were shown to be negligible. Since the reader already flagged this assumption and recommended a conditional verdict, my read does not change the verdict.","tokens_in":19903,"tokens_out":3362,"duration_ms":35777,"concrete_test":"Using the same Brownian dynamics model and system sizes (HOOMD-Blue, N ≥ 55296), prepare a slab with coexisting solid and fluid at the predicted binodal densities for ℓ0/D ≈ 16.9 and 22.3. After reaching steady state, compute Jψ from the microscopic definition in SM Eq. (S12) and the crystallinity generation term sC from the same trajectories, binned in z, and evaluate R = ∫ |∂z Jψ| dz / ∫ |sψ| dz across the interface. If R is not much smaller than 1 (e.g., R > 0.1), the Model A reduction fails and the generalized Maxwell construction is unsupported. As a complementary check, solve the full steady-state equation −∂z Jψ + sC − (ℓ0^2/24) ∂z[U ∂z(U sC)] = 0 with the measured Jψ and compare the resulting ψV(z) profile with the uψ = 0 branch used in the paper.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is that Eqs. (4a)-(4b) reproduce both coexistence binodals and the triple point. Those criteria require the nonconserved crystallinity field to obey Model A dynamics, i.e. |∂z J_ψ| << |s_ψ|, so that the steady-state condition fψ = 0 becomes the bulk pseudopotential condition u_bulk_ψ(ρ,ψ) = 0. The Supplemental Material states this directly (SM Eqs. S11-S12): 'While the flux is clearly not identically zero, we nevertheless approximate it as negligible in comparison to the generation terms.' The main text only says the authors expect this condition to be approximately met for active crystallization; no measurement of Jψ in a coexistence interface, and no scaling estimate, is provided. Since Eq. (5b) and SM Eq. (S30) define uψ entirely from sC and its gradients, any non-negligible ∂z Jψ changes the force balance that fixes ψ*_V(ρ), and the approximate generalized Maxwell construction Eq. (4) no longer follows. This is especially consequential because the authors note that exact nonequilibrium coexistence criteria generally do not exist for ABPs and that Eq. (4) is only well-defined in the high-activity limit, yet it is applied across the full activity range. A direct interfacial measurement of Jψ and sψ would resolve whether the phase diagram is built on a controlled approximation or on an untested cancellation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a statistical-mechanical field theory for three-dimensional active Brownian spheres, describing the system by a conserved density field and a nonconserved crystallinity order parameter. It derives approximate pseudopotentials u_ρ and u_ψ from the microscopic equations of motion using an Irving–Kirkwood procedure, introduces equations of state for the pressure, the renormalized active speed, the preferred crystallinity, and the order-disorder volume fraction (with coefficients fitted to Brownian dynamics simulations), and applies a generalized Maxwell construction from the authors' earlier coexistence framework to compute solid-fluid and liquid-gas coexistence curves and the solid-liquid-gas triple point. The central claim is that the full phase diagram, including the strong activity-induced shift of the solid-fluid binodal and the MIPS binodal, can be obtained without invoking equilibrium thermodynamics, and that the predictions quantitatively match simulation data.","tokens_in":20283,"tokens_out":9956,"duration_ms":96416,"significance":"If the approximations are controlled, this would be a substantial advance: a microscopic dynamical route to the nonequilibrium phase diagram of active Brownian spheres, with explicitly stated criteria rather than an equilibrium free-energy construction. The paper is commendably transparent about where approximations enter, and the comparison to Brownian dynamics provides a concrete, falsifiable target. The strengths are the explicit microscopic derivation, the clear articulation of the domain of validity of the coexistence criteria, and the quantitative comparison to simulation data. The main weaknesses are the untested neglect of the crystallinity flux, the uncontrolled orientational closures, the use of many fitted parameters in the equations of state, and the application of high-activity coexistence criteria across the full activity range. These issues are load-bearing for the central quantitative claim, so the contribution is significant but conditional on additional validation.","major_comments":[{"comment":"The assumption |∂_z J_ψ| << |s_ψ| is load-bearing but never tested. The SM states, \"While the flux is clearly not identically zero, we nevertheless approximate it as negligible in comparison to the generation terms,\" and the main text says only that the authors \"expect\" this condition to be approximately met. Since u_ψ in Eq. (5b) is constructed from s_C alone, and Eq. (4) relies on the bulk condition u_ψ^bulk(ρ,ψ*_V)=0, a non-negligible ∂_z J_ψ would change the force balance that fixes ψ*_V and would invalidate the generalized Maxwell construction. I request a direct measurement of J_ψ and s_ψ across a planar coexistence interface, or at least a scaling estimate in ℓ0/D, to justify this central approximation.","section":"SM, Eqs. (S11)-(S12); main text, paragraph after Eq. (2)"},{"comment":"The closures used to derive the crystallinity dynamics are uncontrolled. In particular, the term proportional to <Σ_{i,j}(q_j q_j - q_i q_i)·∂ψ_i/∂r_ij δ(x-r_i)> is eliminated by an isotropic-tensor ansatz, and the third-moment fields B_ψ^1, B_ψ^2, B_ψ^3 are set to zero. These choices directly determine the gradient contribution to u_ψ in Eq. (5b), but no small parameter or numerical check is provided. I ask for an estimate of the omitted terms, for example by evaluating the relevant moments from Brownian dynamics trajectories, to confirm that the derived form of u_ψ is not an artifact of the closure.","section":"SM, Eqs. (S13), (S22)-(S24); main text, Appendix I"},{"comment":"The range of validity of Eq. (4) is not established. The authors state that exact coexistence criteria generally do not exist for ABPs and that E_ρ^int = p_bulk^C is obtained only in the high-activity limit, yet Eq. (4) is applied throughout the phase diagram, including ℓ0/D = 0.27 in Fig. 2(a). The claimed exact reduction to equilibrium criteria at ℓ0/D → 0 is not demonstrated. As implemented in Fig. 5, Eq. (4) with E_ρ^int = p_bulk^C reduces to an integral over pressure rather than the equilibrium equal-area integral over specific volume, so it is not obvious that the passive hard-sphere binodals are recovered. Please show explicitly how Eq. (4) reduces to Eq. (2) in the passive limit and quantify the error introduced by the approximate interfacial Gibbs-Duhem relation at intermediate activities.","section":"Main text, after Eq. (5) and Appendix III; Fig. 5"},{"comment":"The quantitative agreement in Fig. 3 is obtained with roughly 25 fitted coefficients in the equations of state, and the functions ψ*_N, p_bulk^C, p_bulk^act, and φ_ODT are fitted to the same Brownian dynamics simulation data used for the comparison in Fig. 3. The coexistence densities are not directly fitted, so the Maxwell construction is not wholly circular, but the predictive content is partly inherited from the fitted inputs. A sensitivity analysis (for example, refitting the equations of state to a subset of activities, or perturbing the fitted constants) would help establish that the reported agreement is a robust prediction rather than a consequence of the fitting procedure.","section":"SM, Eqs. (S31)-(S35); main text, Fig. 3"}],"minor_comments":[{"comment":"The inset showing φ_ODT(ℓ0/D) is too small to distinguish the simulation symbols from the fitted line, and no error bars are shown; enlarging the inset and adding error bars would improve the presentation.","section":"Fig. 1"},{"comment":"The switching functions t_100 and t_1 are used in Eq. (S31a) before the general definition t_B is given in Eq. (S31b); reordering the definitions would avoid confusion.","section":"SM, Eq. (S31)"},{"comment":"There is a typo in \"interpactice interactions,\" which should read \"interparticle interactions.\"","section":"Main text, Appendix I"},{"comment":"The description of the red region as \"lighter\" in panel (c) than in panel (b) is hard to appreciate without color information; labeling the metastable region or using a different line style would make the point clearer.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The untested crystallinity-flux assumption is the most serious issue; I would not accept the paper in its current form without a direct test or a scaling estimate. The paper also relies heavily on the authors' pre-existing coexistence framework, so the novelty relative to Ref. [31] and the independence of the EOS fitting should be scrutinized. I believe the central claim is defensible and the requested checks are feasible within the manuscript's scope, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first theory I know of that produces both the MIPS binodal and the active solid-fluid binodal for 3D ABPs from a field-theoretic derivation rather than equilibrium thermodynamics. That claim holds up. The paper does something genuinely new: it derives the crystallinity field dynamics from the ABP Fokker-Planck equation via an Irving-Kirkwood procedure, closes the orientational moment hierarchy, and then uses the authors' own coexistence framework to build a phase diagram that quantitatively tracks the simulation binodals and captures the solid-liquid-gas triple point. The comparison to simulation is honest, including the acknowledged failure to recover the passive hard-sphere limit at low activity.\n\nThe soft spots are real but not fatal. The load-bearing approximation is the neglect of the crystallinity flux Jψ in the interface, stated in SM Eqs. (S11)-(S12) and in the main text right after Eq. (2). The authors say they expect |∂z Jψ| << |sψ|, but they never measure Jψ or give a scaling estimate. If that cancellation fails, the generalized Maxwell construction in Eq. (4) loses its foundation. That is the one assumption I would want tested before taking the quantitative diagram as established. A direct interfacial measurement of the flux and source terms is feasible with their simulation setup and would settle it.\n\nThe second soft spot is the EOS load. The phase diagram rests on a chain of fitted functions — ψ*N(φ; ℓ0/D), φODT, the pressure EOSs with ~20 coefficients, the renormalized speed U. Many of those constants are fit to the same Brownian dynamics data that the predictions are compared to. The coexistence densities themselves are not directly fit; they come out of the Maxwell construction, which is a point in the paper's favor. But the reader should not walk away thinking this is a parameter-free prediction. The paper is transparent about this; the SM lists all coefficients and the equations, so the work is reproducible in that sense. Still, a sensitivity analysis would help, and error bars on the simulation data would make the 'quantitatively recapitulates' claim easier to judge.\n\nThe citation pattern is fine. The coexistence framework is their own prior work (Ref. [31]), and the paper builds on it rather than re-deriving it, which is legitimate. The passive limit issue is acknowledged in the text.\n\nWho is this for? Anyone working on active matter phase behavior, especially the crystallization side. It deserves a serious referee: the derivation is coherent, the comparison is encouraging, and the one untested assumption is precisely the kind of thing a good referee can push on. I'd send it to review, with the expectation of a major revision that measures Jψ or at least provides a controlled argument for its neglect.","headline":"Serious theory paper that produces the active solid-fluid binodal and triple point from a microscopic derivation plus semi-empirical EOSs; the referee should push on the untested crystallinity-flux assumption and the fitted EOS load.","tokens_in":20810,"tokens_out":1736,"would_cite":true,"duration_ms":17163,"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":"The paper claims that a nonequilibrium coexistence theory built from density and crystallinity pseudopotentials, completed with simulation-guided equations of state, quantitatively recapitulates the solid-fluid and liquid-gas binodals and…","keywords":["active Brownian particles","crystallization","nonequilibrium coexistence","motility-induced phase separation","phase diagram","Maxwell construction","order parameter","hard spheres"],"falsifier":"In Brownian dynamics simulations of active hard-sphere solid-fluid coexistence at run lengths such as $\\ell_0/D = 1$, $5$, and $15$, measure the crystallinity flux $J_\\psi$ and the crystallinity source $s_\\psi$ across the interface; if $|\\partial_z J_\\psi|$ is comparable to $|s_\\psi|$, the Model A reduction and the generalized Maxwell construction used to draw the phase diagram lose their theoretical foundation.","tokens_in":19634,"feed_emoji":"🧊","tokens_out":10581,"duration_ms":89517,"temperature":0.7,"pith_summary":"This paper aims to establish that the full phase diagram of three-dimensional active Brownian spheres can be computed from nonequilibrium coexistence theory instead of equilibrium thermodynamics. The central objects are two coupled order-parameter fields: a conserved density $\\phi$ and a nonconserved 'crystallinity' field $\\psi$ tracking local crystalline order, with equations of state guided by Brownian dynamics simulations. Applying a generalized Maxwell construction to the derived pseudopotentials yields solid-fluid and liquid-gas (motility-induced phase separation) binodals and a solid-liquid-gas triple point that match simulation across activities, including the marked shift of the solid branch toward close packing at small run lengths. The paper also argues that the naive equilibrium Maxwell construction is formally incorrect at finite activity and can mispredict the coexisting fluid density by roughly an order of magnitude. If the theory is right, the active hard-sphere phase diagram follows from microscopic particle dynamics plus fitted constitutive equations.","feed_headline":"A nonequilibrium theory maps the active Brownian sphere phase diagram","feed_subtitle":"Predicts solid-fluid and MIPS binodals plus the triple point from dynamics and fitted equations of state.","key_machinery":"The load-bearing device is a pair of pseudopotentials derived from the particle equations of motion through a standard coarse-graining procedure: the 'chemical pseudopotential' $u_\\rho = p_C + p_{\\mathrm{act}}^{\\mathrm{bulk}} - \\frac{\\ell_0^2 U}{20}\\,\\partial_z\\!\\left(U\\,\\partial_z\\,p_C^{\\mathrm{bulk}}\\right)$ and the 'crystallinity pseudopotential' $u_\\psi = s_C - \\frac{\\ell_0^2}{24}\\,\\partial_z\\!\\left[U\\,\\partial_z\\!\\left(U\\,s_C^{\\mathrm{bulk}}\\right)\\right]$, where $p_C$ is the conservative interaction pressure, $p_{\\mathrm{act}}^{\\mathrm{bulk}}$ the bulk active pressure, $U$ a dimensionless equation of state for the activity reduction, and $s_C$ the conservative generation of crystallinity. In each bulk phase the crystallinity pseudopotential vanishes along the preferred-crystallinity branch $\\psi^*_V(\\rho)$, and coexistence is found by requiring equality of $u_\\rho^{\\mathrm{bulk}}$ across phases together with the integral $\\int (u_\\rho^{\\mathrm{bulk}} - u_{\\mathrm{coexist}})\\,dE_\\rho^{\\mathrm{int}} \\approx 0$, evaluated with the interfacial Maxwell construction vector $E_\\rho^{\\mathrm{int}} = p_C^{\\mathrm{bulk}}$. The equations of state for $p_C^{\\mathrm{bulk}}$, $U$, and $\\psi^*_V$ are semi-empirical fits constrained by physical limits such as random close packing, fcc close packing, and the equilibrium hard-sphere pressure at low activity.","core_discovery":"The central discovery claimed is that crystallization of active Brownian spheres is captured by a nonequilibrium description in which a conserved density field and a nonconserved crystallinity field both obey steady-state force-balance conditions, and coexistence is located by a generalized Maxwell construction on the resulting pseudopotentials. The bulk phase diagram obtained this way quantitatively recapitulates the simulation phase diagram: the order-disorder transition shifts to higher packing fractions with activity, the solid branch approaches fcc close packing ($\\phi \\approx 0.74$) by $\\ell_0/D \\approx 1$, and the solid-liquid-gas triple point appears near $\\ell_0/D \\approx 18.3$. The paper further claims that exact coexistence criteria do not generally exist for active Brownian spheres, that a solution exists in the passive limit, and that the approximate interfacial criteria used here become well-defined in the high-activity limit yet remain accurate across intermediate and high activities. It also shows that using the equilibrium Maxwell construction naively at finite activity predicts the wrong coexistence partner for the solid above the triple point, overpredicting the fluid density by about 1000%, whereas the nonequilibrium criteria overpredict by about 1-10%.","pith_inferences":["If the pseudopotential construction is generic, the same machinery should extend to other nonconserved order parameters such as nematic order, polycrystalline orientation, or chiral order, giving phase diagrams for active rods, active liquid crystals, and chiral crystals without equilibrium thermodynamics; the paper only treats cubic single-crystal order.","Because exact coexistence criteria do not exist for active Brownian spheres, the accuracy of the phase diagram depends on the specific interfacial Maxwell vector chosen ($E_\\rho^{\\mathrm{int}} = p_C^{\\mathrm{bulk}}$); a natural stress test is to recompute the binodals with alternative vectors such as the total pressure and check how much the triple point moves.","The semi-empirical equations of state are fitted to homogeneous simulation branches, so the theory could be tested more stringently by replacing the hand-fitted forms with direct simulation data on the homogeneous branches and asking whether the generalized Maxwell construction still closes with the same coexistence densities.","The predicted sign change in solid-fluid interfacial polarization above the triple point is an experimentally accessible signature: orientation-resolved imaging of active colloidal or bacterial crystals should show particle polarization pointing toward the fluid at low activity and toward the solid at high activity."],"forward_implications":["Solid-fluid coexistence shifts to higher packing fractions with activity, with the solid branch reaching near-close-packed densities ($\\phi \\approx 0.74$) at run lengths as small as $\\ell_0/D \\approx 1$.","Below the MIPS critical point ($\\ell_0/D \\approx 16.7$) solid-fluid coexistence is the only stable coexistence; between the critical point and the triple point ($\\ell_0/D \\approx 18.3$) solid-fluid and liquid-gas binodals occupy distinct density ranges; above the triple point solid-gas coexistence engulfs the metastable MIPS region.","The equilibrium Maxwell construction is not merely a poor approximation at finite activity; it formally violates the nonequilibrium Gibbs-Duhem relation and can overpredict the coexisting fluid density by about 1000% near the triple point, while the nonequilibrium criteria overpredict by about 1-10%.","The interfacial insertion work for solid-fluid coexistence is predicted to switch sign above the triple point, implying the active-force polarization of the solid-fluid interface reverses direction, whereas the liquid-gas interface always polarizes toward the liquid phase.","The same pseudopotential machinery naturally recovers the MIPS binodal from the same equations of state, so one theory covers both crystallization and motility-induced phase separation."],"supporting_citations":[{"why":"Supplies the simulation phase diagram, solid-fluid coexistence data, and metastability analysis that the theory is built to recapitulate.","marker":"[10]"},{"why":"Provides the dynamical coexistence theory for coupled conserved and nonconserved order parameters and the generalized Maxwell construction used as the coexistence criteria.","marker":"[31]"},{"why":"Provides the mechanical coexistence theory and the active Brownian fluid equations of state that the paper extends to nonzero crystallinity.","marker":"[36]"},{"why":"Defines the bond-orientational order parameter $q_{12}$ used to construct the per-particle crystallinity $\\psi_N$.","marker":"[32]"},{"why":"Supplies the swim-pressure expression used for the bulk active pressure $p_{\\mathrm{act}}^{\\mathrm{bulk}}$.","marker":"[42]"},{"why":"Provides the equilibrium hard-sphere equation of state used as the low-activity reference for the conservative pressure.","marker":"[45]"},{"why":"Supplies the Model A dynamical classification that underlies treating crystallinity as a nonconserved order parameter with vanishing pseudopotential in bulk.","marker":"[37]"},{"why":"Gives the statistical-mechanical coarse-graining procedure used to derive the density and crystallinity field dynamics from the particle equations of motion.","marker":"[39]"}],"fun_headline_variants":["Active Brownian spheres: new theory maps full phase diagram","Nonequilibrium theory predicts active crystal phase diagram","Solid, fluid, gas: theory predicts active sphere phase diagram","Active spheres: new theory captures solid, fluid, and triple point","Triple point predicted: theory of active Brownian sphere freezing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole construction assumes that, inside the interface between coexisting phases, the crystallinity flux is negligible compared with the local production of crystallinity, so the crystallinity field follows Model A dynamics and the coexistence criteria of the coexistence theory apply; the paper expects but does not measure this directly.","fun_headline_variants_meta":{"raw":{"variants":["Active Brownian spheres: new theory maps full phase diagram","Nonequilibrium theory predicts active crystal phase diagram","Solid, fluid, gas: theory predicts active sphere phase diagram","Active spheres: new theory captures solid, fluid, and triple point","Triple point predicted: theory of active Brownian sphere freezing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000681,"raw_usage":{"total_tokens":3089,"prompt_tokens":938,"completion_tokens":2151,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":2068}},"tokens_in":554,"tokens_out":2151,"duration_ms":16412,"temperature":1.0,"reasoning_tokens":2068,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:09:54.788186+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In Brownian dynamics simulations of active hard-sphere solid-fluid coexistence at run lengths such as $\\ell_0/D = 1$, $5$, and $15$, measure the crystallinity flux $J_\\psi$ and the crystallinity source $s_\\psi$ across the interface; if $|\\partial_z J_\\psi|$ is comparable to $|s_\\psi|$, the Model A reduction and the generalized Maxwell construction used to draw the phase diagram lose their theoretical foundation.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the simulation phase diagram, solid-fluid coexistence data, and metastability analysis that the theory is built to recapitulate."},{"cited_title":"equation of state","cited_arxiv_id":null,"evidence_quote":"Provides the dynamical coexistence theory for coupled conserved and nonconserved order parameters and the generalized Maxwell construction used as the coexistence criteria."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the mechanical coexistence theory and the active Brownian fluid equations of state that the paper extends to nonzero crystallinity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the swim-pressure expression used for the bulk active pressure $p_{\\mathrm{act}}^{\\mathrm{bulk}}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Model A dynamical classification that underlies treating crystallinity as a nonconserved order parameter with vanishing pseudopotential in bulk."}],"review_version":1}