Higher order corrections to the mean-field description of the dynamics of interacting bosons

TL;DR

This paper develops a new method to compute higher order corrections to the mean-field approximation for the dynamics of large interacting boson systems, improving the accuracy of modeling Bose-Einstein condensates.

Contribution

It introduces a novel approach using Duhamel expansions and Bogoliubov time evolution to approximate many-body dynamics with arbitrary precision.

Findings

Derived a sequence of N-body functions approximating true dynamics in L^2 norm

Achieved arbitrary order accuracy in powers of N^{-1}

Extended mean-field theory with higher order corrections

Abstract

In this paper, we introduce a novel method for deriving higher order corrections to the mean-field description of the dynamics of interacting bosons. More precisely, we consider the dynamics of $N$ $d$-dimensional bosons for large $N$. The bosons initially form a Bose-Einstein condensate and interact with each other via a pair potential of the form $(N-1)^{-1}N^{d\beta}v(N^\beta\cdot)$ for $\beta\in[0,\frac{1}{4d})$. We derive a sequence of $N$-body functions which approximate the true many-body dynamics in $L^2(\mathbb{R}^{d N})$-norm to arbitrary precision in powers of $N^{-1}$. The approximating functions are constructed as Duhamel expansions of finite order in terms of the first quantised analogue of a Bogoliubov time evolution.

Equations

$$H^{\beta }(t):=\sum _{j=1}^N\left(-{\Delta }_j+V^{\mathrm{ext}}(t,x_j)\right)+\frac{1}{N-1}\sum _{i<j}{v}^{\beta }(x_i-x_j).$$

$$v^{\beta }(x):=N^{d\beta }v(N^{\beta }x),\beta \in [0,\frac{1}{d}),$$

$$v_{ij}^{\beta }:=v^{\beta }(x_i-x_j).$$

$$\mathrm{i}\frac{\mathrm{d}}{\mathrm{d}t}U(t,s)=H^{\beta }(t)U(t,s),U(s,s)=\text{mathbbm}1.$$

$$\psi (t)=U(t,0)\psi (0)$$

$$\Vert \psi (t)-{\psi }_{\mathrm{approx}}(t){\Vert }_{L^2({\mathbb{R}}^{dN})}^2\le C(t)N^{-\delta }$$

$$\Vert \psi (t)-{\psi }_{\phi }^{(a)}(t){\Vert }_{L^2({\mathbb{R}}^{dN})}^2\le C(t)N^{-a\delta (\beta ,\gamma )}$$

$${⟨\cdot ,\cdot ⟩}_{L^2({\mathbb{R}}^{dN})}=:⟨\cdot ,\cdot ⟩,\Vert \cdot {\Vert }_{L^2({\mathbb{R}}^{dN})}=:\Vert \cdot \Vert ,\Vert \cdot {\Vert }_{\mathcal{L}(L^2({\mathbb{R}}^{dN}))}=:\Vert \cdot {\Vert }_{\mathrm{op}}.$$

$$\lfloor r\rfloor :=\max \left\{z\in \mathbb{Z}:z\le r\right\},\lceil r\rceil :=\min \left\{z\in \mathbb{Z}:z>r\right\}$$

$$\underset{N\to \infty }{\lim }\mathrm{Tr}\left|{\gamma }^{(1)}(0)-|{\phi }_0⟩⟨{\phi }_0|\right|=0,$$

$${\gamma }^{(1)}(t):={\mathrm{Tr}}_{2,\text{makebox}[7.00002pt][c].\text{hfil}.\text{hfil}.,N}|\psi (t)⟩⟨\psi (t)|$$

$$\underset{N\to \infty }{\lim }\mathrm{Tr}\left|{\gamma }^{(1)}(t)-|\phi (t)⟩⟨\phi (t)|\right|=0$$

$$\mathrm{i}\frac{\mathrm{d}}{\mathrm{d}t}\phi (t)=\left(-\Delta +V^{\mathrm{ext}}(t)+{\overline{v}}^{\phi (t)}-{\mu }^{\phi (t)}\right)\phi (t)=:h^{\phi (t)}(t)\phi (t)$$

$${\overline{v}}^{\phi (t)}(x):=\left(v^{\beta }\ast |\phi (t){|}^2\right)(x):={\int }_{{\mathbb{R}}^d}{v}^{\beta }(x-y)|\phi (t,y){|}^2\mathrm{d}y.$$

$${\mu }^{\phi (t)}:=\frac{1}{2}{\int }_{{\mathbb{R}}^d}\mathrm{d}x|\phi (t,x){|}^2{\overline{v}}^{\phi (t)}(x)=\frac{1}{2}{\int }_{{\mathbb{R}}^d}\mathrm{d}x{\int }_{{\mathbb{R}}^d}\mathrm{d}y|\phi (t,x){|}^2|\phi (t,y){|}^2{v}^{\beta }(x-y)$$

$$\psi (t)=\sum _{k=0}^N{\phi (t)}^{\otimes (N-k)}{\otimes }_s{\xi }_{\phi (t)}^{(k)}$$

$${\mathcal{F}}_{\perp \phi }^{\le N}:=\underset{k=0}{\overset{N}{⨁}}\underset{\mathrm{sym}}{\overset{k}{⨂}}\{\phi {\}}^{\perp }$$

$$({\psi }_a{\otimes }_s{\psi }_b)(x_1,\text{makebox}[10.00002pt][c].\text{hfil}.\text{hfil}.,x_{a+b}):=\frac{1}{\sqrt{a!b!(a+b)!}}\sum _{\sigma \in {\mathfrak{S}}_{a+b}}{\psi }_a(x_{\sigma (1)},\text{makebox}[10.00002pt][c].\text{hfil}.\text{hfil}.,x_{\sigma (a)}){\psi }_b(x_{\sigma (a+1)},\text{makebox}[10.00002pt][c].\text{hfil}.\text{hfil}.,x_{\sigma (a+b)}),$$

$${\psi }_0=\sum _{k=0}^N{{\phi }_0}^{\otimes (N-k)}{\otimes }_s{\chi }^{(k)}(0)$$

$${\Vert \psi (t)-\sum _{k=0}^N{\phi (t)}^{\otimes (N-k)}{\otimes }_s{\chi }^{(k)}(t)\Vert }_{L^2({\mathbb{R}}^{dN})}^2\le C{\mathrm{e}}^{C^{\prime }t}{N}^{-\delta },$$

$$\mathrm{i}\frac{\mathrm{d}}{\mathrm{d}t}\chi (t)={\mathbb{H}}_{\mathrm{Bog}}(t)\chi (t).$$

$$\begin{array}{r}{\mathbb{H}}_{\mathrm{Bog}}(t):={\int }_{{\mathbb{R}}^d}{a}_x^{\ast }\left(h^{\phi (t)}(t,x)+K_1(t,x)\right)a_x\mathrm{d}x+\frac{1}{2}{\int }_{{\mathbb{R}}^d}\mathrm{d}x{\int }_{{\mathbb{R}}^d}\mathrm{d}y\left(K_2(t,x,y)a_x^{\ast }{a}_y^{\ast }+\overline{K_2(t,x,y)}{a}_x{a}_y\right),\end{array}$$

$$p^{\phi }:=|\phi ⟩⟨\phi |,q^{\phi }:=\text{mathbbm}1-p^{\phi }$$

$$p_j^{\phi }:=\underset{j-1}{\underbrace{\text{mathbbm}1\otimes \cdots \otimes \text{mathbbm}1}}\otimes p^{\phi }\otimes \underset{N-j}{\underbrace{\text{mathbbm}1\otimes \cdots \otimes \text{mathbbm}1}}\text{and}{q}_j^{\phi }:=\text{mathbbm}1-p_j^{\phi }.$$

$$P_k^{\phi }$$

$$\widehat{f^{\phi }}:=\sum _{k=0}^{N}f(k)P_k^{\phi },\widehat{f_j^{\phi }}:=\sum _{n=-j}^{N-j}f(n+j)P_n^{\phi }.$$

$$n(k):=\sqrt{\frac{k}{N}},m(k):=\sqrt{\frac{k+1}{N}}.$$

$$\psi =\sum _{k=0}^N{P}_k^{\phi }\psi =:\sum _{k=0}^N{\phi }^{\otimes (N-k)}{\otimes }_s{\xi }_{\phi }^{(k)}$$

$$P_k^{\phi }\psi (x_1,\text{makebox}[10.00002pt][c].\text{hfil}.\text{hfil}.,x_N)$$

$${\xi }_{\phi }^{(k)}(x_1,\text{makebox}[10.00002pt][c].\text{hfil}.\text{hfil}.,x_k):=$$

Full text

1 Introduction

We consider a system of $N$ bosons in $\mathbb{R}^{d}$, $d\geq 1$, interacting with each other via pair interactions in the mean field scaling regime. The Hamiltonian of the system is given by

[TABLE]

Here, ${V^{\mathrm{ext}}}$ denotes some possibly time-dependent external potential, and the interaction potential $v^{\beta}$ is defined as

[TABLE]

for some bounded, spherically symmetric and compactly supported function $v:\mathbb{R}^{d}\to\mathbb{R}$. In the following, we will make use of the abbreviation

[TABLE]

Note that the prefactor $(N-1)^{-1}$ in front of $v^{\beta}$ is chosen such that the interaction energy and the kinetic energy per particle are of the same order. The mean inter-particle distance is of order $N^{-\frac{1}{d}}$ and therefore much smaller than the range of the interaction, which scales as $N^{-\beta}$. Hence, on average, every particle interacts with many other particles, and the interactions are weak since $(N-1)^{-1}N^{d\beta}\to 0$ as $N\to\infty$. This implies that we consider a mean-field regime. In particular, the case $\beta=0$ is known as the Hartree scaling regime.

We study the time evolution of the $N$-body system for large $N$ when the bosons initially exhibit Bose–Einstein condensation. We impose suitable conditions on the external potential ${V^{\mathrm{ext}}}(t)$ such that $H^{\beta}(t)$ is self-adjoint on $\mathcal{D}(H^{\beta}(t))=H^{2}(\mathbb{R}^{dN})$ for each $t\in\mathbb{R}$. Consequently, $H^{\beta}(t)$ generates a unique family of unitary time evolution operators $\{U(t,s)\}_{t,s\in\mathbb{R}}$ via the Schrödinger equation

