This paper establishes uniqueness in inverse acoustic scattering problems involving media with unbounded gradient across Lipschitz surfaces, by linking it to Schrödinger operators with singular delta potentials.
Contribution
It introduces a novel uniqueness result for inverse scattering with media having unbounded gradient, extending previous theories to include singular surface-supported potentials.
Findings
01
Proves uniqueness in inverse acoustic scattering with unbounded gradient.
02
Establishes a connection between acoustic scattering and Schrödinger operators with delta potentials.
03
Extends inverse scattering theory to more singular media configurations.
Abstract
We prove uniqueness in inverse acoustic scattering in the case the density of the medium has an unbounded gradient across Σ⊆Γ=∂Ω, where Ω is a bounded open subset of R3 with a Lipschitz boundary. This follows from a uniqueness result in inverse scattering for Schr\"odinger operators with singular δ-type potential supported on the surface Γ and of strength α∈Lp(Γ), p>2.
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
Uniqueness in inverse acoustic scattering with unbounded gradient across
RICAM, Austrian Academy of
Sciences, Altenbergerstr. 69, A-4040 Linz, Austria
Abstract
We prove uniqueness in inverse acoustic scattering in the case the density of
the medium has an unbounded gradient across Σ⊆Γ=∂Ω, where Ω is a bounded open subset of R3
with a Lipschitz boundary. This follows from a uniqueness result in
inverse scattering for Schrödinger operators with singular
δ-type potential supported on the surface Γ and of strength
α∈Lp(Γ), p>2.
The aim of this paper is the study of the uniqueness problem in the inverse scattering for the acoustic wave equation
[TABLE]
in the case ϱ has an unbounded gradient across some surface Σ⊆Γ=∂Ω, where Ω⊂R3 is open and bounded with Lipschitz boundary.
Here u is the pressure field, ϱ is the density and v is the sound speed; we assume that ϱ(x)=v(x)=1 whenever x lies outside some large ball BR⊃Ω.
To introduce our arguments and to allow the reasoning in the following lines, we start assuming that the functions ϱ and v are positive and sufficiently regular (for instance we can take ϱ of class C2 and v bounded). Looking for fixed frequency solutions of the kind u(t,x)=e−iωtuω(x), ω>0, one gets the stationary equation
[TABLE]
Defining u~ω:=ϱ−1uω,
the equation (1.1) transforms into
[TABLE]
where Hφ,v,ω denotes the Schrödinger operator
[TABLE]
[TABLE]
Notice that, since ϱ=v=1 outside BR, the potential Vφ,v,ω is compactly supported.
As well known from stationary scattering theory in quantum mechanics, whenever V is a short-range potential, a generalized eigenfunction for the corresponding Schrödinger operator,
(−Δ+V)ψk=k2ψk, k>0, admits the outgoing representation
[TABLE]
where sV(k,ξ^,ξ^′), ξ^,ξ^′∈S2, denotes the scattering amplitude (see e.g. [3, page 425]).
Since the solution uω of equation (1.1) and the solution uω of the corresponding quantum scattering problem (1.3)-(1.4) identify outside a ball, i.e.: uω(x)=u~ω(x) for ∣x∣>R, the above representation yields the asymptotic formula
[TABLE]
where the far-field pattern uϱ,v∞ is related to the scattering amplitude by the equality
[TABLE]
The inverse acoustic scattering problem consists in recovering the couple of functions (ϱ,v) from the knowledge of the far-field pattern at some fixed frequencies; in particular, to recover the two independent functions ϱ and v, one needs the knowledge of the far-field patters at least for two different frequencies ω and ω~. Clearly, the solvability of such an inverse problem requires a corresponding uniqueness result:
[TABLE]
By (1.5), this uniqueness issue is a consequence of an analogous result concerning Schrödinger operators:
[TABLE]
The justification of the uniqueness property (1.6) goes back to the pioneering works [29, 37, 32]. The idea is based on the orthogonality relation ∫BR(V1(x)−V2(x))u1(x)u2(x)dx=0,
involving the total field solutions uj to the Schrödinger equation with potential Vj, and which can be derived from the equality of the far-field patterns for the two frameworks. Then the strategy consists in constructing a specific set of solutions uj, known as complex geometrical optics solutions or CGO’s in short, and use them to deduce that V1=V2 (here stands for the Fourier transform). Finally, the two equalities Vφ1,v1,ω=Vφ2,v2,ω and Vφ1,v1,ω~=Vφ2,v2,ω~ entail (ϱ1,v1)=(ϱ2,v2).
The aim of our work is to extend the above reasoning and conclusions to the case in which the density function ϱ belongs to Hloc1(R3) and the jump of its normal derivative across some closed set Σ⊂Γ belongs to Lp(Γ), p>2, where Γ is the Lipschitz boundary of some opened and bounded Ω⊂R3 (see Section 7 for the precise hypotheses and statements).
Under these conditions, the corresponding Schrödinger equation is modeled by a potential of the form (1.3)-(1.4) with and additive δ-type potential supported on
Γ with a strength belonging to Lp(Γ),p>2 (see Section 3). Hence, the inverse problem consists in extending the above approach (holding for regular perturbations) to the case of Schrödingers operators
with singular δ-type potentials supported on Γ.
As the setting of the problem is motivated by many applications in sciences and engineering, after those mentioned works, a considerable effort was put to improve and refine these results to deal with potentials in more
general classes of functions and also other models as the electromagnetism and elasticity for instance. The reader can see the following references for more information [11, 21, 33, 38].
A model of particular interest is the EIT (Electrical Impedance Tomography) problem, also called Calderón’s problem, which consists in identifying the conductivity σ using Cauchy data
(u\arrowvert∂Ω,σ∇u⋅ν\arrowvert∂Ω) of the solution of ∇⋅σ∇u=0, in Ω⊂R3.
The uniqueness question of this problem is reduced, in the same way as described above, to the construction of the CGO’s, see [37], where σ is a positive C2-smooth function.
The regularity of σ is reduced to C23+ϵ in [6],
then to C23,∞ in [30] and to C23,p,p≥6 in [7]. Finally in [17, 10] this condition is reduced to W1,∞ and then to W1,3 in [18] where
the CGO’s are constructed allowing potentials of the form ∇⋅f+h, where f∈L3 and h∈L23 with compact supports.
This last result is a key for us as δ-type potentials, with strengths in Lp(Γ)p>2, can be cast in these forms (see Section 6).
In particular, using CGO techniques, the analysis developed in Section 6 and Section 7 provides with a uniqueness result for the case of positive and bounded acoustic densities ϱ which are in Hloc1(R3) and
such that ∣∇Ωin/exϱ∣∈L4(Ωin/ex) and ΔΩin/exϱ∈L2(Ωin/ex), where Ωin/ex denotes the interior or the exterior of Ω, while the normal derivatives across a closed subset Σ of Γ have jumps of regularity Lp(Σ) with p>8/3 (see Theorem 7.4 and Remark 7.3 for the details).
Let us now discuss the forward problem and how we model the acoustic scattering with such regularity of the density.
There are several ways to study and describe the solutions of the forward acoustic scattering and generate the far-field patterns. We mention the variation formulation, see [20, 8] for instance, which
reduces the problem to a bounded domain Ω by introducing a Dirichlet-Neumann map to the exterior problem, i.e. stated in R3∖Ω, where the background is homogeneous.
A second approach consists in using integral equations; this allows to reduce the problem to inverting a Lippmann-Schwinger equation via the Fredholm alternative, see [24]. The approach requires, in addition to the regularity of the coefficients, a positivity of the contrast, i.e. in our case v2ρ=const. and ρ<1, see [24].
In this paper we follow a different strategy and exploit the connection between the acoustic problem and the Schrödinger one, providing the link between (1.1) and (1.2) in the case the density ϱ is no more C2 as supposed in the reasonings above. Due to the lack of regularity of ϱ, we use Schröedinger operators with δ-type potentials and unbounded strengths, thus generalizing previously known results about such kind of operators (see e.g. [5], [27] and references therein); for this class of operators we provide the rigorous construction as self-adjoint extensions of the symmetric operator Δ∣Ccomp∞(R3\Γ). The Schrödinger approach allows the use of techniques from quantum mechanical stationary scattering theory, in particular, by extending some results provided in [27], we get a limiting absorption principle (LAP for short in the following) for our class of Schrödinger operators; as a consequence, the scattering amplitude is derived and used to define the acoustic far-field patterns. Let us remark that, by combining the results contained in [34] with [13, Theorem 16], one could get a non-stationary scattering theory (i.e. the existence of the wave operators) directly for the acoustic model whenever 0<c1≤ϱ,v≤c2<+∞. Nevertheless, using the connection with Schrödinger operators, and the corresponding LAP, our approach has the advantage of easily providing with the acoustic far-field patterns in terms of the (quantum mechanical) scattering amplitude and results better suited for the study of the inverse scattering problem.
The paper is organized as follows. The self-adjoint realizations of such operators are provided in Section 2 and the existence of a limiting absorption principle for them
is given in Section 4. The proof of the connection between Schrödinger operators with δ-type potentials and acoustic operators with densities with unbounded gradients is provided in Section 3.
In Section 5, we give sense to the far-field through the construction of the generalized eigenfunctions. In section 6, we derive the uniqueness result for the Schrödinger model, as Theorem 6.2,
and then we conclude the corresponding result for the acoustic model, as Theorem 7.4, in Section 7.
2 Schrödinger operators with delta interactions of unbounded strength.
Let V∈L2(R3)+L∞(R3); then, by the Kato-Rellich theorem,
[TABLE]
is self-adjoint and bounded from above. Here Hs(R3), s∈R, denotes the scale of Sobolev Spaces on R3, we refer to [16, Chapter 1] for the appropriate definitions of such spaces and for the trace maps defined on them; we also refer to the same book for the definition of the scale Hs(Γ), ∣s∣≤1, of Sobolev spaces on the Lipschitz surface Γ which we use below.
AV can be broadened to an operator in H−2(R3) (by a slight abuse of notation we denote such an operator with the same symbol):
[TABLE]
where now V denotes the linear operator, belonging to B(L2(R3),H−2(R3)) by the Kato-Rellich hypothesis, defined by
[TABLE]
Since AV∈B(H2(R3),L2(R3)), by duality and interpolation one has
[TABLE]
and, setting RzV:=(−AV+z)−1, z∈ρ(AV),
[TABLE]
Lemma 2.1**.**
Let dzV denote the distance of z∈ρ(AV) from σ(AV).
Then there exists cV>0 such that, whenever dzV>cV,
[TABLE]
Proof.
By RzV=Rz0(1+VRz0)−1 and ∥VRz0∥B(L2(R3))→0 as ∣z∣→+∞, one gets, whenever dzV>cV,
[TABLE]
Thus, since
[TABLE]
by interpolation one obtains
[TABLE]
and, by duality,
[TABLE]
The proof is then concluded by interpolation again.
∎
Given Ω⊂R3, open and bounded with Lipschitz boundary Γ, we
introduce the bounded and surjective trace map
[TABLE]
defined as the unique bounded extension of the map
[TABLE]
In the following we also use the extension (denoted by the same symbol) of γ0 to Hlocs+21(R3) defined by γ0u:=γ0(χu), where χ∈Ccomp∞(R3) and χ=1 on an open neighborhood of Γ.
Using the adjoint γ0∗:H−s(Γ)→H−s−21(R3) and RzV we define the bounded operator (the single-layer potential)
[TABLE]
This gives the bounded operator
[TABLE]
Remark 2.2**.**
Given ϕ∈H2(R3) and ξ∈Hs(Γ), ∣s∣≤1, let ψ:=ϕ−SLzVξ. By the definition of SLzV one has (−AV+z)ψ=(−AV+z)ϕ−γ0∗ξ. Thus, notwithstanding neither AVψ nor γ0∗ξ belong to L2(R3), one has
[TABLE]
Lemma 2.3**.**
Let α∈B(Hs(Γ),H−s(Γ)), 0<s<21. Then there exists cα,V>0 such that forall z∈C such that dzV>cα,V one has (1+γ0SLzVα)−1∈B(Hs(Γ)).
Such an inequality show that if 0<s<21 then there exists cα,V>0 such that operator norm ∥γ0SLzVα∥B(Hs(Γ)) is strictly smaller than one whenever dzV>cα,V.
∎
Corollary 2.4**.**
Let α∈B(Hs(Γ),H−s(Γ)), 0<s<21 such that α∗=α. Then there exists a finite set Sα,V⊂R such that (1+αγ0SLzV)−1∈B(H−s(Γ)) for any z∈ρ(AV)\Sα,V. Moreover
[TABLE]
Proof.
Let 0<s<21. By the compact embedding H1−s(Γ)↪Hs(Γ) and by ran(γ0SLzV)⊆H1−s(Γ), the map
γ0SLzV:H−s(Γ)→Hs(Γ) is compact and so γ0SLzVα:Hs(Γ)→Hs(Γ) is compact as well.
Moreover, by the identity SLzV=SLwV+(w−z)RzVSLwV, the map z↦γ0SLzVα is analytic from ρ(AV) to B(Hs(Γ)).
Thus, since the set of z∈ρ(AV) such that (1+γ0SLzVα)−1∈B(Hs(Γ)) is not void by Lemma 2.3, by analytic Fredholm theory (see e.g. [35, Theorem XIII.13]), (1+γ0SLzVα)−1∈B(Hs(Γ)) for any z∈ρ(AV)\Sα,V, where Sα,V is a discrete set. By next Theorem 2.5, Sα,V is contained in the spectrum of a self-adjoint operator and so Sα,V⊂R; hence, by Lemma 2.3, Sα,V⊆[supσ(AV),supσ(AV)+cα,V] and so it is finite being discrete, i.e. without accumulation points.
Since (1+γ0SLzˉVα)∗=(1+γ0RzˉVγ0∗α)∗=1+αγ0RzVγ0∗=1+αSLzV and 1+γ0SLzˉVα is surjective, 1+αγ0SLzV is injective and hence invertible for any z∈ρ(AV)\Sα,V. Moreover
[TABLE]
By the obvious equality (1+αγ0SLzV)α=α(1+γ0SLzVα), one gets (1+αγ0SLzV)−1α=α(1+γ0SLzVα)−1 and so
[TABLE]
∎
By the previous results one has
[TABLE]
Thus
[TABLE]
and
[TABLE]
is a well-defined family of bounded operators in L2(R3).
Taking λ∘∈R∩ρ(AV), in the following we use the shorthand notation SL∘V≡SLλ∘V.
Theorem 2.5**.**
Let V∈L2(R3)+L∞(R3) and α∈B(Hs(Γ),H−s(Γ)), α=α∗, 0<s<21. The family of bounded linear operators RzV,α given in (2.3)
is the resolvent of the self-adjoint operator AV,α in L2(R3) defined, in a λ∘-independent way, by
[TABLE]
[TABLE]
Proof.
We proceed as in the proof of [31, Theorem 2.1]. Setting Λz:=(1+αγ0SLzV)−1α, using the resolvent identity for RzV and definition (2.3), one gets, for any w,z∈ZV,α (see the explicit computation in
[31, page 115])
[TABLE]
By SLzV=RzVγ0∗ and resolvent identity for RzV, it results
Therefore RzV,α is a pseudo-resolvent. Moreover,
RzV,α is injective, since, if ψ∈ker(RzV,α) then
[TABLE]
This gives RzVψ=0 and so ψ=0. Hence, see e.g. [23, Chap. VIII, Section 1.1], RzV,α is the
resolvent of a closed operator A^V,α and the identity (2.1) implies
[TABLE]
so that such an operator is self-adjoint; given z∘∈ZV,α, A^V,α is defined, in a z∘-independent way, by
[TABLE]
Notice that any ψ∈dom(A^V,α) is given
by
[TABLE]
By the mapping properties of SLz and by ran(Λz)⊆H−s(Γ), one gets
dom(A^V,α)=ran(Rz∘V,α)⊆H23−s(R3). Thus
Let ψ=ψz∘−SLz∘VΛz∘γ0ψz∘∈dom(A^V,α). Since
[TABLE]
one has ψ=ψz∘−SLz∘Vαγ0ψ. Then
[TABLE]
and so ψ∈dom(AV,α). Conversely, given ψ∈dom(AV,α), define ψz∘:=ψ+SLz∘γ0ψ.
Then, by (2.10), ψ=ψz∘+SLz∘Λz∘γ0ψz∘ and, by (2.11), ψz∘∈H2(R3). Thus
ψ∈dom(A^V,α) and so dom(A^V,α)=dom(AV,α).
By (2.9),
[TABLE]
Finally,
[TABLE]
shows that the definition of dom(AV,α) is λ∘-independent.
∎
Remark 2.6**.**
A particular case of operator α∈B((Hs(Γ),H−s(Γ)), such that α=α∗ is α∈M(Hs(Γ),H−s(Γ)), α real-valued, where M(Hs(Γ),H−s(Γ)) denotes the set of Sobolev multipliers on Hs(Γ) to H−s(Γ) (here and in the following we use the same notation for a function and for the corresponding multiplication operator).
By proceeding as in the proof of Theorem 2.5.3 in [19], one has
[TABLE]
Then, by Sobolev’s embeddings and Hölder’s inequality, one gets
[TABLE]
Thus we can define AV,α whenever α∈Lp(Γ), p>2.
Remark 2.7**.**
One can check that, in the particular cases where V∈L∞(R3), α∈L∞(Γ) and Γ is smooth, the self-adjoint operators AV,α coincide with the ones studied (and constructed by different methods) in [5, Section 3.2]; also see [26, Section 5.4] for a construction that follows the lines here employed in the case V∈Cb∞(R3), α∈M(H23(Γ)) and Γ is of class C1,1. Similar kind of operators in the case Γ is not necessarily Lipschitz and can have a not integer dimension have been considered in [31, Example 3.6]
Remark 2.8**.**
Here and below we use dualities ⟨⋅,⋅⟩X∗,X which are conjugate linear with respect to the first variable. Let ξ∈H−s(Γ), 0<s≤1. Since
[TABLE]
for all ϕ∈Hs+1/2(R3), the distribution γ0∗ξ has support contained in Γ. In the case ξ∈L2(Γ) one has
[TABLE]
where σΓ denotes the surface measure. In particular γ0∗1=δΓ, where δΓ denotes the Dirac distribution supported on Γ. Introducing the notation γ0∗ξ≡ξδΓ, the operator AV,α is represented as AV,αψ=AVψ−αγ0ψδΓ and this explain why this kind of operators are said to describe quantum mechanical models with singular, δ-type interactions.
Remark 2.9**.**
Notice that AV,α is a self-adjoint extension of the symmetric closed operator AV∣ker(γ0). If α∈M(Hs(Γ),H−s(Γ)) then supp(γ0∗αγ0ψ)⊆Σα, Σα:=supp(α), and so (AV,αψ)∣Σαc=(AVψ)∣Σαc. This shows that AV,α is a self-adjoint extension of the symmetric operator AV∣Ccomp∞(R3\Σα) and so it depends only on Σα and not on the whole Γ: outside Σα we can change Γ at our convenience without modifying the definition of AV,α.
Lemma 2.10**.**
Under the assumptions of Theorem 2.5,
the self-adjoint operator AV,α is bounded from above and σess(AV,α)=(−∞,0]. Moreover, if V is compactly supported and R3\Ω is connected then σp(AV,α)∩(−∞,0)=∅.
Proof.
By V∈L2(R3)+L∞(R3) and by the Kato-Rellich theorem, AV is bounded from above. Thus, by (2.2), there exists λV>sup(σ(AV)) such that: λ∈ZV,α for
all λ>λV. Then, the resolvent formula (2.3)
implies (λV,+∞)⊂ρ(AV,α) and so AV,α is bounded from above.
By Corollary 2.4, by the compact embedding H−s(Γ)↪H−1(Γ) and by (2.3), the resolvent difference RzV,α−RzV is a compact operator. Therefore, since σess(AV)=(−∞,0] (see e.g. [35, Example 6, Section 4, Chapter XIII]), one has σess(AV,α)=σess(AV)=(−∞,0].
Let us now suppose that supp(V) is compact and that exists λ∈σp(AV,α)∩(−∞,0); let ψλ denote a corresponding eigenvector. Let K a compact set containing both Γ and supp(V), so that
(−Δψλ+λψλ)∣Kc=0; by elliptic regularity, ψλ∈C∞(Kc), and, by the Rellich estimate
one gets ψλ∣Kc=0 (see e.g. [25, Corollary 4.8]).
Using the unique continuation property (holding for our exterior problem in R3\Ω according to [22]), we get ψλ∣R3\Ω=0. Since ψλ∈dom(AV,α)⊆H23−s, this gives γ0ψλ=0 and so ψλ∈H2(R3) and (−AV+λ)ψλ=0, i.e. λ∈σp(AV). This contradicts σp(AV)∩(−∞,0)=∅ (which holds for any V∈Lcomp3/2(R3), see [22]).
∎
The next lemma shows that the construction leading to Theorem 2.5 is unaffected by the addition of a bounded potential:
Lemma 2.11**.**
Let V and α be as in Theorem 2.5. If V∞∈L∞(R3) then
[TABLE]
Proof.
According to the representation(2.5), we only need to show that dom(AV+V∞,α)=dom(AV,α). By the definition of SLzV and the second resolvent identity there follows
[TABLE]
Since RzVV∞∈B(L2(R3),H2(R3)), then (2.4) yields the sought domains equality.
∎
3 The connection between acoustic and Schrödinger operators.
We begin the section by reviewing some results about multiplication of distributions and related topics.
Given the couple u∈Hloct(R3), v∈H−s(R3), 0≤s≤t, we can define the product uv∈D′(R3) by
[TABLE]
In particular, the product u(γ0∗ξ)∈D′(R3) is well defined for any
ξ∈H−s(Γ), 0<s≤1, and u∈Hloct(R3), t≥s+21.
Lemma 3.1**.**
If u∈Hloct(R3), v∈H−s(R3), 1≤s+1≤t, then
[TABLE]
Proof.
[TABLE]
∎
Remark 3.2**.**
Notice, that, by the same proof, (3.1) holds true also in the case u∈Hcomp1(R3) and v∈Lloc2(R3).
Lemma 3.3**.**
If u,v∈H1(R3) then uv∈W1,1(R3) and γ0(uv)=γ0uγ0v in L1(Γ).
Proof.
By (3.1) , uv∈W1,1(R3). Since γ0∈B(W1,1(R3),L1(Γ)) one has γ0(uv)∈L1(Γ). Let {un}1∞⊂D(R3), {vn}1∞⊂D(R3) such that un→u and vn→v in H1(R3). Thus, by (3.1) , unvn→uv in W1,1(R3). Since γ0∈B(H1(R3),H21(Γ)), γ0(unvn)=γ0unγ0vn converges in L1(Γ) to both γ0(uv) and γ0uγ0v.
∎
Since W1,∞(Γ)⊆M(Hs(Γ)), 0≤s≤1, we can define the product ζξ∈W1,∞(Γ)′ whenever ζ∈Ht(Γ) and ξ∈H−s(Γ), 0≤s≤t≤1, by
[TABLE]
Notice that the inclusion W1,∞(Γ)⊂H1(Γ) implies H−s(Γ)⊂W1,∞(Γ)′, with 0≤s≤1.
Since γ0ϕ∈W1,∞(Γ) whenever ϕ∈D(R3), given ξ∈W1,∞(Γ)′ one defines γ0∗ξ∈D′(R3) by
[TABLE]
In the case ξ∈H−s(Γ), 0<s≤1, the mapping properties of γ0 imply
γ0∗ξ∈H−s−21(R3); then, from the above identity one recovers the preceding definition in term of the dual map of the trace γ0.
Lemma 3.4**.**
If ξ∈H−s(Γ), 0<s≤1, and u∈Ht+21(R3), t≥s, then
[TABLE]
Proof.
[TABLE]
∎
Lemma 3.5**.**
Let u∈Hloc1(R3) such that u1∈L∞(R3). Then u1∈Hloc1(R3) and
[TABLE]
Proof.
Since u1∈L∞(R3), the definition of the distributional gradient
[TABLE]
shows that ∇u1∈(W1,1(R3))′=W−1,∞(R3). Thus, for any v∈Wloc1,1(R3), we can define the product v∇u1∈D′(R3) by
[TABLE]
Since u∈Hloc1(R3)⊂Wloc1,1(R3), by
[TABLE]
we get
[TABLE]
i.e. u∇u1=−u∇u. Let χ∈Ccomp∞(R3) such that χ=1 on an open neighborhood of Γ; by Lemma 3.3, 1=γ0(χuχu1)=γ0(χu)γ0(χu1)=γ0uγ0u1. Thus γ0u is a.e. different from zero and γ0u1=γ0u1.
∎
Given the real-valued function φ we suppose there exists an open and bounded set Ωφ≡Ω⊂R3 with Lipschitz boundary Γφ≡Γ such that
[TABLE]
where Ωin≡Ω, Ωex≡R3\Ω. Let n(x) denote the exterior unit normal at x∈Γ; the lateral operators defined in Ccomp∞(Ωin/ex) by
[TABLE]
uniquely extend to bounded maps
[TABLE]
Furthermore, by [28, Lemma 4.3 and Theorem 4.4], these extend to
[TABLE]
[TABLE]
as bounded operator with respect to the natural norm
[TABLE]
Therefore the jump across Γ given by
[TABLE]
where χ∈Ccomp∞(R3) is such that χ=1 on an open neighborhood of Γ, is a well-defined distribution in H−21(Γ). Moreover, by Lemma 3.5, γ0u1∈H21(Γ) and its product with [γ^1]φ is well-defined in W1,1(Γ)′. As further assumption, beside (3.2), we suppose
[TABLE]
In particular, by Remark 2.6, hypothesis (3.3) holds true whenever
[TABLE]
Remark 3.6**.**
A more explicit characterization of a class of function φ satisfying hypotheses (3.2) and (3.3) is the following:
[TABLE]
where φ∘∈Hloc2(R3) and ξ∈Lp(Γ), p>2. By the properties of the single layer potential SL (see [14, Theorem 3.1]), one has
[TABLE]
for any ϵ>0 and any χ∈Ccomp∞(R3). Since
W1+1/p−ϵ,p(Ωin/ex)⊂C(Ωin/ex) whenever p>2 and ϵ is sufficiently small, one gets φ∈C(R3) and so φ−1∈L∞(R3) entails (γ0φ)−1∈L∞(Γ). Thus, since [γ^1]φ∘=0, one has αφ=−ξ/γ0φ∈Lp(Γ)⊂M(Hs(Γ),H−s(Γ)).
By hypotheses (3.2), (3.3) and Theorem 2.5, we can introduce the self-adjoint operator in L2(R3) defined by
[TABLE]
The next theorem gives the connection between Aφ and the acoustic operator:
Theorem 3.7**.**
Let φ satisfy hypotheses (3.2) and (3.3), let
ψ∈dom(Aφ) and set u:=φ−1ψ. Then
[TABLE]
and
[TABLE]
Proof.
By the ”half” Green’s formula (see [28, Theorem 4.4], one gets
Since both Δφ and Δψ belong to H−1(R3) (notice that ψ∈dom(Aφ)⊆H1(R3)), the products ψΔφ and φΔψ are well-defined in D′(R3) and from(3.1) there follows
In this section the results provided in [27, Section 4], which in particular apply to A0,α (whenever α∈M(H23(Γ))), are extended to AV,α.
The weighted Sobolev spaces Hwk(R3) are
defined for k=0,1,2 and w∈R by
[TABLE]
[TABLE]
where ⟨x⟩ is a shorthand notation for the function x↦(1+∥x∥2)1/2. In particular, we set Lw2(R3)≡Hw0(R3). Since
[TABLE]
the two conditions ⟨x⟩wu∈L2(R3) and ⟨x⟩w∇u∈L2(R3) are equivalent to
⟨x⟩wu∈H1(R3); hence
[TABLE]
A similar argument applies to Hw2(R3)
[TABLE]
In particular, this provide the equivalent Hw2(R3)-norm
[TABLE]
The above definitions
are generalized to the case of non-integer order s∈R by
[TABLE]
while the corresponding dual spaces (w.r.t. the L2-product) identify with
[TABLE]
For the open subset Ω⊂R3, the spaces Hws(Ω) and Hws(R3\Ωˉ) are defined in a similar way. In particular,
since Ω is bounded, one has: Hws(Ω)=Hs(Ω), the
equalities holding in the Banach space sense; thus
[TABLE]
and
[TABLE]
The trace operators are extended to Hws(R3\Γ), w<0, by
[TABLE]
where χ∈Ccomp∞(Ωc), χ=1 on a neighborhood of Γ.
From now on we suppose that V∈Lcomp2(R3), so that σp(AV)∩(−∞,0)=∅ (see e.g. [22]) and, since V is a short range potential, a limiting absorption principle (LAP for short) holds for AV (see e.g. [2, Theorem 4.2]):
Theorem 4.1**.**
Let V∈Lcomp2(R3). For any k∈R\{0} and for any w>21, the limits
[TABLE]
exist in B(Lw2(R3),H−w2(R3)). Moreover
[TABLE]
and
[TABLE]
Remark 4.2**.**
By duality, the limits (4.4) also exist in B(Hw−2(R3),L−w2(R3)) and so, by interpolation,
[TABLE]
In order to extend LAP to operators of the kind AV,α, we need some preparatory lemmata. In the following BR denotes a sufficiently large ball such that supp(V)⊂BR.
Lemma 4.3**.**
Let V∈Lcomp2(R3). Then, for all z∈ρ(AV) and for all w∈R,
[TABLE]
Proof.
From the resolvent identity RzV=Rz0(1−VRzV), there follows
[TABLE]
Thus, since the thesis hold true in the case V=0 (this a a consequence of [35, Lemma 1, page 170], see the proof of Theorem 4.2 in [27]), we get
[TABLE]
Then the continuous injection H2(BR)↪L∞(BR) yields
[TABLE]
For w≥0, the embedding Lw2(R3)↪L2(R3) and the standard mapping
properties of RzV lead to
[TABLE]
and so, in this case, the statement follows from (4.7) and
(4.8). For w<0 we proceed as in the proof of [35, Lemma 1, page 170] starting from the identity
[TABLE]
An explicit computation leads to
[TABLE]
and so
[TABLE]
If ∣w∣∈[0,1], the functions
Δ⟨x⟩∣w∣ and
∇⟨x⟩∣w∣ are bounded
and smooth; then the functions VRzVΔ⟨x⟩∣w∣RzV
and VRzV∇⟨x⟩∣w∣⋅∇RzV
define bounded maps in B(L2(R3),Hσ2(R3)) for any σ∈R
(since V has compact support). In this case, we get
[TABLE]
and as before, we obtain (4.6) from (4.7). This
result and an induction argument on ∣w∣∈[n,n+1], allow to conclude the proof.
∎
Lemma 4.4**.**
Let V∈Lcomp2(R3) and w>1/2. Then, for all k2>0,
[TABLE]
Proof.
According to our assumptions, it results
[TABLE]
and the injection H2(BR)↪L∞(BR) yields
[TABLE]
Since R−k2V,±∈B(Lw2(R3),H−w2(R3)), the inequalities
[TABLE]
hold for w>1/2. Then the statement follows from (4.12) and
(4.13).
∎
This result yields the following mapping properties.
Lemma 4.5**.**
Let V∈Lcomp2(R3) and w>1. For all compact subsets K⊂(0,+∞) there exists cK>0 such that, for all k2∈K and for all u∈Lw2(R3)∩ker(R−k2V,+−R−k2V,−),
[TABLE]
Proof.
If V=0 the statement follows from [4, Corollary 5.7(b)]; in this
case for all k2∈K and u∈Lw2(R3)∩ker(R−k20,+−R−k20,−),
[TABLE]
for a suitable c~K>0 depending on K. From the identity
(4.5) there follows
[TABLE]
and
[TABLE]
Let u\in L_{{w}}^{2}(\mathbb{R}^{3})\cap\ker\big{(}R_{-k^{2}}^{V,+}-R_{-k^{2}}^{V,-}\big{)}; then
Hence, from the representation (4.16) and the estimates
(4.15), we finally obtain
[TABLE]
∎
The existence of the resolvent’s limits on the continuous spectrum has been
discussed in [36] for a wide class of operators including
singular perturbations. In the particular case of a singularly perturbed
Laplacian described through the general formalism introduced in [26],
a limiting absorption principle has
been given in [27]. In what follows, we use the same strategy used in
these works to establish a limiting absorption principle for the self-adjoint operators given in Theorem 2.5.
Theorem 4.6**.**
Let V∈Lcomp2(R3) and let AV,α defined as in Theorem 2.5. Then the limits
[TABLE]
exist in B(Lw2(R3),L−w2(R3)) for all w>1/2 and
k∈R\{0}.
Proof.
According to Theorem (4.1), the limits R−k2V,±
exists for all k2>0 and w>1/2 in the uniform operator topology of
B(Lw2(R3)),H−w2(R3)). Hence we follow, mutatis mutandis, the same arguments as in the proof on Theorem 4.1 in [27] (corresponding to the case V=0) to which we refer for more details: by [36, Theorem 3.5
and Proposition 4.2], our statement holds whenever there exist c1,
c2 and cK>0 (the last constant depending on K⊂(0,+∞)
compact), such that the following conditions are fulfilled:
[TABLE]
[TABLE]
(here B∞(L2(R3),Lσ2(R3)) denotes the space of compact
operators from L2(R3) to Lσ2(R3)), and, for all compact subset K⊂(0,+∞),
[TABLE]
Recalling that AV is bounded from above, there
exists c1>0 such that z∈ρ(AV) whenever Re(z)>c1; hence (4.23)
holds for RzV by (4.6). Since Γ is compact, by (4.6) and by the mapping properties of γ0, one has γ0RzV∈B(Lσ2(R3),H1(Γ)) and, by duality, SLzV∈B(H−1(Γ),L−σ2(R3)). Thus, formula (2.3)
gives (4.23) for RzV,α.
Since (1+αγ0SLzV)−1α∈B∞(Hs(Γ),H−s(Γ)), 0<s<1/2, by the compact embeddings H1(Γ)↪Hs(Γ) and H−s(Γ)↪H−1(Γ), one has (1+αγ0SLzV)−1α∈B∞(H1(Γ),H−1(Γ)). So, since γ0RzV∈B(L2(R3),H1(Γ)) and SLzV∈B(H−1(Γ),Lσ2(R3)), (4.24) follows from (2.3). Finally, the condition (4.25) holds as a consequence of the Lemma 4.5.
∎
The previous results also allow to prove that the resolvent formula (2.3) survives in the limits
z→−(k2±i0).
Theorem 4.7**.**
Let V∈Lcomp2(R3), k∈R\{0} and let AV,α defined as in Theorem 2.5.
For any w>21, the limits
[TABLE]
exist in B(H−s(Γ),H−w23−s(R3)), 0<s≤1, and
[TABLE]
[TABLE]
The function SL−k2V,±ξ solves, in the distribution space D′(R3\Γ) and for any ξ∈H−1(Γ), the equation
[TABLE]
and there exist ck2±>0 such that
[TABLE]
Moreover, the limits
[TABLE]
exist in B(Hs(Γ),H−s(Γ)), 0<s<1/2, and the operator 1+αγ0SL−k2V,± has a bounded inverse such that
[TABLE]
Finally, the limit resolvent R−k2V,α,± has the
representation
[TABLE]
Proof.
The proof uses exactly the same argumentation of the proofs of Lemma 4.4 and Theorem 4.5 (which give the analogous results in the case V=0) provided in [27] and so is left to the reader.
∎
5 Generalized eigenfunctions.
We say that a function u± which solves, outside some large ball BR, the Helmholtz equation (Δ+k2)u±=0, satisfies the (±)Sommerfeld radiation condition (or u±* is (±) radiating* for short) whenever
[TABLE]
holds uniformly in x^:=x/∣x∣.
Given ψk0=0, a generalized free eigenfunction with eigenvalue k2=0, i.e. ψk0∈Hloc2(R3) and (Δ+k2)ψk0=0, we say that ψkV,+/−=0 is an incoming/outgoing eigenfunction of −AV associated with the free waveψk0 whenever ψkV,±∈Hloc2(R3) solves (A~V+k2)ψkV,±=0 and the scattered fieldψk,scV,±:=ψkV,±−ψk0 satisfies the (±) Sommerfeld radiation condition. Here
A~V:Hloc2(R3)⊂Lloc2(R3)→Lloc2(R3), V∈Lcomp2(R), denotes the broadening of AV defined by A~Vψ:=Δϕ−Vψ.
Let us notice that ψk,scV,± satisfies the Helmholtz equation outside the support of V.
The next result is a consequence of LAP for AV:
Theorem 5.1**.**
The unique incoming and outgoing eigenfunctions of −AV, V∈Lcomp2(R3), associated with the free wave ψk0, k=0, are given by
[TABLE]
Proof.
By definition, ψk+/−∈Hloc2(R3) is an incoming/outgoing eigenfunction of −AV associated with ψk0 if and only if (ψk±−ψk0) is a (±) radiating solution of (A~V+k2)u=Vψk0. Since the potential V is compactly supported, such an equation has an unique (±) radiating solution. Indeed, if u1 and u2 were two different solutions then u:=u1−u2 would be a radiating solution, outside some large ball BR containing the support of V, of (Δ+k2)u=0. Thus u∣BRc=0. Then, by the unique continuation principle for A~V (see [22]), one gets u=0 everywhere. By Theorem 4.1,
ψk,scV,±:=−R−k2V,±Vψk0∈H−w2(R3) solve the equation (A~V+k2)ψk,scV,±=Vψk0. Moreover, by Theorem 4.1,
ψk,scV,±=R−k20,±V(1−R−k2V,±V)ψk0. Since V is compactly supported, ψk,scV,± is (±) radiating by [12, Lemma 7, Subsection 7d, Section 8, Chapter II].
∎
Now we extend the previous result to AV,α. At first we introduce the following broadening of AV,α to the larger space Lloc2(R3):
[TABLE]
[TABLE]
[TABLE]
We say that ψkV,α,+/−=0 is an incoming/outgoing eigenfunction of −AV,α associated with the free wave ψk0, whenever ψkV,α±∈dom(A~V,α) solves the equation (A~V,α+k2)ψkV,α,±=0 and the scattered fieldψk,scV,α±:=ψkV,α±−ψkV,± is (±) radiating, where ψkV,+/− is the unique incoming/outgoing eigenfunction of −AV associated, according to Theorem 5.1, with ψk0. Let us notice that ψk,scV,α,± satisfies the Helmholtz equation outside the supp(V)∪Γ.
Theorem 5.2**.**
Suppose that R3\Ω is connected and V∈Lcomp2(R3); let −AV,α be defined as in Theorem 2.5. Then the unique incoming and outgoing eigenfunctions of −AV,α associated with the free wave ψk0 are given by
[TABLE]
Proof.
By our definitions, ψ~k+/−∈dom(A~V,α) is an incoming/outgoing eigenfunction of −AV,α associated with ψk0 if and only if (ψ~k±−ψkV,±) is a (±) radiating solution of (AV+k2)u−γ0∗αγ0u=γ0∗αγ0ψkV,± belonging to Hloc23−s(R3). Since both the potential V and the distribution γ0∗ξ are compactly supported, such an equation as an unique (±) radiating solution. Indeed, if u1 and u2 were two different solutions then u:=u1−u2 would be a radiating solution, outside some large ball BR containing both supp(V) and supp(γ0∗ξ), of (Δ+k2)u=0. Thus u∣BRc=0. Then, by the unique continuation principle for AV (see [22]), one gets u∣R3\Ω=0. Since u∈Hloc23−s(R3), then γ0u=0 and so u is a radiating solution of (AV+k2)u=0; thus, proceeding as in the proof of Theorem 5.1, u=0 everywhere. To conclude the proof we need to show that ψkV,α,±∈dom(A~V,α), i.e. that ψ∘:=ψkV,α,±+SL∘Vαγ0ψkV,α,±∈Hloc2(R3), that (A~V,α+k2)ψkV,α,±=0 and that SL−k2V,±(1+αγ0SL−k2V,±)−1αγ0ψkV,± in (±) radiating. Since αγ0ψkV,α,±=(1+αγ0SL−k2V,±)−1αγ0ψkV,±, one has, by (4.27),
[TABLE]
Then
[TABLE]
Finally, by (4.5) and by (4.28), SL−k2V,±ξ=R−k20,±(γ0∗ξ−VSL−k2V,±ξ) and so, since both γ0∗ξ and V are compactly supported, SL−k2V,±ξ is (±) radiating by [12, Lemma 7, Subsection 7d, Section 8, Chapter II].
∎
Remark 5.3**.**
By the resolvent identity RzV=Rz0−RzVVRz0 and by (4.27), one gets
SL−k2V,±ξ=SLz0ξ+ϕξ, where ϕξ∈H−w2(R3). Thus, by [γ^1]SLz0ξ=−ξ (see [28, Theorem 6.11]) and H−w2(R3)⊂ker[γ^1], one obtains
[TABLE]
Then, by the identity αγ0ψkV,α,±=(1+αγ0SL−k2V,±)−1αγ0ψkV,±, one obtains the relations
[TABLE]
For any k>0, we define the set
[TABLE]
where ⋅ denotes the euclidean scalar product, equivalently
[TABLE]
Clearly any function of the kind ψρ(x):=eρ⋅x, ρ∈Σk, is a generalized eigenfunction of −Δ with eigenvalue k2.
Corollary 5.4**.**
Given a ball BR∘⊃Ω∪supp(V), the outgoing eigenfunction ψρV,α associated, according to Theorem 5.2, with ψρ(x):=eρ⋅x, ρ∈Σk,
has the asymptotic behavior
[TABLE]
uniformly in all directions x^:=∣x∣x. Moreover
[TABLE]
where
[TABLE]
Proof.
By Theorems 5.1 and 5.2, and by the identity αγ0ψρV,α=(1+αγ0SL−k2V,−)−1αγ0ψρV, where ψρV denotes the outgoing eigenfunction associated, according to Theorem 5.1, with ψρ, one has ψρV,α=ψρ+ψρ,scV,α. Since (−Δ+k2)ψρ,scV,α=0 outside Ω∪supp(V), the thesis is consequence of the asymptotic representation of the radiating solutions of the Helmholtz equation (see e.g. [11, Theorem 2.6]).
∎
Remark 5.5**.**
Since (−Δ+k2)ψρ,scV,α=0 outside Ω∪supp(V), by elliptic regularity ψρ,scV,α is smooth outside Ω∪supp(V) and so relation (5.2) is well defined.
According to Corollary 5.4, the scattering amplitude sV,α for the Schrödinger operator AV,α is then related to the far-field patternψV,α∞ by the simple relation
[TABLE]
Indeed, by Corollary 5.4, the outgoing eigenfunction ψkV,α associated, according to Theorem 5.2, with ψk0(x):=eikξ^⋅x,
has the asymptotic behavior
[TABLE]
The next lemma shows that the scattering amplitude univocally determines both the far-field ψV,α∞ and the scattered field ψρ,scV,α.
Remark 5.6**.**
Here and below, when considering two different self-adjoint operators AV1,α1 and AV2,α2 we mean that they can be eventually be defined in terms of two different subset Ω1 and Ω2, so that (∂Ω1=)Γ1=Γ2(=∂Ω2) is allowed.
Lemma 5.7**.**
Under the same hypotheses as inTheorem 5.2, suppose that, for some k>0,
[TABLE]
Then
[TABLE]
and
[TABLE]
where BR∘⊃(Ω1∪Ω2∪supp(V1)∪supp(V2)).
Proof.
By (5.3) and (5.2), to get (5.4) it suffices to show that, for some R>R∘, if ψikξ^,scV1,α1∣BR=ψikξ^,scV2,α2∣BR, for all ξ^∈S2 then ψρ,scV1,α1∣BR=ψρ,scV2,α2∣BR for all ρ∈Σk.
Since C(S2) is dense in L2(S2), according to [39, Theorem 2] there exists a sequence {fn}1∞⊂C(S2)
such that Hkfn→ψρ in H2(BR), where Hk, k=0, is the Herglotz operator
[TABLE]
Writing the above integral as a limit of a Riemann’s sum, ψρ can be obtained as a H2(BR)-limit of a sequence of functions of the kind ∑m=1nam,neikξ^m,n⋅x. Since, by Theorems 5.1 and 5.2,
[TABLE]
to get (5.4) one needs to show that the linear operator L−k2V,α is continuous on Hloc2(R3) to L−w2(R3). Since V∈Lcomp2(R3), the multiplication operator associated with V belongs to B(H2(BR),Lw2(R3)); thus, by Theorem 4.1, R−k2V,−V∈B(H2(BR),H−w2(R3)). Moreover, by Theorem 4.7, SL−k2V,−(1+αγ0SL−k2V,−)−1αγ0∈B(H21+s(BR),H−w23−s(R3)). So L−k2V,α∈B(H2(BR),H−w23−s(R3)) and (5.4) holds true.
If ψV1,α1∞(k,ρ,x^)=ψV2,α2∞(k,ρ,x^), then, by Corollary 5.4, usc(x):=ψρ,scV1,α1(x)−ψρ,scV2,α2(x)=O(∣x∣−2). Since usc solves the Helmoltz equation (Δ+k2)usc=0 outside BR∘, by Rellich’s lemma (see e.g. [11, Theorem 2.14]), one gets usc∣BR∘=0, i.e. (5.5).
∎
6 Uniqueness in inverse Schrödinger scattering.
Given V∈Lcomp2(R3) and α∈H−s(Γ), 0<s≤1, let us define
V~α∈Hcomp−s−21(R3), by V~α:=V+γ0∗α.
For any ψ∈Hlocs+21(R3) the product ψV~α∈D′(R3) is well defined and, by Lemma 3.4, ψV~α=Vψ+γ0∗αγ0ψ. Thus, since M(Hs(Γ),H−s(Γ))⊆H−s(Γ), we have
[TABLE]
and so ψV~α∈Hcomp−s−21(R3) whenever α∈M(Hs(Γ),H−s(Γ)) and ψ∈Hlocs+21(R3). In particular, whenever V and α are as in the definition of A~V,α and ψ∈dom(A~V,α), one has ψV~α∈Hcomp−1(R3).
We need a preparatory lemma before stating the main result in this section:
Lemma 6.1**.**
Let ψρ1V1,α1 and ψρ2V2,α2 be the outgoing eigenfunctions of AV1,α1 and AV2,α2, associated, according to Theorem 5.2, with ψρ1(x)=eρ1⋅x and ψρ2(x)=eρ2⋅x, where both ρ1 and ρ2 belong to Σk. Then
[TABLE]
Proof.
By the definition of V~α, by Lemma 3.4 and by Remark 5.3, one obtains
[TABLE]
where γ0,m and [γ^1,m] denote the trace operators on Γm and
the jump of the normal derivatives across Γm respectively.
Let ψm,ρ, m=1,2, be outgoing eigenfunctions of A~Vm,αm associated with ψρ, ρ∈Σk. Setting
[TABLE]
where BR⊃(Ω1∪Ω2∪supp(V1)∪supp(V2)),
one gets
[TABLE]
so that
[TABLE]
Then, according to the half Green’s formula (see [28, Theorem 4.4]), one obtains
[TABLE]
where γ0,R and γ1,R denote the trace operator and the normal derivative on ∂BR respectively. Thus, since Δψm,ρ=(V~m−k2)ψm,ρ, by (LABEL:Vt), one gets
Finally we state our uniqueness result for inverse Schrödinger scattering:
Theorem 6.2**.**
Let V1,V2∈Lcomp2(R3), α1,∈M(Hs(Γ1),H−s(Γ1)), α2∈M(Hs(Γ2),H−s(Γ2)), 0<s<1/2, and suppose that R3\Ω1 and R3\Ω2 are connected. Then
[TABLE]
Proof.
Let ψρmVm,αm, m=1,2, be as in Lemma 6.1 and choose ρ1 and ρ2 in Σk in such a way that ρˉ1+ρ2=−iξ, ξ∈R3. Further set ϕρmVm,αm(x):=e−ρ⋅xψρmVm,αm(x)−1, so that ψρmVm,αm(x)=eρ⋅x(1+ϕρmVm,αm(x)). Then, by Lemma 6.1, setting uξ(x):=e−iξ⋅x,
[TABLE]
By our definitions, setting ψρ(x)=eρ⋅x, one obtains
[TABLE]
Thus ϕρmV,α solves the equation
[TABLE]
Decaying estimates of solutions of such an equation have been obtained in various papers concerning Calderón’s uniqueness problem. In particular we use the recent results provided in [18].
Since αm∈M(Hs(Γ),H−s(Γ))⊆H−s(Γ), 0<s<1/2, one has V~αm∈Hcomp−1(R3)⊂Wcomp−1,3/2(R3).
Thus (see e.g. [1, Theorem 3.12 and Corollary 3.23]) V~αm=∑i=13∇ifm,i+hm where fm,i,h∈L3(R3). Then, taking χm∈Ccomp∞(R3) such that χm=1 on supp(V~αm), one has V~αm=χmV~αm=∑i=13∇i(χmfm,i)−∑i=13fm,i∇iχm+χmhm=∑i=13∇if~m,i+h~m where f~m,i,h~∈Lcomp3(R3). Therefore [18, Theorem 5.3]111Such a theorem is stated for q=γ−1/2Δγ1/2, ∇logγ∈Lcomp3(R3); however the proof only uses the decomposition q=∑i=13∇ifi+h where fi∈Lcomp3(R3) and h∈L3/2(R3). applies to V~αm and so, by the same reasoning as in the (second part) of the proof of Theorem 1.1 in [18] (see in particular inequality (31)), for any ξ∈R3 one gets the existence of two suitable sequences {ρm,n}n=1+∞⊂Σk, ρ1,n+ρ2,n=−iξ, such that
[TABLE]
This implies V~α1=V~α2 and so V~α1=V~α2. Then, by the definition of V~αm, one obtains V1−V2=γ0,2∗α2−γ0,1∗α1, where γ0,1 and γ0,2 denote the trace operators on Γ1 and Γ2 respectively.
This entails V1=V2, supp(α1)=supp(α2) and α1=α2.
∎
7 Uniqueness in inverse acoustic scattering.
The next lemma probably contains well-known results but we found no proof in the literature.
Lemma 7.1**.**
Let ϱ≥0 satisfy the hypotheses
[TABLE]
Then
[TABLE]
Let us further suppose that ϱ is constant outside some bounded ball BR∘ and there exists an open and bounded set Ωϱ≡Ω⊂BR∘ with Lipschitz
boundary Γϱ≡Γ such that
[TABLE]
Then
[TABLE]
and
[TABLE]
where γ0ϱ and [γ^1]ϱ denote the trace on Γ and the jump of the normal derivative across Γ respectively.
Proof.
At first we define the sequence ϱn:=eΔ/nϱ, n≥1. Then, since the heat semigroup is positivity-preserving, strongly continuous in both L∞(R3) and L2(R3), commutes with ∇ and etΔ(L∞(R3))⊂C∞(R3) whenever t>0, one gets ϱn≥0, ϱn∈C∞(R3)∩L∞(R3), ∇ϱn∈L2(R3;R3), ϱn→ϱ in L∞(R3) and ∇ϱn→∇ϱ in L2(R3;R3) as n→+∞. Since ϱ(x)≥∥1/ϱ∥L∞(R3)−1 for a.e. x∈R3, ϱn(x) is definitively stricitly positive uniformly in x∈R3 and so ϱn−3/2→ϱ−3/2 in L∞(R3). Thus, by ∇ϱn=−2−1ϱn−3/2∇ϱn→−2−1ϱ−3/2∇ϱ in L2(R3), (7.2) follows.
By Lemma 7.1, one immediately gets χϱ−1/2∣Ωin/ex∈HΔ1(Ωin/ex), χ∈Ccomp∞(R3). Then, by the half Green’s formula (see [28, Theorem 4.4]) one has
[TABLE]
for any v∈H1(R3). Since both ϱ−1 and ϱ−1/2 belong to Hloc1(R3)∩L∞(R3), one has ϱ−3/2∈Hloc1(R3) and so we can use the above Green’s formula in the case v=ϱ−3/2w, w∈Ccomp∞(R3). Then, by (7.6) and (7.2), one obtains
[TABLE]
Therefore [γ^1](ϱ−1/2)=−2−1γ0(ϱ−3/2)[γ^1]ϱ. By Lemma 3.3, γ0(ϱ−3/2)=(γ0ϱ)−3/2 and the proof is concluded.
∎
Now we introduce a further hypothesis on ϱ :
[TABLE]
In particular, by Remark 2.6, hypothesis (7.7) holds true whenever
[TABLE]
Corollary 7.2**.**
If ϱ≥0 satisfies(7.1), (7.3) and (7.7), then φ:=ϱ−1/2 satisfies (3.2) and (3.3).
Proof.
Hypotheses (3.2) are consequence of Lemma 7.1. By Lemma 3.3, one has γ0(ϱ−1/2)=(γ0ϱ)−1/2 and then hypothesis (3.3) follows from (7.5).
∎
Let us now take v≥0 such that v−1∈L∞(R3), suppose that ϱ(x)=v(x)=1 whenever x lies outside some large ball BR∘ and set
[TABLE]
[TABLE]
so that Vφ∈Lcomp2(R3) and Vv,ω∈Lcomp∞(R3).
By Lemma 2.11 and Theorem 3.7, there is a well defined correspondence between the outgoing eigenfunctions ψωφ,v,ω of Aφ+Vv,ω≡AVφ,v,ω,αφ provided in Theorem 5.2 and the acoustic eigenfunctions uωϱ,v such that
[TABLE]
Such a relation is given by uωϱ,v=ϱψωφ,v,ω. Notice that u(t,x):=e−iωtuωϱ,v(x) is a fixed-frequency solution of the acoustic wave equation
[TABLE]
By Corollary 5.4, since uωϱ,v=ψωφ,v,ω outside BR∘, the eigenfunction uωϱ,v has the asymptotic behavior
[TABLE]
uniformly in all directions x^:=∣x∣x, where the far-field pattern uϱ,v∞
is related to the scattering amplitude for the Schrödinger operator AVφ,v,ω,αφ by
[TABLE]
By Theorem 5.2, uωϱ,v is the unique solution of the stationary acoustic equation (7.8) such that the scattered field uω,scϱ,v(x):=uωϱ,v(x)−e−iωξ^⋅x satisfies the outgoing Sommerfeld radiation condition.
Remark 7.3**.**
An more explicit characterization of a class of function φ satisfying hypotheses (7.1), (7.3) and (7.7) is the following:
[TABLE]
where χ∘∈Ccomp∞(R3), χ∘=1 on some large ball containing Ω, ϱ∘∈H2(R3) and ξ∈Lp(Γ), p>8/3. This includes the case where the normal derivative of ϱ has a jump across Γ which is locally supported on a closed subset Σ⊂Γ. By the same reasoning as in Remark 3.6 one has ∣∇ϱ∣∈L2(R3), ΔΩin/exϱ∈L2(Ωin/ex) and [γ^1]ϱ/γ0ϱ∈Lp(Γ)⊂M(Hs(Γ),H−s(Γ)). Moreover, by H2(R3)⊂W1,6(R3), by
SLξ∈W1+1/p−ϵ,p(Ωin/ex) (see [14, Theorem 3.1]) and by
W1+1/p−ϵ,p(Ωin/ex)⊂W1,4(Ωin/ex) whenever p>8/3,
one has ∣∇Ω\/exϱ∣∈L4(Ωin/ex).
Thanks to Theorem 6.2, one gets the following uniqueness result in acoustic scattering:
Theorem 7.4**.**
Let ϱ1≥0, ϱ2≥0 satisfy hypotheses (7.1), (7.3), (7.7) (see for example Remark 7.3) and let v1≥0, v2≥0 such that v1−1,v2−1∈L∞(R3). Suppose R3\Ωϱ1 and R3\Ωϱ2 are connected and that
ϱ1(x)=ϱ2(x)=v1(x)=v1(x)=1 whenever x lies outside some large ball BR∘⊃(Ωϱ1∪Ωϱ2). Then, given ω=ω~=0,
[TABLE]
Proof.
Let φm=ϱm−1/2, m=1,2. By (7.9) and Theorem 6.2, one gets
[TABLE]
Furthermore, considering any open and bounded Ω with Lipschitz boundary Γ such that supp(αφ1)=supp(αφ2)⊆Γ one has, in both L2(Ωin) and L2(Ωex),
[TABLE]
[TABLE]
Thus, since ω=ω~=0, one has v1=v2 and
[TABLE]
Let us set uin/ex:=(φ2−φ1)∣Ωin/ex, so that
[TABLE]
Since uex=0 outside BR∘ and Vex∈Lcomp2(Ωex), the unique continuation principle (see e.g. [22]) leads us to φ1=φ2 on Ωex. In addition, due to Corollary 7.2, both φ1 and φ2 belong to H1(R3); then the previous identity yields γ0φ1=γ0φ2 and γ^1exφ1=γ^1exφ2. Hence, setting u:=uin⊕uex, one has [γ0]u=0 and, by (7.10), [γ^1]u=0. By elliptic regularity, uin∈H2(Ωin); so u belongs to H2(R3) and solves (−Δ+Vin⊕Vex)u=0. Using again the unique continuation principle, this entails u=0. Therefore φ1=φ2, i.e. ϱ1=ϱ2.
∎
[2]
S. Agmon. Spectral properties of Schrödinger operators and
Scattering Theory. Ann. Scuola Sup. Pisa (IV) 11 (1975), 151-218.
[3]
W.O. Amrein, J.M. Jauch, K.B. Sinha. Scattering theory
in quantum mechanics. Physical principles and mathematical methods. W. A.
Benjamin, Reading, Masschusetts, 1977.
[4]
M. Ben-Artzi, A. Devinatz. The limiting absorption principle
for partial differential operators. Memoirs of the AMS vol. 364, 1987.
[5]
J. Behrndt, M. Langer, V. Lotoreichik. Schrödinger operators
with δ- and δ′-potentials supported on hypersurfaces.
Ann. Henri Poincaré14 (2013), 385-423.
[6]
R. Brown. Global uniqueness in the impedance-imaging problem
for less regular conductivities. SIAM J. Math. Anal.27
(1996), 1049-1056.
[7]
R. Brown, R. Torres. Uniqueness in the inverse
conductivity problem for conductivities with 3/2 derivatives in Lp,
p>2n. J. Fourier Anal. Appl.9 (2003), 563-574 .
[8]
F. Cakoni, D. Colton. Qualitative methods in inverse
scattering theory. An introduction. Interaction of Mechanics and
Mathematics. Springer-Verlag, Berlin, 2006.
[9]
A. Calderón. On an inverse boundary value problem, in:
Seminar on Numerical Analysis and its Applications to Continuum
Physics (Rio de Janeiro, 1980), 65-73, Soc. Brasil. Mat., Rio de Janeiro,
[10]
P. Caro, K. Rogers. Global uniqueness for the Calderón
problem with Lipschitz conductivities. Forum Math.4 (2016),
28 pages.
[11]
D. Colton, R. Kress. Inverse acoustic and electromagnetic
scattering theory, 3rd edition. Springer-Verlag, New-York, 2013.
[12]
R. Dautray, J-L. Lions. Mathematical Analysis and Numerical
Methods for Science and Technology. Vol. 1: Physical Origins and Classical
Methods. Springer, Berlin, 1990.
[13]
P.A. Deift. Applications of a commutation formula.
Duke Math. J.45 (1978), 267-310.
[14]
E. Fabes, O. Mendez, M. Mitrea. Boundary layers on Sobolev-Besov
spaces and Poisson’s equation for the Laplacian in Lipschitz domains.
J. Funct. Anal., 159, 323-368, 1998.
[15]
A. Greenleaf, M. Lassas, G. Uhlmann. The Calderón problem
for conormal potentials. I. Global uniqueness and reconstruction.
Comm. Pure Appl. Math.56 (2003), 328-352.
[16]
P. Grisvard. Elliptic problems in nonsmooth domains.
Pitman, Boston, 1985.
[17]
B. Haberman, D. Tataru.: Uniqueness in Calderón
problem with Lipschitz conductivities. Duke Math. J.162
(2013), 497-516.
[18]
B. Haberman. Uniqueness in Calderón’s problem for
conductivities with unbounded gradient. Comm. Math. Phys.340, 639-659, 2015.
[19]
V.G. Maz’ya, T.O. Shaposhnikova. Theory of Sobolev
Multipliers. Springer-Verlag, Berlin, 2009.
[20]
Hähner, Peter. On the uniqueness of the shape of a
penetrable, anisotropic obstacle. J. Comput. Appl. Math.116
(2000), 167-180.
[21]
V. Isakov. Inverse problems for partial
differential equations. Second edition. Springer, New York, 2006.
[22]
D. Jerison, C.E. Kenig. Unique continuation and absence of
positive eigenvalues for Schrödinger operators. With an appendix by E. M.
Stein. Ann. of Math., 121(2), 463-494, 1985.
[23]
T. Kato. Perturbation Theory for Linear Operators.
Springer, Berlin 1976.
[24]
A. Kirsch; A. Lechleiter. The operator equations of
Lippmann-Schwinger type for acoustic and electromagnetic scattering problems
in L2. Appl. Anal.88 (2009), 807-830.
[25]
R. Leis. Initial-boundary value problems in mathematical
physics. B.G. Teubner, Stuttgart, 1986.
[26]
A. Mantile, A. Posilicano. M. Sini. Self-adjoint elliptic
operators with boundary conditions on not closed hypersurfaces. J.
Differential Equations, 261, 1-55, 2016.
[27]
A. Mantile, A. Posilicano. M. Sini. Limiting absorption
principle, generalized eigenfunction and scattering matrix for Laplace
operators with boundary conditions on hypersurfaces. arXiv:1605.03240v3 (to appear in J. Spectr.
Theory ).
[28]
W. McLean. Strongly elliptic systems and boundary
integral equations. Cambridge University Press, Cambrdige, 2000.
[29]
R.G. Novikov.
A multidimensional inverse
spectral problem for the equation −Δψ+(v(x)−Eu(x))ψ=0.
Funct. Anal. Its Appl.22 (1988), 263-272.
[30]
L. Päivärinta, A. Panchenko, G. Uhlmann.: Complex
geometrical optics solutions for Lipschitz conductivities. Rev. Mat.
Iberoamericana, 19 (2003), 57-72.
[31]
A. Posilicano. A Krĕın-like formula for singular
perturbations of self-adjoint operators and applications. J. Funct.
Anal., 183, 109-147, 2001.
[32]
A.G. Ramm.
Recovery of the potential from
fixed-energy scattering data.
Inverse Problems, 4
(1988), 877-886, .
[33]
A. G. Ramm. Inverse problems. Mathematical and
analytical techniques with applications to engineering. Springer, New York, 2005.
[34]
M. Reed, B. Simon. The scattering of classical waves from
inhomogeneous media. Math. Z.155 (1977), 163-180.
[35]
M. Reed, B. Simon. *Methods of modern mathematical physics
vol. IV: Analysis of Operators. *Academic Press, New YorK, 1978.
[36]
W. Renger. Limiting absorption principle for singularly
perturbed operators. Math. Nachr., 228, 163-187, 2001.
[37]
J. Sylvester, J.G. Uhlmann: A global uniqueness theorem for an
inverse boundary value problem. Ann. Math.125 (1987), 153-169.
[38]
G. Uhlmann: Inverse problems: seeing the unseen.
Bull. Math. Sci.4 (2014), 209-279.
[39]
N. Weck. Approximation by Herglotz wave functions.
Math. Methods Appl. Sci., 27, 155-162, 2004.
Bibliography39
The reference list from the paper itself. Each links out to its DOI / PubMed record.
2[2] S. Agmon. Spectral properties of Schrödinger operators and Scattering Theory. Ann. Scuola Sup. Pisa (IV) 11 (1975), 151-218.
3[3] W.O. Amrein, J.M. Jauch, K.B. Sinha. Scattering theory in quantum mechanics. Physical principles and mathematical methods. W. A. Benjamin, Reading, Masschusetts, 1977.
4[4] M. Ben-Artzi, A. Devinatz. The limiting absorption principle for partial differential operators. Memoirs of the AMS vol. 364, 1987.
5[5] 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.
6[6] R. Brown. Global uniqueness in the impedance-imaging problem for less regular conductivities. SIAM J. Math. Anal. 27 (1996), 1049-1056.
7[7] R. Brown, R. Torres. Uniqueness in the inverse conductivity problem for conductivities with 3 / 2 3 2 3/2 derivatives in L p superscript 𝐿 𝑝 L^{p} , p > 2 n 𝑝 2 𝑛 p>2n . J. Fourier Anal. Appl. 9 (2003), 563-574 .
8[8] F. Cakoni, D. Colton. Qualitative methods in inverse scattering theory. An introduction. Interaction of Mechanics and Mathematics. Springer-Verlag, Berlin, 2006.