Blog - network theory (part 12)

This page is a blog article in progress, written by John Baez. To see discussions of this page as it was being written, go to the Azimuth Forum. To see the final polished version, visit the Azimuth Blog.

Last time we proved a version of Noether’s theorem for stochastic mechanics. Now I want to compare that to the more familiar quantum version.

But to do this, I need to say more about the analogy between stochastic mechanics and quantum mechanics. And whenever I try, I get pulled toward explaining some technical issues involving analysis: whether sums converge, whether derivatives exist, and so on. I’ve been trying to avoid such stuff—not because I dislike it, but because I’m afraid *you* might. But the more I put off discussing these issues, the more they fester and make me unhappy. In fact, that’s why it’s taken so long for me to write this post!

So, this time I will gently explore some of these technical issues. But don’t be scared: I’ll *mainly* talk about some simple big ideas. Next time I’ll discuss Noether’s theorem. I hope that by getting the technicalities out of my system, I’ll feel okay about hand-waving whenever I want.

And if you’re an expert on analysis, maybe you can help me with a question.

First, we need to recall the analogy we began sketching in Part 5, and push it a bit further. The idea is that stochastic mechanics differs from quantum mechanics in two big ways:

• First, instead of complex amplitudes, stochastic mechanics uses nonnegative real probabilities. The complex numbers form a ring; the nonnegative real numbers form a mere rig, which is a ‘ri**n**g without **n**egatives’. Rigs are much neglected in the typical math curriculum, but unjustly so: they’re almost as good as rings in many ways, and there are lots of important examples, like the natural numbers $\mathbb{N}$ and the nonnegative real numbers, $[0,\mathrm{\infty})$. For probability theory, we should learn to love rigs.

But there are, alas, situations where we need to subtract probabilities, even when the answer comes out negative: namely when we’re taking the *time derivative* of a probability. So sometimes we need $\mathbb{R}$ instead of just $[0,\mathrm{\infty})$.

• Second, while in quantum mechanics a state is described using a ‘wavefunction’, meaning a complex-valued function obeying

$$\int \mid \psi {\mid}^{2}=1$$

in stochastic mechanics it’s described using a ‘probability distribution’, meaning a nonnegative real function obeying

$$\int \psi =1$$

So, let’s try our best to present the theories in close analogy, while respecting these two differences.

We’ll start with a set $X$ whose points are **states** that a system can be in. Last time I assumed $X$ was a finite set, but this post is so mathematical I might as well let my hair down and assume it’s a measure space. A measure space lets you do integrals, but a finite set is a special case, and then these integrals are just sums. So, I’ll write things like

$$\int f$$

and mean the integral of the function $f$ over the measure space $X$, but if $X$ is a finite set this just means

$$\sum _{x\in X}f(x)$$

Now, I’ve already defined the word ‘state’, but both quantum and stochastic mechanics need a more general concept of state. Let’s call these ‘quantum states’ and ‘stochastic states’:

• In **quantum mechanics**, the system has an **amplitude** $\psi (x)$ of being in any state $x\in X$. These amplitudes are complex numbers with

$$\int \mid \psi {\mid}^{2}=1$$

We call $\psi :X\to \u2102$ obeying this equation a **quantum state**.

• In **stochastic mechanics**, the system has a **probability** $\psi (x)$ of being in any state $x\in X$. These probabilities are nonnegative real numbers with

$$\int \psi =1$$

We call $\psi :X\to [0,\mathrm{\infty})$ obeying this equation a **stochastic state**.

In quantum mechanics we often use this abbreviation:

$$\u27e8\varphi ,\psi \u27e9=\int \overline{\varphi}\psi $$

so that a quantum state has

$$\u27e8\psi ,\psi \u27e9=1$$

Similarly, we could introduce this notation in stochastic mechanics:

$$\u27e8\psi \u27e9=\int \psi $$

so that a stochastic state has

$$\u27e8\psi \u27e9=1$$

But this notation is a bit risky, since angle brackets of this sort often stand for expectation values of observables. So, I’ve been writing $\int \psi $, and I’ll keep on doing this.

In quantum mechanics, $\u27e8\varphi ,\psi \u27e9$ is well-defined whenever both $\varphi $ and $\psi $ live in the vector space

$${L}^{2}(X)=\{\psi :X\to \u2102\phantom{\rule{thickmathspace}{0ex}}:\phantom{\rule{thickmathspace}{0ex}}\int \mid \psi {\mid}^{2}<\mathrm{\infty}\}$$

In stochastic mechanics, $\u27e8\psi \u27e9$ is well-defined whenever $\psi $ lives in the vector space

$${L}^{1}(X)=\{\psi :X\to \mathbb{R}\phantom{\rule{thickmathspace}{0ex}}:\phantom{\rule{thickmathspace}{0ex}}\int \mid \psi \mid <\mathrm{\infty}\}$$

You’ll notice I wrote $\mathbb{R}$ rather than $[0,\mathrm{\infty})$ here. That’s because in some calculations we’ll need functions that take negative values, even though our stochastic states are nonnegative.

A state is a way our system can be. An observable is something we can measure about our system. They fit together: we can measure an observable when our system is in some state. If we repeat this we may get different answers, but there’s a nice formula for average or ‘expected’ answer.

• In quantum mechanics, an **observable** is a self-adjoint operator $A$ on ${L}^{2}(X)$. The **expected value** of $A$ in the state $\psi $ is

$$\u27e8\psi ,A\psi \u27e9$$

Here I’m assuming that we can apply $A$ to $\psi $ and get a new vector $A\psi \in {L}^{2}(X)$. This is automatically true when $X$ is a finite set, but in general we need to be more careful.

• In stochastic mechanics, an **observable** is a real-valued function $A$ on $X$. The **expected value** of $A$ in the state $\psi $ is

$$\int A\psi $$

Here we’re using the fact that we can multiply $A$ and $\psi $ and get a new vector $A\psi \in {L}^{1}(X)$, at least if $A$ is bounded. Again, this is automatic if $X$ is a finite set, but not otherwise.

Besides states and observables, we need ‘symmetries’, which are transformations that map states to states. We use these to describe how our system changes when we wait a while, for example.

• In quantum mechanics, an **isometry** is a linear map $U:{L}^{2}(X)\to {L}^{2}(X)$ such that

$$\u27e8U\varphi ,U\psi \u27e9=\u27e8\varphi ,\psi \u27e9$$

for all $\psi ,\varphi \in {L}^{2}(X)$. If $U$ is an isometry and $\psi $ is a quantum state, then $U\psi $ is again a quantum state.

• In stochastic mechanics, a **stochastic operator** is a linear map $U:{L}^{1}(X)\to {L}^{1}(X)$ such that

$$\int U\psi =\int \psi $$

and

$$\psi \ge 0\phantom{\rule{thickmathspace}{0ex}}\phantom{\rule{thickmathspace}{0ex}}\Rightarrow \phantom{\rule{thickmathspace}{0ex}}\phantom{\rule{thickmathspace}{0ex}}U\psi \ge 0$$

for all $\psi \in {L}^{1}(X)$. If $U$ is stochastic and $\psi $ is a stochastic state, then $U\psi $ is again a stochastic state.

In quantum mechanics we are mainly interested in invertible isometries, which are called **unitary** operators. There are lots of these. There are, however, very few invertible stochastic operators:

**Puzzle 1.** Suppose $X$ is a finite set. Show that every isometry $U:{L}^{2}(X)\to {L}^{2}(X)$ is invertible.

**Puzzle 2.** Suppose $X$ is a finite set. What are the invertible stochastic operators $U:{L}^{1}(X)\to {L}^{1}(X)$?

This is why we usually think of time evolution as being reversible quantum mechanics, but not in stochastic mechanics! In quantum mechanics we often describe time evolution using a ‘1-parameter group’, while in stochastic mechanics we describe it using a 1-parameter *semi*group… meaning that we can run time forwards, but not backwards.

