{"id":"00b71211-8e8d-42b4-a498-c76dfd90b241","arxiv_id":"2411.13893","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"SPD with a novel slip boundary treatment reproduces squirmer motion, flow fields, hydrodynamic interactions, thermal fluctuations, and multiphase droplet behavior.","lead":"This paper shows that a computational fluid dynamics method called smoothed particle dynamics (SPD) can simulate 'squirmer' microswimmers by giving the solid boundary particles an artificial slip velocity. This could be useful for studying how swimming microorganisms and synthetic swimmers behave in complex fluids and droplets.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (36) enforces slip through a regularized extrapolation (denominator dA+αh, not dA+dB) and is never tested for α- or curvature-sensitivity, leaving the central accuracy claim conditional.","rationale":"I focused on Eq. (36) because it is the mechanism by which the squirmer's surface slip is converted into fluid forces; every subsequent validation inherits it. The reader's weakest assumption points to the same region, but I sharpen it: the formula is not even a true two-point interpolation because of the αh regularization, so the linear-profile assumption is not exactly realized in the implementation. The paper's good agreement for the tested parameters and the available resolution study support the central claim, but the absence of an α-sweep and of a high-curvature test leaves a genuine soft spot. I also noticed two reporting inconsistencies that should be corrected: in Sec. III.B the text says 'pushers with β=5' and 'pullers with β=-5', which reverses the convention fixed by Eq. (3); and in Sec. III.C the stated parameters η=15, kBT=1, R=1, U0=2 give Per≈377, not the reported 25.1. These appear to be typographical/editorial issues rather than deep flaws, so they do not overturn the central validation; they are additional reasons to request minor revisions while keeping the overall verdict unchanged.","tokens_in":15324,"tokens_out":16476,"duration_ms":179897,"concrete_test":"Re-run the single-squirmer steady-state benchmark (Re=0.01, B1=0.015) at fixed R/Δx=10 with α = 0.001, 0.01, 0.05, and 0.1, and at α=0.05 with R/Δx=5, 10, 16, and 20. Extract the swimming speed U0 and the surface slip profile using the same procedure as Fig. 2. If U0 and the slip profile remain constant to within about 2% across these runs, the αh term is a benign regularization and the concern is settled. If they shift systematically with α or resolution, then the effective boundary condition is not exactly the prescribed slip velocity in Eq. (6), and the accuracy claims require qualification.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (36) is the only place the squirmer's slip velocity enters the dynamics: v_B = vs + v_C - dB/(dA+αh)(v_A - vs - v_C). This is not the two-point interpolation between fluid particle A and boundary particle B described in the text; a genuine interpolation would use dB/(dA+dB). The αh replacement, with α=0.05 and h=1.2Δx, means the surface condition at point C is enforced only up to a residual that depends on dB, dA, and αh. For fluid particles near the surface, with dA around 0.1–0.5Δx, the correction for individual pairs can be a substantial fraction of (v_A - vs - v_C). The simulations fix R/Δx=10 and α=0.05, and the paper never varies α or reports a curved-surface convergence test. The resolution study in Fig. 5 changes Δx but not α, so it does not separate the intended boundary condition from the regularization error. Because the central claim is that the model accurately represents pushers, neutral swimmers, and pullers across geometries and scales, this unquantified regularization is the weakest load-bearing point in the argument.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a smoothed particle dynamics (SPD) framework for simulating squirmers, in which the tangential slip velocity at the squirmer surface is imposed through artificial velocities assigned to boundary particles during the pairwise dissipative force calculation. The method is validated against analytical Stokes-flow solutions for a single squirmer, perturbation-theory results at small Reynolds number, prior simulations of two-squirmer collisions and near-wall dynamics, and equilibrium statistical-mechanical predictions for velocity and angular-velocity autocorrelation functions. The paper closes with a qualitative demonstration of a squirmer encapsulated in a droplet in a multiphase flow, reporting distinct co-swimming behaviors for pushers, neutral swimmers, and pullers. The central claim is that the SPD-squirmer model accurately represents a range of microswimmer types across scales and flow regimes.","tokens_in":15597,"tokens_out":3306,"duration_ms":35601,"significance":"If the central claim holds, the paper offers a unified Lagrangian solver for squirmer hydrodynamics that handles moving fluid-solid boundaries, multiphase interfaces, and thermal fluctuations within a single framework. The validation is unusually broad: single-squirmer slip-velocity distributions, flow-field decay, swimming speed versus β, pairwise collisions, wall-bounded motion, and Brownian dynamics are all compared with external analytical or numerical benchmarks. A notable strength is that no parameter is fitted to the target swimming quantities; the slip velocity is prescribed, and the swimming speed emerges dynamically. The multiphase section, however, is only qualitative, and the key boundary-condition formulation contains a regularization term whose influence is not quantified. The paper is therefore a potentially valuable methods contribution, but two load-bearing points need additional support.","major_comments":[{"comment":"The text states that the artificial velocity v_B is assigned so that linear interpolation between v_A and v_B satisfies the slip condition at the surface point C, but Eq. (36) uses a denominator d_A + αh, not d_A + d_B. A genuine interpolation would give v_B = v_C - (d_B/d_A)(v_A - v_C - v_s) - v_s, so the regularization with αh means the desired surface velocity is enforced only approximately, with an error that depends on d_A, d_B, and αh. Since the paper's abstract claims that the boundary-condition treatment accurately represents pushers, neutral swimmers, and pullers, the sensitivity of all reported results to α should be established. The resolution study in §III.A changes Δx but keeps α fixed at 0.05, so it does not separate discretization error from this regularization error. Please add an α-sensitivity study (e.g., α from 0.01 to 0.2 at fixed resolution) for the swimming speed and flow field, and report the residual slip-velocity error at the surface as a function of d_A and αh.","section":"§II.D, Eq. (36)"},{"comment":"The multiphase results are presented without any quantitative validation. The paper reports different co-swimming behaviors for pushers, neutral swimmers, and pullers in a droplet, but there is no comparison with an analytical solution, a resolved computational benchmark, or an experimental observation; the only quantitative statement is that 'the results converge as the resolution increases' (just before Fig. 12), which is not shown with error metrics. Given that the abstract and introduction emphasize multiphase interface tracking as a motivation for the SPD approach, the accuracy claim for multiphase squirmer dynamics is not yet supported. Please provide a convergence study with quantified errors for the droplet-squirmer system, or at least a grid-convergence table for the steady-state squirmer velocity and droplet deformation.","section":"§III.D"}],"minor_comments":[{"comment":"The definition 'v_B = (v^τ_A, v^n_A, v^k_A)' appears to be a typo: the right-hand side should be the fluid particle velocity v_A, not v_B.","section":"§II.D, Eq. (34)"},{"comment":"'parawise dissipative force' should be 'pairwise dissipative force'.","section":"§II.D, before Eq. (37)"},{"comment":"The caption says 'For pushers with β = 5', which conflicts with the convention β < 0 for pushers stated in Eq. (3) and used in the text; this is likely a sign error and should be corrected.","section":"Fig. 7 caption"},{"comment":"The VACF and AVACF plots would benefit from showing standard deviations or confidence bands across the 20 independent runs, since the ensemble size is modest and the decay at long times is where statistical noise is most visible.","section":"§III.C, Fig. 9"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid methods contribution with extensive validation, but the boundary-condition regularization in Eq. (36) is a genuine load-bearing concern: the method is claimed to impose the slip velocity exactly, yet the αh term makes the enforcement inexact, and no sensitivity study is provided. The multiphase section is also exclusively qualitative. Both issues are addressable in revision, so I recommend major revision rather than rejection. The authors should also clarify the relationship with their earlier work on slip boundary conditions (ref. 56), since the novelty relative to that paper is not sharply delineated in the present text."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a solid, genuinely new methods paper. It is the first squirmer model in smoothed particle dynamics, and the authors validate it against a lot of external results. The central claim holds up.\n\nWhat is actually new: they extend their earlier slip boundary condition (Cai, Li, Bian, JCP 2023) to the squirmer's prescribed surface slip. The slip velocity is an input, not a fit; the swimming speed, flow field, collision trajectories, and near-wall reorientation all emerge and match analytic or published results. The VACF/AVACF checks at the mesoscale are a nice addition. The multiphase demonstration is qualitative but shows the method's reach.\n\nThe weak spot, and it is a modest one: Eq. (36) is not the linear interpolation the text describes. The denominator is dA + alpha*h, not dA + dB, with alpha = 0.05 and h = 1.2*dx. For a fluid particle sitting 0.1-0.5*dx from the surface, the regularization term is comparable to dA, so the slip condition is enforced only approximately. The paper never varies alpha or checks curvature sensitivity. The resolution study changes dx but keeps alpha fixed, so it does not isolate this error. That said, the validations are accurate enough that the method clearly works. The gap is a missing robustness study, not a demonstrated flaw. A referee should ask for one alpha-variation test and a comment on the interpolation error.\n\nTwo smaller things: the multiphase section is purely qualitative, with no error measures or comparison to other methods. And there is no code or data deposit, which would help a lot for a methods paper like this.\n\nFor a reader: this is for people working on microswimmer simulations, especially SPH/SDPD users, and for the squirmer community looking for a Lagrangian alternative to lattice Boltzmann or MPCD. It is a serious piece of work and should be sent to peer review. I would accept it with minor revisions, conditional on an alpha-sensitivity test and a clearer statement about the regularization.","headline":"First SPD squirmer model with thorough validation; the only real gap is an untested regularization parameter in the slip boundary condition.","tokens_in":16088,"tokens_out":3356,"would_cite":true,"duration_ms":31971,"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":"This paper introduces a boundary treatment that assigns artificial slip velocities to boundary particles, enabling one Lagrangian particle solver to reproduce pusher, neutral, and puller squirmers from clean fluids to droplets with…","keywords":["squirmer","microswimmer","smoothed particle dynamics","slip boundary condition","hydrodynamic interaction","multiphase flow","thermal fluctuations","pusher and puller swimmers"],"falsifier":"Use a squirmer with radius comparable to the particle smoothing length so that the gap from a fluid particle to the tangent plane is only one or two particle spacings, and compare the measured swimming speed and near-field flow to the analytic Stokes solution; a visible departure from the zero-Reynolds speed $\\frac{2}{3}B_1$ would show that the linear interpolation in Eq. (34) no longer enforces the prescribed slip velocity.","tokens_in":15133,"feed_emoji":"🦠","tokens_out":12462,"duration_ms":101325,"temperature":0.7,"pith_summary":"This paper claims that smoothed particle dynamics (SPD), a Lagrangian particle method, can serve as a unified solver for the squirmer model of microswimmers. The essential step is a boundary treatment that assigns an artificial velocity to each solid boundary particle so that the fluid feels the prescribed tangential slip velocity on the sphere's surface. If this works, a single numerical framework can cover a single swimmer's steady velocity and flow field, hydrodynamic interactions with walls and other squirmers, mesoscale Brownian motion, and multiphase droplets. The paper validates each setting against analytic Stokes-flow solutions and published reference trajectories, reporting good agreement at the Reynolds numbers tested.","feed_headline":"One slip-velocity rule lets smoothed particle dynamics model squirmers","feed_subtitle":"The same boundary treatment covers pushers, pullers, walls, collisions, thermal noise, and droplets.","key_machinery":"The central object is the artificial-velocity interpolation expressed in Eq. (34). For every fluid particle $A$ near the squirmer, the paper draws a normal to the tangent plane at the nearest surface point, then assigns each neighboring boundary particle $B$ an artificial velocity $v_B$ by linearly extrapolating the fluid velocity $v_A$ across the gap so that the interpolated value at the surface point equals the prescribed surface velocity $v_C + v_s$. This artificial velocity enters only the dissipative force, leaving the squirmer's rigid-body kinematics untouched; it turns the geometric slip condition into a local pairwise force. Because the tangent plane is defined locally for each fluid particle, the same rule can in principle handle arbitrary boundary geometry, and the paper combines it with pressure interpolation from fluid to boundary particles to maintain a correct pressure gradient at the surface.","core_discovery":"On its own terms, the discovery is that a prescribed surface slip velocity can be injected into an SPD simulation by changing only the pairwise dissipative force between fluid and boundary particles: each boundary particle receives a different artificial velocity for each fluid neighbor, chosen so that the velocity linearly interpolated across the gap matches the rigid-body plus slip velocity at the tangent point. This yields the correct steady swimming speed $U_0 = \\frac{2}{3}|B_1|$, the expected pusher/neutral/puller flow-field structures, and the correct far-field decay. The same construction carries through when random thermal forces are added, and it coexists with the surface-tension force in a squirmer-in-droplet setup, where the authors observe qualitatively different co-swimming behavior depending on squirmer type and droplet size.","pith_inferences":["If the linear-interpolation assumption survives thinner boundary layers, the same tangent-plane construction could implement arbitrary prescribed surface slip for any immersed body in SPH, not only spherical squirmers.","A natural extension is non-spherical active particles such as ellipsoidal squirmers; the local tangent-plane machinery is ready for them, but the interpolation error would need to be checked where curvature is high.","The squirmer-inside-droplet results suggest that the composite object behaves as an effective squirmer whose polarity reverses for pusher and neutral swimmers; measuring the dipole strength of the droplet-swimmer pair in the far field would provide a clean test.","Because the artificial velocity is recomputed per fluid neighbor, the scheme retains the local structure of SPD, so it should scale to dense suspensions of many swimmers in the same framework."],"forward_implications":["A single SPD implementation can now address single-swimmer flow fields, squirmer-wall and squirmer-squirmer interactions, thermally fluctuating squirmers, and a squirmer inside a droplet without changing the boundary method.","The model reproduces the analytic result that the swimming speed is set by the first mode $B_1$ alone, while the squirmer parameter $\\beta = B_2/|B_1|$ selects pusher, neutral, or puller character.","Near a wall at Reynolds number 1, the model reproduces the three published trajectory regimes: escape, oscillate-and-slide along the wall, and forward bouncing.","With thermal fluctuations, the model recovers equipartition values and the long-time decay of the velocity and angular velocity autocorrelation functions at the mesoscale.","In a multiphase setting, the model predicts that pusher and neutral squirmers co-swim with their droplet in a direction opposite to their head orientation, whereas a puller drags the droplet along its heading after puncturing the interface."],"supporting_citations":[{"why":"Introduces the spherical squirmer as a swimming-body model whose surface-velocity picture the paper uses.","marker":"[16]"},{"why":"Refines the squirmer to an effective slip velocity on a sphere and supplies the two-mode form with $B_1$ and $B_2$ used here.","marker":"[17]"},{"why":"Supplies the smoothed dissipative particle dynamics discretization and fluctuation-dissipation structure that the method builds on.","marker":"[47]"},{"why":"Provides the earlier arbitrary-slip boundary treatment from which the present artificial-velocity construction is adapted.","marker":"[56]"},{"why":"Gives the generalized wall boundary condition used to interpolate boundary-particle pressure from the surrounding fluid.","marker":"[61]"},{"why":"Supplies the reference near-wall trajectory regimes at Reynolds number 1 that the model reproduces.","marker":"[20]"},{"why":"Supplies the Stokes-flow two-squirmer collision trajectory used as the pusher reference.","marker":"[22]"},{"why":"Supplies the inertial-regime two-squirmer collision trajectory used as the puller reference.","marker":"[23]"},{"why":"Provides the small-Reynolds perturbation expansion used to compare swimming speeds at nonzero inertia.","marker":"[62]"}],"fun_headline_variants":["SPD with slip-velocity boundary simulates all squirmer types","Slip-velocity rule lets SPD model pushers and pullers","Smoothed particle dynamics captures squirmer hydrodynamics","SPD squirmer model works for collisions and droplets","Slip-velocity boundary trick powers SPD squirmer sims"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the velocity profile between a fluid particle and a nearby boundary particle is linear, so that an artificial boundary-particle velocity can be chosen to enforce the prescribed slip at the surface; a violation of that linearity in thin gaps or strongly curved regions would break the swimming speed and flow field.","fun_headline_variants_meta":{"raw":{"variants":["SPD with slip-velocity boundary simulates all squirmer types","Slip-velocity rule lets SPD model pushers and pullers","Smoothed particle dynamics captures squirmer hydrodynamics","SPD squirmer model works for collisions and droplets","Slip-velocity boundary trick powers SPD squirmer sims"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000243,"raw_usage":{"total_tokens":1547,"prompt_tokens":979,"completion_tokens":568,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":484}},"tokens_in":595,"tokens_out":568,"duration_ms":5832,"temperature":1.0,"reasoning_tokens":484,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:45:26.900787+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Use a squirmer with radius comparable to the particle smoothing length so that the gap from a fluid particle to the tangent plane is only one or two particle spacings, and compare the measured swimming speed and near-field flow to the analytic Stokes solution; a visible departure from the zero-Reynolds speed $\\frac{2}{3}B_1$ would show that the linear interpolation in Eq. (34) no longer enforces the prescribed slip velocity.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the spherical squirmer as a swimming-body model whose surface-velocity picture the paper uses."},{"cited_title":"Espa˜ nol and M","cited_arxiv_id":null,"evidence_quote":"Supplies the smoothed dissipative particle dynamics discretization and fluctuation-dissipation structure that the method builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the earlier arbitrary-slip boundary treatment from which the present artificial-velocity construction is adapted."},{"cited_title":"Adami, X","cited_arxiv_id":null,"evidence_quote":"Gives the generalized wall boundary condition used to interpolate boundary-particle pressure from the surrounding fluid."},{"cited_title":"Li and A","cited_arxiv_id":null,"evidence_quote":"Supplies the reference near-wall trajectory regimes at Reynolds number 1 that the model reproduces."},{"cited_title":"Ishikawa, M","cited_arxiv_id":null,"evidence_quote":"Supplies the Stokes-flow two-squirmer collision trajectory used as the pusher reference."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the inertial-regime two-squirmer collision trajectory used as the puller reference."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the small-Reynolds perturbation expansion used to compare swimming speeds at nonzero inertia."}],"review_version":1}