The programme of reducing N -particle problems to one-particle problems is extended to include the Pauli principle. To calculate the ground-state energy of a system of N fermions interacting by pair forces, a rigorous lower-bound shell model is derived. This shell model follows a building-up principle, and tends to improve with increasing number of particles as well as with increasing strength of interaction. For twenty particles antisymmetric in ordinary space interacting by square-well interaction of strength V 0 a 2 = 200 hslash 2 /2 m the shell model differs by less than 8% from the calculated upper bound, and hence a fortiori from the exact energy. In order to test the quality of the approximation for various interactions, it was necessary to calculate upper bounds.
The object of this thesis is to study the definition and the mathematically rigorous derivation of effective models in many body quantum mechanics. In this field, it is often impossible to compte exactly physically relevant quantities - such as the ground state energy of a system - because of the large number of degrees of freedom. Hence, introducing effective models (depending only on the local density of particles for instance) allows us to compute these quantities. Then, it remains to check that the effective model is indeed an approximation of the exact model in a macroscopic limit. Besides, in some models, some physical quantities (energy, scattering amplitude...), when computed naively, appear to be infinite. One has then to renormalize the theory in order to make these divergences (that are not observed experimentally) disappear. After putting the results of this thesis into context in the introduction, we will focus on the renormalization of a bosonic quadratic Hamiltonian on the Fock space in Chapter 2. This (formal) Hamiltonian is the sum of the second quantization of a one-particle operator and the square of a field operator, and it exhibits an ultraviolet divergence (i.e. coming from local singularities). More accurately, we will try to define a self-adjoint operator, and to show that it is the limit (in the resolvent sense) of the Hamiltonian with ultraviolet cutoff when the cutoff tends to infinity. To do so, we will first diagonalize the Hamiltonian in the regular case, using Bogoliubov transformations. Then, we will prove that the diagonalized Hamiltonian (up to a renormalization of the energy and the charge) still defines a self-adjoint, and we will show the convergence of the Hamiltonian with cut-off. In a second time, we will study the semi-classical limit of a Fermi gas, in two different scalings. In Chapter 3, we will focus on the dilute limite (i.e. with an interaction range very small compared to the typical distance between particles) of an inhomogeneous (trapped with a potential V) Fermi gas, with repulsive interactions and in three dimensions. Due to the Pauli exclusion principle, two fermions of same spin can only rarely be close enough to interact. Thus, we will consider a Fermi gas with two spin components. We will then show that the ground state energy of such a system can be approximated by a Thomas-Fermi energy (in which interactions play the role of a perturbation of the interactionless energy) in the large number of particles limit. Eventually, in Chapter 4, we will study the mean-field limit (with frequent but weak interactions) of an inhomogeneous Fermi gas, with attractive interactions and in one or two dimensions this time. Once more, we will show the convergence of the ground state energy toward a Thomas-Fermi energy. Moreover, we will also prove the convergence of the ground states. More accurately, we will show that the Husimi functions of these ground states converge to a convex combination of minimizers of the Vlasov energy.
The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid
He3, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$.
F is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the
f(rij) embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements
Hmn,
Nmn of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the
Hmn are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (
Hmn) and (
Nmn), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (
Hmn) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (
Hmn). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (
Hmn). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any
F in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{},
Q-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid
He3 has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.
Sums of matrix elements of spin-dependent two-body momentum-independent interactions and sums of their products are calculated analytically in the basis of many-body states with given total spin---the states built from spin and spatial wave functions belonging to multidimensional irreducible representations of the symmetric group, unless the total spin has the maximal allowed value. As in the first part of the series [V. A. Yurovsky, Phys. Rev. A 91, 053601 (2015)], the sum dependence on the many-body states is given by universal factors, which are independent of the Hamiltonians of noninteracting particles. The sum rules are applied to perturbative analysis of energy spectra and to calculation of two-body spin-dependent local correlations.