BROKEN RAY TRANSFORM FOR TWISTED GEODESICS ON SURFACES WITH A REFLECTING OBSTACLE

. We prove a uniqueness result for the broken ray transform acting on the sums of functions and 1-forms on surfaces in the presence of an external force and a reﬂecting obstacle. We assume that the considered twisted geodesic ﬂows have nonpos-itive curvature. The broken rays are generated from the twisted geodesic ﬂows by the law of reﬂection on the boundary of a suitably convex obstacle. Our work generalizes recent results for the broken geodesic ray transform on surfaces to more general families of curves including the magnetic ﬂows and Gaussian thermostats.


Introduction
This article studies generalizations of the geodesic ray transform to general families of curves.Our main focus will be in broken ray tomography where the trajectories of particles may reflect from the boundary of a reflecting obstacle according to the law of reflection.Furthermore, we consider a situation where the trajectories are influenced by an external force such a magnetic field.Our study limits to the two dimensional case.Our main result, stated later in Theorem 1.1, shows that a function is uniquely determined from the collection of all of its line integrals over the twisted broken rays.We also obtain an analogous result corresponding to vector field tomography with its natural gauge.
Let (M, g) be a smooth Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).We say that a curve γ is a twisted geodesics if it satisfies the λ-geodesic equation where i is the rotation by 90 degrees counterclockwise.The term λ(γ, γ)i γ represents an external force that pushes a particle out from its usual geodesic trajectory where a particle without any influence from an external force would continue its motion.The case when λ(x, v) does not depend on the vertical variable v is called magnetic and the case where λ(x, v) is linear in v is called Gaussian thermostatic, both widely studied in dynamics [DP07a,DP07b,MP22,Sma67,Woj00a,Woj00b].
The geodesic ray transform and closely related Radon transforms are studied by many authors and the area has a long history starting from the early twentieth century [Hel11,IM19,PSU23].More recent advances include the solenoidal injectivity of tensor tomography in two dimensions [Muh77,PSU13,Sha07].In higher dimensions, the geodesic ray transforms are fairly well-understood in negative curvature [GPSU16,PS21,PSU15], and when a manifold has a strictly convex foliation [dHUZ19, PSUZ19, SUV18, UV16].The geodesic ray transform is closely related to the boundary rigidity problem [PU05,SUV21] and the spectral rigidity of closed Riemannian manifolds [GL19, GLP23,PSU14].Other recent considerations include generalizations of many existing results to some classes of open Riemannian manifolds [EG22, GGSU19, GLT20, LRS18] and to the matrix weighted ray transforms [IR20,PSUZ19] as well as their statistical analysis [MNP21a,MNP21b].
The ray transforms for twisted geodesics and general families of curves have been studied recently in [AD18,Zha23] and in the appendix of [UV16] by Hanming Zhou.It is shown in [AD18] that the twisted geodesic ray transform is injective on the simple Finslerian surfaces.For the most recent other studies, we refer to [Zha23].In [JKR23], the injectivity of the matrix attenuated and nonabelian ray transforms for the nontrapping magnetic and thermostatic flows on compact surfaces with boundary are studied.We give a more detailed account of the works that study inverse problems for the magnetic or Gaussian thermostat flows in Sections 2.3.1 and 2.3.2,respectively.
The broken ray transform has been studied extensively.In the case of a strictly convex obstacle, the uniqueness result for the broken ray transform of scalar functions on Riemannian surfaces of nonpositive curvature were obtained in [IS16].This result was later generalized to higher dimensions and tensor fields of any order in [IP22].In the case of rotational (or spherical) symmetry, one may sometimes solve these and related problems using local results and data avoiding the obstacle when the manifold satisfies the Herglotz condition [dHIK22,IM23,Sha97].Broken lens rigidity was studied recently in [dHILS21], and a broken non-Abelian ray transform in Minkowski space in [SA22].Other geometric results include boundary determination from a broken ray transform [Ilm14] and a reflection approach using strong symmetry assumptions [Ilm15b], for example letting to solve the broken ray transforms on flat boxes over closed billiard trajectories [Ilm15a,Hub15].
Numerical reconstruction algorithms and stability for the mentioned problem on the flat boxes would follow directly from [IKR20,Rai20].Artifacts appearing in the inversion of a broken ray transform was studied recently in the flat geometry [Zha20].We refer to the work [Esk04] regarding a possibility to have more than just one reflecting obstacle.Finally, some related results without proofs are stated in the setting of curve families in the Euclidean disk [Muk91].
1.1.Main result.We briefly recall the setting of our work, for further details we point to the later sections.Let (M, g) be a compact oriented smooth Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).Assume that ∂M = E ∪ R where E and R are two relatively open disjoint subsets.Let ν denote the inward unit normal.We define the reflection map ρ : ∂SM → ∂SM by ρ(x, v) = (x, v − 2 v, ν g ν) .
We call a dynamical system (M, g, λ, E) admissible (cf.Definition 3.4 with more details) if (i) the emitter E is strictly λ-convex; (ii) the obstacle R has admissible signed λ-curvature; (iii) the Gaussian λ-curvature of (M, g, λ) is nonpositive; (iv) the broken λ-geodesic flow is nontrapping; (v) there exists a > 0 such that every broken λ-ray γ has at most one reflection at R with | ν, γ | < a.The broken λ-ray transform of f ∈ C 2 (SM) is defined by where φ t (x, v) = (γ x,v (t), γx,v (t)) is the broken λ-geodesic flow, τ x,v is the travel time of φ t (x, v) and π : SM → M is the projection π(x, v) = x.Our main theorem is the following uniqueness result for the sums of functions and 1-forms.
Our proof is based on the ideas introduced in [IP22, IS16] and the Pestov identity for the twisted geodesic flows in [AD18,DP07a].The proof could be split into three main parts, each of them having their own technical challenges, resolved in our work for a general λ ∈ C ∞ (SM).
(i) One has to analyze a generalized Pestov identity with boundary terms and decompose the boundary terms into the even and odd parts with respect to the reflection; (ii) One has to reduce the problem into a related transport problem for the broken λ-geodesics and observe if the Pestov identity and the analysis of the first step applies to the solutions of the transport equation, then the problem can be solved; (iii) One has to show sufficient regularity for the solutions of the broken transport equation.This can be done by analyzing carefully the behavior of the broken λ-rays, broken λ-Jacobi fields at the reflection points, and utilizing a time-reversibility property for the pair of flows with respect to λ(x, v) and −λ(x, −v).