[TABLE]

The $N$-body wave function at time $t\in\mathbb{R}$ is determined by

[TABLE]

for some initial datum $\psi(0)=\psi_{0}\in L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$. Due to the interactions, the characterisation of the time evolution $U(t,s)$ is a difficult problem. Even if the system was initially in a factorised state, where all particles are independent of each other, the interactions instantaneously correlate the particles such that an explicit formula for $U(t,s)$ is quite inaccessible.

To describe $U(t,s)$ approximatively, one observes that the dynamics of the many-body system can be decomposed into the dynamics of the condensate wave function $\varphi(t)\in L^{2}(\mathbb{R}^{d})$ and the dynamics of the excitations from the (time-evolved) condensate. The evolution of ${\varphi(t)}$ approximates the $N$-body dynamics $\psi(t)$ in the sense of reduced densities (see Section 2.1). Moreover, if the dynamics of the excitations are suitably approximated and added to the description, one obtains an approximation of $\psi(t)$ with respect to the $L^{2}(\mathbb{R}^{dN})$ norm, in the sense that

[TABLE]

for some power $\delta\in(0,1]$ depending on the choice of $\beta$ (see Section 2.2).

In this paper, we introduce a novel method for deriving a more precise characterisation of the dynamics. This is achieved by constructing a sequence of $N$-body wave functions, which are defined via an iteration of Duhamel’s formula with the time evolution $\widetilde{U}_{\varphi}(t,s)$ generated by an auxiliary Hamiltonian $\widetilde{H}^{\varphi(t)}(t)$ (see (21) for a precise definition). Under the assumption that sufficiently high moments of the number of excitations in the initial state are subleading, we prove higher order corrections to the norm approximation for the scaling regime $\beta\in[0,\frac{1}{4d})$. This is to be understood in the following sense:  we construct a sequence of $N$-body wave functions $\{\psi^{(a)}_{\varphi}(t)\}_{a\in\mathbb{N}}\subset L^{2}(\mathbb{R}^{dN})$ such that, for sufficiently large $N$,

[TABLE]

for some time-dependent constant $C(t)$. The positive exponent $\delta(\beta,\gamma)$ is determined in Theorem 1. It depends on $\beta$ and on a parameter $\gamma$ which is related to the initial number of excitations (see assumption A3).

Let us remark that the approximating functions $\psi^{(a)}_{\varphi}(t)$ are $N$-body wave functions. Since they are explicitly given in terms of the dynamics $\widetilde{U}_{\varphi}(t,s)$ related to the (first order) norm approximation, the functions $\psi^{(a)}_{\varphi}(t)$ are much more accessible than the true dynamics $\psi(t)$. In particular, all higher order corrections can be obtained from the norm approximation $\widetilde{U}_{\varphi}(t,0)\psi_{0}$ by computing an $N$-independent number of integrals.

Moreover, if the initial excitation vector is quasi-free, one can extend the method introduced in this paper. We conjecture that it is possible to approximate all $n$-point correlation functions of the full dynamics $\psi(t)$ to arbitrary precision by expressions depending only on approximations of the 2-point correlation functions, whose computation reduces to solving two coupled linear one-body equations ([28, Equations (17a-b)] and [47, Equation (34)]). This would mean a huge simplification of the full $N$-body problem (3) corresponding to the Hamiltonian (1) for certain initial states, in particular with regard to a numerical analysis, and we plan to show this in a separate paper.

Finally, we note that higher order approximations of the reduced density matrices were obtained by Paul and Pulvirenti in [50] for $\beta=0$ and factorised initial data, based on the method of kinetic errors from the paper by Paul, Pulvirenti and Simonella [51]. For $j\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$, the authors of [50] construct a sequence $\{F^{N,n}_{j}(t)\}_{n\in\mathbb{N}}$ of trace class operators on $L^{2}(\mathbb{R}^{jd})$, which approximate the $j$-particle reduced density matrix $\gamma^{(j)}(t)$ of the system with increasing accuracy up to arbitrary precision. The approximating operators $F^{N,n}_{j}(t)$ can be determined by a number of operations scaling with $n$. They depend on the initial data as well as the knowledge of the solution of the Hartree equation and its linearisation around this solution.

Due to different methods used, it is not straightforward to compare the results of [50] with the results of this paper. However, we list some features of our paper that differ from the operator-based method of kinetic errors [50, 51]. In contrast to the approach in [50], we derive approximations directly for the time-evolved $N$-body wave function, and our construction is implemented as a robust algorithm that requires an $a$-dependent, $N$-independent number of explicit calculations to compute the $a$’th order approximation. Moreover, the results obtained in this paper include positive values of $\beta$ and cover more generic initial data than [50], where the initial state is assumed factorised, i.e., with zero initial excitations.

Notation. In the following, any expression $C$ that is independent of both $N$ and $t$ will be referred to as a constant. Note that constants may depend on all fixed parameters of the model such as ${\varphi_{0}}$, $\psi_{0}$, $v$ and ${V^{\mathrm{ext}}}(0)$. Further, we denote $A\lesssim B$ and $A\gtrsim B$ to indicate that there exists a constant $C>0$ such that $A\leq CB$, resp. $A\geq CB$, and abbreviate

[TABLE]

Finally, we use the notation

[TABLE]

for $r\in\mathbb{R}$.

2 Known results

2.1 Leading order approximation: reduced densities

A first approximation to the $N$-body dynamics is provided by the time evolution of the condensate wave function. Its dynamics yield a macroscopic description of the Bose gas, which, in the limit $N\to\infty$, coincides with the true dynamics in the sense of reduced density matrices. In order to formulate this mathematically, one assumes that the system is initially in a Bose–Einstein condensate with condensate wave function $\varphi_{0}$, i.e.,

[TABLE]

where

[TABLE]

is the one-particle reduced density matrix of $\psi(t)$ at time $t$. Then it has been shown, see e.g. [1, 2, 13, 15, 20, 21, 35, 57], that

[TABLE]

for any $t\in\mathbb{R}$, where ${\varphi(t)}$ is the solution of the Hartree equation

[TABLE]

with initial datum $\varphi(0)=\varphi_{0}$ and with

[TABLE]

Note that for $\beta=0$, the equation (7) is the $N$-independent Hartree (NLH) equation. For $\beta>0$, the evolution is $N$-dependent and converges to the non-linear Schrödinger (NLS) dynamics with $N$-independent coupling parameter $\int v$ in the limit $N\to\infty$. The parameter $\mu^{\varphi(t)}$ is a real-valued phase factor, which we choose as

[TABLE]

for later convenience. For the convergence with respect to reduced densities, this phase is irrelevant since it cancels in the projection $|{\varphi(t)}\rangle\langle{\varphi(t)}|$.

One way to prove the convergence (6), and consequently to derive the NLH/ NLS equation from a system of $N$ bosons, is via the so-called BBGKY111(Bogoliubov-Born-Green-Kirkwood-Yvon) hierarchy, which was prominently used in the works of Lanford for the study of classical mechanical systems in the infinite particle limit [39, 40]. The first derivation of the NLH equation via the BBGKY hierarchy was given by Spohn in [58], and this was further pursued, e.g., in [1, 2, 22, 23]. About a decade ago, Erdős, Schlein and Yau fully developed the BBGKY hierarchy approach to the derivation of the NLH/NLS equation in their seminal works including [20, 21]. Subsequently, a crucial step of this method was revisited by Klainerman and Machedon in [36], based on reformulating combinatorial argument in [20, 21] and a viewpoint inspired by methods of non-linear PDEs. This, in turn, motivated many recent works on the derivation of dispersive PDEs, including [13, 14, 15, 16, 17, 35, 57]. In [55], Rodnianski and Schlein introduced yet another method for proving (6), which uses coherent states on Fock space and was inspired by techniques of quantum field theory and the pioneering work of Hepp [32].

In the context of the current paper, the most relevant works on the derivation of the NLH/NLS equation are due to Pickl [53, 54], who introduced an efficient method for deriving effective equations from the many-body dynamics, transforming the physical idea behind the mean-field description of an $N$-body system into a mathematical algorithm. Instead of describing the condensate as the vacuum of a Fock space of excitations, this approach remains in the $N$-body setting and uses projection operators to factor out the condensate. This strategy was successfully applied to prove effective dynamics for $N$-boson systems in various situations, e.g., [4, 10, 19, 33, 34, 37, 43, 44].

2.2 Next-to-leading order: norm approximation

Whereas closeness in the sense of reduced densities implies that the majority of the particles (up to a relative number that vanishes as $N\to\infty$) is in the state ${\varphi(t)}$, the norm approximation requires the control of all $N$ particles. In particular, this implies that the excitations from the condensate can no longer be omitted from the description. In this sense, the norm approximation of $\psi(t)$ can be understood as next-to-leading order correction to the mean-field description.

A norm approximation for initial coherent states on Fock space was first obtained in [30, 31] by Grillakis, Machedon and Margetis. For initial states $\psi_{0}\in L^{2}(\mathbb{R}^{dN})$ with fixed particle number, Lewin, Nam and Schlein proved in [41] a norm approximation for $\beta=0$ and ${V^{\mathrm{ext}}}=0$ under quite general assumptions on the interaction potential $v$. Nam and Napiórkowski extended this result in [47] to the range $\beta\in[0,\frac{1}{3})$, in [49] to the range $\beta\in[0,\frac{1}{2})$ for the three-dimensional defocusing case, and in [48] to the focusing case in dimensions one and two for $\beta>0$ and $\beta\in(0,1)$, respectively. As proposed in [42], the authors decomposed the $N$-body wave function $\psi(t)$ into condensate and excitations as

[TABLE]

for some \rm=\big{(}\xi^{(k)}_{\varphi(t)}\big{)}_{k=0}^{N}\in{\mathcal{F}^{\leq N}_{\perp{\varphi(t)}}}, where

[TABLE]

is the truncated bosonic Fock space over the orthogonal complement in $L^{2}(\mathbb{R}^{d})$ of the span of $\varphi\in L^{2}(\mathbb{R}^{d})$. A definition of $\xi^{(k)}_{\varphi(t)}$ will be given in (16). Further, $\otimes_{s}$ denotes the symmetric tensor product, which is for $\psi_{a}\in L^{2}(\mathbb{R}^{da})$, $\psi_{b}\in L^{2}(\mathbb{R}^{db})$ defined as

