{"id":"45dce12e-74a1-4d26-993d-0b5994e761d0","arxiv_id":"2603.00281","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In fire ant excavation, the number of active diggers grows as the square root of group size, and a crowding-based model explains the scaling.","lead":"Fire ants digging in narrow tunnels divide work unevenly, and the unevenness grows with group size: the number of ants doing the digging grows roughly like the square root of the total group size. A model where ants stop trying when tunnels are crowded reproduces the pattern and offers a mechanism for this scaling.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The theoretical n is instantaneous tunnel occupancy; the experimental n is a time-integrated active-worker count inferred from the Gini coefficient. The paper asserts but never validates that these are the same quantity, and the central scaling claim depends on it.","rationale":"The reader's weakest assumption is exactly the one I identify: the rate-equation variable n (instantaneous tunnel occupancy) is conflated with the experimentally inferred active-worker count obtained from the Gini coefficient. This is load-bearing because the analytic derivation is otherwise clean: the quadratic blockage probability and the steady-state condition do produce √N scaling for a fixed L, and the CA provides a plausible microscopic mechanism. But the experimental evidence for √N enters only through the transformation G → n with W=1; if that transformation does not correspond to the model's n, the mechanistic explanation is not connected to the measured scaling. The paper itself flags the assumption in a single sentence ('we assume that the number of ants in the tunnel at any given time, n, corresponds to the number of active ants in the collective') but provides no validation. The CA comparison does not help because it compares Gini coefficients, not occupancy versus active-worker counts. A direct re-analysis of the per-ant tracking data—which already exists for trials with ≤10 ants—could settle this. The verdict remains CONDITIONAL: the idea is plausible and valuable, but the central quantity mapping needs empirical support before the claim is accepted as stated.","tokens_in":18266,"tokens_out":8143,"duration_ms":91662,"concrete_test":"From the raw per-ant pellet counts and position tracks, compute for each trial: (a) the actual number of ants that deposited at least one pellet during the first 8 h (direct active count), (b) the mean and maximum instantaneous number of ants inside the tunnel, and (c) the Gini-inferred n = N(1−G) used in Fig. 10B. Plot all three vs N. Also redo the inference without forcing W=1, using the measured fraction of work done by the top n ants. If direct active count and instantaneous occupancy differ materially, or if the two obey different scaling exponents, the rate-equation explanation does not account for the measured inequality scaling.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that measured n ∝ √N is explained by the rate equation dn/dt = c1N − c2 n²/(2L)—hinges on identifying the model's n with the 'active ants' inferred from the Lorenz-curve simplification. In §V.A the model n is explicitly 'the number of ants in the tunnel at any given time.' The experimental n in Fig. 10B, however, is not a direct count of tunnel occupants; it is obtained by inverting G = W − n/N with W=1, i.e., n = N(1−G), where G summarizes the distribution of total grains deposited over 8 h. These are different observables: many distinct ants can take turns digging over 8 h, so the time-integrated set of contributors can be much larger than the instantaneous occupancy that enters the blockage probability. If ants rotate, the model's equilibrium occupancy could obey √N while the measured 'active' count obeys a different scaling (or vice versa). The paper asserts the correspondence in one sentence without validation, and the CA's Gini comparison does not resolve it because the CA also calculates Gini from cumulative per-ant pellet counts, not from instantaneous occupancy.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports experimental measurements of grain-deposition activity in Solenopsis invicta groups of size 2–25. From per-ant pellet counts the authors construct Lorenz curves and Gini coefficients and observe that workload inequality increases with group size. Using a two-segment Lorenz-curve approximation, G = W − n/N, and assuming W = 1, they convert measured G into an inferred number of \"active\" ants n and claim n ∝ √N. They then present a cellular-automaton model with crowding-modulated entry probability and an analytic mean-field rate equation in which the blockage probability scales as n²/(4L²), yielding steady-state n ∝ √N. The paper concludes that local crowding cues mechanistically explain the scaling.","tokens_in":18640,"tokens_out":11485,"duration_ms":104301,"significance":"If the central claims hold, the paper would provide a valuable mechanistic explanation for a scaling law that has been observed in other social systems. The strengths include direct experimental variation of group size, a transparent rate-equation derivation that is not circular, and an explicit CA parameter table. However, the central quantitative claims need additional support: the measured scaling is not fitted, the model variable n is identified with a different experimental observable without validation, and the CA robustness claim is stronger than the one-parameter-set comparison presented.","major_comments":[{"comment":"The model variable n in Eq. (6) is defined as \"the number of ants in the tunnel at any given time,\" but the experimental n plotted in Fig. 10B is obtained from Eq. (1) as n = N(1−G), where G is computed from cumulative 8-h pellet counts. These are different observables: ants that take turns digging can contribute over 8 h without ever being simultaneously present in the tunnel. The CA does not resolve this mismatch because its Gini coefficients are also computed from cumulative per-ant pellet counts, not from instantaneous occupancy. Please validate the identification directly, e.g., by measuring instantaneous tunnel occupancy from the tracking data and testing its scaling, or by modeling the cumulative contributor count; otherwise the analytic derivation does not explain the measured scaling.","section":"§V.A, Fig. 10B"},{"comment":"The headline scaling n ∝ √N is asserted without a regression or uncertainty estimate. With only eight group sizes and 3–4 trials each, a visual comparison to a √N curve is insufficient. Fit log n vs. log N, report the exponent and confidence interval, and assess goodness of fit. In addition, n is inferred under the assumption W = 1; because n = N(W−G), the inferred values and their scaling are sensitive to W. Fit W from the two-segment Lorenz curves or report a sensitivity analysis.","section":"Fig. 10B, §III"},{"comment":"Equation (2) uses c1N as the entry rate while n ants are already in the tunnel. The outside pool is N−n, so the entry rate should be c1(N−n). With this correction the steady-state relation is c1(N−n) = c2 n²/(2L), which reduces to n ∝ √N only asymptotically and has finite-size corrections that are not negligible over the experimental range N = 2–25. Please revise the derivation or explicitly justify the approximation.","section":"Eqs. (2)–(6), §V.A"},{"comment":"The abstract states that the CA \"reproduces experimental Gini coefficients over a wide range of parameters,\" but Fig. 6 shows agreement only for one parameter set (R = 0.1, τ0 = 10, C0 = 1, τs = 2), and Fig. 9 shows strong dependence on R and τs, with only the chosen values matching the experimental trend. Please quantify the parameter region in which the simulated Gini curves are statistically consistent with the data and adjust the claim accordingly.","section":"Abstract, §IV.B, Figs. 6 and 9"}],"minor_comments":[{"comment":"Typo: \"We find that that workload becomes...\" should read \"We find that workload becomes...\".","section":"Abstract"},{"comment":"The derivative of D = n(1 − n²/(2L)) is dD/dn = 1 − 3n²/(2L). The text writes 1 − c n²/(2L) without defining c; if c = 3 is intended, say so explicitly.","section":"§V.B, Eq. (9)"},{"comment":"The automated tracking method is validated only for groups of 10 or fewer; for groups of 15–25 ants activity was tracked manually. This methodological discontinuity could introduce systematic differences in the larger-group Gini values and should be acknowledged in the error analysis.","section":"§VII.A"},{"comment":"The caption reports simulation variability as standard deviation over 10 trials, but the number of experimental trials varies (4 trials for 2–6 ants, 3 trials for 10–25 ants). Please state these numbers in the caption and consider showing individual trial points for the larger groups.","section":"Fig. 6"}],"recommendation":"major_revision","confidential_remarks":"The central concern—the identification of the analytic model's instantaneous tunnel occupancy with the Lorenz-inferred, time-integrated active-worker count—is real and load-bearing. I do not view it as unfixable: the authors have tracking data that could directly test the occupancy scaling, and the model could be extended to predict cumulative contributors. The missing regression and the overstatement of CA parameter robustness are also fixable. I would not reject, but the manuscript should not be accepted in its present form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Main thing to know: this is a genuine empirical result worth taking seriously, but the paper's central mechanistic claim is not yet established. The new observation is that Gini inequality increases with group size in fire ant excavation, and via a simplified Lorenz-curve assumption the inferred number of 'active' ants scales roughly as sqrt(N). The analytic rate-equation explanation is simple and attractive: if two ants abreast cause a blockage, failure rate ~ n^2, so balance c1*N against c2*n^2 gives n ∝ sqrt(N). That derivation is fine on its own terms, and tying it to a measured scaling in a biological collective is genuinely new.\n\nWhat the paper does well: systematic group-size variation (2–25), a real video-tracking innovation, explicit connection to prior excavation work, and a CA model that qualitatively captures the trend with local crowding rules. The CA is not a free-for-all: most parameters come from earlier studies, and only a few are new.\n\nWhere it gets soft. The stress-test concern is real. In V.A, the model n is defined as instantaneous tunnel occupancy. The experimental n in Fig. 10B is obtained from G = W - n/N with W = 1, so it counts the number of ants responsible for essentially all work over the 8-hour trial. Those are different observables unless ants do not take turns. The paper simply asserts the correspondence in one sentence. This matters because the whole point is to give a mechanism for the measured scaling, not just a mechanism for occupancy. The CA does not resolve it either: it computes Gini from cumulative per-ant pellet counts, not from instantaneous occupancy. So the explanation of the Gini scaling is conditional on an unvalidated identification.\n\nSecondary issues: the 'wide range of parameters' claim is too strong. Fig. 9 shows Gini trends depend on R and tau_s; the match is essentially one hand-picked set (R = 0.1, tau_s = 2). That is not disqualifying, but it should be described as a fit, not a robust sweep. Also, the sqrt(N) scaling is asserted from a plot; no regression or exponent uncertainty is reported. And data/code are only 'available upon request,' which is weak for a scaling claim.\n\nOverall, the empirical pattern may well be right and the mechanism is plausible. But the load-bearing equivalence between the two n's needs direct validation — e.g., measuring instantaneous occupancy in the experiments, or deriving Gini from the analytic model instead of assuming W = 1. That is a major revision, not a desk reject.\n\nRecommendation: send it to a serious referee. A good referee will push on the n/n distinction and the parameter selection. If the authors close that gap, this could be a useful, citable paper. I would want to see the revised version before citing it myself.","headline":"Genuine new scaling and a clean mechanism, but the model's n and the measured n are different observables — that gap needs closing before the claim holds.","tokens_in":19057,"tokens_out":3914,"would_cite":false,"duration_ms":41310,"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":"In groups of 2 to 25 fire ants digging a narrow tunnel, the number of ants doing nearly all the work grows as the square root of the group size, a scaling the paper derives from two-ant blockages in the tunnel.","keywords":["fire ants","workload inequality","Gini coefficient","Lorenz curve","square-root scaling","collective excavation","rate equation","cellular automata"],"falsifier":"Track instantaneous tunnel occupancy and per-ant pellet counts in the same 2-to-25-ant trials: if the occupancy count does not rise as √N while the pellet-based active set does, the rate-equation identification is wrong. As a complementary test, widen the tunnel to three ant widths; the quadratic-blockage mechanism predicts the √N scaling should weaken or change exponent.","tokens_in":18219,"feed_emoji":"🐜","tokens_out":7115,"duration_ms":65702,"temperature":0.7,"pith_summary":"This paper asks why work is shared unequally in collective excavation and proposes a mechanism. In experiments with 2 to 25 fire ants digging a two-ant-wide tunnel, the Gini coefficient of workload rises with group size, and the inferred number of ants that do essentially all the digging scales as the square root of the total group size. A cellular-automaton model reproduces the trend when simulated ants reduce their entry probability after being blocked, and a rate equation shows why: in a tunnel of length L with n ants inside, the chance of a traversal-ending blockage is roughly n²/(2L), so at steady state the active count satisfies n ∝ √N. The paper presents this as a mechanistic origin for a scaling law previously seen only observationally in other social systems, and it suggests that local crowding cues, rather than intrinsic differences among ants, create the inactive majority in larger groups.","feed_headline":"Active fire ants scale as square root of group size","feed_subtitle":"Two-ant-wide tunnel blockages, scaling as the square of active ants, turn a social-science law into a physical one.","key_machinery":"The rate equation dn/dt = c1N − c2 n²/(2L), where n is the number of ants in the tunnel, N the group size, L the tunnel length, and c1, c2 constants. The quadratic term is the load-bearing object: it is the probability that a two-ant-wide tunnel is blocked during a traversal, scaling as n²/(4L²) per site and ~n²/(2L) per trip. Setting the flux to zero yields n ∝ √N, and the same quadratic in the digging-rate expression D = n(1 − n²/2L) yields the companion optimum n ∝ √L. The mechanism converts local crowding into a self-limiting participation rate, and it is the piece the cellular automaton needs in order to match the measured inequality.","core_discovery":"The paper's central claim is that the number of ants n that carry out the work scales with the square root of the total group size N, and that this scaling is set by the tunnel geometry. The corridor is only two ant body widths wide; a blockage at a given position occurs when two ants occupy the two cells across the width. With n ants randomly distributed along a tunnel of length L, the probability of a blockage at one site is ~n²/4L², making the probability that a crossing fails ~n²/2L. Modeling the active population as dn/dt = c1N − c2 n²/(2L) — ants entering at a rate set by group size and leaving when blocked — and setting the derivative to zero gives n ∝ √N. The same quadratic term appe","pith_inferences":["Directly measuring instantaneous tunnel occupancy n(t) alongside per-ant pellet totals would test the identification that carries the derivation; if occupancy does not itself scale as √N, the rate equation is not explaining the measured active count.","Varying tunnel width offers a quantitative extension: with a w-ant-wide tunnel, the same occupancy argument predicts a blockage probability ~n^w/(wL)^{w−1} and hence a different scaling exponent, a prediction testable with ants or robots.","If the collision-cue mechanism is general, the same square-root law should appear in any controlled task where agents queue through a fixed-width bottleneck — a prediction that could be checked in human crowds or engineered swarms."],"forward_implications":["Workload inequality in constrained excavation is a by-product of local, geometry-driven crowding: the same rule that makes an ant leave a jammed tunnel produces the global √N scaling.","The experimental Gini coefficients continue rising sublinearly with group size, and the trend extrapolates to the Gini range previously reported for 30-ant groups.","Because the digging rate is maximized when n ∝ √L, and experiment gives n ∝ √N, ants would be excavating at optimal speed whenever tunnel length is proportional to group size — a correlation observed in natural nests.","The mechanism generalizes: any collective in which individuals leave a confined workspace when blocked should show a square-root active subset, which the authors offer as a possible reason similar power laws appear in other biological and social systems."],"fun_headline_variants":["Square-root law for ant workload emerges from tunnel blockages","Why bigger ant groups dig less equally: tunnel width matters","Fire ant work inequality scales as √N, and here's why","Crowding in narrow tunnels explains square-root scaling of ant labor","Gini coefficient for ant digging rises with group size, mechanistically"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The derivation identifies n, the number of ants in the tunnel at any instant, with the n inferred from the final Lorenz curve as the ants responsible for nearly all work; if instantaneous occupancy and the time-integrated active subset are not the same population, the rate equation does not explain the measured scaling.","fun_headline_variants_meta":{"raw":{"variants":["Square-root law for ant workload emerges from tunnel blockages","Why bigger ant groups dig less equally: tunnel width matters","Fire ant work inequality scales as √N, and here's why","Crowding in narrow tunnels explains square-root scaling of ant labor","Gini coefficient for ant digging rises with group size, mechanistically"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001285,"raw_usage":{"total_tokens":5126,"prompt_tokens":819,"completion_tokens":4307,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":563,"completion_tokens_details":{"reasoning_tokens":4220}},"tokens_in":563,"tokens_out":4307,"duration_ms":30355,"temperature":1.0,"reasoning_tokens":4220,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T19:58:29.579729+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track instantaneous tunnel occupancy and per-ant pellet counts in the same 2-to-25-ant trials: if the occupancy count does not rise as √N while the pellet-based active set does, the rate-equation identification is wrong. As a complementary test, widen the tunnel to three ant widths; the quadratic-blockage mechanism predicts the √N scaling should weaken or change exponent.","supporting_citations":[],"review_version":1}