arxiv: v2 [math.ap] 15 Jul 2016

Transcription

1 RECONSTRUCTION AND STABILITY FOR PIECEWISE SMOOTH POTENTIALS IN THE PLANE arxiv: v [math.ap] 15 Jul 016 JORGE TEJERO Abstract. We show that complex-valued potentials with jump discontinuities can be recovered from the Dirichlet-to-Neumann map using Bukhgeim s method. Combining with known formulas, this enables the recovery from the scattering amplitude at a fixed energy. We also provide a priori stability estimates for reconstruction from the Dirichlet -to-neumann map as well as from the scattering amplitude given an approximate knowledge of the location of the discontinuities. 1. Introduction Let Ω be a bounded planar domain that contains a bounded potential V, and consider the Dirichlet problem { u = Vu (1) u Ω = f. Supposing that 0 is not a Dirichlet eigenvalue, there is a unique solution u H 1 (Ω) and the DtN map Λ V can be formally defined by Λ V : f u n Ω. Our goal is then to recover the potential from the information contained in Λ V. This problem has a long history and is closely related to the inverse conductivity problem proposed by Calderón; see [14]. Some relevant work in higher dimensions includes [34, 7, 9, 11, 5,, 1, 16]. The two dimensional question is quite different to the higher dimensional case; for example the inverse problem is no longer overdetermined. In [8] Nachman introduced the -method to prove uniqueness for conductivities in W,p with p > 1, and gave a reconstruction procedure (this work has since been extended to more general cases, see for example [3, 1]). In [7], combining the -method with the theory of quasi-conformal maps, Astala and Päivärinta solved the uniqueness problem for L conductivities. There was little progress for general potentials in the plane until Bukhgeim introduced a new method, proving uniqueness for C 1 -potentials [13]. There he took advantage of solutions of the form u = u,x = e iψ (1+w), ψ(z) = ψ x (z) = 1 ( z1 x 1 +i(z x ) ), 1

2 JORGE TEJERO where w = w,x is small in some sense. Indeed, after proving that the the integral e i(ψ+ψ) V w Ω converges to zero as grows, one uses the stationary phase method to see that the integral e iψ (Λ V Λ 0 )[u] = e i(ψ+ψ) V (1+w) Ω converges to π 1 V(x). This quadratic phase approach has since been extended by many authors; we list a few here. Novikov and Santacesaria obtained stability estimates for C -potentials in [30] andproved that thereconstruction method could be extendedto C 1 matrix-valued potentials in [31]. In[10], Blåsten, Imanuvilov and Yamamoto proved uniqueness for potentials in L p with p >, and gave stability estimates in the L norm for potentials in H s with s > 0. Finally, Astala, Faraco and Rogers gave a reconstruction procedure for potentials in H 1/ and proved that this is best possible in some sense; see [4]. In the recent work [6], Lakshtanov, Novikov and Vainberg provided a reconstruction scheme for real bounded potentials. Their procedure relies on Faddeev s scattering solutions and allows to recover almost each potential, in the sense of their Remark 4.1. Here we will prove that the reconstruction formula given in [4] (mildly different from the original Bukhgeim formula) works for piecewise W,1 - potentials with jump discontinuities on smooth curves. That is to say potentials V that can be written as N V(x) = q j (x)χ Ωj (x), j=1 where q j W,1 (Ω j ) and Ω j are piecewise C,α -domains with α > 1/. By this we mean, the boundary Ω j can be expressed as a finite union of graphs of C,α -functions and that the union is Lipschitz. The most significant novelties of the article are to be found in Section 3, where we will prove the following theorem. There we will also present a potential for which the recovery formula fails at points away from the discontinuities. Theorem 1.1. Let V be a piecewise W,1 -potential with jump discontinuities on smooth curves. Then lim e iψ (Λ V Λ 0 )[u,x ] = V(x), a.e. x Ω. π Ω We refer to Theorem 1.1 of [4] for how to determine the values of the Bukhgeim solutions on the boundary. For a fixed k > 0, we also consider the Schrödinger equation () ( +V)u = k u, Ω

3 PIECEWISE SMOOTH RECOVERY 3 where k is not a Dirichlet eigenvalue for the Hamiltonian + V. For θ S 1, the outgoing scattering solutions satisfy the Lippmann-Schwinger equation u(x,θ) = e ikx θ G 0 (x,y)v(y)u(y,θ)dy, R where G 0 denotes the outgoing Green s function which satisfies ( k )G 0 (x,y) = δ(x y). Then the scattering amplitude A V : S 1 S 1 C at energy k can bewritten A V (η,θ) = e ikη y V(y)u(y,θ)dy. R Given an incident plane wave in direction θ, the scattering amplitude measures the probability of scattering in the direction η. A classical problem is to recover the potential in () from the information contained in A V. By applyingtheformulaobtained in [5], wecan also recover V from A V as long as V is a piecewise W,1 -potential with jump discontinuities on smooth curves. More details will be given in the third section. Stability estimates are a classical theme in inverse problems; see [3, 10] for the Schrödinger equation and [, 8, 9, 17, 15] for the conductivity equation. Notice that in [10] or in [17] there is stability for discontinuous potentials or conductivities butonly in the L sense. A careful analysis of the dependence of the constants in the reconstruction theorem yields an L stability for discontinuous coefficients provided a noise knowledge of the discontinuities of the potential. In the final section we will prove stability estimates for the reconstruction from the DtN map and from the information contained in A V.. Preliminaries.1. Quadratic phase solutions and integrals. Using Wirtinger derivativeswecanwritethethetime-independentschrödingerequationas4 z z u = Vu. Taking solutions of the form u = e iψ (1+w) and multiplying both sides by e iψ we obtain 4e iψ z z e iψ (1+w) = e i(ψ+ψ) V(1+w). Taking into account that z e iψ = z e iψ = 0, this can be rewritten as 4 z e i(ψ+ψ) z w = e i(ψ+ψ) V(1+w). As the derivatives are local operators, and we need only satisfying the equation inside Ω, we can take solutions of the form w = 1 [ [ ]] 4 1 z e i(ψ+ψ) χ Q z 1 e i(ψ+ψ) χ Q V (1+w)

4 4 JORGE TEJERO where Q is an auxiliary axis-parallel square containing Ω. In order to simplify notation we define the multiplication operators and write M ± [F] = e ±i(ψ+ψ) χ Q F, S 1 [F] = z M 1 z M, S V [F] = S 1 [VF]. ForsufficientlysmoothV,theoperatornormofS V issmallforlargeenough; see [4, Lemma.3], so we can invert (I SV ) using Neumann series, yielding w = w,x = (I S V ) 1 S 1 [V]. The following oscillatory integral operator will also be useful in the sequel Tw [F](x) = e i(ψ+ψ) F(z)w,x (z)dz. π R In order to study the behaviour of these operators we will use the homogeneous L Sobolev space, denoted by Ḣs, with norm f Ḣs = s ˆf L, where ˆf is the Fourier transform of f. We have the following bound for M ± ; the proof can be found in [4, Section ]. The key ingredient in the proof is the classical lemma of van der Corput. Lemma.1. Let 0 s 1,s < 1. Then M ± [F] Ḣ s C min{s 1,s } F Ḣs 1, 1. The following two lemmas were essentially proven in [4, Sections and 4]; we present minor modifications, suitable for the stability analysis of the final section. Lemma.. Let 0 < s 1,s < 1. Then there exists a constant C such that C τ where τ = 1 s +min{s 1,s }. S1 Ḣs 1 Ḣ s Proof. Using Lemma.1 twice we get S1 Ḣs M 1 1 Ḣ s z M Ḣs 1 Ḣ s 1 C 1+s z 1 M Ḣs 1 Ḣ 1 s C 1+s M Ḣs 1 Ḣ s and the proof is concluded. C τ,

5 PIECEWISE SMOOTH RECOVERY 5 Lemma.3. Let F,V Ḣs where 0 < s < 1. Then there exists a constant C such that suptw[f](x) C s F Ḣs V Ḣs x Ω when is sufficiently large. Proof. Using Lemma.1 we obtain Tw[F](x) CM [F] Ḣ s w Ḣ s C 1 s F Ḣs (I SV) 1 S1[V] Ḣs. As (I SV ) 1 is bounded for sufficiently large (see [4, Lemma.3]), Tw[F](x) C 1 s F Ḣs S1[V] Ḣs C s F Ḣs V Ḣs, where the last inequality comes from applying Lemma.. Let g : R R C n. We say that x s is a stationary point of g of order m < n if g (k) (x s ) = 0 for 1 k m and g (m+1) (x s ) 0. We say that g has stationary points if such points exist. The following lemma describes the asymptotic behaviour of a one dimensional oscillatory integral with a C phase when there are only a finite number of stationary points of order one. Lemma.4. Let h W 1,1 ([a,b]) and g C ([a,b]) be such that g has only a finite number of stationary points of order at most one. Then there exists a constant C, independent of h and depending continuously on the C norm of g, such that for > 1. b a e ig(x) h(x)dx C 1/ h W 1,1 ([a,b]) Proof. Let {s j } N j=1 denote the stationary points and δ = min j( g (s j ) ). Let ǫ be such that g (x) > δ/ for all x j U j, where U j = (u d j,uu j ) = B ǫ (s j ) [a,b]. Thenwe can make useof a version of Van der Corput s lemma (see [0, Corollary.6.8]), to obtain ( ) ( e ig(x) δ 1/ h(u h(x)dx U j 4 1/ u j ) + h U (x) ) dx. j

6 6 JORGE TEJERO Now let V j = (vj d,vu j ) denote each of the remaining segments of [a,b], such that j V j = [a,b]\ j U j. Integrating by parts in each V j we obtain V j e ig(x) h(x)dx = 1 i [ e ig(x) h(x) g (x) ] v u j v d j 1 e ig(x) d i V j dx ( h(x) g (x) ) dx. By Sobolev embedding [1, Theorem 4.1, Part 1, Case A] we have that h L [a,b] C h W 1,1 [a,b]. Making use of this and Hölder s inequality, altogether we obtain b ( e ig(x) h(x)dx δ 1/ 48N ) 1/ h W 1,1 ([a,b]) where a κ = max j +(N +1) κ 1 h W 1,1 ([a,b]), { (g ) 1 L (V j ) + g g (g ) L (V j ) which is finite, as there are no stationary points in j V j. For > 1 we can take ( ) δ 1/ C = 48N +(N +1)κ. As δ and κ dependcontinuously on the C norm of g, the proof is concluded... Piecewise W s,1 -potentials. We say that a curve C in the plane is contained in C m,α, with m 0 and 0 α 1, if there exists a finite collection ofboundedopensets{u j } N j=1 suchthatc N j=1 U j, andfunctions f j C m,α (R) such that C U j {(x,f j (x)) : x R} or C U j {(f j (y),y) : y R}. For α = 0, we write f C m,0 (R) whenever f and its derivatives up to order m are continuous and bounded; occasionally we will describe the curve as being simply C m. Consider two curves C 1 and C for which there is a finite cover by open sets {U j } N j=1 such that for each j either or C 1 U j {(x,f 1,j (x)) : x R} and C U j {(x,f,j (x)) : x R} C 1 U j {(f 1,j (y),y) : y R} and C U j {(f,j (y),y) : y R}. Then we define the distance between the two curves in C m,α norm as { { } } d(c 1,C ) = inf sup f1,j f,j C m,α j },

7 PIECEWISE SMOOTH RECOVERY 7 where the infimum is taken over all possible common covers. If a common cover does not exists for the curves, we write d(c 1,C ) =. A bounded Lipschitz domain whose boundary is a finite union of C m,α curves will be referred to as a piecewise C m,α domain. The potentials that we will consider exhibit discontinuities over C,α curves, where 1/ < α 1. More precisely, we are concerned with piecewise W s,1 -potentials V that can be expressed as N V(x) = q j (x)χ Ωj (x), j=1 where q j W s,1 (R ) and Ω j are piecewise C,α domains with 1/ < α 1. We will use the following norm for these potentials: N ( V D s,r = inf q j W s,1 1+ ) N χ H Ωj r : V(x) = q j (x)χ Ωj (x) j=1 where q j and Ω j are as previously described. The following lemma provides a bound for the L Sobolev norm for the potentials of our interest. Lemma.5. Let q : R R and let Ω be a bounded Lipschitz domain in the plane. Then there exists a constant C independent of q and Ω such that i) For 0 < r we have ii) For 0 < r < 1/ we have q H r C q W r+1,1. qχ Ω H r C qχ Ω D,r. Proof. Let m be the largest integer less than r, let t = r m and let D t = ( ) t/. For the first case we can use Sobolev embedding; see for example [1, Theorem 4.1, Part 1, Case C], and the fact that W r+1,1 W m+1,1 to obtain q H r ( q H m + D t q ) H 1/ C ( q m H m + D t q ) W 1/ m+1,1 ( ) 1/ C q H m + q W r+1,1 C q W r+1,1. For the second case we can use the generalised Leibniz rule [4, Theorem A.1], which states that j=1 D r (qχ Ω ) qd r (χ Ω ) χ Ω D r (q) L C q L D r (χ Ω ) L, and so by triangle inequality we obtain (3) qχ Ω H r C( q L χ Ω H r + q H r χ Ω L ). For the first term in the right-hand side, we can use Sobolev embedding; see for example [1, Theorem 4.1, Part 1, Case A], to obtain q L C q W,1.