Preliminaries
In this section, we introduce our notation and concisely present many primary definitions employed throughout the article.We closely follow the notation used in the recently published book by Paternain, Salo and Uhlmann [PSU23].We also refer to [Lee18] for the basics of Riemannian geometry.
2.1.Basic notation.Throughout the article, we denote by (M, g) a complete oriented smooth Riemannian manifold with or without boundary.We always assume that M is a surface, i.e. dim(M) = 2.We denote the Levi-Civita connection or covariant derivative of g by ∇, and the determinant of g by |g|.When the covariant derivative is restricted to a smooth curve γ, we simply write D t = ∇ γ .We sometimes emphasize the base point x ∈ M in the notation of the metric as g x and other operators but this is often omitted.We denote the volume form of (M, g) by dV 2 := |g| 2 dx 1 ∧ dx 2 , expressed in any local positively oriented coordinates.For any vector v ∈ T x M, let v ⊥ ∈ T x M denote the unique vector obtained by rotating v counterclockwise by 90 • .This vector satisfies: and forms a positively oriented basis of T x M with v, when v = 0. Note that, often, we may also write v ⊥ as iv.Given a vector v ∈ T x M and a positive orientation on M, we denote by v ⊥ = −v ⊥ the clockwise rotation.We denote by K the Gaussian curvature of (M, g).
The signed curvature of the boundary ∂M is defined as where δ(t) represents an oriented unit-speed curve that parametrizes the boundary ∂M and ν is the inward unit normal along δ(t).For a comprehensive explanation, please refer to [Lee18, Chapter 9, p. 273].Furthermore, we define the second fundamental form of the boundary ∂M as follows: We say that M has a strictly convex boundary if ∂M is strictly convex for any x ∈ ∂M.We say that ∂M is strictly concave at x ∈ ∂M if II x is negative definite, i.e.II x (v, v) < 0 for any v ∈ T x ∂M \ {0}.The relation between the signed curvature and the second fundamental form is given by κ = II( δ, δ) (cf.[Lee18, Chapters 8 and 9]).If the boundary is strictly convex, it implies that the signed curvature is positive, whereas if it is strictly concave, the signed curvature is negative.
2.2.Analysis on the sphere bundle.In this article, we often assume that M is a Riemannian manifold with smooth boundary ∂M.We denote by SM the unit sphere bundle of M, i.e.
The boundary of SM is given by We define the influx and outflux boundaries of SM as the following sets The glancing region is defined as ∂ 0 SM := ∂ + SM ∩ ∂ − SM = S(∂M).We denote by dS x the volume form of (S x M, g x ) for any x ∈ M. The sphere bundle SM is naturally associated with the measure Let us denote by X the geodesic vector field, V the vertical vector field and the orthogonal vector field X ⊥ := [X, V ] (cf.[PSU23, Chapter 3.5] for more details).The following structure equations hold where K is the Gaussian curvature of (M, g) [PSU23, Lemma 3.5.5].The sphere bundle SM is equipped with the unique Riemannian metric G such that {X, −X ⊥ , V } forms a positively oriented orthonormal frame.The metric G is called the Sasaki metric, and it holds that dV G = dΣ 3 [PSU23, Lemma 3.5.11].We denote by dV 1 the volume form of (∂M, g).This leads to the definition of the volume form on ∂SM.We note that dΣ 2 = dV ∂SM where on the right hand side the volume form is induced by the Sasaki metric.We define the following L 2 inner products (u, w) SM := We next recall simple integration by parts formulas.For any u, w ∈ C 1 (SM), the following formulas hold [PSU23, Proposition 3.5.12]: (2.2) Finally, we recall the vertical Fourier decomposition.We define the following spaces of eigenvectors of V and the series converges in C ∞ (SM).We next discuss some important boundary operators following [PSU23, Lemma 4.5.4].We define the tangential vector field T on ∂SM by setting that where µ(x, v) := ν(x), v g .The Pestov identity with boundary terms is due to Ilmavirta and Salo [IS16, Lemma 8]; see also [PSU23, Proposition 4.5.5] and the higher dimensional generalization in [IP22,Lemma 8].For any u ∈ C 2 (SM) it holds that ) is a compact Riemannian surface with smooth boundary.We generalize (2.5) in Proposition 2.1 to the case where X is replaced by the generator of a twisted geodesic flow.
2.3.Twisted geodesic flows.We first recall the concept of λ-geodesic flows from the lectures of Merry and Paternain [MP11,Chapter 7].Let (M, g) be a complete oriented Riemannian surface (with or without boundary) and λ ∈ C ∞ (SM) be a smooth real valued function.We say that a curve γ : When λ ≡ 0, then the λ-geodesics are the usual geodesics of (M, g).One may think that the function λ twists the usual geodesics in order to model trajectories of particles when moving in the presence of external forces.When λ is a smooth function on M or 1-form, then the twisted geodesics correspond to the magnetic and thermostatic geodesics, respectively.For other advances in the context of inverse problems for λ-geodesics, we refer to [AD18,Zha23].In particular, the class of λ-geodesics is large and can be characterized with only three natural properties of a curve family, for details see [AD18, Theorem 1.4].Let γ x,v be the unique λ-geodesic with the initial condition (γ x,v (t), γx,v (t)) = (x, v) and solving (2.6).As in [MP11, Exercise 7.2], one may define the λ-geodesic flow by setting that The infinitesimal generator of the λ-geodesic flow F is given by Using (2.1), one may derive the following commutator formulas [DP07a, p. 537]: We use the following notations Notice that the formal L 2 adjoint of the operator P is P * = (F + V (λ))V = P .Let us define the λ-curvature of (M, g, λ) as a map in C ∞ (SM) by the formula We can now recall the following generalized Pestov identity by Dairbekov and Paternain for closed surfaces M [DP07a, Theorem 3.3]: (2.8) . This also holds on compact surfaces with boundary if one additionally assumes that u| ∂SM = 0. We generalize (2.5) and (2.8) to the surfaces with smooth boundary without making the simplifying assumption that u| ∂SM = 0 (cf.Proposition 2.1).We also remark that similar Pestov identities, but in a slightly less explicit form, were derived by Assylbekov and Dairbekov for the λ-geodesic flows on Finslerian surfaces [AD18, Theorem 2.3].In turn, the generalized Pestov identity with boundary terms is used to study generalized broken ray transforms.We define the signed λ-curvature as Additionally, we introduce a related term: (2.10) which will appear in the Pestov identity (2.13).We remark that κ λ and η λ depend only on the values of λ on ∂SM.For further details, see Subsection A.3.
2.3.1.Magnetic flows.We refer to the articles of Arnold [Arn61] and Anosov-Sinai [AS67] as first mathematical studies of magnetic flows.We will mainly follow the notation used by Ainsworth in the series of works [Ain13, Ain15, AA15], considering the integral geometry of magnetic flows.We further note the following works related to different inverse problems for magnetic flows [AZ15, DPSU07, Her12, HV11, Jol07a, Jol07b], including the boundary, lens and scattering rigidity problems.
Let Ω be a closed 2-form on M modeling a magnetic field.The Lorentz force Y : T M → T M associated with the magnetic field Ω is the unique bundle map such that We say that γ is a magnetic geodesic if it satisfies the magnetic geodesic equation (2.11) Notice now that since M is orientable, there exists a unique function λ : M → R such that Ω = λdV 2 .We may define λ = λ • π.Now it holds that γ solves (2.11) if and only if it solves (2.6).We may define the magnetic flow simply as the corresponding λ-geodesic flow with the fixed energy level 1/2 corresponding to the unit speed curves.One may also view the magnetic flow as the Hamiltonian flow of H(x, v) = 1 2 |v| 2 g under the symplectic form where ω 0 is the symplectic structure of T M generated by the metric pullback of the canonical symplectic form on T * M. The magnetic geodesics are known to be constant speed and different energy levels lead to different curves.We also remark that the magnetic flow is time-reversible if and only if Ω ≡ 0. Therefore, γ x,v (−t) is not a magnetic geodesic of (M, g, Ω).However, we have that the magnetic field with flipped sign −Ω reverses the orientation of geodesics, i.e. γ −Ω x,−v (t) = γ Ω x,v (−t).One may check that −Ω is the dual of Ω in the sense of Section 2.4.2.3.2.Thermostatic flows.The Gaussian thermostats concept was proposed by Hoover for the analysis of dynamical systems in mechanics [Hoo86], but it appears also earlier, for example, in the work of Smale [Sma67].The inverse problem for Gaussian thermostats has been more recently explored in the works of Dairbekov and Paternain [DP07a,DP07b].Other contributions have also been made by Assylbekov and Zhou [AZ17], Assylbekov and Dairbekov [AD18], and Assylbekov and Rea [AR21].In addition, the dynamical and geometrical properties of Gaussian thermostats have been extensively studied, as demonstrated in the contributions by Wojtkowski [Woj00a,Woj00b], Paternain [Pat07], Assylbekov and Dairbekov [AD14], and Mettler and Paternain [MP19,MP20].Gaussian thermostat also arises in Weyl geometry, see for instance [PW08].
Consider a smooth vector field E on M, representing an external field.A thermostatic geodesic satisfies the equation The flow φ t = (γ(t), γ(t)) is called as thermostatic flow.It is worth mentioning that thermostatic geodesics are time-reversible, meaning that φ t (x, −v) = r (φ −t (x, v)), where r : SM → SM is the reversion map defined by r(x, v) = (x, −v).When E = 0, then the thermostatic geodesics are the usual geodesic.Given a 1-form λ defined by λ(x, v) := E(x), iv , the equation (2.12) can be rewritten as the corresponding λ-geodesic equation 2.4.Dual λ-geodesic flow.It is well known that the usual geodesics are time-reversible.The magnetic geodesics are not time reversible unless Ω ≡ 0 (cf.[DPSU07, p. 537]).In [BK22, p. 100], it is mentioned that the magnetic geodesics corresponding to λ and −λ are one-to-one correspondence through time reversal.This means that a curve t → γ x,v (t) is a magnetic λ-geodesic if and only if the curve t → γ x,v (−t) is a magnetic (−λ)-geodesic.However, the thermostatic λ-geodesic flow is time-reversible, see for instance in [Pat07, p. 88]).Therefore, the λ-geodesic flow in the case of λ ∈ C ∞ (SM) is in general not time-reversible.
Next, we will define the corresponding dual λ-geodesic flow.We can define the timereversed dynamical system related to λ by setting