[TABLE]

where $\mathfrak{S}_{a+b}$ denotes the set of all permutations of $a+b$ elements. The addend $k=0$ in (10) describes the condensate, while the terms $k\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$ correspond to the excitations. In the following, we will refer to $\xi^{(k)}_{\varphi}(t)$ as $k$-particle excitation.

In [41, 47, 48, 49], the authors consider initial data of the form

[TABLE]

for some appropriate initial excitation vector $\chi(0):=(\chi^{(k)}(0))^{\infty}_{k=0}\in\mathcal{F}(\{{\varphi_{0}}\}^{\perp})$. It is then shown that there exist constants $C,C^{\prime}>0$ such that

[TABLE]

where $\delta=1-3\beta$ for the three-dimensional defocusing case with $\beta\in[0,\frac{1}{3})$ and $\delta=\frac{1}{2}$ and $\delta<\frac{1}{3}(1-\beta)$ for the one- and two-dimensional focusing case, respectively. The excitations $\chi(t)=(\chi^{(k)}(t))_{k=0}^{\infty}\in\mathcal{F}(\{{\varphi(t)}\}^{\perp})$ at time $t>0$ are determined by the Bogoliubov evolution,

[TABLE]

Here, $\mathbb{H}_{\mathrm{Bog}}(t)$ denotes the Bogoliubov Hamiltonian222 Written in second quantized form, $\mathbb{H}_{\mathrm{Bog}}(t)$ is defined as

\begin{split}\mathbb{H}_{\mathrm{Bog}}(t):=\int_{\mathbb{R}^{d}}a^{*}_{x}\left(h^{{\varphi(t)}}(t,x)+K_{1}(t,x)\right)a_{x}\mathop{}\!\mathrm{d}x+\tfrac{1}{2}\int_{\mathbb{R}^{d}}\mathop{}\!\mathrm{d}x\int_{\mathbb{R}^{d}}\mathop{}\!\mathrm{d}y\left(K_{2}(t,x,y)a^{*}_{x}a^{*}_{y}+\overline{K_{2}(t,x,y)}a_{x}a_{y}\right)\,,\end{split}

where $a^{*}_{x}$ and $a_{x}$ denote the operator-valued distributions corresponding to the usual creation and annihilation operators on $\mathcal{F}(L^{2}(\mathbb{R}^{d}))$. Besides, $K_{1}(t):=Q(t)\widetilde{K}_{1}(t)Q(t)$ with $Q(t):=1-|{\varphi(t)}\rangle\langle{\varphi(t)}|$, where $\widetilde{K}_{1}$ is the Hilbert-Schmidt operator on $L^{2}(\mathbb{R}^{d})$ with kernel $\widetilde{K}_{1}(t,x,y):=\varphi(t,x)v^{\beta}(x-y)\overline{\varphi(t,y)}$. Further, $K_{2}(t):=\left(Q(t)\otimes Q(t)\right)\widetilde{K}_{2}(t,\cdot,\cdot)$, where $\widetilde{K}_{2}(t,x,y):=\varphi(t,x)v^{\beta}(x-y)\varphi(t,y)$ (e.g. [47, Equation (31)]). , an effective Hamiltonian on Fock space which is quadratic in the number of creation and annihilation operators.

For three dimensions and scaling parameter $\beta=0$, a similar result was obtained by Mitrouskas, Petrat and Pickl in [46, 45] via a first quantised approach, where the splitting of $\psi(t)$ into condensate and excitations is realised by means of projections as introduced in [53]. Since we will work in the first quantised setting, let us recall this approach and introduce some notation.

Definition 2.1**.**

Let $\varphi\in L^{2}(\mathbb{R}^{d})$. Define the orthogonal projections on $L^{2}(\mathbb{R}^{d})$

[TABLE]

and the corresponding projection operators on $L^{2}(\mathbb{R}^{dN})$

[TABLE]

For $0\leq k\leq N$, define the many-body projections

[TABLE]

and $P^{\varphi}_{k}=0$ for $k<0$ and $k>N$. Further, for any function $f:\mathbb{N}_{0}\rightarrow\mathbb{R}_{0}^{+}$ and any $j\in\mathbb{Z}$, define the operators $\widehat{f^{\varphi}},\widehat{f^{\varphi}_{j}}\in\mathcal{L}\left(L^{2}(\mathbb{R}^{dN})\right)$ by

[TABLE]

We will in particular need the operators $\widehat{n^{\varphi}}$ and $\widehat{m^{\varphi}}$ corresponding to the weights

[TABLE]

Hence, for any $\psi\in L^{2}(\mathbb{R}^{dN})$, the part of $\psi$ in the condensate $\varphi^{\otimes N}$ is given by $P^{\varphi}_{0}\psi$, and the part of $\psi$ corresponding to $k$ particles being excited from the condensate is precisely $P^{\varphi}_{k}\psi$ for $k\geq 1$. By construction, $P^{\varphi}_{k}P^{\varphi}_{k^{\prime}}=\delta_{k,k^{\prime}}P^{\varphi}_{k}$. Besides, the identity $\sum_{k=0}^{N}P^{\varphi}_{k}=\mathbbm{1}$ implies

[TABLE]

for some $\xi^{(k)}_{\varphi}\in L^{2}(\mathbb{R}^{dk})$. To determine the explicit form of $\xi^{(k)}_{\varphi}$, observe that by Definition 2.1,

[TABLE]

where, by definition of the symmetric tensor product,

[TABLE]

Obviously, $\xi^{(k)}_{\varphi}$ is symmetric under permutations of all of its coordinates, and $\xi^{(k)}_{\varphi}$ is orthogonal to $\varphi$ in every coordinate, i.e.,

[TABLE]

for every $j\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},k\}$. Hence, $\xi^{(k)}_{\varphi}\in\bigotimes_{\mathrm{sym}}^{k}\{\varphi\}^{\perp}$. The excitations $\xi^{(k)}_{\varphi}$, $k\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$, define a vector $\xi_{\varphi}:=\left(\xi^{(0)}_{\varphi},\xi^{(1)}_{\varphi},\makebox[10.00002pt][c]{.\hfil.\hfil.},\xi^{(N)}_{\varphi}\right)$ in the truncated Fock space ${\mathcal{F}^{\leq N}_{\perp\varphi}}$ defined in (11). The relation between the $N$-body state $\psi$ and the corresponding excitation vector $\xi_{\varphi}$ is given by the unitary map

[TABLE]

where $\xi_{\varphi}$ is defined by (16). The vacuum $(1,0,\makebox[10.00002pt][c]{.\hfil.\hfil.},0)$ of ${\mathcal{F}^{\leq N}_{\perp\varphi}}$ corresponds to the condensate $\varphi^{\otimes N}$, and the probability of $k$ particles being outside the condensate equals

[TABLE]

by (16). The number operator $\mathcal{N}_{\varphi}$ on ${\mathcal{F}^{\leq N}_{\perp\varphi}}$, counting the number of excitations, is defined by its action

[TABLE]

The expected number of excitations from the condensate $\varphi^{\otimes N}$ in the state $\psi$ is thus given by

[TABLE]

with $\widehat{n^{\varphi}}$ from Definition 2.1.

In [46], the authors introduce an auxiliary $N$-particle Hamiltonian $\widetilde{H}^{\varphi(t)}(t)$ by subtracting from $H^{\beta}(t)$ in each coordinate the mean-field Hamiltonian $h^{{\varphi(t)}}(t)$ from (7), inserting identities

[TABLE]

on both sides of the difference, and discarding all terms which are cubic, $\mathcal{C}^{\varphi(t)}$, or quartic, $\mathcal{Q}^{\varphi(t)}$, in the number of projections $q^{\varphi(t)}$:

Lemma 2.1**.**

[TABLE]

where

[TABLE]

and with

[TABLE]

Proof.

[TABLE]

Now one inserts identities $\mathbbm{1}=(p^{\varphi(t)}_{i}+q^{\varphi(t)}_{i})(p^{\varphi(t)}_{j}+q^{\varphi(t)}_{j})$ before and after the expression in the brackets and uses the relations

[TABLE]

which concludes the proof. ∎

The auxiliary Hamiltonian $\widetilde{H}^{\varphi(t)}(t)$ has a quadratic structure comparable to that of the Bogoliubov-Hamiltonian $\mathbb{H}_{\mathrm{Bog}}(t)$: all terms in $\widetilde{H}^{\varphi(t)}(t)-\sum_{j}h^{{\varphi(t)}}_{j}(t)$, which form an effective two-body potential, contain exactly two projectors $q^{\varphi(t)}$ onto the complement of the condensate wave function, while $\mathbb{H}_{\mathrm{Bog}}(t)$ is quadratic in the creation and annihilation operators of the excitations. However, $\widetilde{H}^{\varphi(t)}(t)$ is particle number conserving and acts on the $N$-body Hilbert space $L^{2}(\mathbb{R}^{dN})$, i.e., it determines the evolution of the $N$-body wave function consisting of excited particles and particles in the condensate ${\varphi(t)}$, with ${\varphi(t)}$ the solution of (7). In contrast, $\mathbb{H}_{\mathrm{Bog}}(t)$ operates on the excitation Fock space ${\mathcal{F}_{\perp{\varphi(t)}}}$, does not conserve the particle number, and exclusively concerns the dynamics of the excitations with respect to the condensate wave function evolving according to (7).

The time evolution generated by $\widetilde{H}^{\varphi(t)}(t)$ is denoted by $\widetilde{U}_{\varphi}(t,s)$. For an initial datum $\psi_{0}\in L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$, the corresponding $N$-body wave function at time $t\in\mathbb{R}$ is

[TABLE]

Existence and uniqueness of of the time evolution $\widetilde{U}_{\varphi}(t,s)$ are recalled in Lemma 3.1.

Under appropriate assumptions on the initial datum $\psi_{0}$, the time evolution $\widetilde{U}_{\varphi}(t,s)$ approximates the actual time evolution $U(t,s)$. More precisely, there exist constants $C,C^{\prime}>0$ such that

[TABLE]