8 8 JORGE TEJERO On the other hand, we have χ Ω H r for r < 1/; see for example [19]. For the second term in the right-hand side of (3) we can use case i), combined with the embedding W,1 W 1+r,1, concluding the proof..3. A topological property of C 1 graphs. Finally, we provide a simple continuity result that will be useful for characterizing the topological properties of the set of points where the reconstruction is not guaranteed as well as the continuity properties of the error bound of the reconstruction. Lemma.6. Let f C 1 [a,b] and {x j } N j=1 (a,b) be such that f(x j) = 0, f (x j ) 0 and f(x) 0 for all x [a,b] \{x j } N j=1. Then, for any ǫ > 0, there exists δ such that if f g C 1 < δ then there exists {x j }N j=1 (a,b) such that x j x j < ǫ, with g(x j ) = 0, g (x j ) 0 and g(x) 0 for all x [a,b]\{x j }N j=1. Proof. Let η = min j { f (x j ) } and let r be such that f (x) > η/ for all x B r (x j ). Let Then, for all g such that we have and a = inf { f(x) }. x [a,b]\ B r(x j ) f g C 1 < 1 min{a,η/}, g(x) > a/, x [a,b]\ j B r (x j ) g (x) > η/4, x j B r (x j ). Now, as f g C 0 < a/, we know that for all x [a,b] \ B r (x j ) we have f(x)g(x) > 0 (f and g have the same sign outside the balls B r (x j )), and so, by the intermediate value theorem, g must vanish in each of the balls B r (x j ). The fact that g only vanishes at a single point x j in each of the balls is a consequence of the fact that g is monotonous inside them. As g (x) > η/4 inside the balls, then, whenever f g C 0 < ǫη/4, we have that Taking concludes the proof. x j x j < ǫ. δ = min{a/,η/4,ǫη/4}

