Definamos (para $\Re(a)>-1$ ) la función $$ f(a):=\int_{1}^{\infty}x^a\,e^{-x}\left(x\right)\,\mathrm{d}x$$ entonces este es el función gamma incompleta $\,f(a)=\Gamma(a+1,1)$ .
Pero como $\,\displaystyle x^a=e^{a\ln(x)}\,$ tenemos : $$f''(a)=\int_{1}^{\infty}x^a\;\ln^{2}(x)\;e^{-x}\,\mathrm{d}x$$
De modo que su función es la segunda derivada de $\,\Gamma(a+1,1)\,$ en $\,a=0$ (para que $x^a$ desaparece) o la segunda derivada relativa a $a$ de $\,\Gamma(a,1)\,$ en $\,a=1$ : $$f''(0)=\lim_{a\to 1}\Gamma(a,1)''$$ (esta respuesta también la dio el Sr. G en los comentarios)
Añadamos que las derivadas no son mucho más sencillas por lo que probablemente sea mejor mantenerlo así...
Por supuesto, como señala Eric Naslund (+1), la respuesta es más sencilla cuando el límite inferior es $0$ porque en este caso obtenemos simplemente $\;\lim_{a\to 0}\Gamma(a+1,0)''=\Gamma''(1)$ .