Superfluid response of bosonic fluids in composite optical potentials: angular dependence and Leggett's bounds

This paper investigates the superfluid response of dilute bosonic fluids in two-dimensional composite optical potentials by establishing conditions for isotropy, deriving analytical expressions for Leggett's bounds to identify optimal measurement directions, and confirming these findings through numerical simulations.

The Gist

Imagine a crowd of people at a party. If the room is empty and smooth, everyone can glide around effortlessly, moving in perfect unison without bumping into each other. This is like a superfluid: a state of matter where particles flow with zero friction, like a ghostly dance.

Now, imagine you put obstacles in the room—furniture, pillars, or walls. Suddenly, the crowd has to navigate around these things. Some people get stuck, some have to slow down, and the "perfect glide" is broken. In physics, this is what happens when you put a superfluid inside a special kind of light trap called an optical lattice.

This paper is about figuring out exactly how much of that perfect glide remains when the obstacles are arranged in complex, beautiful patterns.

Here is the breakdown of their discovery, using some everyday analogies:

1. The Setup: The "Light Trap"

Scientists use lasers to create invisible grids for atoms. By overlapping different laser beams, they can create grids that look like squares, triangles, honeycombs (Kagomé), or even weird, non-repeating patterns (quasicrystals).

Think of these laser grids as mazes.

• The Question: If I push this crowd of atoms through a square maze, a triangular maze, or a weird 5-pointed star maze, how much of their "super-speed" do they lose? Does the direction I push them matter? (e.g., is it easier to push them North-South than East-West?)

2. The Big Surprise: The "Magic Shield"

Usually, if you have a triangular grid, you expect the physics to look different depending on which way you look at it. It's like a snowflake: it has six points, so it looks different if you rotate it by 30 degrees.

The authors discovered something counter-intuitive: For many of these complex laser mazes, the superfluid doesn't care which way you push it.

• The Analogy: Imagine you are walking through a forest. If the trees are planted in a perfect square grid, you might find it easier to walk between the rows than diagonally through the trees. But the authors found that for certain complex patterns (like the Kagomé lattice or 5-pointed stars), the forest feels perfectly smooth in every direction, even though the trees are arranged in a rigid, angular pattern.
• Why? The "symmetry" of the light pattern in the "frequency world" (Fourier space) acts like a geometric shield. It cancels out all the directional biases, making the superfluid response look perfectly round and smooth, regardless of the jagged edges of the laser grid.

3. The "Leggett Bounds": The Estimator Tool

Measuring exactly how much superfluid you have is hard. It usually requires complex experiments where you drag the fluid and measure the drag.

Enter Anthony Leggett, a Nobel Prize-winning physicist. He came up with a clever trick: The Bounds.

• The Analogy: Imagine you want to guess the weight of a mystery box without lifting it. You can't know the exact weight, but you can say, "It's definitely heavier than 5 lbs (Lower Bound) and definitely lighter than 10 lbs (Upper Bound)."
• Leggett's bounds use the density map of the atoms (where they are crowded and where they are sparse) to give you a "bracket" for the superfluid fraction.
• The Paper's Contribution: The authors realized that these bounds aren't the same in every direction.
  • If you measure along the "grain" of the light pattern, your estimate might be loose (a wide gap between 5 and 10).
  • If you measure at a specific "sweet spot" angle, the gap shrinks, and your estimate becomes very tight (maybe 7.1 to 7.2).
  • They figured out exactly which angles give the best, tightest estimates for different laser patterns.

4. The "Perfect Match"

They found that for square patterns, the upper and lower bounds actually meet in the middle. It's like guessing the weight of the box and getting the exact number immediately. This happens because square patterns are "separable"—they act like independent X and Y walls that don't mess with each other.

For the more complex patterns (like the 5-pointed quasicrystal), the bounds don't meet perfectly, but the authors showed you exactly where to look to get the best possible answer.

Summary: Why Does This Matter?

1. Symmetry is Powerful: Even if the "cause" (the laser grid) is angular and jagged, the "effect" (the superfluid flow) can be perfectly smooth and round. It's a beautiful example of how geometry protects physical properties.
2. Better Measurements: By knowing the optimal angles to measure, experimentalists can get much more accurate data about how "super" their fluids are without needing to do the most difficult, time-consuming experiments.
3. Future Tech: Understanding how atoms flow through these complex light mazes is a stepping stone toward building quantum computers and simulating exotic materials (like neutron stars or high-temperature superconductors) right here in a lab.

In a nutshell: The team proved that certain complex laser mazes act like "magic carpets" for atoms, allowing them to flow smoothly in all directions, and they gave scientists a new map to find the exact spots where they can measure this flow most accurately.

Technical Summary

1. Problem Statement

The paper investigates the superfluid properties of dilute Bose-Einstein condensates (BECs) subjected to two-dimensional composite optical potentials. These potentials are generated by superposing multiple laser beams to create complex lattice structures (e.g., triangular, Kagomé, superlattices, and quasicrystals).

The core scientific questions addressed are:

• Isotropy of Response: Do these potentials, which possess only discrete rotational symmetry (e.g., $N$-fold symmetry), induce an anisotropic superfluid response, or can they yield a fully isotropic superfluid fraction tensor?
• Leggett's Bounds: How do the upper and lower bounds for the superfluid fraction (proposed by A.J. Leggett) depend on the measurement direction relative to the potential's geometry?
• Optimization: What are the optimal measurement axes that minimize the gap between Leggett's bounds, providing the most accurate estimation of the superfluid fraction?