It now follows that γ −
x,−v (t) = γ x,v (−t) where γ − x,−v is a unique λ − -geodesic with initial data (x, −v).We call λ − as the dual of λ.This time-reversibility property can be checked by substituting γ − x,−v (t) := γ x,v (−t) to the λ − -geodesic equation and using the fact that γ x,v solves the λ-geodesic equation.In fact, and in local coordinates where we write simply γ = γ x,v and γ − = γ − x,−v .Since ∇ γ γ = λ(γ, γ)i γ holds for all times s in the maximal domain of the λ-geodesic γ x,v , we can conclude that γ − is a λ − -geodesic.On the other hand, γ− (0) = −v and γ − (0) = x, which justifies writing γ − = γ − x,−v .The generator of the dual λ-geodesic flow φ − t is simply given by F − := X + λ − V .Additionally, it is worth noting that (λ − ) − = λ.We will use the dual flow to establish regularity results for the solutions of a broken transport equation in Sections 3.2 and 3.4.For the sake of completeness, we will discuss in Section A.4 the curvature and signed curvature of λ and λ − .Proposition 2.1 (Generalized Pestov identity).Let (M, g) be a compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).If u ∈ C 2 (SM), then we have We say that a vector field J along a λ-geodesic γ is a λ-Jacobi field if it is a variation field of γ through λ-geodesics, i.e.J(t) = ∂ s γ s (t)| s=0 where γ s (t) is a smooth oneparameter family of λ-geodesics with γ = γ 0 (for further details we refer to Section A.2). We will next state a useful estimate on the growth rate of λ-Jacobi fields.On compact Riemannian surfaces the norm of a Jacobi field and its covariant derivative grow at most exponentially (see e.g.[IS16, Lemma 10]).Such inequalities are useful in studying the stability of geodesics and their relation to the curvature of the manifold.In the following lemma, we establish a similar result for the Jacobi equation associated with λ-geodesics.A detailed proof is given at the end of Section A.2. Lemma 2.2.Let (M, g) be a compact Riemannian surface with or without boundary and λ ∈ C ∞ (SM).Let J be a λ-Jacobi field defined along a unit speed λ-geodesic γ : [a, b] → M. Then the following growth estimate holds for all t ∈ [a, b] where C is a uniform constant depending only on M, g and λ.
2.6.Even and odd decomposition with respect to the reflection map.Given a Riemannian surface (M, g) with smooth boundary.We define the reflection map ρ : ) .We denote by ρ * the pull back of ρ.The even and odd parts of u : SM → R with respect to the reflection map ρ are defined by the formula We will next state a simple lemma related to the reflection and rotation maps.We omit presenting a proof.
Lemma 2.3.Let (M, g) be a Riemannian surface with smooth boundary.The reflection and the rotation maps satisfy The boundary operators κ λ and ν λ , as defined in (2.9) and (2.10), satisfy the following simple identities.These formulas will later on allow us to simplify the boundary terms in (2.1).
Lemma 2.4.Let (M, g) be Riemannian surface with smooth boundary.Then κ λ and η λ satisfy the following properties A proof of Lemma 2.4 is given in Appendix A.5.

