Stationary Solitons in discrete NLS with non-nearest neighbour interactions

This paper employs dynamical systems methods to accurately construct stationary discrete solitons in an extended one-dimensional Discrete NLS model featuring non-nearest neighbour interactions, demonstrating how long-range coupling can induce soliton bistability with potential applications in controllable switching and biological transport.

The Gist

The Big Picture: The "Infinite Chain" of Beads

Imagine a long, infinite chain of beads connected by springs. This is a model for many things in physics, from light traveling through fiber optic cables to energy moving through biological molecules.

Usually, scientists assume that a bead only talks to its immediate neighbors (the one right next to it). If you push one bead, it pulls the one next to it, which pulls the next, and so on. This is like a game of "telephone" where you only whisper to the person standing right beside you.

This paper asks a different question: What happens if a bead can also "whisper" to the bead two spots away? What if the interaction isn't just "nearest neighbor," but extends a bit further? The authors call this "non-nearest neighbour interactions."

The Goal: Finding the "Frozen Wave"

In this chain of beads, you can send a wave down the line. Usually, waves spread out and fade away, like a ripple in a pond. But sometimes, under the right conditions, a wave can stay perfectly shaped and move (or stay still) without changing. This is called a Soliton.

Think of a soliton like a perfectly formed snowball rolling down a hill. Most snowballs melt or break apart as they roll, but this one keeps its shape perfectly.

The authors wanted to find these "frozen" or "stationary" snowballs (solitons) in their specific model where beads talk to neighbors and second-nearest neighbors.

The Challenge: The "Four-Dimensional Maze"

To find these special waves, the authors had to turn the physics problem into a math problem.

1. The Simple Case (A = 0): First, they looked at the standard model where beads only talk to immediate neighbors. They turned this into a simple 2D map (like a flat piece of paper). They found that if the "spring strength" was just right, the system was stable. It was like walking on a flat, calm lake.
2. The Complex Case (A ≠ 0): Then, they turned on the "second neighbor" interaction. Suddenly, the math got much harder. The problem jumped from a flat 2D map to a 4D maze.

The Analogy: Imagine trying to navigate a maze.

• In the simple case, the maze is on a flat floor. You can see the whole path.
• In the complex case, the maze has four dimensions. You can move forward/back, left/right, up/down, and... through time (or some other invisible dimension). It's incredibly hard to visualize or solve.

The Method: The "Parametrization Method" (The Magic Blueprint)

How do you find a path through a 4D maze? You can't just guess. The authors used a powerful mathematical tool called the Parametrization Method.

Think of the "Stationary Soliton" as a secret tunnel that starts at a specific point (the origin), wanders through the 4D maze, and eventually loops back to the exact same starting point.

• The Stable Manifold: Imagine a funnel leading into the starting point. If you drop a ball anywhere in this funnel, it will eventually roll into the center.
• The Unstable Manifold: Imagine a fountain shooting out from the starting point. If you are in this stream, you are pushed away from the center.

To find the soliton, the authors needed to find a spot where the Funnel (Stable) and the Fountain (Unstable) touch each other perfectly. If they touch, you can go out on the fountain and come back in on the funnel, creating a perfect loop.

They used a computer to build a "blueprint" (a polynomial approximation) of these funnels and fountains. They calculated the shape of these invisible surfaces with extreme precision (up to 80 levels of detail!).

The Discovery: The "Sweet Spot"

By running their calculations, they found a very specific "Sweet Spot" in the parameters (the settings of their model):

• Parameter ǫ (Epsilon): Must be positive. (Think of this as the "tightness" of the springs).
• Parameter A: Must be a specific negative number (between -0.145 and -0.115). (Think of this as the "strength of the long-distance whisper").

The Result: When they set the knobs to these specific values, the "Funnel" and the "Fountain" touched. This proved that a Stationary Soliton exists!

They didn't just prove it exists; they built a blueprint of exactly what the wave looks like.

Why Does This Matter?

1. Controllable Switching: The paper mentions that these long-range interactions can cause "bistability." Imagine a light switch that can be stuck in the "ON" position or the "OFF" position, but you can flip it back and forth with a precise nudge. These solitons could be used as tiny, controllable switches in future computers or communication devices.
2. Biological Transport: The authors mention that these models help explain how energy moves through biological molecules (like DNA). If energy can travel as a stable "snowball" (soliton) rather than spreading out and wasting away, it explains how life processes work efficiently.
3. Mathematical Triumph: They showed that even in a complex, 4D world with long-range interactions, you can find perfect, stable structures if you know exactly where to look.

Summary in One Sentence

The authors used advanced math to prove that in a chain of interacting particles where neighbors talk to each other and their neighbors' neighbors, you can create a perfect, unchanging energy wave (a soliton) if you tune the interaction strength to a very specific, narrow range.

Technical Summary

1. Problem Statement

The paper investigates the existence and construction of stationary discrete solitons in an extended one-dimensional Discrete Nonlinear Schrödinger (DNLS) model.

• Context: While standard DNLS models typically assume nearest-neighbor interactions, physical systems (such as energy transport in biological molecules or optical lattices) often exhibit long-range (non-nearest neighbor) interactions.
• The Model: The authors analyze the following equation, which includes a next-nearest-neighbor interaction term:
  $$i\dot{u}_n + |u_n|^2 u_n + \epsilon \left( u_{n+1} - 2u_n + u_{n-1} + A(u_{n+2} + u_{n-2}) \right) = 0$$
  where:
  • $u_n: \mathbb{R}^+ \to \mathbb{C}$ represents the complex amplitude at lattice site $n$.
  • $\epsilon > 0$ is a coupling parameter.
  • $A$ controls the strength of the next-nearest-neighbor interaction.
