🔬 condensed matter

Superfluid dome in the spatially modulated two-dimensional XY model

By combining tensor network methods and Monte Carlo simulations on a spatially modulated two-dimensional XY model, this study reveals a non-monotonic "superfluid dome" in the critical temperature caused by the effective pinning of vortices in modulation valleys, offering new insights into the interplay between superconductivity and charge density waves.

Original authors: Feng-Feng Song, Aditya Chugh, Hanggai Nuomin, Naoki Kawashima, Alexander Wietek

Published 2026-08-11

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Feng-Feng Song, Aditya Chugh, Hanggai Nuomin, Naoki Kawashima, Alexander Wietek

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a world where tiny particles, like electrons, don't just bounce around randomly but dance in perfect unison. When they do this, they create a phenomenon called superconductivity, where electricity flows with zero resistance, like a car zooming down a highway with no friction or traffic jams. This is the dream of many physicists: to build machines that run on this frictionless power. But in the real world, these electrons are often messy. They sometimes get stuck in patterns called "charge density waves," which are like traffic jams where electrons line up in rows. For a long time, scientists thought these two behaviors—superconducting flow and traffic-jam patterns—were enemies that fought for dominance. However, recent discoveries suggest they might actually be partners, intertwining in complex ways that could unlock new types of materials. The big question is: how do these two dance together, and can we tune their partnership to make the superconducting flow stronger?

This paper dives into that mystery by building a digital playground, a simplified model of a two-dimensional grid where these electron dances happen. The researchers, led by Feng-Feng Song and colleagues, wanted to see what happens when you introduce a rhythmic "modulation" to the grid. Imagine the dance floor isn't flat; instead, it has gentle hills and valleys. In some spots, the electrons can hold hands tightly (strong coupling), and in others, they hold hands loosely (weak coupling). By changing the size of these hills and valleys, the team simulated how the temperature at which superconductivity occurs ( $T_c$ ) would change. They used powerful computer methods, including tensor networks (which are like advanced puzzle solvers for quantum states) and Monte Carlo simulations (which are like rolling digital dice millions of times to see what happens), to watch the dance unfold.

The results were surprising and shaped like a hill. The team found that the temperature at which superconductivity starts doesn't just go up or down as they changed the size of the hills. Instead, it formed a "dome." If the hills were too close together, the superconductivity was weak. If they were too far apart, it was also weak. But there was a "sweet spot" in the middle where the superconducting temperature peaked. Even more interesting, as they made the "valleys" deeper (increasing the modulation strength), this sweet spot moved to larger distances between the hills.

Why does this happen? The paper suggests a clever mechanism involving "vortices." Think of a vortex as a tiny whirlpool or a tornado in the electron dance. In a perfect, flat dance floor, these whirlpools can pop up anywhere and ruin the flow. But in this hilly landscape, the whirlpools get stuck in the valleys. When the valleys are just the right size, the whirlpools are effectively "pinned" or trapped there, preventing them from ruining the superconducting dance. This pinning allows the superconductivity to survive at higher temperatures. However, if the valleys get too wide, the whirlpools have too much room to roam and escape, breaking the flow again. The researchers confirmed this by comparing their results to simpler models and using renormalization group analysis (a mathematical way of zooming out to see the big picture), which showed that the interplay between the energy cost of creating a whirlpool and the entropy (disorder) of where it can hide creates this dome shape.

In short, the paper doesn't just say "superconductivity is good"; it reveals a specific, non-monotonic relationship where the structure of the material's landscape dictates the success of the superconducting state. The authors suggest that this "superfluid dome" is a real feature of intertwined orders, offering a new way to think about how to design materials. While these findings come from simulations and mathematical models rather than a physical lab experiment, they provide a strong theoretical framework. The team proposes that future experiments, perhaps using microscopes that can see individual atoms, could look for these pinned whirlpools in real materials to see if nature follows the same rules they simulated. This work doesn't solve the mystery of high-temperature superconductivity overnight, but it adds a crucial piece to the puzzle, showing that sometimes, a little bit of disorder (the hills and valleys) can actually help create order.

Technical Summary: Superfluid Dome in the Spatially Modulated Two-Dimensional XY Model

Problem Statement
In strongly correlated electron systems, the interplay between superconductivity (SC) and charge density waves (CDW) is a central theme, often manifesting as coexisting or intertwined orders. Recent theoretical work suggests that the superconducting order parameter can couple to modulations in the electronic charge density, leading to phenomena such as pair density waves and fragmented superconducting condensates. While numerical studies on Hubbard models have begun to elucidate these states, a minimal model capturing the specific mechanism of how spatial modulation of the coupling strength affects the superconducting transition temperature ( $T_c$ ) remains to be fully characterized. Specifically, the dependence of $T_c$ on the wavelength of such modulations, and the underlying topological mechanisms driving this dependence, require investigation.