Transport equation for the broken λ-geodesics
3.1.Broken ray transforms and the transport equation.Let (M, g) be an orientable, compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).We assume that ∂M can be decomposed to the union of two disjoint relatively open disjoint subsets E and R. In particular, one could think of M as M = M \ O, where M is a larger manifold containing M and O being an open obstacle.In this case, E = ∂ M and R = ∂O are the outer and inner boundaries of M. We call E the emitter and R the reflector of M. In Section 2.4, we denote by φ t the usual λ-geodesic flow and by φ − t its dual flow.By abuse of notation, we now continue to write the same notation φ t for the broken λ-geodesic flow and φ − t to its dual flow.For any (x, v) ∈ SM, we define the forward and dual travel times by x,v for most of (x, v) ∈ SM since the twisted geodesic flows are not reversible.On the other hand, the maximal domain of Let us denote by ρ : π −1 R → π −1 R the reflection map and define by the law of reflection Note that ρ is an involution in the sense that ρ • ρ = Id.Here and subsequently, γ(t 0 −) stands for the left-hand limit of γ at some point t 0 and γ(t 0 +) denotes the right-hand limit of γ at t 0 , which are defined by γ(t 0 −) = lim t→t − 0 γ(t) and γ(t 0 +) = lim t→t + 0 γ(t).Similarly, γ(t 0 −) = lim t→t − 0 γ(t) and γ(t 0 +) = lim t→t + 0 γ(t).Definition 3.1.Let (M, g) be a Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).Assume that ∂M = E ∪ R where E and R are two relatively open disjoint subsets.A curve γ on M is called a broken λ-ray if γ is a λ-geodesic in int(M) and reflects on R according to the law of reflection whenever there is a reflection at γ(t 0 ) ∈ R.
The broken λ-ray transform of f ∈ C 2 (SM) is defined by ) is the broken λ-geodesic flow.We next move towards studying the injectivity of the broken λ-ray transform, that is, is it possible to determine an unknown f from the knowledge of its integrals (3.1) over maximal broken λ-rays?
To proceed, we first reduce the integral equation (3.1) to a certain partial differential equation.Given any f ∈ C 2 (SM), define Notice that the exit time τ x,v is smooth near (x, v) ∈ SM whenever the broken ray γ x,v reflects and exits transversely.A simple application of the fundamental theorem of calculus together with the regularity properties of exit time, we deduce from (3.2) that u satisfies the transport equation We need to make some assumption on the geometry of M and its reflecting boundary parts to make the injectivity problem more approachable.
where (κ λ ) o and (η λ ) o are odd parts of the functions κ λ and η λ respectively.This assumption holds, for example, when λ is independent of the direction v on π −1 R. Definition 3.3.Let (M, g) be a Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).We say that a relatively open subset R ⊂ ∂M of the boundary has an admissible λ-curvature if the following inequality holds: λ is only a function of the base point on R, and R is strictly λ-concave at any x ∈ R, then R has an admissible λ-curvature.
From Corollary A.14 and Remark A.15, we have i.e., an obstacle R has admissible λ-curvature if and only if R has admissible λ − -curvature.If V (λ e )| R = 0, then the condition that the obstacle R has admissible λ-curvature is equivalent to the strict λ-concavity of R.
In [IS16, Theorem 1], it was proved that one can recover an unknown function f from the knowledge of its geodesic broken ray transform If .They assumed that the surface is nontrapping, having nonpositive Gaussian curvature, the reflecting part is strictly concave, and the broken rays allow at most one reflection with | γ, ν | < a. See also [IP22, Definition 1] for similar assumptions used to study the broken ray transforms in three and higher dimensions.We now define a similar class of admissible Riemannian surfaces with broken λ-geodesic flows.Definition 3.4.Let (M, g) be a compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).Assume that ∂M = E ∪ R where E and R are two relatively open disjoint subsets.We call a dynamical system (M, g, λ, E) admissible if (1) the emitter E is strictly λ-convex in the sense of Definition A.11; (2) the obstacle R has admissible λ-curvature in the sense of Definition 3.3; (3) the λ-curvature K λ of (M, g, λ) is nonpositive; (4) the broken λ-geodesic flow is nontrapping: there exists L > 0 such that τ x,v , τ − x,v < L for any (x, v) ∈ SM; (5) there exists a > 0 such that every broken λ-ray γ has at most one reflection at R with | ν, γ | < a.

Dual transport equation.
Let us define the dual transport equation of (3.3) as where f (x, v) := f (x, −v) and I − denote the broken ray transform related to the λ −geodesic flow.For nontrapping broken λ-geodesic flows, we define the scattering relation Lemma 3.5.Let (M, g) be a compact nontrapping Riemannian surface with smooth boundary and where u is the unique solution of the transport equation (3.3) and u − is the unique solution of the dual transport equation (3.4).
Proof.Let us consider (x, v) ∈ SM.Now the union of curves γ x,v and γ − x,−v form a maximal broken λ-geodesic (cf.Section 2.4) with its endpoints on ∂M.By assumption, By the definition u − is the unique solution of (3.4).Note that where α is the scattering relation and r is the reversion map.This implies that u The formulas (3.5) and (3.6) imply that which completes the proof.

