\documentclass[CRMATH,Unicode,biblatex,XML]{cedram}

\TopicFR{Équations aux dérivées partielles, Physique mathématique}
\TopicEN{Partial differential equations, Mathematical physics}

\addbibresource{CRMATH_Carles_20240143.bib}


\usepackage{mathtools}
\usepackage{amssymb}

\newcommand\R{{\mathbb R}}% real numbers
\renewcommand\H{{\mathcal H}}
\newcommand\Sch{{\mathcal S}}% Schwartz space
\renewcommand\O{\mathcal O}
\newcommand\F{\mathcal F}
\renewcommand\le{\leqslant}
\renewcommand\ge{\geqslant}


\newcommand\Eq[2]{\mathop{\sim}\limits_{#1\to#2}}
\newcommand\Tend[2]{\mathop{\to}\limits_{#1\to#2}}

\renewcommand\d{{\partial}}
\newcommand\eps{\varepsilon}
\newcommand\fum{{\hat{u}_-}}
\newcommand\um{{u_-}}
\renewcommand\u{{\mathbf u}}
\newcommand\W{{W_-^{\mathrm{mod}}}}
\newcommand\Scatt{{S^{\mathrm{mod}}}}
\newcommand\cR{{\mathcal R}}

\DeclareMathOperator{\RE}{Re}
\DeclareMathOperator{\IM}{Im}
\DeclareMathOperator{\diver}{div}
\DeclareMathOperator{\expe}{e}



\numberwithin{equation}{section}


\newcommand{\dd}{\mathrm{d}}
\newcommand{\irm}{\mathrm{i}}
\newcommand{\app}{\mathrm{app}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\makeatletter
\let\save@mathaccent\mathaccent
\newcommand*\if@single[3]{%
 \setbox0\hbox{${\mathaccent"0362{#1}}^H$}%
 \setbox2\hbox{${\mathaccent"0362{\kern0pt#1}}^H$}%
 \ifdim\ht0=\ht2 #3\else #2\fi
 }
%The bar will be moved to the right by a half of \macc@kerna, which is computed by amsmath:
\newcommand*\rel@kern[1]{\kern#1\dimexpr\macc@kerna}
%If there's a superscript following the bar, then no negative kern may follow the bar;
%an additional {} makes sure that the superscript is high enough in this case:
\newcommand*\widebar[1]{\@ifnextchar^{{\wide@bar{#1}{0}}}{\wide@bar{#1}{1}}}
%Use a separate algorithm for single symbols:
\newcommand*\wide@bar[2]{\if@single{#1}{\wide@bar@{#1}{#2}{1}}{\wide@bar@{#1}{#2}{2}}}
\newcommand*\wide@bar@[3]{%
 \begingroup
 \def\mathaccent##1##2{%
%Enable nesting of accents:
 \let\mathaccent\save@mathaccent
%If there's more than a single symbol, use the first character instead (see below):
 \if#32 \let\macc@nucleus\first@char \fi
%Determine the italic correction:
 \setbox\z@\hbox{$\macc@style{\macc@nucleus}_{}$}%
 \setbox\tw@\hbox{$\macc@style{\macc@nucleus}{}_{}$}%
 \dimen@\wd\tw@
 \advance\dimen@-\wd\z@
%Now \dimen@ is the italic correction of the symbol.
 \divide\dimen@ 3
 \@tempdima\wd\tw@
 \advance\@tempdima-\scriptspace
%Now \@tempdima is the width of the symbol.
 \divide\@tempdima 10
 \advance\dimen@-\@tempdima
%Now \dimen@ = (italic correction / 3) - (Breite / 10)
 \ifdim\dimen@>\z@ \dimen@0pt\fi
%The bar will be shortened in the case \dimen@<0 !
 \rel@kern{0.6}\kern-\dimen@
 \if#31
 \overline{\rel@kern{-0.6}\kern\dimen@\macc@nucleus\rel@kern{0.4}\kern\dimen@}%
 \advance\dimen@0.4\dimexpr\macc@kerna
%Place the combined final kern (-\dimen@) if it is >0 or if a superscript follows:
 \let\final@kern#2%
 \ifdim\dimen@<\z@ \let\final@kern1\fi
 \if\final@kern1 \kern-\dimen@\fi
 \else
 \overline{\rel@kern{-0.6}\kern\dimen@#1}%
 \fi
 }%
 \macc@depth\@ne
 \let\math@bgroup\@empty \let\math@egroup\macc@set@skewchar
 \mathsurround\z@ \frozen@everymath{\mathgroup\macc@group\relax}%
 \macc@set@skewchar\relax
 \let\mathaccentV\macc@nested@a
%The following initialises \macc@kerna and calls \mathaccent:
 \if#31
 \macc@nested@a\relax111{#1}%
 \else
%If the argument consists of more than one symbol, and if the first token is
%a letter, use that letter for the computations:
 \def\gobble@till@marker##1\endmarker{}%
 \futurelet\first@char\gobble@till@marker#1\endmarker
 \ifcat\noexpand\first@char A\else
 \def\first@char{}%
 \fi
 \macc@nested@a\relax111{\first@char}%
 \fi
 \endgroup
}
\makeatother

\let\oldbar\bar
\renewcommand*{\bar}[1]{\mathchoice{\widebar{#1}}{\widebar{#1}}{\widebar{#1}}{\oldbar{#1}}}

\let\oldtilde\tilde
\renewcommand*{\tilde}[1]{\mathchoice{\widetilde{#1}}{\widetilde{#1}}{\oldtilde{#1}}{\oldtilde{#1}}}

\let\oldhat\hat
\renewcommand*{\hat}[1]{\mathchoice{\widehat{#1}}{\widehat{#1}}{\oldhat{#1}}{\oldhat{#1}}}

\renewcommand*{\to}{\mathchoice{\longrightarrow}{\rightarrow}{\rightarrow}{\rightarrow}}
\let\oldmapsto\mapsto
\renewcommand*{\mapsto}{\mathchoice{\longmapsto}{\oldmapsto}{\oldmapsto}{\oldmapsto}}


\let\oldcheck\check
\renewcommand*{\check}[1]{\mathchoice{\widecheck{#1}}{\widecheck{#1}}{\oldcheck{#1}}{\oldcheck{#1}}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newcommand*{\mk}{\mkern -1mu}
\newcommand*{\Mk}{\mkern -2mu}
\newcommand*{\mK}{\mkern 1mu}
\newcommand*{\MK}{\mkern 2mu}

\hypersetup{urlcolor=purple, linkcolor=blue, citecolor=red}


\newcommand*{\romanenumi}{\renewcommand*{\theenumi}{\roman{enumi}}}
\newcommand*{\Romanenumi}{\renewcommand*{\theenumi}{\Roman{enumi}}}
\newcommand*{\alphenumi}{\renewcommand*{\theenumi}{\alph{enumi}}}
\newcommand*{\Alphenumi}{\renewcommand*{\theenumi}{\Alph{enumi}}}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\title[Long range scattering for cubic NLS]{Dynamics near the origin
 of the long range scattering for the one-dimensional Schr\"odinger equation} 
 
\alttitle{Dynamique près de l'origine du scattering longue portée pour l'équation de Schrödinger unidimensionnelle} 
 
\author{\firstname{R\'emi} \lastname{Carles}}

\address{Univ Rennes, CNRS, IRMAR -- UMR 6625, F-35000
 Rennes, France}

\email{Remi.Carles@math.cnrs.fr}

\keywords{\kwd{Nonlinear Schrödinger equation} \kwd{long range scattering} \kwd{asymptotic expansion} \kwd{error estimate}}

\altkeywords{\kwd{Équation de Schrödinger non linéaire} \kwd{scattering longue portée} \kwd{développement asymptotique} \kwd{estimation d'erreur}}

\subjclass{35Q55, 35P25, 35B40}


\begin{abstract}
 We consider the cubic Schr\"odinger equation on the line, for which
 the scattering theory requires modifications due to long range
 effects. We revisit the construction of the modified wave operator,
 and recall the construction of its inverse, in order to describe the
 asymptotic behavior of these operators near the origin. At leading
 order, these operators, whose definition includes a nonlinear
 modification in the phase compared to the linear dynamics, correspond to the identity. We compute
 explicitly the first corrector in the asymptotic expansion, and
 justify this expansion by error estimates. 
\end{abstract}


\begin{altabstract}
 Nous considèrons l'équation de Schrödinger cubique sur la droite, pour laquelle la théorie du scattering demande des modifications dues aux effets à longue portée. Nous reprenons la construction de l'opérateur d'onde modifié, et rappelons la construction de son inverse, afin de décrire le comportement de ces opérateurs près de l'origine. Au premier ordre, ces opérateurs, dont la définition contient une modification non linéaire de la phase par rapport à la dynamique linéaire, coïncident avec l'identité. Nous calculons explicitement le premier correcteur du développement asymptotique, et justifions ce développement par des estimations d'erreur. 
\end{altabstract}

\thanks{This work was supported by Centre Henri Lebesgue,
program ANR-11-LABX-0020-0.}


\begin{document}
\maketitle




\section{Introduction}
\label{sec:intro}
We consider the cubic Schr\"odinger equation on the line,
\begin{equation}
 \label{eq:NLS}
 \irm\d_t u +\frac{1}{2}\d_{x}^2u =\lambda |u|^2u,\quad x\in \R, 
\end{equation}
with $\lambda \in \R$. It is well known that in terms of scattering,
this equation corresponds to a borderline case, where long range
effects appear (see~\cite{Barab}). The goal of this paper is to
analyze the modified scattering map near the origin, as well as the
modified wave operator and its inverse. We first recall a few aspects
of scattering theory for the nonlinear Schr\"odinger equation in the short
range case, then turn to the specificities of this long range setting. 

\subsection{Short range scattering}




Denote by $U(t)$ the Schr\"odinger group,
\begin{equation*}
 U(t)\coloneqq \expe^{\irm\frac{t}{2}\d_x^2}.
\end{equation*}
In view of the explicit formula for $U(t)$ as a convolution operator,
we have
\begin{equation}\label{eq:asym-lin}
 U(t)f(x)\Eq t {\pm\infty} \expe^{\irm\frac{x^2}{2t}} \frac{1}{(\irm t)^{1/2}}\hat
 f\left(\frac{x}{t}\right),
\end{equation}
where we normalize the Fourier transform as follows,
\begin{equation*}
 {\F} f(\xi)=\hat
 f(\xi)=\frac{1}{\sqrt{2\pi}}\int_{\R}f(x)\expe^{-\irm x \xi} \, \dd x.
\end{equation*}
See Lemma~\ref{lem:R} for a more precise statement regarding
\eqref{eq:asym-lin}. 
Consider, for the
sake of comparison with~\eqref{eq:NLS}, the case of a quintic,
defocusing nonlinearity
\begin{equation}
 \label{eq:NLSquintic}
 \irm\d_t u +\frac{1}{2}\d_{x}^2u = |u|^4 u,\quad x\in \R.
\end{equation}
The discussion for, e.g., the two-dimensional or the three-dimensional
cubic Schr\"odinger equation would be similar. Given 
\begin{equation}\label{eq:Sigma}
 u_-\in \Sigma\coloneqq \left\{f\in H^1(\R);\
 \|f\|_\Sigma\coloneqq \|f\|_{L^2}+\|\d_x f\|_{L^2} + \|x
 f\|_{L^2}<\infty\right\},
\end{equation}
there exists a unique $u_0\in \Sigma=H^1\cap \F(H^1)$ such that the
(unique, global) 
solution $u$ to~\eqref{eq:NLSquintic} with $u_{\mid t=0}=u_0$
satisfies 
\begin{equation*}
 \|U(-t)u(t)-u_-\|_{\Sigma}\Tend t {-\infty}0.
\end{equation*}
We recall that $U(t)$ is unitary on $H^1$, but not on $\F(H^1)$, this
is why the quantity measured above is not $u(t)-U(t)u_-$. 
The map $u_-\mapsto u_0$ is classically referred to as \emph{wave
 operator} (see e.g.~\cite{Ginibre}).

Conversely, given $u_0\in \Sigma$, there exists $u_+\in \Sigma$ such
that 
\begin{equation*}
 \|U(-t)u(t)-u_+\|_{\Sigma}\Tend t {+\infty}0.
\end{equation*}
The map $u_-\mapsto u_+$ is called the \emph{scattering operator}. It
can be defined for other (defocusing) nonlinearities, but strictly
supercubic, in view of~\cite{Barab}, where it is proved that,
typically for~\eqref{eq:NLS}, it is not possible to compare the
nonlinear dynamics with the linear one for large time (see also~\cite{Ginibre}). We will see in the next
subsection that in the cubic case~\eqref{eq:NLS}, nontrivial long range
effects must be taken into account, and that these effects are explicit.
\smallbreak

In general, rather little is known regarding properties of the
scattering map $S$. For instance, it is proven in~\cite{CaGa09} that
for smooth, power-like nonlinearities such that short range scattering
is known, the wave and scattering operators are analytic. The formula
for the associated Taylor series is given, and the formula differs
whether the expansion is considered at the origin or at a nontrivial
state. Such asymptotic expansions (not necessarily in the analytic
case) have proven useful in the context of inverse problems, see
e.g.~\cite{MR4576319,KiMuVi-p} and references therein. 


\subsection{Long range case: modified scattering}
\label{sec:mod}

The picture is different in the case of~\eqref{eq:NLS}: the nonlinear
dynamics cannot be compared to the linear one, unless one considers
the trivial case $u\equiv 0$ (\cite{Barab}), so a modified scattering
theory has been developed in order to describe the large time behavior
of $u$, sharing some similarities with the linear case, presented in
e.g.~\cite{DG}. 
\smallbreak



Given asymptotic states $u_-$ and $u_+$, we introduce the long range phase
corrections (different whether $t\to -\infty$ of $t\to +\infty$),
\begin{equation}
 \label{eq:defS}
 S_\pm(t,x)\coloneqq \mp\lambda \left| \hat{u}_\pm\left( \frac{x}{t}
 \right)\right|^2 \log |t|.
\end{equation}
Loosely speaking, the existence of modified wave operator reads as
follows: given $u_-$ (sufficiently small), there exists a 
solution $u$ to~\eqref{eq:NLS} such that 
\begin{equation*}
 u(t,x)\Eq t {-\infty} \expe^{\irm S_-(t,x)}U(t)u_-(x)\Eq t {-\infty}
 \expe^{\irm S_-(t,x)+\irm \frac{x^2}{2t}} \frac{1}{(\irm t)^{1/2}}\hat 
 u_-\left(\frac{x}{t}\right),
\end{equation*}
where the second approximation stems from~\eqref{eq:asym-lin}. Note
that the phase modification $S_-$ is by no means negligible 
as $t\to -\infty$: it accounts for long range effects, as established
initially in~\cite{Ozawa91}. We denote by $u_{\mid t=0} =\W (u_-)$ the
modified wave operator. 
\smallbreak

The modified asymptotic completeness is similar: given $u_0$
(sufficiently small), there exists $u_+$ such that the solution $u$ to
\eqref{eq:NLS} with $u_{\mid t=0}=u_0$ satisfies
\begin{equation*}
 u(t)\Eq t {+\infty} \expe^{\irm S_+(t)}U(t)u_+.
\end{equation*}
Such a result was proven initially in~\cite{HN98}. 
The modified scattering map is given by $\Scatt (u_-)=u_+$. At this stage,
we have not addressed the function spaces in which the above
asymptotics have been proven. 
\smallbreak

In~\cite{Ozawa91}, the existence of modified wave operators was
established for $u_-\in \F(H^2)$, with $\|\hat u_-\|_{L^\infty}$
sufficiently small, and the solution $u$ to~\eqref{eq:NLS} has an
$L^2$ regularity, $u_0\in L^2(\R)$. We emphasize, as the notation will
be used many times, that the space $\F(H^s)$ is characterized by
\begin{equation*}
 \F(H^s)=\left\{f\in \Sch'(R),\ \|f\|_{\F(H^s)}^2\coloneqq 
 \int_{\R}\left\langle x\right\rangle ^{2s}|f(x)|^2\, \dd x<\infty\right\}, \quad \left\langle x\right\rangle =\sqrt{1+x^2}.
\end{equation*}
The result of~\cite{Ozawa91} also addresses the case where $u_-\in \H$
(defined below, see~\eqref{eq:H}) where, provided again that
$\|\hat u_-\|_{L^\infty}$ 
sufficiently small, the solution $u$ to~\eqref{eq:NLS} has an
$H^1$ regularity, and convergence holds in this space.

In~\cite{HN98}, the asymptotic completeness was proven for $u_0\in
H^\gamma\cap \F(H^\gamma)$ with $\gamma>1/2$ and
$\|u_0\|_{H^\gamma\cap \F(H^\gamma)}$ sufficiently small. The obtained
asymptotic state $u_+$ is such that $\hat u_+\in L^2\cap L^\infty$. 

Denote
\begin{equation}\label{eq:H}
 \begin{aligned}
{\H}{}\coloneqq{} & \{ f\in H^1(\R);\ \left\langle x\right\rangle \d_x f, \left\langle x\right\rangle ^3f \in L^2(\R)\}\\
{}={}& \{f\in {\mathcal S}'(\R); \|f\|_{\H}\coloneqq 
\|\left\langle x\right\rangle \d_xf \|_{L^2}+
\|\left\langle x\right\rangle ^3f \|_{L^2}
<\infty \}.
\end{aligned}
\end{equation}
In~\cite{CaCMP}, the main result of~\cite{Ozawa91} was adapted for
$u_-\in \H$ with $\|\hat u_-\|_{L^\infty}$
sufficiently small, and the regularity of the solution $u$ to
\eqref{eq:NLS} was proven to be at least $\Sigma$,
hence $u_0\in \Sigma$, making it possible to connect this 
result with the asymptotic completeness from~\cite{HN98}, thus
defining a map $u_-\mapsto u_+$ from (a subset of) $\H$ to
$L^2\cap L^\infty$. When invoking the modified scattering operator, we
refer to this notion.

This gap in regularity between $u_-$ and $u_+$ was considerably
diminished in~\cite{HN06}, where, with the same notations as above,
the authors consider the setting (along with smallness
conditions)
\begin{equation*}
 u_-\in \F(H^\alpha),\quad u_0\in \F(H^\beta),\quad u_+\in \F(H^\delta),
\end{equation*}
with the constraints $1/2<\delta<\beta<\alpha<1$, allowing these three
indices to be arbitrarily close one from another. This is
achieved by adapting the notion of (modified) asymptotic completeness
in order to avoid loss of differentiability issues, in space
dimensions two and three, where the borderline nonlinearity in terms
of scattering is $|u|^{2/d}u$.
\smallbreak

We also emphasize that there are many references addressing the theory
of long range scattering for nonlinear Schr\"odinger equations
(e.g.~\cite{MR3724144,MR2047418,MR2864547,LindbladMurphy2006,MaMi18,MaMi19,MTT03}), as
well as for 
Schr\"odinger-like equations, like Hartree 
equation (e.g.~\cite{MR1855975,MR3192651,MR3359525,Wada2000}), derivative nonlinear Schr\"odinger equation
(e.g.~\cite{MR1618664,MR3144794}), Maxwell--Schr\"odinger system (e.g.~\cite{MR2474176,MR2342882}),
wave-Schr\"odinger system (e.g.~\cite{MR2853555}), not to mention other
dispersive equations.
\smallbreak

We note that the very definition of the modified wave and scattering
operators encodes the fact that the nonlinearity is cubic, and
recovering the nonlinearity from the scattering map does not make
sense, contrary to the case of~\cite{MR4576319,KiMuVi-p}. In the case
where $\lambda$ is allowed to depend on $x$ in a somehow perturbative
way, Chen and Murphy~\cite{ChenMurphy} showed that the inverse of the
modified wave operator uniquely determines $\lambda$. One of the tools
of the proof there is to study the behavior of this operator near the
origin, thanks to a rather implicit expansion
(\cite[Proposition~4.1]{ChenMurphy}). We present a 
more explicit formula in the next subsection, when
$\lambda$ is constant. 


To fix the ideas, we summarize the above discussion by introducing the
following definition, where regularity aspects are left out for
simplicity. 
\begin{definition}
 Given $u_-$, the modified wave operator acting on $u_-$ is given by
 $u_{\mid t=0} =\W (u_-)$, where $u$ solves~\eqref{eq:NLS}, and
 satisfies
 \begin{equation*}
 u(t,x)= \expe^{\irm S_-(t,x)}U(t)u_-(x) +o(1) \text{ in }L^2(\R)\text{ as }t\to-\infty,
\end{equation*}
where $S_-$ is defined in~\eqref{eq:defS}.

The modified scattering
operator is defined by $\Scatt (u_-)=u_+$ if the above solution $u$
satisfies in addition 
 \begin{equation*}
 u(t,x)= \expe^{\irm S_+(t,x)}U(t)u_+(x) +o(1) \text{ in }L^2(\R)\text{ as }t\to+\infty,
\end{equation*}
where $S_+$ is defined in~\eqref{eq:defS}.
\end{definition}


\subsection{Main result}
In the same spirit as what has been achieved for the wave and
scattering operators in the short range case, we consider the
asymptotic behavior of the modified wave and scattering operators,
with two restrictions compared to~\cite{CaGa09}. First, we shall
confine ourselves to the asymptotic expansion near the origin. Second, we
compute only the first two terms of this asymptotic expansion. We will
see that this already requires some amount of work, but that the
same method should provide some, if not all,
higher order terms in the expansion at the origin. On the other hand,
the description of these operators near a nontrivial state certainly
requires a different approach. \smallbreak

We emphasize that the question addressed here is different from the
(higher order) asymptotic expansion of the large time 
behavior of $u(t)$, as studied in~\cite{MR1913680}.
Our main results are gathered in the
following statement:

\begin{theorem}\label{theo:main}
 Let $v_-\in \H$ and $v_0\in \Sigma$.
 \begin{itemize}
 \item For any $0<\eta<2$, we have, in $\Sigma$ and as $\eps\to 0$,
 \begin{equation*}
 \W(\eps v_-) = \eps v_- + \eps^3 w_2 +\O\left(\eps^{5-\eta}\right),
 \end{equation*}
 where $w_2\in \Sigma$ is defined by
 \begin{multline*}
 w_2 = -\irm\lambda\int_{-\infty}^{-1}\left( U(-\tau)\left(
 |U(\tau)v_-|^2U(\tau)v_-\right) +\frac{1}{|\tau|} \F^{-1}\left( 
 |\hat v_-|^2\hat v_-\right)\right) \dd \tau\\
\quad -\irm \lambda\int_{-1}^0 U(-\tau) \left( 
 |U(\tau)v_-|^2U(\tau)v_-\right)\dd \tau.
 \end{multline*}
 \item For any $0<\eta<2$, we have, in $L^2$ and as
 $\eps\to 0$,
 \begin{equation*}
 \left(W_+^{\rm mod}\right)^{-1} (\eps v_0) = \eps v_ 0 + \eps^3 \mu_2+
 \O\left(\eps^{5-\eta}\right), 
 \end{equation*}
 where $\mu_2=\mu_2(v_0)\in L^2$ is defined by
 \begin{multline*}
 \mu_2 = -\irm \lambda \int_0^1 \left( U(-s)\left(
 |U(s)v_0|^2U(s)v_0\right)\right)\dd s 
 \\ 
 +\lambda\int_1^\infty 
 \left( M(-\tau)\F^{-1}\left(\left|\hat{M(\tau)v_0}\right|^2
 \hat{M(\tau)v_0} \right)- \F^{-1}\left(|\hat v_0 |^2\hat
 v_0\right)\right) \frac{\dd \tau}{\tau} ,
 \end{multline*}
 and $M(t)$ stands for the multiplication by
 $\expe^{^\irm\frac{x^2}{2t}}$. 
 \item For any $0<\eta<2$, we have, in $L^2$ and as $\eps\to 0$,
 \begin{equation*}
 \Scatt(\eps v_-) = \eps v_- + \eps^3 \nu_2 + \O\left(\eps^{5-\eta}\right),
 \end{equation*}
 where $\nu_2\in L^2$ is defined by
 $ \nu_2 = w_2 + \mu_2(v_-),$
 where $w_2$ is given by the first point, and $\mu_2$ by the
 second. 
 \end{itemize}
\end{theorem}
\begin{remark}
 The functions $\mu_2$ and $\nu_2$ are more regular than merely
 $L^2$, as we will see in Section~\ref{sec:complete-higher} that they
 belong to $H^\gamma\cap \F(H^\gamma)$ for any $0<\gamma<1$.
\end{remark}
\looseness-1
It may be surprising that the leading order behavior of the modified
wave and scattering operators at the origin is the identity: we recall
that from their definition, these operators are already nontrivial,
and account for long range effects. We also emphasize that the first
corrector has a more involved expression than in the short range case,
a case where we would simply have (typically for the two-dimensional
and three-dimensional
cubic Schr\"odinger equations, see e.g.~\cite{CaGa09})
\begin{align*}
 w_2& =-\irm\lambda \int_{-\infty}^0 \left( U(-s)\left(
 |U(s)v_0|^2U(s)v_0\right)\right)\dd s ,\\
 \mu_2& = -\irm\lambda \int_0^\infty \left( U(-s)\left(
 |U(s)v_0|^2U(s)v_0\right)\right)\dd s . 
\end{align*}
Note that in the (one-dimensional) long range case, these integrals diverge. 
\smallbreak

As evoked above, in principle, the method of proof that we present
allows to compute (some) higher order terms in the asymptotic
expansions near the origin. 

\subsection{Outline}

In Section~\ref{sec:prelim}, we recall several properties which are
classical in the context of scattering theory for nonlinear
Schrödinger equations, and which are of constant use in this
paper. In Section~\ref{sec:wave-op}, we revisit the construction of
the modified wave operator $\W$, in such a way that the leading order
behavior of 
this operator at the origin, presented in Section~\ref{sec:wave-zero}, is rather
straightforward. In Section~\ref{sec:wave-zero}, we also consider the
first corrector in the asymptotic expansion of $\W$ at the origin, an
aspect which requires some extra work. In Section~\ref{sec:complete},
we recall the modified asymptotic completeness result established in
\cite{HN98}, and infer the leading of behavior of $\Scatt$ at the
origin. The first corrector is derived in
Section~\ref{sec:complete-higher}, where the last two error estimates
announced in Theorem~\ref{theo:main} are proved. 
\subsection{Notations}
We recall the classical factorization of the Schrödinger group, 
\begin{equation*}
 U(t)= \expe^{\irm \frac{t}{2}\d_x^2}= M(t)D(t)\F M(t),
\end{equation*}
where the multiplication $M(t)$, the dilation $D(t)$ and
the Fourier transform $\F$ are defined by
\begin{equation}\label{eq:MDF}
 \begin{aligned}
 & M(t) = \expe^{\irm \frac{ x^2}{2t}}, \quad D(t) f(x) =
 \frac{1}{(\irm t)^{1/2}}f\left(\frac{x}{t}\right),\\
 &{\F} f(\xi)=\hat
 f(\xi)=\frac{1}{\sqrt{2\pi}}\int_{\R}f(x)\expe^{-\irm x \xi} \dd x.
 \end{aligned}
 \end{equation}
Note that each of these three operators is unitary on $L^2(\R)$, and
that~\eqref{eq:asym-lin} reads $U(t)\approx M(t)D(t)\F $, see
also Lemma~\ref{lem:R}.
We recall that the space $\Sigma$ is defined by
\begin{equation*}
 \Sigma=H^1\cap \F(H^1)=\left\{f\in H^1(\R);\
 \|f\|_\Sigma\coloneqq \|f\|_{L^2}+\|\d_x f\|_{L^2} + \|x
 f\|_{L^2}<\infty\right\}.
\end{equation*}
We note that $\Sigma$ is a Banach algebra, invariant under the Fourier
transform, and such that $\Sigma\hookrightarrow L^1\cap L^\infty$. In
addition, $\Sigma$ is invariant under the action of $U(t)$, for any
$t\in \R$.
\smallbreak

For functions $f^\eps,g^\eps\ge 0$ depending on time $t$ and $\eps$, the
notation
\begin{equation*}
 f^\eps\lesssim g^\eps
\end{equation*}
means that there exists $C$ independent of $t$ and $\eps$ such that
\begin{equation*}
 f^\eps\le C g^\eps. 
\end{equation*}




\section{Technical preliminaries}
\label{sec:prelim}

In this section, we gather classical estimates which can be found in
several references cited in the introduction. 


\begin{lemma}\label{lem:J}
The operator
\[J(t)=x+\irm t\d_x\]
satisfies the following properties:
 \begin{itemize}
 \item $J(t) =U(t)xU(-t)$, and therefore $J$ commutes with the
 linear part of~\eqref{eq:NLS},
\begin{equation}\label{eq:Jcommute}
\left[ J(t),\irm \d_t +\frac{1}{2}\d_x^2\right]=0\, .
\end{equation}
\item It can be factorized as
 \begin{equation*}%\label{eq:Jfactor}
J(t)= \irm t \expe^{\irm\frac{x^2}{2t}}\d_x\Big( \expe^{-\irm\frac{x^2}{2t}}\,
\cdot\Big)\, .
\end{equation*}
As a consequence, $J$ yields weighted Gagliardo--Nirenberg
inequalities. For
$2\le r \le \infty$, there exists
$C(r)$ depending only on $r$ such that
\begin{equation}\label{eq:GNlibre}
\left\| f \right\|_{L^r}\le \frac{C(r)}{|t|^{\delta(r)}} \left\|
f \right\|_{L^2}^{1-\delta(r)} \left\|
J(t) f \right\|_{L^2}^{\delta(r)},\quad \delta(r)\coloneqq 
\frac{1}{2}-\frac{1}{r}\, .
\end{equation}
Also, if $F(z)=G(|z|^2)z$ is $C^1$, then $J(t)$
acts like a derivative on $F(w)$:
\begin{equation}\label{eq:Jder}
J(t)\left(F(w)\right) = \d_z F(w)J(t)w -\d_{\bar z} F(w)\bar{ J(t)w
}\, .
\end{equation}
\end{itemize}
\end{lemma}


\begin{lemma}\label{lem:R}
 Denote by
 \begin{equation*}
 \cR(t)= M(t)D(t)\F \left(M(t)-1\right)\F^{-1},
 \end{equation*}
 where the above terms are defined in~\eqref{eq:MDF}. The following
 estimates hold for $|t|\ge 1$:
\begin{enumerate}
\item\label{lemma5_1} For all $s>1/2$ and $0\le \theta \le 1$,
 \begin{equation*}
 \|\cR(t)f\|_{L^\infty_x}\lesssim
 \frac{1}{|t|^{1/2+\theta}}\|f\|_{H^{s+2\theta}},\quad \forall
 f\in \Sch(\R).
 \end{equation*}
\item\label{lemma5_2} For all $0\le \theta \le 1$,
 \begin{equation*}
 \|\cR(t)f\|_{L^2_x}\lesssim
 \frac{1}{|t|^{\theta}}\|f\|_{H^{2\theta}},\quad \forall
 f\in \Sch(\R).
 \end{equation*} 
\item\label{lemma5_3} For all $0\le \theta \le 1$,
 \begin{equation*}
 \|J(t)\cR(t)f\|_{L^2_x}\lesssim
 \frac{1}{|t|^{\theta}}\|f\|_{H^{1+2\theta}},\quad \forall
 f\in \Sch(\R).
 \end{equation*}
\item\label{lemma5_4} We have
 \begin{equation*}
 \|\d_x \cR(t)f\|_{L^2_x}\lesssim
 \frac{1}{|t|}\|f\|_{H^1}+
 \frac{1}{|t|}\|x f\|_{H^{2}},\quad \forall
 f\in \Sch(\R).
 \end{equation*} 
\end{enumerate}
\end{lemma}
\begin{remark}
 We emphasize that $\cR(t)$ does not map $L^\infty(\R)$ into itself,
 since $\F \left(M(t)-1\right)\F^{-1}= U\big(\frac{1}{t}\big)-1$ and 
 $U\big(\frac{1}{t}\big)\big(\expe^{-\irm t\frac{x^2}{2}}\big)=\delta_0$. 
\end{remark}
Even though such estimates can be found in the existing literature,
we give a proof, as the ideas will be resumed when $M(t)-1$ is present. 
\begin{proof}
 In view of the definition of $\cR$, we readily have
 \begin{equation*}
 \|\cR(t)f\|_{L^\infty_x} =\frac{1}{|t|^{1/2}}\|\F (M(t)-1)\F^{-1}f\|_{L^\infty}.
 \end{equation*}
 Using Hausdorff--Young inequality and the general estimate
 \begin{equation}\label{eq:sinus}
 |M(t)-1| = 2\left|\sin\left(\frac{x^2}{4t}\right)\right|\lesssim
 \left|\frac{x^2}{t}\right|^\theta,\quad \forall \, 0\le \theta\le 1,
 \end{equation}
 we find
 \begin{equation*}
 \|\cR(t)f\|_{L^\infty_x} \lesssim \frac{1}{|t|^{1/2}}\left\|
 \left|\frac{x^2}{t}\right|^\theta \F^{-1}f\right\|_{L^1}. 
 \end{equation*}
 Using the easy property $\F(H^s)\hookrightarrow L^1$ for
 $s>1/2$ (the same is true on $\R^d$ provided that $s>d/2$,
 from Cauchy--Schwarz inequality),
 \begin{equation*}
 \|\cR(t)f\|_{L^\infty_x} \lesssim \frac{1}{|t|^{1/2+\theta}}\left\|
 \left\langle x\right\rangle ^s|x|^{2\theta} \F^{-1}f\right\|_{L^2}\lesssim
 \frac{1}{|t|^{1/2+\theta}} \|f\|_{H^{s+2\theta}}, 
 \end{equation*}
 hence the first inequality. For the second one, since $D(t)$ is
 unitary on $L^2$, Plancherel formula
 and~\eqref{eq:sinus} yield
 \begin{equation*}
 \|\cR(t)f\|_{L^2_x} = \| (M(t)-1)\F^{-1}f\|_{L^2}\lesssim \frac{1}{|t|^{\theta}} \left\|
 |x|^{2\theta} \F^{-1}f\right\|_{L^2} \lesssim
 \frac{1}{|t|^{\theta}} \|f\|_{H^{2\theta}}.
 \end{equation*}
For the third inequality, we use the formula $J(t)=U(t)xU(-t)$, along
with the factorization
\begin{equation}\label{eq:U(-t)}
 U(-t) = iM(-t) \F^{-1}D\left(\frac{1}{t}\right) M(-t),
\end{equation}
to obtain
 \begin{equation*}
 \|J(t)\cR(t)f\|_{L^2_x} = \|x (M(t)-1)\F^{-1}f\|_{L^2}\lesssim
 \frac{1}{|t|^{\theta}} \left\| 
 |x|^{1+2\theta} \F^{-1}f\right\|_{L^2} \lesssim 
 \frac{1}{|t|^{\theta}} \|f\|_{H^{1+2\theta}}.
 \end{equation*}
 To estimate $\d_x \cR(t)f$, two terms appear, whether the derivative
 hits the first factor $M(t)$ or not, and we have
 \begin{equation*}
 \|\d_x\cR(t)f\|_{L^2_x} \le \left\|\frac{x}{t}D(t)\F
 (M(t)-1)\F^{-1}f\right\|_{L^2}+ \frac{1}{|t|}\left\|\d_x\F
 (M(t)-1)\F^{-1}f\right\|_{L^2}.
 \end{equation*}
 The first term on the right hand side is equal to
 \begin{align*}
 \left\|x\F
 (M(t)-1)\F^{-1}f\right\|_{L^2}&= \left\|
 \d_x (M(t)-1)\F^{-1}f\right\|_{L^2}\\
 &\le \left\| \frac{x}{t}\F^{-1}f\right\|_{L^2}
 +\left\|
 (M(t)-1) \d_x\F^{-1}f\right\|_{L^2}\\
 & \lesssim \frac{1}{|t|}
 \|f\|_{H^1} + \frac{1}{|t|^{\theta}}\|xf\|_{H^{2\theta}},
 \end{align*}
 for any $\theta\in [0,1]$. Note that 
 \[ \frac{1}{|t|}\left\|\d_x\F
 (M(t)-1)\F^{-1}f\right\|_{L^2}= \frac{1}{|t|}\left\|x
 (M(t)-1)\F^{-1}f\right\|_{L^2}\lesssim \frac{1}{|t|}\left\|x
 \F^{-1}f\right\|_{L^2},
\]
hence the result, by choosing $\theta=1$.
\end{proof}




\section{Modified wave operator}
\label{sec:wave-op}
A key aspect in the proofs of the existence of modified wave
operators for~\eqref{eq:NLS}, in particular to get minimal regularity
assumptions on the asymptotic state $u_-$, is to consider a suitable
approximate solution near $t=-\infty$. Introduce two such approximate
solutions, $u_1$ and $u_2$, defined for $t\le -10$, given by
\begin{equation*}
 u_1(t,x) \coloneqq \frac{1}{(\irm t)^{1/2}}\expe^{\irm \frac{x^2}{2t}}\fum
 \left(\frac{x}{t}\right) \expe^{\irm S_-(t,x)}=M(t)D(t)\hat w(t,x),
\end{equation*}
where
\begin{equation*}
\hat w(t)\coloneqq \fum \expe^{\irm \lambda |\fum|^2\log|t|},
\end{equation*}
where we resume the same notations as in~\cite{HN06}, 
and
\begin{equation*}
 u_2(t) \coloneqq U(t)\F^{-1} \hat w.
\end{equation*}
As $t\to -\infty$, $u_1$ and $u_2$ are close, in view of
Lemma~\ref{lem:R}, since $u_2(t)-u_1(t)= \cR(t)\hat w$.
Our goal is to keep track of the smallness of the asymptotic state as precisely as
possible, in the sense that we consider $u_-^\eps=\eps v_-$ for some
fixed $v_-$, assume that $\eps>0$ is small, and we want to get larger
powers of $\eps$ in the error terms whenever we can. It turns out that
apparently, $u_2$ is a better candidate than $u_1$ in this
direction. More explicitly, the Duhamel formula for $u-u_1$ (see
\cite{HN06}) contains a linear term ($R\hat w$), which we want to
remove. We first gather some estimates on $u_1$.

\begin{lemma}\label{lem:u1}
 Let $\um\in \H$: $\|u_1(t)\|_{L^2}=\|u_2(t)\|_{L^2}=\|\um\|_{L^2}$ for all $t\le
 -10$. 
 There exists $C>0$ independent of $\um$ such that for all $t\le -10$,
 \begin{align*}
 \| u_1(t)\|_{L^\infty}&\le
 \frac{\|\fum\|_{L^\infty}}{\sqrt{|t|}},\\
 \|\d_x u_1(t)\|_{L^2}&\le C\left(\|\um\|_{\H}
 +\|\um\|_{\H}^3\right) ,\\
 \|J(t) u_1(t)\|_{L^2}&\le C \|\um\|_{\H}
 \left(1+ \|\um\|_{\H}^2\log |t|\right). 
 \end{align*}
\end{lemma}
\begin{proof}
 The conservation of the $L^2$-norm is obvious, as well as 
the estimate 
\begin{equation*}
 \|u_1(t)\|_{L^\infty}\le \frac{\|\fum\|_{L^\infty}}{\sqrt{|t|}}.
\end{equation*}
The expression of $u_1$ yields directly
\begin{equation*}
 \|\d_x u_1(t)\|_{L^2} \le \|\d_x \fum\|_{L^2} +
 |\lambda|\frac{\log |t|}{|t|}\|\fum\|_{L^\infty}^2\|\d_x \fum\|_{L^2} .
\end{equation*}
In view of the factorization for $J$ stated in Lemma~\ref{lem:J},
\begin{align*}
 J(t)u_1(t)&= \irm t M(t) \d_x \left( D(t)\fum
 \expe^{\irm S_-(t)}\right)\\
&= \expe^{\irm \frac{x^2}{2t}+\irm S_-(t)}\left(\irm D(t)(\d_x \fum)-\lambda 
 D(t)\fum \times \log|t| \times\d_x \left|\fum\right|^2\left(\frac{x}{t}\right)\right),
\end{align*}
thus
\begin{equation*}
 \| J(t)u_1(t)\|_{L^2}\lesssim \|\fum\|_{H^1} + 
\log |t|\|\fum\|_{L^\infty}^2 \|\fum\|_{H^1} .
\end{equation*}
The lemma follows. 
\end{proof}

\begin{lemma}\label{lem:Rw}
 Let $\um\in \H$. 
 There exists $C>0$ independent of $\um$ such that for all $t\le -10$,
 \begin{align*}
 \| \cR(t)\hat w\|_{L^\infty}&\le
 \frac{C}{\sqrt{|t|}} \left(\|\fum\|_{H^1}+ \|\fum\|_{H^1}^3\right),\\
 \|\cR(t)\hat w\|_{L^\infty}&\le \frac{C}{|t|^{3/2}}\left(\|
 \um\|_{\H}+ (\log|t|)^{3}\|\um\|_{\H}^7\right).\\
 \| U(-t) \cR(t)\hat w\|_{\Sigma}&\le \frac{C}{|t|}\left( \|\um\|_{\H} +
 (\log|t|)^{3}\| \um\|_{\H}^7\right).
 \end{align*}
\end{lemma}
\begin{proof}
 If we use
the first point of Lemma~\ref{lem:R} with $\theta=0$ and
$s>1/2$, the phase factor causes the appearance of a $\log|t|$ term,
since the $H^1$-norm of $\hat w$ is involved. So we rather take
$\theta>0$ and $s>1/2$ such that $s+2\theta=1$, and so,
for $t\le -10$,
\begin{align*}
 \|\cR(t)\hat w\|_{L^\infty}\lesssim \frac{\|\hat
 w\|_{H^1}}{|t|^{1/2+\theta}}& \lesssim \frac{ 1}{|t|^{1/2+\theta}}
\left(\|\fum\|_{H^1}+(\log|t|)\|\fum\|_{L^\infty}^2\|\fum\|_{H^1}\right)\\ 
&\lesssim \frac{1}{\sqrt{|t|}}\left(\|
 \fum\|_{H^1}+\|\fum\|_{H^1}^3\right),
\end{align*}
where the logarithmic factor was left out since $\theta>0$.
 It is actually at the level of this estimate that
the $H^1$-norm of $\fum$ appears, instead of merely its
$L^\infty$-norm.

To obtain the stronger time decay, we choose $s=\theta=1$ in the first
point of Lemma~\ref{lem:R}, hence, since $H^3(\R)$ is
an algebra, 
\begin{align*}
 \|\cR(t)\hat w\|_{L^\infty}\lesssim \frac{\|\hat
 w\|_{H^3}}{|t|^{3/2}}& \lesssim \frac{ 1}{|t|^{3/2}}
\left(\|\fum\|_{H^3}+(\log|t|)^3\|\fum\|_{H^3}^7\right).
\end{align*}

In view of Lemma~\ref{lem:J},
\begin{equation*}
 \|U(-t) \cR\hat w\|_{\Sigma}= \|\cR\hat w\|_{L^2}+\|\d_x
 \cR\hat w\|_{L^2} + \|J(t)\cR\hat w\|_{L^2},
\end{equation*}
and Lemma~\ref{lem:R} yields
\begin{equation*}
 \|U(-t) \cR(t)\hat w\|_{\Sigma}\lesssim \frac{1}{|t|}\|\hat
 w\|_{H^3} + \frac{1}{|t|}\|x\hat
 w\|_{H^2} .
\end{equation*}
Using the fact that $H^s(\R)$ is an algebra for $s>1/2$, we infer
\begin{align*}
 \|U(-t) \cR(t)\hat w\|_{\Sigma}&\lesssim
 \frac{1}{|t|}\left(\|\fum\|_{H^3} + (\log|t|)^3
 \|\fum\|_{H^3}^7 \right) %\\
%&\quad 
+ \frac{1}{|t|}\left( \|x\fum\|_{H^2} +
 (\log|t|)^2\|x\fum\|_{H^2}^5\right)\\
&\lesssim \frac{1}{|t|}\left( \|\um\|_{\H} +
 \sum_{j=0,1}(\log|t|)^{2+j}\| \um\|_{\H}^{5+2j}\right),
\end{align*}
hence the lemma. 
\end{proof}

\begin{proposition}\label{prop:wave-op}
There exist $\delta_0>0$, $C$, and a polynomial $P$, such that the following holds.
 For any $u_-\in L^2$ with $\um\in \H$, where $\H$ is
 defined in~\eqref{eq:H}, and such that $\|
 \fum\|_{H^1}\le \delta_0$, there exists a unique
 $u_0\in \Sigma$ such that the solution $u\in C(\R,\Sigma)$
 to~\eqref{eq:NLS} with 
 $u_{\mid t=0}=u_0$ satisfies, for $t\le -10$,
 \begin{equation*}
 \left\|U(-t)\left(u(t)-u_2(t)\right)\right\|_{\Sigma}\le
 C\|\um\|_{\H}^3P\left(\|\um\|_{\H}\right) \frac{(\log|t|)^4}{|t|}. 
 \end{equation*}
\end{proposition}
\begin{remark}
 The choice of the approximate solution $u_2$ is tailored to ensure
that, with the proof presented below, the error term is superlinear in
some norm of $u_-$, uniformly in the limit $t\to
-\infty$. 
On the other hand, the smallness assumption, which is a consequence of
the first estimate in 
Lemma~\ref{lem:Rw}, is stronger than in the
cited references, since Sobolev embedding yields $H^1(\R)\subset
L^\infty(\R)$. 
\end{remark}
\begin{proof}
 Like in~\cite{Ozawa91} or~\cite{HN06}, the proof relies on a fixed
 point argument near $t=-\infty$. Resuming the computations from
~\cite{HN06} (and keeping the same notations), 
 \begin{equation*}
 \irm \d_t \hat w =\frac{\lambda}{t}|\hat w|^2\hat w,
 \end{equation*}
and so 
we check that $u_2$ solves the equation
 \begin{equation*}
 \irm \d_t \F U(-t) u_2 = \lambda \F
 U(-t)\left(\frac{1}{t}U(t)\F^{-1}\left(|\hat w|^2\hat w\right)\right).
 \end{equation*}
Therefore, we want to solve
\begin{equation*}
 \irm \d_t \F U(-t)\left( u-u_2\right) = \lambda \F
 U(-t)\left(|u|^2 u-\frac{1}{t}U(t)\F^{-1}\left(|\hat w|^2\hat w\right)\right).
\end{equation*}
We insert the term $|u_2|^2 u_2$ into the above equation,
and use the identity
\begin{align*}
 \frac{1}{t}U(t)\F^{-1}\left(|\hat w|^2\hat w\right)&= \frac{1}{t}MD\F M\F^{-1}
 \left(|\hat w|^2\hat w\right) \\
&= \frac{1}{t}\cR(t)\left( |\hat w|^2\hat w\right)+
 \frac{1}{t}M(t)D(t)\left(|\hat w|^2\hat w\right), 
\end{align*}
with 
\begin{align*}
 \frac{1}{t}M(t)D(t)\left(|\hat w|^2\hat w\right)= M(t)|M(t)D(t)\hat
 w|^2D(t)\hat w = |u_1|^2u_1,
\end{align*}
so Duhamel's formula reads, along with the requirement
$U(-t)\left(u(t)-u_2(t)\right)\to 0 $ as $t\to -\infty$,
\begin{multline}\label{eq:Duhamel1}
 u(t)=u_2(t) - \irm \lambda\int_{-\infty}^t
 U(t-\tau)\left(|u|^2u- |u_2|^2 u_2\right)(\tau)\,\dd \tau
%\\
%\begin{aligned}[b]
+\irm \lambda \int_{-\infty}^t
 U(t-\tau)\left( |u_2|^2 u_2-|u_1|^2u_1\right)(\tau)\, \dd \tau\\
%& 
-\irm \lambda \int_{-\infty}^t
 U(t-\tau)\left(\frac{1}{\tau}\cR(\tau)\left(|\hat w|^2\hat w\right) \right)(\tau)\, \dd \tau.
% \end{aligned}
\end{multline}
The last two lines correspond to source terms, involving only the
various approximate solutions, and can be estimated thanks to
Lemma~\ref{lem:R}, as $u_2 = u_1+\cR\hat w$. 

Denote by $\Phi(u)= u_2+\Phi_1(u)+\Phi_2+\Phi_3$ the right hand side of
\eqref{eq:Duhamel1}, where $\Phi_j$ corresponds to the $j$-th
line (note that $\Phi_2$ and $\Phi_3$ do not depend on $u$). We have
\begin{align*}
 \left|\d_x\left(|u_2|^2 u_2-|u_1|^2u_1\right)\right|&= \left|\d_x\left(|\cR\hat w
 +u_1|^2 (\cR\hat w+u_1)-|u_1|^2u_1\right)\right| \\
&\lesssim \left(|\cR\hat w|^2+|u_1|^2\right) |\d_x \cR\hat w| + \left(|\cR\hat w|+|u_1|\right)
 |\cR\hat w| |\d_x u_1|. 
\end{align*}
Using the formula $J(t)=U(t)xU(-t)$ from Lemma~\ref{lem:J}, as well as
\eqref{eq:Jder} that state that $J(t)$ can be thought of as a
derivative here, we find
\begin{multline*}
 \left\| U(-t)\Phi_2(t)\right\|_{\Sigma}\lesssim \int_{-\infty}^t\left(\|\cR(\tau)\hat
 w\|_{L^\infty}^2+\|u_1(\tau)\|_{L^\infty}^2\right) \|U(-\tau)\cR\hat
 w\|_{\Sigma}\, \dd \tau\\
 + \int_{-\infty}^t\left(\|\cR(\tau)\hat w\|_{L^\infty}+\|u_1(\tau)\|_{L^\infty}\right) 
 \|\cR(\tau)\hat w\|_{L^\infty} \|U(-\tau) u_1(\tau)\|_{\Sigma}\, \dd \tau.
\end{multline*}
We invoke Lemmas~\ref{lem:u1} and~\ref{lem:Rw} to estimate the above
terms. More precisely, then factor $\|\cR(\tau)\hat w\|_{L^\infty}
$ is controlled by the second estimate of Lemma~\ref{lem:Rw}, and we
get, for $\|\fum\|_{H^1}\le 1$ and $t\le -10$:
\begin{equation*}
 \|U(-t)\Phi_2(t)\|_\Sigma\lesssim
 \|\fum\|_{H^1}\|\um\|_{\H}^2P_2\left( 
 \|\um\|_{\H}\right) \int_{-\infty}^t \frac{(\log |\tau|)^4}{\tau^2}\dd \tau,
\end{equation*}
for some polynomial $P_2$. 
The term $\Phi_3$ is estimated thanks to Lemma~\ref{lem:R}, and the
choice of parameter is motivated by the previous estimate, keeping in
mind that a factor $\frac{1}{\tau}$ is already present in the
definition of $\Phi_3$. For the $L^2$-estimate, we thus choose
$\theta=1$ in Lemma~\ref{lem:R}, so that
\begin{align*}
 \|\Phi_3(t)\|_{L^2}&\lesssim \int_{-\infty}^t \left\| |\hat w|^2\hat
 w\right\|_{H^2}\frac{\dd \tau}{\tau^2 }\lesssim \int_{-\infty}^t \|
 \hat w \|_{L^\infty}^2 \| \hat
 w\|_{H^2}\frac{\dd \tau}{\tau^2 }\\
 &\lesssim \|\fum\|_{L^\infty}^2\|\fum \|_{H^2}
 \int_{-\infty}^t
 \left(1+ (\log|\tau|)^2 \|\fum \|_{H^2}^4\right) \frac{\dd \tau}{\tau^2 }.
\end{align*}
The norm $\|xU(-t)\Phi_3(t)\|_{L^2}=\|J(t)\Phi_3(t)\|_{L^2}$ is
estimate similarly, by setting again $\theta=1$ (in the third point of
Lemma~\ref{lem:R}, and we get, since the $H^2$-norm is replaced by an
$H^3$-norm in the above computation,
\begin{equation*}
 \|J(t)\Phi_3(t)\|_{L^2}\lesssim \|\fum\|_{L^\infty}^2\|\fum \|_{H^3}
 \int_{-\infty}^t
 \left(1+ (\log|\tau|)^3\|\fum \|_{H^3}^6\right) \frac{\dd \tau}{\tau^2 }.
 \end{equation*}
 The choice of suitable parameters for the $\dot H^1$-norm was already
 made in the statement of Lemma~\ref{lem:R}, and we have
 \begin{equation*}
 \|\d_x\Phi_3(t)\|_{L^2}\lesssim 
 \int_{-\infty}^t
 \left(\||\hat w|^2 \hat w\|_{H^1} + \|x|\hat w|^2 \hat w\|_{H^2}\right)
 \frac{\dd \tau}{\tau^2 }. 
 \end{equation*}
 The term involving the $H^1$-norm has been estimated above, so we
 focus on the other norm. Direct computations yield
 \begin{equation*}
 \left| \d_x^2\left( x|\hat w|^2 \hat w\right)\right| \lesssim 
 |x||\hat w|^2 |\d_x^2\hat w| + |\hat w|^2 |\d_x\hat w| +
 |x||\hat w||\d_x\hat w|^2,
 \end{equation*}
 and we get
 \begin{multline*}
 \|\d_x\Phi_3(t)\|_{L^2} \lesssim 
 \|\fum\|_{L^\infty}^2 \|\fum\|_{H^1} \int_{-\infty}^t
 \left(1+ (\log|\tau|)\|\fum\|_{H^1}^2\right)
 \frac{\dd \tau}{\tau^2 }
 \\
 +\|\fum\|_{L^\infty}\|\um\|_{\H}^2\int_{-\infty}^t
 \left(1+ (\log|\tau|)^2 \|\um\|_{\H}^4\right) \frac{\dd \tau}{\tau^2 }. 
 \end{multline*}
 The source term is therefore controlled, for $t\le -10$, by:
 \begin{equation}\label{eq:sourcePhi}
 \|U(-t)(\Phi_2(t)+\Phi_3(t))\|_{\Sigma}\lesssim
 \|\fum\|_{H^1}\|\um\|_{\H}^2P\left(\|\um\|_{\H}\right)\int_{-\infty}^t
 \frac{\left(\log|\tau|\right)^4}{\tau^2}\dd \tau,
 \end{equation}
 for some polynomial $P$ whose precise expression is irrelevant.

 The end of the proof consists of a fixed point argument. Mimicking
~\cite{Ozawa91}, for $\alpha\in \mathopen{]}1/2,1[$ arbitrary (but fixed),
 and $T\gg 1$, introduce the
 space
 \begin{multline*}
 X^\alpha(T)=\Bigl\{ u\in L^\infty(\mathopen{]}-\infty,-T],L^2) \text{ such that }t\mapsto U(-t)u(t)\in
 L^\infty(\mathopen{]}-\infty,-T],\Sigma), 
 \\
 \sup_{t\le
 -T}|t|^\alpha\|U(-t)(u(t)-u_2(t))\|_{\Sigma}\le \|\fum\|_{H^1}\Bigr\}.
 \end{multline*}
 and define on $ X^\alpha(T)$ the metric
 \begin{equation*}
 d(u,v) = \sup_{t\le
 -T}|t|^\alpha\|u(t)-v(t)\|_{L^2}.
 \end{equation*}
 Note that here, we choose to measure distance by considering the
 $L^2$ norm only: $X^\alpha(T)$, equipped with this distance, is
 a complete metric space. 
 
 For $u\in X^\alpha(T)$, $\Phi_1(u)$ is estimated by
 \begin{align*}
 \|U(-t)\Phi_1(u(t))\|_{L^2}&\lesssim \int_{-\infty}^t
 \left(\|u(\tau)\|_{L^\infty}^2+
 \|u_2(\tau)\|_{L^\infty}^2\right)\|u(\tau)-u_2(\tau)\|_{L^2}\, \dd \tau \\
 &\lesssim
 \|\fum\|_{H^1}^2\int_{-\infty}^t\frac{\dd \tau}{\tau^{1+\alpha}}\lesssim \frac{ \|\fum\|_{H^1}^2}{|t|^\alpha}.
 \end{align*}
 For $A\in\{\d_x,J(t)\}$, we have
 \begin{equation}\label{eq:est-gen}
 \left|A\left( |u|^2u-|u_2|^2u_2\right)\right| \lesssim |u|^2 |A(u-u_2)|+
 \left(|u|+|u_2|\right)|A u_2||u-u_2|. 
 \end{equation}
 Therefore,
 \begin{align*}
 \|U(-t)\Phi_1(u(t))\|_{\Sigma}&\lesssim \int_{-\infty}^t
 \|u(\tau)\|_{L^\infty}^2\|U(-\tau)(u(\tau)-u_2(\tau))\|_{\Sigma}\, \dd \tau
 \\
&\hspace*{15mm} + \int_{-\infty}^t
 (\|u(\tau)\|_{L^\infty} +\|u_2(\tau)\|_{L^\infty})\|U(-\tau)u_2(\tau)\|_{\Sigma}
 \|u(\tau)-u_2(\tau)\|_{L^\infty}
 \dd \tau\\
 &\lesssim \|\fum\|_{H^1}^2P_3\left(\|\um\|_{\H}\right)\int_{-\infty}^t
 \frac{\log|\tau|}{\tau^{1+\alpha}}\, \dd \tau ,
 \end{align*}
 for some polynomial $P_3$ whose precise expression is irrelevant,
 where we have used Lemmas~\ref{lem:u1} and~\ref{lem:Rw} to estimate 
 $\|U(-\tau)u_2(\tau)\|_{\Sigma}$. The 
 $L^\infty$-norm of $u(\tau)-u_2(\tau)$ is controlled by
 $|\tau|^{-1/2-\alpha}$ in view of Lemma~\ref{lem:J} and the
 definition of $X^\alpha(T)$. Therefore, by choosing $T$
 sufficiently large, we check that $\Phi$ maps $X^\alpha(T)$ to
 itself.

 To conclude by a fixed point argument, we show that
 $\Phi$ is a contraction provided
 that $\|\fum\|_{H^1}$ is sufficiently small. Indeed, for $u,v\in
 X^\alpha(T)$, 
 \begin{equation*}
 \Phi(u)(t)-\Phi(v)(t) =- \irm \lambda \int_{-\infty}^t
 U(t-\tau)\left(|u|^2u-|v|^2v\right)(\tau)\, \dd \tau, 
 \end{equation*}
 and so, for $t\le -10$,
 \begin{align*}
 \left\|\Phi(u)(t)-\Phi(v)(t)\right\|_{L^2}&\lesssim \int_{-\infty}^t
 \left(\|u(\tau)\|_{L^\infty}^2+ \|v(\tau)\|_{L^\infty}^2\right)\|u(\tau)-v(\tau)\|_{L^2}\, \dd \tau
 \\
 &\lesssim \|\fum\|_{H^1}^2 d(u,v) \int_{-\infty}^t\frac{\dd \tau}{\tau^{1+\alpha}},
 \end{align*}
 where we have used
 \begin{align*}
 \|u(\tau)\|_{L^\infty}+ \|v(\tau)\|_{L^\infty}& \le
 2\|u_2(\tau)\|_{L^\infty} + C\|u(\tau)-u_2(\tau)\|_{H^1}+
 C\|v(\tau)-u_2(\tau)\|_{H^1}\\
 &\lesssim \frac{\|\fum\|_{H^1} }{\sqrt{|\tau|}},
 \end{align*}
 since $\alpha>1/2$. Therefore,
 \begin{equation*}
 d\left(\Phi(u),\Phi(v)\right)\lesssim \|\fum\|_{H^1}^2 d(u,v), 
 \end{equation*}
 hence the result provided that $\|\fum\|_{H^1}$ is
 sufficiently small.
 \end{proof}


\section{Behavior of the modified wave operator near the origin}
\label{sec:wave-zero}

\subsection{Leading order asymptotic behavior}


Let $u_-^\eps=\eps v_-$ with $ v_-\in \H$ independent of
$\eps$. Denote by $u^\eps$ the solution provided by
Proposition~\ref{prop:wave-op}. Then like $u_-^\eps$, $u^\eps$ is of order
$\eps$, and the remainder in Proposition~\ref{prop:wave-op} is
$\O(\eps^3)$. The long range phase correction $S_-^\eps$ is 
$\O(\eps^2\log|t|)$, for $t\le -10$: its contribution is negligible
for times $t$ such that $\eps^2|\log|t||\ll 1$, so in particular we can match with a linear
solution at $t^\eps_\gamma = -1/\eps^\gamma$,
\begin{equation*}
 \irm \d_t \u^\eps_1 +\frac{1}{2}\d_x^2 \u^\eps_1 =0 ,
\end{equation*}
with
\begin{equation*}
 \u^\eps_{1\mid
 t=-1/\eps^\gamma} = 
 \frac{1}{(\irm t)^{1/2}}\expe^{\irm \frac{x^2}{2t}}\F(u_-^\eps)\left(\frac{x}{t}\right)
 \Big|_{t=-1/\eps^\gamma}. 
\end{equation*}
The right hand side is $M(t)D(t)\F u_-^\eps$ evaluated at $t = -1/\eps^\gamma$. 
In view of Lemma~\ref{lem:R}, up to a small error term, 
\begin{equation*}
 \u_1^\eps(t) =
 U(t)u_-^\eps, \text{ or, equivalently, }\u_1^\eps=\eps v_1,
 \text{ with }v_1(t)=U(t)v_-.
\end{equation*}
We choose this definition for the first order approximation
$\u_1^\eps$. 
Let $t\ge -1/\eps^\gamma$. Duhamel's formula implies
\begin{multline*}
 U(-t)\left(u^\eps(t)- v^\eps_1(t)\right) \\
 = U\left(-s\right)\left(
 u^\eps(s)- u^\eps_2(s)+u_2^\eps(s)
-\u^\eps_1(s)\right)\Big|_{s=-1/\eps^\gamma} 
 -\irm \lambda
 \int_{-1/\eps^\gamma}^t U(-\tau)\left(|u^\eps|^2u^\eps\right)(\tau)\, \dd \tau,
\end{multline*}
that is
\begin{multline*}
 U(-t)u^\eps(t)- u_-^\eps \\
 = U\left(-s\right)\left(
 u^\eps(s)- u^\eps_2(s)\right)\Big|_{s=-1/\eps^\gamma}
 +U(-s)u_2^\eps(s)\Big|_{s=-1/\eps^\gamma} - u_-^\eps 
 -\irm \lambda
 \int_{-1/\eps^\gamma}^t U(-\tau)\left(|u^\eps|^2u^\eps\right)(\tau)\dd \tau,
\end{multline*}
The first term on the right hand side is estimated thanks to
Proposition~\ref{prop:wave-op}:
\begin{equation*}
 \left\|U\left(-s\right)\left(
 u^\eps(s)- u^\eps_2(s)\right)\Big|_{s=-1/\eps^\gamma}\right\|_{\Sigma}\lesssim
 \eps^{3+\gamma}|\log\eps|.
\end{equation*}
On the other hand,
\begin{equation*}
 \left\|U(-s)u_2^\eps(s)\Big|_{s=-1/\eps^\gamma}
 - u_-^\eps\right\|_{\Sigma} =
 \|\F^{-1}\hat w_{\mid s=-1/\eps^\gamma}- u_-^\eps\|_{\Sigma}= \|\hat
 w_{\mid s=-1/\eps^\gamma}- \F( u_-^\eps)\|_{\Sigma}.
\end{equation*}
Recalling that
\begin{equation*}
 \hat w = \F( u_-^\eps) \exp \left(\irm \,|\F( u_-^\eps) |^2\log|t|\right), 
\end{equation*}
we have directly
\begin{align*}
 \|\hat w- \F( u_-^\eps)\|_{\Sigma}&\lesssim \|\left\langle x\right\rangle \left(\hat w- \F(
 u_-^\eps)\right)\|_{L^2} + \left\|\d_x \hat{u_-^\eps}\left( \exp \left(\irm \,|\F( u_-^\eps)
 |^2\log|t|\right)-1\right)\right\|_{L^2} 
\\ 
&\hspace*{60mm} 
 + |\log\eps|\left\| \F( u_-^\eps)\d_x |\F( u_-^\eps)
 |^2\right\|_{L^2}\\
 &\lesssim \eps^3|\log \eps|. 
\end{align*}
The integral term can be estimated in $\Sigma$ by
\begin{equation*}
 \int_{-1/\eps^\gamma}^t
 \|u^\eps(\tau)\|_{L^\infty}^2\|U(-\tau)u^\eps(\tau)\|_{\Sigma}\, \dd \tau.
\end{equation*}
Since $\| U(T)u^\eps(-T)\|_{\Sigma}\lesssim \eps$
from Lemmas~\ref{lem:u1} and~\ref{lem:Rw}, and
Proposition~\ref{prop:wave-op} (where $T$ 
is independent of $\eps$), we have the following uniform estimate from
\cite[Theorem~1.1]{HN98} (see also~\cite{KatoPusateri2011}), provided
that $\eps>0$ 
is sufficiently small,
\begin{equation}\label{eq:unifLinfty}
 \|u^\eps(\tau)\|_{L^\infty}\lesssim
 \frac{\eps}{\left\langle \tau\right\rangle ^{1/2}},\quad \forall \tau\in \R. 
\end{equation}
Writing
\begin{equation*}
 \|U(-\tau)u^\eps(\tau)\|_{\Sigma}\le
 \|U(-\tau)\left(u^\eps(\tau)-\u_1^\eps(\tau)\right)\|_{\Sigma} +
 \|U(-\tau)\u_1^\eps(\tau)\|_{\Sigma}, 
\end{equation*}
and using the obvious fact that
$\|U(-\tau)\u_1^\eps(\tau)\|_{\Sigma}=\eps\|v_-\|_{\Sigma}$, we have, for $t\ge
-1/\eps^\gamma$, 
\begin{equation*}
 \left\| U(-t)u^\eps(t)- u_-^\eps\right\|_{\Sigma}\lesssim
 \eps^{3}|\log \eps| + \eps^2 \int_{-1/\eps^\gamma}^t \left\|
 U(-\tau)u^\eps(\tau)- u_-^\eps\right\|_{\Sigma}
 \frac{\dd \tau}{\left\langle \tau\right\rangle } + \eps^3 \int_{-1/\eps^\gamma}^t \frac{\dd \tau}{\left\langle \tau\right\rangle }.
\end{equation*}
 Therefore, Gronwall
 lemma yields
 \begin{equation*}
 \sup_{-1/\eps^\gamma\le t\le T}\left\| U(-t)u^\eps(t)-
 u_-^\eps\right\|_{\Sigma} \lesssim 
 \eps^{3}|\log \eps| \expe^{ C \eps^2|\log \eps|+C\eps^2\log\left\langle T\right\rangle },
 \end{equation*}
 for some constant $C$ independent of $\eps$ and $T$. In particular,
 for any $\delta\in \mathopen{]}0,1]$, 
 the right hand side is $\O(\eps^{3-\delta})$ on
 $[-1/\eps^\gamma,1/\eps^\beta]$ for any $\beta>0$. This estimate
 implies the property
 \begin{equation*}
 \W (\eps v_-)=\eps v_- +\O(\eps^{3-\eta}),
 \end{equation*}
 for any $\eta>0$. In the next subsection, we improve this estimate
 by describing the first corrector term. 
\subsection{Next term in the asymptotic expansion}
\label{sec:next1}

Following e.g.~\cite{CaGa09}, a refined asymptotic expansion for
$u^\eps$, at least on bounded time intervals, is given by the solution
$\u_2^\eps$ to
\begin{equation*}
 \irm \d_t \u_2^\eps +\frac{1}{2}\d_x^2 \u_2^\eps = \lambda
 |\u_1^\eps|^2\u_1^\eps, 
\end{equation*}
where the Cauchy data must be chosen carefully.
Formally, as $\u_1^\eps$ is
 of order $\eps$, $\u_2^\eps$ is of order $\eps^3$, provided that its
 Cauchy datum is $\O(\eps^3)$.
 Set
\begin{equation*}
 \u_{\app}^\eps = \u_1^\eps+\u^\eps_2.
\end{equation*}
We now have
\begin{multline*}
 U(-t)\left(u^\eps(t)- \u^\eps_{\app}(t)\right) 
 \\
 = U\left(-s\right)\left(
 u^\eps(s)- u^\eps_2(s)+u_2^\eps(s)
 - \u_{\app}^\eps(s)\right)\Big|_{s=-1/\eps^\gamma} -\irm \lambda
 \int_{-1/\eps^\gamma}^t U(-\tau)\left(|u^\eps|^2u^\eps-
 |\u_1^\eps|^2\u_1^\eps\right)(\tau)\, \dd \tau.
\end{multline*}
We have seen before
that
\begin{equation*}
 \left\|U\left(-s\right)\left(
 u^\eps(s)- u^\eps_2(s)\right)\Big|_{s=-1/\eps^\gamma}\right\|_{\Sigma}\lesssim
 \eps^{3+\gamma}|\log\eps|=o(\eps^3). 
\end{equation*}
So if we want to catch the $\eps^3$ term, this is good. Making $U(-s)(u_2^\eps(s)
 - \u_{\app}^\eps(s))$ small ($o(\eps^3)$ in $\Sigma$) at
 $s=-1/\eps^\gamma$ is what should 
 tell us how to choose the data for $\u_2^\eps$. Indeed, the integral
 term is estimated in $L^2$ ($\Sigma$ will require more care) by
 \begin{equation*}
 \int_{-1/\eps^\gamma}^t \left(\|u^\eps(\tau)\|_{L^\infty}^2 +
 \|\u_1^\eps(\tau)\|_{L^\infty}^2
 \right)\|u^\eps(\tau)-v_1^\eps(\tau)\|_{L^2}\, \dd \tau\lesssim \eps^2
 \int_{-1/\eps^\gamma}^t
 \|u^\eps(\tau)-\u_1^\eps(\tau)\|_{L^2}\frac{\dd \tau}{\left\langle \tau\right\rangle }, 
 \end{equation*}
 and we now write
 \begin{equation*}
 \|u^\eps(\tau)-\u_1^\eps(\tau)\|_{L^2}\le \|u^\eps(\tau)-\u_{\rm
 app}^\eps(\tau)\|_{L^2}+\|\u_2^\eps(\tau)\|_{L^2}, 
 \end{equation*}
so the last term yields a smaller contribution (in terms of $\eps$ at
least) than in the previous
subsection.

Thus, everything seems to boil down to choosing $\u_2^\eps$ at time
$-1/\eps^\gamma$ in an efficient way, so that
\begin{equation*}
 \left\|U(-s)(u_2^\eps(s)-\u_{\rm
 app}^\eps(s))\Big|_{s=-1/\eps^\gamma}\right\|_\Sigma=o(\eps^3). 
\end{equation*}
By construction,
\begin{equation*}
 U(-s)(u_2^\eps(s))\Big|_{s=-1/\eps^\gamma}= \F^{-1}\left( \F(u_-^\eps)
 \exp\left(-\irm \lambda\gamma |\F(u_-^\eps)|^2\log \eps\right)\right). 
\end{equation*}
Taking into account the second term in the asymptotic expansion of the
exponential in $\hat w$, we define $\u_2^\eps$ by
\begin{equation}\label{eq:v2}
% \begin{aligned}
 U(-t) \u_2^\eps (t) %&
 = -\irm \lambda\gamma \eps^3(\log
 \eps)\F^{-1}\left( 
 |\F(v_-)|^2\F(v_-) \right) %\\
% &\quad
 -\irm \lambda\eps^3\int_{-1/\eps^\gamma}^t U(-\tau) \left(
 |U(\tau)v_-|^2U(\tau)v_-\right)\, \dd \tau. 
%\end{aligned}
\end{equation}
\begin{lemma}\label{lem:renormalisation}
 Let $v_-\in \H$. For every $t\in \R$, $U(-t)\u^\eps_2(t)/\eps^3$ converges in $L^2\cap L^\infty(\R)$, and the
 limit is independent of $\gamma$. This limit is given by
 \begin{multline*}
 U(-t) v_2(t) = -\irm \lambda\int_{-\infty}^{-1} \left( U(-\tau)\left(
 |U(\tau)v_-|^2U(\tau)v_-\right) +\frac{1}{|\tau|} \F^{-1}\left( 
 |\F(v_-)|^2\F(v_-)\right)\right) \dd \tau 
 \\
%&\quad 
-\irm \lambda\int_{-1}^t U(-\tau) \left( 
 |U(\tau)v_-|^2U(\tau)v_-\right)\dd \tau.
 \end{multline*}
 The limit also holds in $\Sigma$,
 and for all $t\in \R$,
 \begin{equation*}
 \left\|
 U(-t)\left(\frac{\u^\eps_2(t)}{\eps^3}-v_2(t)\right)\right\|_{\Sigma}\lesssim
 \eps^\gamma.
 \end{equation*}
\end{lemma}
\begin{remark}
 We emphasize that even though the limit $v_2$ is independent of
$\gamma$, the error estimate improves for large values of $\gamma$. In
view of the definition of $\u_2^\eps$, the matching condition at
$t=-1/\eps^\gamma$ 
between $u_2^\eps$ and $\u_{\app}^\eps$ is not better than
$\O\left(\eps^5(\log\eps)^2\right)$, so it is no use (in this paper) to consider
$\gamma>2$. 
\end{remark}
\begin{proof}
 The proof consists in writing a precise asymptotic expansion of the
 argument of the integral involved in~\eqref{eq:v2}. Writing, for any
 $\theta\in \mathopen{]}0,1]$,
 \begin{equation*}
 M(t) = \expe^{\irm \frac{x^2}{2t}} = 1+\O\left(\left|\frac{x^2}{t}\right|^\theta\right) ,
 \end{equation*}
 we have, in the same vein as in Lemma~\ref{lem:R},
 \begin{equation*}
 U(\tau)v_-(x) = M(\tau)D(\tau)\F v_- + R_1(\tau,x),
 \end{equation*}
 where, as $\tau\to -\infty$
 \begin{equation*}
 \|R_1(\tau)\|_{L^2} = \O\left(\frac{1}{|\tau|^{\theta}}\right),\quad
 \|R_1(\tau)\|_{L^\infty} = \O\left(\frac{1}{|\tau|^{\theta+1/2}}\right).
 \end{equation*}
 We infer, as $\tau\to -\infty$,
 \begin{equation*}
 |U(\tau)v_-|^2U(\tau)v_-=
 \expe^{\irm \frac{\pi}{4}}\frac{\expe^{\irm \frac{x^2}{2\tau}}}{|\tau|^{3/2}}
 \left(|\F(v_-)|^2\F(v_-)\right)\left(\frac{x}{\tau}\right)
 + R_2(\tau,x),
 \end{equation*}
 with
 \begin{equation*}
 \|R_2(\tau)\|_{L^2} = \O\left(\frac{1}{|\tau|^{1+\theta}}\right),\quad
 \|R_2(\tau)\|_{L^\infty} = \O\left(\frac{1}{|\tau|^{3/2+\theta}}\right).
 \end{equation*}
 In view of~\eqref{eq:U(-t)}, we next compute
 \begin{align*}
 U(-\tau)\left( |U(\tau)v_-|^2U(\tau)v_-\right) &= \irm \expe^{-\irm \frac{x^2}{2\tau}}
 \F^{-1}D\left(\frac{1}{\tau}\right) \left( \frac{ \expe^{\irm \frac{\pi}{4}}}{|\tau|^{3/2}} \left(|\F(v_-)|^2\F(v_-)\right)\left(\frac{x}{\tau}\right)\right)\\
 &\hspace*{45mm}
 +\irm \expe^{-\irm \frac{x^2}{2\tau}}
 \F^{-1}D\left(\frac{1}{\tau}\right) \left( \expe^{-\irm \frac{x^2}{2\tau}}
 R_2(\tau,x)\right)\\
 &= -\frac{\expe^{-\irm \frac{x^2}{2\tau}}}{|\tau|} \F^{-1}\left(
 |\F(v_-)|^2\F(v_-)\right)(x) + R_3(\tau,x)\\
 &= -\frac{1}{|\tau|} \F^{-1}\left(
 |\F(v_-)|^2\F(v_-)\right)(x) + R_4(\tau,x),
 \end{align*}
 where
\begin{equation*}
 \|R_4(\tau)\|_{L^2\cap L^\infty} =
 \O\left(\frac{1}{|\tau|^{1+\theta}}\right).
\end{equation*}
Set $\theta=1$,
and write
\begin{align*}
 \frac{1}{\lambda \eps^3} U(-t)\u_2^\eps(t) &= -\irm \gamma
 (\log\eps)\F^{-1} \left( |\F(v_-)|^2\F(v_-) \right)
 -\irm \int_{-1/\eps^\gamma}^{-1} U(-\tau) \left(
 |U(\tau)v_-|^2U(\tau)v_-\right)\, \dd \tau\\
&\qquad-\irm \int_{-1}^t U(-\tau) \left(
 |U(\tau)v_-|^2U(\tau)v_-\right)\dd \tau\\
 &= -\irm \gamma
 (\log\eps)\F^{-1} \left( |\F(v_-)|^2\F(v_-) \right)
% \\ &\quad
 \Mk -\Mk \irm \int_{-1/\eps^\gamma}^{-1} \left(-\frac{1}{|\tau|} \F^{-1}\left(
 |\F(v_-)|^2\F(v_-)\right) + R_4(\tau)\right)\dd \tau
 \\
 &\qquad
-\irm \int_{-1}^t U(-\tau) \left(
 |U(\tau)v_-|^2U(\tau)v_-\right)\dd \tau\\
 &= -\irm \int_{-\infty}^{-1} R_4(\tau)\dd \tau +\underbrace{ \O\left(
 \int_{-\infty}^{-1/\eps^\gamma} \frac{\dd \tau}{\tau^2}\right)}_{=\O(\eps^\gamma)}
-\irm \int_{-1}^t U(-\tau) \left(
 |U(\tau)v_-|^2U(\tau)v_-\right)\dd \tau.
\end{align*}
Therefore,
\begin{equation*}
 U(-t)v_2(t) = -\irm \lambda\int_{-\infty}^{-1} R_4(\tau)\dd \tau
 -\irm \lambda\int_{-1}^t U(-\tau) \left( 
 |U(\tau)v_-|^2U(\tau)v_-\right)\dd \tau,
\end{equation*}
and, by construction,
\begin{equation*}
 R_4(\tau) = U(-\tau)\left( |U(\tau)v_-|^2U(\tau)v_-\right) +\frac{1}{|\tau|} \F^{-1}\left(
 |\F(v_-)|^2\F(v_-)\right),
 \end{equation*}
 hence the formula of the lemma.
In particular, by the same arguments as above, the second integral
diverges logarithmically as $t\to +\infty$. We leave out the convergence in
$\Sigma$ (since $v_-\in \H$, so we can pay some 
momentum estimate when controlling $M-1$), which relies on similar ideas,
like the proof of 
Lemma~\ref{lem:R}. Note that the present argument 
\end{proof}
We can now improve the error estimate proven in the previous
subsection with only the leading order approximation $\u_1^\eps$. By
construction,
\begin{equation*}
 \left\|U(-s)(u_2^\eps(s)-\u_{\rm
 app}^\eps(s))\Big|_{s=-1/\eps^\gamma}\right\|_\Sigma=\O\left(\eps^5(\log \eps)^2\right),
\end{equation*}
and we have, if $0<\gamma<2$,
\[
 \left\|U(-t)\left(u^\eps(t)- \u^\eps_{\app}(t)\right)\right\|_{\Sigma} \lesssim
 \eps^{3+\gamma}|\log \eps|
% \\
% &\quad
 +
 \int_{-1/\eps^\gamma}^t \left\| U(-\tau)\left(|u^\eps|^2u^\eps-
 |\u_1^\eps|^2\u_1^\eps\right)(\tau)\right\|_{\Sigma}\dd \tau.
\]
In view of~\eqref{eq:est-gen},
\begin{multline*}
 \| U(-\tau)\left(|u^\eps|^2u^\eps-
 |\u_1^\eps|^2\u_1^\eps\right)(\tau)\|_{\Sigma} \lesssim
 \| u^\eps(\tau)\|_{L^\infty}^2 \| U(-\tau)\left(u^\eps-
 \u_1^\eps\right)(\tau)\|_{\Sigma} \\
 + \left( \|u^\eps(\tau)\|_{L^\infty} +
 \|\u_1^\eps(\tau)\|_{L^\infty} \right) \|U(-\tau)\u_1^\eps(\tau)\|_{\Sigma}
 \|u^\eps(\tau)-\u_1^\eps(\tau)\|_{L^\infty}. 
\end{multline*}
Recalling~\eqref{eq:unifLinfty}, we readily have
\begin{equation*}
 \|u^\eps(\tau)\|_{L^\infty} \lesssim \frac{\eps}{\left\langle \tau\right\rangle ^{1/2}}.
\end{equation*}
On the other hand,
\begin{equation*}
 \|U(-\tau)\u_1^\eps(\tau)\|_{\Sigma}=\eps\|v_-\|_{\Sigma},
\end{equation*}
and Sobolev embedding yields, along with~\eqref{eq:GNlibre},
\begin{equation*}
 \|\u_1^\eps(\tau)\|_{L^\infty} \lesssim \frac{\eps}{\left\langle \tau\right\rangle ^{1/2}},
\end{equation*}
as well as
\begin{equation*}
 \|u^\eps(\tau)-\u_1^\eps(\tau)\|_{L^\infty}\lesssim
 \frac{1}{\left\langle \tau\right\rangle ^{1/2}}\|U(-\tau)\left(u^\eps-
 v_1^\eps\right)(\tau)\|_{\Sigma} . 
\end{equation*}
Therefore,
\begin{align*}
 \| U(-\tau)\left(|u^\eps|^2u^\eps-
 |\u_1^\eps|^2\u_1^\eps\right)(\tau)\|_{\Sigma} &\lesssim
 \frac{\eps^2}{\left\langle \tau\right\rangle } \|U(-\tau)\left(u^\eps-
 \u_1^\eps\right)(\tau)\|_{\Sigma}\\
 &\lesssim
 \frac{\eps^2}{\left\langle \tau\right\rangle } \left(\big\|U(-\tau)(u^\eps-
 \u_{\app}^\eps)(\tau)\big\|_{\Sigma} +
 \|U(-\tau)\u_2^\eps(\tau)\|_{\Sigma}\right). 
\end{align*}
We see that
there exists $C$ such that for all $T\ge 1$,
\begin{equation*}
 \|U(-\tau)\u_2^\eps(\tau)\|_{\Sigma}\lesssim \eps^3\log\left\langle T\right\rangle ,\quad
 \forall \tau \in\left[-\frac{1}{\eps^\gamma},T\right],
 \end{equation*}
 where the logarithmic correction is necessarily present for large
 positive $\tau$, as pointed out above. Gronwall lemma now yields, for
 all $0<\gamma<2$, 
 \begin{equation*}
 \sup_{-1/\eps^\gamma\le t\le T}\left\| U(-t)\left(u^\eps(t)-
 \u_{\app}^\eps(t)\right)\right\|_{\Sigma} \lesssim 
 \eps^{3+\gamma}|\log \eps| \expe^{ C \eps^2|\log \eps|+C\eps^2\log\left\langle T\right\rangle },
 \end{equation*}
 for some constant $C$ independent of $\eps$ and $T$. In particular,
 the right hand side is $o(\eps^{3})$ on
 $[-1/\eps^\gamma,1/\eps^\beta]$ for any $\beta>0$. 

 The first point of Theorem~\ref{theo:main} follows, by considering
 the above error estimate at time $t=0$ and setting $w_2
 = v_{2\mid t=0}$. 



\section{Asymptotic completeness}
\label{sec:complete}

\subsection{Main steps of the construction of \texorpdfstring{$u_+$}{u+}}

We recall the main steps from the proof of~\cite[Theorem~1.2]{HN98}, a
result which we state in the particular case that we shall consider
(the initial data may have a different regularity there, and an extra
short range nonlinearity can be incorporated). In passing, we keep
track of the size of the remainders more precisely.
\begin{theorem}[From~\cite{HN98}]\label{theo:complete}
 Let $u_0\in \Sigma$, with $\|u_0\|_{\Sigma}=\eps'\le \delta$, where
 $\delta>0$ is sufficiently small. There exist unique functions $W\in
 L^\infty\cap L^2$ and $\Phi\in L^\infty$ such that, for $t\ge 1$ and $C\delta<\alpha<1/4$, 
 \begin{equation}
 \label{eq:1.4}
 \left\| \F\left(U(-t)u(t)\right) \exp\left(\irm \lambda\int_1^t|\hat
 u(\tau)|^2\frac{\dd \tau}{\tau}\right)-W\right\|_{L^2\cap L^\infty} 
 \le C(\eps')^3 t^{-\alpha+C(\eps')^2}, 
 \end{equation}
and
\begin{equation}
 \label{eq:1.5}
 \left\| \lambda \int_1^t |\hat u(\tau)|^2\frac{\dd \tau}{\tau}
 -\lambda |W|^2 \log t-\Phi\right\|_{L^\infty} 
 \le C(\eps')^3 t^{-\alpha+C(\eps')^2}.
 \end{equation}
 In particular,
 \begin{equation*}
 u(t,x) =\frac{1}{(\irm t)^{1/2}}W\left(\frac{x}{t}\right)\exp
 \left(\irm \frac{x^2}{2t}-\irm \lambda \left|W\left(\frac{x}{t}\right)\right|^2\log t -
 \irm \Phi\left(\frac{x}{t}\right)\right) + \rho(t,x),
 \end{equation*}
 with
 \begin{equation*}
 \|\rho(t)\|_{L^2} \le C\eps' t^{-\alpha +C(\eps')^2},\quad
 \|\rho(t)\|_{L^\infty} \le C\eps' t^{-1/2-\alpha +C(\eps')^2} .
 \end{equation*}
\end{theorem}
In~\cite{HN98}, the time decay in~\eqref{eq:1.5} is raised to the
power $2/3$, because the setting is more general and includes the case
of dimension three. We explain below why this power can be discarded
in the above statement. The asymptotic state $u_+$ is naturally given
by
\begin{equation*}
 u_+ = \F^{-1} \left(W\expe^{-\irm \Phi}\right).
\end{equation*}
\begin{proof}[Sketch of the proof]
 As announced above, we recall the main steps from~\cite{HN98}, and
 slightly improve some estimates in terms of the powers of
 $\eps'$.

As evoked before in the case of~\eqref{eq:unifLinfty}, (the proof of)~\cite[Theorem~1.1]{HN98} provides global estimates
for $u$ (proven in~\cite[Lemmas~3.2 and~3.3]{HN98}, see also~\cite{KatoPusateri2011}):
\begin{equation}\label{eq:unif-u}
 \|u(t)\|_{L^\infty}\lesssim
 \frac{\eps'}{\left\langle t\right\rangle ^{1/2}},\quad
 \|U(-t)u(t)\|_{\Sigma}\lesssim \eps' \left\langle t\right\rangle ^{C(\eps')^2}.
\end{equation}
Denote
 \begin{equation*}
 v(t)=U(-t)u(t),\quad \hat w = \hat v \underbrace{\exp\left(\irm \lambda
 \int_1^t |\hat v(\tau)|^2\frac{\dd \tau}{\tau}\right)}_{\eqqcolon B(t)}.
 \end{equation*}
 Then~\eqref{eq:NLS} is recasted as
 \begin{equation}\label{eq:w-complete}
 \irm \d_t \hat w =\frac{\lambda}{t}B(t)\left(I_1+I_2\right),
 \end{equation}
 where
 \begin{equation*}
 I_1(t)= \F\left(M(-t)-1\right)\F^{-1}\left(|\hat{M(t) v}|^2\hat{M(t) v}\right),\quad
 I_2(t)= \left|\hat{M(t)v}\right|^2 \hat{M(t)v}- |\hat
 v|^2\hat v.
\end{equation*}
These source terms are controlled as follows, for $t\ge 1$:
\begin{equation*}
 \|I_1(t)\|_{L^2\cap L^\infty} + \|I_2(t)\|_{L^2\cap L^\infty} \lesssim \frac{\|v\|_{\Sigma}^3}{t^\alpha},
\end{equation*}
provided that $\alpha<1/2$. In view of~\eqref{eq:unif-u} and the
definition of $v$, this entails
\begin{equation*}
 \|I_1(t)\|_{L^2\cap L^\infty} + \|I_2(t)\|_{L^2\cap L^\infty}
 \lesssim (\eps')^3 t^{-\alpha+3C(\eps')^2}. 
 \end{equation*}
Renaming $3C$ to $C$, and integrating~\eqref{eq:w-complete}, we
readily find that there exists 
$W\in L^2\cap L^\infty$ such that
\begin{equation}\label{eq:cvw}
 \|\hat w(t)-W\|_{L^2\cap L^\infty} \lesssim (\eps')^3 t^{-\alpha+C(\eps')^2},
\end{equation}
thus proving~\eqref{eq:1.4}, since
\begin{equation*}
 |\hat v(t,\xi) | = \left|\F\left(U(-t)u(t)\right)(\xi)\right|=
 \left|\expe^{-\irm t\frac{\xi^2}{2}} \hat u(t,\xi)\right|=|\hat
 u(t,\xi)|.
\end{equation*}
The phase corrector $\Phi$ is obtained by introducing the function
\begin{equation}\label{eq:defPsi}
 \Psi (t) =\lambda\int_1^t \left(|\hat w(\tau)|^2-|\hat
 w(t)|^2\right)\frac{\dd \tau}{\tau}. 
\end{equation}
We have, for $t>s>2$,
\begin{equation*}
 \Psi(t)-\Psi(s) = \lambda\int_s^t \left(|\hat w(\tau)|^2-|\hat
 w(t)|^2\right)\frac{\dd \tau}{\tau}+\lambda \left( |\hat w(t)|^2-|\hat
 w(s)|^2\right)\log s.
\end{equation*}
In view of the above estimates, for $t_2>t_1>2$,
\begin{align*}
 \left| |\hat w(t_2)|^2-|\hat
 w(t_1)|^2\right|&\lesssim \left( |\hat w(t_2)|+|\hat
 w(t_1)|\right) |\hat w(t_2)-\hat
 w(t_1)|\\
 &\lesssim \left( |\hat w(t_2)|+|\hat
 w(t_1)|\right)(\eps')^3 t_1^{-\alpha+C(\eps')^2}.
\end{align*}
On the other hand, we have
\begin{align*}
 |\hat w(t)|\lesssim \eps',
\end{align*}
from~\eqref{eq:cvw}. Therefore,
\begin{equation}\label{eq:diff-hat-w}
 \left| |\hat w(t_2)|^2-|\hat
 w(t_1)|^2\right|
 \lesssim (\eps')^4 t_1^{-\alpha+C(\eps')^2},\quad t_2>t_1>2,
\end{equation}
and so
\begin{equation*}
 |\Psi(t)-\Psi(s) |\lesssim (\eps')^4 s^{-\alpha+C(\eps')^2}\log s.
\end{equation*}
In particular, there exists $\Phi\in L^\infty$ such that, for $t>2$,
\begin{equation}\label{eq:cvPhi}
 |\Phi-\Psi(t) |\lesssim (\eps')^4 t^{-\alpha+C(\eps')^2}\log t.
\end{equation}
The estimate~\eqref{eq:1.5} then follows from the identity
\begin{equation*}
 \lambda \int_1^t |\hat w(\tau)|^2\frac{\dd \tau}{\tau} =
 \lambda|W|^2\log t + \Phi + \Psi(t)-\Phi+ \lambda\left( |\hat
 w(t)|^2-|W|^2\right)\log t,
\end{equation*}
using~\eqref{eq:cvw} and~\eqref{eq:cvPhi}, possibly modifying the
constant $C$, and using the fact that the condition on $\alpha$ in
Theorem~\ref{theo:complete} is open. 

In view of~\eqref{eq:1.4} and~\eqref{eq:1.5}, we find
\begin{equation*}
 \F\left(U(-t)u(t)\right) = W \expe^{-\irm \lambda |W|^2\log t -\irm \Phi} +r(t),
\end{equation*}
where
\begin{equation*}
 \|r(t)\|_{L^2\cap
 L^\infty} \lesssim (\eps')^3 t^{-\alpha +C(\eps')^2}. 
\end{equation*}
Writing 
\begin{equation*}
 1=U(t)U(-t)= M(t)D(t)\F M(t) U(-t) = M(t)D(t)\F U(-t) + \cR(t) \F
 U(-t) ,
\end{equation*}
we infer
\begin{equation*}
 u(t,x) =M(t)D(t)\F U(-t) u(t)+\rho_1(t,x),
\end{equation*}
with, in view of Lemma~\ref{lem:R},
\begin{equation*}
 \|\rho_1(t) \|_{L^2} \lesssim \frac{\eps'}{\sqrt t}\|\F
 U(-t) u(t)\|_{L^2} \lesssim \frac{\eps'}{\sqrt t},
\end{equation*}
and, since $\alpha<1/4$, using the first point of Lemma~\ref{lem:R}
with $\theta=\alpha$ and $s=1-2\alpha$, 
\begin{equation*} 
\|\rho_1(t) \|_{L^\infty} \lesssim \frac{\eps'}{t^{1/2+\alpha}}\|\F
 U(-t) u(t)\|_{H^1}\lesssim \eps' t^{-1/2-\alpha+C(\eps')^2}.
\end{equation*}
We conclude thanks to the estimates
\begin{align*}
 \|M(t)D(t)r(t)\|_{L^2}&=\|r(t)\|_{L^2} \lesssim (\eps')^3 t^{-\alpha
 +C(\eps')^2},\\
\|M(t)D(t)r(t)\|_{L^\infty}&=\frac{1}{\sqrt
 t}\|r(t)\|_{L^\infty}\lesssim (\eps')^3 t^{-1/2-\alpha
 +C(\eps')^2},
\end{align*}
since $u(t)= M(t)D(t) \big( W \expe^{-\irm \lambda |W|^2\log t -\irm \Phi}\big) +
M(t)D(t)r(t)+\rho_1(t)$. 
\end{proof}

\subsection{Modified scattering operator: leading order
 behavior near the origin}\label{sec:complete1}

In view of Theorem~\ref{theo:complete}, the asymptotic state that we
consider is $u_+^\eps= \F^{-1}\big(W^\eps \expe^{-\irm \Phi^\eps}\big)$, with
$u_0^\eps = u^\eps_{\mid t=0}$. The main remark at this stage is that
the reduction presented in the proof of Theorem~\ref{theo:complete}
boils down the analysis of the asymptotic behavior of $u_+^\eps$ as
$\eps\to 0$ to a regular asymptotic expansion. In order to treat the
last two cases of Theorem~\ref{theo:main}, we suppose that, for $0<\eta<2$,
\begin{equation*}
 u_0^\eps = \eps v_0 + \eps^3 w_2 + \O\left(\eps^{5-\eta} \right),
\end{equation*}
with $v_0,w_2\in \Sigma$ and the remainder term is (measured) in
$\Sigma$. For the last point of Theorem~\ref{theo:main}, we assume 
$w_2=\rho^\eps=0$, while for the second point, $v_0=v_-$ and $w_2 =
v_{2\mid t=0}$ as in the first point.

We keep the same notations as the proof of
Theorem~\ref{theo:complete}, with now $\eps$ instead of $\eps'$. The
definition of $\hat w$, in particular the term $B$, requires as a
first step the asymptotic description of $u^\eps$ at time $t=1$ instead
of only $t=0$. 





In view of~\cite[Theorem~1.1]{HN98}, as in~\eqref{eq:unif-u}, we have
\begin{equation*}
 \|u^\eps(t)\|_{L^\infty}\lesssim
 \frac{\eps}{\left\langle t\right\rangle ^{1/2}},\quad
 \|U(-t)u^\eps(t)\|_{\Sigma}\lesssim \eps \left\langle t\right\rangle ^{C\eps^2}.
 \end{equation*}
 In view of the proof of~\eqref{eq:1.4} and~\eqref{eq:1.5}, we
 directly know that
 \begin{equation*}
 \|\d_t \hat w\|_{L^1([1,\infty[,L^2\cap L^\infty)}\lesssim \eps^3,
 \end{equation*}
 and so
 \begin{equation*}
 \hat w(t) = \hat w(1) + \O(\eps^3) \text{ in }L^\infty([1,\infty[,L^2\cap L^\infty).
 \end{equation*}
 By definition, $\hat w(1) = \hat v(1) = \F\left( U(-1)u^\eps(1)\right)= \eps
 \hat v_0+\O(\eps^3)$ in $\Sigma$.

 Regarding $\Phi$, the definition~\eqref{eq:defPsi}, and~\eqref{eq:diff-hat-w},
 yield
 \begin{equation*}
 \|\Phi\|_{L^\infty}\le
 \|\Phi-\Psi(t)\|_{L^\infty}+\|\Psi(t)\|_{L^\infty}\lesssim \eps^4 ,
 \end{equation*}
where the right hand side is estimated uniformly in $t\ge 1$. 
 Therefore,
 \begin{equation*}
 u_+^\eps = \F^{-1}\left(W\expe^{-\irm \Phi}\right) = \eps v_0+\O(\eps^3)\text{ in
 }L^2. 
 \end{equation*}
 Therefore, at leading order, we have
 \begin{equation*}
 \Scatt(\eps v_-) = \eps v_-+\O(\eps^3),\quad \left(W_+^{\rm mod}\right)^{-1}
 (\eps v_0) = \eps v_ 0 +\O(\eps^3)\quad \text{in }L^2.
 \end{equation*}






\section{Higher order asymptotic expansion of the final state}
\label{sec:complete-higher}


The higher order asymptotic expansion for $u_+^\eps$, involving an
 $\eps^3$ term, requires more work, even at the formal level. We first
 show how to derive this term, then, in a final subsection, prove the
 error estimates announced in Theorem~\ref{theo:main}.
 
 \subsection{Derivation of the first corrector}

 
We first examine the value of $u^\eps$ at time $t=1$, in view of
future connections with 
Theorem~\ref{theo:complete}.
The idea is the same as in Section~\ref{sec:next1}: for finite time,
since we consider small data, the asymptotic expansion is given by
Picard's scheme, and we have directly from Section~\ref{sec:wave-zero},
\begin{equation*}
 u^\eps_{\mid t=1} = \eps U(1)v_0 + \eps^3 \tilde w_2 + \eps^{5-\eta} \hat\rho^\eps,
\end{equation*}
with
\begin{equation*}
 \tilde w_2 = U(1)w_2 -\irm \lambda\int_0^1 U(1-s)\left(|U(s)v_0|^2U(s)v_0\right)\dd s,
\end{equation*}
and $\|\rho^\eps\|_{\Sigma}\lesssim 1$ as $\eps \to 0$. 
\smallbreak

To derive the first corrector, we plug the asymptotic expansion 
\begin{equation*}
 \hat w(t) = \eps \hat v_0 + \eps^3 \hat\nu_2 (t)
 +\eps^{5-\eta}\hat{r^\eps}(t)
\end{equation*}
into~\eqref{eq:w-complete}, with $r^\eps=\O(1)$ in some topology we
shall precise. We first proceed formally, and then check
that the above ansatz indeed provides a corrector of $u^\eps_+$ in
$L^2$. We compute successively
\begin{equation*}
 |\hat v|^2=|\hat w|^2 = \eps^2|\hat v_0|^2 +\O(\eps^4),
\end{equation*}
hence, for bounded $t$,
\begin{equation*}
 B(t) = \expe^{\irm\lambda \int_1^t |\hat v(\tau)|^2\frac{\dd \tau}{\tau}} = 1+
 \irm \lambda \eps^2|\hat v_0|^2\log t + \O(\eps^4),
\end{equation*}
and
\begin{equation*}
 \hat v = \hat w \bar B = \eps \hat v_0 +\eps^3\left( \hat\nu_2
 -\irm \lambda |\hat v_0 |^2\hat v_0 \log t\right)+\O(\eps^{5-\eta}).
\end{equation*}
We then expand the factor $I_1^\eps$ and $I_2^\eps$ from~\eqref{eq:w-complete}:
\begin{align*}
 I_1^\eps(t)& = \eps^3
 \F\left(M(-t)-1\right)\F^{-1}\left(\left|\hat{M(t)v_0}\right|^2 
 \hat{M(t)v_0} \right)+\O(\eps^5),\\
 I_2^\eps(t) & = \eps^3\left( \left|\hat{M(t)v_0}\right|^2
 \hat{M(t)v_0} - |\hat v_0 |^2\hat v_0 \right)+\O(\eps^5).
\end{align*}
Therefore, the natural candidate for $\nu_2$ satisfies:
\begin{equation*}
 \d_t \hat\nu_2 = \frac{\lambda}{t}(J_1+J_2),
\end{equation*}
where
\begin{align*}
 J_1& = \F\left(M(-t)-1\right)\F^{-1}\left(\left|\hat{M(t)v_0}\right|^2
 \hat{M(t)v_0} \right),\\
 J_2 & = \left|\hat{M(t)v_0}\right|^2
 \hat{M(t)v_0} - |\hat v_0 |^2\hat v_0 .
\end{align*}
The Cauchy datum for $\nu_2$ is given at $t=1$: on the one hand, the
above asymptotic expansion implies
\begin{equation*}
 \hat v_{\mid t=1} =\eps \hat v_0 +\eps^3\hat{\nu}_{2\mid t=1}
 +\O(\eps^{5-\eta}),
\end{equation*}
while on the other hand, the analysis of Section~\ref{sec:wave-zero} yields
\begin{equation*}
 U(-t)u^\eps(t)_{\mid t=1} = \eps v_0 + \eps^3 U(-1)\tilde w_2 + \O(\eps^{5-\eta}),
\end{equation*}
so we infer
\begin{equation*}
 \nu_{2\mid t=1}= U(-1)\tilde w_2 = w_2
 -\irm \lambda \int_0^1 U(-s)\left( |U(s)v_0|^2U(s)v_0\right)\dd s.
\end{equation*}
As $v_0\in \Sigma$, the proof of Theorem~\ref{theo:complete} readily
yields $\frac{1}{t}(J_1+J_2)\in L^1([1,\infty\mathclose{[}\,;L^2\cap L^\infty)$.
Since $\Sigma$ is a Banach algebra, invariant under the Fourier
transform, and such that $\Sigma\hookrightarrow L^1\cap L^\infty$, we
have $\nu_{2\mid t=1}\in \Sigma$, and so $\nu_2\in
L^\infty([1,\infty[\,;L^2\cap L^\infty)$ and
\begin{equation*}
 \hat\nu_2(t)\Tend t \infty \hat\nu_{2\mid t=1} + \lambda \int_1^\infty
 \big(J_1(\tau)+J_2(\tau)\big)\frac{\dd \tau}{\tau} \quad\text{in }L^2\cap
 L^\infty. 
\end{equation*}

However, in the asymptotic expansion of $B$, $I_1$ and $I_2$,
the more precise information that we need is that the remainder term
is $\O(\eps^5)$ together with some algebraic decay in time, so we
obtain an asymptotic expansion for $\hat w$ in $L^2$. This requires
more effort, and the proof is presented in the next subsection.


\subsection{Justification of the second order asymptotic expansion}
As a first step, we analyze the regularity of the corrector
$\nu_2$. To do so, we introduce an extra function space: for
$0<\gamma<1$, set
 \begin{equation*}
 \Sigma^\gamma = H^\gamma\cap \F(H^\gamma)=\{f\in H^\gamma(\R),\
 \|f\|_{\Sigma^\gamma}\coloneqq \|f\|_{H^\gamma} +\| \left\langle x\right\rangle ^\gamma f\|_{L^2}<\infty\}.
 \end{equation*}
 We note that like $\Sigma$, $\Sigma^\gamma$ is invariant under the
 Fourier transform, as well as the action of $U(t)$, and, when
 $\gamma>1/2$, it is an
 algebra embedded into $L^1\cap L^\infty$.
\begin{lemma}\label{lem:nu2}
 Let $v_0,w_2\in \Sigma$. Recall that the first corrector is defined by
 \begin{equation*}
 \hat\nu_2(t) = \hat w_2
 -\irm \lambda \int_0^1 \F U(-s)\left( |U(s)v_0|^2U(s)v_0\right)\dd s +\lambda\int_1^t
 (J_1(\tau)+J_2(\tau))\frac{\dd \tau}{\tau} , 
 \end{equation*}
 where $J_1$ and $J_2$ are defined by
 \begin{align*}
 J_1(t)& = \F\left(M(-t)-1\right)\F^{-1}\left(\left|\hat{M(t)v_0}\right|^2
 \hat{M(t)v_0} \right),\\
 J_2(t) & = \left|\hat{M(t)v_0}\right|^2
 \hat{M(t)v_0} - |\hat v_0 |^2\hat v_0 .
 \end{align*}
 Then $J_1,J_2\in L^\infty([1,\infty[,\Sigma)$, with 
 \begin{equation*}
 \sup_{t\ge 1}\|J_1(t)\|_\Sigma + \sup_{t\ge
 1}\|J_2(t)\|_\Sigma\lesssim \|v_0\|_{\Sigma}^3,
 \end{equation*}
 and, for any $1/2<\gamma<1$,
 \begin{equation*}
 \|J_1(t)\|_{\Sigma^\gamma}+ \|J_2(t)\|_{\Sigma^\gamma}\lesssim
 \frac{\|v_0\|_{\Sigma}^3}{t^{(1-\gamma)/2}} ,
 \end{equation*}
 and therefore $\nu_2\in
 L^\infty_{\rm loc}([1,\infty[,\Sigma)\cap
 L^\infty([1,\infty[,\Sigma^\gamma)$. Finally, $\hat \nu_2(t)\to
 \hat\nu_2^\infty$ in $\Sigma^\gamma$ as $t\to \infty$, where
 \begin{equation*}
 \hat\nu_2^\infty = \hat w_2
 -\irm \lambda \int_0^1 \F U(-s)\left( |U(s)v_0|^2U(s)v_0\right)\dd s +\lambda\int_1^\infty
 \big(J_1(\tau)+J_2(\tau)\big)\frac{\dd \tau}{\tau} .
 \end{equation*}
\end{lemma}
\begin{proof}[Proof of Lemma~\ref{lem:nu2}]
 Since $w_2,v_0\in \Sigma$, the first two terms defining $\hat
 \nu_2$ belong to $\Sigma$ (as it is an algebra). For any (fixed)
 $t$, $J_1,J_2\in \Sigma$, and $\nu_2(t)\in \Sigma$: the uniform
 bound in time is straightforward, but to get some time decay, we pay
 a little regularity. Let $0<\gamma<1$: 
 like in the proof of Lemma~\ref{lem:R}, write, for $t\ge 1$,
 \begin{align*}
 \|J_1(t)\|_{\Sigma^\gamma} &= \left\|
 \left(M(-t)-1\right)\F^{-1}\left(\left|\hat{M(t)v_0}\right|^2 
 \hat{M(t)v_0} \right)\right\|_{\Sigma^\gamma}\\
 &\lesssim
 \frac{1}{t^{(1-\gamma)/2}}
 \left\|
 |x|^{1-\gamma}\F^{-1}\left(\left|\hat{M(t)v_0}\right|^2 
 \hat{M(t)v_0} \right)\right\|_{\Sigma^\gamma} \\
 &\lesssim \frac{1}{t^{(1-\gamma)/2}} \left\|
 \F^{-1}\left(\left|\hat{M(t)v_0}\right|^2 
 \hat{M(t)v_0} \right)\right\|_{\Sigma}\lesssim
 \frac{\|v_0\|_{\Sigma}^3}{t^{(1-\gamma)/2}} , 
 \end{align*}
 since $\Sigma$ is an algebra. The assumption
 $\gamma>1/2$ simplifies the computations in the case of $J_2$, as
 $\Sigma^\gamma$ is an algebra, and the first 
 inequality below is straightforward: 
 \begin{align*}
 \|J_2(t)\|_{\Sigma^\gamma} &\lesssim \left(\|M(t)v_0\|_{\Sigma^\gamma}^2 +
 \|v_0\|_{\Sigma^\gamma}^2\right)
 \|(M(t)-1)v_0\|_{\Sigma^\gamma}\\
 &\lesssim
 \|v_0\|_{\Sigma}^2\frac{1}{t^{(1-\gamma)/2}}\||x|^{1-\gamma}
 v_0\|_{\Sigma^\gamma} \lesssim
 \frac{\|v_0\|_{\Sigma}^3}{t^{(1-\gamma)/2}} . 
 \end{align*}
 The lemma follows, since the extra decay in time,
 $t^{(\gamma-1)/2}$, ensures the convergence of the integral in the
 last term defining $\hat \nu_2$. 
\end{proof}
We recall two uniform estimates which will be of constant use in the
course of the proof:
\begin{equation}\label{eq:control-hat-v}
 \|\hat v(t)\|_{L^\infty}+
 \left\|\hat{M(t)v}\right\|_{L^\infty}\lesssim \eps.
\end{equation}
The first quantity is estimated in the proof of
Theorem~\ref{theo:complete}, thanks to~\eqref{eq:cvw}. The second one
is controlled thanks to~\eqref{eq:unif-u}, since, using
\eqref{eq:U(-t)}, 
\begin{equation*}
 \left\|\hat{M(t)v}\right\|_{L^\infty}= \left\|\F
 M(t)U(-t)u(t)\right\|_{L^\infty}
 =\left\|D\left(\frac{1}{t}\right)M(-t)u(t)\right\|_{L^\infty}
= \sqrt t \left\|u(t)\right\|_{L^\infty}.
\end{equation*}

The justification of the formal asymptotic expansion derived in the
previous subsection relies on a Gronwall type argument. The error that
we want to control is $\hat w -\eps\hat v_0 -\eps^3\hat
\nu_2$, which solves, by construction (keeping in mind that $v_0$ does
not depend on time),
\begin{equation*}
 \d_t\left( \hat w -\eps\hat v_0 -\eps^3\hat \nu_2\right)=
 \frac{\lambda}{t}B(I_1^\eps+I_2^\eps) -\frac{\lambda}{t}\eps^3
 (J_1+J_2). 
\end{equation*}
We rewrite the right hand side as
\begin{equation*}
 \frac{\lambda}{t}B\left(I_1^\eps+I_2^\eps-\eps^3J_1-\eps^3J_2\right)
 +\frac{\lambda}{t}(B-1)\eps^3 
 (J_1+J_2), 
\end{equation*}
and the last term will be considered as a source term. Indeed, from
the definition
\begin{equation*}
 B(t) = \expe^{\irm \lambda \int_1^t |\hat v(\tau)|^2\frac{\dd \tau}{\tau}}=
 \expe^{\irm\lambda \int_1^t |\hat w(\tau)|^2\frac{\dd \tau}{\tau}} ,
\end{equation*} 
so from~\eqref{eq:control-hat-v}, 
\begin{equation}
 \label{eq:B-1}
 \|B(t) -1\|_{L^\infty}\lesssim \eps^2\log t,
\end{equation}
and Lemma~\ref{lem:nu2} yields
\begin{equation}\label{eq:source-nu2-1}
 \left\|\eps^3(B-1)
 (J_1+J_2)\right\|_{L^2}\lesssim \eps^5 
\frac{\log t }{t^{(1-\gamma)/2}},
\end{equation}
for any $1/2<\gamma<1$.

We now focus on the term $B(I_1^\eps+I_2^\eps-\eps^3J_1-\eps^3J_2)$,
to estimate it in $L^2$,
and examine successively $B(I_1^\eps-\eps^3J_1)$
and $B(I_2^\eps-\eps^3J_2)$. First, by definition,
\begin{equation*}
 I_1^\eps-\eps^3J_1 = \F\left(M(-t)-1\right)\F^{-1}\left(
\left|\hat{M(t)v}\right|^2
 \hat{M(t)v}
 -\eps^3\left|\hat{M(t)v_0}\right|^2
 \hat{M(t)v_0} \right),
\end{equation*}
so the term $v-\eps v_0$ is naturally factored out, and we must relate
it to the left hand side, involving $w-\eps v_0-\eps^3\nu_2$. Since
$B\bar B=1$, we can write
\begin{equation*}
 \hat \rho \coloneqq \hat w -\eps \hat v_0-\eps^3\hat \nu_2 = B\hat v
 -\eps \hat v_0-\eps^3\hat \nu_2 = B\left(\hat v
 -\eps \bar B \hat v_0-\eps^3\bar B\hat \nu_2\right),
\end{equation*}
so we have
\begin{equation}
 \label{eq:recoll-DA}
 \hat v-\eps\hat v_0 = \bar B \left( \hat w -\eps
 \hat v_0-\eps^3\hat \nu_2\right) + \eps (\bar B -1)\hat
 v_0 +\eps^3 \bar B \hat \nu_2.
\end{equation}
The first term on the right hand side will be treated differently from
the last two terms, which will be considered as source terms. Using
the formula
\begin{equation*}
 |z_2|^2z_2 - |z_1|^2z_1 = |z_2|^2(z_2-z_1) +z_1\RE
 (z_2-z_1)(\bar z_2+\bar z_1),
\end{equation*}
with $z_2 = \F(Mv)$ and $z_1=\eps\F(Mv_0)$, as well as
\eqref{eq:recoll-DA}, yielding
\begin{equation}\label{eq:redress-DA}
 \F\left(M(v-\eps v_0)\right) = \F M \F^{-1}\left( \bar B \left( \hat w -\eps
 \hat v_0-\eps^3\hat \nu_2\right) + \eps (\bar B -1)\hat
 v_0 +\eps^3 \bar B \hat \nu_2\right),
\end{equation}
we write
\begin{equation*}
 \|I_1^\eps-\eps^3J_1\|_{L^2}\le \|G_1\|_{L^2} +
 \|S_1\|_{L^2}+\|S_2\|_{L^2}, 
\end{equation*}
where, simplifying the algebraic structure for the sake of presentation,
\begin{align*}
 G_1& = \F\left(M^{-1}-1\right)\F^{-1}\left( \left(|\F(M v)|^2 + \eps^2|\F(Mv_0)|^2\right) \F
 M \F^{-1} \left( \bar B \hat \rho \right)\right),\\
 S_1& = \eps \F\left(M^{-1}-1\right)\F^{-1}\left( \left(|\F(M v)|^2 +
 \eps^2|\F(Mv_0)|^2\right) \F 
 M \F^{-1} \left( (\bar B-1) \hat v_0 \right)\right),\\
 S_2 & = \eps^3 \F\left(M^{-1}-1\right)\F^{-1}\left( \left(|\F(M v)|^2 +
 \eps^2|\F(Mv_0)|^2\right) \F 
 M \F^{-1} \left( \bar B \hat \nu_2 \right)\right).
\end{align*}
We have left out the time variable in the expression
 of $M$ in order to lighten the notation, and do so below for the
 same reason.
The Gronwall term $G_1$ is estimated by merely using the boundedness
on $L^2$ of $ \F\left(M^{-1}-1\right)\F^{-1}$, and the unitarity of $\F M
\F^{-1}$, and we write
\begin{equation*}
 \|G_1\|_{L^2} \lesssim \left(\|\F(Mv)\|_{L^\infty}^2 +\eps^2
 \|\F(Mv_0)\|_{L^\infty}^2 \right) \|\hat \rho\|_{L^2}\lesssim \eps^2
 \|\hat \rho\|_{L^2}, 
\end{equation*}
by using~\eqref{eq:control-hat-v}. For the source term $S_1$ and
$S_2$, we invoke the same arguments as in the proof of
Lemma~\ref{lem:R}, and the fact that $\F 
 M \F^{-1}$ is unitary on $H^s$ for any $s\ge 0$:
\begin{align*}
 \|S_1\|_{L^2}&\lesssim \frac{\eps}{\sqrt t} \left\| \left(|\F(M v)|^2 +
 \eps^2|\F(Mv_0)|^2\right) \F 
 M \F^{-1} \left( (\bar B-1) \hat v_0 \right)\right\|_{H^1}\\
& \lesssim \frac{\eps}{\sqrt t}\left(\|v\|_{\F (H^1)}^2 + \eps^2
 \|v_0\|_{\F (H^1)}^2 \right) \|(\bar B-1)\hat v_0\|_{H^1}\\
& \lesssim \frac{\eps^3}{\sqrt t}t^{2C\eps^2} \|(\bar B-1)\hat
 v_0\|_{H^1}, 
\end{align*}
where we have used~\eqref{eq:unif-u} for the last inequality. The last
term is controlled by
\begin{equation*}
 \|(\bar B-1)\hat
 v_0\|_{H^1}\le \|\bar B-1\|_{L^\infty} \|v_0\|_\Sigma + \|\hat
 v_0\d_x \bar B \|_{L^2}\lesssim \eps^2\log t + \|\hat
 v_0\d_x \bar B \|_{L^2}, 
\end{equation*}
in view of~\eqref{eq:B-1}. Writing
\begin{equation*}
 \hat
 v_0\d_x \bar B = -2\irm \lambda\hat v_0 \int_1^t
 \RE\left(\overline{\hat v}\d_x \hat v\right)\frac{\dd \tau}{\tau},
\end{equation*}
\eqref{eq:unif-u} now yields
\begin{equation*}
 \|\hat
 v_0\d_x \bar B \|_{L^2}\lesssim \|v_0\|_{L^2}\int_1^t \|
 \hat v(\tau)\|_{L^\infty} \|\d_x \hat
 v(\tau)\|_{L^2}\frac{\dd \tau}{\tau}\lesssim \eps^2 \int_1^t
 \frac{\dd \tau}{\tau^{1-C\eps^2}} \lesssim \eps^2 t^{C\eps^2}.
\end{equation*}
The estimate for $S_2$ is rather similar: 
\begin{align*}
 \|S_2\|_{L^2} & \lesssim \frac{\eps^3}{t^{1/2}}\left\| \left(|\F(M v)|^2 +
 \eps^2|\F(Mv_0)|^2\right) \F 
 M \F^{-1} \left( \bar B \hat \nu_2 \right)\right\|_{H^1}\\
&\lesssim \frac{\eps^3}{t^{1/2}}\left(\|\F(M v)\|_{H^1}^2 +
 \eps^2\|\F(Mv_0)\|_{H^1}^2\right)\|\bar B \hat \nu_2
 \|_{H^1}
 \lesssim \frac{\eps^5}{t^{1/2-2C\eps^2}}\|\bar B \hat \nu_2
 \|_{H^1} . 
\end{align*}
The last term is controlled via Leibniz formula, invoking
Lemma~\ref{lem:nu2}, 
\begin{equation*}
\|\bar B \hat \nu_2
 \|_{H^1}\le \|\hat \nu_2 \|_{H^1} + \|\hat \nu_2\d_x
 \bar B\|_{L^2}\lesssim \log t + \eps^2 t^{C\eps^2},
\end{equation*}
and the source terms are estimated by
\begin{equation*}
 \|S_1\|_{L^2} + \|S_2\|_{L^2} \lesssim \eps^5\frac{\log
 t}{t^{1/2-3C\eps^2}}. 
\end{equation*}
We now turn to the $L^2$ estimate of $I_2^\eps-\eps^3 J_2$, and proceed
along the same spirit. 
\begin{equation*}
 I_2^\eps-\eps^3 J_2 = \left|\hat{M(t)v}\right|^2
 \hat{M(t)v} - |\hat v|^2\hat v-\eps^3
 \left|\hat{M(t)v_0}\right|^2 
 \hat{M(t)v_0} +\eps^3 |\hat v_0 |^2\hat v_0 .
\end{equation*}
We distinguish $\big|\hat{M(t)v}\big|^2
 \hat{M(t)v} -\eps^3 \big|\hat{M(t)v_0}\big|^2 
 \hat{M(t)v_0} $ and $|\hat v|^2\hat v-\eps^3
 |\hat v_0 |^2\hat v_0$.
 Discarding the precise algebraic structure like before,
 \begin{equation*}
 \left|\hat{M(t)v}\right|^2
 \hat{M(t)v} -\eps^3 \left|\hat{M(t)v_0}\right|^2 
 \hat{M(t)v_0} \approx \left( \left|\hat{M(t)v}\right|^2+
 \eps^2 \left|\hat{M(t)v_0}\right|^2 \right) \F M(v-\eps v_0). 
 \end{equation*}
The last factor is again rewritten thanks to~\eqref{eq:redress-DA},
and we estimate $I_2^\eps-\eps^3 J_2$ as
\begin{equation*}
 \|I_2^\eps-\eps^3 J_2\|_{L^2}\le \|G_2\|_{L^2}+\|S_3\|_{L^2} + \|S_4\|_{L^2},
\end{equation*}
where
\begin{align*}
 G_2 &= \left( \left|\hat{M(t)v}\right|^2+
 \eps^2 \left|\hat{M(t)v_0}\right|^2 \right) \F M\F^{-1}\left(\bar
 B \left( \hat w -\eps
 \hat v_0-\eps^3\hat \nu_2\right) \right)
%\\ &\quad
- \left( \left|\hat{v}\right|^2+
 \eps^2 \left|\hat{v_0}\right|^2 \right) \bar
 B \left( \hat w -\eps
 \hat v_0-\eps^3\hat \nu_2\right) ,\\
 S_3& = \eps \left( \left|\hat{M(t)v}\right|^2+
 \eps^2 \left|\hat{M(t)v_0}\right|^2 \right) \F M\F^{-1}\left(\left(\bar
 B -1\right) \hat v_0\right)
%\\ &\quad 
-\eps \left( \left|\hat{v}\right|^2+
 \eps^2 \left|\hat{v_0}\right|^2 \right) \left(\bar
 B -1\right) \hat v_0,\\
 S_4& = \eps^3 \left( \left|\hat{M(t)v}\right|^2+
 \eps^2 \left|\hat{M(t)v_0}\right|^2 \right) \F M\F^{-1}\hat \nu_2 -\eps^3 \left( \left|\hat{v}\right|^2+
 \eps^2 \left|\hat{v_0}\right|^2 \right) \hat \nu_2. 
\end{align*}
For the Gronwall term $G_2$, we proceed like for $G_1$, and write
\begin{equation*}
 \|G_2\|_{L^2}\lesssim \eps^2\| \hat w -\eps
 \hat v_0-\eps^3\hat \nu_2\|_{L^2}=\eps^2\|\hat \rho\|_{L^2}.
\end{equation*}
The source terms $S_3$ and $S_4$ are readily of size $\O(\eps^5)$; we
recover some decay in time by making the quantity $M-1$ appear
systematically. In the case of $S_3$, we write
\begin{multline*}
 S_3 = \eps \left( \left|\hat{M(t)v}\right|^2+
 \eps^2 \left|\hat{M(t)v_0}\right|^2 \right) \F M\F^{-1}\left(\left(\bar
 B -1\right) \hat v_0\right)\\
% &\quad 
 -\eps \left( \left|\hat{v}\right|^2+
 \eps^2 \left|\hat{v_0}\right|^2 \right) \left(\bar
 B -1\right) \hat v_0
 \pm \eps \left( \left|\hat{M(t)v}\right|^2+
 \eps^2 \left|\hat{M(t)v_0}\right|^2 \right) \left(\bar
 B -1\right) \hat v_0,
\end{multline*}
and estimate as follows:
\begin{align*}
 \|S_3\|_{L^2} & \lesssim \eps^3 \| (M-1)\F^{-1}\left(\left(\bar
 B -1\right) \hat v_0\right)\|_{L^2} 
%\\ &\quad
 + \|\left(\bar
 B -1\right) \hat v_0\|_{L^\infty} \|(M-1)\hat v\|_{L^2}\left(
 \|\hat{Mv}\|_{L^\infty} + \|\hat{v}\|_{L^\infty}\right)\\
 & \lesssim \frac{\eps^3}{\sqrt t} \| x\F^{-1}\left(\left(\bar
 B -1\right) \hat v_0\right)\|_{L^2}
 + \eps^5 \frac{\log t}{\sqrt t} \times\|x\hat v\|_{L^2},
\end{align*}
where we have used~\eqref{eq:control-hat-v} and~\eqref{eq:B-1}. We
have already estimated the $H^1$-norm of $\left(\bar
B -1\right) \hat v_0$,
\begin{equation*}
 \|\left(\bar B -1\right) \hat v_0\|_{H^1}\lesssim \eps^2
 \left(\log t + t^{C\eps^2}\right),
\end{equation*}
and therefore
\begin{equation*}
 \|S_3\|_{L^2} \lesssim \eps^5\frac{\log t}{t^{1/2- C\eps^2}}. 
\end{equation*}
The term $S_4$ is controlled similarly, by using the same ideas as
above, and we come up with:
\begin{equation*}
 \frac{\dd }{\dd t}\|\hat \rho\|_{L^2}\lesssim
 \frac{\eps^2}{t}\|\hat \rho\|_{L^2} + \eps^5 \frac{\log t}{t^{3/2- 3C\eps^2}}.
\end{equation*}
Gronwall lemma then implies, provided that $\eps>0$ is
sufficiently small:
 \begin{proposition}\label{prop:approx-nu2}
 Suppose that $u^\eps$ solves~\eqref{eq:NLS}, with $u^\eps_{\mid
 t=0}=u_0^\eps\in \Sigma$ such that 
 \begin{equation*}
 u_0^\eps = \eps v_0+\eps^3 w_2 + \O\left(\eps^{5-\eta}\right),
 \end{equation*}
 for $v_0,w_2\in \Sigma$, and some $0<\eta<2$. Then we have
 \begin{equation*}
 \sup_{t\ge 1} \|\hat w (t)-\eps
 \hat v_0-\eps^3\hat \nu_2(t)\|_{L^2}\lesssim \eps^{5-\eta},
 \end{equation*}
 where $\nu_2$ is defined in Lemma~\ref{lem:nu2}. 
 \end{proposition}
 In particular, letting $t$ go to infinity, we infer
 \begin{equation*}
 W =\eps \hat v_0+\eps^3\hat \nu_2^\infty +
 \O\left(\eps^{5-\eta}\right)\quad \text{in }L^2. 
 \end{equation*}
As $u_+^\eps = \F^{-1}\left(W e^{-i\Phi}\right)$, and we have seen at the end
of Section~\ref{sec:complete} that $\|\Phi\|_{L^\infty}=\O(\eps^4)$,
we infer
\begin{equation*}
 u_+^\eps = \eps v_0 + \eps^3\nu_2 + \O\left(\eps^{5-\eta}\right)\quad \text{in }L^2,
\end{equation*}
thus completing the proof of Theorem~\ref{theo:main}. 


\section*{Acknowledgments} The author wishes to thank L\'eo
Bigorgne for discussions on this topic, as well as Satoshi
Masaki and Pavel Naumkin for useful explanations on
Theorem~\ref{theo:complete}. The author is grateful to the anonymous
referee for their careful reading, and for pointing out some flaw in
the initial submission. 


\section*{Declaration of interests}
The authors do not work for, advise, own shares in, or receive funds from
any organization that could benefit from this article, and have declared no
affiliations other than their research organizations.



\printbibliography



\end{document}