The Calderón problem for the fractional Schrödinger equation with drift

Abstract

We investigate the Calderón problem for the fractional Schrödinger equation with drift, proving that the unknown drift and potential in a bounded domain can be determined simultaneously and uniquely by an infinite number of exterior measurements. In particular, in contrast to its local analogue, this nonlocal problem does not enjoy a gauge invariance. The uniqueness result is complemented by an associated logarithmic stability estimate under suitable apriori assumptions. Also uniqueness under finitely many generic measurements is discussed. Here the genericity is obtained through singularity theory which might also be interesting in the context of hybrid inverse problems. Combined with the results from Ghosh et al. (Uniqueness and reconstruction for the fractional Calderón problem with a single easurement, 2018. arXiv:1801.04449), this yields a finite measurements constructive reconstruction algorithm for the fractional Calderón problem with drift. The inverse problem is formulated as a partial data type nonlocal problem and it is considered in any dimension $n\ge 1$.

1 Introduction

In this article, we consider an inverse problem for a nonlocal Schrödinger equation with drift. Here we seek to study the uniqueness, stability and reconstruction properties in analogy to its local counterpart, which is a type of magnetic Schrödinger equation and had been investigated in Nakamura et al. [34]. As one of our main results, we prove that: In contrast to its local counterpart, for the fractional Calderón problem with drift there is no gauge invariance present (see Theorem 1.1). In particular, this poses an obstruction in possibly extracting information from the nonlocal inverse problem for its local analogue (as $s\rightarrow 1$). As our second main result, we prove a generic, finite measurements reconstruction, which might also be of interest in the context of hybrid inverse problems (see Theorem 1.4 and the following discussions). As a key tool for this, we rely on singularity theory from Whitney [53].

The classical Calderón problem with drift. Before turning to the nonlocal problem, let us recall its local analogue and the known results on this: In the classical Calderón problem for the magnetic Schrödinger equation, the objective is to determine the drift and potential coefficients simultaneously. More precisely, for $n\ge 3$, let $\Omega \subset \mathbb {R}^n$ be a bounded Lipschitz domain, then consider the following Dirichlet boundary value problem

$$\begin{aligned} \begin{aligned} (-\Delta + b\cdot \nabla +c)u&=0 \text { in }\Omega ,\\ u&=f \text { on }\partial \Omega , \end{aligned} \end{aligned}$$

(1.1)

where b and c are sufficiently smooth functions which vanish on $\partial \Omega$. Assuming the well-posedness of the boundary value problem (1.1), one can define boundary measurements given by the Dirichlet-to-Neumann map (abbreviated as the DN map in the rest of this paper)

$$\begin{aligned} \Lambda _{b,c}:H^{1/2}(\partial \Omega ) \rightarrow H^{-1/2}(\partial \Omega ) \text { with }\Lambda _{b,c}:f \mapsto \dfrac{\partial u}{\partial \nu }, \end{aligned}$$

where $u\in H^1(\Omega )$ is the unique solution to (1.1) and $\nu$ is the unit outer normal on $\partial \Omega$.

The Calderón problem for (1.1) consists of trying to recover the unknown coefficients b, c (which are assumed to be in appropriate function spaces) by the information encoded in the boundary measurement operator $\Lambda _{b,c}$ on $\partial \Omega$. In the case of the local magnetic Schrödinger equation, there is however an intrinsic obstruction to the unique identification of these coefficients. To observe this, consider the substitution $v= e^{\phi } u$ which results in an equation for v which is of a similar form as (1.1):

$$\begin{aligned} \begin{aligned} -\Delta v + (b+ 2\nabla \phi ) \cdot \nabla v + (c+\Delta \phi - b\cdot \nabla \phi -|\nabla \phi |^2) v&= 0 \text{ in } \Omega ,\\ v&= e^{\phi } f \text{ on } \partial \Omega . \end{aligned} \end{aligned}$$

(1.2)

If now $\phi = 0$, $\partial _{\nu } \phi = 0$ on $\partial \Omega$, then the DN maps for the old and the new equations (1.1) and (1.2) coincide; there is a hidden gauge invariance. As a consequence, one cannot expect to be able to recover the full information on the coefficients b and c from the knowledge of the Dirichlet-to-Neumann map $\Lambda _{b,c}$ on $\partial \Omega$.

As a matter of fact, one can at most hope to recover b and c up to the described gauge invariance for the inverse boundary value problem with respect to (1.1). This is indeed the case (c.f. [28, Theorem 5.4.1]): If $b_j$ and $c_j$ are compactly supported in a simply connected domain for $j=1,2$, and if their DN maps coincide on the boundary $\partial \Omega$, then one can obtain uniqueness up to the described gauge invariance:

$$\begin{aligned} \mathrm {curl}b_1=\mathrm {curl}b_2 \text { and }4c_1+b_1\cdot b_1-2\mathrm {div}b_1=4c_2+b_2\cdot b_2-2\mathrm {div}b_2 \text { in }\Omega . \end{aligned}$$

(1.3)

We again emphasize that this does not allow us to recover the full fields b, c but only allows one to obtain information up to the above gauge invariance.

Variants of this inverse boundary value problem have also been studied in [43, 48, 49] for the symmetric magnetic Schrödinger operator, in [25, 30] for the low regularity setting and in [16] for the setting in which only partial data are available. The setting of flexible geometries was studied in [8, 31]. The case of systems was considered in [13] and Yang–Mills potentials with arbitrary geometry in [9]. Stability results can be found in [29, 50] and reconstruction results are given in [44]. For more detailed discussions of inverse problems for the magnetic Schrödinger equation, we refer to the book [28, Chapter 5] as well as to the articles [45, 51] and their bibliographies.

The fractional Calderón problem with drift. Keeping the situation of the local problem with $s=1$ in the back of our minds, we turn to the Calderón problem for the analogous fractional Schrödinger equation with drift, which is a nonlocal inverse problem. This inverse problem should be regarded as a generalization of the fractional Calderón problem, which had first been introduced and investigated in Ghosh et al. [20]. In the sequel, we describe this problem more precisely.

Assume that $\Omega \subset \mathbb {R}^n$ is a bounded Lipschitz domain for $n\ge 1$ and let $\frac{1}{2}<s<1$ (so that the fractional nonlocal operator is dominant). Given a drift $b\in W^{1-s,\infty }(\Omega )^n$ and a potential $c\in L^{\infty }(\Omega )$, we consider the following fractional exterior value problem

$$\begin{aligned} \begin{aligned} ((-\Delta )^s + b \cdot \nabla + c)u&= 0 \text{ in } \Omega ,\\ u&= f \text{ in } \Omega _e:=\mathbb {R}^n {\setminus } \overline{\Omega }, \end{aligned} \end{aligned}$$

(1.4)

with some suitable exterior datum f (we will present a rigorous mathematical formulation of this in Sect. 2, c.f. also Remark 2.10). Here the fractional Laplacian $(-\Delta )^s$ is given by

$$\begin{aligned} (-\Delta )^s u:= \mathcal {F}^{-1}\left\{ |\xi |^{2s}\widehat{u}(\xi )\right\} , \text { for }u\in H^s(\mathbb {R}^n), \end{aligned}$$

where $\widehat{u}=\mathcal {F}u$ denotes the Fourier transform of u. We assume that this problem and its “adjoint problem” are well-posed, which is guaranteed by imposing the eigenvalue condition that

$$\begin{aligned} \begin{aligned}&\text {If } w\in H^s(\mathbb {R}^n)\text { is a solution of }(-\Delta )^sw+b\cdot \nabla w +cw=0\text { in }\Omega \text{ with } w = 0 \text{ in } \Omega _e,\\&\text {then we have }w\equiv 0. \end{aligned} \end{aligned}$$

(1.5)

In the sequel, we will always suppose that this condition is satisfied. As in the classical Calderón problem with drift, we are interested in recovering the drift coefficient $b \in W^{1-s,\infty }(\Omega )^n$ and the potential $c \in L^{\infty }(\Omega )$ simultaneously from the associated DN map. With slight abuse of notation (for a precise definition we refer to Sect. 2), the DN map associated with the nonlocal problem can be thought of as the mapping

$$\begin{aligned} \Lambda _{b,c}: \widetilde{H}^{s}(\Omega _e) \rightarrow (\widetilde{H}^{s}(\Omega _e))^{*}, \ f \mapsto (-\Delta )^s u|_{\Omega _e}, \end{aligned}$$

(1.6)

where u is the solution to (1.4) with exterior data f. Here $\widetilde{H}^{s}(\Omega _e)$ denotes the completion of $C_c^{\infty }(\Omega _e)$ with respect to a fractional Sobolev norm. We refer to (2.1) for the precise definition.

The aim of this work is to prove the global uniqueness and stability of the drift coefficient and the potential for this nonlocal inverse problem. This is in strong contrast to the local case, i.e., the case $s=1$, as in the nonlocal setting the gauge invariance (1.3) which presented an obstruction to global uniqueness in the local case disappears.

1.1 The main results

Let us formulate our main results. As a first property we obtain the global uniqueness of the drift coefficient and the potential for the nonlocal fractional Calderón problem with drift:

Theorem 1.1

(Global uniqueness) For $n\ge 1$, let $\Omega \subset \mathbb {R}^n$ be a bounded, open, non-empty Lipschitz domain, and let $\frac{1}{2}<s<1$. Let $b_j\in W^{1-s,\infty }(\Omega )^n$ be two drift fields, and $c_j\in L^\infty (\Omega )$ be potentials for $j=1,2$. Given arbitrary open, non-empty sets $W_1,W_2 \subset \Omega _e$, suppose that the DN maps for the equations

$$\begin{aligned} ((-\Delta )^s+b_j\cdot \nabla +c_j)u_j=0 \text { in }\Omega \end{aligned}$$

satisfy

$$\begin{aligned} \Lambda _{b_1,c_1}f|_{W_2}=\Lambda _{b_2,c_2}f|_{W_2}, \text { for any }f\in C^\infty _c(W_1). \end{aligned}$$

Then $b_1=b_2$ and $c_1=c_2$ in $\Omega$.

Theorem 1.1 can be regarded as a partial data result for our nonlocal inverse problem. Different from the local case, i.e., $s=1$, there is no gauge invariance and thus no intrinsic obstruction to uniqueness in this nonlocal Calderón problem. We expect that it is possible to improve the regularity assumptions on the drift coefficient and the potential. As our main focus in the present article is however on the striking differences between the local and the nonlocal problems in terms of the existence/absence of a gauge, we do not elaborate on this here but postpone this to a future work.

As in previous results on the fractional Calderón problem (c.f. [20]), in this work, we are not using complex geometrical optics solutions. Instead, we rely on the following approximation property.

Theorem 1.2

(Runge approximation) For $n\ge 1$ and $\frac{1}{2}<s<1$, let $\Omega \subset \mathbb {R}^n$ be a bounded open, non-empty Lipschitz set and $\Omega _1\subset \mathbb {R}^n$ be an arbitrary open set such that $\Omega _1 {\setminus } \overline{\Omega }\ne \emptyset$.

• (a) Let $b\in W^{1-s,\infty }(\Omega )^n$ and $c\in L^\infty (\Omega )$. Then, for any $g\in L^2(\Omega )$ and $\epsilon >0$, one can find a solution $u_\epsilon \in H^s(\mathbb {R}^n)$ of

  $$\begin{aligned} \left( (-\Delta )^s+b\cdot \nabla +c\right) u_\epsilon =0 \text { in }\Omega , \text { with }\mathrm {supp}(u_\epsilon )\subset \overline{\Omega _1} \end{aligned}$$

  such that

  $$\begin{aligned} \Vert u_\epsilon - g \Vert _{L^2(\Omega )}<\epsilon . \end{aligned}$$

• (b) If we further assume that $\Omega$ has a $C^\infty$-smooth boundary, $b\in C^\infty _c(\Omega )^n$ and $c\in C^\infty _c(\Omega )$ with ${{\,\mathrm{supp}\,}}(b), {{\,\mathrm{supp}\,}}(c)\Subset \Omega$, given any $g\in C^\infty (\overline{\Omega })$, $\epsilon >0$ and $k\in \mathbb {N}$, then there exists a solution $u_{\epsilon }\in H^{s}(\mathbb {R}^n)$ of

  $$\begin{aligned} \left( (-\Delta )^s+b\cdot \nabla +c\right) u_\epsilon =0 \text { in }\Omega , \text { with }\mathrm {supp}(u_\epsilon )\subset \overline{\Omega _1} \end{aligned}$$

  such that

  $$\begin{aligned} \Vert d^{-s}(x) u_{\epsilon } -g\Vert _{C^k(\Omega )} \le \epsilon . \end{aligned}$$

  The function d(x) is any $C^\infty$-smooth function defined in $\overline{\Omega }$ such that $d>0$ in $\Omega$ and $d(x)=\mathrm {dist}(x,\partial \Omega )$ whenever x is near $\partial \Omega$.

A more quantitative Runge approximation for the fractional Schrödinger equation with drift will be discussed in Sect. 5 in the context of stability estimates.

Remark 1.1

The qualitative Runge approximation property as a key tool for studying fractional Schrödinger type inverse problems had been introduced in [20]. In order to infer such a result, the authors of [20] built on the unique continuation property for fractional Schrödinger equations in the form of Carleman estimates which had been derived in [35], c.f. also [14, 15, 19, 36, 47, 54] for related unique continuation results for fractional Schrödinger equations. For variable coefficient fractional Schrödinger operators, the authors of [18] utilized Almgren’s frequency function to derive such a property. In the context of nonlocal elliptic equations, this had earlier been employed by [15, 54], c.f. also Section 7 in [35] for the derivation of unique continuation properties with variable coefficients.

Let us put these results into the context of the literature on the fractional Calderón problem: The problem was first introduced by [20], where the authors treated the case with $c\in L^{\infty }(\Omega )$, $b=0$ and $s\in (0,1)$, and proved a global uniqueness result for c. For more general nonlocal variable coefficient Schrödinger operators, the fractional Calderón problem was studied in [18]. The techniques based on Runge approximation are strong enough to deal with the case of semilinear equations [32] and low regularity, almost critical function spaces for the potential [38]. Even single measurement results are possible [19] (c.f. the discussion below). Moreover, these techniques have been extended to other nonlocal problems [7, 39] in a slightly different context. Also, for positive and general potentials, monotonicity inversion formulas have been successfully discovered in [26, 27]. Very recently and independently from our work, uniqueness results have been obtained for equations with nonlocal lower order contributions [4].

In addition to uniqueness, stability is of central importance in inverse problems. Stability results for the fractional Calderón problem were first obtained in [38, 40], where optimal logarithmic stability estimates had been derived (c.f. also [41] for improvements of this if structural apriori conditions like the finiteness of the underlying function space are satisfied). It is possible to extend the logarithmic estimates to the setting of the fractional Calderón problem with drift. Here we obtain the following result:

Theorem 1.3

(Logarithmic stability) Let $s\in (\frac{1}{2},1)$, $\Omega \subset \mathbb {R}^n$, $n\ge 1$ be a bounded open, non-empty smooth domain. Let $W_1, W_2$ be open, non-empty sets with $\overline{W}_1, \overline{W}_2 \subset \Omega _e$. Assume that for some constants $M>0$, $\delta >0$

$$\begin{aligned} \Vert b_j\Vert _{W^{1-s+\delta , \infty }(\Omega )} + \Vert c_j\Vert _{H^s(\Omega )} + \Vert c_j\Vert _{W^{1,n+\delta }(\Omega )} \le M, \end{aligned}$$

and that ${{\,\mathrm{supp}\,}}(b_j), {{\,\mathrm{supp}\,}}(c_j) \Subset \Omega$ for $j=1,2$. Then for some constants $\mu >0$ and $C>0$ which depend on $\Omega , W_1, W_2, n, s, M, \delta$, we have

$$\begin{aligned} \Vert c_1-c_2\Vert _{H^{-s}(\Omega )} + \Vert b_1-b_2\Vert _{H^{-s}(\Omega )} \le C\left| \log (\Vert \Lambda _{b_1,c_1}-\Lambda _{b_2,c_2}\Vert _{*})\right| ^{-\mu }, \end{aligned}$$

if $\Vert \Lambda _{b_1,c_1}-\Lambda _{b_2,c_2}\Vert _{*} \le 1$, where $\Vert A\Vert _{*}:= \sup \{ (A f_1, f_2)_{W_2}: \ f_1 \in \widetilde{H}^s(W_1), \ f_2 \in \widetilde{H}^s(W_2) \}$.

As in [38] this relies on quantitative Runge approximation arguments, which are derived from quantitative unique continuation properties. We adapt the arguments from [38] to infer these results for the fractional Schrödinger equation with drift.

Last but not least, based on the higher order Runge approximation property (Theorem 1.2 (b)), we can deduce finite measurements uniqueness results for the fractional Calderón problem with drift.

Theorem 1.4

(Finite measurements uniqueness) Let $\Omega \subset \mathbb {R}^n$ be a bounded, non-empty domain with a $C^\infty$-smooth boundary. Let $W \subset \Omega _e$ be an open, non-empty smooth set such that $\overline{W}\cap \overline{\Omega } = \emptyset$. Let $s\in (\frac{1}{2},1)$ and assume that $b_j\in C^{\infty }_c(\Omega )^n$, $c_j\in C^{\infty }_c(\Omega )$ satisfy (1.5) with ${{\,\mathrm{supp}\,}}(b_j), {{\,\mathrm{supp}\,}}(c_j)\Subset \Omega$ for $j=1,2$. There exist $n+1$ exterior data $f_1,\ldots ,f_{n+1}\in C^{\infty }_c(W)$ such that if

$$\begin{aligned} \Lambda _{b_1,c_1}(f_l) = \Lambda _{b_2,c_2}(f_l) \text{ for } l \in \{1,\ldots ,n+1\}, \end{aligned}$$

then $b_1 = b_2$ and $c_1=c_2$. Moreover, the set of exterior data $f_1,\ldots ,f_{n+1}$, which satisfies this property forms an open and dense subset in $C^{\infty }_{ c}(W)$.

This is analogous to the single measurement results in [19] for the fractional Schrödinger equation, c.f. also [6] for a single measurement result on the detection of an embedded obstacle. However, compared to [19] a word of caution is needed here: In contrast to the result from [19] it is not possible to work with an arbitrary nontrivial set of measurements $f_1,\ldots ,f_{n+1}$. The data $f_1,\ldots ,f_{n+1}$ have to be chosen appropriately from a set which depends on the unknowns b, c. This is similar to results on hybrid inverse problems, c.f. [1, 2].

While the dependence of the admissible exterior data on the unknown drift field and potential seems like a serious restriction at first sight, we emphasise that by proving that the data $f_1,\ldots ,f_{n+1}$ can be chosen in an open and dense set in $C_{c}^{\infty }(W)$, we show that the set of admissible exterior data is very large: Given a (random) exterior measurement $f_1,\ldots ,f_{n+1} \in C_c^{\infty }(W)$, our result states that an arbitrarily small perturbation of this yields an admissible exterior datum from which we can reconstruct the drift field b and potential c. This might also be of interest in the setting of hybrid inverse problems for which we could not find a statement on an openanddense set of admissible measurements. We plan to address this in future research.

For an overview about the fractional Calderón problem, we refer to the surveys [37, 46].

1.2 Outline of the remaining article

The paper is organized as follows. In Sect. 2, we review the notion of a weak solutions of the fractional Schrödinger equation with drift. With this at hand, we define the DN map rigorously. Section 3 demonstrates the $L^2$-Runge approximation property, which proves Theorem 1.2(a). We will prove the global uniqueness result of Theorem 1.1 in Sect. 4, which shows that the nonlocal Calderón problem does not enjoy a gauge invariance in contrast to its local analogue. In Sect. 5, we also prove the stability result of Theorem 1.3 for the fractional Calderón problem with drift and potential with respect to the associated DN maps. In Sect. 6, in the end of this work, we present the proof of Theorem 1.4, where we prove several points on the reconstruction from finitely many exterior measurements. In addition, we study generic unique determination results via singularity theory in the “Appendix”, which is useful to construct the open and dense subset for the exterior data stated in Theorem 1.4.

2 The fractional Schrödinger equation with drift

In this section, we recall the relevant function spaces, prove the well-posedness of the fractional Schrödinger equation with drift and introduce and derive properties of the Dirichlet-to-Neumann map associated with (1.4).

2.1 Preliminaries

We begin by recalling the relevant fractional Sobolev spaces on (bounded) domains. We define the $L^2$-based fractional Sobolev spaces as follows: for $0<s<1$, we consider the fractional Sobolev spaces $H^s(\mathbb {R}^n)=W^{s,2}(\mathbb {R}^n)$ with the norm

$$\begin{aligned} \Vert u\Vert _{H^s(\mathbb {R}^n)}:=\left\| \mathcal {F}^{-1}\left\{ \langle \xi \rangle ^s \widehat{u} \right\} \right\| _{L^2(\mathbb {R}^n)}, \end{aligned}$$

where $\left\langle \xi \right\rangle =(1+|\xi |^{2})^{\frac{1}{2}}$. Let $\mathcal {O} \subset \mathbb {R}^{n}$ be an arbitrary non-empty open set and $0<s<1$, then we define:

$$\begin{aligned} \begin{aligned} H^{s}(\mathcal {O})&:=\{u|_{\mathcal {O}};\,u\in H^{s}(\mathbb {R}^{n})\},\\ \widetilde{H}^{s}(\mathcal {O})&:=\text {closure of } C_{c}^{\infty }(\mathcal {O})\text { in } {H^{s}(\mathbb {R}^{n})},\\ H_{0}^{s}(\mathcal {O})&:=\text {closure of } C_{c}^{\infty }(\mathcal {O})\text { in } {H^{s}(\mathcal {O})}, \end{aligned} \end{aligned}$$

(2.1)

and

$$\begin{aligned} H_{\overline{\mathcal {O}}}^{s}:=\{u\in H^{s}(\mathbb {R}^{n}); \quad \text { with }\mathrm {supp}(u)\subset \overline{\mathcal O}\}. \end{aligned}$$

The norm of $H^{s}(\mathcal {O})$ is denoted by

$$\begin{aligned} \Vert u\Vert _{H^{s}(\mathcal {O})}:=\inf \left\{ \Vert v\Vert _{H^{s}(\mathbb {R}^{n})};v\in H^{s}(\mathbb {R}^{n}) \text{ and } v|_{\mathcal {O}}=u\right\} . \end{aligned}$$

It is known that $\widetilde{H}^{s}(\mathcal {O})\subseteq H_{0}^{s}(\mathcal {O})$, and that $H_{\overline{\mathcal {O}}}^{s}$ is a closed subspace of $H^{s}(\mathbb {R}^{n})$. Further we have for arbitrary non-empty open sets $\mathcal {O}$

$$\begin{aligned} \left( H^{s}(\mathcal {O})\right) ^{*}=\widetilde{H}^{-s}(\mathcal {O}) \text{ and } \left( \widetilde{H}^{s}(\mathcal {O})\right) ^{*}=H^{-s}(\mathcal {O}). \end{aligned}$$

Remark 2.1

When $\mathcal {O}\subset \mathbb {R}^n$ is a bounded Lipschitz domain, we have that for any $s\in \mathbb {R}$,

$$\begin{aligned} \begin{aligned}&\widetilde{H}^{s}(\mathcal {O})=H_{\overline{\mathcal {O}}}^{s}\subsetneq H_{0}^{s}(\mathcal {O}). \end{aligned} \end{aligned}$$

If $s>-\frac{1}{2}$ and $s\notin \{\frac{1}{2},\frac{3}{2},\ldots \}$ the last inclusion also becomes an equality.

We further denote the homogeneous fractional Sobolev spaces as $\dot{H}^s(\mathbb {R}^n)$, where

$$\begin{aligned} \dot{H}^s(\mathbb {R}^n):=\{u: \mathbb {R}^n \rightarrow \mathbb {R}; \ \Vert {{\,\mathrm{\mathcal {F}}\,}}^{-1} \{ |\xi |^s \hat{u} \} \Vert _{L^2(\mathbb {R}^n)} <\infty \}. \end{aligned}$$

We define the associated semi-norm as

$$\begin{aligned} \Vert u\Vert _{\dot{H}^s(\mathbb {R}^n)}:=\Vert {{\,\mathrm{\mathcal {F}}\,}}^{-1} \{ |\xi |^s \hat{u} \} \Vert _{L^2(\mathbb {R}^n)}. \end{aligned}$$

