Determining potentials from the scattering map of the time-dependent Schrödinger equation

In this article, we study the problem of how the scattering map and the Poisson operator of the time-dependent Schrödinger equation in curved spaces determine the potential. More precisely, we consider the Schrödinger operator $P$ on $\mathbb{R}^{n+1}_{z,t}$, where $z\in\mathbb{R}^{n}$ with $n\geq 3$, $t\in\mathbb{R}$,

(1.1)

$$P=D_{t}+\Delta_{g(t)}+V(z,t),$$

where $D_{t}=-i\partial_{t}$, $\Delta_{g(t)}$ is the positive Laplace operator with respect to a smooth family of metrics $g(t)$ on $\mathbb{R}^{n}_{z}$, and $V$ is a smooth real-valued potential function ¹¹1The requirement that $V$ is real only enters towards the end of the article in Section 6. We allow complex-valued potentials for results recalled from [8]. Also, we allow complex valued potentials for results in Section 5 determining the potential from the Poisson operator. Even for results in Section 6, we can still determine the real part of the potential when they are complex valued. See Remark 6.3. . In addition, we make the assumption that $g(t)$ is a compactly supported, in spacetime, perturbation of the flat metric:

(1.2)

$$g(t)=g_{0}+\tilde{g}(t),\text{ where }\tilde{g}(t)\text{ is compactly supported in spacetime},$$

and $V$ is also compactly supported in spacetime. Thus, there exist large constants $R$, $T$ such that if either $|z|\geq R$, or if $|t|\geq T$, then $g(t)=\sum_{i,j}g_{ij}(z,t)dz_{i}dz_{j}$ coincides with the flat metric $g_{0}=\sum_{i}dz_{i}^{2}$, and $V$ vanishes identically. We use $K\subset\mathbb{R}^{n+1}$ to denote a compact set containing the support of $\tilde{g}$ and $V$. So outside $K$, $P$ coincides with the ‘free’ Schrödinger operator $P_{0}:=D_{t}+\Delta_{0}$, where $\Delta_{0}$ is the (positive) Laplacian on $\mathbb{R}^{n}$ with the Euclidean metric. We assume, in addition, that each $g(t)$ is non-trapping in the sense that its geodesics always escape any compact region in $\mathbb{R}^{n}$ in finite time.

The scattering map concerns final state data or scattering data of global solutions to $Pu=0$. Every global solution $u$ has an asymptotic expansion of the form

(1.3)

$$u\sim(4\pi it)^{-n/2}e^{i|z|^{2}/4t}f_{\pm}\big(\frac{z}{2t}\big)+O(|t|^{-n/2-\epsilon}),\quad t\to\pm\infty,\,\epsilon>0$$

for large positive or negative times, where this holds in a pointwise sense if $f_{\pm}$ are sufficiently regular, and distributionally in general. The functions $f_{\pm}$ are by definition the asymptotic data of the solution $u$. Then the scattering map of $P$ is the map

(1.4)

$$S:\;f_{-}\to f_{+}.$$

The main theorem of [8] shows that this map is an elliptic Fourier integral operator—more precisely, a Legendre distribution—of a novel type. See [8, Section 1] for more detailed discussion.

The aim of this article is to answer the inverse problem that arises from this. That is, suppose the scattering maps $S_{i}$ associated with two different operators

(1.5)

$$P_{i}=D_{t}+\Delta_{g(t)}+V_{i}(z,t),\;i=1,2$$

coincide up to a certain order, then does this imply $V_{1}=V_{2}$?

Another object of fundamental importance in the scattering theory is the Poisson operator. In terms of the notation above, forward (for $-$) and backward Poisson operators are the maps

(1.6)

$\displaystyle\mathcal{P}_{\pm}:f_{\pm}\to u.$

One can also ask whether $\mathcal{P}_{\pm}$ associated with $P_{i}$ coinciding up to a certain order implies that $V_{1}=V_{2}$.

We give an affirmative answer to both questions when the metric $g(t)$ is in the following class of metrics, which arises naturally in the geometric inverse problems.

See [25, Lemma 2.1] for conditions ensuring that $g(t)$ admits a convex function. In particular, when $g(t)$ is Euclidean, or has non-positive sectional curvature, or has curvature lower-bounded by a constant depending on $R$, then $g(t)$ admits a convex function.

For the first question concerning $S_{1}-S_{2}$, our answer is the following theorem.

We briefly summarize the strategy of the proof here. The common theme in both inverse results is that the compactness assumption is imposed at the critical order at which the first potential-dependent symbol appears. The compactness condition forces this symbol to vanish. For the Poisson operator this first non-trivial symbol contains the integral of $V_{1}-V_{2}$ along bicharacteristics up to a point. While for the scattering map, it contains an integral of the potential difference along the corresponding bicharacteristic. The latter integral is then interpreted as a weighted geodesic X-ray transform.

