Source author record

Lode Pollet

Lode Pollet appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2025arXiv

Assembling a Bose-Hubbard superfluid from tweezer-controlled single atoms

Quantum simulation relies on the preparation and control of low-entropy many-body systems to reveal the behavior of classically intractable models. The development of new approaches for realizing such systems therefore represents a frontier in quantum science. Here we experimentally demonstrate a new protocol for generating ultracold, itinerant many-body states in a tunnel-coupled two-dimensional optical lattice. We do this by adiabatically connecting a near-ground-state-cooled array of up to 50 single strontium-86 atoms with a Bose-Hubbard superfluid. Through comparison with finite-temperature quantum-Monte-Carlo calculations, we estimate that the entropy per particle of the prepared many-body states is approximately $2 k_B$, and that the achieved temperatures are consistent with a significant superfluid fraction. This represents the first time that itinerant many-body systems have been prepared from rearranged atoms, opening the door to bottom-up assembly of a wide range of neutral-atom and molecular systems.

preprint2022arXiv

Supertransport by Superclimbing Dislocations in $^4$He: When All Dimensions Matter

The unique superflow-through-solid effect observed in solid Helium-4 and attributed to the quasi-one-dimensional superfluidity along the dislocation cores exhibits two extraordinary features: (i) an exponentially strong suppression of the flow by a moderate increase in pressure, and (ii) an unusual temperature dependence of the flow rate with no analogy to any known system and in contradiction with the standard Luttinger liquid paradigm. Based on ab initio and model simulations, we argue that the two features are closely related: Thermal fluctuations of the shape of a superclimbing edge dislocation induce large, correlated, and asymmetric stress fields acting on the superfluid core. The critical flux is most sensitive to strong rare fluctuations and hereby acquires a sharp temperature dependence observed in experiments.

preprint2020arXiv

The view of TK-SVM on the phase hierarchy in the classical kagome Heisenberg antiferromagnet

We illustrate how the tensorial kernel support vector machine (TK-SVM) can probe the hidden multipolar orders and emergent local constraint in the classical kagome Heisenberg antiferromagnet. We show that TK-SVM learns the finite-temperature phase diagram in an unsupervised way. Moreover, in virtue of its strong interpretability, it identifies the tensorial quadrupolar and octupolar orders, which define a biaxial $D_{3h}$ spin nematic, and the local constraint that underlies the selection of coplanar states. We then discuss the disorder hierarchy of the phases, which can be inferred from both the analytical order parameters and a SVM bias parameter. For completeness we mention that the machine also picks up the leading $\sqrt{3} \times \sqrt{3}$ correlations in the dipolar channel at very low temperature, which are however weak compared to the quadrupolar and octupolar orders. Our work shows how TK-SVM can facilitate and speed up the analysis of classical frustrated magnets.

preprint2020arXiv

Z2 parton phases in the mixed-dimensional $t-J_z$ model

We study the interplay of spin- and charge degrees of freedom in a doped Ising antiferromagnet, where the motion of charges is restricted to one dimension. The phase diagram of this mixed-dimensional $t - J_z$ model can be understood in terms of spin-less chargons coupled to a Z2 lattice gauge field. The antiferromagnetic couplings give rise to interactions between Z2 electric field lines which, in turn, lead to a robust stripe phase at low temperatures. At higher temperatures, a confined meson-gas phase is found for low doping whereas at higher doping values, a robust deconfined chargon-gas phase is seen which features hidden antiferromagnetic order. We confirm these phases in quantum Monte Carlo simulations. Our model can be implemented and its phases detected with existing technology in ultracold atom experiments. The critical temperature for stripe formation with a sufficiently high hole concentration is around the spin-exchange energy $J_z$, i.e., well within reach of current experiments.

preprint2019arXiv

Constrained Random Phase Approximation of the effective Coulomb interaction in lattice models of twisted bilayer graphene

Recent experiments on twisted bilayer graphene show the urgent need for establishing a low-energy lattice model for the system. We use the constrained random phase approximation to study the interaction parameters of such models taking into account screening from the moire bands left outside the model space. Based on an atomic-scale tight-binding model, we develop a numerically tractable approximation to the polarization function and study its behavior for different twist angles. We find that the polarization has three different momentum regimes. For small momenta, the polarization is quadratic, leading to a linear dielectric function expected for a two-dimensional dielectric material. For large momenta, the polarization becomes independent of the twist angle and approaches that of uncoupled graphene layers. In the intermediate momentum regime, the dependence on the twist-angle is strong. Close to the largest magic angle the dielectric function peaks at a momentum of $~1/(4 \: nm)$ attaining values of 25, meaning very strong screening at intermediate distances. We also calculate the effective screened Coulomb interaction in real space and give estimates for the on-site and extended interaction terms for the recently developed hexagonal-lattice model. For free-standing TBG the effective interaction decays slower than $1/r$ at intermediate distances $r$, while it remains essentially unscreened at large enough $r$.

