[step:Integrate the squared likelihood and identify the chi-square divergence]
By the definition of the [chi-square divergence](/page/Chi-Square%20Divergence) for $Q_\pi\ll P$ and the [Radon--Nikodym derivative](/page/Radon-Nikodym%20Derivative) identity $L_\pi=dQ_\pi/dP$,
\begin{align*}
1+\chi^2(Q_\pi\|P)
=
\int_{\mathcal X} L_\pi(x)^2\,dP(x).
\end{align*}
Substituting the identity from the previous step, define $G:\Theta\times\Theta\times\mathcal X\to[0,\infty]$ by
\begin{align*}
G(\theta,\theta',x)=L_\theta(x)L_{\theta'}(x).
\end{align*}
The map $G$ is $\mathcal T\otimes\mathcal T\otimes\mathcal A$-measurable because $(\theta,\theta',x)\mapsto(\theta,x)$ and $(\theta,\theta',x)\mapsto(\theta',x)$ are measurable coordinate projections, $(\theta,x)\mapsto L_\theta(x)$ is jointly measurable, and multiplication $[0,\infty]\times[0,\infty]\to[0,\infty]$ is Borel-measurable with the extended non-negative convention. Applying [Tonelli's theorem](/page/Tonelli%20Theorem) to this non-negative measurable function on $(\Theta\times\Theta\times\mathcal X,\mathcal T\otimes\mathcal T\otimes\mathcal A, \pi\otimes\pi\otimes P)$ gives
\begin{align*}
1+\chi^2(Q_\pi\|P)=\int_{\mathcal X}\int_{\Theta\times\Theta} L_\theta(x)L_{\theta'}(x)\,d(\pi\otimes\pi)(\theta,\theta')\,dP(x).
\end{align*}
The same application of Tonelli's theorem exchanges the order of integration, so
\begin{align*}
\int_{\mathcal X}\int_{\Theta\times\Theta} L_\theta(x)L_{\theta'}(x)\,d(\pi\otimes\pi)(\theta,\theta')\,dP(x)=\int_{\Theta\times\Theta}\int_{\mathcal X} L_\theta(x)L_{\theta'}(x)\,dP(x)\,d(\pi\otimes\pi)(\theta,\theta').
\end{align*}
Since independent prior draws $\theta,\theta'\sim\pi$ have joint law $\pi\otimes\pi$, the last integral is $\mathbb E_{\theta,\theta'}\mathbb E_P[L_\theta L_{\theta'}]$.
This proves the stated identity.
[/step]