{"id":"35cba4c9-55cc-4cf7-b6a5-0738631647bd","arxiv_id":"2509.08295","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Differences in membrane surface tension, not just viscosity or bending stiffness, can drive vesicle swinging and tumbling in shear flow, and simulations reveal three previously unobserved motion regimes.","lead":"Using 3D computer simulations of a phase-field membrane model, this paper shows that differences in surface tension across two lipid phases can make a vesicle swing or tumble in shear flow. It also maps a phase diagram with six bending-driven motion regimes, three of them new, which could help design vesicle-based drug carriers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Inextensibility relaxation parameter ξ is never specified and no ξ→0 or grid-convergence check is shown; the central surface-tension-driven swinging/tumbling may be an artifact of membrane extensibility.","rationale":"The reader's CONDITIONAL verdict is appropriate. My stress-test digs deeper: the unstated ξ is more load-bearing than the initial-condition issue because it attacks the physicality of the mechanism itself. Equations (5)/(20) show that local inextensibility is only enforced as ξ→0; a finite ξ permits local area dilatation. The new 'surface tension heterogeneity' energy term is then able to change local area, giving it a driving force absent in an inextensible vesicle. Since ξ never appears with a value in the parameter lists (Sections 3.2, 4.1, 4.2), the central simulations are not reproducible and the new regimes (VII,VIII) cannot be distinguished from regularization artifacts. The reader's initial-condition concern is valid but secondary: even if symmetric tip domains are special, the mechanism could still be real; however, if the effect disappears at ξ→0 or at finer resolution, the mechanism is not physical. I therefore recommend keeping the verdict CONDITIONAL, with the added condition that the authors report ξ and demonstrate convergence to the inextensible limit and to grid independence. The paper deserves credit for a thermodynamically consistent derivation and successful reproduction of experimental phase separation and tank-treading/tumbling at σ_S=1, but these validations do not exercise the new surface-tension-heterogeneity term.","tokens_in":20427,"tokens_out":13480,"duration_ms":159618,"concrete_test":"Rerun the representative cases of §4.2 that define regimes VII (σS=0.85, Ca≈0.02, κB=0.65) and VIII (σS=0.85, Ca=1), first with the currently unstated ξ, then with ξ reduced by factors of 10 down to 10^-5 and with grid 256^3. Monitor max|P:∇u| on the membrane as a direct measure of inextensibility, and record the inclination-angle trajectory and regime. If tumbling/swinging persists as max|P:∇u|→0, the mechanism survives the check; if the regime shifts or disappears as ξ is reduced, or if max|P:∇u| remains O(1), the central claim is a numerical artifact of finite-ξ membrane extensibility and should be rejected or withdrawn until recomputed at strict inextensibility.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—surface-tension heterogeneity alone triggers swinging/tumbling—is supported only by §4.2 simulations in which the local-inextensibility relaxation parameter ξ (Eqs. (5),(20); Table 1; Appendix C Step 4) is never assigned a value. Eq. (5) enforces P:∇u=0 only in the limit ξ→0; for finite ξ the membrane is extensible. The surface-energy term ∫σ_S(c)ζ(φ)dV then couples to area changes, so phase-dependent 'surface tension' can do mechanical work that would be forbidden in a strictly inextensible vesicle. In that regime, the observed tumbling/swinging could be a regularization artifact rather than a physical mechanism. No ξ-sensitivity, ξ→0, or mesh-refinement study is reported, and the phase diagrams (Figs. 13,16) are point-sampled once at 128^3 with a single initial condition (v=0.91, two symmetric tip domains, a_ld=0.4). The experimental validation (§3.2) also sets σ_S=1, so it does not test the heterogeneous-surface-tension term. Thus the headline mechanism is not yet established to be independent of numerical regularization.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a thermodynamically consistent phase-field model for multicomponent vesicles in shear flow, coupling a bulk Navier–Stokes/Cahn–Hilliard system with a diffuse-domain surface Cahn–Hilliard equation for lipid phases. The membrane energy includes bending, surface tension, line tension, and surface-area penalty terms, and local inextensibility is imposed through a harmonic-relaxation Lagrange multiplier. The model is validated against GUV phase-separation experiments [5] and shear-flow experiments [33], including a quantitative comparison of inclination-angle oscillations. The authors then present phase diagrams in the (bending-rigidity contrast, capillary number) and (surface-tension contrast, capillary number) planes, reporting three new bending-driven regimes (lateral ring banding, lateral broken-ring banding, vertical broken-ring banding) and a new surface-tension-driven swinging/tumbling mechanism that operates without viscosity contrast or asymmetric phase distributions.","tokens_in":20787,"tokens_out":5040,"duration_ms":56474,"significance":"If the reported results are robust, the paper offers a useful 3D computational framework for multicomponent vesicles and identifies a previously unrecognized control knob—membrane surface-tension heterogeneity—for vesicle dynamics in shear. The model derivation is a genuine strength: Appendix B supplies a free-energy dissipation identity (B.25) with explicit nonnegative dissipation terms, and the validation against experiments gives nontrivial support for the model's physical content. The inclination-angle comparison in Fig. 12 is particularly valuable. However, the central new claims—the three new bending regimes and the surface-tension-driven swinging/tumbling mechanism—are supported only by simulations whose inextensibility regularization parameter is never reported and for which no mesh or ξ-convergence study is shown, and the phase diagrams are built on a single initial geometry/phase configuration. These gaps must be addressed before the paper's headline conclusions can be considered established.","major_comments":[{"comment":"The inextensibility relaxation parameter ξ is listed in Table 1 but never assigned a value in any numerical setup, including the phase-diagram runs of §4.1–4.2. Eq. (20) enforces P:∇u=0 only in the limit ξ→0; for finite ξ the membrane is extensible, and the surface-energy term ∫σ_S(c_Γ)ζ(ϕ)dV can then do mechanical work that would be forbidden for a strictly inextensible vesicle. The §3.2 validation sets σ_S=1 and therefore does not test the heterogeneous-surface-tension term. Without a reported ξ, a ξ→0 (or sufficiently small) convergence study, and a spatial-resolution check for the swinging/tumbling regimes (VII) and (VIII), the headline surface-tension mechanism cannot be distinguished from a regularization artifact. This is load-bearing for the paper's central claim.","section":"§4.2, Eq. (20), Table 1, Appendix C Step 4"},{"comment":"The phase diagrams and the regime taxonomy are generated from a single initial condition: a prolate vesicle with reduced volume v≈0.91, two symmetric ld domains pre-positioned at the tips, and area fraction a_ld=0.4 (§4.1; §4.2 states 'The numerical setup follows Section 4.1'). No sensitivity study is reported with respect to initial vesicle aspect ratio, initial domain placement or area fraction, or stochastic phase-field noise. Since the claim is that surface-tension heterogeneity (and bending heterogeneity) generically drives these regimes, the authors should show representative checks that the regimes and their boundaries are not artifacts of this one prepared initial state. A single 128^3 point-sampled run per parameter pair is also insufficient to establish convergence of the regime boundaries.","section":"§4.1, §4.2, Figs. 13 and 16"},{"comment":"The abstract and conclusion describe the model as 'quantitatively validated,' but the validation is partly calibrated: in §3.1 the lipid mobility ratio Cn_Γ/Pe_Γ=0.12 and the interface thickness ε_c=0.02 are fitted/adjusted to the experimental data, and in §3.2 ε_c and ε_ϕ are 'fitted to experimental GUVs.' This is not a fatal flaw, but the claim of parameter-free quantitative prediction should be softened, and the fitting ranges and the sensitivity of the validation to those fitted values should be stated explicitly.","section":"§3.1, §3.2, and Abstract"}],"minor_comments":[{"comment":"Notation is inconsistent: the dimensionless surface-tension strength appears as Cs in Eq. (26) and as Cs_Γ in Table 1 and Eq. (C.29). Section 3.2 also uses 'Cs=27' without the subscript. Please unify.","section":"Eqs. (25)–(26), Table 1"},{"comment":"The text says that as σ_S decreases to 0.85, ld domains form a complete ring and its periodic rupture/reconnection 'result[s] in regime (VIII).' But regime (VIII) is later defined as tumbling with ld domains 'stably anchored at high-curvature tips.' This apparent contradiction should be clarified, since the description sounds like regime (VII) swinging.","section":"§4.2, paragraph after Fig. 16"},{"comment":"The line-tension delta function is written as ζ(c_Γ) in Eq. (16), but the variational derivatives in Appendix B and the chemical potential (26) use expressions involving η_c f'(c_Γ)/ε_c − η_c ε_c Δc_Γ. The relationship between δ(c_Γ) and ζ(c_Γ) should be stated precisely, and the signs in the line-tension contributions to ω_c should be checked against the energy (13).","section":"§2.3 and Appendix B"},{"comment":"The local inextensibility correction in Eq. (C.41) uses ∇·u^n in the constraint, while the velocity field has already been updated to u^{n+1} in Step 3. The temporal staggering is likely intentional, but it should be stated explicitly, along with the order of accuracy implied for the constraint.","section":"Appendix C, Step 4"}],"recommendation":"major_revision","confidential_remarks":"The paper has a solid thermodynamic foundation and the validations are meaningful, but the central new claims rest on numerical explorations whose regularization parameter is unspecified and whose initial-condition and grid sensitivity are not demonstrated. The missing ξ values and convergence checks are straightforward to supply and would substantially strengthen the paper. I would also encourage the authors to provide at least one or two variations of the initial phase distribution/aspect ratio to show that the surface-tension-driven swinging/tumbling is not tied to the particular symmetric tip-domain setup. If those robustness checks confirm the reported regimes, the paper would be a strong candidate for acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has real substance. The authors build a thermodynamically consistent phase-field model for multicomponent vesicles in shear flow, derive a clean free-energy dissipation law in Appendix B, and validate against two experiments: phase-separation coarsening on GUVs and shear-driven tank-treading/swinging/tumbling. The inclination-angle comparison in Fig. 12 is genuinely quantitative, not just qualitative. On top of that, they identify three new bending-driven regimes (lateral ring banding, broken-ring variants) and, more interestingly, show that surface-tension heterogeneity alone—without viscosity contrast or asymmetric phase placement—can produce swinging and tumbling. That is a plausible new physical control knob.\n\nThe main soft spot is the one the stress-test flags: the inextensibility relaxation parameter ξ appears in Eqs. (5), (20), and the Step 4 update, but no simulation in the paper assigns it a value, and there is no ξ-sensitivity or ξ→0 convergence check. Since the surface-tension term ∫σ_S(c)ζ(ϕ) does work when the membrane can stretch, the central swinging/tumbling claims could be an artifact of a too-loose inextensibility constraint. The area-penalty term (M_S=2000) presumably keeps area nearly fixed, but that is not a substitute for reporting ξ. This is a concrete omission, not a manufactured concern.\n\nOther soft spots are proportional to the claims. The phase diagrams in Figs. 13 and 16 are computed from one initial condition (v=0.91 prolate, two symmetric tip domains, a_ld=0.4) and sampled once at 128^3. That is thin support for a regime taxonomy; robustness to initial geometry and mesh refinement is not shown. The shear-flow validation in §3.2 sets σ_S=1, so it never tests the heterogeneous surface-tension term being used for the new regimes. And several parameters (ε_c, ε_ϕ, lipid mobility) are fitted, and line tension is varied to reproduce the experimental regimes—so that validation is partly calibration, though the trends are believable.\n\nNone of this is load-bearing enough to warrant rejection. The math in Appendix B looks self-consistent, the experimental comparisons give confidence in the baseline model, and the new mechanism is worth taking seriously. What the paper needs before its claims should be accepted is a clearly stated ξ (with sensitivity runs), a mesh-convergence check on at least one new regime, and ideally a second initial condition for the phase diagram.\n\nI would send this to peer review. It is a substantial modeling contribution with a novel physical hypothesis, and a good referee can demand the missing sensitivity analysis. For my own work, I would cite the model and the phase diagram with caution, pending that verification.\n\nI'd bring it to reading group as a useful example of how to build and benchmark a vesicle simulation, even while debating whether the new regimes survive a stricter inextensibility limit.","headline":"Solid 3D phase-field study with a genuinely new surface-tension contrast mechanism, but the unstated inextensibility relaxation parameter and single-shot phase diagrams mean the headline claims need a sensitivity pass before I'd trust them.","tokens_in":21241,"tokens_out":1910,"would_cite":true,"duration_ms":24028,"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":"The paper claims that membrane surface-tension contrast alone—without viscosity mismatch or asymmetric phase placement—can drive vesicles to swing and tumble under shear, and that bending-rigidity contrast organizes a phase diagram with thr","keywords":["multicomponent vesicles","shear flow","phase-field model","surface tension contrast","bending rigidity contrast","vesicle dynamics","lipid domains","phase diagram"],"falsifier":"Run the same parameters but replace the symmetric tip domains with a random or single off-center ld domain, or change the reduced volume to v=0.8; if surface-tension-driven swinging/tumbling disappears or the six-regime taxonomy collapses into different states, the mechanism is not generic. Alternatively, a microfluidic experiment on a giant unilamellar vesicle with matched bending rigidities but σ_S<1 should show the predicted tumbling at Ca≈1; failure to observe it would contradict the claim.","tokens_in":20319,"feed_emoji":"🫧","tokens_out":4914,"duration_ms":47924,"temperature":0.7,"pith_summary":"The authors build a thermodynamically consistent phase-field model coupling vesicle-fluid flow, membrane shape, and lateral phase separation into co-existing liquid-ordered and liquid-disordered domains, and validate it against experiments on lipid-domain coarsening and on two-phase vesicles in shear. Using this model, they claim that surface tension heterogeneity (the ld phase having lower surface tension) can by itself produce swinging and tumbling—two motions previously attributed to viscosity contrast or asymmetric phase distributions. They also map a phase diagram for bending-rigidity contrast and capillary number, identifying six regimes, three of which (lateral ring banding, lateral broken-ring banding, vertical broken-ring banding) are new. If true, the result turns membrane surface-tension contrast into a controllable design knob for vesicle-based carriers in flow.","feed_headline":"Surface tension contrast alone makes vesicles swing and tumble","feed_subtitle":"A validated 3D model maps new banding and tumbling states, giving a design knob for vesicle carriers.","key_machinery":"A phase-field model with two order parameters: φ marks the vesicle membrane as a diffuse interface, and c_Γ marks the lateral lo/ld composition; the diffuse-domain method extends surface equations into the bulk. The model is derived from a free energy (bending, surface, line, area-penalty) via energy variation, giving thermodynamically consistent forces and fluxes, with local inextensibility enforced by a Lagrange multiplier. The regime map is organized by dimensionless contrasts κ_B = κ_B^ld/κ_B^lo, σ_S = σ_S^ld/σ_S^lo, and the bending capillary number Ca = μ_out U L²/κ_B^lo. The surface-tension mechanism works by lowering the energetic cost of deforming ld domains, anchoring them at the ti","core_discovery":"The central discovery is twofold. First, for a prolate vesicle with reduced volume v≈0.91 and two symmetric ld domains pre-placed at the tips, lowering the ld/lo surface-tension ratio σS below 1 at fixed bending contrast κ_B=0.65 causes the vesicle to enter two new regimes: small-amplitude swinging (VII) and rigid-body-like tumbling (VIII), with no viscosity contrast and no asymmetric domain placement. Second, varying the bending-rigidity contrast κ_B and the bending capillary number Ca sweeps out six banding/treading regimes, including three not previously reported: lateral ring banding (I), lateral broken-ring banding (II), and vertical broken-ring banding (IV); within vertical ring bandin","pith_inferences":["If the surface-tension-only swinging and tumbling mechanism holds generally, it suggests vesicles could be steered by locally modulating tension (e.g., via drugs, light, or other stimuli) rather than by controlling viscosity; a testable extension is to scan reduced volume v and initial domain arrangement to check the robustness of regimes VII and VIII.","The surface-tension-driven swinging bears a formal resemblance to thermocapillary or surfactant-driven migration of droplets; the same diffuse-domain machinery might translate to droplets with heterogeneous surface tension.","The regime taxonomy depends on the diffuse-interface representation of surface tension (through the regularized delta function and the specific surface energy); a sharp-interface benchmark or a convergence study in the interfacial thickness ε_ϕ would clarify whether the new regimes are physical or artifacts of the model's regularization.","Because the phase diagram was computed from a single initial condition (prolate vesicle with v≈0.91 and symmetric tip domains), the regime boundaries are likely to shift with vesicle deflation and domain placement; mapping those shifts would produce a more complete design chart."],"forward_implications":["Surface-tension contrast should be treated as a first-order control parameter for vesicle motion in flow, alongside viscosity contrast and bending contrast.","The three new banding regimes (I, II, IV) expand the known dynamical repertoire and should be observable in experiments with tunable lipid compositions.","The identified 'Low-Rigidity Full-Dynamics Belt' and 'Low-Capillary Full-Dynamics Belt' mark parameter windows where all six regimes coexist, useful for designing vesicle-based carriers.","Subdividing vertical ring banding into temporary, cyclic, and stable rings implies the ring is a metastable structure governed by bending contrast and shear strength.","The model's agreement with experiments on phase separation and shear dynamics offers a validated 3D framework for predicting multicomponent vesicle behavior under flow."],"supporting_citations":[{"why":"Supplies the experimental snapshots and raft-count/perimeter data used to validate the phase-separation dynamics of the model.","marker":"[5]"},{"why":"Supplies the experimental two-phase vesicle shear-flow behaviors (tank treading, swinging, tumbling) used to validate the coupled vesicle-fluid model.","marker":"[33]"},{"why":"Reports previously known phase treading and 2D multicomponent vesicle dynamics that the new regimes (I, II, IV) extend to 3D.","marker":"[30]"},{"why":"Reported vertical ring banding in 3D multicomponent vesicles, which the model reproduces as regime (V) and then subdivides.","marker":"[31]"},{"why":"Reported swinging and tumbling in multicomponent vesicles, the baseline that the surface-tension-only mechanism (regimes VII and VIII) goes beyond.","marker":"[27]"},{"why":"Provides the diffuse-interface formulation of locally inextensible vesicles, including the Lagrange-multiplier inextensibility approach used here.","marker":"[36]"},{"why":"Supplies the block-structured adaptive multigrid solver (staggered finite differences with FAS) that carries out the simulations.","marker":"[48]"}],"fun_headline_variants":["Tension mismatch flips vesicles into swinging and tumbling","Surface tension heterogeneity triggers vesicle swinging and tumbling","New vesicle motion regimes from bending rigidity contrast","Vesicle dynamics: surface tension contrast as a control knob","Shear flow vesicles: tension contrast drives tumbling"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The whole regime map—including the claim that surface tension contrast alone triggers swinging and tumbling—is computed from a single initial condition: a prolate vesicle with reduced volume v≈0.91 and two symmetric ld domains placed at the tips; if the regimes depend sensitively on that geometry or on the diffuse-interface form of the surface-tension energy, the general claim would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Tension mismatch flips vesicles into swinging and tumbling","Surface tension heterogeneity triggers vesicle swinging and tumbling","New vesicle motion regimes from bending rigidity contrast","Vesicle dynamics: surface tension contrast as a control knob","Shear flow vesicles: tension contrast drives tumbling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000638,"raw_usage":{"total_tokens":2742,"prompt_tokens":679,"completion_tokens":2063,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":423,"completion_tokens_details":{"reasoning_tokens":1985}},"tokens_in":423,"tokens_out":2063,"duration_ms":20069,"temperature":1.0,"reasoning_tokens":1985,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T20:50:09.651131+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same parameters but replace the symmetric tip domains with a random or single off-center ld domain, or change the reduced volume to v=0.8; if surface-tension-driven swinging/tumbling disappears or the six-regime taxonomy collapses into different states, the mechanism is not generic. Alternatively, a microfluidic experiment on a giant unilamellar vesicle with matched bending rigidities but σ_S<1 should show the predicted tumbling at Ca≈1; failure to observe it would contradict the claim.","supporting_citations":[],"review_version":1}