[46, Theorem 2.6]. Further, in the limit $N\to\infty$, the excitations in $\widetilde{U}_{\varphi}(t,0)\psi_{0}$ coincide with the solutions of the Bogoliubov evolution equation: let \xi_{\varphi_{0}}=\big{(}\xi^{(k)}_{{\varphi_{0}}}\big{)}_{k=0}^{N} denote the excitations from ${\varphi_{0}}^{\otimes N}$ in the initial state $\psi_{0}$ under the decomposition (10), let \widetilde{\xi}_{\varphi(t)}=\big{(}\widetilde{\xi}_{\varphi(t)}^{\,(k)}\big{)}^{N}_{k=0} denote the excitations from ${\varphi(t)}^{\otimes N}$ in $\widetilde{U}_{\varphi}(t,0)\psi_{0}$, and let $\chi(t)=\left(\chi^{(k)}(t)\right)_{k\geq 0}$ denote the solutions of (14) with initial datum $\chi^{(k)}(0)=\xi_{\varphi_{0}}^{(k)}$ for $0\leq k\leq N$ and $\chi^{(k)}(0)=0$ for $k>N$. Then

[TABLE]

[46, Lemma 2.8]. Hence, the combination of (23) and (24) yields (13), with a different time-dependent constant but the same $N$-dependence.

Finally, let us remark that for larger values of the scaling parameter, beyond the mean field regime, the evolutions of ${\varphi(t)}$ and do not (approximately) decouple any more as a consequence of the short-scale structure related to the two-body scattering process. For $\beta\in(0,1)$, an accordingly adjusted variant of (13) for appropriately modified initial data was proved by Brennecke, Nam, Napiórkowski and Schlein in [11] in the three-dimensional defocusing case. Similar estimates for the many-body evolution of appropriate classes of Fock space initial data have been obtained in [8, 18, 24, 25, 28, 29, 30, 31, 38, 55] for various ranges of the scaling parameter. A related result for Bose gases with large volume and large density was proved in [52].

3 Main results

3.1 Assumptions

Let us state our assumptions on the model (1) and on the initial data.

• A1

Interaction potential. Let $v:\mathbb{R}^{d}\to\mathbb{R}$ be spherically symmetric and bounded uniformly in $N$, i.e., $\lVert v\rVert_{L^{\infty}(\mathbb{R}^{d})}\lesssim 1$. Further, assume that $\operatorname*{\mathrm{supp}}v\subseteq\{x\in\mathbb{R}^{d}:|x|\lesssim 1\}$.

• A2

External potential. Let ${V^{\mathrm{ext}}}:\mathbb{R}\times\mathbb{R}^{d}\to\mathbb{R}$ such that ${V^{\mathrm{ext}}}(\cdot,x)\in\mathcal{C}(\mathbb{R})$ for each $x\in\mathbb{R}^{d}$ and ${V^{\mathrm{ext}}}(t,\cdot)\in L^{\infty}(\mathbb{R}^{d})$ for each $t\in\mathbb{R}$.

• A3

Initial data. Let $\psi_{0}\in H^{2}(\mathbb{R}^{dN})\cap L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$ and $\varphi_{0}\in H^{k}(\mathbb{R}^{d})$, $k=\lceil\tfrac{d}{2}\rceil$, both be normalised. Let $\gamma\in(0,1]$ and $A\in\mathbb{N}$. Assume that for any $a\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},A\}$, there exists a set of non-negative, $a$-dependent constants $\left\{\mathfrak{C}_{\,a}\right\}_{0\leq a\leq A}$ with $\mathfrak{C}_{\,0}=1$ such that, for sufficiently large $N$,

[TABLE]

Our analysis is valid as long as the solution ${\varphi(t)}$ of the non-linear equation (7) exists in $H^{k}(\mathbb{R}^{d})$-sense for $k=\lceil\tfrac{d}{2}\rceil$. The maximal time of $H^{k}(\mathbb{R}^{d})$-existence, $T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}$, is defined as

[TABLE]

and depends on the dimension $d$, the sign of $\overline{v}^{\varphi(t)}$, and the regularity of the external trap ${V^{\mathrm{ext}}}(t)$. Under assumptions A1 and A2 and for times $s,t\in[0,T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}})$, the time evolution $\widetilde{U}_{\varphi}(t,s)$ is well-defined.

Lemma 3.1**.**

Let s,t\in\big{[}0,T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}\big{)}. Then $\widetilde{H}^{\varphi(t)}(t)$ is self-adjoint on $\mathcal{D}(\widetilde{H}^{\varphi(t)}(t))=H^{2}(\mathbb{R}^{dN})$ and generates a unique family of unitary time evolution operators $\widetilde{U}_{\varphi}(t,s)$. $\widetilde{U}_{\varphi}(t,s)$ is strongly continuous jointly in $s,t$ and leaves $H^{2}(\mathbb{R}^{dN})$ invariant.

Proof.

As a consequence of the Sobolev embedding theorem (e.g. [3, Theorem 4.12, Part IA]), $\lVert{\varphi(t)}\rVert_{L^{\infty}(\mathbb{R}^{d})}\lesssim\lVert{\varphi(t)}\rVert_{H^{k}(\mathbb{R}^{d})}$ for $k=\lceil\tfrac{d}{2}\rceil$. Hence, by definition (25) of $T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}$, $\mu^{\varphi(t)}$ and $(N-1)\overline{v}^{\varphi(t)}$ are bounded uniformly in $N$ for t\in\big{[}0,T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}\big{)}. Further, $t\mapsto\widetilde{H}^{\varphi(t)}(t)\psi$ is Lipschitz for all $\psi\in H^{2}(\mathbb{R}^{dN})$ because of (7), since $t\mapsto{V^{\mathrm{ext}}}(t)\in\mathcal{L}(L^{2}(\mathbb{R}^{d}))$ is continuous and as $\tfrac{\mathop{}\!\mathrm{d}}{\mathop{}\!\mathrm{d}t}p^{\varphi(t)}=\mathrm{i}[p^{\varphi(t)},h^{{\varphi(t)}}(t)]$. Hence, the statement of the lemma follows from [27]. ∎

Assumptions A1 and A2 are rather standard in the rigorous treatment of interacting many-boson systems. Note that we make no assumption on the sign of the potential or its scattering length and thus cover both repulsive and attractive interactions. Besides, we admit a large class of time-dependent external traps ${V^{\mathrm{ext}}}$, with basically the only restriction that ${V^{\mathrm{ext}}}(t)$ must not obstruct the self-adjointness of $H^{\beta}(t)$ on $H^{2}(\mathbb{R}^{dN})$.

The third assumption provides a bound on the expected number of excitations from the condensate ${\varphi_{0}}^{\otimes N}$ in the initial state $\psi_{0}$. Note that while $\gamma=0$ is the trivial bound, the condition becomes more restrictive as $\gamma$ increases. We have chosen this particular formulation of A3 for later convenience333 Note that the operators $\widehat{n^{\varphi}}$ and $\widehat{m^{\varphi}}$ are equivalent in the sense that they are related via (35), namely $(\widehat{n^{\varphi}})^{2a}\leq(\widehat{m^{\varphi}})^{2a}\leq 2^{a}(\widehat{n^{\varphi}})^{2a}+N^{-a}$, hence all results in terms of $\widehat{m^{\varphi}}$ can be translated to the corresponding statements in terms of $\widehat{n^{\varphi}}$. We chose to work with $\widehat{m^{\varphi}}$ instead of $\widehat{n}^{\varphi}$ because this makes in particular Proposition 3.3 easier to write. For example, in terms of $\widehat{n^{\varphi}}$, Proposition 3.3b reads

$\lVert(\widehat{n^{\varphi}})^{j}\widetilde{U}_{\varphi}(t,s)\psi\rVert^{2}\lesssim C_{j}^{\,t,s}\sum\limits_{n=0}^{j}N^{n(-1+d\beta)}\left(2^{j-n}\lVert(\widehat{n^{\varphi}})^{j-n}\psi\rVert^{2}+N^{-j+n}\right)\,,$

which contains an additional term $N^{-j+n}$. Since the proof of our main result requires an iteration of this proposition, the version with $\widehat{m^{\varphi}}$ is more practicable. . However, its physical meaning is better understood from one of the following two equivalent versions of A3:

• A3$\,{}^{\prime}$

Let $\psi_{0}\in H^{2}(\mathbb{R}^{dN})\cap L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$ and $\varphi_{0}\in H^{k}(\mathbb{R}^{d})$, $k=\lceil\tfrac{d}{2}\rceil$, both be normalised. Let $\gamma\in(0,1]$ and $A\in\mathbb{N}$. Assume that for any $a\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},A\}$, there exists a set of non-negative, $a$-dependent constants $\{\mathfrak{C}_{\,a}^{\prime}\}_{0\leq a\leq A}$ with $\mathfrak{C}_{\,0}^{\prime}=1$ such that, for sufficiently large $N$,

[TABLE]

• A3$\,{}^{\prime\prime}$

Let $\psi_{0}\in H^{2}(\mathbb{R}^{dN})\cap L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$ and $\varphi_{0}\in H^{k}(\mathbb{R}^{d})$, $k=\lceil\tfrac{d}{2}\rceil$, both be normalised. Let $\gamma\in(0,1]$, $A\in\mathbb{N}$ and $\xi_{\varphi_{0}}=\mathfrak{U}_{N}^{{\varphi_{0}}}\psi_{0}$. Assume that for any $a\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},A\}$, there exists a set of non-negative, $a$-dependent constants $\{\mathfrak{C}_{\,a}^{\prime\prime}\}_{0\leq a\leq A}$ with $\mathfrak{C}_{\,0}^{\prime\prime}=1$ such that, for sufficiently large $N$,

[TABLE]

The equivalence $\mathrm{A3}\Leftrightarrow\mathrm{A3}\,^{\prime}\Leftrightarrow\mathrm{A3}\,^{\prime\prime}$ follows immediately from Lemma 3.2, whose proof is postponed to Section 4.1.

Lemma 3.2**.**

Let $a\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$ and $\varphi\in L^{2}(\mathbb{R}^{d})$. Let $\psi\in L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$ and $\xi_{\varphi}=\mathfrak{U}_{N}^{\varphi}\psi$. Then

1. (a)

\big{\|}q^{\varphi}_{1}\makebox[10.00002pt][c]{\cdot$\hfil$\cdot$\hfil$\cdot}q^{\varphi}_{a}\psi\big{\|}^{2}\;\leq\;\left\|\left(\widehat{m^{\varphi}}\right)^{a}\psi\right\|^{2}\;\leq\;4^{a}a!\,\sum\limits_{j=1}^{a}N^{-a+j}\big{\|}q^{\varphi}_{1}\makebox[10.00002pt][c]{\cdot$\hfil$\cdot$\hfil$\cdot}q^{\varphi}_{j}\psi\big{\|}^{2}+N^{-a}\,,** 2. (b)

