On free Generalized Inverse Gaussian distributions

We study here properties of free Generalized Inverse Gaussian distributions (fGIG) in free probability. We show that in many cases the fGIG shares similar properties with the classical GIG distribution. In particular we prove that fGIG is freely infinitely divisible, free regular and unimodal, and moreover we determine which distributions in this class are freely selfdecomposable. In the second part of the paper we prove that for free random variables $X,Y$ where $Y$ has a free Poisson distribution one has $X\stackrel{d}{=}\frac{1}{X+Y}$ if and only if $X$ has fGIG distribution for special choice of parameters. We also point out that the free GIG distribution maximizes the same free entropy functional as the classical GIG does for the classical entropy.


Introduction
Free probability was introduced by Voiculescu in [Voi85] as a non-commutative probability theory where one defines a new notion of independence, so called freeness or free independence. Non-commutative probability is a counterpart of the classical probability theory where one allows random variables to be non-commutative objects. Instead of defining a probability space as a triplet (Ω, F , P) we switch to a pair (A, ϕ) where A is an algebra of random variables and ϕ : A → C is a linear functional, in classical situation ϕ = E. It is natural then to consider algebras A where random variables do not commute (for example C * or W * -algebras). For bounded random variables independence can be equivalently understood as a rule of calculating mixed moments. It turns out that while for commuting random variables only one such rule leads to a meaningful notion of independence, the non-commutative setting is richer and one can consider several notions of independence. Free independence seems to be the one which is the most important. The precise definition of freeness is stated in Section 2 below.
Free probability emerged from questions related to operator algebras however the development of this theory showed that it is surprisingly closely related with the classical probability theory. First evidence of such relations appeared with Voiculescu's results about asymptotic freeness of random matrices. Asymptotic freeness roughly speaking states that (classically) independent, unitarily invariant random matrices, when size goes to infinity, become free.
Another link between free and classical probability goes via infinite divisibility. With a notion of independence in hand one can consider a convolution of probability measures related to this notion. For free independence such operation is called free convolution and it is denoted by ⊞. More precisely for free random variables X, Y with respective distributions µ, ν the distribution of the sum X + Y is called the free convolution of µ and ν and is denoted by µ ⊞ ν. The next natural step is to ask which probability measures are infinitely divisible with respect to this convolution. We say that µ is freely infinitely divisible if for any n ≥ 1 there exists a probability measure µ n such that µ = µ n ⊞ . . . ⊞ µ n n times .
Here we come across another striking relation between free and classical probability: there exists a bijection between classically and freely infinitely divisible probability measures, this bijection was found in [BP99] and it is called Bercovici-Pata (BP) bijection. This bijection has number of interesting properties, for example measures in bijection have the same domains of attraction. In free probability literature it is standard approach to look for the free counterpart of a classical distribution via BP bijection. For example Wigner's semicircle law plays the role of the Gaussian law in free analogue of Central Limit Theorem, Marchenko-Pastur distribution appears in the limit of free version of Poisson limit theorem and is often called free Poisson distribution.
While BP bijection proved to be a powerful tool, it does not preserve all good properties of distributions. Consider for example Lukacs theorem which says that for classically independent random variables X, Y random variables X + Y and X/(X + Y ) are independent if and only if X, Y have gamma distribution with the same scale parameter [Luk55]. One can consider similar problem in free probability and gets the following result (see [Szp15,Szp16]) for free random variables X, Y random variables X + Y and (X + Y ) −1/2 X(X + Y ) −1/2 are free if and only if X, Y have Marchenko-Pastur (free Poisson) distribution with the same rate. From this example one can see our point -it is not the image under BP bijection of the Gamma distribution (studied in [PAS08,HT14]), which has the Lukacs independence property in free probability, but in this context free Poisson distribution plays the role of the classical Gamma distribution.
In [Szp17] another free independence property was studied -a free version of so called Matsumoto-Yor property (see [MY01,LW00]). In classical probability this property says that for independent X, Y random variables 1/(X + Y ) and 1/X − 1/(X + Y ) are independent if and only if X has a Generalized Inverse Gaussian (GIG) distribution and Y has a Gamma distribution. In the free version of this theorem (i.e. the theorem where one replaces classical independence assumptions by free independence) it turns out that the role of the Gamma distribution is taken again by the free Poisson distribution and the role of the GIG distribution plays a probability measure which appeared for the first time in [Fer06]. We will refer to this measure as the free Generalized Inverse Gaussian distribution or fGIG for short. We give the definition of this distribution in Section 2.
The main motivation of this paper is to study further properties of fGIG distribution. The results from [Szp17] suggest that in some sense (but not by means of the BP bijection) this distribution is the free probability analogue of the classical GIG distribution. It is natural then to ask if fGIG distribution shares more properties with its classical counterpart. It is known that the classical GIG distribution is infinitely divisible (see [BNH77]) and selfdecomposable (see [Hal79,SS79]). In [LS83] the GIG distribution was characterized in terms of an equality in distribution, namely if we take X, Y 1 , Y 2 independent and such that Y 1 and Y 2 have Gamma distributions with suitable parameters and we assume that then X necessarily has a GIG distribution. A simpler version of this theorem characterizes smaller class of fGIG distributions by equality for X and Y 1 as described above.
The overall result of this paper is that the two distributions GIG and fGIG indeed have many similarities. We show that fGIG distribution is freely infinitely divisible and even more that it is free regular. Moreover fGIG distribution can be characterized by the equality in distribution (1.2), where one has to replace the independence assumption by freeness and assume that Y 1 has free Poisson distributions. While there are only several examples of freely selfdecomposable distributions it is interesting to ask whether fGIG has this property. It turns out that selfdecomposability is the point where the symmetry between GIG and fGIG partially breaks down: not all fGIG distributions are freely selfdecomposable. We find conditions on the parameters of fGIG family for which this distributions are freely selfdecomposable. Except from the results mentioned above we prove that fGIG distribution is unimodal. We also point out that in [Fer06] it was proved that fGIG maximizes a certain free entropy functional. An easy application of Gibbs' inequality shows that the classical GIG maximizes the same functional of classical entropy.
The paper is organized as follows: In Section 2 we shortly recall basics of free probability and next we study some properties of fGIG distributions. Section 3 is devoted to the study of free infinite divisibility, free regularity, free selfdecomposability and unimodality of the fGIG distribution. In Section 4 we show that the free counterpart of the characterization of GIG distribution by (1.2) holds true, and we discuss entropy analogies between GIG and fGIG.

