Supersolid phase in two-dimensional soft-core bosons at finite temperature

This study investigates the finite-temperature phase diagram of two-dimensional soft-core bosons using self-consistent Hartree-Fock and quantum Monte Carlo methods, identifying a broad supersolid phase and a potential intermediate hexatic phase while validating mean-field theory as an effective tool for analyzing these transitions.

The Gist

Imagine a magical ballroom filled with thousands of tiny dancers (bosons). Usually, these dancers follow one of two rules:

1. The Superfluid: They all hold hands in a giant, invisible circle, moving in perfect unison without any friction. They flow like water, but they are all one single quantum entity.
2. The Solid: They freeze into a rigid grid, like soldiers standing in perfect rows and columns. They can't flow; they are stuck in place.

But what if they could do both at the same time? What if they formed a rigid crystal grid, yet somehow, they could still flow through each other without friction?

This impossible state is called a Supersolid. It's like a frozen ice cube that can also pour like water.

This paper investigates how this "Supersolid" behaves in a flat, two-dimensional world (like a sheet of paper) when you start to heat it up. The researchers used two different methods to solve this puzzle: a fast, clever mathematical shortcut (Hartree-Fock) and a slow, super-accurate computer simulation (Quantum Monte Carlo).

Here is the story of their findings, explained with some everyday analogies:

1. The Two Maps

The researchers created two "maps" of the ballroom to see what happens as the temperature changes and the dancers get more energetic.

• Map A (The Shortcut): This is like using a weather forecast based on average trends. It's fast and gives a good general idea. It predicted that the Supersolid exists in a wide, comfortable zone between the flowing liquid and the frozen solid.
• Map B (The Satellite): This is like sending a drone to take high-resolution photos of every single dancer. It's computationally expensive (takes a lot of power) but very precise. It confirmed the general shape of Map A but showed that the "Supersolid" zone is actually a bit smaller and the transitions happen at lower temperatures than the shortcut predicted.

2. The "Hexatic" Mystery

In the world of 2D physics, there is a weird, intermediate state called the Hexatic phase.

• Imagine a crowd of people:
  • Solid: Everyone is standing in a perfect grid, facing the same direction.
  • Fluid: Everyone is running around randomly, facing every which way.
  • Hexatic: The people are still running around randomly (no grid), but they are all trying to face the same direction as their neighbors. They have "orientation" but no "position."

The researchers looked for this Hexatic phase. The shortcut map said, "Yes, there's a big room for this!" The high-precision simulation said, "Maybe, but it's a very tiny, narrow hallway between the solid and the fluid." It's so thin that it's hard to see, and they couldn't find a "Hexatic Superfluid" (a flowing Hexatic phase) in this specific system.

3. The "Order-by-Disorder" Surprise

Here is the most mind-bending part of the paper. Usually, when you heat something up, it gets more chaotic. Ice melts into water; a neat stack of papers blows away in the wind.

But in this quantum ballroom, the researchers found a strange phenomenon: Heating the system actually made the dancers stand up straighter.

• The Analogy: Imagine a group of people trying to stand in a line on a wobbly boat. If the boat is perfectly still (very cold), they might be stiff and rigid. But if the boat starts rocking gently (adding a little heat), they might instinctively grab onto each other and align their bodies to stay balanced, creating a more organized structure than before.
• The Result: As they warmed up the Supersolid, the "orientational order" (how well the dancers faced the same way) actually increased for a while before finally melting into chaos. This is called "Order-by-Disorder," where thermal energy (disorder) accidentally helps create order.

4. Why This Matters

This paper is important for a few reasons:

• Validating the Shortcut: They proved that the fast mathematical method (Hartree-Fock) is actually a very good tool for predicting these complex quantum states. This means scientists can use this "shortcut" to explore new materials without needing to run massive, expensive supercomputer simulations every time.
• Understanding Quantum Melting: They showed exactly how a quantum crystal melts. It doesn't just turn into a liquid; it might pass through that weird Hexatic phase, and the transition involves a jump in how the material flows.
• Real-World Applications: While this is theoretical, it helps us understand real experiments with ultra-cold atoms (dipolar gases) that scientists are creating in labs right now. It tells us what to look for when we try to make a Supersolid in the real world.

