Buslaev, Komech, Kopylova and Stuart (Comm. PDE 2008) proved that solutions of the one-dimensional Schrödinger equation coupled to a nonlinear oscillator — a nonlinear Schrödinger equation whose nonlinearity is concentrated at a single point via a Dirac delta — converge, for initial data close to a linearly stable solitary wave and under a spectral condition excluding internal modes, to a solitary wave plus dispersive radiation, with remainder decaying like ∣t∣−1/2|t|^{-1/2}. The authors explicitly left two problems open: (i) whether these soliton-plus-radiation asymptotics hold for every finite-energy solution, i.e. global attraction to the solitary manifold; and (ii) the extension to the parameter regime in which the linearization possesses an oscillatory internal mode. Problem (ii) was substantially resolved by Komech, Kopylova and Stuart (CPAA 2012) for the resonant sub-range of the internal-mode regime under a generic Fermi Golden Rule condition, leaving only the non-resonant sub-range untreated. Problem (i) — global attraction for a nonlinear Schrödinger equation with point-concentrated nonlinearity — remains fully open: the corresponding Klein–Gordon, wave, and Dirac results of the Komech program do not transfer to Schrödinger dispersion.

Background

The nonlinear Schrödinger equation (NLS) with a point-concentrated nonlinearity is a model Hamiltonian dispersive system in which the free one-dimensional Schrödinger equation is coupled, at a single point, to a nonlinear oscillator:

iψ˙(x,t)=−ψ′′(x,t)−δ(x)F(ψ(0,t)),x∈R,i\dot\psi(x,t)=-\psi''(x,t)-\delta(x)F(\psi(0,t)),\qquad x\in\mathbb R,

where δ(x)\delta(x) is the Dirac delta and the oscillator force FF derives from a potential, F(ψ)=−∇U(ψ)F(\psi)=-\nabla U(\psi). When UU is invariant under phase rotations (the U(1)U(1) symmetry), i.e. F(ψ)=a(∣ψ∣2)ψF(\psi)=a(|\psi|^2)\psi with real a∈C1a\in C^1, the equation conserves the charge ∫∣ψ∣2 dx\int|\psi|^2\,dx and the energy 12∫∣ψ′∣2 dx+U(ψ(0))\frac12\int|\psi'|^2\,dx+U(\psi(0)), and admits solitary waves (standing bound states) ψ(x,t)=ψω(x)eiωt\psi(x,t)=\psi_\omega(x)e^{i\omega t}, which here take the explicit form ψω(x)=Ceiθe−ω∣x∣\psi_\omega(x)=Ce^{i\theta}e^{-\sqrt{\omega}|x|} with constant phase θ∈[0,2π]\theta\in[0,2\pi] and, for C>0C>0, amplitude and frequency related by a(C2)=2ωa(C^2)=2\sqrt{\omega}; the family parametrized by (C,θ)(C,\theta) (equivalently (ω,θ)(\omega,\theta)) forms a two-dimensional solitary manifold in the finite-energy phase space. The central question of soliton dynamics is whether a solution converges, as t→±∞t\to\pm\infty, to a single solitary wave plus dispersing radiation — a sum of another solitary wave and a solution W(t)Φ±W(t)\Phi_\pm of the free Schrödinger equation, where W(t)W(t) denotes the free Schrödinger group. Such asymptotic stability was first established by Soffer and Weinstein for weakly nonlinear multichannel scattering models and by Buslaev and Perelman for one-dimensional NLS with data near a stable soliton.

Buslaev, Komech, Kopylova and Stuart (2008) proved this for the point-coupled model: for data sufficiently close to the stable part of the solitary manifold, and under the spectral condition that the linearized operator has no discrete spectrum except the zero eigenvalue forced by the phase symmetry, the solution satisfies

