General phase diagram features of superradiant phase transitions

This paper establishes a general phase diagram feature for superradiant phase transitions, demonstrating that the system originates in a normal phase and undergoes a single transition along the radial direction of coupling parameters, a finding derived via a concise mean-field method that accounts for diverse interactions, multimode effects, and disorder.

The Gist

Imagine a giant dance floor where thousands of tiny dancers (atoms or "qubits") are trying to move in sync with a spotlight (light or "photons"). Usually, they dance randomly, each doing their own thing. This is the Normal Phase (NP)—a state of calm chaos where nothing special is happening.

But sometimes, if the music gets loud enough or the connection between the dancers and the spotlight gets strong enough, something magical happens. Suddenly, everyone locks into a perfect, synchronized rhythm. The spotlight glows intensely, and the dancers move as one giant entity. This is the Superradiant Phase (SP)—a state of collective, super-powered harmony.

This paper is like a master mapmaker trying to draw the rules for when this "dance switch" happens. The authors, Wen Zhao, Junlong Tian, and Jie Peng, looked at many different, complicated versions of this dance floor. Some had multiple spotlights, some had dancers who could talk to each other, and some had weird, twisted rules for how they interact.

Here is the simple breakdown of their findings:

1. The Universal "One-Way Street" Rule

The most surprising thing the authors found is that no matter how complicated the dance floor is, the transition from "chaos" to "harmony" follows a simple, universal pattern.

Imagine you are walking away from the center of a city (the Normal Phase) in any direction you choose. As you walk further out (increasing the strength of the connection between dancers and light), you will eventually cross a specific line. Once you cross that line, you enter the "Harmony Zone" (Superradiant Phase).

The key finding: You can never walk back. Once you are in the Harmony Zone, you stay there. You don't wander back into chaos just because you keep walking. The system has a "one-way street" feature: it starts in the normal state, crosses a boundary, and stays superradiant forever after.

2. The "Teamwork" Shortcut

The paper also discovered a cool trick about how many spotlights (modes) are needed.

• The Old Way: If you only have one spotlight, the dancers need to be incredibly strong and the connection very intense to sync up. It's like trying to get a whole stadium to cheer in unison with just one megaphone; it's hard.
• The New Finding: If you have many spotlights working together, the dancers can sync up much more easily. Even if the connection to each individual spotlight is weak, the combined effort of all the spotlights creates a "teamwork effect." This allows the super-radiant state to happen even when the individual connections are in a "strong coupling" regime (a technical term meaning the interaction is powerful but not yet "ultra-powerful"). It's like getting a crowd to cheer by having ten people with small megaphones instead of one giant one.

3. The "Messy Room" Effect

The authors also looked at what happens if the dancers aren't all identical—maybe some are taller, some are shorter, or the floor is a bit uneven (this is called "disorder").
They found that this messiness doesn't stop the dance from happening, but it does move the line. If the room is messy, the point where the dancers switch from chaos to harmony shifts. You might need to turn the music up a little louder (or a little softer) to get them to sync up, depending on how messy the room is.

4. How They Figured It Out

Instead of getting lost in thousands of complex math equations for every specific type of dance floor, the authors used a clever shortcut. They treated the whole system like a landscape with hills and valleys.

• The Normal Phase is like a ball sitting at the very bottom of a valley at the center of the map.
• The Superradiant Phase is like the ball rolling down to a new, deeper valley on the side.

They proved that no matter how you twist the rules of the dance, the ball always starts in the center valley. As you turn up the "volume" (the coupling strength), the center valley gets shallower, and a new, deeper valley appears on the side. Once the ball rolls into that new valley, it stays there. This simple picture allowed them to predict exactly when the switch happens for all these different complex models.

Summary

In short, this paper says: "Don't worry about how complicated your light-matter system is. If you turn up the connection strength, the system will always start in a normal state, cross a specific boundary, and lock into a super-radiant state, never to return. And if you use multiple light sources, it's actually easier to get them to sync up."

This helps scientists understand and predict these quantum phenomena without needing to solve a different, impossible math problem for every single new experiment they design.

Technical Summary

Technical Summary: General Phase Diagram Features of Superradiant Phase Transitions

Problem Statement
Superradiant phase transitions (SPTs) in light-matter systems have been extensively studied across various models, including the Dicke model, the Quantum Rabi model (QRM), and their generalizations involving anisotropic couplings, multi-photon interactions, inter-cavity hopping, and qubit-qubit interactions. These diverse coupling processes lead to complex and varied phase diagram structures. A central question remains: how does a general light-matter system transition from the normal phase (NP) to the superradiant phase (SP) within a multidimensional coupling parameter space? While specific features have been identified in limited cases (e.g., single-mode dipolar coupling), a unified understanding of the phase diagram topology for general multiqubit, multimode systems with arbitrary coupling processes is lacking.

