We propose a non-linear, hybrid quantum-classical scheme for simulating non-equilibrium dynamics of strongly correlated fermions described by the Hubbard model in a Bethe lattice in the thermodynamic limit. Our scheme implements non-equilibrium dynamical mean field theory (DMFT) and uses a digital quantum simulator to solve a quantum impurity problem whose parameters are iterated to self-consistency via a classically computed feedback loop where quantum gate errors can be partly accounted for. We analyse the performance of the scheme in an example case.
We propose a non-linear, hybrid quantum-classical scheme for simulating non-equilibrium dynamics of strongly correlated fermions described by the Hubbard model in a Bethe lattice in the thermodynamic limit. Our scheme implements non-equilibrium dynamical mean field theory (DMFT) and uses a digital quantum simulator to solve a quantum impurity problem whose parameters are iterated to self-consistency via a classically computed feedback loop where quantum gate errors can be partly accounted for. We analyse the performance of the scheme in an example case.
This work is a theoretical study of fermionic models. We focus on problems where highly controllable quantum simulators of these models have an important role, and we utilise both the analogue and the digital paradigm of quantum simulation. In the part on analogue quantum simulation, we focus on the proposed 'spin-asymmetric' Josephson effect where Cooper-paired spins display frequency synchronized Josephson oscillations with spin-dependent amplitudes. We consider different scenarios where the phenomenon could manifest in ultracold atomic Fermi gases. We study a Fermi gas Josephson junction in the recently realized Josephson plasma oscillation regime with an additional spin-dependent potential and show that the asymmetry in the resulting spin-dependent plasma oscillation amplitudes is on the order of a couple of per cent. We also demonstrate numerically that spin-asymmetric Josephson-like currents occur in a one-dimensional spin-dependent optical superlattice, with amplitude asymmetries up to 39%. Finally, we show that at zero temperature the tunable critical current in ferromagnetic Josephson junctions can be explained by the spin-asymmetric Josephson effect. In the part where digital quantum simulation is used, we propose a hybrid quantum-classical approach to studying strongly correlated fermion models. In this approach, a digital quantum simulator works in conjunction with a classical feedback loop to solve the infinite-dimensional Hubbard model directly in the thermodynamic limit. The scheme implements the well-established dynamical mean-field theory (DMFT) method, such that the digital quantum simulator solves the classically hard DMFT impurity problem and self-consistency is taken care of in a classical computer. We first present a few-qubit proof-of-principle setup for equilibrium systems that implements the simplified 'two-site' DMFT. This few-qubit setup is used for a qualitative description of the Mott transition in the half-filled infinite-dimensional Hubbard model. We then describe a scalable setup for simulating non-equilibrium many-body quantum dynamics by proposing the implementation of the non-equilibrium extension of DMFT with the hybrid device.
We study a proof-of-principle example of the recently proposed hybrid quantum-classical simulation of strongly correlated fermion models in the thermodynamic limit. In a ‘two-site’ dynamical mean-field theory (DMFT) approach we reduce the Hubbard model to an effective impurity model subject to self-consistency conditions. The resulting minimal two-site representation of the non-linear hybrid setup involves four qubits implementing the impurity problem, plus an ancilla qubit on which all measurements are performed. We outline a possible implementation with superconducting circuits feasible with near-future technology.
This work is motivated by the fact that the investigation of non-equilibrium phenomena in strongly correlated electron systems has developed into one of the most active and exciting branches of condensed matter physics as it provides rich new insights that could not be obtained from the study of equilibrium situations. However, a theoretical description of those phenomena is missing. Therefore, in this thesis, we develop a numerical method that can be used to study two minimal models -- the Hubbard model and the Anderson impurity model with general parameter range and time dependence. We begin by introducing the theoretical framework and the general features of the Hubbard model. We then describe the dynamical mean field theory (DMFT), which was first invented by Georges in 1992. It provides a feasible way to approach strongly correlated electron systems and reduces the complexity of the calculations via a mapping of lattice models onto quantum impurity models subject to a self-consistency condition. We employ the non-equilibrium extension of DMFT and map the Hubbard model to the single impurity Anderson model (SIAM). Since the fundamental component of the DMFT method is a solver of the single impurity Anderson model, we continue with a description of the formalism to study the real-time dynamics of the impurity model staring at its thermal equilibrium state. We utilize the non-equilibrium strong-coupling perturbation theory and derive semi-analytical approximation methods such as the non-crossing approximation (NCA) and the one-crossing approximation (OCA). We then use the Quantum Monte-Carlo method (QMC) as a numerically exact method and present proper measurements of local observables, current and Green's functions. We perform simulations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibrium times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in the steady state. However, this bare QMC solver suffers from a dynamical sign problem for long time propagations. To overcome the limitations of this bare treatment, we introduce the ``Inchworm algorithm'', based on iteratively reusing the information obtained in previous steps to extend the propagation to longer times and stabilize the calculations. We show that this algorithm greatly reduces the required order for each simulation and re-scales the exponential challenge to quadratic in time. We introduce a method to compute Green's functions, spectral functions, and currents for inchworm Monte Carlo and show how systematic error assessments in real time can be obtained. We illustrate the capabilities of the algorithm with a study of the behavior of quantum impurities after an instantaneous voltage quench from a thermal equilibrium state. We conclude with the applications of the unbiased inchworm impurity solver to DMFT calculations. We employ the methods for a study of the one-band paramagnetic Hubbard model on the Bethe lattice in equilibrium, where the DMFT approximation becomes exact. We begin with a brief introduction of the Mott metal insulator phase diagram. We present the results of both real time Green's functions and spectral functions from our nonequilibrium calculations. We observe the metal-insulator crossover as the on-site interaction is increased and the formation of a quasi-particle peak as the temperature is lowered. We also illustrate the convergence of our algorithms in different aspects.
The Dynamical Mean Field Theory (DMFT) is a powerful tool for calculating highly correlated systems (both bosonic and fermionic) in a state of thermodynamic equilibrium. However, in the case of non-equilibrium states, the method has significant limitations that do not allow obtaining correct results. The stumbling block here is the impurity solver: a method for calculating the dynamics of an open system. In this work we present the prototype of a universal impurity solver for the DMFT method which can solve both equilibrium and non-equilibrium problems. We analyse spectral functions of Bose-Hubbard model on the Bethe lattice in different regimes and show that our solver gives correct and physical results.
We present a time-domain iteration scheme for solving the Dynamical Mean-Field Theory (DMFT) self-consistent equations using retarded Green's functions in real time. Unlike conventional DMFT approaches that operate in imaginary time or frequency space, our scheme operates directly with real-time quantities. This makes it particularly suitable for near-term quantum computing hardware with limited Hilbert spaces, where real-time propagation can be efficiently implemented via Trotterization or variational quantum algorithms. We map the effective impurity problem to a finite one-dimensional chain with a small number of bath sites, solved via exact diagonalization as a proof-of-concept. The hybridization function is iteratively updated through time-domain fitting until self-consistency. We demonstrate stable convergence across a wide range of interaction strengths for the half-filled Hubbard model on a Bethe lattice, successfully capturing the metal-to-insulator transition. Despite using limited time resolution and a minimal bath discretization, the spectral functions clearly exhibit the emergence of Hubbard bands and the suppression of spectral weight at the Fermi level as interaction strength increases. This overcomes major limitations of two-site DMFT approximations by delivering detailed spectral features while preserving efficiency and compatibility with quantum computing platforms through real-time dynamics.
Dynamical mean-field theory (DMFT) is a widely used, nonperturbative computational framework for strongly correlated electron systems. This thesis extends DMFT in two directions: treating disorder and higher-order correlators away from equilibrium, and using near-term quantum computers as impurity solvers. In the first part, we develop two methodological advances for the nonequilibrium DMFT treatment of disordered, interacting systems. The first overcomes the central computational bottleneck of the recently developed nonequilibrium DMFT plus coherent potential approximation (DMFT+CPA) by replacing its explicit numerical disorder averaging with an exact analytical integration, exploiting a normalization property of the Schwinger-Keldysh contour; applied to a disorder-and-interaction quench of the Anderson-Hubbard model, this reveals that weak random disorder promotes thermalization. The second demonstrates that the out-of-time-order correlator (OTOC) can be successfully computed within the existing DMFT+CPA framework, by generalizing the Schwinger-Keldysh contour to a double-folded structure; as a proof-of-principle demonstration, weak disorder is found to accelerate the OTOC's decay in the half-filled Anderson-Hubbard model in the intermediate-coupling regime. In the second part, we explore the use of near-term quantum computers as DMFT impurity solvers. A stable real-time DMFT iteration scheme, designed for compatibility with quantum hardware, is first developed and validated classically against the metal-to-insulator transition. A symmetry-adapted variational quantum eigensolver (VQE) is then benchmarked as the Anderson impurity solver for the converged DMFT Hamiltonian, and is found to reproduce ground-state energies accurately and to extract the single-particle Green's function reliably at intermediate to strong coupling, while accurate ground-state energies are shown not to guarantee accurate dynamics at weak coupling, a key challenge for extending such approaches to correlated metals. Together, we demonstrate that DMFT can be methodologically extended to study disordered and interacting systems away from equilibrium, and to incorporate near-term quantum hardware as a genuinely new computational resource to enhance its capabilities.