preprint2016arXiv

Accessing many-body localized states through the Generalized Gibbs Ensemble

We show how the thermodynamic properties of large many-body localized systems can be studied using quantum Monte Carlo simulations. To this end we devise a heuristic way of constructing local integrals of motion of very high quality, which are added to the Hamiltonian in conjunction with Lagrange multipliers. The ground state simulation of the shifted Hamiltonian corresponds to a high-energy state of the original Hamiltonian in case of exactly known local integrals of motion. We can show that the inevitable mixing between eigenstates as a consequence of non-perfect integrals of motion is weak enough such that the characteristics of many-body localized systems are not averaged out in our approach, unlike the standard ensembles of statistical mechanics. Our method paves the way to study higher dimensions and indicates that a full many-body localized phase in 2d, where (nearly) all eigenstates are localized, is likely to exist.

preprint2016arXiv

Bosonic self-energy functional theory

We derive the self-energy functional theory for bosonic lattice systems with broken $U(1)$ symmetry by parametrizing the bosonic Baym-Kadanoff effective action in terms of one- and two-point self-energies. The formalism goes beyond other approximate methods such as the pseudoparticle variational cluster approximation, the cluster composite boson mapping, and the Bogoliubov+U theory. It simplifies to bosonic dynamical-mean field theory when constraining to local fields, whereas when neglecting kinetic contributions of non-condensed bosons it reduces to the static mean-field approximation. To benchmark the theory we study the Bose-Hubbard model on the two- and three-dimensional cubic lattice, comparing with exact results from path integral quantum Monte Carlo. We also study the frustrated square lattice with next-nearest neighbor hopping, which is beyond the reach of Monte Carlo simulations. A reference system comprising a single bosonic state, corresponding to three variational parameters, is sufficient to quantitatively describe phase-boundaries, and thermodynamical observables, while qualitatively capturing the spectral functions, as well as the enhancement of kinetic fluctuations in the frustrated case. On the basis of these findings we propose self-energy functional theory as the omnibus framework for treating bosonic lattice models, in particular, in cases where path integral quantum Monte Carlo methods suffer from severe sign problems (e.g. in the presence of non-trivial gauge fields or frustration). Self-energy functional theory enables the construction of diagrammatically sound approximations that are quantitatively precise and controlled in the number of optimization parameters, but nevertheless remain computable by modest means.

preprint2016arXiv

Fulde-Ferrell-Larkin-Ovchinnikov pairing as leading instability on the square lattice

We study attractively interacting spin-1/2 fermions on the square lattice subject to a spin population imbalance. Using unbiased diagrammatic Monte Carlo simulations we find an extended region in the parameter space where the Fermi liquid is unstable towards formation of Cooper pairs with non-zero center-of-mass momentum, known as the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. In contrast to earlier mean-field and quasi-classical studies we provide quantitative and well-controlled predictions on the existence and location of the relevant Fermi-liquid instabilities. The highest temperature where the FFLO instability can be observed is about half of the superfluid transition temperature in the unpolarized system.

preprint2016arXiv

Grassmannization of classical models

Applying Feynman diagrammatics to non-fermionic strongly correlated models with local constraints might seem generically impossible for two separate reasons: (i) the necessity to have a Gaussian (non-interacting) limit on top of which the perturbative diagrammatic expansion is generated by Wick's theorem, and (ii) the Dyson's collapse argument implying that the expansion in powers of coupling constant is divergent. We show that for arbitrary classical lattice models both problems can be solved/circumvented by reformulating the high-temperature expansion (more generally, any discrete representation of the model) in terms of Grassmann integrals. Discrete variables residing on either links, plaquettes, or sites of the lattice are associated with the Grassmann variables in such a way that the partition function (and correlations) of the original system and its Grassmann-field counterpart are identical. The expansion of the latter around its Gaussian point generates Feynman diagrams. A proof-of-principle implementation is presented for the classical 2D Ising model. Our work paves the way for studying lattice gauge theories by treating bosonic and fermionic degrees of freedom on equal footing.

preprint2016arXiv