3.3.
Jacobi fields for broken λ-geodesics.In this subsection, we give a brief exposition of Jacobi fields along broken λ-geodesic flows.In Lemma A.8, we show that Jacobi fields along a λ-geodesic flow can be realized as the variation field of a unit speed λgeodesic variation of γ.Here, we will generalize similar properties of Jacobi fields to the case of broken λ-geodesic flows.The crucial point is to understand how the Jacobi fields behave at reflection points.The Jacobi fields along broken rays have been studied in the case of λ = 0 in [IS16, Section 5] and we follow some of the techniques from this article.Let x 0 ∈ ∂M and ν be the inward unit normal to it.We define a map Φ ζ : for any vector ζ ∈ T x 0 M that is not orthogonal to ν, where the map ϕ ζ : Since ϕ ζ ξ ⊥ ν, the derivative ∇ ϕ ζ ξ ν is well defined.To analyze Jacobi fields at reflection points, we first make an assumption on λ on the reflected part of the boundary.In particular, we require a condition on λ such that where β ∈ L ∞ (π −1 R).We refer to this condition as the β-reflection condition on λ.
Taking (x, v) → ρ(x, v), we have which gives the condition on β that 1 = β(x, v)β(ρ(x, v)).For example, if we consider , then we have a β given as follows However, at the end of this section, we are able to establish suitable growth estimates for λ-Jacobi fields without this additional β-reflection condition.The benefit of the βreflection condition is that it allows one to write down a clean reflection condition for the broken Jacobi fields.Definition 3.7.Let (M, g) be a compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).Assume that ∂M = E ∪ R where E and R are two relatively open disjoint subsets and λ • ρ| π −1 R = βλ| π −1 R for some β ∈ L ∞ (π −1 R).Let γ be a broken λ-ray without tangential reflections.Then a vector field J along γ is a Jacobi field along γ if (i) J is a λ-Jacobi field along the segments of λ-geodesic γ in int(SM) in the usual sense (cf.Section A.2); (ii) if γ has a reflection at γ (t 0 ) ∈ R, then the left and right limits of J at t 0 are related via J (t 0 +) = ρJ (t 0 −) , and where R and hence the identity (3.7) is equivalent to Remark 3.8.In the case of usual geodesics it holds that D t I = D t γ = 0.
In [IS16,Lemma 12], it has been pointed out that if none of the broken geodesic rays γ s have tangential reflections, then J(t) = ∂ s γ s (t)| s=0 is a Jacobi field along γ 0 , where γ s are the variations of γ 0 .Conversely, any Jacobi field can be understood as a variation of the broken geodesic γ 0 .We can now state an analogue of [IS16, Lemma 12] in the case of broken λ-rays, where λ ∈ C ∞ (SM).Lemma 3.9.Let (M, g) be a compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).Assume that ∂M = E ∪ R where E and R are two relatively open disjoint subsets.Let γ s : [0, L] → M be the broken λ-ray starting at (x s , v s ) where the parametrization (−ε, ε) ∋ s → (x s , v s ) ∈ int SM is a smooth map.If the broken λ-rays γ s do not have tangential reflections and | s=0 is a Jacobi field along the broken λ-ray γ 0 .
Proof.By Lemma A.1, it follows that J satisfies the λ-Jacobi equation on each λ-geodesic segment between the reflection points.Therefore, it is sufficient to prove that J satisfies (3.7) at the reflection points.Let γ := γ 0 be a broken λ-ray that does not have a tangential reflection at the reflection point t = t 0 .After possibly shrinking the domain of definition, we may assume that t = t 0 is the only reflection point of γ 0 and each of the broken rays γ s have at most one reflection.
We begin by proving the lemma for a special family of curves corresponding to the tangential Jacobi fields.Let us consider a family of broken λ-geodesics γ s (t) = γ(t + s) with the starting point and velocity given by (γ s (0), γs (0)) = (x s , v s ).We denote I by the variation field corresponding to γ s .Notice that is a λ-Jacobi field except at the reflection point, see Lemma A.1.
We now analyze the behavior of I at reflection point t 0 .By the definition of the broken λ-ray, we have I (t 0 +) = ρI (t 0 −).Since γ s satisfies the λ-geodesic equation, we have D t I(t) = λ(γ(t), γ(t))i γ(t) outside the reflection point, and this leads to Applying Lemma 2.3 to the above identity, we obtain Hence the vector field I satisfies (3.7).Now it remains to prove the lemma in the case of general variations of broken λ-rays.Let J denote the vector field along γ, as given in the statement of the lemma.In view of the proof of [IS16, Lemma 12], we define another vector field J along γ by setting Since γ(t) does not have a tangential reflection, we see that Note that I(t) and J(t) satisfy the λ-Jacobi equation except at the point of reflection, see Lemma A.1.Therefore, by the linearity of the λ-Jacobi equation, the vector field J must satisfy the λ-Jacobi equation except at the reflection points.Similar to Lemma A.8, there exists a corresponding family of broken λ-geodesic variations associated with J, say One can make a change of order s 2 to (x s , v s ) without changing J. Since J(t 0 −) ⊥ ν and γ 0 arrives to R transversely at t 0 , we can introduce a second order change to (x s , v s ) such that γ s (t 0 ) ∈ R for all s ∈ (−ǫ ′ , ǫ ′ ) after choosing a sufficiently small ǫ ′ ∈ (0, ǫ) (cf.[IS16,p. 396]).This variation after the second order change is explicitly given by s → φ τs (x s , v s ) where τ s is the unique element in for a sufficiently small δ > 0 depending upon the choice of ǫ ′ .One can check that τ 0 = 0 and ∂ s τ s | s=0 = 0 using that γ 0 (t 0 ) ∈ R and J(t 0 −) ⊥ ν.Since the family of curves s → γ s (t 0 − δ ′ ), δ ′ > 0 sufficiently small, arrives in the limit δ ′ → 0 tangentially to R at s = 0, it follows that ∂ s τ s | s=0 = 0. We may assume without loss of generality that all γ s have their unique reflection at t 0 .This implies that γ s (t 0 ) = γ xs,vs (t 0 ) is smooth at s = 0 and By (A.6), we have Let us denote by y s = γ s (t 0 ), u s = γs (t 0 −) and ν(y s ) = ν s .Thus we have Using the fact that J(t 0 −) ⊥ ν, we obtain We also have 2 ( (3.11) Taking s = 0 in (3.10), and combining (3.8) with (3.9) and (3.11), we get (3.12) Finally, we have and From (3.13), the above identity becomes which is our desired conclusion.
We now simplify the map Φ γ(t 0 −) J (t 0 −).To do so, we first compute By properties of covariant derivative along curves (cf.[Lee18, Theorem 4.24]), we have Note that Φ γ(t 0 −) is a linear map.Due to the compactness and smoothness of R, similar to [IS16, Remark 13], we establish that where the field of linear maps A is uniformly bounded over the set π −1 R.This follows using compactness and smoothness of R since A at (x, v) ∈ π −1 R depends only ν, its first derivatives and the choice of a unit vector (x, v) = (γ(t 0 −), γ(t 0 −)).This implies that the Jacobi field tends to infinity as the reflection becomes more tangential.In particular, we have obtained the following auxiliary result.
Remark 3.10.The reflection condition on a Jacobi field along a broken λ-ray is given by (3.15) Lemma 3.11.Let (M, g) be a compact Riemannian surface with smooth boundary.Under the assumptions of Lemma 3.9, a Jacobi field J along a broken λ-ray satisfies at every reflection point t 0 , where C is a constant depending on M, g and λ.
Proof.From (3.15) and (3.14), we have which implies Remark 3.12.Consider a family of broken λ-rays on M satisfying | γ, ν | ≥ a at each of the reflection points.Due to the compactness of M and the requirement for traversality | γ, ν | ≥ a, we can assert the existence of a positive real number l that bounds from below by the distance between any two consecutive reflection points for any broken λ-ray in the set.This provides us with a lower bound on the number of reflections.If we denote by N(t) the number of reflections of γ in the time interval (0, t), then from the preceding discussion it follows that (N(t) − 1)l ≤ t, implying N(t) ≤ 1 + t l .Remark 3.13.Consider a unit-speed C 1 curve on the manifold SM, given by the mapping s ∈ (−ε, ε) → (x s , v s ).Let γ xs,vs (t) denote a λ-geodesic such that its initial conditions are (γ xs,vs (0), γxs,vs (0)) = (x s , v s ).Assume that γ := γ x 0 ,v 0 .Now, let us examine the Jacobi field J s (t) = ∂ s γ xs,vs (t).In this context, we have Additionally, consider another Jacobi field I(t) = γx,v (t).In this case, we obtain Without any assumption λ • ρ = λβ on π −1 R, we can still control |J(t)| 2 + |D t J(t)| 2 as follows: Corollary 3.14.Let (M, g) be a compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM), and fix a number a ∈ (0, 1].Let γ : [0, τ ] → M be a broken λ-ray on M such that | γ, ν | ≥ a at every reflection point.Then for any variation field J along γ, we have for all t ∈ [0, τ ], where A and B are constants that depend only on M, g, λ and a. Proof.Let us assume that there is a reflection at each t = t k , where k ∈ {0, • • • , N}.Then for any t ∈ (0, t 0 ), using Lemma 2.2, we have and also by considering Jacobi field I(t) = γ(t), we have . This proves the case for N = 0.
Similar to the proof of Lemma 3.9, we have J such that We have Also we have | J(t 0 +)| = | J(t 0 −)| and from (3.12), we obtain We assume that for t ∈ (0, t k−1 ) i.e. before the k-th reflection, we have We now prove the estimate for any t ∈ (0, t k ), that is before the (k + 1)-th reflection.Using the assumption that | γ, ν | ≥ a where a ∈ (0, 1] and Lemma 3.11, we can deduce that at each reflection point, we have Combining (3.16) with (3.17) and Lemma 3.11, we have .
From this analysis one can see that there is constant A = e Ct 2C 1 a , such that at each reflection point |J(t)| 2 + |D t J(t)| 2 increases by the factor A. Now consider the interval (0, t) where any broken ray has less than 1 + t/l number of reflections by Remark 3.12.We may conclude the estimate , where B = log A l + C. Note that here, A and B are constants that depend only on M, g, λ, and a, which proves the lemma.
3.4.Regularity of solutions to the transport equation.In [IP22, Lemma 3], a regularity result is proven for the primitive function corresponding to f ∈ C 2 (SM), which is simply the solution to the transport equation.This result shows that the primitive function is twice continuously differentiable in the interior of SM and Lipschitz continuous in SM.In this subsection, we extend this result to the solution of the transport equation associated with broken λ-rays.We first aim at establishing the regularity result for both forward and backward exit times.In [PSU23, Lemma 3.2.3], the regularity of exit time has been demonstrated in the context of usual geodesics.Following a similar approach, we extend this result to the case of the λ-geodesic.
Proof.We establish the result for the forward exit time only as the proof is identical for τ − .Let (N, g) be a closed extension of (M, g) (cf.[PSU23, Lemma 3.1.8]).Define ρ as a boundary defining function on M. Consider a λ-geodesic γ with initial conditions γ(0) = x and γ(0) = v.We now analyze the function ρ(γ x,v (t)) for (x, v) ∈ SM \ ∂ 0 SM, where ρ is a boundary defining function on M.
First, we show that τ is continuous on SM.
Therefore, from (3.18) and (3.19), we can find ǫ > 0 such that |τ (x, v) − τ (x 0 , v 0 )| < ǫ for (x, v) in the neighborhood of (x 0 , v 0 ).A similar argument applies for (x 0 , v 0 ) ∈ ∂SM, in which case, τ (x 0 , v 0 ) = 0. Next, we show that τ is smooth in SM \∂ 0 SM.Similar to the proof of [PSU23, Lemma 3.2.3],we have h : Note that γ x,v (τ x,v ) ∈ ∂M.This implies that the tangent vector γx,v (τ x,v ) must lie in ∂ − SM; otherwise if γx,v (τ x,v ) were not in ∂ − SM, the geodesic γ x,v could be extended beyond the point γ x,v (τ x,v ), contradicting the fact that γ x,v (τ x,v ) is the final point on M. By strict λ-convexity, one must have γx,v (τ it follows from the definition of boundary defining function that h(x, v, τ x,v ) = 0 and h is smooth.Finally, by the implicit function theorem, we conclude that τ is smooth in SM \ ∂ 0 SM.
Remark 3.16.Similar to the case of broken rays (cf.[IS16, p. 399]), using Lemma 3.15, τ and τ − are smooth near any point (x, v) such that the broken λ x,v -ray reflects and exits transversely.
Note that the function h (same as defined in Lemma 3.15), satisfies the property that h(x s , v s , τ (x s , v s )) = 0 and ∂h ∂t (x s , v s , τ xs,vs ) = 0.By the implicit function theorem, we have (3.20)Definition 3.18.Let (M, g) be a Riemannian surface and λ ∈ C ∞ (SM).Let us denote the interior λ-scattering relation by α : SM → ∂SM.Given any point and direction (x, v) ∈ SM, we map it via α to the first intersection point and direction of the λ-geodesic γ x,v with the boundary ∂M (i.e.either E or R).
Remark 3.19.Note that in the case a λ-geodesic γ x,v hits ∂M non-tangentially, then in a neighborhood of (x, v) the map α is smooth.
Remark 3.20.When a broken λ-ray γ x,v hits non-tangentially to R (possibly multiple times) and reach a point on E transversally (by strict convexity), then there is a smooth dependence for the end point on E of the broken λ-ray on its initial data (x, v).This follows from the smoothness of α and ρ map.
It is clear that u solves (3.3).Hence, we split the proof into two cases for the regularity of u.
where u − is the solution to the dual transport equation.Hence it suffices to show that either u or u − is C 2 at (x, v) or (x, −v) respectively.Without loss of generality, we may assume γ x,v (t) has no tangential reflections.Now, for some N = 0, 1, 2, . . ., we have Since the broken λ-ray hits R transversely, τ k (x, v) are smooth in some neighborhood of (x, v), λ-geodesic flow is smooth and f ∈ C 2 (SM), we have that all u k (x, v) are C 2 functions in some neighborhood of the point (x, v).As each u k is a C 2 function at the point (x, v), it follows from (3.21) that the function u is also C 2 at (x, v).
Proof of u ∈ Lip(SM).If we show that first order derivatives u are uniformly bounded in int(SM), then this implies u is Lipschitz.To show this, similar to [IP22, p. 1283], we consider a C 1 unit speed curve (−ε, ε) ∋ s → (x s , v s ) ∈ int SM with (x 0 , v 0 ) = (x, v).Using Lemma 3.5 again, we can assume without loss of generality that γ x,v have no tangential reflections and | γ, ν | ≥ a.Now Let us start by examining the integral term where J := ∂ s γ xs,vs is a broken λ-Jacobi field.Since γ x,v contains no reflections with | γ, ν | < a, it follows from Lemma 3.14 and Remark 3.13 that there exists a uniform We now focus on the boundary term.By taking s = 0 in (3.20), we have .
From the proof of Lemma 3.15, for any (x, v) ∈ int SM, γx,v (τ x,v ) , ν < 0. From the expression (3.23), we need to consider two cases: We choose b such that whenever | γx,v (τ x,v ), ν | < b, then the corresponding broken geodesic has no reflections for any (x, v) ∈ int SM.A choice of a very small parameter b > 0 splits int(SM) into two sets corresponding to the broken λ-rays to short ones which are almost tangential to E and all other broken rays since E strictly λ-convex.
By the choice of b, we have that γ x,v never reaches R and corresponds to a short λ-geodesic almost tangential to ∂M.Since the broken ray transform vanishes, we have f (y, w) = 0 holds for y ∈ E and w ∈ S y E. Write (y s , w s ) = (γ xs,vs (τ xs,vs ) , γxs,vs (τ xs,vs )).Since f is Lipschitz, for any w ∈ S ys E, we have Let us express w s in terms of ν and w where we choose the orientation so that w ∈ S ys E and w, w s ≥ 0. Now
4. Uniqueness for scalar functions and 1-forms 4.1.Revisiting the boundary terms in the Pestov identity.The primitive function corresponding to f , denoted by u := u f , is defined as where φ t is the broken λ-geodesic flow.
In the following lemma, we provide a simplified form of the boundary term ∇ T,λ appearing in the Pestov identity (2.13) in terms of the odd and even components of u and the magnetic signed curvature.In [IS16, Lemma 9], a similar identity has been proved for the broken geodesic flow (i.e. when λ = 0).In particular, they showed that where ∇ T = ∇ T,0 , κ := − D t T, ν g is the signed curvature of ∂M, u e and u o are the even and odd components of u| ∂SM with respect to the reflection ρ and u is a primitive function.We now aim to prove the following generalization of [IS16, Lemma 9] to the case of broken λ-geodesic flows.λ ∈ C ∞ (SM).If u ∈ C 2 (SM), then (4.2) Proof.We denote ∇ T,λ u = ∇ T u + Lu where Lu := − v ⊥ , ν λV u − v, ν (V (λ))V u.Note that from [PSU23, p. 119], we have µ(x, v) = v, ν(x) and V (µ)(x, v) = v ⊥ , ν .This implies Lu = −V (λµ)V u.From [IS16, p. 391], we have (ρ * V u) = −V (ρ * u).We compute From (4.3) and (4.4), we have Since ρ is an isometry on S x for each x ∈ ∂M (cf.Remark 4.2), we obtain Combining with (4.1), we have Remark 4.2.Notice that ρ is an isometry on S x for each x ∈ ∂M.The even and odd parts of u are denoted by u e and u 0 respectively with respect to the isometry ρ.Similarly, v e and v o stands for the even and odd parts of v respectively with respect to the isometry ρ.Then where we used the fact that ρ is an isometry.Similarly, where Proof.We have By the assumption on u, we have u e = u on ∂SM and u o = 0 on ∂SM.Combining (4.5) and (4.6), we obtain From (4.2) and Lemma 2.4, we have Proposition 2.1 and Corollary 4.3 now lead to the Pestov identity for all u ∈ C 2 (SM) with u • ρ = u on ∂SM.The important point to note here is that the regularity of the solution u to the transport equation is C 2 (int(SM)) ∩ Lip(SM) by Lemma 3.21.We need to prove the Pestov identity (4.7) for this class of functions.To overcome this difficulty, we use an approximation argument following [IP22, pp. 1289[IP22, pp. -1290]].
Lemma 4.4.Let (M, g) be a compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).If u ∈ C 2 (int(SM)) ∩ Lip(SM), P u ∈ L 2 (SM), and u = u • ρ on ∂SM then Proof.Following the approach taken in the proof of [IP22, Lemma 10], we extend our manifold as follows: Let M be a smooth and compact Riemannian manifold with boundary, such that M ⊂ int M. We extend the function u to a new function ũ : S M → R such that ũ satisfies ũ = u in SM, ũ ∈ C 2 (int SM) ∩ Lip(S M) and has compact support in int S M .
Following the proof of [IP22, Lemma 10], we define a sequence of mollifications (u j ) ∞ j=1 of ũ.By the basic properties of mollifiers, we have u j → u in Lip(SM) and C 2 (int SM).
By applying Lemma 2.1 and Lemma 4.1 to u j | SM , we obtain the following expression (4.9) Here, we have defined u j e = 1 2 (u j + u j • ρ) and u j o = 1 2 (u j − u j • ρ), and we have used the fact that u • ρ = u at ∂SM by assumption.
Note that Lip(∂SM) ⊂ H 1 (∂SM) and u j → u in Lip(∂SM).Therefore, the convergence also holds in H 1 (∂SM).Similar to the proof of [IP22, Lemma 10], we can conclude from H 1 convergence in SM that F u j → F u and V u j → V u in L 2 (SM).We have by the properties of mollification and regularity assumptions on u that P u j → P u in L 2 (SM).Since the all other terms but P u j in (4.9) are known to converge as j → ∞, we may conclude that lim j→∞ P u j exists and is finite.Using the commutator formula (2.7), we have P u j = P u j + X ⊥ u j − V (λ)V u j .We know that P u j → P u in L 2 (SM) by assumption and, on the other hand, X ⊥ u j → X ⊥ u and V (λ)V u j → V (λ)V u in L 2 (SM) since u ∈ Lip(SM) ⊂ H 1 (SM) and λ ∈ C ∞ (SM).This shows that P u j → P u in L 2 (SM).
By combining all of the above facts about the convergence of terms, we obtain (4.8) by taking the limit j → ∞.