Methodology
The authors investigate a two-dimensional XY model on a square lattice with a spatially modulated coupling strength $J_{ij}$ along the $x$ -direction, defined as $J_{ij} = J_0 [1 + \alpha \sin(k x_i)]$ , where $\alpha$ is the modulation amplitude and $\lambda = 2\pi/k$ is the modulation wavelength. The coupling along the $y$ -direction remains uniform ( $J_0$ ). This model is derived from an intertwined Ginzburg-Landau theory where the superconducting amplitude is modulated by a periodic charge density.

To analyze the phase diagram and critical properties, the study employs a combination of:

Tensor Network Methods: Specifically, the Variational Uniform Matrix Product State (VUMPS) algorithm is used to solve the row-to-row transfer operator for the infinite system. The partition function is represented as a contraction of local tensors derived via a dual transformation. Critical temperatures are identified using the entanglement entropy ( $S_E$ ) of the fixed-point matrix product state (MPS), and correlation functions are calculated using impurity tensors.
Monte Carlo (MC) Simulations: Metropolis and Wolff algorithms are utilized for finite-size systems ( $L \times L$ up to $L=128$ ) to study vortex configurations, internal energy, and helicity moduli.
Modulated $q$ -state Clock Models: These are studied to distinguish between domain-wall-driven and vortex-driven transitions, particularly for $q \ge 5$ .
Renormalization Group (RG) Analysis: A hybrid approach combining a single Migdal-Kadanoﬀ (MK) RG step with Kramers-Wannier duality equations is used to estimate critical temperatures and understand the scaling of effective couplings.

Key Results

Non-Monotonic "Superfluid Dome": The critical temperature $T_c$ exhibits a non-monotonic, dome-like dependence on the modulation wavelength $\lambda$ . For a fixed modulation amplitude $\alpha$ , $T_c$ initially increases with $\lambda$ , reaches a peak, and then decreases as $\lambda$ becomes very large.
Shift of Peak Position: As the modulation amplitude $\alpha$ increases, the peak of the $T_c$ dome shifts to larger wavelengths.
Critical Behavior:
- Below $T_c$ , the system exhibits quasi-long-range order (QLRO) with algebraic decay of correlation functions ( $G(r) \sim r^{-\eta}$ ) and a central charge $c=1$ , consistent with the Berezinskii-Kosterlitz-Thouless (BKT) universality class.
- Above $T_c$ , correlations decay exponentially. The decay is anisotropic, with correlations along the modulation direction ( $x$ ) decaying faster than those perpendicular to it ( $y$ ).
- The specific heat peaks broaden with increasing $\lambda$ , reflecting an extended vortex unbinding process.
Vortex Pinning Mechanism: MC simulations reveal that vortices preferentially localize in the "valleys" of the modulation (regions of weak coupling).
- At small $\lambda$ , the valleys are narrow, restricting vortex formation and leading to a lower $T_c$ due to strong anisotropy.
- As $\lambda$ increases, the valleys widen, allowing for more vortex entropy. However, the effective interaction strength also renormalizes. The dome shape arises from a competition: initially, increasing $\lambda$ strengthens the effective interaction (raising $T_c$ ), but further increases facilitate vortex proliferation and entropy (lowering $T_c$ ).
Clock Model Comparison: In modulated $q$ -state clock models, the lower transition temperature $T_{c1}$ (driven by domain walls) increases monotonically with $\lambda$ . In contrast, the higher transition $T_{c2}$ (driven by vortices) for $q \ge 5$ exhibits the same dome-shaped dependence as the XY model, confirming that vortex physics is the origin of the dome.
RG Analysis Limitations: The hybrid MK-RG approach successfully captures the monotonic increase of the effective coupling with $\lambda$ and the suppression of order with increasing $\alpha$ . However, it overestimates $T_c$ for the XY limit and fails to capture the dome structure without explicitly accounting for vortex excitations, highlighting the necessity of topological defects in explaining the phenomenon.

Significance and Claims
The paper claims to provide a microscopic understanding of the "superfluid dome" observed in the context of intertwined superconducting and charge-ordered states. By mapping the problem to a modulated XY model, the authors demonstrate that the non-monotonic dependence of $T_c$ on modulation wavelength is a generic feature arising from the effective pinning of vortices in the weak-coupling regions of the modulation.

The findings suggest that the complex phase behavior in materials where CDW and SC coexist (such as cuprates and kagome metals) can be understood through the thermodynamic competition between vortex energy and entropy in a spatially inhomogeneous environment. The work offers a theoretical framework where the modulation wavelength (related to doping density) acts as a control parameter for optimizing superconductivity, with an optimal wavelength emerging from the interplay of pinning and thermal fluctuations. The authors propose that experimental techniques like scanning tunneling microscopy could verify the predicted vortex pinning to charge density stripes.

More like this