I’m happy to use this bully pulpit to advertise that the following paper has been deemed “probably not terrible”, i.e., published.

Here’s the figure^{ a } and caption:

It is my highly unusual opinion that identifying a definition for the branches in the wavefunction is the most conceptually important problem is physics. The reasoning is straightforward: (1) quantum mechanics is the most profound thing we know about the universe, (2) the measurement process is at the heart of the weirdness, and (3) the critical roadblock to analysis is a definition of what we’re talking about. (Each step is of course highly disputed, and I won’t defend the reasoning here.) In my biased opinion, the paper represents the closest yet anyone has gotten to giving a mathematically precise definition.

On the last page of the paper, I speculate on the possibility that branch finding may have practical (!) applications for speeding up numerical simulations of quantum many-body systems using matrix-product states (MPS), or tensor networks in general. The rough idea is this: Generic quantum systems are exponentially hard to simulate, but classical systems (even stochastic ones) are not. A definition of branches would identify *which degrees of freedom* of a quantum system could be accurately simulated classically, and when. Although classical computational transitions are understood in many certain special cases, our macroscopic observations of the real world strongly suggest that *all* systems we study admit classical descriptions on large enough scales.^{ b } There is reason to think a general (abstract) description of the quantum-classical transition is possible, and would allow us to go beyond special cases.

In the rest this blog post I’m going to construct a simple MPS state featuring branches with records that are arbitrarily redundant. The state’s key features will be that it is translationally invariant and has finite correlation length. Translational invariance makes the state much easier to study and compare with the literature, and is a property shared with the simple inflationary model I’m studying with Elliot Nelson. Finite correlation length eliminates the trivial solution of a generalized GHZ state, guarantees that there are an infinite number of branches (if there are any at all), and is in some sense more natural.

#### An MPS with branches

Our strategy will be to build the state up from a *classical* probability distribution that already has the key features. Let , be a classical state of a finite 1D lattice of bits (), and consider the canonical ensemble starting with a uniform distribution and adding an energetic penalty of for misaligned nearest-neighbors^{ c }:

(1)

with . It’s not hard to check that the classical expectation value obeys and

(2)

If we define our correlation length by the asymptotic behavior , then for large (large ) we have

(3)

The idea is to define our quantum state as a superposition of configurations of a spin chain weighted by this distribution: . Contiguous regions of size will be highly correlated. The coarse-grained variables that are being recorded are something like “whether this local region is mostly or mostly ”, but we don’t necessarily want each qubit to have a perfect record of this information. Rather, we choose

(4)

for some fixed (typically small) angle . For one can’t reliably infer the classical variable from a measurement on a single spin, but one can make an arbitrarily reliable inference by measuring lots of them: . If we define our “recording length” by the asymptotic behavior , then for small we have

(5)

Thus our final (unnormalized) state is

(6)

for arbitrary positive parameters and . The state inherits from the classical probability distribution the properties of finite correlation length and translational invariance. Any contiguous set of spins is very likely to all be in the same state, and by measuring a subset of spins one can distinguish between the cases and with high reliability. Therefore in the limit it makes sense to define the redundancy ; this is the approximate number of disjoint records that are available about each branch, where there is a (binary) branch density of .

Now lets turn this into a matrix-product state.^{ d } That is, it should take the form^{ e }^{ f }

(7)

where the range over the local Hilbert spaces of the qubits and the are the “bonds” representing contraction of matrices of as-yet unspecified dimension (depicted in the figure as lines connecting the ‘s horizontally). The trick is to pretend that our classical probability distribution was achieved by starting with a quantum state of appropriately weighted bras,

(8)

and contracting (i.e., multiplying from the right to inner-product over the fictional bras) with a specially designed state that picks out the correct conditional states of the original Hilbert space:

(9)

The object

(10)

then takes the form of (7) and is equal to so long as we choose to satisfy

(11)

or, more cleanly,

(12)

*[I think Martin Ganahl and Guifre Vidal for discussion.]*

### Footnotes

(↵ returns to text)

- The editor tried to convince me that this figure appeared on the cover for purely aesthetic reasons and this does not mean my letter is the best thing in the issue…but I know better!↵
- Note that whether certain degrees of freedom admit a classical effective description is a
*computational*question. The existence of such a description is incompatible with them exhibiting quantum supremacy, but it is compatible with them encoding circumstantial evidence (at the macro scale) for quantum mechanics (at the micro scale).↵ - We assume periodic boundary conditions: ↵
- This will not be a general algorithm for building MPSs since it exploits the special form of the state (6) we’re trying to build, which allows us to identify the ‘s with the ‘s.↵
- Again, periodic boundary conditions↵
- The tikzpicture code was obtained from Piotr Migdał.↵

I have a stupid(?) question. The criteria you give in the paper doesn’t involve time-evolution in any way. Intuitively, it seems that the reason we can treat branches separately is that they will evolve roughly independently of each other. We have U(B1 + B2) = U(B1) + U(B2), at least on the time scales we care about. Would a satisfactory abstract criteria for identifying branches make use of the form of the time evolution operator?