Ground state phase diagram of Gaussian-Core bosons in two dimensions

The ground state of a two-dimensional (2D) system of Bose particles of spin zero, interacting via a repulsive Gaussian-Core potential, has been investigated by means of Quantum Monte Carlo simulations. The quantum phase diagram is qualitatively identical to that of 2D Yukawa bosons. While the system is a fluid at all densities for weak coupling, in the strong coupling regime it transitions upon compression from a low density superfluid to a crystal, and then into a reentrant superfluid phase. No evidence of a (supersolid) cluster crystal phase is seen.

preprint2016arXiv

Mean-field phase diagram of the Bose-Fermi Hubbard model

We analyze the ground state properties of Bose-Fermi mixtures using a mean-field treatment of the boson-fermion interaction on a simple cubic lattice. In the deep superfluid limit of the bosonic sector and the BCS regime of the fermion sector, we derive BCS-type equations to determine the phase diagram of the system. We find a competition between a charge density wave and a superconducting phase. In the opposite limit, we study the Mott insulator to superfluid transition of the bosonic sector in the presence of a staggered density-induced alternating potential provided by the fermions, and determine the mean-field transition line. In the two-superfluid phase of the mixture we restrict to nearest-neighbor induced interactions between the fermions and consider the extended Hubbard model. We perform a mean-field analysis of the critical temperature for the formation of boson-assisted $s$-, extended $s^-$-, $d$-, and $p$-wave pairs at fermionic half filling. We compare our results with a recent dynamical mean-field study [Anders \emph{et al.} Phys. Rev. Lett. {\bf 109}, 206401].

preprint2016arXiv

Superfluid--Insulator Transition in Strongly Disordered One-dimensional Systems

We present an asymptotically exact renormalization-group theory of the superfluid--insulator transition in one-dimensional disordered systems, with emphasis on an accurate description of the interplay between the Giamarchi--Schulz (instanton--anti-instanton) and weak-link (scratched-XY) criticalities. Combining the theory with extensive quantum Monte Carlo simulations allows us to shed new light on the ground-state phase diagram of the one-dimensional disordered Bose-Hubbard model at unit filling.

preprint2016arXiv

The classical mutual information in mean-field spin glass models

We investigate the classical Renyi entropy S_n and the associated mutual information I_n in the Sherrington-Kirkpatrick (S-K) model, which is the paradigm model of mean-field spin glasses. Using classical Monte Carlo simulations and analytical tools we investigate the S-K model on the n-sheets booklet. This is obtained by gluing together n independent copies of the model, and it is the main ingredient to construct the Renyi entanglement-related quantities. We find a glassy phase at low temperature, whereas at high temperature the model exhibits paramagnetic behavior, consistent with the regular S-K model. The temperature of the paramagnetic-glassy transition depends non-trivially on the geometry of the booklet. At high-temperatures we provide the exact solution of the model by exploiting the replica symmetry. This is the permutation symmetry among the fictitious replicas that are used to perform disorder averages (via the replica trick). In the glassy phase the replica symmetry has to be broken. Using a generalization of the Parisi solution, we provide analytical results for S_n and I_n, and for standard thermodynamic quantities. Both S_n and I_n exhibit a volume law in the whole phase diagram. We characterize the behavior of the corresponding densities S_n/N, I_n/N, in the thermodynamic limit. Interestingly, at the critical point the mutual information does not exhibit any crossing for different system sizes, in contrast with local spin models.

preprint2015arXiv

Diagrammatic Monte Carlo study of mass-imbalanced Fermi-polaron system

We apply the diagrammatic Monte Carlo approach to three-dimensional Fermi-polaron systems with mass-imbalance, where an impurity interacts resonantly with a noninteracting Fermi sea whose atoms have a different mass. This method allows to go beyond frequently used variational techniques by stochastically summing all relevant impurity Feynman diagrams up to a maximum expansion order limited by the sign problem. Polaron energy and quasiparticle residue can be accurately determined over a broad range of impurity masses. Furthermore, the spectral function of an imbalanced polaron demonstrates the stability of the quasiparticle and allows to locate in addition also the repulsive polaron as an excited state. The quantitative exactness of two-particle-hole wave-functions is investigated, resulting in a relative lowering of polaronic energies in the mass-imbalance phase diagram. Tan's contact coefficient for the mass-balanced polaron system is found in good agreement with variational methods. Mass-imbalanced systems can be studied experimentally by ultracold atom mixtures like $^6$Li-$^{40}$K.

preprint2015arXiv

