A Microlocal Analysis of the Lévy Generator with Conjugate Points
Abstract.
We analyze the microlocal structure of the infinitesimal generator of a Lévy process on a closed Riemannian manifold when conjugate points are allowed. We show that if there are no singular conjugate pairs, then the infinitesimal generator can be written as a sum of pseudodifferential operators and Fourier integral operators. This extends and unifies known results for the flat torus, the sphere, and Anosov manifolds.
1. Introduction
In this paper, we will study the infinitesimal generator of a Lévy process on a Riemannian manifold . Loosely speaking, a Lévy process on a manifold looks like Brownian motion interlaced with jumps along geodesics at random times. The infinitesimal generator encodes certain information about the process and reflects the geometry of . For example, the generator of the most famous Lévy process, Brownian motion, is the Laplace-Beltrami operator. Though Brownian motion on manifolds is well studied [Elw88, Hsu02], the theory for more general Lévy processes has received less attention. Hunt introduced Lévy processes on Lie groups [Hun56], while Gangolli initiated the study of symmetric spaces [Gan64, Gan65]. Applebaum and Estrade were the first to construct Lévy processes on arbitrary Riemannian manifolds, under a natural assumption on the Lévy measure [AE00]. Our goal is to analyze the microlocal structure of the infinitesimal generator of this Lévy process.
One motivation for studying Lévy processes on manifolds is to shed light on the Lévy flight foraging hypothesis, which is the foundation of several biological models [Vis+96, BN13] and search algorithms [YD09, YD10, HP17, KKM22]. This controversial hypothesis claims that Lévy processes are a better model of animal foraging behavior than Brownian motion, in the sense that they can optimize search efficiencies [SK86, Vis+99]. Even though the underlying geometry is often curved, until recently this topic had only been studied in Euclidean space. The expected time for a pure jump Lévy process to find a small target on a manifold was first studied in [Cha+22], while [NTT21, NTT22, Nur+23] considered Brownian motion. A numerical comparison of [Cha+22] and [Nur+23] was performed in [TT23], confirming that Brownian motion might be the faster search strategy for small targets on the -torus. This suggests that the underlying geometry could be an important factor in determining whether the Lévy flight foraging hypothesis is valid in a given context.
In [Cha+22], conjugate points are one of the main geometric influences on the expected stopping time. Specifically, the authors show that on the sphere the expected stopping time exhibits singular behavior at antipodal points, but no such anomaly occurs on the flat torus or Anosov manifolds. (We will call an Anosov manifold if its geodesic flow is an Anosov flow on its unit sphere bundle [Ano69, Kni02], or, equivalently, if lies in the interior of the set of metrics on without conjugate points [Rug91].) Conjugate points exert influence on the expected stopping time through the microlocal structure of the Lévy generator (i.e., the infinitesimal generator of the pure jump Lévy process). On the sphere the Lévy generator decomposes into a pseudodifferential operator and a Fourier integral operator, while on the flat torus and Anosov manifolds it is simply a pseudodifferential operator. Inspired by these results, as well as the microlocal analysis of the geodesic X-ray transform done in [SU12, MSU15, HU18], we prove a similar theorem in a more general geometric context.
Our main result, Theorem 3.8, shows that if there are no singular conjugate pairs, then the Lévy generator equals a sum of pseudodifferential operators and Fourier integral operators. Each Fourier integral operator is associated with conjugate pairs of a given order, and the order of the operator depends on the order of those conjugate pairs and the dimension of . The canonical relation of each Fourier integral operator is related to the geometry of the set of conjugate pairs. When is the sphere or an Anosov manifold we recover the results of [Cha+22]. Our theorem also covers all possibilities in two dimensions.
The structure of the paper is as follows. In Section 2 we discuss the necessary background, including an initial decomposition of the Lévy generator into pseudodifferential operators and a remainder term. In Section 3 we prove our main result by showing that the remainder term equals a smoothing operator plus a sum of Fourier integral operators.
Acknowledgements
The author would like to thank Gunther Uhlmann for proposing this problem and for his consistent support, guidance, and patience. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE-2140004.
2. Preliminaries
This section introduces our notation, provides the main definitions, and recalls a result from [Cha+22] that will guide our analysis of the Lévy generator. Throughout the paper, we will assume satisfies the following condition.
Assumption 2.1.
is a smooth, closed (that is, compact without boundary), and connected Riemannian manifold of dimension .
2.1. The Lévy Generator
Before defining the Lévy generator, let us specify our notation for vector bundles and introduce the exponential map.
Unless stated otherwise, we will denote the projection map associated with a vector bundle by with a subscript indicating the total space of the bundle (e.g., is the projection map from the tangent bundle to ). When referring to a point in a vector bundle, one must decide whether to specify the base point. For example, should a point in be denoted by or simply ? We will often omit the base point, but we will include it when there is potential for confusion. We do follow the convention that means and , and similarly for points in other vector bundles.
For each , let be the maximal geodesic with initial data . Since we view as a point in , this initial data includes (loosely speaking) both the initial position and velocity of the geodesic. Each geodesic is defined on all of because is compact. Thus we can define the exponential map by
For each , the restriction of the exponential map to will be denoted by .
Fix . It is well known that the fractional Laplacian of order is the infinitesimal generator of a -stable jump process on Euclidean space. To ensure that the Lévy generator is consistent with the fractional Laplacian, we will include the constant
Definition 2.2.
The Lévy generator is defined for and by
where is the Riemannian density of induced by the metric .
2.2. Averaging Along the Geodesic Flow
This subsection introduces one of the operators which will appear in our initial decomposition of the Lévy generator.
Let be the unit sphere bundle over . Its fiber over a point is
Since is compact and -dimensional, is a compact manifold of dimension , and it is an embedded submanifold of . We will write rather than for the projection map , and we will use for the inclusion map .
The geodesic flow on is the smooth map defined by
It is a smooth submersion because is a diffeomorphism for each . We will denote the differential of the map at a vector by . Also, the geodesic flow on will be denoted by .
Since is a closed Riemannian manifold, its injectivity radius is positive and finite. Choose a bump function satisfying for and for . Let
Then we can define an operator by
(1) |
We may think of as a certain average of along the geodesic flow. This operator will play a crucial role in our analysis of the Lévy generator.
2.3. Pushforward and Pullback by a Smooth Submersion
Next we will define the pushforward and pullback and set our notation for Fourier integral operators.
Definition 2.3.
Let and be smooth manifolds of dimension and , respectively, with smooth positive densities and . Let be a smooth submersion. Then the pushforward and the pullback are defined by the requirement that
for all and .
Explicitly, the pushforward by integrates over the level sets of , while the pullback precomposes with .
We would like to use Definition 2.3 to express as a composition of simpler operators. Let be the operator which multiplies by the smooth function
Then is a properly supported pseudodifferential operator of order on . Let
be the projection map. Since integrates a function over , (1) suggests that equals the composition . This is not quite correct since may not have compact support, even if . But we can use a partition of unity to define the composition of and , and this is how we will view the operator .
Recall that if is a smooth map between smooth manifolds, then for each the differential yields a dual linear map
Since is given in coordinates by the transpose of the matrix of , if is a smooth submersion then is injective for all .
It is known, going back at least to [GS75], that the pushforward and pullback by a smooth submersion are both Fourier integral operators. We will state the version from Lemma of [HU18] to have this result in the precise format we need.
Lemma 2.4.
Suppose we are in the setting of Definition 2.3. Then the pushforward and the pullback are both Fourier integral operators of order . The canonical relation of is
while the canonical relation of is
We will use the notation of [Hör85] for Fourier integral operators. Given smooth manifolds and , the set consists of all operators which map to , and whose Schwartz kernel is a Fourier integral of order with Lagrangian given by the twisted canonical relation . Unlike [Hör85], we will allow to be merely a local canonical relation, meaning it can be immersed rather than embedded.
In this setup, Lemma 2.4 implies that
In the next subsection, we will see that the key to determining the microlocal structure of is to show that a certain composition involving these operators is a Fourier integral operator. Our main tool will be the clean intersection calculus for Fourier integral operators, whose statement may be found in [DG75], [Wei75], and [Hör85].
2.4. Initial Decomposition of the Lévy Generator
The next result, a consequence of Theorem in [Cha+22], is the starting point for our microlocal analysis of . The idea of the proof is to split into two parts: one where is injective, and a remainder term involving , , and . This decomposition is unnecessary if is the flat torus, since is the fractional Laplace-Beltrami operator in that case (Theorem in [Cha+22]). Thus we may safely rule out that case throughout the paper.
Theorem 2.5.
If preserves , then
where (resp. ) is a pseudodifferential operator of order (resp. ) on .
Proof.
Fix and . Write , where
By our choice of , the integrand of vanishes whenever , so we can make the change of variables . Then
where is the Riemannian distance from to , is the Jacobian determinant of the map , and is the Riemannian density of . Hence is a pseudodifferential operator of order on .
For , use polar coordinates to write
where is the Riemannian density of induced by the metric . Since preserves and integrates over the fibers of , we can write
Let be the operator which multiplies by the constant function . Since is a pseudodifferential operator of order on , this completes the proof. ∎
We will see in Theorem 3.1 that is a Fourier integral operator, so it preserves . Therefore the microlocal analysis of boils down to showing that is a Fourier integral operator, a task which we will take up in the next section.
3. Microlocal Structure of the Lévy Generator
This section contains our main result: a decomposition of into pseudodifferential operators and Fourier integral operators. In light of Theorem 2.5, a natural first step is to show that is a Fourier integral operator. Since we will eventually need to split into several pieces, we will actually introduce an arbitrary smooth function into and show this more general operator is a Fourier integral operator.
3.1. is a Fourier Integral Operator
Givan any , let
where is the operator which multiplies by . Since is a properly supported pseudodifferential operator of order on and is in , the operator is also in . Therefore, if we can show that the clean intersection calculus applies to the composition of and , then is a Fourier integral operator. This is the content of the following theorem.
Theorem 3.1.
Let
Then .
Proof.
As noted above, we can use a partition of unity to define the composition of and . Then by localizing and reducing to the case of two properly supported Fourier integral operators, the proof boils down to an analysis of the canonical relations.
Since is the canonical relation of , Lemma 2.4 implies that
Therefore, by the clean intersection calculus, it is enough to prove the following:
-
(i)
The intersection
(2) is clean in the sense that is an embedded submanifold, and at every point the tangent space equals the intersection of the tangent spaces of the two manifolds being intersected.
-
(ii)
The projection map is proper.
-
(iii)
For every , the fiber is connected.
If we are willing to work with local canonical relations, then point (iii) can be omitted. Point (iii) will hold in this case, but in later results it may not.
Consider the map
defined by
If we include the base points, then
(3) |
so is a smooth embedding. The injectivity of implies that is injective. Hence point (iii) holds, and assuming points (i) and (ii) are true the excess of the intersection (2) is zero.
To begin proving point (i), let be the smooth rank- subbundle of whose fiber over each point is
In other words, is the subspace of which is conormal to the kernel of . To see that is a smooth subbundle, just use the rank theorem to write
locally, and note that is a smooth local frame for . Because is a linear isomorphism from to , we can define a smooth map
(4) |
whose restriction to each fiber is the inverse .
The domain of our smooth parametrization of will be the set
We claim that is an embedded submanifold of dimension . Let
(5) |
be the restriction to of the projection map
Since equals the level set , the claim is true if is nonzero for all . Let . Then for some . Choose slice coordinates for near and near , and fix natural coordinates for associated with the latter. Then locally we can write and . Since is a unit-speed geodesic, we may suppose without loss of generality that
Define a curve in as the composition of and the curve
Then is a smooth curve in such that . Moreover,
(6) |
Therefore, as claimed, is an embedded submanifold of dimension .
Now we can use to parametrize . Indeed, if is the inverse of (3) composed with the projection onto the component, then the map
defined by
is a smooth embedding, so is an embedded submanifold of dimension .
To finish proving point (i), fix . Since is necessarily contained in the intersection of the tangent spaces of the manifolds on the right-hand side of (2), it is enough to show the reverse containment. Observe that can be parametrized by the map
defined by
Hence any vector is the velocity of some smooth curve
meaning and
The vector is also in if and only if
Thus, in any local coordinates, and agree to first order at . Then (3) implies that the same is true of and . Hence
which means . This completes the proof of point (i).
For point (ii), suppose we have a compact set . Then there exists a constant such that
Hence is precisely the set
which is compact by continuity. Then is compact, because it is a closed subset of the compact set . Since is a diffeomorphism onto , this implies that is a proper map. Thus point (ii) holds.
Since points (i)-(iii) hold and the excess is zero, the map is a smooth embedding, so its image is an embedded submanifold. Since and are Fourier integral operators of order , we conclude that . ∎
Because is a Fourier integral operator with canonical relation , its twisted wave front set must be contained in . In the next lemma, we will leverage the fact that is a composition of Fourier integral operators to say a bit more.
Lemma 3.2.
The twisted wave front set of is contained in
Proof.
Since equals , we know
Because the Schwartz kernel of is smooth away from , and is supported away from , we also know
Putting these two containments together yields the result. ∎
This simple observation will be useful when we split into several pieces.
3.2. Composition with
Recall that our goal is to understand the microlocal structure of . Theorem 3.1 demonstrated that is a Fourier integral operator for any . The next theorem shows the same is true of the operator
Theorem 3.3.
Let
Then .
Proof.
Using Theorem 3.1 and Lemma 2.4, one can check that equals , so it suffices to show points (i)-(ii) in the proof of Theorem 3.1 hold with (2) replaced by
(7) |
(Since we omit point (iii), in general will only be a local canonical relation.)
To begin proving point (i), let be the smooth rank- subbundle of whose fiber over each vector is
Similar to (4), we can define a smooth map
whose restriction to each fiber is the inverse . Let
Since is a smooth submersion and is an embedded codimension- submanifold of , we know is an embedded codimension- submanifold of .
Using a similar argument as the proof of Theorem 3.1, we will show that
is an embedded submanifold of dimension . Let be the restriction to of the function defined in (5). Then equals the level set , so it suffices to show that is nonzero for all . Let . Then for some , we have
Choose slice coordinates for near , which yield corresponding coordinates for near . Fix natural coordinates on . Then locally we can write . Define a curve in as the composition of and the curve
Then is a smooth curve in such that . Moreover, similar to (6),
This proves that is an embedded submanifold of dimension .
Now we can use to parametrize via the map
defined by
By the last paragraph of the proof of Theorem 3.1, the map
is a smooth embedding, and hence so is its restriction to . Thus is a smooth embedding, so is an embedded submanifold of dimension .
To complete the proof of (i), note that is parametrized by the map
given by
Fix and suppose . Then there is a smooth curve
such that and the velocity of this curve at zero is . As in the proof of Theorem 3.1, if as well, then the derivatives of and agree at in any local coordinates. It follows that
which means . Therefore the intersection (7) is clean and the excess is zero.
The proof of point (ii) is the same as Theorem 3.1, so we omit the details. Thus is a local canonical relation. Since and are Fourier integral operators of order and , respectively, we conclude that . ∎
When , we will write and instead of and . Then by Theorem 2.5,
Hence the microlocal analysis of reduces to understanding the composition .
The main difficulty in this case is that may have multiple connected components. One component corresponds to a smoothing operator, and the others appear only when there are conjugate points. If we rule out certain types of conjugate points, then these additional components give rise to Fourier integral operators whose canonical relations and orders can be determined. This is the content of Theorem 3.8, our main result. In the next subsection, we will provide the additional definitions and lemmas needed to state and prove it.
3.3. Conjugate Pairs
Though conjugate points along a geodesic are often defined in terms of vanishing Jacobi fields, it will be more convenient to work with the corresponding velocity vectors of the geodesic instead. This leads us to the notion of a conjugate pair.
Definition 3.4.
We call a conjugate pair if
If the dimension of is , then is a conjugate pair of order . The set of regular conjugate pairs of order , denoted by , is the set of conjugate pairs which have a neighborhood in such that all other conjugate pairs in have order . The set of singular conjugate pairs, denoted by , is the set of conjugate pairs which are not in for any .
By Lemma in [HU18], Definition 3.4 is equivalent to the traditional definition of conjugate points in terms of vanishing Jacobi fields along a geodesic.
The crucial assumption in Theorem 3.8 is that there are no singular conjugate pairs. This matters because, as the next lemma shows, the set of regular conjugate pairs of order is a smooth manifold, which may not be true of the set of all conjugate pairs.
Proposition 3.5.
For each integer , the set is an embedded -dimensional submanifold of , and the set
is a smooth vector bundle of rank over .
Proof.
For the first point, it is enough to show that each point in has a neighborhood in such that is an embedded submanifold of dimension . To prove this local statement, we will extend the methods of [War65].
Fix . Let be the differential in the fiber variables of the exponential map . By [War65], we can find coordinate neighborhoods of in and of in such that the st elementary symmetric polynomial in the eigenvalues of (denoted by ) has nonzero derivative in the radial direction. Then is an embedded -dimensional submanifold of , and it equals the set of vectors in with conjugate points of order or higher in .
Choose a neighborhood of in such that all other conjugate pairs in have order . Supposing without loss of generality that , we may assume . Consider the smooth map defined by
Then is a smooth immersion and satisfies
Since has a continuous inverse defined on its image by
it follows that is an embedded submanifold of dimension .
To prove the second point, let be the pullback of by the map
Then is the kernel of the smooth bundle homomorphism
This map has constant rank , and we can view it as a bundle homomorphism over by pulling back by . Hence is a smooth rank- subbundle of . ∎
Next we will turn into a symplectic manifold and make some remarks. Let be the canonical symplectic form on , and let be the musical isomorphism induced by the metric . Then we can define a symplectic form on by
Then is a symplectomorphism for each , the kernel of is a Lagrangian subspace of for each , and in natural coordinates on we have
(8) |
In turn, induces a smooth bundle isomorphism defined by
where is interior multiplication by . We will denote the inverse of by .
The next lemma defines a smooth bundle homomorphism that will help us describe the canonical relations of the various pieces of .
Lemma 3.6.
For each integer , there is a smooth bundle homomorphism
defined by the requirement that
(9) |
The proof of this result is essentially the same as that of Lemma in [HU18], so we do not repeat the details here.
Our final lemma is the key geometric tool in the proof of Theorem 3.8. It will allow us to split into different pieces corresponding to different orders of conjugate pairs, each of which is associated with a Fourier integral operator.
Lemma 3.7.
Proof.
To see that (10) and (11) are equivalent, first note that
by (3) and the fact that . Hence (10) holds if and only if
Thus (11) implies (10). For the converse, just apply both sides of (10) to the vector and deduce that .
Before addressing the claims about conjugate pairs, let us make a few observations. It will be useful to work with rather than . To connect the two, observe that
(12) | ||||
(13) |
due to the identities and . Let
(14) |
Equivalently, applying , we have
(15) |
In the next paragraph, we will prove the following analogue of the first condition in (11):
(16) |
Assuming (16) for the moment, the fact that is a symplectomorphism implies
(17) |
Now suppose (11) holds and . We will divide the proof that is a conjugate pair into three steps. The first one is to prove (16). By (12), (13), and (11),
Since is the span of the differential of , this means
for some . Applying both sides to a radial vector , we find
Since by assumption and the vector is parallel to , this implies that and hence completes the proof of (16).
The second step is to show that (14) is in . Since (15) implies that vanishes on , which is a Lagrangian subspace, must be in . Choose normal coordinates on centered at such that
(18) |
and let be natural coordinates on . Then is the span of the vectors . Since is in , we can write
for some . Using (15), (18), and the assumption , we find
Using (8) and the fact that at , this implies . Hence is in , so we can define the vector .
The third step is to show that is in , meaning
Since is in and , we know is in . Hence it is enough to prove that
By the transposes of (12) and (13), this is equivalent to showing
But (17) implies that vanishes on the Lagrangian subspace , so is indeed in . This proves that is a conjugate pair.
Now suppose . Since is in , this means
Using (15) in the first line below and (13) (transposed) and (17) in the second, we find
Hence , which proves that .
Conversely, suppose . Unpacking (9), this means
(19) | ||||
(20) |
Let . Then (20) and the transpose of (13) imply that
Using that is a symplectomorphism together with (19), we find
Hence (16) holds in this direction of the proof as well. Applying to both sides of that equation and rewriting with (12) and (13) establishes the first condition in (11).
3.4. Final Decomposition of the Lévy Generator
Now we can state our main result. It refines Theorem 2.5 by decomposing into a smoothing operator and a sum of Fourier integral operators, assuming there are no singular conjugate pairs.
Theorem 3.8.
Suppose . Then for to , the sets
are either local canonical relations or empty. Let be the connected components of . Let and be the pseudodifferential operators from Theorem 2.5. Then
where is a smoothing operator, and for each either
or and if .
Proof.
We must decompose into a smoothing operator and a sum of Fourier integral operators. Though the clean intersection calculus does not directly apply to , we will cut up this operator so that it applies to each separate piece.
First let us describe . By Lemma 2.4 and Theorem 3.3,
The requirement in is precisely (10), which is equivalent to (11) by Lemma 3.7. We can also use this lemma to cut up into several pieces according to different orders of conjugate pairs. Indeed, if we take and any such that , then , , and satisfy (11), so one piece is the diagonal
If and , then (11) is satisfied for some , and because is injective. Then Lemma 3.7 and the assumption imply that
Our goal is to write as a sum of Fourier integral operators, each having a canonical relation contained in a single set of this union. Since , we can find open subsets of with disjoint closures such that for to . Write
where each open set contains exactly one of the connected components of . Then we can construct a partition of unity on such that
and for to and to . Setting , we have
By Lemma 3.2, the twisted wave front set of is contained in
Therefore we obtain smoothing operators except near points in
Since Lemma 3.7 implies that is empty,
is a smoothing operator. Lemma 3.7 also implies that for to and to . If any of these compositions are empty, we can absorb the corresponding operator into and set .
It remains to show that if is nonempty, then the clean intersection calculus applies to . We again refer to points (i)-(ii) in the proof of Theorem 3.1. (Since we omit point (iii), we only obtain local canonical relations in general.) Define
Let be the restriction of to the th connected component of , and let be the th component function of for . Then Lemma 3.7 implies that is a connected component of
(21) | |||
Thus, to show is an embedded submanifold, it is enough to prove that
(22) |
is a smooth embedding. Let be the projection of the th connected component of onto . Let be the pullback of by . Then is a smooth subbundle of (as in the proof of Proposition 3.5), so it suffices to show that the extension of (22) to is a smooth embedding. Observe that the bundle homomorphism
covers and is a smooth embedding (because it is a proper injective immersion). This implies that is a smooth embedding, so is smoothly isomorphic to the restriction of to . Hence it suffices to show that
is a smooth embedding. But (8) implies that this map is injective in each fiber, and hence a smooth bundle isomorphism onto its image. Therefore (22) is a smooth embedding, which proves that is a connected embedded submanifold of dimension .
Now let be arbitrary, and consider the set
Since has dimension and is contained in , the intersection is clean if the dimension of is at most . Suppose . Then , and an examination of shows that determines . Hence the dimension of is at most (the dimension of ). But for fixed , the set
is a -dimensional vector space contained in
Therefore the dimension of is at most , so the intersection is clean with excess . Because is a component of (21), the projection map
is proper by the same argument as the proof of Theorem 3.1. Since and are Fourier integral operators of order and , respectively, we conclude that
is in whenever the set is nonempty. ∎
Two special cases of Theorem 3.8 are worth mentioning. First, if is Anosov then it has no conjugate points [Rug91], so each Fourier integral operator is zero and we recover Theorem in [Cha+22]. Second, Theorem 3.8 covers all possibilities in two dimensions because singular conjugate pairs cannot exist (since conjugate pairs can only have order ). In higher dimensions, the generic case includes singular conjugate pairs [Arn72, Klo83].
References
- [Ano69] D.V.“@ Anosov “Geodesic flows on closed Riemann manifolds with negative curvature” American Mathematical Society, 1969
- [AE00] D.“@ Applebaum and A.“@ Estrade “Isotropic Lévy processes on Riemannian manifolds” In Ann. Probab. JSTOR, 2000, pp. 166–184
- [Arn72] V.I.“@ Arnol’d “Normal forms for functions near degenerate critical points, the Weyl groups of , , and Lagrangian singularities” In Funct. Anal. App. 6 Springer, 1972, pp. 254–272
- [BN13] P.C.“@ Bressloff and J.M.“@ Newby “Stochastic models of intracellular transport” In Rev. Mod. Phys. 85.1 APS, 2013, pp. 135
- [Cha+22] Y.“@ Chaubet, Y.G.“@ Bonthonneau, T.“@ Lefeuvre and L.“@ Tzou “Geodesic Lévy Flights and Expected Stopping Time for Random Searches” In arXiv preprint arXiv:2211.13973, 2022
- [DG75] J.J.“@ Duistermaat and V.“@ Guillemin “The spectrum of positive operators and periodic geodesies” In Invent. Math 29, 1975, pp. 39–79
- [Elw88] K.D.“@ Elworthy “Geometric aspects of diffusions on manifolds”, 1988 Springer
- [Gan64] R.“@ Gangolli “Isotropic infinitely divisible measures on symmetric spaces” In Acta Math. 111, 1964, pp. 213–46
- [Gan65] R.“@ Gangolli “Sample functions of certain differential processes on symmetric spaces” In Pac. J. Math. 15.2 Mathematical Sciences Publishers, 1965, pp. 477–496
- [GS75] V.“@ Guillemin and D.“@ Schaeffer “Fourier integral operators from the Radon transform point of view” In Proc. Symp. Pure Math. 27, 1975, pp. 297–300
- [HP17] A.A.“@ Heidari and P.“@ Pahlavani “An efficient modified grey wolf optimizer with Lévy flight for optimization tasks” In Appl. Soft Comput. 60 Elsevier, 2017, pp. 115–134
- [HU18] S.“@ Holman and G.“@ Uhlmann “On the microlocal analysis of the geodesic X-ray transform with conjugate points” In J. Differ. Geom. 108.3 Lehigh University, 2018, pp. 459–494
- [Hör85] L.“@ Hörmander “The analysis of linear partial differential operators IV: Fourier integral operators” Springer, 1985
- [Hsu02] E.P.“@ Hsu “Stochastic analysis on manifolds” American Mathematical Society, 2002
- [Hun56] G.A.“@ Hunt “Semi-groups of measures on Lie groups” In Trans. Am. Math. Soc. 81.2, 1956, pp. 264–293
- [KKM22] W.“@ Kaidi, M.“@ Khishe and M.“@ Mohammadi “Dynamic levy flight chimp optimization” In Knowl.-Based Syst. 235 Elsevier, 2022, pp. 107625
- [Klo83] F.“@ Klok “Generic singularities of the exponential map on Riemannian manifolds” In Geom. Dedicata 14.4 Springer, 1983, pp. 317–342
- [Kni02] G.“@ Knieper “Hyperbolic dynamics and Riemannian geometry” In Handbook of dynamical systems 1 Elsevier, 2002, pp. 239–319
- [MSU15] F.“@ Monard, P.“@ Stefanov and G.“@ Uhlmann “The geodesic ray transform on Riemannian surfaces with conjugate points” In Commun. Math. Phys. 337.3 Springer, 2015, pp. 1491–1513
- [Nur+23] M.“@ Nursultanov, W.“@ Trad, J.C.“@ Tzou and L.“@ Tzou “The narrow capture problem on general Riemannian surfaces” In Differ. Integral Equ. 36.11/12 Khayyam Publishing, Inc. West Palm Beach, FL, USA, 2023, pp. 877–906
- [NTT22] M.“@ Nursultanov, W.“@ Trad and L.“@ Tzou “Narrow escape problem in the presence of the force field” In Math. Methods Appl. Sci. 45.16 Wiley Online Library, 2022, pp. 10027–10051
- [NTT21] M.“@ Nursultanov, J.C.“@ Tzou and L.“@ Tzou “On the mean first arrival time of Brownian particles on Riemannian manifolds” In J. Math. Pures Appl. 150 Elsevier, 2021, pp. 202–240
- [Rug91] R.O.“@ Ruggiero “On the creation of conjugate points” In Math. Z. 208 Springer, 1991, pp. 41–55
- [SK86] M.F.“@ Shlesinger and J.“@ Klafter “Lévy walks versus Lévy flights” In On growth and form: Fractal and non-fractal patterns in physics Springer, 1986, pp. 279–283
- [SU12] P.“@ Stefanov and G.“@ Uhlmann “The geodesic X-ray transform with fold caustics” In Anal. PDE 5.2 Mathematical Sciences Publishers, 2012, pp. 219–260
- [TT23] J.C.“@ Tzou and L.“@ Tzou “Challenging the Lévy Flight Foraging Hypothesis-A Joint Monte Carlo and Numerical PDE Approach” In arXiv preprint arXiv:2302.13976, 2023
- [Vis+96] G.M.“@ Viswanathan, V.“@ Afanasyev, S.V.“@ Buldyrev, E.J.“@ Murphy, P.A.“@ Prince and H.E.“@ Stanley “Lévy flight search patterns of wandering albatrosses” In Nature 381.6581 Nature Publishing Group UK London, 1996, pp. 413–415
- [Vis+99] G.M.“@ Viswanathan, S.V.“@ Buldyrev, S.“@ Havlin, M.G.E.“@ Da Luz, E.P.“@ Raposo and H.E.“@ Stanley “Optimizing the success of random searches” In Nature 401.6756 Nature Publishing Group UK London, 1999, pp. 911–914
- [War65] F.W.“@ Warner “The conjugate locus of a Riemannian manifold” In Am. J. Math 87.3 JSTOR, 1965, pp. 575–604
- [Wei75] A.“@ Weinstein “On Maslov’s quantization condition, Fourier Integral Operators and Partial Differential Equations, J. Chazarain, ed.” Springer-Verlag, BerlinNew York, 1975
- [YD09] X.S.“@ Yang and S.“@ Deb “Cuckoo search via Lévy flights” In 2009 World congress on nature & biologically inspired computing (NaBIC), 2009, pp. 210–214 Ieee
- [YD10] X.S.“@ Yang and S.“@ Deb “Eagle strategy using Lévy walk and firefly algorithms for stochastic optimization” In Nature inspired cooperative strategies for optimization (NICSO 2010) Springer, 2010, pp. 101–111