Homoclinic and heteroclinic solutions of the nonlinear Schrödinger equation with a complex Wadati potential

This paper characterizes homoclinic and heteroclinic stationary solutions of the nonlinear Schrödinger equation with a PT-symmetric Wadati potential, utilizing asymptotic analysis and numerical simulations to explore their existence, bifurcations, and role in resonant nonlinear wave generation.

The Gist

Imagine you are standing by a river, but this isn't a normal river. It's a magical, non-Hermitian river.

In a normal river, water flows, and if you throw a stone in, the ripples eventually fade away due to friction. But in this magical river, there are hidden forces at play:

1. Gain: Some parts of the river are secretly being fed by underground springs, making the water flow faster and higher (amplifying the waves).
2. Loss: Other parts are draining into a giant sinkhole, sucking the energy out of the water (damping the waves).
3. The Balance: The river is designed so that the "springs" and "sinks" are perfectly mirrored. If you look at the left side, it's the exact opposite of the right side. This is called PT-symmetry.

The scientists in this paper (Chandramouli, Sprenger, and Hoefer) are trying to figure out what happens when you send a steady stream of water (a "plane wave") through this magical river, especially when it hits a specific kind of bump or barrier in the middle.

The Main Characters: Homoclinic and Heteroclinic Waves

The paper focuses on two special types of "flow patterns" that can form in this river. Think of them as the river's way of reacting to the gain and loss.

1. Homoclinic Solutions: The "Boomerang" Wave

Imagine you throw a ball into a valley. If the valley is shaped just right, the ball might roll down, hit a bump, roll back up the other side, and then roll all the way back to the exact spot where it started, stopping there.

• In the paper: This is a wave that starts as a flat, calm stream far away, gets disturbed by the gain/loss in the middle, creates a big bump (or a dip) in the water level, and then smoothly returns to being a flat, calm stream far away on the other side.
• The Analogy: It's like a surfer who paddles out, catches a wave that goes up and down, and then paddles back to the exact same calm water they started from.
• Two Types:
  • Depression Waves: The water level dips down in the middle (like a hole).
  • Elevation Waves: The water level rises up in the middle (like a hill).

2. Heteroclinic Solutions: The "Bridge" Wave

Now imagine the river on the left side is flowing fast, and the river on the right side is flowing slow. You can't just stop; you have to transition.

• In the paper: This is a wave that starts as one type of flow on the left (e.g., fast and calm) and transitions smoothly into a completely different type of flow on the right (e.g., slow and calm). It acts like a bridge connecting two different states of the river.
• The Analogy: Think of a car driving from a highway (fast) onto a quiet country road (slow). A "heteroclinic" solution is the perfect, smooth ramp that connects the two without a crash or a sudden stop.

The "Magic" of the Wadati Potential

The scientists use a specific mathematical shape for the river's barrier called the Wadati potential. You can think of this as the "shape of the hill" the water has to go over.

• The Shape: It's a smooth, bell-shaped hill (like a seashell curve).
• The Twist: The hill isn't just a physical bump; it's also a "magic zone" where the water gains energy on one side and loses it on the other.
• The Result: Because of this magic, the river can support these special "Boomerang" and "Bridge" waves that would not exist in a normal, boring river. In a normal river, the energy would just dissipate, and the wave would die out. Here, the gain and loss balance each other out to keep the wave alive.

The "Traffic Jam" Analogy (Transcritical Flow)

The paper mentions something called "transcritical flow." Imagine a highway where cars are driving at different speeds.

• Subsonic: Cars are driving slowly (under the speed limit).
• Supersonic: Cars are speeding (over the speed limit).
• Transcritical: The cars are driving exactly at the speed limit.

In normal physics, when traffic hits a bottleneck at the speed limit, chaos ensues (shock waves, traffic jams). The paper shows that in this magical, non-Hermitian river, the water doesn't just crash; it forms these beautiful, stable "Boomerang" and "Bridge" patterns. It's as if the traffic jam organizes itself into a perfect, rhythmic dance instead of a pile-up.