9 PIECEWISE SMOOTH RECOVERY 9 3. Recovery Later we will see that recovery is not guaranteed, even at points that lie far from the discontinuities of the potential. In order to bound the measure of these points, we will require the following key lemmas. Lemma 3.1. Let C be a C 1 curve contained in a bounded planar domain Ω. Then, the union of tangent lines to C with a fixed slope s has zero Lebesgue measure in Ω. Proof. Let C = {(x,γ(x))}, ǫ > 0 and I ǫ = {x : Γ(x) Ω, Γ (x) s < ǫ}. As Γ C 1, I ǫ is open, and as the real line satisfies the countable chain condition, I ǫ must consist of a countable union of disjoint intervals {U j }. For x 0,x 1 U j,x 0 < x 1, by the Fundamental Theorem of Calculus, we have Γ(x 1 ) = Γ(x 0 )+ x1 x 0 Γ (x)dx. As Γ (x) s < ǫ in the domain of integration, then and (x 1 x 0 )(s ǫ) < Γ(x 1 ) Γ(x 0 ) < (x 1 x 0 )(s+ǫ), (4) Γ(x 1 ) Γ(x 0 ) s(x 1 x 0 ) < ǫ(x 1 x 0 ). Now let l x0,s denote the line-segment with slope s that contains the point (x 0,Γ(x 0 )); l x0,s = {(x,y) Ω : y = Γ(x 0 )+s(x x 0 ))}. We write p 0 = (x 0,Γ(x 0 )) l x0,s and p 1 = (x 1,Γ(x 1 )) l x1,s. As p 0 l x0,s, then so is p 0 = (x 1,Γ(x 0 ) + s(x 1 x 0 )). By (4) we know that d(p 0,p 1) < ǫ(x 1 x 0 ) and so it follows that L j = x U j l x,s is contained in a rectangle of width bounded by ǫ U j and length bounded by diam(ω). Thus, L j L j < ǫ U j diam(ω) ǫdiam(ω). j j j Letting ǫ tend to zero, the proof is complete. Lemma 3.. Let Ω be a bounded domain in the plane, let C be a curve contained in Ω which is the graph of a C,α function with 1/ < α 1, and let φ x (z) = (x 1 z 1 ) (x z ). Then the set of points x Ω such that either i) φ x C has an infinite number of stationary points, ii) φ x C has at least one stationary point of order greater than one, has zero Lebesgue measure and is closed.

