Abstract The Green’s function has been an indispensable tool to study many-body systems that remain one of the biggest challenges in modern quantum physics for decades. The complicated calculation of the Green’s function impedes the research of many-body systems. The appearance of the noisy intermediate-scale quantum devices and quantum–classical hybrid algorithm inspire a new method to calculate the Green’s function. Here we design a programmable quantum circuit for photons with utilizing the polarization and the path degrees of freedom to construct a highly-precise variational quantum state of a photon, and first report the experimental realization for calculating the Green’s function of the two-site fermionic Hubbard model, a prototypical model for strongly-correlated materials, in photonic systems. We run the variational quantum eigensolver to obtain the ground state and excited states of the model, and then evaluate the transition amplitudes among the eigenstates. The experimental results present the spectral function of the Green’s function, which agrees well with exact results. Our demonstration provides the further possibility of the photonic system in quantum simulation and applications in solving complicated problems in many-body systems, biological science and so on.
Abstract The Green’s function has been an indispensable tool to study many-body systems that remain one of the biggest challenges in modern quantum physics for decades. The complicated calculation of the Green’s function impedes the research of many-body systems. The appearance of the noisy intermediate-scale quantum devices and quantum–classical hybrid algorithm inspire a new method to calculate the Green’s function. Here we design a programmable quantum circuit for photons with utilizing the polarization and the path degrees of freedom to construct a highly-precise variational quantum state of a photon, and first report the experimental realization for calculating the Green’s function of the two-site fermionic Hubbard model, a prototypical model for strongly-correlated materials, in photonic systems. We run the variational quantum eigensolver to obtain the ground state and excited states of the model, and then evaluate the transition amplitudes among the eigenstates. The experimental results present the spectral function of the Green’s function, which agrees well with exact results. Our demonstration provides the further possibility of the photonic system in quantum simulation and applications in solving complicated problems in many-body systems, biological science and so on.
Computation of the Green's function is crucial to study the properties of quantum many-body systems such as strongly correlated systems. Although the high-precision calculation of the Green's function is a notoriously challenging task on classical computers, the development of quantum computers may enable us to compute the Green's function with high accuracy even for classically-intractable large-scale systems. Here, we propose an efficient method to compute the real-time Green's function based on the local variational quantum compilation (LVQC) algorithm, which simulates the time evolution of a large-scale quantum system using a low-depth quantum circuit constructed through optimization on a smaller-size subsystem. Our method requires shallow quantum circuits to calculate the Green's function and can be utilized on both near-term noisy intermediate-scale and long-term fault-tolerant quantum computers depending on the computational resources we have. We perform a numerical simulation of the Green's function for the one- and two-dimensional Fermi-Hubbard model up to $4\ifmmode\times\else\texttimes\fi{}4$ sites lattice (32 qubits) and demonstrate the validity of our protocol compared to a standard method based on the Trotter decomposition. We finally present a detailed estimation of the gate count for the large-scale Fermi-Hubbard model, which also illustrates the advantage of our method over the Trotter decomposition.
An end-to-end strategy for hybrid quantum-classical computations of Green's functions in many-body systems is presented and applied to the pairing model. The scheme makes explicit use of the spectral representation of the Green's function, which entails the calculation of the N-body ground state as well as eigenstates and associated energies of the (N±1)-body neighbors. While the former is accessed via variational techniques, the latter are constructed by means of the quantum subspace expansion method. Different ansatzes for the ground-state wave function, originating from either classical or quantum approaches, are tested and compared to exact calculations. The resulting one-body Green's functions prove to be accurate approximations of the exact one for a large range of parameters, including across the normal-to-superfluid transition. As a byproduct, this approach yields a good description of odd systems provided that the starting even system is well reproduced by the variational ansatz.
Approximating the ground states of strongly interacting electron systems in quantum chemistry and condensed matter physics is expected to be one of the earliest applications of quantum computers. In this paper, we prepare highly accurate ground states of the Fermi-Hubbard model for small grids up to 6 sites (12 qubits) by using an interpretable, adaptive variational quantum eigensolver(VQE) called ADAPT-VQE. In contrast with non-adaptive VQE, this algorithm builds a system-specific ansatz by adding an optimal gate built from one-body or two-body fermionic operators at each step. We show this adaptive method outperforms the non-adaptive counterpart in terms of fewer variational parameters, short gate depth, and scaling with the system size. The fidelity and energy of the prepared state appear to improve asymptotically with ansatz depth. We also demonstrate the application of adaptive variational methods by preparing excited states and Green functions using a proposed ADAPT-SSVQE algorithm. Lower depth, asymptotic convergence, noise tolerance of a variational approach, and a highly controllable, system-specific ansatz make the adaptive variational methods particularly well-suited for NISQ devices.
In this paper, we present a quantum computational method to calculate the many-body Green's function matrix in a spin orbital basis. We apply our approach to finite-sized fermionic Hubbard models and related impurity models within Dynamical Mean Field Theory, and demonstrate the calculation of Green's functions on Quantinuum's H1-1 trapped-ion quantum computer. Our approach involves a cumulant expansion of the Lanczos method, using Hamiltonian moments as measurable expectation values. This bypasses the need for a large overhead in the number of measurements due to repeated applications of the variational quantum eigensolver (VQE), and instead measures the expectation value of the moments with one set of measurement circuits. From the measured moments, the tridiagonalised Hamiltonian matrix can be computed, which in turn yields the Green's function via continued fractions. While we use a variational algorithm to prepare the ground state in this work, we note that the modularity of our implementation allows for other (non-variational) approaches to be used for the ground state.
Quantum computers have the potential to deliver speed-ups for solving certain important problems that are intractable for classical counterparts, making them a promising avenue for advancing modern computation. However, many quantum algorithms require deep quantum circuits, which are challenging to implement on current noisy devices. To address this limitation, variational quantum algorithms (VQAs) have been actively developed, enabling practical quantum computing in the noisy intermediate-scale quantum (NISQ) era. Among them, the variational quantum eigensolver (VQE) stands out as a leading approach for solving problems in quantum chemistry, many-body physics, and even integer factorization. The VQE algorithm can be implemented on various quantum hardware platforms, including photonic systems, quantum dots, trapped ions, neutral atoms, and superconducting circuits. In particular, photonic platforms offer several advantages: they operate at room temperature, exhibit low decoherence, and support multiple degrees of freedom, making them suitable for scalable, high-dimensional quantum computation. Here we present methodologies for realizing VQE on photonic systems, highlighting their potential for practical quantum computing. We first provide a theoretical overview of the VQE framework, focusing on the procedure for variationally estimating ground state energies. We then explore how photonic systems can implement these processes, showing that a wide variety of problems can be addressed using either multiple qubit states or a single qudit state.
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.