Pith. sign in

REVIEW 3 major objections 5 minor 24 references

Boundary-layer analysis of the partial engulfment of a small particle by a lipid membrane

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Partial engulfment of a small sphere by a Helfrich membrane has a neck energy that vanishes at leading order and first appears as $\pi k\Gamma\sin^4\theta_c\,\epsilon^2\ln(1/\epsilon)$, with $\Gamma=2(\tilde H_0-\sigma_0)^2$.

desk verdict A genuinely new epsilon^2 ln(1/epsilon) neck-energy law backed by strong numerics for the tension part, but the spontaneous-curvature coefficient needs a real derivation and a convergence study before publication. read the letter →

arxiv 2608.11039 v1 pith:PCIWL7AT submitted 2026-08-11 cond-mat.soft math-phmath.MP

Gaetano Napoli This is my paper · ORCID
classification cond-mat.softmath-phmath.MP
keywords HelfrichmembranepartialengulfmentelasticboundarylayercatenoidJacobiequationneckenergymatchedasymptoticexpansionswrappinghysteresis
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper proves that the neck joining a partially wrapped cap to the rest of a lipid membrane is an elastic boundary layer whose energy is invisible to the leading-order catenoid approximation. Because the catenoid is minimal, the first variation of the Helfrich energy vanishes; the cost is carried by the $O(\epsilon)$ correction to the catenoid and appears at order $\epsilon^2\ln(1/\epsilon)$ with a closed-form coefficient that is non-negative and vanishes only when the ambient membrane already carries its preferred curvature. Closing the outer field for a tension reservoir, the author derives a binding threshold independent of tension and spontaneous curvature, a complete-wrapping threshold shifted by bilayer asymmetry and tension, and a hysteresis loop whose width is set entirely by the neck term. The prediction is confirmed against full nonlinear shape equations at the percent level. If right, it says that any wrapping model that stops at the minimal surface omits the whole neck contribution and will misdescribe the envelopment transition.

What carries the argument

The argument runs on matched asymptotic expansions that resolve the neck as an elastic boundary layer. In stretched variables the leading inner surface is a catenoid with throat radius $\sin^2\theta_c$, and the $O(\epsilon)$ correction is a normal displacement $w$ governed by the inhomogeneous Jacobi equation of the catenoid, $\mathcal J[w]=2\sigma_0+2(\tilde H_0-\sigma_0)\tanh u$, whose two homogeneous modes are rigid translations and throat-radius changes. The logarithm enters through the catenoid's logarithmic excess area over the overlap region, and matching the inner slope $\sim a\sin^2\theta_c/r$ to the outer Bessel-$K_0$ field fixes the deflection charge $q=a\sin^2\theta_c$. The layer is then integrated out and replaced by a scalar self-energy at the pole, leaving an outer problem that is closed by a single number.

What would settle it

Solve the exact axisymmetric shape equations for a tension reservoir with $c_0=0$ at fixed wrapping angles and plot $E_{\rm full}/(\Sigma\sin^4\theta_c)$ against $\ln(1/\Sigma)$ over two decades of tension; the paper predicts a single straight line of slope $\pi/2$ for all $\theta_c$. Any systematic separation of the curves by wrapping angle, or any $O(1)$ shift of the binding threshold $\eta^2=2$ when tension or spontaneous curvature is varied at fixed $\epsilon\sqrt{\bar\Sigma}\ll1$, would falsify the central law.

Watch

Extended reading notes

Core claim

The central claim is that the neck energy of partial engulfment is $$E_{\rm neck}=\pi k\Gamma\$sin^{4}$\theta_c\,\$epsilon^{2}$\ln(1/\epsilon)+O(\$epsilon^{2}$)k,\qquad \Gamma=2(\tilde H_0-\sigma_0)^2,$$ where $\epsilon$ is the particle-to-membrane size ratio, $\theta_c$ is the wrapping angle, and $\tilde H_0$ and $\sigma_0$ are the dimensionless mean and spontaneous curvatures of the ambient membrane at the pole. The leading inner surface is a catenoid, whose vanishing mean curvature makes the leading Helfrich density vanish pointwise; the cost is therefore a second-order quantity generated by the first correction to the catenoid. The logarithmic factor is the geometric signature of the catenoid's slow decay, and the coefficient is a perfect square, so the neck is always a penalty, never a gain.

Load-bearing premise

