We present a phase-space method for fermionic systems, which enables simulations of the dynamics and thermal equilibrium states of many-body quantum systems from first principles. As examples, we calculate finite-temperature correlation functions.
We introduce a new class of quantum Monte Carlo methods, based on a Gaussian quantum operator representation of fermionic states. The methods enable first-principles dynamical or equilibrium calculations in many-body Fermi systems, and, combined with the existing Gaussian representation for bosons, provide a unified method of simulating Bose-Fermi systems. As an application relevant to the Fermi sign problem, we calculate finite-temperature properties of the two dimensional Hubbard model and the dynamics in a simple model of coherent molecular dissociation.
We introduce a positive phase-space representation for fermions, using the most general possible multimode Gaussian operator basis. The representation generalizes previous bosonic quantum phase-space methods to Fermi systems. We derive equivalences between quantum and stochastic moments, as well as operator correspondences that map quantum operator evolution onto stochastic processes in phase space. The representation thus enables first-principles quantum dynamical or equilibrium calculations in many-body Fermi systems. Potential applications are to strongly interacting and correlated Fermi gases, including coherent behavior in open systems and nanostructures described by master equations. Examples of an ideal gas and the Hubbard model are given, as well as a generic open system, in order to illustrate these ideas.
Phase-space representations are of increasing importance as a viable and successful means to study exponentially complex quantum many-body systems from first principles. This paper traces the background of these methods, starting from the early work of Wigner, Glauber and Sudarshan. We focus on modern phase-space approaches using non-classical phase-space representations. These lead to the Gaussian representation, which unifies bosonic and fermionic phase-space. Examples treated include quantum solitons in optical fibers, colliding Bose-Einstein condensates, and strongly correlated fermions on lattices.
We demonstrate that the quantum dynamics of a many-body Fermi-Bose system can be simulated using a Gaussian phase-space representation method. In particular, we consider the application of the mixed fermion-boson model to ultracold quantum gases and simulate the dynamics of dissociation of a Bose-Einstein condensate of bosonic dimers into pairs of fermionic atoms. We quantify deviations of atom-atom pair correlations from Wick's factorization scheme, and show that atom-molecule and molecule-molecule correlations grow with time, in clear departures from pairing mean-field theories. As a first-principles approach, the method provides benchmarking of approximate approaches and can be used to validate dynamical probes for characterizing strongly correlated phases of fermionic systems.
We formulate a general multi-mode Gaussian operator basis for fermions, to enable a positive phase-space representation of correlated Fermi states. The Gaussian basis extends existing bosonic phase-space methods to Fermi systems and thus allows first-principles dynamical or equilibrium calculations in quantum many-body Fermi systems. We prove the completeness of the basis and derive differential forms for products with one- and two-body operators. Because the basis satisfies fermionic superselection rules, the resulting phase space involves only c-numbers, without requiring anticommuting Grassmann variables. Furthermore, because of the overcompleteness of the basis, the phase-space distribution can always be chosen positive. This has important consequences for the sign problem in fermion physics. PACS numbers: 03.65.Ca, 42.50.Lc, 03.75.Ss, 71.10.Fd 1.
We give an overview of the Gaussian phase-space methods for fermions, focusing on the example of the Hubbard model, which describes ultracold fermions in a lattice. The phase-space methods do not suffer the sort of sign problem encountered by other quantum Monte Carlo methods, and provide an efficient means of calculating exact zero-temperature Hubbard results. However at finite temperature, there are known issues to do with the boundedness of the underlying distribution. We show how ‘stochastic gauges’ can be used to control the distribution tails and thereby improve the accuracy of the simulations at low temperatures.