Suppose $Z , \epsilon \sim N(0, 1)$ are independent Gaussian random variables. Let $a \ll 1$ be a small positive number. Let $W = aZ + \epsilon$. It can be show that 
\begin{align}
\mathbb{E} [ W^2 (Z^2 - 1)] = 2 a^2.
\end{align}
Now suppose I truncate random variable $W$.
I was wondering what truncation level $R$ should be such that 
\begin{align}
\mathbb{E} [ W^2 \cdot \mathbf{1}\{ | W| > R\} \cdot (Z^2 - 1)] \leq a^2.
\end{align}
Can we set $R$ to be a constant?
This seems an easy problem, but computing the expectation seems not trivial. 



Expectation of Truncated Bivariate Gaussian Random Variables
There is a simple but often efficient trick to facilitate such computations. Let $H\_k$ be the normalized Hermite polynomials and let $\Phi,\Psi$ be any two functions square integrable with respect to the Gaussian measure $\gamma$. Write $\Phi=\sum\_k\varphi\_kH\_k$, $\Psi=\sum\_k\psi\_kH\_k$. Then, if $Y,Z$ are jointly Gaussian with variance $1$ and covariance $\alpha$, we have
$$
E(\Phi(Y)\Psi(Z))=\sum\_k \varphi\_k\psi\_k\alpha^k\,.
$$
In your case $\Psi$ is essentially $H\_2$ (up to normalization) and $\Phi\_r(x)=x^2\chi\_{[-r,r]^c}(x)$, so putting $Y=(a^2+1)^{-1/2}W$, we see that we just need to drop the integral $2\varphi\_2(r)=\int\_{\mathbb R} \Phi\_r(x)(x^2-1)\,d\gamma(x)$ from its initial value $2$ at $r=0$ to $1$. That corresponds to some fixed level $r\_0>0$ and the answer to the original problem is $R=r\_0\sqrt{a^2+1}$. 