We first recall, in Sections 2.1 and 2.2, the two pseudodifferential calculi which describe the two ends of the problem. The parabolic scattering calculus is adapted to the Schrödinger operator on spacetime, while the 1-cusp calculus is adapted to the asymptotic data appearing in (1.3). The relation between them is given by the source-sink structure of the rescaled Hamilton flow of the principal symbol of $P$. After blowing up the radial sets $\mathcal{R}_{\pm}$, the corresponding front faces are canonically identified with $\overline{{}^{\mathrm{1c}}T^{*}\mathbb{R}^{n}}$. Thus a one-cusp covector parametrizes a bicharacteristic of the Schrödinger operator which enters and exits the interaction region. At fiber infinity these bicharacteristics remain in a fixed time slice and project to geodesics for the metric $g(t)$.

Section 3 recalls the calculus of 1c-ps Fourier integral operators and Section 4 recalls the calculus of 1c-1c Fourier integral operators. The microlocalized Poisson operators are 1c-ps Fourier integral operators and the microlocalized scattering map is a 1c-1c Fourier integral operator.

We prove Theorem 1.5 in Section 5. Let $P_{i}=D_{t}+\Delta_{g(t)}+V_{i}$ and write $a_{\pm,j}^{(i)}$ for the symbols of the Poisson parametrices associated to $P_{i}$. Let $H_{p}^{2,0}$ be the rescaled Hamilton vector field associated to either of the $P_{i}$, which coincide at the leading order. The principal symbols of them are obtained by solving transport equations along the $H_{p}^{2,0}$-flow. The leading symbols $a_{\pm,0}^{(i)}$ are independent of the potential.

The first contribution of the potential difference occurs in

(1.8)

$\displaystyle b_{\pm,1}:=a_{\pm,1}^{(1)}-a_{\pm,1}^{(2)},$

which again satisfies a transport equation along the $H_{p}^{2,0}$ with $V_{1}-V_{2}$ appearing in the forcing term. We will then show $V_{1}-V_{2}=0$ by showing $b_{\pm,1}=0$.

The proof of Theorem 1.2 in Section 6 is more involved. We will use the following identity relating the scattering map to Poisson operators:

(1.9)

$\displaystyle S_{i}Q_{\mathrm{1c}}=i(2\pi)^{n}(\mathcal{P}_{+}^{(i)})^{*}[P_{i},Q_{+}]\mathcal{P}_{-}^{(i)}Q_{\mathrm{1c}},$

where $Q_{+},Q_{\mathrm{1c}}$ are suitable microlocal cutoffs. Using this, we will show that the sub-principal symbol of $(S_{1}-S_{2})Q_{\mathrm{1c}}$ is given by an integral along the bicharacteristic line of $V_{1}-V_{2}$ times a certain weight. We will then show that, with carefully chosen microlocal cut-offs, the weight is non-vanishing on the support of $V_{1}-V_{2}$. Then the compactness assumption in Theorem 1.2 will enforce this sub-principal symbol to vanish, which in turn shows $V_{1}-V_{2}=0$.

We conclude this subsection by discussing how an explicit reconstruction of the potential can be achieved using our method. We begin with the Poisson operator, which is relatively more direct. Using (5.26)(5.27) with $V_{1}=V$, $V_{2}=0$, we know that from the principal symbol of $\mathcal{P}_{\pm}^{(1)}-\mathcal{P}_{\pm}^{(2)}$, one obtains a weighted integral of $V$ along a bicharacteristic line up to a certain point. Then differentiating this integral along the rescaled Hamilton vector field recovers $V$ times the weight. In addition, the weight is given by (5.16), which in principle can be computed explicitly in terms of the metric $g$ and this recovers $V$.

If one uses the data from the scattering map, then as sketched above, applying our method to $V_{2}=0,V_{1}=V$, the principal symbol of $S_{1}-S_{2}$ gives the weighted X-ray transform of the potential. Then as pointed out in [28, Section 2], one can construct the inverse of (the normal operator of) the X-ray transform using a Neumann series, which recovers the potential in our setting.

Recovering a potential from scattering data is a classical problem in scattering theory and we can only give a brief review on a limited part of it. On the conceptual level, this is an inverse problem asking whether the interaction can be reconstructed from the observed incoming and outgoing data.

