Stability of new exact solutions of the nonlinear Schrödinger equation in a Pöschl–Teller external potential

Journal Article · 17 November 2017 · Journal of Physics. A, Mathematical and Theoretical

DOI: https://doi.org/10.1088/1751-8121/aa9006 · OSTI ID:1514934

Dawson, John F. ^[1]; Cooper, Fred ^[2]; Khare, Avinash ^[3]; Mihaila, Bogdan ^[4]; Arévalo, Edward ^[5]; Lan, Ruomeng ^[6]; Comech, Andrew ^[7]; Saxena, Avadh  ^[8]

---

1. Univ. of New Hampshire, Durham, NH (United States)
2. The Santa Fe Institute, Santa Fe, NM (United States); Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
3. Savitribai Phule Pune Univ., Pune (India)
4. National Science Foundation, Arlington, VA (United States); Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
5. Pontifical Catholic Univ. of Chile, Santiago (Chile)
6. Texas A & M Univ., College Station, TX (United States)
7. Texas A & M Univ., College Station, TX (United States); St. Petersburg State Univ., St. Petersburg (Russia); Institute for Information Transmission Problems, Moscow (Russia)
8. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

Here, we discuss the stability properties of the solutions of the general nonlinear Schrödinger equation in 1 + 1 dimensions in an external potential derivable from a parity-time ($$ \newcommand{\PT}{\mathcal{PT}} \PT$$ ) symmetric superpotential $W(x)$ that we considered earlier in Kevrekidis et al (2015 Phys. Rev. E 92 042901). In particular we consider the nonlinear partial differential equation $$ \{i \, \partial_t + \partial_x^2 - V(x) + g \vert \psi(x, t) \vert ^{2\kappa} \} \, \psi(x, t) = 0 \>, $$ for arbitrary nonlinearity parameter κ, where $$g= \pm1$$ and V is the well known Pöschl–Teller potential which we allow to be repulsive as well as attractive. Using energy landscape methods, linear stability analysis as well as a time dependent variational approximation, we derive consistent analytic results for the domains of instability of these new exact solutions as a function of the strength of the external potential and κ. For the repulsive potential we show that there is a translational instability which can be understood in terms of the energy landscape as a function of a stretching parameter and a translation parameter being a saddle near the exact solution. In this case, numerical simulations show that if we start with the exact solution, the initial wave function breaks into two pieces traveling in opposite directions. If we explore the slightly perturbed solution situations, a 1% change in initial conditions can change significantly the details of how the wave function breaks into two separate pieces. For the attractive potential, changing the initial conditions by 1% modifies the domain of stability only slightly. For the case of the attractive potential and negative g perturbed solutions merely oscillate with the oscillation frequencies predicted by the variational approximation.

• Research Organization: Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
• Sponsoring Organization: USDOE Laboratory Directed Research and Development (LDRD) Program
• Grant/Contract Number: 89233218CNA000001
• OSTI ID: 1514934
• Report Number(s): LA-UR-17-23542
• Journal Information: Journal of Physics. A, Mathematical and Theoretical, Vol. 50, Issue 50; ISSN 1751-8113
• Publisher: IOP PublishingCopyright Statement
• Country of Publication: United States
• Language: English