<?xml version="1.0" encoding="UTF-8" standalone="yes"?><oembed><version><![CDATA[1.0]]></version><provider_name><![CDATA[Azimuth]]></provider_name><provider_url><![CDATA[https://johncarlosbaez.wordpress.com]]></provider_url><author_name><![CDATA[John Baez]]></author_name><author_url><![CDATA[https://johncarlosbaez.wordpress.com/author/johncarlosbaez/]]></author_url><title><![CDATA[Noether&#8217;s Theorem: Quantum vs&nbsp;Stochastic]]></title><type><![CDATA[link]]></type><html><![CDATA[<p><i>guest post by <b><a href="http://www.azimuthproject.org/azimuth/show/Ville+Bergholm">Ville Bergholm</a></b></i></p>
<p>In 1915 Emmy Noether discovered an important connection between the symmetries of a system and its conserved quantities. Her result has become a staple of modern physics and is known as Noether&#8217;s theorem.   </p>
<div align="center"><img src="https://i0.wp.com/upload.wikimedia.org/wikipedia/commons/e/e5/Noether.jpg" width="250" alt="Photo of Emmy Noether" /></div>
<p>The theorem and its generalizations have found particularly wide use in quantum theory.  Those of you following the <a href="http://math.ucr.edu/home/baez/networks/">Network Theory</a> series here on Azimuth might recall <a href="http://math.ucr.edu/home/baez/networks/networks_11.html">Part 11</a> where John Baez and Brendan Fong proved a version of Noether&#8217;s theorem for stochastic systems.  Their result is now published here:</p>
<p>&bull; John Baez and Brendan Fong, <a href="http://arxiv.org/abs/1203.2035">A Noether theorem for stochastic mechanics</a>, <i>J. Math. Phys.</i> 54:013301 (2013). </p>
<p>One goal of the network theory series here on Azimuth has been to merge ideas appearing in quantum theory with other disciplines.  John and Brendan proved their stochastic version of Noether&#8217;s theorem by exploiting &#8216;stochastic mechanics&#8217; which was formulated in the network theory series to mathematically resemble quantum theory.    Their result, which we will outline below, was different than what would be expected in quantum theory, so it is interesting to try to figure out <i>why</i>.  </p>
<p>Recently <a href="http://www.azimuthproject.org/azimuth/show/Jacob+Biamonte">Jacob Biamonte</a>, <a href="http://www.azimuthproject.org/azimuth/show/Mauro+Faccin">Mauro Faccin</a> and myself have been working to try to get to the bottom of these differences.  What we&#8217;ve done is prove a version of Noether&#8217;s theorem for Dirichlet operators.  As you may recall from Parts <a href="http://math.ucr.edu/home/baez/networks/networks_16.html">16</a> and <a href="http://math.ucr.edu/home/baez/networks/networks_20.html">20</a> of the network theory series, these are the operators that generate both stochastic <i>and</i> quantum processes.   In the language of the series, they lie in the intersection of stochastic and quantum mechanics.  So, they are a subclass of the infinitesimal stochastic operators considered in John and Brendan&#8217;s work. </p>
<p>The extra structure of Dirichlet operators&mdash;compared with the wider class of infinitesimal stochastic operators&mdash;provided a handle for us to dig a little deeper into understanding the intersection of these two theories.  By the end of this article, astute readers will be able to prove that Dirichlet operators generate <a href="http://en.wikipedia.org/wiki/Doubly_stochastic_matrix">doubly stochastic</a> processes.  </p>
<p>Before we get into the details of our proof, let&#8217;s recall first how conservation laws work in quantum mechanics, and then contrast this with what John and Brendan discovered for stochastic systems.  (For a more detailed comparison between the stochastic and quantum versions of the theorem, see <a href="http://math.ucr.edu/home/baez/networks/networks_13.html">Part 13</a> of the network theory series.)</p>
<h3>The quantum case</h3>
<p>I&#8217;ll assume you&#8217;re familiar with quantum theory, but let&#8217;s start with a few reminders.</p>
<p>In standard quantum theory, when we have a closed system with <img src='https://s0.wp.com/latex.php?latex=n&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='n' title='n' class='latex' /> states, the unitary time evolution of a state <img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi%28t%29%5Crangle&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi(t)&#92;rangle' title='|&#92;psi(t)&#92;rangle' class='latex' /> is generated by a self-adjoint <img src='https://s0.wp.com/latex.php?latex=n+%5Ctimes+n&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='n &#92;times n' title='n &#92;times n' class='latex' /> matrix <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> called the Hamiltonian.  In other words, <img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi%28t%29%5Crangle&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi(t)&#92;rangle' title='|&#92;psi(t)&#92;rangle' class='latex' /> satisfies Schr&ouml;dinger&#8217;s equation:</p>
<p><img src='https://s0.wp.com/latex.php?latex=i+%5Chbar+%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%7C%5Cpsi%28t%29+%5Crangle+%3D+H+%7C%5Cpsi%28t%29+%5Crangle.+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='i &#92;hbar &#92;displaystyle{&#92;frac{d}{d t}} |&#92;psi(t) &#92;rangle = H |&#92;psi(t) &#92;rangle. ' title='i &#92;hbar &#92;displaystyle{&#92;frac{d}{d t}} |&#92;psi(t) &#92;rangle = H |&#92;psi(t) &#92;rangle. ' class='latex' /></p>
<p>The state of a system starting off at time zero in the state <img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi_0+%5Crangle&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi_0 &#92;rangle' title='|&#92;psi_0 &#92;rangle' class='latex' /> and evolving for a time <img src='https://s0.wp.com/latex.php?latex=t&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='t' title='t' class='latex' /> is then given by</p>
<p><img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi%28t%29+%5Crangle+%3D+e%5E%7B-i+t+H%7D%7C%5Cpsi_0+%5Crangle.+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi(t) &#92;rangle = e^{-i t H}|&#92;psi_0 &#92;rangle. ' title='|&#92;psi(t) &#92;rangle = e^{-i t H}|&#92;psi_0 &#92;rangle. ' class='latex' /></p>
<p>The observable properties of a quantum system are associated with self-adjoint operators.   In the state <img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi+%5Crangle%2C&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi &#92;rangle,' title='|&#92;psi &#92;rangle,' class='latex' /> the expected value of the observable associated to a self-adjoint operator <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> is</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Clangle+O+%5Crangle_%7B%5Cpsi%7D+%3D+%5Clangle+%5Cpsi+%7C+O+%7C+%5Cpsi+%5Crangle+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;langle O &#92;rangle_{&#92;psi} = &#92;langle &#92;psi | O | &#92;psi &#92;rangle ' title='&#92;langle O &#92;rangle_{&#92;psi} = &#92;langle &#92;psi | O | &#92;psi &#92;rangle ' class='latex' /></p>
<p>This expected value is constant in time for all states if and only if <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> <a href="http://en.wikipedia.org/wiki/Commuting_matrices">commutes</a> with the Hamiltonian <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' />:</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5BO%2C+H%5D+%3D+0+%5Cquad+%5Ciff+%5Cquad+%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O+%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0+%5Cquad+%5Cforall+%5C%3A+%7C%5Cpsi_0+%5Crangle%2C+%5Cforall+t.+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O, H] = 0 &#92;quad &#92;iff &#92;quad &#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = 0 &#92;quad &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t. ' title='[O, H] = 0 &#92;quad &#92;iff &#92;quad &#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = 0 &#92;quad &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t. ' class='latex' /></p>
<p>In this case we say <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> is a &#8216;conserved quantity&#8217;.  The fact that we have two equivalent conditions for this is a quantum version of Noether&#8217;s theorem!</p>
<h3>The stochastic case</h3>
<p>In stochastic mechanics, the story changes a bit.  Now a state <img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi%28t%29%5Crangle&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi(t)&#92;rangle' title='|&#92;psi(t)&#92;rangle' class='latex' /> is a probability distribution: a vector with <img src='https://s0.wp.com/latex.php?latex=n&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='n' title='n' class='latex' /> nonnegative components that sum to 1.  Schr&ouml;dinger&#8217;s equation gets replaced by the <b>master equation</b>:</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%7C%5Cpsi%28t%29+%5Crangle+%3D+H+%7C%5Cpsi%28t%29+%5Crangle+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;frac{d}{d t}} |&#92;psi(t) &#92;rangle = H |&#92;psi(t) &#92;rangle ' title='&#92;displaystyle{&#92;frac{d}{d t}} |&#92;psi(t) &#92;rangle = H |&#92;psi(t) &#92;rangle ' class='latex' /></p>
<p>If we start with a probability distribution <img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi_0+%5Crangle&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi_0 &#92;rangle' title='|&#92;psi_0 &#92;rangle' class='latex' /> at time zero and evolve it according to this equation, at any later time have </p>
<p><img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi%28t%29%5Crangle+%3D+e%5E%7Bt+H%7D+%7C%5Cpsi_0+%5Crangle.+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi(t)&#92;rangle = e^{t H} |&#92;psi_0 &#92;rangle. ' title='|&#92;psi(t)&#92;rangle = e^{t H} |&#92;psi_0 &#92;rangle. ' class='latex' /></p>
<p>We want this always be a probability distribution.  To ensure that this is so, the Hamiltonian <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> must be <b>infinitesimal stochastic</b>: that is, a real-valued <img src='https://s0.wp.com/latex.php?latex=n+%5Ctimes+n&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='n &#92;times n' title='n &#92;times n' class='latex' /> matrix where the off-diagonal entries are nonnegative and the entries of each column sum to zero.  It no longer needs to be self-adjoint!</p>
<p>When <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is infinitesimal stochastic, the operators <img src='https://s0.wp.com/latex.php?latex=e%5E%7Bt+H%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='e^{t H}' title='e^{t H}' class='latex' /> map the set of probability distributions to itself whenever <img src='https://s0.wp.com/latex.php?latex=t+%5Cge+0%2C&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='t &#92;ge 0,' title='t &#92;ge 0,' class='latex' /> and we call this family of operators a <b>continuous-time Markov process</b>, or more precisely a <b>Markov semigroup</b>.  </p>
<p>In stochastic mechanics, we say an <b>observable</b> <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> is a real diagonal <img src='https://s0.wp.com/latex.php?latex=n+%5Ctimes+n&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='n &#92;times n' title='n &#92;times n' class='latex' /> matrix, and its expected value is given by</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Clangle+O%5Crangle_%7B%5Cpsi%7D+%3D+%5Clangle+%5Chat%7BO%7D+%7C+%5Cpsi+%5Crangle+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;langle O&#92;rangle_{&#92;psi} = &#92;langle &#92;hat{O} | &#92;psi &#92;rangle ' title='&#92;langle O&#92;rangle_{&#92;psi} = &#92;langle &#92;hat{O} | &#92;psi &#92;rangle ' class='latex' /></p>
<p>where <img src='https://s0.wp.com/latex.php?latex=%5Chat%7BO%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;hat{O}' title='&#92;hat{O}' class='latex' /> is the vector built from the diagonal entries of <img src='https://s0.wp.com/latex.php?latex=O.&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O.' title='O.' class='latex' />  More concretely, </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Clangle+O%5Crangle_%7B%5Cpsi%7D+%3D+%5Cdisplaystyle%7B+%5Csum_i+O_%7Bi+i%7D+%5Cpsi_i+%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;langle O&#92;rangle_{&#92;psi} = &#92;displaystyle{ &#92;sum_i O_{i i} &#92;psi_i }' title='&#92;langle O&#92;rangle_{&#92;psi} = &#92;displaystyle{ &#92;sum_i O_{i i} &#92;psi_i }' class='latex' /></p>
<p>where <img src='https://s0.wp.com/latex.php?latex=%5Cpsi_i&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;psi_i' title='&#92;psi_i' class='latex' /> is the <img src='https://s0.wp.com/latex.php?latex=i&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='i' title='i' class='latex' />th component of the vector <img src='https://s0.wp.com/latex.php?latex=%7C%5Cpsi%5Crangle.&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;psi&#92;rangle.' title='|&#92;psi&#92;rangle.' class='latex' /></p>
<p>Here is a version of Noether&#8217;s theorem for stochastic mechanics:</p>
<p><b>Noether&#8217;s Theorem for Markov Processes (Baez&ndash;Fong).</b>  Suppose <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is an infinitesimal stochastic operator and <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> is an observable.  Then </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5BO%2CH%5D+%3D0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O,H] =0 ' title='[O,H] =0 ' class='latex' /></p>
<p>if and only if </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O+%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = 0 ' title='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = 0 ' class='latex' /></p>
<p>and </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D%5Clangle+O%5E2+%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;frac{d}{d t}}&#92;langle O^2 &#92;rangle_{&#92;psi(t)} = 0 ' title='&#92;displaystyle{&#92;frac{d}{d t}}&#92;langle O^2 &#92;rangle_{&#92;psi(t)} = 0 ' class='latex' /></p>
<p>for all <img src='https://s0.wp.com/latex.php?latex=t+%5Cge+0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='t &#92;ge 0' title='t &#92;ge 0' class='latex' /> and all <img src='https://s0.wp.com/latex.php?latex=%5Cpsi%28t%29&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;psi(t)' title='&#92;psi(t)' class='latex' /> obeying the master equation.   &nbsp;  &#9608;</p>
<p>So, just as in quantum mechanics, whenever <img src='https://s0.wp.com/latex.php?latex=%5BO%2CH%5D%3D0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O,H]=0' title='[O,H]=0' class='latex' /> the expected value of <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> will be conserved:</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O&#92;rangle_{&#92;psi(t)} = 0' title='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O&#92;rangle_{&#92;psi(t)} = 0' class='latex' /></p>
<p>for any <img src='https://s0.wp.com/latex.php?latex=%5Cpsi_0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;psi_0' title='&#92;psi_0' class='latex' /> and all <img src='https://s0.wp.com/latex.php?latex=t+%5Cge+0.&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='t &#92;ge 0.' title='t &#92;ge 0.' class='latex' />  However, John and Brendan saw that&mdash;unlike in quantum mechanics&mdash;you need more than just the expectation value of the observable <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> to be constant to obtain the equation <img src='https://s0.wp.com/latex.php?latex=%5BO%2CH%5D%3D0.&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O,H]=0.' title='[O,H]=0.' class='latex' />  You really need both</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O&#92;rangle_{&#92;psi(t)} = 0 ' title='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O&#92;rangle_{&#92;psi(t)} = 0 ' class='latex' /></p>
<p>together with </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O%5E2%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O^2&#92;rangle_{&#92;psi(t)} = 0 ' title='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O^2&#92;rangle_{&#92;psi(t)} = 0 ' class='latex' /></p>
<p>for all initial data <img src='https://s0.wp.com/latex.php?latex=%5Cpsi_0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;psi_0' title='&#92;psi_0' class='latex' /> to be sure that <img src='https://s0.wp.com/latex.php?latex=%5BO%2CH%5D%3D0.&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O,H]=0.' title='[O,H]=0.' class='latex' /> </p>
<p>So it&#8217;s a bit subtle, but symmetries and conserved quantities have a rather different relationship than they do in quantum theory.  </p>
<h3>A Noether theorem for Dirichlet operators</h3>
<p>But what if the infinitesimal generator of our Markov semigroup is <i>also self-adjoint?</i>   In other words, what if <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is both an infinitesimal stochastic matrix but also its own transpose: <img src='https://s0.wp.com/latex.php?latex=H+%3D+H%5E%5Ctop&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H = H^&#92;top' title='H = H^&#92;top' class='latex' />?  Then it&#8217;s called a <b>Dirichlet operator</b>&#8230; and we found that in this case, we get a stochastic version of Noether&#8217;s theorem that more closely resembles the usual quantum one:</p>
<p><b>Noether&#8217;s Theorem for Dirichlet Operators.</b> If <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is a Dirichlet operator and <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> is an observable, then</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5BO%2C+H%5D+%3D+0+%5Cquad+%5Ciff+%5Cquad+%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O+%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0+%5Cquad+%5Cforall+%5C%3A+%7C%5Cpsi_0+%5Crangle%2C+%5Cforall+t+%5Cge+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O, H] = 0 &#92;quad &#92;iff &#92;quad &#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = 0 &#92;quad &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t &#92;ge 0 ' title='[O, H] = 0 &#92;quad &#92;iff &#92;quad &#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = 0 &#92;quad &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t &#92;ge 0 ' class='latex' /></p>
<p><b>Proof.</b> The <img src='https://s0.wp.com/latex.php?latex=%5CRightarrow&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;Rightarrow' title='&#92;Rightarrow' class='latex' /> direction is easy to show, and it follows from John and Brendan&#8217;s theorem.  The point is to show the <img src='https://s0.wp.com/latex.php?latex=%5CLeftarrow&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;Leftarrow' title='&#92;Leftarrow' class='latex' /> direction.  Since <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is self-adjoint, we may use a spectral decomposition:</p>
<p><img src='https://s0.wp.com/latex.php?latex=H+%3D+%5Cdisplaystyle%7B+%5Csum_k+E_k+%7C%5Cphi_k+%5Crangle+%5Clangle+%5Cphi_k+%7C%7D+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H = &#92;displaystyle{ &#92;sum_k E_k |&#92;phi_k &#92;rangle &#92;langle &#92;phi_k |} ' title='H = &#92;displaystyle{ &#92;sum_k E_k |&#92;phi_k &#92;rangle &#92;langle &#92;phi_k |} ' class='latex' /></p>
<p>where <img src='https://s0.wp.com/latex.php?latex=%5Cphi_k&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;phi_k' title='&#92;phi_k' class='latex' /> are an orthonormal basis of eigenvectors, and <img src='https://s0.wp.com/latex.php?latex=E_k&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='E_k' title='E_k' class='latex' /> are the corresponding eigenvalues.  We then have:</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O+%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+%5Clangle+%5Chat%7BO%7D+%7C+H+e%5E%7Bt+H%7D+%7C%5Cpsi_0+%5Crangle+%3D+0+%5Cquad+%5Cforall+%5C%3A+%7C%5Cpsi_0+%5Crangle%2C+%5Cforall+t+%5Cge+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = &#92;langle &#92;hat{O} | H e^{t H} |&#92;psi_0 &#92;rangle = 0 &#92;quad &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t &#92;ge 0 ' title='&#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O &#92;rangle_{&#92;psi(t)} = &#92;langle &#92;hat{O} | H e^{t H} |&#92;psi_0 &#92;rangle = 0 &#92;quad &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t &#92;ge 0 ' class='latex' /></p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Ciff+%5Cquad+%5Clangle+%5Chat%7BO%7D%7C+H+e%5E%7Bt+H%7D+%3D+0+%5Cquad+%5Cforall+t+%5Cge+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;iff &#92;quad &#92;langle &#92;hat{O}| H e^{t H} = 0 &#92;quad &#92;forall t &#92;ge 0 ' title='&#92;iff &#92;quad &#92;langle &#92;hat{O}| H e^{t H} = 0 &#92;quad &#92;forall t &#92;ge 0 ' class='latex' /></p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Ciff+%5Cquad+%5Csum_k+%5Clangle+%5Chat%7BO%7D+%7C+%5Cphi_k+%5Crangle+E_k+e%5E%7Bt+E_k%7D+%5Clangle+%5Cphi_k%7C+%3D+0+%5Cquad+%5Cforall+t+%5Cge+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;iff &#92;quad &#92;sum_k &#92;langle &#92;hat{O} | &#92;phi_k &#92;rangle E_k e^{t E_k} &#92;langle &#92;phi_k| = 0 &#92;quad &#92;forall t &#92;ge 0 ' title='&#92;iff &#92;quad &#92;sum_k &#92;langle &#92;hat{O} | &#92;phi_k &#92;rangle E_k e^{t E_k} &#92;langle &#92;phi_k| = 0 &#92;quad &#92;forall t &#92;ge 0 ' class='latex' /></p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Ciff+%5Cquad+%5Clangle+%5Chat%7BO%7D+%7C+%5Cphi_k+%5Crangle+E_k+e%5E%7Bt+E_k%7D+%3D+0+%5Cquad+%5Cforall+t+%5Cge+0+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;iff &#92;quad &#92;langle &#92;hat{O} | &#92;phi_k &#92;rangle E_k e^{t E_k} = 0 &#92;quad &#92;forall t &#92;ge 0 ' title='&#92;iff &#92;quad &#92;langle &#92;hat{O} | &#92;phi_k &#92;rangle E_k e^{t E_k} = 0 &#92;quad &#92;forall t &#92;ge 0 ' class='latex' /></p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Ciff+%5Cquad+%7C%5Chat%7BO%7D+%5Crangle+%5Cin+%5Cmathrm%7BSpan%7D%5C%7B%7C%5Cphi_k+%5Crangle+%5C%2C+%3A+%5C%3B+E_k+%3D+0%5C%7D+%3D+%5Cker+%5C%3A+H%2C+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;iff &#92;quad |&#92;hat{O} &#92;rangle &#92;in &#92;mathrm{Span}&#92;{|&#92;phi_k &#92;rangle &#92;, : &#92;; E_k = 0&#92;} = &#92;ker &#92;: H, ' title='&#92;iff &#92;quad |&#92;hat{O} &#92;rangle &#92;in &#92;mathrm{Span}&#92;{|&#92;phi_k &#92;rangle &#92;, : &#92;; E_k = 0&#92;} = &#92;ker &#92;: H, ' class='latex' /></p>
<p>where the third equivalence is due to the vectors <img src='https://s0.wp.com/latex.php?latex=%7C%5Cphi_k+%5Crangle&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;phi_k &#92;rangle' title='|&#92;phi_k &#92;rangle' class='latex' /> being linearly independent. For any infinitesimal stochastic operator <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> the corresponding transition graph consists of <img src='https://s0.wp.com/latex.php?latex=m&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='m' title='m' class='latex' /> connected components iff we can reorder (permute) the states of the system such that <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> becomes block-diagonal with <img src='https://s0.wp.com/latex.php?latex=m&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='m' title='m' class='latex' /> blocks. Now it is easy to see that the kernel of <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is spanned by <img src='https://s0.wp.com/latex.php?latex=m&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='m' title='m' class='latex' /> eigenvectors, one for each block. Since <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is also symmetric, the elements of each such vector can be chosen to be ones within the block and zeros outside it. Consequently</p>
<p><img src='https://s0.wp.com/latex.php?latex=%7C%5Chat%7BO%7D+%5Crangle+%5Cin+%5Cker+%5C%3A+H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;hat{O} &#92;rangle &#92;in &#92;ker &#92;: H' title='|&#92;hat{O} &#92;rangle &#92;in &#92;ker &#92;: H' class='latex' /></p>
<p>implies that we can choose the basis of eigenvectors of <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> to be the vectors <img src='https://s0.wp.com/latex.php?latex=%7C%5Cphi_k+%5Crangle%2C&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;phi_k &#92;rangle,' title='|&#92;phi_k &#92;rangle,' class='latex' /> which implies</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5BO%2C+H%5D+%3D+0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O, H] = 0' title='[O, H] = 0' class='latex' /></p>
<p>Alternatively,</p>
<p><img src='https://s0.wp.com/latex.php?latex=%7C%5Chat%7BO%7D+%5Crangle+%5Cin+%5Cker+%5C%2C+H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;hat{O} &#92;rangle &#92;in &#92;ker &#92;, H' title='|&#92;hat{O} &#92;rangle &#92;in &#92;ker &#92;, H' class='latex' /> </p>
<p>implies that</p>
<p><img src='https://s0.wp.com/latex.php?latex=%7C%5Chat%7BO%5E2%7D+%5Crangle+%5Cin+%5Cker+%5C%3A+H+%5C%3B+%5Ciff+%5C%3B+%5Ccdots+%5C%3B+%5Ciff+%5C%3B+%5Cdisplaystyle%7B%5Cfrac%7Bd%7D%7Bd+t%7D%7D+%5Clangle+O%5E2+%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+0+%5C%3B+%5Cforall+%5C%3A+%7C%5Cpsi_0+%5Crangle%2C+%5Cforall+t+%5Cge+0%2C+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='|&#92;hat{O^2} &#92;rangle &#92;in &#92;ker &#92;: H &#92;; &#92;iff &#92;; &#92;cdots &#92;; &#92;iff &#92;; &#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O^2 &#92;rangle_{&#92;psi(t)} = 0 &#92;; &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t &#92;ge 0, ' title='|&#92;hat{O^2} &#92;rangle &#92;in &#92;ker &#92;: H &#92;; &#92;iff &#92;; &#92;cdots &#92;; &#92;iff &#92;; &#92;displaystyle{&#92;frac{d}{d t}} &#92;langle O^2 &#92;rangle_{&#92;psi(t)} = 0 &#92;; &#92;forall &#92;: |&#92;psi_0 &#92;rangle, &#92;forall t &#92;ge 0, ' class='latex' /></p>
<p>where we have used the above sequence of equivalences backwards. Now, using John and Brendan&#8217;s original proof, we can obtain <img src='https://s0.wp.com/latex.php?latex=%5BO%2C+H%5D+%3D+0.&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='[O, H] = 0.' title='[O, H] = 0.' class='latex' /> &nbsp;  &#9608;</p>
<p>In summary, by restricting ourselves to the intersection of quantum and stochastic generators, we have found a version of Noether&#8217;s theorem for stochastic mechanics that looks formally just like the quantum version!  However, this simplification comes at a cost.  We find that the only observables <img src='https://s0.wp.com/latex.php?latex=O&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O' title='O' class='latex' /> whose expected value remains constant with time are those of the very restricted type described above, where the observable has the same value in every state in a connected component.</p>
<h3>Puzzles</h3>
<p>Suppose we have a graph whose <a href="http://en.wikipedia.org/wiki/Laplacian_matrix">graph Laplacian matrix</a> <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> generates a Markov semigroup as follows:</p>
<p><img src='https://s0.wp.com/latex.php?latex=U%28t%29+%3D+e%5E%7Bt+H%7D+&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='U(t) = e^{t H} ' title='U(t) = e^{t H} ' class='latex' /> </p>
<p><b>Puzzle 1.</b>  Suppose that also <img src='https://s0.wp.com/latex.php?latex=H+%3D+H%5E%5Ctop%2C&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H = H^&#92;top,' title='H = H^&#92;top,' class='latex' /> so that <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is a Dirichlet operator and hence <img src='https://s0.wp.com/latex.php?latex=i+H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='i H' title='i H' class='latex' /> generates a 1-parameter unitary group.  Show that the <a href="http://en.wikipedia.org/wiki/Directed_graph#Indegree_and_outdegree">indegree and outdegree</a> of any node of our graph must be equal.  Graphs with this property are called <b>balanced</b>.</p>
<p><b>Puzzle 2.</b>  Suppose that <img src='https://s0.wp.com/latex.php?latex=U%28t%29+%3D+e%5E%7Bt+H%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='U(t) = e^{t H}' title='U(t) = e^{t H}' class='latex' /> is <a href="http://en.wikipedia.org/wiki/Doubly_stochastic_matrix">doubly stochastic</a> Markov semigroup, meaning that for all <img src='https://s0.wp.com/latex.php?latex=t+%5Cge+0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='t &#92;ge 0' title='t &#92;ge 0' class='latex' /> each row and each column of <img src='https://s0.wp.com/latex.php?latex=U%28t%29&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='U(t)' title='U(t)' class='latex' /> sums to 1:   </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B+%5Csum_i+U%28t%29_%7Bi+j%7D+%3D+%5Csum_j+U%28t%29_%7Bi+j%7D+%3D+1+%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{ &#92;sum_i U(t)_{i j} = &#92;sum_j U(t)_{i j} = 1 }' title='&#92;displaystyle{ &#92;sum_i U(t)_{i j} = &#92;sum_j U(t)_{i j} = 1 }' class='latex' /></p>
<p>and all the matrix entries are nonnegative.  Show that the Hamiltonian <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> obeys</p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B%5Csum_i+H_%7Bi+j%7D+%3D+%5Csum_j+H_%7Bi+j%7D+%3D+0+%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{&#92;sum_i H_{i j} = &#92;sum_j H_{i j} = 0 }' title='&#92;displaystyle{&#92;sum_i H_{i j} = &#92;sum_j H_{i j} = 0 }' class='latex' /></p>
<p>and all the off-diagonal entries of <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> are nonnegative.  Show the converse is also true.</p>
<p><b>Puzzle 3.</b>  Prove that any doubly stochastic Markov semigroup <img src='https://s0.wp.com/latex.php?latex=U%28t%29&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='U(t)' title='U(t)' class='latex' /> is of the form <img src='https://s0.wp.com/latex.php?latex=e%5E%7Bt+H%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='e^{t H}' title='e^{t H}' class='latex' /> where <img src='https://s0.wp.com/latex.php?latex=H&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='H' title='H' class='latex' /> is the graph Laplacian of a balanced graph.</p>
<p><b>Puzzle 4.</b> Let <img src='https://s0.wp.com/latex.php?latex=O%28t%29&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='O(t)' title='O(t)' class='latex' /> be a possibly time-dependent observable, and write <img src='https://s0.wp.com/latex.php?latex=%5Clangle+O%28t%29+%5Crangle_%7B%5Cpsi%28t%29%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;langle O(t) &#92;rangle_{&#92;psi(t)}' title='&#92;langle O(t) &#92;rangle_{&#92;psi(t)}' class='latex' /> for its expected value with respect to some initial state <img src='https://s0.wp.com/latex.php?latex=%5Cpsi_0&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;psi_0' title='&#92;psi_0' class='latex' /> evolving according to the master equation.  Show that </p>
<p><img src='https://s0.wp.com/latex.php?latex=%5Cdisplaystyle%7B+%5Cfrac%7Bd%7D%7Bd+t%7D%5Clangle+O%28t%29%5Crangle_%7B%5Cpsi%28t%29%7D+%3D+%5Cleft%5Clangle+%5BO%28t%29%2C+H%5D+%5Cright%5Crangle_%7B%5Cpsi%28t%29%7D+%2B+%5Cleft%5Clangle+%5Cfrac%7B%5Cpartial+O%28t%29%7D%7B%5Cpartial+t%7D%5Cright%5Crangle_%7B%5Cpsi%28t%29%7D+%7D&#038;bg=ffffff&#038;fg=000&#038;s=0' alt='&#92;displaystyle{ &#92;frac{d}{d t}&#92;langle O(t)&#92;rangle_{&#92;psi(t)} = &#92;left&#92;langle [O(t), H] &#92;right&#92;rangle_{&#92;psi(t)} + &#92;left&#92;langle &#92;frac{&#92;partial O(t)}{&#92;partial t}&#92;right&#92;rangle_{&#92;psi(t)} }' title='&#92;displaystyle{ &#92;frac{d}{d t}&#92;langle O(t)&#92;rangle_{&#92;psi(t)} = &#92;left&#92;langle [O(t), H] &#92;right&#92;rangle_{&#92;psi(t)} + &#92;left&#92;langle &#92;frac{&#92;partial O(t)}{&#92;partial t}&#92;right&#92;rangle_{&#92;psi(t)} }' class='latex' /></p>
<p>This is a stochastic version of the <a href="http://en.wikipedia.org/wiki/Ehrenfest_theorem">Ehrenfest theorem</a>. </p>
]]></html><thumbnail_url><![CDATA[https://i0.wp.com/upload.wikimedia.org/wikipedia/commons/e/e5/Noether.jpg?fit=440%2C330]]></thumbnail_url><thumbnail_height><![CDATA[330]]></thumbnail_height><thumbnail_width><![CDATA[216]]></thumbnail_width></oembed>