{"id":"93869878-2f23-455f-bad6-00dfb98161bf","arxiv_id":"2506.17653","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In simulated early biofilms, tracer particles show normal diffusion at short times and growth-driven exponential displacement at long times, with escape only for very small tracers in loose colonies.","lead":"This paper simulates how a small tracer particle moves inside a growing bacterial biofilm. It finds that at short times the tracer diffuses normally, while at long times the colony's expansion drags the particle outward, producing motion that grows exponentially with time.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Release-time threshold mstart=100 is explicitly acknowledged as affecting tracer dynamics but never varied; the two-regime picture may be specific to this single release point.","rationale":"The reader's weakest_assumption identifies exactly the concern I find most load-bearing: the untested dependence on mstart, explicitly acknowledged in Section II. The central claim—long-time MSD linear in biomass and independent of tracer size, driven by colony expansion—would have to be robust to release time for it to describe early biofilm diffusion generally. The paper provides no sensitivity analysis for mstart, and the escape data in Fig. 6 show that at Γ=0.1 the smallest tracers are already marginal at mstart=100, so the regime boundaries are plausibly sensitive to this threshold. I considered other potential concerns: the trajectory-conditioning rule (excluding escaped tracers) could bias survivor averages, and the absence of error bars makes collapse claims harder to judge, but these are secondary or would affect interpretation rather than overturn the central mechanism. The kinematic argument for exponential MSD from uniform colony expansion is physically reasonable, and the simulation data appear to support it within the chosen protocol, so the appropriate response is to keep the paper conditional pending the mstart sensitivity check, not to reject it. Since the reader already reached CONDITIONAL for this reason, my read does not change the verdict.","tokens_in":15903,"tokens_out":7604,"duration_ms":87468,"concrete_test":"Repeat the BD protocol for representative cases—e.g., Γ=10 and Γ=0.1 with dt=1σ and dt=10^-2σ—using mstart=25, 50, 200, and 400 while keeping all other parameters and the 1000-trajectory ensemble fixed. Compare (i) the short-time effective diffusion coefficient extracted as in Fig. 4, (ii) the long-time slope d log⟨Δr²⟩/d log m over m∈[1000,4000], and (iii) the escape probability p(t) at fixed biomass values. If the short-time Deff shifts by more than ~20% or if the long-time slope or escape threshold changes materially with mstart, the two-regime description is specific to the chosen release threshold and the paper's generalizations should be explicitly re-scoped.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central two-regime description—short-time Brownian diffusion modulated by colony ageing, followed by long-time growth-induced advection with tracer-size-independent MSD—is established only for tracers released when the colony reaches biomass mstart=100. Section II states this directly: 'We highlight that the tracer’s dynamics are expected to be affected by the value of mstart, as it experiences different pressure environments depending on the colony’s growth stage. While this is an important consideration, a thorough analysis on the influence of mstart on the tracer’s dynamics is beyond the scope of this work.' This is a load-bearing modeling choice, not a peripheral detail. The tracer is held fixed while the colony grows around it, so mstart determines the local packing fraction, pressure, and cavity structure in which the tracer begins to move. Changing mstart could shift the short-time effective diffusion coefficient (already shown in Fig. 4 to depend on colony age via t0) and could change whether small tracers in loose colonies are trapped long enough to exhibit the claimed exponential radial MSD or escape before advection dominates. The fragility is visible in the authors' own Fig. 6: for Γ=0.1 and dt≤10^-2σ, the probability of remaining in the colony falls to zero almost immediately after release at mstart=100, so the long-time 'drag' result for those tracers rests on a rapidly shrinking survivor ensemble. If a different mstart moved the escape threshold or altered the long-time slope, the claimed universality and the biological generalization to 'early growing biofilms' would not hold beyond this arbitrary release condition. The concern is not that the simulations are wrong, but that the headline claim has an acknowledged, untested parameter dependence at its boundary.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports Brownian dynamics simulations of a passive spherical tracer inside a growing two-dimensional bacterial microcolony, using an individual-based model calibrated to Pseudomonas putida. Tracers of four diameters (dt = 10^-3σ to σ) are held fixed until the colony biomass reaches mstart = 100 and then followed in three growth regimes Γ = 10, 1, 0.1. The main findings are: (i) at short times the tracer MSD is linear in time, with an effective diffusion coefficient that depends on tracer size and weakly on colony age; (ii) at long times the radial MSD grows exponentially in time and linearly in biomass, independently of tracer size for compact colonies, while the tangential MSD is subdiffusive; (iii) small tracers in loose colonies (Γ = 0.1) escape and show a different behavior; and (iv) passive microrheology via the Mason relation yields a mostly viscous-like response, with an elastic regime for small tracers in loose colonies. The central mechanistic claim is that colony expansion creates an outward radial drag that dominates long-time transport, with Brownian motion contributing mainly in the tangential direction.","tokens_in":16194,"tokens_out":5827,"duration_ms":61130,"significance":"If the central result holds, the paper makes a useful contribution by identifying growth-induced advection, rather than thermal diffusion, as the dominant long-time transport mechanism for tracers in early biofilms, with implications for nanoparticle and antibiotic transport and for interpreting microrheology in active growing cell communities. Strengths include the direct simulation of 1000 independent trajectories per case, the use of an established agent-based model with parameters fixed in earlier work, and an explicit statement of the model's limitations. The exponential/linear-in-biomass scaling and the exception for small tracers in loose colonies are concrete, falsifiable predictions that could be tested experimentally. However, the main quantitative claims are currently supported mostly by visual collapses without fitted parameters or error bars, and they depend on a single release threshold mstart = 100 whose influence is acknowledged but not tested.","major_comments":[{"comment":"The text explicitly states that tracer dynamics are expected to depend on mstart, yet no mstart variation is reported. Because mstart sets the local packing fraction, pressure, and cavity connectivity at release, it directly controls both the short-time effective diffusion coefficient and whether small tracers in loose colonies escape before the long-time advective regime develops. I request a sensitivity analysis (e.g., mstart = 50, 200, 500) for at least one compact and one loose growth regime and two tracer sizes, reporting Deff and the long-time scaling parameters; without this, the universality of the two-regime description is not established.","section":"Section II, mstart = 100"},{"comment":"The MSD averages are conditioned on the tracer remaining inside the colony at every time step, and Fig. 6 shows that for Γ = 0.1 and dt ≤ 10^-2σ the survival probability p(t) drops to zero almost immediately after release. The long-time MSD for these cases is therefore computed over a rapidly shrinking and increasingly unrepresentative subensemble. Please report the number of surviving trajectories as a function of m(t) or t - t0, provide confidence intervals or bootstrap errors on the MSD, and quantify the conditioning bias, e.g., by comparing with an average that includes trajectories up to their first escape. This is essential because the exponential scaling claim rests on the behavior of the survivors.","section":"Section III.A and Figs. 2-3"},{"comment":"The central quantitative claim is that at long times ⟨Δr²⟩ ∝ m(t) ∝ exp(kt), and that this behavior is independent of tracer size. The evidence currently is the visual alignment of data with the unlabelled dashed lines '~ b m(t)' and '~ e^{kt}'. No fitted values of b or k, no confidence intervals, and no comparison with the independently known growth rate r are given. Please fit the long-time regime for each Γ and tracer size, report the fitted exponents or rates with uncertainties, and define k in terms of the model parameters. This would turn the collapse claim from a visual statement into a quantitative result.","section":"Fig. 2 and its inset"},{"comment":"The tracer-bacterium interaction strength ε_tr = 10 kBT is a free parameter, and no sensitivity analysis or physical justification is provided for this value. Since the escape probability and the short-time Deff both depend on the steric repulsion between the tracer and surrounding rods, a variation of ε_tr (e.g., 1 kBT and 100 kBT) for a representative case would help establish that the two-regime picture does not hinge on this choice.","section":"Section II, Eq. (1)"},{"comment":"The generalized Stokes-Einstein relation used to compute G' and G'' assumes equilibrium linear response and a passive thermal tracer. In the growing colony, the tracer's long-time radial motion is dominated by non-equilibrium advection, so the interpretation of Fig. 5 as the viscoelastic moduli of the biofilm is conceptually strained. The authors do acknowledge that the analysis is qualitative, but the text still refers to 'the elastic and viscous moduli of the growing microcolony' without discussing how the advective contribution is separated from the thermal part in Eq. (14). Please either restrict the MR interpretation to time or frequency windows where the tracer motion is diffusive and thermal, or add an explicit discussion of why the GSE remains approximately valid under growth-induced drag.","section":"Section III.B, Eqs. (14)-(16)"}],"minor_comments":[{"comment":"Several typos should be corrected: 'Manson' should be 'Mason' (Section III.B), 'Boltzmamn constant' should be 'Boltzmann constant', and 'unit ot time' should be 'unit of time' (Section II).","section":"Throughout"},{"comment":"The gamma function Γ[1 + α(ω)] uses the same symbol as the growth-regime parameter Γ; this notational collision should be removed, for example by using a different letter for the gamma function.","section":"Eq. (14)"},{"comment":"The horizontal axis is labeled 'dt', but the reader must infer that it is the tracer diameter normalized by σ; please label it explicitly as dt/σ and add units.","section":"Fig. 6, right panel"},{"comment":"References [22] and [56] are the same paper (Sinha, Samanta, and Thirumalai, Soft Matter 19, 5385 (2023)) and should be consolidated.","section":"References [22] and [56]"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely publishable after a revision that adds the mstart sensitivity analysis and statistical error reporting; the central mechanism is plausible and the simulation setup is transparent. I would not insist on a fully quantitative microrheology validation, but the current MR section needs a clearer caveat about the non-equilibrium contribution. The main risk to the long-time universality claim is survivor bias for Γ = 0.1 small tracers."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a credible, straightforward simulation paper. The new thing is the application of the group's IbM model to tracer diffusion in growing bacterial colonies, and the central result is a clean two-regime picture: short-time Brownian motion that slows as the colony ages, then long-time advection by colony expansion, with the radial MSD growing exponentially in time (linearly with biomass). The collapse of the long-time MSD onto a universal biomass curve across tracer sizes is a genuinely new observation, and the size-dependent escape probabilities for loose colonies (Γ=0.1) are biologically useful. The simulations are heavy: 1000 trajectories per condition, and the authors are careful to check seeding geometry.\n\nThe main soft spot is exactly what the stress-test flags: the tracer is held fixed until the colony reaches biomass mstart=100, and the authors explicitly say in Section II that tracer dynamics are expected to depend on mstart, but they never vary it. That is load-bearing, because mstart sets the local packing fraction and the initial cavity structure around the tracer. If the long-time scaling or the escape thresholds shift with mstart, the paper's generalization to \"early growing biofilms\" is only a specific-case result. Adding even two or three mstart values, or a short paragraph of sensitivity data, would fix this.\n\nThe other weaknesses are real but secondary: no error bars, no code/data release, and no sensitivity analysis for the tracer-rod interaction strength (ε_tr=10 kBT). The microrheology section is explicitly labeled qualitative, so I take it as illustrative. The trajectory-conditioning rule—averaging only trajectories where the tracer stays in the colony—could bias the long-time MSD for escaping tracers, but the escape probability curves in Fig. 6 help put this in context.\n\nThe central argument does hold up: the exponential MSD scaling is an observed simulation result, not a fitted curve, and the model parameters are inherited from prior work. The uniform-density drag explanation is post hoc, but it's a reasonable interpretation, not a circular confirmation.\n\nThis paper deserves a serious referee. It's a solid contribution to biofilm microrheology and nanoparticle transport in growing cell colonies. I'd bring it to a reading group to discuss the release-time issue and the advection mechanism. I'd cite it with a caveat about mstart if I worked in that area. The right recommendation for an editor is to send it out, with a request for a sensitivity analysis on mstart before acceptance.","headline":"Solid simulation study giving a clear two-regime picture of tracer transport in growing biofilms, but the load-bearing mstart=100 release threshold is acknowledged yet untested.","tokens_in":16756,"tokens_out":3054,"would_cite":true,"duration_ms":28256,"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":"Colony expansion, not thermal motion, controls long-time tracer drift in young biofilms.","keywords":["tracer diffusion","biofilm","Brownian dynamics","individual-based model","microrheology","mean squared displacement","subdiffusion","colony growth"],"falsifier":"Run the same simulation with the tracer released at a much smaller or much larger colony size; if the claimed linear-in-biomass scaling does not persist, the release threshold is load-bearing. In the lab, track a tracer in a confined two-dimensional growing colony and decompose its motion into radial and tangential parts; the claim fails if the radial mean squared displacement grows as a power law rather than exponentially at long times.","tokens_in":15731,"feed_emoji":"🦠","tokens_out":7834,"duration_ms":73842,"temperature":0.7,"pith_summary":"This paper tries to establish a two-regime law for how a passive tracer moves inside a growing bacterial colony during its early, quasi-two-dimensional stage. On short time scales the tracer diffuses normally through the spaces between bacteria, with a diffusion coefficient that decreases as the colony ages and compacts. On long time scales, the colony's exponential biomass growth creates an outward drag field, so the tracer's mean squared displacement grows exponentially in time and linearly in biomass, and this long-time curve is independent of tracer size within each growth regime. The notable exception is very small tracers in loose, slowly growing colonies, which slip through the structure and escape. If true, this gives a predictive transport rule for nanoparticles, antibiotics, and viruses in early biofilms and a way to read the mechanical feel of a growing colony from tracer motion.","feed_headline":"Colony growth, not heat, drags tracers in young biofilms","feed_subtitle":"Simulations show long-time tracer motion scales with biomass and is size-independent in compact colonies.","key_machinery":"The load-bearing mechanism is the exponential-growth drag field: if a colony's area obeys $dA/dt = kA$, then every subregion expands by the same law, producing an outward velocity $dR/dt = kR/2$ at radial distance $R$ from the colony center. That single relation converts exponential biomass growth into the prediction $\\langle \\Delta r^2 \\rangle \\propto \\exp(kt) \\propto m(t)$, which the simulations confirm. The computational machinery is an individual-based model (IbM) of spherocylinder bacteria that grow and divide, a spherical tracer interacting through a truncated Kihara potential, Brownian dynamics integration, and the generalized Stokes–Einstein relation used to convert MSD curves into storage and loss moduli.","core_discovery":"The central claim is that in early-stage biofilms, growth-induced advection replaces thermal diffusion as the dominant long-time transport mechanism for embedded tracers. Because biomass grows as $m(t) \\simeq \\exp(rt)$, the colony area grows at a rate proportional to itself, and any interior point experiences an outward radial velocity $dR/dt = kR/2$. This drag makes the long-time mean squared displacement (MSD) scale as $\\langle \\Delta r^2 \\rangle \\propto \\exp(kt) \\propto m(t)$, with the same curve for tracers differing by three orders of magnitude in size within a given growth regime. The radial component of the MSD is exponential while the tangential component is subdiffusive, showing that the tracer is carried outward by colony expansion while Brownian collisions drive sideways motion. Only for the smallest tracers ($d_t = 10^{-2}\\sigma$ and $10^{-3}\\sigma$) in the most open colonies ($\\Gamma = 0.1$) does the tracer escape the drag and behave as a Brownian particle in an elastic structure. The same MSD data, converted through the generalized Stokes–Einstein relation, yield elastic and viscous moduli: the growing colony reads mostly as a viscous fluid, except for small tracers in loose colonies at low frequency, which see an elastic response.","pith_inferences":["I infer that re-plotting existing tracer data from growing tissues as log MSD against biomass, not time, would separate growth drag from other active transport: a straight line would indicate the same mechanism.","I infer a concrete simulation test of the untested release-time choice: vary $m_{\\mathrm{start}}$ and check whether the short-time diffusion coefficient and the exponential prefactor shift; if they do, the two-regime law needs a release-time correction.","I infer that adding the extracellular polymeric substance (EPS) matrix, which the model omits, would slow short-time diffusion and likely suppress escape of small tracers in loose colonies, changing the escape thresholds.","I infer that because the drag velocity grows with distance from the colony center, a tracer starting off-center should escape more readily than the centrally seeded tracers studied here; off-center seeding is a direct way to test the geometry of the drag field."],"forward_implications":["At long times, tracer displacement in a growing colony is set by the biomass growth rate, so transport estimates for antibiotics or nanoparticles in early biofilms can use the colony's growth kinetics rather than only local viscosity.","In compact colonies ($\\Gamma = 10$ and $1$), the long-time MSD is essentially independent of tracer size, so one advection law covers tracers from nanometers to the size of a bacterium.","Very small tracers in loose colonies ($\\Gamma = 0.1$, $d_t \\le 10^{-2}\\sigma$) escape the colony, so size-based filtering of nanoparticles is strongest in open, slowly growing microcolonies.","Passive microrheology of the growing colony reports a predominantly viscous response, with an elastic response only for small tracers in loose colonies at low frequencies; because the underlying equilibrium assumptions may not hold in an active growing system, these moduli should be read qualitatively.","The mechanism implies that reported long-time 'superdiffusive' tracer motion in growing cell collectives such as cancer spheroids may actually be exponential in time; checking the functional form would tell whether the same growth-drag mechanism is at work."],"supporting_citations":[{"why":"Supplies the individual-based model of biofilm growth and the closed/open growth classification that defines the Γ regimes used throughout.","marker":"[40]"},{"why":"Establishes the exponential biomass growth law and the Γ ratio that sets colony compactness, the basis for the drag-field prediction.","marker":"[42]"},{"why":"Provides the passive-microrheology relation connecting tracer MSD to frequency-dependent moduli.","marker":"[8]"},{"why":"Gives the generalized Stokes–Einstein expressions (Eqs. 14–16) used to compute G′ and G″ from the MSD.","marker":"[9]"},{"why":"Provides the minimum-distance algorithm between rods used to evaluate the Kihara interactions in the simulations.","marker":"[48]"},{"why":"Supplies the size-dependent translational and rotational diffusion coefficients for growing spherocylinder bacteria.","marker":"[51]"},{"why":"Experimental microrheology of bacterial biofilms that reports both subdiffusive and superdiffusive behavior, the context the early-growth results are meant to clarify.","marker":"[24]"},{"why":"Reports long-time superdiffusive tracer motion in growing cancer spheroids, the comparison case the paper suggests may share the exponential-drag mechanism.","marker":"[56]"}],"fun_headline_variants":["Growth, not heat, drives tracer motion in biofilms","Biofilm expansion carries tracers like a conveyor belt","In young biofilms, growth powers tracer drift","Colony growth overwhelms thermal diffusion for tracers","Tracers in biofilms are dragged by colony expansion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The tracer is held still until the colony has grown to a fixed size, and the study does not test how that starting point changes the result; if either the short-time diffusion coefficient or the long-time exponential scaling shifts with release time, the two-regime law may be specific to that threshold.","fun_headline_variants_meta":{"raw":{"variants":["Growth, not heat, drives tracer motion in biofilms","Biofilm expansion carries tracers like a conveyor belt","In young biofilms, growth powers tracer drift","Colony growth overwhelms thermal diffusion for tracers","Tracers in biofilms are dragged by colony expansion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000239,"raw_usage":{"total_tokens":1545,"prompt_tokens":1007,"completion_tokens":538,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":623,"completion_tokens_details":{"reasoning_tokens":461}},"tokens_in":623,"tokens_out":538,"duration_ms":5306,"temperature":1.0,"reasoning_tokens":461,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:04:28.154654+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same simulation with the tracer released at a much smaller or much larger colony size; if the claimed linear-in-biomass scaling does not persist, the release threshold is load-bearing. In the lab, track a tracer in a confined two-dimensional growing colony and decompose its motion into radial and tangential parts; the claim fails if the radial mean squared displacement grows as a power law rather than exponentially at long times.","supporting_citations":[{"cited_title":"Volfson, S","cited_arxiv_id":null,"evidence_quote":"Supplies the individual-based model of biofilm growth and the closed/open growth classification that defines the Γ regimes used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the exponential biomass growth law and the Γ ratio that sets colony compactness, the basis for the drag-field prediction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the passive-microrheology relation connecting tracer MSD to frequency-dependent moduli."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the generalized Stokes–Einstein expressions (Eqs. 14–16) used to compute G′ and G″ from the MSD."},{"cited_title":"Taheri-Araghi, S","cited_arxiv_id":null,"evidence_quote":"Provides the minimum-distance algorithm between rods used to evaluate the Kihara interactions in the simulations."},{"cited_title":"Vega and S","cited_arxiv_id":null,"evidence_quote":"Supplies the size-dependent translational and rotational diffusion coefficients for growing spherocylinder bacteria."},{"cited_title":"Weihs, T","cited_arxiv_id":null,"evidence_quote":"Experimental microrheology of bacterial biofilms that reports both subdiffusive and superdiffusive behavior, the context the early-growth results are meant to clarify."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports long-time superdiffusive tracer motion in growing cancer spheroids, the comparison case the paper suggests may share the exponential-drag mechanism."}],"review_version":2}