Next Article in Journal
Sobolev Estimates for the ¯ and the ¯-Neumann Operator on Pseudoconvex Manifolds
Previous Article in Journal
Robust Bi-Level Optimization for Maritime Emergency Materials Distribution in Uncertain Decision-Making Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Monotonicity and Infinite Divisibility Properties for the Mittag-Leffler Function and for Stable Distributions

Department of Statistics & OR, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 2 September 2023 / Revised: 21 September 2023 / Accepted: 25 September 2023 / Published: 30 September 2023
(This article belongs to the Section Probability and Statistics)

Abstract

:
Hyperbolic complete monotonicity property ( HCM ) is a way to check if a distribution is a generalized gamma ( GGC ), hence is infinitely divisible. In this work, we illustrate to which extent the Mittag-Leffler functions  E α , α ( 0 , 2 ] , enjoy the  HCM property, and then intervene deeply in the probabilistic context. We prove that for suitable  α and complex numbers z, the real and imaginary part of the functions  x E α z x , are tightly linked to the stable distributions and to the generalized Cauchy kernel.

1. Introduction and First Results

The Mittag-Leffler function,
E a ( z ) = k 0 z k Γ ( 1 + k a ) , z C , a > 0 ,
is a widely studied special function; see [1] and the references therein. The essence of this work is to unveil several probabilistic interpretation for the functions  E α , α ( 0 , 2 ] . To this end, we need some setting.

1.1. Bernstein Functions and Infinite Divisibility

A function  f : ( 0 , ) ( 0 , ) is completely monotone (we denote  f CM ) if it is infinitely differentiable and it satisfies  ( 1 ) n f ( n ) ( x ) 0 , for all n N , x > 0 . By Bernstein characterization,  f CM if and only if, it is the Laplace transform of some measure
f ( λ ) : = L τ ( λ ) = [ 0 , ) e λ x τ ( d x ) , λ > 0 .
A subset of  CM is the class  S a of generalized Stieltjes transforms of order  a > 0 , viz. of functions f represented by
f ( λ ) = d + q λ a + ( 0 , ) 1 ( λ + x ) a U ( d x ) , λ > 0 ,
where  q , d 0 and U is a Radon measure on  ( 0 , ) such that  ( 0 , ) ( 1 + u ) a U ( d u ) < . Stieltjes transforms of order 1 are usually called Stieltjes functions. The class  BF of Bernstein functions consists of functions of the form
ϕ ( λ ) = q + d λ + ( 0 , ) ( 1 e λ x ) Π ( d x ) = q + d λ + λ ( 0 , ) e λ x Π ( x , ) d x , λ 0 ,
where  q 0 is the so-called killing rate d R is the drift, and  Π is the Lévy measure of  ϕ , i.e., a positive measure on  ( 0 , ) , which satisfies  ( 0 , ) x 1 Π ( d x ) < . Two important subclasses of  BF are the ones of Thorin Bernstein functions and of complete Bernstein functions, denoted by  TBF and  CBF , respectively. A function  ϕ belongs to  TBF (respectively,  CBF ) .... if it is represented, for all  λ 0 , by
ϕ ( λ ) = q + d λ + ( 0 , ) log 1 + λ u U ( d u ) , resp . q + d λ + ( 0 , ) λ λ + u V ( d u ) ,
where the positive measures U and V satisfy
( 0 , 1 ) | log u | U ( d u ) + [ 1 , ) u 1 U ( d u ) < ,
and V ( d u ) u is a   Lévy measure .
Note the equivalences
ϕ BF ( respectively , TBF , CBF ) ϕ 0 and ϕ CM ( respectively , S 1 , S 2 ) .
The book of Hirschman and Widder [2] is a good reference for the classes  S a and, especially for  S 1 , we refer to [3] (Chapter 2). The books of Schilling, Song, and Vondraçek [3] and the one of Steutel and van Harn [4] are our main references for the class  BF and its subclasses, and for infinitely divisible measures as well. The distribution of a nonnegative random variable X is said to be infinitely divisible, and we denote  X I D , if there exists an i.i.d. sequence  ( X i n ) 1 i n , such that  X = d X 1 n + + X n n . The celebrated Lévy-Khintchine formula gives the following characterization through the cumulant function of  X 0 :
X I D E [ e λ X ] t = [ 0 , ) e λ x P ( X d x ) t CM , t > 0 ϕ X ( λ ) : = log E [ e λ X ] BF .

1.2. The HCM Property and GGC Distributions

The class of infinitely divisible distributions behind  TBF (respectively,  CBF ) is known as the generalized gamma convolutions GGC (respectively, the Bondesson class  BO ), and it corresponds to the smallest class of sub-probability measures on  [ 0 , ) which contains mixtures of gamma (respectively, exponential) distributions and which is closed under convolutions and vague limits. See [3] (Theorem 9.7 and Proposition 9.11). These classes were introduced by Olaf Thorin, and were widely developed by Lennart Bondesson in [5]; see also [3,4]. We will now introduce an important subclass of  GGC . A function  f : ( 0 , + ) ( 0 , + ) is said to be hyperbolically completely monotone, and we denote  f HCM if, for every  u > 0 , the function
H u ( w ) : = f ( u v ) f ( u / v ) ,
is completely monotone in the variable w = v + 1 / v 2 ,
(it is easy to see that  H u is always a function of w). In [5] (Theorem 5.3.1), it is shown that the class  HCM corresponds to pointwise limits of functions of the form
x C x a 1 i = 1 N ( x + c i ) t i , N N , C , c i , t i > 0 , 0 < a < i = 1 N c i .
Property [5] ((iv) p. 68) asserts that
f HCM x a f ( x b ) HCM , for all a R , and | b | 1 .
By Bondessons’s definition [3] (Definition 9.10), and by [3] (Theorem 6.1.1), we have
X GGC ϕ X TBF e ϕ X HCM .
The class of  HCM functions are stable under multiplicative convolution [5] (property vii, p. 68), then the Laplace transform of an  HCM function is itself  HCM . If the p.d.f. of a continuous and positive r.v. X is  HCM , we denote  X HCM . Clearly, we have the classification
HCM GGC BO , GGC SD I D ,
where  SD denotes the well-known class of self-decomposable distributions; see [4] for this class. Note that all previous classes of distributions are stable by additive convolution. The importance of the classes  HCM and  GGC is due to their specific stability properties, one can find in [6]: if X and Y are independent and  X , Y GGC (respectively  HCM ), then
X + Y , X × Y GGC ( resp . X q , q > 1 , X × Y , X / Y HCM ) .
Note that the above stability properties are not shared by general distribution of the half real line. Moreover, by (7), the  I D property reads only on the level of the cumulant function, whereas the  HCM one reads on the level of the probability density function, cf. the examples of next subsection. For a better comprehension of the above classes, see Figure 1.

1.3. Gamma and Positive Stable Distributions

In what follows,  G t denotes a standard gamma distributed random variable with shape with parameter  t > 0 , i.e., with Laplace transform, Mellin transform, and p.d.f.
E [ e λ G t ] = 1 1 + λ t , λ 0 , E [ G t s ] = Γ ( t + s ) Γ ( t ) s > t , f G t ( x ) = x t 1 Γ ( t ) e x , x > 0 .
The last p.d.f. form and (8) ensure that all gamma distributions are  HCM . Another example of a  GGC distribution is given by the positive stable r.v.  S α , α ( 0 , 1 ) , which is associated to the Thorin–Bernstein function
ϕ S α ( λ ) = log E [ e λ S α ] = log 0 e λ x f S α ( x ) d x = λ α = α sin ( α π ) π 0 log 1 + λ x d x x 1 α , λ 0 .
Note that the p.d.f.  f S α is not explicit except for  α = 1 / 2 ;
S 1 / 2 = d 1 4 G 1 / 2 and f S 1 / 2 ( x ) = 1 2 π x 3 e 1 / 4 x , x > 0 ,
nevertheless, the Mellin transform of  S α is explicitly given by
E S α s = Γ ( 1 + s α ) Γ ( 1 + s ) , s > α ,
which gives the convergence in distribution
1 ( S α ) α d 1 , as α 1 , and 1 ( S α ) α d G 1 , as α 0 .
It is then natural to adopt the conventions  S 0 = 0 , S 1 . Related literature for the stable distribution is referenced on Nolan’s website [7]. The monograph of Zolotarev [8] is one of the major references for general stable distributions. By representation (13), and by the main result of [9], we know that
S α GGC , for all α ( 0 , 1 ) , and S α HCM α 1 / 2 .
Additionally, in [10], it was shown that if  t > 0 , then
S α t I D t 1 α .
Using (16), and taking  S α an independent version of  S α , let us define the quotient,
T α : = S α S α α , α ( 0 , 1 ) , T 1 = 1 , and T 0 = d G 1 G 1 ,
where  G 1 , G 1 are i.i.d. and standard exponentially distributed. The p.d.f. of  T α is explicitly given by the generalized Cauchy form
f T α ( x ) = sin ( π α ) π α 1 1 + 2 cos ( π α ) x + x 2 , x > 0 ,
see (75) below for instance. Bosch [11] completed (18) by showing
α ( 0 , 1 / 2 ] and | t | 1 α T α 1 / t = S α S α α / t HCM .
With the above information, Bosch and Simon found it natural to raise the following open question in [9]:
S α α / ( 1 α ) HCM S α α / t HCM   if ,   and   only   if ,   α 1 / 2   and   | t | 1 α ?
This open question extends Bondesson’s conjecture [5], which was stated in 1977. This conjecture, which is identical to (22) with  t = α , was investigated in [10] and solved in [9].