Everything rests on the particle being so much smaller than the membrane that the neck is a cleanly separated thin layer, and on the leading minimal neck's excess area supplying the whole logarithmic energy; if the membrane's response length is comparable to the particle, the two regions blend and the formulas lose their accuracy.

Editorial extensions

If this is right

  • The free state loses stability at $\eta^2=2$, independent of tension and spontaneous curvature; the classical onset of wrapping is unchanged by the neck.
  • Complete wrapping is favoured once $\eta^2>2+4\epsilon\sigma_0+\bar\Sigma\epsilon^2+\bar\kappa$, so bilayer asymmetry and tension raise the full-wrapping threshold while a negative Gaussian modulus lowers it.
  • The envelopment transition is an imperfect saddle-node bifurcation: the hysteresis width is $\Delta\eta^2\simeq 1.54\,\Gamma_0\epsilon^2\ln(1/\epsilon)$, so hysteresis is absent at leading order and opens only through the neck term.
  • Any wrapping energy that keeps the catenoidal neck but drops its first correction omits the entire neck contribution, whose $\sin^4\theta_c$ dependence is structurally different from the cap terms.
  • The neck penalty is non-negative and vanishes only at $\tilde H_0=\sigma_0$, so a membrane already relaxed to its preferred curvature can wrap a small particle through a neck that costs nothing.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: the same log-enhanced, second-order neck cost should appear in any near-minimal bridge, including budding and fission necks and non-spherical inclusions, so the mechanism is not specific to spherical particles.
  • Beyond the paper: the reduction of the wrapped particle to a point with a scalar self-energy implies that several wrapped particles on one vesicle should interact through the outer Bessel-$K_0$ fields, with an amplitude controlled by the same curvature-mismatch factor.
  • Beyond the paper: on a membrane with anisotropic ambient curvature, the neck energy depends on the location of the contact line, so a wrapped particle should drift toward points where the two principal curvatures of the ambient membrane are equal.
  • Beyond the paper: for a closed vesicle at fixed area and volume, $\tilde H_0$ and the logarithmic cutoff will differ from the reservoir values, shifting the thresholds without changing the order, sign, or angular dependence of the neck term.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper analyses axisymmetric partial engulfment of a small rigid sphere by a Helfrich membrane in the limit epsilon = a/R_ves << 1. It treats the neck joining the bound cap to the free membrane as an elastic boundary layer, shows that the leading-order inner surface is a catenoid with vanishing Helfrich energy, and obtains the first non-vanishing neck contribution by solving the inhomogeneous Jacobi equation of the catenoid. The central formal result is the closed-form neck-energy law (52), E_neck = pi k Gamma sin^4(theta_c) epsilon^2 ln(1/epsilon) + O(epsilon^2) k with Gamma = 2(tilde(H0) - sigma0)^2. For a tension reservoir the outer field is closed by matching, giving a binding threshold eta_W^2 = 2 independent of tension and spontaneous curvature, a complete-wrapping threshold, and a hysteresis loop whose width is proportional to the neck coefficient. Appendix A reports numerical solution of the full nonlinear shape equations, confirming the tension part of the neck energy to 0.6% and the deflection charge q = a sin^2(theta_c), while the spontaneous-curvature part is only partially confirmed.

Significance. If the central result holds, Eq. (52) provides the first closed-form neck self-energy at order epsilon^2 ln(1/epsilon) that is invisible to leading-order catenoid treatments, and it yields concrete, falsifiable predictions: a binding threshold independent of tension and spontaneous curvature, a complete-wrapping threshold with the bilayer-asymmetry shift, and a hysteresis loop set by the boundary-layer coefficient. The derivation is parameter-free, the tension part is verified against the full nonlinear boundary-value problem to 0.6%, and the geometric cancellation underlying Eq. (63) is confirmed to 0.4%. These are substantial strengths. The main caveat is that the spontaneous-curvature part of the coefficient is not yet quantitatively established, and the derivation of the log coefficient in Section 4 contains a heuristic step that needs to be made rigorous.

major comments (3)
  1. [Section 4, Eqs. (48)-(52)] The passage from the exact layer integral (48) to the coefficient Gamma is not derived from the explicit solution. The text asserts that the leading logarithmic energy is the excess area (51) multiplied by the ambient Helfrich density 2k(H0-c0)^2, but it never substitutes the solution (42) and the logarithmic tail (45) into the integrand rho0 (m1 - 2 sigma0)^2. That integrand is quadratic in the first variation of the mean curvature, and it is not shown that cross terms between the constant source 2 sigma0 and the decaying particular solution either cancel or integrate to exactly the claimed coefficient. Since Eq. (52) is the central result, this gap must be closed by an explicit integration of (48), or by a rigorous derivation of the area-counting equivalence, rather than by the present plausibility argument.
  2. [Appendix A, Fig. 4b] The only numerical test of the spontaneous-curvature part of Gamma reports a fitted slope of 2.8 against the predicted 2 pi q^2 = 3.53, a deficit of about 20%, which is attributed to a 'finite-domain effect' without any convergence study or quantitative estimate. Because the hysteresis width (77) is proportional to Gamma_0, this discrepancy propagates directly into a central quantitative prediction. Please add a systematic study of the fitted slope as a function of the truncation radius (or of ln(lambda/a)), and show that the slope approaches 3.53 in the asymptotic regime, or provide an analytic expression for the finite-domain correction that reconciles the measured 2.8 with the prediction.
  3. [Section 5.2 and Eq. (60)] The paper adds the outer tension contribution (60), E_out = pi Sigma q^2 ln(lambda/a), to the bending neck energy (52) to obtain the reservoir coefficient Gamma_0 = barSigma + 2 sigma0^2. It would be helpful to show explicitly, at the level of the exact energy (78), how the expansion of the full functional separates into the bending part (52), the tension part (60), and a controlled remainder, so that the reader can verify that no cross term between the bending and tension sectors contributes at order epsilon^2 ln(1/epsilon). Without this decomposition, the additivity of the two coefficients is asserted rather than demonstrated.
minor comments (5)
  1. [Section 3.2, Eq. (39)] The matching condition m1 -> 2 tilde(H0) as u -> infinity is used to fix c1, but it is introduced only in the text; it should be stated explicitly as a matching condition before being used.
  2. [Section 4, Eq. (50)] The sign convention for the graph height zeta0(r) and for the normal displacement w is not fixed; since the coefficient B in Eq. (45) changes sign with the orientation of w, a short statement of the convention would remove ambiguity.
  3. [Section 5.1, Eq. (56)] The exclusion of the B ln r mode is justified by the divergence of the tension energy, but only in a sketch; please make the divergence argument explicit, since the same logarithmic mode is later accepted in the neck region where it is cut off.
  4. [Appendix A, Fig. 4b] The fitted slope is quoted as 2.8 in the text and 2.83 in the figure caption; please harmonise these numbers.
  5. [Appendix A, Eq. (79)] The relation ln(lambda/a) = (1/2) ln(1/Sigma) is written in the chosen units a = k = 1; stating this explicitly before the equation would avoid confusion.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the neck-energy coefficient is derived from matched asymptotics and Jacobi-equation matching; the only self-citation is methodological and not load-bearing.

full rationale

The central claim (52) is not circular. The coefficient Γ = 2(\tilde H0 − σ0)^2 follows from the inhomogeneous Jacobi problem (31)–(42), with the source (39) fixed by the far-field matching condition m1 → 2\tilde H0 and the spontaneous-curvature term −2σ0; no parameter is fitted to data to produce this coefficient. The deflection charge q = a sin^2 θc is fixed by matching the 1/r slope of the catenoid tail (50) to the outer near field (57), which is a standard matched-asymptotic reduction, not a fit. The numerical verification in Appendix A solves the full nonlinear shape equations and independently confirms the tension part of the law to 0.6% (Fig. 3) and the deflection charge to 0.4%; the 20% deficit in the spontaneous-curvature part (Fig. 4b) is an acknowledged finite-domain effect and a validation gap, not a circularity. The only self-citation is [7] for the matched-asymptotic method, but the layer equations, catenoid solution, Jacobi operator, and matching are re-derived in the paper; the argument does not reduce to a prior result of the authors. The §4 shortcut—'excess area (51) multiplied by the Helfrich energy density of the ambient membrane'—is a plausible but under-shown derivation step, not an assumption of the answer: the excess area and ambient density are separately defined, and the paper states that expansion of the exact integrand (46) confirms the result term by term. A rigorous evaluation of (48) would strengthen the proof, but that is a correctness/rigor concern, not circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The analysis introduces no new physical entities. It depends on the standard Helfrich model, the axisymmetric and force-free assumptions, and the stated asymptotic scalings. The 'deflection charge' q is a matched-asymptotic bookkeeping quantity, not a new physical object.

assumptions (5)
  • domain assumption The membrane is described by the Helfrich energy (1) with constant bending rigidity k and spontaneous curvature c0.
    Stated in Section 2 and throughout; the entire analysis rests on this continuum model.
  • domain assumption The particle is held by adhesion alone, so the axial resultant f_z = 0 in (13).
    Invoked in Section 2 to make the axial multiplier negligible in the inner problem; a prescribed axial load would alter the boundary layer.
  • domain assumption The tension and pressure scale as Sigma a^2/k = O(epsilon^2) and P a^3/k = O(epsilon^3), and the ambient curvatures as c0 a = epsilon sigma0 and H0 a = epsilon tilde(H0).
    Section 3 defines the small-parameter regime; the asymptotic expansion is only valid in this moderate-tension, small-particle limit.
  • ad hoc to paper The inner expansion is regular: Psi = Psi0 + epsilon Psi1 + O(epsilon^2), with Lambda = O(epsilon).
    Equation (25) in Section 3. This regularity is essential for the Jacobi problem and is not derived from a more general principle.
  • domain assumption The outer membrane can be described by the linearized Monge-gauge equation k grad^4 h - Sigma grad^2 h = 0 (54).
    Section 5.1 linearizes the outer field assuming small slopes and moderate tension; the matching requires the overlap region a << r << lambda.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Boundary-layer analysis of the partial engulfment of a small particle by a lipid membrane." pith.science (2026). https://pith.science/paper/PCIWL7AT

@misc{pith2026260811039,
  author       = {Pith},
  title        = {Pith review of: Boundary-layer analysis of the partial engulfment of a small particle by a lipid membrane},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PCIWL7AT}},
  note         = {Machine review of arXiv:2608.11039}
}
abstract

We study the axisymmetric partial engulfment of a small rigid sphere by a fluid Helfrich membrane. When the size ratio $\eps$ between the particle and the membrane is small, the neck that joins the wrapped cap to the surrounding membrane is an elastic boundary layer, and we analyse it by matched asymptotic expansions. The leading inner surface is a catenoid, a minimal surface that stores no Helfrich energy, so that the energy of partial engulfment is carried by the first correction and appears only at order $\eps^{2}\ln(1/\eps)$. We obtain it by solving the inhomogeneous Jacobi equation of the catenoid, forced by the spontaneous curvature and by the ambient mean curvature that the neck must match, and we give its coefficient in closed form. The boundary layer can then be integrated out, and the neck replaced by a scalar self-energy carried at the pole, so that the outer field can be closed independently of the inner one. For a membrane coupled to a tension reservoir we derive the binding threshold, which turns out to be independent of both the tension and the spontaneous curvature, the complete-wrapping threshold, and the hysteresis of the envelopment transition. The neck self-energy law is confirmed against the full nonlinear shape equations.

Figures

Figures reproduced from arXiv: 2608.11039 by the authors.

