Hostname: page-component-76fb5796d-45l2p Total loading time: 0 Render date: 2024-04-29T17:28:46.127Z Has data issue: false hasContentIssue false

Bogoliubov theory in the Gross-Pitaevskii limit: a simplified approach

Published online by Cambridge University Press:  21 October 2022

Christian Hainzl
Affiliation:
Department of Mathematics, LMU Munich, Theresienstrasse 39, 80333 Munich, Germany; E-mail: hainzl@math.lmu.de
Benjamin Schlein
Affiliation:
Institute of Mathematics, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland; E-mail: benjamin.schlein@math.uzh.ch
Arnaud Triay
Affiliation:
Department of Mathematics, LMU Munich, Theresienstrasse 39, 80333 Munich, Germany; E-mail: triay@math.lmu.de

Abstract

We show that Bogoliubov theory correctly predicts the low-energy spectral properties of Bose gases in the Gross-Pitaevskii regime. We recover recent results from [6, 7]. While our main strategy is similar to the one developed in [6, 7], we combine it with new ideas, taken in part from [15, 25]; this makes our proof substantially simpler and shorter. As an important step towards the proof of Bogoliubov theory, we show that low-energy states exhibit complete Bose-Einstein condensation with optimal control over the number of orthogonal excitations.

Type
Mathematical Physics
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press

1 Introduction

We consider a Bose gas consisting of N particles moving in the box $\Lambda = [-1/2; 1/2]^3$ with periodic boundary conditions. In the Gross-Pitaevskii regime, particles interact through a potential with scattering length of the order $1/N$ . The Hamilton operator acts on the Hilbert space $L^2_s (\Lambda ^N)$ of permutation symmetric complex-valued square integrable functions on $\Lambda ^N$ , and it has the form

(1) $$ \begin{align} H_N = \sum_{i=1}^N -\Delta_i + \sum_{i<j}^N V_N(x_i-x_j) \end{align} $$

where

$$ \begin{align*} V_N(x) := N^{2} V(N x), \end{align*} $$

for a $V\in L^{2}(\mathbb {R}^3)$ nonnegative, radial and compactly supported. We denote the scattering length of V by $\mathfrak {a}> 0$ . Following [Reference Hainzl15, Reference Hainzl and Seiringer16], we define it through the formula

(2) $$ \begin{align} 4\pi \mathfrak{a} = \frac{1}{2} \int_{\mathbb{R}^3} V(x) \mathrm{d}x - \left \langle \frac 12 V, \frac{1} {-\Delta + \frac 1 2 V} \frac 1 2 V \right \rangle. \end{align} $$

As first proven in [Reference Lieb and Yngvason22, Reference Lieb, Seiringer and Yngvason20], the ground state energy $E_N$ of (1) satisfies

(3) $$ \begin{align} E_N / N \to 4\pi \mathfrak{a} \end{align} $$

in the limit $N \to \infty $ . In particular, to leading order, the ground state energy only depends on the interaction potential through its scattering length $\mathfrak {a}$ . In [Reference Lieb and Seiringer18, Reference Lieb and Seiringer21, Reference Nam, Rougerie and Seiringer24], it was also shown that the ground state of the Hamiltonian $H_N$ in equation (1) and, in fact, every normalised sequence $\psi _N \in L_s^2 (\Lambda ^N)$ of approximate ground states with

$$\begin{align*}\frac{1}{N} \, \langle \psi_N , H_N \psi_N \rangle \to 4\pi \mathfrak{a} , \end{align*}$$

exhibits complete Bose-Einstein condensation in the zero-momentum state $\phi _0 (x) = 1$ , for all $x \in \Lambda $ , in the sense that the corresponding one-particle reduced density matrix $\gamma _N$ (normalised so that $\operatorname {\mathrm {Tr}} \, \gamma _N = 1$ ) satisfies

$$\begin{align*}\lim_{N \to \infty} \langle \phi_0, \gamma_N \phi_0 \rangle = 1 \,. \end{align*}$$

Recently, a rigorous version of Bogoliubov theory [Reference Bogoliubov8] has been developed in [Reference Boccato, Brennecke, Cenatiempo and Schlein4, Reference Boccato, Brennecke, Cenatiempo and Schlein5, Reference Boccato, Brennecke, Cenatiempo and Schlein6, Reference Boccato, Brennecke, Cenatiempo and Schlein7] to provide more precise information about the low-energy spectrum of $H_N$ in equation (1), resolving the ground state energy and low-lying excitations up to errors that vanish in the limit $N \to \infty $ ; and about the corresponding eigenvectors, showing Bose-Einstein condensation with optimal control over the number of orthogonal excitations. Analogous results have also been established for Bose gases trapped by external potentials in the Gross-Pitaevskii regime [Reference Brennecke, Schlein and Schraven10, Reference Brennecke, Schlein and Schraven11, Reference Nam, Napiórkowski, Ricaud and Triay23, Reference Nam and Triay25] and for Bose gases in scaling limits interpolating between the Gross-Pitaevskii regime and the thermodynamic limit [Reference Adhikari, Brennecke and Schlein1, Reference Brennecke, Caporaletti and Schlein9]. Very recently, the upper bound for the ground state energy has also been extended to the case of hard-sphere interaction, as announced in [Reference Basti, Cenatiempo, Olgiati, Pasqualetti and Schlein2].

In this paper, we propose a new and substantially simpler proof of the results established in [Reference Boccato, Brennecke, Cenatiempo and Schlein6, Reference Boccato, Brennecke, Cenatiempo and Schlein7]. Our approach follows some of the ideas in the proof of Bose-Einstein condensation with optimal bounds on the number of excitations obtained in [Reference Hainzl15]. Moreover, it makes use of some ideas introduced in [Reference Nam and Triay25] for the case of particles trapped by an external potential. The next theorem is our main result; it describes the low-energy spectrum of $H_N$ in equation (1).

Theorem 1. Let $V \in L^{2} (\mathbb {R}^3)$ be nonnegative, radial and compactly supported, and let $E_N$ denote the ground state energy of $H_N$ in equation (1). Then the spectrum of $H_N - E_N$ below a threshold $\Theta \leq N^{1/17}$ consists of eigenvalues having the form

(4) $$ \begin{align} \sum_{p \in 2\pi \mathbb{Z}^3 \backslash \{ 0 \}} n_p \sqrt{|p|^4 + 16 \pi \mathfrak{a} p^2} + {\mathcal O} (N^{-1/17} \Theta) \end{align} $$

with $n_p \in \mathbb {N}$ , for all $p \in 2\pi \mathbb {Z}^3 \backslash \{ 0 \}$ .

Remark. Our analysis also provides a precise estimate for the ground state energy $E_N$ of the Hamiltonian $H_N$ in equation (1), showing that

(5) $$ \begin{align} E_N = 4\pi \mathfrak{a}_N (N-1) + \frac{1}{2} \sum_{p \in 2\pi\mathbb{Z}^3\backslash \{ 0 \}} \Big[ \sqrt{|p|^4 + 16 \pi \mathfrak{a} p^2} - p^2 - 8\pi \mathfrak{a} + \frac{(8\pi \mathfrak{a})^2}{2p^2} \Big] + {\mathcal O} (N^{-1/17}) \end{align} $$

with a ‘box scattering length’ $\mathfrak {a}_N$ (defined in the next section) satisfying $|\mathfrak {a}_N - \mathfrak {a}| \lesssim N^{-1}$ . This immediately implies that $E_N$ is given by equation (5), with $\mathfrak {a}_N$ replaced by the true scattering length $\mathfrak {a}$ , up to an error that remains bounded as $N \to \infty $ . In [Reference Boccato, Brennecke, Cenatiempo and Schlein6], the order-one correction arising from $N (\mathfrak {a}_N - \mathfrak {a})$ was also computed. Here we skip this step to keep our presentation as simple as possible. Note that an estimate similar to equation (5) has recently been shown to hold in the thermodynamics limit; see [Reference Fournais and Solovej13, Reference Fournais and Solovej14] for the lower bound and [Reference Yau and Yin27, Reference Basti, Cenatiempo and Schlein3] for the upper bound.

The main strategy we use to prove Theorem 1 is similar to the one developed in [Reference Boccato, Brennecke, Cenatiempo and Schlein6]. First we switch to the formalism of second quantisation, expressing the Hamilton operator in momentum space in terms of creation and annihilation operators. Then we renormalise the Hamilton operator, conjugating it first with a generalised Bogoliubov transformation (the exponential of a quadratic expression in the modified creation and annihilation operators $b_p^{\dagger } = a_p^{\dagger } a_0 / \sqrt {N}$ , $b_p = a_0^{\dagger } a_p / \sqrt {N}$ ) and afterwards with the exponential of a cubic expression in (modified) creation and annihilation operators. Effectively, these conjugations regularise the interaction potential. As a last step, we diagonalise the resulting quadratic Hamiltonian; this allows us to establish Bose-Einstein condensation with optimal bounds on the number of excitations and to compute the low-energy spectrum, proving the estimate (4).

Compared with [Reference Boccato, Brennecke, Cenatiempo and Schlein6], our approach has the following advantages. First, we make a different choice for the coefficients ${\varphi }_p$ of the quadratic and cubic transformations used to renormalise the Hamiltonian, which should model correlations among particles. Instead of the ground state of a Neumann problem on a ball of radius $\ell>0$ , we consider here the solution of an appropriate zero-energy scattering equation, describing scattering processes inside the box $\Lambda $ . This simplifies the proof of important properties of ${\varphi }$ and improves cancellations between different terms arising in the many-body analysis. Second, we restrict the quadratic conjugation to momenta $|p|> N^{\alpha }$ for some $0 < \alpha <1$ . Consequently, it is enough to expand its action to first or, in a few cases, second order; higher-order contributions are negligible. This is a substantial advantage compared with [Reference Boccato, Brennecke, Cenatiempo and Schlein6], where no cutoff was imposed and all contributions had to be computed precisely (in contrast to standard Bogoliubov transformations, the action of generalised Bogoliubov transformations is not explicit). The presence of the cutoff means the interaction is regularised only up to length scales $\ell \leq N^{-\alpha }$ ; this needs to be compensated at the end when we diagonalise the quadratic Hamiltonian resulting from the renormalisation procedure. Another important simplification of the analysis concerns the final diagonalisation. As in [Reference Boccato, Brennecke, Cenatiempo and Schlein6], we implement it through a generalised Bogoliubov transformation defined (like the first quadratic transformation) in terms of the modified creation and annihilation operators $b_p^{\dagger }, b_p$ . Here, however, instead of expanding the action of the generalised Bogoliubov transformation to all orders, we compare it directly with the explicit action of the corresponding standard Bogoliubov transformation, using an appropriate interpolation. Finally, we use the tool of localisation in the number of particles not only to show Bose-Einstein condensation (similarly to [Reference Boccato, Brennecke, Cenatiempo and Schlein7]) but also to compute the spectrum and prove Theorem 1. This makes the analysis substantially simpler (but provides a worse estimate of the error).

2 Fock space formalism

We introduce the bosonic Fock space

$$\begin{align*}\mathcal{F} = \bigoplus_{n \geq 0} L^2_s (\Lambda^n). \end{align*}$$

For a momentum $p \in \Lambda ^{\ast } = 2\pi \mathbb {Z}^3$ and denoting $u_p(x) = e^{i p \cdot x}$ , we define $a_p^{\dagger } = a^{\dagger }(u_p)$ and $a_p = a(u_p)$ , where $a^{\dagger }$ and a are the usual creation and annihilation operators. They satisfy the canonical commutation relations

(6) $$ \begin{align} \big[ a_p , a_q^{\dagger} \big] = \delta_{p,q}, \qquad \big[ a_p, a_q \big] = \big[ a_p^{\dagger} , a_q^{\dagger} \big] = 0. \end{align} $$

We denote, in configuration space, the creation and annihilation operator-valued distributions by ${\check a}^{\dagger }_x, {\check a}_x$ ; they satisfy $a^{\dagger }_p = \int e^{ip\cdot x} {\check a}^{\dagger }_x \mathrm {d}x$ , $a_p = \int e^{-ip\cdot x} {\check a}_x \mathrm {d}x$ . The number-of-particles operator $\mathcal {N}$ on $\mathcal {F}$ is given by

$$\begin{align*}\mathcal{N} = \sum_{p \in \Lambda^{\ast}} a_p^{\dagger} a_p. \end{align*}$$

In the formalism of second quantisation, the Hamilton operator in equation (1) takes the form

(7) $$ \begin{align} H_N = \sum_{p \in \Lambda^{\ast}} p^2 a_p^{\dagger} a_p + \frac{1}{2}\sum_{r,p,q \in \Lambda^{\ast}} \hat{V}_N (r) a_{p+r}^{\dagger} a_q^{\dagger} a_p a_{q+r}, \end{align} $$

with

(8) $$ \begin{align} \hat{V}_N (r) = \frac 1{N} \hat V(r/N). \end{align} $$

To recover the expression in equation (1), we have to restrict equation (7) to the sector with $\mathcal {N} = N$ .

Because of the presence of Bose-Einstein condensation, the mode with $p=0$ plays a special role when considering states with low energy. We introduce the notation $\mathcal {N}_0 = a_0^{\dagger } a_0$ and $\mathcal {N}_+ = \mathcal {N} - \mathcal {N}_0$ for the operators measuring the number of particles in the condensate and the number of excitations, respectively. Following Bogoliubov [Reference Bogoliubov8], we decompose equation (7) according to the number of $a_0, a_0^{\dagger }$ operators. Since (on $\{ \mathcal {N} = N \}$ )

$$ \begin{align*} a^{\dagger}_0 a^{\dagger}_0 a_0 a_0 = \mathcal{N}_0(\mathcal{N}_0 - 1) = (N - \mathcal{N}_+)(N - \mathcal{N}_+ -1) = N(N - 1) - \mathcal{N}_+(2N - 1) + \mathcal{N}_+^2, \end{align*} $$

we can rewrite equation (7) as

(9) $$ \begin{align} H_N = H_0+ H_1 + H_2 + Q_2 + Q_3 + Q_4, \end{align} $$

where

$$\begin{align*}\begin{aligned} H_0 &= \frac{\hat{V}_N(0) }{2} N (N - 1 ) , \qquad H_1 = \sum_{p \neq 0} p^2 a_p^{\dagger} a_p, \\[3pt] H_2 &= \sum_{p\neq 0} \hat V_N(p) a^{\dagger}_p a_p(N - \mathcal{N}_+) - \frac {\hat V_N(0) }{2} \mathcal{N}_+(\mathcal{N}_+ -1). \end{aligned} \end{align*}$$

and

(10) $$ \begin{align} \begin{aligned} Q_2 = &\frac{1}{2}\sum_{p \neq 0} \hat V_N(p) [a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0 + \text{h.c.}], \\[3pt] Q_3 = &\sum_{q,r,q+r\neq 0} \hat V_N(r)\left[a_{q+r}^{\dagger} a_{-r}^{\dagger} a_{q}a_0 + \text{h.c.} \right], \\[3pt] Q_4 = &\frac{1}{2}\sum_{p,q \neq 0, r \neq -p, r\neq -q} \hat V_N(r) a_{p+r}^{\dagger} a_q^{\dagger} a_p a_{q+r}. \end{aligned} \end{align} $$

Since we isolated the contributions of the zero modes, from now on we follow the convention that the indices appearing in creation and annihilation operators are always nonzero except when stated otherwise.

Naive power counting, based on the fact that $a_0, a_0^{\dagger } \simeq \sqrt {N}$ due to the presence of Bose-Einstein condensation and on the scaling (8) of the interaction, suggests that the terms $Q_3$ and $Q_4$ are small. For this reason, Bogoliubov neglected these contributions and diagonalised the remaining quadratic terms. This led to expressions similar to equations (4) and (5) for the low-energy spectrum of $H_N$ , but with the scattering length replaced by its first and second Born approximations. In fact, because of the slow decay of the potential in Fourier space, the operators $Q_3$ and $Q_4$ are not small. They instead contain important terms that effectively renormalise the interaction and produce the scattering length appearing in the formulas in equations (4) and (5). To obtain a rigorous proof of Theorem 1, it is therefore crucial that we first extract the large contributions to the energy hidden in the cubic and quartic operators $Q_3$ , $Q_4$ ; only afterwards can we diagonalise the remaining quadratic terms.

Let us give a little more detail about the main ideas of the proof. Following the strategy of [Reference Boccato, Brennecke, Cenatiempo and Schlein6], we will first conjugate equation (9) with a unitary operator of the form $e^{\mathcal {B}_2}$ , where $\mathcal {B}_2$ is a quadratic expression in creation and annihilation operators $a_p, a_p^{\dagger }$ , associated with momenta $p \not = 0$ . The goal of this conjugation is to extract contributions that regularise the off-diagonal term $Q_2$ and, at the same time, reconstruct the leading-order ground state energy $4 \pi {\mathfrak a}_N N$ when combined with $H_0$ . Roughly speaking, neglecting several error terms, we will find

(11) $$ \begin{align} \begin{aligned} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} \simeq & \, \, 4 \pi \mathfrak{a}_N (N-1) + \sum_{|p|\leq N^{\alpha}} \frac{(4\pi \mathfrak{a}_N)^2}{p^2} \\[3pt] &+ \sum_{p \in \Lambda^*} (p^2 + 2\hat V(0) - 8\pi \mathfrak{a}_N) a^{\dagger}_p a_p + \sum_{|p| \leq N^{\alpha}} 4\pi \mathfrak{a}_N[ a_p^{\dagger} a_{-p}^{\dagger} + \text{h.c.} ] + Q_3 + Q_4. \end{aligned} \end{align} $$

As explained in the introduction, an important difference, compared with [Reference Boccato, Brennecke, Cenatiempo and Schlein6], is that here we impose an infrared cutoff in $\mathcal {B}_2$ , defined in equation (13), letting it act only on momenta $|p|> N^{\alpha }$ . On one hand, this choice simplifies the computation of the action of $\mathcal {B}_2$ (it allows us to expand it; important contributions arise only from the first and second commutators). On the other hand, it produces terms, like the sum on the first line and the regularised off-diagonal quadratic term on the second line of equation (11), which contribute to the energy to order $N^{\alpha }$ ; these terms are larger than the precision we are looking for and will need to be compensated for with the second quadratic transformation. Notice that the idea of using an infrared cutoff in the quadratic conjugation already appeared in the proof of complete Bose-Einstein condensation given in [Reference Boccato, Brennecke, Cenatiempo and Schlein4] and, more recently, in the proof of the validity of Bogoliubov theory for Bose gases trapped by an external potential obtained in [Reference Nam and Triay25].

Observing equation (11), it is clear that we still have to renormalise the diagonal quadratic term (proportional to $\hat {V} (0)$ ) and the cubic term $Q_3$ . To this end, we will introduce a unitary transformation $e^{\mathcal {B}_3}$ , with $\mathcal {B}_3$ , defined in equation (59), cubic in the operator $a_p, a_p^{\dagger }$ , with $p \not = 0$ . Up to several negligible errors, conjugation with $e^{\mathcal {B}_3}$ will lead us to

(12) $$ \begin{align} \begin{aligned} e^{-\mathcal{B}_3} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2}e^{\mathcal{B}_3} \simeq \, & \, \, 4 \pi \mathfrak{a}_N (N-1) + \sum_{|p|\leq N^{\alpha}} \frac{(4\pi \mathfrak{a}_N)^2}{p^2}. \\ &+ \sum_{p \in \Lambda^*} (p^2 + 8\pi \mathfrak{a}_N) a^{\dagger}_p a_p + \sum_{|p| \leq N^{\alpha}} 4\pi \mathfrak{a}_N[ a_p^{\dagger} a_{-p}^{\dagger} + \text{h.c.} ] + Q_4 \,. \end{aligned} \end{align} $$

The only term on the right-hand side of the last equation where we still have the original, singular, potential $\hat {V}_N$ is $Q_4$ ; all other terms have been renormalised and are now expressed in terms of the scattering length $\frak {a}_N$ . Fortunately, $Q_4$ is positive; for this reason, we do not need to renormalise it (for lower bounds, it can be neglected; for upper bounds, it only needs to be controlled on special trial states). Finally, in section 5, we will apply a second quadratic transformation $e^{\mathcal {B}_4}$ to diagonalise the remaining quadratic Hamiltonian on the right-hand side of equation (12). This will lead us to

$$ \begin{align*} e^{-\mathcal{B}_4}e^{-\mathcal{B}_3} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2}e^{\mathcal{B}_3} e^{\mathcal{B}_4} \simeq &\; 4 \pi \mathfrak{a}_N (N-1) + \frac{1}{2}\sum_{p} \left[ \sqrt{p^4 + 16 \pi \mathfrak{a}_N p^2} - p^2 - 8 \pi \mathfrak{a}_N + \frac{(8\pi \mathfrak{a}_N)^2}{2p^2}\right] \\ &+ \sum_{p} \sqrt{p^4 + p^2 16 \pi \mathfrak{a}_N } \, a^{\dagger}_p a_p + Q_4, \end{align*} $$

which will allow us to show Theorem 1. To control error terms, we use the tool of localisation in the number of particles to show Bose-Einstein condensation (similarly to [Reference Boccato, Brennecke, Cenatiempo and Schlein7, Reference Nam and Triay25]).

3 Quadratic renormalisation

Starting with the quadratic transformation, we conjugate the Hamiltonian $H_N$ in equation (7) with the unitary $e^{\mathcal {B}_2}$ , where

(13) $$ \begin{align} {\mathcal{B}}_2 = \frac 1{2} \sum_{p} \tilde {\varphi}_p [ a^{\dagger}_p a^{\dagger}_{-p} a_0 a_0 - \textrm{h.c.}]. \end{align} $$

We are going to fix the coefficients $\tilde {\varphi }_p$ so that the commutator $[ H_1 + Q_4 , \mathcal {B}_2 ]$ arising from the action of the unitary $e^{\mathcal{B}_2}$ (13) renormalises the off-diagonal quadratic term $Q_2$ (effectively replacing the singular potential $V_N$ with a regularised interaction having the same scattering length). To this end, we choose ${\varphi }_p$ satisfying the relations

(14) $$ \begin{align} p^2{\varphi}_p + \frac{1}{2}\sum_{q\neq 0} \hat{V}_N \left((p-q)\right){\varphi}_q = -\frac{1}{2} \hat{V}_N(p)\, , \end{align} $$

for all $p \in \Lambda ^{\ast }_+ = \Lambda ^{\ast } \backslash \{ 0 \}$ . Equation (14) is a truncated version of the zero-energy scattering equation for the potential $V_N$ on the whole space $\mathbb {R}^{3}$ .

To prove the existence of a solution for equation (14), we consider the operator

$$\begin{align*}\mathfrak{h} = -\Delta + \frac{1}{2} V_N \end{align*}$$

acting on the one-particle space $L^2 (\Lambda )$ (for N large enough, $V_N$ is supported in $[-1/2; 1/2]^3$ and can be periodically extended to define a function on the torus). Denoting by $P_0^{\perp }$ the orthogonal projection onto the orthogonal complement of the zero-momentum mode ${\varphi }_0$ in $L^2 (\Lambda )$ , we find (since $V_N \geq 0$ ) that $P_0^{\perp } \mathfrak {h} P_0^{\perp } \geq C> 0$ , and therefore that $P_0^{\perp } \mathfrak {h} P_0^{\perp }$ is invertible. Thus, we can define $\check {\varphi } \in L^2 (\Lambda )$ through