Exotic Superconductivity Through Bosons in a Dynamical Cluster Approximation

We study the instabilities towards (exotic) superconductivity of mixtures of spin-$1/2$ fermions coupled to scalar bosons on a two-dimensional square lattice with the Dynamical-Cluster-Approximation (DCA) using a numerically exact continuous-time Monte-Carlo solver. The Bogoliubov bosons provide an effective phononic bath for the fermions and induce a non-local retarded interaction between the fermions, which can lead to (exotic) superconductivity. Because of the sign problem the biggest clusters we can study are limited to $2 \times 2$ in size, but this nevertheless allows us to study the pairing instablilities, and their possible divergence, in the $s$- and $d$ -wave channels as well as the competition with antiferromagnetic fluctuations. At fermionic half-filling we find that $d$-wave is stable when the mediated interaction by the bosons is of the same order as the bare fermionic repulsion. Its critical temperature can be made as high as the maximum one for $s$-wave, which opens perspectives for its detection in a cold atom experiment.

preprint2015arXiv

Identifying a bath-induced Bose liquid in interacting spin-boson models

We study the ground state phase diagram of a one-dimensional hard-core bosonic model with nearest-neighbor interactions (XXZ model) where every site is coupled Ohmically to an independent but identical reservoir, hereby generalizing spin-boson models to interacting spin-boson systems. We show that a bath-induced Bose metal phase can occur in the ground state phase diagram away from half filling. This phase is compressible, gapless, and conducting but not superfluid. At haf-filling, only a Luttinger liquid and a charge density wave are found. The phase transition between them is of Kosterlitz-Thouless type where the Luttinger parameter takes a non-universal value.The applied quantum Monte Carlo method can be used for all open bosonic and unfrustrated spin systems, regardless of their dimension, filling factor and spectrum of the dissipation as long as the quantum system couples to the bath via the density operators.

preprint2015arXiv

p-wave Superfluidity by Spin-Nematic Fermi Surface Deformation

We study attractively interacting fermions on a square lattice with dispersion relations exhibiting strong spin-dependent anisotropy. The resulting Fermi surface mismatch suppresses the s-wave BCS-type instability, clearing the way for unconventional types of order. Unbiased sampling of the Feynman diagrammatic series using Diagrammatic Monte Carlo methods reveals a rich phase diagram in the regime of intermediate coupling strength. Instead of a proposed Cooper-pair Bose metal phase [A. E. Feiguin and M. P. A. Fisher, Phys. Rev. Lett. 103, 025303 (2009)] we find an incommensurate density wave at strong anisotropy and two different p-wave superfluid states with unconventional symmetry at intermediate anisotropy.

preprint2015arXiv

Spin and charge dynamics of a quasi-one-dimensional antiferromagnetic metal

We use quantum Monte Carlo simulations to study a finite-temperature dimensional-crossover-driven evolution of spin and charge dynamics in weakly coupled Hubbard chains with a half-filled band. The low-temperature behavior of the charge gap indicates a crossover between two distinct energy scales: a high-energy one-dimensional (1D) Mott gap due to the umklapp process and a low-energy gap which stems from long-range antiferromagnetic (AF) fluctuations. Away from the 1D regime and at temperature scales above the charge gap, the emergence of a zero-frequency Drude-like feature in the interchain optical conductivity $σ_{\perp}(ω)$ implies the onset of a higher-dimensional metal. In this metallic phase, enhanced quasiparticle scattering off finite-range AF fluctuations results in incoherent single-particle dynamics. The coupling between spin and charge fluctuations is also seen in the spin dynamical structure factor $S({\pmb q},ω)$ displaying damped spin excitations (paramagnons) close to the AF wave-vector ${\pmb q}=(π,π)$ and particle-hole continua near 1D momentum transfers spanning quasiparticles at the Fermi surface. We relate our results to the charge deconfinement in quasi-1D organic Bechgaard-Fabre salts.

preprint2015arXiv

The Massive Goldstone (Higgs) mode in two-dimensional ultracold atomic lattice systems