Free GIG distributions
In this section we recall the definition of free GIG distribution and study basic properties of this distribution. In particular we study in detail the R-transform of fGIG distribution. Some of the properties established in this section will be crucial in the subsequent sections where we study free infinite divisibility of the free GIG distribution and characterization of the free GIG distribution. The free GIG distribution appeared for the first time (not under the name free GIG) as the almost sure weak limit of empirical spectral distribution of GIG matrices (see [Fer06]).

Basics of free probability
This paper deals mainly with properties of free GIG distribution related to free probability and in particular to free convolution. Therefore in this section we introduce notions and tools that we need in this paper. The introduction is far from being detailed, reader not familiar with free probability may find a very good introduction to the theory in [VDN92,NS06,MS17].
2 o Let I be an index set. A family of subalgebras (A i ) i∈I are called free if ϕ(X 1 · · · X n ) = 0 whenever a i ∈ A j i , j 1 = j 2 = . . . = j n and ϕ(X i ) = 0 for all i = 1, . . . , n and n = 1, 2, . . .. Similarly, self-adjoint random variables X, Y ∈ A are free (freely independent) when subalgebras generated by (X, 1 A ) and (Y, 1 A ) are freely independent.
3 o The distribution of a self-adjoint random variable is identified via moments, that is for a random variable X we say that a probability measure µ is the distribution of X if ϕ(X n ) = t n dµ(t), for all n = 1, 2, . . .
Note that since we assume that our algebra A is a C * -algebra, all random variables are bounded, thus the sequence of moments indeed determines a unique probability measure.
4 o The distribution of the sum X + Y for free random variables X, Y with respective distributions µ and ν is called the free convolution of µ and ν, and is denoted by µ ⊞ ν.

