Lieb-Mattis ordering theorem of electronic energy levels in the thermodynamic limit

This paper generalizes the Lieb-Mattis ordering theorem to fermionic mixtures with $N>2$ spinor components in the thermodynamic limit, demonstrating that the lowest-energy states within each permutation symmetry sector are well-approximated by U$(N)$ coherent states and exhibit distinct quantum phase transitions dependent on their symmetry sectors.

The Gist

Imagine a crowded dance floor where thousands of dancers (particles) are trying to find the most comfortable way to move together. In the world of quantum physics, these dancers are "fermions" (like electrons), and they have a strict rule: no two dancers can occupy the exact same spot at the same time.

This paper is about figuring out the lowest energy state (the most relaxed, comfortable arrangement) for these dancers when they have more than just two "types" of moves available.

Here is a breakdown of the paper's ideas using simple analogies:

1. The Players: From Two Colors to Many

Usually, physicists study electrons that have two "flavors" or "colors" (like spin up and spin down, or Red and Blue). This is like a dance floor where everyone is either wearing a Red shirt or a Blue shirt.

However, in modern physics (like in special atomic gases or twisted graphene), electrons can have many more colors (N components). Imagine a dance floor with Red, Blue, Green, Yellow, and even more shirt colors. The paper asks: If we have a huge crowd of these multi-colored dancers, how do they arrange themselves to be the most relaxed?

2. The Sorting Hat: Permutation Symmetry

When you have a crowd of dancers, they naturally group themselves based on how they swap places with each other.

• The "Most Symmetric" Group: Imagine a group where everyone is identical and interchangeable. If you swap any two dancers, the group looks exactly the same. This is the "most symmetric" group.
• The "Mixed" Groups: There are other groups where the dancers are a bit more picky. Swapping two specific dancers might change the "vibe" of the group slightly. These are the "mixed symmetry" groups.

In the past, scientists (using the Lieb-Mattis theorem) knew that for the simple two-color case, the "most symmetric" group always had the lowest energy (was the most comfortable). They also knew that if you took a "mixed" group and made it more symmetric (by moving dancers from the edges to the center, like pouring water from a tall glass into a wide bowl), the energy would go down.

3. The Big Question: What happens with infinite dancers?

The authors wanted to know: Does this rule still hold if we have an infinite number of dancers (the thermodynamic limit) and many more colors (N > 2)?

They used a mathematical tool called Coherent States.

• The Analogy: Imagine trying to describe the movement of a billion dancers. It's impossible to track every single one. Instead, you use a "quasi-classical" average—a smooth, flowing wave that represents the general motion of the crowd. This is what a "Coherent State" is. It's like describing the ocean as a single wave rather than tracking every water molecule.

4. The Discovery: The "Mixed Symmetry" Phase Transition

The paper finds that even with infinite dancers and many colors, the old rules mostly still apply, but with a twist:

• The Hierarchy of Comfort: Just like before, the "most symmetric" arrangement is still the most comfortable (lowest energy). However, the authors proved that even for the "mixed" groups, there is a strict order. If you can "pour" one arrangement into a more symmetric one, the more symmetric one will always have lower energy.
• New Critical Points: In the old two-color world, there was one specific moment (a critical value of interaction strength, $\lambda$) where the dancers suddenly changed their dance style (a Quantum Phase Transition).
  • The authors discovered that every single "mixed" group has its own specific moment where it changes its dance style.
  • Imagine a stadium full of people. In the "Red/Blue" section, everyone stands up at the same time when the music hits a certain beat. But in the "Red/Blue/Green" section, a different group might stand up at a slightly different beat. The paper maps out exactly when each specific group changes its behavior.

5. The Map: A New Phase Diagram

The authors created a new "map" (phase diagram) for this system.

• Old Map: Only showed the transition for the "most symmetric" group.
• New Map: Shows transitions for every possible group arrangement.
• The Result: They proved that even in this complex, infinite world with many colors, the "Lieb-Mattis" ordering rule holds true. The most symmetric groups are always the most stable, and the energy levels follow a predictable, smooth pattern as you change the interaction strength.

Summary

Think of this paper as a guidebook for a massive, multi-colored dance party.

1. The Rule: The most uniform groups of dancers are always the most relaxed.
2. The Twist: Even the less uniform groups have their own specific "moments of change" (phase transitions) depending on how many colors are involved.
3. The Proof: The authors used advanced math (Coherent States) to show that even with an infinite number of dancers, the energy levels follow a predictable, orderly pattern, confirming that the universe prefers symmetry, even in its most complex, multi-colored forms.

They tested this using a specific model (the Lipkin-Meshkov-Glick model) and confirmed that their mathematical predictions match what happens when you simulate these systems on a computer.

Technical Summary

Technical Summary: Lieb-Mattis Ordering Theorem of Electronic Energy Levels in the Thermodynamic Limit

Problem Statement
The Lieb-Mattis theorem traditionally establishes an ordering of the lowest-energy electronic states for a system of $P$ interacting fermions with $N=2$ spinor components (spin-1/2), based on the total spin $s$ and the permutation symmetry of the wavefunction. Specifically, it dictates that within a general class of symmetric Hamiltonians, states with higher symmetry (corresponding to the most symmetric Young diagram) possess lower energy.

