From 61802fe09ed222c69983fe8fef5ce107ba3aa223 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Fri, 5 Mar 2021 17:01:39 +0100
Subject: Initial commit.

---
 stokes.tex | 176 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 176 insertions(+)
 create mode 100644 stokes.tex

(limited to 'stokes.tex')

diff --git a/stokes.tex b/stokes.tex
new file mode 100644
index 0000000..13138f7
--- /dev/null
+++ b/stokes.tex
@@ -0,0 +1,176 @@
+\documentclass[aps,reprint,longbibliography,floatfix]{revtex4-2}
+
+\usepackage[utf8]{inputenc} % why not type "Stokes" with unicode?
+\usepackage[T1]{fontenc} % vector fonts plz
+\usepackage{newtxtext,newtxmath} % Times for PR
+\usepackage[
+  colorlinks=true,
+  urlcolor=purple,
+  citecolor=purple,
+  filecolor=purple,
+  linkcolor=purple
+]{hyperref} % ref and cite links with pretty colors
+\usepackage{amsmath, graphicx, xcolor} % standard packages
+
+\begin{document}
+
+\title{}
+
+\author{Jaron Kent-Dobias}
+\author{Jorge Kurchan}
+
+\affiliation{Laboratoire de Physique de l'Ecole Normale Supérieure, Paris, France}
+
+\date\today
+
+\begin{abstract}
+\end{abstract}
+
+\maketitle
+
+Consider a thermodynamic calculation involving the (real) $p$-spin model for a
+particular instantiation of the coupling tensor $J$
+\begin{equation} \label{eq:partition.function}
+  Z_J(\beta)=\int_{S^{N-1}}ds\,e^{-\beta H_J(s)}
+\end{equation}
+where $S^{N-1}$ is the $N-1$-sphere of radius $N$. Quantities of interest are
+usually related to the quenched free energy, produced by averaging over the
+$J$s the sample free energy $F_J$
+\begin{equation}
+  \overline F(\beta,\kappa)=-\beta^{-1}\int dJ\,p_\kappa(J)\log Z_J
+\end{equation}
+which can depend in general on the inverse temperature $\beta$ and on some
+parameter $\kappa$ which governs the distribution of $J$s. For most
+applications, $\beta$ is taken to be real and positive, and the distribution
+$p_\kappa$ is taken to be Gaussian or discrete on $\pm1$.
+
+We are interested in analytically continuing expressions like $\overline F$
+into the region of complex $\beta$ or distributions $p_\kappa$ involving
+complex $J$. The former has been considered extensively for the Gaussian
+$p$-spin in the past \cite{complex_energy}, while the latter is largely
+unexplored.
+
+When the argument of the exponential integrand in \eqref{eq:partition.function}
+acquires an imaginary component, various numeric and perturbative schemes for
+approximating its value can face immediate difficulties due to the emergence of
+a sign problem, resulting from rapid oscillations coinciding with saddles.
+
+The surface $M\subset\mathbb C^N$ defined by $z^2=N$ is an $N-1$ dimensional
+\emph{Stein manifold}, a type of complex manifold defined by the level set of a
+holomorphic function \cite{Forstneric_2017_Stein}. Suppose that $z:\mathbb
+C^{N-1}\to M$ is a holomorphic map. The Jacobian $J$ of the map is
+\begin{equation}
+  J_{i\alpha}=\frac{\partial z_i}{\partial u_\alpha}=\partial_\alpha z_i
+\end{equation}
+where Greek coefficients run from $1$ to $N-1$ and Latin coefficients from $1$
+to $N$. The hermitian metric is $g=J^\dagger J$. For any smooth function
+$\phi:M\to\mathbb R$, its gradient $\nabla\phi$ is a holomorphic vector field
+given by
+\begin{equation}
+  \nabla\phi=(\partial^*\phi)^\sharp=(\partial^*\phi)g^{-1}
+\end{equation}
+
+For
+coordinates $u\in\mathbb C^{N-1}$, dynamics consists of gradient descent on
+$\operatorname{Re}H$, or
+\begin{equation}
+  \dot u=-\nabla\operatorname{Re}H=-\tfrac12(\partial H)^\dagger g^{-1}
+\end{equation}
+These dynamics preserve $\operatorname{Im}H$ and in fact correspond to
+Hamiltonian dynamics, with the real and imaginary parts of the coordinates
+taking the role of conjugate variables. \cite{Morrow_2006_Complex}
+
+Working with a particular map is inconvenient, and we would like to develop a map-independent dynamics. Using the chain rule, one finds
+\begin{equation}
+  \begin{aligned}
+    \dot z_i
+    &=\dot u_\alpha\partial_\alpha z_i
+    =-\tfrac12(\partial_\beta H)^*g^{-1}_{\beta\alpha}\partial_\alpha z_i\\
+    &=-\tfrac12(\partial_j H)^*(\partial_\beta z_j)^*g^{-1}_{\beta\alpha}\partial_\alpha z_i
+    =-\tfrac12(\partial H)^\dagger(J^\dagger g^{-1}J)\\
+  \end{aligned}
+\end{equation}
+where $J$ is the Jacobian of the coordinate map and $g$ is the metric. In stereographic coordinates this can be worked out directly.
+Consider the coordinates $z_i=u_i$ for $1\leq i\leq N-1$ and
+\begin{equation}
+  z_N=\sqrt{N-u_\alpha u_\alpha}
+\end{equation}
+The Jacobian is
+\begin{equation}
+  J_{\alpha i}=\partial_\alpha z_i=\delta_{\alpha i}-\delta_{Ni}\frac{u_\alpha}{\sqrt{N-u_\beta u_\beta}}
+\end{equation}
+and the corresponding hermitian metric is
+\begin{equation}
+  g_{\alpha\beta}=J_{i\alpha}^*J_{i\beta}
+  =\delta_{\alpha\beta}+\frac{u_\alpha^*u_\beta}{|N-u_\gamma u_\gamma|}
+\end{equation}
+The metric can be inverted explicitly:
+\begin{equation}
+  g^{-1}_{\alpha\beta}
+  =\delta_{\alpha\beta}-\frac{u_\alpha^*u_\beta}{|N-u_\gamma u_\gamma|+|u|^2}.
+\end{equation}
+Putting these pieces together, we find
+\begin{equation}
+  (J^\dagger g^{-1}J)_{ij}
+  =\delta_{ij}-\frac{z_iz_j^*}{|z|^2}
+\end{equation}
+which is just the projector onto the constraint manifold.
+
+Therefore, a map-independent dynamics is given by
+\begin{equation}
+  \dot z
+  =-\frac12(\partial H)^\dagger\left(I-\frac{zz^\dagger}{|z|^2}\right)
+\end{equation}
+
+Stokes lines are trajectories that approach distinct critical points as time
+goes to $\pm\infty$. From the perspective of dynamics, these correspond to
+\emph{heteroclinic orbits}. What are the conditions under which Stokes lines
+appear? Because the dynamics conserves imaginary energy, two critical points
+must have the same imaginary energy if they are to be connected by a Stokes line.
+
+The level sets defined by $\operatorname{Im}H=c$ for $c\in\mathbb R$ are surfaces of
+$2(N-1)-1$ real dimensions. They must be simply connected, since gradient
+descent in $\operatorname{Re}H$ cannot pass the same point twice.
+
+\section{2-spin}
+
+\begin{equation}
+  H_0=\frac12z^TJz
+\end{equation}
+$J$ is generically diagonalizable by a complex orthogonal matrix $P$. With
+$z\mapsto Pz$, $J\mapsto D$ for diagonal $D$. Then $\partial_i H=d_iz_i$.
+Suppose that two critical points have the same imaginary energy; without loss
+of generality, assume these are the first and second components. Since the gradient is proportional to $z$, any components that are zero at some time will be zero at all times. The dynamics for the components of interest assuming all others are zero are
+\begin{equation}
+  \begin{aligned}
+    \dot z_1
+    &=-z_1^*\left(d_1^*-\frac{d_1^*z_1^*z_1+d_2^*z_2^*z_2}{|z_1|^2+|z_2|^2}\right) \\
+    &=-(d_1-d_2)^*z_1^*\frac{|z_2|^2}{|z_1|^2+|z_2|^2}
+  \end{aligned}
+\end{equation}
+and the same for $z_2$ with all indices swapped.  Since $\Delta=d_1-d_N$ is
+real, if $z_1$ begins real it remains real, with the same for $z_2$. Since the critical points are at real $z$, we make this restriction, and find
+\begin{equation}
+  \begin{aligned}
+    \frac d{dt}(z_1^2+z_2^2)=0\\
+    \frac d{dt}\frac{z_2}{z_1}=\Delta\frac{z_2}{z_1}
+  \end{aligned}
+\end{equation}
+Therefore $z_2/z_1=e^{\Delta t}$, with $z_1^2+z_2^2=N$ as necessary.  Depending on the sign of $\Delta$, $z$ flows
+from one critical point to the other over infinite time. This is a Stokes line,
+and establishes that any two critical points in the 2-spin model with the same
+imaginary energy will possess one.
+
+The critical points of the 2-spin model are all adjacent: no critical point is separated from another by the separatrix of a third. This means that when the imaginary energies of two critical points are brought to the same value, their surfaces of constant imaginary energy join.
+
+\section{p-spin}
+
+\section{(2 + 4)-spin}
+
+\begin{acknowledgments}
+  MIT mathematicians have been no help
+\end{acknowledgments}
+
+\bibliography{stokes}
+
+\end{document}
-- 
cgit v1.2.3-70-g09d2


From aec7b1fb31e146d4e192081cc600caa121b692c8 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Mon, 8 Mar 2021 17:16:30 +0100
Subject: More work.

---
 stokes.bib |  15 ++++++++
 stokes.tex | 116 ++++++++++++++++++++++++++++++++++++++-----------------------
 2 files changed, 88 insertions(+), 43 deletions(-)

(limited to 'stokes.tex')

diff --git a/stokes.bib b/stokes.bib
index a1148bc..fb25b3a 100644
--- a/stokes.bib
+++ b/stokes.bib
@@ -22,4 +22,19 @@
  isbn = {9780821840559}
 }
 
+@incollection{Witten_2011_Analytic,
+ author = {Witten, Edward},
+ title = {Analytic continuation of Chern-Simons theory},
+ publisher = {American Mathematical Society},
+ year = {2011},
+ month = {7},
+ volume = {50},
+ pages = {347--446},
+ url = {https://doi.org/10.1090%2Famsip%2F050%2F19},
+ doi = {10.1090/amsip/050/19},
+ booktitle = {Chern-Simons Gauge Theory: 20 Years After},
+ editor = {Andersen, Jørgen E. and Boden, Hans U. and Hahn, Atle and Himpel, Benjamin},
+ series = {AMS/IP Studies in Advanced Mathematics}
+}
+
 
diff --git a/stokes.tex b/stokes.tex
index 13138f7..1ac2298 100644
--- a/stokes.tex
+++ b/stokes.tex
@@ -28,6 +28,8 @@
 
 \maketitle
 
+\cite{Witten_2011_Analytic}
+
 Consider a thermodynamic calculation involving the (real) $p$-spin model for a
 particular instantiation of the coupling tensor $J$
 \begin{equation} \label{eq:partition.function}
@@ -55,68 +57,75 @@ acquires an imaginary component, various numeric and perturbative schemes for
 approximating its value can face immediate difficulties due to the emergence of
 a sign problem, resulting from rapid oscillations coinciding with saddles.
 
-The surface $M\subset\mathbb C^N$ defined by $z^2=N$ is an $N-1$ dimensional
+\section{Geometry}
+
+The surface $M\subset\mathbb C^N$ defined by $N=f(z)=z^2$ is an $N-1$ dimensional
 \emph{Stein manifold}, a type of complex manifold defined by the level set of a
-holomorphic function \cite{Forstneric_2017_Stein}. Suppose that $z:\mathbb
-C^{N-1}\to M$ is a holomorphic map. The Jacobian $J$ of the map is
+holomorphic function \cite{Forstneric_2017_Stein}. One can define a Hermitian
+metric $h$ on $M$ by taking the restriction of the standard metric of $\mathbb
+C^N$ to vectors tangent along $M$. For any smooth function $\phi:M\to\mathbb
+R$, its gradient $\nabla\phi$ is a holomorphic vector field given by
 \begin{equation}
-  J_{i\alpha}=\frac{\partial z_i}{\partial u_\alpha}=\partial_\alpha z_i
+  \nabla^\alpha\phi=h^{\bar\beta\alpha}\bar\partial_{\bar\beta}\phi
 \end{equation}
-where Greek coefficients run from $1$ to $N-1$ and Latin coefficients from $1$
-to $N$. The hermitian metric is $g=J^\dagger J$. For any smooth function
-$\phi:M\to\mathbb R$, its gradient $\nabla\phi$ is a holomorphic vector field
-given by
-\begin{equation}
-  \nabla\phi=(\partial^*\phi)^\sharp=(\partial^*\phi)g^{-1}
+Dynamics consists of gradient descent on $\operatorname{Re}H$, or
+\begin{equation} \label{eq:flow}
+  \dot u^\alpha=-\nabla^\alpha\operatorname{Re}H=-\tfrac12h^{\bar\beta\alpha}\bar\partial_{\bar\beta}\bar H
 \end{equation}
-
-For
-coordinates $u\in\mathbb C^{N-1}$, dynamics consists of gradient descent on
-$\operatorname{Re}H$, or
+Gradient descent on $\operatorname{Re}H$ is equivalent to Hamiltonian dynamics
+with the Hamiltonian $\operatorname{Im}H$. This is because $(M, h)$ is Kähler
+and therefore admits a symplectic structure, but that the flow conserves
+$\operatorname{Im}H$ can be seen from the Cauchy--Riemann equations and
+\eqref{eq:flow}:
 \begin{equation}
-  \dot u=-\nabla\operatorname{Re}H=-\tfrac12(\partial H)^\dagger g^{-1}
+  \begin{aligned}
+    \frac d{dt}\operatorname{Im}H
+    &=\dot u^\alpha\partial_\alpha\operatorname{Im}H+\dot{\bar u}^{\bar\alpha}\bar\partial_{\bar\alpha}\operatorname{Im}H \\
+    &=\tfrac i4\left(\bar\partial_{\bar\beta}\bar Hh^{\bar\beta\alpha}\partial_\alpha H-\partial_\beta H\bar h^{\beta\bar\alpha}\bar\partial_{\bar\alpha}\bar H  \right)\\
+    &=0
+  \end{aligned}
 \end{equation}
-These dynamics preserve $\operatorname{Im}H$ and in fact correspond to
-Hamiltonian dynamics, with the real and imaginary parts of the coordinates
-taking the role of conjugate variables. \cite{Morrow_2006_Complex}
+since $h$ is a Hermitian operator with $\bar h=h^T$.
 
-Working with a particular map is inconvenient, and we would like to develop a map-independent dynamics. Using the chain rule, one finds
+Working with a particular map is inconvenient, and we would like to develop a
+map-independent dynamics. Suppose that $z:\mathbb C^{N-1}\to M$ is a map. Using
+the chain rule, one finds
 \begin{equation}
   \begin{aligned}
-    \dot z_i
-    &=\dot u_\alpha\partial_\alpha z_i
-    =-\tfrac12(\partial_\beta H)^*g^{-1}_{\beta\alpha}\partial_\alpha z_i\\
-    &=-\tfrac12(\partial_j H)^*(\partial_\beta z_j)^*g^{-1}_{\beta\alpha}\partial_\alpha z_i
-    =-\tfrac12(\partial H)^\dagger(J^\dagger g^{-1}J)\\
+    \dot z^i
+    &=\dot u^\alpha\partial_\alpha z^i
+    =-\tfrac12\bar\partial_{\bar\beta}\bar Hh^{\bar\beta\alpha}\partial_\alpha z^i
+    =-\tfrac12\bar\partial_j\bar H\partial_{\bar\beta}\bar z^{\bar\jmath}h^{\bar\beta\alpha}\partial_\alpha z^i \\
+    &=-\tfrac12(J^\dagger h^{-1}J)^{\bar\jmath i}\bar\partial_{\bar\jmath}\bar H\\
   \end{aligned}
 \end{equation}
-where $J$ is the Jacobian of the coordinate map and $g$ is the metric. In stereographic coordinates this can be worked out directly.
-Consider the coordinates $z_i=u_i$ for $1\leq i\leq N-1$ and
+where $J$ is the Jacobian of the coordinate map and $h$ is the Hermitian
+metric. In stereographic coordinates this can be worked out directly.
+Consider the coordinates $z^i=u^i$ for $1\leq i\leq N-1$ and
 \begin{equation}
-  z_N=\sqrt{N-u_\alpha u_\alpha}
+  z^N=\sqrt{N-u^2}
 \end{equation}
 The Jacobian is
 \begin{equation}
-  J_{\alpha i}=\partial_\alpha z_i=\delta_{\alpha i}-\delta_{Ni}\frac{u_\alpha}{\sqrt{N-u_\beta u_\beta}}
+  J_\alpha^{\hphantom\alpha i}=\partial_\alpha z^i=\delta_\alpha^{\hphantom\alpha i}-\delta_N^{\hphantom Ni}\frac{u_\beta}{\sqrt{N-u^2}}
 \end{equation}
-and the corresponding hermitian metric is
+and therefore the Hermitian metric induced by the map is
 \begin{equation}
-  g_{\alpha\beta}=J_{i\alpha}^*J_{i\beta}
-  =\delta_{\alpha\beta}+\frac{u_\alpha^*u_\beta}{|N-u_\gamma u_\gamma|}
+  h_{\alpha\bar\beta}=\bar J_{\alpha}^{\hphantom\alpha i}J_{\bar\beta}^{\hphantom\beta\bar\jmath}\delta_{i\bar\jmath}
+  =\delta_{\bar\alpha\beta}+\frac{\bar u^{\alpha}u^{\bar\beta}}{|N-u^2|}
 \end{equation}
 The metric can be inverted explicitly:
 \begin{equation}
-  g^{-1}_{\alpha\beta}
-  =\delta_{\alpha\beta}-\frac{u_\alpha^*u_\beta}{|N-u_\gamma u_\gamma|+|u|^2}.
+  h^{\bar\beta\alpha}
+  =\delta^{\bar\beta\alpha}-\frac{\bar u^{\bar\beta}u^\alpha}{|N-u^2|+|u|^2}.
 \end{equation}
 Putting these pieces together, we find
 \begin{equation}
-  (J^\dagger g^{-1}J)_{ij}
-  =\delta_{ij}-\frac{z_iz_j^*}{|z|^2}
+  (J^\dagger h^{-1}J)^{\bar\jmath i}
+  =\delta^{\bar\jmath i}-\frac{z^{\bar\jmath}\bar z^i}{|z|^2}
 \end{equation}
-which is just the projector onto the constraint manifold.
-
-Therefore, a map-independent dynamics is given by
+which is just the projector onto the constraint manifold \cite{Morrow_2006_Complex}.
+Therefore, a map-independent dynamics for $z\in M$ is given by
 \begin{equation}
   \dot z
   =-\frac12(\partial H)^\dagger\left(I-\frac{zz^\dagger}{|z|^2}\right)
@@ -132,15 +141,35 @@ The level sets defined by $\operatorname{Im}H=c$ for $c\in\mathbb R$ are surface
 $2(N-1)-1$ real dimensions. They must be simply connected, since gradient
 descent in $\operatorname{Re}H$ cannot pass the same point twice.
 
+
+
 \section{2-spin}
 
+The Hamiltonian of the $2$-spin model is defined for $z\in\mathbb C^N$ by
 \begin{equation}
-  H_0=\frac12z^TJz
+  H_0=\frac12z^TJz.
 \end{equation}
 $J$ is generically diagonalizable by a complex orthogonal matrix $P$. With
-$z\mapsto Pz$, $J\mapsto D$ for diagonal $D$. Then $\partial_i H=d_iz_i$.
+$z\mapsto Pz$, $J\mapsto D$ for diagonal $D$. Then $\partial_i H=d_iz_i$. We will henceforth assume to be working in this basis. To find the critical points, we use the method of Lagrange multipliers. Introducing $\epsilon$, the constrained Hamiltonian is
+\begin{equation}
+  H=H_0+\epsilon(N-z^2)
+\end{equation}
+As usual, $\epsilon$ is equivalent to the energy per spin at any critical point.
+Critical points must satisfy
+\begin{equation}
+  0=\partial_iH=(d_i-2\epsilon)z_i
+\end{equation}
+which is only possible for $z_i=0$ or $\epsilon=\frac12 d_i$. Generically the $d_i$ will all differ, so this can only be satisfied for one $d_i$ at a time, and to be a critical point all other $z_j$ must be zero. In the direction in question,
+\begin{equation}
+  \epsilon=\frac1N\frac12d_iz_i^2=\frac12 d_i,
+\end{equation}
+whence $z_i=\pm\sqrt{N}$. Thus there are $2N$ critical points, each corresponding to $\pm$ the cardinal directions in the diagonalized basis.
+
 Suppose that two critical points have the same imaginary energy; without loss
-of generality, assume these are the first and second components. Since the gradient is proportional to $z$, any components that are zero at some time will be zero at all times. The dynamics for the components of interest assuming all others are zero are
+of generality, assume these are associated with the first and second
+cardinal directions. Since the gradient is proportional to $z$, any components that are
+zero at some time will be zero at all times. The dynamics for the components of
+interest assuming all others are zero are
 \begin{equation}
   \begin{aligned}
     \dot z_1
@@ -148,8 +177,9 @@ of generality, assume these are the first and second components. Since the gradi
     &=-(d_1-d_2)^*z_1^*\frac{|z_2|^2}{|z_1|^2+|z_2|^2}
   \end{aligned}
 \end{equation}
-and the same for $z_2$ with all indices swapped.  Since $\Delta=d_1-d_N$ is
-real, if $z_1$ begins real it remains real, with the same for $z_2$. Since the critical points are at real $z$, we make this restriction, and find
+and the same for $z_2$ with all indices swapped.  Since $\Delta=d_1-d_2$ is
+real, if $z_1$ begins real it remains real, with the same for $z_2$. Since the
+critical points are at real $z$, we make this restriction, and find
 \begin{equation}
   \begin{aligned}
     \frac d{dt}(z_1^2+z_2^2)=0\\
-- 
cgit v1.2.3-70-g09d2


From aa6d18593563f764380859b3f8b522bbb112ac07 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Tue, 9 Mar 2021 16:36:09 +0100
Subject: Removed perhaps suprious statement.

---
 stokes.tex | 4 +---
 1 file changed, 1 insertion(+), 3 deletions(-)

(limited to 'stokes.tex')

diff --git a/stokes.tex b/stokes.tex
index 1ac2298..7dab4da 100644
--- a/stokes.tex
+++ b/stokes.tex
@@ -48,9 +48,7 @@ $p_\kappa$ is taken to be Gaussian or discrete on $\pm1$.
 
 We are interested in analytically continuing expressions like $\overline F$
 into the region of complex $\beta$ or distributions $p_\kappa$ involving
-complex $J$. The former has been considered extensively for the Gaussian
-$p$-spin in the past \cite{complex_energy}, while the latter is largely
-unexplored.
+complex $J$. 
 
 When the argument of the exponential integrand in \eqref{eq:partition.function}
 acquires an imaginary component, various numeric and perturbative schemes for
-- 
cgit v1.2.3-70-g09d2


From 4c1e9aea8011aa6185a351a407594ce020d63ac7 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Thu, 11 Mar 2021 12:26:28 +0100
Subject: Lots of writing, and a figure!

---
 figs/local_flow.pdf | Bin 0 -> 50712 bytes
 stokes.tex          | 164 ++++++++++++++++++++++++++++++++++++++--------------
 2 files changed, 120 insertions(+), 44 deletions(-)
 create mode 100644 figs/local_flow.pdf

(limited to 'stokes.tex')

diff --git a/figs/local_flow.pdf b/figs/local_flow.pdf
new file mode 100644
index 0000000..dc662ad
Binary files /dev/null and b/figs/local_flow.pdf differ
diff --git a/stokes.tex b/stokes.tex
index 7dab4da..41878b7 100644
--- a/stokes.tex
+++ b/stokes.tex
@@ -62,55 +62,38 @@ The surface $M\subset\mathbb C^N$ defined by $N=f(z)=z^2$ is an $N-1$ dimensiona
 holomorphic function \cite{Forstneric_2017_Stein}. One can define a Hermitian
 metric $h$ on $M$ by taking the restriction of the standard metric of $\mathbb
 C^N$ to vectors tangent along $M$. For any smooth function $\phi:M\to\mathbb
-R$, its gradient $\nabla\phi$ is a holomorphic vector field given by
+R$, its gradient is a holomorphic vector field given by
 \begin{equation}
-  \nabla^\alpha\phi=h^{\bar\beta\alpha}\bar\partial_{\bar\beta}\phi
+  \operatorname{grad}\phi=\bar\partial^\sharp\phi
 \end{equation}
-Dynamics consists of gradient descent on $\operatorname{Re}H$, or
-\begin{equation} \label{eq:flow}
-  \dot u^\alpha=-\nabla^\alpha\operatorname{Re}H=-\tfrac12h^{\bar\beta\alpha}\bar\partial_{\bar\beta}\bar H
+Suppose $u:M\to\mathbb C^{N-1}$ is a coordinate system. Then
+\begin{equation}
+  \operatorname{grad}\phi=h^{\bar\beta\alpha}\bar\partial_{\bar\beta}\phi\frac{\partial}{\partial u^\alpha}
 \end{equation}
-Gradient descent on $\operatorname{Re}H$ is equivalent to Hamiltonian dynamics
-with the Hamiltonian $\operatorname{Im}H$. This is because $(M, h)$ is Kähler
-and therefore admits a symplectic structure, but that the flow conserves
-$\operatorname{Im}H$ can be seen from the Cauchy--Riemann equations and
-\eqref{eq:flow}:
+Let $z=u^{-1}$.
 \begin{equation}
-  \begin{aligned}
-    \frac d{dt}\operatorname{Im}H
-    &=\dot u^\alpha\partial_\alpha\operatorname{Im}H+\dot{\bar u}^{\bar\alpha}\bar\partial_{\bar\alpha}\operatorname{Im}H \\
-    &=\tfrac i4\left(\bar\partial_{\bar\beta}\bar Hh^{\bar\beta\alpha}\partial_\alpha H-\partial_\beta H\bar h^{\beta\bar\alpha}\bar\partial_{\bar\alpha}\bar H  \right)\\
-    &=0
-  \end{aligned}
+  \frac\partial{\partial u^\alpha}=\frac{\partial z^i}{\partial u^\alpha}\frac\partial{\partial z^i}
 \end{equation}
-since $h$ is a Hermitian operator with $\bar h=h^T$.
-
-Working with a particular map is inconvenient, and we would like to develop a
-map-independent dynamics. Suppose that $z:\mathbb C^{N-1}\to M$ is a map. Using
-the chain rule, one finds
 \begin{equation}
-  \begin{aligned}
-    \dot z^i
-    &=\dot u^\alpha\partial_\alpha z^i
-    =-\tfrac12\bar\partial_{\bar\beta}\bar Hh^{\bar\beta\alpha}\partial_\alpha z^i
-    =-\tfrac12\bar\partial_j\bar H\partial_{\bar\beta}\bar z^{\bar\jmath}h^{\bar\beta\alpha}\partial_\alpha z^i \\
-    &=-\tfrac12(J^\dagger h^{-1}J)^{\bar\jmath i}\bar\partial_{\bar\jmath}\bar H\\
-  \end{aligned}
+  \bar\partial_{\bar\beta}\phi=\frac{\partial\bar z^i}{\partial\bar u^{\bar\beta}}\frac{\partial\phi}{\partial\bar z^i}
+\end{equation}
+\begin{equation}
+  \operatorname{grad}\phi
+  =\frac{\partial\bar z^{\bar\jmath}}{\partial\bar u^{\bar\beta}}h^{\bar\beta\alpha}\frac{\partial z^i}{\partial u^\alpha}\frac{\partial\phi}{\partial\bar z^{\bar\jmath}}\frac\partial{\partial z^i}
 \end{equation}
-where $J$ is the Jacobian of the coordinate map and $h$ is the Hermitian
-metric. In stereographic coordinates this can be worked out directly.
-Consider the coordinates $z^i=u^i$ for $1\leq i\leq N-1$ and
+At any point $z_0\in M$, $z_0$ is parallel to the normal to $M$. Without loss of generality, take $z_0=e^N$. In a neighborhood of $z_0$, the map $u^\alpha=z^\alpha$ is a coordinate system.
+its inverse is $z^i=u^i$ for $1\leq i\leq N-1$ and
 \begin{equation}
-  z^N=\sqrt{N-u^2}
+  z^N=\sqrt{N-u^2}.
 \end{equation}
 The Jacobian is
 \begin{equation}
-  J_\alpha^{\hphantom\alpha i}=\partial_\alpha z^i=\delta_\alpha^{\hphantom\alpha i}-\delta_N^{\hphantom Ni}\frac{u_\beta}{\sqrt{N-u^2}}
+  \frac{\partial z^i}{\partial u^\alpha}=\delta^i_\alpha-\delta_N^{i}\frac{u^\alpha}{\sqrt{N-u^2}}
 \end{equation}
 and therefore the Hermitian metric induced by the map is
 \begin{equation}
-  h_{\alpha\bar\beta}=\bar J_{\alpha}^{\hphantom\alpha i}J_{\bar\beta}^{\hphantom\beta\bar\jmath}\delta_{i\bar\jmath}
-  =\delta_{\bar\alpha\beta}+\frac{\bar u^{\alpha}u^{\bar\beta}}{|N-u^2|}
+  h_{\alpha\bar\beta}=\frac{\partial\bar z^i}{\partial\bar u^\alpha}\frac{\partial z^{\bar\jmath}}{\partial u^{\bar\beta}}\delta_{i\bar\jmath}
+  =\delta_{\alpha\bar\beta}+\frac{\bar u^{\alpha}u^{\bar\beta}}{|N-u^2|}
 \end{equation}
 The metric can be inverted explicitly:
 \begin{equation}
@@ -119,26 +102,119 @@ The metric can be inverted explicitly:
 \end{equation}
 Putting these pieces together, we find
 \begin{equation}
-  (J^\dagger h^{-1}J)^{\bar\jmath i}
+  \frac{\partial\bar z^{\bar\jmath}}{\partial\bar u^{\bar\beta}}h^{\bar\beta\alpha}\frac{\partial z^i}{\partial u^\alpha}
   =\delta^{\bar\jmath i}-\frac{z^{\bar\jmath}\bar z^i}{|z|^2}
 \end{equation}
-which is just the projector onto the constraint manifold \cite{Morrow_2006_Complex}.
-Therefore, a map-independent dynamics for $z\in M$ is given by
 \begin{equation}
-  \dot z
+  \operatorname{grad}\phi
+  =\left(\delta^{\bar\jmath i}-\frac{z^{\bar\jmath}\bar z^i}{|z|^2}\right)
+  \frac{\partial\phi}{\partial\bar z^{\bar\jmath}}\frac\partial{\partial z^i}
+\end{equation}
+
+\section{Dynamics}
+
+For a holomorphic Hamiltonian $H$, dynamics are defined by gradient descent on $\operatorname{Re}H$, or
+\begin{equation} \label{eq:flow}
+  \dot z=-\operatorname{grad}\operatorname{Re}H
   =-\frac12(\partial H)^\dagger\left(I-\frac{zz^\dagger}{|z|^2}\right)
+  =-\tfrac12(\partial H)^\dagger P,
+\end{equation}
+where $P=I-\hat z\hat z^\dagger$ is the projection matrix on to the tangent space of $M$.
+Gradient descent on $\operatorname{Re}H$ is equivalent to Hamiltonian dynamics
+with the Hamiltonian $\operatorname{Im}H$. This is because $(M, h)$ is Kähler
+and therefore admits a symplectic structure, but that the flow conserves
+$\operatorname{Im}H$ can be seen from the Cauchy--Riemann equations and
+\eqref{eq:flow}:
+\begin{equation}
+  \begin{aligned}
+    \frac d{dt}\operatorname{Im}H
+    &=\dot z_i\partial_i\operatorname{Im}H+\dot{\bar z}_i\bar\partial_i\operatorname{Im}H \\
+    &=\frac i4\left(
+      \bar\partial_j\bar HP_{ji}\partial_i H-\partial_j H\bar P_{ji}\bar\partial_i\bar H
+    \right)
+    =0
+  \end{aligned}
+\end{equation}
+since $P$ is a Hermitian operator. This conservation law indicates that surfaces of constant $\operatorname{Im}H$ will be important when evaluting the possible endpoints of dynamic trajectories.
+
+Let us consider the generic case, where the critical points of $H$ have
+distinct energies. What is the topology of the $C=\operatorname{Im}H$ level
+set? We shall argue its form by construction. Consider initially the situation
+in the absence of any critical point. In this case the level set consists of a
+single simply connected surface, locally diffeomorphic to $\mathbb R^{2(N-1)-1}$. Now, `place' a generic
+(nondegenerate) critical point in the function at $z_0$. In the vicinity of the critical
+point, the flow is locally
+\begin{equation}
+  \begin{aligned}
+    \dot z_i
+    &\simeq-\frac12\left[
+      \partial_j\left(
+        H(z_0)+\frac12\partial_k\partial_\ell H(z_0)\Delta z_k\Delta z_\ell
+      \right)
+    \right]^* P_{ji} \\
+    &=-\frac12\left(
+      \partial_j\partial_kH(z_0)\Delta z_k
+    \right)^* P_{ji} \\
+    &=-\frac12\Delta z_k^*(\partial_k\partial_jH(z_0))^*P_{ji}
+  \end{aligned}
+\end{equation}
+The matrix $(\partial\partial H)^\dagger P$ has a spectrum identical to that of
+$\partial\partial H$ save a single marginal direction corresponding to $z_0$,
+the normal to the constraint surface. Assuming we are working in a diagonal basis, we find
+\begin{equation}
+  \dot z_i=-\frac12\lambda_i\Delta z_i^*+O(\Delta z^2)
 \end{equation}
+Breaking into real and imaginary parts gives
+\begin{equation}
+  \begin{aligned}
+    \frac{d\Delta x_i}{dt}&=-\frac12\left(
+      \operatorname{Re}\lambda_i\Delta x_i+\operatorname{Im}\lambda_i\Delta y_i
+    \right) \\
+    \frac{d\Delta y_i}{dt}&=-\frac12\left(
+      \operatorname{Im}\lambda_i\Delta x_i-\operatorname{Re}\lambda_i\Delta y_i
+    \right)
+  \end{aligned}
+\end{equation}
+Therefore, in the complex plane defined by each eigenvector of $(\partial\partial H)^\dagger P$ there is a separatrix flow of the form in Figure \ref{fig:local_flow}. The effect of these separatrices in each complex direction of the tangent space $T_{z_0}M$ is to separate that space into four quadrants: two disconnected pieces with greater imaginary part than the critical point, and two with lesser imaginary part. This partitioning implies that the level set of $\operatorname{Im} H=C$ for $C\neq\operatorname{Im}H(z_0)$ is splint into two disconnected pieces, one lying in each of two quadrants corresponding with its value relative to that at the critical point.
+
+\begin{figure}
+  \includegraphics[width=\columnwidth]{figs/local_flow.pdf}
+  \caption{
+    Gradient descent in the vicinity of a critical point, in the $z$--$z*$
+    plane for an eigenvector $z$ of $(\partial\partial H)^\dagger P$. The flow
+    lines are colored by the value of $\operatorname{Im}H$.
+  } \label{fig:local_flow}
+\end{figure}
+
+Continuing to `insert' critical points whose imaginary energy differs from $C$,
+one repeatedly partitions the space this way with each insertion. Therefore,
+for the generic case with $\mathcal N$ critical points, with $C$ differing in
+value from all critical points, the level set $\operatorname{Im}H=C$ has
+$\mathcal N+1$ connected components, each of which is simply connected,
+diffeomorphic to $\mathbb R^{2(N-1)-1}$ and connects two sectors of infinity to
+each other.
+
+When $C$ is brought to the same value as the imaginary part of some critical
+point, two of these disconnected surfaces pinch grow nearer and pinch together
+at the critical point when $C=\operatorname{Im}H$, as in the black lines of
+Figure \ref{fig:local_flow}. The unstable manifold of the critical point, which
+corresponds with the portion of this surface that flows away, is known as a
+\emph{Lefshetz thimble}.
 
 Stokes lines are trajectories that approach distinct critical points as time
 goes to $\pm\infty$. From the perspective of dynamics, these correspond to
 \emph{heteroclinic orbits}. What are the conditions under which Stokes lines
 appear? Because the dynamics conserves imaginary energy, two critical points
-must have the same imaginary energy if they are to be connected by a Stokes line.
-
-The level sets defined by $\operatorname{Im}H=c$ for $c\in\mathbb R$ are surfaces of
-$2(N-1)-1$ real dimensions. They must be simply connected, since gradient
-descent in $\operatorname{Re}H$ cannot pass the same point twice.
+must have the same imaginary energy if they are to be connected by a Stokes
+line. This is not a generic phenomena, but will happen often as one model
+parameter is continuously varied. When two critical points do have the same
+imaginary energy and $C$ is brought to that value, the level set
+$C=\operatorname{Im}H$ sees formally disconnected surfaces pinch together in
+two places. We shall argue that generically, a Stokes line will exist whenever
+the two critical points in question lie on the same connected piece of this
+surface.
 
+What are the ramifications of this for disordered Hamiltonians? When some process brings two critical points to the same imaginary energy, whether a Stokes line connects them depends on whether the points are separated from each other by the separatrices of one or more intervening critical points. Therefore, we expect that in regions where critical points with the same value of $\operatorname{Im}H$ tend to be nearby, Stokes lines will proliferate, while in regions where critical points with the same value of $\operatorname{Im}H$ tend to be distant compared to those with different $\operatorname{Im}H$, Stokes lines will be rare.
 
 
 \section{2-spin}
-- 
cgit v1.2.3-70-g09d2


From c08bcb900b7524a1b474062ece3b23dfe76d350f Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Thu, 11 Mar 2021 15:08:04 +0100
Subject: More writing.

---
 stokes.tex | 12 ++++++++++++
 1 file changed, 12 insertions(+)

(limited to 'stokes.tex')

diff --git a/stokes.tex b/stokes.tex
index 41878b7..d07ec6a 100644
--- a/stokes.tex
+++ b/stokes.tex
@@ -271,6 +271,18 @@ The critical points of the 2-spin model are all adjacent: no critical point is s
 
 \section{(2 + 4)-spin}
 
+\section{Numerics}
+
+To confirm the presence of Stokes lines under certain processes in the $p$-spin, we studied the problem numerically.
+\begin{equation}
+  \mathcal L(z(t), z'(t))
+  = 1-\frac{\operatorname{Re}(\dot z^*_iz'_i)}{|\dot z||z'|}
+\end{equation}
+\begin{equation}
+  \mathcal C[\mathcal L]=\int dt\,\mathcal L(z(t), z'(t))
+\end{equation}
+$\mathcal C$ has minimum of zero, which is reached only by functions $z(t)$ whose tangent is everywhere parallel to the direction of the dynamics. Therefore, finding functions which reach this minimum will find time-reparameterized Stokes lines. For the sake of numeric stability, we look for functions whose tangent has constant norm.
+
 \begin{acknowledgments}
   MIT mathematicians have been no help
 \end{acknowledgments}
-- 
cgit v1.2.3-70-g09d2