Note that the norm $\Vert \cdot \Vert _{H^s(\mathbb {R}^n)}$ is equivalent to the norm $\Vert \cdot \Vert _{L^2(\mathbb {R}^n)}+\Vert \cdot \Vert _{\dot{H}^s(\mathbb {R}^n)}$. For a more detailed introduction to fractional Sobolev spaces and related results, we refer the readers to [12, 33].

Since we will use this for our drift fields, we also recall the $L^p$ based fractional Sobolev spaces: we set $\Vert u\Vert _{W^{s,p}(\mathbb {R}^n)} := \Vert \langle D \rangle ^s u \Vert _{L^p(\mathbb {R}^n)}$, where $\langle \xi \rangle = (1+|\xi |^2)^{1/2}$ and $m(D) u = \mathcal {F}^{-1}\{m(\xi ) \widehat{u}(\xi )\}$ for $m\in C^{\infty }(\mathbb {R}^n)$ such that m and all its derivatives are polynomially bounded, and u is a tempered distribution. For a non-empty open set $\mathcal {O}\subset \mathbb {R}^n$ and $p>1$, we then define the space $W^{s,p}(\mathcal {O})$ by

$$\begin{aligned} W^{s,p}(\mathcal {O}) = \{u|_{\mathcal {O}}; \ u \in W^{s,p}(\mathbb {R}^n)\}. \end{aligned}$$

This is equipped with the associated norm

$$\begin{aligned} \Vert u\Vert _{W^{s,p}(\mathcal {O})} = \inf \{\Vert w\Vert _{W^{s,p}(\mathbb {R}^n)}; \ w \in W^{s,p}(\mathbb {R}^n), \ w|_{\mathcal {O}} = u\}. \end{aligned}$$

We also define

$$\begin{aligned} W^{s,p}_0(\mathcal {O}):= \text{ closure } \text{ of } C_c^{\infty }(\mathcal {O}) \text{ in } W^{s,p}(\mathcal {O}). \end{aligned}$$

In the sequel, we will only use these more general $W^{s,p}$ function spaces to quantify the size of the drift field b.

We conclude this section by recalling a fractional Poincaré type inequality:

Lemma 2.2

Let $n\ge 1$ and let $s\in (0,1)$. Assume that $\Omega \subset \mathbb {R}^n$ is non-empty, open and bounded. Then, there exists a constant $C>0$ such that

$$\begin{aligned} \Vert v\Vert _{L^2(\Omega )} \le C ({{\,\mathrm{diam}\,}}(\Omega ))^{s} \Vert (-\Delta )^{s/2} v\Vert _{L^2(\mathbb {R}^n)} \text{ for } v \in \widetilde{H}^s(\Omega ), \end{aligned}$$

where ${{\,\mathrm{diam}\,}}(\Omega )$ denotes the diameter of $\Omega$.

We present the proof for self-containedness but follow the idea from the appendix in [42].

Proof of Lemma 2.2

We first assume that $v\in C_c^{\infty }(\Omega )$. The result will then follow by density of $C_c^{\infty }(\Omega )$ in $\widetilde{H}^s(\Omega )$. Let $x\in \Omega$ be arbitrary and let $x' = x + 2 {{\,\mathrm{diam}\,}}(\Omega ) \frac{x}{|x|}$. Let further $\widetilde{v}$ be the Caffarelli-Silvestre [5] extension of v, i.e. let $\widetilde{v}$ be the solution of

$$\begin{aligned} \nabla \cdot x_{n+1}^{1-2s} \nabla \widetilde{v}&= 0 \text{ in } \mathbb {R}^{n+1}_+,\\ \widetilde{v}&= v \text{ on } \mathbb {R}^n \times \{0\}. \end{aligned}$$

Then the fundamental theorem of calculus and the support condition for v imply that for any $r\in ({{\,\mathrm{diam}\,}}(\Omega ),\infty )$

