 Solution

Use the Fourier transform convention $\hat{f}(s)={\int }_{\mathbb{R}}{e}^{-isx}f(x)dx$. First observe that $f(x)=(2\pi {)}^{-1/2}{e}^{-x^2/2}$ satisfies the differential equation, decay, and normalization, using the Gaussian integral. The uniqueness allowed in the question identifies it with $f$. In particular it is a Schwartz function, so the derivative and multiplication rules below are justified.

Taking the Fourier transform of the differential equation givesThe Riemann-Lebesgue lemma gives decay of $\hat{f}$. Every decaying solution $h$ of this equation is proportional to $e^{-s^2/2}$: setting $w=h^{\prime }+sh$ gives $w^{\prime }=sw$, and integration yieldsThe term with $C\ne 0$ grows like $C{e}^{s^2/2}/(2s)$ as $s\to +\infty$, so decay forces $C=0$. Thus $\hat{f}=\lambda f$ for a scalar $\lambda$. This establishes proportionality before using Fourier inversion to calculate the scalar.