(15) $$ \begin{align} \check{\varphi} = -\frac{1}{2} P_0^{\perp} \left[ P_0^{\perp} \Big(-\Delta + \frac{1}{2} V_N \Big) P_0^{\perp} \right]^{-1} P_0^{\perp} V_N. \end{align} $$

It is then easy to check that the Fourier coefficients of $\check {{\varphi }}$ satisfy the relations in equation (14).

Using the sequence $\{ {\varphi }_p \}_{p \in 2\pi \mathbb {Z}^3 \backslash \{ 0 \}}$ , we can define the ‘box scattering length’ of $V_N$ by

(16) $$ \begin{align} 8 \pi \mathfrak{a}_N := N \Big[ \hat V_N (0) + \sum_p \hat{V}_N (p) {\varphi}_p \Big] = \hat{V} (0) + N \sum_p \hat{V}_N (p) {\varphi}_p. \end{align} $$

As proven in [Reference Hainzl15], we have that $|\mathfrak {a}_N - \mathfrak {a}| \lesssim N^{-1}$ .

As explained earlier, we first renormalise the high-momenta part of $Q_2$ ; for this reason, we use a cutoff version of ${\varphi }_p$ to momenta $|p|> N^{\alpha }$ for some $0< \alpha < 1$ . We therefore define

(17) $$ \begin{align} \tilde {\varphi}_p = {\varphi}_p \chi_{|p|> N^{\alpha}}. \end{align} $$

The next lemma lists some important properties of the sequences ${\varphi }, \tilde {{\varphi }}$ and the scattering length $\mathfrak {a}_N$ that will be useful for our analysis.

Lemma 2. Let $V \in L^{2} (\mathbb {R}^3)$ be nonnegative and compactly supported. Define $\check {\varphi }$ as in equation (15), and denote by ${\varphi }_p$ the corresponding Fourier coefficients. Then ${\varphi }_p \in \mathbb {R}$ , ${\varphi }_{-p} = {\varphi }_p$ and

(18) $$ \begin{align} |{\varphi}_p| \lesssim \frac{1}{N p^2} \end{align} $$

for all $p \in 2\pi \mathbb {Z} \backslash \{ 0 \}$ . Moreover, with equation (17), we have

$$ \begin{align*} & \|\tilde {\varphi} \|_2 \lesssim N^{ - 1 - \alpha/2},\qquad \|\tilde {\varphi}\|_{\infty} \lesssim N^{- 1 - 2\alpha}, \qquad \|\tilde {\varphi} \|_1 \lesssim 1, \end{align*} $$

and

(19) $$ \begin{align} N \sum_{p} \hat{V}_N(p) \tilde {\varphi}_p &= 8\pi \mathfrak{a}_N - \hat V(0) + {\mathcal O}(N^{\alpha-1}). \end{align} $$

Proof. Multiplying equation (14) by ${\varphi }_p$ , summing over p and using that $V \geq 0$ , we obtain

(20) $$ \begin{align} 2 \|p{\varphi}\|_2^2 = - \sum_{p} \hat V_N(p) {\varphi}_p - \sum_{p,q} \hat V_N(p-q) {\varphi}_p {\varphi}_q \leq - \sum_{p} \hat V_N(p) {\varphi}_p. \end{align} $$

On one hand, this implies that $\| p {\varphi } \|_2 \lesssim \| V_N \|_2 \| {\varphi } \|_2 < \infty $ (the last bound is not uniform in N; it follows from equation (15)). On the other hand, the estimate (20) leads to

$$\begin{align*}2 \| p {\varphi} \|^2_2 \leq \| V_N / |p| \|_2 \| p {\varphi} \|_2. \end{align*}$$

Dividing by $\| p {\varphi } \|_2$ and squaring, we obtain

(21) $$ \begin{align} \|p{\varphi}\|_2^2 \lesssim \sum_p \frac{|\hat V_N(p)|^2}{p^2} \lesssim \| \hat{V}_N \|_{\infty}^2 \| |p|^{-2} \chi_{|p| < N} \|_1 + \| \hat{V}_N \|_{\infty} \| \hat{V}_N \|_2 \| |p|^{-2} \chi_{|p|> N} \|_2 \lesssim N^{-1}. \end{align} $$

Using equation (14) again, we obtain the pointwise bound

(22) $$ \begin{align} | p^2 {\varphi}_p| \leq |\hat{V}_N (p)| + \Big[ \sum_q \frac{|\hat{V}_N (p-q)|^2}{q^2} \Big]^{1/2} \| q {\varphi} \|_2 \lesssim N^{-1}, \end{align} $$

where we proceeded as in (21) to bound $\| |\hat {V}_N|^2 * |q|^{-2} \|_{\infty }$ . This proves the bound (18) and immediately implies the bounds for $\| \tilde {\varphi } \|_2 , \| \tilde {\varphi } \|_{\infty }$ . To obtain the bound on $ \| \tilde {\varphi } \|_1$ , we divide equation (14) by $|p|^2$ . Proceeding as in (21), we obtain

$$ \begin{align*} \sum_{p\neq 0} \frac{|\hat{V}_N (p)|}{|p|^2} \lesssim 1 \, , \end{align*} $$

and hence we only have to bound $\||p|^{-2} (\hat {V}_N \ast {\varphi }) \|_{1}$ . Iterating equation (14) and using the regularising estimate $\||p|^{-2} \hat {V}_N \ast g\|_{6p/(6+p) + \varepsilon } \leq C_{\varepsilon } \|\hat {V}_N\|_2 \|g\|_{p}$ for all $\varepsilon>0$ , $p \geq 6/5$ , $g \in \ell ^{p}(\Lambda ^*)$ and some $C_{\varepsilon }> 0$ , we obtain that $\|{\varphi }\|_1 < \infty $ . Separating high and low momenta, we obtain for $A \geq 1$ and $\varepsilon>0$

$$ \begin{align*} \|\varphi\|_{1} &\lesssim 1+ \| \chi_{|p|> A N} |p|^{-2} \|_{2} \|V_N \check {\varphi}\|_{2} + \|\chi_{|p| \leq A N}|p|^{-2}\|_{1} \| \hat{V}_N \ast {\varphi} \|_{\infty} \\ &\lesssim 1 + A^{-\tfrac{1}{2}} \|{\varphi}\|_{1} + A \, , \end{align*} $$

where we used that $\|\check {\varphi }\|_{\infty } \leq \|{\varphi }\|_1$ and the Hölder inequality as in (22) to estimate $\| \hat {V}_N \ast {\varphi } \|_{\infty }$ . Taking A sufficiently large but fixed, we obtain $\|\varphi \|_{1} \lesssim 1$ .

The estimate (19) follows by noticing that, from the definition given in equation (16),

$$\begin{align*}\Big| 8\pi \mathfrak{a}_N - \hat{V} (0) - N \sum_p V_N (p) \tilde{\varphi}_p \Big| \leq N \sum_{|p| \leq N^{\alpha}} |\hat{V}_N (0)| |{\varphi}_p| \lesssim \frac{1}{N} \sum_{|p| \leq N^{\alpha}}\frac{1}{|p|^2} \lesssim N^{-1+\alpha}, \end{align*}$$

where we use equation (18) and $\| \hat {V}_N \|_{\infty } \lesssim N^{-1}$ .

Using the bounds in Lemma 2, we can control the growth of the number of excitations w.r.t. the action of $\mathcal {B}_2$ ; the proof of the next lemma can be found, for example, in [Reference Brennecke and Schlein12, Lemma 3.1].

Lemma 3. For every $n \in \mathbb {N}$ and $|s| \leq 1$ , we have

$$ \begin{align*} \pm (e^{-s{\mathcal{B}}_2}\mathcal{N}_+ e^{s{\mathcal{B}}_2} - \mathcal{N}_+) &\lesssim C N^{-\alpha/2} (\mathcal{N}_++1), \\ e^{-s{\mathcal{B}}_2}(\mathcal{N}_++1)^n e^{s{\mathcal{B}}_2} &\lesssim (\mathcal{N}_++1)^n. \end{align*} $$

In the next proposition, we describe the action of the operator $\mathcal {B}_2$ , defined as in equation (13), on the Hamilton operator in equation (9).

Proposition 4. We have

(23) $$ \begin{align} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} = 4 \pi \mathfrak{a}_N (N-1) + \sum_{|p|\leq N^{\alpha}} \frac{(4\pi \mathfrak{a}_N)^2}{p^2} + H_1 + \tilde H_2 + \tilde Q_2 + Q_3 + Q_4 + \mathcal{E}_{{\mathcal{B}}_2}, \end{align} $$

with

(24) $$ \begin{align} \tilde H_2 &= (2\hat V(0) - 8\pi \mathfrak{a}_N) \mathcal{N}_+, \end{align} $$
(25) $$ \begin{align} \tilde Q_2 &= \sum_{|p| \leq N^{\alpha}} 4\pi \mathfrak{a}_N \Big[ a_p^{\dagger} a_{-p}^{\dagger} \frac{a_0 a_0}{N} + \mathrm{h.c.} \Big] , \end{align} $$

and

$$\begin{align*}\pm\mathcal{E}_{{\mathcal{B}}_2} \lesssim N^{-\alpha/2} Q_4 + \big[N^{-\alpha/2} + N^{-1+5\alpha/2} \big] (\mathcal{N}_+ + 1) + N^{-1+\alpha} \mathcal{N}_+^2 + N^{-2} H_1. \end{align*}$$

To show Proposition 4, we define

(26) $$ \begin{align} \Gamma_2 := [H_1 + Q_4,{\mathcal{B}}_2] + Q_2 - \tilde Q^{\prime}_2 \end{align} $$

with

(27) $$ \begin{align} \tilde Q^{\prime}_2 = \sum_p \hat W (p) a^{\dagger}_p a^{\dagger}_{-p} a_0 a_0 + \mathrm{{h.c.}}, \end{align} $$

and

$$\begin{align*}\hat{W} (p) = \frac{1}{2} \chi_{|p| \leq N^{\alpha}} \Big[ \sum_q \hat{V}_N (p-q) {\varphi}_q + \hat{V}_N (p) \Big] - \frac{1}{2} \sum_{|q| \leq N^{\alpha}} \hat{V}_N (p-q) {\varphi}_q. \end{align*}$$

We observe that

(28) $$ \begin{align} \begin{aligned} e^{-{\mathcal{B}}_2} &H_N e^{{\mathcal{B}}_2} \\ &= H_0+ H_1 + Q_4 + \int_0^1 e^{-t{\mathcal{B}}_2} [H_1+Q_4,{\mathcal{B}}_2] e^{t{\mathcal{B}}_2} \mathrm{d}t + e^{-{\mathcal{B}}_2} Q_2 e^{{\mathcal{B}}_2} + e^{-{\mathcal{B}}_2}(H_2 + Q_3)e^{{\mathcal{B}}_2} \\ &= H_0+ H_1 + Q_4 + \int_0^1 e^{-t{\mathcal{B}}_2} (-Q_2 + \tilde Q^{\prime}_2 + \Gamma_2) e^{t{\mathcal{B}}_2} \mathrm{d}t + e^{-{\mathcal{B}}_2} Q_2 e^{{\mathcal{B}}_2}+e^{-{\mathcal{B}}_2}(H_2 + Q_3)e^{{\mathcal{B}}_2} \\ & = H_0+ H_1 + \tilde{Q}^{\prime}_2 + Q_4 + \int_0^1 \int_s^1 e^{-t{{\mathcal{B}}_2}}[Q_2,{{\mathcal{B}}_2}]e^{t{{\mathcal{B}}_2}} \mathrm{d}t \mathrm{d}s+ \int_0^1 \int_0^s e^{-t \mathcal{B}_2} \tilde{Q}^{\prime}_2 e^{t \mathcal{B}_2} \mathrm{d}t \mathrm{d}s \\ &\qquad + \int_0^1 e^{-t{\mathcal{B}}_2} \Gamma_2 e^{t{\mathcal{B}}_2} \mathrm{d}t + e^{-{\mathcal{B}}_2}(H_2 + Q_3)e^{{\mathcal{B}}_2} , \end{aligned} \end{align} $$

where, in the last step, we use

$$\begin{align*}\begin{aligned} \int_0^1 e^{-t \mathcal{B}_2} \tilde Q^{\prime}_2 e^{t \mathcal{B}_2} \mathrm{d}t &= \tilde Q^{\prime}_2 + \int_0^1 \int_0^s e^{-t \mathcal{B}_2} [\tilde{Q}^{\prime}_2 , \mathcal{B}_2 ] e^{t \mathcal{B}_2} \mathrm{d}t \mathrm{d}s \\ e^{-{\mathcal{B}}_2} Q_2 e^{{\mathcal{B}}_2} - \int_0^1 e^{-t{\mathcal{B}}_2} Q_2 e^{t{\mathcal{B}}_2} \mathrm{d}t &= \int_0^1 \int_s^1 e^{-t{{\mathcal{B}}_2}}[Q_2,{{\mathcal{B}}_2}]e^{t{{\mathcal{B}}_2}} \mathrm{d}t \mathrm{d}s. \end{aligned} \end{align*}$$

The proof of Proposition 4 now follows by controlling the terms on the right-hand side of equation (28). This is accomplished through a series of lemmas. We start by controlling the contribution arising from $H_2$ .

Lemma 5. On $\{\mathcal {N} = N\}$ , we have

(29) $$ \begin{align} e^{-{\mathcal{B}}_2} H_2 e^{{\mathcal{B}}_2} = \hat V (0) \mathcal{N}_+ + \mathcal{E}_{H_2} , \end{align} $$

with

$$ \begin{align*}\pm \mathcal{E}_{H_2} \lesssim N^{-\alpha/2} (\mathcal{N}_++1) + \mathcal{N}_+^2/N + N^{-2} H_1.\end{align*} $$

Proof. We have

(30) $$ \begin{align} e^{-{\mathcal{B}}_2}H_2e^{{\mathcal{B}}_2}=N \sum_{p} \hat V_N(p) &\left(a^{\dagger}_p a_p + \int_0^1 e^{-s{\mathcal{B}}_2} [ a^{\dagger}_p a_p,{\mathcal{B}}_2] e^{s{\mathcal{B}}_2} \mathrm{d}s\right) \notag \\ &\qquad - e^{-{\mathcal{B}}_2} \left( \sum_{p} \hat V_N(p) a^{\dagger}_p a_p\mathcal{N}_+ + \frac {\hat V_N(0)}{2} \mathcal{N}_+(\mathcal{N}_+ -1) \right) e^{{\mathcal{B}}_2}. \end{align} $$

The term on the second line is controlled using Lemma 3 by $N^{-1} \mathcal {N}_+^2$ . For the second term in parenthesis in the first line, we use Lemma 3 to estimate

$$ \begin{align*} \pm \sum_{p} N \hat V_N(p) [ a^{\dagger}_p a_p,{\mathcal{B}}_2] &= \pm \sum_{p} N \hat V_N(p) \tilde {\varphi}_p a^{\dagger}_p a^{\dagger}_{-p} a_0 a_0 +\text{h.c.} \\ &\leq C \|\tilde {\varphi} \|_{2} (\mathcal{N}_0+1) (\mathcal{N}_++1) \leq N^{- \alpha/2} (\mathcal{N}_+ + 1), \end{align*} $$

where we used that $\mathcal {N}_0\lesssim N$ and Lemma 2. Finally, since $V_N$ is even, we obtain

$$ \begin{align*} N |\hat V_N (p) - \hat V_N(0)| \leq \|x^2V\|_{1} N^{-2} p^2 \leq C N^{-2} p^2, \end{align*} $$

which gives

$$ \begin{align*} \pm \left(N \sum_{p} \hat V_N(p) a^{\dagger}_p a_p - \hat V(0) \mathcal{N}_+\right) \leq N^{-2} H_1.\\[-50pt] \end{align*} $$

The estimate of the term involving $Q_3$ is obtained analogously to [Reference Boccato, Brennecke, Cenatiempo and Schlein7, Reference Hainzl15]. We repeat the proof for the sake of completeness.

Lemma 6. We have

$$ \begin{align*} e^{-{\mathcal{B}}_2} Q_3 e^{{\mathcal{B}}_2} = Q_3 + \mathcal{E}_{Q_3} , \end{align*} $$

with

$$ \begin{align*} \pm \mathcal{E}_{Q_3} \lesssim N^{-\alpha/2} (Q_4 + \mathcal{N}_+ +1). \end{align*} $$

Proof. We can rewrite

(31) $$ \begin{align} e^{-{\mathcal{B}}_2} Q_3 e^{{\mathcal{B}}_2} = \sum_{q,r} \hat{V}_N(r)\left[e^{-{\mathcal{B}}_2}a_{q+r}^{\dagger} a_{-r}^{\dagger} e^{{\mathcal{B}}_2} e^{-{\mathcal{B}}_2} a_{q} a_0 e^{{\mathcal{B}}_2} + \text{h.c.} \right]. \end{align} $$

Via Duhamel’s formula, we have

(32) $$ \begin{align} e^{-{\mathcal{B}}_2}a_{q+r}^{\dagger} a_{-r}^{\dagger} e^{{\mathcal{B}}_2} = a^{\dagger}_{q+r} a^{\dagger}_{-r} + \int_0^1 e^{-s{\mathcal{B}}_2} [a^{\dagger}_{q+r} a^{\dagger}_{-r}, {\mathcal{B}}_2] e^{s{\mathcal{B}}_2} \mathrm{d}s \end{align} $$

and

$$ \begin{align*} e^{-{\mathcal{B}}_2} a_{q} a_0 e^{{\mathcal{B}}_2} = a_q a_0 + \int_0^1 e^{-s{\mathcal{B}}_2} (\tilde {\varphi}_q a^{\dagger}_{-q} a_0 a_0 a_0 - a_q \sum_l \tilde {\varphi}_l a_{-l} a_l a_0^{\dagger}) e^{s{\mathcal{B}}_2} \mathrm{d}s. \end{align*} $$

The product of the first terms, combined with its hermitian conjugate, corresponds to $Q_3$ in the statement of the lemma. All other terms will be estimated in three steps.

Step 1. Passing to x-space, we find

(33) $$ \begin{align} \begin{aligned} \Big| \sum_{q,r}& \tilde {\varphi}_q \hat{V}_N(r) \int_0^1 \big\langle \xi , a_{q+r}^{\dagger} a_{-r}^{\dagger} e^{-s{\mathcal{B}}_2} a^{\dagger}_{-q} a_0 a_0 a_0 e^{s{\mathcal{B}}_2} \xi \rangle \, \mathrm{ d}s \Big| \\ &= \int_0^1 \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) \big\langle \xi, {\check a}^{\dagger}_x {\check a}_y^{\dagger} e^{-s{\mathcal{B}}_2} {\check a}^{\dagger}(\check{\tilde{\varphi}}_x) a_0 a_0 a_0 e^{s{\mathcal{B}}_2} \xi \big\rangle \, \mathrm{d}s \\ &\leq \Big( \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) \| {\check a}_x {\check a}_y \xi \|^2 \Big)^{1/2} \\ &\hspace{2cm} \times \Big( \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N (x-y) \int_0^1 \| {\check a}^{\dagger}(\check{\tilde{\varphi}}_x) a_0 a_0 a_0 e^{s{\mathcal{B}}_2} \xi \|^2 \mathrm{d}s \Big)^{1/2} \\ &\leq C N^{3/2} \| V_N \|_1^{1/2} \| \tilde{\varphi} \|_2 \| Q_4^{1/2} \xi \| \| (\mathcal{N}_+ + 1)^{1/2} \xi \| \\ &\leq C N^{-\alpha/2} \langle \xi , (Q_4 + \mathcal{N}_+ + 1) \xi \rangle, \end{aligned} \end{align} $$

where we used Lemma 2, Lemma 3 and the bound $\mathcal {N}_0 \leq N$ .

Step 2. Similarly, we have

(34) $$ \begin{align} \begin{aligned} \Big| \sum_{q,r} &\hat{V}_N(r) \sum_l \tilde {\varphi}_l \int_0^1 \mathrm{d}s \, \big\langle \xi, a_{q+r}^{\dagger} a_{-r}^{\dagger} e^{-s{\mathcal{B}}_2} a_{q} a_{-l} a_l a^{\dagger}_0 e^{s{\mathcal{B}}_2} \xi \big\rangle \Big| \\ &= \Big| \sum_l \tilde {\varphi}_l \int_0^1 \mathrm{d}s \int_{\Lambda^2} \mathrm{ d}x \mathrm{d}y \, V_N(x-y) \langle \xi , a^{\dagger}_x a^{\dagger}_y e^{-s{\mathcal{B}}_2} a_x a_{-l} a_l a^{\dagger}_0 e^{s{\mathcal{B}}_2} \xi \rangle \Big| \\ &\leq C \| \tilde {\varphi} \|_2 \|V_N \|^{1/2}_1 \| Q_4^{1/2} \xi \| \| (\mathcal{N}_++1)^{2} \xi \| \leq C N^{-\alpha/2} \langle \xi , Q_4 + \mathcal{N}_+ + 1) \xi \rangle. \end{aligned} \end{align} $$

Step 3. The remaining term has the form

$$ \begin{align*} \sum_{q,r} \hat{V}_N(r)\int_0^1e^{-s{\mathcal{B}}_2}[a_{q+r}^{\dagger} a_{-r}^{\dagger},{\mathcal{B}}_2] e^{(s-1){\mathcal{B}}_2} a_{q} a_0 e^{{\mathcal{B}}_2} \mathrm{ d}s. \end{align*} $$

A straightforward computation gives

$$ \begin{align*} [a_{q+r}^{\dagger} a_{-r}^{\dagger},{\mathcal{B}}_2] = - \left( \tilde{\varphi}_r a^{\dagger}_{q+r} a_r + \tilde{\varphi}_{-r -q} a^{\dagger}_{-r} a_{-q-r} \right)a^{\dagger}_0a^{\dagger}_0. \end{align*} $$

The contributions of the two terms in parenthesis can be handled similarly. Let us consider, for example, the expectation

(35) $$ \begin{align} \begin{aligned} \Big| \sum_{q,r} \hat{V}_N(r) \tilde {\varphi}_r \int_0^1& \langle \xi, e^{-s{\mathcal{B}}_2}a_{q+r}^{\dagger} a_{r} a^{\dagger}_0 a^{\dagger}_0e^{(s-1){\mathcal{B}}_2} a_{q} a_0 e^{{\mathcal{B}}_2} \xi \rangle \, \mathrm{d}s \Big| \\ &\lesssim \frac{1}{N} \sum_{r,q} |\tilde{\varphi}_r| \int_0^1 \| a_{q+r} a_0 a_0 e^{(s-1){\mathcal{B}}_2} \xi\| \| a_r e^{s {\mathcal{B}}_2} a_q a_0 e^{{\mathcal{B}}_2} \xi \| \mathrm{d}s \\ &\lesssim \| \tilde {\varphi} \|_2 \| \mathcal{N}_+^{1/2} \xi \| \int_0^1 \Big( \sum_{r,q} \| a_r e^{(s-1) {\mathcal{B}}_2} a_q a_0 e^{{\mathcal{B}}_2} \xi \|^2 \Big)^{1/2} \mathrm{d}s \\ &\lesssim \| \tilde {\varphi}\|_2 \| \mathcal{N}_+^{1/2} \xi \| \| (\mathcal{N}_+ +1)^{3/2} \xi \| \leq N^{-\alpha/2} \langle \xi , (\mathcal{N}_+ + 1) \xi \rangle, \end{aligned} \end{align} $$

where we used Lemma 2 and Lemma 3.

