Dynamical mean field theory algorithm and experiment on quantum computers
Ivan Rungger, Nathan Fitzpatrick, Hongxiang Chen, C. Huerta Alderete, Harriet Apel, Alexander Cowtan, Patterson, A., David Muñoz Ramo, Yixuan Zhu, Nhung H. Nguyen, Edward Grant, Stéphane Chrétien, Leonard Wossnig, Norbert M. Linke, Duncan, R.
The developments of quantum computing algorithms and experiments for atomic scale simulations have largely focused on quantum chemistry for molecules, while their application in condensed matter systems is scarcely explored. Here we present a quantum algorithm to perform dynamical mean field theory (DMFT) calculations for condensed matter systems on currently available quantum computers, and demonstrate it on two quantum hardware platforms. DMFT is required to properly describe the large class of materials with strongly correlated electrons. The computationally challenging part arises from solving the effective problem of an interacting impurity coupled to a bath, which scales exponentially with system size on conventional computers. An exponential speedup is expected on quantum computers, but the algorithms proposed so far are based on real time evolution of the wavefunction, which requires high-depth circuits and hence very low noise levels in the quantum hardware. Here we propose an alternative approach, which uses the variational quantum eigensolver (VQE) method for ground and excited states to obtain the needed quantities as part of an exact diagonalization impurity solver. We present the algorithm for a two site DMFT system, which we benchmark using simulations on conventional computers as well as experiments on superconducting and trapped ion qubits, demonstrating that this method is suitable for running DMFT calculations on currently available quantum hardware.
The developments of quantum computing algorithms and experiments for atomic scale simulations have largely focused on quantum chemistry for molecules, while their application in condensed matter systems is scarcely explored. Here we present a quantum algorithm to perform dynamical mean field theory (DMFT) calculations for condensed matter systems on currently available quantum computers, and demonstrate it on two quantum hardware platforms. DMFT is required to properly describe the large class of materials with strongly correlated electrons. The computationally challenging part arises from solving the effective problem of an interacting impurity coupled to a bath, which scales exponentially with system size on conventional computers. An exponential speedup is expected on quantum computers, but the algorithms proposed so far are based on real time evolution of the wavefunction, which requires high-depth circuits and hence very low noise levels in the quantum hardware. Here we propose an alternative approach, which uses the variational quantum eigensolver (VQE) method for ground and excited states to obtain the needed quantities as part of an exact diagonalization impurity solver. We present the algorithm for a two site DMFT system, which we benchmark using simulations on conventional computers as well as experiments on superconducting and trapped ion qubits, demonstrating that this method is suitable for running DMFT calculations on currently available quantum hardware.
The developments of quantum computing algorithms and experiments for atomic scale simulations have largely focused on quantum chemistry for molecules, while their application in condensed matter systems is scarcely explored. Here we present a quantum algorithm to perform dynamical mean field theory (DMFT) calculations for condensed matter systems on currently available quantum computers, and demonstrate it on two quantum hardware platforms. DMFT is required to properly describe the large class of materials with strongly correlated electrons. The computation-ally challenging part arises from solving the effective problem of an interacting impurity coupled to a bath, which scales exponentially with system size on conventional computers. An exponential speedup is expected on quantum computers, but the algorithms proposed so far are based on real time evolution of the wavefunction, which requires high-depth circuits and hence very low noise levels in the quantum hardware. Here we propose an alternative approach, which uses the variational quantum eigensolver (VQE) method for ground and excited states to obtain the needed quantities as part of an exact diagonalization impurity solver. We present the algorithm for a two site DMFT system, which we benchmark using simulations on conventional computers as well as experiments on superconducting and trapped ion qubits, demonstrating that this method is suitable for running DMFT calculations on currently available quantum hardware.
Quantum computers (QC) could harbor the potential to significantly advance materials simulations, particularly at the atomistic scale involving strongly correlated fermionic systems, where an accurate description of quantum many-body effects scales unfavorably with size. While a full-scale treatment of condensed matter systems with currently available noisy quantum computers remains elusive, quantum embedding schemes like dynamical mean-field theory (DMFT) allow the mapping of an effective, reduced subspace Hamiltonian to available devices to improve the accuracy of ab initio calculations such as density functional theory (DFT). Here, we report on the development of a hybrid quantum-classical DFT + DMFT simulation framework which relies on a quantum impurity solver based on the Lehmann representation of the impurity Green’s function. Hardware experiments with up to 14 qubits on the IBM Quantum system are conducted, using advanced error mitigation methods and a novel calibration scheme for an improved zero-noise extrapolation to effectively reduce adverse effects from inherent noise on current quantum devices. We showcase the utility of our quantum DFT + DMFT workflow by assessing the correlation effects on the electronic structure of a real material, Ca 2 CuO 2 Cl 2 , which is mapped to an effective single-band Hubbard Hamiltonian and the subsequently derived Anderson impurity model solved with up to 6 bath sites on available quantum hardware. Further, we carefully benchmark our quantum results with respect to exact reference solutions and experimental spectroscopy measurements. While challenges remain to scale our approach to larger, multi-orbital and multi-site systems with more bath sites, the present work marks an important milestone towards achieving utility-scale quantum computation in materials simulation.
The accurate theoretical description of materials with strongly correlated electrons is a formidable challenge in condensed matter physics and computational chemistry. Dynamical Mean Field Theory (DMFT) is a successful approach that predicts behaviors of such systems by incorporating some of the correlated behavior using an impurity model, but it is limited by the need to calculate the impurity Green’s function. This work proposes a framework for DMFT calculations on quantum computers, focusing on near-term applications. It leverages the structure of the impurity problem, combining a low-rank Gaussian subspace representation of the ground state and a compressed, short-depth quantum circuit that joins state preparation with time evolution to compute Green’s functions. We demonstrate the convergence of the DMFT algorithm using the Gaussian subspace in a noise-free setting, and show the hardware viability of circuit compression by extracting the impurity Green’s function on IBM quantum processors for a single impurity coupled to three bath orbitals (8 qubits, 1 ancilla). We discuss potential paths toward realizing this quantum computing use case in materials science.
Quantitative descriptions of strongly correlated materials pose a considerable challenge in condensed matter physics and chemistry. A promising approach to address this problem is quantum embedding methods. In particular, the dynamical mean-field theory (DMFT) maps the original system to an effective quantum impurity model comprising correlated orbitals embedded in an electron bath. The biggest bottleneck in DMFT calculations is numerically solving the quantum impurity model, i.e., computing the Green's function. Past studies have proposed theoretical methods to compute the Green's function of a quantum impurity model in polynomial time using a quantum computer. So far, however, efficient methods for computing the imaginary-time Green's functions have not been established despite the advantages of the imaginary-time formulation. We propose a quantum-classical hybrid algorithm for computing imaginary-time Green's functions on quantum devices with limited hardware resources by applying the variational quantum simulation. Using a quantum circuit simulator, we verified this algorithm by computing Green's functions for a dimer model as well as a four-site impurity model obtained by DMFT calculations of the single-band Hubbard model, although our method can be applied to general imaginary-time correlation functions.
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.
We propose an algorithm to obtain Green's functions as a continued fraction on quantum computers, which is based on the construction of the Krylov basis using variational quantum algorithms, and included in a Lanczos iterative scheme. This allows the integration of quantum algorithms with first principles material science simulations, as we demonstrate within the dynamical mean-field theory (DMFT) framework. DMFT enables quantitative predictions for strongly correlated materials, and relies on the calculation of Green's functions. On conventional computers the exponential growth of the Hilbert space with the number of orbitals limits DMFT to small systems. Quantum computers open new avenues and can lead to a significant speedup in the computation of expectation values required to obtain the Green's function. We apply our Krylov variational quantum algorithm combined with DMFT to the charge transfer insulator La2CuO4 using a quantum computing emulator, and show that with 8 qubits it predicts the correct insulating material properties for the paramagnetic phase. We therefore expect that the method is ideally suited to perform simulations for real materials on near term quantum hardware.