Variational Quantum Continuous Optimization:
a Cornerstone of Quantum Mathematical Analysis
Pablo Bermejo1, 2 and Rom´an Or´us1, 2, 3
1Multiverse Computing, Paseo de Miram´on 170, E-20014 San Sebasti´an, Spain
2Donostia International Physics Center, Paseo Manuel de Lardizabal 4, E-20018 San Sebasti´an, Spain
3Ikerbasque Foundation for Science, Maria Diaz de Haro 3, E-48013 Bilbao, Spain
Here we show how universal quantum computers based on the quantum circuit model can handle
mathematical analysis calculations for functions with continuous domains, without any digitaliza-
tion, and with remarkably few qubits. The basic building block of our approach is a variational
quantum circuit where each qubit encodes up to three continuous variables (two angles and one
radious in the Bloch sphere). By combining this encoding with quantum state tomography, a vari-
ational quantum circuit of nqubits can optimize functions of up to 3ncontinuous variables in an
analog way. We then explain how this quantum algorithm for continuous optimization is at the basis
of a whole toolbox for mathematical analysis on quantum computers. For instance, we show how to
use it to compute arbitrary series expansions such as, e.g., Fourier (harmonic) decompositions. In
turn, Fourier analysis allows us to implement essentially any task related to function calculus, includ-
ing the evaluation of multidimensional definite integrals, solving (systems of) differential equations,
and more. To prove the validity of our approach, we provide benchmarking calculations for many
of these use-cases implemented on a quantum computer simulator. The advantages with respect to
classical algorithms for mathematical analysis, as well as perspectives and possible extensions, are
also discussed.
I. INTRODUCTION
The field of quantum computing is evolving rapidly.
What looked like science-fiction a few years ago, is now
starting to become a reality. Thanks to recent experi-
mental breakthroughs, we now can enjoy the first com-
mercial quantum processors. Such machines are noisy
and have a limited number of qubits (i.e., they are Noisy
Intermediate-Scale Quantum - NISQ - devices [1]), and
are thus far from ideal. Following the trends in his-
tory, we should expect to have more powerful and noise-
resistant processors in the future. But, as of today, our
quantum computers are made of a small number of noisy
qubits.
Now, here comes a hard question: can we do some-
thing useful with these imperfect quantum machines?
The qualified answer to this question is yes. To put it in
perspective, think of the evolution of classical comput-
ers. One would never use a computer from the 1930’s
to develop, say, a speech recognition system. But those
machines were definitely useful at cracking some crypto-
graphic codes. The situation with current quantum com-
puters is, in a way, similar to that of classical computers
in the 30’s. Of course we cannot hope to use today’s
quantum processors for very complex tasks such as, say,
factoring a 617-digit number (as used in the RSA cyber-
security protocol). But they are certainly useful for other
tasks. And of course, what consumes mental energy is to
find those tasks, where today’s quantum computers can
begin to offer real, practical value. So, what can we al-
ready do with NISQ processors that is actually useful?
A good example where NISQ devices can start offering
value is optimization. By developing hybrid quantum-
classical solutions, present-day quantum annealers can
handle hard optimization problems in real scenarios [2–
10]. Another good example is that of machine learning,
where few-qubit universal quantum systems are, at least,
as powerful as deep neural networks for supervised and
unsupervised learning [11,12], thanks to the use of vari-
ational quantum optimization methods [13,14].
In this paper we follow the above philosophy, and show
how few-qubit quantum computers can already imple-
ment arbitrary multi-dimensional function calculus in a
remarkably efficient way. The basic building block of
our approach is a variational quantum algorithm to op-
timize functions with continuous domains. In our ap-
proach, each qubit in the quantum circuit encodes up
to three continuous variables in the parameters of the
Bloch sphere: two angles and one radius. By combining
this with single-qubit quantum tomography, a variational
optimization of the circuit parameters allows to find the
extreme values of multidimensional functions in a purely-
analog and fast way.
Based on the above, we then show how our quantum
optimization algorithm can be used to implement, essen-
tially, any type of calculus on functions with continuous
domains. As a first step, we show how to compute arbi-
trary series expansions of any function, and in particular
Fourier (harmonic) analysis, with the coefficients of the
expansions calculated via continuous quantum optimiza-
tion. And, then, with this tool at our disposal, we explain
how this allows us to do pretty much any calculation, in-
cluding definite integrals of multidimensional functions,
solving (systems of) differential equations, and more.
The quantum optimization algorithm presented here
is the cornerstone of an efficient quantum toolbox for
mathematical analysis on NISQ processors. To put it in
perspective, a quantum computer of 127 qubits, such as
the IBM-Q System One [15], could use our algorithms to
analyze 381-dimensional functions directly in the contin-
arXiv:2210.03136v1 [quant-ph] 6 Oct 2022