REVIEW 3 major objections 5 minor 42 references
Wavelet based solutions to the Poisson and the Helmholtz equations on the $n$-dimensional unit sphere
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Wavelet inversion gives explicit Green functions for Poisson and Helmholtz equations on the n-dimensional sphere.
desk verdict A genuinely new wavelet derivation and useful new closed forms, but the Green-function series needs a distributional framework and the main proof has an unjustified interchange. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the zonal Green function $G$, a function of the angular distance on $S^n$ represented by its Gegenbauer expansion; the coefficients $1/(a - l(n+l-1))$ are chosen so that convolution inverts the Laplace-Beltrami operator, whose eigenvalues are $-l(n+l-1)$. The argument is carried by the continuous spherical wavelet transform based on approximate identities: applying $\Delta_* + a$ to the analyzing Poisson wavelet produces an admissible wavelet family, and the reconstruction formula converts the differential equation into a convolution with $G$. The closed-form part relies on the Poisson kernel identity $\Sigma_n p_r(\cos\theta) = \sum_{l=0}^\infty r^l \frac{\lambda+l}{\lambda} C_l^\lambda(\cos\theta)$ together with the integral identity $\int_0^1 R^{-(n+2L)} \int_0^R r^{n+L-2} r^l\, dr\, dR = -\frac{1}{(L-l)(n+L+l-1)}$, which reproduces the desired coefficient exactly.
What would settle it
For $n=2$ and $a=0$, set $t=1$ in the series (16): the $l$-th term is $-\frac{2l+1}{l(l+1)} \sim -2/l$, so the series diverges logarithmically and cannot be read as an ordinary pointwise sum. A decisive test is to verify the convolution identity $u = f*G$ in a specified distributional or Sobolev sense for smooth test data with known solutions; if the truncated convolutions do not converge to the exact solution in that norm, the claimed solution formula fails.
Extended reading notes
Core claim
The paper claims that if $f \in C(S^n)$ and $u \in C^2(S^n)$ solve $\Delta_* u + a u = f$ with $a$ not equal to any eigenvalue $l(n+l-1)$, then $u = f * G$, where $G$ is the zonal series $G = \sum_{l=0}^\infty \frac{1}{a - l(n+l-1)} \frac{\lambda+l}{\lambda} C_l^\lambda$, with $\lambda = (n-1)/2$ and $C_l^\lambda$ the Gegenbauer polynomials. For the Poisson equation with zero-mean data the same formula holds starting at $l=1$ with coefficient $-1/(l(n+l-1))$. For $a = L(n+L-1)$ with non-integer $L$, Theorem 4.3 expresses $G$ as a double integral of the Poisson kernel minus the first $L_0$ Gegenbauer terms, plus a finite rational correction; this is the mechanism that yields closed forms. The paper states that for tested cases these closed forms coincide with known Green functions, including the classical logarithmic kernel on $S^2$, and that the negative-$a$ cases have not appeared in closed form before.
Load-bearing premise
The formulas work only if the formal Gegenbauer series for $G$ is given a meaning beyond pointwise convergence, because at the north pole the series diverges logarithmically even in the classical two-dimensional Poisson case.
Editorial extensions
If this is right
- For every dimension $n \ge 2$ and every $a$ that is not an eigenvalue, the Helmholtz solution on $S^n$ is given by a single convolution formula, so the problem is reduced to evaluating one zonal kernel.
- For $a=0$ the Poisson Green function has explicit closed forms in dimensions $2$ through $10$, including the standard logarithmic kernel on the two-sphere.
- For Helmholtz parameters $a = L(n+L-1)$ with integer $L$, the formula with the critical term removed gives the generalized Green function and reproduces previously known closed forms.
- For negative $a$, such as $a = -3/4$ on $S^3$ or $a = -2$ on $S^4$, the closed forms in Table 4 are presented as new, with applications to telegraph and Klein-Gordon type equations.
- The derivation provides an analytical alternative to numerical wavelet solvers, producing exact series and closed integrals rather than approximate solutions.
Reading between the lines
- If the distributional interpretation of the Green series is made explicit and rigorous, the same wavelet-inversion trick should produce Green functions for other rotation-invariant operators on $S^n$, such as polyharmonic operators or fractional powers of $\Delta_*$.
- The double-integral representation in Theorem 4.3 suggests a quadrature-based numerical evaluation that avoids summing the divergent Gegenbauer series, which would be stable away from the singular diagonal.
- Since the integral representation is valid for all non-integer $L$, the tabulated closed forms could be extended systematically to arbitrary rational $a$ by symbolic integration.
- Because the kernel is zonal, the method should extend to non-zonal right-hand sides through the paper's convolution framework, which already handles general $L^2$ functions by rotation-invariant kernels.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a wavelet-transform method for solving the Poisson and Helmholtz equations on the unit n-sphere. It derives convolution representations u = f * G with G given by a Gegenbauer series (Theorem 4.1) and gives a double-integral closed form for G when a = L(n+L-1) with L non-integer (Theorem 4.3), together with extensive tables of explicit formulas for n = 2,...,10 and various values of a. The method is based on the continuous spherical wavelet transform of the authors' earlier work, and the closed forms are claimed to coincide with known results by Szmytkowski for the positive integer resonance cases.
Significance. If the results are made rigorous, they would provide a new wavelet-based derivation of Green functions on spheres and give closed forms for Helmholtz parameters, in particular negative values of a, that have not been published before. The paper is unusual in offering analytical rather than numerical wavelet solutions, and the tabulated formulas are potentially valuable reference material. However, the manuscript's central analytical objects—the Green function series—are not functions in general, so the stated theorems are formal as written and require a distributional or Sobolev-space interpretation. The paper does supply explicit, verifiable formulas for many cases, which is a strength even where the proofs are incomplete.
major comments (3)
- [Eq. (16) and the proof of Theorem 4.1] The series (16) defining G does not converge pointwise for n >= 2: at t = 1 the terms behave like l^{n-3}, and the zonal L^2 norm of partial sums behaves like sum l^{n-5}, which diverges for n >= 5. For n = 2 the series diverges logarithmically at t = 1. Thus G is not generally an ordinary function, and u = f * G in (15) is not a classical convolution between a continuous function and a function kernel. The proof's claim that 'L2-convergence of the triple integral is ensured' does not establish convergence of the final series (16). The manuscript must state the distributional or Sobolev framework in which (15) holds, and prove the identification of G with that object.
- [Proof of Theorem 4.3, Eqs. (24)-(26)] The interchange of the series (24) with the double integration in (23) is not justified. Although (24) converges absolutely for each r in [0,1), the integrated series is bounded by sum_l ((lambda+l)/lambda) |C_l^lambda(cos theta)| / l^2, which behaves like sum l^{n-3} for n >= 3 and diverges logarithmically for n = 2. Dominated convergence therefore fails. The equality of the double integral with the series should be proved in a distributional sense or by an explicit summation method (such as Abel summation), and this is essential to the validity of the closed-form formula (23).
- [Remark 4.2(3) and Corollary 4.4] The paper acknowledges in Remark 4.2(3) that [31] treats convergence of such series in spherical Sobolev spaces 'more precisely,' which indicates the gap is known to the authors. Corollary 4.4 inherits the same issue for integer L: the subtracted series (28) still fails to converge classically for n >= 2. The manuscript should either incorporate a Sobolev-space or distributional framework and prove the convolution statement in that framework, or explicitly state that all Green-function identities are understood distributionally and supply the corresponding convolution theorem for continuous f.
minor comments (5)
- [Introduction, paragraph 2] The name 'Klein–Gordan' should be 'Klein–Gordon'.
- [Tables 3 and 4] The placement of the factor pi in the tabulated expressions is visually ambiguous, because it appears after a large fraction without a clear denominator; please reformat so that the formulas are unambiguous, for example with an explicit multiplication sign or by placing pi in the numerator or denominator consistently.
- [Definition 3.1] In condition (11), the l = 0 statement is written as an integral equal to zero; it would be clearer to say explicitly that the pair has zero mean, as is done in (12).
- [Proof of Theorem 4.1, last paragraph] The sentence 'If integral u is not 0 and a is not 0, it follows directly from (14)' skips the verification that the l = 0 mode of G is 1/a; adding one line showing that the constant mode is handled correctly would improve clarity.
- [Theorem 4.3] The notation L0 := max{[L], [-n-L+1]} with [x] denoting the integer part is functional but slightly terse; a short parenthetical explanation of why this choice of L0 removes the singular modes would help the reader.
Circularity Check
No significant circularity: the Green functions are derived by wavelet inversion and algebraic spectral evaluation, not assumed or fitted; self-citations are background only.
full rationale
The derivation chain for Theorem 4.1 is self-contained in the relevant sense. The paper does not postulate u=f*G; it starts from the eigenfunction identity Δ*Y_l^k=-l(n+l-1)Y_l^k, forms the wavelet family Θρ=(Δ*+ā)Ψρ, and computes the Gegenbauer coefficients of G=∫Ψρ*Ωρ dρ/ρ as 1/(a-l(n+l-1))·(λ+l)/λ. This is a direct algebraic consequence of the chosen reconstruction wavelet, not an input fitted to the target solution. Theorem 4.3 is likewise an identity transformation: it rewrites the same spectral series using the Poisson-kernel expansion (13) and the elementary integral (26); the right-hand side of (23) is shown to equal the series (22), so no target conclusion is smuggled into the hypotheses. The paper's self-citations [22-26] supply background wavelet-transform definitions and inversion results, but those results do not contain the Green-function formula, and the tables are checked against independent sources ([13,12,31,37,38]). There is a genuine rigor gap, not a circularity: the paper does not state the distributional or Sobolev setting in which the series (16)/(22) converges, and the bound (27) does not justify the term-by-term integration in Theorem 4.3, since the absolute terms behave like l^(n-3) at t=1 and diverge for n≥5. The authors themselves acknowledge a related limitation in Remark 4.2(3), citing [31] for more precise convergence considerations. This gap affects correctness and rigor, but it does not make the derivation circular: the proof does not assume the theorem it is proving.
Assumptions & free parameters
assumptions (4)
- standard math Hyperspherical harmonics are eigenfunctions of the Laplace-Beltrami operator: Delta*Y_l^k = -l(n+l-1)Y_l^k (Eq. 6).
- standard math The spherical wavelet transform is invertible for zero-mean L2 functions with the admissible-pair condition (11).
- standard math The Gegenbauer polynomial bound |C_l^lambda(cos theta)| <= (n+l-2)^{n-2} (Eq. 27) holds uniformly.
- ad hoc to paper Convergence of the Green-function series (16)/(18) is not specified for n>=2; the paper assumes it defines an object suitable for convolution with continuous f.
Cite this review
Pith. "Pith review of Wavelet based solutions to the Poisson and the Helmholtz equations on the $n$-dimensional unit sphere." pith.science (2026). https://pith.science/paper/AUIVE6E6
@misc{pith2026250703451,
author = {Pith},
title = {Pith review of: Wavelet based solutions to the Poisson and the Helmholtz equations on the $n$-dimensional unit sphere},
year = {2026},
howpublished = {\url{https://pith.science/paper/AUIVE6E6}},
note = {Machine review of arXiv:2507.03451}
}
abstract
We present a method of solving partial differential equations on the $n$-dimensional unit sphere using methods based on the continuous wavelet transform derived from approximate identities. We give an explicit analytical solution to the Poisson equation and to the Helmholtz equations. For the first one and for some special values of the parameter in the latter one, we derive a closed formula for the generalized Green function.
Reference graph
Works this paper leans on
-
[31]
V. Michel, Geomathematics – modelling and solving mathematical problems in geodesy and geo- physics, Cambridge University Press, Cambridge, 2022,https://doi.org/10.1017/9781108297882
-
[1]
R. Amin, B. Alshahrani, M. Mahmoud, A.-H. Abdel-Aty, K. Shah, W. Deebanif, Haar wavelet method for solution of distributed order time-fractional differential equations, Alexandria Engineering J. 60(3), 3295–3303 (2021), https://doi.org/10.1016/j.aej.2021.01.039
-
[2]
J.-P. Antoine, P. Vandergheynst, Wavelets on the 2-sphere: a group-theoretical approach , Appl. Comput. Harmon. Anal. 7(3), 262–291 (1999), https://doi.org/10.1006/acha.1999.0272
-
[3]
J.-P. Antoine, P. Vandergheynst, Wavelets on the n-sphere and related manifolds , J. Math. Phys. 39(8), 3987–4008 (1998), https://doi.org/10.1063/1.532481
-
[4]
V.M. Babiˇ c, M.B. Kapileviˇ c, S.G. Mihlin, G.I. Natanson, P.M. Riz, L.N. Slobodecki˘i, M.M. Smirnov, Line ˘ ınye uravneniya matematiˇ cesko ˘ ı fiziki(Russian) [The linear equations of mathematical physics] Edited by S.G. Mihlin, Izdat. ”Nauka”, Moscow, 1964
work page 1964
-
[5]
S. Balaji, G. Hariharan, A novel wavelet approximation method for the solution of nonlinear dif- ferential equations with variable coefficients arising in astrophysics , Astrophys Space Sci. 363(16) (2018), https://doi.org/10.1007/s10509-017-3236-3
-
[6]
S. Bernstein, Spherical singular integrals, monogenic kernels and wavelets on the three– dimensional sphere, Adv. Appl. Clifford Algebr. 19(2), 173–189 (2009), https://doi.org/10.1007/ s00006-009-0149-4
work page 2009
-
[7]
S. Ebert, S. Bernstein, P. Cerejeiras, U. K´ ahler,Nonzonal wavelets on S N , 18th International Confer- ence on the Application of Computer Science and Mathematics in Architecture and Civil Engineering, Weimar 2009, https://doi.org/10.25643/bauhaus-universitaet.2840
Show all 42 references
-
[8]
Edelen, Helmhotz equation on the surface of the unit (N + 1)-sphere, J
D.G.B. Edelen, Helmhotz equation on the surface of the unit (N + 1)-sphere, J. Math. Anal. Appl. 23, 99–109 (1968), https://doi.org/10.1016/0022-247X(68)90117-0
1968 doi
-
[9]
Fengler, A nonlinear Galerkin scheme involving vectorial and tensorial spherical wavelets for solving the incompressible Navier-Stokes equation on the sphere , PAMM - Proc
M.J. Fengler, A nonlinear Galerkin scheme involving vectorial and tensorial spherical wavelets for solving the incompressible Navier-Stokes equation on the sphere , PAMM - Proc. Appl. Math. Mech. 5, 457–458 (2005), https://doi.org/10.1002/pamm.200510205
2005 doi
-
[10]
Frazier, An introduction to wavelets through linear algebra , Springer, New York, 1999, https: //doi.org/10.1007/b97841
M.W. Frazier, An introduction to wavelets through linear algebra , Springer, New York, 1999, https: //doi.org/10.1007/b97841
1999 doi
-
[11]
Freeden, T
W. Freeden, T. Gervens, M. Schreiner, Constructive approximation on the sphere. With ap- plications to geomathematics , Clarendon Press, New York, 1998, https://doi.org/10.1017/ S0013091500020836
1998
-
[12]
Freeden, M
W. Freeden, M. Gutting, Special functions of mathematical (geo-)physics , Applied and Numeri- cal Harmonic Analysis, Birkh¨ auser/Springer Basel AG, Basel, 2013, https://doi.org/10.1007/ 978-3-0348-0563-6
2013
-
[13]
Freeden, M
W. Freeden, M. Schreiner, Spherical functions of mathematical geosciences – a scalar, vectorial, and tensorial setup, Springer, Berlin Heidelberg, 2009, https://doi.org/10.1007/978-3-540-85112-7 W A VELET BASED SOLUTIONS TO THE POISSON AND THE HELMHOLTZ EQUATIONS ON THEn-DIMEN...
2009 doi
-
[14]
Freeden, U
W. Freeden, U. Windheuser, Combined spherical harmonic and wavelet expansion – a future concept in Earth’s gravitational determination , Appl. Comput. Harmon. Anal. 4(1), 1–37 (1997), https: //doi.org/10.1006/acha.1996.0192
1997
-
[15]
Freeden, U
W. Freeden, U. Windheuser, Spherical wavelet transform and its discretization, Adv. Comput. Math. 5(1), 51–94 (1996), https://doi.org/10.1007/BF02124735
1996 doi
-
[16]
Gittelson, Adaptive wavelet methods for elliptic partial differential equations with random operators, Numer
C.J. Gittelson, Adaptive wavelet methods for elliptic partial differential equations with random operators, Numer. Math. 126(3), 471–513 (2014), https://www.infona.pl/resource/bwmeta1. element.springer-1c047140-43e6-3ce6-b250-9d984995c541
2014
-
[17]
Grebenitcharsky and M.G
R. Grebenitcharsky and M.G. Sideris, Application of spherical pseudo-differential operators and spherical wavelets for numerical solutions of the fixed altimetry-gravimetry boundary value problem , V Hotine-Marussi Symposium on Mathematical Geodesy, 127, 242–251 (2004), https:...
2004 doi
-
[18]
Gupta, S
A.K. Gupta, S. Saha Ray Wavelet methods for solving fractional order differential equations , Math. Probl. Eng. (2014), https://doi.org/10.1155/2014/140453
2014 doi
-
[19]
Hajji, S
M.A. Hajji, S. Melkonian, R. Vaillancourt, Two-dimensional wavelet bases for partial differential operators and applications , Advances in pseudo-differential operators, Oper. Theory Adv. Appl., 155, 219–233 (2004), https://doi.org/10.1007/978-3-0348-7840-1_13
2004 doi
-
[20]
B. Han, M. Michelle, Y.S. Wong, Wavelet-based methods for numerical solutions of differential equa- tions, https://arxiv.org/abs/1909.12192
1909 arXiv
-
[21]
Hariharan, K
G. Hariharan, K. Kannan, Haar wavelet method for solving some nonlinear parabolic equations , J. Math. Chem. 48(4), 1044–1061 (2010), https://doi.org/10.1007/s10910-010-9724-0
2010 doi
-
[22]
Iglewska-Nowak, A continuous spherical wavelet transform for C(Rn), Appl
I. Iglewska-Nowak, A continuous spherical wavelet transform for C(Rn), Appl. Comput. Harmon. Anal. 47(3), 1033–1039 (2019), https://doi.org/10.1016/j.acha.2019.04.001
2019 doi
-
[23]
Iglewska-Nowak, Continuous wavelet transforms on n-dimensional spheres , Appl
I. Iglewska-Nowak, Continuous wavelet transforms on n-dimensional spheres , Appl. Comput. Har- mon. Anal. 39(2), 248–276 (2015), https://doi.org/10.1016/j.acha.2014.09.006
2015 doi
-
[24]
Iglewska-Nowak, Directional wavelets on n-dimensional spheres , Appl
I. Iglewska-Nowak, Directional wavelets on n-dimensional spheres , Appl. Comput. Harmon. Anal. 44(2), 201–229 (2018), https://doi.org/10.1016/j.acha.2016.04.008
2018 doi
-
[25]
Iglewska-Nowak, Frames of directional wavelets on n–dimensional spheres , Appl
I. Iglewska-Nowak, Frames of directional wavelets on n–dimensional spheres , Appl. Comput. Har- mon. Anal. 43(1), 148–161 (2017), https://doi.org/10.1016/j.acha.2016.01.004
2017 doi
-
[26]
Iglewska-Nowak, Poisson wavelets on n-dimensional spheres, J
I. Iglewska-Nowak, Poisson wavelets on n-dimensional spheres, J. Fourier Anal. Appl.21(1), 206–227 (2015), https://doi.org/10.1007/s00041-014-9366-x
2015 doi
-
[27]
S. Jaffard, Wavelets and analysis of partial differential equations , Probabilistic and stochas- tic methods in analysis, with applications, 372, 3–13 (1992), https://doi.org/10.1007/ 978-94-011-2791-2_1
1992
-
[28]
Lepik, H
¨U. Lepik, H. Hein, Haar wavelets , Mathematical Engineering. Springer Cham, 2014, https://doi. org/10.1007/978-3-319-04295-4 22 ILONA IGLEWSKA-NOW AK AND PIOTR STEF ANIAK
2014 doi
-
[29]
Mehra, N
M. Mehra, N. K.-R. Kevlahan, An adaptive multilevel wavelet solver for elliptic equations on an optimal spherical geodesic grid , SIAM J. Sci. Comput. 30(6), 3073–3086 (2008), https://doi.org/ 10.1137/070689607
2008 doi
-
[30]
Mei, Q.S
S.L. Mei, Q.S. Lu, S.W. Zhang, and L. Jin, Adaptive interval wavelet precise integration method for partial differential equations , Appl. Math. Mech. 26(3), 364–371 (2005), https://doi.org/10. 1007/BF02440087
2005
-
[32]
Perrier, Towards a method for solving partial differential equations using wavelet bases , Wavelets
V. Perrier, Towards a method for solving partial differential equations using wavelet bases , Wavelets. Inverse Probl. Theoret. Imaging, 269–283 (1989),https://doi.org/10.1007/978-3-642-97177-8_ 27
1989 doi
-
[33]
Rehman, D
M. Rehman, D. Baleanu, J. Alzabut, M. Ismail, U. Saeed, Green–Haar wavelets method for gen- eralized fractional differential equations , Adv. Difference Equ. (2020), https://doi.org/10.1186/ s13662-020-02974-6
2020
-
[34]
J. Sun, X. Yi, B. Ye, and Y. Shen, Wavelet-Galerkin Solutions for Differential Equations , Wuhan Univ. J. Nat. Sci. 3(4), 403–406 (1998), https://doi.org/10.1007/BF02830037
1998 doi
-
[35]
Shimakura, Partial differential operators of elliptic type , Translations of Mathematical Mono- graphs, Vol
N. Shimakura, Partial differential operators of elliptic type , Translations of Mathematical Mono- graphs, Vol. 99, Amer. Math. Soc., Providence, Rhode Island, 1992
1992
-
[36]
Szeg¨ o,Orthogonal polynomials, Fourth Edition, AMS Coll
G. Szeg¨ o,Orthogonal polynomials, Fourth Edition, AMS Coll. Publ., Providence, Rhode Island, 1975
1975
-
[37]
Szmytkowski, Closed form of the generalized Green ’s function for the Helmholtz operator on the two-dimensional unit sphere , J
R. Szmytkowski, Closed form of the generalized Green ’s function for the Helmholtz operator on the two-dimensional unit sphere , J. Math. Phys. 47(6) (2006), https://doi.org/10.1063/1.2203430
2006 doi
-
[38]
Szmytkowski, Closed forms of the Green ’s function and the generalized Green ’s function for the Helmholtz operator on the N -dimensional unit sphere , J
R. Szmytkowski, Closed forms of the Green ’s function and the generalized Green ’s function for the Helmholtz operator on the N -dimensional unit sphere , J. Phys. A 40(5), 995–1009 (2007), https: //doi.org/10.1088/1751-8113/40/5/009
2007 doi
-
[39]
Urban, Wavelet methods for elliptic partial differential equations, Numerical Mathematics and Sci- entific Computation
K. Urban, Wavelet methods for elliptic partial differential equations, Numerical Mathematics and Sci- entific Computation. Oxford University Press, Oxford, 2009, https://doi.org/10.1093/acprof: oso/9780198526056.001.0001
2009
-
[40]
N. Ja. Vilenkin, Special functions and the theory of group representations , in Translations of Math- ematical Monographs, Vol. 22, American Mathematical Society, Providence, R. I., 1968
1968
-
[41]
L. Wang, Y. Ma, and Z. Meng, Haar wavelet method for solving fractional partial differential equa- tions numerically , Appl. Math. Comput. 227, 66–76 (2014), https://doi.org/10.1016/j.amc. 2013.11.004
2014 doi
-
[42]
S. Zhi, Y. Xu, and J. Zhao, Haar wavelets method for solving Poisson equations with jump conditions in irregular domain , Adv. Comput. Math. 42(4), 995–1012 (2016), https://doi.org/10.1007/ s10444-015-9450-z W A VELET BASED SOLUTIONS TO THE POISSON AND THE HELMHOLTZ EQUATIONS ...
2016
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.