The Archimedean limit of random sorting networks
A sorting network (also known as a reduced decomposition of the reverse permutation), is a shortest path from $12 \cdots n$ to $n \cdots 21$ in the Cayley graph of the symmetric group $S_n$ generated by adjacent transpositions. We prove that in a uniform random $n$-element sorting network $\sigma^n$, all particle trajectories are close to sine curves with high probability. We also find the weak limit of the time-$t$ permutation matrix measures of $\sigma^n$. As a corollary of these results, we show that if $S_n$ is embedded into $\mathbb{R}^n$ via the map $\tau \mapsto (\tau(1), \tau(2), \dots \tau(n))$, then with high probability, the path $\sigma^n$ is close to a great circle on a particular $(n-2)$-dimensional sphere in $\mathbb{R}^n$. These results prove conjectures of Angel, Holroyd, Romik, and Virag.
\section{Introduction}
\label{S:intro}

Let $\Gamma(S_n)$ be the Cayley graph of the symmetric group  $S_n$ with generators given by the adjacent transpositions $\pi_i = (i, i + 1)$ for $i = 1, \mathellipsis, n -1$. A {\bf sorting network} is a shortest path from the identity permutation $\text{id}_n = 12 \cdots n$ to the reverse permutation $\text{rev}_n = n \cdots 21$ in  $\Gamma(S_n)$. The length of any sorting network is $N = {n \choose 2}$.

\medskip

Sorting networks are known as {\bf reduced decompositions} of the reverse permutation, as any sorting network can equivalently be represented as a minimal length decomposition of the reverse permutation as a product of adjacent transpositions: $\text{rev}_n = \pi_{k_N} \mathellipsis \pi_{k_1}$. The combinatorics of sorting networks have been studied in detail under this name. There are connections between sorting networks and Schubert calculus, quasisymmetric functions, zonotopal tilings of polygons, and aspects of representation theory. For more background on these topics, see Stanley \cite{stanley1984number}; Bj\"orner and Brenti \cite{bjorner2006combinatorics}; Garsia \cite{garsia2002saga}; Tenner \cite{tenner2006reduced}; and Manivel \cite{manivel2001symmetric}.

\begin{figure}
   \centering
   \includegraphics[scale=0.7]{Wiring-3.pdf}
   \caption{A ``wiring diagram" for a sorting network with $n = 4$. In
     this diagram, trajectories are drawn as continuous curves for
     clarity.}
   \label{fig:wiring}
\end{figure}

\medskip

In computer science, sorting networks are viewed as $N$-step algorithms for sorting a list of $n$ numbers. At step $i$, the sorting network algorithm sorts the numbers at positions $k_i$ and $k_{i+1}$ into increasing order. This process sorts any list in $N$ steps. 

\medskip

As in this algorithmic interpretation of sorting networks, we will think of the elements $\{1, \mathellipsis, n \}$ as particles being sorted in time (see Figure \ref{fig:wiring}). We use the notation $\sigma(x, t) = \pi_{k_\floor{t}} \mathellipsis \pi_{k_2}\pi_{k_1}(x)$ for the position of particle $x$ at time $t$. We call $(k_1, \dots, k_N)$ the \textbf{swap sequence} for $\sigma$.

\subsection{Main limit theorems}
\label{SS:main-thms}

Angel, Holroyd, Romik, and Vir\'ag \cite{angel2007random} initiated the study of uniform random $n$-element sorting networks. They studied the limiting behaviour of the space-time swap distribution, rescaled particle trajectories, time-$t$ permutation matrices, and the Cayley graph path itself.

\medskip

They proved a law of large numbers for the space-time swap distribution, and based on strong numerical evidence, made conjectures about the limiting behaviour of the other three objects. In this paper, we prove these conjectures.



\bigskip

\noindent {\bf I. Rescaled particle trajectories.} \qquad  For a sorting network $\sigma$, define the global trajectory
$$
\sigma_G(x, t) = \frac{2\sigma(x, Nt)}n - 1.
$$
The function $\sigma_G(x, \cdot)$ is the trajectory of particle $x$, with time rescaled so that the sorting process finishes at time $1$, and space rescaled so that the trajectory stays in the interval $[-1, 1]$.
In \cite{angel2007random}, the authors conjectured that with high probability, all particle trajectories in a uniform random sorting network are close to sine curves  (see Figure \ref{fig:sinecurves}). 

\medskip

They proved that limiting trajectories are $1/2$-H\"older with H\"older constant $\sqrt{8}$. Dauvergne and Vir\'ag \cite{dauvergne1} improved upon this by showing that limiting trajectories are $\pi$-Lipschitz, and that trajectories of particles starting within $\epsilon$ of the edge are distance $2 \sqrt{\epsilon}$ away from a sine curve in uniform norm.

\medskip

Our first theorem proves the sine curve conjecture from \cite{angel2007random}. Here and throughout the paper we use the notation $\sigma^n$ for a uniform random $n$-element sorting network.

\begin{customthm}{1}[Sine curve limit]
\label{T:sine-curves}
For each $n$ there exist random variables $\{(A^n_i, \Theta^n_i) \in [0, 1] \times [0, 2\pi]\}_{i \in \{1, \mathellipsis, n\}}$ such that for any $\epsilon > 0$, we have that
$$
\mathbb{P} \left( \max_{i \in [1, n] } \max_{t \in [0, 1]} \card{ \sigma^n_G(i, t) - A_i^n \sin(\pi t + \Theta_i^n) } > \epsilon \right) \to 0 \qquad \;\text{as}\; \;\; n \to \infty.
$$
\end{customthm}

\begin{figure}
   \centering
   \includegraphics[scale=0.8]{RSN-sinecurves.pdf}
   \caption{This is a diagram of selected particle trajectories in a 2000 element sorting network. This image is taken from \cite{angel2007random}.}
   \label{fig:sinecurves}
\end{figure}

\bigskip

\noindent {\bf II. Permutation Matrices.} \qquad For a uniform $n$-element sorting network $\sigma^n$, define 
$$
\rho^n_t= \frac{1}n \sum_{i=1}^n \delta \left( \sigma^n_G(i, 0), \sigma^n_G(i, t) \right).
$$
Here $\delta(x, y)$ is a $\delta$-mass at the point $(x, y)$.
The measure $\rho^n_t$ rescales the time-$t$ permutation matrix of $\sigma^n$, placing atoms of weight $1/n$ at the positions of the $1$'s. Define the {\bf Archimedean measure} $\mathfrak{Arch}$ on the square $[-1, 1]^2$ to be the measure with Lebesgue density 
$$
f(x, y) = \frac{1}{2\pi\sqrt{1 - x^2 - y^2}}
$$
on the unit ball $B(0, 1)$, and $0$ elsewhere. This is simply the projected surface area measure of the $2$-sphere. Define the {\bf time-$t$ Archimedean measure} by
$$
(X, X \cos (\pi t) + Y \sin (\pi t)) \stackrel{d}{=} \mathfrak{Arch}_{t}, \qquad \text{ where } \;\; (X, Y) \stackrel{d}{=} \mathfrak{Arch}.
$$
Note that $\mathfrak{Arch} = \mathfrak{Arch}_{1/2}$.
In \cite{angel2007random}, the authors conjectured that $\rho^n_t$ converges weakly to $\mathfrak{Arch}_{t}$ (see Figure \ref{fig:circles}). They proved that for any $t$, the support of $\rho^n_t$ lies in a particular octagon with high probability. In \cite{dauvergne1}, the authors showed that the support of any limit of $\rho^n_t$ lies within the elliptical support of $\mathfrak{Arch}_t$. 
\begin{figure}
   \centering
   \includegraphics[scale=0.9]{RSN-circles.pdf}
   \caption{This is a diagram of the measures $\{\rho^n_t : t \in \{0, 1/10, 2/10, \mathellipsis 1\} \}$ in a $500$-element sorting network. The blue octagons are the octagonal support bounds proved in \cite{angel2007random}, and the red ellipses are the support bounds proved in \cite{dauvergne1}. This figure is from \cite{angel2007random}.}
   \label{fig:circles}
\end{figure}

\medskip

Our second theorem proves the weak convergence of the random measures $\rho^n_t$. We also show that the support of $\rho^n_t$ and $\mathfrak{Arch}_t$ are close. To state this theorem, recall that the Hausdorff distance between two sets $A, B \subset \mathbb{R}^2$ is
$$
d_H(A, B) = \max \left\{\sup_{a \in A} \inf_{b \in B} d(a, b), \;\sup_{b \in B} \inf_{a \in A} d(a, b) \right\}.
$$
\begin{customthm}{2}[Permutation matrix limit]
\label{T:matrices}
For any $t \in [0, 1]$,
$
\rho^n_t \to \mathfrak{Arch}_t
$ 
in probability
as $n \to \infty$. That is, for any weakly open set $U$ in the space of probability measures on $[-1, 1]^2$ containing $\mathfrak{Arch}_t$, we have that
$$
\mathbb{P}(\rho^n_t \in U) \to 1 \qquad \;\text{as}\; \quad n \to \infty.
$$
Moreover, 
$$
d_H(\text{supp}(\rho^n_t), \text{supp}(\mathfrak{Arch}_t)) \to 0 \qquad \text{ in probability } \;\text{as}\; n \to \infty.
$$
\end{customthm}

Angel, Holroyd, Romik, and Vir\'ag \cite{angel2007random} also considered a permutation matrix evolution for $\sigma^n$. Let $j \in \{1, \mathellipsis, n\}$, and consider the random complex-valued function
$$
Z^n_j(t) = e^{\pi i t} \left[\sigma^n_G(j, t) + i\sigma^n_G(j, t + 1/2) \right], \qquad t \in [0, 1/2].
$$
For a fixed $t$, $(Z^n_1(t), \mathellipsis, Z^n_n(t))$ is the set of points in the scaled permutation matrix for 
\[
\sigma^n(\cdot, t + 1/2)(\sigma^n)^{-1}(\cdot, t)
\]
after a counterclockwise rotation by $\pi t$ (see Figure \ref{fig:window}). Theorem \ref{T:sine-curves} guarantees that each of the paths $Z^n_j$ localize.
\begin{customthm}{3}[Path Localization]
\label{T:unif-rotation}
Let $\sigma^n$ be a uniform random $n$-element sorting network.
Then 
$$
\max_{j \in [1, n]} \max_{s, t \in [0, 1]} |Z^n_j(t) - Z^n_j(s)| \to 0 \qquad \text{in probability as } \;\; n \to \infty.
$$
\end{customthm}





\bigskip

\noindent {\bf III. Great Circles.} \qquad We can embed the vertices of $\Gamma(S_n)$ into $\mathbb{R}^n$ by sending the permutation $\tau \in S_n$ to the point $\close{\tau} = (\tau(1), \mathellipsis, \tau(n)) \in \mathbb{R}^n$. For any $\tau \in S_n$, the point $\close{\tau}$ lies on the $(n-2)$-sphere $\mathbb{S}^{n-2} = \mathbb{L}_n \cap \mathbb{K}_n$, where
\begin{align*}
\mathbb{L}_n =& \left\{ (x_1, \mathellipsis x_n) \in \mathbb{R}^n : \sum_{i=1}^n x_i = \frac{n(n+1)}2 \right\} \qquad \;\text{and}\; \\
\mathbb{K}_n =& \left\{ (x_1, \mathellipsis x_n) \in \mathbb{R}^n : \sum_{i=1}^n x_i^2 = \frac{n(n+1)(2n+1)}6\right\}.
\end{align*}
If we also embed the edges of $\Gamma(S_n)$ as straight lines between the embedded vertices, we get an object called the {\bf permutahedron} (see Figure \ref{fig:permutahedron}). 

\medskip



\begin{figure}
   \centering
   \includegraphics[scale=0.3]{Permutahedron.png}
   \caption{The permutahedron for $S_4$.}
   \label{fig:permutahedron}
\end{figure}
In \cite{angel2007random}, the authors conjectured that a uniform random sorting network is close to a great circle in $\mathbb{S}^{n-2}$ under this embedding. They showed that this conjecture implies the other global limiting results for uniform random sorting networks. Our strongest theorem proves this conjecture.

\medskip

For two functions $f, g: [0, 1] \to \mathbb{R}^n$,  define the distance function
$$
d_\infty(f, g) = \sup_{t \in [0, 1]} ||f(t) - g(t)||_\infty.
$$
This is the uniform norm on functions, where the pointwise distance is the $L^\infty$ distance. Also, for an $n$-element sorting network $\sigma$, we define the embedded and time-rescaled path
$
\bar{\sigma}:[0, 1] \to \mathbb{S}^{n-2}
$
by letting 
$$
\bar{\sigma}(t) = (\sigma(1, Nt), \sigma(2, Nt), \dots, \sigma(n, Nt)).
$$

\begin{customthm}{4}[Great circles]
\label{T:geom-limit}
Let $\bar{\sigma}^n$ be the embedding of $\sigma^n$ into $\mathbb{S}^{n-2}$. For every $n$, there exists a random path $C_n:[0, 1] \to \mathbb{S}^{n-2}$ such that $C_n$ is a constant-velocity parametrization of an arc of a great circle in $\mathbb{S}^{n-2}$ starting at $(1, \mathellipsis, n)$ and finishing at $(n, \mathellipsis, 1)$, and such that 
$$
\frac{d_\infty(\bar{\sigma}^n, C_n)}{n} \to 0
\qquad \text{in probability as} \;\;  n \to \infty.
$$
\end{customthm}

Note that it is easy to find sorting networks that aren't close to a great circle in $\mathbb{S}^{n-2}$. For example, the ``bubble sort" sorting network given by the swap sequence 
$$
(1, 2, \mathellipsis n-1, 1, 2, \mathellipsis, n-2, \mathellipsis, 1, 2, 1)
$$
is $d_\infty$-distance $n - 1 - o(1)$ from any great circle. 


\subsection{Random trajectory limits}

To prove the main theorems of Section \ref{SS:main-thms}, we first analyze the limit of a random particle trajectory. This approach was first considered by Rahman, Vir\'ag, and Vizer \cite{rahman2016geometry}, and continued in \cite{dauvergne1}.  

\medskip

Let $\mathcal{D}$ be the closure in the uniform norm $||\cdot||_u$ of the space of all possible sorting network trajectories $\sigma_G(x, \cdot):[0, 1] \to [-1, 1]$. The space $(\mathcal{D}, ||\cdot||_u)$ is a complete separable metric space. The only functions in $\mathcal{D}$ are continuous functions and the sorting network trajectories themselves.

\medskip
Let $Y_n \in \mathcal{D}$ be a uniform $n$-element sorting network trajectory. That is, if $\sigma^n$ is a uniform $n$-element sorting network, and $I_n$ is an independent uniform random variable on $\{1, \mathellipsis n\}$, then
 $$
 Y_n = \sigma^n_G(I_n, \cdot).
 $$
We refer to $Y_n$ as the {\bf trajectory random variable} of $\sigma^n$.
In \cite{dauvergne1}, Dauvergne and Vir\'ag proved that the sequence $\{Y_n\}_{n \in \mathbb{N}}$ is precompact in distribution, and that any subsequential limit is almost surely Lipschitz (we state their results precisely in Section \ref{S:local}). In this paper, we show that $\{Y_n\}_{n \in \mathbb{N}}$ converges in distribution, and identify its limit.

\begin{customthm}{5}[The weak trajectory limit]
\label{T:weak-limit}
Let $(X, Z) \sim \mathfrak{Arch}$, and define the Archimedean path $\mathcal{A} \in \mathcal{D}$ by $\mathcal{A}(t) = X\cos(\pi t) + Z\sin(\pi t)$. Then
 $$
Y_n \stackrel{d}{\to} \mathcal{A} \qquad \;\text{as}\; \; n \to \infty.
 $$
 \end{customthm}
 
Theorem \ref{T:weak-limit} will be used in the proof of all our main theorems from Section \ref{SS:main-thms}. Most of the paper is devoted to its proof. We note here that we can equivalently write 
$$
\mathcal{A}(t) = \sqrt{1 - V^2}\sin(\pi t + 2\pi U),
$$
where $V$ and $U$ are independent uniform random variables on $[0, 1]$.


\bigskip
 
\noindent \textbf{Random $m$-out-of-$n$ sorting networks.} \qquad 
We will also use Theorem \ref{T:weak-limit} to identify the limit of random $m$-out-of-$n$ subnetworks. This answers a question of Angel and Holroyd \cite{angel2010random}. This limit can also be by found by using the stronger great circle theorem (Theorem \ref{T:geom-limit}), as was noted in \cite{angel2010random}.

\medskip

Let $\sigma$ be an $n$-element sorting network. For $A \subset \{1, \dots, n\}$, let $\sigma|_A$ be the $|A|$-element sorting network given by restricting $\sigma$ to the set $A$. Specifically, for $i \in \{1, \dots, N\}$ define $\sigma^*_A(\cdot, i)$ be the relative ordering of the particles in $A$ in the permutation $\sigma(\cdot, i)$. This gives a sequence of $N$ permutations $\{\sigma^*_A(\cdot, i) \in S_{|A|} : i \in \{1, \dots N\} \}$. Removing duplicates gives the permutation sequence for sorting network $\sigma|_A$. For $m < n$, let the \textbf{random $m$-out-of-$n$ subnetwork $\tau^n_m$} be the restriction of $\sigma^n$ to a uniform $m$-element subset of $\{1, \mathellipsis n\}$, chosen independently from $\sigma^n$. 




 \medskip
 
Let $\{x_1, \mathellipsis x_n\}$ be a set of points in $\mathbb{R}^2$ in general position, and such that no two pairs of points determine parallel lines. Label the points in order of increasing $x$-coordinate. For all but finitely many angles $\theta$, listing the labels of the points $\{x_1, x_2, \mathellipsis x_n\}$ in increasing order of their horizontal projections after rotation by $\pi \theta$  gives a permutation $\tau_\theta$. The \textbf{geometric sorting network} associated to $\{x_1, \mathellipsis x_n\}$ is simply the sequence of permutations $\{\tau_\theta, \theta \in [0, 1] \}$ listed in order of increasing $\theta$. 

 \begin{customthm}{6}[The subnetwork limit]
 \label{T:subnetwork}
 Let $\{X_1, X_2, \mathellipsis, X_m \}$ be random points in the unit ball $B(0, 1)$ sampled from the Archimedean distribution $\mathfrak{Arch}$, and let $\tau_m$ be the associated geometric sorting network. Then
 $$
\tau^n_m \stackrel{d}{\to} \tau_m \qquad \;\text{as}\; \; n \to \infty.
 $$ 
 \end{customthm}
 

 
 
 \subsection{The local speed distribution} 
 
As a by-product of the proof of Theorem \ref{T:weak-limit}, we find the distribution of speeds in the local limit of random sorting networks. To state this result, we first give an informal description of this limit (a precise description is given in Section \ref{S:local}). The existence of this limit was established independently by Angel, Dauvergne, Holroyd, and Vir\'ag \cite{angel2017local}, and by Gorin and Rahman \cite{gorin2017}. Define 
$$
U_n (x, t) = \sigma^n(\floor{n/2} + x, nt) - \floor{n/2}.
$$
Each path $U_n(x, \cdot)$ is a locally scaled particle trajectory.
With an appropriate notion of convergence, we have that 
$$
U_n \stackrel{d}{\to} U,
$$
where $U$ is a random function from $\mathds{Z} \times [0, \infty) \to \mathds{Z}$. $U$ is the local limit at the centre of the sorting network. We can also take a local limit centred at particle $\floor{\alpha n}$ for any $\alpha \in (0, 1)$. The result is the process $U$ with time rescaled by a semicircle factor $2\sqrt{\alpha(1 - \alpha)}$. 

\medskip

In \cite{dauvergne1}, Dauvergne and Vir\'ag proved that particles in $U$ have asymptotic speeds. Specifically, they showed that for every $x \in \mathds{Z}$, the limit
$$
S(x) = \lim_{t \to \infty} \frac{U(x, t) - U(x, 0)}{t} \qquad \text{exists }\text{almost surely}.
$$ 
They showed that $S(x)$ has distribution $\mu$ independent of $x$. In this paper, we identify $\mu$.

\begin{customthm}{7}
\label{T:main-2}
The measure $\mu$ is the arcsine distribution on $[-\pi, \pi]$ given by the Lebesgue density
$$
f(x) = \frac{1}{\pi\sqrt{\pi^2 - x^2}}.
$$
\end{customthm}

\subsection*{Related work}

In addition to the papers mentioned above, the local behaviour of random sorting networks has also been studied by Angel, Gorin, and Holroyd \cite{angel2012pattern}. The frequency of particular substrings in the swap sequence of a random sorting network has been analyzed by Reiner \cite{reiner2005note} and Tenner \cite{tenner2014expected}. Fulman and Stein \cite{fulman2014stein} have analyzed the distribution of the first swap in a random sorting network. 

\medskip

All of this work relies on a bijection of Edelman and Greene \cite{edelman1987balanced} between sorting networks and Young tableaux of a particular shape. Little \cite{little2003combinatorial} found another bijection between these two sets, and Hamaker and Young \cite{HY} proved that these bijections coincide.

\medskip

Problems involving limits of sorting networks under different measures have been considered by Angel, Holroyd, and Romik \cite{angel2009oriented}, and also by Young \cite{young2014markov}. Uniform ``relaxed" random sorting networks have been analyzed by Kotowski and Vir\'ag \cite{kotowski2016limits} (see also \cite{rahman2016geometry}).

\subsection*{Outline of the paper}
Most of the paper is devoted to proving Theorem \ref{T:weak-limit}. The key to proving this theorem is the notion of particle flux across a curve $h$. Heuristically, particle flux measures the number of times that particles cross $h$ in a large-$n$ sorting network. It is defined in terms of the local speed distribution $\mu$.

\medskip

 Particle flux is a useful quantity because any limit of $Y_n$ must minimize this quantity among curves $h$ with $h(0) = -h(1)$. This is due to the fact that in any sorting network, every pair of particles swaps \textit{exactly} once, whereas any particle must cross a line with $h(0) = -h(1)$ \textit{at least} once. 

\medskip

In Section \ref{S:local}, we introduce the necessary background for the paper. In Section \ref{S:flux-intro}, we formally define particle flux and establish its properties. In Section \ref{S:unique}, we then use these properties to partially characterize minimal flux paths $h$ with $h(0) = -h(1)$. In Section \ref{S:path-speed}, we use this characterization to relate the local speed distribution $\mu$ to the derivative distribution for subsequential limits of $Y_n$.
\medskip

 This allows us to derive an integral transform formula for $\mu$, which we do in Section \ref{S:integral-formula}. In Section \ref{S:transform}, we invert this integral transform in order to find $\mu$. This allows us to immediately prove Theorem \ref{T:weak-limit}. In this section, we also prove a slightly stronger version of Theorem \ref{T:weak-limit}, allowing us to establish Theorem \ref{T:subnetwork}.
 
 \medskip
 
 In the remainder of the paper, we use Theorem \ref{T:weak-limit} to establish our stronger limit theorems. In Section \ref{S:strong-limit}, we combine Theorem \ref{T:weak-limit} with bounds from \cite{angel2007random} to prove Theorems \ref{T:sine-curves}, \ref{T:matrices}, and \ref{T:unif-rotation}. Finally, in Section \ref{S:geom-limit}, we combine Theorem \ref{T:weak-limit} with Theorem \ref{T:sine-curves} to prove Theorem \ref{T:geom-limit}.


\section{Preliminaries}
\label{S:local} 
In this section, we introduce the necessary background about sorting networks. The most basic observation about uniform $n$-element sorting networks is that they exhibit a type of time-stationarity. This was first observed in \cite{angel2007random}.

\begin{theorem}
\label{T:time-stat}
Let $(K^n_1, K^n_2, \mathellipsis K^n_N)$ be the swap sequence of a random $n$-element sorting network $\sigma^n$. Then 
$$
\left(K^n_1, K^n_2, \mathellipsis K^n_N\right) \stackrel{d}{=} \left(n - K^n_N, K^n_1, K^n_2, \mathellipsis K^n_{N - 1}\right).
$$
\end{theorem}
We will repeatedly use time-stationarity of sorting networks to reduce proofs to statements about the beginning of a sorting network. Using the Edelman-Greene bijection between sorting networks and Young tableaux, Angel, Holroyd, Romik, and Vir\'ag \cite{angel2007random} also found an explicit formula for the distribution of $K^n_1$.

