# Modified Helmholtz Equation in Spherical Coordinates

Writing out the Modified Helmholtz equation in spherically symmetric co-ordinates

Note that $$\nabla^2 \psi(r)\;$$=$$\;\frac{d^{2} \psi}{d r^{2}}+\frac{2}{r} \frac{d \psi}{d r}$$=$$\frac{1}{r} \frac{d^{2}}{d r^{2}}(r \psi)$$ in the spherically symmetric case.

LapSphSym[F[r]] := (1/r)*D[r*F[r], {r, 2}]
eqmHrr1 = LapSphSym[F[r]] - m^2*F[r] == DiracDelta[r - r1]
Assuming[Element[{r, r1, F, m}, Reals], DSolve[eqmHrr1, F[r], r]]


MMa 12.0 output

(** {{F[r] -> C/(E^(m*r)*r) + (E^(m*r)*C)/r -
(E^((-m)*r - m*r1)*(-E^(2*m*r) + E^(2*m*r1))*r1*HeavisideTheta[r - r1])/
(2*m*r)}} **)


Questions

 Is this the correct approach, including DiracDelta[r-r1] , if "r" is a (spherically symmetric) vector?

 How to input the boundary conditions to get the Green's functions?

• Somewhat related: mathematica.stackexchange.com/q/209662/1871 As shown there, GreenFunction knows how to handle the problem in Cartesian coordinates, but doesn't seem to know how to handle the spherical coordinates, at least now. – xzczd Jul 9 at 1:52

A partial attempt is documented below, but it needs further work to get at the solutions.

eqmHrr1 = -LapSphSym[F[r]] - m^2*F[r] == (1/r^2)*DiracDelta[r - r1];
sol = TrigExpand[Assuming[Element[{r, r1, F[r], m}, Reals], DSolve[eqmHrr1, F[r], r]]]
FullSimplify[Simplify[sol, Assumptions -> {0 < r < Infinity && r1 < r < Infinity}]]
(** {{F[r] -> (C/E^(I*m*r) + E^(I*m*r)*C -
Sin[m*(r - r1)]/(m*r1))/r}}**)


How does one enforce limits at infinity, not include the origin, so that C and C can be determined? The correct solution(s) is/are: $$\displaystyle{\frac{e^{\pm \iota~ m~|r - r1|}}{4\pi |r -r1|}}$$

• (1/r^2)*DiracDelta[r - r1] it is not the same as in 1D case. It describes singularity distributed on a sphere r=r1. It can be electric charge on a soap bubble for instance. – Alex Trounev Jul 13 at 9:54
• Thanks. How should the DiracDelta term written correctly, i.e. is there a way to write DiracDlelta[$\vec{r}$ - $\vec{r1}$] using "r" as a scalar in spherically symmetric coordinates? Alternatively, what's the correct way to cast this equation so the solution can be validated in 3D within Mma? Just in case, It's not essential to use "r" as a scalar. – arny Jul 13 at 19:28
• There is a good answer to your question on math.stackexchange.com/questions/2235615/… – Alex Trounev Jul 13 at 20:19
• Also there is a good explanation about 3D Dirac delta on farside.ph.utexas.edu/teaching/jk1/Electromagnetism.pdf – Alex Trounev Jul 13 at 20:33