10 10 JORGE TEJERO Proof. First we see that if φ x C has an infinite number of stationary points, it has at least one stationary point of order greater than one, and so the first case is contained in the second. Let x be such that φ x C has an infinite number of stationary points. Then, by the compactness of C, there exists a sequence of stationary points {z 1,i } i=1 and a point z 1, such that with lim i z 1,i = z 1, φ x C z 1 (z 1, ) = 0. As φ x C is a C function and φx C z 1 vanishes at all the stationary points then we have φ x C φ x C z z1 (z 1, ) = lim 1 (z 1,i+1 ) φx C z 1 (z 1,i ) = 0, i z 1,i+1 z 1,i and therefore z 1, is a stationary point of order greater than one. Now we see that the set of x such that φ x C has a stationary point of order greater than one is null. That is, the set of x such that φ x C z 1 = φ x C z 1 has zero measure. Letting Γ C,α be such that C = {(z 1,Γ(z 1 ))}, the previous condition can be written as = 0 z 1 x 1 Γ (z 1 )(Γ(z 1 ) x ) = 1 Γ (z 1 )(Γ(z 1 ) x ) Γ (z 1 ) = 0 leading to (5) (6) x 1 = z 1 +Γ (z 1 )(x Γ(z 1 )), Γ (z 1 )(x Γ(z 1 )) = Γ (z 1 ) 1. Firstweconsiderthecasewhere Γ (z 1 ) > δ > 0. AsΓ iscontinuous and the real line satisfies the countable chain condition, for this to be satisfied z 1 must lie in one of at most countably many intervals. Taking one such interval U and rearranging (5) and (6), the set of x such that φ x C has a stationary point of order greater than one at z 1 U is given by ( x = (x 1,x ) = G(z 1 ) = z 1 + Γ (z 1 ) 3 Γ (z 1 ) Γ,Γ(z 1 )+ Γ (z 1 ) 1 (z 1 ) Γ (z 1 ) Weseethat theset ofxistheimageofac 0,α function. Toseethat suchaset has zero measure, take {U j } N j=1 a covering of U such that U j = U /N. Now as G(U j ) is contained in a ball of radius ( U /N) α, we obtain G(U) U α N 1 α. As α > 1/, we can let N tend to infinity to conclude that G(U) = 0. This is the only place where we require the Hölder regularity. Now as the countable union of null sets is null, we have concluded the proof in this case. ).

11 PIECEWISE SMOOTH RECOVERY 11 On the other hand, if Γ (z 1 ) δ and x Ω, then by (6) it follows that Γ (z 1 ) must be contained in the interval [1 δ,1+δ ] with δ = δdiam(ω). Let l z1,s be the line that passes through (z 1,Γ(z 1 )) with slope s, and let T(δ ) = l z1,1/γ (z 1 ), z 1 : Γ (z 1 ) [1 δ,1+δ ]. From equation (5) we see that the remaining set of x such that φ x C has a stationary point of order greater than one at (z 1,Γ(z 1 )) is contained in T(δ ). Using Lemma 3.1, we have lim T(1/n) Ω = n (T(1/n) Ω) = T(0) Ω n=1 = l z1,1 Ω + z 1 :Γ (z 1 )=1 z 1 :Γ (z 1 )= 1 l z1, 1 Ω = 0. Therefore, for any ε > 0 we can take δ small enough so that T(δ ) Ω < ε, allowing us to conclude that the set of points x such that φ x C has a stationary point of order greater than one is null. To see that the set is closed, first notice that for any δ > 0, there exists r > 0 such that for any x B r (x) we have φ x C φ x C C < δ. Thus, applying Lemma.6 to φx C z 1 we see that whenever φ x C has a finite number of stationary points of degree at most one then φ x C will have the same number of stationary points and of the same degree for any x close enough to x. This means that the set of points x such that φ x C has a finite number of stationary points of degree at most one is open, and the complement is closed, concluding the proof. Suppose that 0 is not a Dirichlet eigenvalue for the Hamiltonian V. Then, for each f H 1/ there exists a unique solution u H 1 to (1), and the DtN map can be defined by (7) Λ V [f],v Ω = Λ V [f]v Ω = V uv + u v Ω for any v H 1 (Ω). Theorem 1.1 is contained in the following result in which we also obtain decay rates for potentials with slightly more regularity. Theorem 3.3. Let V be a piecewise W s,1 -potential, with 0 < s < r < 1. Then for almost every x Ω, there exists a constant C x = C(x, Ω j ) such that e iψ (Λ V Λ 0 )[u,x ] V(x) π C x 1 s/( ) V D s,r + V D s,r Ω Ω

12 1 JORGE TEJERO whenever is sufficiently large. Moreover, if s = then we have lim e iψ (Λ V Λ 0 )[u,x ] = V(x), a.e. x Ω. π Ω Proof. As the DtN is a self-adjoint operator and e iψx satisfies the Laplace equation, then we can use the DtN map definition (7) to see that π Ωe iψx (Λ V Λ 0 )[u] = e iφx V (1+w), π Ω where φ x = ψ x +ψ x. Recalling that 0 s < r < 1, by Lemma.3 we have sup e iφx V w x Ω π C r V Hr, as, Ω and by part ii) of Lemma.5 this yields (8) sup e iφx V w π C1 s/ V Ds,r, as. x Ω Ω Now, by Lemma 3. we know that for almost every x in Ω, φ x Ωj has only a finite number of stationary points of order at most one. We now prove that the reconstruction formula recovers the potential correctly at these points for piecewise W s,1 -potentials, s > (almost all of them for W,1 -potentials). First we split the integral (9) π Ωe iφx V = π N j=1 Ω j e iφx q j, wherewe write φ for φ x from now on. We will prove that the value of each of these integrals tends to zero sufficiently fast whenever the integration domain does not contain x, the point at which we are reconstructing. Then we show that the value of the integral that contains x converges to V(x). Without loss of generality, we can suppose that x belongs to the interior of Ω 1. For j > 1, we use Green s first identity, with u = eiφ i and v = q j φ, to φ obtain Ω j e iφ q j = 1 i Ω j e iφ q j φ φ n 1 i Ω j e iφ ( ) qj φ φ. Using Green s ( first ) identity again on the second term with u = eiφ i and qj φ v = φ leads to φ φ e iφ q j = 1 e iφ q j φ π Ω j iπ Ω j φ n + 1 ( ) e iφ qj φ φ (10) π Ω j φ φ n 1 ( ( ) ) e iφ qj φ φ π Ω j φ φ.

