1 Introduction
1.1 The Model
The large polaron provides an idealized description for the motion of a slow band electron through a polarizable crystal. The analysis of the polaron is a classic problem in solid-state physics that first appeared in 1933 when Landau introduced the idea of self-trapping of an electron in a polarizable environment [Reference Landau30]. Since it provides a simple model for a particle interacting with a nonrelativistic quantum field, the polaron has been of interest also in field theory and mathematical physics. In particular, the strong coupling theory of the polaron and Pekar’s adiabatic approximation have been the source of interesting and challenging mathematical problems.
Following H. Fröhlich [Reference Fröhlich20], the Hamiltonian of the model acts on the Hilbert space
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn1.png?pub-status=live)
with
${\mathcal F}$
the bosonic Fock space over
$L^2({\mathbb R}^3)$
and is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn2.png?pub-status=live)
Here,
$x\in {\mathbb R}^3$
is the coordinate of the electron,
${\mathbb N}$
denotes the number operator on Fock space and the field operator
$\phi (h_x) =a^{\dagger }(h_x) + a(h_x)$
with coupling function
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn3.png?pub-status=live)
accounts for the interaction between the electron and the quantum field. The creation and annihilation operators satisfy the usual canonical commutation relations
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn4.png?pub-status=live)
Since we set
$\hbar = 1$
and the mass of the electron equal to
$1/2$
, the only free parameter is the coupling constant
$\alpha>0$
.
By rescaling all lengths by a factor
$1/\alpha $
, one can show that
$\alpha ^2 H_{\alpha }$
is unitarily equivalent to the Hamiltonian
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn5.png?pub-status=live)
which is more common in the polaron literature and also explains why
$\alpha \to \infty $
is called the strong coupling limit.
The Fröhlich Hamiltonian defines a translation invariant model, that is, it commutes with the total momentum operator,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn6.png?pub-status=live)
where
$P_f=\mathrm{d} \Gamma (-i\nabla )$
denotes the momentum operator of the phonons. This allows the definition of the energy-momentum relation
$E_{\alpha }(P)$
as the lowest possible energy of
$H_{\alpha }$
when restricted to states with total momentum
$P\in {\mathbb R}^3$
. To this end, it is convenient to switch to the Lee–Low–Pines representation
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn7.png?pub-status=live)
where
$H_{\alpha }(P)$
acts on the Fock space only [Reference Lee, Low and Pines32]. The Fröhlich Hamiltonian
$H_{\alpha }$
is unitarily equivalent to the fiber decomposition
$\int ^{\oplus }_{{\mathbb R}^3} H_{\alpha }(P) \mathrm{d} P$
, which follows easily from transforming
$H_{\alpha }$
with
$e^{iP_fx}$
and diagonalizing the obtained operator in the electron coordinate. The energy-momentum relation is then defined as the ground state energy of the fiber Hamiltonian,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn8.png?pub-status=live)
which by construction satisfies
$E_{\alpha }( R P ) = E_{\alpha } (P)$
for all rotations
$R\in \mathrm{SO}(3)$
. It is known [Reference Gross26] that
$E_{\alpha }(0) \le E_{\alpha }(P)$
and hence
$E_{\alpha }(0)= \inf \sigma (H_{\alpha }) $
(in fact,
$E_{\alpha }(0) < E_{\alpha }(P)$
for all
$P\neq 0$
[Reference Lampart, Mitrouskas and Myśliwy29, Reference Polzer52]). Further properties, such as the domain of analyticity, existence of ground states and the value of the bottom of the essential spectrum, were analyzed in [Reference Miyao44, Reference Fröhlich21, Reference Spohn58, Reference Møller45, Reference Gerlach and Löwen23, Reference Dahlbæk10].
The aim of this work is to analyze the quantitative behavior of the energy-momentum relation for large coupling
$\alpha \to \infty $
. Our main result provides an upper bound for
$E_{\alpha } (P)$
. The upper bound consists of a momentum independent part coinciding with the optimal upper bound for the ground state energy of the strongly coupled polaron at rest and a momentum-dependent part. In more detail, the momentum-independent part is given by the classical Pekar energy and the corresponding quantum fluctuations that are described by the energy of a system of harmonic oscillators with frequencies determined by the Hessian of the corresponding classical field functional. This part agrees with the expected asymptotic form of
$E_{\alpha }(0)$
; see equation (1.25). The momentum-dependent part, on the other hand, describes the energy of a free particle with mass
$M(\alpha ) = \frac {2 \alpha ^4}{3} \int |\nabla \varphi |^2$
, where
$\varphi $
denotes the self-consistent polarization field, which coincides with the classical polaron mass introduced by Landau and Pekar [Reference Landau and Pekar31]; see equation (1.24). As will be explained in Section 1.3, our result confirms the heuristic picture of the polaron (the electron and the accompanying classical field) as a free quasi-particle with largely enhanced mass. To our best knowledge, the upper bound we present in this work is the first rigorous result about the connection between the energy-momentum relation
$E_{\alpha }(P)$
and the classical polaron mass
$M(\alpha )$
.
Starting from the works in the 1930s and 1940s [Reference Landau30, Reference Landau and Pekar31, Reference Fröhlich19], there has been a large number of publications in the physics literature that studied the ground state energy
$E_{\alpha }(0)$
and the effective mass, that is, the inverse curvature of
$E_{\alpha }(P)$
at
$P=0$
. For a comprehensive summary of the earlier results, we refer to [Reference Mitra, Chatterjee and Mukhopadhyay41]. More recent developments are reviewed in [Reference Alexandrov and Devreese1].
Mathematically rigorous results for the leading-order asymptotics of
$E_{\alpha }(0)$
, for
$\alpha $
large, were obtained by Lieb and Yamazaki [Reference Lieb and Yamazaki40] (with nonmatching upper and lower bounds) and by Donsker and Varadhan [Reference Donsker and Varadhan11] as well as Lieb and Thomas [Reference Lieb and Thomas39]. The effective mass has been studied in [Reference Spohn57, Reference Dybalski and Spohn12, Reference Feliciangeli, Rademacher and Seiringer14, Reference Lieb and Seiringer38, Reference Lieb and Seiringer37, Reference Betz and Polzer4]. Other works have considered confined polarons or polaron models with more regular interaction [Reference Frank and Seiringer18, Reference Feliciangeli and Seiringer15, Reference Myśliwy and Seiringer48]. For completeness, let us also mention recent progress in the understanding of the polaron path measure [Reference Mukherjee and Varadhan47, Reference Betz and Polzer3] as well as the increased interest in the analysis of the Schrödinger time evolution of strongly coupled polarons [Reference Griesemer25, Reference Leopold, Mitrouskas, Rademacher, Schlein and Seiringer34, Reference Leopold, Rademacher, Schlein and Seiringer35, Reference Mitrouskas42, Reference Feliciangeli, Rademacher and Seiringer13, Reference Frank and Gang16, Reference Frank and Schlein17].
1.2 Pekar functionals
The semiclassical theory of the polaron has been introduced by Pekar [Reference Pekar51]. It arises naturally in the context of strong coupling, based on the expectation that the electron and the phonons are adiabatically decoupled, similarly as the electrons are adiabatically decoupled from the heavy nuclei in the well-known Born–Oppenheimer theory [Reference Born and Oppenheimer6, Reference Born and Huang5]. With this in mind, one can minimize the Fröhlich Hamiltonian over product states of the form
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn9.png?pub-status=live)
where
$u\in H^1({\mathbb R}^3)$
is a normalized electron wave function,
$\Omega = (1,0,0,\ldots )$
the Fock space vacuum and
$ e^{a^{\dagger } (\alpha v) } \Omega $
the coherent state, up to normalization, that is associated with a classical field
$\alpha v\in L^2({\mathbb R}^3)$
. A simple computation leads to the Pekar energy functional
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn10.png?pub-status=live)
with polarization potential
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn11.png?pub-status=live)
By completing the square, one can further remove the field variable and obtain the energy functional for the electron wave function,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn12.png?pub-status=live)
which is known [Reference Lieb36] to admit a unique rotational invariant minimizer
$\psi> 0$
(the minimizing property is unique only up to translations and multiplication by a constant phase). Alternatively, one can minimize the Pekar energy functional w.r.t. the electron wave function first. This leads to the classical field functional
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn13.png?pub-status=live)
whose unique rotational invariant minimizer is readily shown to be
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn14.png?pub-status=live)
The corresponding classical ground state energy is called the Pekar energy
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn15.png?pub-status=live)
and by the variational principle it provides an upper bound for
$\inf \sigma (H_{\alpha })$
. The validity of Pekar’s ansatz was rigorously verified by Donsker and Varadhan [Reference Donsker and Varadhan11] who proved that
$ \lim _{\alpha \to \infty } \inf \sigma (H_{\alpha }) = e^{\mathrm{Pek}}$
and subsequently by Lieb and Thomas [Reference Lieb and Thomas39] who added a quantitative bound for the error by showing that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn16.png?pub-status=live)
Given the potential
$V^{\varphi }$
for the field
$\varphi $
, one can define the Schrödinger operator
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn17.png?pub-status=live)
with
$\lambda ^{\mathrm{Pek}} = \inf \sigma (-\Delta + V^{\varphi }(x))<0$
and
$\psi $
the corresponding unique ground state. It follows from general arguments for Schrödinger operators that
$h^{\mathrm{Pek}}$
has a finite spectral gap above zero, and thus the reduced resolvent
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn18.png?pub-status=live)
defines a bounded operator (
$P_{\psi }$
denotes the orthogonal projection onto the state
$\psi $
).
The last object to be introduced in this section is the Hessian
$H^{\mathrm{Pek}}$
of the energy functional
${\mathcal F}$
at its minimizer
$\varphi $
, defined by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn19.png?pub-status=live)
In the following lemma, we collect some important properties of
$H^{\mathrm{Pek}}$
.
Lemma 1.1. The linear operator
$H^{\mathrm{Pek}} : L^2({\mathbb R}^3) \to L^2({\mathbb R}^3)$
has integral kernel
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn20.png?pub-status=live)
and satisfies the following properties.
-
(i)
$0\le H^{\mathrm{Pek}} \le 1$
-
(ii)
$\mathrm {Ker} H^{\mathrm{Pek}} = \mathrm {Span}\{\partial _i \varphi \, : \, i = 1,2,3 \}$
-
(iii)
$H^{\mathrm{Pek}} \ge \tau>0$ when restricted to
$(\mathrm {Ker} H^{\mathrm{Pek}})^{\perp }$
-
(iv)
${\mathrm{Tr}}_{L^2}(1- \sqrt {H^{\mathrm{Pek}}}) < \infty $ .
The proof of the lemma, in particular item (ii), is based on the analysis of the Hessian of the energy functional
${\mathcal E}$
[Reference Lenzmann33]. The details are given in Section 4.
1.3 Motivation and goal of this work
In this work, we are interested in the behavior of the energy-momentum relation
$E_{\alpha }(P)$
for large values of the coupling
$\alpha $
. In general,
$E_{\alpha }(P)$
is expected to interpolate between two distinct regimes (see, for instance, [Reference Gerlach and Smondyrev24, Reference Gerlach, Kalina and Smondyrev22, Reference Whitfield and Puff60, Reference Spohn58]): The quasi-particle regime and the radiative regime. The former corresponds to small momenta, and the expectation is that the system behaves effectively like a free particle with energy
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn21.png?pub-status=live)
where the effective mass is determined by the inverse curvature of
$E_{\alpha }(P)$
at
$P=0$
(which is known to be well defined),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn22.png?pub-status=live)
It is easy to verify that
$M^{\mathrm{eff}} (\alpha ) \ge 1 / 2$
(the mass of the electron in our units), and one can further show that the inequality is strict if
$\alpha>0$
so that the emerging quasi-particle is heavier than the bare electron. The heuristic idea is that the electron drags along a cloud of phonons when it moves through the crystal and thus appears to be heavier than it would be without the interaction. The radiative regime, on the other hand, describes a polaron at rest and an unbound/radiative phonon carrying the total momentum P. It is expected to be valid for large momenta and it is characterized by a flat energy-momentum relation that equals or approaches the bottom of the essential spectrum [Reference Møller45] (see also [Reference Lampart, Mitrouskas and Myśliwy29, Lemma 1.1])
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn23.png?pub-status=live)
The two regimes cross at
$|P| = P_{\mathrm c}(\alpha ) : =\sqrt {2M^{\mathrm{eff}}(\alpha )} / \alpha $
which marks a characteristic momentum scale of the polaron. While the quasi-particle picture is expected to be accurate for
$|P|\lesssim P_{\mathrm c}(\alpha )$
, the radiative regime should hold for
$|P|\gtrsim P_{\mathrm c}(\alpha )$
(see also Remark 1.3 below). Between the two regimes there is no concrete prediction for the behavior of
$E_{\alpha }(P)$
. A schematic plot is provided in Figure 1.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_fig1.png?pub-status=live)
Figure 1 The energy-momentum relation
$E_{\alpha }(P)$
is expected to have two characteristic regimes: The parabolic quasi-particle regime for
$| P | \lesssim P_c (\alpha ) = \sqrt {2M^{\mathrm{eff}}(\alpha )} / \alpha $
and the flat radiative regime for larger momenta. For the transition between the two there is no precise prediction. The dashed lines denote the quasi-particle energy (1.21) and the bottom of the essential spectrum (1.23). Their intersection defines the momentum scale
$ P_c( \alpha )$
that is proportional to
$\alpha $
for large coupling. Note that the y-axis is measured in units of order
$ \alpha ^{-2}$
.
One aspect of this work is to show that the quasi-particle picture is mathematically rigorous, insofar as it provides a parabolic upper bound on
$E_{\alpha }(P)$
that coincides with the expected form of the quasi-particle energy in the limit of large coupling. Since the quasi-particle energy (1.21) is determined by the values of
$E_{\alpha }(0)$
and
$M^{\mathrm{eff}}(\alpha )$
, it is instructive to recall two long-standing open conjectures concerning their behavior for
$\alpha \to \infty $
. As explained in the previous section, the phonon field behaves classically for large coupling, and thus it is expected that
$M^{\mathrm{eff}} (\alpha )$
should asymptotically tend to the expression that follows from the corresponding semiclassical counterpart of the problem. This semiclassical theory of the effective mass was introduced by Landau and Pekar in 1948 [Reference Landau and Pekar31], and, based on this work (see also [Reference Spohn57, Reference Feliciangeli, Rademacher and Seiringer14]), it is conjectured that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn24.png?pub-status=live)
Although this problem is many decades old, the best rigorous result available at the time of writing is that
$M^{\mathrm{eff}} (\alpha )$
is divergent [Reference Lieb and Seiringer38] at least as fast as
$\alpha ^{2/5}$
[Reference Betz and Polzer4]. Regarding the ground state energy
$E_{\alpha }(0)$
the prediction from the physics literature (see, e.g., [Reference Allcock2, Reference Miyake43, Reference Tjablikow59, Reference Gross27]) is that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn25.png?pub-status=live)
for some
$\delta>0$
(in fact it is predicted that
$\delta =2$
[Reference Gross27]). Compared to the semiclassical expansion this includes a subleading correction of order
$\alpha ^{-2}$
, which we call the Bogoliubov energy, and which arises from quantum fluctuations of the field around its classical value. For a nice heuristic derivation of this correction, we recommend the study of [Reference Miyake43]. An upper bound of the form (1.25) is an immediate consequence of the results in this paper. We also note that a corresponding lower bound on
$E_{\alpha }(0)$
was recently established in [Reference Brooks and Seiringer8].
Now inserting equations (1.24) and (1.25) into equation (1.21), and based on the expectation that the quasi-particle regime is restricted to
$|P| \lesssim \sqrt {2M^{\mathrm{eff}}(\alpha )} / \alpha \sim \alpha $
, it is clear that the Bogoliubov energy needs to be taken into account in order to see the quasi-particle energy shift given by
$P^2/( 2\alpha ^4M^{\mathrm{LP}}) \le \alpha ^{-2}$
.
To put it concisely, we can summarize the heuristics discussed above in the claim that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn26.png?pub-status=live)
Our main result, Theorem 2.1 below, provides an upper bound for
$E_{\alpha }(\alpha P)$
that is compatible with this claim. To be more precise, our result implies that the left side of equation (1.26), with the limit replaced by the
$\limsup $
, is bounded from above by the expression on the right side. This shows in particular that the corrections to the quasi-particle energy are always negative, a conclusion that is not entirely obvious a priori.
During the publication process of this work, a corresponding lower bound on
$E_{\alpha }(\alpha P)$
was presented in [Reference Brooks and Seiringer9]. When combined with the upper bound presented in this work, the two bounds establish the validity of equation (1.26).
Furthermore, we would like to mention the recent progress made in the analysis of the large coupling limit of the effective mass, as reported in [Reference Brooks and Seiringer9, Reference Sellke55]. These advancements represent a significant step forward towards solving the Landau–Pekar conjecture (1.24).
Remark 1.2. An immediate consequence of equation (1.26) is that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn27.png?pub-status=live)
which is to be compared with equation (1.24) where the limits are taken in reversed order. In light of this, we interpret equation (1.26) as an additional confirmation of the polaron’s quasi-particle nature, which complements the picture suggested by equation (1.24).
Remark 1.3. An unresolved problem of interest is whether the ground state energy
$E_{\alpha }(P)$
enters the essential spectrum for some finite momentum P, which may depend on the dimension and possibly also on the value of
$\alpha $
. It is known that in two dimensions
$E_{\alpha }(P)$
remains an isolated eigenvalue for all P, meaning that the curve approaches
$\inf \sigma _{\mathrm{ess}}(H_{\alpha }(P))$
only in the limit as
$|P|$
goes to infinity [Reference Spohn58]. However, in three dimensions the question is unsettled. In Corollary 2.2 below, we prove that for large
$\alpha $
,
$E_{\alpha }(P)$
remains an isolated eigenvalue for all
$|P| \lesssim \sqrt {2 M^{\mathrm{LP}}} \alpha $
. Nonetheless, there is evidence from results obtained for weak coupling that
$E_{\alpha }(P)$
agrees with the bottom of the essential spectrum when
$|P|$
is sufficiently large [Reference Dahlbæk10].
2 Main result
We are now ready to state the main result.
Theorem 2.1. Let
$E_{\alpha }(P) = \inf \sigma (H_{\alpha }(P))$
and
$M^{\mathrm{LP}} = \frac {2}{3} \vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
with
$\varphi $
defined in equation (1.14). For every
$\varepsilon>0$
, there exists a constant
$C_{\varepsilon }> 0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn28.png?pub-status=live)
for all
$P\in {\mathbb R}^3$
and all
$\alpha $
large enough.
As a consequence of Lemma 1.1, the operator
$\sqrt { H^{\mathrm{Pek}} } - 1 $
is trace-class and nonpositive, implying that the second term on the right-hand side is finite and lowers the energy. This term corresponds to the quantum corrections to the ground state energy of the Fröhlich Hamiltonian, as discussed in Section 1.3. Since
$E_{\alpha }(0) = \inf \sigma (H_{\alpha })$
, our theorem implies a two-term upper bound for the ground state energy of the Fröhlich Hamiltonian. A complementary lower bound for
$E_{\alpha }(0)$
has been recently proved in [Reference Brooks and Seiringer8].
The result for
$|P| / \alpha \ge \sqrt {2 M^{\mathrm{LP}}}$
can be obtained from equations (1.22) and (2.1) for
$P=0$
. The relevant range for the momentum dependent term is
$|P| / \alpha \le \sqrt {2 M^{\mathrm{LP}}}$
. For momenta satisfying
$\alpha ^{-\frac {1}{4}+\frac {\varepsilon }{2}} \ll |P|/\alpha \le \sqrt {2 M^{\mathrm{LP}}} $
, the last term in equation (2.1) is subleading for large
$\alpha $
when compared to the momentum dependent term. In this region, the upper bound describes a quadratic dispersion relation for a free quasi-particle with mass
$\alpha ^4 M^{\mathrm{LP}}$
. The lower restriction
$ |P|/\alpha \gg \alpha ^{-\frac {1}{4}+\frac {\varepsilon }{2}} $
could in principle be improved by deriving a better error term in equation (2.1).
For a long time, the only rigorous lower bound available for nonzero P was the one derived by Lieb and Yamazaki [Reference Lieb and Yamazaki40] in 1958, which states that
$E_{\alpha }(P) \ge c_1 e^{\mathrm{Pek}} + c_2 P^2 / ( 2 \alpha ^4 M^{\mathrm{LP}} ) $
with
$c_1 \approx 3.07$
and
$c_2 \approx 0.11$
(where
$e^{\mathrm{Pek}}$
is negative). After the completion of our paper, a lower bound that matches our upper bound was obtained in [Reference Brooks and Seiringer9]. It should be noted that the approach in [Reference Brooks and Seiringer9] is different from ours, as it does not utilize the fiber decomposition of
$H_{\alpha }$
.
Combining a suitable lower bound on the bottom of the essential spectrum (1.23) with Theorem 2.1 yields an extension of known results regarding the existence of a unique ground state of
$H_{\alpha }(P)$
. Fröhlich [Reference Fröhlich21] showed that a unique ground state exists for
$|P| < \sqrt {2}$
, which was later extended by Spohn [Reference Spohn58] to a larger but unspecified domain. More recently, Polzer [Reference Polzer52] established the existence of the ground state for all
$|P| < \sqrt {2M^{\mathrm{eff}}(\alpha )}/\alpha $
with
$M^{\mathrm{eff}}(\alpha )$
defined by equation (1.22). Our new result demonstrates that the ground state exists for
$|P| < \sqrt {2M^{\mathrm{LP}}}\alpha $
.
Corollary 2.2. For every
$s \in (0, \tfrac {1}{29})$
, there exists a constant
$\alpha (s)>0$
such that
$E_{\alpha }(P)$
is a nondegenerate eigenvalue of
$H_{\alpha }(P)$
for all
$|P| < ( 1 - 2 \alpha ^{-s})^{1/2} \sqrt { 2 M^{\mathrm{LP}}} \alpha $
and
$\alpha \ge \alpha (s)$
.
To prove this statement, we combine equations (2.1) and (1.21) and [Reference Brooks and Seiringer8, Theorem 1.1] to show that
$E_{\alpha }(P) < \inf \sigma _{\mathrm{ess}}(H_{\alpha }(P))$
for the specified values of
$|P|$
and
$\alpha $
. This implies that
$E_{\alpha }(P)$
is part of the discrete spectrum, meaning it is an isolated eigenvalue of finite multiplicity. The nondegeneracy of this eigenvalue can then be established using a Perron–Frobenius type argument, as shown in [Reference Miyao44, Reference Lampart, Mitrouskas and Myśliwy29].
In the next two sections, we provide the definition of a suitable trial state and formulate our main statement as an energy estimate for this trial state. The remainder of the paper is devoted to the proof of this energy estimate. A sketch of the strategy of the proof is given in Section 3.2.
2.1 Bogoliubov Hamiltonian
In this section, we introduce and discuss a quadratic Hamiltonian defined on the Fock space. For its definition, we set
$\Pi _0$
and
$\Pi _1$
to be the orthogonal projectors onto
$\mathrm {Ker}H^{\mathrm{Pek}} = \text {Span}\{ \partial _i \varphi : i=1,2,3\} $
and
$(\mathrm {Ker}H^{\mathrm{Pek}})^{\perp }$
, that is,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn29.png?pub-status=live)
Even though we will not make explicit use of it, it is convenient to keep in mind that the decomposition
$ L^2 ({\mathbb R}^3) \, =\, \text {Ran}(\Pi _0) \oplus \text {Ran}(\Pi _1) $
implies the factorization
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn30.png?pub-status=live)
For technical reasons, which are explained in Section 3.4.3, we introduce the Bogoliubov Hamiltonian
${\mathbb H}_K $
with a momentum cutoff
$K\in (0,\infty ]$
. Setting
${\mathbb N}_1 = \mathrm {d} \Gamma (\Pi _1)$
(the number operator on
${\mathcal F}_1$
) and recalling equation (1.18) we define
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn31.png?pub-status=live)
where the new coupling function
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn32.png?pub-status=live)
results from the coupling function
$h_{x}$
by removing all momenta larger than K and then projecting to
$\text {Ran}(\Pi _1)$
. The second term in equation (2.4) defines the quadratic operator given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn33.png?pub-status=live)
By definition,
${\mathbb H}_K$
acts nontrivially only on the tensor component
${\mathcal F}_1$
. Below we will show that
${\mathbb H}_K$
is bounded from below and diagonalizable by a unitary Bogoliubov transformation. For the precise statement, we need some further preparations.
For
$K\in (0,\infty ]$
, we introduce
$H^{\mathrm{Pek}}_K$
as the operator on
$L^2({\mathbb R}^3)$
defined by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn34.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn35.png?pub-status=live)
where
$T_K$
is defined by the integral kernel
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn36.png?pub-status=live)
By definition
$H^{\mathrm{Pek}}_{\infty } = H^{\mathrm{Pek}}$
; see equation (1.20). Moreover, we set
$\Theta _K = (H^{\mathrm{Pek}}_K)^{1/4}$
and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn37.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn38.png?pub-status=live)
The next lemma, whose proof can be found in Section 4, implies some useful properties of these operators, among others, that there are constants
$C,K_0>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn39.png?pub-status=live)
Lemma 2.3. For
$K_0$
large enough, there exist constants
$\beta \in (0,1)$
and
$C>0$
such that for all
$K\in (K_0,\infty ]$
-
(i)
$0 \le H^{\mathrm{Pek}}_K \le 1$ and
$ (H^{\mathrm{Pek}}_K - \beta ) \restriction \mathrm {Ran} (\Pi _1) \ge 0$
-
(ii)
$(B_K)^2 \le C( 1- H_K^{\mathrm{Pek}})$
-
(iii)
$ {\mathrm {Tr}}_{L^2}(1- H_K^{\mathrm{Pek}}) \le C $ .
Moreover, for all
$K\in (K_0, \infty )$
-
(iv)
${\mathrm {Tr}}_{L^2}((-i\nabla )( 1- H_K^{\mathrm{Pek}} )(-i\nabla )) \le C K $ .
Remark 2.4. Since
$H^{\mathrm{Pek}}_K$
has a real-valued kernel it satisfies
$ H^{\mathrm{Pek}}_K f = \operatorname {Re}( H^{\mathrm{Pek}}_K f) + i \operatorname {Im} (H^{\mathrm{Pek}}_Kf) $
for all
$f\in L^2({\mathbb R}^3)$
, and the same holds for
$\Pi _0$
and
$\Pi _1$
. By the spectral calculus for self-adjoint operators, this property extends to
$\Theta _K$
and
$\Theta _K^{-1}$
.
To make the relation between
${\mathbb H}_K$
and
$H^{\mathrm{Pek}}_K$
precise, we introduce the transformation
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn40.png?pub-status=live)
That this transformation defines a unitary operator
$\mathbb U_K$
for all
$K\in (K_0,\infty ]$
is a consequence of equation (2.10) and
$A_K^2 = 1 + B_K^2$
by the well-known Shale–Stinespring condition (see [Reference Ruijsenaars54, Reference Shale56, Reference Solovej53]). Also, note that
$\mathbb U_K$
does not mix the two factors in
${\mathcal F} = {\mathcal F}_0 \otimes {\mathcal F}_1$
.
Lemma 2.5. For
$K\in (K_0,\infty ]$
with
$K_0$
large enough and
$\mathbb U_K$
, the unitary operator defined by equation (2.11), we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn41.png?pub-status=live)
with
$H^{\mathrm{Pek}}_K$
defined by equations (2.7a) and (2.7b).
The proof is obtained by an explicit computation and postponed to Section 4. From this lemma, we can infer that the ground state energy of
${\mathbb H}_K$
is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn42.png?pub-status=live)
where we also used
$\Pi _1 = 1-\Pi _0$
and
$\text {Tr}_{L^2}(\Pi _0)= 3$
. Moreover, since
$H^{\mathrm{Pek}}_K \le \Pi _1$
we have
$\inf \sigma ( {\mathbb H}_K) < 0$
and from item (iv) of Lemma 2.3 we find that
$\inf \sigma ( {\mathbb H}_K)> - \infty $
uniformly in
$K\to \infty $
.
For the ground state of
, we shall use the notation
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn43.png?pub-status=live)
where it is important to keep in mind that the state
$\Upsilon _K $
has excitations only in
${\mathcal F}_1$
(i.e., no zero-mode excitations) since
$\mathbb U_K^{\dagger }$
acts as the identity on
${\mathcal F}_0$
; see equation (2.9b).
From now on, we shall always assume
$K\ge K_0$
large enough such that Lemmas 2.3 and 2.5 are applicable.
2.2 Trial state and energy estimate
As starting point for the definition of our trial state, consider the Fock space wave function obtained from the fiber decomposition of the classical Pekar product state
$\Psi _{\psi ,\varphi }$
, that is,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn44.png?pub-status=live)
Testing the energy of
$H_{\alpha }(P)$
with
$\Psi ^{\mathrm{Pek}}_{\alpha }(P)$
, one would in fact obtain that
$E_{\alpha }(P)$
is bounded from above by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn45.png?pub-status=live)
For
$E_{\alpha }(0)$
, this provides already a better bound compared to the semiclassical approximation for
$\inf \sigma (H_{\alpha })$
. The improvement comes from taking into account the translational symmetry and can be interpreted as the missing zero-point energy of three quantum oscillators (that turned into translational degrees of freedom). As a side remark, we find it somewhat surprising that fiber decompositions of this form have been employed very rarely in the polaron literature, exceptions being [Reference Höhler and Caianiello28] and [Reference Nagy49]. We think they could be of interest also for other translation-invariant polaron type models.
To obtain the desired bound for
$E_{\alpha }(P)$
, we need to add several modifications to the integrand in equation (2.15). On the one hand, we have to replace the classical field
$\varphi $
by a suitably shifted
$\varphi _P$
in order to get the correct momentum dependent term (note that equation (2.16) is missing a factor
$\frac {1}{2}$
in the quadratic term). The missing part of the rest energy (compare with equation (2.13)), on the other hand, is caused by two types of correlations that need to be added to the Pekar product state. First, we include correlations between the electron and the phonons. This is done in the spirit of first-order adiabatic perturbation theory. Second, we rotate the vacuum by the unitary transformation (2.11) that diagonalizes the Bogoliubov Hamiltonian (2.4). As discussed, the latter describes the quantum fluctuations of the phonons around the classical field. For technical reasons, briefly explained in Section 3.2, we also need to introduce suitable momentum and space cutoffs in the trial state.
Explicitly, we consider the family of Fock space wave functions
$\Psi _{K,\alpha }(P) \in {\mathcal F}$
, depending on the coupling
$\alpha $
, the total momentum
$P\in {\mathbb R}^3$
and the cutoff
$K \in (K_0,\infty )$
, given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn46.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn47.png?pub-status=live)
and (recall equations (1.18) and (2.5))
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn48.png?pub-status=live)
Here,
$u_{\alpha } :{\mathbb R}^3 \to [0,1]$
is a radial function, satisfying
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn49.png?pub-status=live)
for some
$C>0$
. For completeness, we recall that
$\psi> 0$
and
$\varphi $
are the unique rotational invariant minimizers of the Pekar functionals (1.12) and (1.13).
Remark 2.6. Writing
$G_{K,x}^{i}$
, we think of these states as elements in
$L^2({\mathbb R}^3, {\mathcal F})$
and of
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn50.png?pub-status=live)
as an x-dependent Fock space operator. Via the isomorphism
$L^2({\mathbb R}^3, {\mathcal F}) \simeq \mathscr H$
, we can view
$G_{K,x}^{i}$
also as a wave function in
$ \mathscr H$
. In this case, we shall write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn51.png?pub-status=live)
Remark 2.7. Let us note that in equation (2.17), we anticipated the fact that the integrand is in
$ L^1 ({\mathbb R}^3,{\mathcal F})$
and thus
$\Psi _{K,\alpha }(\cdot ) \in C_b({\mathbb R}^3,{\mathcal F})$
. For
$G_{K}^0$
, the integrability follows directly from the exponential decay of
$\psi $
(as shown in Lemma 3.7), while for
$G_{K}^1$
it can be seen from
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn52.png?pub-status=live)
where we used Cauchy–Schwarz and Lemmas 3.8, 3.9 and 3.13. A more precise estimate for the norm of
$\Psi _{K,\alpha }(P)$
for large
$\alpha $
will be given in Proposition 3.17.
For the introduced trial states, we prove the following energy estimate, where
${\mathbb H}_{\infty }$
denotes the Bogoliubov Hamiltonian (2.4) for
$K=\infty $
.
Proposition 2.8. Let
$\Psi _{K,\alpha }(P)\in {\mathcal F}$
as in equation (2.17), choose
$c,\tilde c>0$
and set
$r(K,\alpha ) = K^{-1/2}\alpha ^{-2} + \sqrt K \alpha ^{-3}$
. For every
$\varepsilon>0$
, there exists a constant
$C_{\varepsilon }> 0$
(we omit the dependence on c and
$\tilde c$
) such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn53.png?pub-status=live)
for all
$|P|/\alpha \le c $
,
$ K /\alpha \le \tilde c $
and
$\alpha $
large enough.
The next section, which constitutes the bulk of the paper, is devoted to proving this proposition. Before embarking on the proof, let us now deduce its main consequence and conclude the proof of Theorem 2.1.
Proof of Theorem 2.1.
With equation (2.13) and
$H_{\infty }^{\mathrm{Pek}} = H^{\mathrm{Pek}}$
we can rewrite the term of order
$\alpha ^{-2}$
as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn54.png?pub-status=live)
Choosing K proportional to
$\alpha $
optimizes the asymptotics of the error in equation (2.24) and thus proves equation (2.1) for
$|P|\le \sqrt {2M^{\mathrm{LP}}} \alpha $
by the variational principle. For larger
$|P|$
, we use
$E_{\alpha }(P) \le E_{\alpha }(0) + \alpha ^{-2}$
as a consequence of equation (1.24) and apply equation (2.1) for
$P=0$
.
3 Proof of Proposition 2.8
We recall the definition of the field operators
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn55.png?pub-status=live)
and the Weyl operator
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn56.png?pub-status=live)
The Weyl operator is unitary and satisfies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn57.png?pub-status=live)
3.1 The total energy
The proof of Proposition 2.8 starts with a convenient formula for the energy evaluated in the trial state. For the precise statement, we introduce the y-dependent function in
$L^2({\mathbb R}^3)$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn58.png?pub-status=live)
and the y-dependent Fock space operator
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn59.png?pub-status=live)
Since
$g_{P}(y)$
is real-valued, we have
$(A_{P,y})^{\dagger } = - A_{P,y}$
.
We further consider the shift operator
$T_y = e^{y\nabla }$
on
$L^2({\mathbb R}^3)$
, that is,
$(T_y f)(x) = f(x+y)$
for every
$f\in L^2({\mathbb R}^3)$
, and the Hamiltonian acting on
$\mathscr H$
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn60.png?pub-status=live)
where we recall that
$h^{\mathrm{Pek}} = -\Delta + V^{\varphi } - \lambda ^{\mathrm{Pek}}$
.
Lemma 3.1. For
$\Psi _{K,\alpha }(P)$
defined in equation (2.17), we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn61.png?pub-status=live)
where
$\mathcal N = \vert \hspace {-1pt}\vert \Psi _{K,\alpha }(P)\vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(\mathcal F\)}}}$
and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn62.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn63.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn64.png?pub-status=live)
For the proof, we recall that the Weyl operator shifts the creation and annihilation operators by complex numbers,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn65.png?pub-status=live)
and, as a simple consequence,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn66.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn67.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn68.png?pub-status=live)
Moreover, we need the following identity.
Lemma 3.2. Let
$\varphi _P = \varphi + i \xi _P$
with
$\xi _P$
defined by (2.18). Then
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn69.png?pub-status=live)
Proof of Lemma 3.2.
We first observe that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn70.png?pub-status=live)
which follows from
$\frac {d}{ds} e^{- i s P_f y} a^{\dagger }( e^{(s-1)y\nabla } f ) e^{isP_fy} = 0$
. In combination with equation (3.3), this leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn71.png?pub-status=live)
Recalling
$\varphi _P = \varphi + i \frac {1}{\alpha ^2 M^{\mathrm{LP}}} (P \nabla ) \varphi $
, we compute
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn72.png?pub-status=live)
where we inserted
$e^{-y\nabla } = 1 - (y\nabla ) + \frac {1}{2}( y \nabla )^2 - \int _0^1 \mathrm {d} s\, e^{-s y\nabla } (y\nabla )^3$
and used that, due to rotational invariance of
$\varphi $
,
$ \langle \varphi | (P\nabla ) \varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = 0 = \langle \varphi | (y \nabla )^2 (P\nabla ) \varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} $
. Also, because of rotational invariance,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn73.png?pub-status=live)
and thus,
$ \alpha ^2 \operatorname {Im} \langle \varphi _P | e^{-y \nabla } \varphi _P \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = P y+ g_{P}(y)$
.
Proof of Lemma 3.1.
Throughout this proof, let
$\Xi _i = W ( \alpha \varphi _P ) G^{i}_{K} \in L^2({\mathbb R}^3) \otimes {\mathcal F} $
,
$i=0,1$
, with
$G_{K}^i$
defined in equation (2.19) and set
$\Psi _i = \int \mathrm {d} x \, e^{i (P_f-P)x} \Xi _i(x)$
. First, note that
$\Psi _i \in D(H_{\alpha }(P)^{1/2} ) $
for
$i,j\in \{ 0,1\}$
which follows from
$D(H_{\alpha }(P)^{1/2}) = D( |P_f| + {\mathbb N}^{1/2} ) $
[Reference Lieb and Yamazaki40] together with
$| P_f | \Xi _i \in L^1({\mathbb R}^3,{\mathcal F})$
and
${\mathbb N}^{1/2} \Xi _i \in L^1({\mathbb R}^3,{\mathcal F})$
. The integrability of these states is verified using Lemmas 3.16 and 3.14.
Below, we shall employ the identitiesFootnote 1
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn74.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn75.png?pub-status=live)
where
$H_{\alpha }$
is the Fröhlich Hamiltonian given by equation (1.2). To obtain the first identity, write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqnu1.png?pub-status=live)
and for the second one, use
$e^{-i(P_f-P)x} \phi (v_0) = \phi (v_x) e^{-i(P_f-P)x} $
so that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn76.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn77.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn78.png?pub-status=live)
With equations (3.16a) and (3.16b), the norm and the energy of the trial state are given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn79.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn80.png?pub-status=live)
Inserting
$\Xi _i = W ( \alpha \varphi _P ) G^{i}_{K}$
and applying Lemma 3.2, we find for
$i,j\in \{0,1\}$
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn81.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn82.png?pub-status=live)
Using equations (3.10a) and (3.10b) and
$2\operatorname {Re} \langle \varphi _P | h_{x} \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = 2\operatorname {Re} \langle \varphi | h_{x} \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = V^{\varphi }(x)$
(see equation (1.11)) the Weyl-transformed Hamiltonian becomes
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn83.png?pub-status=live)
with
$\widetilde H_{\alpha ,P}$
defined by equation (3.6). Note that we added and subtracted
$e^{\mathrm{Pek}}= \lambda ^{\mathrm{Pek}} + \vert \hspace {-1pt}\vert \varphi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}^2$
and used equation (3.6). Altogether, this implies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn84.png?pub-status=live)
The claimed result now follows from
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn85.png?pub-status=live)
where we used
$\vert \hspace {-1pt}\vert (P\nabla ) \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}} = \frac {P^2}{3} \vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}} = \frac {P^2}{2} M^{\mathrm{LP}}$
because of rotational invariance of
$\varphi $
.
3.2 A short guide to the proof
3.2.1 Heuristic picture
Given Lemma 3.1, the remaining task is to show that
$({\mathcal E}+{\mathcal G} + {\mathcal K})/\mathcal N$
coincides, up to small errors, with the energy contribution of order
$\alpha ^{-2}$
in equation (2.24). Although our proof is somewhat technical, the main idea is a simple one, and we explain the corresponding heuristics here in order to facilitate the reading. The main point is that the integrals appearing in the terms given in Lemma 3.1 turn out to be, as
$\alpha \rightarrow \infty $
and
$|P|/ \alpha \leq c$
, sharply localized around zero at the length scale of order
$ \alpha ^{-1}$
. In this regime, as formally
$w_{P,y}(z)\approx y\nabla \varphi (z) $
for y small, the Weyl operator
$W(\alpha w_{P,y})$
effectively acts nontrivially only on the
$\mathcal {F}_0$
part of the Fock space (at this point, it is convenient to recall the factorization (2.3)). Moreover, we shall show that
$e^{A_{P,y}}$
can be effectively replaced by the identity operator and it suffices to consider
$T_y \approx 1+ y\nabla $
. Since our trial state coincides with the vacuum on
$\mathcal {F}_0$
, we thus expect for
$|y|$
small that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn86.png?pub-status=live)
with
$\lambda =\frac {1}{6}\vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}^2$
. (Since
$T_y$
acts on the electron coordinate, it commutes with
$e^{A_{P,y}}$
and
$W(\alpha w_{P,y})$
). Taking this approximation for granted, and considering only the term with
$i=j=0$
in equations (3.18a) and (3.19a), would lead to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn87.png?pub-status=live)
With the above replacement and keeping only the terms of order
$\alpha ^{-2}$
(relative to the factor from the norm), the energy terms are found to be given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn88.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn89.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn90.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn91.png?pub-status=live)
From here, the Bogoliubov energy is obtained by setting
$u_{\alpha } = 1$
and
$K=\infty $
in the leading-order terms and using
$R h^{\mathrm{Pek}} R = R$
, since this would imply (omitting the errors)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn92.png?pub-status=live)
The remaining
$-\frac {3}{2\alpha ^2}$
term stems from the part of the interaction involving the zero modes. In equation (3.25b), the term not involving
$y\nabla $
vanishes due to
$\langle \psi |h_{\cdot }\psi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} =-\varphi $
. Moreover,
$\langle \psi | h_{\cdot } \nabla \psi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} =-\frac {1}{2} \nabla \varphi $
using
$\nabla h_{\cdot }=-(\nabla h)_{\cdot }$
via integration by parts (in the sense of distributions). Thus, since
$[a^{\dagger } (y \nabla \varphi ),\mathbb U_{\infty }^{\dagger }]=0$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn93.png?pub-status=live)
Equations (3.26) and (3.27) now add up to the desired energy of order
$\alpha ^{-2}$
; see equation (2.25). Note that for estimating the error induced by replacing
$e^{A_{P,y}}$
by unity we require the momentum cutoff K in the definition of the trial state; see Lemma 3.16.
The main issue in equation (3.23) is that, while for small enough y one can use the first-order approximation
$W(\alpha w_{P,y}) \approx W(\alpha y\nabla \varphi )$
, for y large, on the other hand, the higher-order terms in
$w_{P,y}$
begin to play an important part, ultimately killing the Gaussian factor. Writing
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn94.png?pub-status=live)
we notice that, since
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn95.png?pub-status=live)
the prefactor should lead to a y-independent, exponentially small constant. In order to make use of this exponential decay in
$\alpha $
, however, we need to ensure that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn96.png?pub-status=live)
is polynomially bounded in
$\alpha $
with some integrable function
$g(y)$
, which heuristically can be expected to be true since the average number of particles in the state
$\widetilde H_{\alpha }(P) G_K^i$
is of order one w.r.t.
$\alpha $
. To obtain the required integrability in y is also the reason for introducing the cutoff function
$u_{\alpha }$
in the definition of
$G^1_K$
.
3.2.2 Outline of the proof
Although the replacement (3.23) illustrates the main idea behind extracting the leading order terms, in our proof we do not directly perform this replacement and estimate the resulting error. Instead, when taking inner products, we commute the exponential operators
$e^{a^{\dagger }(\alpha w_P)}$
and
$e^{-a(\alpha w_P)}$
in
$W(\alpha w_{P,y})$
to the left resp. to the right until they hit the vacuum state in
$G_K^i$
. This involves the Bogoliubov transformation, cf. Lemma 3.12 and gives rise to a slight modification of
$w_{P,y}$
, which we denote by
$\widetilde {w}_{P,y}$
. These manipulations naturally lead to a multiplicative factor
$\exp (-\frac {\alpha ^2}{2}\vert \hspace {-1pt}\vert \widetilde {w}_{P,y} \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}} )$
which, as we shall see, indeed behaves like the Gaussian function in equation (3.23) for
$|y|$
small and tends to a constant exponentially small in
$\alpha $
as
$|y| \rightarrow \infty $
. In Lemma 3.5, we prove the large
$\alpha $
asymptotics of integrals of the type
$\int g(y)\exp (-\frac {\alpha ^2}{2}\vert \hspace {-1pt}\vert \widetilde {w}_{P,y}\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}^2 ) \mathrm {d} y$
for a suitable class of functions g. The major part of the proof, apart from extracting the leading order terms, is to establish that the resulting error terms in the integrands are, in fact, functions in this class. This is, for the most part, achieved by use of elementary estimates combined with the commutator method by Lieb and Yamazaki [Reference Lieb and Yamazaki40] in the form stated in Lemma 3.9. As already mentioned, for certain terms this makes the introduction of the space cutoff
$u_{\alpha }$
and the momentum cutoff K necessary, while for others, it is enough to use the well-known regularity properties of
$\psi $
, the relevant consequences of which are summarized in Lemma 3.7.
In the next two sections, we state the remaining necessary lemmas. The main proof is then carried out in Sections 3.5–3.9.
Throughout the remainder of the proof, we will abbreviate constants by the letter C and write
$C_{\tau }$
whenever we want to specify that it depends on a parameter
$\tau $
. As usual, the value of a constant may change from one line to the next.
3.3 The Gaussian lemma
We recall that
$w_{P,y} = (1-e^{-y\nabla })\varphi _P$
and
$\Theta _K = (H^{\mathrm{Pek}}_K)^{1/4} $
and set
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn97.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn98.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn99.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn100.png?pub-status=live)
Remark 3.3. Note that
$(y,z)\mapsto \operatorname {Re}(w_{P,y})(z)$
is even as a function on
$\mathbb {R}^6$
, while
$\operatorname {Im} (w_{P,y})(z)$
is odd on the same space. Since
$\Pi _0$
and
$\Theta _K$
both commute with the reflection operator
$(\pi f) (x) = f(-x)$
, they preserve the parity properties just mentioned. That
$\Pi _0$
has the desired properties follows directly from its explicit form. To see this for
$\Theta _K$
, it is enough to check this for
$H_K^{\mathrm{Pek}}$
, which can be easily done using the fact that the resolvent R commutes with the reflection operator, which, on the other hand, follows from the invariance of
$h^{\mathrm{Pek}}$
and
$P_{\psi }$
under parity, cf. the definition of R (1.18). Thus,
$(y,z) \mapsto \operatorname {Re}(w^i_{P,y})(z)$
is even as a function on
$\mathbb {R}^6$
for
$i=0,1$
while the corresponding imaginary parts are odd on the same space. These facts will be of relevance below where they lead to the vanishing of several integrals.
The following lemma is proven in Section 4.
Lemma 3.4. Let
$\lambda = \frac {1}{6}\vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
. For every
$c>0$
, there exists a constant
$C>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn101.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn102.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn103.png?pub-status=live)
for all
$y\in {\mathbb R}^3$
,
$|P|/\alpha \le c$
and
$\alpha>0$
.
For
$0\le \delta <1$
and
$\eta>0$
, we introduce the weight function
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn104.png?pub-status=live)
where, for ease of notation, the dependence on
$\alpha $
, P and K is omitted. Using the arguments laid down in Remark 3.3, it is easy to see that
$n_{\delta , \eta }(y)$
is even as a function of y. Moreover, in the limit of large
$\alpha $
the dominant part of the weight function when integrated against suitably decaying functions comes from the term in the exponent that is quadratic in y, cf. equation (3.32c). This is a crucial ingredient in our proofs and the content of the next lemma.
Lemma 3.5. Let
$\eta _0> 0$
,
$c>0$
,
$\lambda = \frac {1}{6}\vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
and
$n_{\delta ,\eta }$
defined in equation (3.33). For every
$n\in {\mathbb N}_0$
, there exist constants
$ d, C_n>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn105.png?pub-status=live)
for all nonnegative functions
$g\in L^{\infty }({\mathbb R}^3) \cap L^1({\mathbb R}^3)$
,
$\eta \ge \eta _0$
,
$\delta \in [0,1) $
,
$|P|/\alpha \le c $
and all
$\alpha $
large enough.
At first reading, one should think of
$n=0$
,
$\delta = 0$
,
$\eta = 1$
and g a suitable
$\alpha $
-independent nonnegative function. In this case, the integral involving the Gaussian is of order
$\alpha ^{-3}$
whereas the term on the right-hand side is of order
$\alpha ^{-4}$
and thus contributing a subleading error. The proof of the lemma is given in Section 4. As a direct consequence that will be useful to estimate error terms, we find
Corollary 3.6. Given the same assumptions as in Lemma 3.5, for every
$n \in {\mathbb N}_0$
there exist constants
$d, C_n>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn106.png?pub-status=live)
for all nonnegative functions
$g\in L^{\infty }({\mathbb R}^3) \cap L^1({\mathbb R}^3)$
,
$\eta \ge \eta _0 $
,
$\delta \in [0,1) $
,
$|P|/\alpha \le c$
and all
$K,\alpha $
large enough.
Proof of Corollary 3.6.
Since
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn107.png?pub-status=live)
the statement follows immediately from Lemma 3.5.
3.4 Further preliminaries
3.4.1 Estimates involving the Pekar minimizers
Lemma 3.7. Let
$\psi>0$
be the (normalized) rotational invariant unique minimizer of the Pekar functional (1.12), and let
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn108.png?pub-status=live)
We have that
$\psi $
,
$|\nabla \psi |$
and H are
$L^p(\mathbb {R}^3,(1+|x|^n)\mathrm {d} x)$
functions for all
$1\leq p \leq \infty $
and all
$n\geq 0$
. Moreover, there exists a constant
$C>0$
such that for all x we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn109.png?pub-status=live)
Proof. As follows from [Reference Lieb36],
$\psi (x)$
is monotone decreasing in
$|x|$
; moreover, it is smooth and bounded and vanishes exponentially at infinity, that is, there exists a constant
$C>0$
such that
$\psi (x)\leq Ce^{-|x|/C}$
for all
$|x|$
large enough (for the precise asymptotics see [Reference Moroz and Schaftingen46]). This clearly implies the statement for
$\psi $
. It further implies that all the derivatives of
$\psi $
are bounded. Hence, in order to show the desired result for
$| \nabla \psi |$
, it suffices to show that
$\int \mathrm {d} x |x|^n | \nabla \psi (x)|$
is finite for all
$n\geq 0$
. Since
$\psi $
is radial, that is, there is a function
$\psi ^{\mathrm{rad}}:[0,\infty )\to (0,\infty )$
such that
$\psi (x) = \psi ^{\mathrm{rad}}(|x|)$
, and monotone decreasing, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn110.png?pub-status=live)
for all
$R_0>0$
. Clearly, H is bounded, and hence, by
$|x+y|^n\leq 2^{n-1}\left (|x|^n+|y|^n\right )$
, we can easily bound
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn111.png?pub-status=live)
from which the statement follows also for H. To show (3.38), use the Fourier representation
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn112.png?pub-status=live)
together with
$H(x)\leq 1$
,
$\cos (2\pi kx)\geq 1-\frac {( kx)^2}{2}$
,
$\| \psi \|_{\scriptscriptstyle{\text{\(L^2\)}}} =1 $
and
$\nabla \psi \in L^2$
.
The next lemma contains suitable bounds for the potential
$V^{\varphi }$
and the resolvent R introduced in equations (1.11), (1.14) and (1.18).
Lemma 3.8. There is a constant
$C>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn113.png?pub-status=live)
Proof. For the proof of the first two inequalities, we refer to [Reference Leopold, Rademacher, Schlein and Seiringer35, Lemma III.2]. The bound for the resolvent is obtained through
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn114.png?pub-status=live)
where we made use of the second inequality in equation (3.42).
3.4.2 The commutator method
In the course of the proof, we are frequently faced with bounding field operators like
$\phi (h_x)$
. From the standard estimates for creation and annihilation operators, we would obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn115.png?pub-status=live)
which is not sufficient since
$h_{0}(y)$
is not square-integrable. With the aid of the commutator method introduced by Lieb and Yamazaki [Reference Lieb and Yamazaki40] one obtains suitable upper bounds by using in addition some regularity in the electron variable of the wave function
$\Psi $
. For our purpose, the version summarized in the following lemma will be sufficient.
Lemma 3.9. Let
$h_{K,\cdot }$
for
$K\in (1,\infty ]$
as defined in equation (2.5), let A denote a bounded operator in
$L^2({\mathbb R}^3)$
(acting on the field variable) and
$a^{\bullet } \in \{ a , a^{\dagger } \}$
. Further, let
$X,Y$
be bounded symmetric operators in
$L^2({\mathbb R}^3)$
(acting on the electron variable) that satisfy
$D_0 := \vert \hspace {-1pt}\vert X \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} \vert \hspace {-1pt}\vert Y \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} + \vert \hspace {-1pt}\vert \nabla X \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}}\vert \hspace {-1pt}\vert Y \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} + \vert \hspace {-1pt}\vert X \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} \vert \hspace {-1pt}\vert \nabla Y \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}}< \infty $
. There exists a constant
$C>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn116.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn117.png?pub-status=live)
for all
$y\in {\mathbb R}^3$
,
$\Psi \in \mathscr H$
and
$1 \le K \le \Lambda \le \infty $
.
Remark 3.10. Note that
$A h_{K,\cdot +y}= T_y (Ah_{K,\cdot }) $
and in case that A has an integral kernel,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn118.png?pub-status=live)
Proof of Lemma 3.9.
To obtain the first inequality, write
$h_{K,\cdot } = (h_{K,\cdot }-h_{1,\cdot }) + h_{1,\cdot }$
and apply the second inequality (with
$\Lambda $
and K interchanged) to the term in parenthesis. The bound for the term involving
$h_{1,\cdot }$
follows from equation (3.44), as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn119.png?pub-status=live)
To verify the second inequality, write the difference as a commutator
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn120.png?pub-status=live)
and use that
$\nabla $
and A commute (they act on different variables). Then similarly as in equation (3.47) we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn121.png?pub-status=live)
The desired bound now follows from
$\sup _{\Lambda> K} \vert \hspace {-1pt}\vert j_{K,\Lambda ,0} \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}^2 \le C/K$
.
A simple but useful corollary is given by
Corollary 3.11. Under the same conditions as in Lemma 3.9, with the additional assumption that Y is a rank-one operator, there exists a constant
$C>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn122.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn123.png?pub-status=live)
for all
$y\in {\mathbb R}^3$
and
$1 \le K \le \Lambda \le \infty $
.
Proof. Since Y has rank one, we can use
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn124.png?pub-status=live)
for any
$w\in L^2({\mathbb R}^3)$
, and similarly for equation (3.50b), and apply Lemma 3.9.
3.4.3 Transformation properties of
$\mathbb U_K$
The next lemma collects some useful relations for the Bogoliubov transformation
$\mathbb U_K$
. The proof follows directly from the definition (2.11) and the properties explained in Remark 2.4. We omit the details.
Lemma 3.12. Let
$f \in L^2({\mathbb R}^3)$
,
$f^0 = \Pi _0 f$
,
$f^1 = \Pi _1 f$
with
$\Pi _i$
defined in equation (2.2), and set
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn125.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn126.png?pub-status=live)
The unitary operator
$\mathbb U_K$
defined in equation (2.11) satisfies the relations
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn127.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn128.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn129.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn130.png?pub-status=live)
Note that equation (3.31d) is consistent with the general notation introduced in equation (3.52b). The following statements provide helpful bounds involving the number operator when transformed with the Bogoliubov transformation.
Lemma 3.13. There exists a constant
$b>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn131.png?pub-status=live)
for all
$n\in {\mathbb N} $
.
Proof. With b replaced by
$b_K = 2 \vert \hspace {-1pt}\vert B_K \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{HS}\)}}}^2 + \vert \hspace {-1pt}\vert A_K\vert \hspace {-1pt}\vert _{\mathrm {op}}^2 + 1 $
, both estimates follow from [Reference Boßmann, Petrat, Pickl and Soffer7, Lemma 4.4] together with equations (3.53a) and (3.53b). That
$b_K\le b $
for some K-independent
$b>0$
is inferred from Lemma 2.3.
In the next two statements, we denote by
(resp.
) the orthogonal projector in
${\mathcal F}$
to all states with phonon number larger than (resp. less or equal to) c.
Corollary 3.14. Let
$ \Upsilon _K = \mathbb U_K^{\dagger } \Omega $
and
for
$\delta> 0$
. There exist constants
$b,C_{\delta ,n}> 0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn132.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn133.png?pub-status=live)
for all
$n\in {\mathbb N}_0$
.
Proof. The first bound follows directly from Lemma 3.13. The second one is obtained from
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn134.png?pub-status=live)
with
$m \ge 20 / \delta $
.
Lemma 3.15. For
$\delta>0$
and
$\kappa = 1/ (16 e b \alpha ^{\delta })$
with
$b> 0$
the constant from Lemma 3.13, the operator inequality
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn135.png?pub-status=live)
holds for all
$\alpha $
large enough.
Proof. We write out the Taylor series for the exponential and invoke Lemma 3.13,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn136.png?pub-status=live)
where we used
$1\le 2\alpha ^{\delta } $
in the last step. The stated bound now follows from the elementary inequality
$n! \ge (\frac {n}{e})^n$
.
The reason for introducing the momentum cutoff in
${\mathbb H}_K$
is to obtain a finite upper bound for the norm of the state
$P_f \Upsilon _K $
. This is the content of the next lemma, whose proof is given in Section 4.
Lemma 3.16. Let
$P_f = \int \mathrm {d} k \, k\, a_k^{\dagger } a_k$
and
$K_0$
large enough. There is a
$C>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn137.png?pub-status=live)
for all
$K \in (K_0, \infty )$
.
3.5 Norm of the trial state
In this section, we provide the computation of the norm
$\mathcal N = \vert \hspace {-1pt}\vert \Psi _{K,\alpha }(P)\vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(\mathcal F\)}}}$
.
Proposition 3.17. Let
$\lambda = \frac {1}{6} \vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
and
$c>0$
. For every
$\varepsilon>0$
, there exist a constant
$C_{\varepsilon }> 0$
(we omit the dependence on c) such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn138.png?pub-status=live)
for all
$|P|/\alpha \le c$
and all
$\alpha $
large enough.
Proof. It follows from equations (3.18a) and (3.19a) that
$\mathcal N = \mathcal N_{0} + \mathcal N_1 + \mathcal N_2 $
with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn139.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn140.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn141.png?pub-status=live)
$\underline {\mathrm {Term}\ {\mathcal N}_0}$
. This part contains the leading order contribution
$(\frac {\pi }{\lambda \alpha ^2})^{3/2}$
. With H defined in equation (3.37), let us write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn142.png?pub-status=live)
In the first term, we use
$\Upsilon _K = \mathbb U_K^{\dagger } \Omega $
and apply equation (3.53d) to transform the Weyl operator with the Bogoliubov transformation. This gives
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn143.png?pub-status=live)
with
$\widetilde w_{P,y}$
defined in equation (3.31d). From equations (3.2) and (3.33), we thus obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn144.png?pub-status=live)
Since
$\vert \hspace {-1pt}\vert H\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^1\)}}} + \vert \hspace {-1pt}\vert H\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^{\infty }\)}}} \le C $
, cf. Lemma 3.7, we can apply Lemma 3.5 in order to replace the weight function
$n_{0,1}(y)$
by the Gaussian
$e^{-\lambda \alpha ^2 y^2}$
. More precisely,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn145.png?pub-status=live)
for all
$|P| /\alpha \le c$
and all
$K,\alpha $
large enough. Then we use
$| H(y)-1| \le C y^2$
in order to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn146.png?pub-status=live)
To treat
$\mathcal N_{02}$
, it is convenient to decompose the state
$ \Upsilon _K $
into a part with bounded particle number and a remainder. To this end, we choose a small
$\delta>0$
and write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn147.png?pub-status=live)
Inserting this into
$\mathcal N_{02}$
and using unitarity of
$e^{A_{P,y}}$
and
$\vert \hspace {-1pt}\vert H\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^1\)}}}\le C$
, we can estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn148.png?pub-status=live)
By Corollary 3.14 for
$n=0$
,
$\vert \hspace {-1pt}\vert \Upsilon _K^>\vert \hspace {-1pt}\vert \le C_{\delta } \, \alpha ^{-10}$
. In the remaining expression, we use equation (3.63),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn149.png?pub-status=live)
and insert the identity
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn150.png?pub-status=live)
on the left of the Weyl operator (where
$b>0$
is the constant from Lemma 3.13). After applying the Cauchy–Schwarz inequality, this leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn151.png?pub-status=live)
In the second factor, we then employ
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn152.png?pub-status=live)
and use
$e^{-\kappa {\mathbb N}}a^{\dagger }(f)e^{\kappa {\mathbb N}} = a^{\dagger }(e^{-\kappa } f) $
to write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn153.png?pub-status=live)
Combining the previous two lines, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn154.png?pub-status=live)
for some
$\alpha $
-independent
$\eta> 0 $
and
$\alpha $
large enough. To estimate the first factor in equation (3.71), we apply Lemma 3.15 (note that
)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn155.png?pub-status=live)
On the right side, we use the functional calculus for self-adjoint operators
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn156.png?pub-status=live)
where in the last step we applied Lemma 3.16 and used
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn157.png?pub-status=live)
which is inferred from equation (3.5) using
$\vert \hspace {-1pt}\vert \Delta \varphi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}< \infty $
. Returning to equation (3.71), we have shown that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn158.png?pub-status=live)
and hence we are in a position to apply Corollary 3.6. This implies for all
$\alpha $
large
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn159.png?pub-status=live)
$\underline {\mathrm {Term}\ {\mathcal N}_1}$
. We start by inserting equation (2.22) for
$G_K^1$
in expression (3.61b). Since the Weyl operator commutes with
$u_{\alpha }$
, R and
$P_{\psi } = |\psi \rangle \langle \psi | $
, we can apply equation (3.10a) to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn160.png?pub-status=live)
where we used that
$h_{K,x}$
is real-valued. Note that
$\langle h_{K,\cdot } | \operatorname {Re}( w_{P,y}^1 ) \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} $
is a y-dependent multiplication operator in the electron variable. With
$( T_y e^{A_{P,y}})^{\dagger } = T_{-y} e^{-A_{P,y}}$
and equation (3.67), we can thus write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn161.png?pub-status=live)
where we introduced the operator
$R_{1,y} = R_{1,y}^1 + R_{1,y}^2$
with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn162.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn163.png?pub-status=live)
Using Lemma 3.9 in combination with
$\vert \hspace {-1pt}\vert \nabla P_{\psi } \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} + \vert \hspace {-1pt}\vert \nabla R^{1/2}\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} <\infty $
(see Lemmas 3.7 and 3.8), we can bound the first operator, for any
$\Psi \in \mathscr H$
, by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn164.png?pub-status=live)
To estimate the second operator, we write out the inner product, use Cauchy–Schwarz twice, apply Corollary 3.11 (with
$A=1$
,
$X=R$
and
$Y=P_{\psi }$
) and use equation (3.32a),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn165.png?pub-status=live)
Combining the above estimates we arrive at
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn166.png?pub-status=live)
Since
$\psi (x)$
decays exponentially for large
$|x|$
, the function
$f_{\alpha }(y) := \vert \hspace {-1pt}\vert u_{\alpha } T_{-y} P_{\psi } \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} $
satisfies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn167.png?pub-status=live)
With this at hand, we can estimate the part containing the tail. Invoking Corollary 3.14
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn168.png?pub-status=live)
To estimate the first term in equation (3.81), we proceed similarly as in the bound for
$\mathcal N_{02}$
. We insert the identity (3.70), apply Cauchy–Schwarz and employ equation (3.74). This leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn169.png?pub-status=live)
In the remaining norm, we use the fact that
$R_{1,y}$
changes the number of phonons at most by one, and thus we can apply Lemma 3.15 and equation (3.85), together with equation (3.55a), to get
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn170.png?pub-status=live)
With Corollary 3.6, equation (3.86) and
$\vert \hspace {-1pt}\vert f_{\alpha } \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^{\infty }\)}}} \leq 1 $
, this leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn171.png?pub-status=live)
$\underline {\mathrm {Term}\ {\mathcal N}_{2}}$
. The strategy for estimating this term is similar to the one for
$\mathcal N_1$
. Proceeding as described before equation (3.81), one obtains
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn172.png?pub-status=live)
with
$R_{2,y} = R_{2,y}^1 + R_{2,y}^2$
and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn173.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn174.png?pub-status=live)
It follows in close analogy as for
$R_{1,y}$
in equations (3.82a) and (3.82b) that given any
$\Psi \in \mathscr H$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn175.png?pub-status=live)
and since
, we can use Corollary 3.14 to estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn176.png?pub-status=live)
To bound the first term in equation (3.91), we proceed similarly as for
$\mathcal N_{01}$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn177.png?pub-status=live)
The last integral is estimated again via Corollary 3.6, and thus
$|\mathcal N_{2}^< | \le C \alpha ^{-5+3\delta }$
.
Collecting all relevant estimates and choosing
$\delta>0$
small enough completes the proof of the proposition.
3.6 Energy contribution
${\mathcal E}$
In this section, we prove the following estimate for the energy contribution
${\mathcal E}$
defined in equation (3.8a).
Proposition 3.18. Let
${\mathbb N}_1 = \mathrm {d}\Gamma (\Pi _1)$
and choose
$c>0$
. For every
$\varepsilon>0$
, there is a constant
$C_{\varepsilon }>0$
(we omit the dependence on c) such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn178.png?pub-status=live)
for all
$|P|/ \alpha \le c $
and
$\alpha $
large enough.
Proof. Since
$G_K^0 = \psi \otimes \Upsilon _K$
,
$ h^{\mathrm{Pek}} \psi = 0$
and
${\mathbb N}\Upsilon _K = {\mathbb N}_1 \Upsilon _K $
, one has
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn179.png?pub-status=live)
where both terms provide contributions to the energy of order
$\alpha ^{-2}$
.
$\underline {\mathrm {Term}\ {\mathcal E}_1}$
. Recall that
$H(y) = \langle \psi | T_y \psi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} $
, and use this to write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn180.png?pub-status=live)
With equations (3.63), (3.3) and (3.33), it follows that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn181.png?pub-status=live)
and since
$e^{a^{\dagger }( \alpha w_{P,y}^0) }$
commutes with
$\mathbb U_K {\mathbb N}_1 \mathbb U_K^{\dagger }$
and
$e^{a( \alpha w_{P,y}^0) } \Upsilon _K = \Upsilon _K $
(we use
$\mathbb U_K a^{\dagger }(f^0) \mathbb U_K^{\dagger } = a^{\dagger } (f^0)$
for
$f^0 \in \text {Ran}(\Pi _0)$
), this leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn182.png?pub-status=live)
Because
$\mathbb U_K {\mathbb N}_1 \mathbb U_K^{\dagger } $
is quadratic in creation and annihilation operators, we can expand the exponential in the inner product and use that only the zeroth- and second-order terms give a nonvanishing contribution,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn183.png?pub-status=live)
Next, we add and subtract the Gaussian to separate the leading-order term,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn184.png?pub-status=live)
In
${\mathcal E}_{111}^{\mathrm{lo}}$
, we use
$| H(y) - 1 | \le C y^2$
and Corollary 3.14 to replace
$H(y)$
by unity at the cost of an error of order
$\alpha ^{-7}$
. In the term where
$H(y)$
is replaced by unity, we perform the Gaussian integral and use Proposition 3.17 and again Corollary 3.14. This leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn185.png?pub-status=live)
The error in equation (3.102) is bounded with the help of Lemma 3.5,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn186.png?pub-status=live)
In
${\mathcal E}_{112}$
, we use the Cauchy–Schwarz inequality, Corollary 3.14 and Lemma 3.4, to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn187.png?pub-status=live)
With
$\vert \hspace {-1pt}\vert |\cdot |^n H \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^1\)}}} \le C_n$
, we can now apply Corollary 3.6 to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn188.png?pub-status=live)
In order to bound
${\mathcal E}_{12}$
in equation (3.98), we decompose
$\Upsilon _K = \Upsilon _K^< + \Upsilon _K^>$
for some
$\delta>0$
(see equation (3.67)) and then follow similar steps as described below equation (3.69). This way we can estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn189.png?pub-status=live)
While the second term is bounded via equation (3.55b) by
$C_{\delta }\, \alpha ^{-12}$
, in the first term we apply Lemma 3.15 and use the functional calculus for self-adjoint operators,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn190.png?pub-status=live)
Since
$P_f$
changes the number of phonons in
${\mathcal F}_1$
at most by one, we can proceed by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn191.png?pub-status=live)
where we used
$1\le \alpha ^{\delta }$
, Lemma 3.16 and equation (3.77) in the second step. We conclude via Corollary 3.6 that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn192.png?pub-status=live)
$\underline {\mathrm {Term}\ {\mathcal E}_2}$
. Here, we start with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn193.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn194.png?pub-status=live)
with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn195.png?pub-status=live)
and
$\xi _P$
defined in equation (2.18). We record the following properties of
$l_y$
and its derivative. The proof of the lemma is postponed until the end of the present section.
Lemma 3.19. For
$k=0,1$
and for all
$n\in \mathbb {N}_0$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn196.png?pub-status=live)
Note that, by Lemma 3.7,
$j_y$
clearly has these properties as well. We proceed by writing
${\mathcal E}_{21} = \mathcal {E}_{21}^0 + \mathcal {E}_{21}^P$
with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn197.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn198.png?pub-status=live)
and estimate the two parts separately. Using the canonical commutation relations and (3.53c), we evaluate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn199.png?pub-status=live)
where we used that
$l_y$
is real-valued. Note that
$l_{-y}(-z)=l_y(z)$
. As discussed in Remark 3.3,
$n_{0,1}(y)$
is even, and using the arguments therein one can conclude that
$ \Theta ^{-2}_K\mathrm {Im}(w_{P,y}^1)$
and
$\mathrm {Im}(w_{P,y}^0)$
are odd functions on
$\mathbb {R}^6$
since
$(y,z)\mapsto \mathrm {Im}(w_{P,y})(z)$
is odd on this space, and hence
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn200.png?pub-status=live)
Thus, with
$\mathrm {Re}(w_{P,y})=w_{0,y}$
, and with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn201.png?pub-status=live)
we finally have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn202.png?pub-status=live)
Note that
$v\in L^1\cap L^{\infty }$
since
$y\mapsto \vert \hspace {-1pt}\vert l_y \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}$
is, while
$\vert \hspace {-1pt}\vert w_{0,y}\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}$
is uniformly bounded in y. Because of
$\varphi (z) = -\langle \psi |h_{\cdot }(z)\psi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}}$
and
$\nabla _z h(x-z)\, =\, -\nabla _x h(x-z)$
, we have by integration by parts
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn203.png?pub-status=live)
Thus,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn204.png?pub-status=live)
Since
$\psi $
is a smooth function with uniformly bounded derivatives, there exists a
$C>0$
such that for all y
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn205.png?pub-status=live)
Moreover, for
$k=0,1$
and every
$z\in {\mathbb R}^3$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn206.png?pub-status=live)
The first statement follows easily from Lemma 3.7; to show the second one, use
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn207.png?pub-status=live)
and apply the Hardy–Littlewood–Sobolev inequality. This, together with equation (3.38), shows that there exists a function f in
$L^2({\mathbb R}^3, \mathrm {d} z)$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn208.png?pub-status=live)
Now, let
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn209.png?pub-status=live)
and note that
$\vert \hspace {-1pt}\vert b_y \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}} \leq \frac {1}{4} y^4\vert \hspace {-1pt}\vert \Delta \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
which is finite since
$\Delta \varphi \in L^2$
. This equation, together with equation (3.125), implies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn210.png?pub-status=live)
From this and from
$v\in L^1\cap L^{\infty }$
, it is also easy to deduce that
$|\cdot |^{-2}v \in L^1\cap L^{\infty }$
. We can thus write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn211.png?pub-status=live)
and use Lemma 3.5 for
$g=|\cdot |^{-2}|v|$
to bound
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn212.png?pub-status=live)
Using equation (3.127), the definition of
$\lambda =\frac {1}{6}\vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}^2$
as well as
$\int y^2 e^{-y^2}\mathrm {d} y= \frac {3}{2}\pi ^{3/2}$
, we further have that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn213.png?pub-status=live)
which finally gives the estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn214.png?pub-status=live)
using Proposition 3.17.
In a similar fashion as for
$\mathcal {E}_{21}^0$
, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn215.png?pub-status=live)
Explicit computation, using
$\Pi _0 = \frac {3}{\vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}}\sum _{i=1}^3 |\partial _i \varphi \rangle \langle \partial _i \varphi |$
and
$\langle \varphi | \nabla \varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = 0$
, gives
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn216.png?pub-status=live)
Note that the real part of the above is odd as a function of y and hence
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn217.png?pub-status=live)
and, taking rotational invariance of
$\varphi $
into account, we arrive at
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn218.png?pub-status=live)
Further, note that
$|\vert \hspace {-1pt}\vert \nabla \varphi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}^2 - \left (\nabla \varphi \ast \nabla \varphi \right )(y)|\leq Cy^2$
and thus, by Lemma 3.7 and Corollary 3.6, one obtains
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn219.png?pub-status=live)
This completes the analysis of
${\mathcal E}_{21}$
.
In order to estimate the term
$\mathcal {E}_{22}$
, we proceed as before by splitting
$\Upsilon _K=\Upsilon _K^{<}+\Upsilon _{K}^{>}$
. Using equation (3.44), we can estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn220.png?pub-status=live)
where we used Corollary 3.14 and Lemmas 3.7 and 3.19. The term involving
$\Upsilon _K^{<}$
is split again into two contributions,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn221.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn222.png?pub-status=live)
To bound the first one, we proceed as in equation (3.107), that is, use Lemma 3.15 and the fact that
$\phi (l_y)$
changes the number of phonons at most by one. This leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn223.png?pub-status=live)
Furthermore, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn224.png?pub-status=live)
where we used
$[iP_f y,\phi (f)]=\pi (y\nabla f)$
and
. Note that in order to estimate the remaining expression, it is not sufficient to directly apply Corollary 3.6. To obtain a better bound, we first replace
$n_{\delta ,\eta }(y)$
by
$e^{-\eta \lambda \alpha ^{2(1-\delta )} y^2}$
and then, for the part containing the Gaussian, we use that
$\vert \hspace {-1pt}\vert l_y\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}$
and
$\vert \hspace {-1pt}\vert \nabla l_y \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}$
provide additional factors of
$|y|$
, as is shown below. More precisely, with Lemma 3.19 and the aid of Lemmas 3.5 and 3.16, we bound
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn225.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn226.png?pub-status=live)
Next, we use that by equation (3.125) there exists an
$L^2$
function f such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn227.png?pub-status=live)
Hence, by integration
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn228.png?pub-status=live)
Regarding the second term in equation (3.140), we need to bound
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn229.png?pub-status=live)
where we proceed in a similar way as above, using that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn230.png?pub-status=live)
In fact, since
$\nabla \varphi (z) =-\langle \psi |h_{\cdot }(z) \nabla \psi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} -\langle \nabla \psi |h_{\cdot }(z) \psi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}}$
, we have the identity
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn231.png?pub-status=live)
Again, using that
$\psi $
has bounded derivatives, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn232.png?pub-status=live)
and the desired inequality now follows from
$| H(y)-1 | \le Cy^2$
and equation (3.123). Given equation (3.114), we can use Lemma 3.5 to replace
$n_{\delta ,\eta }(y)$
in equation (3.145) with
$e^{-\lambda \eta \alpha ^{2(1-\delta }) y^2}$
at the energy penalty
$C \alpha ^{-6 + 9\delta /2}$
, and then use equation (3.146) to bound the remaining integral involving the Gaussian factor, which yields an error of the same order. Altogether, this gives the estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn233.png?pub-status=live)
For the term
$\mathcal {E}^P_{22}$
, we proceed in exactly the same way as in equation (3.139):
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn234.png?pub-status=live)
where the last estimate follows from Corollary 3.6 and the assumption
$|P|\le c \alpha $
.
Combining the relevant estimates, that is, equations (3.103), (3.104), (3.106) and (3.110) for
${\mathcal E}_1$
as well as equations (3.131), (3.136), (3.137), (3.149) and (3.150) for
${\mathcal E}_2$
, we arrive at the statement of Proposition 3.18, thus providing an appropriate bound for
$\mathcal {E}$
.
Proof of Lemma 3.19.
Since H has the desired properties, we need to show them for
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn235.png?pub-status=live)
To this end, we introduce
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn236.png?pub-status=live)
and start with the following observation: Suppose
$f_1,f_2,f_3$
and
$f_4$
are functions in
$\mathcal {S}$
. Then
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn237.png?pub-status=live)
In fact,
$|S(y)| \leq C \|f_4\|_{\scriptscriptstyle{\text{\(L^{\infty }\)}}} \|f_3\|_{\scriptscriptstyle{\text{\(L^{\infty }\)}}} \|f_1\|_{\scriptscriptstyle{\text{\(L^p\)}}} \|f_2\|_{\scriptscriptstyle{\text{\(L^q\)}}}$
for all
$1<p<3/2, q=3p/(5p-3)$
by the Hardy–Littlewood–Sobolev inequality. Since
$\int \mathrm {d} y |y|^n f_3(u+y)\leq 2^{n-1}\left (|u|^n \| f_3 \|_{\scriptscriptstyle{\text{\(L^1\)}}} +\||\cdot |^nf_3\|_{\scriptscriptstyle{\text{\(L^1\)}}}\right )$
, we have also
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn238.png?pub-status=live)
from which (3.153) follows. Moreover,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn239.png?pub-status=live)
Indeed, we have for all
$n\geq 0$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn240.png?pub-status=live)
since m can be chosen arbitrarily large by assumption. Thus, it suffices to prove the desired statement for the functions
$\|\nabla ^k l_y^{(1)}\|_{\scriptscriptstyle{\text{\(L^2\)}}}^2$
. For
$k=0$
, we use equation (3.124) to compute
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn241.png?pub-status=live)
The statement now follows easily from equation (3.153) and Lemma 3.7. Arguing again via equation (3.155), for
$k=1$
it suffices to show the statement for
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn242.png?pub-status=live)
(the first equality follows from
$\nabla _z h_x(z) = - \nabla _x h_x(z)$
and integration by parts). Using (3.124), we find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn243.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn244.png?pub-status=live)
We arrive at the desired conclusion by Lemma 3.7 and equation (3.153).
3.7 Energy contribution
${\mathcal G}$
The energy contribution
${\mathcal G}$
, defined in equation (3.8b), is evaluated by the following proposition.
Proposition 3.20. Let
${\mathbb H}_K$
as in equation (2.4),
${\mathbb N}_1= \mathrm {d} \Gamma (\Pi _1)$
and choose
$c>0$
. For every
$\varepsilon> 0$
, there exists a constant
$C_{\varepsilon }>0$
(we omit the dependence on c) such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn245.png?pub-status=live)
for all
$| P |/\alpha \le c$
and all
$\alpha $
large enough.
Proof. Using
$h^{\mathrm{Pek}} G_K^0 = 0$
and
$\mathbb NG_K^0 = {\mathbb N}_1 G_K^0$
, we can decompose
${\mathcal G}$
into two terms
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn246.png?pub-status=live)
where the first term will contribute to the error while the second one provides an energy contribution of order
$\alpha ^{-2}$
. We proceed for each one separately.
$\underline {\mathrm {Term}\ {\mathcal G}_1}$
. With the aid of equations (3.67) and (3.80) and
$(T_y e^{A_{P,y}} )^{\dagger } = T_{-y} e^{-A_{P,y}}$
, one finds
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn247.png?pub-status=live)
where we introduced the operator
$R_{3,y} = R_{3,y}^1 + R_{3,y}^2$
with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn248.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn249.png?pub-status=live)
Proceeding similarly as for
$R_{1,y}^1$
and
$R_{2,y}^2$
in equations (3.82a) and (3.82b), one further verifies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn250.png?pub-status=live)
Recalling the definition
$f_{\alpha }(y) = \vert \hspace {-1pt}\vert u_{\alpha } T_{-y} P_{\psi } \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}}$
and equation (3.86), we can use Corollary 3.14 to find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn251.png?pub-status=live)
In the first term, we proceed with equation (3.74) and Lemma 3.15 to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn252.png?pub-status=live)
which brings us again into a position to apply Corollary 3.6. Hence,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn253.png?pub-status=live)
$\underline {\mathrm {Term}\ {\mathcal G}_2}$
. Here, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn254.png?pub-status=live)
To separate the leading order contribution in
${\mathcal G}_{21}$
, we insert
$1 = \mathbb U_K^{\dagger } \mathbb U_K$
next to
$G_K^0$
and bring
$\mathbb U_K^{\dagger }$
to the right side of the inner product. With
$\mathbb U_K\Upsilon _K = \Omega $
, equation (3.53c) and equation (3.63) this gives
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn255.png?pub-status=live)
where
$\underline {\, \cdot \, }$
is defined in equation (3.52a). Next, we write
$ W(\alpha \widetilde w_{P,y}) = n_{0,1}(y) e^{a^{\dagger }(\alpha \widetilde w_{P,y})} e^{-a(\alpha \widetilde w_{P,y})}$
and move the first exponential to the left side and the second exponential to the right side until they act both on the Fock space vacuum. Using
$e^{-a(f)} a^{\dagger }(g) e^{a(f)} = a^{\dagger }(g) - \langle f |g \rangle $
, we find this way
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn256.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn257.png?pub-status=live)
In the first line, we write
$h_{\cdot } + \varphi _P = h_{\cdot }^0 + h_{\cdot }^1 + \varphi + i \xi _P$
, with
$h_{\cdot }^i = ( \Pi _ih)_{\cdot }$
, and use that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn258.png?pub-status=live)
since
$h_{x}^0 + i \xi _P \in \text {Ran}(\Pi _0)$
whereas
$h_{K,x}^1 \in \text {Ran}(\Pi _1)$
. Finally, we can replace a and
$a^{\dagger }$
by
$\phi $
, and then transform back with
$\mathbb U_K$
, using equation (3.53c), in order to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn259.png?pub-status=live)
To summarize, we have shown that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn260.png?pub-status=live)
with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn261.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn262.png?pub-status=live)
In the first term, we add and subtract the Gaussian,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn263.png?pub-status=live)
and proceed with
${\mathcal G}_{211}^{\mathrm{lo}}$
by inserting
$h_{\cdot }^1 = h_{K,\cdot }^1 +( h_{\cdot }^1 - h_{K,\cdot }^1)$
,
$T_y = 1 + ( T_y-1 ) $
and
$u_{\alpha } = 1 + ( u_{\alpha } -1 ) $
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn264.png?pub-status=live)
Since
$P_{\psi } \phi (\varphi ) R = 0$
, we have
${\mathcal G}_{211}^{{\mathrm{lo}},1} = \frac {2}{\alpha ^2} \big \langle \Upsilon _K | ( {\mathbb H}_K -{\mathbb N}_1 ) \Upsilon _K \big \rangle _{\scriptscriptstyle{\text{\(\mathcal F\)}}} (\frac {\pi }{\lambda \alpha ^2})^{3/2}$
, cf. (2.4), and hence we can use Proposition 3.17 to conclude that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn265.png?pub-status=live)
For the other terms, we shall show the combined error estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn266.png?pub-status=live)
In the last term, we recall
$h_{\cdot }(y)= h_{K=\infty ,\cdot }(y)$
, and apply Lemma 3.9 in combination with
$\vert \hspace {-1pt}\vert R^{1/2} u_{\alpha } T_{-y} \nabla \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} \le C$
. This gives
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn267.png?pub-status=live)
Next, we write
$T_y - 1 = \int _0^1 \mathrm {d} s T_{sy} (y \nabla )$
in the third term to obtain an additional
$|y|$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn268.png?pub-status=live)
where the factor
$\sqrt K$
comes from the
$L^2$
norm of
$h_{K,0}^1$
in the bound on the first field operator (since
$\Delta R^{1/2}$
is unbounded, we cannot apply the commutator method to this part). In the second term, we use
$\psi (x)\le C e^{-|x|/C}$
for some
$C>0$
, and thus
$\vert \hspace {-1pt}\vert (u_{\alpha } -1) \psi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}} \le C e^{-\alpha /C}$
, to estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn269.png?pub-status=live)
This proves equation (3.178).
To bound the remaining contributions in
${\mathcal G}_{211}^{\mathrm{err}}$
and
${\mathcal G}_{212}$
, we shall use
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn270.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn271.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn272.png?pub-status=live)
Using the exponential decay of
$\psi $
and
, for
$k=0,1$
, it is easy to show that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn273.png?pub-status=live)
To verify equations (3.182a) and (3.182b), use
$u_{\alpha } T_{-y}\phi (h_{\cdot })= \phi (h_{\cdot - y } ) u_{\alpha } T_{-y}$
and Cauchy–Schwarz to bound
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn274.png?pub-status=live)
Now, we can use equation (3.44) and Lemma 3.9 to obtain equation (3.182a). To estimate
$\ell _2(y)$
, defined in equation (3.174b), we proceed with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn275.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn276.png?pub-status=live)
and considering the first line, we use Cauchy–Schwarz, write out the two inner products (in the phonon variable) and then use Cauchy–Schwarz again,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn277.png?pub-status=live)
where the last step follows from Lemma 3.4 and Corollary 3.11 together with
$\underline { h_{K,\cdot } } = h_{K,\cdot }^0 + \Theta _K^{-1} h_{K,\cdot }^1$
. Since the second line is estimated similarly, we arrive at equation (3.182b). With equation (3.182a) at hand, we can apply Lemma 3.5 and equation (3.184) to get
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn278.png?pub-status=live)
and further, using equation (3.182b) and Corollary 3.6, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn279.png?pub-status=live)
This completes the analysis of
${\mathcal G}_{21}$
.
Next, we introduce
$ R_{4,y} = R_{4,y}^1 + R_{4,y}^2 $
with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn280.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn281.png?pub-status=live)
Inserting equations (3.67) and (3.80) into equation (3.168) it follows that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn282.png?pub-status=live)
With the aid of Lemma 3.9, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn283.png?pub-status=live)
and proceeding similarly as in equation (3.84), we find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn284.png?pub-status=live)
For
$\Psi = \psi \otimes \Upsilon _K^>$
, a second application of Lemma 3.9 (after using unitarity of
$e^{-A_{P,y}}$
) together with
$\vert \hspace {-1pt}\vert \varphi _P \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}} \le C$
for
$|P|/ \alpha \le c$
and Corollary 3.14 is sufficient to find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn285.png?pub-status=live)
with
$f_{2,\alpha }$
defined in equation (3.183). Using this bound in
$G_{22}^>$
and recalling Corollary 3.14 and equation (3.184) we thus obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn286.png?pub-status=live)
In
${\mathcal G}_{22}^<$
, we proceed by inserting equation (3.70) and use equation (3.74) and Lemma 3.15. This gives
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn287.png?pub-status=live)
The derivation of a suitable bound for the norm in the integrand is more cumbersome, so we go through it step by step. To shorten the notation let
$G_K^{0<} = \psi \otimes \Upsilon _K^<$
. We start from equations (3.192) and (3.193) where we insert
$h_{\cdot } = h_{K,\cdot } + ( h_{\cdot } - h_{K, \cdot } ) $
and use the triangle inequality,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn288.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn289.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn290.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn291.png?pub-status=live)
For the second and fourth line, we apply Lemma 3.9 a second time (after bringing
$({\mathbb N}+1)^{1/2}$
to the right of a and
$a^{\dagger }$
) to find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn292.png?pub-status=live)
In the first and third line, we use the functional calculus and write out
$A_{P,y}= i P_f y + i g_{P}(y)$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn293.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn294.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn295.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn296.png?pub-status=live)
Now, we use
$[iP_f y,\phi (f)]=\pi (y\nabla f)$
such that we can estimate the first line by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn297.png?pub-status=live)
To bound the first line, we use again Lemma 3.9, while in the second line we use
$(\nabla h_K)_{\cdot } = - \nabla (h_{K,\cdot }) = - [\nabla , h_{K,\cdot }]$
and (3.44) together with
$\vert \hspace {-1pt}\vert \nabla \varphi _P \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}} \le C $
for
$|P|/ \alpha \le c$
. Together, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn298.png?pub-status=live)
where the factor
$\sqrt K$
in the first step comes from the
$L^2$
-norm of
$h_{K,0}$
, and the last step follows from Lemma 3.16. In a similar fashion, one shows
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn299.png?pub-status=live)
and, with equation (3.77), one also verifies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn300.png?pub-status=live)
Collecting the estimates (3.198), (3.201), (3.202) and (3.203), we arrive at
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn301.png?pub-status=live)
Now, we can apply Corollary 3.6 together with equation (3.184) to bound the right side of equation (3.196). The result is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn302.png?pub-status=live)
In view of the estimates (3.165), (3.167), (3.177), (3.178), (3.188), (3.189), (3.195) and (3.205), the proof of Proposition 3.20 is now complete.
3.8 Energy contribution
${\mathcal K}$
Recall that
${\mathcal K}$
was defined in (3.8c).
Proposition 3.21. Let
${\mathbb H}_K$
as in equation (2.4),
${\mathbb N}_1 = \mathrm {d} \Gamma (\Pi _1)$
and choose
$c>0$
. For every
$\varepsilon> 0$
, there exists a constant
$C_{\varepsilon }>0$
(we omit the dependence on c) such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn303.png?pub-status=live)
for all
$|P|/\alpha \le c $
and all
$\alpha $
large enough.
Proof. We split this contribution into three terms
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn304.png?pub-status=live)
and note that
${\mathcal K}_1$
provides the energy contribution of order
$\alpha ^{-2}$
.
$\underline {\mathrm {Term}\ {\mathcal K}_1}$
. We start again by writing
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn305.png?pub-status=live)
and proceed for the first term similarly as in the computation of
${\mathcal G}_2$
; see equation (3.168). This leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn306.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn307.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn308.png?pub-status=live)
We go on with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn309.png?pub-status=live)
and in the leading-order term, we insert
$T_y= 1 + (T_y-1)$
and
$u_{\alpha } = 1 + (u_{\alpha } -1 )$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn310.png?pub-status=live)
Since
$R h^{\mathrm{Pek}} R = R$
, one finds
$ {\mathcal K}_{111}^{\mathrm{lo,1}} = - \frac {1}{\alpha ^2} \big \langle \Upsilon _K | ( {\mathbb H}_K - {\mathbb N}_1 ) \Upsilon _K \big \rangle _{\scriptscriptstyle{\text{\(\mathcal F\)}}} (\frac {\pi }{\lambda \alpha ^2})^3$
, cf. equation (2.4), and with the aid of Proposition 3.17, this gives the leading-order contribution
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn311.png?pub-status=live)
For the other terms, we shall show that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn312.png?pub-status=live)
In the second term, we use
$h^{\mathrm{Pek}}R = Q_{\psi } = 1- P_{\psi }$
to write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn313.png?pub-status=live)
which is exponentially small in
$\alpha $
, since
$\vert \hspace {-1pt}\vert (u_{\alpha }-1) \psi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}} \le C e^{-\alpha / C}$
, and thus with Lemma 3.9 one obtains
$|{\mathcal K}_{111}^{\mathrm{lo,2}}| \le C \sqrt K e^{-\alpha / C }$
. In the next term, we use
$ [h^{\mathrm{Pek}}, u_{\alpha }-1] = - [\Delta , u_{\alpha }] $
and again
$h^{\mathrm{Pek}} R = 1-P_{\psi }$
to get
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn314.png?pub-status=live)
Here, the first line is bounded again exponentially in
$\alpha $
, whereas in the second line we use
$[\Delta ,u_{\alpha }] = 2 (\nabla u_{\alpha }) \nabla + (\Delta u_{\alpha })$
and
$\vert \hspace {-1pt}\vert \nabla u_{\alpha } \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^{\infty }\)}}} + \vert \hspace {-1pt}\vert \Delta u_{\alpha } \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^{\infty }\)}}} \le C\alpha ^{-1}$
, see (2.20). Together with Lemmas 3.8 and 3.9, this implies
$ |{\mathcal K}_{111}^{\mathrm{lo,3}}| \le C \alpha ^{- 6 }$
. In the last term, we employ
$T_y-1 = \int _0^1 \mathrm {d} s T_{sy} (y\nabla )$
,
$[h^{\mathrm{Pek}}, u_{\alpha }] = - [\Delta , u_{\alpha }]$
and
$h^{\mathrm{Pek}}R = Q_{\psi }$
to find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn315.png?pub-status=live)
In both lines, there is an additional factor y, and together with equation (2.20), we thus obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn316.png?pub-status=live)
This proves (3.214).
To estimate
${\mathcal K}_{112}$
and
${\mathcal K}_{111}^{\mathrm{err}}$
, we make use of
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn317.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn318.png?pub-status=live)
where
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn319.png?pub-status=live)
Recalling that by definition
for
$k=0,1$
, it follows that
and thus
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn320.png?pub-status=live)
In order to verify equation (3.219a), use
$h^{\mathrm{Pek}} = -\Delta + V^{\varphi } - \lambda ^{\mathrm{Pek}}$
to write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn321.png?pub-status=live)
Since
$\vert \hspace {-1pt}\vert V^{\varphi } R^{1/2} \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} \le C (\vert \hspace {-1pt}\vert R\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} + \vert \hspace {-1pt}\vert \nabla R^{1/2} \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}}) \le C$
(see Lemma 3.8), it thus follows that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn322.png?pub-status=live)
With this at hand, one applies Lemma 3.9 to conclude the bound stated in equation (3.219a). For
$\ell _3(y)$
, we proceed similarly as in equation (3.187), that is,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn323.png?pub-status=live)
Now, we can apply Lemma 3.5 and (3.221) to estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn324.png?pub-status=live)
and further invoke Corollary 3.6 to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn325.png?pub-status=live)
Next, we come to
${\mathcal K}_{12}$
which we rewrite with the aid of equations (3.67) and (3.80) as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn326.png?pub-status=live)
with the operator
$ R_{5,y} = R_{5,y}^1 + R_{5,y}^2$
and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn327.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn328.png?pub-status=live)
Utilizing Lemma 3.9 and (3.32a), we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn329.png?pub-status=live)
and following the same steps as in equation (3.84),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn330.png?pub-status=live)
After using unitarity of
$e^{-A_{P,y}}$
and equation (3.221), we can apply Lemma 3.9 another time to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn331.png?pub-status=live)
Thus, we can estimate the tail with the aid of Corollary 3.14 and equation (3.221),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn332.png?pub-status=live)
Then we use equation (3.63), equaiton (3.74) and apply Lemma 3.15 to get
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn333.png?pub-status=live)
To bound the norm in the integral, we proceed in close analogy to the steps following equation (3.196). We abbreviate again
$G_K^{0<} = \psi \otimes \Upsilon _K^<$
and start from equations (3.229) and (3.230). With equation (3.221), the functional calculus and
$A_{P,y} = iP_fy + i g_{P}(y)$
, one finds
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn334.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn335.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn336.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn337.png?pub-status=live)
In the second and fourth line, we use
$|g_{P}(y)| \le C\alpha |y|^3$
and Lemma 3.9,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn338.png?pub-status=live)
In the first and third line, we employ the commutator
$[iP_f y,\phi (f)]=\pi (y\nabla f)$
to get
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn339.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn340.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn341.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn342.png?pub-status=live)
After another application of Lemma 3.9, we can use equation (3.67) and then Lemma 3.16 for the terms involving
$P_f$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn343.png?pub-status=live)
while in the other two lines, we use
$(\nabla h_{K})_{\cdot } = - [\nabla ,h_{K,\cdot }]$
, to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn344.png?pub-status=live)
Collecting all estimates we have thus shown that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn345.png?pub-status=live)
Using this bound in equation (3.233), we can invoke Corollary 3.6 together with equation (3.221) in order to obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn346.png?pub-status=live)
$\underline {\mathrm {Term}\ {\mathcal K}_2}$
. Using equations (3.67) and (3.80), one finds
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn347.png?pub-status=live)
with the operator
$R_{6,y} = R_{6,y}^1 + R_{6,y}^2$
and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn348.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn349.png?pub-status=live)
With Lemma 3.9 and equation (3.32a) it is not difficult to verify
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn350.png?pub-status=live)
and since
, we can use Corollary 3.14 to estimate the part with the tail by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn351.png?pub-status=live)
To treat
${\mathcal K}_{2}^<$
we proceed as in (3.233), that is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn352.png?pub-status=live)
It now follows from Corollary 3.6 that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn353.png?pub-status=live)
$\underline {\mathrm {Term}\ {\mathcal K}_3}$
. This term is similarly estimated as the previous one. With the aid of equations (3.67) and (3.80), we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn354.png?pub-status=live)
with the operator
$R_{7,y} = R_{7,y}^1 + R_{7,y}^2$
and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn355.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn356.png?pub-status=live)
Utilizing again Lemma 3.9 and equation (3.32a), one shows that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn357.png?pub-status=live)
with
$f_{3,\alpha }$
defined in (3.220). Invoking Corollary 3.14 and equation (3.221), we thus find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn358.png?pub-status=live)
Similarly, as in equation (3.233), we also obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn359.png?pub-status=live)
By Corollary 3.6 and equation (3.220) it follows that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn360.png?pub-status=live)
This completes the analysis of
${\mathcal K}$
. The proof of Proposition 3.21 follows from combining equations (3.213), (3.214), (3.225), (3.226), (3.232), (3.240), (3.244), (3.246), (3.250) and (3.252).
3.9 Concluding the proof of Proposition 2.8
Combining Propositions 3.18, 3.20 and 3.21, we arrive at
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn361.png?pub-status=live)
Now, for
$K\le \widetilde c \alpha $
we know from Proposition 3.17 that
$\mathcal N \ge C \alpha ^{3}$
for some
$C>0$
, such that the right side is bounded by
$C_{\varepsilon } \, \alpha ^{\varepsilon } r(K,\alpha )$
. It remains to show that one can replace
$\alpha ^{-2}\inf \sigma ({\mathbb H}_K)$
by
$\alpha ^{-2} \inf \sigma ({\mathbb H}_{\infty })$
at the cost of an additional error. To this end, recall that
$\inf \sigma ({\mathbb H}_K) = \langle \Upsilon _K | {\mathbb H}_K \Upsilon _K \rangle _{\scriptscriptstyle{\text{\(\mathcal F\)}}}$
and use the variational principle to find
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn362.png?pub-status=live)
Writing
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn363.png?pub-status=live)
and using Lemma 3.9, we can infer that for any
$\Psi \in {\mathcal F}$
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn364.png?pub-status=live)
By Corollary 3.14, we know that
$\big \langle \Upsilon _K | ({\mathbb N}_1 + 1 ) \Upsilon _K \big \rangle _{\scriptscriptstyle{\text{\(\mathcal F\)}}} \le C$
, and thus
$|\inf \sigma ({\mathbb H}_K) - \inf \sigma ({\mathbb H}_{\infty }) | \le CK^{-1/2}$
. In view of equation (3.253) and Lemma 3.1 this completes the proof of Proposition 2.8.
4 Remaining Proofs
Proof of Lemma 1.1.
The form of the kernel is readily found using second order perturbation theory (we omit the details). (i) The lower bound
$H^{\mathrm{Pek}} \ge 0 $
follows from (1.19) whereas
$H^{\mathrm{Pek}} \le 1 $
is a consequence of
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn365.png?pub-status=live)
(ii) That
$\text {Span} \{ \partial _i \varphi :i=1,2,3\} \subseteq \text {Ker}H^{\mathrm{Pek}}$
follows from translation invariance of the energy functional
${\mathcal F}$
(1.15). To show equality, we argue that there is a
$\tau>0$
such that
$\langle v | H^{\mathrm{Pek}} v \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} \ge \tau \vert \hspace {-1pt}\vert v \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
for all
$v\in L^2({\mathbb R}^3)$
with
$\langle v | \nabla \varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} =0$
(note that this also implies (iii)). Since
$H^{\mathrm{Pek}}$
has real-valued kernel, it is sufficient to consider
$v\in L^2_{{\mathbb R}}({\mathbb R}^3)$
. We start by quoting [Reference Feliciangeli, Rademacher and Seiringer13, Lemma 2.7] stating that there exists a constant
$\tau>0$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn366.png?pub-status=live)
for all
$v \in L^2({\mathbb R}^3)$
. (A key ingredient in the proof of this quadratic lower bound are the results about the Hessian of the Pekar energy functional (1.12) that were obtained in [Reference Lenzmann33]; see [Reference Feliciangeli, Rademacher and Seiringer13] for a detailed derivation). Combined with equation (1.19), and using that for small
$\varepsilon $
the infimum over
$y\in {\mathbb R}^3$
exists, this implies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn367.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn368.png?pub-status=live)
Given any v satisfying
$\langle v | \nabla \varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = 0$
, we choose
$y^* ( \varepsilon )$
such that
$f_v (y^*(\varepsilon ), \varepsilon )$
is minimal. Furthermore, note that for every zero sequence
$(\varepsilon _n)_{n\in {\mathbb N}}$
such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn369.png?pub-status=live)
it follows that
$\liminf _{n\to \infty } f_v (y^*(\varepsilon _n) , \varepsilon _n ) = \infty $
, and hence, we can conclude that
$|y^* ( \varepsilon )| \to 0$
as
$\varepsilon \to 0$
. To proceed, let
$\eta (\varepsilon ) := \varphi - \varphi (\cdot - y^*(\varepsilon ) )$
and assume
$\vert y^*(\varepsilon ) \vert> 0$
(for if
$ y^*(\varepsilon ) = 0 $
it follows directly that
$ f_v (y^*(\varepsilon ) , \varepsilon ) = \vert \hspace {-1pt}\vert v \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}} $
). With this, we can estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn370.png?pub-status=live)
To bound the right side, write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn371.png?pub-status=live)
and use, by dominated convergence, that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn372.png?pub-status=live)
Combining the last statement with
$\vert y^*(\varepsilon ) \vert \to 0$
as
$\varepsilon \to 0$
and
$\langle v | \nabla \varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = 0 ,$
we conclude that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn373.png?pub-status=live)
This completes the proof of items (ii) and (iii). Property (iv) follows from
$H^{\mathrm{Pek}} \le (H^{\mathrm{Pek}})^{1/2}$
and
$\text {Tr}_{L^2}(1-H^{\mathrm{Pek}}) <\infty $
; see Lemma 2.3 for
$K=\infty $
.
Proof of Lemma 2.3.
(i) The bound
$H_K^{\mathrm{Pek}} \restriction \text {Ran} (\Pi _1) \le 1$
follows analogously to equation (4.1) and
$H_K^{\mathrm{Pek}} \restriction \text {Ran} (\Pi _0) = 0 $
holds by definition. The lower bound on
$\text {Ran} (\Pi _1)$
is a consequence of
$( H^{\mathrm{Pek}} - \tau ) \restriction \text {Ran} (\Pi _1) \ge 0$
for some
$\tau> 0$
; see Lemma 1.1, in combination with
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn374.png?pub-status=live)
To verify the latter, let
$v\in \text {Ran}(\Pi _1)$
,
$\Pi _v =|v \rangle \langle v|$
and write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn375.png?pub-status=live)
With Cauchy–Schwarz, it follows that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn376.png?pub-status=live)
and from Corollary 3.11, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn377.png?pub-status=live)
(ii) On
$\text {Ran}(\Pi _0)$
the inequality holds trivially, whereas on
$\text {Ran}(\Pi _1)$
, it follows from
$\Theta _K\le 1$
,
$B_K^2 \le \frac {1}{4}( \Theta _K^{-2} - 1 ) $
,
$\Theta _K^{-2} = ( 1 - (1-H_K^{\mathrm{Pek}}) )^{- 1/2}$
and the elementary inequality
$( 1 - x )^{-1/2} \le 1 + \beta ^{-3/2} x$
for all
$x\in (0,1- \beta )$
.
(iii) Here, we use
$\text {Tr}_{\mathrm {Ran}(\Pi _0)}(1-H_K^{\mathrm{Pek}}) = 3 $
, write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn378.png?pub-status=live)
and apply Corollary 3.11.
(iv) Since
$1-H_K^{\mathrm{Pek}} = \Pi _0 + \Pi _1 (1-H_K^{\mathrm{Pek}}) \Pi _1 = \Pi _0 + 4 T_K$
, cf. equations (2.7a) and (2.7b), we can write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn379.png?pub-status=live)
Using the explicit form of
$\Pi _0$
, one shows that the first term is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn380.png?pub-status=live)
which is finite since
$\Delta \varphi \in L^2$
. For the second term, it follows from a short computation that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn381.png?pub-status=live)
Using the Cauchy–Schwarz inequality and
$\vert \hspace {-1pt}\vert \nabla \psi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}} + \vert \hspace {-1pt}\vert R^{1/2} \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}} + \vert \hspace {-1pt}\vert R^{1/2} \nabla \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{op}\)}}}<\infty $
, see Lemmas 3.7 and 3.8, we can estimate the last expression by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn382.png?pub-status=live)
This completes the proof of the lemma.
Proof of Lemma 2.5.
We recall that
$H_K^{\mathrm{Pek}}\restriction \text {Ran}(\Pi _0) = 0$
and
$T_K = \frac {1}{4} (H^{\mathrm{Pek}}_K - \Pi _1) $
, and set
$S_K= \frac {1}{2} (\Pi _1 +H^{\mathrm{Pek}}_K )$
. For
$(u_n)_{n\in {\mathbb N} }$
an orthonormal basis of
$\text {Ran}(\Pi _1)$
, we further set
$a_n = a(u_n)$
and use this to write the Bogoliubov Hamiltonian as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn383.png?pub-status=live)
Applying the transformation (2.11), a straightforward computation leads to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn384.png?pub-status=live)
The statement of the lemma now follows from
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn385.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn386.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn387.png?pub-status=live)
Proof of Lemma 3.4.
To bound
$ \vert \hspace {-1pt}\vert w^1_{P,y} \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
, we expand
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn388.png?pub-status=live)
where we used
$\Pi _1 (y\nabla )\varphi = 0$
. Thus, since
$\Delta \varphi \in L^2$
, we easily arrive at
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn389.png?pub-status=live)
for some constant
$C>0$
, and with
$|P|\le \alpha c$
we obtain the stated estimated. The bound for
$\vert \hspace {-1pt}\vert \widetilde w_{P,y}^1 \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
follows from
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn390.png?pub-status=live)
where we used that
$\Theta _K$
is real-valued and satisfies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn391.png?pub-status=live)
when restricted to
$\text {Ran}(\Pi _1)$
; see Lemma 2.3. To bound
$\vert \hspace {-1pt}\vert w^0_{P,y} \vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}}$
we use
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn392.png?pub-status=live)
since
$\varphi $
,
$\xi _P$
and
$\Pi _0$
are all real-valued. Expanding
$1-e^{-y\nabla }$
as in equation (4.21), it is easy to conclude that
$\vert \hspace {-1pt}\vert \Pi _0(1-e^{-y\nabla })\xi _P\vert \hspace {-1pt}\vert ^2_{\scriptscriptstyle{\text{\(L^2\)}}} \leq CP^2y^2\alpha ^{-4}$
. Using the explicit form of
$\Pi _0$
and
$\langle \nabla \varphi | \varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} = 0$
, we can write
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn393.png?pub-status=live)
Using the Fourier representation and rotation invariance, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn394.png?pub-status=live)
By the elementary inequality
$|\sin z-z|\leq Cz^3$
, the formula
$\vert \hspace {-1pt}\vert (y\nabla )\varphi \vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(L^2\)}}}^2= 2 \lambda y^2$
and the finiteness of
$\|\Delta \varphi \|_{\scriptscriptstyle{\text{\(L^2\)}}},$
we conclude that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn395.png?pub-status=live)
To prove the last bound, we use
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn396.png?pub-status=live)
and hence with equation (4.24),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn397.png?pub-status=live)
The desired bound now follows from equations (4.22) and (4.28).
Proof of Lemma 3.5.
From Lemma 3.4, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn398.png?pub-status=live)
Hence, there is a constant
$\mu>0$
such that for all
$y^2 \le \alpha ^{-1}$
the weight function (3.33) satisfies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn399.png?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn400.png?pub-status=live)
In the remainder, let us abbreviate
$f_n(y) = |y|^n g(y)$
and
$Z(y) = | n_{\delta ,\eta }(y) - e^{-\lambda \eta \alpha ^{2(1-\delta )} y^2} | $
. We then decompose the integral into
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn401.png?pub-status=live)
with
$B_{\alpha } = \{ y \in {\mathbb R}^3 : y^2 \le \alpha ^{-1} \} $
. The bounds (4.32a) and (4.32b) imply that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn402.png?pub-status=live)
and thus by
$|e^z-1| \le z e^z$
for
$z>0$
, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn403.png?pub-status=live)
The last expression is further bounded by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn404.png?pub-status=live)
and since the resulting expression is uniformly bounded in
$\eta \ge \eta _0$
and
$\alpha $
large, we get
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn405.png?pub-status=live)
To bound the second term in equation (4.33), we estimate
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn406.png?pub-status=live)
To see that the first summand is exponentially small as well, we use equation (4.29), equation (4.24) and
$\operatorname {Re}(w_{P,y}^i) = \Pi _i\operatorname {Re}( w_{P,y}) = \Pi _i\operatorname {Re}( w_{0,y}) $
for
$i=0,1$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn407.png?pub-status=live)
and hence
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn408.png?pub-status=live)
Since
$\varphi $
is real-valued, we have
$\langle \varphi |e^{-y\nabla }|\varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}} =\langle \varphi |e^{y\nabla }|\varphi \rangle _{\scriptscriptstyle{\text{\(L^2\)}}}=(\varphi \ast \varphi )(y)$
and thus
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn409.png?pub-status=live)
Recall that, as shown in [Reference Lieb36], the electronic Pekar minimizer
$\psi $
is radial and nonincreasing and hence
$\varphi $
, cf. equation (1.14), is radial and nonincreasing as well, as convolutions of radial nonincreasing functions are themselves radial nonincreasing functions. Consequently,
$q(y)$
is radial and monotone nondecreasing, and thus
$q(y)\ge q(y')$
for all
$y\in B_{\alpha }^c$
,
$y' \in B_{\alpha }$
. On the other hand, by a simple computation, using the regularity of
$\varphi $
, one finds that
$q(y)\geq C_0y^2$
for some
$C_0>0$
and all
$|y|$
small enough, and thus
$q(y)\geq C_0\alpha ^{-1}$
for all
$y\in B_{\alpha }^c$
and
$\alpha $
large. Therefore,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn410.png?pub-status=live)
for some
$d>0$
, which completes the proof of the lemma.
Proof of Lemma 3.16.
Let
$p = -i\nabla $
. By a straightforward computation using equation (2.11), we arrive at
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn411.png?pub-status=live)
for some orthonormal basis
$(u_n)_{n\in {\mathbb N}}$
of
$L^2(\mathbb {R}^3)$
. That
$B_KpB_K$
is trace-class can be seen via
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn412.png?pub-status=live)
where the second step follows from Lemma 2.3, implying
$\vert \hspace {-1pt}\vert B_{K}\vert \hspace {-1pt}\vert _{\scriptscriptstyle{\text{\(\mathrm{HS}\)}}}\le C$
, and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn413.png?pub-status=live)
By rotation invariance
$\mathrm {Tr}_{L^2}(B_KpB_K) = 0$
. The first term in equation (4.43), on the other hand, is seen to be a two-particle wave function
$\Phi _K$
given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn414.png?pub-status=live)
Thus,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20230612135857253-0595:S2050509423000452:S2050509423000452_eqn415.png?pub-status=live)
where we invoked again equation (4.45).
Competing Interest
The authors have no competing interest to declare.
Financial support
This research was supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme grant agreement No. 694227 (R.S.) and the Maria Skłodowska-Curie grant agreement No. 665386 (K.M.).