Skip to content

Quantum Trick or Treat

November 1, 2021


Are crazy quantum walks fact or fiction?

NOAA src

Lesley Smith is a climate researcher at the University of Colorado’s Cooperative Institute for Research in Environmental Sciences. She is also associated to the National Oceanic and Atmospheric Administration’s Physical Science Laboratory and has worked on behalf of several other national organizations. Her PhD was in particle physics at the University of Kansas, and she incorporates quantum physics into some of her non-nonfiction writings.

Today we discuss a quantum walk on a simple graph with behavior wild enough to inspire at least the latter kind of writing.

Smith was most recently second author on a paper in the April, 2021 Journal of Climate whose title immediately speaks a vital subject: “Explaining the Spatial Pattern of U.S. Extreme Daily Precipitation Change.” The years-long drought in the US West has critically depleted lakes and reservoirs and rivers. But just one week ago, much of California had its wettest day ever. Both the drought and the pattern of extremes have spread across the country, as many easily-found stories this year attest. The question is how to distinguish standard fluctuations from deep-seated causal factors.

As for her other writings, this is representative:

Amazon page

The first lines are likewise topical: “Despite it being October, the sun felt warm on my skin. Forest-fire smoke from states west of Colorado gave the air a three-dimensional quality…” She has also written novels titled Quantum Murder, The Quantum Cop, and Quantum Juneteenth, plus just now a short story, “Lucky Halloween,” about a quantum computer scientist. She is writing a series called Kat Cubed—well we will raise Schrödinger cat-walks to powers well beyond cubed below.

Climate and Dynamism

The climate paper’s opening paragraphs review the the general prediction of greater precipitation from the thermal component of climate change and the effects of changes in atmospheric currents. Then it observes:

Yet, the regional pattern of observed heavy precipitation trends across the United States departs from this theoretical expectation, being remarkable for its heterogeneity rather than the more uniform structure that thermodynamic effects alone would predict. … The regionally diverse trends in extreme precipitation are likely related to dynamical factors.

The paper speaks to care needed to distinguish causes and draws two conclusions:

First, the absence of appreciable increases [in annual maximum single day of rain] over the West is attributed to a dynamical effect of the observed centennial-scale trends in [sea-surface temperature], sea ice, and atmospheric composition. … A second factor is sampling variability […from…] a combination of moderate forced increases comingled with a strong articulation of internal atmospheric variability.

The nub of the latter from my perspective is distinguishing what we can infer on the basis of increased dynamism alone. In our post seven years ago on global warming, I opined that greater energy and dynamism should be treated as more-primary effects than temperature. The challenge of inferring further effects from dynamism is recognized by another excerpt from the conclusions:

Recognizing the potency of unforced atmospheric dynamics, [one cannot discount] a scenario in which internal variability could mute if not reverse the observed upward trend in eastern U.S. [extreme rains].

In the rest of this post, we’ll offer an example of wild dynamics from a small quantum system that may speak to both her métiers.

Classical Random and Quantum Walks

Consider the following graph {G}. It looks like a Q-tip. The two self-loops depart from usual examples of undirected graphs.



A classical random walk on {G} starts at some node—say node {1}—and at each step flips a coin to determine the next node. The “walker” stays on node {1} if the result is heads ({H}) but goes to node {2} if tails ({T}). The walker cannot go from {1} to {3} in one step because there is no edge. The probabilistic transition matrix {A_G} is given at right. The vector

\displaystyle  \begin{bmatrix} 0.5 \\ 0.25 \\ 0.25 \end{bmatrix} = A_G\cdot \begin{bmatrix} 0.5 \\ 0.5 \\ 0 \end{bmatrix} = A_G^2 \cdot x^{(0)} \qquad\text{where}\qquad x^{(0)} = \begin{bmatrix} 1 \\ 0 \\ 0 \end{bmatrix}

gives the probabilities for each node after two steps. (Although multiplying row vectors on the left is commonly used for classical walks, we will use column vectors for consistency with quantum usage; that {A_G} is symmetric moots the difference anyway.) For any {t}, {x^{(t)} = A_G^t x_0} can be called the “classical state” of the walk after {t} steps.

