Study of dynamic and ground-state propertiesof dipolar Fermi gases using mean-field and quantum Monte Carlo methods

N.M. Matveeva

Published Feb 18, 2013• 0 citations

Abstract

In this thesis I theoretically study the dynamic and ground state properties of ultracold dipolar Fermi gases. The mean-field approach based on the Thomas-Fermi energy functional is applied to consider the dynamic properties of bilayer harmonically trapped dipolar Fermi gases. The fixed-node Diffusion Monte Carlo method (FNDMC) is used instead to investigate the ground-state properties of two dimensional dipolar Fermi gases. This technique is also applied to the problem of one impurity in a bilayer configuration with dipolar fermions.

Related papers

Study of Ultracold Fermi Gases in the BCS-BEC Crossover: Quantum Monte Carlo Methods, Hydrodynamics and Local Density Approximation.

Gianluca Bertaina

In this Thesis we will theoretically address some issues concerning the Physics of ultracold Fermi gases, all of them with experimental relevance. The field of ultracold gases, and more recently of ultracold Fermi gases, is gathering a lot of experimental and theoretical interest for two main reasons: the great control on the relevant parameters of the problem and the relative simplicity of the minimal theory which is able to correctly describe the system. Building upon this solid background, ultracold gases provided an ideal laboratory for testing more refined theories and also for addressing fundamental issues of quantum mechanics as well as for simulating more complex physical systems such as those encountered in condensed matter. Even if the effective hamiltonian of ultracold gases can be simple, due to the diluteness of the system and the low temperature, which imply low energy physics, the solution of the quantum mechanic equations governing the state of these systems is not always simple. For the static properties usually mean-field solutions exist or perturbative expansions can be produced in some regimes. However Quantum Monte Carlo (QMC) techniques provide more accurate results especially in the strongly interacting regimes. For confined systems it is possible to use QMC only for a few particles, so that, for large number of particles, a fruitful combined use of Density Functional Theory (DFT) and QMC is necessary. The study of the dynamics of ultracold gases has received little attention with QMC techniques, due to the intrinsic computational difficulty of the many-body problem, so that general hydrodynamic equations are often used for studying the propagation of smoothly varying perturbations. In this Thesis we use QMC techniques for studying the problem of ferromagnetism in repulsive or effectively repulsive ultracold gases without a lattice and the problem of the Bardeen-Cooper-Schrieffer to Bose-Einstein-Condensation (BCS-BEC) crossover in two dimensions. We use DFT in the Local Density Approximation (LDA) for calculating the density profiles of ultracold Fermi gases in harmonic magneto-optical traps, starting from QMC equations of state. We study the propagation of first and second sound in ultracold Fermi gases in cylindrical geometry, using the hydrodynamic equations of superfluids.

Jan 1, 20101 citation

Quantum Monte Carlo study of ultracold gases (PhD thesis)

G. E. Astrakharchik

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.

arXiv:1412.4529Dec 15, 20144 citations

Ultracold Bose and Fermi dipolar gases : a quantum Monte Carlo study

Raúl Bombín

arXiv:2002.01378Feb 4, 20201 citation

Effects of the dipole-dipole interaction on the physics of ultracold quantum gases

Marta Abad García

In this thesis we study the effects of the dipole-dipole interaction on the physics of ultracold quantum gases, both bosonic and fermionic, within the theoretical framework provided by the mean-field regime. This kind of interaction takes place in ultracold atomic gases (for instance 52Cr or 164Dy) due to their atomic magnetic dipole moment, and in ultracold molecular gases due to the magnetic or electric dipole moment. \n \nIn the case of quantum gases of bosonic atoms, or Bose-Einstein condensates, the dipole-dipole interaction can be studied within mean-field approximation using the Gross-Pitaevskii equation, which now contains a new non-linear term due to the dipole-dipole interaction. We investigate, on the one hand dipolar condensates confined in harmonic traps, and on the other dipolar condensates confined in toroidal traps. In the harmonic geometry, our focus is on the study of the ground state and the quantized vortex state, where the density profile is characterized as well as some properties leading to the process of vortex formation, such as the critical frequency and the energy barrier that has to be overcome to bring the vortex from the surface to the centre of the gas. We finish the study of dipolar condensates in harmonic traps by dynamically simulating the precession frequency of an off-center vortex in a non-rotating condensate. \n \nIn the toroidal geometry the dipolar effects are strongly magnified when the polarization axis of the dipoles is perpendicular to the trap symmetry axis. In this case, the anisotropic structure of the density can be understood as the response of the system to the double-well effective potential along the ring. We have studied the dynamics of this system when the initial number of atoms in the left and right wells is imbalanced, predicting Josephson and self-trapping oscillations depending on the initial condition. This has led us to name this new system as Self-induced Josephson Junction. We have studied in detail the self-trapping regime and we have seen that the particle flux inversion is closely related to the crossing of vortices across the Josephson junctions. This result opens the door to establishing a more direct connection between the phase-slip regime, widely addressed in superfluid helium, and the self-trapping regime of condensates. \n \nIn the case of quantum gases of fermionic dipolar particles, we have studied how the radial quadrupole mode allows one to distinguish between hydrodynamic and collisionless regimes. We have analytically calculated the frequency of this mode in the mean-field approximation, generalizing the results from the Thomas-Fermi approximation for trapped ideal Fermi gases. On the one hand, we observe that the frequency in the hydrodynamic regime is smaller than in non-dipolar Fermi gases, while in the collisionless regime the frequency is larger or smaller than that corresponding to the non-interacting system depending on the geometry of the harmonic trap. On the other hand, we predict that reducing the trap deformation (aspect ratio) an observable jump in the frequency of the radial quadrupole mode would take place, which would correspond to the transition between the collisionless and hydrodynamic regimes, for instance when the gas undergoes the transition to the superfluid state.

