On the number of spanning trees in random regular graphs
Let $d \geq 3$ be a fixed integer. We give an asympotic formula for the expected number of spanning trees in a uniformly random $d$-regular graph with $n$ vertices. (The asymptotics are as $n\to\infty$, restricted to even $n$ if $d$ is odd.) We also obtain the asymptotic distribution of the number of spanning trees in a uniformly random cubic graph, and conjecture that the corresponding result holds for arbitrary (fixed) $d$. Numerical evidence is presented which supports our conjecture.
\section{Introduction}\label{s:intro}


\global\long\def\zeta{\zeta}

In this paper, $d$ denotes a fixed integer which is at least 2 (and usually
at least 3). All asymptotics are taken as $n\to\infty$, with $n$ restricted to 
even integers when $d$ is odd.

The number of spanning trees in a graph, also called the \emph{complexity}
of the graph,
is of interest for a number of reasons.  The complexity of a graph
is an evaluation of
the Tutte polynomial (see for example~\cite{welsh}).  
The Merino-Welsh conjecture~\cite{MW} relates
the complexity of a graph with two other graph
parameters, namely, the number of acyclic orientations and the number of
totally cyclic orientations of a graph. 
(Noble and Royle~\cite{NR} recently proved that the Merino-Welsh conjecture
is true for series-parallel graphs.) The complexity of a graph also plays
a role in the theory of electrical networks (see for example~\cite{myers}).


We are interested in the number of spanning trees in random regular graphs.
The first significant result in this area is due to McKay~\cite{McKay-0}, 
who proved that for $d\ge3$, the
$n$th root of the number of spanning trees of a random $d$-regular
graph with $n$ vertices converges to 
\begin{equation}
\frac{\left(d-1\right)^{d-1}}{\left(d^{2}-2d\right)^{d/2-1}}\label{eq:mckay}
\end{equation}
as $n\to\infty$, with probability one. An alternative proof of this
was later given
by Lyons~\cite[Example 3.16]{Lyons}.

McKay~\cite[Theorem 4.2]{McKay} gave an asymptotic expression for the
expected
number of spanning trees in a random graph with specified degrees,
up to some unknown constant.  His result holds when the maximum degree
is bounded 
and the average degree is bounded away from 2 (independently of $n$).
When specialised to regular degree sequences,~\cite[Theorem 4.2]{McKay}
states that the
expected number of spanning trees in $\mathcal{G}_{n,\k}$ is asymptotic to 
\begin{equation}
\frac{c_{d}}{n}\left(\frac{\left(d-1\right)^{d-1}}{\left(d^{2}-2d\right)^{d/2-1}}\right)^{n},\label{eq:cd}
\end{equation}
for some unknown constant $c_{d}$.

Other work on asymptotics for the number of spanning trees has focussed on
circulant graphs, grid graphs and tori (see for example~\cite{GYZ} and the
references therein).

Our first result, Theorem~\ref{thm:expectation}, provides the value of the 
constant $c_{d}$ from (\ref{eq:cd}), proving that
\[
c_{d}=
  \exp\left(\frac{6d^{2}-14d+7}{4\left(d-1\right)^{2}}\right)\,
\frac{\left(d-1\right)^{1/2}}{\left(d-2\right)^{3/2}}.
\]
For our second result we investigate
the distribution of the number of spanning trees in random $d$-regular
graphs using the small subgraph conditioning method, and obtain
the asymptotic distribution in the case of cubic graphs, presented in
Theorem~\ref{thm:distribution-3}.  We provide partial calculations for 
arbitrary fixed degrees, which lead us to conjecture that the corresponding
result holds in general (see Conjecture~\ref{cnj:distribution}).

In order to precisely state our main results we must introduce some notation
and terminology.

\subsection{Notation and our main results}\label{ss:notation}

Let $\mathbb{N}$ denote the natural numbers (which includes 0).
For integers $n$, $k$ let $(n)_{k}$ denote the falling factorial
$n(n-1)\cdots (n-k+1)$. Square brackets without subscripts denote
extraction of coefficients of a generating function.
We use $\boldsymbol{1}\left(\cdot \right)$ to denote both the indicator variable of
an event and the characteristic function of a set (the particular set will
appear as a subscript). We use standard asymptotic notation throughout,
with the exception that $\rightsquigarrow$ indicates convergence in
distribution of a sequence of random variables.
(We use this notation rather than $\stackrel{d}{\rightarrow}$
to avoid overloading the symbol ``$d$'',
which we use for the degree of the graph.)


Let $\mathcal{G}_{n,\k}$ denote the uniform model of $d$-regular simple graphs 
on the vertex set $\{ 1,\ldots, n\}$. 
Define the random variable $Y_{\mathcal{G}}$ to be
the number of spanning trees in a random $G\in\mathcal{G}_{n,\k}$. 

Clearly $Y_{\mathcal{G}}$ is identically zero if $n\geq 3$ and
$d < 2$.  
A 2-regular graph has a spanning tree if and only if
it is connected (that is, forms a Hamilton cycle), in which case it has exactly $n$ spanning
trees. Hence the distribution of $Y_{\mathcal{G}}$
can be inferred from~\cite[Equation (11)]{Wormald}.
For the remainder of the paper we assume that $d\geq 3$.

Our first result gives an asymptotic expression for the expectation
of $Y_{\mathcal{G}}$.

\begin{theorem}
\label{thm:expectation}
Let $d\geq 3$ be a fixed integer.  Then
\[
\mathbb{E} Y_{\mathcal{G}} \sim
    \exp\left(\frac{6d^{2}-14d+7}{4\left(d-1\right)^{2}}\right)\,
\frac{\left(d-1\right)^{1/2}}{n\left(d-2\right)^{3/2}}\,
\left(\frac{\left(d-1\right)^{d-1}}{\left(d^{2}-2d\right)^{d/2-1}}
   \right)^{n}\ .
\]
\end{theorem}

This theorem is proved at the end of Section~\ref{sec:E_YX}.

Next, for fixed $d\geq 3$ and for each positive integer $j$, define
\begin{equation}
\lambda_{j}\left(d\right)  =  \frac{\left(d-1\right)^{j}}{2j},\qquad
\zeta_{j}\left(d\right)  = 
  -\frac{2\left(d-1\right)^{j}-1}{\left(d-1\right)^{2j}}
\label{eq:janson-parameters}
\end{equation}

Our second theorem gives the asymptotic distribution of the number of
spanning trees in the case of cubic graphs.

\begin{theorem}
\label{thm:distribution-3}
Let $Z_{j}\sim\mathrm{Poisson}\left(\lambda_{j}\left(3\right)\right)$,
with each $Z_{j}$ independent. Consider the number of spanning trees
in a random cubic graph, normalized by the expectation given in Theorem
\ref{thm:expectation} for $d=3$. The asymptotic distribution of
this quantity is given by
\[
\frac{Y_\mathcal{G}}{\mathbb{E} Y_{\mathcal{G}}} \,\,
\rightsquigarrow\,\,
\prod_{j=3}^{\infty}\, \left(1+\zeta_{j}\left(3\right)\right)^{Z_{j}}
\, e^{-\lambda_{j}\left(3\right)\zeta_{j}\left(3\right)}.
\]
\end{theorem}

This theorem is proved in Section~\ref{sec:EY2}.
We conjecture that an analogous result holds for arbitrary (fixed) $d \geq 3$.

\begin{conjecture}
\label{cnj:distribution}
Let $d\geq 3$ be fixed.  Then
\[
\frac{Y_\mathcal{G}}{\mathbb{E} Y_{\mathcal{G}}} \,\,
\rightsquigarrow
\,\, \prod_{j=3}^{\infty}\,
   \left(1+\zeta_{j}\left(d\right)\right)^{Z_{j}}\,
   e^{-\lambda_{j}\left(d\right)\zeta_{j}\left(d\right)},
\]
where $Z_{j}\sim\mathrm{Poisson}\left(\lambda_{j}\left(d\right)\right)$
and each $Z_{j}$ is independent.
\end{conjecture}

We present numerical evidence which supports this conjecture in
Section~\ref{ss:numerical}.

\subsection{Plan of attack}\label{ss:plan}

From now on, we omit explicit mention
of $d$ in the constants $\zeta_{j}=\zeta_{j}\left(d\right)$ and
$\lambda_{j}=\lambda_{j}\left(d\right)$ from (\ref{eq:janson-parameters}).