13 PIECEWISE SMOOTH RECOVERY 13 As the number of stationary points on Ω j is finite and are of order at most one, and by trace theorem we know that q j Ωj W 1,1 ( Ω j ) (see for example [18, Section 5.5, Theorem 1]), then we can use Lemma.4 on each of the C components of Ω j, together with Hölder s inequality, to see that there exists C x = C(x, Ω j ) such that (11) e iφ q j φ Ω j φ n C x 1/ φ n φ W 1, ( Ω j ) q j W 1,1 ( Ω j ), as. As x belongs to the interior of Ω 1, we have that φ is bounded, and using the trace theorem this yields to (1) e iφ q j φ Ω j φ n Cx 1/ V D s,r, as. For the second term on the right-hand side of (10) we can use Hölder s inequality to obtain ( ) e iφ qj φ φ Ω j φ φ n qj φ φ n φ L 1 ( Ω j ) φ L ( Ω j ) ( φ Ẇ1, ( Ω j ) φ ) L + φ ( Ω j ) φ L ( Ω j ) φ n φ q j W 1,1 ( Ω j ), L ( Ω j ) and by the trace theorem we get ( ) (13) e iφ qj φ φ Ω j φ φ n C x V D s,r. Similarly, for the last term on the right-hand side of (10) we have ( ( ) ) e iφ qj φ φ Ω j φ φ q j φ φ Ẇ1,1 ) φ Ẇ1, (Ω j ) φ L (Ω (Ω j ) j + q j φ φ φ φ L 1 (Ω j ) L (Ω j )

14 14 JORGE TEJERO ( φ W, (Ω j ) φ W + φ 1, (Ω j ) φ φ Ẇ1, (Ω j ) φ L q j (Ω j ) W,1 (Ω j ) ( + φ Ẇ1, (Ω j ) φ L + φ (Ω j ) φ φ φ q j W 1,1 (Ω j ) L (Ω j ) ( φ W, (Ω j ) φ W + φ 1, (Ω j ) φ q j W,1 (Ω j ) Ẇ1, (Ω j ) L (Ω j ) ) W 1, (Ω j ) ) ) yielding (14) e iφ Ω j ( ( ) ) qj φ φ φ φ C x V D s,r. Plugging (1), (13) and (14) into (10) we obtain (15) e iφ q j π Ω j C x 1/ V D s,r as. We now consider Ω 1 e iφ q 1 by decomposing q 1 into q x = q 1 χ, q rem = q 1 (1 χ), where χ(z) is a bump function such that { 1 if z x r1, χ(z) = 0 if z x r, with 0 < r 1 < r < d(x, Ω 1 ). As q rem (y) = 0 for y close enough to x, we can use the same arguments that lead to (15) to obtain (16) e iφ q rem π C x 1/ q rem W,1, as. Ω 1 Ontheother hand, as q 1 W,1 (Ω 1 ), wecan usesobolev embedding(see for example [1, Theorem 4.1, Part 1, Case C]) to see that q x H0 1(Ω 1). Now it was noted in [4] that π Ω 1 e iφ q x can be interpreted as the solution to a nonelliptic time dependent Schrödinger equation, at time 1/. Thus, using the almost everywhere convergence result of [3, Theorem 1] we obtain (17) lim e iφ q x = q x (x) = q 1 (x) = V(x) a.e. x R. π Ω 1

15 PIECEWISE SMOOTH RECOVERY 15 If q 1 W s,1 with s >, then we can recover at all the remaining points and we get a decay rate. Indeed, e iφ q x q x (x) π = ( q x e i(z 1 z)) (x) q x (x) Ω 1 π = 1 ( ) 4π e ix ξ q x (ξ) e i 1 (ξ 1 ξ ) 1 dξ e i 1 (ξ 1 ξ ) 1 1/ q x Ḣs 1 dξ Note that for < s < 3 we have ξ s ( cos(ξ = q x Ḣs 1 1 s/ 1 ξ ) ) 1/ ξ s dξ ( sin = q x Ḣs 1 1 s/ ( 1 (ξ 1 ξ )) ) 1/ ξ s dξ. sin(ξ1 lim ξ ) sin( ξ ) ξ 0 ξ s 1 lim ξ 0 ξ = 1 thus the integral is finite, and we can use part i) of Lemma.5 to obtain (18) e iφ q x q x (x) π C1 s/ q x W s,1, as. Ω 1 Plugging (15), (16), (17), and (18) into (9), together with (8) concludes the proof. Remark 3.4. As noted in [4], π Ω eiφ V can be interpreted as the solution to a nonelliptic time dependent Schrödinger equation at time 1/. Therefore equations (15), (16), (17), and (18) imply almost everywhere convergence to the initial data V, whenever V is piecewise-w s,1 with s < 3. A consequence of Theorem 3.3 is that potentials of this type can also be recovered from the scattering data at a fixed energy. Corollary 3.5. Let V be a piecewise W,1 -potential. Then V can be recovered almost everywhere from the scattering amplitude at a fixed energy k. Proof. Let Q be a square such that N j=1 Ω j Q. In [5] expressions are given for computing Λ V k defined on Q from the scattering amplitude at energy k. Therefore, the recovery from the scattering amplitude follows directly from the fact that if V is piecewise W,1 -potential, then so is V k χ Q, which allows us to recover the potential using Theorem 3.3. For the stability estimates of the sequel we will require some continuity properties of the constant in Theorem 3.3 which we record as a lemma.

16 16 JORGE TEJERO Lemma 3.6. Let N be the set of x such that φ x Ωj has a stationary point of degree greater that one. Then the constant C x = C(x, Ω j ) in Theorem 3.3 has the following continuity properties in Ω\N: i) It is continuous with respect to x. ii) It is continuous with respect to Ω j in the C norm. Proof. Let N, Ω j and r 1 be as in the proof of Theorem 3.3. Let Ω j = Ω j for j =,...,N and let Ω 1 = Ω 1 \B r1 (x). The constant C x is given by C x = N j=1 C x φ n φ W 1, ( Ω j ) + φ n φ φ φ Ẇ1, ( Ω j ) L ( Ω + φ φ W, (Ω j ) W 1, (Ω j ) + j ) + L ( Ω j ) φ φ φ φ L ( Ω j ) W 1, (Ω j ). The constant C x appears in equation (11) by the use of Lemma.4 (taking g = φ Ωj and h = q j φ n Ωj ), so it is continuous with respect to x and φ with respect to Ω j in the C norm, and as φ does not vanish inside any of the Ω j, then so is C x, concluding the proof. Asistobeexpected, theerrorinthereconstructionincreases thecloser we movetothediscontinuities ofthepotential, astheconstant C x blowsup, and we are unable to recover at the discontinuities. It is perhaps more interesting that, for certain potentials, there are points where the reconstruction fails which are far from the discontinuities of the potential. Indeed, consider thepotential given by V = χ Ω1, whereω 1 is the rhombus with vertices at (0,0),(1,1),(,0) and (1, 1); see Figure 1. Consider the problem of recovering the potential inside Ω = [,] [,]. We might expect tobeableto recover at thepoints x = ( t, t), fort (0,), farfrom the potential. However, by Alessandrini s identity [, Lemma 1], we know that the reconstructed potential Ṽ at x is the limit as tends to infinity of π Ωe iψx (Λ V Λ 0 )[u,x ] = π where φ x = ψ x +ψ x, which can be rewritten as ( e iφx V w,x + e iφx V + π Ω Ω\Ω 1 Ω e iφx V (1+w,x ), Ω 1 e iφx V ).

17 PIECEWISE SMOOTH RECOVERY 17 For the firstterm we can uselemma.3 and part ii) of Lemma.5 to obtain e iφx V w,x C r V D,r, as, π Ω for any 0 < r < 1/. As the potential is equal to 1 inside Ω 1 and zero in the rest of the domain, this yields Ṽ(x) = lim e iφx. π Ω 1 Using Green s first identity twice we get e iφx = 1 φ e iφx x π Ω 1 iπ Ω j φ x n + 1 ( ) φx φx e iφx π Ω 1 φ x φ x n 1 ( ( ) ) φx φx e iφx π Ω 1 φ x φ x. As we have seen in the proof of Theorem 3.3, the second and third terms converge to zero as tends to infinity, and for the first term we can write φ x φ x n = φ x φ x n φ x 1 + φ x n Ω j e iφx where l 1 e iφx + l 3 e iφx φ x φ x n 3 + l 1 = (s,s) for s (0,1), l = (1+s,1 s) for s (0,1), l 3 = ( s, s) for s (0,1), l 4 = (1 s,s 1) for s (0,1). As we have φ(z) = (z 1 +t) (z +t), we see that φ x l1 (s) = 0, φ x l (s) = 4s(t+1), φ x l3 (s) = 4(t s+1), φ x l4 (s) = 4t(1 s). l e iφx l 4 e iφx Therefore, we can apply Lemma.4 to three of the sides; φ x φ x n j = O( 1/ ), for j =,3,4, l j e iφx and on the other hand we have φ x φ x n 1 = l 1 e iφx 1 0 φ x φ x n 4 4(s+t) ds = 4 (log(t) log(t+1)).

18 18 JORGE TEJERO Figure 1. Solid lines are the discontinuities of the potential and dashed lines are points far from the discontinuities where the recovery fails. Putting everything together we obtain Ṽ(x) = i 4π log(1+1/t) Stability We begin with some preliminary results that we will require for the proof of the stability estimates. Lemma 4.1. Let V 1,V H s, where 0 < s 1 and let F L. Then there exists a constant C such that suptw 1 w [F](x) C s F L V 1 Ḣs V Ḣs x Ω when is sufficiently large. Proof. By Hölder s inequality and the Hardy-Littlewood-Sobolev lemma we get Tw 1 w [F](x) Fw 1 w L 1 F L w 1 L 4 w L 4 C F L w 1 Ḣ1/ w Ḣ1/.

19 PIECEWISE SMOOTH RECOVERY 19 As (I SV ) 1 is bounded for sufficiently large (see [4, Lemma.3]), this yields Tw 1 w [F](x) C F L S1[V 1 ] Ḣ1/ S1[V ] Ḣ1/. Applying Lemma. twice concludes the proof. We will also require the following crude bound for Bukhgeim solutions. Lemma 4.. Let 0 < r < 1/, and let V be a piecewise W,1 -potential contained in a bounded planar domain Ω with diameter d. Then there exists a constant C depending on Ω such that the Bukhgeim solutions satisfy u,x H 1 (Ω) Ced ( 1+ V D,r ) whenever is sufficiently large. Proof. Writing u = u,x, ( ) u H 1 (Ω) C u L (Ω) + zu L (Ω) + zu L (Ω) = Ce iψ (1+w) +C 1 L z Ve iψ (1+w) (Ω) +C z 1 Ve iψ (1+w) L (Ω) L (Ω) Note that z 1 and z 1 are bounded operators, as V has compact support (see for example [6, Theorem 4.3.1]). As we have L (Ω) C L 4 (Ω), then using Hölder s inequality leads to u H 1 (Ω) Ce iψ L (1+ V L 4) 1+w L (Ω) 4 (Ω) ( ) Ce id (1+ V L 4) 1+ w L 4 (Ω), and by the Hardy-Littlewood-Sobolev inequality we get u H 1 (Ω) Ce id (1+ V L 4) ( 1+ w Ḣ1/). As (I SV ) 1 is bounded for sufficiently large (see [4, Lemma.3]), we have ( ) u H 1 (Ω) Ceid (1+ V L 4) 1+ S1 [V] Ḣ1/. Now, by Lemma. we get u H 1 (Ω) Ceid (1+ V L 4) ( 1+ V Ḣr). Using Sobolev embedding(see for example[1, Theorem 4.1, Part 1, Case A]) we have N V L 4 C q j L 4 V D,r. j=1 Using part ii) of Lemma.5 Lemma concludes the proof.

20 0 JORGE TEJERO We now prove an a priori stability estimate for reconstruction from the DtN map in the L norm. Note that the result is interesting when there is a priori knowledge of where, approximately, the discontinuities lie, as the constant term depends on the point under consideration with respect to the discontinuities. The result has been stated in the following form as in practical situations one could consider where a potential might lie, given a noisy reconstruction of it. If the assumption that the discontinuities of the potentials are closed to each other was dropped, the constant C x would depend on both { Ω 1,j } and { Ω,j }, the discontinuities of V 1 and V respectively. Theorem 4.3. Let 0 < s < r < 1 and let V 1,V be piecewise W s,1 - potentials supported on a Lipschitz domain Ω in R such that their discontinuities are close enough with respect to the C norm. Then, for almost every x Ω, there exists a constant C x = C(x, Ω 1,j ) such that V 1 (x) V (x) C x ln Λ V1 Λ V 1 s/ (κ+κ 3 ) wheneverλ V1 and Λ V are close enough, where κ = max{ V 1 D s,r, V D s,r}. Proof. Let d be the diameter of Ω and let (19) = 1 6d ln Λ V 1 Λ V. Note that whenever Λ V1 and Λ V are sufficiently close, then can be as large as we need. By the triangle inequality, V 1 (x) V (x) V 1 (x) T1+w 1 [V 1 ](x) + V (x) T1+w (0) [V ](x) + T1+w 1 [V 1 ](x) T1+w [V ](x). We can use Theorem 3.3 on the first two terms to obtain (1) V j (x) T1+w j [V j ](x) C x 1 s/ (κ+κ ) wherec x = C(x, Ω 1,j ). We can take the same constant C x for both terms as it is continuous with respect to the discontinuities in the C norm (due to Lemma 3.6). For the last term we have T1+w 1 [V 1 ](x) T1+w [V ](x) T(1+w 1 )(1+w ) [V 1 V ](x) + Tw 1 [V ](x) + Tw [V 1 ](x) + Tw 1 w [V 1 V ](x). By Lemma.3 and part ii) of Lemma.5 we obtain () Tw [V 1 ] C1 s/ κ and T L w 1 [V ] C1 s/ κ. L

