Next Article in Journal
Deep Learning Network Based on Improved Sparrow Search Algorithm Optimization for Rolling Bearing Fault Diagnosis
Previous Article in Journal
Efficient Estimation and Validation of Shrinkage Estimators in Big Data Analytics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differentiability of the Largest Lyapunov Exponent for Non-Planar Open Billiards

1
Department of Mathematics and Statistics, School of Physics, Mathematics and Computing, University of Western Australia, Perth, WA 6009, Australia
2
Department of Mathematics, College of Science and Arts, Najran University, Najran 66262, Saudi Arabia
Mathematics 2023, 11(22), 4633; https://doi.org/10.3390/math11224633
Submission received: 25 July 2023 / Revised: 8 November 2023 / Accepted: 10 November 2023 / Published: 13 November 2023
(This article belongs to the Section Dynamical Systems)

Abstract

:
This paper investigates the behaviour of open billiard systems in high-dimensional spaces. Specifically, we estimate the largest Lyapunov exponent, which quantifies the rate of divergence between nearby trajectories in a dynamical system. This exponent is shown to be continuous and differentiable with respect to a small perturbation parameter. A theoretical analysis forms the basis of the investigation. Our findings contribute to the field of dynamical systems theory and have significant implications for the stability of open billiard systems, which are used to model physical phenomena. The results provide a deeper comprehension of the behaviour of open billiard systems in high-dimensional spaces and emphasise the importance of taking small perturbations into consideration when analysing these systems.
MSC:
37B10; 37D50; 37D20; 34D08

1. Introduction

Billiards represent dynamic systems where a particle travels at a constant speed and collides with the boundary of the billiard’s domain following the principles of geometrical optics, which state that “the angle of incidence equals the angle of reflection”. Open billiards are a specialised category of billiard systems that occur within unbounded regions. In these scenarios, the domain encompasses the area outside a finite number of strictly convex compact obstacles. These obstacles satisfy the no-eclipse condition (H) established by Ikawa [1]. This condition ensures that the convex hull formed by any pair of obstacles does not intersect with any other obstacle, effectively preventing the existence of a straight line that passes through more than two obstacles. This condition preserves the smoothness and predictability of the open billiard system by ensuring that trajectories within the system avoid tangent points and singularities. This simplifies the analysis and understanding of the system’s behaviour, allowing for a clearer study of particle motion within the open billiard system. One way to quantify the behaviour of this system is through the use of Lyapunov exponents. Lyapunov exponents measure the rate at which neighbouring trajectories in a dynamical system diverge or converge, revealing the growth or decay of small perturbations. In the case of open billiards, the hyperbolic non-wandering set of the billiard map indicates the presence of positive and negative Lyapunov exponents. Numerous studies have investigated Lyapunov exponents for billiards, (see [2,3,4,5,6,7]). In this paper, we continue our investigation of the regularity properties of Lyapunov exponents from [8] for open billiards in Euclidean spaces. In [8], we studied the case of open billiards in the plane; here we deal with the higher-dimensional case.
Our primary findings are as follows:
The largest Lyapunov exponent for an open billiard in R n is
λ 1 = lim m 1 m i = 1 m log 1 + d i ( 0 ) i ( 0 ) ,
where d i is the distance between two reflection points along trajectories, and i is related to the curvature of the unstable manifold, as defined in (2).
To comprehend the following theorems, an understanding of non-planar billiard deformations and their concepts is required. A small deformation parameter α [ 0 , b ] , b R characterises billiard deformations, which include position, rotation, and obstacle reshaping. The initial boundary parameterisation and deformation parameter α are used to parameterise the obstacle borders. This was introduced in [9]. Billiard deformations and notations in the subsequent theorems are explained in Section 4. This study assumes that billiard deformation is a differentiable function for both parameters, providing the foundation for the theorems.
Theorem 1. 
Let K ( α ) be a C 4 , 1 billiard deformation in R n , n 3 . Let λ 1 ( α ) be the largest Lyapunov exponent for K ( α ) . Then, the largest Lyapunov exponent is continuous as a function of α.
Theorem 2. 
Let K ( α ) be a C 5 , 2 billiard deformation in R n , n 3 . Let λ 1 ( α ) be the largest Lyapunov exponent for K ( α ) . Then λ 1 ( α ) is C 1 with respect to α.
The proofs of these theorems are in Section 7 and Section 8, respectively.
The demonstration of the continuity and differentiability of the largest Lyapunov exponent for open billiards is highly significant. Open billiards represent a chaotic physical system (see, e.g., [1,10,11]), which appears naturally in scattering theory and tomography. Unlike many systems, which often exhibit complex and unpredictable behaviour in response to perturbations, open billiards show a more regular response, offering insights into the predictability of chaotic systems’ dynamics. This understanding is critical to improving our understanding of complex systems and their applications.

2. Preliminaries

2.1. Open Billiards

Let K 1 , K 2 , , K z 0 be strictly convex, disjoint, and compact subsets of R n , n 3 with smooth boundaries K i , and satisfying condition (H) of Ikawa [1]: for any i j k the convex hull of K i K k does not have any common points with K j . Let Ω be the exterior of K (i.e., Ω = R n K ¯ ), where K = i K i . Let
M = { x = ( q , v ) i n t ( Ω ) × S n 1 o r ( q , v ) Ω × S n 1 : ν K ( q ) , v 0 } ,
where n K ( q ) is the outwards unit normal vector to K at q. For t R and x M , the billiard flow Φ t is defined as Φ t ( x ) = ( q t , v t ) , where q t and v t represent the position and velocity of the x at time t. Let Λ represent the set of all points of M that have bounded trajectories. Let
M = { x = ( q , v ) ( K × S n 1 ) M : ν ( q ) , v 0 } .
Let t j R , where j Z , represent the time of j-th reflection. The billiard map  B : M M is defined as B ( q 0 , v 0 ) = ( q 1 , v 1 ) where q 1 = q 0 + t 1 v 0 K j and v 1 = v 0 2 v 0 , ν ( q 1 ) ν ( q 1 ) . Clearly, B is a smooth diffeomorphism on M. Define the canonical projection map π : M K by π ( q , v ) = q . Let M 0 M be the non-wandering set of the billiard map B, that is, M 0 = { x M :   | t j ( x ) |   < , f o r a l l j Z } . It is clear that M 0 is invariant under B. See [3,4,12,13,14] for general information about billiard dynamical systems.

2.2. Symbolic Coding for Open Billiards

Each particular x M 0 can be coded by a bi-infinite sequence
ξ ( x ) = ( , ξ 1 , ξ 0 , ξ 1 , ) { 1 , 2 , , z 0 } Z ,
in which ξ i ξ i + 1 , for all i Z , and ξ j indicates the obstacle K ξ j such that π B j ( x ) K ξ j . For example, if there are three obstacles K 1 , K 2 , K 3 and K 4 as above and a particular x repeatedly hits K 2 , K 1 , K 4 , K 3 , K 2 , K 1 , K 4 , K 3 , then the bi-infinite sequence is ( , 2 , 1 , 4 , 3 , 2 , 1 , 4 , 3 , ) . Let Σ be the symbol space, which is defined as:
Σ = { ξ = ( , ξ 1 , ξ 0 , ξ 1 , ) { 1 , 2 , , z 0 } Z : ξ i ξ i + 1 , i Z } .
Define the representation map  R : M 0 Σ by R ( x ) = ξ ( x ) . Let σ : Σ Σ be the two-sided subshift map defined by σ ( ξ ) = ξ where ξ i = ξ i + 1 .
It is known that the representation map R : M 0 Σ is a homeomorphism (see, e.g., [12]). See [1,10,11,12,15], for topics related to symbolic dynamics for open billiards.

2.3. Lyapunov Exponents

For the open billiard B : M 0 M 0 in R n we will use the coding R : M 0 Σ from Section 2.2, which conjugates B with the shift map σ : Σ Σ , to define Lyapunov exponents. It follows from the symbolic coding that there are ergodic σ -invariant measures μ on Σ . Let μ be an ergodic σ -invariant probability measure on Σ . The following is a consequence of Oseledets Multiplicative Ergodic Theorem (see, e.g., [16,17]):
Theorem 3 (A Consequence of Oseledets Multiplicative Ergodic Theorem).
There exist real numbers λ 1 > λ 2 > > λ k > 0 > λ k > > λ 1 and vector subspaces E j u ( x ) and E j s ( x ) of T x ( M ) , 1 j k n 1 , x M 0 , depending measurably on R ( x ) Σ such that:
1. 
E u ( x ) = E 1 u ( x ) E k u ( x ) and E s ( x ) = E 1 s ( x ) E k s ( x ) for almost all x M ;
2. 
D x B ( E i u ( x ) ) = E i u ( B ( x ) ) and D x B ( E i s ( x ) ) = E i s ( B ( x ) ) for all x M 0 and all i = 1 , , k , and
3. 
For almost all x M 0 there exists
lim m 1 m log D x B m ( w ) = λ j
whenever w E j u ( x ) E k u ( x ) ( j k ), however, w E j + 1 u ( x ) E k u ( x ) .
And
lim m 1 m log D x B m ( w ) = λ j
whenever w E j s ( x ) E k s ( x ) ( j k ), however, w E j + 1 s ( x ) E k s ( x ) .
Here, “for almost all x” means “for almost all x” with respect to μ . The numbers λ 1 , , λ k are called Lyapunov exponents, while the invariant subspaces E j ( x ) are called Oseledets subspaces.