Methodology
The authors employ a mean-field approximation, validated in the thermodynamic limit ($N \to \infty$) or the classical oscillator limit, to analyze a general Hamiltonian encompassing a wide class of light-matter interactions. The Hamiltonian includes $M$ field modes and $N$ qubits, incorporating:

• Anisotropic couplings (linear and nonlinear).
• One- and two-photon interactions.
• Stark shifts and inter-cavity hopping.
• Qubit-qubit interactions.

The study utilizes the canonical partition function expanded in the Glauber coherent state basis. By rescaling the mean photon number and applying the Laplace method in the limit where $\beta N \Omega \to \infty$ (where $\beta$ is inverse temperature, $N$ is the number of qubits, and $\Omega$ is a generalized qubit energy), the partition function is approximated by the global minimum of a Landau potential, $\phi(\vec{z})$. Here, $\vec{z}$ represents the rescaled coherent state parameters. The rescaled mean photon number serves as the order parameter, corresponding to the location of the global minimum of this potential.

The authors propose a concise method to determine the phase boundary and the order of the transition:

1. Analyze the Landau potential $\phi(\vec{z})$ to find its global minimum.
2. Determine if the origin ($\vec{z}=0$) is the global minimum (NP) or if a non-zero solution becomes energetically favorable (SP).
3. The phase boundary is identified by solving $\phi(\vec{z}_c) = \phi(0)$, where $\vec{z}_c$ is an extremum.
4. The order of the transition is determined by the behavior of the order parameter at the boundary: continuous change implies a second-order transition, while a discontinuous jump implies a first-order transition.

Key Contributions and Results
The paper establishes a general phase diagram feature for the described class of systems:

• Radial Transition: In the coupling parameter space defined by a vector $\vec{\lambda}$, the system is naturally in the NP at the origin ($\vec{\lambda} \approx 0$). As the magnitude of the coupling vector $|\vec{\lambda}|$ increases radially, the system undergoes a transition to the SP and remains in the SP thereafter; it does not revert to the NP.
• Multimode Cooperative Behavior: For multimode systems, the condition for SPT depends on the norm of the coupling vector, $|\vec{\lambda}| = |(\lambda_1, \lambda_2, \dots)|$. This implies that the coupling strength required for each individual mode can be significantly reduced compared to the single-mode case, potentially allowing SPTs to occur even in the strong coupling regime through collective multimode behavior.
• Disorder Effects: The introduction of disorder (inhomogeneities in qubit frequencies or coupling strengths) shifts the phase boundary. The authors derive analytical conditions showing how the critical coupling values change in the presence of disorder, generally lowering the threshold for the transition in specific configurations.

The authors illustrate these general features using three specific models:

1. Multimode Dicke Model (One- and Two-Photon Terms): The analysis confirms that the phase boundary follows $|\vec{\lambda}| = 1$ (in dimensionless units). The presence of two-photon terms ($\gamma' \neq 0$) leads to a first-order transition, whereas the pure one-photon case yields a second-order transition. The multimode nature reduces the critical coupling strength per mode.
2. Rabi-Stark-Hubbard Model: This model includes Stark shifts and inter-cavity hopping. The phase boundary is found to be $\gamma^2 + \tilde{J} + \tilde{U} = 1$, where $\gamma$ is the coupling, $\tilde{J}$ is the hopping, and $\tilde{U}$ is the Stark shift. The transition is second-order. The results align with self-consistent numerical methods in the classical oscillator limit.
3. Anisotropic Rabi-Stark Model: This model exhibits two distinct SPs (u-type and v-type) corresponding to different symmetry breakings. The transition from NP to either SP is second-order, while the transition between the two SPs is first-order. Despite the complexity of multiple SPs, the system still enters the SP region radially from the NP in the coupling parameter space.

Significance and Claims
The paper claims to provide a unified framework for understanding SPTs in general light-matter systems. By identifying the universal feature that the origin is always in the NP and the transition occurs radially outward in coupling space, the authors offer a simplified method to calculate phase boundaries and determine the order of transitions without needing to solve complex differential equations for every specific model.

The significance lies in:

• Simplification: The ability to determine SPT properties (existence, order, boundary) via a concise method based on the global minimum of a Landau potential.
• Experimental Relevance: The finding that multimode collective behavior can lower the coupling strength requirements for SPT suggests that these transitions might be achievable in regimes previously considered difficult (e.g., strong coupling rather than ultrastrong coupling).
• Robustness: The demonstration that disorder shifts but does not necessarily destroy the phase transition, providing insight into the stability of SPTs in realistic, imperfect systems.

The authors conclude that these general features facilitate the study of SPTs and their diverse applications, particularly in quantum information processing, by offering a predictive tool for the phase behavior of complex light-matter Hamiltonians.