Consequences of the Pauli exclusion principle for the Bose-Einstein condensation of atoms and excitons

The bosonic atoms used in present day experiments on Bose-Einstein condensation are made up of fermionic electrons and nucleons. In this Letter we demonstrate how the Pauli exclusion principle for these constituents puts an upper limit on the Bose-Einstein-condensed fraction. Detailed numerical results are presented for hydrogen atoms in a cubic volume and for excitons in semiconductors and semiconductor bilayer systems. The resulting condensate depletion scales differently from what one expects for bosons with a repulsive hard-core interaction. At high densities, Pauli exclusion results in significantly more condensate depletion. These results also shed a new light on the low condensed fraction in liquid helium II.Comment: 4 pages, 2 figures, revised version, now includes a direct comparison with hard-sphere QMC results, submitted to Phys. Rev. Let

Polynomial complexity despite the fermionic sign

It is commonly believed that in quantum Monte Carlo approaches to fermionic many- body problems, the infamous sign problem generically implies prohibitively large computational times for obtaining thermodynamic-limit quantities. We point out that for convergent Feynman diagrammatic series evaluated with the Monte Carlo algorithm of [Rossi, arXiv:1612.05184], the computational time increases only polynomially with the inverse error on thermodynamic-limit quantities

Diagrammatic Monte Carlo algorithm for the resonant Fermi gas

We provide a description of a diagrammatic Monte Carlo algorithm for the resonant Fermi gas in the normal phase. Details are given on diagrammatic framework, Monte Carlo moves, and incorporation of ultraviolet asymptotics. Apart from the self-consistent bold scheme, we also describe a non-self-consistent scheme, for which the ultraviolet treatment is more involved.Comment: Revised and extended versio

Schwinger-Dyson equations in large-N quantum field theories and nonlinear random processes

We propose a stochastic method for solving Schwinger-Dyson equations in large-N quantum field theories. Expectation values of single-trace operators are sampled by stationary probability distributions of the so-called nonlinear random processes. The set of all histories of such processes corresponds to the set of all planar diagrams in the perturbative expansions of the expectation values of singlet operators. We illustrate the method on the examples of the matrix-valued scalar field theory and the Weingarten model of random planar surfaces on the lattice. For theories with compact field variables, such as sigma-models or non-Abelian lattice gauge theories, the method does not converge in the physically most interesting weak-coupling limit. In this case one can absorb the divergences into a self-consistent redefinition of expansion parameters. Stochastic solution of the self-consistency conditions can be implemented as a "memory" of the random process, so that some parameters of the process are estimated from its previous history. We illustrate this idea on the example of two-dimensional O(N) sigma-model. Extension to non-Abelian lattice gauge theories is discussed.Comment: 16 pages RevTeX, 14 figures; v2: Algorithm for the Weingarten model corrected; v3: published versio

High-precision numerical solution of the Fermi polaron problem and large-order behavior of its diagrammatic series

We introduce a simple determinant diagrammatic Monte Carlo algorithm to compute the ground-state properties of a particle interacting with a Fermi sea through a zero-range interaction. The fermionic sign does not cause any fundamental problem when going to high diagram orders, and we reach order $N=30$. The data reveal that the diagrammatic series diverges exponentially as $(-1/R)^{N}$ with a radius of convergence $R<1$. Furthermore, on the polaron side of the polaron-dimeron transition, the value of $R$ is determined by a special class of three-body diagrams, corresponding to repeated scattering of the impurity between two particles of the Fermi sea. A power-counting argument explains why finite $R$ is possible for zero-range interactions in three dimensions. Resumming the divergent series through a conformal mapping yields the polaron energy with record accuracy

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 Feynman's diagrammatic series for the single-particle self-energy we can study moderate values of the on-site repulsion ($U/t \sim 4$) and temperatures down to $T/t=1/40$. We compare our results with high temperature series expansion and with single-site and cluster dynamical mean-field theory.Comment: 4 pages, 5 figures, stylistic change

Integrable models for asymmetric Fermi superfluids: Emergence of a new exotic pairing phase

We introduce an exactly-solvable model to study the competition between the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) and breached-pair superfluid in strongly interacting ultracold asymmetric Fermi gases. One can thus investigate homogeneous and inhomogeneous states on an equal footing and establish the quantum phase diagram. For certain values of the filling and the interaction strength, the model exhibits a new stable exotic pairing phase which combines an inhomogeneous state with an interior gap to pair-excitations. It is proven that this phase is the exact ground state in the strong coupling limit, while numerical examples demonstrate that also at finite interaction strength it can have lower energy than the breached-pair or LOFF states.Comment: Revised version accepted for publicatio

From Popov-Fedotov trick to universal fermionization

We show that Popov-Fedotov trick of mapping spin-1/2 lattice systems on two-component fermions with imaginary chemical potential readily generalizes to bosons with a fixed (but not limited) maximal site occupation number, as well as to fermionic Hamiltonians with various constraints on the site Fock states. In a general case, the mapping---fermionization---is on multi-component fermions with many-body non-Hermitian interactions. Additionally, the fermionization approach allows one to convert large many-body couplings into single-particle energies, rendering the diagrammatic series free of large expansion parameters; the latter is essential for the efficiency and convergence of the diagrammatic Monte Carlo method.Comment: 4 pages, no figures (v2 contains some improvements; the most important one is the generic complex chemical potential trick for spins/bosons

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 \psi \vert^4$ theory.Comment: 4 pages, 3 figure