pith. sign in

arxiv: 2605.15994 · v2 · pith:K2OKFXJ4new · submitted 2026-05-15 · ❄️ cond-mat.quant-gas · cond-mat.stat-mech

Large-N scaling of Tan's contact for the harmonically trapped Tonks--Girardeau gas at finite temperature

Felipe Taha Sant'Ana This is my paper

Pith reviewed 2026-05-19 19:09 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas cond-mat.stat-mech
keywords Tan's contactTonks-Girardeau gaslarge-N scalingcanonical ensembleharmonic trapfinite temperatureFermi factor integralsensemble difference
0
0 comments X

The pith

In the large-N limit the subleading coefficient of Tan's contact for trapped Tonks-Girardeau bosons is given by universal phase-space integrals of the Fermi factor and equals the canonical-grand-canonical ensemble difference at fixed mean N

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives the finite-temperature scaling of Tan's contact in the canonical ensemble for N harmonically trapped Tonks-Girardeau bosons when N is large. The leading term recovers the local-density-approximation result via a saddle-point reduction of a contour-integral representation of the partition function. The central new result is an explicit expression for the subleading coefficient as phase-space integrals over the Fermi factor, together with its closed-form Sommerfeld and virial limits. This coefficient is identified with the difference between canonical and grand-canonical values of the contact at fixed average particle number. In the high-temperature regime the ratio of subleading to leading coefficients collapses to a universal constant traceable to Poissonian number fluctuations in the grand-canonical ensemble.

Core claim

The subleading coefficient in the large-N expansion of Tan's contact admits an explicit representation in terms of universal phase-space integrals of the Fermi factor, possesses closed-form Sommerfeld and virial limits, and is identified with the canonical-versus-grand-canonical ensemble difference at fixed mean particle number.

What carries the argument

Saddle-point reduction of the contour-integral representation of the canonical partition function to a phase-space integral with self-consistent scaled chemical potential; the subleading coefficient extracted from that integral.

If this is right

  • Padé approximants constructed from the Sommerfeld and virial limits remain uniformly accurate over the full temperature range and asymptotically correct at both ends.
  • In the high-temperature Boltzmann regime the ratio of subleading to leading coefficients collapses to a universal constant set by Poissonian particle-number statistics.
  • The entire scaling law is confirmed by direct comparison with canonical contour-integration data over the full temperature range.
  • The subleading term supplies the leading finite-N correction to the local-density result for Tan's contact.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same saddle-point technique may furnish subleading finite-N corrections for other local observables in trapped one-dimensional gases.
  • Measuring the contact while varying particle number at fixed trap frequency could isolate the predicted linear-in-N correction term.
  • The identification with ensemble difference suggests that grand-canonical calculations can be corrected to canonical accuracy by adding this explicit subleading piece.

Load-bearing premise

The saddle-point reduction of the contour-integral representation of the canonical partition function to a phase-space integral with a self-consistent scaled chemical potential remains accurate across the full temperature range for large but finite N.

What would settle it

Systematic deviation between the predicted subleading coefficient and direct numerical evaluation of the canonical contour integral for increasing N at intermediate temperatures where the saddle-point approximation is expected to be least accurate.

Figures

Figures reproduced from arXiv: 2605.15994 by Felipe Taha Sant'Ana.

