-
Quantum algorithms for the variational optimization of correlated electronic states with stochastic reconfiguration and the linear method
Authors:
Mario Motta,
Kevin J. Sung,
James Shee
Abstract:
Solving the electronic Schrodinger equation for strongly correlated ground states is a long-standing challenge. We present quantum algorithms for the variational optimization of wavefunctions correlated by products of unitary operators, such as Local Unitary Cluster Jastrow (LUCJ) ansatzes, using stochastic reconfiguration and the linear method. While an implementation on classical computing hardw…
▽ More
Solving the electronic Schrodinger equation for strongly correlated ground states is a long-standing challenge. We present quantum algorithms for the variational optimization of wavefunctions correlated by products of unitary operators, such as Local Unitary Cluster Jastrow (LUCJ) ansatzes, using stochastic reconfiguration and the linear method. While an implementation on classical computing hardware would require exponentially growing compute cost, the cost (number of circuits and shots) of our quantum algorithms is polynomial in system size. We find that classical simulations of optimization with the linear method consistently find lower energy solutions than with the L-BFGS-B optimizer across the dissociation curves of the notoriously difficult N$_2$ and C$_2$ dimers; LUCJ predictions of the ground-state energies deviate from exact diagonalization by 1 kcal/mol or less at all points on the potential energy curve. While we do characterize the effect of shot noise on the LM optimization, these noiseless results highlight the critical but often overlooked role that optimization techniques must play in attacking the electronic structure problem (on both classical and quantum hardware), for which even mean-field optimization is formally NP hard. We also discuss the challenge of obtaining smooth curves in these strongly correlated regimes, and propose a number of quantum-friendly solutions ranging from symmetry-projected ansatz forms to a symmetry-constrained optimization algorithm.
△ Less
Submitted 3 August, 2024;
originally announced August 2024.
-
Chemistry Beyond Exact Solutions on a Quantum-Centric Supercomputer
Authors:
Javier Robledo-Moreno,
Mario Motta,
Holger Haas,
Ali Javadi-Abhari,
Petar Jurcevic,
William Kirby,
Simon Martiel,
Kunal Sharma,
Sandeep Sharma,
Tomonori Shirakawa,
Iskandar Sitdikov,
Rong-Yang Sun,
Kevin J. Sung,
Maika Takita,
Minh C. Tran,
Seiji Yunoki,
Antonio Mezzacapo
Abstract:
A universal quantum computer can be used as a simulator capable of predicting properties of diverse quantum systems. Electronic structure problems in chemistry offer practical use cases around the hundred-qubit mark. This appears promising since current quantum processors have reached these sizes. However, mapping these use cases onto quantum computers yields deep circuits, and for for pre-fault-t…
▽ More
A universal quantum computer can be used as a simulator capable of predicting properties of diverse quantum systems. Electronic structure problems in chemistry offer practical use cases around the hundred-qubit mark. This appears promising since current quantum processors have reached these sizes. However, mapping these use cases onto quantum computers yields deep circuits, and for for pre-fault-tolerant quantum processors, the large number of measurements to estimate molecular energies leads to prohibitive runtimes. As a result, realistic chemistry is out of reach of current quantum computers in isolation. A natural question is whether classical distributed computation can relieve quantum processors from parsing all but a core, intrinsically quantum component of a chemistry workflow. Here, we incorporate quantum computations of chemistry in a quantum-centric supercomputing architecture, using up to 6400 nodes of the supercomputer Fugaku to assist a Heron superconducting quantum processor. We simulate the N$_2$ triple bond breaking in a correlation-consistent cc-pVDZ basis set, and the active-space electronic structure of [2Fe-2S] and [4Fe-4S] clusters, using 58, 45 and 77 qubits respectively, with quantum circuits of up to 10570 (3590 2-qubit) quantum gates. We obtain our results using a class of quantum circuits that approximates molecular eigenstates, and a hybrid estimator. The estimator processes quantum samples, produces upper bounds to the ground-state energy and wavefunctions supported on a polynomial number of states. This guarantees an unconditional quality metric for quantum advantage, certifiable by classical computers at polynomial cost. For current error rates, our results show that classical distributed computing coupled to quantum processors can produce good approximate solutions for practical problems beyond sizes amenable to exact diagonalization.
△ Less
Submitted 8 May, 2024;
originally announced May 2024.
-
Subspace methods for electronic structure simulations on quantum computers
Authors:
Mario Motta,
William Kirby,
Ieva Liepuoniute,
Kevin J. Sung,
Jeffrey Cohn,
Antonio Mezzacapo,
Katherine Klymko,
Nam Nguyen,
Nobuyuki Yoshioka,
Julia E. Rice
Abstract:
Quantum subspace methods (QSMs) are a class of quantum computing algorithms where the time-independent Schrodinger equation for a quantum system is projected onto a subspace of the underlying Hilbert space. This projection transforms the Schrodinger equation into an eigenvalue problem determined by measurements carried out on a quantum device. The eigenvalue problem is then solved on a classical c…
▽ More
Quantum subspace methods (QSMs) are a class of quantum computing algorithms where the time-independent Schrodinger equation for a quantum system is projected onto a subspace of the underlying Hilbert space. This projection transforms the Schrodinger equation into an eigenvalue problem determined by measurements carried out on a quantum device. The eigenvalue problem is then solved on a classical computer, yielding approximations to ground- and excited-state energies and wavefunctions. QSMs are examples of hybrid quantum-classical methods, where a quantum device supported by classical computational resources is employed to tackle a problem. QSMs are rapidly gaining traction as a strategy to simulate electronic wavefunctions on quantum computers, and thus their design, development, and application is a key research field at the interface between quantum computation and electronic structure. In this review, we provide a self-contained introduction to QSMs, with emphasis on their application to the electronic structure of molecules. We present the theoretical foundations and applications of QSMs, and we discuss their implementation on quantum hardware, illustrating the impact of noise on their performance.
△ Less
Submitted 30 November, 2023;
originally announced December 2023.
-
Simulating Majorana zero modes on a noisy quantum processor
Authors:
Kevin J. Sung,
Marko J. Rančić,
Olivia T. Lanes,
Nicholas T. Bronn
Abstract:
The simulation of systems of interacting fermions is one of the most anticipated applications of quantum computers. The most interesting simulations will require a fault-tolerant quantum computer, and building such a device remains a long-term goal. However, the capabilities of existing noisy quantum processors have steadily improved, sparking an interest in running simulations that, while not nec…
▽ More
The simulation of systems of interacting fermions is one of the most anticipated applications of quantum computers. The most interesting simulations will require a fault-tolerant quantum computer, and building such a device remains a long-term goal. However, the capabilities of existing noisy quantum processors have steadily improved, sparking an interest in running simulations that, while not necessarily classically intractable, may serve as device benchmarks and help elucidate the challenges to achieving practical applications on near-term devices. Systems of non-interacting fermions are ideally suited to serve these purposes. While they display rich physics and generate highly entangled states when simulated on a quantum processor, their classical tractability enables experimental results to be verified even at large system sizes that would typically defy classical simulation. In this work, we use a noisy superconducting quantum processor to prepare Majorana zero modes as eigenstates of the Kitaev chain Hamiltonian, a model of non-interacting fermions. Our work builds on previous experiments with non-interacting fermionic systems. Previous work demonstrated error mitigation techniques applicable to the special case of Slater determinants. Here, we show how to extend these techniques to the case of general fermionic Gaussian states, and demonstrate them by preparing Majorana zero modes on systems of up to 7 qubits.
△ Less
Submitted 24 January, 2023; v1 submitted 1 June, 2022;
originally announced June 2022.
-
Observation of separated dynamics of charge and spin in the Fermi-Hubbard model
Authors:
Frank Arute,
Kunal Arya,
Ryan Babbush,
Dave Bacon,
Joseph C. Bardin,
Rami Barends,
Andreas Bengtsson,
Sergio Boixo,
Michael Broughton,
Bob B. Buckley,
David A. Buell,
Brian Burkett,
Nicholas Bushnell,
Yu Chen,
Zijun Chen,
Yu-An Chen,
Ben Chiaro,
Roberto Collins,
Stephen J. Cotton,
William Courtney,
Sean Demura,
Alan Derk,
Andrew Dunsworth,
Daniel Eppens,
Thomas Eckl
, et al. (74 additional authors not shown)
Abstract:
Strongly correlated quantum systems give rise to many exotic physical phenomena, including high-temperature superconductivity. Simulating these systems on quantum computers may avoid the prohibitively high computational cost incurred in classical approaches. However, systematic errors and decoherence effects presented in current quantum devices make it difficult to achieve this. Here, we simulate…
▽ More
Strongly correlated quantum systems give rise to many exotic physical phenomena, including high-temperature superconductivity. Simulating these systems on quantum computers may avoid the prohibitively high computational cost incurred in classical approaches. However, systematic errors and decoherence effects presented in current quantum devices make it difficult to achieve this. Here, we simulate the dynamics of the one-dimensional Fermi-Hubbard model using 16 qubits on a digital superconducting quantum processor. We observe separations in the spreading velocities of charge and spin densities in the highly excited regime, a regime that is beyond the conventional quasiparticle picture. To minimize systematic errors, we introduce an accurate gate calibration procedure that is fast enough to capture temporal drifts of the gate parameters. We also employ a sequence of error-mitigation techniques to reduce decoherence effects and residual systematic errors. These procedures allow us to simulate the time evolution of the model faithfully despite having over 600 two-qubit gates in our circuits. Our experiment charts a path to practical quantum simulation of strongly correlated phenomena using available quantum devices.
△ Less
Submitted 15 October, 2020;
originally announced October 2020.
-
Using models to improve optimizers for variational quantum algorithms
Authors:
Kevin J. Sung,
Jiahao Yao,
Matthew P. Harrigan,
Nicholas C. Rubin,
Zhang Jiang,
Lin Lin,
Ryan Babbush,
Jarrod R. McClean
Abstract:
Variational quantum algorithms are a leading candidate for early applications on noisy intermediate-scale quantum computers. These algorithms depend on a classical optimization outer-loop that minimizes some function of a parameterized quantum circuit. In practice, finite sampling error and gate errors make this a stochastic optimization with unique challenges that must be addressed at the level o…
▽ More
Variational quantum algorithms are a leading candidate for early applications on noisy intermediate-scale quantum computers. These algorithms depend on a classical optimization outer-loop that minimizes some function of a parameterized quantum circuit. In practice, finite sampling error and gate errors make this a stochastic optimization with unique challenges that must be addressed at the level of the optimizer. The sharp trade-off between precision and sampling time in conjunction with experimental constraints necessitates the development of new optimization strategies to minimize overall wall clock time in this setting. In this work, we introduce two optimization methods and numerically compare their performance with common methods in use today. The methods are surrogate model-based algorithms designed to improve reuse of collected data. They do so by utilizing a least-squares quadratic fit of sampled function values within a moving trusted region to estimate the gradient or a policy gradient. To make fair comparisons between optimization methods, we develop experimentally relevant cost models designed to balance efficiency in testing and accuracy with respect to cloud quantum computing systems. The results here underscore the need to both use relevant cost models and optimize hyperparameters of existing optimization methods for competitive performance. The methods introduced here have several practical advantages in realistic experimental settings, and we have used one of them successfully in a separately published experiment on Google's Sycamore device.
△ Less
Submitted 11 August, 2020; v1 submitted 22 May, 2020;
originally announced May 2020.
-
Creating and manipulating a Laughlin-type $ν=1/3$ fractional quantum Hall state on a quantum computer with linear depth circuits
Authors:
Armin Rahmani,
Kevin J. Sung,
Harald Putterman,
Pedram Roushan,
Pouyan Ghaemi,
Zhang Jiang
Abstract:
Here we present an efficient quantum algorithm to generate an equivalent many-body state to Laughlin's $ν=1/3$ fractional quantum Hall state on a digitized quantum computer. Our algorithm only uses quantum gates acting on neighboring qubits in a quasi-one-dimensional setting, and its circuit depth is linear in the number of qubits, i.e., the number of Landau orbitals in the second quantized pictur…
▽ More
Here we present an efficient quantum algorithm to generate an equivalent many-body state to Laughlin's $ν=1/3$ fractional quantum Hall state on a digitized quantum computer. Our algorithm only uses quantum gates acting on neighboring qubits in a quasi-one-dimensional setting, and its circuit depth is linear in the number of qubits, i.e., the number of Landau orbitals in the second quantized picture. We identify correlation functions that serve as signatures of the Laughlin state and discuss how to obtain them on a quantum computer. We also discuss a generalization of the algorithm for creating quasiparticles in the Laughlin state. This paves the way for several important studies, including quantum simulation of nonequilibrium dynamics and braiding of quasiparticles in quantum Hall states.
△ Less
Submitted 7 November, 2020; v1 submitted 5 May, 2020;
originally announced May 2020.
-
Quantum Approximate Optimization of Non-Planar Graph Problems on a Planar Superconducting Processor
Authors:
Matthew P. Harrigan,
Kevin J. Sung,
Matthew Neeley,
Kevin J. Satzinger,
Frank Arute,
Kunal Arya,
Juan Atalaya,
Joseph C. Bardin,
Rami Barends,
Sergio Boixo,
Michael Broughton,
Bob B. Buckley,
David A. Buell,
Brian Burkett,
Nicholas Bushnell,
Yu Chen,
Zijun Chen,
Ben Chiaro,
Roberto Collins,
William Courtney,
Sean Demura,
Andrew Dunsworth,
Daniel Eppens,
Austin Fowler,
Brooks Foxen
, et al. (61 additional authors not shown)
Abstract:
We demonstrate the application of the Google Sycamore superconducting qubit quantum processor to combinatorial optimization problems with the quantum approximate optimization algorithm (QAOA). Like past QAOA experiments, we study performance for problems defined on the (planar) connectivity graph of our hardware; however, we also apply the QAOA to the Sherrington-Kirkpatrick model and MaxCut, both…
▽ More
We demonstrate the application of the Google Sycamore superconducting qubit quantum processor to combinatorial optimization problems with the quantum approximate optimization algorithm (QAOA). Like past QAOA experiments, we study performance for problems defined on the (planar) connectivity graph of our hardware; however, we also apply the QAOA to the Sherrington-Kirkpatrick model and MaxCut, both high dimensional graph problems for which the QAOA requires significant compilation. Experimental scans of the QAOA energy landscape show good agreement with theory across even the largest instances studied (23 qubits) and we are able to perform variational optimization successfully. For problems defined on our hardware graph we obtain an approximation ratio that is independent of problem size and observe, for the first time, that performance increases with circuit depth. For problems requiring compilation, performance decreases with problem size but still provides an advantage over random guessing for circuits involving several thousand gates. This behavior highlights the challenge of using near-term quantum computers to optimize problems on graphs differing from hardware connectivity. As these graphs are more representative of real world instances, our results advocate for more emphasis on such problems in the developing tradition of using the QAOA as a holistic, device-level benchmark of quantum processors.
△ Less
Submitted 30 January, 2021; v1 submitted 8 April, 2020;
originally announced April 2020.
-
Hartree-Fock on a superconducting qubit quantum computer
Authors:
Frank Arute,
Kunal Arya,
Ryan Babbush,
Dave Bacon,
Joseph C. Bardin,
Rami Barends,
Sergio Boixo,
Michael Broughton,
Bob B. Buckley,
David A. Buell,
Brian Burkett,
Nicholas Bushnell,
Yu Chen,
Zijun Chen,
Benjamin Chiaro,
Roberto Collins,
William Courtney,
Sean Demura,
Andrew Dunsworth,
Daniel Eppens,
Edward Farhi,
Austin Fowler,
Brooks Foxen,
Craig Gidney,
Marissa Giustina
, et al. (57 additional authors not shown)
Abstract:
As the search continues for useful applications of noisy intermediate scale quantum devices, variational simulations of fermionic systems remain one of the most promising directions. Here, we perform a series of quantum simulations of chemistry the largest of which involved a dozen qubits, 78 two-qubit gates, and 114 one-qubit gates. We model the binding energy of ${\rm H}_6$, ${\rm H}_8$,…
▽ More
As the search continues for useful applications of noisy intermediate scale quantum devices, variational simulations of fermionic systems remain one of the most promising directions. Here, we perform a series of quantum simulations of chemistry the largest of which involved a dozen qubits, 78 two-qubit gates, and 114 one-qubit gates. We model the binding energy of ${\rm H}_6$, ${\rm H}_8$, ${\rm H}_{10}$ and ${\rm H}_{12}$ chains as well as the isomerization of diazene. We also demonstrate error-mitigation strategies based on $N$-representability which dramatically improve the effective fidelity of our experiments. Our parameterized ansatz circuits realize the Givens rotation approach to non-interacting fermion evolution, which we variationally optimize to prepare the Hartree-Fock wavefunction. This ubiquitous algorithmic primitive corresponds to a rotation of the orbital basis and is required by many proposals for correlated simulations of molecules and Hubbard models. Because non-interacting fermion evolutions are classically tractable to simulate, yet still generate highly entangled states over the computational basis, we use these experiments to benchmark the performance of our hardware while establishing a foundation for scaling up more complex correlated quantum simulations of chemistry.
△ Less
Submitted 18 September, 2020; v1 submitted 8 April, 2020;
originally announced April 2020.
-
Supplementary information for "Quantum supremacy using a programmable superconducting processor"
Authors:
Frank Arute,
Kunal Arya,
Ryan Babbush,
Dave Bacon,
Joseph C. Bardin,
Rami Barends,
Rupak Biswas,
Sergio Boixo,
Fernando G. S. L. Brandao,
David A. Buell,
Brian Burkett,
Yu Chen,
Zijun Chen,
Ben Chiaro,
Roberto Collins,
William Courtney,
Andrew Dunsworth,
Edward Farhi,
Brooks Foxen,
Austin Fowler,
Craig Gidney,
Marissa Giustina,
Rob Graff,
Keith Guerin,
Steve Habegger
, et al. (52 additional authors not shown)
Abstract:
This is an updated version of supplementary information to accompany "Quantum supremacy using a programmable superconducting processor", an article published in the October 24, 2019 issue of Nature. The main article is freely available at https://www.nature.com/articles/s41586-019-1666-5. Summary of changes since arXiv:1910.11333v1 (submitted 23 Oct 2019): added URL for qFlex source code; added Er…
▽ More
This is an updated version of supplementary information to accompany "Quantum supremacy using a programmable superconducting processor", an article published in the October 24, 2019 issue of Nature. The main article is freely available at https://www.nature.com/articles/s41586-019-1666-5. Summary of changes since arXiv:1910.11333v1 (submitted 23 Oct 2019): added URL for qFlex source code; added Erratum section; added Figure S41 comparing statistical and total uncertainty for log and linear XEB; new References [1,65]; miscellaneous updates for clarity and style consistency; miscellaneous typographical and formatting corrections.
△ Less
Submitted 28 December, 2019; v1 submitted 23 October, 2019;
originally announced October 2019.
-
Learning to learn with quantum neural networks via classical neural networks
Authors:
Guillaume Verdon,
Michael Broughton,
Jarrod R. McClean,
Kevin J. Sung,
Ryan Babbush,
Zhang Jiang,
Hartmut Neven,
Masoud Mohseni
Abstract:
Quantum Neural Networks (QNNs) are a promising variational learning paradigm with applications to near-term quantum processors, however they still face some significant challenges. One such challenge is finding good parameter initialization heuristics that ensure rapid and consistent convergence to local minima of the parameterized quantum circuit landscape. In this work, we train classical neural…
▽ More
Quantum Neural Networks (QNNs) are a promising variational learning paradigm with applications to near-term quantum processors, however they still face some significant challenges. One such challenge is finding good parameter initialization heuristics that ensure rapid and consistent convergence to local minima of the parameterized quantum circuit landscape. In this work, we train classical neural networks to assist in the quantum learning process, also know as meta-learning, to rapidly find approximate optima in the parameter landscape for several classes of quantum variational algorithms. Specifically, we train classical recurrent neural networks to find approximately optimal parameters within a small number of queries of the cost function for the Quantum Approximate Optimization Algorithm (QAOA) for MaxCut, QAOA for Sherrington-Kirkpatrick Ising model, and for a Variational Quantum Eigensolver for the Hubbard model. By initializing other optimizers at parameter values suggested by the classical neural network, we demonstrate a significant improvement in the total number of optimization iterations required to reach a given accuracy. We further demonstrate that the optimization strategies learned by the neural network generalize well across a range of problem instance sizes. This opens up the possibility of training on small, classically simulatable problem instances, in order to initialize larger, classically intractably simulatable problem instances on quantum devices, thereby significantly reducing the number of required quantum-classical optimization iterations.
△ Less
Submitted 11 July, 2019;
originally announced July 2019.
-
Quantum algorithms to simulate many-body physics of correlated fermions
Authors:
Zhang Jiang,
Kevin J. Sung,
Kostyantyn Kechedzhi,
Vadim N. Smelyanskiy,
Sergio Boixo
Abstract:
Simulating strongly correlated fermionic systems is notoriously hard on classical computers. An alternative approach, as proposed by Feynman, is to use a quantum computer. Here, we discuss quantum simulation of strongly correlated fermionic systems. We focus specifically on 2D and linear geometry with nearest neighbor qubit-qubit couplings, typical for superconducting transmon qubit arrays. We imp…
▽ More
Simulating strongly correlated fermionic systems is notoriously hard on classical computers. An alternative approach, as proposed by Feynman, is to use a quantum computer. Here, we discuss quantum simulation of strongly correlated fermionic systems. We focus specifically on 2D and linear geometry with nearest neighbor qubit-qubit couplings, typical for superconducting transmon qubit arrays. We improve an existing algorithm to prepare an arbitrary Slater determinant by exploiting a unitary symmetry. We also present a quantum algorithm to prepare an arbitrary fermionic Gaussian state with $O(N^2)$ gates and $O(N)$ circuit depth. Both algorithms are optimal in the sense that the numbers of parameters in the quantum circuits are equal to those to describe the quantum states. Furthermore, we propose an algorithm to implement the 2-dimensional (2D) fermionic Fourier transformation on a 2D qubit array with only $O(N^{1.5})$ gates and $O(\sqrt N)$ circuit depth, which is the minimum depth required for quantum information to travel across the qubit array. We also present methods to simulate each time step in the evolution of the 2D Fermi-Hubbard model---again on a 2D qubit array---with $O(N)$ gates and $O(\sqrt N)$ circuit depth. Finally, we discuss how these algorithms can be used to determine the ground state properties and phase diagrams of strongly correlated quantum systems using the Hubbard model as an example.
△ Less
Submitted 30 April, 2018; v1 submitted 14 November, 2017;
originally announced November 2017.
-
OpenFermion: The Electronic Structure Package for Quantum Computers
Authors:
Jarrod R. McClean,
Kevin J. Sung,
Ian D. Kivlichan,
Yudong Cao,
Chengyu Dai,
E. Schuyler Fried,
Craig Gidney,
Brendan Gimby,
Pranav Gokhale,
Thomas Häner,
Tarini Hardikar,
Vojtěch Havlíček,
Oscar Higgott,
Cupjin Huang,
Josh Izaac,
Zhang Jiang,
Xinle Liu,
Sam McArdle,
Matthew Neeley,
Thomas O'Brien,
Bryan O'Gorman,
Isil Ozfidan,
Maxwell D. Radin,
Jhonathan Romero,
Nicholas Rubin
, et al. (10 additional authors not shown)
Abstract:
Quantum simulation of chemistry and materials is predicted to be an important application for both near-term and fault-tolerant quantum devices. However, at present, developing and studying algorithms for these problems can be difficult due to the prohibitive amount of domain knowledge required in both the area of chemistry and quantum algorithms. To help bridge this gap and open the field to more…
▽ More
Quantum simulation of chemistry and materials is predicted to be an important application for both near-term and fault-tolerant quantum devices. However, at present, developing and studying algorithms for these problems can be difficult due to the prohibitive amount of domain knowledge required in both the area of chemistry and quantum algorithms. To help bridge this gap and open the field to more researchers, we have developed the OpenFermion software package (www.openfermion.org). OpenFermion is an open-source software library written largely in Python under an Apache 2.0 license, aimed at enabling the simulation of fermionic models and quantum chemistry problems on quantum hardware. Beginning with an interface to common electronic structure packages, it simplifies the translation between a molecular specification and a quantum circuit for solving or studying the electronic structure problem on a quantum computer, minimizing the amount of domain expertise required to enter the field. The package is designed to be extensible and robust, maintaining high software standards in documentation and testing. This release paper outlines the key motivations behind design choices in OpenFermion and discusses some basic OpenFermion functionality which we believe will aid the community in the development of better quantum algorithms and tools for this exciting area of research.
△ Less
Submitted 27 February, 2019; v1 submitted 20 October, 2017;
originally announced October 2017.