Ordering in statistical systems on the way to the thermodynamic limit

This paper introduces the concept of "order indices" to quantitatively describe the growth of preordering in finite statistical systems as they approach the thermodynamic limit, illustrating the approach through mean-field models of phenomena such as Bose-Einstein condensation, superconductivity, magnetization, and crystallization.

The Gist

Imagine you are trying to understand how a crowd of people suddenly starts marching in perfect unison. In the world of physics, this "marching" is called order, and it's what happens when things like magnets align, water freezes into ice, or atoms clump together to form a superfluid.

For a long time, physicists had a rule: you can only truly say this "marching" has started if the crowd is infinite. If the crowd is finite (even if it's a million people), the math says they can't be perfectly synchronized. This is the "thermodynamic limit"—a theoretical state where the number of particles is infinite.

But here's the problem: In the real world, we never have infinite crowds. We have huge, but finite, systems. So, how does the order grow as we add more and more people to the crowd? Does it just pop into existence instantly when we hit a certain number, or does it build up gradually?

This paper by Yukalov and Yukalova says: It builds up gradually. And they have invented a new way to measure exactly how much "marching" is happening at any stage of growth.

The New Tool: The "Order Index"

Think of the Order Index as a "Synchronization Score."

• Score of 0 (or negative): The crowd is chaotic. Everyone is walking in random directions. There is no order.
• Score of 1: The crowd is perfectly synchronized. Everyone is marching in lockstep. This is the "thermodynamic limit" (perfect order).
• Scores between 0 and 1: The crowd is starting to get organized. Some people are looking at each other and copying steps, but it's not perfect yet.

The authors show that as you increase the size of the system (add more particles), this score doesn't jump from 0 to 1. It climbs steadily. The "Order Index" tells you exactly how high the score is for a specific system size.

How They Measure It: The "Echo" Analogy

To measure this score, the authors look at correlations. Imagine you shout in a large hall.

• In a small, chaotic room, your shout dies out immediately. The "echo" (correlation) is short.
• In a perfectly ordered, giant hall, your shout might bounce around and be heard clearly across the entire room. The "echo" is long.

The authors use a mathematical tool called a Reduced Density Operator. Think of this as a device that measures how far the "echo" of one particle reaches to influence its neighbors.

• If the echo is short, the Order Index is low.
• If the echo stretches across the whole system, the Order Index is high (close to 1).

They apply this same logic to different types of "echoes":

1. Single particles: How does one atom influence another?
2. Pairs of particles: How do two atoms dance together?

The Four Examples They Tested

To prove their idea works, they ran simulations on four different physical phenomena, treating them like different types of crowds:

1. Bose-Einstein Condensation (The "Super-Flow" Crowd)

• The Scene: Atoms cooled down so much they all decide to move as a single giant wave.
• The Finding: As you add more atoms, the "synchronization score" rises. However, if the atoms interact too strongly (like a rowdy crowd pushing each other), it takes more people to get the score up. Strong interactions make it harder to organize.

2. Superconductivity (The "Dancing Pairs" Crowd)

• The Scene: Electrons usually run around chaotically. But in a superconductor, they pair up and dance in perfect sync.
• The Finding: Here's a twist. If you look at individual electrons, they still look chaotic (Score ~ 0). But if you look at the pairs, the score shoots up! The "Order Index" for pairs reaches 0.5 (halfway to perfection) as the system grows. This explains why superconductivity is a "pairing" phenomenon, not a single-particle one.

3. Magnetization (The "Compass" Crowd)

• The Scene: Tiny magnets (spins) that want to point in the same direction.
• The Finding: As the system grows, the "Compass Score" climbs. Even if the magnetism is weak (a small fraction of people pointing the right way), the score grows steadily with the size of the system until it hits the maximum.

4. Crystallization (The "Grid" Crowd)

• The Scene: Liquid turning into a solid crystal.
• The Finding: In a liquid, particles are everywhere. In a crystal, they are locked in a grid. The authors measured how much the density fluctuates from the average. As the system grows, the "Grid Score" rises, showing the transition from a messy liquid to an ordered solid.

The Big Picture

The main takeaway is simple: Order is not a switch that flips on; it's a dimmer that slowly turns up.

Before the system becomes "infinite" (which is impossible in reality), it passes through a stage where it is "large but finite." In this stage, order is already forming, and the Order Index is the ruler we use to measure exactly how much order exists.

• Small systems: Low score, chaotic.
• Medium systems: The score climbs, order begins to appear.
• Huge systems: The score gets very close to 1, approaching perfect order.

This paper provides the mathematical "ruler" to measure that growth, proving that we don't need to wait for an infinite universe to see order; we can see it growing right now in large, finite systems.

Technical Summary

Technical Summary: Ordering in Statistical Systems on the Way to the Thermodynamic Limit

Problem Statement
The paper addresses a fundamental conceptual gap in statistical mechanics regarding the emergence of order and symmetry breaking. While the rigorous mathematical definition of phase transitions and order parameters (such as magnetization or condensate fraction) strictly requires the thermodynamic limit ($N \to \infty, V \to \infty, \rho = \text{const}$), physical systems are finite. The authors note that order does not appear instantaneously upon reaching infinity; rather, a "preordering" process must occur as the system size increases. A principal question remains unclear: how does order quantitatively grow in finite systems as they approach the thermodynamic limit? Existing methods, such as finite-size scaling used in numerical extrapolations or the Bogolubov quasi-average method (which requires taking the limit first before removing symmetry-breaking terms), do not provide a direct quantitative description of the growth of order in a sequence of finite systems.

Methodology
The authors propose a quantitative framework based on order indices derived from the properties of trace-class operators, specifically reduced density operators and correlation operators.

1. Definition of Order Indices: For a semi-positive trace-class operator $\hat{A}$, the order index $\omega(\hat{A})$ is defined as:
   $$\omega(\hat{A}) \equiv \frac{\log ||\hat{A}||}{\log |\text{Tr}\,\hat{A}|}$$
   where $||\hat{A}||$ is the operator norm (the largest eigenvalue) and $\text{Tr}\,\hat{A}$ is the trace.

   • Interpretation: The index characterizes the relationship between the operator norm and the system size (particle number $N$).
     • If $\omega \le 0$, there is no order (correlation length $\ll$ system size).
     • If $0 < \omega < 1$, order is developing (correlation length is comparable to but smaller than system size).
     • If $\omega \to 1$, long-range order is established (correlation length $\sim$ system size).
   • The authors generalize this from reduced density matrices (standard in statistical mechanics) to arbitrary correlation operators.

2. Physical Picture: The order index is linked to the correlation length $l_{\text{cor}}$. For a first-order density operator, $||\hat{\rho}_1|| \sim (l_{\text{cor}}/a)^3$, where $a$ is the inter-particle distance. As the system size $L$ increases, if $l_{\text{cor}}$ grows to become comparable to $L$, the norm scales with $N$, driving $\omega$ toward 1.

3. Application to Models: The authors apply this formalism to four specific phenomena using the mean-field approximation to ensure analytical tractability and clarity:

   • Bose-Einstein Condensation (BEC)
   • Superconducting Transition (BCS theory)
   • Magnetic Transition (Spin systems)
   • Crystal-Liquid Transition

Key Results

• Bose-Einstein Condensation:

  • The first-order order index $\omega(\hat{\rho}_1)$ is calculated for a uniform Bose gas with contact interactions.
  • For small system sizes $N$, $\omega$ can be negative (indicating no order). As $N$ increases, $\omega$ grows toward 1.
  • Interactions suppress the rate of order growth; stronger interactions (higher gas parameter $\gamma$) result in a slower approach to $\omega \approx 1$.
  • The second-order index $\omega(\hat{\rho}_2)$ also approaches 1 in the thermodynamic limit, but its growth is similarly hindered by interactions.

• Superconducting Transition:

  • A crucial distinction is found between single-particle and pair correlations.
  • The first-order index $\omega(\hat{\rho}_1)$ remains near zero ($\approx 0$) even in the thermodynamic limit because the single-particle occupation numbers do not scale with $N$ (no macroscopic occupation of a single state in the same sense as BEC).
  • However, the second-order index $\omega(\hat{\rho}_2)$, which characterizes pair correlations, grows from zero to 1/2 in the thermodynamic limit. This reflects that the order is carried by pairs, not single particles.

• Magnetic Transition:

  • Using spin correlation operators, the first-order index $\omega(\hat{C}_1)$ is derived.
  • For any non-zero magnetization $M$, the index grows from negative values to 1 as $N \to \infty$.
  • The second-order index $\omega(\hat{C}_2)$ coincides with the first-order index, indicating that the ordering of spins is captured at the single-particle correlation level.

• Crystal-Liquid Transition:

  • Order is defined via density deviation operators.
  • Similar to magnetization, the first-order index $\omega(\hat{C}_1)$ grows to 1 in the thermodynamic limit if the density is non-uniform ($D > 0$).

Significance and Claims
The paper claims to provide a rigorous mathematical tool to describe the process of ordering in finite systems, bridging the gap between microscopic finite systems and the macroscopic thermodynamic limit.

• Quantification of Preordering: The order index allows for the quantification of "preordering" states where the system is large but finite. It demonstrates that order is not an all-or-nothing property but a continuous growth dependent on system size.
• Distinction of Order Types: The framework successfully distinguishes between different types of ordering. It correctly identifies that BEC and magnetization are characterized by a first-order index approaching 1, while superconductivity is characterized by a vanishing first-order index but a second-order index approaching 1/2.
• Finite-Size Effects: The authors emphasize that for sufficiently large $N$ ($N \gg 1$), boundary effects become negligible, and the order index provides a clear measure of how close the system is to the thermodynamic limit.
• Conjectures: The paper concludes with four conjectures regarding the hierarchy of order indices, suggesting that if order exists in a lower-order operator, it generally exists in higher orders, but the converse is not necessarily true (e.g., order may exist in pairs but not single particles).

The authors explicitly state that the paper does not aim to prove the existence of phase transitions (a separate difficult problem) nor to study critical regions or nucleation dynamics. Instead, the focus is strictly on the quantitative description of the growth of equilibrium order as system size increases.