This paper addresses the generalization of these predictions to fermionic mixtures with $N > 2$ spinor components/species. Furthermore, it investigates the behavior of these systems in the thermodynamic limit ($P \to \infty$), a regime where quantum phase transitions (QPTs) typically occur. The authors aim to determine if the Lieb-Mattis ordering holds for all irreducible representations (IRs) of the unitary group $U(N)$ arising in the $P$-fold tensor product decomposition, not just the fully symmetric sector where the global ground state resides.

Methodology
The authors employ a variational approach using $U(N)$ coherent states (CSs) defined on Grassmannian phase spaces. The methodology proceeds as follows:

1. Model Formulation: The study considers a system of $P$ fermions with $N$ components distributed on a lattice (or Landau sites). The Hamiltonian is formulated in terms of collective $U(N)$-spin operators $S_{ij}$, focusing on pairing correlations and exchange interactions.
2. Symmetry Classification: The Hilbert space is decomposed into irreducible representations (IRs) of $U(N)$ using Young diagrams. The permutation symmetry sectors are labeled by partitions $h = [h_1, \dots, h_N]$ of $P$. The authors utilize the "dominance order" (or pouring principle) $\succeq$ between Young diagrams to compare different symmetry sectors.
3. Coherent State Approximation: For each permutation symmetry sector $h$, the lowest-energy state is approximated by a $U(N)$ coherent state $|h, Z\rangle$, constructed as a rotation of the highest-weight (HW) vector $|h, 0\rangle$. In the thermodynamic limit ($P \to \infty$), these CSs are shown to faithfully reproduce the ground state energy of the critical quantum system.
4. Energy Surface Calculation: The expectation value of the Hamiltonian density, $\langle h, Z | H | h, Z \rangle$, is computed. In the limit $P \to \infty$, quantum fluctuations vanish, allowing quadratic operator terms to be replaced by products of expectation values ($s_{ij}s_{kl}$).
5. Specific Examples: The general formalism is applied to a generalized Lipkin-Meshkov-Glick (LMG) model for $N=2$ and $N=3$ components. The energy surfaces are minimized with respect to the coherent state parameters to find the lowest energy density $E^{(0)}_{\mu,\nu}(\lambda)$ for each sector.

Key Contributions

• Generalization of Lieb-Mattis to $N>2$: The paper extends the Lieb-Mattis ordering theorem to fermionic mixtures with arbitrary $N$ components, demonstrating that the energy ordering based on permutation symmetry holds even when the system is not restricted to the $N=2$ spin case.
• Thermodynamic Limit Analysis: The authors extend the analysis to the thermodynamic limit ($P \to \infty$). They show that the discrete dominance order between Young diagrams translates into monotonic properties of the energy density function with respect to continuous parameters $(\mu, \nu)$ that characterize the IRs in the large-$P$ limit.
• Mixed Symmetry Quantum Phase Transitions (MSQPT): The study introduces and characterizes "Mixed Symmetry QPTs." While standard QPTs are typically analyzed in the fully symmetric sector (ground state), this work demonstrates that QPTs occur within every permutation symmetry sector $h$ at specific critical values of the coupling constant $\lambda_c(h)$.
• Variational Framework: The work establishes $U(N)$ coherent states on Grassmannian manifolds as effective variational tools for describing the low-energy physics of $U(N)$ quantum Hall ferromagnets and interacting fermion mixtures in the thermodynamic limit.

Results

• N=2 Case: For $N=2$ (spin-1/2), the authors recover the standard result where the energy density decreases as the total spin $s$ increases (or as the parameter $\mu$ increases). A second-order QPT occurs at a critical coupling $\lambda_c(\mu) = 1/(2\mu - 1)$, which depends on the symmetry sector.
• N=3 Case: For $N=3$, the phase diagram is extended to a hyperplane defined by the interaction strength $\lambda$ and two continuous symmetry parameters $(\mu, \nu)$. The analysis reveals four distinct quantum phases separated by critical curves.
  • The energy density $E^{(0)}_{\mu,\nu}(\lambda)$ is proven to be a decreasing function of $\mu$ and an increasing function of $\nu$. This confirms that if a representation $h'$ can be "poured" into $h$ (i.e., $h \succeq h'$), then $E(h) \leq E(h')$, validating the Lieb-Mattis ordering in the thermodynamic limit.
  • The critical values $\lambda_c$ for phase transitions depend explicitly on the symmetry sector $(\mu, \nu)$.
• Numerical Verification: Numerical diagonalization of the LMG Hamiltonian for finite $P$ (e.g., $P=25, 50$) confirms the convergence of the energy density to the analytic thermodynamic limit results derived from the coherent state approximation.

Significance and Claims
The paper claims to provide a rigorous extension of the Lieb-Mattis theorem to multicomponent fermionic systems in the thermodynamic limit. By utilizing coherent states, the authors demonstrate that the ordering of energy levels by permutation symmetry is robust even as $P \to \infty$.

The primary significance lies in the concept of Mixed Symmetry QPTs. The authors argue that the standard view of QPTs, which focuses solely on the global ground state (the most symmetric sector), is incomplete. Instead, every permutation symmetry sector undergoes its own phase transition at a sector-dependent critical coupling. This enriches the phase diagram from a single parameter space ($\lambda$) to an extended space $(\lambda, h)$, where $h$ represents the mixed symmetry sector. The work suggests that this framework is applicable to $U(N)$ Quantum Hall Ferromagnets and ultracold atomic gases with high $SU(N)$ symmetry, providing a theoretical basis for understanding the low-energy structure of these complex quantum systems.