Dynamical quantum-cluster approaches, such as different cluster extensions of the dynamical mean-field theory (cluster DMFT) or the variational cluster approximation (VCA), combined with efficient cluster solvers, such as the quantum Monte Carlo (QMC) method, provide controlled approximations of the single-particle Green's function for lattice models of strongly correlated electrons. To access the thermodynamics, however, a thermodynamical potential is needed. We present an efficient numerical algorithm to compute the grand potential within cluster-embedding approaches that are based on novel continuous-time QMC schemes. It is shown that the numerically exact cluster grand potential can be obtained from a quantum Wang-Landau technique to reweight the coefficients in the expansion of the partition function. The lattice contributions to the grand potential are computed by a proper infinite summation over Matsubara frequencies. A proof of principle is given by applying the VCA to antiferromagnetic (short-range) order in the two-dimensional Hubbard model at finite temperatures.
In the framework of the hybrid quantum-classical variational cluster approach to strongly correlated electron systems one of the goals of a quantum subroutine is to find single-particle correlation functions of lattice fermions in polynomial time. Previous works suggested to use variants of the Hadamard test for this purpose, which requires an implementation of controlled single-particle fermionic operators. However, for a number of locality-preserving mappings to encode fermions into qubits, a direct construction of such operators is not possible. In this work, we propose a quantum algorithm that uses an analog of the Kubo formula adapted to a quantum circuit simulating the Hubbard model. It allows to access the Green's function of a cluster directly using only bilinears of fermionic operators and circumvents the usage of the Hadamard test. We test our algorithm in practice by using open-access simulators of noisy IBM superconducting chips.
Quantum computers are the ideal platform for quantum simulations. Given enough coherent operations and qubits, such machines can be leveraged to simulate strongly correlated materials, where intricate quantum effects give rise to counterintuitive macroscopic phenomena such as high-temperature superconductivity. In this paper, we provide a gate decomposition and an architecture for a quantum simulator used to simulate the Fermi-Hubbard model in a hybrid variational quantum-classical algorithm. We propose a simple planar implementation-independent layout of qubits that can also be used to simulate more general fermionic systems. By working through a concrete application, we show the gate decomposition used to simulate the Hamiltonian of a cluster of the Fermi-Hubbard model. We briefly analyze the Trotter-Suzuki errors and estimate the scaling properties of the algorithm for more complex applications.
Preparing Gibbs states on a quantum computer offers a way to computationally probe how microscopic interactions give rise to macroscopic phenomena at finite temperature. However, an uncomfortable tension has set in: our strongest general guarantees for Gibbs state preparation hold in regimes where the resulting states can also be simulated classically. We propose investigating physical models with competing energy scales. The Fermi--Hubbard model is a natural target: the competition between electron hopping, on-site repulsion and chemical potential encapsulates many mysteries of condensed matter physics, including Mott insulators, striping and superconductivity.
Our main result identifies a simple guiding principle: the temperature scale for efficient preparation is set by the couplings between sites, while the on-site energy scale is unrestricted. For the Fermi--Hubbard model on an
n-site nearest-neighbor lattice of fixed dimension, we construct a quantum Gibbs sampler that mixes in
O(log(n/ε)) time whenever the temperature is at least a constant factor larger than the hopping strength, but is allowed to be independent of the repulsion and chemical potential. We obtain analogous guarantees for local qubit systems with bounded intersite interactions and arbitrary external fields. As a consequence, our algorithms can efficiently prepare Fermi--Hubbard Gibbs states that are provably outside the convex hull of fermionic Gaussian states and qubit Gibbs states for locally interacting Hamiltonians outside the set of product states.
We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.
We formulate a finite-temperature scheme of the variational cluster approximation (VCA) particularly suitable for an exact-diagonalization cluster solver. Based on the analytical properties of the single-particle Green's function matrices, we explicitly show the branch-cut structure of logarithm of the complex determinant functions appearing in the self-energy-functional theory (SFT) and whereby construct an efficient scheme for the finite-temperature VCA. We also derive the explicit formulas for entropy and specific heat within the framework of the SFT. We first apply the method to explore the antiferromagnetic order in a half-filled Hubbard model by calculating the entropy, specific heat, and single-particle excitation spectrum for different values of on-site Coulomb repulsion
U and temperature
T. We also calculate the
T dependence of the single-particle excitation spectrum in the strong coupling region, and discuss the overall similarities to and the fine differences from the spectrum obtained by the spin-density-wave mean-field theory at low temperatures and the Hubbard-I approximation at high temperatures. Moreover, we show a necessary and sufficient condition for the third law of thermodynamics in the SFT. On the basis of the thermodynamic properties, such as the entropy and the double occupancy, calculated via the
T and/or
U derivative of the grand potential, we obtain a crossover diagram in the
(U,T) plane, which separates a Slater-type insulator and a Mott-type insulator. Next, we demonstrate the finite-temperature scheme in the cluster-dynamical-impurity approximation (CDIA), i.e., the VCA with noninteracting bath orbitals attached to each cluster, and study the paramagnetic Mott metal-insulator transition in the half-filled Hubbard model. Formulating the finite-temperature CDIA, we first address a subtle issue regarding the treatment of the artificially introduced bath degrees of freedom which are absent in the originally considered Hubbard model. We then apply the finite-temperature CDIA to calculate the finite-temperature phase diagram in the
(U,T) plane. Metallic, insulating, coexistence, and crossover regions are distinguished from the bath-cluster hybridization-variational-parameter dependence of the grand-potential functional. We find that the Mott transition at low temperatures is discontinuous, and the coexistence region of the metallic and insulating states persists down to zero temperature. The result obtained here by the finite-temperature CDIA is complementary to the previously reported zero-temperature CDIA phase diagram.