The Takeaway

The universe is full of surprises. Sometimes, adding heat doesn't just break things down; it can rearrange them into new, strange, and beautiful patterns. This paper maps out the territory of a "Supersolid," showing us that in the quantum world, you can be a rigid crystal and a flowing river at the same time, and sometimes, a little bit of warmth helps you hold your shape even better.

Technical Summary

1. Problem Statement

The paper investigates the finite-temperature phase diagram of a two-dimensional (2D) system of bosons interacting via a soft-core repulsive potential. While the existence of supersolidity (simultaneous breaking of $U(1)$ gauge symmetry and translational symmetry) has been established in 3D dipolar gases and at zero temperature in 2D models, the behavior of these systems at finite temperatures remains complex.

Key challenges addressed include:

• 2D Fluctuations: In 2D, true long-range order (LRO) is forbidden by thermal fluctuations (Mermin-Wagner theorem). Instead, phases are characterized by quasi-long-range order (power-law decay of correlations) and Berezinskii-Kosterlitz-Thouless (BKT) transitions.
• Melting Scenarios: The melting of a 2D crystal can follow different pathways, potentially involving an intermediate hexatic phase (short-range positional order, quasi-long-range orientational order) as predicted by KTHNY theory, or a direct first-order transition.
• Quantum vs. Thermal Effects: It is unclear how quantum exchange effects influence the freezing/melting transition and whether a hexatic superfluid phase (superfluidity coexisting with hexatic order) can exist.
• Methodological Gap: Existing studies often rely on computationally expensive Quantum Monte Carlo (QMC) or simplified mean-field approximations (like Local Density Approximation) that may fail to capture the full self-consistent nature of the condensate and thermal excitations.

2. Methodology

The authors employ a dual-method approach to construct and validate the phase diagram:

A. Self-Consistent Hartree-Fock (HF) Theory

• Approach: They develop a fully self-consistent HF theory for bosons at finite temperature beyond the Local Density Approximation (LDA).
• Formalism:
  • Based on a finite-temperature variational principle using the Bogoliubov inequality.
  • Solves the extended Gross-Pitaevskii (GP) equation coupled with a self-consistency equation for the Hartree-Fock potential.
  • Treats the condensate wavefunction ($\psi$) and thermal excitations (fluctuation operators) self-consistently.
• Calculations:
  • Determines the free energy of uniform (fluid) and spatially modulated (solid/supersolid) phases to identify first-order transitions.
  • Calculates the superfluid density ($\rho_s$) using the phase-twist method (twisted boundary conditions).
  • Computes the Young modulus ($Y$) by applying shear and hydrostatic deformations to the lattice.
  • Uses universal BKT relations to estimate critical temperatures for superfluid and melting transitions:
    • Superfluid transition: $\frac{\hbar^2 \rho_s}{m} = \frac{2 k_B T_{BKT}}{\pi}$
    • Solid-Hexatic transition: $Y a^2 = 16 \pi k_B T_{BKT}$

B. Path-Integral Monte Carlo (PIMC) Simulations

• Approach: Exact numerical simulations of the microscopic Hamiltonian using the continuous-space PIMC algorithm with periodic boundary conditions.
• Scale: Simulations performed with large system sizes ($N = 4032$ and $N = 16128$ particles) to minimize finite-size effects.
• Observables:
  • Superfluid Fraction: Calculated via the winding number estimator.
  • Orientational Order: Defined by the six-fold director field $\Psi_6$ and the global order parameter $|\psi_6|^2$.
  • Correlation Functions: Analysis of density-density ($g_2$) and orientational ($G_6$) correlations to distinguish between fluid, solid, supersolid, and hexatic phases based on decay laws (exponential vs. power-law).

