Coexistence Regime and Thermal Crystallization in the cavity-mediated extended Bose-Hubbard Model

Using path integral Monte Carlo simulations, this study reveals that in the cavity-mediated extended Bose-Hubbard model, the coexistence regime between superfluid and charge-density-wave phases exhibits strong metastability and a counterintuitive thermally induced crystallization when starting from a superfluid state, contrasting with the smooth melting observed from a charge-density-wave initialization.

The Gist

Imagine a crowded dance floor where thousands of tiny, invisible dancers (atoms) are trying to move in perfect sync. This paper is about a special kind of dance floor where the dancers can influence each other not just by touching, but by shouting across the room through a giant, magical megaphone (the optical cavity).

The scientists wanted to understand how these dancers behave when the room gets hot. They discovered some very surprising things about how the dancers switch between different "modes" of dancing, and how the history of the dance (where they started) changes the outcome.

Here is the story of their discovery, broken down into simple concepts:

1. The Two Main Dance Styles

At the very bottom of the temperature scale (absolute zero), the dancers can only do two main things:

• The Superfluid (The Smooth Flow): Everyone holds hands and moves together in a giant, fluid wave. No one bumps into anyone; it's a perfect, frictionless glide.
• The Crystal (The Checkerboard): Everyone freezes into a rigid grid, like soldiers standing in perfect rows and columns. They are stuck in place, alternating between "full" spots and "empty" spots.

Usually, you expect a system to be one or the other. But in this specific setup, the scientists found a Battle Zone.

2. The Battle Zone (Coexistence Regime)

In a specific area of their experiment, the "Smooth Flow" and the "Crystal" are equally strong. It's like a tug-of-war where both teams are pulling with exactly the same force.

• If you start the simulation with the dancers in a Smooth Flow, they stay in a Smooth Flow.
• If you start them in a Crystal, they stay in a Crystal.
• Neither side can easily knock the other over. This is called metastability—a state that is stuck in place, waiting for a nudge.

3. The Weird Heat Effect: "Thermal Crystallization"

This is the most surprising part. Usually, when you heat something up, it melts. Ice turns to water; a solid turns to liquid. You expect heat to destroy order.

But here, the scientists found a weird trick:

• Scenario A (Starting from Smooth Flow): They heated up the Smooth Flow dancers. First, the dancers got confused and stopped holding hands (becoming a "Normal Fluid"). But then, as they got even hotter, something magical happened. The heat actually helped them organize into a Crystal!
  • Analogy: Imagine a chaotic crowd of people running around. If you turn up the music (heat) just right, they suddenly stop running and start marching in perfect lockstep. The heat forced them to organize.
• Scenario B (Starting from Crystal): They heated up the Crystal dancers. They just slowly melted into a chaotic "Normal Fluid." They never went back to being a Smooth Flow.

4. Why Does This Happen?

The "megaphone" (the cavity) is the key. It connects every dancer to every other dancer instantly.

• When the system is in the Smooth Flow, the dancers are constantly swapping places (exchanging positions).
• The heat makes it harder for them to swap places.
• Because the megaphone makes them care about the pattern of who is where, when the swapping stops due to heat, the dancers get "stuck" in the rigid Crystal pattern. The heat actually helped them lock into the crystal shape because it stopped them from messing up the pattern by swapping around.

5. The Takeaway

The paper teaches us two main lessons:

1. History Matters: In this "Battle Zone," what happens next depends entirely on where you started. If you start smooth, you might get a crystal later. If you start as a crystal, you just melt.
2. Heat Can Build Order: We usually think heat destroys order. But in this specific quantum dance, heat actually helped build a crystal structure by stopping the dancers from swapping places.

In summary: The scientists found a strange quantum dance floor where heating up a fluid can accidentally turn it into a solid crystal, but only if you started with the fluid. If you started with the crystal, it just melts. It's a reminder that in the quantum world, the path you take to get somewhere is just as important as the destination.

Technical Summary

1. Problem Statement

The paper investigates the finite-temperature behavior of the extended Bose-Hubbard model on a two-dimensional square lattice, specifically focusing on systems with cavity-mediated long-range interactions at unit filling ($\langle n \rangle = 1$).

While previous studies have explored the thermal melting of supersolid phases (where diagonal and off-diagonal orders coexist), this work addresses a gap in understanding the first-order phase transition regime between the Superfluid (SF) and Charge-Density-Wave (CDW) phases. Specifically, the authors aim to:

• Characterize the ground-state phase diagram, particularly the width and nature of the SF-CDW coexistence region.
• Analyze how competing SF and CDW orders evolve under thermal fluctuations.
• Determine if metastability persists at finite temperatures and if thermal effects can induce crystallization (thermocrystallization) in a system initialized as a superfluid.

2. Methodology

The authors employ Path-Integral Quantum Monte Carlo (PIMC) simulations based on the worm algorithm. Key methodological aspects include:

• Model Hamiltonian: An extended Bose-Hubbard model including nearest-neighbor hopping ($t$), on-site repulsion ($U_s$), chemical potential ($\mu$), and a global cavity-mediated interaction ($U_\ell$) that favors particle number imbalance between even and odd sublattices.
• Order Parameters:
  • Superfluid Density ($\rho_s$): Calculated via winding number fluctuations (Pollock–Ceperley formula).
  • Compressibility ($\kappa$): Derived from particle number fluctuations.
  • Structure Factor ($S(\pi, \pi)$): Used to detect CDW order (checkerboard pattern).
• Finite-Size Scaling: Simulations were performed on lattice sizes ranging from $L=10$ to $L=30$. Critical points were determined using data collapse techniques and universal scaling functions for Ising and XY universality classes.
• Hysteresis Protocols: To identify first-order transitions and coexistence regions, the authors performed "forward sweeps" (starting from SF and increasing interaction) and "backward sweeps" (starting from CDW and decreasing interaction) to observe path dependence and metastability.
• Thermal Trajectories: Simulations were initialized in specific states (pure SF or pure CDW) and heated to observe the evolution of order parameters.

3. Key Contributions and Results

A. Ground-State Phase Diagram Refinement

• Phases Identified: The study confirms four stable phases: Superfluid (SF), Mott Insulator (MI), Supersolid (SS), and Charge-Density-Wave (CDW).
• Transition Nature:
  • SF $\leftrightarrow$ SS and CDW $\leftrightarrow$ SS transitions are continuous (second-order), belonging to the $(2+1)$-dimensional Ising and XY universality classes, respectively.
  • SF $\leftrightarrow$ CDW and MI $\leftrightarrow$ CDW transitions are strongly first-order.
• Broad Coexistence Region: A major finding is the discovery of a significantly broader coexistence region between SF and CDW phases than previously reported. This region is characterized by strong metastability where the final state depends entirely on the initial microscopic configuration.

B. Thermal Evolution and "Thermocrystallization"

The most significant discovery concerns the behavior of the system within the first-order coexistence region as temperature ($T$) increases:

• SF-Initialized Trajectory: When starting from a Superfluid state:
  1. $\rho_s$ vanishes at low $T$ ($T/t \sim 1.0$), transitioning to a Normal Fluid (NF).
  2. Upon further heating ($T/t \approx 3.3$), a thermally induced crystalline order emerges. The structure factor $S(\pi, \pi)$ peaks, indicating a stable CDW phase. This phenomenon is termed thermocrystallization.
  3. At very high $T$ ($T/t \sim 5.5$), the CDW order melts into the Normal Fluid.
  • Sequence: $SF \to NF \to CDW \to NF$.
• CDW-Initialized Trajectory: When starting from a CDW state:
  1. The system remains in the CDW phase as temperature increases.
  2. No reentrant superfluidity is observed.
  3. The CDW order melts directly into the Normal Fluid at high $T$.
  • Sequence: $CDW \to NF$.
• Implication: This path dependence confirms that at intermediate temperatures, the CDW phase represents the global free-energy minimum, while the SF phase is metastable. The emergence of CDW order upon heating a superfluid is driven by thermal suppression of particle exchange cycles, which favors solid ordering.

C. Finite-Temperature Phase Diagram

• The authors constructed a $(U_\ell/U_s, T/t)$ phase diagram.
• Metastability Window: A gray-shaded region exists where both SF and CDW (or NF and CDW) are locally stable. The realized phase depends on initial conditions.
• Supersolid Melting: In the supersolid regime, thermal melting occurs in two steps: either phase coherence is lost first (KT transition) followed by density modulation loss, or vice versa, depending on the interaction strength.
• Robustness: CDW order induced by strong cavity interactions is highly robust against thermal fluctuations, with melting temperatures increasing rapidly with interaction strength.

4. Significance

• Theoretical Insight: The paper challenges the notion that thermal fluctuations only melt ordered states. It demonstrates that in systems with cavity-mediated long-range interactions, thermal fluctuations can actually induce crystallization (thermocrystallization) by suppressing quantum exchange processes that favor superfluidity.
• Experimental Relevance: The findings are directly applicable to cavity-QED experiments with ultracold atoms. The ability to tune cavity-mediated interactions and control heating rates allows experimentalists to access these metastable domains and observe thermally assisted ordering processes.
• Resolution of Discrepancies: The work resolves discrepancies with earlier studies (e.g., Ref. 12) by revealing that the SF-CDW transition is not a sharp line but a broad coexistence region, fundamentally changing the understanding of the phase diagram's topology.
• Metastability: It establishes that metastability in these quantum systems persists up to finite temperatures, creating a complex landscape where the history of the system (initial state) dictates its thermodynamic path.

In conclusion, this study provides a comprehensive map of the quantum and thermal phases of cavity-mediated bosons, highlighting a unique mechanism where heating a superfluid can paradoxically stabilize a crystalline charge-density-wave state before eventual thermal melting.