{"id":"b0c0268e-7d42-41f6-924a-928c8722fe9f","arxiv_id":"1909.02066","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"An interfacial regularized Stokeslet method computes advection, rotation, and effective viscosity of arbitrarily shaped rigid inclusions in lipid membranes, reproducing known limits and predicting higher viscosity for linear oligomers.","lead":"This paper presents a numerical method to predict how proteins and gel domains drift and rotate inside flowing lipid membranes, and to estimate how much they thicken the membrane. The method works for arbitrarily shaped objects, where simple analytical formulas fail, and suggests that measuring membrane viscosity could reveal whether proteins cluster into chains.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified","rationale":"The reader's weakest assumption points to the dilute-limit multipole/correlation truncation in Appendix A. That is also the most natural place to probe the argument, but it is not a load-bearing concern: the paper explicitly restricts claims to the dilute limit, and the standard first-order-in-phi effective-viscosity derivation is precisely a single-particle stresslet calculation with correlations and higher multipoles relegated to higher order. The RS method itself is validated against known analytical limits, and the new oligomer results are supported by an independent analytical-orientation-averaging route in Appendix B with very small numerical difference. The remaining caveats, such as the unexplained 2/3 factor relative to Henle-Levine and the absence of error bars in Fig. 5, are acknowledged in the manuscript and do not undermine the central conclusions. Therefore the reader's ACCEPT verdict should stand unchanged.","tokens_in":16065,"tokens_out":22142,"duration_ms":269081,"concrete_test":"Run the public code for the a/L_sd = 0.1 cylinder case with blob spacings s = 0.025a and s = 0.0125a (epsilon = s/2), re-extrapolating alpha from three-point fits, and recompute the N = 4 oligomer point with 40 orientation samples instead of 10. If either alpha shifts by more than about 1% from Figs. 4 and 5, the numerical-convergence claims need revision; otherwise, the conclusions stand.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After reading the derivation and numerics, I do not find a load-bearing defect in the central claim. The RS constraint formulation is validated by the stationary-cylinder test and the comparison against the Oppenheimer-Diamant Faxén relations. The effective-viscosity claim rests on the dilute-limit assumptions in Appendix A: single-particle stresslet, no particle-particle correlations, and uniform orientation averaging. These are explicitly stated and are the standard first-order-in-phi conditions for an intrinsic-viscosity calculation; the dipole truncation is appropriate for the long-wavelength, dilute response that defines eta_eff^m. The unexplained 2/3 factor relative to Henle-Levine is definitional and is acknowledged as such, and the lack of error bars in Fig. 5 does not threaten the OD-matching cylinder result or the qualitative oligomer trend. The internal grand-resistance-matrix cross-check in Appendix B agrees to 2e-5, providing independent support.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the interfacial regularized Stokeslet (RS) method developed by the authors to compute the translational and rotational velocities of force- and torque-free solid objects embedded in a lipid bilayer membrane and advected by an ambient surface flow (Eqs. 3–5). The method is validated by recovering, in the limit a ≪ L_sd, the approximate Faxén relations of Oppenheimer and Diamant (Eqs. 11–12) for cylinders and by the stationary-cylinder test in extensional flow. The same stresslet machinery is used to compute the dilute-limit effective membrane viscosity η_eff^m = η_m(1 + αφ) from the orientationally averaged single-particle force dipole (Appendix A). For cylinders, α → 2 as a/L_sd → 0, consistent with Oppenheimer and Diamant; α grows with a/L_sd, and with a rescaling factor of 2/3 the results agree with Henle and Levine. For rigid linear oligomers of circular monomers, α increases with oligomer length. The grand resistance matrix formulation in Appendix B provides an independent cross-check, agreeing with the direct stresslet calculation to 2×10^-5.","tokens_in":16232,"tokens_out":10161,"duration_ms":91319,"significance":"If the results hold, the paper provides a general numerical tool for membrane hydrodynamics beyond simple cylinders, applicable to arbitrary shapes and spatially varying flows. It gives a route to compute membrane Einstein corrections for non-circular inclusions, and predicts a measurable dependence of membrane effective viscosity on protein oligomerization state. Strengths include: the method is validated against the known Oppenheimer-Diamant limit; the effective-viscosity calculation rests on explicitly stated dilute-limit assumptions (single-particle stresslet, no correlations) that are the standard first-order-in-φ conditions; an independent grand-resistance-matrix cross-check agrees to 2×10^-5; and the code is publicly available. The paper is honest about limitations (range of area fractions not quantified, 2/3 rescaling not understood).","major_comments":[],"minor_comments":[{"comment":"The statement that the Henle-Levine predictions are 'confirmed' is stronger than the evidence supports, given that agreement requires an unexplained factor of 2/3; I suggest rewording to 'reproduced up to an overall factor of 2/3'.","section":"Section IV / Discussion"},{"comment":"There is a typo in the sentence 'extending the original Oppenheimer Diamant calculation to determine the the effective viscosity'; delete the duplicated 'the'.","section":"Section IV"},{"comment":"The oligomer intrinsic-viscosity results in Fig. 5 are presented without error bars; adding uncertainties from the spacing extrapolation and orientation averaging would help the reader judge the significance of the increase with oligomer length.","section":"Figure 5"},{"comment":"The convergence study is described only in figure captions; a representative plot or table showing the extrapolation to zero blob spacing for one test case would strengthen the numerical claims.","section":"Section II / Figures 2–5"},{"comment":"The statement 'The origin of the factor of 2/3 is not understood' (Section IV) is in tension with the Discussion's remark that this factor 'was previously attributed to the different definitions of effective viscosity'; please clarify whether a definitional explanation is accepted.","section":"Section IV / Discussion"},{"comment":"In Eq. (B7), the notation '(G H)' is not defined; writing the block matrices explicitly would improve readability.","section":"Appendix B, Eq. (B7)"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid contribution to membrane hydrodynamics. The only substantive reservation is the unexplained 2/3 rescaling in the Henle-Levine comparison, but this does not affect the central results, which are independently validated by the Oppenheimer-Diamant limit and the grand-resistance-matrix cross-check. The manuscript is suitable for publication after minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth your time if you care about membrane hydrodynamics. This is a clean incremental extension of Camley and Brown's earlier interfacial regularized Stokeslet work, and it delivers what it promises: a numerical route to compute advection, rotation, and effective viscosity for arbitrarily shaped membrane inclusions. The new piece is the force- and torque-free formulation in an ambient flow, plus the stresslet-based intrinsic viscosity calculation. The validation is the strong part: they recover the Oppenheimer-Diamant alpha=2 limit for small cylinders, and the grand-resistance-matrix cross-check in Appendix B agrees to 2e-5. Code is on GitHub, which makes the numerics reproducible rather than just asserted.\n\nSoft spots are real but not load-bearing. The 2/3 rescaling of Henle-Levine is left unexplained; the authors say it's definitional, and the comparison is honest, but a referee should push for a clearer statement about why a simple factor should collapse the two theories across the whole range. The dilute-limit derivation in Appendix A assumes single-particle stresslets, no correlations, and orientation averaging; those are standard first-order-in-phi assumptions, and the paper states them explicitly. Fig. 5 lacks error bars, which is minor given the convergence extrapolation described. The advection results are interesting but not dramatic; the deviations from Faxen are modest unless the flow varies on the particle scale, which is fine.\n\nOverall the physics looks sound. The intrinsic viscosity framework is the most valuable contribution, especially the oligomer trend that suggests a possible experimental handle on oligomerization state. This is not a breakthrough, but it is a competent, useful paper that fills a real gap. I would send it to peer review rather than desk reject, and I'd expect it to be accepted after minor revision.","headline":"A solid, useful extension of the authors' regularized Stokeslet method; the effective-viscosity results are credible and the paper deserves a serious referee.","tokens_in":16713,"tokens_out":956,"would_cite":true,"duration_ms":11893,"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":"An interfacial regularized Stokeslet scheme predicts the motion of arbitrary rigid bodies in flowing lipid-bilayer membranes and yields the dilute-limit effective viscosity $\\eta_m(1+\\alpha\\phi)$.","keywords":["lipid bilayer membranes","Saffman-Delbrück model","regularized Stokeslets","effective membrane viscosity","intrinsic viscosity","membrane hydrodynamics","Faxén relationships","rigid oligomers"],"falsifier":"A many-body numerical experiment—randomly dispersing rigid disks or rods in a sheared membrane at several small area fractions and measuring how the averaged stresslet and effective viscosity grow with $\\phi$—would settle whether the single-particle dipole law holds; if the slope $\\alpha$ changes with $\\phi$ or with multipole truncation, the derivation fails.","tokens_in":15907,"feed_emoji":"🧬","tokens_out":11071,"duration_ms":109044,"temperature":0.7,"pith_summary":"This paper extends the interfacial regularized Stokeslet method, which previously computed drag and diffusion of membrane-embedded objects, to predict how such objects are advected and rotated by an externally imposed membrane flow. The central claim is that the same calculation yields the dilute-limit effective membrane viscosity, $\\eta_{\\rm eff}^m = \\eta_m(1+\\alpha\\phi)$, with the intrinsic viscosity $\\alpha$ obtained from the orientationally averaged force dipole (stresslet) exerted by a single object. In the small-object limit the method confirms the analytical prediction that circular inclusions give $\\alpha\\to 2$, and it extends the calculation to arbitrary shapes and to rigid linear oligomers, whose intrinsic viscosity increases with chain length. A sympathetic reader would care because this provides a numerical route to compute protein motion in complex flows where approximate Faxén relations fail, and it suggests that membrane viscosity measurements could report on protein oligomerization state.","feed_headline":"Rigid protein chains raise membrane viscosity more than monomers","feed_subtitle":"Interfacial regularized Stokeslets predict advection, rotation, and dilute-limit viscosity of lipid-bilayer inclusions.","key_machinery":"The central object is the interfacial regularized Stokeslet: the Saffman-Delbrück membrane's Oseen tensor, regularized by replacing point forces with Gaussian 'blobs' so that a solid body can be discretized as a cluster of constrained fluid regions. The argument is carried by the linear system that enforces rigid-body motion at every blob, the vanishing of total constraint force and torque, and the subsequent construction of the stresslet from the blob forces. For anisotropic objects, the intrinsic viscosity is obtained from the orientational average of the stresslet, computed either by repeated rotations or analytically through the grand resistance matrix, which also provides an order-of-magnitude speedup.","core_discovery":"The paper's central discovery is that force- and torque-free rigid bodies embedded in a flowing bilayer can be fully characterized by solving a linear system for the constraint forces on a cluster of regularized blobs: the membrane velocity at each blob must equal $\\mathbf{U} + \\boldsymbol{\\Omega}\\times\\mathbf{R}_m$, while the sum of blob forces and torques vanishes. The solution gives the body's translational velocity $\\mathbf{U}$ and angular velocity $\\boldsymbol{\\Omega}$ in any ambient field, reducing to the known Faxén relations when the flow is smooth and the body is small compared to the Saffman-Delbrück length, and deviating from them when the flow varies on the scale of the body. Under a pure shear, the same forces define the stresslet $S_{ij} = \\frac12\\sum_n (R_{n,i} g_j[\\mathbf{R}_n] + R_{n,j} g_i[\\mathbf{R}_n])$, whose orientational average satisfies $S_{ij} = -\\alpha\\eta_m A_p(\\partial_i v_j + \\partial_j v_i)$; inserting this into the averaged membrane response yields the Einstein-type correction $\\eta_{\\rm eff}^m = \\eta_m(1+\\alpha\\phi)$. Numerically, $\\alpha\\to 2$ for small cylinders, matching the analytical result of Ref. 6, and grows with $a/L_{\\rm sd}$; rigid linear oligomers display larger $\\alpha$ than monomers at the same area fraction.","pith_inferences":["The persistence of the unexplained 2/3 factor between the two prior analytical predictions across all $a/L_{\\rm sd}$ suggests it is a systematic convention difference in how the induced force dipole is defined, not a small-particle artifact; identifying that convention would reconcile the two definitions.","Because the method computes the full blob force distribution, the same machinery could handle externally forced or self-propelled inclusions by relaxing the zero-force/zero-torque constraints, and finite-concentration systems by adding pair interactions and higher multipoles—extensions the paper does not make.","If chain-length dependence of $\\alpha$ is robust, membrane shear-viscosity measurements in systems with controlled protein clustering could serve as a shape-sensitive probe of oligomerization state once the experimental precision the paper calls for is reached.","The paper's caveat that elongated particles may leave the linear regime at smaller area fractions than compact ones is testable: a finite-$\\phi$ simulation of sheared rigid rods should show where $\\eta_{\\rm eff}^m/\\eta_m - 1$ becomes nonlinear in $\\phi$."],"forward_implications":["In smooth, slowly varying ambient flows the scheme reproduces the approximate Faxén relations, so it provides a quantitative check on when those approximations are valid.","In rapidly varying flows, the Faxén truncation can produce unphysical oscillatory trajectories, while the regularized-Stokeslet trajectory remains smooth, so particle paths can be computed reliably in complex flow fields.","For dilute suspensions of small cylindrical inclusions the effective membrane viscosity is $\\eta_m(1+2\\phi)$, in agreement with the analytical prediction of Ref. 6.","For rigid linear oligomers at fixed area fraction, the intrinsic viscosity increases with chain length, with a stronger relative effect when $a/L_{\\rm sd}$ is larger.","The grand resistance matrix formulation permits analytic orientational averaging, reducing the cost of the effective-viscosity computation by roughly an order of magnitude."],"supporting_citations":[{"why":"Supplies the approximate Faxén relations and the force-dipole derivation of effective viscosity that this paper extends numerically and confirms in the small-particle limit.","marker":"[6]"},{"why":"Gives the alternative analytical effective-viscosity prediction that the numerics reproduce after a 2/3 rescaling across the full range of $a/L_{\\rm sd}$.","marker":"[26]"},{"why":"Introduced the interfacial regularized Stokeslet method for drag and diffusion of arbitrary membrane-embedded shapes, which this paper extends to advection and viscosity.","marker":"[24]"},{"why":"Supplies the regularized-Stokeslet Green's function formulas and implementation details used in the discretization.","marker":"[25]"},{"why":"Establishes the method of regularized Stokeslets in three dimensions, the representation of solid bodies as constrained fluid regions on which the scheme relies.","marker":"[29]"},{"why":"Provides the grand resistance matrix formalism used in Appendix B to compute the orientationally averaged stresslet analytically.","marker":"[36]"},{"why":"Supports the expectation that the linear effective-viscosity correction remains reasonable up to area fractions around 0.1, which underlies the dilute-limit interpretation.","marker":"[37]"}],"fun_headline_variants":["Rigid membrane inclusions: advection, rotation, viscosity","Chains beat monomers in raising bilayer viscosity","Regularized Stokeslets predict motion of membrane objects","How rigid bodies drift, spin, and thicken bilayers","Rod-shaped inclusions boost membrane viscosity more"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything about the effective-viscosity formula rests on the dilute-limit assumption that many inclusions act only as independent single-particle force dipoles, with no higher multipoles and no particle-particle correlations; if those contribute at the area fractions of interest, $\\eta_{\\rm eff}^m = \\eta_m(1+\\alpha\\phi)$ would be inaccurate.","fun_headline_variants_meta":{"raw":{"variants":["Rigid membrane inclusions: advection, rotation, viscosity","Chains beat monomers in raising bilayer viscosity","Regularized Stokeslets predict motion of membrane objects","How rigid bodies drift, spin, and thicken bilayers","Rod-shaped inclusions boost membrane viscosity more"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000605,"raw_usage":{"total_tokens":2824,"prompt_tokens":948,"completion_tokens":1876,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":1802}},"tokens_in":564,"tokens_out":1876,"duration_ms":12420,"temperature":1.0,"reasoning_tokens":1802,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:00:48.662476+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A many-body numerical experiment—randomly dispersing rigid disks or rods in a sheared membrane at several small area fractions and measuring how the averaged stresslet and effective viscosity grow with $\\phi$—would settle whether the single-particle dipole law holds; if the slope $\\alpha$ changes with $\\phi$ or with multipole truncation, the derivation fails.","supporting_citations":[{"cited_title":"Correlated diffusion of membrane proteins and their effect on membrane viscosity","cited_arxiv_id":null,"evidence_quote":"Supplies the approximate Faxén relations and the force-dipole derivation of effective viscosity that this paper extends numerically and confirms in the small-particle limit."},{"cited_title":"Henle and A.J","cited_arxiv_id":null,"evidence_quote":"Gives the alternative analytical effective-viscosity prediction that the numerics reproduce after a 2/3 rescaling across the full range of $a/L_{\\rm sd}$."},{"cited_title":"Diffusion of complex objects embedded in free and supported lipid bilayer membranes: role of shape anisotropy and leaflet structure","cited_arxiv_id":null,"evidence_quote":"Introduced the interfacial regularized Stokeslet method for drag and diffusion of arbitrary membrane-embedded shapes, which this paper extends to advection and viscosity."},{"cited_title":"Calculating hydrodynamic interactions for membrane-embedded objects","cited_arxiv_id":null,"evidence_quote":"Supplies the regularized-Stokeslet Green's function formulas and implementation details used in the discretization."},{"cited_title":"The method of regularized Stokeslets in three dimensions: Analysis, validation, and application to helical swimming","cited_arxiv_id":null,"evidence_quote":"Establishes the method of regularized Stokeslets in three dimensions, the representation of solid bodies as constrained fluid regions on which the scheme relies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the grand resistance matrix formalism used in Appendix B to compute the orientationally averaged stresslet analytically."},{"cited_title":"Fluctuating hydrodynamics of multicomponent membranes with embedded proteins","cited_arxiv_id":null,"evidence_quote":"Supports the expectation that the linear effective-viscosity correction remains reasonable up to area fractions around 0.1, which underlies the dilute-limit interpretation."}],"review_version":1}