1.4. First Results on the Mittag-Leffler Functions

From (15), one obtains a first link of the Mittag-Leffler function with the positive stable distribution:
E e z / ( S α ) α = E α ( z ) , α ( 0 , 1 ) , z C .
By (13), it is clear that
E [ e λ T α 1 / α ] = E [ e λ S α / S α ] = E [ e ( λ / S α ) α ] = E α ( λ α ) , λ 0 .
In [10], if was shown that
0 < α 1 2 T α 1 / α HCM E α ( λ α ) HCM .
Note that the r.v.
R α : = 1 + T | α 1 | 1 / α , α [ 0 , 2 ] ,
appeared in [12] ((184) and (223)), in a  GGC context without noticing its  GGC property: by (21), we have
1 2 α 1 T | α 1 | 1 / α = R α 1 HCM R α GGC .
Another monotonicity property of the Mittag-Leffler functions can be found in [13]:
0 < α 1 2 λ E α ( λ ) BF 1 Γ ( 1 α ) λ E α ( λ ) CM ,
and results in case  α ( 1 / 2 , 1 ) are also shown there. Note that the function  E α ( λ ) is completely monotone, contrary to  E α ( λ ) . In [14] (Theorem 1.2, (a) and (c)), Simon showed that
λ F α ( λ ) : = e λ α E α ( λ α ) 1 α CM , for all α ( 0 , 1 ) ( 1 , 2 ] .
The choice of the function  F α is intuitive because  e λ = E 1 ( λ ) = lim α 1 α E α ( λ α ) .
The main objective of this work is to exhibit new monotonicity properties of the  GGC and  HCM type for the Mittag-Leffler functions. To give a taste of our results, we start by improving Simon’s results (27).
Theorem 1.
For  α ( 0 , 1 ) ( 1 , 2 ] , the following holds.
(1) 
With  T | α 1 | given by (19), and  F α by (27), we have the Laplace representation:
F α ( λ ) : = E [ e λ T | α 1 | 1 / α ] , λ 0 .
(2) 
Assume  α ( 0 , 1 ) .
(a) 
The functions  1 α e λ E α ( λ α ) and  e λ / E α ( λ α ) are completely monotone, and  e λ E α ( λ α ) 1 is Bernstein.
(b) 
The function  F α is Stieltjes. Moreover,  F α HCM , if and only if,  1 / 2 α < 1 . In this case,  1 α e λ E α ( λ α ) HCM and  F α / F α S 1 .
(3) 
Assume  α ( 1 , 2 ] .
(a) 
The function  α e λ E α ( λ α ) 1 is completely monotone and  1 e λ E α ( λ α ) is Bernstein.
(b) 
The function  F α is not  HCM if  α > 4 / 3 .
(c) 
If  1 < α 3 / 2 and  0 < γ ( 2 α ) / α , then  F α ( λ γ ) is  HCM , and both functions  F α ( λ γ ) and  λ γ 1 F α ( λ γ ) / F α ( λ γ ) are in  S 1 .
In the same direction as (27), Simon considered in [14,15] the function
G α ( x ) : = E α ( x α ) d d x E α ( x α ) = 1 α α F α ( x ) F α ( x ) , x > 0 .
From (25) and (28), we deduce the expressions
G α ( x ) = α 1 α E 1 + T | α 1 | 1 / α e x T | α 1 | 1 / α , α ( 0 , 1 ) ( 1 , 2 ] .
and
α α 1 e x G α ( x ) = E R α e x R α
is the p.d.f. of the independent quotient G 1 / R α .
Using (15), we see that
α ( 1 , 2 ] E T α 1 1 / α = 1 α 1 ,
thus
G α ( 0 + ) = α 1 α E 1 + T α 1 1 / α = 1 = 0 G α ( x ) d x ,
which means that  G α , α ( 1 , 2 ] , is a p.d.f. and is also the Laplace transform of a probability distribution; see (69) below. In case  α ( 0 , 1 ) , we have  E [ T 1 α 1 / α ] = , and the function  G α is not a p.d.f. We complete our previous discussion with the following result.
Theorem 2. 
The r.v.  R α and the function  G α , given in (29) and (25), satisfy the following.
(1) 
If  α ( 1 , 2 ] , then the p.d.f.  G α is the one of a  GGC distribution but not a  HCM .
(2) 
If  α ( 0 , 1 ) , then the function  G α is a widened  GGC , if and only if,  α 1 / 2 . In this case,  G 1 / R α GGC .
In [16], Simon was interested in the first passage time of the normalized spectrally positive  α -stable Lévy process  ( S t ( α ) ) t 0 , α ( 1 , 2 ] , viz.
E [ e λ S t ( α ) ] = e t λ α , λ , t 0 ,
and
τ ^ x = inf { t > 0 ; S t ( α ) = x } and τ x = inf { t > 0 ; S t ( α ) = x } , x > 0 .
The self-similarity property of  ( S t ( α ) ) t 0 with index  1 = 1 / α entails
S t ( α ) = d t 1 / α S 1 ( α ) , ( τ ^ x , τ x ) = d x α ( τ ^ 1 , τ 1 ) .
Additionally, it is well known that
τ ^ 1 = d inf { t > 0 ; S t ( α ) > 1 } = d S 1 / α ,
and this explains our focus on  τ 1 . To formalize our next results, we need to recall the biasing procedure for the distribution of a non-negative r.v. Z. For  u R such that  E [ Z u ] < , we denote by  Z [ u ] a version of the length-biased distribution of order u of Z, viz.
P ( Z [ u ] d x ) = d x u E [ Z u ] P ( Z d x ) .
We complete Simon’s main result in [16] by explaining (in (37) below) his factorization of  τ 1 obtained in his main theorem, and by the use of the function
H α ( λ ) : = ( 1 α ) F α ( λ 1 / α ) = e λ 1 α α 2 λ 1 1 α E α ( λ ) , α ( 1 , 2 ] .
Corollary 1.
Let  α ( 1 , 2 ] . The r.v.  τ 1 given by (33), and  H α by (35), satisfy the following.
(1) 
We have the Laplace representation:
E [ e λ τ 1 ] = H α ( λ ) , λ 0 .
(2) 
With the convention (19), we have the independent factorizations in law:
τ 1 = d S 1 / α ( T α 1 ) [ 1 / α ] a n d G 1 τ 1 1 / α = d G 1 ( T α 1 1 / α ) [ 1 ] ,
where the distributions of  ( T α 1 ) [ 1 / α ] and  ( T α 1 1 / α ) [ 1 ] are given by the biasing procedure (34).
(3) 
( G 1 / τ 1 ) 1 / ( 2 α ) HCM .
(4) 
If  α 3 / 2 , then  τ 1 GGC , hence  H α HCM .
In addition to the above results, which are a continuation of Simon’s ones [14,15,16], we obtain the following ones that we divide, for sake of coherence, into Section 2 and Section 3:
  • Proposition 1 is a key to illustrate several other links between the r.v.s  T α (19), and the Mittag-Leffler functions;
  • A consequence of Proposition 1 is Corollary 3, which builds a new class of  HCM functions, similar to the shapes (8) and (20);
  • In Corollary 2 (and also in Theorem 3), we obtain that for any complex number z in the first quadrant, and for  0 < t 1 arg z , the real and imaginary part of  λ E t z λ t also enjoy  CM GGC and  HCM properties;
  • With the help of Proposition 1, we characterize in Theorem 3 the distributions of a peculiar family of distributions, which might be helpful in solving the open question (22);
  • We conclude with Corollary 4, which provides more information than
To the best of our knowledge, the results of Section 2 and Section 3 have not been addressed in the literature in a form even close to ours, and they constitute a contribution to infinite divisibility, especially in the Bondesson’s  HCM and Mittag-Leffler contexts.
In Section 4, we give comments and perspectives on some of the above results and also settle a prerequisite on the positive stable distributions on the Bernstein functions. The proofs are postponed to Section 5.

2. A New Class of HCM Distributions and  CM Property for  λ Re E t ( z λ t ) , z C

2.1. The Biasing and the Gamma-Mixture Procedure

Clearly,
E [ ( Z [ u ] ) λ ] = E [ Z λ + u ] E [ Z u ] and ( Z v ) [ u ] = d ( Z [ u v ] ) v , if E [ Z u v ] < .
A useful result is the stochastic interpretation of property (9):
Z HCM ( Z q ) [ u ] = d ( Z [ q u ] ) q HCM , for | q | 1   and   all   u   s.t. E [ Z u q ] < .
A random variable  B a , b , a , b > 0 has the beta distribution if it has the p.d.f.
Γ ( a + b ) Γ ( a ) Γ ( b ) x a 1 ( 1 x ) b 1 , 0 < x < 1 .
Observe that  G t [ u ] = d G t + u and  B a , b [ v ] = d B a + v , b , if  u > t and  v > a . By the beta-gamma algebra, we know that if  B a , b and  G a + b are independent, then
( B a , b , G a + b ) = d G u G a + G b , G a + G b ,
where G a and G b are independent .
The following fact will be used in the sequel to exhibit several factorizations in law. A positive r.v. X has p.d.f. of the form
f X ( x ) = x t 1 g ( x ) , t > 0 and g CM ,
if and only if, the distribution of X is a gamma-mixture of order t (shortly X is a  G t -mixture). The latter is equivalent to the existence of a positive r.v. Y such that we have the independent factorization
X = d G t Y ,
and in (42), necessarily
g ( x ) = 1 Γ ( t ) E Y t e x Y = E [ Y t ] Γ ( t ) E e x Y [ t ] , if E [ Y t ] < .
An important result due to Kristiansen [17] asserts that
X has a   G 2 mixture distribution X I D ,
and the beta-gamma algebra ensures that the same holds if the distribution of X is a  G t -mixture,  t < 2 .

