{"id":"ff896817-c227-4382-953c-8330e9d9cd3a","arxiv_id":"2506.19407","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A single integrated pair correlation function, measured across interaction strengths, suffices to reconstruct pressure, entropy, heat capacity, and other thermodynamic variables of a quantum many-body system.","lead":"Physicists show that all thermodynamic quantities, such as pressure, entropy, and heat capacity, can be extracted from measurements of a single two-particle correlation function. The method turns the standard Maxwell relations of thermodynamics into a practical recipe for quantum simulators, where correlations are easy to measure but thermodynamics is not.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Experimental universality requires c to enter the Hamiltonian exactly as c G2; real control knobs that shift H0 or the G2 operator make the measured G2 a biased conjugate variable, and the paper does not quantify this failure mode.","rationale":"The mathematical core of the paper is sound: under the explicit assumption H(c)=H0+cG2 with H0 and G2 c-independent, the Hellmann–Feynman theorem gives ∂F/∂c=G2, and the generalized Maxwell relations follow from equality of mixed partial derivatives. The benchmarks against Bethe ansatz and DMRG support the identity, and the paper contains no free parameters or circular fitting. The load-bearing weakness is not in the derivation but in the breadth of the 'universal/readily accessible' claim: the construction requires that the experimentally tuned parameter c appear exactly linearly in the Hamiltonian. The paper does not address what happens when the tuning mechanism also changes H0 or the operator G2, which is a real possibility in several quantum simulator platforms. The reader's weakest_assumption identifies precisely this issue, and I agree. Because the central theoretical result is correct and the experimental overclaim is a condition that could be met by careful engineering, the appropriate verdict remains CONDITIONAL; no change from the reader's verdict is needed.","tokens_in":15117,"tokens_out":6531,"duration_ms":76181,"concrete_test":"Take the transverse-field Ising Hamiltonian with a coupled control: H(c) = (h0 + λ c) ∑_j S^x_j + c ∑_j S^z_j S^z_{j+1}. For several λ values (including λ=0 and a realistically nonzero λ), compute the exact finite-temperature magnetization or entropy from the full Hamiltonian (Jordan-Wigner or exact diagonalization), and compare it with the value inferred from Eq. (6) using the bare G2(c) = ⟨∑ S^z_j S^z_{j+1}⟩ but ignoring the λ term (i.e., the procedure the paper prescribes). If the inferred value deviates from the exact value by an amount that grows with λ, the experimental universality claim is not supported; if the deviation is negligible for λ within the range of current ion-trap experiments, the concern is mitigated. This directly tests the correction terms in ∂F/∂c=G2+⟨∂H0/∂c⟩.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (17) — ∂F/∂c = G2 — is derived from the Hellmann–Feynman theorem under the assumption H(c) = H0 + c G2 with H0 and G2 independent of c. If the experimental control called 'c' also modifies H0 or the operator G2, then ∂F/∂c = G2 + ⟨∂H0/∂c⟩ + c ⟨∂G2/∂c⟩ (operator derivative). Equations (5) and (7) then acquire extra terms such as ∂/∂X ⟨∂H0/∂c⟩ + ∂/∂X [c ⟨∂G2/∂c⟩], and the integral relations (6) and (8) are biased by integrals of these corrections. The paper's benchmarks keep H0 and G2 strictly c-independent (e.g., Eq. (26) with h_x held fixed), so they do not test this failure mode. Yet the abstract claims the approach is 'readily accessible in quantum simulators'; in several leading platforms the tuning knob does not satisfy the exact linear-coupling condition. Trapped-ion Ising simulators, for example, generate the spin–spin coupling with laser beams that also cause AC-Stark shifts (entering H0), and lattice-depth variation in optical lattices changes both t and U as well as the Wannier functions that define the G2 operator. The authors do not provide a parametric analysis of the resulting error, so the universal experimental claim rests on an unvalidated exactness condition.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript derives generalized Maxwell relations for Hamiltonians of the form H = H0 + c G2, where G2 is an integrated pair-correlation operator. The central identities are Eq. (5), relating the derivative of a thermodynamic quantity Y with respect to the interaction strength c to a derivative of G2 with respect to a standard thermodynamic variable X, and Eq. (7) for second derivatives, together with the integral forms Eqs. (6) and (8). The authors argue that this lets one reconstruct pressure, entropy, chemical potential, magnetization, heat capacity, and compressibility from measurements of G2, given a known reference value. They illustrate the method on the Lieb-Liniger model, the Yang-Gaudin model, the Fermi- and Bose-Hubbard models, and the transverse-field Ising model, benchmarking against exact Bethe-ansatz results, DMRG calculations, and the exact Ising solution. They also derive an approximate analytic entropy for the high-temperature Yang-Gaudin gas and show that ∂G2/∂hz displays a lambda anomaly near the Ising quantum critical point.","tokens_in":15241,"tokens_out":8863,"duration_ms":106239,"significance":"The theoretical content is elementary but broadly useful: the generalized Maxwell relations are exact identities with no fitted parameters, and the numerical benchmarks are extensive and convincing, including strong-interaction regimes and a quantum critical point. If the experimental claims are appropriately qualified, the method could provide a practical bridge between correlation-function measurements and thermodynamics in quantum simulators. The paper also gives credit to prior work by citing the existing Maxwell relation between entropy and atom-atom pair correlation and extends it to a full family of thermodynamic quantities. The main weakness is that the experimental universality claim rests on a linear-coupling assumption that is not satisfied by every control knob in current platforms, and this failure mode is not quantified.","major_comments":[{"comment":"The central reconstruction formula, Eq. (17), relies on the exact linear-coupling assumption H(c) = H0 + c G2 with H0 and G2 independent of c. The paper's claim that the approach is 'readily accessible in quantum simulators' (abstract and Conclusions) is not supported for common tuning schemes where this condition fails: Feshbach sweeps alter the effective range of the interaction, optical-lattice depth changes both t and the Wannier functions defining the onsite operators, and trapped-ion laser-induced spin couplings generically produce AC-Stark shifts in H0. In such cases Eq. (17) acquires extra terms involving ∂H0/∂c and ∂G2/∂c, and Eqs. (5)-(8) receive corrections that the paper neither derives nor bounds. The benchmarks, e.g., Eq. (26) for the transverse-field Ising model with hx held fixed, deliberately keep H0 and G2 c-independent and therefore do not test this failure mode. Please either restrict the experimental universality claim to control knobs that realize the exact linear coupling, or provide a quantitative error estimate for representative experimental platforms.","section":"Methods, Maxwell relations; Conclusions"},{"comment":"The heat-capacity benchmark in Fig. 2f is obtained after employing a Savitzky-Golay filter to smooth G2 at low temperatures, but the manuscript does not report the filter window, polynomial order, or any comparison of the smoothed second derivative with the raw numerical second derivative. Since CV is a second derivative of G2 (Eq. (13)), the filter can introduce a systematic bias that is not visible in the figure. Please provide the filter parameters and a convergence or robustness study, for example by varying the window width, before the agreement with the energy-variance calculation can be taken as quantitative.","section":"Methods, Transverse field Ising model"},{"comment":"With m defined as -V^{-1}∇_h F in the main text and G2 defined as the expectation value of the extensive operator in Eq. (3), the Maxwell relation should read ∂m/∂c = -V^{-1} ∂G2/∂h, not -∂G2/∂h. The derivation in Supplementary Information S6 appears to use the intensive nearest-neighbour correlation per bond (Eq. S11) instead, so Eq. (12) in Table I is correct only if G2 is redefined as an intensive quantity. Please make the extensive/intensive convention explicit and consistent throughout, since a reader implementing Eq. (12) with the paper's Eq. (3) definition would obtain results off by a factor of the system volume.","section":"Table I, Eq. (12)"}],"minor_comments":[{"comment":"The phrase 'from G2 alone' overstates the input required: one also needs the reference value Y(c0) and derivatives of G2 with respect to X at every intermediate c. Suggest softening this wording.","section":"Abstract; Eqs. (6) and (8)"},{"comment":"Ref. [4] duplicates Ref. [1] (Kalmutzki, Hanikel, and Yaghi) and is presumably meant to be a different paper on atomically thin materials; please correct.","section":"References"},{"comment":"The pressure benchmark is described as 'setting c0 → ∞', but Eq. (6) is written as an integral from c0 to c; please state the limiting procedure and the Tonks-Girardeau reference value explicitly.","section":"Methods, Lieb-Liniger model"},{"comment":"The extensivity argument for G2 infers extensivity of G2 from extensivity of cG2 and the kinetic energy; this inference assumes short-range interactions, so the statement that 'G2 is a thermodynamic variable' should be qualified to short-range pairwise interactions.","section":"Supplementary Information S1"},{"comment":"The critical-exponent extraction uses the exact order-parameter form m_z = (1-(|c|/hx)^{-2})^β to construct Eq. (S14); this is a useful consistency check, but the text should state that the exponent is recovered from the known form, not determined independently from the G2 data alone.","section":"Fig. 2e inset; Supplementary Information S7"}],"recommendation":"major_revision","confidential_remarks":"The main theoretical derivation is correct and the numerical benchmarks are strong. My recommendation is driven by the gap between the experimental universality claim and the exact linear-coupling assumption, together with the need for reproducible details on the heat-capacity smoothing. I see no grounds for rejection, but the experimental claims and the Table I conventions should be fixed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a methods paper with a genuinely useful core. The generalised Maxwell relations, Eqs. (5)-(8), really do imply that any thermodynamic quantity for a Hamiltonian of the form H = H0 + c G2 can be reconstructed from the integrated pair correlation alone. The derivation is almost trivial—Hellmann-Feynman plus commuting second derivatives—but the framework, especially the second-derivative relations for heat capacity and compressibility, is not in the earlier PRL and is worth having in one place. The benchmarks are honest: five model systems, exact Bethe ansatz or DMRG comparisons, no fitted parameters, and the Mott gap and critical behaviour come out correctly. The analytic entropy expression for the Yang-Gaudin gas is a nice byproduct.\n\nThe soft spots are real but not fatal. One is that the claim of being 'readily accessible in quantum simulators' needs qualification. The relation relies on the experimental control knob being exactly the linear coupling c G2 with H0 and G2 fixed. In many platforms that is only approximately true—tuning a Feshbach resonance can change the effective range, lattice-depth sweeps change t and U together, and laser-induced Stark shifts in ion traps enter H0. The paper does not estimate how large the resulting bias is for any real platform. That is a gap between the abstract and the benchmarks, which keep H0 and G2 strictly c-independent.\n\nThe second soft spot is the numerics reporting. No error bars on the benchmark comparisons, no code or data deposited yet, and the heat capacity uses a Savitzky-Golay filter on G2. None of this undermines the central derivation, but it does mean a referee would need to see error analysis and, ideally, release of the data to fully trust the accuracy claims.\n\nThird, the novelty is moderate. The entropy-pair-correlation special case is in the authors' own PRL [23], and the generalisation is mathematically immediate. But the paper is not pretending otherwise; the value is in the systematic packaging and the demonstrations.\n\nOverall: this deserves serious peer review. The core result is correct, the benchmarks are convincing as internal consistency checks, and the method is likely to be useful for theory even before the experimental calibration issue is sorted out. I would ask the authors to soften the experimental universality claim and add a sensitivity analysis for realistic control knobs.","headline":"A neat calculus identity made into a useful framework, with benchmarks that hold up; the experimental universality claim outruns what the math actually allows.","tokens_in":15937,"tokens_out":1715,"would_cite":true,"duration_ms":18658,"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":"This paper proves that every first- or second-derivative thermodynamic quantity of a many-body system can be reconstructed from the integrated pair-correlation function $G_2$ alone, via generalized Maxwell relations.","keywords":["Maxwell relations","quantum simulation","pair correlation functions","thermodynamics","quantum many-body systems","Lieb-Liniger model","Hubbard model","transverse-field Ising model"],"falsifier":"Perform a high-precision measurement of $G_2(T,V,N,c)$ in an ultracold Bose gas across a Feshbach resonance, reconstruct $P(c)$ via Eq. (6), and compare with direct pressure measurement; systematic deviations growing with $c$ would show where the exact-linear-coupling assumption fails.","tokens_in":14776,"feed_emoji":"⚛️","tokens_out":11663,"duration_ms":101433,"temperature":0.7,"pith_summary":"This paper establishes a universal route from a single measured correlation function to the full thermodynamics of a quantum many-body system. It derives generalized Maxwell relations in which any thermodynamic quantity $Y=(\\partial F/\\partial X)$ varies with interaction strength $c$ according to the derivative of the integrated pair correlation $G_2$ with respect to $X$. Because $Y(c)$ can be integrated from a known reference value $Y(c_0)$, quantities such as pressure, entropy, chemical potential, magnetization, heat capacity, and compressibility become accessible wherever $G_2$ can be measured or calculated. The authors validate the relations on the Lieb-Liniger gas, the Fermi- and Bose-Hubbard models, and the transverse-field Ising model, including at a quantum critical point. If correct, the method lets quantum simulators, which measure correlations far more easily than thermodynamics, supply direct thermodynamic comparison with materials experiments.","feed_headline":"Every thermodynamic quantity from a single pair correlation","feed_subtitle":"Generalized Maxwell relations let quantum simulators read off pressure, entropy, heat capacity from G2.","key_machinery":"The central object is the integrated pair-correlation function $G_2 = \\langle \\hat{G}_2 \\rangle$, the expectation value of the two-body operator $\\hat{G}_2$ that enters the Hamiltonian linearly as $c\\hat{G}_2$. That operator is the thermodynamic conjugate of the interaction strength $c$, and the identity $\\partial F/\\partial c = G_2$ is the load-bearing mechanism: it converts the commutativity of mixed partial derivatives of the free energy into Maxwell relations that exchange derivatives with respect to $c$ for derivatives with respect to $T$, $V$, $N$, or field components. The integral forms of these relations use a known reference value $Y(c_0)$ and a sweep of $c$, which is experimentally available in quantum simulators, to reconstruct the full thermodynamic quantity $Y(c)$.","core_discovery":"The central claim is a pair of identities. For a Hamiltonian of the form $\\hat{H} = \\hat{H}_0 + c \\hat{G}_2$, the Helmholtz free energy satisfies $\\partial F/\\partial c = G_2$ by the Hellmann-Feynman theorem, so mixed second derivatives of $F$ commute. This yields, for $Y = (\\partial F/\\partial X)_{c,\\ldots}$, the generalized Maxwell relation $(\\partial Y/\\partial c)_{X,\\ldots} = (\\partial G_2/\\partial X)_{c,\\ldots}$, and for second derivatives $(\\partial Y/\\partial c)_{X,X',\\ldots} = (\\partial^2 G_2/\\partial X\\partial X')_{c,\\ldots}$. Integrating these from $c_0$ to $c$ reconstructs $Y(c)$ from $Y(c_0)$ and the measured or calculated derivative of $G_2$. The paper benchmarks this scheme by reproducing exact results for Lieb-Liniger pressure, Fermi- and Bose-Hubbard chemical potentials with Mott gaps, transverse-field Ising magnetization and heat capacity, and by deriving a new high-temperature entropy expression for the Yang-Gaudin gas. It also shows that $G_2$ itself is a thermodynamic variable, whose derivative displays a $\\lambda$ anomaly at the quantum critical point and yields the critical exponent of the order parameter.","pith_inferences":["Beyond the paper, a natural test is thermometry: because entropy follows from $\\partial G_2/\\partial T$, a calibrated measurement of $G_2$ could assign absolute temperatures in quantum gas microscopes where standard thermometry is difficult.","For experimental tuning knobs that do not exactly realize $\\hat{H} = \\hat{H}_0 + c\\hat{G}_2$, such as Feshbach resonances that also shift the effective range, the relations will acquire correction terms; quantifying these corrections per platform would extend the method to real devices.","The lambda anomaly in $\\partial G_2/\\partial h_z$ near the critical point suggests that correlation measurements could serve as a thermodynamic probe in two-dimensional materials, where substrate-dominated heat capacity blocks conventional thermodynamic measurements."],"forward_implications":["In any quantum simulator where $c$ can be tuned continuously, one dataset of $G_2(T,V,N,c)$ suffices to reconstruct pressure, entropy, chemical potential, magnetization, heat capacity, and compressibility.","Near a continuous phase transition, derivatives of $G_2$ inherit the singular behavior of free-energy derivatives, giving a correlation-based phase-transition signature and a route to critical exponents.","The same relations reproduce Mott-insulator gaps in both Fermi- and Bose-Hubbard models from $G_2$ alone, confirming the method in strongly correlated lattice systems.","The relations work in reverse: in materials where thermodynamic measurements are straightforward but correlation measurements are not, integrated correlation functions can be inferred from standard thermodynamic data.","For Hamiltonians with higher-body interactions, the pair correlation $G_2$ is replaced by the corresponding higher-order correlation function in the same Maxwell-relation construction."],"supporting_citations":[{"why":"Supplies the standard free-energy formalism and Maxwell-relation derivation that the generalized relations extend.","marker":"[21]"},{"why":"Establishes the prior entropy-pair-correlation Maxwell relation that this work generalizes to all thermodynamic quantities.","marker":"[23]"},{"why":"Provides the Hellmann-Feynman theorem used to prove that the free-energy derivative with respect to c equals G2.","marker":"[47]"},{"why":"Defines the Lieb-Liniger model and its exact ground-state solution, the first benchmark system.","marker":"[26]"},{"why":"Provides the thermodynamic Bethe ansatz results used as exact reference for the strongly interacting Bose gas.","marker":"[31]"},{"why":"Introduces the DMRG method used to compute G2 in the benchmark models.","marker":"[32]"},{"why":"Supplies the exact transverse-field Ising solution used to benchmark magnetization, heat capacity, and the critical exponent.","marker":"[42]"},{"why":"Gives pair-correlation results for the one-dimensional Bose gas used in the numerical evaluation of G2.","marker":"[48]"}],"fun_headline_variants":["Thermodynamics from a single pair correlation","Generalized Maxwell relations for quantum simulators","Pair correlations unlock thermodynamic properties","Quantum simulators: extract entropy from G2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the interaction strength $c$ enters the Hamiltonian only as one linear term $c\\hat{G}_2$, with $c$ appearing nowhere else; if a real control knob changes anything else while tuning $c$, the central equality $\\partial F/\\partial c = G_2$ no longer holds exactly.","fun_headline_variants_meta":{"raw":{"variants":["Thermodynamics from a single pair correlation","Generalized Maxwell relations for quantum simulators","Pair correlations unlock thermodynamic properties","Quantum simulators: extract entropy from G2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000563,"raw_usage":{"total_tokens":2683,"prompt_tokens":967,"completion_tokens":1716,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":1663}},"tokens_in":583,"tokens_out":1716,"duration_ms":13125,"temperature":1.0,"reasoning_tokens":1663,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T23:07:39.520658+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform a high-precision measurement of $G_2(T,V,N,c)$ in an ultracold Bose gas across a Feshbach resonance, reconstruct $P(c)$ via Eq. (6), and compare with direct pressure measurement; systematic deviations growing with $c$ would show where the exact-linear-coupling assumption fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the standard free-energy formalism and Maxwell-relation derivation that the generalized relations extend."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the prior entropy-pair-correlation Maxwell relation that this work generalizes to all thermodynamic quantities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Lieb-Liniger model and its exact ground-state solution, the first benchmark system."},{"cited_title":"Yang and C","cited_arxiv_id":null,"evidence_quote":"Provides the thermodynamic Bethe ansatz results used as exact reference for the strongly interacting Bose gas."},{"cited_title":"Pfeuty, The one-dimensional Ising model with a transverse field, Ann","cited_arxiv_id":null,"evidence_quote":"Supplies the exact transverse-field Ising solution used to benchmark magnetization, heat capacity, and the critical exponent."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives pair-correlation results for the one-dimensional Bose gas used in the numerical evaluation of G2."}],"review_version":1}