{"id":"e2c40e83-48ec-483e-bcf4-d9ff39ec39aa","arxiv_id":"2607.14664","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Epigenetic memory, thresholds, and hysteresis are recast as rate-independent dissipative evolution: an energy landscape plus a resistance-to-change cost yields closed-form switch thresholds and a convergent, energy-consistent numerical scheme.","lead":"This paper builds a mathematical framework, Rate-Independent Epigenetics, that models epigenetic marks as rate-independent dissipative variables: an energy landscape plus a resistance-to-change threshold reproduces memory, threshold activation, and hysteresis, using the mathematics of plasticity and magnetism. A generalist reader might care because it is the first variational, thermodynamically-consistent treatment of epigenetic switching, with few parameters and explicit fal","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Rate-independence is the load-bearing empirical premise; the paper's own falsifier—frequency-dependent scarring—is untested, so the biological claim remains conditional.","rationale":"The reader's weakest assumption is the right one. I found no internal mathematical error in the derivation: the energetic formulation, Biot inclusion, saturation identity, and energy-consistency estimates are coherent, and the paper is unusually explicit about its idealisations (Secs. 7.1, 7.5) and about the distinction between energetic and balanced-viscosity switch selection (Remark 6.1). The BV uniqueness proof is a sketch and the bistable benchmark computes the energetic solution rather than the BV-selected path, but these are secondary to the fact that the framework's biological relevance depends on rate-independence, which is asserted and falsifiable yet not tested. Since the paper itself identifies the falsifier, the honest verdict is CONDITIONAL: accept the mathematical framework, condition the biological claim on the rate-independence test. This does not change the reader's verdict.","tokens_in":32369,"tokens_out":6428,"duration_ms":73135,"concrete_test":"Take a concrete epigenetic system (e.g., hypoxia-driven DNA methylation, Ref. [36]) and apply an intermittent stimulus of fixed amplitude with two or more distinct frequencies—equivalently, reparametrize the same loading path S(τ) at speeds differing by a factor of at least 10 while keeping endpoints fixed. Measure the residual scar q(T) and the hysteresis loop area. If either depends on frequency or sweep rate, RIE's rate-independent predictions fail and the viscous extension of §7.4 becomes mandatory; if both are invariant under time reparametrization, the premise is supported. A lower-cost first pass is to fit the play-operator or bistable model to existing time-course data and test whether inferred thresholds and scars are rate-independent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 7.1 states the premise and its falsifier: the pure rate-independent model has no intrinsic time scale, cannot represent relaxation, and predicts no ratcheting and no frequency dependence; 'an experimentally observed dependence of the epigenetic scar on the frequency of an intermittent stimulus at fixed amplitude would signal a genuine viscous (rate-dependent) contribution.' No such experiment or dataset is presented. All the quantitative biological outputs—the play-operator lock-in criterion ρ>ℓmax/2 (Eq. 27), the bistable jump at ℓ=ρ (Prop. D.1), the geometric-area dissipation law (§3.3), and the no-ratcheting/frequency-independence prediction (§7.1)—are rate-independent statements. If chromatin-state dynamics is not quasi-static at the timescales of interest, those thresholds and scars will not be observed, and the central claim that epigenetic change can be modelled as rate-independent dissipative evolution lacks empirical support. The mathematics remains self-consistent; the biological transfer does not.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes Rate-Independent Epigenetics (RIE), a variational framework in which an epigenetic system is specified by a triple (Q,E,Ψ): a configuration space, a stored energy E(q,S)=F(q)−q·ℓ(S), and a 1-homogeneous dissipation potential Ψ. From the energetic evolution principle (global stability (S) plus energy balance (E), Definition 2.1), the authors derive the Biot-type subdifferential inclusion (Proposition 2.2), read the energy balance and dissipation inequality as a first and second law, and prove existence of energetic solutions (Theorem 4.3), uniqueness under convexity (Theorem 4.4), and existence and scalar uniqueness of balanced-viscosity solutions (Theorems 4.6, 4.7). They construct an incremental-minimisation integrator, prove convergence (Theorem 5.1) and a first-order global energy-consistency estimate (Theorem 5.2), and test the method on the scalar play operator, a piecewise-quadratic double well, and a nonlinear convex example. The paper is explicit that rate independence is an idealisation (Section 7.1) and that the potentials and thermodynamic statements are effective (Section 7.5).","tokens_in":32515,"tokens_out":10737,"duration_ms":125034,"significance":"If the advertised results hold, this is a useful and largely rigorous transfer of rate-independent-system theory to epigenetics. The paper's main strengths are its honest delineation of what is proved versus what is assumed; the explicit falsifier in Section 7.1 (frequency-dependent scarring would signal viscous behaviour); the careful separation of the proven energy-consistency estimate (Theorem 5.2) from the merely empirical state-error rates (Proposition C.4, Remark 5.3); and the closed-form energetic jump threshold ℓ=ρ for the bistable benchmark (Proposition D.1). These components give the framework a solid applied-mathematics core even though its biological validation is deferred.","major_comments":[{"comment":"The biological claim is conditional on a premise that is not tested. The abstract says RIE 'models epigenetic change', but the only empirical content is the falsifiable prediction in Section 7.1 that frequency-dependent scarring would signal a viscous contribution; no frequency-resolved or otherwise discriminating experiment is reported. All quantitative outputs (Equation (27), Proposition D.1, Section 3.3) are rate-independent statements. The authors are honest about this, but the manuscript should either present/analyse such an experiment or explicitly restrict the claims in the abstract and title to a quasi-static mathematical skeleton rather than a validated model of epigenetic response.","section":"Section 7.1 / Abstract"},{"comment":"The proof of the scalar uniqueness theorem is not a proof. It asserts without derivation that the balanced-viscosity jump condition follows 'the path of steepest descent of q↦E(q,S(t∗))+ρ|q−q−|' and selects the 'first local minimiser of E(·,S(t∗))'. These two descriptions are inconsistent: stationary points of E+ρ|·| satisfy E_q=±ρ, not E_q=0. Since Theorem 4.7 is advertised in the abstract, the jump-selection mechanism needs a rigorous derivation from the vanishing-viscosity construction, or a precise reference to a theorem that covers this scalar setting; as written, the uniqueness claim is unsupported.","section":"Appendix B.3, Theorem 4.7"},{"comment":"Existence of the ε-regularised viscous solution is attributed to Brézis [6], which concerns maximal monotone operators. For the bistable landscapes of central interest, ∂_q E is not monotone, so the cited theorem is not directly applicable. A finite-dimensional argument using the single-valued Lipschitz inverse of ∂Ψ+εI would close the gap, but as written the construction of BV solutions is not fully justified.","section":"Appendix B.3, Theorem 4.6, Step 1"},{"comment":"The bistable numerical benchmark does not validate the balanced-viscosity solution. The incremental scheme (20) minimises globally and therefore computes the energetic solution, with the energetic jump threshold ℓ=ρ (Proposition D.1), whereas the physically motivated BV solution would have the local threshold ℓ=ka+ρ (Remark 6.1). The paper explicitly does not carry out the BV construction numerically. Consequently Section 6.2 validates only the energetic-solution branch, not the uniqueness/selection mechanism of Theorem 4.7 that is central to the framework's physical interpretation. This limitation should be stated more prominently than a remark, and ideally a local/BV implementation should be added or the scope reduced.","section":"Section 6.2 / Remark 6.1"}],"minor_comments":[{"comment":"The notation S(t) is used both for the loading path and for the stability set, which is confusing. Rename the stability set, e.g. Σ(t).","section":"Definition 4.2"},{"comment":"The phrase 'global energy-consistency estimate' in the abstract is stronger than Theorem 5.2, which assumes q∈AC or q∈W^{1,∞}. Across jumps, the estimate is only observed numerically (Section 6.2.3). Please state this regularity restriction explicitly in the abstract or in the statement of Theorem 5.2.","section":"Section 5.3 / Abstract"},{"comment":"In the unloading branch, the formula q(t)=min{qmax, u(t)+r} suppresses the lock value at the turning point. The closed form is correct but the piecewise ranges could be written more transparently, especially for readers unfamiliar with the play operator.","section":"Equation (26)"},{"comment":"The sup-norm state error is reported only away from the jump; the full sup-norm error is O(1) because of the one-step jump location error. This is consistent with Proposition C.4, but should be stated in the caption of Figure 9 to avoid misreading.","section":"Section 6.2.3"},{"comment":"The paper contains no data- or code-availability statement. Given the reproducible benchmarks and closed-form solutions, a small reproducibility remark would increase the value of the numerical sections.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is mathematically serious and unusually honest about its own limitations. My main reservations are: (i) the abstract-level biological claim outruns the evidence, and (ii) one advertised theorem (Theorem 4.7) lacks a rigorous proof, while the proof of Theorem 4.6 cites an inapplicable monotone-operator result. These are fixable within the manuscript's scope. The paper is more applied mathematics than empirical biology; whether that fits physics.bio-ph depends on the journal's tolerance for theoretical framework papers with no biological data. I would not reject on mathematical grounds."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper is best read as a variational mechanics paper with an epigenetic application, not as an empirical biology paper. It transfers energetic/balanced-viscosity rate-independent theory to a scalar internal-variable model of chromatin state, and it does that transfer carefully. The novelty is real: the closed-form thresholds (Prop D.1, ℓ=ρ vs local ℓ=ka+ρ), the node-exactness analysis (Prop C.4), and the two-sided return-map reduction are genuinely new in this setting. The derivations are clean, the proofs I checked in the appendices are rigorous (Prop 2.2, energy estimates, convergence, consistency), and the paper is unusually candid about what it doesn't do—Sec 7.1 spells out the frequency-dependence falsifier, Sec 7.5 warns that the thermodynamics is effective, Remark 5.3 warns that nodal exactness is not general. That honesty earns credit.\n\nSoft spots, in proportion. First, the central premise—rate independence as quasi-static limit—rests on no biological data. The paper itself supplies the falsifier: if epigenetic scar depends on stimulus frequency at fixed amplitude, the skeleton fails as a biological statement. No such experiment is presented. That makes the biological significance conditional, not nonexistent. For a theory paper, that's acceptable if framed as a skeleton; the abstract slightly oversells by claiming validation on the bistable benchmark when the body carefully says the incremental scheme computes the energetic solution, not the BV-selected path. The switch-selection mechanism at jumps is therefore not demonstrated numerically.\n\nSecond, Theorem 4.7 (scalar uniqueness of BV solutions) is the load-bearing uniqueness result for the non-convex case, and its proof is a sketch: the \"first local minimiser encountered\" selection is asserted rather than proved. I don't see a clear error, but the crucial jump-selection step needs a real proof or a citation that covers this exact setting.\n\nThe sub-step loading experiment (Sec 6.1.4) is a nice negative control: it shows the state error plateau at O(1) when h≫τ while the energy defect stays bounded. That is the kind of honesty that makes the numerical section trustworthy.\n\nOverall: a solid mathematical paper, honest about its limits, with a conditional biological payoff. A serious referee should engage with it, but should push on Theorem 4.7 and on the missing empirical test. I would send it to peer review, and would cite it if I worked on variational methods for epigenetic or hysteretic biological systems.\n\nBest,\n[You]","headline":"A mathematically careful transfer of rate-independent systems to epigenetics, with honest caveats, but the biological premise is untested and the BV uniqueness proof is sketched.","tokens_in":33135,"tokens_out":1788,"would_cite":true,"duration_ms":18626,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["49J40","74C05","34G25","74N30","65M12","92C40","80A17","37D35"],"pacs":[],"model":"deepseek-v4-flash","headline":"Epigenetic change can be modelled as rate-independent dissipative evolution, with thresholds, memory, and hysteresis following from one variational principle.","keywords":["rate-independent systems","epigenetics","hysteresis","dissipation potential","energetic solutions","balanced-viscosity solutions","variational time integration","thermodynamic consistency"],"falsifier":"Measure the epigenetic scar (residual mark) after applying a fixed-amplitude, intermittent stimulus at different frequencies. The rate-independent model predicts the scar is independent of frequency; an observed dependence of the scar on frequency would signal a genuine viscous (rate-dependent) contribution and refute the pure rate-independent description of epigenetic response. Alternatively, measuring the hysteresis-loop area under loading cycles of different speeds and finding that the area changes with sweep rate would similarly falsify the rate-independence hypothesis.","tokens_in":32123,"feed_emoji":"🧬","tokens_out":2590,"duration_ms":30487,"temperature":0.7,"pith_summary":"This paper tries to establish that epigenetic changes—heritable, reversible modifications of chromatin state—can be captured by the theory of rate-independent systems. A model is fixed by a triple: a space of epigenetic configurations, a stored-energy landscape depending on an external micro-environmental loading, and a 1-homogeneous dissipation potential. Postulating an energetic evolution principle (global stability plus energy balance) generates all governing equations without further hypotheses, and the resulting solutions satisfy exact energy conservation and non-negative, minimal dissipation by construction. The paper further proves existence, uniqueness under convexity, and uniqueness of 'balanced-viscosity' solutions that resolve how a state jumps at an epigenetic switch, and it constructs a numerical integrator that converges and is first-order consistent in the energy balance. In the scalar bistable case, the framework yields closed-form thresholds—such as a jump at ℓ=ρ and permanent lock-in when the dissipation threshold exceeds ρ>ℓ_max/2—so a sympathetic reader would care because the paper supplies a thermodynamically consistent skeleton that turns epigenetic phenomenology into computable, testable predictions.","feed_headline":"Epigenetic thresholds arise from one variational principle","feed_subtitle":"A new rate-independent framework builds thermodynamic consistency into every epigenetic model, yielding exact switch and lock-in predictions","key_machinery":"The central object is the rate-independent system triple (Q,E,Ψ): a closed convex state space, a stored energy E(q,S)=F(q)−q·ℓ(S) that splits into a configurational free energy and a coupling to the loading, and a positively 1-homogeneous dissipation potential Ψ whose subdifferential at zero defines an elastic domain E. The argument is carried by the energetic evolution principle (global stability (S) and energy balance (E)), the derived Biot inclusion 0∈∂Ψ(q̇)+∂E/∂q, and the saturation identity f·q̇=Ψ(q̇), which together encode exact energy conservation and minimal dissipation. Balanced-viscosity solutions, obtained as the vanishing-viscosity limit of a regularised inclusion, select the phy","core_discovery":"The central claim is that an epigenetic system can be specified by a triple (Q,E,Ψ) — a state space of epigenetic configurations, a stored-energy functional E(q,S)=F(q)−q·ℓ(S) with a micro-environmental loading S, and a convex, 1-homogeneous dissipation potential Ψ encoding resistance to remodelling. Under the energetic evolution principle (Definition 2.1), the state evolves so that the energy is globally stable at each time and the energy balance holds exactly. From these two conditions alone, the governing subdifferential inclusion (the Biot equation) follows, and the saturation identity f·q̇=Ψ(q̇) forces dissipation to be non-negative and minimal. The laws of thermodynamics are therefore","pith_inferences":["Extension: The sharpest falsifiable prediction is frequency independence: an experiment that shows the epigenetic scar after an intermittent stimulus depends on stimulus frequency at fixed amplitude would directly indicate a viscous contribution beyond the rate-independent skeleton, exactly as the paper's own Section 7.1 states.","Extension: The framework could be extended to vectorial epigenetic marks (n≥2), where the balanced-viscosity selection becomes genuinely path-dependent; such an extension would likely require additional numerical machinery for the jump construction.","Extension: The rate-independent description may be the leading-order macroscopic reduction of a faster, stochastic gene-regulatory dynamics; if so, the potentials F, ℓ, and Ψ could be learned from fine-scale data rather than posited, connecting RIE to structure-preserving coarse-graining approaches.","Extension: The exact bistable threshold ℓ=ρ suggests a direct experimental test: measure the loading level at which a reversible mark irreversibly flips between two phenotypes, and compare with the model's predicted dependence on the dissipation resistance ρ."],"forward_implications":["Any model built on the triple (Q,E,Ψ) is thermodynamically consistent a priori: the first and second laws hold for every choice of landscape and dissipation, so no ad-hoc evolution laws are needed.","Threshold activation, memory, and hysteresis—the defining phenomenology of epigenetic marks—emerge directly from the variational structure, and the energy dissipated over a closed loading cycle equals the area of the hysteresis loop regardless of how fast the cycle is traversed.","In the scalar bistable case, the framework yields exact analytical predictions: the state jumps between wells at the threshold ℓ=ρ, and the mark is permanently locked in whenever the dissipation resistance exceeds ℓ_max/2.","The numerical integrator converges to energetic solutions and carries a proven global first-order energy-consistency estimate, so computed epigenetic trajectories come with a quantified energy-balance error.","The model predicts the absence of ratcheting: repeating an identical loading cycle leaves the same residual state, and the maximal damage and residual scar depend only on the extreme values of the loading, not on its rhythm.","If the paper's framework is correct, it provides a route to calibrating the effective energy and dissipation potentials from macroscopic loading-unloading experiments, making the ingredients identifiable rather than purely formal."],"fun_headline_variants":["One variational principle predicts epigenetic thresholds and memory","Thermodynamically consistent framework for epigenetic switches","Rate-independent epigenetics: a rigorous mathematical model","A single variational principle yields epigenetic hysteresis"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that real epigenetic dynamics operates in the rate-independent (quasi-static) regime, where the chromatin state either equilibrates quickly relative to the variation of its micro-environment or stays effectively frozen except when a threshold is crossed; if chromatin-state change carries intrinsic rate dependence at the timescales of interest, the closed-form thresholds, hysteresis areas, and no-ratcheting predictions fail as biological statements.","fun_headline_variants_meta":{"raw":{"variants":["One variational principle predicts epigenetic thresholds and memory","Thermodynamically consistent framework for epigenetic switches","Rate-independent epigenetics: a rigorous mathematical model","A single variational principle yields epigenetic hysteresis"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000663,"raw_usage":{"total_tokens":2913,"prompt_tokens":842,"completion_tokens":2071,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":586,"completion_tokens_details":{"reasoning_tokens":2015}},"tokens_in":586,"tokens_out":2071,"duration_ms":14925,"temperature":1.0,"reasoning_tokens":2015,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T01:27:10.451027+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the epigenetic scar (residual mark) after applying a fixed-amplitude, intermittent stimulus at different frequencies. The rate-independent model predicts the scar is independent of frequency; an observed dependence of the scar on frequency would signal a genuine viscous (rate-dependent) contribution and refute the pure rate-independent description of epigenetic response. Alternatively, measuring the hysteresis-loop area under loading cycles of different speeds and finding that the area changes with sweep rate would similarly falsify the rate-independence hypothesis.","supporting_citations":[],"review_version":1}