{"id":"558bd3f8-33a3-4f8f-9cfe-18bb2c9a0d0b","arxiv_id":"1908.06704","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the critical 2D Ising model on a 2N by 2M rectangle with periodic-free boundaries, the total energy, after centering by (4/π)N ln N and scaling by sqrt(4MN ln N), converges to a Gaussian with variance 8/π as N and M grow together.","lead":"This paper proves that the total energy of the critical 2D Ising model on a large rectangle, with horizontal periodic and vertical free boundary conditions, has Gaussian fluctuations when both dimensions grow together. The main advance is a simultaneous-limit central limit theorem with an explicit boundary correction term, going beyond earlier iterated-limit results.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 9's key bound relies on g_θ ≥ 0, but (17) at β_c gives g_θ < 0 for small θ, so the proof of (29) is not established.","rationale":"I agree with the reader that the inequality gθ ≥ 0 is the weakest assumption, but the situation is more severe than 'plausible and probably true': it is contradicted by the paper's own formulas. Direct computation from (17) at β_c gives gθ < 0 for small θ, so the bound (72) in Lemma 9 is false as stated. Moreover, the paper's (74) and (54) are algebraically inconsistent with (17), indicating a sign or transcription error in a central formula. Since Lemma 9 is required to show that the f''_θ remainder in the Taylor expansion (62) vanishes, the proof of Proposition 1, and hence of Theorem 1, is not established. The additional gap from moment-generating-function convergence on t ≥ 0 to weak convergence (via the cited Problem 30.4) is secondary: that standard result requires convergence in a neighborhood of 0 or additional control of the lower tail, which is not supplied. Both issues could potentially be repaired, but as written the manuscript does not provide a valid proof of the central claim. Therefore the verdict should move from CONDITIONAL to REJECT (or at least UNVERDICTED until the false estimate is replaced).","tokens_in":12655,"tokens_out":13600,"duration_ms":121504,"concrete_test":"Evaluate gθ exactly at β = β_c from (17) for θ = π/(2N) with N = 10^2 and N = 10^4: compute cosh(2β_c) = √2, coth(2β_c) = √2, and sinh γθ = √((1−cosθ)(3−cosθ)); the resulting gθ is negative for small θ, contradicting the hypothesis of Lemma 9. Also compare the denominator in (54) with 1+η^{-4M}+(1−η^{-4M})gθ from (50) at β = β_c to verify that the value of gθ used in (54) is not the value given by (17); symbolically substituting θ = π/2 into both expressions shows they differ.","verdict_should_be":"REJECT","load_bearing_attack":"The decisive step is the proof of Lemma 9 (equation (71)), which is needed to control the second-order Taylor remainder of L4 in (62). Its proof begins with the assertion “Since g_θ ≥ 0 for each β ∈ (0, β_c]” and uses it to claim |1 + e^{-4Mγθ} + (1 − e^{-4Mγθ})gθ| ≥ 1 in (72). This assertion is false under the paper's own definition (17). At β = β_c, Lemma 2 gives coth(2β_c) = √2, cosh(2β_c) = √2, and sinh γθ = √((1−cosθ)(3−cosθ)), so gθ = √2 − √2 cosθ / √((1−cosθ)(3−cosθ)). For small θ, say θ = π/(2N), we have sinh γθ ~ θ and hence gθ ≈ √2 − √2/θ → −∞. Thus gθ is strictly negative for small θ, and (72) is not valid. A further internal inconsistency appears in (74), where |gθ| is written as cosh(2β)|[csch(2β) − cosθ] csch(γθ)|; this expression differs from (17) by an extra factor csch(γθ) in the first term. Also (54) implicitly uses gθ at β_c equal to √2(1−cosθ)^{1/2}(3−cosθ)^{−1/2}, which is not what (17) gives. Therefore Lemma 9 is not a minor missing justification; a load-bearing estimate in the proof is based on a false inequality. Without a replacement bound for the denominator in (51), the estimate (78) is unsupported, and Proposition 1(29) does not follow.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proves a central limit theorem for the total energy of the two-dimensional critical Ising model on a (2N) x (2M) rectangle, with periodic boundary condition in the horizontal direction and free boundary condition in the vertical direction. Theorem 1 states that, under the growth condition N(ln ln N)^2/(M ln N) -> 0, the normalized energy (E_{M,N} + 4√2 MN - (4/π)N ln N)/√(4MN ln N) has moment generating function converging to e^{4t^2/π} for each t ≥ 0, and hence converges weakly to a Gaussian with variance 8/π. The proof follows the classical partition-function-ratio method: the moment generating function is written as a ratio of partition functions at shifted inverse temperatures, and the exact McCoy-Wu formula (15) is expanded in Taylor series around β_c. The main work is in four estimates (26)-(29) for the components L_1,...,L_4 of the log partition function.","tokens_in":12966,"tokens_out":15492,"duration_ms":155660,"significance":"If the theorem is correct, it provides the first simultaneous-limit CLT for the critical Ising energy with a precise free-boundary contribution (4/π)N ln N, going beyond the earlier iterated limit of De Coninck. The derivation is coherent in outline and uses no fitted parameters: the constants 4/π and 8/π arise from the asymptotics of the exact partition function, and the paper connects the result to the expected non-existence of a usual probabilistic scaling limit for the energy field. These are genuine strengths. However, the proof as written contains a false inequality and an algebraic inconsistency in the two most delicate estimates, so the central claim is not currently established.","major_comments":[{"comment":"The assertion that g_θ ≥ 0 for every β ∈ (0, β_c] and every θ ∈ (0, π] is false. From (17) and Lemma 2, at β = β_c one has g_θ = √2 − √2 cos θ / √((1−cos θ)(3−cos θ)), which tends to −∞ as θ ↓ 0. Moreover, the same negativity occurs for β immediately below β_c and sufficiently small θ: at θ = 0, g_0 = coth(2β) − cosh(2β)/√(A(A−2)) with A = coth(2β) cosh(2β) > 2, and this tends to −∞ as β ↑ β_c. Since (72) is used to obtain |1 + e^{−4Mγ} + (1 − e^{−4Mγ})g_θ| ≥ 1, the bound (73) and the final estimate (78) are not established. Lemma 9, and therefore the proof of (29), rests on an invalid inequality and needs a replacement argument.","section":"Lemma 9, Eq. (72)"},{"comment":"Formula (54) for f'_θ at β_c is inconsistent with (17), (50), and (52). At β_c, (52) gives g'_θ = −2(1 + cos θ) csch(γ_θ), so (50) yields f'_θ|β_c = −2(1 − e^{−4Mγ_θ})(1 + cos θ) csch(γ_θ) / [1 + e^{−4Mγ_θ} + (1 − e^{−4Mγ_θ}) g_θ], with g_θ = √2 − √2 cos θ csch(γ_θ). The denominator written in (54) instead contains √2 (1 − cos θ)^{1/2} (3 − cos θ)^{-1/2}, which is not equal to g_θ (for small θ the former tends to 0 while g_θ tends to −∞). Therefore the decomposition in (64)-(69) and the proof of Lemma 8 do not establish the claimed limit (63). This is a load-bearing algebraic error, not a mere typo, because Lemma 8 supplies the free-boundary shift in (29).","section":"Eq. (54), Lemma 8"},{"comment":"The deduction of weak convergence from (10) is not justified as written. The paper cites Problem 30.4 of [2], but the standard result requires convergence of moment generating functions on an interval around 0 (or for all real t). Here (10) is proved only for t ≥ 0, and convergence of positive-t moment generating functions does not control left tails and does not by itself imply tightness or weak convergence. The author should either extend the moment generating function computation to t in a neighborhood of 0 or provide a separate tightness/characteristic-function argument.","section":"Theorem 1, proof after (10)"}],"minor_comments":[{"comment":"In the sentence defining the continuous function, the denominator is printed as 3 − 4 cos θ − cos² θ; it should be 3 − 4 cos θ + cos² θ, as in (47) and Lemma 2.","section":"Proof of (28), around Eq. (47)"},{"comment":"The last term in (78), C21 M² csch²(γ_θ)/(M N ln N), is more clearly written as C21 M csch²(γ_θ)/(N ln N); the displayed form is confusing though mathematically equivalent.","section":"Eq. (78)"},{"comment":"The proof of (33) says it follows from the monotonicity of csch and the mean value theorem; for clarity, the reader needs the explicit bound on |β − β_c| in terms of 1/√(4M N ln N), which is implicit and could be stated.","section":"Lemma 4, inequality (33)"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a natural and worthwhile problem, and the partition-function-ratio strategy is plausible. However, the proof of Lemma 9 relies on an inequality that is false in exactly the regime used, and Lemma 8 contains an algebraic error in its central formula. These are not presentation issues; they invalidate the proof as it stands. I still recommend major revision rather than rejection because the overall method may be repairable, but the revision must supply a correct denominator estimate for the second Taylor remainder and a correct derivation of the f'_θ asymptotics, and it must address the t ≥ 0 versus weak-convergence gap."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jianping Jiang proves a simultaneous-limit CLT for the total energy of the critical 2D Ising model on a 2N×2M rectangle, with the free-boundary N ln N term made explicit. The math looks right to me, and the stress-test's main objection does not survive reading the actual formula. Equation (17) is plainly meant to be (coth(2β) − cosh(2β) cosθ)/sinh γθ, not coth − (cosh cos)/sinh. With the correct reading, at β_c one gets gθ = √2(1−cosθ)/√((1−cosθ)(3−cosθ)) ≥ 0, and since the numerator decreases in β, gθ ≥ 0 on (0, β_c]. So the 'gθ < 0 for small θ' example is an artifact of misparsing the fraction. The same correct reading makes (54) and (74) consistent.\n\nNow the substance. The new content: De Coninck had the iterated limit M→∞, N→∞; here the limits are simultaneous under M ≥ N/(ln N)^α (α<1), and the expected energy has the boundary term (4/π)N ln N subtracted. That term is derived from the partition function, not inserted. The method is a clean use of the McCoy–Wu product formula: Taylor expand the three parts L1, L3, L4, control sums with Lemma 7's exponential-in-M estimates and harmonic bounds. I did not find fitted parameters or circular reasoning.\n\nSoft spots, in proportion: The passage from (10) to weak convergence is one sentence, citing Problem 30.4 of Billingsley. That is terse; the argument requires continuity of the limiting MGF at 0 to rule out mass escaping to −∞. I believe the cited problem covers it, but a referee will want the one-line tightness explanation. Lemma 4's proof is compressed and the algebra behind (34)–(35) is easy to stumble over, but it is correct. Nothing else looks broken.\n\nThe only real presentation issue: the paper only proves the MGF for t≥0, and Remark 5 acknowledges that t<0 needs more work. For a CLT, one-sided convergence plus continuity at 0 is enough, but the authors should say so explicitly.\n\nWho is it for: people working on critical Ising fluctuations, energy field scaling limits, and numerical studies of finite-size corrections. It deserves a serious referee and probably acceptance after minor revisions.","headline":"The stress-test's 'false inequality' is a misreading of formula (17); with the intended fraction, gθ ≥ 0 on (0,β_c] and the paper's CLT argument holds together.","tokens_in":13529,"tokens_out":16682,"would_cite":true,"duration_ms":160855,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60F05","82B20","82B27"],"pacs":[],"model":"deepseek-v4-flash","headline":"The total energy of the critical two-dimensional Ising model on a growing rectangle converges to a Gaussian with mean zero and variance 8/π.","keywords":["Ising model","critical temperature","central limit theorem","total energy","partition function","free boundary condition","periodic boundary condition","two-dimensional lattice"],"falsifier":"Compute $g_\\theta=(\\coth(2\\beta)-\\cosh(2\\beta)\\cos\\theta)/\\sinh\\gamma_\\theta$ on a fine grid of $\\theta\\in(0,\\pi]$ for $\\beta$ just below $\\beta_c$; a negative value would disprove the unproved inequality on which Lemma 9 rests. Alternatively, simulate the critical Ising model on $2N\\times2N$ periodic-free rectangles for increasing $N$ and compare the empirical distribution of $(E_{M,N}+4\\sqrt{2}MN-(4/\\pi)N\\ln N)/\\sqrt{4MN\\ln N}$ with a standard Gaussian; a systematic departure in variance from $8/\\pi$ would refute the theorem.","tokens_in":12397,"feed_emoji":"🧲","tokens_out":6487,"duration_ms":64920,"temperature":0.7,"pith_summary":"This paper proves a central limit theorem for the total energy (Hamiltonian) of the critical two-dimensional Ising model on a $[1,2N]\\times[1,2M]$ rectangle, with periodic boundary condition horizontally and free boundary condition vertically, when both $M$ and $N$ tend to infinity simultaneously. It provides a simultaneous-limit CLT for this observable and isolates the free-boundary contribution $(4/\\pi)N\\ln N$ to the mean energy. After centering by $4\\sqrt{2}MN-(4/\\pi)N\\ln N$ and scaling by $\\sqrt{4MN\\ln N}$, the energy has moment generating function tending to $e^{4t^2/\\pi}$, hence converges weakly to a Gaussian with mean $0$ and variance $8/\\pi$. The result matters because it shows precisely how boundary conditions shape the energy fluctuations at criticality and because the same mechanism is a step toward understanding the scaling limit of the energy field.","feed_headline":"Critical 2D Ising energy is Gaussian with variance 8/π","feed_subtitle":"Simultaneous growth in both directions yields the free-boundary term (4/π)N ln N and a standard normal limit.","key_machinery":"The central object is the exact partition function of the periodic-free rectangle, written as a product over Fourier modes $\\theta=\\pi(2n-1)/(2N)$, with $\\cosh\\gamma_\\theta=\\coth(2\\beta)\\cosh(2\\beta)-\\cos\\theta$ and $\\gamma_\\theta\\ge 0$. The moment generating function is the ratio $Z_{M,N}(\\beta-s)/Z_{M,N}(\\beta)$, and the proof expands $\\ln Z$ around the critical inverse temperature $\\beta_c$. The essential mechanism is that the mode sum for $L_3$ produces the variance $8/\\pi$ through a sum of $16(3-4\\cos\\theta+\\cos^2\\theta)^{-1/2}$ behaving like $4N\\ln N$, while $L_4$ contributes the boundary term $-(4/\\pi)N\\ln N$ to the centering.","core_discovery":"Theorem 1 asserts that for each $t\\ge 0$, $\\lim_{N\\to\\infty}\\langle e^{t\\hat E_{M,N}}\\rangle_{\\beta_c}=e^{4t^2/\\pi}$, where $\\hat E_{M,N}=(E_{M,N}+4\\sqrt{2}MN-(4/\\pi)N\\ln N)/\\sqrt{4MN\\ln N}$, provided $N(\\ln\\ln N)^2/(M\\ln N)\\to 0$. Consequently $\\hat E_{M,N}$ converges weakly to a Gaussian with mean $0$ and variance $8/\\pi$. The centering contains both the bulk term $4\\sqrt{2}MN$ and the free-boundary term $(4/\\pi)N\\ln N$; the proof obtains the Gaussian limit by expanding the logarithm of the moment generating function into four parts $L_1,L_2,L_3,L_4$ coming from the log partition function.","pith_inferences":["If the denominator bound used in Lemma 9 can be extended to both signs of $t$, the same Taylor-expansion scheme would likely yield the CLT for all real $t$ and, as the paper conjectures, for other boundary conditions such as free or all-plus.","A natural testable extension is to subtract the boundary term and study the distributional scaling limit of the energy field after renormalization; the paper only conjectures a white-noise limit, so a precise theorem would require new mode-sum estimates.","The ratio-of-partition-functions method may transfer to other exactly solvable lattice models whose partition function factors over Fourier modes, converting a CLT for additive observables into a Taylor-expansion estimate for mode sums.","The unproved sign condition on $g_\\theta$ is checkable by direct numerical evaluation; if it fails, the proof has a gap, although the theorem itself might still be true by a sharper two-sided estimate."],"forward_implications":["The same CLT holds whenever $M\\ge N/(\\ln N)^\\alpha$ for $\\alpha\\in[0,1)$, so the result covers squares and moderately rectangular boxes; in particular $M=N$ is allowed.","The variance grows as $(32/\\pi)MN\\ln N$, a factor $\\ln N$ larger than the $MN$ bulk scale, showing that the free vertical boundaries dominate the energy fluctuations.","The proof identifies the exact leading free-boundary contribution $(4/\\pi)N\\ln N$ to the mean energy, consistent with conformal-covariant behavior of the energy density near the boundary.","The Gaussian limiting fluctuation, together with the $|z_1-z_2|^{-2}$ decay of energy correlations, suggests that the critical energy field does not have a finite white-noise scaling limit in the usual probabilistic sense.","The theorem gives the first simultaneous-limit statement of this type for the critical Ising energy, removing the need to send $M\\to\\infty$ before $N\\to\\infty$ as in earlier work."],"supporting_citations":[{"why":"Supplies the exact partition function formula for the periodic-free rectangle that is the starting point of the moment-generating-function computation.","marker":"[11]"},{"why":"Establishes the earlier iterated-limit CLT for the same model that the paper extends to simultaneous limits.","marker":"[4]"},{"why":"Provides the ratio-of-partition-functions and Taylor-expansion strategy adapted in the proof.","marker":"[5]"},{"why":"Gives the standard probability theorem converting convergence of moment generating functions into weak convergence.","marker":"[2]"}],"fun_headline_variants":["2D critical Ising energy has Gaussian limit with variance 8/π","Critical Ising energy fluctuations are Gaussian: variance 8/π","Ising energy at criticality: Gaussian limit with variance 8/π","CLT for 2D critical Ising energy: variance 8/π"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof assumes without proof that the function $g_\\theta=(\\coth(2\\beta)-\\cosh(2\\beta)\\cos\\theta)/\\sinh\\gamma_\\theta$ is nonnegative for every $\\beta\\in(0,\\beta_c]$ and every $\\theta\\in(0,\\pi]$, because this keeps the denominator of the second-derivative estimate at least one in absolute value.","fun_headline_variants_meta":{"raw":{"variants":["2D critical Ising energy has Gaussian limit with variance 8/π","Critical Ising energy fluctuations are Gaussian: variance 8/π","Ising energy at criticality: Gaussian limit with variance 8/π","CLT for 2D critical Ising energy: variance 8/π"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000763,"raw_usage":{"total_tokens":3352,"prompt_tokens":881,"completion_tokens":2471,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":2389}},"tokens_in":497,"tokens_out":2471,"duration_ms":16334,"temperature":1.0,"reasoning_tokens":2389,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:39:57.617428+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $g_\\theta=(\\coth(2\\beta)-\\cosh(2\\beta)\\cos\\theta)/\\sinh\\gamma_\\theta$ on a fine grid of $\\theta\\in(0,\\pi]$ for $\\beta$ just below $\\beta_c$; a negative value would disprove the unproved inequality on which Lemma 9 rests. Alternatively, simulate the critical Ising model on $2N\\times2N$ periodic-free rectangles for increasing $N$ and compare the empirical distribution of $(E_{M,N}+4\\sqrt{2}MN-(4/\\pi)N\\ln N)/\\sqrt{4MN\\ln N}$ with a standard Gaussian; a systematic departure in variance from $8/\\pi$ would refute the theorem.","supporting_citations":[{"cited_title":"Mccoy and T.T","cited_arxiv_id":null,"evidence_quote":"Supplies the exact partition function formula for the periodic-free rectangle that is the starting point of the moment-generating-function computation."},{"cited_title":"De Coninck (1984)","cited_arxiv_id":null,"evidence_quote":"Establishes the earlier iterated-limit CLT for the same model that the paper extends to simultaneous limits."},{"cited_title":"De Coninck (1987)","cited_arxiv_id":null,"evidence_quote":"Provides the ratio-of-partition-functions and Taylor-expansion strategy adapted in the proof."},{"cited_title":"Billingsley (1995)","cited_arxiv_id":null,"evidence_quote":"Gives the standard probability theorem converting convergence of moment generating functions into weak convergence."}],"review_version":1}