{"id":"0f656bc2-28c3-46f2-91be-75eea5854f13","arxiv_id":"2607.16735","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Closed-form elliptic-integral approximations for the field-dependent thermodynamics of the 2D Ising model, exact at H=0, valid at high temperature and weak field.","lead":"This paper derives approximate closed-form formulas for the free energy, magnetization, susceptibility, internal energy, and specific heat of the two-dimensional Ising magnet in an external magnetic field. The formulas are exact when the field is zero, match Monte Carlo data at high temperature and weak field, but miss the physics at and below the critical point.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"ψ→4ψ truncation drops a term whose H-derivative is nonzero at H=0; Eq. (60) has the wrong O(βJ) coefficient, so the small-H response is uncontrolled.","rationale":"The paper's central claim is that derivatives of Eq. (51) give field-dependent thermodynamics matching BKL in the small-angle regime. The crux is how the H-dependence of the truncated transfer matrix is handled. The exact T contains a field term tanφ(Z_k+Z_{k+1}) in Eq. (32) whose derivative at H=0 is first order in βe^{2βJ}; by setting tanφ=0 before differentiating, the approximation removes the leading contribution to χ_0. This is not merely an absent error bound: it is a demonstrable error in the first nontrivial derivative. My expansion of Eq. (60) shows the O(βJ) coefficient is −3, while the exact high-temperature series gives +4. The paper's Fig. 3 at J=1, T/J=6 should show this discrepancy; reporting absolute differences hides it because β is small. I still do not think this warrants REJECT: the paper is candid about limitations, reduces to Onsager at H=0, and reproduces BKL at very high temperature/weak field. The reader's CONDITIONAL verdict is the right one, but the condition should be sharpened: the author should supply a next-order error estimate and demonstrate relative, not just absolute, agreement. This is the same load-bearing weakness the reader identified, sharpened to a concrete, testable failure.","tokens_in":18577,"tokens_out":24355,"duration_ms":210851,"concrete_test":"Analytically expand Eq. (60) in v=tanh βJ and compare the O(v) coefficient with the exact high-temperature susceptibility expansion βχ = 1 + 4v + 12v² + ... for the square-lattice Ising model (or Eq. (65) of ref [24] at small βJ). Then evaluate both at a point inside the claimed validity window, e.g. J=1, T=6 (βJ=1/6): the exact βχ ≈ 1 + 4v ≈ 1.66, while Eq. (60) gives βχ_a ≈ 1 − 3βJ ≈ 0.5, a large relative discrepancy. If the O(v) coefficient is not +4, the H-derivative of the truncated free energy is not controlled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the exact transfer matrix, the field enters Eq. (32) through tanφ (Z_k+Z_{k+1}). The small-angle step in Sec. III.B sets tanφ=0 and X'_k=X_k, then rescales ψ→4ψ. The discarded term vanishes at H=0, but its H-derivative at H=0 is nonzero: ∂tanφ/∂H|_0 = (β/2)e^{2βJ}. It therefore contributes directly to ∂²(βF)/∂H²|_0, i.e. to the zero-field susceptibility. Dropping it before differentiating is not a controlled truncation. A concrete symptom is Eq. (60): expanding around βJ=0 gives χ_a/β = 1 − 3βJ + O((βJ)^3), whereas the exact square-lattice high-temperature series is βχ = 1 + 4v + O(v²) with v=tanh βJ. The approximate formula misses the leading interaction correction and even reverses its sign; it approaches the exact χ only at infinite temperature. Consequently, in the H→0 limit at fixed βJ, the approximate free energy's H² coefficient does not converge to the exact one, even though βF_a at H=0 is exactly Onsager. The BKL agreement at J=0.01 reflects that the absolute error O(βJ) is numerically small, but the error is relative and the stated validity condition φa≪1 does not restore the discarded term. The same uncontrolled derivative underlies the spurious κ'=0 divergence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes approximate closed-form thermodynamics for the 2D Ising model in a magnetic field. Using the quaternion/Kaufman transfer-matrix representation, the author applies a 'small-angle' truncation (Sec. III.B) and a field rescale psi -> 4 psi; the resulting centrosymmetric transfer matrix is then diagonalized by Kaufman's method, giving the approximate free energy beta F_a in Eq. (51). Derivatives yield approximate expressions for M, chi, U, and C (Eqs. 52, 55, 57, 58), which reduce to the Onsager results at H=0 and are compared with BKL Monte Carlo data. A supplementary section applies the same construction to the 1D Ising model. The paper states the approximation is valid when phi_a << 1 and acknowledges failures near criticality.","tokens_in":18983,"tokens_out":23527,"duration_ms":224159,"significance":"If the derivation were controlled, this would be a useful contribution: explicit field-dependent formulas in a regime without an exact solution, with an exact H=0 limit, analytic derivatives, and direct comparison to simulations in a high-temperature weak-field window. The algebraic path from Eqs. (13)-(45) is coherent, and the Monte Carlo comparisons at J=0.01 show that the formulas capture the magnitudes of M, chi, U, and C in the tested regime. However, as detailed below, the small-angle truncation is uncontrolled, the field rescale is calibrated rather than derived, and the 1D companion section contains an algebraic error. The paper is therefore better viewed as a promising phenomenological approximation than as a derivation of the model's thermodynamics; in its current form it needs substantial revision.","major_comments":[{"comment":"The small-angle step is not a controlled truncation. Setting tan phi = 0 in Eq. (32) removes a term whose H-derivative at H=0 is nonzero; rescaling psi -> 4 psi in Eq. (36) cannot restore it. Consequently, the zero-field susceptibility Eq. (60) expands to beta chi_a = 1 + beta J + O((beta J)^3), while the exact square-lattice high-temperature series is beta chi = 1 + 4 v + ... = 1 + 4 beta J + ... . The leading interaction correction is missed by a factor of four. Moreover, phi_a = 0 identically at H=0 for any beta J, so the stated criterion phi_a << 1 does not imply beta J << 1; yet Eq. (60) is not exact at finite beta J. No error bound or next-order term is provided, so the field derivatives of Eq. (51) are not controlled in the claimed validity region.","section":"Sec. III.B, Eqs. (32), (36), (60)"},{"comment":"The factor 4 in psi -> 4 psi (and 2 in 1D) is inserted after the truncation, and Sec. VI states it was 'necessary for the approximate expressions to agree with the numerical simulations.' This makes the factor a fitted parameter rather than a consequence of the approximation. The same numerics are also used to define the validity region, so agreement cannot independently validate the central claim. The paper should either derive this factor from a systematic expansion, or be reframed explicitly as a calibrated/phenomenological approximation.","section":"Sec. III.B, Eq. (36); Sec. VI"},{"comment":"Eq. (C35) is not the derivative of Eq. (C34). Differentiation gives chi = 2 beta e^(2 beta J) [cosh(2 beta H)+e^(2 beta J)] / [1+e^(2 beta J) cosh(2 beta H)]^2, not the displayed expression. At J=0 the correct derivative is 2 beta / [1+cosh(2 beta H)], whereas Eq. (C35) gives 2 beta / [1+cosh(2 beta H)]^3; the exact value is beta sech^2(beta H). The sentence after Eq. (C37) claiming exact J=0 agreement is therefore false. This error in the companion validation should be corrected.","section":"Supplementary Material, Eqs. (C34)-(C37)"}],"minor_comments":[{"comment":"The Monte Carlo comparisons give absolute-difference thresholds but no statistical errors or number of independent runs; finite-size effects between n=20 and n=50 are not quantified. Please add error bars or state uncertainties.","section":"Secs. IV-V"},{"comment":"The kappa' = 0 divergence is described as a 'sharp transition between ordered and disordered phase,' but for H != 0 the exact model has no phase transition; the text should call this an artifact of the approximation.","section":"Sec. VI"},{"comment":"'phi_a small implies 2 beta J << 1 and 2 beta H << 1' is not a valid implication when H=0 (phi_a=0 for any beta J); the condition T >> 2(H+J) should be justified separately.","section":"Sec. III.B, Eq. (36)"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the manuscript is a borderline case. The central derivation is uncontrolled and the rescale is calibrated, not derived; the supplementary 1D susceptibility contains an algebraic error. I would require a major revision: either supply an error estimate for the neglected O(phi) terms and a derivation of the rescale, or explicitly rewrite the paper as a fitted approximation with restricted claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives closed-form elliptic-integral expressions for the free energy and its derivatives for the 2D Ising model in a field, and the H=0 limit reproduces Onsager exactly. That is real work, and the algebra through the quaternion transfer matrix is internally consistent. The BKL comparisons in the high-temperature, small-φa regime look fine, and the author is candid about the regime's limits. This is a legitimate attempt at a hard problem, not a crank derivation.\n\nThe soft spot is the one the stress-test note lands on: the small-angle truncation drops the tanφ term before differentiating, and that term's H-derivative at H=0 is nonzero. Concretely, Eq. (60) expands to χ_a/β = 1 − 3βJ + O((βJ)²), whereas the exact high-temperature series is βχ = 1 + 4v with v=tanh(βJ). So the approximate susceptibility at H=0 misses the leading interaction correction and even gets the sign wrong. The MC agreement at J=0.01 is just because the absolute error is small there; the error is relative, and the stated validity condition φa≪1 does not restore the discarded term. This means the field dependence of the free energy is not actually controlled by the approximation, only by the rescale ψ→4ψ, which the author admits was chosen to match numerics.\n\nOther weaknesses are proportionate: the κ′=0 divergence at H≠0 is spurious, and the near-critical susceptibility is explicitly out of reach, as the author says. The paper is honest about these. But the core issue is that the approximation's error is not estimated, and the one clean check at H=0 fails for the susceptibility, which is the quantity that matters.\n\nWho gets value from this? Someone working on approximate transfer-matrix methods for the Ising model in a field, or wanting a worked example of the Kaufman-Onsager machinery. Not someone needing reliable field-dependent thermodynamics near criticality. It deserves a serious referee, because the derivation is concrete and the claims are falsifiable, but the referee should ask for an error bound or next-order term, and a proper expansion of χ at H=0 that recovers the exact first-order coefficient.\n\nRecommendation: send it to peer review, but expect major revision before it's trustworthy as more than a curve-fit.","headline":"A genuine but uncontrolled approximation: the field-dependent closed forms are new and the H=0 limit is exactly Onsager, but the ψ→4ψ rescale is fit to numerics and the zero-field susceptibility is wrong at first order in βJ.","tokens_in":19486,"tokens_out":2044,"would_cite":false,"duration_ms":24003,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B20","82B26","82B80"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that the two-dimensional Ising model in an external magnetic field, at high temperature and weak field, has its free energy given by a closed expression built from complete elliptic integrals and hyperbolic functions, with","keywords":["two-dimensional Ising model","external magnetic field","quaternion representation","small-angle approximation","transfer matrix","free energy","magnetization","susceptibility"],"falsifier":"Compute the exact largest eigenvalue of the full 2D transfer matrix for a finite lattice (e.g., N = 20) at values of βJ and βH where φa is small but not tiny (say φa = 0.1), and compare it to the largest eigenvalue of T_a; if the relative difference does not vanish at a rate consistent with the claimed O(φ) truncation, the approximation's validity region is not controlled in the way the paper states.","tokens_in":18391,"feed_emoji":"🧲","tokens_out":3518,"duration_ms":36614,"temperature":0.7,"pith_summary":"The two-dimensional Ising model in a magnetic field has no exact solution, and the symmetry-broken transfer matrix resists the methods used for the zero-field case. This paper proposes that a small-angle rotation in the quaternion representation of the model produces a centrosymmetric transfer matrix whose largest eigenvalue can be written down explicitly, yielding a closed-form approximate free energy. From that free energy the paper derives analytic expressions for magnetization, susceptibility, internal energy, and specific heat that depend on the field strength. The paper reports that these expressions match BKL simulations in the regime T ≫ 2(H + J) and reduce exactly to the Onsager results when the field is zero. The importance lies in providing the first field-dependent closed-form thermodynamics for the model in any region, however restricted.","feed_headline":"Closed-form thermodynamics found for 2D Ising model in a field","feed_subtitle":"Small-angle quaternion approximation yields analytic free energy and derived quantities at high temperature, matching simulations in the val","key_machinery":"The quaternion (Pauli-matrix) representation of the transfer matrix, in which the field-dependent single-spin operators are rotated by an angle φ relative to the coupling operators. The approximation takes φ → 0, dropping the tan(φ) coupling terms while rescaling ψ → 4ψ, so the transfer matrix becomes centrosymmetric and can be diagonalized with Kaufman's plane-rotation machinery. The resulting largest eigenvalue depends on the field through the effective parameters η and ωa, yielding the free energy and all derived thermodynamic quantities.","core_discovery":"The central claim is that, for φa = arctan[sinh(2βH) exp(2βJ)] ≪ 1, the free energy of the 2D Ising model in a magnetic field is βFa = −log[2 cosh(2η) cosh(2ψ)] − (1/π)∫₀^{π/2} dθ log[½(1 + √(1 − κ² sin²θ))], where tanh η = [cosh(2ψ)e^{2βJ} − 1]/[cosh(2ψ)e^{2βJ} + 1], κ = 2 sinh(2η)/cosh²(2η), and ψ = βH. The paper derives the magnetization, susceptibility, internal energy, and specific heat as derivatives of this expression, correcting the magnetization by a factor of 1/4 to account for the small-angle rescaling ψ→4ψ. It claims these formulas agree with BKL simulations in the stated regime and reduce exactly to the known H = 0 results.","pith_inferences":["The validity region is empirical: no error bound is derived for the φ → 0 truncation, so the claimed regime is defined by observed simulation agreement rather than by a controlled expansion, and the method is unlikely to be reliable near the critical temperature or at moderate fields.","The empirical factor of 4 in the ψ rescaling (and 2 in one dimension) suggests a coordination-number-dependent normalization; testing whether φa = arctan[sinh(2βH) exp(2βJ)] remains the correct angle for lattices with coordination numbers other than 4 would be a natural next step.","The closest existing benchmark for the susceptibility, the high-temperature perturbative series, is compared in the paper and found to be more accurate in some ranges; a head-to-head test of the two approximations against high-precision simulations at φa just below the stated validity threshold could separate truncation error from field-rescaling error.","The machinery might be extended to produce approximate spin-spin correlation functions, since the plane-rotation diagonalization used for the eigenvalue also yields eigenvectors, though the paper does not attempt this."],"forward_implications":["If correct, the four thermodynamic functions (magnetization, susceptibility, internal energy, specific heat) are available in closed form, explicitly dependent on the external field, in the high-temperature weak-field regime.","The expressions reduce exactly to the Onsager free energy, internal energy, and specific heat at H = 0, providing a smooth interpolation between the known zero-field exact results and the field-dependent approximation.","The approximation defines a precise validity criterion, φa ≪ 1, equivalently T ≫ 2(H + J), which the paper checks against BKL simulations for J = H and J ≠ H.","The mechanism, applied to the one-dimensional model in the supplementary material, yields closed-form approximations that also reduce to the exact paramagnetic limit when J = 0.","The paper suggests the same construction can be extended to anisotropic couplings, other two-dimensional lattice types, and—if a quaternion representation exists—possibly the three-dimensional model."],"fun_headline_variants":["Analytic Ising thermodynamics found for 2D in a field","Small-angle quaternion gives 2D Ising free energy in field","Ising in a magnetic field: approximate analytic formulas","Approximate thermodynamics for 2D Ising in external field","Closed-form Ising thermodynamics via quaternion small-angle trick"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the small-angle truncation (φ → 0, tan φ = 0) together with the rescaling ψ → 4ψ leaves the largest eigenvalue of the exact field-dependent transfer matrix essentially unchanged, an assumption supported only by numerical agreement and never by a derived error bound.","fun_headline_variants_meta":{"raw":{"variants":["Analytic Ising thermodynamics found for 2D in a field","Small-angle quaternion gives 2D Ising free energy in field","Ising in a magnetic field: approximate analytic formulas","Approximate thermodynamics for 2D Ising in external field","Closed-form Ising thermodynamics via quaternion small-angle trick"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000967,"raw_usage":{"total_tokens":3949,"prompt_tokens":739,"completion_tokens":3210,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":483,"completion_tokens_details":{"reasoning_tokens":3124}},"tokens_in":483,"tokens_out":3210,"duration_ms":22928,"temperature":1.0,"reasoning_tokens":3124,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T20:07:31.572285+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact largest eigenvalue of the full 2D transfer matrix for a finite lattice (e.g., N = 20) at values of βJ and βH where φa is small but not tiny (say φa = 0.1), and compare it to the largest eigenvalue of T_a; if the relative difference does not vanish at a rate consistent with the claimed O(φ) truncation, the approximation's validity region is not controlled in the way the paper states.","supporting_citations":[],"review_version":1}