{"id":"4553d044-c1c9-4313-b019-9c17656a4205","arxiv_id":"2501.05367","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Rotational drag on a spherical-cap inclusion in a spherical membrane is solved semi-analytically as a function of inclusion size and Saffman-Delbrück length, revealing five mobility regimes.","lead":"This paper computes how fast a rigid patch of a spherical cell membrane, such as a protein or lipid domain, can rotate when the membrane sits in a viscous fluid. The result gives a new theoretical tool for measuring membrane viscosity from the rotation of large inclusions, where the classic flat-membrane formula fails.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Finite-η Legendre truncation is not checked; the numerical Λ_R values and regime boundaries in Figs 5-6 could shift if the k=300 solution is not converged at ε=O(1).","rationale":"The analytical core is solid: the reduction of torque to c1 (Eq. 2.21) follows from orthogonality, and the matched asymptotics for regions (1), (3), and (5) are self-consistent and reproduce the planar [30] and spherical-cap [48] results. The numerical truncation at finite ε is the weakest load-bearing point because Figs. 5-6 and the experimental proposal depend on the quantitative values of Λ_R, not just on the asymptotic scaling. A convergence study at fixed finite ε with higher k would resolve this. Separately, the manuscript contains several pages of extraneous text from other papers (colloidosomes, DNA-membrane experiments, liquid-domain diffusion) that must be removed before publication; this does not affect the hydrodynamic argument but is a serious manuscript hygiene issue. The absence of code and data is secondary. The reader's conditional verdict remains appropriate pending the finite-ε convergence check.","tokens_in":22579,"tokens_out":45060,"duration_ms":444876,"concrete_test":"For θp=0.1 and ε ∈ {10^{-4}, 10^{-3}, 10^{-2}, 0.1, 1, 10, 100}, recompute Λ_R with Legendre truncations k=400 and k=600 (and k=800 for a subset) using the same ODE solver; if the maximum relative change from k=300 exceeds 2% at any finite ε, the finite-η numerical results in Figs. 5-6 and the inferred regime boundaries are not converged. Cross-check the k→∞ extrapolation of c1 against the analytical asymptotics in Table 2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative form of the central claim—the values of Λ_R(ε,θp) in Figs 5-6 and the locations of the five regimes—rests on the semi-analytical solution truncated at k=100-300 modes. The only convergence check reported is in the limit η→0 (Section 3(a)), where the flow reduces to a rigid spherical cap in an unbounded solvent; no truncation or residual error estimate is given for the finite ε values plotted (ε=0.1,1,10 in Fig 6A; ε=0.1-100 in Fig 5A). The truncation error can depend on ε: at intermediate ε the membrane velocity has a boundary layer of width O(θ_p^2) near the cap edge (Eqs. 3.3-3.4), requiring many Legendre modes to resolve, while the singular corner is smoothed by membrane viscosity, so it is not obvious that the η→0 check is the worst case. If the error grows at finite ε, the reported Λ_R values and the inferred boundaries of regions (2) and (4) shift, directly affecting the proposed use of rotational diffusion to measure membrane viscosity. The asymptotic scalings in regions (1), (3), and (5) are derived analytically and are robust; the concern is the numerical support for the quantitative crossover and prefactors.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a semi-analytical theory for the rotational mobility of a rigid spherical-cap inclusion embedded in a spherical incompressible membrane surrounded by Stokes solvents. The membrane and solvent flows are expanded in a Legendre-derived rotlet multipole basis, reducing the problem to a coupled system, Eqs. (2.16)-(2.19), for the expansion coefficients and the membrane angular velocity. The total torque is shown to be controlled solely by the rotlet coefficient c1, Eq. (2.21), so the dimensionless drag Lambda_R(epsilon, theta_p) depends only on the Saffman-Delbruck ratio epsilon=2 mu R_p/eta and the cap half-angle theta_p=R_p/R_m. The paper computes Lambda_R numerically for a range of parameters and identifies five asymptotic regimes for small inclusions, namely rigid whole-vesicle rotation, a crossover, two-dimensional membrane drag, the planar Saffman-Delbruck regime, and three-dimensional solvent-dominated drag, together with a large-inclusion limit that recovers the spherical-cap result of Dorrepaal.","tokens_in":22839,"tokens_out":16627,"duration_ms":166430,"significance":"If the numerical results are converged, this is a substantive contribution. It gives a parameter-free model for rotational drag on finite-size inclusions in spherical vesicles, includes curvature effects that are absent from planar Saffman-Delbruck theory, and offers a concrete route to measuring membrane viscosity from rotational diffusion of large inclusions. The clean reduction of the torque to the rotlet coefficient c1 and the analytic derivation of the asymptotic scalings are clear strengths, and the comparisons with the planar result of Ref. [30] and the spherical-cap result of Ref. [48] are appropriate external benchmarks. The main outstanding issue is numerical: the truncation accuracy is verified only in the eta->0 limit, so the quantitative values in Figs. 5 and 6 and the inferred crossover positions need confirmation before the quantitative claims can be accepted.","major_comments":[{"comment":"The stated truncation at k=100-300 modes is validated only in the eta->0 limit, where the spherical-cap mobility [48] is recovered with a relative error no larger than 2%; no convergence test or error estimate is reported for the finite values of epsilon plotted in Fig. 5A (epsilon=0.1-100) and Fig. 6A (epsilon=0.1, 1, 10). Since the quantitative positions of the crossovers in regions (2) and (4) and the values of Lambda_R in Figs. 5 and 6 are numerical, this is a load-bearing gap: if the Legendre truncation error grows at finite epsilon, the regime boundaries and the proposed viscosity-measurement method would shift. Please add a convergence study at representative finite epsilon, for example at theta_p=0.1 with epsilon=0.1, 1, and 10 and k=300, 600, and 1200, and report either the relative change in Lambda_R or an a posteriori residual of Eqs. (2.16)-(2.18).","section":"Section 3(a), Figs. 5 and 6"},{"comment":"The five-regime picture for small inclusions is tested numerically only at theta_p=0.1. While this is a small cap, it is not asymptotically small for all purposes, and the planar-limit comparison in region (4) and the crossover between regions (3) and (5) would be considerably more convincing if a second, smaller value, such as theta_p=0.05, were shown or if the theta_p-dependence of Lambda_R were demonstrated to have converged. Without such evidence, the stated universal boundaries, which are obtained by scaling rather than by exact asymptotics, may not yet represent the theta_p->0 limit accurately.","section":"Section 3(b)(i), Fig. 6A"}],"minor_comments":[{"comment":"The manuscript contains several paragraphs that appear to be reprinted from Refs. [6], [7], and [39], including passages titled 'Unbinding of DNA from Cationic Membranes', 'Coverage of the emulsified droplets by colloidal particles', and 'Diffusion of Liquid Domains in Lipid Bilayer Membranes'. These paragraphs do not belong to the paper's narrative and should be removed before submission; as printed, they disrupt the introduction and create a provenance concern.","section":"Section 1, after Fig. 1"},{"comment":"Equation (A7) has a spurious trailing '=0' after the explicit non-zero traction sum; this is inconsistent with Eq. (2.17) and with the preceding line (A6), and should be removed.","section":"Appendix C, Eq. (A7)"},{"comment":"The region-(1) boundary condition is written inconsistently: the text gives epsilon << theta_p^3, which is equivalent to l_SD >> R_m^3/R_p^2, while Table 2 and the following sentence write l_SD >> R_m^3/R_p^3. The exponents should be harmonized, with R_m^3/R_p^2 being the version consistent with the text.","section":"Section 3(b)(i) and Table 2"},{"comment":"The phrase 'rotlet 2n-pole' should presumably be 'rotlet 2^n-pole', since the n-th derivative of the rotlet produces a 2^n-pole singularity.","section":"Table 1 caption"},{"comment":"The phrase 'the solid angle formed by the particle' is used to mean the polar half-angle theta_p; consider using 'solid angle' only for the actual quantity 2*pi*(1-cos(theta_p)), and referring to theta_p as the cap half-angle.","section":"Abstract and Section 2(a)"},{"comment":"Several reference entries have questionable or mismatched DOIs; for example, Ref. [14] gives a DOI from a Biophysical Journal article that does not appear to correspond to Brochard and Lennon 1975. The reference list should be checked systematically.","section":"References"},{"comment":"The sentence 'These values of k ensure that, in the limit eta->0, we recover the mobility of a rotating spherical cap with a relative error no larger than 2% (see § ii)' points to Section 3(b)(ii), which is the large-particle asymptotic section, not the truncation check; the pointer should be to Eq. (3.7) and Fig. 6C, where the spherical-cap comparison is made.","section":"Section 3(a)"}],"recommendation":"major_revision","confidential_remarks":"The main scientific hurdle is the missing finite-epsilon convergence study; if the authors supply it and it confirms the plotted values, the paper should be suitable for publication. Also, the extraneous reference-manager text inserted in Section 1 should be cleaned before acceptance; it is a manuscript-integrity issue separate from the scientific content. The paper appears within the journal's scope, and the novelty disclosure is adequate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a real contribution, not a repackaging. The paper solves the rotational drag on a rigid spherical cap in a spherical membrane coupled to Stokes solvents, treating cap size and Saffman-Delbrück length as two independent parameters. The reduction of the torque to the rotlet coefficient c1 is clean (Eqs. 2.20-2.21), the governing equations are derived carefully from surface Stokes balances, and the known limits all check out: rigid-body rotation for large Saffman length, the ε^{-1} two-dimensional drag, the planar Hughes-Pailthorpe-White result, and the Dorrepaal spherical-cap torque in the solvent-dominated limit. The five-regime asymptotic map for small inclusions is the most useful output, and the analytic scalings in regimes (1), (3), and (5) do not depend on the numerics. That is solid, careful work.\n\nThe soft spots are proportionate. The stress-test concern is fair: the 2% convergence check is only reported for η→0, and the finite-ε truncation error is not assessed. The boundary-layer structure near the cap edge at intermediate ε suggests the k=300 truncation could in principle be less accurate there than in the solvent-only limit, so the quantitative crossover positions and prefactors in regions (2) and (4) rest on an unverified numerical assumption. I would not call this fatal, because the main physical conclusions would not change if those prefactors moved by a few percent, but the authors should either provide a convergence study at finite ε or release code so readers can check. Second, no code or data is included, which is a practical barrier for a semi-analytical method that others will want to use. Third, there is a page of unrelated text from another paper inserted into the arXiv PDF (the DNA/colloidosome material around section 1). That is plainly a production error and should be fixed before any publication, but it does not affect the science.\n\nWho is this for? Anyone working on membrane hydrodynamics, especially experimentalists wanting to extract membrane viscosity from rotational diffusion of large inclusions. The paper gives them a concrete prediction and a clear warning about where planar Saffman-Delbrück theory fails. I would cite it if I were writing on curved-membrane mobility. It deserves a serious referee: the derivation is original, the asymptotics are informative, and the numerical concern is addressable. My recommendation is to send it to review, with the request that the authors add a finite-ε convergence check and share their code.","headline":"A genuinely new finite-size curved-membrane mobility calculation with a clean derivation and a useful five-regime map; the main weakness is an unverified truncation error at finite membrane viscosity, plus an obvious manuscript-production glitch.","tokens_in":23367,"tokens_out":1100,"would_cite":true,"duration_ms":13092,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76D07","76Z99"],"pacs":[],"model":"deepseek-v4-flash","headline":"A rotating cap-shaped inclusion in a spherical membrane has its entire rotational drag set by the first rotlet coefficient, so the dimensionless mobility depends only on inclusion size and the Saffman–Delbrück ratio.","keywords":["membrane hydrodynamics","rotational mobility","Saffman–Delbrück length","spherical membrane","spherical-cap inclusion","rotlet","Legendre expansion","vesicle"],"falsifier":"Recompute $\\Lambda_R(\\varepsilon,\\theta_p)$ at a fixed finite $\\varepsilon$ (for example $\\varepsilon=1$, $\\theta_p=0.1$) with increasing truncation $k=100,200,300,600$; if the values do not converge to within a few percent, the claimed mobility curve and the five-regime boundaries are not established. A complementary check would be to measure the rotational diffusivity of a micrometer-scale domain in a giant unilamellar vesicle with an independently known membrane viscosity and compare the inferred $\\Lambda_R$ with Fig. 5.","tokens_in":22344,"feed_emoji":"🔄","tokens_out":9433,"duration_ms":85207,"temperature":0.7,"pith_summary":"This paper computes how hard it is to rotate a rigid cap-shaped inclusion embedded in a spherical fluid membrane surrounded by viscous solvent on both sides. It claims that the entire rotational drag is captured by a single Legendre coefficient, the rotlet strength $c_1$, so the dimensionless drag $\\Lambda_R(\\varepsilon,\\theta_p)=\\Omega^{-1}\\theta_p^{-3}c_1$ depends only on the inclusion's angular size $\\theta_p=R_p/R_m$ and the Saffman–Delbrück ratio $\\varepsilon=2\\mu R_p/\\eta$. As the Saffman–Delbrück length is varied, a small cap passes through five asymptotic drag regimes, from rigid whole-vesicle rotation through nearly planar membrane drag to bare three-dimensional solvent drag. The flat Saffman–Delbrück result is recovered only in the intermediate window $R_p \\ll \\ell_{\\mathrm{SD}} \\ll R_m^3/R_p^2$, and the paper identifies when planar theory fails for large inclusions or highly viscous membranes. If correct, this provides a route to measuring membrane viscosity from the rotational diffusion of large inclusions.","feed_headline":"Five drag regimes for spinning caps inside curved membranes","feed_subtitle":"A single rotlet coefficient predicts how membrane and solvent viscosity set the torque on a cap-shaped particle.","key_machinery":"The load-bearing object is the rotlet coefficient $c_1$, the first term in a Legendre/rotlet-multipole expansion of the fluid velocity in the outer solvent, inner solvent, and membrane (Eqs. 2.12–2.14). Orthogonality of the Legendre derivatives $P_n'$ reduces the momentum balance to a linear recurrence for the coefficients coupled to an ordinary differential equation for the membrane angular velocity $v(x)$ (Eqs. 2.16–2.19); truncating at $k$ modes gives a finite linear system whose solution yields $c_1$. The identity $G_p=8\\pi\\mu R_m^3 c_1$ (Eq. 2.21) is what carries the argument, because it shows that all details of the flow and inclusion shape enter the torque only through this single coefficient. The rest of the machinery turns the calculation of $c_1$ into a two-parameter problem depending only on $\\varepsilon$ and $\\theta_p$.","core_discovery":"The paper's central claim is that for a spherical membrane of radius $R_m$ containing a rigid spherical cap of half-angle $\\theta_p$ rotating with angular velocity $\\Omega$, the torque needed to maintain rotation is $G_p = 8\\pi\\mu R_m^3 c_1$, where $c_1$ is the rotlet (first Legendre) coefficient of the surrounding Stokes flow (Eq. 2.21). All higher-order multipoles reshape the flow but contribute nothing to the torque, so the dimensionless drag coefficient $\\Lambda_R(\\varepsilon,\\theta_p)=\\Omega^{-1}\\theta_p^{-3}c_1$ is a function of only two parameters: $\\varepsilon=2\\mu R_p/\\eta$ and $\\theta_p=R_p/R_m$. In the small-particle limit the paper identifies five asymptotic regimes as $\\ell_{\\mathrm{SD}}=\\eta/\\mu$ grows: rigid rotation of the whole vesicle, a transition region, solvent-free two-dimensional disk drag, the planar Saffman–Delbrück regime, and finally three-dimensional disk drag when the membrane effectively disappears. The planar Saffman–Delbrück mobility and the exact rotational drag of a spherical cap in an unbounded fluid are recovered as special limits of the same calculation. The result is obtained by solving the coupled membrane–solvent Stokes equations semi-analytically with truncated Legendre expansions.","pith_inferences":["A natural extension, which the paper leaves for future work, is the translational analogue: because purely azimuthal flow on a sphere is automatically divergence-free, the regular $\\eta\\to0$ limit found here is special to rotation, and translational mobility should still feel membrane incompressibility and curvature.","The unsteady version needed for rotational Brownian motion of liquid domains, where random molecular torques act on fast timescales, may not obey the same five-regime ordering because oscillatory flows introduce a new length scale into the membrane–solvent coupling.","Because the torque depends only on $c_1$, a measurement of rotational drag (or rotational diffusion) gives a direct handle on the rotlet amplitude, potentially allowing the prediction of the far-field solvent flow without reconstructing the full membrane velocity field.","The same Legendre machinery could be applied to asymmetric inclusions on a sphere by using vector spherical harmonics, though the loss of axisymmetry would replace the scalar ODE with a more complex coupled system."],"forward_implications":["For a small inclusion in a vesicle, planar Saffman–Delbrück theory is valid only while $\\ell_{\\mathrm{SD}} \\ll R_m^3/R_p^2$; outside this window the spherical geometry or bare solvent drag dominates and planar fits misestimate membrane viscosity.","Rotational diffusion measurements of large inclusions, where planar theory fails, can be converted into a membrane-viscosity estimate using the computed relation between $\\Lambda_R$ and $\\varepsilon$.","In the limit of vanishing membrane viscosity the spherical membrane imposes no constraint on purely azimuthal rotational flow, so it simply disappears and the drag tends to that of a rotating spherical cap in an unbounded fluid.","For very viscous membranes the whole vesicle rotates almost rigidly with the cap, producing a torque proportional to $R_m^3$ and a drag coefficient that scales as $\\theta_p^{-3}$, a regime absent from flat-membrane theories.","The five-regime sequence is controlled by two dimensionless parameters, so a single plot such as Fig. 5 can be used to read off the expected rotational drag for any vesicle radius, inclusion size, and viscosity ratio."],"supporting_citations":[{"why":"supplies the Saffman–Delbrück model of a membrane as a 2D viscous fluid coupled to Stokesian solvents and defines the length scale $\\ell_{\\mathrm{SD}}=\\eta/\\mu$.","marker":"[1]"},{"why":"introduces the curved-membrane mobility setting and the co-rotating drag definition that the paper extends to finite-size inclusions.","marker":"[26]"},{"why":"provides the planar rotational and translational drag of a cylinder in a membrane, the baseline recovered in region (4) and used for comparison in Fig. 6A.","marker":"[30]"},{"why":"gives the three-dimensional drag of a rotating disk, the solvent-dominated limit $\\Lambda_R\\sim 4/3\\pi$ in region (5).","marker":"[27]"},{"why":"provides the exact rotational Stokes resistance of a spherical cap in an unbounded fluid, the $\\eta\\to0$ limit of the large-particle analysis.","marker":"[48]"},{"why":"supplies experimental giant-unilamellar-vesicle parameters (domain sizes and viscosity range) used to choose representative values of $\\varepsilon$ and $\\theta_p$.","marker":"[39]"},{"why":"establishes the thin-sheet Brownian-motion drag used in the two-dimensional membrane-dominated regime.","marker":"[28]"}],"fun_headline_variants":["Five drag regimes for rotors in spherical membranes","Torque on a spinning cap predicted by one coefficient","Five regimes from flat to curved membrane drag","SD length and cap size set torque in curved membranes","Spherical cap rotation mapped by SD length and size"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical truncated-Legendre solution is trusted at the finite membrane viscosities used in the results, even though the paper's stated 2% accuracy check (Section 3(a), with $k=100$–$300$ modes) was performed only in the limit $\\eta\\to0$; if the truncation error grows at finite $\\eta$, the reported drag values and the inferred boundaries between the five regimes would shift.","fun_headline_variants_meta":{"raw":{"variants":["Five drag regimes for rotors in spherical membranes","Torque on a spinning cap predicted by one coefficient","Five regimes from flat to curved membrane drag","SD length and cap size set torque in curved membranes","Spherical cap rotation mapped by SD length and size"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000638,"raw_usage":{"total_tokens":2987,"prompt_tokens":1041,"completion_tokens":1946,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":1873}},"tokens_in":657,"tokens_out":1946,"duration_ms":13190,"temperature":1.0,"reasoning_tokens":1873,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:14:59.541198+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute $\\Lambda_R(\\varepsilon,\\theta_p)$ at a fixed finite $\\varepsilon$ (for example $\\varepsilon=1$, $\\theta_p=0.1$) with increasing truncation $k=100,200,300,600$; if the values do not converge to within a few percent, the claimed mobility curve and the five-regime boundaries are not established. A complementary check would be to measure the rotational diffusivity of a micrometer-scale domain in a giant unilamellar vesicle with an independently known membrane viscosity and compare the inferred $\\Lambda_R$ with Fig. 5.","supporting_citations":[{"cited_title":"1975 Brownian motion in biological membranes","cited_arxiv_id":null,"evidence_quote":"supplies the Saffman–Delbrück model of a membrane as a 2D viscous fluid coupled to Stokesian solvents and defines the length scale $\\ell_{\\mathrm{SD}}=\\eta/\\mu$."},{"cited_title":"2008 The effect of curvature and topology on membrane hydrodynamics","cited_arxiv_id":null,"evidence_quote":"introduces the curved-membrane mobility setting and the co-rotating drag definition that the paper extends to finite-size inclusions."},{"cited_title":"1981 The translational and rotational drag on a cylinder moving in a membrane","cited_arxiv_id":null,"evidence_quote":"provides the planar rotational and translational drag of a cylinder in a membrane, the baseline recovered in region (4) and used for comparison in Fig. 6A."},{"cited_title":"1965 Low Reynolds Number Hydrodynamics","cited_arxiv_id":null,"evidence_quote":"gives the three-dimensional drag of a rotating disk, the solvent-dominated limit $\\Lambda_R\\sim 4/3\\pi$ in region (5)."},{"cited_title":"1978 The Stokes resistance of a spherical cap to translational and rotational motions in a linear shear flow","cited_arxiv_id":null,"evidence_quote":"provides the exact rotational Stokes resistance of a spherical cap in an unbounded fluid, the $\\eta\\to0$ limit of the large-particle analysis."},{"cited_title":"2007 Diffusion of liquid domains in lipid bilayer membranes","cited_arxiv_id":null,"evidence_quote":"supplies experimental giant-unilamellar-vesicle parameters (domain sizes and viscosity range) used to choose representative values of $\\varepsilon$ and $\\theta_p$."},{"cited_title":"1976 Brownian motion in thin sheets of viscous fluid","cited_arxiv_id":null,"evidence_quote":"establishes the thin-sheet Brownian-motion drag used in the two-dimensional membrane-dominated regime."}],"review_version":1}