text:
\section{Introduction}
\label{SecIntro}

Let $(M,g)$ be a Lorentzian manifold with a non-empty boundary. We consider the problem of reconstructing the topological, differentiable, and conformal structure of subsets $S\subset M^\circ$ by boundary observations of light cones emanating from points in $S$, with light rays being reflected at $\partial M$ according to Snell's law. We accomplish this under a convexity assumption on $\partial M$ and assuming that broken (reflected) null-geodesics from $S$ have no conjugate points lying on $\partial M$. The present paper is similar in spirit to the work by Kurylev, Lassas, and Uhlmann \cite{KurylevLassasUhlmannSpacetime}: they consider a related reconstruction problem using light observation sets in the \emph{interior} of globally hyperbolic spacetimes without boundary. The presence of a boundary leads to a much richer structure of the broken null-geodesic flow, and observing only \emph{at} the boundary limits the available leeway when light cones are singular (conjugate points or self-intersections) at $\partial M$.

To state a simple example to which our main result, stated below, applies, consider the manifold $M=\{(t,x)\in\mathbb{R}^{1+2}\colon |x|<1\}$, equipped with the Minkowski metric $g=-dt^2+dx^2$, and let the set $S$ of sources be an open subset $S\subseteq\{(t,x)\colon|t|<1/2-|x|\}\subset M$. The boundary light observation set from a point $q=(t_0,x_0)\in S$ within the subset $\mc U:=\{(t,x)\colon 0<t<2,\ |x|=1\}\subset\partial M$ is the intersection
\[
  \mc L^+_q \cap \mc U = \{ (t,x)\in\mc U \colon t\geq t_0,\ t-t_0=|x-x_0| \}
\]
of the future light cone from $q$ with $\mc U$. See Figure~\ref{FigIntroEx}. Let $\mc S=\{\mc L^+_q\cap\mc U\colon q\in S\}$ denote the family (as an unlabelled set) of boundary light observation sets. Then from $\mc S$, we can reconstruct $S$ as a smooth manifold, as well as the conformal class of the metric $g|_S$.

This example generalizes in a straightforward manner to higher dimensions; in $1+3$ dimensions, this would be a very simple model for wave propagation in the interior $|x|<1$, $x\in\mathbb{R}^3$, of the Earth, with observations taking place for some limited period of time on the surface of the Earth. More generally, our main theorem allows the wave speed to be inhomogeneous, anisotropic, and time-dependent.

\begin{figure}[!ht]
\centering
  \inclfig{IntroEx}
  \caption{One can recover the topological, differentiable, and conformal structure of $S$ from the collection of boundary light observation sets.}
\label{FigIntroEx}
\end{figure}

In general, the future light cone $\mc L^+_q$ from a point $q\in M$ is defined as the union of all future-directed broken null-geodesics. (See Figure~\ref{FigGeoGBroken} for an illustration, and Definition~\ref{DefGeoG} for the precise definition.) Our main theorem applies to rather general Lorentzian manifolds, and allows for the reconstruction of $S$ from boundary light observation sets involving multiple reflections. (See Remark~\ref{RmkRUsefulRefl}.) To set this up, we define the class of manifolds we will work with:

\begin{definition}
\label{DefIntroAdm}
  Let $n\geq 1$. Let $(M,g)$ be a smooth connected $(n+1)$-dimensional Lorentzian manifold with non-empty boundary; thus, $g$ has signature $(-,+,\ldots,+)$. We call $(M,g)$ \emph{admissible} if
  \begin{enumerate}
    \item there exists a proper, surjective function $t\colon M\to\mathbb{R}$ such that $dt$ is everywhere timelike;
    \item the boundary $\partial M$ is timelike, i.e.\ the induced metric $g_\partial:=g|_{\partial M}$ is Lorentzian;
    \item\label{ItIntroAdmNullConvex} $\partial M$ is \emph{null-convex}: if $\nu$ denotes the outward pointing unit normal vector field on $\partial M$, then
    \begin{equation}
    \label{EqIntroAdmNull}
      II(V,V) = g(\nabla_V\nu, V) \geq 0
    \end{equation}
    for all null vectors $V\in T_p\partial M$.
  \end{enumerate}
\end{definition}

We recall that a vector $V\in T_p M$ in a Lorentzian manifold $(M,g)$ is called \emph{timelike}, \emph{spacelike}, or \emph{lightlike} (\emph{null}) whenever $g_p(V,V)<0$, $g_p(V,V)>0$, or $g_p(V,V)=0$, respectively. An admissible manifold $(M,g)$ is time orientable, as we can declare $dt$ to be \emph{past} timelike. (We refer the reader to \cite{ONeillSemi} for further background on Lorentzian geometry.) If $n=1$, then condition~\itref{ItIntroAdmNullConvex} is vacuous.

For the purposes of this introduction, we will work with manifolds $(M,g)$ with \emph{strictly null-convex} boundaries, that is, \eqref{EqIntroAdmNull} holds with \emph{strict} inequality for $V\neq 0$. In this case, all broken null-geodesics are well-defined globally on $M$, see \S\ref{SubsecGeoTame}.

\begin{thm}
\label{ThmIntro}
  Let $(M_j,g_j)$, $j=1,2$, be two admissible Lorentzian manifolds with strictly null-convex boundaries, let $S_j\subset M_j^\circ$ be open with compact closure in $M_j$, and let $\mc U_j'\Subset\mc U_j\subset\partial M_j$ be open. Let
  \[
    \mc S_j:=\{\mc L^+_q\cap\mc U_j\colon q\in S_j\}.
  \]
  Assume that for $q_1,q_2\in\overline{S_j}$, the equality of boundary light observation sets $\mc L_{q_1}^+\cap\mc U'_j=\mc L_{q_2}^+\cap\mc U'_j$ implies $q_1=q_2$. Assume moreover that for $q\in S_j$, no point in $\mc U_j$ which lies on a future-directed broken null-geodesic starting at $q$ is conjugate to $q$.

  Suppose there exists a diffeomorphism $\Phi\colon\mc U_1\xrightarrow{\cong}\mc U_2$ which identifies the families of boundary light observation sets, that is, $\mc S_2=\{\Phi(L)\colon L\in\mc S_1\}$. Then there exists a conformal diffeomorphism $\Psi\colon(S_1,g_1|_{S_1})\xrightarrow{\cong}(S_2,g_2|_{S_2})$.

  If in addition $\Phi$ is conformal for the metrics $g_j|_{\mc U_j}$ on $\mc U_j$ and time orientation preserving, then $\Psi$ preserves the time orientation as well.
\end{thm}

Thus, if the smooth structure of the observation set $\mc U_j$ is given, then the collection of light observation sets --- carrying no structure other than that of a set! --- uniquely determines the topological, differentiable, and conformal structure of the set of sources; given a conformal structure and time orientation on $\mc U_j$, one can in addition recover the time orientation of the set of sources. See Theorem~\ref{ThmR} for the full statement which replaces the strict null-convexity condition with a certain non-degeneracy condition (called \emph{tameness} in \S\ref{SubsecGeoTame}) on broken null-geodesics.