Next, we recall the definition of $\Gamma_2$ given in equation (26) and consider the term containing $\Gamma _2$ appearing on the right-hand side of equation (28).

Lemma 7. We have

(36) $$ \begin{align} \Gamma_2 = \sum_{r,p,q} \hat{V}_N(r) \tilde {\varphi}_p[a_{p+r}^{\dagger} a_q^{\dagger} a_{-p}^{\dagger} a_{q+r}a_0a_0 + \text{h.c.}] \end{align} $$

and

(37) $$ \begin{align} \int_0^1 e^{-t{\mathcal{B}}_2} \Gamma_2 e^{t{\mathcal{B}}_2} \mathrm{d}t \lesssim N^{-\alpha/2} (Q_4 + \mathcal{N}_+ + 1). \end{align} $$

Proof. Straightforward calculations yield

$$ \begin{align*} [H_1,{\mathcal{B}}_2] = \sum_{p } p^2 \tilde {\varphi}_p [a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0 + \text{h.c.}] \end{align*} $$

and

(38) $$ \begin{align} [Q_4, {\mathcal{B}}_2] = \frac{1}{2} \sum_{p,q} \hat{V}_N(p-q) \tilde {\varphi}_q a_{p}^{\dagger} a_{-p}^{\dagger} a_0a_0 + \sum_{r,p,q} \hat{V}_N(r) \tilde {\varphi}_p a_{p+r}^{\dagger} a_q^{\dagger} a_{-p}^{\dagger} a_{q+r}a_0 a_0 + \text{h.c.} \end{align} $$

Hence,

$$\begin{align*}\begin{aligned} [H_1 + Q_4,{\mathcal{B}}_2] + Q_2 = \; &\sum_{|p|> N^{\alpha}} \Big(p^2 {\varphi}_p + \frac{1}{2} \sum_{q} \hat{V}_N(p-q) {\varphi}_q + \frac{1}{2} \hat{V}_N (p)\Big) (a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0+ \text{h.c.} ) \\ &- \frac{1}{2} \sum_{|p| > N^{\alpha} , |q| \leq N^{\alpha}} \hat{V}_N (p-q) {\varphi}_q (a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0+ \text{h.c.} ) \\ &+ \frac{1}{2} \sum_{|p| \leq N^{\alpha}} \Big( \sum_{|q| > N^{\alpha}} \hat{V}_N (p-q) {\varphi}_q + \hat{V}_N (p) \Big) (a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0+ \text{h.c.} ) \\ &+ \sum_{r,p,q} \hat{V}_N(r) \tilde {\varphi}_p a_{p+r}^{\dagger} a_q^{\dagger} a_{-p}^{\dagger} a_{q+r}a_0 a_0 \\ = \: &\tilde Q^{\prime}_2 + \sum_{r,p,q} \hat{V}_N(r) \tilde {\varphi}_p a_{p+r}^{\dagger} a_q^{\dagger} a_{-p}^{\dagger} a_{q+r}a_0 a_0, \end{aligned} \end{align*}$$

where we used the scattering equation (14) and the definition of $\tilde Q'_2$ given in equation (27). Comparing with equation (26), we find equation (36).

To prove the estimate (37), we write

$$\begin{align*}\int_0^1 e^{-t \mathcal{B}_2} \Gamma_2 e^{t \mathcal{B}_2} \mathrm{d}t = \sum_{p,q,r} \hat{V}_{N}(r)\tilde{\varphi}_p \int_0^1 e^{-t{\mathcal{B}}_2} a_{p+r}^{\dagger} a_q^{\dagger} e^{t{\mathcal{B}}_2}e^{-t{\mathcal{B}}_2} a_{-p}^{\dagger} a_{q+r} a_0 a_0 e^{t{\mathcal{B}}_2} \mathrm{d}t, \end{align*}$$

and we proceed similarly as in Lemma 6. We omit further details.

Next, we focus on the contribution with the commutator $[Q_2, \mathcal {B}_2]$ on the right-hand side of equation (28).

Lemma 8. On $\{ \mathcal {N} = N \}$ , we have

(39) $$ \begin{align} \begin{aligned} \int_0^1 & \int_s^1 e^{-t{\mathcal{B}}_2}[Q_2,{\mathcal{B}}_2]e^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{d}s \\ &= - \frac {(N-1)}{2} (\hat V(0) -8\pi \mathfrak{a}_N ) -\frac { N(N-1)}{2} \sum_{|p|\leq N^{\alpha}} \hat{V}_N(p) {\varphi}_p + \mathcal{N}_+ (\hat V(0) -8\pi \mathfrak{a}_N ) + \mathcal{E}_{[Q_2,{\mathcal{B}}_2]} , \end{aligned} \end{align} $$

and

(40) $$ \begin{align} \pm \mathcal{E}_{[Q_2,{\mathcal{B}}_2]} \lesssim (N^{-\alpha/2} + N^{\alpha -1}) \mathcal{N}_+ + N^{-1} (\mathcal{N}_+ + 1)^2 + N^{-\alpha/2} Q_4. \end{align} $$

Proof. First, we claim that

(41) $$ \begin{align} \int_0^1 \int_s^1 e^{-t{\mathcal{B}}_2}[Q_2,{\mathcal{B}}_2]e^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{d}s &= \frac { N(N-1)}{2} \sum_p \hat{V}_N(p) \tilde {\varphi}_p -N \mathcal{N}_+ \sum_p \hat{V}_N(p) \tilde {\varphi}_p + \mathcal{E}^{\prime}_{[Q_2,{\mathcal{B}}_2]} \end{align} $$

with

(42) $$ \begin{align} \mathcal{E}^{\prime}_{[Q_2,{\mathcal{B}}_2]} = - 2N &\sum_p \hat{V}_N (p) \tilde{\varphi}_p \int_0^1 \int_s^1 \left[ e^{-t \mathcal{B}_2} \mathcal{N}_+ e^{t \mathcal{B}_2} - \mathcal{N}_+ \right] \mathrm{d}t \mathrm{ d}s \nonumber \\ &\quad + \sum_p \hat{V}_N(p) \tilde{\varphi}_p \int_0^1 \int_s^1 e^{-t{\mathcal{B}}_2} \mathcal{N}_+( \mathcal{N}_+ + 1 )e^{t{\mathcal{B}}_2}\mathrm{d}t \mathrm{d}s \nonumber \\ &\quad +2\sum_p \hat{V}_N(p) \tilde {\varphi}_p \int_0^1 \int_s^1 e^{-t{\mathcal{B}}_2} \mathcal{N}_0(\mathcal{N}_0 -1) a^{\dagger}_p a_pe^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{d}s \nonumber \\ &\quad -\sum_{p,q}\hat{V}_N(p) \tilde {\varphi}_q \int_0^1 \int_s^1 e^{-t{\mathcal{B}}_2} a^{\dagger}_{p} a^{\dagger}_{-p} a_{-q} a_{q} (2 \mathcal{N}_0+1 ) e^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{d}s. \end{align} $$

To prove equation (42), we calculate

(43) $$ \begin{align}\nonumber [Q_2,{\mathcal{B}}_2] &= \frac{1}{4} \sum_{p,q} \hat{V}_N(p) \tilde{\varphi}_q [a_p^{\dagger} a_{-p}^{\dagger} a_0a_0 + a_{-p}a_pa_0^{\dagger} a_0^{\dagger} , a_q^{\dagger} a_{-q}^{\dagger} a_0a_0 - a_{-q}a_q a_0^{\dagger} a_0^{\dagger}] \\ &= \frac{1}{4} \sum_{p,q} \hat{V}_N(p) \tilde{\varphi}_q \left( [a_{-p}a_pa_0^{\dagger} a_0^{\dagger} , a_q^{\dagger} a_{-q}^{\dagger} a_0a_0] - [a_p^{\dagger} a_{-p}^{\dagger} a_0a_0,a_{-q}a_q a_0^{\dagger} a_0^{\dagger}] \right). \end{align} $$

The two terms in brackets are hermitian conjugates. Hence, it suffices to compute the second one

$$ \begin{align*} - [a_p^{\dagger} a_{-p}^{\dagger} a_0a_0,a_{-q}a_q a_0^{\dagger} a_0^{\dagger}] &= [a_{-q}a_q, a_p^{\dagger} a_{-p}^{\dagger}] a_0^{\dagger} a_0^{\dagger} a_0a_0 - a_p^{\dagger} a_{-p}^{\dagger} a_{-q} a_q [a_0a_0,a_0^{\dagger} a_0^{\dagger}], \end{align*} $$

where

(44) $$ \begin{align} [a_{-p}a_p,a_q^{\dagger} a_{-q}^{\dagger} ] = (\delta_{p,q} + \delta_{p,-q})( 1 + a^{\dagger}_p a_p + a^{\dagger}_{-p} a_{-p} ), \end{align} $$

and

(45) $$ \begin{align} \begin{aligned} a_0^{\dagger} a_0^{\dagger} a_0 a_0 &= \mathcal{N}_0(\mathcal{N}_0-1)= N(N -1) - 2 N \mathcal{N}_+ + \mathcal{N}_+(\mathcal{N}_+ +1), \\ [a_0a_0,a^{\dagger}_0a^{\dagger}_0] &= 2(2 \mathcal{N}_0 + 1). \end{aligned} \end{align} $$

Inserting these identities on the right-hand side of equation (43), conjugating with $e^{t\mathcal {B}_2}$ and integrating over $t\ \mathrm{and}\ s$ , we obtain equation (42).

With the definition given in equation (16), we write

$$\begin{align*}\begin{aligned} \frac{N(N-1)}{2} \sum_p \hat{V}_N (p) \tilde{{\varphi}}_p &= \frac{(N-1)}{2} (\hat{V} (0) - 8 \pi \mathfrak{a}_N) - \frac{N(N-1)}{2} \sum_{|p| \leq N^{\alpha}} \hat{V}_N (p) {\varphi}_p , \end{aligned} \end{align*}$$

and we use the bound (19) to estimate

$$\begin{align*}\pm \Big[ - N \mathcal{N}_+ \sum_p \hat{V}_N (p) \tilde{\varphi}_p + (8\pi \mathfrak{a}_N - \hat{V} (0)) \mathcal{N}_+ \Big] \lesssim N^{\alpha-1} \mathcal{N}_+. \end{align*}$$

Thus, Lemma 8 follows from equation (41), if we can prove that $\mathcal {E}^{\prime }_{[Q_2,{\mathcal {B}}_2]}$ satisfies the estimate given in (40).

Using the bound

$$\begin{align*}\Big| \sum_p \hat{V}_N (p) \tilde{\varphi}_p \Big| \lesssim N^{-1} \end{align*}$$

and Lemma 3, we can bound the first term on the right-hand side of equation (42) by

$$\begin{align*}\pm 2N \sum_p \hat{V}_N (p) \tilde{\varphi}_p \int_0^1 \int_s^1 \left[ e^{-t \mathcal{B}_2} \mathcal{N}_+ e^{t \mathcal{B}_2} - \mathcal{N}_+ \right] \mathrm{d}t \mathrm{d}s \lesssim N^{-\alpha/2} (\mathcal{N}_+ + 1). \end{align*}$$

Also, the second term on the right-hand side of equation (42) can be bounded with Lemma 3; we find

$$\begin{align*}\pm \sum_p \hat{V}_N (p) \tilde{\varphi}_p \int_0^1 \int_s^1 e^{-t \mathcal{B}_2} \mathcal{N}_+ ( \mathcal{N}_+ + 1) e^{t \mathcal{B}_2} \mathrm{d}t \mathrm{d}s \lesssim N^{-1} (\mathcal{N}_+ + 1)^2. \end{align*}$$

For the third term, we use $\| \hat {V}_N \|_{\infty } \lesssim N^{-1}$ , $\| \tilde {\varphi } \|_{\infty } \lesssim N^{-1-2\alpha }$ together with $\mathcal {N}_0 (\mathcal {N}_0 -1) \leq N^2$ and Lemma 3 again to conclude that

$$\begin{align*}\pm 2 \sum_p \hat{V}_N (p) \tilde{\varphi}_p \int_0^1 \int_s^1 e^{-t \mathcal{B}_2} \mathcal{N}_0 ( \mathcal{N}_0 + 1) e^{t \mathcal{B}_2} \mathrm{d}t \mathrm{d}s \lesssim N^{-2\alpha} \mathcal{N}_+. \end{align*}$$

To control the last term on the right-hand side of equation (42), we write

(46) $$ \begin{align} \sum_{p,q} \hat{V}_N(p) \tilde {\varphi}_q \int_0^1 \int_s^1 e^{-t{\mathcal{B}}_2} & a^{\dagger}_{p} a^{\dagger}_{-p} a_{-q} a_{q} (2 \mathcal{N}_0+1) e^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{d}s \notag \\ &\quad =\sum_{p} \hat{V}_N(p) \int_0^1 \int_s^1 e^{-t{\mathcal{B}}_2} a^{\dagger}_{p} a^{\dagger}_{-p} e^{t{\mathcal{B}}_2} e^{-t{\mathcal{B}}_2} \Phi (2 \mathcal{N}_0+1) e^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{ d}s, \end{align} $$

where we define $\Phi = \sum _q \tilde {\varphi }_q a_{-q} a_q$ so that, by Lemma 2,

(47) $$ \begin{align} \| \Phi \xi \| \lesssim \| \tilde{\varphi} \|_2 \| \mathcal{N}_+ \xi \| \lesssim N^{-1-\alpha/2} \| \mathcal{N}_+ \xi \|. \end{align} $$

Next, we expand

(48) $$ \begin{align} e^{- t{\mathcal{B}}_2} a^{\dagger}_{p} a^{\dagger}_{-p} e^{t{\mathcal{B}}_2} = a^{\dagger}_{p} a^{\dagger}_{-p} + \int_0^t e^{-\tau {\mathcal{B}}_2} [a^{\dagger}_{p} a^{\dagger}_{-p} ,{\mathcal{B}}_2]e^{\tau {\mathcal{B}}_2}d\tau. \end{align} $$

Inserting this identity into equation (46), we obtain two contributions. The first contribution can be controlled by passing to position space. We find

$$\begin{align*}\begin{aligned} \pm \sum_{p} \hat{V}_N(p) &\int_0^1 \int_s^1 a^{\dagger}_{p} a^{\dagger}_{-p} e^{-t{\mathcal{B}}_2} \Phi (2 \mathcal{N}_0+1) e^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{d}s \\ &= \pm \int_0^1 \int_s^1 \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y \kappa V_N(x-y) {\check a}_x^{\dagger} {\check a}_y^{\dagger} e^{-t{\mathcal{B}}} \Phi (2 \mathcal{N}_0 +1) e^{t{\mathcal{B}}} \mathrm{d}t \mathrm{d}s + \text{h.c.} \\ &\lesssim \delta Q_4 + \delta^{-1} N^2 \|\varphi\|_{L^2}^2 \|V_N\|_{1} (\mathcal{N}_++1)^2 \lesssim N^{-\alpha/2} Q_4 + N^{-1-\alpha/2} (\mathcal{N}_++1)^2 .\end{aligned} \end{align*}$$

On the other hand, the contribution arising from the second term on the right-hand side of equation (48) can be controlled by

(49) $$ \begin{align} \pm \sum_{p} \hat{V}_N(p) \tilde {\varphi}_p \int_0^1 \int_s^1\int_0^t e^{-\tau {\mathcal{B}}_2} a^{\dagger}_0 a^{\dagger}_0 (2 a^{\dagger}_p a_p + 1 )e^{\tau {\mathcal{B}}_2} e^{-t{\mathcal{B}}} \Phi &(2 \mathcal{N}_0 +1) e^{t{\mathcal{B}}} \mathrm{d}t \mathrm{d}s + \text{h.c.} \notag \\ &\qquad \lesssim N^{-1 -\alpha/2} (\mathcal{N}_++1)^2. \end{align} $$

This concludes the proof of (39) and (40).

Finally, we control the contribution with the commutator $[ \tilde {Q}_2 , \mathcal {B}_2 ]$ in equation (28).

Lemma 9. We have

(50) $$ \begin{align} \int_0^1 \int_0^s e^{-t{\mathcal{B}}_2}[\tilde Q_2,{\mathcal{B}}_2] e^{t{\mathcal{B}}_2} \mathrm{d}t \mathrm{d}s= - \frac {N(N-1)} 2 \sum_{|p|> N^{\alpha}, |q|\leq N^{\alpha}} \hat{V}_N(p-q) {\varphi}_p {\varphi}_q + \mathcal{E}_{[\tilde Q_2,{\mathcal{B}}_2]}, \end{align} $$

with

$$ \begin{align*} \pm \mathcal{E}_{[\tilde Q_2 , {\mathcal{B}}_2]} \lesssim N^{-\alpha/2} Q_4 + N^{-1 + \alpha} (\mathcal{N}_++1)^2. \end{align*} $$

Proof. Using that $ \chi _{|p| \leq N^{\alpha }} \tilde \varphi _p = 0$ , similar computations as in the proof of Lemma 41 yield

(51) $$ \begin{align} \begin{aligned} [\tilde Q_2 &,{\mathcal{B}}_2] \\ = \; &\frac{1}{4} \sum_{|p|>N^{\alpha}, |q|\leq N^{\alpha},r} \hat V_N(p-q) {\varphi}_q \tilde {\varphi}_r \\ &\times (a^{\dagger}_p a^{\dagger}_{-p} a_r a_{-r}[a_0 a_0,a^{\dagger}_0 a^{\dagger}_0] +[a^{\dagger}_p a^{\dagger}_{-p},a_r a_{-r}]a^{\dagger}_0 a^{\dagger}_0 a_0 a_0)+\text{h.c.} \\ = \; &\sum_{|p|>N^{\alpha}, |q| \leq N^{\alpha}} \frac{1}{4}\hat V_N(p-q) {\varphi}_q \\ & \times \big( (2 a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1) + \text{h.c.}) - 4 \tilde {\varphi}_p (a^{\dagger}_p a_p +a^{\dagger}_{-p} a_{-p}) \mathcal{N}_0 (\mathcal{N}_0-1) -4\tilde{\varphi}_p \mathcal{N}_0 (\mathcal{N}_0-1)\big), \end{aligned} \end{align} $$

with the notation $\Phi = \sum _r \tilde {\varphi }_r a_r a_{-r}$ . To bound the contribution arising from the first term in parentheses, we decompose

(52) $$ \begin{align} \begin{aligned} \frac{1}{2} &\sum_{|p|>N^{\alpha}, |q|\leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1)\\ = \; &\frac{1}{2} \sum_{p,|q|\leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1) - \frac{1}{2} \sum_{|p|, |q|\leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1). \end{aligned} \end{align} $$

The first term can be controlled by switching to position space. With the notation $\check {{\varphi }}^<$ for the Fourier series of $\chi _{|q| \leq N^{\alpha }} {\varphi }_q$ , we find

$$\begin{align*}\begin{aligned} \pm \frac{1}{2} &\sum_{p,|q|\leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1) \\ &= \pm \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N (x-y) \check {\varphi}^< (x-y) a^{\dagger}_x a^{\dagger}_y \Phi (2\mathcal{N}_0+1) + \text{h.c.} \nonumber \\ &\lesssim \delta Q_4 + \delta^{-1} N^{-\alpha} \|V_N^{1/2} \check{\varphi}^< \|_2^2 (\mathcal{N}_++1)^2, \end{aligned} \end{align*}$$

where we used the bound (47) for $\Phi $ and $\mathcal {N}_0 \leq N$ . With

(53) $$ \begin{align} \| V_N^{1/2} \check{\varphi}^< \|_2^2 = \sum_{|p|, |q| \leq N^{\alpha}} \hat{V}_N (p-q) {\varphi}_p {\varphi}_q \lesssim \frac{1}{N^3} \Big[ \sum_{|p| \leq N^{\alpha}} \frac{1}{|p|^2} \Big]^2 \lesssim N^{2\alpha - 3} \end{align} $$

and choosing $\delta = N^{-\alpha /2}$ , we conclude (since $\alpha < 1$ ) that

$$\begin{align*}\pm \frac{1}{2} \sum_{p,|q|\leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1) \lesssim N^{-\alpha/2} Q_4 + N^{-1-\alpha/2} (\mathcal{N}_+ + 1)^2. \end{align*}$$

For the second term on the right-hand side of equation (52), we estimate

$$\begin{align*}\pm \frac{1}{2} \sum_{|p|, |q|\leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1) + \text{h.c.} \lesssim N^{-\alpha/2} \| \chi_{|p|\leq N^{\alpha}} \hat V_N \ast {\varphi}^{<} \|_{2} \| (\mathcal{N}_+ + 1)^2, \end{align*}$$

where, again, we used the bound (47) and $\mathcal {N}_0 \leq N$ . With

$$\begin{align*}\| \chi_{|p|\leq N^{\alpha}} \hat V_N \ast {\varphi}^{<}\|_{2} \leq \|\chi_{|p|\leq N^{\alpha}}\|_{L^2} \|V_N \check{\varphi}^{<}\|_{1} \lesssim N^{3\alpha/2}\| V_N^{1/2}\|_{2} \| V_N^{1/2} \check {\varphi}^{<}\|_{2} \lesssim N^{-2+5\alpha/2} \end{align*}$$

we conclude that

