- Research
- Open access
- Published:
Asymptotic stability of standing waves for the coupled nonlinear Schrödinger system
Boundary Value Problems volume 2015, Article number: 3 (2015)
Abstract
This paper considers the asymptotic stability of standing waves of a coupled nonlinear Schrödinger system with an attractive potential. Meanwhile, the existence of a center manifold is obtained.
1 Introduction and main results
The two-component system of time-dependent nonlinear Schrödinger equations arises in the binary mixture of Bose-Einstein condensates with two different hyperfine states (see [1]):
where \(u(x,t)\) and \(v(x,t)\) denote the wave functions, ħ is the Planck constant divided by 2π, m is atom mass, \(V_{i}\) is the trapping potential for the ith hyperfine state, \(a_{i}\geq0\), \(i=1,2,3,4\), is the associated axial frequency.
By rescaling and some simple assumptions, the above system could be viewed as the coupled nonlinear Schrödinger system:
with initial data
In the present paper, we assume that \(a_{2}=a_{4}\).
In the last decades, there has been a lot of interest in the study of localized modes in the coupled nonlinear Schrödinger system. In particular, for the existence of standing wave (periodic in time and exponentially localized in space) for various coupled nonlinear Schrödinger system, one may refer to [2–6] for more details. As we know, there are few results on the asymptotic stability of standing waves for the coupled nonlinear Schrödinger system (1.1)-(1.2).
Assume that the standing wave of system (1.1)-(1.2) has the form as follows:
where \(E_{1},E_{2}\in\mathbb{R}\) and \(\psi_{E_{1}},\phi_{E_{2}}\in\mathbf {H}^{2}(\mathbb{R}^{2})\) satisfy the time independent equations:
Here, we give some notations. Assume that \(\langle x\rangle=(1+|x|^{2})^{\frac {1}{2}}\). For \(\sigma\in\mathbb{R}\), the Sobolev space \(L_{\sigma }^{2}\) with the norm \(\|f(x)\|_{L^{2}_{\sigma}}=\|\langle x\rangle ^{\sigma }f(x)\|_{L^{2}}\) denotes the space of functions \(f(x)\) such that \(\langle x\rangle ^{\sigma}f(x)\) are square integrable. \(C_{a,b,\ldots}\) denotes a constant depending on \(a,b,\ldots\) .
Define
and
Then the system (1.3)-(1.4) can be rewritten as
It is obvious that \(X\equiv0\) is a trivial solution of (1.5). In order to make a bifurcation with a nontrivial, one parameter family of solutions, we need (iii) in the following assumptions:
-
(A)
Assume that:
-
(i)
There exist C and α such that
$$\begin{aligned} \bigl|V_{i}(x)\bigr|\leq C\langle x\rangle ^{-\alpha},\quad i=1,2, \forall x\in\mathbb{R}^{2}. \end{aligned}$$Note that potential \(V_{1}\), \(V_{2}\) could be different in our paper.
-
(ii)
0 is a regular point of the spectrum of the linear operator \(-\triangle+V_{i}\) acting on \(\mathbf{L}^{2}\), \(i=1,2\).
-
(iii)
For \(i=1,2\), \(-\triangle+V_{i}\) acting on \(\mathbf{L}^{2}\) has exactly two negative eigenvalues \(\omega_{i}\) with corresponding normalized eigenvectors \(\psi_{0}\) and \(\phi_{0}\), respectively. It is well known that \(\psi_{0}\) and \(\phi_{0}\) are exponentially decaying as \(|x|\rightarrow \infty\), and they could be chosen strictly positive.
-
(i)
The main result of this paper is the following.
Theorem 1.1
Assume that hypothesis (A) holds. Then there exists an \(\varepsilon_{0}>0\) such that for any fixed \(\sigma>1\) and the initial data \(u_{0}\), \(v_{0}\) satisfying
the initial value problem (1.1)-(1.2) is globally well-posed in \(\mathbf{H}^{1}\times\mathbf{H}^{1}\).
Moreover, the solution of the initial value problem (1.1)-(1.2) has the form
with
where \(2\leq p\leq p_{0}\).
This paper is organized as follows. In the next section, we show the existence of a center manifold of the coupled nonlinear Schrödinger system (1.1)-(1.2). We investigate the asymptotic stability of the system (1.1)-(1.2) in Section 3.
2 The center manifold
This section is devoted to the proof of the existence of a center manifold. The method of constructing the center manifold is based on the standard bifurcation argument in Banach spaces for (1.3)-(1.4) at \((\omega_{1},\omega_{2})^{T}\) (see [7]). Since the spectrum of the operator \(-\triangle+V_{i}\) has a discrete and continuous part, we follow the idea of [8, 9] and decompose the solution of (1.5) in its projection onto the discrete and continuous part:
where \(P_{c}\) denotes the projector onto the continuous spectrum of \(-\triangle+V_{i}\) in \(L^{2}\).
Now, projecting (1.5) onto \(X_{0}=(\psi_{0},\phi_{0})^{T}\) and its orthogonal complement, i.e. range \(P_{c}\), we have
where
and
Now we have the following result on the center manifold of the system (1.1)-(1.2):
Theorem 2.1
Let \(\delta_{E_{1}},\delta_{E_{2}},\delta>0\). Then for \(|E_{1}-\omega _{1}|<\delta_{E_{1}}\), \(|E_{2}-\omega_{2}|<\delta_{E_{2}}\), \(\|\psi_{E_{1}}\| _{\mathbf{L}^{2}_{\sigma}\cap\mathbf{H}^{2}}<\delta\), and \(\|\phi_{E_{2}}\| _{\mathbf{L}^{2}_{\sigma}\cap\mathbf{H}^{2}}<\delta\), there exist \(\mathbf{C}^{1}\) functions
such that the eigenvalue problem (1.3)-(1.4) has a unique solution up to multiplication with \(e^{i\theta}\), \(\theta\in(0,2\pi)\), which can be represented as
Proof
We follow the idea of [8, 9] and use the implicit function theorem to solve (2.1)-(2.2). We define a \(\mathbf{C}^{1}\) function \(\mathcal{F}:(-\infty,0)^{2}\times \mathbb{C}^{2}\times(\mathbf{L}^{2}_{\sigma}\cap\mathbf{H}^{2})^{2}\) as
It is easy to see that
Therefore, from the implicit function theorem, we know that (2.1) has a unique solution \(X=\bar{X}(\mathcal{E},a,b)\) with \(\|X\| _{(\mathbf{L}^{2}_{\sigma}\cap\mathbf{H}^{2})^{2}}<\delta_{1}\), where \(\bar{X}(\mathcal{E},a,b)\) is a \(\mathbf{C}^{1}\) function from \((\omega _{1}-\delta_{1},\omega_{1}+\delta_{1})\times(\omega_{2}-\delta_{1},\omega _{2}+\delta_{1})\times\{(a,b)\in\mathbb{C}^{2}:|a|,|b|<\delta_{1}\}\) to \((\mathbf{L}^{2}_{\sigma}\cap\mathbf{H}^{2})^{2}\) for \(\delta_{1}>0\). By direct computation, we know that \((e^{i\theta}a,e^{i\theta}b,e^{i\theta}X)\) is also a solution of (2.1) for \(\theta\in(0,2\pi)\). By uniqueness, we have \(\bar{X}(\mathcal{E},a,b)=\tilde{\mathcal{A}}\bar{X}(\mathcal {E},|a|,|b|)\), where \(\tilde{\mathcal{A}}=\operatorname{diag} (\frac{a}{|a|},\frac {b}{|b|} )\).
Substituting \(X=\bar{X}(\mathcal{E},a,b)\) into (2.2). From \(\bar {X}(\mathcal{E},a,b)=\tilde{\mathcal{A}}\bar{X}(\mathcal{E},|a|,|b|)\), we have
where
Define a \(\mathbf{C}^{1}\) function \(\mathcal{G}\) on \((\omega_{1}-\delta _{1},\omega_{1}+\delta_{1})\times(\omega_{2}-\delta_{1},\omega_{2}+\delta _{1})\times(-\delta_{1},\delta_{1})^{2}\) as follows:
We can see that
From the implicit function theorem, we find that (2.2) has a unique solution \(\mathcal{E}=\tilde{\mathcal{E}}(|a|,|b|)\) with \(X=\bar {X}(\mathcal{E},a,b)\), where \(\tilde{\mathcal{E}}(|a|,|b|)\) is a \(\mathbf{C}^{1}\) function from \((-\delta,\delta)\times(-\delta,\delta)\) to \((\omega_{1}-\delta _{E_{1}},\omega_{1}+\delta_{E_{1}})\times(\omega_{2}-\delta_{E_{2}},\omega _{2}+\delta_{E_{2}})\) for \(0<\delta\leq\delta_{1}\) and \(0<\delta _{E_{1}},\delta_{E_{2}}\leq\delta_{1}\). □
3 Proof of Theorem 1.1
We consider the asymptotic stability of the coupled Schrödinger system (1.1)-(1.2) in this section. The global well-posedness for the system (1.1)-(1.2) with small initial data is obtained by Cazenave [10].
Define
By Theorem 2.1, we can define \(h_{1}(a(t))\) and \(h_{2}(b(t))\in\mathbb{C}\) by choosing \(\varepsilon_{0}<\delta\). Therefore,
We can see that the solution is described by the scalar \(a(t),b(t)\in \mathbb{C}\) and \(r_{1}(t),r_{2}(t)\in\mathbf{C}(\mathbb{R},\mathbf{H}^{1})\). More precisely, from (1.6)-(1.7), we have
Note that \(h_{1}\), \(h_{2}\), and \(Dh_{1}\), \(Dh_{2}\) have range orthogonal to \(\psi_{0}\), \(\phi_{0}\). \(r_{1}\), \(r_{2}\) and \(\frac {dr_{1}}{dt}\), \(\frac{dr_{2}}{dt}\) are orthogonal to \(\psi_{0}\), \(\phi_{0}\), respectively. Therefore, from (3.1)-(3.2), we get
Using \(Dh_{1}|_{a}(iE_{1}a)=-E_{1}ih_{1}(a)\) and \(Dh_{2}|_{b}(iE_{2}b)=-E_{2}ih_{2}(b)\), we have
The linear part of (3.5)-(3.6) is
Define the operator \(\mathcal{S}(t,s)Y\) as the solution of the linear equation (3.7)-(3.8):
The following lemma can be obtained by a small modification of the proof of Theorem 4.1 in [8], so we omit the proof.
Lemma 3.1
There exists \(\varepsilon>0\) such that if
then there exist \(C, C_{p}>0\) with the property that for any \(t,s\in \mathbb{R}\), we have
where \(p\geq2\) and \(\frac{1}{p'}+\frac{1}{p}=1\).
By Duhamel’s principle, it follows from (3.5)-(3.6) that
where \(R=(r_{1},r_{2})^{T}\), \(G_{1}=(g_{1},g_{2})^{T}\), \(G_{2}=(g_{3},g_{4})^{T}\), and \(DH=\operatorname{diag}(Dh_{1}|_{a},Dh_{2}|_{b})\) with
For fixed \(p\geq6\), we define the Banach space
endowed with the norm
Consider the nonlinear part in (3.12):
Lemma 3.2
-
(1)
For \(R\in\mathbb{B}\times\mathbb{B}\), the nonlinear operator \(\mathcal{N}:\mathbb{B}\times\mathbb{B}\rightarrow\mathbb{B}\times \mathbb{B}\) is well defined.
-
(2)
We have
$$\begin{aligned} \|\mathcal{N}R_{1}-\mathcal{N}R_{2}\|_{\mathbb{B}\times\mathbb{B}} \leq &C_{a_{1},a_{2},a_{3},a_{4},p}\bigl(\|R_{1}\|_{\mathbb{B}\times\mathbb{B}}+\|R_{2}\| _{\mathbb{B}\times\mathbb{B}}+\|R_{1}\|_{\mathbb{B}\times\mathbb{B}}^{2}+\| R_{2}\|_{\mathbb{B}\times\mathbb{B}}^{2}\bigr) \\ &{}\times\bigl(\|R_{1}-R_{2}\|_{\mathbb{B}\times\mathbb{B}}\bigr). \end{aligned}$$
Proof
It is obvious that the part (1) can be obtained by part (2) choosing \(R_{2}=0\). Now, we only need to prove (2). This proof is based on these estimates of \(\mathcal{N}R_{1}-\mathcal{N}R_{2}\): the \(\mathbf{L}_{-\sigma }^{2}\times\mathbf{L}_{-\sigma}^{2}\) estimate, the \(\mathbf{L}^{p}\times \mathbf{L}^{p}\) estimate and the \(\mathbf{L}^{2}\times\mathbf{L}^{2}\) estimate. In fact, with a similar argument of getting the \(\mathbf {L}_{-\sigma}^{2}\times\mathbf{L}_{-\sigma}^{2}\) estimate, we could obtain the \(\mathbf{L}^{p}\times\mathbf{L}^{p}\) estimate and the \(\mathbf {L}^{2}\times\mathbf{L}^{2}\) estimate from (3.10) and (3.11). Here, we consider the \(\mathbf{L}_{-\sigma}^{2}\times\mathbf{L}_{-\sigma }^{2}\) estimate.
Let \(R_{1}=(r_{1},r_{2})^{T}, R_{2}=(\tilde{r}_{1},\tilde{r}_{2})^{T}\in\mathbb{B}\times \mathbb{B}\). One obtains
where \(G_{1}-\tilde{G}_{1}=(g_{1}-\tilde{g}_{1},g_{2}-\tilde{g}_{2})^{T}\) and \(G_{2}-\tilde{G}_{2}=(g_{3}-\tilde{g}_{3},g_{4}-\tilde{g}_{4})^{T}\) with
Let \(4< p<\infty\) and \(\frac{1}{p'}+\frac{1}{p}=1\). We have
where \(A=(A_{1},A_{2})^{T}\) and \(B=(B_{1},B_{2})^{T}\) with
In what follows, we consider the integral terms on A, B, and \(|G_{2}-\tilde{G}_{2}|\) in (3.23). We first estimate the integral term on A. For \(\frac{1}{\alpha}+\frac{2}{p}=\frac{1}{2}\), we have
From (3.9) in Lemma 3.1, we have
where the constants
and
Using the interpolate inequality \(\|r\|_{\mathbf{L}^{\alpha}}\leq\|r\| _{\mathbf{L}^{2}}^{1-b}\|r\|^{b}_{\mathbf{L}^{p}}\) for \(\frac{1}{\alpha}=\frac {1-b}{2}+\frac{b}{p}\) with \(p\geq4\) and \(2\leq\alpha\leq p\), one obtains
where \(\frac{2}{\alpha}+\frac{1}{p}=\frac{1}{p'}\) and \(\frac {2(1-b)}{2}+\frac{2b}{p}+\frac{1}{p}=\frac{1}{p'}\) with \((1-\frac {2}{p})(1+2b)=1+\frac{2}{p}>1\).
Note that \((1-\frac{2}{p})(1+2b)>1\). Therefore, from (3.10), we have
where the constant
It is easy to obtain
for \(\frac{1}{\alpha}+\frac{2}{p}=\frac{1}{2}\). From (3.9) in Lemma 3.1, we get
where we have the constant
Hence, we conclude that
This completes the proof. □
Proof of Theorem 1.1
Define a closed ball \(B(R_{0},\mathbf {r})\subset\mathbb{B}\times\mathbb{B}\) with center \(R_{0}=\mathcal {S}(t,0)R(0)\) and radius \(\mathbf{r}=\frac{L\|R_{0}\|_{\mathbb{B}\times \mathbb{B}}}{2-\operatorname{Lip}}\). By Lemma 3.1, there exists a constant \(C_{6}\) such that
Choose \(\varepsilon_{0}\) such that \(C_{6}\varepsilon_{0}<\frac{1}{2}(\sqrt {1+2C^{-1}_{a_{1},a_{2},a_{3},a_{4},p}}-1)\). Then there exists a constant \(0<\operatorname{Lip}<1\) such that
It is easy to conclude that the right hand side of (3.12) leaves \(B(R_{0},\mathbf{r})\) invariant, and it is a contraction with Lipschitz constant Lip on \(B(R_{0},\mathbf{r})\). From the contraction mapping theorem, (3.12) has a unique solution in \(B(R_{0},\mathbf{r})\). If we have two solutions of (3.12), one in \(\mathbf{C}(\mathbb {R},\mathbf{H}^{1}\times\mathbf{H}^{1})\) from classical well-posedness theory and one in \(\mathbf{C}(\mathbb{R},(\mathbf{L}^{2}_{-\sigma}\cap \mathbf{L}^{2}\cap \mathbf{L}^{p})\times(\mathbf{L}^{2}_{-\sigma}\cap\mathbf{L}^{2}\cap\mathbf {L}^{p}))\) from the above argument for \(p\geq6\). By uniqueness and the continuous embedding of \(\mathbf{H}^{1}\) in \(\mathbf{L}^{2}_{-\sigma}\cap \mathbf{L}^{2}\cap\mathbf{L}^{p}\), we infer that the two solutions must coincide. Therefore, the time decaying estimates also hold for the \(\mathbf{H}^{1}\times\mathbf{H}^{1}\) solutions. This completes the proof of Theorem 1.1. □
References
Pitaevskii, L, Stingari, S: Bose-Einstein Condensation. Oxford University Press, Oxford (2003)
Ambrosetti, A, Colorado, E: Standing waves of some coupled nonlinear Schrödinger equations. J. Lond. Math. Soc. 75, 67 (2007)
Belmonte-Beitia, VM, Pérez-García, J, Torres, PJ: Solitary waves for linearly coupled nonlinear Schrödinger equations with inhomogeneous coefficients. J. Nonlinear Sci. 19, 437 (2009)
Hioe, FT, Carroll, CE: Coupled Gross-Pitaevskii equations with a square-well and a Coulomb potential. Phys. Lett. A 299, 189 (2002)
Roskes, GJ: Some nonlinear multiphase interactions. Stud. Appl. Math. 55, 231 (1976)
Yang, J: Multiple permanent-wave trains in nonlinear systems. Stud. Appl. Math. 100, 127 (1998)
Nirenberg, L: Topics in Nonlinear Functional Analysis. Chapter 6 by E Zehnder. Notes by RA Artino. Revised reprint of the 1974 original. Courant Lecture Notes in Mathematics, vol. 6. New York University, Courant Institute of Mathematical Sciences, New York (2001)
Kirr, E, Zarnescu, A: On the asymptotic stability of bound states in 2D cubic Schrödinger equation. Commun. Math. Phys. 272, 443 (2007)
Pillet, CA, Wayne, CE: Invariant manifolds for a class of dispersive, Hamiltonian partial differential equations. J. Differ. Equ. 141, 310 (1997)
Cazenave, T: Semilinear Schrödinger Equations. Courant Lecture Notes in Mathematics, vol. 10. New York University, Courant Institute of Mathematical Sciences, New York (2003)
Acknowledgements
The authors express their sincere thanks to the anonymous referees for the comments.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.
Rights and permissions
Open Access This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.
About this article
Cite this article
Liao, Y., Sun, Q., Zhao, X. et al. Asymptotic stability of standing waves for the coupled nonlinear Schrödinger system. Bound Value Probl 2015, 3 (2015). https://doi.org/10.1186/s13661-014-0266-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-014-0266-4