$\left\langle\xi_{\varphi},\mathcal{N}_{\varphi}^{a}\,\xi_{\varphi}\right\rangle_{\mathcal{F}^{\leq N}_{\perp\varphi}}\;\leq\;N^{a}\left\|\left(\widehat{m^{\varphi}}\right)^{a}\psi\right\|^{2}\;\leq\;1+2^{a}\left\langle\xi_{\varphi},\mathcal{N}_{\varphi}^{a}\,\xi_{\varphi}\right\rangle_{\mathcal{F}^{\leq N}_{\perp\varphi}}$* .*

Hence, A3 can be understood as follows: Let $A\in\mathbb{N}$ and consider sufficiently large $N$ such that $A=\mathcal{O}(1)$ with respect to $N$, i.e. $A\lesssim 1$. Then we assume that for any $a\leq A$, the part of the wave function with any $a$ particles outside the condensate is at most of order $N^{-\gamma a}$.

Equivalently, A3 states that the first $A\lesssim 1$ moments of the number of excitations must be sub-leading with respect to the particle number; for $\gamma=1$, they must even be bounded uniformly in $N$. Here, “sub-leading” means that the moments of the relative number of excitations, i.e., the expectation values of $(\mathcal{N}_{\varphi(t)}/N)^{A}$, vanish as $N\to\infty$. This, in turn, provides a bound on the high components of the excitation vector: for example, $\sum_{k=0}^{N}k^{A}\lVert\xi^{(k)}_{{\varphi_{0}}}\rVert^{2}_{L^{2}(\mathbb{R}^{dk})}\lesssim N^{(1-\gamma)A}$ implies $\lVert\xi^{(N)}_{\varphi_{0}}\rVert^{2}_{L^{2}(\mathbb{R}^{dk})}\lesssim N^{-\gamma A}$. In other words, it must be very unlikely that significantly many particles are outside the condensate, whereas we impose no restriction on excitations involving only few particles (with respect to $N$).

As soon as $a$ becomes comparable to $N$, i.e., $a\gtrsim N$, the constants $\mathfrak{C}_{a}^{(\prime,\prime\prime)}$ are $N$-dependent and the assumption is trivially satisfied. However, note that we demand that $N$ be large enough that $A\lesssim 1$.

The simplest example of an $N$-body state satisfying A3 is the product state $\psi={\varphi_{0}}^{\otimes N}$. Whereas the ground state of non-interacting bosons ($v=0$) is of this form, the ground state as well as the lower excited states of interacting systems are not close to an exact product with respect to the $L^{2}(\mathbb{R}^{dN})$-norm due to the correlation structure related to the interactions.

Regarding interacting bosons, A3 is fulfilled for quasi-free states with subleading expected number of excitations, since for any quasi-free state $\xi\in\mathcal{F}$ and any $a\geq 1$ there exists a constant $C_{a}>0$ such that

[TABLE]

(e.g. [47, Lemma 5]). In [42, Theorem A.1], it was shown that the ground state of $\mathbb{H}_{\mathrm{Bog}}$ is a quasi-free state, which, via the map $\mathfrak{U}_{N}^{\varphi}$, defines an $N$-body state $\psi_{\mathrm{Bog}}$ that converges to the actual ground state $\psi_{0}$ in norm as $N\to\infty$ [42, Theorem 2.2]. Note that we require a certain minimal size of $\gamma$, which is strictly greater than $\tfrac{2}{3}$.

Further, let us remark that assumption A3 for $A=1$ means complete BEC, which was shown to be a sufficient condition for the validity of the Bogoliubov approximation [42]. It was shown in [56, Lemma 1] and [26, Lemma 1] that A3 with $A=1$ and $\gamma=1$ is satisfied for low-energy states of a $d$-dimensional Bose gas for $\beta=0$ in the homogeneous and inhomogeneous setting, respectively. For the Gross–Pitaevskii scaling $\beta=1$ in three dimensions, a comparable bound was proved in [6, 7].

Besides, A3 with parameter $\gamma<1$ is comparable to part of the assumption made in [49]. In this work, the authors assume that the initial excitation vector $\xi_{\varphi_{0}}$ be a quasi-free state in $\mathcal{F}_{\perp{\varphi(t)}}$ such that

[TABLE]

for all $\varepsilon>0$ and with $\kappa_{\varepsilon}>0$ independent of $N$, and prove a norm approximation for $\beta\in[0,\frac{1}{2})$ with parameter $\delta=(1-\varepsilon-2\beta)/2$. Since $\xi_{\varphi_{0}}$ is quasi-free, the first part of this assumption is comparable to A3, whereas we do not impose any condition on the initial energy of the excitations.

Finally, Mitrouskas showed in [45, Chapter 3] that assumption A3 with $\gamma=1$ (and consequently for all $\gamma\in(0,1]$) is fulfilled by the ground state and lower excited states of a homogeneous Bose gas on the $d$-dimensional torus for $\beta=0$. More precisely, let ${\varphi_{0}}$ be the minimizer of the Hartree functional on the torus with ground state energy $E_{0}$, and let $\psi_{n}$ denote the $n$’th excited state with energy $E_{n}$. Then the author proves that there exist constants $C,D>0$ such that $\lVert P^{\varphi_{0}}_{a}\psi_{n}\rVert^{2}\leq C\mathrm{e}^{-Da}$ for all $(E_{n}-E_{0})\leq a\leq N$. As a corollary of this statement, it is shown that there exists $C_{a}>0$ such that

[TABLE]

which implies that assumption A3$\,{}^{\prime}$ is satisfied.

3.2 Control of higher moments of the number of excitations

In our first result, we estimate the growth of the first $A$ moments of the number of excitations when the system evolves under the dynamics $U(t,s)$ or $\widetilde{U}_{\varphi}(t,s)$. Estimates of this kind are often needed to derive effective descriptions of the dynamics of interacting bosons, e.g., in [5, 8, 12, 46, 52, 55]. Our proof extends comparable statements for $\beta=0$ and $d=3$ obtained in [46, Lemma 2.1] and [55, Proposition 3.3], and for Bose gases with large volume and large density in [52, Corollary 4.2]. The estimates are stated for $\lVert(\widehat{m^{\varphi}})^{a}\psi\rVert^{2}$ as these expressions are required for the proof of our main theorem. By Lemma 3.2, they easily translate to bounds on the corresponding quantities $\lVert q_{1}\cdots q_{a}\psi\rVert^{2}$ and $\left\langle\xi_{\varphi},\mathcal{N}_{\varphi}^{a}\xi_{\varphi}\right\rangle$. The proofs of Proposition 3.3 and Corollary 3.4 are postponed to Section 4.2.

Proposition 3.3**.**

Let $\beta\in[0,\frac{1}{d})$, assume A1 and A2 and let $\psi\in L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$. Let $s\in\mathbb{R}$, ${\varphi(s)}\in H^{k}(\mathbb{R}^{d})$ for $k=\lceil\tfrac{d}{2}\rceil$, and let ${\varphi(t)}$ be the solution of (7) with initial datum ${\varphi(s)}$. Then it holds for t\in\big{[}s,s+T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}\big{)} and $j\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$ that

1. (a)

for any $b\in\mathbb{N}_{0}$,

[TABLE] 2. (b)

[TABLE]

where $C_{j}^{\,t,s}:=j!\,3^{j(j+1)}\mathrm{e}^{9^{j}\int_{s}^{t}\lVert\varphi(s_{1})\rVert^{2}_{H^{k}(\mathbb{R}^{d})}\mathop{}\!\mathrm{d}s_{1}}$.

Under the additional assumption A3 on the initial data, this implies that at any time $t$ and for sufficiently large $N$, the first $A$ moments of the number of excitations remain sub-leading:

Corollary 3.4**.**

Assume A1 – A2 and A3 with $\gamma\in(0,1]$ and $A\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$. Let $\psi(t)$, $\widetilde{\psi}_{\varphi}(t)$ and ${\varphi(t)}$ denote the solutions of (4), (22) and (7) with initial data $\psi_{0}$ and $\varphi_{0}$ from A3. Let

[TABLE]

Then for t\in\big{[}0,T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}\big{)}, sufficiently large $N$ and $a\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},A\}$, it holds that

1. (a)

for the time evolution $U(t,s)$ that

[TABLE]

or, equivalently, that

[TABLE] 2. (b)

for the time evolution $\widetilde{U}_{\varphi}(t,0)$ and $\beta\in[0,\tfrac{1}{d})$ that

[TABLE]

or, equivalently, that

[TABLE]

with $C_{a}^{\,t}:=C_{a}^{t,0}$ and where we estimated $a,\mathfrak{C}_{\,a},\mathfrak{C}_{\,a}^{\prime\prime}\lesssim 1$ for the sake of readability.

At the threshold $\gamma=1-d\beta$, the leading order terms in the sums in Proposition 3.3 change, hence we obtain two different estimates. The additional restrictions on $\beta$ and $\gamma$ in part (a) stem from the second sum in Proposition 3.3a. Only if either $\beta<\frac{1}{2d}$ or $\gamma>d\beta$, it is possible to choose $b$ sufficiently large that the first sum dominates for large $N$. For $\beta=0$, both time evolutions preserve the property A3 ′′ exactly, i.e., with the same power $\gamma$ of $N$, up to a constant growing rapidly in $t$ and $a$. For $\beta>0$, the conservation is exact only for small $\gamma$, whereas one looses some power of $N$ for larger $\gamma$. Further, note that for the range $\gamma\in(0,d\beta]$, we do not obtain a non-trivial estimate for the excitations in $U(t,0)\psi_{0}$.

3.3 Higher order corrections to the norm approximation

Based on the estimates obtained in Proposition 3.3, our main result establishes corrections of any order to the norm approximations (13) and (23): under assumption A3 on the initial data, we construct a sequence $\{\psi^{(a)}_{\varphi}\}_{a\in\mathbb{N}}\subset L^{2}(\mathbb{R}^{dN})$ such that

[TABLE]

