Role of Bose statistics in crystallization and quantum jamming.

The indistinguishability of particles is a crucial factor destabilizing crystalline order in Bose systems. We describe this effect in terms of damped quasiparticle modes and in the dual language of Feynman paths, and illustrate it by first-principles simulations of dipolar bosons and bulk condensed 4He. The first major implication is that, contrary to conventional wisdom, zero-point motion alone cannot prevent 4He crystallization at near zero pressure. Second, Bose statistics leads to quantum jamming at finite temperature, dramatically enhancing the metastability of superfluid glasses. Only studies of indistinguishable particles can reliably address these issues.

Higgs mode in a two-dimensional superfluid.

We present solid evidence for the existence of a well-defined Higgs amplitude mode in two-dimensional relativistic field theories based on analytically continued results from quantum Monte Carlo simulations of the Bose-Hubbard model in the vicinity of the superfluid-Mott insulator quantum critical point, featuring emergent particle-hole symmetry and Lorentz invariance. The Higgs boson, seen as a well-defined low-frequency resonance in the spectral density, is quickly pushed to high energies in the superfluid phase and disappears by merging with the broad secondary peak at the characteristic interaction scale. Simulations of a trapped system of ultracold (87)Rb atoms demonstrate that the low-frequency resonance is lost for typical experimental parameters, while the characteristic frequency for the onset of a strong response is preserved.

Observation of correlated particle-hole pairs and string order in low-dimensional Mott insulators.

Quantum phases of matter are characterized by the underlying correlations of the many-body system. Although this is typically captured by a local order parameter, it has been shown that a broad class of many-body systems possesses a hidden nonlocal order. In the case of bosonic Mott insulators, the ground state properties are governed by quantum fluctuations in the form of correlated particle-hole pairs that lead to the emergence of a nonlocal string order in one dimension. By using high-resolution imaging of low-dimensional quantum gases in an optical lattice, we directly detect these pairs with single-site and single-particle sensitivity and observe string order in the one-dimensional case.

Criticality in trapped atomic systems.

We discuss generic limits posed by the trap in atomic systems on the accurate determination of critical parameters for second-order phase transitions, from which we deduce optimal protocols to extract them. We show that under current experimental conditions the in situ density profiles are barely suitable for an accurate study of critical points in the strongly correlated regime. Contrary to recent claims, the proper analysis of time-of-fight images yields critical parameters accurately.

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 used single atom-single lattice site imaging to investigate the Bose-Hubbard model on a microscopic level. 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 time scales. 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.

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.

Supersolid phase with cold polar molecules on a triangular lattice.

We study a system of heteronuclear molecules on a triangular lattice and analyze the potential of this system for the experimental realization of a supersolid phase. The ground state phase diagram contains superfluid, solid, and supersolid phases. At finite temperatures and strong interactions there is an additional emulsion region, in contrast with similar models with short-range interactions. We derive the maximal critical temperature T{c} and the corresponding entropy S/N=0.04(1) for supersolidity and find feasible experimental conditions for its realization.

Absence of a direct superfluid to mott insulator transition in disordered bose systems.

We prove the absence of a direct quantum phase transition between a superfluid and a Mott insulator in a bosonic system with generic, bounded disorder. We also prove the compressibility of the system on the superfluid-insulator critical line and in its neighborhood. These conclusions follow from a general theorem of inclusions, which states that for any transition in a disordered system, one can always find rare regions of the competing phase on either side of the transition line. Quantum Monte Carlo simulations for the disordered Bose-Hubbard model show an even stronger result, important for the nature of the Mott insulator to Bose glass phase transition: the critical disorder bound Delta(c) corresponding to the onset of disorder-induced superfluidity, satisfies the relation Delta(c)>Eg/2, with Eg/2 the half-width of the Mott gap in the pure system.

Underlying mechanism for the giant isochoric compressibility of solid (4)He: superclimb of dislocations[ corrected].

In the experiment on superfluid transport in solid 4He [Phys. Rev. Lett. 100, 235301 (2008)], Ray and Hallock observed an anomalously large isochoric compressibility: the supersolid samples demonstrated a significant and apparently spatially uniform response of density and pressure to chemical potential, applied locally through Vycor "electrodes." We propose that the effect is due to superclimb: edge dislocations can climb because of mass transport along superfluid cores. We corroborate the scenario by ab initio simulations of an edge dislocation in solid 4He at T = 0.5 K. We argue that at low temperature the effect must be suppressed due to a crossover to the smooth dislocation.