Feb 16, 20120 citations

Quantum Monte Carlo study of low-dimensional Fermi fluids of dipolar atoms

Clio Johnson, N. D. Drummond, J. P. Hague +1 more

Fermionic cold atoms in optical traps provide viable quantum simulators of correlation effects in electronic systems. For dressed Rydberg atoms in two-dimensional traps with out-of-plane dipole moments, a realistic model of the pairwise interaction is of repulsive dipolar 1 / r 3 form at long range, softened to a constant at short range. This study provides parametrizations of fixed-node diffusion Monte Carlo energy data for ferromagnetic (one-component) and paramagnetic (two-component) two-dimensional homogeneous Fermi fluids of interacting dipolar atoms. We find itinerant ferromagnetism to be unstable within our parameter spaces for dipolar interactions both with and without softening. Our parametrization of the energy as a function of density will enable density functional theory to support experimental studies of inhomogeneous fermionic cold-atom systems.

arXiv:2510.19533Physical Review ADec 3, 20250 citations

Ultracold Bose and Fermi dipolar gases : a quantum Monte Carlo study

Raúl Bombín

The object of study of this thesis are dipolar systems in the quantum degenerate regime. In general, dealing with many-body systems and evaluating their properties requires to deal with the the Schrödinger equation. In the present study we employ different Monte Carlo methods that allow to find numerical solutions to it by employing a set of stochastic techniques. The simplest one that we introduce corresponds to the Variational Monte Carlo (VMC) method, that despite its simplicity, allows to obtain variational solutions to the many-body problem. A more accurate description is provided by Diffusion Monte Carlo (DMC), that provides exact solutions for the ground state of the system when dealing with bosons. We continue presenting two methods that rely on the Feynman's path integral formalism of quantum mechanics: Path Integral Monte Carlo (PIMC) and Path Integral Ground State (PIGS), that provide exact solutions for the bosonic problem at finite and zero temperature respectively. In order to work with fermionic systems, as we do in chapter 4 of this thesis, the DMC algorithm has to be modified with the Fixed-Node (FN) approximation, what alows to avoid the sign problem. Doing so, the results obtained with DMC correspond to variational upper bounds to the energy. In chapter 3 we study the superfluid properties of a system of dipolar bosons that are fully polarized and in which the atoms are restricted to move in the plane. We also consider that all the dipolar moments form a certain tilting angle with the axis perpendicular to the plane, what allows to introduce anisotropy in the system. The phase diagram at zero temperature of this system reveals the existence of three different phases: gas, stripe and solid. Here we focus on the characterization of the superfluid properties across that phase diagram. Our calculations allow to address the question of whether the stripe phase of this system could be a candidate for the supersolid, a system that simultaneously exhibit spatial long-range order and superfluidity. By the employment of DMC and PIGS, we report finite supefluid and condensate fractions, both in the gas and the stripe phases. Then, the study is completed by performing finite temperature calculations, where the use of PIMC allows to characterize the BKT transition and to report the critical temperature at which it occurs in the different phases. Finally, by direct comparison with the predictions of the Luttinger Liquid theory, we explicitly show that the stripe phase can not be described as an ensemble of 1D isolated systems. In chapter 4, we study the fermionic dipolar system in two dimensions, focusing in the case in which all dipoles are polarized along the direction, that in this case is chosen to be the one perpendicular to the plane containing their movement. We compute the equation of state of the system in a wide range of interaction parameters. In the low density regime, the comparison of our results for the dipolar model with those of a hard-disks one allows to determine the regime of universality. On the other hand, at higher densities ( and before crystallization), we discuss the issue of itinerant ferromagnetism, that is, the possibility of having a polarized phase as the ground state of the system. The repulsive Fermi polaron with dipolar interaction, that corresponds to the limit of ultralow concentration of impurities embebed in a fermionic bath is also studied. Here we determine the regime of universality for this problem and compute observables that allow to discuss the validity of the quasi-particle picture. In the last part of the thesis, the formation of ultra-dilute dipolar droplets is studied. Our results are in agreement with experimental measurements performed with dysprosium atoms. On the other hand, the evaluation of their differences with the prediction of the extended Gross-Pitaevskii equation makes it possible to determine the limits of the mean-field approach to this problem. El objeto de estudio de esta tesis son los sistemas dipolares en el régimen cuántico degenardo. Usualmente, tratar con sistemas de muchos cuerpos y para evaluar sus propiedades requiere ser capaz de resolver la ecuación de Schrödinger. En el presente estudio, empleamos diferentes métodos de Monte Carlo, que permiten encontrar soluciones numéricas de forma estocásticas. La primera y más simple de estas técnicas es el método Variational Monte Carlo (VMC), que da una solución variacional. Una mejora sobre lo anterior consiste en emplear el método Diffusion Monte Carlo (DMC) que permite obtener soluciones exactas para el estado fundamental del sistema (cuando se estudian sistemas bosónicos). Continuamos presentando dos métodos que se basan en el formalismo Feynman de la mecánica cuántica: Path Integral Monte Carlo (PIMC) y Path Integral Ground State (PIGS), que proporcionan soluciones exactas para el problema bosónico a temperatura finita y en el límite de temperatura cero respectivamente. Para trabajar con sistemas fermionicos, como es el caso del capítulo 4 de esta tesis, el algoritmo DMC tiene que ser modificado según la prescripción Fixed-Node para evitar el problema del signo. Al hacerlo, los resultados obtenidos con DMC se corresponden a soluciones variacionales a la energía. En el capítulo 3 estudiamos las propiedades superfluidas de un sistema de bosones dipolares completamente polarizados y en el que el movimiento de los dipolos está restringido al plano. También consideramos que los momentos dipolares forman un cierto ángulo con el eje perpendicular al plano, lo que permite introducir anisotropía en el sistema. El diagrama de fases a temperatura cero de este sistema revela la existencia de tres fases diferentes: gas, stripe y sólido. Aquí nos centramos en la caracterización de las propiedades superfluidas en ese diagrama de fases. Nuestros cálculos permiten abordar la cuestión de si la fase stripe de este sistema podría ser un candidato para el supersólido: un sistema en el que dos simetrías U (1) se rompen simultáneamente, permitiendo al sistema exhibir orden espacial de largo alcance y a la vez ser superfluido. Mediante el empleo de DMC y PIGS, evaluamos la fracción superfluída y el condensado, tanto en las fases de gas como en el stripe. Este estudio se completa con la extensión a temperatura finita, donde el uso de PIMC permite caracterizar la transición superfluida y obtener la temperatura crítica a la que ésta ocurre en las fases gas y stripe. Finalmente, por comparación directa con las predicciones Líquido de Luttinger, mostramos explícitamente que la fase de stripe no puede describirse como un conjunto de sistemas 1D aislados. En el capítulo 4, estudiamos el sistema dipolar fermiónico en dos dimensiones, enfocándonos en el caso en que todos los dipolos están polarizados a lo largo de la dirección que es perpendicular al plano que contiene su movimiento. Calculamos la ecuación de estado del sistema en un amplio rango de parámetros de interacción: a baja densidad, la comparación de nuestro modelo dipolar con uno de discos duros permite determinar el régimen de universalidad, mientras que a densidades más altas (antes de la cristalización), discutimos la posibilidad de una fase polarizada como el estado fundamental del sistema (ferromagnetismo itinerante). El polaron fermiónico dipolar, correspondiente al límite de impurezas ultradiluídas en un baño fermiónico también es estudiado, determinando el régimen de universalidad y los límites de validez de la aproximación de quasi-partícula. En la última parte de la tesis, la formación de gotas dipolares ultradiluídas es estudiada. Nuestros resultados están en acuerdo con medidas experimentales con átomos de disprosio. Por otro lado, la evaluación de las diferencias entre éstos y la predicción dada por la ecuación de Gross-Pitaevskii extendida, permite evaluar los límites de la aproximación

Dec 13, 20192 citations