for some $\delta(\beta,\gamma)>0$, which may depend on $\beta$ as well as on the parameter $\gamma$ from assumption A3. For reasons given below, our analysis is restricted to the scaling regime $\beta\in[0,\frac{1}{4d})$.

As explained in the introduction, it is well known that the actual time evolution $\psi(t)$ is close to the evolution $\widetilde{\psi}_{\varphi}(t)$ from (22) in norm. Hence, the first element of the approximating sequence $\{\psi^{(a)}_{\varphi}\}_{a\in\mathbb{N}}$ is determined by

[TABLE]

Using Duhamel’s formula, the difference between $U(t,s)\psi$ and $\widetilde{U}_{\varphi}(t,s)\psi$ can be expressed as

[TABLE]

for any $\psi\in L^{2}(\mathbb{R}^{dN})$. Consequently,

[TABLE]

by the triangle inequality and as a consequence of the unitarity of $U(t,s)$. The leading order contribution in (27) is the term containing $\mathcal{C}^{\varphi(s)}$ because the cubic interaction terms are larger than the quartic ones in the following sense:

Lemma 3.5**.**

Let $\psi\in L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$ and denote by ${\varphi(t)}$ the solution of (7) with initial datum ${\varphi_{0}}\in H^{k}(\mathbb{R}^{d})$, $k=\lceil\tfrac{d}{2}\rceil$. Then for any $j\in\mathbb{N}_{0}$ and t\in\big{[}0,T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}\big{)},

1. (a)

\lVert\left(\widehat{m^{\varphi(t)}}\right)^{j}\mathcal{Q}^{\varphi(t)}\psi\rVert^{2}\lesssim N^{2+2d\beta}\lVert\big{(}\widehat{m^{\varphi(t)}}\big{)}^{4+j}\psi\rVert^{2}, 2. (b)

\lVert\left(\widehat{m^{\varphi(t)}}\right)^{j}\mathcal{C}^{\varphi(t)}\psi\rVert^{2}\lesssim 4^{j}\lVert{\varphi(t)}\rVert^{2}_{H^{k}(\mathbb{R}^{d})}N^{2+d\beta}\lVert\big{(}\widehat{m^{\varphi(t)}}\big{)}^{3+j}\psi\rVert^{2}.

The proof of this lemma is postponed to Section 4.3. For $j=0$, it gives a bound on the cubic and quartic terms; the more general statement $j\geq 0$ is included for later convenience.

When applying Lemma 3.5 to (27), we obtain expressions of the form $\lVert(\widehat{m^{\varphi(s)}})^{j}\widetilde{U}_{\varphi}(s,0)\psi_{0}\rVert^{2}$. To be able to use assumption A3 on the initial data, we need to interchange, in a sense, the order of $\widetilde{U}_{\varphi}(s,0)$ and $(\widehat{m^{\varphi(s)}})^{j}$. This is where Proposition 3.3 comes into play: from part 3.3b, it follows for sufficiently large $N$ that

[TABLE]

To enhance readability, we do not keep track of the different constants $C(t)$ for now, but we specify it in more detail in Theorem 1. As in Corollary 3.4, the size of $\gamma$ determines the leading order term in the sum: for $\gamma\geq 1-d\beta$, the dominant contribution issues from $n=3$, whereas otherwise the addend corresponding to $n=0$ is of leading order. Consequently,

[TABLE]

To ensure that (28) converges to zero as $N\to\infty$, we restricted the range of parameters $\gamma$ admitted by assumption A3 to $\gamma\in(\tfrac{2+d\beta}{3},1]$. Besides, in the first case, the bound is only small for $\beta<\frac{1}{4d}$, and the second case is anyway only possible for $\beta<\tfrac{1}{4d}$. Hence, we can only cover the parameter regime $\beta\in[0,\tfrac{1}{4d})$. Analogously to (28), we also obtain

[TABLE]

Note that $\beta<\tfrac{1}{4d}$ implies that $-2+6d\beta<-1+4d\beta$, and besides, it follows from $\gamma>\tfrac{2+3d}{3}$ and $\beta<\tfrac{1}{4d}$ that $2+2d\beta-4\gamma<2+d\beta-3\gamma$. Consequently, the contribution with $\mathcal{C}^{\varphi(s)}$ dominates in (27) for sufficiently large $N$, which leads to the estimate

[TABLE]

with

[TABLE]

This yields (5) for $n=1$.

To construct the second element $\psi^{(2)}_{\varphi}(t)$ of the approximating sequence, we need to extract from (26) the relevant contributions such that $\lVert\psi(t)-\psi^{(2)}_{\varphi}(t)\rVert^{2}\leq C(t)N^{-2\delta(\beta,\gamma)}$. As a consequence of Lemma 3.5, we define

[TABLE]

which equals the leading order contribution in (26) but with the true time evolution $U(t,s)$ replaced by $\widetilde{U}_{\varphi}(t,s)$. Put differently, the leading order contribution is cancelled but for the difference between $U(t,s)$ and $\widetilde{U}_{\varphi}(t,s)$. Since this difference is evaluated on $\mathcal{C}^{\varphi(s)}\widetilde{U}_{\varphi}(s,0)\psi_{0}$, which is small in norm, this is an improvement compared to the first order approximation $\psi^{(1)}_{\varphi}(t)$. To verify this, let us compute the difference between $\psi(t)$ and $\psi^{(2)}_{\varphi}(t)$. Using Duhamel’s formula twice, we obtain

[TABLE]

Due to the unitarity of $U(t,s)$, we obtain with the triangle inequality

[TABLE]

The leading order term in (32) can be estimated as

[TABLE]

As before, considering the two ranges of $\gamma$ separately yields for sufficiently large $N$

[TABLE]

with $\delta(\beta,\gamma)$ from (31). Analogously, the second term can be estimated as

[TABLE]

and the third term was already treated in (29). Combining all bounds, we obtain

[TABLE]

which yields (5) for $a=2$.

Iterating Duhamel’s formula $(a-1)$ times, we construct $\psi^{(a)}_{\varphi}(t)$ as an expansion with $a-1$ terms, where the last term contains the true time evolution $U(t,s)$ and all others exclusively contain $\widetilde{U}_{\varphi}(t,s)$. Consequently, to construct $\psi^{(3)}_{\varphi}(t)$, we iterate (26) once more, which yields

[TABLE]

The leading order contributions issue from the first integral and from the expression with two cubic interaction terms. Analogously to above, they determine the next element $\psi^{(3)}_{\varphi}$ of the sequence $\{\psi^{(a)}_{\varphi}\}_{a\in\mathbb{N}}$ as

[TABLE]

and similar calculations as before yield $\lVert\psi(t)-\psi^{(3)}_{\varphi}(t)\rVert^{2}\lesssim C(t)N^{-3\delta(\beta,\gamma)}$. Continuing the iteration of (26), we obtain for any $a\geq 1$ and $s_{0}=0$ the expansion

[TABLE]

All products are to be understood as ordered, i.e. $\prod_{\ell=0}^{L}P_{\ell}:=P_{0}P_{1}\cdots P_{L}$ for $L\in\mathbb{N}$ and any expressions $P_{\ell}$. Extracting the leading contributions in each order, we construct the sequence $\{\psi^{(a)}_{\varphi}(t)\}_{a\in\mathbb{N}}$ as follows:

Definition 3.1**.**

Let $I_{1}^{\varphi(t)}:=\mathcal{C}^{\varphi(t)}$ and $I_{2}^{\varphi(t)}:=\mathcal{Q}^{\varphi(t)}$. Define the set

[TABLE]

i.e., the set of $n$-tuples with elements in $\{1,2\}$ such that the elements of each tuple add to $k$. Define for $n\in\mathbb{N}$ and $n\leq k\leq 2n$

[TABLE]

where $s_{0}:=0$. As above, the products are ordered. For $n=k=0$, let $T_{0}^{(0)}:=\widetilde{U}_{\varphi}(t,0)\psi_{0}$, and $T_{n}^{(k)}:=0$ for $k<n$ and $k>2n$. Hence, $T_{n}^{(k)}$ is an $n$-dimensional integral where the integrand contains all possible combinations of $I^{\varphi(s_{l})}_{j_{l}}$ such that $\sum_{l=1}^{n}j_{l}=k$.

Finally, the elements of the sequence $\{\psi^{(a)}_{\varphi}\}_{a\in\mathbb{N}}$ are defined as

[TABLE]

Theorem 1**.**

Let $\beta\in[0,\frac{1}{4d})$ and assume A1 – A3 with $A\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$ and $\gamma\in(\tfrac{2+d\beta}{3},1]$. Let $\psi(t)$ and ${\varphi(t)}$ denote the solutions of (4) and (7) with initial data $\psi_{0}$ and $\varphi_{0}$ from A3, respectively, and let $\psi^{(a)}_{\varphi}(t)$ be defined as in Definition 3.1. Then for sufficiently large $N$, t\in\big{[}0,T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}\big{)} and $a\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},\lfloor\tfrac{A}{6}\rfloor\}$, there exists a constant $c(a)$ such that

[TABLE]

where

[TABLE]

Hence, given any desired precision of the approximation, there exists some $a\in\mathbb{N}$ such that the corresponding function $\psi^{(a)}_{\varphi}(t)$ approximates the actual $N$-body dynamics $\psi(t)$ to this order for large $N$. To compute $\psi^{(a)}_{\varphi}(t)$, an $a$-dependent number of steps is required, as well as the knowledge of the first quantised Bogoliubov time evolution. We cover initial states where the first $A$ moments of the number of excitations are sub-leading, where $A$ depends on $a$ but is independent of $N$.

4 Proofs

4.1 Preliminaries

Lemma 4.1**.**

Let ${\varphi_{0}}\in H^{k}(\mathbb{R}^{d})$ for $k=\lceil\tfrac{d}{2}\rceil$, t\in\big{[}0,T^{\mathrm{ex}}_{d,v,{V^{\mathrm{ext}}}}\big{)} and ${\varphi(t)}$ the solution of (7) with inital datum ${\varphi_{0}}$.

1. (a)

Let $f:\mathbb{R}^{d}\times\mathbb{R}^{d}\rightarrow\mathbb{R}$ be a measurable function such that $|f(z_{j},z_{k})|\leq F(z_{k}-z_{j})$ almost everywhere for some $F:\mathbb{R}^{d}\rightarrow\mathbb{R}$. Then

[TABLE] 2. (b)