The state {x} whose entries {x[i]} are given by the degree of node {i} divided by the sum of the degrees makes {A_G x = x}, and thus gives a stable distribution for the walk. A central theorem states that for every {G} that is not bipartite, the powers of {A_G} converge entrywise in every column to a unique stable distribution. When {G} is regular, as here, this is simply the uniform distribution. The convergence in our case is quite quick, as signified by the exact values

\displaystyle  A_G^{10} = \begin{bmatrix} 0.333984375 & 0.3330078125 & 0.3330078125 \\ 0.3330078125 & 0.333984375 & 0.3330078125 \\ 0.3330078125 & 0.3330078125 & 0.333984375 \end{bmatrix}.

The quantum case needs to bind the graph and the “coin” into one system, so its domain has six members, which we list in order

\displaystyle  \mathcal{S} = \{(1,H),(1,T),(2,H),(2,T),(3,H),(3,T)\}.

In quantum notation, this is the tensor product of the “node space” {V = \{1,2,3\}} and the “coin space” {\mathcal{C} = \{H,T\}}. The meaning of {(v,\mathfrak{c})} is that the walker just arrived at node {v} via the coin outcome {\mathfrak{c}}. It is incumbent that every node {v} can be reached on either outcome. Our {H}{T} labeling above complies. Then we get a permutation {P} on {S} from the action

\displaystyle  P(v,\mathfrak{c}) = (v',\mathfrak{c}),

where {v'} is the vertex reached from {v} by the outcome {\mathfrak{c}} from the classical walk, and {\mathfrak{c}} is kept the same. The penultimate trick is to expand this to a mapping {P': \mathcal{S} \times \mathcal{C} \rightarrow \mathcal{S}} by

\displaystyle  P'((v,\mathfrak{b}),\mathfrak{c}) = (v',\mathfrak{c}),

