Articulo 3
Articulo 3
IonQ, Inc.
Abstract
Quantum computers are expected to contribute more efficient and accurate ways of modeling economic
processes. Quantum hardware is currently available at a relatively small scale, but effective algorithms are
limited by the number of logic gates that can be used, before noise from gate inaccuracies tends to dominate
results. Some theoretical algorithms that have been proposed and studied for years do not perform well
yet on quantum hardware in practice. This encourages the development of suitable alternative algorithms
that play similar roles in limited contexts.
This paper implements this strategy in the case of quantum counting, which is used as a component
for keeping track of position in a quantum walk, which is used as a model for simulating asset prices
over time. We introduce quantum approximate counting circuits that use far fewer 2-qubit entangling
gates than traditional quantum counting that relies on binary positional encoding. The robustness of these
circuits to noise is demonstrated.
We compare the results to price change distributions from stock indices, and compare the behavior of
quantum circuits with and without mid-measurement to trends in the housing market. The housing data
shows that low liquidity brings price volatility, as expected with the quantum models.
1
model a continuous time walk on a connected graph with 4 states. Stamatopoulos et al. (2020) used 3
superconducting qubits for option pricing. and Zhu et al. (2022b) used 6 trapped-ion qubits to perform
generative modeling for correlated stock prices.
The scale of such experiments has been particularly limited by qubit availability and quantum gate
accuracy. For example, the 3-qubit circuit of Stamatopoulos et al. (2020) was optimized down to 18 2-
qubit entangling gates and 33 single-qubit gates, but even with this small circuit, error rates in the results
ranged from 62% raw, to 21% using Richardson extrapolation for error-correction. This is not surprising,
since the accuracies of the single- and 2-qubit gates were estimated at 99.7% and 97.8% respectively, and
0.99733 × 0.97818 ≈ 0.587, so the compounded gate error rate is at least 40%.
A safe implementation strategy might be to wait for large-scale fault-tolerant quantum computers to
become available, but this runs the risk of missing opportunities in the meantime. Instead, researchers such
as Stamatopoulos et al. (2020) and Zhu et al. (2022b) try to use currently-available quantum hardware, and
ask whether implementations can be made robust enough to provide value sooner. In the current NISQ era
(Noisy Intermediate-Scale Quantum), the scarce resources include the number of qubits, and also, as seen
above, the number of gates, and especially the number of 2-qubit entangling gates. Circuits are sometimes
described in terms of width (number of qubits) and depth (number of dates, or layers of gates), and both
need to be minimized.
Quantum developers sometimes have many suggested designs to start from: quantum information
processing has been explored as an academic field for some decades, and established literature provides
many circuit recipes (Nielsen and Chuang, 2002). A natural strategy is to take such designs, consider their
NISQ era limitations, and see if there are alternatives that provide some of the same functionality using
fewer qubits or gates.
This paper develops some new examples of this approach, with the basic example of quantum counting.
The central novel contribution of the paper are the approximate quantum counting circuits, introduced in
Section 5. The motivation is that quantum counting is used as a component in the implementation of quan-
tum random walks, which are proposed as a model for stock prices, and also for beliefs about the future
value of stock prices, for the pricing of stock options. Beliefs are less exact than prices: it is not very impor-
tant to make sure that an estimate of $1,000 comes $1 after an estimate of $999 and $1 before an estimate of
$1,001; but it is important to make sure that these are all treated similarly, and that doubling any of them
gives something in the region of $2,000. The circuits proposed in this paper demonstrate such properties,
albeit approximately, but much more accurately than is currently possible on quantum computers that use
positional binary representations for numbers that strictly follow the axioms of arithmetic.
A distinct feature of quantum systems including quantum walks is that they behave differently when
they are measured, compared to when they are left to evolve dynamically. Such behavior has been demon-
strated with humans in psychology experiments (Kvam et al., 2015; Yearsley and Pothos, 2016), and is an
important feature in quantum cognition and economics (Busemeyer and Bruza, 2012; Orrell, 2020). Section
7 investigates the simulated behavior of the approximate counting circuits with mid-measurement, and
shows that they exhibit desirable behavior (in this case, that more frequent measurement tends to reduce
the chances of large changes). In a departure from many quantum cognition models, the macro effects of
beliefs on transactions and prices cannot be described as the decision of particular cognitive agent, and it
may be more appropriate to think of quantum models as representing the beliefs of whole groups of buyers
and sellers, and measurements as decisions observed by the whole market. This theme is considered as
part of introducing quantum walks in finance in Section 3, and the effects of measurement in Section 7,
especially with reference to the housing market.
To begin with, the next few sections review some of the basic quantum logic gates and how they are
put together into quantum circuits, the use of random walks and quantum walks in finance, and how these
2
<latexit sha1_base64="zJhMk4p1qYWXPXBZBWXd7jKPKs4=">AAACdXicbZBPb9NAEMUnLn9K+JcWbghpRaDilNqogh4juHAsEmkjxVEYryfuKuu1uzsuSq18Pb4D34ErcEOsnSLRNiNZ+um9N5r1S0qtHIfh906wdev2nbvb97r3Hzx89Li3s3vsispKGslCF3acoCOtDI1YsaZxaQnzRNNJsvjQ+CfnZJ0qzGdeljTNMTNqriSyl2a9L3FCmTI1q8WFTFciPvsq9kScIVM9XjXohZhM+i/R3bzAaLN6U37W64eDsB1xE6JL6A+fQjtHs96fOC1klZNhqdG5SRSWPK3RspKa/P3KUYlygRlNPBrMyU3rtomVeOWVVMwL6z/DolX/36gxd26ZJz6ZI5+6614jbvImFc8Pp7UyZcVk5PrQvNKCC9HUKlJlSbJeekBplX+rkKdoUbIvv3vlTFPM+jca0iqxaJf1WYWmdXxh0fV6bsLxm0H0dnDw6aA/fL9uDrbhGbyA1xDBOxjCRziCEUj4Bj/gJ/zq/A6eBy+DvXU06FzuPIErE+z/BRaWw1M=</latexit>
0 1
Pauli-X (NOT) X
1 0
1 1
<latexit sha1_base64="yYk4hGQTAux1voauwt2FW+84k5U=">AAACRXicbZDPbtNAEMbHLdA0pRAKNy4rIlBPkY0q4BiVS45FIm2lOIrG60m6ynrt7o6LgpXX6NNwBZ6Bh+CGemxZO0Wif0Za6afvm9HsfEmhleMw/BWsrT94+Gijtdneerz95Gnn2c6hy0sraShzndvjBB1pZWjIijUdF5YwSzQdJfOPtX90Rtap3HzmRUHjDGdGTZVE9tKkE8YJzZSpWM2/ynQp4tMv4o2IZ8hUDZY1eiEmk/7rmHS6YS9sStyF6Bq6/RfQ1MGkcxmnuSwzMiw1OjeKwoLHFVpWUtOyHZeOCpRznNHIo8GM3LhqLluK115JxTS3/hkWjfr/RIWZc4ss8Z0Z8om77dXifd6o5OmHcaVMUTIZuVo0LbXgXNQxiVRZkqwXHlBa5f8q5AlalOzDbN9YUwezOqMmrRKLdlGdlmgaxwcW3Y7nLhy+7UXvenuf9rr9/VVy0IKX8Ap2IYL30IcBHMAQJJzDN/gOP4Kfwe/gT3Cxal0Lrmeew40Krv4CtCGzyQ==</latexit>
Hadamard (H) H √1
2 1 −1
cos θ2 i sin θ2
RX rotation
<latexit sha1_base64="Bc4gbJR/MsEQg8ZiX2fvBRh/W+I=">AAACT3icbZBNb9NAEIbH4asNX6Fw47IiApVLZKMKUE8RXHosiLSR4igaryfJKuu1uzumSq38F34NV+ixv6S3irUTJNoy0kqP3ndWM/MmhVaOw/AiaN25e+/+g63t9sNHj5887TzbOXJ5aSUNZK5zO0zQkVaGBqxY07CwhFmi6ThZfK794+9kncrNN14WNM5wZtRUSWQvTTr7cUIzZSpWizOZrkR8cireiHiGTNXXyXA35jkxvl3VordiMunf3kmnG/bCpsRtiDbQ7b+Apg4nnas4zWWZkWGp0blRFBY8rtCykppW7bh0VKBc4IxGHg1m5MZVc+NKvPZKKqa59c+waNR/f1SYObfMEt+ZIc/dTa8W/+eNSp5+HFfKFCWTketB01ILzkUdmEiVJcl66QGlVX5XIedoUbKPtX1tTB3M+oyatEos2mV1UqJpHB9YdDOe23D0rhe97+192ev2P62Tgy14Ca9gFyL4AH04gEMYgIQf8BN+we/gPLgMrlqb1lawgedwrVrbfwC3+LWz</latexit>
RX (✓)
exp(−i θ2 X) i sin θ2 cos θ2
Figure 1: Some standard single-qubit gates and their corresponding matrices, which operate on the super-
position state α |0⟩ + β |1⟩ written as a column vector (α, β)T .
3
<latexit sha1_base64="a1DAmmYPemi+THLqKCCl9SHhfoA=">AAACXXicbZDBbtNAEIYnLoViSkmhEgcuKyIQp8iuKugxgkuPrUTaSnEUjdeTdJX12t0dg4KVV+JpeuHQ8iaItZNKtGWklX59/4xm509LrRxH0a9OsPFo8/GTrafhs+3nOy+6uy9PXVFZSUNZ6MKep+hIK0NDVqzpvLSEearpLJ1/afyzb2SdKsxXXpQ0znFm1FRJZI8m3aMkpZkyNav5D5ktRZhcfhfvRSLZ6jpeNtKDJLnljHZWr3GYkMluJyfdXtSP2hIPRbwWvcEetHU86f5JskJWORmWGp0bxVHJ4xotK6lpGSaVoxLlHGc08tJgTm5ctxcvxTtPMjEtrH+GRUv/nagxd26Rp74zR75w970G/s8bVTw9HNfKlBWTkatF00oLLkQTn8iUJcl64QVKq/xfhbxAi5J9yOGdNU0wqzMapVVq0S7qywpN6/jA4vvxPBSn+/34Y//g5KA3+LxKDrbgDbyFDxDDJxjAERzDECT8hCu4hpvO72Az2A52Vq1BZz3zCu5U8PovuKe35g==</latexit>
1 0 0 0
Controlled Not 0 1 0 0
(CNOT, CX) 0 0 0 1
0 0 1 0
Figure 2: The CNOT gate is a 2-qubit entangling gate, that acts upon the state α |00⟩ + β |01⟩ + γ |10⟩ + δ|11⟩.
In the standard basis, its behavior can be described as “performing a NOT operation on the target qubit if
the control qubit is in state |1⟩”.
<latexit sha1_base64="MEDC3vV1k2CPz1dBN5+E0IgBpdA=">AAACbnicbZDbbtNAEIY3boFiTqEgbqqKFRGoV5GNKsplBDe9LBJpi+IoGq8n7irrtbs7pgqrPBvPwQNwCzwBYm2nUg+MtKtP/z+zhz+tlLQURT96wcbmnbv3tu6HDx4+evyk/3T72Ja1ETgWpSrNaQoWldQ4JkkKTyuDUKQKT9LFx8Y/+YrGylJ/pmWF0wJyLedSAHlp1v+SpJhL7UguvolsxcPk/IK/4UkOhO5w1aAgo1zcoveS5LKl2wlM7tZmmKDOLo+a9QfRMGqL34Z4DYPRc9bW0az/N8lKUReoSSiwdhJHFU0dGJJC4SpMaosViAXkOPGooUA7dW0EK/7aKxmfl8YvTbxVr044KKxdFqnvLIDO7E2vEf/nTWqav586qauaUIvuonmtOJW8yZNn0qAgtfQAwkj/Vi7OwIAgn3p47ZommO4bDSmZGjBLd16Dbh0fWHwznttw/HYYvxvuf9ofjD50ybEttsNesT0WswM2YofsiI2ZYN/ZT/aL/e79CV4Eu8HLrjXorWeesWsV7P0D7yG93Q==</latexit>
Figure 3: Hadamard and CNOT gates in sequence make a quantum circuit that prepares the Bell state
√1 (|00⟩ + |11⟩).
2
distinguishes quantum computing algorithmically, because predicting the probability distributions that re-
sult from quantum operations with entanglement can become exponentially hard for classical computers.
In simpler terms, quantum computing is special because it offers special kinds of interference, not because
it offers special kinds of in-between-ness.
A quantum circuit consists of a register of qubits, and a sequence of logic gates that act on these qubits.
For example, the circuit in Figure 3 prepares the famous Bell state (named after physicist John Bell, whose
pioneering theorem motivated experiments that demonstrated real entanglement). It maps the input state
|00⟩ to the state √12 (|00⟩+|11⟩), which has the crucial ‘entangled’ behavior whereby if one qubit is measured
to be in the |0⟩ state, the other qubit must also be in the |0⟩ state, and vice versa.
Quantum circuits finish with measurements that record the 0 or 1 state of at least some of the qubits,
and output this as classical information. (Some platforms also support measuring qubits before the end of
a circuit.) The measurement outcomes are probabilistic, following the Born rule. A run of a single quantum
circuit including outputting one sample of measurements is typically called a shot. Circuits are usually
repeated several times, and the output counts from many individual shots are summarized into an output
distribution. The process of running all the shots and gathering the outcomes is typically called a job.
There are many standard gate recipes and equivalences. In particular, larger operators are often thought
of as distinct gates in their own right, an important example being the 3-qubit Toffoli gate of Figure 4. This
is like an extended CNOT gate — it has 2 control qubits instead of 1, and performs an X-rotation / NOT
operation on the target qubit if both the control qubits are in the |1⟩ state. The decomposition in Figure
4 shows that 5 CNOT gates are needed for each Toffoli gate. There are variants of this, but as a general
rule-of-thumb, the error-rate of a Toffoli gate will be at least 4 times the error-rate of the 2-qubit gates form
which it is assembled. Toffoli gates are particularly important for binary arithmetic, as seen in Section 5.
In the NISQ era, such considerations are pervasive: there is a ubiquitous tradeoff between circuit com-
plexity (the number of gates needed to execute a given algorithm) and expected circuit accuracy (the more
gates we use, the more inaccurate our results become).
4
<latexit sha1_base64="eY0tWFxVn+L74jRQRGFt9jB8VXY=">AAADrHiclVJdb9MwFHUTPkb46oA3XiwqEA9QJdPEeKzgZY9DardJTahuHCez6jiZfcNUov4N/hs/BoSTDilNJ8SuZOvonHt97Osbl1IY9P2fA8e9c/fe/b0H3sNHj588He4/OzVFpRmfsUIW+jwGw6VQfIYCJT8vNYc8lvwsXn5u9LNvXBtRqCmuSh7lkCmRCgZoqcX+4EcY80yoGsXyO0vW1Asvr+gbesPOUMv6YP2fYksELZEB8nq67rE2LQy93UOCfznsip3DEXRWdxy/hglkGdc9ccu4zTy+TfntDTsVHbveC5srcZX8/YfFcOSP/TboLgiuwWjygrRxshj+DpOCVTlXyCQYMw/8EqMaNAom+doLK8NLYEvI+NxCBTk3Ud0O0Jq+tkxC00LbpZC2bLeihtyYVR7bzBzwwvS1hrxJm1eYfoxqocoKuWIbo7SSFAvaTCNNhOYM5coCYFrYu1J2ARoY2pn1tmyaxmye0SApYg16VV9WoFrFNizot2cXnB6Mgw/jwy+Ho8mnTefIHnlJXpG3JCBHZEKOyQmZETb45Yycd857d+xO3bkbbVKdwXXNc7IVbvoH4uEiXA==</latexit>
<latexit sha1_base64="Nq3KI4ImyEkMTXPe9r+FU5fwrlw=">AAACeHicbZBNa9tAEIbH6lfqfsRteutlqenXxUgmtD2a9pJjAnESsIwZrcbO4tVK2R0luMI/sD+hv6LXJoeQlZxCHHdg4eV5Z3Z236TQynEY/m4FDx4+evxk62n72fMXL7c7r14fuby0koYy17k9SdCRVoaGrFjTSWEJs0TTcTL/UfvH52Sdys0hLwoaZzgzaqokskeTjowTmilTsZr/lOlStOOzC/FBxJKtrvrLWnoQx2s82uSMdlbdxWTSf5dOOt2wFzYlNkV0K7qDN9DU/qRzHae5LDMyLDU6N4rCgscVWlZS07Idl44KlHOc0chLgxm5cdWEsRTvPUnFNLf+GBYNvTtRYebcIkt8Z4Z86u57NfyfNyp5+m1cKVOUTEauFk1LLTgXdbIiVZYk64UXKK3ybxXyFC1K9vm319bUway+USutEot2UZ2VaBrHBxbdj2dTHPV70Zfe7sFud/B9lRxswVt4B58ggq8wgD3YhyFI+AV/4C9ctq4CEXwMPq9ag9btzA6sVdC/AZu8wT0=</latexit>
T
Toffoli
T T†
(CCX,
CCNOT) H T† T T† T H
Figure 4: The Toffoli gate diagram, showing that it performs a NOT operation on the target qubit if both
the control qubits are in the |1⟩ state. On the right is one of its standard decompositions into CNOT and
single-qubit gates, in which 5 CNOT gates are needed to implement one Toffoli gate.
The stock price follows a random walk in continuous time with a variance rate proportional to
the square of the stock price. (Scholes and Black, 1973, §2(b)).
A classical random walk with unit steps up-or-down leads to a binomial distribution, which at large
scales is approximated by a corresponding normal distribution. Thus, for large simulations, the simplifying
assumption of a fixed size for each step is immaterial, because the overall distribution is normal. However,
the most standard formulation for the Black-Scholes model assume that the price change for each unit of
time is not fixed, but (log-)normally distributed. The Black-Scholes formula has been widely used as a
pricing tool: indeed, over-reliance on the model, and the amounts of money entrusted to it, have been
found to be significant contributors to the 2008 market crash (Cady, 2015; Wilmott and Orrell, 2017). One
particular observation is that the assumption of a constant rate of volatility is not borne out by the long-tail
of variations in strike-price, leading to the claim that a ‘volatility smile’ distribution is a more faithful model
in practice (Orrell and Richards, 2023).
Quantum random walks were introduced in the 1990s (Aharonov et al., 1993) as a counterpart for clas-
sical random walks, and have become a rich and established area of quantum modeling (Venegas-Andraca,
2012). In economics, quantum walks have been proposed as an alternative that takes into account key fac-
tors including varying beliefs or opinions about the future, and the transactions between different traders
(Orrell, 2020, Ch 7). Another anticipated benefit of these quantum walk models is that they will work
natively on quantum computers, when large fault-tolerant quantum hardware is available (Orrell, 2021).
Quantum walks are thus expected to be a powerful component in pricing models: for example, they may
be used to model the input distributions on which the Monte Carlo methods proposed by Stamatopoulos
et al. (2020) depend.
Often the term ‘quantum walk’ is preferred to the term ‘quantum random walk’, not only for brevity, but
because the internal state of a quantum walk is typically an entirely deterministic superposition of different
states. For example, a walk that starts in position 0 with a 50-50 chance of going in either direction will,
after one step, be in a superposition of the states representing positions −1 and +1, with equal amplitudes
in the superposition. The quantum state vector representing this superposition can be predicted exactly: it
is only the measurement outcome that is probabilistic, when one of these distinct possibilities is randomly
selected.
5
In the most standard presentation, a quantum walk consists of a quantum walker and quantum coin. At
each turn, the coin is tossed, and the walker’s position moves depending on the coin’s resulting state. A
canonical example is an unrestricted discrete walk, where the positions correspond to integers, and each
move is a single step, represented by incrementing the position integer by ±1.
This leads to an elegant expression for the shift or translation operator (Venegas-Andraca, 2012, Eq. 9)
(Orrell, 2020, §7.1): X X
|0⟩c ⟨0| ⊗ |i + 1⟩p ⟨i| + |1⟩c ⟨1| ⊗ |i − 1⟩p ⟨i| . (1)
i i
The c and p subscripts refer to the coin and position registers. The positions are represented by integer
states |n⟩ for n ∈ Z.
The quantum walk model can be used in two different but related ways. Firstly, it can be used to model
evolving beliefs or opinions, for example, subjective estimates of what a particular asset will be worth at
a given future time. Different and even contradictory beliefs can be modeled in the quantum walk as a
coherent superposition of different possible values. Secondly, it can be used to model actual asset prices, as
observed in recorded transactions. A transaction behaves like a measurement on the quantum walk state,
which forces the coherent superposition into a particular state. Some ways of inserting and tuning the level
of decoherence are discussed by Orrell (2021), noting in particular that if the quantum walk is measured at
every time increment, decoherence is complete and the quantum walk collapses to the classical version.
An interesting feature of long, coherent quantum walks is that the distributions become two-tailed, as do
those of the quantum harmonic oscillator at higher energy levels (Orrell, 2020, Ch 10). Quantum dynamic
models have been used to represent several cognitive scenarios (Busemeyer and Bruza, 2012, Ch 9), and
oscillator models in particular have been used to model fluctuations in the stock market (Orrell, 2022). It is
intriguing that the distributions produced by quantum walks with and without decoherence are similar to
those produced by oscillator models at high and low energy levels.
This has provided considerable motivation for implementing quantum walk and oscillator models.
However, experiments in simulating quantum walks and harmonic oscillators on real quantum comput-
ers have been very small so far, restricted to just 2 or 3 qubits, and have reported very noisy results using
superconducting hardware (Qiang et al., 2016; Kadian et al., 2021; Puengtambol et al., 2021). The reasons
for this are explained in the next section.
6
a a
<latexit sha1_base64="zUGtFotOymYgIzpVGX6s1MzC2qs=">AAADh3icbVJNb9NAEF3HQNPw0RS4cVmRFnEKdqnaHgNcOBaJtJXiKFqvN+nK67WzO6Yylv8n/JmK9QdunHQkS0/vzb4ZP42fCK7Bcf5aPfvJ02d7/f3B8xcvXx0MD19f6ThVlE1pLGJ14xPNBJdsChwEu0kUI5Ev2LUffiv1619MaR7Ln5AlbB6RleRLTgkYanForT2frbjMgYe/aVDggSc0cBrmpMDe+g5/wB4FJfLPRQvdCtbaRkdLe+rBwWv9/K7fSdUNRK3yYovccGv1jrHfMaZEqQxzWXR32tnsYVbb0BjixiNOoeN85IUMcqc4Kpon2xaPWHYj0GlUOzIZ/I94MRw5Y6cqvAvcBowmb1FVl4vhvRfENI2YBCqI1jPXSWCeE2VGCFYMvFSzhNCQrNjMQEkipud5dRsFPjZMgJexMp8EXLGbL3ISaZ1FJtLjiMCt3tZK8jFtlsLyYp5zmaTAJK0HLVOBIcbloeGAK0ZBZAYQqrjZFdNboggFc46DzpgymPo3SiS4r4jK8nVKZKWYwNzteHbB1cnYPRuf/jgdTb7WyaE+eofeo4/IRedogr6jSzRF1Ppj3ff2en173/5kn9kXdWvPat68QZ2yv/wD13oclg==</latexit>
b b
|0i sum
Figure 5: A quantum full adder circuit, first introduced by Feynman (1986), uses 2 Toffoli and 3 CNOT
gates, which is at least 11 2-qubit gates.
Figure 6: Ideal simulated quantum walk after 14 and 15 steps (left), compared with results simulated with
expected noise (right). The ideal distribution has the two-tailed peaks characteristic of a quantum walk or
harmonic oscillator, but this quickly gets lost with noise (right).
in NISQ-era machines limits the number of steps we can reliably count (Orts et al., 2020).
Choosing a binary positional encoding, as used in classical computing, makes quantum counting very
susceptible to 2-qubit gate errors, because manipulating binary encodings takes a lot of entanglement and
coordination between qubits. To compute the sum of two binary numbers A and B of bitlength n using
classical Boolean algebra, we add (XOR) the least significant bits, and then at each other position we add
the corresponding bits along with a ‘carry’ from the previous stage, setting the output and passing a ‘carry’
on to the next position. Feynman (1986) explained the quantum gate operations needed for each such
step in the quantum full adder circuit of Figure 5, and modern versions are optimized variants on this
theme (Orts et al., 2020). If it takes 11 2-qubit gates for each full-adder, then adding two single-byte (8-bit)
numbers using this approach uses ∼ 80 2-qubit gates, so by the time such a register has successively added
10 numbers, the chances of an error are over 50% even with a two-qubit gate fidelity of 99.9%, which is on
the high-end of performance estimates at the time of writing (IonQ Aria, 2022).
Error rates with quantum counting can thus undermine the simulation of quantum walks, and block this
application of quantum finance. The problem is demonstrated in Figure 6, which compares quantum walks
with ideal outcomes and with noise. This shows the vulnerability of the quantum counting process to noise
after a handful of steps. It should be noted that some noise or decoherence can be useful in quantum models,
for example, in deriving price estimates from many quantum walk simulations (Orrell, 2020, Ch 6). The
predictions later in this paper are also averages over many somewhat-noisy circuits: the key engineering
challenge here is to understand the noisy behaviors, and find those that work well enough to enable the
7
task at hand.
There are many optimizations and alternatives. We expect progress in quantum hardware to enable
greater fidelity and stability, and eventually mid-circuit quantum error-correction should make the current
problems with quantum counting obsolete — but this comes at the cost of waiting for fault-tolerant quan-
tum computing. Quantum addition algorithms can be optimized (Cuccaro et al., 2004; Gidney, 2018), and
the incremental operation of counting can be made simpler than repeated full-register addition (Li et al.,
2014). An interesting benchmark challenge could be to design and evaluate quantum circuits and see how
far they can count with > 50% fidelity, but that task is not undertaken here, because none of these meth-
ods simplify counting enough for many successive counting operations on nontrivial quantum registers to
be performed accurately yet. Instead, this paper proposes alternative circuits that can be used to simulate
steps and positions in a walk, without requiring exact counting.
8
Figure 7: Arc counter circuit
1 160 20
500 2 21
3 140 22
400 4 120 23
5 24
100 25
300
80
200 60
40
100
20
0 0
0 20 40 60 80 100 120 0 50 100 150 200
Figure 8: Quantum walk results after different numbers of steps with arc counter circuit and an 8-qubit
register. 1000 shots for each number of steps.
that qubit’s significance, as in Figure 7. This means that the n + 1th qubit rotates at half the rate of the nth
qubit.
A good analogy for this representation is to think of each qubit as one of the hands on a traditional
analog clock. On an analog clock, the minute hand cycles at 12 times the speed of the hour hand, and the
second hand 60 times faster still, whereas in our binary clock, the ratio between the speed of rotation of
each successive pair of ‘hands’ is 2:1. This analogy works well for the standard binary positional notation
for integers, which can be thought of as a binary digital clock. In a digital register (or an abacus), the digits
logically depend on one another for correct incrementing, because we need to know that one register is full
before we increment the next. By contrast, the hour hand on a clock does not ‘carry’ information when
the minute hand completes a cycle — it just rotates at its own slower pace. Thus the rotation encoding
clock-based design requires much less coordination (and hence entanglement) between the qubits.
This comes at a representational cost — the register does not represent exact integers, and random
variations in the outputs are expected, because many fractional angles are used throughout the circuit.
(This is true for the basic counting operation, irrespective of whether the counting is coupled with a ‘coin
toss’ operation.)
Sample results from a quantum walk with an 8 qubit register are shown in Figure 8. Statistically note-
worthy properties include:
9
• The mean distance from the starting point generally increases with the number of steps.
• In some cases, the position appears to jump ahead, because a high-order qubit is measured in the |1⟩
state. This can happen (with low probability) after just a single step.
• There are sometimes peaks in the distributions after specific powers of two or their combinations (e.g.,
peaks at 48, 64, 96). It may be possible and desirable to find ways to smooth out these peaks.
Since there are no 2-qubit entangling gates, error rates are lower, but there’s also no physical quantum
advantage from this design — it is easy to model this distribution on a classical computer. It’s possible
that such distributions might be useful models for random processes, but this would not require quantum
computers to simulate.
10
Figure 9: Two-directional arc walk circuit from combining positive and negative distributions. 1000 shots
for each number of steps.
11
250
10 10
11 80 11
200 12 12
13 13
14 60 14
150 15 15
40
100
50 20
0 0
0 50 100 150 150 100 50 0 50 100 150 200
Figure 11: Example results for arc walk counter circuit (8 counter qubits, 10 to 15 steps). On the left is a
purely incrementing circuit, on the right is a two-way walk using the reversal and superposition technique
above. 1000 shots for each number of steps.
12
12 3 9 3
4 8 4
10 5 5
6 7 6
7 7
8 8 6 8
9 9
10 5 10
6
4
4 3
2 2
1
0 20 40 60 80 100 50 0 50 100
Figure 13: Random jump circuit results (positive only on the left, two-way walk using reversal and super-
position on the right). 30 random jump circuits per step, 30 shots per circuit.
<latexit sha1_base64="T1zFrnKHuRl/OrOIX3MX7xDYjkY=">AAAC+XicdVJLb9NAEF6bVwmPhseNy4oIxCmyqwo4RnBBnIpE2khxFI3Xk3SV9drdHVO5K/8Ybogrv4Y7vwPhR0CpS0fa1adv5pvZmdk4V9JSEPz0/Bs3b92+s3d3cO/+g4f7w0ePj21WGIFTkanMzGKwqKTGKUlSOMsNQhorPIk37xv/yRc0Vmb6M5U5LlJYa7mSAqimlsNfUYxrqR3JzYVIKj6Izs75Sx6tgdDNqgYKMsodVD22C/t3R9F1yrCvtOeQb9nr5NEGyQVVG0Ng1q6foyFnO+yuYJsOdfK3q+VwFIyD1vhVEG7BaPKUtXa0HP6OkkwUKWoSCqydh0FOCweGpFBYDaLCYg5iA2uc11BDinbh2m1U/EXNJHyVmfpo4i27q3CQWlumcR2ZAp3avq8h/+ebF7R6u3BS5wWhFl2hVaE4ZbxZLU+kQUGqrAEII+u3cnEKBgTVH2BwqUwzmK6NBikZGzClOytAt556YGF/PFfB8cE4fD0+/HQ4mrzrJsf22DP2nL1iIXvDJuwDO2JTJryPXu6V3oXv/K/+N/97F+p7W80Tdsn8H38AGX/sdQ==</latexit>
X X
X X
|0i X |0i
Figure 14: Boolean in-place disjunction circuit, which sets the state of the higher order qubit to the Boolean
OR of the input states of the lower and higher order qubits. The X gates surrounding the Toffoli gate
implement the usual NOT operations to turn an AND conjunction operator into an OR disjunction operator
(A∨B = ¬(¬A∧¬B)). Finally, the swap gate and the reset to |0⟩ operations ensure that the output is written
back into the higher order qubit, and the ancilla qubit is reset to its initial state. Using this construction
several times requires mid-circuit reset, so that the ancilla qubit can be reused.
Results for up to 9 steps, using an 8 qubit counter register, are shown in Figure 13. As expected, each
walk can go both ways (because randomly flipping a bit can reduce as well as increase a number), though
the average tendency for an individual walk is to increase. This is because the registers are initialized to
zero, so the process randomly diffuses out from zero.
13
3 16 3
10 4 4
5 14 5
6 12 6
8 7 7
8 10 8
9 9
6 10 8 10
4 6
4
2 2
0 10 20 30 40 50 0 50 100
Figure 15: Cascading disjunction circuit results (positive only on the left, two-way walk using reversal and
superposition on the right). 30 random jump circuits with cascading disjunctions per step, 30 shots per
circuit.
reversal and superposition on the right). The walks are still random, but propagating more reliably with
the cascading disjunction components.
Another variant of this technique would be to add a conjunction between the coin qubit and the target
qubit, that sets the next higher qubit to |1⟩ before reversing the target qubit to |0⟩. This would behave
like a limited-carry operation: it performs some of the coordination between qubit values in the traditional
bit-counter circuits in the literature, but much less, in a much more targeted fashion.
14
Random
Binary Arc Random jump with
Steps Arc walk cascading
counter counter jump disjunc-
tions
0 0.000 0.000 0.000 0.000 0.000
1 1.000 1.108 0.491 4.386 1.133
2 2.000 3.276 1.522 3.383 5.281
3 3.000 4.742 2.079 8.488 5.783
4 4.000 6.764 2.826 11.450 6.341
5 5.000 8.233 3.833 9.317 13.539
6 6.000 10.564 4.860 12.468 15.510
7 7.000 10.955 6.476 13.964 16.366
8 8.000 12.664 6.458 9.804 21.479
9 9.000 15.097 6.958 14.261 16.837
10 10.000 17.836 7.810 12.261 21.919
Binary Arc
Binary Arc Random
Steps counter counter Arc walk
counter counter jump
(QPU) (QPU)
0 0.000 0.00 0.000 0.002 0.000 0.000
1 15.606 21.946 1.312 1.94 2.562 1.742
2 25.123 32.771 3.461 4.412 3.859 2.618
3 30.504 32.802 5.488 5.782 6.223 5.919
4 29.921 33.603 6.824 6.236 8.721 9.279
5 31.046 32.212 8.931 10.208 10.773 10.250
6 31.732 33.317 10.442 12.984 12.627 10.981
7 30.837 32.978 11.183 12.389 14.978 12.401
8 31.894 33.775 13.148 13.682 15.898 14.568
9 30.99 32.989 14.614 15.504 17.300 12.684
10 31.912 34.101 18.459 19.666 20.907 17.216
Table 2: Average distances traveled after n steps, noisy simulation and QPU
deterministically, so these results include 1000 shots for a single circuit, and all the randomness is quantum.
Key findings include:
• The binary counter results are perfect with ideal simulation, but are rendered useless in the noisy
simulation. They quickly converge to a random number around 32, which is the average
• For all the other circuits, the difference between ideal and noisy results is much less.
• The average results for the arc counter and arc walk circuits are the most reliable for simulating a
monotonically-increasing position, with or without noise.
• The random jump results also tend to increase, but tend to plateau and then move up and down
randomly. (This randomness is smaller with a larger register.)
The QPU results for binary and arc counting are compared graphically in Figure 16. This shows how
quickly the binary counter becomes useless on a real QPU, whereas by contrast, the arc counter QPU results
stay close to the ideal simulated results.
15
Figure 16: Ideal and actual QPU results for binary counter and arc walk circuits. The QPU are much closer
to the ideal monotonically-increasing results for the arc counter, whereas they are useless after 2 steps for
the binary counter
16
Dow Jones Index Arc Walk Circuit Random Jump Circuit
5 yrs of data from 8 qubits 8 qubits
Yahoo! Finance 1000 shots per job 8 circuits per sample
1000 shots per job
Figure 17: Distributions of relative changes in the Dow Jones Industrial Average, and in quantum approx-
imate counting simulations. The dashed curve (red) shows the best-fit normal distribution, that is, the
normal with the same mean and variance.
In related work, IonQ quantum computers have also been used to model the normal distribution itself,
using a matrix product state technique that can readily be adapted to other distributions, because it relies on
piecewise polynomial approximation (Iaconis et al., 2023). One of the longer-term promises of such work
is that such distributions can be used as inputs for models such as the Monte Carlo simulations advocated
by Egger et al. (2020). If we have a reliable circuit for preparing a particular distribution, then such a circuit
could be used as input for Monte Carlo modeling by entangling its output with the simulated variables,
rather than by sampling an individual number from the distribution and using this as a single ‘classical’
random input value.
The probability of transitioning to the state |1⟩ is thus proportional to sin2 (θ) which tends to zero for small
θ, and it is easy to see that if a larger angle is divided into smaller and smaller increments, the prob-
ability of observing a transition to the state |1⟩ in any of these increments also tends to zero, because
limn→∞ (n) sin2 n → 0.
This phenomenon is sometimes called the quantum Zeno effect, after Zeno’s classical paradox of motion.
Of crucial interest for this paper, such effects have also been observed in psychology. Kvam et al. (2015)
demonstrated that participants are likely to form less extreme judgments of moving scenes if asked to
judge the motion in smaller time-frames, and Yearsley and Pothos (2016) demonstrated that participants
evaluating evidence in a criminal trial are more likely to change their minds if several pieces of evidence
are presented before asking for a decision.
17
Figure 18: Arc counter circuit results, simulating the results of quantum walks with 20 steps, with and
without mid-circuit measurement at the given positions.
It is easy to add mid-circuit measurement to our quantum approximate counting circuits and to evaluate
the results, at least in simulation. (The availability of mid-circuit measurement varies across quantum plat-
forms currently, partly because the accuracy of the measurement and reset operations is hard to guarantee.)
Example results are shown in Figure 18, simulating walks with 20 steps, with no mid-measurement, mea-
surement every 7 steps, and measurement every step. The average positions reached by these walks were
35.6, 14.3, and 5.7 respectively, so as expected, the use of mid-measurement reduces the average distance
traveled in the quantum walk. (It is not always this simple, particularly due to periodicity.) A further step
for this research would be to experiment with parameters including the number of steps, distribution of
step sizes, and frequency of measurements, to see if different levels of deliberate decoherence can bring the
quantum results of Figure 18 closer to the distributions observed with the Dow Jones Index and other asset
prices.
18
published by Zillow.
When a house sells for less than its original listing price, we assume that this indicates a difference
between the seller’s and the buyer’s estimate of the house’s value. Larger uncertainties in the market
would support larger differences of opinion. Even when considering monthly averages of data, we would
expect that a smaller number of sales in a given area would contribute to greater market uncertainty, and
this should correlate with a greater difference between the list price and the sale price.
By contrast, when a house sells for more than its original listing price, we assume that there are other fac-
tors involved: in particular, this situation is common when there are other bids on the property from other
potential buyers, so a minimum value is already established without the need for previous comparable
transactions.
Thus, we assume that the markets where lack of comparables is a primary factor in price uncertainty
are those where the average sale price is less than the average list price. It follows that, if we restrict our
attention to markets where the average sale price is less than the average list price, we should see evidence
that lower transaction volumes are correlated with greater disparities between list price and sale price.
Data used to test these hypotheses was gathered from the Zillow Housing Data portal1 . The datasets
are summary statistics: counts and averages. These are only comparable within a given metro area: for
example, 2000 sales in a month would be very low for New York, NY, and very high for Wichita, KS. Thus
we compute correlations by comparing monthly statistics within each metro area.
The algorithmic steps are as follows:
– Collect the sales count and the average list-to-sale price ratio.
– If the average list-to-sale price ratio is greater than 1, skip this month.
• This gives a set of (count, ratio) pairs, e.g. [(822, 0.98), (785, 0.96), (803, 0.97)], etc.
• Compute the Pearson correlation coefficient between these sales counts and list-to-sale price ratios.
• Gather the Pearson correlation coefficients into a histogram to see if there is a general trend.
The result is in Figure 19. Nearly all the correlations are strongly positive. This shows that, in cases
where a house is sold for less than its asking price, there is a very strong correlation between the translation
volume, and the closeness of the list and the sale prices.
In addition to the number of sales in a region, it is instructive to look at the correlation between home
prices and the list-to-sale price ratio since the price of a home is the driving factor in a home purchase. As
a proxy for home prices we use the Zillow Home Value Index (ZHVI) which reflects the typical value for
homes in the 35th to 65th percentile range.
Figure 20 shows both the raw monthly values (above) for the San Francisco region as well as monthly
changes (below) comparing the ZHVI and list-to-sale ratio values for the last 6 years. The bottom chart
shows more clearly the positive correlation between the two if we consider monthly changes. Monthly
changes give a better sense of market ups and downs. As prices fall, we see that buyers are more likely to
bid lower than the list price. Conversely, in a rising market, buyers are more likely to bid more than the list
price.
1 https://www.zillow.com/research/data/, accessed 2023-10-05.
19
Figure 19: Histogram showing correlations between larger numbers of transactions and smaller list-to-sale
price differences. Data from Zillow Housing Data.
Figure 20: Monthly values of Zillow Home Value Index (ZHVI) and list-to-sale ratio as well as monthly
changes that show positive correlation in both rising and falling markets.
20
That fewer transactions correlates with larger list-to-sale price ratios is in line with the trends expected
from quantum economics models, in which various beliefs and opinions about value can evolve and di-
verge more when there are fewer transactions or measurements. The correspondence between overall price
changes and list-to-sale ratios could be described by a quantum walk with a clear directional momentum.
However, it is also easy to propose simple non-quantum models for these behaviors. Fewer comparable
samples leads to greater sampling error and thus greater price uncertainty: thus lower liquidity brings
higher volatility. One potential strategy for evaluating and distinguishing which approaches are better
would be to consider the dynamics / evolution of prices in such models: for example, to try applying
the quantum Monte Carlo sampling reduction described by Egger et al. (2020) to the problem of making
accurate price estimates with fewer comparable sales.
21
8 Conclusions and Future Work
This paper has introduced and explored quantum approximate counting circuits, as fault-tolerant alterna-
tives to the traditional quantum walk design, particularly for the way position is tracked and incremented.
The new designs presented here lack some of the mathematical elegance, and the theoretical results, that
accompany the traditional quantum walk design: and in particular, there are no longer unit increment
and decrement operators that correspond to the ladder operators of a quantum oscillator. However, the
enormous advantage for the simpler models presented here is that they behave much more accurately on
NISQ-era quantum hardware, which could contribute to commercially advantageous applications of quan-
tum computers in economics.
These are just prototype designs so far. The main next steps for this work are to evaluate the proposals
more quantitatively, answering the following two questions:
1. How do results on NISQ-era quantum computers correspond to ideal or simulated results for small
circuits, and what does this indicate about the expected behavior on quantum hardware for systems
that are too big to simulate on classical hardware?
2. How do results compare with the distributions observed with real market behaviors?
The ideal outcome of this research is that we would find circuit walk designs that are robust enough to
given better models of market behavior that include some of the benefits of quantum approaches noted by
Orrell (2020), while being able to run on today’s quantum hardware without waiting for error-correction.
Given the crucial and explicit role that measurement plays in quantum models, it is possible that some
of the earliest such quantum advantages will be apparent in markets where a small number of significant
transactions can dramatically influence the price of a particular asset. An initial analysis suggests that the
housing market may be an appropriate area to test this hypothesis.
This work can be seen as part of a larger program to bring value in economic modeling on quantum com-
puters. Other successes for quantum circuit designs include modeling and sampling from key distributions
(Iaconis et al., 2023), and demonstrating particularly effective time-series models using copula functions
implemented using entanglement (Zhu et al., 2022a). Related work in cognitive science has demonstrated
that simple quantum circuits can also be used to model decision-making processes (Widdows and Rani,
2022). In the next few years, it is likely that several such small components, being developed today, will be
used as key building blocks in the first profitable applications of quantum computing in economics.
9 Acknowledgements
The author would like to thank Emmanuel Pothos and David Orrell for interesting conversations and en-
couragement.
10 Funding
This work was funded by IonQ, Inc.
References
Aharonov Y, Davidovich L and Zagury N (1993) Quantum random walks. Physical Review A 48(2): 1687.
22
Bachelier L (1900) Théorie de la spéculation. In: Annales scientifiques de l’École normale supérieure, volume 17. pp. 21–86.
Busemeyer JR and Bruza PD (2012) Quantum models of cognition and decision. : Cambridge University Press.
Chen W, Gan J, Zhang JN, Matuskevich D and Kim K (2021) Quantum computation and simulation with vibrational
modes of trapped ions. Chinese Physics B 30(6): 060311.
Cuccaro SA, Draper TG, Kutin SA and Moulton DP (2004) A new quantum ripple-carry addition circuit. arXiv preprint
quant-ph/0410184 .
Egger DJ, Gambella C, Marecek J, McFaddin S, Mevissen M, Raymond R, Simonetto A, Woerner S and Yndurain E (2020)
Quantum computing for finance: State-of-the-art and future prospects. IEEE Transactions on Quantum Engineering 1:
1–24.
Haven E, Khrennikov AY and Robinson TR (2017) Quantum methods in social science: A first course. World Scientific
Publishing Company.
Iaconis J, Johri S and Zhu EY (2023) Quantum state preparation of normal distributions using matrix product states.
arXiv preprint arXiv:2303.01562 .
IonQ Aria (2022) IonQ Aria furthers lead as world’s most powerful quantum computer. URL https://ionq.com/
news/february-23-2022-ionq-aria-furthers-lead. Accessed 2022-05-28.
Jain VK, Behera BK and Panigrahi PK (2021) Quantum simulation of discretized harmonic oscillator. Quantum Studies:
Mathematics and Foundations 8: 375–390.
Kadian K, Garhwal S and Kumar A (2021) Quantum walk and its application domains: A systematic review. Computer
Science Review 41: 100419.
Kitto K and Widdows D (2016) Ideologies and their points of view. In: Quantum Interaction: 9th International Conference,
QI 2015, Filzbach, Switzerland, July 15-17, 2015, Revised Selected Papers 9. Springer, pp. 216–227.
Kvam PD, Pleskac TJ, Yu S and Busemeyer JR (2015) Interference effects of choice on confidence: Quantum characteris-
tics of evidence accumulation. Proceedings of the National Academy of Sciences 112(34): 10645–10650.
Li X, Yang G, Torres Jr CM, Zheng D and Wang KL (2014) A class of efficient quantum incrementer gates for quantum
circuit synthesis. International Journal of Modern Physics B 28(01): 1350191.
Mandelbrot B (1963) New methods in statistical economics. Journal of political economy 71(5): 421–440.
Nielsen MA and Chuang I (2002) Quantum computation and quantum information. Cambridge University Press Edition,
2016: American Association of Physics Teachers.
Orrell D (2020) Quantum Economics and Finance: An Applied Mathematics Introduction. New York: Panda Ohana Publish-
ing. ISBN 9781916081611.
Orrell D (2021) A quantum walk model of financial options. Wilmott 2021(112): 62–69.
Orrell D (2022) A quantum oscillator model of stock markets. Available at SSRN 3941518 .
Orrell D and Richards L (2023) Keep on smiling: Market imbalance, option pricing, and the volatility smile. Wilmott
124: 58–64. URL https://ssrn.com/abstract=4205729.
23
Orts F, Ortega G, Combarro EF and Garzón EM (2020) A review on reversible quantum adders. Journal of Network and
Computer Applications 170: 102810.
Pagourtzi E, Assimakopoulos V, Hatzichristos T and French N (2003) Real estate appraisal: a review of valuation meth-
ods. Journal of property investment & finance 21(4): 383–401.
Puengtambol W, Prechaprapranwong P and Taetragool U (2021) Implementation of quantum random walk on a real
quantum computer. In: Journal of Physics: Conference Series, volume 1719. IOP Publishing, p. 012103.
Qiang X, Loke T, Montanaro A, Aungskunsiri K, Zhou X, O’Brien JL, Wang JB and Matthews JC (2016) Efficient quantum
walk on a quantum processor. Nature communications 7(1): 11511.
Rebentrost P, Gupt B and Bromley TR (2018) Quantum computational finance: Monte carlo pricing of financial deriva-
tives. Physical Review A 98(2): 022321.
Scholes M and Black F (1973) The pricing of options and corporate liabilities. Journal of political Economy 81(3): 637–654.
Schuld M and Petruccione F (2021) Machine Learning with Quantum Computers. Springer.
Stamatopoulos N, Egger DJ, Sun Y, Zoufal C, Iten R, Shen N and Woerner S (2020) Option pricing using quantum
computers. Quantum 4: 291.
Standage T (1998) The Victorian Internet: The remarkable story of the telegraph and the nineteenth century’s online pioneers.
Phoenix London.
Venegas-Andraca SE (2012) Quantum walks: a comprehensive review. Quantum Information Processing 11(5): 1015–1106.
Waddup OJ, Yearsley JM, Blasiak P and Pothos EM (2023) Temporal bell inequalities in cognition. Psychonomic Bulletin
& Review : 1–8.
Widdows D and Rani J (2022) Quantum circuits modeling cognitive order and disjunction effects. Internal Preprint at
https://drive.google.com/file/d/1DJQdcvXuty2KdBdQHBXm0W1K6COQ9T9P.
Wilmott P and Orrell D (2017) The money formula: Dodgy finance, pseudo science, and how mathematicians took over the
markets. John Wiley & Sons.
Wright K, Beck KM, Debnath S, Amini J, Nam Y, Grzesiak N, Chen JS, Pisenti N, Chmielewski M, Collins C et al. (2019)
Benchmarking an 11-qubit quantum computer. Nature communications 10(1): 1–6.
Xia F, Liu J, Nie H, Fu Y, Wan L and Kong X (2019) Random walks: A review of algorithms and applications. IEEE
Transactions on Emerging Topics in Computational Intelligence 4(2): 95–107.
Yearsley JM and Pothos EM (2016) Zeno’s paradox in decision-making. Proceedings of the Royal Society B: Biological
Sciences 283(1828): 20160291.
Zhu D, Shen W, Giani A, Majumder SR, Neculaes B and Johri S (2022a) Copula-based risk aggregation with trapped
ion quantum computers. arXiv preprint arXiv:2206.11937 .
Zhu EY, Johri S, Bacon D, Esencan M, Kim J, Muir M, Murgai N, Nguyen J, Pisenti N, Schouela A et al. (2022b) Genera-
tive quantum learning of joint probability distribution functions. Physical Review Research 4(4): 043092.
Zi-Yi G (2017) Heavy-tailed distributions and risk management of equity market tail events. Journal of Risk and Control
4(1).
24