This paper considers deriving new exact solutions of a nonlinear complex generalized Zakharov dynamical system for two different definitions of derivative operators called conformable and M− truncated. The system models the spread of the Langmuir waves in ionized plasma. The extended rational sine−cosine and sinh−cosh methods are used to solve the considered system. The paper also includes a comparison between the solutions of the models containing separately conformable and M− truncated derivatives. The solutions are compared in the 2D and 3D graphics. All computations and representations of the solutions are fulfilled with the help of Mathematica 12. The methods are efficient and easily computable, so they can be applied to get exact solutions of non-linear PDEs (or PDE systems) with the different types of derivatives.
Citation: Melih Cinar, Ismail Onder, Aydin Secer, Mustafa Bayram, Abdullahi Yusuf, Tukur Abdulkadir Sulaiman. A comparison of analytical solutions of nonlinear complex generalized Zakharov dynamical system for various definitions of the differential operator[J]. Electronic Research Archive, 2022, 30(1): 335-361. doi: 10.3934/era.2022018
[1] | Antonio Iannizzotto, Giovanni Porru . Optimization problems in rearrangement classes for fractional p-Laplacian equations. Mathematics in Engineering, 2025, 7(1): 13-34. doi: 10.3934/mine.2025002 |
[2] | Catharine W. K. Lo, José Francisco Rodrigues . On an anisotropic fractional Stefan-type problem with Dirichlet boundary conditions. Mathematics in Engineering, 2023, 5(3): 1-38. doi: 10.3934/mine.2023047 |
[3] | Pier Domenico Lamberti, Michele Zaccaron . Spectral stability of the curlcurl operator via uniform Gaffney inequalities on perturbed electromagnetic cavities. Mathematics in Engineering, 2023, 5(1): 1-31. doi: 10.3934/mine.2023018 |
[4] | David Arcoya, Lucio Boccardo, Luigi Orsina . Hardy potential versus lower order terms in Dirichlet problems: regularizing effects. Mathematics in Engineering, 2023, 5(1): 1-14. doi: 10.3934/mine.2023004 |
[5] | Hugo Tavares, Alessandro Zilio . Regularity of all minimizers of a class of spectral partition problems. Mathematics in Engineering, 2021, 3(1): 1-31. doi: 10.3934/mine.2021002 |
[6] | Gennaro Ciampa, Gianluca Crippa, Stefano Spirito . Propagation of logarithmic regularity and inviscid limit for the 2D Euler equations. Mathematics in Engineering, 2024, 6(4): 494-509. doi: 10.3934/mine.2024020 |
[7] | Lorenzo Brasco . Convex duality for principal frequencies. Mathematics in Engineering, 2022, 4(4): 1-28. doi: 10.3934/mine.2022032 |
[8] | Youchan Kim, Seungjin Ryu, Pilsoo Shin . Approximation of elliptic and parabolic equations with Dirichlet boundary conditions. Mathematics in Engineering, 2023, 5(4): 1-43. doi: 10.3934/mine.2023079 |
[9] | Filippo Gazzola, Gianmarco Sperone . Remarks on radial symmetry and monotonicity for solutions of semilinear higher order elliptic equations. Mathematics in Engineering, 2022, 4(5): 1-24. doi: 10.3934/mine.2022040 |
[10] | Boumediene Abdellaoui, Pablo Ochoa, Ireneo Peral . A note on quasilinear equations with fractional diffusion. Mathematics in Engineering, 2021, 3(2): 1-28. doi: 10.3934/mine.2021018 |
This paper considers deriving new exact solutions of a nonlinear complex generalized Zakharov dynamical system for two different definitions of derivative operators called conformable and M− truncated. The system models the spread of the Langmuir waves in ionized plasma. The extended rational sine−cosine and sinh−cosh methods are used to solve the considered system. The paper also includes a comparison between the solutions of the models containing separately conformable and M− truncated derivatives. The solutions are compared in the 2D and 3D graphics. All computations and representations of the solutions are fulfilled with the help of Mathematica 12. The methods are efficient and easily computable, so they can be applied to get exact solutions of non-linear PDEs (or PDE systems) with the different types of derivatives.
Boundary value problems for fractional-order pseudodifferential operators P, in particular where P is a generalization of the fractional Laplacian (−Δ)a (0<a<1), have currently received much interest in applications, such as in financial theory and probability (but also in mathematical physics and differential geometry), and many methods have been used, most often probabilistic or potential-theoretic methods.
The author has studied such problems by pseudodifferential methods in [8,9,10,11,12,13], under the assumption that the operators satisfy a μ-transmission condition at the boundary of the domain Ω⊂Rn, which allows to show regularity results for solutions of the Dirichlet problem in elliptic cases, to show integration by parts formulas, and much else.
In the present paper we consider translation-invariant pseudodifferential operators (ψdo's) P=Op(p(ξ)) of order 2a>0 with homogeneous symbol p(ξ), which are only taken to satisfy the top-order equation in the μ-transmission condition (relative to the domain Ω=Rn+), we call this the principal μ-transmission condition. It is shown that they retain some of the features: The solution spaces for the homogeneous Dirichlet problem in the elliptic case equal the μ-transmission spaces from [8] (in a setting of low-order Sobolev spaces), having a factor xμn. The integration by parts formula holds (even when P is not elliptic):
∫Rn+Pu∂nˉu′dx+∫Rn+∂nu¯P∗u′dx=Γ(μ+1)Γ(μ′+1)∫Rn−1s0γ0(u/xμn)γ0(ˉu′/xμ′n)dx′, |
when u and u′ are in xμnC∞(¯Rn+) resp. xˉμ′nC∞(¯Rn+) (μ′=2a−μ) and compactly supported.
We also treat nonhomogeneous local Dirichlet problems with Dirichlet trace γ0(u/xμ−1n), and show how the above formula implies a "halfways" Green's formula where one factor has nonzero Dirichlet trace. P can be of any positive order, and μ can be complex.
The results apply in particular to the operator L=Op(A(ξ)+iB(ξ)) with A real, positive and even in ξ, B real and odd in ξ, which satisfies the principal μ-transmission equation for a suitable real μ. Hereby we can compensate for an error made in the recent publication [13] (see also [14]), where it was overlooked that L may not satisfy the full μ-transmission condition when B≠0 (it does so for B=0). The general L are now covered by the present work. They were treated earlier by Dipierro, Ros-Oton, Serrra and Valdinoci [5] under some hypotheses on a and μ; they come up in applications as infinitesimal generators of α-stable n-dimensional Lévy processes, see [5]. (The calculations in [13] are valid when applied to operators satisfying the full μ-transmission condition.)
The study of x-independent ψdo's P on the half-space Rn+ serves as a model case for operators on domains Ω⊂Rn with curved boundary and possible x-dependence, and can be expected to be a useful ingredient in the general treatment, as carried out for the operator L in [5].
Plan of the paper: In Section 2 we give an overview of the aims and results of the paper with only few technicalities. Section 3 introduces the principal transmission condition in detail for homogeneous ψdo symbols. In Section 4, the Wiener-Hopf method is applied to derive basic decomposition and factorization formulas for such symbols. This is used in Section 5 to establish mapping properties for the operators, and regularity properties for solutions of the homogeneous Dirichlet problem in strongly elliptic cases; here μ-transmission spaces (known from [8]) defined in an L2-framework play an important role. Section 6 gives the proof of the above-mentioned integration by parts formula on Rn+. Section 7 treats nonhomogeneous local Dirichlet conditions, and a halfways Green's formula is established.
The study is concerned with the so-called model case, where the pseudodifferential operators have x-independent symbols, hence act as simple multiplication operators in the Fourier transformed space (this frees us from using the deeper composition rules needed for x-dependent symbols), and the considered open subset Ω of Rn is simplest possible, namely Ω=Rn+={x∈Rn∣xn>0}. We assume n≥2 and denote x=(x1,…,xn)=(x′,xn), x′=(x1,…,xn−1). Recall the formulas for the Fourier transform F and the operator P=Op(p(ξ)):
Fu=ˆu(ξ)=∫Rne−ix⋅ξu(x)dx,F−1v=(2π)−n∫Rneix⋅ξv(ξ)dξ,Pu=Op(p(ξ))u=F−1(p(ξ)(Fu)(ξ)). | (2.1) |
We work in L2(Rn) and L2(Rn+) and their derived L2-Sobolev spaces (the reader is urged to consult (5.1) below for notation). On L2(Rn), the Plancherel theorem
‖u‖L2(Rn)=c‖ˆu‖L2(Rn),c=(2π)−n/2, | (2.2) |
makes norm estimates of operators easy. (There is more on Fourier transforms and distribution theory e.g., in [7].) The model case serves both as a simplified special case, and as a proof ingredient for more general cases of domains with curved boundaries, and possibly x-dependent symbols.
The symbols p(ξ) we shall consider are scalar and homogeneous of degree m=2a>0 in ξ, i.e., p(tξ)=tmp(ξ) for t>0, and are C1 for ξ≠0, defining operators P=Op(p).
A typical example is the squareroot Laplacian with drift:
L1=(−Δ)12+bb⋅∇,withsymbolL1(ξ)=|ξ|+ibb⋅ξ, | (2.3) |
where bb=(b1,…,bn) is a real vector. Here m=1, a=12. It satisfies the condition for strong ellipticity, which is:
Rep(ξ)≥c0|ξ|mwithc0>0,allξ∈Rn; | (2.4) |
this is important in regularity discussions. Some results are obtained without the ellipticity hypothesis; as an example we can take the operator L2 with symbol
L2(ξ)=|ξ1+⋯+ξn|+ibb⋅ξ, | (2.5) |
whose real part is zero e.g., when ξ=(1,−1,0,…,0).
The operators are well-defined on the Sobolev spaces over Rn: When p is homogeneous of degree m≥0, there is an inequality
|p(ξ)|≤C|ξ|m≤C⟨ξ⟩m,⟨ξ⟩=(1+|ξ|2)12 |
(we say that p is of order m); then
‖Pu‖L2(Rn)=c‖p(ξ)ˆu(ξ)‖L2(Rn)≤cC‖⟨ξ⟩mˆu‖L2(Rn)=C′‖u‖Hm(Rn), | (2.6) |
so P maps Hm(Rn) continuously into L2(Rn). Similarly, it maps Hs+m(Rn) continuously into Hs(Rn) for all s∈R.
But for these pseudodifferential operators it is not obvious how to define them relative to the subset Rn+, since they are not defined pointwise like differential operators, but by integrals (they are nonlocal). The convention is here to let them act on suitable linear subsets of L2(Rn+), where we identify L2(Rn+) with the set of u∈L2(Rn) that are zero on Rn−, i.e., have their support suppu⊂¯Rn+. (The support suppu of a function or distribution u is the complement of the largest open set where u=0. The operator that extends functions on Rn+ by zero on Rn− is denoted e+.) Then we apply P and restrict to Rn+ afterwards; this is the operator r+P. (r+ stands for restriction from Rn to Rn+.)
Aiming for the integration by parts formula mentioned in the start, we have to clarify for which functions u,u′ the integrals make sense. It can be expected from earlier studies ([5,10,20]) that the integral will be meaningful for solutions of the so-called homogeneous Dirichlet problem on Rn+, namely the problem
r+Pu=fonRn+,u=0onRn− | (2.7) |
(where the latter condition can also be written suppu⊂¯Rn+). This raises the question of where r+P lands; which f can be prescribed? Or, if f is given in certain space, where should u lie in order to hit the space where f lies?
Altogether, we address the following three questions on P:
(1) Forward mapping properties. From which spaces does r+P map into an Hs-space for f?
(2) Regularity properties. If u solves (2.7) with f in an Hs-space for a high s, will u then belong to a space with a similar high regularity?
(3) Integration by parts formula for functions in spaces where r+P is well-defined.
It turns out that the answers to all three points depend profoundly on the introduction of so-called μ-transmission spaces. To explain their importance, we turn for a moment to the fractional Laplacian which has a well-established treatment:
For the case of (−Δ)a, 0<a<1, it was shown in [8] that the following space is relevant:
Ea(¯Rn+)=e+xanC∞(¯Rn+). | (2.8) |
It has the property that (−Δ)a maps it to C∞(¯Rn+); more precisely,
r+(−Δ)amapsEa(¯Rn+)∩E′(Rn)intoC∞(¯Rn+). | (2.9) |
Here E′(Rn) is the space of distributions with compact support, so the intersection with this space means that we consider functions in Ea that are zero outside a compact set.
For Sobolev spaces, it was found in [8] that the good space for u is the so-called a-transmission space Ha(t)(¯Rn+); here
r+(−Δ)amapsHa(t)(¯Rn+)∩E′(Rn)into¯Ht−2a(Rn+), | (2.10) |
for all t≥a (say). Ea(¯Rn+)∩E′(Rn) is a dense subset of Ha(t)(¯Rn+). The definition of the space Ha(t)(¯Rn+) is recalled below in (2.15) and in more detail in Section 5.3; let us for the moment just mention that it is the sum of the space ˙Ht(¯Rn+) and a certain subspace of xan¯Ht−a(Rn+). This also holds when a is replaced by a more general μ.
For (−Δ)a, the a-transmission spaces provide the right answers to question (1), and they are likewise right for question (2) (both facts established in [8]), and there are integration by parts formulas for (−Δ)a applied to elements of these spaces, [10,11].
The key to the proofs is the so-called a-transmission condition that (−Δ)a satisfies; it is an infinite list of equations for p(ξ) and its derivatives, linking the values on the interior normal to Rn+ with the values on the exterior normal. We formulate it below with a replaced by a general μ.
Definition 2.1. Let μ∈C, and let p(ξ) be homogeneous of degree m. Denote the interior resp. exterior normal to the boundary of Rn+ by (0,±1)={(ξ′,ξn)∣ξ′=0,ξn=±1}.
1∘ p (and P=Op(p)) is said to satisfy the principal μ-transmission condition at Rn+ if
p(0,−1)=eiπ(m−2μ)p(0,1). | (2.11) |
2∘ p (and P=Op(p)) is said to satisfy the μ-transmission condition at Rn+ if
∂αξp(0,−1)=eiπ(m−2μ−|α|)∂αξp(0,1),forallα∈Nn0. | (2.12) |
Note that μ is determined from p in (2.11) up to addition of an integer, when p(0,1)≠0.
The operators considered on smooth domains Ω in [8] were assumed to satisfy (2.12) (for the top-order term p0 in the symbol) at all boundary points x0∈∂Ω, with (0,1) replaced by the interior normal ν at x0, and (0,−1) replaced by −ν. The lower-order terms pj in the symbol, homogeneous of degree m−j, should then satisfy analogous rules with m−j instead of m.
The principal μ-transmission condition (2.11) is of course much less demanding than the full μ-transmission condition (2.12). What we show in the present paper is that when (2.11) holds, the μ-transmission spaces are still relevant, and provide the appropriate answers to both questions (1) and (2), however just for t (the regularity parameter) in a limited range. This range is large enough that integration by parts formulas can be established, answering (3).
By simple geometric considerations one finds:
Proposition 2.2. 1∘ When p(ξ) is homogeneous of degree m, there is a μ∈C, uniquely determined modulo Z if p(0,1)≠0, such that (2.11) holds.
2∘ If moreover, p is strongly elliptic (2.4) and m=2a>0, μ can be chosen uniquely to satisfy μ=a+δ with |Reδ|<12.
This is shown in Section 3. From here on we work under two slightly different assumptions. The symbol p(ξ) is in both cases taken homogeneous of degree m=2a>0 and C1 for ξ≠0. We pose Assumption 3.1 requiring that p is strongly elliptic and μ is chosen as in Proposition 2.2 2∘. We pose Assumption 3.2 just requiring that μ is defined according to Proposition 2.2 1∘. In all cases we write μ=a+δ, and define μ′=a−δ=2a−μ.
Example 2.3. Consider L1=|ξ|+ibb⋅ξ defined in (2.3). The order is 1=m=2a, so a=12. Here L1(0,1)=1+ibn and L1(0,−1)=1−ibn. The angle θ in C=R2 between the positive real axis and 1+ibn is θ=Arctanbn. Set δ=θ/π, then
L1(0,1)=eiπδ|L1(0,1)|=eiπδ(1+|bb|2)12, similarly L1(0,−1)=e−iπδ|L1(0,−1)|=e−iπδ(1+|bb|2)12. |
Moreover,
L1(0,−1)/L1(0,1)=e−2iπδ=eiπ(2a−2(a+δ)), when a=12, |
so (2.11) holds with m=2a=1, μ=12+δ, where δ=1πArctanbn, and Assumption 3.1 is satisfied. Note that δ∈]−12,12[.
For L2 in (2.5), the values at (0,1) and (0,−1) are the same as the values for L1, so (2.11) holds with the same values, and Assumption 3.2 is satisfied. But not Assumption 3.1 since L2 is not strongly elliptic.
When bn≠0, hence δ≠0, neither of these symbols satisfy the full μ-transmission condition Definition 2.1 2∘, since second derivatives remove the (ibb⋅ξ)-term so that the resulting symbol is even (with μ=a+δ replaced by μ=a). $
Our answer to (1) is now the following (achieved in Section 5.4):
Theorem 2.4. Let P satisfy Assumption 3.2. For Reμ−12<t<Reμ+32, r+P defines a continuous linear mapping
r+P:Hμ(t)(¯Rn+)→¯Ht−2a(Rn+). | (2.13) |
It is important to note that r+P then also makes good sense on subsets of Hμ(t)(¯Rn+). In particular, since Eμ(¯Rn+)∩E′(Rn) is a subset of Hμ(t)(¯Rn+) for all t, the operator r+P is well-defined on Eμ(¯Rn+)∩E′(Rn), mapping it into ⋂t<Reμ+32¯Ht−2a(Rn+)⊂¯HReδ+32−a−ε(Rn+), any ε>0, by (2.13). When Reδ>−12 (always true under Assumption 3.1), this is assured to be contained in ¯H1−a(Rn+).
Our answer to (2) is (cf. Section 5.4):
Theorem 2.5. Let P satisfy Assumption 3.1. Then P=ˆP+P′, where P′ is of order 2a−1, and r+ˆP is a bijection from Hμ(t)(¯Rn+) to ¯Ht−2a(Rn+) for Reμ−12<t<Reμ+32. In other words, there is unique solvability of (2.7) with P replaced by ˆP, in the mentioned spaces.
For r+P itself, there holds the regularity property: Let Reμ−12<t<Reμ+32, let f∈¯Ht−2a(Rn+), and let u∈˙Hσ(¯Rn+) (for some σ>Reμ−12) solve the homogeneous Dirichlet problem (2.7). Then u∈Hμ(t)(¯Rn+).
The last statement shows a lifting of the regularity of u in the elliptic case, namely if it solves (2.7) lying in a low-order space ˙Hσ(¯Rn+), then it is in the best possible μ-transmission space according to Theorem 2.4, mapping into the given range space ¯Ht−2a(¯Rn+). In other words, the domain of the homogeneous Dirichlet problem with range in ¯Ht−2a(¯Rn+) equals Hμ(t)(¯Rn+).
The strategy for both theorems is, briefly expressed, as follows: The first step is to replace P=Op(p(ξ)) by ˆP=Op(ˆp(ξ)), where ˆp(ξ) is better controlled at ξ′=0 and p′(ξ)=p(ξ)−ˆp(ξ) is O(|ξ|2a−1) for |ξ|→∞. The second step is to reduce ˆP to order 0 by composition with "plus/minus order-reducing operators" Ξt±=Op((⟨ξ′⟩±iξn)t) ((3.11), (5.2)) geared to the value μ (recall μ′=2a−μ):
ˆQ=Ξ−μ′−ˆPΞ−μ+. | (2.14) |
Then the homogeneous symbol q associated with ˆQ satisfies the principal 0-transmission condition. The third step is to decompose ˆQ into a sum (when Assumption 3.2 holds) or a product (when Assumption 3.1 holds) of operators whose action relative to the usual Sobolev spaces ˙Hs(¯Rn+) and ¯Hs(Rn+) can be well understood, so that we can show forward mapping properties and (in the strongly elliptic case) bijectiveness properties for ˆQ. The fourth step is to carry this over to forward mapping properties and (in the strongly elliptic case) bijectiveness properties for ˆP. The fifth and last step is to take P′=P−ˆP back into the picture and deduce the forward mapping resp. regularity properties for the original operator P.
It is the right-hand factor Ξ−μ+ in (2.14) that is the reason why the μ-transmission spaces, defined by
Hμ(t)(¯Rn+)=Ξ−μ+e+¯Ht−Reμ(Rn+), | (2.15) |
enter. Here e+¯Ht−Reμ(Rn+) has a jump at xn=0 when t>Reμ+12, and then the coefficient xμn appears.
The analysis of ˆQ is based on a Wiener-Hopf technique (cf. Section 4) explained in Eskin's book [6], instead of the involvement of the extensive Boutet de Monvel calculus used in [8].
An interesting feature of the results is that the μ-transmission spaces have a universal role, depending only on μ and not on the exact form of P.
Finally, we answer (3) by showing an integration by parts formula, based just on Assumption 3.2.
Theorem 2.6. Let P satisfy Assumption 3.2, and assume moreover that Reμ>−1, Reμ′>−1. For u∈Eμ(¯Rn+)∩E′(Rn), u′∈Eˉμ′(¯Rn+)∩E′(Rn), there holds
∫Rn+Pu∂nˉu′dx+∫Rn+∂nu¯P∗u′dx=Γ(μ+1)Γ(μ′+1)∫Rn−1s0γ0(u/xμn)γ0(ˉu′/xμ′n)dx′, | (2.16) |
where s0=e−iπδp(0,1). The formula extends to u∈Hμ(t)(¯Rn+), u′∈Hˉμ′(t′)(¯Rn+), for t>Reμ+12, t′>Reμ′+12.
The integrals over Rn+ in (2.16) are interpreted as dualities when needed. The basic step in the proof is the treatment of one order-reducing operator in Proposition 6.1, by an argument shown in detail in [10,Th. 3.1,Rem. 3.2], and recalled in [13,Th. 4.1].
In the proof of (2.16) in Section 6, the formula is first shown for the nicer operator ˆP, and thereafter extended to P. (The formula (2.16) for (−Δ)a in Ros-Oton and Serra [20,Th. 1.9] should have a minus sign on the boundary contribution; this has been corrected by Ros-Oton in the survey [19,p. 350].)
The theory will be carried further, to include "large" solutions of a nonhomogeneous local Dirichlet problem, and to show regularity results and a "halfways Green's formula", see Section 7, but we shall leave those aspects out of this preview.
The example L1 in (2.3) is a special case of the operator L=Op(L(ξ)), where L(ξ)=A(ξ)+iB(ξ) with A(ξ) real, even in ξ and positive, and B(ξ) real and odd in ξ. There are more details below in (3.5)ff. (this stands for (3.5) and the near following text) and Examples 5.9, 6.5, 7.4. L was first studied in [5] (under certain restrictions on a and μ), and our results apply to it. Theorem 2.6 gives an alternative proof for the same integration by parts formula, established in [5,Prop. 1.4] by extensive real function-theoretic methods.
The result on the integral over Rn+ is combined in [5] with localization techniques to get an interesting result for curved domains, and it is our hope that the present results for more general strongly elliptic operators can be used in a similar way.
Let p(ξ) be a complex function on Rn that is homogeneous of degree m in ξ, and let ν∈Rn be a unit vector. For a complex number μ, we shall say that p satisfies the principal μ-transmission condition in the direction ν, when
p(−ν)=eiπ(m−2μ)p(ν). | (3.1) |
When p(ν)≠0, we can rewrite (3.1) as
eiπ(m−2μ)=p(−ν)p(ν),i.e.,μ=m2−12πilogp(−ν)p(ν), |
where log is a complex logaritm. This determines the possible μ up to addition of an integer.
The (full) μ-transmission property defined in [8] demands much more, namely that
∂αξp(−ν)=eiπ(m−2μ−|α|)∂αξp(ν),allα∈Nn0. | (3.2) |
Besides assuming infinite differentiability, this is a stronger condition than (3.1) in particular because of the requirements it puts on derivatives of p transversal to ν.
To analyse this we observe that when a (sufficiently smooth) function f(t) on R∖{0} is homogeneous of degree m∈R, then it has the form, for some c1,c2∈C,
f(t)={c1tmfort>0,c2(−t)mfort<0, |
and its derivative outside t=0 is a function homogeneous of degree m−1 satisfying
∂tf(t)={c1mtm−1fort>0,−c2m(−t)m−1fort<0. |
In particular, if c1≠0, m≠0,
f(−1)/f(1)=c2/c1,∂tf(−1)/∂tf(1)=−c2/c1. |
In the case m=0, f is constant for t>0 and t<0, and the derivative is zero there.
Thus, when p(ξ) is a (sufficiently smooth) function on Rn∖{0} that is homogeneous of degree m≠0, and we consider it on a two-sided ray {tν∣t∈R} where ν is a unit vector and p(ν)≠0, then
p(−ν)=c0p(ν)⟹∂tp(tν)|t=−1=−c0∂tp(tν)|t=1. | (3.3) |
So for example, when ν is the inward normal (0,1)={(ξ′,ξn)∣ξ′=0,ξn=1} to Rn+,
p(0,−1)=c0p(0,1)⟹∂ξnp(0,−1)=−c0∂ξnp(0,1). |
For p(ξ) satisfying (3.1), this means that when p(ν)≠0, it will also satisfy
∂tp(tν)|t=−1=eiπ(m−2μ−1)∂tp(tν)|t=1, |
in view of (3.3). This argument can be repeated, showing that
∂ktp(tν)|t=−1=eiπ(m−2μ−k)∂ktp(tν)|t=1, | (3.4) |
as long as the derivatives at t=1 do not vanish. That can happen when m is a nonnegative integer (namely from the (m+1)'st step on); then (3.4) is trivially satisfied. On the other hand, we cannot infer that derivatives of ∂αξp for arbitrary α have the property (3.2); this will be illustrated in examples below.
In general, μ takes different values for different ν. When Ω is a sufficiently smooth subset of Rn with interior normal ν(x) at boundary points x∈∂Ω, we say that p satisfies the principal μ-transmission condition at Ω if μ(x) is a function on ∂Ω such that (3.1) holds with this μ(x) at boundary points x∈∂Ω. For Ω=Rn+, the normal ν equals (0,1) at all boundary points and μ is a constant; this is the situation considered in the present paper.
In [13] we have studied a special class of symbols first considered by Dipierro, Ros-Oton, Serra and Valdinoci in [5]:
L(ξ)=A(ξ)+iB(ξ), | (3.5) |
the functions being C∞ for ξ≠0 and homogeneous in ξ of degree 2a>0 (a<1), and where A(ξ) is real and even in ξ (i.e., A(−ξ)=A(ξ)), B(ξ) is real and odd in ξ (i.e., B(−ξ)=−B(ξ)), and L is strongly elliptic (i.e., A(ξ)>0 for ξ≠0). As shown in [13,Sect. 2], L satisfies (3.1) on each unit vector ν, for m=2a and
μ(ν)=a+δ(ν),withδ(ν)=1πArctanb,b=B(ν)/A(ν); | (3.6) |
this follows straightforwardly (as in Example 2.3) from the observation that L(−ν)/L(ν)=(1−ib)/(1+ib), b=B(ν)/A(ν). It then also satisfies (3.4) with this μ.
But the full μ-transmission condition need not hold. For example, the symbol L1(ξ)=|ξ|+ibb⋅ξ in (2.3) (with bb∈Rn) satisfies the principal μ-transmission condition for ν=(0,1) with μ=12+δ, δ≠0 if bn≠0, whereas
∂2ξ1L1=(ξ22+⋯+ξ2n)/|ξ|3 |
and its derivatives satisfy the conditions in (3.2) for ν=(0,1) with μ replaced by 12.
The statement in [13,Th. 3.1] that solutions of the homogeneous Dirichlet problem have a structure with the factor xμn, was quoted from [8] based on the full μ-transmision condition, and therefore applies to L=Op(L) when B=0 (a case belonging to [8]), but not in general when B≠0. Likewise, the integration by parts formulas for L derived in [13] using details from the Boutet de Monvel calculus are justified when B=0 or when other operators P satisfying the full μ-transmission condition are inserted, but not in general when B≠0. Fortunately, there are cruder methods that do lead to such results, on the basis of the principal μ-transmission condition alone, and that is what we show in this paper.
The treatment of L will be incorporated in a treatment of general strongly elliptic homogeneous symbols in the following. This requires that we allow complex values of μ.
Let P=Op(p(ξ)) be defined by (2.1) from a symbol p(ξ) that is C1 for ξ≠0, homogeneous of order m=2a>0, and now also strongly elliptic (2.4). To fix the ideas, we shall consider the operator relative to the set Rn+, with interior normal ν=(0,1). Denote p(ξ)|ξ|−2a=p1(ξ); it is homogeneous of degree 0. Both p and p1 take values in a closed subsector of {z∈C∣Reξn>0}∪{0}. For any ξ′∈Rn−1, one has for +1 and −1 respectively,
limξn→±∞p1(ξ′,ξn)=limξn→±∞p1(ξ′/|ξn|,±1)=p1(0,±1)=p(0,±1). |
With the logarithm logz defined to be positive for real z>1, with a cut along the negative real axis, denote logp(0,±1)=α±; here Reα±=log|p(0,±1)| and Imα± is the argument of p(0,±1). With this notation,
p(0,−1)/p(0,1)=eα−/eα+=eα–α+, |
so (3.1) for m=2a holds with ν=(0,1) when α–α+=iπ(2a−2μ), i.e.,
μ=a+δwithδ=(α+−α−)/2πi; | (3.7) |
this μ is the factorization index. These calculations were given in [8,Sect. 3] (with m=2a), and are in principle consistent with the determination of the factorization index by Eskin in [6,Ex. 6.1] (which has different plus/minus conventions because of a different definition of the Fourier transform).
Since p(ξ) takes values in {Rez>0} for ξ≠0, both p(0,1) and p(0,−1) lie there and the difference between their arguments is less than π, so |Im(α+−α−)/2π|<12; in other words
|Reδ|<12. | (3.8) |
Note that δ is real in the case (3.5).
We collect the information on P in the following description:
Assumption 3.1. The operator P=Op(p(ξ)) is defined from a symbol p(ξ) that is C1 for ξ≠0, homogeneous of order m=2a>0, and strongly elliptic (2.4). It satisfies the principal μ-transmission condition in the direction (0,1):
p(0,−1)=eiπ(m−2μ)p(0,1), |
with μ equal to the factorization index μ=a+δ derived around (3.7), and |Reδ|<12. Denote μ′=2a−μ=a−δ.
In Eskin's book [6], the case of constant-coefficient pseudodifferential operators considered on Rn+ is studied in §§–17, and the calculations rely on the principal transmission condition up to and including §. From § on, additional conditions on transversal derivatives are required (the symbol class D(0)α+iβ seems to correspond to our full 0-transmission condition, giving operators preserving smoothness up to the boundary). In the following, we draw on some of the points made in §§–7 there.
For an operator A defined from a homogeneous symbol a(ξ), the behavior at zero can be problematic to deal with. In [6,§7], there is introduced a technique that leads to a nicer operator, in the context of operators relative to Rn+: One eliminates the singularity at ξ′=0 by replacing the homogeneous symbol a(ξ′,ξn) by
ˆa(ξ′,ξn)=a(⟨ξ′⟩ξ′/|ξ′|,ξn), | (3.9) |
the corresponding operator denoted ˆA. (In comparison with [6] we have replaced the factor 1+|ξ′| used there by ⟨ξ′⟩=(1+|ξ′|2)12.) It is shown there that when a(ξ) is homogeneous of degree α+iβ, then
a′(ξ)=a(ξ)−ˆa(ξ)isO(|ξ|α−1)for|ξ|≥2, | (3.10) |
hence is of lower order in a certain sense. Many results with Sobolev estimates are then shown primarily for the "hatted" version ˆA=Op(ˆa), and supplied afterwards with information on A′=Op(a′). Indeed, we shall see that the results we are after for our operators P=Op(p), can be obtained in a manageable way for ˆP=Op(ˆp), and then extended to P by a supplementing analysis of P′. The important thing is that special properties with respect to ξn, such as holomorphic extendability into C+ or C−, are not disturbed when a is replaced by ˆa.
Some of the results that we shall show do not require ellipticity of P. We therefore introduce also a weaker assumption:
Assumption 3.2. The operator P=Op(p(ξ)) is defined from a symbol p(ξ) that is C1 for ξ≠0, homogeneous of order m=2a>0, and satisfies the principal μ-transmission condition in the direction (0,1) with μ=a+δ for some δ∈C. Denote a−δ=μ′.
For the symbols p considered in the rest of the paper, we assume at least that Assumption 3.2 holds. As noted earlier, when P satisfies (3.1) for some μ, it also does so with μ replaced by μ+k, k∈Z. The precision in Assumption 3.1, that μ should equal the factorization index, is needed for elliptic solvability statements.
Consider the symbols of "order-reducing" operators (more on them in Section 5):
χt0,±(ξ)=(|ξ′|±iξn)t;consequentlyˆχt0,±(ξ)=(|⟨ξ′⟩ξ′/|ξ′||±iξn)t=(⟨ξ′⟩±iξn)t=χt±(ξ); | (3.11) |
the last entry is the usual notation. Together with our symbol p(ξ) of order 2a, we shall consider its reduction to a symbol q of order 0 defined by:
q(ξ)=χ−μ′0,−p(ξ)χ−μ0,+,herebyp(ξ)=χμ′0,−q(ξ)χμ0,+. | (3.12) |
The "hatted" version is:
ˆq(ξ)=χ−μ′−ˆp(ξ)χ−μ+,herebyˆp(ξ)=χμ′−ˆq(ξ)χμ+. | (3.13) |
Here q is continuous and homogeneous of degree 0 for ξ≠0; it is C1 in ξn there, and C1 in ξ′ for ξ′≠0 with bounded first derivatives on |ξ|=1. Since i=eiπ/2,
q(0,1)=(−i)μ−2ap(0,1)i−μ=i2a−2μp(0,1)=eiπ(a−μ)p(0,1),q(0,−1)=(+i)μ−2ap(0,−1)(−i)−μ=i2μ−2aeiπ(2a−2μ)p(0,1)=eiπ(a−μ)p(0,1)=q(0,1), |
so q satisfies the principal 0-transmission condition in the direction ν=(0,1):
q(0,−1)=q(0,1). | (3.14) |
In view of (3.1)–(3.4), we have moreover when p(0,1)≠0 that
∂ξnq(0,−1)=−∂ξnq(0,1). | (3.15) |
Note that since μ−a=δ, q(0,1)=e−iπδp(0,1). We shall denote
s0=q(0,1)=e−iπδp(0,1). | (3.16) |
In the case p=L in (3.5)–(3.6), L(0,1)=eiπδ|L(0,1)| with δ real, so
s0=e−iπδL(0,1)=|L(0,1)|=(A(0,1)2+B(0,1)2)12then. | (3.17) |
Since p(ξ) is only assumed to satisfy the principal μ-transmission condition, q(ξ) will in general only satisfy the principal 0-transmission condition, not the full one, so the techniques of the Boutet de Monvel calculus brought forward in [8] are not available. Instead we go back to a more elementary application of the original Wiener-Hopf method [22].
When b(ξn) is a function on R, denote
b+(ξn+iτ)=i2π∫Rb(ηn)ηn−ξn−iτdηnforτ<0,b−(ξn+iτ)=−i2π∫Rb(ηn)ηn−ξn−iτdηnforτ>0, | (4.1) |
when the integrals have a sense. When b is suitably nice, b+ is holomorphic in ξn+iτ for τ<0 and extends to a continuous function on ¯C− (also denoted b+), b− has these properties relative to ¯C+, and b(ξn)=b+(ξn)+b−(ξn) on R. With the notation of spaces H, H± introduced by Boutet de Monvel in [4], denoted H, H± in our subsequent works, the decomposition holds for b∈H with b±∈H± on R. Since we are presently dealing with functions with cruder properties, we shall instead apply a useful lemma shown in [6,Lemma 6.1]:
Lemma 4.1. Suppose that b(ξ′,ξn) is homogeneous of degree 0 in ξ, is C1 for ξ′≠0, and satisfies
|b(ξ′,ξn)|≤C|ξ′||ξ|−1,|∂jb(ξ′,ξn)|≤C|ξ|−1forj≤n−1. | (4.2) |
Then the function defined for τ<0 by
b+(ξ′,ξn+iτ)=i2π∫Rb(ξ′,ηn)ηn−ξn−iτdηn | (4.3) |
is holomorphic with respect to ξn+iτ in C−, is homogeneous of degree 0, extends by continuity with respect to (ξ′,ξn+iτ)∈¯C− for |ξ|+|τ|>0, τ≤0, and satisfies the estimate
|b+(ξ′,ξn+iτ)|≤Cε|ξ′|1−ε(|ξ|+|τ|)ε−1,anyε>0. | (4.4) |
There is an analogous statement for b− with C− replaced by C+.
The symbol q derived from p by (3.12) satisfies
q(ξ)=s0+f(ξ), |
where f is likewise homogeneous of degree 0, and has f(0,1)=f(0,−1)=0. We make two applications of Lemma 4.1. One is, under Assumption 3.2, to apply it directly to f to get a sum decomposition f=f++f− where the terms extend holomorphically to C− resp. C+ with respect to ξn; this will be convenient in establishing the forward mapping properties and integration by parts formula for the present operators. The other is, under Assumption 3.1, to apply the lemma to the function b(ξ)=logq(ξ) to get a sum decomposition of b and hence a factorization of q; this is used to show that P has appropriate solvability properties (the solutions exhibiting a singularity xμn at the boundary).
We show that f has the properties required for Lemma 4.1 as follows: To see that (4.2) is verified by f, note that the second inequality follows since ∂jf is bounded on the unit sphere {|ξ|=1} and homogeneous of degree −1. For the first inequality we have, when ξn>|ξ′| (hence |ξ′/ξn|<1),
|f(ξ′,ξn)|=|q(ξ′ξn,1)−q(0,1)|≤∑j<n|ξjξn|sup|η′|≤1|∂jq(η′,1)|≤C|ξ′||ξn|≤C′|ξ′||ξ|, | (4.5) |
using the mean value theorem and the fact that |ξn|∼|ξ| when |ξn|≥|ξ′|. A similar estimate is found for ξn<−|ξ′|. For |ξn|≤|ξ′|, we use that q is bounded, so that |q(ξ)−s0||ξ|/|ξ′|≤c|q(ξ)−s0||ξ′|/|ξ′|≤c′. We have obtained:
Proposition 4.2. When p satisfies Assumption 3.2 and q is derived from p by (3.12), then there is a sum decomposition of f=q−s0:
q(ξ)−s0=f+(ξ)+f−(ξ), |
where f+(ξ′,ξn) is holomorphic with respect to ξn+iτ in C−, and continuous with respect to (ξ′,ξn+iτ)∈¯C− for |ξ|+|τ|>0, τ≤0, and satisfies estimates
|f+(ξ′,ξn+iτ)|≤Cε|ξ′|1−ε(|ξ|+|τ|)ε−1,anyε>0, | (4.6) |
and f− has the analogous properties with C− replaced by C+.
For the corresponding hatted symbol, we then have ˆq=s0+ˆf++ˆf−, with ˆf± defined from f±. They have similar holomorphy properties, and satisfy estimates as in (4.6) with |ξ′| replaced by ⟨ξ′⟩.
In order to obtain a factorization for symbols satisfying Assumption 3.1, we shall study logq. By the strong ellipticity, q(ξ)≠0 for ξ≠0. Moreover, p(ξ)|ξ|−2a=χ−a0,−p(ξ)χ−a0,+ takes values in a subsector of {z∈C∣Rez>0} and the multiplication by χδ0,− and χ−δ0,+ gives the function q taking values in the sector {z∈C∣|argz|≤π(12+|Reδ|)} disjoint from the negative real axis. So the logarithm is well-defined with inverse exp.
Assume first that s0=1; this can simply be obtained by dividing out q(0,1). The function b(ξ)=logq(ξ) is homogeneous of degree 0 and has b(0,1)=b(0,−1)=0 and the appropriate continuity properties, and bounds on first derivatives, so the same proof as for f applies to b to give the decomposition b=b++b−. Then we define q±=exp(b±), they are homogeneous of degree 0. For example,
q+=1+g+,whereg+=∑k≥1(b+)k. |
Here |b+(ξ)|≤Cε|ξ′|1−ε|ξ|−1+ε, and there is a constant C′ε such that Cε|ξ′|1−ε|ξ|−1+ε≤12 for |ξn|≥C′ε|ξ′|. On this set the series for g+ converges with |g+|≤|b+|, hence g+ satisfies an estimate of the form (4.4) there. It likewise does so on the set |ξn|≤C′ε|ξ′| since |ξ′|∼|ξ| there. There are similar results for q−=exp(b−)=1+g− with C− replaced by C+. This shows:
Proposition 4.3. When p satisfies Assumption 3.1 and q is derived from p by (3.12) and satisfies s0=1, then there is a factorization of q:
q(ξ)=q−(ξ)q+(ξ), |
where q+(ξ′,ξn) is holomorphic with respect to ξn+iτ in C−, and continuous with respect to (ξ′,ξn+iτ)∈¯C− for |ξ|+|τ|>0, τ≤0. Moreover, g+=q+−1 satisfies estimates
|g+(ξ′,ξn+iτ)|≤Cε|ξ′|1−ε(|ξ|+|τ|)ε−1,allε>0, | (4.7) |
and q−, g−=q–1 have the analogous properties with C− replaced by C+. The symbols are homogeneous of degree 0, and q+ and q− are elliptic.
For general s0, we apply the factorization to q0=s−10q, so that q0=q−0q+0; then q=q−q+ with q−=s0q−0=s0(1+g−) and q+=q+0=1+g+.
The ellipticity follows from the construction as exp(b±), or one can observe that the product q+q−=q is elliptic (i.e., nonzero for ξ≠0).
The notation with upper index ± is chosen here to avoid confusion with the lower + used later to indicate truncation, P+=r+Pe+.
Turning to the corresponding hatted symbols, we have obtained ˆq=ˆq−ˆq+, with ˆq±, ˆg± defined from q±, g±, respectively. They have similar holomorphy properties, the ˆq± are elliptic, and the ˆg± satisfy estimates as in (4.7) with |ξ′| replaced by ⟨ξ′⟩:
|ˆg+(ξ′,ξn+iτ)|≤Cε⟨ξ′⟩1−ε(⟨ξ⟩+|τ|)ε−1,allε>0. | (4.8) |
First recall some terminology: E′(Rn) is the space of distributions on Rn with compact support, S(Rn) is the Schwartz space of C∞-functions f on Rn such that xβDαf is bounded for all α,β, and S′(Rn) is its dual space of temperate distributions. ⟨ξ⟩ stands for (1+|ξ|2)12. We denote by r+ the operator restricting distributions on Rn to distributions on Rn+, and by e+ the operator extending functions on Rn+ by zero on Rn∖Rn+. Then r+S(Rn) is denoted S(¯Rn+). The following notation for L2-Sobolev spaces will be used, for s∈R:
Hs(Rn)={u∈S′(Rn)∣⟨ξ⟩sFu∈L2(Rn)},¯Hs(Rn+)=r+Hs(Rn),therestrictedspace,˙Hs(¯Rn+)={u∈Hs(Rn)∣suppu⊂¯Rn+},thesupportedspace, | (5.1) |
as in our earlier papers on fractional-order operators. An elaborate presentation of Lp-based spaces was given in [8]. (The notation with dots and overlines stems from Hörmander [17,App. B.2] and is practical in formulas where both types of spaces occur. There are other notations without the overline, and where the dot is replaced by a ring or twiddle.)
Here ¯Hs(Rn+) identifies with the dual space of ˙H−s(¯Rn+) for all s∈R (the duality extending the L2(Rn+) scalar product). When |s|<12, there is an identification of ˙Hs(¯Rn+) with ¯Hs(Rn+) (more precisely with e+¯Hs(Rn+)). The trace operator γ0:u↦limxn→0+u(x′,xn) extends to a continuous mapping γ0:¯Hs(Rn+)→Hs−12(Rn−1) for s>12.
The order-reducing operators Ξt± are defined for t∈C by Ξt±=Op(χt±), where χt±=(⟨ξ′⟩±iξn)t, cf. (3.11). These operators have the homeomorphism properties:
Ξt+:˙Hs(¯Rn+)∼→˙Hs−Ret(¯Rn+),r+Ξt−e+:¯Hs(Rn+)∼→¯Hs−Ret(Rn+),alls∈R,t∈C; | (5.2) |
r+Ξt−e+ is often denoted Ξt−,+ for short. For each t∈C, the operators Ξt+ and Ξ¯t−,+ identify with each other's adjoints over ¯Rn+ (more comments on this in [8,Rem. 1.1]). Recall also the simple composition rules (as noted e.g., in [15,Th. 1.2]):
Ξs+Ξt+=Ξs+t+,Ξs−,+Ξt−,+=Ξs+t−,+fors,t∈C. |
We define
Eμ(¯Rn+)=e+xμnC∞(¯Rn+)whenReμ>−1, | (5.3) |
and Eμ(¯Rn+) is defined successively as the linear hull of first-order derivatives of elements of Eμ+1(¯Rn+) when Reμ≤−1 (then distributions supported in the boundary can occur). The spaces were introduced in Hörmander's unpublished lecture notes [16] and are presented in [8] (and with a different notation in [17,Sect. 18.2]), and they satisfy for all μ (cf. [8,Props. 1.7,4.1]):
Eμ(¯Rn+)∩E′(Rn)⊂Ξ−μ+e+⋂s¯Hs(¯Rn+). | (5.4) |
A sharper statement follows from [13,Lemma 6.1] (when Reμ>−1):
e+xμnS(¯Rn+)=Ξ−μ+e+S(¯Rn+). | (5.5) |
Let P satisfy Assumption 3.1, and consider ˆQ±=Op(ˆq±), defined from the symbols q±(ξ) introduced in Proposition 4.3. Since ˆq± are bounded symbols with bounded inverses, and extend holomorphically in ξn into C− resp. C+,
ˆQ+:˙Hs(¯Rn+)∼→˙Hs(¯Rn+)andˆQ−+=r+ˆQ−e+:¯Hs(Rn+)∼→¯Hs(Rn+),foralls∈R; | (5.6) |
the latter follows since r+ˆQ−e+ is the adjoint of Op(¯ˆq−) over Rn+, where Op(¯ˆq−) defines homeomorphisms in ˙Hs(¯Rn+) (since ¯ˆq− has similar properties as ˆq+). The inverses (ˆQ±)−1=Op((ˆq±)−1) have similar homeomorphism properties. Since ¯Hs(Rn+)=˙Hs(¯Rn+) for |s|<12, it follows that we also have for |s|<12:
ˆQ++=r+ˆQ+e+:¯Hs(Rn+)∼→¯Hs(Rn+),ˆQ−+ˆQ++:¯Hs(Rn+)∼→¯Hs(Rn+). |
If q satisfies the full 0-transmission condition, we are in the case studied in [8], and the bijectiveness in ¯Hs(Rn+) can be lifted to all higher s by use of elements of the Boutet de Monvel calculus, as accounted for in the proof of [8,Th. 4.4]. The symbol q presently considered is only known to satisfy the principal 0-transmission condition (and possibly a few more identities). We shall here show that a lifting is possible in general up to s<32.
Proposition 5.1. Let P satisfy Assumption 3.1, and consider ˆQ+=Op(ˆq+) derived from it in Section 4.
For any −12<s<32, ˆQ++=r+ˆQ+e+ is continuous
r+ˆQ+e+:¯Hs(Rn+)→¯Hs(Rn+), | (5.7) |
and the same holds for the operator ((ˆQ+)−1)+ defined from its inverse (ˆQ+)−1.
In fact, (5.7) is a homeomorphism, and the inverse of ˆQ++ is ((ˆQ+)−1)+.
Proof. We already have the mapping property (5.7) for |s|<12, because ˆq+ is a bounded symbol, and e+¯Hs(Rn+) identifies with ˙Hs(¯Rn+) then. Now let s=32−ε for a small ε>0. Here we need to show that when u∈¯H32−ε(Rn+), then r+∂jˆQ+e+u∈¯H12−ε(Rn+) for j=1,…,n. For j<n, this follows simply because ∂j can be commuted through r+, ˆQ+ and e+ so that we can use that ∂ju∈¯H12−ε(Rn+). For j=n, we proceed as follows:
Since u∈¯H32−ε(Rn+), the extension by zero e+u has a jump at xn=0, and a rule for distributions applies:
∂ne+u=e+∂nu+(γ0u)(x′)⊗δ(xn),γ0u∈H1−ε(Rn−1). | (5.8) |
(The rule is obvious when u∈C∞(¯Rn+), and extends by continuity to Sobolev spaces.) Therefore, since ˆQ+=I+ˆG+ where ˆG+=Op(ˆg+(ξ)) from Proposition 4.3,
∂nˆQ+e+u=ˆQ+∂ne+u=ˆQ+e+∂nu+(I+ˆG+)(γ0u⊗δ(xn)). |
In the restriction to Rn+, r+I(γ0u⊗δ(xn)) drops out, so we are left with
r+∂nˆQ+e+u=r+ˆQ+∂ne+u=r+ˆQ+e+∂nu+Kˆg+γ0u,Kˆg+φ=r+ˆG+(φ(x′)⊗δ(xn)). |
Here Kˆg+ is a potential operator (in the terminology of Eskin [6] and Rempel-Schulze [18], generalizing the concept of Poisson operator of Boutet de Monvel [3,4]), which acts as follows:
Kˆg+φ=r+F−1[ˆg+(ξ)ˆφ(ξ′)]. |
By (4.8),
|ˆg+(ξ)|≤C⟨ξ′⟩1−ε/2⟨ξ⟩ε/2−1, |
hence
‖Kˆg+φ‖2¯H12−ε(Rn+)≤‖ˆG+(φ⊗δ)‖2H12−ε(Rn)=c∫Rn|ˆg(ξ)|2|ˆφ(ξ′)|2⟨ξ⟩1−2εdξ≤C∫Rn|ˆφ(ξ′)|2⟨ξ⟩1−2ε−2+ε⟨ξ′⟩2−εdξ=C∫Rn|ˆφ(ξ′)|2⟨ξ⟩−1−ε⟨ξ′⟩2−εdξ=C′∫Rn−1|ˆφ(ξ′)|2⟨ξ′⟩2−2εdξ′=C"‖φ‖2H1−ε(Rn−1), |
since ∫R⟨ξ⟩−1−εdξn=⟨ξ′⟩−ε∫R⟨ηn⟩−1−εdηn. Inserting φ=γ0u, we thus have
‖Kˆg+γ0u‖¯H12−ε(Rn+)≤C1‖γ0u‖H1−ε(Rn−1)≤C2‖u‖¯H32−ε(Rn+). |
Thus
‖r+∂nˆQ+e+u‖¯H12−ε≤‖r+ˆQ+e+∂nu‖¯H12−ε+‖Kˆg+γ0u‖¯H12−ε≤C3‖u‖¯H32−ε. |
Altogether, this shows the desired mapping property for s=32−ε, and the property for general 12≤s<32 follows by interpolation with the case s=0.
The mapping property (5.7) holds for the inverse (ˆQ+)−1, since its symbol (q+)−1 equals 1+∑k≥1(−b+)k with essentially the same structure.
The identity ((ˆQ+)−1)+ˆQ++=I=ˆQ++((ˆQ+)−1)+ valid on L2(Rn+), holds a fortiori on ¯Hs(Rn+) for 0<s<32, and extends by continuity to ¯Hs(Rn+) for −12<s<0.
When P merely satisfies Assumption 3.2, we can still show a useful forward mapping property of ˆQ, based on the decomposition in Proposition 4.2.
Proposition 5.2. Let P satisfy Assumption 3.2, and consider ˆQ and ˆF±=Op(ˆf±) derived from it in Section 4.
The operator ˆF+,+=r+ˆF+e+ is continuous
r+ˆF+e+:¯Hs(Rn+)→¯Hs(Rn+)forany−12<s<32. | (5.9) |
The operator ˆF−,+=r+ˆF−e+ is continuous from ¯Hs(Rn+) to ¯Hs(Rn+) for any s∈R.
The operator ˆQ+=r+ˆQe+ is continuous
r+ˆQe+:¯Hs(Rn+)→¯Hs(Rn+)forany−12<s<32. | (5.10) |
Proof. Since ˆF+ has bounded symbol, it maps ˙Hs(¯Rn+) into Hs(Rn) for all s, so for |s|<12, (5.9) follows since ˙Hs(¯Rn+)=e+¯Hs(Rn+) then. For 12<s<32, we proceed as in the proof of Proposition 5.1, using that
r+∂nˆF+e+u=r+ˆF+∂ne+u=r+ˆF+e+∂nu+Kˆf+γ0u,Kˆf+φ=r+ˆF+(φ(x′)⊗δ(xn)), |
where Kˆf+ satisfies similar estimates as Kˆg+ by Proposition 4.2.
For r+ˆF−e+, the statement follows since it is on Rn+ the adjoint of Op(¯ˆf−), which preserves support in ¯Rn+ and therefore maps ˙Hs(¯Rn+) into itself for all s∈R. For ˆQ, the statement now follows since it equals s0+ˆF−+ˆF+.
This is as far as we get by applying Lemma 4.1 to f. To obtain the mapping property for higher s would require a control over the potential operators
φ↦r+Op(ξjnˆf+(ξ))(φ(x′)⊗δ(xn)) |
for j≥1 as well. At any rate, the property shown in Proposition 5.2 will be sufficient for the integration by parts formulas we are aiming for.
In the elliptic case, we conclude from Proposition 5.1 for the operator ˆQ:
Corollary 5.3. Let P satisfy Assumption 3.1, and consider the operators ˆQ, ˆQ+, ˆQ− with symbols ˆq, ˆq+, ˆq− derived from it in Section 4. The operator ˆQ+≡r+ˆQe+ acts like r+ˆQ−e+r+ˆQ+e+=ˆQ−+ˆQ++, mapping continuously and bijectively
ˆQ+=r+ˆQe+:¯Hs(Rn+)∼→¯Hs(Rn+)for−12<s<32, | (5.11) |
and the inverse (continuous in the opposite direction) equals
(r+ˆQe+)−1=r+(ˆQ+)−1e+r+(ˆQ−)−1e+. | (5.12) |
Proof. We have for u∈˙Hs(¯Rn+)≃¯Hs(Rn+), |s|<12, that
r+ˆQe+u=r+ˆQ−ˆQ+e+u=r+ˆQ−(e+r++e−r−)ˆQ+e+u=r+ˆQ−e+r+ˆQ+e+u, |
since r−ˆQ+e+u=0; this identity is also valid on the subspaces ¯Hs(Rn+) with s≥12. Combining the homeomorphism property of r+ˆQ+e+ shown in Proposition 5.1 with the known homeomorphism property of r+ˆQ−e+ on ¯Hs(Rn+)-spaces (cf. (5.6)), we get (5.11). The inverse is pinned down by using that r+ˆQ−e+ has inverse r+(ˆQ−)−1e+ on ¯Hs(Rn+) for all s, and r+ˆQ+e+ has inverse r+(ˆQ+)−1e+ on ¯Hs(Rn+) for −12<s<32 in view of Proposition 5.1.
Now turn the attention to ˆP, which is related to ˆQ by
ˆP=Ξμ′−ˆQΞμ+,ˆQ=Ξ−μ′−ˆPΞ−μ+, | (5.13) |
cf. (3.12)–(3.13).
We shall describe the solutions of the homogeneous Dirichlet problem (in the strongly elliptic case)
r+ˆPu=f,suppu⊂¯Rn+, | (5.14) |
with f given in a space ¯Hs(Rn+), and u assumed a priori to lie in a space ˙Hσ(¯Ω) for low σ, e.g., with σ=a.
First we observe for Ξ−μ′−,+=r+Ξ−μ′−e+ that
Ξ−μ′−,+r+ˆP=r+Ξ−μ′−ˆP, | (5.15) |
since, as accounted for in [8,Rem. 1.1,(1.13)], the action of r+Ξ−μ′− is independent of how r+ˆP is extended into ¯Rn−. Thus, in view of the mapping properties (5.2) of Ξ−μ′−,+,
‖r+ˆPu‖¯Hs(Rn+)≃‖Ξ−μ′−,+r+ˆPu‖¯Hs+Reμ′(Rn+)=‖r+Ξ−μ′−ˆPu‖¯Hs+Reμ′(Rn+). | (5.16) |
Composing the equation in (5.14) with Ξ−μ′−,+ to the left, we can therefore write it as
r+Ξ−μ′−Pu=g,whereg=Ξ−μ′−,+f∈¯Hs+Reμ′(Rn+). | (5.17) |
Next, we shall also replace u. Because of the right-hand factor Ξ−μ+ in the expression for ˆQ in (5.13), we need to introduce the μ-transmission spaces
Hμ(t)(¯Rn+)≡Ξ−μ+e+¯Ht−Reμ(Rn+)fort>Reμ−12, | (5.18) |
defined in [8]; they are Hilbert spaces. (For t≤Reμ−12, the convention is to take Hμ(t)(¯Rn+)=˙Ht(¯Rn+), but this is rarely used.) The following properties were shown in [8]:
Theorem 5.4. [8] Let t>Reμ−12.
1∘ The mapping r+Ξμ+ is a homeomorphism of Hμ(t)(¯Rn+) onto ¯Ht−Reμ(Rn+) with inverse Ξ−μ+e+.
2∘ For |t−Reμ|<12, Hμ(t)(¯Rn+)=˙Ht(¯Rn+).
3∘ Assume Reμ>−1 and t>Reμ+12. Then
Hμ(t)(¯Rn+)⊂˙Ht(¯Rn+)+xμne+¯Ht−Reμ(Rn+), | (5.19) |
where ˙Ht(¯Rn+) is replaced by ˙Ht−ε(¯Rn+) if t−Reμ−12∈N. Moreover, the trace of u/xμn is well-defined on Hμ(t)(¯Rn+) and satisfies
Γ(1+μ)γ0(u/xμn)=γ0Ξμ+u∈Ht−Reμ−12(Rn−1). | (5.20) |
Rule 1∘ is shown in [8,Prop. 1.7]. Rule 2∘, shown in [8,(1.26)], holds because of the mapping property (5.2) for Ξ−μ+ and the identification of e+¯Ht−Reμ(Rn+) with ˙Ht−Reμ(¯Rn+) when t−Reμ∈]−12,12[. Rule 3∘ is shown in [8,Th. 5.1,Cor. 5.3,Th. 5.4]; it deals with a higher t, where e+¯Ht−Reμ(Rn+) has a jump at xn=0, and the coefficient xμn appears. Let us just mention the key formula
F−1ξn→xn[(⟨ξ′⟩+iξn)−μ(⟨ξ′⟩+iξn)−1]=1Γ(μ+1)e+r+xμne−⟨ξ′⟩xn, |
which indicates how Ξ−μ+=Op((⟨ξ′⟩+iξn)−μ) is connected with the factor xμn. Besides in [8,Sect. 5], explicit calculations are carried out e.g., in [12,Lemma 3.3] (and [9,Appendix]).
We note in passing that in the definition (5.18), one can equivalently replace the order-reducing operator family Ξt+=Op((⟨ξ′⟩+iξn)t) by Op(([ξ′]+iξn)t), or by Λt+, as defined in [8].
Now continue the discussion of (5.17): In view of Theorem 5.4 1∘, we can set v=r+Ξμ+u, where u=Ξ−μ+e+v, and hereby
r+Ξ−μ′ˆPu=r+Ξ−μ′ˆPΞ−μ+e+v=r+ˆQe+v=ˆQ+v. |
Then the Eq (5.17) reduces to an equivalent equation
ˆQ+v=g, |
with g given in ¯Hs+Reμ′(Rn+) and v a priori taken in ˙Hσ−Reμ(¯Rn+). We shall denote s=t−2a, so s+Reμ′=t−2a+2a−Reμ=t−Reμ. The equation was solved in Corollary 5.3 and we find for r+ˆP:
Theorem 5.5. Let P satisfy Assumption 3.1. For Reμ−12<t<Reμ+32, r+ˆP defines a homeomorphism (continuous bijective operator with continuous inverse)
r+ˆP:Hμ(t)(¯Rn+)∼→¯Ht−2a(Rn+). | (5.21) |
Furthermore, if u is in ˙Hσ(¯Rn+) for some σ>Reμ−12 (this includes the value σ=a) and solves (5.14) with f∈¯Ht−2a(Rn+), then u∈Hμ(t)(¯Rn+).
Here Reμ>−12 since a>0 and |Reδ|<12, so the rules in Theorem 5.4 3∘ apply.
Proof. In view of (5.15) and (5.16), and the mapping property of ˆQ established in Corollary 5.3, r+ˆP has the forward mapping property in (5.21).
To solve (5.14), let σ=Reμ−12+ε for a small ε, set g=Ξ−μ′−,+f∈¯Ht−2a+Reμ′(Rn+)=¯Ht−Reμ(Rn+) and v=r+Ξμ+u∈˙H−12+ε(¯Rn+). Then (5.14) reduces to solving
ˆQ+v=g, | (5.22) |
with g given in ¯Ht−Reμ(Rn+) and v a priori lying in ˙H−12+ε(¯Rn+). By Corollary 5.3, (5.22) has a unique solution v∈¯Ht−Reμ(Rn+), so u must lie in Hμ(t)(¯Rn+), and the mapping (5.21) is bijective.
Remark 5.6.This theorem differs from the strategy pursued in [6], and gives a new insight. The technique in [6,Th. 7.3] for showing solvability in a higher-order Sobolev space, say with 12<t−Reμ<32, f given in ¯Ht−2a(Rn+), is to supplement ˆP with a potential operator KˆP constructed from ˆP such that the solutions are of the form u=u++KˆPφ with u+∈˙Ht(¯Rn+), φ a generalized trace derived from f. Our aim is to show that there is a universal description of the space of solutions u of (5.14) with right-hand side in ¯Ht−2a(Rn+), that depends only on μ, and applies to any P of the given type. The μ-transmission spaces (5.18) serve this purpose. In [8], they are shown to have this role for arbitrarily high t when the full μ-transmission condition holds. $
One more important property of μ-transmision spaces is that the spaces with C∞-ingredients Eμ(¯Rn+)∩E′(Rn) and e+xμnS(¯Rn+) are dense subsets of Hμ(t)(¯Rn+) for all t>Reμ−12, Reμ>−1 (cf. [8,Prop. 4.1] and [13,Lemma 7.1]). Recall also (5.5), which makes the statement for e+xμnS(¯Rn+) rather evident, since S(¯Rn+) is dense in ¯Hs(Rn+) for all s∈R. Hence r+ˆP applies nicely to these spaces.
When P merely satisfies Assumption 3.2, we have at least the forward mapping part of (5.21):
Theorem 5.7. Let P satisfy Assumption 3.2. For Reμ−12<t<Reμ+32, r+ˆP maps continuously
r+ˆP:Hμ(t)(¯Rn+)→¯Ht−2a(Rn+). | (5.23) |
Proof. This follows as in the preceding proof, now using the mapping property of r+ˆQe+ established in Proposition 5.2.
The following consequences can be drawn for the original operator P:
Theorem 5.8. 1∘ Let P satisfy Assumption 3.2. Then P=ˆP+P′, where ˆP is defined by (3.9) and P′ is of order 2a−1. For Reμ−12<t<Reμ+32, r+P maps continuously
r+P:Hμ(t)(¯Rn+)→¯Ht−2a(Rn+). | (5.24) |
2∘ Let P satisfy Assumption 3.1. Then in the decomposition P=ˆP+P′, r+ˆP is invertible, as described in Theorem 5.5.
Let Reμ−12<t<Reμ+32, let f∈¯Ht−2a(Rn+), and let u∈˙Hσ(¯Rn+) (for some σ>Reμ−12) solve the homogeneous Dirichlet problem
r+Pu=fonRn+,suppu⊂¯Rn+. | (5.25) |
Then u∈Hμ(t)(¯Rn+).
Proof. The original operator P equals Op(p(ξ)) with p(ξ) homogeneous on Rn of degree 2a>0; in particular it is continuous at 0. It is decomposed into
p(ξ)=ˆp(ξ)+p′(ξ). | (5.26) |
where p′(ξ) is O(⟨ξ⟩2a−1) for |ξ|≥2 by (3.10) and continuous, hence
|p′(ξ)|≤C′⟨ξ⟩2a−1forξ∈Rn. |
This implies that P′=Op(p′) maps Hs(Rn) continuously into Hs−2a+1(Rn) for all s∈R, and hence
r+P′:˙Hs(¯Rn+)→¯Hs−2a+1(Rn+)foralls∈R. | (5.27) |
1∘. The forward mapping property (5.23) holds for r+ˆP by Theorem 5.7. To show that it holds for r+P′, let Reμ−12<t<Reμ+32.
If t−Reμ<12, then Hμ(t)(¯Rn+)=˙Ht(¯Rn+), and r+P′˙Ht(¯Rn+)⊂¯Ht−2a+1(Rn+)⊂¯Ht−2a(Rn+) by (5.27), matching the mapping property of ˆP.
If 12≤t−Reμ<32, we use the definition of Hμ(t)(¯Rn+) to see that for small ε>0,
r+P′Hμ(t)(¯Rn+)=r+P′Ξ−μ+e+¯Ht−Reμ(Rn+)⊂r+P′Ξ−μ+˙H12−ε(¯Rn+)=r+P′˙H12−ε+Reμ(¯Rn+)⊂¯H32−ε+Reμ−2a(¯Rn+)⊂¯Ht−2a(Rn+), |
also matching the mapping property of ˆP.
Now (5.24) follows by adding the statements for P′ and ˆP. This shows 1∘.
2∘. The first statement registers what we already know about r+ˆP. Proof of the regularity statement: With u and f as defined there, denote σ=Reμ−12+ε; here ε>0. Then
r+ˆPu=r+Pu−r+P′u∈¯Ht−2a(Rn+)+¯HReμ+12+ε−2a(Rn+). |
If t≤Reμ+12+ε, r+ˆPu∈¯Ht−2a(Rn+), and we conclude from Theorem 5.5 that u∈Hμ(t)(¯Rn+).
If t>Reμ+12+ε, r+ˆPu∈¯HReμ+12+ε−2a(Rn+); here Theorem 5.5 applies to give the intermediate information that u∈Hμ(Reμ+12+ε)(¯Rn+). From this follows that
u∈Ξ−μ+e+¯H12+ε(Rn+)⊂Ξ−μ+˙H12−ε′(¯Rn+)=˙HReμ+12−ε′(¯Rn+), |
for any ε′>0. Then r+P′u∈¯HReμ+32−ε′−2a(Rn+). Choosing ε′ so small that Reμ+32−ε′≥t, we have that r+P′u∈¯Ht−2a(¯Rn+); hence r+ˆPu∈¯Ht−2a(¯Rn+), so it follows from Theorem 5.5 that u∈Hμ(t)(¯Rn+). This ends the proof of 2∘.
Example 5.9. Theorem 5.8 applies to the operator L=Op(L(ξ)) described in (3.5)ff., showing that it maps Hμ(t)(¯Rn+) to ¯Ht−2a(Rn+) for −12<t−μ<32, and that solutions of the homogeneous Dirichlet problem with f∈¯Ht−2a(Rn+) are in Hμ(t)(¯Rn+) for these t. The appearance of the factor xμn (cf. (5.19)) is consistent with the regularity shown in terms of Hölder spaces in [5].
In particular, the result provides a valid basis for applying r+L to Eμ(¯Rn+)∩E′(Rn) or e+xμnS(¯Rn+), mapping these spaces into ⋂ε>0¯H32−a+δ−ε(Rn+). $
Remark 5.10. The domain spaces Hμ(t)(¯Rn+) entering in Theorem 5.8 can be precisely described: For |t−Reμ|<12, we already know from Theorem 5.4 2∘ that Hμ(t)(¯Rn+)=˙Ht(¯Rn+). For 12<t−Reμ<32, we have by [12,Lemma 3.3] that u∈Hμ(t)(¯Rn+) if and only if
u=v+w, where w∈˙Ht(¯Rn+) and v=e+xμnK0γ0(u/xμn); |
here K0 is the Poisson operator K0:φ↦z solving the Dirichlet problem for 1−Δ,
(1−Δ)z=0 on Rn+,γ0z=φ at xn=0, |
with φ∈Ht−Reμ−12(Rn−1). For t−Reμ=12, we have the information u∈⋂ε>0˙Ht−ε(¯Rn+). As a concrete example, the elements u of H12(32)(¯Rn+) are the functions u=v+w, where w∈˙H32(¯Rn+) and v=x12nK0φ for some φ∈H12(Rn−1); this φ equals γ0(u/x12n). $
It will now be shown that the operators P satisfying merely the principal μ-transmission condition (Assumption 3.2) have an integration by parts formula over Rn+, involving traces γ0(u/xμn). The study will cover the special operator L in Example 5.9 (regardless of whether a full μ-transmission condition might hold, as assumed in [13]). It also covers more general strongly elliptic operators, and it covers operators that are not necessarily elliptic.
The basic observation is:
Proposition 6.1. Let μ∈C. Let w∈⋂s¯Hs(Rn+), and let u′∈Eˉμ(¯Rn+)∩E′(Rn). Denote w′=r+Ξˉμ+u′∈⋂s¯Hs(Rn+); correspondingly u′=Ξ−ˉμ+e+w′ in view of Theorem 5.4 1∘. Then
(I≡)∫Rn+Ξμ−e+w∂nˉu′dx=(γ0w,γ0w′)L2(Rn−1)+(w,∂nw′)L2(Rn+). | (6.1) |
The left-hand side is interpreted as in (6.2) below when Reμ≤0.
The formula extends to w∈¯H12+ε(Rn+) and u′∈Hˉμ(t)(¯Rn+) with t≥Reμ+12−ε (for small ε>0), using the representation (6.2).
Proof. This was proved in [10,Th. 3.1] for μ=a>0 (see also Remark 3.2 there with the elementary case a=1), and in [13,Th. 4.1] for real μ>−12, so the main task is to check that the larger range of complex μ is allowed. We write ¯u′ as ˉu′ for short.
Note that when Imμ≠0, Eμ(¯Rn+) is different from EReμ(¯Rn+), e.g., since xμn/xReμn=xiImμn=eiImμlogxn has absolute value 1 and is C∞ for xn>0, but oscillates when xn→0.
By the mapping properties of Ξμ−,+ (cf. (5.2)), r+Ξμ−e+w∈⋂s¯Hs(Rn+), hence is integrable. When Reμ>0, the function ∂nu′ is O(xReμ−1n) and compactly supported, so the left-hand side of (6.1) makes sense as an integral of an L1-function. When μ is general, we observe that for any small ε>0,
∂nu′∈Eˉμ−1(¯Rn+)∩E′(Rn)⊂Ξ1−ˉμ+e+⋂s¯Hs(Rn+)⊂Ξ1−ˉμ+˙H12−ε(¯Rn+)=˙H−12+Reμ−ε(¯Rn+), |
so the integral I makes sense as the duality
I=⟨r+Ξμ−e+w,∂nu′⟩¯H12−Reμ+ε(Rn+),˙H−12+Reμ−ε(¯Rn+). | (6.2) |
Since the adjoint of r+Ξμ−e+ equals Ξˉμ+, I is by transposition turned into
I=⟨w,Ξˉμ+∂nu′⟩¯H12+ε,˙H−12−ε=⟨w,∂nΞˉμ+Ξ−ˉμ+e+w′⟩¯H12+ε,˙H−12−ε=⟨w,∂ne+w′⟩¯H12+ε,˙H−12−ε. |
Note that ∂ne+w satisfies an equation like (5.8), which fits in here since the space ˙H−12−ε(¯Rn+) contains distributions of the form φ(x′)⊗δ(xn). The expression is analysed as in [10,Th. 3.1] (and [13,Th. 4.1]), leading to
I=(γ0w,γ0w′)L2(Rn−1)+(w,e+∂nw′)L2(Rn+), | (6.3) |
which shows (6.1).
For the whole analysis, it suffices that w∈¯Hs(Rn+) with s=12+ε, since Ξμ−,+w∈¯H12−Reμ+ε(Rn+) then. For u′, it then suffices that u′∈Hˉμ(t)(¯Rn+)=Ξ−ˉμ+e+¯Ht−Reμ(Rn+) with t≥Reμ+12−ε (assuming 0<ε<1), since
∂nu′∈Ξ1−ˉμ+e+¯Ht−Reμ(Rn+)⊂Ξ1−ˉμ+˙H12−ε(¯Rn+)=˙H−12+Reμ−ε(¯Rn+) |
then, so that the duality in (6.2) is well-defined.
We shall now show:
Theorem 6.2. Let P satisfy Assumption 3.2; it is of order 2a and satisfies the principal μ-transmission condition in the direction (0,1) for some μ=a+δ∈C, and we denote a−δ=μ′. Assume moreover that Reμ>−1, Reμ′>−1. Consider ˆP=Op(ˆp(ξ)), as defined by (3.9). For u∈Eμ(¯Rn+)∩E′(Rn), u′∈Eˉμ′(¯Rn+)∩E′(Rn), there holds
∫Rn+ˆPu∂nˉu′dx+∫Rn+∂nu¯ˆP∗u′dx=Γ(μ+1)Γ(μ′+1)∫Rn−1s0γ0(u/xμn)γ0(ˉu′/xμ′n)dx′, | (6.4) |
where s0=e−iπδp(0,1). The formula extends to u∈Hμ(t)(¯Rn+), u′∈Hˉμ′(t′)(¯Rn+), for t>Reμ+12, t′>Reμ′+12.
The integrals over Rn+ are interpreted as dualities (as in Proposition 6.1) when Reμ or Reμ′≤0, and when extended to general u,u′.
Proof. Since integration over Rn+ in itself indicates that the functions behind the integration sign are restricted to Rn+, we can leave out the explicit mention of r+. Recall that
ˆp=χμ′−ˆqχμ+,ˆP=Ξμ′−ˆQΞμ+, |
cf. (3.13). The adjoint is ˆP∗=Ξˉμ−ˆQΞˉμ′+. Recall from Proposition 4.2 that
q(ξ)=s0+f+(ξ)+f−(ξ),henceˆQ=s0+ˆF++ˆF−, |
where ˆf±(ξ) extend holomorphically in ξn+iτ into C− resp. C+, estimated as in (4.8).
Accordingly, ˆP splits up in three terms
ˆP=ˆP1+ˆP2+ˆP3,whereˆP1=s0Ξμ′−Ξμ+,ˆP2=Ξμ′−ˆF+Ξμ+,ˆP3=Ξμ′−ˆF−Ξμ+. | (6.5) |
Consider the contribution from ˆP1:
∫Rn+ˆP1u∂nˉu′dx+∫Rn+∂nu¯ˆP∗1u′dx=s0∫Rn+Ξμ′−Ξμ+u∂nˉu′dx+s0∫Rn+∂nu¯Ξˉμ−Ξˉμ′+u′dx. |
Recall that s0=q(0,1)=e−iπδp(0,1) by (3.16); this constant is left out of the next calculations.
When u∈Eμ(¯Rn+)∩E′(Rn), then w=r+Ξμ+u∈⋂s¯Hs(Rn+). Similarly as in (5.15), r+Ξμ′−Ξμ+u=r+Ξμ′−e+r+Ξμ+u, which equals r+Ξμ′−e+w, hence lies in ⋂s¯Hs(Rn+) by (5.2). An application of Proposition 6.1 with μ replaced by μ′ gives:
∫Rn+Ξμ′−Ξμ+u∂nˉu′dx=∫Rn+r+Ξμ′−e+w∂nˉu′dx=(γ0w,γ0w′)L2(Rn−1)+(w,∂nw′)L2(Rn+), |
where w′=r+Ξˉμ′+u′.
We can apply the analogous argument to show that the conjugate of ∫Rn+∂nu¯Ξˉμ−Ξˉμ′+u′dx satisfies
∫Rn+Ξˉμ−Ξˉμ′+u′∂nˉudx=(γ0w′,γ0w)L2(Rn−1)+(w′,∂nw)L2(Rn+); |
here w′=r+Ξˉμ′+u′ and w=r+Ξμ+u are the same as the functions defined in the treatment of the first integral.
It follows by addition that
∫Rn+Ξμ′−Ξμ+u∂nˉu′dx+∫Rn+∂nu¯Ξˉμ−Ξˉμ′+u′dx=2(γ0w,γ0w′)L2(Rn−1)+(w,∂nw′)L2(Rn+)+(∂nw,w′)L2(Rn+)=(γ0w,γ0w′)L2(Rn−1); |
in the last step we used that ∫Rn+(w∂nˉw′+∂nwˉw′)dx=−∫Rn−1γ0wγ0ˉw′dx′. Insertion of γ0w=γ0Ξμ+u=Γ(1+μ)γ0(u/xμn) (cf. (5.20)), and similarly γ0w′=γ0Ξˉμ′+u′=Γ(1+ˉμ′)γ0(u′/xˉμ′n), gives (6.4) with ˆP replaced by ˆP1 (using also that ¯Γ(1+ˉμ′)=Γ(1+μ′)).
As for extension of the formula to larger spaces, we note that by Proposition 6.1, the calculations for the first integral allow w∈¯H12+ε(Rn+), corresponding to u∈Hμ(t)(¯Rn+) with t=Reμ+12+ε, and u′∈Hˉμ′(t′)(¯Rn+) with t′≥Reμ′+12−ε. With the analogous conditions for the calculations of the second integral, we find altogether that t>Reμ+12, t′>Reμ′+12, is allowed.
The contributions from ˆP2 and ˆP3 will be treated by variants of this proof, where we show that their boundary integrals give zero.
Consider ˆP2. As in (5.15), we have:
r+ˆP2u=r+Ξμ′−ˆF+Ξμ+u=r+Ξμ′−e+r+(ˆF+Ξμ+u),r+ˆP∗2u′=r+Ξˉμ−ˆF∗+Ξˉμ′+u′=r+(Ξˉμ−ˆF∗+)e+r+Ξˉμ′+u′, |
where ˆF∗+=Op(¯ˆf+). Set
w=r+Ξμ+u,w1=r+ˆF+Ξμ+u,w′=r+Ξˉμ′+u′. | (6.6) |
Here when u∈Hμ(t)(¯Rn+), w∈¯Ht−Reμ(Rn+), and when u′∈Hˉμ′(t′)(¯Rn+), w′∈¯Ht′−Reμ′(Rn+). For w1 we have since u=Ξ−μ+e+w (by Theorem 5.4 1∘), that
w1=r+ˆF+Ξμ+u=r+ˆF+Ξμ+Ξ−μ+e+w=r+ˆF+e+w∈¯Ht−Reμ(Rn+), |
when −12<t−Reμ<32, by the mapping property for ˆF+ established in Proposition 5.2.
We can then apply Proposition 6.1 to the first integral for ˆP2, with μ replaced by μ′, giving when t−Reμ>12:
∫Rn+ˆP2u∂nˉu′dx=∫Rn+Ξμ′−ˆF+Ξμ+u∂nˉu′dx=∫Rn+Ξμ′−e+w1∂nˉu′dx=(γ0w1,γ0w′)L2(Rn−1)+(w1,∂nw′)L2(Rn+). | (6.7) |
There is a general formula for the trace, entering in Vishik and Eskin's calculus as well as that of Boutet de Monvel,
γ0v=(2π)−n∫Rn−1eix′⋅ξ′∫RF(e+v)dξndξ′, |
where the integral over R is read either as an ordinary integral or, if necessary, as the integral ∫+ defined e.g., in [7,(10.85)] (also recalled in [11,(A.1),(A.15)]). Applying this to w1, we find:
γ0w1=γ0(ˆF+,+w)=(2π)−n∫Rn−1eix′⋅ξ′∫Rˆf+(ξ′,ξn)F(e+w)dξndξ′. | (6.8) |
This integral gives 0 for the following reason: It suffices to take w in the dense subspace of ¯Ht−Reμ(Rn+) of compactly supported functions in C∞(¯Rn+). Both ˆf+ and F(e+w) are holomorphic in C− as functions of ξn, f+ being O(⟨ξn⟩−1+ε) and F(e+w) being O(⟨ξn⟩−1) on ¯C−, whereby the integrand is O(⟨ξn⟩−2+ε) there (and is in L1 on R); then the integral over R can be transformed to a closed contour in C− and gives 0.
We can then conclude:
∫Rn+ˆP2u∂nˉu′dx=∫Rn+Ξμ′−ˆF+Ξμ+u∂nˉu′dx=(w1,∂nw′)L2(Rn+). | (6.9) |
The other contribution from ˆP2 is, in conjugated form,
∫Rn+ˆP∗2u′∂nˉudx=∫Rn+Ξˉμ−ˆF∗+e+r+Ξˉμ′+u′∂nˉudx=⟨r+Ξˉμ−ˆF∗+e+r+Ξˉμ′+u′,∂nu⟩¯H12−Reμ+ε,˙H−12+Reμ−ε=⟨r+Ξˉμ′+u′,ˆF+Ξμ+∂nu⟩¯H12+ε,˙H−12−ε=⟨r+Ξˉμ′+u′,∂nˆF+Ξμ+u⟩¯H12+ε,˙H−12−ε=⟨r+Ξˉμ′+u′,∂nˆF+e+w⟩¯H12+ε,˙H−12−ε=⟨w′,∂ne+w1⟩¯H12+ε,˙H−12−ε=(γ0w′,γ0w1)L2(Rn−1)+(w′,∂nw1)L2(Rn+)=(w′,∂nw1)L2(Rn+), |
where we used Proposition 6.1 in a similar way, and at the end used that γ0w1=0, cf. (6.8)ff. Finally, taking the contributions from ˆP2 together, we get
∫Rn+ˆP2u∂nˉu′dx+∫Rn+∂nu¯ˆP∗2u′dx=(w1,∂nw′)L2(Rn+)+(∂nw1,w′)L2(Rn+)=−(γ0w1,γ0w′)L2(Rn−1)=0, |
using again that γ0w1=0.
It is found in a similar way, using that ˆF∗− is of plus-type, that ˆP3 contributes with zero.
To extend the formula to the original operator P, we shall show that P′=P−ˆP (cf. Theorem 5.8 1∘) gives a zero boundary contribution.
Lemma 6.3. Let a>0 and let S=Op(s(ξ)), where s(ξ) is O(⟨ξ⟩2a−1). Then
∫Rn+Su∂nˉu′dx+∫Rn+∂nu¯S∗u′dx=0, | (6.10) |
for any u,u′∈˙Ha(¯Rn+).
Proof. Since u∈˙Ha(¯Rn+), r+Su∈¯H1−a(Rn+); moreover ∂nu′∈˙Ha−1(¯Rn+), so we can write the first integral as
⟨r+Su,∂nu′⟩¯H1−a(Rn+),˙Ha−1(¯Rn+). |
Approximate u′ in ˙Ha(¯Rn+) by a sequence of functions φk∈C∞0(Rn+), k∈N; then
⟨r+Su,∂nφk⟩¯H1−a,˙Ha−1=−⟨r+∂nSu,φk⟩¯H−a,˙Ha→−⟨r+∂nSu,u′⟩¯H−a,˙Ha. |
With a similar argument for the second integral, we have
⟨r+Su,∂nu⟩¯H1−a,˙Ha−1+⟨∂nu,r+S∗u′⟩˙Ha−1,¯H1−a=−⟨r+∂nSu,u′⟩¯H−a,˙Ha−⟨u,r+∂nS∗u⟩˙Ha,¯H−a=−⟨r+∂nSu,u′⟩¯H−a,˙Ha+⟨u,r+(∂nS)∗u⟩˙Ha,¯H−a=0, |
since ∂nS∗=S∗∂n=−(∂nS)∗, and it is well-known that the operator S1=∂nS of order 2a satisfies ⟨r+S1u,u′⟩=⟨u,r+S∗1u′⟩ for u,u′∈˙Ha(¯Rn+).
We can then conclude:
Theorem 6.4. Let P, μ, μ′ be as in Theorem 6.2. For u∈Eμ(¯Rn+)∩E′(Rn), u′∈Eˉμ′(¯Rn+)∩E′(Rn), there holds
∫Rn+Pu∂nˉu′dx+∫Rn+∂nu¯P∗u′dx=Γ(μ+1)Γ(μ′+1)∫Rn−1s0γ0(u/xμn)γ0(ˉu′/xμ′n)dx′, | (6.11) |
where s0=e−iπδp(0,1). The formula extends to u∈Hμ(t)(¯Rn+), u′∈Hˉμ′(t′)(¯Rn+), for t>Reμ+12, t′>Reμ′+12, with t,t′≥a.
The integrals over Rn+ are interpreted as dualities (as in Proposition 6.1 and Lemma 6.3) when Reμ or Reμ′≤0, and when extended to general u,u′.
Proof. Recall that P=ˆP+P′, where P′=Op(p′), |p′(ξ)|≤C⟨ξ⟩2a−1 (cf. Theorem 5.8 1∘). We have the identities (6.10) with S=P′ and (6.4) for u∈Hμ(t)(¯Rn+), u′∈Hˉμ′(t′)(¯Rn+) with t,t′≥a, t>Reμ+12, t′>Reμ′+12. Adding the identities for ˆP and P′ we obtain (6.11). It holds a fortiori for u∈Eμ(¯Rn+)∩E′(Rn), u′∈Eˉμ′(¯Rn+)∩E′(Rn).
Example 6.5. The theorem applies in particular to L=Op(L(ξ)) studied in (3.5)–(3.6) and Example 5.9, showing that
∫Rn+Lu∂nˉu′dx+∫Rn+∂nu¯L∗u′dx=Γ(μ+1)Γ(μ′+1)∫Rn−1|L(0,1)|γ0(u/xμn)γ0(ˉu′/xμ′n)dx′, | (6.12) |
The value s0=|L(0,1)|=(A(0,1)2+B(0,1)2)12 is found in (3.17).
This result was proved in [5,Prop. 1.4] by completely different, real methods, for μ∈ ]0,2a[∩]2a−1,1[.
The result is one of the key ingredients in the proof of integration by parts formulas for operators L on bounded domains Ω⊂Rn in [5], where μ(ν) varies as the normal ν varies along the boundary. It would be interesting to extend this knowledge to general strongly elliptic operators P on bounded domains by similar applications of Theorem 6.4.
Example 6.6. Here is an example of an application to a nonelliptic operator satisfying Assumption 3.2. Let
P=|∂1+∂2|m+|∂3|m−1∂3, with symbol p(ξ)=|ξ1+ξ2|m+isignξ3|ξ3|m |
on R3, for some 1<m<2. For R3+={x3>0} we have the normal ν=(0,0,1), where
p(0,0,1)=i,p(0,0,−1)=−i, so (3.1) holds with m−2μ=1, |
i.e., μ=(m−1)/2, μ′=(m+1)/2. Then by Theorem 6.4,
∫x3>0(Puˉv−u¯P∗v)dx=Γ(m2+12)Γ(m2+32)∫x3=0γ0(ux(m−1)/23)γ0(ˉvx(m+1)/23)dx′, |
for functions u∈x(m−1)/23S(¯R3+), v∈x(m+1)/23S(¯R3+).
The halfspace {x2>0} has the normal ν′=(0,1,0) and
p(ν′)=1,p(−ν′)=1, so (3.1) holds with m−2μ=0, |
i.e., μ=m/2, μ′=m/2. Here by Theorem 6.4,
∫x2>0(Puˉv−u¯P∗v)dx=Γ(m2+1)2∫x2=0γ0(uxm/22)γ0(ˉvxm/22)dx1dx3, |
for functions with a factor xm/22. $
Let P satisfy Assumption 3.1, and assume Reμ>0. Along with the homogeneous Dirichlet problem (5.25), one can consider a nonhomogeneous local Dirichlet problem if the scope is expanded to allow so-called "large solutions", behaving like xμ−1n near the boundary of Rn+; such solutions blow up at the boundary when Reμ<1. Namely, one can pose the nonhomogeneous Dirichlet problem
r+Pu=fonRn+,γ0(u/xμ−1n)=φonRn−1,suppu⊂¯Rn+. | (7.1) |
Problem (7.1) was studied earlier for operators satisfying the a-transmission property in [8,9] (including the fractional Laplacian (−Δ)a), and a halfways Green's formula was shown in [11,Cor. 4.5]. The problem (7.1) for the fractional Laplacian, and the halfways Green's formula——with applications to solution formulas——were also studied in Abatangelo [1] (independently of [8]); the boundary condition there is given in a less explicit way except when Ω is a ball. There have been further studies of such problems, see e.g., Abatangelo, Gomez-Castro and Vazquez [2] and its references.
Note that the boundary condition in (7.1) is local. There is a different problem which is also regarded as a nonhomogeneous Dirichlet problem, namely to prescribe nonzero values of u in the exterior of Ω; it has somewhat different solution spaces (a link between this and the homogeneous Dirichlet problem is described in [9]).
For the general operators P considered here, we shall now show that problem (7.1) has a good sense for u∈H(μ−1)(t)(¯Rn+) with suitable t.
More precisely, since P also satisfies the principal (μ−1)-transmission condition (as remarked after Definition 2.1), Theorem 5.8 1∘ can be applied with μ replaced by μ−1, implying that r+P maps
r+P:H(μ−1)(t)(¯Rn+)→¯Ht−2a(Rn+)forReμ−32<t<Reμ+12. | (7.2) |
This is also valid in the case where P is only assumed to satisfy Assumption 3.2.
From Theorem 5.4 we have (note that Reμ−1>−1)
H(μ−1)(t)(¯Rn+){=˙Ht(¯Rn+)when−32<t−Reμ<−12,⊂˙Ht(¯Rn+)+xμ−1ne+¯Ht−Reμ+1(Rn+)when−12<t−Reμ<12. | (7.3) |
When t−Reμ>−12, the weighted boundary value is well-defined, cf. (5.20):
γμ−10u≡Γ(μ)γ0(u/xμ−1n)=γ0(Ξμ−1+u)∈Ht−Reμ+12(Rn−1). | (7.4) |
The following regularity result holds for the nonhomogeneous Dirichlet problem:
Theorem 7.1. Let P satisfy Assumption 3.1 with Reμ>0, and let −12<t−Reμ<12. When f∈¯Ht−2a(Rn+) and φ∈Ht−Reμ+12(Rn−1) are given, and u solves the nonhomogeneous Dirichlet problem (7.1) with u∈H(μ−1)(σ)(¯Rn+) for some −12<σ−Reμ<12, then in fact u∈H(μ−1)(t)(¯Rn+).
Proof. It is known from [8,Th. 6.1] that Hμ(σ)(¯Rn+) is a closed subspace of H(μ−1)(σ)(¯Rn+), equal to the set of v∈H(μ−1)(σ)(¯Rn+) for which γ0(v/xμ−1n)=0. From the given φ we define
w=Γ(μ)Ξ−μ+1+e+K0φ∈H(μ−1)(t)(¯Rn+), |
where K0 is the standard Poisson operator φ↦K0φ=F−1ξ′→x′[ˆφ(ξ′)e−⟨ξ′⟩xn], xn>0. Then in view of (7.4),
γ0(w/xμ−1n)=Γ(μ)−1γ0(Ξμ−1+w)=γ0(Ξμ−1+Ξ1−μ+e+K0φ)=γ0K0φ=φ, |
so that v=u−w solves (7.1) with f replaced by f−r+Pw∈¯Ht−2a(Rn+), φ replaced by 0. This is a homogeneous Dirichlet problem as in (5.25). Since v∈H(μ−1)(σ)(¯Rn+) with γ0(v/xμ−1n)=0, it is in Hμ(σ)(¯Rn+). It then follows from Theorem 5.8 that v∈Hμ(t)(¯Rn+), and hence u=v+w∈H(μ−1)(t)(¯Rn+).
For the hatted version ˆP there is even an existence and uniqueness result in these spaces.
Theorem 7.2. Let P satisfy Assumption 3.1 with Reμ>0, and let −12<t−Reμ<12. Then r+ˆP together with γμ−10 defines a homeomorphism:
{r+ˆP,γμ−10}:H(μ−1)(t)(¯Rn+)∼→¯Ht−2a(Rn+)×Ht−Reμ+12(Rn−1). | (7.5) |
Proof. The forward mapping properties are accounted for above. The existence of a unique solution u∈H(μ−1)(t)(¯Rn+) of
r+ˆPu=fonRn+,γμ−10u=φonRn−1,suppu⊂¯Rn+, | (7.6) |
for given f∈¯Ht−2a(Rn+), φ∈Ht−Reμ+12(Rn−1), is shown as in Theorem 7.1, now referring to Theorem 5.5 instead of Theorem 5.8.
These theorems show that H(μ−1)(t)(¯Rn+) is the correct domain space for the nonhomogeneous Dirichlet problem, at least in the small range −12<t−Reμ<12. Recall that Eμ−1(¯Rn+)∩E′(Rn) and e+xμ−1nS(¯Rn+) are dense subsets of H(μ−1)(t)(¯Rn+) for all t>Reμ−32.
We now show a "halfways Green's formula", where one factor u is in the domain of the nonhomogeneous Dirichlet problem for P and the other factor v is in the domain of the homogeneous Dirichlet problem for P∗:
Theorem 7.3. Let P satisfy Assumption 3.2, and assume moreover that 0<Reμ<a+12, Reμ′>0.
For u∈Eμ−1(¯Rn+)∩E′(Rn) and v∈Eˉμ′(¯Rn+)∩E′(Rn), there holds
∫Rn+Puˉvdx−∫Rn+u¯P∗vdx=−Γ(μ)Γ(μ′+1)∫Rn−1s0γ0(u/xμ−1n)γ0(ˉv/xμ′n)dx′, | (7.7) |
where s0=e−iπδp(0,1). The formula extends to u∈H(μ−1)(t)(¯Rn+) with t>Reμ−12, v∈Hˉμ′(t′)(¯Rn+) with t′>Reμ′+12.
The left-hand side is interpreted as follows, for small ε>0:
⟨r+Pu,v⟩¯H−Reμ′−12+ε(Rn+),˙HReμ′+12−ε(¯Rn+)−⟨u,P∗v⟩˙HReμ−12−ε(¯Rn+),¯H−Reμ+12+ε(Rn+). | (7.8) |
Proof. We shall show how the result can be derived from Theorem 6.4. Let u∈Eμ−1(¯Rn+)∩E′(Rn) and v∈Eˉμ′(¯Rn+)∩E′(Rn). As shown in [8,p. 494], there exist functions U and u1 in Eμ(¯Rn+)∩E′(Rn) such that u=∂nU+u1.
In terms of the Hilbert spaces: When u∈H(μ−1)(t)(¯Rn+) with |t−Reμ|<12, let z=r+Ξμ−1+u∈¯Ht−Reμ+1(Rn+), then (denoting Op(⟨ξ′⟩)=⟨D′⟩)
u=Ξ−μ+1+e+z=(⟨D′⟩+∂n)Ξ−μ+e+z=u1+∂nU,withu1=⟨D′⟩Ξ−μ+e+z∈⟨D′⟩Hμ(t+1)(¯Rn+)⊂Hμ(t)(¯Rn+),U=Ξ−μ+e+z∈Hμ(t+1)(¯Rn+),∂nU∈Hμ(t)(¯Rn+). | (7.9) |
Here Hμ(t)(¯Rn+)=˙Ht(¯Rn+) since |t−Reμ|<12. Moreover, when t=Reμ−12+ε for a small ε>0, then
r+Pu=r+Pu1+r+P∂nU=r+Pu1+∂nr+PU | (7.10) |
where both terms are in ¯Ht−2a(Rn+)=¯HReμ−12+ε−2a(Rn+)=¯H−Reμ′−12+ε(Rn+); we here use Theorem 5.8 1∘.
For v, we note that when v∈Hˉμ′(t′)(¯Rn+) with t′=Reμ′+12+ε, then
v∈Hˉμ′(t′)(¯Rn+)=Ξ−ˉμ′+e+¯Ht′−Reμ′(Rn+)⊂Ξ−ˉμ′+˙H12−ε(¯Rn+)=˙HReμ′+12−ε(¯Rn+),r+P∗v∈¯Ht′−2a(Rn+)=¯HReμ′+12+ε−2a(Rn+)=¯H−Reμ+12+ε(Rn+). |
Then the dualities in (7.8) are well-defined and serve as an interpretation of the left-hand side in (7.7).
The formula (7.7) will first be proved for u∈Eμ−1(¯Rn+)∩E′(Rn) and v∈Eˉμ′(¯Rn+)∩E′(Rn), and afterwards be extended by continuity to general u,v. We use the decomposition (7.9), that leads to elements of Eμ(¯Rn+) for t→∞. When u is supported in a ball {|ξ|≤R}, we can cut u1 and U down to have support in {|ξ|≤2R}.
Consider the contribution from u1. Here there holds
⟨r+Pu1,v⟩¯H−a,˙Ha−⟨u1,r+P∗v⟩˙Ha,¯H−a=0, | (7.11) |
when u1 and v are in ˙Ha(¯Rn+), since P is of order 2a. This gives the contribution 0 to (7.7) since t=a is allowed in the definition of u1 (recall that a>Reμ−12 by hypothesis), and t′≥a holds for the values of t′ allowed in the definition of v (where t′>Reμ′+12=2a−Reμ+12>2a−a=a). Thus u1 contributes to the boundary integral with 0.
For the contribution from ∂nU, we note that, writing U=xμnw for xn>0, w∈C∞(¯Rn+),
∂nU=∂n(xμnw)=μxμ−1nw+xμn∂nwforxn>0, |
so the weighted boundary value for xn→0+ satisfies (since xμn∂nw/xμ−1n=xn∂nw→0)
γ0(∂nU/xμ−1n)=μγ0w=μγ0(U/xμn). | (7.12) |
Moreover, by a simple integration by parts,
⟨r+P∂nU,v⟩=⟨r+∂nPU,v⟩=−⟨r+PU,∂nv⟩, |
since γ0v=0 because of Reμ′>0. Thus, by use of Theorem 6.4 and (7.12),
⟨r+P∂nU,v⟩−⟨∂nU,r+P∗v⟩=−⟨r+PU,∂nv⟩−⟨∂nU,r+P∗v⟩=−Γ(μ+1)Γ(μ′+1)s0∫Rn−1γ0(U/xμn)γ0(ˉv/xμ′n)dx′=−Γ(μ)Γ(μ′+1)s0∫Rn−1γ0(∂nU/xμ−1n)γ0(ˉv/xμ′n)dx′. |
Since u1∈Eμ(¯Rn+), γ0(u1/xμ−1n)=0, so u1 can be added to ∂nU in the last integral. Adding also (7.11) to the left-hand side, we find (7.7).
Since the expressions depend continuously on u,v in the presented norms, the formula extends to the indicated spaces.
Example 7.4. Theorems 7.1 and 7.2 apply in particular to the operator L considered in (3.5)–(3.6) and Examples 5.9 and 6.5, when μ>0 (this holds automatically if a≥12, since |δ|<12). Theorem 7.3 applies to L when μ and μ′>0 (again automatically satisfied when a≥12). $
Remark 7.5. The transmission spaces can also be defined in terms of other scales of function spaces. The case of Bessel-potential spaces Hsp, 1<p<∞, is a main subject in our preceding papers. There is also the Hölder-Zygmund scale Cs∗(Rn), coinciding with the Hölder scale Cs(Rn) when s∈R+∖N, with spaces over Rn+ defined as in (5.1). Here since Cs+ε∗(Rn)⊂Hs(Rn), also Cμ(s+ε)∗(¯Rn+)⊂Hμ(s)(¯Rn+) for ε>0. (More details on such spaces in our earlier papers, e.g., in [12].) So the results dealing with forward mapping properties of r+P have useful consequences involving these spaces as well. Namely, Theorem 5.8 1∘ implies that r+P maps
r+P:Cμ(t+ε)∗(¯Rn+)→¯Ht−2a(Rn+) for Reμ−12<t<Reμ+32, |
and the integration by parts formulas in Sections 6 and 7 hold for functions in Cμ(t)∗-type spaces, for the same t.
In the opposite direction, an inclusion of an Hs-space in a Hölder spaces loses n/2 in the regularity parameter, hence does not give very good results. For better regularity results, it would be interesting to extend the above theory to Hsp-spaces with general 1<p<∞, possibly under further hypotheses; this remains to be done. More smoothness than C1 is needed for a symbol q(ξ) to be a Fourier multiplier in Lp (some well-known conditions are recalled in [15,Sect. 1.3]). There is an extension of Vishik and Eskin's work to Lp-based spaces by Shargorodsky [21], which should be useful. It is there pointed out that [6,Lemma 17.1] shows how smoothness properties carry over to the factors in the Wiener-Hopf factorization. $
The author declares no conflict of interest.
[1] | J. K. Hale, S. M. V. Lunel, Introduction to Functional Differential Equations, Springer-Verlag, New York, 1993. https://doi.org/10.1007/978-1-4612-4342-7 |
[2] | I. Podlubny, Fractional differential equations : an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications, Academic Press, New York, 1999. |
[3] |
M. Caputo, Mean fractional-order-derivatives differential equations and filters, Annali dell'Università di Ferrara, 41 (1995), 73–84. https://doi.org/10.1007/BF02826009 doi: 10.1007/BF02826009
![]() |
[4] |
M. Caputo, M. Fabrizio, Progress in fractional differentiation and applications a new definition of fractional derivative without singular kernel, Prog. Fractional Differ. Appl., 1 (2015), 73–85. https://doi.org/10.12785/pfda/010201 doi: 10.12785/pfda/010201
![]() |
[5] |
A. Atangana, D. Baleanu, New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model, Ther. Sci., 20 (2016), 763–769. https://doi.org/10.2298/TSCI160111018A doi: 10.2298/TSCI160111018A
![]() |
[6] |
R. Khalil, M. Al Horani, A. Yousef, M. Sababheh, A new definition of fractional derivative, J. Comput. Appl. Math., 264 (2014), 65–70. https://doi.org/10.1016/j.cam.2014.01.002 doi: 10.1016/j.cam.2014.01.002
![]() |
[7] |
J. Sousa, E. V. De Oliviera, A new truncated m-fractional derivative type unifying some fractional derivative types with classical properties, arXiv: Classical Anal. ODEs, 16 (2017), 83–96. https://doi.org/10.28924/2291-8639-16-2018-83 doi: 10.28924/2291-8639-16-2018-83
![]() |
[8] |
A. Atangana, D. Baleanu, A. Alsaedi, New properties of conformable derivative, Open Math., 13 (2015), 889–898. https://doi.org/10.1515/math-2015-0081 doi: 10.1515/math-2015-0081
![]() |
[9] |
R. Khalil, M. A. Horani, M. A. Hammad, Geometric meaning of conformable derivative via fractional cords, J. Math. Comput. Sci, 19 (2019), 241–245. https://doi.org/10.22436/jmcs.019.04.03 doi: 10.22436/jmcs.019.04.03
![]() |
[10] |
E. Bas, B. Acar, R. Ozarslan, The price adjustment equation with different types of conformable derivatives in market equilibrium, AIMS Math., 4 (2019), 805–820. https://doi.org/10.3934/math.2019.3.805 doi: 10.3934/math.2019.3.805
![]() |
[11] |
B. Xin, W. Peng, Y. Kwon, Y. Liu, Modeling, discretization, and hyperchaos detection of conformable derivative approach to a financial system with market confidence and ethics risk, Adv. Differ. Equations, 2019 (2019), 138. https://doi.org/10.1186/s13662-019-2074-8 doi: 10.1186/s13662-019-2074-8
![]() |
[12] |
M. Yavuz, B. Yaskiran, Conformable derivative operator in modelling neuronal dynamics, Appl. Appl. Math., 13 (2018), 803–817. https://doi.org/10.1186/s13662-019-2074-8 doi: 10.1186/s13662-019-2074-8
![]() |
[13] |
S. Uçar, N. Y. Özgür, B. B. İ. Eroğlu, Complex conformable derivative, Arabian J. Geosci., 12 (2019), 1–6. https://doi.org/10.1007/s12517-019-4396-y doi: 10.1007/s12517-019-4396-y
![]() |
[14] |
B. Ghanbari, D. Baleanu, New optical solutions of the fractional gerdjikov-ivanov equation with conformable derivative, Front. Phys., 8 (2020). https://doi.org/10.3389/fphy.2020.00167 doi: 10.3389/fphy.2020.00167
![]() |
[15] |
M. Senol, O. Tasbozan, A. Kurt, Numerical solutions of fractional Burgers' type equations with conformable derivative, Chinese J. Phys., 58 (2019), 75–84. https://doi.org/10.1016/j.cjph.2019.01.001 doi: 10.1016/j.cjph.2019.01.001
![]() |
[16] |
B. Ghanbari, M. S. Osman, D. Baleanu, Numerical solutions of fractional Burgers' type equations with conformable derivative, Mod. Phys. Lett. A, 34 (2019). https://doi.org/10.1142/S0217732319501554 doi: 10.1142/S0217732319501554
![]() |
[17] |
A. A. Hyder, A. H. Soliman, Exact solutions of space-time local fractal nonlinear evolution equations: A generalized conformable derivative approach, Results Phys., 17 (2020), 103135. https://doi.org/10.1016/j.rinp.2020.103135 doi: 10.1016/j.rinp.2020.103135
![]() |
[18] |
A. Yusuf, M. Inc, A. A. Isa, Fractional solitons for the nonlinear Pochhammer-Chree equation with conformable derivative, J. Coupled Syst. Multiscale Dyn., 6 (2018), 158–162. https://doi.org/10.1166/jcsmd.2018.1149 doi: 10.1166/jcsmd.2018.1149
![]() |
[19] |
M. Inc, A. Yusuf, A. A. Isa, D. Baleanu, Dark and singular optical solitons for the conformable space-time nonlinear Schrödinger equation with Kerr and power law nonlinearity, Optik, 162 (2018), 65–75. https://doi.org/10.1016/j.ijleo.2018.02.085 doi: 10.1016/j.ijleo.2018.02.085
![]() |
[20] |
G. Yel, T. A. Sulaiman, H. M. Baskonus, On the complex solutions to the (3 + 1)-dimensional conformable fractional modified KdV-Zakharov-Kuznetsov equation, Mod. Phys. Lett. B, 34 (2020), 2050069. https://doi.org/10.1142/S0217984920500694 doi: 10.1142/S0217984920500694
![]() |
[21] |
H. Rezazadeh, D. Kumar, T. A. Sulaiman, H. Bulut, New complex hyperbolic and trigonometric solutions for the generalized conformable fractional Gardner equation, Mod. Phys. Lett. B, 33 (2019), 1950196. https://doi.org/10.1142/S0217984919501963 doi: 10.1142/S0217984919501963
![]() |
[22] |
H. Bulut, T. A. Sulaiman, H. M. Baskonus, H. Rezazadeh, M. Eslami, M. Mirzazadeh, Optical solitons and other solutions to the conformable space–time fractional Fokas–Lenells equation, Optik, 172 (2018), 20–27. https://doi.org/10.1016/j.ijleo.2018.06.108 doi: 10.1016/j.ijleo.2018.06.108
![]() |
[23] |
Z. Korpinar, M. Inc, A. S. Alshomrani, D. Baleanu, The deterministic and stochastic solutions of the schrodinger equation with time conformable derivative in birefrigent fibers, AIMS Math., 5 (2020), 2326–2345. https://doi.org/10.3934/math.2020154 doi: 10.3934/math.2020154
![]() |
[24] |
H. Yépez-Martinez, J. F. Gómez-Aguilar, Local M-derivative of order α and the modified expansion function method applied to the longitudinal wave equation in a magneto electro-elastic circular rod, Opt. Quantum Electron., 50 (2018), 375. https://doi.org/10.1007/s11082-018-1643-5 doi: 10.1007/s11082-018-1643-5
![]() |
[25] |
H. M. Baskonus, J. F. Gómez-Aguilar, New singular soliton solutions to the longitudinal wave equation in a magneto-electro-elastic circular rod with M-derivative, Mod. Phys. Lett. B, 33 (2019), 1950251. https://doi.org/10.1142/S0217984919502518 doi: 10.1142/S0217984919502518
![]() |
[26] |
H. M. Baskonus, New singular soliton solutions to the longitudinal wave equation in a magneto-electro-elastic circular rod with M-derivative, Eur. Phys. J. Plus, 134 (2019), 322. https://doi.org/10.1140/epjp/i2019-12680-4 doi: 10.1140/epjp/i2019-12680-4
![]() |
[27] |
H. Yépez-Martinez, J. F. Gómez-Aguilar, M-derivative applied to the soliton solutions for the Lakshmanan-Porsezian-Daniel equation with dual-dispersion for optical fibers, Opt. Quantum Electron., 51 (2019), 31. https://doi.org/10.1007/s11082-018-1740-5 doi: 10.1007/s11082-018-1740-5
![]() |
[28] |
M. S. Osman, A. Zafar, K. K. Ali, W. Razzaq, Novel optical solitons to the perturbed Gerdjikov-Ivanov equation with truncated M-fractional conformable derivative, Optik, 222 (2020), 165418. https://doi.org/10.1016/j.ijleo.2020.165418 doi: 10.1016/j.ijleo.2020.165418
![]() |
[29] |
K. U. Tariq, M. Younis, S. T. R. Rizvi, H. Bulut, M-Truncated fractional optical solitons and other periodic wave structures with Schrödinger-Hirota equation, Mod. Phys. Lett. B, 34 (2020), 2050427. https://doi.org/10.1142/S0217984920504278 doi: 10.1142/S0217984920504278
![]() |
[30] |
E. Bas, B. Acay, The direct spectral problem via local derivative including truncated Mittag-Leffler function, Appl. Math. Comput., 367 (2020), 124787. https://doi.org/10.1016/j.amc.2019.124787 doi: 10.1016/j.amc.2019.124787
![]() |
[31] |
T. A. Sulaiman, G. Yel, H. Bulut, M-fractional solitons and periodic wave solutions to the Hirota-Maccari system, Mod. Phys. Lett. B, 33 (2019), 1950052. https://doi.org/10.1142/S0217984919500520 doi: 10.1142/S0217984919500520
![]() |
[32] |
A. Zafar, A. Bekir, M. Raheel, W. Razzaq, Optical soliton solutions to Biswas-Arshed model with truncated M-fractional derivative, Optik, 222 (2020), 165355. https://doi.org/10.1016/j.ijleo.2020.165355 doi: 10.1016/j.ijleo.2020.165355
![]() |
[33] |
A. Yusuf, M. Inc, D. Baleanu, Optical solitons with M-truncated and beta derivatives in nonlinear optics, Front. Phys., 7 (2019), 126. https://doi.org/10.3389/fphy.2019.00126 doi: 10.3389/fphy.2019.00126
![]() |
[34] | B. Guo, Z. Gan, L. Kong, J. Zhang, The Zakharov System and its Soliton Solutions, Springer, Singapore, 2016. https: //doi.org/10.1007/978-981-10-2582-2 |
[35] |
N. Mahak, G. Akram, Extension of rational sine-cosine and rational sinh-cosh techniques to extract solutions for the perturbed NLSE with Kerr law nonlinearity, Eur. Phys.l J. Plus, 134 (2019), 159. https://doi.org/10.1140/epjp/i2019-12545-x doi: 10.1140/epjp/i2019-12545-x
![]() |
[36] |
N. Mahak, G. Akram, Exact solitary wave solutions by extended rational sine-cosine and extended rational sinh-cosh techniques, Phys. Scr., 94 (2019), 115212. https://doi.org/10.1088/1402-4896/ab20f3 doi: 10.1088/1402-4896/ab20f3
![]() |
[37] |
M. Cinar, I. Onder, A. Secer, A. Yusuf, T. A. Sulaiman, M. Bayram, H. Aydin, The analytical solutions of Zoomeron equation via extended rational sin-cos and sinh-cosh methods, Phys. Scr., 96 (2021), 094002. https://doi.org/10.1088/1402-4896/ac0374 doi: 10.1088/1402-4896/ac0374
![]() |
[38] |
B. Malomed, D. Anderson, M. Lisak, M. L. Quiroga-Teixeiro, L. Stenflo, Dynamics of solitary waves in the Zakharov model equations, Phys. Rev. E, 55 (1997), 962–968. https://doi.org/10.1103/PhysRevE.55.962 doi: 10.1103/PhysRevE.55.962
![]() |
[39] | V. Zakharov, Collapse of langmuir waves, Sov. Phys. JETP, 35 (1972). |
[40] |
J. Ginibre, Y. Tsutsumi, G. Velo, On the Cauchy problem for the Zakharov system, J. Funct. Anal., 151 (1997), 384–436. https://doi.org/10.1006/jfan.1997.3148 doi: 10.1006/jfan.1997.3148
![]() |
[41] |
M. Marklund, On the Cauchy problem for the Zakharov system, Phys. Plasmas, 12 (2005), 082110. https://doi.org/10.1063/1.2012147 doi: 10.1063/1.2012147
![]() |
[42] |
J. Bourgain, J. Colliander, On wellposedness of the Zakharov system, Int. Math. Res. Not., 1996 (1995), 515–546. https://doi.org/10.1155/S1073792896000359 doi: 10.1155/S1073792896000359
![]() |
[43] |
A. P. Misra, D. Ghosh, A. R. Chowdhury, A novel hyperchaos in the quantum Zakharov system for plasmas, Phys. Lett., Sect. A: Gen., At. Solid State Phys., 372 (2008), 1469–1476. https://doi.org/10.1016/j.physleta.2007.09.054 doi: 10.1016/j.physleta.2007.09.054
![]() |
[44] |
W. Bao, F. Sun, G. W. Wei, Numerical methods for the generalized Zakharov system, J. Comput. Phys., 190 (2003), 201–228. https://doi.org/10.1016/S0021-9991(03)00271-7 doi: 10.1016/S0021-9991(03)00271-7
![]() |
[45] |
S. T. Demiray, H. Bulut, Some exact solutions of generalized Zakharov system, Wave Rand. Complex Media, 25 (2015), 75–90. https://doi.org/10.1080/17455030.2014.966798 doi: 10.1080/17455030.2014.966798
![]() |
[46] |
A. Borhanifar, M. M. Kabir, L. M. Vahdat, New periodic and soliton wave solutions for the generalized Zakharov system and (2 + 1)-dimensional Nizhnik-Novikov-Veselov system, Chaos, Solitons Fractals, 42 (2009), 1646–1654. https://doi.org/10.1016/j.chaos.2009.03.064 doi: 10.1016/j.chaos.2009.03.064
![]() |
[47] |
Y. Xin, Y. Xu, C. W. Shu, Local discontinuous Galerkin methods for the generalized Zakharov system, J. Comput. Phys., 229 (2010), 1238–1259. https://doi.org/10.1016/j.jcp.2009.10.029 doi: 10.1016/j.jcp.2009.10.029
![]() |
[48] |
Q. Shi, Q. Xiao, X. Liu, Extended wave solutions for a nonlinear Klein-Gordon-Zakharov system, Appl. Math. Comput., 218 (2012), 9922–9929. https://doi.org/10.1016/j.amc.2012.03.079 doi: 10.1016/j.amc.2012.03.079
![]() |
[49] |
D. Lu, A. R. Seadawy, M. M. A. Khater, Structure of solitary wave solutions of the nonlinear complex fractional generalized Zakharov dynamical system, Adv. Differ. Equations, 2018 (2018), 266. https://doi.org/10.1186/s13662-018-1734-4 doi: 10.1186/s13662-018-1734-4
![]() |
[50] | E. W. Weisstein, CRC Concise Encyclopedia of Mathematics, Chapman and Hall/CRC, New York, 2002. https://doi.org/10.1201/9781420035223 |
1. | Gerd Grubb, Resolvents for fractional-order operators with nonhomogeneous local boundary conditions, 2023, 284, 00221236, 109815, 10.1016/j.jfa.2022.109815 | |
2. | B.-Wolfgang Schulze, Jörg Seiler, Truncation quantization in the edge calculus, 2023, 14, 1662-9981, 10.1007/s11868-023-00564-0 |