Free GIG distribution
In this paper we are concerned with a specific family of probability measures which we will refer to as free GIG (fGIG) distributions.
Definition 2.1. The free Generalized Inverse Gaussian (fGIG) distribution is a measure µ = µ(α, β, λ), where λ ∈ R and α, β > 0 which is compactly supported on the interval [a, b] with the density where 0 < a < b are the solution of Observe that the system of equations for coefficients for fixed λ ∈ R and α, β > 0 has a unique solution 0 < a < b. We can easily get the following Remark 2.2. Let λ ∈ R. Given α, β > 0, the system of equations (2.1), (2.2) has a unique solution (a, b) such that Conversely, given (a, b) satisfying (2.3), the set of equations (2.1)-(2.2) has a unique solution (α, β), which is given by Thus we may parametrize fGIG distribution using parameters (a, b, λ) satisfying (2.3) instead of (α, β, λ). We will make it clear whenever we will use a parametrization different than (α, β, λ).

R-transform of fGIG distribution
The R-transform of the measure µ(α, β, λ) was calculated in [Szp17]. Since the R-transform will play a crucial role in the paper we devote this section for a detailed study of its properties. We also point out some properties of fGIG distribution which are derived from properties of the R-transform. Before we present the R-transform of fGIG distribution let us briefly recall how the Rtransform is defined and stress its importance for free probability.
Remark 2.4. 1 o For a probability measure µ one defines its Cauchy transform via It is an analytic function on the upper-half plane with values in the lower half-plane. Cauchy transform determines uniquely the measure and there is an inversion formula called Stieltjes inversion formula, namely for h ε (t) = − 1 π Im G µ (t + iε) one has where the limit is taken in the weak topology.
2 o For a compactly supported measure µ one can define in a neighbourhood of the origin so called R-transform by where by G −1 µ we denote the inverse under composition of the Cauchy transform of µ. The relevance of the R-transform for free probability comes form the fact that it linearizes free convolution, that is R µ⊞ν = R µ + R ν in a neighbourhood of zero.
The R-transform of fGIG distribution is given by in a neighbourhood of 0, where the square root is the principal value, .
Note that z = 0 is a removable singular point of r α,β,λ . Observe that in terms of A, B defined by (2.6) we have It is straightforward to observe that (2.7) implies A(λA + B) < 2AB < 2B 2 , thus we have γ < 0.
The following remark was used in [Szp17, Remark 2.1] without a proof. We give a proof here.
Since all roots of f α,β,λ are real, the square root f α,β,λ (z) may be defined continuously on C − ∪ R so that f α,β,λ (0) = α. As noted above δ < 0, and continuity of f α,β,λ implies that we have where we take the principal value of the square root in the expression 4β(η − z). Thus finally we arrive at the following form of the R-transform which is analytic in C − and continuous in (C − ∪ R) \ {α} as required.
Next we describe the behaviour of the R-transform around the singular point z = α.
Corollary 2.8. In the case λ < 0 one can extend r α,β,λ to an analytic function in C − and continuous in C − ∪ R.

Some properties of fGIG distribution
We study here further properties of free GIG distribution. Some of them motivate Section 4 where we will characterize fGIG distribution in a way analogous to classical GIG distribution. The next remark recalls the definition and some basic facts about free Poisson distribution, which will play an important role in this paper. Remark 2.9. 1 o Marchenko-Pastur (or free-Poisson) distribution ν = ν(γ, λ) is defined by the formula where γ, λ > 0 and the measureν, supported on the interval (γ(1 − √ λ) 2 , γ(1 + √ λ) 2 ), has the density (with respect to the Lebesgue measure) The next proposition was proved in [Szp17, Remark 2.1] which is the free counterpart of a convolution property of classical Gamma and GIG distribution. The proof is a straightforward calculation of the R-transform with the help of Remark 2.5.
We also quote another result from [Szp17, Remark 2.2] which is again the free analogue of a property of classical GIG distribution. The proof is a simple calculation of the density.
The two propositions above imply some distributional properties of fGIG distribution. In the Section 4 we will study characterization of the fGIG distribution related to these properties.
3 o Both identities above can be iterated finitely many times, so that one obtains that X d = Y 1 + (Y 2 + · · · ) −1 −1 , where Y 1 , Y 2 , . . . are free, for k odd Y k has the free Poisson distribution ν(1/β, λ) and for k even Y k has the distribution ν(1/α, λ). For the case described in 1 o one simply has to take α = β. We are not sure if infinite continued fractions can be defined.
Remark 2.14. This result reflects the fact that GIG matrix generalizes the Wishart matrix for λ ≥ 1, but not for λ < 1 (see [Fer06] for GIG matrix and [HP00] for the Wishart matrix).
Corollary 2.15. Considering the continuous dependence of roots on parameters shows the following asymptotic behaviour of the double root δ < 0 and the simple root η ≥ α.

