The Huang–Yang formula for the low-density Fermi gas: upper bound
Abstract
We study the ground state energy of a gas of spin fermions with repulsive short-range interactions. We derive an upper bound that agrees, at low density , with the Huang–Yang conjecture. The latter captures the first three terms in an asymptotic low-density expansion, and in particular the Huang–Yang correction term of order . Our trial state is constructed using an adaptation of the bosonic Bogoliubov theory to the Fermi system, where the correlation structure of fermionic particles is incorporated by quasi-bosonic Bogoliubov transformations. In the latter, it is important to consider a modified zero-energy scattering equation that takes into account the presence of the Fermi sea, in the spirit of the Bethe–Goldstone equation.
Contents
- 1 Introduction and main result
- 2 Main ingredients in the proof
- 3 Preliminary bounds and effective correlation operator
- 4 First quasi-bosonic Bogoliubov transformation
- 5 Second quasi-bosonic Bogoliubov transformation
- 6 Conclusion of Theorem 1.2
- A Useful bounds on scattering solution and cut-off functions
- B Evaluation of
- C Estimate of the integral in (6.5)
- References
1 Introduction and main result
Establishing asymptotic formulas for the ground state energy of dilute gases has been an active area in Mathematical Physics in recent years, partly motivated by the advances in the physics of cold atoms and molecules. In the bosonic case, the validity of the Lee–Huang–Yang formula [33] for the first two terms in a low-density expansion was recently established [1, 19, 20, 47], following earlier work on the first term in [16, 36]. In this work, we are concerned with the analogous question for fermions, for which the Huang–Yang formula [28] gives the first three terms in an asymptotic expansion at low density. We shall establish the validity of this formula, at least as an upper bound.
We consider a system of fermions with spin in a box with periodic boundary conditions. The interaction between particles is described by the periodization of a potential , i.e., , where we assume that is nonnegative, radial and compactly supported. The Hamiltonian describing the system is given by
(1.1) |
acting on , the Hilbert space of antisymmetric square-integrable functions of space-spin variables , with and . Since is spin-independent, if , then leaves invariant the subspace consisting of -body wave functions with exactly particles of spin .
The Hamiltonian is bounded from below and can be defined as a self-adjoint operator by Friedrichs’ method. The ground state energy in the subspace is given by
(1.2) |
and the ground state energy density in the thermodynamic limit is
(1.3) |
It is well-known that the thermodynamic limit exists and is independent of the boundary conditions [41, 42].
By using a pseudopotential method, Huang and Yang [28] predicted that in the low-density limit ,
(1.4) |
where is the -wave scattering length of the interaction potential , defined by with the solution of the zero-scattering equation
(1.5) |
The Huang–Yang conjecture was derived from a consideration of hard-sphere interactions, but it is expected that the formula extends to a large class of repulsive and short-range interaction potentials. The formula (1.4) is of great interest in the theory of dilute quantum gases, since it exhibits a remarkable universality, with the energy depending on the interaction potential only via its scattering length, at least to the order considered. We refer to [38] for an experimental verification.
On the mathematical side, for a large class of repulsive interactions, including hard spheres, the validity of the first two terms of (1.4) was proved in [34]; more precisely, [34, Thm. 1] states that
(1.6) |
The leading order term comes from the kinetic energy of a non-interacting Fermi gas, while the second-order term provides an effective description of the interaction between components of different spins. (In the case of fermions without spin, this second term vanishes; see [31, 32] for the corresponding result in this case.) An analogue of (1.6) for the thermodynamic pressure at positive temperature was derived in [45].
An alternative proof of (1.6) was recently given in [17], by applying elements of bosonic Bogoliubov theory to Fermi systems. The method in [17] was later extended in [21] to derive an optimal error estimate in the upper bound, and will also be of crucial importance in our analysis here. (See also [30] for an alternative approach.) In fact, in the present paper we shall prove an upper bound containing the precise Huang–Yang correction of order . We impose the following conditions on the interaction potential.
Assumption 1.1.
The interaction potential is nonnegative, radial and compactly supported.
In order to state our main result, we shall need the function defined as
(1.7) |
It is a continuous, increasing function, with .
Theorem 1.2 (Upper bound).
As we will see in the proof, the contribution in (1.8) comes from the expression
(1.10) |
with the Fermi momenta of the two spin components. In particular, the roles of and are symmetric. In fact, since for . The explicit form of in (1) has appeared in the physics literature in [29] (see also [12, 39]).
The Huang–Yang correction in (1.4) can be interpreted as the fermionic analogue of the Lee–Huang–Yang correction for the energy of a dilute Bose gas, whose validity was proved in [19, 20] (lower bound) and [47, 1] (upper bound) (see also [25, 24, 18] for extensions to the free energy at positive temperature). In the analysis of the Bose gas, an important role is played by Bogoliubov’s theory [8], where the correlation between particles is captured by unitary transformations containing quadratic expressions in the creation and annihilation operators on Fock space. The solution of the zero-energy scattering equation (1.5) naturally enters in these transformations.
Our proof of Theorem 1.2 is based on an adaptation of the bosonic Bogoliubov theory to Fermi systems. We will interpret suitable pairs of fermions as bosons, and then construct a trial state using the corresponding quasi-bosonic Bogoliubov transformations. This bosonization method goes back to heuristic arguments in [43, 44] in the attempt to explain the random phase approximation proposed in the 1950s [9, 10, 11, 40]. In the context of the high-density Fermi gas, the bosonization method has been recently made rigorous to study the correlation energy in the mean-field regime [3, 4, 5, 13, 14, 15].
For the low-density Fermi gas, the bosonization method was used in [17] to give a new proof of (1.6), and was improved further in [21] to obtain an optimal error estimate in the upper bound. However, to obtain the precise Huang–Yang correction of order , we need to go substantially beyond the existing analysis. We will in fact introduce two quasi-bosonic Bogoliubov transformations to extract the correlation contribution of the particles. The first Bogoliubov transformation is closely related to the construction in [21], and involves the solution to the zero-energy scattering equation in (1.5). In the present paper, this transformation carries an additional momentum cut-off to ensure that all relevant error terms are of order . The main new tool is the second Bogoliubov transformation, which contains the solution of a modified scattering equation taking into account the presence of the Fermi sea at low momentum. The main ingredients of the proof will be explained in the next section.
Our analysis can be extended in a straightforward way to higher spin fermions or, equivalently, more than two species of spinless fermions. For simplicity, we stick to the case of spin here.
Finally, we remark that the method of the present paper can be used to derive a lower bound with the optimal error ; see [22] for details. The matching lower bound for (1.8) remains open, however.
Organization of the paper. In the next section, we shall explain the main structure of the proof. In Sections 2.1 and 2.2 we introduce the fermionic Fock space and the particle-hole transformation, which gives a convenient representation of the correlation Hamiltonian in Proposition 2.3. In Section 2.3 we construct the trial state in terms of two quasi-bosonic unitary transformations and , defined in Definition 2.7. A heuristic discussion of the effect of the unitary transformations on the correlation Hamiltonian, and a sketch of the key estimates required, is then given in Section 2.4. Section 3 contains some preliminary bounds that in particular will allow to reduce the correlation Hamiltonian to a simpler, effective one, as formulated in Proposition 3.6. The main technical analysis needed in the proof of Theorem 1.2 is done in Sections 4 and 5, respectively, where the effect of the two unitary transformations and on the various operators of relevance is investigated. The main results of these two sections are summarized in Propositions 4.1 and 5.1, respectively. Finally, in Section 6 we combine the results of the previous sections, and complete the proof of Theorem 1.2. The appendix contains various auxiliary results needed in the proof.
2 Main ingredients in the proof
In this section we shall explain in detail the construction of the trial state and the main strategy of the proof of Theorem 1.2. We will briefly recall the Fock space formalism and the particle-hole transformation, which will be convenient in order to focus on the excitations around the Fermi sea. Then we will introduce quasi-bosonic transformations, which are our main tool to impose the precise correlations among particles needed to capture the energy to the desired order. Finally we will explain heuristically some key estimates in the computation of the energy, leading to the upper bound in Theorem 1.2.
2.1 The fermionic Fock space
It will be convenient to work in second quantization. The fermionic Fock space is given by
Any is of the form with . The vacuum state will be denoted by .
For any , we denote by and the creation and annihilation operators, respectively. These operators satisfy the canonical commutation relations
for any , where . As a consequence, these operators are bounded, with
(2.1) |
For and , we denote and
(2.2) |
where we introduced the operator-valued distributions with . Furthermore, we denote by , , the number operators counting the number of particles of spin , which can be written as
and the total number operator .
The Hamiltonian in (1.1) can be extended to the operator on Fock space
That is, is the restriction of to . Equivalently, we can write
where are the Fourier coefficients of . Here and in the following, we use the convention
for the Fourier transform in a box . (The Fourier transform in will instead be denoted by in the following.)
Finally, we can identity with the subspace where for (which is left invariant by ), and write the ground state energy in (1.2) as
2.2 The particle-hole transformation
For a non-interacting system, the ground state is given by the free Fermi gas state , which is simply the Slater determinant of the plane waves with momentum inside the Fermi sphere(s), namely
where is the Fermi ball of spin , with . As in [17, 21], we can suppose the Fermi ball(s) to be completely filled, i.e., . A simple computation (see e.g. [17, 21]) shows that the energy per unit volume of the free Fermi gas state , , satisfies
(2.3) |
for small , where the shorthand notation stands for the thermodynamic limit.
We will denote by the one-particle reduced density matrix of , which is an orthogonal projection with integral kernel given by
(2.4) |
Note that is real-valued. We shall also define , with integral kernel
(2.5) |
satisfying . We will write and for simplicity. Furthermore, we define and .
Definition 2.1 (Particle-hole transformation).
Remark 2.2.
As explained in [6], the particle-hole transformation is a particular example of a fermionic Bogoliubov transformation, which allows to focus on the excitations around the free Fermi gas. Equivalently, the transformation can be defined by
(2.8) |
which justifies the name particle-hole transformation.
Note that we define the particle-hole transformation slightly differently from the one used in the previous work [17, 21], in order to preserve translation-invariance. Since we have . With the aid of (2.6) we can compute explicitly. Proceeding as in [17] (using (2.6) and normal-ordering all the terms), we obtain
Proposition 2.3 (Conjugation of by ).
In the following, it will be convenient to further decompose into a part involving interactions of particles with the same spin, and involving interactions of particles of opposite spin, i.e.,
(2.11) |
For our trial state, the main contribution to the energy will come from the effective correlation Hamiltonian
(2.12) |
(see Section 3.2 for more details).
2.3 Trial state and quasi-bosonic Bogoliubov transformations
The trial state we are going to use in order to prove Theorem 1.2 is of the form
(2.13) |
where is the particle hole transformation defined in Definition 2.1, is the vacuum state, and , are certain unitary transformations that will be defined below.
The main idea of our approach is to describe the low energy excitations around the Fermi ball by pairs of fermionic particles that display an approximate bosonic behavior. To do that we introduce the quasi-bosonic annihilation operators
(2.14) |
and their adjoints, the corresponding creation operators, where are the Fourier coefficients of the kernels introduced in (2.5). In particular,
(2.15) |
Remark 2.4.
The unitary transformations and are defined as exponentials of expressions quadratic in the operators. The two transformations correspond to two different regimes of the momentum of the bosonic quasi-particle: for high-momenta (with respect to ) we use a “less refined” expression which helps to extract the leading term of the interaction energy, and at the same time renormalizes the interaction potential, while for low momenta we use a “more refined” expression which allows to extract the correct energy of order . In order to introduce the momentum splitting, we choose two smooth functions , denoted by and , respectively, satisfying
(2.16) |
Here is a parameter which will be optimized over at the end; we shall in fact choose . From we can construct the periodic function using the Fourier coefficients , which satisfies the uniform bound (see Lemma A.5). In the following we will also use the notation .
The operators and will implement the desired correlation structure between particles by using two different expressions both related to the zero-energy scattering equation (1.5).
Definition 2.5 (Scattering solutions).
We define the periodic function by for , , where is the solution of the zero-scattering equation (1.5). Moreover, for with a parameter , we define
(2.17) |
where denotes the scattering length of .
Remark 2.6.
Definition 2.7 (Quasi-bosonic transformations).
For a parameter , we define
(2.18) | ||||
(2.19) |
where and are given in Definition 2.5. We also write and in place of and , respectively.
The cut-off helps to regularize since the function with is uniformly bounded in (see Lemma A.2) as . In the expression for in (2.19), we only use for momenta , , , , where and The parameter in the denominator of is introduced for technical reasons to avoid logarithmic divergences otherwise encountered in integrals. At the end we shall choose to satisfy .
Remark 2.8 (Relation to Bethe–Goldstone equation).
From our analysis, it may seem more natural to consider a trial state of the form , with the unitary given by
(2.20) |
where is chosen to satisfy the equation
(2.21) |
(in the infinity volume limit). Note that (2.21) reduces to (1.5) for . Equation (2.21) can be interpreted as describing the scattering (at zero energy) of two particles at relative momentum , with initial momenta within the Fermi sphere which is fully occupied. From the definition of and , the function appearing in (2.20) is naturally supported on the set where . Note that the effective kinetic energy
describes the difference between the kinetic energy of the two particles with initial momenta and excited momenta , and is positive. Equation (2.21) plays a prominent role in the physics literature. With the help of , it can be rewritten in the form
(2.22) |
which is known as the Bethe–Goldstone equation, see [7, Equation 2.2].
For technical reasons, we find it simpler not to work with , but rather split the unitary into two parts, and . Note that the in the definition of depends only on the momentum ; in this sense is “less refined” but it is sufficient to renormalize the interaction in into a softer one with integral given by . On the other hand, the transformation is “more refined” since the in depends on , but it is also simpler than since the function is explicit.
Note that on the support of , the term is subleading compared to , hence it is natural to replace by the usual scattering solution in . On the other hand, on the support of and for , on the right-hand side of (2.21) can be replaced by to leading order, naturally reducing to in this regime.
Remark 2.9 (Configuration space representation).
In the calculation below it will be convenient to estimate the error terms in configuration space. For this reason, we will often write
(2.23) |
where
(2.24) |
Useful bounds on the function will be given in Lemma A.2.
Writing the unitary in configuration space is a bit more complicated, since depends on all the momenta involved in the definition of the transformation. First of all we notice that in the definition of in (2.19), the coefficients and appearing in (2.14) can be replaced by and with
(2.25) |
since and both and are bounded by (a constant times) , hence for small. Moreover, since both and appearing in the denominator in the definition of in (2.17) are positive, we can write
for , , and . Therefore, introducing the functions
(2.26) |
with and defined in (2.25) and (2.4), respectively, we can write in configuration space as
(2.27) |
A big advantage of choosing the explicit function (2.17) (as opposed to the solution of (2.21)) is to allow for this representation in configuration space.
Remark 2.10.
The transformation is related to the unitary transformation used in [21, Definition 4.3] (see also Remark 4.4 and Remark 4.5 in [21]). In [21, Definition 4.3], contains a cut-off around in configuration space, which is essentially equivalent to a cut-off around in momentum space. In the present paper we choose a weaker localization with a momentum cut-off around in order to obtain better error estimates of order .
We conclude this section by showing that the trial state introduced in (2.13) is admissible, i.e., is a state with particles with spin and particles with spin . To see this, note that is clearly admissible, and commutes with for and , which can be deduced from the definitions in a straightforward way. Hence
2.4 Key estimates
Using the trial state (2.13) and Proposition 2.3, we obtain as an upper bound to the ground state energy
where was defined in (2.12). We shall now explain some key estimates in the computation of . In the first step, the unitary operator introduced in (2.18) is responsible for extracting the leading order of the interaction energy given by . Moreover, via we also renormalize the interaction, effectively replacing by the softer interaction
(2.28) |
At this point the correlation operator to leading order reduces to To extract the next order correction to the correlation energy we then conjugate by the unitary operator in the second step.
Step 1: We conjugate the effective correlation operator with the unitary by writing
recalling the notation introduced in (2.18). To leading order, we shall see that
where denotes the “discrete convolution” for functions defined on momentum space, namely
This explains the choice of the function appearing in : From the zero-energy scattering equation (1.5) we deduce that . The operator in (2.3) can equivalently be written as
Since for small , we conclude that
(2.29) |
where is defined in (2.28) and
(2.30) |
We emphasize that is a renormalized version of , where after taking the 2-body scattering process into account, we have the softer interaction potential instead of . The operator will be important only in the next order of the Duhamel expansion, where it leads via a commutator to a term relevant for the correct Huang–Yang correction of order , as will be demonstrated below. In fact, applying (2.29) we can extend the above Duhamel expansion as
(2.31) |
Computing the commutator in the last integrand and normal-ordering the expression leads to constant term, which is the leading contribution. The integration over then just gives a factor . We shall see that
The first term on the right-hand side naturally combines with from (2.3) to yield the scattering length, since
Therefore, we obtain the upper bound
(2.32) |
The rigorous justification of this statement is the content of Proposition 4.1.
Step 2: To complete the proof, we still have to compute the last term in (2.32). We will see that is negligible, i.e., that (see Proposition 5.3). It therefore remains to conjugate the final effective correlation Hamiltonian under . As above we can write
where . From the definition of the unitary operator (see (2.19) and (2.17)), we find that
(2.33) |
As a consequence, we obtain
Computing the commutator and normal-ordering leads to a constant contribution in the last expression, which is the leading term and given by
(2.34) |
where and are given in (2.15). The rigorous justification of this formula is contained in Proposition 5.5. Combining both transformations we find that the energy per unit volume satisfies
To conclude we can remove the and the momentum cut-off in the third term (see Lemma 6.1) to get, in the thermodynamic limit,
The evaluation of the last integral yields the expression in Theorem 1.2, as will be demonstrated in Appendix B.
3 Preliminary bounds and effective correlation operator
In this section, we collect some useful preliminary bounds that we will use frequently in our analysis. We also identify the effective correlation operator from which we will extract the desired expression for the energy in the upper bound in Theorem 1.2. As usual, we will denote by generic constants in the following, which may have different values at different occurrences.
3.1 Some operator bounds
In this subsection we collect some general operator inequalities. For , define
(3.1) |
where we recall the notation (2.7). Similarly, we define
(3.2) |
Lemma 3.1 (Bounds for , ).
Proof.
The proof is analogous to the one of [21, Lemma 4.8]. We have
(3.4) | |||||
where we used that the conditions and imply that for small , and hence . We multiply and divide the expression on the right hand side of (3.4) by and decompose as a sum of three terms. Writing them in configuration space, we obtain
where
(3.5) |
Using (2.1) we thus obtain
where in last step we used (A.3) to bound the various norms involving , as well as and .
For the proof of the second estimate in (3.3) we can proceed in the same way. The extra factor arises from the additional derivative hitting . ∎
Lemma 3.2.
For any , we have
(3.6) |
Proof.
Next, we turn to the number operator on Fock space, which, after the particle-hole transformation , measures the number of excitations above the Fermi sea. We have the following estimates concerning its behavior under the quasi-bosonic transformations and defined in Section 2.3.
Proposition 3.3 (Estimate for – Part I).
Let and . For any normalized , we have
(3.7) |
Proof.
Before considering the analogous problem for we shall prove the following Lemma.
Lemma 3.4 (Integration over ).
The proof shows that the factor can be replaced by any other number. The choice guarantees that the final error bound is independent of it. In the following we will often use
(3.10) |
which follows from Lemma 3.4 by the Cauchy–Schwarz inequality.
Proof.
We start with proving (3.8). Writing the -norm as a sum in momentum space, we have
Since is supported on (see (2.25)), the contribution to the sum from is bounded by . For the remaining part, we bound the sum by the corresponding integral and obtain
The last factor grows only logarithmically with , hence we can bound it by . This term is negligible compared to the one of order above, and we arrive at (3.8). For (3.9) we have analogously
The sum can be bounded very similarly as above, with the result that
∎
Proposition 3.5 (Estimate for – Part II).
Let for some , and . For any normalized , we have
(3.11) |
Furthermore,
(3.12) |
3.2 Simplified correlation operator
In this section, we collect some bounds already discussed in [17, 21], which allow to reduce the correlation operator defined in (2.3) to a simplified one as far as the energy of our trial state is concerned. We shall see that we can ignore the contributions of , , and .
Estimates for and
Estimate for
We note that the state is such that in the -particle sector of the Fock space unless for . This is a consequence of the fact that both and either create or annihilate particles. Since changes the particle number by , we immediately conclude that .
Estimate for
We shall now argue, similarly as above, that also . This follows from the fact that and are orthogonal, since creates or annihilates particles of the same spin, while in the state the number of particles in each spin component is the same, since both and create particle (and annihilate) pairs of opposite spin simultaneously.
In summary, we obtain the following simplification of the correlation Hamiltonian.
Proposition 3.6 (Effective correlation operator).
4 First quasi-bosonic Bogoliubov transformation
In this section we perform the first step discussed in Section 2.4, i.e., we conjugate the effective correlation operator in Proposition 3.6 with the unitary introduced in (2.18). The main result of this section is the following rigorous version of (2.32). Throughout this section, we shall denote by terms that vanish in the thermodynamic limit (and likewise for those that vanish after dividing by ). We shall not specify their -dependence since we take before considering the low density limit.
Proposition 4.1 (Conjugation by ).
The operator can be interpreted as a renormalized version of , where the effect of the 2-body scattering process have been taken into account, leading to the softer interaction potential instead of . This operator, together with the kinetic energy operator , allows us to extract the full constant term of order after conjugation with . The error term will be estimated in the next section.
To prove Proposition 4.1, we first consider the conjugation of the operators , and with respect to the unitary separately in Sections 4.1–4.3. Then in Section 4.4 we use the zero-energy scattering equation (1.5) to justify (2.29), effectively leading to as a renormalization of . In Section 4.5 we prove propagation bounds for the operators and , which help to control some of the error terms. Finally, in Section 4.6 we collect all the estimates and complete the proof of Proposition 4.1.
In the following we will often use the expression (2.23) for in configuration space, together with several estimates on from Lemma A.2.
4.1 Conjugation of
First, let us consider the conjugation of by . Our goal is to extract a (quasi-bosonic) quadratic term which helps to renormalize the interaction potential in .
Proposition 4.2 (Conjugation of by ).
Remark 4.3.
Remark 4.4 (Comparison with [21]).
Proof of Proposition 4.2.
The equality in (4.3) is a straightforward consequence of
inserting the definitions of and in (2.3) and (2.18), respectively, and computing the commutator. We shall skip the details. To prove (4.5), we compute
For simplicity, we only consider the first term involving in the definition (4.2) of , i.e.,
(4.6) |
the estimate for the other term involving is analogous. For the calculation of the commutator , we need to compute
(4.7) | ||||
Correspondingly, we have , with
(4.8) |
(4.9) |
(4.10) |
(4.11) |
(4.12) |
and
(4.13) |
In the following, we shall estimate all these terms. We shall repeatedly use the bound
(4.14) |
for the fermionic creation and annihilation operators. Moreover, we shall employ the bounds in Lemma 3.1 in Lemma A.2. We will estimate the expectation value of all the terms in a general state . We first consider the error terms coming from the second line on the right hand side in (4.7), i.e., , and . It is convenient to divide them all into two parts, which corresponds to an integration by parts in configuration space: this allows us to get a factor when the derivative hits a , and to use instead of to estimate the error terms when the derivative hits a . We start by estimating the term in (4.8). Writing , we split into two terms, which we denote by and , respectively. In it is convenient to replace by with defined as
(4.15) |
This can be done since and for all the terms in the sum. The same argument applies to and .
To estimate , we rewrite both and in configuration space. We have
For fixed and , we can use (4.14) and to bound
where we used the notation , , as in Lemma 3.1 and the discussion after (2.5). Similarly,
(4.16) |
With the aid of the Cauchy–Schwarz’s inequality, we hence find that
Thus, using also Lemma A.2, we obtain
(4.17) |
The estimate of the term is similar. We have
where . Proceeding as above, we can bound
(4.18) |
Using also (4.16) and that
(4.19) |
we obtain the bound
(4.20) |
We now estimate the term defined in (4.9). Similarly as above, we write and correspondingly we split as a sum of two different terms, and . In order to bound , we rewrite it partially in configuration space. Recalling the definition of in (3.1)), we find
We can bound for all . From the bounds in Lemma 3.1, we find
and hence we obtain
(4.21) |
The estimate of can be done analogously as the one of , using in addition (4.19). We omit the details, and directly write the final estimate as
(4.22) |
Next we consider in (4.10). Again we split the term in two, and , by writing . The term can then be written as
where we used the operators introduced in (3.1) and (3.2). Using the bounds in Lemma 3.1 and that , we get with the aid of the Cauchy–Schwarz inequality
(4.23) |
The analysis of can be done similarly as the one for , using also the inequality in (4.19). The result is
(4.24) |
The next error term we estimate is in (4.11). As above, it is convenient to rewrite it in configuration space. Since the constraints and imply that , we have that in (4.11), and hence
Using Lemma A.2 and (4.14), we can bound
(4.25) | |||||
The estimate for can be done similarly, we omit the details. Also for this term we obtain
(4.26) |
Finally we consider the last term in (4.13). Similarly as above, we have for all the terms in the sum. To estimate it, we shall write
and correspondingly split in two terms, denoted by and . The last one vanishes,
since the sum over is zero by symmetry. The term we can write in configuration space as
(4.27) |
We use and integrate by parts in . The derivate then hits either or , and using (4.14) we can bound
(4.28) |
where we used Lemma A.2 in the last step. ∎
4.2 Conjugation of
Proposition 4.5 (Conjugation of by ).
Let and . Under the same assumptions of Theorem 1.2
(4.29) |
where
(4.30) |
and is such that for any
(4.31) |
Proof.
We start by computing
with
Recall from (2.14) that . Hence we need to compute
Using that , and performing a suitable change of variables, we can combine the first two terms in the commutator above, as well as the third and the forth and the last two. Therefore, we find that
with
(4.32) |
In we have due to the constraints on and, with another change of variables, we find that defined in (4.30). In other words, .
In the following we only consider , the bound on works in exactly the same way. Using again that in (4.32), we can write in configuration space
(4.33) |
where we recall the definition of the operator in (3.1). Using , Lemma 3.1 as well as the Cauchy–Schwarz inequality, we obtain
(4.34) |
This yields (4.31). ∎
4.3 Conjugation of
In the next proposition we conjugate the operator introduced in (2.11) with respect to the unitary operator defined in (2.18).
Proposition 4.6 (Conjugation of by ).
Proof.
As in the previous propositions, we have
(4.36) |
For the computation of we shall use that
and hence we need
(4.37) |
Correspondingly, from the six terms on the right-hand side we obtain , and we shall now estimate these terms. The main term comes from normal-ordering of . As in the previous proofs, in our estimates we often use that , together with the bounds proved in Lemma 3.1, Lemma A.2 and (4.14) with . We start by considering . It is convenient to rewrite the term partially in configuration space, as
(4.38) |
where we used again the operator introduced in (3.1). Using that and the Cauchy–Schwarz inequality, we get
We can further bound
where we used (4.14) and . Using Lemma 3.1 we can also estimate
In combination, we thus obtain
(4.39) |
The term in can be estimated similarly, with the same outcome, and we omit the details.
Next we consider . Similarly as for the other terms, we rewrite it in configuration space as
(4.40) |
By the Cauchy–Schwarz’s inequality, we have
The first factor is bounded by . To bound the second, we proceed as in (4.16), i.e., using (4.14) and , to obtain
Here we have used Lemma A.2 in the last step.
The next term we consider is . Before estimating it, we note that from the constraint and , we get that . Hence . Rewriting in configuration space, we then get
(4.41) |
Hence, using again Lemma A.2 and the fact that commutes with ,
(4.42) | |||||
The term can be estimated in the same way, obtaining the same bound.
Finally we consider the term , from which we extract the leading contribution. We start by normal-ordering the terms, using
(4.43) |
This way, we obtain
(4.44) |
where we have used again that in all the summands. The main contribution is the constant first term. Since
the term in the second line is bounded by . To bound the term in the last line in (4.44), we rewrite it in configuration space as
(4.45) |
Hence also this term is bounded by , using that because of (A.2).
Combining all the bounds the proof of Proposition 4.6 is complete. ∎
4.4 Renormalization of by scattering equation
In this subsection we will combine the “ operators” from Propositions 4.2 and 4.5 with in order to obtain a renormalized interaction. Recall the definitions of and in (4.4) and in (4.30), as well as the ones of and in (2.28) and (2.30). Note also that, as all these operators, in (2.11) is quadratic in the ’s, namely,
Proposition 4.7 (Scattering equation cancellation).
Let . Under the same assumptions of Theorem 1.2,
(4.46) |
where is such that for any normalized
(4.47) |
Proof.
Using , a simple computation yields
with
(4.48) |
In the following we will show that, as a consequence of the zero-energy scattering equation (1.5), is suitably small (for large and small ). We shall write , with
(4.49) |
According to (1.5), , hence takes into account the finite-size effect of the periodic box, and vanishes in the thermodynamic limit. In Lemma A.4, we shall in fact show that the periodic function with Fourier coefficients satisfies and as .
In configuration space we can write
where we used the notation introduced in (3.1). With the aid of Lemma 3.2 and the Cauchy–Schwarz inequality, we can estimate
(4.50) |
This term is thus in the thermodynamic limit.
We now consider the remaining term
(4.51) |
To start, we note that
i.e., we can replace and by and , with defined in (2.25), due to the support properties of and . Using also that
we can write
and hence, by the Cauchy–Schwarz inequality,
The last factor can simply be bounded by . In the first factor, we first normal and write , leading to the bound . Since , this yields
We are left with bounding . Note that . Since is assumed to have compact support, . Using also (A.4), we see that
(4.52) |
and hence (using that is bounded, see Lemma A.1). Altogether, this shows that , and thus completes the proof. ∎
Proposition 4.8 (Propagation estimates for and ).
Let and . Under the same assumptions as in Theorem 1.2
(4.53) | ||||
(4.54) |
where the error term is such that for any
Furthermore, for any normalized , we have
(4.55) |
(4.56) |
Proof of Proposition 4.8.
By Duhamel’s formula,
(4.57) |
We start by computing . The terms have the same structure as in Proposition 4.6, the only difference is that has to be replaced by everywhere. We again write them as
As in Proposition 4.6, putting in normal order, we get the constant term. All the other terms are error terms, which we estimate in expectation in a state . In some error terms it will again be convenient to replace by , with defined in (2.25), which is possible by the compact support of . In the estimates below, we often use the bounds in Lemma 3.1, Lemma A.2 and Lemma A.5 together with (2.1) and (4.14).
For , we have, analogously to (4.38),
(4.58) |
Using and the Cauchy–Schwarz inequality, we can bound
(4.59) |
The term can be estimated in the same way, with the same result.
Next we consider . In analogy with (4.40), it is given by
Using (4.16) and the Cauchy–Schwarz inequality, we obtain the bound
(4.60) |
where we used the bounds mentioned at the beginning of the proof.
The term can be written, similarly to (4.41), as
Hence
(4.61) |
where used that commutes with , as well as Lemmas A.2 and A.5 in the last step. The term can be estimated in the same way, with the same result.
Finally we consider the last term coming from the commutator , which equals (4.44) with replaced by . The term in the first line is the desired constant. For the term in the second line, we use
(4.62) |
which follows from the fact that is bounded, as a consequence of (1.5). Hence the expectation value of this term is bounded by . In the last term corresponding to the third line of (4.44), we observe that the summation over is restricted to , due to the constraints on and in the summands. In other words, we can add for free the cut-off function given by
arriving at
Introducing the function with Fourier coefficients and rewriting the term in configuration space as in (4.45), we obtain the bound . We have using the same bound as in (4.62). Altogether, we have thus shown that
with such that
We now consider . Again all the terms in the commutator are of the same type as the ones in in Proposition 4.6, now with replaced by . They can then be treated exactly in the same way as the corresponding error terms in Proposition 4.6, using that
(4.63) |
as in (4.62). We thus get that
with such that
This proves (4.54).
Finally we discuss the bounds in (4.55) and (4.56). Rewriting in configuration space, we have
where is the function with Fourier coefficients . The bound in (4.55) is then a consequence of as shown in (4.63). Indeed, via the Cauchy–Schwarz inequality, we get
Similarly, we can write in configuration space as
where we again could replace by defined in (2.25). We then find
(4.64) |
This completes the proof of Proposition 4.8. ∎
4.5 Propagation of the estimates for and
In this subsection we obtain propagation estimates for and , which will be useful to estimate some of the error terms.
Proposition 4.9 (Propagation of the estimates for and ).
Let , , with . Under the same assumptions as in Theorem 1.2
Note that since and are nonnegative, the above estimate gives two separate bounds for and . The term on the right-hand side is actually small compared to , as will be proved later in Propositions 5.2 and 5.3.
Proof.
Let . We prove the propagation estimate by Grönwall’s Lemma. By Propositions 4.2 and 4.5, we have
(4.65) |
where the four terms on the right-hand side are defined in (4.4), (4.30), (4.2) and (4.29), respectively. Eq. (4.31) in Proposition 4.5 gives a bound on . For the term in (4.2), we first consider the terms involving in (4.6) and write, as in the proof of Proposition 4.2, . Correspondingly, writing the resulting expression in configuration space and recalling the definition of in (3.1), we find
(4.66) |
where in the first term we again replaced by defined in (4.15), which is possible due the presence of the cutoff . Using Lemma 3.1 together with the bounds in (4.14), we can estimate via the Cauchy–Schwarz inequality
(4.67) |
For the first term on the right-hand side of (4.66), we can proceed as in (4.19) to bound
(4.68) |
For the terms involving in in (4.2) we can proceed in the same way, obtaining the same bound.
We now consider the expectation of in (4.65). From Proposition 4.7, we know that , with bounded as in (4.47), and Eqs. (4.55) and (4.56) from Proposition 4.8 give a bound on . It thus remains to give a bound on . By the Cauchy–Schwarz inequality and (4.14),
where we used and . Therefore,
(4.69) |
Inserting the bounds (4.31), (4.67), (4.68) and (4.69) in (4.65), we find
Using the estimate for the number operator in (3.12) together with the constraints on the parameters , and , we conclude that
The desired result then follows from Grönwall’s Lemma. ∎
4.6 Conclusion of Proposition 4.1
Now we have all the necessary ingredients in order to give the proof of Proposition 4.1. Recall the definition of in (2.12). We start by writing
(4.70) |
From Proposition 4.2 and Remark 4.3,
(4.71) |
where is defined in (4.4), and the last error term can be controlled by combining the bound (4.5) from Proposition 4.2 with Propositions 3.5 and 4.9 as
(4.72) |
From Proposition 4.5 we get
(4.73) |
where is defined in (4.30) and the error term can be controlled, again with Proposition 3.5 and 4.9, as
(4.74) |
The terms and can be combined as in Proposition 4.7 to with and defined in (2.28) and (2.30), and satisfying
(4.75) |
where we used the estimate in (4.47) together with Proposition 3.5. Therefore, combining (4.70)–(4.75), we obtain
(4.76) |
with
(4.77) |
The terms involving in the second line on the right-hand of (4.76) above can be written as
which, according to Proposition 4.6, equals
(4.78) |
with
(4.79) |
for all , where we used again Proposition 3.5 and Proposition 4.9. Next we consider the contribution coming from in the first line on the right-hand side of (4.76). By Proposition 4.8,
with
where we used again the bounds in Propositions 3.5 and 4.9. We can further bound the term on the right hand side above using (4.55). Therefore, combining all the estimates, we proved that
(4.80) |
To complete the proof of Proposition 4.1, we need the following Lemma, evaluating the constant contribution in the first line on the right-hand side of (4.80).
Lemma 4.10 (Constant contribution to ).
Proof.
Using , we can rewrite the left-hand side of (4.81) (multiplied by ) as
(4.82) |
where we used the definition of in (4.51). As in (4.52), we have . At the end of the proof of Proposition 4.7, we also showed that . The zero-energy scattering equation (1.5) implies that is bounded, hence
Using that , the third term in (4.82) is bounded as
(4.83) |
Finally, an argument as at the end of the proof of Proposition 4.7 shows at , hence
with . This completes the proof. ∎
5 Second quasi-bosonic Bogoliubov transformation
In this section we complete the second step discussed in Section 2.4. After the conjugation with , the new effective correlation operator is . By conjugating this operator with , we will extract the full constant of order . The main result of this section is a rigorous version of (2.34), plus estimates on the conjugations of , and , which are needed to control the error terms in (4.1) from Proposition 4.1.
Proposition 5.1 (Conjugation by ).
Let , with . Then
(5.1) |
with
Moreover,
(5.2) | ||||
(5.3) | ||||
(5.4) |
We will study the individual terms , , and in Sections 5.1–5.4, and then complete the proof of Proposition 5.1 in Section 5.5. In the proof, it will be very helpful to use the configuration space representation of in (2.27).
5.1 Conjugation of
In this subsection we investigate , and in particular validate (2.33). Moreover, we also prove a rough estimate for , which is helpful to control the corresponding error term in Proposition 4.1.
Proposition 5.2 (Conjugation of by ).
Let , , with . Under the same assumptions as in Theorem 1.2,
(5.5) |
where is such that for any normalized
(5.6) |
Moreover, for any normalized
(5.7) |
Applying the bounds to and , respectively, we obtain from (3.12) (with )
(5.8) |
Proof.
We have
(5.9) |
An explicit computation, using that for , , ,
by the definition of in (2.17), gives . We have already estimated the expectation value of in the proof of Proposition 3.5. Recalling that and applying this bound, we directly obtain (5.6). To prove (5.7) we use Grönwall’s Lemma together with the bound (4.56) on . ∎
5.2 Conjugation of
In this subsection we conjugate with respect to . We will show that does not contribute to the energy to order .
Proposition 5.3 (Conjugation of by ).
Let , , with . Under the same assumptions as in Theorem 1.2
(5.10) |
Proof.
For the proof we find it convenient to split the operator according to the support of the momenta of each of the . To make this precise, we introduce smooth functions with such that
(5.11) |
We denote we define the periodic function with Fourier coefficients for , and also , . Accordingly, we split
where
and contains all the other possible terms. Specifically, , where contains all the terms that have functions and functions .
Let . Using that , it is easy to see that
(5.12) |
where denotes the part of with one and one . The other contributions to can be bounded via the Cauchy–Schwarz inequality as
(5.13) |
In particular, from (5.12) and (5.13) we find that
(5.14) |
We are thus left with estimating , and we will use Grönwall’s Lemma. We start by computing
The terms we find from the commutator above are of the same type as the ones in Proposition 4.5. Explicitly, with
(5.15) |
(5.16) |
In what follows we will often use that , the bounds in Lemma A.5 and Lemma A.6 and that
(5.17) |
with defined in (2.26).
We start by considering . Since , the definition of implies that , and hence since and . Moreover, also and hence we can replace in (5.15) by with the characteristic function of the set . Proceeding then as in Remark 2.9 we can rewrite in configuration space as
(5.18) |
with
(5.19) |
We shall write the right-hand side of (5.18) as . With the aid of (5.17), Lemma A.6 and the Cauchy–Schwarz inequality, we find
For the integration over , we can use the Cauchy–Schwarz’s inequality to get
The second integral above is bounded in (3.9). The first equals
(5.20) |
In combination, we have thus shown that
(5.21) |
The term can be bounded in the same way.
To conclude the analysis of , we now consider in (5.16). Arguing as above, we can replace by . In configuration space, we then obtain
with as in (5.19). Writing and using (5.17), we can estimate
With
and the Cauchy–Schwarz inequality we then get
The integral is again bounded as in (5.20). Inserting this bound, we obtain
(5.22) |
Combining (5.21) and (5.22), we get
where we used the bound on the number operator from Proposition 3.5 (Eq. (3.12) for ) and the constraints and in the last estimate. Therefore, by Grönwall’s Lemma, for any , we have,
Inserting this estimate in (5.14) and using again the bound for the number operator from Proposition 3.5, we find (5.10). ∎
5.3 Conjugation of
The purpose of this subsection is to show that , with defined in (4.2) the kinetic error term coming from Proposition 4.2 from the conjugation of by .
Proposition 5.4 (Final estimate for ).
Let , with . Under the same assumptions of Theorem 1.2
(5.23) |
Proof.
Again we use Grönwall’s Lemma, and compute
The terms in are of the same type as those in Proposition 4.2, see (4.8)–(4.13). The only difference is that each of them has in place of , which distinguishes from . As in the proof of Proposition 4.2, it is enough to prove the bound in (5.23) with replaced by , defined in 4.6. We shall again write and consider the terms individually. We will estimate the expectation value of all the terms in the state . In the estimates we will often use the bounds in Lemma A.2, Lemma 3.1, Lemma A.5 as well as (4.14) and (5.17). Similarly to (3.14), we shall also need that
(5.24) |
Moreover, we will often use the constraints on to replace and by and , respectively, with defined as in (2.25). Similarly as in Proposition 4.2, it is convenient to split , and into a sum of two terms each.
Proceeding as in the proof of Proposition 4.2, taking into account Remark 2.9, we have for the analogue of (4.8)
(5.25) |
with
and
with defined in (4.15). With the aid of Cauchy–Schwarz inequality, together with bounds as in (4.16) and (4.18), using that and , we obtain for the first term
To bound the terms, we shall use Lemma A.2 and , as well as (3.14). This way, we obtain
The Cauchy–Schwarz inequality and (3.9) for the integral with respect to then yields
(5.26) |
The term can be bounded very similarly, using in addition (4.19). We obtain
(5.27) |
Next we consider the term . Proceeding as above, we find again (5.25) with replaced by everywhere, and
where we used the notation in (3.1) and introduced in addition
(5.28) |
Using and the Cauchy–Schwarz inequality, we obtain
Using the bound from Lemma 3.1 for , we find
Similarly, from
(5.29) |
together with (5.24), we find that
Using (3.10) for the integral with respect to , we obtain
The estimate for is similar. It is given by
Using (4.19) and proceeding as above, we obtain
The term can be estimated in the same way as , obtaining the same final bound. We omit the details. We now consider the term . In contrast to the corresponding one in Proposition 4.2 in (4.11), to estimate this term we rewrite it in normal order using that
(5.30) |
Correspondingly, we have . The first term equals
where we used again the notations introduced in (3.1) and (5.28). Since , the Cauchy–Schwarz inequality yields
The last factor can be bounded with the aid of (5.29) and (5.24). For the first factor, we use
and hence
In combination with Lemma A.2 and (3.10), this yields
(5.31) |
Next we consider , which is given by
and can be bounded as
With (3.14), , , Lemma A.2 and Cauchy–Schwarz we thus obtain
In combination with (3.10) this yields the same bound as in (5.31) after integration over . The same applies to and we omit the details. The last term in is
Since and
for all , we can again obtain a bound of the same kind. Altogether, we have thus shown that
(5.32) |
The term can be estimated very similarly, obtaining the same result.
To conclude, we consider the last error term , which corresponds to (4.13) in Proposition 4.2. We shall again start by writing it in normal order, using
and correspondingly obtain terms, which we write as . The term is a constant, given by
To arrive at this expression, we used that support properties of the various functions imply that . If we replace by , the resulting sum over is zero by symmetry. Hence we can also write
On the support of , i.e., for , one readily checks that for and (and ). Moreover, as already used earlier, as a consequence of (1.5). Hence
(5.33) |
The term is given by
Proceeding as with above, we obtain
The integral over can then be bounded with the aid of (3.8), with the result that
(5.34) |
The same bounds holds for .
The last term to study is , which in configuration space can be written as
Using (4.14), (5.17) and (5.24) together with the Cauchy–Schwarz inequality we can bound
where we also used the bounds from Lemma A.2 and Lemma A.5 in the last step. Again we can bound the integral over using (3.8) to conclude that (5.34) also holds for . Combining all the bounds above, we obtain
(5.35) |
Combining all these bounds with (3.12) (for ) and the second bound in (5.8), and using Grönwall’s Lemma, we obtain (5.23). ∎
5.4 Conjugation of
In the next proposition we finally take into account the conjugation of by . The contribution from , when combined with the main contribution of in Proposition 5.2, will give the correct energy of order .
Proposition 5.5 (Conjugation of by ).
Proof.
Since , we have to compute the commutator . The structure of the various terms is the same as those in the commutator in (4.36) in Proposition 4.6. Following the same strategy, we have
The leading order constant term comes from normal ordering , all the other terms are error terms, which we now estimate individually. In all the error terms, we can replace each by defined in (2.25) as a consequence of the constraints on , which will be convenient for the estimates.
The first term to consider is the analogue of (4.38),
where we used the definition (5.28) and introduced
(5.39) |
By the Cauchy–Schwarz inequality and ,
From the bounds (5.29) and , we find with (3.14) that
(5.40) |
The integral over is bounded in (3.10), resulting in
(5.41) |
The estimate of works in the same way, with the same result.
Next we consider the analogue of (4.40), given by
Using the Cauchy–Schwarz inequality together with the bound (3.14), we can estimate
where we also used Lemma A.5 in the last step. Applying again (3.10) for the integration over shows that (5.41) also holds for .
The term is the analogue of (4.41). Similarly as for the term in the proof of Proposition 5.4, it is convenient to write in normal order. Using (5.30), we can write , and the individual terms can be analyzed in a similar way as the corresponding terms in Proposition 5.4. The first term equals
where we used again the operators introduced in (5.28) and (5.39). Using and the Cauchy–Schwarz inequality, we find
We can further bound and similarly for , with the result that
where we also used again (3.14). Applying now (3.9) for the integration over , we arrive at
The term is given by
Similarly as above, we can bound it with the aid of the Cauchy–Schwarz inequality as
With (3.10) to bound the integral over we therefore obtain
(5.42) |
The estimate for the term , which is given by
can be done in a similar way, obtaining the same result. Finally, the term
can be estimated similarly as above, using also that , with the result that
Again (3.10) then yields the validity of the bound (5.42) also for . Combining the bounds, we conclude that for ,
(5.43) |
The estimate for the term can be done in the same way, obtaining the same result.
Finally we analyze , which is the analogue of (4.44). As already mentioned at the beginning of the proof, we shall put it in normal order in order to extract the constant contribution. Using (4.43), we obtain , where the first is the desired constant,
The second term equals
and is thus bounded by , using . The third term we can write in configuration space as
Proceeding similarly as for above, we can bound it as
For the integration over , we can again use Lemma 3.4 to conclude that
(5.44) |
Combining all the estimates we proved (5.37). ∎
5.5 Conclusion of Proposition 5.1
Proof of Proposition 5.1.
We shall now combine the results of Propositions 5.2–5.5 in order to prove Proposition 5.1. The bounds (5.2), (5.3) and (5.4) were already proved in Proposition 5.2, Proposition 5.3 and Proposition 5.4, respectively. Moreover, combining Proposition 5.2 and Proposition 5.5 we obtain
Inserting the definition of in (2.17), the desired estimate (5.1) follows from (5.8) and (5.38). ∎
6 Conclusion of Theorem 1.2
We shall now complete the proof of Theorem 1.2. As already explained in Section 2.4, the starting point is the bound
using the trial state (2.13) and Proposition 2.3. From Proposition 3.6 we further obtain
Combining now Propositions 4.1 and 5.1 and taking the thermodynamic limit we find that
(6.1) |
with
To obtain (6), we also used (2.3) and that . Recall the conditions on the parameters, , , and . By choosing and , we fulfill these conditions and find that .
To complete the proof of Theorem 1.2, we still have to remove the parameter in (6), which was used to have strictly positive lower bound on the relevant denominator, and the cut-off which restricts the integration to . This is the content of the following lemma.
Lemma 6.1 (Removing the gap and the cut-off).
We have
(6.2) |
Proof.
Recall the notation . We split the proof in two parts. We first prove that
(6.3) |
The difference between the two terms in (6.3) is given by
(6.4) |
Without loss of generality let us assume that . Rescaling the variables respect to , we get
(6.5) |
where the boundedness of the integral follows from Lemma C.1 in Appendix C. To conclude it remains to show that
The integration over runs over only, hence we have . Moreover, using that
(6.6) |
we can write
since the first term in (6.6) integrates to zero by symmetry. On the support of , is bounded from below by , hence
This completes the proof of the lemma. ∎
Appendix A Useful bounds on scattering solution and cut-off functions
In this section we prove bounds on the solution of the zero-energy scattering equation (1.5), and on the cut-off functions we use in our analysis.
A.1 Bounds on scattering solution
Let us start this section by collecting some well-known estimates for the scattering solution .
Lemma A.1 (Properties of ).
For the proof of Lemma A.1 we refer to [37] or [35]. Next we investigate defined in (2.24), a periodized version of regularized by a momentum cut-off. Recall also the definition of in Definition 2.5.
Lemma A.2 (Bounds for ).
Proof.
We start by proving (A.2). Let . By the definition (2.24) and the zero-energy scattering equation (1.5) we have that
where we used that because of Lemma A.1. Assumption 1.1 implies that also and the function has compact support. Hence, for large enough , the -norm of its periodization equals the norm, and in particular,
(A.5) |
Therefore, by Cauchy–Schwarz and using that , we can conclude that . For the -norm, we have
where we used again that .
Remark A.3.
The proof above actually shows that in Definition 2.5 is uniformly bounded, i.e., . However, diverges as .
In the following, we bound the norms of and its derivatives. With the aid of (1.5) and (2.16) we can write
where has Fourier transform . The function is thus the periodization of , and hence
It is a simple exercise, using the smoothness and support properties of , to check that , yielding the first bound in (A.3). The second bound on follows in the same way, using . We omit the details. For the last bound in (A.3), we can directly use that with at most one non-vanishing term in the sum, due to the compact support of . This implies that on . Using that , it then follows that
Hence
where we used again Assumption 1.1, (1.5), Lemma A.1 and the fact that by Lemma A.5 below.
It remains to prove (A.4). We have
The function is square integrable and smooth, uniformly in for bounded . Approximating the integral by a Riemann sum, it is then a simple exercise to show that the right-hand side is for large , uniformly in for bounded . We leave the details to the reader. ∎
Lemma A.4 (Estimates of the norms of ).
Let be the periodic function with the Fourier coefficients defined in (4.49). Then and as .
Proof.
Using (1.5) and , we can write
The function with Fourier coefficients is thus the periodization of , and because of the compact support of , we have for (and large enough)
From Assumption 1.1 and A.4, it readily follows that and as . Now
hence the same holds for since , as will be shown in Lemma A.5 below. This completes the proof. ∎
A.2 Bounds on the cut-off functions
The following Lemma is an easy exercise that we leave to the reader. (Compare with [17, Prop. 4.2].)
Lemma A.5 ( norm of the cut-off functions , ).
The next lemma is somewhat more subtle, as it requires uniformity in .
Lemma A.6 ( norm of the function ).
Let be as in (5.19) for . Then
(A.6) |
Appendix B Evaluation of
In this section we explicitly evaluate the integral in (1.10), and show that it equals with given in (1). In fact, rescaling the integrations over in (1.10) by , we find
(B.1) |
and
(B.2) |
Lemma B.1.
For every , can be written explicitly as in (1).
Proof.
Let us rewrite as
where stands for the characteristic function of the unit ball, and we have used that
because of symmetry. In particular, with
We claim that
which follows from
for fixed and . In fact, the Fourier transform of the integrand above is proportional to
Taking the real part and , this evaluates to at .
We conclude that
By simple scaling, , hence also
(B.3) |
We have
Let us compute
and
which can further be rewritten as
For the integration over , we use that (for )
as well as
Straightforward manipulations then lead to
Appendix C Estimate of the integral in (6.5)
In this section we prove the bound claimed in (6.5).
Lemma C.1.
Proof.
We start by observing that it suffices to integrate only over . For , the integrand is bounded by , hence the resulting integral is finite and bounded uniformly in .
In the following, let us use the notation . Using that and , the integral to consider is bounded above by
Let us split the domain of integration into two domains, depending on whether or , and denote the corresponding integrals by and . In the first case, we have . Dropping the positive term in the denominator, we thus get the upper bound
where we have used that for . We further have, for any ,
(C.1) |
Hence . For the second term we proceed similarly. If we have , with denoting the positive part. We change variables to . Using and we get
which is finite due to (C.1). ∎
References
- [1] G. Basti, S. Cenatiempo, and B. Schlein. A new second-order upper bound for the ground state energy of dilute Bose gases. Forum of Mathematics, Sigma, Vol. 9 (2021), e74.
- [2] N. Benedikter, V. Jaksic, M. Porta, C. Saffirio, and B. Schlein. Mean-field evolution of fermionic mixed states. Comm. Pure Appl. Math. 69 (2016), 2250–2303.
- [3] N. Benedikter, P. T. Nam, M. Porta, B. Schlein, and R. Seiringer. Optimal upper bound for the correlation energy of a Fermi gas in the mean-field regime. Commun. Math. Phys. 374 (2020), 2097–2150.
- [4] N. Benedikter, P. T. Nam, M. Porta, B. Schlein, and R. Seiringer. Correlation energy of a weakly interacting Fermi gas. Invent. Math. 225 (2021), pp. 885–979.
- [5] N. Benedikter, M. Porta, B. Schlein, and R. Seiringer. Correlation Energy of a Weakly Interacting Fermi Gas with Large Interaction Potential. Arch. Rational Mech. Anal. 247 (2023), Art. 65.
- [6] N. Benedikter, M. Porta, and B. Schlein. Mean-field Evolution of Fermionic Systems. Comm. Math. Phys. 331 (2014), 1087–1131.
- [7] H. A. Bethe and J. Goldstone. Effect of a Repulsive Core in the Theory of Complex Nuclei. Proc. R. Soc. Lond. A 238 (1957), 551–567 .
- [8] N. N. Bogoliubov. On the theory of superfluidity. J. Phys. (USSR) 11 (1947), p. 23.
- [9] D. Bohm and D. Pines. A collective description of electron interactions. I. Magnetic interactions. Phys. Rev. 82 (1951), 625–634.
- [10] D. Bohm and D. Pines. A collective description of electron interactions: II. Collective vs. Individual particle aspects of the interactions. Phys. Rev. 85 (1952), 338–353.
- [11] D. Bohm and D. Pines. A collective description of electron interactions: III. Coulomb interactions in a degenerate electron gas. Phys. Rev. 92 (1953), 609–625.
- [12] P. Chankowski and J. Wojtkiewicz. Ground-state energy of the polarized dilute gas of interacting spin- fermions. Phys. Rev. B 104 (2021), 144425.
- [13] M. R. Christiansen, C. Hainzl, and P. T. Nam. The Random Phase Approximation for Interacting Fermi Gases in the Mean-Field Regime. Forum Math. Pi 11 (2023), e32, 1–131.
- [14] M. R. Christiansen, C. Hainzl, and P. T. Nam. The Gell-Mann–Brueckner Formula for the Correlation Energy of the Electron Gas: A Rigorous Upper Bound in the Mean-Field Regime. Commun. Math. Phys. 401 (2023), 1469–1529, .
- [15] M. R. Christiansen, C. Hainzl, and P. T. Nam. The Correlation Energy of the Electron Gas in the Mean-Field Regime. Preprint 2024. arXiv:2405.01386.
- [16] F. J. Dyson. Ground-State Energy of a Hard-Sphere Gas. Phys. Rev. 106 (1957), 20–26.
- [17] M. Falconi, E. L. Giacomelli, C. Hainzl, and M. Porta. The dilute Fermi gas via Bogoliubov theory. Ann. Henri Poincaré 22 (2021), 2283–2353.
- [18] S. Fournais, L. Junge, T. Girardot, L. Morin, M. Olivieri, A. Triay. The free energy of dilute Bose gases at low temperatures interacting via strong potentials. Preprint 2024. arXiv:2408.14222.
- [19] S. Fournais and J. P. Solovej. The energy of dilute Bose gases. Annals of Math. 192 (2020), 893–976.
- [20] S. Fournais and J. P. Solovej. The energy of dilute Bose gases II: the general case. Invent. Math. 232 (2023), 863–994.
- [21] E. L. Giacomelli. An optimal upper bound for the dilute Fermi gas in three dimensions. J. Funct. Anal. 285 (2023), 110073.
- [22] E. L. Giacomelli. An optimal lower bound for the dilute Fermi gas. In preparation.
- [23] E. L. Giacomelli. Bogoliubov Theory for the Dilute Fermi Gas in Three Dimensions, in Quantum Mathematics II, Ed. by M. Correggi and M. Falconi. Springer INdAM Series. Singapore: Springer Nature, 2023, pp. 313–329.
- [24] F. Haberberger, C. Hainzl, B. Schlein and A. Triay. Upper Bound for the Free Energy of Dilute Bose Gases at Low Temperature. Preprint 2024. arXiv:2405.03378.
- [25] F. Haberberger, C. Hainzl, P.T. Nam, R. Seiringer, and A. Triay. The free energy of dilute Bose gases at low temperatures. Preprint 2024. arXiv:2304.02405.
- [26] C. Hainzl. Another Proof of BEC in the GP-limit. J. Math. Phys. 62 (2021), 459-485 .
- [27] C. Hainzl, B. Schlein, and A. Triay. Bogoliubov theory in the Gross-Pitaevskii limit: A simplified approach. Forum of Mathematics, Sigma 10 (2022), E90.
- [28] K. Huang and C. N. Yang. Quantum-mechanical many-body problem with hard-sphere interaction. Physical Review 105 (1957), 767–775.
- [29] S. Kanno. Criterion for the Ferromagnetism of Hard Sphere Fermi Liquid. II. Prog. Theor. Phys. 44 (1970), 813–815.
- [30] A. B. Lauritsen. Almost optimal upper bound for the ground state energy of a dilute Fermi gas via cluster expansion. Preprint 2023. arXiv:2301.08005.
- [31] A. B. Lauritsen and R. Seiringer. Ground state energy of the dilute spin-polarized Fermi gas: Upper bound via cluster expansion. J. Funct. Anal. 286 (2024), 110320.
- [32] A. B. Lauritsen and R. Seiringer. Ground state energy of the dilute spin-polarized Fermi gas: lower bound. Preprint 2024. arXiv:2402.17558.
- [33] T. D. Lee, K. Huang, and C. N. Yang. Eigenvalues and eigenfunctions of a Bose system of hard spheres and its low-temperature properties. Phys. Rev. 106 (1957), 1135–1145.
- [34] E. Lieb, R. Seiringer, and J. P. Solovej. Ground state energy of the low density Fermi gas. Physical Review A 71 (2005), 053605.
- [35] E.H. Lieb, R. Seiringer, J.P. Solovej, and J. Yngvason. The Mathematics of the Bose Gas and its Condensation. Birkhäuser, Basel (2005).
- [36] E. H. Lieb and J. Yngvason. Ground State Energy of the Low Density Bose Gas. Phys. Rev. Lett. 80 (1998), 2504–2507.
- [37] E.H. Lieb and J. Yngvason. The Ground State Energy of a Dilute Two-Dimensional Bose Gas. J. Stat. Phys. 103 (2001), 509–526.
- [38] N. Navon, S. Nascimbene, F. Chevy, and C. Salomon. The equation of state of a low temperature Fermi gas with tunable interactions. Science 328 (2010), 729–732.
- [39] J. Pera, J. Casulleras, and J. Boronat. Itinerant ferromagnetism in dilute Fermi gases. SciPost Phys. 14 (2023), 038.
- [40] D. Pines. A Collective Description of Electron Interactions: IV. Electron Interaction in Metals. Phys. Rev. 92 (1953), 626.
- [41] D. Ruelle. Statistical Mechanics. Rigorous Results. World Scientific (1999).
- [42] D. Robinson. The Thermodynamic Pressure in Quantum Statistical Mechanics. Springer Lecture Notes in Physics, Vol. 9 (1971).
- [43] K. Sawada. Correlation energy of an electron gas at high density. Phys. Rev. 106 (1957), 372.
- [44] K. Sawada, K. A. Brueckner, N. Fukuda, and R. Brout. Correlation energy of an electron gas at high density: Plasma oscillations. Phys. Rev. 108 (1957), 507.
- [45] R. Seiringer. The Thermodynamic Pressure of a Dilute Fermi Gas. Commun. Math. Phys. 261 (2006), 729–758.
- [46] J. P. Solovej. Many Body Quantum Mechanics. Lecture Notes (2007). Available online at: http://web.math.ku.dk/~solovej/MANYBODY/mbnotes-ptn-5-3-14.pdf
- [47] H.-T. Yau and J. Yin. The second order upper bound for the ground energy of a Bose gas. J. Stat. Phys. 136 (2009), 453–503.