3. Key Contributions

1. Beyond-LDA HF Theory: Demonstrated the feasibility and accuracy of a self-consistent HF scheme for 2D bosons without LDA, showing it qualitatively reproduces complex phase diagrams and serves as an efficient tool for exploring parameter spaces.
2. Finite-Temperature Phase Diagram: Mapped the complete phase diagram of 2D soft-core bosons, identifying regions for Normal Fluid, Superfluid, Normal Solid, Supersolid, and potential Hexatic phases.
3. Hexatic Superfluid Investigation: Systematically investigated the existence of a hexatic superfluid phase. While HF theory predicts a hexatic phase separating the fluid and solid, PIMC results suggest this region is extremely narrow or potentially non-existent as a distinct superfluid phase.
4. Counter-Intuitive Ordering: Discovered a regime where orientational order increases with temperature within the solid phase before melting, attributed to enhanced quantum fluctuations (an "order-by-disorder" mechanism).

4. Key Results

Phase Diagram Structure

• Supersolid Region: A broad region exists at low temperatures and high interaction strengths ($W$) where the system exhibits both superfluidity and a triangular density modulation (supersolid).
• Transitions:
  • Superfluid $\to$ Normal Fluid: A BKT transition characterized by a universal jump in superfluid density.
  • Supersolid $\to$ Normal Solid: PIMC simulations indicate a discontinuous jump in the orientational order parameter $|\psi_6|^2$ and superfluid density, suggesting a first-order transition or a very sharp BKT transition.
  • Solid $\to$ Fluid: The melting transition is complex. HF theory predicts a hexatic phase between the solid and fluid. PIMC simulations show a narrow "unknown" region where correlation functions decay slowly, hinting at a hexatic phase, but finite-size effects make a definitive conclusion difficult.

Comparison of Methods

• HF vs. PIMC: The HF phase diagram is in qualitative agreement with PIMC but overestimates the stability regions of the superfluid and supersolid phases (predicting transitions at higher temperatures).
• Critical Points: HF predicts the onset of supersolidity at interaction strengths consistent with PIMC. However, HF fails to capture the BKT suppression of superfluidity, predicting finite $\rho_s$ above the critical temperature.

Anomalous Behavior

• Temperature-Induced Ordering: In the solid phase at intermediate interaction strengths, the orientational order parameter $|\psi_6|^2$ increases as temperature rises, reaches a maximum, and then drops as the system melts. This contradicts classical intuition (where thermal disorder reduces order) and is attributed to the interplay between thermal fluctuations and quantum exchange effects.

Hexatic Phase

• HF theory predicts a stable hexatic phase separating the normal fluid and normal solid.
• PIMC simulations find a very thin region of "unknown" points between the fluid and solid where $G_6(r)$ decays algebraically. However, the existence of a hexatic superfluid (simultaneous superfluidity and hexatic order) is not clearly supported; the region is likely too narrow to be distinct or is absent.

5. Significance

• Theoretical Tool Validation: The study establishes self-consistent HF theory (without LDA) as a computationally efficient and reliable method for exploring the phase diagrams of quantum many-body systems, complementing expensive QMC simulations.
• Understanding 2D Supersolids: It provides a microscopic characterization of the melting of supersolids in 2D, clarifying the role of thermal fluctuations and the nature of the transitions (BKT vs. First Order).
• Quantum Fluctuations: The identification of temperature-induced ordering highlights the non-trivial role of quantum fluctuations in 2D quantum crystals, a phenomenon not captured by classical theories.
• Experimental Relevance: The findings are relevant for ongoing experiments with ultracold dipolar gases and 2D quantum materials, guiding the search for hexatic superfluids and the characterization of finite-temperature supersolid phases.

In conclusion, the paper successfully bridges the gap between mean-field approximations and exact simulations, providing a comprehensive view of the rich phase behavior of 2D soft-core bosons and revealing subtle quantum effects in the melting of supersolids.