Green's function approach using Lebesgue integration

In summary, it is difficult to use the Green's function approach rigorously when considering the Dirac Delta function. While the integral solution for Laplace's Equation is well-defined as both a Riemann-Darboux and Lebesgue integral, it is not possible to exchange the limiting operations in the Lebesgue integral. One possible method to solve this issue is to use a variable change and apply calculus tricks with the divergence theorem and integration by parts. However, the use of Dirac measure and delta distribution may only be useful for defining properties and stating results, rather than actually proving them.
  • #1
bdforbes
152
0
I can't figure out how to use the Green's function approach rigorously, i.e., taking into account the fact that the Dirac Delta function is not a function on the reals.

Suppose we have Laplace's Equation:

[tex]\nabla^2 \phi(\vec{x})=f(\vec{x})[/tex]

The solution, for "well-behaved" [itex]f(\vec{x})[/itex] is

[tex]\phi(\vec{x})=\frac{-1}{4\pi}\int \frac{f(\vec{x}')}{\left|\vec{x}-\vec{x}'\right|}d^3\vec{x}'[/tex]

It is my understanding that this integral is well-defined as both a Riemann-Darboux and Lebesgue integral. If we treat it as a Lebesgue integral, I believe the limiting operations can be exchanged, i.e., we can apply the Laplacian to the integrand:

[tex]\nabla^2 \phi(\vec{x})=\frac{-1}{4\pi}\int\nabla^2\left(\frac{1}{\left|\vec{x}-\vec{x}'\right|}\right)f(\vec{x}')d^3\vec{x}'[/tex]

But now it looks like this Lebesgue integral is NOT well-defined! How do we deal with [itex]\nabla^2(1/|x-x'|)[/itex] at the singular point?
If we naively apply the divergence theorem, we can arrive at the desired result, but that is not good enough for me.

How can we do this rigorously? Is there a way to use the Dirac measure, or the Dirac Delta as a linear functional?
 
Physics news on Phys.org
  • #2
bdforbes said:
I can't figure out how to use the Green's function approach rigorously, i.e., taking into account the fact that the Dirac Delta function is not a function on the reals.

Suppose we have Laplace's Equation:

[tex]\nabla^2 \phi(\vec{x})=f(\vec{x})[/tex]

The solution, for "well-behaved" [itex]f(\vec{x})[/itex] is

[tex]\phi(\vec{x})=\frac{-1}{4\pi}\int \frac{f(\vec{x}')}{\left|\vec{x}-\vec{x}'\right|}d^3\vec{x}'[/tex]

It is my understanding that this integral is well-defined as both a Riemann-Darboux and Lebesgue integral.

According to my experience, we are usually interested in such functions [itex]f[/itex], that it doesn't matter what definition of the integral we are using. They all agree on the relevant functions [itex]f[/itex].

If we treat it as a Lebesgue integral, I believe the limiting operations can be exchanged, i.e., we can apply the Laplacian to the integrand:

That's a mistake. If you actually tried to find some dominating function so that you could use the Lebesgue's dominated convergence, you would notice that you cannot find suitable dominating function. Take a look at the simpler example, [itex]\partial_x \theta(x-x') = \delta(x-x')[/itex]. If you try to do this

[tex]
\partial_x \int dx'\; \theta(x-x') f(x') = \int dx'\; \partial_x \theta(x-x') f(x'),
[/tex]

you would find yourself trying to find such integrable function [itex]h[/itex] that

[tex]
\Big|\frac{\theta(x+\Delta x - x') - \theta(x - x')}{\Delta x}\Big| = \frac{1}{|\Delta x|} \chi_{[x,x+\Delta x]}(x') \leq h(x')
[/tex]

for all [itex]\Delta x[/itex], and that's not possible.

[tex]\nabla^2 \phi(\vec{x})=\frac{-1}{4\pi}\int\nabla^2\left(\frac{1}{\left|\vec{x}-\vec{x}'\right|}\right)f(\vec{x}')d^3\vec{x}'[/tex]

But now it looks like this Lebesgue integral is NOT well-defined! How do we deal with [itex]\nabla^2(1/|x-x'|)[/itex] at the singular point?

Actually the integrals are well defined, because in Lebesgue integration you can always ignore values of functions at single points (or other sets of zero measure). The integral in the right side is simply zero. On the other hand the integral on the left side is usually not zero. So both sides of the equation are well defined, but the equation it self is not correct.

If we naively apply the divergence theorem, we can arrive at the desired result, but that is not good enough for me.

How can we do this rigorously?

One way is do a variable change like this

[tex]
\nabla^2_x \int \frac{f(x')}{|x-x'|} d^3x' = \nabla^2_x \int \frac{f(x-u)}{|u|} d^3u
[/tex]

and then use niceness properties of [itex]f[/itex] for commutation of derivation and integration. After commutation some calculus trickery with divergence theorem and integration by parts will be needed.

Is there a way to use the Dirac measure, or the Dirac Delta as a linear functional?

According to my experience, Dirac measure and delta distribution are useful for defining some properties by definition, or for stating results once they have been proven, but not useful for actually proving anything.
 
  • #3
Thanks for those insights, that's very helpful. I'll work through the variable change method and see what I get.
 
  • #4
jostpuur said:
One way is do a variable change like this

[tex]
\nabla^2_x \int \frac{f(x')}{|x-x'|} d^3x' = \nabla^2_x \int \frac{f(x-u)}{|u|} d^3u
[/tex]

and then use niceness properties of [itex]f[/itex] for commutation of derivation and integration. After commutation some calculus trickery with divergence theorem and integration by parts will be needed.

I've been trying to do this for the simpler one dimensional case, but I don't get the desired result. Here is my working:

[tex]\frac{\partial^2}{\partial x^2}\phi(x)=f(x)[/tex]

[tex]\phi(x)=\int_a^b \frac{1}{2}|x-x'|f(x')dx'[/tex] with a<x<b

Let u=x-x'

[tex]\frac{\partial^2}{\partial x^2}\phi(x)=-\int_{x-a}^{x-b}\frac{1}{2}|u|\frac{\partial^2}{\partial x^2}f(x-u)du

=\left[\frac{1}{2}|u|\frac{\partial}{\partial x}f(x-u)\right]^{x-a}_{x-b}

=\frac{1}{2}|x-a|f'(a)-\frac{1}{2}|x-b|f'(b)[/tex]

The final result almost looks right, but it would only work for a and b very close to x, wouldn't it?
 
  • #5
Those are calculation mistakes, made in too quick calculation. Notice that you need pay attention how to switch operators [tex]\frac{\partial}{\partial x}[/tex] and [tex]\frac{\partial}{\partial u}[/tex]. (It could be you would notice this soon anyway...)

Not any [itex]a,b,f[/itex] are going to be fine. [itex]f[/itex] ([itex]f'[/itex] too) will have to be sufficiently zero somewhere for the integration by parts to work.
 
  • #6
I'm not sure what you mean. If we view u as depending on x, why are we still allowed to commute the integration and differentiation? I thought the point was switch the x dependence to the function f because we can assume it is a nice function.
 
  • #7
It could be a good idea to set [itex]a=-\infty[/itex] and [itex]b=\infty[/itex] in the beginning. It will save from some trouble when the order of derivation and integration are supposed to be changed in

[tex]
\frac{d^2}{dx^2} \int\limits_{x-b}^{x-a} \cdots du
[/tex]

The boundaries will give some extra terms if [itex]|a|,|b|<\infty[/itex]. Actually I've never tried to carry out calculations like the in full rigor. It is not awfully difficult to get right answers, though.
 
Last edited:
  • #8
I'll now set [itex]a=-\infty[/itex] and [itex]b=\infty[/itex], and assume that [itex]f[/itex] has sufficient properties so that [itex]f(\pm\infty)[/itex] and [itex]f'(\pm\infty)[/itex] will bring all else to zero in integration by parts, and also so that [itex]\frac{d^2}{dx^2}[/itex] and [itex]\int du[/itex] can be commuted after the change of variable [itex]u=x-x'[/itex].

[tex]
\frac{d^2}{dx^2} \int\limits_{-\infty}^{\infty} \frac{1}{2}|x-x'| f(x')dx' = \int\limits_{-\infty}^{\infty} \frac{1}{2}|u| \frac{d^2}{dx^2} f(x-u) du = -\frac{1}{2}\int\limits_{-\infty}^0 u \frac{d^2}{du^2} f(x-u) du + \frac{1}{2}\int\limits_0^{\infty} u \frac{d^2}{du^2} f(x-u)du = \cdots
[/tex]

Now substitute

[tex]
u \frac{d^2}{du^2} f(x-u) = \frac{d}{du}\Big( u \frac{d}{du} f(x-u)\Big) - \frac{d}{du} f(x-u).
[/tex]

[tex]
\cdots = \infty f'(\infty) - \infty f'(-\infty) + \underbrace{\frac{1}{2}\int\limits_{-\infty}^0 \frac{d}{du}f(x-u) du}_{=\frac{1}{2}f(x) - \frac{1}{2}f(-\infty)} - \underbrace{\frac{1}{2}\int\limits_0^{\infty} \frac{d}{du} f(x-u) du}_{=-\frac{1}{2}f(\infty) + \frac{1}{2}f(x)} = f(x)
[/tex]
 
Last edited:
  • #9
Thanks, it does look like this can be done rigorously. One question, why do we turn [itex]d^2/dx^2[/itex] into [itex]d^2/du^2[/itex]?

I believe my result
[tex]
\frac{1}{2}|x-a|f'(a)-\frac{1}{2}|x-b|f'(b)
[/tex]

is also good, if we only consider a very small interval around the singularity. The other contributions would go to zero anyway since if we avoid the singularity, we can take the double derivative of |x-x'| which is zero. My answer approaches f(x) in the limit a->x<-b.
 
  • #10
bdforbes said:
One question, why do we turn [itex]d^2/dx^2[/itex] into [itex]d^2/du^2[/itex]?

For the purpose of using the fundamental theorem of calculus, and integration by parts. We are integrating over [itex]u[/itex], so we want derivatives with respect to [itex]u[/itex] inside the integral too.
 
  • #11
But why is it valid to do it? It's not immediately obvious to me.
 
  • #12
[tex]
\frac{d}{dx} f(x-u) = f'(x-u)
[/tex]

[tex]
\frac{d}{du} f(x-u) = -f'(x-u)
[/tex]

Then take second derivatives the same way.
 
  • #13
Thanks, that makes sense.
 
  • #14
What if I now wanted to use a generalised function approach? Ie, write the Green's function as a limit of ordinary functions, which would enable me to commute the differentiation and integration. If I chose the sequence cleverly, perhaps the second derivative would give a sequence equivalent to the Dirac delta, and I would thus arrive at the result. It feels like the usual physicist approach to Green's functions, whereby one flails around wildly until a reasonable result if found, is basically shorthand for using generalised functions. The method you have given above appears to be qualitatively different.
 
  • #15
bdforbes said:
What if I now wanted to use a generalised function approach? Ie, write the Green's function as a limit of ordinary functions, which would enable me to commute the differentiation and integration. If I chose the sequence cleverly, perhaps the second derivative would give a sequence equivalent to the Dirac delta, and I would thus arrive at the result. It feels like the usual physicist approach to Green's functions, whereby one flails around wildly until a reasonable result if found, is basically shorthand for using generalised functions. The method you have given above appears to be qualitatively different.

Good thing that you asked. For some reason I didn't bother trying to mention the other way I already knew about. Notice that I was careful to say

jostpuur said:
One way is to do a variable change like this...

in my original response. I didn't claim it would be "the" way.

I learned the change of variable trick [itex]u=x-x'[/itex] from the https://www.amazon.com/dp/0821807722/?tag=pfamazon01-20 some years ago. Evans wasn't particularly speaking about Green's functions, but instead only stated that solutions to some PDEs could be written as some integral expressions. When looking at the proofs, I recognized a solution also to the Green's function problem that had been bothering me already earlier.

I was reading the https://www.amazon.com/dp/047130932X/?tag=pfamazon01-20 this spring, and to my positive surprise I noticed that Jackson too comments this same problem. On page 35 (of 3th edition) he shows a following calculation

[tex]
\nabla_x^2 \int d^3x'\; \frac{\rho(x')}{\sqrt{|x-x'|^2}} = \lim_{a\to 0} \int d^3x'\; \nabla^2_x \frac{\rho(x')}{\sqrt{|x-x'|^2 + a^2}} = \lim_{a\to 0} \int d^3x'\; \nabla_x\cdot\Big(-\frac{(x-x')\rho(x')}{(|x-x'|^2 + a^2)^{3/2}}\Big)
[/tex]
[tex]
= -\lim_{a\to 0} \int d^3x'\; \frac{3a^2 \rho(x')}{(|x-x'|^2 + a^2)^{5/2}} = -4\pi \rho(x)
[/tex]

If you are interested in rigor, this is not necessarily easier than the trick I got from Evans' book. You would first need to justify the commutation of [itex]\int d^3x'[/itex] and [itex]\lim_{a\to 0}[/itex], then the commutation of [itex]\nabla^2_x[/itex] and [itex]\lim_{a\to 0}[/itex], and then the commutation of [itex]\nabla^2_x[/itex] and [itex]\int d^3x'[/itex]. So there's lot to do. I have never tried to figure out what arguments one should use to justify these (obviously Jackson doesn't speak about justifying them either, because his is a physics book), but I believe that they can be justified, because these steps do not lead to paradoxes. Like for example direct commutation of [itex]\nabla^2_x[/itex] and [itex]\int d^3x'[/itex] does lead to a paradox.

The last steps uses a delta function identity

[tex]
\frac{3a^2}{(|x-x'|^2 + a^2)^{5/2}} \underset{a\to 0}{\to} 4\pi \delta^3(x-x')
[/tex]

It can be showed with a variable change [itex]x'=x+au[/itex].

[tex]
\int d^3x'\; \frac{3a^2\rho(x')}{(|x-x'|^2 + a^2)^{5/2}} = \int d^3u\; a^3 \frac{3a^2\rho(x+au)}{(|au|^2 + a^2)^{5/2}} = 3 \int d^3u\; \frac{\rho(x+au)}{(|u|^2 + 1)^{5/2}}
[/tex]
[tex]
\underset{a\to 0}{\to} 3\int d^3u\; \frac{\rho(x)}{(|u|^2 + 1)^{5/2}} = 12\pi \rho(x) \underbrace{\int\limits_0^{\infty} \frac{r^2 dr}{(r^2 + 1)^{5/2}}}_{\cdots = 1/3} = 4\pi \rho(x)
[/tex]
 
Last edited by a moderator:
  • #16
bdforbes said:
How can we do this rigorously? Is there a way to use the Dirac measure, or the Dirac Delta as a linear functional?

You can. First, you need the identity

[tex]
\nabla_x^2\left(\frac{1}{|x-x^\prime|}\right)=-4\pi\delta(x-x^\prime).
\eqno{(1)}
[/tex]

Both sides of this equality are to be regarded as distributions. That is, they are linear functionals mapping the smooth real valued functions of compact support [itex]u\colon\mathbb{R}^3\to\mathbb{R}[/itex] to the real numbers, with the mapping being written as the integral against u.

[tex]
\int u(x)\nabla_x^2\left(\frac{1}{|x-x^\prime|}\right)\,dx=4\pi\int u(x)\delta(x-x^\prime)\,dx.
[/tex]

The right hand side is simply [itex]u(x^\prime)[/itex] by definition of the Dirac delta. The derivative 'in the sense of distributions' on the left hand side is defined via integration by parts

[tex]
\int \left(\nabla_x^2u(x)\right)\frac{1}{|x-x^\prime|}\,dx=4\pi u(x^\prime).\eqno{(2)}
[/tex]

Equation (2) is then an equivalent statement of (1). Can prove it by using the divergence theorem.

You want to prove

[tex]
\nabla^2_{x^\prime}\,\phi(x^\prime)=f(x^\prime)
[/tex]

which, in distributional form, is equivalent to

[tex]
\int \left(\nabla^2_{x^\prime}\,u(x^\prime)\right)\phi(x^\prime)\,dx^\prime=\int u(x^\prime)f(x^\prime)\,dx^\prime.
[/tex]

To prove this, substitute in the expression for [itex]\phi[/itex], commute the order of integration, and apply (1) or equivalently (2).
 
Last edited:
  • #17
jostpuur said:
I learned the change of variable trick [itex]u=x-x'[/itex] from the https://www.amazon.com/dp/0821807722/?tag=pfamazon01-20 some years ago.

I will definitely take a look at that!

jostpuur said:

I have been reading it too by coincidence.

jostpuur said:
If you are interested in rigor, this is not necessarily easier than the trick I got from Evans' book. You would first need to justify the commutation of [itex]\int d^3x'[/itex] and [itex]\lim_{a\to 0}[/itex], then the commutation of [itex]\nabla^2_x[/itex] and [itex]\lim_{a\to 0}[/itex], and then the commutation of [itex]\nabla^2_x[/itex] and [itex]\int d^3x'[/itex]. So there's lot to do. I have never tried to figure out what arguments one should use to justify these (obviously Jackson doesn't speak about justifying them either, because his is a physics book), but I believe that they can be justified, because these steps do not lead to paradoxes. Like for example direct commutation of [itex]\nabla^2_x[/itex] and [itex]\int d^3x'[/itex] does lead to a paradox.
Perhaps we need to appeal to absolute continuity or uniform convergence?
 
Last edited by a moderator:
  • #18
gel said:
You can. First, you need the identity

[tex]
\nabla_x^2\left(\frac{1}{|x-x^\prime|}\right)=\delta(x-x^\prime).
\eqno{(1)}
[/tex]

Both sides of this equality are to be regarded as distributions. That is, they are linear functionals mapping the smooth real valued functions of compact support [itex]u\colon\mathbb{R}^3\to\mathbb{R}[/itex] to the real numbers, with the mapping being written as the integral against u.

I can believe this, it is the result obtained formally by ignoring the singularity and applying the divergence theorem, which must be equivalent to treating the objects as distributions.

gel said:
[tex]
\int u(x)\nabla_x^2\left(\frac{1}{|x-x^\prime|}\right)\,dx=\int u(x)\delta(x-x^\prime)\,dx.
[/tex]

The right hand side is simply [itex]u(x^\prime)[/itex] by definition of the Dirac delta. The derivative 'in the sense of distributions' on the left hand side is defined via integration by parts

[tex]
\int \left(\nabla_x^2u(x)\right)\frac{1}{|x-x^\prime|}\,dx=u(x^\prime).\eqno{(2)}
[/tex]

Equation (2) is then an equivalent statement of (1). Can prove it by using the divergence theorem.

Do you mean that we can prove statement (2) is true, or that it is equivalent to (1)? It seems to me that the latter has already been proven above.
Also, could you demonstrate how to use the divergence theorem rigorously here? Every time I see the divergence theorem used in this context, the author ignores the second bounding surface which arises from the integral being improper. This second bounding surface always cancels out the contribution from the first. Their mistake was really to commute differentiation and integration with a singular integrand.

gel said:
You want to prove

[tex]
\nabla^2_{x^\prime}\,\phi(x^\prime)=f(x^\prime)
[/tex]

which, in distributional form, is equivalent to

[tex]
\int \left(\nabla^2_{x^\prime}\,u(x^\prime)\right)\phi(x^\prime)\,dx^\prime=\int u(x^\prime)f(x^\prime)\,dx^\prime.
[/tex]

To prove this, substitute in the expression for [itex]\phi[/itex], commute the order of integration, and apply (1) or equivalently (2).

This seems very reasonable. Are there any conditions required for commuting the order of integration? Also, did the Laplacian find itself acting on u(x') via integration by parts, as above? I assume compact support is essential for this to work.
 
  • #19
bdforbes said:
Equation (2) is then an equivalent statement of (1). Can prove it by using the divergence theorem.
Do you mean that we can prove statement (2) is true, or that it is equivalent to (1)?

I hope he meant that the equation (2) can be proven using divergence theorem (and possibly some other calculus stuff). You should not attempt to prove that equations (1) and (2) are equivalent. gel's equation (2) is the definition of the equation (1). The equation (1) should be considered to be notation for the equation (2).
 
  • #20
jostpuur said:
I hope he meant that the equation (2) can be proven using divergence theorem (and possibly some other calculus stuff). You should not attempt to prove that equations (1) and (2) are equivalent. gel's equation (2) is the definition of the equation (1). The equation (1) should be considered to be notation for the equation (2).

That makes sense. I guess using the divergence theorem here would be similar to the method you showed earlier, right jostpuur?
 
  • #21
jostpuur said:
I hope he meant that the equation (2) can be proven using divergence theorem (and possibly some other calculus stuff). You should not attempt to prove that equations (1) and (2) are equivalent. gel's equation (2) is the definition of the equation (1). The equation (1) should be considered to be notation for the equation (2).

Yes, (2) is just the definition of (1). Can break it into two steps. First, we can prove
[tex]
\nabla\left(\frac{1}{|x|}\right)=-\frac{\hat x}{|x|^2}. \eqno{(3a)}
[/tex]
This is easily calculated for x != 0, but we need to show that it holds everywhere in the sense of distributions. In the sense of distributions, (3a) is an equivalent statement to
[tex]
\int \left(\nabla u(x)\right)\frac{1}{|x|}\,dx=\int u(x)\frac{\hat x}{|x|^2}\,dx\eqno{(3b)}
[/tex]
for all smooth u with compact support.
Choose r>0 and let Br be a ball of radius r, with Bcr being its complement.
[tex]
\begin{align*}
\int \left(\nabla u(x)\right)\frac{1}{|x|}\,dx &=\lim_{r\rightarrow 0}\int_{B_r^c}\left(\nabla u(x)\right)\frac{1}{|x|}\,dx\\
&= \lim_{r\rightarrow 0}\int_{B_r^c}\nabla\left(u(x)/|x|\right)\,dx - \lim_{r\to 0}\int_{B_r^c}u(x)\nabla\left(\frac{1}{|x|}\right)\,dx\\
&= -\lim_{r\to 0}\int_{S_r} \hat x u(x)/|x|\,d\sigma(x) + \lim_{r\to 0}\int_{B_r^c} u(x)\frac{\hat x}{|x|^2}\,dx
\end{align*}
[/tex]

Here, the divergence theorem has been applied to the first integral on the RHS to write it as the surface integral [itex]d\sigma(x)[/itex] over the sphere Sr of radius r.
Note that there is no 'second surface' because u is chosen to have compact support. The surface integral is bounded by [itex]\Vert u\Vert 4\pi r[/itex] which vanishes as r->0, so we get (3b).

The second part is to prove
[tex]
\nabla \cdot\left(\frac{\hat x}{|x|^2}\right)=4\pi\delta(x).\eqno{(4a)}
[/tex]
Again, this is differentiation in the sense of distributions and is equivalent to
[tex]
\int \nabla(u(x))\cdot\left(\frac{\hat x}{|x|^2}\right)\,dx=-4\pi u(0).\eqno{(4b)}
[/tex]
Can prove this in a similar way as above.
[tex]
\begin{align*}
\int \nabla(u(x))\cdot\left(\frac{\hat x}{|x|^2}\right)\,dx &=\lim_{r\to 0}\int_{B_r^c} \nabla(u(x))\cdot\left(\frac{\hat x}{|x|^2}\right)\,dx\\
&=\lim_{r\to 0}\int_{B_r^c} \nabla\cdot \left(u(x)\frac{\hat x}{|x|^2}\right)\,dx-\lim_{r\to 0}\int_{B_r^c}u(x)\nabla\cdot\left(\frac{\hat x}{|x|^2}\right)\,dx
\end{align*}
[/tex]
The second integral on the RHS is 0, since [itex]\nabla\cdot(\hat x/|x|^2)=0[/itex] for x!=0. Use the divergence theorem for the first integral (again, no second surface as u has compact support).
[tex]
\begin{align*}
\int \nabla(u(x))\cdot\left(\frac{\hat x}{|x|^2}\right)\,dx
&=-\lim_{r\to 0}\int_{S_r}\hat x\cdot\left(u(x)\frac{\hat x}{|x|^2}\right)\,d\sigma(x)\\
&=-\lim_{r\to 0}\int_{S_r}u(x)\frac{1}{|x|^2}\,d\sigma(x)
\end{align*}
[/tex]
As |x|->0 then u(x)->u(0), so this integral converges to [itex]-u(0)\int_{S_r}(1/r^2)\,d\sigma(x)=-4\pi u(0)[/itex], giving (4b).

Then, (3a)+(4a) gives

[tex]
\nabla^2\left(\frac{1}{|x|}\right) = -\nabla\cdot\left(\frac{\hat x}{|x|^2}\right)=-4\pi\delta(x).
[/tex]
Finally, its just a translation to replace x by x-x'.
 
Last edited:
  • #22
bdforbes said:
This seems very reasonable. Are there any conditions required for commuting the order of integration? Also, did the Laplacian find itself acting on u(x') via integration by parts, as above? I assume compact support is essential for this to work.

More explicitly
[tex]
\begin{align*}
\int \left(\nabla_{x^\prime}^2\,u(x^\prime)\right)\phi(x^\prime)\,dx^\prime
&=-\frac{1}{4\pi}\int\int \left(\nabla_{x^\prime}^2\,u(x^\prime)\right)f(x)\frac{1}{|x-x^\prime|}\,dx\,dx^\prime\\
&=-\frac{1}{4\pi}\int\int \left(\nabla^2_{x^\prime}\,u(x^\prime)\right)f(x)\frac{1}{|x-x^\prime|}\,dx^\prime\,dx\textrm{\ \ (commute order of integration)}\\
&=-\frac{1}{4\pi}\int f(x)\int \left(\nabla_{x^\prime}^2\,u(x^\prime)\right)\frac{1}{|x-x^\prime|}\,dx^\prime\,dx\\
&=-\frac{1}{4\pi}\int f(x) (-4\pi u(x))\,dx = \int u(x)f(x)\,dx
\end{align*}
[/tex]

Commuting the order of integration works, because the integrand is integrable (Fubini's theorem). Its easy to show this if u and f have compact support. If f isn't compact support then you'll have to enforce some other boundary condition to ensure integrability.
The result of this is to prove that [itex]\nabla^2\phi=f[/itex] in the sense of distributions.
 
Last edited:
  • #23
That's one way, but it doesn't look optimal. It is easier to begin like this

jostpuur said:
[tex]
\nabla^2_x \int \frac{f(x')}{|x-x'|} d^3x' = \nabla^2_x \int \frac{f(x-u)}{|u|} d^3u
[/tex]

[tex]
\cdots = \int \nabla^2_x \frac{f(x-u)}{|u|} d^3u
[/tex]

and then carry out the calculation to the end mostly like was shown in the post #21.

Also recall the beginning of the problem. We want to solve a PDE

bdforbes said:
Suppose we have Laplace's Equation:

[tex]\nabla^2 \phi(\vec{x})=f(\vec{x})[/tex]

and here [itex]f[/itex] is usually a function, not a distribution, and [itex]\nabla^2[/itex] is a derivative operator defined by limit. The claim is that if we set

The solution, for "well-behaved" [itex]f(\vec{x})[/itex] is

[tex]\phi(\vec{x})=\frac{-1}{4\pi}\int \frac{f(\vec{x}')}{\left|\vec{x}-\vec{x}'\right|}d^3\vec{x}'[/tex]

where the integral is a Riemann integral, then this [itex]\phi[/itex] will satisfy the original PDE. If the proof goes through a step that looks like this

gel said:
[tex]
\int \left(\nabla^2_{x^\prime}\,u(x^\prime)\right)\phi(x^\prime)\,dx^\prime=\int u(x^\prime)f(x^\prime)\,dx^\prime.
[/tex]

it creates an appearance that the nabla-operator in the original PDE problem would have been something else than an ordinary derivative operator defined by limit.
 
  • #24
Yeah, my proof shows that [itex]\nabla^2\phi = f[/itex] in the sense of distributions. The 'in the sense of distributions' can be dropped if you know that phi is indeed twice differentiable. Haven't looked through your proof in enough detail to tell if it proves twice differentiability or assumes it (does it?).
Doing it in terms of distributions is a bit more optimal than it looks as I wrote it out, because I was careful to write it out in 'dual' form (integrated vs u). Really, all that is happening is that you are applying
[tex]
\begin{align*}
\nabla_{x}^2\phi(x)&=-\frac{1}{4\pi}\int f(x^\prime)\nabla_x^2\left(\frac{1}{|x-x^\prime|}\right)\,dx^\prime\\
&= -\frac{1}{4\pi}\int f(x^\prime)(-4\pi)\delta(x-x^\prime)\,dx^\prime\\
&= f(x).
\end{align*}
[/tex]
Then, I was integrating vs a smooth function of compact support u(x) in order to make both steps rigorous.Edit: Also, the identity [itex]\nabla^2(1/|x|)=-4\pi\delta(x)[/itex] only makes sense in terms of distributions. So, if you want to use that, then you have introduced distributions from the start.
 
Last edited:
  • #25
I did not follow the whole discussion, but it seems to me that you approach the problem from the wrong direction. I recommend my book 'Introduction to Boundary Elements', Springer-Verlag, 1989 p. 108 the section with the title: The influence function for the deflection u(x) (that is the potential).

There you see how by considering Green's second identity on a punctured domain, you spare the source point x (where the Dirac delta is located), the influence function for the potential is easily derived.

The same simple approach works in 3-D. No need to discuss Lebesgue or Riemann etc.

Friedel Hartmann
 
  • #26
I'm not about to buy the book mentioned in Friedel's post, but from Wikipedia, Green's 2nd identity is

[tex] \int_U \left( \psi \nabla^2 \varphi - \varphi \nabla^2 \psi\right)\, dV = \oint_{\partial U} \left( \psi {\partial \varphi \over \partial n} - \varphi {\partial \psi \over \partial n}\right)\, dS. [/tex]

which you can use with U=Brc to get something like I had using divergence theorem above, if you put my 2 steps together
 
  • #27
Thinking about this a bit more. My approach showed that [itex]\nabla^2\phi = f[/itex] in the sense of distributions. To be rigorous, you would also have to show that [itex]\phi[/itex] is twice differentiable. Maybe this isn't true if f is only assumed to be continuous (and compact support, for simplicity).
If you assume that f is continuously differentiable, then you can show

[tex]
\nabla_i\nabla_j\phi(x) = \frac{-1}{4\pi}\int \left(\nabla_i f(x+y)\right)\frac{\hat y}{|y|^2}\,dy.
[/tex]

This can be proven in the same way as my workings out above (i.e. prove 'in the sense of distributions' by integrating vs a smooth function). The right hand side is continuous if f is continuously differentiable, so [itex]\phi[/itex] is twice continuously differentiable.

If f isn't continuously differentiable, then maybe [itex]\phi[/itex] isn't twice continuously differentiable?
 
Last edited:
  • #28
Consider the example [itex]f(x) = 1_{\{|x|<1\}}\sqrt{(1-|x|^3)}[/itex]. Here f is continuous, but not differentiable. As it is spherically symmetric, we can solve by concentrating the mass at the origin (http://en.wikipedia.org/wiki/Shell_theorem" [Broken])
[tex]
\nabla \phi(x) = \frac{M_{|x|}\hat x}{4\pi|x|^2}.
[/tex]

Here, Mr is the total mass inside radius r. For r <= 1,
[tex]
M_r = 4\pi\int_0^r\sqrt{1-s^3}s^2\,ds=\frac{8\pi}{9}-\frac{8\pi}{9}(1-r^3)^{\frac{3}{2}}
[/tex]
and, for r>=1,
[tex]
M_r=M_1=\frac{8\pi}{9}
[/tex]
This is continuous, but not differentiable at |x|=1, so strictly speaking [itex]\nabla^2\phi=f[/itex] fails except in the sense of distributions.

Edit: Mr is continuously differentiable at r=1. interesting.
 
Last edited by a moderator:
  • #29
gel said:
If f isn't continuously differentiable, then maybe [itex]\phi[/itex] isn't twice continuously differentiable?

That must be true, since the "second derivative" of phi IS f. What would be the implication of this situation? Could it be that the closed integral form for phi I quoted in the OP would not work then?

PS: I still have to work through your posts in detail, but it all looks good.

EDIT: You answered me while I was drafting my response.
 
  • #30
gel said:
This is continuous, but not differentiable at |x|=1, so strictly speaking [itex]\nabla^2\phi=f[/itex] fails except in the sense of distributions.

Could this occur in a physical situation? We often get point sources and surface discontinuities in real applications. How should we interpret the fact that the differential equation only holds in the sense of distributions?
 
  • #31
bdforbes said:
That must be true, since the "second derivative" of phi IS f. What would be the implication of this situation? Could it be that the closed integral form for phi I quoted in the OP would not work then?

First, I don't think it is obvious without looking at an example. If f wasn't continuous then clearly phi couldn't be twice continuously differentiable.
If f is continuous, but not continuously differentiable, then it is not clear to me except by looking at an example. My example above shows that in fact phi might not be twice differentiable. In this case [itex]\nabla^2\phi=f[/itex] only makes sense in a distributional sense, or restricted to the regions where f is differentiable.
 
  • #32
bdforbes said:
Could this occur in a physical situation? We often get point sources and surface discontinuities in real applications. How should we interpret the fact that the differential equation only holds in the sense of distributions?

A point source is tricky. Then you have the density [itex]f(x)=\delta(x)[/itex] being a distribution.
If f is continuously differentiable everywhere except for on some surface, where it is continuous, then you can say that [itex]\nabla^2\phi=f[/itex] away from this surface, and that [itex]\phi[/itex] is once continuously differentiable (across the surface). (Maybe even twice diff - see my Edit above). If f is discontinuous at the surface, then phi will be continuous but not continuously differentiable here (Edit: maybe it is continuously differentiable. Not sure).

Anyway, way past my bedtime, so that's enough from me for today.
 
Last edited:
  • #33
Here is my working for Jackson's approach to the one-dimensional Laplace equation.

[tex]\frac{\partial^2}{\partial x^2}\phi(x)=f(x)[/tex]

[tex]\phi(x)=\int_{-\infty}^\infty \frac{1}{2}|x-x'|f(x')dx'[/tex]

[tex]\phi_a(x)=\int_{-\infty}^\infty \frac{1}{2}\sqrt{(x-x')^2+a^2}f(x')dx'[/tex]

[tex]
\begin{align*}
\frac{\partial^2}{\partial x^2}\phi_a(x)&=\int_{-\infty}^\infty \frac{1}{2}\frac{a^2}{\left[(x-x')^2+a^2\right]^{3/2}}f(x')dx'\\
&=\int_{-\infty}^\infty \frac{1}{2}\frac{a^2}{(r^2+a^2)^{3/2}}f(x+r)dr\\
&=\int_{|r|>R} + \int_{-R}^{R} \frac{1}{2}\frac{a^2}{(r^2+a^2)^{3/2}}f(x+r)dr\\
&=\int_{-R}^{R} \frac{1}{2}\frac{a^2}{(r^2+a^2)^{3/2}}\left[\Sigma \frac{f^{(n)}(x)}{n!}r^n\right]dr\\
&=\left[\frac{r}{2\sqrt{r^2+a^2}}f(x)+\cdots\right]_{-R}^R\\
&=\frac{R}{\sqrt{R^2+a^2}}f(x)+\cdots\\
&\rightarrow f(x) \mbox{ as }a\rightarrow 0
\end{align*}
[/tex]

The contribution from |r|>R goes like [itex]a^2[/itex], so it vanishes as a->0.
The other terms in the series are [itex]O(a^2)[/itex], so they vanish.
I am slightly uncomfortable with assuming the function has a Taylor series expansion. Is it general enough?

Here is the key part I guess:

[tex]\phi(x)=\lim_{a\rightarrow 0}\phi_a(x)[/tex]

[tex]
\frac{\partial^2}{\partial x^2}\phi(x)=\frac{\partial^2}{\partial x^2}\lim_{a\rightarrow 0}\phi_a(x)=\lim_{a\rightarrow 0}\frac{\partial^2}{\partial x^2}\phi_a(x)=f(x)[/tex]

I think this exchange of limiting operations requires that [itex]\phi_a(x)\rightarrow \phi(x)[/itex] uniformly. Provided f(x) either has compact support or falls off nicely for large x, and the increasing powers of a dominate any divergent behaviour in the derivatives of f(x), shouldn't we expect uniform convergence?

Does this method of solution imply that the key issue about proving the integral form for phi is related to our ability to exchange these limiting operations? We used the trick of substituting a nicer function phi_a in order to alter our task from taking derivatives of an integral, to exchanging derivatives and limits of sequences of functions.
 
  • #34
bdforbes said:
[tex]
\frac{\partial^2}{\partial x^2}\phi_a(x)=\int_{-\infty}^\infty \frac{1}{2}\frac{a^2}{\left[(x-x')^2+a^2\right]^{3/2}}f(x')dx'
[/tex]

At this point it is nicer to do a variable change [itex]x'=x+ar[/itex], so that the next integral is

[tex]
\int\limits_{-\infty}^{\infty} \frac{1}{2}\frac{1}{(r^2 + 1)^{3/2}} f(x+ar) dr
[/tex]

Then there will be no need for Taylor series of [itex]f[/itex].
 
  • #35
Been thinking about this problem a bit more, and I think I can give a fairly complete answer, using the method of distributions. Recall that we have [itex]\phi(x)=(-1/4\pi)\int f(y)/|x-y|\,dy[/itex] and want to show that it solves the pde [itex]\nabla^2\phi=f[/itex]. I assume that f has compact support for simplicity (we have to assume some boundary conditions, otherwise the integral won't be defined).

In my previous posts I showed that the PDE is satisfied in the sense of distributions. Jostpuur/bdforbes showed that the PDE is satisfied if we can assume that some limit commutes with the differentiation. In fact, limits always commute with differentiation if they are done in the sense of distributions, so we have the same result.

The question remains whether [itex]\phi[/itex] is twice continuously differentiable, so that the PDE will be satisfied in the standard sense and not just 'in distribution'. If f is continuously differentiable, then this will be the case (see my previous post). What about if f is continuous, but not necessarily differentiable?
One method is to smooth f by convolving it with a smooth function which is close to a dirac delta. Inside the integral, this is the same as convolving 1/|x-y| with a smooth function and reduces to the same method as jostpuur/bdforbes, where the smooth dirac delta approximation is [itex]a(r^2+a^2)^{-3/2}/2[/itex] (in 1d, at least), and you then need to know whether the limit a->0 commutes with differentiation.

So the question still remains as to whether a continuous f means that [itex]\phi[/itex] is twice continuously differentiable. In fact, this is false.

We know the following cases.
i) For a ball of uniform density, it is a standard result that the gravitational field ([itex]=-\nabla\phi[/itex]) varies linearly with distance to the centre inside the sphere, then follows the inverse square rule outside. So its derivative jumps (and is undefined) at the boundary. This is not too bad, as the derivative is at least bounded and we already know that it must be discontinuous wherever f is.
ii) In a previous post I tried to construct a counterexample example with a non-differentiable spherically symmetric mass density, but [itex]\phi[/itex] was still twice continuously differentiable.

Now consider a solid object with uniform density f inside and f=0 outside. If the surface is smooth, then the behaviour is like in (i) above with discontinuous but bounded second derivatives for [itex]\phi[/itex].
However, if the boundary has any ridges or spikes, then at these points the second derivative of [itex]\phi[/itex] will diverge to infinity. Can prove this in a separate post (this one is getting a bit long).

Now suppose that you have a solid body and again that the density f is 0 outside, but inside we have the density decreasing suddenly (but continuously) to 0 as you approach the boundary. Decreasing as one over the log will do the trick,

[tex]
f(x) = 1/\max(1,\log(-r))\textrm{\ \ (r=distance to edge of body)}
[/tex]

Then it is still the case that the derivative of the gravitational field (second derivative of [itex]\phi[/itex]) will diverge to infinity along any ridges or spikes on the surface. Have to finish here, but can show this in a separate post.
 
Last edited:
<h2>1. What is the Green's function approach using Lebesgue integration?</h2><p>The Green's function approach using Lebesgue integration is a mathematical method used in solving differential equations. It involves using the Lebesgue integral to calculate the Green's function, which is a function that satisfies a particular differential equation with a given set of boundary conditions.</p><h2>2. How is the Green's function calculated using Lebesgue integration?</h2><p>The Green's function is calculated by solving the differential equation with a delta function as the input, and then taking the Lebesgue integral of the resulting solution. This integral gives the Green's function, which can then be used to solve the original differential equation with any input function.</p><h2>3. What are the advantages of using Lebesgue integration in the Green's function approach?</h2><p>Lebesgue integration allows for a more general and flexible approach to calculating the Green's function. It can handle a wider range of functions and boundary conditions compared to other methods, such as the Fourier transform or Laplace transform.</p><h2>4. What are some applications of the Green's function approach using Lebesgue integration?</h2><p>The Green's function approach using Lebesgue integration has many applications in physics, engineering, and other fields. It is commonly used in solving problems involving heat transfer, fluid dynamics, and electrical circuits, among others.</p><h2>5. Are there any limitations to using the Green's function approach with Lebesgue integration?</h2><p>One limitation is that it can be more computationally intensive compared to other methods, such as the Laplace transform. It also requires a good understanding of Lebesgue integration and its properties, which may be challenging for some users.</p>

1. What is the Green's function approach using Lebesgue integration?

The Green's function approach using Lebesgue integration is a mathematical method used in solving differential equations. It involves using the Lebesgue integral to calculate the Green's function, which is a function that satisfies a particular differential equation with a given set of boundary conditions.

2. How is the Green's function calculated using Lebesgue integration?

The Green's function is calculated by solving the differential equation with a delta function as the input, and then taking the Lebesgue integral of the resulting solution. This integral gives the Green's function, which can then be used to solve the original differential equation with any input function.

3. What are the advantages of using Lebesgue integration in the Green's function approach?

Lebesgue integration allows for a more general and flexible approach to calculating the Green's function. It can handle a wider range of functions and boundary conditions compared to other methods, such as the Fourier transform or Laplace transform.

4. What are some applications of the Green's function approach using Lebesgue integration?

The Green's function approach using Lebesgue integration has many applications in physics, engineering, and other fields. It is commonly used in solving problems involving heat transfer, fluid dynamics, and electrical circuits, among others.

5. Are there any limitations to using the Green's function approach with Lebesgue integration?

One limitation is that it can be more computationally intensive compared to other methods, such as the Laplace transform. It also requires a good understanding of Lebesgue integration and its properties, which may be challenging for some users.

Similar threads

  • Differential Equations
Replies
1
Views
2K
  • Differential Equations
Replies
1
Views
1K
Replies
6
Views
2K
  • Differential Equations
Replies
2
Views
2K
  • Differential Equations
Replies
33
Views
5K
  • Advanced Physics Homework Help
Replies
3
Views
287
  • Differential Equations
Replies
4
Views
1K
  • Differential Equations
Replies
1
Views
976
Replies
8
Views
428
  • Differential Equations
Replies
3
Views
1K
Back
Top