{"id":"d099fd17-2365-43f2-aa57-a31dc49d0819","arxiv_id":"1908.06872","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A 0/1 spin Ising model on Barabasi-Albert networks shows first-order transitions and hysteresis, but reduces by the authors' own variable transformation to a classical Ising model with local fields.","lead":"The authors model gene activity as binary 0/1 spins on a scale-free Barabasi-Albert network and study phase transitions with Monte Carlo and mean-field theory. Their own mapping shows this is exactly a classical Ising model with node-dependent local fields, so the claimed novelty is a coordinate transformation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The derivation of h_c = Jm (Sec. IV B, Eqs. 26-29) is internally inconsistent: Eq. 27 contradicts the paper's own mapping in Eq. 22, and Eq. 28 with the stated mean degree k_bar = 2m gives 2Jm, not Jm; the headline quantitative claim is not derived as written.","rationale":"I independently re-derived the paper's two central analytical steps. The mapping in Eqs. 20-23 is exact: with s_i = (s'_i+1)/2, the 0/1 Hamiltonian becomes the classical Ising Hamiltonian with coupling J'_ij = JA_ij/4, local field h'_i = h/2 + (J/4) sum_j A_ij = h/2 + Jk_i/4, and a state-independent constant. So the 'modified Ising model' is the standard lattice gas, i.e., the Ising model with a degree-dependent local field on a BA network. The paper's own Eq. 23 states this equivalence, which directly undercuts the novelty claim ('To our knowledge, such an analysis ... has not been investigated so far'), while the low-temperature behavior (discontinuous magnetization jump, hysteresis) is standard for such lattice-gas-type systems. However, the paper fails to use its own mapping consistently. Eq. 22 states h' = h/2 + (1/2) sum_j J_ij (factor-2 error, plus sign); Eq. 27 states <h'_i> = h/2 - (J/2)k_i (minus sign, contradicting Eq. 22, and still the wrong prefactor relative to the exact mapping). Eq. 28, h_c/2 - (J/2)k_bar = 0, with the paper's own k_bar = 2m gives h_c = 2Jm; the stated h_c = Jm is arithmetic with the wrong degree. Two compensating errors (sign and degree substitution) are required for the paper's equations to produce Jm. That is the load-bearing concern: the headline quantitative result is not derived as written. In fairness, the result itself is not clearly wrong. A T->0 solution of the paper's own central mean-field equation (Eq. 15) shows the M = 1 branch is a fixed point only for h > -Jm (using the BA minimum degree k_min = m), so the reverse-sweep loss of the ordered branch occurs at |h| = Jm; the corrected mapping with the exact prefactor J/4 similarly gives h_c = -Jm. So an independent, correct derivation supports the numerical value. This is why I do not escalate the mathematics alone to a stronger claim: the defect is in the argument, not evidently in the result. The verdict REJECT still stands because the paper presents itself as having derived h_c = Jm (abstract, Sec. IV B, conclusions) when its derivation is invalid; the central model is a one-line reparametrization of a known model, contradicting the novelty framing; and the Monte Carlo evidence (Figs. 4-8) has no error bars, no convergence diagnostics, no code or data availability, and an incoherent sampling statement ('3x10^4 MC steps ... average 10 spin flips per spin' is inconsistent with N = 5x10^3 by orders of magnitude under either reading of 'MC step'). These are more than cosmetic; the paper cannot be checked or built upon as submitted. The reader's weakest_assumption (annealed versus quenched disorder) is real but secondary: at T->0 the h_c = Jm result depends on the degree sequence (k_min), which is realization-independent for BA networks, so the annealed <A_ij> approximation is not the main vulnerability. The reader's rationale, which flags the sign error in Eq. 27 and the 2Jm arithmetic of Eq. 28, is the concern I confirm and sharpen.","tokens_in":14338,"tokens_out":30429,"duration_ms":277685,"concrete_test":"Recompute h_c from the paper's own equations: (i) replace Eq. 27 by the exact-mapping result <h'_i> = h/2 + (J/4)k_i and substitute k_bar = 2m into Eq. 28; (ii) solve the T->0 limit of Eq. 15, M = (1/N) sum theta(h + JMk_i), on the actual degree sequence of the simulated BA network (N = 5000, m = 5), locating the fields at which the M = 0 and M = 1 branches appear and disappear. Then fix one operational definition of h_c (e.g., the midpoint between forward and reverse jump fields) and compare the candidate values Jm, 2Jm, Jm/2 with the MC data behind Fig. 8A-B for J = 1-5 and m = 3-7. If the branch-boundary field equals Jm under the same definition, the result survives but Sec. IV B must be rewritten; if it equals 2Jm or Jm/2, the central quantitative claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative result, h_c = Jm, is not established by the derivation in Sec. IV B. Three internal inconsistencies. (1) The exact mapping in Eqs. 20-23 gives the local field h'_i = h/2 + (J/4) sum_j A_ij = h/2 + Jk_i/4 (expand (s'_i+1)(s'_j+1)/4 and collect the linear s'_i terms). The paper's Eq. 22 instead states h' = h/2 + (1/2) sum_j J_ij, and Eq. 27 then flips the sign: <h'_i> = h/2 - (J/2)k_i. So Eq. 27 contradicts Eq. 22 in sign and contradicts the exact mapping in both sign and prefactor. (2) Eq. 28 sets h_c/2 - (J/2)k_bar = 0; with k_bar = 2m, which the paper itself computes one line below, this gives h_c = 2Jm, not Jm. The conclusion h_c = Jm in Eq. 29 follows only if the minimum degree m is silently substituted for k_bar, contradicting the text that defines k_bar = 2m. (3) The paper never defines how h_c is read from the hysteresis sweeps behind Fig. 8 (forward jump, reverse jump, or loop midpoint), so the claimed agreement with simulations is not checkable. This does not prove the result false: a T->0 solution of the paper's own Eq. 15 shows the M = 1 branch is lost at h = -Jm (using k_min = m), and the corrected mapping also gives |h_c| = Jm. But the derivation as written is invalid, and the headline claim rests on equations that do not chain together.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a modified Ising model with 0/1 spins (gene-type spins) on Barabasi-Albert networks and studies its phase behavior via Monte Carlo simulations and mean-field theory. The authors report a discontinuous first-order transition with hysteresis under an external field, derive a mean-field self-consistency equation, and claim a critical field h_c = Jm, where J is the coupling and m the number of preferential-attachment links. They also map the 0/1 model to a classical +/-1 Ising model on the same network. The stated goal is to provide a framework for translating statistical-mechanics tools to biological networks.","tokens_in":14664,"tokens_out":10962,"duration_ms":96547,"significance":"If the central claims hold, the paper would provide a simple, analytically tractable example of a binary-state network model with abrupt, history-dependent transitions on scale-free topologies, and a concrete prediction h_c = Jm that is testable in simulations. The model is a natural extension of the Ising model to asymmetric states, and the mean-field self-consistency equation and the exact spin mapping are useful ingredients. The qualitative Monte Carlo observations of a first-order transition and hysteresis are plausible and align with the T to 0 solution of the paper's own mean-field equation. However, the quantitative derivation of h_c = Jm contains internal sign, prefactor, and arithmetic inconsistencies, and the numerical validation lacks error bars, multiple network realizations, and a clearly defined critical-field estimator. The contribution is therefore not yet established as written.","major_comments":[{"comment":"Equation (27) is inconsistent with the mapping established in Eqs. (20)-(22). Expanding Eq. (21) with s_i=(s'_i+1)/2 gives the exact local field h'_i = h/2 + (1/4) sum_j J_ij = h/2 + J k_i/4, whereas Eq. (22) writes h' = h/2 + (1/2) sum_j J_ij and Eq. (27) changes the sign to h/2 - (J/2) k_i. Since Eq. (27) is the input to the critical-field condition in Eq. (28), the derivation of h_c does not follow from the mapping as written.","section":"IV B, Eq. (22) and Eq. (27)"},{"comment":"The step from Eq. (28) to Eq. (29) is arithmetically inconsistent with the definition of the mean degree given in the same paragraph. The paper states that k_bar = 2m, so h_c/2 - J k_bar/2 = 0 gives h_c = 2Jm, not h_c = Jm. The stated conclusion requires silently replacing the mean degree 2m by the minimum degree m, which is not justified. In addition, the paper does not specify whether h_c in Fig. 8 is read from the forward jump, the reverse jump, or the midpoint of the hysteresis loop, so the claimed agreement with Monte Carlo data is not checkable. The final result may be repairable (a zero-temperature analysis of the paper's own Eq. (15), using the minimum degree, also gives |h_c| = Jm), but the derivation as written is invalid.","section":"IV B, Eqs. (28)-(29)"},{"comment":"The Monte Carlo evidence for the first-order transition and for the h_c = Jm scaling is statistically weak. The simulations use a single Barabasi-Albert realization of size N = 5e3, equilibration of 2e4 Monte Carlo steps, and a production run of 3e4 steps that the text itself estimates to allow only about 10 spin flips per spin. No error bars, no averaging over network realizations, and no thermalization or autocorrelation diagnostics are reported. Because the mean-field calculation uses the annealed (ensemble-averaged) adjacency matrix while the simulations use one quenched realization, the quantitative comparison in Fig. 8 is not established. The authors should provide results averaged over multiple network realizations, with a clearly defined protocol for extracting the critical field from the hysteresis sweeps.","section":"III (Simulations) and Fig. 8"}],"minor_comments":[{"comment":"The effective field in Eq. (10) is written as h + Jm k_i with a lower-case m, but the same quantity appears as h + JM k_i in Eqs. (13)-(15); since m already denotes the number of preferential-attachment links, the symbol in Eq. (10) should be the order parameter M.","section":"IV A, Eq. (10)"},{"comment":"The constant term in the mapped Hamiltonian is missing a factor of N: it should be -sum_ij J_ij/8 - hN/2, not -h/2.","section":"IV B, Eq. (22)"},{"comment":"The sentence 'when h > Jm then m > 1/2, however, for h < Jm, m > 1/2' appears to use the lower-case m for the order parameter and is self-contradictory; the symbols and inequalities should be corrected.","section":"IV A, discussion after Eq. (18)"},{"comment":"The axis label |B_c| is used although the text and equations denote the critical field by h_c; B_c is never defined.","section":"Fig. 8"},{"comment":"The reference list contains malformed entries (e.g., ', Kumar, R.' with a leading comma) and inconsistent spellings (e.g., 'Dorogovtsev' vs 'Godtsev'); the list should be cleaned before publication.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The qualitative phenomenon reported here is plausible and the h_c = Jm result may well be correct, but the manuscript currently does not support it: the derivation in Sec. IV B has concrete algebraic errors, and the Monte Carlo validation needs to be redone with proper averaging and a defined critical-field estimator. These are substantial but fixable issues, so I recommend major revision rather than rejection. I would also ask the editor to ensure that the revised manuscript includes reproducible simulation details (number of realizations, sweep protocol, error bars)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper presents a 0/1 spin Ising model on Barabási-Albert networks as a new framework for gene-state transitions, claims a first-order transition with hysteresis and a critical field h_c = Jm. The qualitative picture is plausible, but the paper's own mapping (Sec. IV B) shows the model is just the classical Ising model with a node-dependent field—i.e., the lattice gas / random-field Ising model. That is not new; the 0/1 to ±1 transform is textbook. What is new, the specific prediction h_c = Jm, is not actually derived.\n\nThe mean-field equation Eq. 15 is correct and standard. But the critical-field derivation in Sec. IV B falls apart. The local field in the mapped model should be h_i' = h/2 + J k_i/4; Eq. 22 gives h_i' = h/2 + J k_i/2. Then Eq. 27 flips the sign and uses h/2 - J k_i/2. Either way, Eq. 28 with kbar = 2m gives h_c = 2Jm, not Jm. The conclusion relies on silently replacing the mean degree with the minimum degree. The paper never says how h_c is read from the hysteresis sweeps, so the agreement with Fig. 8 is not checkable.\n\nThe numerics are also weak: roughly 10 flips per spin, single network realization, no error bars, no convergence diagnostics. For a first-order transition with hysteresis, that's far too little. No code or data are provided.\n\nCredit where due: the paper is clearly written, the qualitative hysteresis observation is consistent with known random-field behavior, and the idea of using such models for gene networks is reasonable. But the central quantitative claim is not established, and the novelty claim does not survive the mapping.\n\nI would not send this to peer review in its current form. It deserves a desk reject with instructions to fix the mapping and the h_c derivation, redo the simulations with proper thermalization and error bars, and clarify how h_c is defined. If they do that, a resubmission might be worth a look.","headline":"The paper's 0/1 spin model is just the lattice gas / random-field Ising model, and its central result h_c = Jm is not actually derived as written; the Monte Carlo evidence is far too thin to save it.","tokens_in":15276,"tokens_out":11845,"would_cite":false,"duration_ms":105403,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A 0/1 gene-spin Ising model on a Barabási-Albert network undergoes a first-order phase transition with hysteresis, at a critical field $h_c=Jm$.","keywords":["modified Ising model","gene-type spins","Barabási-Albert network","first-order phase transition","hysteresis","mean-field approximation","Monte Carlo simulation","systems biology"],"falsifier":"Perform the hysteresis protocol of Fig. 5 on several independent Barabási-Albert realizations with $N=5\\times10^3$, $m=5$, $J=1$, and record the field at which the order parameter jumps. If the jump is located at $h\\approx 2Jm=10$ instead of $h\\approx Jm=5$, or if its location varies noticeably between realizations of the same $(m,J)$, the annealed mean-field critical-field prediction is contradicted.","tokens_in":14083,"feed_emoji":"🧬","tokens_out":8966,"duration_ms":80274,"temperature":0.7,"pith_summary":"The paper proposes a modified Ising model for gene regulatory networks: each node carries a 'gene-type' spin that is either inactive (0) or active (1), placed on a Barabási-Albert scale-free network. Working in the canonical ensemble, the authors argue that this asymmetric two-state system behaves like a classical ferromagnetic Ising model in an external field: as the field $h$ is swept at low temperature, the order parameter $M$ jumps discontinuously between the inactive and active states, a first-order transition, and the system exhibits hysteresis, so the state depends on the history of $h$. The central quantitative claim is that the critical field is approximately $h_c=Jm$, where $J$ is the coupling strength and $m$ the number of preferential-attachment links per node, a result derived from a mean-field solution and matched by Monte Carlo simulations. If correct, the model provides a statistical-mechanics picture of normal-to-disease transitions in gene networks as abrupt, history-dependent collective switches rather than gradual changes.","feed_headline":"Binary gene networks flip abruptly at a critical field","feed_subtitle":"A 0/1 Ising model on scale-free wiring shows first-order transitions and history-dependent hysteresis.","key_machinery":"The load-bearing object is the modified Ising Hamiltonian with spins restricted to $\\{0,1\\}$ on a Barabási-Albert network, together with two analytic tools. The first is the annealed mean-field approximation, which replaces the adjacency matrix by its realization average $\\langle A_{ij}\\rangle=k_i k_j/(2mN)$ and neglects fluctuations around the mean order parameter $M$; this reduces the interacting problem to independent spins in an effective local field $h_i^{\\rm eff}=h+Jm k_i$ and yields the central self-consistency equation for $M$. The second is the exact spin mapping $s_i'=2(s_i-\\tfrac12)$, which rewrites the 0/1 model as a classical $\\pm1$ Ising model with coupling $J/4$ and a degree-dependent local field $h_i'=h/2+(J/2)k_i$; the degree dependence is the trace of the asymmetry between the inactive and active states. Together these tools produce the critical-field estimate $h_c\\approx Jm$, explain the first-order jump and hysteresis seen in simulation, and provide a partition-function identity linking the modified model to established results on classical Ising models on scale-free networks.","core_discovery":"The paper's central claim is that the modified Ising Hamiltonian with spins $s_i\\in\\{0,1\\}$ and ferromagnetic coupling $J>0$, $$H_{[0,1]}=-\\frac12 \\sum_{i,j} J A_{ij}s_i s_j - h \\sum_i s_i,$$ placed on a Barabási-Albert network, undergoes a discontinuous first-order phase transition in the external field $h$ and shows hysteresis at low temperature. Using the annealed average $\\langle A_{ij}\\rangle = k_i k_j/(2mN)$ over network realizations, the mean-field reduction gives an effective local field $h_i^{\\rm eff}=h+Jm k_i$ and the self-consistency equation $M=\\frac1N\\sum_i [1+e^{-\\beta(h+JM k_i)}]^{-1}$, from which the authors derive the approximate critical field $h_c\\approx Jm$. The paper further shows that the 0/1 model maps exactly onto the classical $\\pm1$ Ising model on the same network through $s_i'=2(s_i-\\tfrac12)$, with renormalized coupling $J/4$, site-dependent field $h_i'=h/2+(J/2)k_i$, and a constant energy offset, so all classical Ising results transfer through $Z_{[0,1]}(J,h)=e^{\\beta E_0}Z_{[-1,1]}(J/4,h/2+Jk_i/2)$. Monte Carlo simulations on networks of $N=5\\times10^3$ nodes with $m=5$ confirm the discontinuous jump and the roughly rectangular hysteresis loop for ferromagnetic coupling.","pith_inferences":["The degree-dependent local field in the mapped classical model implies that, during a sweep, low-degree nodes should flip before high-degree hubs; this nucleation ordering is not stated in the paper but is directly testable in the same Monte Carlo setup.","Because the mean-field derivation anneals over graph realizations, the sharp $h_c=Jm$ prediction may fail on individual fixed networks; a practical test is to repeat the hysteresis protocol on several independent Barabási-Albert realizations and check whether the jump location changes between realizations.","The same effective-field construction should carry over to any network with a specified degree distribution, so the workflow could be applied to measured gene-regulatory connectivity matrices; the paper's own caveat that the connectivity matrix is assumed binary is the main obstacle for such an application."],"forward_implications":["For ferromagnetic coupling, sweeping the external field through $\\pm Jm$ at low temperature flips the whole network between the inactive ($M\\approx 0$) and active ($M\\approx 1$) states discontinuously, so a gene network can switch phenotype abruptly rather than gradually.","The critical field grows linearly with both the coupling $J$ and the preferential-attachment parameter $m$: stronger interactions or denser local wiring push the transition further away from zero field, making the active state harder to destroy.","Because the transition is first-order, the system is bistable near the critical field: the same value of $h$ can support two different activity levels depending on whether the field was increasing or decreasing (hysteresis).","Through the exact mapping, any quantity computed for the classical Ising model on a Barabási-Albert network (partition function, correlations, critical behavior) can be translated into the gene-type spin model with the renormalized coupling and degree-dependent field."],"supporting_citations":[{"why":"Defines the Barabási-Albert network model (growth and preferential attachment) whose scale-free degree distribution the simulations use.","marker":"Albert 2002"},{"why":"Supplies the annealed adjacency average $\\langle A_{ij}\\rangle=k_i k_j/(2mN)$ and the mean-field treatment of the Ising model on BA networks on which the analytic solution is built.","marker":"Bianconi 2002"},{"why":"Establishes the ferromagnetic phase transition of the classical Ising model on Barabási-Albert networks, the reference behavior the modified model is compared with.","marker":"Aleksiejuk 2002"},{"why":"Provides the single-spin-flip Monte Carlo algorithm used for thermal equilibration and hysteresis-loop sampling.","marker":"Metropolis 1953"}],"fun_headline_variants":["Gene-like 0/1 spins flip abruptly on scale-free networks","First-order jump in gene-state Ising on Barabasi-Albert nets","Modified Ising with 0/1 spins: discontinuous transition on scale-free networks","0/1 spins on scale-free nets: first-order jump with hysteresis"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central calculation assumes a single network can be represented by the annealed average $\\langle A_{ij}\\rangle=k_i k_j/(2mN)$ with negligible fluctuations about the mean spin, and it obtains $h_c=Jm$ by using the minimum degree $m$ rather than the mean degree $2m$; on a fixed finite network neither approximation is guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["Gene-like 0/1 spins flip abruptly on scale-free networks","First-order jump in gene-state Ising on Barabasi-Albert nets","Modified Ising with 0/1 spins: discontinuous transition on scale-free networks","0/1 spins on scale-free nets: first-order jump with hysteresis"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001589,"raw_usage":{"total_tokens":6453,"prompt_tokens":1182,"completion_tokens":5271,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":798,"completion_tokens_details":{"reasoning_tokens":5190}},"tokens_in":798,"tokens_out":5271,"duration_ms":34663,"temperature":1.0,"reasoning_tokens":5190,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:32:37.169601+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform the hysteresis protocol of Fig. 5 on several independent Barabási-Albert realizations with $N=5\\times10^3$, $m=5$, $J=1$, and record the field at which the order parameter jumps. If the jump is located at $h\\approx 2Jm=10$ instead of $h\\approx Jm=5$, or if its location varies noticeably between realizations of the same $(m,J)$, the annealed mean-field critical-field prediction is contradicted.","supporting_citations":[{"cited_title":"Statistical mechanics of complex networks Rev","cited_arxiv_id":null,"evidence_quote":"Defines the Barabási-Albert network model (growth and preferential attachment) whose scale-free degree distribution the simulations use."},{"cited_title":"Mean field solution of the Ising model on a Barab\\' a si–Albert network Physics Letters A 303 (2002) 166–168","cited_arxiv_id":null,"evidence_quote":"Supplies the annealed adjacency average $\\langle A_{ij}\\rangle=k_i k_j/(2mN)$ and the mean-field treatment of the Ising model on BA networks on which the analytic solution is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the ferromagnetic phase transition of the classical Ising model on Barabási-Albert networks, the reference behavior the modified model is compared with."},{"cited_title":"W, Rosenbluth, M","cited_arxiv_id":null,"evidence_quote":"Provides the single-spin-flip Monte Carlo algorithm used for thermal equilibration and hysteresis-loop sampling."}],"review_version":1}