Why Does This Matter?

You might ask, "Who cares about magical rivers?"

1. Lasers and Light: This math describes how light behaves in special materials (optics) that have gain (amplifiers) and loss (absorbers). This could help us build better lasers or optical computers.
2. Superfluids: It helps explain how super-cold fluids (like those in quantum experiments) behave when they flow past obstacles.
3. New Physics: It proves that when you mix "gain" and "loss" in a balanced way, nature creates entirely new types of waves that we've never seen before.

The Big Takeaway

The scientists used math and computer simulations to map out exactly where these special waves can exist. They found that:

• If the "river" is too fast or too slow, the waves disappear.
• There are specific "sweet spots" (like a Goldilocks zone) where these Boomerang and Bridge waves can exist perfectly.
• They discovered that these waves can change shape, split, or merge depending on the speed of the flow and the shape of the barrier.

In short, they found the recipe for creating stable, self-sustaining waves in a world where energy is constantly being added and removed, revealing a hidden order in what looks like chaos.

Technical Summary

1. Problem Statement

The paper investigates stationary solutions of the Nonlinear Schrödinger (NLS) equation subject to a PT-symmetric, complex linear potential known as the Wadati potential. The governing equation is:
$$i\psi_t = -\frac{1}{2}\psi_{xx} + V(x)\psi + |\psi|^2\psi$$
where the potential is $V(x) = \frac{1}{2}(-w(x)^2 + iw'(x))$. Here, $w(x)$ is a real-valued, smooth, even function (e.g., $\text{sech}^n(x)$) representing a spatially balanced gain-loss profile (gain where $w$ increases, loss where it decreases) and a repulsive real potential.

The primary objective is to characterize homoclinic (solutions asymptoting to the same nonlinear plane wave at $\pm \infty$) and heteroclinic (solutions asymptoting to different plane waves) solutions. These solutions model steady flow patterns in non-Hermitian dispersive media, analogous to transcritical flow in fluid dynamics where background flow speeds approach the long-wave speed, potentially leading to resonant wave generation (e.g., Dispersive Shock Waves).

2. Methodology

The authors employ a combination of asymptotic analysis, hydrodynamic formulation, and numerical simulations:

• Hydrodynamic Formulation: The NLS equation is transformed using the Madelung transformation ($\psi = \sqrt{\rho}e^{i\phi}$) into a system of equations for density ($\rho$), velocity ($u = \phi_x$), and retarded momentum ($m = \rho(u-w)$). This reveals the system as a dispersive regularization of compressible Euler fluid dynamics with gain/loss and a repulsive force.
• Steady State Reduction: Assuming steady solutions ($\partial_t = 0$), the system reduces to a third-order non-autonomous Ordinary Differential Equation (ODE) system. A conserved "zero-energy" integral ($E=0$) is derived, constraining the dynamics to a specific energy surface in the $(\rho, \rho', m)$ phase space.
• Asymptotic Regimes:
  • Hydraulic Limit: Analysis of slowly varying potentials ($w(\epsilon x)$) where dispersion is neglected, reducing the problem to algebraic equations.
  • Far-Field Linearization: Investigation of decay rates and resonance conditions (supersonic vs. subsonic regimes) to determine solution polarity (depression vs. elevation).
  • Large Density/Large Velocity Limits: Perturbation expansions to connect solutions to dark solitons or forced harmonic oscillators.
• Numerical Methods:
  • Spectral Methods: Fourier collocation and Levenberg-Marquardt solvers to solve the boundary value problem (BVP) for homoclinic orbits.
  • Shooting Methods: Used to compute resonant heteroclinic solutions with oscillatory tails by enforcing symmetry conditions at the origin.
  • Continuation: Parameter continuation in density ($\rho_0$) and velocity ($u_0$) to map existence regions and bifurcation points.

3. Key Contributions and Results

A. Homoclinic Solutions (Same Far-Field State)

The study identifies two distinct families of homoclinic solutions bifurcating from constant intensity (CI) waves, separated by a critical density threshold $\rho_{cr} \approx 0.25$ (for $w=\text{sech}(x)$):

1. Depression Waves ($\rho_0 > \rho_{cr}$):

   • Occur in the subsonic regime ($u_0 < \sqrt{\rho_0}$).
   • Characterized by a dip in density at the potential peak.
   • Exist in a two-parameter region in the $(\rho_0, u_0)$ plane bounded by the sonic curve and a specific existence boundary $u_0 = d(\rho_0)$.
   • Asymptotically approach the dark soliton limit in the large density regime.

2. Elevation Waves ($\rho_0 < \rho_{cr}$):

   • Also subsonic but characterized by a peak in density.
   • The shifted momentum profile is sign-indefinite, a feature not captured by standard hydraulic theory.
   • Resonant Elevation Waves: In the supersonic regime ($u_0 > \sqrt{\rho_0}$), solutions do not decay to a constant but exhibit oscillatory tails (delocalized density, localized momentum). These are "generalized elevation waves" homoclinic to a periodic orbit.
   • Analytical formulas for tail amplitudes are derived, showing exponential dependence on the inverse velocity and potential scale.

B. Heteroclinic Solutions (Different Far-Field States)

The paper classifies heteroclinic solutions into two types based on boundary conditions:

1. Type-I Heteroclinics:

   • Connect states with equal density ($\rho_- = \rho_+$) but opposite velocities ($u_- = -u_+$).
   • These solutions emerge via a symmetry-breaking bifurcation from the constant intensity limit at the critical density $\rho_{cr}$.
   • They exist in a 2D phase plane and are unique to the non-conservative (gain-loss) system; they have no analog in conservative NLS.

2. Type-II Heteroclinics (Kinks and Antikinks):

   • Connect states with different densities and velocities ($\rho_- \neq \rho_+$).
   • In the hydraulic limit, these arise from a saddle-node bifurcation where a pair of homoclinic solutions (one depression, one elevation) merge and annihilate.
   • Numerical results show these solutions bifurcate from Type-I solutions in the small-density limit.
   • They form one-parameter families where one far-field state is supersonic and the other subsonic.

C. Phase Space Geometry

The authors map the existence regions on the zero-energy surface.

• For low Mach numbers ($M_0 < 2$), the energy surface is connected.
• For high Mach numbers ($M_0 > 2$), the surface becomes disconnected.
• Homoclinic and heteroclinic orbits are constrained to lie on these surfaces, with their existence dictated by the intersection of the zero-energy manifold and the specific boundary conditions.

4. Significance

• Non-Hermitian Hydrodynamics: The work establishes a rigorous framework for "dispersive hydrodynamics" in systems with gain and loss, extending classical fluid concepts (like transcritical flow and shock waves) to non-Hermitian optics and quantum systems.
• Resonant Wave Generation: The identification of resonant elevation waves with oscillatory tails provides a mechanism for the spontaneous generation of nonlinear wave trains (DSWs) in non-conservative media, relevant to polariton condensates and optical systems.
• Bifurcation Theory: The paper details novel bifurcation scenarios (symmetry breaking, saddle-node merging) specific to PT-symmetric potentials, distinguishing them from conservative systems.
• Experimental Relevance: The findings offer theoretical predictions for experimental setups in non-Hermitian photonics (e.g., waveguide arrays with balanced gain/loss) and exciton-polariton condensates, where such stationary flow patterns and resonant transitions can be observed.

In summary, this paper provides a comprehensive characterization of stationary wave structures in PT-symmetric NLS systems, revealing a rich landscape of homoclinic and heteroclinic solutions that bridge the gap between conservative soliton dynamics and non-conservative resonant wave generation.