\begin{theorem}
\label{T:dist-k1}
Let $K^n_1$ be the location of the first swap of $\sigma^n$. For any $i \in \{1, \dots, n-1\}$, we have that
$$
\mathbb{P}(K^n_1 = i) = \frac{1}N \frac{[3\cdot5\cdot7\cdots(2i - 1)][3\cdot 5 \cdots (2(n-i) - 1)]}{[2\cdot4\cdot6\cdots(2i - 2)][2\cdot 4 \cdots (2(n-i) - 2)]} \le \frac{3}n.
$$
Moreover, if $\{i_n\}_{n \in \mathbb{N}}$ is any sequence of integers such that $i_n/n\to (\alpha + 1)/2$ for some $\alpha \in (-1, 1)$, then
$$
n\mathbb{P}(K^n_1 = i_n) \to \frac{4\sqrt{1 - \alpha^2}}{\pi} \qquad \;\text{as}\; n \to \infty.
$$
\end{theorem}

\bigskip

\noindent \textbf{The local limit.} \qquad  Define a {\bf swap function} as a map $V:\mathbb{Z} \times [0, \infty) \to \mathds{Z}$ with the following properties:

\smallskip

 \begin{enumerate}[nosep,label=(\roman*)]
  \item For each $x$, $V(x,\cdot)$ is cadlag with nearest neighbour jumps.
  \item For each $t$,  $V(\cdot,t)$ is a bijection from $\mathds{Z}$  to $\mathds{Z}$.
  \item Define $V^{-1}(x, t)$ by $V(V^{-1}(x, t),t) = x$.  Then for each $x$, $V^{-1}(x, \cdot)$
    is a cadlag path with nearest neighbour jumps.
  \item For any time $t \in (0, \infty)$ and any $x \in \mathds{Z}$, 
  $$
  \lim_{s \to t^-} V^{-1}(x, s) = V^{-1}(x + 1, t) \qquad \text{if and only if} \qquad  \lim_{s \to t^-} V^{-1}(x +1, s) = V^{-1}(x, t).
  $$
  \end{enumerate}
  We think of a swap function as a collection of particle trajectories $\{V(x, \cdot) : x \in \mathds{Z}\}$.
Condition (iv) guarantees that the only way that a particle at position $x$ can move up at time $t$ is if the particle at position $x+1$ moves down. That is, particles move by swapping with their neighbours.   

\medskip
  
  Let $\mathcal{A}$ be the space of swap functions endowed with the following topology. A sequence of swap functions $V_n \to V$ if each of the cadlag paths $V_n(x, \cdot) \to V(x, \cdot)$ and $V^{-1}_n(x, \cdot) \to V^{-1}(x, \cdot)$. Convergence of cadlag paths is convergence in the Skorokhod topology. We refer to a random swap function as a \textbf{swap process}.
  
  \medskip
  
For a swap function $V$ and a time $t \in (0, \infty)$, define
$$
V(\cdot, t, s) = V(V^{-1}(t, \cdot), t + s).
$$
The function $V(\cdot, t, s)$ is the increment of $V$ from time $t$ to time $t + s$. 

\medskip

Now for $i \in \{1, \dots n\}$, define 
$$
r_n(i) = \sqrt{1 - (2i/n - 1)^2},
$$
and consider the shifted, time-scaled swap process 
$$
U_{n}^{i}(x, s) = \sigma^n \left(i + x, \frac{ns}{r_n(i)} \right) - i.
$$
To ensure that $U_{n}^{i}$ fits the definition of a swap process, we can extend it to a random function from $\mathds{Z} \times [0, \infty) \to \mathds{Z}$ by letting $U_{n}^{i}$  be constant after time $\frac{n-1}{2r_n(i)}$, and with the convention that  $U_{n}^{i}(x, s)= x$ whenever $x \notin \{1 -i, \mathellipsis n - i\}$. In the swap processes $U_{n}^{i}$, all particles are labelled by their initial positions.
The following is shown in \cite{angel2017local}, and also essentially in \cite{gorin2017}. 

\begin{theorem}
\label{T:local}
There exists a swap process $U$ such that the following holds. For any $\alpha \in (-1, 1)$, and any sequence of integers $\{i_n\}_{n \in \mathbb{N}}$ such that $i_n/n \to (\alpha + 1)/2$, we have that
  $$
U_{n}^{i_n}  \stackrel{d}{\to} U \qquad \;\text{as}\; \; n \to \infty.
  $$
The swap process $U$ has the following properties:
\begin{enumerate}[nosep,label=(\roman*)]
 \item $U$ is stationary and mixing of all orders with respect to the spatial shift $\tau U(x, t) = U(x + 1, t) - 1$.
  \item $U$ has stationary increments in time: for any $t \ge 0$, the
  process $U(\cdot, t, s)_{s\ge 0}$ has the same law
  as $U(\cdot,s)_{s\geq 0}$.
  \item $U$ is symmetric: $U(\cdot, \cdot) \stackrel{d}{=} - \; U(- \; \cdot, \cdot)$.
  \item For any $t \in [0, \infty)$, $\mathbb{P}($There exists $x \in \mathds{Z}$ such that $U(x, t) \neq \lim_{s \to t^-} U(x, t)) = 0$.
  \end{enumerate}
  \smallskip
  
  Moreover, for any sequence of times $\{t_n : n \in \mathbb{N}\}$ such that $\frac{n-1}{2r_n(i_n)} - t_n \to \infty$ as $n \to \infty$,
  $$
U^{i_n}_n (\cdot, t_n, \cdot) \stackrel{d}{\to} U \qquad \;\text{as}\; n \to \infty.
$$ 
\end{theorem}

Now, for a swap function $V$, let $W(V, t)$ be the number of times that the particles at locations $0$ and $1$ swap in the interval $[0, t]$. That is,
$$
W(V, t) = \card {\left\{s \in (0, t] : 
  \lim_{r \to t^-} V^{-1}(0, r) = V^{-1}(1, s)  \right\}}.
 $$

As a by-product of the proof of convergence of $U_n^{i_n}$ to $U$, the authors of \cite{angel2017local} also found the expectation of $W(U, t)$.

\begin{theorem}[Proposition 7.10, \cite{angel2017local}]
\label{T:swaps}
Let $\alpha \in (-1, 1)$, and let $\{i_n\}_{n \in \mathbb{N}}$ be any sequence of integers converging to $(\alpha + 1)/2.$ Then for any $t \in [0, \infty)$, we have that
$$
\mathbb{E} W(U_n^{i_n}, t) \to \mathbb{E} W(U, t) = \frac{4t}{\pi}.
$$
\end{theorem}

Now let $Q(V, t)$ be the number of swaps that particle $0$ makes by time $t$ in a swap function $V$. That is,
$$
Q(V, t) = \card {\Big\{ s \in (0, t] : \lim_{r \to s^-} V(0, r) \neq V(0, s) \Big\}}.
$$
Dauvergne and Vir\'ag \cite{dauvergne1} used a stationarity argument to prove an analogous result to Theorem \ref{T:swaps} for $Q(u, t)$. 
\begin{theorem}[Lemma 3.2, \cite{dauvergne1}]
\label{T:swaps-2}
For any $t \in [0, \infty)$, we have that
$$
\mathbb{E} Q(U, t) = \frac{8t}{\pi}.
$$
\end{theorem}


The fact that $U$ is stationary in both time and space implies that the point process of swaps of a given particle $x$ in $U$ is also stationary. This realization combined with the ergodic theorem allows us to conclude that all particles in $U$ have asymptotic speeds. This observation was used in \cite{dauvergne1} to prove results about the relationship between the local and global limit. We use their results as a starting point in our proofs.

\begin{theorem}[Theorem 1.7, \cite{dauvergne1}]
 \label{T:local-speeds}
 For every $x \in \mathds{Z}$, the limit
$$
S(x) = \lim_{t \to \infty} \frac{U(x, t) - U(x, 0)}{t}  \qquad \text{exists } \text{almost surely}.
$$
$S(x)$ is a symmetric random variable with distribution $\mu$ independent of $x$. The support of $\mu$ is contained in the interval $[-\pi, \pi]$.
\end{theorem}
We refer to $\mu$ as the {\bf local speed distribution}. Dauvergne and Vir\'ag \cite{dauvergne1} also used Theorem \ref{T:local-speeds} to find limiting swap rates in $U$. Define 
\begin{equation*}
D_\mu^+(c) = \int(y - c)^+d \mu(y), \;\;\;\; D_\mu^-(c) = \int(y - c)^-d \mu(y) \;\;\;\;\;\text{and}\; \;\;\;\; D_\mu(c) = \int|y - c|d \mu(y).
\end{equation*}

\begin{theorem}[Theorem 1.8, \cite{dauvergne1}]
\label{T:particle-swaps} Let $S(0)$ be the asymptotic speed of particle $0$ in $U$. For any $x \in \mathds{Z}$, 
$$
\frac{Q(U, t)}t \to D_\mu(S(0)) \qquad \text{ almost surely and in } L^1.
$$
In particular, the random variables $Q(U, t)/t$ are uniformly integrable and $\mathbb{E} D_\mu(S(0)) = 8/\pi$.
\end{theorem}
We also need an analogous result regarding crossings of lines in the local limit. Let $L(t) = ct + d$ be a line of constant slope $c$. For a swap function $V$, define
$$
C^+(V, L, t) = \card{\Big\{x \in \mathds{Z}: V(x, 0) \le L(0), V(x, t) > L(t)\Big\}}.
$$
$C^+(V, L, t)$ is the total number of particles that are below $L$ at time $0$ and above $L$ at time $t$. We symmetrically define $C^-(V, L, t)$ as the total number of particles that are above $L$ at time $0$ and below $L$ at time $t$, and let $C(V, L, t) = C^-(V, L, t) + C^+(V, L, t)$.

\begin{theorem}[Theorem 5.7, \cite{dauvergne1}]
 \label{T:line-rate}
Let $L(t) = ct + d$. Then almost surely and in $L^1$, we have that
\begin{align*}
\lim_{t \to \infty} \frac{C^+(U, L, t)}t = D_\mu^+(c), \qquad \lim_{t \to \infty} \frac{C^-(U, L, t)}t = D_\mu^-(c), \qquad \lim_{t \to \infty}  \frac{C(U, L, t)}t =  D_\mu(c).
\end{align*}
\end{theorem}

We also record here a few basic facts about the functions $D_\mu$ and $D_\mu^+$ that will be used throughout the paper. These properties can be proven using basic facts about integrals, and the fact that $\mu$ is symmetric and supported in $[-\pi, \pi]$.

\begin{lemma}
\label{L:convex}
\begin{enumerate}[nosep, label=(\roman*)]
\item Both $D_\mu$ and $D_\mu^+$ are convex, $1$-Lipschitz functions.
\item For all $x$, we have that $D_\mu^+(x) \le D_\mu(x) \le |x| \vee \pi$.
\item Suppose that $L(t) = at + b$ is tangent to either $D_\mu$ or $D_\mu^+$. Then $b \in [0, \pi]$.
\item $D_\mu$ is a symmetric function, and hence minimized at $0$.
\item $D_\mu^+(\pi) = 0$, and $D_\mu(\pm \pi) = \pi$.
\end{enumerate}

\end{lemma}

\bigskip

\noindent \textbf{Subsequential limits of $Y_n$.} \qquad Recall that $Y_n$ is the trajectory random variable of $\sigma^n$. We record here the main result of \cite{dauvergne1} regarding subsequential limits of $Y_n$. Here and throughout the paper, the phrase ``subsequential limit of $Y_n$'' always refers to a subsequential limit of $Y_n$ in distribution. 

\medskip

We say that a path $y \in \mathcal{D}$ is $g(y)$-Lipschitz if $y$ is absolutely continuous and if for almost every $t$, $|y'(t)| \le g(y(t))$. 
\begin{theorem}[Theorem 1.4, \cite{dauvergne1}]
\label{T:bounded-speed}
\begin{enumerate}[nosep, label=(\roman*)]
\item The sequence $\{Y_n\}$ is precompact in distribution.
\item Suppose that $Y$ is a subsequential limit of $Y_n$ (in distribution). Then
$$
\mathbb{P}\bigg(Y \text{ is } \pi\sqrt{1-y^2}\text{-Lipschitz},\;Y(0) = -Y(1) \bigg) = 1.
$$
Moreover, $Y(t)$ is uniformly distributed on $[-1, 1]$ for every $t$.
\end{enumerate}
\end{theorem}

In addition, we observe here that any subsequential limit $Y$ of $Y_n$ inherits certain symmetries from $\sigma^n$.

\begin{prop}
\label{P:Y-symmetries} Let $Y$ be any subsequential limit of $Y_n$.
\smallskip
\begin{enumerate}[nosep, label=(\roman*)]
\item Define $Y_t \in \mathcal{D}$ by 
$$
Y_t(s) = \begin{cases}
\;\;Y(s + t), \quad &s \le t.\\
\;\; -Y(s + t - 1), \quad &s > t.
\end{cases}
$$
For any $t \in [0, 1]$, we have that $Y_t \stackrel{d}{=} Y$.
\item $Y \stackrel{d}{=} - Y$.
\item Define $Z \in \mathcal{D}$ by $Z(t) = Y(1 - t)$. Then $Z \stackrel{d}{=} Y$.
\end{enumerate}
\end{prop}

\begin{proof}
Property (i) follows from time stationarity of random sorting networks (Theorem \ref{T:time-stat}). Properties (ii) and (iii) follow from the corresponding properties of the swap sequence $(K^n_1, \dots, K^n_N)$ of $\sigma^n$:
\[
\left(K^n_1, K^n_2, \mathellipsis K^n_N\right) \stackrel{d}{=} \left(n - K^n_1, n - K^n_2,  \mathellipsis, n - K^n_{N}\right) \stackrel{d}{=} \left(K^n_N, \mathellipsis, K^n_2, K^n_1\right). \qquad \qedhere
\]
\end{proof}


\section{Particle flux for Lipschitz paths}
\label{S:flux-intro}

