4
$\begingroup$

I have an integral of the form

$$\int\nolimits^\infty_{-\infty}\mathrm d \omega \, \frac{\omega^2}{k^2 + \gamma^2 \omega^2}$$

which diverges. This integral should have a finite value, as it must related to some physical measurement. I am trying to assign a value to the integral, kind of like how one does using regularisation. In a few papers on theoretical physics (Which is the field I am in), I have seen people use the Cauchy principal value in the form

$$-\!\!\!\!\!\!\!\int^\infty_{-\infty}\mathrm dx \, f(x) = \lim_{L \to \infty} \, \frac1{L} \int^L_{-L}\mathrm dx \, f(x)$$

but I am not sure how one deduces that from

http://en.wikipedia.org/wiki/Cauchy_principal_value

  • 1
    Dividing by L seems to be an invention of your own and makes the result, if it exists, behave very differently than any (reasonably defined) integral on the real line. In any case this is not Cauchy principal value (and the only reasonable value for the integral you are interested in is $+\infty$).2011-09-17
  • 0
    I don't know if you access to this article, but http://prl.aps.org.ez.sun.ac.za/pdf/PRL/v89/i14/e144101 does something to this effect. Sorry, maybe I got confused by terminology, they use "principal value" in quotes, so it seems to be a theoretical physics invention perhaps? In any case, is there a sensible way to assign a value to said integral? $\omega$ is a frequency, so of course it cannot really become infinite, but there needs to be some finite cutoff...2011-09-17
  • 1
    Karl, *Off-campus access is restricted to Stellenbosch University staff and currently registered students only*... Anyway, the limit you want to consider is more a *mean value* of $f$ than an integral.2011-09-17
  • 0
    Sorry, copy/paste the wrong thing: http://prl.aps.org/pdf/PRL/v89/i14/e144101 . Anyways, a mean value might be what I want. Thanks. I will think a bit more about this.2011-09-17
  • 0
    Karl, flunked again: *Authorization Required. Individual Subscribers. Please log in with your APS Journals Account below.*2011-09-17
  • 0
    Ah I guess if you don't have site access from a university then you won't be able to access it.2011-09-17
  • 0
    @Karl: Do you mind specifying which line of the article you are referring too?2011-09-17
  • 0
    @eric-naslund I am referring to equation 16. Maybe I am just misinterpreting the article.2011-09-17
  • 0
    @Karl: Ok, then in that case I think my answer should be what you are looking for.2011-09-17
  • 0
    For completeness: here is a [free version](http://arxiv.org/pdf/quant-ph/0107083v1.pdf) of the article Karl refers to.2014-06-24

4 Answers 4

1

Average value is not that it usually meant as a regularized value of (say) series. One of the standard ways of regularization is to construct the analytical continuation. If it diverges at the point of interest, then the constant term of the expansion at this point is taken as the regularized value. In this case it seems to work like follows. For $s>0$ consider convergent integrals $$ \int_{-\infty }^{\infty } \frac{w^2}{\left(s^2 w^2+1\right) \left(k^2+\gamma ^2 w^2\right)} \, dw=\frac{\pi }{\gamma k s^2+\gamma ^2 s}= \frac{\pi }{\gamma ^2 s}-\frac{\pi k}{\gamma ^3}+O(s), $$ $$ \int_{-\infty }^{\infty } \frac{w^2 e^{-s w^2}}{k^2+\gamma ^2 w^2} \, dw= \frac{\sqrt{\pi }}{\gamma ^2 \sqrt{s}}-\frac{\pi k}{\gamma ^3}+O\left(\sqrt{s}\right), $$ $$ \int_{-\infty }^{\infty } \frac{w^2 e^{-s|w|}}{k^2+\gamma ^2 w^2} \, dw= \frac{2}{\gamma ^2 s}-\frac{\pi k}{\gamma ^3}+O(s). $$ Note that the constant term $-\pi k/\gamma ^3$ is the same in all three cases. So it's a candidate for the regularized value.

  • 0
    It is worth notice that this same value occurs if we just remove the $\frac{2L}{\gamma^2}$ term from the integral. That is, if we subtract off the main term. In other words $$\lim_{L\rightarrow\infty}\left(\int_{-L}^{L}\frac{\omega^{2}}{k^{2}+\gamma^{2}\omega^{2}}d\omega-\frac{2L}{\gamma^2}\right)=\frac{-\pi k}{\gamma^3}.$$2011-09-17
  • 0
    @Andrew do you perhaps have any advice where I can read up more about regularization? I come from a theoretical physics background, but not really quantum field theories (even though I make heavy use of path integrals), so I would really like to learn more. EDIT: Also, when moving the integral into the complex plane, I found also $-\pi k/\gamma^3$ + a contour integral that follows a semi circle. That contour integral seemed to diverge so I just gave up using that method.2011-09-17
  • 0
    @Karl There are different ways to regularize series/integrals. For series there is a classic book of Hardy *Divergent Series* http://www.archive.org/details/divergentseries033523mbp. As for using the analytical continuation it is well-known method. In Terry Tao blog http://terrytao.wordpress.com/2010/04/10/the-euler-maclaurin-formula-bernoulli-numbers-the-zeta-function-and-real-variable-analytic-continuation it is discussed how to obtain such results staying on real line, without use of complex analysis instruments like analytic continuation. I think the same ideas are applicable to integrals.2011-09-18
2

We can evaluate $$ \int_{-L}^{L}\frac{\omega^{2}}{k^{2}+\gamma^{2}\omega^{2}}d\omega$$ completely. Our integral is $$\frac{1}{k^{2}}\int_{-L}^{L}\frac{\omega^{2}}{1+\frac{\gamma^{2}}{k^{2}}\omega^{2}}d\omega=\frac{k}{\gamma^{3}}\int_{-\frac{\gamma L}{k}}^{\frac{\gamma L}{k}}\frac{u^{2}}{1+u^2}d\omega=\frac{2L}{\gamma^2}-\frac{k}{\gamma^3}\arctan\left(\frac{\gamma}{k}L\right).$$ The last equality follows since the anti derivative of $\frac{u^{2}}{1+u^{2}}$ is $u-\arctan(u)$.

The Average Value: From the above we can evaluate $\frac{1}{L}\lim_{L\rightarrow \infty}\frac{1}{2L} \int_{-L}^L \frac{\omega^{2}}{k^{2}+\gamma^{2}\omega^{2}}d\omega$. In particular we have $$\lim_{L\rightarrow \infty}\frac{1}{2L} \int_{-L}^L \frac{\omega^{2}}{k^{2}+\gamma^{2}\omega^{2}}d\omega =\frac{1}{\gamma^2},$$ which means that the average value of the function on the real line is $\frac{1}{\gamma^2}$. (notice I divided by $2L$ rather than $L$ because the interval is of length $2L$.)

  • 1
    I would have thought $\lim_{L\rightarrow\infty}\left(u \biggr|_{-L}^L\right) = \lim_{L\rightarrow\infty}\left(2L\right)$.2011-09-17
  • 0
    @Henry, Didier: Thanks. My mistake was the classic $L-(-L)=L-L=0$. It has been corrected.2011-09-17
  • 0
    @eric-naslund Thanks that seems to be the correct answer for my question.2011-09-17
  • 0
    @Karl: It may be worth noting for future reference that for this case the name is the "Expected Value" rather than the "Cauchy Principle Value."2011-09-17
  • 0
    @eric-naslund Noted. Physicists have a tendency to abuse mathematics... Either way, I guess this is the "best" (not that there is any way to assign a measure of good/worse...) to solve my problem. Unfortunately I still want the value of the integral and not just the average value of the function I am integrating over. I guess this will be a case of "the important behaviour is captured by the mean value" as one usually encounters in fairly gross approximations...2011-09-17
  • 0
    @Eric: A bit of nitpicking: It's Cauchy *principal* value. ("Principle" is a noun.)2011-09-17
2

Hint: Assume that $f(x)\to \ell_+$ when $x\to+\infty$ and $f(x)\to \ell_-$ when $x\to-\infty$. Then $$ \lim\limits_{L\to+\infty}\frac1{L}\int\limits_{-L}^Lf(x)\mathrm dx=\ell_++\ell_-. $$ In your case $f(x)=\dfrac{x^2}{k^2+\gamma^2x^2}$ hence $\ell_+=\ell_-=\dfrac1{\gamma^2}$ and the limit is $\dfrac2{\gamma^2}$.

As you can see the result is robust in the sense that it has nothing to do with whether we know exactly a primitive of $f$ or not, nor even with the exact form of $f$.

1

I am not sure what $1/L$ is doing in your expression for the Cauchy principal value.

In any case, your integral is as clearly infinite as $\dfrac{\int^\infty_{-\infty} d \omega}{\gamma^2} $ is as the limit of the average value is $\dfrac{1}{\gamma^2}$.

  • 0
    Sorry, My Mistake!2011-09-17
  • 1
    @Eric: I suspect the question might be asking for the Cauchy principal value of $\int^\infty_{-\infty}\mathrm d \omega \, \frac{\omega}{k^2 + \gamma^2 \omega^2}$ (note my change to the numerator) but that it is not what it says.2011-09-17
  • 0
    I want $\omega^2$ in the numerator..either way, yes the integral diverges...I just want to "fix" it in the theoretical physics sense. (renormalization/regularization). I expect great flak from proper mathematicians...2011-09-17