Differences between fundamental solutions of general higher order elliptic operators and of products of second order operators

We study fundamental solutions of elliptic operators of order 2m≥4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2m\ge 4$$\end{document} with constant coefficients in large dimensions n≥2m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\ge 2m$$\end{document}, where their singularities become unbounded. For compositions of second order operators these can be chosen as convolution products of positive singular functions, which are positive themselves. As soon as n≥3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\ge 3$$\end{document}, the polyharmonic operator (-Δ)m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(-\Delta )^m$$\end{document} may no longer serve as a prototype for the general elliptic operator. It is known from examples of Maz’ya and Nazarov (Math. Notes 39:14–16, 1986; Transl. of Mat. Zametki 39, 24–28, 1986) and Davies (J Differ Equ 135:83–102, 1997) that in dimensions n≥2m+3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\ge 2m+3$$\end{document} fundamental solutions of specific operators of order 2m≥4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2m\ge 4$$\end{document} may change sign near their singularities: there are “positive” as well as “negative” directions along which the fundamental solution tends to +∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$+\infty $$\end{document} and -∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-\infty $$\end{document} respectively, when approaching its pole. In order to understand this phenomenon systematically we first show that existence of a “positive” direction directly follows from the ellipticity of the operator. We establish an inductive argument by space dimension which shows that sign change in some dimension implies sign change in any larger dimension for suitably constructed operators. Moreover, we deduce for n=2m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2m$$\end{document}, n=2m+2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2m+2$$\end{document} and for all odd dimensions an explicit closed expression for the fundamental solution in terms of its symbol. From such formulae it becomes clear that the sign of the fundamental solution for such operators depends on the dimension. Indeed, we show that we have even sign change for a suitable operator of order 2m in dimension n=2m+2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2m+2$$\end{document}. On the other hand we show that in the dimensions n=2m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2m$$\end{document} and n=2m+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2m+1$$\end{document} the fundamental solution of any such elliptic operator is always positive around its singularity.


General constant coefficients elliptic operators
We focus our attention to uniformly elliptic operators of order 2m with constant coefficients which involve only the highest order derivatives, namely where the 2m-homogeneous characteristic polynomial is called (possibly up to a sign) the symbol of the operator. Uniform ellipticity means then that Q is strictly positive on the unit sphere, i.e. there exists a constant λ > 0 such that ∀ξ ∈ R n : Q(ξ ) ≥ λ|ξ | 2m .

Fundamental solutions
In order to construct and to understand solutions u to the differential equation Lu = f for a given right-hand side f , one introduces the concept of a fundamental solution K L (x, . ) for any "pole" x ∈ R n which is defined as a solution to the equations L * K L (x, . ) = δ x and L K L ( . , x) = δ x in the distributional sense where δ x is the δ-distribution located at x. This means that for any test function ψ ∈ C ∞ 0 (R n ) one has R n K L (x, y)Lψ(y) dy = ψ(x), being the adjoint operator of L. Because L has only constant coefficients and only of the highest even order 2m, we have that L = L * . Moreover, we may achieve that K L (x, y) = K L (0, x − y) = K L (0, y − x) = K L (y, x).
For given f ∈ C ∞ 0 (R n ), any fundamental solution yields a solution to the differential equation Lu = f in R n by putting u(x) := R n K L (x, y) f (y) dy.
One should also notice that, if a fundamental solution exists, it is not unique: one may add any smooth solution of Lv = 0, namelyK L (x, y) = K L (x, y) + v(x − y) yields another fundamental solution.

Green functions
When the problem Lu = f is considered in a sufficiently smooth bounded domain, one may still obtain solution and even representation formulae by means of suitable fundamental solutions. Indeed, let ⊂ R n be a bounded smooth domain and consider the problem where f ∈ C 0,γ ( ) and the boundary conditions verify a complementing condition, see [1]. As a typical and most frequently studied prototype one may think of Dirichlet boundary conditions Then the unique solution of (3) is given by Notice that in general it is not straightforward to infer the existence of such h L, ,B . However, exploiting the general elliptic theory of Agmon, Douglis, and Nirenberg [1] this is always possible in our special case when the operator L has only constant coefficients of highest order, if Dirichlet boundary conditions B = B D are imposed and the domain is C 2m,γ -smooth.
In this case, one also infers by standard estimates that the function h L, ,B is regular in . Since in large dimensions fundamental solutions have a singularity near the pole, it becomes clear that, in order to understand G L, ,B , we need first to investigate the behaviour of fundamental solutions.

Positivity questions
Positivity properties for G L, ,B concern the question whether a positive right-hand side yields a positive solution: if u is a solution of (3), does it hold that f ≥ 0 ⇒ u ≥ 0? One often expects such a behaviour for physical or geometrical reasons. However, for equations of order at least 4, such a positivity preserving property will fail in general, see [7] for historical remarks and detailed references. This question concerns a nonlocal behaviour of the full boundary value problem and often the influence of boundary conditions spoils the expected positivity. However, physically, one would hope that when applying an extremely concentrated right-hand side-a δ-distribution-then close to this point the solution should respond in the same direction. This leads to the related but relaxed local question: Is a suitable fundamental solution to the differential equation positive, at least close to its pole? This question is reasonable only for large dimensions n ≥ 2m because only here, fundamental solutions become unbounded and they are unique only up to locally bounded regular solutions of the homogeneous equation. If n > 2m one may achieve uniqueness of the fundamental solution by imposing zero (Dirichlet) boundary conditions at infinity. In this case K L may be considered as the Green function G L,R n ,B D in the whole space. This means that one considers just the behaviour of the differential equation and disregards the influence of possible boundary conditions (being infinitely far apart).

Previous results
In the context of second order equations (m = 1), both local and nonlocal behaviours are well established. Indeed, within the class of constant coefficients operators, the Laplacian − is, up to a change of coordinates, the only such operator. Its fundamental solutions are known explicitly and in particular they are positive (if n = 2, at least close to the pole). Moreover, the maximum principle holds for such operators, so positive data yield positive solutions (see [8]). In other words, the Green function is always positive.
When one moves to the higher order setting (m ≥ 2), several differences arise, even for (− ) m or, equivalently, for powers of second order operators with constant coefficients.
Indeed, if one investigates the positivity preserving property in bounded domains, then the answer is largely affected by the choice of boundary conditions. As an example, on the one hand, with Navier boundary conditions (u = u = · · · = m−1 u = 0) one may rewrite the problem as a second order system and thus the maximum principle implies positivity. On the other hand, this tool is in general not available when dealing with Dirichlet boundary conditions (u = ∂ ν u = · · · = ∂ m−1 ν u = 0) and one cannot expect positivity, in general not even in convex bounded smooth domains (see [5]). Nevertheless, positivity holds in balls and their small smooth deformations (see [2,9]). We refer to [7] for an extensive survey of the topic.
However, within that class of powers of second order operators, if one restricts to a "local" question, meaning the positivity of Green functions under Dirichlet boundary conditions near the pole, the answer is still affirmative. Indeed, a uniform local positivity can be proved, namely the existence of constants r m, > 0, δ m, > 0 such that This means that the negative part and the singularity of the Green function are uniformly apart.
A consequence of that result is that the size of negative part of the Green function, if present at all, is small compared to its positive part. Indeed, concerning Dirichlet problems, positivity for a rank-1-correction of the polyharmonic Green function is retrieved, namely where d denotes the distance to the boundary and c m, is a sufficiently large positive constant, see [9,10].
These results have been extended later on by Pulst in his PhD-dissertation [15] for formally selfadjoint positive definite operators of order 2m with the polyharmonic operator (− ) m or an m-th power of a second order elliptic operator with constant coefficients as the leading term. Lower order terms are permitted provided they can be written in divergence form and have sufficiently smooth and uniformly bounded coefficients.
In two dimensions, i.e. n = 2, the symbol Q with real coefficients can be split into 2m linear terms. Combining mutually conjugate pairs (ξ 1 + a k ξ 2 ) and (ξ 1 + a k ξ 2 ) of these linear terms with nonreal a k we see that is a product of second order symbols. However, in dimensions n > 2 powers of second order operators are not the prototype of a general operator L of order 2m, not even in the case of constant coefficients. Moreover, it is in general not possible to rewrite L as an m-fold composition of (possibly different) second order operators. Indeed, let us simply consider the case of a homogeneous fourth order operator with a symbol of the kind Q(x, y, z) = x 4 + y 4 + z 4 + i+ j+k=4 0≤i, j,k≤3 c i, j,k x i y j z k and suppose that it is the product of two second order polynomials q 1 , q 2 . One may assume that both polynomials have their coefficients in front of x 2 equal to 1 and then, they would necessarily be of the kind q 1 (x, y, z) = x 2 + cy 2 + dz 2 + a 1 x y + a 2 xz + a 3 yz The smooth map from (0, ∞) 2 × R 6 into the 12-dimensional vector space of such symbols Q which maps is not surjective.
Concerning explicit formulae and (local) positivity properties of fundamental solutions of such general elliptic operators only little is known. Existence of fundamental solution is shown in [12] in a very general framework, and rather involved formulae are obtained. In the particular case of a 2m-homogeneous higher order uniformly elliptic operator with constant coefficients, different implicit expressions have been found according to the parity of the dimension n. In what follows we always assume that n ≥ 2m.
For odd n, the general formula for a fundamental solution [12, (3.44)] simplifies as (from [12, (3.54)]), while for even n one has (see [12, (3.62)]) We recall that Q denotes the symbol (possibly up to a sign) of the operator L. On the other hand, motivated by questions in potential and Schrödinger semigroup theory, respectively, and without referring to (4), (5) or even [12], Maz'ya-Nazarov [14] and Davies [3] found examples of elliptic operators of order 2m ≥ 4 in dimensions n ≥ 2m+3 with sign changing fundamental solutions. The precise range of dimensions where this phenomenon may be observed remained open as well as a systematic study, see [3, p. 85]: "It seems to be difficult to find a useful characterization of the symbols of those constant coefficient elliptic operators with this property."

Aim and results
The aim of this paper is a systematic investigation of the behaviour of fundamental solutions -and in particular whether or not they are positive close to the pole -for this class of uniformly elliptic operators of order 2m with constant coefficients. We find the above mentioned examples of sign changing fundamental solutions somehow unexpected because this means that even when applying a right hand side, which is concentrated at some point and points into one direction, the response of any solution to the differential equation will be sign changing and so -in some regions arbitrarily close to this point -in opposite direction to the right hand side. Indeed, we show in Theorem 2.3 that "positivity" is somehow the expected behaviour related to ellipticity.
In Sect. 2.3 we establish an inductive argument by space dimension. Roughly speaking this says that for understanding in any dimension whether one finds operators with sign changing fundamental solutions it suffices to understand the behaviour in "small" dimensions.
In Sect. 3 we calculate in almost closed form fundamental solutions for some specific fourth order elliptic operators in any dimension n ≥ 5. For a specific direction we even find a very simple closed expression. From n = 6 on, we observe sign change. While for n = 6, 7 we need a nonconvex symbol, for n ≥ 8 even convex symbols are admissible. The examples presented here use similar symbols as in [3] and [14], but they are constructed with the help of a different method based on the Fourier transform and the residue theorem. For n = 6 even the observation is new-to the best of our knowledge.
These examples yield the first important result. In order to answer the question asked by Davies for a systematic understanding of this phenomenon, we find in Sect. 4 explicit formulae for fundamental solutions from which it becomes clear which kind of elliptic symbols yield positive and sign changing fundamental solutions, respectively.
In odd dimensions we find the following general elegant formula. Theorem 1.2 Let n ≥ 2m + 1 be odd. Then, the fundamental solution K L is given by Here, if T is a j-multilinear form and v is a vector, we use the compact tensorial notation so, in particular, To avoid redundant parenthesis we write y |y| In even dimensions, due to the presence of the logarithm in (5), we are not able to achieve a comparable compact result, computations being much more involved. However, we show a related formula for the first "critical" dimension n = 2m + 2. Theorem 1.3 Let n = 2m + 2. Then, the fundamental solution K L is given by The proof is given in Sect. 4.2.
The difference between even and odd dimensions here reminds us somehow of the same distinction for the wave equation. In Theorem 1.3 (even dimensional) the integration is carried out over a one-codimensional surface with a weight function, which becomes infinite at its boundary. On the other hand, in Theorem 1.2 (odd dimensional) the integration is carried out over the boundary of this surface, i.e. a 2codimensional surface. The method of descent, as outlined in Sect. 2.3, gives further support to this observation.  Theorem 1.2 shows further which kind of symbol Q will e.g. in space dimension n = 2m + 3 yield K (0, e 1 ) < 0. For this one needs ∂ 2 1 1 Q (0, ξ ) > 0 for |ξ | = 1. This would follow e.g. from ∂ 1 Q(0, ξ ) = 0 and ∂ 2 1 Q(0, ξ ) < 0, i.e. (only) from a nonconvex shape of the level set of Q in these points.
For an example of this kind see [3, p. 100] and pp. 20-23 of a preliminary preprint version of this article which can be found at arXiv:1902.06503v1.
The situation is similar but not this clear in the even dimension n = 2m + 2, due to the relatively higher dimensional domain of integration and to the presence of further terms. However, thanks to the logarithmic singularity one may expect that also here, a nonconvex symbol may work. Indeed we prove in Sect. 5 the following result. In view of Sect. 2.3 and Theorem 2.3 we may immediately conclude the following general result.

Corollary 1.6
For any m ≥ 2 and n ≥ 2m + 2 there exists an elliptic symbol Q such that the fundamental solution of the associated operator L Q is sign changing for |x − y| → 0.
Together with Corollary 1.4 we have so obtained a complete picture of positivity and change of sign, respectively, in all "singular" dimensions n ≥ 2m.

Notation
We denote the partial derivative as D α or ∂ α or ∂ ∂α , where α is a multi-index, with the convention that if D 0 u = u for any function. Moreover, if j ∈ N, ∇ j u stands for the tensor of the j-th derivatives. Finally, we denote by H k the k-th dimensional Hausdorff measure.

Basic observations
In this section we consider only the case n > 2m of large dimensions. In what follows L always denotes a uniformly elliptic operator as in (1) with constant coefficients and only of highest order 2m. K L denotes John's fundamental solution as it is given in (4) and (5), respectively. By (2), without loss of generality, we may consider 0 as the pole of K L .

