Uniqueness in the inverse conductivity problem for complex-valued Lipschitz conductivities in the plane
Evgeny Lakshtanov, Jorge Tejero, Boris Vainberg

TL;DR
This paper proves that in the plane, the conductivity in impedance tomography can be uniquely identified from boundary measurements when the conductivity is complex-valued and Lipschitz continuous.
Contribution
It introduces a novel approach using Bukhgeim's scattering data for the Dirac problem to establish uniqueness in the inverse conductivity problem.
Findings
Conductivity is uniquely determined by the Dirichlet-to-Neumann map.
The method applies to complex-valued Lipschitz conductivities in the plane.
The approach advances the theoretical understanding of inverse boundary value problems.
Abstract
We consider the impedance tomography problem in the plane. Using Bukhgeim's scattering data for the Dirac problem, we prove that the conductivity is uniquely determined by the Dirichlet-to-Neuman map
Peer Reviews
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.
Taxonomy
TopicsNumerical methods in inverse problems · Mathematical Analysis and Transform Methods · Advanced Mathematical Modeling in Engineering
Uniqueness in the inverse conductivity problem for complex-valued Lipschitz conductivities in the plane
Evgeny Lakshtanov Department of Mathematics, Aveiro University, Aveiro 3810, Portugal. This work was supported by Portuguese funds through the CIDMA - Center for Research and Development in Mathematics and Applications and the Portuguese Foundation for Science and Technology (“FCT–Fundção para a Ciência e a Tecnologia”), within project UID/MAT/0416/2013 ([email protected]).
Jorge Tejero Universidad Autónoma de Madrid and ICMAT. The work was partially supported by the ERC grants ERC-301179 and ERC-277778 and the MINECO grants SEV-2015-0554 and MTM-2013-41780-P (Spain) ([email protected]).
Boris Vainberg Department of Mathematics and Statistics, University of North Carolina, Charlotte, NC 28223, USA. The work was partially supported by the NSF grant DMS-1410547 ([email protected]).
Abstract
We consider the inverse impedance tomography problem in the plane. Using Bukhgeim’s scattering data for the Dirac problem, we prove that the conductivity is uniquely determined by the Dirichlet-to-Neuman map.
Key words: Bukhgeim’s scattering problem, inverse Dirac problem, inverse conductivity problem, complex conductivity.
1 Introduction
Let be a bounded domain in . The electrical impedance tomography problem (e.g., [6]) concerns determining the impedance in the interior of , given simultaneous measurements of direct or alternating electric currents and voltages at the boundary . If the magnetic permeability can be neglected, then the problem can be reduced to the inverse conductivity problem (ICP), i.e., to the problem of reconstructing function , from the set of data , dense in an adequate topology, where
[TABLE]
Here is the unit outward normal to , , where is the electric conductivity and is the electric permittivity. If the frequency is negligibly small, then one can assume that is a real-valued function, otherwise it is supposed to be a complex-valued function.
An extensive list of references on the tomography problem can be found in the review [6]. Here we will mention only the papers that seem to be particularly related to the present work.
For real , the inverse conductivity problem has been reduced to the inverse problem for the Schrödinger equation. The latter was solved by Nachman in [14] in the class of twice differentiable conductivities. Later, Brown and Uhlmann [7] reduced the ICP to the inverse problem for the Dirac equation, which has been solved in [4], [15]. This approach requires the existence of only one derivative of . The authors of [7] proved the uniqueness for the ICP. Later, Knudsen and Tamasan [11] extended this approach and obtained a method to reconstruct the conductivity. Finally, the ICP has been solved by Astala and Paivarinta in [3] for real conductivities when both and are in .
If a complex conductivity has at least two derivatives, then one can reduce equation (1) to the Schrödinger equation and apply the method of Bukhgeim [8] (or some of the works extending this method, such as [5], [12] or [16]). This approach does not work in the case of only one time differentiable complex valued conductivities. On the other hand, the work of Francini [10], where the ideas of [7] were extended to deal with complex conductivities with small imaginary part, are not applicable to general complex conductivities due to possible existence of the so called exceptional points. In [13], Lakstanov and Vainberg extended the ideas of [12] to apply the -method in the presence of exceptional points and reconstructed generic conductivities under the assumption that and (here is the Fourier transform).
In this paper, we will prove that complex-valued Lipschitz conductivities are uniquely determined by information on the boundary. Since we use the standard reduction of (1) to the Dirac equation followed by the solution of the inverse problem for the Dirac equation, the condition on can be restated in the form , where is the potential in the Dirac equation. Our present result is based on a development of the Bukhgeim approach, combined with some of the arguments of Brown and Uhlmann from [7]. The statement of our main theorem is the following.
Theorem 1.1**.**
Let be a bounded Lipschitz domain in the plane and let be complex-valued Lipschitz conductivities. Then
[TABLE]
where is the Dirichlet-to-Neumann map for the conductivity .
The Dirichlet-to-Neumann (DtN) map is defined by
[TABLE]
where is a solution to (1) and is the normal derivative of at the boundary of . Function is defined as such an element of the space dual to that
[TABLE]
for each .
In section 2, we will describe our approach, stating the most relevant results. All the proofs will be given in section 3.
2 Main steps
2.1 Reduction to the Dirac equation
From now on, we will consider as a point of a complex plane: , and will be considered as a domain in . The following observation made in [7] plays an important role. Let be a solution of (1) and let . Then the pair satisfies the Dirac equation
[TABLE]
where
[TABLE]
Thus the inverse Dirac scattering problem is closely related to the ICP. If is found and the conductivity is known at one point , then in can be immediately found from (5).
From now on, we will use a different form of equation (2): instead of Beals-Coifmann notations , we will rewrite the equation in Sung notations: . We will consider the equation in the whole plane by extending the potential outside by zero. Then the vector is a solution of the following system
[TABLE]
where
[TABLE]
2.2 Solving the Dirac equation for large
Let be a matrix solution of (6) that depends on parameter and has the following behavior at infinity
[TABLE]
Note that the unperturbed wave
[TABLE]
depends on the spacial parameter and the spectral parameter , and grows at infinity exponentially in some directions. The same is true for the elements of the matrix . Let us stress that, contrary to the standard practice, we consider function (and other functions defined by ) for all complex values of , not just for . This allows us to generalize the Bukhgeim method to the case of potentials in . From the technical point of view, this allows us to use the Hausdorff-Young inequality.
Problem (6)-(10) can be rewritten using a bounded function
[TABLE]
i.e., (6)-(10) is equivalent to
[TABLE]
Using the fact that , equation (13) can be reduced to the Lippmann-Schwinger equation
[TABLE]
where and as .
Denote
[TABLE]
Then equation (14) implies that
[TABLE]
In particular, for the component of the matrix , we have , with leading to
[TABLE]
By inverting , we can obtain . Other components of can be found similarly.
Denote by the space of bounded functions of with values in a Banach space . The following two lemmas show that is a contractive operator in the space if is large enough, and that also belongs to this space. After these lemmas are proved, one can find the solution of (14) (using, for example, the Neumann series for the inversion of ). Then formula (12) provides the solution of (6)-(10).
Lemma 2.1**.**
Let . Then
[TABLE]
Lemma 2.2**.**
Let . Then there exists such that
[TABLE]
Note that (17) together with Lemmas 2.1 and 2.2 allows one to solve the direct but not the inverse problem, since operator depends on . The following inclusion is an immediate consequence of (17) and Lemmas 2.1 and 2.2:
[TABLE]
for large enough .
2.3 Determination of the potential
Let the matrix be the (generalized) scattering data, given by the formula
[TABLE]
One can use Green’s formula
[TABLE]
to rewrite as
[TABLE]
Thus, one does not need to know the potential in order to find . Function can be evaluated if the Dirichlet data is known for equation (6), since in (20) can be expressed via using (12).
The spectral parameter with real was used in the standard approach to recover the potential from scattering data (19), and the potential was recovered by the limit of the scattering data as . Instead, in the present work, we have , and the potential is determined by integrating the scattering data over a large annulus in the complex -plane.
Let be the operator defined by
[TABLE]
where can be a matrix- or scalar-valued function. Then
[TABLE]
We will show that the following statement is valid.
Theorem 2.3**.**
Let be a complex-valued bounded potential. Then
[TABLE]
and
[TABLE]
for every smooth with a compact support in . Thus
[TABLE]
Therefore, if the scattering data is uniquely determined by the DtN map, then so is the potential .
In order to prove (23), we use the two lemmas stated below and (17) rewritten as follows
[TABLE]
(other entries of the matrix can be handled in a similar way). Relation (24) follows from the stationary phase approximation.
Lemma 2.4**.**
Let . Then there exists such that
[TABLE]
Lemma 2.5**.**
Let . Then there exists such that
[TABLE]
3 Proofs
In order to make the calculations more compact, we introduce the following notation for the -space on the complement of the ball:
[TABLE]
We will also use the real-valued function
[TABLE]
where the dependence on and will be omitted in some cases.
3.1 Preliminary results
Lemma 3.1**.**
Let . Then the following estimate is valid for an arbitrary and some constants and :
[TABLE]
Remark. A more accurate estimate will be proved below with if , and with the right-hand side replaced by when , or by a constant when .
Proof. The statement is obvious if . If , then the left-hand side in the inequality above takes the following form after the substitution :
[TABLE]
Without loss of the generality, one can assume that . We split the function into two terms obtained by multiplying by and , respectively, where is the indicator function of the disk of radius two. The norm of can be estimated from above by an -independent constant. The second function can be estimated from above by . The norm of the latter function can be easily evaluated, and it does not exceed a constant if . It does not exceed if , and it does not exceed if . Since , we can replace in (26) by , and this implies the statement of the lemma. ∎
Lemma 3.2**.**
Let , and . Then
[TABLE]
where constant depends only on the support of and on .
Proof. Denote by the integral in the left-hand side of the inequality above. We change variables in and take into account that . Then
[TABLE]
Using the Hausdorff-Young inequality with and Lemma 3.1, we obtain that
[TABLE]
∎
3.2 Proof of Lemma 2.1
Let
[TABLE]
so that
[TABLE]
Then, from the Minkowski’s integral inequality, we have
[TABLE]
Thus it remains to show that, uniformly in and , we have
[TABLE]
Let be given by (27) with the extra factor in the integrand, where outside of a neighborhood of the origin, and vanishes in a smaller neighborhood of the origin. Since
[TABLE]
for each there exists such that
[TABLE]
for all the values of . Denote by the function with potentials replaced by their -approximations . Since the other factors in the integrand of are bounded (they are infinitely smooth), we can choose these approximations in such a way that
[TABLE]
for all the values of . Now it is enough to show that
[TABLE]
uniformly in . The latter relation follows immediately from the stationary phase method, since the amplitude function in the integral and all the derivatives in of the amplitude function are uniformly bounded with respect to all the arguments. ∎
3.3 Proof of Lemma 2.2
Recall that
[TABLE]
Let be a constant that may depend on and . Then, by Minkowski’s integral inequality and Lemma 3.2, we have
[TABLE]
since . ∎
3.4 Proof of Lemma 2.4
Let be a constant that may depend on and . Then, applying successively Minkowski’s integral inequality, Holder’s inequality, and Lemma 3.2, we see that
[TABLE]
as . ∎
3.5 Proof of Lemma 2.5
Let and let be a constant that may depend on and . Then the same arguments as in the proof of Lemma 2.4 imply that
[TABLE]
since and (18) holds for . ∎
3.6 Proof of Theorem 2.3
Let us prove (23). We fix . From (25) and Lemmas 2.4 and 2.5, it follows that there exists such that . Other entries of matrix can be treated similarly, i.e.,
[TABLE]
Since , Holder’s inequality implies that
[TABLE]
as Relation (23) is proved.
The stationary phase approximation implies that
[TABLE]
This immediately justifies (24). The last statement of the theorem follows from (22)-(24). ∎
3.7 Proof of Theorem 1.1
Due to Theorem 2.3, one only needs to show that the scattering data for is uniquely determined by the Dirichlet-to-Neumann operator . This will be done by repeating the arguments used in [7, Theorem 4.1] and [10, Theorem 5.1].
Let be two Lipshitz conductivities in such that . Since is Lipschitz continuous, it is differentiable almost everywhere, and the derivatives are bounded [9]. Since and , we have (see [1]). We extend outside in such a way that in and . Let be a bounded domain with a smooth boundary that contains supports of functions . All the previous results will be used below with replaced by and extended as described above. Let be the potential and the solution in (6), the function in (12), and the scattering data in (19) associated with the extended conductivity . Let us note that functions defined by the conductivity problem in are not extensions of the functions defined by the problem in .
Due to equation (20), we have
[TABLE]
Thus it is enough to prove that
[TABLE]
Let be the first column of and , . Since , and equation (2) holds for , it follows that in , and therefore there exists such that
[TABLE]
which is a solution to
[TABLE]
Now we define by
[TABLE]
where is the solution to the Dirichlet problem
[TABLE]
Let . Then
[TABLE]
Hence in . Then
[TABLE]
is the solution of (2) with , and
[TABLE]
is the solution of (6) with .
Lemmas 2.1 and 2.2 imply the unique solvability of the Lippmann-Schwinger equation when and is large enough. Thus, is equal to the first column of when . On the other hand, in coincides with the first column of . Thus the first columns of and are equal on when . Repeating the same steps with the second columns of , we obtain that when , and therefore (28) holds.
The uniqueness of and Theorem 2.3 imply that the potential in the Dirac equation (6) is defined uniquely, and therefore is defined uniquely. Now the conductivity can be found from (5) uniquely up to an additive constant. Finally, this constant can be defined uniquely since is defined uniquely by .
∎
Acknowledgments. The authors are thankful to Daniel Faraco and Keith Rogers for useful discussions.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] G. Alessandrini (1990), Singular Solutions of Elliptic Equations and the Determination of Conductivity by Boundary Measurements. Journal of Differential Equations.
- 2[2] Astala, K., Faraco, D., Rogers, K. M. (2016), Unbounded potential recovery in the plane. Annales scientifiques de l´ École Normale Supérieure. 49, 1023–1047.
- 3[3] Astala K., P a ¨ ¨ 𝑎 \ddot{a} iv a ¨ ¨ 𝑎 \ddot{a} rinta L. (2006). Calderón’s inverse conductivity problem in the plane, Annals of Mathematics, 265-299.
- 4[4] Beals R., Coifman R.R., (1985). Multidimensional inverse scatterings and nonlinear partial differential equations. In F. Treves, editor, Pseudodifferential operators and applications, volume 43 of Proceedings of symposia in pure mathematics, Amer. Math. Soc., 45-70.
- 5[5] 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 (2015).
- 6[6] Borcea L., (2002). Electrical impedance tomography, Inverse Problems 18, 99–136.
- 7[7] Brown R. M., Uhlmann G. A. (1997). Uniqueness in the inverse conductivity problem for nonsmooth conductivities in two dimensions, Communications in partial differential equations, 22(5-6), 1009-1027.
- 8[8] Bukhgeim A. L. (2008). Recovering a potential from Cauchy data in the two-dimensional case, J. Inverse Ill-Posed Probl., 16(1), 19-33.