As is standard in this area, most of our calculations
will be performed in the uniform probability space $\mathcal{P}_{n,\k}$ of
\textit{pairings} (also called the 
\textit{configuration model}~\cite{Bollobas,McKay-pairings,Wormald}.
Let $d$ and $n$ be positive integers such that $d n$ is even.
Consider a set of $d n$ \emph{prevertices} distributed evenly into 
$n$ sets, called
\emph{buckets}.  
(We prefer the terminology ``prevertices'' to ``points''.)  
A \emph{pairing} is a partition of the prevertices into
$d n/2$ sets of size 2, called \emph{pairs}.  
Then
\begin{equation}
|\mathcal{P}_{n,\k}| = \#P\left(d n\right) = \frac{(d n)!}{(d n/2)!\, 2^{d n/2}} 
  \sim\sqrt{2}\left(\frac{d n}{e}\right)^{d n/2},
\label{eq:num-pairings-asymptotic}
\end{equation}
using Stirling's formula.

By contracting the prevertices in each bucket to a vertex, each pairing
projects to a labelled $d$-regular multigraph, with loops permitted.
Let $\Omega_{n,\k}$ denote the set of such multigraphs, and denote the projection
of a pairing $P$ by $G\left(P\right)$.
(We will occasionally informally refer to ``partial'' pairings, where
only a subset of the prevertices are paired. The projection of a partial
pairing is defined in the same way.)

Each $G\in\mathcal{G}_{n,\k}$ is the projection of $\left(d!\right)^{n}$ different
pairings (permuting the prevertices in each bucket), so we can recover
the uniform model $\mathcal{G}_{n,\k}$ from $\mathcal{P}_{n,\k}$ by conditioning on the event that
the projected multigraph of a random pairing is simple.

We will apply the small subgraph conditioning method in the
form given by Janson~\cite[Theorem 1]{Janson}.

\begin{theorem}
\label{thm:janson}Let $\lambda_{j}>0$ and $\zeta_{j}\ge-1$, $j=1,2,\dots,$
be constants. Suppose that for each $n$ we have a sequence 
$\boldsymbol{X} =(X_1,X_2,X_3,\ldots)$ of
non-negative integer valued random variables and a random variable
$Y$ with $\mathbb{E} Y\ne0$ (at least for large $n$). Further suppose the
following conditions are satisfied:

\begin{enumerate}[label=\textup{(A\arabic*)}]

\item\label{cond:A1}
For $m\ge1$, $\left(X_{1},\dots,X_{m}\right)\rightsquigarrow\left(Z_{1},\dots,Z_{m}\right)$,
where $Z_{j}$ are independent Poisson random variables with means
$\lambda_{j}$; 

\item\label{cond:A2}
For any $m\ge0$, $\delta^{1}\in\mathbb{N}^{m}$,
\[
\frac{\mathbb{E}\left[Y|X_{1}=\delta^{1}_{1},\dots,X_{m}=\delta^{1}_{m}\right]}{\mathbb{E} Y} \longrightarrow \prod_{j=1}^{m}\left(1+\zeta_{j}\right)^{\delta^{1}_{j}}e^{-\lambda_{j}\zeta_{j}};
\]


\item\label{cond:A3}
${\displaystyle \sum_{j=1}^{\infty}\lambda_{j}\zeta_{j}^{2}<\infty}$;

\item\label{cond:A4}
$
{\displaystyle \frac{\mathbb{E} Y^{2}}{\left(\mathbb{E} Y\right)^{2}}\rightarrow\exp\left(\sum_{j=1}^{\infty}\lambda_{j}\zeta_{j}^{2}\right)}$.\end{enumerate}
Then
\begin{equation}
\frac{Y}{\mathbb{E} Y}\rightsquigarrow W:=\prod_{j=1}^{\infty}\left(1+\zeta_{j}\right)^{Z_{j}}e^{-\lambda_{j}\zeta_{j}}.\label{eq:W}
\end{equation}
Moreover, this and the convergence in \ref{cond:A1} hold jointly.
\end{theorem}

We also need a related lemma:

\begin{lemma}
\label{lem:janson}\textup{(\cite[Lemma 1]{Janson})} Let $\lambda_{j}'>0$,
$j=1,2,\dots,$ be constants. Suppose that \ref{cond:A1} holds and
that $Y\ge0$.\textup{ S}uppose:

\begin{enumerate}[label=\textup{(A2$'$)}]

\item\label{cond:A2'}
For any $m\ge0$, $\delta^{1}\in\mathbb{N}^{m}$,
\[
\frac{\mathbb{E}\left[Y\prod_{j=1}^{m}\falling{X_{j,n}}{\delta^{1}_{j}}\right]}{\mathbb{E} Y}\longrightarrow\prod_{j=1}^{m}\left(\lambda_{j}'\right)^{\delta^{1}_{j}}.
\]
\end{enumerate}
Then \ref{cond:A2} holds with $\lambda_{j}\left(1+\zeta_{j}\right)=\lambda_{j}'$
for all positive integers $j$.
\end{lemma}

We now define the random variables $X_{j}$ and $Y$ to which these
results will be applied.

For each $j\geq 1$, let $\gamma_{j}:\Omega_{n,\k}\to\mathbb{N}$ give the number
of cycles of length $j$ in a multigraph. (A loop is a 1-cycle, and a pair
of edges on the same two vertices is a 2-cycle.)
Then the random variable $X_{j}=\gamma_{j}\circ G$
is the number of $j$-cycles in the projection of a random pairing $P\in\mathcal{P}_{n,\k}$. Write
$\boldsymbol{X} = \left(X_{j}\right)_{j\geq 1}$ for the
sequence of all cycle counts.
It is well known~\cite{Bollobas} that for any positive integer $m$,
the random variables $X_{1},\ldots,X_{m}$ 
are asymptotically independent Poisson random variables, and that
the mean of $X_{j}$ tends to the quantity $\lambda_{j} = \lambda_{j}(d)$ 
given in (\ref{eq:janson-parameters}).
Hence Condition~\ref{cond:A1} of Theorem~\ref{thm:janson} holds.

Let $\tau:\Omega_{n,\k}\to\mathbb{N}$ be the function which counts 
spanning trees in $d$-regular multigraphs. Define $Y_{\mathcal{G}}$
as the restriction of $\tau$ to $\mathcal{G}_{n,\k}$, and define $Y=\tau\circ G$.
Then, $Y_{\mathcal{G}}$ is the number of spanning trees in a random
$G\in\mathcal{G}_{n,\k}$, as in Section~\ref{ss:notation}. $Y$ is accordingly
the number of spanning trees in the projection of a random pairing 
$P\in\mathcal{P}_{n,\k}$. We will
investigate the asymptotic distribution of $Y_{\mathcal{G}}$
through analysis of $Y$. 

In Section~\ref{sec:EY} we obtain an asymptotic formula for the expected
value of $Y$.  In Section~\ref{sec:E_YX} we analyse the interaction of the
number of spanning trees with short cycles, establishing that \ref{cond:A2}
holds for $\lambda_{j}$ and $\zeta_{j}$ as given in (\ref{eq:janson-parameters}).
This enables us to prove Theorem~\ref{thm:expectation} and to prove
that \ref{cond:A3} holds.  Then in Section~\ref{sec:EY2} we investigate
the second moment of $Y$.  We can prove that Condition~\ref{cond:A4} holds
when $d=3$, leading to a proof of Theorem~\ref{thm:distribution-3}.
Using our partial calculations for general degrees, we provide numerical
evidence that strongly supports Conjecture~\ref{cnj:distribution}.



\section{Expected number of spanning trees }\label{sec:EY}

In this section we compute $\mathbb{E} Y$. Let $\mathcal{T}_{n}$ denote the set of labelled
trees on $n$ vertices, so that $\left|\mathcal{T}_{n}\right|=n^{n-2}$ by Cayley's
formula (see for example, \cite{Moon}).

Recalling the definition of $Y$, 
we have
\[
\mathbb{E}Y = \sum_{P\in\mathcal{P}_{n,\k}}\frac{1}{\left|\mathcal{P}_{n,\k}\right|}\,\tau\left(G\left(P\right)\right),
\]
and hence
\begin{equation}
\left|\mathcal{P}_{n,\k}\right|\:\mathbb{E}Y = \sum_{P\in\mathcal{P}_{n,\k}}\sum_{T\in\mathcal{T}_{n}}M_{T,P},\label{eq:|P|EY-sum-0}
\end{equation}
where $M_{T,P}$ is the number of ways to embed $T$ into the multigraph
$G\left(P\right)$. (When $G\left(P\right)$ is simple, $M_{T,P}$
is zero or one).

Now, we want to condition on the degree of each of the $n$ vertices
in $T$. Define the set of possible degree sequences 
\[
\mathcal{D}_{n}=\left\{ \delta\in\mathbb{N}^{n}:\quad\sum_{j=1}^{n}\delta_{j}=2\left(n-1\right)\right\} .
\]
We can decompose $\left|\mathcal{P}_{n,\k}\right|\:\mathbb{E}Y$ as
\begin{equation}
\sum_{\delta\in\mathcal{D}_{n}}\sum_{T\in\mathcal{T}_{n}}\boldsymbol{1}\left(T\sim\delta\right)\,
  \sum_{P\in\mathcal{P}_{n,\k}}\, M_{T,P},\label{eq:|P|EY-sum-1}
\end{equation}
where $T\sim\delta$ denotes the event that vertex $j$ has degree $\delta_{j}$
in $T$, for all $j=1,\ldots n$.

To evaluate the innermost sum in (\ref{eq:|P|EY-sum-1}), fix some
$\delta\in\mathcal{D}_{n}$ and some $T\in\mathcal{T}_{n}$ with $T\sim\delta$. We need to count the
number of pairings that include $T$, with the embedding of $T$ identified.
That is, if for some pairing $P$, the tree $T$ can be embedded in
$G\left(P\right)$ in multiple ways, then we count each different
way separately.

Now, exactly $\delta_{j}$ of the prevertices in bucket $j$ must contribute
to $T$, and there are $\falling{d}{\delta_{j}}$ ways to choose and
order these prevertices. So, there are
\[
\prod_{j=1}^{n}\falling{d}{\delta_{j}}
\]
ways to pair up the $n-1$ edges corresponding to a copy of $T$.
Then, there are 
\[
d n-2\left(n-1\right)=\left(d-2\right)n+2
\]
prevertices remaining, which can be paired in $\#P\left(\left(d-2\right)n+2\right)$
ways. This yields
\begin{equation}
\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y=
  \#P\left(\left(d-2\right)n+2\right)\,
 \sum_{\delta\in\mathcal{D}_{n}}\prod_{j=1}^{n}\, \falling{d}{\delta_{j}} \,
 \sum_{T\in\mathcal{T}_{n}}
   \, \boldsymbol{1}\left(T\sim\delta\right).
\label{eq:|P|EY-sum-2}
\end{equation}
The inner sum in (\ref{eq:|P|EY-sum-2}) is the number of trees with
degree sequence $\delta$, which is the multinomial
\begin{equation}
\binom{n-2}{\delta_{1}-1,\ldots,\delta_{n}-1}
  =\frac{\left(n-2\right)!}{\prod_{j=1}^{n}\left(\delta_{j}-1\right)!}.
\label{eq:moon}
\end{equation}
(See for example Moon~\cite[Theorem 3.1]{Moon}.)

Hence
\[
\left|\mathcal{P}_{n,\k}\right|\:\mathbb{E}Y = \left(n-2\right)!\,\#P\left(nd-2\left(n-1\right)\right)\sum_{\delta\in\mathcal{D}_{n}}\prod_{j=1}^{n}\frac{\falling{d}{\delta_{j}}}{\left(\delta_{j}-1\right)!}.
\]
It follows that the total number of ways to choose a spanning tree on $n$
vertices and choose a partial pairing that projects to that tree is 
\begin{align}
(n-2)!\, \sum_{\delta\in\mathcal{D}_{n}}\, \prod_{j=1}^{n}\, 
  \frac{\falling{d}{\delta_{j}}}{\left(\delta_{j}-1\right)!} 
  &=  (n-2)!\, \left[x^{2\left(n-1\right)}\right]\,
  \left(\sum_{j=1}^{\infty}\, 
   \frac{\falling{d}{j}}{\left(j-1\right)!}\, x^{j}\right)^{n}\nonumber \\
 &=  (n-2)!\, \left[x^{2\left(n-1\right)}\right]\, 
   \left(d x\left(1+x\right)^{d-1}\right)^{n}\nonumber\\
 &=  (n-2)!\, d^{n}\, \binom{(d-1)n}{n-2}. \label{useful}
\end{align}
Hence, by Stirling's approximation and (\ref{eq:num-pairings-asymptotic})
we conclude that
\begin{align}
\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E} Y & =  \left(n-2\right)!\,\#P\left(nd-2\left(n-1\right)\right)\,d^{n}\, \binom{(d-1)n}{n-2}\nonumber \\
 & \sim  \frac{\sqrt{2}\left(d-1\right)^{1/2}}{n\left(d-2\right)^{3/2}}\left(d\left(d-2\right)\left(d-1\right)^{d-1}\left(\frac{n}{\left(d-2\right)e}\right)^{d/2}\right)^{n}.\label{eq:|P|EY-asymptotic}
\end{align}

It follows that
\begin{align}
\mathbb{E} Y & = \frac{\left(n-2\right)!\,\#P\left(nd-2\left(n-1\right)\right)}{\#P\left(nd\right)}\,d^{n}\, \binom{(d-1)n}{n-2} \nonumber \\
 & \sim \frac{\left(d-1\right)^{1/2}}{n\left(d-2\right)^{3/2}}\left(\frac{\left(d-1\right)^{d-1}}{\left(d^{2}-2d\right)^{d/2-1}}\right)^{n}.
\label{eq:EY-asymptotic}
\end{align}
Hence for $d\ge3$ and $n$ sufficiently large, we have $\mathbb{E}Y\ne0$.


\section{Interaction with short cycles}\label{sec:E_YX}

\global\long\defX_{\r}{X_{\delta^{1}}}
\global\long\def\mathcal{P}_{n,\r}{\mathcal{P}_{n,\delta^{1}}}
\global\long\defs{s}
\global\long\defn_{\r}{n_{\delta^{1}}}
\global\long\def\left|\r\right|{\left|\delta^{1}\right|}
\global\long\def\delta'{\delta'}
\global\long\def\mathcal{D}_{Q}{\mathcal{D}_{Q}}
\global\long\deft{t}
\global\long\defE{E}
\global\long\def\mathcal{D}_{\dI}{\mathcal{D}_{\delta'}}
\global\long\def\mu{\mu}
\global\long\def\kappa{\kappa}
\global\long\def\Lambda{\Lambda}
\global\long\defg{g}
\global\long\def\mathcal{R}_{n,\r}{\mathcal{R}_{n,\delta^{1}}}
\global\long\def\mathcal{I}_{\i}{\mathcal{I}_{j}}
\global\long\def\mathcal{Q}_{\r}{\mathcal{Q}_{\delta^{1}}}

Recall that $X_{j}=\gamma_{j}\circ G$ is the number of cycles of
length $j$ in the projection of a random pairing $P\in\mathcal{P}_{n,\k}$.
For some fixed $m\ge0$, $\delta^{1}\in\mathbb{N}^{m}$, let $X_{\r}=\prod_{j=1}^{m}\falling{X_{j}}{\delta^{1}_{j}}$.
In this section we will compute an asymptotic formula for $\mathbb{E}\left[YX_{\r}\right]/\mathbb{E}Y$,
in the form required by Condition \ref{cond:A2'}.

We have
\[
\mathbb{E}\left[YX_{\r}\right]=\sum_{P\in\mathcal{P}_{n,\k}}\frac{1}{\left|\mathcal{P}_{n,\k}\right|}\,\tau\left(G\left(P\right)\right)\prod_{j=1}^{m}\falling{\gamma_{j}\left(G\left(P\right)\right)}{\delta^{1}_{j}}.
\]
Note that
\[
\prod_{j=1}^{m}\falling{\gamma_{j}\left(G\left(P\right)\right)}{\delta^{1}_{j}}
\]
is the number of ways to choose, for each $j\in\{ 1,\ldots, m\}$,
 an ordered set of $\delta^{1}_{j}$ cycles
of length $j$.
This will result in an ordered set of 
\[
\left|\r\right|=\sum_{j=1}^{m}\delta^{1}_{j}
\]
cycles.

We make the decomposition 
\[
\prod_{j=1}^{m}\falling{\gamma_{j}\left(G\left(P\right)\right)}{\delta^{1}_{j}}
  =\gamma_{\delta^{1}}^{\left(0\right)}+\gamma_{\delta^{1}}',
\]
where $\gamma_{\delta^{1}}^{\left(0\right)}$ is the number of ordered sets of cycles
in which each cycle is disjoint, and $\gamma_{\delta^{1}}'$ is the number
of ordered sets in which some vertices are shared between multiple cycles.
We can further decompose $\gamma_{\delta^{1}}'$ by the structure of the interaction
between the cycles. That is, according to the multigraph that is the
union of the cycles, and the specification of which edges of this
union belong to which cycle. This expresses $\gamma_{\delta^{1}}'$ as a sum
of terms $\gamma_{\delta^{1}}^{\left(j\right)}$. The number of terms $J$
in this decomposition depends on $\rho$, but is $O\left(1\right)$
as $n\to\infty$.

Define 
\[
E^{\left(j\right)}=\sum_{P\in\mathcal{P}_{n,\k}}\,
  \tau\left(G\left(P\right)\right)\, \gamma_{\delta^{1}}^{\left(j\right)},
\]
so that we have $\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}\left[YX_{\r}\right]=\sum_{j=1}^{J}E^{\left(j\right)}$.

We proceed to calculate $E^{\left(0\right)}$. As is standard when
applying this method (see for example, \cite[Theorem 9.6]{JLR}),
$\mathbb{E}\left[YX_{\r}\right]$ is asymptotically dominated by $E^{\left(0\right)}$
(the contribution due to disjoint cycles). See Lemma \ref{lem:Ej}
for some justification for this fact.

Let $\mathcal{C}_{n,j}$ be the set of all $j$-cycles on the vertex
set $\left\{ 1,\dots,n\right\} $, and define the Cartesian product
\[
\mathcal{C}_{n,\delta^{1}}=\prod_{j=1}^{m}\mathcal{C}_{n,j}^{\delta^{1}_{j}}.
\]
Each $R\in\mathcal{C}_{n,\delta^{1}}$ is an ordered set of $\left|\delta^{1}\right|$
cycles. We use the notation $R_{j,k}$ for the $k$th cycle of
length $j$ in $R$. Next, define
\[
\mathcal{R}_{n,\r} = \left\{ R\in\mathcal{C}_{n,\delta^{1}}:\mbox{ the cycles in $R$ are pairwise
 disjoint} \right\} .
\]
Similarly to (\ref{eq:|P|EY-sum-0}), we have
\[
E^{\left(0\right)} = \sum_{R\in\mathcal{R}_{n,\r}}\sum_{T\in\mathcal{T}_{n}}\sum_{P\in\mathcal{P}_{n,\k}}M_{\left(T,R\right),P},
\]
where $M_{\left(T,R\right),P}$ is the number of ways to embed the
tree $T$ and the cycles in $R$ into the multigraph $G\left(P\right)$.

We further condition on the edge intersection between the embedding
of $T$ and the cycles in $R$. We use a binary sequence of length
$j$ to encode the intersection of an $j$-cycle with a spanning
tree. Picking an arbitrary start vertex and direction for the cycle,
if the $k$th edge of the cycle is to be included in the intersection,
then the $k$th element of the corresponding sequence is one; otherwise
it is zero. All sequences $q\in\left\{ 0,1\right\} ^{j}$ represent
possible intersections, except the sequence $\left(1,\dots,1\right)$,
because a tree contains no cycles.

Define the set of all possible intersection sequences for a cycle
of length $j$, by 
\[
\mathcal{I}_{\i}=\left\{ 0,1\right\} ^{j}\setminus \left(1,\dots,1\right).
\]
Also, define the Cartesian product
\[
\mathcal{Q}_{\r}=\prod_{j=1}^{m}\mathcal{I}_{\i}^{\delta^{1}_{j}}.
\]
So, for each $R\in \mathcal{R}_{n,\r}$, specifying some $Q\in\mathcal{Q}_{\r}$ fully specifies the
intersection between the cycles in $R$ and a tree $T$.

We have 
\begin{equation}
E^{\left(0\right)}=\sum_{Q\in\mathcal{Q}_{\r}}\sum_{R\in\mathcal{R}_{n,\r}}\sum_{T\in\mathcal{T}_{n}}\sum_{P\in\mathcal{P}_{n,\k}}M_{\left(T,R,Q\right),P},\label{eq:|P|E[YX]}
\end{equation}
where $M_{\left(T,R,Q\right),P}$ is the number of ways to embed $T$
and the cycles in $R$ in $P$, such that the intersection between
the embedding of $T$ and the cycle $R_{j,k}$ is consistent with
$Q_{j,k}$, for $j = 1,\ldots, m$ and $k = 1,\ldots, \delta^{1}_{j}$.

Fixing $Q\in\mathcal{Q}_{\r}$, we will now evaluate the innermost triple sum in
(\ref{eq:|P|E[YX]}). Consider the following process:

\renewcommand{\theenumi}{\arabic{enumi}}
\renewcommand{\labelenumi}{\theenumi.}
\setenumerate{itemsep=5pt}
\begin{enumerate}
\item \label{step:3}Choose some $R\in\mathcal{R}_{n,\r}$.
\item \label{step:4} Choose a partial pairing that projects to $R$.
\item \label{step:6} Extend this to a pairing of a spanning tree consistent
with $Q$.
\item \label{step:7}Pair the remaining prevertices arbitrarily.
\end{enumerate}

We will find that the number of ways to complete each step is independent
of the other steps. Then, $E^{\left(0\right)}$ is a product of the
number of ways to complete each step, summed over all $Q\in\mathcal{Q}_{\r}$.

First let $n_{\r}=\sum_{j=1}^{m}j\delta^{1}_{j}$ be the total number of vertices
in each $R\in\mathcal{R}_{n,\r}$. The number of ways to choose the vertices for some
$R\in\mathcal{R}_{n,\r}$ is
\[
\binom{n}{n_{\r}}
\]
and the number of different arrangements of disjoint cycles on those
vertices is 
\[
\frac{n_{\r}!}{\prod_{j=1}^{m}\left(j\chi\left(j\right)\right)^{\delta^{1}_{j}}},
\]
where $j\chi\left(j\right)$ is the size of the automorphism group
of a $j$-cycle:
\[
\chi\left(j\right)=\begin{cases}
1 & \mbox{for }j\le2\\
2 & \mbox{for }j>2.
\end{cases}
\]
That is, the number of ways to complete Step \ref{step:3} is 
\[
s_{\ref{step:3}}^{\left(0\right)}=\frac{n!}{\left(n-n_{\r}\right)!\prod_{j=1}^{m}\left(j\chi\left(j\right)\right)^{\delta^{1}_{j}}}.
\]
Next, the number of ways to complete Step \ref{step:4} is
\[
s_{\ref{step:4}}^{\left(0\right)}=\prod_{j=1}^{m}\left(\frac{\chi\left(j\right)\left(d\left(d-1\right)\right)^{j}}{2}\right)^{\delta^{1}_{j}}.
\]
Note for future reference that we have 
\begin{equation}
s_{\ref{step:3}}^{\left(0\right)}s_{\ref{step:4}}^{\left(0\right)}\sim n^{n_{\r}}\prod_{j=1}^{m}\left(\frac{\left(d\left(d-1\right)\right)^{j}}{2j}\right)^{\delta^{1}_{j}}.\label{eq:f_3f_4}
\end{equation}
Next, we count the number of ways to extend this pairing to
a tree $T$ consistent with $Q$. We do this by constructing a new
irregular pairing model $\mathcal{P}_{n,\r}$ from the prevertices still unpaired after
Step \ref{step:4}. Recall that $Q$ describes a union of disjoint
paths; for each of these paths, combine the unpaired prevertices remaining
in each constituent vertex of the path to form a \emph{super-bucket}. If the
path has $k$ vertices then the resulting super-bucket
has $k\left(d-2\right)$ prevertices. Let $\left|Q\right|$ be the number
of super-buckets formed in this way, so the total number of buckets in
$\mathcal{P}_{n,\r}$ is $n':=n-n_{\r}+\left|Q\right|$.

Now, consider an extension of a pairing of cycles from Step \ref{step:4},
to a (partial) pairing of the edges of a tree $T$ consistent with
$Q$, as per Step \ref{step:6}. The pairs from this extension correspond
uniquely to a (partial) pairing $P'$ in the pairing model $\mathcal{P}_{n,\r}$.
By the construction of $\mathcal{P}_{n,\r}$, the projection $T'=G\left(P'\right)$
of this pairing is simply $T$ with some subpaths contracted
to single vertices. Since contracting edges of a tree cannot create
cycles, $G\left(P'\right)$ is itself a (spanning) tree. Similarly,
every pairing of a tree in $\mathcal{P}_{n,\r}$ corresponds to an extension of
a pairing of cycles to a pairing of a tree in $\mathcal{P}_{n,\k}$ consistent with
$Q$. So the number of ways $s_{\ref{step:6}}^{\left(0\right)}$
to complete Step \ref{step:6} equals the number of ways to
choose and pair up a spanning tree in $\mathcal{P}_{n,\r}$.

We will perform this count as in Section \ref{sec:EY}, by conditioning
on the degree in $T'$ of each bucket in $\mathcal{P}_{n,\r}$. Put an arbitrary
ordering on the $\left|Q\right|$ super-buckets, and let $d_{j}$
be the number of prevertices in the $j$th super-bucket. For a degree
sequence $\delta$, let $\left|\delta\right|$ be its degree sum. Define the
sets
\begin{align*}
\mathcal{D}_{Q} & = \left\{ \delta'\in\mathbb{N}^{\left|Q\right|}:\quad\delta'_{j}\le d_{j}\mbox{ for all }j\right\} ,\\
\mathcal{D}_{\dI} & = \left\{ \delta\in\mathbb{N}^{n-n_{\r}}:\quad\sum_{j=1}^{n-n_{\r}}\delta_{j}=2\left(n'-1\right)-\left|\delta'\right|\right\} .
\end{align*}
The set $\mathcal{D}_{Q}$ contains all possible degree-in-$T'$ sequences for
the $\left|Q\right|$ super-buckets. For some $\delta'\in\mathcal{D}_{Q}$, the set
$\mathcal{D}_{\dI}$ contains all possible degree sequences for the $n-n_{\r}$ remaining
ordinary buckets. So, we have
\[
s_{\ref{step:6}}^{\left(0\right)}=
  \sum_{\delta'\in\mathcal{D}_{Q}} \sum_{\delta\in\mathcal{D}_{\dI}}\, \sum_{T'\in\mathcal{T}_{n'}}\,
  \boldsymbol{1}\left(T'\sim\left(\delta',\delta\right)\right)\, \sum_{P'\in\mathcal{P}_{n,\r}}M_{T',P},
\]
where $T'\sim\left(\delta',\delta\right)$ denotes the event that the super-buckets
have degree-in-$T'$ sequence $\delta'$ and the remaining vertices have
degree-in-$T'$ sequence $\delta$. Proceeding as before, after fixing
some $\left(T',\delta',\delta\right)$, there are
\[
\left(\prod_{j=1}^{n-n_{\r}}\falling d{\delta_{j}}\right)\left(\prod_{j=1}^{\left|Q\right|}\falling{d_{j}}{\delta'_{j}}\right)
\]
ways to pair the edges of $T'$. There are 
\[
\frac{\left(n'-2\right)!}{\left(\prod_{j=1}^{n-n_{\r}}\left(\delta_{j}-1\right)!\right)\left(\prod_{j=1}^{\left|Q\right|}\left(\delta'_{j}-1\right)!\right)}
\]
trees with $T'\sim\left(\delta',\delta\right)$, by (\ref{eq:moon}). So,
we have
\begin{equation}
s_{\ref{step:6}}^{\left(0\right)}=\sum_{\delta'\in\mathcal{D}_{Q}}A_{\delta'}\prod_{j=1}^{\left|Q\right|}\frac{\falling{d_{j}}{\delta'_{j}}}{\left(\delta'_{j}-1\right)!},\label{eq:f_6}
\end{equation}
where
\begin{align}
A_{\delta'} & = \left(n'-2\right)!\sum_{\delta\in\mathcal{D}_{\dI}}\prod_{j=1}^{n-n_{\r}}\frac{\falling d{\delta_{j}}}{\left(\delta_{j}-1\right)!}\nonumber \\
 & = \left(n'-2\right)!\,\left[x^{2\left(n'-1\right)-\left|\delta'\right|}\right]\left(d x\left(1+x\right)^{d-1}\right)^{n-n_{\r}}\nonumber \\
 & = \left(n'-2\right)!\,d^{n-n_{\r}}\,
  \binom{(d-1)(n-n_{\r})}{2(n'-1)-|\delta'|-(n-n_{\r})}\nonumber \\
 & \sim \frac{\left(d-2\right)^{2\left|Q\right|-\left|\delta'\right|-5/2}\,\left(d-1\right)^{1/2}}{n^{n_{\r}-\left|Q\right|+2}}\left(\frac{\left(d-2\right)^{d-2}}{d\left(d-1\right)^{d-1}}\right)^{n_{\r}}\left(\frac{d\left(d-1\right)^{d-1}n}{e\left(d-2\right)^{d-2}}\right)^{n}.\label{eq:A_d'-asymptotic}
\end{align}
Finally, for Step \ref{step:7} there are $d n-2n_{\r}-2\left(n'-1\right)=\left(d-2\right)n-2\left(\left|Q\right|-1\right)$
prevertices remaining, which can be paired in 
\begin{equation}
s_{\ref{step:7}}^{\left(0\right)}=\#P\left(\left(d-2\right)n-2\left(\left|Q\right|-1\right)\right)\sim\sqrt{2}\left(\frac{\left(d-2\right)n}{e}\right)^{\left(d/2-1\right)n}\left(\left(d-2\right)n\right)^{1-\left|Q\right|}\label{eq:f_7}
\end{equation}
 ways.

Combining (\ref{eq:|P|EY-asymptotic}), (\ref{eq:f_3f_4}), (\ref{eq:f_6}),
(\ref{eq:A_d'-asymptotic}) and (\ref{eq:f_7}), we have
\begin{align}
\frac{E^{\left(0\right)}}{\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y} 
 & = \sum_{Q\in\mathcal{Q}_{\r}}\, 
  \frac{s_{\ref{step:3}}^{\left(0\right)}s_{\ref{step:4}}^{\left(0\right)}s_{\ref{step:6}}^{\left(0\right)}s_{\ref{step:7}}^{\left(0\right)}}
 {\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y}\nonumber \\
 & \to \left(\prod_{j=1}^{m}\left(\frac{\left(d\left(d-1\right)\right)^{j}}{2j}\right)^{\delta^{1}_{j}}\right)\nonumber \\
 &  \quad\times \sum_{Q\in \mathcal{Q}_{\r}}\, \sum_{\delta'\in\mathcal{D}_{Q}}\left(d-2\right)^{\left|Q\right|-\left|\delta'\right|}\left(\frac{\left(d-2\right)^{d-2}}{d\left(d-1\right)^{d-1}}\right)^{n_{\r}}\prod_{j=1}^{\left|Q\right|}\frac{\falling{d_{j}}{\delta'_{j}}}{\left(\delta'_{j}-1\right)!}.\label{eq:E0}
\end{align}

As is standard in these arguments, the only significant contribution to
$\mathbb{E}\left[Y X_{\r}\right]$ comes from $E^{\left(0\right)}$, where the cycles
do not overlap.  For completeness we sketch a proof of this below.

\begin{lemma}
\label{lem:Ej}$\mathbb{E}\left[YX_{\r}\right]$ is dominated by the contribution
from $E^{\left(0\right)}$. That is, 
\[
\frac{\mathbb{E}\left[YX_{\r}\right]}{\mathbb{E}Y}\sim\frac{E^{\left(0\right)}}{\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y}.
\]
\end{lemma}

\begin{proof}
We can estimate general $E^{\left(j\right)}$ with some slight modifications
to the above calculations. We would need a different $\mathcal{R}_{n,\r}'\subseteq\mathcal{C}_{n,\delta^{1}}$
that contains all possible ways to embed an ordered set of cycles
with a particular union $U$ into the vertex set $\left\{ 1,\dots,n\right\} $.
The intersection between a spanning tree and the cycles in some $R\in\mathcal{R}_{n,\r}$
would then be a subforest of $U$, so it would be more complicated
to explicitly define a set $\mathcal{Q}_{\r}$ that encodes all possibilities for
the intersection. However, the number of possible intersections is
still independent of $n$.

Using the same 4 steps, the decomposition $E^{\left(j\right)}=s_{\ref{step:3}}^{\left(j\right)}s_{\ref{step:4}}^{\left(j\right)}s_{\ref{step:6}}^{\left(j\right)}s_{\ref{step:7}}^{\left(j\right)}$
is still valid. Let $\left|U\right|$ and $\left\Vert U\right\Vert $
be the number of vertices and edges in the multigraph $U$, respectively.
Carefully adjusting the calculations for $E^{\left(0\right)}$, we
have $s_{\ref{step:3}}^{\left(j\right)}/s_{\ref{step:3}}^{\left(0\right)}=O\left(n^{\left|U\right|-n_{\r}}\right)$,
$s_{\ref{step:4}}^{\left(j\right)}/s_{\ref{step:4}}^{\left(0\right)}=O\left(1\right)$
and $s_{\ref{step:6}}^{\left(j\right)}/s_{\ref{step:6}}^{\left(0\right)}=O\left(1\right)$.
For Step \ref{step:7} there would be $d n-2\left\Vert U\right\Vert -2\left(n'-1\right)$
prevertices remaining, so $s_{\ref{step:7}}^{\left(j\right)}/s_{\ref{step:7}}^{\left(0\right)}=O\left(n^{n_{\r}-\left\Vert U\right\Vert }\right)$.

We conclude that $E^{\left(j\right)}/E^{\left(0\right)}=O\left(n^{\left|U\right|-\left\Vert U\right\Vert }\right)$.
Any non-disjoint union between distinct cycles has more edges than
vertices, so the lemma is proved.
\end{proof}
We now want to to express (\ref{eq:E0}) in the form required by \ref{cond:A2'}.
So, we consider each cycle independently. For a sequence $q\in\mathcal{I}_{\i}$,
let $q\left[k\right]$ be the number of paths with $k$ vertices
in the intersection encoded by $q$. So, for $Q\in\mathcal{Q}_{\r}$ we have 
\[
\left|Q\right|=\sum_{j=1}^{m}\sum_{\ell=1}^{\delta^{1}_{j}}\sum_{k=1}^{j}Q_{j,\ell}\left[k\right].
\]
Also, let $\left|q\right|=\sum_{k=1}^{j}q\left[k\right]$ be
the total number of paths in the intersection encoded by $q$.

Now, recall that if the $j$th super-bucket was collapsed from a
path of length $k$, then $d_{j}=k\left(d-2\right)$. Also recall
that $n_{\r}=\sum_{j=1}^{m}j\delta^{1}_{j}$ and note that 
\[
\left(d-2\right)^{\left|Q\right|-\left|\delta'\right|}=\prod_{j=1}^{\left|Q\right|}\left(d-2\right)^{1-\delta'_{j}}.
\]
As a result, we have
\[
\frac{\mathbb{E}\left[YX_{\r}\right]}{\mathbb{E}Y}\rightarrow \prod_{j=1}^{m}\left(\lambda_{j}'\right)^{\delta^{1}_{j}},
\]
where
\[
\lambda_{j}' = \frac{\left(d\left(d-1\right)\right)^{j}}{2j}\left(\frac{\left(d-2\right)^{d-2}}{d\left(d-1\right)^{d-1}}\right)^{j}\sum_{q\in\mathcal{I}_{\i}}\prod_{k=1}^{j}\left(\sum_{\ell=1}^{k\left(d-2\right)}\frac{\falling{k\left(d-2\right)}{\ell}}{\left(\ell-1\right)!\,\left(d-2\right)^{\ell-1}}\right)^{q\left[k\right]}.
\]
Note that $\ell$ takes the role of $\delta'_{j}$ for the $j$th super-bucket.
We have proved that Condition \ref{cond:A2'} is satisfied. It remains
to simplify our expression for $\lambda_{j}'$. We have
\begin{align*}
\lambda_{j}' & = \frac{1}{2j}\,
  \left(\frac{d-2}{d-1}\right)^{j\left(d-2\right)}\,
  \sum_{q\in\mathcal{I}_{\i}}\, \prod_{k=1}^{j}\, \left(k\left(d-2\right)\,
  \sum_{\ell=0}^{k\left(d-2\right)-1}\,
  \binom{k(d-2)-1}{\ell}\, 
  \left(d-2\right)^{-\ell}\right)^{q\left[k\right]}\\
 & = \frac{1}{2j}\left(\frac{d-2}{d-1}\right)^{j\left(d-2\right)}\sum_{q\in\mathcal{I}_{\i}}\prod_{k=1}^{j}\left(k\left(d-2\right)\left(\frac{1}{d-2}+1\right)^{k\left(d-2\right)-1}\right)^{q\left[k\right]}\\
 & = \frac{1}{2j}\sum_{q\in\mathcal{I}_{\i}}\left(\frac{\left(d-2\right)^{2}}{d-1}\right)^{\left|q\right|}\prod_{k=1}^{j}k^{q\left[k\right]}.
\end{align*}
Now, recall that $\left(1,\dots,1\right)\notin\mathcal{I}_{\i}$. So, to evaluate
the sum over $q$, we may identify a particular element in the sequence
to be zero. By symmetry, we arbitrarily choose the last. We also condition
on $\left|q\right|$: define
\begin{equation}
\mu  =  \frac{\left(d-2\right)^{2}}{d-1},\qquad
\Lambda_{j,t}  =  \sum_{\substack{q\in\mathcal{I}_{\i}\\
\left|q\right|=t\\
q_{j}=0
}
}\mu^{t}\prod_{k=1}^{j}k^{q\left[k\right]}.
\end{equation}
Note that $\Lambda_{j,1}=j\mu$ for all $j$, because the only sequence
$q\in\mathcal{I}_{\i}$ with $\left|q\right|=1$ and $q_{j}=0$ is $\left(1,\dots,1,0\right)$.
For $\left|q\right|>1$, the first path in (the intersection encoded
by) a sequence can contain anywhere between 1 and $j-1$ vertices.
Ranging over the possibilities, we have
\[
\Lambda_{j,t}=\sum_{k=1}^{j-1}k\mu\Lambda_{j-k,t-1}
\]
for $t>0$. To solve this recurrence, define the generating function
\[
\Lambda\left(x,y\right) = \sum_{j=1}^{\infty}\sum_{t=1}^{\infty}\Lambda_{j,t}x^{j}y^{t}.
\]
We have
\[
\Lambda\left(x,y\right)-\sum_{j=1}^{\infty}\Lambda_{j,1}x^{j}y=\sum_{j=1}^{\infty}\sum_{t=2}^{\infty}\sum_{k=1}^{j-1}k\mu\Lambda_{j-k,t-1}x^{j}y^{t},
\]
so
\begin{align*}
\Lambda\left(x,y\right) & = \sum_{k=1}^{\infty}k x^{k}y\mu\sum_{j=k+1}^{\infty}\sum_{t'=1}^{\infty}\Lambda_{j-k,t'}x^{j-k}y^{t'}+\sum_{j=1}^{\infty}j x^{j}y\mu\\
 & = \left(\Lambda\left(x,y\right)+1\right)\sum_{k=1}^{\infty}k x^{k}y\mu.
\end{align*}
Now, defining
\[
g\left(x\right)=\sum_{k=1}^{\infty}k x^{k}\mu=\frac{x\mu}{\left(1-x\right)^{2}},
\]
we have 
\[
\Lambda\left(x,y\right)=\frac{g\left(x\right)y}{1-g\left(x\right)y}.
\]
If $q\in\mathcal{I}_{\i}$, then there are $j$ positions to place a zero, and
if $\left|q\right|=t$ then there are $t$ zeros in $q$.  Hence
\begin{align*}
\lambda_{j}' & =  \frac{1}{2j}\sum_{t=1}^{j}\frac{j\Lambda_{j,t}}{t}\\
  &=  \frac{1}{2}\left[x^{j}\right]\sum_{t=1}^{j}\frac{1}{t}\left[y^{t-1}\right]\frac{g\left(x\right)}{1-g\left(x\right)y}\\
 & =  \frac{1}{2}\left[x^{j}\right]\int_{0}^{1}\frac{g\left(x\right)}{1-g\left(x\right)y}\mathrm{d}y\\
 & =  -\frac{1}{2}\left[x^{j}\right]\log\left(1-g\left(x\right)\right).
\end{align*}
Now, defining $\kappa=\sqrt{d-1}$ we
have
\[
1-g\left(x\right)=\frac{1+x^{2}-\left(2+\mu\right)x}{\left(1-x\right)^{2}}=\frac{\left(1-\kappa^{2}x\right)\left(1-\kappa^{-2}x\right)}{\left(1-x\right)^{2}},
\]
so 
\begin{align}
\lambda_{j}' & = \frac{1}{2}\left[x^{j}\right]\left(2\log\left(1-x\right)-\log\left(1-\kappa^{2}x\right)-\log\left(1-\kappa^{-2}x\right)\right) \nonumber\\
 & = \frac{1}{2}\left[x^{j}\right]\sum_{k=1}^{\infty}\frac{-2x^{k}+\left(\kappa^{2}x\right)^{k}+\left(\kappa^{-2}x\right)^{k}}{k}\nonumber\\
 & = \frac{1}{2j}\left(\kappa^{j}-\kappa^{-j}\right)^{2}\nonumber\\
 & = \frac{\left(\left(d-1\right)^{j}-1\right)^{2}}{2j\left(d-1\right)^{j}}
\label{eq:lambda'}
\end{align}
for $j \geq 1$.

To complete this section we will 
establish that conditions \ref{cond:A2} and \ref{cond:A3} of 
Theorem~\ref{thm:janson} hold, and prove
Theorem~\ref{thm:expectation}.

\begin{lemma}
Let $d\geq 3$ be a fixed integer.  Then 
Conditions~\ref{cond:A2} and~\ref{cond:A3} of Theorem~\ref{thm:janson} are 
satisfied, and 
\[
\exp\left(\sum_{j=1}^{\infty}\lambda_{j}\zeta_{j}^{2}\right)=\frac{d^{2}}{\sqrt{\left(d-1\right)\left(d-2\right)\left(d^{2}-d+1\right)}}.
\]
\label{EY2/EY2-target}
\end{lemma}

\begin{proof}
The calculations of this section show that Condition~\ref{cond:A2'}
of Lemma~\ref{lem:janson} is satisfied with $\lambda_{j}'$ given by
(\ref{eq:lambda'}).  Then Lemma~\ref{lem:janson} guarantees that \ref{cond:A2}
is satisfied.
Using the Taylor expansion of $\log\left(1-z\right)$,
it follows from (\ref{eq:janson-parameters}) that
\begin{align*}
& \sum_{j=1}^{\infty}\lambda_{j}\zeta_{j}^{2} \\
&=
\sum_{j=1}^{\infty} \, \dfrac{1}{2j}\left(4\left(d-1\right)^{-j}-4\left(d-1\right)^{-2j}+\left(d-1\right)^{-3j}\right)\\
&=\dfrac{1}{2}\left(-4\log\left(1-\left(d-1\right)^{-1}\right)+4\log\left(1-\left(d-1\right)^{-2}\right)-\log\left(1-\left(d-1\right)^{-3}\right)\right).
\end{align*}
Taking the exponential of both sides and rearranging establishes
the stated expression for 
$\exp\left(\sum_{j = 1}^\infty \lambda_{j}\zeta_{j}^2\right)$,
which is finite for $d\geq 3$.  Hence Condition~\ref{cond:A3} holds, 
as required.
\end{proof}


\begin{proof}[Proof of Theorem \ref{thm:expectation}]
We know that condition~\ref{cond:A2} holds (as proved above), and hence
\[
\mathbb{E} Y_{\mathcal{G}}  =  \mathbb{E}\left[Y|X_{1}=X_{2}=0\right]\\
  \rightarrow   \mathbb{E} Y\,\exp\left(-\lambda_{1}\zeta_{1}-\lambda_{2}\zeta_{2}\right).
\]
Substituting using (\ref{eq:janson-parameters}) and (\ref{eq:EY-asymptotic}) 
completes the proof.
\end{proof}

\section{The second moment}\label{sec:EY2}

\global\long\def\delta^{1}{\delta^{1}}
\global\long\def\delta^{2}{\delta^{2}}
\global\long\def\eta^{1}{\eta^{1}}
\global\long\def\eta^{2}{\eta^{2}}
\global\long\def\eta^{3}{\eta^{3}}
\global\long\def\floor#1{\left\lfloor #1\right\rfloor }
\global\long\def\ceil#1{\left\lceil #1\right\rceil }
\global\long\defC{C}

\global\long\def\mathcal{X}{\mathcal{X}}

\global\long\def\Xk#1{\mathcal{X}^{\left(#1\right)}}

\global\long\defF_n{F_n}

\global\long\def\Hk#1{F_n^{\left(#1\right)}}

\global\long\defx{x}

\global\long\defy{y}

\global\long\deff{f}

\global\long\def\mathcal{S}{\mathcal{S}}
\global\long\defa_{n}{a_{n}}
\global\long\def\hat{a}_{n}{\hat{a}_{n}}
\global\long\def\mathcal{O}{\mathcal{O}}

We now want to calculate $\mathbb{E}Y^{2}$. 
As a first step we transform this problem into one of evaluating
the coefficient of a certain generating function.

\begin{lemma}
Let $d\geq 3$ be fixed and define
\[ N(n,d) = \begin{cases} n & \text{ if $d\geq 4$,}\\
            n/2 + 2 & \text{ if $d=3$}.
\end{cases}\]
Then
\begin{align*}
& \left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y^{2} \\
&= \frac{n!\, \left((d-2)n\right)!\, d^n}{2^{(d/2-1)n+2}}\,
   \sum_{b=1}^{N(n,d)}\, \frac{2^b}
  {b!\,\left(\left(d/2-1\right)n-b+2\right)!}\,
  \left[z^{n}\right]\,\left(\sum_{j=1}^{\infty}
 \binom{(d-1)j}{j}\, z^{j}\right)^{b}.
\end{align*}
\label{lem:|P|EY2}
\end{lemma}

\begin{proof}
We write 
\[
\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y^{2} = \sum_{P\in\mathcal{P}_{n,\k}}\, \sum_{T_{1}\in\mathcal{T}_{n}}\,
  \sum_{T_{2}\in\mathcal{T}_{n}}\, M_{\left(T_{1},T_{2}\right),P},
\]
where $M_{(T_{1},T_{2}),P}$ is the number of ways to embed the ordered
pair of trees $\left(T_{1},T_{2}\right)$ into the multigraph $G\left(P\right)$.
We will estimate this sum by choosing some $T_{1},T_{2}\in\mathcal{T}_{n}$ and
counting the ways to pair up their edges, then counting the ways to
complete the pairing. We break up this process in a similar way to
Section \ref{sec:E_YX}:

\begin{enumerate}
\item Choose $b\in\{1,\ldots, n\}$, which will be the number of 
connected components in the 
intersection of the embeddings of $T_{1}$ and $T_{2}$.
(As we will see later, when $d=3$ we must restrict to $b\in \{1,\ldots,
n/2 + 2\}$.)
\item \label{step:EY2-1}
Choose a partition $(\nu_1,\ldots, \nu_b)$ of $n$ into positive parts.
That is, $\nu_{j}$ is a positive integer for $j=1,\ldots, b$ and
$\sum_{j=1}^b \nu_{j} = n$.  Here $\nu_{j}$ will be the number of
vertices in the $j$th connected component of the interesection.
(We should divide by
$b!$ to account for our assumption that the connected components
are labelled).
\item \label{step:EY2-2}
Choose a partition of the $n$ vertices into $b$ groups, where the
size of the $j$th group is $\nu_{j}$.
\item \label{step:EY2-3} In each group, choose a spanning tree on that
group and choose a partial pairing that projects to that tree.
This specifies a component of the intersection.
\end{enumerate}

Now, collapse the buckets in each group into a single super-bucket,
giving exactly $b$ super-buckets.  The
$j$th super-bucket has $d\nu_{j}-2\left(\nu_{j}-1\right)$
unpaired prevertices. We now want to pair up two pair-disjoint spanning
trees $T_{1}',T_{2}'$ in the collapsed pairing model. These will
extend to $T_{1}$ and $T_{2}$ using the intersection subtrees chosen
in Step \ref{step:EY2-3}.\vspace{-5pt}

\begin{enumerate}
\setcounter{enumi}{4}
\item 
\label{step:EY2-4}
For $j=1,\ldots, b$, choose $\delta^{1}_{j}$ and
$\delta^{2}_{j}$, the degree of vertex $j$ in $T_{1}'$ and $T_{2}'$
respectively, in such a way that 
$\delta^{1}_{j}+\delta^{2}_{j}\led\nu_{j}-2\left(\nu_{j}-1\right)$.
These must also satisfy
\[ \sum_{j=1}^b \delta^{1}_{j} = \sum_{j=1}^b \delta^{2}_{j} = 2(b-1),\]
as they are the degree sequence of a spanning tree on $b$ vertices.
\item \label{step:EY2-5}
Choose two trees $T_{1}',T_{2}'$ on the $b$ vertices
that are consistent with the degree sequences chosen in Step \ref{step:EY2-5}.
\item \label{step:EY2-6}
Pair up these two trees in a pair-disjoint way.
\item \label{step:EY2-7}
Pair all remaining prevertices to complete a $d$-regular pairing.
\end{enumerate}

Given $\nu$, the number of ways to complete Step \ref{step:EY2-2}
is 
\[
s_{\ref{step:EY2-2}}= \binom{n}{\nu_{1},\ldots,\nu_{b}}.
\]
By (\ref{useful}), 
the number of ways to complete Step \ref{step:EY2-3} is 
\[
s_{\ref{step:EY2-3}}=
d^{n}\, \prod_{j=1}^{b}\, \left(\nu_{j}-2\right)!\,
\binom{(d-1)\nu_{j}}{\nu_{j}-2}.
\]
The number of ways to complete Step \ref{step:EY2-5} is 
\[
s_{\ref{step:EY2-5}}= \binom{b-2}{\delta^{1}_{1}-1,\ldots,\delta^{1}_{b}-1}\,
  \binom{b-2}{\delta^{2}_{1}-1,\ldots,\delta^{2}_{b}-1},
\]
 by (\ref{eq:moon}), and the number of ways to complete Step \ref{step:EY2-6}
is 
\[
s_{\ref{step:EY2-6}}=\prod_{j=1}^{b}\falling{(d-2)\nu_{j}+2}{\delta^{1}_{j}+\delta^{2}_{j}}.
\]
Finally, for Step \ref{step:EY2-7} there are 
\[
\sum_{j=1}^{b}\left(d\nu_{j}-2\left(\nu_{j}-1\right)-\delta^{1}_{j}-\delta^{2}_{j}\right) 
  =  \left(d-2\right)n-2\left(b-2\right)
\]
prevertices remaining, so the number of ways to complete Step \ref{step:EY2-7}
is 
\[
s_{\ref{step:EY2-7}}=\#P\left(\left(d-2\right)n-2\left(b-2\right)\right).
\]
For this construction to make sense, the quantity $(d-2)n-2(b-2)$ must be 
nonnegative. This is certainly true when $d\geq 4$, but when $d=3$ this
imposes the constraint that $b\leq n/2+2$.  (This explains the definition
of $N(n,d)$ in the lemma statement.)

It will be convenient to work with the nonnegative variables
\begin{equation}
\label{translate} \eta^{1}_{j}=\delta^{1}_{j}-1, \quad \eta^{2}_{j}=\delta^{2}_{j}-1 \quad \text{ and }
 \quad \eta^{3}_{j}=\left(d-2\right)\nu_{j}-\eta^{1}_{j}-\eta^{2}_{j}
\end{equation}
defined for $j=1,\ldots, b$.
Let 
\[ \mathcal{S}_{\ref{step:EY2-1}}(b) = \{ \nu \in \{ 1,\ldots, n\}^b : 
       \sum_{j=1}^b \nu_{j} = n \}
\]
be the set of possible sequences $\nu$
from Step \ref{step:EY2-1} and let 
\begin{align*} \mathcal{S}_{\ref{step:EY2-4}}(\nu) 
 = \left\{ \left(\eta^{1},\eta^{2},\eta^{3}\right)\in\left(\mathbb{N}^{b}\right)^{3}: \right. &
   \left.
  \quad\eta^{1}_{j}+\eta^{2}_{j}+\eta^{3}_{j}=\left(d-2\right)\nu_{j} \, \mbox{ for }\,
  j=1,\ldots, b,\right.\\
 &  \left. {} \quad 
   \sum_{j=1}^b \eta^{1}_{j} = \sum_{j=1}^b \eta^{2}_{j} = b-2\right\} 
\end{align*}
be the set of sequences arising from Step~\ref{step:EY2-4} using (\ref{translate}).

Combining all of the above gives
\begin{align*}
\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y^{2} & = \sum_{b=1}^{N(n,d)}\, 
 \frac{1}{b!}\, \sum_{\nu\in\mathcal{S}_{\ref{step:EY2-1}}(b)}\, 
  s_{\ref{step:EY2-2}}s_{\ref{step:EY2-3}}\,
  \sum_{\left(\eta^{1},\eta^{2},\eta^{3}\right)\in\mathcal{S}_{\ref{step:EY2-4}}(\nu)}\,
  s_{\ref{step:EY2-5}}s_{\ref{step:EY2-6}}s_{\ref{step:EY2-7}}\\
 & = \sum_{b=1}^{N(n,d)}\, \frac{n!\, d^{n}}
  {b!\,\left(\left(d/2-1\right)n-b+2\right)!\,2^{\left(d/2-1\right)n-b+2}}\, 
\sum_{\nu\in\mathcal{S}_{\ref{step:EY2-1}}(b)}\,
 \prod_{j=1}^{b}\, \frac{\left(\left(d-1\right)\nu_{j}\right)!}{\nu_{j}!}\\
 &  \quad\times\sum_{\left(\eta^{1},\eta^{2},\eta^{3}\right)\in\mathcal{S}_{\ref{step:EY2-4}}(\nu)}
 \, \binom{b-2}{\eta^{1}_{1},\ldots,\eta^{1}_{b}}\, \binom{b-2}{\eta^{2}_{1},\ldots,\eta^{2}_{b}}
  \,\binom{(d-2)n-2(b-2)}{\eta^{3}_{1},\dots,\eta^{3}_{b}}.
\end{align*}
Now
\begin{align*}
 &  \sum_{\left(\eta^{1},\eta^{2},\eta^{3}\right)\in\mathcal{S}_{\ref{step:EY2-4}}(\nu)}\,
  \binom{b-2}{\eta^{1}_{1},\ldots,\eta^{1}_{b}}\,
  \binom{b-2}{\eta^{2}_{1},\ldots,\eta^{2}_{b}}\,
  \binom{(d-2)n-2(b-2)}{\eta^{3}_{1},\ldots,\eta^{3}_{b}}\\
 & = \sum_{\left(\eta^{1},\eta^{2},\eta^{3}\right)\in\mathcal{S}_{\ref{step:EY2-4}}(\nu)}\,
  \left[z_{1}^{\eta^{1}_{1}}\dots z_{b}^{\eta^{1}_{b}}\right]\,
   \left(\sum_{j=1}^{b}z_{j}\right)^{b-2}\\
 &  \quad\times\left[z_{1}^{\eta^{2}_{1}}\dots z_{b}^{\eta^{2}_{b}}\right]\left(\sum_{j=1}^{b}z_{j}\right)^{b-2}\left[z_{1}^{\eta^{3}_{1}}\dots z_{b}^{\eta^{3}_{b}}\right]\left(\sum_{j=1}^{b}z_{j}\right)^{\left(d-2\right)n-2\left(b-2\right)}\\
 & = \left[ z_1^{(d-2)\nu_1}\, z_2^{(d-2)\nu_2} \cdots z_b^{(d-2)\nu_b}
  \right]\, \left(\sum_{j=1}^{b}z_{j}\right)^{\left(d-2\right)n}\\
 & = \binom{(d-2)n}{(d-2)\nu_{1},\ldots,(d-2)\nu_{b}}.
\end{align*}
It follows that
\begin{align*}
& \left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y^{2} \\
  &= \frac{n!\, \left((d-2)n\right)!\, d^n}{2^{(d/2-1)n+2}}\,
  \sum_{b=1}^{N(n,d)} \, \frac{2^b}{b!\,\left(\left(d/2-1\right)n-b+2\right)!}
 \, \sum_{\nu\in\mathcal{S}_{\ref{step:EY2-1}}(b)}\,
   \prod_{j=1}^{b}\, \binom{(d-1)\nu_{j}}{\nu_{j}}. 
\end{align*}
This is equal to the expression in the lemma statement,
by definition of $\mathcal{S}_{\ref{step:EY2-1}}(b)$. 
\end{proof}


We now seek to evaluate
\[ [z^n] \left( \sum_{j=1}^\infty \binom{(d-1)j}{j} z^{j}\right)^b.\]
By Stirling's approximation and the ratio test, the radius of
convergence of the series $\sum_{j=1}^{\infty}\, \binom{(d-1)j}{j}\, z^{j}$
equals $\frac{\left(d-2\right)^{d-2}}{\left(d-1\right)^{d-1}}$.
Hence, 
\[
f\left(z\right):=\sum_{j=1}^{\infty}\,
 \binom{(d-1)j}{j}\, \left(\frac{\left(d-2\right)^{d-2}}{\left(d-1\right)^{d-1}}\right)^{j}z^{j}
\]
 is analytic in the disk $\left\{ z:\left|z\right|<1\right\} $. Define
$\beta=b/n$ and let $r_{\beta}\in\left(0,1\right)$ be fixed
for each $\beta$ (we will determine this later). Then, with the contour
$\Gamma:\left[-\pi,\pi\right]\to\mathbb{C}$ defined by 
$\theta\mapsto r_{\beta}e^{i\theta}$, we have
\begin{align}
\left[z^{n}\right]f\left(z\right)^{b} & = 
  \frac{1}{2\pi i}\, \int_{\Gamma}\frac{f\left(z\right)^{b}}{z^{n+1}}\,
    \mathrm{d}z\nonumber \\
 & = \frac{1}{2\pi}\int_{-\pi}^{\pi}\left(\frac{f\left(r_{\beta}e^{i\theta}\right)^{\beta}}{r_{\beta}e^{i\theta}}\right)^{n}\mathrm{d}\theta,\label{eq:contour-integral}
\end{align}
by Cauchy's coefficient formula.
Let
\[ \mathcal{X}_{n}=\left\{ \dfrac{1}{n},\dfrac{2}{n},\dots,\dfrac{N(n,d)}{n}\right\} 
   \times\left[-\pi,\pi\right]
\]
be the sample space for pairs $\left(\beta,\theta\right)$, 
and define
$a_{n}:\mathcal{X}_{n}\to\mathbb{C}$ by
\begin{equation}
\label{andef}
a_{n}\left(\beta,\theta\right)=\frac{2^{\beta n}}{\left(\beta n\right)!\,\left(\left(d/2-1\right)n-\beta n+2\right)!}\left(\frac{f\left(r_{\beta}e^{i\theta}\right)^{\beta}}{r_{\beta}e^{i\theta}}\right)^{n}.
\end{equation}
Finally, let
\begin{equation}
\label{eq:E}
F_n=\sum_{b=1}^{N(n,d)}\int_{-\pi}^{\pi}a_{n}\left(b/n,\theta\right)\,\mathrm{d}\theta.
\end{equation}
Then, by Lemma~\ref{lem:|P|EY2} and (\ref{eq:contour-integral}), 
\begin{equation}
\left|\mathcal{P}_{n,\k}\right|\,\mathbb{E}Y^{2} = 
\frac
{\left(d-1\right)^{n\left(d-1\right)}n!\,\left(\left(d-2\right)n\right)!\,d^{n}}
{2\pi\left(d-2\right)^{n\left(d-2\right)}2^{\left(d/2-1\right)n+2}}\, F_n .
\label{eq:E2}
\end{equation}

We now apply the saddle point method to estimate the sum in
(\ref{eq:E}) in the case that $d=3$. 
Our proof is adapted from that of~\cite[Theorem 2.3]{GJR}.

When $d=3$ the function $f$ satisfies 
\[
f\left(z\right)=\sum_{j=1}^{\infty}\, \binom{2j}{j}\, 
  \left(\frac{z}{4}\right)^{j}=\left(1-z\right)^{-1/2}-1.
\]
We note for later that if $\theta\in [-\pi,\pi]$ is nonzero then
\begin{equation}
 \left| f\left( r_\beta e^{i\theta}\right)\right| = 
   \left| \, \sum_{j=1}^\infty 
     \binom{2j}{j}\, \left(\frac{r_\beta}4\right)^j
 e^{ij\theta}\, \right|
   <
\sum_{j=1}^\infty 
     \binom{2j}{j}\, \left(\frac{r_\beta}4\right)^j
 = |f(r_\beta)|,
\label{lem:technical}
\end{equation} 
using the triangle inequality.
Hence for each $\beta$ the function $\theta \mapsto |f( r_\beta e^{i\theta})|$
on $[-\pi,\pi]$ is uniquely maximised at $\theta =0$.


Define $\mathcal{X}=\left(0,1\right]\times\left[-\pi,\pi\right]$ 
and let $\mathcal{X}^{*}\subset\mathcal{X}$ be a set (to be determined) such that for $\left(\beta,\theta\right)\in\mathcal{X}^{*}$,
both $\beta$ and $d/2-1-\beta=1/2-\beta$ are bounded below by some
positive constant. 
Then Stirling's approximation gives, for $\beta\in\mathcal{X}^*\cap\mathcal{X}_n$,
\begin{align}
& a_{n}\left(\beta,\theta\right) \nonumber\\
 & \sim \frac{e^{n/2}}{2\pi\, n^{n/2+3}\, \sqrt{\beta\left(1/2-\beta\right)}\, \left(1/2-\beta\right)^{2}}\left(\frac{\left(2f\left(r_{\beta}e^{i\theta}\right)\right)^{\beta}}{r_{\beta}e^{i\theta}\beta^{\beta}\left(1/2-\beta\right)^{\left(1/2-\beta\right)}}\right)^{n}.\label{eq:g-asymptotic-working}
\end{align}
Next, define the half-spaces $\mathcal{X}^{1/2}=\left(0,1/2\right]\times\left[-\pi,\pi\right]$
and $\bar{\mathcal{X}}^{1/2}=\left[0,1/2\right]\times\left[-\pi,\pi\right]$.
Define the real-valued sequence $\left(c_{n}\right)_{n\in\mathbb{N}}$
and the functions $\psi:\bar{\mathcal{X}}^{1/2}\to\mathbb{R}$ and $\phi:\mathcal{X}^{1/2}\to\mathbb{C}$
by 
\begin{align*}
c_{n} & = \frac{e^{n/2}}{2\pi n^{n/2+3}},\\
\psi\left(\beta,\theta\right) & = \beta^{-1/2}\left(1/2-\beta\right)^{-5/2},\\
\phi\left(\beta,\theta\right) & = \beta\log\left(2f\left(r_{\beta}e^{i\theta}\right)\right)-\log r_{\beta}-i\theta-\beta\log\beta-\left(1/2-\beta\right)\log\left(1/2-\beta\right),
\end{align*}
so that we have 
\begin{equation}
a_{n}\left(\beta,\theta\right)\sim c_{n}\, \psi\left(\beta,\theta\right)\,
 e^{n\phi\left(\beta,\theta\right)}\label{eq:g-asymptotic-statement}
\end{equation}
uniformly for $\left(\beta,\theta\right)\in\mathcal{X}_{n}\cap\mathcal{X}^{*}$.

Let $D$ denote the differential operator 
\[ \left(D\phi\left(x\right)\right)_{j}=
   \frac{\partial\phi\left(x\right)}{\partial x_{j}}.\]
We seek a stationary point of $\phi$.  The condition
$\frac{\partial\phi\left(\beta,0\right)}{\partial\theta}=0$
is equivalent to the condition $\beta r_\beta f'(r_\beta) = f(r_\beta)$.
Solving for $r_\beta$ gives
\[
r_{\beta}=\frac{1}{8}\left(8-4\beta-\beta^{2}\pm\sqrt{\beta^{3}\left(8+\beta\right)}\right).
\]
We choose
$r_{\beta}=\frac{1}{8}\left(8-4\beta-\beta^{2}-
   \sqrt{\beta^{3}\left(8+\beta\right)}\right)\in\left(0,1\right)$,
which ensures that
$\frac{\partial\phi\left(\beta,0\right)}{\partial\theta}=0$.
Next, we calculate that with this choice of $r_\beta$,
\[ \frac{\partial\phi}{\partial \beta}(\beta,0) = 
  \log\left(\frac{\left(4-\beta-\sqrt{\beta(8+\beta)}\right)(1-2\beta)}
                 {\beta\left(\beta + \sqrt{\beta(8+\beta)}\right)}\right).
\]
Setting this equal to 0 and solving for $\beta$ gives the equation
$(3\beta-1)(\beta^2-4\beta+2)=0$.
The only solution with $\beta\in\left(0,\frac12\right]$ is $\beta=\frac13$ so we choose
$x^{*}=\left(\frac{1}{3},0\right)$ and check that 
$D\phi\left(x^{*}\right)=0$.

Note that $\phi\left(x^{*}\right)=\log\left(4\sqrt{\frac{2}{3}}\right)$,
and
\[
H=-\begin{pmatrix}\frac{63}{5} & 0\\
0 & \frac{5}{2}
\end{pmatrix}
\]
is the Hessian matrix of $\phi$ at $x^{*}$. Define $C_{1}=5/8$,
so that $-4C_{1}$ is the largest eigenvalue of $H$.

Now, define $\hat{\phi}$ by $\hat{\phi}\left(x\right)=\phi\left(x\right)-\phi\left(x^{*}\right)$,
and define $\hat{a}_{n}:\mathcal{X}_{n}\to\mathbb{C}$ by 
\[
\hat{a}_{n}\left(x\right)=
 c_n^{-1}\, e^{-n\phi\left(x^{*}\right)}\, 
  a_{n}\left(x\right).
\]
With a Taylor expansion about $x^{*}$, for $x\in\mathcal{X}^{1/2}$ we have
\begin{equation}
\hat{\phi}\left(x\right) = \frac{1}{2}\left(x-x^{*}\right)^{T}H\left(x-x^{*}\right)+h\left(x\right)\left|x-x^{*}\right|^{2},\label{eq:taylor}
\end{equation}
where $h\left(x\right)$ is complex and $h\left(x\right)\to0$ as
$x\tox^{*}$. For all $v\in\mathbb{R}^{2}$ we have $v^{T}Hv\le-2C_{1}\left|v\right|^{2}$,
so we can choose $\xi<\frac{1}{6}$ such that $\Re\hat{\phi}\left(x\right)\le-C_{1}\left|x-x^{*}\right|^{2}$
for $\left|x-x^{*}\right|<\xi$. Define $\mathcal{X}^{*}=\left\{ x\in\mathcal{X}^{1/2}:\left|x-x^{*}\right|<\xi\right\} $,
satisfying the requirement for (\ref{eq:g-asymptotic-working}).

Next, define the sets 
\begin{align*}
\Xk 1 & = \left\{ x\in\mathcal{X}^{*}:\left|x-x^{*}\right|<n^{-1/3}\right\} ,\\
\Xk 2 & = \mathcal{X}^{*}\setminus \Xk 1,\\
\Xk 3 & = \mathcal{X}^{1/2}\setminus \mathcal{X}^{*},\\
\Xk 4 & = \mathcal{X}\setminus \mathcal{X}^{1/2},
\end{align*}
so that with
\[
\Hk{j}=\sum_{b=1}^{n/2+2}\int_{-\pi}^{\pi}\hat{a}_{n}\left(b/n,\theta\right)\,
  \boldsymbol{1}_{\Xk{j}}(b/n,\theta)\, \mathrm{d}\theta
\]
we have 
\begin{equation}
\label{split}
 c_n^{-1}\, e^{-n\phi\left(x^{*}\right)} F_n=\Hk 1+\Hk 2+\Hk 3+\Hk 4.
\end{equation}

\begin{lemma}
With notation as above, we have
\[ \Hk 1 + \Hk 2 + \Hk 3 + \Hk 4 \sim \Hk 1
   \sim \frac{144\pi}{\sqrt{7}}.\]
\label{lem:dominates}
\end{lemma}

\begin{proof}
Note that $\psi$ is a continuous function defined on a compact set.
So, $\psi$ is absolutely bounded on its domain, by $C_{2}$ say.
By (\ref{eq:g-asymptotic-statement}), it follows that $\left|\hat{a}_{n}\left(x\right)\right|=O\left(e^{n\Re\hat{\phi}\left(x\right)}\right)$
uniformly for $x\in\mathcal{X}^{*}$. For $x\in\Xk 2$ we have 
\[
n\Re\hat{\phi}\left(x\right)\le-nC_{1}\left|x-x^{*}\right|^{2}\le-C_{1}n^{1/3}\to-\infty
\]
and consequently
\begin{equation}
\left|\Hk 2\right|  \, \leq \, 
  \sum_{b=1}^{n} \int_{-\pi}^{\pi}\, \left|\hat{a}_{n}\left(b/n,\theta\right)\right|\,
  \boldsymbol{1}_{\Xk 2}(b/n,\theta)\, \mathrm{d}\theta \,
  = \, O\left(ne^{-C_{1}n^{1/3}}\right) = o(1).
\label{F2bound}
\end{equation}

Now (\ref{lem:technical}) implies that
for each $\beta$, $\Re\hat{\phi}\left(\beta,\theta\right)$
is uniquely maximized when $\theta=0$. 
Also
$\frac{\partial\Re\hat{\phi}}{\partial\beta}\left(\beta,0\right)=\frac{\partial\hat{\phi}}{\partial\beta}\left(\beta,0\right)=0$ only for $\left(\beta,0\right)=x^*$, since $\hat{\phi}$ is real along the line $\theta=0$.
Checking the values of $\Re\hat{\phi}(\beta,0)$ in the limit as $\beta\to 0$ 
and $\beta\to \nfrac{1}{2}$, 
it follows that $\Re\hat{\phi}$ attains a unique
maximum on $\mathcal{X}^{1/2}$ at $x^{*}$. 
Let $-C_{3}<0$ be the maximum value of $\Re\hat{\phi}$ on $\Xk 3$.

Let $u\lor w=\max\left\{ u,w\right\} $ for real numbers $u,w$. 
We now redo the calculations of (\ref{eq:g-asymptotic-working}) using an
alternate form of Stirling's inequality which holds for all $k\geq 0$,
namely $\sqrt{k\lor1}\left(\frac{k}{e}\right)^{k}\le k!$.
For $\left(\beta,\theta\right) \in \mathcal{X}^{1/2}\cap \mathcal{X}_n$,
\begin{align*}
\left|\hat{a}_{n}\left(\beta,\theta\right)\right| & \le 
\frac{ne^{2}}{\sqrt{\left(\beta n\lor 1\right)
  \left(\left(n/2-\beta n+2\right)\lor 1\right)}
  \left(1/2-\beta\right)^{2}}\, e^{n\Re\hat{\phi}\left(x\right)}\\
 & = e^{n\Re\hat{\phi}\left(x\right)+o\left(n\right)}.
\end{align*}
It follows that
\begin{equation}
\label{F3bound}
\left|\Hk 3\right|=O\left(ne^{-C_{3}n/2}\right) = o(1).
\end{equation}

Next suppose that $\left(\beta,\theta\right)\in\Hk 4\cap \mathcal{X}_n$.
Then we have $\frac{1}{2}<\beta\le\frac{1}{2}+o(1)$ and $\left(n/2-b+2\right)!=1$. 
By the alternate form of Stirling's inequality and (\ref{lem:technical}),
\begin{align*}
\left|a_{n}\left(\beta,\theta\right)\right| & \le \frac{n^{3}e^{2}}{\sqrt{\left(\beta n\lor 1\right)}}\, c_{n}\left|\frac{\left(2f\left(r_{\beta}e^{i\theta}\right)\right)^{\beta}}{r_{\beta}e^{i\theta}\beta^{\beta}}\right|^{n}\\
& \le e^{o\left(n\right)}c_{n}\left(\frac{\left(2f\left(r_{1/2}\right)\right)^{1/2}}{r_{1/2}\left(\frac12\right)^{1/2}}\right)^{n}.
\end{align*}
By direct computation,
\[
\log\frac{\left(2f\left(r_{1/2}\right)\right)^{1/2}}{r_{1/2}\left(\frac12\right)^{1/2}}=\phi(x^*)+C_4
\]
for some $C_4>0$. It follows that
\begin{equation}
\label{F4bound}
\left|\Hk 4\right|=O\left(ne^{-C_{4}n/2}\right) = o(1).
\end{equation}
It remains to consider $\Hk 1$. Define $\ceil{\left(\beta,\theta\right)}=\left(\frac{\ceil{\beta n}}{n},\theta\right)$,
so that $F_n=n\int_{\mathcal{X}}\hat{a}_{n}\left(\ceil{x}\right)\mathrm{d}x$. For
any $y\in\mathbb{R}^{2}$, define $x_{y}=x^{*}+y/\sqrt{n}$ and
$B_{n}=\left\{ y:\ceil{x_{y}}\in\Xk 1\right\} $, so that we can
make the change of variables
\[
\Hk 1=\int_{B_{n}}\hat{a}_{n}\left(\ceil{x_{y}}\right)\mathrm{d}y.
\]
Note that 
\[
\left|y/\sqrt{n}\right|=\left|x_{y}-x^{*}\right|\le\left|x_{y}-\ceil{x_{y}}\right|+\left|\ceil{x_{y}}-x^{*}\right|=O\left(n^{-1/3}+n^{-1}\right)=O\left(n^{-1/3}\right)
\]
for $y\in B_{n}$, so that $B_{n}$ is approximately a ball of radius
$O\left(n^{1/6}\right)$.

Next, a first-order Taylor expansion of $D\phi$ about $x^{*}$ gives
\[
\left|D\phi\left(x_{y}\right)\right|=O\left(\left|y/\sqrt{n}\right|\right)=O\left(n^{-1/3}\right).
\]
Another first-order Taylor expansion of $\phi$ about $x_{y}$ gives
\[
\phi\left(\ceil{x_{y}}\right)-\phi\left(x_{y}\right)=O\left(\left|D\phi\left(x_{y}\right)\right|\left|\ceil{x_{y}}-x_{y}\right|\right)=O\left(n^{-4/3}\right),
\]
so that 
\begin{equation}
e^{n\phi\left(x_{y}\right)}\sim e^{n\phi\left(\ceil{x_{y}}\right)}\label{eq:ceil-asymptotic}
\end{equation}
uniformly. Now, for each $y\in\mathbb{R}^{2}$, we have $\ceil{x_{y}}\tox^{*}$.
For $n$ large enough so that $y\in B_{n}$, we have $\psi\left(\ceil{x_{y}}\right)\to\psi\left(x^{*}\right)$
by continuity and $e^{n\phi\left(x_{y}\right)}\to e^{\frac{1}{2}y^{T}Hy}$
by (\ref{eq:taylor}). We therefore have $\boldsymbol{1}_{B_{n}}\left(y\right)\hat{a}_{n}\left(\ceil{x_{y}}\right)\to\psi\left(x^{*}\right)e^{\frac{1}{2}y^{T}Hy}$
for all $y$.

Recalling that $C_{2}$ and $2C_{1}$ are bounds involving $\psi$
and $\phi$ respectively, with (\ref{eq:ceil-asymptotic}) we have
$\left|\boldsymbol{1}_{B_{n}}\left(y\right)\hat{a}_{n}\left(\ceil{x_{y}}\right)\right|\le2C_{2}e^{-C_{1}\left|y\right|^{2}}$
for sufficiently large $n$. 
Since
\[
\left(\det\left(-H\right)\right)^{-1/2}  = \frac{1}{3}\sqrt{\frac{2}{7}},\quad
\psi\left(x^{*}\right)  =  108\sqrt{2}
\]
we obtain, by the dominated convergence theorem,
\[
\Hk 1\to\psi\left(x^{*}\right)\int_{\mathbb{R}^{2}}\, e^{\frac{1}{2}y^{T}Hy}\, 
  \mathrm{d}y=2\pi\, \psi\left(x^{*}\right)\, 
     \left(\det\left(-H\right)\right)^{-1/2} = \frac{144\pi}{\sqrt{7}}.
\]
Combining this with (\ref{F2bound}--\ref{F4bound})
completes the proof.
\end{proof}

We now pull these calculations together to prove the following.

\begin{lemma}
Let $d=3$. Then
\[ \mathbb{E}[Y^2] \sim \frac{18}{\sqrt{14}}\, \left(\frac{16}{3}\right)^n,\]
and hence
\[ \frac{\mathbb{E}[Y^2]}{[\mathbb{E} Y]^2} \rightarrow \frac{9}{\sqrt{14}}.
\]
It follows that Condition~\ref{cond:A4} holds when $d=3$.
\label{lem:secondmoment}
\end{lemma}

\begin{proof}
Lemma~\ref{lem:dominates} and (\ref{split}) prove that
\[ F_n \sim \frac{72}{n^3\sqrt{7}}\, \left( 4\sqrt{\frac{2e}{3n}}\right)^n.\]
Substituting $d=3$ into (\ref{eq:E2}) and applying
(\ref{eq:num-pairings-asymptotic}) gives 
\[
\mathbb{E}Y^{2} = \frac{(6\sqrt{2})^n\, (n!)^2}{4\pi\, \#P(3n)}\, F_n
      \sim\frac{18}{\sqrt{14}}\left(\frac{16}{3}\right)^{n},
\]
using Stirling's approximation. 
Then, with (\ref{eq:EY-asymptotic}) and Lemma~\ref{EY2/EY2-target},
we conclude that
\[
\frac{\mathbb{E}Y^{2}}{(\mathbb{E}Y)^{2}}\to\frac{9}{\sqrt{14}}=\exp\left(\sum_{j=1}^{\infty}\lambda_{j}\gamma_{j}^{2}\right).
\]
This establishes Condition~\ref{cond:A4}, as required. 
\end{proof}


We can now complete the proof of Theorem~\ref{thm:distribution-3}.

\begin{proof}[Proof of Theorem \ref{thm:distribution-3}]
We will prove that for general $d \geq 3$, if Condition~\ref{cond:A4}
holds then Conjecture~\ref{cnj:distribution} is true.
In particular, this will prove Theorem~\ref{thm:distribution-3},
using Lemma~\ref{lem:secondmoment}.

Suppose that Condition \ref{cond:A4} is satisfied for some fixed
integer $d\geq 3$. Then by Lemma~\ref{EY2/EY2-target} we may apply
Theorem~\ref{thm:janson} to conclude that
(\ref{eq:W}) holds for $Y$.  Therefore, for all real numbers $y$ we have
\begin{alignat*}{2}
\mathbb{P}\left(Y_{\mathcal{G}}/\mathbb{E} Y_{\mathcal{G}}<y\right) 
  &=\enspace&& \mathbb{P}\left(Y/\mathbb{E} Y_{\mathcal{G}}<y|X_{1}=X_{2}=0\right)\\
 & \rightarrow &&
  \mathbb{P}\left(W\exp\left(\lambda_{1}\zeta_{1}+\lambda_{2}\zeta_{2}\right)<y|
           Z_{1}=Z_{2}=0\right)\\
 &=&& \mathbb{P}\left(\prod_{j=3}^{\infty}\left(1+\zeta_{j}\right)^{Z_{j}}
        e^{-\lambda_{j}\zeta_{j}}<y\right).
\end{alignat*}
Hence Conjecture \ref{cnj:distribution} is a consequence
of \ref{cond:A4}. 
\end{proof}

\subsection{Support for Conjecture~\ref{cnj:distribution}}\label{ss:numerical}
\global\long\defp{p}


Let $p_{d}\left(n\right)$ denote the quotient 
\[
\frac{\mathbb{E} Y^{2}}{\left(\mathbb{E} Y\right)^{2}}\left/\exp\left(\sum_{j=1}^{\infty}\lambda_{j}\zeta_{j}^{2}\right)\right. .
\]
For any fixed integer $d\geq 4$, Conjecture~\ref{cnj:distribution} holds 
if and only Condition~\ref{cond:A4} from Theorem~\ref{thm:janson} is satisfied;
that is, if and only if $p_{d}(n)\sim 1$.
Using Lemma~\ref{lem:|P|EY2} we can compute
$p_{d}\left(100\right)$ for various values of $d$:
\begin{center}
\begin{tabular}{|c|c|c|c|c|c|}
\hline 
$d$ & 3 & 4 & 5 & 6 & 100\tabularnewline
\hline 
$p_{d}\left(100\right)$ & 0.9761 & 0.9881 & 0.9921 & 0.9942 & 0.9998\tabularnewline
\hline 
\end{tabular}
\par\end{center}
Figure \ref{fig:EY2-plot} is a plot of
$p_{d}\left(n\right)$ for $d\in \{ 3,4,5,6,100\}$ and $n\leq 50$.  

\begin{figure}[ht!]
\begin{center}
\centerline{\includegraphics{spanningtree-arxiv-figure0}}
\caption{\label{fig:EY2-plot} A plot of $p_d(n)$ for $d \in \{ 3,4,5,6,100\}$}
\end{center}
\end{figure}



  
	
 

	
	
	
		
		
	
		
		
	
		
		
	
		
		
	
		
		

This plot supports our conjecture that
$p_{d}(n)\sim 1$ for all $d\geq 4$. 
Indeed, the rate of convergence to 1 appears to increase as $d$ increases.

We now give an asymptotic result which is equivalent to 
Conjecture~\ref{cnj:distribution}.
Combining (\ref{eq:num-pairings-asymptotic}), \ref{cond:A4}, (\ref{eq:EY-asymptotic}), Lemma~\ref{EY2/EY2-target}, Lemma~\ref{lem:|P|EY2} and applying Stirling's 
formula shows that for a fixed integer $d\geq 4$,
Conjecture~\ref{cnj:distribution} holds if and only if
\begin{align}
\sum_{b=1}^{n}\,&
  \frac{2^{b}}{b!\,\left(\left(d/2-1\right)n-b+2\right)!}\,
  \left[z^{n}\right]\,\left(\sum_{j=1}^{\infty}\binom{(d-1)j}{j}\, z^{j}\right)^{b} \notag\\
&\sim
  \frac{2d^{2}}{\pi\left(d-2\right)^{4} \, n^3}\,\,
  \sqrt{\frac{2d-2}{d^{2}-d+1}}\,\,
   \left(\frac{\left(d-1\right)^{2(d-1)}}{\left(d-2\right)^{2(d-2)}}\left(\frac{2e}{d n}\right)^{d/2-1}\right)^{n}.
\label{sufficient}
\end{align}


\section*{Acknowledgements} 
We would like to thank the referee for their helpful comments.


2014-02-19T02:02:34.000
1309
1309.6710
en
https://arxiv.org/abs/1309.6710
Let $d \geq 3$ be a fixed integer. We give an asympotic formula for the expected number of spanning trees in a uniformly random $d$-regular graph with $n$ vertices. (The asymptotics are as $n\to\infty$, restricted to even $n$ if $d$ is odd.) We also obtain the asymptotic distribution of the number of spanning trees in a uniformly random cubic graph, and conjecture that the corresponding result holds for arbitrary (fixed) $d$. Numerical evidence is presented which supports our conjecture.
Combinatorics (math.CO); Probability (math.PR)
On the number of spanning trees in random regular graphs
Qwen/Qwen-72B
1. YES
2. YES