Homogeneity, decay and uniqueness of John's fundamental solution
In particular this yields for all x ∈ R n \ {0} and all multi-indices α ∈ N n Proof In case of odd n > 2m, (4) shows that K L can be obtained as (− ) (n+1−2m)/2 of a 1-homogeneous function. In case of even n ≥ 2m + 2, the proof of Theorem 1.3 in Sect. 4.2 shows that K L can be obtained (− ) (n+1−2m)/2 of a (−2)-homogeneous function. K L (0, x) = K L (0, −x) follows from the corresponding property of the symbol.
Proposition 2.2 Let K L andK L be two fundamental solutions for L which both obey (6). Then we find a solution of Lu = 0 in R n which thanks to elliptic regularity theory satisfies u ∈ C ∞ (R n ). To see this one may combine local elliptic L p -estimates (see [1,Theorem 15.1"]), the difference quotient method as outlined in [8,Section 7.11] and a bootstrapping argument. Since both K L andK L satisfy (6) we find that for any x ∈ R n \ {0} and any σ ∈ R \ {0}: Since n > 2m we conclude by continuity of u in 0 from letting σ → 0 that u(x) ≡ 0.

Ellipticity and positive directions
We prove the existence of "positive" directions (observe our sign convention for ellipticity) for the fundamental solutions which is somehow the simpler case and which one expects from the notion of ellipticity.

Theorem 2.3
Assume that n > 2m, L is a uniformly elliptic operator with constant coefficients of order 2m as introduced in (1) and consider the (2m − n)-homogeneous fundamental solution K L according to (4) and (5), respectively. Then there exists y ∈ R n \ {0} such that Proof We assume by contradiction that ∀y ∈ R n \ {0} : K L (0, y) ≤ 0.
Certainly, K L (0, y) ≡ 0. By continuity there exists a nonempty open set ⊂ S n−1 such that we have ∀y ∈ C : K L (0, y) < 0 on the corresponding cone We consider a fixed radially symmetric We introduce a corresponding solution (defined in the whole space) of the differential equation by Since for x ∈ B 1/2 (0) the intersection (x + C ) ∩ B 1/2 (0) is nonempty, U is strictly negative there. By compactness we find a constant C 0 > 0 such that Next, we introduce a scaling parameter σ ∈ (0, 1] and consider for the solution of the corresponding Dirichlet problem in B 1 (0) Here, denotes the corresponding Green function and its decomposition into fundamental solution and regular part. By continuous dependence on parameters and general elliptic theory (see [1]) we find that In what follows we consider only x ∈ B 1/2 (0). By the (2m − n)-homogeneity of the fundamental solution we obtain: provided that σ ∈ (0, 1] is chosen small enough. We fix such a suitable parameter and keep the corresponding u σ and ϕ σ fixed. We recall that we have shown: This yields (we recall that λ denotes the ellipticity constant of L) In the last step we used the elementary form of Gårding's inequality (see [6]) for operators, which have only constant coefficients and only of highest order, which follows from the ellipticity condition by employing the Fourier transform.
An alternative proof would follow from Corollary 1.4 and the inductive argument of Proposition 2.7 below.

An inductive argument
For simplicity we write in the remainder of this section In what follows we always assume that n − 1 > 2m, i.e. that n > 2m + 1.

A method of descent with respect to space dimension
Here we use the notation Let L n be an elliptic operator in R n as in (1) From this we obtain an operator in R n−1 by simply "forgetting" the x n -coordinate or by considering only functions, which do not depend on x n : with the corresponding symbol L n−1 is an elliptic operator. Next we define and aim at showing that this is John's (unique) fundamental solution for L n−1 .
We prove first that we have the expected homogeneity and hence also the expected decay at ∞.
Next we prove that K n−1 is in fact a fundamental solution for L n−1 .

Lemma 2.5 For all
and find that In order to proceed we need to overcome the difficulty thatφ / ∈ C ∞ 0 (R n ) by a suitable approximation. To this end we choose We find because we assume that n > 2m + 1. With this we conclude from (10) as claimed.
Combining Lemmas 2.4 and 2.5 with the uniqueness result of fundamental solutions with suitable degree of homogeneity from Proposition 2.2 we conclude: Proposition 2.6 K n−1 as defined in (9) is John's fundamental solution for L n−1 as it is given in (4) and (5), respectively.

Understanding "small" dimensions is sufficient
Proposition 2.7 Assume that n > 2m +1 and that for all elliptic operators L n−1 of the form (8) Proof Let L n be an arbitrary elliptic operator of the form (1) in R n with corresponding John's fundamental solution K n . We define L n−1 and K n−1 as in Sect. 2.3.1. Then Prop. 2.6 shows that K n−1 is John's corresponding fundamental solution. Making use of (9), the assumption yields the existence of a vector x ∈ R n−1 \ {0} such that This shows that there exists a point x n ∈ R which satisfies K n (x , x n ) > 0. Remark 1 Theorem 1.2 (i) and Proposition 2.7 yield a different proof of Theorem 2.3 by means of the inductive procedure.

Proposition 2.8
Assume that n > 2m + 1 and that there exists one elliptic operator L n−1 of the form (8) in R n−1 with John's corresponding fundamental solution K n−1 for which one finds a vector x ∈ R n−1 \ {0} such that K n−1 (x ) < 0. Then there exists one elliptic operator L n of the form (1) in R n with John's corresponding fundamental solution K n for which one finds a vector x = (x , x n ) ∈ R n \{0} such that K n (x ) < 0.
Proof Let L n−1 be an elliptic operator of the form (8) in R n−1 with symbol Q n−1 and corresponding John's fundamental solution K n−1 for which one finds a vector The operator L n is connected to L n−1 by the procedure described in Sect. 2.3.1. In particular John's fundamental solution K n−1 corresponding to L n−1 is given by (9).
The assumption yields the existence of a vector x ∈ R n−1 \ {0} such that This shows that there exists a point x n ∈ R which satisfies which completes the proof.

Sign changing fundamental solutions for m = 2
On S (R n ), the space of rapidly decreasing functions, one may define the Fourier- The inverse on So formally one would obtain the following expression for the corresponding fundamental solution: with δ 0 the delta-distribution in 0 and Fδ 0 = (2π) −n/2 , cf. (16) below. Due to the homogeneity of Q the integral in (12) is however not defined in L 1 (R n ) but at most as an oscillatory integral.
The Malgrange-Ehrenpreis Theorem, see [16,Theorem IX.23], states that a distributional solution F exists for L F = δ 0 , whenever L is a differential operator with constant coefficients. For elliptic operators the zero sets of Q are small in R n , which may allow one to give a classical meaning to (12) and gives a route to the fundamental solution. In some special cases this formula even allows one to derive an (almost) explicit fundamental solution. One such case is the following class of fourth order elliptic operators: where x = (x 1 , . . . , x n−1 ) and Although L is only interesting in the present setting whenever n ≥ 5, allow us to classify L for all dimensions. Proof To prove that ellipticity holds if and only if α > −2, is elementary. For |α| ≥ 2 one may split the symbol Q for L in (13) into real quadratic polynomials by: Whenever n = 2 and α ∈ (−2, 2] the operator L can be split into a product of two real second order elliptic operators following: This last splitting in dimensions n ≥ 3 with α ∈ (−2, 2), that is, replacing ξ 1 by ξ , would lead to a Fourier multiplier operator of order 2 with nonsmooth symbol, i.e. not even to a pseudodifferential operator.
The interesting case is hence α ∈ (−2, 2) and then it is convenient to use α = 2 cos γ with γ ∈ (0, π). So we write with corresponding symbol Q ξ , ξ n = ξ 4 +2 cos γ ξ 2 ξ 2 n +ξ 4 n . The fundamental solution for (14) is a regular distribution, so a function, which is C ∞ on R n \ {0} and moreover, homogeneous of degree 4 − n. We will recall that fact as the first step, when we prove the following result. Proposition 3.2 Let n ≥ 5 and γ ∈ (0, π). The fundamental solution F n,γ for (14) satisfies: • when x = 0 and x n = 0 one has • when x = 0 and x n = 0 one obtains Proof The proof is divided into 4 steps.
(i) Fundamental solution as distribution through an inverse Fourier transform. Let us first discuss the extensions of the Fourier transform F in (11). Both F and its inverse are well defined on S (R n ). By definition, a sequence (ϕ ) ∈N ⊂ S (R n ) converges to ϕ ∈ S (R n ) iff for all k ∈ N and for all multi-indices α ∈ N n 0 one has sup x∈R n (1 + |x|) k |D α (ϕ − ϕ)(x)| → 0 as → ∞. The natural extension to the space of tempered distributions S (R n ) is then [11, Definition 7.1.9] as follows: with a similar version for the inverse; ·, · denotes the duality between distribution and test function. (11) is well defined and one finds Fu ∈ L ∞ (R n ) and even the estimate Fu ∞ ≤ 1 (2π ) n/2 u 1 , but generically Fu / ∈ L 1 (R n ). In general The Fourier-transformation can also be extended to L 2 (R n ) by Plancherel and hence [11,Theorem 7.1.13] With the definition of the (inverse) Fourier transform in (16) The dot in (17) The distribution Q −1 on D (R n \ {0}), when extended to D (R n ), can only add a combination of the δ 0 -distribution and its distributional derivatives. Since in R n each such a distribution is homogeneous of degree −n or less, one finds that the extension is the regular distribution, that is, the function ξ → Q −1 (ξ ) on R n and we may skip the dot. By [11, Theorem 7.1.16] F n,γ with n ≥ 5 is then homogeneous of degree 4 − n and by [11,Theorem 7.1.18] one finds that F n,γ ∈ S (R n ) and that F n,γ | R n \{0} ∈ C ∞ (R n \ {0}) is a function. Also here the extension in 0 of this function can only add a combination of the δ 0 -distribution and its distributional derivatives and again, in R n each such distribution is homogeneous of degree −n or less. So indeed, one finds that also F n,γ is given by a function satisfying: With M = sup |ω|=1 F n,γ (ω) one finds from this (ii) Approximation as distribution through a summability kernel. Since we have established that F n,γ is a function, we will try next to derive a more explicit formula. Since Q −1 is not an L 1 -function the direct definition of the inverse Fourier-transform just after (11) is not applicable. We will use an approximation through a special positive summability kernel k ε , see [13, Section VI. 1.9]. A positive summability kernel on R n is defined as a family (k ε ) ε∈(0,ε 0 ] ⊂ C (R n ) with k ε ≥ 0 satisfying: The summability kernel that we use is a combination of a Gauss kernel in x ∈ R n−1 and 1 2ε e −|x n |/ε . We set and define Recall that So one finds (Fk ε ) (ξ ) = (Fk 1 ) (εξ ) and One obtains for all ϕ ∈ S (R n ), exploiting k ε * ϕ → ϕ in S (R n ), that and from the properties of distributions: In other words F −1 Fk ε 1 Q → F n,γ for ε ↓ 0 in the sense of distributions.
(iii) Approximation as a function through the summability kernel. Since F n,γ is a regular distribution and we may also write Setting we find from (18) and (20) that Since k 1 (x) ≤ ce −|x|/2 for some c > 0, we may estimate the last expression in (25) by splitting the corresponding integral in two parts: |x − y| < 1 2 |x|, implying |y| ≥ 1 2 |x|, and |x − y| > 1 2 |x|. Indeed, one finds for some C 1 > 0 that with σ n = S n−1 dω = 2 π n/2 (n/2) , the surface area of the unit sphere in R n . The result is that This allows us to use Lebesgue's dominated convergence theorem to find with (21) and (23) that with the last identity for any measurable ϕ : R n → R such that R n |ϕ (x)| |x| 4−n dx < ∞ .
(iv) An almost explicit formula by a contour integral. Next we will compute f ε using the formula F −1 Fk ε 1 Q from (22). The symbol Q = Q γ for (14) satisfies and the approximation f ε of the fundamental solution F n,γ (x) for (14) becomes Note that the integral converges near 0 for n ≥ 5. Near ∞ the integral converges for all n.
The integrand in (27) contains an analytic function of ξ n ∈ C and we find for γ ∈ (0, π) and x n > 0 by a contour integral in Since γ ∈ (0, π) holds and by the negative exponents in the exponential, the integral in (28) converges. The last term of (28) contains exp (−x n /ε), which shows that Hence by taking k > n − 4 it follows, whenever x n > 0, for ε ↓ 0 that It remains to consider which is well-defined for x n > 0. For x n < 0 one replaces x n by |x n | in (29). For the remaining integral in (29) we proceed by using for ξ ∈ R n−1 with n ≥ 4 the coordinates ξ = r cos ϕ r ω sin ϕ with r ≥ 0, ω ∈ S n−3 and ϕ ∈ [0, π] , where S n−3 = v ∈ R n−2 ; |v| = 1 . By the rotational symmetry in x ∈ R n−1 one may assume that x = x e 1 and we find through e i x ·ξ = e i|x |ξ1 = e ir|x | cos(ϕ) and dξ = (sin ϕ) n−3 r n−2 dϕdrdω , that Notice that We find for x = 0 and x n = 0 from (30) that which shows (15).
Together with Theorem 2.3, this proposition yields the proof of Theorem 1.1. Notice that whenever n ≥ 8, the fundamental solution changes sign even for γ < 1 2 π , where the level hypersurfaces of the symbol are still convex. , which is negative for γ ∈ 2 5 π, 4 5 π . To get an impression for which combinations of γ and β := arctan |x 8 | |x | the fundamental solution is negative, Fig. 1 contains a graph of which indeed only depends on β, γ ∈ 0, 1 2 π × (0, π).

Odd dimensions n > 2m: explicit fundamental solutions
In this section we prove Theorem 1.2 starting from John's formula (4). In other words, recalling that K L (x, y) = K L (0, x − y), we shall compute the iterated Laplacian for

The first iteration: 1F
Proposition 4.1 With the notation above, we have Before going into details of the proof, let us remark an important consequence of Proposition 4.1.

Remark 2
In the case n = 2m + 1, e.g. for a fourth order elliptic operator in R 5 , (4) and (35) imply which is thus a positive fundamental solution having the expected order of singularity. The only difference with the polyharmonic case is that an "angular dependent" positive factor appears. Notice that for the model polyharmonic case, as Q m (ξ ) := |ξ | 2m , then such factor is identically 1 and we of course retrieve its well-known fundamental solution.
We first recall a classical result about integrations of differential forms (see for instance [4,Satz 3]).

be a vector field. Then we have
with the convention that dx i means that this factor is missing.

Proof of Proposition 2
Step 1 Let us first compute ∇ F on points y = (y 1 , 0, . . . , 0) = 0. Here and everywhere in what follows we assume y 1 > 0. Then, by means of a rigid motion in R n , we will extend the result to the general case.
Writing y = r η, where |η| = 1, we have Noticing that the derivative in the first direction is indeed a normal derivative, therefore we infer Let us now compute all other first and second derivatives, which thus involve tangential directions. Without loss of generality, we may consider just ∂ 3 F and ∂ 2 33 F, the general case being similar. To this aim, we introduce the rotation matrix B ϕ which roughly speaking exchanges the first direction with the third: On the other hand, by definition Therefore, using Lemma 4.2, We observe that for the (n − 2) -form Hence, by (40)-(42) and Stokes' Theorem, noticing that ω = 0 on {ξ 1 = 0}, we infer Analogously one may compute ∂ k F(y 1 , 0, . . . , 0) for k = 1 and then obtain Indeed, it is sufficient to consider instead of B ϕ a similar matrix corresponding to a rotation in the plane y 1 , y k . Hence, (38) and (43) yield Step 2 Now we want to extend this identity to a generic point y ∈ R n \ {0}. To this aim, let us write y = |y|b 1 , where |b 1 | = 1, and complete this unit vector to a matrix The change of variable ξ = Bξ yields finally (34).
Step 3 Now it is the turn of second derivatives. We compute them with the same method we applied so far, so first we consider the easier case y = (y 1 , 0, . . . , 0) T with y 1 > 0 and then we extend this to a general y ∈ R n \ {0}.
Once again, we may for simplicity consider just ∂ 3 F, the other cases being similar as already mentioned. Let B ϕ be as in (39) and for y = (y 1 , 0, . . . , 0) Exactly as for F in (40), we have We denote R(z) := z 3 Q(z) for all z ∈ R n . Then we have having applied Stokes' Theorem. Letν denote the exterior normal of the half-sphere on its boundary, which gives its induced orientation. Hence, using Lemma 4.2 and the definition of R, Therefore, we may conclude that for all k ∈ {2, . . . , n} we have Recalling that ∂ 2 11 F(y 1 , 0, . . . , 0) = 0, we thus infer F(y 1 , 0, . . . , 0) = 1 |y| |ξ |=1 Step 4 Let now y ∈ R n \ {0} which we write as y = |y|b 1 , with |b 1 | = 1, and let us complete this unit vector to a matrix B := b 1 | . . . | b n ∈ SO(n). Moreover, recall F := F • B. Then, one has and therefore by (44),

The k-th iteration 1 k F and the proof of Theorem 1.2
The following theorem provides a general formula for the iterated Laplacian of F. We will see that tangential derivatives of the symbol play a fundamental role in the formula.
where d k,1 (n, m) = 1 and for j = 2, . . . , k: with (using the convention that the product is 1 whenever j = 1) Note that Theorem 1.2 is an immediate consequence of Theorem 4.3. This follows from putting k = n−2m+1 The rest of the subsection is devoted to proving Theorem 4.3, the strategy being the following. Firstly, we show that each term of the sum, namely once the Laplacian is applied, produces only terms of the same kind (so only even derivatives of 1 Q are involved) with order at most 2 j +2, each of them multiplied by the same suitable power of 1 |y| . This is achieved in Proposition 4.4. As a consequence, we obtain some recurrence formulae for the coefficients d k, j in the proof of Theorem 4.3. These relations will be important to finally prove the theorem by induction. Let us fix k ∈ N and j ∈ {0, . . . , k − 1} and define J k, j (y 1 , 0, . . . , 0) where Therefore, we obtain J k, j (y) = c (n, m, j)J k+1, j−1 + c 0 (n, m, k, j)J k+1, j + J k+1, j+1 .
Proof In order to simplify the notation, as k, j are fixed, we write J (y) instead of J k, j (y).
Step 1 Let y = (y 1 , 0, . . . , 0), y 1 > 0. First of all, writing J in polar coordinates and We may rewrite the argument as with the shorter notation Since the terms which remain must have only cosines, everything vanishes except for the term with h = 0 in the first sum and the one with h = 1 in the third. Therefore, Moreover, we compute on Inserting in (54), we infer Of course, an analogous formula holds for any h ∈ {2, . . . , n}, namely Step 2 Let us now consider y ∈ R n \ {0}, so y = |y|b 1 , where |b 1 | = 1, and let J (B(|y|, 0, . . . , 0) ∇J (|y|, 0, . . . , 0) Returning therefore to the variable ξ = Bξ , we get that is, (48).
Step 3 Let us again consider y = (y 1 , 0, . . . , 0), y 1 > 0, and compute J (y). DefiningH (ϕ) := ∂ y 3 J (B ϕ y), according to the splitting in (56), we havẽ Let us differentiate with respect to ϕ term by term. Firstly, Therefore we obtain As in Step 1, the terms which remain are the first with h = 0 and the third with h = 1, so Differentiating the first term as in (55) we obtain hence Let us now address to the second term in (57): Hence, differentiating with respect to ϕ, with similar computations as forH 0 , we obtain: With similar computations as in (55), we infer Finally, we have to consider the third term in (57): Hence, The first two terms vanish for any choice of h, while the last one remains only for Hence, recalling the splitting (57), by (58)-(60) we obtain (omitting from now on in each integral its differential dH n−2 (ξ )): Therefore, the same being valid for any variable y h with h ∈ {2, . . . , n}, and recalling (53) if h = 1, we may compute the Laplacian of J : Here, we denote By homogeneity of the symbol, one has (see Lemma 4.5 below) Moreover, in order to handle the term ∂ 2 j−2 1 2 j−2 1 Q , we may apply the following wellknown identity with u = ∂ Noticing that the normal derivatives in (63) may be handled as in (62), it remains the term with the tangential part of the Laplacian. However, it vanishes when integrated on S. Hence, Inserting (62) and (64) in (61) and summing the constants, we finally end up with (49) and thus with our formula (52).

Lemma 4.5
Let Q(·) be positive and p-homogeneous. Then, one has for any multiindex α ∈ N n 0 and any x ∈ R n \ {0}: Proof By assumption we have for r > 0 that Differentiation with respect to x yields: Differentiating now with respect to r gives: The claim follows by putting r = 1.

Proof of Theorem 4.3
Notice that for k = 1, F is already in the form (45) by Proposition 4.1. We thus proceed by induction and let us suppose that k F has the form (45) for some k ∈ N, k > 1, with coefficients as in (46)-(47), namely Applying the recursive formula (52), we thus have for j = 3, . . . , k + 1. We have d k,k+1 according to (47) and put d k,0 := d k,−1 := 0. Hence, for j = 1 and j = 2 we have the recurrence relations and for those the formulae (46) and (47) are easily checked. We show the shape of d k+1, j (n, m) for j = 3, . . . , k + 1, the cases j = 1, 2 being analogous but simpler. In order to prove that d k+1, j (n, m) has the shape as in (45) with k + 1, which will prove Theorem 4.3, we show that the following term is equal to 0: This concludes the proof of Theorem 4.3.

The even dimension n = 2m + 2: explicit fundamental solutions
Here we provide the proof of Theorem 1.3. The starting point is John's formula (5), according to which we have For n = 2m we immediately see that the fundamental solution of L satisfies which is similar to the case n = 2m + 1 above. Hence, in order to obtain an explicit expression for K L , we thus have to compute the Laplacian of Due to the logarithmic term, the calculations for even dimensions cannot be simplified similarly to the previous section. An application of Stokes' theorem would change log(ξ 1 ) into 1 ξ 1 , a non-integrable singularity. For this reason we restrict ourselves to the case n = 2m + 2.
The strategy is similar to the one applied in the proofs of Propositions 4.1 and 4.4.
Step 1 Let y = (y 1 , 0, . . . , 0) with y 1 > 0. Splitting G as we first infer as Q(ξ ) = Q(−ξ) by homogeneity of the symbol. Moreover, In order to compute the other first derivatives of G, let us split it as dH n−1 (ξ ) =: log |y|G 1 (y) +G 2 (y). Hence, We use the matrix B ϕ defined in (39) and definẽ Concerning the first term, exploiting Q(ξ ) = Q(−ξ), we find for any ỹ Concerning the second term, reasoning as in (40), we get Therefore, by (69)-(71) we conclude and, analogously, for any k ∈ {2, . . . , n}, one has Because this formula is consistent with (67), we can write it in a more compact way as Step 2 Let now y ∈ R n \ {0} so that y = |y|b 1 with |b 1 | = 1, and let B = b 1 | · · · | b n ∈ SO(n). Defining we have Therefore, Observing that (B T ξ) 1 = y·ξ |y| and multiplying by B T shows that (72) holds also for any y ∈ R n \ {0}.

A change of variables and (39) imply
As the same holds for any k ∈ {2, . . . , n}, we infer This together with (68) yields (omitting from now on the differentials) Step 4 Let now y ∈ R n \ {0}. Reasoning as in Step 2 and using the same notation as there we obtain Recalling that B · e 1 = y |y| , (73) implies Finally, due to the 2m-homogeneity of Q by means of Lemma 4.5 we have Hence the previous formula simplifies to and, recalling (2), the proof is concluded.

Further examples: sign change of suitable fundamental solutions for n ≥ 2m + 2
Proposition 2.8 shows that in any space dimension n ≥ 2m + 2 there exists an elliptic operator L of the form (1) in R n with corresponding fundamental solution K L , where one finds a vector y ∈ R n \ {0} such that K L (0, y) < 0, provided we are able to construct such an operator of order 2m in space dimension n = 2m + 2. Together with Theorem 2.3 this will show that K L is sign changing near the origin, i.e. near its singularity. In view of Theorem 2.3 this will prove Theorem 1.1. Fig. 2 The shape of sublevel sets of Q α for α < 0 (left) and α > 0 (right)
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.