Lode Pollet

The ALPS project release 2.0: open source software for strongly correlated systems

We present release 2.0 of the ALPS (Algorithms and Libraries for Physics Simulations) project, an open source software project to develop libraries and application programs for the simulation of strongly correlated quantum lattice models such as quantum magnets, lattice bosons, and strongly correlated fermion systems. The code development is centered on common XML and HDF5 data formats, libraries to simplify and speed up code development, common evaluation and plotting tools, and simulation programs. The programs enable non-experts to start carrying out serial or parallel numerical simulations by providing basic implementations of the important algorithms for quantum lattice models: classical and quantum Monte Carlo (QMC) using non-local updates, extended ensemble simulations, exact and full diagonalization (ED), the density matrix renormalization group (DMRG) both in a static version and a dynamic time-evolving block decimation (TEBD) code, and quantum Monte Carlo solvers for dynamical mean field theory (DMFT). The ALPS libraries provide a powerful framework for programmers to develop their own applications, which, for instance, greatly simplify the steps of porting a serial code onto a parallel, distributed memory machine. Major changes in release 2.0 include the use of HDF5 for binary data, evaluation tools in Python, support for the Windows operating system, the use of CMake as build system and binary installation packages for Mac OS X and Windows, and integration with the VisTrails workflow provenance tool. The software is available from our web server at http://alps.comp-phys.org/.

Superfluidity of Grain Boundaries in Solid 4He

By large-scale quantum Monte Carlo simulations we show that grain boundaries in 4He crystals are generically superfluid at low temperature, with a transition temperature of the order of approximately 0.5 K at the melting pressure; nonsuperfluid grain boundaries are found only for special orientations of the grains. We also find that close vicinity to the melting line is not a necessary condition for superfluid grain boundaries, and a grain boundary in direct contact with the superfluid liquid at the melting curve is found to be mechanically stable and the grain-boundary superfluidity observed by Sasaki et al. [Science 313, 1098 (2006)10.1126/science.1130879] is not just a crack filled with superfluid.

Quantitative determination of temperature in the approach to magnetic order of ultracold fermions in an optical lattice.

We perform a quantitative simulation of the repulsive Fermi-Hubbard model using an ultracold gas trapped in an optical lattice. The entropy of the system is determined by comparing accurate measurements of the equilibrium double occupancy with theoretical calculations over a wide range of parameters. We demonstrate the applicability of both high-temperature series and dynamical mean-field theory to obtain quantitative agreement with the experimental data. The reliability of the entropy determination is confirmed by a comprehensive analysis of all systematic errors. In the center of the Mott insulating cloud we obtain an entropy per atom as low as 0.77k(B) which is about twice as large as the entropy at the Néel transition. The corresponding temperature depends on the atom number and for small fillings reaches values on the order of the tunneling energy.

Competition between pairing and ferromagnetic instabilities in ultracold Fermi gases near Feshbach resonances.

We study the quench dynamics of a two-component ultracold Fermi gas from the weak into the strong interaction regime, where the short time dynamics are governed by the exponential growth rate of unstable collective modes. We obtain an effective interaction that takes into account both Pauli blocking and the energy dependence of the scattering amplitude near a Feshbach resonance. Using this interaction we analyze the competing instabilities towards Stoner ferromagnetism and pairing.

Local Stress and Superfluid Properties of Solid 4He

We provide a semiquantitative tool, derived from first-principles simulations, for answering the question of whether certain types of defects in solid 4He support mass superflow. Although ideal crystals of 4He are not supersolid, the gap for vacancy creation closes when applying a moderate stress. While a homogeneous system becomes unstable at this point, the stressed core of crystalline defects (dislocations and grain boundaries) can turn superfluid.

Phase diagram of the disordered Bose-Hubbard model

We establish the phase diagram of the disordered three-dimensional Bose-Hubbard model at unity filling which has been controversial for many years. The theorem of inclusions, proven by Pollet et al. [Phys. Rev. Lett. 103, 140402 (2009)] states that the Bose-glass phase always intervenes between the Mott insulating and superfluid phases. Here, we note that assumptions on which the theorem is based exclude phase transitions between gapped (Mott insulator) and gapless phases (Bose glass). The apparent paradox is resolved through a unique mechanism: such transitions have to be of the Griffiths type when the vanishing of the gap at the critical point is due to a zero concentration of rare regions where extreme fluctuations of disorder mimic a regular gapless system. An exactly solvable random transverse field Ising model in one dimension is used to illustrate the point. A highly nontrivial overall shape of the phase diagram is revealed with the worm algorithm. The phase diagram features a long superfluid finger at strong disorder and on-site interaction. Moreover, bosonic superfluidity is extremely robust against disorder in a broad range of interaction parameters; it persists in random potentials nearly 50 (!) times larger than the particle half-bandwidth. Finally, we comment on the feasibility of obtaining this phase diagram in cold-atom experiments, which work with trapped systems at finite temperature.

Thermodynamics of the 3D Hubbard Model on Approaching the Neel Transition

We study the thermodynamic properties of the 3D Hubbard model for temperatures down to the Néel temperature by using cluster dynamical mean-field theory. In particular, we calculate the energy, entropy, density, double occupancy, and nearest-neighbor spin correlations as a function of chemical potential, temperature, and repulsion strength. To make contact with cold-gas experiments, we also compute properties of the system subject to an external trap in the local density approximation. We find that an entropy per particle S/N ≈ 0.65(6) at U/t = 8 is sufficient to achieve a Néel state in the center of the trap, substantially higher than the entropy required in a homogeneous system. Precursors to antiferromagnetism can clearly be observed in nearest-neighbor spin correlators.

Diagrammatic Monte Carlo for correlated fermions

We show that Monte Carlo sampling of the Feynman diagrammatic series (DiagMC) can be used for tackling hard fermionic quantum many-body problems in the thermodynamic limit by presenting accurate results for the repulsive Hubbard model in the correlated Fermi liquid regime. Sampling Feynmans diagrammatic series for the single-particle self-energy, we can study moderate values of the on-site repulsion (U/t~4) and temperatures down to T/t=1/40. We compare our results with high-temperature series expansions and with single-site and cluster dynamical mean-field theory.

Temperature changes when adiabatically ramping up an optical lattice

When atoms are loaded into an optical lattice, the process of gradually turning on the lattice is almost adiabatic. In this paper, we investigate how the temperature changes when going from the gapless superfluid phase to the gapped Mott phase along isentropic lines. To do so we calculate the entropy in the single-band Bose–Hubbard model for various densities, interaction strengths and temperatures in one and two dimensions for homogeneous and trapped systems. Our theory is able to reproduce the experimentally observed visibilities and therefore strongly supports the view that current experiments remain in the quantum regime for the considered lattice depths with low temperatures and minimal heating.

Dynamical mean-field theory for bosons

We discuss the recently developed bosonic dynamical mean-field theory (B-DMFT) framework, which maps a bosonic lattice model onto the self-consistent solution of a bosonic impurity model with coupling to a reservoir of normal and condensed bosons. The effective impurity action is derived in several ways: (i) as an approximation to the kinetic energy functional of the lattice problem, (ii) using a cavity approach and (iii) using an effective medium approach based on adding a one-loop correction to the self-consistently defined condensate. To solve the impurity problem, we use a continuous-time Monte Carlo algorithm based on the sampling of a perturbation expansion in the hybridization functions and the condensate wave function. As applications of the formalism, we present finite-temperature B-DMFT phase diagrams for the bosonic Hubbard model on a three-dimensional (3D) cubic and a 2D square lattice, the condensate order parameter as a function of chemical potential, critical exponents for the condensate, the approach to the weakly interacting Bose gas regime for weak repulsions and the kinetic energy as a function of temperature.