Upper Bound on a en exponential function I am trying to upper bound the following function and find its growth rate:
\begin{equation}
\psi(y) \stackrel{\triangle}{=} \int_{0}^{\infty}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx,~y>0,
\end{equation}
where $f(x)>0$ satisfies $\int_{0}^{+\infty}f(x)\,dx = 1$ (it is a pdf of the random variable $X$). I would like to find the growth rate of $\psi(y)$ for large values of $y$, i.e., would like to know if $\psi(y)$ grows like $O(y^{-k})$, for some $k>1$.  
Here are my attempts so far: 
It is straightforward to show that the function $\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)\in[0,1]$ is increasing in $x$ for $y>x>0$ and is decreasing in $x$ for $0<y<x$.
Then one can break the bounds of the integral to $[0,y/2], \, [y/2,3y/2], \, [3y/2, y^{3/2}],\,[y^{3/2},\infty)$ and then evaluate each integral. Doing so will give the following:
\begin{align}
\int_{0}^{y/2}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &= O(y^{-k}), \,k>1\\
\int_{y/2}^{3y/2}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &\stackrel{?}{=}O(y^{-k}),\, k>1\\
\int_{3y/2}^{y^{3/2}}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &= O(y^{-k}), \,k>1\\
\int_{y^{3/2}}^{\infty}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)f(x)\,dx &= O(y^{-k}), \,k>1
\end{align}
Therefore, if we can prove that the second integral also grows like $O(y^{-k}),\,k>1$, then we are done. However, I have not been able to show this and am stuck. Any help would be highly appreciated!

First, let us see an example in which $\lim_{y\to \infty} y\psi(y) = \infty$.
Let (log-Cauchy distribution) 
$$f(x) = \frac{1}{\pi x (1 + (\ln x)^2)}, \ x > 0.$$
We have, for sufficiently large $y$, 
\begin{align}
\psi(y) &= \int_{0}^{\infty}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)
\frac{1}{\pi x (1 + (\ln x)^2)}\,\mathrm{d} x\\
&\ge \int_{y - \sqrt[4]{y}}^{y + \sqrt[4]{y}}\exp\left(-\frac{(y-x)^2}{2(\sigma_0^2 + \sigma_1^2 x)}\right)
\frac{1}{\pi x (1 + (\ln x)^2)}\,\mathrm{d} x \\
&\ge \int_{y - \sqrt[4]{y}}^{y + \sqrt[4]{y}}\exp\left(-\frac{\sqrt{y}}{2(\sigma_0^2 + \sigma_1^2 y/2)}\right)
\frac{1}{\pi x (1 + (\ln x)^2)}\,\mathrm{d} x \\
&\ge \int_{y - \sqrt[4]{y}}^{y + \sqrt[4]{y}} \frac{1}{2}\cdot
\frac{1}{\pi x (1 + (\ln x)^2)}\,\mathrm{d} x \\
&= \frac{\arctan(\ln(y + \sqrt[4]{y})) - \arctan(\ln(y - \sqrt[4]{y}))}{2\pi}\\
&= \frac{1}{2\pi}
\arctan \frac{\ln(y + \sqrt[4]{y}) - \ln(y - \sqrt[4]{y})}{1 + \ln(y + \sqrt[4]{y})\ln(y - \sqrt[4]{y})}
\end{align}
where we have used $\tan (x-y) = \frac{\tan x - \tan y}{1 + \tan x \tan y}$.
Thus, we have
\begin{align}
\lim_{y\to \infty} y\psi(y) &= \lim_{y\to \infty} \frac{1}{2\pi} y
\arctan \frac{\ln(y + \sqrt[4]{y}) - \ln(y - \sqrt[4]{y})}{1 + \ln(y + \sqrt[4]{y})\ln(y - \sqrt[4]{y})}\\
&= \lim_{y\to \infty} \frac{1}{2\pi} y \frac{\ln(y + \sqrt[4]{y}) - \ln(y - \sqrt[4]{y})}{1 + \ln(y + \sqrt[4]{y})\ln(y - \sqrt[4]{y})}\\
&= \infty
\end{align}
where we have used $\lim_{z \to 0}\frac{\arctan z}{z} = 1$.
Second, if $\mathbb{E}[X]$ and $\mathbb{E}[X^2]$ are both finite, one can prove that $\psi(y) = O(y^{-2})$. 