$$\begin{aligned} |v(x)|&= |\widetilde{v}(x,0)| \le |\widetilde{v}(x,r)| + \int \limits _{0}^{r} |\partial _{n+1} \widetilde{v}(x,t)| dt\\&\le |\widetilde{v}(x',r)| + \int \limits _{0}^{1}|\nabla ' \widetilde{v}(tx + (1-t)x',r) | |x-x'| dt + \int \limits _{0}^{r} |\partial _{n+1} \widetilde{v}(x,t)| dt\\&\le \int \limits _{0}^r |\partial _{n+1} \widetilde{v}(x',t)|+ \int \limits _{0}^{r} |\partial _{n+1} \widetilde{v}(x,t)| dt + \int \limits _{0}^{1}|\nabla ' \widetilde{v}(tx + (1-t)x',r) | |x-x'| dt. \end{aligned}$$

Next we insert the weights $t^{\frac{1-2s}{2}}$, apply Hölder’s inequality, estimate $|x'-x| \le Cr$ and fix $r \in (r_0, 2r_0)$, where $r_0 = {{\,\mathrm{diam}\,}}(\Omega )$:

$$\begin{aligned} |v(x)|&\le \int \limits _{0}^r t^{\frac{2s-1}{2}}t^{\frac{1-2s}{2}}|\partial _{n+1} \widetilde{v}(x',t)| dt + \int \limits _{0}^{r} t^{\frac{2s-1}{2}}t^{\frac{1-2s}{2}} |\partial _{n+1} \widetilde{v}(x,t)| dt \\&\quad + \int \limits _{0}^{1}r^{\frac{1-2s}{2}} |\nabla ' \widetilde{v}(tx + (1-t)x',r) | |x-x'| r^{\frac{2s-1}{2}} dt\\&\le Cr^s\left( \left( \int \limits _{0}^{r} x_{n+1}^{1-2s} |\partial _{n+1} \widetilde{v}(x',x_{n+1})|^2 dx_{n+1} \right) ^{\frac{1}{2}} + \left( \int \limits _{0}^r x_{n+1}^{1-2s} |\partial _{n+1} \widetilde{v}(x,x_{n+1})|^2 dx_{n+1} \right) ^{\frac{1}{2}} \right. \\&\quad \left. + r^{\frac{1}{2}}\left( \int \limits _{0}^1 r^{1-2s} |\nabla ' \widetilde{v}(t x + (1-t)x',r)|^2 dt \right) ^{\frac{1}{2}} \right) \\&\le Cr^s_0\left( \left( \int \limits _{0}^{2 r_0} x_{n+1}^{1-2s} |\partial _{n+1} \widetilde{v}(x',x_{n+1})|^2 dx_{n+1} \right) ^{\frac{1}{2}} + \left( \int \limits _{0}^{2 r_0} x_{n+1}^{1-2s} |\partial _{n+1} \widetilde{v}(x,x_{n+1})|^2 dx_{n+1} \right) ^{\frac{1}{2}} \right. \\&\quad \left. + r^{\frac{1}{2}}_0\left( \int \limits _{0}^1 r^{1-2s} |\nabla ' \widetilde{v}(t x + (1-t)x',r)|^2 dt \right) ^{\frac{1}{2}} \right) . \end{aligned}$$

Here the constant $C>0$ in particular depends on s. Next we square the estimate and integrate it in the normal direction in the interval $r\in (r_0/2, 2r_0)$. This yields

$$\begin{aligned} |v(x)|^2&\le Cr^{2s}_0\left( \int \limits _{0}^{2 r_0} x_{n+1}^{1-2s} |\partial _{n+1} \widetilde{v}(x',x_{n+1})|^2 dx_{n+1} + \int \limits _{0}^{2 r_0} x_{n+1}^{1-2s} |\partial _{n+1} \widetilde{v}(x,x_{n+1})|^2 dx_{n+1} \right. \\&\quad \left. + \int \limits _{0}^1 \int \limits _{0}^{2r_0} x_{n+1}^{1-2s} |\nabla ' \widetilde{v}(t x + (1-t)x',x_{n+1})|^2 dt dx_{n+1} \right) . \end{aligned}$$

Finally, integrating in tangential directions and using the support condition for v, we obtain

$$\begin{aligned} \Vert v\Vert _{L^2(\Omega )}^2 \le C r_0^{2s} \Vert x_{n+1}^{\frac{1-2s}{2}} \nabla \widetilde{v}\Vert _{L^2(\mathbb {R}^{n+1}_+)}^2. \end{aligned}$$

Since by the work of Caffarelli-Silvestre [5] we have

$$\begin{aligned} \Vert x_{n+1}^{\frac{1-2s}{2}} \nabla \widetilde{v}\Vert _{L^2(\mathbb {R}^{n+1}_+)}^2 = \Vert v\Vert _{\dot{H}^{s}(\mathbb {R}^n)}, \end{aligned}$$

and since for $v\in C_c^{\infty }(\Omega )$ we have (in the sense of norm equivalences)

$$\begin{aligned} \Vert (-\Delta )^{s/2} v\Vert _{L^2(\mathbb {R}^n)} \sim \Vert v\Vert _{\dot{H}^{s}(\Omega )} , \end{aligned}$$

this concludes the proof. $\square$

2.2 Well-posedness

In this section we discuss the well-posedness of the Eq. (1.4) and its dual equation.

To this end, let $\Omega \subset \mathbb {R}^n$ be a non-empty bounded domain (open and connected), $b\in W^{1-s,\infty }(\Omega )^n$ be a drift coefficient, $c\in L^\infty (\Omega )$ a potential and let $\frac{1}{2}<s<1$ be a constant. For $F\in (\widetilde{H}^s(\Omega ))^{*}$ (the dual space of $H^s(\Omega )$), $f\in H^{s}(\mathbb {R}^n)$, let us consider the following Dirichlet problem

$$\begin{aligned} \begin{aligned} ((-\Delta )^s+b\cdot \nabla +c)u=F&\text { in }\Omega ,\\ u-f&\in \widetilde{H}^s(\Omega ). \end{aligned} \end{aligned}$$

(2.2)

Given an arbitrary non-empty open set $\mathcal O\subset \mathbb {R}^n$ and $v,w\in L^2(\mathcal O)$, we use the notation

$$\begin{aligned} (v,w)_{\mathcal O}:=\int _{\mathcal O}vw dx. \end{aligned}$$

For $v,w \in C^{\infty }_c(\mathbb {R}^n)$ we define the bilinear form $B_{b,c}(\cdot ,\cdot )$ by

$$\begin{aligned} B_{b,c}(v,w):=((-\Delta )^{s/2}v,(-\Delta )^{s/2}w)_{\mathbb {R}^n}+(b\cdot \nabla v,w)_{\Omega }+(cv,w)_\Omega . \end{aligned}$$

(2.3)

Notice that the bilinear form $B_{b,c}(\cdot , \cdot )$ is not symmetric, so we also introduce the adjoint bilinear form as

$$\begin{aligned} B_{b,c}^{*}(v^{*},w^{*}):=((-\Delta )^{s/2}v^{*},(-\Delta )^{s/2}w^{*})_{\mathbb {R}^n} +(b v^{*}, \nabla w^{*} )_\Omega +(c v^{*},w^{*})_\Omega , \end{aligned}$$

(2.4)

for $v^{*},w^{*} \in C^{\infty }_c(\mathbb {R}^n)$.

We remark that the term “adjoint” is used with a slight abuse of notation here, e.g. as we did not specify the underlying function spaces. We however think of the adjoint bilinear form (2.4), as the bilinear form associated with the adjoint exterior value problem

$$\begin{aligned} \begin{aligned} (-\Delta )^s u^{*} - \nabla \cdot (b u^{*})+cu^{*}=F^{*}&\text { in }\Omega , \\ u^{*} - f^{*}&\in \widetilde{H}^s(\Omega ), \end{aligned} \end{aligned}$$

(2.5)

for some suitable source $F^{*}$ and exterior datum $f^{*}$.

Remark 2.3

We further stress that there is a slight discrepancy between the bilinear form (2.4) and the adjoint Dirichlet problem (2.5) in that we have ignored the boundary contribution originating from the integration by parts of (2.5) in the definition of the bilinear form (2.3). In the sequel, this will for instance be reflected in the (symmetry) properties of the operator $\Lambda _{b,c}^{*}$, c.f. Lemma 2.13. More precisely, the DN map associated with the Eq. (2.5) would contain a boundary contribution on $\partial \Omega$ if b (or $u^{*}$) does not vanish there. This is a consequence of an integration by parts which is used to obtain the weak form of the DN map associated with the Eq. (2.5). For the operator $\Lambda _{b,c}^{*}$ which is defined through the “adjoint bilinear form” (2.4) this contribution is not present, as the integration by parts has already been carried out.

Before continuing in our discussion, we observe that under our regularity assumptions the bilinear forms (2.3), (2.4) are well-defined on $H^{s}(\mathbb {R}^n)$, for $\frac{1}{2}<s<1$.

Lemma 2.4

(Boundedness of bilinear forms) For $n\ge 1$, let $\Omega \subset \mathbb {R}^n$ be a bounded Lipschitz domain and $\frac{1}{2}<s<1$. Assume that $b\in W^{1-s,\infty }(\Omega )^n$, $c\in L^{\infty }(\Omega )$. Let $B_{b,c}(\cdot , \cdot )$ and $B^{*}_{b,c}(\cdot ,\cdot )$ be the bilinear and the adjoint bilinear form defined by (2.3) and (2.4), respectively. Then, $B_{b,c}(\cdot ,\cdot )$ and $B_{b,c}^{*}(\cdot ,\cdot )$ extend as bounded bilinear forms on $H^s(\mathbb {R}^n)\times H^{s}(\mathbb {R}^n)$.

Proof

It suffices to discuss the extension of $B_{b,c}(\cdot ,\cdot )$ as the argument for $B_{b,c}^{*}(\cdot ,\cdot )$ is analogous. First, we directly have

$$\begin{aligned}&\left| ((-\Delta )^{s/2}v,(-\Delta )^{s/2}w)_{\mathbb {R}^n}\right| +\left| (cv,w)_\Omega \right| \nonumber \\&\quad \le C(1+\Vert c\Vert _{L^{\infty }(\Omega )}) \left( \Vert v\Vert _{\dot{H}^{s}(\mathbb {R}^n)}\Vert w\Vert _{\dot{H}^{s}(\mathbb {R}^n)} + \Vert v\Vert _{L^2(\Omega )}\Vert w\Vert _{L^2(\Omega )}\right) \nonumber \\&\quad \le C(1+\Vert c\Vert _{L^{\infty }(\Omega )}) \Vert v\Vert _{H^{s}(\mathbb {R}^n)}\Vert w\Vert _{H^{s}(\mathbb {R}^n)}, \end{aligned}$$

(2.6)

for some constant $C>0$ independent of v, w. It remains to discuss the contribution of the drift term. For this we note that as $1-s \in (0,\frac{1}{2})$ and as $\Omega$ is a Lipschitz domain, we have $\widetilde{H}^{1-s}(\Omega ) = H^{1-s}(\Omega )$ and $\Vert w\Vert _{\widetilde{H}^{1-s}(\Omega )} \le C \Vert w\Vert _{H^{1-s}(\mathbb {R}^n)}$. Next, we choose $B \in W^{1-s,\infty }(\mathbb {R}^n)$ such that

$$\begin{aligned} B|_{\Omega } = b, \ \Vert B\Vert _{W^{1-s,\infty }(\mathbb {R}^n)} \le 2\Vert b\Vert _{W^{1-s,\infty }(\Omega )}. \end{aligned}$$

and estimate

$$\begin{aligned} \begin{aligned} \left| \int _\Omega w (b\cdot \nabla v) dx\right|&\le \Vert b w\Vert _{\widetilde{H}^{1-s}(\Omega )} \Vert \nabla v\Vert _{H^{s-1}(\Omega )} \\&\le C \Vert B w \Vert _{H^{1-s}(\mathbb {R}^n)} \Vert \nabla v\Vert _{H^{s-1}(\mathbb {R}^n)}\\&\le C \Vert w \Vert _{H^{1-s}(\mathbb {R}^n)} \Vert b \Vert _{W^{1-s,\infty }(\Omega )} \Vert \nabla v\Vert _{H^{s-1}(\mathbb {R}^n)} \\&\le C \Vert w \Vert _{H^{1-s}(\mathbb {R}^n)} \Vert b \Vert _{W^{1-s,\infty }(\Omega )} \Vert v\Vert _{H^{s}(\mathbb {R}^n)}\\&\le C \Vert w \Vert _{H^{s}(\mathbb {R}^n)} \Vert b \Vert _{W^{1-s,\infty }(\Omega )} \Vert v\Vert _{H^{s}(\mathbb {R}^n)}. \end{aligned} \end{aligned}$$

(2.7)

Here we used the Kato-Ponce inequality [22] in order to obtain a suitable multiplier estimate

$$\begin{aligned} \begin{aligned} \Vert B w\Vert _{H^{1-s}(\mathbb {R}^n)}&\le C \Vert J^{1-s}(B w)\Vert _{L^2(\mathbb {R}^n)} \\&\le C \left( \Vert B\Vert _{L^{\infty }(\mathbb {R}^n)} \Vert J^{1-s} w \Vert _{L^2(\mathbb {R}^n)} + \Vert J^{1-s}B \Vert _{L^{\infty }(\mathbb {R}^n)} \Vert w \Vert _{L^2(\mathbb {R}^n)}\right) \\&\le C \Vert B\Vert _{W^{1-s,\infty }(\mathbb {R}^n)} \Vert w\Vert _{H^{1-s}(\mathbb {R}^n)}\\&\le 2C \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \Vert w\Vert _{H^{1-s}(\mathbb {R}^n)}\\&\le 2C \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \Vert w\Vert _{H^{s}(\mathbb {R}^n)}, \end{aligned} \end{aligned}$$

(2.8)

where $J^{1-s}:=(\Delta -1)^{\frac{1-s}{2}}$ and $0<1-s< \frac{1}{2}$. Finally, combining (2.6), (2.7) and (2.8), one has the desired result

$$\begin{aligned} \left| B_{b,c}(v,w)\right| \le C\Vert v\Vert _{H^s(\mathbb {R}^n)}\Vert w\Vert _{H^s(\mathbb {R}^n)}, \end{aligned}$$

(2.9)

for $\frac{1}{2}<s<1$ and for some constant $C>0$ independent of v, w. A similar argument also yields the boundedness for the adjoint bilinear form, which reads

$$\begin{aligned} \left| B_{b,c}^*(v^*,w^*)\right| \le C\Vert v^*\Vert _{H^s(\mathbb {R}^n)}\Vert w^*\Vert _{H^s(\mathbb {R}^n)}. \end{aligned}$$

(2.10)

$\square$

Remark 2.5

For the above well-definedness argument, we point out the following observations:

• (a) The requirement that $\Omega$ is a bounded Lipschitz domain only entered in the identity $\widetilde{H}^{1-s}(\Omega ) = H^{1-s}(\Omega )$. If this identity holds for a less regular domain $\Omega$ (e.g. if $\partial \Omega$ is Hölder continuous with sufficiently high Hölder exponent depending on s) the well-definedness of the bilinear forms (2.3), (2.4) persists for this domain. In order to avoid technicalities, we do not address this issue in the sequel, but will always assume that $\Omega$ is Lipschitz. Nevertheless, we phrase our results (e.g. the Dirichlet-to-Neumann map in Definition 2.11) such that they remain valid for a more general class of domains which satisfy the condition $\widetilde{H}^{1-s}(\Omega ) = H^{1-s}(\Omega )$.

• (b) As an alternative condition to imposing Lipschitz regularity, we might also have asked for $b\in W^{1-s,\infty }_0(\Omega )^n$. As the drift field b however is one of the main objects of interest in our argument, we opted for rather imposing regularity on $\Omega$ than on restricting the class of admissible drift fields.

We in particular notice that by the above definitions, (2.3), (2.4) and the estimates (2.9) and (2.10), we obtain

$$\begin{aligned} B_{b,c}(u,w) = B_{b,c}^{*}(w,u), \text { for any }u,w\in H^s(\mathbb {R}^n). \end{aligned}$$

(2.11)

With these bilinear forms at hand, we define $u\in H^s(\mathbb {R}^n)$ with $u-f\in \widetilde{H}^s(\Omega )$ to be a solution of (2.2) if for all $w\in \widetilde{H}^{s}(\mathbb {R}^n)$ we have $B_{b,c}(u,w)=F(w)$. A solution to the adjoint problem is defined analogously using $B^{*}_{b,c}(\cdot ,\cdot )$.

Relying on energy estimates, we can prove the existence and uniqueness of solutions to (2.2) and (2.5), outside of a discrete set of eigenvalues:

Proposition 2.6

Let $\Omega \subset \mathbb {R}^n$, $n\ge 1$, be a bounded Lipschitz domain and $\frac{1}{2}<s<1$. Assume that $b\in W^{1-s,\infty }(\Omega )^n$, $c\in L^{\infty }(\Omega )$. Let $B_{b,c}(\cdot , \cdot )$ and $B^{*}_{b,c}(\cdot ,\cdot )$ be the bilinear and the adjoint bilinear form defined by (2.3) and (2.4), respectively. Then the following properties hold:

• (a) There exists a countable set $\Sigma = \{ \lambda _j \}_{j=1}^{\infty } \subset \mathbb {R}$ with $\lambda _1 \le \lambda _2 \le \cdots \rightarrow \infty$ such that if $\lambda \in \mathbb {R}{\setminus } \Sigma$, for any $f\in H^{s}(\mathbb {R}^n)$ and $F\in (\widetilde{H}^s(\Omega ))^{*}$ there exists a unique function $u\in H^{s}(\mathbb {R}^n)$ such that

  $$\begin{aligned} B_{b,c}(u,w) - \lambda (u,w)_{\Omega } = F(w) \text{ for } \text{ all } w \in \widetilde{H}^s(\Omega ) \text{ and } \ u-f \in \widetilde{H}^{s}(\Omega ). \end{aligned}$$

  We then have

  $$\begin{aligned} \Vert u\Vert _{H^s(\mathbb {R}^n)} \le C(\Vert f\Vert _{H^{s}(\mathbb {R}^n)} + \Vert F\Vert _{(H^{s}(\mathbb {R}^n))^{*}}), \end{aligned}$$

  where the constant $C>0$ depends on $n,s, \Vert b\Vert _{W^{1-s,\infty }(\Omega )}, \Vert c\Vert _{L^{\infty }(\Omega )}, \lambda$. Moreover, $\Sigma \subset (-\Vert c\Vert _{L^{\infty }(\Omega )}-C \Vert b\Vert _{W^{1-s,\infty }(\Omega )}, \infty )$ where the constant $C>0$ is sufficiently large and only depends on s, n.

• (b) Similarly, for the adjoint bilinear form (2.4), there exists a countable set $\Sigma ^{*}\subset \mathbb {R}$ such that if $\lambda ^{*}\in \mathbb {R}{\setminus } \Sigma ^{*}$, for all $f^* \in H^s(\mathbb {R}^n)$ there is a unique $u^* \in H^s(\mathbb {R}^n)$ with

  $$\begin{aligned} B_{b,c}^{*}(u^{*},w^{*}) - \lambda ^{*}(u^{*},w^{*})_{\Omega } = F^{*}(w^{*}), \text { for all }w^{*}\in \widetilde{H}^s(\Omega ) \text { and } u^{*}-f^{*} \in \widetilde{H}^{s}(\Omega ). \end{aligned}$$

Remark 2.7

We point out that the regularity of $b \in W^{1-s,\infty }(\Omega )^n$ is determined by the properties of multipliers in the function space $H^{1-s}(\Omega )$ for $s\in (\frac{1}{2},1)$ (through the application of the Kato-Ponce inequality, c.f. the proof of the well-posedness result below and [22]). This is consistent with the fact that for $s=1$ it is easily seen that $b\in L^{\infty }(\Omega )$ is a sufficient condition for the well-posedness of the problem from Proposition 2.6.

Remark 2.8

For $\Omega$ Lipschitz, we recall that $\widetilde{H}^s(\Omega ) = H^{s}_{\overline{\Omega }}$. In particular the above well-posedness theory then also transfers to these functions spaces immediately.

Remark 2.9

We emphasize that we did not attempt to optimize the function spaces for the drift and the potential contributions here. It is an interesting question which we plan to investigate in future work under which regularity conditions the properties of the inverse problem under investigation persist.

Proof of Proposition 2.6

We first give the proof of (a); it will follow from the spectral theorem for compact operators. Indeed, we first assume that $f=0$, which can always be achieved by considering $v=u-f\in \widetilde{H}^s(\Omega )$. Arguing similarly as in (2.8) for the resulting drift term, we note that this yields an appropriately modified functional $\widetilde{F} \in (\widetilde{H}^{s}(\Omega ))^{*}$. In the sequel, we only consider the function v. Since the boundedness of the bilinear form was already proven by Lemma 2.4, we only need to prove the coercivity of $B_{b,c}(\cdot ,\cdot )$.

To this end, we first note that for $\phi \in C_c^{\infty }(\mathbb {R}^n)$ the following interpolation estimate holds

$$\begin{aligned} \Vert \phi \Vert _{H^{1-s}(\mathbb {R}^n)} \le C \Vert \phi \Vert _{L^2(\mathbb {R}^n)}^{\frac{2s-1}{s}} \Vert \phi \Vert _{H^{s}(\mathbb {R}^n)}^{\frac{1-s}{s}} \text{ for } \frac{1}{2}<s<1. \end{aligned}$$

(2.12)

By virtue of Remark 2.1 and the interpolation estimate (2.12), one has for $\phi , v \in C_c^{\infty }(\mathbb {R}^n)$,

$$\begin{aligned} \begin{aligned} \left| \int _\Omega \phi (\nabla v\cdot b) dx\right|&\le \Vert b \phi \Vert _{\widetilde{H}^{1-s}(\Omega )} \Vert \nabla v\Vert _{H^{s-1}(\Omega )} \\&\le C \Vert b \phi \Vert _{H^{1-s}(\mathbb {R}^n)} \Vert \nabla v\Vert _{H^{s-1}(\mathbb {R}^n)}\\&\le C \Vert \phi \Vert _{H^{1-s}(\mathbb {R}^n)} \Vert b \Vert _{W^{1-s,\infty }(\Omega )} \Vert \nabla v\Vert _{H^{s-1}(\mathbb {R}^n)} \\&\le C \Vert \phi \Vert _{H^{1-s}(\mathbb {R}^n)} \Vert b \Vert _{W^{1-s,\infty }(\Omega )} \Vert v\Vert _{H^{s}(\mathbb {R}^n)}\\&\le C \Vert v \Vert _{H^{s}(\mathbb {R}^n)} \Vert b \Vert _{W^{1-s,\infty }(\Omega )} \Vert \phi \Vert _{L^2(\mathbb {R}^n)}^{\frac{2s-1}{s}} \Vert \phi \Vert _{H^{s}(\mathbb {R}^n)}^{\frac{1-s}{s}}\\&\le C_{\epsilon } \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \Vert v \Vert _{H^{s}(\mathbb {R}^n)} \Vert \phi \Vert _{L^{2}(\mathbb {R}^n)} + \epsilon \Vert v \Vert _{H^{s}(\mathbb {R}^n)} \Vert \phi \Vert _{H^{s}(\mathbb {R}^n)}, \end{aligned} \end{aligned}$$

(2.13)

which holds for any $\epsilon >0$ and for (generic) constants $C,C_\epsilon >0$ which are independent of v and $\phi$. Here we have utilized Young’s inequality in the last line; and used the Kato-Ponce inequality as in (2.7). Notice that by the density of $C_c^{\infty }(\Omega )$ in $\widetilde{H}^s(\Omega )$, the estimate (2.13) also extends to $v, \phi \in \widetilde{H}^s(\Omega )$. A matching coercivity estimate follows from the Poincaré inequality (c.f. Lemma 2.2)

$$\begin{aligned} \Vert v\Vert _{L^2(\mathbb {R}^n)} \le C \Vert (-\Delta )^{s/2} v\Vert _{L^2(\mathbb {R}^n)}, \ v \in \widetilde{H}^s(\Omega ), \end{aligned}$$

in combination with (2.13) and Young’s inequality:

$$\begin{aligned} B_{b,c}(v,v)&\ge \Vert (-\Delta )^{s/2} v\Vert _{L^2(\mathbb {R}^n)}^2 - \Vert c\Vert _{L^{\infty }(\mathbb {R}^n)}\Vert v\Vert _{L^2(\mathbb {R}^n)}^2 - \left| \int \limits _{\Omega } v (\nabla v \cdot b) dx \right| \\&\ge \Vert (-\Delta )^{s/2} v\Vert _{L^2(\mathbb {R}^n)}^2 - \Vert c\Vert _{L^{\infty }(\mathbb {R}^n)}\Vert v\Vert _{L^2(\mathbb {R}^n)}^2 \\&\quad - C_{\epsilon } \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \Vert v\Vert _{L^2(\mathbb {R}^n)}^2 - \epsilon \Vert v\Vert _{H^s(\mathbb {R}^n)}^2 \\&\ge c_0(\Vert (-\Delta )^{s/2} v\Vert _{L^2(\mathbb {R}^n)}^2 + \Vert v\Vert _{L^2(\mathbb {R}^n)}^2) - \Vert c\Vert _{L^{\infty }(\mathbb {R}^n)}\Vert v\Vert _{L^2(\mathbb {R}^n)}^2 \\&\quad - C_{\epsilon } \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \Vert v\Vert _{L^2(\mathbb {R}^n)}^2 - \epsilon \Vert v\Vert _{H^s(\mathbb {R}^n)}^2 \\&\ge \frac{c_0}{2}\Vert v\Vert _{H^{s}(\mathbb {R}^n)}^2 - \Vert c\Vert _{L^{\infty }(\mathbb {R}^n)}\Vert v\Vert _{L^2(\mathbb {R}^n)}^2 - C_{c_0} \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \Vert v\Vert _{L^2(\mathbb {R}^n)}^2, \end{aligned}$$

for some constant $c_0>0$ and $C>0$ independent of v. Here we have chosen $\epsilon = \frac{c_0}{2}$ in the above calculation.

As a consequence, for $\mu = C_{c_0}\Vert b\Vert _{W^{1-s,\infty }(\Omega )} + \Vert c\Vert _{L^{\infty }(\Omega )} \ge 0$ we have

$$\begin{aligned} B_{b,c}(v,v) + \mu (v,v)_{L^2(\mathbb {R}^n)} \ge \frac{c_0}{2} \Vert v\Vert _{H^s(\mathbb {R}^n)}^2, \text { for any }v\in \widetilde{H}^s(\Omega ). \end{aligned}$$

Thus, for $\mu$ as above, $B_{b,c}(\cdot , \cdot ) + \mu (\cdot ,\cdot )_{L^2(\Omega )}$ is a scalar product on $\widetilde{H}^s(\Omega )$. Applying the Riesz representation theorem, we then infer that for all $\widetilde{F}\in (\widetilde{H}^s(\Omega ))^{*}$ there exists $v\in \widetilde{H}^s(\Omega )$ such that

$$\begin{aligned} B_{b,c}(v,\phi ) + \mu (v,\phi )_{L^2(\Omega )}&=\widetilde{F}(\phi ) \text{ for } \text{ all } \phi \in \widetilde{H}^s(\Omega ). \end{aligned}$$

In fact, we may write $v = G_{\mu }(\widetilde{F})$, where $G_{\mu }$ is a bounded, linear operator from $(\widetilde{H}^s(\Omega ))^{*} \rightarrow \widetilde{H}^s(\Omega )$. By the compact Sobolev embedding of $\widetilde{H}^s(\Omega ) \hookrightarrow L^2(\Omega )$ the operator $G_{\mu }$ is a compact operator, if interpreted as a map from $L^2(\Omega )$ to $L^2(\Omega )$.

Since the equation

$$\begin{aligned} B_{b,c}(v,\cdot ) - \lambda (v,\cdot ) = \widetilde{F}(\cdot ) \text{ on } \widetilde{H}^s(\Omega ) \end{aligned}$$

(2.14)

is equivalent to the equation $v= G_{\mu }((\mu +\lambda )v + \widetilde{F})$, we may now invoke the spectral theorem for compact operators to infer that there exists a countable, decreasing sequence $\{\gamma _j\}_{j=1}^{\infty } \subset [0,\infty )$ with $\gamma _j \rightarrow 0$ such that if $\gamma \notin \{\gamma _j\}_{j=1}^{\infty }$ the operator $G-\gamma Id$ is invertible. In particular, we may rewrite $\gamma _j = (\lambda _j+\mu )^{-1}$ for a countable, increasing sequence with $\lambda _j \rightarrow \infty$. Recalling that (2.14) can be written as a corresponding operator equation, then concludes the proof of the solvability result for $B_{b,c}(\cdot ,\cdot )$.

For (b), the proof is similar to (a). As in the previous arguments, one only needs to check the boundedness of $\left| \int _{\Omega } bv^{*}\cdot \nabla \phi ^{*}dx \right|$, where $v^{*}=u^{*}-f^{*}\in \widetilde{H}^s(\Omega )$ and $\phi ^{*}\in \widetilde{H}^s(\Omega )$ is an arbitrary test function. Here $u^{*}$ and $f^{*}$ are the functions which appeared in the Eq. (2.5). In particular, by using the Kato-Ponce inequality, the interpolation inequality (2.12) and Young’s inequality again, for any $\epsilon >0$, we have

$$\begin{aligned} \left| \int _\Omega bv^{*}\cdot \nabla \phi ^{*} dx\right| \le C_{\epsilon } \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \Vert \phi ^{*}\Vert _{H^{s}(\mathbb {R}^n)} \Vert v^{*}\Vert _{L^{2}(\mathbb {R}^n)} + \epsilon \Vert \phi ^{*}\Vert _{H^{s}(\mathbb {R}^n)} \Vert v^{*}\Vert _{H^{s}(\mathbb {R}^n)}, \end{aligned}$$

for some positive constants $C_\epsilon$ independent of $v^{*}$ and $\phi ^{*}$. The remainder of the proof follows along the same lines as the proof of (a). This completes the proof of the proposition. $\square$

Remark 2.10

We point out that:

• (a) Relying on our assumption (1.5), and combining this with the results of Proposition 2.6 and the Fredholm alternative, then implies the well-posedness of our problem (1.4). In particular, the Fredholm alternative also yields that (1.5) is equivalent to the following condition (for example, see [33, Theorem 2.27] or [21, Chapter 5.3])

  $$\begin{aligned}&w^{*}\in H^s(\mathbb {R}^n) \text { is a solution of }(-\Delta )^s w^{*}-\nabla \cdot (b w^{*} ) +cw^{*}=0\text { in }\Omega \text{ with } w^{*} = 0 \hbox { in } \Omega _e, \\&\text {then we have }w^{*} \equiv 0. \end{aligned}$$

• (b) The proofs of Lemma 2.4 and Proposition 2.6 relied on the fact that the fractional Laplacian was the leading order operator of our equations which was ensured by the condition $\frac{1}{2}<s<1$. In particular, the above arguments and results do not persist in the regime $0<s\le \frac{1}{2}$.

Now, let $\Omega \subset \mathbb {R}^n$ be a bounded, non-empty open set and $\frac{1}{2}<s<1$. We consider the Dirichlet problem (2.2) with a zero source function, i.e.,

$$\begin{aligned} \begin{aligned} ((-\Delta )^s+b\cdot \nabla +c)u=0&\text { in }\Omega ,\\ u-f&\in \widetilde{H}^s(\Omega ), \end{aligned} \end{aligned}$$

(2.15)

for some $f\in H^s(\Omega _e)$. Recall that the eigenvalue condition (1.5) in combination with Proposition 2.6 shows that there exists a unique solution $u\in H^s(\mathbb {R}^n)$ of (2.15). We would like to point out that the solution u here depends only on f modulo $\widetilde{H}^s(\Omega )$. In effect, we emphasize that the solution of the Dirichlet problem (2.15) is only affected by the exterior datum.

Therefore, following [18, 20], we introduce the quotient space

$$\begin{aligned} \mathbb {X} = H^s(\mathbb {R}^n) / \widetilde{H}^s(\Omega ). \end{aligned}$$

We will denote the equivalence class of $f \in H^s(\mathbb {R}^n)$ by [f]. We also recall that for $\Omega$ a bounded, non-empty open Lipschitz set, we have $\mathbb {X}=H^s(\Omega _e)$ by Remark 2.1. Based on the well-posedness of (1.4), one can define the DN map rigorously as follows.

Definition 2.11

(DN map) Let $\Omega \subset \mathbb {R}^n$, $n\ge 1$, be a bounded Lipschitz domain, let $s\in (\frac{1}{2},1)$ and let $b\in W^{1-s,\infty }(\Omega )^n$, $c\in L^{\infty }(\Omega )$. Let $B_{b,c}(\cdot ,\cdot )$ be given as (2.3). Then we define the Dirichlet-to-Neumann operator associated with the Eq. (1.4) as

$$\begin{aligned} \Lambda _{b,c}: \mathbb {X} \rightarrow \mathbb {X}^{*},\ (\Lambda _{b,c}[f], [g]) = B_{b,c}(u_f, g), \end{aligned}$$

where $f, g\in H^{s}(\mathbb {R}^n)$ and $u_f$ is the weak solution to (1.4) with exterior datum f.

We first remark that this is well-defined, since by the definition of a solution to (1.4) we have $B_{b,c}(u_f,\widetilde{g})=B_{b,c}(u_f, \widetilde{g}+ \psi )$ for any $\psi \in \widetilde{H}^s(\Omega )$. Also $B_{b, c}(u_{f + \varphi }, g) = B_{b, c}(u_f, g)$ for $\varphi \in \widetilde{H}^s(\Omega )$, since $u_{f + \varphi } = u_f$ (by the uniqueness of solutions). Finally, the boundedness of the DN map $\Lambda _{b, c}: \mathbb {X} \rightarrow \mathbb {X}^*$ follows easily, once we have the multiplier estimate (2.8) and part (a) of Proposition 2.6.

Remark 2.12

Furthermore, we remark that a formal calculation using Definition 2.11 yields

$$\begin{aligned} (\Lambda _{b,c}[f],[g])&=B_{b,c}(u_{f},g)\nonumber \\&=\int _{\mathbb {R}^{n}}(-\Delta )^{s/2}u_f(-\Delta )^{s/2}gdx+\int _\Omega b\cdot \nabla u_f gdx+\int _\Omega cu_f gdx\nonumber \\&=\int _{\Omega _{e}}g (-\Delta )^s u_f dx, \end{aligned}$$

(2.16)

where $u_f$ is the weak solution to (1.4) with exterior datum f. Then from (2.16), one can formally obtain that

$$\begin{aligned} \Lambda _{b,c}[f]=\left. (-\Delta )^s u_{f}\right| _{\Omega _{e}}, \end{aligned}$$

which gives evidence of (1.6). We remark that in order to make this calculation rigorous, higher regularity assumptions have to be imposed. We refer to Section 3, Lemma 3.1, in [20] for details on this.

We may also consider the adjoint problem with a zero source term:

$$\begin{aligned} \begin{aligned} (-\Delta )^s u^{*} - \nabla \cdot (b u^{*}) + c u^{*} = 0&\text{ in } \Omega ,\\ u^{*} - f&\in \widetilde{H}^s(\Omega ), \end{aligned} \end{aligned}$$

(2.17)

Then by using the eigenvalue condition (1.5) together with Proposition 2.6, this problem has a unique solution $u^*$. With slight abuse of notation, one can analogously define the DN map $\Lambda _{b,c}^{*}$ of the adjoint bilinear form (2.4), i.e.,

$$\begin{aligned} \Lambda _{b,c}^{*}: \mathbb {X} \rightarrow \mathbb {X}^{*},\ (\Lambda _{b,c}^{*}[f], [g]) = B_{b,c}^{*}(u_f^*, g). \end{aligned}$$

Here $u_f^* \in H^s(\mathbb {R})^n$ is a weak solution of (2.17) with exterior data f.

Next, let $\Lambda '_{b,c}: \mathbb {X} \rightarrow \mathbb {X}^*$ be the adjoint operator with respect to the DN map $\Lambda _{b,c}$, i.e., the adjoint DN map $\Lambda '_{b,c}$ is defined via

$$\begin{aligned} (\Lambda _{b,c}f,g) = (f, \Lambda '_{b,c}g), \text { for }f,g\in H^s(\mathbb {R}^n). \end{aligned}$$

Then we further observe the following relation between the DN map associated with (1.4) and the DN map of the adjoint bilinear form, which states that the adjoint DN map $\Lambda '_{b,c}$ and the DN map $\Lambda _{b,c}^*$ with respect to the adjoint equation coincide. In order to simplify notation, here and in the sequel, we drop the brackets $[\cdot ]$ for the elements of $\mathbb {X}$.

Lemma 2.13

Let $\Omega \subset \mathbb {R}^n$, $n\ge 1$, be a bounded, non-empty open Lipschitz set, $s\in (\frac{1}{2},1)$, $b\in W^{1-s,\infty }(\Omega )^n$, $c\in L^{\infty }(\Omega )$ and assume that (1.5) holds. Let $B_{b,c}(\cdot ,\cdot )$, $B_{b,c}^{*}(\cdot ,\cdot )$ and $\Lambda _{b,c}, \Lambda _{b,c}^{*}$ be given as above. Then, for all $f,g \in \mathbb {X}$

$$\begin{aligned} (\Lambda _{b,c}f,g) = (f, \Lambda _{b,c}^{*}g). \end{aligned}$$

(2.18)

Proof

The claim follows from the independence of the DN map and the adjoint DN map of the extension of g into $\Omega$ in the duality pairing from Definition 2.11. Let $u_f$ be a solution to (2.2) with exterior data $f\in \mathbb {X}$ and let $u^{*}_g$ be a solution of the adjoint equation (2.17) with exterior data $g \in \mathbb {X}$. Then, by (2.11) and the definition of $\Lambda _{b,c}$, $\Lambda _{b,c}^{*}$

$$\begin{aligned} (\Lambda _{b,c}f,g)&= B_{b,c}(u_f, u^{*}_g) = B_{b,c}^{*}(u^{*}_g, u_f) = (\Lambda _{b,c}^{*}g,f) = (f,\Lambda _{b,c}^{*}g). \end{aligned}$$

$\square$

With this at hand, we seek to derive a corresponding Alessandrini type identity. It will play a key role in our uniqueness and stability arguments.

Lemma 2.14

(Alessandrini identity) Let $\Omega \subset \mathbb {R}^n$, $n\ge 1$, be a bounded Lipschitz domain and $\frac{1}{2}<s<1$. Let $b_j\in W^{1-s,\infty }(\Omega )^n$ and $c_j\in L^\infty (\Omega )$ be such that (1.5) holds for $j=1,2$. For any $f_1, f_2\in \mathbb {X}$, we have

$$\begin{aligned} \left( (\Lambda _{b_1,c_1}-\Lambda _{b_2,c_2})f_1,f_2\right) _{\mathbb {X}^{*}\times \mathbb {X}}=((b_1-b_2)\cdot \nabla u_1,u_2^{*})_\Omega +((c_1-c_2)u_1,u_2^{*})_\Omega , \end{aligned}$$

(2.19)

where $u_1 \in H^s(\mathbb {R}^n)$ is the solution to $((-\Delta )^s+b_1\cdot \nabla +c_1)u_1=0$ in $\Omega$ and $u_2^{*}\in H^s(\mathbb {R}^n)$ is the solution to $(-\Delta )^s u_2^{*}-\nabla \cdot (b_2 u_2^{*}) + c_2 u_2^{*}=0$ in $\Omega$ with $u_1 - f_1 \in \widetilde{H}^s(\Omega )$ and $u_2 - f_2 \in \widetilde{H}^s(\Omega )$.

Proof

By (2.18), one has

$$\begin{aligned} ((\Lambda _{b_1,c_1}-\Lambda _{b_2,c_2})f_1,f_2)&=(\Lambda _{b_1,c_1}f_1,f_2)-(f_1,\Lambda _{b_2,c_2}^{*}f_2) \\&= B_{b_1,c_1}(u_1,u_2^{*})-B_{b_2,c_2}^{*}(u_2^{*},u_1)\\&= ((b_1-b_2)\cdot \nabla u_1,u_2^{*})_\Omega +((c_1-c_2)u_1,u_2^{*})_\Omega . \end{aligned}$$

$\square$

Last but not least, we define the Poisson operator associated with the Eq. (1.4):

Definition 2.15

Let $\Omega \subset \mathbb {R}^n$, $n\ge 1$, be a bounded, non-empty open Lipschitz set. Let $s\in (\frac{1}{2},1)$ and assume that (1.5) holds. Then, the Poisson operator $P_{b,c}$ associated with (1.4) is defined as

$$\begin{aligned} P_{b,c}: \mathbb {X} \rightarrow H^{s}(\mathbb {R}^n), \ f \mapsto u_f, \end{aligned}$$

(2.20)

where $u_f \in H^s(\mathbb {R}^n)$ with $u_f-f\in \widetilde{H}^s(\Omega )$ denotes the unique solution to (1.4).

3 Approximation property

In this section we discuss the Runge approximation property for solutions to (1.4). To this end, we first recall the strong uniqueness property for the fractional Laplacian, which was proved in [20, Theorem 1.2].

Proposition 3.1

(Global weak unique continuation) Let $n\ge 1$, $s\in (0,1)$ and $u\in H^{-r}(\mathbb {R}^n)$ for some $r\in \mathbb {R}$. Assume that for some non-empty open set $W\subset \mathbb {R}^n$ we have

$$\begin{aligned} u=(-\Delta )^su =0 \text{ in } W. \end{aligned}$$

Then $u\equiv 0$ in $\mathbb {R}^n$.

By using the above strong uniqueness property and a duality argument, one can derive the following Runge approximation [(which then immediately entails Theorem 1.2 (a)].

Lemma 3.2

For $n\ge 1$ and $\frac{1}{2}<s<1$, let $\Omega \subset \mathbb {R}^n$ be a bounded, non-empty open Lipschitz set, $b\in W^{1-s,\infty }(\Omega )^n$ and $c\in L^\infty (\Omega )$ and let $P_{b,c}$ denote the Poisson operator from Definition 2.15. Further let $W\subset \Omega _e$ be an arbitrary non-empty open set. Let

$$\begin{aligned} \mathcal {D}:=\left\{ P_{b,c}f-f: \ f\in C^\infty _c (W)\right\} . \end{aligned}$$

Then we have the following results:

• (a) The set $\mathcal {D}$ is dense in $L^2(\Omega )$.

• (b) The set $\mathcal {D}$ is dense in $\widetilde{H}^s(\Omega )$.

Proof

The proofs of (a) and (b) are similar and follow from a duality argument. Hence, we only need to prove the case (b). First, it is easy to see that $\mathcal D \subset \widetilde{H}^s(\Omega )$. By the Hahn-Banach theorem, it suffices to show that for any $F\in (\widetilde{H}^s(\Omega ))^*$ such that $F(v)=0$ for any $v\in \mathcal D$, we must have $F\equiv 0$. Since $F(v)=0$ for any $v\in \mathcal D$, we have

$$\begin{aligned} F(P_{b,c}f-f)=0 \text { for }f\in C^\infty _c (W). \end{aligned}$$

(3.1)

Next, we claim that

$$\begin{aligned} F(P_{b,c}f-f)=-B_{b,c}(f,\varphi ), \text { for }f\in C^\infty _c(W), \end{aligned}$$

(3.2)

where $\varphi \in \widetilde{H}^s(\Omega )$ is the unique solution of

$$\begin{aligned} (-\Delta )^s \varphi -\nabla \cdot (b \varphi ) +c \varphi =F \text { in }\Omega \text { with }\varphi =0 \text { in }\Omega _e. \end{aligned}$$

We remark that by Proposition 2.6 and the assumption (1.5) this problem is well-posed. In its weak form, it becomes

$$\begin{aligned} B_{b,c}^{*}(\varphi , w)=F(w) \text{ for } w\in \widetilde{H}^s(\Omega ). \end{aligned}$$

We next address the proof of (3.2). Let $f\in C^\infty _c(W)$ and $u_f:=P_{b,c}f\in H^s(\mathbb {R}^n)$, then $u_f-f\in \widetilde{H}^s(\Omega )$ and

$$\begin{aligned} F(P_{b,c}f-f)=B_{b,c}^{*}(\varphi ,u_f-f)= B_{b,c}(u_f-f,\varphi )=-B_{b,c}(f,\varphi ), \end{aligned}$$

where we have utilized (2.11) and the fact that $u_f$ is a solution to (1.4) and $\varphi \in \widetilde{H}^s(\Omega )$. By means of the relations (3.1) and (3.2), we hence conclude that

$$\begin{aligned} B_{b,c}(f,\varphi )=0 \text { for }f\in C^\infty _c (W). \end{aligned}$$

Since $f =0$ in $\mathbb {R}^n {\setminus } \overline{W}$ and $\varphi \in \widetilde{H}^s(\Omega )$, $B_{b,c}(f,\varphi )=0$ implies that

$$\begin{aligned} 0=((-\Delta )^{s/2}\varphi ,(-\Delta )^{s/2}f)_{\mathbb {R}^n}=((-\Delta )^{s}\varphi ,f)_{\mathbb {R}^n},\text { for }f\in C^\infty _c(W). \end{aligned}$$

In particular, $\varphi \in H^s(\mathbb {R}^n)$ satisfies

$$\begin{aligned} \varphi = (-\Delta )^s \varphi =0 \text { in }W. \end{aligned}$$

By invoking Proposition 3.1, this implies $\varphi \equiv 0$ in $\mathbb {R}^n$ and therefore $F\equiv 0$ as well. $\square$

The above lemma proves Theorem 1.2 (a). The proof of Theorem 1.2 (b) is postponed to the last section of this paper. Without major modifications of the above argument, we also obtain the Runge approximation for the adjoint equation.

Corollary 3.3

For $n\ge 1$ and $\frac{1}{2}<s<1$, let $\Omega \subset \mathbb {R}^n$ be a bounded, non-empty open Lipschitz set, $b\in W^{1-s,\infty }(\Omega )^n$ and $c\in L^\infty (\Omega )$ and let $P_{b,c}^{*}$ denote the Poisson operator for the Eq. (2.5), i.e. let

$$\begin{aligned} P_{b,c}^{*}:\mathbb {X} \rightarrow H^s(\mathbb {R}^n), \ f^*\mapsto u_{f^*}^*, \end{aligned}$$

where $u_{f^*}^*\in H^s(\mathbb {R}^n)$ is the solution of the adjoint equation

$$\begin{aligned} (-\Delta )^su_{f^*}^*-\nabla \cdot (bu_{f^*}^*)+cu_{f^*}^*=0 \text { in }\Omega \text { with }u_{f^*}^*=f^*\text { in }\Omega _e. \end{aligned}$$

Let $W\subset \Omega _e$ be an arbitrary non-empty open set. Then the sets $\mathcal {D}^*:=\left\{ P_{b,c}^{*}f^*-f^*: f^*\in C^\infty _c (W)\right\}$ is dense in $L^2(\Omega )$ and $\widetilde{H}^s(\Omega )$.

4 Global uniqueness

In this section, we prove the global uniqueness result from Theorem 1.1 which follows from the knowledge of the DN map $\Lambda _{b,c}$ and its adjoint.

Proof of Theorem 1.1

Since $\Lambda _{b_1,c_1}f|_{W_2}=\Lambda _{b_2,c_2}f|_{W_2}$ for any $f\in C^\infty _c(W_1)$, where $W_1,W_2$ are arbitrary, but fixed non-empty open sets in $\Omega _e$, by using the Alessandrini identity (2.19), we have

$$\begin{aligned} \int _\Omega \left( (b_1-b_2)\cdot \nabla u_1 u_2^{*}+(c_1-c_2)u_1u_2^{*}\right) dx=0, \end{aligned}$$

(4.1)

where $u_1 \in H^s(\mathbb {R}^n)$ is the solution to $((-\Delta )^s+b_1\cdot \nabla +c_1)u_1=0$ in $\Omega$ and $u_2^{*}\in H^s(\mathbb {R}^n)$ is the solution to $(-\Delta )^s u_2^{*}- \nabla \cdot (b_2 u_2^{*}) + c_2 u_2^{*}=0$ in $\Omega$ with exterior data $u_1|_{\Omega _e}=f_1\in C^\infty _c(W_1)$ and $u_2^{*}|_{\Omega _e}=f_2^{*}\in C^\infty _c(W_2)$. In the sequel, we seek to recover the potential and the drift coefficients separately.

Step 1: Recovery of the potentialc. Let $\psi _2 \in C_c^\infty (\Omega )$ be arbitrary. Then choose $\psi _1 \in C_c^\infty (\Omega )$ such that $\psi _1 = 1$ on the set ${{\,\mathrm{supp}\,}}(\psi _2) \Subset \Omega$. By the Runge approximation of $(-\Delta )^s+b\cdot \nabla +c$ and its adjoint (see Lemma 3.2 and Corollary 3.3), there exist sequences of solutions $\{u_j^1\}_{j=1}^\infty$ and $\{u_j^{2,*}\}_{j=1}^\infty$ in $H^s(\mathbb {R}^n)$ such that

$$\begin{aligned} \begin{aligned}&((-\Delta )^s+b_1\cdot \nabla +c_1)u_j^1= (-\Delta )^s u_j^{2,*} - \nabla \cdot (b_2 u_j^{2,*}) + c_2 u_j^{2,*}=0\text { in }\Omega ,\\&\mathrm {supp}(u_j^1)\subset \overline{\Omega _1}\text { and }\mathrm {supp}(u_j^{2,*})\subset \overline{\Omega _2},\\&u_j^1|_\Omega = \psi _1 +r_j^1 \text { and }u_j^{2,*}|_\Omega = \psi _2+r_j^{2,*}, \end{aligned} \end{aligned}$$

(4.2)

where $\Omega _1, \Omega _2$ are non-empty open sets in $\mathbb {R}^n$ containing $\overline{\Omega }$, and $r_j^1, r_j^{2,*}\rightarrow 0$ strongly in $\widetilde{H}^s(\Omega )$ as $j\rightarrow \infty$.

With these solutions at hand, we observe that as $r_j^1, r_j^{2,*} \in \widetilde{H}^s(\Omega )$

$$\begin{aligned} \begin{aligned} \left| \int \limits _{\Omega } (b_1-b_2) \cdot \nabla u_j^1 u_j^{2,*} dx \right|&\le \left| \int \limits _{\Omega } (b_1 -b_2) \cdot \nabla r_j^{1} r_j^{2,*} dx \right| + \left| \int \limits _{\Omega } (b_1 -b_2) \cdot \nabla r_j^{1} \psi _2 dx \right| \\&\quad + \left| \int \limits _{\Omega } (b_1 -b_2) \cdot \nabla \psi _1 \psi _2 dx \right| + \left| \int \limits _{\Omega } (b_1 -b_2) \cdot \nabla \psi _1 r_j^{2, *} dx \right| \\&\le \Vert \nabla r_j^1\Vert _{H^{s-1}(\Omega )} \Vert (b_1-b_2) r_j^{2,*}\Vert _{\widetilde{H}^{1-s}(\Omega )}\\&\quad + \Vert \nabla r_j^{1}\Vert _{H^{s-1}(\Omega )} \Vert (b_1-b_2) \psi _2\Vert _{\widetilde{H}^{1-s}(\Omega )}\\&\quad + \Vert \nabla \psi _1\Vert _{H^{s-1}(\Omega )} \Vert (b_1-b_2) r_j^{2,*}\Vert _{\widetilde{H}^{1-s}(\Omega )}\\&\le C\Vert b_1-b_2\Vert _{W^{1-s,\infty }(\Omega )} \Vert \nabla r_j^1\Vert _{H^{s-1}(\Omega )} \Vert r_j^{2,*}\Vert _{\widetilde{H}^{1-s}(\Omega )}\\&\quad + C\Vert b_1 - b_2\Vert _{W^{1-s,\infty }(\Omega )}\Vert \nabla r_j^{1}\Vert _{H^{s-1}(\Omega )} \Vert \psi _2\Vert _{\widetilde{H}^{1-s}(\Omega )}\\&\quad + C\Vert b_1-b_2\Vert _{W^{1-s,\infty }(\Omega )} \Vert r_j^{2, *}\Vert _{\widetilde{H}^{1-s}(\Omega )} \Vert \nabla \psi _1\Vert _{H^{s-1}(\Omega )}\\&\le C\Vert b_1-b_2\Vert _{W^{1-s,\infty }(\Omega )} \Big (\Vert r_j^{1}\Vert _{\widetilde{H}^{s}(\Omega )} \Vert r_j^{2,*}\Vert _{\widetilde{H}^{s}(\Omega )}\\&\quad + \Vert r_j^{1}\Vert _{\widetilde{H}^{s}(\Omega )}\Vert \psi _2\Vert _{\widetilde{H}^{1-s}(\Omega )} + \Vert r_j^{2, *}\Vert _{\widetilde{H}^{s}(\Omega )} \Vert \psi _1\Vert _{\widetilde{H}^{s}(\Omega )}\Big )\\&\rightarrow 0, \text { as }j \rightarrow \infty . \end{aligned} \end{aligned}$$

(4.3)

Here we have used that $b_1-b_2\in W^{1-s,\infty }(\Omega )$ is a bounded multiplier from $\widetilde{H}^{1-s}(\Omega )$ into itself (which follows from the same argument as in the proof of Proposition 2.6) and the estimate $\Vert u\Vert _{\widetilde{H}^{1-s}(\Omega )} \le C \Vert u\Vert _{\widetilde{H}^{s}(\Omega )}$ for $s\in (\frac{1}{2},1)$. Further, we made strong use of the support assumptions for $\psi _1$ and $\psi _2$, which in particular allow us to drop the third term in the first estimate in (4.3).

Inserting these solutions $\{u_j^1\}$, $\{u_j^{2,*}\}$ and the estimate (4.3) into (4.1) and taking $j\rightarrow \infty$, together with the assumptions on $\psi _1$ and $\psi _2$, we derive

$$\begin{aligned} \int _\Omega (c_1-c_2) \psi _2 dx =0. \end{aligned}$$

Since $\psi _2 \in C_c^\infty (\Omega )$ was arbitrary, by density of $C_c^\infty (\Omega )$ in $L^2(\Omega )$ and the previous identity, we obtain that $c_1 = c_2$ in $\Omega$.

Step 2: Recovery of the driftb. Since we have $c_1 = c_2$ in $\Omega$, (4.1) becomes

$$\begin{aligned} \int _\Omega (b_1-b_2)\cdot \nabla u_1 u_2^{*} dx =0. \end{aligned}$$

(4.4)

Fix an arbitrary $\psi _2 \in C_c^\infty (\Omega )$. Then choose $\psi _{x_k}\in C_c^\infty (\Omega )$ equal $x_k$ on ${{\,\mathrm{supp}\,}}(\psi _2) \Subset \Omega$, where for $k=1,2,\ldots ,n$ the function $x_k$ denotes the restriction to the k-th component of x. By using the Runge approximation again as in (4.2), we can find sequences of solutions $u_j^1 = \psi _{x_k} + r_j^{1}$ and $u_j^{2, *} = \psi _2 + r_j^{2, *}$, with $r_j^{1}, r_j^{2, *} \rightarrow 0$ strongly in $\widetilde{H}^s(\Omega )$ as $j\rightarrow \infty$. Plugging the Runge approximations of the functions $\psi _{x_k}$ and $\psi _2$ into (4.4), we obtain

$$\begin{aligned}&\int _\Omega (b_1-b_2)_k \psi _2 dx + \int \limits _{\Omega } (b_1-b_2) \cdot \nabla r_j^{1} \psi _2 dx + \int \limits _{\Omega } (b_1-b_2) \cdot \nabla r_j^{1} r_j^{2,*} dx\\&\quad + \int \limits _{\Omega } (b_1-b_2) \cdot \nabla {\psi _{x_k}} r_j^{2,*} dx=0, \end{aligned}$$

for $k=1,2,\ldots ,n$, where $(b_1-b_2)_k$ denotes the k-th component of the vector valued function $b_1-b_2$. By arguing similarly as in (4.3), in the limit $j\rightarrow \infty$ we arrive at

$$\begin{aligned} \int _\Omega (b_1-b_2)_k \psi _2 dx = 0 \text { for }k=1,2,\ldots ,n. \end{aligned}$$

Since $\psi _2 \in C_c^\infty (\Omega )$ is arbitrary and as $C_c^\infty (\Omega )$ is dense in $L^2(\Omega )$, we also conclude that $(b_1-b_2)_k=0$ for all $k=1,2,\ldots ,n$. Therefore, $b_1=b_2$, which completes the proof. $\square$

5 Stability

In this section, we study the stability result for the fractional Schrödinger equation with drift.

5.1 Auxiliary results

We begin by proving several auxiliary results, which will be used in deducing a quantitative Runge approximation result in the next section. To this end, we will mainly be studying the dual equation to (1.4)

$$\begin{aligned} \begin{aligned} (-\Delta )^s w - \nabla \cdot (b w) + c w&= v \text{ in } \Omega ,\\ w&= 0 \text{ in } \Omega _e. \end{aligned} \end{aligned}$$

(5.1)

Throughout this section, we assume that the drift field b and the potential c satisfy (1.5).

Lemma 5.1

Let $\Omega \subset \mathbb {R}^n$ be an open, non-empty bounded Lipschitz domain. Let $s\in (\frac{1}{2},1)$, $b\in W^{1-s,\infty }(\Omega )^n$, $c\in L^{\infty }(\Omega )$, $v\in H^{-s}(\Omega )$ and assume that $w\in H^s_{\overline{\Omega }}$ is the solution to (5.1). Then there exists a constant $C>1$ independent of v, w such that

$$\begin{aligned} C^{-1}\Vert v\Vert _{H^{-s}(\Omega )} \le \Vert w\Vert _{H^{s}_{\overline{\Omega }}} \le C \Vert v\Vert _{H^{-s}(\Omega )}. \end{aligned}$$

Proof

The upper bound follows from the well-posedness result of Proposition 2.6 and the assumption (1.5). It hence remains to discuss the lower bound. To this end, we use the triangle inequality and the equation (5.1), which lead to

$$\begin{aligned} \begin{aligned} \Vert v\Vert _{H^{-s}(\Omega )}&\le \Vert (-\Delta )^s w\Vert _{H^{-s}(\Omega )} + \Vert \nabla \cdot (b w)\Vert _{H^{-s}(\Omega )} + \Vert c w\Vert _{H^{-s}(\Omega )}\\&\le \Vert (-\Delta )^s w\Vert _{H^{-s}(\mathbb {R}^n)} + \Vert \nabla \cdot (b w)\Vert _{H^{-s}(\Omega )} + \Vert c w\Vert _{H^{-s}(\Omega )}\\&\le \Vert w\Vert _{H^{s}_{\overline{\Omega }}} + \Vert \nabla \cdot (b w)\Vert _{H^{-s}(\Omega )} + \Vert c w\Vert _{H^{-s}(\Omega )}. \end{aligned} \end{aligned}$$

(5.2)

We estimate the terms with the drift and the potential separately. On the one hand, by integration by parts, we observe that for the drift, we can find a constant $C>0$ independent of b and w such that

$$\begin{aligned} \Vert \nabla \cdot (b w)\Vert _{H^{-s}(\Omega )}&= \sup \limits _{\Vert \phi \Vert _{\widetilde{H}^{s}(\Omega )}=1}\left| (\nabla \cdot (b w), \phi )_{\Omega }\right| \\&\le \sup \limits _{\Vert \phi \Vert _{\widetilde{H}^{s}(\Omega )}=1} \Vert b w\Vert _{\widetilde{H}^{1-s}(\Omega )} \Vert \nabla \phi \Vert _{H^{s-1}(\Omega )}\\&\le C \Vert w\Vert _{\widetilde{H}^{1-s}(\Omega )} \Vert b\Vert _{W^{1-s,\infty }(\Omega )} \sup \limits _{\Vert \phi \Vert _{\widetilde{H}^{s}(\Omega )}=1} \Vert \phi \Vert _{\widetilde{H}^{s}(\Omega )}\\&\le C \Vert w\Vert _{H^s_{\overline{\Omega }}} \Vert b\Vert _{W^{1-s,\infty }(\Omega )} , \end{aligned}$$

where we used that for $s\in (1/2,1)$ it holds that $\Vert w \Vert _{\widetilde{H}^{1-s}(\Omega )} \le C\Vert w \Vert _{\widetilde{H}^{s}_{\overline{\Omega }}}$.

On the other hand, in order to estimate the potential c, we use Hölder’s and Poincaré’s inequalities to observe

$$\begin{aligned} \Vert c w\Vert _{H^{-s}(\Omega )}&= \sup \limits _{\Vert \phi \Vert _{H^{s}_{\overline{\Omega }}}=1}\left| (c w, \phi )_{\Omega }\right| \\&= \sup \limits _{\Vert \phi \Vert _{H^{s}_{\overline{\Omega }}}=1}\Vert c\Vert _{L^{\infty }(\Omega )} \Vert w\Vert _{L^2(\Omega )} \Vert \phi \Vert _{L^2(\Omega )} \le C\Vert c\Vert _{L^{\infty }(\Omega )} \Vert w\Vert _{H^s_{\overline{\Omega }}}. \end{aligned}$$

Inserting the estimates for the contributions involving b, c into (5.2) then concludes the proof of Lemma 5.1. $\square$

Next, we prove Vishik–Eskin type estimates for the fractional Schrödinger equation with drift, c.f. [52] and also [23, Section 3].

Lemma 5.2

(Vishik–Eskin) Let $\delta \in (-\frac{1}{2},\frac{1}{2})$. Let $\Omega \subset \mathbb {R}^n$ be an open, non-empty bounded Lipschitz domain. Let $s\in (\frac{1}{2},1)$, $b\in W^{1-s+\delta ,\infty }(\Omega )^n$, $c\in L^{\infty }(\Omega )$ satisfy (1.5), $v\in H^{-s}(\Omega )$ and assume that $w\in H^s_{\overline{\Omega }}$ is the solution to (5.1). Then there exists a constant $C>1$ such that

$$\begin{aligned} \Vert w\Vert _{H^{s+\delta }_{\overline{\Omega }}} \le C \Vert v\Vert _{H^{-s+\delta }(\Omega )}. \end{aligned}$$

The argument for this follows from a perturbation of the original estimates due to Vishik and Eskin [52].

Proof

We rewrite the Eq. (5.1) as

$$\begin{aligned} (-\Delta )^s w&= G \text{ in } \Omega ,\\ w&= 0 \text{ in } \Omega _e, \end{aligned}$$

where $G= \nabla \cdot (b w) - cw+v$. Then the estimates of Vishik and Eskin (see [24, Theorem 3.1]) yield that

$$\begin{aligned} \Vert w\Vert _{H^{s+\delta }_{\overline{\Omega }}} \le C \Vert G\Vert _{H^{-s+\delta }(\Omega )} \le C (\Vert cw\Vert _{H^{-s+\delta }(\Omega )} + \Vert \nabla \cdot (bw)\Vert _{H^{-s+\delta }(\Omega )} + \Vert v\Vert _{H^{-s+\delta }(\Omega )}). \end{aligned}$$

(5.3)

We estimate the drift and the potential contributions separately: for $\delta \in (0,\min \{2s-1,\frac{1}{2}\})$ by integration by parts and duality

$$\begin{aligned} \Vert \nabla \cdot (b w)\Vert _{H^{-s+\delta }(\Omega )}&= \sup \limits _{\Vert \phi \Vert _{\widetilde{H}^{s-\delta }(\Omega )}=1}\left| (\nabla \cdot (bw), \phi )_{L^2(\Omega )}\right| \\&= \sup \limits _{\Vert \phi \Vert _{\widetilde{H}^{s-\delta }(\Omega )}=1}\left| (bw , \nabla \phi )_{L^2(\Omega )}\right| \\&\le \sup \limits _{\Vert \phi \Vert _{\widetilde{H}^{s-\delta }(\Omega )}=1} \Vert b w\Vert _{\widetilde{H}^{1-s+\delta }(\Omega )} \Vert \nabla \phi \Vert _{H^{s-1-\delta }(\Omega )}\\&\le \sup \limits _{\Vert \phi \Vert _{\widetilde{H}^{s-\delta }(\Omega )}=1}\Vert b\Vert _{W^{1-s+\delta ,\infty }(\Omega )} \Vert w\Vert _{\widetilde{H}^{1-s+\delta }(\Omega )} \Vert \phi \Vert _{\widetilde{H}^{s-\delta }(\Omega )}\\&\le C \Vert b\Vert _{W^{1-s+\delta ,\infty }(\Omega )} \Vert w\Vert _{H^{s}_{\overline{\Omega }}}. \end{aligned}$$

The potential term is estimated by Hölder’s and Poincaré’s inequalities

$$\begin{aligned} \Vert c w\Vert _{H^{-s+\delta }(\Omega )}&= \sup \limits _{\Vert \phi \Vert _{H^{s}_{\overline{\Omega }}}=1}\left| (c w, \phi )_{\Omega }\right| \\&= \sup \limits _{\Vert \phi \Vert _{H^{s-\delta }_{\overline{\Omega }}}=1}\Vert c\Vert _{L^{\infty }(\Omega )} \Vert w\Vert _{L^2(\Omega )} \Vert \phi \Vert _{L^2(\Omega )}\\&\le C\Vert c\Vert _{L^{\infty }(\Omega )} \Vert w\Vert _{H^s_{\overline{\Omega }}}. \end{aligned}$$

Combining these bounds with the estimate $\Vert w\Vert _{H^{s}_{\overline{\Omega }}} \le C \Vert v\Vert _{H^{-s}(\Omega )}$, which follows from the well-posedness result of Proposition 2.6 and returning to (5.3) concludes the argument. $\square$

5.2 A quantitative approximation result

As a final preparation for the stability proof, in this section we deduce a quantitative Runge approximation result for fractional Schrödinger equations with drift terms:

Proposition 5.3

For $n\ge 1$, let $\Omega \subset \mathbb {R}^n$ be a bounded, non-empty open set with a $C^\infty$-smooth boundary and $s\in (\frac{1}{2},1)$. Let $W \Subset \Omega _e$ be non-empty and open such that $\overline{\Omega } \cap \overline{W} = \emptyset$. Further suppose that $\delta \in (0,\frac{2s-1}{2})$ and that $b\in W^{1-s+\delta ,\infty }(\Omega )^n$, $c \in L^{\infty }(\Omega )$. Then, for each $\epsilon >0$ and for each $\overline{v}\in H^{s}_{\overline{\Omega }}$ there exists $f_{\epsilon } \in H^{s}_{\overline{W}}$ such that the following approximation estimates hold true

$$\begin{aligned} \Vert P_{b,c} f_{\epsilon } - f_{\epsilon } - \overline{v}\Vert _{H^{s-\delta }_{\overline{\Omega }}} \le \epsilon \Vert \overline{v}\Vert _{H^{s}_{\overline{\Omega }}}, \ \Vert f_{\epsilon }\Vert _{H^{s}_{\overline{W}}} \le C e^{C \epsilon ^{-\mu (\delta )}} \Vert \overline{v}\Vert _{H^{s-\delta }_{\overline{\Omega }}}. \end{aligned}$$

As in [38] the approximation property follows from a quantitative unique continuation result:

Proposition 5.4

For $n\ge 1$, let $\Omega \subset \mathbb {R}^n$ be a bounded, non-empty open set with a $C^\infty$-smooth boundary and $s\in (\frac{1}{2},1)$. Let $W \Subset \Omega _e$ be non-empty and open such that $\overline{\Omega } \cap \overline{W} = \emptyset$. Further suppose that $\delta \in (0,\frac{2s-1}{2})$ and that $b\in W^{1-s+\delta ,\infty }(\Omega )^n$, $c \in L^{\infty }(\Omega )$.

Assume that for each $v\in H^{s-\delta }_{\overline{\Omega }}$ it holds that

$$\begin{aligned} \Vert v\Vert _{H^{s-2\delta }_{\overline{\Omega }}} \le \frac{C}{\left| \log \left( C \frac{\Vert v\Vert _{H^{s-\delta }_{\overline{\Omega }}}}{\Vert (-\Delta )^s w\Vert _{H^{-s}(\Omega )}} \right) \right| ^{\sigma (\delta )} } \Vert v\Vert _{H^{s-\delta }_{\overline{\Omega }}}, \end{aligned}$$

(5.4)

where $w\in H^{s}_{\overline{\Omega }}$ is the solution of (5.1). Then, for each $\epsilon >0$ and for each $\overline{v}\in H^{s}_{\overline{\Omega }}$ there exists $f_{\epsilon } \in H^{s}_{\overline{W}}$ such that the following approximation estimate holds true

$$\begin{aligned} \Vert P_{b,c} f_{\epsilon } - f_{\epsilon } - \overline{v}\Vert _{H^{s-\delta }_{\overline{\Omega }}} \le \epsilon \Vert \overline{v}\Vert _{H^{s}_{\overline{\Omega }}}, \quad \quad \Vert f_{\epsilon }\Vert _{H^{s}_{\overline{W}}} \le C e^{C \epsilon ^{-\mu (\delta )}} \Vert \overline{v}\Vert _{H^{s-\delta }_{\overline{\Omega }}}. \end{aligned}$$

This statement is the exact analogue of Lemma 8.2 in [38] with the argument for Proposition 5.4 following verbatim as in the proof of Lemma 8.2 in [38]: Indeed, the only property of the Eq. (5.1) which is used, is its mapping property. By virtue of the regularity assumptions on b, c and the well-posedness results of Proposition 2.6, solutions to (5.4) also enjoy exactly the same regularity and compactness estimates as the ones from [38].

Proof of Proposition 5.4

We consider the operator

$$\begin{aligned} A:H^{s}_{\overline{W}} \rightarrow H^{s}_{\overline{\Omega }} \hookrightarrow H^{s-\delta }_{\overline{\Omega }}, \ f \mapsto j(P_{b,c}(f)-f), \end{aligned}$$

where $j:H^{s}_{\overline{\Omega }} \hookrightarrow H^{s-\delta }_{\overline{\Omega }}$ is a compact embedding. Thus, A is a compact, injective operator. Here injectivity follows from the strong uniqueness result of Proposition 3.1. In addition, by (a slight adaptation of) Lemma 3.2, it has a dense range. Thus, we may apply the spectral theorem for compact operators and obtain sequences $\{\mu _j\}_{j=1}^{\infty }\subset \mathbb {R}_+$ decreasing, and $\{w_j\}_{j=1}^{\infty } \subset H^{s}_{\overline{W}}$ such that $A^{*}A w_j = \mu _j w_j$.

The set $\{w_j\}_{j=1}^{\infty }$ forms an orthonormal basis with respect to the $H^{s}_{\overline{W}}$ scalar product. By the density of the range of A, it also follows that the set $\{\varphi _j\}_{j=1}^{\infty }:=\left\{ \frac{1}{\sigma _j}A w_j\right\} _{j=1}^{\infty }$ with $\sigma _j := \mu _j^{\frac{1}{2}}$ is an orthonormal basis of $H^{s-\delta }_{\overline{\Omega }}$. As a consequence, $\left\{ (\sigma _j, w_j, \varphi _j)\right\} _{j=1}^\infty \subset \mathbb {R}_+\times H^{s}_{\overline{W}}\times H^{s-\delta }_{\overline{\Omega }}$ is the singular value decomposition of A. By the characterization of $A^{*}$, the assumption (5.4) can be rephrased as the estimate

$$\begin{aligned} \Vert v\Vert _{H^{s-2\delta }_{\overline{\Omega }}} \le \frac{C}{\left| \log \left( C \frac{\Vert v\Vert _{H^{s-\delta }_{\overline{\Omega }}}}{\Vert A^{*}v\Vert _{H^{s}_{\overline{W}}}} \right) \right| ^{\sigma (\delta )}} \Vert v\Vert _{H^{s-\delta }_{\overline{\Omega }}}. \end{aligned}$$

(5.5)

Using this and the singular value decomposition from above, we deduce the approximation property along the same lines as in [38, Lemma 8.2]: Let $\bar{v}\in H^{s}_{\overline{\Omega }}$. For $\alpha \in (0,1)$, let $r_{\alpha }:= \sum \nolimits _{\sigma _j\le \alpha }(\bar{v},w_j)w_j \in H^{s-\delta }_{\overline{\Omega }}$ and $R_{\alpha } \bar{v}:=\sum \limits _{\sigma _j> \alpha }\sigma _j^{-1}(\bar{v},w_j)\varphi _j \in H^{s}_{\overline{W}}$.

Therefore, on the one hand, we have

$$\begin{aligned} \Vert R_{\alpha } (\bar{v})\Vert _{H^{s}_{\overline{W}}} \le \frac{C}{\alpha }\Vert \bar{v}\Vert _{H^{s-\delta }_{\overline{\Omega }}}. \end{aligned}$$

(5.6)

On the other hand,

$$\begin{aligned} \begin{aligned}&\Vert \bar{v}-A R_{\alpha }(\bar{v})\Vert _{H^{s-\delta }_{\overline{\Omega }}}^2 \\&\quad = \sum \limits _{\sigma _j \le \alpha }|(\bar{v},w_j)_{H^{s-\delta }_{\overline{\Omega }}}|^2 \le (\bar{v},r_{\alpha })_{H^{s-\delta }_{\overline{\Omega }}}\\&\quad =(\bar{v},r_{\alpha })_{H^{s-\delta }(\mathbb {R}^n)} \le \Vert \bar{v}\Vert _{H^{s}(\mathbb {R}^n)} \Vert r_{\alpha }\Vert _{H^{s-2\delta }(\mathbb {R}^n)} =\Vert \bar{v}\Vert _{H^{s}_{\overline{\Omega }}} \Vert r_{\alpha }\Vert _{H^{s-2\delta }_{\overline{\Omega }}}\\&\quad \le \Vert \bar{v}\Vert _{H^{s}_{\overline{\Omega }}} \frac{C}{\left| \log \left( C \frac{\Vert r_{\alpha }\Vert _{H^{s-\delta }_{\overline{\Omega }}}}{\Vert A^{*}r_{\alpha }\Vert _{H^{s}_{\overline{W}}}} \right) \right| ^{\sigma (\delta )}} \Vert r_{\alpha }\Vert _{H^{s-\delta }_{\overline{\Omega }}} \le \Vert \bar{v}\Vert _{H^{s}_{\overline{\Omega }}} \frac{C}{\left| \log \left( C \frac{1}{\alpha } \right) \right| ^{\sigma (\delta )}} \Vert r_{\alpha }\Vert _{H^{s-\delta }_{\overline{\Omega }}}. \end{aligned} \end{aligned}$$

(5.7)

Optimizing (5.6), (5.7) in $\alpha$ then implies the claim. $\square$

As a second main ingredient in the proof of Proposition 5.3, we rely on Theorem 5.1 from [38]. This is a propagation of smallness estimate from the boundary into the bulk for the Caffarelli–Silvestre extension of a general function. It does not use the specific equation at hand. For completeness, we recall the statement:

Proposition 5.5

([38], Proposition 5.1) For $n\ge 1$, let $\Omega \subset \mathbb {R}^n$ be an open, non-empty bounded and smooth domain. Let $W \Subset \Omega _e$ be open, non-empty bounded and Lipschitz with $\overline{\Omega }\cap \overline{W} = \emptyset$. Suppose that $s\in (0,1)$ and that $\widetilde{w}$ is a solution to

$$\begin{aligned} \begin{aligned} \nabla \cdot x_{n+1}^{1-2s}\nabla \widetilde{w}&= 0 \text{ in } \mathbb {R}^{n+1}_+,\\ \widetilde{w}&= w \text{ on } \mathbb {R}^n \times \{0\}, \end{aligned} \end{aligned}$$

(5.8)

where $w \in H^{s}(\mathbb {R}^n)$ is a function which vanishes in an open neighbourhood of W. Assume further that for some constants $C_1>0$ and $\delta >0$ one has the a priori bounds

$$\begin{aligned}&\left\| x_{n+1}^{1-2s} \partial _{n+1} \widetilde{w}\right\| _{H^{-s}(W)}\le \eta ,\\&\left\| x_{n+1}^{\frac{1-2s}{2}} \widetilde{w}\right\| _{L^2(\mathbb {R}^n \times [0,C_1])} + \left\| x_{n+1}^{\frac{1-2s}{2}} \nabla \widetilde{w}\right\| _{L^2(\mathbb {R}^{n+1}_+)} + \left\| x_{n+1}^{\frac{1-2s}{2}-\delta } \nabla \widetilde{w}\right\| _{L^2(\mathbb {R}^{n+1}_+)} \le E, \end{aligned}$$

for some constants $\eta , E$ with $\frac{E}{\eta }>1$. Then, there exist constants $C>1$, $\mu >0$ which depend on $n,s,C_1,\delta , \Omega , W$ such that

$$\begin{aligned} \left\| x_{n+1}^{\frac{1-2s}{2}} \widetilde{w}\right\| _{L^2(2\Omega \times [0,1])} + \Vert x_{n+1}^{\frac{1-2s}{2}} \nabla \widetilde{w}\Vert _{L^2(2\Omega \times [0,1])} \le C \frac{E}{\left| \log \left( C\frac{E}{\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}} \right) \right| ^{\mu }}. \end{aligned}$$

Here $2\Omega := \{x\in \mathbb {R}^n: \ {{\,\mathrm{dist}\,}}(x, \Omega ) \le \min \{\frac{1}{2} {{\,\mathrm{dist}\,}}(\overline{W},\overline{\Omega }),2 \}\}$.

We remark that although Proposition 5.1 in [38] is formulated for the Caffarelli–Silvestre extension of the solution to the fractional Schrödinger equation which is studied there (in [38] the situation $b=0$ is considered), only the vanishing of w in a neighbourhood of W is used in the argument which leads to Proposition 5.1. We thus do not present the proof of this result, but refer the reader to Section 5 in [38]. We will apply it to solutions to (5.1) in the sequel.

With Propositions 5.4 and 5.5 at hand, we address the proof of the approximation result of Proposition 5.3:

Proof of Proposition 5.3

By Proposition 5.4, it suffices to argue that (5.4) holds true.

Step 1: Reduction. As in [38] the estimate (5.4) is derived by interpolation from the following quantitative unique continuation result

$$\begin{aligned} \Vert v\Vert _{H^{-s}(\Omega )} \le \frac{C}{ \left| \log \left( \frac{CE }{\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}} \right) \right| ^{\sigma }} E, \end{aligned}$$

(5.9)

where $E \ge \Vert v\Vert _{L^2(\Omega )}$. We recall the argument that (5.9) implies (5.4) (c.f. Step 1 in the proof of Theorem 1.2 in [38]): For $\theta \in (0,1)$ such that $0\le s\theta < \frac{1}{2}$ we first interpolate (5.9) with some $E \ge \Vert v\Vert _{L^2(\Omega )}$. This yields

$$\begin{aligned} \Vert v\Vert _{H^{-\theta s}(\overline{\Omega })} \le C^{\theta }\frac{1}{\left| \log \left( \frac{CE }{\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}} \right) \right| ^{\sigma (\delta )}} E. \end{aligned}$$

Using that $H^{-\theta s}(\Omega ) = H^{-\theta s}_{\overline{\Omega }}$ for $0<\theta s < \frac{1}{2}$, we then interpolate once more with the trivial bound $\Vert v\Vert _{H^{s-\delta }_{\overline{\Omega }}} \le \Vert v\Vert _{H^{s-\delta }_{\overline{\Omega }}}$. Choosing $E \ge \Vert v\Vert _{H^{s-\delta }_{\overline{\Omega }}} \ge \Vert v\Vert _{L^2(\Omega )}$ then yields

$$\begin{aligned} \Vert v\Vert _{H^{-s+2\delta }_{\overline{\Omega }}} \le \frac{C}{ \left| \log \left( \frac{CE }{\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}} \right) \right| ^{\widetilde{\sigma }(\delta )} }E. \end{aligned}$$

We again emphasize that this reduction does not rely on any properties of the Eq. (5.1) but only follows from general interpolation arguments.

Step 2: Proof of (5.9). Also in the proof of (5.9) only three ingredients involving the solution of (5.1) are exploited:

• (i) It is used that $\Vert w\Vert _{H^s(\mathbb {R}^n)} \le C \Vert v\Vert _{L^2(\Omega )}$,

• (ii) that $\Vert w\Vert _{H^{s+\widetilde{\delta }}(\mathbb {R}^n)} \le C \Vert v\Vert _{H^{-s+\widetilde{\delta }}(\Omega )} \le C \Vert v\Vert _{L^2(\Omega )}$ for some $\widetilde{\delta }\in (0,1/2)$,

• (iii) and that $\Vert v\Vert _{H^{-s}(\Omega )} \le C \Vert w\Vert _{H^s_{\overline{\Omega }}}$.

Indeed, choosing a constant $E>0$ such that $\Vert v\Vert _{L^2(\Omega )}\le E$ and assuming that (i), (ii) hold, the properties of the Caffarelli–Silvestre extension yield

• (i’) $\left\| x_{n+1}^{\frac{1-2s}{2}} \widetilde{w}\right\| _{L^2(\mathbb {R}^n \times [0,C_1])} + \Vert x_{n+1}^{\frac{1-2s}{2}} \nabla \widetilde{w}\Vert _{L^2(\mathbb {R}^{n+1}_+)} \le C E$,

• (ii’) $\left\| x_{n+1}^{\frac{1-2s}{2}-\delta }\widetilde{w}\right\| _{L^2(\mathbb {R}^{n+1}_+)} \le C E$.

Here $\widetilde{w}$ denotes the Caffarelli–Silvestre extension of w; for the details of the estimates we refer to the proof of Theorem 1.3 in [38]. If (i’), (ii’) are available and setting $\eta = C\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}$ (so that by the properties of the Caffarelli–Silvestre extension [5], also $\Vert \lim \limits _{x_{n+1}\rightarrow 0} x_{n+1}^{1-2s}\partial _{n+1}\widetilde{w}\Vert _{H^{-s}(W)} \le \eta$), Theorem 5.1 in [38] is applicable and results in

$$\begin{aligned} \left\| x_{n+1}^{\frac{1-2s}{2}} \widetilde{w}\right\| _{L^2(2\Omega \times [0,1])} + \Vert x_{n+1}^{\frac{1-2s}{2}} \nabla \widetilde{w}\Vert _{L^2(2\Omega \times [0,1])} \le C \frac{E}{\left| \log \left( C\frac{E}{\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}} \right) \right| ^{\mu }}. \end{aligned}$$

Recalling localized trace estimates (c.f. Lemma 4.4 in [38]) then gives

$$\begin{aligned} \Vert w\Vert _{H^s_{\overline{\Omega }}} \le C \frac{E}{\left| \log \left( C\frac{E}{\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}} \right) \right| ^{\mu }}. \end{aligned}$$

Finally, applying (iii), we can further bound $\Vert w\Vert _{H^{s}_{\overline{\Omega }}}$ from below and obtain

$$\begin{aligned} \Vert v\Vert _{H^{-s}(\Omega )} \le C \frac{E}{\left| \log \left( C\frac{E}{\Vert (-\Delta )^s w\Vert _{H^{-s}(W)}} \right) \right| ^{\mu }}. \end{aligned}$$

Since Lemmas 5.1, 5.2 imply the conditions (i), (ii), (iii), this concludes the proof of Proposition 5.3. $\square$

5.3 Proof of Theorem 1.3

With the approximation result of Proposition 5.3 at hand, we present the proof of the stability estimate from Theorem 1.3.

Proof of Theorem 1.3

As in the uniqueness proof we argue in two steps:

Step 1: Stability estimate forc. First, we seek to obtain an estimate on $\Vert c_1-c_2\Vert _{H^{-s}(\Omega )}$. Let $\psi _1$ be a smooth function which is compactly supported in $\Omega$ and which satisfies $\psi _1 = 1$ in ${{\,\mathrm{supp}\,}}(c_1), {{\,\mathrm{supp}\,}}(c_2), {{\,\mathrm{supp}\,}}(b_1), {{\,\mathrm{supp}\,}}(b_2)$. Let further $g\in \widetilde{H}^{s}(\Omega )$ be arbitrary. Using the quantitative Runge approximation result from Proposition 5.3 and considering approximate solutions

$$\begin{aligned} u^1_j = \psi _1 + r_j^1, \ u^{2,*}_j = g + r_j^2, \end{aligned}$$

with associated exterior data $f_{j}^1, f_{j}^2$, we infer from Alessandrini’s identity

$$\begin{aligned} \begin{aligned} ((c_1-c_2), g)_{\Omega }&= ((\Lambda _{b_1,c_1}-\Lambda _{b_2,c_2})f_{j}^1,f_{j}^2)_{W_2} -((c_1-c_2)r_j^1, g)_{\Omega } \\&\quad - ((c_1-c_2)r_j^1, r_j^2)_{\Omega } - ((b_1-b_2) \cdot \nabla r_{1}^j, u_{j}^{2,*})_{\Omega }. \end{aligned} \end{aligned}$$

(5.10)

We estimate the terms involving the drift and the potential as follows

$$\begin{aligned} \begin{aligned} |((b_1-b_2)\cdot \nabla r_j^1, u_j^{2,*})_{\Omega }|&\le \Vert \nabla r_j^1\Vert _{H^{s-1-\delta }(\Omega )} \Vert (b_1-b_2) u_j^{2,*}\Vert _{\widetilde{H}^{1-s+\delta }(\Omega )}\\&\le C\Vert r_j^1\Vert _{H^{s-\delta }(\Omega )} \Vert b_1-b_2\Vert _{W^{1-s+\delta ,\infty }(\Omega )} \Vert u_j^{2,*}\Vert _{\widetilde{H}^{1-s+\delta }(\Omega )}\\&\le C M \Vert r_j^1\Vert _{H^{s-\delta }(\Omega )} \Vert u_j^{2,*}\Vert _{\widetilde{H}^{s-\delta }(\Omega )}, \end{aligned} \end{aligned}$$

(5.11)

where we used that $s\in (\frac{1}{2},1)$ and chose $\delta \in (0,\frac{2s-1}{2})$ sufficiently small, and

$$\begin{aligned} \begin{aligned} |((c_1-c_2) r_j^1, \psi _1)_{\Omega }|&\le \Vert c_1-c_2\Vert _{L^{\infty }(\Omega )} \Vert r_j^1\Vert _{L^2(\Omega )} \Vert \psi _1\Vert _{L^2(\Omega )}\\&\le CM \Vert r_j^1\Vert _{\widetilde{H}^{s-\delta }(\Omega )}\Vert \psi _1\Vert _{L^2(\Omega )}. \end{aligned} \end{aligned}$$

(5.12)

Applying the bounds from Proposition 5.3 and combining the bounds from (5.11), (5.12) with (5.10), we therefore infer

$$\begin{aligned} |((c_1-c_2), g)_{\Omega }| \le C(\Vert \Lambda _{b_1,c_1}-\Lambda _{b_2,c_2}\Vert _{*} e^{C \epsilon ^{-\mu }} + \epsilon M ) \Vert \psi _1\Vert _{H^{s}_{\overline{\Omega }}} \Vert g\Vert _{H^{s}_{\overline{\Omega }}}. \end{aligned}$$

Optimizing in $\epsilon >0$ by choosing $\epsilon = \left| \log (\Vert \Lambda _{b_1,c_2}-\Lambda _{b_2,c_2}\Vert _{*})\right| ^{-\nu }$ for some $\nu >0$ then implies the estimate

$$\begin{aligned} \Vert c_1-c_2\Vert _{H^{-s}(\Omega )}&\le \sup \limits _{\Vert g\Vert _{H^s_{\overline{\Omega }}=1}} \left| ((c_1-c_2),g)_{\Omega }\right| \\&\le \sup \limits _{\Vert g\Vert _{H^s_{\overline{\Omega }}=1}} C\left| \log (\Vert \Lambda _{b_1,c_2}-\Lambda _{b_2,c_2}\Vert _{*})\right| ^{-\nu } \Vert \psi _1\Vert _{H^s_{\overline{\Omega }}} \Vert g\Vert _{H^s_{\overline{\Omega }}}\\&\le C\left| \log (\Vert \Lambda _{b_1,c_2}-\Lambda _{b_2,c_2}\Vert _{*})\right| ^{-\nu } \Vert \psi _1\Vert _{H^s_{\overline{\Omega }}} . \end{aligned}$$

Since $\psi _1$ is fixed, we arrive at the estimate for $\Vert c_1-c_2\Vert _{H^{-s}(\Omega )}$.

Step 2: Stability estimate forb. As a preparation for the stability estimate for b, we next note that by interpolation, for $\widetilde{\delta }>0$ and some $\theta \in (0,1)$, which only depends on $n,s,\widetilde{\delta }$, we have

$$\begin{aligned} \begin{aligned}&\Vert c_1-c_2\Vert _{L^{\frac{n}{2s}+\widetilde{\delta }}(\Omega )} \\&\quad \le C \Vert c_1-c_2\Vert _{L^2(\Omega )}^{\theta } \Vert c_1-c_2\Vert _{L^{\infty }(\Omega )}^{1-\theta } \\&\quad \le C \Vert c_1-c_2\Vert _{L^2(\Omega )}^{\theta } \Vert c_1-c_2\Vert _{W^{1,n+\delta }(\Omega )}^{1-\theta }\\&\quad \le C \Vert c_1-c_2\Vert _{H^{-s}(\Omega )}^{\frac{\theta }{2}} \Vert c_1-c_2\Vert _{\widetilde{H}^{s}(\Omega )}^{\frac{\theta }{2}} \Vert c_1-c_2\Vert _{W^{1,n+\delta }(\Omega )}^{1-\theta }\\&\quad \le C \left( \left| \log (\Vert \Lambda _{b_1,c_2}-\Lambda _{b_2,c_2}\Vert _{*})\right| ^{-\nu } \Vert \psi _1\Vert _{H^s_{\overline{\Omega }}} \right) ^{\frac{\theta }{2}} \Vert c_1-c_2\Vert _{\widetilde{H}^{s}(\Omega )}^{\frac{\theta }{2}} \Vert c_1-c_2\Vert _{W^{1,n+\delta }(\Omega )}^{1-\theta }. \end{aligned} \end{aligned}$$

(5.13)

With this estimate at hand, for $k\in \{1,\dots ,n\}$ we again consider approximate solutions

$$\begin{aligned} u_j^1 = \psi _1 x_k + r_j^1, \ u_j^{2,*} = g + r_j^2, \end{aligned}$$

with corresponding exterior data $f_{j}^1, f_{j}^{2,*}$. Here $\psi _1$ is chosen as above, i.e. in particular such that $\psi _1 = 1$ on the support of the functions $b_{1},b_2, c_1,c_2$. The solutions $u_j^1, u_{j}^{2,*}$ are chosen such that the estimates of Proposition 5.3 hold.

As before, we exploit Alessandrini’s identity. Now we have to treat the full term involving $((c_1-c_2) u_{j}^1, u_j^{2,*})_{\Omega }$ as an error: For all $k\in \{1,\dots ,n\}$

$$\begin{aligned} ((b_1-b_2)_k, g)_{\Omega }&= ((\Lambda _{b_1,c_1}-\Lambda _{b_2,c_2})f_{j}^1,f_{j}^2)_{W_2} \nonumber \\&\quad -((c_1-c_2)u_j^1, u_j^{2,*})_{\Omega } - ((b_1-b_2) \cdot \nabla r_{1}^j, u_{2}^{j,*})_{\Omega }. \end{aligned}$$

(5.14)

The terms involving the DN map and the drift term on the right hand side are estimated as in (5.11), (5.12). For the term involving the difference of the potentials, we invoke the interpolation estimate from above; we estimate it as follows

$$\begin{aligned}&|((c_1-c_2) u_{j}^1, u_j^{2,*})_{\Omega }| \\&\quad \le \Vert c_1-c_2\Vert _{L^{\frac{n}{2s}+\widetilde{\delta }}(\Omega )} \Vert u_{j}^1\Vert _{L^{2^{\star }-\widetilde{\delta }}(\Omega )} \Vert u_{j}^{2,*}\Vert _{L^{2^{\star }-\widetilde{\delta }}(\Omega )}\\&\quad \le C\Vert c_1-c_2\Vert _{H^{-s}(\Omega )}^{\frac{\theta }{2}} \Vert c_1-c_2\Vert _{\widetilde{H}^{s}(\Omega )}^{\frac{\theta }{2}} \Vert c_1-c_2\Vert _{W^{1,n+\delta }(\Omega )}^{1-\theta } \Vert u_{j}^1\Vert _{\widetilde{H}^{s-\delta }(\Omega )} \Vert u_{j}^{2,*}\Vert _{\widetilde{H}^{s-\delta }(\Omega )}\\&\quad \le C \left| \log (\Vert \Lambda _{b_1,c_1}-\Lambda _{b_2,c_2}\Vert _{*})\right| ^{-\nu \frac{\theta }{2}} \Vert c_1-c_2\Vert _{\widetilde{H}^{s}(\Omega )}^{\frac{\theta }{2}} \Vert c_1-c_2\Vert _{W^{1,n+\delta }(\Omega )}^{1-\theta } \Vert u_{j}^1\Vert _{\widetilde{H}^{s-\delta }(\Omega )} \Vert u_{j}^{2,*}\Vert _{\widetilde{H}^{s-\delta }(\Omega )}. \end{aligned}$$

Here $2^{\star }= \frac{2n}{n-2s}$ denotes the corresponding (fractional) Sobolev embedding exponent and $\widetilde{\delta }>0$ is sufficiently small. Inserting this into (5.14) and using the bounds from Proposition 5.3, we infer

$$\begin{aligned}&|((b_1-b_2)_k, g)_{\Omega }|\\&\quad \le C (1+M) \left( \Vert \Lambda _{b_1,c_1}-\Lambda _{b_2,c_2}\Vert _{*} e^{C \epsilon ^{-\mu }}\right. \\&\qquad \left. + \left| \log (\Vert \Lambda _{b_1,c_1}-\Lambda _{b_2,c_2}\Vert _{*})\right| ^{-\nu \frac{\theta }{2}} + \epsilon \right) \Vert \psi _1 x_k\Vert _{H^s_{\overline{\Omega }}} \Vert g\Vert _{H^s_{\overline{\Omega }}}. \end{aligned}$$

for any $k\in \{1,\dots ,n\}$. Optimizing once more in $\epsilon$ and taking the supremum over $g\in \widetilde{H}^{s}(\Omega )$ with $\Vert g\Vert _{\widetilde{H}^s(\Omega )}=1$ and over $k\in \{1,\dots ,n\}$, hence also leads to the desired logarithmic stability estimate for the difference of the drifts. This concludes the proof of the stability estimate. $\square$

6 Reconstruction and finite measurements uniqueness

Last but not least, we present a few results on reconstruction procedures and finite measurement statements for the fractional Calderón problem with drift. More precisely, we show that uniqueness for the fractional Calderón problem with $C^\infty _c$ drift and potential can be guaranteed from $n+1$ measurements only.

6.1 Higher order approximation

Before addressing the finite measurement results, we recall the higher order Runge approximation property and present the proof of Theorem 1.2 (b). The structure of the proof for the higher order Runge approximation is similar as the arguments presented in Sect. 3. However, since we seek to approximate solutions in high regularity function spaces, by using a duality argument, we need to consider the corresponding Dirichlet problem in Sobolev spaces of negative orders. The argument for this follows along the same lines as the proofs in [20, Section 7], which in the sequel we recall for self-containedness for the fractional Schrödinger equation with drift.

Let us consider the fractional Laplacian $(-\Delta )^s$ with $s \in (\frac{1}{2},1)$. Let $\Omega \subset \mathbb {R}^n$ be a bounded domain with a $C^\infty$-smooth boundary, $b\in C^\infty _c(\Omega )^n$, $c\in C^\infty _c(\Omega )$ with ${{\,\mathrm{supp}\,}}(b),{{\,\mathrm{supp}\,}}(c)\Subset \Omega$ satisfy (1.5). Here we impose the described compact support as well as the high regularity conditions for b and c in order to satisfy the assumptions from the theory which is presented in [24]. In the sequel, we consider the function space

$$\begin{aligned} \mathcal {E}^s(\overline{\Omega }):=e_\Omega d(x)^s C^\infty (\overline{\Omega }), \end{aligned}$$

where $e_\Omega$ denotes extension by zero from $\Omega$ to $\mathbb {R}^n$, and $d=d(x)$ is a $C^\infty$ function in $\overline{\Omega }$ with $d(x)>0$ for $x\in \Omega$ and $d(x)={{\,\mathrm{dist}\,}}(x,\partial \Omega )$ near $\partial \Omega$.

Further for $\mu > s- \frac{1}{2}$ we work in the Banach space $H^{s(\mu )}(\overline{\Omega })$ which is the space which is introduced in [24] as the Banach space tailored for solutions u solving the problem

$$\begin{aligned} r_\Omega ((-\Delta )^s + b\cdot \nabla +c)u \in H^{\mu -2s}(\Omega ) \text { with }u=0 \text { in }\Omega _e, \end{aligned}$$

where $r_\Omega$ is the restriction map from $\mathbb {R}^n$ to $\Omega$ such that $r_\Omega u=u|_\Omega$.

In order to deduce the desired higher order approximation property, we recall the following result from [24], which was also used in [20, Lemma 7.1] for deducing higher order approximation for the fractional Schrödinger equation.

Proposition 6.1

(Lemma 7.1 in [20]) For $\mu >s-\frac{1}{2}$ and a smooth bounded domain $\Omega \subset \mathbb {R}^n$, there exists a Banach space $H^{s(\mu )}(\overline{\Omega })$ with the following properties:

• (a) $H^{s(\mu )}(\overline{\Omega })\subset H^{s-\frac{1}{2}}_{\overline{\Omega }}$ with a continuous inclusion;

• (b) $H^{s(\mu )}(\overline{\Omega })=H^{\mu }_{\overline{\Omega }}$ for $\mu \in (s-\frac{1}{2},s+\frac{1}{2})$;

• (c) $r_\Omega ((-\Delta )^{s}+b\cdot \nabla +c)$ is a homeomorphism from $H^{s(\mu )}(\overline{\Omega })$ onto $H^{\mu -2s}(\Omega )$;

• (d) $H^{\mu }_{\overline{\Omega }}\subset H^{s(\mu )}(\overline{\Omega })\subset H^{\mu }_{loc}(\Omega )$ with continuous inclusions, or the multiplication by any smooth cut-off $\chi \in C^\infty _c(\Omega )$ is bounded from $H^{s(\mu )}(\overline{\Omega })$ to $H^{\mu }(\Omega )$;

• (e) $\mathcal {E}^s (\overline{\Omega })=\cap _{\mu >s-\frac{1}{2}}H^{s(\mu )}(\overline{\Omega })$ and the set $\mathcal {E}^s (\overline{\Omega })$ is dense in $H^{s(\mu )}(\overline{\Omega })$.

For the proof of this result we refer to [20, 23, 24]. We remark that equipped with the topology induced by $\{\Vert \cdot \Vert _{H^{s(k)}}\}_{k=1}^\infty$, the space $\mathcal {E}^s (\overline{\Omega })$ is a Fréchet space.

Building on these properties of the spaces $H^{s(\mu )}(\overline{\Omega })$, following the argument of [20], we prove a higher order approximation property for solutions to (1.4) in $\mathcal {E}^s (\overline{\Omega })$. The following result was proved by [20, Lemma 7.2 ] for the case $b=0$.

Lemma 6.2

Let $\Omega \subset \mathbb {R}^n$ be a bounded domain with a $C^\infty$-smooth boundary, and $\frac{1}{2}<s <1$. Let $W\subset \Omega _e$ be a non-empty, open set, and let $b\in C^\infty _c(\Omega )^n$, $c\in C^\infty _c(\Omega )$ with ${{\,\mathrm{supp}\,}}(b),{{\,\mathrm{supp}\,}}(c)\Subset \Omega$ be such that (1.5) holds. Let $P_{b,c}$ be the Poisson operator given by (2.20) and

$$\begin{aligned} \mathcal {D}:=\{e_\Omega (r_\Omega P_{b,c}f): f\in C^\infty _c(W)\}. \end{aligned}$$

Then the set $\mathcal {D}$ is dense in the Fréchet space $\mathcal {E}^s (\overline{\Omega })$ with the topology induced by $\{\Vert \cdot \Vert _{H^{s(k)}}\}_{k=1}^\infty$.

Proof

We follow the argument of [20, Lemma 7.2]. First, notice that for any $f\in C^\infty _c(W)$, by the definition of the Poisson operator (2.20), one has $P_{b,c}f=f+v$, where $v \in \widetilde{H}^s(\Omega )$ satisfies

$$\begin{aligned} r_\Omega ((-\Delta )^s+b\cdot \nabla +c)v\in C^\infty (\overline{\Omega }). \end{aligned}$$

By Proposition 6.1, we have $v\in \mathcal {E}^s(\overline{\Omega })$, which implies that $\mathcal {D}\subset \mathcal {E}^s (\overline{\Omega })$. Next, let $\mathcal {L}$ be a continuous linear functional defined on $\mathcal {E}^s (\overline{\Omega })$ satisfying

$$\begin{aligned} \mathcal {L}(e_\Omega (r_\Omega P_{b,c}f))=0, \text { for all }f\in C^\infty _c(W). \end{aligned}$$

By the Hahn-Banach theorem, it suffices to show that $\mathcal {L}\equiv 0$. By using the definition of the topology of the Fréchet space $\mathcal {E}^s (\overline{\Omega })$, one can find an integer $\ell$ so that

$$\begin{aligned} |\mathcal {L}(u)|\le C\sum _{m=1}^{\ell } \Vert u\Vert _{H^{s(m)}(\overline{\Omega })}\le C\Vert u\Vert _{H^{s(\ell )}(\overline{\Omega })}, \text { for }u\in \mathcal {E}^{s}(\overline{\Omega }), \end{aligned}$$

for some constant $C>0$ which is independent of u. By virtue of Proposition 6.1 (e), $\mathcal {E}^{s}(\overline{\Omega })$ is dense in $H^{s(\ell )}(\overline{\Omega })$. Thus, $\mathcal {L}$ has a unique bounded extension $\widetilde{\mathcal {L}}\in (H^{s(\ell )}(\overline{\Omega }))^{*}$.

Let us consider the same homeomorphism in Proposition 6.1 (c),

$$\begin{aligned} \mathcal {T}=r_\Omega ((-\Delta )^s+b\cdot \nabla +c):H^{s(\ell )}(\overline{\Omega })\rightarrow H^{\ell -2s}(\Omega ). \end{aligned}$$

The adjoint of $\mathcal T$ is a bounded map between the dual Banach spaces with

$$\begin{aligned} \mathcal {T}^{*}:(H^{\ell -2s}(\Omega ))^{*}\rightarrow (H^{s(\ell )}(\overline{\Omega }))^*. \end{aligned}$$

Note that the adjoint map $\mathcal {T}^*$ is also homeomorphism with the inverse $(\mathcal {T}^{-1})^{*}$. Moreover, by Remark 2.1, we have $(H^{\ell -2s}(\Omega ))^*=H^{-\ell +2s}_{\overline{\Omega }}$ such that

$$\begin{aligned} \mathcal {T}^{*}v(w)=(v,\mathcal {T}w)_{H^{-\ell +2s}_{\overline{\Omega }} \times H^{\ell -2s}(\Omega )}, \text { for }v\in H^{-\ell +2s}_{\overline{\Omega }} \text { and }w\in H^{s(\ell )}(\overline{\Omega }). \end{aligned}$$

Let $v\in H^{-\ell +2s}_{\overline{\Omega }}$ be the unique function satisfying $\mathcal {T}^{*}v=\widetilde{\mathcal {L}}$ and choose a sequence $\{v_k\}_{k\in \mathbb {N}}\subset C^\infty _c(\Omega )$ with $v_k \rightarrow v$ in $H^{-\ell +2s}(\Omega )$ as $k\rightarrow \infty$. Now, let $f\in C^\infty _c(W)$, recalling that $e_\Omega (r_\Omega P_{b,c}f)=P_{b,c}f-f$, then we have

$$\begin{aligned} 0&=\mathcal {L}(e_\Omega (r_\Omega P_{b,c}f))=\widetilde{\mathcal {L}}(P_{b,c}f-f)=\mathcal {T}^*v(P_{b,c}f-f)\nonumber \\&= (v,\mathcal {T}(P_{b,c}f-f)) = -(v,\mathcal {T}f)=\lim _{k\rightarrow \infty }(v_k,((-\Delta )^s +b\cdot \nabla +c)f)\nonumber \\&=-\lim _{k\rightarrow \infty }((-\Delta )^s v_k -\nabla \cdot (bv_k)+cv_k, f), \end{aligned}$$

(6.1)

where we have utilized that $\mathcal {T}P_{b,c}f=0$ and $v_k\in C^\infty _c(\Omega )$. Finally, since $f\in C^\infty _c(W)$ with $\overline{W}\cap \overline{\Omega }=\emptyset$, then the last equation of (6.1) reads

$$\begin{aligned} ((-\Delta )^s v,f)=\lim _{k\rightarrow \infty }((-\Delta )^s v_k,f)=0, \text { for }f\in C^\infty _c(W). \end{aligned}$$

Thus, we obtain that $v\in H^{-\ell +2s}(\mathbb {R}^n)$ satisfies

$$\begin{aligned} v|_W=(-\Delta )^s v|_W=0, \end{aligned}$$

and the strong uniqueness (Proposition 3.1) implies that $v\equiv 0$ in $\mathbb {R}^n$. Therefore, we obtain $\widetilde{\mathcal {L}}=0$ and hence $\mathcal {L}=0$, which completes the proof. $\square$

Now, we are ready to prove the higher regularity Runge approximation property.

Proof of Theorem1.2(b) As $\Omega \Subset \Omega _1$ with $\mathrm {int}(\Omega _1 {\setminus } \overline{\Omega })\ne \emptyset$, it is possible to find a small ball W with $\overline{W}\subset \Omega _1{\setminus } \overline{\Omega }$. Let $g\in C^\infty (\overline{\Omega })$ and $h:=e_{\Omega } d(x)^{s} g \in \mathcal {E}^{s}(\overline{\Omega })$, then Lemma 6.2 shows that one can find a sequence of solutions $\{u_j\} \subset H^s(\mathbb {R}^n)$ satisfying

$$\begin{aligned} ((-\Delta )^s+b\cdot \nabla +c)u_j=0 \text { in }\Omega \text { with }{{\,\mathrm{supp}\,}}(u_j)\subset \overline{\Omega _1}, \end{aligned}$$

so that $e_\Omega r_\Omega u_j \in \mathcal {E}^s(\overline{\Omega })$ and

$$\begin{aligned} e_\Omega r_\Omega u_j \rightarrow h \text { in }\mathcal {E}^{s}(\overline{\Omega }) \text { as }j\rightarrow \infty . \end{aligned}$$

The higher order approximation will hold if we can show that

$$\begin{aligned} \mathcal {M}:C^\infty (\Omega ) \rightarrow \mathcal {E}^{s}(\overline{\Omega }) \text { with }\mathcal {M}g=e_\Omega d(x)^s g \end{aligned}$$

is a homeomorphism, as it is then possible to apply $\mathcal {M}^{-1}=d(x)^{-s}r_\Omega$. This then gives

$$\begin{aligned} d(x)^{-s}r_\Omega u_j \rightarrow g \text { in }C^\infty (\overline{\Omega }). \end{aligned}$$

Note that the map $\mathcal {M}$ is a bijective linear map between Fréchet spaces and has a closed graph, i.e., if $g_j \rightarrow g$ in $C^\infty$ and $\mathcal {M}g_j\rightarrow h$ in $\mathcal {E}^s$, then also $\mathcal {M}g_j \rightarrow \mathcal {M}g$ in $L^\infty$. Then by the uniqueness of the limit, one obtains that $\mathcal {M}g=h$ as distributional limits. Hence, $\mathcal {M}$ is a homeomorphism by the closed graph and the open mapping theorems. This finishes the proof. $\square$

6.2 Finite measurements reconstruction without openness

In this section, we discuss a first result towards the proof of Theorem 1.4 by using higher order Runge approximation [Theorem 1.2 (b)]. However, before proving the full result of Theorem 1.4, we prove a weaker (but technically considerably easier) result, which still proves finite measurement reconstruction but only asserts that the set of measurement data contains a non-empty open set (a priori this argument does not prove the density of the set of good data). The technically more involved statement on the openness and density of the set of good data will be proved in the subsequent sections.

Proof of Theorem1.4without the density result We show that for any drift $b \in C^{\infty }_c(\Omega )^n$ and any potential $c\in C^{\infty }_c(\Omega )$ with ${{\,\mathrm{supp}\,}}(b),{{\,\mathrm{supp}\,}}(c)\Subset \Omega$, there exist exterior Dirichlet data $f_1,\dots ,f_{n+1}$ such that the b and c can be uniquely reconstructed from the knowledge of $f_1,\dots ,f_{n+1}$ and $\Lambda _{b,c}(f_1),\dots , \Lambda _{b,c}(f_{n+1})$.

By Runge approximation in $C^{k}$ spaces [see Theorem 1.2 (b)], we have that for any $g\in C^{\infty }(\overline{\Omega })$ there exists a sequence of solutions $\{u_j\}_{j \in \mathbb {N}}$ to the fractional Schrödinger equation (1.4) with drift (and compactly supported coefficients) such that for any $k \in \mathbb {N}$

$$\begin{aligned} \Vert g - d^{-s} u_j \Vert _{C^k(\Omega )} \rightarrow 0 \text { as }j\rightarrow \infty , \end{aligned}$$

where $d(x)= {{\,\mathrm{dist}\,}}(x,\partial \Omega )$ if $x\in \Omega$ is sufficiently close to the boundary of $\Omega$ and d(x) is extended to a positive function smoothly into the interior of $\Omega$. Next, we choose $n+1$ smooth functions $g_1,\dots ,g_{n+1}$ defined in $\Omega$ with the property that

$$\begin{aligned} h(g_1,g_2,\ldots ,g_{n+1})(x) := \det \begin{pmatrix} \partial _1 g_1 &{} \ldots &{} \partial _n g_1 &{} g_1 \\ \vdots &{} \vdots &{} \vdots &{} \vdots \\ \partial _1 g_{n} &{} \ldots &{} \partial _n g_{n} &{} g_{n}\\ \partial _1 g_{n+1} &{} \ldots &{} \partial _n g_{n+1} &{} g_{n+1} \end{pmatrix} (x)\ne 0. \end{aligned}$$

(6.2)

An example for this would be the functions $g_j= x_j$ for $j\in \{1,\dots ,n\}$ and $g_{n+1}=1$. We further set $\widetilde{g}_l= d^{-s} (\chi g_l)$, where $\chi \in C_c^\infty (\Omega )$ and $\chi = 1$ on K, for $l\in \{1,\dots ,n+1\}$ and apply Theorem 1.2 (b). As a consequence, for each $l\in \{1,\dots ,n+1\}$, and in any compact subset $K \Subset \Omega$ there exists a sequence of solutions $\{u_{j,K}^l\}_{j\in \mathbb {N}}$ such that for any $k\in \mathbb {N}$

$$\begin{aligned} \Vert d^{-s}( g_l-u_{j,K}^l)\Vert _{C^k(K)} \rightarrow 0 \text{ as } j \rightarrow \infty . \end{aligned}$$

As $K \subset \Omega$ and as $d(x)>0$ in K, we then also have

$$\begin{aligned} \Vert g_l-u_{j,K}^l\Vert _{C^k(K)} \rightarrow 0 \text{ as } j \rightarrow \infty . \end{aligned}$$

Hence, choosing $j\ge j_0$ large enough and K such that ${{\,\mathrm{supp}\,}}(b)\cup {{\,\mathrm{supp}\,}}(c)\Subset K$, we obtain that

$$\begin{aligned} h( u_{j,K}^1, \ldots , u_{j,K}^{n+1})(x)= \det \begin{pmatrix} \partial _1 u_{j,K}^1 &{} \ldots &{} \partial _n u_{j,K}^1 &{} u_{j,K}^1 \\ \vdots &{} \vdots &{} \vdots &{} \vdots \\ \partial _1 u_{j,K}^n &{} \ldots &{} \partial _n u_{j,K}^n &{} u_{j,K}^n\\ \partial _1 u_{j,K}^{n+1} &{} \ldots &{} \partial _n u_{j,K}^{n+1} &{} u_{j,K}^{n+1} \end{pmatrix} (x)\ne 0. \end{aligned}$$

(6.3)

As a consequence, for these values of $u_{j,K}^l$ the (linear) system for $b_1,\dots ,b_n$ and c

$$\begin{aligned} \begin{pmatrix} \partial _1 u_{j,K}^1 &{} \ldots &{} \partial _n u_{j,K}^1 &{} u_{j,K}^1 \\ \vdots &{} \vdots &{} \vdots &{} \vdots \\ \partial _1 u_{j,K}^n &{} \ldots &{} \partial _n u_{j,K}^n &{} u_{j,K}^n\\ \partial _1 u_{j,K}^{n+1} &{} \ldots &{} \partial _n u_{j,K}^{n+1} &{} u_{j,K}^{n+1} \end{pmatrix} \begin{pmatrix} b_1 \\ \vdots \\ b_n\\ c \end{pmatrix} = -\begin{pmatrix} (-\Delta )^s u_{j,K}^1\\ \vdots \\ (-\Delta )^s u_{j,K}^n \\ (-\Delta )^s u_{j,K}^{n+1} \end{pmatrix} \end{aligned}$$

(6.4)

is solvable (since the $(n+1)\times (n+1)$ matrix in the left hand side of (6.4) is invertible). Thus, from the knowledge of $u_{j,K}^l$, $l\in \{1,\dots ,n+1\}$ for some $j\ge j_0$, it is possible to uniquely recover the drift and the potential simultaneously.

As by the global (nonlocal) unique continuation arguments in [19] (c.f. also Proposition 3.1 from above) it is possible to recover $u_{j,K}^{l}$ given the measurements $f_{j,K}^l$ and $\Lambda _{b,c}(f_{j,K}^l)$ we infer the finite measurement recovery statement.

Finally, in order to infer the openness of the set of possible exterior data, we note that for any $\epsilon >0$ there exists $\delta >0$ such that for any $f=(f^1,\dots , f^{n+1})\in C_c^{\infty }(W)^{n+1}$ with

$$\begin{aligned} \Vert f-f_{j,K}\Vert _{C^{k}_{c}(W)} < \delta , \ k \in \mathbb {N}, \end{aligned}$$

we have by boundedness of the mapping $C_c^k(W)^{n+1} \ni f \mapsto u \in C^{k}(K)$

$$\begin{aligned} \Vert u-u_{j,K}\Vert _{C^k(K)} < \epsilon . \end{aligned}$$

Here $f_{j,K}:=(f_{j,K}^1,\dots ,f_{j,K}^{n+1})\in C_c^{\infty }(W)^{n+1}$ are the exterior data from above, $u=(u^1,\dots ,u^{n+1})$ are the solutions to (1.4) corresponding to the data f and $u_{j,K}:=(u_{j,K}^1,\dots ,u_{j,K}^{n+1})$ are the solutions to (1.4) corresponding to the data $f_{j,K}=(f_{j,K}^1,\dots ,f_{j,K}^{n+1})$. In particular, assuming that

$$\begin{aligned} \det \begin{pmatrix} \partial _1 u_{j,K}^1 &{} \ldots &{} \partial _n u_{j,K}^1 &{} u_{j,K}^1 \\ \vdots &{} \vdots &{} \vdots &{} \vdots \\ \partial _1 u_{j,K}^n &{} \ldots &{} \partial _n u_{j,K}^n &{} u_{j,K}^n\\ \partial _1 u_{j,K}^{n+1} &{} \ldots &{} \partial _n u_{j,K}^{n+1} &{} u_{j,K}^{n+1} \end{pmatrix} \ge c>0 \text{ in } K, \end{aligned}$$

the triangle inequality implies that if $\epsilon >0$ (and correspondingly $\delta$) is chosen sufficiently small, also

$$\begin{aligned} \det \begin{pmatrix} \partial _1 u^1 &{} \ldots &{} \partial _n u^1 &{} u^1 \\ \vdots &{} \vdots &{} \vdots &{} \vdots \\ \partial _1 u^n &{} \ldots &{} \partial _n u^n &{} u^n\\ \partial _1 u^{n+1} &{} \ldots &{} \partial _n u^{n+1} &{} u^{n+1} \end{pmatrix} >0 \text{ in } K. \end{aligned}$$

This concludes the argument. $\square$

Remark 6.3

We conclude this section by some comments on the assumptions of the theorem:

• (a) The compact support condition for the functions $b_j$, $c_j$, $j=1,2$, is assumed here in order to be able to apply the theory of Grubb [24].

• (b) In order to obtain our result we in principle do not need the full strength of the $C^k$, $k\in \mathbb {N}$, approximation result from [20]. It would for instance be sufficient to use a $C^1$ approximation result for which only lower regularity on the coefficients is needed. Since the theory of Grubb is however formulated in the smooth set-up, we do not optimize the regularity dependences here.

• (c) We again emphasize that although the variant of Theorem 1.4 which is proved in this section is interesting from a theoretical point of view, a word of caution is needed as follows: In contrast to the single measurement result from [19], the exterior data $f_1,\dots ,f_{n+1}$ are not arbitrary. In general the specific choice of these data depends on $b_1, b_2$, $c_1,c_2$ and hence the explicit choice of the functions $f_1,\dots , f_{n+1}$ is not known in general. This has a similar background (see the following proof) as many results on hybrid inverse problems, where it is important that constraints are satisfied, see [1, 2]. The result from this section will be improved considerably in the argument leading to the proof of Theorem 1.4.

• (d) In addition to the previous point, there are examples of matrices with entries satisfying elliptic equations, for which the determinant vanishes on an open set, see [10]. This indicates that the zero set of the determinant (6.3) can indeed be large.

6.3 Variations and extensions of the reconstruction result

We discuss a slight variation of the reconstruction result from Sect. 6.2 by relaxing the condition that the fields b, c are compactly supported in $\Omega$. Recall that a $C^\infty$-smooth function f is vanishing to infinite order at a point $x_0$ provided that $\partial ^\alpha f(x_0)=0$ holds for any multi-index $\alpha =(\alpha _1,\ldots ,\alpha _n) \in (\mathbb {N}\cup \{0\})^n$.

Proposition 6.4

Let $\Omega \subset \mathbb {R}^n$ be a bounded domain with a $C^\infty$-smooth boundary. Let $W \subset \Omega _e$ be a non-empty open, smooth domain containing an open neighbourhood of $\partial \Omega$. Let $\frac{1}{2}<s<1$ and assume that $b_j\in C^{\infty }(\overline{\Omega })^n$, $c_j\in C^{\infty }(\overline{\Omega })$ satisfy (1.5) and

$$\begin{aligned} b_1-b_2, \ c_1-c_2 \text{ vanish } \text{ to } \text{ infinite } \text{ order } \text{ on } \partial \Omega . \end{aligned}$$

Then, there exist $n+1$ exterior Dirichlet data $f_1,\dots ,f_{n+1}\in C^{\infty }_c(W)$ such that if

$$\begin{aligned} \Lambda _{b_1,c_1}(f_l) = \Lambda _{b_2,c_2}(f_l) \text{ for } l \in \{1,\dots ,n+1\}, \end{aligned}$$

then $b_1 = b_2$ and $c_1=c_2$ in $\Omega$. Moreover, for any $k\in \mathbb {N}$ the set of exterior data $f_1,\dots ,f_{n+1}$, which satisfies this property forms an open subset in $C^{k}_{ c}(W)$.

This follows from an auxiliary result, which states that under geometric restrictions on the set where we measure the Dirichlet data, we may enlarge our domain and that the DN map on the larger domain is determined by the DN map on the smaller set. This is well-known in the study of local inverse problems – see e.g. [43, Lemma 4.2].

In what follows, we denote by $r_X \widetilde{H}^s(Y)$ the set of all restrictions $f|_{X}$ of functions $f \in \widetilde{H}^s(Y)$ for open $X, Y \subset \mathbb {R}^n$.

Lemma 6.5

Assume that $\Omega$ and $\Omega '$ with $\Omega \subset \Omega '$ are bounded domains with Lipschitz boundaries and $W_1, W_2 \subset \mathbb {R}^n$ are two non-empty open sets, such that $\Omega '{\setminus } \Omega \Subset W_1$ (so, in particular, $W_1 \cap \Omega '_e$ is non-empty). Then let $b_1, b_2 \in W^{1-s, \infty }(\Omega ')$ and $c_1, c_2 \in L^\infty (\Omega ')$, and $\frac{1}{2}< s < 1$, and assume $c_1 = c_2$ and $b_1 = b_2$ in $\Omega ' {\setminus } \Omega$.

Additionally, assume that the coefficients satisfy the eigenvalue condition (1.5) both on $\Omega '$ and on $\Omega$. Assume the equality of DN maps $\Lambda _{b_j,c_j}$ with respect to $\Omega$

$$\begin{aligned} \Lambda _{b_1, c_1}f|_{W_2} = \Lambda _{b_2, c_2}f|_{W_2} \text { for all } f \in r_{W_1 \cap \Omega _e}\widetilde{H}^s(W_1), \end{aligned}$$

Then we have the equality of DN maps $\Lambda '_{b_j,c_j}$ with respect to $\Omega '$:

$$\begin{aligned} \Lambda '_{b_1, c_1}f|_{W_2} = \Lambda '_{b_2, c_2}f|_{W_2} \text { for all }f \in r_{W_1 \cap \Omega '_e}\widetilde{H}^s(W_1). \end{aligned}$$

Proof

Assume $u_1' \in H^s(\mathbb {R}^n)$ solves the Dirichlet problem for $f \in \widetilde{H}^s(W_1 \cap \Omega '_e)$

$$\begin{aligned} ((-\Delta )^s + b_1 \cdot \nabla + c_1)u_1'&= 0 \, \text { in } \, \Omega ',\\ u_1'|_{\Omega '_e}&= f|_{\Omega '_e}. \end{aligned}$$

Then solve the analogous Dirichlet problem with respect to the smaller domain $\Omega$:

$$\begin{aligned} ((-\Delta )^s + b_2 \cdot \nabla + c_2)u_2&= 0 \, \text { in } \, \Omega ,\\ u_2|_{\Omega _e}&= u_1'|_{\Omega _e}. \end{aligned}$$

By using the computation from Remark 2.12 for $\phi \in C_c^\infty (\Omega _e \cap W_2)$ (we need the support condition here) and the hypothesis of equality of the DN maps, we get:

$$\begin{aligned}&\Lambda _{b_1, c_1}f|_{W_2 \cap \Omega _e} = \Lambda _{b_2, c_2}f|_{W_2 \cap \Omega _e}, \\&\quad \implies B_1(u_1', \phi ) = B_2(u_2, \phi ),\\&\quad \implies \int _{\mathbb {R}^n} \phi (-\Delta )^s u_1' dx= \int _{\mathbb {R}^n} \phi (-\Delta )^s u_2 dx, \\&\quad \implies (-\Delta )^s u_1' \equiv (-\Delta )^s u_2 \text { on } W_2 \cap \Omega _e. \end{aligned}$$

Here, for convenience, we have used the abbreviation $B_{j}(\cdot , \cdot ) = B_{b_{j},c_j}(\cdot , \cdot )$ for $j\in \{1,2\}$. Furthermore, note that here we also use that $u_1'|_{\Omega _e} \in r_{\Omega _e \cap W_1} \widetilde{H}^s(W_1)$, since $\Omega ' {\setminus } \Omega \Subset W_1$. Therefore, we have the following relations:

$$\begin{aligned} (-\Delta )^s(u_1' - u_2) = 0 \text { on } W_2 \cap \Omega _e \quad \text {and} \quad u_1' - u_2 = 0 \text { on } W_2 \cap \Omega _e . \end{aligned}$$

(6.5)

By the strong uniqueness Proposition from 3.1, we conclude $u_1' \equiv u_2$ on the whole of $\mathbb {R}^n$.

Now we proceed to compare the DN maps on the domain $\Omega '$. We have the following chain of equalities for $w \in C_c^\infty (W_2 \cup \Omega )$ and $u_1', u_2$ as above:

$$\begin{aligned} B_2'(u_2, w)&= \int _{\mathbb {R}^n} (-\Delta )^{\frac{s}{2}} u_2 \cdot (-\Delta )^{\frac{s}{2}} w dx + \int _{\Omega '} b_2 \cdot \nabla u_2 w dx + \int _{\Omega '} c_2 u_2 w dx\\&= B_2(u_2, w) + \int _{\Omega ' {\setminus } \Omega } b_2 \cdot \nabla u_2 w dx + \int _{\Omega ' {\setminus } \Omega } c_2 u_2 w\\&= B_1(u_1', w) + \int _{\Omega ' {\setminus } \Omega } b_1 \cdot \nabla u_1' w dx + \int _{\Omega ' {\setminus } \Omega } c_1 u_1' w\\&= B_1'(u_1', w). \end{aligned}$$

Here we split the integrals over $\Omega$ and $\Omega ' {\setminus } \Omega$, used the notation $B'_1, B_2'$ for the bilinear form associated to the equation on $\Omega '$, the definition of the DN map (see Definition 2.11), the fact that $u_2 = u_1'$, and $b_1 = b_2$ and $c_1 = c_2$ on $\Omega ' {\setminus } \Omega$.

We conclude that $u_2$ solves $(-\Delta )^{s} u_2 + b_2 \cdot \nabla u_2 + c_2 u_2 = 0$ in $\Omega '$, with $u_2|_{\Omega '_e} = f|_{\Omega '_e}$. Also, we conclude that for all w as above

$$\begin{aligned}(\Lambda '_{b_2, c_2}f, w) = B_2'(u_2, w) = B_1'(u_1', w) = (\Lambda '_{b_1, c_1} f, w).\end{aligned}$$

Here we use $\Lambda '$ notation for the DN map on $\Omega '$. The main claim follows by observing that we may pick arbitrary $w \in C_c^\infty (W_2)$. Note that we need to assume that the operators satisfy condition (1.5) on the bigger set $\Omega '$ to assume well-definedness of the DN maps. $\square$

Sketch of the argument for Proposition6.4 The proof of Proposition 6.4 is similar to Section 6.2, so we only sketch the arguments here. With the auxiliary result of Lemma 6.5 at hand, we deduce the claim of Proposition 6.4 by extending $\Omega$ to $\Omega '$ suitably such that the geometric conditions on the domains are satisfied. As our operator is not self-adjoint, in the extended domain $\Omega '$ we might possibly work with Cauchy data, as we could catch a finite number of Dirichlet eigenvalues (c.f. Remark 6.6 below on how to avoid this in some cases). However, using arguments as in [41], it is possible to deduce analogous results. $\square$

Remark 6.6

(Domain monotonicity) We remark that when the drift term $b=0$, then it is possible to avoid dealing with Cauchy data by using perturbation of domain arguments. Indeed, for the fractional Laplacian (and self-adjoint fractional Schrödinger operators), it is possible to characterize the Dirichlet spectrum through min-max formulations [17] (see also [11] for similar settings for the classical Laplacian). Relying on the weak unique continuation property, it is possible to show the monotonicity of eigenvalues

$$\begin{aligned} \lambda _k(\Omega _1) > \lambda _k(\Omega _2)\text { for }\Omega _1 \subsetneq \Omega _2\text { and }k\in \mathbb {N}, \end{aligned}$$

where $\lambda _k$ is the k-th Dirichlet eigenvalue of the fractional Laplacian. Thus, perturbing the domain suitably and only considering a finite number of eigenvalues in the case of self-adjoint fractional Schrödinger operators, one might work with the DN map instead of having to resort to Cauchy data.

6.4 Proof of Theorem 1.4 and generic properties of determinants via singularity theory

In this section, we prove the full result of Theorem 1.4, i.e. we show that the set of exterior data from which we can choose $n + 1$ measurements in order to recover the coefficients on a compact set K as in previous sections is open and dense. This significantly improves the result from Sect. 6.2 in that the data still depend on the unknown potentials b, c, but in a precise sense they form a large set, i.e. given (random) exterior data, we know that an arbitrarily small perturbation of them will render them admissible in our reconstruction scheme.

The main idea of our argument is to relax the condition that for admissible exterior data $f_1,\dots ,f_{n+1}\in C_c^{\infty }(W)$ we require

$$\begin{aligned} h(P_{b,c}(f_1),\dots ,P_{b,c}(f_{n+1})) \ne 0 \text{ in } K \subset \Omega , \end{aligned}$$

where h was the function from (6.2). Instead, we consider data $f_1,\dots ,f_{n+1}\in C_c^{\infty }(W)$ such that $h(P_{b,c}(f_1),\dots ,P_{b,c}(f_{n+1}))$ is only allowed to vanish to a finite (dimension-dependent) order. Then, by known results from [3, Lemma 3], it follows that the set

$$\begin{aligned} \left\{ x\in K; \ h(P_{b,c}(f_1),\dots ,P_{b,c}(f_{n+1}))(x)=0\right\} \end{aligned}$$

is of measure zero in $K \Subset \Omega$.¹ As a consequence, due to the continuity of b, c, it is then possible to reconstruct both coefficients (see Lemma 6.8). Simultaneously, the set of exterior data $f_1,\dots ,f_{n+1}\in C_c^{\infty }(W)$ for which $h(P_{b,c}(f_1),\dots ,P_{b,c}(f_{n+1}))$ vanishes only of finite (dimension-dependent) order immediately by definition is open in $C_c^{\infty }(W)^{n+1}$. The density of such data will be obtained via small perturbations, relying on ideas of Whitney’s work [53], which had been developed in the context of singularity theory. Technically, this is the most involved part of our arguments.

Let us introduce the set of our admissible exterior conditions.

Definition 6.7

Let b, c satisfy the conditions in Theorem 1.4. Let $n\in \mathbb {N}$ and $k(n)= \Big \lceil {\sqrt{n + 1}}\Big \rceil \in \mathbb {N}$ (i.e. k(n) is the smallest positive integer greater or equal to $\sqrt{n + 1}$). Let $\Omega \subset \mathbb {R}^n$ be as in Theorem 1.4 and $K\Subset \Omega$ be a compact set as in Sect. 6.2. Then, we define the set $\mathcal {F} \subset C_c^{\infty }(W)^{n + 1}$ to be the set of exterior data, such that $(f_1, \cdots , f_{n + 1}) \in {{\,\mathrm{\mathcal {F}}\,}}$, if $h(P_{b,c} f_1, \cdots , P_{b, c} f_{n + 1})$ has at most order of vanishing $k(n)-2$ at each point of $K \Subset \Omega$.

We claim that the set ${{\,\mathrm{\mathcal {F}}\,}}$ yields the desired set of exterior data for the proof of Theorem 1.4. To this end, we first show that given exterior data $(f_1, \cdots , f_{n + 1})\in {{\,\mathrm{\mathcal {F}}\,}}$, it is possible to recover b and c:

Lemma 6.8

Assume that the conditions of Theorem 1.4 hold. Let $(f_1, \cdots , f_{n + 1})\in {{\,\mathrm{\mathcal {F}}\,}}$. Then it is possible reconstruct b, c from the exterior measurements of $f_l$ and $\Lambda _{b,c}(f_l)$, for $l \in \{1,\cdots ,n+1\}$.

Proof

We first recall that the zero set of a smooth function not vanishing to infinite order is contained in a countable union of codimension one submanifolds (see e.g. [3, Lemma 3]). In particular, by definition of the set ${{\,\mathrm{\mathcal {F}}\,}}$, we thus infer that $h(P_{b,c}(f_1),\dots ,P_{b,c}(f_{n+1}))$ vanishes only on a set of Lebesgue measure zero. Let us denote this measure zero set by $B \subset K$. Therefore, the argument in (6.4) goes through on the set $K {\setminus } B$. But b, c satisfy ${{\,\mathrm{supp}\,}}(b)\cup {{\,\mathrm{supp}\,}}(c)\Subset K$ and are smooth, so have unique continuous extensions to K, which can be determined from $b|_{K{\setminus } B}$ and $c|_{K {\setminus } B}$. This concludes the argument for the reconstruction of b, c from $f_l$ and $\Lambda _{b,c}(f_l)$ for $l\in \{1,\cdots ,n+1\}$. $\square$

Hence, it remains to prove the openness and density of the set ${{\,\mathrm{\mathcal {F}}\,}}\subset C^\infty _0(W)^{n+1}$. While the openness is a direct consequence of the definition of the set ${{\,\mathrm{\mathcal {F}}\,}}$, the density of the set ${{\,\mathrm{\mathcal {F}}\,}}$ requires careful arguments. This will be the content of the remaining subsections.

6.4.1 Generic properties of determinants via singularity theory

In order to deduce the density of the set ${{\,\mathrm{\mathcal {F}}\,}}$, we seek to argue by perturbation: The main idea is that for an m-tuple of functions $f = (f_1, \cdots , f_m)$ on $\mathbb {R}^n$, and some differential relation $P(x,D)(f)(x)=0$ on $\mathbb {R}^n$, we may generate a parametric family $f_{\alpha }$, in such a way that on a compact set, near any $\alpha _0$ there is an $\alpha$ arbitrarily close, such that $P(x,D)(f_\alpha ) \ne 0$ on K. Here $\alpha \in \mathbb {R}^N$ for some (large) $N \in \mathbb {N}$ and $f_\alpha$ is given by adding a polynomial of degree r (related to N) to each of the entries $f_i$ of f, with coefficients given by reading off indices of $\alpha$ in a suitable order.

A famous example, due to Morse, of this fact is just a $C^2$ function $f: \mathbb {R}^n \rightarrow \mathbb {R}$. Then by looking at $f_\alpha (x) = f(x) + \langle {\alpha , x}\rangle$, where $\langle {\cdot , \cdot }\rangle$ is the inner product and $\alpha \in \mathbb {R}^n$, by Sard’s theorem the set of $\alpha$ for which $f_\alpha$ has a degenerate critical point is of measure zero. These ideas were generalised by Whitney [53] and others in the area of singularity theory to study generic maps $\mathbb {R}^n \rightarrow \mathbb {R}^m$.

In our case, the idea is that by adding generic polynomials of degree $k(n)\in \mathbb {N}$ with small coefficients to $P_{b, c} f_1, \cdots , P_{b, c} f_{n + 1}$, we may obtain perturbations such that the determinant function h from (6.2) only vanishes of order at most $k(n)-2$. Here $k(n)\in \mathbb {N}$ is the constant from Definition 6.7. Since the perturbation is just by polynomials of order k(n), i.e. by a linear combination of one of

$$\begin{aligned} N_0 = \sum \limits _{j=0}^{k(n)} {n + j - 1 \atopwithdelims ()j} = {n +k(n) \atopwithdelims ()k(n)} \end{aligned}$$

(6.6)

linearly independent polynomials of the form $1, x_i, x_i x_j, ...$, by Runge approximation we may approximate these by $P_{b, c} (f_{i,m})$ for $i\in \{1 , \cdots , N_0\}$ arbitrarily close, for some suitable exterior data $\left\{ f_{i,m}\right\}$’s. By adding a linear combination of $f_{i,m}$ with coefficients $\alpha$ to the exterior data, one can obtain an arbitrary close measurement for which the zero set of h is “good”, i.e. is just given by a stratification of smooth hypersurfaces (see Propositions 6.10 and 6.14 below). In the sequel, we present the details of this argument.

6.4.2 Preliminaries

We need to import some (old) technology from [53, Parts A and B], which allows us to modify functions in a favorable way by simply applying dimension-counting arguments. We consider a mapping $f = (f_1, \cdots , f_m): \mathbb {R}^n \rightarrow \mathbb {R}^m$. We consider derivatives of order up to $r \in \mathbb {N}$ of f and a map $\bar{f}: \mathbb {R}^n \rightarrow \mathbb {R}^N$ given by arranging the partial derivatives of f in some fixed order. Then there is a “bad” set $S \subset \mathbb {R}^N$ that we would like to avoid. In general, the space S is stratified, i.e. there is a splitting $S = \cup _{i \le \mu } S_i$, where $S_i$ are smooth manifolds of dimension $\dim S_i$ for $i\in \{1,\cdots ,\mu \}$. More precisely, we say S is a manifold collection of defect$\delta$, if $\cup _{i \le j} S_i$ is closed for all $j\in \{1,\cdots ,\mu \}$ and ${{\,\mathrm{codim}\,}}S_i \le \delta$ for all i.

We alter f by adding to it a polynomial of degree $\le r$ whose coefficients form a set $\alpha$ of very small numbers. If $f_{\alpha }$ is the resulting mapping $\mathbb {R}^n \rightarrow \mathbb {R}^N$, we may prove that for compact subsets $K \subset \mathbb {R}^n$ and $T \subset S$ there exists $\alpha$ arbitrarily small such that $f_{\alpha }(K) \cap T = \emptyset$. If we fix K, by taking $\bar{f}(K) \subset B_L \subset \mathbb {R}^N$ for some large L, we prove that for any K, there is an $\alpha$ such that $f_{\alpha }(K) \cap S = \emptyset$.

Let $N \in \mathbb {N}$ be the number given as above. Assume that we have a smooth map for $(p, \alpha ) \in \Omega \times R_1 \subset \mathbb {R}^n \times \mathbb {R}^N$

$$\begin{aligned}F(p, \alpha ) = f_\alpha (p) = f_p^*(\alpha )\end{aligned}$$

into $\mathbb {R}^N$. The family $f_{\alpha }$ is called an N-parameter family of mappings of $\Omega$ into $\mathbb {R}^N$ if the matrix $\nabla _\alpha f_p^*(\alpha )$ has full rank for all $p \in \Omega$.

Next, we say a subset Q of a Euclidean space is of finite$\mu$-extent if there is a number A with the following property. For any integer $\kappa$, there are sets $Q_1, \cdots , Q_a$ such that

$$\begin{aligned} Q = Q_1 \cup \cdots \cup Q_a, \quad {{\,\mathrm{diam}\,}}(Q_j) < 1/2^\kappa \ (\text {for all } j),\quad a \le 2^{\mu \kappa }A. \end{aligned}$$

(6.7)

Lemma 6.9

(Lemma 9a in [53]) Let the $f_{\alpha }$ form an N-parameter family of mappings of $\Omega \subset \mathbb {R}^n$ into $\mathbb {R}^N$ and let $\Omega _1\subset \Omega$ and Q be compact subsets of $\mathbb {R}^n$ and of $\mathbb {R}^N$ of finite $\omega$-extent and of finite q-extent, respectively, and suppose $\omega +q < N$. Then for any $\alpha _0\in \mathbb {R}^N$ and for any $\epsilon >0$ there exists $\alpha$ with $|\alpha -\alpha _0|<\epsilon$ such that

$$\begin{aligned} f_\alpha (\Omega _1) \cap Q = 0. \end{aligned}$$

By [53, Section 10], the family of functions $f_\alpha$ given by adding polynomials of order $\le r$ is an N-parameter family of mapping of $\mathbb {R}^n$ into $\mathbb {R}^N$. We call $\delta := N - q$ the defect of S in $\mathbb {R}^N$, so the condition in Lemma 6.9 can be restated as simply $\delta > \omega$. It can be easily seen that the conclusion of the above lemma holds if Q is a manifold collection of defect $\delta > \omega$.

6.4.3 Density argument

Let $\mathcal {F} \subset C_c^{\infty }(W)^{n + 1}$ be the set, which contains exterior measurements from Definition 6.7. By Lemma 6.8, we may reconstruct the drift b and potential c for such exterior measurements.

Notice that $\mathcal {F}$ is an open set, since the set of all functions g with finite order of vanishing at each point is open and the operator $P_{b, c}$ is continuous in given topologies. In the sequel, we seek to prove that $\mathcal {F}$ is also dense.

We first illustrate the argument for the case $n=1$ in which case it is rather transparent. We will then present the proof for the general case below.

Proposition 6.10

Assume $n = 1$. Then $\mathcal {F} \subset C_c^{\infty }(W)^2$ is open and dense.

Proof

Take $f = (f_1, f_2) \in C_c^{\infty }(W)^2$ to be any exterior data and consider $g_i := P_{b, c}f_i \in C^\infty (K)$ for $i = 1, 2$. Write $g = (g_1, g_2)$. We will construct an approximation of $(f_1, f_2)$ in the topology of $C_c^{\infty }(W)^2$ lying in $\mathcal {F}$. In one dimension, we have

$$\begin{aligned} \begin{aligned} h(g_1, g_2) = \begin{pmatrix} g_1' &{} g_1\\ g_2' &{} g_2 \end{pmatrix} = g_1' g_2 - g_1 g_2'. \end{aligned} \end{aligned}$$

(6.8)

In other words, h is the Wronskian of $g_1$ and $g_2$. Now, in this case, we have

$$\begin{aligned} \bar{g}(x) = \big (g_1(x), g_2(x), g_1'(x), g_2'(x), g_1''(x), g_2''(x)\big ) \end{aligned}$$

and $N = 2N_0= 6$ [where $N_0$ was defined in (6.6)]. We consider perturbations of the form

$$\begin{aligned} g_\alpha (x) = g(x) + (\alpha _1 + \alpha _1' x + \alpha _1'' x^2)e_1 + (\alpha _2 + \alpha _2' x + \alpha _2'' x^2)e_2 \in \mathbb {R}^2, \end{aligned}$$

where $e_1, e_2$ are the canonical coordinate vectors in $\mathbb {R}^2$.

We compute the defect of the bad set given by $h=0$ and $\nabla h = 0$, i.e.

$$\begin{aligned} S := \{J_1 = \alpha _1' \alpha _2 - \alpha _1 \alpha _2' = 0\} \cap \{J_2 = \alpha _1'' \alpha _2 - \alpha _1 \alpha _2'' =0 \} \subset \mathbb {R}^6. \end{aligned}$$

We need to show that the defect $\delta = 6 - q$, where q is the extent of S, is bigger that $n = 1$ to apply Lemma 6.9. To this end, we compute the gradients

$$\begin{aligned} \nabla \bar{J}_1 = (-\alpha _2', \alpha _1', \alpha _2, -\alpha _1, 0, 0),\\ \nabla \bar{J}_2 = (-\alpha _2'', \alpha _1'', 0, 0, \alpha _2, -\alpha _1). \end{aligned}$$

These are clearly linearly independent for $\alpha _1 \ne 0$ or $\alpha _2 \ne 0$ or $\det \begin{pmatrix} \alpha _1' &{} \alpha _2'\\ \alpha _1'' &{} \alpha _2'' \end{pmatrix} \ne 0$, so S is a manifold collection of defect $\delta = 2$. So we apply Lemma 6.9 to obtain arbitrarily small values of $\alpha = (\alpha _1, \alpha _1', \alpha _1'', \alpha _2, \alpha _2', \alpha _2'')$ such that $g_\alpha$ satisfies the property that $h(g_\alpha )$ has an empty critical zero set on $K \subset \overline{\Omega }$.

Next, by Runge approximation (see Lemma 6.2), there exists $f_{0, m}, f_{1, m}, f_{2, m} \in C_c^{\infty }(W)$ with $P_{b, c} f_{i, m} \rightarrow x^i$ in $C^\infty (K)$ for $i = 0,1, 2$ and as $m \rightarrow \infty$. Therefore, we define

$$\begin{aligned} f_{\alpha , m} = f + (\alpha _1 f_{0, m} + \alpha _1' f_{1, m} + \alpha _1'' f_{2, m}) e_1 + (\alpha _2 f_{0, m} + \alpha _2' f_{1, m} + \alpha _2'' f_{2, m}) e_2 \in \mathbb {R}^2. \end{aligned}$$

Now fix m large enough, so $P_{b, c} f_{i, m}$ is close to $x^i$, such that the perturbation by elements $f_{i, m}$ of f makes a 6-parameter family of mappings $\Omega ' \rightarrow \mathbb {R}^6$, for $\alpha \in B_1(0) \subset \mathbb {R}^6$ in the unit ball (say), where $K \subset \Omega ' \Subset \Omega$, by compactness. Then we again apply Lemma 6.9 and get that $g_{\alpha ,m} := P_{b, c} f_{\alpha , m} \rightarrow g$ in $C^\infty (K)$ on a sequence of $\alpha$ converging to zero. By construction $f_{\alpha , m} \in \mathcal {F}$, so this finishes the proof. $\square$

Remark 6.11

To prove the desired genericity property, we need to approximate $P_{b, c} f$ by either polynomials or other nice functions (analytic, generic etc.), by use of a linear approximation operator $T_m f$, but such that

$$\begin{aligned} T_m P_{b, c} f = P_{b, c} T_m f \quad \text {and} \quad \lim _{m \rightarrow \infty } T_m f = f. \end{aligned}$$

This is the reason why the usual approximation operators, such is the Bernstein polynomials operator, or the general Weierstrass approximation theorem approach are not good for this purpose. The above approximation argument however proves the existence of such $T_m$, which is obtained by adding a finite linear combination of suitable functions with coefficients going to zero as $m \rightarrow \infty$.

Remark 6.12

There is an alternative proof by hand of the above statement for $n = 1$, not using the Whitney machinery, but only the genericity of Morse functions.

We seek to extend the previous argument to dimension $n=2$ and higher. Unfortunately, in this context, it does not suffice to consider the critical zero set of h, i.e. the set $h=0$ and $\nabla h=0$. The computations below show that we have to include higher order derivatives of h to obtain genericity in the sense of the previous proposition.

Let $N = (n + 1) \times {n + k(n) \atopwithdelims ()k(n)}$ be the number of polynomials of degree $\le k(n)$ with which we perturb (we multiply by $n + 1$ as this is our number of functions). Then we consider the map $\bar{g}:\mathbb {R}^n \rightarrow \mathbb {R}^N$ (see previous subsection) of evaluating the derivatives of order $\le k(n)$ at each point. We define the bad set to be $S = S_{k(n)} \subset \mathbb {R}^N$ consisting of points given by the condition that h and its derivatives up to order $k(n) - 1$ vanish. For a given function $g = (g_1, \cdots , g_{n + 1})$, we denote this set by $Z_0(g) = h^{-1}(0)$, $Z_1(g) = Z_0 \cap (\nabla h)^{-1} (0)$, and inductively we define $Z_j(g) = Z_{j-1}(g) \cap (\nabla ^j h)^{-1}(0)$.

Similarly as in one dimension, we then also have the following lemma, which we prove in the “Appendix”:

Lemma 6.13

(Determinant genericity) The bad set $S_{k(n)} \subset \mathbb {R}^N$ is a manifold collection of defect $n + 1$. In particular, there are arbitrarily small perturbations $g_\alpha$ with $\alpha \in \mathbb {R}^N$, such that $Z_{k(n) - 1}(g_{\alpha }) = \emptyset$ on an arbitrary compact set.

As a corollary, we deduce the main result.

Proposition 6.14

For any $n \in \mathbb {N}$, the set $\mathcal {F} \subset C_c^{\infty }(W)^{n + 1}$ is open and dense.

Proof

The proof is an immediate corollary of Lemma 6.13 and the method of proof of Proposition 6.10. Indeed, openness again follows from the definition of the set ${{\,\mathrm{\mathcal {F}}\,}}$. In order to infer the density of ${{\,\mathrm{\mathcal {F}}\,}}$, we argue along the lines of Proposition 6.10: Let $f=(f_1,\dots ,f_{n+1})\in C_c^{\infty }(W)^{n+1}$ be arbitrary but fixed. By Runge approximation (see Theorem 1.2 (b)), for each $\beta \in \mathbb {N}^{n}$ with $|\beta |\le \lceil {\sqrt{n + 1}}\Big \rceil$ there exists $f_{\beta ,m} \in C_c^{\infty }(W)$ such that

$$\begin{aligned} P_{b,c}(f_{\beta ,m}) \rightarrow x^{\beta } \text{ in } C^{\infty }(K). \end{aligned}$$

As the set of polynomials up to degree $|\beta |\le \Big \lceil {\sqrt{n + 1}\Big \rceil }$ forms a $\nu$-parameter family with $\nu = N_0$ and $N_0$ as in (6.6), for $m\ge m_0$ sufficiently large, also the set $P_{b, c} (f_{\beta ,m})$ with $|\beta | \le \Big \lceil {\sqrt{n + 1}}\Big \rceil$ forms a $N_0$-parameter family. As a result, Lemma 6.9 can be applied to

$$\begin{aligned} P_{b,c}(f_{\alpha }):= P_{b,c}(f) + \sum \limits _{j=1}^{ n + 1} \sum \limits _{\beta \in \mathbb {N}^n, \ |\beta | \le \Big \lceil {\sqrt{n + 1}}\Big \rceil } \alpha _{\beta j} P_{b,c}(f_{\beta , m}) e_j, \end{aligned}$$

where $\{e_1, \dots , e_{n+1}\}$ denotes the canonical basis of $\mathbb {R}^{n + 1}$, $\alpha _{\beta j} \in \mathbb {R}$ and $x^{\beta }:= \prod \nolimits _{\ell =1}^{n} x_{\ell }^{\beta _{\ell }}$. Thus, for any $\epsilon >0$ there exists $\alpha _{\epsilon } \in \mathbb {R}^{N_0 (n +1)}$ with $|\alpha _{\epsilon }|\le \epsilon$ such that $Z_{ k(n) - 1}(P_{b,c}(f_{\alpha _{\epsilon }})) = \emptyset$ on $K \subset \Omega$. By construction, we have $f_{\alpha _{\epsilon }} \in \mathcal {{{\,\mathrm{\mathcal {F}}\,}}}$. This concludes the density proof. $\square$

Combining Lemma 6.8 and Propositions 6.10, 6.14 then implies the result of Theorem 1.4.

Remark 6.15

We remark that Theorem 1.4 together with the results from [19] also yields a constructive reconstruction algorithm for the fractional Calderón problem with drift.

Remark 6.16

Last but not least, we point out that similar openness and density results can also be obtained for the Jacobian by arguing along the same lines as in the “Appendix”. More specifically, genericity results in Lemma 6.13 can be shown to hold with the same critical index $k(n) - 1$, if instead of the determinant in Eq. (6.2) we consider the Jacobian determinant of n functions, which might be of independent interest.

Appendix A. Genericity of determinants

The aim of this appendix is to prove Lemma 6.13. In order to determine the conditions which are imposed by the sets $Z_i$ from above, we first rewrite the condition $\nabla h = 0$ in closed form. Using that $\frac{\partial \det (M)}{\partial M_{ij}} = {{\,\mathrm{Cof}\,}}(M)_{ij}$ for a matrix M, where ${{\,\mathrm{Cof}\,}}(M)$ denotes the matrix of cofactors, we obtain

$$\begin{aligned} \partial _{x_j} h(g_1,\dots ,g_{n+1}) = \frac{\partial h(M)}{\partial M_{kl}} \frac{\partial M_{kl}}{\partial x_{j}} = {{\,\mathrm{Cof}\,}}(M)_{kl} \frac{\partial M_{kl}}{\partial x_{j}}, \end{aligned}$$

where we use summation convention of repeated indices, and

$$\begin{aligned} M=M(g_1,\dots ,g_{n+1}) = \begin{pmatrix} \frac{\partial g_1}{\partial x_1} &{} \ldots &{} \frac{\partial g_1}{\partial x_n} &{} g_1\\ \vdots &{} \ldots &{} \vdots &{} \vdots \\ \frac{\partial g_{n+1}}{\partial x_1} &{} \ldots &{} \frac{\partial g_{n+1}}{\partial x_{n}} &{} g_{n+1} \end{pmatrix} \end{aligned}$$

abbreviates the entries of h. Hence, by the column-wise expansion of the determinant the condition $\partial _{j}h=0$ can be reformulated as

$$\begin{aligned} \sum \limits _{l=1}^{n+1} \det (\widetilde{M}_l^j) = 0 \text{ for } j\in \{1,\dots , n\}, \end{aligned}$$

(A.1)

where

$$\begin{aligned} \widetilde{M}_l^j = M\left( g_1, \dots , g_{l-1}, \frac{\partial g_l}{\partial x_j}, g_{l+1}, \dots , g_{n+1}\right) . \end{aligned}$$

Higher order derivatives can be computed analogously (see “Appendix A.2”).

Hence, in the formal variables $\alpha \in \mathbb {R}^{N_0}$ the two constraints associated with $Z_1$ turn into the formal constraints

$$\begin{aligned} \det (M(\alpha ))&= 0, \end{aligned}$$

(A.2)

$$\begin{aligned} \sum \limits _{l=1}^{n+1} \det (\widetilde{M}_{l}^j(\alpha ))&= 0 \text{ for } j =1,\dots ,n, \end{aligned}$$

(A.3)

where $\alpha \in \mathbb {R}^{N_0}$ denotes the vector of the component mapping and

$$\begin{aligned} M(\alpha ) = \begin{pmatrix} \alpha ^1_1 &{} \alpha _2^1 &{} \ldots &{} \alpha _{n}^1 &{} \alpha ^1\\ \alpha _1^2 &{} \alpha _2^2 &{} \ldots &{} \alpha _n^2 &{} \alpha ^2\\ \vdots &{} \vdots &{} \vdots &{} \vdots &{} \vdots \\ \alpha _1^{n+1} &{} \alpha _2^{n+1} &{} \ldots &{} \alpha _{n}^{n+1} &{} \alpha ^{n+1} \end{pmatrix}. \end{aligned}$$

Here we have used the convention that in the arrangement of the derivatives defining $\alpha \in \mathbb {R}^{N_0}$ in (6.6) we have set $\alpha ^{j}$ to correspond to the function $g_j$, $j\in \{1,\dots ,n+1\}$, and $\alpha ^{j}_k$ corresponding to the partial derivatives $\frac{\partial g_j}{\partial x_k}$, $j\in \{1,\dots ,n+1\}$, $k\in \{1,\dots ,n\}$.

As already indicated above, we note that in general the conditions imposed by $Z_1$ do not suffice to apply Lemma 6.9: Indeed, in order to apply Lemma 6.9 we have to prove that the bad set is of co-dimension $n+1$. When $n \ge 4$ the set of matrices with co-rank equal to two for instance is non-negligible (i.e. it has larger co-dimension than $n+1$; where we recall that the set of $m_1 \times m_2$ matrices of rank r is a manifold of co-dimension equal to $(m_1 - r) (m_2 - r)$). However, if ${{\,\mathrm{codim}\,}}M(\alpha ) = 2$, then the equations in (A.3) are void, i.e. are automatically implied by the corank condition (by definition of the cofactor matrix). Thus, in order to arrive at a stratification of the bad set with a sufficiently large co-dimension, for $n \ge 4$ we must include higher order derivatives of the determinant h; these are encoded in the sets $Z_i$.

The proof of Lemma 6.13 will be an induction on the corank of $M(\alpha )$, which we discuss in the next three sections.

1.1 Corank one case

We analyse the case when the corank of the matrix $M(\alpha )$ is equal to one and prove a stratification into a collection of codimension at least $n + 1$ submanifolds in that case. We will use the notation $\alpha ^i$ to represent the parameter corresponding to the function $f_i$. The derivatives $\partial _{jk} f_i$ correspond to $\alpha ^i_{jk} = \alpha ^i_{kj}$ etc. Assume ${{\,\mathrm{corank}\,}}M(\alpha ) = 1$.

Step 1. Assume ${{\,\mathrm{Cof}\,}}(M)_{11} \ne 0$. By the expansion in (A.3) for $j = 1$ and the Laplace expansion of the determinant, we get an equation for $\alpha ^{1}_{11}$. Since this is the unique equation for $\alpha ^{1}_{11}$ (the relations (A.3) for $j = 2, \cdots , n$ do not contain information on $\alpha _{11}^1$), it is automatically linearly independent of any other relation. We denote the gradient of equation (A.3) for $j = 1$ by $G_1$. Denote by $G_j$ the gradient of equation (A.3) applied to j. We claim that the gradients $G_1, \cdots , G_n$ are linearly independent.

From (A.3) for $j = 2$, we get an equation for $\alpha ^{1}_{12}$ (as ${{\,\mathrm{Cof}\,}}(M)_{11} \ne 0$). The gradient $G_2$ of this equation is linearly independent from $G_1$ since this equation doesn’t contain $\alpha _{11}^1$.

Similarly, from (A.3) for $j = 3$, we get an equation for $\alpha ^{1}_{13}$ (as ${{\,\mathrm{Cof}\,}}(M)_{11} \ne 0$). The gradient $G_3$ of this equation is linearly independent of $G_1$ and $G_2$, since $G_1$ is the only one containing a non-zero component at the place of $\alpha ^1_{11}$. Then, the component of $G_2$ corresponding to $\alpha ^{1}_{13}$ is zero, so we obtain the claimed independence.

Iteratively, we obtain that for any $j=1, \cdots , n$, we can extract $\alpha ^{1}_{1j}$ from (A.3) applied to j, and the so obtained gradient $G_j$ is linearly independent from $G_1, \cdots , G_{j-1}$ by a similar argument as in previous two paragraphs.

It is left to observe that the gradient G of the equation (A.2) is independent of $G_1, \cdots , G_n$, since (A.2) doesn’t contain information about $\alpha ^{1}_{1j}$ for any $j=1,\cdots , n$. Therefore, we are contained in a codimension $n + 1$ set.

Step 2. Assume ${{\,\mathrm{Cof}\,}}(M)_{11} = 0$, but ${{\,\mathrm{Cof}\,}}(M)_{12} \ne 0$. Similarly as before, by (A.3) for $j = 1$ and the Laplace expansion of the determinant, we get an equation for $\alpha ^{1}_{12}$ (as ${{\,\mathrm{Cof}\,}}(M)_{12} \ne 0$). Denote the gradient of this equation by $G_1$. Also, denote by $G_j$ the gradient of (A.3) for j. We claim that $G_1, \cdots , G_n$ are linearly independent.

From (A.3) for $j = 2$, we get an equation for $\alpha ^{1}_{22}$ (as ${{\,\mathrm{Cof}\,}}(M)_{12} \ne 0$). The gradient $G_2$ of this equation is linearly independent of $G_1$, since the component of $G_1$ at the place of $\alpha ^{1}_{22}$ is equal to zero. Similarly, from (A.3) for $j = 3$, we get an equation for $\alpha ^{1}_{23}$ (as ${{\,\mathrm{Cof}\,}}(M)_{12} \ne 0$). The gradient $G_3$ of this equation is linearly independent of $G_1$ and $G_2$. This is because $G_1$ and $G_3$ do not contain information about $\alpha ^1_{22}$; then because the equation corresponding to $G_1$ contains no information about $\alpha ^{1}_{32}$.

Iterating (as before) we obtain linearly independent gradients $G_1, \cdots , G_n$. It is left to observe that the gradient G of equation (A.2) is independent of $G_1, \cdots , G_n$, since (A.2) does not contain information about $\alpha ^{1}_{2,j}$ for any $j=1,\cdots , n$, by looking at the entry $\alpha ^1_2$. Therefore, we are contained in a codimension $n + 1$ set.

Step 3. We assume ${{\,\mathrm{Cof}\,}}(M)_{ij} = 0$ for $i=1,\cdots ,n + 1$ and $j=1,\cdots , n$ (otherwise the arguments are reduced to previous two steps). Then without loss of generality (w.l.o.g.), assume ${{\,\mathrm{Cof}\,}}(M)_{1, n + 1} \ne 0$ (we know that ${{\,\mathrm{Cof}\,}}(M)$ has rank equal to 1 by the condition ${{\,\mathrm{corank}\,}}(M) = 1$). Then from (A.3), we get n independent relations by looking at the coefficient next to $\alpha ^1_{j}$ for $j = 1, \cdots , n$, i.e. the gradients $G_1, \cdots , G_n$ of these equations are linearly independent.

Now we are just left to observe that the components of the gradient G of (A.2) corresponding to $\alpha ^1_i$ are equal to zero by the cofactor condition. Thus G is linearly independent of $G_1, \cdots , G_n$, by looking at the coefficient of $\alpha ^1$.

1.2 Corank two case

Assume now ${{\,\mathrm{corank}\,}}(M) = 2$. We need to include the second order derivatives, unless $n<4$. This is because matrices of ${{\,\mathrm{corank}\,}}M = 2$ have codimension 4, and the equations (A.3) are automatically satisfied in the co-rank two case. The equations for the second order derivatives read

$$\begin{aligned} \sum _{k, l = 1}^{n+1} \det \big (\widetilde{M}_{kl}^{ij}(\alpha )\big ) = 0, \quad \text { for each } \quad i, j =1,\cdots , n. \end{aligned}$$

(A.4)

Here we write $\widetilde{M}^{ij}_{kl}(\alpha )$ for the matrix

$$\begin{aligned} \widetilde{M}^{ij}_{kl} = M\left( g_1, \dots , g_{k-1}, \frac{\partial g_k}{\partial x_i}, g_{k+1}, \dots , g_{l-1}, \frac{\partial g_l}{\partial x_j}, g_{l+1}, \dots , g_{n+1}\right) . \end{aligned}$$

Here w.l.o.g. we may assume $k < l$. If $k = l$, then we get

$$\begin{aligned} \widetilde{M}^{ij}_{kk} = M\left( g_1, \dots , g_{k-1}, \frac{\partial ^2 g_k}{\partial x_i \partial x_j}, g_{k+1}, \dots , g_{n+1}\right) . \end{aligned}$$

Since ${{\,\mathrm{corank}\,}}M(\alpha ) = 2$, it is easy to see that $\det \widetilde{M}^{ij}_{kk}(\alpha ) = 0$ and we only need to consider matrices with first order derivatives (even in the higher corank case) in the Eq. (A.4).

Main argument. W.l.o.g., we assume that the matrix $M_0$ is obtained by erasing the first row and the first column of $M(\alpha )$ has corank equal to one (other cases are dealt with similarly as in Step 2 above). Consider Eq. (A.4) for $i = 1$ and $j=2,\cdots , n$. By the standard cofactor expansion of the determinant, the coefficient next to $\alpha ^1_{11}$ is equal to

$$\begin{aligned} C_j = \sum _{k = 2}^{n+1} \det \widetilde{M}_{1, k}^{j}(\alpha ). \end{aligned}$$

(A.5)

Here the matrix $\widetilde{M}_{1, k}^{j}(\alpha )$ is obtained by erasing the first row, i.e.

$$\begin{aligned} \widetilde{M}_{1, k}^j = M\left( g_2, \dots , g_{k-1}, \frac{\partial g_k}{\partial x_j}, g_{k+1}, \dots , g_{n+1}\right) . \end{aligned}$$

Case 1. Assume we have $C_j = 0$ for each $j=2,\cdots , n$. In this case the n equations in (A.5) together with the condition that $\det (M_0)=0$ are exactly of the form as in the corank one case, which had been discussed in the previous section. Hence, by the inductive hypothesis (i.e. by the ${{\,\mathrm{corank}\,}}M(\alpha ) = 1$ case), we have that gradients of $C_j$, together with the gradient of $\det M_0$ span an n-dimensional space.

Now putting the first row and column of $M(\alpha )$ back to $M_0$, by the corank condition on $M(\alpha )$, we get additionally at least one more independent relation and so this totals to $n + 1$ independent gradients. This means that the co-dimension of the bad set is at least $n + 1$.

Case 2. Assume w.l.o.g. $C_2 \ne 0$. Similarly to Steps 1 and 2 above, we consider the gradients of the Eq. (A.4) for a range of indices $i = 1, \cdots , n$ and $j = 2$. Denote these gradients by $G_1, \cdots , G_n$. We claim that $G_1, \cdots , G_n$ are linearly independent. Now the gradient $G_n$ is non-zero at the component of $\alpha ^1_{n1}$, with the coefficients $C_2$. Note that the component of $\alpha ^1_{n1}$ in $G_2, \cdots , G_{n - 1}$ is equal to zero. Also, the component of $\alpha _{11}^1$ is equal to zero in $G_2, \cdots , G_n$ and so $G_n$ is linearly independent of $G_1, \cdots , G_{n-1}$.

Similarly, for each $k=1,\cdots ,n$, the component of $\alpha ^{1}_{k1}$ in $G_k$ is equal to $C_2$ and the gradients $G_2, \cdots , G_{k - 1}$ have the component of $\alpha ^{1}_{k1}$ equal to zero. But the component of $\alpha _{11}^1$ in $G_2, \cdots , G_k$ is equal to zero. We therefore inductively conclude that $G_k$ is linearly independent of $G_1, \cdots , G_{k - 1}$. This proves the claim.

As a consequence, we obtain n linearly independent conditions from the gradients $G_1,\dots ,G_n$ which put constraints on the components $\alpha _{j1}^1$ with $j =1,\cdots , n$. Combined with the condition from $\det (M(\alpha ))=0$ (which does not involve any component of $\alpha$ with two indeces) this yields that the bad set is at least of co-dimension $n+1$.

1.3 Corank at least three

The argument from the corank two case is now easy to generalise. For the case that corank of M is equal to k, we need to include k derivatives of $\det M$. The extra equations obtained are analogous to (A.4) and now we have a similar equation for each k-tuple of integers $1 \le i_1, i_2, \cdots , i_k \le n$. Again these extra equations only involve $\alpha ^{i}_{jk}$ (and not $\alpha$’s corresponding to higher derivatives), due to the corank condition on $M(\alpha )$.

The cases 1 and 2 above are then obtained in a similar fashion. The case 1 always uses the induction hypothesis in the setting of $n-1$ variables and $k-1$ (first order) derivatives [the equation for $Z_k$ always involves k first order derivatives of $g_1,\dots ,g_{n+1}$ of which we always delete one of the rows with one of the first order derivatives defining quantities analogous to $C_j$ from (A.5)]. Thus, the fact that the bad set is of co-dimension at least one is always a consequence of the co-dimension estimate which follows from the setting with $n-1$ variables and $k-1$ (first order) derivatives (which was proved in the step before) and yields n co-dimensions. Adding the co-dimension from the deletion of the additional row thus gives $n+1$ co-dimensions as desired.

The case 2 is the induction step and is carried out in a very similar way as outlined above, i.e. we use the condition obtained from $Z_k$ involving k first order derivatives of $g_j$. Arguing in this way yields at least $n + 1$ independent relations. We can see that, under the assumption ${{\,\mathrm{corank}\,}}M(\alpha ) = k$, the co-dimension of the bad set should be roughly $n + k^2$, but this is not important for the main argument, as long as we have co-dimension at least $n + 1$.

This concludes the argument for the estimate on the co-dimension of the bad set if one considers $Z_1,\dots ,Z_k$ with

$$\begin{aligned} k = \Big \lceil {\sqrt{n + 1}}\Big \rceil -1, \end{aligned}$$

as for ${{\,\mathrm{corank}\,}}M(\alpha ) \ge \Big \lceil {\sqrt{n + 1}}\Big \rceil$ already by the linear dependence relations we obtain the desired co-dimension $n+1$ condition.