Let $f:\mathbb{N}_{0}\rightarrow\mathbb{R}^{+}_{0}$. Then P^{\varphi(t)}_{k},\;\widehat{f^{\varphi(t)}}\in\mathcal{C}^{1}\big{(}\mathbb{R},\mathcal{L}\left(L^{2}(\mathbb{R}^{dN})\right)\big{)} for $0\leq k\leq N$ and

[TABLE]

where $h^{{\varphi(t)}}_{j}(t)$ denotes the one-particle operator $h^{{\varphi(t)}}(t)$ from (7) acting on the $j$th coordinate.

Proof.

For part (a), see, e.g., [54, Lemma 4.1] and note that $\lVert{\varphi(t)}\rVert_{L^{\infty}(\mathbb{R}^{d})}\lesssim\lVert{\varphi(t)}\rVert_{H^{k}(\mathbb{R}^{d})}$ by the Sobolev embedding theorem. Part (b) can be shown as in the proof of [54, Lemma 6.2]. ∎

Lemma 4.2**.**

Let $\psi\in L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$, $\varphi\in L^{2}(\mathbb{R}^{d})$ and $f:\mathbb{N}_{0}\rightarrow\mathbb{R}_{0}^{+}$.

1. (a)

$\left(\widehat{n^{\varphi}}\right)^{2}=\frac{1}{N}\sum\limits_{j=1}^{N}q^{\varphi}_{j}\,.$** 2. (b)

Let $a\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$. Then for $j\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},a\}$,

[TABLE] 3. (c)

In particular, this implies

[TABLE]

Proof.

For simplicity, let us drop all superscripts $\varphi$. Part (a) is shown e.g. in [54, Lemma 4.1]. For part (b), observe that for any $1\leq j\leq N$,

[TABLE]

by part (a). Since $\widehat{n}\psi$ is again symmetric, the statement follows by iteration. ∎

Lemma 4.3**.**

Denote by $T_{ij}$ an operator acting non-trivially only on coordinates $i$ and $j$.

1. (a)

Let $\varphi\in L^{2}(\mathbb{R}^{d})$, let $f,g:\mathbb{N}_{0}\rightarrow\mathbb{R}_{0}^{+}$ be any weights and $i,j\in\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$. Let $Q^{\varphi}_{0}:=p^{\varphi}_{i}p^{\varphi}_{j}$, $Q^{\varphi}_{1}\in\{p^{\varphi}_{i}q^{\varphi}_{j},q^{\varphi}_{i}p^{\varphi}_{j}\}$ and $Q^{\varphi}_{2}:=q^{\varphi}_{i}q^{\varphi}_{j}$. Then, for $\mu,\nu\in\{0,1,2\}$,

[TABLE] 2. (b)

Let $\Gamma,\Lambda\in L^{2}(\mathbb{R}^{dN})$ be symmetric under the exchange of coordinates in a subset $\mathcal{M}\subseteq\{1,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$ such that $j\notin\mathcal{M}$ and $k,l\in\mathcal{M}$. Then

[TABLE]

Proof.

[54, Lemma 4.1] and [9, Lemma 4.7]. ∎

Proof of Lemma 3.2. Let us for simplicity drop all superscripts $\varphi$. First, observe that

[TABLE]

hence

[TABLE]

in the sense of operators. The first part of (a) follows from Lemma 4.2b and the first line in (34). For the second part, Lemma 4.2a implies

[TABLE]

for $a_{1},\makebox[10.00002pt][c]{.\hfil.\hfil.},a_{N}\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},a\}$. Due to the symmetry of $\psi$, since there are $\binom{a-1}{j-1}$ possibilities to write $a$ as the sum of $j$ positive integers and with $\binom{a}{a_{1},\makebox[7.00002pt][c]{.\hfil.\hfil.},a_{N}}\leq a!$, this yields

[TABLE]

Further, note that

[TABLE]

and $\tbinom{N}{j}\leq N^{j}$, hence

[TABLE]

Part (b) follows from (20) and (35). ∎

4.2 Proof of Proposition 3.3

Proof of Proposition 3.3. The proof of this proposition is essentially an adaptation of the proof of [52, Corollary 4.2]. We begin with part (a). Let $\psi\in L^{2}(\mathbb{R}^{dN})$ symmetric, $s\in\mathbb{R}$ and $f:\mathbb{N}_{0}\to\mathbb{R}_{0}^{+}$ some weight function. Define

[TABLE]

and

[TABLE]

Let us for the moment abbreviate $U(t,s)\psi=:\psi_{t}$. By Lemma 4.1b,

[TABLE]

where we have inserted $\mathbbm{1}=(p^{\varphi(t)}_{1}+q^{\varphi(t)}_{1})(p^{\varphi(t)}_{2}+q^{\varphi(t)}_{2})$ on both sides of the commutator and used Lemma 4.3a. Since $q^{\varphi(t)}_{1}p^{\varphi(t)}_{2}Z^{\beta}_{12}p^{\varphi(t)}_{1}p^{\varphi(t)}_{2}=0$, we conclude that (41) equals zero. From now on, we will for simplicity drop the superscripts ${\varphi(t)}$. Let

[TABLE]

Since, for example, $\left(\widehat{f}-\widehat{f_{-2}}\right)^{\frac{1}{2}}q_{1}q_{2}=\left(\sum\limits_{k=2}^{N}(f(k)-f(k-2))P^{\varphi(t)}_{k}\right)^{\frac{1}{2}}q_{1}q_{2}$, this yields

[TABLE]

By Lemmas 4.1 and 4.2 and since $\lVert v^{\beta}\rVert_{L^{2}(\mathbb{R}^{d})}^{2}\lesssim N^{d\beta}$, the first term in (43) leads to

[TABLE]

To obtain the estimate in the last line, note first that

[TABLE]

Now we decompose $v^{\beta}=v^{\beta}_{+}-v^{\beta}_{-}$ into its positive and negative part such that $v^{\beta}_{\pm}\geq 0$, hence $v^{\beta}_{\pm}(x)=\sqrt{v^{\beta}_{\pm}(x)}\sqrt{v^{\beta}_{\pm}(x)}$, which leads to

[TABLE]

by Lemma 4.1. The second term in (43) can be estimated as

[TABLE]

Now we choose for $f$ the family of weight functions $w_{\lambda}^{j}:k\mapsto\left(w_{\lambda}(k)\right)^{j}$ given by

[TABLE]

for some $0<\lambda\leq 1-d\beta$ and $j\in\{0,\makebox[10.00002pt][c]{.\hfil.\hfil.},N\}$. The set corresponding to $L_{f}$ from (42) is called $L_{w^{j}_{\lambda}}$. To bound the operators in $L_{w^{j}_{\lambda}}$, note that for any $a,b\in\mathbb{N}_{0}$, $a>b$,

[TABLE]

where we have used in the second line that for every $1\leq m\leq j-1$,

[TABLE]

and that $a^{j}\geq a^{\ell}-b^{\ell}$ for any $1\leq\ell\leq j$ and $j\geq 1$ (the statement is trivial for $j=0$). Since $w_{\lambda}(k)\leq\tfrac{k+1}{N^{\lambda}}$ for all $k$, especially also if $k>N^{\lambda}-1$, we conclude that

[TABLE]

Besides, one computes analogously to above that $(k+1)^{j}-(k-1)^{j}\leq 2j(k+1)^{j-1}$, hence

[TABLE]

Finally, $w_{\lambda}(k)=1$ for $k>N^{\lambda}-1$, hence the above estimates imply

[TABLE]

For all other values of $k$, the differences yield zero. Thus, every element of $L_{w_{\lambda}^{j}}$ can be bounded, in the sense of operators, by the operator corresponding to the weight function

[TABLE]

Besides, since $l_{\lambda}^{j}(k)=0$ for $k>N^{\lambda}+1$, one obtains

[TABLE]

Inserting (47) to (49) into (44) and (45) and using that $\lambda\leq 1-d\beta$ implies $N^{\frac{d\beta+\lambda-1}{2}}\leq 1$, we conclude that

[TABLE]

Now we apply Grönwall’s inequality, for now on using the abbreviations $\alpha_{\psi,\varphi,s}(w_{\lambda}^{j};t)=:\alpha_{j}(t)$ and $I_{t}:=\int_{s}^{t}\lVert\varphi(s_{1})\rVert_{H^{k}(\mathbb{R}^{d})}^{2}\mathop{}\!\mathrm{d}s_{1}$. This yields

[TABLE]

where we have used that all integrands are non-negative and thus the upper boundary of all integrals could be replaced by $t$. Written explicitly, this gives

[TABLE]

with

[TABLE]

where we have estimated $I_{t}^{j}\mathrm{e}^{2j3^{j}I_{t}}<\mathrm{e}^{9^{j}I_{t}}$. To relate this estimate to $\lVert\widehat{m}^{j}\psi\rVert^{2}$, observe that for $0\leq k\leq N$,

[TABLE]

and

[TABLE]

Consequently, $m^{2j}(k)\leq N^{-j(1-\lambda)}w_{\lambda}^{j}(k)+w_{\lambda}^{b}(k)$, and we conclude

[TABLE]

for any $b\in\mathbb{N}$. Inserting these estimates into (51) yields

[TABLE]

To minimise the second term, we choose the maximal $\lambda=1-d\beta$, which concludes the proof of part (a).

The proof of part (b) is much simpler since we now consider the time evolution $\widetilde{U}_{\varphi}(t,s)$. The term corresponding to (41) vanishes, which implies that we may directly consider the weights $m^{2j}(k)$ instead of taking the detour via $w_{\lambda}^{j}(k)$. Analogously to (37), we define

[TABLE]

We will now abbreviate $\widetilde{U}_{\varphi}(t,s)\psi=:\widetilde{\psi}_{t}$. In this notation,

[TABLE]

We now evaluate this expression for the weight $m^{2j}(k)$, i.e.

[TABLE]

This corresponds to $w_{\lambda}^{j}(k)$ with the choice $\lambda=1$ in (46). Consequently, we define $l^{j}(k):=j3^{j}N^{-1}m^{2(j-1)}(k)$ analogously to (47) and conclude that $m^{2j}(k)-m^{2j}(k-2)\leq l^{j}(k)$ and $m^{2j}(k+2)-m^{2j}(k)\leq l^{j}(k)$. Analogously to the estimate of the first term in (43) and using the relation (48) for $\lambda=1$, we obtain

[TABLE]

The same Grönwall argument which led to (51) concludes the proof. ∎

Proof of Corollary 3.4. From Proposition 3.3a and the assumptions on the initial data, we conclude that for every $b\in\mathbb{N}$ and sufficiently large $N$,

[TABLE]

If $\gamma\geq 1-d\beta$, the leading order terms in both sums are the ones with maximal $n$, hence

[TABLE]

If one chooses $b>a\frac{1-d\beta}{1-2d\beta}$ for fixed $\beta<\frac{1}{2d}$, the second term is for sufficiently large $N$ dominated by the first one. For $\gamma<1-d\beta$, the leading order terms are those with $n=0$, hence

[TABLE]

which yields a non-trivial bound only for $\gamma>d\beta$. Part (b) follows analogously from part (b) of Proposition 3.3 without the restrictions on $\beta$ and $\gamma$ that are due to the second sum. ∎

4.3 Proof of Theorem 1

Proof of Lemma 3.5. We use the abbreviation $Z^{\beta}_{ij}=v^{\beta}_{ij}-\overline{v}^{\varphi(t)}(x_{i})-\overline{v}^{\varphi(t)}(x_{j})+2\mu^{\varphi(t)}$ as in (38), and drop all superscripts ${\varphi(t)}$ in $p^{\varphi(t)}$, $q^{\varphi(t)}$ and $\widehat{m^{\varphi(t)}}$ for simplicity. By Lemma 4.3a, $\mathcal{Q}^{\varphi(t)}\widehat{m}^{a}=\widehat{m}^{a}\mathcal{Q}^{\varphi(t)}$, hence

[TABLE]

where we have used that $\lVert Z^{\beta}_{ij}\rVert_{L^{\infty}(\mathbb{R}^{d})}\lesssim N^{d\beta}$ by Young’s inequality. Now observe that

[TABLE]

hence

[TABLE]

by (34), and analogously

[TABLE]

This implies part (a). For part (b), note that by Lemma 4.3a,

[TABLE]

Consequently,

[TABLE]

similarly to the estimate of $\lVert\widehat{m}^{a}\mathcal{Q}^{\varphi(t)}\psi\rVert$. The last inequality follows because by Lemma 4.1a, $\lVert p_{1}Z^{\beta}_{12}\rVert_{\mathrm{op}}^{2}\lesssim N^{d\beta}\lVert{\varphi(t)}\rVert^{2}_{H^{k}(\mathbb{R}^{d})}$ due to Young’s inequality and since $\lVert v^{\beta}\rVert_{L^{2}(\mathbb{R}^{d})}^{2}\lesssim N^{d\beta}$. Further, note that

[TABLE]

in the sense of operators. As in the estimate of $\mathcal{Q}^{\varphi(t)}$, we thus obtain for $\ell\in\{-1,1\}$

[TABLE]

and analogously $\lVert q_{1}q_{2}\widehat{m}^{a}_{\ell}\psi\rVert<4^{a}N\lVert\widehat{m}^{a+3}\psi\rVert^{2}$ and $\lVert q_{1}q_{2}q_{3}\widehat{m}^{a}_{\ell}\psi\rVert<4^{a}\lVert\widehat{m}^{a+3}\psi\rVert^{2}$. Together, this implies part (b). ∎

Proof of Theorem 1. Let $a\in\mathbb{N}_{0}$ such that $6a\leq A$. Recall that by Definition 3.1,

[TABLE]

for any $a\geq 0$, where $T_{n}^{(k)}$ is given by

[TABLE]

where

[TABLE]

with $I^{\varphi(t)}_{1}=\mathcal{C}^{\varphi(t)}$ and $I^{\varphi(t)}_{2}=\mathcal{Q}^{\varphi(t)}$ and $(j_{1},\makebox[10.00002pt][c]{.\hfil.\hfil.},j_{n})\in\mathcal{S}_{n}^{(k)}$. In this notation,

[TABLE]

hence the Duhamel expansion (33) of $\psi(t)$ reads

[TABLE]

Here, $\widetilde{T}_{n}^{(k)}$ is obtained from $T_{n}^{(k)}$ by replacing the first $\widetilde{U}_{\varphi}(t,s_{n})$ by the full time evolution $U(t,s_{n})$, i.e., for $n<k<2n$,

[TABLE]

Consequently,

[TABLE]

since the first double sum contributes only if $2n\geq a+1$, and in this case $\min\{2n,a\}=a$. Note that for $k=n$, $j_{1}=\dots=j_{k}=1$, hence $T_{k}^{(k)}$ and $\widetilde{T}_{k}^{(k)}$ exclusively contain $\mathcal{C}^{\varphi(s_{l})}$. Using Duhamel’s formula, the last expression can thus be expanded as

[TABLE]

By unitarity of $U(t,s)$ and $\widetilde{U}_{\varphi}(t,s)$,

[TABLE]

With this, (52) and (53) imply for $a=0,1$

[TABLE]

which coincides with (27) and (32). For $a\geq 2$, we find

[TABLE]

where we used that $a+2\leq 2a$ for $a\geq 2$. To estimate $\lVert t_{(j_{1},\makebox[7.00002pt][c]{.\hfil.\hfil.},j_{n})}^{(k)}\rVert^{2}$ for $a+1\leq k\leq 2a$ and $n\leq k$, note first that Lemma 3.5 and Proposition 3.3b can be combined into the single statement

[TABLE]

for $j\in\{1,2\}$ and any $\psi\in L^{2}_{\mathrm{sym}}(\mathbb{R}^{dN})$. Hence, with \delta_{\mu}:=2(n-\mu+1)+(j_{n}+j_{n-1}+\makebox[10.00002pt][c]{\cdot$\hfil$\cdot$\hfil$\cdot}+j_{\mu}) and $\eta_{\mu}:=\prod_{\ell=0}^{\mu}\lVert\varphi(s_{n-\ell})\rVert^{2}_{H^{k}(\mathbb{R}^{d})}$, we obtain for $n\leq k$

[TABLE]

Since $j_{1}+,\makebox[10.00002pt][c]{.\hfil.\hfil.},+j_{n}=k$ and $n\leq k\leq 2a$, we find $\delta_{1}=2n+k\leq 3k\leq 6a\leq A$, hence assumption A3 yields

[TABLE]

Let us for the moment focus on the $N$-dependent factors in (58), thereby neglecting all other contributions to the sum. This yields

[TABLE]

For $\gamma\geq 1-d\beta$, the leading order term in the sum $\sum_{\nu_{n}}$ is the term corresponding to the choice \nu_{n}=\delta_{1}-(\nu_{1}+\makebox[10.00002pt][c]{\cdot$\hfil$\cdot$\hfil$\cdot}+\nu_{n-1})=2n+k-(\nu_{1}+\makebox[10.00002pt][c]{\cdot$\hfil$\cdot$\hfil$\cdot}+\nu_{n-1}), which yields the total factor $N^{k(-1+d\beta)}N^{d\beta\delta_{1}}=N^{-k+2d\beta(n+k)}$. This factor is maximal for $n=k$. For $\gamma<1-d\beta$, the leading term corresponds to the choice $\nu_{1}=\dots=\nu_{n}=0$, which yields $N^{2n(1-\gamma)+k(d\beta-\gamma)}$. Also here, the maximal contribution issues from $n=k$. In fact, the leading contributions for both ranges of $\gamma$ can be summarised as $N^{-k\delta(\beta,\gamma)}$, where

[TABLE]

as defined in (31). Hence, for sufficiently large $N$, the dominating terms is the one with $n=k$, which comes from $t^{(k)}_{(j_{1},\makebox[7.00002pt][c]{.\hfil.\hfil.},j_{k})}=t^{(k)}_{(1,\makebox[7.00002pt][c]{.\hfil.\hfil.},1)}$.

[TABLE]

and (54) to (56) can be summarised as

[TABLE]

It remains to evaluate the estimate (58) for $n=k$. In this case, $j_{1}=\dots=j_{k}=1$ and $\delta_{\mu}=3(k-\mu+1)$. Note also that the constants $C_{a}^{\,t}$ are increasing in $a$ and $t$, hence C_{\delta_{k-\mu}-(\nu_{1}+\makebox[7.00002pt][c]{\cdot$\hfil$\cdot$\hfil$\cdot}+\nu_{\mu-1})}^{\,s_{k-\mu}-s_{k-\mu-1}}\leq C_{3(\mu+1)}^{\,s_{k-\mu}}. Further, observe that \delta_{k}+\makebox[10.00002pt][c]{\cdot$\hfil$\cdot$\hfil$\cdot}+\delta_{2}=\tfrac{3}{2}k(k-1)\leq\tfrac{3}{2}k^{2}. Consequently,

[TABLE]

where we have used that each sum $\sum_{\nu_{\mu}}$ in (58) contains at most $\delta_{k-\mu+1}+1=3\mu+1$ addends, and that the prefactor of the leading order term for $\gamma\geq 1-d\beta$ is $\mathfrak{C}_{\,0}=1$, whereas it is $\mathfrak{C}_{\,3k}$ for $\gamma<1-d\beta$. Consequently, for sufficiently large $N$, the maximum in (59) is attained for $k=a+1$. Inserting the explicit formula $C_{j}^{\,t,s}=j!\,3^{j(j+1)}\mathrm{e}^{9^{j}I_{t}}$ with $I_{t}=\int_{s}^{t}\lVert\varphi(s_{1})\rVert^{2}_{H^{k}(\mathbb{R}^{d})}\mathop{}\!\mathrm{d}s_{1}$ yields

[TABLE]

where we have bounded all $a$-dependent, time-independent expressions by a constant $c\lesssim 1$. ∎

Acknowledgements

L.B. gratefully acknowledges the support by the German Research Foundation (DFG) within the Research Training Group 1838 “Spectral Theory and Dynamics of Quantum Systems”, and wishes to thank Stefan Teufel, Sören Petrat and Marcello Porta for helpful discussions. This project has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie Grant Agreement No. 754411. N.P. gratefully acknowledges support from NSF grant DMS-1516228 and DMS-1840314. P.P.’s research was funded by DFG Grant no. PI 1114/3-1. Part of this work was done when N.P. and P.P. were visiting CCNU, Wuhan. N.P. and P.P. thank A.S. for his hospitality at CCNU.