Time-harmonic solutions for Maxwell's equations in anisotropic media and Bochner-Riesz estimates with negative index for non-elliptic surfaces

We solve time-harmonic Maxwell's equations in anisotropic, spatially homogeneous media in intersections of $L^p$-spaces. The material laws are time-independent. The analysis requires Fourier restriction-extension estimates for perturbations of Fresnel's wave surface. This surface can be decomposed into finitely many components of the following three types: smooth surfaces with non-vanishing Gaussian curvature, smooth surfaces with Gaussian curvature vanishing along one-dimensional submanifolds, but without flat points, and surfaces with conical singularities. Our estimates are based on new Bochner-Riesz estimates with negative index for non-elliptic surfaces.


Introduction
The purpose of this article is to prove the existence of solutions to the timeharmonic Maxwell's equations and estimating the solutions (electromagnetic fields) in terms of the input data (currents) in L p -spaces. Let (E, H) : R × R 3 → R 3 × R 3 denote the electric and magnetic field, (D, B) : R × R 3 → R 3 × R 3 the displacement field and magnetic induction, and (J e , J m ) : R × R 3 → R 3 × R 3 the electric and magnetic current. Maxwell's equations in the absence of charges are given by (1) ∂ t D = ∇ × H − J e , ∇ · D = ∇ · B = ∇ · J e = ∇ · J m = 0, We suppose that displacement and magnetic field are related with electric field and magnetic induction through time-independent and spatially homogeneous material laws. This leads to supplementing (1) with (2) D(t, x) = εE(t, x), B(t, x) = µH(t, x), ε ∈ R 3×3 , µ ∈ R 3×3 .
ε is referred to as permittivity, and µ is referred to as permeability. Permittivity and permeability are positive-definite in classical physical applications. We suppose in the following that ε and µ are diagonal matrices and write Maxwell's equations are invariant under change of basis, i.e., the transformations X (t, x) = M X(t, M t x) for the involved vector fields with M ∈ SO(3), and timeparity symmetry (t, x) → (−t, −x). Hence, the more general case when ε and µ are commuting positive-definite matrices, or equivalently: simultaneously orthogonally diagonalizable, reduces to (3). For physical explanations, we refer to [15,34]. The assumption ∇ · D = 0 corresponds to the absence of electrical charges and ∇ · B = 0 translates to the absence of magnetic monopoles. Due to conservation of charges, the currents are likewise divergence-free. Since magnetic monopoles are hypothetical, J m is vanishing for most applications. Here, we consider the more general case, which will highlight symmetry between E and H. In this paper we focus on the fully anisotropic case Upon considering the time-harmonic, monochromatic ansatz With (2) we arrive at the equations Below W m,p (R d ) denotes the L p -based Sobolev space defined by We prove the following: Theorem 1.1. Let 1 ≤ p 1 , p 2 , q ≤ ∞, ε, µ ∈ R 3 as in (3), (4) and and 0 ≤ 1 then, for any given ω ∈ R \ {0}, there exists a distributional time-harmonic solution to fully anisotropic Maxwell's equations (6) that satisfies with locally uniform dependence with respect to ω ∈ R\{0}. If additionally J e , J m ∈ L q (R 3 ), q < ∞, then E, H ∈ W 1,q (R 3 ) is a weak solution satisfying (E, H) W 1,q (R 3 ) p,q,ω (J e , J m ) L p 1 (R 3 )∩L q (R 3 ) .
We shall see that the Fourier multiplier derived by inverting (6) for ω ∈ R is not well-defined in the sense of distributions. A common regularization is to consider ω ∈ C\R and derive estimates independent of dist(ω, R). This program was carried out in our previous works [11,39], which were concerned with isotropic, possibly inhomogeneous, respectively, partially anisotropic, but homogeneous media. The necessity of considering (J e , J m ) within intersections of L p -spaces and the connection with resolvent estimates for the Half-Laplacian was discussed in [39]. In the present work we need to regularize differently due to a more complicated behaviour of the involved Fourier symbols with respect to the change ω → ω + iε. In other words, we do not prove a Limiting Absorption Principle in the classical sense.
In the proof we will reduce the analysis to the case µ 1 = µ 2 = µ 3 = 1 as in [36] in order to simplify the notation. We will justify this step in Section 3. In the partially anisotropic case ε 1 = ε 2 = ε 3 the matrix-valued Fourier multiplier associated with Maxwell's equations can be diagonalized easily and a combination of Riesz transform estimates and resolvent estimates for the Half-Laplacian are used to prove uniform bounds. In our fully anisotropic case (4) this does not work at all. Instead of diagonalizing the symbol, we take the more direct approach of inverting the matrix Fourier multiplier associated with (6). Taking the Fourier transform in R 3 , denoting with ξ ∈ R 3 the dual variable of x ∈ R 3 , and the vector-valued Fourier transform of E withÊ, likewise for the other vector-valued quantities, we find that (6) is equivalent to (9) ib(ξ)Ê(ξ) + iωµĤ(ξ) =Ĵ m , ξ ·Ĵ m = ξ ·Ĵ e = 0, ib(ξ)Ĥ(ξ) − iωεÊ(ξ) =Ĵ e .
From Proposition 1.2 we infer that solutions to anisotropic Maxwell's equations can be found provided that the mapping properties of the Fourier multiplier with symbol p −1 (ω, ξ) or, actually, an adequate regularization of this, can be controlled. The first step of this analysis is to develop a sound understanding of the geometry of S := {ξ ∈ R 3 : p(ω, ξ) = 0}, with an emphasis on its principal curvatures. This has essentially been carried out by Darboux [12] and Liess [36,Appendix]. We devote Section 3 to recapitulate these facts along with some computational details that were omitted in [36]. S is known as Fresnel's wave surface, which was previously described, e.g., in [12,36,32,16]. We refer to Figure 2 for visualizations. Despite its seemingly complicated structure, this surface can be perceived as nonsmooth deformation of the doubly covered sphere in R 3 . For the involved algebraic computations we provide a MAPLE TM sheet for verification.
We turn to a discussion of the regularization of p(ω, ξ) −1 . Motivated by Cramer's rule, we multiply (10) with the adjugate matrices and divide by p(ω, ξ) + iδ. This leads us to approximate solutions (E δ , H δ ). We postpone the precise definition to Section 2. The main part of the proof of Theorem 1.1 is then to show uniform bounds in δ = 0: for q, p 1 , p 2 as in Theorem 1.1. In Section 2 we shall see how this allows us to infer the existence of distributional solutions to (6) and how the limits can be understood as principal value distribution and delta distribution for Fresnel's wave surface in Fourier space. Moreover, the distributional solutions are weak solutions provided that the currents have sufficiently high integrability.
We point out the connection to Bochner-Riesz operators of negative index and seemingly digress for a moment to explain key points for these operators. For 0 < α < 1, consider the Bochner-Riesz operator with negative index given by C d denotes a dimensional constant, Γ denotes the Gamma function, and x + = max(x, 0). For 1 ≤ α ≤ (d + 1)/2, S α is explained by analytic continuation. The body of literature concerned with Bochner-Riesz estimates with negative index is huge, see, e.g., [40,5,24,10,33]. In Section 4 we give a more exhaustive overview. For α = 1, we find because the distribution in (13) for α = 1 coincides with the delta distribution up to a factor. Estimates for such Fourier restriction-extension operators are the backbone of the Limiting Absorption Principle for the Helmholtz equation (cf. [25]). It turns out that we need more general Fourier restriction-extension estimates than the ones associated with elliptic surfaces because the Gaussian curvature of the Fresnel surface S changes sign, as we shall see in Section 3. We take the opportunity to prove estimates for generalized Bochner-Riesz operators of negative index for non-elliptic surfaces as the associated Fourier restriction-extension operators will be important in the proof of Theorem 1.1.
To describe our results in this direction, let d ≥ 3 and S = {(ξ , ψ(ξ )) : ξ ∈ [−1, 1] d−1 } be a smooth surface with k ∈ {1, . . . , d−1} principal curvatures bounded from below. The case d = 2 was disclosed by Bak [2] and Gutiérrez [24]. Let In the following theorem, we show L p -L q -bounds within a pentagonal region (see Figure 1) For two points X, Y ∈ [0, 1] 2 , let At its inner endpoints B α,k , B α,k , we show restricted weak bounds and on part of its boundary, we show weak bounds (14) holds true for Furthermore, we find estimates (17) to hold for For any α the constant in (14)-(18) depends on the lower bounds of the principal curvatures and χ C N and ψ C N for N = N (p, q, d). In particular it is stable under smooth perturbations of χ and ψ.
The proof is based on the decay of the Fourier transform of the surface measure on S (cf. [37], [41,Section VIII.5.8]) and convenient decompositions of the distribution . We also show that the strong bounds are sharp for α ≥ 1 2 . In the elliptic case the currently best results were shown by Kwon-Lee [33,Section 2.6]. This also shows that our strong bounds are not sharp for α < 1 2 . We refer to Section 4 for further discussion. To describe the remainder of our analysis, we recall important properties of the Fresnel surface. Up to arbitrary neighbourhoods of four singular points, the surface is a smooth compact manifold with two connected components. The Gaussian curvature vanishes precisely along the so-called Hamiltonian circles on the outer sheet. However, the surface is never flat, i.e., there is always a principal section away from zero. Around the singular points, the surface looks conical and ceases to be a smooth manifold.
We briefly explain how this leads to an analysis of the Fourier multiplier (p(ω, ξ)+ iδ) −1 , ω ∈ R\{0}, 0 < |δ| 1. We recall that solutions to time-harmonic Maxwell's equations are constructed by considering δ → 0 with bounds independent of δ. The non-resonant contribution of {ξ ∈ R 3 : |p(ω, ξ)| ≥ t 0 }, t 0 > 0 away from Fresnel's wave surface is estimated by Mikhlin's theorem and standard estimates for Bessel potentials. This high-frequency part of the solutions is responsible for the condition 0 ≤ 1 p2 − 1 q ≤ 1 3 in (7). This contribution was called global in [39]. We refer to [39,Section 3] for further explanation how this contribution does not allow for an . After smoothly cutting away the global contribution, we focus on estimates for the multiplier (p(ω, ξ) + iδ) −1 in a neighbourhood {|p(ω, ξ)| ≤ t 0 } near the surface. It turns out that around the smooth elliptic part with Gaussian curvature bounded away from zero, we can use the estimates for the Bochner-Riesz operator from Theorem 1.3 for d = 3, k = 2, α = 1. However, there is also a smooth non-elliptic part where the modulus of the Gaussian curvature is small and vanishes precisely along the Hamiltonian circles. Here, Theorem 1.3 applies for d = 3, k = 1, α = 1. In the corresponding analyis of the multiplier (p(ω, ξ) + iδ) −1 we foliate the neighbourhoods of the Fresnel surface by level sets of p(ω, ξ). The contributions of the single layers are estimated with the Fourier restriction-extension theorem. In the analysis we use decompositions in Fourier space generalizing arguments of Kwon-Lee [33,Section 4], where the decompositions were adapted to the sphere.
For the contribution coming from neighbourhoods of the four isolated conical singularities, we will apply Theorem 1.3 once more for d = 3, k = 1, α = 1. On a technical level, a major difference compared to the other regions comes from the fact that the cone is not a smooth manifold: we use an additional Littlewood-Paley decomposition and scaling to uncover its mapping properties. Jeong-Kwon-Lee [30] previously applied related arguments to analyze Sobolev inequalities for second degree non-elliptic operators.
We further mention the very recent preprint by Castéras-Földes [9] (see also [4]). In [9] L p -mapping properties of Fourier multipliers (Q(ξ) + iε) −1 for fourth order polynomials Q were analyzed in the context of traveling waves for nonlinear equations. The analysis in [9] does not cover surfaces {Q(ξ) = 0} containing singular points, and the L p -L q -boundedness range stated in [9,Theorem 3.3] is strictly smaller than in the corresponding results given in Theorem 1.3.
Outline of the paper. In Section 2 we carry out reductions to prove Theorem 1.1. We anticipate the uniform estimates of the regularized solutions that we will prove in Sections 5 -6, by which we finish the proof of Theorem 1.1. In Section 3 we recall the relevant geometric properties of the Fresnel surface and reduce our analysis to the case ω = µ 1 = µ 2 = µ 3 = 1. In Section 4 we recall results on Bochner-Riesz estimates with negative index for elliptic surface and extend those to estimates for a class of more general nondegenerate surfaces by proving Theorem 1.3. In Section 5 we use these estimates to uniformly bound solutions to (5) corresponding to the smooth part of the Fresnel surface. In Section 6 we finally estimate the contribution with Fourier support close to the four singular points.

