We review phase-space simulation techniques for fermions, showing how a Gaussian operator basis leads to exact calculations of the evolution of a many-body quantum system in both real and imaginary ...
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 introduce a phase-space representation for qubits and spin models. The technique uses an SU(n) coherent-state basis and can equally be used for either static or dynamical simulations. We review previously known definitions and operator identities, and show how these can be used to define an off-diagonal, positive phase-space representation analogous to the positive-P function. As an illustration of the phase-space method, we use the example of the Ising model, which has exact solutions for the finite-temperature canonical ensemble in two dimensions. We show how a canonical ensemble for an Ising model of arbitrary structure can be efficiently simulated using SU(2) or atomic coherent states. The technique utilizes a transformation from a canonical (imaginary-time) weighted simulation to an equivalent unweighted real-time simulation. The results are compared to the exactly soluble two-dimensional case. We note that Ising models in one, two, or three dimensions are potentially achievable experimentally as a lattice gas of ultracold atoms in optical lattices. The technique is not restricted to canonical ensembles or to Ising-like couplings. It is also able to be used for real-time evolution and for systems whose time evolution follows a master equation describing decoherence and coupling to external reservoirs. The case of SU(n) phase space is used to describe n-level systems. In general, the requirement that time evolution be stochastic corresponds to a restriction to Hamiltonians and master equations that are quadratic in the group generators or generalized spin operators.
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.