You’re right that the criteria for defining branches I give in the paper doesn’t make explicit reference to time evolution. Rather, the criteria for branches at a given time is based only on the wavefunction at that time, so that, in principle, the branch decomposition could change dramatically from one moment in time to another.

Our

expectation(or hope) is that the number of branches increases monotonically in time and that, furthermore, the branches at an earlier time, when evolved forward, are just a coarse-graining of the branches at a later time. Here, “coarse-graining” reflects that fact that individual branches may subdivide (e.g., when a measurement is made), but if you add the sub-branches together you should recover the parent branches (suitably time-evolved). More precisely, we expect that if the branch decomposition is at time , and is at some later time , then , where the form a partition of the set , i.e., and for . Let us call this property “branch fine-graining under time evolution” (BFGUTE)BFGUTE is merely an unproven desiderata. Indeed, one way to describe the way in which our understanding of quantum mechanics is incomplete is that we have never proven BFGUTE, which is basically the statement that we haven’t proven the Copenhagen and Everett interpretations equivalent. The first step in proving BFGUTE is to find a precise definition of wavefunction branches. Although it would be very nice if we could find a definition of branches that automatically implied BFGUTE, a little thought shows that the BFGUTE property is not enough, on its own, to define branches.

I’m not sure if that answered your questions, but maybe it clarified things enough that you could re-ask?

Right, so to argue that the branches evolve independently, one must first define what branches are. That makes sense.

I guess my question might be, is it necessarily the case that the branch decomposition can be derived from the tensor product structure/ entanglement alone? It seems that more structure is needed, as in the preferred length scale in the paper. The time evolution operator is a natural source of such a structure, e.g. if the qubits were fixed in some grid such that only neighbours can influence each other, that would be reflected in the Hamiltonian. So then the ‘correct’ branch decomposition would be a function of the Hamiltonian/additional structure, not just the state itself. Is this basically correct?

Great question. It is not obvious that we should define branches from the tensor structure alone, but here are two considerations:

(1) Strategically, we want to start with the accepted recipe for measurement. As shown by Zurek (discussed here), measurements are really about

amplification, and the best way I know how to formalize amplification is in terms ofcopies of information. Although I’m very interested in finding another formalization involving time (e.g., divergent information flow?), the simplest mathematical criteria seems to involve just identifying correlated information at some fixed time after the measurement.(2) Intuitively, if you hand me the macroscopic wavefunction describing Schrödinger’s cat, it seems to me that we can identify the branches without any reference to the historical evolution that generated that state. All I need to do is just look at the entanglement structure of a single time-slice and it’s obvious.

That said, it’s conceivable a toy model with the following features could be found: On a single late-time slice, there are two incompatible candidate branch decompositions with equivalent spatial entanglement structure, but earlier time evolution unambiguously picks out one branch decomposition as the “correct” one. (This might be based on one decomposition having the BFGUTE property.)

I have spent some time trying to find an additional criterion — grounded in the Hamiltonian, the lattice spacing, or the time evolution — that would unambiguously pick out a unique decomposition (unlike in the paper), but everything I tried was ugly/ad-hoc. Obviously, you can just declare a criterion, but I’d like something that was as compelling as the idea (used in the paper) that anything that deserves to be called a measurement must make at least three copies.

Thanks. Do you know of any systems, toy or real, where we can identify a naturally emerging branch structure?

I would be interested in trying to run a simulation of a simple system with some of the features that seem important to branching in the real world, such as locality. A QCA for instance. Seeing such a system evolve, and hopefully form branch-like structures, might be helpful for trying to figure out what ‘branches’ really are (although of course we could only do it for a handful of particles/sites) Has anybody tried to do anything like this? Do you know what sort of software/methods could be used to implement such a simulation?

This prompted me to take an old email where I listed some models of branching and convert it into a blog post here.

I think the most promising model would be a 1D chain of spins with some sort of local interactions. See especially the Brun-Halliwell model [arXiv:quant-ph/9601004]. The idea is that hydrodynamic variables (i.e., locally conserved densities) are excellent candidates for variables that would be redundantly recorded, and would follow quasiclassical trajectories on short timescales.

I don’t know anyone who has tried to do a numerical simulation to find branches based on spatial entanglement structure (rather than by imposing a system-environment distinction by hand). This is something I’ve wanted to do for a while. Email me if you want to discuss more: jessriedel@gmail.com

My technical experience here doesn’t go beyond some simple Mathematica models, but I imagine this could only feasibly be done for a reasonably large spin chain by using a tensor network (in this case a matrix-product state) .

A quantum cellular automata would be very interesting. One drawback would (I think) be that there wouldn’t be quasi-classical evolution, which would be ideal for convincingly deriving the appearance of classical trajectories.