We discuss how to reveal the massive Goldstone mode, often referred to as the Higgs amplitude mode, near the Superfluid-to-Insulator quantum critical point (QCP) in a system of two-dimensional ultracold bosonic atoms in optical lattices. The spectral function of the amplitude response is obtained by analytic continuation of the kinetic energy correlation function calculated by Monte Carlo methods. Our results enable a direct comparison with the recent experiment [M. Endres, T. Fukuhara, D. Pekker, M. Cheneau, P. Schauß, C. Gross, E. Demler, S. Kuhr, and I. Bloch, Nature 487, 454-458 (2012)], and demonstrate a good agreement for temperature shifts induced by lattice modulation. Based on our numerical analysis, we formulate the necessary conditions in terms of homogeneity, detuning from the QCP and temperature in order to reveal the massive Goldstone resonance peak in spectral functions experimentally. We also propose to apply a local modulation at the trap center to overcome the inhomogeneous broadening caused by the parabolic trap confinement.

preprint2015arXiv

Thermodynamics of the Bose-Hubbard model in a Bogoliubov+U theory

We derive the Bogoliubov+U formalism to study the thermodynamical properties of the Bose-Hubbard model. The framework can be viewed as the zero-frequency limit of bosonic dynamical mean-field theory (B-DMFT), but equally well as an extension of the mean-field decoupling approximation in which pair creation and annihilation of depleted particles is taken into account. The self-energy on the impurity site is treated variationally, minimizing the grand potential. The theory containing just three parameters that are determined self-consistently reproduces the T=0 phase diagrams of the three-dimensional and two-dimensional Bose-Hubbard model with an accuracy of 1% or better. The superfluid to normal transition at finite temperature is also reproduced well and only slightly less accurately than in B-DMFT.

preprint2014arXiv

Asymptotically Exact Scenario of Strong-Disorder Criticality in One-Dimensional Superfluids

We present a controlled rare-weak-link theory of the superfluid-to-Bose/Mott glass transition in one-dimensional disordered systems. The transition has Kosterlitz-Thouless critical properties but may occur at an arbitrary large value of the Luttinger parameter $K$. In contrast to the scenario by Altman {\it et al.} [Phys. Rev. B {\bf 81}, 174528 (2010)], the hydrodynamic description is valid under the correlation radius and defines criticality via the renormalization of microscopically weak links, along the lines of Kane and Fisher [Phys. Rev. Lett. {\bf 68}, 1220 (1992)]. The hallmark of the theory is the relation $K^{(c)}=1/ζ$ between the critical value of the Luttinger parameter at macroscopic scales and the microscopic (irrenormalizable) exponent $ζ$ describing the scaling $\propto 1/N^{1-ζ}$ for the strength of the weakest link among the $N/L \gg 1$ disorder realizations in a system of fixed mesoscopic size $L$.

preprint2014arXiv

Diagrammatic Monte Carlo study of quasi-two-dimensional Fermi-polarons

We apply a diagrammatic Monte Carlo method to the problem of an impurity interacting resonantly with a homogeneous Fermi bath for a quasi-two-dimensional setup. Notwithstanding the series divergence, we can show numerically that the three particle-hole diagrammatic contributions are not contributing significantly to the final answer, thus demonstrating a nearly perfect destructive interference of contributions in subspaces with higher-order particle-hole lines. Consequently, for strong enough confinement in the third direction, the transition between the polaron and the molecule ground state is found to be in good agreement with the pure two-dimensional case and agrees very well with the one found by the wave-function approach in the two-particle-hole subspace.

preprint2014arXiv

Ground state phase diagram of the 2d Bose-Hubbard model with anisotropic hopping

We compute the ground state phase diagram of the 2d Bose-Hubbard model with anisotropic hopping using quantum Monte Carlo simulations, connecting the 1d to the 2d system. We find that the tip of the lobe lies on a curve controlled by the 1d limit over the full anisotropy range while the universality class is always the same as in the isotropic 2d system. This behavior can be derived analytically from the lowest RG equations and has a form typical for the underlying Kosterlitz-Thouless transition in 1d. We also compute the phase boundary of the Mott lobe for strong anisotropy and compare it to the 1d system. Our calculations shed light on recent cold gas experiments monitoring the dynamics of an expanding cloud.

preprint2014arXiv

Quantum Plasticity and Supersolid Response in Helium-4

We argue that the three key phenomena recently observed in solid $^4$He---mass supertransport, anomalous isochoric compressibility (syringe effect), and giant plasticity---are closely linked to each other through the physics of an interconnected network of tilted quantum-rough dislocations. As immediate implications of this connection several predictions follow: In the absence of $^3$He impurities, the syringe effect and giant plasticity persist down to $T=0$; the dynamical low-frequency syringe and giant-plasticity responses are dispersionless; and similarly to giant plasticity but without direct relationship to the supertransport along the dislocation cores, $^3$He impurities should suppress the syringe effect partially or completely at appropriately low temperatures.

