Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDensity functional theory versus quantum Monte Carlo simulations of Fermi gases in the optical-lattice arena
93
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We benchmark the ground state energies and the density profiles of atomic repulsive Fermi gases in optical lattices computed via Density Functional Theory (DFT) against the results of diffusion Monte Carlo (DMC) simulations. The main focus is on a half-filled one-dimensional optical lattices, for which the DMC simulations performed within the fixed-node approach provide unbiased results. This allows us to demonstrate that the local spin-density approximation (LSDA) to the exchange-correlation functional of DFT is very accurate in the weak and intermediate interactions regime, and also to underline its limitations close to the strongly-interacting Tonks-Girardeau limit and in very deep optical lattices. We also consider a three dimensional optical lattice at quarter filling, showing also in this case the high accuracy of the LSDA in the moderate interaction regime. The one-dimensional data provided in this study may represent a useful benchmark to further develop DFT methods beyond the LSDA and they will hopefully motivate experimental studies to accurately measure the equation of state of Fermi gases in higher-dimensional geometries.
rate research
Read More
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.
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.
The dynamics of a one-dimensional two-component Fermi gas in the presence of a quasi-periodic optical lattice (OL) is investigated by means of a Density Functional Theory approach. Inspired by the protocol implemented in recent cold-atom experiments, designed to identify the many-body localization transition, we analyze the relaxation of an initially prepared imbalance between the occupation number of odd and of even sites. For quasi-disorder strength beyond the Anderson localization transition, the imbalance survives for long times, indicating the inability of the system to reach local equilibrium. The late time value diminishes for increasing interaction strength. Close to the critical quasi-disorder strength corresponding to the noninteracting (Anderson) transition, the interacting system displays an extremely slow relaxation dynamics, consistent with sub-diffusive behavior. The amplitude of the imbalance fluctuations around its running average is found to decrease with time, and such damping is more effective with increasing interaction strengths. While our study addresses the setup with two equally intense OLs, very similar effects due to interactions have been observed also in recent cold-atom experiments performed in the tight-binding regime, i.e. where one of the two OLs is very deep and the other is much weaker.
We present in detail two variants of the lattice Monte Carlo method aimed at tackling systems in external trapping potentials: a uniform-lattice approach with hard-wall boundary conditions, and a non-uniform Gauss-Hermite lattice approach. Using those two methods, we compute the ground-state energy and spatial density profile for systems of N=4 - 8 harmonically trapped fermions in one dimension. From the favorable comparison of both energies and density profiles (particularly in regions of low density), we conclude that the trapping potential is properly resolved by the hard-wall basis. Our work paves the way to higher dimensions and finite temperature analyses, as calculations with the hard-wall basis can be accelerated via fast Fourier transforms, the cost of unaccelerated methods is otherwise prohibitive due to the unfavorable scaling with system size.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