$$\begin{align*}\pm \frac{1}{2} \sum_{|p|, |q|\leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q a^{\dagger}_p a^{\dagger}_{-p} \Phi (2\mathcal{N}_0+1) + \text{h.c.} \lesssim N^{-2-2\alpha} (\mathcal{N}_+ + 1)^2. \end{align*}$$

The contribution arising from the second term in parentheses on the right-hand side of equation (51) can be bounded by

$$ \begin{align*} \pm \sum_{|p|>N^{\alpha}, |q| \leq N^{\alpha}} \hat V_N(p-q) {\varphi}_q \tilde {\varphi}_p a^{\dagger}_p a_p \mathcal{N}_0 (\mathcal{N}_0-1) \lesssim N^{-1- \alpha} \mathcal{N}_+, \end{align*} $$

using that $ \|\tilde {\varphi } (\hat V_N \ast {\varphi }^{<})\|_{L^{\infty }} \lesssim \|\tilde {\varphi } \|_{\infty } \| V_N \check {\varphi }^{<}\|_1 \lesssim N^{-3-\alpha }$ . For the contribution arising from the last term on the right-hand side of equation (51), we write $\mathcal {N}_0(\mathcal {N}_0-1)= N(N -1) - 2 N \mathcal {N}_+ + \mathcal {N}_+ (\mathcal {N}_+ +1)$ . The contribution proportional to $N(N-1)$ produces the main term on the right-hand side of equation (50). The other contributions can be bounded, noticing that

$$ \begin{align*} \Big| \sum_{|p|> N^{\alpha}, |q|\leq N^{\alpha}} \hat{V}_N(p-q) {\varphi}_p {\varphi}_q \Big| \leq \|\check {\tilde {\varphi}} V_N \check {\varphi}^{<}\|_{1} \leq \|\check {\tilde {\varphi}}\|_{\infty} \|V_N \check{\varphi}^< \|_1 \lesssim N^{-2+\alpha}, \end{align*} $$

where we used $\| \check {\tilde {\varphi }} \|_{\infty } \lesssim \| \tilde {\varphi } \|_1 \lesssim 1$ , by Lemma 2.

We can now finish the proof of Proposition 4.

Proof of Proposition 4.

Combining equation (28) with the bounds proven in Lemma 5, Lemma 6, Lemma 7, Lemma 8 and Lemma 9, we conclude that

(54) $$ \begin{align} e^{-\mathcal{B}_2} H_N e^{\mathcal{B}_2} = 4 \pi \mathfrak{a}_N (N-1) + A_{\alpha} + H_1 + \tilde{H}_2 + \tilde{Q}^{\prime}_2 + Q_3 + Q_4 + \mathcal{E}, \end{align} $$

where

(55) $$ \begin{align} \pm \mathcal{E} \lesssim N^{-\alpha/2} Q_4 + \big [ N^{-\alpha/2} + N^{\alpha-1} \big] (\mathcal{N}_+ + 1) + N^{-1+\alpha} \mathcal{N}_+^2 + N^{-2} H_1 \end{align} $$

and where we defined

$$\begin{align*}\begin{aligned} A_{\alpha} &= - \frac{N(N-1)}{2} \Big[ \sum_{|p| \leq N^{\alpha}} \hat{V}_N (p) {\varphi}_p + \sum_{|p|> N^{\alpha}, |q| \leq N^{\alpha}} \hat{V}_N (p-q) {\varphi}_p {\varphi}_q \Big] \\ &= - \frac{N(N-1)}{2} \Big[ \sum_{|p| \leq N^{\alpha}} (\hat{V}_N (p) + \hat{V}_N * {\varphi}) {\varphi}_p - \sum_{|p|, |q| \leq N^{\alpha}} \hat{V}_N (p-q) {\varphi}_p {\varphi}_q \Big]. \end{aligned} \end{align*}$$

The second term in parentheses can be estimated as in (53). Setting, in position space, $f = 1 + \check {{\varphi }}$ , we find

(56) $$ \begin{align} A_{\alpha} = - \frac{N(N-1)}{2} \sum_{|p| \leq N^{\alpha}} (\hat{V}_N * \hat{f}) (p) {\varphi}_p + \mathcal{O} (N^{2\alpha-1}). \end{align} $$

From equation (16), we have $\widehat {V_N f} (0) = 8\pi \mathfrak {a}_N$ . Hence

(57) $$ \begin{align} | (\hat{V}_N * \hat{f}) (p) - 8 \pi \mathfrak{a}_N / N | \leq \int_{\Lambda} V_N (x) f (x) |e^{-ip \cdot x} -1| \mathrm{d}x \leq C |p| / N^2. \end{align} $$

Moreover, from the scattering equation (14), we find

$$\begin{align*}{\varphi}_p = - \frac{1}{2p^2} (\hat{V}_N * \hat{f}) (p), \end{align*}$$

which implies, by the bound (57),

$$\begin{align*}\Big| {\varphi}_p + \frac{4\pi \mathfrak{a}_N}{N p^2} \Big| \leq \frac{1}{2p^2} \big| (\hat{V}_N * \hat{f}) (p) - 8 \pi \mathfrak{a}_N \big| \leq \frac{C}{|p|N^2}. \end{align*}$$

Inserting in equation (56), we obtain

$$\begin{align*}A_{\alpha} = \sum_{|p| \leq N^{\alpha}} \frac{(4\pi \mathfrak{a}_N)^2}{p^2} + \mathcal{O} (N^{2\alpha-1}). \end{align*}$$

To conclude the proof of Proposition 4, we still have to compare the operator $\tilde Q^{\prime }_2$ appearing on the right-hand side of equation (54) with the operator $\tilde Q_2$ defined in equation (25). From equation (27), we can write, using again the notation $f = 1+ \check {{\varphi }}$ ,

(58) $$ \begin{align} \tilde Q^{\prime}_2 - \tilde Q_2 = \frac{1}{2} \sum_{|p| \leq N^{\alpha}} \Big[ (\hat{V}_N * \hat{f})(p) - \frac{8\pi \mathfrak{a}_N}{N} \Big] \, a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0 - \frac{1}{2} \sum_{p, |q| \leq N^{\alpha}} \hat{V}_N (p-q) {\varphi}_q a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0 + \text{h.c.} \end{align} $$

The first term on the right-hand side of equation (58) can be bounded using the bound (57) by

$$\begin{align*}\begin{aligned} \pm \sum_{|p| \leq N^{\alpha}} \Big[ (\hat{V}_N * \hat{f})(p) &- \frac{8\pi \mathfrak{a}_N}{N} \Big] \, \big( a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0 +\text{h.c.}\big) \\ &\lesssim \frac{1}{N} \| |p| \chi_{|p| \leq N^{\alpha}} \|_2 (\mathcal{N}_+ + 1) \lesssim N^{-1+5\alpha/2} (\mathcal{N}_+ + 1). \end{aligned} \end{align*}$$

For the second term on the right-hand side of equation (58), we set ${\varphi }^<_p = {\varphi }_p \chi _{|p| \leq N^{\alpha }}$ and estimate, switching to position space,

$$\begin{align*}\begin{aligned} \pm \sum_{p} &(\hat{V}_N * {\varphi}^< )(p) \, a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0 + \text{h.c.} \\ &= \pm \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y \, V_N (x-y) \check{{\varphi}}^< (x-y) a_x^{\dagger} a_y^{\dagger} a_0 a_0 + \text{h.c.} \\ &\lesssim \delta Q_4 + \delta^{-1} N^2 \| V_N \|_1 \| \check{{\varphi}}^< \|_{\infty}^2 \lesssim \delta Q_4 + \delta^{-1} N^{-1+2\alpha} \leq N^{-\alpha/2} Q_4 + N^{-1+5\alpha/2}, \end{aligned} \end{align*}$$

since $\| \check {{\varphi }}^< \|_{\infty } \lesssim \| {\varphi }^< \|_1 \lesssim N^{\alpha }$ , from Lemma 2 (in the last step, we chose $\delta = N^{-\alpha /2}$ ). The last two estimates show that the difference $\tilde Q^{\prime }_2 - \tilde Q_2$ can be added to the error (55) and therefore conclude the proof of Proposition 4.

4 Cubic renormalisation

While conjugation with $e^{{\mathcal {B}}_2}$ allowed us to renormalise the quadratic part of the Hamiltonian $H_N$ , regularising the off-diagonal term $Q_2$ , it did not significantly change the cubic operator $Q_3$ . To renormalise $Q_3$ , we proceed with a second conjugation, with a unitary operator $e^{{\mathcal {B}}_3}$ , where

(59) $$ \begin{align} {\mathcal{B}}_3 = \sum_{p,q} \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} \, a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 - \text{h.c.},\notag\\[-24pt] \end{align} $$

with the same $0 < \alpha <1$ used in the definition given in equation (17) of $\tilde {\varphi }$ . Similarly to what we did in Lemma 3 for the action of ${\mathcal {B}}_2$ , it is important to notice that conjugation with $e^{{\mathcal {B}}_3}$ does not substantially change the number of excitations.

Lemma 10. For all $s \in [-1;1]$ and all $k \in \mathbb {N}$ , we have

(60) $$ \begin{align} \pm \Big[ e^{-s{\mathcal{B}}_3} \mathcal{N}_+ e^{s {\mathcal{B}}_3} - \mathcal{N}_+ \Big] &\lesssim N^{-\alpha/2} (\mathcal{N}_+ + 1) \, , \end{align} $$
(61) $$ \begin{align} e^{-s{\mathcal{B}}_3} (\mathcal{N}_++1)^k e^{s{\mathcal{B}}_3} &\lesssim (\mathcal{N}_++1)^k. \end{align} $$

Proof. We proceed similarly as in the proof of [Reference Boccato, Brennecke, Cenatiempo and Schlein7, Prop. 5.1]. For $\xi \in \mathcal {F}$ , we set $f(s) = \langle \xi , e^{-s{\mathcal {B}}_3} (\mathcal {N}_++1) e^{s{\mathcal {B}}_3} \xi \rangle $ . For $s \in (0;1)$ , we find

(62) $$ \begin{align} \begin{aligned} f'(s) &= \langle \xi, e^{-s{\mathcal{B}}_3} [(\mathcal{N}_++1),{\mathcal{B}}_3] e^{s{\mathcal{B}}_3} \xi \rangle \\ &= \sum_{p,q} \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} \, \langle \xi, e^{-s{\mathcal{B}}_3} a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 e^{s\mathcal{B}_3} \xi \rangle + \text{h.c.} \\ &\lesssim \delta \sum_{p,q} \langle \xi, e^{-s{\mathcal{B}}_3} a^{\dagger}_{p+q} a^{\dagger}_{-p} a_{-p} a_{p+q} e^{s{\mathcal{B}}_3} \xi \rangle + \delta^{-1} \sum_{p,q} |\tilde {\varphi}_p|^2 \, \langle \xi, e^{-s{\mathcal{B}}_3} a^{\dagger}_q a_q a^{\dagger}_0 a_0 e^{s{\mathcal{B}}_3} \xi \rangle \\ &\lesssim \delta \langle e^{-s{\mathcal{B}}_3} \mathcal{N}_+^2 e^{s{\mathcal{B}}_3} \xi \rangle + C \delta^{-1} N \| \tilde {\varphi} \|_2^2 \langle \xi , e^{-s{\mathcal{B}}_3} \mathcal{N}_+ e^{s {\mathcal{B}}_3} \xi \rangle \lesssim N^{-\alpha/2} f(s), \end{aligned} \end{align} $$

where we put $\delta = N^{-1-\alpha /2}$ and used that $\mathcal {N}_+, \mathcal {N}_0 \leq N$ , $\| {\varphi } \|_2 \lesssim N^{-1-\alpha /2}$ , by Lemma 2. With Gronwall’s lemma [Reference Pachpatte26, Theorem 1.2.2], we obtain $f(s) \lesssim \langle \xi , (\mathcal {N}_+ + 1) \xi \rangle $ for all $s \in [-1;1]$ , proving the bound (61). Inserting this estimate on the right-hand side of (62) and integrating over s, we obtain the desired bound (60). For $k> 1$ , the bound given in (61) can be shown similarly.

The operator ${\mathcal {B}}_3$ is chosen (similarly as we did with ${\mathcal {B}}_2$ in Section 3) so that the commutator $[H_1 + Q_4 , {\mathcal {B}}_3]$ arising from conjugation with $e^{{\mathcal {B}}_3}$ cancels the main part of $Q_3$ . The goal of this section is to use this cancellation to prove the following proposition.

Proposition 11. We have

(63) $$ \begin{align} \begin{aligned} e^{-{\mathcal{B}}_3} &e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} \\ &= 4 \pi \mathfrak{a}_N (N-1) + \frac{1}{4}\sum_{|p|\leq N^{\alpha}} \frac{(8\pi \mathfrak{a}_N)^2}{p^2} \\ &\quad + \sum_{p} (p^2 + 8\pi \mathfrak{a}_N \chi_{|p| \leq N^{\alpha}}) a^{\dagger}_p a_p + \frac{1}{2}\sum_{|p|\leq N^{\alpha}} 8 \pi \mathfrak{a}_N [a^{\dagger}_p a^{\dagger}_{-p} \frac{a_0 a_0 }N+\text{h.c.}] + Q_4 + \mathcal{E}_{{\mathcal{B}}_3}, \end{aligned} \end{align} $$

with

(64) $$ \begin{align} \begin{aligned} \pm \mathcal{E}_{{\mathcal{B}}_3} \lesssim \; & N^{-3\alpha/2} H_1 + N^{-\alpha/2} Q_4 + N^{-\alpha/2} (\mathcal{N}_+ + 1) \\ &+ N^{(3\alpha-1)/2} (\mathcal{N}_+ + 1)^{3/2} + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{aligned} \end{align} $$

To prove Proposition 11, we define

(65) $$ \begin{align} \Gamma_3 := \big[ H_1 + Q_4 , \mathcal{B}_3 \big] + Q_3. \end{align} $$

Starting from equation (23), we compute

$$ \begin{align*}\nonumber e^{-{\mathcal{B}}_3} e^{-{\mathcal{B}}_2} &H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} - 4 \pi \mathfrak{a}_N (N-1) - \sum_{|p|\leq N^{\alpha}} \frac{(4\pi \mathfrak{a}_N)^2}{p^2} - e^{-{\mathcal{B}}_3}(\tilde H_2 + \tilde Q_2 + \mathcal E_{{\mathcal{B}}_2})e^{{\mathcal{B}}_3} \\ &= H_1 + Q_4 + \int_0^1 e^{-t{\mathcal{B}}_3} [H_1+Q_4,{\mathcal{B}}_3] e^{t{\mathcal{B}}_2} \mathrm{ d}t + e^{-{\mathcal{B}}_3} Q_3 e^{{\mathcal{B}}_3} \nonumber \\ \nonumber & = H_1 + Q_4 + \int_0^1 e^{-t{\mathcal{B}}_3} (-Q_3 + \Gamma_3) e^{t{\mathcal{B}}_3} \mathrm{d}t + e^{-{\mathcal{B}}_3} Q_3 e^{{\mathcal{B}}_3}, \end{align*} $$

which leads to

(66) $$ \begin{align} \begin{aligned} e^{-{\mathcal{B}}_3} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} = \; & 4 \pi \mathfrak{a}_N (N-1) + \sum_{|p|\leq N^{\alpha}} \frac{(4\pi \mathfrak{a}_N)^2}{p^2} + H_1 + Q_4 + e^{-{\mathcal{B}}_3}(\tilde H_2 + \tilde Q_2 + \mathcal E_{{\mathcal{B}}_2} )e^{{\mathcal{B}}_3} \\ &+ \int_0^1 e^{-t{\mathcal{B}}_3} \Gamma_3 e^{t{\mathcal{B}}_3} \mathrm{d}t + \int_0^1 \int_s^1 e^{-t{{\mathcal{B}}_3}}[Q_3,{{\mathcal{B}}_3}]e^{t{{\mathcal{B}}_3}} \mathrm{d}t \mathrm{d}s. \end{aligned} \end{align} $$

To show Proposition 11, we are going to control all terms on the right-hand side of equation (66). We start by computing and estimating the commutator in equation (65), defining the error term $\Gamma _3$ .

Lemma 12. We have

(67) $$ \begin{align} [H_1,{\mathcal{B}}_3] &= \sum_{p,q,r} \hat V_N(p-r) (\delta_{0,r}+{\varphi}_r) \chi_{|p|> N^{\alpha}} \chi_{|q| \leq N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.} + \mathcal{E}_{[H_1,{\mathcal{B}}_3]}, \end{align} $$
(68) $$ \begin{align} [Q_4,{\mathcal{B}}_3] &= \sum_{r,p,q} \hat V_N(p-r)\tilde {\varphi}_r \chi_{|q| \leq N^{\alpha}} (a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.})+ \mathcal{E}_{[Q_4,{\mathcal{B}}_3]}, \end{align} $$

where

(69) $$ \begin{align} \begin{aligned} \pm\mathcal{E}_{[H_1,{\mathcal{B}}_3]} &\lesssim N^{-3\alpha/2} H_1 + N^{-1+5\alpha/2} (\mathcal{N}_++1)^2, \\ \pm \mathcal{E}_{[Q_4,{\mathcal{B}}_3]} &\lesssim N^{-\alpha/2} Q_4 + N^{-1+5\alpha/2} (\mathcal{N}_++1)^2. \end{aligned} \end{align} $$

Proof. A simple computation shows that

$$ \begin{align*} [H_1,{\mathcal{B}}_3] &= \sum_{p,q} \left[(p+q)^2 + p^2 - q^2 \right] \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.}, \\ &= 2 \sum_{p,q} p^2 \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 +\text{h.c.}+ \mathcal{E}_{[H_1,{\mathcal{B}}_3]}, \end{align*} $$

with

$$ \begin{align*} \mathcal{E}_{[H_1,{\mathcal{B}}_3]}& = 2 \sum_{p,q} p\cdot q \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 +\text{h.c.} \end{align*} $$

Using the scattering equation (14) yields (67). We now estimate $\mathcal {E}_{[H_1,{\mathcal {B}}_3]}$ . Using $|q| \leq N^{\alpha }$ , we find, for any $\delta>0$ ,

$$ \begin{align*} \pm \mathcal{E}_{[H_1,{\mathcal{B}}_3]} &= \pm 2 \sum_{p,q} p\cdot q \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p} (\mathcal{N}_++1)^{-1/2} (\mathcal{N}_++1)^{1/2} a_q a_0 + \text{h.c.} \\ &\lesssim \delta \sum_{p,q} p^2 a^{\dagger}_{p+q} a^{\dagger}_{-p} (\mathcal{N}_++1)^{-1} a_{-p} a_{p+q} + \delta^{-1} \sum_{p,|q| \leq N^{\alpha}} q^2 |\tilde {\varphi}_p |^2 \, a^{\dagger}_{q} (\mathcal{N}_++1) a_{q}(a^{\dagger}_{0} a_{0}) \\ &\lesssim \delta H_1 + \delta^{-1} N^{1+2\alpha} \|\tilde {\varphi}\|_2^2 (\mathcal{N}_++1)^2. \end{align*} $$

Choosing $\delta = N^{-3\alpha /2}$ , we conclude that

$$\begin{align*}\pm \mathcal{E}_{[H_1,{\mathcal{B}}_3]} \lesssim N^{-3\alpha/2} H_1 + N^{-1+5\alpha/2} (\mathcal{N}_++1)^2. \end{align*}$$

Let us now turn to (68). Recalling (10) and (59), we find

$$ \begin{align*} [Q_4,{\mathcal{B}}_3] &= \frac{1}{2} \sum_{r,p,q} \sum_{m,n} \hat{V}_N (r) \tilde {\varphi}_m \chi_{|n| \leq N^{\alpha}} \big[ a_{p+r}^{\dagger} a_q^{\dagger} a_p a_{q+r}, a_{m+n}^{\dagger} a_{-m}^{\dagger} a_n a_0 \big] + \text{h.c.} \\ &= \frac{1}{2} \sum_{r,p,q} \sum_{m,n} \hat{V}_N (r) \tilde {\varphi}_m \chi_{|n| \leq N^{\alpha}} \\ &\hspace{1.5cm} \times \Big\{ a_{p+r}^{\dagger} a_q^{\dagger} \big[ a_p a_{q+r} , a_{m+n}^{\dagger} a_{-m}^{\dagger} \big] a_n a_0 + a_{m+n}^{\dagger} a_{-m}^{\dagger} \big[ a_{p+r}^{\dagger} a_q^{\dagger} , a_n a_0 \big] a_p a_{q+r} \Big\}. \end{align*} $$

Using equation (6) and rearranging all terms in normal order, we arrive at

$$ \begin{align*} [Q_4,{\mathcal{B}}_3] &= \sum_{r,p,q} \hat V_N(p-r)\tilde {\varphi}_r \chi_{|q| \leq N^{\alpha}} (a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.})+ \mathcal{E}_{[Q_4,{\mathcal{B}}_3]}, \end{align*} $$

where

$$\begin{align*}\begin{aligned} 2 \mathcal{E}_{[Q_4,{\mathcal{B}}_3]}= \; &- \sum_{p,q,m, r} \hat V_N(r) \tilde {\varphi}_m \chi_{|p+r| \leq N^{\alpha}} \, a^{\dagger}_{m+p+r} a^{\dagger}_{-m} a^{\dagger}_q a_{q+r} a_p a_0 \\[4pt] & - \sum_{p,q,m, r} \hat V_N(r) \tilde {\varphi}_m \chi_{|q| \leq N^{\alpha}} \, a^{\dagger}_{m+q} a^{\dagger}_{-m} a^{\dagger}_{p+r} a_{q+r} a_p a_0 \\[4pt] & + \sum_{p,q,m, r} \hat V_N(r) \tilde {\varphi}_m \chi_{|q+r-m| \leq N^{\alpha}} \, a^{\dagger}_{p+r} a^{\dagger}_{-m} a^{\dagger}_{q} a_{q+r-m} a_p a_0 \\[4pt] & + \sum_{p,q,m, r} \hat V_N(r) \tilde {\varphi}_m \chi_{|p-m| \leq N^{\alpha}} \, a^{\dagger}_{p+r} a^{\dagger}_{-m} a^{\dagger}_{q} a_{q+r} a_{p-m} a_0 \\[4pt] & + \sum_{p,q,m, r} \hat V_N(r) \tilde {\varphi}_{-q-r} \chi_{|m| \leq N^{\alpha}} \, a^{\dagger}_{p+r} a^{\dagger}_{q} a^{\dagger}_{-q-r+m} a_{p} a_{m} a_0 \\[4pt] & + \sum_{p,q,m, r} \hat V_N(r) \tilde {\varphi}_{-p} \chi_{|p-m| \leq N^{\alpha}} \, a^{\dagger}_{p+r} a^{\dagger}_{q} a^{\dagger}_{-m} a_{q+r} a_{p-m} a_0 + \text{h.c.} =: \; \sum_{i=1}^6 \mathcal{E}_i. \end{aligned} \end{align*}$$

For a parameter $\delta> 0$ , we find

$$ \begin{align*} \pm \mathcal{E}_1 &\lesssim \delta \|\hat V_N\|_{\infty}^2 \|\chi_{|\cdot|\leq N^{\alpha}}\|_{1} \mathcal{N}_+^{3} + \delta^{-1} \|\tilde {\varphi}\|_2^2 \| \chi_{|\cdot|\leq N^{\alpha}}\|_{1} \mathcal{N}^2_+ \mathcal{N}_0 \, , \\[4pt] &\lesssim N^{-1+5\alpha/2} \mathcal{N}_+^2, \end{align*} $$

where, in the last step, we chose $\delta = N^{-\alpha /2}$ and used $\mathcal {N}_0 \leq N$ (and Lemma 2). To estimate $\mathcal {E}_2, \dots , \mathcal {E}_6$ , we switch to position space. For arbitrary $\delta> 0$ , we find

$$ \begin{align*} \pm \mathcal{E}_2 &= \pm \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) \check{\chi}_{|\cdot|\leq N^{\alpha}}(z-y) a^{\dagger}_z a^{\dagger}(\check{\tilde {\varphi}}^z) a^{\dagger}_x a_x a_y a_0 + \text{h.c.} \\[4pt] &\leq \delta \|\chi_{|\cdot|\leq N^{\alpha}}\|_2^2 \, Q_4 a^{\dagger}_0 a_0 + \delta^{-1} \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) a^{\dagger}_z a^{\dagger}(\check{\tilde {\varphi}}^z) a^{\dagger}_x a_x a(\check{\tilde {\varphi}}^z) a_z \\[4pt] &\lesssim \delta N^{1+3\alpha} Q_4 + \delta^{-1} N^{-3-\alpha} \mathcal{N}_+^3 \\[4pt] &\lesssim N^{-\alpha/2} Q_4 + N^{-2+5\alpha/2} \mathcal{N}_+^3, \end{align*} $$

where, in the last line, we fixed $\delta = N^{-1-7\alpha /2}$ . Similarly, we find

$$ \begin{align*} \pm \mathcal{E}_3 &= \pm \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) a^{\dagger}_xa^{\dagger}_y a^{\dagger}(\check{\tilde {\varphi}}^z) a(\check{\chi}_{|\cdot|\leq N^{\alpha}}^y) a_x a_0 + \text{h.c.} \\[4pt] &\lesssim \delta Q_4 + \delta^{-1} N \|V_N\|_{1} \|\check \chi_{|\cdot|\leq N^{\alpha}}\|_{2}^2 \|\tilde {\varphi}\|_2^2 (\mathcal{N}_++1)^3 \\[4pt] &\lesssim N^{-\alpha/2} Q_4 + N^{-2+5\alpha/2} (\mathcal{N}_+ + 1)^3, \end{align*} $$

