Next Article in Journal
Information Dynamics of a Nonlinear Stochastic Nanopore System
Next Article in Special Issue
On f-Divergences: Integral Representations, Local Behavior, and Inequalities
Previous Article in Journal
Robust Covariance Estimators Based on Information Divergences and Riemannian Manifold
Previous Article in Special Issue
Some Inequalities Combining Rough and Random Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Logarithmic Sobolev Inequality and Exponential Convergence of a Markovian Semigroup in the Zygmund Space

Department of Mathematics, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
Entropy 2018, 20(4), 220; https://doi.org/10.3390/e20040220
Submission received: 29 December 2017 / Revised: 19 March 2018 / Accepted: 19 March 2018 / Published: 23 March 2018
(This article belongs to the Special Issue Entropy and Information Inequalities)

Abstract

:
We investigate the exponential convergence of a Markovian semigroup in the Zygmund space under the assumption of logarithmic Sobolev inequality. We show that the convergence rate is greater than the logarithmic Sobolev constant. To do this, we use the notion of entropy. We also give an example of a Laguerre operator. We determine the spectrum in the Orlicz space and discuss the relation between the logarithmic Sobolev constant and the spectral gap.

1. Introduction

Let ( M , B , m ) be a measure space with m ( M ) = 1 . Suppose we are given a symmetric Dirichlet form E in L 2 ( m ) . The associated Markovian semigroup is denoted by { T t } and we assume that T t 1 = 1 . Here, 1 stands for a constant function of M, taking the value 1. For any f L 1 , we use the notation
f = M f d m .
We also assume that 1 is the unique invariant function for the semigroup { T t } . Then, as t , we have
T t f f
in L 2 ( m ) . The semigroup { T t } is called ergodic when Equation (2) holds. We define the index γ 2 2 by
γ 2 2 = lim ¯ 1 t log T t m 2 2
which is often called the spectral gap (see e.g., Theorem 4.2.5 of [1]). Here, m stands for a linear operator f m ( f ) = f and 2 2 stands for the operator norm from L 2 ( m ) to L 2 ( m ) . In connection to this index, we are interested in another index γ Z Z defined by
γ Z Z = lim ¯ 1 t log T t m Z Z .
Here, Z is the Zygmund space (sometimes denoted by L log L ). The space Z is defined as follows. Set ϕ ( x ) = log ( 1 + x ) and Φ ( x ) = 0 x ϕ ( y ) d y . Then
Z = { f ; M Φ ( | f | ) d m < } .
We introduce norms in Z later (see Section 2).
On the other hand, the logarithmic Sobolev inequality is a powerful tool in the analysis of Markovian semigroups. The inequality takes the following form:
M f 2 ( x ) log ( | f ( x ) | / f 2 ) d m 1 γ LS E ( f , f ) .
Here, 2 stands for the L 2 -norm and the constant γ LS is chosen to be maximal and is called the logarithmic Sobolev constant. The form of the inequality reminds us the notion of entropy:
Ent ( f ) = E [ f log ( f / f ) ] .
An important application of the logarithmic Sobolev inequality is the following estimate of the entropy (see e.g., Chapter 6.1 of [2]):
Ent ( T t f ) e 2 γ LS t Ent f .
We are interested in the relation between γ Z Z and γ LS . In fact, we show the inequality γ Z Z γ LS . This kind of estimate of γ Z Z is given in [3], but in this paper we give a direct connection to the constant γ LS .
The organization of this paper is as follows. In Section 2, we give several kinds of norms in the Zygmund space Z. Using these notions, we show relations between the entropy and the norm in the space Z and give a proof of the main result. As an example, we discuss the Laguerre operator in Section 3. We give a precise expression of the resolvent kernel. In Section 4, we introduce Orlicz spaces ( L log β L ). We also discuss how to show the boundedness of operators in Orlicz spaces. Using these, we investigate the spectrum of the Laguerre operator in Orlicz spaces in Section 5. We can completely determine the spectrum and can see the relation between the spectral gap and the logarithmic Sobolev constant.

2. Entropy and the Zygmund Space

2.1. The Zygmund Space

We start with the Zygmund space. Let ( M , B , m ) be a measure space, and we assume that m ( M ) = 1 , i.e., m is a probability measure. All functions in the paper are assumed to be B -measurable. We denote the integration with respect to m by f . Of course, we assume the integrability of a function f. We also use the notation E [ f ] for f .
The Zygmund space is the set of all measurable functions f with E [ | f | log | f | ] < . We denote it by Z or L log L . We can define a norm in this space. To do this, we introduce a function ϕ on [ 0 , ) defined by
ϕ ( x ) = log ( 1 + x )
and further, we define
Φ ( x ) = 0 x ϕ ( y ) d y = ( 1 + x ) log ( 1 + x ) x .
Φ is a convex function. Now, define N Φ by
N Φ ( f ) = inf { λ ; E [ Φ ( | f | / λ ) ] 1 } .
This norm is sometimes called the Luxemburg norm (see e.g., [4]). The norm of the constant function 1 can be computed as
N Φ ( 1 ) = inf { λ ; E [ Φ ( 1 / λ ) ] 1 } = inf { λ ; 1 / λ Φ 1 ( 1 ) } = inf { λ ; 1 / λ e 1 } ( Φ 1 ( 1 ) = e 1 ) = 1 e 1 .
Z becomes a Banach with the norm N Φ .
The dual space of Z is given as follows. Let ψ be the inverse function of ϕ , i.e.,
ψ ( x ) = e x 1 .
Using this, define Ψ by
Ψ ( x ) = 0 x ψ ( y ) d y = 0 x ( e y 1 ) d y = e x x 1 .
The dual space of Z can be identified with the space of all measurable functions f with E [ Ψ ( ε | f | ) ] < for some ε > 0 (see [4]).
The following inequality is fundamental:
x y Φ ( x ) + Ψ ( y ) .
Using this, we can show that
f 1 ( e 1 ) N Φ ( f ) .
In fact, if N Φ ( f ) = 1 , we have
E [ | f | y ] E [ Φ ( | f | ) ] + E [ Ψ ( y ) ] = 1 + Ψ ( y ) = e y y .
Hence,
E [ | f | ] e y y 1 .
The right-hand side takes its minimum e 1 when y = 1 . Hence, we obtain Equation (13). This shows that Z L 1 .
Further, we have
N Φ ( f f ) N Φ ( f ) + N Φ ( f ) = N Φ ( f ) + | f | N Φ ( 1 ) N Φ ( f ) + f 1 1 e 1 N Φ ( f ) + N Φ ( f ) ( Equation ( 13 ) ) = 2 N Φ ( f ) .

