I have an equation which is similar to Hermite's differential equation $y''(x)+xy'(x) +\lambda y(x)=0$ with decaying boundary conditions at infinity $y\rightarrow 0$ as $x\rightarrow\pm\infty$. I believe the eigenvalues take the form $\lambda_n =-n$ where $n=0,1,2,...$ from the power series solution.
The eigenfunction which satisfies the boundary condition is $\phi_n(x) = A_n e^{-y^2/2}H_{\lambda_n-1}\left(\frac{y}{\sqrt 2}\right)$. All of which are Hermite polynomials for negative integers. The recurrence relation $2 n H_{n-1}(x) = 2x H_n(x)-H_{n+1}(x)$ is no helpful since it gives indeterminate values.
But I have seen them that these polynomial of non-negative orders are expressed in terms of error function, but I could not find any reference to this. Is there any way to derive them in terms of error function? It would be helpful if someone clarifies about the orthogonality condition for negative integers. Is it the same as Hermite's polynomials of positive integers?
For $H(-n,z)$ there is the rather complicated expression with the error function and $H(n,z)$ from http://functions.wolfram.com/07.01.03.0014.01 . It seems simpler to use Kummer’s confluent hypergeometric function, see http://functions.wolfram.com/07.01.02.0001.01
$$H(\nu, z) = 2^{\nu} \sqrt{\pi}\Big(\frac{1}{\Gamma((1 - \nu)/2)} {_1F_1}(-\nu/2, 1/2, z^2) -\frac{2z}{\Gamma(-\nu/2)}{_1F_1}((1 - \nu)/2, 3/2, z^2))\Big)$$
I dont know about orthogonality for negative indices.