35
$\begingroup$

I haven't been able to find a good answer to this searching around online. There is a related old question here, but it never received much attention.

Suppose I have some physical property that I believe depends on $\int_{-\infty}^{\infty}xdx$. Ordinarily, we would say the integral is undefined, but we could also hypothetically take the Cauchy principal value and say the integral evaluates to zero.

My question is under what conditions should I let the integral remain undefined and when should it equal zero? It seems odd to me that based on different circumstances of the problem we could say that the same integral has different values. For pure math, as long as you work within a consistent framework, it doesn't really matter which you decide is true. But for physics, chemistry, etc if the integral relates to some property, it would seem there would be some definitive, empirical solution to the problem.

This question arises from some problems I have been doing with the Cauchy distribution. It seems that for certain physical examples, there is certain camp that treats the location parameter as the mean of the distribution, which is effectively saying the integral to determine the mean is equal to its principal value. I don't like this because it seems to suggest to me that all the higher order odd moments could also be argued to exist, at least in case where the distribution is symmetric about zero,

$\endgroup$
3
  • 1
    $\begingroup$ Basically it makes sense if symmetry of the integration domain is imposed by the model or by the math. For example, if the problem is due to insufficient decay at infinity, it may be that you are really not integrating from $-\infty$ to $\infty$, physically, but are rather integrating from $-M$ to $M$ and are estimating the result. If that symmetry of the integration domain comes from the physics, then PV is appropriate for your problem. $\endgroup$
    – Ian
    Commented Sep 29, 2017 at 20:44
  • $\begingroup$ This example of the Cauchy distribution IMO is easily resolved by plotting sample means. Go into Matlab and run plot(cumsum(tan(pi*rand(10000,1)-pi/2))./(1:10000)') and compare the results to plot(cumsum(randn(10000,1))./(1:10000)'). There is something fundamentally different going on there. $\endgroup$
    – Ian
    Commented Sep 29, 2017 at 20:54
  • 3
    $\begingroup$ This is a good question and I don't think it received definitive answers yet. I read on the book of Folland "Real Analysis" that considering an integral in the principal value sense is akin to subtracting infinities in quantum field theory. Indeed, he speaks of a "regularized", or "renormalized" integral. Besides, why use the term "principal", in "principal value"? It's a kind of "main order term". But I still cannot make sense of these intuitions. +1 to your question. $\endgroup$ Commented Jan 30, 2020 at 19:01

4 Answers 4

13
$\begingroup$

The Cauchy principal value is very important, especially in cases where the Lebesgue integral (which it seems you refer to as the improper integral) does not exist. The issue is that the Lebesgue integral doesn't deal too well with really big oscillations. Indeed, a measurable function $f$ is Lebesgue integrable if and only if $|f|$ is, so oscillations don't matter - only the magnitude does. This becomes problematic when dealing with functions such as $\frac{\sin x}{x}$ on $(0,\infty)$ or $\frac{1}{x}$ on $(-1,1)\setminus\{0\}$.

A purpose of the Cauchy principal value is to rectify this problem, to take into account oscillations like the Riemann integral does and give a meaningful number that represents the integral (i.e. scaled average) of the function in question. The Cauchy p.v. of $\frac{1}{x}$ is $\lim\limits_{\epsilon \to 0^+} \int_{-1}^{-\epsilon} \frac{1}{x}dx+\int_{\epsilon}^1 \frac{1}{x}dx = 0$, which coincide with our intuition for what the average value of $\frac{1}{x}$ should be.

The most prominent use of the Cauchy principal value is the Hilbert transform, in which we study $\int \frac{f(x-y)}{y}dy$, which of course, needs to be defined properly. It is critical here that we don't just say the improper integral exists, but rather get a quantitative sense of the oscillation present.

$\endgroup$
10
  • 4
    $\begingroup$ I don't think I was entirely clear with my question. I recognize the use of the Cauchy principal value. My question is more in regards to what justifies the use of it or what is its meaning. Like you say the Cauchy principal value for that integral of 1/x is 0, but when should I just go ahead and say that is the value and when should I just leave it undefined. If there aren't rigorous criteria, it would seem I can arbitrarily decide the value of the integral to suit my need. $\endgroup$
    – Tyberius
    Commented Sep 29, 2017 at 20:58
  • 1
    $\begingroup$ Of course you arbitrarily decide the value of the integral to suit your need. That's what all of mathematics is. The integral is not something that is just lurking in the air. We define it (and try to define it naturally) but of course I can't "prove" to you that the Cauchy p.v. is the correct value. What would you expect such a "proof" to entail $\endgroup$ Commented Sep 29, 2017 at 21:03
  • 3
    $\begingroup$ @mechanodroid The issue is essentially the same, it is all about the order in which positive and negative contributions to the integral are summed. $\endgroup$
    – Ian
    Commented Sep 29, 2017 at 21:15
  • 1
    $\begingroup$ @mathworker21 I'm coming at this more from the perspective of a physical scientist. Like I said in the question statement, obviously it is fine to just choose a value in pure math. But when these integrals occur in physical systems, there is an actual value for that property and so the integral does have a "correct" value in that sense. $\endgroup$
    – Tyberius
    Commented Sep 29, 2017 at 21:22
  • 2
    $\begingroup$ IMHO a good summary $\endgroup$ Commented Sep 30, 2017 at 5:33
12
$\begingroup$

A physics related example may also be useful. Consider a potential of the form $V=-Log(\left|x\right|)$. Then the force is given by $$ F=-{d V \over d x}={1\over x}. $$ The work done on a particle which moves from $-a$ to $a$ (where $a>0$) is $$ W=\int_{-a}^a{1\over x} dx. $$ From conservation of energy, we would expect $W=0$ as $V(a)=V(-a)$ and so, in that case, the Cauchy principal value gives the correct answer. However, in general, we would expect that it is not physically possible to have an infinite potential and so some new physics would come in near to $x=0$ which would change the shape of the potential in that region. But, provided that energy is still conserved we would still get $W=0$ when moving the particle from $-a$ to $a$ and so the Cauchy principle value would still have given the correct answer even without accounting directly for the new physics.

Another interesting example is when $V=-\frac{1}{x}$ which is plotted below:enter image description here

If you imagine a ball being rolled from $x<0$ it would need infinite kinetic energy to get to $x>0$. However, we expect some new physics to come in and keep $V$ finite and continuous around $V=0$. In that case the difference in potential energy between $x=-1$ and $x=1$ would simply be $V(-1)-V(1)=2$. So we would want $$ W=-\int_{-1}^1 \frac{1}{x^2}\,dx=2 $$ As explained in this answer, for this case the Cauchy principle value for the integral gives $-\infty$ while contour integration gives $2$. So in this case, contour integration is better than Cauchy's principle value from a physics perspective.

$\endgroup$
4
  • $\begingroup$ a bit of confusion here, is there a difference between the Cauchy principal value and the contour integration approach? I thought the latter was a method to find the former. $\endgroup$ Commented Apr 22, 2021 at 20:28
  • $\begingroup$ I think you are right. I have removed the last paragraph as I think it was not correct. $\endgroup$
    – Virgo
    Commented Apr 23, 2021 at 0:44
  • $\begingroup$ I think that is true, but the link you had before the edit had yet another integration method (change the $x \to x + i \epsilon$ and define the limit of $\epsilon \to 0$ to be outside of the integral) and that got the answer -2. It was really interesting to me that for the potential $V=1/x$, there is no notion of integration which gives the correct amount of change in energy!! What is the deal! I checked the $\epsilon$ method and it really gives $\int_{-1}^{1} \frac{1}{x^2} dx = -2$, whereas both the improper integral and the cauchy principal value were $\infty \neq 0$. $\endgroup$ Commented Apr 23, 2021 at 3:02
  • $\begingroup$ OK, I put it back in with some more explanation. $\endgroup$
    – Virgo
    Commented Apr 23, 2021 at 5:31
8
$\begingroup$

I like Abhimanyu Pallavi Sudhir answer: we can assign the integral $\int_{-1}^1\tfrac1x\,dx$ the value $\infty$, $-\infty$ or any other number, and we would be equally right. There is no criterion to select a result over the other. Sometimes this is presented as a flaw of the Lebesgue integral, which I think is nonsense; this integral just does not have an unambiguous meaning and it is not the fault of Lebesgue or of Riemann.

If, however, we consider the integral operator $$ T(f):=\int_{-1}^1 \frac{f(x)}{x}\, dx, $$ where $f$ is an arbitrary smooth function, then things are different. Let us develop $f$ in a Taylor series; we get, formally, $$ \int_{-1}^1 \frac{f(x)}{x}\, dx = f(0)\int_{-1}^1\frac{dx}{x} + f'(0)\int_{-1}^1\, dx + \ldots$$ The only ambiguous term is the zeroeth-order one; all the others make perfectly sense as standard integrals.

Conclusion. We need a renormalization procedure that removes the zeroeth-order term. This procedure amounts to exploit symmetry, prescribing that $\int_{-1}^1 \frac{1}{x}\, dx$ should be interpreted as $$ \lim_{\epsilon \to 0}\int_{\epsilon<\lvert x \rvert <1} \frac{dx}{x}=0,$$ that is, considering integrals in the principal value sense.

This is the starting point of the Calderón-Zygmund theory of singular integrals. In general, given $$ K(x)=\frac{\Omega(\tfrac{x}{|x|})}{|x|^n}, $$ where $\Omega$ is a function on the unit sphere $\mathbb S^{n-1}$ with zero integral, that is, $$ \int_{\mathbb S^{n-1}}\Omega\, dS =0, $$ we can define an integral operator $$ Sf(x):=\int_{\mathbb R^n} K(x-y)f(y)\, dy, $$ where the integral is interpreted in the principal value sense, exactly as before. Some conditions must be put on $K$ to ensure that this operator is continuous in appropriate function spaces. But most importantly, this definition based on the principal value is not crazy, or arbitrary, as it might seem at first sight. It is the only possible choice that removes the singular term from the integral, while leaving all the other terms untouched.


A final remark. I think that "principal value" is a gallicism, which obscures the true meaning. "Renormalization" would be better, as Gerry Folland suggests in his book "Real Analysis".

$\endgroup$
1
  • 1
    $\begingroup$ This is the most reasonable answer, actually giving a reason why one could consider such integrals. $\endgroup$
    – shuhalo
    Commented Dec 10, 2020 at 4:03
3
$\begingroup$

It isn't the "correct value" for the integral any more than the principal root is the correct value for a root or the principal logarithm is the correct value for a logarithm, or setting $C=0$ gives you the correct value of the antiderivative. There are plenty of other values the integral can take, depending on how you take the limit. See my answer to Why can't $\int_{-1}^1\frac{dx}x$ be evaluated?

$\endgroup$
1

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .