2
$\begingroup$

How to extract the divergent part of the following integral simply as $u \rightarrow \infty$

$$g(u) = \frac{\sqrt{2u}}{\pi} \int^1_{\frac{1}{u}} dz \frac{\sqrt{z-1}}{\sqrt{z^2-u^{-2}}} $$

$\endgroup$
1
  • 1
    $\begingroup$ Have you tried expanding in power series? $\endgroup$
    – S. Carnahan
    Nov 30, 2013 at 12:32

3 Answers 3

4
$\begingroup$

it's an elliptic integral; a series expansion gives

$$g(u)=\frac{\sqrt{2u}}{\pi} \int^1_{1/u} dz \frac{\sqrt{z-1}}{\sqrt{z^2-u^{-2}}}=i\frac{1}{\pi}(2u)^{1/2}\;[\ln (8u)-2]+{\cal O}(u^{-1/2})$$

so the integral diverges as $\sqrt{u}\ln u$

here is a plot of $-i(\pi/\sqrt{2u})g(u)$, evaluated numerically, and $\ln(8u)-2$ versus $u$, just as a check:

quite a fast convergence

$\endgroup$
1
  • $\begingroup$ I checked it numerically, it seems correct as stated (so $\ln(8u)-2$). $\endgroup$ Dec 1, 2013 at 14:19
3
$\begingroup$

$\newcommand{\+}{^{\dagger}}% \newcommand{\angles}[1]{\left\langle #1 \right\rangle}% \newcommand{\braces}[1]{\left\lbrace #1 \right\rbrace}% \newcommand{\bracks}[1]{\left\lbrack #1 \right\rbrack}% \newcommand{\ceil}[1]{\,\left\lceil #1 \right\rceil\,}% \newcommand{\dd}{{\rm d}}% \newcommand{\down}{\downarrow}% \newcommand{\ds}[1]{\displaystyle{#1}}% \newcommand{\equalby}[1]{{#1 \atop {= \atop \vphantom{\huge A}}}}% \newcommand{\expo}[1]{\,{\rm e}^{#1}\,}% \newcommand{\fermi}{\,{\rm f}}% \newcommand{\floor}[1]{\,\left\lfloor #1 \right\rfloor\,}% \newcommand{\half}{{1 \over 2}}% \newcommand{\ic}{{\rm i}}% \newcommand{\iff}{\Longleftrightarrow} \newcommand{\imp}{\Longrightarrow}% \newcommand{\isdiv}{\,\left.\right\vert\,}% \newcommand{\ket}[1]{\left\vert #1\right\rangle}% \newcommand{\ol}[1]{\overline{#1}}% \newcommand{\pars}[1]{\left( #1 \right)}% \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}} \newcommand{\pp}{{\cal P}}% \newcommand{\root}[2][]{\,\sqrt[#1]{\,#2\,}\,}% \newcommand{\sech}{\,{\rm sech}}% \newcommand{\sgn}{\,{\rm sgn}}% \newcommand{\totald}[3][]{\frac{{\rm d}^{#1} #2}{{\rm d} #3^{#1}}} \newcommand{\ul}[1]{\underline{#1}}% \newcommand{\verts}[1]{\left\vert\, #1 \,\right\vert}$ $$ {\rm g}\pars{u} = {\root{2u} \over \pi}\int_{1/u}^{1}{\root{z - 1} \over \sqrt{z^{2} - u^{-2}}}\,\dd z = {\root{2u} \over \pi}\int_{1/u}^{1}\root{z - 1}\varphi'\pars{z}\,\dd z $$ where $\ds{\varphi\pars{z} \equiv \int_{1/u}^{z}{\dd t \over \root{t^{2} - u^{-2}}}}$

Then, \begin{align} {\rm g}\pars{u} &= -\,{\root{2} \over 2\pi}\,u^{1/2}\int_{1/u}^{1}\, {\varphi\pars{z} \over \root{z - 1}}\,\dd z \end{align} With the change of variables $\ds{t \equiv u^{-1}\sec\pars{x}}$: \begin{align} \varphi\pars{z} &=\int_{0}^{\arccos\pars{u^{-1}z^{-1}}}\sec\pars{x}\,\dd x =\left.\ln\pars{\sec\pars{x} + \tan\pars{x}}\right\vert_{0}^{\arccos\pars{u^{-1}z^{-1}}} \\[3mm]&=\left.\ln\pars{ut + \root{u^{2}t^{2} - 1}}\right\vert_{u^{-1}}^{z} \\[3mm]&=\ln\pars{uz + \root{u^{2}z^{2} - 1}} =\ln\pars{u} + \ln\pars{z + \root{z^{2} - u^{-2}}} \end{align}

The 'leading term' when $u \gg 1$ becomes: $$ {\rm g}\pars{u}\sim -\,{\root{2} \over 2\pi}\,u^{1/2}\ln\pars{u} \int_{u^{-1}}^{1}{\dd z \over \root{z - 1}} = {\root{2} \over \pi}\,u^{1/2}\ln\pars{u}\root{u^{-1} - 1} $$

$\endgroup$
2
$\begingroup$

The Maple command $$series(sqrt(z-1)/sqrt(z^2-1/u^2), z = 1/u, 2) assuming u>0$$ produces $$ 1/2\,\sqrt {-{\frac {-1+u}{u}}}\sqrt {2}\sqrt {u}{\frac {1}{\sqrt {z-{ u}^{-1}}}}+ \left( -1/8\,\sqrt {-{\frac {-1+u}{u}}}\sqrt {2}{u}^{3/2}- 1/4\,\sqrt {-{\frac {-1+u}{u}}}\sqrt {2}{u}^{3/2} \left( -1+u \right) ^{-1} \right) \sqrt {z-{u}^{-1}}+O \left( \left( z-{u}^{-1} \right) ^ {3/2} \right) .$$

$\endgroup$
3
  • $\begingroup$ The Maple command $$asympt(int(sqrt(-z+1)/sqrt(z^2-1/u^2), z = 1/u .. 1), u, 1)assuming u>0 $$ produces $$3\,\ln \left( 2 \right) +\ln \left( u \right) +O \left( {u}^{-1} \right) .$$ $\endgroup$
    – user64494
    Dec 1, 2013 at 16:04
  • $\begingroup$ so it misses the extra term $-2$ ? $\endgroup$ Dec 1, 2013 at 19:39
  • $\begingroup$ This is included in $O(u^{-1})$. $\endgroup$
    – user64494
    Dec 1, 2013 at 20:32

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.