{"id":"370f7699-2d21-490e-b6b6-360f3823ac20","arxiv_id":"1908.08842","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a two-dimensional chemotaxis model with density-dependent pressure, the aggregation transition is discontinuous with hysteresis, and aggregates form near the boundary.","lead":"Bacteria that attract each other but also repel when crowded can undergo a sudden, first-order aggregation transition on a Petri dish, according to this model study. The authors show the transition is discontinuous, has hysteresis, and forms clusters near the boundary, which may bear on how bacterial aggregates form in blood vessels.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Both numerical methods support the discontinuity only at P=M=3 mode truncation, with no convergence test; higher modes could alter the clustered branch or hysteresis.","rationale":"The reader's weakest assumption is exactly the absence of mode-truncation convergence, and I find that this is the most load-bearing concern. The Lyapunov functional derivation (Appendix A) and the linear stability analysis (Appendix C) are internally consistent; the two-mode approximation in Eq. (8) is explicitly illustrative and is not the main quantitative evidence. The decisive evidence for the discontinuous-transition claim is the coexistence of branches in the gradient descent and the hysteresis in the direct simulation. But both calculations truncate the expansion at P=M=3, so they cannot distinguish a genuine first-order transition in the PDE from an artifact of a severely restricted Galerkin subspace. The boundary-localized aggregate is precisely the kind of solution that high-order modes can strongly affect. Since the reader already conditions acceptance on a convergence study, I do not change the verdict. My recommendation is to keep the paper CONDITIONAL pending the explicit P,M-convergence test.","tokens_in":13960,"tokens_out":4365,"duration_ms":47762,"concrete_test":"Repeat the spectral simulation of Fig. 3 and the gradient-descent minimization of Fig. 2 with P=M=5, 7, and 9, keeping all other parameters fixed, and also halve the Euler time step at the largest truncation. For each truncation, record the order-parameter jump size and location, the up-sweep and down-sweep spinodal values of f0, and the amplitudes of the newly added modes near the aggregate. If the hysteresis loop persists with roughly unchanged width and thresholds, and the high-order mode amplitudes remain small relative to R11, the discontinuous-transition claim is supported in the PDE limit. If the loop narrows, shifts by more than the f0 resolution, or disappears as the truncation grows, the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the aggregation transition is discontinuous in the full two-dimensional PDE. The numerical evidence for this claim is the hysteresis loop in the spectral simulation (Fig. 3) and the two branches in the Lyapunov-functional minimization (Fig. 2), but both use the same Fourier-Bessel truncation P=M=3 (Eqs. 10-12). These two methods are not independent with respect to truncation error: they project onto the same low-order subspace. The aggregate that forms sits near the boundary (Fig. 2d), so resolving it should require many high-order angular and radial modes; with only three radial modes per angular order, the apparent bistability could be a projection artifact. No convergence study in P or M, no error bound, and no comparison with larger truncations is reported. The paper's own radial-symmetry analysis in Sec. IV shows that when the allowed mode structure is changed, the transition becomes continuous, underscoring that the discontinuity is sensitive to the modal truncation. Therefore the load-bearing assumption, that neglected modes do not change the stability of the clustered state or the shape of the Lyapunov landscape, is unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a two-dimensional Patlak-Keller-Segel model with density-dependent pressure, where the effective pressure is truncated at quadratic order in density. The authors construct an exact Lyapunov functional and analyze it by expanding the density and chemoattractant fields in Fourier-Bessel modes on a disk. For the linear-pressure case (κ=0) the functional decouples into independent modes and the transition to aggregation is discontinuous but without hysteresis. For nonlinear pressure (κ>0), a two-mode approximation reveals bistability, and numerical minimization of the truncated Lyapunov functional together with direct spectral simulations, both using P=M=3 modes, show a discontinuous transition with hysteresis. A radially symmetric restriction (p=0) is also analyzed and found to yield a continuous transition. The paper concludes that the full aggregation transition is discontinuous and interprets the Lyapunov functional as a free energy.","tokens_in":14124,"tokens_out":6526,"duration_ms":63535,"significance":"If the central claim holds, the paper provides an interesting example of a discontinuous, hysteretic aggregation transition in a chemotaxis model with a mechanical regularization, and it demonstrates the usefulness of Lyapunov-functionals as free-energy analogues for active matter. The derivations in the appendices are careful and transparent: the Lyapunov functional is derived from the equations rather than assumed, and the linear stability threshold is a derived quantity. The numerical methods are clearly described and the hysteresis signal is visually consistent. However, the significance is tempered by the fact that the main quantitative evidence for the full-PDE claim rests entirely on a low-order Fourier-Bessel truncation (P=M=3) without any convergence test, and the paper's own radial-symmetry analysis shows that the transition order can change when the allowed mode structure is changed.","major_comments":[{"comment":"The central claim that the aggregation transition is discontinuous in the full two-dimensional PDE is supported only by calculations with P=M=3 modes. No convergence study in P and M is reported, and the two numerical methods (gradient-descent minimization of the Lyapunov functional and the spectral simulation) are not independent with respect to truncation because they project onto the same low-order subspace. Since the aggregate forms near the boundary, where the first unstable mode J_1(j'_{11}r/l) peaks, an adequate representation may require many higher-order radial and angular modes. Please add a systematic convergence test (e.g., increasing P and M to 4, 5, or 6) for both the Lyapunov-minimization branches and the hysteresis loop, or provide a rigorous error bound showing that neglected modes do not alter the stability of the clustered state or the shape of the Lyapunov landscape.","section":"§III, Figs. 2 and 3"},{"comment":"The paper's own radially symmetric analysis shows that when the allowed mode structure is restricted to p=0 and m=1,2, the transition becomes continuous, with the order parameter increasing linearly near threshold. This is explicitly acknowledged in the Discussion as a case where 'the transition behavior may change if we suppress aggregation at the boundary.' This example demonstrates that the order of the transition is sensitive to the modal truncation, which reinforces the need for a convergence check in the full p>0 case. Absent such a check, the discontinuous-transition claim is not established beyond the truncated model.","section":"§IV, Eq. (13)"},{"comment":"The dispersion relation is misprinted: the line reads 'η^2 + [ ... ] + k^2(...)(...)' with no linear term in η and no '= 0'. This makes it difficult to verify the derivation of the stability condition in Eq. (C4). Please correct the equation.","section":"Appendix C, Eq. (C3)"}],"minor_comments":[{"comment":"The caption states the minimization uses P=M=3, but the rendered panel (c) appears to include a legend with 'P=M=5' alongside '3'. If panel (c) shows mode coefficients for two different truncations, the text and caption should state this explicitly; if not, the stray 'P=M=5' label should be removed.","section":"Fig. 2 caption"},{"comment":"The matrix notation in the linearized equation is slightly unconventional; the diffusion term is written as a matrix acting on ∇^2 of the perturbation vector. This is understandable, but a sentence explaining the matrix structure would improve clarity.","section":"Appendix C, Eq. (C2)"},{"comment":"The two-mode approximation is introduced as an illustration, and the conclusion is explicitly limited to that truncation. This is fine, but the transition from 'the answer is negative, at least within the two-mode approximation' to the later claim about the full system would benefit from a clearer statement that the full-system conclusion relies on the numerical results.","section":"Sec. II, Eq. (8)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central idea is potentially interesting. The main concern is the missing convergence study for the mode truncation; this is fixable and should be requested in revision. The self-citation for the parameter choice χ0=4 is not an issue. The authors may also want to address the apparent ambiguity in Fig. 2(c) regarding P=M=5."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main new result here is that the two-dimensional PKS model with quadratic density-dependent pressure undergoes a discontinuous aggregation transition with hysteresis on a disk, and that the aggregate forms near the boundary. The Lyapunov-functional machinery is handled well: the functional is derived, not assumed, the κ=0 case is solved exactly in the Fourier-Bessel basis, and the two-mode approximation gives a transparent mechanism for bistability. The linear stability threshold matches the loss of stability of the homogeneous state, and the two numerical approaches (gradient descent on the Lyapunov functional and spectral simulation of the PDE) agree with each other and with the analytic picture. The hysteresis loop in Fig. 3 is a genuine predictive signal, and the paper is honest about the thermodynamic interpretation being a heuristic analogy.\n\nWhere I waver is the mode truncation. Both numerical methods use P=M=3, i.e., they project onto the same low-order subspace. That means they are not truly independent checks: they share the same truncation error. The aggregate sits near the boundary, and resolving that kind of localized structure typically requires higher angular and radial modes, so the apparent bistability could in principle be a projection artifact. The paper reports no convergence study in P or M, no error bound, and no test with a larger truncation. That is a real gap, and it is load-bearing because the central claim is about the full PDE, not about a three-mode toy model.\n\nThe radial-symmetry analysis in Sec. IV does not worry me as much as the stress-test suggests. Restricting to p=0 is a different physical sector, not simply a lower truncation, and the paper explicitly notes that suppressing boundary aggregation changes the transition to continuous. That is a known limitation, not a contradiction. But it does reinforce that transition order is sensitive to which modes are allowed, which makes the missing convergence test more salient.\n\nOverall: the analytic work is solid, the numerical evidence is consistent within the truncation, and the question is clearly posed and answered. The paper deserves a serious referee, but a referee should push for a convergence check in P and M before accepting the discontinuous-transition claim as a statement about the full PDE. I would cite this for the Lyapunov-functional treatment and the two-mode picture, and would bring it to a reading group interested in chemotaxis or active matter.","headline":"Careful Lyapunov-functional analysis of a chemotaxis model claims a discontinuous aggregation transition with hysteresis, but the numerical evidence rests on a single low-order mode truncation that needs a convergence check.","tokens_in":14674,"tokens_out":1641,"would_cite":true,"duration_ms":17288,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q92","92C17","35B36"],"pacs":["64.60.-i"],"model":"deepseek-v4-flash","headline":"The aggregation transition in a two-dimensional chemotaxis model with density-dependent pressure is discontinuous: the order parameter jumps, the Lyapunov landscape has two coexisting minima, and direct simulation exhibits hysteresis.","keywords":["chemotaxis","Patlak-Keller-Segel","density-dependent pressure","Lyapunov functional","discontinuous phase transition","hysteresis","Fourier-Bessel expansion","aggregation"],"falsifier":"Run the full PDE Eq. (1) with many Fourier-Bessel modes or a high-resolution finite-element scheme across $f_0=2.15$ through $2.26$ and measure the order parameter $A$ on slow up- and down-ramps; if the aggregate dissolves immediately below $f_0^{\\mathrm{unstable}}\\approx2.195$, or if $A$ grows continuously instead of jumping, the discontinuous-transition claim fails.","tokens_in":13704,"feed_emoji":"🦠","tokens_out":6846,"duration_ms":67217,"temperature":0.7,"pith_summary":"The paper studies a two-dimensional Patlak-Keller-Segel model in which chemotactic attraction is balanced by a density-dependent pressure that is quadratic in the local density. Its central claim is that the transition from a uniform population to an aggregated one is discontinuous: the order parameter $A = \\max\\rho - \\min\\rho$ jumps at a threshold, the Lyapunov landscape develops two coexisting minima, and direct numerical simulation shows hysteresis, so an existing aggregate survives when the attraction strength drops below the linear instability threshold $f_0^{\\mathrm{unstable}}\\approx2.195$. A two-mode Fourier-Bessel approximation shows why: the quadratic pressure term couples modes, and once the lowest mode is excited, the coupling can create a metastable aggregate branch before the homogeneous solution loses stability. If the claim is right, bacterial aggregation in confined geometries is an abrupt, history-dependent event rather than a smooth condensation.","feed_headline":"Chemotactic clumping is an abrupt, history-dependent jump","feed_subtitle":"Lyapunov analysis and simulation show clumps persist even after attraction falls below the instability threshold.","key_machinery":"The central object is the Lyapunov functional $W$ of the system, equation (3), which decreases along every trajectory and is interpreted as a free energy. The argument is carried by expanding $\\rho$ and $c$ in Fourier-Bessel modes $J_p(j'_{pm} r/l)\\cos[p(\\theta-\\eta_{pm})]$ on the disk, where $j'_{pm}$ are zeros of the derivative of the Bessel function; this makes $W_{\\mathrm{st}}$ a sum over modes plus a coupling term. The crucial identity is the mode-coupling coefficient $I_{112}=\\int_0^1 x J_1^2(j'_{11}x)J_2(j'_{21}x)\\,dx\\approx0.0474258$, which appears in the cubic term $\\kappa I_{112}z_{11}^2 z_{21} R_{11}^2 R_{21}\\cos 2(\\eta_{11}-\\eta_{21})$ in Eq. (8). This term is what produces a second local minimum before the homogeneous solution becomes linearly unstable, and the same Bessel orthogonality relations are used to derive the quadratic form for $\\kappa=0$ and the radial-symmetry quartic analysis.","core_discovery":"On the paper's own terms, the main discovery is that adding a quadratic term to the effective pressure changes the character of the aggregation instability expected from the classical PKS model. For $\\kappa=D_1/D_0=0$, each Fourier-Bessel mode contributes independently to the Lyapunov functional and the first excited mode $j'_{11}$ produces a jump with no hysteresis. For $\\kappa>0$, the cubic coupling term $\\propto R_{11}^2 R_{21}$ makes the landscape bistable: the aggregate branch coexists with the uniform branch, so minimization of $W_{\\mathrm{st}}$ gives a crossing of branches and a discontinuous jump of $A$, while spectral simulation of the PDE shows that an aggregate formed at $f_0=2.26$ dissolves only at $f_0\\lesssim2.15$, well below $f_0^{\\mathrm{unstable}}\\approx2.195$. The authors also show the boundary matters: an aggregate forms near the disk edge, and under imposed radial symmetry the same model gives a continuous transition with order parameter growing linearly with slope $\\propto \\kappa^{-1}$.","pith_inferences":["If the discontinuity survives the full infinite-mode limit, then a microfluidic or swarm experiment that slowly ramps chemoattractant up and down should observe a hysteresis loop in population clumpiness, not just in the model; this is directly testable with bacterial or synthetic chemotactic swimmers.","The paper's own lack of convergence checks for the truncation sizes $P$ and $M$ means the strongest plausible threat is that high-order modes round off the jump or destabilize the aggregate branch; a high-resolution spectral or finite-element study of Eq. (1) is the natural check.","For active-matter thermodynamics, the two-branch structure of $W_{\\min}$ suggests defining a latent-like quantity from the derivative of the minimized functional, which could support a Maxwell equal-area construction for the aggregation transition.","Because Appendix A derives the Lyapunov functional for an arbitrary $\\rho$-dependent pressure, the same two-mode criterion could classify supercritical versus subcritical aggregation for other pressure virial forms, not just the quadratic truncation studied here."],"forward_implications":["Aggregation is not gradual: as chemotactic strength $f_0$ rises through the threshold, a compact aggregate forms at the boundary with a sudden jump in the order parameter $A$.","The transition is hysteretic: an aggregate created at high $f_0$ remains intact down to about $f_0\\simeq2.15$, so reversing the control parameter does not retrace the forward path.","The quadratic pressure term is essential: for linear pressure ($\\kappa=0$) the mode analysis gives a coexistence-free jump, while $\\kappa>0$ creates bistability and hysteresis.","If aggregation is forced to be radial about the disk center, the transition becomes continuous with $A\\propto \\kappa^{-1}$ above threshold, so boundary geometry can change the order of the transition.","The Lyapunov functional can serve as a free energy; a derivative of the minimized functional with respect to disk area acts like a thermodynamic pressure, suggesting a thermodynamic reading of the aggregation transition."],"supporting_citations":[{"why":"Supplies the original derivation of the Lyapunov functional that is the paper's main analytical tool.","marker":"[15]"},{"why":"Introduces the density-dependent-pressure PKS variant in which the authors set their model.","marker":"[16]"},{"why":"Provides the experimental polynomial fit for pressure in insect swarms that motivates the quadratic truncation $\\Pi = D_0\\rho + D_1\\rho^2/2$.","marker":"[17]"},{"why":"Supplies the parameter values $\\chi_0=4$, $\\rho_0=l=D_0=\\nu_0=g_0=\\kappa=1$ used in the numerical calculations.","marker":"[27]"},{"why":"Supplies the Fourier-Bessel expansion and orthogonality relations used to derive the mode-decomposed Lyapunov functional.","marker":"[32]"},{"why":"Gives the classical PKS model whose discontinuous transition under radial symmetry is contrasted with the continuous radial result of this paper.","marker":"[11]"}],"fun_headline_variants":["Discontinuous transition in bacterial aggregation","Abrupt clumping with density-dependent pressure","Bistable bacterial clumps: abrupt jump with hysteresis","Boundary-driven jump in chemotactic aggregation","Nonlinear pressure makes bacterial clumping discontinuous"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on truncating the infinite Fourier-Bessel expansion to a small number of modes ($P=M=3$ in the numerics, two modes in the analytic picture), and the paper never checks that the jump and hysteresis survive when more modes are included.","fun_headline_variants_meta":{"raw":{"variants":["Discontinuous transition in bacterial aggregation","Abrupt clumping with density-dependent pressure","Bistable bacterial clumps: abrupt jump with hysteresis","Boundary-driven jump in chemotactic aggregation","Nonlinear pressure makes bacterial clumping discontinuous"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000514,"raw_usage":{"total_tokens":2473,"prompt_tokens":901,"completion_tokens":1572,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":517,"completion_tokens_details":{"reasoning_tokens":1504}},"tokens_in":517,"tokens_out":1572,"duration_ms":12455,"temperature":1.0,"reasoning_tokens":1504,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:27:18.007965+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the full PDE Eq. (1) with many Fourier-Bessel modes or a high-resolution finite-element scheme across $f_0=2.15$ through $2.26$ and measure the order parameter $A$ on slow up- and down-ramps; if the aggregate dissolves immediately below $f_0^{\\mathrm{unstable}}\\approx2.195$, or if $A$ grows continuously instead of jumping, the discontinuous-transition claim fails.","supporting_citations":[{"cited_title":"Hillen and K","cited_arxiv_id":null,"evidence_quote":"Supplies the original derivation of the Lyapunov functional that is the paper's main analytical tool."},{"cited_title":"Gamba, D","cited_arxiv_id":null,"evidence_quote":"Introduces the density-dependent-pressure PKS variant in which the authors set their model."},{"cited_title":"Kowalczyk, J","cited_arxiv_id":null,"evidence_quote":"Provides the experimental polynomial fit for pressure in insect swarms that motivates the quadratic truncation $\\Pi = D_0\\rho + D_1\\rho^2/2$."},{"cited_title":"Fatkullin, Nonlinearity 26, 81 (2013)","cited_arxiv_id":null,"evidence_quote":"Supplies the parameter values $\\chi_0=4$, $\\rho_0=l=D_0=\\nu_0=g_0=\\kappa=1$ used in the numerical calculations."},{"cited_title":"Madzvamuse, H","cited_arxiv_id":null,"evidence_quote":"Supplies the Fourier-Bessel expansion and orthogonality relations used to derive the mode-decomposed Lyapunov functional."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the classical PKS model whose discontinuous transition under radial symmetry is contrasted with the continuous radial result of this paper."}],"review_version":1}