{"id":"65b14c77-8eb0-4447-9f58-64fa2b7ed776","arxiv_id":"2505.12033","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new Monte Carlo method measures non-local bending modulus and spontaneous curvature of asymmetric lipid bilayers, and near-gel simulations reveal a density-curvature buckling instability.","lead":"This paper uses Monte Carlo simulations of a coarse-grained lipid bilayer with unequal leaflet lipid numbers to measure membrane elastic constants, including the non-local bending modulus that ordinary periodic-boundary simulations cannot access. It shows the standard quadratic elasticity theory works for fluid asymmetric membranes but breaks down near the gel transition, where a compressed leaflet phase-separates and drives curvature instabilities.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The non-local bending measurement assumes the two arcs are independent open replicas; curvature-dependent seam energy or a constraint-induced shift in local κ would directly bias κnl, which has no independent benchmark.","rationale":"I read the paper as proposing a constrained Monte Carlo method to measure κ, κnl, C*, and the area modulus for asymmetric bilayers, with validation for κ on a symmetric CD membrane (κ=7±0.5 kBT vs 8±1 kBT in ref. [37]) and a claim of theory agreement far from the gel transition. The near-gel section is explicitly operational, using 'effective' free energies. The most load-bearing step is the identification of κnl: the no-exchange setting changes the ensemble, and Eq. (15) assigns the measured curvature stiffness to κ+κnl. Everything downstream—the fluid-asymmetry agreement, the statement that κnl≈25 kBT, and the interpretation of the near-gel stiffening—depends on that assignment. The paper provides no independent check of κnl, no error bars on the δ-dependent plots, and no code or data, so the reader's CONDITIONAL verdict is appropriate. My concern sharpens the reader's weakest assumption: it is not merely that interface energy could exist; rather, the no-exchange constraint is itself what generates the area-difference elasticity, so if the seam or the constraint changes the quadratic coefficient, the measured κnl is an ensemble artifact rather than the vesicle property defined in Eq. (8). The four-arc test isolates seam contributions by changing the number of seams while keeping the bulk material and arc geometry fixed. I therefore leave the verdict UNCHANGED: CONDITIONAL remains the right call, with the condition being a quantitative demonstration that Eq. (15) is uncontaminated by seam and constraint effects.","tokens_in":943,"tokens_out":1234,"duration_ms":147084,"concrete_test":"Run the no-exchange protocol with four arcs instead of two (two positive and two negative segments in the periodic cell, same total area, same N±, same θ|c| move set). If the fitted κtot changes by more than the reported ±0.3 kBT, the two seams in the original cell are not inert and Eq. (15) does not isolate the intrinsic κ+κnl.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The method's headline deliverable is κnl, obtained by subtracting κ (measured in the exchange setting) from κtot (measured in the no-exchange setting). Equation (15) presupposes that in the no-exchange setting the curvature free energy of one arc is (A0/2)(κ+κnl)(c−C*)². That presupposition requires the two arcs to behave as independent closed-vesicle-like replicas, with the two seams contributing only an irrelevant constant. This is not demonstrated. The arcs are open patches joined periodically; the global average curvature is zero, so the non-local term in Eq. (8) is constant for the whole box. The extra stiffness in the no-exchange setting therefore arises entirely from the constraint that leaflet numbers are fixed separately in each arc. The coefficient of that constrained quadratic term equals κ+κnl only if (i) the seam energy is curvature-independent and negligible, and (ii) the local κ under the no-exchange constraint equals the κ=7±0.5 kBT measured with exchange. The manuscript provides no check of either condition; it even concedes in Sec. V that inter-leaflet and inter-segment coupling is non-negligible when the membrane is heterogeneous. Since κnl≈25.4 kBT is the paper's novel quantity and has no independent benchmark, a curvature-dependent seam contribution or a constraint-induced shift in local κ would directly change the central claim. This is the same load-bearing assumption identified by the reader, and it is the principal soft spot in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents Monte Carlo simulations of an ultra coarse-grained (Cooke-Deserno) lipid bilayer with unequal leaflet populations, using a two-arc geometry with opposite curvatures under periodic boundary conditions. A key methodological innovation is the use of rejection rules that prevent flip-flop and optionally prevent lipid exchange between the two arcs, allowing the system to sample constrained ensembles. From the Gaussian free energies obtained by fitting the distributions of area and curvature, the author extracts the area elasticity K, the local bending rigidity κ, and, in the no-exchange setting, the total bending rigidity κ_tot = κ + κ_nl and the spontaneous curvature C*. For fluid membranes (ε = 1.05 kBT) the symmetric validation gives κ = 7.0 ± 0.5 kBT, consistent with the earlier Fourier-spectrum value of 8 ± 1 kBT, and the asymmetric results agree with the quadratic monolayer-additive theory. Near the gel transition (ε = 1.3 kBT), increasing asymmetry triggers phase separation in the compressed leaflet, density variations in the dilated leaflet, and a buckling instability, with an abrupt rise in the effective bending rigidity. The paper claims that the method can measure the non-local bending modulus and spontaneous curvature, which are not accessible in standard simulations with periodic boundary conditions.","tokens_in":17006,"tokens_out":8053,"duration_ms":84617,"significance":"If the proposed method is valid, it would fill a real gap: the non-local bending modulus κ_nl is difficult to measure in simulations of periodic membranes, and a practical route to it would be valuable for studying vesicle processes such as budding and fusion. The paper also offers a plausible mechanism linking density asymmetry, gel-phase separation, and curvature instability, which is relevant to raft-forming membranes. Strengths of the work include the validation of the symmetric bending rigidity against published Fourier-spectrum data, the clean Gaussian fits used to extract elastic parameters, and the operational definition of an effective bending rigidity in the two-phase regime. However, the headline new observable κ_nl is not checked against any independent benchmark, and its extraction rests on assumptions about the independence of the two arcs and about the invariance of the local bending rigidity under the no-exchange constraint. The significance of the paper therefore hinges on whether these assumptions can be validated.","major_comments":[{"comment":"The measurement of κ_nl rests on treating the two cylindrical arcs as independent closed-vesicle-like replicas. In the periodic simulation box the two arcs have equal area and opposite curvature, so the global average curvature is identically zero and the non-local term in Eq. (8) is constant for the whole system. The additional curvature stiffness observed in the no-exchange setting must therefore come from the constraint that leaflet counts are fixed separately in each arc. Interpreting that stiffness as (κ + κ_nl)/2 per arc, as in Eq. (15), requires two conditions: (i) the seam energies at the two interfaces are negligible and independent of curvature, and (ii) the local bending rigidity under the no-exchange constraint equals the value measured with exchange. Neither condition is demonstrated. Since κ_nl ≈ 25.4 kBT is the paper's central new result and has no independent benchmark, this is a load-bearing gap. I recommend adding tests such as varying the box length or seam geometry, computing the seam energy as a function of c, or measuring κ_nl through an independent route, for example a closed vesicle with imposed area difference.","section":"Section III, Eq. (15) and Fig. 3"},{"comment":"The value κ_nl = κ_tot − κ is obtained by subtracting the local bending rigidity measured in the exchange setting (Fig. 2) from the total bending rigidity measured in the no-exchange setting (Fig. 3). The paper verifies that the area parameters A0 and K are similar in the two settings (Figs. 4 and 5), but it does not verify that the local bending rigidity is unchanged by the no-exchange constraint. A constraint-induced shift in the local κ would directly bias κ_nl by the same amount. This possibility should be checked, for instance by comparing results for systems of different sizes or by designing an independent way to isolate the local curvature response in the no-exchange ensemble.","section":"Section III, Figs. 2–3"},{"comment":"The manuscript explicitly concedes that inter-leaflet and inter-segment coupling is non-negligible when the membrane is heterogeneous, stating that this coupling has been 'completely ignored thus far, both in sections II and V.' This admission is made for the phase-separated regime near the gel transition, but it raises a quantitative question for the fluid regime as well: how large is the same coupling at ε = 1.05 kBT, where Eqs. (10) and (15) are used to extract κ and κ_nl? A control calculation or an estimate of the neglected coupling in the fluid regime would substantially strengthen the central claim.","section":"Section V, after Fig. 10"}],"minor_comments":[{"comment":"The statement that area and curvature are independent quadratic degrees of freedom is only approximate because the curvature term in Eq. (10) contains a factor A. The text should state explicitly that A is replaced by A0 and should estimate the size of the resulting A–c coupling.","section":"Section III, Eq. (10)"},{"comment":"The notation for the densities in Fig. 10 appears to list ρ^2_1 twice; it should presumably be ρ^1_1, ρ^2_1, ρ^1_2, and ρ^2_2. Please correct this.","section":"Fig. 10 and accompanying text"},{"comment":"The comparison of the area elasticity modulus K between the exchange and no-exchange settings is described as 'fairly good agreement,' but no quantitative statement of the difference relative to the error bars is given. Please add this information.","section":"Section IV, Figs. 4–5"},{"comment":"The fit shown in Fig. 8(b) for δ = 0.10 is made to one minimum of the double-well free energy. Please state the fitting range and the criterion used to exclude the unstable branch between the two minima.","section":"Section V, Eq. (22)"},{"comment":"There are several typographical errors, for example 'membran es' in the abstract and title area, and 'they can are joined smoothly' in Section III. These should be corrected before publication.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague —\n\nThis paper is worth reading. The MC method is new and clever: by rejecting flip-flops and, in the no-exchange setting, also rejecting lipid transfer between the two arcs, Farago can keep leaflet asymmetry in the fast-flip-flop Cooke–Deserno model. That alone is a solid technical advance. The symmetric validation is good: κ = 7.0 ± 0.5 kBT from the curvature distribution matches the published Fourier-spectrum value of 8 ± 1 kBT. The fluid asymmetric results (weak δ-dependence of K, κ, a0; linear C*(δ)) agree cleanly with the quadratic theory. The near-gel part is also interesting: the compressed leaflet appears to phase-separate and the effective κ jumps, which matches experimental reports of stiffening in asymmetric vesicles.\n\nThe concern is the headline κnl ≈ 25.4 kBT. The extraction assumes the two arcs are independent closed-vesicle replicas, with seam energies contributing only a constant. That assumption is not demonstrated. In the no-exchange setting the global average curvature is zero, so the non-local term in Eq. (8) is constant for the whole box; the observed extra stiffness comes entirely from the per-segment leaflet constraints. Calling it κ+κnl requires that the seams are curvature-independent and that the local κ is unaffected by the no-exchange constraint. Neither is checked, and κnl has no independent benchmark. This is the principal soft spot.\n\nOther issues are minor: the δ-dependent plots in Secs. IV–V lack error bars, and no code/data are provided, so the numbers cannot be independently checked. The near-gel identification is indirect (diffusion drop, density discontinuity), but it is framed honestly, including the admission that inter-leaflet coupling matters there.\n\nWho is this for? Researchers who need a practical simulation route to spontaneous curvature and non-local bending in asymmetric, periodically bounded membranes. The method is promising enough that a serious referee should engage. I would send it out, with a request to benchmark κnl — for instance by varying the spanning angle or box dimensions to show κtot is geometry-independent, or by comparing against a vesicle simulation. I would not cite the κnl value itself until that is done.","headline":"A promising MC method for measuring spontaneous curvature and, in principle, the non-local bending modulus of asymmetric membranes, but the headline κnl rests on an unverified independence assumption and needs a benchmark.","tokens_in":17472,"tokens_out":4871,"would_cite":false,"duration_ms":48561,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A Monte Carlo setup with two oppositely curved membrane arcs measures the non-local bending modulus and shows that leaflet asymmetry triggers buckling near the gel transition.","keywords":["asymmetric lipid bilayer","Monte Carlo simulation","bending rigidity","non-local bending modulus","spontaneous curvature","membrane elasticity","gel transition","curvature instability"],"falsifier":"Measure $\\kappa$ and $\\kappa_{\\rm nl}$ with the same two-arc protocol at several spanning angles; if the extracted moduli drift with arc length while the membrane stays fluid, the independent-segment assumption fails.","tokens_in":16478,"feed_emoji":"🧪","tokens_out":8963,"duration_ms":81791,"temperature":0.7,"pith_summary":"This paper develops a Monte Carlo protocol for measuring the elastic parameters of lipid bilayers whose two leaflets hold different numbers of lipids. The membrane is arranged as two oppositely curved cylindrical arcs in a periodic box, and rejection rules freeze the leaflet asymmetry; in a second mode, exchange between the arcs is blocked so each arc relaxes independently. From the probability distributions of area and curvature, the simulations extract the area compression modulus $K$, the local bending rigidity $\\kappa$, the spontaneous curvature, and the non-local bending modulus $\\kappa_{\\rm nl}$ that standard fluctuation methods cannot measure. For fluid membranes far from the gel transition the measured values agree with the quadratic monolayer-additive theory, while near the transition the method reveals an asymmetry-driven gel phase and a buckling instability with an abrupt rise in effective stiffness.","feed_headline":"Two curved arcs reveal a membrane's hidden bending stiffness","feed_subtitle":"Monte Carlo rules freeze leaflet asymmetry and expose the non-local bending modulus, plus buckling near the gel transition.","key_machinery":"The central object is the two-arc geometry: a periodically repeated box containing two cylindrical arcs of equal spanning angle and opposite curvatures $\\pm c$, with a rejection rule that confines each lipid to its own leaflet and prevents escape. In the free-exchange version, lipids can diffuse between the arcs so the average curvature stays zero and the non-local term is constant; in the blocked-exchange version, each arc relaxes independently and the curvature free energy takes the form $(\\kappa+\\kappa_{\\rm nl})(c-C^*)^2/2$. The logarithms of the sampled area and curvature distributions are fitted to these quadratic forms, and the fitted coefficients give the elastic moduli.","core_discovery":"On its own terms, the paper shows that Monte Carlo sampling restricted to a partially constrained configuration space—two cylindrical arcs of opposite curvature joined smoothly under periodic boundary conditions—can act as a mechanical balance for an asymmetric bilayer. Fitting the logarithms of the measured area and curvature distributions to quadratic free energies yields $K$, $\\kappa$, the spontaneous curvature, and the non-local bending modulus $\\kappa_{\\rm nl}$, which is ordinarily invisible at fixed average curvature. In fluid membranes far from the liquid–gel transition, the extracted parameters obey the quadratic monolayer-additive theory: the bilayer bending modulus is the sum of the leaflet values, the spontaneous curvature grows linearly with the asymmetry parameter, and $\\kappa_{\\rm nl}$ is large compared with $\\kappa$. Near the gel density the quadratic theory fails: the compressed leaflet phase-separates, the two arcs buckle to opposite curvatures, and the effective bending rigidity rises from about 12.5 to 18 $k_{\\rm B}T$ across the transition.","pith_inferences":["The interface between the two arcs is never separately accounted for; a direct check would be to repeat the protocol with different spanning angles and see whether the extracted $\\kappa$ and $\\kappa_{\\rm nl}$ stay constant.","Because $\\kappa_{\\rm nl}$ turns out to be several times larger than $\\kappa$, area-difference elasticity should dominate whenever a vesicle's overall curvature changes substantially, such as in budding or tether pulling; the same method could parameterize those continuum models.","The curvature–density coupling mechanism may operate continuously in liquid-ordered raft domains surrounded by a disordered liquid matrix, not only at a first-order gel transition; this could be tested with coexisting-liquid simulations or composition gradients.","Near the transition the fitted free energies are operational rather than true thermodynamic potentials, so effective bending rigidities extracted from one minimum should be interpreted with care and compared with single-segment measurements."],"forward_implications":["A flat, periodically repeated bilayer can be used to measure the non-local bending modulus $\\kappa_{\\rm nl}$, which thermal-fluctuation spectra cannot determine.","For fluid membranes away from the gel transition, the full set of elastic parameters—$K$, $\\kappa$, spontaneous curvature, and $\\kappa_{\\rm nl}$—follows the quadratic monolayer-additive theory, with $K$ and $\\kappa$ depending only weakly on leaflet asymmetry and the spontaneous curvature growing linearly with $\\delta$.","Near the gel transition, increasing asymmetry makes the compressed leaflet phase-separate and drives the two arcs to buckle to opposite curvatures, with the effective bending rigidity jumping from about 12.5 to 18 $k_{\\rm B}T$.","Phase separation in the compressed leaflet can induce density differences between the oppositely curved segments of the dilated liquid leaflet, demonstrating cross-leaflet mechanical coupling.","The density–curvature buckling mechanism offers a physical route by which raft-like domains could generate local curvature changes in cellular membranes."],"supporting_citations":[{"why":"Supplies the solvent-free ultra-coarse-grained lipid model whose three-bead lipids are simulated with the new Monte Carlo moves and rejection rules.","marker":"[23]"},{"why":"Provides the benchmark bending rigidity of the same model from Fourier fluctuation analysis, used to validate the two-arc method.","marker":"[37]"},{"why":"Derives the monolayer-additive expressions for the bilayer area, bending, and non-local moduli that the simulations are compared against.","marker":"[29]"},{"why":"Introduces the area-difference-elasticity term that the non-local bending modulus $\\kappa_{\\rm nl}$ is identified with.","marker":"[31]"},{"why":"Earlier simulation study of asymmetric membranes reporting spontaneous curvature and bending modulus changes, providing the comparison for the near-gel stiffening.","marker":"[26]"},{"why":"Experimental tether-formation measurement that defines local and non-local curvature elasticity and shows why the non-local modulus matters for geometry changes.","marker":"[35]"},{"why":"Experiment observing increased bending rigidity in asymmetric vesicles, matched by the simulation's abrupt rise in effective $\\kappa$ near the transition.","marker":"[41]"}],"fun_headline_variants":["Asymmetric membrane arcs expose non-local bending modulus","New Monte Carlo method measures asymmetric bilayer elasticity","Two curved arcs probe membrane bending and spontaneous curvature","Monte Carlo arcs reveal asymmetric membrane bending stiffness","Arcs measure hidden bending modulus in asymmetric membranes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The two curved membrane segments are treated as independent quadratic springs whose area and curvature energies add, with no significant energy stored at the joints where the arcs meet.","fun_headline_variants_meta":{"raw":{"variants":["Asymmetric membrane arcs expose non-local bending modulus","New Monte Carlo method measures asymmetric bilayer elasticity","Two curved arcs probe membrane bending and spontaneous curvature","Monte Carlo arcs reveal asymmetric membrane bending stiffness","Arcs measure hidden bending modulus in asymmetric membranes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001371,"raw_usage":{"total_tokens":5590,"prompt_tokens":1011,"completion_tokens":4579,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":627,"completion_tokens_details":{"reasoning_tokens":4509}},"tokens_in":627,"tokens_out":4579,"duration_ms":28943,"temperature":1.0,"reasoning_tokens":4509,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:43:34.138375+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure $\\kappa$ and $\\kappa_{\\rm nl}$ with the same two-arc protocol at several spanning angles; if the extracted moduli drift with arc length while the membrane stays fluid, the independent-segment assumption fails.","supporting_citations":[{"cited_title":"Blumer, S","cited_arxiv_id":null,"evidence_quote":"Supplies the solvent-free ultra-coarse-grained lipid model whose three-bead lipids are simulated with the new Monte Carlo moves and rejection rules."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the benchmark bending rigidity of the same model from Fourier fluctuation analysis, used to validate the two-arc method."},{"cited_title":"Takiue, Heterogeneity and deformation behavior of l ipid vesicles, Current Opinion in Colloid Interface Scienc e 62, 101646 (2022)","cited_arxiv_id":null,"evidence_quote":"Derives the monolayer-additive expressions for the bilayer area, bending, and non-local moduli that the simulations are compared against."},{"cited_title":"Svetina and B","cited_arxiv_id":null,"evidence_quote":"Introduces the area-difference-elasticity term that the non-local bending modulus $\\kappa_{\\rm nl}$ is identified with."},{"cited_title":"Foley and M","cited_arxiv_id":null,"evidence_quote":"Earlier simulation study of asymmetric membranes reporting spontaneous curvature and bending modulus changes, providing the comparison for the near-gel stiffening."},{"cited_title":"Deserno, Fluid lipid membranes: From differential ge ometry to curvature stresses, Chemistry and Physics of Lipi ds 185, 11 (2015)","cited_arxiv_id":null,"evidence_quote":"Experimental tether-formation measurement that defines local and non-local curvature elasticity and shows why the non-local modulus matters for geometry changes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experiment observing increased bending rigidity in asymmetric vesicles, matched by the simulation's abrupt rise in effective $\\kappa$ near the transition."}],"review_version":1}