2.2. A Generalization of Property (21)

In this work, the functions, defined for  c 0 and  t > 0 , by
ϕ c , t ( λ ) : = log 1 + 2 c λ t + λ 2 t , λ 0 ,
turn out to be crucial, and their properties will unblock several of our problems. Notice that the necessity of the condition  c 0 is essential since we need  ϕ c , t to be a Bernstein function in the sequel. Notice that  ϕ c , t is the composition  ϕ c , t ( λ ) = φ ( ψ c , t ( λ ) ) , where  φ ( λ ) : = log ( 1 + λ ) TBF , whereas
ψ c , t ( λ ) : = 2 c λ t + λ 2 t TBF
if , and only if , t 1 / 2 .
On the other hand, the class  CBF being stable by composition, we see that
λ λ 2 t 2 c λ t + λ 2 t
is a complete Bernstein function if t 1 / 2 .
By [3] (Theorem 7.3), we deduce that the logarithmic derivative
ψ c , t ( λ ) ψ c , t ( λ ) = t λ 1 + λ 2 t 2 c λ t + λ 2 t
is a Stieljes function if t 1 / 2 ,
i.e.,  ψ c , t has the representation by (2). By (78) below, we easily deduce that  ϕ c , t TBF , for any  c 0 and  t 1 / 2 . A larger range for t is provided by the following result.
Proposition 1.
Let  t > 0 and  c 0 . Then, the following holds.
(1) 
The function
ϕ c , t ( λ ) : = log ( 1 + 2 c λ t + λ 2 t ) , λ 0 ,
is Thorin–Bernstein if and only if  t 1 and  c + cos ( π t ) 0 .
(2) 
Let  α [ 0 , 1 ] and  c = cos ( π α ) . Then, the following assertions are equivalent.
(i) 
α 1 / 2 and  t 1 α ;
(ii) 
The function  ( 1 + 2 cos ( π α ) λ t + λ 2 t ) 1 is completely monotone;
(iii) 
The function  ϕ c , t is Bernstein;
(iv) 
The function  ϕ c , t is Thorin–Bernstein.
A consequence of Proposition 1 is the following result, which gives more information on the function  ϕ c , t ( λ ) , and provides a monotonicity property for the Mittag-Leffler function in case  c [ 0 , 1 ] .
Corollary 2.
Let  c [ 0 , 1 ] , j c : = c + i 1 c 2 . For  t ( 0 , 1 ) , let  E t be the Mittag-Leffler function and  S t be a positive stable random variable. Then, the following holds.
(1) 
The series
C c , t ( x ) : = n = 0 ( 1 ) n cos n arccos ( c ) Γ n t + 1 x n t , x 0 ,
is represented by
C c , t ( x ) = Re E t j c x t = E χ c x S t t , w h e r e χ c ( x ) : = e c x cos ( 1 c 2 x ) .
(2) 
The following assertions are equivalent.
(i) 
0 t 1 arccos ( c ) ;
(ii) 
the function  C c , t is completely monotone;
(iii) 
the function  1 C c , t is Bernstein;
(iv) 
the function  ϕ c , t ( λ ) = log ( 1 + 2 c λ t + λ 2 t ) is Thorin–Bernstein.
(3) 
Under any of the conditions in (2), the functions  ϕ c , t and  C c , t are represented by
ϕ c , t ( λ ) = 2 t 0 ( 1 e λ x ) C c , t ( x ) x d x , C c , t ( y ) = E [ e y E c , t ] , λ , y 0 ,
where  E c , t is a positive r.v. whose distribution gives no mass in 0. The Thorin measure  U c , t of  ϕ c , t , obtained by representation (4), is  U c , t ( d u ) = 2 t P ( E c , t d u ) .
(4) 
Let z be a complex number such that  ( z ) , I m ( z ) 0 . Then,  λ Re E t z λ t CM , if, and only if,  t 1 arg ( z ) .
Remark 1.
Using (23), (24), and (46), we obtain
C 1 , t ( x ) = E e ( x / S t ) t = E e x T t 1 / t = E t ( x t )
and
C 0 , t ( x ) = E cos x t / S t t = Re E t ( i x t ) ,
then the completely monotonicity of  C 1 , t is not a surprise. On the other hand, Proposition 1 gives
0 < t 1 ϕ 1 , t ( λ ) = 2 log ( 1 + λ t ) = 2 t 0 ( 1 e λ x ) E t ( x t ) x d x TBF E t ( x t ) CM , 0 < t 1 / 2 ϕ 0 , t ( λ ) = log ( 1 + λ 2 t ) = 2 t 0 ( 1 e λ x ) Re E t ( i x t ) x d x TBF Re E t ( i x t ) CM .
The trivial relation
Re E a ( i x ) = n 0 ( 1 ) n x 2 n Γ ( 2 n a + 1 ) = E 2 a ( x 2 ) , x 0 , a > 0 ,
complies with the equality  ϕ 1 , t ( λ ) = 2 ϕ 0 , t / 2 , t 1 .
We are now able to introduce a new class of  HCM p.d.f.s that are reminiscent of representation (8).
Corollary 3.
Probability density functions of the form
f c , a , t ( x ) = C x a 1 1 + 2 c x t + x 2 t , C , t , a , x > 0 ,
are  HCM , if and only if,  t 1 c 0 and  c + cos ( π t ) 0 . Pointwise limits of functions of the form  HCM :
C x a 1 N 1 + 2 c i x t i + x 2 t i γ i ,
where  N N , C , γ i > 0 , a R , 0 < t i 1 , c i 0 , a n d c i + cos ( π t i ) 0 . are  HCM .
Thanks to (20), the p.d.f. of  T α 1 / t = ( S α / S α ) α / t , t > 0 is explicitly given by
f T α 1 / t ( x ) = t x t 1 f T α ( x t ) = t sin ( π α ) π α x t 1 1 + 2 cos ( π α ) x t + x 2 t , x > 0 .
Therefore, the main result of Bosch [11] (Theorem 1.2) is a particular case of Corollary 3 with  c = cos ( π α ) and  a = t there. Indeed, the condition  0 < t 1 c 0 and  c + cos ( π t ) 0 is equivalent to
0 α 1 2 , 0 t 1 α ,
and then (21) holds true.

3. Stochastic Interpretation of the p.d.f.s in (48) and  CM Property for  λ Im E t ( z λ t ) , z C