taking $\delta = N^{-\alpha /2}$ . Furthermore, for an arbitrary $\xi \in \mathcal {F}$ , we have

$$ \begin{align*} |\langle \xi, \mathcal{E}_4 \xi \rangle| &= \Big| \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) \check{\chi}_{|\cdot|\leq N^{\alpha}}(x-z) \langle \xi, a^{\dagger}_xa^{\dagger}_z a^{\dagger}(\check{\tilde {\varphi}}^x) a_y a_z a_0 \xi \rangle \Big| \\ &\leq \| \check{\chi}_{|.| \leq N^{\alpha}} \|_{\infty} \int _{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N (x-y) \| a_x a_z a (\check{\tilde{\varphi}}^x) \xi \| \| a_y a_z a_0 \xi \| \\ &\lesssim \| \chi_{|.| \leq N^{\alpha}} \|_1 \| \tilde{\varphi} \|_2 \| V_N \|_1 \left[ \langle \xi , \mathcal{N}^3_+ \xi \rangle + \langle \xi, \mathcal{N}_+^2 \mathcal{N}_0 \xi \rangle \right] \\ &\lesssim N^{-1+5\alpha/2} \langle \xi, \mathcal{N}_+^2 \xi \rangle. \end{align*} $$

For $\mathcal {E}_5$ , we estimate

$$ \begin{align*} \pm \mathcal{E}_5 &= \pm \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) \check{\tilde {\varphi}}(y-z) a^{\dagger}_xa^{\dagger}_y a^{\dagger}_z (\mathcal{N}_++1)^{-1/2}(\mathcal{N}_++1)^{1/2} a_x a(\check{\chi}_{|\cdot|\leq N^{\alpha}}^z) a_0 + \text{h.c.} \\ &\leq \delta Q_4 + \delta^{-1} \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N(x-y) |\check{\tilde {\varphi}}(z-y)|^2 a_x^{\dagger} a^{\dagger} (\check{\chi}^z_{|.| \leq N^{\alpha}}) a_0^{\dagger} (\mathcal{N}_+ + 1) a_0 a (\check{\chi}^z_{|.| \leq N^{\alpha}}) a_x \\ &\lesssim \delta Q_4 + \delta^{-1} N \|V_N\|_{1} \|\check \chi_{|\cdot|\leq N^{\alpha}}\|_{2}^2 \|\tilde {\varphi}\|_2^2 (\mathcal{N}_++1)^3 \\ &\lesssim N^{-\alpha/2} Q_4 + N^{-2+5\alpha/2} \mathcal{N}_+^3, \end{align*} $$

again choosing $\delta = N^{-\alpha /2}$ . By a simple change of variable, it is easy to check that $\mathcal {E}_6 = \mathcal {E}_5$ . This concludes the proof of the lemma.

With the bounds from the last lemma, we can estimate the operator $\Gamma _3$ defined in equation (65).

Lemma 13. We have

(70) $$ \begin{align} \begin{aligned} \pm \Gamma_3 \lesssim \; &N^{-3\alpha/2} H_1 + N^{-\alpha/2} Q_4 + N^{3\alpha/2-1/2} (\mathcal{N}_+ + 1)^{3/2} + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{aligned} \end{align} $$

Proof. With Lemma 12, we find, using the scattering equation (14),

$$\begin{align*}\Gamma_3 = \big[ H_1 + Q_4 , {\mathcal{B}}_3 \big] + Q_3 = \tilde{Q}_{3,1} + \tilde{Q}_{3,2} + Q_3^> + \mathcal{E}_{[H_1,{\mathcal{B}}_3]} + \mathcal{E}_{[Q_4,{\mathcal{B}}_3]} \end{align*}$$

with

(71) $$ \begin{align} \tilde Q_{3,1} & = - \sum_{p,q} p^2 {\varphi}_p \chi_{|p| \leq N^{\alpha}} \chi_{|q| \leq N^{\alpha}} \, a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.}, \end{align} $$
(72) $$ \begin{align} \tilde Q_{3,2} &= - \sum_{p,q, r} \hat{V}_N (p-r) {\varphi}_r \chi_{|r| \leq N^{\alpha}} \chi_{|q| \leq N^{\alpha}} \, a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.}, \end{align} $$
(73) $$ \begin{align} Q_3^{>} & = \sum_{p,q} \hat V_N(p) \chi_{|q| > N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 +\text{h.c.} \end{align} $$

It follows easily from Lemma 2 that $\| \chi _{|p| \leq N^{\alpha }} p^2 {\varphi }_p\|_{2}\lesssim N^{-1+3\alpha /2}$ ; thus

$$\begin{align*}\pm \tilde Q_{3,1} \lesssim N^{-1/2 + 3\alpha/2} (\mathcal{N}_++1)^{3/2}. \end{align*}$$

Denoting ${\varphi }^{<}_p = {\varphi }_p \chi _{|p| \leq N^{\alpha }}$ , we write $\tilde {Q}_{3,2}$ in position space as

$$ \begin{align*} \pm \tilde Q_{3,2} &= \pm \int_{\Lambda^3} \mathrm{d}x \mathrm{d}y \mathrm{d}z \, V_N (x-y) \check{{\varphi}}^{<} (x-y) \check{\chi}_{|\cdot|\leq N^{\alpha}}(x-z) a^{\dagger}_x a^{\dagger}_y a_z a_0 + \text{h.c.}\\ &\leq \delta Q_4 + \delta^{-1} N \| V_N^{1/2} \check{\varphi}^< \|_2^2 \|\chi_{|\cdot|\leq N^{\alpha}}\|_{\infty}^2 (\mathcal{N}_++1) \\ &\lesssim N^{-\alpha/2} Q_4 + N^{-2+5\alpha/2} (\mathcal{N}_+ + 1) \, , \end{align*} $$

where we chose $\delta = N^{-\alpha /2}$ and used the estimate

(74) $$ \begin{align} \int_{\Lambda} a^{\dagger} (\check{\chi}_{|.| \leq N^{\alpha}}^x) \, a (\check{\chi}_{|.| \leq N^{\alpha}}^x) \, \mathrm{d}x = d\Gamma (\chi^2_{|p| \leq N^{\alpha}}) \leq \| \chi_{|.| \leq N^{\alpha}} \|_{\infty}^2 \, \mathcal{N}_+ \leq \mathcal{N}_+. \end{align} $$

Proceeding similarly, we find

$$ \begin{align*} \pm Q_3^{>} &= \pm \int_{\Lambda^3} \mathrm{d}x \mathrm{d}y \mathrm{d}z \, V_N (x-y) \check{\chi}_{|\cdot|> N^{\alpha}}(x-z) a^{\dagger}_x a^{\dagger}_y a_z a_0 + \text{h.c.}\\ &\lesssim \delta Q_4 + \delta^{-1} N^{1-2\alpha} \| V_N \|_1 H_1 \lesssim N^{-\alpha/2} Q_4 + N^{-3\alpha/2} H_1, \end{align*} $$

where we took $\delta = N^{-\alpha /2}$ and used that

$$\begin{align*}\int_{\Lambda} a^{\dagger}(\check{\chi}_{|\cdot|> N^{\alpha}}^{x}) a(\check{\chi}_{|\cdot|> N^{\alpha}}^x) \mathrm{d}x = \sum_{|p| > N^{\alpha}} a_p^{\dagger} a_p \leq N^{-2\alpha} H_1. \end{align*}$$

Combining the bounds for $\tilde {Q}_{3,1} , \tilde {Q}_{3,2}, Q_3^>$ with the estimates for $\mathcal {E}_{[H_1,{\mathcal {B}}_3]}, \mathcal {E}_{[Q_4,{\mathcal {B}}_3]}$ from Lemma 12, we obtain the bound (70).

To obtain similar bounds for the integral, we also need a priori control over the growth of $H_1, Q_4$ .

Lemma 14. We have

(75) $$ \begin{align} e^{-s {\mathcal{B}}_3} Q_4 e^{s {\mathcal{B}}_3} &\lesssim Q_4 + \mathcal{N}_+ + 1 + N^{-1+5\alpha/2} (\mathcal{N}_++1)^2, \end{align} $$
(76) $$ \begin{align} e^{-s {\mathcal{B}}_3} H_1e^{s {\mathcal{B}}_3} &\lesssim H_1 + Q_4 + \mathcal{N}_+ +1 + N^{-1+3\alpha} (\mathcal{N}_++1)^2. \end{align} $$

Proof. For arbitrary $\xi \in \mathcal {F}$ , we define $f(s) = \langle \xi , e^{-s{\mathcal {B}}_3} Q_4 e^{s {\mathcal {B}}_3} \xi \rangle $ so that

$$\begin{align*}f' (s) = \langle \xi, e^{-s {\mathcal{B}}_3} [ Q_4 , \mathcal{B}_3 ] e^{s {\mathcal{B}}_3} \xi \rangle. \end{align*}$$

From Lemma 12, we find

$$\begin{align*}[Q_4,{\mathcal{B}}_3] = \sum_{r,p,q} \hat V_N(p-r)\tilde {\varphi}_r \chi_{|q| \leq N^{\alpha}} (a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.})+ \mathcal{E}_{[Q_4,{\mathcal{B}}_3]}, \end{align*}$$

where

$$\begin{align*}\pm \mathcal{E}_{[Q_4,{\mathcal{B}}_3]} \lesssim N^{-\alpha/2} Q_4 + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{align*}$$

Switching to position space, we have

$$\begin{align*}\begin{aligned} \sum_{r,p,q} \hat V_N(p-r) &\tilde {\varphi}_r \chi_{|q| \leq N^{\alpha}} (a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.}) \\ &= \sum_{r,p,q} \hat V_N(p-r)\tilde {\varphi}_r \chi_{|q| \leq N^{\alpha}} (a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0 + \text{h.c.})\\ &= \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y \, V_N (x-y) \check{\tilde{\varphi}} (x-y) a^{\dagger}_x a^{\dagger}_y a( \check{\chi}_{|\cdot| \leq N^{\alpha}}^x) a_0 + \text{h.c.} \\ &\lesssim Q_4 + \int \mathrm{d}x \mathrm{d}y \, V_N(x-y) |\check{\tilde{\varphi}} (x-y)|^2 a_0^{\dagger} a^{\dagger} (\check{\chi}^x_{|.| \leq N^{\alpha}}) a (\check{\chi}^x_{|.| \leq N^{\alpha}}) a_0 \\ &\lesssim Q_4 + N \| V_N \|_1 \| \tilde {\varphi} \|_{\infty}^2 \| \chi_{|.| \leq N^{\alpha}}^2 \|_{\infty} \lesssim Q_4 + \mathcal{N}_+ , \end{aligned} \end{align*}$$

where we used Lemma 2 and argued as in (74). We conclude that

$$\begin{align*}\pm \big[ Q_4 , {\mathcal{B}}_3 \big] \lesssim Q_4 + \mathcal{N}_+ + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{align*}$$

Therefore, using Lemma 3, we find

$$\begin{align*}f' (s) \lesssim f(s) + \langle \xi, \mathcal{N}_+ \xi \rangle + N^{-1+5\alpha/2} \langle \xi , (\mathcal{N}_+ + 1)^2 \xi \rangle. \end{align*}$$

By the Gronwall lemma, we obtain the bound (75).

To prove the estimate (76), we proceed similarly. For $\xi \in \mathcal {F}$ , we define $g(s) = \langle \xi , e^{-s{\mathcal {B}}_3} H_1 e^{s {\mathcal {B}}_3} \xi \rangle $ for any $|s| \leq 1$ , which leads to

$$\begin{align*}g' (s) = \langle \xi , e^{-\mathcal{B}_3} H_1 e^{s\mathcal{B}_3} \xi \rangle. \end{align*}$$

From Lemma 12, we have

$$\begin{align*}\big[ H_1, {\mathcal{B}}_3 \big] = \sum_{p,q,r} \hat{V}_N (p-r) (\delta_{0,r} + {\varphi}_r) \chi_{|p|> N^{\alpha}} \chi_{|q| \leq N^{\alpha}} a_{p+q}^{\dagger} a_{-p}^{\dagger} a_q a_0 + \text{h.c.} + \mathcal{E}_{[H_1,{\mathcal{B}}_3]}, \end{align*}$$

where

$$\begin{align*}\pm \mathcal{E}_{[H_1,{\mathcal{B}}_3]} \leq N^{-3\alpha/2} H_1 + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{align*}$$

Writing $\chi _{|p|> N^{\alpha }} = 1 - \chi _{|p| \leq N^{\alpha }}$ , we decompose

$$\begin{align*}\sum_{p,q,r} \hat{V}_N (p-r) (\delta_{0,r} + {\varphi}_r) \chi_{|p|> N^{\alpha}} \chi_{|q| \leq N^{\alpha}} a_{p+q}^{\dagger} a_{-p}^{\dagger} a_q a_0 + \text{h.c.} = \mathcal{E}_1 + \mathcal{E}_2, \end{align*}$$

where

$$\begin{align*}\begin{aligned} \mathcal{E}_1 &= \pm \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N (x-y) (1 + \check{{\varphi}}) (x-y) a_x^{\dagger} a_y^{\dagger} a(\check{\chi}_{|.| \leq N^{\alpha}}^x) a_0 + \text{h.c.} \\ &\lesssim Q_4 + N \|V_N\|_{1} \|(1+\check {\varphi})\|_{\infty}^2 \|\chi_{|\cdot|\leq N^{\alpha}}\|_{\infty}^2 \mathcal{N}_+ \lesssim Q_4 + \mathcal{N}_+ \end{aligned} \end{align*}$$

and

$$\begin{align*}\begin{aligned} \pm \mathcal{E}_2 &= \pm \sum_{p,q,r} \hat{V}_N (p-r) (\delta_{0,r} + {\varphi}_r ) \chi_{|p| \leq N^{\alpha}} \chi_{|q| \leq N^{\alpha}} a_{p+q}^{\dagger} a_{-p}^{\dagger} a_q a_0 + \text{h.c.} \\ & \lesssim \| V_N \|_1 \| (1+ \check{{\varphi}}) \|_{\infty} \Big[ \delta \mathcal{N}_+^2 + \delta^{-1} N \| \chi_{|.| \leq N^{\alpha}} \|_2^2 \mathcal{N}_+ \Big] \\ &\lesssim N^{-1} \Big[ \delta \mathcal{N}_+^2 + \delta^{-1} N^{1+3\alpha} \mathcal{N}_+ \Big] \lesssim \mathcal{N}_+ + N^{-1+3\alpha} \mathcal{N}_+^2\, , \end{aligned} \end{align*}$$

choosing in the last step $\delta = N^{3\alpha }$ . Thus, with Lemma 3, we find

$$\begin{align*}g' (s) \lesssim f(s) + g(s) + \langle \xi, \mathcal{N}_+ \xi \rangle + N^{-1+3\alpha} \langle \xi, (\mathcal{N}_+ + 1)^2 \xi \rangle. \end{align*}$$

With the estimate (75) and applying Gronwall’s lemma, we obtain the desired estimate (76).

In the next lemma, we control the contribution on the right-hand side of equation (66) arising from the commutator $[Q_3, \mathcal {B}_3]$ .

Lemma 15. We have

$$ \begin{align*} \int_0^1\int_s^1 e^{-t{\mathcal{B}}_3} [Q_3,B_3]e^{t{\mathcal{B}}_3} \mathrm{d}t\mathrm{d}s = 2 (8\pi \mathfrak{a}_N - \hat{V} (0)) \mathcal{N}_+ + \mathcal{E}_{[Q_3,{\mathcal{B}}_3]}, \end{align*} $$

with

$$ \begin{align*}\pm \mathcal{E}_{[Q_3,{\mathcal{B}}_3]} \lesssim N^{-2\alpha} H_1 + N^{-\alpha/2} Q_4 + \big[ N^{-\alpha/2} + N^{-1+\alpha} \big] (\mathcal{N}_+ + 1) + N^{-1+5\alpha/2} (\mathcal{N}_++1)^2. \end{align*} $$

Proof. We compute

(77) $$ \begin{align} \begin{aligned} [Q_3,{\mathcal{B}}_3] &= \sum_{p,q,r,s } \hat V_N(p) \tilde {\varphi}_r \chi_{|s|\leq N^{\alpha}} [a^{\dagger}_{p+q} a^{\dagger}_{-p}a_q a_0 + a^{\dagger}_0 a^{\dagger}_q a_{-p} a_{p+q}, a^{\dagger}_{r+s}a^{\dagger}_{-r} a_sa:0 - a^{\dagger}_0 a^{\dagger}_sa:{-r} a_{r+s}] \\ &=\sum_{p,q,r,s } \hat V_N(p) \tilde {\varphi}_r \chi_{|s|\leq N^{\alpha}} [a^{\dagger}_{p+q} a^{\dagger}_{-p}a_q a_0, a^{\dagger}_{r+s}a^{\dagger}_{-r} a_sa:0] + \text{h.c.} \\ &\quad + \sum_{p,q,r,s } \hat V_N(p) \tilde {\varphi}_r \chi_{|s|\leq N^{\alpha}} [a^{\dagger}_0 a^{\dagger}_sa:{-r} a_{r+s},a^{\dagger}_{p+q} a^{\dagger}_{-p}a_q a_0] + \text{h.c.} =: (\textrm{I}) + (\textrm{II}). \end{aligned} \end{align} $$

We start by estimating term $(\textrm {I})$ . With the canonical commutation relations, we obtain

(78) $$ \begin{align} \begin{aligned} (\textrm{I}) = \; &\sum_{p,q,r}\hat V_N(p) \left[ \tilde {\varphi}_r \chi_{|q-r|\leq N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p}a^{\dagger}_{-r} a_{q-r} + \tilde {\varphi}_q \chi_{|r|\leq N^{\alpha}} a^{\dagger}_{p+q} a^{\dagger}_{-p}a^{\dagger}_{-q+r} a_{r} \right. \\ &\hspace{1cm} - \left. \tilde {\varphi}_r \chi_{|q+p|\leq N^{\alpha}} a^{\dagger}_{r+p+q} a^{\dagger}_{-r}a^{\dagger}_{-p} a_{q} - \tilde {\varphi}_r \chi_{|p|\leq N^{\alpha}} a^{\dagger}_{r-p} a^{\dagger}_{-r}a^{\dagger}_{p+q} a_{q}\right] a_0a_0 + \text{h.c.} \\ =: \; &(\textrm{I})_a + (\textrm{I})_b + (\textrm{I})_c + (\textrm{I})_d. \end{aligned} \end{align} $$

To estimate the first term, we rewrite it in position space. We find

$$ \begin{align*} \pm (\textrm{I})_a &= \pm \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y \, V_N(x-y) a^{\dagger}_x a^{\dagger}_y a^{\dagger}(\check {\tilde {\varphi}}^x) a(\check{ \chi}^x_{|\cdot| \leq N^{\alpha}}) a_0 a_0 + \text{h.c.} \\[3pt] &\leq \delta Q_4 + \delta^{-1} \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y \, V_N(x-y) a^{\dagger} (\check{ \chi}^x_{|\cdot| \leq N^{\alpha}}) a(\check {\tilde {\varphi}}^x)a^{\dagger}(\check {\tilde {\varphi}}^x) a(\check{ \chi}^x_{|\cdot| \leq N^{\alpha}}) a^{\dagger}_0a^{\dagger}_0 a_0 a_0 \\[3pt] &\lesssim \delta Q_4 + \delta^{-1} N^2 \|\tilde {\varphi}\|_2^2 \|V_N\|_1 \|\chi_{|\cdot| \leq N^{\alpha}}\|_{\infty}^2 (\mathcal{N}_++1)^2 \\[3pt] &\lesssim \delta Q_4 + \delta^{-1} N^{-\alpha-1} (\mathcal{N}_++1)^2 \leq N^{-\alpha/2} Q_4 + N^{-1-\alpha/2} (\mathcal{N}_+ + 1)^2, \end{align*} $$

where we used that $\mathcal {N}_0 \leq N$ and the bound (74) and, in the last step, set $\delta = N^{-\alpha /2}$ . The second term in equation (78) is dealt with similarly. We obtain

$$ \begin{align*} \pm (\textrm{I})_b &= \pm \int_{\Lambda^3} \mathrm{d}x \mathrm{d}y \mathrm{d}z \, V_N(x-y) \check {\tilde {\varphi}}(x-z) \, a^{\dagger}_x a^{\dagger}_y a^{\dagger}_z a(\check{ \chi}^z_{|\cdot| \leq N^{\alpha}}) a_0 a_0 + \text{h.c.} \\[3pt] &\leq \delta Q_4 + \delta^{-1}\int_{\Lambda^3} \mathrm{d}x \mathrm{d}y \mathrm{d}z \, V_N(x-y) |\check {\tilde {\varphi}}(x-z)|^2 a^{\dagger}(\check{ \chi}^z_{|\cdot| \leq N^{\alpha}}) (\mathcal{N}_++1)a(\check{ \chi}^z_{|\cdot| \leq N^{\alpha}}) a^{\dagger}_0a^{\dagger}_0a_0a_0 \\[3pt] & \lesssim \delta Q_4 + \delta^{-1} N^2 \|\tilde {\varphi}\|_2^2 \|V_N\|_{1} \|\chi_{|\cdot| \leq N^{\alpha}}\|_{\infty}^2 (\mathcal{N}_++1)^2 \\[3pt] &\lesssim N^{-\alpha/2} Q_4 + N^{-1-\alpha/2} (\mathcal{N}_+ + 1)^2, \end{align*} $$

choosing again $\delta = N^{-\alpha /2}$ . For the third term in equation (78), we bound it, for an arbitrary $\delta> 0$ , with Cauchy-Schwarz by

$$ \begin{align*} \pm (\textrm{I})_c&\lesssim \delta \sum_{p,q,r}|\hat V_N(p)| a^{\dagger}_{r+p+q} a^{\dagger}_{-r}a^{\dagger}_{-p} (\mathcal{N}_+ +1)^{-1}a_{-p} a_{-r} a_{r+p+q} \\ &\qquad + \delta^{-1} \sum_{p,q,r} |\hat V_N(p)| |\tilde {\varphi}_r|^2 \chi_{|q+p|\leq N^{\alpha}} a^{\dagger}_{q}(\mathcal{N}_++1)a_{q}a^{\dagger}_0a^{\dagger}_0a_0a_0 \\ &\lesssim (\delta \|\hat V_N\|_{\infty} + \delta^{-1} N^2 \|\tilde {\varphi}\|_{2}^2 \|\hat V_N \ast \chi_{|\cdot|\leq N^{\alpha}}\|_{\infty})(\mathcal{N}_++1)^2 \\ &\lesssim N^{- 1+5\alpha/2} (\mathcal{N}_+ + 1)^2, \end{align*} $$

where, at the end, we took $\delta = N^{-\alpha }$ and used $\|\hat V_N \ast \chi _{|\cdot |\leq N^{\alpha }}\|_{\infty } \leq \|\hat V_N\|_{\infty } \|\chi _{|\cdot |\leq N^{\alpha }}\|_{1} \lesssim N^{-1+3\alpha }$ . The last term in equation (78) can be bounded, again by Cauchy-Schwarz, by

$$ \begin{align*} \pm (\textrm{I})_d &\lesssim \delta (\mathcal{N}_+ + 1)^2 + \delta^{-1} N^2 \| \hat{V}_N \|_{\infty}^2 \| \tilde{\varphi} \|_2^2 \| \chi_{|.| \leq N^{\alpha}} \|_1 (\mathcal{N}_+ + 1)^2 \lesssim N^{-1+\alpha} (\mathcal{N}_+ + 1)^2, \end{align*} $$