But let’s see how this works in detail!

In quantum mechanics there’s a beautiful relation between observables and symmetries, which goes like this. Suppose that for each time $t$ we want a unitary operator $U(t):{L}^{2}(X)\to {L}^{2}(X)$ that describes time evolution. Then it makes a lot of sense to demand that these operators form a 1-parameter group:

**Definition.** A collection of linear operators U(t) ($t\in \mathbb{R}$) on some vector space forms a **1-parameter group** if

$$U(0)=1$$

and

$$U(s+t)=U(s)U(t)$$

for all $s,t\in \mathbb{R}$.

Note that these conditions force all the operators $U(t)$ to be invertible.

Now suppose our vector space is a Hilbert space, like ${L}^{2}(X)$. Then we call a 1-parameter group a **1-parameter unitary group** if the operators involved are all unitary.

It turns out that 1-parameter unitary groups are either continuous in a certain way, or so pathological that you can’t even prove they exist without the axiom of choice! So, we always focus on the continuous case:

**Definition.** A 1-parameter unitary group is **strongly continuous** if $U(t)\psi $ depends continuously on $t$ for all $\psi $, in this sense:

$${t}_{i}\to t\phantom{\rule{thickmathspace}{0ex}}\phantom{\rule{thickmathspace}{0ex}}\Rightarrow \phantom{\rule{thickmathspace}{0ex}}\phantom{\rule{thickmathspace}{0ex}}\parallel U({t}_{i})\psi -U(t)\psi \parallel \to 0$$

Then we get a classic result proved by Marshall Stone back in the early 1930s. You may not know him, but he was so influential at the University of Chicago during this period that it’s often called the “Stone Age”. And here’s one reason why:

**Stone's Theorem.** There is a one-to-one correspondence between strongly continuous 1-parameter unitary groups on a Hilbert space and self-adjoint operators on that Hilbert space, given as follows. Given a strongly continuous 1-parameter unitary group $U(t)$ we can always write

$$U(t)=\mathrm{exp}(-itH)$$

for a unique self-adjoint operator $H$. Conversely, any self-adjoint operator determines a strongly continuous 1-parameter group this way. For all vectors $\psi $ for which $H\psi $ is well-defined, we have

$${\frac{d}{dt}U(t)\psi \mid}_{t=0}=-iH\psi $$

Moreover, for any of these vectors, if we set

$$\psi (t)=\mathrm{exp}(-itH)\psi $$

we have

$$\frac{d}{dt}\psi (t)=-iH\psi (t)$$

When $U(t)=\mathrm{exp}(-itH)$ describes the evolution of a system in time, $H$ is is called the **Hamiltonian**, and it has the physical meaning of ‘energy’. The equation I just wrote down is then called **Schrödinger's equation**.

So, simply put, in quantum mechanics we have a correspondence between observables and nice one-parameter groups of symmetries. Not surprisingly, our favorite observable, energy, corresponds to our favorite symmetry: time evolution!

However, if you were paying attention, you noticed that I carefully avoided explaining how we define $\mathrm{exp}(-itH)$. I didn’t even say what a self-adjoint operator is. This is where the technicalities come in: they arise when $H$ is unbounded, and not defined on all vectors in our Hilbert space.

Luckily, these technicalities evaporate for *finite-dimensional* Hilbert spaces, such as ${L}^{2}(X)$ for a finite set $X$. Then we get:

**Stone’s Theorem (Baby Version).** Suppose we are given a *finite-dimensional* Hilbert space. In this case, a linear operator $H$ on this space is self-adjoint iff it’s defined on the whole space and

$$\u27e8\varphi ,H\psi \u27e9=\u27e8H\varphi ,\psi \u27e9$$

for all vectors $\varphi ,\psi $. Given a strongly continuous 1-parameter unitary group $U(t)$ we can always write

$$U(t)=\mathrm{exp}(-itH)$$

for a unique self-adjoint operator $H$, where

$$\mathrm{exp}(-itH)\psi =\sum _{n=0}^{\mathrm{\infty}}\frac{(-itH{)}^{n}}{n!}\psi $$