2.2. Entropy

Now, we recall the notion of entropy. In this section, all functions are taken from Z. For any non-negative function f, the entropy Ent ( f ) is defined by
Ent ( f ) = E [ f log ( f / f ) ] .
We will discuss the relation between N Φ ( f ) and Ent ( f ) . First, we show the following.
Proposition 1.
For any non-negative function f, we have
f E [ Φ ( | ( f f ) / f | ) ] Ent ( f ) .
Proof. 
We note the following inequality
Φ ( | x 1 | ) x log x x + 1 .
Using this, we can get
E [ Φ ( | ( f / f ) 1 | ) E [ ( f / f ) log ( f / f ) ( f / f ) + 1 ] = 1 f Ent ( f ) ,
which is the desired result. ☐
If, in addition, we assume f 1 , we can get another estimate.
Proposition 2.
If a non-negative function f satisfies f 1 , then we have
E [ Φ ( | f f | ) ] f Ent ( f ) .
Proof. 
Let us show the inequality
Φ ( | x f | ) f ( x log x x log f x + f )
for any x 0 . Set
F ( x ) = f ( x log x x log f x + f ) Φ ( | x f | ) .
(1) The case x f .
By the definition,
F ( x ) = f ( x log x x log f x + f ) ( 1 + x f ) log ( 1 + x f ) + x f .
Hence, F ( f ) = 0 . By differentiating the function F, we have
F ( x ) = f ( log x log f ) log ( 1 + x f )
and so we easily see that F ( f ) = 0 . The second-order derivative is given by
F ( x ) = f x 1 1 + x f = ( f 1 ) ( x f ) x ( 1 + x f ) 0 .
Thus, we have F ( x ) 0 for x f .
(2) The case x f .
In this case, we have
F ( x ) = f ( x log x x log f x + f ) ( 1 + f x ) log ( 1 + f x ) + f x .
So F ( f ) = 0 is clear. The derivative of F is
F ( x ) = f ( log x log f ) + log ( 1 + f x )
and so we easily see that F ( f ) = 0 . Furthermore, we have
F ( x ) = f x 1 1 + f x = f + f 2 f x x x ( 1 + f x ) = ( f + 1 ) ( f x ) x ( 1 + x f ) 0 .
Thus, we have F ( x ) 0 and F ( x ) 0 for x f .
Using the inequality Equation (18), we have
E [ Φ ( | f f | ) ] f E [ f log f f log f f + f ) ] f Ent ( f )
which completes the proof. ☐
Now we are ready to show that the N Φ -norm is dominated by the entropy.
Proposition 3.
For any non-negative function f, we have
N Φ ( f f ) max { f , Ent ( f ) } Ent ( f ) .
Proof. 
We note that since Φ is convex and Φ ( 0 ) = 0 , Φ satisfies that for J 1
Φ ( J x ) J Φ ( x )
and, for ε 1
Φ ( ε x ) ε Φ ( x ) .
The proof of (19) is divided into two cases.
(1) The case f N Φ ( f f ) .
Set N = N Φ ( f f ) . Applying Proposition 1 to the function f N , we have
Ent ( f N ) f / N E [ Φ ( | ( ( f / N ) f / N ) | f / N ) ] f / N 1 f / N E [ Φ ( | ( ( f / N ) f / N ) | ) ] = E [ Φ ( | f f | N ) ] = 1 .
We used Equation (20) in the second line. Thus we have, in this case,
N Φ ( f f ) Ent ( f ) .
(2) The case f N Φ ( f f ) .
Set N = N Φ ( f f ) . Since f N 1 , we can apply Proposition 2 to f N and obtain
E [ Φ ( | ( f / N ) f / N ) | ) ] f N Ent ( f N ) .
Now, N 2 f Ent ( f ) follows since the left-hand side equals 1. Hence,
N Φ ( f f ) f 1 / 2 Ent ( f ) 1 / 2 .
Combining both of them, we have
N Φ ( f f ) max { Ent ( f ) 1 / 2 , f 1 / 2 } Ent ( f ) 1 / 2
which completes the proof. ☐
In turn, we prove the inequality of the reversed direction.
Proposition 4.
We have the following inequality:
Ent ( f ) f log ( 4 / e ) E [ Φ ( | ( f f ) / f | ) ] .
Proof. 
Note that
n x log x x + 1 = Φ ( | x 1 | ) for x 1 ,
x log x x + 1 C Φ ( | x 1 | ) for x [ 0 , 1 ] .
Here, C = 1 log ( 4 / e ) . In fact, Equation (23) is clear and so we only show Equation (24). For x [ 0 , 1 ] , we have
Φ ( | x 1 | ) = Φ ( 1 x ) = ( 1 + 1 x ) log ( 1 + 1 x ) ( 1 x ) = ( 2 x ) log ( 2 x ) + x 1 .
We set
f ( x ) = C { ( 2 x ) log ( 2 x ) + x 1 } x log x + x 1 .
Then
f ( 0 ) = C ( 2 log 2 1 ) 1 = 0 , f ( 1 ) = 0 .
It is not hard to show that f ( x ) 0 for x [ 0 , 1 ] . Hence, we have Equation (24).
Since C > 1 , we have, for all x 0 ,
x log x x + 1 C Φ ( | x 1 | ) .
Substituting x = f f in it and integrating both hands, we have
E f f log f f f f + 1 C E Φ f f 1 .
Hence,
E [ f log ( f / f ) ] C f E [ Φ ( | ( f f ) / f | ] .
This completes the proof. ☐
When f 1 , we can show another inequality.
Proposition 5.
If we assume f 1 , then we have
Ent ( f ) E [ Φ ( | f f | ) ] + 2 .
Proof. 
Since f 1 , we easily see
f log f ( 1 + | f f | ) log ( 1 + | f f | ) .
Using this, we have
Φ ( | f f | ) = ( 1 + | f f | ) log ( 1 + | f f | ) f f f log f f f f log f f log f + f log f f f .
Integrating both hands, we have
E [ Φ ( | f f | ) ] Ent ( f ) + f log f 2 f .
Now set g ( x ) = x log x + 2 x on [ 0 , 1 ] . Then, g ( x ) = log x 1 + 2 = 1 log x 0 for x [ 0 , 1 ] . Hence, g takes its maximum at x = 1 . Therefore,
f log f 2 f g ( 1 ) = 2 .
Thus we get the desired result. ☐
We are ready to show that the entropy is dominated by the N Φ -norm.
Proposition 6.
We have the following inequality:
Ent ( f ) 3 N Φ ( f f ) .
Proof. 
Since the function Φ is convex, Φ satisfies following inequality. For J 1 , we have
Φ ( J x ) J Φ ( x ) ,
and, for ε 1 ,
Φ ( ε x ) ε Φ ( x ) .
The proof is divided into two cases.
(1) The case f N Φ ( f f ) 1 .
For notational simplicity, we denote N Φ ( f f ) by N. Using Proposition 4 for f N ,
Ent ( f N ) f / N log ( 4 / e ) E [ Φ ( | ( ( f / N ) f / N ) | f / N ) ] f / N log ( 4 / e ) 1 f / N E [ Φ ( | ( ( f / N ) f / N ) | ) ] 1 log ( 4 / e ) .
Here we used Equation (29) in the second line.
(2) The case f N Φ ( f f ) 1 .
This time we use Proposition 5 and obtain
Ent ( f N ) E [ Φ ( | ( ( f / N ) f / N ) | ) ] + 2 = 1 + 2 = 3 .
Since 1 log ( 4 / e ) 3 , we have Equation (27).
Let us recall the logarithmic Sobolev inequality:
M f 2 ( x ) log ( f ( x ) 2 / f 2 2 ) d m 2 γ LS E ( f , f ) ,
which yields the following entropy estimate:
Ent ( T t f ) e 2 γ LS t Ent f .
Now, we are in a position to prove the following main theorem.
Theorem 1.
We have the following inequality.
γ LS γ Z Z
Proof. 
We may assume γ LS > 0 . Let f be a non-negative function. If N Φ ( f ) 1 , then we have for sufficiently large t
N Φ ( T t f f ) Ent ( T t f ) ( Ent ( T t f ) f ) e γ LS t Ent ( f ) ( e γ LS t Ent ( f ) f ) ( Equation ( 31 ) ) e γ LS t 3 N Φ ( f f ) ( e γ LS t 3 N Φ ( f f ) f ) e γ LS t 6 N Φ ( f ) ( e γ LS t 6 N Φ ( f ) ( e 1 ) N Φ ( f ) ) e γ LS t 6 ( e 1 ) .
Next, we take a general f. If N Φ ( f ) 1 , then N Φ ( f + ), NΦ(f ) ≤ NΦ(|f|) = NΦ(f) ≤ 1 and so
N Φ ( T t f f ) = N Φ ( T t f + f + T t f + f ) N Φ ( T t f + f ) + N Φ ( T t f f ) e γ LS t 24 ( e 1 ) .
Therefore, we have
T t m Z Z 24 ( e 1 ) e γ LS t .
Hence, this completes the proof. ☐
In Theorem 1, we have shown that γ Z Z γ LS . We now connect the Logarithmic Sobolev constant γ LS and the spectral gap. Let us denote the set of spectrum of A in the Zygmund space Z by σ ( A Z ) . Then, the following inequality is known (see e.g., Chapter IV, Proposition 2.2 of Engel-Nagel [5])
sup { λ ; λ σ ( A Z ) } lim t 1 t log T t Z Z .
If we restrict ourselves to the mean 0 functions, we have
sup { λ ; λ σ ( A Z ) \ { 0 } } lim t 1 t log T t m Z Z = γ Z Z .
Now we set
ζ SG = sup { λ ; λ σ ( A Z ) \ { 0 } }
and call it the spectral gap in Z. Hence, we have the following inequalities:
γ LS γ Z Z ζ SG .
Example 1.
Let us consider the Ornstein–Uhlenbeck operator A = d 2 d x 2 x d d x on R . The reference measure is m ( d x ) = 1 2 π e x 2 / 2 d x and the Dirichlet form is given by
E ( f , g ) = R f ( x ) g ( x ) m ( d x ) .
In this case, it is known that γ LS = 1 . Moreover f ( x ) = x is an eigenfunction for the eigenvalue 1 . Hence, we have
1 ζ SG γ Z Z γ LS = 1
which shows γ LS = γ Z Z = ζ SG . In Section 5, we will give an example that ζ SG > γ LS holds.

3. Spectrum of the Laguerre Operator

In this section, we give an example.

3.1. The Laguerre Operator

We consider the following operator:
A = x d d x 2 + ( α + 1 x ) d d x
Since eigenfunctions of A are Laguerre polynomials (see e.g., Lebedev [6]), we call the diffusion process generated by A the Laguerre process as in [7]. It is also an interest rate model called the Cox-Ingersoll-Ross process in mathematical finance.
We assume that α > 1 . This is necessary to ensure that the invariant measure becomes a probability measure. Actually, the invariant probability measure is given by
m ( d x ) = 1 Γ ( α + 1 ) x α e x d x ,
which is the gamma distribution of the parameters α + 1 , 1.
There is another characterization of a diffusion process by a speed measure and a scale function. In our case, setting
ρ ( x ) = 1 Γ ( α + 1 ) x α e x ,
the speed measure is m = ρ ( x ) d x and the scale function s is determined by d s = 1 x ρ ( x ) d x . Following Feller, the boundary 0 is classified as a non-exit, an entrance when α 0 and exit, and an entrance when 1 < α < 0 . We impose the Neumann boundary condition when 1 < α < 0 to ensure that the associated diffusion process is conservative.
We can give the associated Dirichlet form E as
E ( u , v ) = 0 d u d s d v d s d s .
Here, d u d s = x ρ ( x ) d u d x . Therefore,
0 d u d s d v d s d s = 0 x ρ ( x ) d u d x x ρ ( x ) d v d x 1 x ρ ( x ) d x = 0 x d u d x d v d x ρ ( x ) d x = 0 x d u d x d v d x m ( d x ) .
This means that
E ( u , v ) = 0 x u ( x ) v ( x ) m ( d x ) .
It is well-known that the set of the spectrum of A in L 2 ( m ) is Z + and eigenfunctions are Laguerre polynomials. We also give an expression of the resolvent. To do this, we need confluent hypergeometric functions.

3.2. Confluent Hypergeometric Functions

We recall confluent hypergeometric functions (see, e.g., Beals-Wong [8] or Lebedev [6]). They are defined by
1 F 1 ( a ; c ; x ) = n = 0 ( a ) n ( c ) n n ! x n .
Here, ( a ) n is the Pochhammer symbol, i.e.,
( a ) n = Γ ( a + n ) Γ ( a ) = a ( a + 1 ) ( a + n 1 ) n 1 1 n = 0
A function defied by Equation (38) converges for all x C and is an analytic function. This function satisfies the following differential equation:
x u + ( c x ) u = a u .
This equation is called the Kummer equation (or the confluent hypergeometric equation), and, of course, is closely related to our generator A in Equation (35). Our interest is in the spectrum of A , and so confluent hypergeometric functions are candidates of eigenfunctions. If 1 F 1 belongs to L 2 , it is an eigenfunction and it is so when a = n , n Z + . In this case, 1 F 1 (−n; c; x) is nothing but a Laguerre polynomial (up to constant) and is an eigenfunction. For simplicity, we introduce the following notation:
M ( a , 1 + α ; x ) = 1 F 1 ( a ; 1 + α ; x ) .
By the way, Equation (42) is a second-order differential equation; there is another independent solution, which is given by
U ( a , 1 + α ; x ) = Γ ( α ) Γ ( a α ) M ( a , 1 + α ; x ) + Γ ( α ) Γ ( a ) x α M ( a α , 1 α ; x ) .
This function is called a confluent hypergeometric function of the second kind. Their Wronskian is
W ( M ( a , 1 + α ; · ) U ( a , 1 + α ; · ) ) ( x ) = Γ ( 1 + α ) Γ ( a ) x α 1 e x .
The Laguerre polynomial is written as
L n α ( x ) = ( α + 1 ) n n ! M ( n , α + 1 ; x ) .
Our parameter α is chosen to be consistent with the parameter of the Laguerre polynomial. The asymptotic behavior of these function is given as follows (see e.g., Lebedev [6]). When x 0 ,
M ( a , 1 + α ; x ) 1 ,
U ( a , 1 + α ; x ) Γ ( α ) Γ ( a ) x α .
However, when α = 0 , x α should be replaced by log x .
When x ,
M ( a , 1 + α ; x ) Γ ( 1 + α ) Γ ( a ) e x x a 1 α ,
U ( a , 1 + α ; x ) x a .
Here, we assumed a, 1 + α 0 , 1 , 2 , .
Now we can give an expression of the resolvent. Recall that we assumed α > 0 and a 0 , 1 , 2 , . The resolvent G a = ( a A ) 1 has the following kernel expression:
G a f ( x ) = 0 G a ( x , y ) f ( y ) d y
where
G a ( x , y ) = M ( a , 1 + α ; y ) U ( a , 1 + α , x ) 1 p ( y ) W ( y ) y < x , M ( a , 1 + α ; x ) U ( a , 1 + α , y ) 1 p ( y ) W ( y ) y > x .
Here, W stands for the Wronskian in Equation (45) and p ( y ) = y . Hence, we have
G a ( x , y ) = Γ ( a ) Γ ( 1 + α ) M ( a , 1 + α ; y ) U ( a , 1 + α , x ) e y y α y < x , Γ ( a ) Γ ( 1 + α ) M ( a , 1 + α ; x ) U ( a , 1 + α , y ) e y y α y > x .
G a is a bounded operator in L 2 ( m ) if a 0 , 1 , 2 , . We will discuss later what happens in the Zygmund space.

3.3. The Logarithmic Sobolev Inequality

We show that the logarithmic Sobolev inequality holds for the Laguerre operator A . You can also see the result in Chapter 2.7.3 of [1] from the view point of the curvature dimension condition. Recall that the Dirichlet form associated with A is given by Equation (39).
Theorem 2.
We assume that α > 1 2 . Then, the following logarithmic Sobolev inequality holds for the Dirichlet form E in (36):
0 u 2 log ( | u | / u 2 ) ν ( d x ) 2 E ( u , u ) .
Proof. 
It is enough to check Bakry-Emery’s Γ 2 -criterion. It is as follows. From Equation (36), the square field Γ is given by
Γ ( f , g ) = x f ( x ) g ( x ) .
The generator is A u = x u + ( α x ) u . Hence, the Γ 2 is computed as
2 Γ 2 ( u , u ) = A Γ ( u , u ) 2 Γ ( A u , u ) = A ( x u 2 ) 2 x ( A u ) u = x ( x u 2 ) + ( α x ) ( x u 2 ) 2 x ( x u + ( α x ) u ) u = x ( u 2 + 2 x u u ) + ( α x ) ( u 2 + 2 x u u ) 2 x u ( u + x u u + ( α x ) u ) = x ( 2 u u + 2 u u + 2 x u 2 + 2 x u u ) + ( α x ) ( u 2 + 2 x u u ) 2 x ( u u + x u u u 2 + ( α x ) u u ) = 2 x u u + 2 x 2 u 2 + ( α + x ) u 2 = 2 ( x u + 1 2 u ) 2 1 2 u 2 + ( α + x ) u 2 = 2 ( x u + 1 2 u ) 2 + ( 2 α 1 2 x + 1 ) x u 2 .
Thus, we have
Γ 2 ( u , u ) = ( x u + 1 2 u ) 2 + 1 2 ( 2 α 1 2 x + 1 ) Γ ( u , u ) .
From this we have Γ 2 ( u , u ) 1 2 Γ ( u , u ) under the condition α 1 2 . Due to Bakry-Emery’s Γ 2 -criterion, this implies that γ LS 1 2 . ☐
Taking f ( x ) = e ξ x , we can see that γ LS = 1 2 is the best constant.
Remark 1.
This result was shown in Korzeniowski-Stroock [7] when α = 1 . In that paper, it was emphasized that the logarithmic Sobolev constant differs from the spectral gap.

4. Orlicz Space L log β L

We start with the definition of the Orlicz space. Take any β > 0 and fix it. We introduce a norm in the space of all functions f with E [ | f | log β ( 1 + | f | ) ] < . Define a function ϕ on [ 0 , ) by
ϕ ( x ) = log ( 1 + x ) .
Then, further define
Φ β ( x ) = 0 x ϕ β ( y ) d y = 0 x log β ( 1 + y ) d y .
Φ β is a concave function. To get the behavior of Φ β at , we use the l’Hospital theorem and get
lim x ( x ) x log β ( 1 + x ) = lim x Φ β ( x ) ( x log β ( 1 + x ) ) = lim x log β ( 1 + x ) log β ( 1 + x ) + β x log β 1 ( 1 + x ) 1 1 + x = lim x 1 1 + β x ( 1 + x ) log ( 1 + x ) = 1 .
Therefore, when x , we can see
Φ β ( x ) x log β ( 1 + x ) .
We define the space L log β L by
L log β L = { f ; E [ Φ β ( | f | ) ] < } .
Then, L log β L becomes a Banach space with the norm N Φ β defined by
N Φ β ( f ) = inf { λ ; E [ Φ β ( | f | / λ ) ] 1 } .
For instance, the norm of the constant function 1 is
N Φ β ( 1 ) = inf { λ ; E [ Φ ( 1 / λ ) ] 1 } = inf { λ ; 1 / λ Φ 1 ( 1 ) } = 1 Φ 1 ( 1 ) .
If β = 1 , then Φ 1 ( 1 ) = e 1 . In the sequel, the operator norm of linear operators from L log β L into L log β L is defined by using the norm N Φ β .

4.1. Dual Space

The dual space of L log β L is characterized as follows. Let ψ β be the inverse function of log β ( 1 + x ) :
ψ β ( x ) = e x 1 / β 1 .
Further, we define
Ψ β ( x ) = 0 x ψ β ( y ) d y .
The Orlicz space associated with Ψ β is the dual space of L log β L . Let us study the asymptotic behavior of Ψ β at x = .
Proposition 7.
We have the following:
Ψ β ( x ) β e x 1 / β x ( β 1 ) / β a s x .
Proof. 
We use the l’Hôspital theorem.
lim x Ψ β ( x ) e x 1 / β x ( β 1 ) / β = lim x Ψ β ( x ) ( e x 1 / β x ( β 1 ) / β ) = lim x e x 1 / β 1 e x 1 / β ( 1 / β ) x 1 / β 1 x ( β 1 ) / β + e x 1 / β β 1 β x 1 / β = lim x e x 1 / β 1 ( e x 1 / β / β ) + e x 1 / β β 1 β x 1 / β = lim x 1 e x 1 / β 1 β + β 1 β x 1 / β = β .
Equation (63) easily follows from this. ☐
The following Hausdorff-Young inequality plays a fundamental role in the later computation.
x y Φ β ( x ) + Ψ β ( y ) .
For example, if N Φ ( f ) = 1 , then we can show that for y > 0
E [ | f | y ] E [ Φ β ( | f | ) ] + E [ Ψ β ( y ) ] = 1 + Ψ β ( y ) .
Hence,
E [ | f | ] 1 + Ψ β ( y ) y .
This shows that L log β L L 1 and there exists a constant κ β > 0 so that
f 1 κ β N Φ β ( f ) .

4.2. Linear Operators in Orlicz Spaces

Orlicz space L log β L is a Banach space with the norm N Φ β . The operator norm can also be defined in terms of this norm. However, since this norm is hard to calculate concretely, we take another way.
We introduce a new norm Φ , which is called the Orlicz norm, by
f Φ β = sup { E [ g | f | ] ; E [ Ψ β ( g ) ] 1 } .
Here, g runs over all functions satisfying E [ Ψ β ( g ) ] 1 . Replacing f with f / 2 ,
f / 2 Φ β = sup { E [ g | f / 2 | ] ; E [ Ψ β ( g ) ] 1 } = sup { E [ ( g / 2 ) | f | ] ; E [ Ψ β ( 2 g / 2 ) ] 1 } = sup { E [ g | f | ] ; E [ Ψ β ( 2 g ) ] 1 } .
Hence, we can rewrite Equation (65) as follows:
f Φ β = 2 sup { E [ g | f | ] ; E [ Ψ β ( 2 g ) ] 1 } .
We will rewrite the condition E [ Ψ β ( 2 g ) ] 1 .
Proposition 8.
We have
sup { E [ e g 1 / β ] ; g 0 and E [ Ψ β ( 2 g ) ] 1 } < .
Proof. 
From Proposition 7, we have
Ψ β ( 2 x ) β e ( 2 x ) 1 / β ( 2 x ) ( β 1 ) / β .
We can take large constant C β so that
e x 1 / β Ψ β ( 2 x ) + C β .
Therefore, if E [ Ψ β ( 2 g ) ] 1 , then we have
E [ e g 1 / β ] E [ Ψ β ( 2 g ) ] + C β 1 + C β ,
which is the desired result. ☐
We set
K β = sup { E [ e g 1 / β ] ; g 0 and E [ Ψ β ( 2 g ) ] 1 } .
Then, by Proposition 8, we can see that E [ e g 1 / β ] K β if E [ Ψ β ( 2 g ) ] 1 .
It is well-known that two norms, N Φ β and Φ β , are equivalent (see e.g., p. 61, Chapter III. 3.3, Proposition 4 of Rao-Ren [4]):
N Φ β ( f ) f Φ β 2 N Φ β ( f ) .
From this, we have the following
Proposition 9.
A linear operator T on L log β L is bounded if there exist positive constants A and B so that
T f Φ β A E [ Φ β ( | f | ) ] + B .
Proof. 
If we assume Equation (68), then we have
T f / N Φ β ( f ) Φ β A E [ Φ β ( | f | / N Φ β ( f ) ) ] + B = A + B ,
which implies T f Φ β ( A + B ) N Φ β ( f ) . The rest is easy from Equation (67).
Corollary 1.
Let K β be a constant defined by (65). Then a linear operator T on L log β L is bounded if there exist constants A and B so that for any non-negative function g with E [ e g 1 / β ] K β and any non-negative function f L log β L , we have
E [ g | T f | ] A E [ | f | log + β | f | ] + B .
Proof. 
From Equation (58), we have Φ ( x ) x log β ( x + 1 ) and hence we can find constants a and b so that
x log + β x a Φ β ( x ) + b .
Then
T f Φ β = 2 sup { E [ g | T f | ] ; E [ Ψ β ( 2 g ) ] 1 } 2 sup { E [ g | T f | ] ; E [ e g 1 / β K β } 2 { A E [ | f | log + β | f | ] + B } ( Equation ( 69 ) ) 2 A E [ | f | log + β | f | ] + 2 B 2 A E [ a Φ β ( | f | ) + b ] + 2 B 2 a A E [ Φ β ( | f | ) ] + 2 b A + 2 B .
Now, from Proposition 9, T is bounded. ☐
We list up some inequalities which are necessary later. For x, y 0 ,
( x + y ) p 2 p 1 ( a p + b p ) , p 1 , ( x + y ) p a p + b p , p 1 .
There exists a positive constant A β so that
x y A β ( x log + β x + e y 1 / β ) .
This inequality is a modification of the following Hausdorff-Young inequality:
x y x log x x + e y .

5. The Spectrum of the Laguerre Operator in L log β L

The kernel representation of the resolvent of Laguerre operator was given in Equation (53). It is bounded in L 2 . We will examine whether it is bounded in L log β L . Recall that our reference measure is m = 1 Γ ( α + 1 ) x α e x d x . We assume α > 1 . From now on, we ignore the constant and consider the measure x α e x d x .
We take any f L log β L . We also take a non-negative function g satisfying E [ e g 1 / β ] K β . Our aim is to show that there exist constants A and B so that E [ g | G a f | ] A E [ | f | log + β | f | ] + B . The integrability is important, and we do not need the precise constant. Hence, we use the following notation:
x y .
This means that there exist constants k and l so that
x k y + l .
Here, constants k and l are independent of functions f and g. This is important but we do not mention this each time.
We starts with an estimate of the defective Gamma function.
Proposition 10.
Take any k R . If k < 0 , then
y x k e x d x y k e y , y 0 .
If k 0 , then there exists a constant c k so that
y x k e x d x c k ( y + 1 ) k e y , y 0 .
These are easily obtained by seeing that y x k e x d x y k e y d x as y .
Proposition 11.
Assume that κ > 0 and λ > 1 . Then, there exists a constant C depending on κ and λ so that
0 y ( log x ) κ x λ d x C ( log y + 1 ) κ y λ + 1 , y 1 .
Proof. 
This inequality can be reduced to the previous proposition. By the change of variable formula, we have
0 y ( log x ) κ x λ d x = ( λ + 1 ) log y u λ + 1 κ e u λ + 1 λ e u λ + 1 d u λ + 1 u = ( λ + 1 ) log x , x = e u / ( λ + 1 ) , d x = e u / ( λ + 1 ) d u λ + 1 = 1 λ + 1 κ + 1 ( λ + 1 ) log y u κ e u d u C 1 λ + 1 κ + 1 ( ( λ + 1 ) log y + 1 ) κ e ( λ + 1 ) log y ( Equation ( 72 ) ) c 1 λ + 1 log y + 1 λ + 1 κ y λ + 1 c 2 λ + 1 ( log y + 1 ) κ y λ + 1 .
This completes the proof. ☐
We study integrals involving function g.
Proposition 12.
For any k R , α > 1 and β > 0 , there exist constants C 1 , C 2 so that
y g ( x ) x k x α e x d x C 1 e y y k y e g ( x ) 1 / β x α e x d x + C 2 y k + β + α e y , y 1 .
We have assumed that E [ e g 1 / β ] K β , so we have that there exists a constant c depending only on k, β, and α so that
y g ( x ) x k e α e x d x c y k y ( β + α ) 0 e y , y 1 .
Proof. 
Set F ( x ) = x k e x , x 0 . Then
F ( x ) = k x k 1 e x e x x k = ( k x ) x k 1 e x
F takes its maximum at x = k and for x k , F is decreasing. Hence, if k y x , then
x k e x y k e y
and if 1 y k , then for y x , we have
x k e x F ( k ) F ( k ) F ( 1 ) 1 F ( 1 ) F ( k ) F ( 1 ) 1 F ( y ) F ( k ) F ( 1 ) 1 y k e y .
Hence, there exists a constant c k depending on k such that for 1 y x ,
x k e x c k y k e y .
Therefore
y g ( x ) x k x α e x d x = y ( g ( x ) e x ) e x x k x α e x d x A β y ( e g ( x ) 1 / β + e x log + β e x ) e x x k x α e x d x ( Equation ( 70 ) ) A β y e g ( x ) 1 / β ( x k e x ) x α e x d x + A β y x k + β x α e x d x A β c k y k e y y e g ( x ) 1 / β x α e x d x + A β c k + β + α y k + β y α e y . ( Proposition 10 and Equation ( 76 ) )
This completes the proof. ☐
Proposition 13.
Assume k + β + α + 1 > 0 . Then, there exists a constant C so that
1 y g ( x ) x k x α d x y k + β + α + 1 C 1 + C 2 1 y e g ( x ) 1 / β x α e x d x , y 1 .
Recalling that E [ e g 1 / β ] K β , we have
1 y g ( x ) x k x α d x C 3 y k + β + α + 1 , y 1 .
When k + β + α + 1 = 0 , we have
1 y g ( x ) x k x α d x C 1 log y + C 2 1 y e g ( x ) 1 / β x α e x d x , y 1 .
Again by E [ e g 1 / β ] K β ,
1 y g ( x ) x k x α d x C 1 log y + C 2 , y 1 .
Proof. 
By using Equation (70), we have
1 y g ( x ) x k x α d x = 1 y ( e x g ( x ) ) e x x k x α d x A β 1 y { e g ( x ) 1 / β + e x log + β e x } x k x α e x d x ( Equation ( 70 ) ) A β 1 y { x k + β + α + e g ( x ) 1 / β x k + α e x } d x A β 1 k + β + α + 1 [ x k + β + α + 1 ] 1 y + A β 1 y e g ( x ) 1 / β x k + β + α + 1 x β α 1 x α e x d x A β 1 k + β + α + 1 y k + β + α + 1 + A β y k + β + α + 1 1 y e g ( x ) 1 / β x α e x d x .
When k + β + α + 1 = 0 , in the computation above, we just need to note that the primitive function of x 1 is log x . ☐
Of course, when k + β + α + 1 < 0 , the left-hand side of Equation (77) is bounded.
We have seen the asymptotic behavior of integrals as y . We can also get the asymptotics as y 0 . We will prove this by reducing to the previous case.
Proposition 14.
Suppose that α > 1 , β > 0 . Then, there exist constants C 1 , C 2 so that
0 y g ( x ) x α d x C 1 y α + 1 0 y e g ( x ) 1 / β x α d x + C 2 y α + 1 ( log y + 1 ) β , y 1 .
By using E [ e g 1 / β ] A β , we have
0 y g ( x ) x α d x C y α + 1 ( log y + 1 ) β y 1 .
Proof. 
For y 1 , we have
0 y g ( x ) x α d x = 0 y ( g ( x ) x α 1 ) x α + 1 x α d x A β 0 y ( e g ( x ) 1 / β + x α 1 log + β x α 1 ) x α + 1 x α d x ( Equation ( 70 ) ) = A β 0 y e g ( x ) 1 / β x α + 1 x α d x + A β 0 y x α 1 ( α + 1 ) ( log x ) β x α + 1 x α d x A β y α + 1 0 y e g ( x ) 1 / β x α d x + A β ( α + 1 ) 0 y ( log x ) β x α d x A β y α + 1 0 y e g ( x ) 1 / β x α d x + A β ( α + 1 ) C ( log y + 1 ) β y α + 1 , ( Equation ( 73 ) )
which is the desired result. ☐
Lastly, we show the estimate involving the function f.
Proposition 15.
We have the following inequality for f.
1 x β f ( x ) x α e x d x C 1 + C 2 1 f ( x ) log + β f ( x ) x α e x d x .
Proof. 
From Young’s inequality
1 x β f ( x ) x α e x d x = 2 β 1 ( x / 2 ) β f ( x ) x α e x d x 2 β A β 1 { e ( ( x / 2 ) β ) 1 / β + f ( x ) log + β f ( x ) } x α e x d x ( Equation ( 70 ) ) 2 β A β 1 e x / 2 x α e x d x + 2 β A β 1 f ( x ) log + β f ( x ) x α e x d x 2 β A β 1 x α e x / 2 d x + 2 β A β 1 f ( x ) log + β f ( x ) x α e x d x C 1 + C 2 1 f ( x ) log + β f ( x ) x α e x d x .
 ☐
We now investigate the spectrum. Let us start with the point spectrum.
Theorem 3.
The point spectrum of A is { z ; z < β } { 0 , 1 , 2 , , [ β ] } . Here, [ β ] stands for the integer part of β.
Proof. 
We show that a is an eigenvalue of A if a < β . We have seen that M ( x ) = M ( a , α + 1 ; x ) satisfies the differential equation A M = a M (see Equation (42)). We only need to show that M L log β L . The integrability of Φ β ( M ) on [ 0 , 1 ] is trivial since M is bounded on [ 0 , 1 ] . We see the integrability on [ 1 , ) :
1 | M ( x ) | log + β | M ( x ) | x α e x d x 1 e x x a 1 α ( x + | a 1 α | log + x ) β x α e x d x 1 x a 1 x β d x 1 x a 1 + β d x < .
The finiteness in the last line follows from a 1 + β < 1 .
It remains to be shown that n ( 0 n < β ) is an eigenvalue. In fact, M ( n , α + 1 ; x ) is an polynomial of order n (a Laguerre polynomial up to normalization) and hence the integrability of Φ β ( M ) follows easily.
Theorem 4.
If a > β , then a belongs to the resolvent set except when a = 0 , 1 , 2 , , [ β ] .
Proof. 
We show that G a in (51) is bounded in L log β L . To do this, we use Corollary 1. We recall the kernel G a ( x , y ) .
(1) The case y < x .
Let us consider the resolvent kernel M ( y ) U ( x ) y α e y (recall Equation (53) but we ignore a constant multiplication). In the region of x 1 , we have
0 1 g ( x ) x α e x d x 0 x | M ( y ) U ( x ) | | f ( y ) | y α e y d y 0 1 g ( x ) x α e x d x 0 x x α | f ( y ) | y α e y d y ( M is bounded and | U | x α on [ 0 , 1 ] ) 0 1 g ( x ) x α x α e x d x 0 1 | f ( y ) | y α e y d y 0 1 { e g ( x ) 1 / β + x α α β ( log x ) β } x α e x d x 0 1 | f ( y ) | log + β | f ( y ) | y α e y d y ( x x log + β x and Equation ( 70 ) ) 0 1 | f ( y ) | log + β | f ( y ) | y α e y d y .
When α = 0 , | U | log x in the above computation should be changed as follows:
0 1 g ( x ) e x d x 0 x | M ( y ) U ( x ) | | f ( y ) | e y d y 0 1 g ( x ) e x d x 0 x | log x | | f ( y ) | e y d y 0 1 g ( x ) | log x | e x d x 0 1 | f ( y ) | e y d y 0 1 { e g ( x ) 1 / β + | log x | log + β | log x | } e x d x 0 1 | f ( y ) | log + β | f ( y ) | e y d y ( x x log + β x and Equation ( 70 ) ) 0 1 | f ( y ) | log + β | f ( y ) | e y d y .
In the region of x 1 , we have
1 g ( x ) x α e x d x 0 x | M ( y ) U ( x ) | | f ( y ) | y α e y d y 1 x a g ( x ) x α e x d x 0 x | M ( y ) | f ( y ) | y α e y d y ( | U | x a on [ 1 , ) ] ) = 1 x a g ( x ) x α e x d x 0 1 | M ( y ) | | f ( y ) | y α e y d y + 1 x a g ( x ) x α e x d x 1 x | M ( y ) | | f ( y ) | y α e y d y 1 x a g ( x ) x α e x d x 0 1 | f ( y ) | y α e y d y + 1 x a g ( x ) x α e x d x 1 x e y y a 1 α | f ( y ) | y α e y d y ( | M | is bounded on [ 0 , 1 ] and | M | e y y a 1 α on [ 1 , ) ) 1 x a g ( x ) x α e x d x 0 1 | f ( y ) | y α e y d y + 1 y a 1 | f ( y ) | d y y x a g ( x ) x α e x d x ( Fubini s theorem ) 0 1 | f ( y ) | y α e y d y + 1 y a 1 | f ( y ) | y a y ( α + β ) 0 e y d y ( Equation ( 75 ) ) 0 1 | f ( y ) | y α e y d y + 1 | f ( y ) | y 1 y ( α + β ) 0 e y d y 0 1 | f ( y ) | y α e y d y + 1 | f ( y ) | y α + β e y d y ( ( α + β ) 0 < α + β + 1 ) 0 1 | f ( y ) | log + β | f ( y ) | y α e y d y + 1 | f ( y ) | log + β | f ( y ) | y α e y d y . ( x x log + β x and Equation ( 83 ) )
(2) The case y > x .
We consider the resolvent kernel M ( x ) U ( y ) y α e y . In the region x 1 , we have
0 1 g ( x ) x α e x d x x | M ( x ) U ( y ) | | f ( y ) | y α e y d y 0 1 g ( x ) x α e x d x x 1 | U ( y ) | | f ( y ) | y α e y d y + 0 1 g ( x ) x α e x d x 1 | U ( y ) | | f ( y ) | y α e y d y 0 1 g ( x ) x α e x d x x 1 y α | f ( y ) | y α e y d y + 0 1 g ( x ) x α e x d x 1 y a | f ( y ) | y α e y d y ( | U | y α on [ 0 , 1 ] , | U | y a on [ 1 , ) ] ) 0 1 | f ( y ) | e y d y 0 y g ( x ) x α e x d x + 1 y a | f ( y ) | y α e y d y ( Fubini s theorem and Equation ( 82 ) ) for y = 1 ) 0 1 | f ( y ) | e y y α + 1 ( log y + 1 ) β d y + 1 y β | f ( y ) | y α e y d y ( Equation ( 82 ) ) and a > β ) 0 1 | f ( y ) | y α e y d y + 1 | f ( y ) | log + β | f ( y ) | y α e y d y ( y ( log y + 1 ) β is bounded on [ 0 , 1 ] and Equation ( 83 ) ) 0 1 | f ( y ) | log + β | f ( y ) | y α e y d y + 1 | f ( y ) | log + β | f ( y ) | y α e y d y ( x x log + β x ) .
In the region x 1 , we have
1 g ( x ) x α e x d x x | M ( x ) U ( y ) | | f ( y ) | y α e y d y 1 g ( x ) x α e x d x x e x x a 1 α y a | f ( y ) | y α e y d y ( | M | e x x a 1 α and | U | y a on [ 1 , ) ) 1 y a y α | f ( y ) | e y d y 1 y g ( x ) x a 1 d x ( Fubini s theorem ) 1 y a y α | f ( y ) | e y y a 1 + β + 1 d y ( Equation ( 78 ) ) 1 y β | f ( y ) | y α e y d y 1 | f ( y ) | log + β | f ( y ) | y α e y d y ( Equation ( 83 ) )
Thus, we have shown that G a is bounded in L log β L . Hence, the spectrum of A is completely determined. ☐
The spectrum is shown as in Figure 1. The case of β = 1 is the Zygmund space. Hence, the spectral gap equals 1. So we have ζ SG = 1 where ζ SG is defined by Equation (33). Therefore, by Equation (34), we have
1 2 = γ LS γ Z Z ζ SG = 1
and so it shows that γ LS ζ SG may happen. This is a well-known result in the case of L 2 . Furthermore, in [9], we have shown that assuming the logarithmic Sobolev inequality, the spectra in L p ( 1 < p < ) are all the same. In our case, the spectrum in L 2 is 0, 1 , 2 , The spectrum in L log β L is quite different from that. Moreover, the logarithmic Sobolev inequality may not give a uniform estimate for spectral gaps among the Orlicz spaces L log β L ( β > 0 ).

