[step:Separate the null edge from the planted outlier above the BBP threshold]
Assume $\theta>\sqrt{\gamma}$. By the null Marchenko--Pastur upper-edge conclusion in the real Wishart version of the [BBP transition for spiked covariance matrices](/page/BBP%20Transition), equivalently the zero-spike case of that theorem, the hypotheses are satisfied under $\mathbb{P}_{0,n}$ because the columns are independent real Gaussian vectors with covariance $I_{p_n}$ and $p_n/n\to\gamma\in(0,\infty)$. This formulation covers every positive aspect ratio, including $\gamma>1$, where $S_n$ may have zero eigenvalues because $p_n>n$ asymptotically; the upper spectral edge is still $(1+\sqrt{\gamma})^2$. Hence, under $\mathbb{P}_{0,n}$,
\begin{align*}
\lambda_{\max}(S_n) \to (1+\sqrt{\gamma})^2
\end{align*}
in probability.
The supercritical rank-one outlier conclusion in the same real Wishart [BBP transition for spiked covariance matrices](/page/BBP%20Transition) applies conditionally on $v_n$ because the covariance matrix $I_{p_n}+\theta v_nv_n^\top$ has one population eigenvalue $1+\theta$ in the direction $v_n$ and all remaining population eigenvalues equal to $1$; the spike strength $\theta$ is fixed, the spike rank is one, the aspect ratio satisfies $p_n/n\to\gamma\in(0,\infty)$, and the supercritical hypothesis is $\theta>\sqrt{\gamma}$. Let $O(p_n):=\{Q\in\mathbb{R}^{p_n\times p_n}:Q^\top Q=I_{p_n}\}$ denote the real orthogonal group, let $e_1:=(1,0,\dots,0)\in\mathbb{R}^{p_n}$ denote the first coordinate vector, and choose an orthogonal matrix $Q_n\in O(p_n)$ such that $Q_nv_n=e_1$. The conditional law of the eigenvalues of $S_n$ does not depend on the particular value of $v_n$: under this change of coordinates, $Q_nX_n$ has covariance $I_{p_n}+\theta e_1e_1^\top$ and $S_n(Q_nX_n)=Q_nS_n(X_n)Q_n^\top$, which has the same eigenvalues as $S_n(X_n)$. Hence the conditional convergence holds with the same deterministic limit for every $v_n$, and integrating over the spherical prior transfers the convergence to the mixture law $\mathbb{P}_{1,n}$. Therefore, under $\mathbb{P}_{1,n}$,
\begin{align*}
\lambda_{\max}(S_n) \to (1+\theta)\left(1+\frac{\gamma}{\theta}\right)
\end{align*}
in probability.
The two deterministic limits are strictly separated. Indeed,
\begin{align*}
(1+\theta)\left(1+\frac{\gamma}{\theta}\right)-(1+\sqrt{\gamma})^2 = \frac{(\theta-\sqrt{\gamma})^2}{\theta}>0.
\end{align*}
Choose a threshold $\tau$ satisfying
\begin{align*}
(1+\sqrt{\gamma})^2 < \tau < (1+\theta)\left(1+\frac{\gamma}{\theta}\right).
\end{align*}
[/step]