The classical way to treat the scattering map is to decompose it with respect to energy. In one dimension, the Gel’fand–Levitan–Marchenko theory (see e.g. [5, 1]) shows that the reflection/transmission data together with bound-state information determine the potential; see [2] and references therein. In higher dimensions the fixed-energy problem is more subtle. The works of Ramm, Nachman and Novikov [26, 23, 24] showed that, under compact-support or exponential decay hypotheses, a time-independent potential can be uniquely determined and in principle reconstructed from a fixed-energy scattering matrix or amplitude. See the lecture notes of Melrose [20] for the geometric microlocal analysis approach to this problem, which is closer to this article. In particular, in [20, Section 3.6], Melrose discussed how the leading order scattering amplitude associated to the Helmholtz operator recovers the X-ray transform of the potential, which in turn determines the potential. For time-dependent potentials, see the works of Kitada and Yajima [17][18]. See also the work of Enss-Weder [4] for the problem in the N-body setting, and Weder [32] in such a setting with a time-dependent potential. For a similar problem for the wave equation, see the works of Stefanov [27] and Uhlmann [29] and references therein.

The main advantage of our microlocal framework, compared with more classical scattering theory, is to handle curved spaces equipped with a time dependent metric and time-dependent potentials, on which the literature remains relatively sparse. In the time-independent setting, this is treated by Joshi and Joshi-Sá Berreto in a series of works [14, 15, 13, 16] using microlocal tools developed by Melrose-Zworski [21].

For a similar inverse problem determining the metric instead of the potential, we refer readers to the discussion in [12, Section 1.3]. The calculus of Fourier integral operators we are using was introduced by Hassell and the author in [8]. It borrows ideas from Hörmander [10], Duistermaat-Hörmander [3], Melrose-Zworski [21], Vasy [31], Hassell-Wunsch [9].

In this section we give a brief introduction of the parabolic scattering pseudodifferential algebra. As the name indicates, it is the ‘parabolic version’ of the scattering pseudodifferential algebra, where ‘parabolic’ refers to the way to compactify the fiber infinity. More precisely, let $(t,z,\tau,\zeta)$ be coordinates on $T^{*}\mathbb{R}^{n+1}$, with $\tau dt+\zeta\cdot dz$ being the canonical form. Then the (compactified) parabolic scattering cotangent bundle, denoted by $\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$, is defined by compactifying $T^{*}\mathbb{R}^{n+1}$ in the following way. The ‘base’ $\mathbb{R}^{n+1}$ is compactified radially to be a ball with boundary defining function $x_{\mathrm{ps}}=(1+t^{2}+|z|^{2})^{-1/2}$, while each fiber is compactified ‘parabolically’ to be a ball with boundary defining function

(2.1)

$\displaystyle\rho_{\mathrm{ps}}=(1+\tau^{2}+|\zeta|^{4})^{-1/4},$

with $\mathrm{ps}$ standing for ‘parabolic scattering’.

Then the corresponding parabolic scattering pseudodifferential operators $\Psi_{\mathrm{ps}}^{m,l}(\mathbb{R}^{n+1})$ are operators that are quantizations of symbols in $S_{\mathrm{ps}}^{m,l}(\mathbb{R}^{n+1})$, which means they have Schwartz kernels of the form

$\displaystyle\mathrm{Op}(a)=(2\pi)^{-(n+1)}\int e^{i(t-t^{\prime})\tau+(z-z^{\prime})\cdot\zeta}a(t,z,\tau,\zeta)d\zeta d\tau,$

in the distributional sense.

The principal symbol of $A=\mathrm{Op}(a)\in\Psi_{\mathrm{ps}}^{m,l}(\mathbb{R}^{n+1})$ is defined to be the equivalence class of $a$ in $S_{\mathrm{ps}}^{m,l}(\mathbb{R}^{n+1})$ quotient by $S_{\mathrm{ps}}^{m-1,l-1}(\mathbb{R}^{n+1})$:

$\displaystyle[a]\in S_{\mathrm{ps}}^{m,l}(\mathbb{R}^{n+1})/S_{\mathrm{ps}}^{m-1,l-1}(\mathbb{R}^{n+1}).$

A symbol $a\in S_{\mathrm{ps}}^{m,l}(\mathbb{R}^{n+1})$ (and corresponding operator $A=\mathrm{Op}(a)$) is said to be elliptic if it satisfies:

$\displaystyle|\rho_{\mathrm{ps}}^{m}x_{\mathrm{ps}}^{l}a|\geq C,\text{ when }\rho_{\mathrm{ps}}\leq\epsilon\text{ or }x_{\mathrm{ps}}\leq\epsilon,$

for some $\epsilon>0,C>0$. And we say that $a$ (and corresponding operator $A=\mathrm{Op}(a)$) is elliptic at $q\in\partial\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$ if there is a neighborhood of $q$ on which the above inequality is satisfied. And the set of all such $q\in\partial\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$ is denoted by $\mathrm{Ell}_{\mathrm{ps}}^{m,l}(a)$ (or $\mathrm{Ell}_{\mathrm{ps}}^{m,l}(A)$).

Finally we recall the concept of (parabolic) wavefront sets, which is also called the micro-support. For $A=\mathrm{Op}(a)$, its parabolic scattering operator wavefront set $WF^{\prime}_{\mathrm{ps}}(A)$ is defined, as a subset of $\partial\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$, as follows. For $q\in\partial\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$, we say $q\notin\mathrm{WF}^{\prime}_{\mathrm{ps}}(A)$ if and only if there is $\chi\in C^{\infty}(\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}})$ with $\chi(q)=1$ such that $\chi a\in\mathcal{S}(\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}})$. In particular, $\mathrm{WF}^{\prime}_{\mathrm{ps}}(A)=\emptyset$ if and only if $A\in\Psi_{\mathrm{ps}}^{-\infty,-\infty}(\mathbb{R}^{n+1})$.

In this section, we briefly introduce the 1-cusp pseudodifferential algebra, and refer readers to [33][30, Section 2][11, Section 2] for more details. For its connection to time-dependent Schrödinger equations, see [8, Section 3.3].

Let $M$ be an $m-$dimensional manifold with boundary, with boundary defining function $x_{\mathrm{1c}}$, and suppose $y=(y_{1},..,y_{m-1})$ is a coordinate system of $\partial M$, which together with $x_{\mathrm{1c}}$ forms a local coordinate system of $M$ near the boundary. Then the space of 1-cusp vector fields, denoted by $\mathcal{V}_{\mathrm{1c}}$, is locally spanned over $C^{\infty}(M)$ by

(2.2)

$\displaystyle x_{\mathrm{1c}}^{3}\partial_{x_{\mathrm{1c}}},x_{\mathrm{1c}}\partial_{y_{j}},j=1,2,...,m-1.$

$\mathcal{V}_{\mathrm{1c}}$ gives rise to a vector bundle with (2.2) being its local frame, which is called 1-cusp tangent bundle, and denoted by ${}^{\mathrm{1c}}TM$.

The class of 1-cusp differential operators of order (at most) $k$, denoted by $\mathrm{Diff}_{\mathrm{1c}}^{k}(M)$, consists of polynomials of these vector fields of degree at most $k$:

(2.3)

$\displaystyle\begin{split}\mathrm{Diff}_{\mathrm{1c}}^{k}(M)=\{\sum_{\alpha+|\beta|\leq k}a_{\alpha\beta}(x_{\mathrm{1c}},y)(x_{\mathrm{1c}}^{3}\partial_{x_{\mathrm{1c}}})^{\alpha}(x_{\mathrm{1c}}\partial_{y})^{\beta}:\\ a_{\alpha\beta}\in C^{\infty}(M),\alpha\in\mathbb{N},\beta\in\mathbb{N}^{m-1}\}.\end{split}$

The 1-cusp cotangent bundle, denoted by ${}^{\mathrm{1c}}T^{*}M$ is the dual bundle of ${}^{\mathrm{1c}}TM$. It is locally spanned over $C^{\infty}(M)$ by

(2.4)

$\displaystyle\frac{dx_{\mathrm{1c}}}{x_{\mathrm{1c}}^{3}},\frac{dy_{j}}{x_{\mathrm{1c}}},j=1,2,...,m-1.$

Then $T^{*}M$ embeds into ${}^{\mathrm{1c}}T^{*}M$ canonically in the interior of $M$, giving it a symplectic structure naturally. In particular, one may write the canonical one form as

(2.5)

$\displaystyle\xi_{\mathrm{1c}}\frac{dx_{\mathrm{1c}}}{x_{\mathrm{1c}}^{3}}+\eta_{\mathrm{1c}}\cdot\frac{dy}{x_{\mathrm{1c}}},$

where $(\xi_{\mathrm{1c}},\eta_{\mathrm{1c}})$ are coordinates of fibers on ${}^{\mathrm{1c}}T^{*}M$. We use $\overline{{}^{\mathrm{1c}}T^{*}M}$ to denote the compactification of ${}^{\mathrm{1c}}T^{*}M$, obtained by compactifying each fiber radially to be a ball with boundary defining function

(2.6)

$\displaystyle\rho_{\mathrm{1c}}=(1+\xi_{\mathrm{1c}}^{2}+|\eta_{\mathrm{1c}}|^{2})^{-1/2}.$

Notice the minor difference with $\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$, which is compactified on both base and fiber level, is because the 1-cusp construction is happening on the manifold with boundary $M$, which is already ‘compactified a priori’.

The 1-cusp symbol class of differential order $m$ and decay (in fact growth) order $l$, denoted by $S_{\mathrm{1c}}^{m,l}(M)$, is defined to be smooth functions on ${}^{\mathrm{1c}}T^{*}M$ satisfying that for any $j,\alpha,k,\beta$ there exists a constant $C_{j\alpha k\beta}$ such that

(2.7)