Reduction to multiplier estimates related to the Fresnel surface
The purpose of this section is to carry out the reductions indicated in the Introduction. We first define suitable approximate solutions (E δ , H δ ) and present estimates for those related to the different parts of the Fresnel surface and away from the Fresnel surface. With these estimates at hand, to be shown in the upcoming sections, we finish the proof of Theorem 1.1. At the end of the section we give explicit formulae for the solution.
We work with the following convention for the Fourier transform: For f ∈ S(R d ) the Fourier transform is defined bŷ and as usually extended by duality to S (R d ). The Fourier inversion formula reads for f ∈ S(R d ) 2.1. Approximate solutions. By Proposition 1.2 the original anisotropic Maxwell system leads to the following second order 3 × 3-system for E and H (19) ( for the polynomials q 0 , q 1 as defined there. Inverting M E (ξ) − ω 2 using Cramer's rule, we find for all ξ ∈ R 3 such that p(ω, ξ) = 0: Here, adj(M ) denotes the adjugate matrix of M . Sarrus's rule and lengthy computations yield that the components of Z = Z ε,µ are given as follows: A crucial observation is that the associated matrix-valued Fourier multiplier will be applied to divergence-free functions. This is a consequence of (19) and (20). For that reason the fourth order terms in the entries can be ignored (if convenient), which becomes important when estimating the large frequency parts of our approximate solutions. Let Z eff (ξ) = Z eff ε,µ (ξ) denote the unique matrix-valued polynomial of degree 2 such that with v · ξ = 0. In view of (19) and (20) it is natural to define the approximate solutions (E δ , H δ ) for |δ| = 0 as follows: e (ξ) + ωĴ m (ξ)). To prove Theorem 1.1, we show estimates for these functions that are uniform with respect to δ. The global part away from the Fresnel surface is considered in the next subsection, the remaining estimates will be done later. Then, taking these estimates for granted, we show how to conclude the argument.