where we used $\delta = N^{\alpha }$ .

Let us now consider term $(\textrm {II})$ in equation (77). We write

(79) $$ \begin{align} \begin{aligned} (\textrm{II}) &= \sum_{p,q,r,s} \hat{V}_N (p) \tilde{\varphi}_r \chi_{|s| \leq N^{\alpha}} \\ &\hspace{1cm} \times \Big\{ a_0^{\dagger} a_s^{\dagger} \big[ a_{-r} a_{r+s} , a^{\dagger}_{p+q} a^{\dagger}_{-p} \big] a_q a_0 + a_{p+q}^{\dagger} a_{-p}^{\dagger} \big[ a_0^{\dagger} a_s^{\dagger} , a_q a_0 \big] a_{-r} a_{r+s} \Big\} + \text{h.c.} \\ &=: (\textrm{II})_a + (\textrm{II})_b. \end{aligned} \end{align} $$

With

$$\begin{align*}a_{p+q}^{\dagger} a_{-p}^{\dagger} \big[ a_0^{\dagger} a_s^{\dagger} , a_q a_0 \big] a_{-r} a_{r+s} = - \delta_{sq} \, a_{p+q}^{\dagger} a_{-p}^{\dagger} a_0^{\dagger} a_0 a_{-r} a_{r+s} - a_{p+q}^{\dagger} a_{-p}^{\dagger} a_q a_s^{\dagger} a_{-r} a_{r+s}, \end{align*}$$

we obtain

$$\begin{align*}\begin{aligned} (\textrm{II})_b = \; &-\sum_{p,q,r} \hat{V}_N (p) \tilde{\varphi}_r \chi_{|q| \leq N^{\alpha}} \, a_{p+q}^{\dagger} a_{-p}^{\dagger} a_0^{\dagger} a_0 a_{-r} a_{r+q} \\ &- \sum_{p,q,r,s} \hat{V}_N (p) \tilde{\varphi}_r \chi_{|s| \leq N^{\alpha}} \, a_{p+q}^{\dagger} a_{-p}^{\dagger} a_s^{\dagger} a_q a_{-r} a_{r+s} \\ &- \sum_{p,r,q} \hat{V}_N (p) \tilde{\varphi}_r \chi_{|q| \leq N^{\alpha}} \, a_{p+q}^{\dagger} a^{\dagger}_{-p} a_{-r} a_{r+q} =: (\textrm{II})_{b1} + (\textrm{II})_{b2} + (\textrm{II})_{b3}. \end{aligned} \end{align*}$$

We can bound $(\textrm {II})_{b3}$ by switching to position space. We find

$$\begin{align*}\begin{aligned} \pm (\textrm{II})_{b3} &= \pm \int_{\Lambda^4} \mathrm{d}x \mathrm{d}y \mathrm{d}u \mathrm{d}v \, V_N (x-y) \check{\chi}_{|.| \leq N^{\alpha}} (x-u) \check{\tilde {\varphi}} (u-v) a_x^{\dagger} a_y^{\dagger} a_u a_v + \text{h.c.} \\ &\lesssim \delta \| \check{\tilde {\varphi}} \|_2^2 \, Q_4 + \delta^{-1} \| \check{\chi}_{|.| \leq N^{\alpha}} \|_2^2 \, (\mathcal{N}_+ + 1)^2 \lesssim N^{-\alpha/2} Q_4 + N^{-3 +5 \alpha/2} (\mathcal{N}_+ + 1)^2. \end{aligned} \end{align*}$$

Term $(\textrm {II})_{b1}$ can be bounded analogously, but it contains an additional factor $\mathcal {N}_0 = a_0^{\dagger } a_0 \leq N$ . Thus

$$\begin{align*}\pm (\textrm{II})_{b1} \leq N^{-\alpha/2} Q_4 + N^{-1 +5 \alpha/2} (\mathcal{N}_+ + 1)^2. \end{align*}$$

Term $(\textrm {II})_{b1}$ can also be bounded in position space. We obtain

$$\begin{align*}\begin{aligned} \pm (\textrm{II})_{b2} &= \pm \int_{\Lambda^3} \mathrm{d}x \mathrm{d}y \mathrm{d}u \, V_N (x-y) \, a_x^{\dagger} a_y^{\dagger} a^{\dagger} (\check{\chi}_{|.| \leq N^{\alpha}}^u) a_x a (\check{\tilde {\varphi}}^u) a_u + \text{h.c.} \\ &\leq \delta \| \check{\chi}_{|.| \leq N^{\alpha}} \|_2^2 \, Q_4 \mathcal{N}_+ + \delta^{-1} \| V_N \|_1 \| \check{\tilde {\varphi}} \|_2^2 (\mathcal{N}_+ + 1)^3 \lesssim N^{-\alpha/2} Q_4 + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2, \end{aligned} \end{align*}$$

where we chose $\delta = N^{-1-7\alpha /2}$ .

Let us now consider term $(\textrm {II})_a$ , defined on the right-hand side of equation (79). With

$$\begin{align*}\begin{aligned} a_0^{\dagger} &a_s^{\dagger} \big[ a_{-r} a_{r+s} , a^{\dagger}_{p+q} a^{\dagger}_{-p} \big] a_q a_0 \\ &= a_0^{\dagger} a_s^{\dagger} \Big\{ \delta_{-r,p+q} a_{-p}^{\dagger} a_{r+s} + \delta_{r,p} a^{\dagger}_{p+q} a_{r+s} + \delta_{r+s, p+q} a_{-r} a_{-p}^{\dagger} + \delta_{r+s,-p} a_{-r} a_{p+q}^{\dagger} \Big\} a_q a_0, \end{aligned} \end{align*}$$

we obtain, rearranging the terms in normal order (with appropriate changes of variables),

(80) $$ \begin{align} \begin{aligned} (\textrm{II})_a = \; &\sum_{p,q} \big( \hat{V}_N (p) + \hat{V}_N (p-q) \big) \tilde{\varphi}_p \chi_{|q| \leq N^{\alpha}} a_0^{\dagger} a_q^{\dagger} a_q a_0 \\ &+ 4 \sum_{p,q,s} \hat{V}_N (p) \tilde{\varphi}_{p+q} \chi_{|s| \leq N^{\alpha}} \, a_0^{\dagger} a_s^{\dagger} a_{-p}^{\dagger} a_{-p-q+s} a_q a_0 + \text{h.c.} \\ =: \; &2 \sum_{p,q} \big( \hat{V}_N (p) + \hat{V}_N (p-q) \big) \tilde{\varphi}_p \chi_{|q| \leq N^{\alpha}} a_0^{\dagger} a_q^{\dagger} a_q a_0 + (\textrm{II})_{a1}, \end{aligned} \end{align} $$

where we can bound, with $\mathcal {N}_0 = a_0^{\dagger } a_0 \leq N$ ,

$$\begin{align*}\begin{aligned} \pm (\textrm{II})_{a1} \lesssim N \| \hat{V}_N \|_{\infty} \Big[ \delta \| \tilde {\varphi} \|_2^2 + \delta^{-1} \| \chi_{|.| \leq N^{\alpha}} \|_2^2 \Big] \mathcal{N}_+^2 \lesssim N^{-1+\alpha} (\mathcal{N}_+ + 1)^2, \end{aligned} \end{align*}$$

choosing $\delta = N^{1-\alpha }$ . Collecting all the estimates we have proved so far, we conclude from equation (77) that

$$\begin{align*}\big[ Q_3, \mathcal{B}_3 \big] = 2 \sum_{p,q} \big( \hat{V}_N (p) + \hat{V}_N (p-q) \big) \tilde{\varphi}_p \chi_{|q| \leq N^{\alpha}} a_0^{\dagger} a_q^{\dagger} a_q a_0 + \mathcal{E}_1, \end{align*}$$

where

(81) $$ \begin{align} \pm \mathcal{E}_1 \lesssim N^{-\alpha/2} Q_4 + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{align} $$

At the expense of adding an additional small error to the right-hand side of the estimate (81), in the main term, we can replace $a_0^{\dagger } a_0 = N - \mathcal {N}_+$ by a factor of N, since

$$\begin{align*}\pm \sum_{p,q} (\hat{V}_N (p) + \hat{V}_N (p-q) ) \tilde{\varphi}_p \chi_{|q| \leq N^{\alpha}} a_q^{\dagger} \mathcal{N}_+ a_q \lesssim N^{-1} \| \tilde {\varphi} \|_1 (\mathcal{N}_++1)^2 \lesssim N^{-1} (\mathcal{N}_+ + 1)^2. \end{align*}$$

Moreover, from

(82) $$ \begin{align} \pm N\sum_{p,q} &\left( \hat V_N(p-q)- \hat V_N(p)\right)\tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} a^{\dagger}_q a_q \notag \\ &\qquad\qquad\lesssim \sum_{p,q} | q/N| \| \nabla \hat V\|_{\infty} |\tilde {\varphi}_p| \chi_{|q| \leq N^{\alpha}} a^{\dagger}_q a_q\lesssim \|\tilde {\varphi}\|_1 N^{-1+\alpha} \mathcal{N}_+ \lesssim N^{-1+\alpha} \mathcal{N}_+, \end{align} $$

we arrive at

$$\begin{align*}\big[ Q_3, \mathcal{B}_3 \big] = 4N \sum_{p} \hat{V}_N (p) \tilde{\varphi}_p \sum_{|q| \leq N^{\alpha}} a_q^{\dagger} a_q + \mathcal{E}_2, \end{align*}$$

where

$$\begin{align*}\pm \mathcal{E}_2 \lesssim N^{-\alpha/2} Q_4 + N^{-1+\alpha} \mathcal{N}_+ + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{align*}$$

Conjugating with $e^{t {\mathcal {B}}_3}$ and integrating over $t\ \mathrm{and}\ s$ , we obtain, with the help of Lemma 10 and Lemma 14 (and the observation that the first estimate in Lemma 10 also holds, if we replace $\mathcal {N}_+$ on the left-hand side by $\sum _{|p| \leq N^{\alpha }} a_p^{\dagger } a_p$ ),

$$\begin{align*}\int_0^1 \int_s^1 e^{-t{\mathcal{B}}_3} [Q_3, {\mathcal{B}}_3] e^{t {\mathcal{B}}_3} \mathrm{d}t \mathrm{ d}s = 2N \sum_q \hat{V}_N (q) \tilde{\varphi}_q \sum_{|p| \leq N^{\alpha}} a_p^{\dagger} a_p + \mathcal{E}_3, \end{align*}$$

where

$$\begin{align*}\pm \mathcal{E}_3 \lesssim N^{-\alpha/2} (Q_4 + \mathcal{N}_+ + 1) + N^{-1+\alpha} \mathcal{N}_+ + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{align*}$$

The claim now follows from equation (19) and the observation that

$$\begin{align*}\mathcal{N}_+ - \sum_{|p| \leq N^{\alpha}} a_p^{\dagger} a_p = \sum_{|p| \geq N^{\alpha}} a_p^{\dagger} a_p \leq N^{-2\alpha} H_1.\\[-47pt] \end{align*}$$

Finally, we consider the conjugation of the operator $\tilde {Q}_2$ , defined in equation (25).

Lemma 16. We have

(83) $$ \begin{align} e^{-{\mathcal{B}}_3} \tilde Q_2 e^{{\mathcal{B}}_3}= \tilde{Q}_2 + \mathcal{E}_{\tilde Q_2} \end{align} $$

with

$$\begin{align*}\pm \mathcal{E}_{\tilde Q_2} \lesssim N^{-1/2+\alpha} (\mathcal{N}_++1)^{3/2} + N^{-3/2+\alpha} (\mathcal{N}_++ 1)^{5/2}. \end{align*}$$

Proof. We have

$$\begin{align*}\begin{aligned} e^{-{\mathcal{B}}_3} \tilde Q_2 e^{{\mathcal{B}}_3} - \tilde Q_2 &= \int_0^1 e^{-s{\mathcal{B}}_3}[ \tilde Q_2,{\mathcal{B}}_3] e^{s{\mathcal{B}}_3}\mathrm{d}s \\ &= \frac{4\pi \mathfrak{a}_N}{N} \sum_{|r| \leq N^{\alpha}} \int_0^1 e^{-s{\mathcal{B}}_3}[ a^{\dagger}_r a^{\dagger}_{-r} a_0 a_0 + \text{h.c.}, {\mathcal{B}}_3] e^{s{\mathcal{B}}_3}\mathrm{d}s. \end{aligned} \end{align*}$$

We compute the commutator

(84) $$ \begin{align} [a^{\dagger}_r a^{\dagger}_{-r} a_0 a_0 &+ \text{h.c.} , a_{p+q}^{\dagger} a_{-p}^{\dagger} a_q a_0 - \text{h.c.} ] \end{align} $$
(85) $$ \begin{align} &=[a^{\dagger}_r a^{\dagger}_{-r}a_0a_0,a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0]+[a^{\dagger}_0 a^{\dagger}_0 a_r a_{-r},a^{\dagger}_{p+q} a^{\dagger}_{-p} a_q a_0]+\text{h.c.}\nonumber \\ &= 2\Big\{ a^{\dagger}_0 a^{\dagger}_0 \big( \delta_{p+q,r} a^{\dagger}_{-p} a_{-p-q} +\delta_{p,r} a^{\dagger}_{p+q} a_{p} \big) a_q a_0 \nonumber\\ &\qquad\qquad\qquad -\delta_{r,q}a^{\dagger}_{p+q} a^{\dagger}_{-p} a^{\dagger}_{-q} a_0 a_0 a_0- a^{\dagger}_{p+q} a^{\dagger}_{-p} a^{\dagger}_{0} a_{q} a_r a_{-r}\Big\} +\text{h.c.} \end{align} $$

Hence, we obtain

$$\begin{align*}\begin{aligned} \frac{4\pi \mathfrak{a}_N}{N} \sum_{|r| \leq N^{\alpha}} \big[ a_r^{\dagger} &a_{-r}^{\dagger} a_0 a_0 + \text{h.c.} , {\mathcal{B}}_3 \big] \\ =\; &\frac{8\pi \mathfrak{a}_N}{N} \sum_{p,q : |p+q| \leq N^{\alpha}} \tilde{\varphi}_p \chi_{|q| \leq N^{\alpha}} a_0^{\dagger} a_0^{\dagger} a_{-p}^{\dagger} a_{-p-q} a_q a_0 \\ &- \frac{4\pi \mathfrak{a}_N}{N} \sum_{p, |q| \leq N^{\alpha}} \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} a_{p+q}^{\dagger} a_{-p}^{\dagger} a_{-q}^{\dagger} a_0 a_0 a_0 \\ &- \frac{4\pi \mathfrak{a}_N}{N} \sum_{p, q, |r| \leq N^{\alpha}} \tilde {\varphi}_p \chi_{|q| \leq N^{\alpha}} a_{p+q}^{\dagger} a_{-p}^{\dagger} a_{0}^{\dagger} a_q a_r a_{-r} =: (\textrm{I}) + (\textrm{II}) + (\textrm{III}). \end{aligned} \end{align*}$$

With Cauchy-Schwarz and using the bounds from Lemma 2, we can bound

$$\begin{align*}\begin{aligned} \pm (\textrm{I}) &\lesssim N^{-1/2-\alpha/2} (\mathcal{N}_+ + 1)^{3/2} , \\ \pm (\textrm{II}) &\lesssim N^{-1/2+\alpha} (\mathcal{N}_+ + 1)^{3/2} , \\ \pm (\textrm{III}) &\lesssim N^{-3/2+\alpha} (\mathcal{N}_+ + 1)^{5/2}. \end{aligned} \end{align*}$$

The claim now follows with Lemma 10.

We are now ready to conclude the proof of Proposition 11.

Proof of Proposition 11.

Recall from equation (66) that

$$\begin{align*}\begin{aligned} e^{-{\mathcal{B}}_3} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} = \; & 4 \pi \mathfrak{a}_N (N-1) + \sum_{|p|\leq N^{\alpha}} \frac{(4\pi \mathfrak{a}_N)^2}{p^2} + H_1 + Q_4 + e^{-{\mathcal{B}}_3}(\tilde H_2 + \tilde Q_2 + \mathcal E_{{\mathcal{B}}_2})e^{{\mathcal{B}}_3} \\ &+ \int_0^1 e^{-t{\mathcal{B}}_3} \Gamma_3 e^{t{\mathcal{B}}_3} \mathrm{d}t + \int_0^1 \int_s^1 e^{-t{{\mathcal{B}}_3}}[Q_3,{{\mathcal{B}}_3}]e^{t{{\mathcal{B}}_3}} \mathrm{d}t \mathrm{d}s, \end{aligned} \end{align*}$$

where

$$\begin{align*}\pm \mathcal{E}_{{\mathcal{B}}_2} \leq N^{-\alpha/2} Q_4 + \big[ N^{-\alpha/2} + N^{-1+5\alpha/2} \big] (\mathcal{N}_+ + 1) + N^{-1+\alpha} \mathcal{N}_+^2 + N^{-2} H_1. \end{align*}$$

With Lemma 10 and Lemma 14, this also implies that

$$\begin{align*}\begin{aligned} \pm e^{-{\mathcal{B}}_3} \mathcal{E}_{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} \lesssim \; &N^{-\alpha/2} Q_4 + N^{-2} H_1 + \big[ N^{-\alpha/2} + N^{-1+5\alpha/2} \big] (\mathcal{N}_+ + 1) \\ &+ N^{-1+2\alpha} (\mathcal{N}_+ + 1)^2. \end{aligned} \end{align*}$$

Applying the first bound in Lemma 10 to the operator $\tilde H_2 = (2\hat {V} (0) - 8 \pi \mathfrak {a}_N) \mathcal {N}_+$ , defined in equation (24), we obtain

$$\begin{align*}\pm \Big[ e^{-{\mathcal{B}}_3} \tilde H_2 e^{{\mathcal{B}}_3} - \tilde H_2 \Big] \lesssim N^{-\alpha/2} (\mathcal{N}_+ + 1). \end{align*}$$

Combining Lemma 13 with Lemma 14, we obtain

$$\begin{align*}\begin{aligned} \int_0^1 e^{-t {\mathcal{B}}_3} &\Gamma_3 e^{t {\mathcal{B}}_3} \mathrm{d}t \\ \lesssim \; & N^{-3\alpha/2} H_1 + N^{-\alpha/2} Q_4 + N^{-\alpha/2} (\mathcal{N}_+ + 1) \\ &+ N^{3\alpha/2-1/2} (\mathcal{N}_+ + 1)^{3/2} + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{aligned} \end{align*}$$

Together with the bounds in Lemmas 15 and 16, and with the observation that

$$\begin{align*}8\pi \mathfrak{a}_N \sum_{|p|> N^{\alpha}} a_p^{\dagger} a_p \leq N^{-2\alpha} H_1, \end{align*}$$

we conclude the proof of Proposition 11.

5 Diagonalisation of a quadratic Hamiltonian

From Proposition 11, we find

(86) $$ \begin{align} \begin{aligned} e^{-{\mathcal{B}}_3} &e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} \\ &= 4 \pi \mathfrak{a}_N (N-1) + \frac{1}{4}\sum_{|p|\leq N^{\alpha}} \frac{(8\pi \mathfrak{a}_N)^2}{p^2} + \sum_{|p|> N^{\alpha}} p^2 a_p^{\dagger} a_p + Q_4 \\ &\quad + \sum_{|p| \leq N^{\alpha}} (p^2 + 8\pi \mathfrak{a}_N ) a^{\dagger}_p \frac{a_0 a_0^{\dagger}}{N} a_p + \frac{1}{2}\sum_{|p| \leq N^{\alpha}} 8 \pi \mathfrak{a}_N [a^{\dagger}_p a^{\dagger}_{-p} \frac{a_0 a_0 }N+\text{h.c.}] + \mathcal{E}, \end{aligned} \end{align} $$

with an error $\mathcal {E}$ satisfying the bound (64). Here we used the observation that, on the sector $\{ \mathcal {N} = N \}$ , we can write

$$\begin{align*}\begin{aligned} \sum_{|p| \leq N^{\alpha}} (p^2 + 8\pi \mathfrak{a}_N) a_p^{\dagger} a_p &= \sum_{|p| \leq N^{\alpha}} (p^2 + 8\pi \mathfrak{a}_N) a_p^{\dagger} \frac{a_0^{\dagger} a_0 + \mathcal{N}_+ + 1}{N} a_p \\ &= \sum_{|p| \leq N^{\alpha}} (p^2 + 8\pi \mathfrak{a}_N) a_p^{\dagger} \frac{a_0 a_0^{\dagger}}{N} a_p + \frac{1}{N} \sum_{|p| \leq N^{\alpha}} (p^2 + 8\pi \mathfrak{a}_N) a_p^{\dagger} \, \mathcal{N}_+ a_p \end{aligned} \end{align*}$$

where the term

$$\begin{align*}\frac{1}{N} \sum_{|p| \leq N^{\alpha}} (p^2 + 8\pi \mathfrak{a}_N) a_p^{\dagger} \, \mathcal{N}_+ a_p \lesssim N^{2\alpha-1} (\mathcal{N}_+ + 1)^2 \end{align*}$$

can be absorbed on the right-hand side of the estimate (64).

In this section, we will diagonalise the operator on the last line of equation (86). Inspired by Bogoliubov theory (on states with $a_0, a_0^{\dagger } \sim \sqrt {N}$ , this operator is approximately quadratic), we define, for $|p| \leq N^{\alpha }$ , the coefficients

$$\begin{align*}\tau_p = - \frac{1}{4} \log \left[ 1 + \frac{16\pi \mathfrak{a}_N}{p^2} \right] \end{align*}$$

so that

$$\begin{align*}\tanh (2 \tau_p) = - \frac{8\pi\mathfrak{a}_N}{p^2 + 8\pi \mathfrak{a}_N}. \end{align*}$$

We also introduce the notation $\gamma _p = \cosh \tau _p$ and $\nu _p = \sinh \tau _p$ .

Lemma 17. We have the pointwise bound $\gamma _p \lesssim 1$ and $\tau _p, \nu _p \lesssim \chi _{|p| \leq N^{\alpha }} /p^2$ . Moreover,

$$\begin{align*}\begin{aligned} \| \tau \|_{\infty} &\leq \| \tau \|_2 \lesssim 1 , \qquad \| \nu \|_{\infty} \leq \| \nu \|_2 \lesssim 1, \qquad \| \gamma - 1 \|_{\infty} \leq \| \gamma - 1 \|_2 \lesssim 1 \end{aligned} \end{align*}$$

and

$$\begin{align*}\begin{aligned} \| \check{\tau} \|_{\infty} &\leq \| \tau \|_1 \lesssim N^{\alpha} , \qquad \| \check{\nu} \|_{\infty} \leq \| \nu \|_1 \lesssim N^{\alpha}. \end{aligned} \end{align*}$$

With these coefficients, we can write