$\displaystyle|x_{\mathrm{1c}}^{l}\rho_{\mathrm{1c}}^{m-k-|\beta|}(x_{\mathrm{1c}}\partial_{x_{\mathrm{1c}}})^{j}\partial_{y}^{\alpha}\partial_{\xi_{\mathrm{1c}}}^{k}\partial_{\eta_{\mathrm{1c}}}^{\beta}a(x_{\mathrm{1c}},y,\xi_{\mathrm{1c}},\eta_{\mathrm{1c}})|\leq C_{j\alpha k\beta}.$

And locally they quantize to be operators acting by

(2.8)

$\displaystyle\mathrm{Op}(a)u(x_{\mathrm{1c}},y)=(2\pi)^{-n}\int e^{i\xi_{\mathrm{1c}}\frac{x_{\mathrm{1c}}-x_{\mathrm{1c}}^{\prime}}{x_{\mathrm{1c}}^{3}}+\eta_{\mathrm{1c}}\cdot\frac{y-y^{\prime}}{x_{\mathrm{1c}}}}a(x_{\mathrm{1c}},y,\xi_{\mathrm{1c}},\eta_{\mathrm{1c}})u(x_{\mathrm{1c}}^{\prime},y^{\prime})d\xi_{\mathrm{1c}}d\eta_{\mathrm{1c}}\frac{dx_{\mathrm{1c}}^{\prime}dy^{\prime}}{(x_{\mathrm{1c}}^{\prime})^{n+2}}.$

The collection of all such operators is called 1-cusp pseudodifferential operators with differential order $m$ and decay order $l$, and denoted by $\Psi_{\mathrm{1c}}^{m,l}(M)$. The operator wavefront set $\mathrm{WF}_{\mathrm{1c}}^{\prime}(A)$ is the part of $\partial\overline{{}^{\mathrm{1c}}T^{*}M}$ defined by

(2.9)

$\displaystyle\begin{split}\mathrm{WF}^{\prime}_{\mathrm{1c}}(A)=\Big(\{q\in\partial\big(\overline{{}^{\mathrm{1c}}T^{*}M})\;|\text{There exists }\chi\in C^{\infty}(\overline{{}^{\mathrm{1c}}T^{*}M})\text{ such that }\chi a\in\mathcal{S}(\overline{{}^{\mathrm{1c}}T^{*}M}),\,\chi(q)=1\}\Big)^{\complement}.\end{split}$

Next we define the ellipticity of symbols and operators.

Under the elliptic condition in Definition 2.2, see [30, Section 2.5], its quantization $A$ has a parametrix $B\in\Psi_{\mathrm{1c}}^{-m,-l}(M)$ such that

$\displaystyle AB-\mathrm{Id},\;BA-\mathrm{Id}\in\Psi_{\mathrm{1c}}^{-\infty,-\infty}(M).$

One can now define 1-cusp Sobolev spaces $H^{s,r}_{\mathrm{1c}}(M)$ (see [30, Section 2.5]) for $s\geq 0$ by choosing $A\in\Psi_{\mathrm{1c}}^{s,0}(M)$ elliptic, and demanding

(2.10)

$$u\in H^{s,r}_{\mathrm{1c}}(M)\Leftrightarrow u\in x_{\mathrm{1c}}^{r}L_{\mathrm{1c}}^{2}(M)\ \text{and}\ Au\in x_{\mathrm{1c}}^{r}L_{\mathrm{1c}}^{2}(M);$$

here $L_{\mathrm{1c}}^{2}(M)$ is the $L^{2}$ space relative to the 1-cusp density $\frac{dx\,dy}{x_{\mathrm{1c}}^{n+2}}$. They are equipped with norms:

$$\|u\|^{2}_{H^{s,r}_{\mathrm{1c},h}}=\|x_{\mathrm{1c}}^{-r}u\|_{L_{\mathrm{1c}}^{2}}^{2}+\sum_{j+|\alpha|\leq s}\|(hx_{\mathrm{1c}}^{3}D_{x_{\mathrm{1c}}})^{j}(h^{1/2}x_{\mathrm{1c}}D_{y})^{\alpha}u\|_{L_{\mathrm{1c}}^{2}}^{2}.$$

The 1-cusp Sobolev spaces for other $s$ are defined via interpolation and duality. Then 1-cusp pseudodifferential operators are bounded on these Sobolev spaces, namely for $A\in\Psi_{\mathrm{1c}}^{m,l}(M)$ and all $s,r$, $A$ is a bounded linear operator from $H^{s,r}_{\mathrm{1c}}$ to $H^{s-m,r-l}_{\mathrm{1c}}$.

We recall the phase space dynamics associated with the time-dependent Schrödinger equation and discuss how $\overline{{}^{\mathrm{1c}}T^{*}\mathbb{R}^{n}}$ parametrizes its bicharacteristic lines.

Let $p$ be the left symbol of $P$, the interesting part is where the Schrödinger operator is non-elliptic, i.e., the characteristic variety:

(3.1)

