A reflection approach to the broken ray transform
We reduce the broken ray transform on some Riemannian manifolds (with corners) to the geodesic ray transform on another manifold, which is obtained from the original one by reflection. We give examples of this idea and present injectivity results for the broken ray transform using corresponding earlier results for the geodesic ray transform. Examples of manifolds where the broken ray transform is injective include Euclidean cones and parts of the spheres . In addition, we introduce the periodic broken ray transform and use the reflection argument to produce examples of manifolds where it is injective. We also give counterexamples to both periodic and nonperiodic cases. The broken ray transform arises in Calderón’s problem with partial data, and we give implications of our results for this application.
Key words and phrases:Broken ray transform, X-ray transform, Calderón’s problem, inverse problems
2010 Mathematics Subject Classification:Primary 53C65, 78A05; Secondary 35R30, 58J32
Suppose we have an unknown compactly supported continuous function in the upper half plane and we know its integrals over all broken lines in the upper half plane, which reflect at according to the usual law of geometrical optics. We can deduce the function from these integrals by reflecting the half plane to fill the entire plane and unfolding the broken rays into straight lines. If we let for , then we may reconstruct the integral of over any straight line in the plane. By injectivity of the Radon transform in the plane, we can deduce the original function from this information. In this article we generalize this reflection argument to show injectivity of the broken ray transform in various domains.
Let be an -dimensional compact Riemannian manifold with boundary. We assume that the boundary is a disjoint union of , , and such that and are open and in the topology of .
We consider broken rays to be piecewise geodesic paths on such that
starts and ends in the set ,
is a geodesic in , and
is reflected on according to the usual reflection law: the angle of incidence equals the angle of reflection.
If convenient, we may also allow reflections on ; such paths can be constructed by concatenating a finite number of broken rays as defined above. Since all broken rays have endpoints in the set , we call it the set of tomography.
We ask the following questions: If the integral of an unknown real valued function on is known over all broken rays, can be reconstructed? If yes, is the reconstruction stable? How do answers to these questions depend on the regularity assumptions on , and ?
To answer these questions, we reduce the problem to injectivity and regularity of the geodesic ray transform on Riemannian manifolds via reflections. This can be done most naturally on manifolds with corners as discussed and proven in Section 4. For manifolds with smooth boundary, more steps have to be taken, and they are outlined in Section 5.
We define the geodesic ray transform and the broken ray transform as follows:
For a manifold with boundary we denote the set of all geodesics joining boundary points by . For two classes of functions and any we define the attenuated (geodesic) ray transform by
when is a geodesic in with unit speed.
For a set of tomography , we denote the set of broken rays from to by (allowing reflections on ). We define similarly the attenuated broken ray transform by
If attenuation nor the word ‘attenuated’ is not mentioned, the attenuation is assumed to vanish identically.
We also study the periodic broken ray transform where the entire boundary is reflecting and integrals of the unknown function are known over all periodic broken rays. The precise definition is the following:
Let be a Riemannian manifold with boundary. Let be the set of all periodic broken rays in . The mapping , , is the periodic broken ray transform.
Using the reflection approach, we show that the broken ray transform is injective on the following manifolds (regularity requirements for functions vary, but suffices in each case):
Quarter of the sphere for where the set of tomography if half of the boundary. (See Proposition 26.)
The two dimensional hemisphere where the set of tomography is slightly larger than half of the boundary. (See Proposition 26.)
An octant of the sphere for the periodic transform. (See Proposition 30.)
The cube , , for the periodic transform. (See Proposition 31.)
Eskin  reduced the recovery of an electromagnetic potential from partial data to injectivity of the broken ray transform, and showed the transform to be injective in a Euclidean domain with convex reflecting obstacles. The broken ray transform has recently been studied in its own right [14, 12]. This research has been motivated by the fact that Kenig and Salo  reduced Calderón’s problem with partial data to the injectivity of the broken ray transform. We will discuss this in more detail in Section 1.1 below.
The recovery of a function from its ray transform is a well understood problem in a Euclidean domain (see textbooks [11, 22]), and there are also a number of results on Riemannian manifolds (of which we mention [4, 26, 30, 18]) and also in greater generality (see e.g. [9, 3]). The broken ray transform, however, is much less studied, which makes it appealing to try to reduce broken ray problems to the usual ray tomography.
It should be noted that while the geodesic ray transform is a good model for measuring the attenuation coefficient in a material with light, the broken ray transform is not a very good model if the light ray is allowed to reflect. After a few reflections the signal is essentially lost, and reconstruction methods using broken rays with one reflection only are more appropriate for this application. The model with one reflection (in the interior of the domain) is known as the V-line Radon transform [21, 2]. In the case of multiple reflections of light it is more appropriate to use the radiative transfer equation to model propagation of light.
This paper is organized as follows. In Section 1.1 we recall the relation between Calderón’s problem and the broken ray transform and show how to translate the results in this paper to results for Calderón’s problem. Section 2 gives examples of the reflection construction by proving injectivity of the broken ray transform in Euclidean cones. More general examples on manifolds are given in Section 3. To prove these more general examples, reflected manifolds are constructed and studied in Section 4, and a generalization of the result for Euclidean cones is given in Theorem 16. The examples in Section 3 are based on this theorem. The results up to this point require that the manifold has a corner at . A method for removing this restriction is presented in Section 5. In Section 6 we give examples (and counterexamples) of specific manifolds where the broken ray transform is injective. In Section 7 we demonstrate by example that a similar reflection approach can also be used for the periodic broken ray transform, where the entire boundary is reflecting and one integrates over periodic broken rays.
1.1. Relation to Calderón’s problem
As mentioned above, our main motivation for the study of the broken ray transform comes from the Calderón problem with partial data. The recent result by Kenig and Salo [17, Theorem 2.4] states roughly the following if the broken ray transform on with set of tomography is injective: the partial Cauchy data for the Schrödinger equation on a manifold determines the potential uniquely, provided that contains a tubular part and the inaccessible part of the boundary of is contained in . The Calderón problem can be reduced to the corresponding problem for the Schrödinger equation. For basic results for the Calderón problem we refer to the review article  and references therein.
As an example, we state the result for Calderón’s problem arising from injectivity of the broken ray proven in Proposition 6.
Let be any constant, and define the gutter , where . Let be a bounded domain such that for some we have . Suppose that is outside the line . (One such domain is sketched in Figure 1 for .)
Let be the inaccessible part of the boundary and denote the accessible part by . Then the partial Cauchy data
determines the potential uniquely in . (The inaccessible part of the boundary is shadowed in Figure 1.)
This implies that the similarly partial Cauchy data for the Calderón problem determines a conductivity uniquely.
By the result [17, Theorem 2.4] it suffices to show that the broken ray transform is injective on the transversal manifold with all constant attenuations, where the set of tomography is . But this is done in Proposition 6(2) below; the parameters and are related by . Although we only prove Proposition 6(2) without attenuation, the same proof is valid for any constant attenuation as noted in Remark 19. ∎
Theorem 3 was proven for (half space) by Isakov  by another reflection method. For it was proven by Kenig and Salo ; in this case it suffices to study broken rays without reflections. Our result generalizes the previous ones to . Injectivity results for the broken ray transform can be turned into partial data results for Calderón’s problem in corresponding tubular domains; Theorem 3 is an example of this.
A partial version of Theorem 3 remains true if the corner of the transversal domain is smoothed out as follows. Let be a transversal domain satisfying the assumptions of the theorem. Suppose then is a subdomain such that for some we have . Now . As demonstrated at the end of Section 2.1, if the broken ray transform of is known, one can determine outside . Thus for a domain with transversal domain as in the theorem the partial Cauchy data with determines outside the tube . The domain can have a smooth boundary, unlike .
Full recovery is possible if the smoothened tip of the cone is in not reflective but available for measurements; see Lemma 5 below.
2. First examples
We present some examples in the following proposition which demonstrate the idea that we wish to generalize. We begin with a lemma that contains a general observation. Here we denote by the piecewise continuous functions from to .
Suppose the broken ray transform on a domain is injective on with some set of tomography . Then it is also injective on any subdomain (now on ) with a new set of tomography .
The lemma is true also for manifolds and functions and the proof is the same.
Proof of Lemma 5.
Suppose integrates to zero over all broken rays. Define then by letting on and on . We clearly have .
Take now any broken ray in . Intersecting it with gives segments which are broken rays in (endpoints on ). Using this decomposition and the definition of , we observe that the broken ray transform of vanishes. By assumption this implies and hence . ∎
The broken ray transform is injective in the following Euclidean domains for functions in (and thus also on ):
A domain with on a cone with opening angle , , and any set of tomography . For example in polar coordinates
for continuous functions , where .
The previous example works for all without the restriction . (That is, the opening angle may be anything in the range .)
Also higher dimensional cones
for continuous functions and parameters such that .
The previous example works for all even without the restriction that . (That is, the opening angle may be anything in the range .)
These injectivity results are also true for any subdomain with new reflecting set . In particular, the cone need not contain a nonsmooth tip.
The last remark follows from Lemma 5.
(1) Reflect over one side of the cone. Then reflect this new copy of over the other side of the cone, and carry on until there are copies of the original domain. Because the opening angle is , the total angle adds up to and the copies of form a domain . This reflection process is illustrated in Figs. 2 and 3; the domain in Fig. 2 with is copied and reflected to constitute the domain in Fig. 3. Note that if one continues the construction by reflecting the last (th) copy of , one ends up with the original domain in its original position.
Let be the natural projection map that undoes the copying, rotating and reflecting done in the construction of . Let be any function. We define a reflected version of by letting .
Take any line in that does not meet the origin (for example the one in Fig. 3). Let ; then is a broken ray in (the broken ray in Fig. 2 corresponds to the line in Fig. 3 in this way). Because of this correspondence between and we write . This correspondence of lines and broken rays is illustrated in Figs. 2 and 3. (Note that for each broken ray in that does not hit the tip of the cone there are lines in .)
Since we have
we may construct the Radon transform of from the broken ray transform of . In particular, vanishing broken ray transform of in implies that the Radon transform of vanishes.
Since the Radon transform is injective, vanishing broken ray transform of implies that and so .
(2) Just like above, reflect and copy the domain in the plane. The plane cannot be filled as nicely, but it is enough to construct a cone of opening angle at least . Now does not cover all angles as in part (1) above, but this problem can be bypassed.
Let . We cover the angle left out by by a compact cone with radius centered at the same point as copies of (the origin). This construction is demonstrated in Figs. 4 and 5 (analogously to Figs. 2 and 3).
For , define in as above and extend by zero to . Suppose the broken ray transform of vanishes.
As in the proof of part (1), the vanishing broken ray transform of in implies that the integral of over a line is zero whenever . (One such line is drawn in Fig. 5.) By Helgason’s support theorem [11, Theorem 2.6] (and a little mollification argument, see e.g. [14, proof of proposition 5]) vanishes outside and thus especially in . Thus the original function vanishes in .
(3) The case is literally part (1) and others can be reduced to it.
Let be the Grassmannian of all two dimensional subspaces of . For symmetry reasons every broken ray in the cone is confined to some . For any part (1) shows injectivity in the planar cone ; guarantees that the opening angle is . Since , the broken ray transform is injective in .
(4) Follows from part (2) just like in part (3) follows from part (1). ∎
In the above proposition parts (1) and (3) follow trivially from (2) and (4), but we present them separately since their proofs are somewhat different. In particular, (2) and (4) are based on Helgason’s support theorem.
The broken ray transform is also injective on unbounded cones of the types (1) and (3) with , since the Radon transform is injective in the whole plane. The set of tomography may be taken to be “at infinity” in the sense that broken rays are allowed to tend to infinity. Integrability assumptions are then needed for the unknown functions.
It is important to notice that the gluing in Proposition 6 was done along flat parts of the boundary (line segments in the plane). A particular case of Proposition 6(3) is the half space (with ); the reflection is simply obtained by and there are no corners (or the corners have angle ). We will reflect and glue together manifolds in Section 4 below, and flatness of the gluing boundary will lead to regularity of the reflected manifold (see Lemma 13). We recall that in Isakov’s reflection method for the Calderón problem  reflection is made along a part of the boundary which is flat (hyperplane) or can be conformally flattened (sphere).
The recent result by Hubenthal in the square  heavily relies on the geometry of the square: reflections are done at straight lines and corners have angle . By Proposition 6(1) the broken ray transform is injective in the square, provided that the set of tomography contains two adjacent edges of the square. Similarly the broken ray transform is injective in any polygon if the reflecting part of the boundary contains at most two adjacent edges. Although these results are different in their formulation and methods of proof, the underlying geometrical structure of the square is heavily relied on.
The shape of the domain in dimensions three and higher can be other than the cone in Proposition 6(5). For example, if is a cube with three adjacent faces as the set of tomography, eight copies of it can be glued together to form a bigger cube in a fashion similar to gluing four squares to form a bigger square in dimension two. (Such a construction is used to prove Proposition 31.) The correspondence between lines and broken rays is the same. We do not elaborate on all the possibilities here; we only wish to present the idea in a fair amount of generality.
In the discussion below we will focus on the analogue of the half space. Corners can be allowed, but for the sake of simplicity we shall not allow them. It is the author’s belief that if the corners add up nicely as in Proposition 6(1), Theorem 16 remains true. The technical difficulty lies in the fact that one needs some kind of “corner normal coordinates” at a corner point of a manifold. For more general corners one needs something to replace Helgason’s support theorem in the proof of Proposition 6(2). Support theorems as simple and powerful as Helgason’s seem not to be available for general manifolds.
The examples above were such that the domain was constructed as a submanifold of a particularly nice domain in so that the construction reduces the problem to a planar one. The examples in Section 6 are also of this type. If we start with an arbitrary manifold, the resulting manifold is not generally any simpler than the one we started with.
If the angle in Proposition 6(2) is a reflex angle (between and ), there is no need for a reflection construction. Helgason’s support theorem immediately gives injectivity for the broken ray transform, and reflected rays need not be considered at all. This observation holds true whenever the reflector is concave.
The cone in Proposition 6(2) need not have an angle. If the domain looks like the domain of Proposition 6(2) outside some neighborhood of the origin, we can reconstruct outside some (possibly larger) neighborhood of the origin. This can be done with the same method; Helgason’s support theorem tells that vanishes outside the convex hull of the set (and its copies) where the boundary of is not conical.
3. Two applications on manifolds
We use Theorem 16 to give two theorems of injectivity of the broken ray transform on in a fairly large class of manifolds. In brief, Theorem 16 tells that injectivity of the broken ray transform on a manifold can be reduced to the injectivity of the geodesic ray transform on a reflected manifold in analogue to the half plane example given in the beginning of this article. The notation and necessary results are given in Section 4 below. The purpose of this section is to motivate the general construction. More examples are given in Section 6.
Let satisfy the following assumptions:
is a smooth Riemannian surace with corners, and the boundary is a disjoint union of the sets , , and .
The open smooth boundary components and meet orthogonally at .
is strictly convex and is -flat.
The local boundary defining functions of near can be chosen to be -even at .
For any two points on (but not both of them on ) and a chosen parity (odd or even), there is a unique broken ray in with the chosen parity (even or odd number of reflections) joining the points, and this geodesic depends smoothly on the endpoints.
If both endpoints lie in , the geodesic is contained in .
The normal derivatives of odd orders with respect to endpoints of the geodesic vanish at .
Then the broken ray transform in is injective with set of tomography for smooth functions in which are -even at .
By Theorem 16 it is enough to show that the geodesic ray transform is injective on the Riemannian manifold in the class .
Let be a compact Riemannian manifold with boundary and suppose that . Assume satisfies the following:
is -even at .
If , the set has zero measure.
The level set is strictly convex in for all .
Then the broken ray transform is injective in the class in , when the set of tomography is .
To prove the theorem, we need the following result.
Theorem 9 (Corollary in ).
Let be a compact Riemannian manifold of dimension at least 3 with boundary embedded in a manifold . Assume there is a function such that , , on . Let . Assume furthermore that is strictly convex for all and has zero measure.
Then the geodesic ray transform is injective in the class on the manifold .
Proof of Theorem 8.
4. Reflected manifolds
The key idea in the proof of Proposition 6 was to glue together copies of the original conical domain . The most simple case was when was part of a half space and the reflecting part of lay on the boundary of the half space. In this case two copies of the original domain could be glued together to form .
Similar reflecting and gluing can be done for Riemannian manifolds. We focus here on the case analogous with the Euclidean half space. The construction for more complicated Euclidean domains presented above can be generalized in the same fashion, but for the sake of simplicity we omit them here. The analogue of Proposition 6 for Riemannian manifolds can be used to show injectivity of the broken ray transform on some Riemannian manifolds.
4.1. Construction of reflected manifolds
In the following is a smooth compact Riemannian manifold with boundary and is a closed subset with and the boundary of meets orthogonally. We reflect the set with respect to . The manifold constructed below is this reflection of . The construction is illustrated if Figure 6.
Let , , , and . The “boundary” (it is not the boundary in the topology of ) of is the (disjoint) union . Now is a topological manifold with boundary with and the boundary is . Furthermore, is a smooth manifold with corners with interior , smooth part of the boundary , and nonsmooth part of the boundary . The sets , and are the strata of of depths 0, 1 and 2, respectively, in the sense of .
By we always mean the boundary of in the topology of . Thus, if we identify with , we have .
The higher depth strata of are empty (there can only be a corner in one direction), but in the example of Proposition 31 strata of all possible depths appear. We bound the depth of strata by 2 for technical simplicity. As manifolds with smooth boundary are more convenient to work with, we will reduce the depth bound to 1 in Section 5.
We define as an identical copy of with all labels changed, and glue and together along and to form a manifold with (smooth) boundary.
Let be the natural bijection. We define the relation on the disjoint union by letting if , or and , or and . This is obviously an equivalence relation, and the quotient space is well defined. We denote by the image of (or ; the image is the same) under the quotient map. We consider , , to be a subspace of in the natural way. We define to be the natural injection.
It is geometrically rather obvious that is an -dimensional topological manifold with boundary. The boundary of is a disjoint union of , , and , and the interior is a disjoint union of , , and . Using Riemannian boundary normal coordinates at we also turn it into a smooth manifold with boundary in a natural way. In we may proceed similarly, but the model space is which after reflection and gluing becomes ; we have as expected.
Due to the use of Riemannian boundary normal coordinates the transition maps are smooth at , and is indeed a smooth manifold with boundary.
The natural projection map defined by and for each is a covering map. Identifying with and and writing as the projection, we have the obvious property
The projection can be used to pull back (scalar, vector, and tensor) functions from to the reflected manifold ; we define for any the reflected version and similarly for higher rank tensors.
Special care must be taken since is not smooth but only continuous on . This is related to the fact that not all choices of boundary charts at give a smooth structure. Some additional conditions thus need to be satisfied to guarantee that is smooth if is. This issue is considered in Section 4.2 in more detail.
The above construction does not use the smoothness of the manifold and its metric . If we instead equip with differentiable structure and take with , then is naturally equipped with a structure. The condition ensures that geodesics on do not branch, and the boundary normal coordinates actually provide coordinates.
We remark that the set is not a manifold with boundary, since it is not locally diffeomorphic to at points in . At these points the proper model space is , which makes a manifold with corners in the sense defined in . Theorem 16 for manifolds with corners will be generalized to manifolds without corners in Section 5.
We wish to point out that both precomposition (for functions on a manifold) and inverse of postcomposition (for curves) of an object with the projection are denoted by a tilde. An object with tilde should therefore be understood as the natural corresponding object (or one of them in the case of curves) on the reflected manifold .
4.2. Regularity of the reflected manifold
We keep the assumptions made in the beginning of Section 4.1 regarding , and but the smoothness requirement is only that with . Lemma 13 below demonstrates the correspondence between the properties of and its reflected version .
Let and . A function is -even at if for all odd .
A rank two tensor (written in boundary normal coordinates near the boundary) of class is -even at if the functions and for are -even at in the above sense and for all even . In particular, is called -flat if the metric is -even at .
The above definition can be easily extended to tensors of any rank, but we do not need to consider ranks other than zero and two here. The definitions are given so that a tensor field on is -even at if and only if is a tensor field on . This correspondence is the basis of Lemma 13 below.
Let be an closed subset of as in Section 4.1. A set (boundary in the topology of ) is evenly (resp. oddly) strictly broken ray convex in , if for any two points on there is a broken ray with an even (resp. odd) number of reflections (possibly zero) connecting the two points such that the interior of the broken ray is in .
Regularity of functions and metrics on the original manifold and the reflecting manifold correspond in the following way (here ):
is 1-flat if and only if the second fundamental form vanishes on .
If and is -flat, then .
If and is -even at , then .
If , then .
If is (strictly) evenly and oddly broken ray convex, then is geodesically convex.
If is strictly convex , then is strictly convex.
If is convex, then is convex.
Suppose and are . If is strictly convex, there are no geodesics tangent to .
If is strictly concave, geodesics tangent to branch.
is simple if the assumptions listed in Theorem 7 hold (with the word ‘surface’ replaced by ‘manifold’).
In the parts (8–9) a geodesic means a locally length minimizing curve, since the geodesic equation does not make sense on when is not 1-flat.
(1) Use boundary normal coordinates with the first coordinate as the normal direction to the boundary. In these coordinates is constant and when . Thus it suffices to study for all .
In these coordinates is the outward unit normal vector. Since is constant in these coordinates, for two vectors and tangent to the boundary at a boundary point we have for the second fundamental form
Thus vanishes on if and only if for all .
(8) Let be the map which reflects with respect to in the natural way.
Suppose is a geodesic which meets tangentially at . The intersection points of and cannot accumulate at unless ; the points of reflection of a broken ray of finite length near a strictly convex part of the boundary cannot accumulate. The proof of this statement is too long to be included here; see  for proof and explanation of the assumption.
There are three options left: (a) intersects only at and stays on one side of (say, ), (b) intersects only at and changes side there (we may choose for and for ), or (c) lies in .
In case (a) cannot be a geodesic because of strict convexity of . In case (b) define a curve as for and for . By construction of the reflected manifold, the curve is also a geodesic. But now falls in the case (a), which is impossible. Also case (c) is impossible, since a curve lying at a strictly convex subset of the boundary cannot be a geodesic.
We conclude that a geodesic tangential to at a point where is strictly convex cannot exist.
(9) Suppose is a geodesic which meets tangentially at . Consider the case when for . Now construct another geodesic by letting for and for . By the construction of the reflected manifold also is a geodesic. Thus branches at .
Then consider the case when for or the points where intersects accumulate at . Now define as for as above and let for be the unique geodesic in with initial direction . This geodesic exists if is small enough. Again, the curve is a geodesic and branches.
We conclude that any geodesic tangent to at a strictly concave point always has nonunique continuation.
(10) The assumptions imply that is smooth and has smooth and strictly convex boundary. Also for any two boundary points there is a unique geodesic joining them and the geodesic depends smoothly on its endpoints. Thus is simple by definition. ∎
To have unique geodesics on , we want to be (or ). By the above lemma, for this we need that the original metric is or and the reflector is flat in the sense that the second fundamental form vanishes. For higher regularity of we need higher order flatness of the reflector.
For shorthand, we give the following definition so that the various cases of Lemma 13 need not be listed again when it is used.
If is class of functions from to (e.g. or ), we define the reflected class of functions by
4.3. From broken ray transform to geodesic ray transform
The main result we present is Theorem 16. It is a direct generalization of the ideas behind the proofs in Proposition 6, but we state it as a theorem to highlight the generality of the reflection construction. We gave two applications of this theorem in Section 3 to show injectivity of the broken ray transform on a fairly large class of manifolds. Simpler and more concrete examples are given in Section 6. The geodesic ray transform and the broken ray transform were defined in definition 1.
We remind the reader that the set is a manifold with corners. The case of manifolds with smooth boundary requires more work and will be discussed in Section 5.
Let be as in Section 4.1. Let be some classes of functions on and let be the set of tomography. Then:
If determines both and , then determines both and .
If determines for a fixed (known) , then for a fixed (known) .
We only prove part (1); part (2) results by letting . Suppose indeed determines both and , and that is given.
Any stability result for the geodesic ray transform on immediately yields a stability result for the broken ray transform on . Since stability is inherited in such a way, we do not discuss the stability of the broken ray transform on different manifolds.
The above theorem only considers scalar functions . We can similarly define the broken ray transform for a tensor field of any order just like one defines the geodesic ray transform for a tensor field. The theorem holds true for tensor fields as well, and the proof is the same; a tensor field is reflected to instead of simply reflecting all the component functions. (A tensor function can only be recovered up to the natural gauge freedom; see e.g. .) The theorem also remains true if one introduces a weight in the broken ray transform. Replacing real numbers with complex numbers is also a trivial generalization.
The examples in Proposition 6 were concerned with zero attenuation. The attenuated broken ray transform is injective provided the corresponding attenuated ray transform in the plane is injective. The analogue of Helgason’s support theorem holds true with constant attenuation [19, Theorem 4.2]. For more results on attenuated ray transforms in Euclidean spaces, we refer to [23, 5, 25, 8]. Attenuated transforms have also been considered on manifolds (see e.g. [28, 27]), but we set our focus on the nonattenuated setting.
5. Support theorems and manifolds without corners
Theorem 16 above was stated for a manifold with corners. It is appealing to consider the broken ray transform on a manifold with smooth boundary (that is, without corners). To achieve this we take the following two steps: First, using geodesics (broken rays without reflections) with endpoints in one can in favorable situations recover the unknown function in a neighborhood of . (We refer to results of this nature as “support theorems” since they in a way generalize Helgason’s support theorem.) Second, if this neighborhood is nice enough, is a manifold with corners and Theorem 16 is applicable.
We do not formulate this two step procedure as a theorem since the geometry of the set is difficult to control in terms of assumptions on and . A specific example where this idea works is given in Proposition 26(3).
Theorem 20 (Theorem in ).
Let be a manifold with boundary with dimension 3 or greater. If is strictly convex at , then there is a neighborhood of such that the geodesic ray transform is injective in in the following sense: If the integral of a function in vanishes over all geodesics with interior in and endpoints in , then vanishes in .
Theorem 21 ().
Let be a simple Riemannian manifold embedded in a slightly larger manifold and assume that the metric is real analytic. Let be an open set of geodesics in such that that each geodesic can be deformed to a point on the boundary by geodesics in . Let be the set of points lying on the intersection of these geodesics with . If is a function such that the integral of is zero over every geodesic in that has an extension in , then on .
Let be a manifold with boundary and the set of tomography. Suppose any one of the following:
and is open and strictly convex.
The metric is analytic and the manifold simple, the set is open and strictly convex, and can be extended to slightly larger manifold . Any geodesic in with endpoints in can be extended to a geodesic in and can be deformed to a point on by geodesics that do not intersect .
If the broken ray transform of a function vanishes, then in some neighborhood of . Furthermore, in the case (2) vanishes on each geodesic with endpoints in .
(1) Fix any and use Theorem 20 near it. The set in Theorem 20 is constructed so that it can be shrinked to be inside any given neighborhood of in . (For details, see .) Thus in the present case we may choose so that .
If the broken ray transform of a function vanishes, its integral over any geodesic with endpoints in is zero. Therefore vanishes in by Theorem 20. The conclusion holds for all , whence vanishes in a neighborhood of .
(2) Follows from Theorem 21. ∎
Let be the neighborhood of in which vanishes in part (1) of the above theorem. If (in the topology of ) is strictly convex from the side of , the theorem may be used on it again. Such layer stripping might show that vanishes in a relatively large neighborhood (like in part (2)), whose geometry can be controlled more strongly. The global injectivity result of Uhlmann and Vasy  is based on such an argument.
If we know that the support a function is a positive distance away from (the connected set) in case (2) of the above theorem, we do not need to extend . First, we replace by zero outside the convex hull of ; this does not alter the integral of over geodesics with endpoints in , but makes sure that has compact support in . If now , we can use part (2) of the above theorem on the manifold .
In the Euclidean case one can simply use Helgason’s support theorem. This support theorem can be viewed as a special case of part (2) of the above theorem.
6. Examples and counterexamples
Proposition 26 (Examples).
The broken ray transform is injective in the following manifolds (with or without corners):
Consider the quadrant of a sphere
when and the set of tomography . The broken ray transform is injective in the class .
The previous example with in the class .
with the set of tomography for some in the class .
(1) We use Theorem 8. Suppose in the given class has vanishing broken ray transform. Define and let be so small that . Now has vanishing broken ray transform in with the set of tomography .
Let and , where is the intrinsic (Riemannian) metric on , , and is a constant chosen so that whenever . By Theorem 8 the broken ray transform is injective on , whence .
Proposition 27 (Counterexamples).
The broken ray transform fails to be injective on the following kinds of manifolds and sets of tomography :
The manifold is such that the geodesic ray transform is not injective, e.g. a one dimensional manifold. may be anything.
The manifold contains a reflecting tubular part: for , , a bounded set and , the manifold with boundary embeds isometrically to such that is mapped to the complement of .
The manifold contains a reflecting generalized tubular part: for and manifolds with boundary such that the geodesic ray transform on is not injective, the manifold embeds isometrically to such that is mapped to the complement of . We must have but can have .
Parts (1) and (3) hold for the function classes where the geodesic ray transform is non-injective. Part (2) holds for the class .
(1) In the case the broken ray transform is the geodesic ray transform. If the broken ray transform with set of tomography is not injective, it is not injective with any set of tomography either.
(2) Take a function such that but does not vanish identically. Define by . Using the embedding, extend by zero to . We claim that the broken ray transform of vanishes.
It suffices to show that integrates to zero on any (unit speed) broken ray in starting at and ending at . Let be the unit vector normal to the hypersurfaces . Possible reflections at are such that is preserved. Thus the integral over the broken ray becomes (up to a multiplicative constant) the integral of over , which vanishes.
(3) There is a function in such that it integrates to zero over all maximal geodesics in . Define by . For a unit speed broken ray in with both endpoints in the component of the gradient is conserved in reflections and along geodesics. Thus integrates to zero over any such broken rays just like in part (2). ∎
Part (2) of the above proposition is related to the fact that the geodesic ray transform on a one dimensional manifold is not injective; a function on the real line cannot be recovered from its integral. Part (3) naturally generalizes this observation.
As an example of part (3) with we mention .
There is a counterexample to the counterexample given in Proposition 27 which warns us that some counterexamples may fail when attenuation is introduced. We give this as the following proposition. The result could be given for more general manifolds and broken rays, but we only state it here for the simple cylindrical case.
Consider the manifold with boundary. Let be a constant attenuation coefficient. For a function define by .
If , there is a nonzero function such that the ray transform vanishes.
If , there is no nonzero function such that the ray transform vanishes.
(1) Choose a smooth which integrates to zero as in the proof of Proposition 27(2).
(2) Consider geodesics from to of the form , , where . (All geodesics are of this form up to trivial transformations.) We wish to show that if for all , then .
After extending by zero to we find
where is the Laplace transform of . Thus, if for all , we have that for all .
Since is bounded and has compact support, is real analytic on . But vanishes in , so . It follows from the properties of the Laplace transform that . ∎
7. The periodic broken ray transform
In analogue to the broken ray transform introduced in the beginning of Section 1, we now turn to the periodic broken ray transform. In this case the entire boundary is reflecting and the integrals of the unknown function are known over periodic broken rays.
Periodic broken rays are analogous to periodic geodesics on a closed manifold, and this analogy is made precise in the proof of the following two propositions. Guillemin and Kazhdan  reduced spectral rigidity of negatively curved closed Riemannian surfaces to determining a function from its integrals over all periodic geodesics. We therefore expect spectral rigidity of negatively curved surfaces with boundary to be related to determining a function from its integrals over all periodic broken rays. Lengths of periodic broken rays (or periodic billiard orbits) play an important role in spectral geometry (see ). Since linearizing lengths of geodesics with respect to the metric leads to X-ray transforms, the periodic broken ray transform can be expected to have applications in spectral geometry.
In the introductory examples of Section 2 and more generally in Theorem 16 the injectivity of broken ray transforms was reduced to injectivity of certain related geodesic ray transforms via reflections. The geodesics and broken rays considered there joined two points on the boundary or the set of tomography.
The same idea can be carried over to the case of the periodic broken ray transform. We study this idea briefly in this section. The periodic broken ray transform were defined in definition 2.
It is clear that the periodic broken ray transform fails to be injective if there are too few periodic broken rays on the manifold. We consider below specific examples, where the geometry allows for a large number of periodic broken rays.
The periodic broken ray transform is injective for the Riemannian manifold with boundary