whereby the previous coin outcome {\mathfrak{b}} is “forgotten” while the current flip is preserved in the state. We can represent this by the directed graph {G'} shown next.



Incidentally, {G'} is planar (switch nodes {2H} and {2T} to see). The last trick is to label its edges by quantum amplitudes for the coin outcomes.

Realm of the Coin

Our quantum coin can be “flipped” by applying any chosen {2 \times 2} matrix {\mathbf{U}} to the current state of the coin, provided {\mathbf{U}} is a unitary matrix. In general, if {G} is a {d}-regular graph and we have a suitable labeling for results of a {d}-sided die, then {\mathbf{U}} can be a {d \times d} matrix acting on a single qudit, or we can use {\lceil \log_2 d \rceil} qubits to encode the coin. But here we just need one qubit. Before fixing the choice of coin matrix, we can represent it abstractly as

\displaystyle  \mathbf{U} = \begin{bmatrix} a & b \\ c & d \end{bmatrix}

To apply the coin on {\mathcal{S}}, we form the {6 \times 6} matrix {\mathbf{U'} = \mathbf{I} \otimes \mathbf{U}}, where {\mathbf{I}} is the {3 \times 3} identity matrix. To move the walker after the coin result, we apply the {6 \times 6} matrix {\mathbf{P}} of the permutation {P}. Since we use column vectors for states, the walk matrix is given by

\displaystyle  \mathbf{W} = \mathbf{P}\cdot\mathbf{U'} = \begin{bmatrix} 1 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 \\ \end{bmatrix} \cdot \begin{bmatrix} a & b & 0 & 0 & 0 & 0 \\ c & d & 0 & 0 & 0 & 0 \\ 0 & 0 & a & b & 0 & 0 \\ 0 & 0 & c & d & 0 & 0 \\ 0 & 0 & 0 & 0 & a & b \\ 0 & 0 & 0 & 0 & c & d \end{bmatrix} = \left[\begin{array}{cc|cc|cc} a & b & 0 & 0 & 0 & 0 \\ 0 & 0 & c & d & 0 & 0 \\ \hline 0 & 0 & 0 & 0 & a & b \\ c & d & 0 & 0 & 0 & 0 \\ \hline 0 & 0 & a & b & 0 & 0 \\ 0 & 0 & 0 & 0 & c & d \end{array}\right]

We have drawn lines to show how the coin acts between each pair of nodes. Notice that {\mathbf{W}} is no longer symmetric, so specifying column-vector inputs matters. For instance, the {b} in the top row comes in from the column at {1T} and exits at {1H}. This means the previous coin result was tails and took the walker (from node 2) to node {1}; now the coin gives heads so the walker stays at node {1}. We can diagram this action also on the expanded graph {G'}, so that the {b} goes on the arrow from {1T} to {1H}, and so on:

The first step has no “previous coin result” but we still have to specify one to initialize the coin state as well as the walker’s state. For starting on node {1} there is still an infinitude of combinations of the basis states {1H} and {1T}. If we fix the start as {1H}, this is like saying the coin was heads-up before the initial act of flipping it.

Walking Into Chaos

A natural and popular choice of coin is the Hadamard matrix

\displaystyle  \mathbf{H} = \frac{1}{\sqrt{2}}\begin{bmatrix} 1 & 1 \\ 1 & -1 \end{bmatrix}.

We need only mark the edges with {d = -1}, keeping the normalizing {\frac{1}{\sqrt{2}}} in mind. Here are {G'} and the three steps of the quantum walk after the first:

To visualize the {t}-step walk, say starting from {\mathfrak{s}_0 = 1H}, we take each possible path of length {t} and multiply the numbers on its edges. Over all paths that end at the same {\mathfrak{s}=(v,\mathfrak{c})} we add those products to give the amplitude {a_{u,v}}, which in our column-first notation equals the entry {\mathbf{W}^t[\mathfrak{s},\mathfrak{s}_0]}.

For example, there are two paths of {t=3} steps from {1H} back to itself. One takes the loop at {1H} three times and contributes the phase value {1}, but the other does {1H\rightarrow 2T \rightarrow 1T \rightarrow 1H} and picks up the value {-1}. Adding those values cancels them, leaving the red {0} shown at upper left in {\mathbf{W}^3}. If our walker is a Schrödinger cat, that part of it is not just “dead” but completely annihilated.

It is still possible for the cat to end on node {1} after being cubed, because the path {1H \rightarrow 1H \rightarrow 2T \rightarrow 1T} is not canceled. That the cat ends with phase {-1} does not matter, because to get the probabilities—given we started at {1H}—we take the squared magnitude of each entry in column {1} and add the adjacent pairs denoting the same vertex of the original graph {G}. Thus the cat has probability only

\displaystyle  p_1^{(3)} = \frac{0^2 + (-1)^2}{(2\sqrt{2})^2} = \frac{1}{8} = 0.125

of being still on node {1} after three time steps. Given the ease of looping at node {1}, this is counterintuitive. After nine and ten steps, we have:

\displaystyle  \mathbf{W}^9 = \frac{1}{2^{4.5}}\begin{bmatrix} 18 & 4 & -1 & 13 & -1 & -1\\ 13 & 1 & 4 & -18 & -1 & 1\\ -1 & 13 & -1 & -1 & 18 & 4\\ 4 & -18 & -1 & 1 & 13 & 1\\ -1 & -1 & 18 & 4 & -1 & 13\\ -1 & 1 & 13 & 1 & 4 & -18 \end{bmatrix} .\qquad \mathbf{W}^{10} = \frac{1}{32}\begin{bmatrix} 31 & 5 & 3 & -5 & -2 & 0\\ -5 & 31 & 0 & -2 & 5 & 3\\ -2 & 0 & 31 & 5 & 3 & -5\\ 5 & 3 & -5 & 31 & 0 & -2\\ 3 & -5 & -2 & 0 & 31 & 5\\ 0 & -2 & 5 & 3 & -5 & 31 \end{bmatrix}.

The corresponding node probabilities are {p^{(9)} = [ 0.96289, 0.0332, 0.00391]} and {p^{(10)} = [ 0.96289, 0.02832, 0.00879]}. The probability of the first node is the same between every odd and even step, echoing how starting at {1H} presumes a previous coin flip of heads, which kept the cat on node {1} at time {t = -1}. At {t=10} the cat has almost come back fully to node {1}, indeed to the initial quantum state {|1H\rangle}, but not quite. The swing is even closer at {t=38}:

\displaystyle  \mathbf{W}^{38} = \frac{1}{524288}\begin{bmatrix} 522919 & -23939 & -11285 & 23939 & 12654 & 0\\ 23939 & 522919 & 0 & 12654 & -23939 & -11285\\ 12654 & 0 & 522919 & -23939 & -11285 & 23939\\ -23939 & -11285 & 23939 & 522919 & 0 & 12654\\ -11285 & 23939 & 12654 & 0 & 522919 & -23939\\ 0 & 12654 & -23939 & -11285 & 23939 & 522919 \end{bmatrix}

with state {\mathfrak{s}_{38} = [0.99739, 0.04566, 0.02414, -0.04566, -0.02152, 0]^T} and node probabilities {p^{(38)} = [ 0.99687, 0.00267, 0.00046]}. The next probabilities, however, are {p^{(39)} = [ 0.54641, 0.45313, 0.00046]}, which is not close to the first-step distribution {p^{(1)} = [0.5,0.5,0]}.

Most treatments of quantum walks, including a seminal paper by Dorit Aharonov, Andris Ambainis, Julia Kempe, and Umesh Vazirani, and an influential 2003 survey by Kempe, emphasize the dispersion time being linear rather than quadratic as in classical random walks on undirected graphs. Here, the quadratic nature effects quick divergence between close states like {\mathfrak{s}_0} and {\mathfrak{s}_{38}}, which is the signature of dynamical chaos.

The coin {\mathsf{J} = \frac{1}{\sqrt{2}}\begin{bmatrix} 1 & i \\ i & 1\end{bmatrix}} is reputed to make quantum walks smoother. That works to some extent for our three-node graph, but there appears to be no periodic stability, let alone convergence as in the classical case. Simple Python code for trying other coins is available here.

Real or Fiction?

Our main question is simple:

Are these crazy walks real?

The answer is yes if we can engineer a qutrit+qubit system that evolves by repeating {\mathbf{W}}. Simulating such a system throws up a second issue: No one would say that the classical computer executing my Python code is behaving chaotically. We can program {\mathbf{W}} via basic gates available to run on real hardware at IBM Quantum or some other services.

Is quantum hardware that applies powers of {\mathbf{W}} behaving chaotically? Does that pose a concretely physical impediment to maintaining its coherence?

On expecting a yes answer to the walks’ physical existence in the microworld, our questions next try to connect to Smith’s paper:

Does the mathematical existence of chaotic quantum systems on such small scales have macroscopic effects in our world? Such as on the chaotic dynamics of our atmosphere?

Again, the answer should be yes at the level of connecting Brownian motion to quantum processes, or of quantum-walk interpretations of physical systems such as the quantum harmonic oscillator (QHO). So the most specific form of our question is:

Does our (not-so-)simple three-node quantum walk appear as a non-negligible term in a Feynman-style summation needed to understand a macroscopic physical system?

Open Problems

Can the quantum chaos of the Q-tip graph be tamed? See this for success on other graphs including the 3-cycle. There are are affinities between the 3-cycle and Q-tip walks, especially at even timesteps, and the 3-cycle is an option in the Python code.

Smith also has many activities promoting STEM, including her Phyiscs Is Fun website and a book on particle physics. She maintains a blog of her writing and reading, plus much else on her author site. Dick and I were already wondering about women in the community of math/theory bloggers, such as Tanya Khovanova, Sophia Wood, Vi Hart, and Anna Weltman.


[fixed link for Python code, some word changes, fixed entry of J matrix]

5 Comments leave one →
  1. Ali Dasdan permalink
    November 4, 2021 12:45 pm

    the “physics is fun” URL has an extra ” at the end that needs to be removed to work.

    • November 4, 2021 11:36 pm

      Thanks. Actually it was a missing ” at the other end, an error WordPress is IMHO overforgiving of.

  2. January 19, 2022 9:41 pm

    I suspect that the quantum walk “chaos” happens to relate to John Horton Conway’s “fractal” programming using his language called FRACTRAN.
    See the pp. 4, ff. of
    https://www.amazon.co.jp/Problems-Communication-Computation-Thomas-Cover/dp/0387966218
    (I ask you to Google-translate my book-review on it in the link, too.)

Trackbacks

  1. POPL 2022—Et Tu, Brute? | Gödel's Lost Letter and P=NP
  2. Quantum Circuits in the New York Times | Gödel's Lost Letter and P=NP

Leave a Reply

Discover more from Gödel's Lost Letter and P=NP

Subscribe now to keep reading and get access to the full archive.

Continue reading