{"id":"1d7c7b83-7bb0-4686-b876-0f65a05b2375","arxiv_id":"2411.13359","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In DPD simulations, sedimented squirmer suspensions form hexagonal monolayers whose defect annealing is faster for pullers than for pushers.","lead":"Simulations of sedimenting microswimmer suspensions show that as gravity increases, the bottom layer forms a hexagonal crystal, and active swimmers, especially puller-type, fix defects in that crystal faster than passive colloids do. The result helps understand how swimming style could control biofilm formation and particle settling in bioreactors.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The puller/pusher ordering comparison is not isolated from the raspberry shape: high-gravity orientation peaks are attributed to the polyhedral facets (SM fig. 7), yet the claim that shape effects are 'the same in both' is asserted without a spherical-colloid control.","rationale":"I considered the reader's identified weakest assumption and found it to be the most load-bearing concern. The paper's symmetry argument ('the effects ... are the same in both') is valid for the static shape potential, but it does not cover the hydrodynamic force application, which is mediated by the same discrete filler particles. Since the sign of B2 changes the solvent force field, the reaction forces sampled at the polyhedral filler positions can produce β-dependent torques; this is a concrete mechanism by which the shape could couple differently to pullers and pushers. The paper's own attribution of the high-gravity P(α) peaks to the raspberry structure (SM fig. 7, §3.3) shows that shape effects are not negligible in the regime where the ordering claim is made. A spherical-colloid control, which the authors state is achievable, would settle this. I also note secondary concerns (no error bars on P(ψ6)/G6(r), and the un-matched first-layer packing fraction between pullers and pushers, ~0.82 vs ~0.78), but these would weaken the statistical or mechanistic interpretation without invalidating the comparison as fundamentally as a shape bias would. Therefore the reader's CONDITIONAL verdict is appropriate; requiring the spherical-shape check as a condition of acceptance is consistent with UNCHANGED.","tokens_in":24008,"tokens_out":9755,"duration_ms":110031,"concrete_test":"Re-run the bottom-layer simulations at Fg/Fp = 1.5 and 2.25 using the spherical-colloid variant the paper itself describes in §2.1.2: switch off conservative interactions of the 18 filler particles and instead give the central thruster a conservative DPD cutoff large enough to reproduce the same effective sphere (or use a colloid with 33 or 66 surface fillers). Recompute P(ψ6), G6(r), and P(α) for pullers and pushers. If the puller/pusher difference in G6(r) at high gravity persists, the conclusion is robust; if it vanishes or reverses, the reported 'pullers preserve order better' is a raspberry-shape artifact rather than a stresslet-sign effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that pullers preserve hexagonal order better than pushers rests on the assumption, stated in §3.2, that 'the effects due to weak departure from sphericity are the same in both' swimmers, so differences are attributable to hydrodynamics. This assumption is not tested and is contradicted in spirit by the paper's own use of SM fig. 7: the 19-bead raspberry has preferred facet normals near 30.4°, 54.7°, and 69.1° (with supplements at 149.6°, 125.3°, and 110.9°), which align closely with the observed P(α) peaks, ≈125° and ≈155° for pullers, ≈125° for pushers, and the ≈60° and ≈125° peaks at high gravity. Because the self-propulsion force in eq. (2) is exerted on solvent in a spherical shell while the reaction is applied to the nearest filler particle (eq. 4), the discrete polyhedral surface determines the local torque on the colloid. The B2 term changes sign between pullers and pushers, so the reaction-force distribution over the faceted surface is not guaranteed to be identical for the two swimming modes; a β-dependent effective torque or wall–colloid coupling would make the ordering difference a shape artifact rather than a property of stresslet sign. The paper's own §3.3 attributes the high-gravity orientation angles to the raspberry structure, further weakening the 'same in both' assertion. A spherical-colloid control is therefore required to establish the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript uses dissipative particle dynamics simulations of 19-bead raspberry squirmer colloids immersed in an explicit DPD solvent, under a gravitational force applied to the colloids, to study sedimentation and the structure of the sedimented bottom layer for passive colloids, pullers (β = +10), and pushers (β = −10). The passive sedimentation velocity is validated against Stokes law with a fitted effective radius, and the squirmer sedimentation profiles are compared with earlier MPCD and LB studies, reproducing known trends such as stronger wall accumulation for pullers and longer sedimentation lengths for pushers. The main new claim is that, as the gravitational field increases, the bottom layer undergoes a transition to a hexagonal crystal, that activity helps anneal defects relative to the kinetically trapped passive case, and that pullers preserve hexagonal order better than pushers, reaching an almost defect-free crystalline state at lower Fg. The structural analysis uses P(ψ6), G6(r), g(r), the static structure factor, polar order parameters, and orientation distributions P(α). The authors are appropriately cautious about the model's limitations, explicitly stating that lubrication forces are not included and that quantitative predictions are not the goal.","tokens_in":24242,"tokens_out":6432,"duration_ms":67710,"significance":"If the central claim survives scrutiny, the paper makes a useful contribution by demonstrating, within a single DPD model, that activity and specifically puller-type swimming accelerate the formation of nearly defect-free hexagonal monolayers in sedimented squirmer suspensions, and by connecting this to known wall-alignment hydrodynamics. The manuscript has real strengths: it validates the model against Stokes law and against independent MPCD/LB results, it describes the model and its parameters in detail, and it explicitly quantifies the expected size of lubrication effects. However, the load-bearing comparison between pullers and pushers currently rests on a single simulation protocol and on an untested assumption that the raspberry polyhedral shape affects both swimmers identically. The significance is therefore conditional on additional controls and statistical support.","major_comments":[{"comment":"The statement in §3.2 that 'the effects due to weak departure from sphericity are the same in both' is load-bearing for the central claim of §3.3 that pullers are better than pushers at preserving hexagonal order, but it is only asserted, not demonstrated. SM Fig. 7 shows that the 19-bead raspberry has preferred facet normals at 30.4°, 54.7°, and 69.1°, with supplements at 149.6°, 125.3°, and 110.9°, values that are suspiciously close to the observed P(α) peaks (≈125° and ≈155° for pullers, ≈125° for pushers, and ≈60° and ≈125° at high Fg). Because the propulsion reaction force is applied to the nearest filler particle (Eq. 4) and the B2 term changes sign between pullers and pushers, the discrete polyhedral surface can in principle produce a β-dependent effective wall torque or wall–colloid coupling. The paper's own §3.3 attributes the high-gravity orientation angles to 'the swimmers' raspberry structures', which further weakens the 'same in both' assumption. A control simulation with a spherical colloid—which the manuscript itself notes in §2.1.2 is achievable by switching off the filler conservative interactions—or an explicit calculation showing that the shape-induced free-energy landscape is identical for β and −β is required before the ordering difference can be attributed to the stresslet sign.","section":"§3.2, §3.3, SM Fig. 7"},{"comment":"The claim that pullers undergo the hexagonal-order transition at smaller Fg than pushers is based on visual comparison of P(ψ6) distributions and G6(r) curves from what appears to be a single independent run per state point. The text says averages are taken over 600 independent configurations, but no information is given about the time separation between them or the number of independent initial conditions; if the configurations come from a single trajectory, the differences at Fg/Fp = 1.5 and 2.25 could be within run-to-run fluctuations. With only 500 colloids (and fewer in the bottom layer at low Fg), finite-size fluctuations are non-negligible. I recommend reporting ⟨ψ6⟩ or the fraction of particles with ψ6 above a threshold as a function of Fg, with standard errors or confidence intervals from independent runs, and stating a quantitative criterion for 'hexagonal crystal' (e.g., a plateau in G6(r) or a decay exponent) to replace the qualitative 'virtually overlapping' comparison.","section":"§3.3, Figs. 5 and 6"}],"minor_comments":[{"comment":"In the sentence about colloidal beads, 'not only onserved the same non-linear behaviour' should read 'observed'.","section":"Introduction"},{"comment":"'zoom avaible in SM' contains a typo; it should be 'available'.","section":"§3.2"},{"comment":"The phrase 'previously reported bioconvenction' contains a typo; it should be 'bioconvection'.","section":"§4 Conclusions"},{"comment":"The caption defines the 'Static Structure Factor (SFF)', but the acronym used throughout the text is SSF; please make this consistent.","section":"Fig. 4 caption"},{"comment":"The table caption begins 'T able 1' and the footnote contains 'rencently'; both are typos.","section":"Table 1 and footnote in §2.1.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the model description is detailed. My main concern is the untested shape effect on the puller/pusher comparison; requesting a spherical-colloid control is reasonable given the authors' own statement that such a control is achievable. I see no other issues with novelty or attribution of prior work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a careful, honest simulation paper with a genuinely new observation: under strong gravity, puller squirmers form a nearly defect-free hexagonal monolayer at lower gravitational fields than pushers, and both active types anneal defects better than passive colloids. The DPD raspberry-squirmer model is validated against Stokes law for passive sedimentation and against earlier MPCD/LB results for active sedimentation lengths and wall accumulation, which gives me confidence in the model. The structural characterization is thorough—P(psi6), G6(r), g(r), SSF, and snapshots all point in the same direction—and the authors explicitly flag the missing lubrication forces and the kinetic trapping of the passive baseline. That is good practice.\n\nThe soft spot is the shape-bias assumption. In §3.2 they assert that the weak departure from sphericity affects pullers and pushers the same way, so the differences are hydrodynamic. But the propulsion reaction force is applied to the nearest filler particle, so the faceted 19-bead raspberry experiences different force distributions for the two swimmer types; the B2 term changes sign, and the polyhedral surface sets the local lever arms. The paper itself later attributes the high-gravity orientation peaks to the raspberry facets (SM fig. 7), which sits awkwardly with the \"same in both\" claim. Without a spherical-colloid control—which the authors describe as straightforward in §2.1.2—the puller/pusher ordering difference could be a model artifact rather than a property of stresslet sign. This is the main thing a referee should push on.\n\nTwo smaller issues: the structural claims lack error bars, and there is no quantitative order-parameter curve versus Fg or phase-boundary analysis, so the \"transition\" language is more qualitative than the text sometimes suggests. The data are not shipped, only available on request, which is common but not ideal.\n\nAll that said, the paper deserves peer review. The qualitative direction is plausible, the model validation is solid, and the central claim is worth testing. I would condition acceptance on addressing the shape-control point and adding error bars, not on a rewrite. This is a paper for people working on active sedimentation and DPD microswimmer models; a serious referee should take it.","headline":"A credible, honest DPD simulation study with a new qualitative claim about pullers annealing defects in sedimented monolayers, but the headline comparison rests on an untested shape-neutrality assumption.","tokens_in":689,"tokens_out":888,"would_cite":false,"duration_ms":41030,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In sedimenting suspensions of model microswimmers, increasing gravity drives the bottom monolayer into a hexagonal crystal, and puller-type swimmers reach a nearly defect-free crystal at lower gravitational forcing than pushers or passive…","keywords":["squirmer suspensions","sedimentation","dissipative particle dynamics","hexagonal order","active colloids","puller and pusher swimmers","wall orientation","monolayer crystallization"],"falsifier":"Run the same bottom-layer simulations with effectively spherical colloids, for example with conservative interactions only on the central thruster bead or with a raspberry made of many more fillers, and compare $\\psi_6$, the wall-angle peaks, and the $F_g/F_p$ at which crystallinity appears. If the puller advantage in preserving hexagonal order shrinks or the preferred wall angles move away from the polyhedron face angles reported in the supplementary material, the central claim is an artifact of raspberry shape rather than a property of the stresslet sign.","tokens_in":23686,"feed_emoji":"🦠","tokens_out":6895,"duration_ms":71104,"temperature":0.7,"pith_summary":"The paper uses coarse-grained dissipative-particle-dynamics simulations of 500 \"raspberry\" colloids to ask how gravity and self-propulsion jointly shape where microswimmers settle and what structure the sedimented layer takes. It shows that as the gravitational force increases, the bottom monolayer of passive colloids, pullers, and pushers all move toward a hexagonal crystal, but the path differs: passive colloids get kinetically trapped with defects, while active swimmers keep rearranging and anneal those defects. Among swimmers, pullers preserve hexagonal order better than pushers, forming an almost perfect hexagonal crystal at lower gravitational fields. The authors take this as evidence that the sign of the swimmer's stresslet, whether it pulls fluid in at the equator or pushes it out, controls both wall orientation and how easily a sedimented monolayer heals. These results matter because settling of biological microswimmers underlies biofilm formation and bioconvection, so swimmer-type-dependent order could affect how dense, structured microbial layers form.","feed_headline":"Puller microswimmers form near-perfect crystal beds under gravity","feed_subtitle":"Simulations show puller-type swimmers crystallize at lower gravity than pushers or passive colloids.","key_machinery":"The argument runs on the raspberry-DPD squirmer: each colloid is a rigid body of 19 DPD beads (18 fillers on a sphere plus a central thruster), and self-propulsion is imposed through a force field on solvent particles in a shell around the colloid, $\\mathbf{F}_H(r,\\theta) = (B_1\\sin\\theta + B_2\\sin\\theta\\cos\\theta)\\,\\hat{\\mathbf{e}}_\\theta$, with $\\beta = B_2/B_1$ positive for pullers and negative for pushers. The structural analysis is carried by the local hexagonal-order parameter $\\psi_6$, its correlation function $G_6(r)$, the static structure factor, and the distribution of wall-orientation angles $\\alpha$. The mechanism that distinguishes swimmer types is the wall-induced hydrodynamic torque: pullers reorient toward the wall and stay persistently tilted, while pushers prefer orientations parallel to the wall, making them more motile and more disruptive to the hexagonal layer.","core_discovery":"On the paper's own terms, the central discovery is that a sedimented bed of active colloids is not passively frozen: activity lets the bottom layer escape kinetic trapping and reorganize into an ordered state, and the type of activity selects how perfect that order is. At the two strongest gravitational fields, puller suspensions show a pronounced peak in the hexagonal-order parameter at $\\psi_6 \\approx 0.95$ already at $F_g/F_p = 1.50$, and their hexagonal correlation function $G_6(r)$ does not decay, as expected for a crystal; pushers need the highest field, $F_g/F_p = 2.25$, to reach the same crystalline signature, and at $F_g/F_p = 1.50$ they show a power-law decay consistent with a hexatic phase. Passive colloids form a hexagonal structure full of defects because thermal motion cannot overcome the trapped states. The explanation offered is hydrodynamic: pullers tend to point toward the wall with stable orientations, which hinders their lateral mobility and protects the lattice, whereas pushers align more parallel to the wall, stay more motile, and are more prone to disrupt the layer.","pith_inferences":["If stresslet sign alone drives the defect-annealing gap, then a purely steric active colloid without hydrodynamics should lose the puller/pusher difference; this is a direct numerical lever to test the proposed mechanism.","By the same logic, pusher-like bacteria such as E. coli should produce more defective sedimented monolayers than puller-like algae, which could give a physical explanation for heterogeneous biofilm structure.","The near-perfect hexagonal order achieved by pullers suggests a possible route to gravity-assisted colloidal crystal templating using active colloids, where defects self-heal before the bed freezes.","A natural next step the paper leaves implicit is to measure how long a single defect takes to anneal as a function of $\\beta$ and $F_g$; if the annealing rate is maximal at intermediate puller strength, the effect would be tunable."],"forward_implications":["If puller suspensions crystallize at lower gravity, gravitational settling of puller-like microalgae will produce sedimented monolayers with fewer defects at the same forcing, which changes how one predicts biofilm and bio-sediment microstructure from swimmer type.","Since activity removes kinetic trapping, active sedimented layers can anneal defects on simulation timescales where passive layers cannot; the quality of the crystal is set by the competition between gravitational compression and swimmer mobility.","The exponential sedimentation regime at moderate gravity reproduces the known trend that sedimentation lengths are larger for pushers than pullers and shrink with stronger gravity, so the model is consistent with earlier LB and MPCD results.","At high gravity, wall orientations converge to a common roughly 120 degrees toward-wall and 60 degrees away-from-wall population, meaning hydrodynamic swimmer-type differences are strongest in the transition regime rather than in the deep-crystal regime."],"supporting_citations":[{"why":"Supplies the raspberry-DPD squirmer model and the tangential force field ($B_1$, $B_2$) on which all simulations are built.","marker":"[52]"},{"why":"The MPCD study of collective squirmer sedimentation whose sedimentation profiles and lengths the paper validates and compares against.","marker":"[36]"},{"why":"The LB study of squirmers under gravity used for comparison of sedimentation and of wall-orientation distributions.","marker":"[35]"},{"why":"Showed hexagonal order in a gravity-confined squirmer monolayer, the phenomenon this paper extends to a defect-free crystal.","marker":"[40]"},{"why":"Provides the single-squirmer sedimentation regimes and the wall-distance-dependent friction estimate used to interpret the passive sedimentation velocity.","marker":"[25]"},{"why":"Describes pairwise far-field hydrodynamic interactions near a wall that the paper uses to explain why pushers access higher regions and leave the first layer.","marker":"[48]"},{"why":"Reports wall-orientation angles near 60 degrees and 120 degrees in confined ellipsoidal squirmers, a reference for the tilt angles observed here.","marker":"[51]"},{"why":"Gives the far-field theory that pullers align perpendicular to walls while pushers align parallel, the central mechanism invoked for puller/pusher differences.","marker":"[89]"}],"fun_headline_variants":["Puller swimmers crystallize best under gravity","Gravity crystals: pullers outperform pushers","Active colloids anneal defects in sediment beds","Pullers form near-perfect crystal beds under gravity","Squirmer suspensions: pullers win the crystal race"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion that pullers preserve hexagonal order better than pushers rests on the assumption that the small departures from perfect sphericity in the 19-bead colloid affect pullers and pushers identically, so every difference between them is caused by their swimming mechanism rather than by shape.","fun_headline_variants_meta":{"raw":{"variants":["Puller swimmers crystallize best under gravity","Gravity crystals: pullers outperform pushers","Active colloids anneal defects in sediment beds","Pullers form near-perfect crystal beds under gravity","Squirmer suspensions: pullers win the crystal race"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000159,"raw_usage":{"total_tokens":1225,"prompt_tokens":940,"completion_tokens":285,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":212}},"tokens_in":556,"tokens_out":285,"duration_ms":3454,"temperature":1.0,"reasoning_tokens":212,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:31:06.006071+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same bottom-layer simulations with effectively spherical colloids, for example with conservative interactions only on the central thruster bead or with a raspberry made of many more fillers, and compare $\\psi_6$, the wall-angle peaks, and the $F_g/F_p$ at which crystallinity appears. If the puller advantage in preserving hexagonal order shrinks or the preferred wall angles move away from the polyhedron face angles reported in the supplementary material, the central claim is an artifact of raspberry shape rather than a property of the stresslet sign.","supporting_citations":[{"cited_title":"Schaar, A","cited_arxiv_id":null,"evidence_quote":"Supplies the raspberry-DPD squirmer model and the tangential force field ($B_1$, $B_2$) on which all simulations are built."},{"cited_title":"Chat \\'e , F","cited_arxiv_id":null,"evidence_quote":"The MPCD study of collective squirmer sedimentation whose sedimentation profiles and lengths the paper validates and compares against."},{"cited_title":"Vicsek, A","cited_arxiv_id":null,"evidence_quote":"The LB study of squirmers under gravity used for comparison of sedimentation and of wall-orientation distributions."},{"cited_title":"Scagliarini and I","cited_arxiv_id":null,"evidence_quote":"Showed hexagonal order in a gravity-confined squirmer monolayer, the phenomenon this paper extends to a defect-free crystal."},{"cited_title":"Delfau, J","cited_arxiv_id":null,"evidence_quote":"Describes pairwise far-field hydrodynamic interactions near a wall that the paper uses to explain why pushers access higher regions and leave the first layer."},{"cited_title":"Ishimoto and E","cited_arxiv_id":null,"evidence_quote":"Reports wall-orientation angles near 60 degrees and 120 degrees in confined ellipsoidal squirmers, a reference for the tilt angles observed here."},{"cited_title":"Happel and H","cited_arxiv_id":null,"evidence_quote":"Gives the far-field theory that pullers align perpendicular to walls while pushers align parallel, the central mechanism invoked for puller/pusher differences."}],"review_version":1}