No Arabic abstract
This Dissertation presents results of a thorough study of ultracold bosonic and fermionic gases in three-dimensional and quasi-one-dimensional systems. Although the analyses are carried out within various theoretical frameworks (Gross-Pitaevskii, Bethe ansatz, local density approximation, etc.) the main tool of the study is the Quantum Monte Carlo method in different modifications (variational Monte Carlo, diffusion Monte Carlo, fixed-node Monte Carlo methods). We benchmark our Monte Carlo calculations by recovering known analytical results (perturbative theories in dilute limits, exactly solvable models, etc.) and extend calculations to regimes, where the results are so far unknown. In particular we calculate the equation of state and correlation functions for gases in various geometries and with various interatomic interactions.
We present two Diffusion Monte Carlo (DMC) algorithms for systems of ultracold quantum gases featuring synthetic spin-orbit interactions. The first one is a discrete spin generalization of the T- moves spin-orbit DMC, which provides an upper bound to the fixed-phase energy. The second is a spin-integrated DMC method which recovers the fixed-phase property by avoiding the definition of the effective Hamiltonian involved in the T-moves approach. The latter is a more accurate method but it is restricted to spin-independent two-body interactions. We report a comparison between both algorithms for different systems. As a check of the efficiency of both methods, we compare the DMC energies with results obtained with other numerical methods, finding agreement between both estimation
The ground state properties of spin-polarized deuterium (D$downarrow$) at zero temperature are obtained by means of the diffusion Monte Carlo calculations within the fixed-node approximation. Three D$downarrow$ species have been investigated (D$downarrow_1$, D$downarrow_2$, D$downarrow_3$), corresponding respectively to one, two and three equally occupied nuclear spin states. Influence of the backflow correlations on the ground state energy of the systems is explored. The equilibrium densities for D$downarrow_2$ and D$downarrow_3$ liquids are obtained and compared with ones obtained in previous approximate prediction. The density and the pressure at which the gas-liquid phase transition occurs at $T$=0 is obtained for D$downarrow_1$.
The ground-state properties of two-component repulsive Fermi gases in two dimensions are investigated by means of fixed-node diffusion Monte Carlo simulations. The energy per particle is determined as a function of the intercomponent interaction strength and of the population imbalance. The regime of universality in terms of the s-wave scattering length is identified by comparing results for hard-disk and for soft-disk potentials. In the large imbalance regime, the equation of state turns out to be well described by a Landau-Pomeranchuk functional for two-dimensional polarons. To fully characterize this expansion, we determine the polarons effective mass and their coupling parameter, complementing previous studies on their chemical potential. Furthermore, we extract the magnetic susceptibility from low-imbalance data, finding only small deviations from the mean-field prediction. While the mean-field theory predicts a direct transition from a paramagnetic to a fully ferromagnetic phase, our diffusion Monte Carlo results suggest that the partially ferromagnetic phase is stable in a narrow interval of the interaction parameter. This finding calls for further analyses on the effects due to the fixed-node constraint.
This is a review of recent developments in Monte Carlo methods in the field of ultra cold gases. For bosonic atoms in an optical lattice we discuss path integral Monte Carlo simulations with worm updates and show the excellent agreement with cold atom experiments. We also review recent progress in simulating bosonic systems with long-range interactions, disordered bosons, mixtures of bosons, and spinful bosonic systems. For repulsive fermionic systems determinantal methods at half filling are sign free, but in general no sign-free method exists. We review the developments in diagrammatic Monte Carlo for the Fermi polaron problem and the Hubbard model, and show the connection with dynamical mean-field theory. We end the review with diffusion Monte Carlo for the Stoner problem in cold gases.
We study a resonant Bose-Fermi mixture at zero temperature by using the fixed-node diffusion Monte Carlo method. We explore the system from weak to strong boson-fermion interaction, for different concentrations of the bosons relative to the fermion component. We focus on the case where the boson density $n_B$ is smaller than the fermion density $n_F$, for which a first-order quantum phase transition is found from a state with condensed bosons immersed in a Fermi sea, to a Fermi-Fermi mixture of composite fermions and unpaired fermions. We obtain the equation of state and the phase diagram, and we find that the region of phase separation shrinks to zero for vanishing $n_B$.