Probably not how you want this done, but a useful and elegant method nonetheless:
Note that $$f(x)=\frac{2x}{1+x^2}$$ has the property that $f(x)<1$ except at $f(1)=1$. Since exponentiation to the power of $n\to\infty$ kills off anything less than $1$, we expect that the dominant contribution is near $1$.
$f^n \equiv \exp(n \log f)$ and $n\log(f(1+u))=-n\frac{1}{2}u^2+n\mathcal O(u^3)$. Then letting $v=\sqrt{n} u$,
$$\int_0^\infty n^{1/2}\mathrm d x \exp(-\frac{1}{2}v^2+\mathcal O(n^{-1/2})) \approx \int_{-\infty}^\infty \mathrm d v\, e^{-v^2/2} = \sqrt{2\pi}$$
Verification: WolframAlpha.
Edit: This is called Laplace's method amongst other things, and is a core part of the general theory of asymptotics of integrals. (It's also how Stirling's approximation is most easily derived! Here one applies the method to the $\Gamma$ function.) It can be tightened up and made rigorous straightforwardly by estimating errors. You can also derive an asymptotic expansion by considering the correction terms.
One useful insight it offers is why the $\sqrt{n}$ is really there: it's because the peak contributing the dominant term has a width decaying as $1/\sqrt{n}$.