with the sum converging for all $\psi $. Conversely, any self-adjoint operator on our space determines a strongly continuous 1-parameter group this way. For all vectors $\psi $ in our space we then have

$${\frac{d}{dt}U(t)\psi \mid}_{t=0}=-iH\psi $$

and if we set

$$\psi (t)=\mathrm{exp}(-itH)\psi $$

we have

$$\frac{d}{dt}\psi (t)=-iH\psi (t)$$

We’ve seen that in quantum mechanics, time evolution is usually described by a 1-parameter group of operators that comes from an observable: the Hamiltonian. Stochastic mechanics is different!

First, since stochastic operators aren’t usually invertible, we typically describe time evolution by a mere ‘semigroup’:

**Definition.** A collection of linear operators $U(t)$ ($t\in [0,\mathrm{\infty})$) on some vector space forms a **1-parameter semigroup** if

$$U(0)=1$$

and

$$U(s+t)=U(s)U(t)$$

for all $s,t\ge 0$.

Now suppose this vector space is ${L}^{1}(X)$ for some measure space $X$. We want to focus on the case where the operators $U(t)$ are stochastic and depend continuously on $t$ in the same sense we discussed earlier.

**Definition.** A 1-parameter strongly continuous semigroup of stochastic operators $U(t):{L}^{1}(X)\to {L}^{1}(X)$ is called a **Markov semigroup**.

What’s the analogue of Stone’s theorem for Markov semigroups? I don’t know a fully satisfactory answer! If you know, please tell me.

Later I’ll say what I *do* know—I’m not *completely* clueless—but for now let’s look at the ‘baby’ case where $X$ is a finite set. Then the story is neat and complete:

**Theorem.** Suppose we are given a *finite* set $X$. In this case, a linear operator $H$ on ${L}^{1}(X)$ is **infinitesimal stochastic** iff it’s defined on the whole space,

$$\int H\psi =0$$

for all $\psi \in {L}^{1}(X)$, and the matrix of $H$ in terms of the obvious basis obeys

$${H}_{ij}\ge 0$$

for all $j\ne i$. Given a Markov semigroup $U(t)$ on ${L}^{1}(X)$, we can always write

$$U(t)=\mathrm{exp}(tH)$$

for a unique infinitesimal stochastic operator $H$, where

$$\mathrm{exp}(tH)\psi =\sum _{n=0}^{\mathrm{\infty}}\frac{(tH{)}^{n}}{n!}\psi $$

with the sum converging for all $\psi $. Conversely, any infinitesimal operator on our space determines a Markov semigroup this way. For all $\psi \in {L}^{1}(X)$ we then have

$${\frac{d}{dt}U(t)\psi \mid}_{t=0}=H\psi $$

and if we set

$$\psi (t)=\mathrm{exp}(tH)\psi $$

we have the **master equation**:

$$\frac{d}{dt}\psi (t)=H\psi (t)$$

In short, time evolution in stochastic mechanics is a lot like time evolution in quantum mechanics, except it’s typically not invertible, and the Hamiltonian is typically *not an observable*.

Why not? Because we defined an observable to be a function $A:X\to \mathbb{R}$. We can think of this as giving an operator on ${L}^{1}(X)$, namely the operator of multiplication by $A$. That’s a nice trick, which we used to good effect last time. However, at least when $X$ is a finite set, this operator will be diagonal in the obvious basis consisting of functions that equal 1 one point of $X$ and zero elsewhere. So, it can only be infinitesimal stochastic if it’s zero!

**Puzzle 3.** If $X$ is a finite set, show that any operator on ${L}^{1}(X)$ that’s both diagonal and infinitesimal stochastic must be zero.

I’ve now told you everything you really need to know… but not everything I want to say. What happens when $X$ is not a finite set? What are Markov semigroups like then? I can’t abide letting this obvious question go unresolved! Unfortunately I only know a partial answer.

We can get a certain distance using the Hille-Yosida theorem, which is much more general.