Regularity of fGIG distribution under free convolution
In this section we study in detail regularity properties of the fGIG distribution related to the operation of free additive convolution. In the next theorem we collect all the results proved in this section. The theorem contains several statements about free GIG distributions. Each subsection of the present section proves a part of the theorem.
2 o The free Levy measure is of the form 3 o It is free regular with zero drift for all α, β > 0 and λ ∈ R.

Free infinite divisibility and free Lévy measure
As we mentioned before, having the operation of free convolution defined, it is natural to study infinite divisibility with respect to ⊞. We say that µ is freely infinitely divisible if for any n ≥ 1 there exists a probability measure µ n such that µ = µ n ⊞ . . . ⊞ µ n n times . It turns out that free infinite divisibility of compactly supported measures can by described in terms of analytic properties of the R-transform. In particular it was proved in [Voi86,Theorem 4.3] that the free infinite divisibility is equivalent to the inequality Im(r α,β,λ (z)) ≤ 0 for all z ∈ C − .
As in the classical case, for freely infinitely divisible probability measures, one can represent its free cumulant transform with a Lévy-Khintchine type formula. For a probability measure µ on R, the free cumulant transform is defined by Then µ is FID if and only if C ⊞ µ can be analytically extended to C − via the formula where ξ ∈ R, ζ ≥ 0 and τ is a measure on R such that The triplet (ξ, ζ, τ ) is called the free characteristic triplet of µ, and τ is called the free Lévy measure of µ. The formula (3.3) is called the free Lévy-Khintchine formula.
Remark 3.2. The above form of free Lévy-Khintchine formula was obtained by Barndorff-Nielsen and Thorbjørnsen [BNT02b] and it has a probabilistic interpretation (see [Sat13]). Another form was obtained by Bercovici and Voiculescu [BV93], which is more suitable for limit theorems.
In order to prove that all fGIG distributions are freely infinitely divisible we will use the following lemma.
Proof. Since f is analytic the function Imf is harmonic and thus satisfies the maximum principle. Fix ε > 0. Since f (z) → 0 uniformly with z → ∞, let R > 0 be such that Imf (z) < ε. Consider a domain D ε with the boundary Observe that on ∂D ε Imf (z) < ε by assumptions, and hence by the maximum principle we have Imf (z) < ε on whole D ε . Letting ε → 0 we get that Imf (z) ≤ 0 on C − .
Next we proceed with the proof of free infinite divisibility of fGIG distributions.
In this case similar argument shows that µ(α, β, λ) is FID. Moreover by (2.16) point z = α is a removable singularity and r α,β,λ extends to a continuous function on C − ∪ R. Thus one does not need to take care of the behaviour around z = α.
Case 3 λ = 0. For λ = 0 one can adopt a similar argumentation using (2.17). It also follows from the fact that free GIG family µ(α, β, λ) is weakly continuous with respect to λ. Since free infinite divisibility is preserved by weak limits, then the case λ = 0 may be deduced from the previous two cases.
Next we will determine the free Lévy measure of free GIG distribution µ(α, β, λ).

Free regularity
In this subsection we will deal with a property stronger than free infinite divisibility, so called free regularity. Let µ be a FID distribution with the free characteristic triplet (ξ, ζ, τ ). When the semicircular part ζ is zero and the free Lévy measure τ satisfies a stronger integrability property R min{1, |x|}τ (dx) < ∞, then the free Lévy-Khintchine representation reduces to The distribution µ is said to be free regular [PAS12] if ξ ′ ≥ 0 and τ is supported on (0, ∞). A probability measure µ on R is free regular if and only if the free convolution power µ ⊞t is supported on [0, ∞) for every t > 0, see [AHS13]. Examples of free regular distributions include positive free stable distributions, free Poisson distributions and powers of free Poisson distributions [Has16]. and free selfadjoint elements X 1 , . . . , X n which have symmetric FID distributions, the polynomial n i,j=1 a ij X i X j has a free regular distribution with zero drift. Proof of Theorem 3.1 3 o . For the free GIG distributions, the semicircular part can be found by ζ α,β,λ = lim z→∞ z −1 r α,β,λ (z) = 0. The free Lévy measure (3.1) satisfies and so we have the reduced formula (3.10). The drift is given by ξ ′ α,β,λ = lim u→−∞ r α,β,λ (u) = 0.