Figure 1
Figure 1. Figure 1: Canonical-ensemble contact (5) as a function of the reduced temperature [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Low-τ verification of the canonical scaling law (12). The ratio RN defined in (68), with A(τ ) and B(τ ) obtained from the universal integral representations (32) and (67) using the self-consistent scaled chemical potential ξ(τ ) from (26), is plotted on τ ∈ [0.005, 0.1] for N = 10, . . . , 100. Convergence to unity (dashed line) confirms the leading and subleading scaling at low temperature, including the… view at source ↗
Figure 3
Figure 3. Figure 3: High-τ verification of the canonical scaling law: the same ratio RN as in [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Full-τ scaling verification of Eq. (12) across τ ∈ [0, 10] and N = 10, . . . , 100. (a) Ratio RN from (68), evaluated with the numerical scaling functions A(τ ) and B(τ ). (b) Ratio R (P) N from (185), evaluated with the Pad´e approximants AP(τ ) and BP(τ ) of (177) and (181). In both panels the curves collapse onto unity (dashed line), confirming the canonical scaling law from the degenerate (τ ≪ 1) to th… view at source ↗
Figure 5
Figure 5. Figure 5: Universal scaling functions A(τ ) (left panel) and B(τ ) (right panel). Solid blue: numerical evaluation of the integral representations (32) and (67) using the self-consistent ξ(τ ) from (26). Dashed black: Pad´e approximants (177) and (181). Dotted red: low-τ Sommerfeld expansions (39) and (91). Dotted green: high-τ virial expansions (51) and (103). The Pad´e approximants interpolate between the two asym… view at source ↗
read the original abstract

We derive the canonical-ensemble scaling of Tan's contact for $N$ harmonically trapped Tonks--Girardeau bosons at finite temperature in the large-$N$ limit. The leading scaling coefficient reproduces the local-density-approximation result and is obtained from a contour-integral representation of the canonical partition function followed by a saddle-point reduction to a phase-space integral with a self-consistent scaled chemical potential. The subleading coefficient is the central new object of this work: it admits an explicit representation in terms of universal phase-space integrals of the Fermi factor, has closed-form Sommerfeld and virial limits, and is identified with the canonical-versus-grand-canonical ensemble difference at fixed mean particle number. In the high-temperature Boltzmann regime the ratio of subleading to leading coefficients collapses to a universal value, traceable to the Poissonian particle-number statistics of the dilute grand-canonical gas. We construct Pad\'e approximants for both scaling functions that interpolate uniformly between the low-temperature Sommerfeld and high-temperature virial regimes; for the subleading coefficient we report a form that is uniformly accurate on our working range of temperatures and asymptotically correct beyond. The scaling law is verified against canonical contour-integration data across the full temperature range.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript derives the large-N scaling of Tan's contact for N harmonically trapped Tonks-Girardeau bosons at finite temperature in the canonical ensemble. The leading coefficient is obtained via saddle-point reduction of a contour-integral representation of the canonical partition function to a phase-space integral over the Fermi factor with self-consistent scaled chemical potential, reproducing the local-density approximation. The subleading coefficient is expressed as universal phase-space integrals of the Fermi factor, identified with the canonical-grand-canonical ensemble difference at fixed mean N, with closed-form Sommerfeld and virial limits; Padé approximants are constructed and the scaling is verified against canonical contour-integration data.

Significance. If the central derivation holds, the work supplies an explicit analytic form for the subleading correction to Tan's contact in this integrable 1D model, linking it directly to ensemble inequivalence and providing uniformly accurate interpolants between limiting regimes. The numerical verification against exact canonical data and the collapse to a universal ratio in the high-T Boltzmann regime are notable strengths that could guide experiments with ultracold atoms.

major comments (2)
  1. [§3.2, Eq. (25)] §3.2, around Eq. (25): the extraction of the subleading coefficient as the exact canonical-grand-canonical difference at fixed mean N via saddle-point reduction assumes that O(1/sqrt(N)) Gaussian and higher fluctuation corrections around the saddle are either negligible or cancel precisely in the difference; the manuscript does not provide an explicit error bound or asymptotic analysis showing this cancellation holds uniformly across the intermediate-temperature regime where neither Sommerfeld nor virial expansions apply.
  2. [§4, Eq. (32)] §4, Eq. (32): the self-consistent determination of the scaled chemical potential is used both for the leading LDA term and to define the subleading ensemble-difference integrals; it is not shown whether residual N-dependence in this self-consistency condition propagates into the claimed universal phase-space integrals at the subleading order.
minor comments (2)
  1. [Abstract] The abstract states that the Padé approximant for the subleading coefficient is 'uniformly accurate on our working range of temperatures' but does not report the maximum relative error or the precise temperature interval used for fitting; adding this quantitative information would allow readers to assess the approximant independently.
  2. [§5] In §5 the comparison plots versus canonical data would benefit from explicit indication of the particle numbers N employed and any temperature points excluded from the fit; this would clarify the range of validity of the large-N scaling.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We respond to each major comment below.

read point-by-point responses

  1. Referee: [§3.2, Eq. (25)] §3.2, around Eq. (25): the extraction of the subleading coefficient as the exact canonical-grand-canonical difference at fixed mean N via saddle-point reduction assumes that O(1/sqrt(N)) Gaussian and higher fluctuation corrections around the saddle are either negligible or cancel precisely in the difference; the manuscript does not provide an explicit error bound or asymptotic analysis showing this cancellation holds uniformly across the intermediate-temperature regime where neither Sommerfeld nor virial expansions apply.

    Authors: The saddle-point reduction of the contour-integral representation yields the leading term as the phase-space integral of the Fermi factor evaluated at the self-consistent scaled chemical potential. The subleading coefficient is obtained directly as the difference between the canonical and grand-canonical expressions for Tan's contact at fixed mean particle number; this difference isolates the ensemble-inequivalence correction at order 1/N. Gaussian fluctuations around the saddle contribute O(1/N) corrections to the free energy, but for the contact (a derivative with respect to the interaction strength) these enter only at higher order or cancel between the two ensembles. While the manuscript does not contain a rigorous uniform error bound, the scaling is verified numerically against exact canonical contour-integration results over the full temperature range, including the intermediate regime. We will add a short paragraph discussing the expected order of the neglected fluctuation terms. revision: partial

  2. Referee: [§4, Eq. (32)] §4, Eq. (32): the self-consistent determination of the scaled chemical potential is used both for the leading LDA term and to define the subleading ensemble-difference integrals; it is not shown whether residual N-dependence in this self-consistency condition propagates into the claimed universal phase-space integrals at the subleading order.

    Authors: The scaled chemical potential is fixed by the leading-order condition that the grand-canonical particle number equals N; in the scaled variables this condition is N-independent. The subleading phase-space integrals are therefore evaluated at this fixed scaled value. Any residual N-dependence in the chemical potential appears only at O(1/N) and affects terms beyond the subleading coefficient. Consequently the integrals remain universal at the order considered. We will insert a clarifying sentence in §4. revision: yes

Circularity Check

1 steps flagged

Subleading coefficient extracted via self-consistent saddle-point chemical potential and identified with ensemble difference

specific steps
  1. fitted input called prediction [Abstract]
    "The leading scaling coefficient reproduces the local-density-approximation result and is obtained from a contour-integral representation of the canonical partition function followed by a saddle-point reduction to a phase-space integral with a self-consistent scaled chemical potential. The subleading coefficient is the central new object of this work: it admits an explicit representation in terms of universal phase-space integrals of the Fermi factor, has closed-form Sommerfeld and virial limits, and is identified with the canonical-versus-grand-canonical ensemble difference at fixed mean粒子数."

    The subleading coefficient is obtained after determining a self-consistent scaled chemical potential via the saddle-point reduction (which enforces the mean particle number). Identifying this coefficient directly with the canonical-grand-canonical difference at fixed mean N makes the subleading term a direct consequence of the same self-consistency condition used to define the leading term, rather than an independent derivation.

full rationale

The derivation obtains the leading LDA term from saddle-point reduction of the contour-integral canonical partition function. The subleading coefficient is then represented as universal phase-space integrals of the Fermi factor and explicitly identified with the canonical-grand-canonical difference at fixed mean N. This identification relies on the self-consistent scaled chemical potential obtained from the same saddle-point procedure. While the explicit integral representation provides independent content, the extraction of the subleading term as the ensemble difference at fixed N is tied to the fitting of the chemical potential, raising the possibility that the claimed universality partially reflects the construction of the saddle-point approximation rather than a fully independent result. No self-citation load-bearing or ansatz smuggling is evident from the provided text.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract-only review limits visibility into explicit assumptions; main structural premises are the validity of the saddle-point approximation for large N and the identification of the subleading term with ensemble difference.

free parameters (1)
  • scaled chemical potential
    Self-consistent scaled chemical potential introduced in the saddle-point reduction of the contour integral.
axioms (1)
  • domain assumption Saddle-point approximation is valid for large but finite N across all temperatures
    Invoked to reduce the canonical partition function contour integral to a phase-space integral.

pith-pipeline@v0.9.0 · 5756 in / 1277 out tokens · 40880 ms · 2026-05-19T19:09:14.810172+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.