In this section, we introduce particle flux and prove that it measures the amount of particles that cross a line in a typical sorting network. Let $\text{\fontfamily{ppl}\selectfont Lip}$ be the set of Lipschitz paths from $[0,1] \to [-1, 1]$ (we will use this notation throughout the paper).
Define the \textbf{local speed} of a function $h$ at time $t$ by
$$
s(t) = \frac{d(\arcsin(h))}{dt} = \begin{cases}
\frac{h'(t)}{\sqrt{1 - h^2(t)}}, \quad &|h(t)| < 1 \\
0, \quad &|h(t)| = 1.
\end{cases} 
$$
When $h \in \text{\fontfamily{ppl}\selectfont Lip}$, the local speed $s(t)$ exists for almost every time $t$.
Recalling the definition of $D_\mu$ from the previous section, we then define the {\bf particle flux} of $h$ over a set $A$ by
\begin{equation}
\label{E:path-flux}
J(h; A) = \frac{1}2 \int_A D_\mu(s(t)) \sqrt{1-h^2(t)}dt.
\end{equation}
We define $J(h) = J(h; [0, 1])$. Note that $J(h) < \infty$ for any Lipschitz function $h$. This follows from Lemma \ref{L:convex} (ii), which implies that 
$$
D_\mu(s(t)) \sqrt{1-h^2(t)} \le [\pi \vee |s(t)|] \sqrt{1-h^2(t)} \le \pi \vee |h'(t)|.
$$
We will also consider {\bf positive particle flux} $J^+(h; A)$ and {\bf negative particle flux} $J^-(h; A)$ defined by
$$
J^+(h; A) = \frac{1}2 \int_A D^+_\mu(s(t)) \sqrt{1-h^2(t)}dt, \qquad J^-(h; A) = \frac{1}2 \int_A D^-_\mu(s(t)) \sqrt{1-h^2(t)}dt.
$$ 
Again, we let $J^+(h) = J^+(h; [0, 1])$ and $J^-(h) = J^-(h; [0, 1])$.

\medskip

We now connect flux to random sorting networks.
If a random sorting network resembles the local limit in a local window around the global space-time position $(t, h(t))$, then by Theorem \ref{T:line-rate}, the number of distinct particles that cross $h$ in this window should be proportional to 
$$
D_\mu\left(\frac{h'(t)}{\sqrt{1 -h^2(t)}}\right) \sqrt{1-h^2(t)}.
$$
The scaling factors of $\sqrt{1-h^2}$ come from the semicircle rescaling of time in the local limit away from the center.

\medskip

Therefore in a typical large-$n$ sorting network, where most local windows resemble the local limit, $J(h)$ should be proportional to the number of particles that cross the line $h$, counting multiple crossings for a given particle if and only if the crossings happen at globally distinguishable locations. 

\medskip

The factor of $1/2$ is to account for the difference between the ${n \choose 2}^{-1}$ scaling in the global limit and the $n^{-1}$ scaling in the local limit. Similarly, $J^+(h)$ and $J^-(h)$ should be proportional to the number of upcrossings and downcrossings of the line $h$ in a large $n$ sorting network. 

\medskip

Now let $\text{\fontfamily{ppl}\selectfont Lip}_r$ be the set of paths $h \in \text{\fontfamily{ppl}\selectfont Lip}$ with $h(0) = -h(1)$. If $h \in \text{\fontfamily{ppl}\selectfont Lip}_{r}$, then in any sorting network, every particle must cross $h$ at least once unless the particle starts at $h(0)$. Therefore $J(h)$ should be bounded below for such paths. Since any two particles in a sorting network cross each other exactly once, $J(h)$ should achieve this lower bound when $h$ is a trajectory limit. 

\medskip
 
The next theorem makes rigorous this intuition behind the definition of particle flux. To state the theorem, for a random variable $Y \in \mathcal{D}$ and a path $h \in \text{\fontfamily{ppl}\selectfont Lip}$, we define
$$
P^+_Y(h; [a, b]) = \mathbb{P} \Big( \exists s < t \in [a, b] \text{ such that } Y(s) < h(s) \;\text{and}\; Y(t) > h(t) \Big).
$$
In other words, $P^+_Y(h; [a, b])$ is the probability that $Y$ upcrosses $h$ in the interval $[a, b]$. We similarly define the downcrossing probability
$$
P^-_Y(h; [a, b]) = \mathbb{P} \Big( \exists s < t \in [a, b] \text{ such that } Y(s) > h(s) \;\text{and}\; Y(t) < h(t) \Big).
$$

\begin{theorem} 
\label{T:particle-crossings} Let $Y$ be any subsequential limit of $Y_n$.
\begin{enumerate}[label=(\roman*)]
\item Let $h \in \text{\fontfamily{ppl}\selectfont Lip}$ and $[a, b] \subset [0, 1]$. Then $P^+_Y(h; [a, b]) \le J^+(h; \; [a, b])$ and $P^-_Y(h; [a, b]) \le J^-(h; [a, b])$.
\item Let $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$. Then $J(h) \ge 1$.
\item $\mathbb{P}(J(Y) = 1) = 1$.
\end{enumerate}
\end{theorem}

To prove this theorem, we first need two key lemmas about convergence to the local limit.
Recall that $C^+(V, L, t)$ is the number of upcrossings of a line $L(s) = cs + d$ in the interval $[0, t]$ in a swap function $V$. Recall also the definition of $U_n^{i}$ from Section \ref{S:local}.\begin{lemma}
\label{L:L1-swap-rate}
Let $\alpha \in (-1, 1)$, and suppose that $\{i_n\}_{n \in \mathbb{N}}$ is a sequence of integers such that $i_n/n \to (1 + \alpha)/2$. Let $X$ be a uniform random variable on $[0, 1]$ that is independent of $U, U_n^{i_n}$, let $\{c_n\}_{n \in \mathbb{N}}$ be a sequence of real numbers converging to $c \in \mathbb{R}$, and let $\{d_n\}_{n \in \mathbb{N}}$ be a sequence of real numbers in $[0, 1]$. Define 
$$
L_n(s) = c_ns + d_n + X, \quad \text{and} \quad L(s) = cs + X.
$$
Then for any time $t \in (0, \infty)$, we have that 
\begin{enumerate}[label=(\roman*)]
\item $C^+(U_n^{i_n}, L_n, t) \stackrel{d}{\to} C^+(U, L, t)$ as $n \to \infty$.
\item $\mathbb{E} C^+(U_n^{i_n}, L_n, t) \to \mathbb{E} C^+(U, L, t)$ as $n \to \infty$.
\item $\mathbb{E} C^+(U_n^{i_n}, L_n, t) \le 3t + |c_n t| + 2$ for all $n$.
\end{enumerate}
\end{lemma}


\begin{proof} 
First note that $d_n + X \Mod 1 \stackrel{d}{=} X$ for all $n$. Therefore by the stationarity of $U$ with respect to integer-valued spatial shifts (Theorem \ref{T:local} (i)), we have that
\begin{equation}
\label{E:dn-shift}
C^+(U^{i_n}_n, L_n, t) \stackrel{d}{\to} C^+(U, L, t) \quad \text{if and only if} \quad C^+(U^{i_n}_n, c_ns + X, t) \stackrel{d}{\to} C^+(U, L, t).
\end{equation}
Now if $V_n \in \mathcal{A}$ is a sequence of swap functions converging to swap function $V$, then
$$
C^+(V_n, c_ns + X, t) \to C^+(V, L, t)
$$ 
unless $V$ either has a swap at time $t$, or $ct + X \in \mathds{Z}$. By Theorem \ref{T:local} (iv), for any time $t$, the event where $U$ has a swap at time $t$ has probability $0$. Moreover, the probability that $ct + X \in \mathds{Z}$ is also $0$. Therefore 
$$
C^+(U^{i_n}_n, c_ns + X, t) \stackrel{d}{\to} C^+(U, L, t),
$$
and hence (i) follows by statement \eqref{E:dn-shift}. Now recall that $W(V, t)$ is the number of swaps at location $0$ in a swap function $V$ in the interval $[0, t]$. For any swap function $V$ and any line $H(s) = as+b$, we have that
\begin{equation}
\label{E:C-W}
C^+(V, H, t) \le W(V, t) + |a t| + 2.
\end{equation}
To see why this is true, observe that every particle $x$ with $x \le 0$ and $V(x, t) > 1$ must move from position $0$ to position $1$ at least once in the interval $[0, t]$, therefore contributing to $W(V, t)$. Every particle $x$ that upcrosses $H$ in the interval $[0, t]$ fits this description, unless
$$
at \le V(x, t) \le 1.
$$
There are at most $|a t| + 2$ such values of $x$, proving \eqref{E:C-W}.

\medskip

 Now again since $U$ has no swap at time $t$ almost surely (Theorem \ref{T:local} (iv)), $W(U_n^{i_n}, t) \stackrel{d}{\to} W(U, t)$. Also, by Theorem \ref{T:swaps}, we have that
$$
\mathbb{E} W(U_n^{i_n}, t) \to \mathbb{E} W(U, t).
$$
Hence, the random variables $W(U_n^{i_n}, t)$ are uniformly integrable (see Proposition 3.12, \cite{kallenberg2006foundations}). Therefore by \eqref{E:C-W} applied to the swap processes $U_n^{i_n}$ and the lines $L_n$, the random variables $C^+(U_n^{i_n}, L_n, t)$ are also uniformly integrable, and hence converge in expectation since they converge in distribution (again by Proposition 3.12, \cite{kallenberg2006foundations}).

\medskip

Finally, the bound on the probability distribution for the first swap location in a random sorting network (Theorem \ref{T:dist-k1}) and time stationarity (Theorem \ref{T:time-stat}) allows us to conclude that for any $n, i_n$ and $t$, that
$$
\mathbb{E} W(U_n^{i_n}, t) \le 3t,
$$
which in turn proves the bound on the expectation of $C^+(U_n^{i_n}, L_n, t)$ via \eqref{E:C-W}.
\end{proof}

This next lemma is similar to Lemma \ref{L:L1-swap-rate}, but deals with particle speeds rather than upcrossing rates.
For a swap function $V$ we define
$$
S(V, t) = \frac{V(0, t)}{t}.
$$
This is the average speed of particle $0$ in the interval $[0, t]$. 

\begin{lemma}
\label{L:L1-flux}
Let $\{U_n^{i} : i \in \{i, \dots n\}, n \in \mathbb{N}\}$ be the array of locally scaled random sorting networks defined in Section \ref{S:local}, and let $U$ be the local limit. Let $I$ be a uniform random variable on $(-1, 1)$, independent of everything else.
For each $n$, define $I_n = \ceil{n(I + 1)/2}$, let 
$$
R_n = \sqrt{1 - \left[\frac{2I_n}n - 1\right]^2}, \qquad \;\text{and}\; \qquad R = \sqrt{1 - I^2}.
$$
Then the following statements hold.
\begin{enumerate}[label=(\roman*)]
\item For any $t \in (0, \infty)$, we have that $R_nD_\mu(S(U_n^{I_n}, R_nt)) \xrightarrow[n \to \infty]{d} RD_\mu(S(U, R t))$.
\item For any $t \in (0, \infty)$, we have that $\mathbb{E} R_n D_\mu(S(U_n^{I_n}, R_nt)) \xrightarrow[n \to \infty]{} \mathbb{E} R D_\mu \left(S(U, R t) \right)$.
\item $\mathbb{E} R D_\mu \left(S(U, R t) \right) \to 2$ as $t \to \infty$.
\end{enumerate}
\end{lemma}

\begin{proof}
Fix $t \in (0, \infty)$. If $V_n$ is a sequence of swap functions converging to a swap function $V$ with no swap at time $t$, then $D_\mu(S(V_n, t)) \to D_\mu(S(V, t))$. Now condition on $I$. This fixes $R_n, I_n,$ and $R$. For any fixed time $t \in (0, \infty)$, the process $U$ has no swap at time $t$ almost surely (Theorem \ref{T:local} (iv)). Therefore for any bounded continuous test function $f$, we have that
$$
\mathbb{E} \bigg[f\Big(R_nD_\mu(S(U_n^{I_n}, R_nt))\Big) \; \Big | \; I \bigg] \to \mathbb{E} \bigg[f\Big(RD_\mu(S(U,Rt))\Big) \; \Big| \; I \bigg].
$$
Taking expectations proves the distributional convergence in (i).
Now, by Lemma \ref{L:convex}, we have that $D_\mu(s) \le |s| + \pi$ for any $s \in \mathbb{R}$. Recalling that $Q(V, t)$ is the number of swaps made by particle $0$ in the interval $[0, t]$ in a swap process $V$, this implies that
\begin{equation}
\label{E:Q-bd}
R_n\left(D_\mu(S_{I_n}(U_n^{I_n}, R_nt))\right)\le R_n(|S_{I_n}(U_n^{I_n}, R_nt)| + \pi) \le \frac{Q(U_n^{I_n}, R_nt)}{t} + \pi.
\end{equation}

Now, we similarly have that $Q(U_n^{I_n}, R_nt) \stackrel{d}{\to} Q(U, Rt)$, again since $U$ has no swap at time $t$ almost surely. Moreover, 
$$
\mathbb{E} Q(U_n^{I_n}, R_nt) = \frac{2\floor{nt}}n,
$$
since this expectation is simply the expected number of swaps made by a random particle in a sorting network after $2nt$ steps. By Theorem \ref{T:swaps-2}, we also have that
$$
\mathbb{E} Q(U, Rt) = \frac{4t}{\pi} \int_{-1}^1 \sqrt{1 - x^2}dx = 2t,
$$
so $\mathbb{E} Q(U_n^{I_n}, R_nt) \to \mathbb{E} Q(U, Rt)$. Again, by Proposition 3.12 from \cite{kallenberg2006foundations}, this implies that the random variables $\{Q(U_n^{I_n}, R_nt)\}$ are uniformly integrable, and hence so are the random variables $R_n D_\mu(S(U_n^{I_n}, R_nt))$. Since these random variables converge in distribution to $R D_\mu(S(U, R t))$, they must also converge in expectation, proving (ii).

\medskip

Now we prove (iii). First, Theorem \ref{T:local-speeds} implies that 
\begin{equation}
\label{E:D-mu-lim}
R D_\mu(S(U, Rt)) \xrightarrow[t \to \infty]{} R D_\mu(S(0)) \quad \text{almost surely},
\end{equation}
where $S(0)$ is the speed of particle $0$ in $U$. Analogously to \eqref{E:Q-bd}, we also have that
$$
R D_\mu(S(U, R t)) \le \frac{Q(U, t)} {t} + \pi.
$$
By Theorem \ref{T:particle-swaps}, the right hand side above is uniformly integrable, and hence so is the left hand side. Therefore since $RD_\mu(S(U, R t))$ has an almost sure limit by \eqref{E:D-mu-lim}, it also converges in expectation. Finally, Theorem \ref{T:particle-swaps} implies that
\begin{equation*}
\mathbb{E} R D_\mu(S(0)) = 2. \qedhere
\end{equation*}
\end{proof}

To prove Theorem \ref{T:particle-crossings}, we also need two auxiliary results. The first is a simple lemma about uniform convergence of functions. This will be used in the proof of Theorem \ref{T:particle-crossings} (i).

\begin{lemma}
\label{L:mesh}
Let $f_n:[0, 1] \to [-1, 1]$ be a sequence of functions converging uniformly to a continuous function $f$, and let $h:[0, 1] \to [-1, 1]$ be any continuous function. Let 
$$
\{\Pi_n = \{t_{n, 0} =0< t_{n, 1} < \dots < t_{n, m(n)}=1\}\}_{n \in \mathbb{N}}
$$
be a sequence of partitions of $[0, 1]$ such that
$$
\text{mesh}(\Pi_n) = \max_{i \in \{0, 1, \dots, m(n) - 1\}} |t_{n, i + 1} - t_{n, i}| \to 0 \qquad \;\text{as}\; n \to \infty.
$$
Let $[a, b] \subset [0, 1]$. If there exist times $s < t \in [a, b]$ such that $f(s) < h(s)$ and $f(t) > h(t)$, then for all large enough $n$, there exists a time $t_{n, i} \in [a, b]$ such that $f_n(t_{n, i}) \le h(t_{n, i})$ and $f_n(t_{n, i + 1}) > h(t_{n, i + 1})$.
\end{lemma}

\begin{proof}
By the continuity of $f$ and $h$, there exists an $\epsilon > 0$ and disjoint intervals $[s, s_+]$ and $[t_-, t]$ such that $f(r) < h(r) - \epsilon$ for all $r \in[s, s_+]$ and $f(r) > h(r) + \epsilon$ for all $r \in [t_-, t]$. Therefore by uniform convergence, for all large enough $n$ we have that $f_n(r) < h(r) - \epsilon/2$ for all $r \in[s, s_+]$ and $f(r) > h(r) + \epsilon/2$ for all $r \in [t_-, t]$. Now, for large enough $n$ we also have that 
$$
\text{mesh}(\Pi_n) < \min(s_+ - s, t - t_-).
$$
Therefore for such $n$, there must exists $i_1 < i_2 \in \Pi_n \cap [a, b]$ such that  $f_n(t_{n, i_1}) < h(t_{n, i_1})$ and $f_n(t_{n, i_2}) > h(t_{n, i_2})$. Thus for some $i \in \{i_1, \dots i_2 - 1\}$, we must have that $f_n(t_{n, i}) \le h(t_{n, i})$ and $f_n(t_{n, i + 1}) > h(t_{n, i + 1})$.
\end{proof}

To prove part (iii), we also need that $J(\cdot)$ is lower semicontinuous.

\begin{prop}
\label{P:lsc-E}
Let $h_n \in \text{\fontfamily{ppl}\selectfont Lip}$ be a sequence converging uniformly to $h$. Then for any set $A \subset [0, 1]$, we have that
$$
J(h; A) \le \liminf_{n \to \infty} J(h_n ; A).
$$
\end{prop}

\begin{proof}
Note that $J(h; A)$ is of the form $\int_A G(h(t), h'(t))dt$, where $G$ is a continuous, positive function, such that for any fixed value of $a$, $G(a, \cdot)$ is convex. This follows from the convexity of $D_\mu$ (Lemma \ref{L:convex}). Functionals of this form are lower semicontinuous in the uniform norm by Theorem 1.6 from \cite{struwe1996variational}.
\end{proof}
\begin{proof}[Proof of Theorem \ref{T:particle-crossings}.]
\noindent \textbf{Proof of (i):} \qquad Note first that by the symmetry of $Y$ (Proposition \ref{P:Y-symmetries}), that
$$
\mathbb{P}^-_Y(h; [a, b]) = \mathbb{P}^+_{-Y}(-h; [a, b]) = \mathbb{P}^+_{Y}(-h; [a, b])
$$
for any path $h \in \text{\fontfamily{ppl}\selectfont Lip}$ and any interval $[a, b]$. Also, $J^+(-h) = J^-(h)$ by the symmetry of $\mu$ (Theorem \ref{T:local-speeds}). Therefore the assertion that
$
\mathbb{P}^-_Y(h; [a, b]) \le J^-(h)
$
is equivalent to the assertion that $\mathbb{P}^+_Y(-h; [a, b]) \le J^+(-h)$, and so to prove Theorem \ref{T:particle-crossings} (i), it suffices to prove that $\mathbb{P}^+_Y(h; [a, b]) \le J^+(h)$ for every path $h \in \text{\fontfamily{ppl}\selectfont Lip}$.

\medskip
We first prove this for $h \in \text{\fontfamily{ppl}\selectfont Lip}$ with range in the open interval $(-1, 1)$. Let $t \in \mathds{Z} \cap (0, \infty)$, and for $s \in [0, 1]$, define 
$$
s_{n, t} = \frac{2t}{n-1} \floor{\frac{(n-1)s}{2t}}, \quad \;\text{and}\; \quad s^+_{n, t} = \min\left(s_{n, t} + \frac{2t}{n-1}, 1\right).
$$
Let $X$ be a uniform random variable on $[0, 1]$, independent of the sequence $\{Y_n\}$, and let
$$
A^+_{n, t, s} = \left\{ Y_n(s_{n, t}) < h(s_{n, t}) + \frac{2X}n \;\; \;\text{and}\; \;\; Y_n(s^+_{n, t}) \ge h(s^+_{n, t}) + \frac{2X}n \right\}.
$$
Now observe that when $s < b_{n, t}$, then $s_{n, t}^+ \le b_{n, t}$, and so $s^+_{n, t} =  s_{n, t} + 2t/(n-1)$. Therefore for such $s$, time stationarity of random sorting networks (Theorem \ref{T:time-stat}) implies that $A^+_{n, t, s}$ has the same probability as the event
$$
\left\{ Y_n(0) < h(s_{n, t}) + \frac{2X}n \;\; \;\text{and}\; \;\; Y_n\left(\frac{2t}{n-1}\right) \ge h(s^+_{n, t}) + \frac{2X}n \right\}.
 $$
 Here have used that $t \in \mathds{Z}$ to apply time-stationarity. 
We can then express the upcrossing probability $\mathbb{P}(A^+_{n, t, s})$ in terms of the expected number of upcrossings in the local scaling. For $u \in [-1, 1]$, define
$$
[u]_n = \floor{\frac{n(u + 1)}2}, \quad \;\; \{u\}_n = \frac{n(u + 1)}2 - \floor{\frac{n(u + 1)}2}, \quad \;\; g_n(u) = \sqrt{1 - \left(\frac{2[u]_n}n - 1\right)^2}.
$$
Then for $s < b_{n, t}$, we have that
$$
n \mathbb{P}(A^+_{n, t, s}) = \mathbb{E} C^+\left(U_n^{[h(s_{n, t})]_n}, L_{n, s}, g_n(h(s_{n, t})) t \right),
$$
where
$$
L_{n, s}(r) = \frac{n\left(h(s^+_{n, t}) - h(s_{n, t})\right)}{2tg_n(h(s_{n, t}))} r +  \{h(s_{n, t})\}_n +  X.
$$
Also, let 
$$
L_{s}(r) = \frac{h'(s)}{\sqrt{1 -h^2(s)}}r + X.
$$
Since $h$ is Lipschitz and hence differentiable almost everywhere, Lemma \ref{L:L1-swap-rate} (ii) implies that for almost every $s \in [0, b)$, that
\begin{equation}
\label{E:n-PA}
\lim_{n \to \infty} n\mathbb{P}(A^+_{n, t, s}) = \lim_{n \to \infty} \mathbb{E} C^+\left(U_n^{[h(s_{n, t})]_n}, L_{n, s}, g_n(h(s_{n, t})) t \right) =\mathbb{E} C^+(U, L_s, \sqrt{1 - h^2(s)}t).
\end{equation}
Here we have used that the range of $h$ is in $(-1, 1)$ to ensure convergence. Moreover, there exists a constant $C$ such that 
\begin{equation}
\label{E:finite-bd}
n\mathbb{P}(A^+_{n, t, s}) \le Ct
\end{equation}
for all $n \in \mathbb{N}, s < b_{n, t}$.  This follows from Lemma \ref{L:L1-swap-rate} (iii) and the fact that $h$ is Lipschitz. Now let $Z_{n, t}$ be the number of times $s$ of the form $s_{n, t}$ in the interval $[a_{n, t}, b_{n, t})$ such that the upcrossing event $A^+_{n, t, s}$ occurs. We have that
$$
Z_{n, t} = \int_{a_{n, t}}^{b_{n, t}} \frac{(n-1)\mathbbm{1} (A^+_{n, t, s})}{2t} ds.
$$
Therefore by the bounded convergence theorem, we have that
\begin{equation}
\label{E:cvg-ab-flux}
\mathbb{E} Z_{n, t} = \int_{a_{n, t}}^{b_{n, t}} \frac{(n-1)\mathbb{P} (A^+_{n, t, s})}{2t} ds \xrightarrow[n \to \infty]{} \frac{1}{2} \int_{a}^b \frac{\mathbb{E} C^+(U, L_s, \sqrt{1 - h^2(s)}t)}t ds.
\end{equation}
Now, $\mathbb{E} Z_{n, t}$ bounds $\mathbb{P}(Z_{n, t} \ge 1)$. Therefore for any subsequential limit $Y$ of $Y_n$, Lemma \ref{L:mesh} applied to a subsequence $Y_{n_i} - 2X/n \to Y$ implies that
$$
\lim_{n \to \infty} \mathbb{E} Z_{n, t} \ge \liminf_{i \to \infty} \mathbb{P}(Z_{n_i, t} \ge 1) \ge P^+_Y(h; [a, b]).
$$
The integrand on the right hand side of \eqref{E:cvg-ab-flux} is bounded above by $C$ for all $t \in \mathds{Z} \cap (0, \infty)$ by \eqref{E:n-PA} and \eqref{E:finite-bd}. Therefore by Theorem \ref{T:line-rate} and the bounded convergence theorem, the right hand side of \eqref{E:cvg-ab-flux} converges to 
$$
\frac{1}{2} \int_{a}^b D^+_\mu\left(\frac{h'(s)}{\sqrt{1 - h^2(s)}}\right) \sqrt{1 - h^2(s)} ds = J^+(h; \; [a, b]) \qquad \;\text{as}\; \;\; t \to \infty.
$$
This proves Theorem \ref{T:particle-crossings} (i) for $h \in \text{\fontfamily{ppl}\selectfont Lip}$ with range in $(-1, 1)$. Now we extend this to all $h \in \text{\fontfamily{ppl}\selectfont Lip}$. Define $h_m = h \vee(-1 + 1/m) \wedge (1 - 1/m)$, and let 
 $$
 A_m = \{t \in [0, 1] : h(t) \ge 1- 1/m \text{ or } h(t) \le -1 + 1/m\}.
 $$
 Letting $\mathcal{L}$ be Lebesgue measure on $[0, 1]$, we have that
$$
J^+(h_m ; A_m) = D^+_\mu(0)\sqrt{1 - (1- 1/m)^2}\mathcal{L}(A_m), \qquad \;\text{and}\; \qquad J^+(h_m ; A_m^c) = J(h ; A_m^c).
$$
The flux $J^+(h_m ; A_m) \to 0$ as $m \to \infty$, so
$$
\limsup_{m \to \infty} J^+(h_m; [a, b]) \le J^+(h; [a, b]).
$$
Moreover, $h_m$ converges uniformly to $h$, so if $Y$ upcrosses $h$, then $Y$ will eventually upcross $h_m$. Therefore
$$
\liminf_{m \to \infty} P^+_Y(h_m; [a, b]) \ge P^+_Y(h; [a, b]).
$$
Putting these two inequalities together completes the proof.

\medskip

\noindent \textbf{Proof of (ii):} \qquad 
Let $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$, and let $Y$ be a subsequential limit of $Y_n$. Therefore by Theorem \ref{T:particle-crossings} (i),
\begin{align*}
\mathbb{P} \Big(\exists r < t \in [0, 1] \text{ such that either } Y(r) < h(r), \;Y(t) > h(t),\text{ or } &Y(r) > h(r), \;Y(t) < h(t) \Big) \\
&\le J^+(h) + J^-(h) = J(h).
\end{align*}
The event above holds unless $Y(0) = h(0)$ since $Y(0) = - Y(1)$ almost surely by Theorem \ref{T:bounded-speed}. Since $Y(0)$ is uniformly distributed (Theorem \ref{T:bounded-speed}), $Y(0) \ne h(0)$ with probability one, so the left hand side above is equal to $1$.
\medskip

\noindent \textbf{Proof of (iii):} \qquad
Let $Y$ be any subsequential limit of $Y_n$. Fix $t \in (0, \infty) \cap \mathds{Z}$, and define 
$$
t_{n, j} = \frac{2jt}{n-1} \quad \text{for  } j \in \{0, \mathellipsis \floor{(n-1)/(2t)}\} \quad \text{ and let } t_{n, \floor{\frac{n-1}{2t}} + 1} = 1.
$$
Define $Y_{n, t}$ so that 
$
Y_{n, t}\left(t_{n, j}\right) = Y_{n}\left(t_{n, j}\right)$ for $j \in \{0, \mathellipsis \floor{(n-1)/(2t)} + 1\}$, and so that $Y_{n, t}$ is linear at times in between. By time stationarity of random sorting networks (Theorem \ref{T:time-stat}), we can write
\begin{equation}
\label{E:JY-ineq}
\mathbb{E} J(Y_{n, t}) = \sum_{j=0}^\floor{\frac{n-1}{2t}} \mathbb{E} J(Y_{n, t} ;  \; [t_{n, j}, t_{n, j+1}] ) \le \left(\floor{\frac{n-1}{2t}} + 1\right) \mathbb{E} J(Y_{n, t} ;  \; [0, t_{n, 1}]). 
\end{equation}
Here we have used that $t \in \mathds{Z}$ to apply time stationarity of sorting networks. The final term in the sum above may be slightly smaller than the previous terms since the length of interval may be less than $2t/(n-1)$; this gives rise to the inequality.

\medskip

Now we have that
\begin{align}
\nonumber
J(Y_{n, t} ;  \; [0, t_{n, 1}]) &= \frac{1}{2} \int_0^{t_{n, 1}} \int \left|Y_{n, t}'(0)- \sqrt{1 - Y_{n, t}^2(r)}x\right|d\mu(x) dr \\
\label{E:frac-en}
&=  \frac{1}{2} \int_0^{t_{n, 1}} \int \left|Y_{n, t}'(0)- \sqrt{1 - Y_{n, t}^2(0)}x\right|d\mu(x) dr + \epsilon_n
\end{align}
for some error term $\epsilon_n$. Here we have used that $Y'_{n, t}(r) = Y'_{n, t}(0)$ for all $r \in [0, t_{n, 1}]$ by piecewise linearity of $Y'_{n, t}$. Now if $|Y_{n, t}'(0)| \ge \pi$, then since $\mu$ is symmetric and supported in $[-\pi, \pi]$ (Theorem \ref{T:local-speeds}), the two inner integrals are the same. Therefore $\epsilon_n = 0$ in this case. Also, when 
$
|Y_{n, t}'(0)| < \pi 
$
and $x \in [-\pi, \pi]$, the difference of the integrands is bounded above by
$$
\left|\sqrt{1 - Y_{n, t}^2(0)}x - \sqrt{1 - Y_{n, t}^2(r)}x\right| \le |x|\sqrt{1 - \left(1 - \frac{2\pi t}{n-1}\right)^2} \le \pi\sqrt{\frac{4\pi t}{n-1}}.
$$
Hence $|\epsilon_n| \le k(t/n)^{3/2}$ for some constant $k$. Now, letting $I_n = n(Y_{n, t}(0) + 1)/2$, and using the notation of Lemma \ref{L:L1-flux}, \eqref{E:frac-en} is equal to
$$
\frac{t}{n-1}R_nD_\mu\left(S(U_n^{I_n}, R_n t)\right) + \epsilon_n.
$$
Therefore by Lemma \ref{L:L1-flux} (ii) and the bound on $\epsilon_n$, \eqref{E:JY-ineq} implies that
\begin{align*}
\mathbb{E} J(Y_{n, t}) \le \left(\floor{\frac{n - 1}{2t}} + 1 \right)\left[\frac{t}{n-1} \mathbb{E}  R_nD_\mu\left(S(U_n^{I_n}, R_n t)\right)  + \mathbb{E} \epsilon_n \right]\xrightarrow[n \to \infty]{} \frac{1}2 \mathbb{E} R D_\mu \left(S(U, Rt) \right).
\end{align*}
Letting $t \to \infty$, Lemma \ref{L:L1-flux} (iii) then implies that
\begin{equation}
\label{E:flux-lim}
\lim_{t \to \infty} \lim_{n \to \infty} \mathbb{E} J(Y_{n, t}) \le 1.
\end{equation}
Since subsequential limits of $Y_n$ are Lipschitz by Theorem \ref{T:bounded-speed}, $Y$ is a subsequential limit of $Y_n$ if and only if $Y$ is a subsequential limit of $Y_{n, t}$. Therefore, by Fatou's lemma and the lower semicontinuity of $J(\cdot)$ (Proposition \ref{P:lsc-E}), 
\begin{equation*}
\lim_{n \to \infty} \mathbb{E} J(Y_{n, t}) \ge \mathbb{E} J(Y)
\end{equation*}
for all $t$, so $\mathbb{E} J(Y) \le 1$ by \eqref{E:flux-lim}. Moreover, $J(Y) \ge 1$ almost surely by Theorem \ref{T:particle-crossings} (i), so $J(Y) = 1$ almost surely.
\end{proof}

\section{Characterization of minimal flux paths}
\label{S:unique}

In this section, we show that any subsequential limit $Y$ of $Y_n$ is uniquely determined by its initial position, maximum absolute value, and whether it is initially increasing or decreasing. By Theorem \ref{T:particle-crossings} (iii), it is enough to characterize paths $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with $J(h) = 1$. 

\medskip

Let $\text{\fontfamily{ppl}\selectfont Lip}_0$ be the set of path $h \in \text{\fontfamily{ppl}\selectfont Lip}$ with $h(0) = h(1) = 0$. We can first recognize that to characterize minimal flux paths $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$, it is enough to characterize minimal flux paths $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$. This is a simple consequence of the following simple fact.

\begin{lemma}
\label{L:path-shift}
Let $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$ and $t_0 = \inf \{t : h(t) = 0\}$. Define the path
\begin{equation*}
h_0(t) = \begin{cases}
&h(t + t_0), \qquad \quad \quad t \le 1- t_0 \\
&-h(t + t_0 - 1), \qquad t > 1 - t_0.
\end{cases}
\end{equation*}
Then $J(h) = J(h_0)$. In particular, every path $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with flux $J(h) = 1$ can be shifted to a path $h_0 \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with $J(h_0) = r$.
\end{lemma}

We build up to the following characterization of minimal flux paths in $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$. Define the maximum height $m(h)$ of a continuous path $h:[0, 1] \to [-1, 1]$ by
$$
m(h) = \max \{ |h(t)| : t \in [0, 1] \}.
$$ 
Recall also the definition of $g(y)$-Lipschitz from Section \ref{S:local}.
\begin{theorem}
\label{T:path-set}
For every $m \in [0, 1]$, there exists exactly one $\pi\sqrt{1- h^2}$-Lipschitz path $h_m \in \text{\fontfamily{ppl}\selectfont Lip}_0$ such that $J_m(h_m) = 1$, $m(h_m) = m$, and $h_m \ge 0$. The paths $h_m$ satisfy the following properties.
\smallskip
\begin{enumerate}[nosep,label=(\roman*)]
  \item $h_m(1/2) = m$, and $h_m(t) = h_m(1 - t)$ for $t \in [0, 1/2]$.
  \item $h_m$ is strictly increasing on the interval $[0, 1/2]$.
  \item For any $\ell \in [0, m]$, we have that $h_{\ell}(s) \le h_{m}(s)$ for all $s \in [0, 1]$.
\end{enumerate}
\smallskip
Now define $h_{-m} = -h_m$. If $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$ is a $\pi\sqrt{1 -h^2}$-Lipschitz path with $J_m(h) = 1$ and $m(h) = m$, then either $h = h_m$ or $h = h_{-m}$.
  \end{theorem}
  
Theorem \ref{T:path-set} will be proven as Proposition \ref{P:monotone-traj}, Proposition \ref{P:unique-path}, Corollary \ref{C:sym}, Lemma \ref{L:no-cross}, and Proposition \ref{P:exist-m}. We also state a analogue of Theorem \ref{T:path-set} for minimal flux paths $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$. First, for a path $h \in \text{\fontfamily{ppl}\selectfont Lip}$, define
$$
S(h) = \sup_{t \in [0, 1]} h(t) \qquad \;\text{and}\; \qquad I(h) = \inf_{t \in [0, 1]} h(t).
$$
\begin{theorem}
\label{T:unique-2}
Fix $a \in [-1, 1]$ and $k \in [|a|, 1]$. Then the following hold:

\begin{enumerate}[nosep,label=(\roman*)]
  \item There is exactly one $\pi\sqrt{1 - h^2}$-Lipschitz path $h_{a, k} \in \text{\fontfamily{ppl}\selectfont Lip}_r$ such that $h_{a, k}(0) = - h_{a, k}(1) = a$ and $m(h) = S(h) = k$.
 \item There is a exactly one path  $\pi\sqrt{1 - h^2}$-Lipschitz path $h_{a, -k} \in \text{\fontfamily{ppl}\selectfont Lip}_r$ such that $h_{a, -k}(0) = - h_{a, -k}(1) = a$ and $-m(h) = I(h) = - k$. We have that $h_{a, k}(t) = -h_{a, k}(1 - t)$ for all $t$.
  \item If $a \neq 0$, there is a unique time $t \in (0, 1)$ such that $h_{a, k}(t) = 0$. Moreover, the path
$$
g(s) = 
\begin{cases}
&h(t + s), \qquad s \le 1 - t, \\
&-h(t + s - 1), \qquad s > 1 - t 
\end{cases}
$$
is equal to $h_k$ if $a < 0$ and $h_{-k}$ if $a > 0$.
  \item For $k_1 \le k_2 \in [-1, -|a|] \cup [|a|, 1]$, we have that $h_{a, k_1}(t) \le h_{a, k_2}(t)$ for all $t \in [0, 1]$. 
  \item $h_{a, a} = h_{-a, a}$.
\end{enumerate}
\end{theorem}

All parts of this theorem follow from applying Lemma \ref{L:path-shift} to Theorem \ref{T:path-set}, except part (iv). This will be proven in Lemma \ref{L:no-cross}.


\subsection{Basic bounds on $D_\mu$}
\label{SS:basic-bounds}

In order to work with the functional $J(\cdot)$, we first prove a few basic bounds on $D_\mu$.


\begin{lemma}
\label{L:max-heights}
Let $Y$ be a subsequential limit of $Y_n$. For all $a \in [0, 1]$, we have that
$$
\mathbb{P}(m(Y) > a) \leq \frac{D_\mu(0)\sqrt{1- a^2}}{2}.
$$
\end{lemma}

\begin{proof}
Since $Y(0) = -Y(1)$ and $Y(0)$ is uniformly distributed (Theorem \ref{T:bounded-speed}), the left hand side of this inequality can be written as
\begin{align*}
&\mathbb{P}\bigg(Y(0) < a \;\text{and}\; Y(t) > a \text{ for some } t > 0, 
\;\text{or}\; Y(0) > - a \;\text{and}\; Y(t) < - a \text{ for some } t > 0 \bigg).
\end{align*}
By Theorem \ref{T:particle-crossings}, this is bounded above by $J^+(g_a) + J^-(g_{-a})$, where $g_a$ is the path of constant height $a$. Finally, 
\begin{equation*}
J^+(g_a) + J^-(g_{-a}) = \frac{D_\mu(0)\sqrt{1- a^2}}2. \qedhere
\end{equation*}
\end{proof}
\begin{lemma}
\label{L:expect-2}
$D_\mu(0) = 2.$ That is, if $X$ is a random variable with distribution $\mu$, then $\mathbb{E}|X| = 2$.
\end{lemma}

\begin{proof}  Let $Y$ be a subsequential limit of $Y_n$, and let $\alpha = D_\mu(0)$.  By Lemma \ref{L:max-heights}, we have that
$$
1 = \mathbb{P}(m(Y) > 0) \leq \frac{\alpha}2,
$$
so $\alpha \geq 2$. The equality above follows since $m(Y) \ge |Y(0)|$ and $|Y(0)|$ is uniformly distributed on $[-1, 1]$ (Theorem \ref{T:bounded-speed}). Now, for every height $a$ such that 
$$
a > \sqrt{1 - \left(\frac{2}{\alpha}\right)^2},
$$
Lemma \ref{L:max-heights} guarantees that $m(Y) \le a$ with positive probability. Therefore by Theorem \ref{T:particle-crossings} (iii), for any $\epsilon > 0$ there is a $\pi\sqrt{1-y^2}$-Lipschitz path $h_\epsilon \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with $J(h_\epsilon) = 1$ such that
\begin{equation}
\label{E:m-hep}
m(h_\epsilon) \leq m_\epsilon := \sqrt{1 - \left(\frac{2}{\alpha} - \epsilon\right)^2}.
\end{equation}
Using that $D_\mu$ is minimized at $0$ (Lemma \ref{L:convex}), we have the bound
\begin{equation}
\label{E:flux-he}
\begin{split}
1 = J(h_\epsilon) &\ge \frac{\alpha}2 \int_0^1\sqrt{1-h_\epsilon^2(t)}dt.
\end{split}
\end{equation}
The above integrand is always bounded by $2/\alpha - \epsilon$. Also, since $h_\epsilon(0) = -h_\epsilon(1)$ and $h_\epsilon$ is $\pi$-Lipschitz, the amount of time that $h_\epsilon$ spends in the interval $[-m_\epsilon/2, m_\epsilon/2]$ is at least $m_\epsilon/(2\pi)$. Therefore
$$
\frac{2}{\alpha} \ge \int_0^1\sqrt{1-h_\epsilon^2(t)}dt \ge \left(\frac{2}\alpha - \epsilon\right)\left(1 - \frac{m_\epsilon}{2\pi}\right) + \frac{m_\epsilon}{2\pi}\sqrt{1 - \frac{m_\epsilon^2}{4}}.
$$
Letting $\epsilon \to 0$, we get an inequality in $\alpha$ which implies that $\alpha = 2$.
\end{proof}

\begin{lemma}
\label{L:crude-deviation-bd} There is a sequence $u_n \to 0$ such that 
$$
D_\mu(u_n) \leq \frac{2}{\cos(u_n/2)}.
$$
\end{lemma}

\begin{proof}
For any subsequential limit $Y$ of $Y_n$ and any $\epsilon > 0$, Lemmas \ref{L:max-heights} and \ref{L:expect-2} imply that
$$
\mathbb{P}(m(Y) \le \epsilon) \geq 1 - \sqrt{1-\epsilon^2} > 0.
$$
Also, $Y(0)$ is uniformly distributed by Theorem \ref{T:bounded-speed}, so $\mathbb{P}(m(Y) = 0) = 0$. Therefore by Theorem \ref{T:particle-crossings} (iii), we can find a sequence of positive numbers $\alpha_n \to 0$ and a sequence of paths $h_n \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with  $m(h_n) = \alpha_n$ and $J(h_n) = 1$. 
\medskip

Let $s_n$ be the local speed of $h_n$. The total variation of each of the paths $h_n$ is at least $2\alpha_n$, and hence the average absolute local speed of each $h_n$ is at least $2 \arcsin(\alpha_n)$.  The convexity and symmetry of $D_\mu$ (Lemma \ref{L:convex}) then gives the following bound.
\begin{align*}
J(h_n) &\geq \sqrt{1- \alpha_n^2}\int_0^1  D_\mu(s_n(t)) dt \\
&\geq \sqrt {1-\alpha_n^2} D_\mu(2 \arcsin(\alpha_n)).
\end{align*}
Letting $u_n = 2\arcsin(\alpha_n)$ and rearranging completes the proof.
\end{proof}

\subsection{Monotonicity and uniqueness for minimal flux paths}
\label{SS:monotone}
In this subsection, we prove that minimal flux paths $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$ with a particular maximum height $m(h)$ are unique up to sign, and that they satisfy a particular monotonicity relation.

\medskip

We start with two simple lemmas. The first lemma shows that minimal flux paths minimize flux on every subinterval of $[0, 1]$.
  
\begin{lemma}
\label{L:min-path}
Let $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$ be a path with $J(h) = 1$. For any interval $[a, b] \subset [0, 1]$ and any path $g \in \text{\fontfamily{ppl}\selectfont Lip}$ with $g(a) = h(a)$ and $g(b) = h(b)$, we have that
$$
J(h; [a, b]) \le J(g; [a, b]).
$$
Moreover, if $f \in \text{\fontfamily{ppl}\selectfont Lip}_r$ is another path with $J(f) = 1$, $f(a) = h(a)$, and $f(b) = h(b)$, then we can form a new path $p \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with $J(p) = 1$ by letting
\begin{equation*}
p(t) = \begin{cases}
&f(t) \qquad t \leq t_1, t \ge t_2 \\
&h(t) \qquad t \in [t_1, t_2].
\end{cases}
\end{equation*}
\end{lemma}

\begin{proof}
If there is some $g$ with $J(h; [a, b]) > J(g; [a, b])$, then we can make a new path $p \in \text{\fontfamily{ppl}\selectfont Lip}_r$ which is equal to $g$ on $[a, b]$ and equal to $h$ on $[a, b]^c$. This path $p$ will have $J(p) < 1$, contradicting Theorem \ref{T:particle-crossings} (ii). The second part of the lemma is a consequence of the fact that 
\[
J(p) = J(p ; [a, b]) + J(p ; [a, b]^c). \qquad \qquad \qedhere
\]
\end{proof}

The next lemma uses the bounds established in Section \ref{SS:basic-bounds} to eliminate plateaus in minimal flux paths. 

\begin{lemma}
\label{L:no-plat}
For any height $\alpha \in (0, 1)$ and any interval $[t_1, t_2] \subset [0, 1]$, we have that 
$$
\inf_{h \in \text{\fontfamily{ppl}\selectfont Lip}} \big\{J(h; [t_1, t_2]) : h(t_1) = \alpha, h(t_2) = \alpha \big\} < (t_2-t_1)\sqrt{1-\alpha^2}.
$$
\end{lemma}

\begin{proof}
Without loss of generality, we may assume that  $t_1 = 0$. Let $\{u_n\}_{n \in \mathbb{N}}$ be as in Lemma \ref{L:crude-deviation-bd}, and define a sequence of paths $h_n \in \text{\fontfamily{ppl}\selectfont Lip}$ by letting

\begin{equation*}
h_n(t) = \begin{cases} 
\;\sin (\arcsin(\alpha) + u_n t), \qquad &t \leq t_2/2 \\
\;h_n(t_2 - t), \qquad \qquad \quad \;\;\;\;\; &t \in [t_2/2, t_2] \\
\;\alpha, \qquad \qquad \qquad \quad \quad \; &t > t_2.
\end{cases} 
\end{equation*}
Lemma \ref{L:crude-deviation-bd} gives the following bound on the flux of $h_n$ on the interval $[0, t_2]$.
\begin{align*}
J(h_n; [0, t_2]) &\leq \frac{2}{\cos(u_n/2)} \int_0^{t_2/2} \cos(\arcsin(\alpha) + u_n t)dt \\
&= t_2 \sqrt{1-\alpha^2} - \frac{t_2^2 \alpha u_n}{4} + O(u_n^2).
\end{align*}
For small enough $u_n$,  this calculation shows that $J(h_n) < t_2 \sqrt{1-\alpha^2}$.
\end{proof}

We can now prove that minimal flux paths $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$ are unimodal.

\begin{prop}
\label{P:monotone-traj}
Let $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$ be such that $J(h) = 1$, and $m(h) \in (0, 1)$. Then $|h(1/2)| = m(h)$. Moreover, if $h(1/2) = m(h)$, then $h$ is strictly increasing on $[0, 1/2]$ and strictly decreasing on $[1/2, 1]$. If $-h(1/2) = m(h)$, then $h$ is strictly decreasing on $[0, 1/2]$ and strictly increasing on $[1/2, 1]$. 
\end{prop}

\begin{proof}
First consider the case where $h \ge 0$. Let $t_m$ be any time when $h(t_m) = m(h)$. Suppose that there exist times $s_1 < s_2 \in [0, t_m]$ such that $h(s_1) = h(s_2)$, and $h(t) \le h(s_1)$ on the interval $[s_1, s_2]$. Define a new path
$$
r(t) = 
\begin{cases}
h(t), \qquad &t \le s_1, t \ge t_m\\
h(t + (s_2 - s_1)), \qquad &t \in (s_1, t_m - (s_2 -s_1)]\\
m(h) \qquad &t \in [t_m -(s_2 - s_1), t_m].
\end{cases}
$$
The path $r$ replaces the segment of $h$ in the interval $[s_1, s_2]$ with a plateau at height $t_m$ at a later time. This operation cannot increase flux since $D_\mu$ is minimized at $0$ (Lemma \ref{L:convex}), so $r$ must be a minimal flux path in $\text{\fontfamily{ppl}\selectfont Lip}_r$. However 
$$
J(r; [t_m -(s_2 - s_1), t_m]) = (s_2 - s_1)\sqrt{1 - m^2(h)},
$$
which contradicts Lemmas \ref{L:min-path} and \ref{L:no-plat}. Therefore there is no interval $[s_1, s_2] \subset [0, t_m]$ where $h(s_1) = h(s_2)$ and $h(t) \le h(s_1)$ for $t \in [s_1, s_2]$, so $h$ must be strictly increasing on $[0, t_m]$. Similarly, $h$ is strictly decreasing on $[t_m, 1]$. The point $t_m$ is the unique time when $h$ achieves its maximum.

\medskip

We now show that $t_m = 1/2$. Without loss of generality, assume that $t_m \le 1/2$. Define a new path $g(t) = h(1 - t)$. The path $g \in \text{\fontfamily{ppl}\selectfont Lip}_0$ also satisfies $J(g) = 1$ and $m(g) = m$. We have that $g(1/2) = h(1/2)$, and so by Lemma \ref{L:min-path} we can create a path $p \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with $J(p) = 1$ and $m(p) = m$ by letting
\begin{equation*}
p(t) =\begin{cases}
&h(t) \qquad t \leq 1/2\\
&g(t) \qquad t \geq 1/2.
\end{cases}
\end{equation*}
The path $p(t)$ achieves its maximum height at both $t_m$ and $1- t_m$, so we must have that $t_m = 1 - t_m$, and hence $t_m = 1/2$.

\medskip

Now if $h$ is not a strictly non-negative path, then we can create a non-negative path $p(t) = |h(t)|$. By the symmetry of $\mu$ (Theorem \ref{T:local-speeds}), the path $p$ must again have minimal flux, so by the above argument $p$ is strictly increasing on $[0, 1/2]$ and strictly decreasing on $[1/2, 1]$. Therefore either $h = p$ or $h = -p$, completing the proof.
\end{proof}

We can now use Proposition \ref{P:monotone-traj} to prove uniqueness of minimal flux paths with a given maximum height.
\begin{prop}
\label{P:unique-path}
For every $m \in [0, 1]$, there is at most one $\pi\sqrt{1- y^2}$-Lipschitz path $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$ with $J(h) = 1$, $m(h) = m$, and $h \ge 0$. If such a path $h$ exists, then the only other $\pi\sqrt{1- y^2}$-Lipschitz path $g \in \text{\fontfamily{ppl}\selectfont Lip}_0$ with $J(g) = 1$ and $m(g) = m$ is $g = -h$.
\end{prop}

\begin{proof}
First observe that the only $\pi\sqrt{1- y^2}$-Lipschitz paths in $\text{\fontfamily{ppl}\selectfont Lip}_0$ with $m(h) = 1$ are $h = \pm \sin(\pi t)$. Similarly, the only path $h \in \text{\fontfamily{ppl}\selectfont Lip}$ with $m(h) = 0$ is $h = 0$. This proves the proposition for $m \in \{0, 1\}$. Now we assume $m \in (0, 1)$.

\medskip

Suppose that $h, g \in \text{\fontfamily{ppl}\selectfont Lip}_0$ are two distinct non-negative paths with $J(h) = J(g) = 1$ and $m(h) = m(g) = m$. By Proposition \ref{P:monotone-traj}, $h(1/2) = g(1/2) = 1/2$. Without loss of generality, there exists a time $t_2 \in [0, 1/2)$ such that $g(t_2) < h(t_2)$. Let $t_1 < t_2$ be such that $h(t_1) = g(t_2)$. Define the shifted path
\begin{equation*}
g^*(t) = \begin{cases}
\;\;g(t + (t_2 - t_1)), \qquad &t \leq 1 - (t_2 - t_1) \\
\;\;- g(t + (t_2 - t_1) - 1), \qquad& t \geq 1 - (t_2 - t_1).
\end{cases}
\end{equation*}
By Proposition \ref{P:monotone-traj}, we have that 
$$
m = g^*(1/2 - (t_2 - t_1)) > h(1/2 - (t_2 - t_1)) \quad \;\text{and}\; \quad g^*(1/2) < h(1/2) = m.
$$
Therefore there is some time $t_3 \in (1/2 - (t_2 - t_1), 1/2)$ such that $g^*(t_3) = h(t_3)$.  Moreover, $g^* \in \text{\fontfamily{ppl}\selectfont Lip}_r$, $J(g^*)  = 1$, and $g^*(t_1) = h(t_1)$. Therefore by Lemma \ref{L:min-path} we can create a path $p \in \text{\fontfamily{ppl}\selectfont Lip}_0$ with $J(p) = 1$ by letting
\begin{equation*}
p(t) = \begin{cases} 
\;\;h(t), \qquad& t \leq t_1 \;\text{or}\; t \geq t_3 \\
\;\;g^*(t), \qquad& t \in (t_1, t_3).
\end{cases}
\end{equation*}
This new path $p$ is a non-negative minimal flux path in $\text{\fontfamily{ppl}\selectfont Lip}_0$ which does not uniquely achieve its maximum value at $1/2$, contradicting Proposition \ref{P:monotone-traj}.

\medskip

Now, if $g \in \text{\fontfamily{ppl}\selectfont Lip}_0$ is another minimal flux path with $J(g) = 1$ and $m(g) = m$, then $|g|$ is a non-negative minimal flux path in $\text{\fontfamily{ppl}\selectfont Lip}_0$, so $|g| = h$. Since $h(t) > 0$ for $t \in (0, 1)$ by Proposition \ref{P:monotone-traj}, either $g = h$ or $g = -h$.
\end{proof}

Proposition \ref{P:unique-path} gives us the following easy corollary.

\begin{corollary} 
\label{C:sym}
Any $\pi\sqrt{1-h^2}$-Lipschitz path $h \in \text{\fontfamily{ppl}\selectfont Lip}_0$ with $J(h) = 1$ has $h(t) = h(1-t)$ for all $t \in [0, 1]$.
\end{corollary}
\begin{proof}
Without loss of generality, assume that $h \ge 0$. 
Both $h(t)$ and $h(1-t)$ are non-negative minimal flux paths satisfying all the conditions of Proposition \ref{P:unique-path} and with the same maximum height, so by that proposition, $h(t) = h(1-t)$ for all $t$.
\end{proof}

\begin{remark}
\label{R:shift}
Note that the uniqueness proofs in this section automatically give the uniqueness of the paths $h_{a, k}$ in Theorem \ref{T:unique-2}. The fact that uniqueness immediately carries over to shifted paths follows from the strict monotonicity in Proposition \ref{P:monotone-traj}.
\end{remark}


\subsection{Existence of minimal flux paths}
\label{SS:unique}

We now show that for every $m \in [0, 1]$, that there exists a minimal flux path in $\text{\fontfamily{ppl}\selectfont Lip}_0$ with maximum height $m$. This is a consequence of the following proposition, which shows that the inequality in Theorem \ref{T:particle-crossings} is an equality for minimal flux paths. Recall the definition of the upcrossing probability $P^+_Y(h; [a, b])$ and the downcrossing probability $P^-_Y(h; [a, b])$ from Section \ref{S:flux-intro}.

\begin{prop}
\label{P:flux-upcross}
Let $h \in \text{\fontfamily{ppl}\selectfont Lip}_r$ be a path with $J(h) = 1$. Let $Y$ be a subsequential limit of $Y_n$. Then for any interval $[a, b]$, we have that
$$
P^+_Y(h; [a, b])  = J^+(h; [a, b]) \quad \;\text{and}\; \quad P^-_Y(h; [a, b]) = J^-(h; [a, b]).
$$
\end{prop}



\begin{proof}
We will only prove the first equality, as the second one follows by the symmetry of $Y$ and $\mu$ (Theorem \ref{T:local-speeds}, Proposition \ref{P:Y-symmetries}).
Since $Y(0) = -Y(1)$ almost surely by Theorem \ref{T:bounded-speed}, we have that
\begin{equation*}
1 = P^+_Y(h; [0, 1]) + P^-_Y(h; [0, 1])= J^+(h) +J^-(h).
\end{equation*}
Therefore by Theorem \ref{T:particle-crossings} (i), we have that $P^+_Y(h; [0, 1]) = J^+(h)$.
Now we have that
\begin{equation}
\label{E:big-J+}
\begin{split}
P^+_Y(h; [0, 1])& \le
P^+_Y(h; [0, a]) + P^+_Y(h; [a, b]) + P^+_Y(h; [b, 1])\\
 &\le J^+(h; [a, b]^c) + P^+_Y(h; [a, b]) \\
 &\le  J^+(h; [a, b]^c) + J^+(h; [a, b]) = J^+(h).
\end{split}
\end{equation}
To see the first inequality above, note that the union of the three events on the right hand side contains the event on the left hand side minus the set $\{Y(a) = h(a) \;\text{or}\; Y(b) = h(b)\}$. This set has probability $0$ since $Y(t)$ is uniformly distributed for all $t$ (Theorem \ref{T:bounded-speed}). The second and third inequalities follow from Theorem \ref{T:particle-crossings} (i). 

\medskip

Since $P^+_Y(h; [0, 1]) = J^+(h)$, all the inequalities in \eqref{E:big-J+} must in fact be equalities, proving the lemma.
\end{proof}

To apply the above proposition in order to prove the existence of minimal flux paths at every height, we need three lemmas. The first lemma shows that any two minimal flux paths in $\text{\fontfamily{ppl}\selectfont Lip}_r$ cannot cross each other more than once. This lemma also proves part (iii) of Theorem \ref{T:path-set}.

\begin{lemma}
\label{L:no-cross}
Suppose that $h, g \in \text{\fontfamily{ppl}\selectfont Lip}_r$ are non-negative paths with $J(h) = J(g) = 1$. Then there cannot exist times $t_0 < t_1 < t_2 \in [0, 1]$ such that $h(t_0) < g(t_0)$, $h(t_1) > g(t_1)$, and $h(t_2) = g(t_2)$.
In particular, if $g(0) = h(0)$, then either $h \le g$ or $g \le h$.
\end{lemma}

\begin{proof}
Suppose there exist times $t_0 < t_1 < t_2 \in [0, 1]$ such that $h(t_0) < g(t_0)$, $h(t_1) > g(t_1)$, and $h(t_2) = g(t_2)$. First, by possibly shifting the paths as in Lemma \ref{L:path-shift}, we may assume that $t_2 = 1$. Hence $g(1) = h(1) = -g(0) = -h(0)$. Now let $s \in (t_0, t_1)$ be such that $h(s) = g(s)$. Letting $a = g(0)$, Remark \ref{R:shift} implies that $h = h_{a, k_1}$ and $g = h_{a, k_2}$ for some $k_1, k_2 \in [-1, -|a|] \cup [|a|, 1]$. Here the paths $h_{a, k}$ are as in Theorem \ref{T:unique-2}. Without loss of generality, assume that $k_2 > k_1$ and that $k_2 \ge a$; the other cases follow by symmetric arguments.

\medskip

Define $p = \max(h, g)$. By Lemma \ref{L:min-path}, $J(p) = 1$. Moreover, $p(0) = -p(1) = a$, and since $k_2 > k_1$ and $k_2 \ge |a|$, we have that $m(p) = S(p) = k_2$. This contradicts the uniqueness of the path $h_{a, k_2}$ established in Remark \ref{R:shift}.
\end{proof}

The next lemma establishes a strong concavity property for minimal flux paths.

\begin{lemma}
\label{L:concave-paths} 
For every $k \in [0, 1]$, the path $\arcsin(h_k)$ is concave. In other words, the local speed $s_k$ of $h_k$ is a non-increasing function of $t$.
\end{lemma}

\begin{proof} 
By the symmetry of $h_k$ (Corollary \ref{C:sym}), it suffices to prove that $\arcsin(h_k)$ is concave on $[0, 1/2]$. Suppose that there exist times $t_1 < t_2 < t_3 \le 1/2$ such that $h(t_2) < p(t_2)$, where $p$ is the path of constant local speed $\close{s} \in [0, \pi]$ with $p(t_1) = h(t_1)$ and $p(t_2) = h(t_2)$. Let 
$$
t_4 = \sup \{t \in [t_1, t_2] : p(t) \le h(t)\} \qquad \;\text{and}\; \qquad t_5 = \inf \{t \in [t_2, t_3] : p(t) \le h(t)\}.
$$
Then $h(t_4) = p(t_4)$ and $h(t_5) = p(t_5)$, and for every $t \in (t_4, t_5)$, we have that $h(t) < p(t)$. By Lemma \ref{L:convex}, we can find a line $L(s) = as + b$ with $b \ge 0$ such that $L(s) \leq D_\mu(s)$ for all $s$, and such that $L(\close{s}) = D_\mu(\close{s}).$ Therefore 
\begin{align*}
J(h; [t_4, t_5]) 
& \geq \int_{t_4}^{t_5} (as(t) + b) \sqrt{1-h^2(t)} dt \\
&\ge \int_{t_4}^{t_5} ah'(t) dt + \int_{t_4}^{t_5} b \sqrt{1-p^2(t)} dt.
\end{align*}
The inequality in the second line follows since $p(t) > h(t) \ge 0$ for all $t \in (t_4, t_5)$. By the fundamental theorem of calculus, since $p(t_4) = h(t_4)$ and $p(t_5) = h(t_5)$, we have
\begin{align*}
\int_{t_4}^{t_5} ah'(t) dt + \int_{t_4}^{t_5} b \sqrt{1-p^2(t)} dt &= 
\int_{t_4}^{t_5} ap'(t) dt + \int_{t_4}^{t_5} b \sqrt{1-p^2(t)} dt \\
& = J(p; [t_4, t_5]).
\end{align*}
We can then create a non-negative  $\pi\sqrt{1 - h^2}$-Lipschitz path $g \in \text{\fontfamily{ppl}\selectfont Lip}_0$ with $m(g) = k$ by letting $g(t) = p(t)$ for $t \in [t_4, t_5]$ and $g(t) = h(t)$ otherwise. $J(g) \le 1$, contradicting the uniqueness of $h_k$. 
\end{proof}

The next lemma is a consequence of the symmetry and concavity of $Y$.

\begin{lemma}
\label{L:Y-sym-type}
Let $Y$ be a subsequential limit of $Y_n$, let $t \in [0, 1)$, and let $[a, b] \subset [0, 1]$. Then $Y'(t)$ exists almost surely and 
\begin{enumerate}[label=(\roman*)]
\item $\mathbb{P}(Y'(t) = 0) = 0.$
\item $\mathbb{P}(Y'(t) > 0 , \; Y(t) \in [a, b]) = \mathbb{P}(Y'(t) < 0 , \; Y(t) \in [a, b]) = 1/2.$
\end{enumerate}
\end{lemma}

\begin{proof}
By time-stationarity of $Y$ (Proposition \ref{P:Y-symmetries} (i)) it suffices to consider $t = 0$. This also implies that
$$
Y'(0) \stackrel{d}{=} Y'(U),
$$
for an independent uniform random variable $U$ on $[0, 1]$. Since $Y$ is Lipschitz and hence almost everywhere differentiable, this proves that $Y'(0)$ exists almost surely. By the concavity and strict monotonicity of minimal flux paths (Lemma \ref{L:concave-paths} and Proposition \ref{P:monotone-traj}), we have that
\begin{equation}
\label{E:prob-Y-0}
\mathbb{P} (Y'(U) = 0) = \mathbb{P}( m(Y) = |Y(U)|) = 0.
\end{equation}
Define $Z(t) = -Y(1- t)$. Since $Y(0) = -Y(1)$ (Theorem \ref{T:bounded-speed}), we have that
$
Z(0) = Y(0)
$
and $Z'(0) = -Y'(0)$. By Proposition \ref{P:Y-symmetries}, $Y \stackrel{d}{=} Z$, so 
$$
\mathbb{P}(Y'(0) > 0 , Y(0) \in [a, b]) = \mathbb{P}(Y'(0) < 0 , Y(0) \in [a, b]).
$$
Putting this together with \eqref{E:prob-Y-0} completes the proof.
\end{proof}

We are now ready to prove the existence of minimal flux paths at every height.

\begin{prop}
\label{P:exist-m}
For every $m \in [0, 1]$, there exists a $\pi\sqrt{1- h^2}$-Lipschitz path $h_m \in \text{\fontfamily{ppl}\selectfont Lip}_0$ such that $m(h_m) = m$, $J(h_m) = 1$, and $h \ge 0$.
\end{prop}


\begin{proof} First observe that we can set $h_0(t) = 0$ and $h_1(t) = \sin (\pi t)$. Both of these paths satisfy all of the above properties. For this, we use that $D_\mu(\pm \pi) = \pi$ and $D_\mu(0) = 2$ (Lemma \ref{L:convex} and Lemma \ref{L:expect-2}).

\medskip

Let $A \subset [0, 1]$ be the set of all $m$ such there is a $\pi\sqrt{1- h^2}$-Lipschitz path $h_m \in \text{\fontfamily{ppl}\selectfont Lip}_0$ with $h_m \ge 0$, $J(h_m) = 1$ and $m(h_m) = m$.  Let $m \in \close{A}$, and let $m_k \in A$ be a sequence converging to $m$. 

\medskip

By the Arzel\`a-Ascoli Theorem, there exists a $\pi\sqrt{1-h^2}$-Lipschitz subsequential limit $h_m \in \text{\fontfamily{ppl}\selectfont Lip}_0$ of $h_{m_k}$ with $m(h_m) = m$ and $h_m \ge 0$. By the lower semicontinuity of $J$ (Proposition \ref{P:lsc-E}), $J(h_m) = 1$. Therefore $A$ is closed.

 \medskip
 
 Now suppose that $A \ne [0, 1]$. Then there exists an interval $(\ell, m) \subset [0, 1]$ such that $m \in A$ and $(\ell, m) \cap A = \emptyset$. Let $k = (\ell + m)/2$, and choose $t_1 < t_2 \in [0, 1/2]$ so that $
 h_m(t_1) = \ell $ and $h_m(t_2) = k.$
We will show that
 \begin{equation}
 \label{E:P=P}
P^+_Y(h_m ; [t_1, t_2]) + P^-_Y(h_m ; [t_1, t_2]) =  P^+_Y\left(k ; [t_1, t_2]\right) + P^-_Y\left(k; [t_1, t_2]\right).
 \end{equation}
 Here we use $k$ as shorthand for the line $h(t) = k$.
 
 \medskip
 
Let $g \in \text{\fontfamily{ppl}\selectfont Lip}_r$ be any minimal flux $\pi\sqrt{1- y^2}$-Lipschitz path. By the monotonicity established in Proposition \ref{P:monotone-traj}, $g$ upcrosses $k$ at most once in the interval $[0, 1]$, and downcrosses $k$ at most once. Moreover, by Lemma \ref{L:no-cross}, $g$ crosses $h_m$ at most once.

\medskip

If $g(t_1) \notin [\ell, k]$, and $g(t_2) \ne k$, then the total number of crossings of the lines $h_m$ and $k$ is even. Hence, either $g$ crosses $h_m$ the same amount of times that $g$ crosses $k$, or $g$ first crosses $k$ going down, and then crosses $k$ going up. This second possibility cannot happen by the monotonicity of $g$ (Proposition \ref{P:monotone-traj}).

\medskip

If $g(t_1) \in (\ell, k)$ and $g'(t_1) < 0$, then again by the monotonicity properties of $g$, $g$ downcrosses $h_m$, does not upcross $h_m$, and neither upcrosses nor downcrosses $k$. 

\medskip

If $g(t_1) \in (\ell, k)$, $g'(t_1) > 0$, and $g(t_2) \ne k$, then $m(g) \ge m$ since $(\ell, m) \cap A = \emptyset$. By the monotonicity of $g$ (Proposition \ref{P:monotone-traj}), $g(t^*) = m(g)$ for some $t^* > t_1$. If $g(t_2) < h(t_2)$, then since both $g$ and $h_m$ are minimal flux paths, the restrictions on crossings imposed by Lemma \ref{L:no-cross} imply that $t^* < t_2$. In this case, $g$ both upcrosses and downcrosses $k$ in the interval $[t_1, t_2]$, and downcrosses but does not upcross $h_m$.

\medskip

 On the other hand, if $g(t_2) > h(t_2)$, then $g$ upcrosses $k$ in the interval $[t_1, t_2]$ but does not downcross $k$, and does not cross $h$ at all in this interval.


\medskip


Now let $Y$ be any subsequential limit of $Y_n$. By Theorem \ref{T:bounded-speed} and Lemma \ref{L:Y-sym-type} (i), we have that
$$
\mathbb{P} \Big(Y(t_1) \notin \{\ell, k\}, Y(t_2) \ne k, Y'(t_1) \text{ exists and is not equal to $0$}\Big) = 1.
$$
Therefore by the above analysis of minimal flux $\pi\sqrt{1- y^2}$-Lipschitz paths $g \in \text{\fontfamily{ppl}\selectfont Lip}_r$, we have that
\begin{align*}
P^+_Y(h_m ; [t_1, t_2]) - P^+_Y\left(k ; [t_1, t_2]\right) &= - \mathbb{P}(Y'(t_1) > 0, Y(t_1) \in [\ell, k]), \quad \;\text{and}\; \\
P^-_Y(h_m ; [t_1, t_2]) - P^-_Y\left(k ; [t_1, t_2]\right) &= \mathbb{P}(Y'(t_1) < 0, Y(t_1) \in [\ell, k]).
\end{align*}
Equation \eqref{E:P=P} then follows from the symmetry established in Lemma \ref{L:Y-sym-type} (ii). Now by Proposition \ref{P:flux-upcross}, the left hand side of \eqref{E:P=P} is equal to $J(h_m ; [t_1, t_2])$, and by Theorem \ref{T:particle-crossings} (i), the right hand side of \eqref{E:P=P} is bounded above by $J\left(k; [t_1, t_2] \right)$.
However, since $D_\mu$ is minimized at $0$ by Lemma \ref{L:convex}, we can also easily calculate that
\begin{align*}
J(h_m ; [t_1, t_2]) &\ge \frac{D_\mu(0)}{2} \int_{t_1}^{t_2} \sqrt{1 - h^2_m(t)}dt \\
&> \frac{D_\mu(0)}2(t_2 - t_1)\sqrt{1 - {k^2}} = J\left(k; [t_1, t_2] \right).
\end{align*}
This is a contradiction, so $A$ must be the whole interval $[0,1]$.
\end{proof}


\section{The derivative distribution}
\label{S:path-speed}

Let
$
\mathbb{P}_{a}(Y \in \cdot)
$
be the regular conditional distribution of $Y$ given that $Y(0) = a$. In this section, we use the structure of minimal flux paths to find $\mathbb{P}_a(Y'(0) \in \cdot)$.

\begin{prop}
\label{P:local-global-deriv}
Let $Y$ be any subsequential limit of $Y_n$. With probability one, $Y'(0)$ exists. Moreover,
\begin{equation}
\label{E:X-Y'}
\mathbb{P}_{a} \left(\frac{Y'(0)}{\sqrt{1 - a^2}} \in \cdot\right) = \mu,
\end{equation}
for Lebesgue-a.e. $a \in (-1, 1)$.
\end{prop}

As an immediate corollary of Proposition \ref{P:local-global-deriv}, we will be able to find the distribution of the maximum height function $m(Y)$ for any subsequential limit $Y$ of $Y_n$. This will allow us to show that the sequence $\{Y_n\}$ has a unique limit point. 

\medskip

Moreover, Proposition \ref{P:local-global-deriv} will be used later to prove an integral transform formula for $\mu$, which will allow us to determine $\mu$, and then $Y$. This is done in Section \ref{S:integral-formula}.


\medskip
 
To prove Proposition \ref{P:local-global-deriv}, we first show that minimal flux paths fill space. For this lemma, we use the notation of Theorem \ref{T:unique-2}.

\begin{lemma}
\label{L:closed-int}
Let $a \in [0, 1]$, and let 
$$
M_a = \left\{h_{a, k} : k \in  [-1, -a] \cup [a, 1] \right\}.
$$
 Then for any $t \in [0, 1]$, we have that
 $$
\{h(t) : h \in M_a\} = [\sin(-\pi t + \arcsin(a)), \sin(\pi t + \arcsin(a))].
$$
\end{lemma}

\begin{proof}
Let $K_{a, t} = \{h(t) : h \in M_a\}.$
$K_{a, t}$ contains the endpoints of the above interval, since $M_a$ contains the two extremal functions $h_{a, -1}(t) = \sin(-\pi t + \arcsin(a))$ and $h_{a, 1}(t) = \sin(\pi t + \arcsin(a))$. Moreover, the ordering on minimal flux paths (Theorem \ref{T:unique-2} (iv)) implies that 
$$
K_{a, t} \subset [\sin(-\pi t + \arcsin(a)), \sin(\pi t + \arcsin(a))].
$$
Now suppose that $c \in (\sin(-\pi t + \arcsin(a)), \sin(\pi t + \arcsin(a)))$. Let 
$$
A = \{h_{a, k} \in M_a : h_{a, k}(t) \le c\} \qquad \;\text{and}\; \qquad B = \{h_{a, k} \in M_a : h_{a, k}(t) > c\}.
$$
Both of these sets are non-empty.
 By Theorem \ref{T:unique-2} (iv), $M_a$ is a linearly ordered subset of the set of $\pi \sqrt{1-h^2}$-Lipschitz paths in $\text{\fontfamily{ppl}\selectfont Lip}$ with the usual partial order on functions. Therefore the sets $A$ and $B$ have a supremum and infimum in $\text{\fontfamily{ppl}\selectfont Lip}$ by the Arzel\`a-Ascoli Theorem. By the lower semicontinuity of flux (Proposition \ref{P:lsc-E}), $\sup A$ and $\inf B$ lie in the set $M_a$. This implies that
 \begin{align*}
  \sup A &= h_{a, k_A}, \quad \text{ where } \quad k_A = \sup \{ k : h_{a, k} \in A\}, \quad \;\text{and}\; \\
   \inf A &= h_{a, k_B}, \quad \text{ where } \quad k_B = \inf \{ k : h_{a, k} \in B\}.
 \end{align*}
Since $A \cup B = M_a$, we either have that $k_A = k_B$, or else $k_A = -a$ and $k_B = a$. Since $h_{a, a} = h_{-a, a}$ by Theorem \ref{T:unique-2} (v), this implies that $h_{a, k_A}(t) = h_{a, k_B}(t) =  c$, so $c \in K_{a, t}$.
\end{proof}

\begin{proof}[Proof of Proposition \ref{P:local-global-deriv}] First note that $Y'(0)$ exists almost surely by Lemma \ref{L:Y-sym-type}.
We now show that for every $q \in (-\pi, \pi)$ and $b < c \in (-1, 1)$, that
\begin{equation}
\label{E:deriv-integrated}
\int_{b}^c \mathbb{E}_{a} (Y'(0) - \sqrt{1 - a^2}q)^+ da = D_\mu^+(q) \int_b^c \sqrt{1-a^2}da.
\end{equation}
Here $\mathbb{E}_a$ is the expectation taken with respect to $\mathbb{P}_a$. Lemma \ref{L:closed-int} guarantees that there exists a time $t_0>0$ such that for any $r \in [0, t_0]$ and $a \in [b, c]$, there is a minimal flux $\pi\sqrt{1-g^2}$-Lipschitz path $g \in \text{\fontfamily{ppl}\selectfont Lip}_r$ with $g(0) = a$ and $g(r) = \sin(\arcsin(a) + qr)$. Call this path $g_{a, r}$. Letting
$$
s_{a, r}(t) = \frac{g'_{a, r}(t)}{\sqrt{1 - g^2_{a, r}(t)}}
$$
be the local speed of $g_{a, r}$, we have that
\begin{align}
\nonumber
J^+(g_{a, r} ; [0, r]) &\ge \frac{1}2 \min_{t \in [0, r]} \sqrt{1 - g^2_{a, r}(t)} \int_0^t D^+_\mu(s_{a, r}(t)) dt \\
\label{E:J-lower}
&\ge  \frac{1}2 \min_{t \in [0, r]} {\sqrt{1 - g^2_{a, r}(t)}} t D^+_\mu(q).
\end{align}
Here the second inequality follows by Jensen's inequality, using that $D^+_\mu$ is convex (Lemma \ref{L:convex}) and that the average local speed of $g_{a, r}$ is $q$. Now letting $g_a(t) = \sin(\arcsin(a) + qt)$, Lemma \ref{L:min-path} implies that
\begin{equation}
\label{E:g-J-h}
J^+(g_{a, r} ; [0,r]) \le J^+(g_a ; [0, r]) \le \frac{1}2  \max_{t \in [0, r]} {\sqrt{1 - g^2_{a}(t)}} r D^+_\mu(q).
\end{equation}
Combining this inequality with \eqref{E:J-lower} implies that
$$
\lim_{r \to 0} \frac{J^+(g_{a, r} ; [0, r])}{r} = \frac{D_\mu^+(q)}2 \sqrt{1 - a^2}.
$$
Moreover, \eqref{E:g-J-h} implies that $J^+(g_{a, r} ; [0, r]) \le rD^+_\mu(q)$ for all $a \in [b, c], r \in [0, t_0]$. Therefore
by the bounded convergence theorem, we have that
\begin{equation}
\label{E:want-1}
\lim_{t \to 0} \int_{b}^c \frac{J^+(g_{a, r} ; [0, r])}r da = \frac{D^+_\mu(q)}2\int_b^c \sqrt{1- a^2}da.
\end{equation}
Now recall that $P^+_Y(g_{a, r} ; [0, r])$ is the probability that $Y$ upcrosses $g_{a, r}$ in the interval $[0, r]$. Since minimal flux paths cross at most once (Lemma \ref{L:no-cross}),
$$ 
\mathbb{P}(Y(0) < a, Y(r) > g_{a, r}(r))= P^+_Y(g_{a, r} ; [0, r]).
$$
Therefore by Proposition \ref{P:flux-upcross}, we can write 
\begin{align}
\nonumber
\int_{b}^c \frac{J^+(g_{a, r} ; [0, r])}r da &= \int_{b}^c \frac{\mathbb{P}(Y(0) < a, Y(r) > g_{a, r}(r))}r da.
\end{align}
Then defining $P_x(a, r) = \mathbb{P}_x(Y(r) > g_{a, r}(r))/r$, the right hand side above is equal to 
\begin{align}
\label{E:J-+-da}
\int_{b}^c \frac{1}2 \int_{-1}^a P_x(a, r)dxda = \frac{1}2 \int_{-1}^c \int_{x \vee b}^c  P_x(a, r) dadx.
\end{align}
Now define $S_{x, y} = \sin(\arcsin(x) + (\pi - q)y)$. Since $Y$ is $\pi\sqrt{1-h^2}$-Lipschitz almost surely (Theorem \ref{T:bounded-speed}), for almost every $x \in [-1, 1]$, we have that
$$
\mathbb{P}_x(Y(t) > g_{a, r}(r)) = 0
$$
whenever $a > S_{x, r}.$ Therefore we can rewrite the right hand side of \eqref{E:J-+-da} as
\begin{equation}
\label{E:want-2}
\begin{split}
\frac{1}2 \left[\int_b^{c} \int_{x}^{S_{x, r}} P_x(a, r) da dx - \int_{S_{c, -r}}^c \int_c^{S_{x, r}} P_x(a, r) da dx + \int_{S_{b, -r}}^b \int_b^{S_{x, r}} P_x(a, r) da dx\right].
\end{split}
\end{equation}
In the second and third terms above, the integrand is of size $O(r^{-1})$, whereas the region of integration is of size $O(r^{2})$. Therefore these terms go to zero as $r \to 0$. Now letting
$$
f(a, x, s, r) = \frac{\sin(\arcsin(a) + sr) - x}r,
$$
and making the substitution $u = f(a, x, q, r)$, we can write the first term of \eqref{E:want-2} as
\begin{equation}
\label{E:b-c}
\begin{split}
\int_b^c \int_{f(x, x, q, r)}^{f(x, x, \pi, r)} \mathbb{P}_x\left(\frac{Y(r) - Y(0)}r > u \right) \frac{\cos(\arcsin(ru + x) - qr)}{\sqrt{1 - (ru + x)^2}}du dx.
\end{split}
\end{equation}
Moreover, for every $x \in [-1, 1]$ such that $Y'(0)$ exists $\mathbb{P}_x$-almost surely, and for $u \in [q\sqrt{1 - x^2}, \pi\sqrt{1 - x^2}]$ for which $u$ is not an atom of the distribution $\mathbb{P}_x(Y'(0) \in \cdot)$, we have that
$$
\mathbb{P}_x\left(\frac{Y(r) - Y(0)}r > u \right) \xrightarrow[r \to 0]{} \mathbb{P}_x(Y'(0) > u).
$$
Therefore by the dominated convergence theorem, \eqref{E:b-c} converges to 
$$
\int_b^c \int_{q\sqrt{1 - x^2}}^{\pi\sqrt{1 - x^2}} \mathbb{P}_x(Y'(0) > u) du dx = \int_{b}^c \mathbb{E}_{x} (Y'(0) - \sqrt{1 - x^2}q)^+ dx \qquad \;\text{as}\; \quad r \to 0.
$$
The equality above again uses that $|Y'(0)| \le \pi\sqrt{1- Y^2(0)}$ almost surely, since $Y$ is $\pi\sqrt{1- h^2}$-Lipschitz.
Combining this with \eqref{E:want-1}, \eqref{E:J-+-da}, and \eqref{E:want-2} proves \eqref{E:deriv-integrated}.
Now define 
\[
d(Y) = \frac{Y'(0)}{\sqrt{1 - Y^2(0)}}.
\]
By \eqref{E:deriv-integrated}, for almost every $a \in (-1, 1)$, for every $q \in \mathbb{Q} \cap (\pi, \pi)$, we have that 
$$
\sqrt{1 - a^2}\mathbb{E}_{a} \left(d(Y) -  q\right)^+ = \mathbb{E}_{a} \left(Y'(0) - \sqrt{1 - a^2} q\right)^+ = \sqrt{1- a^2}D_\mu^+(q).
$$
Therefore by continuity in $s$ of the functions $\mathbb{E}_{a} \left(d (Y) - s\right)^+$ and $D_\mu^+(s)$, we have that
\begin{equation}
\label{E:daY}
\mathbb{E}_{a} \left(d (Y) - s\right)^+ = D_\mu^+(s) \qquad \text{ for all } s \in (-\pi, \pi).
\end{equation}
Now, $D_\mu(s)$ is Lipschitz and hence differentiable almost everywhere (Lemma \ref{L:convex}), so we can differentiate both sides of the above equation to get that for almost every $s \in (-\pi, \pi)$, that
$$
\mathbb{P}_a(d(Y) > s) = \mu(s, \infty).
$$
Finally, for almost every $a \in (-1, 1)$, we have that 
$$
\mathbb{P}_a(d(Y) \in [-\pi, \pi]) = \mu[\pi, \pi] = 1
$$
since $Y$ is almost surely $\pi\sqrt{1 - h^2}$-Lipschitz (Theorem \ref{T:bounded-speed}).
Therefore $\mathbb{P}_a(d(Y) \in \cdot) = \mu$.
\end{proof}

\subsection{The max-height distribution and the uniqueness of $Y$}

As an application of Proposition \ref{P:local-global-deriv}, we can find $\mathbb{P}(m(Y) > a)$ for all $a \in [0, 1]$. This allows us to deduce the uniqueness of $Y$.

\begin{theorem}
\label{T:max-height}
Let $Y$ be any subsequential limit of $Y_n$. For all $a \in [0, 1]$, we have that 
$ \mathbb{P}(m(Y) > a) = \sqrt{1 - a^2}$. That is, $m(Y) \stackrel{d}{=} \sqrt{1- U^2}$ for a uniform random variable $U$ on $[0, 1]$.
\end{theorem}

\begin{proof}
Let $V(Y)$ be the total variation of $Y$ on the interval $[0, 1]$. We have that
\begin{align*}
\mathbb{E} V(Y) &= \mathbb{E} \int_0^1 |Y'(t)|dt = \mathbb{E} \left(\frac{|Y'(0)|}{\sqrt{1-Y^2(0)}}\right) \sqrt{1 - Y^2(0)} = 2 \mathbb{E} \sqrt{1 - U^2},
\end{align*}
where $U$ is a uniform random variable on $[0, 1]$. Here the second equality follows from time stationarity of $Y$ (Proposition \ref{P:Y-symmetries} (i)). The third equality follows by Proposition \ref{P:local-global-deriv}, the fact that the first moment of $\mu$ is $2$ (Lemma \ref{L:expect-2}), and the fact that $Y(0)$ is uniformly distributed (Theorem \ref{T:bounded-speed}). 

\medskip

Moreover, the characterization of minimal flux paths (Theorem \ref{T:path-set}) implies that 
$
V(Y) = 2 \mathbb{E} m(Y),
$
so $\mathbb{E} m(Y) = \mathbb{E} \sqrt{1 - U^2}$.
Finally, the bound in Lemma \ref{L:max-heights} implies that 
$$
\mathbb{P}(m(Y) > a) \le \frac{D_\mu(0)}2 \sqrt{1 - a^2} = \sqrt{1- a^2} = \mathbb{P}(\sqrt{1 - U^2} > a),
$$
so the random variable $m(Y)$ is stochastically dominated by $\sqrt{1 - U^2}$. Hence $m(Y) \stackrel{d}{=} \sqrt{1 - U^2}$.
\end{proof}

Now we can prove the existence of a unique limit $Y$ of $Y_n$. First extend the paths $h_m$ defined in Theorem \ref{T:path-set} to paths $h_m:[0, \infty) \to [-1, 1]$ by letting $h_m(t) = -h_m(t + 1)$ for all $t > 1$.

\begin{theorem}
\label{T:Y-unique}
The sequence $Y_n$ has a distributional limit $Y$ given by
$$
Y(t) = h_{\sqrt{1 - U^2}}(V + t).
$$
Here $U$ is a uniform random variable on $[0, 1]$ and $V$ is uniform on $[0, 2]$, independent of $U$.
\end{theorem}

\begin{proof}
By Theorem \ref{T:unique-2} and Theorem \ref{T:particle-crossings} (iii), any subsequential limit $Y$ of $Y_n$ is supported on the set of paths 
$$
\{h_m (\cdot + u) : m \in [0, 1], u \in [0, 2]\}, 
$$ 
so we can write 
$
Y(\cdot) = h_M(\cdot + V),
$
for a pair of random variables $(M, V)$. By Theorem \ref{T:max-height}, $M = \sqrt{1 - U^2}$, for a uniform random variable $U$ on $[0, 1]$. By the time stationarity of $Y$ (Proposition \ref{P:Y-symmetries} (i)), for any $t \in [0, 2]$, we have that
$$
h_M(\cdot + t +  V) \stackrel{d}{=} h_M(\cdot +  V),
$$
which in turn implies that $(M, V) \stackrel{d}{=} (M, V + t \Mod 2 )$ for any $t \in [0, 2]$. This implies that $V$ has uniform distribution, and that $M$ is independent of $V$.
\end{proof}

\section{The integral transform formula}
\label{S:integral-formula}
Let $\mu_+$ be the pushforward of the local speed distribution $\mu$ under the map $x \mapsto |x|$.
In this section we prove an integral transform formula for $\mu_+$. This transform formula allows us to identify $\mu$ as the arcsine distribution on $[-\pi, \pi]$. Once we know $\mu$, we can find the set of minimal flux paths $h_m$ introduced in Theorem \ref{T:path-set}, and then in turn identify $Y$.

\medskip

 To find the integral transform formula, we will calculate $\mathbb{P}(m(Y) > k)$ in two different ways. We first give a heuristic explanation of how to do this when the local speed distribution $\mu$ has no atoms. By Theorem \ref{T:max-height}, we have that
$$
\mathbb{P}(m(Y) > k) = \sqrt{1 - k^2}.
$$
We can also calculate $\mathbb{P}(m(Y) > k)$ by integrating the marginal probabilities $\mathbb{P}_a(m(Y) > k)$. This gives that
\begin{equation}
\label{E:heur1}
\sqrt{1 - k^2} = \frac{1}2 \int_{-1}^1 \mathbb{P}_a(m(Y) > k)da = 1 - k + \frac{1}2 \int_{-k}^k \mathbb{P}_a(m(Y) > k)da.
\end{equation}
Now we want to find an expression for $\mathbb{P}_a(m(Y) > k)$ when $|a| < k$.  By the ordering on minimal flux paths, if $Y(0) = a$, then $m(Y) > k$ if and only if $|Y'(0)|$ is greater than some threshold value. Therefore for some $s_{a, k} \in \mathbb{R}$, we have that
\begin{equation}
\label{E:heur2}
\mathbb{P}_a(m(Y) > k) = \mathbb{P}_a (|Y'(0)| > \sqrt{1 - a^2}s_{a, k}) = \mu_+(s_{a, k}, \infty).
\end{equation}
The final equality here follows from Proposition \ref{P:local-global-deriv}. To find $s_{a, k}$, we calculate $\mathbb{P}(m(Y) > k | \; m(Y) > a)$. If we let 
\[
T_a = \inf \{t \in [0, 1] : Y(t) = a\},
\]
then we should have that
$$
\mathbb{P}(m(Y) > k | \; m(Y) > a) = \mathbb{P}(|Y'(T_a)| > \sqrt{1 - a^2}s_{a, k}).
$$
Now, the ``amount of time" that $Y$ spends at height $a$ is inversely proportional to its speed at that location. Therefore the distribution of $|Y'(T_a)|$ should be a size-biased version of the distribution of $Y'(0)$ given that $Y(0) = a$. Hence,  
$$
\mathbb{P}(|Y'(T_a)| > \sqrt{1 - a^2}s_{a, k}) = \hat{\mu}_+(s_{a, k}, \infty),
$$
where $\hat{\mu}_+$ is the size-biased distribution of $\mu_+$ (we define this formally in the next paragraph). Finally, we can also calculate $\mathbb{P}(m(Y) > k | \; m(Y) > a)$ using Theorem \ref{T:max-height}. This gives that
$$
\hat{\mu}_+(s_{a, k}, \infty) = \mathbb{P}(m(Y) > k | \; m(Y) > a) = \frac{\mathbb{P}(m(Y) > k)}{\mathbb{P}(m(Y) > a)} = \frac{\sqrt{1 - k^2}}{\sqrt{1 - a^2}}.
$$
We can combine this with \eqref{E:heur2} and \eqref{E:heur1} to get an integral transform formula involving the function $r_{\mu_+}(x) = S(\hat{S}^{-1}(x))$, where $S$ and $\hat{S}$ are the survival functions of $\mu_+$ and $\hat{\mu}_+$, respectively. 

\medskip

We now precisely define everything that is needed to state the integral transform formula. For a probability measure $\nu$ on $(0, \infty)$ with finite mean, define the {\bf size-biased distribution} $\hat{\nu}$ on $(0, \infty)$ by the Radon-Nikodym derivative formula
 $$
 \frac{d\hat{\nu}}{d \nu}(x) = \frac{x}{\int_0^\infty x d \nu(x)}.
 $$
In order to define the integral transform when $\mu$ has atoms, we define the \textbf{extended survival function} $S:\mathbb{R} \times [0,1] \to [0, 1]$ of a probability measure $\nu$ by
$$
S(x, q) = \nu(x, \infty) + (1 - q) \nu(x).
$$
$S$ is a non-increasing, continuous function in the lexicographic ordering on $\mathbb{R} \times [0,1]$. $S$ can be thought of as the survival function in the lexicographic ordering for $\nu \times \mathcal{L}$, where $\mathcal{L}$ is uniform measure on $[0, 1]$. Now define the \textbf{ size-bias ratio function} $r_\nu(x):[0, 1] \to [0, 1]$ of a probability measure $\nu$ on $(0, \infty)$ with finite mean by
\begin{equation*}
r_\nu(x) = S(\hat{S}^{-1}(x)),
\end{equation*}
where $\hat{S}$ is the extended survival function of $\hat{\nu}$ and $S$ is the extended survival function of $\nu$. Here the inverse function is given by
$$
\hat{S}^{-1}(x) = \sup \{(y, q) \in \mathbb{R} \times [0,1] : \hat{S}(y, q) = x \},
$$
where the supremum is taken with respect to the lexicographic ordering. When $\nu$ has no atoms, we can define $r_\nu$ in terms of the usual survival functions. 
\begin{prop}[The integral transform] 
\label{P:integral-equality} Let $\mu_+$ be the pushforward of the measure $\mu$ under the map $f(x) = |x|$. 
For every $k \in (0, 1)$, we have that
\begin{equation}
\label{E:integral-H}
 1 - k + \int_0^k r_{\mu_+} \left(\frac{\sqrt{1 - k^2}}{\sqrt{1- x^2}}\right) dx = \sqrt{1-k^2}.
\end{equation}
\end{prop}

To prove Proposition \ref{P:integral-equality}, we first establish that $\mu(0) = 0$. 
\begin{lemma}
\label{L:X-atom}
$\mu(0) = 0$.
\end{lemma}

\begin{proof} 
By Lemma \ref{L:Y-sym-type}, $\mathbb{P} \left(\frac{Y'(0)}{\sqrt{1 - Y^2(0)}} = 0 \right) = 0$. By Proposition \ref{P:local-global-deriv}, this is equal to $\mu(0)$. \end{proof}

Now let $m < 1$, and let $h_m$ be as in Theorem \ref{T:path-set}. For any $a \in [0, m]$, define
$$
s_{a, m} = \frac{1}{\sqrt{1 - a^2}} \liminf_{r \to 0} \frac{|h_m(t^* + r) - h_m(t^*)|}{r},
$$
where $t^*$ is any point where $h_m(t^*) = a$. Note that $s_{a, m}$ is independent of $t^*$ by the symmetry of minimal flux paths (Theorems \ref{T:path-set} (i)).
We now prove an integral formula relating the speeds $s_{a, m}$ to the local speed distribution.

\begin{lemma}
\label{L:bias-unbias}
Let $\nu$ be the law of $m(Y)$. For almost every $a \in [0, 1)$, for every $k \in (a, 1)$ there exists a constant $q_{a, k} \in [0, 1]$ such that
\begin{equation}
\begin{split}
\mathbb{P}_a(m(Y) > k) &= \frac{1}{\sqrt{1 - a^2}} \int \mathbbm{1} \Big(s_{a, m} > s_{a, k} \;\; \;\text{or}\; \;\; s_{a, m} = s_{a, k}, \; m > k\Big) \frac{2}{s_{a, m}} d\nu(m) \\
\label{E:s-a-k}
&= \mu_+(s_{a, k}, \infty) + q_{a, k}\mu_+(s_{a, k}).
\end{split}
\end{equation}
Moreover, $\mathbb{P}_a(m(Y) > k)$ is a continuous function of $k \in (a, 1)$ for almost every $a \in [0, 1)$.
\end{lemma}

\begin{proof}
Let $h_m$ be as in Theorem \ref{T:path-set}, and let $a < b \in [0, 1)$. We first compute the amount of time that $h_m$ spends in the interval $[a, b]$. Define $h_m^{-1}:[0, m] \to [0, 1/2]$ by
$$
h_m^{-1}(x) = \inf \{t : h_m(t) = x \}.
$$
By the strict monotonicity and symmetry of $h_m$ (Theorem \ref{T:path-set} (i), (ii)), we can write
\begin{equation}
\label{E:L-h}
\mathcal{L}\{t: |h_m(t)| \in [a, b]\} = 2[h_m^{-1}(b) -  h_m^{-1}(a)],
\end{equation}
where $\mathcal{L}$ is Lebesgue measure on $[0, 1]$.
Now by the concavity of $h_m$ (Lemma \ref{L:concave-paths}), the inverse $h_m^{-1}(x)$ is almost everywhere differentiable  with derivative $2/[s_{x, m} \sqrt{1 - x^2}]$. Therefore the left hand side of \eqref{E:L-h} is equal to 
$$
\int_a^{b} \frac{2}{s_{x, m} \sqrt{1 - x^2}}dx.
$$
Now letting $U$ be a uniform random variable on $[0, 1]$ that is independent of $Y$, for any $a < b < k \in [0, 1)$, we have that
\begin{equation*}
\begin{split}
\mathbb{P}\big(m(Y) > k \;\text{and}\; Y(0) \in [a, b]\big) &= \frac{1}2\mathbb{P}\big(m(Y) > k \;\text{and}\; |Y(U)| \in [a, b]\big) \\
&= \frac{1}2 \int_k^1 \mathcal{L}\{t: |h_m(t)| \in [a, b]\} d\nu(m) \\
&= \int_a^{b} \frac{1}{\sqrt{1 - x^2}} \int_k^1 \frac{1}{s_{x, m}} d\nu(m)dx.
\end{split}
\end{equation*}
The first equality above follows by the time-stationarity and symmetry of $Y$ (Proposition \ref{P:Y-symmetries} (i) and (ii)). The second equality follows since $Y$ is supported on shifts of the minimal flux paths $h_m$ (Theorem \ref{T:unique-2}). This implies that for almost every pair $a < k \in [0, 1)$, that
\begin{equation}
\label{E:f-amid}
\mathbb{P}_a(m(Y) > k) = \frac{1}{\sqrt{1 - a^2}} \int_k^1 \frac{2}{s_{a, m}} d\nu(m).
\end{equation}
Now by the ordering on minimal flux paths (Theorem \ref{T:unique-2}(iv)), we have that
\begin{equation}
\label{E:speed-mono}
\text{if} \;\; s_{a, m(Y)} < s_{a, k}, \qquad \text{then } m(Y) < k.
\end{equation}
This allows us to rewrite \eqref{E:f-amid} to get that
\begin{equation}
\label{E:f-amid-2}
\mathbb{P}_a(m(Y) > k) = \frac{1}{\sqrt{1 - a^2}} \int \mathbbm{1} \Big(s_{a, m} > s_{a, k} \;\; \;\text{or}\; \;\; s_{a, m} = s_{a, k}, m > k\Big) \frac{2}{s_{a, m}} d\nu(m)
\end{equation}
for almost every $a < k \in [0, 1)$.
Moreover, for almost every $a \in [0, 1)$, we have that
$$
|Y'(0)| = \sqrt{1 - a^2}s_{a, m(Y)} \qquad \mathbb{P}_a\text{-almost surely}.
$$
This follows by time stationarity of $Y$, since $Y$ is almost everywhere differentiable. Also, $Y'(0)/\sqrt{1- a^2}$ has distribution $\mu$ for almost every $a$ (Proposition \ref{P:local-global-deriv}). Therefore \eqref{E:speed-mono} implies that for almost every $a \in [0, 1)$, for every $k \in (a, 1)$ there exists a constant $q_{a, k} \in [0, 1]$ such that
\begin{equation}
\label{E:Pa-mY}
\mathbb{P}_a(m(Y) > k) = \mu_+(s_{a, k}, \infty) + q_{a, k}\mu_+(s_{a, k}).
\end{equation}
Now let $a$ be such that \eqref{E:Pa-mY} holds for every $k \in (a, 1)$ and \eqref{E:f-amid} and \eqref{E:f-amid-2} hold for almost every $k \in (a, 1)$. 

\medskip

By concavity of minimal flux paths (Lemma \ref{L:concave-paths}), $s_{a, m} > 0$ whenever $m > a$. Therefore since $\nu$ has a Lebesgue density by Theorem \ref{T:max-height}, the right hand side of \eqref{E:f-amid} is continuous and non-increasing. Since both sides of \eqref{E:f-amid-2} are also non-increasing, and are equal to the right hand side of \eqref{E:f-amid} for almost every $k \in (a, 1)$, they must be equal for every $k \in (a, 1)$. Combining this with \eqref{E:Pa-mY} proves the lemma.
\end{proof}


\begin{lemma}
\label{L:equal-measure} Let $\mathfrak{s}_a$ be the law of $\mathbb{P}(s_{a, m(Y)} \in \cdot \; | m(Y) > a)$, and define the measure $\bar{\mathfrak{s}}_a$ by the Radon-Nikodym formula 
$$
\frac{2}{s} d\mathfrak{s}_a(s) = d\bar{\mathfrak{s}}_a(s).
$$
Then for almost every $a \in [0, 1)$, we have that $\bar{\mathfrak{s}}_a = \mu_+$.
In particular, for such $a$, we have that
\begin{equation}
\label{E:q-ak}
q_{a, k}\mu_+(s_{a, k}) = \frac{2}{s_{a, k}\sqrt{1 - a^2}} \int \mathbbm{1} \Big(s_{a, m} = s_{a, k}, m > k\Big)d\nu(m).
\end{equation}
\end{lemma}

\begin{proof}
Let $a$ be such that \eqref{E:s-a-k} holds for every $k \in (0, 1)$ and such that $\mathbb{P}_a(m(Y) > k)$ is continuous. Further assume that
$$
\mathbb{P}_a\left( Y'(0) \text{ exists and is non-zero }, m(Y) < 1, J(Y) = 1 \right) = 1.
$$
These conditions hold for almost every $a \in [0, 1)$ (Proposition \ref{P:local-global-deriv}, Lemma \ref{L:bias-unbias}). Noting that $\sqrt{1 - a^2} = \mathbb{P}(m(Y) > a)$ by Theorem \ref{T:max-height}, equation \eqref{E:s-a-k} implies that for every $k \in (a, 1)$, there exists a $p_{a, k} \in [0, 1]$ such that
\begin{equation}
\label{E:comp}
\mathbb{P}_a(m(Y) > k) = \bar{\mathfrak{s}}_a(s_{a, k}, \infty) + p_{a, k}\bar{\mathfrak{s}}_a(s_{a, k}) = \mu_+(s_{a, k}, \infty) + q_{a, k}\mu_+(s_{a, k}).
\end{equation}
Now, since $Y'(0) \ne 0$, $\mathbb{P}_a$-almost surely, the concavity of minimal flux paths (Lemma \ref{L:concave-paths}) implies that $\mathbb{P}_a(m(Y) = a) = 0$. Moreover, since $m(Y) < 1$, $\mathbb{P}_a$-almost surely, we have that $\mathbb{P}_a(m(Y) \in (a, 1)) = 1$, and so
\begin{equation}
\label{E:mY}
\lim_{k \to 1} \mathbb{P}_a(m(Y) > k) = 0 \quad \;\text{and}\; \quad \lim_{k \to a} \mathbb{P}_a(m(Y) > k) = 1.
\end{equation}
Therefore, \eqref{E:comp} and the continuity of $\mathbb{P}_a(m(Y) > k)$ implies that 
\begin{equation}
\label{E:full}
\mu_+ \left( s_{a, k} : k \in (0, 1) \right) = \bar{\mathfrak{s}}_a \left( s_{a, k} : k \in (0, 1) \right) = 1.
\end{equation}
Now fix $k \in (a, 1)$. Let 
$$
k^* = \inf \{\ell  \in [k, 1] :  s_{a, k} < s_{a, m} \;\; \text{ for all } \;\;m > \ell \}.
$$
Equation \eqref{E:comp} and the continuity of $\mathbb{P}_a(m(Y) > k)$ then implies that
$$
\mathbb{P}_a(m(Y) > k^*) = \bar{\mathfrak{s}}_a(s_{a, k}, \infty) = \mu_+(s_{a, k}, \infty).
$$
Combining this with \eqref{E:full} proves that $\mu_+ = \bar{\mathfrak{s}}_a$. Equation \eqref{E:q-ak} follows by using that $\mu_+ = \bar{\mathfrak{s}}_a$ to simplify equation \eqref{E:s-a-k}.
\end{proof}

\begin{proof}[Proof of Proposition \ref{P:integral-equality}]
Fix $a$ so that the conclusion of Lemma \ref{L:equal-measure} holds, and let $k > a$. By Theorem \ref{T:max-height} and the ordering on minimal flux paths, we can write 
\begin{align*}
\nonumber\mathbb{P}(m(Y) > k | \; m(Y) > a)  &= \frac{1}{\sqrt{1 - a^2}} \left[\int \mathbbm{1}(s_{a, m} > s_{a, k} ) d\nu(m)  + \int \mathbbm{1}(s_{a, m} = s_{a, k}, m > k) d\nu(m) \right].
\end{align*}
By Lemma \ref{L:equal-measure}, we can rewrite the first integral above in terms of $\mathfrak{s}_a$, and then in terms of $d\mu_+$. This gives that
$$
\frac{1}{\sqrt{1 - a^2}}\int \mathbbm{1}(s_{a, m} > s_{a, k} ) d\nu(m) = \int \mathbbm{1}(s > s_{a, k}) d\mathfrak{s}_a(s) = \int \mathbbm{1}(s > s_{a, k})\frac{s}2d\mu_+(s).
$$
We can also rewrite the second integral using \eqref{E:q-ak}. This implies that
$$
\mathbb{P}(m(Y) > k | \; m(Y) > a) = \int \mathbbm{1}(s > s_{a, k})\frac{s}2d\mu_+(s) + \frac{s_{a, k}}2 q_{a, k} \mu_+(s_{a, k}).
$$
By Lemma \ref{L:expect-2}, we can recognize $\frac{s}{2}$ as the Radon-Nikodym derivative of the size-biased random variable $\hat{\mu}_+$ with respect to $\mu_+$, proving that
$$
\frac{\sqrt{1 - k^2}}{\sqrt{1 - a^2}} = \mathbb{P}(m(Y) > k | m(Y) > a) = \hat{\mu}(s_{a, k}, \infty) + q_{a, k}\hat{\mu}(s_{a, k})
$$
for almost every $a < k \in [0, 1)$. Now, Lemma \ref{L:bias-unbias} allow us to conclude that 
\begin{align*}
\mathbb{P}_a(m(Y) > k) &= S(s_{a, k}, 1 - q_{a, k}) = S\left(\hat{S}^{-1}\left(\frac{\sqrt{1- k^2}}{\sqrt{1- a^2}}\right)\right)
\end{align*}
for almost every $a < k \in [0, 1)$. Combining this with the symmetry of $Y$ (Proposition \ref{P:Y-symmetries}) and Theorem \ref{T:max-height} implies that
\[
\sqrt{1 - k^2} = \mathbb{P}(m(Y) > k) = \int_{0}^1 \mathbb{P}_a(m(Y) > k) da = 1 - k + \int_0^k r_{\mu_+}\left(\frac{\sqrt{1- k^2}}{\sqrt{1- a^2}}\right) da. \qedhere
\]
\end{proof}
\section{The weak trajectory limit}
\label{S:transform}

In this section we show that the integral transform in Proposition \ref{P:integral-equality} determines the local speed distribution. We begin with two basic lemmas about size-bias ratio functions.

\begin{lemma} 
\label{L:loc-lip}
For any probability distribution $\nu$ on $(0, \infty)$ with finite first moment $m$, the size-bias ratio function $r_\nu$ is locally Lipschitz on $[0, 1)$ and continuous on $[0, 1]$.
\end{lemma}

\begin{proof}
Let $\pi_1(x, y) = x.$ By calculus, when $x \in (0, 1)$ we have that
$$
\partial_+ r_\nu(x) = \lim_{h \to 0+} \frac{m}{\pi_1(\hat{S}^{-1}(x + h))} \qquad \;\text{and}\; \qquad
\partial_- r_\nu(x) = \lim_{h \to 0+} \frac{m}{\pi_1(\hat{S}^{-1}(x - h))}.
$$
The first equation also holds as $x = 0$. As $\pi_1(\hat{S}^{-1}(y))$ is a decreasing function of $y$, and strictly positive for all $y \in [0, 1)$, this shows that $r_\nu$ is locally Lipschitz on $[0, 1)$. It is straightforward to check that $r_\nu$ is continuous at $1$.
\end{proof}

\begin{lemma}
\label{L:rX-det}
Suppose that $\nu_1$ and $\nu_2$ are probability measures on $(0, \infty)$ with the same first moment, such that $r_{\nu_1} = r_{\nu_2}$. Then $\nu_1 = \nu_2$.
\end{lemma}

\begin{proof} We prove the contrapositive.
Let $\nu_1$ and $\nu_2$ be measures with the same first moment, and suppose that $\nu_1 \ne \nu_2$. Since $\nu_1$ and $\nu_2$ have the same first moment, $\hat{\nu}_1 \ne \hat{\nu}_2$, so  for some value of $(y, q) \in (0, \infty) \times [0, 1]$, we have that $\hat{S}_{\nu_1} (y, q) \ne \hat{S}_{\nu_2} (y, q)$. Since the functions $\hat{S}_{\nu_1}(y, \cdot)$ and $\hat{S}_{\nu_2}(y, \cdot)$ are linear for any fixed value of $y$, this implies that $\hat{S}_{\nu_1} (y, r) \ne \hat{S}_{\nu_2} (y, r)$ for some $r \in \{0, 1\}$. Moreover, since 
$$
\hat{S}_{\nu_i}(y, 1) = \sup \{\hat{S}_{\nu_i}(z, 0) : z > y \}, \qquad i = 1, 2,
$$
we can conclude that $\hat{S}_{\nu_1} (z, 0) \ne \hat{S}_{\nu_2} (z, 0)$ for some $z \in (0, \infty)$. Without loss of generality, assume that $a_1 := \hat{S}_{\nu_1} (z, 0) > a_2 :=\hat{S}_{\nu_2} (z, 0)$. Therefore using the notation of the previous lemma, letting $b = (a_1 + a_2)/2$, we have that
$$
\lim_{h \to 0^+} \pi_1(\hat{S}_{\nu_1}^{-1}(b + h)) \ge z, \qquad \text{whereas} \qquad \lim_{h \to 0^+} \pi_1(\hat{S}_{\nu_1}^{-1}(b + h)) < z.
$$
Since $b \in (0, 1)$, the derivative computation in the previous lemma combined with the fact that $\nu_1$ and $\nu_2$ have the same first moment implies that $r_{\nu_1} \ne r_{\nu_2}$.
\end{proof}
Now let $\mathcal{X}$ be the space of continuous functions from $[0, 1] \to \mathbb{R}$ that are locally Lipschitz on $[0, 1)$. Define an integral transform $H$ on $\mathcal{X}$ by 
\begin{equation*}
H(r)(k) = \int_0^k r \left(\frac{\sqrt{1 - k^2}}{\sqrt{1- x^2}}\right) dx,
\end{equation*}
for $k \in [0, 1]$. By Lemma \ref{L:loc-lip}, any linear combination of size-bias ratio functions is in $\mathcal{X}$, so if the integral transform $H$ is injective on $\mathcal{X}$, then $H(r_{\mu_+})$ determines $r_{\mu_+}$.

\begin{lemma}
\label{L:injective-transform-lip}
$H$ is injective on $\mathcal{X}$.
\end{lemma}

\begin{proof}
Let $r \in \mathcal{X}, r \ne 0$. We will show that $H(r) \ne 0$. Without loss of generality, we may assume that 
$$
\max_{x \in [0, 1]} r(x) = \delta > 0.
$$
Letting $u = \frac{\sqrt{1 - k^2}}{\sqrt{1- x^2}}$ and $y = \sqrt{1 -k^2}$, we have that
\begin{equation*}
H(r)(\sqrt{1-y^2}) = \int_y^1 r(u) \frac{y^2}{u^2\sqrt{u^2 - y^2}}du.
\end{equation*}
For $y \in (0, 1]$, we have that
\begin{equation}
\label{E:H-bar}
\close{H}(r)(y) := \frac{H(r)(\sqrt{1 - y^2})}{y^2} = \int_y^1 \frac{r(u)}{u^2\sqrt{u^2 - y^2}}du.
\end{equation}
It suffices to show that $\close{H}(r)(y) \ne 0$ for some $y \in (0, 1]$. 
Observe that if $r(1) > 0$, then there would exist a $\gamma > 0$ such that $r(x) > 0$ for all $x > \gamma$, so $\close{H}(r)(\gamma) > 0$. Also, if $r(0) > 0$, then $\close{H}(r)(y) \to \infty$ as $y \to 0$. Therefore there must exist $y \in (0,1)$ be such that  $r(y) = \delta$.

 Since $r$ is locally Lipschitz on $[0, 1)$, we can find $k, \epsilon > 0$ such that $r(x) \ge \delta - k \epsilon$ for all $x \in [y, y + \epsilon]$. 
Therefore
\begin{align}
\label{E:de-k-bd}
\int_y^{y+\epsilon}\frac{r(u)}{u^2\sqrt{u^2 - y^2}}du \ge \frac{(\delta - k \epsilon)\sqrt{\epsilon^2 + 2y\epsilon}}{y^2(y + \epsilon)}.
\end{align}

Now consider the difference between $\close{H}(r)(y)$ and $\close{H}(r)(y + \epsilon)$. We have 
\begin{align}
\nonumber \close{H}&(r)(y) - \close{H}(r)(y + \epsilon) \\
\nonumber &= \int_y^{y+\epsilon}\frac{r(u)}{u^2\sqrt{u^2 - y^2}}du + \int_{y+ \epsilon}^1 \frac{r(u)}{u^2\sqrt{u^2 - y^2}}du - \int_{y+ \epsilon}^1 \frac{r(u)}{u^2\sqrt{u^2 - (y+ \epsilon)^2}}du  \\
\nonumber &\ge \int_y^{y+\epsilon}\frac{r(u)}{u^2\sqrt{u^2 - y^2}}du + \delta \int_{y + \epsilon}^1 \left[\frac{1}{u^2\sqrt{u^2 - y^2}} - \frac{1}{u^2\sqrt{u^2 - (y + \epsilon)^2}}\right]du \\
\label{E:int-alm} &\ge \frac{(\delta - k \epsilon) \sqrt{\epsilon^2 + 2y\epsilon}}{y^2(y + \epsilon)} + \delta \left[\frac{(y + \epsilon)^2\sqrt{1-y^2} -y^2\sqrt{1-(y + \epsilon)^2}- (y + \epsilon)\sqrt{\epsilon^2 + 2y\epsilon}}{(y + \epsilon)^2y^2}\right] .
\end{align}
In above calculation the first inequality comes from the fact that $r(x) \le \delta$ for all $x$, and the observation that
$$
\frac{1}{u^2\sqrt{u^2 - y^2}} < \frac{1}{u^2\sqrt{u^2 - (y + \epsilon)^2}}
$$ 
for all $u \in (y + \epsilon, 1)$. The second equality follows by integration and plugging in the bound in \eqref{E:de-k-bd}. Now expanding in $\epsilon$ about $\epsilon = 0$, we get that \eqref{E:int-alm} is equal to \begin{align*}
\frac{(2-y^2)\delta\epsilon}{y^3\sqrt{1-y^2}} + O(\epsilon^{3/2}).
\end{align*}
This is strictly greater than 0 for small enough $\epsilon$, so
$
\close{H}(r)(y) - \close{H}(r)(y + \epsilon) \ne 0
$
for such $\epsilon$. Hence $\close{H}(r)(x) \ne 0$ for some $x \in [0, 1)$. 
\end{proof}

\begin{prop}
\label{P:arcsin-en} Let $\mathfrak{arc}$ be the arcsine distribution on $[-\pi, \pi]$, and let $\mathfrak{arc}_+$ be the pushforward of $\mathfrak{arc}$ under the map $x \mapsto |x|$. Then for every $k \in [0, 1]$, we have that
\begin{equation}
\label{E:integral-H-2}
 1 - k + \int_0^k r_{\mathfrak{arc}_+} \left(\frac{\sqrt{1 - k^2}}{\sqrt{1- x^2}}\right) dx = \sqrt{1-k^2}.
\end{equation}
\end{prop}

\begin{proof} The distribution $\mathfrak{arc}$ has density $(\pi \sqrt{\pi^2 - x^2})^{-1}$ on $[-\pi, \pi]$. From this, we can calculate that
\begin{align*}
r_{\mathfrak{arc}_+}(y) = 1 - \frac{2}{\pi}\arcsin(\sqrt{1-y^2}).
\end{align*}
We now use the connection between the arcsine distribution and the Archimedean measure $\mathfrak{Arch}$ to help evaluate the integral in \eqref{E:integral-H-2}. Let $(X_1, X_2) \sim \mathfrak{Arch}$ and let 
$$
\mathcal{A}(t) = X_1 \cos (\pi t) + X_2 \sin(\pi t)
$$
be the Archimedean path. Writing $\mathbb{P}(m(\mathcal{A}) > k) = \mathbb{P}(X_1^2 + X_2^2 > k^2)$ in both polar and Cartesian coordinates, we get that
$$
\int_k^1 \frac{rdr}{\sqrt{1 - r^2}} = 4 \int_{0}^1 \int_{\sqrt{k^2 - y^2} \vee 0}^{\sqrt{1 - y^2}} \frac{dxdy}{2 \pi \sqrt{1 - x^2 - y^2}}.
$$
The left hand side can easily be evaluated as $\sqrt{1 - k^2}$. The right hand side is equal to \begin{align*}
(1 - k) + \int_{0}^k \frac{2}{\pi} \left(\frac{\pi}{2} - \arcsin\left(\frac{\sqrt{k^2 - y^2}}{\sqrt{1- y^2}}\right)\right)dy
= (1- k) + \int_{0}^k r_\mathfrak{arc}\left(\frac{\sqrt{1 - k^2}}{\sqrt{1- y^2}}\right) dy. \qquad \qedhere
\end{align*}

\end{proof}

We can now prove Theorem \ref{T:main-2}, and in turn use that to prove Theorem \ref{T:weak-limit}.

\begin{proof}[Proof of Theorem \ref{T:main-2}]
By Proposition \ref{P:integral-equality}, Lemma \ref{L:injective-transform-lip} and Proposition \ref{P:arcsin-en}, we have that  $r_{\mu_+} = r_{\mathfrak{arc}_+}$. Moreover, the first moment of $\mathfrak{arc}_+$ is $2$. By Lemma \ref{L:expect-2}, this matches the first moment of $\mu_+$. Therefore by Lemma \ref{L:rX-det}, $\mu_+ = \mathfrak{arc}_+$. Finally, symmetry of both measures implies that $\mu = \mathfrak{arc}$.
\end{proof}


\begin{proof}[Proof of Theorem \ref{T:weak-limit}] Fix $m > 0$.
Let $g_m(t) = m \sin (\pi t)$, and let $g_m^{-1}$ be the inverse of $g_m$ on the interval $[0, 1/2]$. Let 
$$
s_m(t) = \frac{g_m'(t)}{\sqrt{1 - g_m^2(t)}}
$$
be the local speed of $g_m$. We can calculate that
\begin{align}
\nonumber J(g_m) = 2J(g_m; [0, 1/2]) &= \int_0^{1/2} D_\mu(s_m(t)) \sqrt{1 - g_m^2(t)}dt \\
\label{E:s-h}
&= \int_0^m \frac{D_\mu(s_m(g_m^{-1}(x))}{s_m(g_m^{-1}(x))}dx.
\end{align}
Here we have made the substitution $x = g_m(t)$ to go from the first to the second line. Now using Theorem \ref{T:main-2}, we can calculate
$$
D_\mu(c) = D_\mathfrak{arc}(c) = \frac{2}{\pi} \left(c \arcsin \left(\frac{c}{\pi}\right) + \sqrt{\pi^2 - c^2}\right).
$$
From here we can use that $s_m(g_m^{-1}(x)) = \frac{\pi \sqrt{m^2 - x^2}}{\sqrt{1 - x^2}}$ to compute that
\begin{equation*}
J(g_m) = \int_0^m \frac{2}{\pi} \left[ \frac{\sqrt{1-m^2}}{\sqrt{m^2 - x^2}} + \arcsin\left(\frac{\sqrt{m^2 - x^2}}{\sqrt{1- x^2}}\right)\right]dx
\end{equation*}
The first part of this integral can be easily evaluated, and the second part can be evaluated by comparing with the integral in the final line of the proof of Proposition \ref{P:arcsin-en}. Putting this all together yields that $J(g_m) = 1$, as desired. Therefore by Theorem \ref{T:Y-unique}, we can write
$$
Y(t) = \sqrt{1 - V^2}\sin(\pi t + 2 \pi U),
$$
where $U$ and $V$ are independent uniform random variables on $[0, 1]$. This is the Archimedean path.
\end{proof}

\subsection{The empirical distribution of trajectories}
\label{SS:slightly}

By a compactness argument, we can immediately prove a slightly stronger version of Theorem \ref{T:weak-limit}. This will allow us to conclude Theorem \ref{T:subnetwork}. This theorem will also be necessary for establishing the stronger limits in Sections \ref{S:strong-limit} and \ref{S:geom-limit}.

\medskip

For a Polish space $S$, let $\mathcal{M}(S)$ be the space of probability measures on $S$ with the topology of weak convergence. Note that $\mathcal{M}(S)$ is itself a Polish space.
Recall the notation $\sigma_G$ introduced in Section \ref{S:intro}. For a fixed sorting network $\sigma$, let $\nu_\sigma \in \mathcal{M}(D)$ be uniform measure on the set $\{\sigma_G(i, \cdot)\}_{i \in \{1, \dots, n\}}$. Now let $\Omega_n$ be the space of all $n$-element sorting networks, and define 
$$
\nu_n = \frac{1}{\card{\Omega_n}} \sum_{\sigma \in \Omega_n} \delta(\nu_\sigma).
$$
For each $n$, $\nu_n \in \mathcal{M}(\mathcal{M}(\mathcal{D}))$. We now extend Theorem \ref{T:weak-limit} to give a limit theorem for the sequence $\nu_n$. 
\begin{theorem}
\label{T:weak-limit'}
Let $\nu_n \in \mathcal{M}(\mathcal{M}(\mathcal{D}))$ be as above, and let $\mathbb{P}_\mathcal{A}$ be the law of the Archimedean path $\mathcal{A}(t) = X_1\cos(\pi t) + X_2 \sin (\pi t)$, where $(X_1, X_2) \sim \mathfrak{Arch}$. Then $\nu_n \to \delta_{\mathbb{P}_\mathcal{A}}$.
\end{theorem}

Note that the law of $Y_n$ is the expectation of the random measure $\nu_n$. Therefore Theorem \ref{T:weak-limit} gives convergence in expectation for $\nu_n$, whereas Theorem \ref{T:weak-limit'} gives pointwise convergence.

\begin{proof}
By Remark 2.4 from \cite{dauvergne1}, the sequence $\nu_n$ is precompact. For any subsequential limit $\nu$ of $\nu_n$, Theorem \ref{T:weak-limit} implies that $\nu$-almost every $\rho$ is supported on curves of the form $a\sin(\pi t) + b \cos(\pi t)$. Hence if $Z$ is a random path with law $\rho$, then
$$
Z(t) = X_1 \cos(\pi t) + X_2 \sin(\pi t)
$$
for some random variables $(X_1, X_2)$. Moreover, $Z(t)$ is uniform for every $t$ since each measure $\nu_n$ is supported on the set $\{\nu_\sigma\}_{\sigma \in \Omega_n}$.  Therefore $(X_1, X_2) \sim \mathfrak{Arch}$, so $Z$ is the Archimedean path, and hence $\nu = \delta_{\mathbb{P}_\mathcal{A}}$.
\end{proof}

\begin{proof}[Proof of Theorem \ref{T:subnetwork}]
Fix $m$. With the notation of Theorem \ref{T:subnetwork}, we have that
\begin{equation}
\label{E:2-cvg}
\tau^n_m \stackrel{d}{\to} \tau_m \quad \;\text{as}\; n \to \infty \quad \text{ if and only if } \quad
\frac{1}{|\Omega_n|} \sum_{\sigma \in \Omega_n} \nu_n^m \to \delta_{\mathbb{P}_\mathcal{A}^m} \;\text{as}\; n \to \infty.
\end{equation}
in the space $\mathcal{M}(\mathcal{M}(\mathcal{D}^n))$. Here $\nu_n^m$ is the $m$-fold product measure $\nu_n \times \nu_n \times \dots \nu_n$. The reason for this is that is if $I_1, \dots, I_m$ are uniform random variables on $\{1, \dots n\}$, then $I_j \ne I_k$ for all $j, k \in \{1, \dots n\}$ with high probability. Hence to find the limit of an $m$-out-of-$n$ sorting network, it suffices to find the limit of $m$ independently chosen trajectories. Finally, the second convergence statement in \eqref{E:2-cvg} follows from Theorem \ref{T:weak-limit'}.
\end{proof}

\begin{section}{The strong sine curve limit}
\label{S:strong-limit}
Theorem \ref{T:weak-limit} shows that for any $\epsilon > 0$, with high probability $(1- \epsilon)n$ particle trajectories in a random sorting network are close to sine curves. In this section, we extend this result to all particle trajectories, thus proving Theorem \ref{T:sine-curves}. By combining Theorem \ref{T:sine-curves} and Theorem \ref{T:weak-limit}, we also prove Theorems \ref{T:matrices} and \ref{T:unif-rotation}.

\medskip

The idea behind the proof is as follows. By Theorem \ref{T:weak-limit}, we know that most trajectories in a typical large-$n$ sorting network are close to sine curves. Since any two trajectories must cross exactly once, this restricts the type of behaviour that the remaining trajectories can have. Specifically, this forces all other trajectories to be either sine curves themselves, or to spend a lot of time at the edge of the sorting network. We can eliminate this second case by using the ``octagon" bound from \cite{angel2007random}. To state this bound, let $A_{n, \gamma}$ be the event where
\begin{align*}
\left|\sigma^n_G(i, t) - \sigma^n_G(i, 0)\right| < 2\sqrt{2t - t^2} + \gamma \;\; &\;\text{and}\; \;\; \left|\sigma^n_G(i, t) - \sigma^n_G(i, 1)\right| < 2\sqrt{1 - t^2} + \gamma \\
&\text{for all } \;\; t \in [0, 1], i \in \{1, \mathellipsis n\}.
\end{align*}


\begin{theorem}[Octagon bound, \cite{angel2007random}]
\label{T:octagon}
For any $\gamma > 0$, we have that 
$$
\lim_{n \to \infty} \mathbb{P}(A_{n, \gamma}) = 1.
$$
\end{theorem}

Now define 
$$
L^n_{i}(\delta) = \mathcal{L} \{ t: |\sigma^n_G(i, t)| \ge 1 - \delta\},
$$
where $\mathcal{L}$ is Lebesgue measure on $[0, 1]$. This is the amount of time that particle $i$ spends within $\delta$ of the edge in the random sorting network $\sigma^n$.
\begin{lemma}
\label{L:edge-time}
For every $\epsilon > 0$, we have that
$$
\lim_{n \to \infty} \mathbb{P} \left( \max_{i \in [1, n] } L^n_i(\epsilon^2/16) > \epsilon \right) =0.
$$
\end{lemma}

\begin{proof}
Fix $\epsilon > 0, n \in \mathbb{N}$, and let $i$ be such that $|\sigma_G(i, 0)| \ge 1 - \epsilon^2/16$. On the event $A_{n, \epsilon^4/64}$, a simple calculation shows that
$$
-1 + \epsilon^2/16 < \sigma_G(i, t) < 1 - \epsilon^2/16 \qquad \text{ for all } t \in [\epsilon/2, 1 - \epsilon/2].
$$
Therefore by Theorem \ref{T:octagon}, we have that
\begin{equation}
\label{E:time-shift}
\lim_{n \to \infty} \mathbb{P} \bigg( \max \left\{L^n_i(\epsilon^2/16)  : \; i \in [1, n], \;|\sigma^n_G(i, 0)| \ge 1 - \epsilon^2/16 \right\} > \epsilon \bigg) = 0.
\end{equation}
By time stationarity of random sorting networks (Theorem \ref{T:time-stat}), for each $n$ the above probability is greater than or equal to 
$$
\epsilon \mathbb{P} \left( \max_{i \in [1, n] } L^n_i(\epsilon^2/16) > \epsilon \right).
$$
Therefore \eqref{E:time-shift} implies the lemma.
\end{proof}

Recall that $\nu_\sigma$ is uniform measure on the trajectories of a sorting network $\sigma$, and that $\mathbb{P}_\mathcal{A}$ is the law of the Archimedean path.
To prove Theorem \ref{T:sine-curves} from here, we must show that if $\nu_\sigma$ is close to $\mathbb{P}_\mathcal{A}$ in the weak topology, then for any particle $i$, either $\sigma_G(i, \cdot)$  spends a lot of time at the edge of the sorting network, or else $\sigma_G(i, \cdot)$ is close to a sine curve. 

\medskip

Recall that $\mathcal{D}$ is the closure of all sorting network trajectories in the uniform norm. To metrize weak convergence on the space of probability measures on $\mathcal{D}$, we use the L\'evy-Prokhorov metric $d_{LP}$. For a set $A \subset \mathcal{D}$, define 
$$
A^\epsilon = \{ f \in \mathcal{D} : ||f - g||_u < \epsilon \;\; \text{for some} \; g \in A \}.
$$
Here and throughout the remaining proofs, $||\cdot||_u$ is the uniform norm. For two probability measures $\nu_1, \nu_2$ on $\mathcal{D}$, define
$$
d_{LP}(\nu_1, \nu_2) = \inf \big\{ \epsilon > 0 : \nu_1(A) \le \nu_2(A^\epsilon) + \epsilon \text{  for all Borel sets } A \subset \mathcal{D} \big\}.
$$
We now prove two lemmas characterizing particle behaviour when $\nu_\sigma$ is close to ${\mathbb{P}_\mathcal{A}}$. The first gives conditions under which a particle spends time close to the edge of a sorting network. Let $a_{n, i} = 2i/n - 1$, and define 
$$
C_{n, i}(t) = \sin(\arcsin(a_{n, i}) + \pi t) \qquad \;\text{and}\; \qquad c_{n, i}(t) = \sin(\arcsin(a_{n, i}) - \pi t).
$$
Note that these are simply the paths $h_{a_{n, i}, \pm 1}$ in Theorem \ref{T:unique-2}. Recalling the definition of the Archimedean measure $\mathfrak{Arch}$ from Section \ref{S:intro}, for $\epsilon \in (0, 1]$ define
$$
L(\epsilon) =\mathfrak{Arch}( r \ge 1 - \epsilon, \theta \in [2\pi x, 2\pi (x + \epsilon)]) =\epsilon\sqrt{2\epsilon - \epsilon^2},
$$
where $(r, \theta)$ are polar coordinates. Note that $L(\epsilon)$ is independent of $x$ by the rotational symmetry of $\mathfrak{Arch}$.

\begin{lemma}
\label{L:edge-cond}
Let $\gamma \in (0, 2]$ and $\epsilon \in (0, \gamma^2/100)$. 
Suppose that for a fixed $n$-element sorting network $\sigma$ and a particle $i \in \{1, \mathellipsis n \}$, that $d_{LP}(\nu_\sigma, {\mathbb{P}_\mathcal{A}}) < L(\epsilon)/2$, and either
$$
\max_{t \in [0, 1]} [\sigma_G(i, t)  - C_{n,i}(t)]\ge \gamma \quad \text{or} \quad \max_{t \in [0, 1]} [c_{n, i}(t) - \sigma_G(i, t)] \ge \gamma.
$$
Then 
$$
\mathcal{L}\left\{ t: |\sigma_G(i, t)| \ge 1- 3\epsilon \right\} \ge \gamma/6.
$$
\end{lemma}

\begin{proof}
Without loss of generality, we can assume that there exists an $s \in [0, \arccos(a_{n, i})/\pi]$ such that
\begin{equation}
\label{E:sig-ga}
\sigma_G(i, s)  - C_{n, i}(s) \ge \gamma.
\end{equation}
The other cases follow by symmetric arguments.
Since $d_{LP}(\nu_\sigma, {\mathbb{P}_\mathcal{A}}) < L(\epsilon)/2$, for every $x \in [0, 1]$, there exists a particle $j(x)$ such that $$
||\sigma_G(j(x), t) - r\sin(\pi t + 2 \pi (x + \alpha)) ||_u < L(\epsilon)/2
$$
for some $r \ge 1 - \epsilon$ and $\alpha \in [0, \epsilon]$. Now since
$$
||r\sin(\pi t + 2 \pi (x + \alpha)) - \sin(\pi t + 2 \pi x)||_u \le 2 \epsilon,
$$
this implies that 
\begin{equation}
\label{E:sig-jx}
||\sigma_G(j(x), t) -  \sin(\pi t + 2 \pi x)||_u < L(\epsilon)/2 + 2 \epsilon
\end{equation}
for all $x \in [0, 1]$. Let $\delta = \gamma - \frac{L(\epsilon)}2 - 2 \epsilon$. Note that $\delta$ is positive. For all $\alpha \in [0, \delta]$, the inequality \eqref{E:sig-ga} implies that
\begin{align}
\label{E:cond-1}
\sigma_G(i, s)- \sin(\arcsin(a_{n, i}) + \alpha + \pi s) > \gamma - \delta = L(\epsilon)/2 + 2 \epsilon.
\end{align}
Now define 
$j_{\alpha} = j(\arcsin(a_{n, i}) + \alpha)/(2\pi)).$ Combining \eqref{E:cond-1} with \eqref{E:sig-jx}, we have that
\begin{align}
\nonumber
\sigma_G(i, s) > \sigma_G\left(j_\alpha, s\right).
\end{align}
 Moreover, if $1 - \cos(\alpha) > L(\epsilon)/2 + 2 \epsilon$, then 
 $$
 \sigma_G(i, 0) = a_{n, i} <  \sin(\arcsin(a_{n, i}) + \alpha) - L(\epsilon)/2 + 2 \epsilon,
 $$
 and so 
 $$
 \sigma_G(i, 0) < \sigma_G\left(j_\alpha, 0\right).
 $$

Thus for every $\alpha \in [\arccos(1 - L(\epsilon)/2 - 2 \epsilon), \delta]$, the particles $i$ and $j_\alpha$ must cross during the interval $[0, s]$. Therefore
\begin{equation}
\label{E:sig-G}
\sigma_G(i, t) > \sigma_G\left(j_\alpha, t\right)\qquad \text{ for all } t > s, \;\; \alpha \in [\arccos(1 - L(\epsilon)/2 - 2 \epsilon), \delta].
\end{equation}
 We now show that this forces the particle $i$ to spend a large amount of time close to the edge of the sorting network.  For all $\alpha \in [0, \delta]$, there must be some time $t_\alpha \in [0, 1]$ such that 
$$
\sin(\pi t_\alpha +\arcsin(a_{n, i}) + \alpha) = 1.
$$
The time  $t_\alpha \notin [0, s]$ since for every $\alpha \in [0, \delta]$, we have that
$$
||\sin(\pi t +\arcsin(a_{n, i}) + \alpha) - C_{n, i}(t)||_u < \gamma, \quad \;\text{and}\; \;\; C_{n, i}(t) \le C_{n, i}(s) \le 1 - \gamma \text{ for all } t \in [0, s].
$$
We have used \eqref{E:sig-ga} to get the second statement above. Therefore for all $\alpha \in [\arccos(1 - L(\epsilon)/2 - 2 \epsilon), \delta]$, \eqref{E:sig-G} implies that
$$
\sigma_G(i, t_\alpha) > \sigma(j_\alpha, t_\alpha) \ge 1 - L(\epsilon)/2 - 2\epsilon \ge 1 - 3\epsilon.
$$ 
 Using the fact that $\arccos(1 - x) \le 2\sqrt{x}$ for $x \in [0, \pi/2)$, we have that
 \[
\mathcal{L} \big\{t_\alpha : \alpha \in [\arccos(1 - L(\epsilon)/2 - 2 \epsilon), \delta] \big\} \ge \frac{\delta - 2\sqrt{L(\epsilon)/2 + 2 \epsilon}}{\pi} \ge \frac{\gamma - 3\epsilon - 2\sqrt{3\epsilon}}{\pi} \ge \frac{\gamma}{6}.
\]
Here the final bound follows from the fact that $\epsilon < \gamma^2/100$.
\end{proof}

The second lemma shows that if a curve $\sigma_G(i, \cdot)$ stays close to the region between the curves $c_{n, i}$ and $C_{n, i}$, then it must be close to a sine curve. 


\begin{lemma}
\label{L:sine-close}
Let $\sigma$ be an $n$-element sorting network, $i \in \{1, \mathellipsis n \}$, and $\gamma \in (0, 1)$. Suppose that
$$
\max_{t \in [0, 1]} [\sigma_G(i, t)  - C_{n,i}(t)] < \gamma \quad \text{and} \quad \max_{t \in [0, 1]} [c_{n, i}(t) - \sigma_G(i, t)] < \gamma,
$$
and that $d_{LP}(\nu_\sigma, {\mathbb{P}_\mathcal{A}}) < \frac{\gamma^4}{128}$. Then there exist constants $a \in [0, 1]$ and $\theta \in [0, 2\pi]$ such that 
$$
||\sigma_G(i, t) - a\sin(\pi t + \theta) ||_u < 2\gamma + 2/n.
$$
\end{lemma}

\begin{proof} 
Suppose first that for some $s \in [0,1]$, that 
\begin{equation}
\label{E:sig-GG}
\sigma_G(i, s) \in [c_{n, i}(s) + \gamma, C_{n, i}(s) - \gamma]. 
\end{equation}
Observe that $s \in [\arcsin(\gamma)/\pi, 1- \arcsin(\gamma)/\pi]$, since otherwise the above interval 
is empty. By \eqref{E:sig-GG}, we can find a point $(x_0, y_0) \in B(0, 1 - \gamma)$ such that 
$$
x_0 = \sigma_G(i, 0) \qquad \;\text{and}\; \qquad x_0\cos(\pi s) + y_0 \sin(\pi s) = \sigma_G(i, s).
$$
Since $s \le 1- \arcsin(\gamma)/\pi$, we can find another point 
$
(x_1, y_1) \in
B((x_0, y_0), \gamma) \subset B(0, 1)
$
such that 
 \begin{equation}
 \label{E:x1-x0}
x_1 > x_0 + \gamma^2/3 \quad \;\text{and}\; \quad x_1\cos(\pi s) + y_1\sin(\pi s) > x_0\cos(\pi s) + y_0\sin(\pi s) + \gamma^2/3.
 \end{equation}
Now observe that
$$
\gamma^4/64\le \inf \left\{ \mathfrak{Arch}(B(x, \gamma^2/6)) : x \in B(0, 1) \right\}.
$$
Therefore since $d_{LP}(\nu_\sigma, \mathbb{P}_{\mathbb{P}_\mathcal{A}}) < \gamma^4/128$, there must be a point $(x_2, y_2) \in  B(0, 1) \cap B((x_1, y_1),\gamma^2/6)$ and a particle $j \in \{1, \dots, n\}$ such that 
\begin{equation}
\label{E:G-jj}
\begin{split}
||\sigma_G(j, t) - (x_1\cos(\pi t) &+ y_1\sin(\pi t))||_u \\
&\le ||\sigma_G(j, t) - (x_2\cos(\pi t) + y_2\sin(\pi t))||_u + \gamma^2/6 < \gamma^2/3.
\end{split}
\end{equation}
By \eqref{E:x1-x0}, this implies that 
$$
\sigma_G(j, 0) > \sigma_G(i, 0) \qquad \;\text{and}\; \qquad \sigma_G(j, s) > \sigma_G(i, s),
$$
and hence $\sigma_G(j, t) > \sigma_G(i, t)$ for all $t \in [0, s]$. Combining this with \eqref{E:G-jj} and the fact that $(x_1, y_1) \in B((x_0, y_0), \gamma)$ implies that
$$
x_0\sin(\pi t) + y_0 \cos(\pi t) - \sigma_G(i, t)  \le \gamma^2/3 + \gamma < 2\gamma \qquad \text{ for all } t \in [0, s].
$$
Symmetric arguments give the same upper bound on this difference over the interval $[s, 1]$, and on the difference $\sigma_G(i, t) - x_0\sin(\pi t) + y_0 \cos(\pi t)$ over the interval $[0, 1]$. This implies that
$$
||\sigma_G(i, t) - [x_0\sin(\pi t) + y_0 \cos(\pi t)]||_u \le 2\gamma.
$$
 Now suppose that there does not exist an $s \in [0, 1]$ such that $\sigma_G(i, s) \in [c_{n, i}(s) + \gamma, C_{n, i}(s) - \gamma]$. Then either
$$
||\sigma_G(i, \cdot) - C_{n, i}(\cdot)||_u < 2\gamma + \frac{2}n \qquad \;\text{or}\; \qquad ||\sigma_G(i, \cdot) - c_{n, i}(\cdot)||_u < 2\gamma + \frac{2}n.
$$
This can be seen by simply analyzing the paths $C_{n, i}$ and $c_{n, i}$, noting that the trajectory $\sigma_G(i, \cdot)$ can ``jump" distances of size $2/n$. 

\medskip
As all the functions $C_{n, i}$, $c_{n, i}$, and $x_0\sin(\pi t) + y_0 \cos(\pi t)$ are of the form $a\sin(\pi t + \theta)$, for some $(a, \theta) \in [0, 1] \times[0, 2\pi]$, this proves the lemma.
\end{proof}

We can now put together all the ingredients to prove Theorem \ref{T:sine-curves}.

\begin{proof} Fix $\gamma \in (0, 1)$, and recall from Section \ref{SS:slightly} that $\nu_n$ is uniform measure on the set $\{\nu_\sigma\}_{\sigma \in \Omega_n}$. For small enough $\epsilon > 0$, Lemma \ref{L:edge-time} and Theorem \ref{T:weak-limit'} imply that 
$$
\mathbb{P} \left ( \max_{i \in [1, n]} \mathcal{L}\left\{t: |\sigma^n_G(i, t)| \ge 1 -3\epsilon \right\} < \frac{\gamma}{6}, \;\; d_{LP}(\nu_n, \mathbb{P}_\mathcal{A}) < \frac{\gamma^4}{128} \wedge \frac{L(\epsilon)}2 \right) \to 1 \qquad \;\text{as}\; n \to \infty.
$$
Combining Lemmas \ref{L:edge-cond} and \ref{L:sine-close}, this implies that there exist random variables $A_{n, i, \gamma} \in [0, 1]$ and $\Theta_{n, i, \gamma} \in [0, 2\pi]$ such that 
$$
P_{n, \gamma} := \mathbb{P} \left(\max_{i \in [1, n]} ||\sigma^n_G(i, t) - A_{n, i, \gamma} \sin(\pi t + \Theta_{n, i, \gamma})||_u < 2\gamma + \frac{1}n \right) \to 1 \qquad \;\text{as}\; n \to \infty.
$$
As constructed, the random variables $A_{n, i, \gamma}$  and $\Theta_{n, i, \gamma}$ depend on $\gamma$. To remove this dependence, let $\gamma_n \to 0$ be a sequence such that $P_{n, \gamma_n} \to 1$ as $n \to \infty$. Let $A_{n, i} = A_{n, i, \gamma_n}$ and $\Theta_{n, i} = \Theta_{n, i, \gamma_n}$, 
and define
\begin{equation}
\label{E:Bnga}
B_{n, \gamma} = \left\{ \max_{i \in [1, n]} ||\sigma^n_G(i, t)- A_{n, i} \sin(\pi t + \Theta_{n, i})||_u < \gamma \right\},
\end{equation}
Then for any $\gamma > 0$, we have that $\mathbb{P}(B_{n, \gamma}) \to 1$ as $n \to \infty$.
\end{proof}

We can also prove Theorem \ref{T:matrices} and Theorem \ref{T:unif-rotation}. For these we use the notation $B_{n, \gamma}$ from \eqref{E:Bnga}.

\begin{proof}[Proof of Theorem \ref{T:matrices}] Fix $t \in [0, 1]$.
Since $\nu_{n} \to \delta_{\mathbb{P}_\mathcal{A}}$ (Theorem \ref{T:weak-limit'}), the random measure $\rho_t^n$ converges in probability to the law of $(X, X \cos (\pi t )+ Y \sin (\pi t))$, where $(X, Y) \sim \mathfrak{Arch}$. This law is simply $\mathfrak{Arch}_t$. 

\medskip

Now on the event $B_{n, \gamma}$, the support of the measure $\rho_t^n$ is contained in the set 
$$
\text{supp}(\mathfrak{Arch}_t)^\gamma = \{ x \in [-1, 1]^2 : d(x, \text{supp}(\mathfrak{Arch}_t)) < \gamma \}.
$$
Moreover, since $\mathfrak{Arch}_t$ has a Lebesgue density that is bounded below, the weak convergence of $\rho^n_t$ to $\rho$ implies that with high probability,
$
\text{supp}(\mathfrak{Arch}_t) \subset \text{supp}(\rho_n^t)^\gamma
$
for any $\gamma > 0$. Therefore since $\mathbb{P}(B_{n, \gamma}) \to 1$ as $n \to \infty$ by Theorem \ref{T:sine-curves}, we have that
\[
\lim_{n \to \infty} d_H(\text{supp}(\rho^n_t), \text{supp}(\mathfrak{Arch}_t)) = 0 \qquad \text{ in probability. }  \quad \qedhere
\]
\end{proof}

\begin{proof}[Proof of Theorem \ref{T:unif-rotation}]
First observe that for any $a \in [0, 1]$, $t \in [0, 1/2]$, and $\theta \in [0, 2\pi]$, that
$$
e^{\pi i t} \left[a\sin(\pi t + \theta) + ia\sin(\pi t + \pi/2 + \theta) \right] = a \sin(\theta).
$$
Therefore on the event $B_{n, \gamma}$, we have that
$$
\max_{j \in [1, n]} \max_{s, t \in [0, 1]} |Z^n_j(t) - Z^n_j(s)| \le 2\gamma.
$$
Since $\mathbb{P}(B_{n, \gamma}) \to 1$ as $n \to \infty$ for any $\gamma > 0$ by Theorem \ref{T:sine-curves}, this proves the theorem.
\end{proof}



\end{section}
\section{The geometric limit}
\label{S:geom-limit}

In this section, we use Theorem \ref{T:sine-curves} and Theorem \ref{T:weak-limit} to prove Theorem \ref{T:geom-limit}. Recall from Section \ref{S:intro} that $d_\infty(f, g)$ is the uniform norm between two $\mathbb{R}^n$-valued functions $f$ and $g$, where the pointwise distance is the $L^\infty$-distance. Recall also that $\bar{\sigma}$ is the embedding of a sorting network $\sigma$ into the $(n-2)$-dimensional sphere $\mathbb{S}^{n-2} \subset \mathbb{R}^n$.

\medskip

By Theorem \ref{T:sine-curves} and a change of variables,  there exist random vectors $\mathbf{X^n} = (X^n_1, \mathellipsis X^n_n)$ and $\mathbf{V^n} = (V^n_1, \mathellipsis V^n_n)$ such that $d_\infty(f_n, \close{\sigma}^n)/n \to 0$ in probability, where 
$$
f_n(t) = \mathbf{X^n} \cos(\pi t) + \mathbf{V^n} \sin(\pi t) + \mathbf{c}, \qquad \text{where} \qquad 
\mathbf{c} = \left( \frac{n + 1}2, \mathellipsis , \frac{n + 1}2\right).
$$
Moreover, we can assume that $X^n_i = i - (n + 1)/2$ and that $V^n_i =  \sigma^n(i, N/2) -  (n + 1)/2$, as these changes only shift the curve  $\mathbf{X^n} \cos(\pi t) + \mathbf{V^n} \sin(\pi t)$ by $d_\infty$-distance $o(n)$ in probability. 

\medskip

The point $\mathbf{c}$ is the center of $\mathbb{S}^{n-2}$. It remains to show that we can shift the curve $f_n$ by $d_\infty$-distance $o(n)$ to obtain a great circle in $\mathbb{S}^{n-2}$. For this we need the following lemma.

\begin{lemma}
\label{L:dot-prod-small}
Let the vectors $\mathbf{X^n}$ and $\mathbf{V^n}$ be as above. Then 
$$
\lim_{n \to \infty} \frac{\left\langle\mathbf{X^n} , \mathbf{V^n} \right\rangle}{n^3} = 0 \qquad \text{in probability.}
$$
\end{lemma}

\begin{proof}
Let $I_n$ be a uniform random variable on $\{1, \dots, n\}$, independent of $\mathbf{X^n}$ and $\mathbf{V^n}$, and define
$$
(\tilde{X}^n, \tilde{V}^n) =  \left(\frac{2X^n_{I_n}}n, \frac{2V^n_{I_n} }n\right).
$$
We have that 
$$
(\tilde{X}^n, \tilde{V}^n) \stackrel{d}{=} (Y_n(0) - 1/n, Y_n(1/2) - 1/n),
$$
where $Y_n$ is the trajectory random variable of $\sigma^n$.  Therefore by Theorem \ref{T:weak-limit}, 
$$
(\tilde{X}^n, \tilde{V}^n) \stackrel{d}{\to} (X, V) \sim \mathfrak{Arch}.
$$
By the bounded convergence theorem, this implies that $\mathbb{E} \tilde{X}^n\tilde{V}^n \to \mathbb{E} XV$ in probability. Observing that $n^3 \mathbb{E} \tilde{X}^n\tilde{V}^n = 4\left\langle\mathbf{X^n} , \mathbf{V^n} \right\rangle$ and that $\mathbb{E} XV = 0$ completes the proof.
\end{proof}

\begin{proof}[Proof of Theorem \ref{T:geom-limit}.] Fix $n \in \mathbb{N}$. For ease of bounding errors, we assume that $n \ge 27$. Let $\sigma$ be a fixed $n$-element sorting network. Our goal is to perturb the vector $\mathbf{V^n}$ to a new vector $\mathbf{W^n}$ so that the path
\begin{equation}
\label{E:C-nep}
C_{n}(t) := \mathbf{X^n} \cos(\pi t) + \mathbf{W^n} \sin(\pi t) + \mathbf{c}
\end{equation}
follows a great circle on $\mathbb{S}^{n-2}$, and so that $d_\infty(C_{n}, f_n)$ is small. In order to do this, we just need to find $\mathbf{W^n}$ such that
$$
\sum_{i=1}^n W^n_i = 0, \;\left\langle\mathbf{X^n} , \mathbf{W^n} \right\rangle = 0, \quad \;\text{and}\; \quad ||W^n||_2^2 = (n^3 - n)/12.
$$
The first of these conditions guarantees that $C_n$ lies in the correct hyperplane $\mathbb{L}_n \subset \mathbb{R}^n$, and the second and third conditions guarantee that $C_n$ traces out a great circle on the sphere $\mathbb{S}^{n-2}$ within that hyperplane. We first perturb $\mathbf{V^n}$ to a vector $\mathbf{Z^n}$ satisfying the first two properties. Define the random variable
$$
A_n = \frac{\left|\left\langle\mathbf{X^n} , \mathbf{V^n} \right\rangle\right|}{n^3}.
$$
Without loss of generality, assume that $\left\langle\mathbf{X^n} , \mathbf{V^n} \right\rangle > 0$. 

\medskip

Choose any $i \ge 3(n + 1)/4$. Let $k > 0$, and decrease the value of $V^n_i$ by $kn$ and increase the value of $V^n_{n + 1 - i}$ by $kn$. Call this new vector $\mathbf{V^n_1}$. Note that $X^n_i \ge n/4$, and that $X^n_i = -X^n_{n + 1 - i}$. Therefore $\left\langle\mathbf{X^n} , \mathbf{V^n_1} \right\rangle \le\left\langle\mathbf{X^n} , \mathbf{V^n} \right\rangle - kn^2/2$.
  
 \medskip
 
By iterating the above procedure to repeatedly lower the dot product $\left\langle\mathbf{X^n} , \mathbf{V^n} \right\rangle$, we can obtain a vector $\mathbf{K_n} = (K_1n, \mathellipsis K_nn)$ with the following properties (here is where we require that $n \ge 27$).
\smallskip
\begin{enumerate}[nosep,label=(\roman*)]
\item For all $i$, $|K_i| \le 9A_n$, $- K_i = K_{n + 1 - i}$, and $k_i = 0$ if $i \in \left(\frac{n+1}4, \frac{3(n+1)}4\right)$.
\item $\left\langle\mathbf{X^n} , \mathbf{V^n + K^n} \right\rangle = 0.$
\end{enumerate}
\smallskip

Let $\mathbf{Z^n} = \mathbf{V^n + K^n}$. Observe that $\sum_{i=1}^n Z^n_i = 0$ since both $\sum_{i=1}^n K_i = 0$ and $\sum_{i=1}^n V^n_i = 0$. The fact that $\sum_{i=1}^n V^n_i = 0$ follows since $\mathbf{V^n} + \mathbf{c}$ is a permutation of the vector $(1, 2, \mathellipsis n)$.
For any $t$, we have that 
\begin{align*}
||f_n(t) - (\mathbf{X^n} \cos(\pi t) + \mathbf{Z^n} \sin(\pi t) + \mathbf{c})||_\infty = \max_{i \in [1, n]} |V^n_i \sin(\pi t) - Z^n_i \sin(\pi t)| \le 9 A_n n.
\end{align*}
We can now define $\mathbf{W^n} = M_n\mathbf{Z^n}$, where $M_n$ is a random constant chosen so that $||\mathbf{W^n}||^2_2 = (n^3 - n)/12$. Using the definition of $C_n$ in \eqref{E:C-nep}, for every $t \in [0, 1]$, we have that
 \begin{equation}
 \label{E:fC}
 \begin{split}
||f_n(t) - C_n(t)||_\infty &\le 9\epsilon n +  \max_{i \in [1, n]} |W^n_i \sin(\pi t) - Z^n_i \sin(\pi t)| \\
&< 9 A_n n + |M_n - 1|\left[\frac{n + 1}2 + 9 A_n n\right].
\end{split}
\end{equation}
In the last inequality, we have used that $|Z^n_i - V^n_i| \le 9 A_n n$. The inequality \eqref{E:fC} implies that
\begin{equation}
\label{E:d-inf-comp}
\frac{d_\infty(C_n, \close{\sigma}^n)}n \le 9 A_n + |M_n - 1|\left[\frac{n + 1}{2n} + 9 A_n \right] + \frac{d_\infty(f_n, \close{\sigma}^n)}n.
\end{equation}
Now again using that $|Z^n_i - V^n_i| \le 9 A_n n$, we have that  
\begin{align*}
\big|\;||\mathbf{Z^n}||^2_2 - ||\mathbf{V^n}||^2_2 \; \big| \le \sum_{i=1}^n |Z^n_i - V^n_i||Z^n_i + V^n_i| \le 9A_n n^2(n + 1 + 9 A_n n).
\end{align*}
Therefore since $||\mathbf{V^n}||^2_2 = (n^3 - n)/12$, we have that
$$
M_n \in  \left[\sqrt{1 - 217A_n}, \sqrt{1 + 217A_n} \right]
$$
whenever $A_n < 1 /217$. Lemma \ref{L:dot-prod-small} implies that $A_n \to 0$ as $n \to \infty$ in probability, and hence so does $M_n$. Therefore since $d_\infty(f_n, \close{\sigma}^n)/n \to 0$ in probability as $n \to \infty$, \eqref{E:d-inf-comp} implies that $d_\infty(C_n, \close{\sigma}^n)/n$ does as well.
\end{proof}

\section*{Acknowledgements} The author would like to thank B\'alint Vir\'ag for many fruitful discussions about the problem, and for many constructive comments about previous drafts.


\bibliographystyle{alpha}

2018-03-20T01:13:07.000
1802
1802.08934
en
https://arxiv.org/abs/1802.08934
A sorting network (also known as a reduced decomposition of the reverse permutation), is a shortest path from $12 \cdots n$ to $n \cdots 21$ in the Cayley graph of the symmetric group $S_n$ generated by adjacent transpositions. We prove that in a uniform random $n$-element sorting network $\sigma^n$, all particle trajectories are close to sine curves with high probability. We also find the weak limit of the time-$t$ permutation matrix measures of $\sigma^n$. As a corollary of these results, we show that if $S_n$ is embedded into $\mathbb{R}^n$ via the map $\tau \mapsto (\tau(1), \tau(2), \dots \tau(n))$, then with high probability, the path $\sigma^n$ is close to a great circle on a particular $(n-2)$-dimensional sphere in $\mathbb{R}^n$. These results prove conjectures of Angel, Holroyd, Romik, and Virag.
Probability (math.PR); Combinatorics (math.CO)
The Archimedean limit of random sorting networks
Qwen/Qwen-72B
1. YES
2. YES