A discrete Liouville identity for numerical reconstruction of Schrödinger potentials

We propose a discrete approach for solving an inverse problem for Schrödinger's equation in two dimensions, where the unknown potential is to be determined from boundary measurements of the Dirichlet to Neumann map. For absorptive potentials, and in the continuum, it is known that by using the Liouville identity we obtain an inverse conductivity problem. Its discrete analogue is to find a resistor network that matches the measurements, and is well understood. Here we show how to use a discrete Liouville identity to transform its solution to that of Schrödinger's problem. The discrete Schrödinger potential given by the discrete Liouville identity can be used to reconstruct the potential in the continuum in two ways. First, we can obtain a direct but coarse reconstruction by interpreting the values of the discrete Schrödinger potential as averages of the continuum Schrödinger potential on a special sensitivity grid. Second, the discrete Schrödinger potential may be used to reformulate the conventional nonlinear output least squares optimization formulation of the inverse Schrödinger problem. Instead of minimizing the boundary measurement misfit, we minimize the misfit between the discrete Schrödinger potentials. This results in a better behaved optimization problem that converges in a single Gauss-Newton iteration, and gives good quality reconstructions of the potential, as illustrated by the numerical results.