ψ(⋅,t)=ψω±eiω±t+W(t)Φ±+r±(t),∥r±(t)∥Cb(R)∩L2(R)=O(∣t∣−1/2),\psi(\cdot,t)=\psi_{\omega_\pm}e^{i\omega_\pm t}+W(t)\Phi_\pm+r_\pm(t),\qquad \|r_\pm(t)\|_{C_b(\mathbb R)\cap L^2(\mathbb R)}=O(|t|^{-1/2}),

where Cb(R)C_b(\mathbb R) is the space of bounded continuous functions. The point model is simple enough that all spectral properties of the linearization can be computed exactly. Writing ν=a′(C2)C2/a(C2)\nu=a'(C^2)C^2/a(C^2) for the normalized derivative of the nonlinearity at the solitary-wave amplitude CC, the linearly stable range is ν<1\nu<1: for ν>1\nu>1 the wave is linearly unstable; for ν∈(−∞,0)∪(0,1/2)\nu\in(-\infty,0)\cup(0,1/\sqrt2) the linearization has no additional discrete spectrum (the regime of the 2008 theorem); and for ν∈(1/2,1)\nu\in(1/\sqrt2,1) it acquires a pair of purely imaginary eigenvalues ±iμ\pm i\mu with μ<ω\mu<\omega — an oscillatory internal mode, i.e. a discrete non-dispersing oscillatory component of the linearized dynamics. The paper concludes by posing two open problems: whether the above asymptotics hold for any finite-energy solution (global attraction to the solitary manifold), which the authors hoped to attack with techniques then being developed for the Klein–Gordon equation with the same point nonlinearity; and the extension to the internal-mode regime. The subsequent history is asymmetric. The internal-mode regime was largely settled for the resonant part of the range by Komech, Kopylova and Stuart (2012) via a nonlinear Fermi golden rule condition — a non-degeneracy of the coupling between the second harmonic of the internal mode and the continuous spectrum, where the golden rule is the standard mechanism by which a discrete mode is damped through resonant energy leakage into radiation. Global attraction, in contrast, was proved in the Komech program for the Klein–Gordon, wave, and Dirac analogs with concentrated nonlinearities but has resisted all attempts in the Schrödinger case, whose one-dimensional radiation decays only like t−1/2t^{-1/2}, too slowly for the local-energy-decay and tauberian arguments (methods of harmonic analysis that convert averaged decay of the solution's Fourier transform into pointwise decay in time) that drive the relativistic proofs.

Problem Statement

The overarching objective is to complete the soliton-plus-radiation asymptotics of the one-dimensional Schrödinger equation coupled to a nonlinear oscillator,

iψ˙(x,t)=−ψ′′(x,t)−δ(x) a(∣ψ(0,t)∣2) ψ(0,t),x∈R,i\dot\psi(x,t)=-\psi''(x,t)-\delta(x)\,a(|\psi(0,t)|^2)\,\psi(0,t),\qquad x\in\mathbb R,

with a U(1)U(1)-invariant C2C^2 potential, i.e. F(ψ)=a(∣ψ∣2)ψF(\psi)=a(|\psi|^2)\psi with real a∈C1a\in C^1 and U(z)≥A−B∣z∣2U(z)\ge A-B|z|^2 for some constants A∈RA\in\mathbb R and B>0B>0, in the two directions left open by Buslaev–Komech–Kopylova–Stuart (2008).

(1) Global attraction (principal part). Prove, or disprove, that for every admissible nonlinearity aa, every finite-energy solution ψ∈C(R,H1(R))\psi\in C(\mathbb R,H^1(\mathbb R)) admits the long-time asymptotics

ψ(⋅,t)=ψω±eiω±t+W(t)Φ±+r±(t),t→±∞,\psi(\cdot,t)=\psi_{\omega_\pm}e^{i\omega_\pm t}+W(t)\Phi_\pm+r_\pm(t),\qquad t\to\pm\infty,

