Scale invariance of the polaron energy at the Mott-superfluid critical point

This paper demonstrates through ground-state quantum Monte Carlo calculations that the energy of a mobile impurity in a lattice Bose gas exhibits scale invariance at the Mott-superfluid critical point, establishing impurity spectroscopy as a viable method for probing the critical properties of quantum phase transitions.

The Gist

Imagine a crowded dance floor where everyone is moving in perfect, synchronized harmony. This is a superfluid—a state of matter where particles flow without any friction, like a super-cooled liquid that never stops moving. Now, imagine a few dancers suddenly stop moving and lock into rigid, individual spots, refusing to dance with the crowd. This is a Mott insulator—a state where particles are stuck in place.

The moment the music changes and the dancers switch from "stuck" to "flowing" (or vice versa) is called a Quantum Phase Transition. It's a dramatic shift in the rules of the game, driven not by temperature, but by the quantum nature of the particles themselves.

The Problem: Finding the "Switch"

Scientists have long known how to describe these transitions, but finding the exact moment the switch flips is incredibly hard. Usually, you have to measure complex, invisible patterns (like how every dancer is connected to every other dancer) to know where the transition is. It's like trying to find the exact second a crowd turns into a riot by listening to the entire stadium at once—it's noisy and difficult.

The Solution: The "Impurity" Detective

This paper introduces a clever trick: The Polaron.

Imagine dropping a single, slightly different dancer (an impurity) onto the floor.

• If the crowd is stuck (insulator), this new dancer bumps into rigid walls and feels heavy.
• If the crowd is flowing (superfluid), the dancer gets swept up in the current, dragging a little cloud of excited neighbors with them. This new, heavier package is called a Polaron.

The researchers asked: Can we just measure how heavy this "Polaron dancer" feels to figure out exactly when the crowd switches from stuck to flowing?

The Discovery: The Magic Moment

Using powerful supercomputers (simulating the dance floor with millions of calculations), the team dropped this impurity into a 2D grid of bosons (the dancers) and watched what happened as they changed the "music" (the hopping strength).

They found something magical:

1. The Scale-Invariant Energy: At the exact moment the transition happens, the energy (or "heaviness") of the impurity becomes scale-invariant.
   • The Analogy: Imagine looking at a fractal (like a fern leaf). No matter how much you zoom in or out, the pattern looks the same. Similarly, at the critical point, the impurity's energy looks the same whether you are looking at a tiny dance floor (small system) or a massive stadium (large system). It doesn't matter how big the room is; the "feeling" of the transition is identical.
2. The Flattening Effect: They also looked at how the impurity interacts with dancers far away.
   • The Analogy: Usually, if you push a dancer, the ripple effect fades quickly. But right at the transition point, the ripple effect becomes so strong and long-lasting that the interaction looks "flat" and uniform across the whole floor. It's as if the whole dance floor holds its breath at the same time.

The Surprise: A New Mystery

The team expected the impurity to perfectly mimic the crowd's behavior. However, they found a twist. The mathematical "speed" at which the impurity's energy changes near the transition (a number called the exponent) was different from the speed at which the crowd itself changes.

• The Analogy: It's like the crowd is running a race at a specific speed, but the impurity is running alongside them at a slightly different, unexplained speed. Theoreticians don't yet have a rulebook that explains why the impurity behaves this way. It's a new puzzle for physics to solve.

Why This Matters

This research is a game-changer for two reasons:

1. Easier Measurement: Instead of trying to measure the complex, invisible connections of the whole crowd, scientists can now just measure the "weight" of a single impurity. It's a much simpler, more direct way to find the critical point.
2. New Tools for Quantum Tech: As we build better quantum computers and simulators using cold atoms, having a simple "probe" (the impurity) to tell us exactly what state the system is in is invaluable. It's like having a simple thermometer that tells you exactly when water turns to ice, without needing to analyze the molecular structure of the ice.

In short: The paper shows that by dropping a single "spy" into a quantum system, we can easily spot the exact moment the system changes its fundamental nature, revealing a hidden, scale-invariant beauty at the heart of quantum transitions.

Technical Summary

1. Problem Statement

Quantum phase transitions (QPTs) are typically characterized by an order parameter and long-range correlation functions. However, these quantities are often difficult to access experimentally or via direct numerical simulation due to the need for high-resolution, phase-sensitive measurements.

• The Challenge: Identifying and measuring the order parameter for specific transitions (like the Mott insulator–superfluid transition) is nontrivial.
• The Proposed Solution: The authors investigate whether the energy of a mobile impurity (a polaron) embedded in a many-body system can serve as a sensitive, accessible probe of critical behavior. Specifically, they ask if the polaron energy exhibits scale invariance at the critical point of the Bose-Hubbard model, similar to the host system's bulk properties.

2. Methodology

The study employs Ground-State Full Configuration Interaction Quantum Monte Carlo (FCIQMC) to simulate a two-dimensional Bose-Hubbard model with a single impurity.

• Model System:
  • Host: Bosons on a 2D square lattice ($L \times L$) with periodic boundaries, described by the Bose-Hubbard Hamiltonian $\hat{H}_B$ with hopping $t$ and repulsive interaction $U$. The system is at unit filling ($n_B = 1$).
  • Impurity: A single mobile impurity is added via Hamiltonian $\hat{H}_I$, interacting with host bosons with strength $U_{IB}$.
  • Total Hamiltonian: $\hat{H} = \hat{H}_B + \hat{H}_I + U_{IB} \sum_i \hat{n}_{i,I}\hat{n}_{i,B}$.