$$\begin{align*}\begin{aligned} \sum_{|p| \leq N^{\alpha}} (p^2 &+ 8\pi \mathfrak{a}_N ) a^{\dagger}_p \frac{a_0 a_0^{\dagger}}{N} a_p + \frac{1}{2}\sum_{|p| \leq N^{\alpha}} 8 \pi \mathfrak{a}_N [a^{\dagger}_p a^{\dagger}_{-p} \frac{a_0 a_0 }N+\text{h.c.}] \\ = \; & \sum_{|p| \leq N^{\alpha}} \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} \, \big(\gamma_p a_p^{\dagger} \frac{a_0}{\sqrt{N}} + \nu_p \frac{a_0^{\dagger}}{\sqrt{N}} a_{-p} \big) \big( \gamma_p \frac{a_0^{\dagger}}{\sqrt{N}} a_p + \nu_p a_{-p}^{\dagger} \frac{a_0}{\sqrt{N}} \big) \\ &-\frac{1}{2} \sum_{|p| \leq N^{\alpha}} \big[ p^2 + 8\pi \mathfrak{a}_N - \sqrt{|p|^4 + 16\pi \mathfrak{a}_N p^2} \big] + \delta \end{aligned} \end{align*}$$

with an error $\delta $ satisfying

$$\begin{align*}\pm \delta \lesssim N^{\alpha-1} (\mathcal{N}_+ + 1). \end{align*}$$

Here, we used the relations

$$\begin{align*}\begin{aligned} \gamma^2_p + \nu_p^2 &= \cosh (2\tau_p) = \frac{1}{\sqrt{1- \tanh^2 (2\tau_p)}} = \frac{p^2 + 8\pi \mathfrak{a}_N}{\sqrt{|p|^4 + 16\pi \mathfrak{a}_N p^2}}\, , \\ 2\gamma_p \nu_p &= \sinh (2\tau_p) = \frac{\tanh (2\tau_p)}{\sqrt{1- \tanh^2 (2\tau_p)}} = \frac{8\pi \mathfrak{a}_N}{\sqrt{|p|^4 + 16\pi \mathfrak{a}_N p^2}} \, , \\ \nu^2_p &= \frac{1}{2} \big[ \cosh (2\tau_p) -1 \big] = \frac{1}{2} \frac{p^2 + 8\pi \mathfrak{a}_N - \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2}}{\sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2}} \, \end{aligned} \end{align*}$$

and the commutator

$$\begin{align*}\big[ \frac{a_0^{\dagger}}{\sqrt{N}} a_{-p} , a_{-p}^{\dagger} \frac{a_0}{\sqrt{N}} \big] = \frac{1}{N} a_0^{\dagger} a_0 - \frac{1}{N} a_{-p}^{\dagger} a_{-p} = 1 - \frac{1}{N} (\mathcal{N}_+ + a_{-p}^{\dagger} a_{-p}). \end{align*}$$

The contribution proportional to $N^{-1}$ on the right-hand side of the last equation produces (using Lemma 17) the error $\delta $ . Inserting in equation (86), we conclude that

(87) $$ \begin{align} \begin{aligned} &e^{-{\mathcal{B}}_3} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} \\ &= 4 \pi \mathfrak{a}_N (N-1) - \frac{1}{2} \sum_{|p| \leq N^{\alpha}} \Big[ p^2 + 8\pi \mathfrak{a}_N - \sqrt{|p|^4 + 16\pi \mathfrak{a}_N p^2} - \frac{(8 \pi \mathfrak{a}_N)^2}{2p^2} \Big] \\ &\quad + \sum_{|p| \leq N^{\alpha}} \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} \, \big(\gamma_p b_p^{\dagger} + \nu_p b_{-p} \big) \big( \gamma_p b_p + \nu_p b_{-p}^{\dagger} \big) + \sum_{|p|> N^{\alpha}} p^2 a_p^{\dagger} a_p + Q_4 + \mathcal{E} , \end{aligned} \end{align} $$

where $\mathcal {E}$ still satisfies the estimate (64) and where we introduced the modified creation and annihilation operators

(88) $$ \begin{align} b_p = \frac{a_0^{\dagger}}{\sqrt{N}} a_p, \quad b_p^{\dagger} = a_p^{\dagger} \frac{a_0}{\sqrt{N}} \end{align} $$

satisfying the commutation relations

(89) $$ \begin{align} \big[ b_p, b_q\big] = \big[ b_p^{\dagger} , b_q^{\dagger} \big] = 0 , \qquad \big[ b_p , b_q^{\dagger} \big] = \delta_{p,q} \big( 1- \frac{\mathcal{N}_+}{N} \big) - \frac{1}{N} a_q^{\dagger} a_p \end{align} $$

and $[ a_p^{\dagger } a_r , b^{\dagger }_q ] = \delta _{r,q} b_p^{\dagger }$ , $[a_p^{\dagger } a_r, b_q ] = - \delta _{p,q} b_r$ . On states with few excitations $a_0, a_0^{\dagger } \simeq \sqrt {N}$ , we have $b_p^{\dagger } \simeq a_p^{\dagger }$ , $b_p \simeq a_p$ . According to Bogoliubov theory, we can therefore expect that the operators $(\gamma _p b_p^{\dagger }. + \nu _p b_{-p})$ and $(\gamma _p b_p + \nu _p b_{-p}^{\dagger })$ can be rotated back to $b_p^{\dagger }$ and, respectively, $b_p$ through conjugation of the Hamiltonian with the unitary transformation generated by the antisymmetric operator

$$\begin{align*}{\mathcal{B}}_4 =\frac{1}{2} \sum_{|p| \leq N^{\alpha}} \tau_p \big( b_p^{\dagger} b_{-p}^{\dagger} - b_p b_{-p} \big) = \frac{1}{2} \sum_{|p| \leq N^{\alpha}} \tau_p \big( a_p^{\dagger} a_{-p}^{\dagger} \frac{a_0 a_0}{N} - \text{h.c.} \big). \end{align*}$$

Notice that ${\mathcal {B}}_4$ has the same form as the operator ${\mathcal {B}}_2$ defined in equation (13) (with a different choice of the coefficients, of course; here it is more convenient to keep the factor $N^{-1}$ out of $\tau _p$ ). To control the action of ${\mathcal {B}}_4$ , we will need rough a priori bounds on the growth of the number and the energy of the excitations.

Lemma 18. For every $k \in \mathbb {N}$ , we haveFootnote 1

(90) $$ \begin{align} e^{-{\mathcal{B}}_4} (\mathcal{N}_++1)^k e^{{\mathcal{B}}_4} &\lesssim (\mathcal{N}_++1)^k. \end{align} $$

Moreover,

(91) $$ \begin{align} e^{-{\mathcal{B}}_4} H_1 e^{{\mathcal{B}}_4} &\lesssim H_1 + N^{\alpha}, \end{align} $$
(92) $$ \begin{align} e^{-{\mathcal{B}}_4} Q_4 e^{{\mathcal{B}}_4} &\lesssim Q_4 + N^{-1+2\alpha} + N^{-1} (\mathcal{N}_+ + 1)^2. \end{align} $$

Proof. The proof of the bound (90) is standard (based on Gronwall’s lemma and the bounds in Lemma 17). To prove equation (91), we define $g(s) = e^{-s{\mathcal {B}}_4} H_1 e^{s{\mathcal {B}}_4}$ and compute (using the commutation relations after equation (89))

$$ \begin{align*} g'(s) &= e^{-s{\mathcal{B}}_4} [H_1,{\mathcal{B}}_4] e^{s{\mathcal{B}}_4} = \sum_{|p| \leq N^{\alpha}} p^2 \tau_p e^{-s{\mathcal{B}}_4} \big[ b^{\dagger}_p b^{\dagger}_{-p} + \text{h.c.} \big] e^{s{\mathcal{B}}_4} \lesssim e^{-s {\mathcal{B}}_4} H_1 e^{s {\mathcal{B}}_4} + \sum_{|p| \leq N^{\alpha}} p^2 \tau_p^2. \end{align*} $$

From Lemma 17, we have $\sum _p p^2 \tau _p^2 \lesssim N^{\alpha }$ ; with Gronwall’s lemma, we obtain the bound (91).

To show the bound (92), we set $h(s) = e^{-s{\mathcal {B}}_4} Q_4 e^{s{\mathcal {B}}_4}$ , and then

(93) $$ \begin{align} h' (s) = e^{-s{\mathcal{B}}_4} [Q_4, {\mathcal{B}}_4] e^{s{\mathcal{B}}_4}. \end{align} $$

Proceeding as in (38), we find (we use here the convention that $\tau _q = 0$ , for $|q|> N^{\alpha }$ )

$$\begin{align*}\begin{aligned} [Q_4, {\mathcal{B}}_4 ] = \; & \frac{1}{2} \sum_{q} (\hat{V}_N * \tau) (q) b_q^{\dagger} b_{-q}^{\dagger} + \sum_{p,q,s} \hat{V}_N (q-p) \tau_s b^{\dagger}_{p+s-q} b^{\dagger}_q a^{\dagger}_{-s} a_p +\text{h.c.} \end{aligned} \end{align*}$$

Switching to position space, we write

$$\begin{align*}[Q_4, {\mathcal{B}}_4 ] = \frac{1}{2} \int_{\Lambda^2} \mathrm{d}x \mathrm{d}y V_N (x-y) \check{\tau} (x-y) b_x^{\dagger} b_y^{\dagger} + \int_{\Lambda^3} \mathrm{d}x \mathrm{d}y \mathrm{d}z V_N (x-y) \check{\tau} (x-z) b_x^{\dagger} b_y^{\dagger} a_z^{\dagger} a_y + \text{h.c.} \end{align*}$$

With the bounds from Lemma 17, we conclude

$$\begin{align*}\pm [Q_4, {\mathcal{B}}_4] \lesssim Q_4 + N^{-1+2\alpha} + N^{-1} (\mathcal{N}_+ + 1)^2. \end{align*}$$

Inserting this in equation (93), applying the bound (90) and then Gronwall’s lemma, we obtain the desired bound (92).

We are now ready to state the main result of this section, which shows that conjugation with $e^{\mathcal {B}_4}$ diagonalises the quadratic part of the Hamilton operator.

Proposition 19. We have

(94) $$ \begin{align} \begin{aligned} e^{-{\mathcal{B}}_4} e^{-{\mathcal{B}}_3} &e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3}e^{{\mathcal{B}}_4} \\ = &\; 4 \pi \mathfrak{a}_N (N-1) + \frac{1}{2}\sum_{p} \left[ \sqrt{p^4 + 16 \pi \mathfrak{a}_N p^2} - p^2 - 8 \pi \mathfrak{a}_N + \frac{(8\pi \mathfrak{a}_N)^2}{2p^2}\right] \\ &+ \sum_{p} \sqrt{p^4 + 16 \pi \mathfrak{a}_N p^2} \, a^{\dagger}_p a_p + e^{-{\mathcal{B}}_4}Q_4e^{{\mathcal{B}}_4} + \mathcal E_{{\mathcal{B}}_4}, \end{aligned} \end{align} $$

where

(95) $$ \begin{align} \begin{aligned} \pm \mathcal{E}_{{\mathcal{B}}_4} \lesssim \; &N^{-3\alpha/2} (H_1 +N^{\alpha}) + N^{-\alpha/2} e^{-{\mathcal{B}}_4} Q_4 e^{{\mathcal{B}}_4} + N^{-\alpha/2} (\mathcal{N}_+ + 1) \\ &+ N^{(3\alpha-1)/2} (\mathcal{N}_+ + 1)^{3/2} + N^{-1+5\alpha/2} (\mathcal{N}_+ + 1)^2. \end{aligned} \end{align} $$

Proof. For $s \in [0;1]$ , we define

$$\begin{align*}E(s) = \sum_{|p| \leq N^{\alpha}} \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} \, \big(\gamma^{s}_p b_p^{\dagger} + \nu_p^{s} b_{-p} \big) \big( \gamma_p^{s} b_p + \nu^{s}_p b_{-p}^{\dagger} \big) \end{align*}$$

with the operators $b_p, b^{\dagger }_p$ defined in equation (88) and with the notation $\gamma _p^{s} = \cosh (s \tau _p)$ and $\nu _p^{s} = \sinh (s \tau _p)$ . In particular, for $s=1$ , this is exactly the operator appearing on the third line in equation (87). For $\psi $ in the sector $\{ \mathcal {N} = N \}$ , we define $f_{\psi } :[0;1] \to \mathbb {R}$ by

$$\begin{align*}f_{\psi} (s) = \langle \psi , e^{-s{\mathcal{B}}_4} E(s) e^{s{\mathcal{B}}_4} \psi \rangle. \end{align*}$$

The idea is that the generalised Bogoliubov transformation $e^{s{\mathcal {B}}_4}$ approximately cancels (on states with few excitations) the symplectic rotations determined by the coefficients $\gamma _p^{s}$ , $\nu _p^{s}$ (it would precisely cancel them if the operators $b_p^{\dagger }, b_p$ satisfied canonical commutation relations); hence, on states with few excitations, we expect $f_{\psi }$ to be approximately constant in s. More precisely, we claim that

(96) $$ \begin{align} \big| f^{\prime}_{\psi} (s) \big| \lesssim N^{2\alpha-1} \langle \psi, (\mathcal{N}_+ + 1)^2 \psi \rangle. \end{align} $$

Assuming for a moment that the bound (96) holds true, we could conclude, integrating over $s \in [0;1]$ , that

$$\begin{align*}e^{-{\mathcal{B}}_4} E(1) e^{{\mathcal{B}}_4} = E(0) + \delta \end{align*}$$

with

$$\begin{align*}\pm \delta \lesssim N^{2\alpha-1} (\mathcal{N}_+ + 1)^2. \end{align*}$$

With the bounds from Lemma 18 (and noticing that the action of ${\mathcal {B}}_4$ on the high-momenta part of the kinetic energy is trivial), this would imply that

(97) $$ \begin{align} \begin{aligned} e^{-{\mathcal{B}}_4} &e^{-{\mathcal{B}}_3} e^{-{\mathcal{B}}_2} H_N e^{{\mathcal{B}}_2} e^{{\mathcal{B}}_3} e^{{\mathcal{B}}_4} \\ &= 4 \pi \mathfrak{a}_N (N-1) - \frac{1}{2} \sum_{|p| \leq N^{\alpha}} \Big[ p^2 + 8\pi \mathfrak{a}_N - \sqrt{|p|^4 + 16\pi \mathfrak{a}_N p^2} - \frac{(8 \pi \mathfrak{a}_N)^2}{2p^2} \Big] \\ &\quad + \sum_{|p| \leq N^{\alpha}} \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} \, a_p^{\dagger} \frac{a_0 a_0^{\dagger}}{N} a_p + \sum_{|p|> N^{\alpha}} p^2 a_p^{\dagger} a_p + e^{-{\mathcal{B}}_4} Q_4 e^{{\mathcal{B}}_4} + \mathcal{E} , \end{aligned} \end{align} $$

where $\mathcal {E}$ satisfies the bound (95). Writing $a_0 a_0^{\dagger } = a_0^{\dagger } a_0 + 1 = N - \mathcal {N}_+ + 1$ , we could then replace

$$\begin{align*}\sum_{|p| \leq N^{\alpha}} \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} \, a_p^{\dagger} \frac{a_0 a_0^{\dagger}}{N} a_p = \sum_{|p| \leq N^{\alpha}} \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} \, a_p^{\dagger} a_p + \delta \end{align*}$$

with

$$\begin{align*}\pm \delta \leq \frac{1}{N} \sum_{|p| \leq N^{\alpha}} \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} a_p^{\dagger} (\mathcal{N}_+ + 1) a_p \lesssim N^{-1+2\alpha} (\mathcal{N}_++1)^2. \end{align*}$$

Furthermore, since

$$ \begin{align*} \Big| \, p^2 + 8 \pi \mathfrak{a}_N - \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} - \frac{(8\pi \mathfrak{a}_N)^2}{2p^2} \Big| \lesssim (1+p^2)^{-2}, \end{align*} $$

we could write

$$ \begin{align*} \sum_{|p|\leq N^{\alpha}} &\left[ p^2 + 8\pi \mathfrak{a}_N - \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} - \frac{(8\pi \mathfrak{a}_N)^2}{2p^2}\right] \\ &\qquad\qquad\qquad\qquad = \sum_{p} \left[p^2 + 8\pi \mathfrak{a}_N - \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} - \frac{(8\pi \mathfrak{a}_N)^2}{2p^2}\right] + {\mathcal O}(N^{-\alpha}). \end{align*} $$

Similarly, from $\big | p^2 - \sqrt {|p|^4 + 16 \pi \mathfrak {a}_N p^2} \big | \lesssim 1$ , we could bound

$$\begin{align*}\pm \sum_{|p|> N^{\alpha}} \big[ p^2 - \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} \big] a_p^{\dagger} a_p \lesssim N^{-2\alpha} H_1. \end{align*}$$

Inserting all these estimates in equation (97), we would end up with (94) and (95).

It remains to show the bound (96). To this end, we observe that

(98) $$ \begin{align} f^{\prime}_{\psi} (s) = \frac{d}{ds} \big\langle \psi, e^{s {\mathcal{B}}_4} E(s) e^{-s{\mathcal{B}}_4} \psi \big\rangle = \big\langle \psi, e^{s {\mathcal{B}}_4} \big\{ \big[ {\mathcal{B}}_4, E(s) \big] + \frac{\partial E(s)}{\partial s} \big\} e^{-s{\mathcal{B}}_4} \psi \big\rangle. \end{align} $$

We have, denoting $\varepsilon _p = \sqrt {|p|^4 + 16 \pi \mathfrak {a}_N p^2}$ ,

$$\begin{align*}\begin{aligned} \big[ {\mathcal{B}}_4, E(s) \big] = \frac{1}{2} \sum_{|p| \leq N^{\alpha}} \varepsilon_p \sum_{|q| \leq N^{\alpha}} \tau_q \big\{ &\big[ b_q^{\dagger} b_{-q}^{\dagger}, (\gamma_p^s b_p^{\dagger} + \nu_p^s b_{-p}) \big] (\gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger}) \\ &+ (\gamma_p^s b_p^{\dagger} + \nu_p^s b_{-p}) \big[ b_q^{\dagger} b_{-q}^{\dagger}, (\gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger}) \big] \big\} + \text{h.c.} \end{aligned} \end{align*}$$

A long but straightforward computation, based on the commutation relations (89), leads to

(99) $$ \begin{align} \begin{aligned} \big[ &{\mathcal{B}}_4, E(s) \big] \\ = &\, -\frac{1}{2}\kern-1.2pt \sum_{|p| \leq N^{\alpha}} \varepsilon_p \tau_p \nu_p^s b_p^{\dagger} \big(1-\frac{\mathcal{N}_+}{N}\big) (\gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger} ) - \frac{1}{2} \kern-1.2pt\sum_{|p| \leq N^{\alpha}} \kern-1.2pt\varepsilon_p \tau_p \nu_p^s \big(1-\frac{\mathcal{N}_+}{N} \big) b_p^{\dagger} (\gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger} ) \\ &-\frac{1}{2} \sum_{|p| \leq N^{\alpha}} \varepsilon_p \tau_p \gamma_p^s (\gamma_p^s b^{\dagger}_p + \nu_p^s b_{-p} ) b_{-p}^{\dagger} \big(1-\frac{\mathcal{N}_+}{N}\big) -\frac{1}{2} \sum_{|p| \leq N^{\alpha}} \varepsilon_p \tau_p \gamma_p^s (\gamma_p^s b^{\dagger}_p + \nu_p^s b_{-p} )\big(1-\frac{\mathcal{N}_+}{N}\big) b_{-p}^{\dagger} \\ &+ \frac{1}{2N} \sum_{|p|, |q| \leq N^{\alpha}} \varepsilon_p \tau_q \nu_p^s \big[ b_q^{\dagger} a_{-q}^{\dagger} a_{-p} + a_q^{\dagger} a_{-p} b_{-q}^{\dagger} \big] (\gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger} ) \\ &+ \frac{1}{2N} \sum_{|p|, |q| \leq N^{\alpha}} \varepsilon_p \tau_q \gamma_p^s (\gamma_p^s b^{\dagger}_p + \nu_p^s b_{-p}) \big[ b_q^{\dagger} a_{-q}^{\dagger} a_{p} + a_q^{\dagger} a_{p} b_{-q}^{\dagger} \big] + \text{h.c.} \end{aligned} \end{align} $$

To compute the explicit time derivative of the observable $E(s)$ , on the other hand, we notice that $ d\gamma _p^s /ds = \tau _p \nu _p^s$ and $d\nu _p^s / ds = \tau _p \gamma _p^s$ . Thus, we obtain

$$\begin{align*}\begin{aligned} \frac{\partial E(s)}{\partial s} = &\!\sum_{|p| \leq N^{\alpha}} \!\!\varepsilon_p \tau_p \big(\nu_p^s b_p^{\dagger} + \gamma_p^s b_{-p} \big) \big( \gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger} \big) +\! \sum_{|p| \leq N^{\alpha}} \!\!\varepsilon_p \tau_p (\gamma_p^s b_p^{\dagger} + \nu_p^s b_{-p}) (\nu_p^s b_p + \gamma_p^s b_{-p}^{\dagger}). \end{aligned} \end{align*}$$

Combining the last equation with equation (99), we observe that (as expected) all large contributions cancel. We find

$$\begin{align*}\begin{aligned} \big[ {\mathcal{B}}_4, E(s) \big] + \frac{\partial E(s)}{\partial s} = \; &\frac{1}{N} \sum_{|p| \leq N^{\alpha}} \varepsilon_p \tau_p \nu_p^s \, \mathcal{N}_+ b_p^{\dagger} (\gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger}) \\ &+\frac{1}{N} \sum_{|p| \leq N^{\alpha}} \varepsilon_p \tau_p \gamma_p^s \, (\gamma_p^s b^{\dagger}_p + \nu_p^s b_{-p}) \mathcal{N}_+ b_{-p}^{\dagger} \\ &+\frac{1}{N} \sum_{|p|, |q| \leq N^{\alpha}} \varepsilon_p \tau_q \nu_p^s \, b_q^{\dagger} a_{-q}^{\dagger} a_{-p} (\gamma_p^s b_p + \nu_p^s b_{-p}^{\dagger}) \\ &+ \frac{1}{N} \sum_{|p|, |q| \leq N^{\alpha}} \varepsilon_p \tau_q \gamma_p^s \, (\gamma_p^s b_p^{\dagger} + \nu_p^s b_{-p}) b_q^{\dagger} a_{-q}^{\dagger} a_p + \text{h.c.} \end{aligned} \end{align*}$$

Using the bounds in Lemma 17, with the estimate $\varepsilon _p \lesssim p^2$ and the restrictions $|p|, |q| \leq N^{\alpha }$ , we arrive at

$$\begin{align*}\pm \left\{ \big[ {\mathcal{B}}_4, E(s) \big] + \frac{\partial E(s)}{\partial s} \right\} \lesssim N^{2\alpha-1} (\mathcal{N}_+ + 1)^2. \end{align*}$$

From equation (98), this implies that

$$\begin{align*}|f^{\prime}_{\psi} (s)| \lesssim N^{2\alpha-1} \langle \psi, e^{-s{\mathcal{B}}_4} (\mathcal{N}_+ + 1)^2 e^{s{\mathcal{B}}_4} \psi \rangle. \end{align*}$$

Applying Lemma 18, we obtain the desired bound (96).

6 Optimal BEC and proof of Theorem 1

Let us denote

(100) $$ \begin{align} \widetilde H_N = H_N - 4 \pi \mathfrak{a}_N (N-1) + \frac{1}{2}\sum_{p} \left[ \sqrt{|p|^4 + 16 \pi \mathfrak{a}_N p^2} - p^2 - 8 \pi \mathfrak{a}_N + \frac{(8\pi \mathfrak{a}_N)^2}{2p^2}\right], \end{align} $$

and