21 PIECEWISE SMOOTH RECOVERY 1 and by Lemma 4.1 and part ii) of Lemma.5 we obtain (3) Tw 1 w [V 1 V ] C1 s/ κ 3. L Let u j be Bukhgeim solutions to u j = V j u j. Then we have T(1+w 1 )(1+w ) [V 1 V ] = L L π 1 V )u 1 u Ω(V = π f 1 Λ V [f ] f Λ V1 [f 1 ] Ω L where f j = u j Ω. As the DtN is a self-adjoint operator, we have that π f 1 Λ V [f ] f Λ V1 [f 1 ] = L L π (Λ V1 Λ V )[f 1 ]f. Ω For f H 1/ the DtN map satisfies Λ v [f] H 1/, where H 1/ is the dual of H 1/. Thus, for any x Ω, we have π (Λ V1 Λ V )[f 1 ]f Λ V1 Λ V f 1 H 1/ ( Ω) f H 1/ ( Ω) Ω Λ V1 Λ V u 1 H 1 (Ω) u H 1 (Ω) and we can use Lemma 4. to obtain (4) π (Λ V1 Λ V )[f 1 ]f Ce d Λ V1 Λ V ( 1+κ 4). Ω Inserting (1), (), (3) and (4) into (0), and noting that < e d for sufficiently large, leads to V 1 (x) V (x) C x 1 s/ (κ+κ 3 )+ C Λ V1 Λ V e 3d (1+κ 4 ). Taking as in (19) we obtain ( ) 1 1 s/ V 1 (x) V (x) C x 6d ln Λ V 1 Λ V (1+κ 3 ) Ω +C Λ V1 Λ V 1/ (1+κ 4 ) where the second term can be omitted for Λ V1 Λ V small enough, concluding the proof. We now provide a link to the scattering amplitude. We adapt the proof of Stefanov (see [33]) to the two-dimensional case. Due to the severe illposedness of the problem, we need to introduce a norm for the scattering amplitude which penalizes the higher components of the frequency spectrum. Let V be a potential supported on the unit disk, then we define the norm for its scattering amplitude at a fixed energy k as A V k = n,m Z ( ) 3+3 n n ( ) 3+3 m m a (n,m) V k k 1/,,

22 JORGE TEJERO where a (n,m) V are the Fourier coefficients of A V A V (η,θ) = a (n,m) V e inη+imθ. n,m Z Before passing to the proof, we define the single layer potential operator S V [f](x) = G V (x,y)f(y)dy, where G V is the outgoing Green s function which satisfies Ω ( +V k )G V (x,y) = δ(x y). Lemma 4.4. Let V 1,V be two potentials supported in the unit disk. Then there exists a constant C k = C(k) such that Λ V1 Λ V H 1/ (S 1 ) H 1/ (S 1 ) C k A V1 A V k. Proof. Using Nachman s formula [7, Theorem 1.6] we have Λ V1 k Λ V k = S 1 V 1 S 1 V = S 1 V 1 (S V S V1 )S 1 V. As S V is a bounded and invertible mapping from H 1/ (S 1 ) to H 1/ (S 1 ) (see [3, Proposition A.1]), we have ΛV1 k Λ H V k 1/ (S 1 ) H 1/ (S 1 ) C S V 1 S V H 1/ (S 1 ) H 1/ (S 1 ). Letting B x = (1 x ) 1/4, we write S V1 S V H 1/ (S 1 ) H 1/ (S 1 ) = sup (G V1 (x,y) G V (x,y))f(y)dy f =1 S 1 H 1/ (S 1 ) ( ) 1/ = sup (B x (G V1 (x,y) G V (x,y))f(y)dy dx) f =1 S 1 S 1 ( ( ) 1/ = sup B x (G V1 (x,y) G V (x,y))f(y)dy dx) f =1 S 1 S 1 by Pareseval s identity we have ( ( = sup (B y B x (G V1 (x,y) G V (x,y)))by f(y)dy) 1 dx f =1 S 1 S 1 ) 1/ using Minkowski s integral inequality we get ( ) ( 1/ sup (By B x (G V1 (x,y) G V (x,y)))by 1 f(y)) dx dy f =1 S 1 S 1 = sup f =1 S 1 B 1 y f(y) ( ) 1/ (B y B x (G V1 (x,y) G V (x,y))) dx dy S 1

23 and using the Cauchy-Schwarz inequality PIECEWISE SMOOTH RECOVERY 3 sup f H 1/ (S 1 ) G V 1 G V H 1/ (S 1 ) H 1/ (S 1 ) f =1 = G V1 G V H 1/ (S 1 ) H 1/ (S 1 ). From [5, Theorem.] we know that G V1 (x,y) G V (x,y) = n,m Z ( 1) n 16 in+m( a (n,m) V 1 ) a (n,m) V H (1) n (k x )H (1) m (k y )e inφx+imφy where H (1) denotes the Hankel function of the first kind. Now, using Parseval s identity and the bound for the Hankel function in [5, Lemma.3], there exists C k = C(k) such that G V1 G V H 1/ (S 1 ) H 1/ (S 1 ) (1+n ) 1/ (1+m ) 1/ a (n,m) C k C k n,m Z n,m Z n,m Z concluding the proof. V 1 (1+n ) 1/ (1+m ) 1/ a (n,m) ( ) 3+3 n n ( ) 3+3 m m a (n,m) k k V 1 a (n,m) V H (1) n (k) H (1) m (k) ( ) a (n,m) V n! m! 3 n + m k V 1 a (n,m) V Corollary 4.5. Let 0 < s < r < 1 and let V 1,V be two piecewise W s,1 potentials supported on a bounded domain in R such that their discontinuities are close enough in the C norm. Then, for almost every x Ω, there exist constants C x = C(x, Ω j ),C k = C(k), such that V 1 (x) V (x) C x ln(ck A V1 A V k ) 1 s/ (κ+κ 3 ) wheneverλ V1 and Λ V are close enough, where κ = max{ V 1 D s,r, V D s,r}. Acknowledgements: The author thanks Daniel Faraco and Keith Rogers for their valuable comments and corrections and thanks Evgeny Lakshtanov for bringing [6] to his attention. This work has been partially supported by ERC and by ICMAT Severo Ochoa project SEV (MINECO). References [1] R. A. Adams and John J. F. Fournier. Sobolev Spaces. Pure and Applied Mathematics, Elsevier (003). [] G. Alessandrini. Stable determination of conductivity by boundary measurements. Applicable Analysis (1988). [3] G. Alessandrini. Global stability for a coupled physics inverse problem. Inverse Problems (014).

24 4 JORGE TEJERO [4] K. Astala, D. Faraco and K. M. Rogers. Unbounded potential recovery in the plane. Annales scientifiques de l École Normale Supérieure, to appear, arxiv: (015). [5] K. Astala, D. Faraco and K. M. Rogers. Recovery of the Dirichlet-to-Neumann map from scattering data in the plane. RIMS Kokyuroku Bessatsu (014). [6] K. Astala, T. Iwaniec and G. Martin. Elliptic Partial Differential Equations and Quasiconformal Mappings in the Plane. Princeton University Press (008). [7] K. Astala and L. Päivärinta. Calderón Inverse Conductivity problem in plane. Annals of Mathematics (006). [8] J. A. Barceló, T. Barceló and A. Ruiz. Stability of the Inverse Conductivity Problem in the Plane for Less Regular Conductivities. Journal of Differential Equations (001). [9] T. Barceló, D. Faraco and A. Ruiz. Stability of Calderón inverse conductivity problem in the plane. Journal de Mathématiques Pures et Appliquées (007). [10] E. Blasten, O. Yu. Imanuvilov and M. Yamamoto. Stability and uniqueness for a two-dimensional inverse boundary value problem for less regular potentials. Inverse Problems and Imaging (015). [11] R. M. Brown and R. H. Torres. Uniqueness in the inverse conductivity problem for conductivities with 3/ derivatives in L p,p > n. Journal of Fourier Analysis and Applications (003). [1] R. M. Brown and G. A. Uhlmann. Uniqueness in the inverse conductivity problem for nonsmooth conductivities in two dimensions. Communications in Partial Differential Equations (1997). [13] A. L. Bukhgeim. Recovering a potential from Cauchy data in the two-dimensional case. Journal of Inverse and Ill-Posed Problems (008). [14] A. P. Calderón. On an inverse boundary value problem. Seminar on Numerical Analysis and its Applications to Continuum Physics, Rio de Janeiro (1980). [15] P. Caro, A. García and J. M. Reyes. Stability of the Calderón problem for less regular conductivities. Journal of Differential Equations (013). [16] P. Caro and K. M. Rogers. Global uniqueness for the Calderón problem with Lipschitz conductivities. Forum of Mathematics Pi (016). [17] A. Clop, D. Faraco and A. Ruiz. Stability of Calderón s inverse conductivity problem in the plane for discontinuous conductivities. Inverse Problems and Imaging (010). [18] L. C. Evans. Partial Differential Equations. Graduate Studies in Mathematics (010). [19] D. Faraco and K. M. Rogers. The Sobolev norm of characteristic functions with applications to the Calderón inverse problem. The Quarterly Journal of Mathematics (013). [0] L. Grafakos. Classical Fourier analysis. Springer, Graduate Texts in Mathematics (008). [1] B. Haberman. Uniqueness in Calderón s problem for conductivities with unbounded gradient. Communications on Mathematical Physics (015). [] B. Haberman, D. Tataru. Uniqueness in Calderón s problem with Lipschitz conductivities. Duke Mathematical Journal (013). [3] V. Isakov and A. I. Nachman. Global uniqueness for a two-dimensional semilinear elliptic inverse problem. Transactions of the American Mathematical Society (1995). [4] C. Kenig, G. Ponce and L. Vega. Well-posedness and scattering results for the generalized Korteweg-de Vries equation via the contraction principle. Communications in Pure and Applied Mathematics (1993). [5] C. Kenig, J. Sjöstrand and G. Uhlmann. The Calderón problem with partial data. Annals of Mathematics (007). [6] E. L. Lakshtanov, R. G. Novikov and B. R. Vainberg. A global Riemann-Hilbert problem for two-dimensional inverse scattering at fixed energy. arxiv: (015). [7] A. I. Nachman. Reconstructions from boundary measurements. Annals of Mathematics (1988).

25 PIECEWISE SMOOTH RECOVERY 5 [8] A. I. Nachman. Global uniqueness for a two-dimensional inverse boundary value problem. Annals of Mathematics (1996). [9] R. G. Novikov. Multidimensional inverse spectral problem for the equation ψ + (v(x) Eu(x))ψ = 0. Functional Analysis and Its Applications (1988). [30] R. G. Novikov and M. Santacesaria. A global stability estimate for the Gel fand- Calderón inverse problem in two dimensions. Journal of Inverse and Ill-Posed Problems (010). [31] R. G. Novikov and M. Santacesaria. Global uniqueness and reconstruction for the multi-channel Gel fand-calderón inverse problem in two dimensions. Bulletin des Sciences Mathématiques (011). [3] K. M. Rogers, A. Vargas and L. Vega. Pointwise convergence of solutions to the nonelliptic Schrödinger equation. Indiana University Mathematical Journal (006). [33] P. Stefanov. Stability of the inverse problem in potential scattering at fixed energy. Annales de l institue Fourier (1990). [34] J. Sylvester and G. Uhlmann. A global uniqueness theorem for an inverse boundary value problem. Annals of Mathematics (1987).