{"id":"5d3000b5-b2b2-4bbc-88d8-a7434c48668d","arxiv_id":"1909.00217","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Charge regulation can make adjacent membranes in a stack take opposite charges, producing attraction strong enough to dominate van der Waals forces and potentially driving thylakoid grana formation.","lead":"This paper models stacks of chargeable biological membranes and shows that the membranes can settle into alternating positive and negative charge states. The resulting electrostatic attraction can dominate over van der Waals forces, offering a possible explanation for how photosynthetic thylakoid membranes stay stacked.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central equilibrium claim rests on an unverified global minimization: Appendix §2's steepest-descent loop lacks initial-condition and local-minimum checks, so the alternating patterns and attractive pressures may be metastable artifacts.","rationale":"The reader's weakest assumption — that the steepest-descent iteration actually finds the global minimum — is exactly the flaw that matters. All subsequent conclusions, including the 'radical modification' of electrostatics and the dominance over vdW, are computed from σ* states obtained by that iteration. The analytic derivations in the Appendix (e.g., the one-surface integral reduction leading to Eq. (A15)) are self-consistent, and the capacitor-limit disjoining pressure matches the numerical slope, but these checks validate the minimization machinery only conditional on reaching the intended minimum. Since the paper supplies no convexity proof, no multi-start search, and no report of the initial condition used, the equilibrium interpretation of the patterns is unverified. This does not invalidate the model; it makes the central claim conditional on a feasible numerical check. The word 'quasiperiodic' in the abstract is a separate imprecision, but it is not load-bearing: the figures show finite alternating or near-alternating sequences rather than a rigorously quasiperiodic sequence, and the main claims do not depend on that word. Similarly, the uncertainty in biological parameters (Hamaker constant, ionic strength, site density) affects the quantitative comparison but not the qualitative existence of the symmetry-broken state. The appropriate disposition remains conditional acceptance, with the global-minimization test as the condition.","tokens_in":10745,"tokens_out":4792,"duration_ms":47198,"concrete_test":"For representative points in Figs. 3–4 (e.g., the three (α,χ) values used in Figs. 5–6, plus one point in each tail), rerun the outer minimization from (i) all 2^N sign-ordered endpoint profiles with σ*_j = ±1/2, and (ii) at least 1000 random initial vectors in [−1/2,1/2]^N for N = 3, 4, and 11, at several κL values. Keep the inner PB solve unchanged and compare the final Ω values and σ*_mca. If any lower minimum is found, or if different initial conditions converge to distinct σ* with energy differences above numerical tolerance, the reported diagrams are metastable and the equilibrium attraction claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline claim — that the equilibrium state of a stack can be a symmetry-broken alternating charge pattern with an attraction that dominates van der Waals forces — requires that the σ* profiles reported in Figs. 3–6 are global minima of the grand potential Eq. (1). The Appendix ('Minimization of the grand potential') describes only an outer steepest-descent loop, 'The iteration is done until convergence is reached,' with no stated initial condition, no multi-start protocol, and no test against other stationary profiles. The functional is not shown to be convex, and the sign-flip symmetry plus the N-dependent pattern changes (e.g., (−,+,−,−) → (−,+,−,+) → (+,+,−,+) in Fig. 4) indicate multiple candidate states. If a run converges to a metastable pattern, the dark-blue regions of the phase diagrams, the σ*_mca(L) curves in Fig. 6, and especially the binding-energy curves in Fig. 5 would no longer describe the equilibrium stack; the magnitude of the predicted attraction and its comparison with vdW would be unsupported. Thus the absence of a global-minimization check is load-bearing for the central claim, not a cosmetic omission.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a stack of N parallel, charge-regulating membranes in a symmetric monovalent electrolyte. The model is a Poisson-Boltzmann grand potential, Eq. (1), combined with a surface free energy containing adsorption parameters α and χ. By numerically minimizing this functional, the authors compute the mean charge asymmetry σ*_mca (Eq. (2)) in the α–χ plane for several N values, identify symmetry-broken alternating charge patterns, and calculate stack binding energies and disjoining pressures (Fig. 5). They argue that the resulting electrostatic attraction dominates van der Waals forces in thylakoid stacks. The Appendix contains an analytic reduction of the outside-field contributions and a capacitor-limit expression for the short-stack-width disjoining pressure, which is checked against the numerical slope in Fig. 5.","tokens_in":10987,"tokens_out":6616,"duration_ms":65568,"significance":"If the equilibrium-state claims hold, the paper provides a concrete charge-regulation mechanism for attractive forces in membrane stacks: symmetry-broken alternating charge patterns produce an attractive electrostatic pressure whose short-distance asymptote is parameter-free (σ²/(8ε0εr)). This is relevant to thylakoid grana and other biological stacks. The analytic Appendix is a clear strength, and the small-width disjoining-pressure prediction is a falsifiable, parameter-free result that matches the numerics. The phase diagrams for N>2 extend the authors' earlier two-surface work and are not fitted to any target data. The main risks are the lack of a documented global-minimization protocol and the unsupported quantitative comparison with van der Waals forces.","major_comments":[{"comment":"The central equilibrium claim requires that the σ* profiles in Figs. 3–6 are global minima of Eq. (1), but the minimization protocol is described only as an outer steepest-descent loop that runs 'until convergence is reached,' with no statement of initial conditions, no multi-start search, and no comparison with other stationary profiles. The functional is not shown to be convex, and the paper itself exhibits several distinct charging patterns for neighboring parameter values (e.g., the sequence (−,+,−,−) → (−,+,−,+) → (+,+,−,+) in Fig. 4), so the iteration could converge to a metastable branch. If so, the phase diagrams, the σ*_mca(L) curves in Fig. 6, and the binding energies in Fig. 5 would not describe the equilibrium stack, and the predicted dominance over van der Waals forces would be unsupported. Please add a multi-start or continuation protocol and identify the global minimum, or explicitly restrict the claims to locally stable states.","section":"Appendix 2 and Figs. 3–6"},{"comment":"The statement that the charge-regulation-induced attraction 'dominates quantitatively over the attractive van der Waals force' is asserted but not demonstrated. The only quantitative input is a Hamaker constant A≈4.8×10−20 J; no curve or expression shows how the computed disjoining pressure (the κL≲4 asymptote or the Fig. 5 curves) compares with the van der Waals pressure at the relevant separations. Because this comparison is the basis for the thylakoid conclusion, please include the actual calculation or a graph with both contributions.","section":"Sec. III, thylakoid paragraph; Conclusions"},{"comment":"The abstract advertises a 'quasiperiodic effective charge sequence,' but the term is never defined in the body and no example of a quasiperiodic sequence is given. The reported patterns are short alternating sequences such as (−,+,−,+) and (+,+,−,+); these are not shown to have any quasiperiodic structure. Please either define and demonstrate the quasiperiodic sequences or revise the abstract to describe the symmetry-broken alternating charge states actually documented in the text.","section":"Abstract"}],"minor_comments":[{"comment":"The text says that up to κL≲10 the disjoining pressure is attractive 'for all asymmetrically charged configurations,' but Fig. 5 shows only three parameter sets; please rephrase to refer to the computed cases or provide a systematic scan.","section":"Sec. III, Fig. 5 discussion"},{"comment":"The parameter θ=2 is introduced as the ratio of neutral sites to fixed charges, but the relation to the allowed interval σ*_j∈[−1/2,1/2] is not explained; a sentence defining the total site density would help.","section":"Sec. II"},{"comment":"The phrase 'non-electrostatic adsorption free energy penalty per ion' for the term −αη_j is confusing when α is negative in the relevant regions; consider calling it a binding or adsorption energy parameter.","section":"Sec. II"},{"comment":"The statement that the N=2 results are in 'perfect agreement' with Ref. [12] is followed by a description of quantitative changes (thinning of the tails); please specify whether the agreement is qualitative or quantitative and, if quantitative, show the comparison.","section":"Sec. III, N=2 paragraph"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the paper is a solid numerical extension of the authors' own two-surface charge-regulation work, and the analytic capacitor-limit result is a genuine strength. The two main risks are the undocumented global-minimization procedure and the unsupported 'quasiperiodic' claim in the abstract; both appear fixable in revision, so I do not see a need for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, the genuinely new physics is in the N>2 stack behavior: even/odd N differences in the charge patterns, the crossover across chi = -2alpha, and the binding-energy curves showing an electrostatic attraction that can dominate van der Waals. The N=2 case is mostly a check of the authors' earlier Soft Matter 2018 paper, with the outside reservoir added, and they say so plainly. Second, the central caveat is real: the paper claims equilibrium states but never demonstrates that the steepest-descent iteration finds global minima of the grand potential. No initial conditions, no multi-start, no check against other stationary profiles. Given the sign-flip symmetry and the pattern changes reported in Fig. 4, local minima are a live possibility. If the reported patterns are metastable, the phase diagrams and the predicted attraction could change. That concern is load-bearing, not cosmetic.\n\nWhat the paper does well: the analytic reduction in the Appendix is careful and the capacitor-limit disjoining pressure matches the numerical slopes. The work is not fitted to data; it is a parameter-free derivation in the mean-field charge-regulation model. The biological motivation (thylakoid grana) is reasonable, and the claim that simple DLVO balance is insufficient for membrane stacks is worth taking seriously. The paper also credits its own prior work appropriately instead of hiding it.\n\nSoft spots, in proportion. The global-minimization gap is the main one. A second, smaller issue: the abstract's phrase 'quasiperiodic effective charge sequence' is never demonstrated in the body; the patterns shown are periodic or simple alternations. 'Quasiperiodic' is an overreach. Also, the biological conclusions depend on parameters (alpha, chi, site spacing, neutral-site ratio) that are not pinned down by experiment, so the thylakoid relevance is suggestive rather than established.\n\nWho is this for? Soft-matter theorists and biophysicists working on charge regulation, membrane stacking, or non-DLVO forces. A referee can reasonably ask for the minimization to be made rigorous or at least explored, and for the 'quasiperiodic' wording to be fixed. The framework is coherent and the results are plausible; this is a conditional accept, not a reject.\n\nMy recommendation: send it to peer review. The paper deserves referee time, and the authors can address the concerns without changing the framework.","headline":"A solid extension of a two-surface charge-regulation model to stacks, with a real though fixable gap: the numerical minimization is not shown to find global minima.","tokens_in":11504,"tokens_out":940,"would_cite":true,"duration_ms":10786,"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":"Charge-regulating membrane stacks can settle into alternating charges that attract.","keywords":["charge regulation","membrane stack","Poisson-Boltzmann","symmetry breaking","surface charge","thylakoid","disjoining pressure","van der Waals"],"falsifier":"Run the same minimization from many random initial charge configurations (or an exhaustive grid) for the reported parameters and find a lower grand potential than the steepest-descent result; that would show the alternating state is not the equilibrium. Alternatively, measure the disjoining pressure of a stack of pH-responsive membranes at low salt and small separation: the model predicts an attractive pressure near $-11.3$ bar in the alternating regime, so a measured repulsion there would rule out the claim.","tokens_in":10545,"feed_emoji":"⚡","tokens_out":10898,"duration_ms":98047,"temperature":0.7,"pith_summary":"The paper argues that a stack of charge-regulating membranes—surfaces whose charge is set by proton (de)adsorption in equilibrium with a bathing salt solution—can settle into an alternating positive/negative charge pattern. In that symmetry-broken state the electrostatic interaction between neighboring membranes is attractive, and for parameters representative of thylakoid stacks this attraction dominates the van der Waals attraction. The authors establish this by minimizing a grand-potential functional over all surface charge densities for stacks of $N$ membranes, and they map out when the alternating state occurs as a function of reaction strength, salt concentration, and membrane separation. If the claim holds, modeling grana formation as fixed-charge surfaces balanced by van der Waals forces is insufficient; the stack's charge state itself must be treated as an equilibrium variable.","feed_headline":"Membrane stacks switch to alternating charges that attract","feed_subtitle":"Charge regulation makes thylakoid-like stacks bind electrostatically, beating van der Waals at short separations.","key_machinery":"The load-bearing object is the grand potential per unit area $\\beta\\Omega[\\sigma^*]$, Eq. (1): a mean-field functional of the surface charge densities $\\sigma^*_j$ that combines the Poisson-Boltzmann field energy in the electrolyte with surface terms for proton adsorption, including the adsorption penalty $-\\alpha\\eta_j$, the in-plane ion-ion interaction $-\\chi\\eta_j^2/2$, and the entropy of occupied and empty sites. Equilibria are the charge profiles that minimize this functional. The mechanism that produces alternating stacks is a near-cancellation of the $\\alpha$ and $\\chi$ terms on the line $\\chi=-2\\alpha$: with those terms balanced, the remaining electrostatic attraction between oppositely charged neighboring surfaces becomes the dominant contribution, selecting patterns of alternating sign. The paper evaluates the outside-region contributions analytically using the known single-surface Poisson-Boltzmann solution and minimizes numerically by nested loops over the potential and a steepest-descent step over the charges.","core_discovery":"The central discovery is that the equilibrium of a stack of charge-regulated membranes is not generally a uniform or symmetric charge profile. For $N\\ge 2$ membranes there are regions of the parameter plane (roughly around $\\chi=-2\\alpha$ with sufficiently negative $\\alpha$) where the lowest grand potential belongs to a quasiperiodic sequence of strongly charged surfaces alternating in sign, with mean charge asymmetry $\\sigma^*_{mca}\\approx 1$. Odd-$N$ stacks in this state carry a small net charge that flips sign as $\\chi$ crosses $-2\\alpha$, while even-$N$ stacks remain globally neutral but reorganize their pattern, e.g. for $N=4$ from $(-,+,-,-)$ to $(-,+,-,+)$ to $(+,+,-,+)$, to keep the in-plane ion-ion interaction favorable. In the alternating state the disjoining pressure (the force per unit area between the membranes) is attractive; at stack widths much smaller than the Debye length the attraction is $11.3$ bar per surface area, independent of $N$, and for thylakoid parameters it exceeds the van der Waals attraction even for a Hamaker constant one order of magnitude above standard estimates.","pith_inferences":["If the reported states are the true global minima, then the effective interaction between membranes in a stack is inherently collective: measuring the force between two isolated charge-regulated surfaces would not predict the stack's behavior, because the alternating pattern can only pay off when several surfaces share the charge imbalance.","A testable extension would be to expose a stack of pH-responsive surfaces to decreasing salt and watch for an abrupt crossover from repulsive to attractive disjoining pressure at the predicted $\\alpha$-$\\chi$ boundary; the crossover's position would probe the adsorption parameters directly.","The same mechanism may apply beyond biological membranes—for example to multilayered polyelectrolyte films or clay tactoids with dissociable surface groups, where alternating charge states could stabilize multilayer stacking without added multivalent ions.","Since the short-separation attraction is independent of $N$, this mechanism offers a reason why grana-like stacks prefer several layers: binding energy per area does not decrease as layers are added, while the collective charge pattern keeps every interface attractive."],"forward_implications":["At short stack widths $\\kappa L\\lesssim 4$ and high charge asymmetry, the electrostatic disjoining pressure approaches $-11.3$ bar and is independent of the number of membranes and of the reaction parameters, so adding layers does not dilute the per-area binding.","For thylakoid-relevant parameters, the charge-regulation attraction dominates the van der Waals attraction over the full salinity range considered, from a few mM up to 200 mM monovalent salt.","The predicted charge patterns are quasiperiodic rather than simply alternating for stacks of three or more membranes: odd-$N$ stacks have a net charge that flips at $\\chi=-2\\alpha$, while even-$N$ stacks stay neutral and instead reorganize their in-plane proton occupation.","Because the charge asymmetry persists up to large $\\kappa L$ for parameters on the $\\chi=-2\\alpha$ line, the electrostatic attraction remains active at salt concentrations and separations far beyond the range where fixed-charge Poisson-Boltzmann interactions would already be screened."],"supporting_citations":[{"why":"Supplies the grand-potential functional, the charge-regulation boundary conditions, and the two-surface symmetry-broken states that this paper extends to stacks.","marker":"[12]"},{"why":"Provides the thylakoid stack parameters (interdisk separation, ionic strength ranges) and the Hamaker constant used to compare the electrostatic attraction with van der Waals.","marker":"[20]"},{"why":"Gives the Poisson-Boltzmann equation and the analytic single-surface potential used to evaluate the electrolyte contributions outside the stack in the grand potential.","marker":"[29]"}],"fun_headline_variants":["Charge-regulated stacks alternate charges to attract","Membrane stacks flip charges in quasiperiodic pattern to attract","Alternating charges make membrane stacks attract beyond van der Waals","Charge regulation shapes membrane stacks into attractive alternators"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the numerical minimization always reaches the true global minimum of the grand potential; if the reported alternating states are metastable local minima rather than equilibrium states, the phase diagrams and the predicted attractions would change.","fun_headline_variants_meta":{"raw":{"variants":["Charge-regulated stacks alternate charges to attract","Membrane stacks flip charges in quasiperiodic pattern to attract","Alternating charges make membrane stacks attract beyond van der Waals","Charge regulation shapes membrane stacks into attractive alternators"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000195,"raw_usage":{"total_tokens":1330,"prompt_tokens":890,"completion_tokens":440,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":377}},"tokens_in":506,"tokens_out":440,"duration_ms":14778,"temperature":1.0,"reasoning_tokens":377,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:58:12.207417+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same minimization from many random initial charge configurations (or an exhaustive grid) for the reported parameters and find a lower grand potential than the steepest-descent result; that would show the alternating state is not the equilibrium. Alternatively, measure the disjoining pressure of a stack of pH-responsive membranes at low salt and small separation: the model predicts an attractive pressure near $-11.3$ bar in the alternating regime, so a measured repulsion there would rule out the claim.","supporting_citations":[{"cited_title":"Adži ´c and R","cited_arxiv_id":null,"evidence_quote":"Supplies the grand-potential functional, the charge-regulation boundary conditions, and the two-surface symmetry-broken states that this paper extends to stacks."},{"cited_title":"Markovich, D","cited_arxiv_id":null,"evidence_quote":"Provides the thylakoid stack parameters (interdisk separation, ionic strength ranges) and the Hamaker constant used to compare the electrostatic attraction with van der Waals."},{"cited_title":"Mustárdy, K","cited_arxiv_id":null,"evidence_quote":"Gives the Poisson-Boltzmann equation and the analytic single-surface potential used to evaluate the electrolyte contributions outside the stack in the grand potential."}],"review_version":1}