One of the most promising applications for near term quantum computers is the simulation of physical quantum systems, particularly many-electron systems in chemistry and condensed matter physics. In solid state physics, finding the correct symmetry broken ground state of an interacting electron system is one of the central challenges. The Variational Hamiltonian Ansatz (VHA), a variational hybrid quantum-classical algorithm especially suited for finding the ground state of a solid state system, will in general not prepare a broken symmetry state unless the initial state is chosen to exhibit the correct symmetry. In this work, we discuss three variations of the VHA designed to find the correct broken symmetry states close to a transition point between different orders. As a test case we use the two-dimensional Hubbard model where we break the symmetry explicitly by means of external fields coupling to the Hamiltonian and calculate the response to these fields. For the calculation we simulate a gate-based quantum computer and also consider the effects of dephasing noise on the algorithms. We find that two of the three algorithms are in good agreement with the exact solution for the considered parameter range. The third algorithm agrees with the exact solution only for a part of the parameter regime, but is more robust with respect to dephasing compared to the other two algorithms.
One of the most promising applications for near term quantum computers is the simulation of physical quantum systems, particularly many-electron systems in chemistry and condensed matter physics. In solid state physics, finding the correct symmetry broken ground state of an interacting electron system is one of the central challenges. The Variational Hamiltonian Ansatz (VHA), a variational hybrid quantum-classical algorithm especially suited for finding the ground state of a solid state system, will in general not prepare a broken symmetry state unless the initial state is chosen to exhibit the correct symmetry. In this work, we discuss three variations of the VHA designed to find the correct broken symmetry states close to a transition point between different orders. As a test case we use the two-dimensional Hubbard model where we break the symmetry explicitly by means of external fields coupling to the Hamiltonian and calculate the response to these fields. For the calculation we simulate a gate-based quantum computer and also consider the effects of dephasing noise on the algorithms. We find that two of the three algorithms are in good agreement with the exact solution for the considered parameter range. The third algorithm agrees with the exact solution only for a part of the parameter regime, but is more robust with respect to dephasing compared to the other two algorithms.
Abstract One of the most promising applications for near term quantum computers is the simulation of physical quantum systems, particularly many-electron systems in chemistry and condensed matter physics. In solid state physics, finding the correct symmetry broken ground state of an interacting electron system is one of the central challenges. To help finding the correct broken symmetries in the thermodynamic limit methods that allow to determine the groundstate of large but finite interacting electron systems are very useful. The variational Hamiltonian ansatz (VHA), a variational hybrid quantum-classical algorithm especially suited for finding the ground state of a solid state system, will in general not prepare a broken symmetry state unless the initial state is chosen to exhibit the correct symmetry. In this work, we discuss three variations of the VHA designed to find the symmetry-breaking groundstate of a finite system close to a transition point between different orders. As a test case we use the two-dimensional Hubbard model where we break the symmetry explicitly by means of external fields coupling to the Hamiltonian and calculate the response to these fields. For the calculation we simulate a gate-based quantum computer and also consider the effects of dephasing noise on the algorithms. We find that two of the three algorithms are in good agreement with the exact solution for the considered parameter range. The third algorithm agrees with the exact solution only for a part of the parameter regime, but is more robust with respect to dephasing compared to the other two algorithms.
Abstract A key goal of digital quantum computing is the simulation of fermionic systems such as molecules or the Hubbard model. Unfortunately, for present and near-future quantum computers the use of quantum error correction schemes is still out of reach. Hence, the finite error rate limits the use of quantum computers to algorithms with a low number of gates. The variational Hamiltonian ansatz (VHA) has been shown to produce the ground state in good approximation in a manageable number of steps. Here we study explicitly the effect of gate errors on its performance. The VHA is inspired by the adiabatic quantum evolution under the influence of a time-dependent Hamiltonian, where the—ideally short—fixed Trotter time steps are replaced by variational parameters. The method profits substantially from quantum variational error suppression, e.g. unitary quasi-static errors are mitigated within the algorithm. We test the performance of the VHA when applied to the Hubbard model in the presence of unitary control errors on quantum computers with realistic gate fidelities.
Variational algorithms are promising candidates to be implemented on near-term quantum computers. In our work, we investigate the variational Hamiltonian ansatz (VHA), where a parametrized trial state of the quantum mechanical wave function is optimized to obtain the ground-state energy. In the VHA, the trial state is given by a noninteracting reference state modified by unitary rotations using generators that are part of the Hamiltonian describing the system. The lowest energy is obtained by optimizing the angles of those unitary rotations. A standard procedure to optimize the variational parameters is to use gradient-based algorithms. However, shot noise and the intrinsic noise of the quantum device affect the evaluation of the required gradients. We study how different methods for obtaining the gradient, specifically the finite-difference and the parameter-shift rule, are affected by shot noise and the noise of the quantum computer. To this end, we simulate a simple quantum circuit, as well as the two-site and six-site Hubbard models.
Abstract We present a detailed analysis and optimization of the variational quantum algorithms required to find the ground state of a correlated electron model, using several types of variational ansatz. Specifically, we apply our approach to the atomic limit of the Anderson model, which is widely studied in condensed matter physics since it can simulate fundamental physical phenomena, ranging from magnetism to superconductivity. The method is developed by presenting efficient state preparation circuits that exhibit total spin, spin projection, particle number and time-reversal symmetries. These states contain the minimal number of variational parameters needed to fully span the appropriate symmetry subspace allowing to avoid irrelevant sectors of Hilbert space. Then, we show how to construct quantum circuits, providing explicit decomposition and gate count in terms of standard gate sets. We test these quantum algorithms looking at ideal quantum computer simulations as well as implementing quantum noisy simulations. We finally perform an accurate comparative analysis among the approaches implemented, highlighting their merits and shortcomings.
The Variational Quantum Eigensolver (VQE), as a hybrid quantum-classical algorithm, is an important tool for effective quantum computing in the current noisy intermediate-scale quantum (NISQ) era. However, the traditional hardware-efficient ansatz without taking into account symmetries requires more computational resources to explore the unnecessary regions in the Hilbert space. The conventional Subspace-Search VQE (SSVQE) algorithm, which can calculate excited states, is also unable to effectively handle degenerate states since the loss function only contains the expectation value of the Hamiltonian. In this study, the energy eigenstates of the one-dimensional Fermi-Hubbard model with two lattice sites and the two-dimensional Hubbard model with four lattice sites are calculated. By incorporating symmetries into the quantum circuits and loss function, we find that both the ground state and excited state calculations are improved greatly compared to the case without symmetries. The enhancement in excited state calculations is particularly significant. This is because quantum circuits that conserve the particle number are used, and appropriate penalty terms are added to the loss function, enabling the optimization process to correctly identify degenerate states. The results are verified through repeated simulations.
Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.