2. Methodology

The authors employ a combination of analytical perturbation theory and numerical simulations:

• Theoretical Framework:

  • Linear Response Theory: The superfluid fraction tensor $f_{\alpha\beta}$ is defined via the response of the system to a weak dragging potential moving with velocity $v_0$.
  • Gross-Pitaevskii Equation (GPE): The system is modeled using the GPE in the frame co-moving with the drag. The wavefunction is expanded perturbatively in terms of the external potential strength ($V$) and the drag wavevector ($k_0$).
  • Composite Potentials: The external potential is modeled as a sum of "shells" of plane waves: $V(r) = \sum_{l=1}^M V_l \sum_{j=1}^{N_l/2} \cos(q_{lj} \cdot r)$, where each shell $l$ consists of $N_l$ waves arranged in a regular polygon in momentum space.
  • Leggett's Bounds: The authors derive analytical expressions for Leggett's upper ($f^+$) and lower ($f^-$) bounds within the perturbative regime, analyzing their angular dependence relative to the measurement axis.

• Numerical Validation:

  • The authors solve the GPE numerically for specific experimental configurations: a Kagomé lattice and a five-fold quasicrystal.
  • Simulations cover both the perturbative regime and strong potential strengths to test the robustness of the analytical predictions.

3. Key Contributions and Results

A. Isotropic Superfluid Response from Discrete Symmetry

• Finding: The authors prove that any composite potential constructed from superpositions of regular polygons in Fourier space (e.g., square, triangular, hexagonal, Kagomé, and specific quasicrystals) yields a fully isotropic superfluid response ($f_{xx} = f_{yy}$ and $f_{xy} = 0$), even though the potential itself only has discrete rotational symmetry.
• Mechanism: Using perturbation theory, they show that the anisotropic terms in the superfluid fraction tensor involve sums of the form $\sum \hat{q}_{x}^2 - \hat{q}_{y}^2$ and $\sum \hat{q}_x \hat{q}_y$ over the vertices of a regular polygon. These sums identically vanish due to the symmetry of the roots of unity.
• Robustness: This isotropy holds to all orders in the potential strength, not just the perturbative limit. The "geometrical protection" provided by the Fourier structure ensures the response remains rotationally invariant regardless of how strong the lattice potential is.
• Verification: Numerical results for Kagomé and five-fold quasicrystal lattices confirm that the superfluid fraction remains isotropic even in the non-perturbative regime.

B. Angular Dependence of Leggett's Bounds

• Finding: Unlike the superfluid fraction itself, Leggett's bounds exhibit strong angular dependence. The tightness of the bounds (the gap between $f^+$ and $f^-$) varies significantly with the angle $\phi$ between the measurement axis and the potential's lattice vectors.
• Optimal Directions:
  • Upper Bound ($f^+$): Is tightest (closest to the true value) when the measurement axis is aligned with the dominant Fourier components of the potential.
  • Lower Bound ($f^-$): Is tightest when the measurement axis is perpendicular to the dominant Fourier components.
• Tightest Bracketing: The authors derive a condition for the tightest possible window between the bounds. They find that for square lattices (and superlattices composed of them, where $N_l=4$), the bounds coincide ($f^+ = f^-$) in the perturbative limit. This is because the potential is separable ($V(x,y) = V(x) + V(y)$), eliminating cross-terms that create the gap.

C. Analytical Expressions

The paper provides explicit perturbative formulas for the bounds as a function of angle $\phi$:

• $f^+(\phi)$ and $f^-(\phi)$ are expressed in terms of the weighted potential amplitudes $\tilde{V}_l$ and the number of waves $N_l$ in each shell.
• These expressions allow experimentalists to predict exactly which direction to measure to minimize uncertainty in the superfluid fraction estimation.

4. Significance and Implications

• Resolution of Symmetry Paradox: The work demonstrates a counter-intuitive physical phenomenon where the "effect" (isotropic superfluidity) possesses higher symmetry than the "cause" (discrete symmetric potential). This validates Pierre Curie's principle that effects can be more symmetric than their causes in specific physical contexts.
• Experimental Guidance: The findings provide a practical roadmap for ultracold atom experiments. Researchers can now:
  1. Confidently assume isotropic superfluid response in complex lattices (like Kagomé or quasicrystals) without needing to measure all tensor components.
  2. Optimize the measurement of superfluid fractions by aligning the drag force along specific crystallographic directions to minimize the error bars provided by Leggett's bounds.
• Theoretical Tool: The derivation of the angular dependence of Leggett's bounds offers a rigorous method to estimate superfluidity in density-modulated systems without performing complex dynamical simulations or measurements.
• Future Directions: The authors conjecture that these geometric protections and optimal bound directions hold even for strongly interacting regimes and finite temperatures, suggesting a new avenue for research in quantum many-body physics.

In summary, this paper bridges the gap between the geometric structure of optical lattices and the macroscopic transport properties of quantum fluids, establishing that specific composite potentials guarantee isotropic superfluidity while providing a precise analytical framework for optimizing the measurement of superfluid fractions.