![Arxiv:2008.07195V3 [Math.PR] 26 Nov 2020 (1) As: Functions Smooth on Xoeta Ucinl;Bueo’ Identity](http://data.docslib.org/img/f43586387783278a544201f025130853-1.webp)
arXiv:2008.07195v3 [math.PR] 26 Nov 2020 (1) as: functions smooth on xoeta ucinl;Bueo’ identity. Bougerol’s functionals; Exponential generalisa a of context the if in [3] Besides, and [2] in studied was which arcs(e lo[][].U oorbs nweg,tesp the knowledge, best our to Up [8],[9]). so-calle the also analogue by (see matrix-variate given matrices the distribution that stationary fact with pseudo-J the ergodic (or by Routh-Romanovski justified to be may relation misle in be may ago and [6] long-time from borrowed is ‘Hua-Pickrell’ name The nue h ttoaiyo h Pdffso u loisinvar its also but diffusion HP the r u (s of The [34] stationarity equation. work the Pearson Wong’s ensures are the E. solve generators in distributions whose time stationary diffusions first stationary the classified for author appeared diffusion HP (2) edo rmtecneto thswt xoeta function exponential with change variable has the it precisely, connection More the from off read − → h u-ikel(P rcs sadffso audi h whol the in valued diffusion a is process (HP) Hua-Pickrell The e od n phrases. and words Key u XLCTEPESOSO H U-IKELSEMI-GROUP HUA-PICKRELL THE OF EXPRESSIONS EXPLICIT ftesm-ru est ftemlidmninlHua-Pic multi-dimensional the Karlin-McG the of Assiotis. use density we semi-group section, the last the of e In and Laplacian field. Maass magnetic the Pro the of (S´eminaire density de semi-group Shiraishi the and to appeal Matsumoto Alili, Hua-Pickrell by the of raised Af representations integral cases. novel general The and viding Girsanov’s stationary and the stati between hand necessarily one connection the the (not on integral general Beta the Cauchy cor Using on diffusions t Hua-Pickrell focus through between we relation density intertwining Next, semi-group its which cut. of for the expression Wong E. unified by a considered down case stationary the visiting Abstract. hnsatn tzr.Mr eeal,teln-iebehavi long-time the generally, More zero. at starting when Y t ( A,K nti ae,w td h n-iesoa u-ikeldi Hua-Pickrell one-dimensional the study we paper, this In ) t , ≥ L u-ikeldffso;RuhRmnvk oyoil;As polynomials; Routh-Romanovski diffusion; Hua-Pickrell 0 A,K stedffso orsodn o()satn at starting (2) to corresponding diffusion the is d =( + (1 := d 2 2 y OA RSAADNZRDEMNI NIZAR AND ARISTA JONAS + y (2 u sinh( = 2 A 1. ) du − d Introduction 2 )tanh( 1) 2 u (2 + rnfrs()it h olwn generator: following the into (1) transforms ) 1 Au y + ) + epnigt ieetst fparameters. of sets different to responding K eigopdniy nwrn question a answering density, semi-group cosh( Pmaueo h pc fHermitian of space the on measure HP d ) ewrs epoeormi eutpro- result main our prove we terwards, rl atcesse nrdcdb T. by introduced system particle krell du srcint hs aaeesntonly not parameters these to estriction aii´s 5 01.T hsed we end, this To 2001). babilit´es, 35, K tn tt ueyiaiayvle of values purely-imaginary to it xtend d eo oml odrv nexpression an derive to formula regor ino h lsia ogrlsidentity. Bougerol’s classical the of tion espl mte eal n write and details omitted supply we nr)cs o hc epoean prove we which for case onary) ca instance ecial tr-iuil prtr n whose and operators Sturm-Liouville eascae eedefnto in function Legendre associated he rmo h te ad ediscuss we hand, other the on orem ,K A, , dn ic h prtr()appeared (1) operator the since ading y l fBona oin ihdrifts. with motions Brownian of als cb)plnmas oehls,it Nonetheless, polynomials. acobi) fti iuinitoue n[]is [6] in introduced diffusion this of ac ndsrbto ne h flip the under distribution in iance ) dy d elln hs eeao acts generator whose line real e ro h Pdffso a be may diffusion HP the of or , eEape().Tee the There, (E)). Example ee uin esatb re- by start We ffusion. ∈ R . oitdLgnr function; Legendre sociated < A Y 0 ( A,K 1 / 2 ) K , ∈ R 0 = hnthe then , , fthe of HP process sinh(Y (A,K)), t 0, has the same distribution as ([2], Proposition 1): t ≥ − 2t − B(A 1/2) (A,K) −B(A 1/2) (K/2) e 2t sinh(Y )+ e s dγ , t 0, 0 s ≥ Z0 where B(A−1/2) and γ(K/2) are independent standard Brownian motions with drifts A 1/2 and − K/2, respectively. In particular, this connection answers a question raised in [34] (see Section 3 there). Moreover, the time-reversal property of the Brownian motion shows that for fixed t > 0, (A,K) sinh(Yt ) is distributed as 2t (A,K) (A−1/2) (A−1/2) B2t Bs (K/2) (3) sinh(Y0 )e + e dγs . Z0 Consequently, the Dubins-Scharwtz Theorem implies that the HP process converges in distribution as t + if A< 1/2 and oscillates otherwise. → ∞ At the geometrical side, the operator d2 (4) + (N + 1) tanh(y), d2y where N 0 is a positive integer, is the radial part of the Laplace-Beltrami operator of the ≥ one-sheeted hyperboloid (de Sitter space): x2 + + x2 x2 = 1 SO (N + 2, 1)/SO (N + 1, 1). { 1 · · · N+2 − N+3 }≈ 0 0 Here, SO0(N + 1, 1) is the identity component of the pseudo-orthogonal group SO(N, 1) (see [16], part II, and references therein). In particular, if N = 1 then the corresponding de Sitter space is realised as: SO(3, 1)/SO(2, 1) Sl(2, C)/Sl(2, R), ≈ and the radial part of its Laplace-Beltrami operator is the generator of the process (B + ǫt,t 0), t ≥ where ǫ is a symmetric Bernoulli random variable which is independent of B ([2], see also Example 2 in [29]). Another striking and interesting observation made in [2] reveals that the operator (4) coincides also with the Laplace-Beltrami operator in the real hyperbolic space HN+2 RN+3 ⊂ written in equidistant coordinates. In the mathematical physics realm, the Sturm-Liouville operator LA,K may be mapped into the Hamiltonian with hyperbolic scarf potential: d2 sinh(u) 1 (5) + C1(A, K) + C2(A, K) , u R, − du2 cosh2(u) cosh2(u) ∈ where C1(A, K),C2(A, K), are constants depending on (A, K) (see e.g. [26], [30]). In representation-theoretical words, this Hamiltonian admits a SU(1, 1)-symmetry ([11], [18]). Despite the multiple occurrences of the HP generator LA,K in its various forms and though it is a Sturm-Liouville operator, still many rich clarifications related to its semi-group density are signi- ficant being highlighted. To start with, there is a typo in Wong’s expression of the eigenfunctions corresponding to the continuous spectrum and the present work develops the details leading to the correct one. Next, we prove that the expansion of the semi-group density in the stationary case given in [34] may be written in a unified way using the associated Legendre function in the cut. The newly-obtained expression has the merit to extend to the HP particle system after the use of the Karlin-McGregor and Andreief formulas. The main ingredient used to prove the unified expression is a connection we prove between ‘symmetric’ Romanovski polynomials and associated Legendre function in the cut. This connection provides also the details of the derivation of another formula 2 given in [34] for the semi-group density in the case when A is half of integer. For general parameters (A, K) (not necessarily stationary case), we exhibit an intertwining property between stationary and non stationary HP processes with parameters (A, K) and (2 A, K) respectively. This prop- − − erty extends the one given in [22] when K = 0. Using Cauchy Beta integral on the one hand and some absolute continuity relations given by Girsanov’s Theorem on the other hand, we discuss the connection between the stationary and general cases. Afterwards, we prove our main result and derive integral representations of the HP semi-group density valid for any A = 1, K R, avoiding 6 ∈ the Sturm-Liouville theory and rather inverting the Fourier transform derived in [3], Proposition 6.1. To this end, we use the analytic extension of the semi-group density of the Maass Laplacian to purely-imaginary values of the magnetic field and appeal to the connection between the Maass and the Morse semi-groups (see [17]). In particular, when A = 1, K = 0, we get a new expression of 6 the HP semi-group density as a partial Mellin transform of the semi-group density of the Brownian motion in the hyperbolic plane. Finally, we give some interest in the multivariate analogue of the HP process introduced and studied by T. Assiotis ([6]). This is a one-dimensional particle system depending on a complex parameter and realised as a Doob transform of independent stationary HP processes. Using the Karlin-McGregor and Andreief formulas, we express the semi-group density of this particle system with real parameter through multivariate symmetric Legendre functions.
1.1. Organisation of the paper. In Section 2 we revisit the stationary case of Wong’s paper [34]. There, we correct the eigenfunction expression corresponding to the continuous spectrum and prove the connection between ‘symmetric’ Romanovski polynomials and associated Legendre function in the cut. We close that section by describing the main steps toward Wong’s integral representation of the HP semi-group density proved for half-integer values of A. In Section 3, we consider the HP diffusion with general parameters for which we first recall the spectral decomposition of its generator then prove the intertwining property described above. We also prove in the same section an identity transforming a ‘symmetric’ Romanovski polynomial to a ‘non-symmetric’ one which may be seen as another type of intertwining since it involves the Cauchy Beta kernel. These intertwinings are next related to Girsanov’s Theorem connecting the stationary and general cases. In Section 4, we prove our main result and derive an integral representation of the HP semi-group density. In the last section, we recall basic facts on the HP particle system and derive an expression for its semi-group density when the parameter it depends on is real.
1.2. Acknowledgements. The authors appreciated stimulating discussions with R. Chhaibi, F. Baudoin, L. Miclo, H. Matsumoto and N. O’Connell. The first author gratefully acknowledges the support of the European Research Council (Grant number 669306).
2. The stationary case: E. Wong’s paper revisited In [34], example (E), the following Fokker-Planck equation was considered:
∂2[(1 + u2)p] + (2α 1)∂ (up)= ∂ p, u − u t where α > 0 and u R. Using the separation of variables technique, the following singular ∈ Sturm-Liouville equation was obtained (see eq. (12) in [34]): d 1 dφ λ (u) = φ(u), du (1 + u2)α−1/2 du −(1 + u2)α+1/2 3 or equivalently, d2φ dφ (6) (1 + u2) + (1 2α)u = λφ. d2u − du −
The diffusion operator displayed in the LHS of (6) is an instance of (1): L(1/2)−α,0. Its discrete spectrum consists of the finite sequence of eigenvalues ([30]): λ = n(2α n), 0 n M := [α 1], n − ≤ ≤ − and the corresponding eigenfunctions are ‘symmetric’ Routh-Romanovski polynomials: dn (7) R(α)(u) := (1 + u2)α+(1/2) [(1 + u2)n−α−(1/2)]. n dun They are finitely-orthogonal with respect to the symmetric (Student) weight ([30]): 1 (8) W (u) := , (1 + u2)α+1/2 (−α−1/2,−α−1/2) and may be expressed through Jacobi polynomials Jn (iu) with imaginary variable and negative parameters as (whence the name pseudo-Jacobi polynomials): (α) n (−α−1/2,−α−1/2) Rn (u)= cni Jn (iu) ((1/2) α) 1 1 iu = c in − n F n,n 2α, α; − , n M. n n! 2 1 − − 2 − 2 ≤ n Here cn = ( 2) n! is a normalisation constant which may be computed by comparing the leading − (α) coefficients of both sides (see for e.g. [30], eq. (31), for the leading coefficient of Rn ) and 2F1 L is the Gauss hypergeometric function ([1]). As to the continuous spectrum of (1/2)−α,0, it is parametrized by µ R: ∈ λ = α2 + µ2 and the corresponding eigenfunctions were given in [34] by:
2 iµ 2 1/2 1 u φα,µ := (u + 1+ u ) (1 + u ) 2F1 α, α + 1, 1+ iµ; + . − 2 2√1+ u2 p However, E. Wong did not provide any detail about the derivation of φα,µ and the correct expression is rather given by the following:
Proposition 1. We have 1 u φ (u) = (u + 1+ u2)iµ(1 + u2)α/2 F α, α + 1, 1+ iµ; + α,µ 2 1 − 2 √ 2 2 1+ u p 1 1+ α + iµ iµ α 1 (9) = 2F1 , − , 1+ iµ; . 2iµ(1 + u2)(iµ−α)/2 2 2 1+ u2 Proof. First of all, note that the quadratic transformation (see e.g. [14], p.112): c a c + a 1 (10) F (a, 1 a, c; u) = (1 u)c−1 F − , − , c; 4u(1 u) , 2 1 − − 2 1 2 2 − transforms the expression of φα,µ given by Wong into: 1 1+ α + iµ iµ α 1 φα,µ(u)= 2F1 , − , 1+ iµ; . 2iµ(1 + u2)(iµ−1)/2 2 2 1+ u2 4 L 1 As a matter of fact, it is natural to seek an eigenfunction of (1/2)−α,0 of the form :
1 1 φ (u)= f , y = , α,µ 1+ u2 1+ u2 for some smooth function f. In this respect, straightforward computations show that f satisfies:
4y2(1 y)f ′′(y)+ 4(α + 1)y (4α + 6)y2 f ′(y)= λf(y). − − − Now, assume further that f(y)= yβg(y),y (0, 1], ∈ for some power β R and a smooth function g. Then ∈ 4y2(1 y)g′′(y) + [8βy(1 y)+ y(4(α + 1) (4α + 6)y)]g′(y)= [λ+ − − − − (11) β(4(α + 1) (4α + 6)y) + 4β(β 1)(1 y)] g(y). − − − Next we take λ = α2 + µ2 and choose β as a root of
4β(β 1) + 4(α + 1)β + λ = 0, − that is β = ( iµ α)/2. If β = (iµ α)/2, then ± − − β(4α + 6) + 4β(β 1) = (iµ α)(α +1+ iµ) − − and 3 8β(1 y) + (4(α + 1) (4α + 6)y) = 4(iµ + 1) 4y iµ + . − − − 2 Keeping in mind (11) and the hypergeometric equation ([14]):
y(1 y)g′′(y) + (c (a + b + 1)y)g′(y)= abg(y), − − we get f(y)= φ (y). If β = (iµ + α)/2, then we similarly get α,µ − β(4α + 6) + 4β(β 1) = (iµ + α)(iµ α 1) − − − and 3 8β(1 y) + (4(α + 1) (4α + 6)y) = 4(1 iµ) + 4y iµ , − − − 2 − leading to φ = φ . Note that if β = (iµ α)/2, then the second linearly-independent solution α,−µ α,µ − of the hypergeometric equation which is regular around y = 0 is given by ([1], [14]):
1+ α iµ iµ α g(y)= y−iµ F − , − − , 1 iµ; y , 2 1 2 2 − (iµ−α)/2 so that f(y)= y g(y)= φα,−µ(y)= φα,µ(y).
1We discard the factor 2iµ. 5 2.1. The semi-group density and associated Legendre function on the cut. By the gen- eral expansion theorem for singular Sturm-Liouville operators ([31]), Chapter III), the semi-group density of the Hua-Pickrell process considered in [34] admits the following expansion: (12) M (α) (α) (α) 1 −n(2α−n)t Rn (v)Rn (u) 1 −(α2+µ2)t qt (v, u)= e + e φα,−µ(v)φα,µ(u)dµ , (1 + u2)α+(1/2) (α) 2 2π R (n=0 ( Rn 2) ) X || || Z where (α) 2 (α) (α) du ( Rn 2) := Rn (u)Rn (u) . || || R (1 + u2)α+(1/2) Z This squared L2-norm may be computed using the Rodrigues formula (7) and integration by parts as (see also Corollary 1 in [20]): πn! Γ(2α + 1 n) (13) ( R(α) )2 = − , || n ||2 22α−2n(α n) [Γ(α n + (1/2))]2 − − which agrees with equation (34) in [34]. (α) The semi-group density qt may be written as a single integral of associated Legendre functions on the cut ( 1, 1), known also as Ferrer functions ([14], eq. (6), p.143): − 1 1+ z µ/2 (14) P µ(z) := F ( α, 1+ α, 1 µ; (1 z)/2), z ( 1, 1). α Γ(1 µ) 1 z 2 1 − − − ∈ − − − This idea is actually motivated by the fact that (9) readily implies: u (15) φ (u)=(1+ u2)α/2Γ(1 + iµ)P −iµ α,µ α −√ 2 1+ u and similarly for φα,−µ(v). However, the connection between ‘symmetric’ Routh-Romanovski poly- nomials and Ferrer functions is not so obvious and we have not been able to find any reference where it is clearly stated. For the reader’s convenience, we prove it below using variable transformations of the Gauss hypergeometric function. Proposition 2. Let α 1. Then, for any n [α 1], ≥ ≤ − 2n−α√πΓ(1 + 2α n) v R(α)(v) = ( 1)n − (1 + v2)α/2P n−α . n − Γ(α + 1/2 n) α −√1+ v2 − Proof. Appealing again to (10), the Ferrer function may be written as: µ µ v 2 2 µ/2 1 µ + α µ + α 1 P = (1 + v ) 2F1 − , , 1 µ; , v R. α −√1+ v2 Γ(1 µ) 2 − 2 − 1+ v2 ∈ − Secondly, we apply the following quadratic transformation ([14], p. 112) : a a 1 z2 F , b , b + , = (1 z)a/2 F (a, b, 2b; z), 2 1 2 − 2 2 4(z 1) − 2 1 − with a = (µ + α)= n, b = (1/2) µ, z = 2/(1 iv) − − − − n−α 2 n/2 2 −n/2 2 α/2 n−α v 2 (1 + v ) 1+ v (1 + v ) Pα = −√1+ v2 Γ(1 + α n) (v + i)2 − 1 2 F + α n, n, 1 + 2α 2n; , 2 1 2 − − − 1 iv − 6 where for determinacy purposes we can assume that v 0 so that [(v + i)2]n/2 = (v + i)n. This ≥ (α) assumption is not a loss of generality since the orthogonality measure (8) is even so that Rn ( v)= n (α) n−α − ( 1) Rn (v) and since the same symmetry relation is satisfied by P ([14], p.144, (14)). Finally, − α we appeal to the following transformation: (b) 1 (16) n ( z)n F n, 1 c n, 1 b n; = F ( n, b, c, z) , (c) − 2 1 − − − − − z 2 1 − n which is readily checked by expanding the LHS and using the identity:
n−k (c)n (1 c n)n−k = ( 1) , 0 k n. − − − (c)k ≤ ≤ We apply it with b = (1/2) + α n, c =1+2α 2n to get: − − n 2n−α 2 α/2 n−α v n ( 1) 2 (1/2+ α n)n 1 1 iv (1 + v ) Pα = i − − 2F1 n,n 2α, α; − −√1+ v2 (1 + 2α 2n)nΓ(1 + α n) − − 2 − 2 − − ( 1)n22n−αΓ(1/2+ α)Γ(1 + 2α 2n) = in − − Γ(1 + 2α n)Γ(1 + α n)Γ(α n + 1/2) − − − 1 1 iv ( 1)n2αΓ(1/2+ α)n! F n,n 2α, α; − = − inJ (−α−1/2,−α−1/2)(iv), 2 1 − − 2 − 2 √π(1/2 α) Γ(1 + 2α n) n − n − where the last equality follows from Legendre duplication’s formula: 22z−1 1 Γ(2z +1) = Γ(z + )Γ(z + 1). √π 2 As a result, R(α)(v) = ( 2n)n!inJ (−α−1/2,−α−1/2)(iv) n − n 2n−α√π(1/2 α) Γ(1 + 2α n) v = − n − (1 + v2)α/2P n−α , Γ(α + 1/2) α −√ 2 1+ v and the lemma follows from the identity Γ(α + 1/2) (1/2 α) = ( 1)n . − n − Γ(α n + 1/2) −
Remark 1. When α 1 is an integer, the Ferrer function satisfies ([14], eq. (17), p. 144): ≥ v v Γ(2α n + 1)P n−α = ( 1)α−nΓ(α + 1)P α−n − α −√ 2 − α −√ 2 1+ v 1+ v so that ( 1)α2n−α√πΓ(α + 1) v R(α)(v)= − (1 + v2)α/2P α−n . n Γ(α + 1/2 n) α −√1+ v2 − Keeping in mind (13) and using Proposition 2, it follows that: (α) (α) Rn (v)Rn (u) Γ(2α + 1 n)(α n) 2 2 α/2 n−α u n−α v = − − [(1+u )(1+v )] Pα Pα . (α) 2 2 2 ( Rn ) n! −√1+ u −√1+ v || ||2 Besides, the discrete and continuous spectral values may be factorized as: α2 + µ2 = (α + iµ)(α iµ) n(2α n) = (α + n α)(α + α n). − 7 − − − Consequently, (α) Corollary 1. For any α> 0, the semi-group density qt may be written as: 2 α/2 (1 + v ) 2 2 v u q(α)(v, u)= e(α −µ )tP −µ P µ κ(dµ) t (1 + u2)(α+1)/2 α −√1+ v2 α −√1+ u2 Z where [α−1] Γ(2α + 1 n)(α n) Γ(1 µ)Γ(1 + µ) κ(dµ)= − − δ + − 1 R(µ)dµ. n! α−n 2iπ i nX=0 Remark 2. Since qα(v, ) is a probability density and since t · M R(α)(v)R(α)(u) du e−n(2α−n)t n n = 1, R (α) 2 (1 + u2)α+(1/2) n=0 ( Rn 2) Z X || || then we must have:
−µ2t du (17) e φα,−µ(v)φα,µ(u)dµ = 0, R R (1 + u2)α+(1/2) Z Z for all v R. In terms of Ferrer functions, the double integral (17) may be expressed as: ∈ 1 2 α/2 −µ2t −iµ v 2 iµ 2 α/2−1 (1 + v ) e Pα Γ(1 + iµ) dµ Pα (u)(1 u ) du, R −√1+ v2 | | − Z Z−1 which vanishes by the virtue of Lemma 2.13 in [10]. If µ = 0, α N 0 , then P 0 reduces to ∈ \ { } α a Legendre polynomial. In this case, (17) vanishes simply due to the orthogonality of Legendre polynomials with respect to Lebesgue measure on [ 1, 1] when α is even and due to their parity − otherwise. 2.2. Integer values of α: Wong’s second formula. In this case, Proposition (2) allows also to prove the ‘peculiar’ formula displayed in [34], equation (36), whose proof was ommitted. However, the computations are tedious and we do not detail them here since they do not lead to an easier formula for the semi-group density compared to the one obtained in Corollary 1. Nonetheless, the major steps towards Wong’s formula are as follows: Use (9) to expand φ (u) as: • α,µ α ( α) (1 + α) 1 + tanh(y) n φ (u) = [cosh(y)]αeiµy − n n , α,µ n!(1 + iµ) 2 n=0 n X where we set u := sinh(y), and do similarly for φα,µ(v) with v := sinh(z): α ( α) (1 + α) 1 + tanh(z) m φ (v) = [cosh(z)]αe−iµz − m m . α,−µ m!(1 iµ) 2 m=0 m X − Consider the integral: • −α2t 1 2 2 dµ e 2 dµ e−(α +µ )t = e−µ t 2π R (1 iµ) (1 + iµ) 2π R (m iµ) (1 iµ)(1 + iµ) (n + iµ) Z − m n Z − · · · − · · · and decompose m n 1 a (m,n) b (m,n) = k + k (m iµ) (1 iµ)(1 + iµ) (n + iµ) k iµ k + iµ − · · · − · · · k=1 − k=1 8 X X where ( 1)k−1k ( 1)k−1k a (m,n)= − , b (m,n)= − . k (m k)!(n + k)! k (n k)!(m + k)! − − Appeal to the integrals • ∞ 1 2 π 2 = e−(k±iµ)vdv, e−µ teiwµdµ = e−w /(4t), w R, k iµ R t ∈ ± Z0 Z r and to Fubini Theorem. Doing so gives rise to the incomplete Gamma integral ∞ 2 e−kve−(v±(y−z)) /(4t)dv, Z0 appearing in Wong’s formula. Reverse the summation order and arrange the different terms of the (n,m)-sums to get • a product of two Ferrer functions P k−α( tanh(y))P k−α( tanh(z)) to which Lemma 2 α − α − applies.
3. The HP diffusion with general parameters In this section, we consider the HP diffusion with general parameters A, K R (not necessarily ∈ stationary) whose generator is given by LA,K. We start with a review of its spectral resolution then prove an intertwining property between LA,K and L2−A,−K. In particular, this intertwining is deterministic in the sense that it does not involve a Markov kernel. Afterwards, we prove an identity (α) which maps a ‘symmetric’ Romanovski polynomial Rn to a ‘non-symmetric’ one recalled below. This identity follows from the Cauchy Beta integral and may be seen as a ‘random’ intertwining between LA,0 and LA,K. In the last section we notice that these intertwinings are closely related to some absolute continuity relations between the law of the HP diffusion with different parameters given by Girsanov’s Theorem.
3.1. Spectral resolution of the HP generator. The ‘non-symmetric’ Routh-Romanovski poly- nomials may be defined by (there is a misprint in [30], eq. (97)):
(α,K) n (−α−1/2−iK/2,−α−1/2+iK/2) Rn (u) := cni Jn (iu) ((1 iK)/2) α) 1 iK 1 iu = c in − − n F n,n 2α, − α; − . n n! 2 1 − − 2 − 2 They are eigenfunctions of L with eigenvalues n(2α n) and are finitely orthogonal with A,K − − respect to the weight:
−1 1 −1 W (u)e−K cot (u) = e−K cot (u). (1 + u2)α+1/2 Note in passing that the eigenvalues do not depend on K since the operator d K du lowers the polynomial degrees. Similar to LA,0, the HP generator LA,K admits a continuous spectrum of the form (α2 + µ2) corresponding to ([18], [19], [25]): − 1 iK 1 iu φ (u) := F α iµ, α + iµ, − α; − . α,K,µ 2 1 − − − 2 − 2 9 Therefore, the semi-group density of the HP diffusion admits an eigenfunction expansion similar to (12) which may be written as a single integral through the hypergeometric function: 1 iK 1 iu F α µ, α + µ, − α; − , µ n α, n [α 1] iR. 2 1 − − − 2 − 2 ∈ { − ≤ − }∪ 3.2. Intertwining. Recall the notations A 1/2= α and recall from [22], (paragraph 5.1) that − − the unitary operator: 1 1 f fW, W (u)= = , 7→ (1 + u2)(1/2)+α (1 + u2)1−A intertwines LA,0 and L2−A,0. This intertwining property extends to general parameters as follows: Lemma 1. For any A, K R and any smooth function f, ∈ −1 −1 eK cot eK cot (u) LA,K f (u)= [L2−A,−K(f)(u) + 2(1 A)] . W ! W (u) − Moreover, this is the only intertwining relation of this form.
Proof. Let f be a smooth function. Then Straightforward computations show that (1) acts on functions of the form −1 u f(u)(1 + u2)aeb cot (u), a, b R, 7→ ∈ as:
−1 (1 + u2)aeb cot (u) (1 + u2)∂2f + [(2A + 4a)u + K 2b]∂ f + 2af u − u 4a(A + a 1)u2 + 2[aK + b(1 A 2a)]u + b2 Kb + − − − − f 1+ u2 Hence, we readily see that the couple (a, b) = (1 A, K) annihilates the last term in which case − 2a = 2(1 A), 2A + 4a = 2(2 A), K 2b = K. One may also seek parameters (a, b) such that − − − − b2 Kb = 4a(A + a 1) = 0, aK = b(A + 2a 1). − − 6 − This forces in particular a = 0, b =0. If K = 0 then a = (1 A)/2 which lead to the contradiction: 6 6 − b2 = (A 1)2, − − since b takes real values. Otherwise, K = 0 and in turn A + 2a 1 = 0 so that 6 − 6 aK b = . A + 2a 1 − Consequently, (1 A)2 (A + 2a 1)2 a(A + a 1) b2 Kb = K2 − − − = K2 − = 4a(A + a 1), − 4(A + 2a 1)2 − (A + 2a 1)2 − − − which has no real solution as well.
According to this lemma, we may restrict the range of the parameters (A, K) to [1, [ R in ∞ × + our study of the HP diffusion. In Section 3.4 we shall see that the intertwining of Lemma 1 is ultimately related to the absolute continuity property between the laws of the corresponding HP processes by means of Girsanov’s Theorem. 10 3.3. From ‘symmetric’ to ‘non-symmetric’ Routh-Romanovski polynomials. Observe that (α,K) (α,0) (α) Rn is a deformation of Rn = Rn which adds ( iK/2, iK/2) to the parameters (1/2 α, 1/2 − − − α) of the Jacobi polynomials. On the one hand, this deformation reminds, yet in different from, the following identity valid for (c) > 0 (see e.g. (1.18) in [5]): ℜ J (a+c,a−c)(u) Γ(a + c + 1) 1 J (a,b)(u) (18) (1 u)a+c n = (1 y)a n (y u)c−1dy, u ( 1, 1), − (a+c,a−c) Γ(a + 1)Γ(c) − (a,b) − ∈ − Jn (1) Zu Jn (1) which follows from Euler’s first Beta integral. Remarkably, a similar fractional integral was the key ingredient used in [13] (see eq. (4.16)) for the derivation of the density of the exponential functional of Brownian motion (µ) (see Section 3.4 below). On the other hand, it is tempting At and quite interesting to seek an intertwining operator between HP generators corresponding to (A, K) and (A, 0) (see equation (20) in the next paragraph). In the next proposition, we prove an analogue of (18) relying this time on Cauchy’s Beta integral recalled below. Note that a variant of this integral was already used in the study of the Hua-Pickrell measure ([7]).
Proposition 3. For any n [α 1] and u R, ≤ − ∈
Γ(α n + 1/2) (α) dy (1 + iK/2)α−n (α,K) − lim Rn (y) = Rn (u). 2π w→−iu R (1 iy)α+1/2(w + iy)1+iK/2 (1 iu)α+(1+iK)/2 ℜ(w)>0 Z − −
(α) Proof. Apply the transformation (16) with b = n 2α, c = (1/2) α to write Rn as: − − (α) n (−α−1/2,−α−1/2) Rn (y)= cni Jn (iy) (n 2α) (1 iy)n 1 2 = c ( i)n − n − F n, + α n, 1 + 2α 2n; n − n! 2n 2 1 − 2 − − 1 iy n − c inΓ(2α n + 1) n (α n + 1/2) 2m−n = n − ( 1)m − m . n!Γ(2α 2n + 1) − m (2α 2n + 1) (1 iy)m−n − m=0 − m − Xn c inΓ(2α n + 1) n Γ(α n + 1/2+ m) 2m−n = n − ( 1)m − . n!Γ(α n + 1/2) − m Γ(2α 2n +1+ m) (1 iy)m−n m=0 − X − − Similarly,
(α,K) n (−α−1/2−iK/2,−α−1/2+iK/2) Rn (u)= cni Jn (iy) (n 2α) (1 iy)n 1+ iK 2 = c ( i)n − n − F n, + α n, 1 + 2α 2n; n − n! 2n 2 1 − 2 − − 1 iy n − c inΓ(2α n + 1) n (α n +(1+ iK)/2) 2m−n = n − ( 1)m − m n!Γ(2α 2n + 1) − m (2α 2n + 1) (1 iy)m−n − m=0 − m − X n c inΓ(2α n + 1) n Γ(α n + m +(1+ iK)/2) 2m−n = n − ( 1)m − . n!Γ(α n +(1+ iK)/2) − m Γ(2α 2n +1+ m) (1 iy)m−n m=0 − X − − Next, we apply Cauchy’s Beta integral ([4]): 1 dy Γ(a + b 1) = − , (w) > 0, (a + b) > 1, 2π R (1 iy)a(w + iy)b Γ(a)Γ(b)(1 + w)a+b−1 ℜ ℜ Z − 11 with a = α n + m + 1/2, b =1+ iK/2, to get − n Γ(α n + 1/2) (α) dy cni Γ(2α n + 1) − Rn (iy) = − 2π R (1 iy)α+1/2(w + iy)α n!Γ(1 + iK/2)(1 + w)α+(1+iK)/2 Z − n n Γ(α n + m +(1+ iK)/2) 2m−n ( 1)m − . − k Γ(2α 2n + m + 1) (1 + w)m−n m=0 X − Letting w iu, we are done. →−
3.4. Remarks on Girsanov’s Theorem and exponential functionals. In relation to a certain generalisation of Bougerol’s identity [32], the authors in [2] and [3] studied the diffusion whose infinitesimal generator acts on smooth functions as2: 1 d2 ν d + µ tanh(y)+ µ,ν R. 2 d2y cosh(y) dy ∈ With regard to (2), the corresponding diffusion is 1 K Y (A,K),t 0 , µ = A ,ν = , t/2 ≥ − 2 2 that we shall simply denote by Y (µ,ν), t 0. Recall from [2] that if P µ,ν denotes the probability t ≥ law on the canonical space C(R , R) of the diffusion Y (µ,ν), t 0, then the following absolute + t ≥ continuity relation holds
µ,ν µ,ν 0,0 (19) P Y = Dt (Y )P Y , |Ft |Ft where Y is the canonical coordinate process, P 0,0 denotes the Wiener measure and
t ν 1 t ν 2 Dµ,ν(Y ) = exp µ tanh (Y )+ dY ds µ tanh(Y )+ . t s cosh(Y ) s − 2 s cosh(Y ) Z0 s Z0 s ! We deduce from Itˆo’s formula that an intertwining operator between HP generators corresponding to (µ,ν) and (µ, 0) (equiv. (A, K) and (A, 0)) is related to the following absolute continuity relation (equation (20) corrects a typo in Section 3.1 of [2]): t 2 t µ,ν 1 2µ sinh(Ys) ν ds µ,0 (20) P Y = exp ν h(Yt) νh(Y0)+ ν − ds P Y , |Ft − 2 cosh2(Y ) − 2 cosh2(Y ) |Ft Z0 s Z0 s where x dy π h(x)= = tan−1(sinh(x)) = cot−1(sinh(x)). cosh(y) 2 − Z0 Note the occurrence of the hyperbolic scarf potential (5) in the RHS of (20). The constants ν(1 2µ)/2 and ν2/2 in equation (20) are invariant under the relation (µ,ν) − ↔ (1 µ, ν) and, therefore, for an intertwining operator between HP generators with parameters − − (µ,ν) and (1 µ, ν) (equiv. (A, K) and (2 A, K)) we further obtain the simplified relation: − − − − µ,ν (1/2−µ)t 1−µ,−ν (21) P Y = e W (sinh(Yt)) exp (2ν h(Yt)) P Y , |Ft |Ft
2 The parameter µ has no relation to the previous one encoding the continuous spectrum of LA.K . 12 where W (u) is given as in Section 3.2: 1 1 W (u)= = . (1 + u2)1/2−µ (1 + u2)1−A One readily sees that equation (21) is equivalent to the intertwining property of Lemma 1 and we can restrict the range of parameters (µ,ν) to [1/2, [ R . ∞ × + On the other hand, when µ < 1/2 and ν = 0 (equiv. A < 1 and K = 0), it is shown in [22], equation (4.2), that (21) is also equivalent to the identity:
µ2t/2E 1 1 (1−µ)2 t/2E 1 1 (22) e ψ = e ψ + 2γ1/2−µ , (µ) (µ) (1−µ) (1−µ) t ! t ! At A At A q q for every nonnegative Borel function ψ and t > 0, where (µ) is the exponential functional of At Brownian motion given by t (µ) (µ) := e2Bs ds, At Z0 and γ1/2−µ is a Gamma variable (independent of B) with density: e−u u−(1/2+µ) 1 . Γ(1/2 µ) {u>0} − Note that the equation (22) is reminiscent of Dufresne’s identity (see e.g. [13]); however, the additional square root factor of 1/ (µ) is related to the intertwining (µ, 0) (1 µ, 0) (equiv. At ↔ − (A, 0) (2 A, 0)), instead of the “opposite drift” case (µ, 0) ( µ, 0) considered e.g. in [21]. ↔ − ↔ − We remark that it would be interesting to find a (two-dimensional) analogue of the identity (22) corresponding to the general case ν = 0 (equiv. K = 0) in equation (21). In this regard, since for 6 6 fixed t> 0 we have ([2], Proposition 1)
t (µ) (µ,ν) d (µ) (µ,ν) (µ) Bs sinh(Yt ) = β (µ) + νat , with Y0 = 0, at := e ds, At Z0 where B and β are independent Brownian motions, one easily obtains from equation (21) that:
2 − 2 − 1 − u−νa(µ) /2A(µ) 1 − u+νa(1 µ) /2A(1 µ) (23) E e t t = e(1/2−µ)tη(u) E e t t , (µ) (1−µ) 2π 2π At At q −1 q a.s. u R, where η(u) := W (u)eν(π−2 cot (u)) (see also Lemma 1). The analysis of equation (23) ∈ is related to the joint distribution of the vector (see [35], Section 6)
t (µ) (µ) d Bs (µ) β (µ) , a = e dγ , a , t> 0, A t s t t 0 Z where γ is another Brownian motion independent of B. Moreover, if µ< 0 (equiv. A< 1/2), note that the HP process sinh(Y µ,ν), t 0, converges in law, as t , to the so-called ‘subordinated t ≥ → ∞ perpetuity’:
∞ (µ) Bs (ν) d (µ) (ν) e dγ = β (µ) + νa , where γ = γs + νs, s A∞ ∞ s Z0 and whose (unnormalized) probability density is given by the function η(u) defined above (see [21], Section 4). 13 4. Integral representation of the HP semi-group density The left-hand side of Equation (23) gives a representation formula for the HP semi-group density; however, this is not very tractable. In this section we prove an explicit integral representation through an extension of the Maass semi-group density to imaginary values of the magnetic field. With the same notation as in Section 3.4, the following Fourier transform of the HP process was derived in [3], Proposition 6.1. If Y (µ,ν) = v R: 0 ∈ ∞ ∞ E (µ,ν) −µ2t/2 y (24) v exp iλ sinh Yt = e dz dy exp(iλ sinh(v)e + µy + iλνz) 0 −∞ h i Z Z
λ t exp( λ(1 + ey) coth(λz/2))θ , 4 sinh(λz/2) − φ(z,y;λ) 4 where ∞ r 2 2 πw (25) θ (t) := eπ /2t e−w /(2t)e−r cosh(w) sinh(w)sin dw, r 3 1/2 (2π t) 0 t Z and 2λey/2 φ(z,y; λ) := . sinh(λz/2) N Note that the double integral in (24) is absolutely convergent for all µ since θr(t) = O(r ) as r gets close to zero, for any N > 0 (see [22], p. 188 and paragraph 6.2 there). In order to state our result, recall from [3] (see also [17]) the Maass Laplacian: y2 k2 H := (∂2 + ∂2)+ iky∂ + , k − 2 w y w 2 where k R is interpreted in the physics realm as the intensity of a constant magnetic field. This ∈ is a densely-defined and essentially-self-adjoint operator in dy L2 H, dw , y2 where H is the Poincar´eupper-half-plane, and reduces to the Laplace-Beltrami operator in H when k = 0: y2 H = (∂2 + ∂2). 0 − 2 w y The semi-group density of its self-adjoint closure with initial point z = i admits the following expression ([3], [17]):
2 w + i(y + 1) k √2e−t/8−k t/2 (26) QH (i, w + iy)= t,k i(y + 1) w (2πt)3/2 − 2 ∞ ze−z /(2t) cosh(z/2) dz cosh 2k cosh−1 . r cosh(z) cosh(r) cosh(r/2) Z − with respect to the hyperbolic volume measurep dwdy/y2. Here the principal determination of the power is taken and r = dH(i, w + iy) is the hyperbolic distance: w2 + (y 1)2 cosh(r)=1+ − . 2y 14 With these notations, we have the following:
Theorem 1. The semi-group density of the HP process sinh Y (µ,ν) , t 0, admits the following t ≥ integral representation: ∞ 2 2 g(µ,ν)(sinh(v), w)= e(1−4µ )t/8e−ν t/2 dy yµ−3/2QH (i, w sinh(v)y + iy), w, v R, t t,iν − ∈ Z0 H where Qt,iν is the extension of the RHS of (26) to purely imaginary values k = iν. Proof. Perform the variable change z 2z/λ in the double integral displayed in (24) and change 7→ the order of integration there to write: ∞ E (µ,ν) −µ2t/2 y (27) v exp(iλ sinh(Yt )) = e dy exp(iλ sinh(v)e + µy) −∞ h i Z ∞ e2iνz t dz exp( λ(1 + ey) coth(z))θ . 2 sinh(z) − φ(2z/λ,y;λ) 4 Z0 On the other hand, Proposition 4.1 and equation (5.2) in [3] together with Fourier inversion Theorem show that for any λ> 0, k R, ∈ ∞ e2kz t dz exp( λ(1 + ey) coth(z))θ = qM (0,y) 2 sinh(z) − φ(2z/λ,y;λ) 4 t,λ Z0 t/8+k2t/2 −y/2 iλw H y (28) = e e dwe Qt,k(i, w + ie ), R Z M where qt,λ(0,y) is the semi-group density of the Schr¨odinger operator with Morse potential: 1 d2 λ2 HM := + e2u, λ −2 du2 2 starting at the origin. Now, the LHS of (27) is analytic in k since the integral converges absolutely for any k R. Besides, we can rewrite QH (i, w + iy) as ∈ t,k 2 √2e−t/8−k t/2 QH (i, w + iy)= eik(2 arg(w+(y+1)i)−π) t,k (2πt)3/2 2 ∞ ze−z /(2t) cosh(z/2) dz cosh 2k cosh−1 , r cosh(z) cosh(r) cosh(r/2) Z − which extends analytically to k C. Since thep integral displayed in the RHS of (28) still converges ∈ absolutely for k C, then (28) extends analytically to the whole complex plane. In particular, for ∈ k = iν, the inner integral in (27) is expressed as:
∞ 2iνz e y t t/8+k2t/2 −y/2 iλw H y dz exp( λ(1+e ) coth(z))θφ(2z/λ,y;λ) = e e dwe Qt,iν(i, w+ie ) 2 sinh(z) − 4 R Z0 Z whence ∞ E (µ,ν) (1−4µ2)t/8 −ν2t/2 iλ sinh(v)y µ−3/2 iλw H v exp(iλ sinh(Yt )) = e e dye y dwe Qt,iν(i, w + iy). 0 R h i Z Z The double integral in the right-hand side converges absolutely for any µ R. Indeed, ∈ H ν2t/2 2νπ H (29) Qt,iν(i, w + iy) e e Qt,0(i, w + iy) | |≤ 15 and QH (i, w + iy) is the semi-group density of the hyperbolic Brownian motion (N , O ) H t,0 t t t≥0 ∈ solution of the following SDE:
W 1−t/2 Ot = e t , t 2 Nt = OsdWs , Z0 1 2 for a planar (Euclidean) Brownian motion (Wt ,Wt )t≥0. Thus, ∞ dy µ+1/2 H ν2t/2 2νπE µ+1/2 y dw Qt,iν (i, w + iy) e e (Ot ) < . y2 R | |≤ ∞ Z0 Z As a matter, we can invert order of integration and get the semi-group density of the HP diffusion sinh Y (µ,ν) , t 0. The theorem is proved. t/2 ≥ Remark 3. From Remark 3.1. in [17], it follows that for any k R: 3 ∈ s H dwdy [s(s−1)−k2]t/2 y Qt,k(i, w + iy) = e . H y2 Z By the same argument used in the previous proof, this equality extends to all complex values of k and shows in particular, with s = µ + 1/2 and k = iν, that ∞ 2 2 w e(1−4µ )t/8e−ν t/2 dyyµ−3/2QH (i, w sinh(v)y + iy), 7→ t,iν − Z0 is indeed a probability density. On the other hand, the intertwining relation proved in Lemma 1 yields another integral representations:
2ν cot−1(sinh(v)) 2 1/2−µ (µ,ν) (1−4µ2)t/8 −ν2t/2 e (1 + (sinh (v)) g (sinh(v), w)= e e − t e2ν cot 1(w)(1 + w2)1/2−µ ∞ dy y−µ−1/2QH (i, w sinh(v)y + iy). t,−iν − Z0 In particular, both representations coincide for µ = 1/2,ν = 0.
5. The Hua-Pickrell particle system In this final section we give some interest in the multivariate analogue of the HP process: the HP particle system. More precisely, let s be a complex parameter such that (s) > 1/2, the HP N ℜ − particle system X = (Xn)n=1 was introduced in [6] as the N-dimensional diffusion satisfying the stochastic differential system:
2 2 1 + (Xn(t)) dXn(t)= 2(1 + (Xn(t)) dWn(t) + 2 [(1 N (s)]Xn(t)+ (s)] + dt, − −ℜ ℑ Xn(t) Xj(t) j6=n p X − 1 n N, ≤ ≤
3The eigenvalue is λ = [−s(s − 1) + B2]/2 and not [−s(s − 1)/2+ B2]/2. 16 N where (Wn)n=1 is a N-dimensional Brownian motion. It was also shown there (see Lemma 4.2) (s,N) that the particles never collide almost surely and that the multivariate HP measure τHP : N τ (s,N)(dx ,...,dx ) := C [V (x)]2 (1 + x2 )−ℜ(s)−N e2ℑ(s) arg(1+ixn)dx , x = (x )N RN , HP 1 N s,N n n n n=1 ∈ nY=1 where Cs,N is a normalizing constant and V (x)= (x x ) n − j 1≤n
Proposition 4. For any real s> 1/2, the semi-group density of the HP-particle system is given − by: n 2 (N+s−(1/2))/2 (s,N) [(N+s−1/2)2−λ ]t (1 + v ) 2 (−µ2+···+µ2 )t G (x,y)= e s,N V (y) e 1 N t (1 + u2)(N+s+(1/2))/2 n=1 µ1>···>µN Y Z N N x y det P −µj n det P −µj n N+s−1/2 2 N+s−1/2 2 N − 1+ xn !! − 1+ yn !! j,n=1 j,n=1 κ(dµ ). V (x) V (y) n p p n=1 17 Y Proof. It suffices to write n 2 (N+s−(1/2))/2 N (s,N) −λs,N t (1 + v ) V (y) (N+s−1/2) Gt (x,y)= e det qt (xn,yj (1 + u2)(N+s+(1/2))/2 V (x) n,j=1 n=1 Yn (1 + v2)(N+s−(1/2))/2 V (y) = e−λs,N t (1 + u2)(N+s+(1/2))/2 V (x) nY=1 N 2 2 x y det e[(N+s−1/2) −µ ]tP −µ n P µ j κ(dµ) . N+s−1/2 − 1+ x2 N+s−1/2 − 2 Z n ! 1+ yj n,j=1 p q
Remark 4. The function N x det P −µj n N+s−1/2 2 − 1+ xn !! x j,n=1 7→ V (xp) is a symmetric multivariable extension of the associated Legendre function on the cut.
References
[1] G. E. Andrews, R. Askey, R. Roy. Special Functions. Cambridge University Press, Cambridge. (1999). [2] L. Alili, D. Dufresne, M. Yor. Sur l’identiti´ede Bougerol pour les fonctionnelles exponentielles du mouvement Brownien avec drift. [3] L. Alili, H. Matsumoto, T. Shiraishi. On a triplet of exponential Brownian functionals. S´eminaire de Probabilit´es. 35, (2001), 396-415. [4] R. Askey. An integral of Ramanujan and orthogonal polynomials, J. Indian Math. Soc. 51, (1987), 27-36. [5] R. Askey, J. Fitch. Integral representations for Jacobi polynomials and some applications. J. Math. Anal. Appl. 26. (1969), 411-437. [6] T. Assiotis. Hua-Pickrell diffusions and Feller processes on the boundary of the graph of spectra. Ann. de l’Institut Henri Poincare (B): Probabilites et Statistiques, Vol. 56, no. 2, 1251-1283 (2020). [7] A. Borodin, G. Olshanski. Infinite Random Matrices and Ergodic Measures. Comm. Math. Phys. 223, (2001), 87-123. [8] A. Bufetov, Y. Qiu. The explicit formulae for scaling limits in the ergodic decomposition of infinite Pickrell measures, Arkiv f¨or Matematik, 54 (2016), no.2, 403-435. [9] P. Bourgade, A. Nikeghbali, A. Rouault. Circular Jacobi ensembles and deformed Verblunsky coefficients. Int. Math. Res. Not. IMRN (2009), no. 23, 4357-4394. [10] H. S. Cohl, T. H. Dang, T. M. Dunster. Fundammental solution and Gegenbauer expansions of Helmholtz operators in Riemannian spaces of constant curvature. SIGMA. [11] J. W. Dabrowska, A. Khare, U. P. Sukhatme. Explicit wavefunctions for shape-invariant potentials by operator techniques. J. Phys. A: Math. Gen. 21, (1988), 195-200. [12] P. Deift, D. Gioev. Random matrix theory: invariant ensembles and universality. Courant Lecture Notes in Math- ematics, 18. Courant Institute of Mathematical Sciences, New York; American Mathematical Society, Providence, RI, 2009. [13] D. Dufresne. The integral of geometric Brownian motion. Adv. Appl. Prob. 33, (2001), 223-241. [14] A. Erdelyi, W. Magnus, F. Oberhettinger, F. G. Tricomi. Higher transcendental functions. Vol. I. McGraw-Hill Book Company, Inc., New York-Toronto-London, 1953. xvii+396 pp. [15] J. C. Gruet. Semi-groupe du mouvement Brownien hyperbolique. Stochastics Stochastic Rep. 56, (1996), 53-61. [16] G. Heckman, G. Schlichkrull. Harmonic Analysis and Special Functions on Symmetric spaces. Perspectives in Mathematics. 16, (1994). 18 [17] N. Ikeda, H. Matsumoto. Brownian motion of the hyperbolic plane and Selberg trace formula. J. Funct. Anal. 163, (1999), 63-110. [18] A. Khare, U. P. Sukhatme. Scattering amplitudes for supersymmetric shape-invariant potentials by operator methods. J. Phys. A: Math. Gen. 21, (1988), 501-508. [19] G. L´evai, F. Cannata, A. Ventura. Algebraic and scattering aspects of a PT-supersymmetric solvable potential. J. Phys. A: Math. Gen. 34, (2001), 839-844. [20] M. Masjed-Jamei, F. Marcellan, E. J. Huertas. A finite class of orthogonal functions generated by Routh- Romanovski polynomials. Complex Var. Elliptic Equ. 59 (2014), no. 2, 162-171. [21] H. Matsumoto, M. Yor. A relationship between Brownian motions with opposite drifts via certain enlargements of the Brownian filtration. Osaka J. Math.. 38, (2001), 383-398. [22] H. Matsumoto, M. Yor. On Dufresne’s relation between the probability laws of exponential functionals of Brow- nian motions with different drifts. Adv. Apple. Probab. 35, No. 1. (2003), 184-206. [23] H. Matsumoto, M. Yor. Exponential functionals of Brownian motion, I: Probability laws at fixed time . Probability Surveys. 2, (2005), 312-347. [24] I. G. MacDonald. Symmetric Functions and Hall Polynomials. Second edition, Mathematical Monographs, Oxford. 1995. [25] Yu. A. Neretin. Some Continuous Analogs of the Expansion in Jacobi Polynomials and Vector-Valued Orthogonal Bases. Funct. Anal. Appl. 39, (2005), 106-119. [26] J. Peetre. Correspondence principle for the quantized annulus, Romanovski polynomials and Morse potentials. J. Funct. Analysis. 117, (1993), 377-400. [27] M. V. Romanovski. Sur quelques classes nouvelles de polynˆomes orthogonaux. C. R. A. S. 188. (1929), 1023-1025. [28] E. J. Routh. On some properties of certain solutions of a differential equation of second order. Proc. London Math. Soc. 16: 245. (1884). [29] J. W. Pitman, L. C. G. Rogers. Markov functions. Ann. Probab. Vol. 9, No. .4. (1981). 573-582. [30] A. Raposi, H.J. Weber, D. Alvarez-Castillo, M. Kirchbach. Romanovski polynomials in selected physics problems. Cent. Eur. J. Phys. 5 (3), (2007), 253-284. [31] E. C. Titchmarsh. Eigenfunction Expansions associated with Second-order Differential Equations. Part I. Second Edition, Oxford University Press, 1962. [32] S. Vakeroudis. Bougerol’s identity in law and extensions. Probab. Surv. 9. (2012), 411-437. [33] J. F. van Diejen. A combinatorial formula for the associated Legendre functions of integer degree. Adv. Maths. 149. (2000), 61-88. [34] E. Wong. The construction of a class of stationary Markoff processes. (1964), Proc. Sympos. Appl. Math., Vol. XVI pp. 264-276. Amer. Math. Soc., Providence, R.I. [35] M. Yor, Interpretations in terms of Brownian and Bessel meanders of the distribution of a subordinated perpe- tuity, L´evy processes. Birkh¨auser Boston, Boston, MA, (2001). 361–375. [36] M. Yor, On some exponential functionals of Brownian motion, Adv. Appl. Prob. 24, (1992). 509-531.
School of Mathematics and Statistics, University College Dublin, Dublin 4, Ireland Email address: [email protected]
IRMAR, Universite´ de Rennes 1, Campus de Beaulieu, 35042 Rennes cedex, France Email address: [email protected]
19