Proof of Theorem
It follows from the definition of u in (3.2) that F u = −f in the interior of SM as stated in (3.3).By Lemma 3.21 and the identity V F u = −V f ∈ C 1 (SM), we know that u satisfies the assumptions of Lemma 4.4.
From (2.3) and the orthogonality (2.4), we obtain Combining the Pestov identity (4.8) with (4.10), we have Using the integration by parts formulas (2.2), we have ˆSM Applying the integrating by parts formula once again, we have An identity similar to [DP07a, p. 538] can be obtained using the commutator relations (as illustrated in (2.7)), as follows We now integrate the above equation over SM to get Notice that the identity above is a generalization of [DP07a, eq.( 14)] with the boundary term.
Hence, we have This is precisely the assertion of the proposition.
Proof.Define φ −t := r •φ − t •r, where φ − t is the dual λ-geodesic flow and r is reversion map i.e., r(x, v) = (x, −v).Since composition of smooth maps is smooth and φ −t • φ t = Id in SM, we have φ −t is a smooth map.Therefore, for each t ∈ [0, ∞), φ t : SM → SM is a diffeomorphism.This implies that the differential dφ t is invertible for every t ∈ [0, ∞).Consequently, the differential dφ −t is invertible for all t ∈ [0, ∞) and hence dφ −1 t = dφ −t .Similar to [MP11, Proof of Proposition 4.14], we apply dφ −t to both sides of (A.14).Then we obtain Differentiating both sides of the above identity with respect to t and applying the derivative formula for Lie brackets (A.13), we have Applying the commutator formulas, linearity and grouping like terms, we obtain Since dφ −t is invertible and {F (t), X ⊥ (t), V (t)} is a basis of each tangent space T φt(x,v) SM, the coefficients of F (t), X ⊥ (t) and V (t) must vanish for all t, and hence the result follows.
Lemma A.10.Let (M, g) be a Riemannian surface with or without boundary and λ ∈ C ∞ (SM).For any (x, v) ∈ SM with ξ ∈ T (x,v) SM, the differential of the λ-geodesic flow can be decomposed as where J ξ is the unique Jacobi field along the λ-geodesic π(φ t (x, v)) with the initial condition (J ξ (0), D t J ξ (0)) = ξ.
Proof.Let us consider σ : (−ǫ, ǫ) → SM be a smooth curve such that σ(0) = (x, v) and σ ′ (0) = ξ.Now we consider a variation of the λ-geodesic Using the definition of connection map, we have Combining (A.6) with (A.17), we deduce that From (A.11), we obtain This result establishes the lemma.
We will next proof an estimate needed for the proof of our main theorem.
Proof of Lemma 2.2.If M is a surface with boundary, then we extend (M, g) into a closed surface N and extend λ to a smooth function on SN (cf.[PSU23, Lemma 3.1.8]).Then the λ-geodesic flow is defined for all t ∈ R. If (2.14) holds for the closed extension N, then it also holds for M. Therefore we can assume without loss of generality, we are in a closed setting.
Let us first define Z = (Z 1 , Z 2 , Z 3 ) such that Z 1 := a, Z 2 := b, Z 3 := c.Then the equation (A.15) can be written as This is equivalent to is a bounded linear map for each t.By compactness, there is a positive constant C = C(M, g, λ) such that A γ t ≤ C/2 for all geodesics γ and all times t.Thus, we have D t |Z| 2 = 2 Z, A γ t Z ≤ C|Z| 2 .By Grönwall's inequality (cf.[Eva10, p. 624]), the above implies that |Z(t)| 2 ≤ e Ct |Z(0)| 2 for all t ≥ 0. Now if J is a λ-Jacobi field, then from (A.14) and (A.16), we have Z = (J, D t J) for some initial data, from which the claim follows.
A.3.Convexity, concavity and signed λ-curvature.Recall that II denotes the second fundamental form of the boundary ∂M, and ν(x) the inward unit normal vector to ∂M at x ∈ ∂M.Definition A.11.Let (M, g) be a Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).We say that the boundary ∂M is strictly λ-convex at a point x ∈ ∂M if the following inequality holds: for all v ∈ S x (∂M).
Similarly, we say that the boundary ∂M is strictly λ-concave at a point x ∈ ∂M if the following inequality holds: II x (v, v) < λ(x, v)iv, ν(x) for all v ∈ S x (∂M).
In [PSU23, Lemma 3.1.12],it is demonstrated that any geodesic tangent to the boundary ∂M stays in the exterior to the surface M for both small positive and negative time intervals.Moreover, the surface M is strictly magnetic convex at x ∈ ∂M, if II(x, v) > Y x (v), ν(x) , holds for all v ∈ S x (∂M) where Y x denotes the Lorentz force, see [DPSU07,Lemma A.6].An analogous property will be discussed in the next lemma in the case of a strictly λ-convex boundary which generalizes both of these results.Before stating the next lemma, let us consider the boundary defining map following [PSU23, Lemma 3.1.10].Let (N, g) be a closed extension of (M, g).Then there is a function ρ ∈ C ∞ (N), called a boundary defining function, such that ρ(x) = d(x, ∂M) near ∂M in M, and M = {x ∈ N : ρ ≥ 0}, ∂M = {x ∈ N : ρ = 0}, and N \ M = {x ∈ N : ρ < 0}.Moreover, ∇ρ(x) = ν(x) holds for all x ∈ ∂M.
Lemma A.12.Let (M, g) be a compact Riemannian surface with smooth boundary and λ ∈ C ∞ (SM).Suppose (N, g) is a closed extension of (M, g).Then ∂M is strictly λconvex if and only if any λ-geodesic in N starting from some point (x, v) ∈ ∂ 0 SM satisfies d 2 dt 2 [ρ • γ x,v (t)]| t=0 < 0. Furthermore, any λ-geodesic tangent to ∂M stays outside M for small positive and negative time intervals.Also any maximal λ-geodesic going from ∂M into M stays in the interior of M excepts for its endpoints.
In the next subsection, we will see how the signed λ-curvature for λ-geodesics is related to its dual λ-geodesic.
A.4. Curvature of the dual λ-geodesic flow.Let us define the reversion map r : SM → SM by r(x, v) = (x, −v).Let u ∈ C ∞ (SM), and consider its extension as a homogeneous function of degree zero in C ∞ (T M \ 0), denoted by u(x, y/|y|).We define the horizontal and vertical derivatives as follows: Here, Γ l ik denotes the Christoffel symbols of the metric g (cf.[IS16,p. 386]).Lemma A.13.Let (M, g) be a Riemannian surface with or without boundary.For any function f ∈ C 1 (SM), the following formulas hold: We omit the derivations of the formulas in A.13. Their proofs are straightforward computations.
Note that λ − = −λ • r by the definitions.We define K λ − (x, v) as the Gaussian λ − curvature corresponding to the λ − -geodesic flow.We prove in the next corollary that the global negative curvature assumption for the λ-geodesic flow is equivalent to the same property of the λ − -geodesic flow.
Corollary A.14.Let (M, g) be a Riemannian surface with or without boundary and λ ∈ C ∞ (SM).Then the (Gaussian) λ-curvature of the dual system satisfies for all (x, v) ∈ SM.
Proof.Let (x, v) ∈ SM.We may compute using Lemma A.13 that

2. 5 .
Lemmas for twisted geodesic flows.In the following proposition, we provide a generalized version of the Pestov identity for the generators of twisted geodesic flows.Similar identities are proved earlier in [DP07a, Theorem 3.3] for closed Riemannian surfaces, in [IS16, Lemma 8] for surfaces with boundary under the condition λ ≡ 0, and for Finslerian surfaces in terms of Lie derivatives on the boundary [AD18, Theorem 2.3].A detailed proof is given in Appendix A.1.