preprint2014arXiv

Spontaneous emissions and thermalization of cold bosons in optical lattices

We study the thermalization of excitations generated by spontaneous emission events for cold bosons in an optical lattice. Computing the dynamics described by the many-body master equation, we characterize equilibration timescales in different parameter regimes. For simple observables, we find regimes in which the system relaxes rapidly to values in agreement with a thermal distribution, and others where thermalization does not occur on typical experimental timescales. Because spontaneous emissions lead effectively to a local quantum quench, this behavior is strongly dependent on the low-energy spectrum of the Hamiltonian, and undergoes a qualitative change at the Mott Insulator-superfluid transition point. These results have important implications for the understanding of thermalization after localized quenches in isolated quantum gases, as well as the characterization of heating in experiments.

preprint2014arXiv

Thermalization of strongly interacting bosons after spontaneous emissions in optical lattices

We study the out-of-equilibrium dynamics of bosonic atoms in a 1D optical lattice, after the ground-state is excited by a single spontaneous emission event, i.e. after an absorption and re-emission of a lattice photon. This is an important fundamental source of decoherence for current experiments, and understanding the resulting dynamics and changes in the many-body state is important for controlling heating in quantum simulators. Previously it was found that in the superfluid regime, simple observables relax to values that can be described by a thermal distribution on experimental time-scales, and that this breaks down for strong interactions (in the Mott insulator regime). Here we expand on this result, investigating the relaxation of the momentum distribution as a function of time, and discussing the relationship to eigenstate thermalization. For the strongly interacting limit, we provide an analytical analysis for the behavior of the system, based on an effective low-energy Hamiltonian in which the dynamics can be understood based on correlated doublon-holon pairs.

preprint2014arXiv

Topological quantum phases of ${^4}$He confined to nanoporous materials

The ground state of $^4$He confined in a system with the topology of a cylinder can display properties of a solid, superfluid and liquid crystal. This phase, which we call compactified supersolid (CSS), originates from wrapping the basal planes of the bulk hcp solid into concentric cylindrical shells, with several central shells exhibiting superfluidity along the axial direction. Its main feature is the presence of a topological defect which can be viewed as a disclination with Frank index $n=1$ observed in liquid crystals, and which, in addition, has a superfluid core. The CSS as well as its transition to an insulating compactified solid with a very wide hysteresis loop are found by ab initio Monte Carlo simulations. A simple analytical model captures qualitatively correctly the main property of the CSS -- a gradual decrease of the superfluid response with increasing pressure.

preprint2014arXiv

Universal Conductivity in a Two-Dimensional Superfluid-to-Insulator Quantum Critical System

We compute the universal conductivity of the (2+1)-dimensional XY universality class, which is realized for a superfluid-to-Mott insulator quantum phase transition at constant density. Based on large-scale Monte Carlo simulations of the classical (2+1)-dimensional $J$-current model and the two-dimensional Bose-Hubbard model, we can precisely determine the conductivity on the quantum critical plateau, $σ(\infty)=0.359(4)σ_Q$ with $σ_Q$ the conductivity quantum. The universal conductivity curve is the standard example with the lowest number of components where the bottoms-up AdS/CFT correspondence from string theory can be tested and made to use [R. C. Myers, S. Sachdev, and A. Singh, Phys. Rev. D83, 066017 (2011)]. For the first time, the shape of the $σ(iω_n)- σ(\infty)$ function in the Matsubara representation is accurate enough for a conclusive comparison and establishes the particle-like nature of charge transport. We find that the holographic gauge/gravity duality theory for transport properties can be made compatible with the data if temperature of the horizon of the black brane is different from the temperature of the conformal field theory. The requirements for measuring the universal conductivity in a cold gas experiment are also determined by our calculation.

preprint2013arXiv

A review of Monte Carlo simulations for the Bose-Hubbard model with diagonal disorder

We review the physics of the Bose-Hubbard model with disorder in the chemical potential focusing on recently published analytical arguments in combination with quantum Monte Carlo simulations. Apart from the superfluid and Mott insulator phases that can occur in this system without disorder, disorder allows for an additional phase, called the Bose glass phase. The topology of the phase diagram is subject to strong theorems proving that the Bose Glass phase must intervene between the superfluid and the Mott insulator and implying a Griffiths transition between the Mott insulator and the Bose glass. The full phase diagrams in 3d and 2d are discussed, and we zoom in on the insensitivity of the transition line between the superfluid and the Bose glass in the close vicinity of the tip of the Mott insulator lobe. We briefly comment on the established and remaining questions in the 1d case, and give a short overview of numerical work on related models.

preprint2013arXiv

Demixing in symmetric supersolid mixtures

The droplet crystal phase of a symmetric binary mixture of Rydberg-blockaded dipolar Bose gases is studied by computer simulation. At high temperature each droplet comprises on average equal numbers of particles of either component, but the two components demix below the supersolid transition temperature, i.e., droplets mostly consist of particles of one component. Droplets consisting of the same component will also favor clustering. Demixing is driven by quantum tunnelling of particles across droplets over the system, and does not take place in a non-superfluid crystal. This effect should be easily detectable in a cold gas experiment.

preprint2013arXiv

Phase diagram of $^{4}$He on graphene

The low temperature phase diagram of $^4$He adsorbed on a single graphene sheet is studied by computer simulation of a system comprising nearly thousand helium atoms. In the first layer, two commensurate solid phases are observed, with fillings 1/3 and 7/16 respectively, separated by a domain wall phase, as well as an incommensurate crystal at higher coverage. No evidence of a thermodynamically stable superfliuid phase is found for the first adlayer. Second layer promotion occurs at a coverage of 0.111(4) $Å^{-2}$. In the second layer two phases are observed, namely a superfluid and an incommensurate solid, with no commensurate solid intervening between these two phases. The computed phase diagram closely resembles that predicted for helium on graphite.

preprint2013arXiv

Single-site- and single-atom-resolved measurement of correlation functions

Correlation functions play an important role for the theoretical and experimental characterization of many-body systems. In solid-state systems, they are usually determined through scattering experiments whereas in cold-gases systems, time-of-flight and in-situ absorption imaging are the standard observation techniques. However, none of these methods allow the in-situ detection of spatially resolved correlation functions at the single-particle level. Here we give a more detailed account of recent advances in the detection of correlation functions using in-situ fluorescence imaging of ultracold bosonic atoms in an optical lattice. This method yields single-site and single-atom-resolved images of the lattice gas in a single experimental run, thus gaining direct access to fluctuations in the many-body system. As a consequence, the detection of correlation functions between an arbitrary set of lattice sites is possible. This enables not only the detection of two-site correlation functions but also the evaluation of non-local correlations, which originate from an extended region of the system and are used for the characterization of quantum phases that do not possess (quasi-)long-range order in the traditional sense.

preprint2013arXiv

Universal Properties of the Higgs Resonance in (2+1)-Dimensional U(1) Critical Systems

We present spectral functions for the magnitude squared of the order parameter in the scaling limit of the two-dimensional superfluid to Mott insulator quantum phase transition at constant density, which has emergent particle-hole symmetry and Lorentz invariance. The universal functions for the superfluid, Mott insulator, and normal liquid phases reveal a low-frequency resonance which is relatively sharp and is followed by a damped oscillation (in the first two phases only) before saturating to the quantum critical plateau. The counterintuitive resonance feature in the insulating and normal phases calls for deeper understanding of collective modes in the strongly coupled (2+1)-dimensional relativistic field theory. Our results are derived from analytically continued correlation functions obtained from path-integral Monte Carlo simulations of the Bose-Hubbard model.

preprint2012arXiv

From the Cooper problem to canted supersolids in Bose-Fermi mixtures

We calculate the phase diagram of the Bose-Fermi Hubbard model on the 3d cubic lattice at fermionic half filling and bosonic unit filling by means of single-site dynamical mean-field theory. For fast bosons, this is equivalent to the Cooper problem in which the bosons can induce s-wave pairing between the fermions. We also find miscible superfluid and canted supersolid phases depending on the interspecies coupling strength. In contrast, slow bosons favor fermionic charge density wave structures for attractive fermionic interactions. These competing instabilities lead to a rich phase diagram within reach of cold gas experiments.

preprint2012arXiv

Recent developments in Quantum Monte-Carlo simulations with applications for cold gases

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.

preprint2011arXiv

Dynamical mean-field theory for bosons

We discuss the recently developed bosonic dynamical mean-field (B-DMFT) framework, which maps a bosonic lattice model onto the selfconsistent 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) by using an effective medium approach based on adding a one-loop correction to the selfconsistently defined condensate. To solve the impurity problem, we use a continuous-time Monte Carlo algorithm based on a 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 3d cubic and 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.

preprint2011arXiv

Localization of Toric Code Defects

