Towards the complete description of stationary states of a Bose-Einstein condensate in a one-dimensional quasiperiodic lattice: A coding approach

This paper establishes and numerically verifies sufficient conditions for a one-to-one correspondence between stationary states of a Bose-Einstein condensate in a one-dimensional quasiperiodic lattice and bi-infinite sequences over a finite alphabet, thereby providing a coding framework to describe these states.

The Gist

Imagine you are trying to describe every possible shape a wave of frozen atoms (a Bose-Einstein Condensate) can take when it's trapped in a special kind of cage.

Usually, if the cage is a perfect, repeating pattern (like a standard brick wall), the waves can slide around easily. If you find one wave shape, you can just slide it over to the left or right, and it's essentially the same shape. It's like having a wallpaper pattern; once you know the pattern, you know everything.

But what if the cage is quasiperiodic? Imagine a wall made of bricks, but the pattern never quite repeats. It's ordered, but it never becomes the same twice. It's like a musical rhythm that follows a complex rule (like the Fibonacci sequence) but never loops back to the start. In this chaotic-looking cage, every single wave shape is unique. There is no "sliding" allowed. If you have a million different wave shapes, they are all distinct.

The Problem:
Because every shape is unique and the pattern never repeats, it seems impossible to write a complete list or a "rulebook" for all the possible wave shapes. It feels like trying to catalog every single grain of sand on a beach that keeps changing its shape.

The Solution (The "Coding" Approach):
The authors of this paper found a clever way to organize this chaos. They realized that even though the cage is complex, the waves inside it behave like a code.

Here is the analogy:
Imagine the wave is a traveler walking through a landscape of hills and valleys (the potential energy).

1. The "Singular" Travelers: Most travelers get stuck. They try to climb a hill that is too steep and fall off the edge into infinity. In physics, these are "singular" solutions—they don't exist in the real world because the atoms would fly apart. We ignore them.
2. The "Regular" Travelers: Only a few travelers manage to stay on the path forever, never falling off. These are the "regular" stationary states we care about.

The authors discovered that these "lucky" travelers can be described by a string of symbols, like a secret code.

How the Code Works

Think of the landscape as having a few specific "safe zones" or islands where the traveler can stand without falling.

• Let's say there are three safe islands. We can label them A, B, and C.
• As the traveler moves forward in time (or space), they hop from one island to another.
• The entire history of the wave is just a sequence of these hops: ... A, C, B, A, B, B, C ...

The paper proves that:

1. Every valid, stable wave shape corresponds to exactly one of these infinite sequences of letters.
2. Every possible sequence of these letters corresponds to exactly one real wave shape.

It's like a Rosetta Stone for quantum waves. Instead of trying to draw the complex, jagged shape of the wave, you just write down the code: "A-B-C-A..." and a computer can instantly draw the exact wave for you.

The "Donut" and the "Map"

To prove this works, the authors used a mathematical tool called a Poincaré map.

• Imagine taking a snapshot of the traveler every time they pass a specific checkpoint.
• If you plot these snapshots on a map, the "safe" travelers form a shape that looks like a donut (or a stack of donuts).
• The authors showed that if you stretch and fold this donut in a specific way (like kneading dough), it creates a pattern that guarantees the code works. They called this a "Smale horseshoe" (a classic shape in chaos theory that looks like a U-shape being stretched and folded).

They developed a numerical "checklist" (an algorithm) to verify if the landscape is "kneaded" correctly. If the checklist passes, you know you can use the code.

Why This Matters

1. Complete Description: Before this, we could only find a few specific wave shapes by guessing. Now, we have a method to describe every single possible stable shape in this complex environment.
2. Predictability: If you want a wave that looks a certain way, you just write the code and generate it.
3. Real World: This applies to Bose-Einstein Condensates (super-cold atoms) trapped in laser grids. It helps physicists design experiments to create specific quantum states for computing or sensing.

The "What If" Scenarios

The paper also shows when this code fails:

• Example 2: If the landscape is too shallow, the "islands" merge together. The traveler can't tell which island they are on, and the code breaks down.
• Example 3: If the hills are too steep or twisted, the traveler might jump in a way that doesn't follow the simple A-B-C rules. The code no longer matches the reality.

Summary

In simple terms, the authors took a messy, non-repeating quantum problem and realized it could be solved by treating the solutions like a language. They proved that for certain conditions, the complex behavior of atoms in a quasiperiodic lattice can be reduced to a simple, infinite string of letters (like ...010110...).

It turns a chaotic, infinite puzzle into a structured, codable system, allowing scientists to "read" and "write" the states of matter in ways that were previously thought impossible.

Technical Summary

1. Problem Statement

The paper addresses the challenge of classifying and describing the stationary states of a one-dimensional Bose-Einstein Condensate (BEC) trapped in a quasiperiodic optical lattice.

• Context: In periodic lattices, stationary states (gap solitons) possess translational symmetry, meaning they have countable exact copies. In quasiperiodic lattices (characterized by incommensurate frequencies), translational symmetry is broken, making every stationary state unique.
• Difficulty: The interplay between quasiperiodicity and nonlinearity (described by the Gross-Pitaevskii equation, GPE) creates a complex landscape. Unlike periodic systems where spectral gaps are well-defined, quasiperiodic systems exhibit fractal spectra, Anderson localization, and mobility edges.
• Goal: The authors aim to determine if the set of all regular (non-singular) stationary states in a quasiperiodic potential can be completely classified via a coding approach, establishing a one-to-one correspondence between physical solutions and bi-infinite sequences of symbols from a finite alphabet.

2. Methodology

The authors extend a symbolic dynamics method originally developed for periodic lattices to the quasiperiodic regime. The core methodology involves:

A. Mathematical Formulation

The system is modeled by the dimensionless Gross-Pitaevskii equation:
$$\ddot{u} + (\mu - V(x))u + P(x)u^3 = 0$$
where $u(x)$ is the real-valued stationary wavefunction, $\mu$ is the chemical potential, $V(x)$ is the trap potential, and $P(x)$ is the pseudopotential (nonlinearity coefficient). The study focuses on cases where $V(x)$ and/or $P(x)$ are quasiperiodic (e.g., $V(x) = A_1 \cos x + A_2 \cos(\theta x + \delta)$ with irrational $\theta$).

B. Singular vs. Regular Solutions

• Singular Solutions: Most solutions to the Cauchy problem for this equation tend to infinity (collapse) at a finite point on the real axis. These are physically irrelevant.
• Regular Solutions: Solutions that exist for all $x \in (-\infty, \infty)$. These form a "scarce" set that the authors aim to characterize.

C. The Poincaré Map and Symbolic Dynamics

The authors define a map $P$ acting on the space of initial conditions $(\Delta, u, u')$, where $\Delta$ represents the phase shift of the quasiperiodic potential.

• $\gamma$-Islands and Donuts: They identify regions in the phase space (called $\gamma$-islands) where solutions remain regular over a specific interval. By varying the phase parameter $\Delta$, these islands deform continuously to form $\gamma$-donuts (sets $T$).
• Hypotheses for Coding:
  1. Hypothesis 1: The domain of regular solutions contains a $\gamma$-donut set $T$ composed of a finite number of disjoint donuts.
  2. Hypothesis 2 (Hyperbolicity): The map $P$ acts on these donuts in a hyperbolic manner. Specifically, it contracts "h-strips" (horizontal strips) and "v-strips" (vertical strips) by a factor $\rho > 1$, ensuring that the intersection of forward and backward iterates converges to a unique point.

D. Numerical Verification (Strip Mapping Theorems)

Directly verifying the hyperbolicity of the map is computationally difficult. The authors introduce Strip Mapping Theorems (Theorems 2 and 3) which reduce the verification to checking:

1. Monotonicity: The boundaries of the islands maintain specific monotonicity properties.
2. Linearization Signs: The signs of the entries in the Jacobian matrix ($DP$) of the map follow specific patterns (e.g., all positive or specific sign configurations).
3. Contraction: The magnitude of specific matrix entries ($|a_{11}|$) is strictly greater than 1.

3. Key Contributions

1. Extension to Quasiperiodicity: The paper successfully generalizes the symbolic dynamics coding method from periodic to quasiperiodic potentials, a significant theoretical leap given the lack of translational symmetry.
2. Rigorous Conditions for Coding: It formulates sufficient conditions (Hypotheses 1 and 2) under which a one-to-one correspondence exists between regular solutions and bi-infinite sequences.
3. Practical Numerical Algorithm: It provides a concrete algorithm (Algorithm 1) to numerically verify these conditions using linearization matrices, making the theory testable for specific physical parameters.
4. Alternative Interpretation: It offers a heuristic physical interpretation where the quasiperiodic potential is viewed as an infinite chain of non-identical wells, and the coding represents the "composition" of bound states from individual wells.

4. Results

The authors applied their method to a specific case: a BEC in a bichromatic potential $V(x) = A_1 \cos x + A_2 \cos(\theta x + \delta)$ with repulsive nonlinearity ($P(x) = -1$).

• Successful Coding (Example 1): For parameters $\mu=1, A_1=3, A_2=2$, and $\theta$ (golden ratio), the numerical verification confirmed both hypotheses.

  • The set of regular solutions corresponds exactly to the set of bi-infinite sequences over a 3-symbol alphabet $\mathcal{A}_3 = \{-1, 0, 1\}$.
  • Each symbol corresponds to a specific "island" in the phase space.
  • The authors demonstrated that distinct codes produce distinct, localized soliton solutions, even though the potential lacks translational symmetry.
  • The method was robust to rational approximations of the golden ratio (Fibonacci numbers), showing consistency with periodic limits.

• Failure Cases (Examples 2 & 3):

  • Example 2: At $\mu=0.4$, the islands in the phase space merged, violating Hypothesis 1 (no disjoint donut set). Coding failed.
  • Example 3: At $\mu=1.8$, the islands remained disjoint, but the linearization matrix signs violated the hyperbolicity conditions (Theorem 2). Coding failed.
  • Significance: This proves the coding is not automatic; it depends critically on the chemical potential and lattice amplitudes.

5. Significance

• Complete Classification: The work provides a "complete description" of the stationary states for a specific class of quasiperiodic systems, transforming an infinite-dimensional problem into a discrete symbolic one.
• Predictive Power: The coding approach allows researchers to predict the existence and structure of complex nonlinear modes (e.g., multi-soliton states) simply by constructing the corresponding symbol sequence.
• Stability and Topology: The authors suggest that the coding sequence may correlate with the stability of the solitons. Furthermore, they propose a link to Thouless pumping, where adiabatically varying the phase shift $\Delta$ corresponds to traversing the "donut" in phase space, potentially enabling topological transport of solitons.
• Experimental Relevance: Since quasiperiodic lattices are experimentally realizable with cold atoms (using superimposed laser beams with incommensurate wavelengths), this theory offers a roadmap for controlling and classifying BEC states in these complex environments.

In summary, the paper demonstrates that despite the complexity introduced by quasiperiodicity, the nonlinear stationary states of a BEC can be rigorously classified using symbolic dynamics, provided specific hyperbolicity conditions are met. This bridges the gap between the theory of quasicrystals and nonlinear wave dynamics.