Proposition 1 enables us to introduce the positive r.v.  X c , t GGC associated with the  TBF -function  ϕ c , t in (45), i.e.,
E [ e λ X c , t ] = 1 1 + 2 c λ t + λ 2 t = e ϕ c , t ( λ ) , 0 < t 1 , c 0 and c + cos ( π t ) 0 .
With the convention  S 1 = 1 , the random variable  X 1 , t (respectively for  X 0 , t ) is well defined for  t 1 (respectively,  0 < t 1 / 2 ). These r.v.s enjoy a simple independent factorization in law:
E [ e λ X 1 , t ] = 1 ( 1 + λ t ) 2 = E [ e λ S t G 2 1 / t ] X 1 , t = d S t G 2 1 / t ,
and
E [ e λ X 0 , t ] = 1 1 + λ 2 t = E [ e λ S 2 t G 2 1 / t ] X 0 , t = d S 2 t G 1 1 / t .
For  0 < t < s 1 and c as in Proposition 1, we have this algebra: if  S t / s is independent of  X c , s , then
E [ e λ X c , t ] = E [ e λ t / s X c , s ] = E [ e λ S t / s ( X c , s ) s / t ] X c , t = d S t / s ( X c , s ) s / t .
The latter is a reminiscent of the subordination relation for stable distributions: if  S α / γ , 0 < α < γ < 1 , is a stable r.v. independent of  S γ , then
E [ e λ S α ] = e λ α = e ( λ α / γ ) α = E [ e λ α / γ S α ] = E [ e λ S α / γ ( S γ ) γ / α ] ,
which yields
S α = d S α / γ ( S γ ) γ / α .
The ordinary generating function for the Chebychev’s polynomial of the second kind is
n 0 U n ( x ) u n = 1 1 2 x u + u 2 , | x | , | u | < 1 .
When  c [ 0 , 1 ] , we may take
c = cos ( π α ) , 0 α 1 2 , 0 t 1 α ,
and obtain the Laplace transform
E [ e λ X c , t ] = 1 1 + 2 cos ( π α ) λ t + λ 2 t = n 0 U n ( cos π α ) λ n t = n 0 ( 1 ) n sin π ( n + 1 ) α sin π α λ n t , | λ | < 1 .
Using the convention (56) and the form (49), we obtain
f T α 1 / t ( x ) = t sin ( π α ) π α x t 1 E [ e x X c , t ] , x > 0 .
Integrating the latter on  ( 0 , ) , we obtain
E [ ( X c , t ) t ] = π α sin ( π α ) Γ ( t + 1 ) .
By (34) and (42), we deduce
f T α 1 / t ( x ) = x t 1 Γ ( t ) E X c , t [ t ] t e x X c , t [ t ] , x > 0 .
Hence,
α ( 0 , 1 / 2 ] , t ( 0 , 1 α ] , and c = cos ( π α ) T α 1 / t = d G t ( X c , t ) [ t ] .
Then taking the Mellin transform in both sides of the latter, we obtain
E [ ( X c , t [ t ] ) λ ] = E [ ( X c , t ) ( λ + t ) ] E [ X c , t t ] = Γ ( t ) Γ ( 1 λ t ) Γ ( 1 + λ t ) Γ ( t + λ ) Γ ( 1 α λ t ) ) Γ ( 1 + α λ t ) , | λ | < t ,
or equivalently,
E [ X c , t x ] = π α sin ( π α ) Γ 2 x t Γ 1 + x t Γ ( 1 + x ) Γ 1 α ( 1 x t ) Γ 1 + α ( 1 x t ) = π α sin ( π α ) Γ ( t + 1 ) , if x = t Γ ( 1 x t ) Γ ( 1 + x t ) Γ ( 1 + x ) sin π α ( 1 x t ) sin ( π α ) , if x [ 0 , t ) ( t , 2 t ) .
Motivated by the duplication Formula (54), the link (57) with Chebychev’s polynomials, and (60), we explain the distribution of the r.v.s  X c , t , c ( 0 , 1 ) ; in the following result, the case  c = 1 , 0 , is described in (52) and (53).
Theorem 3. 
Let  α ( 0 , 1 / 2 ] , t ( 0 , 1 α ] c = cos ( π α ) , j c = e i π α and  j ¯ c = e i π α   .   T h e GGC  d i s t r i b u t i o n   o f   r . v .   X c , t   i s   d e s c r i b e d   b y   t h e   f o l l o w i n g .
(1) 
T h e   c . d . f .   a n d   t h e   p . d . f .   o f   t h e   X c , t   a r e
P ( X c , t x ) = 1 sin ( π α ) Im E t ( j c x t ) 1 j c , x 0 ,
f X c , t ( x ) = t sin ( π α ) x t 1 Im E t ( j c x t ) = t sin ( π α ) x t 1 Im E e j ¯ c ( x / S t ) t S t t , x > 0 .
(2) 
There exists a positive r.v.  D c , t , such that
E [ D c , t 2 t ] = 1 , D c , t [ t ] = d 1 D c , t [ t ] ,
whose Mellin transform is
E D c , t [ t ] λ = Γ ( t ) 2 Γ ( 1 λ t ) Γ ( 1 + λ t ) Γ ( t λ ) Γ ( t + λ ) Γ ( 1 α λ t ) Γ ( 1 + α λ t ) , | λ | < t ,
and such that we have the independent factorization in law
X c , t = d G 2 t D c , t GGC , a n d ( X c , t ) [ t ] = d G t D c , t [ t ] GGC ,
(recall the size biasing notation (34) for  Z [ u ] ). In particular, we have the Laplace transform representation:
E [ e λ D c , t [ 2 t ] ] = t Γ ( 2 t ) sin ( π α ) Im E t ( j c λ t ) λ t = t Γ ( 2 t ) sin ( π α ) λ t E e cos ( π α ) ( λ / S t ) t S t t sin sin ( π α ) ( λ / S t ) t , λ 0 .
(3) 
If z is a complex number such that  ( z ) , I m ( z ) 0 and  0 < t 1 arg ( z ) , then  λ Im E t ( z λ t ) / λ t CM .
The next result completes Theorem 3 and Bosch’s characterization (21).
Corollary 4.
For  α ( 0 , 1 ) and  t > 0 , we have the equivalences:
(1) 
α 1 / 2 and  t 1 α ;
(2) 
x ( 1 + 2 cos ( π α ) x t + x 2 t ) 1 CM ;
(3) 
x ( 1 + 2 cos ( π α ) x t + x 2 t ) 1 S 2 t ;
(4) 
the distribution of  T α 1 / t is a  G t -mixture;
(5) 
T α 1 / t HCM .
With  c = cos ( π α ) , we have the independent factorizations in law
T α 1 / t = S α S α α t = d G t ( X c , t ) [ t ] = d G t G t D c , t [ t ] HCM .

4. Comments and Prerequisite for the Proofs

4.1. Comments on Theorem 1

Our results in Theorem 1 merit some comments.
(a)
By (42) and by (28), it is clear that the function  F α is the p.d.f. of the independent quotient  G 1 / T | α 1 | 1 / α .
(b)
By (16), we have the convergence in law
T | α 1 | 1 / α d G 1 G 1 , as α 1 ,
where the independent ratio  G 1 / G 1 has the Pareto with the  HCM p.d.f.  ( x + 1 ) 2 , x > 0 . Note that the latter can be obtained using (21) and (66), with  t = 1 α 0 , and the stability of the  HCM class under pointwise limits. By (28), the latter also reads
lim α 1 F α ( λ ) = α E α ( λ α ) α | α = 1 = e λ + E α ( λ α ) α | α = 1 = 0 e λ x ( x + 1 ) 2 d x , λ 0 .

4.2. Comments on Theorem 2

