[step:Average the Hermitian form and choose the Jucys-Murphy eigenbasis]
Let $(\cdot,\cdot)_0:S^\lambda\times S^\lambda\to\mathbb C$ be any positive definite Hermitian form. Define a Hermitian form $(\cdot,\cdot):S^\lambda\times S^\lambda\to\mathbb C$ by
\begin{align*}
(u,v):=\frac{1}{|S_n|}\sum_{\sigma\in S_n}(\sigma u,\sigma v)_0.
\end{align*}
For $u\ne 0$, the summand with $\sigma=1_{S_n}$ is $(u,u)_0>0$ and every summand $(\sigma u,\sigma u)_0$ is non-negative, so $(\cdot,\cdot)$ is positive definite. If $\tau\in S_n$, then right multiplication by $\tau$ is a bijection $S_n\to S_n$, and therefore
\begin{align*}
(\tau u,\tau v)=\frac{1}{|S_n|}\sum_{\sigma\in S_n}(\sigma\tau u,\sigma\tau v)_0=(u,v).
\end{align*}
Thus every group element acts unitarily on $S^\lambda$.
For each transposition $(j\ k)$, the associated operator is self-adjoint because it is unitary and equal to its inverse. Hence each $X_k$ is self-adjoint, being a real linear combination of self-adjoint transposition operators. By [[Commutativity of the Jucys-Murphy Elements](/theorems/8446)][citetheorem:8446], the operators $X_1,\dots,X_n$ commute. The finite-dimensional spectral theorem for commuting [self-adjoint operators](/page/Self-Adjoint%20Operators) therefore applies to the Hermitian space $S^\lambda$ and gives an orthonormal simultaneous eigenbasis.
We now use the standard Young-Jucys-Murphy spectral theorem for complex Specht modules in the following exact form: on $S^\lambda$, the joint eigenspaces of $X_1,\dots,X_n$ are precisely the one-dimensional spaces
\begin{align*}
E_T:=\{u\in S^\lambda:X_ku=c_k(T)u\text{ for every }1\le k\le n\}
\end{align*}
indexed by $T\in\operatorname{SYT}(\lambda)$, and these spaces exhaust $S^\lambda$. This is the content of [[Content Eigenvalues of the Jucys-Murphy Elements](/theorems/8447)][citetheorem:8447]. Its hypotheses are satisfied here because the module is the complex Specht module $S^\lambda$ of shape $\lambda$, the elements $X_k$ are the Jucys-Murphy elements with convention $X_k=\sum_{1\le j<k}(j\ k)$, and the content convention is $c_k(T)=\operatorname{col}_T(k)-\operatorname{row}_T(k)$. Choose a unit vector $v_T\in E_T$ for every $T\in\operatorname{SYT}(\lambda)$. The spaces $E_T$ are mutually orthogonal because the $X_k$ are commuting self-adjoint operators, and they exhaust $S^\lambda$; hence $(v_T)_{T\in\operatorname{SYT}(\lambda)}$ is an [orthonormal basis](/page/Orthonormal%20Basis) satisfying
\begin{align*}
X_kv_T=c_k(T)v_T.
\end{align*}
[/step]