Mapping Between Nonlinear Schödinger Equations with Real and Complex Potentials

This paper constructs a mapping between stationary solutions of nonlinear Schrödinger equations with real and complex potentials to derive exact solutions with real energies, specifically applying this framework to damped quantum harmonic oscillators and dissipative periodic solitons.

The Gist

Imagine you are trying to navigate a boat through a stormy sea. In the world of physics, this "sea" is described by an equation called the Nonlinear Schrödinger Equation (NLS). Usually, this equation helps us predict how waves (like light in a fiber optic cable or atoms in a Bose-Einstein condensate) move.

However, in many real-world situations, the sea isn't just stormy; it's also leaky or magical. Some parts of the sea absorb your boat's energy (dissipation), while others pump extra energy into it (amplification). In physics, we call this a complex potential.

The problem is that equations with these "leaky/magical" features are incredibly hard to solve. They are like trying to solve a maze while the walls are constantly moving. Usually, when you add these complex, energy-changing elements, the math gets so messy that the solutions become unstable or impossible to predict.

The Paper's Big Idea: The "Magic Translator"

Mario Salerno, the author of this paper, has built a translator.

He discovered a way to take a problem that looks impossible (a wave in a leaky, magical sea) and translate it into a problem that is already solved (a wave in a calm, normal sea).

Here is the analogy:

• The Hard Problem: You want to know how a boat behaves in a storm where the water is also trying to suck the boat down or push it up (Complex Potential).
• The Easy Problem: You already know exactly how a boat behaves in calm, flat water (Real Potential).
• The Translator: Salerno's method is a set of rules that says, "If you take the solution for the calm water and twist it just right (adding a specific phase and adjusting the shape), it will perfectly match the behavior in the stormy, leaky water."

How the Translator Works

The paper describes a two-step process:

1. Solve the Simple Version First:
   First, you ignore the "leaks" and "magic." You solve the equation for a standard, real-world scenario. You find the shape of the wave (the amplitude) and its energy. Let's call this the Master Wave.

2. Apply the "Twist":
   Once you have the Master Wave, Salerno's method tells you exactly how to "dress it up" to fit the complex, leaky environment.

   • It calculates a phase shift (like rotating the wave slightly).
   • It calculates the complex potential (the specific pattern of leaks and pumps) required to keep that wave stable.

The magic trick is that this translation guarantees the resulting wave has a real energy. In physics, "real energy" is crucial because it means the system is stable and physically possible. Without this translation, the energy might become imaginary (a mathematical ghost), meaning the wave would explode or vanish instantly.

What They Found (The Examples)

Salerno didn't just build the translator; he used it to find new, specific shapes of waves that were previously unknown.

• The "Elliptic" Solitons: He found that in these complex, leaky environments, waves can form stable, repeating patterns (like a train of waves) that look like elliptic functions (mathematical shapes that are more complex than simple sine waves).
• Dissipative Solitons: These are special "packets" of energy that stay together even though the environment is trying to tear them apart or drain them. Think of a whirlpool that stays in one place even though the river is flowing and losing water. Salerno showed exactly how to build these whirlpools.

Why Does This Matter?

This isn't just abstract math. It has real-world applications:

• Fiber Optics: Imagine sending internet data through glass fibers. If the fiber has imperfections that absorb or amplify light, the signal gets messy. This paper helps engineers design fibers where the signal stays strong and stable despite those imperfections.
• Quantum Computing: In the quantum world, atoms can be trapped in "optical lattices" (grids of light). If these grids absorb energy, the atoms might fall out of place. This method helps scientists figure out how to keep those atoms stable.
• PT-Symmetry: Recently, physicists discovered a special kind of symmetry (called PT-symmetry) where a system can have real energy even with complex potentials. This paper shows that you don't need that special symmetry to find stable waves; you can find them in many other complex situations too.

The Bottom Line

Mario Salerno's paper is like giving a sailor a map. Before, sailors (physicists) were afraid to sail into the "complex potential" waters because the maps were blank or the storms looked too dangerous.

Salerno said, "Don't worry. If you know how to sail in calm water, I can give you a set of instructions (a mapping) that tells you exactly how to sail in the stormy, leaky waters, and I can promise you that your boat will stay afloat with a real, stable energy."

This allows scientists to design better lasers, faster internet, and more stable quantum computers by understanding how to control waves in messy, real-world environments.

Technical Summary

1. Problem Statement

The paper addresses the challenge of finding stationary solutions for the Nonlinear Schrödinger Equation (NLS) with complex potentials. Complex potentials ($V(x) + iW(x)$) are essential for modeling physical systems involving dissipation and amplification, such as:

• Light propagation in nonlinear optical fibers with periodic modulations of the complex refractive index.
• Matter-wave solitons in Bose-Einstein condensates (BEC) trapped in absorbing optical lattices.
• Quantum dissipative systems.

While recent interest has focused on $PT$-symmetric potentials (which guarantee real spectra), the author notes that real-energy solutions can exist for generic complex potentials independent of $PT$-symmetry. The core problem is to systematically construct these exact stationary solutions with real energies (or real chemical potentials) for a broad class of complex potentials, rather than relying on specific symmetry constraints.

2. Methodology: The Mapping Approach

The author proposes a systematic mapping technique that transforms the problem of solving a complex NLS equation into a nonlocal eigenvalue problem involving only real potentials.

The Model:
The starting point is the generalized NLS equation with linear and nonlinear complex potentials:
$$i\psi_t = -\frac{1}{2}\psi_{xx} + (V_l + iW_l)\psi + (\sigma + V_{nl} + iW_{nl})|\psi|^2\psi$$
The author seeks stationary solutions of the form $\psi(x, t) = A(x)e^{i\theta(x)}e^{-i\omega t}$, where $A(x)$ is the real amplitude and $\theta(x)$ is the real phase.

The Derivation:

1. Separation of Variables: Substituting the ansatz into the NLS equation yields a coupled system for $A(x)$ and $\theta(x)$.
2. Phase Integration: By integrating the phase equation, the phase gradient $\theta_x$ is expressed in terms of an integral function $F(x)$ dependent on the complex gain/loss terms ($W_l, W_{nl}$) and the amplitude $A$.
3. The Mapping: The author defines a specific functional relationship where the complex potential terms are determined by the amplitude profile $A(x)$. Specifically, they propose an ansatz where the integral function $F(x)$ is an analytical function of $A^2$ (e.g., $F(x) = \frac{1}{2}C_n A^{n+2}$).
4. Reduction to Real Eigenvalue Problem: By substituting this ansatz back into the amplitude equation, the complex NLS is reduced to a real nonlinear eigenvalue problem:
   $$\left[ -\frac{1}{2}\frac{\partial^2}{\partial x^2} + V_l + (\sigma + V_{nl})A^2 + \frac{1}{2}\left(\frac{F(x)}{A^2}\right)^2 \right] A = \omega A$$
   This equation contains only real potentials and the real amplitude $A$.

The Inverse Process:
Once a solution $A(x)$ and eigenvalue $\omega$ are found for the real equation (analytically or numerically), the corresponding complex potentials ($W_l, W_{nl}$) and the phase $\theta(x)$ are self-consistently determined by the mapping relations. This guarantees that the resulting complex solution possesses a real energy $\omega$.

3. Key Contributions

• General Mapping Framework: The paper establishes a general method to map stationary solutions of NLS equations with real potentials to those with complex potentials. This bypasses the need for $PT$-symmetry to ensure real spectra.
• Systematic Construction of Exact Solutions: The method allows for the construction of exact analytical solutions for complex potentials using known solutions (elliptic functions) of real NLS equations.
• Handling Higher-Order Nonlinearities: The framework is flexible enough to handle cases where the mapping introduces higher-order nonlinearities (e.g., quintic terms), which can be absorbed by redefining the real linear or nonlinear potentials.
• Extension to General Functions: The mapping is generalized to include infinite series of terms ($F(x) = \sum C_n A^{n+2}$), allowing for a wide variety of complex potential profiles.

4. Results and Specific Examples

The author illustrates the method using elliptic functions (Jacobi elliptic functions: $sn, cn, dn$) as the real amplitude solutions.

• Case $n=1$ (Linear Complex Potentials):

  • Assuming $W_{nl}=0$ and $V_{nl}=0$, the mapping relates linear complex potentials to the derivative of the amplitude.
  • Result: Exact periodic dissipative solitons are constructed where the complex potential $W_l(x)$ is proportional to $A_x$ (the derivative of the amplitude). For example, if $A(x) = A_0 cn(x, k)$, the complex potential becomes $W_l(x) \propto -sn(x)dn(x)$.
  • The phase $\theta(x)$ is explicitly calculated as an integral of the amplitude.

• Case $n=2$ (Higher-Order Nonlinearities):

  • The mapping introduces quintic nonlinearities. The author shows how to balance these by adjusting the real potential $V_l$ or $V_{nl}$.
  • Result: Exact solutions are derived for pure nonlinear optical lattices ($V_l=W_l=0$) where the complex part of the nonlinear potential $W_{nl}$ is determined by the amplitude profile.

• General Case ($n>2$ and Series Expansion):

  • A more complex ansatz involving a sum of terms ($C_0 + C_2 A^2$) is used to generate solutions for mixed linear and nonlinear complex potentials.
  • Stability: The paper notes (without detailed proof) that these constructed solutions are stable under time evolution.

5. Significance

• Physical Relevance: The method provides a rigorous way to design optical or matter-wave systems with specific gain/loss profiles (complex potentials) that support stable, localized, periodic states with real energies. This is crucial for engineering photonic devices and BEC traps.
• Beyond $PT$-Symmetry: By demonstrating that real spectra exist for generic complex potentials, the work broadens the scope of non-Hermitian physics beyond the strict constraints of $PT$-symmetry.
• Analytical Utility: It transforms a difficult problem (solving NLS with complex potentials) into a well-understood problem (solving NLS with real potentials), leveraging existing analytical solutions (like elliptic functions) to generate new classes of exact solutions.
• Versatility: The approach is applicable not only to the standard NLS but also to linear Schrödinger equations (quantum dissipative oscillators) and equations with arbitrary higher-order nonlinearities.

In summary, Salerno presents a powerful mathematical tool that decouples the complexity of the potential from the solvability of the equation, allowing researchers to "reverse engineer" complex potentials that support specific, physically realizable stationary states.