The proof of Theorem~\ref{ThmIntro} proceeds in three steps. First, we define a topology on $\mc S_j$ by declaring collections of boundary light observation sets to be open if they intersect, resp.\ miss, a fixed open, resp.\ compact, subset of $\mc U$: this topology is shown to be equal to the subspace topology of $S_j$ via the bijection $S_j\ni q\mapsto\mc L_q^+\cap\mc U_j$; see \S\ref{SubsecRT}. Second, we show how to construct (intrinsically within $\mc S_j$ and $\mc U_j$) a large class of functions which are smooth on $S_j$: these functions $x^\mu$ assign to a point $q'$ close to a fixed point $q$ the unique parameter $x^\mu(q')$ along suitable curves $\mu\subset\mc U_j$ at which $\mu$ intersects $\mc L_{q'}^+$. (In \cite{KurylevLassasUhlmannSpacetime}, a similar construction was used \emph{globally}.) We show that \emph{all} smooth functions on $S_j$ are, locally, $\cC^\infty$ functions of these $x^\mu$ for varying $q$ and $\mu$; see \S\ref{SubsecRS}. In order to reconstruct the conformal class of $g_j$ on $S_j$, we show how to identify a large number of null-geodesics $s\mapsto q(s)$ in $S_j$ in terms of the boundary light observation sets of the points $q(s)$; see \S\ref{SubsecRC}. Since light cones are well-defined given merely the conformal class of a Lorentzian metric, one can in general not recover the metric itself. (Under additional assumptions, this may be possible, see~\cite[Corollary~1.3]{KurylevLassasUhlmannSpacetime}.) Finally, the time orientation on $S_j$ can be determined by analyzing the behavior of $\mc L_q^+\cap\mc U_j$ as $q$ moves along a timelike curve in $S_j$; see \S\ref{SubsecRO}.

It would be interesting to reconstruct suitable subsets of $(M,g)$ from \emph{active measurements}, namely from the Dirichlet-to-Neumann map of initial boundary value problems for non-linear wave equations. (In the boundary-less setting, the analogous inverse problem was first solved in the context of the quasilinear Einstein equation \cite{KurylevLassasUhlmannEinstein}, see also \cite{KurylevLassasUhlmannSpacetime}, with improvements by Lassas, Uhlmann, and Wang \cite{LassasUhlmannWangSemi,LassasUhlmannWangEinsteinMaxwell}.) The idea is to generate singular small amplitude distorted plane waves by imposing suitable singular Dirichlet data: these can be engineered so that their non-linear interaction generates point sources at points $q\in M^\circ$, allowing one to identify the boundary light observation set $\mc L_q^+\cap\partial M$ by measuring singularities of the Neumann trace; this puts one into the setting of Theorem~\ref{ThmR}. We hope to address this problem in future work. See also~\cite{BelishevKurylevBCMethod,EskinOpticalBlackHoles,LassasOksanenDNDisjoint} for results in related contexts.

For further results on the reconstruction of Lorentzian manifolds, we mention Larsson's work \cite{LarssonSkyShadow} using broken causal lens data or sky shadow data (see also the related \cite{KurylevLassasUhlmannBrokenGeodesicRigidity}), and the work by Lassas, Oksanen, and Yang \cite{LassasOksananYangTimeMeasurements} on the reconstruction of the jet of a Lorentzian metric on a timelike hypersurface from time measurements. There is a large amount of literature on inverse problems on Riemannian manifolds with boundary; we refer to \cite{PestovUhlmann2DSimple,StefanovUhlmannVasyBoundaryRidigity} and the references therein.

The plan of the paper is as follows: in \S\ref{SubsecGeoA}, we analyze the properties of admissible Lorentzian manifolds and give an equivalent formulation of the null-convexity assumption; in \S\ref{SubsecGeoG}, we define the broken null-geodesic flow and discuss its basic properties. We introduce the important notion of \emph{tameness} in \S\ref{SubsecGeoTame}; on admissible manifolds with strictly null-convex boundary, all broken null-geodesics are tame. In \S\ref{SecR} finally, we prove the main result, Theorem~\ref{ThmR}, following the steps outlined above.

\subsection*{Acknowledgments}

PH gratefully acknowledges support from the Miller Institute at the University of California, Berkeley. GU is partially supported by the NSF grant DMS-1265958, a Si-Yuan Professorship at HKUST, a Simons Fellowship, and the Academy of Finland.

\section{Geometric preliminaries}
\label{SecGeo}

\subsection{Structure of admissible manifolds}
\label{SubsecGeoA}

We begin by elucidating the smooth structure of admissible manifolds, see Definition~\ref{DefIntroAdm}:

\begin{lemma}
\label{LemmaGeoATime}
  Let $(M,g)$ be an admissible Lorentzian manifold. Then $X:=\{t=0\}$ is a compact submanifold with boundary $\partial X\subset\partial M$, and there exists a diffeomorphism $M\cong\mathbb{R}_t\times X$. Furthermore, there exists a global future timelike vector field $T\in\Vf_\bop(M)$, i.e.\ $T$ is tangent to $\partial M$, such that $T t=1$.
\end{lemma}
\begin{proof}
  Since $t$ is proper with $dt\neq 0$, the first claim is immediate. Moreover, the time orientation on $M$ induces a time orientation on $\partial M$, since the latter is assumed to be Lorentzian; with this time orientation, $dt|_{\partial M}$ is past timelike.

  Since for $O\subset M$ open the set of future timelike vector fields $V\in\Vf_\bop(O)$ with $V t=1$ is convex, it suffices to construct $T$ locally. In the interior of $M$, this is straightforward. In a neighborhood $O$ of a point $p\in\partial M$, one first constructs $T'\in\mc V(O\cap\partial M)$ with $T' t=1$; one then extends $T'$ arbitrarily to a vector field $\widetilde T\in\Vf_\bop(O)$, which thus satisfies $\widetilde T t>1/2$ in a smaller neighborhood $O'\subset O$ of $p$, thus $T=(\widetilde T t)^{-1}\widetilde T\in\Vf_\bop(O')$ is the desired vector field near $p$.

  The flow $\phi\colon\mathbb{R}\times X\in(s,x)\mapsto\exp_x(s T)\in M$ exists globally; indeed, $t(\exp_x(s T))=s$ for all $(s,x)$, since this holds for $s=0$, and the $s$-derivative of both sides is equal to $1$ by construction. The inverse of $\phi$ is given by $\phi^{-1}(p)=(a,\exp_p(-a T))$ when $p\in t^{-1}(a)$. Thus, $\phi$ establishes a diffeomorphism $\mathbb{R}\times X\cong M$.
\end{proof}

It will be useful to embed $(M,g)$ into a larger spacetime without boundary.
\begin{lemma}
\label{LemmaGeoAEmbed}
  There exists a time-oriented smooth Lorentzian manifold $(\widetilde M,\widetilde g)$ into which $M$ embeds isometrically as a submanifold with boundary.
\end{lemma}
\begin{proof}
  Let $M'$ be any open manifold into which $M$ embeds as a submanifold with boundary, e.g.\ take $M'$ to be the double of $M$. Extend $g$ to a symmetric 2-tensor $\widetilde g$ on $M'$, and extend $t$ to an arbitrary smooth function, still denoted $t$, on $M'$. Since the set of Lorentzian metrics on a fixed vector space is open, and since the condition that $dt$ is timelike (in particular $dt\neq 0$) is open, there exists an open neighborhood $\widetilde M$ of $M$ on which $\widetilde g$ is Lorentzian and $dt$ timelike; declaring $dt$ to be \emph{past} timelike endows $\widetilde M$ with a time orientation.
\end{proof}

Write $\widetilde\exp$ for the exponential map on $(\widetilde M,\widetilde g)$. Denote by $g^+$ a fixed smooth Riemannian metric on $(\widetilde M,\widetilde g)$, and write
\[
  |V|_{g^+}:=g^+(V,V)^{1/2}.
\]
(All our arguments will take place in compact subsets of $M$, hence the concrete choice of $g^+$ will be irrelevant.)

We now analyze the null-convexity condition. (We encourage the reader to keep the simpler case in mind that the boundary is \emph{strictly} null-convex.) We introduce the outward ($+$) and inward ($-$) pointing tangent bundles
\[
  T^\pm_{\partial M}M = \{ V\in T_{\partial M}M \colon \pm g(V,\nu) > 0 \},
\]
where $\nu$ is the outward pointing unit normal. Thus, $dx(\nu)<0$ for any boundary defining function $x$ (that is, $x=0$ and $dx\neq 0$ at $\partial M$, while $x>0$ in $M^\circ$), and we therefore also have
\[
  T^\pm_{\partial M}M = \{ V\in T_{\partial M}M \colon \mp dx(V) > 0 \}.
\]
Define the future/past light cones
\[
  L^\pm_p M = \{ V\in T_p M \colon V\ \text{is future ($+$), resp.\ past ($-$), lightlike} \},
\]
and the light cone
\[
  L_p M = \{ V\in T_p M \colon V\ \text{is lightlike} \} \cup \{ 0 \} \subset T_p M.
\]
As a first step, we show:

\begin{lemma}
\label{LemmaGeoABdy}
  Let $(M,g)$ be a Lorentzian manifold with null-convex timelike boundary $\partial M$ and outward pointing unit normal $\nu$. Let $p\in\partial M$. Then there exists $s_0>0$ such that for all lightlike $V\in L_p M$, $|V|_{g^+}=1$, the following holds for the null-geodesic $\gamma(s):=\widetilde\exp_p(s V)$:
  \begin{enumerate}
    \item\label{ItGeoABdyOut} If $V\in T^\pm_{\partial M}M$, then $\gamma(s)\in\widetilde M\setminus M$ for $0<\pm s\leq s_0$.
    \item\label{ItGeoABdyTgt} If $V\in T_p\partial M$ is tangent to $\partial M$, then $\gamma(s)\in\widetilde M\setminus M^\circ$ for $|s|\leq s_0$.
  \end{enumerate}
\end{lemma}
\begin{proof}
  Pick a boundary defining function $x\in\cC^\infty(\widetilde M)$, so $x^{-1}(0)=\partial M$ and $dx\neq 0$ on $\partial M$, and $x>0$ in $M^\circ$, while $x<0$ in $\widetilde M\setminus M$. Since the outward pointing unit normal to $\partial M$ is then given by $\nu=-|\nabla x|^{-1}\nabla x$, one computes
  \begin{equation}
  \label{EqGeoABdyII}
    II(V,W) = -|\nabla x|^{-1}(H x)(V,W),\ \ V,W\in T\partial M,
  \end{equation}
  where $H x=\nabla^2 x$ is the Hessian of $x$ with respect to $\widetilde g$. Therefore, the null-convexity condition is equivalent to $(H x)(V,V)\leq 0$ for all $V\in L\partial M=\bigsqcup_{p\in\partial M}L_p\partial M$.

  Denote by $y_1,\ldots,y_n$ smooth coordinates on a neighborhood $U_\partial\subset\partial M$ of $p$, with $y_j=0$ at $p$ for $1\leq j\leq n$. Using a collar neighborhood of $\partial M$, identify the set $U:=U_\partial\times(-x_0,x_0)_x$ (with $x_0>0$ small) with a neighborhood of $p$ in $\widetilde M$. We will construct a foliation of a small neighborhood of $p$ intersected with $U\cap\{x<0\}$ by \emph{strictly} null-convex hypersurfaces which will act as barriers for the geodesic $\gamma$, roughly speaking preventing it from crossing $\partial M$ into $M^\circ$ too quickly.

   To construct the foliation, let $\delta\in(0,x_0)$ and define the function
   \[
     x_\epsilon := x + \epsilon(1-\delta^{-2}Y^2)\ \ \text{for}\ \  Y:=\Bigl(\sum_{i=1}^n y_i^2\Bigr)^{1/2}<\delta^2,\ 0\leq\epsilon<\delta.
   \]
   We claim that for $\delta>0$ sufficiently small, the level sets $D_\epsilon:=x_\epsilon^{-1}(0)$ are strictly null-convex for $\epsilon>0$. To see this, note that the conormal $dx_\epsilon=dx-2\epsilon\delta^{-2}\sum_{i=1}^n y_i dy_i$ of $D_\epsilon$ is $\epsilon\delta^{-2}Y$-close (with respect to $g^+$) to $dx$; furthermore, on $D_\epsilon$, we have $x\in[-\epsilon,0]$. Given the bound we are imposing on $Y$, we conclude that null vectors $W\in L_{(x,y_1,\ldots,y_n)}D_\epsilon$ with $|W|_{g^+}=1$ are $\epsilon$-close to the boundary light cone $L_{(0,y_1,\ldots,y_n)}\partial M$. Since $\{x=0\}$ is null-convex, this implies that
   \[
     (H x)(W,W) \leq C_1\epsilon
   \]
   for some constant $C_1$. Furthermore, we have $\sum_{i=1}^n dy_i(W)^2\geq C_2>0$ for such $W$ provided $\delta>0$ is sufficiently small. Therefore,
   \begin{align*}
     (H x_\epsilon)(W,W) &= (H x)(W,W) - 2\epsilon\delta^{-2}\sum_{i=1}^n dy_i(W)^2 - 2\epsilon\delta^{-2}\sum_{i=1}^n y_i(H y_i)(W,W) \\
       &\leq C_1\epsilon - 2\epsilon\delta^{-2}C_2 - C_3\epsilon\delta^{-2}Y \\
       &\leq -C_2\epsilon\delta^{-2}
   \end{align*}
   for sufficiently small $\delta>0$, proving the strict null-convexity of $D_\epsilon$. Fixing such a $\delta>0$, define
   \[
     \mc B:=\bigl\{-\delta/2<x\leq 0,\ Y<\delta^2\bigr\}\subset U,
   \]
   and consider the function
   \[
     f := \frac{x}{1-\delta^{-2}Y^2}
   \]
   on $\mc B$, so $D_\epsilon=f^{-1}(-\epsilon)$; since $df\neq 0$ is inward pointing at $p$ (indeed, $df=dx$ there), formula~\eqref{EqGeoABdyII} shows that $(H f)(W,W)<0$ for $0\neq W\in L_p D_\epsilon$, $\epsilon>0$. See also Figure~\ref{FigGeoABdy}.

   \begin{figure}[!ht]
   \centering
     \inclfig{GeoABdy}
     \caption{Foliation of a neighborhood of $p\in\partial M$ in $\widetilde M\setminus M$ by hypersurfaces $D_\epsilon$, $\epsilon>0$, which are \emph{strictly} null-convex. These hypersurfaces are barriers for null-geodesics in the set $\mc B$: a null-geodesic $\gamma$ in $\mc B$ which emanates from a point in $\partial M$ and has outward pointing initial velocity, that is, $(x\circ\gamma)'(0)<0$, cannot cross $D_\epsilon$ in the inward direction while in $\mc B$.}
   \label{FigGeoABdy}
   \end{figure}
   Consider now $V\in L_p M\cap T^+_{\partial M}M$, $|V|_{g^+}=1$, $\gamma(s)=\widetilde\exp_p(s V)$. The point of the above construction is that the function $d(s):=f(\gamma(s))$ is negative and strictly decreasing for $s>0$ as long as $\gamma(s)\in\mc B$. Indeed, note first that we have $d(0)=0$ and $d'(0)<0$, hence $d(s),d'(s)<0$ for small $s>0$. Suppose now that $d'(s)$ vanishes for some $s>0$ with $\gamma(s)\in\mc B$, and let $s'>0$ be the first zero of $d'(s')=0$. Then, letting $\epsilon:=d(s')$, we have $\gamma'(s')\in L_{\gamma(s')}D_\epsilon$. The strict null-convexity of $D_\epsilon$ forces $d''(s')<0$, so $d'(s)$ is strictly decreasing near $s'$; since $d'(s)<0$ for $s<s'$, this contradicts the assumption that $d'(s')=0$.

   Therefore, we have
   \[
     s_0:=\inf_{\genfrac{}{}{0pt}{}{V\in L_p M\cap T^+_{\partial M}M}{|V|_{g^+}=1}}\sup\{s>0\colon \widetilde\exp_p((0,s]V)\subset\mc B^\circ\}>0.
   \]
   (In fact, our arguments show $s_0\gtrsim\delta^2$.) The conclusion of part~\itref{ItGeoABdyOut} then holds for this value of $s_0$.

   Part~\itref{ItGeoABdyTgt} follows from part~\itref{ItGeoABdyOut} by a simple limiting argument: let $V_\epsilon:=V+\epsilon\nu$, $|\epsilon|<1$, which is outward pointing for $\epsilon>0$ and inward pointing for $\epsilon<0$. By part~\itref{ItGeoABdyOut}, there exists $s_0>0$ such that $\widetilde\exp_p(s V_\epsilon)\in\widetilde M\setminus M$ for $0<(\operatorname{sgn}\epsilon)s\leq s_0$. Letting $\epsilon\to 0$, this implies $\gamma(s)\in\overline{\widetilde M\setminus M}=\widetilde M\setminus M^\circ$ for $0\leq|s|\leq s_0$, as claimed.
\end{proof}

We can now give a useful equivalent formulation of the null-convexity condition.

\begin{prop}
\label{PropGeoANull}
  Let $(M,g)$ be a Lorentzian manifold with timelike boundary $\partial M$ and outward pointing unit normal $\nu$. Then the following are equivalent:
  \begin{enumerate}
    \item\label{ItGeoANullConvex} $\partial M$ is null-convex, i.e.\ the inequality~\eqref{EqIntroAdmNull} holds.
    \item\label{ItGeoANullOutward} If $\gamma\colon(-\epsilon,0]\to M$ is a null-geodesic segment with $\gamma(0)\in\partial M$ and $\gamma(s)\in M^\circ$ for $s\leq 0$, then $\gamma'(0)\in T^+_{\partial M}M$. Likewise, if $\gamma\colon[0,\epsilon)\to M$ is a null-geodesic segment with $\gamma(0)\in\partial M$ and $\gamma(s)\in M^\circ$ for $s>0$, then $\gamma'(0)\in T^-_{\partial M}M$.
  \end{enumerate}
\end{prop}
\begin{proof}
  \underline{\itref{ItGeoANullConvex} $\Longrightarrow$ \itref{ItGeoANullOutward}}: for a null-geodesic $\gamma\colon(-\epsilon,0]\to M$ as in \itref{ItGeoANullOutward}, the conclusion $\gamma'(0)\in T^+_{\partial M}M\cup T\partial M\setminus\{0\}$ is clear. But by Lemma~\ref{LemmaGeoABdy}, which uses condition~\itref{ItGeoANullConvex}, $\gamma'(0)\in T\partial M$ would imply that $\gamma(s)\in\widetilde M\setminus M^\circ$ for small $s$. Hence $\gamma'(0)\not\in T\partial M$.

  \underline{\itref{ItGeoANullOutward} $\Longrightarrow$ \itref{ItGeoANullConvex}}: suppose that condition~\itref{ItGeoANullConvex} is \emph{violated}, hence there exists $V\in L_p\partial M$, $p\in\partial M$, with $II(V,V)<0$, in particular $V\neq 0$. Define $\gamma(s)=\exp_p(s V)$ for $s\in[0,\epsilon)$, $\epsilon>0$ small, and let $f=x\circ\gamma\colon[0,\epsilon)\to\mathbb{R}$, with $x$ a boundary defining function as in the proof of Lemma~\ref{LemmaGeoABdy}. Then
  \[
    f(0)=0,\ \ f'(0)=0,\ \ f''(0)=(H x)(V,V)>0.
  \]
  Therefore, $\gamma(s)\in M^\circ$ for $s\in(0,s_1)$ for sufficiently small $s_1\in(0,\epsilon)$. Since $\gamma'(0)=V\in T\partial M$, this contradicts condition~\itref{ItGeoANullOutward}.
\end{proof}

We end this section with a geometric lemma linking boundary light observation sets with spacetime light cones on an infinitesimal level. We denote by $\rho(V)$, $V\in T_p M$, $p\in\partial M$, the reflection of $V$ across $\partial M$, that is,
\begin{equation}
\label{EqGeoARefl}
  \rho(V) := V - 2 g(V,\nu)\nu,
\end{equation}
with $\nu$ the outward pointing unit normal. One easily checks $\rho\colon L M\to L M$. Moreover, if $T\in T\partial M$, then $g(\rho(V),T)=g(V,T)$; this in particular applies to future timelike $T$, hence $\rho\colon L^\pm M\to L^\pm M$ preserves the time orientation of lightlike vectors.

\begin{lemma}
\label{LemmaGeoANull}
  Suppose $(M,g)$ is a time-oriented manifold with timelike boundary $\partial M$. Let $p\in\partial M$. Then there exists an isomorphism $\phi$ between the space $\ms S$ of spacelike hypersurfaces $S\subset T_p\partial M$ and the space $\ms V$ of rays $\mathbb{R}_+ V\subset T_p M$ along future-directed outward pointing null vectors, given by mapping $S\in\ms S$ to the unique future-directed outward pointing null ray $\phi(S)$ contained in $S^\perp$. The inverse map is given by $\ms V\ni\mathbb{R}_+ V\mapsto T_p\partial M\cap V^\perp\in\ms S$.

  Moreover, there exists an isomorphism between $\ms S$ and the space $\ms N$ of null hypersurfaces $N\subset T_p M$ which contain a future-directed outward pointing null vector, given by $\ms S\ni S\mapsto S\oplus\operatorname{span} \phi(S)\in\ms N$.
\end{lemma}

See Figure~\ref{FigGeoANull}.

\begin{figure}[!ht]
\centering
  \inclfig{GeoANull}
  \caption{Illustration of Lemma~\ref{LemmaGeoANull}; everything takes place in $T_p M$. $N$ is a null hypersurface containing an outward pointing null vector $V$, while $S$ is a spacelike hypersurface in $T_p\partial M$. We have $\mathbb{R}_+V=\phi(S)$, $N^\perp=\mathbb{R} V$, and $S=N\cap T_p\partial M$.}
\label{FigGeoANull}
\end{figure}

\begin{proof}[Proof of Lemma~\ref{LemmaGeoANull}]
  Given a spacelike hypersurface $S\subset T_p\partial M$, the orthocomplement $S^\perp$ is a time-oriented 2-dimensional vector space with signature $(1,1)$, hence there exists a non-zero null vector $W\in S^\perp$; the four distinct rays of null vectors contained in $S^\perp$ are then the positive scalar multiples of $W$, $-W$, $\rho(W)$, $-\rho(W)$. Since multiplication by $-1$ exchanges future- and past-directed null as well as outward and inward pointing vectors, and since application of $\rho$ exchanges outward and inward pointing vectors but preserves the time orientation, exactly one of these four vectors, which we call $V$, is future-directed and outward pointing; and $\phi^+(S)=\mathbb{R}_+ V$.

  On the other hand, if $0\neq V\in T_p M$ is null (thus $V^\perp/\mathbb{R} V$ is spacelike) and outward pointing, in particular $V\not\in T_p\partial M$, then the composition $V^\perp\cap T_p\partial M\hookrightarrow V^\perp\to V^\perp/\mathbb{R} V$ is an isometric isomorphism, hence $S:=V^\perp\cap T_p\partial M$ is a spacelike hypersurface. This establishes the isomorphism $\ms S\cong\ms V^+$ (as smooth manifolds).

  For the last claim, we note that $\ms N^+\ni N\mapsto N^\perp\cap L^+_p M\in\ms V^+$ maps a null hypersurface $N$ into the unique ray along a future-directed outward pointing null generator of $N$. The inverse of this map is given by $\ms V^+\ni\mathbb{R}_+ V\mapsto V^\perp=V\oplus(V^\perp\cap T_p\partial M)\in\ms N^+$. Composition of these maps with $\phi^+$ gives the desired isomorphism $\ms S\xrightarrow{\cong}\ms N^+$. The inverse of this isomorphism is given by $N\mapsto N\cap T_p\partial M$.
\end{proof}

\subsection{Examples of admissible manifolds}
\label{SubsecGEx}

Small perturbations of admissible Lo\-rentz\-ian manifolds with \emph{strictly} null-convex boundaries are admissible:

\begin{lemma}
\label{LemmaGExOp}
  Suppose $(M,g)$ is admissible and strictly null-convex, with an embedding $(M,g)\hookrightarrow(\widetilde M,\widetilde g)$ as in Lemma~\ref{LemmaGeoAEmbed}. Let $K\Subset\widetilde M$, and define $\mc C^k$ spaces using the Riemannian metric $g^+$ on $\widetilde M$.
  \begin{enumerate}
    \item Let $x\in\cC^\infty(\widetilde M)$ denote a defining function of $\partial M$. If $x'\in\cC^\infty(\widetilde M)$ is equal to $x$ outside of $K$ and sufficiently close in $\mc C^2$ to $x$ in $K$, then $M':=\{x'\geq 0\}$ is admissible and strictly null-convex.
    \item If $g'$ is a smooth Lorentzian metric on $M$, equal to $g$ on $M\setminus K$ and sufficiently close in $\mc C^1(K)$ to $g$, then $(M,g')$ is admissible.
  \end{enumerate}
\end{lemma}
\begin{proof}
  This follows from the observation that the assumption of strict null-convexity involves up to first derivatives of the metric and up to second derivatives of the boundary defining function, see \eqref{EqGeoABdyII}.
\end{proof}

If more is known about the global structure of $(M,g)$, one can allow non-compact perturbations as well. For example, the cylinder
\begin{equation}
\label{EqGExM0}
  M_0 := \{ (t,x) \in \mathbb{R}\times\mathbb{R}^n \colon |x|<R \}, \ \ R>0,
\end{equation}
with the Minkowski metric $g=-dt^2+dx^2$, is admissible with strictly null-convex boundary; indeed,
\[
  II(V,V) = R^{-1}dx(V)^2
\]
is strictly positive for non-zero null vectors $V\in T\partial M$. If now $f\colon\mathbb{R}^{1+n}\to\mathbb{R}$ has small $\mc C^2$ norm, then
\begin{equation}
\label{EqGExMf}
  M_f := \{ (t,r\omega) \in \mathbb{R}\times\mathbb{R}^n \colon r<(1+f(t,\omega))R \}
\end{equation}
is admissible, with strictly null-convex boundary; see Figure~\ref{FigGExAdm}.

\begin{figure}[!ht]
\centering
  \inclfig{GExAdm}
  \caption{An admissible manifold $(M_f,g)\hookrightarrow(\mathbb{R}^{1+n}_{t,x},-dt^2+dx^2)$, with $f$ having small $\mc C^2$ norm.}
\label{FigGExAdm}
\end{figure}

Another interesting class of examples, which includes the cylinder~\eqref{EqGExM0}, is obtained as follows: let $(X,h)$ be a compact Riemannian manifold with convex boundary, so $II(V,V)=h(\nabla_V \nu,V)\geq 0$ for all $V\in T\partial X$, where $\nu$ is the outward pointing unit normal. (We thus allow for the possibility that parts of the boundary are totally geodesic.) Then the product manifold $M:=\mathbb{R}_t\times X$, $g=-dt^2+h$, is admissible, with $\partial M$ strictly null-convex if and only if $\partial X$ is strictly convex. See Figure~\ref{FigGExProd}; another example is (partially) shown in Figure~\ref{FigRSTrans}.

\begin{figure}[!ht]
\centering
  \inclfig{GExProd}
  \caption{An admissible manifold, obtained as the product of the real line with a Riemannian manifold $X$ with convex boundary. We can allow parts of $\partial X$ to be flat.}
\label{FigGExProd}
\end{figure}

\subsection{Broken null-geodesics}
\label{SubsecGeoG}

Throughout this section, $(M,g)$ will be a fixed admissible Lorentzian manifold. Motivated by the fact that singularities of solutions of wave equations on $(M,g)$ propagate along null-geodesics in $M^\circ$ and undergo reflection according to Snell's law at the boundary $\partial M$, see Taylor~\cite{TaylorReflection}, we rigorously define and study such \emph{broken null-geodesics} in this section. Define the open submanifold
\begin{equation}
\label{EqGeoGRedLM}
  L^{\mathrm{b}} M := L M \setminus ( T^+_{\partial M}M \cup T\partial M )
\end{equation}
of $L M$, so $(p,V)\inL^{\mathrm{b}} M$ if and only if $V\in L_p M$ and $V\in T^-_{\partial M}M$ when $p\in\partial M$. We then introduce:

\begin{definition}
\label{DefGeoG}
  Let $(p,V)\inL^{\mathrm{b}} M$. We call a piecewise smooth curve $\gamma\colon I\to M$, with $0\in I\subset\mathbb{R}$ an open connected interval, $\gamma(0)=p$, $\gamma'(0)=V$, a \emph{broken null-geodesic}, if
  \begin{enumerate}
  \item for all open intervals $J\subset I$ with $\gamma(J)\cap\partial M=\emptyset$, $\gamma|_J$ is an affinely parameterized null-geodesic in $(M^\circ,g)$;
  \item if $s\in I$, $\gamma(s)\in\partial M$, then for small $\epsilon>0$, $\gamma|_{(s-\epsilon,s]}$ and $\gamma|_{[s,s+\epsilon)}$ are null-geodesics with $\gamma(s\pm(0,\epsilon))\subset M^\circ$, and $\gamma'(s+0)=\rho(\gamma'(s-0))$, where $\rho$ is the reflection map~\eqref{EqGeoARefl}.
  \end{enumerate}
\end{definition}

Thus, broken null-geodesics are null-geodesics which undergo reflection at the boundary $\partial M$ preserving their velocity tangent to $\partial M$; see Figure~\ref{FigGeoGBroken}.

\begin{figure}[!ht]
\centering
  \inclfig{GeoGBroken}
  \caption{A broken null-geodesic undergoing a reflection at $\gamma(s)\in\partial M$. Here, $V=\gamma'(s)$.}
\label{FigGeoGBroken}
\end{figure}

A broken null-geodesic with $(\gamma(0),\gamma'(0))=(p,V)$ as in this definition always exists on sufficiently small intervals $I=(-\epsilon,\epsilon)$, $\epsilon>0$: when $p\in M^\circ$, $\gamma|_I$ is an interior null-geodesic, while for $p\in\partial M$, one takes $\gamma(s)=\exp_p(s V)$ for $s\geq 0$ and $\gamma(s)=\exp_p(s\rho(V))$ for $s\leq 0$. Also note that if $\gamma_j\colon I_j\to M$, $j\in J$, are broken null-geodesics which all have the same initial condition, then the prescription $\gamma|_{I_j}=\gamma_j$ defines a broken null-geodesic $\gamma\colon \bigcup_{j\in J}I_j\to M$. Thus, for $(p,V)$ as in the above definition, there always exists an \emph{inextendible broken null-geodesic} with initial position $p$ and speed $V$.

\begin{definition}
\label{DefGeoGExp}
  For $(p,V)\inL^{\mathrm{b}} M$, let $\gamma\colon I\to M$ denote the unique inextendible broken null-geodesic with $(\gamma(0),\gamma'(0))=(p,V)$. Suppose $1\in I$. We then define the \emph{broken exponential map} by $\exp^{\mathrm{b}}_p(V):=\gamma(1)$. Denote the domain of definition of $\exp^{\mathrm{b}}$ by $\mc D\subsetL^{\mathrm{b}} M$.
\end{definition}

Thus, for $(p,V)$ with $\gamma((0,1))\subset M^\circ$, we simply have $\exp^{\mathrm{b}}_p(V)=\exp_p(V)$. We proceed to analyze the properties of inextendible broken null-geodesics. For convenience, we make our choice of the Riemannian metric $g^+$ on $\widetilde M$ more specific by demanding
\begin{equation}
\label{EqGeoGReflIso}
  g^+(\rho(V),\rho(V))=g^+(V,V),\ \ V\in T_p M,\ p\in\partial M.
\end{equation}
This is easily arranged by taking any Riemannian metric $g^+_0$ on $T_{\partial M}M$, then letting $g^+_1(V,W):=g^+_0(V,W)+g^+_0(\rho(V),\rho(W))$ for $V,W\in T_{\partial M}M$, and finally taking $g^+$ to be a Riemannian metric on $\widetilde M$ extending $g^+_1$ smoothly to the rest of $\widetilde M$. The consequence of \eqref{EqGeoGReflIso} is that the $g^+$-length of the tangent vector of a broken null-geodesic $\gamma$ is continuous when $\gamma$ hits the boundary.

\begin{prop}
\label{PropGeoGMax}
  Let $\gamma\colon I\to M$ be a broken null-geodesic with $\gamma'(0)\in L^+_p M$, and let $I_+:=\sup I\in\mathbb{R}\cup\{\infty\}$. Then $\gamma$ is future inextendible\footnote{By this we mean that the parameter interval for which the maximal broken null-geodesic with the same initial data as $\gamma$ is defined has supremum equal to $\sup I$.} if and only if one of the following happens:
  \begin{enumerate}
  \item\label{ItGeoGMaxT} $t(\gamma(s))\to\infty$ as $s\to I_+$.
  \item\label{ItGeoGMaxL} $I_+<\infty$, $t(\gamma(s))\to t_\infty<\infty$ as $s\to I_+$, and $I_+\in\overline{\gamma^{-1}(\partial M)}$.
  \end{enumerate}
  There exists an analogous characterization of past inextendibility.
\end{prop}

In other words, a broken null-geodesic is future-inextendible if and only if it leaves every region $\{t\leq t_0\}$, $t_0<\infty$ (this may happen even in the case $I_+<\infty$, e.g.\ for $M=\{(\tilde t,x)\colon|\tilde t|<\pi/2,|x|\leq 1\}$ with the metric $g=-dt^2+dx^2$, $t=\tan\tilde t$), or it undergoes infinitely many reflections as $s\to I_+<\infty$; similarly for past inextendibility. We remark that the latter scenario can indeed occur in certain cases when $\partial M$ is flat to infinite order; see \cite[\S6]{TaylorGrazing}.

\begin{proof}[Proof of Proposition~\ref{PropGeoGMax}]
  We note that $dt(\gamma'(s))>0$ for all $s\in I$ since $dt$ is past timelike and $\gamma'(s)$ is future causal; hence $t\circ\gamma$ is strictly increasing.

  If $I_+=\infty$, then $\gamma$ is clearly future inextendible. Suppose $I_+<\infty$ and $t(\gamma(s))\to\infty$ as $s\to I_+$. If there were an extension $\gamma_1\colon I_1\to M$, $I_+\in I_1$, of $\gamma$, then
  \[
    \infty>t(\gamma_1(I_+))=\lim_{\epsilon\to 0+} t(\gamma_1(I_+-\epsilon)) = \lim_{\epsilon\to 0+} t(\gamma(I_+-\epsilon)) = \infty,
  \]
  a contradiction. We next claim that
  \begin{equation}
  \label{EqGeoGMaxEscape}
    I_+=\infty\ \Longrightarrow\ t(\gamma(s))\to\infty,\ s\to\infty.
  \end{equation}
  Taking this for granted, the assumption $t(\gamma(s))\to t_\infty<\infty$ implies $I_+<\infty$; moreover, $\gamma(s)$ stays in a fixed compact set as $s\to I_+$ since $t$ is proper. If there exists a broken null-geodesic $\gamma_1\colon I_1\to M$ extending $\gamma$, $I_+\in I_1$, then since $\gamma_1^{-1}(\partial M)\subset I_1$ is discrete by definition, we infer that $I_+$ is not a limit point of $\gamma_1^{-1}(\partial M)$, hence not of $\gamma^{-1}(\partial M)$. Conversely, if $I_+\not\in\overline{\gamma^{-1}(\partial M)}$, then $s_0:=\max\{s\in I\colon \gamma(s)\in\partial M\}<I_+$. Let $(p_0,V_0):=(\gamma(s_0),\gamma'(s_0+0))$, then
  \[
    \gamma_1(s) := \widetilde\exp_{p_0}\bigl((s-s_0)V_0\bigr)
  \]
  satisfies $\gamma_1(s)=\gamma(s)\in M^\circ$ for $s\in(s_0,\infty)\cap I$. If $\gamma_1(I_+)\not\in\partial M$, then $\gamma_1|_{I\cup(s_0,I_++\epsilon)}$, defined for small $\epsilon>0$, is an extension of $\gamma$. Otherwise, $\gamma_1$ intersects $\partial M$ at $s=I_+$, and it necessarily does so \emph{transversally} according to Proposition~\ref{PropGeoANull}; hence we can continue $\gamma_1(s)$ past $s=s_1$ as a broken null-geodesic by defining
  \[
    \gamma_1(s_1+s'):=\widetilde\exp_{\gamma_1(s_1)}\bigl(s'\rho(\gamma_1'(s_1-0))\bigr),
  \]
  $s'>0$ small. This construction shows that $\gamma$ is future extendible past $I_+$.

  It remains to prove~\eqref{EqGeoGMaxEscape}. Assume to the contrary that
  \begin{equation}
  \label{EqGeoGTInfty}
    t(\gamma(s))\to t_\infty<\infty\ \ \text{as}\ \ s\to\infty.
  \end{equation}
  Since $K:=t^{-1}([t(\gamma(0)),t_\infty])$ is compact, there exists $c>0$ such that
  \begin{equation}
  \label{EqGeoGMaxBound}
    dt(V) \geq c|V|_{g^+}, \ \ V\in L^+_K M.
  \end{equation}
  Define $\ell(s):=|\gamma'(s)|_{g^+}$. Since the difference of connections $D=\nabla^{g^+}-\nabla^g$ induces a bilinear map $T_p M\otimes T_p M\to T_p M$, $D_p(X,Y)=\nabla^{g^+}_X Y-\nabla^g_X Y$, we can write for $s$ with $\gamma(s)\in M^\circ$:
  \[
    \frac{1}{2}\frac{d}{ds}\bigl(\ell(s)^2\bigr) = g^+|_{\gamma(s)}\bigl(\nabla^{g^+}_{\gamma'(s)}\gamma'(s),\gamma'(s)\bigr) = g^+|_{\gamma(s)}\bigl(D_{\gamma(s)}(\gamma'(s),\gamma'(s)),\gamma'(s)\bigr)
  \]
  On the other hand, if $\gamma(s)\in\partial M$, then $\ell(s+0)=\ell(s-0)$ in view of \eqref{EqGeoGReflIso}. Since $\gamma(s)\in K$ remains in a compact set, this implies
  \[
    \frac{1}{2}\frac{d}{ds}\bigl(\ell(s)^2\bigr) \geq -C\ell(s)^3
  \]
  where $C>0$ is a uniform constant only depending on $K$. Rewriting this differential inequality as $(1/\ell)'\leq C$, we obtain
  \begin{equation}
  \label{EqGeoGMaxLengthBound}
    \ell(s) \geq \frac{1}{\frac{1}{\ell(0)}+C s}.
  \end{equation}
  Therefore, the bound~\eqref{EqGeoGMaxBound} implies
  \[
    t(\gamma(s)) \geq t(\gamma(0)) + \int_0^s c\ell(u)\,du \geq t(\gamma(0)) + \frac{c}{C}\log\bigl(1+C s\ell(0)\bigr),
  \]
  which exceeds $t_\infty$ for sufficiently large $s$, contradicting \eqref{EqGeoGTInfty}. The proof is complete.
\end{proof}

We next study the regularity properties of the broken exponential map. For $(p,V)\in\mc D$, the domain of definition of $\exp^{\mathrm{b}}$, consider the maximal broken null-geodesic $\gamma(s)=\exp^{\mathrm{b}}_p(s V)$, $s\in I$, let
\[
  \mathsf{J}(p,V):=\#\{s>0\colon\gamma(s)\in\partial M\}\in\mathbb{N}_0\cup\{\infty\}
\]
denote the number of reflections at $\partial M$, and enumerate the affine parameters for which $\gamma$ intersects the boundary:
\[
  \gamma^{-1}(\partial M)\cap(0,\infty) =: \{\mathsf{s}_j(p,V) \colon j=1,\ldots,\mathsf{J}(p,V)\},\ \ 0<\mathsf{s}_j(p,V)<\mathsf{s}_{j+1}(p,V).
\]
Further, let
\[
  \mathsf{p}_j(p,V):=\gamma(\mathsf{s}_j(p,V)),\ \ \mathsf{V}_j(p,V):=\gamma'(\mathsf{s}_j(p,V)+0) \in L^{\mathrm{b}}_{\mathsf{p}_j(p,V)}M
\]
denote the position and the velocity of the broken null-geodesic leaving the boundary at a reflection point. For $k\in\mathbb{N}_0$, define
\[
  \mc D_k^\circ := \bigl\{ (p,V)\in\mc D \colon \#\{s\in(0,1)\colon \gamma(s)\in\partial M\}=k,\ \exp^{\mathrm{b}}_p(V)\not\in\partial M\bigr\},
\]
i.e.\ $k$ is the number of reflections of the broken null-geodesic segment $\gamma((0,1))$. Let $\mc D_k$ denote the closure of $\mc D_k^\circ$ in $\mc D$. See Figure~\ref{FigGeoGExpRefl}.

\begin{figure}[!ht]
\centering
  \inclfig{GeoGExpRefl}
  \caption{A broken null-geodesic with initial data $(p,V)\in\mc D_2^\circ$. Here $\mathsf{p}_j=\mathsf{p}_j(p,V)=\exp^{\mathrm{b}}_p(\mathsf{s}_j(p,V)V)$ and $\mathsf{V}_j=\mathsf{V}_j(p,V)$.}
\label{FigGeoGExpRefl}
\end{figure}

\begin{prop}
\label{PropGeoGExpCont}
  The broken null-geodesic flow on an admissible Lorentzian manifold $(M,g)$ has the following properties:
  \begin{enumerate}
  \item\label{ItGeoGExpContOpenk} For every $k\in\mathbb{N}_0$, the set $\mc D_k^\circ\subsetL^{\mathrm{b}} M$ is open. The functions $\mathsf{s}_j(p,V)$ as well as the points $(\mathsf{p}_j(p,V),\mathsf{V}_j(p,V))\inL^{\mathrm{b}} M$ depend smoothly on $(p,V)\in\mc D_k^\circ$, and so does $\exp^{\mathrm{b}}_p(V)\in M^\circ$.
  \item\label{ItGeoGExpCl} We have the decomposition $\partial\mc D_k=B_{k,-}\sqcup B_{k,+}$ into a disjoint union of the closed sets
    \begin{align*}
      B_{k,-} &= \{ (p,V)\in\mc D \colon \mathsf{J}(p,V)\geq k,\ \mathsf{s}_k(p,V)=1 \}, \\
      B_{k,+} &= \{ (p,V)\in\mc D \colon \mathsf{J}(p,V)\geq k+1,\ \mathsf{s}_{k+1}(p,V)=1 \}.
    \end{align*}
    See Figure~\ref{FigGExpCl}. Furthermore, $\mathsf{s}_j$, $\mathsf{p}_j$, and $\mathsf{V}_j$ for $1\leq j\leq k$ as well as $\exp^{\mathrm{b}}$ extend from $\mc D_k^\circ$ to smooth functions on $\mc D_k$.
  \item\label{ItGeoGExpContOpen} $\mc D$ is open in $L^{\mathrm{b}} M$; more precisely, $\mc D_k^\circ\cup\bigcup_{j<k}\mc D_j$ is open for all $k\in\mathbb{N}_0$. The map $\exp^{\mathrm{b}}$ is continuous on $\mc D$.
  \end{enumerate}
\end{prop}

\begin{figure}[!ht]
\centering
  \inclfig{GeoGExpCl}
  \caption{Illustration of part~\itref{ItGeoGExpCl} of Proposition~\ref{PropGeoGExpCont}. \textit{Left:} initial data $(p,V)\in B_{1,+}$. \textit{Right:} initial data $(p,V)\in B_{1,-}$, with $(\mathsf{p}_1,\mathsf{V}_1)$ defined by smooth extension from $\mc D_1^\circ$.}
\label{FigGExpCl}
\end{figure}

\begin{proof}[Proof of Proposition~\ref{PropGeoGExpCont}]
  \underline{\itref{ItGeoGExpContOpenk}}: for $k=0$ and $(p,V)\in\mc D_0^\circ$, we have $\exp^{\mathrm{b}}_p(s V)=\widetilde\exp_p(s V)$ for $0\leq s\leq 1$, and $\widetilde\exp_p(V)\in M^\circ$; hence the smooth dependence of $\exp^{\mathrm{b}}_p(V)$ on $(p,V)$ follows from that of the standard exponential map $\widetilde\exp_p(V)$.

  Consider now $k\geq 1$. Denote the map acting by dilation by $c\in\mathbb{R}$ in the fibers by $R_c\colon L M\to L M$, $(p,V)\mapsto (p, c V)$. Let $(p,V)\in\mc D_k^\circ$, and let $\bar\mathsf{s}_j=\mathsf{s}_j(p,V)$, $\bar\mathsf{p}_j=\mathsf{p}_j(p,V)$, and $\bar\mathsf{V}_j=\mathsf{V}_j(p,V)$ for $j=1,\ldots,k$. We start by defining neighborhoods of $(\bar\mathsf{p}_j,\bar\mathsf{V}_j)$ of initial conditions of null-geodesics for which the next intersection with $\partial M$ is controlled. Thus, for $j=0$, let
  \[
    Z_0 := (\widetilde\exp\circ R_{\bar\mathsf{s}_1})^{-1}(\partial M).
  \]
  Note that $(p,V)\in Z_0$. Moreover, in a neighborhood of $(p,V)$, $Z_0$ is a smooth codimension 1 submanifold of $L^{\mathrm{b}} M$ which is transversal to $\mathbb{R}_+ V\subsetL^{\mathrm{b}}_p M$; this follows from the implicit function theorem applied to the map $x\circ\widetilde\exp\circ R_{\bar s_1}$ where $x\in\cC^\infty(M)$ is a boundary defining function: this map has a non-zero differential at $(p,V)$ due to part~\itref{ItGeoANullOutward} of Proposition~\ref{PropGeoANull}. See Figure~\ref{FigGeoGExpContZ0}.

  \begin{figure}[!ht]
  \centering
    \inclfig{GeoGExpContZ0}
    \caption{Illustration of the proof of part~\itref{ItGeoGExpContOpenk} of Proposition~\ref{PropGeoGExpCont}. Shown are the preimage $Z_0$ of $\partial M$ under $\widetilde\exp\circ R_{\bar\mathsf{s}_1}$ within $L_p M$ (only the future half of which is drawn) near $V$, and the intersection of the neighborhood $\mc U_0$ of $(p,V)$ with $L_p M$.}
  \label{FigGeoGExpContZ0}
  \end{figure}

  For a small neighborhood $\mc U_0\subsetL^{\mathrm{b}} M$ of $(p,V)$ such that $\mc U_0\subset\bigcup_{c\in(1-\epsilon,1+\epsilon)}R_c Z_0$, with $\epsilon>0$ small, and such that $Z_0\cap\overline{\mc U_0}$ is a smooth connected submanifold transversal to all dilation (in the fiber) orbits intersecting $\mc U_0$, define the function $d_0\in\cC^\infty(\mc U_0)$ by
  \begin{equation}
  \label{EqGeoGExpContScaling}
    d_0(q,W) := \bar\mathsf{s}_1 c, \ \ \ R_c(q,W)\in Z_0,\ c\in(\tfrac{1}{1+\epsilon},\tfrac{1}{1-\epsilon}),
  \end{equation}
  so $d_0(p,V)=\bar\mathsf{s}_1$, and $\widetilde\exp_q(R_{d_0(q,W)}W)\in\partial M$ for $(q,W)\in\mc U_0$. Similarly, but now working within $\partial M$, we define for $1\leq j\leq k-1$
  \[
    Z_j^\partial := \bigl(\widetilde\exp\circ R_{\bar\mathsf{s}_{j+1}-\bar\mathsf{s}_j}|_{L^{\mathrm{b}}_{\partial M}M}\bigr)^{-1}(\partial M) \subset L^{\mathrm{b}}_{\partial M}M,
  \]
  so $(\bar\mathsf{p}_j,\bar\mathsf{V}_j)\in Z_j^\partial$, and near this point, $Z_j^\partial$ is a smooth codimension 1 submanifold of $L^{\mathrm{b}}_{\partial M}M$ transversal to $\mathbb{R}_+\bar\mathsf{V}_j\subsetL^{\mathrm{b}}_{\bar\mathsf{p}_j}M$. For a small neighborhood $\mc U_j^\partial\subsetL^{\mathrm{b}}_{\partial M}M$ of $(\bar\mathsf{p}_j,\bar\mathsf{V}_j)$ such that $\mc U_j^\partial\subset\bigcup_{c\in(1-\epsilon,1+\epsilon)}R_c Z_j^\partial$, with $Z_j^\partial\cap\overline{\mc U_j^\partial}$ smooth and connected, define $d_j\in\cC^\infty(\mc U_j^\partial)$ by
  \[
    d_j(q,W) = (\bar\mathsf{s}_{j+1}-\bar\mathsf{s}_j)c,\ \ \ R_c(q,W)\in Z_j^\partial,
  \]
  so $d_j(\bar\mathsf{p}_j,\bar\mathsf{V}_j)=\bar\mathsf{s}_{j+1}-\bar\mathsf{s}_j$, and $\widetilde\exp_q(R_{d_j(q,W)}W)\in\partial M$ for $(q,W)\in\mc U_j^\partial$. Lastly, let $\mc U_k^\partial\subsetL^{\mathrm{b}}_{\partial M}M$ denote a small neighborhood of $(\bar\mathsf{p}_k,\bar\mathsf{V}_k)$ such that $\mc U_k^\partial\subset(\widetilde\exp\circ R_{1-\bar\mathsf{s}_k})^{-1}(M^\circ)$; in particular $(\bar\mathsf{p}_k,\bar\mathsf{V}_k)\in\mc U_k^\partial$.

  We now construct a neighborhood of $(p,V)$ for which the $j$-th reflection point and velocity lie in $\mc U_j^\partial$. Thus, encoding point and velocity of the extended manifold by the map
  \[
    \widetilde\exp'(q,W) := \bigl(\widetilde\exp_q(W), \tfrac{d}{ds}\widetilde\exp_q(s W)\bigr|_{s=1}\bigr)
  \]
  we inductively define $\mc V_k^\partial:=\mc U_k^\partial$ and, for $1\leq j\leq k-1$,
  \[
    V_j^\partial := \mc U_j^\partial \cap \bigl(\rho\circ R_{d_j}^{-1}\circ\widetilde\exp'\circ R_{d_j}\bigr)^{-1}(V_{j+1}^\partial),
  \]
  where we used the reflection map $\rho$ defined in equation~\eqref{EqGeoARefl}, and finally
  \[
    \mc V_0 := \mc U_0 \cap \bigl(\rho\circ R_{d_0}^{-1}\circ\widetilde\exp'\circ R_{d_0}\bigr)^{-1}(V_1^\partial).
  \]
  Then $\mc V_0$ is the desired neighborhood of $(p,V)$. Indeed, if $(q,W)\in\mc V_0$, we inductively define $(q_0,W_0)=(q,W)$, and for $j=0,\ldots,k-1$
  \[
    (q_{j+1},W_{j+1}) := \rho\Bigl(R_{d_j(q_j,W_j)}^{-1}\bigl(\widetilde\exp'\bigl(R_{d_j(q_j,W_j)}(q_j,W_j)\bigr)\bigr)\Bigr) \in \mc V_{j+1}^\partial.
  \]
  Then we have $\exp^{\mathrm{b}}_q(\mathsf{s}_j(q,W) W)=q_j$ for $1\leq j\leq k$, where
  \[
    \mathsf{s}_j(q,W) = \sum_{i<j} d_i(q_i,W_i).
  \]
  Therefore $\exp^{\mathrm{b}}_q(W)=\widetilde\exp_{q_k}\bigl((1-s_k(q,W))W_k\bigr)$, and by construction, we also have $(\mathsf{p}_j(q,W),\mathsf{V}_j(q,W))=(q_j,W_j)$, with smooth dependence on $(q,W)$.

  \underline{\itref{ItGeoGExpCl}}: Suppose $\mc D\supset\partial\mc D_k\ni(p,V)=\lim_{i\to\infty}(p_i,V_i)$ with $(p_i,V_i)\in\mc D_k^\circ$, and denote $\bar\mathsf{J}:=\mathsf{J}(p,V)$, $\bar\mathsf{s}_j:=\mathsf{s}_j(p,V)$. The above arguments imply $\bar\mathsf{J}\geq k-1$, and
  \[
    \bar\mathsf{s}_j = \lim_{i\to\infty} \mathsf{s}_j(p_i,V_i),\ \ j\leq\min(k,\bar\mathsf{J}),
  \]
  likewise $(\bar\mathsf{p}_j,\bar\mathsf{V}_j):=(\mathsf{p}_j(p,V),\mathsf{V}_j(p,V))=\lim_{i\to\infty} (\mathsf{p}_j(p_i,V_i),\mathsf{V}_j(p_i,V_i))$ for these $j$. Let $\gamma(s)=\exp^{\mathrm{b}}_p(s V)$ and $\gamma_i(s)=\exp^{\mathrm{b}}_{p_i}(s V_i)$.

  Suppose first that $\bar\mathsf{J}\geq k$, then $\bar\mathsf{s}_j\in(0,1]$ for $j\leq k$. If $\bar\mathsf{s}_k=1$, then $\gamma(s)$ undergoes $(k-1)$ reflections and ends at $\gamma(1)\in\partial M$, and $(p,V)\in B_{k,-}$. If $\mathsf{s}_k<1$, then both the case $\mathsf{J}=k$ and the case $\mathsf{J}\geq k+1$, $\mathsf{s}_{k+1}>1$, would imply $(p,V)\in\mc D_k^\circ$. Hence, we must have $\mathsf{J}\geq k+1$, and $\mathsf{s}_{k+1}\leq 1$. If $\mathsf{s}_{k+1}<1$, then the arguments for \itref{ItGeoGExpContOpenk} would imply that $\gamma_i((0,1))$ intersects $\partial M$ at least $k+1$ times for large $i$. Thus necessarily $\mathsf{s}_{k+1}=1$, and $(p,V)\in B_{k,+}$.

  In order to exclude the case that $\bar\mathsf{J}=k-1$ ($k\geq 1$), note that $(p,V)\in\mc D$ implies that we can define $\gamma(s)$ as a broken null-geodesic for $s\in[0,1+\epsilon]$ for some $\epsilon>0$. We claim that $\gamma$ necessarily has a $k$-th intersection point with $\partial M$, contradicting $\bar\mathsf{J}<k$. If $k=1$, this is straightforward, as $\gamma|_{[0,1+\epsilon]}$ not intersecting $\partial M$ would imply (by continuity of $\widetilde\exp$) the same statement for $\gamma_i|_{[0,1+\epsilon]}$, contradicting $V_i\in\mc D_1^\circ$. For $k\geq 2$, we note that
  \[
    \mathsf{p}_k(p_i,V_i) = \widetilde\exp_{\mathsf{p}_{k-1}(p_i,V_i)}\bigl((\mathsf{s}_k(p_i,V_i)-\mathsf{s}_{k-1}(p_i,V_i))\mathsf{V}_{k-1}(p_i,V_i)\bigr);
  \]
  passing to a subsequence, we may assume that $\mathsf{s}_k(p_i,V_i)\to\bar\mathsf{s}\in[\bar\mathsf{s}_{k-1},1]$. Since $\lim_{i\to\infty}\mathsf{V}_{k-1}(p_i,V_i)=\bar\mathsf{V}_{k-1}\in T^-_{\partial M}M$ is \emph{strictly} inward pointing (as follows from the definition of $\exp^{\mathrm{b}}$ and $\mc D$), there exists $c>0$ such that for all $i$, $\mathsf{s}_k(p_i,V_i)-\mathsf{s}_{k-1}(p_i,V_i)\geq c$; therefore $\bar\mathsf{s}>\bar\mathsf{s}_{k-1}$, and we obtain
  \[
    \widetilde\exp_{\bar\mathsf{p}_{k-1}(p,V)}\bigl((\mathsf{s}-\mathsf{s}_{k-1}(p,V))\mathsf{V}_{k-1}(p,V)\bigr) \in \partial M,
  \]
  so indeed $\bar\mathsf{J}\geq k$ (and $\mathsf{s}_k(p,V)=\bar\mathsf{s}$).

  This proves the inclusion $\partial\mc D_k\subset B_{k,-}\sqcup B_{k,+}$. (The disjointness of the two sets on the right is evident.) For the reverse inclusion, we note that $(p,V)\in B_{k,-}$ is the limit as $\epsilon\to 0$ of $(p,(1+\epsilon)V)\in\mc D_k^\circ$ (this uses that $\mathsf{s}_j(p,c V)=c^{-1}\mathsf{s}_j(p,V)$ for $c>0$), while $(p,V)\in B_{k,+}$ is the limit of $(p,(1-\epsilon)V)\in\mc D_k^\circ$. The smooth extendibility of $\mathsf{s}_j$ etc.\ follows easily from the construction used in the proof of part~\itref{ItGeoGExpContOpenk}.

  \underline{\itref{ItGeoGExpContOpen}}: Note that $B_{k,+}=B_{k+1,-}$, so the family $\exp^{\mathrm{b}}|_{\mc D_k}$, $k\in\mathbb{N}_0$, of smooth maps does glue to a continuous function on $\bigcup\mc D_k$; furthermore, $\mc D=\bigcup_{k\in\mathbb{N}_0}\mc D_k$ by definition of broken null-geodesics. In view of \itref{ItGeoGExpContOpenk}, and noting that $B_{0,-}=\emptyset$, it remains to show that every $(p,V)\in B_{k-1,+}$, $k\in\mathbb{N}$, has an open neighborhood in $L^{\mathrm{b}} M$ which is contained in $\mc D_{k-1}\cup\mc D_k^\circ$; but this follows again from the proof of part~\itref{ItGeoGExpContOpenk}.
\end{proof}

\subsection{Tame broken null-geodesics}
\label{SubsecGeoTame}

We define the class of tame null-geodesics for which the possibility~\itref{ItGeoGMaxL} in Proposition~\ref{PropGeoGMax} does not occur for a given range of values of $t$:

\begin{definition}
\label{DefGeoTame}
  We call an inextendible broken null-geodesic $\gamma\colon I\to M$ \emph{tame for $-\infty\leq a<t<b\leq\infty$} if for all $a<a',b'<b$, we have $t(\gamma(I))\cap(a,a'),t(\gamma(I))\cap(b',b)\neq\emptyset$. If $\gamma$ is tame for $-\infty<t<\infty$, we simply say that $\gamma$ is \emph{tame}.
\end{definition}

By Proposition~\ref{PropGeoGMax} and its proof, this can be rephrased as follows: an inextendible broken null-geodesic $\gamma$ is tame for $a<t<b$ if and only if the only possible accumulation points of $\gamma^{-1}(\partial M)\cap(a,b)\subset\overline\mathbb{R}=\mathbb{R}\cup\{\pm\infty\}$ are $a$ and $b$; that is, $\gamma$ only undergoes a finite number of reflections whenever $t\circ\gamma$ stays in a fixed compact subset of $(a,b)$. Tame geodesics are precise those for which $t(\gamma(I))=\mathbb{R}$.

From the point of view of solving boundary value problems for wave equations, we have precise control over the singularities of geometric optics solutions along tame broken null-geodesics \cite{TaylorReflection}. (There are much more general results about the propagation of singularities for boundary value problems, see for example \cite{TaylorGrazing,MelroseSjostrandSingBVPI,MelroseSjostrandSingBVPII,MelroseTaylorParametrix}, which would become relevant if one dropped the null-convexity assumption on $\partial M$.)

To illustrate Definition~\ref{DefGeoTame} and to provide a natural class of examples, we show:

\begin{prop}
\label{PropGeoTameStr}
  If $(M,g)$ is admissible with \emph{strictly} null-convex boundary, then all inextendible broken null-geodesics $\gamma\colon I\to M$ are tame.
\end{prop}

This is a generalization of \cite[Lemma~6.1]{TaylorGrazing}.

\begin{proof}[Proof of Proposition~\ref{PropGeoTameStr}]
  Assume the conclusion is false, then we must have $I_+:=\sup I<\infty$, and $\gamma^{-1}(\partial M)\cap[0,I_+)=\{b_j\colon j\in\mathbb{N}_0\}\subset I$, with
  \begin{equation}
  \label{EqGeoTameStrTimes}
    0\leq b_{j-1}<b_j\to I_+, \ \ t(\gamma(b_j))\to t_\infty<\infty.
  \end{equation}
  Denote $p_j:=\gamma(b_j)\in\partial M$. By the proof of Proposition~\ref{PropGeoGMax}, in particular the estimate~\eqref{EqGeoGMaxLengthBound}, there exists a constant $C_+>1$ such that $C_+^{-1}\leq |\gamma'(s)|_{g^+}\leq C_+$ for all $s\in[0,I_+)$; thus, $\gamma(s)$ is uniformly continuous, which implies that the limit $\lim_{s\to I_+}\gamma(s)=:p_\infty\in\partial M$ exists.

  Letting
  \[
    V_j:=\gamma'(b_j+0)\inL^{\mathrm{b,+}}_{p_j}M,
  \]
  we claim that $V_j$ converges to some $0\neq V_\infty\in L_{p_\infty}\partial M$, i.e.\ $V_\infty$ is tangent to the boundary; note that $C_+^{-1}\leq|V_j|_{g^+}\leq C_+$ for all $j$, proving that any subsequential limit of the $V_j$ must be a \emph{non-zero} element of $L\partial M$. To prove the convergence, denote by $\nu$ the outward unit normal to $\partial M$, and assume to the contrary that there is a subsequence $V_{j_k}$ such that $|g(V_{j_k},\nu)|\geq C_\partial>0$ for some fixed constant $C_\partial$. Using a finite number of local coordinate charts covering the compact set $K:=\partial M\cap t^{-1}([t(\gamma(0)),t_\infty])$, one easily sees that
  \begin{align*}
    C_R := \inf_{q\in K} \sup \{ t \colon \widetilde\exp_q((0,t)W) &\subset M^\circ\ \text{for all}\ W\in (T^-_{\partial M}M)_q, \\
      & C_+^{-1}\leq |W|_{g^+}\leq C_+,\ |g(W,\nu)|\geq C_\partial \}
  \end{align*}
  is \emph{positive}, as follows from the fact that in a local coordinate chart and for such $W$, we have $\widetilde\exp_q(s W)=q+s W+\mc O(s^2)$, which does not return to $\partial M$ for a uniform amount of time (depending on $C_\partial,C_+,K$, and the $\mc C^1(K)$ norm of the metric $g$). But then $b_{j_k+1}-b_{j_k}\geq C_R$, contradicting \eqref{EqGeoTameStrTimes}. A similar argument shows more generally that
  \begin{equation}
  \label{EqGeoTameStrGammaPrime}
    \lim_{s\to I_+}\gamma'(s)=V_\infty \in L_{p_\infty}\partial M.
  \end{equation}
  By affinely reparameterizing $\gamma$, we may assume $|V_\infty|_{g^+}=1$.

  Fix a boundary defining function $x$, and let
  \[
    f:= x\circ\gamma \geq 0,
  \]
  then $f$ is continuous on the closed interval $[0,I_+]$, with $f(b_j)=0$ for all $j$; therefore $\lim_{s\to I_+}f(s)=f(I_+)=0$. Let further
  \[
    \theta_j:=f'(b_j+0)=dx(V_j)>0,
  \]
  then $\lim_{j\to\infty}\theta_j=0$. We aim to prove estimates on the `chord lengths' $b_{j+1}-b_j$ and the `reflection angles' $\theta_j$ as $j\to\infty$ which will contradict the convergence \eqref{EqGeoTameStrTimes}; our arguments will slightly more generally prove that reflection points cannot accumulate near a strict null-convex boundary point.

  The strict null-convexity of $\partial M$ at $p_\infty$ implies, by continuity, that $(H x)(V,V)\leq-k<0$ for some constant $k>0$ whenever $V\in L_p\partial M$, $|V|_{g^+}=1$, $p$ near $p_\infty$. For large $j$ then, by~\eqref{EqGeoTameStrGammaPrime}, we have
  \[
    f''(s) = (H x)(\gamma'(s),\gamma'(s)) \leq -k/2,\ \ s\in(b_j,b_{j+1}),
  \]
  which gives an estimate for how close $\gamma$ stays to $\partial M$:
  \begin{equation}
  \label{EqGeoStrTameClose}
    f(s)\leq\theta_j(s-b_j)-k(s-b_j)^2/4\leq \theta_j^2/k.
  \end{equation}
  Furthermore, $f(b_{j+1})=0$ implies the estimate
  \begin{equation}
  \label{EqGeoStrTameClose2}
    b_{j+1}-b_j\leq 4\theta_j/k.
  \end{equation}
  Consider now a reflection point $p_j$, $j$ large, then $V_j$ is $\theta_j$-close to a null vector $V'_j\in L_{p_j}\partial M$. Let $k_j:=-(H x)(V'_j,V'_j)>0$, then the smoothness of $H x$ and the estimates~\eqref{EqGeoStrTameClose}--\eqref{EqGeoStrTameClose2} give
  \[
    f''(s) = -k_j + \mc O(\theta_j + f(s) + |s-b_j|)=-k_j+\mc O(\theta_j),\ \ s\in(b_{j-1},b_j).
  \]
  We also record that $k_j\geq k/2$ for large $j$. Therefore, for such $s$, we have
  \begin{align*}
    f'(s) &= \theta_j - (s-b_j)(k_j+\mc O(\theta_j)), \\
    f(s) &= (s-b_j)\theta_j - (s-b_j)^2(k_j+\mc O(\theta_j))/2.
  \end{align*}
  Hence, $f(b_{j+1})=0$ implies
  \begin{equation}
  \label{EqGeoTameStrBDiff}
    b_{j+1}-b_j = \frac{2\theta_j}{k_j+\mc O(\theta_j)} \geq C_B\theta_j,
  \end{equation}
  and thus
  \[
    \theta_{j+1} = -f'(b_{j+1}) = -\theta_j + \frac{2\theta_j(k_j+\mc O(\theta_j)}{k_j+\mc O(\theta_j)} \geq \theta_j - C\theta_j^2,
  \]
  with $C$ and $C_B>0$ constants independent of $j$.

  Fix $j_0$ such that $\theta_j<\frac{1}{2 C}$ for $j\geq j_0$. Since $x\mapsto x-C x^2$ is increasing for $x<\frac{1}{2 C}$, we conclude that $\theta_j\geq a_j/C$, where $a_{j_0}=C \theta_{j_0}\in(0,1/2)$ and
  \[
    a_{j+1} = a_j - a_j^2.
  \]
  Since $a_j\geq\frac{C_A}{j-(j_0-1)}$ by Lemma~\ref{LemmaGeoTameSeq} below, the estimate~\eqref{EqGeoTameStrBDiff} implies that $b_j\geq C_0 + C_1\log j$ for some constants $C_0$ and $C_1>0$, contradicting \eqref{EqGeoTameStrTimes}. The proof is complete.
\end{proof}

\begin{lemma}
\label{LemmaGeoTameSeq}
  If $a_1\in(0,1/2)$ and $a_{j+1}=a_j-a_j^2$, then $a_j \geq C/j$ for some $C>0$.
\end{lemma}
\begin{proof}
  Clearly, $a_j>0$ for all $j$. Write $a_j = b_j / j$, then $b_j>0$, and
  \[
    b_{j+1} = b_j\left(1+\frac{1}{j}\Bigl(1-\frac{j+1}{j}b_j\Bigr)\right)
  \]
  If $b_j\leq\frac{j}{j+1}$ (this holds for  $j=1$), this gives $b_j\leq b_{j+1}\leq 1$. If on the other hand $\frac{j}{j+1}\leq b_j\leq 1$, then $\frac{j-1}{j}\leq b_{j+1}\leq 1$. Thus, $b_j\geq C$ for some $C>0$, as claimed.
\end{proof}

\section{Reconstruction from boundary light observation sets}
\label{SecR}

In this section, we prove (a generalization of) Theorem~\ref{ThmIntro}, showing how one can reconstruct the topological, smooth, and conformal structure of suitable precompact subsets $S\Subset M$ from the observation of light cones on (subsets of) the null-convex boundary $\partial M$, following the arguments outlined in the introduction.

There are substantial differences compared to the arguments in \cite{KurylevLassasUhlmannSpacetime} due to the presence of a boundary which we will explain in more detail below: the boundary allows for the reconstruction of $S$ using (multiply) reflected rays; it necessitates certain restrictions on $S$ due to possible strong refocusing after reflection; and the codimension $1$ nature of $\partial M$ causes complications when there are null conjugate points on $\partial M$ --- we circumvent the latter by \emph{assuming} that there are no such conjugate points in the set $\mc U\subset\partial M$ where we observe the future light cones from points in $S$.

Let $(M,g)$ denote an admissible manifold.

\begin{definition}
\label{DefRConj}
  Let $(q,V)\inL^{\mathrm{b}} M$, see~\eqref{EqGeoGRedLM}, and suppose $p:=\exp^{\mathrm{b}}_q(V)\in\partial M$. Then we say that $(q,V)$ and $p$ are \emph{not conjugate} if $\exp^{\mathrm{b}}_q|_{L_q M}$ has injective differential at $V$, where we define the differential as the limit $D_{(1-\epsilon)V}\exp^{\mathrm{b}}_q$ as $\epsilon\to 0+$.
\end{definition}

If $\exp^{\mathrm{b}}_q(s V)\not\in\partial M$ for $0<s<1$, this can be phrased equivalently as the condition that the exponential map $\widetilde\exp_q|_{L_p M}$ has injective differential at $V$.

The existence of the limit follows from part~\itref{ItGeoGExpCl} of Proposition~\ref{PropGeoGExpCont}, since $(q,V)\in B_k^+$ for some $k\in\mathbb{N}_0$. Since broken null-geodesics are transversal to $\partial M$, we can rephrase the definition as follows: denote $Z=(\exp^{\mathrm{b}}_q)^{-1}(\partial M)\subsetL^{\mathrm{b}}_q M$, which is a smooth codimension $1$ submanifold near $V$ (see the proof of Proposition~\ref{PropGeoGExpCont}). Then $(q,V)$ and $p$ are not conjugate if and only if the map $Z\ni W\mapsto\exp^{\mathrm{b}}_q(W)\in\partial M$ has injective differential at $V$; that is, the boundary point near $p$ depends non-degenerately on the initial velocity. See also Figure~\ref{FigGeoGExpContZ0} for a closely related setting.

The implicit function theorem immediately gives:
\begin{lemma}
\label{LemmaRConj}
  If $(q,V)$ and $p$ are not conjugate, then, in the above notation, there exists a neighborhood $U\subset Z$ of $V$ such that $\exp^{\mathrm{b}}_q|_U\colon U\to\exp^{\mathrm{b}}_q(U)$ is a diffeomorphism onto its image, which is thus a 1-codimensional smooth submanifold of $\partial M$.
\end{lemma}

Recalling~\eqref{EqGeoGRedLM}, denote by
\[
  L^{\mathrm{b,+}} M := L^{\mathrm{b}} M \cap L^+ M = L^+M \setminus (T^+_{\partial M}M\cup T\partial M)
\]
the set of future-directed light-like vectors which are inward pointing at the boundary. We then define by
\begin{equation}
\label{EqRFutLight}
  \mc L_q^+ := \overline{\exp^{\mathrm{b}}_q(L^{\mathrm{b,+}}_q M\cap\mc D)}, \ \ q\in M,
\end{equation}
the future light cone from $q$. Thus, if $q\in M^\circ$, we simply have $\mc L_q^+=\{q\}\cup\exp^{\mathrm{b}}_q(L_q^+ M)$. (The set on the right hand side is already closed since $\exp^{\mathrm{b}}_q$ is proper; this uses that there exists a global timelike function on $M$.)

\begin{thm}
\label{ThmR}
  Let $(M_j,g_j)$, $j=1,2$, be admissible Lorentzian manifolds, let $S_j\subset M_j^\circ$ be open with compact closure in $M_j$, and let $\mc U_j'\Subset\mc U_j\subset\partial M_j$ be open. Denote the collection of light observation sets by
  \begin{equation}
  \label{EqRObsSets}
    \mc S_j := \{ \mc L_q^+\cap\mc U_j \colon q\in S_j \}.
  \end{equation}
  Assume:
  \begin{enumerate}
    \item \label{ItRDistinct} for any two points $q_1\neq q_2\in\overline{S_j}$, we have $\mc L_{q_1}^+\cap\mc U'_j\neq\mc L_{q_2}^+\cap\mc U'_j$.
    \item \label{ItRTame} all inextendible broken null-geodesics passing through a point in $\overline{S_j}$ are tame, see Definition~\ref{DefGeoTame};
    \item \label{ItRConj} for all $q\in S_j$ and $V\inL^{\mathrm{b,+}}_q M_j$ such that $p=\exp^{\mathrm{b}}_q(V)\in\mc U_j$, $(q,V)$ and $p$ are not conjugate.
  \end{enumerate}
  Suppose there exists a diffeomorphism $\Phi\colon\mc U_1\xrightarrow{\cong}\mc U_2$ such that
  \[
    \mc S_2=\{\Phi(L)\colon L\in\mc S_1\}.
  \]
  Then there exists a conformal diffeomorphism $\Psi\colon(S_1,g_1|_{S_1})\xrightarrow{\cong}(S_2,g_2|_{S_2})$.

  If in addition $\Phi$ is conformal, i.e.\ $\Phi^*(g_2|_{\mc U_2})=f g_1|_{\mc U_1}$ for some function $f\neq 0$, and preserves the time orientation, then $\Psi$ preserves the time orientation as well.
\end{thm}

In fact, we will show that the map $\Psi\colon S_1\to S_2$ given by the composition of $S_1\ni q\mapsto\mc L_q^+\cap\mc U_1\in\mc S_1$, $\Phi$, and the inverse of $S_2\ni q\mapsto\mc L_q^+\cap\mc U_2\in\mc S_2$ is a conformal (and time orientation preserving) diffeomorphism.

\begin{rmk}
\label{RmkRTame}
  For a general admissible manifold $(M,g)$, the constructions below allow for the reconstruction of $S$ from light observation sets if the closure $\bar S$ of the set $S$ of light sources as well as the subset of the boundary on which we observe are contained in a fixed slab $M':=t^{-1}((I_-,I_+))$, $-\infty\leq I_-<I_+\leq\infty$, with the property that all inextendible broken null-geodesics passing through a point in $\bar S$ are tame for $I_-<t<I_+$. One can then define a new time function $t'$, proper as a map $M'\to\mathbb{R}$, such that $t'\to\pm\infty$ as $t\to I_\pm$. Replacing $M$ by $M'$, condition~\itref{ItRTame} is satisfied.
\end{rmk}

\begin{rmk}
\label{RmkRUsefulRefl}
  Theorem~\ref{ThmR} allows for the reconstruction of subsets $S\subset M^\circ$ even in certain situations in which the first intersection point of future null-geodesics from sources in $S$ with $\partial M$ is not contained in $\mc U$; that is, the theorem crucially uses possibly multiply reflected broken null-geodesics. As an example, in Figure~\ref{FigRSTrans}, one can take $S\subset M^\circ$ and $\mc U\subset\partial M$ to be small neighborhoods of $q$ and $p$, respectively; if $\mc U$ is sufficiently small, then the shown \emph{once broken} null-geodesics are the \emph{only} broken null-geodesics starting at $q$ and intersecting $\mc U$.
\end{rmk}

Assumption~\itref{ItRDistinct} is very natural: we illustrate this with two examples.

\begin{example}
\label{ExRDistinct}
  Consider the cylinder $M_0=\{(t,r\omega)\colon r<1\}\subset\mathbb{R}\times\mathbb{R}^n$, $n\geq 1$, of radius $1$, see also equation~\eqref{EqGExM0}. Let $S_1=\{(t,r\omega)\colon |t|<1/2-r,\ r<1/2\}$ and $\mc U=(0,2)\times\mathbb{S}^{n-1}$, $\mc U'=[1/2,3/2]\times\mathbb{S}^{n-1}$. Then Theorem~\ref{ThmR} applies: the topological, differentiable, and conformal structure of $S_1$ can be recovered from the light observation sets $\mc L_q^+\cap\mc U$, $q\in S_1$.

  Denoting by $T_a\colon(t,x)\mapsto(t+a,x)$ the time translation operator, let $S_2:=S_1\cup T_{-2}(S_1)$. Using the notation \eqref{EqRObsSets}, we then have $\mc S_1=\mc S_2$, hence observers in $\mc U$ cannot distinguish $S_1$ and $S_2$, even though the sets $S_1$ and $S_2$ are not homeomorphic ($S_1$ is connected, $S_2$ is not); assumption~\itref{ItRDistinct} is violated. See the left panel of Figure~\ref{FigRDistinct}.
\end{example}

\begin{figure}[!ht]
\centering
  \inclfig{RDistinct}
  \caption{Illustration of examples which violate assumption~\itref{ItRDistinct} of Theorem~\ref{ThmR}. \textit{Left:} example~\ref{ExRDistinct}. The sets $S_1$ and $S_2$ of sources have the same light observation sets in $\mc U$, for instance $\mc L_{q_1}^+\cap\mc U=\mc L_{q_2}^+\cap\mc U$ for $q_1=(0,0)$ and $q_2=(-2,0)$. \textit{Right:} example~\ref{ExRDistinctCl}. All light observation sets from points in $S$ are distinct in $\mc U$, but as $q\to(-2,0)$, the observation set $\mc L_q^+\cap\mc U$ converges to $\mc L_{(0,0)}^+\cap\mc U$. (The light cone based at $q$ does \emph{not} refocus near $(0,0)$ in three and more spacetime dimensions due to its distorted form, contrary to the appearance in this 2-dimensional picture.)}
\label{FigRDistinct}
\end{figure}

\begin{example}
\label{ExRDistinctCl}
  Consider $M_f\subset\mathbb{R}\times\mathbb{R}^n$, $n\geq 2$, defined in equation~\eqref{EqGExMf} with $R=1$, for the function $f(t,\omega)=\chi_0(t)\chi(\omega)$, where $\chi_0(t)\equiv 0$ for $t\geq 0$ and $\chi_0(t)=\delta e^{1/t}$ for $t\leq 0$, and where $\chi\in\cC^\infty(\mathbb{S}^{n-1})$ is identically $1$ in the neighborhood $|\omega-\omega_0|<1/2$ of some fixed $\omega_0\in\mathbb{S}^{n-1}\subset\mathbb{R}^n$, and $\chi(\omega)=0$ for $|\omega-\omega_0|>1$. See the right panel of Figure~\ref{FigRDistinct}. For $\delta>0$ sufficiently small, $M_f$ has a strictly null-convex boundary by Lemma~\ref{LemmaGExOp}. Let
  \begin{align*}
    S &= S' \cup S'', \\
      &\qquad S' = \{ (t,r\omega)\in\mathbb{R}^{1+n} \colon |t|<1/2-r,\ r<1/2 \}, \\
      &\qquad S''= \{ (t,r\omega)\in\mathbb{R}^{1+n} \colon |t+9/4|<1/4-r,\ r<1/4 \}.
  \end{align*}
  We use the observation set $\mc U=(0,2)\times\mathbb{S}^{n-1}\subset\partial M_f$ and $\mc U'=[1/2,3/2]\times\mathbb{S}^{n-1}$. Theorem~\ref{ThmR} applies to the set $S'$, and in fact yields a conformal diffeomorphism (which in this case is just the identity map on $\mathbb{R}^{1+n}$) between $(S',-dt^2+dx^2)$ and $(S_1,-dt^2+dx^2)$ from Example~\ref{ExRDistinct}. Theorem~\ref{ThmR} can also be shown (by a perturbative argument off the case $\delta=0$) to apply to $S''$ and $\mc U$ for small $\delta>0$. If one attempts to recover $S$, all light observation sets $\mc L_q^+\cap\mc U$, $q\in S$, are distinct. However, we have
  \[
    \lim_{t\to -2} \mc L_{(t,0)}\cap\mc U \to \mc L_{(0,0)}\cap\mc U
  \]
  as smooth submanifolds of $\mc U$. That is, separated points can have very similar light observation sets. This motivates the stronger hypothesis that light observation sets from points in the \emph{closure} $\bar S$ are distinct.
\end{example}

Fix an admissible Lorentzian manifold $(M,g)$, and sets $S\subset M$ and $\mc U'\Subset\mc U\subset\partial M$ satisfying the assumptions of Theorem~\ref{ThmR}, and denote $\bar\mc S=\{\mc L_q^+\cap\mc U\colon q\in\bar S\}$. By assumption~\itref{ItRDistinct}, the map
\begin{equation}
\label{EqRLMap}
  L\colon \bar S\ni q\mapsto\mc L_q^+\cap\mc U\in\bar\mc S
\end{equation}
is bijective, as is its restriction to $S$ as a map $S\to\mc S$. There exists a unique topological, smooth, and conformal structure on $\mc S$, defined by pushing these structures forward from $S\subset M$ to $\mc S$ using $L$, which makes this map a conformal diffeomorphism. In order to prove Theorem~\ref{ThmR}, we need to show that we can uniquely recover these structures merely from the knowledge of the collection $\mc S$ of subsets of the manifold $\mc U$ and the conformal class of $g|_{\mc U}$. \emph{From now on, we identify the set $S$ of sources and the set $\mc S$ of light observation sets using the map~\eqref{EqRLMap}, and use the two interchangeably.}

The proof of Theorem~\ref{ThmR} will occupy the remainder of this section: in \S\ref{SubsecRT}, we show how to recover the topology of $\mc S$, in \S\ref{SubsecRS} we recover the smooth structure, in \S\ref{SubsecRC} the conformal structure, and finally in \S\ref{SubsecRO} the time orientation of $S$.

\subsection{Topology}
\label{SubsecRT}

We define a topology $\mc T$ on $\mc S$ by using the collection of sets of the form
\begin{align*}
  U_O &:= \{ L\in\mc S \colon L\cap O\neq\emptyset \}, \ \ O\subset\mc U\ \text{open},\\
  U^K &:= \{ L\in\mc S \colon L\cap K=\emptyset \}, \ \ K\subset\mc U\ \text{compact},
\end{align*}
as a subbasis. Note that the definition of $\mc T$ only involves the \emph{set} $\mc S$ and the a priori known topology of $\mc U$.

\begin{prop}
\label{PropRT}
  The topology $\mc T$ is equal to the subspace topology $\mc T_M$ of $S\subset M$.
\end{prop}
\begin{proof}
  \underline{$\mc T\subset\mc T_M$}: We show that sets of the form $U_O$ and $U^K$ and open in $S\subset M$. If $O=\emptyset$, then $U_O=\emptyset$ is open. If on the other hand $O\neq\emptyset$, let $L=\mc L_q^+\cap\mc U\in\mc S$, $q\in S$, and suppose $V\inL^{\mathrm{b,+}}_q M$ is such that $\exp^{\mathrm{b}}_q(V)\in O$; in the notation of Proposition~\ref{PropGeoGExpCont}, we have $V\in B_{k,-}$ for some $k\in\mathbb{N}$, i.e.\ $p$ is the $k$-th intersection of the broken null-geodesic with initial data $(q,V)$ with the boundary $\partial M$, and $\mathsf{p}_k(V)=p$. By part~\itref{ItGeoGExpCl} of that proposition, $\mathsf{p}_k(W)$ depends continuously on $q'\in M$, $W\inL^{\mathrm{b}}_{q'} M$, hence $\mathsf{p}_k(W)\in O$ when $(q',W)$ is close to $(q,V)$. This shows that $\mc L_{q'}^+\cap O\neq\emptyset$, as desired.

  For $K\subset\mc U$ compact, we claim that $S\setminus U^K$ is closed in the subspace topology of $M$: if $q_j\in S$, $\lim_{j\to\infty}q_j=:\bar q\in S$, and $V_j\inL^{\mathrm{b,+}}_{q_j}M$, $p_j:=\exp^{\mathrm{b}}_{q_j}(V_j)\in K$, then, passing to a subsequence if necessary, we may assume that $p_j\to\bar p\in K$ as $j\to\infty$. Moreover, it follows from the proof of Proposition~\ref{PropGeoGMax}, see in particular the estimates~\eqref{EqGeoGMaxBound} and \eqref{EqGeoGMaxLengthBound}, that $|V_j|_{g^+}$ remains in a compact subset of $(0,\infty)$, hence we may assume that $(q_j,V_j)\to(\bar q,\bar V)\inL^{\mathrm{b,+}} M$. But by Proposition~\ref{PropGeoGExpCont}, we then have $\bar p=\exp^{\mathrm{b}}_{\bar q}(\bar V)\in K$, hence $\bar q\not\in U^K$, as claimed.

  \underline{$\mc T_M\subset\mc T$}: we need to prove that for any $\mc T_M$-open set $U\subset S$, every $q\in U$ has a $\mc T$-open neighborhood which is contained in $U$. To see this, denote $L:=\mc L_q^+\cap\overline{\mc U'}$, and fix a compact set $K'$ with $\mc U'\Subset K'\subset\mc U$; for any $\epsilon>0$, let $K_\epsilon:=K'\setminus\bigl\{p\in\partial M\colon d_{g^+}(p,L)<\epsilon\bigr\}$, where $d_{g^+}(p,L)=\inf_{p'\in L}d_{g^+}(p,p')$ is defined using the Riemannian distance function of $g^+$. Further, pick a countable dense subset $\{p_i\colon i\in\mathbb{N}\}\subset\mc L_q^+\cap\overline{\mc U'}$, and let $O_{i,\epsilon}=\{p\in\mc U\colon d_{g^+}(p,p_i)<\epsilon\}$. By the compactness of $L$, for each $\epsilon>0$, there exists a finite number $N(\epsilon)$ such that $L\subset\bigcup_{i=1}^{N(\epsilon)} O_{i,\epsilon}$. Consider now the nested sequence of $\mc T$-open neighborhoods
  \[
    U_j := U^{K_{1/j}} \cap \bigcap_{i=1}^{N(1/j)} U_{O_{i,1/j}}
  \]
  of $\mc L^+_q\cap\mc U$.

  Suppose that $U_j\not\subset U$ for all $j$, then we can pick a sequence $q_j\in U_j\setminus U\subset S$, and we may assume without loss that $q_j\to\bar q\in\bar S$. It then follows that $\mc L_{\bar q}^+\cap\overline{\mc U'}$ is equal to the set of limit points of the sequence of sets $\mc L_{q_j}^+\cap\overline{\mc U'}$; for $\bar q\not\in\partial M$, this is a consequence of Proposition~\ref{PropGeoGExpCont}, while for $\bar q\in\partial M$, recalling the definition~\eqref{EqRFutLight}, this follows from a simple approximation argument. By definition of the sets $U^{K_{1/j}}$, we infer that $\mc L_{\bar q}^+\cap\overline{\mc U'}\subset L$. If this were a strict inclusion (of \emph{closed} sets), we could find $i_0\in\mathbb{N}$ with $p_{i_0}\in L\setminus\mc L_{\bar q}^+$ and $j_0\in\mathbb{N}$ such that $O_{i_0,1/j}\subset L\setminus\mc L_{\bar q}^+$ for all $j\geq j_0$. However, by definition of $U_{O_{i_0,1/j}}$, there exists, for all $j\geq j_0$, a point $x_j\in\mc L_{q_j}^+\cap O_{i_0,1/j}$; hence $p_{i_0}=\lim_{j\to\infty}x_j$ is a limit point of the sets $\mc L_{q_j}^+$, hence contained in $\mc L^+_{\bar q}$, which is a contradiction. Therefore, $\mc L_{\bar q}^+\cap\overline{\mc U'}=L=\mc L_q^+\cap\overline{\mc U'}$. By assumption~\itref{ItRDistinct} of Theorem~\ref{ThmR}, this implies $\bar S\setminus U\ni\bar q=q\in U$. This contradiction shows that $U_j\subset U$ for sufficiently large $j$, and the proof is complete.
\end{proof}

\begin{example}
\label{ExREarliest}
  A key construction in \cite{KurylevLassasUhlmannSpacetime} is the earliest observation time along timelike curves in the observation region. We give an example to indicate why, without modifications as in \S\ref{SubsecRS} below, this is not as useful in the present setting. Consider the cylinder $M_0\subset\mathbb{R}^{1+n}$, $n\geq 1$, with radius $1$, see equation~\eqref{EqGExM0}, and consider the set
  \[
    S = \{ (t,r\omega)\in\mathbb{R}^{1+n} \colon |t+1|<1/2-r,\ r<1/2 \}.
  \]
  We observe in the set $\mc U=(0,3)\times\mathbb{S}^{n-1}$. Thus,
  \[
    \mc L^+_{(t,0)}\cap\mc U
      =\begin{cases}
         \{t+3\}\times\mathbb{S}^{n-1}, & t\leq -1, \\
         \{t+1,t+3\}\times\mathbb{S}^{n-1}, & t>-1.
       \end{cases}
  \]
  Correspondingly, the earliest observation time of $\mc L^+_{(t,0)}$ along the timelike curve $\gamma(s)=(s,\omega_0)$ (with $\omega_0\in\mathbb{S}^{n-1}$ fixed) within $\partial M_0$, defined by $s_\gamma(t):=\inf\{s\colon\gamma(s)\in\mc L^+_{(t,0)}\}$, is discontinuous, namely $s_\gamma(t)=t+3$ for $t\leq -1$, and $s_\gamma(t)=t+1$ for $t>-1$.
\end{example}

\subsection{Smooth structure}
\label{SubsecRS}

With the topology of $\mc S$ at our disposal, the space of continuous maps from $\mc S$ into any topological space is well-defined. In order to recover the structure of $\mc S$ as an (open) smooth manifold, we will, in a neighborhood of any point $q\in S$, define a coordinate system by using `earliest observation times' along suitable curves passing through points where $\mc L^+_q\cap\mc U$ is a smooth submanifold.

\begin{lemma}
\label{LemmaRSFinite}
  For $q\in M$ and $p\in\partial M$, the number of different vectors $V\inL^{\mathrm{b}}_q M$ for which $\exp^{\mathrm{b}}_q(V)=p$ is finite.
\end{lemma}
\begin{proof}
  Note that all such $V$ have to be non-zero. If $V$ is such a vector, then $\exp^{\mathrm{b}}_q(s V)=p$ only holds for $s=1$. Thus, it suffices to prove that there are only finitely many rays in $L^{\mathrm{b}}_q M\setminus\{0\}$ whose image under $\exp^{\mathrm{b}}_q$ passes through $p$. Since $(L^{\mathrm{b}}_q M\setminus\{0\})/\mathbb{R}_+$ is a compact space, it suffices to prove that every ray $\mathbb{R}_+ V$ with $\exp^{\mathrm{b}}_q(V)=p$ has a punctured neighborhood consisting of rays whose image under $\exp^{\mathrm{b}}_q$ does not contain $p$. But this follows from Lemma~\ref{LemmaRConj} (using assumption~\itref{ItRConj} of Theorem~\ref{ThmR}).
\end{proof}

\begin{lemma}
\label{LemmaRSTransversal}
  Suppose that $q\in S$ and $p\in\partial M$ are such that for all $V\inL^{\mathrm{b,+}}_q M$ satisfying $\exp^{\mathrm{b}}_q(V)=p$, $(q,V)$ and $p$ are not conjugate. Then there exist a neighborhood $O\ni p$, an integer $N<\infty$, and $N$ pairwise transversal smooth codimension $1$ submanifolds $L_j$ ($j=1,\ldots,N$) of $O$ such that $\mc L_q^+\cap O=\bigcup_{j=1}^N L_j$.
\end{lemma}

See Figure~\ref{FigRSTrans}.

\begin{figure}[!ht]
\centering
  \inclfig{RSTrans}
  \caption{Two different future-directed light rays from $q$ intersecting the boundary at the same point $p$. Under the assumption that $(q,V_j)$ and $p$ are not conjugate for $j=1,2$, the light observation set $\mc L_q^+\cap\partial M$ is, near $p$, the union of two transversally intersecting codimension $1$ submanifolds $L_1,L_2\subset\partial M$.}
\label{FigRSTrans}
\end{figure}

\begin{proof}[Proof of Lemma~\ref{LemmaRSTransversal}]
  Let $\{V_1,\ldots,V_N\}:=(\exp^{\mathrm{b}}_q)^{-1}(p)\capL^{\mathrm{b,+}}_q M$. As in Lemma~\ref{LemmaRConj}, there exists a smooth codimension $1$ submanifold $Z_j\subsetL^{\mathrm{b}}_q M$ containing $V_j$ such that $L_j:=\exp^{\mathrm{b}}_q(Z_j)$ is a smooth codimension $1$ submanifold of $\partial M$, which moreover is spacelike by Lemma~\ref{LemmaGeoANull}. Furthermore, by construction, $\bigcup_{j=1}^N L_j=\mc L_q^+\cap O'$ for a sufficiently small neighborhood $O'$ of $p$.

  If $N=1$, we can take $O=O'$, and the proof is complete. If $N\geq 2$, we first establish the transversality of $L_j$ and $L_k$, $1\leq j\neq k\leq N$, at $p$. Let $\gamma_i(s)=\exp^{\mathrm{b}}_q(s V_i)$, $i=j,k$; we then observe that $T_p L_i$ uniquely determines an outward lightlike ray through $p$, which is necessarily equal to the ray $\mathbb{R}_+W_i$, where $W_i:=\gamma_i'(1-0)\in (T_{\partial M}^+ M)_p$ for $i=j,k$; this uses Lemma~\ref{LemmaGeoANull}. Thus, if $T_p L_j=T_p L_k$, then $c W_j=W_k$ for some $c>0$. But then
  \begin{equation}
  \label{EqRSTransversalEq}
    \exp^{\mathrm{b}}_p(-c s W_j)=\exp^{\mathrm{b}}_p(-s W_k)\in\gamma_j([0,1])\cap\gamma_k([0,1]),\ \ s\in[0,1].
  \end{equation}
  Now for $s=1$, we have $\exp^{\mathrm{b}}_p(-W_k)=q$, but we also have $\exp^{\mathrm{b}}_p(-W_j)=q$ by construction. Thus, $c=1$, and by differentiating the equality in~\eqref{EqRSTransversalEq} in $s$ at $s=1$, we find $-V_j=-V_k$, contradicting $j\neq k$. The conclusion of the Lemma follows if we take $O\subset O'$ to be a sufficiently small neighborhood of $p$.
\end{proof}

Thus, away from a finite union of smooth codimension $2$ submanifolds of $\partial M$, $\mc L_q^+\cap O$ is a smooth codimension $1$ submanifold of $\partial M$. Define the smooth part
\begin{equation}
\label{EqRSLReg}
\begin{split}
  \mc L_q^{\mathrm{reg}} := \big\{ p \in \mc L_q^+\cap\mc U &\colon \text{there exists}\ p\in O\subset\partial M\ \text{open such that}\ \mc L_q^+\cap O\ \text{is a} \\
  &\qquad \text{smooth connected codimension 1 submanifold of}\ O \bigr\}.
\end{split}
\end{equation}
(In the notation of the Lemma~\ref{LemmaRSTransversal}, we have $\mc L_q^{\mathrm{reg}}\cap U=(\mc L_q^+\cap U)\setminus\bigcup_{j\neq k}(L_j\cap L_k)$.)

Fix now any $q\in S$, and denote by $\mu\colon[-1,1]\to\mc U$ a smooth curve in $\partial M$ which is transversal to $\mc L_q^{\mathrm{reg}}$, with $\mu'(s)\neq 0$ for $s\in[-1,1]$, and such that $\mu(0)\in\mc L_q^{\mathrm{reg}}$ and $\mu(s)\not\in\mc L^+_q$ for $s\neq 0$. Consider the set
\[
  R'(\mu):=\{q'\in S\colon \#(\mc L^{\mathrm{reg}}_{q'}\cap\mu([-1,1]))=1,\ \text{with transversal intersection}\}.
\]
While $R'(\mu)$ is neither open nor closed in general, it does contain an open neighborhood of $q$. Therefore, the set
\[
  R(\mu) := \bigcup_{\genfrac{}{}{0pt}{}{R\subset R'(\mu)}{\text{open in}\ S}} R
\]
is a non-empty open neighborhood of $q$. By part~\itref{ItGeoGExpCl} of Proposition~\ref{PropGeoGExpCont}, the earliest observation time
\[
  x^\mu \colon R(\mu)\ni q \mapsto s \in [-1,1],\ \ \text{where}\ \mu(s)\in\mc L^+_q,
\]
is a smooth function on $R(\mu)$, and $x^\mu(q)=0$. (We stress that $R(\mu)$ and $x^\mu$ are well-defined given the topology of $S$ and the smooth structure of $\mc U$.) We aim to show that suitable families of such functions $x^\mu$ give local coordinates near $q$. The key step is to show that there is always a large supply of curves $\mu$ for which $x^\mu$ is non-degenerate at $q$; more precisely:

\begin{lemma}
\label{LemmaRSGrad}
  Fix $q\in S$, and let
  \begin{align*}
    \mathfrak{M}:=\{\mu\colon[-1,1]\xrightarrow{\cC^\infty}\mc U \colon \mu\ &\text{is transversal to}\ \mc L_q^+,\ \mu'(s)\neq 0\ \text{for}\ s\in[-1,1], \\
      & \qquad\qquad \mu(0)\in\mc L_q^{\mathrm{reg}},\ \mu(s)\not\in\mc L_q^+\ \text{for}\ s\in[-1,1]\}.
  \end{align*}
  Then
  \[
    \bigcap_{\mu\in\mathfrak{M}} \ker(dx^\mu|_q) = \{0\} \subset T_q M.
  \]
\end{lemma}

We give an analytic proof here, arguing by contradiction. The arguments in \S\ref{SubsecRO} below provide a different, more geometric proof.

\begin{proof}[Proof of Lemma~\ref{LemmaRSGrad}]
  Let $(-1,1)\ni r\mapsto q(r)\in S$ be a smooth path with $q(0)=q$ and $V:=q'(0)\neq 0\in T_q M$. Suppose that
  \begin{equation}
  \label{EqRSGrad0}
    dx^\mu(V)=0\ \text{for all}\ \mu\in\mathfrak{M};
  \end{equation}
  equivalently, for all $\mu\in\mathfrak{M}$, the curve $r\mapsto\tilde\mu(r)$ defined by $\{\tilde\mu(r)\}=\mc L_{q(r)}^+\cap\mu([-1,1])$ for small $r$, so $\tilde\mu(0)=\mu(0)$, satisfies $\tilde\mu'(0)=0$.

  Let now $O\subset\partial M$ be an open neighborhood of a point in $\mc L^{\mathrm{reg}}_q$, as in \eqref{EqRSLReg}. Pick any non-empty $O'\Subset O$, and denote $L:=\mc L_q^{\mathrm{reg}}\cap O'$. Since $\mc L_q^{\mathrm{reg}}\cap O$ is smooth of codimension $1$, we can pick a smooth open map
  \[
    \underline\mu \colon L \times (-2,2) \to O
  \]
  such that $s\mapsto\underline\mu(p,s)$ is a curve with $\underline\mu(p,0)=p$, transversal to $\mc L_q^{\mathrm{reg}}\cap O$, and so that $\underline\mu$ is a diffeomorphism onto its image $O''\subset\mc U$.

  For small $r$, the preimage $\underline\mu^{-1}(\mc L^+_{q(r)}\cap O'')$ is the graph of a smooth function $f(r,\cdot)\colon L\to(-2,2)$ and in fact $f\colon(-\epsilon,\epsilon)\times L\to(-2,2)$ (shrinking $\epsilon>0$ if necessary) is smooth. (These are consequences of Proposition~\ref{PropGeoGExpCont}.) Furthermore, $f(0,\cdot)\equiv 0$. Since we are assuming that \eqref{EqRSGrad0} holds, so $\partial_r f(0,p)\equiv 0$, the tangent space
  \[
    T(r,p) := T_{\underline\mu(p,f(r,p))}\mc L^+_{q(r)}
  \]
  is $r^2$-close to $T(0,p)=T_p\mc L^+_q$, uniformly for all $p\in O'$, hence the same is true for the unique future lightlike, outward pointing ray $\ell(r,p)\subset T(r,p)^\perp$, see Lemma~\ref{LemmaGeoANull}.

  Let now $V_1,V_2\inL^{\mathrm{b,+}}_q  M$ be two \emph{distinct} non-zero tangent vectors such that $p_j:=\exp^{\mathrm{b}}_q(V_j)\in L$, $j=1,2$, then $\mathbb{R}_+\cdot\frac{d}{ds}\exp^{\mathrm{b}}_q(s V_j)|_{s=1-0}=\ell(0,p_j)$. Denote by $W(r,p)$ a generator of $\ell(r,p)$ which depends smoothly on $(r,p)\in(-\epsilon,\epsilon)\times L$, and which is $r^2$-close to $W(0,p)$. Then the images of the two broken null-geodesics $s\mapsto\exp^{\mathrm{b}}_{p_j}(-s W(0,p_j))$ for $j=1,2$ intersect \emph{cleanly} at $q$. But this implies that the point $q(r)$ is the \emph{unique} element near $q$ of the set of intersections of the broken null-geodesics $\exp^{\mathrm{b}}_{\underline\mu(p_j,f(r,p_j))}(-s W(r,p_j))$, $j=1,2$, and moreover $q(r)$ depends smoothly on $f(r,p_j)$ and $W(r,p_j)$. The properties of $f$ and $W$ therefore imply that $q(r)$ is $r^2$-close to $q=q(0)$, contradicting the assumption $q'(0)\neq 0$ and completing the proof.
\end{proof}

In particular, for every $q\in S$, there exist $(n+1)$ curves $\mu_j\in\mathfrak{M}$ such that the set $\{dx^{\mu_j}\colon j=0,\ldots,n\}$ is linearly independent at $q$, and therefore $(x^{\mu_j})_{j=0,\ldots,n}$ is a smooth local coordinate system near $q$. However, only knowing the collection $\mc S$ of light observation sets, it is not immediately clear how to determine if a family $\mu_j$, $j=0,\ldots,n$, has this property. We thus argue indirectly: define a subalgebra
\[
  \mc C\subset\mc C^0(S)
\]
by declaring that $f\in\mc C^0(S)$ is an element of $\mc C$ if and only if for every $q\in S$, there exist an open neighborhood $U\ni q$ and curves $\mu_i\in\mathfrak{M}$ (in the notation of Lemma~\ref{LemmaRSGrad}) for $0\leq i\leq n$ such that $U\subset\bigcap_{i=0}^n R(\mu_i)$, and a smooth function $F\colon\mathbb{R}^{n+1}\to\mathbb{R}$ so that
\begin{equation}
\label{EqRSFunc}
  f(q')=F(x^{\mu_0}(q'),\ldots,x^{\mu_n}(q')), \ \ q'\in U.
\end{equation}
By the arguments presented in this section, $\mc C=\cC^\infty(S)$, and hence we have recovered the algebra of smooth functions on $S$ from the family of sets $\mc S$.

Lastly then, a set of $n+1$ curves $\mu_i$, $0\leq i\leq n$, for which \emph{every} element of $\mc C$ can be expressed in a neighborhood $U$ of $q$ in the form~\eqref{EqRSFunc} gives rise to a local coordinate system $(x^{\mu_j})_{j=0,\ldots,n}\colon U\to\mathbb{R}^{n+1}$. This completes the reconstruction of $S$ as a smooth manifold.

\subsection{Conformal structure}
\label{SubsecRC}

The reconstruction of the conformal structure of $S$ is straightforward: if $q\in S$, let $V\inL^{\mathrm{b,+}}_q M$ be such that $p=\exp^{\mathrm{b}}_q(V)\in\mc L_q^{\mathrm{reg}}$, and put $L:=T_p(\mc L_q^+\cap\mc U)$. Consider the set
\[
  Q = \bigl\{ \mu\colon(-1,1)\xrightarrow{\cC^\infty} S \colon \mu(0)=q,\ p\in\mc L_{\mu(r)}^+\ \text{for all}\ r,\ T_p(\mc L_{\mu(r)}^+\cap\mc U)=L \bigr\}
\]
of all paths which have the same outgoing future null ray at $p$, see Lemma~\ref{LemmaGeoANull}. Then $\{\mu'(0)\colon\mu\in Q\}=\mathbb{R} V\in T_q M$ recovers a 1-dimensional lightlike subspace of $T_q S$. Repeating this procedure for all points $p\in\mc L_q^{\mathrm{reg}}$, and noting that $\mc L_q^{\mathrm{reg}}\subset\mc L_q\cap\mc U$ is dense by Lemma~\ref{LemmaRSTransversal}, we can reconstruct an open subset of the light cone $L_q M\subset T_q M$. But $L_q M$ is a real-analytic submanifold of $T_q M$, hence this determines $L_q M$ uniquely. Since $q\in S$ was arbitrary, this proves that we can recover $L_S M$, hence the conformal structure of $S$. This finishes the proof of the first part of Theorem~\ref{ThmR}.

\subsection{Time orientation}
\label{SubsecRO}

In order to recover the time orientation of $S$ when we are given the conformal structure of $\mc U$ as well as its time orientation, we analyze the dependence of the boundary intersection point of a broken null-geodesic on its initial point:

\begin{lemma}
\label{LemmaRODep}
  Suppose $(-1,1)\ni r\mapsto (q(r),V(r))\inL^{\mathrm{b,+}} M$ is a smooth path such that $p(r):=\exp^{\mathrm{b}}_{q(r)}V(r)\in\partial M$. Let $\gamma(r,s):=\exp^{\mathrm{b}}_{q(r)}(s V(r))$ and $\gamma(s):=\gamma(0,s)$. Then
  \begin{equation}
  \label{EqRODep}
    g(q'(0),V(0)) = g(p'(0),\gamma'(1)).
  \end{equation}
\end{lemma}
\begin{proof}
  The values $0<\mathsf{s}_1(r)<\cdots<\mathsf{s}_k(r)=1$ of $s$ for which $\gamma(r,s)=0$ are smooth functions of $r$ for $r$ small, likewise the boundary intersection points $\mathsf{p}_j(r)=\gamma(r,\mathsf{s}_j(r))$; see also the discussion preceding Proposition~\ref{PropGeoGExpCont}. Define $\mathsf{s}_0(r):=0$ and $\mathsf{p}_0(r):=q(r)$. For $j=0,1,\ldots,k-1$, we then have $\int_{\mathsf{s}_j(r)}^{\mathsf{s}_{j+1}(r)} |\partial_s\gamma(r,s)|_g^2\,ds=0$ for all $r$, hence by differentiation in $r$, using that $\partial_s\gamma(r,s)$ is null for all $s$, and further using that $\gamma(r,s)$ is a null-geodesic for $s\in(\mathsf{s}_j(r),\mathsf{s}_{j+1}(r))$,
  \begin{align*}
    0 &= \int_{\mathsf{s}_j(r)}^{\mathsf{s}_{j+1}(r)} g(\partial_s\gamma(0,s),D_s\partial_r\gamma(r,s)|_{r=0})\,ds \\
      &= g\bigl(\gamma'(\mathsf{s}_{j+1}(r)-0), \mathsf{p}_{j+1}'(0)\bigr) - g\bigl(\gamma'(\mathsf{s}_j(r)+0), \mathsf{p}_j'(0)\bigr).
  \end{align*}
  Summing these identities and using that $\gamma'(\mathsf{s}_j(r)+0)-\gamma'(\mathsf{s}_j(r)-0)\perp T\partial M\ni\mathsf{p}_j'(0)$, all but the first and last terms cancel, and we obtain~\eqref{EqRODep}.
\end{proof}

Let now $(-1,1)\ni r\mapsto q(r)\in S$ be a timelike path; we show that one can determine whether $q$ is \emph{future} timelike:

\begin{prop}
\label{PropRO}
  Let $p(r)\in\mc L_{q(r)}^{\mathrm{reg}}\cap\mc U$ be a smooth path, and denote by $N\in T_{p(0)}\partial M$ the future-directed unit normal to the spacelike hypersurface $T_{p(0)}\mc L_{q(0)}^+$. Then $q$ is future timelike if and only if $g(p'(0),N)<0$.
\end{prop}

We stress that this criterion only uses the conformal structure and time orientation of $(\mc U,g|_{\mc U})$.

\begin{proof}[Proof of Proposition~\ref{PropRO}]
  We claim that $p(r)=\exp^{\mathrm{b}}_{q(r)}(V(r))$ with $V(r)\inL^{\mathrm{b,+}}_{q(r)}M$ smooth in $r$. By definition of $\mc L_{q(0)}^{\mathrm{reg}}$, there exists a unique $V(0)\inL^{\mathrm{b,+}}_{q(0)}M$ such that $p(0)=\gamma(1)$ for $\gamma(s):=\exp^{\mathrm{b}}_{q(0)}(s V(0))$. Similarly to the proof of Lemma~\ref{LemmaRSGrad}, let $W(r)$ denote a generator of the future-directed outward pointing null ray orthogonal to $\mc L^+_{q(r)}\cap\mc U$, see Lemma~\ref{LemmaGeoANull}, so that $W(0)=\gamma'(1)$. Since $q$ is timelike, the intersection of the broken null-geodesic $\mu_r(s):=\exp^{\mathrm{b}}_{p(r)}(-s W(r))$ with $q$ is unique (if necessary shrinking the interval that $r$ takes values in) and clean; therefore we can choose a smooth function $s(r)$ such that $\mu_r(s(r))=q(r)$, with $s(0)=1$. But then $V(r)=-\mu_r'(s(r))/s(r)$ is smooth, as claimed.

  We can now apply Lemma~\ref{LemmaRODep} and use that the orthogonal projection of $\gamma'(1)=W(0)\in T_{p(0)}M$ to $T_{p(0)}\partial M$ is a positive multiple of $N$; since $V(0)$ is future-directed, we conclude that $q'(0)$ is future timelike iff $g(p'(0),N)<0$, proving the proposition.
\end{proof}

This finishes the proof of the second part of Theorem~\ref{ThmR}.

\bibliographystyle{alpha}
\input{bdylo.bbl}

\end{document}

redpajama_set_name: RedPajamaArXiv


text:
Kevin Harvick wins Stage 1 of Bank of America 500
By Daniel McFadinOct 8, 2017, 2:11 PM EDT
CONCORD, North Carolina — Kevin Harvick won Stage 1 of the Bank of America 500 at Charlotte Motor Speedway for his fourth stage win of the year.
Harvick led the final 39 laps of the stage after taking the lead from Chase Elliott on a restart.
The top 10 through 90 laps is Harvick, Chase Elliott, Jamie McMurray, Kyle Larson, Kyle Busch, pole-sitter Denny Hamlin, Matt Kenseth, Kurt Busch, Daniel Suarez and Jimmie Johnson.
Martin Truex Jr. finished 12th in the stage. It is the first time he has finished outside the top 10 in Stage 1 since the July race at Daytona.
The stage was slowed by two cautions. The first was a competition caution on Lap 35.
Kyle Larson missed his pit stall under caution and had to pit again. He restarted 31st before charging his way to the top five.
The second caution waved on Lap 44 for a crash involving Clint Bowyer and Ryan Newman. The crash was triggered when Bowyer bounced off the wall on the backstretch and into Newman, who was turned into the wall. Newman's No. 31 Chevrolet briefly erupted in flames, but Newman was able to exit the car quickly unharmed.
"I know I got tight underneath him and washed up and checked up," Newman told NBC. "I don't know if he turned me on purpose. He had a right too. It was either my mistake or his mistake. We were both going to same piece of real estate."
Harvick took the lead on the ensuing Lap 52 restart.
Hamlin led 40 laps and Elliott led 11.
The second stage of the race will end on Lap 180.

redpajama_set_name: RedPajamaCommonCrawl


I hope that your 2018 has been awesome so far and that January brought you joy!
As some of you may have noticed, I had different plans for January: I've lost 2 weeks to massive water damage in my apartment and the aftermath of that catastrophe, which wrecked my planned coverage of Starfinder-materials. A lot of the coverage of this month either was based on already finished tests and/or my unpublished backlog. Rest assured that there will be more Starfinder-material in February, now that I slowly start getting my bearing once more.
That being said, if you're one of the many, many folks that check out my site and reviews on a regular base, please take a look at my patreon.
Right now, this whole reviewing gig takes up almost the equivalent on a proper dayjob in the time that I devote to it (Who am I kidding? More than that, I just don't consider reading the books as work…just the testing, writing, etc.); the water damage has made me very cognizant of the fact that it, right now, amounts to a 0-sum game that does not allow me to accumulate any savings. While I am not wholly reliant on reviewing, my day-job doesn't pay well either. I the long run, I'll have to reevaluate how much time I can devote to reviewing if there is not an upswing in support. I absolutely believe that this is possible – the gaming community is absolutely awesome!
So yeah, if you want to continue seeing this level of commitment and coverage, please contemplate supporting my patreon – even a small donation will help and if everyone who checks out this site donated just a single buck per month, I could do this full-time.
Straight Path Games – What Is…Incorporeal?
Straight Path Games – How Do I…Use Magic Items?