Free selfdecomposability
Classical GIG distribution is selfdecomposable [Hal79,SS79] (more strongly, hyperbolically completely monotone [Bon92, p. 74]), and hence it is natural to ask whether free GIG distribution is freely selfdecomposable. A distribution µ is said to be freely selfdecomposable (FSD) [BNT02a] if for any c ∈ (0, 1) there exists a probability measure µ c such that µ = (D c µ)⊞µ c , where D c µ is the dilation of µ, namely is non-increasing on (0, 1/η). The derivative is (3.14) Hence FSD is equivalent to Using η ≥ α > 0 > δ, one can show that 2αη − 2ηδ + αδ > 0, a straightforward calculation shows that the function g takes a minimum at a point in (0, 1/η). Thus FSD is equivalent to In order to determine when the above inequality holds, it is convenient to switch to parameters A, B defined by (2.6). Using formulas derived in Section 2.3 we obtain (3.17) Calculating λ for which D is non-positive we obtain that Corollary 3.4. One can easily find that the maximum of the function − B

Unimodality
Since relations of free infinite divisibility and free self decomposability were studied in the literature, we decided to determine whether measures from the free GIG family are unimodal.
A measure µ is said to be unimodal if for some c ∈ R is non-decreasing on (−∞, c) and non-increasing on (c, ∞). In this case c is called the mode. Hasebe and Thorbjørnsen [HT16] proved that FSD distributions are unimodal. Since some free GIG distributions are not FSD, the result from [HT16] does not apply. However it turns out that free GIG measures are unimodal.
Proof of Theorem 3.1 5 o . Calculating the derivative of the density of µ(α, β, λ) one obtains Denoting by f (x) the quadratic polynomial in the numerator, one can easily see from the shape of the density that f (a) > 0 > f (b) and hence the derivative vanishes at a unique point in (a, b) (since f is quadratic).

Characterizations the free GIG distribution
In this section we show that the fGIG distribution can be characterized similarly as classical GIG distribution. In [Szp17] fGIG was characterized in terms of free independence property, the classical probability analogue of this result characterizes classical GIG distribution. In this section we find two more instances where such analogy holds true, one is a characterization by some distributional properties related with continued fractions, the other is maximization of free entropy.