2.2.
Local and global contributions. We turn to the description of the different contributions of (E δ , H δ ). We split the local and global contribution. Let β 1 , β 2 ∈ C ∞ (R 3 ) satisfy β 1 (ξ) + β 2 (ξ) = 1 with where t 0 > 0 denotes a small constant. t 0 will be chosen later when carrying out the estimates close to the surface. Also, for m ∈ C ∞ (R d ) we write Proposition 2.1. Let E δ , H δ be given by (22). Then, we find the following estimate to hold uniformly in |δ| > 0: with |p(ω, ξ)| ≥ c > 0 on supp(χ) and χ(ξ) = 1 on supp(β 2 ). We first consider the case q = ∞. Then By the choice of χ we have the following uniform estimates with respect to δ: Since 1 < q < ∞, Mikhlin's theorem (cf. [19,Chapter 6]) applies and Bessel potential estimates (see for instance [11,Theorem 30]) yield for the claimed range of exponents. If q = ∞, we first use Sobolev embedding to lower q < ∞, and applying the previous argument gives (23) for 0 < 1 p < 1 3 , which is all we had to show in this case. This gives the claim concerning L q (R 3 )-integrability. For the Sobolev regularity we obtain in a similar fashion The proof is complete.
The paper is mainly devoted to estimate the local contribution close to the Fresnel surface S = {p(ω, ξ) = 0}. In Section 3 we shall see that the Fresnel surface has components of the following type: • smooth components with non-vanishing Gaussian curvature, • smooth components with curvature vanishing along a one-dimensional submanifold (Hamiltonian circles), but without flat points, • neighbourhoods of conical singularities. This fact is established in Corollary 3.8. Precisely, it suffices to consider six components of the first kind, and four components of the second and third type.
Corresponding to the three types listed above, we split with smooth compactly supported functions localizing to neighbourhoods of the components of the above types. The estimate for the smooth components with non-vanishing Gaussian curvature is a consequence of estimates for Bochner-Riesz operators with negative index that we will prove in Section 4: Proposition 2.2. Let 1 ≤ p, q ≤ ∞ and (E δ , H δ ) as in (22). We find the following estimate to hold uniformly in |δ| = 0: By similar means, we show the inferior estimate for components with vanishing Gaussian curvature along the Hamiltonian circles: (22). We find the following estimate to hold uniformly in |δ| = 0: At last, the estimate around the singular points is shown in Section 6: (22). We find the following estimates to hold uniformly in |δ| = 0: For these estimates, due to bounded frequencies, the precise form of Z ε,µ (or Z eff ε,µ ) is not important. It suffices to show the above estimates for the multiplier Again due to bounded frequencies, the as a consquence of Young's inequality.

2.3.
Proof of Theorem 1.1. By Propositions 2.1 -2.4 we have uniform bounds in δ = 0: for q, p 1 , p 2 as in the assumptions of Theorem 1.1. Hence, there is a weak limit (E, H) ∈ L q (R 3 ; C 6 ), which satisfies the same bound by the Banach-Alaoglu-Bourbaki theorem. We have to show that the approximate solutions weakly converge to distributional solutions of Indeed, (22) gives

From (21) one infers after lengthy computations
As a consequence we obtain By Proposition 2.1 -2.4, and Remark 2.5 we have and, when assuming J m ∈ L q (R 3 ), so that the only δ-dependent term vanishes as δ → 0. This implies in the distributional sense and even in the weak sense for J m ∈ L q (R 3 ). Similarly, one proves the validity of the second equation in (24), and the proof is complete.

Explicit representations of solutions.
At last, we give explicit representations of the constructed solutions. By Sokhotsky's formula (cf. Sections 3.2 and 6.1 in [26]): Furthermore, in the sense of distributions.
In the context of the easier Helmholtz equation this allows to write for so-called outgoing solutions Proposition 2.6 suggests that the solutions to anisotropic Maxwell's equations can again be written as principal value and delta distribution in Fourier space. However, Proposition 2.6 only allows to make sense of the principal value and delta distribution if S = {ξ ∈ R 3 : p(ω, ξ) = 0} is a smooth manifold. But there are four isolated singular points ζ 1 , . . . , ζ 4 ∈ S as we will prove in Proposition 3.2. Still, we shall see how v.p. 1 p(ω,ξ) and δ S (ξ) can be understood as Fourier multipliers with certain L p -mapping properties. For a dense set, e.g., J ∈ S(R 3 ), ζ i / ∈ supp(Ĵ), we can explain δ S as a Fourier multiplier The density follows by Littlewood-Paley theory. As a consequence of Sections 5 and 6, we have S e ix.ξĴ (ξ)dσ(ξ) for p and q as in Proposition 2.4 with a bound independent of the support ofĴ. This allows to extend . This explains the formula for solutions to anisotropic Maxwell's equations. Notice that in these formulae we may replace the matrices Z ε,µ (ξ), Z µ,ε (ξ) by the corresponding effective matrices.