$$\Sigma(P):=p^{-1}(0)\subset\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}.$$

The set

(3.2)

$$\mathrm{Char}(P)=\Sigma(P)\cap\partial(\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}})$$

is where the microlocal propagation takes place.

One of the key ingredients of our analysis is the global source-sink structure of the rescaled Hamilton flow associated to $P$, which is denoted by $H_{p}^{2,0}$ in (3.4) below. Its integral curves (more precisely, their closures) are called bicharacteristic lines. As the discussion in [6, Section 3] shows, the bicharacteristic lines starting on the boundary with finite frequency will remain on the boundary for all time and only those bicharacteristics at fiber infinity of $\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$ will leave $\partial\overline{\mathbb{R}^{n+1}}$ and meet the metric perturbation, so we will restrict ourselves to be near the characteristic set in such a region. In addition, such a bicharacteristic line will stay in the same time slice, or in other words, it reaches the endpoint ‘instantly’ in terms of $t$. In this region, we know $|\tau|^{1/2}$ and $|\zeta|_{g}$ are comparable and we can use ${\rho_{\mathrm{ps}}}=\big(\sum_{j,k}g^{jk}(z,t)\zeta_{j}\zeta_{k}\big)^{-1/2}$.

Let $\hat{\tau}=\rho_{\mathrm{ps}}^{2}\tau\in\mathbb{R}$, and $\hat{\zeta}\in\mathbb{R}^{n-1}$ be coordinates parametrizing $\rho_{\mathrm{ps}}\zeta$ in the sphere with respect to $g$. A valid coordinate system with $(t,z)$ in a compact region is

(3.3)

$$(t,z,\rho_{\mathrm{ps}},\hat{\tau},\hat{\zeta}).$$

Let $\rho_{\mathrm{base}}=(1+t^{2}+|z|^{2})^{-1/2}$ be the defining function of the spacetime infinity, then the rescaled Hamilton vector field of $p$

(3.4)

$$H_{p}^{2,0}={\rho_{\mathrm{ps}}}\rho_{\mathrm{base}}^{-1}H_{p}$$

is a smooth vector field on $\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$ that is tangent to its boundary. The flow of $H_{p}^{2,0}$ has a global source-sink structure with the following ‘radial sets’ being the source and sink.

Since the flow of $H^{2,0}_{p}$ has a source-sink structure, we have the decomposition of $\mathcal{R}$:

(3.6)

$$\mathcal{R}=\mathcal{R}_{-}\cup\mathcal{R}_{+}\subset\overline{{}^{\mathrm{ps}}T_{\partial\overline{\mathbb{R}^{n+1}}}^{*}\mathbb{R}^{n+1}},$$

where $\mathcal{R}_{-}$ is the source and $\mathcal{R}_{+}$ is the sink. In the region $|t/z|\leq C$ for a constant $C>0$, they are given by

(3.7)

$$\mathcal{R}_{\pm}\supset\{x_{\mathrm{ps}}=0,\ t/|z|=\pm\rho_{\mathrm{ps}}/2,\ \hat{\zeta}/|\zeta|=\pm z/|z|,\ \tau/|\zeta|^{2}=-1\}\cap\{\frac{|t|}{|z|}\leq C\}.$$

So one can see that $\mathcal{R}_{-}$ is a graph over the southern hemisphere (i.e., the part of $\partial\overline{\mathbb{R}^{n+1}}$ with $t/\langle z\rangle\leq 0$) and $\mathcal{R}_{+}$ is a graph over the northern hemisphere (i.e., the part of $\partial\overline{\mathbb{R}^{n+1}}$ with $t/\langle z\rangle\geq 0$). Both of them turn vertical (in the sense of tending to fiber infinity) when approaching the equator $\{t/|z|=0,1/|z|=0\}$. Though they overlap at the equator on the base level, they tend to the opposite direction in terms of the frequency, hence they remain disjoint on the phase space level.

As aforementioned, for each point on $\partial\overline{\mathbb{R}^{n+1}}$, there is an $n$-dimensional family of geodesics tending to it. On a phase space level, this corresponds to different bicharacteristic lines tending to or emanating from the same point in $\mathcal{R}_{\pm}$. Then they are distinguished via blowing up $\mathcal{R}_{\pm}$ within $\Sigma$:

(3.8)

$$[\Sigma;\mathcal{R}_{\pm}].$$

We denote the front face created by the blow up by $W_{\pm}$, then $W_{\pm}$ precisely parametrizes all bicharacteristic lines. Using $W_{\pm}$, we have the following correspondence between a point on a bicharacteristic line and the ‘endpoint’ of this bicharacteristic line.

Since we have a canonical identification between $\mathcal{R}_{\pm}$ and $\overline{\mathbb{R}^{n}}$ via projection to the base, $\overline{{}^{\mathrm{1c}}T^{*}\mathcal{R}_{\pm}}$ above can be replaced by $\overline{{}^{\mathrm{1c}}T^{*}\mathbb{R}^{n}}$. See [7, Section 3.3] for a more explicit characterization of this correspondence.