• Objective: The goal is to find steady-state solutions (where $\dot{u}_n = 0$) that correspond to stationary solitons. This requires finding orbits that are homoclinic to the origin (starting and ending at the zero state) in the associated discrete dynamical system.

2. Methodology

The authors employ a dynamical systems approach, converting the differential-difference equation into a discrete mapping. The methodology is divided into two main cases based on the parameter $A$.

Case 1: Nearest-Neighbor Only ($A = 0$)

• Reduction: Setting $A=0$ reduces the system to a 2-dimensional mapping (a generalized Hénon map):
  $$f_0(x, y) = (y, -x + 2y - \frac{1}{\epsilon}y^3)$$
• Analysis: The authors prove this map is area-preserving. They analyze the stability of the fixed point at the origin, showing it is stable if and only if $\epsilon > 0$. In this regime, no homoclinic orbits (solitons) exist because the origin is a stable sink, not a saddle.

Case 2: Long-Range Interactions ($A \neq 0$)

• Reduction: Setting $A \neq 0$ and looking for steady states ($\dot{u}_n=0$) yields a 4-dimensional mapping:
  $$f(x, y, z, w) = (y, z, w, -x - \frac{1}{A}y + \frac{2}{A}z - \frac{1}{\epsilon A}z^3 - \frac{1}{A}w)$$
  This map is volume-preserving and possesses specific symmetry properties (reversibility and involution symmetries).
• Fixed Point Analysis: The origin $(0,0,0,0)$ is a fixed point. The authors analyze its stability by examining the eigenvalues of the Jacobian matrix. They identify a parameter region where the origin becomes a hyperbolic saddle with a 2D stable manifold ($W^s$) and a 2D unstable manifold ($W^u$).
• Parametrization Method: To locate homoclinic orbits (intersections of $W^s$ and $W^u$), the authors use the Parametrization Method:
  1. Series Expansion: They represent the stable and unstable manifolds as analytic power series ($S_s$ and $S_u$) in local coordinates.
  2. Invariance Equation: They solve the functional equations $f \circ S_s = S_s \circ \Lambda_s$ and $f \circ S_u = S_u \circ \Lambda_u$, where $\Lambda$ represents the linear action of the map on the manifolds.
  3. Coefficient Calculation: This generates a linear system for the coefficients of the power series, solved recursively up to high orders (order 80).
  4. Intersection Search: They numerically solve the equation $P^u(u_1, v_1) = P^s(u_2, v_2)$, where $P$ are polynomial approximations of the manifolds. A solution indicates a homoclinic point.
  5. Transversality Check: They verify that the intersection is transverse by computing the determinant of the Jacobian matrix formed by the tangent vectors of the manifolds at the intersection point.

3. Key Contributions

• Construction of Solitons: Unlike previous studies that might only prove existence theoretically, this paper provides a constructive method to generate stationary discrete solitons with high numerical accuracy.
• Parametrization Method Application: The paper successfully applies the parametrization method to a 4-dimensional reversible map, a complex task requiring the handling of high-order polynomial coefficients and symmetry constraints.
• Parameter Space Mapping: The study rigorously defines the region in the $(\epsilon, A)$ parameter space where these solitons exist.

4. Key Results

• Existence Region: Homoclinic orbits (and thus stationary solitons) exist for:
  • $\epsilon \in (0, 1]$
  • $A \in [-0.145, -0.115]$ (specifically within the range where the origin is a hyperbolic saddle with real eigenvalues).
• Stability Conditions:
  • If $\epsilon < 0$, the origin is stable, and no solitons exist.
  • If $A=0$, the system reduces to the 2D case where the origin is stable for $\epsilon > 0$, preventing soliton formation.
  • The introduction of the non-nearest neighbor term ($A \neq 0$) is crucial for creating the saddle structure necessary for homoclinic orbits.
• Transversality: The authors verified that the intersection of the stable and unstable manifolds is transverse (the determinant of the intersection matrix is non-zero). This implies the existence of complex dynamical behavior (chaos) in the phase space and confirms the robustness of the soliton solutions.
• Numerical Accuracy: Using terms up to order 80 in the power series expansion, the authors achieved error norms of approximately $10^{-13}$ for specific test cases (e.g., $A = -0.125, \epsilon = 0.0004$).

5. Significance

• Physical Relevance: The results are significant for modeling physical systems where long-range interactions are non-negligible, such as energy and charge transport in biological molecules and optical waveguide arrays.
• Bistability and Switching: The paper notes that long-range interactions can lead to soliton bistability. The ability to construct these solutions accurately suggests potential applications in controllable switching devices in nonlinear optics and information processing.
• Methodological Advancement: The work demonstrates that dynamical systems techniques (specifically the parametrization method) are powerful tools for solving nonlinear lattice equations that are not integrable, providing a rigorous alternative to purely numerical shooting methods.

In conclusion, the paper establishes that extending the DNLS model to include non-nearest neighbor interactions fundamentally alters the phase space structure, enabling the existence of stable, stationary discrete solitons in a specific parameter regime, and provides a precise mathematical framework for their construction.