This paper develops a unified approach combining Mathematical Scattering Theory and the Factorization Method to solve inverse scattering problems involving the Laplace operator with boundary conditions on Lipschitz surfaces, determining boundaries from scattering data.
Contribution
It introduces a general scheme for inverse scattering with singular perturbations of the Laplacian, applicable to boundary conditions on Lipschitz surfaces, using a combined theoretical framework.
Findings
01
Boundary surfaces are uniquely determined by single-frequency scattering data.
02
The method applies to standard boundary conditions on Lipschitz domains.
03
Results extend to boundary conditions on subsets of surfaces.
Abstract
We provide a general scheme, in the combined frameworks of Mathematical Scattering Theory and Factorization Method, for inverse scattering for the couple of self-adjoint operators (Δ,Δ), where Δ is the free Laplacian in L2(R3) and Δ is one of its singular perturbations, i.e., such that the set {u∈H2(R3)∩dom(Δ):Δu=Δu} is dense. Typically Δ corresponds to a self-adjoint realization of the Laplace operator with some kind of boundary conditions imposed on a null subset; in particular our results apply to standard, either separating or semi-transparent, boundary conditions at Γ=∂Ω, where Ω⊂R3 is a bounded Lipschitz domain. Similar results hold in the case the boundary conditions are assigned…
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Full text
Inverse Scattering for the Laplace operator with boundary conditions on Lipschitz surfaces
Andrea Mantile
and
Andrea Posilicano
Laboratoire de Mathématiques, Université de Reims -
FR3399 CNRS, Moulin de la Housse BP 1039, 51687 Reims, France
DiSAT, Sezione di Matematica, Università dell’Insubria, via Valleggio 11, I-22100
Como, Italy
We provide a general scheme, in the combined frameworks of Mathematical Scattering Theory and Factorization Method, for inverse scattering for the couple of self-adjoint operators (Δ,Δ), where Δ is the free Laplacian in L2(R3) and Δ is one of its singular perturbations, i.e., such that the set {u∈H2(R3)∩dom(Δ):Δu=Δu} is dense. Typically Δ corresponds to a self-adjoint realization of the Laplace operator with some kind of boundary conditions imposed on a null subset; in particular our results apply to standard, either separating or semi-transparent, boundary conditions at Γ=∂Ω, where Ω⊂R3 is a bounded Lipschitz domain. Similar results hold in the case the boundary conditions are assigned only on Σ⊂Γ, a relatively open subset with a Lipschitz boundary. We show that either the obstacle Ω or the screen Σ are determined by the knowledge of the Scattering Matrix, equivalently of the Far Field Operator, at a single frequency.
1. Introduction
In the recent paper [22] (also see [24] for the case of smooth boundaries and [5] for similar results in the case of smooth boundaries and under additional trace-class conditions) we obtained a representation formula for the scattering matrix SλΛ:L2(S2)→L2(S2) relative to the scattering couple (Δ,ΔΛ), where Δ is the self-adjoint free Laplacian in L2(R3) and ΔΛ is a self-adjoint realization of the Laplacian with boundary conditions at Γ, the Lipschitz boundary of the bounded domain Ω⊂R3. Here Λ:z↦Λz is an operator-valued map which univocally defines ΔΛ and fixes the boundary conditions realized by the corresponding operator (see Sections 4.1 and 5.1 below for various explicit examples). Our representation formula gives SλΛ=1−2πiLλΛλ+Lλ∗, where Λλ+ is the limit of Λλ+iϵ as ϵ↓0 (which, under suitable hypotheses, exists in operator norm through a Limiting Absorption Principle, see [22]), and Lλ is defined in term of the trace (either Dirichlet or Neumann or both) at the boundary Γ of the free waves with wavenumber ∣λ∣1/2. Introducing the Far-Field operator FλΛ:=2πi1(1−SλΛ) (see [18, relation (1.31)]), one gets FλΛ=LλΛλ+Lλ∗; such a factorized form suggests to study the inverse scattering problem (concerning the reconstruction of the shape of Ω by the knowledge of the scattering data at a fixed frequency) by means of Kirsch’s Factorization Method (see [18] and references therein).
Our result is the following (see Theorem 4.14 for the complete statement): let Λλ+=(Mλ+)−1, where the bounded bijective operator Mλ+ has the decomposition Mλ+=M∘+Kλ, M∘ sign-definite and Kλ compact; then
[TABLE]
where ϕλx(ξ):=ei∣λ∣1/2ξ⋅x and the sequences {zλ,kΛ}1∞⊂C\{0} and {ψλ,kΛ}1∞⊂L2(S2) provide the spectral resolution of the compact normal operator FλΛ. While such a result conforms to the standard ones (the inf-criterion and the (F∗F)1/4-method) provided in [18, Section 1.4], its advantage is due to the fact that
we use a factorization where all the informations regarding the boundary conditions are encoded in the operator Λλ+, whereas Lλ, for which one needs to characterize the range, is model-independent; this enhances the flexibility of our approach. Moreover, with a minimal effort (which in essence consists in compressing the operator Λz onto subspaces of functions with supports contained in Σ⊂Γ) one gets similar results in the case the boundary conditions are imposed not on the whole Γ but only on a relatively open subset Σ with a Lipschitz boundary. In this case the result is of the same kind, only the family of testing functions changes (see Theorem 5.6 for the precise statement): let Σ∘⊂Γ∘, Γ∘ a Lipschitz boundary; then
[TABLE]
where
ϕλΣ∘(ξ):=∫Σ∘ϕλx(ξ)dσΓ∘(x).
We provide several examples where Theorems 4.14 and 5.6 apply. In particular, we consider obstacles and screens reconstruction for the following boundary conditions:
∙
Dirichlet γ0u=0 (see Subsections 4.1.1 and 5.1.1);
∙
Neumann γ1u=0 (see Subsections 4.1.2 and 5.1.2);
∙
semi-transparent
[TABLE]
either α>0 or α<0 (see Subsections 4.1.3 and 5.1.3);
b11<0, b22 real-valued (see Subsections 4.1.5 and 5.1.5).
A huge literature is devoted to obstacle reconstruction from scattering data; we just recall some papers where the Factorization Method is used in connection with the models here treated. Dirichlet and Neumann obstacles have been considered in [17] (see also [18, Chap. 1]); Dirichlet screens have been studied firstly, in a 2-dimensional setting, in [20]. Semi-transparent interface conditions appear, apart in quantum mechanical models (see, e.g., [8], [4] and references therein), in connections with acoustic models with gradient singularities, see [25]. Conditions of the type αγ0=[γ1]u appear in [19] and [6] in a non self-adjoint setting (i.e. when α is complex-valued): this compels the use of different data operators.
An appropriate choice of the functions bij in (1.1) gives the classical Robin boundary conditions; the latter have been considered in [11] (see also [18, Chap. 2]) and [7]. In these papers, as in the previous case, a non self-adjoint setting is used and different data operators enters in the reconstruction formulae.
In this paper, as regards scattering, we use a quantum mechanics point of view (see Section 3.2); however, as recalled in Section 3.3 below (see also [35] for the case of Neumann boundary conditions), the scattering theory for Schrödinger-type equations is equivalent to the one for wave-type equations. Hence our reconstruction results apply to diffusions of both classical and quantum waves.
In order to simplify the exposition, our results are stated in dimension d=3; however they hold in any dimension d≥2. Finally, we presume that, by the same techniques, our approach can be extended to the case in which the Laplace operator is replaced by a more general 2nd order elliptic differential operator.
Acknowledgements. The authors are indebted to Mourad Sini for the
fruitful discussions which largely inspired this work.
2. Notations and preliminaries.
2.1. Notations.
∙∥⋅∥X denotes the norm on the complex Banach space X; in case X is a Hilbert space, ⟨⋅,⋅⟩X denotes the (conjugate-linear w.r.t. the first argument) scalar product.
∙⟨⋅,⋅⟩X∗,X denotes the duality (assumed to be conjugate-linear w.r.t. the first argument) between the dual couple (X∗,X).
∙L∗:dom(L∗)⊆Y∗→X∗ denotes the dual of the densely defined linear operator L:dom(L)⊆X→Y; in a Hilbert spaces setting L∗ denotes the adjoint operator.
∙ρ(A) and σ(A) denote the resolvent set and the spectrum of the self-adjoint operator A; σp(A), σpp(A), σac(A), σsc(A), σess(A), σdisc(A), denote the point, pure point, absolutely continuous, singular continuous, essential and discrete spectra.
∙B(X,Y), B(X)≡B(X,X), denote the Banach space of bounded linear operator on the Banach space X to the Banach space Y; ∥⋅∥X,Y denotes the corresponding norm.
∙X↪Y means that X⊆Y and for any u∈X there exists c>0 such that
∥u∥Y≤c∥u∥X; we say that X is continuously embedded into Y.
∙u∣Γ denotes the restriction of the function u to the set Γ; L∣V denotes the restriction of the linear operator L to the subspace V.
∙Hs(R3), s∈R, denotes the scale of Hilbert space of Sobolev functions on R3, i.e. u∈Hs(R3) if and only if k↦(1+∥k∥2)s/2u(k) is square integrable, u denoting Fourier transform.
∙Ω≡Ωin⊂R3 denotes a bounded open set with a Lipschitz boundary Γ; Ωex:=R3\Ω.
∙γ0 and γ1 denote the Dirichlet and Neumann traces on the boundary Γ.
∙ΔΩin/exD denotes the self-adjoint operator in L2(Ωin/ex) representing the Laplace operator with homogeneous Dirichlet boundary conditions at Γ.
∙ΔΩin/exN denotes the self-adjoint operator in L2(Ωin/ex) representing the Laplace operator with homogeneous Neumann boundary conditions at Γ.
∙Hs(Ωin/ex), s∈R, denotes the scale of Hilbert space of Sobolev functions on Ωin/ex.
∙Cκ(Γ) denotes the space of Hölder-continuous functions of order κ on Γ.
∙Hs(Γ), ∣s∣≤1, denotes the Hilbert space of Sobolev functions of order s on Γ.
∙M(Hs(Γ),Ht(Γ)), M(Hs(Γ),Hs(Γ))≡M(Hs(Γ)), denotes the space of Sobolev multipliers from Hs(Γ) to Ht(Γ).
∙s♯, ♯=D,N, denote the indices
sD=1/2, sN=−1/2.
∙φn⇀φ means that the sequence {φn}1∞ weakly converges to φ.
∙V⊥⊆X∗, denotes the annihilator {\mathsf{V}}^{\perp}=\{x^{*}\in\mathsf{X}^{*}:\langle x^{*},x\rangle_{\mathsf{X}^{*}\!,\mathsf{X}}=0\ \text{for all x\in{\mathsf{V}}}\} of the subspace V⊆X.
2.2. Trace maps and layer operators on Lipschitz manifolds.
Let Γ be the compact Lipschitz manifold given by the boundary of Ω⊂R3. Let γ0 be the map defined by the restriction of u∈Ccomp∞(R3) along the set Γ: γ0u:=u∣Γ. Then, by [13, Theorem 1, Chapter VII], such a map has a bounded and surjective extension to Hs+1/2(R3) for any s>0:
[TABLE]
Here the Hilbert space B2,2s(Γ) is a Besov-like space (see [13, Section 2, Chapter V] for the precise definitions); B2,2s(Γ) identifies with Hs(Γ) whenever
0<s<1 (see [13, Section 1.1, chap. V]), where Hs(Γ) denotes the usual fractional Sobolev space on Γ (see e.g. [26, Chapter 3]). If Γ is a manifold of class Cκ,1, κ≥0, then B2,2s(Γ)=Hs(Γ) for any s≤κ+1. We use the following notations for the dual (with respect to the L2(Γ)-pairing) spaces: (B2,2s(Γ))∗≡B2,2−s(Γ).
By [33, Proposition 20.5], the embeddings
B2,2s2(Γ)↪B2,2s1(Γ), s2>s1, and B2,2s(Γ)↪L2/(1−s)(Γ), 0<s<1, are compact.
Let Δ:Hs+2(R3)→Hs(R3) be the distributional Laplacian; in the following the resolvent Rz0≡(−Δ+z)−1, z∈C\(−∞,0], is viewed as a map in B(Hs(Rn),Hs+2(Rn)), s∈R. Given s>0, by the mapping properties (2.1) one gets, for the dual of the trace map,
[TABLE]
and so we can define the bounded operator (the single-layer potential)
[TABLE]
By resolvent identity one has
[TABLE]
By (2.1) and (2.2), one obtains the bounded operator
[TABLE]
In the following ΔΩin/ex denote the distributional Laplacians on Ωin/ex.
The one-sided, zero and first order, trace operators γ0in/ex and γ1in/ex=ν⋅γ0in/ex∇ (ν denoting the outward normal vector at the boundary) defined on
smooth functions in Ccomp∞(Ωin/ex) extend to bounded and surjective linear operators (see e.g. [26, Theorem
3.38])
[TABLE]
and
[TABLE]
(we refer to [26, Chapter 3] for the definition of the Sobolev spaces Hs(Ωin/ex) and Hs(Γ)). Using these maps and setting Hs(R3\Γ):=Hs(Ωin)⊕Hs(Ωex), the two-sided bounded and surjective trace operators are defined according
to
[TABLE]
[TABLE]
while the corresponding jumps are
[TABLE]
[TABLE]
Let us notice that in the case u=uin⊕uex∈Hs+1/2(Rn), 0<s<1, γ0 in (2.6) coincides with the map defined in (2.1) and so there is no ambiguity in our notations; this also entails that γ0 remains surjective even if restricted to H2(R3). Similarly the map γ1 is surjective onto Hs(Γ) even if restricted to Hs+3/2(R3).
By [26, Lemma 4.3], the trace maps γ1in/ex can be extended to the spaces
[TABLE]
[TABLE]
This gives the analogous extensions of the maps γ1 and [γ1] defined on HΔ1(R3\Γ):=HΔ1(Ωin)⊕HΔ1(Ωex) with values in H−1/2(Γ).
By using a cut-off function χ∈Ccomp∞(Rn) such that χ=1 in a neighborhood of Ωin, all the maps defined above can be extended (and we use the same notation) to functions u such that χu is in the right function space.
The single-layer operator SLz has been already introduced above; now we recall the definition of double-layer operator DLz, z∈C\(−∞,0]: by the dual map
[TABLE]
and by the resolvent Rz0∈B(Hs(R3),Hs+2(R3)), one defines the bounded operator
[TABLE]
By resolvent identity one has
[TABLE]
By the mapping properties of the layer operators, one gets (see [26, Theorem 6.11])
[TABLE]
for any χ∈Ccomp∞(R3); by (−(ΔΩin⊕ΔΩex)+z)SLzϕ=(−(ΔΩin⊕ΔΩex)+z)DLzφ=0, one gets χSLzϕ∈HΔ1(Rn\Γ), ϕ∈H−1/2(Γ), and χDLzφ∈HΔ1(Rn\Γ), φ∈H1/2(Γ). Thus
[TABLE]
These mapping properties can be extended to a larger range of Sobolev spaces (see, e.g., [26, Theorem 6.12 and successive remarks]):
[TABLE]
By the Limiting Absorption Principle for the free Laplacian (see, e.g., [21, Section 18]), duality and interpolation, one has that the limits
[TABLE]
exist in ∈B(Hw−s(R3),H−w−s+2(R3)), w>1/2, 0≤s≤2 (here Hws(R3) denotes the weighted Sobolev space of order s with weight φ(x)=(1+∥x∥2)w/2). Thus, since Γ is bounded,
the limits
[TABLE]
exist in B(B2,2−s(Γ),H−w3/2−s(R3)), 0<s≤3/2, and
B(H−s(Γ),H−w1/2−s(R3)), 0<s≤1/2, respectively. Moreover, by the identities (2.3),(2.11) and by SLz∈B(B2,2−3/2(Γ),Lw2(Rn)), DLz∈B(H−1/2(Γ),Lw2(Rn)) (see [24, relation (4.10)]) one has
[TABLE]
3. Direct Scattering Theory for Singular Perturbations.
3.1. Singular Perturbations of the Laplace operator.
Let Δ:H2(R3)⊆L2(R3)→L2(R3) be the self-adjoint operator given by the free Laplacian on the whole space. Another self-adjoint operator Δ:dom(Δ)⊆L2(R3)→L2(R3) is said to
be a singular perturbation of Δ if the set
[TABLE]
is dense in
L2(R3). Our aim is the study of direct and inverse scattering for the couple (Δ,Δ). Notice that Δ is a self-adjoint extension of the symmetric operator Δ∘:=Δ∣D≡Δ∣D; in typical situations Δ represents the Laplace operator with some kind of boundary condition holding on a null subset.
3.2. Wave Operators.
Given the two self-adjoint operators Δ and Δ, let eitΔ and eitΔ be the corresponding unitary groups of evolution providing solutions of the Cauchy problems for the Schrödinger equations
[TABLE]
As usual in Quantum Mechanics (see, e.g., [31]), we define the Wave Operators for the scattering couple (Δ,Δ) as
[TABLE]
One says that W±(Δ,Δ) exist whenever the limits exist for any vector u∈L2(R3) and then that are complete whenever
[TABLE]
where L2(R3)ac denotes the absolutely continuous subspace of Δ.
It is known that the existence of both the wave operators W±(Δ,Δ) and W±(Δ,Δ) gives completeness. From the point of view of physical interpretation, a more relevant definition is the following: W±(Δ,Δ) are said to be asymptotically complete whenever they are complete and
[TABLE]
where L2(R3)pp denotes the pure point subspace of Δ; equivalently, whenever
they are complete and the singular continuous spectrum of Δ is empty: σsc(Δ)=∅. In this case L2(R3) decomposes into the direct sum of scattering states and bound states.
3.3. Scattering theory for wave equations.
Suppose that Δ is real (i.e., it maps real-valued functions to real-valued functions), not positive and injective (these hypotheses can be weakened, it suffices to require Δ upper semi-bounded, see [15, Sections 8 and 9], [3, Section 10.3]). Let Hhom1(R3) be the homogeneous Sobolev space of order one and let Hhom1(R3) the completion, with respect to the norm ∥u∥:=∥(−Δ)1/2u∥L2(R3), of dom(−Δ)1/2).
Then the unitary group of evolutions providing the solutions of the Cauchy problems with real initial conditions
[TABLE]
are unitary equivalent, by the maps
[TABLE]
to the Schrödinger unitary groups in the complex Hilbert space L2(R3) given by
e−it(−Δ)1/2 and e−it(−Δ)1/2 respectively. By the Kato-Birman invariance principle (see, e.g., [3, Section 11.3.3]), if both the wave operators W±(Δ,Δ) and
W±(−(−Δ)1/2,−(−Δ)1/2) exist, then they are equal (by the Kato-Birman criterion, see [16, Theorem 4.8, Chapter X], equality holds whenever the difference of some power of the resolvents is trace-class; for the models discussed below this is true under some additional regularity hypotheses on Γ, see [23, Theorems 4.11 and 4.12]). In this case
the scattering theory for the couple of Schrödinger equations (3.1) is equivalent to the one for the couple of wave equations
[TABLE]
3.4. A resolvent formula for singular perturbations.
Given an auxiliary Hilbert space K, we introduce a linear application τ:H2(R3)→K
which plays the role of an abstract trace (evaluation) map. We assume that 1. τ is continuous;
τ is surjective (so that K plays the role of the trace space);
ker(τ) is dense in L2(R3).
In the following we do not identify K with its dual K∗; however we use K∗∗≡K. Tipically K↪K0↪K∗ and the K-K∗ duality ⟨⋅,⋅⟩K∗,K (conjugate-linear with respect to the first variable) is defined in terms of the scalar product of the Hilbert space K0. For any z∈ρ(A0) we define the bounded operators
[TABLE]
and
[TABLE]
Then, given a reflexive Banach space X such that K↪X, we consider, for some not empty set ZΛ⊆C\(−∞,0] which is symmetric with respect to the real axis (i.e., z∈ZΛ⇒z∗∈ZΛ), a map
[TABLE]
such that
[TABLE]
Remark 3.1**.**
Notice that whenever there exists a family of bijections Mz∈B(X∗,X), z∈ZΛ, such that Λz=Mz−1, then (3.3) is equivalent to
[TABLE]
The following result is a useful ingredient in the successive discussion about inverse scattering:
By (3.4), one has Im⟨ϕ,Mzϕ⟩X∗,X=Im(z)∥Gzϕ∥L2(R3)2. Since Gz∗=τRz∗0 is surjective onto K, Gz has closed range by the closed range theorem. Hence, see [16, Theorem 5.2, page 231], there exists c>0 such that ∥Gzϕ∥L2(R3)2≥c∥ϕ∥K∗2. Therefore, whenever Im(z)=0,
[TABLE]
and the proof is done.
∎
Now we recall the key result about singular perturbations of Δ (see [27, Theorem 2.1], [28, Corollary 3.2], [29, Corollary 3.2], [22, Theorem 2.4]):
Theorem 3.3**.**
Let τ and Λ be as above. Then the family of bounded linear maps in L2(R3)
[TABLE]
is the resolvent of a self-adjoint operator ΔΛ which is a singular perturbation of Δ. Moreover, ΔΛ is a self-adjoint extension of the closed symmetric operator Δ∣ker(τ) and all its self-adjoint extensions (and any singular perturbation of Δ as well) are of this kind.
Remark 3.4**.**
The map Λ:z↦Λz introduced in (3.2) and (3.3) encodes the boundary conditions that the functions belonging to the self-adjointness domain of the corresponding ΔΛ have to satisfy. We refer to the successive Sections 4.1 and 5.1 below for various explicit examples. Notice that the properties required in (3.3) are necessary for the operator family z↦RzΛ in (3.5) to satisfy the first resolvent identity and (RzΛ)∗=Rz∗Λ
(see [27, page 113]).
Then, building on some results by Schechter conceived for perturbations by a regular potential (see [30, Section 9.4]), one gets a completeness criterion for the scattering couple (ΔΛ,Δ) (see [22, Theorem 2.8]):
Theorem 3.5**.**
Suppose that there exists an open subset E⊆R of full measure such that for any open and bounded I, I⊂E,
[TABLE]
and
[TABLE]
Then both the wave operators W±(ΔΛ,Δ) and W±(Δ,ΔΛ) exists and are complete.
3.5. The Scattering Matrix.
According to Theorem 3.3, whenever (3.6) and (3.7) hold, the scattering operator
[TABLE]
is a well defined unitary map. Given the direct integral representation of L2(R3) with respect to the spectral measure of Δ, i.e. the unitary map (here S2 denotes the 2-dimensional unitary sphere in R3)
[TABLE]
which diagonalizes Δ,
we define the scattering matrix
[TABLE]
by the relation
[TABLE]
The scattering matrix is better studied using Limiting Absorption Principle and stationary scattering theory (see, e.g., [34]). However, for typical scattering couples (ΔΛ,Δ), the hypotheses required in [34] are not satisfied. Thus at first one considers the scattering matrix for the resolvent couple (RμΛ,Rμ0), μ∈ρ(ΔΛ)∩(0,+∞), so to exploit the factorized form of the resolvent difference RμΛ−Rμ0 provided by formula (3.5), and then uses the Birman-Kato invariance principle (see [22, Section 4]). At the end, one obtains the following (see [22, Theorem 5.1]; notice that in reference [22], due to a repeated misprint, the t→±∞
limits has to be replaced by the t→∓∞ ones)
Theorem 3.6**.**
Let ΔΛ denote the self-adjoint operator corresponding to Λ={Λz}z∈ZΛ, Λz∈B(X,X∗), K↪X. Suppose that:
[TABLE]
[TABLE]
[TABLE]
Then asymptotic completeness holds for the scattering couple (ΔΛ,Δ).
Moreover,
[TABLE]
the scattering matrix SλΛ is given by
[TABLE]
where σp−(ΔΛ):=(−∞,0)∩σp(ΔΛ) is a (possibly empty) discrete set,
[TABLE]
and
[TABLE]
Here uλξ(x)=ei∣λ∣1/2ξ⋅x denotes the plane wave with direction ξ∈S2 and wavenumber ∣λ∣1/2.
Remark 3.7**.**
Let Λz=Mz−1 as in Remark 3.1 and suppose that the limit Mλ+:=limϵ↓0Mλ+iϵ exists in B(X∗,X). Then, by Theorem 3.6,
the inverse (Mλ+)−1 exists in B(X,X∗) and Λλ+=(Mλ+)−1.
4. Inverse Scattering for the Laplace operator with boundary conditions on Lipschitz surfaces.
With reference to Theorem 3.6 and given an open, bounded set Ω≡Ωin⊂R3 with a Lipschitz boundary Γ and such that Ωex:=R3\Ω is connected, we consider models where the map τ:H2(R3)→K corresponds to one of the following three different cases:
[TABLE]
[TABLE]
[TABLE]
These settings, with suitable choice of the map Λ, allow to obtain all the self-adjoint extensions of the closed symmetric operator Δ∣Ccomp∞(R3\Γ). In particular, any self-adjoint realization of the Laplace operator with boundary conditions prescribed either on the surface Γ or on a relatively open subset Σ⊂Γ can be defined in one of the above schemes, see [23, Theorem 4.4] for the case of smooth hypersurfaces. In the present framework, Theorem 3.6 allows the boundary Γ to be Lipschitz; in the applications we give in Sections 4.1 and 5.1 hypothesis (3.10) is always satisfied since Ω is bounded; hypotheses (3.8) and (3.9) also hold, (3.8) by a direct checking and (3.9) by compact Sobolev embeddings.
The results we provide in this section apply to the cases where the boundary conditions are assigned on the whole boundary Γ. Then ΔΛ can be interpreted as a model either of an extended obstacle or of a semi-transparent interface supported on Γ, whose physical properties are encoded by Λ.
Defining the Far Field operator
[TABLE]
the inverse scattering problem consists in recovering the shape of the obstacle Ω from the knowledge of FλΛ, or, equivalently, from knowledge of the scattering matrix SλΛ.
Notation 4.1**.**
In the following we refer to the different settings 1) - 3) above by introducing the index ♯, with ♯=D,N,DN according to the possible different choices, to label the operators
[TABLE]
associated to one of the traces τD=γ0, τN=γ1, τDN=γ0⊕γ1, and the spaces X=X♯s, where
[TABLE]
Furthermore the adopt the short-hand notations s♯, ♯=D,N, to denote the indices
sD=1/2, sN=−1/2.
Remark 4.2**.**
Since X♯0↪X♯s, and hence X♯s∗↪X♯0∗, we do not put any index s in the notation for Lλ♯, since we can always suppose that Lλ♯ acts on X♯0∗ and is then restricted to the proper space according to the case.
Lemma 4.3**.**
Let λ∈(−∞,0)\σdisc(ΔΩ♯), ♯=D,N, and set
[TABLE]
Then
[TABLE]
Proof.
Given λ∈(−∞,0), let uλ,ϕ♯ be the radiating solution (i.e satisfying the Sommerfeld radiating condition) in Ωex:=R3\Ω of Helmholtz equation (−Δ+λ)uλ,ϕ=0 with either Dirichlet (whenever ♯=D) or Neumann (whenever ♯=N) boundary condition ϕ∈Hs♯(Γ). Such a solution is unique in
[TABLE]
where uB:=u∣Ωex∩B (see, e.g., [26, Theorem 9.11] for the Dirichlet case and [26, Exercise 9.5] for the Neumann case). Then (see, e.g., [18, Theorem 1.4], [26, Exercise 9.4(iv)]) there exists a unique uλ,ϕ♯,∞∈C∞(S2) such that
[TABLE]
This defines the data-to-pattern operator
[TABLE]
Introducing the Herglotz operators Hλ♯:L2(S2)→Hs♯(Γ) defined by
[TABLE]
one has
[TABLE]
Since, (see [18, proofs of Theorems 1.15 and 1.26])
[TABLE]
one gets
[TABLE]
Since, for any s∈[0,1/2],
[TABLE]
and
[TABLE]
are bijections (by [22, relations (5.32) and (5.33)] and the regularity results in [10, Theorem 3]), one has
[TABLE]
Finally, by [18, Theorems 1.12 and 1.27] (it is easy to check that the proofs, there given for s=0, hold for any s∈[0,1/2]), one has
[TABLE]
and the thesis is proven.∎
Corollary 4.4**.**
Let λ∈(−∞,0)\(σdisc(ΔΩD)∩σdisc(ΔΩN)). Then
[TABLE]
Proof.
Let λ∈(−∞,0). Since (−Δ+λ)SLλ+(x)=(−Δ+λ)DLλ+(x)=0, x∈Ωex, one gets the identities KλDγ0SLλ+=KλNγ1SLλ+ and
KλDγ0DLλ+=KλNγ1DLλ+. Thus, given ϕ⊕φ∈Hs−1/2(Γ)⊕Ht+1/2(Γ), one has
[TABLE]
Therefore the thesis is consequence of (4.6), (4.7) and Lemma 4.3.
∎
Let us recall the following definitions:
Definition 4.5**.**
Let Y be a reflexive Banach space. C∈B(Y∗,Y) is said to be:
coercive, whenever there exists c>0 such that
[TABLE]
positive, whenever C=C∗ and there exists c>0 such that
[TABLE]
sign-definite, whenever either C or −C is positive.
Remark 4.6**.**
Let C∈B(Y∗,Y) be coercive. Then C∗ is injective and so ran(C) is dense by ran(C)=ker(C∗)⊥=Y. Since (4.9) implies ∥Cφ∥Y≥c∥φ∥Y∗, ran(C) is closed by [16, Theorem 5.2, page 231]. Hence C is a continuous bijection and therefore C−1∈B(Y,Y∗) by the inverse mapping theorem.
We also recall the following useful coercivity criterion (see [18, Lemma 1.17]; since our statement is slightly different from the original one, for the reader convenience we give a sketch of the proof there provided):
Lemma 4.7**.**
Let C∈B(Y∗,Y) be such that Im⟨φ,Cφ⟩Y∗,Y=0 for any φ∈Y∗\{0}. Suppose C has the decomposition C=C∘+K, where C∘=C∘∗ is coercive and K is compact. Then C is coercive.
Proof.
Supposing that C does not satisfy (4.9), one gets a sequence {φn}1∞, ∥φn∥Y∗=1, φn⇀φ, such that ⟨φn,Cφn⟩Y∗,Y→0. Since
[TABLE]
and ∥K(φn−φ)∥Y→0, one gets
[TABLE]
i.e., Im⟨φ,Cφ⟩Y∗,Y=0, which gives φ=0.
Thus φn⇀0 and the inequality
[TABLE]
is violated for n sufficiently large.
∎
Notation 4.8**.**
[TABLE]
The factorized form of the operator FλΛ, Lemma 4.3 and Corollary 4.4 suggest to take into account Kirsch’s inf-criterion:
Theorem 4.9**.**
Let λ∈E♯−∩EΛ−, ♯=D,N,DN, and suppose that the Far Field Operator can be factorized as
[TABLE]
where C∈B(Y∗,Y), Y a reflexive Banach space, is coercive and B∈B(Y,L2(S2)) is such that
By (4.9) and by [18, Theorem 1.16],
for any ϕ∈L2(S2)\{0}, one has
[TABLE]
The proof is then concluded by (4.11), Lemma 4.3 and Corollary 4.4.
∎
The next results is a key ingredient for obtaining a different identification criterion for the shape of Ω.
Theorem 4.10**.**
Let λ∈EΛ−. Then FλΛ is a normal compact operator.
Proof.
Since the scattering matrix SλΛ is unitary,
[TABLE]
and so FλΛ is normal. By
[TABLE]
and
[TABLE]
[TABLE]
one gets (here the constant c changes from line to line)
[TABLE]
Therefore Lλ is a bounded map with values in the space Lip(S2) of Lipschitz functions and so Lλ in Theorem 3.6 is a compact operator by the compact embedding Lip(S2)↪L2(S2). In conclusion, FλΛ=LλΛλ+Lλ∗ is compact since Λλ+ is bounded.
∎
Remark 4.11**.**
As consequence of Theorem 4.10 (and since 1−2πiFλΛ is unitary), by spectral theory for compact normal operators (see, e.g., [14, Section 6]), one has
[TABLE]
and there exists an orthonormal sequence {ψλ,kΛ}1∞⊂L2(S2) such that for every ψ∈L2(S2),
[TABLE]
and
[TABLE]
Remark 4.12**.**
Notice that, by Remark 4.11, {ψλ,kΛ}1∞⊂ker(FλΛ)⊥ and so ran(FλΛ)⊆ker(FλΛ)⊥.
Theorem 4.13**.**
Let FλΛ=BCB∗, where B satisfies (4.11) and C, with Im⟨φ,Cφ⟩Y∗,Y=0 for any φ∈Y∗\{0}, has the decomposition C=C∘+K, where C∘ is sign-definite and K is compact. Then
Let P0:L2(S2)→L2(S2) be the orthogonal projection
such that ran(P0)=L⊥2(S2):=ker(FλΛ)⊥. Then, by Remark 4.12, FλΛ=P0FλΛP0; hence FλΛ=P0BCB∗P0=(P0B)C(P0B)∗, and so, by [18, Theorem 1.16], ran(B)=ran(P0B). Let FλΛ:L⊥2(S2)→L⊥2(S2) be the injective normal compact operator given by the compression of FλΛ to L⊥2(S2). By Remark 4.11, {ψλ,kΛ}1∞ is an orthonormal basis in L⊥2(S2) and FλΛ=∑k=1∞zλ,kΛψλ,kΛ⊗ψλ,kΛ. By functional calculus for normal operators, using the factorization of z∈C\{0} given by z=∣z∣1/2sgn(z)∣z∣1/2, sgn(z):=∣z∣−1z, one gets
[TABLE]
Since FλΛ=BCB∗, where B:=P0B (here P0 means the surjection P0:L2(S2)→L⊥2(S2)), by [18, Theorem 1.23], ran(∣FλΛ∣1/2)=ran(B)=ran(P0B)=ran(B).
Hence ran(∣FλΛ∣1/2)=ran(Lλ♯∣X♯s∗) and so, by Lemma 4.3 and Corollary 4.4, x∈Ω if and only if ϕλx∈ran(∣FλΛ∣1/2), equivalently if and only if ϕλx∈dom(∣FλΛ∣−1/2).
Since ∣FλΛ∣−1/2=∑k=1∞∣zλ,kΛ∣−1/2ψλ,kΛ⊗ψλ,kΛ,
ϕλx∈dom(∣FλΛ∣−1/2) if and only if the series ∑k=1∞∣zλ,kΛ∣−1∣⟨ϕλx,ψλ,kΛ⟩L2(S2)∣2 converges.
∎
In applications to concrete models, the following consequence of Theorems 4.9 and 4.13 turns out to be useful:
Theorem 4.14**.**
Let
[TABLE]
and suppose that Λλ+=(Mλ+)−1, where the bijection Mλ+∈B(X♯s∗,X♯s), s∈[0,1/2], has the decomposition Mλ+=M∘++Kλ+, with M∘+ sign-definite and Kλ+ compact. Then
[TABLE]
where the sequences {zλ,kΛ}1∞⊂C\{0} and {ψλ,kΛ}1∞⊂L2(S2) provide the spectral resolution of FλΛ as in Remark 4.11 and ϕλx is defined in (4.2).
Proof.
Let us consider the factorization F^{\Lambda}_{\lambda}=\big{(}L^{\sharp}_{\lambda}(M^{+}_{\lambda})^{-1}\big{)}(M^{+}_{\lambda})^{*}\big{(}L^{\sharp}_{\lambda}(M^{+}_{\lambda})^{-1}\big{)}^{*}. Then the thesis is consequence of Lemma 4.7, Theorems 4.9 and 4.13 once one shows that
[TABLE]
Equivalently, let us prove that Im⟨ϕ,(Mλ+)∗ϕ⟩X♯s∗,X♯s=0 implies ϕ=0 (our reasonings below are inspired by the ones given in [18, page 51]).
By the definition of FλΛ and since SλΛ is unitary, one gets
[TABLE]
Setting Bλ:=Lλ♯(Mλ+)−1, this gives the identity
[TABLE]
Let ♯=N,D; then by (4.5), ker(Bλ)=ker(Kλ♯); hence, by [18, Lemma 1.13 and Theorem 1.26(b)], one has ker(Bλ)={0} and so ran(Bλ∗) is dense. Let ♯=DN; then, by (4), ker(Bλ)=ker(LND)=ker(KλDγ0Gλ+), where Gλ+(ϕ⊕φ):=(SLλ+ϕ+DLλ+φ). Since Gλ+(ϕ⊕φ) is a radiating solution of the Helmholtz equation in Ωex, γ0Gλ+(ϕ⊕φ)=0 implies Gλ+(ϕ⊕φ)=0. Hence ker(Bλ)=ker(Gλ+). Since Gλ+iϵ converges to Gλ+ in B(B2,2−3/2(Γ)⊕H1/2(Γ),Lw2(R3)) (see (2.13)) and
there exists c>0 such that, for any ϵ>0,
[TABLE]
(see [22, proof of Lemma 3.6]), Gλ+ is injective and so ran(Bλ) is dense whenever ♯=DN as well.
Let ϕ∈X♯s∗ be such that Im⟨ϕ,(Mλ+)∗ϕ⟩X♯s∗,X♯s=0; let {ψn}1∞⊂L2(S2) be a sequence such that Bλ∗ψn→ϕ. Then, by (4.12), FλΛψn→0 and so, for any ψ∈L2(S2),
[TABLE]
Therefore (Mλ+)∗ϕ∈ran(Bλ∗)⊥={0}. Since Mλ+ is a bijection, (Mλ+)∗ is injective and so ϕ=0.∎
Remark 4.15**.**
If Mλ+ in Theorem 4.14 is merely coercive, then the “inf” criterion still holds.
4.1. Applications
4.1.1. Dirichlet obstacles.
Let ΔΩin/exD denote the self-adjoint operators in L2(Ωin/ex) corresponding to the Laplace operator with Dirichlet boundary conditions. One has ΔΩinD⊕ΔΩexD=ΔΛD, where ΛzD=−(γ0SLz)−1∈B(H1/2(Γ),H−1/2(Γ)), z∈C\(−∞,0], and Theorem 3.6 holds in this case (see [22, Section 5.2]). By first resolvent identity, γ0SLλ+:H−1/2(Γ)→H1/2(Γ) can be additively decomposed as γ0SLλ+=γ0SLμ+(λ−μ)γ0Rλ0,+SLμ, μ>0,
with γ0SLμ positive (see [23, Lemma 3.2]) and γ0Rλ0,+SLμ compact (see [22, Section 5.1.3]).
Thus Theorem 4.14 applies to FλΛD, λ∈ED−.
4.1.2. Neumann obstacles.
Let ΔΩin/exN denote the self-adjoint operators in L2(Ωin/ex) corresponding to the Laplace operator with Neumann boundary conditions. One has ΔΩinN⊕ΔΩexN=ΔΛN, where ΛzN=−(γ1DLz)−1∈B(H−1/2(Γ),H1/2(Γ)), z∈C\(−∞,0], and Theorem 3.6 holds in this case (see [22, Section 5.3]). By first resolvent identity, γ1DLλ+:H1/2(Γ)→H−1/2(Γ) can be additively decomposed as γ1DLλ+=γ1DLμ+(λ−μ)γ1Rλ0,+DLμ, μ>0,
with −γ1DLμ positive (see [23, Lemma 3.2]) and γ1Rλ0,+DLμ compact (see [22, Section 5.1.3]).
Thus Theorem 4.14 applies to FλΛN, λ∈EN−.
4.1.3. Obstacles with semitransparent boundary conditions αγ0u=[γ1]u.
Here α is a real-valued function and we use the same symbol to denote the corresponding multiplication operator.
Lemma 4.16**.**
1) If α∈L6(Γ) and α1∈L∞(Γ) then (α1+γ0SLz)−1∈B(L2(Γ)), z∈C\R. 2) If both α and α1 belong to
L∞(Γ) and sgn(α) is constant, then α1+γ0SLλ+, λ∈(−∞,0), is the sum of a sign-definite operator plus a compact one.
Proof.
Since L6(Γ)⊆M(H2/3(Γ),L2(Γ)) and γ0SLz∈B(Ht−1/2(Γ),Ht+1/2(Γ)), 0<t≤1/2 (see [22, equation (5.27)], one has that 1+αγ0SLz∈B(L2(Γ)) and it is injective since it is invertible (and hence injective) as a map in H−1/3(Γ) (use [22, Lemma 5.8]). Let us now suppose that it is not surjective from L2(Γ) onto itself, i.e. we suppose that there exists ψ∈L2(Γ) such that ψ=ϕ+αγ0SLzϕ with ϕ∈H−1/3(Γ), and ϕ∈/L2(Γ).
Hence αγ0SLzϕ∈/L2(Γ), which is not possible since SLzϕ∈H2/3(Γ) and α∈M(H2/3(Γ),L2(Γ)). In conclusion 1+αγ0SLz∈B(L2(Γ)) is a bounded bijection in L2(Γ) and so (1+αγ0SLz)−1∈B(L2(Γ)) by the inverse mapping theorem. Since α is a.e. finite, α1:L2(Γ)→L2(Γ) is a continuous bijection. Hence α1+γ0SLz=α1(1+αγ0SLz) is a continuous bijection and so (α1+γ0SLz)−1∈B(L2(Γ)) by the inverse mapping theorem.
Since γ0SLλ+ maps L2(Γ) onto H1(Γ), by the compact embedding H1(Γ)↪L2(Γ), it is compact. Since
⟨φ,∣α∣1φ⟩L2(Γ)≥∥α∥L∞(Γ)−1∥φ∥L2(Γ)2 and sgn(α) is constant, α1 is sign-definite. ∎
We consider the self-adjoint operator ΔΛα, where
[TABLE]
Λzα is well-defined, i.e., Mzα has a bounded inverse, by Lemma 4.16. By [9, Theorem 2.19], the map z↦Λzα and the resolvent formula
(3.5) extend to ZΛα:=ρ(ΔΛα)∩C\(−∞,0].
ΔΛα provides a self-adjoint realization of the (bounded form above) Laplacian on R3\Γ with the semi-transparent boundary conditions at Γ given by αγ0u=[γ1]u, [γ0]u=0; moreover Theorem 3.6 holds in this case (see [22, Corollary 5.12]). By point 2 in Lemma 4.16, Theorem 4.14 applies to FλΛα, λ∈ED− (here EΛα−=(−∞,0) by [24, Remark 3.8]).
4.1.4. Obstacles with semitransparent boundary conditions γ1u=θ[γ0]u.
Here θ is a real-valued function and we use the same symbol to denote the corresponding multiplication operator.
Lemma 4.17**.**
Let θ∈Lp(Γ), p>2. Then 1) (θ−γ1DLz)−1∈B(H−1/2(Γ),H1/2(Γ)), z∈C\R ; 2) θ−γ1DLλ+, λ∈(−∞,0), can be decomposed as the sum of a compact operator plus a sign-definite one.
Proof.
Point 1 is consequence of [22, Lemma 5.14]. Since L1/s(Γ)⊆M(Hs(Γ),H−s(Γ)), s∈[0,1], the map θ:H1/2(Γ)→H−1/2(Γ) is compact by the compact embedding H−1/p(Γ)↪H−1/2(Γ). The difference γ1DLλ+−γ1DLμ is compact for any μ>0 (see [22, Section 5.1.3]) and −γ1DLμ is positive (see [23, Lemma 3.2]). ∎
We consider the self-adjoint operator
ΔΛθ,
[TABLE]
Λzθ is well-defined, i.e., Mzθ has a bounded inverse, by Lemma 4.17. By [9, Theorem 2.19], the map z↦Λzθ and the resolvent formula
(3.5) extend to ZΛθ:=ρ(ΔΛθ)∩C\(−∞,0]. ΔΛθ provides a self-adjoint realization of the (bounded form above) Laplacian on R3\Γ with the semi-transparent boundary conditions at Γ given by γ1u=θ[γ0]u, [γ1]u=0; moreover Theorem 3.6 holds in this case (see [22, Section 5.5]). By point 2 in Lemma 4.17, Corollary 4.7 applies to FλΛθ, λ∈EN− (here EΛθ−=(−∞,0) by [24, Remark 3.8]).
4.1.5. Obstacles with local boundary conditions.
Lemma 4.18**.**
Let b11 and b22 real-valued, b11<0, b11∈L∞(Γ), b11−1∈L∞(Γ), b22∈Lp(Γ), p>2, b12∈Cκ(Γ) for some κ∈(0,1). Then
[TABLE]
[TABLE]
is coercive for any z∈C\R.
Proof.
Given μ>0, let us consider the decomposition Mzb=M(1)+M(2)+M(3), where
and M(1) is sign-definite. M(2) is compact
since both its diagonal elements are compact (here one argues as in the proofs of Lemmata 4.16 and 4.17). Since Cκ(Γ)⊆M(Hs(Γ))⊆L∞(Γ), 0<s<κ, and γ0DLz∈B(H1/2(Γ)), γ1SLz∈B(L2(Γ)) (see, e.g., [26, Theorem 6.12 and successive remarks]), one has that M(3) maps L2(Γ)⊕H1/2(Γ) into
Hs(Γ)⊕L2(Γ) for any s∈[0,1/2]∩[0,κ]; hence it is compact by the compact embeddings Hs(Γ)↪L2(Γ), s>0, and L2(Γ)↪H−1/2(Γ). Therefore Mzb decomposes as the sum of a sign-definite operator plus a compact one. Since, by resolvent identity, Mzb satisfies (3.4), the proof is then concluded by Lemmata 3.2 and 4.7.∎
By Lemma 4.18 and Remark 4.6, the operator-valued map z↦Λzb,
Λzb:=(Mzb)−1, z∈C\R, is well defined and, by (3.4), satisfies (3.3). Therefore, by Theorem 3.3, we can define the self-adjoint operator ΔΛb; it provides a self-adjoint realization of the Laplacian on R3\Γ with boundary conditions
[TABLE]
(see [23, Corollary 4.9]). By [9, Theorem 2.19], the map z↦Λzb and the resolvent formula
(3.5) extend to ZΛb:=ρ(ΔΛb)∩C\(−∞,0]. The choice
[TABLE]
gives ΔΛb=ΔΩinR⊕ΔΩexR, where ΔΩin/exR denotes the Laplacian in L2(Ωin/ex) with Robin boundary conditions γ1in/exuin/ex=bin/exγ0in/exuin/ex (see [23, Section 5.3]). Notice that, since γ1in/ex are both defined in terms of the outward normal vector, the case describing the same Robin boundary conditions at both sides of Γ corresponds to the choice bin=b=−bex (thus b11=2b1, b12=0, b22=−2b).
Arguing as in [24, page 1480], one shows that ΔΛb is bounded from above; moreover ran(Λzb)=L2(Γ)⊕H1/2(Γ) is compactly embedded in K∗=B2,2−3/2(Γ)⊕H−1/2(Γ). Thus Theorem 3.5 applies and the limit operator Λλb,+ exists for any λ∈EΛb− and Λλb,+=(Mλb,+)−1, where
[TABLE]
Proceeding exactly in the same way as in the proof of Lemma 4.18, one shows that Mλb,+
is the sum of a sign-definite operator plus a compact one. Therefore Theorem 4.14 applies to FλΛb, λ∈EDN−∩EΛb−.
5. Inverse Scattering for the Laplace operator with boundary conditions on non-closed Lipschitz surfaces.
We focus now on the case of boundary conditions assigned on a relatively open subset Σ of the boundary Γ of the domain Ω. In this framework ΔΛ provides models of obstacles supported on the non-closed interface Σ; our aim is to determine Σ from the knowledge of the Scattering Matrix by implementing the Factorization Method. An important difference with respect to the previous case appears: in fact the crucial coercivity hypothesis in Theorem 4.14 (by Lemma 4.7, Mλ+ there needs to be coercive) fails to hold in the spaces X♯s, which are made of functions defined on the whole Γ (see Notation 4.1).
To avoid such a problem one introduces (as in [23] and [24]) projectors onto subspaces of functions supported on Σ.
In the following, given X⊂Γ closed, we use the definition
[TABLE]
Given Σ⊂Γ relatively open with a Lipschitz boundary, we denote by ΠΣ the orthogonal projector in the Hilbert space Hs(Γ), ∣s∣≤1, such that ran(ΠΣ) is the subspace orthogonal to HΣcs(Γ).
Lemma 5.1**.**
The orthogonal projection ΠΣ identifies with the restriction map
is an unitary isomorphism. Therefore we can regard Hs(Σ) as a closed subspace of Hs(Γ). Using the decomposition ϕ=(1−ΠΣ)ϕ⊕UΣ−1(ϕ∣Σ), the restriction operator RΣϕ:=0⊕UΣΠΣϕ=0⊕(ϕ∣Σ) is the orthogonal projection from Hs(Γ)≃HΣcs(Γ)⊕Hs(Σ) onto Hs(Σ). Thus, using the identifications ran(ΠΣ)≃Hs(Σ) and HΣ−s(Γ)≃Hs(Σ)∗ (see, e.g., [12, Lemma 4.3.1]), the orthogonal projection ΠΣ identifies with RΣ
and its dual ΠΣ∗ identifies with RΣ∗.
∎
Remark 5.2**.**
Let us notice that if a bounded linear operator M:H−s(Γ)→Hs(Γ) is coercive then RΣMRΣ∗:HΣ−s(Γ)→Hs(Σ) is coercive as well by
[TABLE]
Therefore (see Remark 4.6) (RΣMRΣ∗)−1∈B(Hs(Σ),HΣ−s(Γ)). Moreover, if M=M∘+K with M∘ sign-definite and K compact, then RΣMRΣ∗=RΣM∘RΣ∗+RΣKRΣ∗, with RΣM∘RΣ∗ sign-definite and RΣKRΣ∗ compact. Analogously, if Im⟨ϕ,Mϕ⟩H−s(Γ),Hs(Γ)=0 implies ϕ=0, then Im⟨ϕ,RΣMRΣ∗ϕ⟩HΣ−s(Γ),Hs(Σ)=0 implies RΣ∗ϕ=0 and hence ϕ=0.
The same considerations apply to M:H−s(Γ)⊕H−t(Γ)→Hs(Γ)⊕Ht(Γ) and (RΣ⊕RΣ)M(RΣ∗⊕RΣ∗):HΣ−s(Γ)⊕HΣ−t(Γ)→Hs(Σ)⊕Ht(Σ).
In the following Γ∘ is the Lipschitz boundary of an open bounded set Ω∘⊂R3 and Σ∘⊂Γ∘ is relatively open with Lipschitz boundary.
Lemma 5.3**.**
Let Σ⊂Γ and Σ∘⊂Γ∘ such that R3\(Σ∘∪Σ) is connected. Then
[TABLE]
where
[TABLE]
Proof.
Let uλ,ϕ♯ be the radiating (i.e satisfying the Sommerfeld radiating condition) solution in R3\Σ of Helmholtz equation (−Δ+λ)uλ,ϕ♯=0 with either Dirichlet (whenever ♯=D) or Neumann (whenever ♯=N) boundary condition ϕ∈Hs♯(Σ). Such a solution exists and is unique in
[TABLE]
where uB:=u∣B∩R3\Σ (see [1, Theorems 3.1 and 3.3], see also [2, Section 12.8] and, for the case with smooth boundaries, [32]). Then (see, e.g., [26, Exercise 9.4(iv)]) there exists a unique uλ,ϕ♯,∞∈C∞(S2) such that
[TABLE]
This defines the data-to-pattern operator
[TABLE]
Introducing the Herglotz operators
[TABLE]
where Hλ♯ is defined in (4.3), one has, for any ϕ∈HΣ−s♯(Γ) and f∈L2(S2),
[TABLE]
Proceeding as in [18, proofs of Theorems 1.15 and 1.26] leading to (4.4), one gets
[TABLE]
and so
[TABLE]
By the mapping properties of SLλ+ and DLλ+ and by Remark 5.2, one has RΣγ0SLλ+RΣ∗∈B(HΣs−1/2(Γ),Hs+1/2(Σ))
and RΣγ1DLλ+RΣ∗∈B(HΣs+1/2(Γ),Hs−1/2(Σ)), s∈[0,1/2]. These maps
are bijections (by (5.4), (5.5) in next Subsections 5.1.1 and 5.1.2 and by the regularity results in [10, Theorem 3]; see also [32] for the case of smooth boundaries), and so
[TABLE]
Therefore to conclude the proof we need to show that
[TABLE]
Here we follows the same kind of reasonings as in [19, Section 3.2]. Assume that Σ∘⊂Σ; let uλ♯,Σ∘ be defined according to
[TABLE]
It solves the Helmoltz equation (−Δ+λ)uλ♯,Σ∘=0 in R3\Σ∘ and hence in R3\Σ as well.
Let ϕΣ∘D:=RΣγ0uλD,Σ∘∈H1/2(Σ), ϕΣ∘N:=RΣγ1uλN,Σ∘∈H−1/2(Σ). Then Kλ♯ϕΣ∘♯=ϕλΣ∘. Suppose now that Σ∘∩Σc=∅. Let B⊂R3 be an open ball such that B∩Σ=∅, B∩Σ∘=∅. Assume that ϕλΣ∘=Kλ♯ϕ♯ for some ϕ♯∈Hs+s♯(Σ) and consider the corresponding radiating solution uλ,ϕ♯♯. Then, since Kλ♯ϕ♯=Kλ∘,♯ϕΣ∘∘,♯ (here the apex ∘ denotes objects defined by using the surface Γ∘), one has, by Rellich’s Lemma and unique continuation, uλ,ϕ♯♯∣R3\(Σ∘∪Σ)=uλ∘,♯,Σ∘∣R3\(Σ∘∪Σ). By elliptic regularity,
(−Δ+λ)uλ,ϕ♯∣B=0 implies
uλ,ϕ♯∣B∈H2(B); this leads to a contradiction, since uλ∘,♯,Σ∘∣B∈/H2(B).
∎
By the same kind of proof provided for Corollary 4.4, one gets the following:
Corollary 5.4**.**
Let Σ⊂Γ and Σ∘⊂Γ∘ such that R3\(Σ∘∪Σ) is connected. Then
[TABLE]
Notation 5.5**.**
We introduce the spaces
[TABLE]
so that
[TABLE]
The following results is the analogue for screens of Theorem 4.14 :
Theorem 5.6**.**
Let
[TABLE]
and suppose that Λλ+=RΣ∗(Mλ+,Σ)−1RΣ, where the bijection Mλ+,Σ∈B((X♯s)∗,X♯s), s∈[0,1/2], has the decomposition Mλ+,Σ=M∘+,Σ+Kλ+,Σ, where M∘+,Σ is sign-definite and Kλ+,Σ is compact. Let Σ∘⊂Γ∘ such that R3\(Σ∘∪Σ) is connected; then
[TABLE]
where the sequences {zλ,kΛ}1∞⊂C\{0} and {ψλ,kΛ}1∞⊂L2(S2) provide the spectral resolution of FλΛ as in Remark 4.11 and ϕλΣ∘ is defined in (5.1).
Proof.
We use the factorization F^{\Lambda}_{\lambda}=\big{(}L^{\sharp}_{\lambda}R^{*}_{\Sigma}(M^{+,\Sigma}_{\lambda})^{-1}\big{)}(M^{+,\Sigma}_{\lambda})^{*}\big{(}L^{\sharp}_{\lambda}R^{*}_{\Sigma}(M^{+,\Sigma}_{\lambda})^{-1}\big{)}^{*}.
By proceeding as in the proof of Theorem 4.14 (where now Bλ=Lλ♯RΣ∗(Mλ+,Σ)−1), one gets Im⟨ϕ,Mλ+,Σϕ⟩(X♯s)∗,X♯s=0 for any ϕ=0. Since Mλ+,Σ is a bijection, one has
ran(Lλ♯RΣ∗(Mλ+,Σ)−1)=ran(Lλ♯∣(X♯s)∗). Hence, by Lemma 5.3, Corollary 5.4 and by [18, Theorem 1.16],
[TABLE]
By proceeding as in the proof of Theorem 4.13, ran(Lλ♯RΣ∗(Mλ+,Σ)−1)=ran(∣FλΛ∣1/2) and then one concludes by the same arguments.
∎
Remark 5.7**.**
If Mλ+,Σ in Theoren 5.6 is merely coercive, then the “inf” criterion still holds.
5.1. Applications
Here we apply Theorem 5.6 the analogue of models in the examples considered in Section 4.1 where now the boundary conditions holds only on Σ. Before considering the specific examples, let us explain our strategy.
At first, notice that all the examples in Section 4.1 consider self-adjoint operators ΔΛ with Λz=Mz−1, where the map z↦Mz satisfies (3.4) (see Remark 3.1). Hence, by Lemma 3.2, Im⟨ϕ,Mzϕ⟩X♯s∗,X♯s=0 for any ϕ=0. Furthermore, all such Mz’s have a decomposition Mz=M∘+Kz with M∘ sign-definite, Kz compact; this can be checked by proceeding as in the proof of Lemma 4.18 using identities (2.14). Then, by Lemma 4.7, Mz is coercive. Now, the dual couple of projectors RΣ, RΣ∗ in Lemma 5.1 come into play: by Remark 5.2, these properties of Mz transfer to MzΣ:=RΣMzRΣ∗ (here and in the following lines, RΣ has to be replaced by RΣ⊕RΣ when one considers example in Subsection 4.1.5), and so, in particular, MzΣ is coercive and (RΣMzRΣ∗)−1∈B(X♯s,(X♯s)∗). Then, setting
[TABLE]
it is immediate to check that z↦Λz∈B(X♯s,X♯s∗) satisfies (3.3), and so, by Theorem 3.3, it defines a self-adjoint operator ΔΛ. Such an operator describes the model corresponding to the same kind of boundary conditions associated to ΔΛ, now assigned only on Σ (see [23, Section 6], [24, Section 7]). Since the limit operator Mλ+ exists (use (2.13)) and, by Theorem 3.6, the limit Λλ+ exists as well, one gets Λλ+=RΣ∗(RΣMλ+RΣ∗)−1RΣ. Now, since all Mλ+ appearing in the examples in Section 4.1 decompose as the sum of a sign-definite operator plus a compact one, by Remark 5.2 the same is true for RΣMλ+RΣ∗.
In conclusion, the assumptions in Theorem 5.6 hold for any Λ defined as in (5.3) where Mz is any of the operators given in the examples in Section 4.1; hence the reconstruction formula (5.2) applies to FλΛ. In what follows this scheme is implemented case by case.
5.1.1. Dirichlet screens.
One considers ΔΛD with
[TABLE]
ΔΛD is a (bounded from above) self-adjoint representation of the Laplacian on R3\Σ with homogeneous Dirichlet boundary conditions at Σ (see [24, Example 7.1]).
By [9, Theorem 2.19], the map z↦ΛzD and the corresponding resolvent formula (3.5) extends to ZΛD:=ρ(ΔΛD)∩C\(−∞,0]=C\(−∞,0]. By [24, Theorem 3.7], σp−(ΔΛD) is empty and so, by Theorem 3.6,
ΔΛN is a (bounded from above) self-adjoint representation of the Laplacian on R3\Σ with homogeneous Neumann boundary conditions at Σ (see [24, Example 7.2]). By [9, Theorem 2.19], the map z↦ΛzN and the corresponding resolvent formula (3.5) extends to ZΛN:=ρ(ΔΛN)∩C\(−∞,0]=C\(−∞,0]. By [24, Theorem 3.7],
σp−(ΔΛN) is empty and so, by Theorem 3.6,
5.1.3. Screens with semitransparent boundary conditions αΣγ0u=[γ1]u.
Let α∈L∞(Γ) real-valued such that sgn(α) is constant and α1∈L∞(Γ); let us define αΣ:=α∣Σ and αΣ−1∈B(LΣ2(Γ),L2(Σ))
by αΣ−1ϕ:=(α−1ϕ)∣Σ. Since −(αΣ−1+RΣγ0SLzRΣ∗)=RΣMzαRΣ∗, where Mzα is defined in (4.13), one considers ΔΛα, where
[TABLE]
ΔΛα is a self-adjoint representation of the (bounded from above) Laplacian on R3\Σ with boundary conditions at Σ given by αΣRΣγ0u=RΣ[γ1]u, RΣ[γ0]u=0 (see [24, Example 7.4]). By [9, Theorem 2.19], the map z↦Λzα and the resolvent formula (3.5) extend to ZΛα:=ρ(ΔΛα)∩C\(−∞,0]. By [24, Theorem 3.7], σp−(ΔΛα) is empty and so, by Theorem 3.6,
5.1.4. Screens with semitransparent boundary conditions θΣγ1u=[γ0]u.
Let θ∈Lp(Γ), p>2; set θΣ:=θ∣Σ and define the corresponding operator θΣ∈B(HΣ1/2(Γ),H−1/2(Σ)) by θΣφ:=(θφ)∣Σ. Since (θΣ−RΣγ1DLzRΣ∗)=RΣMzθRΣ∗, where Mzθ is defined in (4.14), one considers ΔΛθ, where
[TABLE]
ΔΛθ is a self-adjoint representation of the (bounded from above) Laplacian on R3\Σ with boundary conditions at Σ given by θΣRΣγ1u=RΣ[γ0]u, RΣ[γ1]u=0 (see [24, Example 7.5]). By [9, Theorem 2.19], the map z↦Λzθ and the resolvent formula (3.5) extend to ZΛθ:=ρ(ΔΛθ)∩C\(−∞,0]. By [24, Theorem 3.7], σp−(ΔΛθ) is empty and so, by Theorem 3.6,
Let b11∈L∞(Γ), b11−1∈L∞(Γ), b22∈Lp(Γ), p>2, b12∈Cκ(Γ), 0<κ<1, with b11 and b22 real-valued, b11<0. Set bijΣ:=bij∣Σ and define the corresponding multiplication operator by bijΣϕ:=(bijϕ)∣Σ, where supp(ϕ)⊆Σ.
Since
[TABLE]
where Mzb is defined in Lemma (4.18), one considers ΔΛb,
where
[TABLE]
ΔΛb is a self-adjoint representation of the (bounded from above, this follows proceeding as in [24, page 1480]) Laplacian on R3\Σ with boundary conditions at Σ given by
[TABLE]
By [9, Theorem 2.19], the map z↦Λzb and the resolvent formula (3.5) extend to ZΛb:=ρ(ΔΛb)∩C\(−∞,0]. Since σp−(ΔΛb) is empty (see [24, Theorem 3.7]),
[TABLE]
exists for any λ∈(−∞,0) by Theorem 3.6. Therefore Theorem 5.6 applies to FλΛb, λ∈(−∞,0).
Bibliography35
The reference list from the paper itself. Each links out to its DOI / PubMed record.
1[1] M. S. Agranovich: Strongly elliptic second order systems with spectral parameter in transmission conditions on a nonclosed surface. In: P. Boggiatto, L. Rodino, J. Toft, M. W. Wang (eds.), Pseudo-differential operators and related topics , Operator Theory: Advances and Applications, vol. 164, 1-21, Birkhäuser, Basel, 2006.
2[2] M. S. Agranovich: Sobolev Spaces, Their Generalizations and Elliptic Problems in Smooth and Lipschitz Domains. Springer, Berlin, 2015.
3[3] H. Baumgärtel, M. Wollenberg: Mathematical Scattering Theory , Akademie-Verlag, Berlin, 1983.
4[4] J. Behrndt, M. Langer, V. Lotoreichik: Schrödinger operators with δ 𝛿 \delta - and δ ′ superscript 𝛿 ′ \delta^{\prime} -potentials supported on hypersurfaces. Ann. Henri Poincaré 14 (2013), 385-423.
5[5] J. Behrndt, M. Malamud, H. Neidhardt: Scattering matrices and Dirichlet-to-Neumann maps. J. Funct. Anal. 273 (2017), 1970-2025.
6[6] O. Bondarenko, X. Liu: The factorization method for inverse obstacle scattering with conductive boundary condition. Inverse Problems 29 (2013), 095021 (25pp).
7[7] Y. Boukari, H. Haddar: The factorization method applied to cracks with impedance boundary conditions. Inverse Probl. Imaging 7 (2013), 1123-1138.
8[8] J.F. Brasche, P. Exner, Y.A. Kuperin, P. Šeba: Schrödinger operators with singular interactions. J. Math. Anal. Appl. 184 (1994), 112-139.