Figure 1
Figure 1. Axisymmetric geometry of the model. A rigid sphere of radius [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Reservoir wrapping of a small particle ( [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Full nonlinear shape equations versus the effective law ( [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Spontaneous-curvature tests against the full nonlinear shape equations ( [PITH_FULL_IMAGE:figures/full_fig_p022_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [1]

    Lipowsky

    R. Lipowsky. The conformation of membranes.Nature, 349:475, 1991

    work page 1991

  2. [2]

    U. Seifert. Configurations of fluid membranes and vesicles.Adv. Phys., 46:13, 1997

    work page 1997

  3. [3]

    Agudo-Canalejo and R

    J. Agudo-Canalejo and R. Lipowsky. Critical particle sizes for the engulfment of nanopar- ticles by membranes and vesicles with bilayer asymmetry.ACS Nano, 9:3704, 2015. 20 2 3 4 5 6 7 8 9 ln(1/Σ) = 2ln(1/ϵ) + O(1) 6 8 10 12 14 16Efull / (Σ sin4θc) (a) neck energy: full nonlinear vs effective law θc = 45 ◦ (full BVP) θc = 60 ◦ (full BVP) θc = 75 ◦ (full B...

    work page 2015

  4. [4]

    Deserno and T

    M. Deserno and T. Bickel. Wrapping of a spherical colloid by a fluid membrane.Europhys. Lett., 62:767, 2003

    work page 2003

  5. [5]

    M. Deserno. Elastic deformation of a fluid membrane upon colloid binding.Phys. Rev. E, 69:031903, 2004

    work page 2004

  6. [6]

    Napoli and A

    G. Napoli and A. Goriely. Elastocytosis.J. Mech. Phys. Solids, 145:104133, 2020

    work page 2020

  7. [7]

    Biscari and G

    P. Biscari and G. Napoli. Inclusion-induced boundary layers in lipid vesicles.Biomech. Model. Mechanobiol., 6:297, 2007

    work page 2007

  8. [8]

    Fessler, P

    F. Fessler, P. Muller, and A. Stocco. Energetics and dynamics of membrane necks in particle wrapping.J. Colloid Interface Sci., 700:138524, 2025

    work page 2025

Show all 24 references
  1. [9]

    P. B. Canham. The minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell.J. Theor. Biol., 26:61, 1970

    1970

  2. [10]

    Helfrich

    W. Helfrich. Elastic properties of lipid bilayers: theory and possible experiments.Z. Naturforsch. C, 28:693, 1973

    1973

  3. [11]

    Seifert, K

    U. Seifert, K. Berndl, and R. Lipowsky. Shape transformations of vesicles: Phase diagram for spontaneous-curvature and bilayer-coupling models.Phys. Rev. A, 44:1182, 1991

    1991

  4. [12]

    Capovilla and J

    R. Capovilla and J. Guven. Stresses in lipid membranes.J. Phys. A: Math. Gen., 35:6233, 2002

    2002

  5. [13]

    E. J. Hinch.Perturbation Methods. Cambridge Texts in Applied Mathematics. Cambridge University Press, Cambridge, 1991

    1991

  6. [14]

    Osserman.A Survey of Minimal Surfaces

    R. Osserman.A Survey of Minimal Surfaces. Dover Publications, New York, 2nd edition, 1986

    1986

  7. [15]

    M. M. M¨ uller, M. Deserno, and J. Guven. Geometry of surface-mediated interactions. Europhys. Lett., 69:482, 2005. 21 30 45 60 75 90 wrapping angle θc (deg) 0.999 1.000 1.001 1.002 1.003 1.004 1.005 1.006 ∫ 2H dA / 2πasin2 θc 0.03% 0.10% 0.21% 0.33% 0.40% (a) total mean curva...

    2005

  8. [16]

    M. Deserno. Fluid lipid membranes: From differential geometry to curvature stresses. Chem. Phys. Lipids, 185:11, 2015

    2015

  9. [17]

    Abramowitz and I

    M. Abramowitz and I. A. Stegun.Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Dover Publications, New York, 1972

    1972

  10. [18]

    Agudo-Canalejo and R

    J. Agudo-Canalejo and R. Lipowsky. Uniform and janus-like nanoparticles in contact with vesicles: energy landscapes and curvature-induced forces.Soft Matter, 13:2155, 2017

    2017

  11. [19]

    J¨ ulicher and R

    F. J¨ ulicher and R. Lipowsky. Domain-induced budding of vesicles.Phys. Rev. Lett., 70:2964, 1993

    1993

  12. [20]

    B. J. Reynwar, G. Illya, V. A. Harmandaris, M. M. M¨ uller, K. Kremer, and M. Deserno. Aggregation and vesiculation of membrane proteins by curvature-mediated interactions. Nature, 447:461, 2007

    2007

  13. [21]

    Midya, T

    J. Midya, T. Auth, and G. Gompper. Membrane-mediated interactions between nonspher- ical elastic particles.ACS Nano, 17:1935, 2023

    1935

  14. [22]

    N. Li, N. Sharifi-Mood, F. Tu, D. Lee, R. Radhakrishnan, T. Baumgart, and K. J. Stebe. Curvature-driven migration of colloids on tense lipid bilayers.Langmuir, 33:600, 2017

    2017

  15. [23]

    Kralj-Igliˇ c, V

    V. Kralj-Igliˇ c, V. Heinrich, S. Svetina, and B. ˇZekˇ s. Free energy of closed membrane with anisotropic inclusions.Eur. Phys. J. B, 10:5, 1999

    1999

  16. [24]

    Kierzenka and L

    J. Kierzenka and L. F. Shampine. A bvp solver based on residual control and the MATLAB PSE.ACM Trans. Math. Software, 27:299, 2001. 22

    2001

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.