Covariance between squared and exponential of Gaussian random variables

Assuming the random vector $[X \ \ Y]’$ follows a bivariate Gaussian distribution with mean $[\mu_X \ \ \mu_Y]’$, and covariance matrix $ \left[ \begin{array}{cc}
\sigma_X ^2 & \sigma_Y \ \sigma_X \ \rho \\
\sigma_Y \ \sigma_X \ \rho & \sigma_Y^2
\end{array} \right]$, I am looking for the expression of $\text{Cov}(X^2, \exp{Y})$.

I know there exists a formula (Stein’s lemma) for $\text{Cov}(X, \exp{Y})$ and $\text{Cov}(X^2, {Y})$, but I did not manage to find a closed form formula for $\text{Cov}(X^2, \exp{Y})$ given the two transformations.

Solutions Collecting From Web of "Covariance between squared and exponential of Gaussian random variables"