\documentclass[12pt,titlepage]{article}
\usepackage{amsmath}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsthm}
\usepackage{mathtools}
\usepackage{graphicx}
\usepackage{color}
\usepackage{ucs}
\usepackage[utf8x]{inputenc}
\usepackage{xparse}
\usepackage{hyperref}
%----Macros----------
%
% Unresolved issues:
%
% \righttoleftarrow
% \lefttorightarrow
%
% \color{} with HTML colorspec
% \bgcolor
% \array with options (without options, it's equivalent to the matrix environment)
% Of the standard HTML named colors, white, black, red, green, blue and yellow
% are predefined in the color package. Here are the rest.
\definecolor{aqua}{rgb}{0, 1.0, 1.0}
\definecolor{fuschia}{rgb}{1.0, 0, 1.0}
\definecolor{gray}{rgb}{0.502, 0.502, 0.502}
\definecolor{lime}{rgb}{0, 1.0, 0}
\definecolor{maroon}{rgb}{0.502, 0, 0}
\definecolor{navy}{rgb}{0, 0, 0.502}
\definecolor{olive}{rgb}{0.502, 0.502, 0}
\definecolor{purple}{rgb}{0.502, 0, 0.502}
\definecolor{silver}{rgb}{0.753, 0.753, 0.753}
\definecolor{teal}{rgb}{0, 0.502, 0.502}
% Because of conflicts, \space and \mathop are converted to
% \itexspace and \operatorname during preprocessing.
% itex: \space{ht}{dp}{wd}
%
% Height and baseline depth measurements are in units of tenths of an ex while
% the width is measured in tenths of an em.
\makeatletter
\newdimen\itex@wd%
\newdimen\itex@dp%
\newdimen\itex@thd%
\def\itexspace#1#2#3{\itex@wd=#3em%
\itex@wd=0.1\itex@wd%
\itex@dp=#2ex%
\itex@dp=0.1\itex@dp%
\itex@thd=#1ex%
\itex@thd=0.1\itex@thd%
\advance\itex@thd\the\itex@dp%
\makebox[\the\itex@wd]{\rule[-\the\itex@dp]{0cm}{\the\itex@thd}}}
\makeatother
% \tensor and \multiscript
\makeatletter
\newif\if@sup
\newtoks\@sups
\def\append@sup#1{\edef\act{\noexpand\@sups={\the\@sups #1}}\act}%
\def\reset@sup{\@supfalse\@sups={}}%
\def\mk@scripts#1#2{\if #2/ \if@sup ^{\the\@sups}\fi \else%
\ifx #1_ \if@sup ^{\the\@sups}\reset@sup \fi {}_{#2}%
\else \append@sup#2 \@suptrue \fi%
\expandafter\mk@scripts\fi}
\def\tensor#1#2{\reset@sup#1\mk@scripts#2_/}
\def\multiscripts#1#2#3{\reset@sup{}\mk@scripts#1_/#2%
\reset@sup\mk@scripts#3_/}
\makeatother
% \slash
\makeatletter
\newbox\slashbox \setbox\slashbox=\hbox{$/$}
\def\itex@pslash#1{\setbox\@tempboxa=\hbox{$#1$}
\@tempdima=0.5\wd\slashbox \advance\@tempdima 0.5\wd\@tempboxa
\copy\slashbox \kern-\@tempdima \box\@tempboxa}
\def\slash{\protect\itex@pslash}
\makeatother
% math-mode versions of \rlap, etc
% from Alexander Perlis, "A complement to \smash, \llap, and lap"
% http://math.arizona.edu/~aprl/publications/mathclap/
\def\clap#1{\hbox to 0pt{\hss#1\hss}}
\def\mathllap{\mathpalette\mathllapinternal}
\def\mathrlap{\mathpalette\mathrlapinternal}
\def\mathclap{\mathpalette\mathclapinternal}
\def\mathllapinternal#1#2{\llap{$\mathsurround=0pt#1{#2}$}}
\def\mathrlapinternal#1#2{\rlap{$\mathsurround=0pt#1{#2}$}}
\def\mathclapinternal#1#2{\clap{$\mathsurround=0pt#1{#2}$}}
% Renames \sqrt as \oldsqrt and redefine root to result in \sqrt[#1]{#2}
\let\oldroot\root
\def\root#1#2{\oldroot #1 \of{#2}}
\renewcommand{\sqrt}[2][]{\oldroot #1 \of{#2}}
% Manually declare the txfonts symbolsC font
\DeclareSymbolFont{symbolsC}{U}{txsyc}{m}{n}
\SetSymbolFont{symbolsC}{bold}{U}{txsyc}{bx}{n}
\DeclareFontSubstitution{U}{txsyc}{m}{n}
% Manually declare the stmaryrd font
\DeclareSymbolFont{stmry}{U}{stmry}{m}{n}
\SetSymbolFont{stmry}{bold}{U}{stmry}{b}{n}
% Manually declare the MnSymbolE font
\DeclareFontFamily{OMX}{MnSymbolE}{}
\DeclareSymbolFont{mnomx}{OMX}{MnSymbolE}{m}{n}
\SetSymbolFont{mnomx}{bold}{OMX}{MnSymbolE}{b}{n}
\DeclareFontShape{OMX}{MnSymbolE}{m}{n}{
<-6> MnSymbolE5
<6-7> MnSymbolE6
<7-8> MnSymbolE7
<8-9> MnSymbolE8
<9-10> MnSymbolE9
<10-12> MnSymbolE10
<12-> MnSymbolE12}{}
% Declare specific arrows from txfonts without loading the full package
\makeatletter
\def\re@DeclareMathSymbol#1#2#3#4{%
\let#1=\undefined
\DeclareMathSymbol{#1}{#2}{#3}{#4}}
\re@DeclareMathSymbol{\neArrow}{\mathrel}{symbolsC}{116}
\re@DeclareMathSymbol{\neArr}{\mathrel}{symbolsC}{116}
\re@DeclareMathSymbol{\seArrow}{\mathrel}{symbolsC}{117}
\re@DeclareMathSymbol{\seArr}{\mathrel}{symbolsC}{117}
\re@DeclareMathSymbol{\nwArrow}{\mathrel}{symbolsC}{118}
\re@DeclareMathSymbol{\nwArr}{\mathrel}{symbolsC}{118}
\re@DeclareMathSymbol{\swArrow}{\mathrel}{symbolsC}{119}
\re@DeclareMathSymbol{\swArr}{\mathrel}{symbolsC}{119}
\re@DeclareMathSymbol{\nequiv}{\mathrel}{symbolsC}{46}
\re@DeclareMathSymbol{\Perp}{\mathrel}{symbolsC}{121}
\re@DeclareMathSymbol{\Vbar}{\mathrel}{symbolsC}{121}
\re@DeclareMathSymbol{\sslash}{\mathrel}{stmry}{12}
\re@DeclareMathSymbol{\bigsqcap}{\mathop}{stmry}{"64}
\re@DeclareMathSymbol{\biginterleave}{\mathop}{stmry}{"6}
\re@DeclareMathSymbol{\invamp}{\mathrel}{symbolsC}{77}
\re@DeclareMathSymbol{\parr}{\mathrel}{symbolsC}{77}
\makeatother
% \llangle, \rrangle, \lmoustache and \rmoustache from MnSymbolE
\makeatletter
\def\Decl@Mn@Delim#1#2#3#4{%
\if\relax\noexpand#1%
\let#1\undefined
\fi
\DeclareMathDelimiter{#1}{#2}{#3}{#4}{#3}{#4}}
\def\Decl@Mn@Open#1#2#3{\Decl@Mn@Delim{#1}{\mathopen}{#2}{#3}}
\def\Decl@Mn@Close#1#2#3{\Decl@Mn@Delim{#1}{\mathclose}{#2}{#3}}
\Decl@Mn@Open{\llangle}{mnomx}{'164}
\Decl@Mn@Close{\rrangle}{mnomx}{'171}
\Decl@Mn@Open{\lmoustache}{mnomx}{'245}
\Decl@Mn@Close{\rmoustache}{mnomx}{'244}
\makeatother
% Widecheck
\makeatletter
\DeclareRobustCommand\widecheck[1]{{\mathpalette\@widecheck{#1}}}
\def\@widecheck#1#2{%
\setbox\z@\hbox{\m@th$#1#2$}%
\setbox\tw@\hbox{\m@th$#1%
\widehat{%
\vrule\@width\z@\@height\ht\z@
\vrule\@height\z@\@width\wd\z@}$}%
\dp\tw@-\ht\z@
\@tempdima\ht\z@ \advance\@tempdima2\ht\tw@ \divide\@tempdima\thr@@
\setbox\tw@\hbox{%
\raise\@tempdima\hbox{\scalebox{1}[-1]{\lower\@tempdima\box
\tw@}}}%
{\ooalign{\box\tw@ \cr \box\z@}}}
\makeatother
% \mathraisebox{voffset}[height][depth]{something}
\makeatletter
\NewDocumentCommand\mathraisebox{moom}{%
\IfNoValueTF{#2}{\def\@temp##1##2{\raisebox{#1}{$\m@th##1##2$}}}{%
\IfNoValueTF{#3}{\def\@temp##1##2{\raisebox{#1}[#2]{$\m@th##1##2$}}%
}{\def\@temp##1##2{\raisebox{#1}[#2][#3]{$\m@th##1##2$}}}}%
\mathpalette\@temp{#4}}
\makeatletter
% udots (taken from yhmath)
\makeatletter
\def\udots{\mathinner{\mkern2mu\raise\p@\hbox{.}
\mkern2mu\raise4\p@\hbox{.}\mkern1mu
\raise7\p@\vbox{\kern7\p@\hbox{.}}\mkern1mu}}
\makeatother
%% Fix array
\newcommand{\itexarray}[1]{\begin{matrix}#1\end{matrix}}
%% \itexnum is a noop
\newcommand{\itexnum}[1]{#1}
%% Renaming existing commands
\newcommand{\underoverset}[3]{\underset{#1}{\overset{#2}{#3}}}
\newcommand{\widevec}{\overrightarrow}
\newcommand{\darr}{\downarrow}
\newcommand{\nearr}{\nearrow}
\newcommand{\nwarr}{\nwarrow}
\newcommand{\searr}{\searrow}
\newcommand{\swarr}{\swarrow}
\newcommand{\curvearrowbotright}{\curvearrowright}
\newcommand{\uparr}{\uparrow}
\newcommand{\downuparrow}{\updownarrow}
\newcommand{\duparr}{\updownarrow}
\newcommand{\updarr}{\updownarrow}
\newcommand{\gt}{>}
\newcommand{\lt}{<}
\newcommand{\map}{\mapsto}
\newcommand{\embedsin}{\hookrightarrow}
\newcommand{\Alpha}{A}
\newcommand{\Beta}{B}
\newcommand{\Zeta}{Z}
\newcommand{\Eta}{H}
\newcommand{\Iota}{I}
\newcommand{\Kappa}{K}
\newcommand{\Mu}{M}
\newcommand{\Nu}{N}
\newcommand{\Rho}{P}
\newcommand{\Tau}{T}
\newcommand{\Upsi}{\Upsilon}
\newcommand{\omicron}{o}
\newcommand{\lang}{\langle}
\newcommand{\rang}{\rangle}
\newcommand{\Union}{\bigcup}
\newcommand{\Intersection}{\bigcap}
\newcommand{\Oplus}{\bigoplus}
\newcommand{\Otimes}{\bigotimes}
\newcommand{\Wedge}{\bigwedge}
\newcommand{\Vee}{\bigvee}
\newcommand{\coproduct}{\coprod}
\newcommand{\product}{\prod}
\newcommand{\closure}{\overline}
\newcommand{\integral}{\int}
\newcommand{\doubleintegral}{\iint}
\newcommand{\tripleintegral}{\iiint}
\newcommand{\quadrupleintegral}{\iiiint}
\newcommand{\conint}{\oint}
\newcommand{\contourintegral}{\oint}
\newcommand{\infinity}{\infty}
\newcommand{\bottom}{\bot}
\newcommand{\minusb}{\boxminus}
\newcommand{\plusb}{\boxplus}
\newcommand{\timesb}{\boxtimes}
\newcommand{\intersection}{\cap}
\newcommand{\union}{\cup}
\newcommand{\Del}{\nabla}
\newcommand{\odash}{\circleddash}
\newcommand{\negspace}{\!}
\newcommand{\widebar}{\overline}
\newcommand{\textsize}{\normalsize}
\renewcommand{\scriptsize}{\scriptstyle}
\newcommand{\scriptscriptsize}{\scriptscriptstyle}
\newcommand{\mathfr}{\mathfrak}
\newcommand{\statusline}[2]{#2}
\newcommand{\tooltip}[2]{#2}
\newcommand{\toggle}[2]{#2}
% Theorem Environments
\theoremstyle{plain}
\newtheorem{theorem}{Theorem}
\newtheorem{lemma}{Lemma}
\newtheorem{prop}{Proposition}
\newtheorem{cor}{Corollary}
\newtheorem*{utheorem}{Theorem}
\newtheorem*{ulemma}{Lemma}
\newtheorem*{uprop}{Proposition}
\newtheorem*{ucor}{Corollary}
\theoremstyle{definition}
\newtheorem{defn}{Definition}
\newtheorem{example}{Example}
\newtheorem*{udefn}{Definition}
\newtheorem*{uexample}{Example}
\theoremstyle{remark}
\newtheorem{remark}{Remark}
\newtheorem{note}{Note}
\newtheorem*{uremark}{Remark}
\newtheorem*{unote}{Note}
%-------------------------------------------------------------------
\begin{document}
%-------------------------------------------------------------------
\section*{Blog - quantum network theory (part 2)}
This page is a [[Blog articles in progress|blog article in progress]], written by [[Tomi Johnson]]. Part 1 is found [[Blog - quantum network theory (part 1)|here]]. To see discussion of the pair of articles while they are being written, visit the \href{http://forum.azimuthproject.org/discussion/1228/blog-quantum-complex-networks/}{Azimuth Forum}.
guest post by
I make a couple of references to part 1. Currently I don't add any hyperlinks to the text, as I wouldn't know where to link to. Let me know if I should do anything different. - Tomi
Last time I told you how a stochastic walk (termed the \textbf{uniform escape walk}) could be used to analyze a network. In particular, Google uses it to rank nodes. For the case of an undirected network, the steady state of the stochastic walk tells us the degrees of the nodes.
Here, I'm going to prove this to you. I'm also going to exploit the connection between this stochastic walk and a corresponding quantum walk, also introduced last time. In particular, I'll connected the properties of this quantum walk to the degrees of a network by exploiting its relationship with the stochastic walk.
This is pretty useful considering how tricky these quantum walks can be to study. As the parts of the world that we model using quantum mechanics get bigger and have more complicated structures, like biological network, we need all the help in understanding quantum walks that we can get. Then I better start!
Previously, we introduced a stochastic and quantum walk, generated by the similar matrices $S$ and $Q$ respectively. The first describes the uniform escape stochastic walk who's steady state tells us so much about the degrees of a network. The second is the related quantum walk we're going to solve in this article.
Both walks follow from a graph and all we learnt last time is summarized by the following diagram
where
$\bullet$ $G$ is a simple graph that specifies
$\bullet$ $A$ the adjacency matrix (generator of a quantum walk) with elements $A_{i j}$ equal to unity if nodes $i$ and $j$ are connected and zero otherwise ($A_{i i} = 0$), which subtracted from
$\bullet$ $D$ the diagonal matrix of degrees $D_{i i} = \sum_j A_{i j}$ gives
$\bullet$ $L = D - A$ the symmetric Laplacian (generator of stochastic and quantum walks), which when normalized by $D$ returns both
$\bullet$ $S = L D^{-1}$ the generator of the uniform escape stochastic walk and
$\bullet$ $Q = D^{-1/2} L D^{-1/2}$ the quantum walk generator to which it is similar!
Now we remember where we are, next I’ll talk you through the mathematics of the uniform escape stochastic walk $S$ and how it connects to the degrees of the nodes in the long-time limit. Then I’ll show you how this helps us solve aspects of the quantum walk generated by $Q$.
The uniform escape stochastic walk generated by $S$ is popular because it has a really useful .
This standout should appear OK on the blog, sorry for its current appearance. - Tomi
To find this unique stationary state, consider the Laplacian $L$, which is both infinitesimal stochastic and symmetric. Among other things, this means the rows of $L$ sum to zero
\begin{displaymath}
\displaystyle{ \sum_j L_{i j} = 0 }
\end{displaymath}
which means the all ones vector $\mathbf{1}$ is an eigenvector of $L$ with zero eigenvalue
\begin{displaymath}
L \mathbf{1} = 0
\end{displaymath}
Inserting the identity $I = D^{-1} D$ into this equation we then find $D \mathbf{1}$ is a zero eigenvector of $S$
\begin{displaymath}
L \mathbf{1} = ( L D^{-1} ) ( D \mathbf{1} ) = S ( D \mathbf{1} ) = 0
\end{displaymath}
Therefore we just need to normalize this to get the long-time stationary state of the walk
\begin{displaymath}
\pi = \frac{D \mathbf{1}}{\sum_i D_{i i}}
\end{displaymath}
Each element $\langle i , \pi \rangle$ of this state, the long-time probability of finding a walker at node $i$, is proportional to the degree $D_{i i}$ of node $i$!
We can check this holds for the following graph
We find that $\pi$ is
\begin{displaymath}
\left( \begin{matrix} 1/6 \\ 1/6 \\ 1/4 \\ 1/4 \\ 1/6 \end{matrix} \right)
\end{displaymath}
which implies long-time probability $1/6$ for nodes $1$, $2$ and $5$, and $1/4$ for nodes $3$ and $4$. Comparing this to the original graph, this exactly reflects the arrangement of degrees, as we knew it must. Math works!
Next up is the quantum walk generated by $Q$. Not a lot is known about quantum walks on networks of arbitrary geometry, but below we’ll see some analytical results are obtained by exploiting the similarity of $S$ and $Q$.
Where to start? Well, let's start from the bottom, what quantum physicists call the . In contrast to stochastic walks, for a quantum walk every eigenvector $\phi_i^k$ of $Q$ is a stationary state of the quantum walk (in puzzle 5, at the bottom of this page, I ask you to prove this). The stationary state $\phi_0$ is of particular interest physically and mathematically. Physically, since eigenvectors of the $Q$ correspond to states of well-defined energy equal to the associated eigenvalue, $\phi_0$ is the state of lowest energy $\epsilon_0 = 0$, hence the name ground state (in puzzle 3 we ask you to prove that all eigenvalues of $Q$ are non-negative, so zero really does correspond to the ground state).
Mathematically, the relationship between eigenvectors implied by the similarity of $S$ and $Q$ means $\phi_0 \propto D^{-1/2} \pi \propto D^{1/2} \mathbf{1}$. So, if in the ground state, the probability of being measured to be at node $i$ is $| \langle i , \phi_0 \rangle |^2 \propto | \langle i , D^{1/2} \rangle |^2 = D_{i i}$. Amazingly, this probability is proportional to the degree and so is exactly the same as $\langle i , \pi \rangle$, the probability in the stationary state $\pi$ of the stochastic walk.
A zero energy quantum walk $Q$ leads to exactly the same distribution of the walker over the nodes as in the long-time limit of the uniform escape stochastic walk $S$.
The normally classical notion of degree distribution plays a role in quantum walks!
This is already pretty exciting. But what else can we say? If you are someone who feels faint at the sight of quantum mechanics, well done for getting this far, but be prepared for what's coming next.
What if the walker starts in some other initial state? Is there some quantum walk analogue of the unique long-time state of a stochastic walk?
In fact, the quantum walk in general does not converge to a stationary state. But there is a probability distribution that can be thought to characterize the quantum walk in the same way as the long-time state characterizes the stochastic walk. It's the $P$.
So we have a distribution that characterizes a quantum walk for a general initial state but it's a complicated beast. What can we say about it?
Our best hope of understanding the long-time average probability is through the term $| \langle i, \Phi_0 \psi (0) \rangle |^2$ associated with the zero energy eigenspace, since we know everything about this space.
For example, we know the zero energy eigenspace is one-dimensional and spanned by the eigenvector $\phi_0$. This means that the projector is just the usual outer product
\begin{displaymath}
\Phi_0 = | \phi_0 \rangle \langle \phi_0 | = \phi_0 \phi_0^\dagger
\end{displaymath}
where we have normalized $\phi_0$ according to the inner product $\langle \phi_0, \phi_0\rangle = 1$. (If you're wondering why I'm using all these angled brackets, well, it's a type of notation named after that is adored by quantum physicists.)
The zero eigenspace contribution to the long-time average probability then breaks nicely into two
\begin{displaymath}
| \langle i, \Phi_0 \psi (0) \rangle |^2 = | \langle i, \phi_0\rangle \langle \phi_0, \psi (0) \rangle |^2 = | \langle i, \phi_0\rangle |^2 | \langle \phi_0 , \psi (0) \rangle |^2 = \langle i , \pi \rangle | \langle \phi_0 , \psi (0) \rangle |^2
\end{displaymath}
This is just the product of two probabilities!
First the probability $\langle i , \pi \rangle$ for a quantum state in the zero energy eigenspace to be at node $i$ (as we found above) and second the probability $| \langle \phi_0, \psi (0)\rangle |^2$ of being in this eigenspace to begin with (in quantum mechanics the probability of measuring the system to have an energy is the modulus squared of the projection of the state onto the associated eigenspace, which for the one dimensional zero energy eigenspace means just the inner product with the ground state.)
This is all we need to say something interesting about the long-time average probability for all states. We've basically shown that we can break the long-time probability vector $P$ into a sum of two normalized probability vectors
\begin{displaymath}
P = (1- \eta) \pi + \eta \Omega
\end{displaymath}
the first $\pi$ being the degree dependent stochastic stationary state associated with the zero energy eigenspace and the second associated with the higher energy eigenspaces $\Omega$, with
\begin{displaymath}
\displaystyle{ \langle i , \Omega \rangle = \frac{ \sum_{k\neq 0} | \langle i, \Phi_k \psi (0) \rangle |^2 }{ \eta} }
\end{displaymath}
The weight of each term is governed by the parameter
\begin{displaymath}
\eta = 1 - | \langle \phi_0, \psi (0)\rangle |^2
\end{displaymath}
which you could think of as the of the result. This is unity minus the probability of the walker being in the zero energy eigenspace, or equivalently the probability of the walker being outside the zero energy eigenspace.
So even though we don't know anything about $\Omega$ we know its importance is controlled by a parameter $\eta$ that governs how close the long-time average distribution $P$ of the quantum walk is to the corresponding stochastic stationary distribution $\pi$. What do we mean by close? Find out for yourself:
There are several occasions on which I wish to start a line with {\colorbox[rgb]{1.00,0.93,1.00}{\tt \char60b\char62}}. This causes an error. As a temporary work around, on each occasion I have added ``UGH - A BUG '' to the start of the line. - Tomi
UGH - A BUG Show, using a triangle inequality, that the between the two characteristic stochastic and quantum distributions $\{ \langle i , P \rangle \}_i$ and $\{ \langle i , \pi \rangle \}_i$ is upper-bounded by $2 \eta$.
Can we say anything physical about when the quantumness $\eta$ is big or small?
Because the eigenvalues of $Q$ have a physical interpretation in terms of energy, the answer is yes. The quantumness $\eta$ is the probability of being outside the zero energy state. Call the next lowest eigenvalue $\Delta = \min_{k \neq 0} \epsilon_k$ the energy gap. If the quantum walk is not in the zero energy eigenspace then it must be in an eigenspace of energy greater or equal to $\Delta$. Therefore the expected energy $E$ of the quantum walker must bound the quantumness $E \ge \eta \Delta$.
This tells us that a quantum walk with a low energy is similar to a stochastic walk in the long-time limit (we already knew this exactly in the zero energy limit).
So little is known about quantum walks on networks of arbitrary geometry, that we were very pleased to find this result. It says there is a special case in which the walk is characterized by the degree distribution of the network and gives us a clear physical parameter that bounds how far the walk is from this special case.
Also, in finding it we learnt that the difficulties of the initial state dependence, enhanced by the lack of convergence to a stationary state, could be overcome for a quantum walk, and that the relationships between quantum and stochastic walks extend beyond those with shared generators.
That’s all for the latest bit of idea sharing at the interface between stochastic and quantum systems.
I hope I’ve piqued your interest about quantum walks. There’s so much still left to work out about this topic, and your help is needed!
Other questions we have include: What holds analytically about the form of the quantum correction? Numerically it is known that the so-called quantum correction $\Omega$ tends to enhance the probability of being found on nodes of low degree compared to $\pi$, can someone explain why? What happens if a small amount of stochastic noise is added to a quantum walk? Or a lot of noise?
It’s difficult to know who is best placed to answer these questions: experts in quantum physics, graph theory, complex networks or stochastic processes? I suspect it’ll take a bit of help from everyone.
A couple of textbooks with comprehensive sections on non-negative matrices and continuous-time stochastic processes are:
$\bullet$ Peter Lancaster and Miron Tismenetsky, , 2nd Ed., Academic Press, San Diego, 1985.
$\bullet$ James R. Norris, , Cambridge University Press, Cambridge, 1997.
There is, of course, the book that arose from the Azimuth , which considers several relationships between quantum and stochastic processes on networks:
$\bullet$ John Baez and Jacob Biamonte, , 2012.
Another couple on complex networks are:
$\bullet$ Mark Newman, , Oxford University Press, Oxford, 2010.
$\bullet$ Ernesto Estrada, , Oxford University Press, Oxford, 2011. Note that the is available free online.
There are plenty more useful references in our article on this topic:
$\bullet$ Mauro Faccin, Tomi Johnson, Jacob Biamonte, Sabre Kais and Piotr Migda, .
Sadly I didn't have space to show proofs of all the theorems I used. So here are a few puzzles that guide you to doing the proofs for yourself:
UGH - A BUG (For the hard core) Prove there is always a unique positive eigenvector for a stochastic walk generated by $S$. You’ll need the assumption that the graph $G$ is connected. It's not simple, and you’ll probably need help from a book, perhaps one of those above by Lancaster and Tismenetsky, and Norris.
UGH - A BUG Show that the eigenvalues of $S$ (and therefore $Q$) are non-negative i.e. $\epsilon_k \ge 0$. A good way to start this proof is to apply the to the non-negative matrix $M = - S + I \max_i S_{i i}$. This means $M$ has a positive eigenvalue $r$ equal to its spectral radius ($r = \max_k | \lambda_k |$ where $\{ \lambda_k \}_k$ are the eigenvalues of $M$), and the associated eigenvector $v$ is positive. Since $S = - M + I \max_i S_{i i}$, it follows that $S$ shares the eigenvectors of $M$ and the associated eigenvalues are related by inverted translation $\epsilon_k = - \lambda_k + \max_i S_{i i}$.
UGH - A BUG Prove that regardless of the initial state $\psi(0)$, the zero eigenvector $\pi$ is obtained in the long-time limit $\lim_{t \rightarrow \infty} \psi(t) = \pi$ of the walk generated by $S$. This breaks down into two parts:
UGH - A BUG Using the approach from puzzle 5, to show that $S v = \epsilon_v v$, the positivity of $v$ and the infinitesimal stochastic property $\sum_i S_{i j} = 0$ imply that $\epsilon_v = \epsilon_0 = 0$ and thus $v = \pi$ is actually the unique zero eigenvector and stationary state of $S$ (its uniqueness follows from puzzle 4, you don’t need to re-prove it).
UGH - A BUG By inserting the decomposition $S = \sum_k \epsilon_k \Pi_k$ into $e^{-S t}$ and using the result of puzzle 5, complete the proof.
(Though I ask you to use the diagonalizability of $S$, the main results still hold if the generator is irreducible but not diagonalizable.)
Here are a couple of extra puzzles for those of you interested in quantum mechanics:
UGH - A BUG In quantum mechanics, probabilities are given by the moduli squared of amplitudes, so multiplying a state by a number of modulus unity has no physical effect. By inserting $Q= \sum_k \epsilon_k \Phi_k$ into the quantum evolution matrix $e^{-Q t}$ show that if $\psi(0) = \phi_i^k$ then $\psi(t) = e^{ - i \epsilon_k t} \psi(0)$, hence $\phi_i^k$ is a stationary state as probabilities don't change in time.
UGH - A BUG By expanding the initial state $\psi(0)$ in terms of the complete orthogonal basis vectors $\{ \phi_k^i \}_i$, show that for a quantum walk $\psi(t)$ never converges to a stationary state unless it began in one.
\end{document}