(101) $$ \begin{align} E_{\infty} = \sum_{p} \sqrt{|p|^4 + p^2 16 \pi \mathfrak{a}_N } \, a_p^{\dagger} a_p. \end{align} $$

Moreover, let $\mathcal {U} = e^{{\mathcal {B}}_2} e^{{\mathcal {B}}_3} e^{{\mathcal {B}}_4}$ . Observe that $\mathcal {U}$ is a unitary operator. From Proposition 19, we have

(102) $$ \begin{align} {\mathcal U}^{\dagger} \widetilde H_N {\mathcal U} = E_{\infty} + e^{-{\mathcal{B}}_4}Q_4e^{{\mathcal{B}}_4} + \mathcal E_{{\mathcal{B}}_4} , \end{align} $$

where $\mathcal {E}_{{\mathcal {B}}_4}$ satisfies the bound (95). To prove that the error term $\mathcal {E}_{{\mathcal {B}}_4}$ is small, we show first that low-energy states exhibit complete Bose-Einstein condensation.

Proposition 20 (Optimal BEC).

On $\{\mathcal {N} = N\}$ , we have

(103) $$ \begin{align} H_N \geq 4\pi \mathfrak{a}_N N + C^{-1} \mathcal{N}_+ - C, \end{align} $$

for some constant $C>0$ independent of N.

Proof. To take care of the terms on the second line of the bound (95), we use localisation in the number of particles, a tool developed in [Reference Lewin, Nam, Serfaty and Solovej17] and, in the present setting, in [Reference Boccato, Brennecke, Cenatiempo and Schlein7]. Here, we make use of the results of [Reference Lieb and Seiringer18, Reference Lieb and Seiringer19, Reference Nam, Rougerie and Seiringer24], which imply that, if $\psi _N \in L^2_s (\Lambda ^N)$ is a normalised sequence of approximate ground states of the Hamilton operator $H_N$ satisfying

$$\begin{align*}\Big| \frac{1}{N} \langle \psi_N, H_N \psi_N \rangle - 4\pi \mathfrak{a}_N \Big| \to 0 \end{align*}$$

as $N \to \infty $ , then $\psi _N$ exhibit condensation, in the sense that

(104) $$ \begin{align} \lim_{N \to \infty} \frac{1}{N} \langle \psi_N , \mathcal{N}_+ \psi_N \rangle = 0. \end{align} $$

Now let $f,g : \mathbb {R} \to [0,1]$ be smooth functions such that $f(s)^2 + g(s)^2 = 1$ for all $s \in \mathbb {R}$ , and $f(s) = 1$ for $s \leq 1/2$ , $f(s) =0$ for $s\geq 1$ . For $M_0 \geq 1$ , we define $f_{M_0}(\mathcal {N}_+) = f(\mathcal {N}_+ / M_0)$ and $g_{M_0}(\mathcal {N}_+) = g(\mathcal {N}_+/M_0)$ . Then we have

(105) $$ \begin{align} H_N &= f_{M_0} H_N f_{M_0} + g_{M_0} H_N g_{M_0} + \mathcal E_{M_0}, \end{align} $$

with

(106) $$ \begin{align} \mathcal E_{M_0} = \frac{1}{2} \left( [f_{M_0},[f_{M_0},H_N]] + [g_{M_0},[g_{M_0},H_N]]\right). \end{align} $$

In view of equation (9), we can write (with $h =f,g$ )

$$ \begin{align*} [h_{M_0} , [h_{M_0},H_N]] &= \big[ h (\mathcal{N}_+ / M_0) - h ((\mathcal{N}_+ -2)/M_0)\big]^2 \sum_{p \neq 0} \hat V_N(p) a_p^{\dagger} a_{-p}^{\dagger} a_0 a_0 + \text{h.c.}, \\ &\quad + \big[ h (\mathcal{N}_+ / M_0) - h ((\mathcal{N}_+ -1)/M_0)\big]^2 \sum_{q,r,q+r\neq 0} \hat V_N(r) a_{q+r}^{\dagger} a_{-r}^{\dagger} a_{q}a_0 + \text{h.c.} \end{align*} $$

This easily implies that

(107)

We choose $M_0 = \varepsilon N$ , for some $\varepsilon>0$ independent of N, to be fixed later. We introduce the notation $\mathcal {N}_{+}^{\, {\mathcal U}} = {\mathcal U}^{\dagger } \mathcal {N}_+ {\mathcal U}$ . We use equation (102) with the bound (95) for the error term $\mathcal {E}_4$ ; we pick $\alpha = - \log \ell / \log N$ so that $N^{\alpha } = 1/\ell $ , for some $\ell> 0$ independent of N, to be specified below. For N large enough, we obtain from Proposition 19

(108) $$ \begin{align} f_{M_0}(\mathcal{N}_+) H_N &f_{M_0}(\mathcal{N}_+) = {\mathcal U} f_{M_0}(\mathcal{N}_{+}^{\, {\mathcal U}}) {\mathcal U}^{\dagger} H_N {\mathcal U} \, f_{M_0}(\mathcal{N}_{+}^{\, {\mathcal U}}) {\mathcal U}^{\dagger} \nonumber \\ &\geq {\mathcal U} f_{M_0}(\mathcal{N}_{+}^{\, {\mathcal U}}) \left( 4\pi \mathfrak{a}_N N + H_1 -C + e^{-{\mathcal{B}}_4}Q_4e^{{\mathcal{B}}_4} + \mathcal E_4 \right)f_{M_0}(\mathcal{N}_{+}^{\, {\mathcal U}}) {\mathcal U}^{\dagger}\nonumber \\ &\geq {\mathcal U} f_{M_0}(\mathcal{N}_{+}^{\, {\mathcal U}}) \big( 4\pi \mathfrak{a}_N N + (1 - C \ell^{1/2} - C\ell^{-5/2} N^{-1}- C \ell^{-7/2}\varepsilon) (H_1+1) \nonumber \\ &\qquad \qquad\qquad\qquad\qquad\qquad\qquad + (1-C \ell^{1/2}) e^{-{\mathcal{B}}_4}Q_4e^{{\mathcal{B}}_4}\big)f_{M_0}(\mathcal{N}_{+}^{\, {\mathcal U}}) {\mathcal U}^{\dagger} \nonumber \\ &\geq f_{M_0}(\mathcal{N}_{+})^2 (4\pi \mathfrak{a}_N N - C + C^{-1} \mathcal{N}_+), \end{align} $$

choosing first $\ell> 0$ small enough and then $\varepsilon> 0$ sufficiently small. Here, we used $(\mathcal {N}_++1)^j \lesssim (\mathcal {N}_+^{\, {\mathcal U}}+1)^j$ (as follows from Lemma 3, Lemma 10 and Lemma 18), to estimate the error terms on the second line of (95). Moreover, we used the bounds $\mathcal {N}_+ , \mathcal {N}^{{\mathcal U}}_+ \lesssim H_1$ .

On the other hand, following an argument from [Reference Boccato, Brennecke, Cenatiempo and Schlein7, Prop. 6.1], we find

(109) $$ \begin{align} g_{M_0}(\mathcal{N}_+) \left(H_N -4\pi \mathfrak{a}_N N\right) &g_{M_0}(\mathcal{N}_+) \geq C^{-1} \mathcal{N}_+ g_{M_0}(\mathcal{N}_+)^2. \end{align} $$

Indeed, otherwise, we could find a normalised sequence $\Psi _N$ , supported on $\{\mathcal {N}_+> \varepsilon N\}$ , satisfying

$$\begin{align*}\Big| \frac{1}{N} \langle \psi_N , H_N \psi_N \rangle_{\Psi_N} - 4\pi \mathfrak{a}_N \Big| \to 0 \end{align*}$$

as $N \to \infty $ , in contradiction with (104).

Finally, we deal with the error term $\mathcal E_{M_0=\varepsilon N}$ . For $\psi _N$ with $\langle \psi _N, H_N \psi _N \rangle \leq C N$ , we immediately find, from the bound (107), that

$$\begin{align*}\langle \psi_N, \mathcal E_{M_0 = \varepsilon N} \psi_N \rangle \lesssim \varepsilon^{-2} N^{-1}. \end{align*}$$

Since the bound (103) holds trivially on states with $\langle \psi _N , H_N \psi _N \rangle \geq CN$ , this, together with the estimates (108) and (109), concludes the proof of Proposition 20.

With Proposition 19 and Proposition 20, we are now ready to show Theorem 1, determining the low-energy spectrum of the operator Hamilton operator $H_N$ .

Proof of Theorem 1.

We continue to use the notation $\widetilde {H}_N$ and $E_{\infty }$ introduced in equations (100) and (101). Moreover, we denote by $\lambda _1 (\widetilde {H}_N) \leq \lambda _2 (\widetilde {H}_N) \leq \dots $ and $\lambda _1 (E_{\infty }) \leq \lambda _2 (E_{\infty }) \leq \dots $ the ordered eigenvalues of $\widetilde {H}_N$ and, respectively, $E_{\infty }$ . We now choose $L \in \mathbb {N}$ , with $\lambda _L (\widetilde H_N) \leq \Theta $ for some $1 \leq \Theta \leq N^{1/17}$ . Then we claim that

(110) $$ \begin{align} \lambda_L (\widetilde H_N) = \lambda_L(E_{\infty}) + {\mathcal O} (\Theta N^{-1/17}). \end{align} $$

Since $\lambda _0 (E_{\infty }) = 0$ , the estimate (110) shows that the ground state energy $E_N$ of $H_N$ satisfies the estimate (5). It is then easy to check, using (110), that the excitations of $H_N - E_N$ satisfy the claim in (4).

To prove equation (110), we show first a lower bound and then a matching upper bound. We again use Proposition 19, but this time we choose the exponents $\alpha = 2/17$ .

Lower bound on $\lambda _L (\widetilde {H}_N)$ . We again use the localisation identity given by equation (105), but this time we take $M_0 = N^{1/2 + 1/34}$ . Let Y denote the subspace generated by the first L eigenfunctions of $\widetilde H_N$ , and let us denote $Z = {\mathcal U}^{\dagger } Y$ , which is of dimension L. From the decomposition (105), we have

(111) $$ \begin{align} \lambda_{L}(\widetilde H_N) \geq P_Y \left(f_{M_0}(\mathcal{N}_+) \widetilde H_N f_{M_0}(\mathcal{N}_+) + g_{M_0}(\mathcal{N}_+)\widetilde H_N g_{M_0}(\mathcal{N}_+) + \mathcal E_{M_0}\right) P_Y. \end{align} $$

From Proposition 20, we have

$$ \begin{align*} g_{M_0}(\mathcal{N}_+) \widetilde H_N g_{M_0}(\mathcal{N}_+) \geq C g^2_{M_0}(\mathcal{N}_+) (C^{-1} M_0 - C) \geq 0 \end{align*} $$

for N large enough (recall the choice $M_0 = N^{1/2+1/34}$ ). Here, we used that $g_{M_0}$ is supported on $\mathcal {N}_+> M_0/3$ . Moreover, with the bound (107), we find

$$ \begin{align*} P_Y \mathcal E_{M_0} P_Y &\geq - C M_0^{-2} P_Y (Q_4 + N) P_Y \geq - C M_0^{-2} P_Y (H_N + N) P_Y \geq - C M_0^2 N \geq - C N^{-1/17} \end{align*} $$

because (from the upper bound), we know that $H_N \leq C N$ on Y. From the estimate (111), we obtain

$$ \begin{align*} \lambda_{L}(\widetilde H_N) \geq P_Y f_{M_0}(\mathcal{N}_+) \widetilde H_N f_{M_0}(\mathcal{N}_+) P_Y - CN^{-1/17}. \end{align*} $$

We now use Proposition 19 to estimate

Using that $\mathcal {N}_+ \lesssim H_1 \leq E_{\infty }$ and the choices $M_0 = N^{1/2 + 1/34}$ , $\alpha = 2/17$ , we find

$$ \begin{align*} P_Y f_{M_0}(\mathcal{N}_+) \,&\widetilde H_N \,f_{M_0}(\mathcal{N}_+) P_Y \\ &\geq \left(1- C N^{-1/17}\right) {\mathcal U} P_Z f_{M_0}(\mathcal{N}_+^{{\mathcal U}}) e^{{\mathcal{B}}_4} \, E_{\infty} \, e^{-{\mathcal{B}}_4} f_{M_0}(\mathcal{N}_+^{{\mathcal U}})P_Z {\mathcal U}^{\dagger}. \end{align*} $$

Now it turns out that for N large enough, $\dim f_{M_0}(\mathcal {N}_+^{{\mathcal U}})P_Z = L$ because

$$ \begin{align*} \max_{\xi \in P_Z} \frac{\bigg\|\sqrt{1-f_{M_0}(\mathcal{N}_+^{{\mathcal U}})^2}\xi\bigg\|^2}{\|\xi\|^2} \leq C M_0^{-1} \max_{\xi \in P_Y} \frac{\|\mathcal{N}_+^{1/2}\xi\|^2}{\|\xi\|^2} \leq C \lambda_{L}(\widetilde H_N) M_0^{-1} \underset{N \to \infty}{\longrightarrow} 0; \end{align*} $$

see, for instance, [Reference Lewin, Nam, Serfaty and Solovej17, Prop. 6.1 ii)]. Thus

$$ \begin{align*} \lambda_{L}(\widetilde H_N) &\geq \max_{\xi \in Y} \frac{\langle e^{-{\mathcal{B}}_4} f_{M_0}(\mathcal{N}_+^{{\mathcal U}}) P_Z \mathcal{U}^{\dagger} \xi , E_{\infty} e^{-{\mathcal{B}}_4} f_{M_0} (\mathcal{N}_+^{{\mathcal U}}) P_Z {\mathcal U}^{\dagger}\xi \rangle}{\|\xi\|^2} - C \Theta N^{-1/17} \\ &\geq \max_{\xi \in e^{-{\mathcal{B}}_4} f_{M_0}(\mathcal{N}_+^{{\mathcal U}}) Z} \frac{\langle \xi, E_{\infty} \xi \rangle}{\|\xi\|^2}(1- C \lambda_{L}(\widetilde H_N) M_0^{-1}) - C \Theta N^{-1/17} \\ &\geq \min_{\dim X = L} \max_{\xi \in X} \frac{\langle \xi , E_{\infty} \xi \rangle}{\|\xi\|^2}(1- C \lambda_{L}(\widetilde H_N) M_0^{-1})- C\Theta N^{-1/17} \\ &\geq \lambda_L (E_{\infty}) (1 -C\Theta N^{-1/2}) -C \Theta N^{-1/17}, \end{align*} $$

which implies $\lambda _L(\widetilde H_N) \geq \lambda _{L}(E_{\infty }) + {\mathcal O}(\Theta N^{-1/17})$ .

Upper bound on $\lambda _L (\widetilde {H}_N)$ . Let Z denote the subspace generated by the first L eigenfunctions of $E_{\infty }$ and $P_Z$ be the orthogonal projection onto Z. The normalised eigenfunctions of $E_{\infty }$ have the form

(112) $$ \begin{align} \xi = \prod_{j=1}^k \frac{a^{\dagger}(p_j)^{n_j}}{\sqrt{n_j!}} \Omega \end{align} $$

for some $k\geq 1$ , $p_j \in \Lambda _+^*$ , $n_j \geq 1$ and where $\Omega $ is the vacuum. Note that

$$ \begin{align*} P_Z \mathcal{N}_+ P_Z \lesssim P_Z H_1 P_Z \leq P_Z E_{\infty} P_Z \leq \lambda_{L}(E_{\infty}) \leq C \Theta \leq C N^{1/17}, \end{align*} $$

where we used the lower bound that we proved above. Since $[P_Z, \mathcal {N}_+]=0$ , this bound can also be applied to powers of $\mathcal {N}_+$ . Note also that $e^{{\mathcal {B}}_4} E_{\infty } e^{-{\mathcal {B}}_4}$ almost commutes with $\mathcal {N}_+$ in the sense that

Hence, we have

Together with Proposition 19, we find

(113)

Again because

we have

for N large enough. With Lemma 18, we obtain

To estimate $P_Z Q_4 P_Z$ , we use an argument from [Reference Boccato, Brennecke, Cenatiempo and Schlein6, Lemma 6.1]: from $P_Z E_{\infty } P_Z \leq \Theta $ , we must have $a_p \xi = 0$ for all $|p|> \Theta ^{1/2}$ and $\xi \in Z$ . This implies that

$$ \begin{align*} \langle \xi, Q_4 \xi \rangle &\leq \sum_{p,q,r} \hat V_N(r) \chi_{|r|\leq \Theta^{1/2}} \|a_{p+r} a_q \xi \| \|a_{p} a_{q+r} \xi \| \\ & \leq C \Theta^{3/2} N^{-1} \|(\mathcal{N}_++1)\xi\|^2 \leq C \Theta^{7/2} N^{-1} \|\xi\|^2 \leq C N^{-1/17} \|\xi\|^2, \end{align*} $$

for all $\xi \in Z$ . Applying the min-max principle, we conclude from the estimate (113) that

Acknowledgements

We thank P.T. Nam for insightful discussions.

Conflict of Interest

None.

Financial Support

B.S. gratefully acknowledges support from the European Research Council through the ERC Advanced Grant CLaQS. Additionally, B.S. acknowledges partial support from the NCCR SwissMAP and the Swiss National Science Foundation through the grant ‘Dynamical and energetic properties of Bose-Einstein condensates’.

Footnotes

1 The estimate for $Q_4$ will only be used in the next section to show upper bounds on the eigenvalues of $H_N$ ; for the lower bounds, we will only use the fact that $Q_4 \geq 0$ .

References

Adhikari, A., Brennecke, C., Schlein, B.. ‘Bose-Einstein condensation beyond the Gross-Pitaevskii regime’. Ann. Henri Poincaré 22, (2021) 11631233.CrossRefGoogle ScholarPubMed
Basti, G., Cenatiempo, S., Olgiati, A., Pasqualetti, G., Schlein, B.. ‘Ground state energy of a Bose gas in the Gross-Pitaevskii regime’. J. Math. Phys. 63 (2022), 041101. https://doi.org/10.1063/5.0087116.CrossRefGoogle Scholar
Basti, G., Cenatiempo, S., Schlein, B.. ‘A new second-order upper bound for the ground state energy of dilute Bose gases’. Forum of Mathematics, Sigma 9 (2021), e74.CrossRefGoogle Scholar
Boccato, C., Brennecke, C., Cenatiempo, S., Schlein, B.. ‘Complete Bose-Einstein condensation in the Gross-Pitaevskii regime’. Comm. Math. Phys. 359(3) (2018), 9751026.CrossRefGoogle Scholar
Boccato, C., Brennecke, C., Cenatiempo, S., Schlein, B.. ‘The excitation spectrum of Bose gases interacting through singular potentials’. J. Eur. Math. Soc. 22(7) (2020), 23312403.CrossRefGoogle Scholar
Boccato, C., Brennecke, C., Cenatiempo, S., Schlein, B.. ‘Bogoliubov Theory in the Gross-Pitaevskii limit’. Acta Math. 222(2) (2019), 219335.CrossRefGoogle Scholar
Boccato, C., Brennecke, C., Cenatiempo, S., Schlein, B.. ‘Optimal rate for Bose-Einstein condensation in the Gross-Pitaevskii regime’. Comm. Math. Phys. 376 (2020), 13111395.CrossRefGoogle Scholar
Bogoliubov, N. N.. ‘On the theory of superfluidity’. Izv. Akad. Nauk. USSR 11 (1947), 77. Engl. Transl. J. Phys. (USSR) 11 (1947), 23.Google Scholar
Brennecke, C., Caporaletti, M., Schlein, B.. ‘Excitation spectrum for Bose gases beyond the Gross-Pitaevskii regime’. Preprint arxiv:2104.13003; https://doi.org/10.1142/S0129055X22500271.CrossRefGoogle Scholar
Brennecke, C., Schlein, B., Schraven, S.. ‘Bose-Einstein condensation with optimal rate for trapped bosons in the Gross-Pitaevskii regime’. Math. Phys. Anal. Geom. 25 (2022), 12. https://doi.org/10.1007/s11040-022-09424-7.CrossRefGoogle Scholar
Brennecke, C., Schlein, B., Schraven, S.. ‘Bogoliubov theory for trapped bosons in the Gross-Pitaevskii regime’. Ann. Henri Poincaré (2022) https://doi.org/10.1007/s00023-021-01151-z.CrossRefGoogle Scholar
Brennecke, C., Schlein, B.. ‘Gross-Pitaevskii dynamics for Bose-Einstein condensates’. Anal. PDE 12(6) (2019), 15131596.CrossRefGoogle Scholar
Fournais, S., Solovej, J.P.. ‘The energy of dilute Bose gases’. Annals of Mathematics, 192 (2020).CrossRefGoogle Scholar
Fournais, S., Solovej, J.P.. ‘The energy of dilute Bose gases II: The general case’. Preprint arXiv:2108.12022.Google Scholar
Hainzl, C.. ‘Another proof of BEC in the GP-limit’. J. Math. Phys. 62(5) (2021), 459485.CrossRefGoogle Scholar
Hainzl, C., Seiringer, R.. ‘The BCS critical temperature for potentials with negative scattering length’. Lett. Math. Phys. 84(2-3) (2008), 99107.CrossRefGoogle Scholar
Lewin, M., Nam, P. T., Serfaty, S., Solovej, J.P.. ‘Bogoliubov spectrum of interacting Bose gases’. Comm. Pure Appl. Math. 68(3) (2014), 413471.CrossRefGoogle Scholar
Lieb, E. H. and Seiringer, R.. ‘Proof of Bose-Einstein condensation for dilute trapped gases’. Phys. Rev. Lett. 88 (2002), 170409.CrossRefGoogle ScholarPubMed
Lieb, E. H. and Seiringer, R.. ‘Derivation of the Gross-Pitaevskii equation for rotating Bose gases’. Comm. Math. Phys. 264(2) (2006), 505537.CrossRefGoogle Scholar
Lieb, E. H., Seiringer, R., and Yngvason, J.. ‘Bosons in a trap: A rigorous derivation of the Gross-Pitaevskii energy functional’. Phys. Rev. A 61 (2000), 043602.CrossRefGoogle Scholar
Lieb, E. H. and Seiringer, R.. ‘Derivation of the Gross-Pitaevskii equation for rotating Bose gases’. Comm. Math. Phys. 264 (2006), 505537.CrossRefGoogle Scholar
Lieb, E. H., Yngvason, J.. ‘Ground state energy of the low density Bose gas’. Phys. Rev. Lett. 80 (1998), 25042507.CrossRefGoogle Scholar
Nam, P. T., Napiórkowski, M., Ricaud, J., Triay, A.. ‘Optimal rate of condensation for trapped bosons in the Gross–Pitaevskii regime’. Anal. PDE (in press)Google Scholar
Nam, P. T., Rougerie, N., Seiringer, R.. ‘Ground states of large bosonic systems: The Gross-Pitaevskii limit revisited’. Anal. PDE 9(2) (2016), 459485.CrossRefGoogle Scholar
Nam, P. T., Triay, A.. ‘Bogoliubov excitation spectrum of trapped Bose gases in the Gross-Pitaevskii regime’. Preprint arXiv:2106.11949.Google Scholar
Pachpatte, B. G., Inequalities for Differential and Integral Equations. Elsevier (1997).Google Scholar
Yau, H.-T., Yin, J.. ‘The second order upper bound for the ground state energy of a Bose gas’. J. Stat. Phys. 136(3) (2009), 453503.CrossRefGoogle Scholar