**On asymptotic stability of moving kink for relativistic Ginzburg-Landau equation** E. A. Kopylova 111Supported partly by the grants of DFG, FWF and RFBR. *Institute for Information Transmission Problems RAS B.Karetnyi 19, Moscow 101447,GSP-4, Russia *e-mail: elena.kopylova@univie.ac.at A. I. Komech${}^{\mbox{\scriptsize 1,}}\!\!$ 222Supported partly by the Alexander von Humboldt Research Award. *Fakultät für Mathematik, Universität Wien and Institute for Information Transmission Problems RAS *e-mail: alexander.komech@univie.ac.at We prove the asymptotic stability of the moving kinks for the nonlinear relativistic wave equations in one space dimension with a Ginzburg-Landau potential: starting in a small neighborhood of the kink, the solution, asymptotically in time, is the sum of a uniformly moving kink and dispersive part described by the free Klein-Gordon equation. The remainder decays in a global energy norm. Crucial role in the proofs play our recent results on the weighted energy decay for the Klein-Gordon equations. *Keywords*: Relativistic nonlinear wave equation, asymptotic stability, kink, weighted energy decay, symplectic projection, modulation equations. *2000 Mathematics Subject Classification*: 35Q51, 37K40. ## 1 Introduction There has been wide spread interest in the dynamics of topological excitations of classical relativistic field theories [[2](#bib.bib2), [3](#bib.bib3)]. These excitations are finite energy solutions which do not decay to one of the true ground states because of topological constraints, said differently, these excitations are separated by an infinitely high potential barrier from the ground state. In our contribution we will study in mathematical detail one of the simplest examples. The field, $\psi$, is real valued and defined on the line, $\psi:{\mathbb{R}}\to{\mathbb{R}}$. The Hamiltonian function reads [TABLE] with $\pi$ the momentum canonically conjugate to $\psi$ and a smooth potential $U$. This leads to the equation of motion [TABLE] where $F(\psi)=-U^{\prime}(\psi)$. For an introduction, let us consider the Ginzburg-Landau quartic double well potential of the form $U(\psi)=(\psi^{2}-a^{2})^{2}/(4a^{2})$. Then the topological excitations are defined through ${\cal H}<\infty$ and the boundary conditions [TABLE] with a fixed $a>0$. Amongst them there are soliton-like solutions which travel with constant velocity, [TABLE] where $\gamma=1/\sqrt{1-v^{2}}$ is the Lorentz contraction. The solitons ([1.6](#S1.E6)) are related by a Lorentz boost, since equation ([1.2](#S1.E2)) is relativistically invariant. We will consider more general double well potentials for which [TABLE] and [TABLE] similarly to the quartic potential. In this case the soliton-like solutions also exist, [TABLE] where $s(\cdot)$ is a ”kink” solution to the corresponding stationary equation [TABLE] In general our goal is to clarify the special role of the soliton-like solutions ([1.6](#S1.E6)) as long time asymptotics for any finite energy topological excitations satisfying ([1.3](#S1.E3)). Namely, if one chooses some arbitrary finite energy initial state satisfying ([1.3](#S1.E3)), one would expect that for $t\to\infty$ the solution separates into two pieces: one piece is a finite collection of travelling solitons of the form ([1.6](#S1.E6)) and their negatives with some velocities $v_{j}\in(-1,1)$ and the shifts $q_{j}$ depending in a complicated way on the initial data, and the second radiative piece which is a dispersive solution to the free Klein-Gordon equation which propagates to infinity with the velocity $1$. Our aim here is to elucidate this general picture by mathematical arguments for initial data sufficiently close to a soliton ([1.6](#S1.E6)). Let us discuss our choice of the smooth potentials $U$. The condition ([1.5](#S1.E5)) is necessary and sufficient for the existence of a finite energy static solution $s(x)$ to ([1.7](#S1.E7)) when ([1.4](#S1.E4)) holds. Indeed, the condition is obviously sufficient. On the other hand, the ”energy conservation” [TABLE] and $s(\pm\infty)=\pm a$ imply that $E=0$. Therefore, $U(\psi)>0$ for $\psi\in(-a,a)$ since otherwise the boundary conditions $s(\pm\infty)=\pm a$ would fail. As a byproduct, our kink solution is monotone increasing, and [TABLE] Let us note that only the behavior of $U$ near the interval $[-a,a]$ is of importance since the solution is expected to be close to a soliton. However, we will assume additionally the potential to be bounded from below [TABLE] to have a well posed Cauchy problem for all finite energy initial states. Summarising, we formulate our first basic condition on the potential, for technical reasons adding a flatness condition. **Condition U1**. *The potential $U$ is a real smooth function which satisfies ([1.4](#S1.E4)), ( [1.5](#S1.E5)), ([1.10](#S1.E10)), and the following condition holds with some $m>0$* [TABLE] Let us comment on the condition ([1.11](#S1.E11)) (see also Remark [4.11](#S4.Thmtheorem11)). First, the condition means that $U^{\prime\prime}(-a)=U^{\prime\prime}(a)$ though we do not need the potential to be reflection symmetric. We consider the solutions close to the kink, $\psi(x,t)=s(\gamma(x-vt-q))+\phi(x,t)$, with small perturbations $\phi(x,t)$. For such solution the condition ([1.11](#S1.E11)) and the asymptotics ([1.3](#S1.E3)) mean that the equation ([1.2](#S1.E2)) is almost linear Klein-Gordon equation for large $|x|$ which is helpful for application of the dispersive properties. Finally, we expect that the degree $14$ in ([1.11](#S1.E11)) is technical, and a smaller degree should be sufficent. Let us note that a similar condition has been introduced in [[4](#bib.bib4), [5](#bib.bib5)] in the context of the Schrödinger equation. Further we need some assumptions on the spectrum of the linearised equation. Let us rewrite the equation ([1.2](#S1.E2)) in the vector form, [TABLE] Now the soliton-like solutions ([1.6](#S1.E6)) become [TABLE] for $q,v\in{\mathbb{R}}$ with $|v|<1$, where [TABLE] The states $S_{q,v}:=Y_{q,v}(0)$ form the solitary manifold [TABLE] The linearized operator near the soliton solution $Y_{q,v}(t)$ is (see Section [4](#S4), formula ([4.20](#S4.E20))) [TABLE] where [TABLE] By ([1.21](#S1.E21)), we have [TABLE] since [TABLE] by the condition **U1**. In Section [4](#S4) we show that the spectral properties of the operator $A_{v}$ are determined by the corresponding properties of its determinant, which is the Schrödinger operator [TABLE] The spectral properties of $H_{v}$ are identical for all $v\in(-1,1)$ since the relation $V_{v}(x)=V_{0}(\gamma x)$ implies [TABLE] This equivalence manifests the relativistic invariance of the equation ([1.15](#S1.E15)). The continuous spectrum of the operator $H_{v}$ coincides with $[m^{2},\infty)$. The point [math] belongs to the discrete spectrum with corresponding eigenfunction $\psi^{\prime}_{v}$. By ([1.17](#S1.E17)) and ([1.9](#S1.E9)) we have $\psi^{\prime}_{v}(x)=\gamma s^{\prime}(\gamma x)>0$ for $x\in{\mathbb{R}}$. Hence, $\psi^{\prime}_{v}$ is the groundstate, and all remaining discrete spectrum is contained in $(0,m^{2}]$. For $\alpha\in{\mathbb{R}}$, $p\geq 1$, and $l=0,1,2,...$ let us denote by $W^{l,p}_{\alpha}$, the weighted Sobolev space of the functions with the finite norm [TABLE] and $H^{l}_{\alpha}:=W^{l,2}_{\alpha}$, so $L^{2}_{\alpha}:=H^{0}_{\alpha}$ are the Agmon’s weighted spaces. ###### **Definition 1.1****.** *(cf. [[9](#bib.bib9), [16](#bib.bib16)]) A nonzero solution $\psi\in L^{2}_{-1/2-0}({\mathbb{R}})\setminus L^{2}({\mathbb{R}})$ to $H_{v}\psi=m^{2}\psi$ is called a resonance.* Now we can formulate our second basic condition on the potential. **Condition U2**. *For any $v\in(-1,1)$ i) [math] is only eigenvalue of $H_{v}$. ii) $m^{2}$ is not a resonance of $H_{v}$. * We show that Condition **U2** implies the boundedness of the resolvent of the operator $A_{v}$ in the corresponding weighted Agmon spaces at the edge points $\pm im/\gamma$ of its continuous spectrum. Both conditions **U1**, **U2** can be satisfied though it is non-obvious. Let us note that the quartic Ginzburg-Landau potential does not satisfy nor ([1.11](#S1.E11)) nor the conditions **U2** i) nor **U2** ii). We will prove elsewhere that the corresponding examples of potentials satisfying both **U1** and **U2** can be constructed as smoothened piece-wise quadratic potentials. We now can formulate the main result of our paper. Namely, we will prove the following asymptotics [TABLE] for solutions to ([1.15](#S1.E15)) with initial states close to a soliton-like solution ([1.6](#S1.E6)). Here $W_{0}(t)$ is the dynamical group of the free Klein-Gordon equation, ${\Phi}_{\pm}$ are the corresponding asymptotic states, and the remainder converges to zero $\sim t^{-1/2}$ in the global energy norm of the Sobolev space $H^{1}({\mathbb{R}})\oplus L^{2}({\mathbb{R}})$. Let us comment on previous results in this field. $\bullet$ *Orbital stability of the kinks* For 1D relativistic nonlinear Ginzburg-Landau equations ([1.2](#S1.E2)) the orbital stability of the kinks has been proved in [[10](#bib.bib10)]. $\bullet$ *The Schrödinger equation* The asymptotics of type ([1.24](#S1.E24)) were established for the first time by Soffer and Weinstein [[23](#bib.bib23), [24](#bib.bib24)] (see also [[19](#bib.bib19)]) for nonlinear $U(1)$-invariant Schrödinger equation with a potential for small initial states and sufficiently small nonlinear coupling constant. The results have been extended by Buslaev and Perelman [[4](#bib.bib4)] to the translation invariant 1D nonlinear $U(1)$-invariant Schrödinger equation. The novel techniques [[4](#bib.bib4)] are based on the ”separation of variables” along the solitary manifold and in the transversal directions. The symplectic projection allows to exclude the unstable directions corresponding to the zero discrete spectrum of the linearized dynamics. Similar techniques were developed by Miller, Pego and Weinstein for the 1D modified KdV and RLW equations, [[17](#bib.bib17), [18](#bib.bib18)].The extensions to higher dimensions were obtained in [[6](#bib.bib6), [12](#bib.bib12), [22](#bib.bib22), [27](#bib.bib27)]. $\bullet$ *Nonrelativistic Klein-Gordon equations* The asymptotics of type ([1.24](#S1.E24)) were extended to the nonlinear 3D Klein-Gordon equations with a potential [[25](#bib.bib25)], and for translation invariant system of the 3D Klein-Gordon equation coupled to a particle [[11](#bib.bib11)]. $\bullet$ *Wave front of 3D Ginzburg-Landau equation* The asymptotic stability of wave front was proved for 3D relativistic Ginzburg-Landau equation with initial data which differ from the wave front on a compact set [[7](#bib.bib7)]. The wave front is the solution which depends on one space variable only, so it is not a soliton. The equation differs from the 1D equation ([1.2](#S1.E2)) by the additional 2D Laplacian. The additional Laplacian improves the dispersive decay for the corresponding linearized Klein-Gordon equation in the continuous spectral space that provides the needed decay for the transversal dynamics. The proving of the asymptotic stability of the solitons and kinks for relativistic equations remained an open problem till now. The investigation crucially depends on the spectral properties for linearized equation which are completely unknown for higher dimensions. For 1D case main obstacle was the slow decay $\sim t^{-1/2}$ for the free 1D Klein-Gordon equation (see the discussion in [[7](#bib.bib7), Introduction]). Let us comment on our approach. We follow general strategy of [[4](#bib.bib4), [5](#bib.bib5), [6](#bib.bib6), [7](#bib.bib7), [11](#bib.bib11), [25](#bib.bib25)]: symplectic projection onto the solitary manifold, modulation equations, linearization of the transversal equations and further Taylor expansion of the nonlinearity etc. We develop for relativistic equations general scheme which is common in almost all papers in this area: dispersive estimates for high energy and low energy components of the solutions to linearized equation, virial and $L^{1}-L^{\infty}$ estimates and method of majorants. However, the corresponding statements and their proofs in the context of relativistic equations are completely new. Let us comment on our novel techniques. I. The decay $\sim t^{-3/2}$ from Theorem [4.8](#S4.Thmtheorem8) for the linearized transversal dynamics relies on our novel approach [[13](#bib.bib13), [14](#bib.bib14)] to 1D Klein-Gordon equation. II. The novel ”virial type” estimate ([4.43](#S4.E43)) is relativistic version of the bound [[5](#bib.bib5), (1.2.5)] used in [[5](#bib.bib5)] in the context of the nonlinear Schrödinger equation (see Remark [4.11](#S4.Thmtheorem11)). III. We establish an appropriate relativistic version ([4.32](#S4.E32)) of $L^{1}\to L^{\infty}$ estimates. Both estimates ([4.43](#S4.E43)) and ([4.32](#S4.E32)) play crucial role in obtaining the bounds for the majorants. IV. Finally, we give the complete proof of the soliton asymptotics ([1.24](#S1.E24)). In the context of the Schrödinger equation, the proof of the corresponding asymptotics were sketched in [[5](#bib.bib5)]. Our paper is organized as follows. In Section [2](#S2) we formulate the main theorem. In Section [3](#S3) we introduce the symplectic projection onto the solitary manifold. The linearized equation is defined in Section [4](#S4). In Section [5](#S5) we split the dynamics in two components: along the solitary manifold and in the transversal directions. In Section [6](#S6) the modulation equations for the parameters of the soliton are displayed. The time decay of the transversal component is established in Sections [7](#S7)-[11](#S11). Finally, in Section [12](#S12) we obtain the soliton asymptotics ([1.24](#S1.E24)). ## 2 Main results ### 2.1 Existence of dynamics We consider the Cauchy problem for the Hamilton system ([1.15](#S1.E15)) which we write as [TABLE] Here $Y(t)=(\psi(t),\pi(t))$, $Y_{0}=(\psi_{0},\pi_{0})$, and all derivatives are understood in the sense of distributions. To formulate our results precisely, let us first introduce a suitable phase space for the Cauchy problem ([2.1](#S2.E1)). ###### **Definition 2.1****.** *i) $E_{\alpha}:=H^{1}_{\alpha}\oplus L^{2}_{\alpha}$ is the space of the states $Y=(\psi,\pi)$ with finite norm* [TABLE] *ii) The phase space ${\cal E}:={{\cal S}}+E$, where $E=E_{0}$ and ${\cal S}$ is defined in ([1.18](#S1.E18)). The metric in ${\cal E}$ is defined as* [TABLE] *iii) $W:=W^{2,1}_{0}\oplus W^{1,1}_{0}$ is the space of the states $Y=(\psi,\pi)$ with the finite norm* [TABLE] Obviously, the Hamilton functional ([1.1](#S1.E1)) is continuous on the phase space ${\cal E}$. The existence and uniqueness of the solutions to the Cauchy problem ([2.1](#S2.E1)) follows by methods [[15](#bib.bib15), [20](#bib.bib20), [26](#bib.bib26)]: ###### **Proposition 2.2****.** *i) For any initial data $Y_{0}\in{\cal E}$ there exists the unique solution $Y(t)\in C({\mathbb{R}},\cal E)$ to the problem ([2.1](#S2.E1)). (ii) For every $t\in{\mathbb{R}}$, the map $U(t):Y_{0}\mapsto Y(t)$ is continuous in ${\cal E}$. (iii) The energy is conserved, i.e.* [TABLE] ### 2.2 Solitary manifold and main result Let us consider the solitons ([1.17](#S1.E17)). The substitution to ([1.15](#S1.E15)) gives the following stationary equations, [TABLE] ###### **Definition 2.3****.** *A soliton state is $S(\sigma):=(\psi_{v}(x-b),\pi_{v}(x-b))$, where $\sigma:=(b,v)$ with $b\in{\mathbb{R}}$ and $v\in(-1,1)$.* Obviously, the soliton solution ([1.16](#S1.E16)) admits the representation $S(\sigma(t))$, where [TABLE] ###### **Definition 2.4****.** *A solitary manifold is the set ${\cal S}:=\{S(\sigma):\sigma\in\Sigma:={\mathbb{R}}\times(-1,1)$.* The main result of our paper is the following theorem ###### **Theorem 2.5****.** *Let the conditions ***U1*** and ***U2*** hold, and $Y(t)$ be the solution to the Cauchy problem ([2.1](#S2.E1)) with an initial state $Y_{0}\in{\cal E}$ which is close to a kink $S(\sigma_{0})=S_{q_{0},v_{0}}$:* [TABLE] *where $\beta>5/2$. Then for $d_{0}$ sufficiently small the solution admits the asymptotics:* [TABLE] *where $v_{\pm}$ and $q_{\pm}$ are constants, $\Phi_{\pm}\in E$, and $W_{0}(t)$ is the dynamical group of the free Klein-Gordon equation , while* [TABLE] It suffices to prove the asymptotics ([2.9](#S2.E9)) for $t\to+\infty$ since the system ([1.15](#S1.E15)) is time reversible. ## 3 Symplectic projection ### 3.1 Symplectic structure and hamiltonian form The system ([2.1](#S2.E1)) reads as the Hamilton system [TABLE] where ${\cal D}{\cal H}$ is the Fréchet derivative of the Hamilton functional ([1.1](#S1.E1)). Let us identify the tangent space of $E$, at every point, with the space $E$. Consider the symplectic form $\Omega$ on $E$ defined by [TABLE] where [TABLE] and $\langle\psi_{1},\psi_{2}\rangle=\displaystyle\int\psi_{1}(x)\psi_{2}(x)dx$ etc. It is clear that the form $\Omega$ is non-degenerate, i.e. [TABLE] ###### **Definition 3.1****.** *i) The symbol $Y_{1}\nmid Y_{2}$ means that $Y_{1}\in E$, $Y_{2}\in E$, and $Y_{1}$ is symplectic orthogonal to $Y_{2}$, i.e. $\Omega(Y_{1},Y_{2})=0$.* *ii) A projection operator ${\bf P}:E\to E$ is said to be symplectic orthogonal if $Y_{1}\nmid Y_{2}$ for $Y_{1}\in\mbox{\rm Ker}\,{\bf P}$ and $Y_{2}\in\mbox{\rm Range}\,{\bf P}$.* ### 3.2 Symplectic projection onto solitary manifold Let us consider the tangent space ${\cal T}_{S(\sigma)}{\cal S}$ of the manifold ${\cal S}$ at a point $S(\sigma)$. The vectors [TABLE] form a basis in ${\cal T}_{S(\sigma)}{\cal S}$. Here $y:=x-b$ is the “moving frame coordinate”. Let us stress that the functions $\tau_{j}$ are always regarded as functions of $y$ rather than those of $x$. Formula ([1.17](#S1.E17)) implies that [TABLE] ###### **Lemma 3.2****.** *The symplectic form $\Omega$ is nondegenerate on the tangent space ${\cal T}_{S(\sigma)}{\cal S}$ i.e. ${\cal T}_{S(\sigma)}{\cal S}$ is a symplectic subspace.* ###### Proof. Let us compute the vectors $\tau_{1}$ and $\tau_{2}$. Recall that $\psi_{v}(y)=s(\gamma y)$ and $\pi_{v}=-v\psi^{\prime}_{v}(y)=-v\gamma s^{\prime}(\gamma y)$ with $\gamma=1/\sqrt{1-v^{2}}$. Then [TABLE] [TABLE] Therefore [TABLE] ∎ Now we show that in a small neighborhood of the soliton manifold ${\cal S}$ a “symplectic orthogonal projection” onto ${\cal S}$ is well-defined. Let us introduce the translations $T_{q}:(\psi(x),\pi(x))\mapsto(\psi(x-q),\pi(x-q))$, $q\in{\mathbb{R}}$. Note that the manifold ${\cal S}$ is invariant with respect to the translations. ###### **Definition 3.3****.** *For any $\overline{v}<1$ denote by $\Sigma(\overline{v})=\{\sigma=(b,v):b\in{\mathbb{R}},|v|\leq\overline{v}\}$.* Let us note that ${\cal S}\in E_{\alpha}$ with $\alpha<-1/2$. ###### **Lemma 3.4****.** *Let $\alpha<-1/2$ and $\overline{v}<1$. Then i) there exists a neighborhood ${\cal O}_{\alpha}({\cal S})$ of ${\cal S}$ in $E_{\alpha}$ and a mapping ${\bf\Pi}:{\cal O}_{\alpha}({\cal S})\to{\cal S}$ such that ${\bf\Pi}$ is uniformly continuous on ${\cal O}_{\alpha}({\cal S})$ in the metric of $E_{\alpha}$,* [TABLE] *ii) ${\cal O}_{\alpha}({\cal S})$ is invariant with respect to the translations $T_{q}$, and* [TABLE] *iii) For any $\overline{v}<1$ there exists an $r_{\alpha}(\overline{v})>0$ s.t. $S(\sigma)+X\in{\cal O}_{\alpha}({\cal S})$ if $\sigma\in\Sigma(\overline{v})$ and $\|X\|_{E_{\alpha}}5/2$. Then for any $X\in E_{\beta}$, the weighted energy decay holds* [TABLE] ###### **Corollary 4.9****.** *For $\sigma>5/2$ and for $X\in E_{\beta}\cap W$* [TABLE] *Here $(\cdot)_{1}$ stands for the first component of the vector function.* ###### Proof. Let us apply the projector ${\bf P}^{c}_{v}$ to both sides of ([4.29](#S4.E29)): [TABLE] where [TABLE] Hence, the Duhamel representation gives, [TABLE] Let us note, that $e^{A_{v}^{0}t}Z=e^{A_{0}^{0}t}T_{v}Z$, where $T_{v}Z(x,t)=Z(x+vt,t)$. Then ([4.34](#S4.E34)) reads [TABLE] Applying estimate (265) from [[21](#bib.bib21)], the Hölder inequality and Proposition [4.8](#S4.Thmtheorem8) we obtain [TABLE] ∎ ### 4.3 Taylor expansion for nonlinear term Now let us expand $N(v,\Psi)$ from ([4.9](#S4.E9)) in the Taylor series [TABLE] where [TABLE] and $N_{R}$ is the remainder. By condition **U1** we have [TABLE] Hence, the functions $F^{(j)}(\psi_{v}(y))$, $2\leq j\leq 12$ decrease exponentially as $|y|\to\infty$ by ([1.21](#S1.E21)) and ([1.17](#S1.E17)). Therefore, [TABLE] For the remainder $N_{R}$ we have [TABLE] where ${\cal R}(A)$ is a general notation for a positive function which remains bounded as $A$ is sufficiently small. ###### **Lemma 4.10****.** *The bounds hold* [TABLE] [TABLE] ###### Proof. *Step i)* By the Cauchy formula, [TABLE] Therefore, [TABLE] since $\|\Psi\|_{L^{2}}\leq C(d_{0})$ by the results of [[10](#bib.bib10)]. Differentiating ([4.42](#S4.E42)), we obtain [TABLE] Hence, [TABLE] since $\displaystyle\int|\Psi(x)\Psi^{\prime}(x)|dx\leq\|\Psi\|_{L^{2}}\|\Psi^{\prime}\|_{L^{2}}\leq C(d_{0})$. Finally, note that [TABLE] Then ([4.40](#S4.E40)) follows. *Step ii)* The bound ([4.39](#S4.E39)) implies [TABLE] We will prove in Appendix B that [TABLE] Then ([4.41](#S4.E41)) follows. ∎ ###### **Remark 4.11****.** Our choice of the degree $14$ in the condition ([1.11](#S1.E11)) is due to the competition between the factors in the estimate ([4.41](#S4.E41)) for the remainder. Namely, the factor $(1+t)^{4+\nu}$ with $\nu<1/2$ comes from the virial type estimate ([4.43](#S4.E43)) describing the expansion of the support for the perturbation of the kink. On the other hand, $\|\Psi\|^{12}_{L^{\infty}}\sim t^{-6}$ by the crucial decay estimate ([7.1](#S7.E1)). Hence, the right hand side ([4.41](#S4.E41)) decays like $\sim t^{-2+\nu}$ where $-2+\nu<-3/2$ which is sufficient for the method of majorants (in integral inequalities ([9.2](#S9.E2)) and ([9.3](#S9.E3))). ## 5 Symplectic decomposition of the dynamics Here we decompose the dynamics in two components: along the manifold ${\cal S}$ and in transversal directions. The equation ([4.5](#S4.E5)) is obtained without any assumption on $\sigma(t)$ in ([4.1](#S4.E1)). We are going to choose $S(\sigma(t)):={\bf\Pi}Y(t)$, but then we need to know that [TABLE] with some ${\cal O}_{\alpha}({\cal S})$ defined in Lemma [3.4](#S3.Thmtheorem4). It is true for $t=0$ by our main assumption ([2.8](#S2.E8)) with sufficiently small $d_{0}>0$. Then $S(\sigma(0))={\bf\Pi}Y(0)$ and $X(0)=Y(0)-S(\sigma(0))$ are well defined. We will prove below that ([5.1](#S5.E1)) holds with $\alpha=-\beta$ if $d_{0}$ is sufficiently small. First, we choose $\overline{v}<1$ such that [TABLE] Denote by $r_{-\beta}(\overline{v})$ the positive number from Lemma [3.4](#S3.Thmtheorem4) iii) which corresponds to $\alpha=-\beta$. Then $S(\sigma)+X\in{\cal O}_{-\beta}({\cal S})$ if $\sigma=(b,v)$ with $|v|<\overline{v}$ and $\|X\|_{E_{-\beta}}0$, we have for $t_{1}0.$ In order to prove the bound ([10.5](#S10.E5)), we will need the formula ([4.30](#S4.E30)) and the following relation which follows from ([4.30](#S4.E30)): [TABLE] where $\tilde{\tau}_{j}(v(t))$ are the vectors $\tau_{j}(v(t))$ expressed in the variables $y_{1}$. In detail (cf. ([3.3](#S3.E3))), [TABLE] where $v=v(t)$. Since $\tau^{\prime}_{j}$ are smooth and rapidly decaying at infinity functions, then Lemma [9.1](#S9.Thmtheorem1) implies [TABLE] Furthermore, [TABLE] and therefore [TABLE] [TABLE] since $|\partial_{v}p_{jl}(v(s))|$ is uniformly bounded by ([5.2](#S5.E2)). Further, [TABLE] Hence, the bounds ([10.5](#S10.E5)) will follow from ([4.30](#S4.E30)), ([10.6](#S10.E6)) and ([10.8](#S10.E8))-([10.10](#S10.E10)) if we choose $\varepsilon>0$ small enough. The proof is completed. ∎ ## 11 Decay of transversal component Here we prove Proposition [7.1](#S7.Thmtheorem1). *Step i)* We fix $\varepsilon>0$ and $t_{*}=t_{*}(\varepsilon)$ for which Lemma [10.1](#S10.Thmtheorem1) holds. Then the bounds of type ([9.7](#S9.E7)) and ([9.8](#S9.E8)) holds with $\|{\bf P}_{1}^{d}\tilde{X}(t)\|_{E_{-\beta}}$ and $\|({\bf P}_{1}^{d}\tilde{X}(t))_{1}\|_{L^{\infty}}$ in the left hand sides replaced by $\|X(t)\|_{E_{-\beta}}$ and $\|\Psi(t)\|_{L^{\infty}}$ : [TABLE] [TABLE] for $0\leq t\leq t_{1}$ and $t_{1}0$ if $d_{0}$ is small enough. ## 12 Soliton asymptotics Here we prove our main Theorem [2.5](#S2.Thmtheorem5) under the assumption that the decay ([7.1](#S7.E1)) holds. The estimates ([6.11](#S6.E11)) and ([7.1](#S7.E1)) imply that [TABLE] Therefore, $c(t)=c_{+}+{\cal O}(t^{-2})$ and $v(t)=v_{+}+{\cal O}(t^{-2})$, $t\to\infty$. Similarly, [TABLE] We have obtained the solution $Y(x,t)=(\psi(x,t),\pi(x,t))$ to ([1.15](#S1.E15)) in the form [TABLE] where we define now $v(t)=\dot{b}(t)=v_{+}+\dot{\alpha}(t)$. Since [TABLE] it remains to extract the dispersive wave $W_{0}(t)\Phi_{+}$ from the term $X(x-b(t),t)$. Substituting ([12.3](#S12.E3)) into ([1.15](#S1.E15)) we obtain by ([2.6](#S2.E6)) the inhomogeneous Klein-Gordon equation for the $X(x-b(t),t)$: [TABLE] where $y=x-b(t)$, and [TABLE] Now we change the variable $y\mapsto y_{1}=y+\alpha(t)+q_{+}$. Then we obtain the “frozen” equation [TABLE] where $\tilde{X}(t)$ and $\tilde{R}(t)$ are $X(t)$ and $R(t)$ of $y=y_{1}-\alpha(t)-q_{+}$, and [TABLE] Equation ([12.5](#S12.E5)) implies [TABLE] where $W_{+}(t)=e^{A_{+}t}$ is the integral operator with integral kernel [TABLE] since by ([12.2](#S12.E2)) [TABLE] Hence, equation ([12.7](#S12.E7)) implies [TABLE] Let us rewrite ([12.8](#S12.E8)) as [TABLE] To establish the asymptotics ([2.9](#S2.E9)), it suffices to prove that [TABLE] Assumption ([2.8](#S2.E8)) implies that $\tilde{X}(0)\in E$. Let us split $\tilde{R}(s)$ as the sum [TABLE] By ([12.1](#S12.E1)), we obtain [TABLE] Let us consider $\tilde{R}^{\prime\prime}=(0,\tilde{R}^{\prime\prime}_{2})$. We have [TABLE] By ([1.20](#S1.E20)) and ([7.1](#S7.E1)), we obtain [TABLE] since $|q_{+}+\alpha(t)|\leq C$. Finally, ([7.1](#S7.E1)), ([7.4](#S7.E4)), and ([8.6](#S8.E6)) imply [TABLE] Hence, ([12.11](#S12.E11))-([12.12](#S12.E12)) imply [TABLE] and ([12.9](#S12.E9)) follows by ([12.10](#S12.E10)) and ([12.13](#S12.E13)). ## A Virial type estimates Here we prove the weighted estimate ([4.43](#S4.E43)). Let us recall that we split the solution $Y(t)=(\psi(\cdot,t),\pi(\cdot,t))=S(\sigma(t))+X(t)$, and denote $X(t)=(\Psi(t),\Pi(t))$, $(\Psi_{0},\Pi_{0}):=(\Psi(0),\Pi(0))$. Our basic condition ([2.8](#S2.E8)) implies that for some $\nu>0$ [TABLE] ###### **Proposition A.1****.** *Let the potential $U$ satisfy conditions ***U1***, and $\Psi_{0}$ satisfy ([A. 1](#S1.E1a)). Then the bounds hold* [TABLE] We will deduce the proposition from the following two lemmas. The first lemma is well known. Denote [TABLE] ###### **Lemma A.2****.** *For the solution $\psi(x,t)$ of Klein-Gordon equation ([1.2](#S1.E2)) the local energy estimate holds* [TABLE] ###### Proof. The estimate follows by standard arguments: multiplication of the equation ([1.2](#S1.E2)) by $\dot{\psi}(x,t)$ and integration over the trapezium $ABCD$, where $A=(a_{1}-t,0)$, $B=(a_{1},t)$, $C=(a_{2},t)$, $D=(a_{2}+t,0)$. Then ([A. 3](#S1.E3a)) is obtained after partial integration using that $U(\psi)\geq 0$. ∎ ###### **Lemma A.3****.** *For any $\sigma\geq 0$* [TABLE] ###### Proof. By ([A. 3](#S1.E3a)) [TABLE] Hence, [TABLE] Obviously, [TABLE] with some $c(\sigma)>0$. On the other hand, [TABLE] since $\sigma\geq 0$. Finally, ([A. 5](#S1.E5a))-([A. 7](#S1.E7a)) imply ([A. 4](#S1.E4a)). ∎ **Proof of Proposition [A.1](#S1.Thmtheorem1a)** First, we verify that [TABLE] Indeed, $\psi_{0}(x)=\psi_{v_{0}}(x-q_{0})+\Psi_{0}(x)$ is bounded since $\Psi_{0}\in H^{1}({\mathbb{R}})$. Hence **U1** implies that [TABLE] and then ([A. 8](#S1.E8a)) follows by ([1.17](#S1.E17)), ([1.21](#S1.E21)) and ([A. 1](#S1.E1a)). Further, we have [TABLE] Using the bounds ([6.11](#S6.E11)) we obtain [TABLE] [TABLE] Hence ([4.2](#S4.E2)) implies that [TABLE] Substituting ([A](#S1.Ex59)) into ([A. 9](#S1.E9a)) and changing variables we obtain by ([A. 4](#S1.E4a)), ([A. 10](#S1.E10a)) and ([A. 8](#S1.E8a)) [TABLE] ## References - [1] S. Agmon, Spectral properties of Schrödinger operator and scattering theory, *Ann. Scuola Norm. Sup. Pisa*, Ser. IV **2**, 151-218 (1975). - [2] F. A. Bais, Topological excitations in gauge theories; An introduction from the physical point of view, Springer Lecture Notes in Mathematics, vol. 926 (1982). - [3] F. Bj$\o$rn, Geometry, Particles, and Fields, Springer, NY, 1998. - [4] V.S. Buslaev, G.S. Perelman, Scattering for the nonlinear Schrödinger equations: states close to a soliton, *St. Petersburg Math. J.* **4** (1993), no.6, 1111-1142. - [5] V.S. Buslaev, C. Sulem, On asymptotic stability of solitary waves for nonlinear Schrödinger equations,*Ann. Inst. Henri Poincaré, Anal. Non Linéaire* **20**(2003), no.3, 419-475. - [6] S. Cuccagna, Stabilization of solutions to nonlinear Schrödinger equations, *Comm. Pure Appl. Math.* **54** (2001), 1110-1145. - [7] S. Cuccagna, On asymptotic stability in 3D of kinks for the $\phi^{4}$ model, *Transactions of AMS* **360** (2008), no. 5, 2581-2614. - [8] A. Jensen, T. Kato, Spectral properties of Schrödinger operators and time-decay of the wave functions, *Duke Math. J.* **46**, 583-611 (1979). - [9] A. Jensen, G. Nenciu, A unified approach to resolvent expansions at thresholds, *Reviews in Math. Phys.* **13** (2001), no. 6, 717-754. - [10] D.B. Henry, J.F. Perez, W.F. Wreszinski, Stability theory for solitary-wave solutions of scalar field equations, *Comm. Math. Phys.* **85** (1982), 351-361. - [11] V. Imaikin, A.I. Komech, B. Vainberg, On scattering of solitons for the Klein-Gordon equation coupled to a particle, *Comm. Math. Phys.* **268** (2006), no.2, 321-367. - [12] E. Kirr, A. Zarnesku, On the asymptotic stability of bound states in 2D cubic Schro”dinger equation, *Comm. Math. Phys.* **272** (2007), no. 2, 443-468. - [13] Komech A., Kopylova E., Weighted energy decay for 1D Klein-Gordon equation, *Comm. PDE* **35**, no.2, 353-374 (2010). - [14] E. Kopylova, On long-time decay for Klein-Gordon equation, ArXiv: 1009.2649 - [15] J.L. Lions, “Quelques Mèthodes de Rèsolution des Problémes aux Limites non Linéaires”, Paris, Dunod, 1969. - [16] M. Murata, Asymptotic expansions in time for solutions of Schrödinger-type equations, *J. Funct. Anal.***49**, 10-56 (1982). - [17] J. Miller, M. Weinstein, Asymptotic stability of solitary waves for the regularized long-wave equation, *Comm. Pure Appl. Math.* **49** (1996), no. 4, 399-441. - [18] R.L. Pego, M.I. Weinstein, Asymptotic stability of solitary waves, *Commun. Math. Phys.* **164** (1994), 305-349. - [19] C.A. Pillet, C.E. Wayne, Invariant manifolds for a class of dispersive, Hamiltonian, partial differential equations, *J. Differ. Equations* **141** (1997), No.2, 310-326. - [20] M. Reed, “Abstract Non-Linear Wave Equations”, Lecture Notes in Mathematics 507 (1976), Springer, Berlin. - [21] M. Reed, B. Simon, Methods of Modern Mathematical Physics, III, Academic Press, 1979. - [22] I. Rodnianski, W. Schlag, A. Soffer, Dispersive analysis of charge transfer models, *Commun. Pure Appl. Math.* **58** (2005), no. 2, 149-216. - [23] A. Soffer, M.I. Weinstein, Multichannel nonlinear scattering for nonintegrable equations, *Comm. Math. Phys.* **133** (1990), 119-146. - [24] A. Soffer, M.I. Weinstein, Multichannel nonlinear scattering for nonintegrable equations. II. The case of anisotropic potentials and data, *J. Differential Equations* **98** (1992), no. 2, 376-390. - [25] A. Soffer, M.I. Weinstein, Resonances, radiation damping and instability in Hamiltonian nonlinear wave equations, *Invent. Math.* **136** (1999), 9-74. - [26] W.A. Strauss, Nonlinear invariant wave equations, Lecture Notes in Physics 73 (1978), Springer, Berlin, 197-249. - [27] T.-P. Tsai, H.-T. Yau, Asymptotic dynamics of nonlinear Schrödinger equations: resonance-dominated and dispersion-dominated solutions, *Commun. Pure Appl. Math.* **55** (2002), no.2, 153-216.