This is a perfect example of a question-begging mathematical "proof": The kind of "proof" you are trying to construct here is really not terribly useful. If you don't know how to prove the integral:
$$\int \limits_0^\infty y^{-1/2} e^{-y} dy = \sqrt{\pi},$$
then you need to ask yourself why are you willing to assume that you already know the integral:
$$\int \limits_{-\infty}^\infty y^2 e^{-\tfrac{1}{2} y^2} dy = \sqrt{2 \pi},$$
which is what you are assuming when you appeal to the fact that a standard normal random variable has a variance of one. In this particular context, the properties of the normal density function first had to be derived by elementary methods, that are essentially equivalent to the result you are trying to prove. Indeed, even the fact that the normal density is a density function (i.e., that it integrates to one) is essentially just a restatement of an integral that is a simple variation of the integral you are trying to prove.
Appealing to the properties of the normal density to prove the gamma integral is an example of a "proof" that tries to demonstrate result $A$ by appealing to result $B$, where the latter theorem is a result that is itself proved by the same technique as an elementary proof of $A$. In these cases you need to be very careful because the latter result may even have a "proof" that hinges on the former result, in which case you get a completely circular argument.
Proof by elementary methods: Fortunately, in this particular case, the gamma integral and the resulting properties of the normal distribution are all provable by elementary methods. Rather than trying to "prove" the gamma integral by appealing to what are essentially variations of the same result, you can use Liouville's method (hat tip to whuber for mentioning it in comments), which uses polar coordinates to establish that:
$$\begin{equation} \begin{aligned}
\Bigg( \int \limits_{-\infty}^\infty e^{- \frac{1}{2} y^2} dy \Bigg)^2
&= \Bigg( \int \limits_{-\infty}^\infty e^{- \frac{1}{2} x^2} dx \Bigg) \Bigg( \int \limits_{-\infty}^\infty e^{- \frac{1}{2} y^2} dy \Bigg) \\[6pt]
&= \int \limits_{-\infty}^\infty \int \limits_{-\infty}^\infty e^{- \frac{1}{2} (x^2+y^2)} dx dy \\[6pt]
&= \int \limits_0^{2 \pi} \int \limits_{0}^\infty r e^{- \frac{1}{2} r^2} dr d\theta \\[6pt]
&= \int \limits_0^{2 \pi} \Bigg[ - e^{- \frac{1}{2} r^2} \Bigg]_{r=0}^{ r \rightarrow\infty} d\theta \\[6pt]
&= \int \limits_0^{2 \pi} \Bigg[ 0 - (-1) \Bigg] d\theta \\[6pt]
&= \int \limits_0^{2 \pi} d\theta \\[6pt]
&= 2 \pi. \\[6pt]
\end{aligned} \end{equation}$$
Using the substitution $r = \tfrac{1}{2} y^2$ gives $dr = y dy$, so you then have:
$$\begin{equation} \begin{aligned}
\Gamma(\tfrac{1}{2})
&= \int \limits_0^\infty r^{-1/2} e^{-r} dr \\[6pt]
&= \int \limits_0^\infty \sqrt{\frac{2}{y^2}} e^{- \tfrac{1}{2} y^2} y \ dy \\[6pt]
&= \sqrt{2} \int \limits_0^\infty e^{- \tfrac{1}{2} y^2} \ dy \\[6pt]
&= \frac{1}{\sqrt{2}} \int \limits_{-\infty}^\infty e^{- \tfrac{1}{2} y^2} \ dy \\[6pt]
&= \frac{1}{\sqrt{2}} \cdot \sqrt{2 \pi} = \sqrt{\pi}. \\[6pt]
\end{aligned} \end{equation}$$
This method uses only substituions and familiarity with polar coordinates, and most importantly, it does not require you to assume any premise that is effectively just a variation of the result.