If a heuristic consideration is acceptable, we can do the following.
$$I(n)=\int_1^\infty e^{-t}\ln^nt\,dt\overset{t=nx}{=}n\int_{1/n}^\infty e^{-n\big(x-\ln(\ln n+\ln x)\big)}dx=n\int_{1/n}^\infty e^{-nf(x,n)}dx$$
Though the function $f(x,n)=x-\ln(\ln n+\ln x)$ depends on $n$, we can try to apply the Laplace method. Logarithm is a very slowly changing function, and the extremum of $f(x,n)$ lies nearly $x=\frac1{\ln n};$ (can be easily checked by iterations).
For $n\gg1$ the point $x=\frac1{\ln n}$ lies far from the end of the interval of integration ($x=\frac1n$), so we may suppose that the Laplace' method is applicable (checking afterwards that we get a reasonable approximation).
$$\frac{d}{dx}f(x,n)=1-\frac1{x\ln(xn)}=0\,\,\Rightarrow\,\,x\ln(xn)=1\,\,\Rightarrow\,\,\frac1xe^\frac1x=n$$
so the extremum (maximum) point is $x_0=\frac1{W(n)}$, where $W$ is the Lambert's function, defined as $\,We^W=n$
$$f''(x,n)=\frac{1+\ln(xn)}{x^2\ln^2(xn)}; \,f(x_0;n)=W(n)+1$$
$$I(n)\approx ne^{-n\left(\frac1{W(n)}-\ln W(n)\right)}\int_{1/n}^\infty e^{-\frac n2\big(W(n)+1\big)\left(x-\frac1{W(n)}\right)^2}dx$$
$$\approx ne^{-n\left(\frac1{W(n)}-\ln W(n)\right)}\int_{-\infty}^\infty e^{-\frac n2\big(W(n)+1\big)t^2}dt$$
$$\boxed{\,\,I(n)\approx\sqrt{\frac{2\pi n}{W(n)+1}}e^{n\left(\ln W(n)-\frac1{W(n)}\right)}=\sqrt{\frac{2\pi n}{W(n)+1}}(W(n))^ne^{-e^{W(n)}}\,\,}$$
For the rough approximation we can take $W(n)\approx \ln n$ (for example, for $W(n)=5,\, n=742.066$, and $\ln (742.066)=6.609...$), though to get higher accuracy we should use the exact value of $W(n)$ for every $n$ (or vice versa - evaluate $n$ for every $W$).
Numeric check with WolframAlpha
- $W=2, \, n=14.778,\, I=98.4065,\, approximation =95.5954$
- $W=3, \, n=60.26,\, I=1.044\times 10^{21},\, approximation =1.037\times 10^{21}$
- $W=4, \, n=218.2926,\, I=9.86\times 10^{108},\, approximation =9.84\times 10^{108}$
- $W=5, \, n=742.066,\, I=4.7033\times 10^{455},\, approximation =4.6997\times 10^{455}$