{"id":"fff259ed-e471-409f-97b5-3a0399e6be2a","arxiv_id":"1908.11739","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"A continuous field theory model predicts that protein beta strands have a universal length of about 12 Å, with their abundance controlled by a single parameter.","lead":"This paper models protein chains as continuous curves using a physics theory called the Abelian Higgs model, suggesting that helices are the natural ground state and beta strands are metastable segments. It derives a predicted average length for beta strands of about 12 Å, but the claim rests on parameters fitted to the same data used for validation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 12 Å beta-strand prediction is not independently tested: F and epsilon are fitted from the same beta-strand Figure 2, R_beta is never compared with measured strand lengths, and the derivation is in an unavailable companion paper.","rationale":"The Reader correctly highlighted the fragility of identifying beta strands with zero-curvature configurations and of measuring ribbon torsion rather than backbone torsion. My stress test points to a closely related but more direct gap: the paper's central quantitative prediction, R_beta ~ 12 Å, is not compared with independent beta strand length data, and the derivation of the potential and parameter mapping is delegated to an unavailable companion paper. This circularity and missing support reinforce the REJECT verdict. The proposed check is a direct way to decide whether the \"universal size\" claim has real empirical content or is an artifact of fitting the same data that the prediction is supposed to explain.","tokens_in":5950,"tokens_out":21069,"duration_ms":191862,"concrete_test":"Assemble an independent set of beta strand lengths from a protein structure database, using the same inclusion or exclusion rule as Fig. 2, and compare the distribution to R_beta ~ 12 Å with uncertainty propagation. If the independent distribution does not peak near 12 Å, or if the spread is comparable to the difference, the universal prediction is unsupported. A complementary check is to refit F and epsilon using only alpha-helix points and recompute Eq. (7); if the value shifts by more than a few angstroms, the circular use of strand data is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (7) is presented as \"a prediction of the universal size of the length of the beta strand in the model.\" The prediction is not independent of the data it is supposed to explain. F = 0.70 Å^-1 and epsilon = 1.5 Å^-1 are obtained by fitting Eq. (2) to the Fig. 2 distribution of helix and strand points (Eq. 6); the same Fig. 2 is then cited as evidence that relation (2) is phenomenologically valid. The resulting R_beta ~ 12 Å is computed from these fitted parameters, but no measured distribution of beta strand lengths is shown for comparison, and no error bars are given. If the fit already absorbs the strand geometry through F and epsilon, a \"prediction\" that comes out near a typical strand length is circular. The missing derivation of potential (4), the parameter mapping, and the fit methodology are explicitly deferred to ref. 10, which the footnote reports is not available on arXiv. Thus the paper's own text flags that the central quantitative claim lacks the support needed to check it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a continuous one-dimensional Abelian Higgs model with a Chern-Simons term and a Proca mass term as an effective field theory for protein backbone geometry. After integrating out the torsion, it claims a reduced energy functional whose derivative-free part can be rewritten as a two-minimum potential (4), with helices as the ground state, beta strands as metastable zero-curvature configurations, and loops as kinks interpolating between helical minima. The authors fit the model parameters to curvature-torsion data for helices and beta strands, and from the fitted parameters they derive a universal beta-strand length R_beta ≈ 12 Å, stated as a prediction of the model. They further argue that a single parameter kappa1 controls the abundance of beta strands, acting like a chemical potential.","tokens_in":6107,"tokens_out":21313,"duration_ms":184478,"significance":"If the central derivation were supplied and the 12 Å prediction independently tested, the paper would offer a genuinely attractive statement: a two-minimum effective potential with a single controlling parameter that unifies helices, beta strands, and loops in one universality class. The qualitative mapping of kinks, sphalerons, and zero-curvature segments to secondary-structure motifs is concrete and pedagogically useful, and the claimed R_beta ≈ 12 Å is a crisp, falsifiable quantitative prediction. However, in its current form the letter does not establish the derivation of the effective potential, and the empirical validation is not independent of the fitted parameters. The paper therefore cannot yet support its main quantitative claim.","major_comments":[{"comment":"The reduction from the gauge-invariant energy (1) to the algebraic relation (2) and the reduced functional (3) is not shown and does not follow from Eq. (1) as written. The torsion-dependent part of Eq. (1) is -∫ F τ̂ ds + (1/2ε²)∫(τ̂ - η')² ds. Varying with respect to the physical torsion τ = τ̂ - η' gives τ = F ε², a κ-independent constant, rather than the κ-dependent relation τ = F/(κ² + ε²) needed to produce the F²/(κ² + ε²) term in Eq. (3). Either Eq. (1) is missing a κ-dependent torsion mass or kinetic term, or Eq. (2) is simply asserted. Since Eq. (3) and all subsequent analysis use Eq. (2), this is a load-bearing gap.","section":"Eqs. (1)–(3)"},{"comment":"The paper states that Eq. (4) is a rewrite of the derivative-free part of Eq. (3), but this is not algebraically correct as printed. Multiplying the potential in Eq. (3), U(κ) = ½(-m²κ² + λκ⁴) - F²/[2(κ² + ε²)], by 2(κ² + ε²) and comparing with Eq. (4) forces the constant-term condition -F² = λ κ₀⁴ κ₁², which is incompatible with real parameters and with the stated assumption 0 ≤ κ₁² ≤ ε². An unstated additive constant and additional parameter relations would be required to make the two potentials coincide. The 'two different parameterizations' relation between λ and F invoked in Eq. (7) is also deferred to the companion Ref. 10, which the footnote reports is not available on arXiv. Because Eq. (4) is the basis of the kink solutions, Fig. 3, and Eq. (7), the central derivation is not checkable in this letter.","section":"Eq. (4)"},{"comment":"The R_β ≈ 12 Å result is not an independent prediction. The values F = 0.70 Å⁻¹ and ε = 1.5 Å⁻¹ in Eq. (6) are obtained by fitting relation (2) to the helix and beta-strand points in Fig. 2, the same data whose compatibility with Eq. (2) is presented as validation. Equation (7) then evaluates R_β from these fitted values together with κ₀ from Eq. (5), but no measured distribution of beta-strand lengths is shown and no error bars or sensitivity analysis are given for R_β ≈ 12 Å. As a result, the paper does not demonstrate that this is a falsifiable prediction or that it agrees with protein phenomenology beyond the data already used to fix the parameters.","section":"Eqs. (6)–(7) and Fig. 2"},{"comment":"The empirical support for Eq. (2) is weakened by the use of a different geometric observable for beta strands. The fields κ and τ in Eqs. (1)–(3) are properties of the backbone curve, and Eq. (2) is a relation between backbone curvature and backbone torsion. For beta strands the authors state that they measured 'the torsion of the ribbon, rather than that of the backbone chain.' At κ → 0 the Frenet torsion of a nearly straight backbone is not a well-defined continuous quantity, whereas the ribbon torsion is a different quantity; no justification is given for using ribbon torsion to test a relation derived for the backbone. Consequently, the fitted parameters in Eq. (6) and the compatibility claim for Fig. 2 are not established as tests of the model.","section":"Identification of β strands and Fig. 2"}],"minor_comments":[{"comment":"The companion paper is explicitly reported to be 'still on hold by the arXiv moderators' as of the submission date; the authors should either make the derivation available as supplementary material or state clearly which results in the letter depend on it.","section":"Reference 10"},{"comment":"The displayed formula for R_β is garbled in the typesetting; it should be clarified whether the integrand is 1/[√λ κ₀ ε √(k² + κ₁²)] and how the final rational expression in F, ε, and κ₀ is obtained.","section":"Eq. (7)"},{"comment":"The figure caption and text do not specify how the beta-strand points were obtained from protein structures, which proteins were used, or how the ribbon torsion was computed; these details are needed for reproducibility.","section":"Fig. 2"},{"comment":"The term 'universality class' is used without a precise definition; the authors should specify what can and cannot vary within the class and what observables are universal.","section":"General"}],"recommendation":"reject","confidential_remarks":"The central quantitative claim rests on derivations deferred to an unavailable companion paper, and the paper itself contains an apparent inconsistency between Eq. (1) and Eq. (2) as well as a non-equivalence between Eq. (3) and Eq. (4). The empirical validation is not independent. I would be willing to reconsider a substantially revised version that supplies the derivation, corrects the model equations, and tests R_β against measured strand lengths."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the headline claim—a universal beta-strand length of about 12 Å—is not independently tested. F and ε are fitted to the same Figure 2 that is then offered as evidence that relation (2) holds, and the fitted values feed directly into Rβ via Eq. (7). The paper even says the derivation is in a companion paper (Ref. 10) not available at submission. So the quantitative centerpiece is unsupported as a prediction.\n\nWhat the paper does well: the effective field theory picture is coherent, and the reinterpretation of κ1 as a chemical potential controlling strand abundance is a nice idea. The kink solutions and loop size dependence are worked out explicitly, with credit given to prior kink literature (Refs. 19–21). The paper is honest about what it did: it explicitly notes that the ribbon torsion, not the backbone torsion, was measured for strands, and it states that Ref. 10 contains the details. That transparency is worth something.\n\nSoft spots: the circularity is real and load-bearing. Figure 2 is used both to fit F and ε and to claim compatibility with relation (2); the 12 Å value is computed from those fitted parameters and is never compared to measured strand-length distributions. No error bars are given, so the \"prediction\" has no stated uncertainty. The exclusion of the 2pne protein's strands is post hoc—the authors say those strands do not fit the universality class, which is a defensible modeling decision but weakens the empirical test. And the identification of beta strands with zero-curvature configurations is assumed rather than argued, which matters because the fit depends on measuring the ribbon's torsion. These are not minor caveats; they bear on whether the central claim is established.\n\nWhere I'd push back on the skeptical notes: the model itself is not nonsense—it is a specific, falsifiable framework, and the math up to Eq. (7) is straightforward. The problem is the evidence, not the idea. Also, the paper's self-citation to Ref. 10 is not a flaw per se; it becomes a flaw only because the companion paper is central and unavailable to the reader.\n\nBottom line: this paper deserves a serious referee. It is a plausible and interesting framework, but the current letter does not support the universal-prediction claim. A referee should require either an independent test (e.g., measured strand lengths with error bars) or a self-contained derivation and fit methodology. If that is added, this could be a useful contribution to coarse-grained protein modeling.","headline":"An interesting effective-field-theory sketch with a circular quantitative centerpiece: the 12 Å beta-strand 'prediction' comes from parameters fitted to the same data it claims to describe.","tokens_in":6718,"tokens_out":2198,"would_cite":false,"duration_ms":21013,"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":"The paper predicts that beta strands have a universal length of about 12 Å, produced by a two-minimum field theory of protein geometry.","keywords":["protein secondary structure","effective field theory","Abelian Higgs model","beta strand length","protein geometry","Chern-Simons term","kink solitons"],"falsifier":"Take a large, structurally diverse set of high-resolution protein structures and measure, for every $\\beta$ strand, the ribbon torsion and the local curvature of the fitted continuous curve. If the $(\\kappa,\\tau)$ pairs do not follow $\\tau = F^2\\kappa/(\\kappa^2+\\epsilon^2)$ with the published values of $F$ and $\\epsilon$, or if the typical isolated strand length is not near 12 Å, the model's central prediction is falsified.","tokens_in":5650,"feed_emoji":"🧬","tokens_out":5382,"duration_ms":50404,"temperature":0.7,"pith_summary":"This paper argues that the large-scale geometry of proteins can be described by a continuous one-dimensional field theory with only a few parameters. In this theory, helices emerge as the stable ground state, beta strands are nearly straight, twisted configurations, and the abundance of strands is set by one parameter acting like a chemical potential. The model yields a quantitative prediction: a universal beta-strand length of about 12 Å, independent of the remaining parameters. If the model holds, many protein secondary-structure features would follow from a single symmetry-breaking mechanism rather than from the detailed chemistry of each amino acid.","feed_headline":"Protein field theory predicts beta strands near 12 Å","feed_subtitle":"If the prediction holds, helices are the ground state and one parameter controls whether strands form.","key_machinery":"The load-bearing object is the one-dimensional Abelian Higgs model with a Chern-Simons term and a Proca mass term, Eq. (1), written for a complex curvature field $\\hat\\kappa = \\kappa e^{i\\eta}$ and a torsion field $\\hat\\tau$. Gauge invariance removes $\\eta$ by shifting the torsion, giving the effective potential $V(\\kappa) = \\lambda(\\kappa^2-\\kappa_0^2)^2(\\kappa^2+\\kappa_1^2)/(2(\\kappa^2+\\epsilon^2))$ with two minima. The special relation $\\tau = F^2\\kappa/(\\kappa^2+\\epsilon^2)$ emerges from integrating out torsion and is the empirical tie to real proteins. Static solutions of this potential, namely constant-curvature helices, zero-curvature strands, kinks as loops, and sphalerons as hairpins, carry the paper's interpretations, and the size of the flat step inside the kink gives the universal strand length.","core_discovery":"The paper's central claim is that a continuous Abelian Higgs model with a two-minimum potential captures universal geometric features of protein secondary structure. In the model, a complex curvature field has a symmetry-broken ground state with constant curvature and torsion, meaning a helix, while zero-curvature configurations correspond to $\\beta$ strands. Integrating out the torsion field produces an effective potential for curvature whose minima and kink solutions reproduce helices, loops, strands, and hairpin-like sphalerons. Fitting the curvature-torsion relation to real structures gives $F = 0.70\\ \\text{\\AA}^{-1}$ and $\\epsilon = 1.5\\ \\text{\\AA}^{-1}$; with the helix parameters $\\kappa_0 \\simeq 1.6\\ \\text{\\AA}^{-1}$ and $\\tau_0 \\simeq 0.15\\ \\text{\\AA}^{-1}$, the model predicts a universal $\\beta$-strand length $R_\\beta \\simeq 12\\ \\text{\\AA}$. The paper further claims that $\\beta$-strand abundance is controlled by the parameter $\\kappa_1$, which can either suppress strands completely or make them abundant.","pith_inferences":["If the 12 Å prediction is correct, beta-strand length distributions in structural databases should show a floor near this value across fold families; this could be checked directly without any new model assumptions.","The chemical-potential interpretation of $\\kappa_1$ suggests a concrete test: protein families in different cellular environments or with different amino-acid compositions should differ statistically in beta-strand abundance in the direction the parameter prescribes.","The model's curvature-torsion relation implies an anticorrelation between ribbon curvature and torsion along strands; re-measuring this relation on a large independent structure set would be a sharper test than the single fitted figure.","Extending the model to position-dependent parameters would turn it into a generative model for secondary-structure patterns, allowing comparison of predicted loop-size distributions with observed ones."],"forward_implications":["Helices are the generic ground state of the model, while beta strands appear only when the chemical-potential parameter $\\kappa_1$ permits them, making strand abundance a tunable rather than fixed property.","The predicted beta-strand length $R_\\beta \\approx 12\\ \\text{\\AA}$ is parameter-independent in the model, so it should appear as a universal scale across unrelated proteins.","Long loops interpolating between helices are kink solutions whose size grows logarithmically as $\\kappa_1 \\to 0$, matching the expectation that low-$\\kappa_1$ proteins have extended near-straight segments.","Sphaleron solutions give a natural field-theoretic counterpart of beta hairpins, predicting higher-curvature loop regions connecting beta strands.","The model's few parameters define a universality class; discrete and inhomogeneous generalizations should preserve the geometric relations rather than reproduce exact atomic coordinates."],"supporting_citations":[{"why":"Establishes the effective field theory approach to proteins that this paper extends to the continuous model.","marker":"[1,2]"},{"why":"Companion paper supplying the curvature-torsion fits, the ribbon-torsion measurement for beta strands, and the parameter estimates used here.","marker":"[10]"},{"why":"Provides the form of the effective energy functional for curvature and torsion fields used in Eq. (1).","marker":"[11,12]"},{"why":"Ramachandran-plot literature that motivates casting the curvature-torsion distribution as a diagnostic plot for secondary structure.","marker":"[16-18]"},{"why":"Earlier work on similar kink solutions, supporting the interpretation of loops as kinks in the model.","marker":"[19-21]"}],"fun_headline_variants":["Field theory predicts universal beta-strand length","Helix ground state, single parameter controls strand abundance","Protein curves: helical ground state, 12 Å beta strands","One knob makes beta strands abundant or absent","Continuous protein model: helices win, strands optional"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a beta strand is correctly represented by the model's nearly straight, twisted-ribbon solution, with the ribbon's torsion, not the backbone's torsion, used as the measured quantity; if this identification fails, the fitted parameters and the 12 Å prediction lose their empirical support.","fun_headline_variants_meta":{"raw":{"variants":["Field theory predicts universal beta-strand length","Helix ground state, single parameter controls strand abundance","Protein curves: helical ground state, 12 Å beta strands","One knob makes beta strands abundant or absent","Continuous protein model: helices win, strands optional"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000609,"raw_usage":{"total_tokens":2833,"prompt_tokens":938,"completion_tokens":1895,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":1821}},"tokens_in":554,"tokens_out":1895,"duration_ms":12026,"temperature":1.0,"reasoning_tokens":1821,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:07:57.844587+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a large, structurally diverse set of high-resolution protein structures and measure, for every $\\beta$ strand, the ribbon torsion and the local curvature of the fitted continuous curve. If the $(\\kappa,\\tau)$ pairs do not follow $\\tau = F^2\\kappa/(\\kappa^2+\\epsilon^2)$ with the published values of $F$ and $\\epsilon$, or if the typical isolated strand length is not near 12 Å, the model's central prediction is falsified.","supporting_citations":[{"cited_title":"Melnikov and A","cited_arxiv_id":null,"evidence_quote":"Companion paper supplying the curvature-torsion fits, the ribbon-torsion measurement for beta strands, and the parameter estimates used here."}],"review_version":1}