• Computational Approach:
  • FCIQMC: A projector Monte Carlo method that iterates a stochastic vector to converge to the ground state eigenvector.
  • Bias Control: Importance sampling and population control techniques are used to minimize statistical bias, allowing simulations of Hilbert space dimensions up to $\approx 10^{60}$.
  • System Sizes: Simulations were performed for lattice sizes $L$ ranging from 5 to 10 for the scaling analysis, and exact diagonalization was used for smaller systems ($L=3, 4$) to study correlation functions.
• Scaling Analysis:
  • The polaron energy $E_p$ is defined as the energy difference between the system with and without the impurity: $E_p = E(U_{IB}) - E(0)$.
  • A finite-size scaling ansatz was applied:
    $$\frac{E_p - U_{IB}}{U} = L^{-b} f\left( L^a \frac{t - t_c}{U} \right)$$
    where $a$ is the crossover exponent, $b$ is the leading exponent, and $t_c$ is the critical hopping strength.
  • A cost function $C$ (similar to reduced $\chi^2$) was optimized to find parameters that cause data from different $L$ to collapse onto a single universal curve.

3. Key Contributions

1. Scale Invariance of Polaron Energy: The authors demonstrate that the energy of a weakly coupled impurity is itself scale-invariant at the quantum critical point, making it a viable predictor for locating the transition.
2. Novel Scaling Exponent: They extract a crossover exponent for the polaron energy ($a \approx 0.74$) that is distinct from the known critical exponent of the host system ($1/\nu \approx 1.49$). This discrepancy presents a new theoretical challenge.
3. Connection to Correlations: The paper establishes a direct link between the polaron energy and impurity-boson density-density correlations. It shows that the flattening of these correlations at the critical point serves as an experimental signature accessible via quantum gas microscopy.
4. Experimental Feasibility: The work suggests that impurity spectroscopy (measuring $E_p$) is a powerful alternative to bulk correlation measurements for probing QPTs in cold atom experiments.

4. Key Results

A. Finite-Size Scaling of Polaron Energy

• Data Collapse: For weak impurity-boson interactions ($|U_{IB}/U| \lesssim 0.5$), the polaron energy curves for different system sizes $L$ cross near the critical point $t_c/U \approx 0.061$.
• Scaling Parameters:
  • Critical Point: $t_c/U = 0.061 \pm 0.0013$, consistent with literature values.
  • Leading Exponent ($b$): $b = 0.05 \pm 0.078$. Since this is consistent with zero, the polaron energy $E_p$ is scale invariant at the critical point (it does not vanish or diverge with system size).
  • Crossover Exponent ($a$): $a = 0.74 \pm 0.26$.
• Discrepancy: The host system's charge gap scales with $1/\nu \approx 1.49$. The polaron energy's exponent ($a \approx 0.74$) is significantly different, suggesting the impurity probes a different aspect of the critical fluctuations or that the single-impurity limit introduces unique scaling behavior not captured by standard bulk theories.

B. Impurity-Boson Correlations

• Flattening Effect: In small systems ($L=4$), the impurity-boson density-density correlation function $G^{(2)}_{IB}(d)$ exhibits a distinct "flattening" (minimal standard deviation $\sigma$) near the critical point ($t/U \approx 0.065$).
• Physical Interpretation:
  • In the Mott insulator phase (low $t/U$), the host boson density is enhanced near the impurity.
  • In the superfluid phase (high $t/U$), the density is depressed near the repulsive impurity.
  • At the critical point, the correlation profile becomes flat, hinting at a divergent correlation length in the thermodynamic limit.
• Theoretical Link: Using the Hellmann-Feynman theorem, the authors show $\frac{dE_p}{dU_{IB}} = G^{(2)}_{IB}(0)$. Furthermore, a perturbative expansion relates $E_p$ to the dynamic structure factor $S(k, \omega)$ of the host, confirming that $E_p$ probes non-local correlations.

C. Validity Limits

• The scaling hypothesis holds well for weak interactions ($|U_{IB}/U| \lesssim 0.5$).
• As interactions strengthen, the impurity significantly disturbs the host system, breaking the "passive probe" assumption and causing the data collapse to fail (increased cost function $C$).

5. Significance and Outlook

• New Probe for QPTs: The results validate impurity spectroscopy as a robust method to detect quantum critical points, particularly in regimes where traditional order parameters are hard to measure.
• Experimental Relevance: The predicted "flattening" of non-local correlations is directly observable in quantum gas microscope experiments, offering a concrete experimental signature for the Mott-superfluid transition.
• Theoretical Challenge: The discovery of a distinct crossover exponent ($a \approx 0.74$) for the polaron energy compared to the host system ($1/\nu \approx 1.49$) is an open theoretical problem. It suggests that the critical behavior of a single impurity in a quantum critical bath may belong to a different universality class or involve complex many-body effects not yet fully understood.
• Broader Impact: This work opens new research directions for using impurities to map critical properties in high-precision quantum simulations with ultracold atoms.