Excited-state quantum phase transitions and chaos in a three-level Lipkin model

This paper investigates excited-state quantum phase transitions in a three-level Lipkin-Meshkov-Glick model by combining chaos-sensitive measures like the Kullback-Leibler divergence with standard diagnostics to establish a robust framework for analyzing ESQPTs in systems exhibiting mixed regular and chaotic dynamics.

The Gist

Imagine you have a giant, complex machine made of thousands of tiny, dancing balls. In the world of quantum physics, these "balls" are particles, and how they dance together determines the state of matter (like whether it's a solid, a liquid, or something stranger).

This paper is about studying a specific, slightly more complicated version of this machine: a three-level Lipkin model. Think of it as a dance floor with three different heights (levels) where the particles can stand, rather than just two.

Here is the breakdown of what the researchers did, using simple analogies:

1. The Big Picture: Finding the "Tipping Points"

In physics, when you change a control knob (like temperature or pressure), a system can suddenly switch from one behavior to another. This is called a Phase Transition.

• Ground State Transitions: This happens when the system is at its lowest energy (like ice melting into water).
• Excited-State Transitions (ESQPTs): This is the paper's main focus. It's like asking: "What happens if we shake the machine really hard?" At certain specific energy levels, the way the particles dance changes abruptly, even though the machine is still "hot" or "excited."

2. The Problem: Order vs. Chaos

The researchers found that in this three-level machine, things get messy.

• The Regular Dance: Sometimes, the particles move in perfect, predictable patterns (like soldiers marching in step).
• The Chaotic Dance: Other times, they move randomly and unpredictably (like a mosh pit at a rock concert).
• The Messy Middle: The hardest part is the "mixed" zones where the dance is half-organized and half-messy. It's very hard to tell exactly where the "marching" stops and the "mosh pit" begins.

3. The Tools: How They Mapped the Machine

To understand this chaos, the authors used a toolbox of different "cameras" and "maps":

• The Mean Field Map (The Blueprint): First, they looked at the machine as if it were infinite and smooth (ignoring individual particles). This gave them a blueprint showing where the "tipping points" (called separatrices) are located. Think of these as invisible walls on the dance floor that separate different types of dancing.
• Poincaré Sections (The Stroboscope): They used a technique to take "snapshots" of the particles' paths.
  • If the dance is regular, the snapshots form neat, closed loops (like drawing a circle).
  • If the dance is chaotic, the snapshots look like a scattered cloud of dust.
• Peres Lattices (The Grid): They plotted the energy of the particles against how much they occupy each level.
  • In a regular system, the dots line up in a perfect grid (like a chessboard).
  • In a chaotic system, the dots scatter randomly, breaking the grid.
• The "Chaos Meter" (KL Divergence): This is their secret weapon. They invented a way to measure how much the "mosh pit" looks like a true random mess versus a structured pattern. They found that a specific math tool called Kullback-Leibler divergence was the best at spotting the exact moment the system switched from ordered to chaotic.

4. The Discovery: The "Zones" of the Dance Floor

By combining all these tools, they successfully mapped out the three-level machine for high-energy states. They found that the dance floor is divided into four distinct zones as you go from low energy (calm) to high energy (wild):

1. Quasi-Integrable (The Calm Zone): Near the bottom, everything is orderly.
2. Chaotic (The Mosh Pit): As you add energy, the system hits a "separatrix" wall and suddenly becomes completely chaotic.
3. Quasi-Chaotic (The Confused Zone): As you add even more energy, it gets messy again, but not fully chaotic. It's a weird, mixed state.
4. Quasi-Integrable (The Calm Returns): Surprisingly, at the very highest energies, the system starts to organize itself again!

5. Why Does This Matter?

You might ask, "Who cares about a theoretical dance floor?"

• Future Computers: Today's quantum computers mostly use "qubits" (two-level systems, like a light switch: on or off). This research explores "qutrits" (three-level systems, like a dimmer switch: off, dim, bright). Understanding how these three-level systems behave is crucial for building the next generation of super-powerful quantum computers.
• Sensors: It helps us design better sensors that can detect tiny changes in the environment.

The Takeaway

The authors successfully built a "rulebook" for understanding how complex, three-level quantum systems behave. They showed that even in a chaotic, messy environment, there are hidden structures and specific "tipping points" (ESQPTs) that dictate when the system will organize or disorganize. They proved that by using the right mathematical "lenses," we can predict exactly where the chaos begins and ends, paving the way for better quantum technologies.

Technical Summary

1. Problem Statement

While Excited-State Quantum Phase Transitions (ESQPTs) have been extensively characterized in two-level models (e.g., the standard Lipkin-Meshkov-Glick model), their analysis in systems with mixed regular and chaotic dynamics remains challenging.

• The Challenge: In complex systems, the spectral signatures of ESQPTs (typically level clustering) are often obscured by avoided crossings characteristic of chaotic dynamics.
• The Gap: There is a lack of a robust framework to distinguish between regular, chaotic, and quasi-chaotic dynamical regimes in excited states of multi-level systems, particularly where multiple separatrices partition the energy spectrum.
• Objective: To investigate ESQPTs within the three-level Lipkin-Meshkov-Glick (LMG) model, utilizing an enlarged Hilbert space to study the interplay between ESQPTs and quantum chaos.

2. Methodology

The authors employ a combination of analytical mean-field theory, numerical diagonalization, and dynamical system analysis.

A. The Model

• Hamiltonian: A three-level LMG model with $N$ identical particles and infinite-range interactions, governed by the $U(3)$ Lie algebra.
• Control Parameter: The interaction strength is renormalized into a dimensionless parameter $\lambda \in [0, 1]$.
• Symmetry: The system possesses a discrete $Z_2 \times Z_2 \times Z_2$ parity symmetry, dividing the Hilbert space into independent sectors based on particle number parity in each level.
• Basis: Calculations are performed in the fully symmetric irreducible representation of $U(3)$ (Bose-Einstein-Fock basis), reducing the dimension to $d = (N+1)(N+2)/2$.

B. Analytical Framework (Mean Field)

• Coherent States: The authors utilize $U(3)$-spin coherent states (3SCSs) to derive the energy surface $E(z, \lambda)$ in the thermodynamic limit ($N \to \infty$).
• Separatrices: By minimizing the energy surface and finding stationary points ($\nabla E = 0$), they identify 14 distinct ESQPT separatrices. Four of these are highlighted as critical for defining dynamical regions:
  • $f_1, f_2$: Associated with ground-state phase transitions and low-energy excitations.
  • $f_3, f_4$: Defining boundaries for higher-energy chaotic and quasi-integrable regions.

C. Numerical and Dynamical Analysis

To characterize the spectrum for finite $N$ (specifically $N=15$ and $N=100$), the authors employ:

1. Poincaré Sections: To visualize classical trajectories in phase space, distinguishing between closed curves (integrable/quasi-integrable) and scattered points (chaotic).
2. Peres Lattices: Plotting expectation values of conserved observables ($\langle S_{ii} \rangle$) against energy to detect regularity (lattice structure) vs. chaos (scattering).
3. Spectral Statistics:
   • Unfolding: Smoothing the Density of States (DOS) to remove global trends.
   • Nearest-Neighbor Spacing Distribution (NNSD): Comparing level spacings to Wigner (chaotic) and Poisson (integrable) distributions.
   • Chaos Measures: Calculating the Degree of Chaos ($\eta$) and Kullback-Leibler (KL) Divergence to quantify the deviation from integrability.
   • Participation Ratio (PR): Measuring the delocalization of eigenstates in the Fock basis.
   • Overlap Distribution: Analyzing the scalar product between eigenstates and basis states to distinguish Gaussian Orthogonal Ensemble (GOE) behavior from integrable delta-like distributions.

3. Key Results

A. Phase Structure and Separatrices

• The mean-field analysis reveals three ground-state phases separated by second-order QPTs at $\lambda = 1/3$ and $\lambda = 3/5$.
• In the excited spectrum, four main separatrices ($f_1, f_2, f_3, f_4$) partition the energy spectrum into distinct dynamical regimes, particularly for high interaction strengths ($\lambda \approx 0.98$):
  1. Quasi-integrable: Low energy ($E < f_2$) and very high energy ($E > f_3$).
  2. Chaotic: Intermediate energy ($f_2 < E < f_1$).
  3. Quasi-chaotic: A transitional region ($f_1 < E < f_3$).

B. Chaos and ESQPT Signatures

• Poincaré & Peres: Visualizations confirm that regions between separatrices $f_2$ and $f_1$ exhibit fully chaotic dynamics (scattered points), while regions outside these boundaries show regular, lattice-like structures.
• Spectral Statistics:
  • The NNSD in the chaotic region fits the Wigner distribution, while quasi-integrable regions fit the Poisson distribution.
  • The Degree of Chaos ($\eta$) and KL Divergence serve as effective quantitative markers. The KL divergence, in particular, shows a sharp drop at separatrix $f_2$, clearly demarcating the transition from quasi-integrable to chaotic regimes.
• Participation Ratio (PR): The PR exhibits clustering in quasi-integrable regions and spreads significantly in chaotic regions. Local minima in the PR correspond to the ESQPT separatrices.

C. Limitations

• The clear separation of dynamical regimes is most pronounced at high interaction strengths ($\lambda \approx 1$).
• For intermediate $\lambda$, the separatrices merge or cross, making the classification of dynamical regions difficult due to the complexity of the energy surface stationary points.

4. Key Contributions

1. Framework for 3-Level Systems: The paper establishes a systematic method to analyze ESQPTs in three-level systems, extending previous work limited to two-level models.
2. Chaos-ESQPT Connection: It explicitly links the location of ESQPT separatrices to the onset of quantum chaos, demonstrating that separatrices act as boundaries between regular and chaotic dynamical phases.
3. Robust Diagnostic Tools: The authors validate the Kullback-Leibler divergence and the Degree of Chaos ($\eta$) as superior static indicators for detecting ESQPTs in chaotic environments compared to traditional methods alone.
4. Analytical-Numerical Bridge: The work successfully uses mean-field stationary points to predict the location of chaotic regions in the finite-$N$ quantum spectrum.

5. Significance and Future Outlook

• Fundamental Physics: The study deepens the understanding of how critical phenomena (ESQPTs) coexist with and influence quantum chaos in complex many-body systems.
• Quantum Technologies: The findings have direct implications for quantum sensing (e.g., NV-centers) and quantum computation using qutrits (three-level systems). As the field moves beyond qubits to more flexible level structures, understanding the stability and chaotic behavior of these systems is crucial.
• Scalability: The analytical and numerical techniques developed here provide a template for studying other multi-level models, such as extensions of the Jaynes-Cummings, Dicke, and Agassi models, and other permutation symmetry sectors.

In summary, this paper provides a comprehensive characterization of the interplay between excited-state phase transitions and chaos in a three-level Lipkin model, offering a validated toolkit for future investigations into complex quantum many-body dynamics.