The case  α ( 1 , 2 ) in Theorem 2 is specific.
(a)
The complete monotonicity of the function  G α in (67) is also a direct consequence of the one of  F α in (28). Indeed, since  F α 0 , then
F α CM F α CM G α = F α F α CM .
The functions
( α 1 ) F α ( λ ) = e λ α d d λ E α ( λ α ) = E [ e λ α τ 1 ] , α ( 1 , 2 ) ,
and  F α , α ( 0 , 1 ) ( 1 , 2 ) , were shown to be completely monotone in [16] (Equation (1.4)) and in [14] (Theorem 1.1), respectively. The connection with the r.v.  T 1 α , α ( 0 , 1 ) was not noticed there, nor was the  HCM property of the function  F α , α ( 0 , 1 ) ( 1 , 2 ) .
(b)
One has  D 1 / 2 ( x ) = ( π x ) 1 HCM . One could ask if there exists some  α ( 1 / 2 , 1 ) such that  G α HCM . Observe that, contrary to the case  α ( 1 , 2 ] , the function
x 1 + x x 2 α 2 cos ( π α ) x α + 1 , x 0 , α ( 1 / 2 , 1 ) ,
is completely monotone (one could even show that it is a Stieltjes function).
(c)
In [15] (Equation (3)) and [14] (Equation (2.1)), one can find the following representations valid for  α ( 0 , 1 ) ( 1 , 2 ] :
G α ( x ) = sin ( π α ) π 0 e x u u α 1 ( 1 + u ) u 2 α 2 cos ( π α ) u α + 1 d u , x > 0 ,
and
L G α ( λ ) = 0 e λ x G α ( x ) d x = λ α 1 1 λ α 1 , λ 0 .
In (32), we have seen that  G α ( 0 + ) = 1 , in case  α ( 1 , 2 ] , consequently, the function
sin ( π α ) π × u α 1 ( 1 + u ) u 2 α 2 cos ( π α ) u α + 1 , u > 0 ,
is a p.d.f. These expressions were not used in our approach.
(d)
Using (68), we see that Theorem 2 could be restated as follows: for  α ( 0 , 1 ) ( 1 , 2 ] , we have
λ α | α 1 | λ α 1 1 λ α 1 1 l λ 1 + 1 l λ = 1 HCM α 1 2 .
In [5] (Theorem 5.7.1), we have the following computation, valid for all  c ( 0 , 1 ) :
E exp λ G 1 × G 1 c G c = 1 λ c 1 λ if λ 1 and E exp G 1 × G 1 c G c = c .
Note that the latter function is  HCM , since it is the Laplace transform of the product and quotient of independent  HCM random variables. Using property (9), we see that
λ 1 c 1 λ θ c 1 λ θ 1 l λ 1 + 1 l λ = 1 HCM , θ [ 1 , 1 ] .
Taking  c = | α 1 | / α for  α [ 1 / 2 , 1 ) ( 1 , 2 ] , we retrieve the “if part” in (70), and also point (1) in Theorem 2. When  α > 2 , the function in (70) does not extend anymore to an analytic function on  C ( , 0 ] , and by [5] ( (property (ix) p. 68), it follows that this function cannot be  HCM .

4.3. Comments on Theorem 3

In (64), the l.h.s. term is certainly the Laplace transform of a distribution on the positive line, whereas the r.h.s. term is the imaginary part of the Laplace transform (with complex argument) of a signed function. We were not able to invert (64) to obtain the explicit distribution of  D c , t [ t ] with elementary computations. On the other hand, formula [18] ((FI II 812, BI (361)(9), p. 498) asserts that
0 x μ 1 e a x sin ( b x ) d x = Γ ( μ ) ( a 2 + b 2 ) μ 2 sin μ arctan a b , Re ( μ ) > 1 , Re ( a ) > | Im ( b ) | .
Using (64) and last formula, we may write: for all  λ ( t , t ) ,
E [ D c , t t + λ ] = 1 Γ ( t λ ) 0 u t λ 1 E [ D c , t 2 t e u D c , t ] d u = t Γ ( 2 t ) sin ( π α ) Γ ( t λ ) 0 u λ 1 E e cos ( π α ) ( u / S t ) t S t t sin sin ( π α ) ( u / S t ) t d u = Γ ( 2 t ) sin ( π α ) Γ ( t λ ) 0 E S t t λ v λ t 1 e cos ( π α ) ( v ) sin sin ( π α ) v d v = Γ ( 2 t ) sin ( π α ) Γ ( t λ ) Γ 2 + λ t Γ ( 1 + t + λ ) Γ λ t sin λ t π α = Γ ( 2 t ) sin ( π α ) Γ ( t λ ) Γ 1 + λ t Γ ( t + λ ) Γ 1 λ t sin λ t π α λ .
Letting  λ 0 , we obtain
E [ D c , t t ] = π α Γ ( 2 t ) t sin ( π α ) Γ ( t ) 2 .
Then, dividing (72) by the latter, we retrieve another computation of the complicated gamma ratio form of the Mellin transform (62). After searching in several books specialized in integral representations, such as [18], we were unsuccessful in finding this complicated expression. Finally, using identity (63), the fact that  D c , t [ t ] 1 / D c , t [ t ] , and the beta prime p.d.f. of  G t / G t explicitly given by  Γ ( 2 t ) Γ ( t ) 2 x t 1 ( 1 + x ) 2 t , x > 0 , we obtain the alternative representation
f T α 1 / t ( x ) = t sin ( π α ) π α x t 1 1 + 2 cos ( π α ) x t + x 2 t = Γ ( 2 t ) Γ ( t ) 2 x t 1 E D c , t [ t ] t x + D c , t [ t ] 2 t , x > 0 ,
where those parameters are as in Theorem 3. In particular,
x 1 1 + 2 cos ( π α ) x t + x 2 t S 2 t .

4.4. Comments on Corollary 4

Point (2) in Corollary (4) completes (21) and comforts conjecture (22). We aim to deepen the investigation of this conjecture in a forthcoming work.

4.5. Some Account of Stable Distributions

We will need the several identities for the positive stable distributions. Shanbhag and Sreehari [19] (Theorem 1) exhibited following independent factorization:
G r 1 / α = d G α r S α [ α r ] , r > 0 , 1 > α > 0 .
Next independent factorization more classical:
G 1 = d G 1 S α α ,
where in the l.h.s, the exponentially distributed r.v.  G 1 is independent of  S α , and can be easily seen from
P ( G 1 1 / α > λ ) = e λ α = E [ e λ S α ] = P ( G 1 > λ S α ) = P G 1 S α > λ , λ 0 .
Note that the last identity yields the Mellin transform computation in (15), and that the explicit expression of the p.d.f. of  T α is easily retrieved through the expression of the Mellin transform of  T α obtained from (15), the Euler’s reflection formula for the gamma function, and the residue theorem. Additionally, if  | s | < 1 , then
E T α s = Γ ( 1 s ) Γ ( 1 + s ) Γ ( 1 α s ) Γ ( 1 + α s ) = sin ( π α s ) sin ( π s ) = sin ( π α ) π 0 x s 1 + 2 cos ( π α ) x + x 2 d x ,
and the injectivity of the Mellin transform provides the p.d.f. of  T α in (20). Also note that a possible approach to conjecture (22) is Kanter’s factorization, which can be found in [20] (Corollary 4.1):
1 S α α / ( 1 α ) = d G 1 s α 1 / ( 1 α ) ( U ) ,
where  U denotes an uniform random variable on  ( 0 , 1 ) independent of  G 1 and  s α is the function defined on  ( 0 , 1 ) by
s α ( u ) = sin α π α u sin 1 α π ( 1 α ) u sin ( π u ) , u ( 0 , 1 ) .
Combining identities (73) and (76) with the beta-gamma algebra (40), we can complete (18) by the following: for  α ( 0 , 1 ) , it holds that
0 < t 1 α S α α t I D S α α t i s a G 1 m i x t u r e T α 1 t i s a G 1 m i x t u r e T α 1 t I D .
Thus, (18), (21), and (77) appear to give a certain credit to conjecture (22).

4.6. Some Account of Thorin and Complete Bernstein Functions

Note that function ϕ belongs to  TBF (respectively,  CBF ) if  ϕ BF if its Lévy measure Π in (3) has a density function of the form  L U ( x ) / x (respectively,  L V ( x ) ),  x > 0 , where the U and V and are positive measures on  ( 0 , ) satisfying (5). Both classes  CBF and  TBF are convex cones that are closed under pointwise limits; the class  CBF is stable by composition, and by [3] (Theorem 8.4), for a Thorin–Bernstein function ϕ, we have
φ ϕ TBF , φ TBF ϕ ϕ S 1 .
It is immediate that
ϕ CBF λ ϕ 1 λ S 1 λ ϕ ( λ ) λ S 1 .
The following fact is much less evident. By [3] (Theorem 7.3), we also have
ϕ CBF 1 ϕ S 1 .
An important representation for the logarithmic derivative of  ϕ CBF is provided by [3] (Theorem 6.17): there exists (a unique) pair  γ R and a measurable function  η : ( 0 , ) [ 0 , 1 ] such that
ϕ ( λ ) ϕ ( λ ) = ( 0 , ) η ( u ) ( λ + u ) 2 d u , λ > 0 .
For instance, the Thorin–Bernstein function  ϕ α ( λ ) = λ α , α ( 0 , 1 ) , has an η-function equal to  α 1 l u > 0 .

5. The Proofs

Proof of Theorem 1.
(1) From [14] (Equation (3.4) and (3.6)), we have
e λ α E α ( λ α ) = α sin ( π α ) π 0 e λ t t α 1 t 2 α 2 cos ( π α ) t α + 1 d x , λ 0 .
After some identification in (49), we obtain the Laplace representation (28):
e λ α E α ( λ α ) = ( 1 α ) E e λ T 1 α 1 / α , if α ( 0 , 1 ) , E e λ T α 1 1 / α , if α ( 1 , 2 ] .
(2a), we have
1 α e λ E α ( λ α ) = ( 1 α ) E [ e λ R α ] , e λ E α ( λ α ) 1 = 1 α α 1 E [ e λ R α ] .
Let  R n ( α ) be the random walk generated by  R α given by (25), i.e.,  R 0 ( α ) = 0 , R n ( α ) = R 1 , α + + R n , α , n 1 , where  R 1 , α , , R n , α , are n independent copies of  R α . Let  Z α be geometrically distributed on  N with parameter α, and independent of the sequence  ( R n ( α ) ) n 0 . Using (25), (28), and the subordinated r.v.  R Z α ( α ) , we have the Laplace representation
e λ E α ( λ α ) = α 1 ( 1 α ) E [ e λ R α ] = n = 0 α ( 1 α ) n E [ e λ R α ] n = n = 0 P ( Z α = n ) E [ e λ R n ( α ) ] = E [ e λ R Z α ( α ) ] , λ 0 .
(2b) From Kanter’s factorization (76), we see that  S 1 α ( 1 α ) / α has a completely monotone density, and so does  T 1 α 1 / α . We deduce that  F α is a Stieltjes function. The equivalence between  F α HCM , and  1 / 2 α < 1 , is due to (21). The function  F α / F α is Stieltjes because  log F α TBF , and because of (6). The function  1 α e λ E α ( λ α ) is  HCM due to (26) and (82).
(3a) We perform as in (2a).
(3b) By (21), we know that if  t > 2 ( 2 α ) , then the distribution of  T α 1 1 / t is not  GGC . Since  α > 2 ( 2 α ) in case  α > 4 / 3 , then  F α could not be  HCM .
(3c) For the last assertion, use (11), (17) and (21), and observe that  0 < γ ( 2 α ) / α < 1 . Then,
1 < α 3 2 0 < α 1 < 1 2 S γ T α 1 1 / γ α GGC F α ( λ γ ) = E e λ γ T α 1 1 / α = E e λ S γ T α 1 1 / γ α HCM .
As in (2), the mixture property (77) entails  T α 1 1 / γ α has a completely monotone density; hence, so does  S γ T α 1 1 / γ α . The latter proves that  F α ( λ γ ) is Stieltjes; for its logarithmic derivative, we conclude as in the end of the proof of point (2b).
Proof of Theorem 2.
(1) Assume  α ( 1 , 2 ] . The function  L D 2 ( λ ) = ( λ + 1 ) 1 is clearly an  HCM function, and by (10), we deduce  D 2 GGC . If  α ( 1 , 2 ) , then the function
1 + x x 2 α 2 cos ( π α ) x α + 1 , x > 0 ,
is locally increasing in a neighborhood of  0 + and cannot be completely monotone. Suppose that  G α HCM . By [5] (Theorem 5.4.1), the function  G α is then the Laplace transform of a  GGC and by [5] (Theorems 4.1.1 and 4.1.4), the Thorin mass in (69) equals α; thus, the function in (83) would be completely monotone, a contradiction. We deduce that  G α is not  HCM when  α ( 1 , 2 ] . To show that  G α GGC , it suffices to show that  G 1 , α : = log L G α is a Thorin–Bernstein function (equivalently,  G 1 , α is a Stieltjes function), using formula (68). The function  G 1 , α (and then  G 1 , α ) extends to an analytic function on  C R . We will use the characterization of the Stieltjes function given [3] (by Corollary 7.4), namely, we aim to prove that
Im ( z ) > 0 Im ( G ( z ) ) < 0 .
Observe that  Im ( G 1 , α ) is harmonic on the upper half-plane as the imaginary part of an analytic function. Moreover,  G 1 , α ( z ) 0 uniformly as  | z | + . Then, from a compactness argument and from the minimum principle, it suffices to show that
lim sup z 0 , Im ( z ) > 0 Im ( G 1 , α ( z ) ) 0 .
Elementary computations give
G 1 , α ( z ) = α z α 1 z α 1 ( α 1 ) z α 2 z α 1 1 , z C R ,
and for all  x > 0 , we have
G 1 , α ( x + ) : = lim z x Im ( z ) > 0 G 1 , α ( z ) = e i π α x α 2 α 1 x α 1 e i π α + 1 α x x α e i π α 1 = x α 2 x α e i 2 π α + ( α 1 ) e i π α + α x e i π α x 2 α 1 e i 2 π α + ( x α x α 1 ) e i π α 1 , Im ( G 1 , α ( x + ) ) = sin ( π α ) x α 2 A 1 , α ( x ) B 1 , α ( x ) ,
where
A 1 , α ( x ) = ( α 1 ) x 2 α + α x 2 α 1 + 2 cos ( π α ) x α + α x + α 1 ,
and
B 1 , α ( x ) = x 2 α 1 e i 2 π α + ( x α x α 1 ) e i π α 1 2 .
The function  A 1 , α is a positive on  [ 0 , ) . Indeed,
A 1 , α ( x ) = α 2 ( α 1 ) x 2 α 1 + ( 2 α 1 ) x 2 α 2 + 2 cos ( π α ) x α 1 + 1 α x 2 α 2 2 x α 1 + 1 = α x α 1 1 2 0 ,
thus,  A 1 , α is a non-decreasing function on  [ 0 , ) , and since  A 1 , α ( 0 ) = α 1 > 0 , we deduce that  A 1 , α is positive and then  Im ( G 1 , α ( x + ) ) 0 for all  x > 0 . Next, since  1 < α 2 < 0 , then
lim sup z 0 Im ( z ) > 0 Im ( G 1 , α ( z ) ) = lim sup z 0 Im ( z ) > 0 ( α 1 ) z α 2 = ( α 1 ) lim sup r 0 , r > 0 θ ( 0 , 1 ) r α 2 sin ( π ( α 2 ) θ ) = 0 ,
Finally,  Im ( G 1 , α ( x ) ) = 0 , for all  x > 0 . All in all, we have proved that the function  G 1 , α is Stieltjes, hence  G 1 , α is Thorin–Bernstein.
(2) Assume  α ( 0 , 1 ) . Using [5] (Theorem 5.4.1) again and the definition of a widened  GGC in [5] (Section 3.5), we obtain
G α is   a   widened GGC λ L ( G α ) ( λ ) = 1 λ α 1 λ α 1 HCM λ λ 1 α 1 λ α 1 HCM ,
and it suffices to prove that the function
G 2 , α ( λ ) : = log λ 1 α 1 λ α 1 , λ 0 ,
is Thorin–Bernstein (equivalently,  G 2 , α is a Stieltjes function) if and only if  α 1 / 2 . Elementary computations give
G 2 , α ( z ) = α z α 1 z α 1 ( 1 α ) z α z 1 α 1 , z C ( , 0 ] ,
and then,
Im ( G 2 , α ( x + ) ) = sin ( π α ) A 2 , α ( x ) B 2 , α ( x ) , x > 0 ,
where
A 2 , α ( x ) = ( 1 α ) ( x α + x α ) α ( x 1 α + x α 1 ) 2 cos ( π α )
and
B 2 , α ( x ) = x 1 + x α e i π α x 1 α e i π α 2 .
We proceed as we did for  G 1 , α to show that  G 2 , α is a Stieltjes function. For this, we have to check the sign of  A 2 , α . The function  H 1 / 2 = 0 is a trivial Thorin–Bernstein function. Since for  α 1 / 2
A 2 , α ( 1 ) = 2 ( 1 2 α cos ( π α ) ) < 0 , if α ( 0 , 1 / 2 ) > 0 , if α ( 1 / 2 , 1 ) ,
we see that  α ( 1 / 2 , 1 ) is a necessary condition for  A 2 , α to be positive. Next, let
C α ( u ) = A 2 , α ( e u ) = 2 ( 1 α ) cosh ( α u ) α cosh ( α u ) cos ( π α ) , u R .
With the expression
C α ( u ) = 2 α ( 1 α ) sinh ( α u ) sinh ( ( 1 α ) u ) ,
one deduces that if  α ( 1 / 2 , 1 ) , then the function  A 2 , α decreases on  ( 0 , 1 ) and increases on  ( 1 , + ) ; thus,  A 2 , α ( x ) A 2 , α ( 1 ) > 0 for all  x > 0 . Finally, we deduce that if  α ( 1 / 2 , 1 ) , then
lim sup z 0 Im ( z ) > 0 Im ( G 2 , α ( z ) ) = lim sup r 0 , r > 0 θ ( 0 , 1 ) sin ( π α θ ) ( α r α 1 ( 1 α ) e α ) = 0 ,
and this shows that  G 2 , α is a Stieltjes function.
The  GGC property of  G 1 / R α is straightforward thanks to (31). □
Proof of Corollary 1.
If  α = 2 , then  E 2 ( λ ) = cosh ( λ ) ( X t ( 2 ) ) t 0 is Brownian motion, and  τ 1 = d S 1 / 2 . The rest of the assertions are straightforward due to the convention  T 0 = 1 . We then treat the case  α ( 1 , 2 ) .
(1) Using [16] (Equation (1.4)) and the definition of  F α in (28), one has the Laplace representation
E [ e λ α τ 1 ] = ( 1 α ) F α ( λ ) = H α ( λ α ) ,
and this gives the Laplace transform representation (36).
(2) From (28) and (85), we obtain
E [ e λ τ 1 ] = ( α 1 ) E T α 1 1 / α e λ 1 / α T α 1 1 / α = ( α 1 ) E T α 1 1 / α e λ S 1 / α T α 1 ,
where  S 1 / α is a positive r.v. independent of  T α . Further, (86) or (32) imply
1 = ( α 1 ) E [ T α 1 1 / α ] ,
hence,
E [ e λ τ 1 ] = E [ T α 1 1 / α e λ S 1 / α T α 1 ] E [ T α 1 1 / α ] = E e λ S 1 / α ( T α 1 ) [ 1 / α ] ,
which shows the first factorization in (37). The second factorization in (37) is due to (74) and to the second identity in (38) applied to  T α 1 .
(3) The second factorization in (37) gives
G 1 τ 1 1 / ( 2 α ) = d G 1 α / ( 2 α ) T α 1 1 / ( 2 α ) [ ( 2 α ) / α ] .
Since  q = α / ( 2 α ) > 1 , then the power property in (11) ensures that  G 1 α / ( 2 α ) HCM . Property (21) with  t = ( 2 α ) asserts that  T α 1 1 / ( 2 α ) HCM , and by property (39) we deduce that  T α 1 1 / ( 2 α ) [ ( 2 α ) / α ] HCM . Finally, the stability property by independent products in (11) yields  ( G 1 / τ 1 ) 1 / ( 2 α ) HCM .
(4) We know that  S 1 / α GGC . Then,
0 < α 1 1 2 ( 21 ) T α 1 HCM ( 39 ) ( T α 1 ) [ 1 / α ] HCM ( 11 ) τ 1 = d S 1 / α ( T α 1 ) [ 1 / α ] GGC ,
the  HCM property of  H α follows from (10). □
Proof of Proposition 1.
(1) By (78)–(80), we have
ϕ c , t TBF λ ϕ c , t ( λ ) = 2 t λ c λ t + λ 2 t 1 + 2 c λ t + λ 2 t S 1 λ c λ t + λ 2 t 1 + 2 c λ t + λ 2 t CBF λ 1 + 2 c λ t + λ 2 t c λ t + λ 2 t = 1 + 1 + c λ t c λ t + λ 2 t S 1 λ 1 + c λ t c λ t + λ 2 t S 1 λ φ c , t ( λ ) : = c λ t + λ 2 t 1 + c λ t CBF .
(1)
The case  t = 1 . We have  φ 1 , 1 ( λ ) = 1 and the function  φ c , 1 ( λ ) / λ = ( c + λ ) / ( 1 + c λ ) is not completely monotone if  c 1 , because its derivative  c ( 1 c ) ( 1 λ ) / ( 1 + c λ ) 2 has a change sign. Thus,  φ c , 1 BF if  c 1 .
(2)
The case  c = 1 . Trivially,  φ 1 , t ( λ ) = λ t CBF t 1 .
(3)
The case  c 1 . The function  φ c , t is  CBF if and only if  φ c , t satisfies (80). We then study the logarithmic derivative of  φ c , t :
φ c , t φ c , t ( λ ) = t λ c λ t + 2 λ 2 t c λ t + λ 2 t c λ t 1 + c λ t = t λ λ 2 t c λ t + λ 2 t + 1 c λ t 1 + c λ t = t λ λ t c + λ t + 1 1 + c λ t .
(3a) The case  c 1 , t > 1 . Here, we have a first conclusion: if  t > 1 , then  φ c , t does not belong to  CBF , and not even to  BF , since
φ c , t ( λ ) λ λ t 1 + , when λ + .
(3b) The case  c 1 , 0 < t < 1 . Recall the r.v.  T t given in (20). Observe that if  G 1 is standard exponentially distributed and is independent of  T t , then
1 1 + λ t = E [ e λ G 1 T t 1 / t ] , λ t 1 + λ t = 1 E [ e λ G 1 T t 1 / t ] .
Using the fact that used that  T t 1 / t = d T t 1 / t and the latter, we obtain
1 λ 1 1 + c λ t = 1 λ 1 1 + ( c 1 / t λ ) t = 1 λ E e λ G 1 ( c T t ) 1 / t = 0 e λ x P ( G 1 ( c T t ) 1 / t x ) d x = 0 e λ x x 1 E [ e x ( T t / c ) 1 / t ] x d x = 0 e λ x x L f ( x ) d x = 0 f ( u ) ( λ + u ) 2 d u ,
where
f ( u ) = P ( T t 1 / t > u c 1 / t ) = P ( T t > c u t ) , u 0 .
Similarly, we have
1 λ λ t λ t + c = 1 λ ( λ / c 1 / t ) t 1 + ( λ / c 1 / t ) t = 1 λ 1 E e λ G 1 ( T t / c ) 1 / t = 0 e λ x P ( G 1 ( T t / c ) 1 / t > x ) d x = 0 e λ x x E e x ( c T t ) 1 / t x d x = 0 e λ x x L g ( x ) d x = 0 g ( u ) ( λ + u ) 2 d u ,
where
g ( u ) = P c T t ) 1 / t < u = P ( T t < u t / c ) , u 0 .
We finally obtain the representation of the logarithmic derivative of  φ c , t :
φ c , t φ c , t ( λ ) = 0 η t ( u ) ( λ + u ) 2 d u ,
where
η t ( u ) : = t f ( u ) + g ( u ) = t P ( T t > c u t ) + P ( T t < u t / c ) .
Due to (81), we can now assert that
ϕ c , t TBF η c , t ( u ) [ 0 , 1 ] , u > 0 .
We arrive to the last step. Since  T t has the explicit density
f T t ( u ) = sin ( π t ) π t 1 + 2 cos ( π t ) u + u 2 , u > 0 ,
then the derivative of  η ˜ c , t ( u ) : = η c , t ( u t ) is
η ˜ c , t ( u ) = t 1 c f T t u c c f T t ( c u ) = sin ( π t ) c π 1 c 2 + 2 c cos ( π t ) u + u 2 1 1 + 2 c cos ( π t ) u + c 2 u 2 ,
and finally,
η ˜ c , t ( u ) 0 1 + 2 c cos ( π t ) u + c 2 u 2 c 2 + 2 c cos ( π t ) u + u 2 1 c 2 ( 1 c 2 ) u 2 u [ 1 , ) .
(3b)(i) The case  c > 1 , 0 < t < 1 . By (88), the function  η c , t decreases on  [ 0 , 1 ] , then increases on  [ 1 , ) , and its maximum is equal to  t = η c , t ( 0 ) = η c , t ( ) . We deduce that (87) is true.
(3b)(ii) The case  c < 1 , 0 < t < 1 . By (88), the function  η c , t increases on  [ 0 , 1 ] , then decreases on  [ 1 , ) , and its maximum is equal to
η c , t ( 1 ) = t P ( T t > c ) + P ( T t < 1 / c ) = 2 t P ( T t > c ) = 2 t c sin ( π t ) π t ( 1 + 2 cos ( π t ) u + u 2 ) = 2 π π 2 arctan c + cos ( π t ) sin ( π t ) .
Note that we have performed the obvious change of variable  u = sin ( π t ) v cos ( π t ) to obtain the last expression. We deduce that (87) is true if and only if
arctan c + cos ( π t ) sin ( π t ) 0 c + cos ( π t ) 0 .
(2) By (1), we have  ( i ) ( i v ) . Trivially,  ( i v ) ( i i i ) ( i i ) . It suffices to check  ( i i ) ( 1 ) . Assume
g α , t ( λ ) : = 1 1 + 2 cos ( π α ) λ t + λ 2 t CM .
If  α ( 1 / 2 , 1 ) , then the denominator has a sign change, and  g α , t could not be  CM . Assume  α ( 0 , 1 / 2 ] and  t > 1 α , then  g α , t has two poles  e ± i π ( 1 α ) / t . In this case,  g α , t cannot extend to an analytic function on  { z C / Re ( z ) > 0 } , then, it could not be  CM . □
Proof of Corollary 2.
(1) Using (24), note that
χ c ( u ) = Re e j c u a n d C c , t ( x ) = Re E t ( j c x t ) = Re E χ c ( x / S t ) t .
(2) The equivalences are obtained by the form of the derivative of  ϕ c , t :
ϕ c , t ( λ ) = 2 t c λ t 1 + λ 2 t 1 1 + 2 c λ t + λ 2 t = 2 t λ 1 s c , t ( λ ) , s c , t ( λ ) = c λ t + 1 1 + 2 c λ t + λ 2 t λ > 0 .
By the latter, Proposition 1, and by (6) and (79), we have the equivalences
c + cos ( π t ) 0 ϕ c , t TBF ϕ c , t S 1 1 s c , t CBF s c , t S 1 .
Using the representations
χ c ( u ) = e c u c cos ( 1 c 2 u ) + 1 c 2 sin 1 c 2 u ,
x 1 + x 2 = 0 e x u cos ( u ) d u , 1 1 + x 2 = 0 e x u sin ( u ) d u , x > 0 ,
and (13), we may write
s c , t ( λ ) = c λ t + 1 1 + 2 c λ t + λ 2 t = c 1 c 2 λ t + c 1 c 2 1 + λ t + c 1 c 2 2 + 1 1 + λ t + c 1 c 2 2 = 0 e ( λ t + c ) u c cos ( 1 c 2 u ) + 1 c 2 sin 1 c 2 u d u = 0 e λ t u χ c ( u ) d u = 0 E e λ u 1 / t S t χ c ( u ) d u .
Applying Fubini–Tonelli’s theorem, and performing the change of variable  v = u 1 / t S t under the expectation, we obtain the Laplace transform representation of  s c , t :
s c , t ( λ ) = 0 e λ v t v t 1 E 1 S t t χ c v S t t d v = 0 e λ v d d v E χ c v S t t d v = 0 e λ v C c , t ( v ) d v .
Thus,  s c , t S 1 , if, and only if,  C c , t CM and  C c , t (hence decreases from  C c , t ( 0 ) = 1 to  C c , t ( ) = 0 , thus  C c , t [ 0 , 1 ] ). The latter is equivalent to  C c , t CM or to  1 C c , t BF .
(3) Using (89), (91), and then performing an integration by parts, we obtain the expression
ϕ c , t ( λ ) = 2 t λ 0 ( 1 e λ v ) C c , t ( v ) d v = 2 t 0 e λ x C c , t ( x ) d x .
Integrating the last expression from 0 to λ, we retrieve (47). The representation of  C c , t as the Laplace transform or a positive r.v.  E c , t is evident since  C c , t is completely monotone and  C c , t ( 0 ) = 1 . The last assertion is also evident due to the alternative Frullani integral form of (47).
ϕ c , t ( λ ) = 2 t E log 1 + λ E c , t = 2 t 0 log 1 + λ u P ( E c , t d u ) .
(4) It suffices to write  z = | z | j c , where  c = cos ( arg z ) . □
Proof of Corollary 3.
The first assertion is a straightforward consequence of the Thorin property of the Bernstein function  ϕ c , t in Proposition 1. Indeed, the definition of  GGC distributions gives that
ϕ c , t TBF x e ϕ c , t ( x ) = 1 1 + 2 c x t + x 2 t HCM .
For the second assertion, use the fact that the  HCM class is stable by product, closed by pointwise limits [5] (property (ii) p. 68), and property (9). □
Proof of Theorem 3.
(1) We proceed as in the proof of Corollary 2. Let us define
σ c ( x ) : = Im e j c x 1 j c a n d S c , t ( x ) : = 1 1 c 2 E σ c x S t t , x 0 .
and observe that
σ c ( x ) = e c x sin 1 c 2 x = Im ( e j c ¯ x ) .
Using representations (51) and (90), then performing the change of variable  u = ( x / S t ) t under the expectation, we obtain
E [ e λ X c , t ] = 1 1 + 2 c λ t + λ 2 t = 1 1 c 2 1 1 + λ t + c 1 c 2 2 = 1 1 c 2 0 e ( λ t + c ) u sin 1 c 2 u d u = 1 1 c 2 0 E [ e λ u 1 / t S t ] σ c ( u ) d u = t 1 c 2 0 e λ x x t 1 E S t t σ c ( x / S t ) t d x
The latter gives the expression of  f X c , t ( x ) , and by integration over  ( 0 , x ] , we obtain the one of  S c , t ( x ) = P ( X c , t x ) .
(2) By (51), and Corollary 3, recall that the distribution  X c , t is  GGC , and is associated with the Thorin–Bernstein function  ϕ c , t in (47) and to the Thorin measure  U c , t given by (4). Using Corollary 2, then [5] (Theorem 4.1.1), or equivalently [5] (Theorem 4.1.4), we see that the total mass of  U c , t equals
U c , t ( 0 , ) = 2 t = sup s ; lim x 0 + x 1 s f X c , t ( x ) = 0 .
and we conclude that  X c , t GGC if and only if identity (63) holds. The Laplace transform in (64) is identified by the equality  E [ D c , t 2 t ] = 1 , which is obtained by taking  lim x 0 + x 1 t f X c , t ( x ) in (61). By (42), (61) and (63), it is straightforward that
E [ D c , t t ] E D c , t [ t ] t e λ D c , t [ t ] = E [ D c , t 2 t e λ D c , t ] = Γ ( 2 t ) λ 1 2 t f X c , t ( λ ) = t Γ ( 2 t ) λ t 1 c 2 Im E e j c ¯ ( λ / S t ) t S t t = t Γ ( 2 t ) λ t sin ( π α ) E e cos ( π α ) λ / S t t S t t sin sin ( π α ) λ / S t t , λ 0 .
(3) This performed as in proof of point (3) of Corollary 2. □
Proof of Corollary 4.
Equivalences  ( 1 ) ( 2 ) ( 3 ) are due to point (2) of Proposition 1 and to (72),  ( 2 ) ( 4 ) corresponds to (41) ⟺ (42), and  ( 1 ) ( 5 ) is due to point (2) of Proposition 1 and to (51). For the last assertion, assume  α ( 0 , 1 / 2 ) 0 < t 1 α c = cos ( π α ) in (51), then consider the  GGC r.v.  X c , t , which is linked to  T α 1 / t by the expression
f T α 1 / t ( x ) = t sin ( π α ) π α x t 1 1 + 2 cos ( π α ) x t + x 2 t = t sin ( π α ) π α x t 1 E [ e x X c , t ] .
Integrating both sides in (92), we see that
E [ ( X c , t ) t ] = π α sin ( π α ) Γ ( t + 1 ) < ,
which by the procedure (34) allows one to introduce the r.v.  X c , t [ t ] , whose distribution is also a  GGC , by (9) and (10), we obtain the equivalences
X c , t GGC E [ e x X c , t ] HCM f T α 1 / t ( x ) HCM .
Using (63) and the effect of size-biasing in (38), we obtain
X c , t [ t ] = d G 2 t D c , t [ t ] = d G t D c , t [ t ] ,
hence,  X c , t [ t ] is also a  G t -mixture. Then, (42) and (93) yield
f T α 1 / t ( x ) = t sin ( π α ) π α x t 1 E [ e x X c , t ] = t sin ( π α ) π α x t 1 E [ ( X c , t ) t ( X c , t ) t e x X c , t ] = t sin ( π α ) π α E [ ( X c , t ) t ] x t 1 E [ ( X c , t [ t ] ) t e x X c , t [ t ] ] = x t 1 Γ ( t ) E [ ( X c , t [ t ] ) t e x X c , t [ t ] ] = f G t / X c , t [ t ] ( x ) ,
which, combined with (94), gives (65). □

6. Conclusions and Perspectives

In this paper, we showed that the Mittag-Leffler function (with an eventually complex argument) is tightly linked to the stable distributions by various aspects: we exhibited its non-trivial infinite divisibility,  GGC and  HCM properties, and we provided its explicit intervention in the distributional properties for the first passage time of the spectrally positive stable process. We also introduced new classes of  HCM distributions, and gave a possible direction to solve the open question (22) on the power of the positive stable r.v.s. Indeed, (65) gives
S α S α α t = d G t ( X c , t ) [ t ] = d G t G t D c , t [ t ] HCM , i f α ( 0 , 1 / 2 ] , c = cos ( π α ) a n d t = 1 α .
The latter indicates that two independent factorizations are feasible, namely,
S α α t = d G t Z c , t , D c , t [ t ] = d Z c , t Z c , t , Z c , t , Z c , t i . i . d .
In other terms, the distribution of  log D c , t [ t ] is symmetric. More investigation of the distribution of  D c , t [ t ] is then necessary to solve the open question (22).

Author Contributions

Writing—review & editing, N.A. and W.J. All authors have read and agreed to the published version of the manuscript.

Funding

Researchers Supporting Project (RSP2023R162), King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

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

Acknowledgments

We thank the two anonymous reviewers for their valuable comments and suggestions which substantially improved the presentation and the results of this article. The work of the second author was supported by the “Researchers Supporting Project (RSP2023R162), King Saud University, Riyadh, Saudi Arabia”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag–Leffler Functions (Related Topics and Applications); Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  2. Hirschman, I.I.; Widder, D.V. The Convolution Transform; Princeton University Press: Princeton, NJ, USA, 1955; Volume 20. [Google Scholar]
  3. Schilling, R.L.; Song, R.; Vondraček, Z. Bernstein Functions. Theory and Applications, 2nd ed.; de Gruyter Studies in Mathematics; Walter de Gruyter & Co.: Berlin, Germany, 2012; Volume 37. [Google Scholar]
  4. Steutel, L.; Van Harn, K. Infinite Divisibility of Probability Distributions on the Real Line; Marcel Dekker: New York, NY, USA; Basel, Switzerland, 2004. [Google Scholar]
  5. Bondesson, L. Generalized Gamma Convolutions and Related Classes of Distributions and Densities; Lecture Notes in Statistics; Springer: New York, NY, USA, 1992; Volume 76. [Google Scholar]
  6. Bondesson, L. A class of probability distributions that is closed with respect to addition as well as multiplication of independent random variables. J. Theor. Probab. 2015, 28, 1063–1081. [Google Scholar] [CrossRef]
  7. Nolan, J. Bibliography on Stable Distributions, Processes and Related Topics. Available online: https://edspace.american.edu/jpnolan/wp-content/uploads/sites/1720/2022/05/StableBibliography.pdf (accessed on 25 September 2023).
  8. Zolotarev, V.M. One Dimensional Stable Laws; Translations of Mathematical Monographs; American Mathematical Society: Providence, RI, USA, 1986. [Google Scholar]
  9. Bosch, P.; Simon, T. A proof of Bondesson’s conjecture on stable densities. Ark. Mat. 2016, 54, 31–38. [Google Scholar] [CrossRef]
  10. Jedidi, W.; Simon, T. Further examples of GGC and HCM densities. Bernoulli 2013, 19, 1818–1838. [Google Scholar] [CrossRef]
  11. Bosch, P. HCM property and the half-cauchy distribution. Probab. Math. Stat. 2015, 35, 191–200. [Google Scholar]
  12. James, L.F.; Roynette, B.; Yor, M. Generalized Gamma Convolutions, Dirichlet means, Thorin measures, with explicit examples. Probab. Surv. 2008, 5, 346–415. [Google Scholar] [CrossRef]
  13. Bridaa, S.; Jedidi, W.; Sendov, H. Generalized unimodality and subordinators, with applications to stable laws and to the Mittag-Leffler function. J. Theor. Probab. 2023. [Google Scholar] [CrossRef]
  14. Simon, T. Mittag-Leffler functions and complete monotonicity. Integral Transforms Spec. Funct. 2015, 26, 36–50. [Google Scholar] [CrossRef]
  15. Simon, T. Fonctions de Mittag-Leffler et processus de Lévy stables sans sauts négatifs. EXPO Math. 2010, 28, 290–298. [Google Scholar] [CrossRef]
  16. Simon, T. Hitting densities for spectrally positive stable processes. Stochastics 2011, 83, 203–214. [Google Scholar] [CrossRef]
  17. Kristiansen, G.K. A proof of Steutel’s conjecture. Ann. Probab. 1994, 22, 442–452. [Google Scholar] [CrossRef]
  18. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 7th ed.; Elsevier: Amsterdam, The Netherlands, 2007. [Google Scholar]
  19. Shanbhag, D.N.; Sreehari, M. An extension of Goldie’s result and further results in infinite divisibility. Z. Wahrscheinlichkeitstheorie Verwandte Geb. 1979, 47, 19–25. [Google Scholar] [CrossRef]
  20. Kanter, M. Stable densities under change of scale and total variation inequalities. Ann. Probab. 1975, 3, 697–707. [Google Scholar] [CrossRef]
Figure 1. Classes of infinitely divisible distributions, and motivations.
Figure 1. Classes of infinitely divisible distributions, and motivations.
Mathematics 11 04141 g001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Altaymani, N.; Jedidi, W. New Monotonicity and Infinite Divisibility Properties for the Mittag-Leffler Function and for Stable Distributions. Mathematics 2023, 11, 4141. https://0-doi-org.brum.beds.ac.uk/10.3390/math11194141

AMA Style

Altaymani N, Jedidi W. New Monotonicity and Infinite Divisibility Properties for the Mittag-Leffler Function and for Stable Distributions. Mathematics. 2023; 11(19):4141. https://0-doi-org.brum.beds.ac.uk/10.3390/math11194141

Chicago/Turabian Style

Altaymani, Nuha, and Wissem Jedidi. 2023. "New Monotonicity and Infinite Divisibility Properties for the Mittag-Leffler Function and for Stable Distributions" Mathematics 11, no. 19: 4141. https://0-doi-org.brum.beds.ac.uk/10.3390/math11194141

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

Article Metrics

Back to TopTop