The phase space for 1c-ps Lagrangian distributions (a class of operators that includes suitably microlocalized Poisson operators) is obtained by blowing up the corner of

(3.9)

$\displaystyle\mathcal{M}_{0}=\overline{{}^{\mathrm{ps}}{T^{*}\mathbb{R}^{n+1}}}\times\overline{{}^{\mathrm{1c}}{T^{*}\mathbb{R}^{n}}}$

at base infinity of $\overline{{}^{\mathrm{1c}}{T^{*}\mathbb{R}^{n}}}$ and fiber-infinity of $\overline{{}^{\mathrm{ps}}{T^{*}\mathbb{R}^{n+1}}}$. Here the $\mathbb{R}^{n}$ factor represents the interior of either $\mathcal{R}_{+}$ or $\mathcal{R}_{-}$. We refer readers to [22][19] for more details about blow ups. Concretely, we define

(3.10)

$\displaystyle\mathcal{M}:=[\mathcal{M}_{0};\{\rho_{\mathrm{ps}}=0,x_{\mathrm{1c}}=0\}],\quad{\rho_{\mathrm{ps}}}=\big(\sum_{j,k}g^{jk}(z,t)\zeta_{j}\zeta_{k}\big)^{-1/2},$

and denote the blow down map by

(3.11)

$\displaystyle\beta_{\mathrm{1c}-\mathrm{ps}}:\mathcal{M}\rightarrow\mathcal{M}_{0}.$

The front face created by this blow-up shall be denoted ${\mathrm{ff}_{\mathrm{1c}-\mathrm{ps}}}$. The new smooth coordinate on ${\mathrm{ff}_{\mathrm{1c}-\mathrm{ps}}}$ introduced by the blow up is

(3.12)

$$\varsigma=\frac{x_{\mathrm{1c}}}{\rho_{\mathrm{ps}}},$$

or its reciprocal. In the interior of ${\mathrm{ff}_{\mathrm{1c}-\mathrm{ps}}}$, either ${x_{\mathrm{1c}}}$ or ${\rho_{\mathrm{ps}}}$ can be taken as a boundary defining function.

The manifolds $\mathcal{M}_{0}$ and $\mathcal{M}$ have codimension 4 corners. However, we shall only be interested in a neighbourhood of a compact subset $K$ of the interior of ${\mathrm{ff}_{\mathrm{1c}-\mathrm{ps}}}$; in particular, we shall stay away from all the other boundary hypersurfaces. So, in effect, we are dealing with a manifold with boundary.

The manifold $\mathcal{M}$ is endowed with a canonical symplectic structure²²2In this article, we allow symplectic structures to blow up or degenerate at the boundary $\omega$ from $\mathcal{M}_{0}$ by lifting the symplectic form on $\mathcal{M}_{0}$, which in turn is equipped with the product symplectic structure from its two factors. We are particularly interested in the symplectic/contact structures on ${\mathrm{ff}_{\mathrm{1c}-\mathrm{ps}}}$ induced by this symplectic structure in the interior. To prepare for this, we shall specify coordinates to use in a neighbourhood of $K\subset{\mathrm{ff}_{\mathrm{1c}-\mathrm{ps}}}$. These will be ${x_{\mathrm{1c}}}$ (a boundary defining function) and ${y_{\mathrm{1c}}}$, which are base coordinates on $\mathbb{R}^{n}$ near base infinity; ${\xi_{\mathrm{1c}}}$, ${\eta_{\mathrm{1c}}}$, their one-cusp dual coordinates as defined in Section 2.2; $z$, $t$, Euclidean space and time coordinates; $\varsigma$ as defined in (3.12); and fibre coordinates near fibre-infinity, which we take to be $\tilde{\tau}=\tau{x_{\mathrm{1c}}}^{2}$, $\hat{\zeta}=\zeta/|\zeta|$. We remark that we are also mostly interested in $(z,t)$ near the perturbation of the metric, which by assumption is a compact set in spacetime. On the other hand, $(\zeta,\tau)$ will be near infinity, since ${\rho_{\mathrm{ps}}}=0$ at ${\mathrm{ff}_{\mathrm{1c}-\mathrm{ps}}}$. To summarize, our coordinates are

(3.13)

$${x_{\mathrm{1c}}},\quad{y_{\mathrm{1c}}},\quad{\xi_{\mathrm{1c}}},\quad{\eta_{\mathrm{1c}}},\quad z,\quad t,\quad\varsigma,\quad\tilde{\tau},\quad\hat{\zeta}.$$

