Also known as QAOA · VQE · variational quantum algorithm · digital-analog · quantum annealing
A variational quantum algorithm (VQA) runs a parametrized evolution \(U(\boldsymbol\theta)\) on a quantum device,
measures a cost such as an energy \(C(\boldsymbol\theta) = \langle\psi_0|U^\dagger(\boldsymbol\theta)\,H\,
U(\boldsymbol\theta)|\psi_0\rangle\), and updates the parameters with a classical optimizer
[1, 2]. The best-known examples are the variational quantum eigensolver (VQE)
[3] and the quantum approximate optimization algorithm (QAOA) [4]. Read as control
problems, the parameters are controls, the cost is a figure of merit, and the hybrid loop is
closed-loop optimization on hardware. VQAs and optimal control are variational optimization at two
levels of abstraction, circuits and pulses. Results move between them in both directions [5].
where \(H_C\) encodes the cost function and \(H_M\) is a “mixer” [4]. This is the evolution under
\(H(t) = u(t) H_C + [1-u(t)] H_M\) with a control \(u(t) \in \{0, 1\}\) that alternates between its two values in \(2p\)
segments, whose durations \(\gamma_k, \beta_k\) are the parameters. In other words, it is bang-bang control. Three control-theoretic results frame it:
Optimal protocols. Pontryagin’s principle has been used to argue that the optimal protocols at fixed total time
are bang-bang [6]. A later analysis found that generic optimal protocols are bang–anneal–bang, with a
smooth annealing segment between pulsed ends [7]. See Pontryagin's maximum principle.
Annealing and beyond. A smooth \(u(t)\) that varies slowly is quantum annealing [8] and
adiabatic quantum computation [9], which QAOA can approximate at large depth. Optimized QAOA can also
learn non-adiabatic mechanisms that get around the small spectral gaps where annealing fails [10].
Universality. Whether alternating two fixed Hamiltonians generates every unitary is a
Quantum controllability question. For one-dimensional cost Hamiltonians with nearest-neighbour \(ZZ\) terms, the precise
conditions for universality have been proven [11].
Landscapes and trainability
Barren plateaus. For a wide class of randomly initialized circuits, the gradient is exponentially small in the
number of qubits [12]. The dynamical Lie algebra of the circuit’s generators, the same object that
decides controllability, controls the variance of the cost [13]. A general Lie-algebraic theory
gives this variance exactly for deep circuits and unifies the known sources of barren plateaus: expressiveness,
entangled inputs, local observables and noise [14]. Full controllability and trainability therefore
pull in opposite directions.
Overparametrization. Above a critical number of parameters, bounded by the dimension of that Lie algebra,
spurious local minima start to disappear and training becomes much easier [15].
Problem-inspired ansätze such as the Hamiltonian variational ansatz show mild barren plateaus and become almost trap
free when overparametrized [16]. This is the circuit counterpart of the trap-free
Control landscapes of optimal control.
From circuits to pulses, and back
Gradients from measurements. The parameter-shift rule obtains exact gradients from the same circuit run at
shifted parameters [17, 18]. A stochastic version extends it to general
parametrized evolutions, such as pulse amplitudes [19]. These are hardware-ready gradient
estimators of the kind Closed-loop & in-situ optimization control needs.
Optimizing pulses directly. In ctrl-VQE, the pulse shapes of the device, not gate angles, are the variational
parameters. In simulations of molecules on transmons, this shortened the state-preparation time by roughly three
orders of magnitude [20]. See Pulse-level compilation.
Ansätze from control. Hardware-efficient ansätze use the device’s native interactions [21].
Ansätze inspired by optimal control add symmetry-breaking terms that improve convergence, for example for the
Fermi–Hubbard model [22]. Counterdiabatic terms enrich QAOA layers [23] (see
Shortcuts to adiabaticity). Digital–analog schemes combine single-qubit gates with an always-on analog Ising
interaction [24].
Optimal control brings controllability and Lie-algebra tools, landscape theory, Quantum speed limits and
robustness, and the lesson that enough control resources (parameters, time) make
optimization easy [5]. Variational algorithms bring gradient estimation from measurements,
optimizers that cope with shot noise, and a precise theory of when the cost concentrates. Both use the same
mathematics for problems on many-body Hamiltonians.
Limitations
Every cost evaluation is a finite-shot estimate. Gradients below the shot-noise level cannot be resolved, which
turns barren plateaus into a practical limit.
Hardware noise flattens landscapes as well and biases the cost.
Pulse-level ansätze optimized in simulation need an accurate device model. Without one, they must be optimized on
the hardware, with every evaluation costing measurements.
Key papers
Edward Farhi, Jeffrey Goldstone and Sam GutmannA Quantum Approximate Optimization Algorithmpreprint (2014)arXiv:1411.4028
Introduces QAOA — alternating cost and mixer Hamiltonians, with the switching times optimized classically.
Alberto Peruzzo, Jarrod McClean, Peter Shadbolt, Man-Hong Yung, Xiao-Qi Zhou, et al.A variational eigenvalue solver on a photonic quantum processorNat Commun5, 4213 (2014)doi
Introduces the variational quantum eigensolver, demonstrated on a photonic processor.
Applies Pontryagin's principle to variational algorithms and argues for bang-bang (QAOA-like) protocols.
Alicia B. Magann, Christian Arenz, Matthew D. Grace, Tak-San Ho, Robert L. Kosut, et al.From Pulses to Circuits and Back Again: A Quantum Optimal Control Perspective on Variational Quantum AlgorithmsPRX Quantum2, 010101 (2021)doi
Perspective that puts variational algorithms and optimal control in one hierarchy of variational optimization, from circuits to pulses.
Martin Larocca, Piotr Czarnik, Kunal Sharma, Gopikrishnan Muraleedharan, Patrick J. Coles and M. CerezoDiagnosing Barren Plateaus with Tools from Quantum Optimal ControlQuantum6, 824 (2022)doi
Uses the dynamical Lie algebra, a controllability tool, to diagnose barren plateaus.
M. Cerezo, Andrew Arrasmith, Ryan Babbush, Simon C. Benjamin, Suguru Endo, et al.Variational quantum algorithmsNat Rev Phys3, 625-644 (2021)doi
Broad review of variational quantum algorithms, their applications and their trainability problems.
Software
PennyLane — differentiable variational circuits, including pulse-level programs with parametrized Hamiltonians
References
M. Cerezo, Andrew Arrasmith, Ryan Babbush, Simon C. Benjamin, Suguru Endo, et al.Variational quantum algorithmsNat Rev Phys3, 625-644 (2021)doi
Kishor Bharti, Alba Cervera-Lierta, Thi Ha Kyaw, Tobias Haug, Sumner Alperin-Lea, et al.Noisy intermediate-scale quantum algorithmsRev. Mod. Phys.94, 015004 (2022)doi
Alberto Peruzzo, Jarrod McClean, Peter Shadbolt, Man-Hong Yung, Xiao-Qi Zhou, et al.A variational eigenvalue solver on a photonic quantum processorNat Commun5, 4213 (2014)doi
Edward Farhi, Jeffrey Goldstone and Sam GutmannA Quantum Approximate Optimization Algorithmpreprint (2014)arXiv:1411.4028
Alicia B. Magann, Christian Arenz, Matthew D. Grace, Tak-San Ho, Robert L. Kosut, et al.From Pulses to Circuits and Back Again: A Quantum Optimal Control Perspective on Variational Quantum AlgorithmsPRX Quantum2, 010101 (2021)doi
Lucas T. Brady, Christopher L. Baldwin, Aniruddha Bapat, Yaroslav Kharkov and Alexey V. GorshkovOptimal Protocols in Quantum Annealing and Quantum Approximate Optimization Algorithm ProblemsPhys. Rev. Lett.126, 070505 (2021)doi
Tadashi Kadowaki and Hidetoshi NishimoriQuantum annealing in the transverse Ising modelPhys. Rev. E58, 5355-5363 (1998)doi
Tameem Albash and Daniel A. LidarAdiabatic quantum computationRev. Mod. Phys.90, 015002 (2018)doi
Leo Zhou, Sheng-Tao Wang, Soonwon Choi, Hannes Pichler and Mikhail D. LukinQuantum Approximate Optimization Algorithm: Performance, Mechanism, and Implementation on Near-Term DevicesPhys. Rev. X10, 021067 (2020)doi
M. E. S. Morales, J. D. Biamonte and Z. ZimborásOn the universality of the quantum approximate optimization algorithmQuantum Inf Process19, 291 (2020)doi
Jarrod R. McClean, Sergio Boixo, Vadim N. Smelyanskiy, Ryan Babbush and Hartmut NevenBarren plateaus in quantum neural network training landscapesNat Commun9, 4812 (2018)doi
Martin Larocca, Piotr Czarnik, Kunal Sharma, Gopikrishnan Muraleedharan, Patrick J. Coles and M. CerezoDiagnosing Barren Plateaus with Tools from Quantum Optimal ControlQuantum6, 824 (2022)doi
Michael Ragone, Bojko N. Bakalov, Frédéric Sauvage, Alexander F. Kemper, Carlos Ortiz Marrero, et al.A Lie algebraic theory of barren plateaus for deep parameterized quantum circuitsNat Commun15, 7172 (2024)doi
Martín Larocca, Nathan Ju, Diego García-Martín, Patrick J. Coles and Marco CerezoTheory of overparametrization in quantum neural networksNat Comput Sci3, 542-551 (2023)doi
Roeland Wiersema, Cunlu Zhou, Yvette de Sereville, Juan Felipe Carrasquilla, Yong Baek Kim and Henry YuenExploring Entanglement and Optimization within the Hamiltonian Variational AnsatzPRX Quantum1, 020319 (2020)doi
K. Mitarai, M. Negoro, M. Kitagawa and K. FujiiQuantum circuit learningPhys. Rev. A98, 032309 (2018)doi
Maria Schuld, Ville Bergholm, Christian Gogolin, Josh Izaac and Nathan KilloranEvaluating analytic gradients on quantum hardwarePhys. Rev. A99, 032331 (2019)doi
Leonardo Banchi and Gavin E. CrooksMeasuring Analytic Gradients of General Quantum Evolution with the Stochastic Parameter Shift RuleQuantum5, 386 (2021)doi
Oinam Romesh Meitei, Bryan T. Gard, George S. Barron, David P. Pappas, Sophia E. Economou, et al.Gate-free state preparation for fast variational quantum eigensolver simulationsnpj Quantum Inf7, 155 (2021)doi
Abhinav Kandala, Antonio Mezzacapo, Kristan Temme, Maika Takita, Markus Brink, et al.Hardware-efficient variational quantum eigensolver for small molecules and quantum magnetsNature549, 242-246 (2017)doi
Alexandre Choquette, Agustin Di Paolo, Panagiotis Kl. Barkoutsos, David Sénéchal, Ivano Tavernelli and Alexandre BlaisQuantum-optimal-control-inspired ansatz for variational quantum algorithmsPhys. Rev. Research3, 023092 (2021)doi
P. Chandarana, N. N. Hegade, K. Paul, F. Albarrán-Arriagada, E. Solano, et al.Digitized-counterdiabatic quantum approximate optimization algorithmPhys. Rev. Research4, 013141 (2022)doi
Adrian Parra-Rodriguez, Pavel Lougovski, Lucas Lamata, Enrique Solano and Mikel SanzDigital-analog quantum computationPhys. Rev. A101, 022305 (2020)doi
Alicia B. Magann, Kenneth M. Rudinger, Matthew D. Grace and Mohan SarovarFeedback-Based Quantum OptimizationPhys. Rev. Lett.129, 250502 (2022)doi
Alicia B. Magann, Kenneth M. Rudinger, Matthew D. Grace and Mohan SarovarLyapunov-control-inspired strategies for quantum combinatorial optimizationPhys. Rev. A106, 062414 (2022)doi
Recent papers
304 tracked papers since 2000 are classified
under this topic. See all ·
Feed
Preparing quantum states in bosonic qudits with SNAP--Displacement sequences requires both optimizing the gate parameters and deciding which Fock levels should be controlled by each SNAP gate. A fully parameterized ansatz assigns an independent phase to every Fock level in every gate, leading to a number of phases that grows steeply. This constitutes a problem for short to medium term systems as controlling many phases requires complex control pulses which are subject to increased control errors. Previous work reduced this cost by restricting the trainable phases to predefined sparse patterns, such as a diagonal or a grid, which however may not be optimal for a specific target state. In this work, we address this limitation by introducing a method that learns which phases are needed together with their values, using a differentiable mask. A sparsity penalty encourages the removal of unnecessary phases, while temperature annealing gradually turns the masks into binary decisions. We compare the proposed method with fully parameterized and predefined sparse ansatzes on Fourier and Haar-random target states. The evaluation considers both the minimum number of phases required to reach different infidelity thresholds and the overall trade-off between fidelity and phases required, measured through Pareto frontiers and hypervolume. The results show that learning the location of the phases together with their values provides competitive or improved fidelity--resource trade-offs, although the advantage depends on the target state and qudit dimension.
We introduce an inertial counterdiabatic (CD) protocol that reformulates shortcut-to-adiabaticity design as a moving-frame diagonalization problem. The method employs the Magnus expansion to construct a variational propagator generated by an accessible CD ansatz. In the new rotating frame, the controls are optimized to suppress the off-diagonal Hamiltonian terms, and the CD terms are defined to counteract the corresponding inertial contribution. The propagator generated with the Magnus expansion naturally exploits additional control directions generated through its nested commutator, while the sequential design of the propagator introduces temporal nonlocality into the optimization. We benchmark the method on a two-qubit system and on the one-dimensional quantum Ising chain, where it consistently outperforms variational CD protocols based on action minimization. The largest improvements are observed near gap-closing regions, demonstrating the advantage of combining temporal and algebraic nonlocality for suppressing diabatic excitations.
Preparing stabilizer-code states with shallow circuits is an important primitive for near-term quantum error correction in superconducting circuits, where logical-state initialization requires four-body stabilizer correlations from native one- and two-body controls. We introduce a digital-analog method for the preparation of surface-code ground-state manifolds. The method exploits a structural property of the adiabatic gauge potential, mapping stabilizer-local counterdiabatic terms to digital-analog blocks on the superconducting layout using fixed-angle two-qubit XY gates, single-qubit rotations, and analog exchange-interaction evolutions. The dressed blocks generate higher-order operators used as the variational ansatz to determine the counterdiabatic terms within the Sels-Polkovnikov adiabatic-gauge-potential framework. We benchmark checkerboard plaquette-stabilizer grids up to 4 x 4; the square-grid instances correspond to rotated surface codes, while the rectangular cases probe the scaling of the method. For all grids, the proposed method generates the dominant four-body counterdiabatic basis and, at short evolution times, substantially improves ground-state preparation compared with bare adiabatic evolution. We further generalize the construction to n x m stabilizer lattices and show that digital-analog synthesis can reduce the entangling depth by an order of magnitude. Together, these results establish a direct connection between stabilizer geometry, the structure of the counterdiabatic gauge potential, and digital-analog control for measurement-free preparation of stabilizer-code states on near-term superconducting architectures.
Frequency crowding, the close spacing of qubit transitions, limits superconducting quantum processors through crosstalk, leakage, and spectator dependent errors. We introduce a Hamiltonian based framework that maps two- and three-body collision mechanisms onto a Max-k-Cut problem on qubit frequency allocation. In order to make larger lattices tractable, we introduce a cluster-stitching protocol that decomposes larger lattices into overlapping subgraphs and reconciles their local frequency assignments into globally consistent configurations. We tested our proposed framework across multiple qubit architectures to suppress dominant collision penalties by revealing a hardware residual crosstalk. We benchmark digitized counterdiabatic QAOA (DC-QAOA), standard QAOA, and simulated annealing. Counterdiabatic terms improve QAOA optimization near steep collision penalties, while simulated annealing achieves the lowest costs on the tested instances. These results indicate that our hardware-aware optimization connects microscopic interaction physics to processor-level frequency planning. These results enable a systematic solver for benchmarking of scalable software hardware co design.
The Quantum Approximate Optimization Algorithm (QAOA), whose design is underpinned by the adiabatic theorem, is one of the leading variational quantum algorithms for preparing the ground states of gauge theories. A recently proposed variant, DC-QAOA, incorporates counterdiabatic driving to accelerate the adiabatic process, reducing both circuit depth and runtime. In this work, we present a systematic analysis, from both theoretical and computational perspectives, of DC-QAOA and related counterdiabatic variants for preparing the ground state of the Schwinger model, the paradigmatic (1+1)-dimensional lattice gauge theory. Benchmarking these methods against standard QAOA, we show that counterdiabatic driving improves ground-state preparation while lowering circuit depth requirements. Our results establish counterdiabatic protocols as a practical route towards efficient quantum simulation of gauge theories.