where ψω±(x)=C±eiθ±e−ω±∣x∣\psi_{\omega_\pm}(x)=C_\pm e^{i\theta_\pm}e^{-\sqrt{\omega_\pm}|x|}, C±≥0C_\pm\ge0, θ±∈[0,2π]\theta_\pm\in[0,2\pi], is the solitary wave of frequency ω±>0\omega_\pm>0 — the constant phase θ±\theta_\pm being part of the solitary wave by the U(1)U(1) symmetry, and, for C±>0C_\pm>0, amplitude and frequency being related by a(C±2)=2ω±a(C_\pm^2)=2\sqrt{\omega_\pm} (the two-parameter family of these waves is the solitary manifold; following the source's convention, C=0C=0 gives the zero function, which is a solitary wave for every admissible nonlinearity since F(0)=0F(0)=0) — W(t)W(t) denotes the free one-dimensional Schrödinger group, the asymptotic scattering states are Φ±∈Cb(R)∩L2(R)\Phi_\pm\in C_b(\mathbb R)\cap L^2(\mathbb R), and the remainder r±(t)→0r_\pm(t)\to0 in the global norm Cb(R)∩L2(R)C_b(\mathbb R)\cap L^2(\mathbb R) (the rate O(∣t∣−1/2)O(|t|^{-1/2}) of the local theorem is expected, but the essential claim is the decomposition with vanishing remainder), under no restriction on the initial data beyond finite energy — no closeness to the solitary manifold and no smallness assumption. Under this convention, the zero solution ψ≡0\psi\equiv0, and the degenerate linear case a≡0a\equiv0 (in which every solution is a free solution), satisfy the decomposition trivially with the zero solitary wave and hence do not constitute counterexamples; a disproof requires exhibiting an admissible nonlinearity and a finite-energy solution that admits no such decomposition for any C≥0C\ge0, θ±∈[0,2π]\theta_\pm\in[0,2\pi], ω±>0\omega_\pm>0, and Φ±∈Cb(R)∩L2(R)\Phi_\pm\in C_b(\mathbb R)\cap L^2(\mathbb R).

(2) The non-resonant internal-mode regime (complementary part). Prove, or disprove, that the local asymptotic-stability theorem (data in a fixed neighborhood of a linearly stable solitary wave) extends to the remaining parameter window ν∈(2(1+3)/4, 1)\nu\in(\sqrt2(1+\sqrt3)/4,\,1), where ν=a′(C2)C2/a(C2)\nu=a'(C^2)C^2/a(C^2) is the normalized derivative of the nonlinearity at the solitary-wave amplitude CC: in this window the linearization possesses an oscillatory internal mode — a pair of purely imaginary eigenvalues ±iμ\pm i\mu — whose second harmonic 2μ2\mu lies below the continuous-spectrum edge ω\omega, so that the Fermi-golden-rule mechanism of Komech–Kopylova–Stuart (2012), the damping of the internal mode by resonant energy leakage into radiation, is unavailable; a disproof requires an admissible nonlinearity in this window together with data in a fixed neighborhood of a linearly stable solitary wave for which the local asymptotics fail; the source's announced goal was coverage of the entire linearly stable range ν<1\nu<1.

Part (1) is the principal open question; part (2) is complementary to it and can be addressed independently. Part (2) extends the parameter coverage of the local theory to the remaining non-resonant open window — the threshold values ν=0\nu=0, ν=1/2\nu=1/\sqrt2, and ν=2(1+3)/4\nu=\sqrt2(1+\sqrt3)/4, at which the corresponding spectral conditions degenerate to equalities, being borderline cases. The two parts are logically linked: an affirmative proof of part (1) supplies the qualitative decomposition but not the rate refinement that part (2) asserts, whereas a non-resonant-window counterexample — an admissible nonlinearity in the window together with data near a linearly stable solitary wave for which the local asymptotics fail — refutes both parts, its data lying within part (1)'s scope. The two parts share one scientific objective: the complete soliton-radiation description of all finite-energy dynamics of the model over its full linearly stable parameter range.

References

  1. V. S. Buslaev, A. I. Komech, E. A. Kopylova, D. M. A. Stuart, On Asymptotic Stability of Solitary Waves in Schrödinger Equation Coupled to Nonlinear Oscillator, Communications in Partial Differential Equations 33 (2008), no. 4, 669–705, DOI 10.1080/03605300801970937, https://doi.org/10.1080/03605300801970937; preprint arXiv:math-ph/0702013, https://arxiv.org/abs/math-ph/0702013.
  2. A. I. Komech, E. A. Kopylova, D. M. A. Stuart, On asymptotic stability of solitons in a nonlinear Schrödinger equation, Communications on Pure and Applied Analysis 11 (2012), no. 3, 1063–1079, DOI 10.3934/cpaa.2012.11.1063, https://doi.org/10.3934/cpaa.2012.11.1063; preprint arXiv:0807.1878, https://arxiv.org/abs/0807.1878.
  3. A. Comech, A. Komech, E. Kopylova, Attractors of Hamiltonian nonlinear partial differential equations, preprint arXiv:2212.14152 (2022), https://arxiv.org/abs/2212.14152 (survey version of the monograph Attractors of Hamiltonian Nonlinear Partial Differential Equations, Cambridge University Press, 2021).
  4. A. I. Komech, A. A. Komech, On the global attraction to solitary waves for the Klein–Gordon equation coupled to a nonlinear oscillator, Comptes Rendus. Mathématique 343 (2006), no. 2, 111–114, DOI 10.1016/j.crma.2006.06.009, https://doi.org/10.1016/j.crma.2006.06.009.
  5. A. I. Komech, A. A. Komech, Global well-posedness for the Schrödinger equation coupled to a nonlinear oscillator, Russian Journal of Mathematical Physics 14 (2007), no. 2, 164–173, DOI 10.1134/S1061920807020057, https://doi.org/10.1134/S1061920807020057.
  6. E. A. Kopylova, On global attraction to solitary waves for the Klein–Gordon equation with concentrated nonlinearity, Nonlinearity 30 (2017), no. 11, 4191–4207, DOI 10.1088/1361-6544/aa84bf, https://doi.org/10.1088/1361-6544/aa84bf.
  7. E. A. Kopylova, Soliton type asymptotics for the Klein-Gordon equation coupled to nonlinear oscillator, Nonlinear Differential Equations and Applications NoDEA 32 (2025), no. 4, DOI 10.1007/s00030-025-01078-y, https://doi.org/10.1007/s00030-025-01078-y.
  8. R. Adami, D. Noja, C. Ortoleva, Orbital and asymptotic stability for standing waves of a nonlinear Schrödinger equation with concentrated nonlinearity in dimension three, Journal of Mathematical Physics 54 (2013), no. 1, 013501, DOI 10.1063/1.4772490, https://doi.org/10.1063/1.4772490.
  9. R. Adami, D. Noja, C. Ortoleva, Asymptotic stability for standing waves of a NLS equation with subcritical concentrated nonlinearity in dimension three: Neutral modes, Discrete and Continuous Dynamical Systems 36 (2016), no. 11, 5837–5879, DOI 10.3934/dcds.2016057, https://doi.org/10.3934/dcds.2016057.
  10. C. Collot, P. Germain, Asymptotic stability of solitary waves for one-dimensional nonlinear Schrödinger equations, Journal of the European Mathematical Society (2025), DOI 10.4171/jems/1642, https://doi.org/10.4171/jems/1642.
  11. V. S. Buslaev, G. S. Perelman, On the stability of solitary waves for nonlinear Schrödinger equations, American Mathematical Society Translations: Series 2, 164 (1995), 75–98, DOI 10.1090/trans2/164/04, https://doi.org/10.1090/trans2/164/04.
  12. A. Soffer, M. I. Weinstein, Resonances, radiation damping and instability in Hamiltonian nonlinear wave equations, Inventiones Mathematicae 136 (1999), no. 1, 9–74, DOI 10.1007/s002220050303, https://doi.org/10.1007/s002220050303.