See [8, Proposition 4.4] for the reason for the term ‘fibred’ from a symplectic fibration. The 1c-ps Lagrangian submanifolds that are used for our analysis are the twisted forward and backward sojourn relations which arise from the bulk-boundary duality we discussed in Section 3.1. As the calculation in [8, Section 4.4] shows, the rescaled Hamilton vector field $H_{p}^{2,0}$ is a smooth vector field tangent to the boundary of $[\Sigma;\mathcal{R}]$ except at $W_{\pm}$, where it is transverse: inward-pointing at $W_{-}$ and outward-pointing at $W_{+}$. Using the non-trapping assumption we see that each point of $W_{\pm}$ gives rise to a smooth integral curve of $H_{p}^{2,0}$, that travels from $W_{-}$ to $W_{+}$ in finite parameter time $s$. Let $q_{-}$ be a point of $W_{-}$, and let $\gamma_{q_{-}}(s)$ be the integral curve of $H_{p}^{2,0}$ emanating from $q_{-}$ at time $s=0$, and arriving at $W_{+}$ at time $T(q_{-})>0$. Similarly, let $q_{+}$ be a point of $W_{+}$, and let $\mu_{q_{+}}(s)$ be the integral curve of $H_{p}^{2,0}$ emanating from $q_{+}$ at time $s=0$, and arriving at time $T^{\prime}(q_{+})<0$. The forward sojourn relation is the subset of $\mathcal{M}_{0}=\overline{{}^{\mathrm{ps}}{T^{*}\mathbb{R}^{n+1}}}\times\overline{{}^{\mathrm{1c}}T^{*}{\mathcal{R}_{-}}}$ defined by

(3.14)

$$\overline{\Lambda_{-}^{\prime}}=\{(\gamma_{q_{-}}(s),q_{-})\mid q_{-}\in W_{-},\ s\in[0,T(q_{-})]\}.$$

Similarly, the backward sojourn relation is defined by

(3.15)

$$\overline{\Lambda_{+}^{\prime}}=\{(\mu_{q_{+}}(s),q_{+})\mid q_{+}\in W_{+},\ s\in[T^{\prime}(q_{+}),0]\}.$$

Then our Lagrangian will be a lifted and microlocalized version of them, after twisting the sign of the 1-cusp frequencies. Let $U_{\pm}\subset W_{\pm}\sim\overline{{}^{\mathrm{1c}}T^{*}\mathbb{R}^{n}}$ be open sets disjoint from fibre-infinity, and let $G_{\pm}\subset\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$ be open sets disjoint from spacetime infinity. Consider the microlocalized Lagrangians

(3.16)

$$\Lambda_{\pm}=\beta_{\mathrm{1c}-\mathrm{ps}}^{*}\big(\{(\gamma_{q_{\pm}}(s),-q_{\pm})\in\overline{\Lambda_{\pm}}\mid q_{\pm}\in U_{\pm},\ \gamma_{q_{\pm}}(s)\in G_{\pm}\}\big)$$

where $-q_{\pm}$ means changing the sign of the fiber part and we now view it as being in $\mathcal{M}$ rather than $\mathcal{M}_{0}$. That is, we take these sets to be the closure, in $\mathcal{M}$, of their interiors lifted to $\mathcal{M}$ from $\mathcal{M}_{0}$ via the blowdown map $\beta_{\mathrm{1c}-\mathrm{ps}}$. This $\Lambda_{\pm}$ depends on the choice of $U_{\pm}$ and $G_{\pm}$ but we do not indicate this in the notation, regarding these choices of open sets as fixed. As shown in [8, Proposition 4.12], this $\Lambda_{\pm}$ is admissible in the sense of Definition 3.3.

The discussion above also gives rise to the following map $\mathsf{I}_{\pm}$ that we will use later. For $\tilde{q}\in\mathrm{Char}(P)$, there is a unique $q_{\pm}\in W_{\pm}$ such that

(3.17)

$$(\tilde{q},q_{\pm})\in\overline{\Lambda^{\prime}_{\pm}}.$$

That is, $q_{\pm}$ is the endpoint of the $H_{p}^{2,0}$-flow starting from $\tilde{q}$ in the forward/backward direction. Then we view $q_{\pm}$ as a point in $\overline{{}^{\mathrm{1c}}T^{*}\mathbb{R}^{n}}$ and define

(3.18)

$$\mathsf{I}_{\pm}(\tilde{q})=q_{\pm}\in\overline{{}^{\mathrm{1c}}T^{*}\mathbb{R}^{n}}.$$

In addition, for $\tilde{q}$ at fiber infinity of $\overline{{}^{\mathrm{ps}}T^{*}\mathbb{R}^{n+1}}$, we have $\mathsf{I}_{\pm}(\tilde{q})\in\overline{{}^{\mathrm{1c}}T_{\partial\overline{\mathbb{R}^{n}}}^{*}\mathbb{R}^{n}}$.