Acknowledgments

This research was partially supported by the Japanese Ministry of Education, Culture, Sports, Science, and Technology, Grant-in-Aid for Scientific Research (B), No. 15H03624, 2017. The author would like to thank referees for their careful reading of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bakry, D.; Gentil, I.; Ledoux, M. Analysis and Geometry of Markov Diffusion Operators; Grundlehren der Mathematischen Wissenschaften 348; Springer: Cham, Switzerland, 2014. [Google Scholar]
  2. Deuschel, J.-D.; Stroock, D.W. Large Deviations; Academic Press: San Diego, CA, USA, 1989. [Google Scholar]
  3. Roberto, C.; Zegarliński, B. Orlicz-Sobolev inequalities for sub-Gaussian measures and ergodicity of Markov semi-groups. J. Funct. Anal. 2007, 243, 28–66. [Google Scholar] [CrossRef]
  4. Rao, M.M.; Ren, Z.D. Theory of Orlicz Spaces; Monographs and Textbooks in Pure and Applied Mathematics, 146; Marcel Dekker: New York, NY, USA, 1991. [Google Scholar]
  5. Engel, K.-J.; Nagel, R. One-Parameter Semigroups for Linear Evolution Equations; Springer: New York, NY, USA, 2000. [Google Scholar]
  6. Lebedev, N.N. Special Functions and Their Applications; Dover: Mineola, NY, USA, 1972. [Google Scholar]
  7. Korzeniowski, A.; Stroock, D.W. An example in the theory of hypercontractive semigroups. Proc. Am. Math. Soc. 1985, 94, 87–90. [Google Scholar] [CrossRef]
  8. Beals, R.; Wong, R. Special Functions; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  9. Kusuoka, S.; Shigekawa, I. Exponential convergence of Markovian semigroups and their spectra on Lp-spaces. Kyoto J. Math. 2014, 54, 367–399. [Google Scholar] [CrossRef]
Figure 1. The spectrum in the space L log β L .
Figure 1. The spectrum in the space L log β L .
Entropy 20 00220 g001

Share and Cite

MDPI and ACS Style

Shigekawa, I. Logarithmic Sobolev Inequality and Exponential Convergence of a Markovian Semigroup in the Zygmund Space. Entropy 2018, 20, 220. https://doi.org/10.3390/e20040220

AMA Style

Shigekawa I. Logarithmic Sobolev Inequality and Exponential Convergence of a Markovian Semigroup in the Zygmund Space. Entropy. 2018; 20(4):220. https://doi.org/10.3390/e20040220

Chicago/Turabian Style

Shigekawa, Ichiro. 2018. "Logarithmic Sobolev Inequality and Exponential Convergence of a Markovian Semigroup in the Zygmund Space" Entropy 20, no. 4: 220. https://doi.org/10.3390/e20040220

APA Style

Shigekawa, I. (2018). Logarithmic Sobolev Inequality and Exponential Convergence of a Markovian Semigroup in the Zygmund Space. Entropy, 20(4), 220. https://doi.org/10.3390/e20040220

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