**Definition.** A Banach space is vector space with a norm such that any Cauchy sequence converges.

Examples include Hilbert spaces like ${L}^{2}(X)$ for any measure space, but also other spaces like ${L}^{1}(X)$ for any measure space!

**Definition.** If $V$ is a Banach space, a 1-parameter semigroup of operators $U(t):V\to V$ is called a **contraction semigroup** if it’s strongly continuous and

$$\parallel U(t)\psi \parallel \le \parallel \psi \parallel $$

for all $t\ge 0$ and all $\psi \in V$.

Examples include strongly continuous 1-parameter unitary groups, but also Markov semigroups!

**Puzzle 4.** Show any Markov semigroup is a contraction semigroup.

The Hille–Yosida theorem generalizes Stone’s theorem to contraction semigroups. In my misspent youth, I spent a lot of time carrying around Yosida’s book *Functional Analysis*. Furthermore, Einar Hille was the advisor of my thesis advisor, Irving Segal. Segal generalized the Hille–Yosida theorem to *nonlinear* operators. I used this generalization a lot back when I studied nonlinear partial differential equations. So, I feel compelled to tell you this theorem:

**Hille-Yosida Theorem.** Given a contraction semigroup $U(t)$ we can always write

$$U(t)=\mathrm{exp}(tH)$$

for some closed and densely defined operator $H$ such that $H-\lambda I$ has a bounded inverse for all $\lambda >0$ and

$$\parallel (H-\lambda I{)}^{-1}\psi \parallel \le \frac{1}{\lambda}\parallel \psi \parallel$$

for all $\psi \in v$. Conversely, any such operator determines a strongly continuous 1-parameter group. For all vectors $\psi $ for which $H\psi $ is well-defined, we have

$${\frac{d}{dt}U(t)\psi \mid}_{t=0}=H\psi $$

Moreover, for any of these vectors, if we set

$$\psi (t)=U(t)\psi $$

we have

$$\frac{d}{dt}\psi (t)=H\psi (t)$$

If you like, you can take the stuff at the end of this theorem to be what we mean by saying $U(t)=\mathrm{exp}(tH)$.

But now suppose $V={L}^{1}(X)$. What *extra conditions* on $H$ are necessary and sufficient for $\mathrm{exp}(tH)$ to be a Markov semigroup? In other words, what’s a definition of ‘infinitesimal stochastic operator’ that’s suitable not only when $X$ is a finite set, but an arbitrary measure space?

I asked this question on Mathoverflow a few months ago, and so far the answers have not been completely satisfactory. Some people mentioned the Hille–Yosida theorem, which is surely a step in the right direction, but not the full answer.

Others discussed the special case when $\mathrm{exp}(tH)$ extends to a bounded self-adjoint operator on ${L}^{2}(X)$. When $X$ is a finite set, this special case happens precisely when the matrix ${H}_{ij}$ is *symmetric*: the probability of hopping from $j$ to $i$ equals the probability of hopping from $i$ to $j$. This is a fascinating special case, not least because when $H$ is both infinitesimal stochastic *and* self-adjoint, we can use it as a Hamiltonian for both stochastic mechanics *and* quantum mechanics! However, it’s just a special case.

After grabbing people by the collar and insisting that I wanted to know the answer to the question I actually asked—not some vaguely similar question— the best answer seems to be Martin Gisser’s reference to this book:

• Zhi-Ming Ma and Michael Röckner, *Introduction to the Theory of (Non-Symmetric) Dirichlet Forms*, Springer, Berlin, 1992.

However, as best as I can tell, this does *not* answer my question in general, but only when the skew-symmetric part of $H$ is dominated (in a certain sense) by the symmetric part.

So, I’m stuck on this front, but that needn’t bring the whole project to a halt. We’ll just sidestep this question.

For a good well-rounded introduction to Markov semigroups and what they’re good for, try:

• Ryszard Rudnicki, Katarzyna Pichór and Marta Tyran-Kamínska, Markov semigroups and their applications.

category: blog