Properties of the Fresnel surface
As explained above, the set {ξ ∈ R 3 : p(ω, ξ) = 0} plays a decisive role for our analysis. This classical quartic surface is known as Fresnel's surface initially discovered by Augustin-Jean Fresnel in 1822 to describe the phenomenon of double refraction. In an optically anisotropic medium, e.g., a biaxial crystal, Fresnel's surface corresponds to Huygen's elementary spherical wave surfaces in isotropic media. This surface was already studied in the 19th century by Darboux [12]. For an account on classical references we refer to the survey by Knörrer [32]. In the present context the curvature properties will be most important, which were collected by Liess [36,Appendix]. We think it is worthwhile to elaborate on Liess's presentation, as we shall also discuss first and second fundamental form in suitable coordinates.
We recall the key properties of Fresnel's wave surface and Recall that we assume full anisotropy (4). We first notice that we can reduce our analysis to the case µ 1 = µ 2 = µ 3 = ω = 1. This results from the change of coordinates ξ → η given by Notice that this change of coordinates results from a suitable dilation of the coordinates, which corresponds to an appropriate dilation in physical space. To see the equivalence, let us introduce the corresponding quantities for Then one observes ω 4 N (η) = p(ω, ξ), hence the qualitative properties of Fresnel's surface in the special case µ 1 = µ 2 = µ 3 = ω = 1 carry over to the general case.
For this reason we focus on the analysis of . Notice that (4) then reads In the following we write We first show that S * is a smooth manifold away from four singular points. To see this, we compute The set of singular points is denoted by Σ.
The reason for this definition is that S * \ Σ is a smooth manifold, whereas the neighbourhood of the singular points require a separate analysis. It turns out that there are precisely four singular points. This is a consequence of the following result.
In particular, the Gaussian curvature is well-defined and smooth on S * \ Σ, i.e., away from the four singular points. We now introduce the explicit parametrization of S * by Darboux and Liess ([36,A3]). Our parameters (s, t) correspond to (β, α ) in Liess' work. As in [36], this parametrization is given away from the four singular points and the principal sections S ∩ {η i = 0}.
Proof. If we define η := Φ(s, t) ∈ R 3 , then one can subsequently verify for all s, t such that the argument of the square root is positive. On the other hand, every point of ( can be written in this way. To see this, one solves the linear system In this way one finds is well-defined (the arguments of all square roots are positive) if and only if either ε j < s < ε j+1 < t < ε j+2 or ε j < t < ε j+1 < s < ε j+2 holds provided that ε j < ε j+1 < ε j+2 .
We note that the two parameter regions ε j < s < ε j+1 < t < ε j+2 and ε j < t < ε j+1 < s < ε j+2 give rise to the inner, respectively, outer sheet of the wave surface, cf. Figure 2. Both sheets meet at the singular points that formally correspond to s = t = ε j+1 where, in accordance with Proposition 3.2, one has η j+1 = Φ j+1 (s, t) = 0. We now turn towards the computation of the Gaussian curvature on S * \Σ. This will first be done away from the principal sections, but the formula will prevail also in the principal sections since S is a smooth manifold in that region as we showed above. We start with computing the relevant derivatives for the first and second fundamental form of S * : From these formulae one gets the following.
Proposition 3.4. The first fundamental form of S * \ Σ is given by Proof. This follows from lengthy, but straightforward computations based on To write down the second fundamental form, we introduce the following functions: Proposition 3.5. The second fundamental form of S * \ Σ is given by Proof. By definition, the functions L, M, N are given by where ν(s, t) denotes the outer unit normal on S \ Σ at the point Φ(s, t). In Euclidean coordinates, a normal at η = Φ(s, t) is given by and obtain after normalization Using this formula for the unit normal field ν, and plugging in the formulae for Φ ss , Φ st , Φ tt , one obtains the above expressions for L(s, t), M (s, t), N (s, t).
We continue with the formulae for the Gaussian and mean curvature, which were given in (A.1),(A.2) in Liess' work [36].
Proof. The determinant of the first fundamental form is given by The determinant of the second fundamental form is .
So the Gaussian curvature at the point Φ(s, t) is Following Liess, we define α(s, t) to be the squared distance of the origin to the tangent plane through Φ(s, t) ∈ S * \ Σ. Then From this we deduce Proposition 3.7. The mean curvature at Φ(s, t) ∈ S * \ Σ is given by (α = α(s, t)) Proof. This is a consequence of the formula We remark that our result deviates by the factor − 1 2 from Liess' formula [36, (A.2), p. 91]. This does not change the curvature properties, which we describe in the following: The Gaussian curvature K vanishes precisely in those points where α(s, t) attains one of the values ε 1 , ε 2 , ε 3 . We assume ε 1 < ε 2 < ε 3 for simplicity. Then one has ε 1 < α(s, t) < ε 3 so that the Gaussian curvature vanishes precisely at those points where α(s, t) = ε 2 . Those are given by t = T (s) where This is the parametrization of a one-dimensional submanifold that is called a Hamiltonian circle. Notice that each of the four singular point has its own Hamiltonian circle. (They are distinguished by σ i , σ i+2 ∈ {−1, +1} in Proposition 3.2). By Proposition 3.7 the mean curvature is non-zero along the Hamiltonian circles. We thus conclude that in the smooth regular part of Fresnel's wave surface, there is at least one principal curvature bounded away from zero. The Gaussian curvature is To summarize the geometric properties, we can perceive S as union of two sheets A and B, linked together at the singular points, when A is completely encased by B. A is convex, but B is not. Close to the singular points, B is not convex, and the Gaussian curvature is negative. Increasing geodesic distance from the singular points on B, we reach the Hamiltonian circles: the curvature vanishes. Beyond the Hamiltonian circles, B is locally convex, too, and has again positive Gaussian curvature.
Corollary 3.8. The wave surface S = {ξ ∈ R 3 : p(ω, ξ) = 0} admits a decomposition S = S 1 ∪ S 2 ∪ S 3 , where (i) S 1 is a compact smooth regular manifold with two non-vanishing principal curvatures in the interior, (ii) S 2 is a compact smooth regular manifold with one non-vanishing principal curvature in the interior, (iii) S 3 is the union of (small) neighbourhoods of the singular points described in Proposition 3.2.
For later sections, it will be important to have these curvature properties likewise for level sets {p(ω, ξ) = t} t∈[−t0,t0] for some 0 < t 0 1 with uniform bounds in t. For this purpose, recall that for an implicitly defined surface {ξ ∈ R 3 : F (ξ) = 0} the Gaussian curvature is given by (cf. [ Hence, again due to smoothness of p and |∇F | ≥ d > 0, along the curves on the level sets, where the Gaussian curvature vanishes, we have one principal curvature bounded from below. Choosing the level sets close to the original surface, we find one principal curvature bounded from below likewise on all the layers.

Generalized Bochner-Riesz estimates with negative index
The purpose of this section is to show Theorem 1.3. In the following let d ≥ 2 and S = {(ξ , ψ(ξ )) : ξ ∈ [−1, 1] d−1 } be a smooth surface with k ∈ {1, . . . , d − 1} principal curvatures bounded from below. Let We show strong estimates for a range of p and q with weak endpoint estimates as stated in Theorem 1.3. We start with recapitulating Bochner-Riesz estimates in the elliptic case, which is understood best.
4.1. Bochner-Riesz estimates with negative index for elliptic surfaces.
If ψ is elliptic, i.e.,the Hessian ∂ 2 ψ has eigenvalues of a fixed sign on [−1, 1] d−1 , then T α is a Bochner-Riesz operator of negative index. As explained above, we shall show bounds also for possibly degenerate ψ, which will be useful in the next sections. For solutions to time-harmonic Maxwell's equations we are interested in the case d = 3, α = 1, corresponding to restriction-extension operator: We take a more general point of view to show that the considerations in the next section also apply in higher dimensions and general α. To put our results into context, we digress for a moment and recapitulate results on the classical Bochner-Riesz problem.
For α > 0 recall The Bochner-Riesz conjecture (for elliptic surfaces) with negative index states: The necessity of these conditions was proved by Börjeson [5]. We refer to [33, Section 2.6] for a survey, where the currently widest range is covered. In the special case α > 1 2 contributions are due to Bak-McMichael-Oberlin [1, Theorem 3] and Gutiérrez [24,Theorem 1], see also [40,2]. In [10, Remark 2.3] was also pointed out that T α : 2 . In the following we recall arguments from [10,33], which were needed for the proofs and will be used in the next section for more general surfaces. In the first step we decompose the multiplier distribution.
For α > 0, let D α ∈ S (R d ) be defined by which is again extended by analytic continuation to the range 1 ≤ α < d+1 2 . We recall the following lemma to decompose the Fourier multiplier: . For α > 0, there is a smooth function φ α satisfying supp (φ α ) ⊆ {t : |t| ∼ 1} such that for all g ∈ S(R d ), The importance for our analysis comes from T α f (x) = D α , g x (S ,S) where g x (ξ) = e ix.ξf (ξ). We are thus reduced to study the operators where φ α ∈ S satisfies suppφ α ⊆ {t : |t| ∼ 1} and δ = 2 −j > 0. Fourier restriction estimates can be applied to T δ , and interpolation with a kernel estimate takes advantage of the decomposition given by Lemma 4.1. The Tomas-Stein restriction theorem (cf. [43,41]) suffices already for the sharp estimates for the restrictionextension operator (α = 1) due to Gutiérrez [25]. Cho et al. [10] made further progress building on Tao's bilinear restriction theorem [42]. The most recent result is due to Kwon-Lee [33] additionally using sharp oscillatory integral estimates by Guth-Hickman-Iliopoulou [23].

4.2.
Bochner-Riesz estimates with negative index for general non-flat surfaces. In this section we extend the analysis to compact pieces of smooth regular hypersurfaces S ⊂ R d with k non-vanishing principal curvatures and k ∈ {1, . . . , d − 1}, d ∈ N, d ≥ 3. Notice that the case d = 2, α > 0 was entirely settled by Bak [2] and Gutiérrez [24,Theorem 1]. Our argument is based on decompositions in Fourier space as in [10,33]. By further localization in Fourier space we may suppose S = {(ξ , ψ(ξ )) : ξ ∈ [−1, 1] d−1 }. Notice that the case of k = 0 corresponds to possibly flat surfaces, where no decay of the Fourier transform can be expected. The case k = d − 1 means that the Gaussian curvature is non-vanishing.
In the special case that all principal curvatures have the same sign, the surface is elliptic and so is ψ.
In the following we show L p -L q boundedness from Theorem 1.3 of the operator for 0 < α < k+2 2 and p, q ∈ [1, ∞] depending on the decay of the Fourier transform of the surface measure and thus by the number of non-vanishing principal curvatures. As above the operator T α for 1 ≤ α < k+2 2 is defined by analytic continuation (cf. [10,33]). We comment on α = k+2 2 after the proof of Lemma 4.4. By Lemma 4.1, we decompose the operator T α as distribution: where φ := φ α ∈ S satisfies suppφ ⊆ {t : |t| ∼ 1}. In view of (26) we have so that it suffices to consider the operators for δ > 0. The contribution away from the surface corresponding to δ 1 or j ≤ 0 in the above display can be estimated by Young's inequality, see below for a precise kernel estimate. This gives summability for j ≤ 0. We focus on the main contribution from j ≥ 0.
We start with using an L 2 -restriction theorem for the surface S. To begin with, we recall the classical result due to Littman [37]; see also [41, Section VIII.5.8], which gives the following decay of the Fourier transform of the surface measure µ: By the T T * -argument (cf. [43,17,31]) this can be recast into an L 2 -L q estimate as already recorded by Greenleaf [22]. The decay of the Fourier transform is crucial for the verification of assumption (ii) in the following special case of the abstract result from [31, Theorem 1.2]: Theorem 4.2 (Keel-Tao). Let (X, dx) be a measure space and H a Hilbert space. Suppose that for each t ∈ R we have an operator U (t) : H → L 2 (X) which satisfies the following assumptions for σ > 0: (i) For all t and f ∈ H we have the energy estimate: (ii) For all t = s and g ∈ L 1 (X) we have the decay estimate Then, for q ≥ 2(1+σ) σ , the estimate The following two lemmas are the key ingredients in the proof of Theorem 1.3. Both rely on (27), which in turn depends on the lower bounds of the k non-vanishing curvatures and ψ C N , χ C N for some large enough N ∈ N. This leads to the claimed stability in Theorem 1.3 of the estimates on ψ and χ.
In the following lemma we apply Theorem 4.2 to T δ and σ = k 2 : Then we have Proof. We perform a linear change of variables to rewrite For the kernel in the inner integral we find by the assumptions on ψ From this and Theorem 4.2, applied to U (t)g(x ) = R d−1 e ix .ξ +tψ(ξ ) χ(ξ )g(ξ ) dξ , we infer

By (29) and Minkowski's inequality, we find
The ultimate estimate follows from the Cauchy-Schwarz inequality, Plancherel's theorem, and inverting the change of variables.
Further estimates for T δ are derived from Integration by parts leads to the following kernel estimate: } and the following estimates hold: Proof. Changing variables ξ d → ξ d + ψ(ξ ) and integrating in ξ d , we have

Sinceφ is supported in {t
So the method of non-stationary phase gives for |x | |x d | Notice that we used δ|x| ≥ δ|x d | 1 holds in this case. On the other hand, (30) We remark that the kernel estimate shows that T α : 2 by the same argument as in [10, Remark 2.3]. With Lemma 4.4 at hand, we may now localize f to cubes of size δ −1 by the following argument, originally due to Fefferman [14]; see also [41, p. 422-423], and [35,10]: Let (Q j ) j∈Z d denote a finitely overlapping covering of R d with cubes of sidelength 2δ −1 centered at jδ −1 and aligned parallel to the coordinate axes. Let C d > 0 be such that |j − k| > C d implies dist(Q j , Q k ) δ −1 |j − k| uniformly with respect to j, k, δ. Then, we obtain If |k − j| > C d , we use the first kernel estimate in (31) and obtain for all N ∈ N Hence, choosing N ∈ N large enough, these terms allow for summation by Young's inequality for series: The penultimate estimate follows from the embedding p → q , p ≤ q, and the last line from the finite overlapping property. For the "diagonal" set, |k − j| ≤ C d , we use (28) as well as Hölder's inequality: Here we have used that the support of f k has measure ∼ δ −d and p ≥ 2. We conclude like above due to the embedding p → q for p ≤ q and the finite overlapping property. Combining the off-diagonal and the diagonal estimates for large enough N , we get for q ≥ 2(2+k) k and 2 ≤ p ≤ q.
By the kernel estimate (31), we find |K δ (x)| δ k+2 2 for all x ∈ R d and thus Next we interpolate (32) and (33) to prove our bounds. To this end we distinguish the cases 1 2 < α < k+2 2 and 0 < α ≤ 1 2 . We obtain weak endpoint estimates using a special case of Bourgain's summation argument (cf. [6,8]). The present version is taken from [10, Lemma 2.5], see also [35,Lemma 2.3] for an elementary proof: For every j ∈ Z let T j be a linear operator, which satisfies Then, for θ, p, q defined by θ = ε2 ε1+ε2 , 1 q = θ q1 + 1−θ q2 and 1 p = θ p1 + 1−θ p2 , the following hold: Interpolating the estimates at the points ( 1 2 , 1 q1 ), 1 q1 ∈ 0, k 2(k+2) from (32) and A := (1, 0) from (33) gives We use this bound for p 1 , p 2 , q 1 , q 2 given by Here, ε > 0 is chosen so small that 1 p1 , 1 p2 ∈ [ 1 2 , 1] holds, which is possible thanks to our assumption 1 2 < α < k+2 2 . So (34) from Lemma 4.5 gives (37) Furthermore, since T 2 −j coincides with its dual, we have under the same conditions on p, q as above: So (36) gives for p 1 = p 2 = 1, and hence, by duality, Since T α coincides with its dual, the estimates (37), (38) imply Finally, we have the trivial strong estimate We refer to Figure 3 for a visualization of the situation. From the estimates (37)- (40) we now derive our claim using the real interpolation identity (cf. [3, Theorem 5.3.1]) as well as the Lorentz space embeddings Lp( In this way, we obtain strong estimates for the operator T α in the interior of the pentagon conv(A, C α,k , B α,k , B α,k , C α,k ) as well on (B α,k , B α,k ): Real interpolation with parameters (θ,q) gives the estimate for (1/p, 1/q) ∈ (B α,k , B α,k ). We have shown strong bounds for p, q such that All these estimates are valid for α > 1 2 . The strong bounds for α = 1 2 can be obtained using Stein's interpolation theorem for analytic families of operators and the estimates for α > 1 2 just proved and for α < 1 2 that we prove below. Proof of Theorem 1.3 (ii): 0 < α < 1 2 . We use the estimates from (32) and the same interpolation procedure as above to find By duality, Again we have the trivial strong estimate (40). Interpolating these estimates as above, we get strong bounds precisely for p, q such that This finishes the proof of Theorem 1.3.

4.3.
Necessary conditions for generalized Bochner-Riesz estimates with negative index. In this subsection we discuss necessary conditions for estimates We shall see that for α ≥ 1/2 the established strong estimates are sharp, but for 0 < α < 1/2 these are in general not. For this purpose, we compare to the estimates for elliptic surfaces in lower dimensions where the bounds are known to be sharp, see [33, p.1419]. Suppose that for d ≥ 3, there is 1 ≤ k ≤ d−1 and (p,q) such that (41) holds true for all regular hypersurfaces with k non-vanishing principal curvatures. Then, let d 1 := k + 1 and let S = {(ξ , ψ(ξ )) ∈ R d1 : ξ ∈ B(0, c)} be an elliptic surface with k = d 1 − 1 positive principal curvatures. This can be trivially embedded into R d considering S = {(ξ , ξ , ψ(ξ )) ∈ R d1+d2 : ξ ∈ B(0, c)}. We consider the operator Apparently, As L α is the kernel of a Bochner-Riesz operator with negative index for an elliptic surface in R k+1 , we know that for 1 2 ≤ α < k+2 2 the corresponding operator Hence, R α : Lp(R d1 ) → Lq(R d1 ) is bounded. By the sharpness of our conditions for elliptic surfaces we infer (1/p, 1/q) ∈ P α (d 1 − 1) = P α (k), which is all we had to show.
On the other hand, we see that the estimates proved in Theorem 1.3 are not sharp for 0 < α < 1/2 as in the elliptic case better estimates are known to hold. Apparently, for 0 < α < 1/2 the geometry of the surface becomes more important. We believe that the optimal estimates will also depend on the difference between positive and negative curvatures as for oscillatory integral operators (cf. [45,7,23]).

Estimates for the regular part
In this section we estimate the contribution of (E, H) with Fourier support close to smooth and regular component of the Fresnel surface by proving Proposition 2.2 and Proposition 2.3. We recall that the first proposition deals with those parts where two principal curvatures are non-zero, whereas the latter proposition deals with frequencies close to the Hamiltonian circles where only one principal curvature is bounded away from zero. As explained in the Introduction, our estimates result from uniform estimates for the Fourier multipliers (P (ξ) + iδ) −1 as δ → 0 with P (ξ) = p(ω, ξ). We stress that ω ∈ R \ {0} is fixed from now on.
We first use our estimates for the Bochner-Riesz operator T α from the previous section to prove a Fourier restriction-extension estimate related to the two parts of the Fresnel surface mentioned above. To carry out the estimates for both parts, we change from implicit to graph representation and apply Theorem 1.3 for (α, k) = (1, 2), respectively (α, k) = (1, 1). The L p -L q -estimates are not affected by this change of representation, see Corollary 5.1. Then we use this result to prove uniform estimates for (P (D) + iδ) −1 by a foliation with level sets of P and the Fourier restriction-extension theorem for the single layer. [10, p. 152] it was stated that a compact convex surface with curvature bounded from below can be written locally as the graph of an elliptic function. Moreover, it was stated that these parametrizations do not affect Bochner-Riesz estimates. To see that this is also true in the non-elliptic case, we explain this in a nutshell.

Parametric representation. Already in
So let M ⊂ R d be a compact part of a smooth regular hypersurface with k non-vanishing curvatures where k ∈ {1, . . . , d − 1}. After finite decompositions and rigid motions, which leave the L p − L q -estimates invariant, we find finitely many local graph representations of M of the form where at least k eigenvalues of the Hessian matrices ∂ 2 ψ(x), x ∈ B(0, c) are bounded away from zero. Taylor's formula gives for ∆ := ξ d − ψ(ξ ) p loc (ξ) = p loc (ξ , ψ(ξ ) + ∆) By the properties of p loc , we find m ∈ C ∞ (B ) with the properties 0 < c 1 ≤ m ≤ c 2 and |∂ γ m| γ 1 for γ ∈ N 2 0 . The Fourier multiplier m α defined by for a suitable cutoff β ∈ C ∞ c (B ), defines a bounded mapping L p (R d ) → L p (R d ), 1 ≤ p ≤ ∞ via Young's convolution inequality. Real interpolation of these estimates also yields the boundedness L p,r (R d ) → L p,r (R d ) for 1 < p < ∞, 1 ≤ r ≤ ∞. Accordingly, choosing a suitable finite partition of unity we find that the operators are well-defined for 0 < α < k+2 2 through analytic continuation and satisfy the same (weak) L p − L q -estimates as the Bochner-Riesz operators that we analyzed in Theorem 1.3. For α = 1 this gives the following: for M t := {ξ ∈ K : P (ξ) = t} and ( 1 p , 1 q ) ∈ {B 1,k , B 1,k }. We have (L p,1 (R d ), L q (R d ))bounds for ( 1 p , 1 q ) ∈ (B 1,k , C 1,k ], (L p (R d ), L q,∞ (R d ))-bounds for ( 1 p , 1 q ) ∈ (B 1,k , C 1,k ] and strong (L p (R d ), L q (R d ))-bounds for ( 1 p , 1 q ) ∈ P 1 (k). As described at the end of Section 3, the principal curvatures of M t vary continuously with respect to t so that the curvature properties of M t for small |t| are inherited from those for t = 0. The estimates leading to the proof of Proposition 2.2 will result from an application of Corollary 5.1 for d = 3, K = supp(β 11 ), k = 2 whereas Proposition 2.3 corresponds to the choice d = 3, K = supp(β 12 ), k = 1. To prove both results simultaneously, we therefore assume that K ⊂ R 3 and k ∈ {1, 2} satisfy the conditions of the corollary.

5.2.
Uniform estimates for the singular multiplier. To prove the desired uniform resolvent estimates for (P (ξ) + iδ) −1 , we consider It is actually enough to show the restricted weak type bound for (1/p 0 , 1/q 0 ) = ( 2(k+1)(k+2) k 2 +6k+4 , k 2(k+1) ) = B and for the remaining tuples (1/p, 1/q) ∈ (B , C ] where C = (1, k 2(k+1) ).  We focus on (42) in the following. To reduce our analysis to the region {ξ ∈ K : |P (ξ)| < t 0 } for t 0 as in Corollary 5.1, we introduce a cut-off function χ ∈ C ∞ c (R d ) such that |P (ξ)| < t 0 for χ(ξ) = 0 and P (ξ) > t 0 /2 for χ(ξ) = 1. We then have Since P is smooth and bounded away from zero on supp(1 − χ), the Fourier multiplier in the latter expression is Schwartz and the claimed estimates (in fact even much stronger ones) hold for this second part. For this reason we may from now on concentrate on the first part. We change to generalized polar coordinates via the coarea formula: In the following we estimate this expression with the aid of Corollary 5.1 by decomposition in Fourier space as in [30, p.346].
The estimate for I(D) is based on the coarea formula, Corollary 5.1, and Young's inequality in Lorentz spaces.
We turn to the estimate of R(D), which requires an additional decomposition: For the existence of φ we refer to the proof of [30, Lemma 2.2], where it is denoted by ψ. We split The coarea formula, Minkowski's inequality, and Corollary 5.1 yield as above Here we used the estimate |φ(s)| s −1 , which holds becauseφ is a Schwartz function. By similar means, we find Here, the estimate from the third to the fourth line uses |φ(2 −j t)| = |φ(2 −j t)|2 −j t 2 −j t. For the most involved estimate of C j , we need the following lemma: and that the level sets {ξ ∈ supp(χ) : P (ξ) = t} have k principal curvatures uniformly bounded from below in modulus for all |t| ≤ t 0 . Then, for 1 ≤ p, q ≤ ∞ with q ≥ 2 and 1 q ≥ k+2 k 1 − 1 p , we find the following estimate to hold for all λ > 0: Proof. By interpolation, it suffices to prove the endpoint estimates for (p, q) = ( 2(k+2) k+4 , 2) and (p, q) = (1, ∞), (p, q) = (1, 2). Since the multiplier is regular for λ ≥ 1, we may henceforth suppose λ ≤ 1. For q = 2 we use Plancherel's theorem, the coarea formula and the L 2(k+2) k+4 -L 2 restriction-extension estimate from Corollary 5.1: .
Using the trivial estimate |f (ξ)| ≤ f L 1 (R d ) instead (from the third to the fourth line), we find the endpoint estimate for (p, q) = (1, 2).
For the endpoint (p, q) = (1, ∞) it suffices to show the kernel estimate The function a is smooth, all its derivatives are bounded functions and its support is bounded with respect to t. So the principle of non-stationary phase yields for |x| In particular, this holds for M = k+2 2 . For |x| λ −1 we can use the dispersive estimate |a(t, x)| (1 + |x|) −k/2 , which holds due to method of stationary phase and the presence of k non-vanishing principal curvatures. We thus get for |x| λ −1 The proof is complete.
The lemma allows to bound the C j -terms as follows: Combining the estimates (43)-(46), we get the claimed estimate This proves Proposition 2.2 (k = 2) and Proposition 2.3 (k = 1).

5.
3. An improved Fourier restriction-extension estimate for the Fresnel surface close to Hamiltonian circles. The purpose of this section is to point out how the special degeneracy along the Hamiltonian circles might allow for improved estimates in Proposition 2.3. In our proof in the previous section we exploited that one principal curvature is bounded away from zero close to these circles. But actually we have more: The other principal curvature does not vanish identically in that region, but only vanishes at the Hamiltonian circle, which is a curve on the Fresnel surface. We refer to Figure 2 for an illustration of the situation. For surfaces with vanishing Gaussian curvature, but no flat points, improved results were established in special cases. For in a sense generic surfaces in R 3 with Gaussian curvature vanishing along a one-dimensional sub-manifold, the decay was shown by Erdős-Salmhofer [13]. (In our proof we used Bochner-Riesz estimates resulting from the weaker bound |μ(ξ)| ξ − 1 2 , which is (27) for k = 1.) We shall show the corresponding L p -L q estimates for these surfaces in future work. However, these results are not applicable in our case. Indeed, one can still show that the gradient of the curvature ∇ Σ K(p) = 0 does not vanish along the Hamiltonian circles (cf. [13,Assumption 2]) and any unit vector has only finitely many preimages under the normal ν : Σ → S 2 (cf. [13,Assumption 3]). It turns out that the transversality assumption [13,Assumption 4] regarding the Hamiltonian circle and the direction of the non-vanishing principal curvature fails: On the curves Γ = {K = 0}, exactly one of the principal curvatures vanishes. We define a (local) unit vectorfield Z ∈ T Σ along Γ in the tangent plane of Σ. On Γ, Z is supposed to point into the direction of the vanishing principal curvature, and Z can be extended to a neighbourhood of Γ as the direction of the principal curvature that is small and vanishes on Γ.
To apply the arguments from [13], it is required that Z is transversal to Γ up to finitely many points, and the angle between Z and Γ increases linearly.
This allows to solve for t = t(s) (cf. (25)): and further, ∂t In the following we shall see that the direction of the vanishing principal curvature is tangential to the Hamiltonian circles. This violates the transversality assumption. For this purpose, consider ξ = ξ(s, t) with t = ε1ε3 ε1+ε3−s . This yields a parametrization of the Hamiltonian circles. For a tangent vector we find dξ ds = ∂ξ ∂s + ∂ξ ∂t ∂t ∂s = e s + ∂t ∂s e t .
A straight-forward computation shows that this is in the kernel of the second fundamental form II, which was computed in (s, t) coordinates in Proposition 3.5: Note that Consequently, From this follows and thus, ∂t ∂s with P L defined in Section 3. Furthermore, and therefore, by plugging the definition of α into P N , Thus, Still, there is hope that one can show better decay for some δ > 0 using stationary phase estimates for functions with degenerate Hessian as in Ikromov-Müller [28, Corollary 1.6] applied by Greenblatt [21]; see also [44,38,20,27,29] and references therein. Since the singular points of our Fresnel surface (to be discussed in the following section) give rise to the worse total decay |μ(ξ)| ξ − 1 2 of the Fourier transform, the analysis is not detailed here.

Estimates for neighbourhoods of the singular points
The purpose of this section is to prove the estimate with β 13 defined in Section 2 as smooth cutoff localizing to a neighbourhood of the singular points. We shall also take the opportunity to derive estimates for perturbed cone multipliers in R d . These naturally arise for surfaces S = {ξ ∈ R d : p(ξ) = 0} at singular points ξ ∈ S with ∇p(ξ) = 0, and ∂ 2 p with signature (1, d − 1). In the first step, to clarify the nature of S, we shall change to parametric representation in Section 6.1. We will see that it suffices to analyze two perturbed half-cones This yields that for a small, but fixed distance from the origin, we have the curvature properties of the cone and can apply Theorem 1.3 with α = 1, k = d − 2 to derive Fourier restriction-extension estimates for the layers. Then, the arguments of Section 5.2 apply again. We derive the estimates for the generalized cone multiplier and (47) by an additional Littlewood-Paley decomposition and a scaling argument in Subsection 6.2.
Coming back to Fresnel's surface, we first prove that S looks like a cone around the singular points. We recall that we assumed without loss of generality µ 1 = µ 2 = µ 3 = ω = 1 so that the results from Section 3 apply for S = S * . Proposition 6.1. Set ζ ∈ S be one of the four singular points given by Proposition 3.2. Then and D 2 p(ω, ζ) has two positive and one negative eigenvalue.
Accordingly, after suitable rotations, translations and multiplication by −1, we may suppose in the following that the analyzed singular point lies in the origin and that the Taylor expansion of the Fourier symbol around the singular point is given byp (ξ) = ξ 2 3 − |ξ | 2 + g(ξ), |∂ α g(ξ)| α |ξ| 3−|α| (α ∈ N 3 0 ). We will discuss the corresponding Fourier multiplier given by where β ∈ C ∞ c (R 3 ). The support of β will later be assumed to be close to zero so that the mapping properties of A δ are determined by the Taylor expansion ofp around zero. The aim is to show estimates for p,q as in Proposition 2.4 as previously independent of δ. This will be proved in Subsection 6.3.
The proposition generalizes Lee's result [35, Theorem 1.1] for α > 1/2: Fourier supports and perturbations of the cone including the singular point are covered and the space dimension is not restricted to d = 3. As we obtain the same conditions on (p, q) as Lee, which he showed to be sharp in the case d = 3, the conditions in Proposition 6.3 are clearly sharp. It seems likely that by bilinear restriction the result can be improved as in [35] for α < 1 2 . To reduce the estimates to Theorem 1.3, we apply a Littlewood-Paley decomposition. Let β l (ξ) = β 0 (2 l ξ) with supp(β 0 ) ⊆ B(0, 2)\B(0, 1/2) and l≥0 β l · β = β.
We define We have the following consequence of Littlewood-Paley theory: Lemma 6.4. Assume that there are 1 < p < 2 < q < ∞ and r 1 ∈ {1, p} and r 2 ∈ {q, ∞} such that C α l f L q,r 2 (R d ) ≤ C f L p,r 1 (R d ) holds for all l ∈ N 0 . Then C α f L q,r 2 (R d ) C f L p,r 1 (R d ) .
Proof. We write by the square function estimate, which also holds in Lorentz spaces, see, e.g. [30,Lemma 3.2], and Minkowski's inequality (note that L q 2 ,∞ is normable because q > 2) C α f L q,r 2 (R d )
By hypothesis and noting that C α l f = C α l |l −l|≤2 β l (D)f , we find l≥0 C α l f 2 L q,r 2 (R d ) Notice that the ultimate estimate is dual to the previous display.
We are ready for the proof of Proposition 6.3.
With m(ξ) = ξ d + |ξ | + r 2,l (ξ ) ξ d for ξ ∈ supp(β 1 ) ∩ supp(β) and |∂ α m(ξ)| 1, by Young's inequality it is enough to considerS α l,1 given by To this operator, we can apply the estimates of Theorem 1.3 for k = d − 2 since in each point of the perturbed cone d − 2 principal curvatures are bounded from below in modulus uniformly with respect to k. Moreover, the rescaled surfaces {ζ d = ∓|ζ | + r i,l (ζ )} can be approximated with the cone in any C N -norm. As a consequence, S α k has the mapping properties described in Theorem 1.3 (i) for 1 2 < α < d 2 with a uniform mapping constant. From C α l f (x) = 2 2αl−3l S α l f l (2 −l x) we conclude Given that the conditions on p, q imply 2α+ d q − d p ≤ 0, we obtain the desired uniform estimates for any fixed l ∈ N 0 . Hence, an application of Lemma 6.4 finishes the proof for p = 1, q = ∞ because of p < 2 < q. If p > 1, q = ∞, we can find q * < ∞ such that the conditions hold for (p, q * ). This is true because 1 p − 1 q = 1 − 1 q * > 2α d for large enough q * . Take χ a cut-off function with χ = 1 on supp(β). Then, by bounded frequencies and Young's inequality, The case p = 1, q < ∞ is dual and thus proved as well. The proof is complete.
These multipliers will be estimated by Fourier restriction-extension estimates for the level sets of p given by Theorem 1.3 for k = d − 2 = 1, α = 1. To carry out the program of Section 5.2, we need to change to generalized polar coordinates ξ = ξ(p, q) in supp (χ). We can suppose that this is possible as |∇p(ξ)| c > 0 for p(ξ) = 0, |ξ| ∼ c, after making the support of β 01 closer to the characteristic set, if necessary. Furthermore, with graph parametrizations (ξ , ψ(ξ )) of {ξ ∈ supp(χ) : p(ξ) = t} uniform in t ∈ (−t 0 , t 0 ), t 0 chosen small enough, Theorem 1.3 yields uniform bounds. Also note that Lemma 5.2 applies with k = 1. This finishes the proof of (47) for = 0. We show the bounds for ≥ 1 by rescaling. A change of variables gives Let p (ζ) = ζ 2 3 − ζ 2 1 − ζ 2 2 + 2 2 g(2 − ζ), δ = 2 2 δ. Recall that |∂ α g(ξ)| |ξ| 3−|α| , which previously allowed to carry out the proof for = 0 for c chosen small enough depending on finitely many C α in (49). Furthermore, we find with implicit constant independent of ≥ 1 choosing c small enough depending only on finitely many C α . Hence, taking (51) and (52) together, gives Hence, Lemma 6.5 applies for p = 1 and q = ∞ because for our choice of p and q we find 2 3 ≤ 1 p − 1 q . For q = ∞ or p = 1, we use that frequencies are compactly supported to reduce to p = 1 and q = ∞ like at the end of the proof of Proposition 6.3. The proof of Proposition 2.4 is complete.