We explore the possibility of passive error correction in the toric code model. We first show that even coherent dynamics, stemming from spin interactions or the coupling to an external magnetic field, lead to logical errors. We then argue that Anderson localization of the defects, arising from unavoidable fluctuations of the coupling constants, provides a remedy. This protection is demonstrated using general analytical arguments that are complemented with numerical results which demonstrate that self-correcting memory can in principle be achieved in the limit of a nonzero density of identical defects.

preprint2010arXiv

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.

preprint2010arXiv

Dynamical mean field solution of the Bose-Hubbard model

We present the effective action and self-consistency equations for the bosonic dynamical mean field (B-DMFT) approximation to the bosonic Hubbard model and show that it provides remarkably accurate phase diagrams and correlation functions. To solve the bosonic dynamical mean field equations we use a continuous-time Monte Carlo method for bosonic impurity models based on a diagrammatic expansion in the hybridization and condensate coupling. This method is readily generalized to bosonic mixtures, spinful bosons, and Bose-Fermi mixtures.

preprint2010arXiv

Measuring the equation of state of trapped ultracold bosonic systems in an optical lattice with in-situ density imaging

We analyze quantitatively how imaging techniques with single-site resolution allow to measure thermodynamical properties that cannot be inferred from time-of-light images for the trapped Bose-Hubbard model. If the normal state extends over a sufficiently large range, the chemical potential and the temperature can be extracted from a single shot, provided the sample is in thermodynamic equilibrium. When the normal state is too narrow, temperature is low but can still be extracted using the fluctuation-dissipation theorem over the entire trap range as long as the local density approximation remains valid, as was recently suggested by Qi Zhou and Tin-Lun Ho [arXiv:0908.3015]. However, for typical present-day experiments, the number of samples needed is of the order of 1000 in order to get the temperature at least $10 \%$ accurate, but it is possible to reduce the variance by 2 orders of magnitude if the density-density correlation length is short, which is the case for the Bose-Hubbard model. Our results provide further evidence that cold gases in an optical lattices can be viewed as quantum analog computers.

preprint2010arXiv

Probing the Superfluid to Mott Insulator Transition at the Single Atom Level

Quantum gases in optical lattices offer an opportunity to experimentally realize and explore condensed matter models in a clean, tunable system. We investigate the Bose-Hubbard model on a microscopic level using single atom-single lattice site imaging; our technique enables space- and time-resolved characterization of the number statistics across the superfluid-Mott insulator quantum phase transition. Site-resolved probing of fluctuations provides us with a sensitive local thermometer, allows us to identify microscopic heterostructures of low entropy Mott domains, and enables us to measure local quantum dynamics, revealing surprisingly fast transition timescales. Our results may serve as a benchmark for theoretical studies of quantum dynamics, and may guide the engineering of low entropy phases in a lattice.

preprint2010arXiv

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.77kB which is about twice as large as the entropy at the Neel transition. The corresponding temperature depends on the atom number and for small fillings reaches values on the order of the tunneling energy.

preprint2010arXiv

Regularization of Diagrammatic Series with Zero Convergence Radius

The divergence of perturbative expansions for the vast majority of macroscopic systems, which follows from Dyson's collapse argument, prevents Feynman's diagrammatic technique from being directly used for controllable studies of strongly interacting systems. We show how the problem of divergence can be solved by replacing the original model with a convergent sequence of successive approximations which have a convergent perturbative series. As a prototypical model, we consider the zero-dimensional $\vert ψ\vert^4$ theory.

preprint2009arXiv

Disorder-induced superfluidity

We use quantum Monte Carlo simulations to study the phase diagram of hard-core bosons with short-ranged {\it attractive} interactions, in the presence of uniform diagonal disorder. It is shown that moderate disorder stabilizes a glassy superfluid phase in a range of values of the attractive interaction for which the system is a Mott insulator, in the absence of disorder. A transition to an insulating Bose glass phase occurs as the strength of the disorder or interactions increases.

preprint2008arXiv

Influence of the trap shape on the superfluid-Mott transition in ultracold atomic gases

The coexistence of superfluid and Mott insulator, due to the quadratic confinement potential in current optical lattice experiments, makes the accurate detection of the superfluid-Mott transition difficult. Studying alternative trapping potentials which are experimentally realizable and have a flatter center, we find that the transition can be better resolved, but at the cost of a more difficult tuning of the particle filling. When mapping out the phase diagram using local probes and the local density approximation we find that the smoother gradient of the parabolic trap is advantageous.

preprint2008arXiv

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 that current experiments remain in the quantum regime for all considered lattice depths with low temperatures and minimal heating.