Discerning incompressible and compressible phases of cold atoms in optical lattices.

Experiments with cold atoms trapped in optical lattices offer the potential to realize a variety of novel phases but suffer from severe spatial inhomogeneity that can obscure signatures of new phases of matter and phase boundaries. We use a high temperature series expansion to show that compressibility in the core of a trapped Fermi-Hubbard system is related to measurements of changes in double occupancy. This core compressibility filters out edge effects, offering a direct probe of compressibility independent of inhomogeneity. A comparison with experiments is made.

Binding of a 3He impurity to a screw dislocation in solid 4He.

Using first-principles simulations for the probability density of finding a 3He atom in the vicinity of the screw dislocation in solid 4He, we determine the binding energy to the dislocation nucleus E(B)=0.8+/-0.1 K and the density of localized states at larger distances. The specific heat due to 3He features a peak similar to the one observed in recent experiments, and our model can also account for the observed increase in shear modulus at low temperature. We further discuss the role of 3He in the picture of superfluid defects.

Expansion of a quantum gas released from an optical lattice.

We analyze the interference pattern produced by ultracold atoms released from an optical lattice, commonly interpreted as the momentum distributions of the trapped quantum gas. We show that for finite times of flight the resulting density distribution can, however, be significantly altered, similar to a near-field diffraction regime in optics. We illustrate our findings with a simple model and realistic quantum Monte Carlo simulations for bosonic atoms and compare the latter to experiments.

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.

Luttinger liquid in the core of a screw dislocation in helium-4.

On the basis of first-principles Monte Carlo simulations we find that the screw dislocation along the hexagonal axis of an hcp 4He crystal features a superfluid (at T-->0) core. This is the first example of a regular quasi-one-dimensional supersolid--the phase featuring both translational and superfluid orders, and one of the cleanest cases of a Luttinger-liquid system. In contrast, the same type of screw dislocation in solid H2 is insulating.

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.

Fate of vacancy-induced supersolidity in 4He.

The supersolid state of matter, exhibiting nondissipative flow in solids, has been elusive for 35 years. The recent discovery of a nonclassical moment of inertia in solid 4He by Kim and Chan provided the first experimental evidence, although the interpretation in terms of supersolidity of the ideal crystal phase remains a subject to debate. Using quantum Monte Carlo methods we investigate the long-standing question of vacancy-induced superflow and find that vacancies in a 4He crystal phase separate instead of forming a supersolid. On the other hand, nonequilibrium vacancies relaxing on defects of polycrystalline samples could provide an explanation for the experimental observations.

Quantum Monte Carlo simulation in the canonical ensemble at finite temperature.

A quantum Monte Carlo method with a nonlocal update scheme is presented. The method is based on a path-integral decomposition and a worm operator which is local in imaginary time. It generates states with a fixed number of particles and respects other exact symmetries. Observables like the equal-time Green's function can be evaluated in an efficient way. To demonstrate the versatility of the method, results for the one-dimensional Bose-Hubbard model and a nuclear pairing model are presented. Within the context of the Bose-Hubbard model the efficiency of the algorithm is discussed.

Loop updates for quantum Monte Carlo simulations in the canonical ensemble.

We present a new nonlocal updating scheme for quantum Monte Carlo simulations, which conserves particle number and other symmetries. It allows exact symmetry projection and direct evaluation of the equal-time Green's function and other observables in the canonical ensemble. The method is applicable to a wide variety of systems. We show results for bosonic atoms in optical lattices, neutron pairs in atomic nuclei, and electron pairs in ultrasmall superconducting grains.

Ultracold atoms in one-dimensional optical lattices approaching the Tonks-Girardeau regime.

Recent experiments on ultracold atomic alkali gases in a one-dimensional optical lattice have demonstrated the transition from a gas of soft-core bosons to a Tonks-Girardeau gas in the hard-core limit, where one-dimensional bosons behave like fermions in many respects. We have studied the underlying many-body physics through numerical simulations which accommodate both the soft-core and hard-core limits in one single framework. We find that the Tonks-Girardeau gas is reached only at the strongest optical lattice potentials. Results for slightly higher densities, where the gas develops a Mott-like phase already at weaker optical lattice potentials, show that these Mott-like short-range correlations do not enhance the convergence to the hard-core limit.