Continued fraction characterization
In this section we study a characterization of fGIG distribution which is analogous to the characterization of GIG distribution proved in [LS83]. Our strategy is different from the one used in [LS83]. We will not deal with continued fractions, but we will take advantage of subordination for free convolutions, which allows us to prove the simpler version of "continued fraction" characterization of fGIG distribution.
Theorem 4.1. Let Y have the free Poisson distribution ν(1/α, λ) and let X be free from Y , where α, λ > 0 and X > 0, then we have if and only if X has free GIG distribution µ(α, α, −λ).
Remark 4.2. Observe that the "if" part of the above theorems is contained in the remark 2.12. We only have to show that if (4.1) holds where Y has free Poisson distribution ν(1/α, λ), then X has free GIG distribution.
As mentioned above our proof of the above theorem uses subordination of free convolution. This property of free convolution was first observed by Voiculescu [Voi93] and then generalized by Biane [Bia98]. Let us shortly recall what we mean by subordination of free additive convolution.
has the desired property, which in particular means that for X distributed µ(α, α, λ) equation (4.7) is satisfied. Thus if there is a unique solution it has to be the Cauchy transform of the free GIG distribution.
To prove uniqueness of the Cauchy transform of X, we will prove that coefficients of the expansion of G X at a special "good" point, are uniquely determined by α and λ.
First we will determine the point at which we will expand the function. Observe that with our assumptions G X −1 is well defined on the negative half-line, moreover G X −1 (x) < 0 for any x < 0, and we have G X −1 (x) → 0 with x → −∞. On the other hand the function f (x) = 1/x−x is decreasing on the negative half-line, and negative for x ∈ (−1, 0). Thus there exist a unique point c ∈ (−1, 0) such that Let us denote , (4.10) where the last equality follows from (4.8).
One has N(c) = c, and our functional equation (4.7) may be rewritten (with the help of (4.8)) as − M(z) + z = z 2 M(N(z)). Observe that β 0 = c since N(c) = c. Differentiating (4.10) we observe that any β n , n ≥ 1 is a rational function of α, λ, c, α 0 , α 1 , . . . , α n . Moreover any β n , n ≥ 1 is a degree one polynomial in α n . We have where R n is a rational function of n + 3 variables evaluated at (α, λ, c, α 0 , α 1 , . . . , α n−1 ), which does not depend on the distribution of X. For example β 1 is given by (4.13) Next we investigate some properties of c, α 0 and α 1 . Evaluating both sides of (4.11) at z = c yields −M(c) + c = c 2 M(N(c)) = c 2 M(c), (4.14) Observe that α 0 = M(c) = G X (1/c) and α 1 = M ′ (c) = −c −2 G ′ X (1/c) hence we have where µ X is the distribution of X. Using the Schwarz inequality for the first estimate and a simple observation that 0 ≤ 1/(1 − cx) ≤ 1 for x > 0, for the latter estimate we obtain (4.15) The equation (4.9) together with (4.8) gives (4.16) Substituting (4.14) to (4.16) after simple calculations we get We start by showing that α 0 is determined only by α and λ. We will show that c, which we showed before is a unique number, depends only on α and λ and thus (4.14) shows that α 0 is determined by α and λ.
Since the polynomial c 4 − (1 + λ)c 3 is non-negative for c < 0 and has a root at c = 0, and the polynomial (λ − 1)c + α equals α > 0 at c = 0 it follows that there is only one negative c, such that the two polynomials are equal and thus the number c is uniquely determined by (α, λ). From (4.14) we see that α 0 is also uniquely determined by (α, λ).
Next we will prove that α 1 only depends on α and λ. Differentiating (4.11) and evaluating at z = c we obtain 1 − α 1 = 2cα 0 + c 2 α 1 β 1 . (4.18) Substituting α 0 and λ from the equations (4.14) and (4.17) we simplify (4.13) and we get and then equation (4.18) may be expressed in the form The above is a degree 2 polynomial in α 1 , denote this polynomial by f , we have then Where the first inequality follows from the fact that c < 0. Since the coefficient c(1 + c 2 ) 2 is negative we conclude that f has one root in the interval (0, 1/(1 + c 2 )) and the other in (1/(1 + c 2 ), ∞). The inequality (4.15) implies that α 1 is the smaller root of f , which is a function of α and c and hence of α and λ.
Thus we have shown that, if a random variable X > 0 satisfies the functional equation (4.7) for fixed α > 0 and λ > 0, then the point c and all the coefficients α 0 , α 1 , α 2 , . . . of the series expansion of M(z) at z = c are determined only by α and λ. By analytic continuation, the Cauchy transform G X is determined uniquely by α and λ, so there is only one distribution of X for which this equation is satisfied.
Here we point out the classical analogue. The (classical) GIG distribution is the probability measure on (0, ∞) with the density (α/β) λ/2 2K λ (2 √ αβ) x λ−1 e −(αx+β/x) , α, β > 0, λ ∈ R, (4.23) where K λ is the modified Bessel function of the second kind. Note that this density is proportional to exp(−V α,β,λ (x)). Kawamura and Iwase [KI03] proved that the GIG distribution is a unique probability measure which maximizes the classical entropy with the same potential H α,β,λ (p) = − p(x) log p(x) dx − V α,β,λ (x)p(x) dx among all the probability density functions p on (0, ∞). This statement is slightly different from the original one [KI03, Theorem 2], and for the reader's convenience a short proof is given below. The proof is a straightforward application of the Gibbs' inequality − p(x) log p(x) dx ≤ − p(x) log q(x) dx, (4.24) for all probability density functions p and q, say on (0, ∞). Taking q to be the density (4.23) of the classical GIG distribution and computing log q(x), we obtain the inequality H α,β,λ (p) ≤ − log (α/β) λ/2 2K λ (2 √ αβ) . (4.25) Since the Gibbs inequality (4.24) becomes equality if and only if p = q, the equality in (4.25) holds if and only if p = q, as well.
Remark 4.4. From the above observation, it is tempting to investigate the map Ce −V (x) dx → the maximizer µ V of the free entropy functional I V with potential V , where C > 0 is a normalizing constant. Under some assumption on V , the free entropy functional I V is known to have a unique maximizer (see [ST97]) and so the above map is well defined. Note that the density function Ce −V (x) is the maximizer of the classical entropy functional with potential V , which follows from the same arguments as above. This map sends Gaussian to semicircle, gamma to free Poisson (when λ ≥ 1), and GIG to free GIG. More examples can be found in [ST97].