Order-by-disorder and emergent Kosterlitz-Thouless phase in triangular Rydberg array

This study employs numerically exact quantum Monte Carlo simulations to reveal that triangular Rydberg atom arrays exhibit $\sqrt{3}\times\sqrt{3}$ anti-ferromagnetic order at specific fillings and an emergent Kosterlitz-Thouless phase driven by an order-by-disorder mechanism at half-filling, providing theoretical predictions for future experimental verification.

The Gist

The Big Picture: A Quantum Dance Floor

Imagine a giant, perfectly round dance floor made of a triangular grid (like a honeycomb). On this floor, we have a bunch of dancers (atoms). These aren't normal dancers; they are Rydberg atoms. Think of them as "super-dancers" that can be in two states:

1. Ground State: They are sitting quietly on the floor.
2. Rydberg State: They are jumping up high, but if two of them jump too close to each other, they repel each other violently (like magnets with the same pole facing each other).

Scientists can control how many dancers jump and how hard they push each other. This setup is called a Quantum Simulator. It's like a video game where the physics engine is real, allowing us to watch how these atoms behave in ways that are impossible to calculate with a normal computer.

The Problem: The "Frustrated" Triangle

The dance floor is a triangular lattice. This creates a problem known as geometric frustration.

• The Analogy: Imagine three friends (atoms) standing at the corners of a triangle. They all want to be "different" from their neighbors (like one sitting, one jumping, one sitting).
  • If Friend A sits and Friend B jumps, Friend C is stuck. If C sits, they are the same as A. If C jumps, they are the same as B.
  • There is no perfect arrangement where everyone is happy. This is "frustration."

Usually, when things are frustrated, they just get messy and chaotic. But this paper asks: Can chaos actually create order?

The Discovery 1: The "Order-by-Disorder" Magic

The researchers used a super-advanced computer simulation (Quantum Monte Carlo) to watch what happens at different "filling levels" (how many atoms are jumping).

1. The Easy Wins (1/3 and 2/3 Full):
When the floor is 1/3 or 2/3 full, the atoms naturally fall into a neat, repeating pattern (like a checkerboard but triangular). This is what experiments had already seen. It's like a well-organized marching band.

2. The Magic Trick (1/2 Full):
Here is the exciting part. When the floor is exactly half-full, the atoms should be a mess because of the frustration. There are millions of ways they could arrange themselves, and they are all equally "happy" (energetically).

However, the simulation revealed a miracle: Order-by-Disorder (OBD).

• The Analogy: Imagine a room full of people trying to decide where to stand. There are 100 spots, and everyone is equally happy in any spot. But, if they all start shuffling around randomly (thermal fluctuation), they might accidentally find a specific pattern that lets them move the most freely without bumping into each other.
• The Result: The "noise" and "chaos" of the atoms actually force them into a specific, beautiful pattern ($\sqrt{3} \times \sqrt{3}$). The disorder creates the order. This is a rare and beautiful phenomenon in physics.

The Discovery 2: The "Liquid Crystal" Phase (Kosterlitz-Thouless)

The paper also looked at what happens when you heat up the dance floor (increase the temperature).

At 1/3 and 2/3 full:
As you heat it up, the neat marching band breaks down quickly and turns into a chaotic crowd. This is a standard "phase transition."

At 1/2 full (The Special Case):
Because of the "Order-by-Disorder" magic, something weird happens as you heat it up.

• The Analogy: Imagine a group of dancers who are holding hands in a circle. As the music gets faster (temperature rises), they don't just break apart immediately. Instead, they enter a "liquid crystal" state. They aren't frozen in a rigid grid, but they aren't a total mess either. They are swirling in a fluid, continuous circle.
• The Science: In physics terms, the atoms gain U(1) symmetry. This means they can rotate their "phase" (their position in the cycle) freely, like a spinning top that can point in any direction, rather than being stuck pointing North, East, or South.
• The Transition: Eventually, if it gets hot enough, this fluid swirl breaks apart into total chaos. The transition from the "swirling fluid" to "total chaos" is called a Kosterlitz-Thouless (KT) transition. It's a very special type of change that only happens in 2D systems and is famous in physics.

Why Does This Matter?

1. It's Real: The researchers didn't just guess; they used a "numerically exact" method to prove this happens.
2. It's Predictable: They predict that if experimentalists tweak their Rydberg atom arrays to be exactly half-full, they will see this "Order-by-Disorder" pattern and the special "swirling" phase transition.
3. New Physics: It shows that in the quantum world, chaos isn't always the enemy. Sometimes, the noise is the glue that holds a new, exotic state of matter together.

Summary in a Nutshell

• The Setup: A triangular grid of atoms that hate being close to each other.
• The Surprise: At half-capacity, the atoms' natural chaos forces them into a perfect, repeating pattern (Order-by-Disorder).
• The Twist: When heated, this pattern doesn't just break; it turns into a fluid, spinning state (Emergent U(1) symmetry) before finally dissolving.
• The Takeaway: Nature is full of surprises where disorder creates order, and we can now simulate and predict these exotic dances of atoms.

Technical Summary

1. Problem Statement

The paper addresses the challenge of deciphering exotic physics in frustrated quantum magnets, specifically those realized on triangular lattices. While recent experiments using Rydberg atom arrays have successfully observed classical $\sqrt{3} \times \sqrt{3}$ antiferromagnetic (AFM) order at specific fillings, several theoretically predicted quantum phenomena remain elusive in experimental settings. Key open questions include:

• Can the Order-by-Disorder (OBD) mechanism, where quantum or thermal fluctuations select a specific ordered state from a macroscopically degenerate manifold, be realized in Rydberg arrays at 1/2-filling?
• Does an emergent U(1) symmetry arise at 1/2-filling, leading to a Kosterlitz-Thouless (KT) phase transition, as opposed to the standard discrete symmetry breaking seen at other fillings?
• A rigorous, numerically exact theoretical framework is needed to validate the plausibility of these phenomena before they are targeted in future experiments.

2. Methodology

The authors employ numerically exact Quantum Monte Carlo (QMC) simulations based on the Stochastic Series Expansion (SSE) algorithm.

• Model: They simulate a 2D interacting Rydberg atom array on a triangular lattice described by the Hamiltonian:
  $$\hat{H} = \sum_{i<j} U_{ij}\hat{n}_i\hat{n}_j + \frac{\hbar\Omega}{2}\sum_i (\hat{b}^\dagger_i + \hat{b}_i) - \hbar\delta\sum_i \hat{n}_i$$
  where $U_{ij}$ represents long-range van der Waals interactions ($1/r^6$), $\Omega$ is the Rabi frequency, and $\delta$ is the detuning.
• Algorithmic Innovation: To handle large system sizes ($L \times L$ up to $L=18$) and low temperatures efficiently, they implement a combination of line updates and multi-branch updates.
• Advantage: The model is free of the sign problem, allowing for precise calculations of ground-state and finite-temperature properties without statistical bias.
• Analysis: They utilize finite-size scaling analysis, calculating order parameters ($m$), static structure factors ($S(\vec{Q})$), Binder ratios ($U_B$), and correlation length ratios ($C$) to determine phase boundaries and universality classes.

3. Key Contributions and Results

A. Ground State Phase Diagram

The authors map the ground state phase diagram as a function of detuning ($\delta$) and Rydberg filling ($\bar{n}$):

• 1/3 and 2/3 Filling: They confirm the existence of a robust $\sqrt{3} \times \sqrt{3}$ Triangular Antiferromagnetic (TAF) long-range order. This result is consistent with recent experimental observations (e.g., Scholl et al., Nature 2021).
• Phase Transitions: The transitions between the disordered phase and the TAF phases at 1/3 and 2/3 fillings are identified as first-order transitions, occurring at critical detuning values $\delta_c \approx 1.2$ and $\delta_c \approx 17.8$.

B. Order-by-Disorder (OBD) at 1/2 Filling

The most significant finding is the discovery of a $\sqrt{3} \times \sqrt{3}$ TAF ordered phase at 1/2-filling, which is driven by the Order-by-Disorder mechanism.

• Unlike the 1/3 and 2/3 fillings where the order is classical, the order at 1/2-filling emerges from the selection of a specific ground state from a macroscopically degenerate manifold due to quantum fluctuations.
• The simulation confirms the existence of finite order parameters and structure factors at thermodynamic limits for this filling.

C. Emergent U(1) Symmetry and KT Phase Transition

The study reveals a profound difference in symmetry breaking between the 1/3/2/3 fillings and the 1/2 filling at finite temperatures:

• 1/3 and 2/3 Fillings: The system breaks $Z_3$ symmetry. The transition to the disordered phase belongs to the 3-state Potts universality class.
• 1/2 Filling: The system possesses an extra $Z_2$ particle-hole symmetry, resulting in a six-fold degenerate ground state.
  • Irrelevance of Anisotropy: The authors demonstrate that the $Z_6$ anisotropic term (which would break continuous symmetry) is irrelevant in the presence of strong thermal fluctuations.
  • Emergent U(1): Consequently, the system exhibits an emergent continuous U(1) symmetry over a wide temperature range.
  • KT Transition: The transition from this quasi-long-range ordered phase to the high-temperature disordered phase belongs to the Kosterlitz-Thouless (KT) universality class.
  • Evidence:
    • The Binder ratio crossing points are not distinct (characteristic of KT transitions).
    • The six-fold anisotropy parameter $C_6$ vanishes as system size increases.
    • Histograms of the order parameter show an isotropic distribution in the complex plane (U(1) symmetric) rather than discrete peaks.
    • Data collapse analysis of susceptibility confirms the scaling behavior consistent with KT theory ($\eta = 0.25$).

4. Significance

• Theoretical Validation: This work provides the first numerically exact confirmation that Rydberg atom arrays on triangular lattices can host OBD-driven phases and emergent U(1) symmetries.
• Experimental Guidance: The results offer concrete theoretical support and specific targets (e.g., specific filling factors and temperature regimes) for future experiments to observe the elusive KT phase transition and OBD mechanisms in quantum simulators.
• Platform Potential: It reinforces the Rydberg atom array as a premier platform for simulating complex frustrated magnetism, capable of accessing phases that are difficult to realize in natural materials.
• Methodological Benchmark: The successful application of SSE-QMC with multi-branch updates to large triangular lattice systems sets a benchmark for studying other frustrated quantum systems without sign problems.

In summary, the paper bridges the gap between theoretical predictions of exotic frustrated magnetism and experimental realization, demonstrating that triangular Rydberg arrays can host a rich phase diagram including OBD-driven order and topological KT transitions.