{"id":"714b3e23-965e-4cfe-91dc-4a5a4be48bc3","arxiv_id":"1908.10658","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Aligning searchers around a common home undergo a soft transition to collective circular rotation at μ_crit = 2σ²/v0², with closed-form expressions for velocities, variances, and spatial density.","lead":"The paper models groups of animals or robots that share a home and align their motion with neighbors, and it shows that above a critical alignment strength the group spontaneously starts circling the home together. The authors derive formulas for the transition and for the collective circular motion, then check them against computer simulations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Overcritical von Mises closure is imposed: the conditional bracket ⟨1/r−1/r_c+µu_r/v0²⟩_z is assumed to vanish for every z, and Eqs. (34), (52), (53) inherit this unverified ansatz.","rationale":"The reader's CONDITIONAL verdict is appropriate. The central qualitative result, a pitchfork bifurcation of the mean orbital velocity at µ_crit=2σ²/v0², is supported by the independent linear-stability calculation (38) and by the simulation data in Fig. 5(a). My stress-test pass confirms that the weakest point is the closure in Section IV C. The bracket condition is the coefficient of the sin z drift in the marginal angular FPE; dropping it is necessary for the stationary solution to be the symmetric von Mises distribution (34). The authors justify it by the observed z-symmetry and weak z-dependence of P(r|z), but they never measure the conditional average or show it vanishes exactly. Since the decoupling identities (41)–(43) and hence the quantitative overcritical formulas (52)–(53) rest on the von Mises form, the overcritical predictions inherit this empirical closure. This is not fatal: the paper openly states the assumption, and the simulations used as comparison are the same ones that motivate it, so the agreement in Figs. 7–9 is partly circular but still nontrivial. The simulations do show the expected von Mises shape and the predicted scaling of variances, which gives independent support. A secondary issue is a typographical inconsistency in the printed exponent of Eq. (34): as written, exp(µu_z/(σ²v0) sin z) disagrees by a factor v0² with the Bessel argument in Eq. (35), which is the dimensionally consistent expression; this should be corrected but does not change the central claim. The proposed measurement of B(z) would settle whether the main closure is exact or only approximate. If B(z) is nonzero, the verdict should remain conditional with a request to quantify the error; if it is zero, the concern is resolved.","tokens_in":25907,"tokens_out":28864,"duration_ms":301712,"concrete_test":"Reuse the existing simulation data (N=1000, v0=1, κ=1, σ²=0.1, rsens=0.2): in the overcritical steady state for µ=0.3, 0.4, and 0.7, bin particles by z (about 20 bins) and compute B(z)=⟨1/r−1/rc+µu_r(r)/v0²⟩_z with bootstrap error bars. Check whether B(z) is consistent with zero in every bin. If it deviates by more than the statistical uncertainty, solve the stationary marginal FPE (33) with the measured B(z) instead of setting it to zero and compare the resulting u_z² and σ_zz² to Eqs. (52)–(53); a shift larger than about 10% would show that the von Mises overcritical formulas are not quantitatively reliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section IV C derives the marginal angular density by integrating the FPE (23) over r. The coefficient of sin z in the resulting Eq. (33) is ⟨1/r−1/r_c+µu_r/v0²⟩_r, where the average is conditional on z. The authors state that, 'to obey the numerically found asymptotic z-symmetry,' this bracket should vanish for arbitrary z, and they then drop it. This is not derived from the microscopic dynamics or from the stationary transport equations; it is an ansatz motivated by the observed z-symmetry and the statement that P(r|z) depends weakly on z. The von Mises solution (34) and the self-consistent relation (35) exist only if this condition holds exactly. The decoupling identities (41)–(43), which are used to close the mixed-moment balance and to obtain the overcritical mean orbital velocity (52) and orbital variance (53), rely on the von Mises form. If the conditional bracket is nonzero or z-dependent, the angular density is not von Mises, and those quantitative overcritical predictions are approximate. The critical coupling µ_crit=2σ²/v0² is on firmer ground: Eq. (38) comes from linear stability of the undercritical state with uniform z and does not require the bracket condition. The paper does not report any direct measurement of the conditional bracket, so the magnitude of the error in (52)–(53) is unknown.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript analyzes a mean-field model of N identical active searchers that move at constant speed v0, are attracted to a common home by a torque of strength κ, align their headings with neighbors within a sensing radius, and experience angular white noise of strength σ. The central analytical results are a pitchfork transition of the mean orbital velocity at the critical coupling µ_crit = 2σ²/v0² (Eq. 36), the stationary marginal spatial density ρ(r) = C [r exp(-κr/v0)]^(v0²µ/σ² - 1) (Eq. 48), and closed-form expressions for the overcritical radial variance σ²_rr = σ²/µ (Eq. 46) and orbital variance σ²_zz = 2σ⁴/(µ²v0²) (Eq. 53). These are derived from a nonlinear mean-field Fokker-Planck equation (Eq. 23) using a von Mises closure for the angular distribution, and are compared with agent-based simulations (Figs. 5, 7-9). The paper also derives overdamped Smoluchowski equations for the marginal density in both regimes.","tokens_in":26232,"tokens_out":9415,"duration_ms":94531,"significance":"If the reported results hold, the paper provides a rare analytically tractable description of a collective-motion transition in a central-place foraging context. The critical coupling µ_crit = 2σ²/v0² is parameter-free and follows from a simple linear-stability analysis of the orientationally uniform state (Eq. 38), independently of the von Mises closure. The explicit formulas for the stationary densities and velocity variances are falsifiable and are tested against simulations, with generally reasonable agreement. The comparison between von Mises and Gaussian decoupling schemes (Appendix B) is instructive and shows which closure better captures the periodicity of the angular variable. The manuscript is clearly written and the simulation data are presented in a reproducible manner.","major_comments":[{"comment":"The derivation of the von Mises angular density (34) and the self-consistent mean orbital velocity (35) rests on the condition that the conditional average ⟨1/r − 1/r_c + µu_r/v0²⟩_r vanishes for every z, which the authors impose 'to obey the numerically found asymptotic z-symmetry.' This condition is not derived from the stationary Fokker-Planck equation or from the microscopic dynamics, and the manuscript reports no direct numerical measurement of this conditional bracket. Since the decoupling identities (41)-(43) and the overcritical expressions for the mean orbital velocity (52) and orbital variance (53) inherit the von Mises form, the error incurred by this ansatz is uncontrolled. The authors should either derive the condition from the assumed z-symmetry and the stationary FPE, or provide a direct simulation test showing that the conditional average is small in the overcritical regime and quantify the resulting error in (52)-(53).","section":"Section IV C, Eqs. (33)-(35)"},{"comment":"The paper presents two distinct analytical predictions for the overcritical orbital velocity: the self-consistent von Mises result (35) and the transport-equation result (52). These two expressions visibly differ for moderate µ (e.g., µ ≈ 0.4 in Fig. 7), and the simulation data lie between them. The authors attribute the difference to finite particle number, but no finite-size scaling or error estimates are provided to support this attribution. Because the difference between the two analytical predictions is of the same order as the simulation scatter, the manuscript should quantify the finite-size corrections or discuss which closure error is responsible, rather than invoking finite N qualitatively.","section":"Section IV D 2 and Fig. 7"}],"minor_comments":[{"comment":"The definitions of the variances in Eq. (27) contain a typo: the integrals should involve v0² cos²(z) and v0² sin²(z), respectively, for dimensional consistency with the subtracted u²_r ρ and u²_z ρ terms; the current expressions with a single factor v0 are dimensionally inconsistent.","section":"Eq. (27)"},{"comment":"The notation 'ρz(z.t)' in the sentence before Eq. (33) should read 'ρz(z,t)'.","section":"Section IV C, Eq. (33)"},{"comment":"The abstract contains a typo: 'rhytmically' should be 'rhythmically'.","section":"Abstract"},{"comment":"The caption of Fig. 7(b) is confusing: 'according to Eq. (46) as red dotted line as dashed line' mixes line styles and equations; please rephrase to clearly identify each curve.","section":"Fig. 7 caption"},{"comment":"Reference [38] appears malformed ('arXiv:q-bio 325, 0404018 (2004)') and should be corrected.","section":"Reference [38]"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the main phase-transition result is solid, with µ_crit = 2σ²/v0² derived cleanly from linear stability. The primary obstacle is the von Mises closure in Section IV C: the conditional average is imposed rather than derived, and the overcritical quantitative predictions (52)-(53) inherit this assumption. Please ask the authors to provide a numerical test of the conditional bracket or a derivation from the stationary FPE, and to quantify the discrepancy between Eqs. (35) and (52). If these points are resolved, the paper is likely acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this paper is a solid mean-field treatment of aligning central-place searchers, and its main transition survives scrutiny; the quantitative overcritical predictions rest on a closure assumption that is plausible but not derived.\n\nWhat's new: the combination of a homing torque, constant-speed motion, and pairwise alignment, analyzed through a nonlinear Fokker-Planck equation. They get a clean pitchfork condition µ_crit = 2σ²/v0², closed-form stationary densities and velocity moments, and a Smoluchowski limit. The critical coupling comes from two independent routes: linear stability of the uniform state (Eq. 38) and the self-consistent von Mises equation (Eq. 35). That the two agree and match simulations is genuinely good evidence.\n\nThe paper also does useful housekeeping. The comparison with Gaussian decoupling clarifies why the von Mises closure is better. The finite-size appendix addresses the N→∞, rsens→0 assumptions honestly. The Daphnia connection is motivating but not oversold.\n\nSoft spots, in proportion. The main one is in Section IV C. To get the von Mises angular density, they require the conditional bracket ⟨1/r −1/r_c + µu_r/v0²⟩_r to vanish for every z. That is an imposed ansatz, motivated by weak z-dependence seen in simulations, not derived from the dynamics. Later they verify that the spatial average of that bracket vanishes in both regimes, but that is a weaker statement than the conditional one. So Eqs. (34), (52), and (53) inherit this assumption. I don't think it breaks the central claim—µ_crit is independent of it, and the simulation agreement is decent—but the paper should either derive the condition from the transport equations or add a direct simulation measurement of the conditional bracket and a statement of how big the error is. Minor points: no error bars in Figs. 5, 7–9, no code/data deposit, and a few typos. None of these affect the qualitative picture.\n\nWho this is for: people working on active matter, collective motion, and mean-field kinetics of swarms. It is a subfield contribution, not a major advance, but it is a clean one and deserves a serious referee.\n\nRecommendation: send to peer review. Ask the authors to justify or measure the closure condition, and to add reproducibility details. I would accept it after those revisions.","headline":"Clean mean-field treatment of aligning central-place searchers with a robust critical coupling, though the quantitative overcritical predictions rest on an imposed closure condition that should be justified or tested.","tokens_in":26735,"tokens_out":2901,"would_cite":true,"duration_ms":29659,"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":"Above a critical alignment, searchers spontaneously rotate as a group around their home.","keywords":["active searchers","central place foraging","alignment interaction","collective rotation","mean-field Fokker-Planck equation","von Mises distribution","pitchfork bifurcation","Daphnia"],"falsifier":"Measure, in simulation or tracking, the squared mean orbital velocity as a function of alignment $\\mu$ at fixed $v_0$ and $\\sigma$; the theory fails if the onset is not at $\\mu=2\\sigma^2/v_0^2$ or if the overcritical orbital variance does not scale as $2\\sigma^4/(\\mu^2 v_0^2)$. Alternatively, condition on the angle $z$ and compute $\\langle 1/r-1/r_c+\\mu u_r/v_0^2\\rangle_z$: if it is not near zero across $z$, the assumed symmetry behind the von Mises density is broken.","tokens_in":2016,"feed_emoji":"🌀","tokens_out":2212,"duration_ms":92633,"temperature":0.7,"pith_summary":"This paper asks what happens when many active searchers that each drift home with constant speed also align their headings with neighbours. It claims a sharp, second-order transition at coupling $\\mu_{\\rm crit}=2\\sigma^2/v_0^2$: below it, searchers move in independent random loops and the angular distribution of headings is uniform; above it, the ensemble spontaneously rotates collectively around the home, clockwise or anticlockwise. The claim matters because it turns a minimal model of central-place foraging into a solvable kinetic theory, with explicit predictions for mean orbital velocity, velocity variances, and the stationary spatial density around the home. These predictions match particle simulations over a broad range of coupling strengths.","feed_headline":"Alignment makes homing foragers circle their home in unison","feed_subtitle":"Below a sharp threshold each searcher loops alone; above it the whole group rotates one way around the home.","key_machinery":"The machine is the nonlinear mean-field Fokker-Planck equation for the joint density of the distance $r$ and the angle $z=\\varphi-\\beta$ between a searcher's heading and its position vector. With a small sensing radius and rotational symmetry, the alignment force reduces to $\\mu(u_z\\cos z-u_r\\sin z)$, which closes the reduced $(r,z)$ dynamics. The stationary angular density is a von Mises distribution $\\rho_z(z)\\propto\\exp\\left(\\mu u_z v_0\\sin z/\\sigma^2\\right)$, whose self-consistency equation yields the pitchfork bifurcation at $\\mu_{\\rm crit}=2\\sigma^2/v_0^2$. This von Mises form is then used to decouple third-order moments in the transport equations, producing the overcritical velocities, variances, and spatial density.","core_discovery":"For an ensemble of $N$ searchers moving at speed $v_0$, attracted to a common home with strength $\\kappa$, and aligning with neighbours inside a sensing radius with strength $\\mu$, the paper establishes that collective circular motion sets in only above $\\mu_{\\rm crit}=2\\sigma^2/v_0^2$. In the ordered state the mean orbital velocity obeys $u_z^2=v_0^2-(\\sigma^2/\\mu)(1+2\\sigma^2/(\\mu v_0^2))$, the radial velocity variance is $\\sigma_{rr}^2=\\sigma^2/\\mu$, and the orbital variance is $\\sigma_{zz}^2=2\\sigma^4/(\\mu^2 v_0^2)$; the stationary spatial density takes the form $\\rho(r)=C\\left[r\\exp(-\\kappa r/v_0)\\right]^{v_0^2\\mu/\\sigma^2-1}$. Below the threshold the variances are equal to $v_0^2/2$ and the density reduces to the single-searcher exponential $\\rho(r)=(r/r_c^2)\\exp(-r/r_c)$ with $r_c=v_0/\\kappa$. The paper verifies these predictions against simulations of up to 2000 particles and finds good agreement.","pith_inferences":["The paper does not draw this out, but the effective diffusion coefficient $D_{\\rm eff}=v_0^2/|\\mu-\\mu_{\\rm crit}|$ from the Smoluchowski limit implies critical slowing of spatial spreading near the transition, which could be measured independently of the velocity moments.","Because the overcritical moments and the threshold do not contain $\\kappa$, one can test the theory by varying the homing strength: the onset of rotation should stay put while the spatial density narrows or widens.","The reduced angle $z=\\varphi-\\beta$ between heading and position vector is a natural single-particle observable; tracking data could extract $\\mu u_z v_0/\\sigma^2$ from the fitted von Mises shape and then predict the onset from the noise level alone.","For a finite sensing radius the approximation $\\beta_i\\approx\\beta_j$ used in the alignment force will break down near the home, so finite-$N$ swarms should show apparent critical couplings that shift with the sensing radius, suggesting a finite-size scaling check."],"forward_implications":["Below $\\mu_{\\rm crit}=2\\sigma^2/v_0^2$, the angular distribution is uniform, the two velocity variances equal $v_0^2/2$, and the spatial density is the same exponential form as for an isolated searcher.","Above the threshold, the squared mean orbital velocity grows continuously from zero and saturates at $v_0^2$; the radial velocity variance is $\\sigma^2/\\mu$ and the orbital variance is $2\\sigma^4/(\\mu^2 v_0^2)$, so orbital fluctuations are suppressed more strongly than radial ones.","In the rotating state the stationary spatial density becomes $\\propto\\left[r\\exp(-\\kappa r/v_0)\\right]^{v_0^2\\mu/\\sigma^2-1}$: it vanishes at the home for strong coupling, and its peak stays at $r_{\\max}=v_0/\\kappa$ independent of $\\mu$.","The transition is a second-order pitchfork: at $\\mu=\\mu_{\\rm crit}$ the overcritical density coincides with the undercritical one, and the chosen rotation direction depends on the initial state, as reported for Daphnia swarms around a light shaft.","In the overdamped limit both regimes obey a Smoluchowski equation with an effective diffusion coefficient that depends on distance from the threshold; undercritical spreading speeds up with $\\mu$ while overcritical spreading slows."],"supporting_citations":[{"why":"supplies the single-searcher search-and-return model and its stationary exponential spatial density, which the undercritical ensemble reproduces.","marker":"[26]"},{"why":"companion derivation of the single-searcher kinetics whose transport-equation approach the paper generalizes to the interacting ensemble.","marker":"[27]"},{"why":"reports Daphnia experiments and simulations showing unidirectional rotation around a light shaft, the motivating empirical scenario.","marker":"[36]"},{"why":"provides the fluid-drag argument that aligning interactions can arise from the comoving medium around swarming Daphnia.","marker":"[38]"},{"why":"defines alignment as the arrayal force among animal groups, the interaction adopted here.","marker":"[47]"},{"why":"introduces the Vicsek-type aligning rule that the heading dynamics uses.","marker":"[48]"},{"why":"supplies the mean-field Fokker-Planck and moment-equation framework for active particles used in the analytical derivation.","marker":"[52]"},{"why":"the self-consistent equation for the orbital velocity is of Kuramoto type, whose bifurcation structure fixes the critical coupling.","marker":"[55]"},{"why":"provides the von Mises distribution used for the stationary angular density and moment closure.","marker":"[78]"}],"fun_headline_variants":["Aligning foragers circle home together past a sharp threshold","Critical alignment makes homing searchers rotate as one","Homebound searchers shift to collective circling at critical strength","Below threshold foragers loop alone; above it they circle in sync","Searchers align to orbit home together above a critical threshold"],"cache_read_input_tokens":28800,"weakest_assumption_plain":"The load-bearing premise is that the averaged combination $\\langle 1/r-1/r_c+\\mu u_r/v_0^2\\rangle_r$ vanishes for every angle $z$, imposed to match the numerically observed $z$-symmetry rather than derived, so if that bracket is only approximately zero the von Mises density and all overcritical moments inherit the error, although the critical coupling itself survives because it also follows from linear stability.","fun_headline_variants_meta":{"raw":{"variants":["Aligning foragers circle home together past a sharp threshold","Critical alignment makes homing searchers rotate as one","Homebound searchers shift to collective circling at critical strength","Below threshold foragers loop alone; above it they circle in sync","Searchers align to orbit home together above a critical threshold"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00066,"raw_usage":{"total_tokens":3093,"prompt_tokens":1098,"completion_tokens":1995,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":714,"completion_tokens_details":{"reasoning_tokens":1911}},"tokens_in":714,"tokens_out":1995,"duration_ms":15814,"temperature":1.0,"reasoning_tokens":1911,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:37:52.468275+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, in simulation or tracking, the squared mean orbital velocity as a function of alignment $\\mu$ at fixed $v_0$ and $\\sigma$; the theory fails if the onset is not at $\\mu=2\\sigma^2/v_0^2$ or if the overcritical orbital variance does not scale as $2\\sigma^4/(\\mu^2 v_0^2)$. Alternatively, condition on the angle $z$ and compute $\\langle 1/r-1/r_c+\\mu u_r/v_0^2\\rangle_z$: if it is not near zero across $z$, the assumed symmetry behind the von Mises density is broken.","supporting_citations":[{"cited_title":"Noetel, V","cited_arxiv_id":null,"evidence_quote":"supplies the single-searcher search-and-return model and its stationary exponential spatial density, which the undercritical ensemble reproduces."},{"cited_title":"Noetel, V","cited_arxiv_id":null,"evidence_quote":"companion derivation of the single-searcher kinetics whose transport-equation approach the paper generalizes to the interacting ensemble."},{"cited_title":"Ordemann, G","cited_arxiv_id":null,"evidence_quote":"reports Daphnia experiments and simulations showing unidirectional rotation around a light shaft, the motivating empirical scenario."},{"cited_title":"Erdmann, W","cited_arxiv_id":null,"evidence_quote":"provides the fluid-drag argument that aligning interactions can arise from the comoving medium around swarming Daphnia."},{"cited_title":"Okubo and S","cited_arxiv_id":null,"evidence_quote":"defines alignment as the arrayal force among animal groups, the interaction adopted here."},{"cited_title":"Romanczuk, M","cited_arxiv_id":null,"evidence_quote":"supplies the mean-field Fokker-Planck and moment-equation framework for active particles used in the analytical derivation."},{"cited_title":"Kuramoto, Chemical Oscillations, Waves, and Turbulence (Springer-Verlag, New York, 1984)","cited_arxiv_id":null,"evidence_quote":"the self-consistent equation for the orbital velocity is of Kuramoto type, whose bifurcation structure fixes the critical coupling."},{"cited_title":"Stratonovich, Topics in the Theory of Random Noise , Vol","cited_arxiv_id":null,"evidence_quote":"provides the von Mises distribution used for the stationary angular density and moment closure."}],"review_version":1}