2.4. Propagation of Unstable Manifolds for Open Billiards

This part explains the relationship between unstable manifolds for the billiard ball map and the billiard flow in R n . Recall that M 0 represents the non-wandering set of the billiard ball map, while Λ represents the non-wandering set of the billiard flow. For the billiard map, the unstable manifolds of size ϵ are
W ϵ u ( x ) = { y M : d ( B n ( x ) , B n ( y ) ) ϵ for all n N , d ( B n ( x ) , B n ( y ) ) n 0 } .
Similarly, for x = ( q , v ) Λ the unstable manifolds W ˜ ϵ u ( x ) for the billiard ball flow are
W ˜ ϵ u ( x ) = { y M : d ( Φ t ( x ) , Φ t ( y ) ) ϵ for all t 0 , d ( Φ t ( x ) , Φ t ( y ) ) t 0 } .
It is known that the unstable manifolds for the billiard ball map and the billiard flow naturally are related. This correspondence can be described as follows geometrically. Given a point x 0 = ( q 0 , v 0 ) in the non-wandering set M 0 , and a small number 0 < t ( 1 ) < t 1 , let y 0 = ( q 0 + t ( 1 ) v 0 , v 0 ) . Then, there exists a one-to-one correspondence map φ ˜ 1 between the unstable manifold W u ( x 0 ) and the unstable manifold W ˜ u ( y 0 ) . In addition, it follows from Sinai [13,14], that the unstable manifold W ˜ u ( y 0 ) has the form
W ˜ u ( y 0 ) = { ( p , n Y 0 ( p ) ) : p Y 0 } ,
where Y 0 is a smooth ( n 1 ) -dimensional hypersurface in R n containing the point y 0 and is strictly convex with respect to the unit normal field n Y 1 ( p ) .
Likewise, for all j = 1 , 2 , , m , there exists one-to-one correspondence φ ˜ j between the unstable manifolds W u ( x j ) and W ˜ u ( y j ) , where B j ( x 0 ) = x j = ( q j , v j ) and y j = ( q j + t ( j ) v j , v j ) , for a small positive t j 1 < t ( j ) < t j t j 1 . Moreover, W ˜ u ( y j ) takes the form Y j ˜ = { ( p j , n Y j ( p j ) ) ; j Y j } , where Y j is also a smooth ( n 1 ) -dimensional hypersurface in R n containing the point y j and is strictly convex with respect to the unit normal field n Y j ( p j ) .
The following are the commutative diagrams involving the unstable manifolds for open billiard maps and for the billiard flows, as shown above:
W u ( x 0 ) B ( x 0 ) W u ( x 1 ) B 2 ( x 0 ) W u ( x 2 ) W u ( x m 1 ) B m ( x 0 ) W u ( x m ) φ ˜ 0 φ ˜ 1 φ ˜ 2 φ ˜ m 1 φ ˜ m W ˜ u ( y 0 ) Φ τ 1 ( y 0 ) W ˜ u ( y 1 ) Φ τ 2 ( y 0 ) W ˜ u ( y 2 ) W ˜ u ( y m 1 ) Φ τ m ( y 0 ) W ˜ u ( y m )
where τ j = t j + t ( j + 1 ) . Similarly, the following are the commutative diagrams involving the corresponding tangent spaces of unstable manifolds under the derivative of the billiard ball maps and the derivative of the billiard flows:
E u ( x 0 ) D B E u ( x 1 ) D B 2 E u ( x 2 ) E u ( x m 1 ) D B m E u ( x m ) D φ ˜ 0 D φ ˜ 1 D φ ˜ 2 D φ ˜ m 1 D φ ˜ m E ˜ u ( y 0 ) D Φ τ 1 E ˜ u ( y 1 ) D Φ τ 2 E ˜ u ( y 2 ) E ˜ u ( y m 1 ) D Φ τ m E ˜ u ( y m )
where E u ( x j ) = T x j ( W u ( x j ) ) and E ˜ u ( y j ) = T y j ( W ˜ u ( y j ) ) , and E ˜ u ( y j ) v j . Since the derivatives D φ ˜ j 1 and D φ ˜ j , for all j, are uniformly bounded [13], then there exist global constants C > c > 0 such that
c D Φ τ m D B m C D Φ τ m .
This will be used to calculate the largest Lyapunov exponent λ 1 for the open billiard map in Section 3.
Based on the prior discussion, we now apply the concept of unstable manifold propagation to write the main theorem, which involves the propagation of an appropriate convex curve on a convex hypersurface. This theorem was proved in reference [18] and is utilised to calculate the largest Lyapunov exponent, λ 1 .
Let x = ( q , v ) M 0 and let W ϵ u ( x 0 ) be the local unstable manifold for x 0 for sufficiently small ϵ > 0 . Then, W ϵ u ( x 0 ) = { ( x , n X ( x ) ) : x X } , where X is a convex curve on a smooth hypersurface X ˜ containing q such that X ˜ is strictly convex with respect to the unit normal field n X ˜ . This follows from (cf. [13,14]), see also [18]. Let K x : T q X ˜ T q X ˜ be the curvature operator (second fundamental form). Since X ˜ is strictly convex, then the curvature K x is positive definite with respect to the unit normal field n X ˜ .
Let X be parameterised by q ( s ) , s [ 0 , a ] , such that q ( 0 ) = q 0 + r v 0 for a small r > 0 , and by the unit normal field n X ( q ( s ) ) . Let q j ( s ) , j = 1 , 2 , , m be the j-th reflection points of the forward billiard trajectory γ ( s ) generated by x ( s ) = ( q ( s ) , n X ( q ( s ) ) . We assume that a > 0 is sufficiently small so that the j-th reflection points q j ( s ) belong to the same boundary component K ξ j for every s [ 0 , a ] . Let 0 < t 1 ( x ( s ) ) < < t m ( x ( s ) ) be the times of the reflections of the ray γ ( s ) at K . Let κ j ( s ) be the curvature of K at q j ( s ) and ϕ j ( s ) be the collision angle between the unit normal ν j ( s ) to K and the reflection ray of γ ( s ) at q j ( s ) . Also, let d j ( s ) be the distance between two reflection points, i.e., d j ( s ) = q j + 1 ( s ) q j ( s ) , j = 0 , 1 , , m .
Given a large m 1 , let t m ( x ( s ) ) < t < t m + 1 ( x ( s ) ) . Set X ^ = { ( q ( s ) , n X ( s ) ) : s [ 0 , a ] } , and Φ t ( X ^ ) = X ^ t . Let π ( Φ t ( x ( s ) ) ) = p ( s ) . Then p ( s ) , s [ 0 , a ] , is a parameterisation of the C 3 curve X t = π ( Φ t ( X ^ ) ) .
Let k j ( s ) be the normal curvature of X t j ( s ) = lim t t j ( s ) X t at q j ( s ) in the direction w ^ j ( s ) ( w ^ j ( s ) = 1 ) of lim t t j ( s ) ( d / d s ) w t ( s ) | s = s where w ^ j ( s ) = q ˙ t ( s ) / q ˙ t ( s ) . For j 0 let
K j ( s ) : T q j ( s ) ( X t j ( s ) ) T q j ( s ) ( X t j ( s ) )
be the curvature operator of X t j ( s ) at q j ( s ) , and define l j ( s ) > 0 by
1 + d j ( s ) j ( s ) 2 = 1 + 2 d j ( s ) k j ( s ) + ( d j ( s ) ) 2 K j ( s ) w ^ j ( s ) 2 .
Set
δ j ( s ) = 1 1 + d j ( s ) j ( s ) , 1 j m .
Theorem 4 ([18]). 
For all s [ 0 , a ] we have
q ˙ ( s ) = p ˙ ( s ) δ 1 ( s ) δ 2 ( s ) δ m ( s ) .

2.5. Curvature of Unstable Manifolds

Following [13,14,19] (see also [20]), here we express the curvature operator K j of the convex front X j at q j ( s ) using certain related objects:
  • X ˜ j represents the convex front passing q j ( s ) before collision, that is X ˜ j = lim τ t j X ˜ τ , where t j 1 < τ < t j .
  • X ˜ j + represents the convex front passing q j ( s ) after collision, that is X ˜ j + = lim τ ˜ t j X ˜ τ ˜ , where t j < τ ˜ < t j + 1 . We write X ˜ j to indicate to X ˜ j + .
  • J is the hyperplane of X ˜ j at q j ( s ) , i.e., J = T q j ( X ˜ j ) , which is perpendicular to v j 1 .
  • J + is the hyperplane of X ˜ j + at q j ( s ) , i.e., J + = T q j ( X ˜ j + ) , which is perpendicular to v j .
  • T is the hyperplane of K at q j ( s ) , i.e., T = T q j ( K ) , which is perpendicular to ν j ( q j ( s ) ) .
Now, the curvature operator K j : J + J + is given by
K j = U j 1 K j U j + 2 cos ϕ j V j N j V j ,
where
  • K j is the curvature operator of X ˜ j , which defined as K j = K j 1 ( I + d j 1 K j 1 ) 1 or we can write K j = ( K j 1 1 + d j 1 I ) 1 .
  • The unitary operator U j : J j + J j is a projection parallel to ν j ( q j ( s ) ) , defined as; for all w j J j +
    U j w j = w j w j , v j 1 cos ϕ j ν j .
  • U j 1 : J j J j + is a projection parallel to ν j ( q j ( s ) ) , defined as; for all w j 1 J j
    U j 1 w j 1 = w j 1 w j 1 , v j cos ϕ j ν j .
  • V j : J j + T j is a projection parallel to v j , defined as; for all w j J j +
    V j w j = w j w j , ν j cos ϕ j v j .
  • V j : T j J j + is a projection parallel to ν j ( q ( s ) ) , defined as; for all u ˜ j T j
    V j u ˜ j = u ˜ j u ˜ j , v j cos ϕ j ν j .
  • N j is the curvature operator s.f.f. of K at q j ( s ) .
The operator V N V is bounded by κ min V N V κ max cos 2 ϕ j where κ min and κ max are the minimum and maximum eigenvalues of the normal curvature N and again ϕ j ( s ) is the collision angle at q j ( s ) and ϕ j ( s ) [ 0 , π / 2 ] . Let μ j ( s ) = μ j ( x ( s ) ) and η j ( s ) = η j ( x ( s ) ) be the eigenvalues of the curvature operator K j . Then by using these and (5), we obtain
μ min 2 κ min μ j + 1 ( s ) η j + 1 ( s ) η j ( s ) 1 + d j ( s ) η j ( s ) + 2 κ max cos ϕ j ( s ) 1 d min + 2 κ max cos ϕ max η max .

3. Estimation of the Largest Lyapunov Exponent for Non-Planar Open Billiards

Here, we want to estimate the largest Lyapunov exponent λ 1 for non-planar open billiards. We use Oseledets multiplicative ergodic Theorems 3 and 4.
Assume that μ is an ergodic σ -invariant measure on Σ , and let x = ( q , v ) M 0 correspond to a typical point in Σ with respect to μ via the representation map R. As in Theorem 3, there exists a subset A 0 of Σ with μ ( A 0 ) = 1 such that
λ 1 = lim m 1 m log D x B m ( w ) ,
with w E 1 u ( x ) E 2 u ( x ) E k u ( x ) E 2 u ( x ) E 3 u ( x ) E k u ( x ) .
Let w = q ˙ ( s ) as in Section 4. And then by using (1), there exist some global constants c 1 > c 2 > 0 , independent of x 0 , X, m, etc., such that
c 2 p ˙ ( s )     D x 0 B m ( w )   c 1 p ˙ ( s )
for all s [ 0 , a ] . So, by (4),
c 2 δ 1 ( 0 ) δ 2 ( 0 ) δ m ( 0 ) D x 0 B m ( w ) c 1 δ 1 ( 0 ) δ 2 ( 0 ) δ m ( 0 )
for all s [ 0 , a ] . Using this for s = 0 , taking logarithms and limits as m , we obtain
lim m 1 m log δ 1 ( 0 ) δ 2 ( 0 ) δ m ( 0 ) lim m 1 m D x 0 B m ( w ) lim m 1 m log δ 1 ( 0 ) δ 2 ( 0 ) δ m ( 0 ) .
Hence,
λ 1 = lim m 1 m i = 1 m log δ i ( 0 ) = lim m 1 m i = 1 m log 1 + d i ( 0 ) i ( 0 ) .
We can use (2) and the (10) to estimate the upper and lower of the largest Lyapunov exponent λ 1 as follows:
1 + d j ( s ) j ( s ) 2 1 + 2 d j ( s ) η j ( s ) + ( d j ( s ) ) 2 ( η j ( s ) ) 2 = 1 + d j ( s ) η j ( s ) 2 .
This implies that
1 + d j ( s ) j ( s ) 1 + d j ( s ) η j ( s ) 1 + d max η max .
In the same way, we obtain
1 + d j ( s ) j ( s ) 1 + d j ( s ) μ j ( s ) 1 + d min μ min .
Therefore,
log ( 1 + d min μ min ) λ 1 log ( 1 + d max η max ) .

4. Billiard Deformations in R n

Let α I = [ 0 , b ] , for some b R + , be a deformation parameter and let K j ( α ) be parameterised counterclockwise by φ j ( u j ( 1 ) , u j ( 2 ) , , u j ( n 1 ) , α ) . Let q j = φ j ( u j ( 1 ) , u j ( 2 ) , , u j ( n 1 ) , α ) be a point that lies on K j ( α ) . Denote the perimeter of K i ( α ) by L j ( α ) , and let R j = { ( u j ( 1 ) , u j ( 2 ) , , u j ( n 1 ) : α I , u j ( t ) [ 0 , L j ( α ) ] } .
Definition 1 ([9]). 
For any α I = [ 0 , b ] , let K ( α ) be a subset of R n , n 3 . For integers r 2 , r 1 , we call K ( α ) a C r , r -billiard deformation if the following conditions hold for all α I :
1. 
K ( α ) = i = 1 z 0 K i ( α ) satisfies the no-eclipse condition ( H ) .
2. 
Each K i ( α ) is a compact, strictly convex set with C r boundary, and K i ( α ) K j ( α ) = ϕ for i j .
3. 
For each i = 1 , 2 , , m and all p K i ( α ) , there is a rectangle R p R n 1 and a C r , r function φ j : R p R n , which is an orthonormal parametrisation of K i ( α ) at p.
4. 
For all integers 0 l r , 0 l r (apart from l = l = 0 ), there exist constants C φ ( l , l ) depending only on the choice of the billiard deformation and the parameterisation φ j , such that for all integers j = 1 , 2 , 3 , , z 0 ,
l α l i l φ j C φ ( l , l ) .
We consider the open billiard deformation map, denoted as B α , defined on the non-wandering set M α for K ( α ) . As in Section 2.2, we define Σ α and R α as the mapping from M α to Σ α such that R α ( x ( α ) ) = ξ ( x ( α ) ) . Using the parameterisation defined earlier, we can express the point corresponding to deformed billiard trajectories as q ξ j ( α ) = φ ξ j ( u j ( 1 ) ( α ) , u j ( 2 ) ( α ) , , u j ( n 1 ) ( α ) , α ) K ξ j ( α ) , where u ξ j ( t ) ( α ) [ 0 , L ξ j ( α ) ] . We write q j ( α ) = φ j ( u j ( 1 ) ( α ) , u j ( 2 ) ( α ) , , u j ( n 1 ) ( α ) , α ) for brevity.
It was shown in [9] that u j ( t ) ( α ) = u ξ j ( t ) ( α ) , where t = 1 , 2 , , n 1 , for a fixed ξ Σ α , is differentiable with respect to α and its higher derivative is bounded by a constant independent of α and j.
Theorem 5. 
[9] Let K ( α ) be a C r , r billiard deformation with r 2 , r 1 . Then for all t = 1 , 2 , , n 1 , u j ( t ) ( α ) is C min { r 1 , r } with respect to α, and there exist constants C u ( s ) > 0 such that
s u j ( t ) ( α ) α s C u ( s ) .

5. Propagation of Unstable Manifolds for Non-Planar Billiard Deformations

Here, we want to re-describe the propagation of the unstable manifold mentioned in Section 4 for the billiard deformation K ( α ) with respect to the deformation parameter α I = [ 0 , b ] . Let K ( α ) be a C r , r billiard deformation in R n , with r 4 , r 1 . Let X ˜ be a smooth hypersurface passing through x 0 ( α ) = ( q 0 ( α ) , v 0 ( α ) ) M α and let X be a C 3 convex curve, with respect to the unit normal field ν X , on X ˜ . Let X be parameterised by q ( s , α ) for all s [ 0 , a ] and set X ^ α = { ( q ( s , α ) , ν X ( q ( s , α ) ) ) , s [ 0 , a ] } . For j = 1 , 2 , , m , let q j ( s , α ) be the reflection points generated by x ( s , α ) = ( q ( s , α ) , ν ( q ( s , α ) ) ) , d j ( s , α ) be the distance between q j ( s , α ) and q j + 1 ( s , α ) , and ϕ j ( s , α ) be the angle of reflection at q j ( s , α ) .
Let t m 1 ( x ( s , α ) ) < t < t m ( x ( s , α ) ) for some large m 2 . Let π ( Φ t ( X ^ α ) ) = X α t . Then, X α t is a C 3 curve on X ˜ t parameterised by u t ( s , α ) .
Let K j ( s , α ) be the curvature operator of X α t j ( s ) = lim t t j ( x ( s ) ) X α t at q j ( s , α ) in the direction w ^ j ( s , α ) of lim t t j ( ( x ( s ) ) ( d / d s ) w ^ t ( s , α ) ) | ( s , α ) . Let K 0 ( s ) be the curvature operator of X at q ( s , α ) , which is independent of α . And, for j 1 let
K j ( s , α ) : T q j ( s , α ) ( X α t j ( x ( s ) ) ) T q j ( s , α ) ( X α t j ( x ( s , α ) ) )
be the curvature operator of X t j ( x ( s , α ) ) at q j ( s , α ) . Define j ( s , α ) by
1 + d j ( s , α ) j ( s , α ) 2 = w j ( s , α ) + d j ( s , α ) K j ( s , α ) 2 , for all j = 0 , 1 , , m .
Then
1 + d j ( s , α ) j ( s , α ) 2 = 1 + 2 d j ( s , α ) k j ( s , α ) + ( d j ( s , α ) ) 2 K j ( s ) w ^ j ( s ) 2 ,
where k j ( s , α ) is the normal curvature of X α j ( s , α ) at q j ( s , α ) in the direction w ^ j ( s , α ) , which is given by k j ( s , α ) = k j ( w ^ j ( s , α ) ) = K j ( s , α ) w ^ j ( s , α ) , w ^ j ( s , α ) .
Now, we want to re-write the curvature operator in (5) with respect to the billiard deformation parameter α . So, we have for j = 1 , 2 , , m
K j ( s , α ) = U j 1 ( s , α ) ( K j 1 1 ( s , α ) + d j 1 ( s , α ) I ) 1 U j ( s , α ) + 2 cos ϕ j ( s , α ) V j ( s , α ) N j ( s , α ) V j ( s , α )
For brevity, we will write all previous characteristics as, e.g., K ( 0 , α ) = K ( α ) , d ( 0 , α ) = d ( α ) , …, in the case s = 0 . We can write (16) as follows
K j ( α ) = U j 1 ( α ) ( K j 1 1 ( α ) + d j 1 ( α ) I ) 1 U j ( α ) + 2 Θ j ( α ) ,
where Θ j ( α ) = cos ϕ j V j N j V j . Note that all terms in the last formula are functions of α , and are defined as in Section 2.5 with respect to α .

5.1. Estimates of the Higher Derivative of Billiard Characteristics in R n

This section aims to demonstrate the differentiability of the billiard deformation characteristics in high dimensions with respect to α . These characteristics are described in Section 5. Furthermore, we establish that the derivatives of these characteristics are bounded by constants independent of deformation parameter α [ 0 , b ] and the number of reflections j Z + . In particular, we show that the first and second derivatives are bounded, which holds significant relevance for the subsequent Section 7 and Section 8. The higher derivatives are bounded via induction. All corollaries, which are provided here, are based on Definition 1 and Theorem 5.
Corollary 1. 
Let K ( α ) be a C r , r billiard deformation with r 2 , r 1 , in R n with n 3 . Let q j ( α ) belong to K ξ j . Then q j ( α ) is C s , where s = min { r 1 , r } , with respect to α, and there exist constants C q ( s ) > 0 such that
d o r d e r = s q j ( α ) d α o r d e r = s C q ( s ) .
Proof. 
Let q j ( α ) = φ j u j ( 1 ) ( α ) , u j ( 2 ) ( α ) , , u j ( n 1 ) ( α ) , α K ξ j . Then, from Definition 1 and Theorem 5, q j ( α ) is C min { r 1 , r } with respect to α . For the first derivatives of q j ( α ) we have
d q j d α = i = 1 n 1 φ j u j ( i ) u j ( i ) α + φ j α .
From Definition 1 and Theorem 5, there exists C q ( 1 ) > 0 such that
d q j d α ( n 1 ) C u ( 1 ) + C φ ( 0 , 1 ) = C q ( 1 ) .
For the second derivative, we have
d o r d e r = 2 q j d α o r d e r = 2 = i = 1 n 1 k = 1 n 1 2 φ j u j ( k ) u j ( i ) u j ( i ) α u j ( k ) α + i = 1 n 1 φ j u j ( i ) o r d e r = 2 u j ( i ) α o r d e r = 2 + 2 i = 1 n 1 φ j u j ( i ) α u j ( i ) α + o r d e r = 2 φ j α o r d e r = 2 .
As before, there exists C q ( 2 ) > 0 such that
d o r d e r = 2 q j d α o r d e r = 2 ( n 1 ) 2 C φ ( 2 , 0 ) ( C u ( 1 ) ) 2 + ( n 1 ) C u ( 2 ) + 2 ( n 1 ) C φ ( 1 , 1 ) C u ( 1 ) + C φ ( 2 , 0 ) = C q ( 2 ) .
This constant is independent of j and α . Continuing by induction, we can see that the s -th derivative of q j ( α ) is bounded by a constant that depends only on s and n. Thus, the statement is proved. □
Corollary 2. 
Let K ( α ) be a C r , r billiard deformation in R n with r 2 , r 1 . Then d j ( α ) is C s , where s = min { r 1 , r } , with respect to α, and there exist constants C v ( s ) > 0 depending only on s , and n such that
| d o r d e r = s d j ( α ) d α o r d e r = s | C d ( s ) .
Proof. 
Let d j = q j + 1 ( α ) q j ( α ) , where q j ( α ) = φ j u j ( 1 ) ( α ) , u j ( 2 ) ( α ) , , u j ( n 1 ) ( α ) , α belongs to K ξ j . Then, from Corollary 1, d j ( α ) is C min { r 1 , r } and its derivative with respect to α is
d d j d α = φ j + 1 φ j φ j + 1 φ j , φ j + 1 α + i = 1 n 1 φ j + 1 u j + 1 ( i ) u j + 1 ( i ) α φ j α i = 1 n 1 φ j u j ( i ) u j ( i ) α .
By using the estimations in Definition 1 and Theorem 5, we obtain
| d d j d α | 2 C φ ( 0 , 1 ) + 2 ( n 1 ) C u ( 1 ) = C d ( 1 ) ,
where C d ( 1 ) > 0 is a constant depending only on n, and its second derivative is
d o r d e r = 2 d j d α o r d e r = 2 = φ j + 1 α + i = 1 n 1 φ j + 1 u j + 1 ( i ) u j + 1 ( i ) α φ j α i = 1 n 1 φ j u j ( i ) u j ( i ) α 2 φ j + 1 φ j ( φ j + 1 φ j ) 2 φ j + 1 α + i = 1 n 1 φ j + 1 u j + 1 ( i ) u j + 1 ( i ) α φ j α i = 1 n 1 φ j u j ( i ) u j ( i ) α 2 φ j + 1 φ j 3 + φ j + 1 φ j φ j + 1 φ j · ( 2 φ j + 1 α 2 + k = 1 n 1 i = 1 n 1 2 φ j + 1 u j + 1 ( k ) u j + 1 ( i ) u j + 1 ( i ) α u j + 1 ( k ) α + i = 1 n 1 φ j + 1 u j + 1 ( i ) 2 u j + 1 ( i ) α 2 + i = 1 n 1 2 φ j + 1 u j + 1 ( i ) α u j + 1 ( i ) α 2 φ j α 2 k = 1 n 1 i = 1 n 1 2 φ j u j ( k ) u j ( i ) u j ( i ) α u j ( k ) α i = 1 n 1 φ j u j ( i ) 2 u j ( i ) α 2 i = 1 n 1 2 φ j u j ( i ) α u j ( i ) α ) .
This is bounded by a constant C d ( 2 ) depending only on n such that
| d o r d e r = 2 d j d α o r d e r = 2 | 2 ( C d ( 1 ) ) 2 d min + 2 C φ ( 0 , 2 ) + 2 ( n 1 ) 2 C φ ( 2 , 0 ) ( C u ( 1 ) ) 2 + 2 ( n 1 ) C u ( 2 ) + 2 ( n 1 ) C φ ( 1 , 1 ) C u ( 1 ) = C d ( 2 ) .
Continuing by induction, we can see that the s -th derivative of d j ( α ) is bounded by a constant which depends only on n and s . This proves the statement. □
Corollary 3. 
Let K ( α ) be a C r , r billiard deformation in R n with r 2 , r 1 . Let v j ( α ) be the unit speed vector from q j ( α ) to q j + 1 ( α ) . Then, v j ( α ) is C s , where s = min { r 1 , r } , with respect to α, and there exist constants C v ( s ) > 0 depending only on n and s , such that
d o r d e r = s v j ( α ) d α o r d e r = s C v ( s ) .
Proof. 
We can write v j ( α ) = q j + 1 ( α ) q j ( α ) d j ( α ) . And then by using Corollaries 1 and 2, the statement is proved. □
Corollary 4. 
Let K ( α ) be a C r , r billiard deformation in R n with r 3 , r 1 , n 3 . Let ν j ( α ) be the normal vector field to the boundary of obstacle K ξ j at the point q j ( α ) . Then, ν j ( α ) is at least C s , where s = min { r 2 , r } , and there exist constants C ν ( s ) > 0 depending only on n and s such that
d o r d e r = s ν j d α o r d e r = s C ν ( s ) .
Proof. 
Let ν j ( α ) the normal vector field to the boundary of obstacle K ξ j at the point q j ( α ) . We can write
ν j ( α ) = φ j u j ( 1 ) × φ j u j ( 2 ) × × φ j u j ( n 1 ) .
Then, from Definition 1 and Theorem 5, ν j is at least C min { r 2 , r } with respect to α . To show its derivatives are bounded, we have
d ν j d α = i = 1 n 1 2 φ j u j ( i ) u j ( 1 ) u j ( i ) α + 2 φ j u ( 1 ) α × φ j u j ( 2 ) × × φ j u j n 1 + + φ j u j ( 1 ) × × i = 1 ( n 1 ) 2 φ j u j ( i ) u j ( n 1 ) u j ( i ) α + 2 φ j u j ( n 1 ) α .
Then, there exists a constant C ν ( 1 ) depending only on n such that
d ν j d α ( n 1 ) ( n 1 ) C φ ( 2 , 0 ) C u ( 1 ) + C φ ( 1 , 1 ) = C ν ( 1 ) .
Next, we want to show that the second derivative of ν j ( α ) is bounded by a constant.
To simplify the last derivative, we can write d ν j d α = Ψ 1 + Ψ 2 + + Ψ n 1 , where Ψ i corresponds to one square bracket […]. The first derivative of Ψ 1 with respect to α is
d Ψ 1 d α = [ ( k = 1 n 1 i = 1 n 1 3 φ j u j ( k ) u j ( i ) u j ( 1 ) u j ( i ) α u j ( k ) α + i = 1 n 1 2 φ j u j ( i ) u j ( 1 ) 2 u j ( i ) α 2 + 2 i = 1 n 1 3 φ j u j ( i ) u j ( 1 ) α u j ( i ) α + 3 φ j u j ( 1 ) α 2 ) × × φ j u j ( n 1 ) ] + i = 1 n 1 2 φ j u j ( i ) u j ( 1 ) u j ( i ) α + 2 φ j u ( 1 ) α × i = 1 n 1 2 φ j u j ( i ) u j ( 2 ) u j ( i ) α + 2 φ j u ( 1 ) α × × φ j u j ( n 1 ) + + i = 1 n 1 2 φ j u j ( i ) u j ( 1 ) u j ( i ) α + 2 φ j u ( 1 ) α × φ j u j ( 2 ) × × i = 1 n 1 2 φ j u j ( i ) u j ( n 1 ) u j ( i ) α + 2 φ j u ( n 1 ) α .
And
d Ψ 1 d α ( n 1 ) 2 C φ ( 3 , 0 ) ( C u ( 1 ) ) 2 + ( n 1 ) C φ ( 2 , 0 ) C u ( 2 ) + 2 ( n 1 ) C φ ( 2 , 1 ) C u ( 1 ) + C φ ( 1 , 2 ) + ( n 1 ) ( n 1 ) C φ ( 2 , 0 ) C u ( 1 ) + C φ ( 1 , 1 ) 2 = C Ψ .
We can obtain similar estimates for d Ψ 2 d α , , d Ψ n 1 d α . Therefore, there exists a constant C ν ( 2 ) > 0 depends only on n such that
d o r d e r = 2 ν j d α o r d e r = 2 ( n 1 ) C Ψ = C ν ( 2 ) .
So, we can see by induction that the s -th derivative, where s = min { r 2 , r } , is bounded by a constant C ν ( s ) > 0 , which depends only on n and s . This proves the statement. □
Corollary 5. 
Let K ( α ) be a C r , r billiard deformation in R n with r 3 , r 1 , n 3 . Let ϕ j ( α ) the collision angle at q j ( α ) . Then, cos ϕ j ( α ) is C s , where s = min { r 2 , r } , and there exist constants C ϕ ( s ) > 0 depending only on n and s such that
| d o r d e r = s cos ϕ j d α o r d e r = s | C ϕ ( s ) .
Proof. 
Recall that cos ϕ j ( α ) = v j ( α ) , ν j ( α ) . Then, by using Corollaries 3 and 4, the statement is proved. □
Corollary 6. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 , n 3 . Let N j ( α ) be the shape operator of K ξ j at the point q j ( α ) . Then N j ( α ) is C s , where s = min { r 3 , r } , and there exist constants C κ ( s ) > 0 , which depend only on s and n such that
d o r d e r = s N j d α o r d e r = s C N ( s ) .
Proof. 
The shape operator N j : T ( K ξ j ( α ) ) T ( K ξ j ( α ) ) at the point q j ( α ) , is defined by N j = ν j , where ν j is the normal vector field to the boundary of obstacle K ξ j at the point q j ( α ) . First, from the expression of ν j , we have
ν j = φ j u j ( 1 ) × φ j u j ( 2 ) × × φ j u j ( n 1 ) . = φ j u j ( 1 ) × φ j u j ( 2 ) × × φ j u j ( n 1 ) + . + φ j u j ( 1 ) × φ j u j ( 2 ) × × φ j u j ( n 1 ) .
And
φ j u j ( 1 ) = i = 1 n 1 2 φ j ( 1 ) u j ( i ) u j ( 1 ) , , i = 1 n 1 2 φ j ( n ) u j ( i ) u j ( 1 ) = i = 1 n 1 2 φ j ( k ) u j ( i ) u j ( 1 ) k = 1 n .
This is bounded by a constant C 1 ( 1 ) depending only on n such that
φ j u j ( 1 ) n ( n 1 ) C φ ( 2 , 0 ) = C 1 .
Also, φ j u j ( 1 ) is C min { r 3 , r } with respect to α and its first derivative is
d d α φ j u j ( 1 ) = d d α i = 1 n 1 2 φ j ( k ) u j ( i ) u j ( 1 ) k = 1 ( n ) = i = 1 n 1 p = 1 n 1 3 φ j ( k ) u j ( p ) u j ( i ) u j ( 1 ) u j ( p ) α + i = 1 n 1 3 φ j ( k ) u j ( i ) u j ( 1 ) α k = 1 ( n ) .
Thus, there exists a constant C 1 ( 1 ) depending only on n such that
d d α φ j u j ( 1 ) n ( n 1 ) 2 C φ ( 3 , 0 ) C u ( 1 ) + ( n 1 ) C φ ( 2 , 1 ) = C 1 ( 1 ) .
And then, by deriving (19) with respect to α and using the estimations in (18), (20) and (22), we obtain
d N j d α = d d α ν j ( n 1 ) C 1 + C 1 ( 1 ) C ν ( 1 ) = C N ( 1 ) ,
which is a constant depending only on n. By induction, we can see that the s -derivative of N j ( α ) with respect to α is bounded by a constant C N ( s ) > 0 depending only on n and s . □
In the next Corollary, we want to show that
u ˜ j = φ j ( k ) u j ( t ) k = 1 n T q j ( K )
is differentiable with respect to α and its derivative is bounded by a constant independent of j and α .
Corollary 7. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 . For all j = 1 , 2 , m , let u ˜ j ( α ) as in (24). Then u ˜ j ( α ) is C s , where s = min { r 2 , r } , and there exist constants C u ˜ ( s ) > 0 , which are independent of j and α such that
d o r d e r = s u ˜ j d α o r d e r = s C u ˜ ( s ) .
Proof. 
Let u ˜ j ( α ) be as in (24). From Definition 1 and Theorem 5, u ˜ is C min { r 2 , r } with respect to α and its first derivative is
d u ˜ j d α = d d α φ j ( k ) u j ( t ) k = 1 n = i = 1 n 1 2 φ j ( k ) u j ( i ) u j ( t ) u j ( i ) α + 2 φ j ( k ) u j ( t ) α k = 1 n .
This is bounded by a constant C u ˜ ( 1 ) > 0 as follows
d u ˜ j ( α ) d α n ( n 1 ) C φ ( 2 , 0 ) C u ( 1 ) + C φ ( 1 , 1 ) = C u ˜ ( 1 ) .
It is clear that this constant depends only on n. And the second derivative of u ˜ j ( α ) with respect to α is
d o r d e r = 2 u ˜ j d α o r d e r = 2 = ( p = 1 n 1 i = 1 n 1 3 φ j ( k ) u j ( p ) u j ( i ) u j ( t ) u j ( i ) α u j ( p ) α + i = 1 n 1 2 φ j ( k ) u j ( t ) α u j ( i ) α + i = 1 n 1 2 φ j ( k ) u j ( i ) u j ( t ) 2 u j ( i ) α 2 p = 1 n 1 3 φ j ( k ) u j ( p ) u j ( t ) α u j ( p ) α + 2 φ j ( k ) u j ( t ) α ) k = 1 n .
Again from Definition 1 and Theorem 5, there exists a constant C u ˜ ( 2 ) > 0 such that
d o r d e r = 2 u ˜ j ( α ) d α o r d e r = 2 n ( n 1 ) ( n 1 ) C φ ( 3 , 0 ) ( C u ( 1 ) ) 2 + 2 C φ ( 2 , 1 ) C u ( 1 ) + C φ ( 2 , 0 ) C u ( 2 ) + C φ ( 1 , 2 ) = C u ˜ ( 2 ) .
This constant depends only on n. By induction, we can see that the s -th derivative of u ˜ j with respect to α is bounded by a constant that depends only on s and n. This proves the statement. □
The next Corollary follows from Corollaries 6 and 7.
Corollary 8. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 . Let κ j ( α ) be the normal curvature of the boundary K ξ j at the point q j ( α ) in the direction u j T q j ( K j ) , which is given by κ j ( α ) = N j ( u ˜ j ) , u j ˜ . Then κ j ( α ) is C s , where s = min { r 3 , r } and there exist constants C κ ( s ) > 0 which depend only on s and n such that
| d o r d e r = s κ j d α o r d e r = s | C κ ( s ) .
Corollary 9. 
Let K ( α ) be a C r , r billiard deformation in R n with r 3 , r 1 . Let U j and U j 1 be as in (6) and (7). Then, U j and U j 1 are at least C s , where s = min { r 2 , r } and there exist constants C U ( s ) > 0 depending only on n and s such that
d o r d e r = s U j d α o r d e r = s C n ( s ) a n d d o r d e r = s U j 1 d α o r d e r = s C n ( s ) .
Proof. 
From (6) and (7), we can write
U j = I ν j v j 1 cos ϕ j and U j 1 = I ν j v j + 1 cos ϕ j ,
where the operators U j 1 and U j depend on α . From Corollaries 3–5, U j and U j 1 are C min { r 2 , r } . Moreover, we can see that their first derivatives are bounded by the same constant C n ( 1 ) such that,
C n ( 1 ) = C ν ( 1 ) + C v ( 1 ) cos ϕ max ,
where cos ϕ max [ 0 , π / 2 ) . C n ( 1 ) is only depending on n. And the second derivatives of U j and U j 1 with respect to α are bounded by C n ( 2 ) such that
C n ( 2 ) = 1 ( cos ϕ max ) 4 C ϕ ( 1 ) 4 C ν ( 1 ) + 4 C v ( 1 ) + C ϕ ( 1 ) + C ν ( 2 ) + C ν ( 1 ) C v ( 1 ) + C v ( 2 ) + C ϕ ( 2 ) .
This constant depends only on n. By induction, we can see that the s -th derivatives of U j and U j 1 are bounded by a constant C n ( s ) > 0 that depends only on n and s . This proves the statement. □
Corollary 10. 
Let K ( α ) be a C r , r billiard deformation in R n with r 3 , r 1 . Let V j and V j be as in (8) and (9). Then, V j and V j are at least C s , where s = min { r 2 , r } and there exist constants C n ( s ) > 0 depending only on n and s such that
d o r d e r = s V j d α o r d e r = s C n ( s ) a n d d o r d e r = s V j d α o r d e r = s C n ( s ) .
Proof. 
From (8) and (9), we can write
V = I v j ν j cos ϕ j and V = I ν j v j cos ϕ j .
Then, the rest of the proof is the same as the proof of Corollary 9. □
The next Corollary follows from Corollaries 5, 6 and 10.
Corollary 11. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 . Let Θ j ( α ) = cos ϕ j V j N j V j . Then, Θ j ( α ) is C s , where s = min { r 3 , r } and there exist constants C Θ ( s ) > 0 depending only on s and n such that
d o r d e r = s Θ j d α o r d e r = s C Θ ( s ) .
In the next Corollary, we want to show that w ^ j ( α ) T q j ( X j ( α ) ) , such that
w ^ j = w j / w j and
w j = u ˜ j u ˜ j , v j cos ϕ j ν j ,
where u ˜ j is as in (24).
Corollary 12. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 , n 3 . For all j = 1 , 2 , m , let w j ( α ) be as in (31). Then, w ^ j ( α ) is C s , where s = min { r 2 , r } , and there exist constants C w ^ ( s ) > 0 , which are independent of j and α such that
d o r d e r = s w ^ j d α o r d e r = s C w ^ ( s ) .
Proof. 
Let w j ( α ) J j as in (31), where J j is the hyperplane to the convex front X ˜ j at the point q j ( α ) . From Corollaries 3, 4 and 7, w ^ j is at least C min { r 2 , r } with respect to α , and there exists a constant C w ^ ( 1 ) > 0 such that
d w ^ j d α 2 ( 1 + cos ϕ + ) cos 2 ϕ + 2 C u ˜ ( 1 ) + C v ( 1 ) + C ν ( 1 ) + C ϕ ( 1 ) = C w ^ ( 1 ) .
This constant depends only on n. By induction, we can see that the s -th derivative of w ^ j with respect to α is bounded by a constant that depends only on s and n. This proves the statement. □
Corollary 13. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 . Then for all j = 1 , 2 , , m , K j ( α ) is C s , where s = min { r 3 , r } and there exists a constant C K ( s ) > 0 , which depends on s and n such that
d o r d e r = s K j ( α ) d α o r d e r = s C K ( s ) .
Proof. 
From (16), recall that
K j ( α ) = U j 1 K j U j + 2 Θ j ,
where K j ( α ) = K j 1 ( I + d j 1 K j 1 ) 1 . Given that every term in the right-hand side of the equation is differentiable with respect to α , and we can derive K j 1 in terms of K 0 on the right-hand side, where K 0 is independent of α , we can conclude that K j is differentiable with respect to α . Specifically, K j is at least C min { r 3 , r } , where r 4 and r 1 . This result follows from the Corollaries 11 and 9.
Next, we want to show that the first derivative of K j with respect to α is bounded by a constant independent of α and j. First, we have
d K j d α = d K j 1 d α ( I + d j 1 K j 1 ) 1 ( I + d j 1 K j 1 ) 1 ( d j 1 d K j 1 d α + d d j 1 d α K j 1 ) ( I + d j 1 K j 1 ) 1 = I d j 1 K j 1 ( I + d j 1 K j 1 ) 1 d K j 1 d α ( I + d j 1 K j 1 ) 1 ( I + d j 1 K j 1 ) 1 d d j 1 d α K j 1 2 ( I + d j 1 K j 1 ) 1 = ( I + d j 1 K j 1 ) 1 d K j 1 d α ( I + d j 1 K j 1 ) 1 ( I + d j 1 K j 1 ) 1 d d j 1 d α K j 1 2 ( I + d j 1 K j 1 ) 1 = D j 1 d K j 1 d α D j 1 + E j 1
where D j 1 ( α ) = ( I + d j 1 K j 1 ) 1 and E j 1 = d d j 1 d α D j 1 K j 1 2 D j 1 . Second, we derive
d K j d α = U j 1 d K j 1 d α U j U j 1 d U j 1 d α U j 1 K j 1 U j + U j 1 K j 1 d U j d α + 2 d Θ j d α .
From (32), we obtain
d K j d α = U j 1 D j 1 d K j 1 d α D j 1 U j + U j 1 E j 1 U j d U j 1 1 d α K j 1 + K j 1 d U j 1 d α + 2 d Θ j d α .
Let P ˜ j = U j 1 D j 1 , P j = D j 1 U j , and R j = U j 1 E j 1 U j d U j 1 d α K j 1 + K j 1 d U j d α + 2 d Θ j d α . So, we can write
d K j d α = P ˜ j d K j 1 d α P j + R j .
To obtain the estimation, we have
d K j d α = P ˜ j P ˜ j 1 d K j 2 d α P j 1 + R j 1 P j + R j = P ˜ j P ˜ j 1 P ˜ j 2 d K j 3 d α P j 2 + R j 2 P j + R j 1 P j 1 + R j = P ˜ j P ˜ j 1 P ˜ 1 d K 0 d α P 1 P j 1 P j + P ˜ j R j 1 P j + + P ˜ j P ˜ j 1 P ˜ 2 R 1 P 2 P j 1 P j + R j .
Since K 0 is independent of α , then
d K j d α = P ˜ j R j 1 P j + + P ˜ j P ˜ 2 R 1 P 2 P j + R j .
From (10) and Corollaries 2, 9 and 11, we have P j D max and P ˜ j D max , where D max = 1 1 + d min μ min , and
R j R max = ( D max ) 2 C d ( 1 ) ( η max ) 2 + C w ( 1 ) η max + C w ( 1 ) d min + 2 C Θ ( 1 ) .
And then
d K j d α 1 + ( D max ) 2 + ( D max ) 4 + + ( D max ) 2 ( j 1 ) R max = ( 1 ( D max ) 2 j 1 ( D max ) 2 R max .
For a large j, ( D max ) 2 j converges to 0. Thus, there exists a constant C K ( 1 ) > 0 independent of j and α , such that
d K j d α C K ( 1 ) .
By using the same approach we can see that the s -th derivative of K j is bounded by a constant C K ( s ) > 0 , which depends only on n and s . □
The next Corollary follows from Corollaries 12 and 13.
Corollary 14. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 , and n 3 . Let k j ( α ) be the normal curvature of the unstable manifold X ˜ j at the point q j ( α ) in the direction w ^ j T q j ( X ˜ j ) , which is given by k j ( α ) = K j ( w ^ j ) , w ^ j . Then, k j ( α ) is C s , where s = min { r 3 , r } and there exist constants C k ( s ) > 0 , which depend only on s and n such that
| d o r d e r = s k j d α o r d e r = s | C k ( s ) .
Corollary 15. 
Let K ( α ) be a C r , r billiard deformation in R n with r 4 , r 1 . Let j ( α ) be as defined in (14). Then, j ( α ) is C s , where s = min { r 3 , r } , and there exist constants C ( s ) > 0 , which depend on s and n such that
| d o r d e r = s j d α o r d e r = s | C ( s ) .
Proof. 
Recall that
( 1 + d j ( α ) j ( α ) ) 2 = w ^ j ( α ) + ( d j ( α ) ) K j ( w ^ j ( α ) ) 2 .
From this and Corollaries 2, 12 and 13, j ( α ) is C s , where s = min { r 3 , r } . And its first derivative with respect to α is
d j d α = 1 d j w ^ j + d j K j ( w ^ j ) w ^ j + d j K j ( w ^ j ) ( d d α w ^ j + d d α d j K j ( w ^ j ) + d j d d α K j ( w ^ j ) + d j K d d α w ^ j d d α d j j .
And again by Corollaries 2, 12 and 13, we have
| d j d α | 1 d min C w 1 ( 1 ) + 2 C d ( 1 ) η max + d max C K 1 ( 1 ) .
Then, there exists a C 1 ( 1 ) > 0 depending on n, such that | d j d α | C 1 ( 1 ) . By induction, we can see that | d o r d e r = s j d α o r d e r = s | C j ( s ) , where C j ( s ) is constant depending only on s and n. □

6. Estimation the Largest Lyapunov Exponent for Open Billiard Deformation

Recall from Section 4 that B α : M α M α and R α : M α Σ , for a small α [ 0 , a ] where a R + . According to Theorem 3 there exists a subset A α of Σ with μ ( A α ) = 1 such that
λ 1 ( α ) = lim m 1 m log D x B α m ( w )
for μ -almost x ( α ) M α with R α ( x ) A α , and all w E 1 u ( x ) E 2 u ( x ) E k u ( x ) E 2 u ( x ) E 3 u ( x ) E k u ( x ) . In case α = 0 , we have A 0 with μ ( A 0 ) = 1 , which we showed in Section 3.
The next Lemma shows that for μ -almost x ( α ) M α , we can choose ξ Σ such that the last Formula (33) holds for α = α p and x ( α p ) for all p = 1 , 2 , , and also for α = 0 and x ( 0 ) , which refers to the initial open billiard.
Lemma 1. 
Given an arbitrary sequence
α 1 , α 2 , , α p ,
of elements of [ 0 , b ] , for μ-almost all ξ Σ the formula (33) is valid for α = α p and x = R α 1 ( ξ ) for all p = 1 , 2 , and also for α = 0 and x = R 1 ( ξ ) .
Proof. 
See [8]. □
From the formula for the largest Lyapunov exponent (11), we can write the Lyapunov exponents for K ( α ) and K ( 0 ) as follows:
λ 1 ( α ) = lim m 1 m j = 1 m log 1 + d j ( α ) j ( α ) = lim m λ 1 ( m ) ( α ) , λ 1 ( 0 ) = lim m 1 m j = 1 m log 1 + d j ( 0 ) j ( 0 ) = lim m λ 1 ( m ) ( 0 ) ,
where
λ 1 ( m ) ( α ) = 1 m j = 1 m log 1 + d j ( α ) j ( α ) and λ 1 ( m ) ( 0 ) = 1 m j = 1 m log 1 + d j ( 0 ) j ( 0 ) .

7. The Continuity of the Largest Lyapunov Exponent for Non-Planar Billiard Deformation

This section provides a rigorous proof that the largest Lyapunov exponent of the billiard deformation in R n , where n 3 , is a continuous function at every perturbation parameter α [ 0 , b ] . Our proof closely follows the argument used in the proof of the continuity of the largest Lyapunov exponent for planar billiard deformation [8], as well as the results established in Section 5.1 and previous sections.
Here we prove Theorem 1:
Proof of Theorem 1. 
We will show that λ 1 ( α ) is continuous at α = 0 . From this, the continuity is followed at every α [ 0 , b ] . To prove that, we will show that for any sequence of points α p [ 0 , b ] as in Lemma 1 with α p 0 when p , λ 1 m ( α p ) λ 1 m ( 0 ) .
From (34), we have
λ 1 m ( α p ) λ 1 m ( 0 ) = 1 m i = 1 m log ( 1 + d j ( α p ) j ( α p ) ) log ( 1 + d j ( 0 ) j ( 0 ) ) .
Let f j ( α p ) = log ( 1 + d j ( α p ) j ( α p ) ) . From Corollaries 2 and 15, f j ( α p ) is C 2 with respect to α p , and its first derivative is
d f j ( α p ) d α p = d d j ( α p ) d α p j ( α p ) + d j ( α p ) d j ( α p ) d α p 1 + d j ( α p ) j ( α p ) .
In addition, there exists a constant C f ( 1 ) > 0 independent of j and α p such that,
| d f j ( α p ) d α p | C d ( 1 ) η max + d max C ( 1 ) 1 + d min μ min = C f ( 1 ) .
From this and applying the Mean Value Theorem to (35), we obtain
| λ 1 m ( α p ) λ 1 m ( 0 ) | 1 m i = 1 m α p C f ( 1 ) = α p C f ( 1 ) .
By letting α p approach 0 as p approaches infinity, and then letting m approach infinity, we obtain λ 1 ( α p ) approaching λ 1 ( 0 ) , for any sequence α 1 > α 2 > > α p > 0 as p in [ 0 , b ] as in Lemma 1. This proves the statement. □

8. The Differentiability of the Largest Lyapunov Exponent for the Non-Planar Billiard Deformation

In this section, we present a proof of the differentiability of the largest Lyapunov exponent λ 1 for the billiard deformation in R n , with n 3 , with respect to a small perturbation α . This proof is based on the findings in the preceding sections and resembles the proof in [8] for the case of planar billiard deformation.
Here is the proof of Theorem 2:
Proof of Theorem 2. 
We prove differentiability at α = 0 , and this implies differentiability for any α [ 0 , b ] . To prove the differentiability at α = 0 , we have to show that there exists
lim α 0 λ 1 ( α ) λ 1 ( 0 ) α .
Equivalently, there exists a number F such that
lim p λ 1 ( α p ) λ 1 ( 0 ) α p = F ,
for any sequence α 1 > α 2 > > α p > 0 as p in [ 0 , b ] as in Lemma 1. By using this Lemma and the expressions of λ 1 m ( α ) for α = α p and λ 1 m ( 0 ) for α = 0 in (34), we have λ 1 ( m ) ( α p ) λ 1 ( α p ) and λ 1 ( m ) ( 0 ) λ 1 ( 0 ) when m . Let f j ( α p ) = log 1 + d j ( α p ) j ( α p ) and f j ( 0 ) = log 1 + d j ( 0 ) j ( 0 ) . And Taylor–Lagrange formula gives
f j ( α p ) = f j ( 0 ) + α p f ˙ j ( 0 ) + α p 2 2 f ¨ j ( r j ( α p ) )
for some r j ( α p ) [ 0 , α p ] . Then,
f j ( α p ) f j ( 0 ) α p f ˙ j ( 0 ) = α p 2 f ¨ j ( r j ( α p ) ) .
After multiplying by m 1 and summing up for j = 1 to m, we have
1 m j = 1 m f j ( α p ) f j ( 0 ) α p 1 m j = 1 m f ˙ j ( 0 ) = 1 m j = 1 m α p 2 f ¨ j ( r j ( α p ) ) .
Thus,
λ 1 ( m ) ( α p ) λ 1 ( m ) ( 0 ) α p F m = 1 m j = 1 m α p 2 f ¨ j ( r j ( α p ) ) ,
where F m = 1 m j = 1 m f ˙ j ( 0 ) . Next, we showed the first derivative of f j ( α p ) in (36). As f j ( α p ) is C 2 , we can now proceed to calculate the second derivative as follows:
f ¨ j ( α p ) = d ¨ j ( α p ) j ( α p ) + 2 d ˙ j ( α p ) ˙ j ( α p ) + d j ( α p ) ¨ j ( α p ) 1 + d j ( α p ) j ( α p ) 1 + d j ( α p ) j ( α p ) 2 d ˙ j ( α p ) ( α p ) + d j ( α p ) ˙ j ( α p ) 2 1 + d j ( α p ) j ( α p ) 2 .
According to Corollaries 2 and 15, there exists a constant C f ( 2 ) > 0 that does not depend on j and α p such that,
| f ¨ j ( α p ) | C d ( 2 ) η max + 2 C d ( 1 ) C ( 1 ) + d max C ( 2 ) 1 + d max η max 1 + d min μ min 2 + C d ( 1 ) η max + d max C ( 1 ) 2 1 + d min μ min 2 = C f ( 2 ) .
Therefore, | f ¨ j ( r j ( α p ) ) | C f ( 2 ) , implying that
| λ 1 ( m ) ( α p ) λ 1 ( m ) ( 0 ) α p F m | 1 m j = 1 m α p 2 | f ¨ j ( t j ( α p ) ) | C f ( 2 ) 2 α p .
From (37), we have | F m | 1 m j = 1 m | f ˙ j ( 0 ) | C f ( 1 ) . So, by the Bolzano–Weierstrass Theorem, there exists a convergent subsequence { F m h } of { F m } such that lim h F m h = F . Then
| λ 1 ( m h ) ( α p ) λ 1 ( m h ) ( 0 ) α p F m h | C f ( 2 ) 2 α p .
By letting h , and letting lim p α p = 0 , we obtain for every sequence α 1 > α 2 > > α p > 0 as p in [ 0 , b ] as in Lemma 1,
lim p λ 1 ( α p ) λ 1 ( 0 ) α p = F .
Thus, there exists F = lim m 1 m j = 1 m f ˙ j ( 0 ) . This holds for any subsequence { m h } . Therefore, for each subsequence, we can conclude that F m h approaches F. Therefore, F m also converges to F. This means that
lim α 0 λ 1 ( α ) λ 1 ( 0 ) α = F ,
exists. Consequently, λ 1 is differentiable at α = 0 and λ ˙ 1 ( 0 ) = F . This proves the statement. □

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

Prof. Luchezar Stoyanov is gratefully acknowledged by the author for his insightful advice, helpful criticism, and assistance. The author would also like to thank Dr. Des Hill and the referees for their careful reading, remarks, and suggestions. This work was supported by a scholarship from Najran University, Saudi Arabia.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Ikawa, M. Decay of solutions of the wave equation in the exterior of several strictly convex bodies. Ann. Inst. Fourier 1988, 38, 113–146. [Google Scholar] [CrossRef]
  2. Barreira, L.; Pesin, Y. Lyapunov Exponents and Smooth Ergodic Theory; American Mathematical Society: Providence, RI, USA, 2002. [Google Scholar]
  3. Chernov, N. Entropy, Lyapunov exponents, and mean free path for billiards. J. Stat. Phys. 1997, 88, 1–29. [Google Scholar] [CrossRef]
  4. Chernov, N.; Markarian, R. Introduction to the Ergodic Theory of Chaotic Billiards; IMCA: Lima, Peru, 2001. [Google Scholar]
  5. Markarian, R. Billiards with Pesin Region of Measure one. Commun. Math. Phys. 1988, 118, 87–97. [Google Scholar] [CrossRef]
  6. Markarian, R. New ergodic Billiards: Exact results. Nonlinearity 1993, 6, 819–841. [Google Scholar] [CrossRef]
  7. Wojtkowski, M.P. Principles for the design of billiards with nonvanishing Lyapunov exponents. Commun. Math. Phys. 1986, 105, 391–414. [Google Scholar] [CrossRef]
  8. Al Dowais, A. Differentiability of the largest Lyapunov exponent for planar open billiards. Dyn. Syst. 2023, 38, 556–575. [Google Scholar] [CrossRef]
  9. Wright, P. Dimensional Characteristics of the Non-Wandering Sets of Open Billiards. Ph.D. Dissertation, The University of Western Australia, Perth, Australia, 2014. [Google Scholar]
  10. Lopes, A.; Markarian, R. Open billiards: Invariant and conditionally invariant probabilities on Cantor sets. SIAM J. Appl. Math. 1996, 56, 651–680. [Google Scholar] [CrossRef]
  11. Stoyanov, L. Exponential instability and entropy for a class of dispersing billiards. Ergod. Theory Dyn. Syst. 1999, 19, 201–226. [Google Scholar] [CrossRef]
  12. Petkov, V.; Stoyanov, L. Geometry of the Generalized Geodesic Flow and Inverse Spectral Problems, 2nd ed.; John Wiley & Sons: Chichester, UK, 2017. [Google Scholar]
  13. Sinai, Y. Dynamical systems with elastic reflections. Russ. Math. Surv. 1970, 25, 137–190. [Google Scholar] [CrossRef]
  14. Sinai, Y. Development of Krylov’s Ideas, An Addendum to: N.S.Krylov. Works on the Foundations of Statistical Physics; Princeton Univ. Press: Princeton, NJ, USA, 1979; pp. 239–281. [Google Scholar]
  15. Morita, T. The symbolic representation of billiards without boundary condition. Trans. Am. Math. Soc. 1991, 325, 819–828. [Google Scholar] [CrossRef]
  16. Katok, A.; Strelcyn, J.M. Invariant Manifolds, Entropy and Billiards; Smooth Maps with Singularities; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 1986; Volume 1222. [Google Scholar]
  17. Viana, M. Lectures on Lyapunov Exponents; Cambridge Studies in Adv. Math.; Cambridge Univ. Press: Cambridge, UK, 2014; Volume 145. [Google Scholar]
  18. Stoyanov, L. Non-integrability of open billiard flows and Dolgopyat-type estimates. Ergod. Theory Dyn. Syst. 2012, 32, 295–313. [Google Scholar] [CrossRef]
  19. Chernov, N. Structure of transversal leaves in multidimensional semidispersing billiards. Funct. Anal. Appl. 1982, 16, 35–46. [Google Scholar] [CrossRef]
  20. Bálint, P.; Chernov, N.; Szász, D.; Tóth, I.P. Geometry of multi-dimensional dispersing billiards. Asterisque 2003, 286, 119–150. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al Dowais, A. Differentiability of the Largest Lyapunov Exponent for Non-Planar Open Billiards. Mathematics 2023, 11, 4633. https://doi.org/10.3390/math11224633

AMA Style

Al Dowais A. Differentiability of the Largest Lyapunov Exponent for Non-Planar Open Billiards. Mathematics. 2023; 11(22):4633. https://doi.org/10.3390/math11224633

Chicago/Turabian Style

Al Dowais, Amal. 2023. "Differentiability of the Largest Lyapunov Exponent for Non-Planar Open Billiards" Mathematics 11, no. 22: 4633. https://doi.org/10.3390/math11224633

APA Style

Al Dowais, A. (2023). Differentiability of the Largest Lyapunov Exponent for Non-Planar Open Billiards. Mathematics, 11(22), 4633. https://doi.org/10.3390/math11224633

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop