{'fraction_non_alphanumeric': 0.06492870965589341, 'fraction_numerical': 0.01637119914585048, 'mean_word_length': 3.337900424546507, 'pattern_counts': {'":': 0, '<': 29, '<?xml version=': 0, '>': 51, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 20, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
SHRINKING TARGETS FOR THE GEODESIC FLOW ON GEOMETRICALLY FINITE HYPERBOLIC MANIFOLDS

Let M be a geometrically finite hyperbolic manifold. We present a very general theorem on the shrinking target problem for the geodesic flow, using its exponential mixing. This includes a strengthening of Sullivan’s logarithm law for the excursion rate of the geodesic flow. More generally, we prove logarithm laws for the first hitting time for shrinking cusp neighborhoods, shrinking tubular neighborhoods of a closed geodesic, and shrinking metric balls, as well as give quantitative estimates for the time a generic geodesic spends in such shrinking targets.


INTRODUCTION
Let M be a complete hyperbolic manifold of dimension n ≥ 2. Denote by G t the geodesic flow on the unit tangent bundle T 1 (M ). If M is of finite volume, but non-compact, Sullivan [27] showed in 1982 the following logarithm law for the rate of the excursion of the geodesic flow: for any o ∈ M , and for almost all x ∈ T 1 (M ), where d (G t (x), o) is the hyperbolic distance between the basepoint of G t (x) and o.
This result can be viewed as a special case of the so-called shrinking target problem for the geodesic flow, which asks the behavior of a generic geodesic ray with respect to a given sequence of shrinking subsets. Indeed, if we consider the family of shrinking cuspidal neighborhoods h t := {z ∈ M : d (o, z) > t }, t > 1, then (1.1) is equivalent to the following logarithm law for the first hitting time: for almost all x, where τ h t (x) := inf{s > 0 : G s (x) ∈ h t }.
In this paper, we investigate shrinking target problems for the geodesic flow on a geometrically finite hyperbolic manifold M , and prove results which are far reaching strengthening and generalizations of (1.2), and hence of (1.1).
Let H n denote the n-dimensional hyperbolic space and let G := Isom + (H n ) be the group of all orientation preserving isometries. We may present a complete hyperbolic manifold M as the quotient Γ\H n where Γ is a torsion-free discrete subgroup of G. We assume that Γ is Zariski dense and geometrically finite in the whole paper. Denote by Λ ⊂ ∂H n the limit set of Γ and by 0 < δ ≤ n − 1 the critical exponent of Γ. The maximal entropy of the geodesic flow on T 1 (M ) is given by δ, and there exists a unique ergodic probability measure of maximal entropy, called the Bowen-Margulis-Sullivan measure, which we denote by m. The support of m is precisely the non-wandering set for the geodesic flow, and hence the shrinking target problem in this setting is interesting only for those shrinking subsets in the support of m and for m-almost all points. Now since G t is ergodic for m, the Birkhoff ergodic theorem says that for a given Borel subset B ⊂ T 1 (M ), we have the following: for m-almost all x ∈ T 1 (M ), for an infinite sequence of times t tending to ∞ or for all sufficiently large t 1? 3. (Quantitative estimates) How fast can B t shrink so that 1 for an infinite sequence of times t tending to ∞, or for all sufficiently large t 1?
In order to address the above questions, we need to impose certain regularity conditions on the shrinking targets. Let K < G be a maximal compact subgroup and identify M with Γ\G/K . There exists a one parameter diagonalizable subgroup A = {a t } so that if M denotes the centralizer of A in K , then the unit tangent bundle T 1 (M ) can be identified with Γ\G/M in the way that the geodesic flow G t on T 1 (M ) corresponds to the right translation action of a t 1 The notation f t g t means that for all t > 1, f t ≤ c g t for some absolute constant c > 0, and we write f t g t if f t g t and g t f t . We sometimes indicate the dependence of the implied constant in subscripts. on Γ\G/M . We fix dim(M ) and the Sobolev norm S = S ∞, on C ∞ (Γ\G) given by where the sum is taken over all monomials in a fixed basis of Lie(G) of order at most . A family of shrinking targets in T 1 (M ) means a collection B = {B t ⊂ T 1 (M ) : t > 1} such that m(B t ) > 0, B t ⊃ B s for s > t , and lim t →∞ m(B t ) = 0. A family {B t } of shrinking targets is said to be inner regular (resp. outer regular) if there exist α > 0 and a family of functions Ψ − t ∈ C ∞ (T 1 (M )) (resp. Ψ + t ∈ C ∞ (T 1 (M ))) such that where the implied constants are independent of t . A family {B t } is said to be regular if it is both inner and outer regular. We note that this regularity condition is rather mild, and is satisfied by most families of naturally occurring shrinking targets. Such examples include shrinking cusp neighborhoods, shrinking tubular neighborhoods of a closed geodesic and shrinking metric balls, as will be shown later.
In the rest of the introduction, we assume that B = {B t : t 1} is a family of shrinking targets in T 1 (M ).
1.1. Logarithm laws. For discrete time dynamical systems, it is expected that the first hitting time would be inversely proportional to the measure of the shrinking target; it is indeed the case for the discretized geodesic flow. For the continuous geodesic flow, it turns out that it is inversely proportional to the measure of a thickened setB t := ∪ |s|<1/2 G s (B t ): THEOREM 1.1.

REMARK 1.2.
When | log m(B t )| | log m(B t )|, the first hitting time for the discrete flow {G n : n ∈ N} behaves in the same way for the continuous flow. This is indeed the case for shrinking cusp neighborhoods or tubular neighborhoods of a closed geodesic. However, there are also cases when | log m(B t )| is much larger than | log m(B t )|, such as the case of shrinking metric balls.
We note that logarithm laws for the first hitting time were studied for certain families of shrinking targets in many examples of discrete time dynamical systems with fast mixing, see, e.g., [6,7,8].
1.2. Shrinking rate threshold. In order to ensure that a generic orbit G s (x) hits B t before time t for an infinite sequence of t tending to ∞, the easy half of the Borel-Canteli lemma implies that it is necessary to have k m(B k ) = ∞, from which lim sup t →∞ log 2 (t )t · m(B t ) = ∞ follows. The first part of the following theorem says that this condition is also sufficient, up to logarithmic factors. The second part says that a generic orbit G s (x) hits B t before time t , for all sufficiently large t , under a slightly stronger assumption on the rate of shrinking (see Theorem 4.11).

Quantitative estimates.
In order to answer a more refined question regarding the amount of time that a generic geodesic ray spends in a shrinking target, we require our family of targets to be regular and their measures do not change too fast in the sense that m(B t ) m(B 2t ).
With these additional regularity assumptions, we have the following (see Theorem 4.8 below for a more general result).
We observe that unlike Theorems 1.1 and 1.3, the amount of time that the geodesic flow spends in the targets is governed by the measure of the original targets rather than by their thickenings.   [27], or spherical balls [21] (or general spherical targets if one considers discrete time dynamics [13]). There are no known sharp results regarding Theorem 1.4 (2). We refer to [15] where this kind of problem is studied for systems with almost perfect mixing. 3. All the results described above still hold as stated if we replace the unit tangent bundle T 1 (M ) with the frame bundle Γ\G, provided δ > n−2. We note if M contains a co-dimension one properly immersed totally geodesic submanifold of finite volume, then δ > n − 2, so this stronger condition still holds in many examples.
For some concrete applications of these results, we discuss three families of shrinking targets to which our theorems apply. In order to define these families, we fix a left G-invariant and right K -invariant metric d on G which descends to the hyperbolic metric on H n = G/K . This metric then naturally defines a distance function, dist(·, ·) on T 1 (M ) = Γ\G/M .

Cusp excursion.
The convex core of M is defined by core(M ) = Γ\ hull(Λ), where hull(Λ) defines the convex hull of the limit set Λ. As M is geometrically finite, there are finitely many disjoint cuspidal regions whose complement in core(M ) is a compact submanifold. Let h i , 1 ≤ i ≤ k, denote the pre-images in T 1 (M ) of these cuspidal regions under the base point projection π : T 1 (M ) → M .
For each i , we denote by κ i the rank of h i , that is, the rank of the maximal free abelian subgroup of the stabilizer Stab Γ (h i ). It is known that κ i < 2δ.
For each i and t > 1, consider the following cusp neighborhood (1.5) For each i , we show that the shrinking family {h i ,t : t > 1} is regular and that (1.6) (see Section 5.1). Applying our results to this family, we get the following: REMARK 1.7. As mentioned before, it is not hard to show that where h t = 1≤i ≤k h i ,t . Stratmann and Velani showed that (1.7) is equal to 2δ − max i κ i [28], and hence extended Sullivan's logarithm law (1.1) to geometrically finite manifolds. Theorem 1.6(1) presents a stronger version, as we consider excursion to individual cusps as well as obtain an actual limit rather than lim inf.
For the sake of a concrete application, we give a reformulation of Theorem 1.6(1) in the case of Apollonian manifolds. An Apollonian gasket P = C i is a countable union of circles obtained by repeatedly inscribing circles into the triangular interstices of four mutually tangent circles with disjoint interiors in the complex plane (where lines are considered as circles). The symmetry group {g ∈ PSL 2 (C) : g (P ) = P } is a discrete subgroup of PSL 2 (C) which acts on C by Möbius transformations and its torsion-free subgroup of finite index is called an Apollonian group, which we denote by Γ. Via the Poincaré extension theorem, we can identify PSL 2 (C) with Isom + (H 3 ) for the upper-half space model H 3 of the hyperbolic space. The quotient manifold Γ\H 3 is called an Apollonian manifold, which is known to be geometrically finite with all cusps having rank one. Its limit set is equal to the closure P , and supports a locally finite Hausdorff measure H of dimension δ = 1.30568(8) [11].
Fix a tangent point ξ = C i ∩ C j for i = j and consider a sufficiently small Euclidean ball B in H 3 based at ξ, so that B = Γ(B ) is a disjoint collection of Euclidean balls.
Fix o ∈ H 3 outside of B , let B (t ) ⊂ B be the Euclidean ball based at ξ and The following is a consequence of Theorem 1.6: COROLLARY 1.8. Let P be an Apollonian gasket. For H -almost all initial direction v toward P , where v s denotes the base point of G s (v).
1.5. Tubular neighborhoods. Another natural family of shrinking targets is given by tubular neighborhoods of a closed geodesic. For a closed geodesic C ⊂ T 1 (M ) and > 0, we consider the -tubular neighborhood of C : The family {C 1/t : t > 1} forms a family of shrinking neighborhoods of C . We show that {C 1/t : t > 1} is a regular family with m(C 1/t ) m(C 1/t ) t −2δ . Applying our results to this family of shrinking targets gives the following result on the amount of time a generic geodesic spirals near a fixed closed geodesic (cf. [10, Theorem 1.1] for a similar result in a negatively curved compact manifold). THEOREM 1.9. Let C ⊂ T 1 (M ) be a closed geodesic. Then for m-a.e. x ∈ T 1 (M ), we have the following: 2. For any 0 < η < 1 2δ and for all t > 1, REMARK 1.10. Since for any x ∈ T 1 (M ) we have that which was previously shown in [4,Theorem 4] to hold for the special case of convex co-compact hyperbolic surfaces.
x ∈ T 1 (M ), When M has cusps, the situation is more complicated as m(B t (x 0 )) can fluctuate, with the fluctuation depending on x 0 (or more precisely on the cusp excursions of the geodesic emanating from x 0 ∈ T 1 (M )). Combining our previous results on cusp excursions, we can show the following THEOREM 1.12. Suppose that M has cusps.

1.
For m-a.e. x 0 ∈ T 1 (M ), and for m-a.e. x ∈ T 1 (M ), 2. For any pair of distinct cusps of ranks κ 1 , κ 2 , we can find x 0 ∈ T 1 (M ) such that for m-a.e. x ∈ T 1 (M ), REMARK 1. 13. We note that if M has finite volume, then δ = n − 1 and m(B t (x 0 )) t −(2n−2) . Hence, in this case, the same arguments imply that for m-a.e. x ∈ T 1 (M ), we have lim t →∞ log τ B t (x 0 ) (x) log t = 2(n − 1). We note that here the shrinking targets are in T 1 (M ), unlike the results of [21] which considered shrinking balls inside M , in which case the limit is n−1 (see also [17], for related result for the discrete time geodesic flow).
If Ψ is the characteristic function of B , we simply write λ T (B ) instead of λ T (1 B ). The Birkhoff ergodic theorem implies that for a.e. x ∈ X , We note that if we had a rate control in this convergence such as An estimate like (1.10) is too strong to be true for a.e. individual points x. So, instead, we prove its mean-version for all smooth functions Ψ ∈ L 2 (T 1 (M ), m), that is, for some uniform constant C > 0. The regularity conditions imposed on the thickeningsB t of our shrinking targets are precisely so that we could apply (1.12) to smooth functions which approximates 1B t and deduce This effective mean ergodic theorem forB t 's enables us to obtain that for a.e. x, (1.14) for all sufficiently large t . Using that |τ d This is the non-trivial direction of the logarithm law Theorem 1.1; the other direction holds for general shrinking targets in any dynamical system (see, e.g., [14,Lemma 2.2]). Theorems 1.3 and 1.4 are also proved in a similar spirit using the effective mean ergodic theorem.
The use of quantitative mixing of geodesic flow in the shrinking target problem in the homogeneous setting goes back to the work of Kleinbock and Margulis [16], and the idea of using an effective mean ergodic theorem was first introduced in [9] and more explicitly in [13,14], where these ideas were used to prove the analogous results for finite volume hyperbolic manifold. In this paper we follow closely the arguments of [13], in particular, for the discretized geodesic flows; once we establish the relevant result on exponential decay of matrix coefficients as in Theorem 1.14, the arguments are essentially identical, although establishing the needed regularity conditions for explicit examples is technically more difficult in the infinite volume case (see Section 5). On the other hand, our results for continuous flows are new even for the finite volume case. It is worthwhile mentioning that, for spherical shrinking targets, the results of [13] also dealt with unipotent flows having only polynomial decay of matrix coefficients, and established a dynamical Borel-Cantelli Lemma, which is stronger than the logarithm law. While it would be very interesting to extend such results to the infinite volume setting, our current methods do not seem to apply.
Here we will use the following exponential decay of matrix coefficients for geometrically finite hyperbolic manifolds: THEOREM 1.14. There exists η 0 > 0 such that for any Ψ 1 , Ψ 2 ∈ C ∞ (T 1 (M )) with support in one-neighborhood of supp(m), for all t ≥ 1, Moreover, η 0 is explicitly computable when δ > n− 1 2 , depending only on the spectral gap for the Laplacian on L 2 (M ). If Γ is convex cocompact or δ > n − 2, (1.15) with m replaced by its M -invariant lift on Γ\G holds for any Ψ 1 , Ψ 2 ∈ C ∞ (Γ\G).
This theorem was obtained in [22,5] for compactly supported functions under the assumption δ > n−1 2 and in [25] for any convex cocompact Γ (see also [24] for the same result for the frame flow). In order to study shrinking target problem for cusp neighborhoods as described in Theorem 1.6, removing the compact support condition is crucial as we need to study functions that are positive on cusps. We use the quantitative decay of the matrix coefficient of the functions L 2 (Γ\G) with respect to the Haar measure m Haar in [22], and exploit the product structures of m and m Haar to transfer the exponential rate information on the transversal intersections of G t (B (x)) for the flow box B (x), that we get from the behavior of the correlation function with respect to m Haar , to the behavior of the correlation function with respect to m. Here depends on the injectivity radius of x, and as we need to control the exponential rate independent of the injectivity radius for Theorem 1.14, which is required to deal with functions which are not compactly supported, the whole procedure turns out to be technically quite subtle. The remaining cases of geometrically finite manifolds with cusps are proved in a recent work of Li-Pan [20].
After some preliminaries given in Section 2, we devote Section 3 to the proof of Theorem 1.14. With this result in hand, we prove effective mean ergodic theorem in this setting (see Theorem 4.2), and use it in Section 4 to establish results on shrinking target problems for both the discrete and continuous time flow. While the results we obtain for the discrete time flow are essentially optimal, this is not the case for some of the results for continuous time flow. Nevertheless, in Section 4.5, we show how one can obtain optimal results for the continuous flow by translating it into a discrete time flow problem for a thickened target. In Section 5, we deduce Theorems 1.6, 1.9, 1.11 and 1.12 by proving the regularity of the corresponding shrinking sets and by computing their volumes using Sullivan's shadow lemma and the structure of cusps for geometrically finite manifolds.

Notations and conventions.
Let G ∼ = SO(n, 1) o be the group of orientation preserving isometries of H n , and Γ < G a geometrically finite, torsion-free, Zariski dense, discrete subgroup of G. We denote by Λ the limit set of Γ, and by 0 < δ ≤ n − 1 the Hausdorff dimension of Λ, which is equal to the critical exponent of Γ. Let M = Γ\H n . Let K < G be a maximal compact subgroup and identify M with Γ\G/K . There exists a one parameter diagonalizable subgroup A = {a t } so that if M denotes the centralizer of A in K , then the unit tangent bundle T 1 (M ) can be identified with Γ\G/M in the way that the geodesic flow G t on T 1 (M ) corresponds to the right translation action of a t on Γ\G/M . With this identification we can work in the homogeneous space Γ\G and think of subsets and functions on T 1 (M ) and M respectively as M -invariant (resp. K invariant) subsets and functions on Γ\G.
We say that two families {B t } and {A t } of shrinking sets are Lipschitz equivalent and write B t A t , if there are some positive constants c 1 , c 2 such that We fix a left G-invariant and right K -invariant metric d on G which descends to the hyperbolic metric on H n = G/K . This induces a unique metric on G/M which we will also denote by d by abuse of notation. The metric d defines a distance function on T 1 (M ) = Γ\G/M given by dist(Γg , Γh) = inf γ∈Γ d (γg , h).

Invariant measures.
For ξ ∈ ∂H n , let β ξ : H n × H n → R denote the Busemann function for the geodesic flow, defined by with ξ(t ) a unit speed geodesic ray toward ξ. A family of measures {µ x : x ∈ H n } is called a Γ-invariant conformal density of dimension δ µ > 0 on ∂H n , if each µ x is a non-zero finite Borel measure on ∂H n satisfying for any x, y ∈ H n , ξ ∈ ∂H n and γ ∈ Γ, In particular, the Patterson-Sullivan density {ν x } is a Γ-invariant conformal density supported on the limit set Λ of dimension δ and the Lebesgue density {m x } is a G-invariant conformal density of dimension (n − 1) (both are unique up to scalar multiplications).
Let π : T 1 (H n ) → H n be the basepoint projection. For u ∈ T 1 (H n ), we denote by u ± ∈ ∂H n the forward and the backward endpoints of the geodesic determined by u.
In these coordinates, the BMS measure m = m BMS , the Haar measure m Haar , and the Burger-Roblin measure m BR on T 1 (H n ) are given by These measures are all left Γ-invariant, and hence descend to corresponding measures on T 1 (M ) . Using T 1 (H n ) = G/M , we can lift the above measures to right M -invariant measures on Γ\G, which we still denote by m, m Haar and m BR by abuse of notation. The measure m is finite and ergodic with respect to the geodesic flow [27]. We will normalize the Patterson-Sullivan density {ν x } so that m(T 1 (M )) = m(Γ\G) = 1.
Let N = N + and N − denote the expanding and the contracting horospherical subgroups respectively, i.e., N ± = {g ∈ G : a s g a −s → e as s → ±∞}.

Note that
The BMS measure m has a natural foliation corresponding to the decomposition P N = G (modulo a Zariski closed subset) with P = N − AM . Explicitly, for any g ∈ G, we define the PS-measure and the Lebesgue measure on the coset g N , by respectively. We also define the measureν g P on the coset g P by g p). Using the decomposition G = g P N and noting that (g pn) − = (g p) − , we have that for any Ψ ∈ C c (G), Finally, for x = [g ] ∈ Γ\G and > 0 smaller than the injectivity radius at x, we denote by d µ PS xN and d ν xP the measures induced by dμ PS g N and dν g P on xN and xP respectively.

Cusp decomposition.
Let X 0 be the pre-image of the convex core of M under the base point projection map π : Γ\G → Γ\G/K = M and let X be the unit neighborhood of X 0 . Then Ω ⊆ X 0 ⊆ X and since M is geometrically finite, X has finite Haar-measure. When M is convex cocompact, X is compact, and otherwise it can be decomposed into a compact part and finitely many cusp neighborhoods, as we describe below.
Let Λ p ⊂ Λ denote the set of parabolic fixed points (i.e., points fixed by some parabolic element of Γ). Since Γ is geometrically finite, Λ p consists of finitely many Γ-orbits represented by {ξ 1 , . . . , ξ k } which are called cusps of M . A cuspidal neighborhood of ξ i ∈ Λ p is a set of the form where H ξ ⊆ H n is some fixed horoball tangent to ξ such that γH ξ ∩ H ξ = if and only if γ fixes ξ. For each i , the stabilizer Stab Γ (ξ i ) is a free abelian subgroup and we denote its rank by κ i . We set κ max := max κ i and κ min := min κ i . Note that 2δ > k max (see [3,Lemma 3.5]).
For x ∈ Γ\G, we denote by r x the injectivity radius at x. For all sufficiently small > 0, let X ( ) = {x ∈ X : r x < }, so that is compact, and the family X ( ) with < 0 forms a shrinking family of cusp neighborhoods.
More explicitly, we show in Section 5.1 that for all sufficiently small > 0, where the sum is taken over all monomials D of order at most l in the basis elements, and Ψ Haar p denotes the L p (Γ\G, m Haar )-norm of Ψ. While this norm depends on the choice of basis, changing the basis will only change the norm by some bounded factor.
We will mostly use the norms S ∞,l , which we will denote by S l to simplify notation. Since supp(m) ⊂ X , it is sufficient for our purpose to consider functions supported on X , and since X has finite Haar measure we can, and will, use the bound where the implied constant is independent of Ψ ∈ C ∞ (X ).

DECAY OF MATRIX COEFFICIENTS
A crucial ingredient in our proof is the exponential mixing of the geodesic flow with respect to the BMS-measure. We use the inner product notation: By the remarks following Theorem 1.14, the following theorem is the only missing part of it, given the works [24] and [20].
Then there exist an explicit η 0 > 0 (depending only on the spectral gap of L 2 (M )) and l ∈ N, such that for any bounded In the rest of this section, we assume δ > (n − 1)/2. Theorem 3.1 with an explicit η 0 depending only on the spectral gap of L 2 (M ) is then proved in ([22, Theorem 6.16], [5]) under the assumption that the test functions are compactly supported. In order to complete the proof of the theorem we need to remove the assumption on the support of the test functions.
To do this, we will approximate Ψ as the sum Ψ + (Ψ − Ψ ) where Ψ is a smooth function supported on Y ( ), and similarly for Φ. In view of (2.7), the main term will be reduced to 〈a t Ψ , Φ 〉, for which the result follows from [22,Theorem 6.16]. However, since the dependence on the supports of Ψ and Φ was not made explicit in terms of in [22], we need to redo their arguments while keeping track of the dependence on as well as on all implied constants along the proof.

Control of BR measures.
Since m BR (Γ\G) = ∞ when Γ < G is not a lattice, and some of the implied constants in [22,Thm. 6.16] depend on m BR (supp(Ψ)), we need the following result to control the dependence on these measures. Proof. Recall that by [26] and [19], there exists a positive eigenfunction φ 0 ∈ C ∞ (Γ\G) K for the Laplace operator such that Under the assumption δ > n−1 2 , we have φ 0 Since X is K -invariant with m Haar (X ) < ∞, the following is a consequence of Lemma 3.2:

Test function supported on small balls.
For a subset S ⊆ G and > 0, S denotes the -neighborhood of e in S, that is, S = {g ∈ S : d (g , e) ≤ }. Set B := P N ; and note that G B for all sufficiently small > 0. In this subsection, we will prove the following.
). There exist l ∈ N depending only on dim(G) and η > 0 (depending only on the spectral gap of Γ) such that for any ∈ (0, 1) small and any where the implied constant is absolute.
Proof. Fix Φ, Ψ ∈ C ∞ (xB ). In the case when n−1 2 < δ ≤ n − 2, we assume that Fix a small 0 < 0 < 2 and consider the functions Ψ ± 0 on Γ\G defined by and let We then have that , where we used that m(X ) < ∞. We will also use the notation and similarly get that µ PS . Now by ([22, Lemma 6.2], [5]), there exists some absolute constant c > 0, such that the integral is bounded from above and below, respectively, by where P x (t ) is the finite set defined by Moreover, by the proof of [22,Theorem 6.7], there are positive constants η > 0 (depending only on the spectral gap of Γ) and α > 0 such that A BR Ψ := S ∞,1 (Ψ)m BR (supp(Ψ)) S ∞,1 (Ψ), (by Lemma 3.2), and . Notice that the injectivity radii of the supports of φ and ψ are at least and since we chose 0 2 much smaller, all the implied constants are absolute and independent of and 0 . Combining these results and estimating Since all implied constants are independent of 0 , taking the limit as 0 → 0 gives Now, integrating over xP , and noting that xP µ PS the integral of the first remainder term is bounded by O(e −ηt S l (Ψ)S l (Φ)). For the second remainder term, we bound We now use Lemma 3.2 to bound and since there is a uniform constant c > 0 such that P N K ⊆ P c K , noting that P K = N − A K , we can bound Combining the two remainder terms, and bounding all norms by S l (Ψ)S l (Φ) we get that where the implied constant is absolute.

General test functions.
We now use a partition of unity to reduce the case of a general test function to the case of functions with small support.
For ∈ (0, 1) sufficiently small, let Q be a maximal family of points in X ∩ Y such that the sets yB 3 , y ∈ Q , are disjoint and meet Y 2 , and let Note that the collection {yB 2 : y ∈ Q } covers X ∩ Y 2 and that {yB 3  Proof. Let y ∈ Q . If x ∈ yB 3 B 3 , then β y, (x) = 1 and hence z∈Q β z, (x) ≥ 1. Now suppose that x ∈ yB 2 yB 3 B 3 , in which case xB 3 ∩yB 3 = . Since y ∈ Y 3 and x ∈ yB 2 we have that x ∈ Y 2 ∩X . By the maximality of Q , there exists z ∈ Q such that xB 3 ∩zB 3 = . This implies that x ∈ zB 3 B 3 and hence β z, (x) ≥ 1.
s (x) . Since α y, is supported on yB 2 , we only need to bound its derivatives there in which case we have that  Proof. Since τ = y∈Q β y, y∈Q β y, , it is clear that 0 ≤ τ ≤ 1. Note that if y ∈ Q Q , then yB 2 ∩ Y 4 = . Hence if x ∈ X ∩ Y 4 satisfies x ∈ yB 2 , then β y, (x) = 0. This shows that y∈Q β y, (x) = y∈Q β y, (x). Moroever, since X ∩ Y 4 is covered by {yB 3 B 3 : y ∈ Q }, we have that y∈Q β y, (x) = 0 on X ∩ Y 4 and hence indeed τ = 1 there. Next, since for any y ∈ Q , we have that yB 2 ⊆ Y ; so τ (x) = 0 outside of Y . Finally we can bound S l (τ ) ≤ y∈Q S l α y, −p+3 dim(G) .

SHRINKING TARGET PROBLEMS
We now use the results on the exponential decay of matrix coefficients to prove an effective mean ergodic theorem and apply it to various shrinking target problems. As before, we assume that Γ is a geometrically finite, Zariski dense subgroup of G = SO(n, 1) • . For n ≥ 5, in the case where Γ has a cusp and δ ≤ n − 2, all functions and shrinking targets on X we consider below are assumed to be M -invariant so that Theorem 1.14 applies to them. All functions below are also assumed to be real-valued functions. where Ψ denotes the L 2 -norm of Ψ. In the entire section, we will take λ to be either the Lebesgue measure on R (when considering a continuous time flow) or the counting measure on Z (for a discrete time flow). For T ≥ 1, consider the averaging operator λ T on L 2 (X , m) given by

THEOREM 4.2.
For any non-zero Ψ ∈ C ∞ (X ), and for all T 1, where we used that λ is translation invariant and λ([0, T ) ∩ [t , t + T )) = T − |t | (where in the discrete case we may and will assume that t and T are integers). Now fix a large parameter M to be determined later. For |t | ≥ M large we use Theorem 1.14 to get that for some η 0 ∈ (0, 1). On the other hand, for |t | < M small, we bound 〈a t Ψ, Ψ〉 ≤ Ψ 2 , to get that using m(Ψ) ≤ Ψ . Using these estimates, we get that It remains to set M = 2 log(S * (ψ))η 0 −1 to finish the proof.
Following [13], for a non-negative function Ψ on X , we define Note that C o T,Ψ ⊆ C T,Ψ . As a direct consequence of the effective mean ergodic theorem, we get the following bounds: Proof. On one hand, On the other hand, by Theorem 4.2, Putting these two together gives the result.
Having control on the measures of these subsets has immediate consequences to several shrinking target problems. Indeed, a simple adaptation of [13,Lemmas 13 and 14] gives the following result.
for all t x 1.

Hitting along a subsequence.
In the rest of this section, let B = {B t } be a family of shrinking targets in X . Recall that a family B is inner regular (resp. outer regular) if there exist c > 0, α > 0 and smooth positive functions A family B is regular if it is inner and outer regular. When we want to emphasize the parameters c and α, we say that a family is (c, α)-regular. Our first application of the effective mean ergodic theorem is the following.

PROPOSITION 4.5. Assume that B is inner regular and satisfies
Then there is a subsequence t j → ∞ such that for m-a.e. x ∈ X , If B is also outer regular, then for m-a.e. x ∈ X , Id B t such that log(S * (Ψ t )) log(m(B t )) and m(Ψ t ) m(B t ). The mean ergodic theorem (Theorem 4.2) applied to Ψ t implies that where we used that Ψ t 2 ≤ m(B t ). From our assumption, there is some subsequence t j such that log(m(B t j )) m(B t j )·t j → 0. Hence λ t j (Ψ t j ) → 1 in L 2 (Γ\G, m) and, after JOURNAL OF MODERN DYNAMICS VOLUME 17, 2021, 401-434 perhaps passing to another subsequence, we get λ t j (Ψ t j )(x) → 1 for m-a.e x ∈ X .
For any x in this full measure subset, the inequality implies that λ({t ≤ t j : xa t ∈ B t j }) t j m(B t j ) as claimed. Assuming that {B t } is also outer regular, repeating the same argument for functions approximating Id B t from above gives the other inequality.
In particular, taking λ to be the Lebesgue measure gives the first part of Theorem 1.4. Moreover, by taking λ to be the counting measure, we get the following consequence implying a discrete version of Theorem 1.3(1).

COROLLARY 4.6. If B is inner regular and
Proof. Applying the above result with λ the counting measure shows that for along some subsequence t j . Since B t j ⊆ B k for any k ≤ t j , it follows that the subset {k : xa k ∈ B k } is unbounded as well.

Orbits eventually always hitting.
The results of the previous section allow us to control how orbits hit the shrinking targets along a subsequence of times, however, under the same hypothesis we could also have different subsequences for which this asymptotic fails, and for which the set {k ≤ k j : xa k ∈ B k j } may even be empty (see, e.g. [13,Proposition 12]). A more subtle question is to ask what conditions on the shrinking sets guarantee that the truncated orbit {xa j : j ≤ k} is eventually always hitting the targets B k , and moreover, how large is their intersection? This is the content of the following Theorem 4.7, which is a discrete version of Theorem 1.3 (2).

THEOREM 4.7. Assume that B is inner regular and that
for some sequence t j → ∞. Then for m-a.e. x ∈ X and for all t x 1, we have Proof. From the inner regularity, we can find smooth functions 0 ≤ Ψ t ≤ Id B t satisfying log(S l (Ψ t )) log(m(B t )) and m(B t ) m(Ψ t ). By Proposition 4.3, we can estimate thatfor any s, t > 1 , we obtain j m C t j −1 ,Ψ t j < ∞. Hence by the first part of Lemma 4.4, we have that for m-a.e. x ∈ X , λ t Ψ t (x) = 0 for all sufficiently large t . Taking λ to be the counting measure on N, this implies that {k ∈ N : k ≤ t , xa k ∈ B t } = for all sufficiently large t .

THEOREM 4.8. Assume that B is regular and that
then, for m-a.e. x, and for all t x 1, Proof. Let Ψ ± t be functions which approximate Id B t from above and below such that For each of these functions we can use Proposition 4.3 as before to So by the second and third part of Lemma 4.4 we get that for m-a.e. x ∈ X and for all sufficiently large t , we have This implies that λ t (Id B t ) m(B t ). Finally, taking λ to be the counting measure on N (resp. the Lebesgue measure) gives the result for discrete (resp. continuous) time flow.

Logarithm law for the first hitting time.
Using similar arguments utilizing the effective mean ergodic theorem, we can prove the logarithm law for the first hitting time for the discrete flow. Recall the discrete first hitting time function Proof. We first note that the bound holds for m-a.e.x; indeed, this holds in general for any monotone sequence of shrinking targets in a measure preserving dynamical system (see [14,Lemma 2.2]). It is thus sufficient to show that for m-a.e. x, Fix a small > 0 and set Note that if x ∈ A + , then there are arbitrarily large values of t for which Now for any j ∈ N, we choose y j ∈ 1 2 j +1 , 1 2 j such that either t j = sup{t : is continuous, we may simply take y j = 2 − j . In general, since the function t → m(B t ) is monotone decreasing, it has at most countably many points of discontinuity and hence we can always find such points). We partition [0, ∞) into intervals I j = {t : m(B t ) ∈ [y j +1 , y j )} and write For all sufficiently large j and any t with m(B t ) ∈ [y j +1 , y j ), we have that k (t ) ∈ We can thus further bound From our choice of y j and t j , we have that m(Ψ t j ) = y j ∈ 1 j . Using Proposition 4.3 for the smooth functions as before, we bound Hence m(A + ) ≤ j >k j 2 − j k2 − k for all k ∈ N. Therefore m(A + ) = 0 and This holds for any > 0. Hence, by taking a sequence of j → 0, we finish the proof.

4.5.
Thickening along the flow. We note that if {k ∈ N : xa k ∈ B k } is unbounded (resp. { j ≤ k : xa j ∈ B k } = ), then {t ∈ R : xa t ∈ B t } is unbounded (resp. {t ≤ k : Hence the same assumptions on the shrinking rate of m(B t ) as in Proposition 4.5 give the same conclusions also for the continuous flow. However, it is possible for the set {t ∈ (0, ∞) : xa t ∈ B t } to be unbounded even when it is bounded for the discrete time flow. In order to get the correct thresholds for the continuous flow, one needs to consider the thickened targets. For any set B ⊆ X we define its thickeningB to bẽ In the following lemma we observe that the shrinking target problems for the continuous flow can be translated to similar problems for the discrete flow hitting the thickened targets.

LEMMA 4.10.
For any B ⊆ X and x ∈ X , we have: The proof of these observations is easy once stated and we omit the details. Using this, we get the following sharper results for the continuous time flow, which imply Theorems 1.1 and 1.3. THEOREM 4.11. Suppose that the family {B t } t ≥1 of thickened targets is inner regular.
Proof. The first condition (with k replaced by k + 1) implies that the set {k ∈ N : xa k ∈B k+1 } is unbounded. For each k in this set, there is some t k ∈ [k − 1/2, k + 1/2] with xa t k ∈ B k+1 ⊆ B t k , proving the first part.
For the second part, the summability condition implies that for m-a.e. x, we have that {xa j : j ≤ k)} ∩B k = for all sufficiently large k > k 0 . Now for t ≥ k 0 + 1 and k := t , there is some j ≤ k with xa j ∈B k ; hence there is s ≤ t with xa s ∈ B k ⊆ B t .
Finally for the last part, since . REMARK 4.12. The problem of estimating {t ≤ k : xa t ∈ B k }, for the continuous time flow, does not easily reduce to the discrete time problem for the thickened targets. Here, knowing that xa k ∈B k only tells us that xa t ∈ B k for some t close to k but not on the amount of time spent there. Hence, to get asymptotics we need the stronger condition that ∞ j =1 | log(m(B 2 j ))| 2 j m(B 2 j ) < ∞ for the original sets and not the thickened sets. In particular, if m(B k ) k −a for some a ≥ 1 and m(B k ) k −b for some b < 1, then by reducing to the thickened case, we know that for all sufficiently large k, {t ≤ k : xa t ∈ B k } = , but we do not get an asymptotic estimate for the size of these sets.

EXPLICIT EXAMPLES
In this section, we consider explicit examples of shrinking targets given by shrinking cusp neighborhoods, shrinking metric balls and shrinking tubular neighborhoods, and show that they are regular and approximate their measure.

Cusp neighborhoods.
Let h 1 , . . . , h k and h i ,t be the cusp neighborhoods defined in (2.5). In order to apply our results for these sets we need to verify that the family {h i ,t } t ≥1 is regular and satisfies m(h i ,t ) e −t (2δ−κ i ) where κ i is the rank of the parabolic fixed points associated to h i . While the upper bound m(h i ,t ) e −t (2δ−κ i ) is proved in [2] and [23], we could not find a reference where the lower bound is established; so we include a proof for the convenience of readers.
The important feature of a geometrically finite group is that all of its parabolic fixed points are bounded, i.e., the stabilizer of ξ in Γ acts cocompactly on Λ−{ξ} for each parabolic fixed point ξ. This is the main ingredient of the argument below. We refer to [1] for the description of horoballs in geometrically finite manifolds that will be used below.
We will work here with the upper half space model and fix our base point to be o = (0, 1). Since we will work with one fixed cusp, we may assume without loss of generality that it is the infiniy ∞. Set Γ ∞ := Stab Γ ∞ and κ to be the rank of ∞. Without loss of generality, we assume that Γ ∞ = Z κ . Fix a horoballH (0) ⊂ H n such that In fact,H (0) is of the form {(x, y) : y = y 0 } for some y 0 > 0. For the notational simplicity, we assume y 0 = 1. Set Without loss of generality, we may assume π Choose a fundamental domain F ∞ ⊆ R n−1 for the action of Γ ∞ on R n−1 containing the origin so that the sets, int(γF ∞ ), are mutually disjoint for γ ∈ Γ ∞ . Note that H (t ) = {z = (x, y) : x ∈ F ∞ : y ≥ e t } is a fundamental domain for π(h t ) = Γ ∞ \H (t ). We can choose a compact fundamental parallelepiped P containing F ∞ ∩ Λ such that Γ ∞ P covers Λ {∞} and int(γP )s are mutually disjoint for all γ ∈ Γ ∞ . We may choose P to contain the origin so that if H (t ) := H (t ) ∩ hull(Λ), then As P is compact, we have for any z ∈ H (t ), we have d (Γo, z) = d (o, z), and for z ∈ ∂H (t ), . The following is also clear from (5.1): PROPOSITION 5.1. The injectivity radius r z at any point z ∈ ∂H (t ) satisfies r z e −t , where the implied constants are uniform for all t 1.
We will use the following well-known fact: There exists c > 0 such that for all t ≥ 0, Next we want to estimate the measure m(h t ) for large t . For any ξ − = ξ + ∈ ∂H n − {∞} and s ∈ R, we denote by ξ s the unit speed geodesic from ξ − to ξ + (where s is the signed distance from the highest point of the geodesic), and recall that this gives us the coordinates (ξ − , ξ + , s) parametrizing T 1 (M ). Let Λ = Λ {∞} and let P 0 = F ∞ ∩ Λ. We first show the following: Proof. Let F Γ ⊆ H n be a fundamental domain for Γ\H n containing o such that for t ≥ 0 sufficiently large, we have that F Γ ∩H (t ) = H (t ), so that m(h t ) = Now decomposing Λ as a union over translates γP 0 with γ ∈ Γ ∞ , we can rewrite where for the second line we made a change of variables ξ → γξ and in the last line we used thatH (t ) = γ∈Γ ∞ γH (t ).
In order to evaluate this, we need the following geometric estimate.
LEMMA 5.4. Let ξ − ∈ P 0 and ξ + ∈ γP 0 with γ ∈ Γ ∞ . Then there exists c > 0 such that Proof. Recall that o = (0, 1) and note that γo = (v, 1) for some v ∈ R n−1 . Since P 0 is a compact set containing the origin, then γP 0 is a compact set (of the same diameter) containing v and hence ξ − − ξ Now assume that ξ ∩ H (t ) = and let z 1 , z 4 ∈ H n be the first and second intersections of the geodesic ξ s with ∂H (0) and z 2 , z 3 the first and second intersections with ∂H (t ). Writing z i = (x i , y i ), we have that x 1  Recall the notationh t = ∪ |s|<1/2 G s h t .
With this observation together with Lemma 5.4, we get that Next, to estimate ν o (γP 0 ) = ν γo (P 0 ), we use the Γ-conformality to get that To estimate β ξ (γo, o), let z 1 , z 2 be the two points in the intersection of ∂H (0) and the geodesic connecting ξ to γξ. Then d (z 1 , o) and d (z 2 , γo) are uniformly bounded and β ξ (z 1 , (1). Plugging in this estimate gives as claimed.
For the thickened target, for any x ∈ h t and |s| ≤ 1/2, if xa t ∈ h t −1 , then h t ⊆ h t ⊆ h t −1 . Hence m(h t ) e −t (2δ−κ) as well.
Next we show regularity.
Since m(h t ) m(h t ±1 ), we also get that m(Ψ ± t ) m(h t ), implying that {h t } is regular.

Shrinking balls in Γ\G.
In this subsection, our goal is to show that for x ∈ supp(m), the family {xG : 0 < < 1} is regular as stated in Proposition 5.9. We may assume that o ∈ hull Λ and fix v o ∈ T o (H n ) so that M = Stab(v o ). For ξ ∈ ∂H n and > 0, let B ξ ( ) denote the Euclidian ball of radius around ξ. When Γ is convex co-compact, Sullivan's shadow lemma implies that ν o (B ξ ( )) δ , but when Γ has cusps, the measure ν o (B ξ ( )) fluctuates as → 0. Nevertheless, we have the following: LEMMA 5.7. For any ξ ∈ Λ, the following holds for all sufficiently small > 0: Proof. For ξ ∈ Λ let ξ t ⊂ hull(Λ) denote the unit speed geodesic ray connecting o to ξ. Let b(ξ t ) ⊂ ∂(H n ) denote the shadow at infinity of the hyperbolic hyperplane meeting ξ t orthogonally. Then If {ξ 1 , · · · , ξ k } denotes the set of all representatives of Γ-orbits in the set of parabolic limit points, and H ξ i ⊂ H n is a sufficiently deep horoball based at ξ i , then H := k i =1 Γ(H ξ i ) forms a family of disjoint horoballs. Now by [28, Theorem 2], we have where κ(ξ t ) is the rank of ξ i if ξ t ∈ Γ(H ξ i ) for some i , and κ(ξ t ) = δ otherwise. Now, to prove (1), let = e −t . First, if ξ t is not in H , the claim follows easily. Next, if ξ t ∈ Γ(H ξ i ), then κ(ξ t ) = κ i and ν o (b(ξ t )) e −δt +d (ξ t ,Γo)(δ−κ i )) . If κ i ≤ δ, then so that e −δt ν o (b(ξ t )) e −(2δ−κ max )t . This proves (1).
PROPOSITION 5.8. Let K ⊆ X be a compact subset. Let δ − = min{δ, 2δ − k max } and δ + = max{δ, 2δ − k min }. For any x ∈ K ∩ Ω, we have that for all 0 < < r x , where all the implied constants above are uniform over all x ∈ K .
We finish this section with the proofs of Theorems 1.11 and 1.12.

Shrinking tubular neighborhoods.
For a fixed closed geodesic C ⊂ T 1 (M ) and > 0, we set C = {x ∈ T 1 (M ) : d (C , x) < }. The proof of Theorem 1.9 follows as above from the following. Proof. Recall the notations X , X ( ) and Y ( ) = X − X ( ) from Section 2.3. They are all M -invariant subsets of Γ\G, and in the following proof, we will regard them as subsets in Γ\G/M . We can present C = [g 0 ]AM /M and an element of C is represented by [g 0 ]a t M for a unique 0 ≤ t < L where L is the length of C . Let 0 be sufficiently small so that C ⊆ Y ( 0 ) and let 0 < ≤ 0 . Let Q denote a maximal set of points x i ∈ C such that the sets x i G M are pairwise disjoint. Writing x i = x 0 a t i the condition that x i G ∩ x j G = imply that |t i − t j | ≥ and the maximality condition implies that |t i − t i +1 | ≤ 3 . Hence #Q L −1 . Since we can estimate Now Proposition 5.11 (2) implies that m(g i G M ) 2δ+1 where the implied constant does not depend on i . Summing over all x i ∈ Q , we get that indeed m(C ) 2δ .
Next, to show regularity, for each point x i ∈ Q , let Ψ ± ,i be smooth nonnegative functions approximating x i G M from below and x i G 3 M from above respectively, with S l (Ψ ± ,i ) −l , and define Ψ ± = i Ψ ± ,i . Since the sets x i G M are pairwise disjoint, we have that Ψ − ≤ Id C ≤ Ψ + and moreover and similarly that m(C ) m(Ψ − ). Since #Q −1 , we can bound S l (Ψ ± ) −(l +1) m(C ) −α with α = l +1 2δ , showing that the family {C } is (c, α) regular for some c > 1, and α = l +1 2δ . Finally, note that there is c ≥ 1 such that a −s G a s ⊆ G c for all |s| ≤ 1/2. Then any point x ∈C is of the form x = x 0 a t g a s M with 0 ≤ t ≤ L, g ∈ G and |s| ≤ 1/2. We can write g a s = a s a −s g a s ∈ a s G c , to get that x ∈ x 0 a t +s G c ∈ C c . Therefore C ⊆C ⊆ C c , implying that {C } is also regular with m(C ) m(C ).