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

This paper derives the large-$N$ scaling of Tan's contact for harmonically trapped Tonks--Girardeau bosons at finite temperature by identifying a new subleading coefficient that quantifies the canonical-versus-grand-canonical ensemble difference, providing explicit universal representations and accurate Padé approximants that interpolate between low- and high-temperature regimes.

The Gist

Imagine a crowded dance floor where the dancers are tiny, invisible particles called bosons. In this specific scenario, these particles are in a "Tonks–Girardeau" state, which is a fancy way of saying they are extremely grumpy and refuse to touch each other. If two try to occupy the same spot, they bounce off with infinite force, like hard-core billiard balls.

The paper investigates a specific property of this crowd called Tan's Contact. Think of this "Contact" as a measure of how often these grumpy dancers bump into each other. In the quantum world, these bumps aren't just physical collisions; they create a specific "tail" in the way the particles move, a signature that tells us everything about their interactions.

The author, Felipe Taha Sant'Ana, is trying to figure out exactly how this "bumping rate" changes based on two things:

1. How many dancers are on the floor ($N$): The paper looks at the "Large-N" limit, meaning a very large crowd.
2. How hot the dance floor is ($T$): From freezing cold (where quantum rules dominate) to hot and chaotic (where classical rules take over).

The Main Discovery: A Two-Part Formula

The paper derives a mathematical recipe (a scaling law) to predict the "bumping rate" for a huge crowd. The recipe has two main ingredients, like a cake with a base layer and a frosting layer:

1. The Big Layer (The Leading Term):
This is the main part of the answer. It scales with the number of particles to the power of 2.5 ($N^{5/2}$).

• The Analogy: Imagine the size of the dance floor. As you add more dancers, the total number of potential collisions grows very fast. This part of the formula is what you would expect if you just looked at the average density of the crowd. It matches what scientists have known for a long time using a method called the "Local Density Approximation" (essentially, treating the crowd as a smooth fluid).

2. The Small Layer (The Subleading Term):
This is the new discovery of the paper. It is a smaller correction that scales with $N^{1.5}$ ($N^{3/2}$).

• The Analogy: This is the "fine print." While the big layer tells you the average behavior, this small layer accounts for the fact that the number of dancers is fixed.
• The "Fixed vs. Floating" Problem: In physics, you can calculate things in two ways:
  • Grand-Canonical: You imagine the dance floor is connected to a giant reservoir. Dancers can wander in and out freely. The number of dancers fluctuates.
  • Canonical: You lock the door. The number of dancers is fixed exactly at $N$.
  • The paper shows that the "Small Layer" is exactly the difference between these two scenarios. Because the door is locked in the real experiment (Canonical), the particles have to "adjust" their behavior slightly compared to the floating scenario. This adjustment creates a specific, predictable correction to the bumping rate.

The Temperature Journey

The paper maps out how this formula works across different temperatures:

• The Freezing Cold (Low Temperature):
  The dancers are very organized, almost like a perfect crystal. The "Small Layer" correction is negative and grows linearly with temperature. It's like a subtle shiver in the crowd that changes how they bump.
• The Hot Chaos (High Temperature):
  The dancers are moving wildly and rarely bumping. In this "Boltzmann" regime, the paper finds a surprising universal truth: the "Small Layer" becomes exactly the negative of the "Big Layer."
  • The Metaphor: It's as if the correction cancels out the main effect in a specific ratio. This happens because, in the hot, dilute gas, the number of particles behaves like a random coin flip (Poissonian statistics). The math shows that the "locked door" effect is exactly equal and opposite to the main crowd size effect in this extreme heat.

The "Universal" Bridge

One of the paper's most practical achievements is creating Padé approximants.

• The Analogy: Imagine you have a map of the terrain at the very bottom of a valley (cold) and at the very top of a mountain (hot), but you don't have a map for the middle. The author builds a smooth, curved bridge (a mathematical function) that connects the bottom and the top perfectly.
• This bridge allows scientists to calculate the "bumping rate" for any temperature in between, without needing to run complex, slow computer simulations every time. The paper provides these formulas so experimentalists can use them immediately.

Why This Matters (According to the Paper)

The paper doesn't claim to cure diseases or build new engines. Its value is purely in precision physics.

• Recent experiments have finally been able to measure this "Tan's Contact" directly in 1D gases.
• Before this paper, scientists had a good guess for the main part of the answer, but they lacked the precise correction for the "fixed number of particles" scenario.
• This paper provides the exact "correction factor" needed to match theory with those new, high-precision experiments. It tells experimentalists: "If you have $N$ particles at temperature $T$, here is the exact number you should see, including the subtle difference caused by locking the particle count."

In short, the paper takes a complex quantum crowd, breaks down its "bumping rate" into a main effect and a subtle correction, explains exactly why that correction exists (the difference between a fixed crowd and a floating one), and provides a smooth mathematical map to predict it at any temperature.

Technical Summary

Technical Summary: Large-N Scaling of Tan's Contact for the Harmonically Trapped Tonks–Girardeau Gas at Finite Temperature

Problem Statement
The paper addresses the theoretical determination of Tan's contact, $C$, for a one-dimensional gas of $N$ bosons with infinite repulsive interactions (the Tonks–Girardeau or TG gas) confined in a harmonic trap at finite temperature. While Tan's contact is a universal observable governing short-range correlations and thermodynamic identities in strongly interacting quantum gases, its precise scaling with particle number $N$ and temperature $T$ in the canonical ensemble (fixed $N$) has been a subject of investigation, particularly regarding the transition from few-body to many-body regimes. Previous work established the contact's behavior in the grand-canonical ensemble (GCE) using the local-density approximation (LDA), but the specific finite-$N$ corrections arising from the canonical constraint were not fully characterized in the large-$N$ limit.

Methodology
The authors derive the large-$N$ scaling of the contact by employing a rigorous saddle-point analysis of the canonical partition function. The methodology proceeds through the following steps:

1. Kernel Representation: The contact is expressed as a canonical thermal average of a specific integrand $F_N(x)$, which is constructed from the diagonal and near-diagonal elements of the fermionic kernel (via the Bose-Fermi mapping). This integrand is represented as a contour integral over the complex fugacity $z$, weighted by the grand partition function $\Xi(z)$.
2. Saddle-Point Reduction: In the limit of large $N$ at fixed reduced temperature $\tau = T/T_F$, the contour integrals are evaluated using the saddle-point method. The leading order is obtained by evaluating the kernel functional at the saddle point $z^*$, which corresponds to the grand-canonical solution with a self-consistent scaled chemical potential $\xi(\tau)$.
3. Fluctuation Analysis: To capture subleading corrections, the authors expand the integrand to include Gaussian fluctuations around the saddle point. This involves calculating particle-number cumulants ($\kappa_2, \kappa_3$) and their relation to the derivatives of the kernel functional with respect to the chemical potential.
4. Asymptotic and Numerical Evaluation: The derived scaling functions are evaluated analytically in the low-temperature (Sommerfeld, $\tau \ll 1$) and high-temperature (Boltzmann/Virial, $\tau \gg 1$) limits. For intermediate temperatures, the universal phase-space integrals are computed numerically.
5. Verification: The derived scaling law is verified against direct numerical evaluation of the canonical contour integrals for particle numbers $N$ up to 100 across the full temperature range $\tau \in [0, 10]$.

Key Contributions and Results
The primary result is a two-term large-$N$ expansion for the canonical contact:
$$C_N(\tau) = A(\tau) N^{5/2} + B(\tau) N^{3/2} + O(N^{1/2})$$
where $A(\tau)$ and $B(\tau)$ are universal functions of the reduced temperature $\tau$.

• Leading Term ($A(\tau)$): The coefficient $A(\tau)$ reproduces the known grand-canonical LDA result. The authors re-derive it using the canonical contour representation and provide closed-form asymptotic expansions for both the Sommerfeld and virial limits.
• Subleading Term ($B(\tau)$): This is the central new object of the work. The authors derive an explicit first-principles representation for $B(\tau)$ in terms of universal phase-space integrals of the Fermi factor.
  • Physical Interpretation: $B(\tau)$ is identified as the precise difference between the canonical and grand-canonical contacts at fixed mean particle number. It arises from the constraint of fixed $N$, which suppresses particle-number fluctuations present in the GCE.
  • Asymptotic Behavior:
    • At low temperature ($\tau \ll 1$), $B(\tau)$ is linear in $\tau$ with a negative slope.
    • At high temperature ($\tau \gg 1$), $B(\tau)$ approaches $-A(\tau)$. This leads to a universal ratio $B/A \to -1$ in the Boltzmann regime, a consequence of the Poissonian particle-number statistics of the dilute GCE.
• Padé Approximants: The authors construct closed-form Padé approximants for both $A(\tau)$ and $B(\tau)$ that uniformly interpolate between the low- and high-temperature asymptotic regimes. These approximants are reported to be uniformly accurate across the working temperature range and asymptotically correct beyond it.
• Numerical Verification: The scaling law is verified against canonical contour-integration data. The ratio of the numerical contact to the scaling prediction collapses to unity as $N$ increases, confirming the $N^{5/2}$ and $N^{3/2}$ scaling exponents and the accuracy of the derived coefficients.

Significance
The paper claims to provide a quantitative target for finite-temperature experiments on trapped one-dimensional Bose gases, specifically in light of recent direct measurements of the contact in Lieb-Liniger gases. By identifying the subleading $N^{3/2}$ term as a finite-$N$ ensemble-correspondence effect, the work extends previous few-body canonical analyses to the many-body limit. The results clarify the origin of ensemble differences in strongly interacting 1D systems and offer a precise theoretical framework for comparing canonical and grand-canonical descriptions in the presence of a harmonic trap. The authors note that their analysis is valid for $\tau \gg N^{-2/3}$, distinguishing it from the strictly zero-temperature regime where different finite-$N$ edge corrections (scaling as $N^{3/4}$) dominate.