{"id":"45eb53c2-64e7-496f-bc72-a925b5729f93","arxiv_id":"2608.11039","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The neck energy of a partially engulfed particle scales as epsilon^2 ln(1/epsilon) with a closed-form coefficient 2(H0 - c0)^2 sin^4(theta_c), leading to tension-independent binding and hysteresis.","lead":"This paper derives a closed-form expression for the elastic energy of the narrow membrane neck that forms when a tiny particle is partially engulfed by a lipid membrane. The result gives quantitative predictions for when wrapping begins, when it becomes complete, and why the transition can be hysteretic.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (52)'s log coefficient is inferred from an area-counting shortcut rather than derived from the Jacobi solution; the c0≠0 numerical check is off by ~20%, so the neck-energy law is not yet established.","rationale":"The reader's conditional verdict is appropriate, and my stress-test does not move it to accept or reject. The paper has real independent support: the binding threshold η²W = 2, being an exact statement about the second derivative at θc = 0, is robust; the tension-driven part of the neck energy is confirmed to 0.6% against full nonlinear solves (Fig. 3); and the deflection charge q = sin²θc matches to 0.4%. The soft spot is precisely the coefficient Γ of the log term, which is the quantity that controls the hysteresis width and the general closed-vesicle prediction. My concern sharpens the reader's weakest-assumption point in two ways. First, the path from (48) to (52) is a shortcut: the log coefficient is obtained by multiplying the leading-catenoid excess area by the ambient Helfrich density, rather than by evaluating the second-order layer energy on the actual solution of the Jacobi problem (42). Second, the numerical test of the spontaneous-curvature part shows a 20% slope deficit, and the paper's finite-domain explanation is not backed by a convergence study. The cross-term issue for \\tilde{H}0 ≠ 0 is a further unexamined risk: the area difference between the actual neck and the ambient curved membrane contains a term of order ε^{-1} that would threaten the claimed order of the expansion unless the mean-curvature terms cancel it. None of this demonstrates the result is wrong; it demonstrates that the central law is not yet established beyond reasonable doubt. A direct derivation from (42) into (48), or a numerical convergence study at smaller tension, would settle it. Hence the verdict remains conditional: the paper should be revised to supply that derivation or verification before the coefficient in (52) is used as a quantitative prediction.","tokens_in":21020,"tokens_out":14225,"duration_ms":125609,"concrete_test":"Substitute the explicit normal displacement w(u) of Eq. (42), with b1, b2 fixed by w(uc) = 0 and w'(uc) = 0, into the layer energy integral (48) written in the conformal variable u, and asymptotically extract the coefficient of ε^2 ln(1/ε). If it differs from πΓ sin^4θc with Γ = 2(\\tilde{H}0 - σ0)^2, then Eq. (52) is not the correct log coefficient. As a cross-check, repeat the Fig. 4b second-difference calculation at Σ = 10^{-5} and 10^{-6} with the outer truncation at ≥ 100λ and verify whether the fitted slope converges to 2π sin^4θc (3.53 at θc = 60°) or remains at the observed ≈ 2.8.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the neck energy (52): E_neck = πk Γ sin^4θc ε^2 ln(1/ε) with Γ = 2(\\tilde{H}0 - σ0)^2. The derivation in Section 4 moves from the exact layer integral (48), which involves the first-order mean-curvature deviation (m1 - 2σ0)^2, to the final coefficient by asserting that the leading logarithmic energy is 'the excess area (51) multiplied by the Helfrich energy density of the ambient membrane.' No substitution of the explicit solution w(u) from Eq. (42) into (48) is shown; the coefficient is obtained by an overlap-region area-counting argument that is plausible but not a derivation. For a curved ambient (\\tilde{H}0 ≠ 0), the area excess of the actual neck profile relative to the ambient contains a cross term of order ε^{-1} in physical units (from the product of the catenoid tail q/r and the ambient slope H0 r), which would enter at lower order than ε^2 ln unless canceled by the mean-curvature variation; the paper does not analyze this cancellation. The numerical support is also incomplete: the tension part (Fig. 3) matches to 0.6%, but the spontaneous-curvature part (Fig. 4b) shows a fitted slope 2.8 versus the predicted 2πq^2 = 3.53, a 20% deficit attributed to a 'finite-domain effect' without a convergence study. Since the hysteresis width (77) is proportional to Γ0 = \\barΣ + 2σ0^2, an error in this coefficient propagates directly to a central quantitative prediction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":21372,"tokens_out":8399,"duration_ms":82684,"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":[{"comment":"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.","section":"Section 4, Eqs. (48)-(52)"},{"comment":"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.","section":"Appendix A, Fig. 4b"},{"comment":"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.","section":"Section 5.2 and Eq. (60)"}],"minor_comments":[{"comment":"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.","section":"Section 3.2, Eq. (39)"},{"comment":"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.","section":"Section 4, Eq. (50)"},{"comment":"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.","section":"Section 5.1, Eq. (56)"},{"comment":"The fitted slope is quoted as 2.8 in the text and 2.83 in the figure caption; please harmonise these numbers.","section":"Appendix A, Fig. 4b"},{"comment":"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.","section":"Appendix A, Eq. (79)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and the asymptotic programme is promising, but the central neck-energy law is not yet fully supported: the Section 4 derivation of the log coefficient is heuristic, and the spontaneous-curvature numerical check in Fig. 4b shows a 20% deficit that is not explained by a convergence study. Both issues are fixable within the manuscript's scope, so I do not see grounds for rejection, but they need to be addressed before the paper can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's central claim — E_neck = πkΓ sin^4θ_c ε^2 ln(1/ε) with Γ = 2(Ĥ_0 − σ_0)^2 — is genuinely new, and the tension-driven part is solidly confirmed by full nonlinear solves. But the spontaneous-curvature part of the coefficient is derived by an area-counting shortcut rather than by evaluating the layer integral, and its one numerical test is 20% off. That is a real gap, and the referee should push on it.\n\nWhat is actually new and good: previous treatments stop at the catenoid and either zero the neck energy or keep only an O(ε) line tension. This paper correctly identifies that the cost lives in the O(ε) correction to the catenoid and appears at ε^2 ln(1/ε). The matched-asymptotic structure is coherent: leading catenoid, inhomogeneous Jacobi equation for the first correction, log from the neck's 1/r tail. The deflection charge q = a sin^2θ_c is fixed by matching, not fitting, and the full BVP confirms it to 0.4%. The binding threshold η_W^2 = 2, independent of tension and spontaneous curvature, is exact at leading order, and the cap–neck cancellation in ∫2H (Eqs. 62–63) that defers the bilayer-asymmetry shift to complete wrapping is elegant and numerically verified. The hysteresis loop is a testable prediction. The numerical support for the tension part is genuinely strong: slope 1.564 versus π/2, a 0.6% discrepancy, across three wrapping angles, with no free parameters.\n\nSoft spots, in proportion. First, the coefficient in Sec. 4 is obtained by the area-counting argument — excess area of the leading catenoid tail multiplied by the ambient Helfrich density — not by a direct evaluation of the layer integral (48) with the explicit solution (42). The paper asserts that the integrand confirms the result term by term but shows no such evaluation. This matters most for a curved ambient (Ĥ_0 ≠ 0), where the cross term between the neck's 1/r slope and the ambient slope is not analyzed and could enter at a competitive order. The reservoir closure actually solved in the paper has Ĥ_0 = 0, so the main quantitative results are protected; the general formula (52) is not yet fully backed. Second, the only test of the spontaneous-curvature part (Fig. 4b) gives a fitted slope of 2.8 versus the predicted 3.53, a 20% deficit dismissed as a finite-domain effect without a convergence study. Since the hysteresis width (77) is proportional to this coefficient, that needs to be resolved before publication.\n\nThe circularity concern does not land: the matched results are parameter-free and independently verified by the nonlinear solves. This paper is for membrane-mechanics and nanoparticle-wrapping readers. It deserves a serious referee, with the request to fill in the direct layer-integral evaluation and the convergence study.\n\nRecommendation: send to peer review. Expect a revision that tightens Sec. 4 and the c_0 ≠ 0 numerics.","headline":"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.","tokens_in":21882,"tokens_out":13083,"would_cite":true,"duration_ms":108161,"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":"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$.","keywords":["Helfrich membrane","partial engulfment","elastic boundary layer","catenoid","Jacobi equation","neck energy","matched asymptotic expansions","wrapping hysteresis"],"falsifier":"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.","tokens_in":20798,"feed_emoji":"🫧","tokens_out":11052,"duration_ms":97063,"temperature":0.7,"pith_summary":"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.","feed_headline":"Nanoparticle wrapping pays a neck penalty leading-order theory misses","feed_subtitle":"The cost is set by curvature mismatch and appears only at ε² ln(1/ε), making the wrapping transition discontinuous.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"It supplies the matched-asymptotic boundary-layer treatment of inclusion-induced membrane deformations that the paper adapts to the wrapping neck.","marker":"[7]"},{"why":"It gives the classical elastic deformation analysis of a bound colloid and the catenoidal neck that the present expansion refines.","marker":"[5]"},{"why":"It supplies the tension-clamped neck energetics and the elementary line-tension picture that the paper extends to order $\\epsilon^2\\ln(1/\\epsilon)$.","marker":"[8]"},{"why":"It establishes the wrapping energy landscape and the threshold $\\eta^2=2$ against which the binding result is compared.","marker":"[4]"},{"why":"It provides the bilayer-asymmetry shift of the engulfment threshold that the paper recovers at complete wrapping but not at onset.","marker":"[3]"},{"why":"It sets out the matched asymptotic expansion method used throughout Sections 3 and 5.","marker":"[13]"},{"why":"It supplies the Jacobi operator and stability theory of minimal surfaces used to solve the first correction.","marker":"[14]"},{"why":"It motivates the effective-source description of the neck as a deflection charge seen by the outer membrane.","marker":"[15]"},{"why":"It provides the collocation solver used to verify the neck-energy law and deflection charge against the full nonlinear equations.","marker":"[24]"}],"fun_headline_variants":["Neck penalty for nanoparticle wrapping emerges only at ε² ln(1/ε)","Leading-order theory misses the neck penalty for wrapping","The hidden cost of partial engulfment: a logarithmic neck energy","Neck energy is always a penalty, never a gain, in engulfment","Why partial engulfment has a discontinuous wrapped-transition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Neck penalty for nanoparticle wrapping emerges only at ε² ln(1/ε)","Leading-order theory misses the neck penalty for wrapping","The hidden cost of partial engulfment: a logarithmic neck energy","Neck energy is always a penalty, never a gain, in engulfment","Why partial engulfment has a discontinuous wrapped-transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000918,"raw_usage":{"total_tokens":3947,"prompt_tokens":960,"completion_tokens":2987,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":2912}},"tokens_in":576,"tokens_out":2987,"duration_ms":21674,"temperature":1.0,"reasoning_tokens":2912,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:32:23.572119+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Biscari and G","cited_arxiv_id":null,"evidence_quote":"It supplies the matched-asymptotic boundary-layer treatment of inclusion-induced membrane deformations that the paper adapts to the wrapping neck."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the classical elastic deformation analysis of a bound colloid and the catenoidal neck that the present expansion refines."},{"cited_title":"Fessler, P","cited_arxiv_id":null,"evidence_quote":"It supplies the tension-clamped neck energetics and the elementary line-tension picture that the paper extends to order $\\epsilon^2\\ln(1/\\epsilon)$."},{"cited_title":"Deserno and T","cited_arxiv_id":null,"evidence_quote":"It establishes the wrapping energy landscape and the threshold $\\eta^2=2$ against which the binding result is compared."},{"cited_title":"Agudo-Canalejo and R","cited_arxiv_id":null,"evidence_quote":"It provides the bilayer-asymmetry shift of the engulfment threshold that the paper recovers at complete wrapping but not at onset."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It sets out the matched asymptotic expansion method used throughout Sections 3 and 5."},{"cited_title":"Osserman.A Survey of Minimal Surfaces","cited_arxiv_id":null,"evidence_quote":"It supplies the Jacobi operator and stability theory of minimal surfaces used to solve the first correction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It motivates the effective-source description of the neck as a deflection charge seen by the outer membrane."},{"cited_title":"Kierzenka and L","cited_arxiv_id":null,"evidence_quote":"It provides the collocation solver used to verify the neck-energy law and deflection charge against the full nonlinear equations."}],"review_version":1}