Arxiv:2108.10514V1 [Math-Ph] 24 Aug 2021 Iers Oblzangbssaitc,Bs-Isensaitc,A of States Statistics, of Number Bose-Einstein Stat the Statistics, Counting Respectively

arXiv:2108.10514v1 [math-ph] 24 Aug 2021 iers oBlzanGbssaitc,Bs-isensaitc,a of states statistics, of number Bose-Einstein stat The statistics, counting respectively. statistics Boltzmann-Gibbs of Fermi-Dirac enumerat to to rules related rise different closely give three is Traditionally, it that expected combinatorics. naturally is it states, hoishv o enapidi hi study. their can in problems applied combinato advanced been counting more not these hence have However methods, theories elementary by (8.40)-(8.42)].) solved [23, be e.g. (See atce ftype of particles by eas ttsia ehnc eeal novscutn numbers counting involves generally mechanics statistical Because NSM AHMTC EAE OTHE TO RELATED MATHEMATICS SOME ON h on fve fqatmmcaisadsaitclmechanics. statistical from and mechanics calculus quantum umbral of reexamine view to th of point On suggests the spec also etc. this trees, recursions, hand, topological counting other Eynard-Orantin operads, complex in curves laws, genera, tral group complex formal differ- theory, many including to bordism effects objects, Hall makes mathematical quantum the connection fractional ent this study relate hand to to de- one possible by bridge the it entropy On a maximum functions. as principle entropy it formed the use by and statistics calculus interpolation interpolating umbral the to connect Bose-Einste We statistics interpolating statistics. statistics Fermi-Dirac and of statistics space universal a troduce Abstract. NEPLTN STATISTICS INTERPOLATING i W W W oiae yfatoa unu aleet,w in- we effects, Hall quantum fractional by Motivated cuyn ru of group a occupying D F BG BE ( ( ( { { { n n n 1. i i i } } } = ) = ) = ) Introduction INZHOU JIAN Y 1 g n n i i n i i !( ! i . ( g g n i g g i i − i ! i !( + ttsi,rsetvl,given respectively, is, states g n n i i )! i − − , 1)! 1)! , n i in identical e - rial ive nd es of 2 JIAN ZHOU

In 1982, the fractional quantum Hall effect [45] was discovered by Tsui, Stormer and Gossard. This leads to many new progresses in statistical physics. Following the original work of Laughlin [29] describing 1 the ground state of the FQHE for a filling fraction ν = m (where m is odd), Halperin [21] conjectured and Arovas, Shrieffer and Wilczek [2] showed that the quasiparticles and quasiholes in quantum Hall systems not only have fractional charges, they also obey fractional statistics in the sense that they are anyons in the sense of Wilczek [46] (see also Leinnaas-Myrheim [30]). Suppose that we have two identical particles in two dimensions. Then when one particle is exchanged in a counter- clockwise manner with the other, the wavefunction can change by an arbitrary phases (1) ψ(~r , ~r ) eiθψ(~r , ~r ), 1 2 → 1 2 The special cases θ =0, π correspond to bosons and fermions, respectively. Particles with other values of the angle θ = απ are called anyons. Fractional statistics in this sense involves geometric phases [5], and it indicates the topological nature of the fractional quantum Hall system. In fact, fractional quantum Hall effect is an example of topological order, a notion introduced by Wen [47] in 1989. In 1991, Haldane [21] reformulated the concept of “fractional statistics” as a generalization of the Pauli exclusion principle. He considered the situations in which the number gi of independent states of a single particle of type i is finite and extensive. This means the quasiparticles and qausiholes are confined to a finite region, and their number is proportional to the size of region. The generalized Pauli exclusion principle proposed by Haldane is the following linear relation: (2) ∆g = g ∆n . i − ij j j X Let di be the value of gi when ni = 1, then the number of gi for general n are { } { i} (3) g = d + n 1 g (n δ ). i i i − − ij j − ij j X For state-counting purposes at fixed particle numbers the particles can be regarded as bosons with a Fock space of dimension gi, or fermions with a Fock space of dimension gi + ni 1, so the number of states for given n is given by the following formula− in Wu [49, (4)]: { i} (gi +(ni 1) gij(nj δij)]! (4) W ( n )= − − j − . { i} n !(g 1 g (n δ )]! i i i Pj ij j ij Y − − − P MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 3

When gij = βδij, it becomes (g +(n 1)(1 β)]! (5) W ( n )= i i − − . { i} n !(g 1 β(n 1)]! i i i i Y − − − According to the fundamental principles of quantum statistical mechanics, the grand partition function is given by

(6) Z = W ( n ) exp (µ ǫ )/kT , { i} i − i ni i X{ } X where T is the temperature, k is the Boltzmann constant, ǫi is the energy level of a particle of species i, and µi is the chemical potential for species i. Then one has

(7) Z = Zi, i Y where Zi is given by (g +(n 1)(1 β)]! (8) Z = i i − − exp (µ ǫ )/kT . i n !(g 1 β(n 1)]! i − i n i i i Xi − − − This leads to the notion of statistics interpolating between Bose and Fermi statistics [8]. We will refer to them as interpolating statistics. For a comparative study of various interpolating statistics, see [8]. As our point of departure, we will consider the universal space of all possible interpolating statistics and study some mathematics re-S lated to it. This is consistent with a common mathematical practice: Form a set by collecting all mathematical objects with the same properties, and then study the mathematical structures on this set. This is also inspired by Wilson’s renormalization theory in which he considered the space of all coupling constants and the dynamics on it. The basic mathematical objects in our study are formal power series, so ultimately the mathematical setting should be that of formal algebraic geometry. However, on the outset we will first establish a link between interpolating statistics to umbral calculus in combinatorics, and we will use this as a stepping stone to applications from the point of formal algebraic geometry in subsequent work. Before 1970s, umbral calculus refered to some powerful techniques to prove combinatorial identities introduced by Blissard in 1861. In the 1970s, Rota and his collaborators developed the mathematical foundations of umbral calculus by means of linear functionals on the space of polynomials [36, 35, 34]. Another approach to this theory is through the techniques of Hopf algebras [26]. Clearly a connection between umbral calculus and quantum mechanics was known to Rota and his 4 JIAN ZHOU collaborators, but this has not been fully elaborated. The link with interpolating statistics provides more stimulus to interpret umbral calculus from the point of view of quantum mechanics. In the text we will reexamine umbral calculus using representations of Heisenberg commutation relations, coherent states, Segal-Bargmann transform, Wick quantization, etc. An application of the link between the interpolating statistics and the umbral calculus is the study of interpolating statistics by the principle of maximum entropy. This principle was first expounded by Jaynes in two papers in 1957 [24, 25] where he emphasized a natural correspondence between statistical mechanics and information theory. Ac- cording to this principle, Boltzmann distribution can be derived by maximizing the Shannon entropy under the constraint of a fixed expectation value of energy. To derive the Bose-Eisnstein distribution and Fermi-Dirac distribution, or more generally, all the interpolating statistics, one needs to find suitable generalized entropy functions tobe maximized. In 1967, Havrda-Charvát introduced structural α-entropy [22] within information theory. In 1988, Tsallis [41] rediscovered this entropy and proposed to use it to study nonextensive statistical mechanics. For more information, see his monograph [44] published in 2009. The Tsallis q-entropy is defined by q 1 q 1 p p − 1 (9) S = − i i = pq i − . q 1 − i · 1 q i −P X − In 1994, Tsallis [42] proposed to use the q-logarithm function

1 q x − 1 (10) log x := − q 1 q − and the q-exponential function

1 (11) exp (x) := [1 + (1 q)x] 1−q q − to unify the expression in traditional Boltzmann-Gibbs statistical mechanics and his nonextensive statistical mechanics. This idea was generalized by Naudts [32] in 2002. He introduced more general deformed exponentials and logarithms in generalized statistical physics. See also his monograph [33] published in 2011. We make a connection between umbral calculus and the generalized entropy based on the following observation: In umbral calculus, one deals with the deformed power functions, and so their exponential generating functions can be regarded as deformed exponential functions. As a consequence of the connections of umbral calculus to both the MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 5 interpolating statistics and the generalized entropy, we establish a connection between interpolating statistics and generalized entropy. Our main result is that the generalized entropy function is essentially the Legendre transform of the free energy (cf. Theorem 6.1), as in ordinary thermodynamics. This demonstrates the surprising power of applying umbral calculus to study various aspects of interpolating statistics. Such applications will be further pursued in subsequent work. The rest of this paper is arranged as follows. In Section 2 we recall some deﬁnitions related to the statistics interpolating Bose-Einstein statistics and Fermi-Dirac statistics. In a very short Section 3 we introduce the space of interpolating statistics and a sequence of group structures on it. We also introduce two diﬀerent but related ways to associated spectral curves to interpolating statistics. We relate interpolating statistics to umbral calculus in Section 4. We give umbral interpretation of the group structure of the space of interpolating statistics. Next we relate umbral calculus to quantumS mechanics in Section 5. In particular, generating series of associated series in umbral calculus are interpreted as coherent states in quantum mechanics. In Section 6 we interpret such generating series in umbral calculus as deformed exponential functions and this motivates us to relate interpolating statistics to generalized entropy functions. Some concluding remarks are made in Section 7. We present many examples in Appen- dix A. Most of them are from Roman [34, Chapter 4], but we present some computations related to their relationship with generalized entropy. Some new examples are also presented there. Many sequences on The On-Line Encyclopedia of Integer Sequences [37] appear in these examples, indicating very rich combinatorial interpretations to be ex- plored.

2. Statistics Interpolating Bose and Fermi Statistics In this Section we recall the general framework for interpolations between Bose-Einstein and Fermi-Dirac statistics at the level of grand canonical partition functions as developed by Chaturvedi and Srini- vasan [8]. We will also deﬁne the entropy associated with the one- particle free enrrgy by Legendre transformation.

2.1. Grand canonical partition functions of noninteracting particles. Assume that for a system of N identical particles each of which can occupy n states corresponding to energies E < < E , the 1 ··· n 6 JIAN ZHOU canonical partition function has the following structure:

(12) Z (x)= W W xN1 xNn , N N1 ··· Nn 1 ··· n Pn N =N i=1Xi where x = x1,...,xn; xi = exp( βEi). This assumption means that there is no interactions between particles− at diﬀerent energy levels, and there are WNi particles at energy level Ei for each i = 1,...,n. Here W0 = 1 is always assumed. In most of our example, we will also assume that W1 = 1. Under the above assumption, the grand partition canonical function has the following factorized form:

n (13) (X ,...,X )= eNµβZ (x)= z(X ), Z 1 n N i N 0 i=1 X≥ Y where X = exp( β(E µ)) and i − i − ∞ n (14) z(X)= WnX n=0 X is called the one-particle partition function.

2.2. Examples. The Bose-Einstein, Maxwell-Boltzmann and Fermi- Dirac statistics are given by 1 (15) zBE (X)= , zFD(X)=1+ X, zBG(X)= eX , 1 X − respectively.

2.3. Occupation numbers. In the following expansion:

N k (16) z(X) = Wk(N)X , k 0 X≥ the coeﬃcients Wk(N) are the numbers of N particles occupying k states, and they are called the occupation numbers. From the fact that z(X)N1 z(X)N2 = z(X)N1+N2 , one gets the following recursion relations for occupation· numbers:

(17) Wk(N1 + N2)= Wi(N1)Wk i(N2). − i=0 X MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 7

2.4. Mean number of particles. The expression for the mean number of particles in the energy state Ei for this class of statistics reads ∂ ∂ (18) w(Xi)= Xi log (W )= Xi log z(Xi). ∂Xi Z ∂Xi It will be called the weight function in this paper. For Bose-Einstein, Maxwell-Boltzmann and Fermi-Dirac statistics we have X X (19) wBE(X)= , wFD(X)= , wBG(X)= X, 1 X 1+ X − β(E µ) respectively. When X = e− − , one gets the familiar expressions:

BE 1 FD 1 BG β(E µ) (20) w = , w = , w = e− − . eβ(E µ) 1 eβ(E µ) +1 − − − 2.5. Cluster coeﬃcients. We now expand the one-particle free energy as follows:

∞ w (21) F (X) := log z(X)= n Xn, n n=1 X and so

∞ n (22) w(X)= wnX . n=1 X The coefficients wn’s are called the cluster coefficients. One can easily find the following combinatorial formula expressing the Wn’s as polynomials in the wn’s:

mi wi (23) Wn = . imi m ! i i PXmii=n Y For example,

W1 = w1, w w2 W = 2 + 1 , 2 2 2 w w w w3 W = 3 + 2 1 + 1 . 3 3 2 3!

Conversely, the sequence wn n 1 is obtained from the sequence Wn n 1by: { } ≥ { } ≥ ∞ w ∞ (24) n Xn = log(1 + W Xn). n n n=1 n=1 X X 8 JIAN ZHOU

For example,

w1 = W1, W 2 w = W 1 , 2 2 − 2 W 3 w = W W W + 1 . 3 3 − 2 1 3

2.6. Entropy associated with the one-particle partition function. The entropy function associated to the one-particle partition function z = eF (x) is obtained by the Legendre transformation of the free energy F (X) as follows: dF (X) (25) H(X) := F (X) log X = F (X) log X w(X). − · d log X − · The entropy is also called the eﬀective action or one-particle irreducible partition correlation function. As usual, if one can invert the function w = w(X) to get X = X(w), then H can be expressed in terms of w:

(26) H = w log X(w)+ F (X(w)). − For example, we have

HBG = X log X + X = wBG log wBG + wBG, − − X 1 HBE = log X + log −1 X · 1 X − − wBE = wBE log + log(1 + wBE), − 1+ wBE X HFD = log X + log(1 + X) −1+ X · wFD 1 = wFD log + log . − · 1 wFD 1 wFD − −

3. Space of Interpolating Statistics and Its Group Structures As our point of departure, we will introduce in this Section the space of interpolating statistics, on it an involution that interchanges the Bose-EinsteinS statistics with the Fermi-Dirac statistics, and a sequence of group structures on . Furthermore, we deﬁne two versions of spectral curves associated withS elements in . S MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 9

3.1. Space of interpolating statistics. Since we have w1 = W1 = 1, w(X) is of the form

n (27) w(X)= X + wnX . n 2 X≥ Given any series of this form, one gets a statistics whose grand canonical partition function is given by (13). So we now introduce a space of interpolating statistics: (28) = w(X)= X + w Xn . S { n } n 2 X≥ 3.2. Dual cluster coefficients and dual statistics. Given any for- n mal sequence w = X + n 2 wnX , by Lagrangian inversion, one can show that ≥ ∈ S P n (29) X = w + wˆnw n 2 X≥ for some weighted homogeneous polynomialsw ˆ in w ,w ,w of de- n { 2 2 n} gree n 1, where deg wj = j 1. The coefficientsw ˆn will be called the dual cluster− coefficients. Therefore,− one can define the dual sequence of w by:

n (30)w ˆ(X)= X + wˆnX . n 2 X≥ The map w wˆ deﬁnes an involution σ on . It is clear that 7→ S wBG(X)= X is the only ﬁxed point of this involution. One can easily check that X X (31) wˆ (X)= = w (X), wˆ (X)= = w (X). BE 1+ X FD FD 1 X BE − In other words, the involution σ : interchanges the Bose- Einstein statistics with the Fermi-DiracS statistics,→ S and so we will refer to it as the generalized boson-fermion duality.

3.3. A sequence of group structures on . The construction in last subsection inspires us to introduce a sequenceS of group structures n on . Given two elements w(X) = X + n>1 wnX and v(X) = S n X + vnX in , one can deﬁne v w by composition of series: n>1 S ◦ P P(v w)(X) = w(X)+ v w(X)n ◦ n n 1 X≥ k k n = (X + wkX )+ vn(X + wkX ) . n>1 Xk>1 X Xk>1 10 JIANZHOU

Clearly, ( , ) is a group. For m 1, and w = X + w Xn, define S ◦ ≥ N>1 n (m) m m w (X)= X + n wnX P, n>1 X and define v w = v(m) w(m). ◦m ◦ We call it the the m-th group multiplication on . These group structures will not play a role in this paper. We expectS they are part of bigger symmetry structure on to be discovered in the future. S 3.4. Spectral curves associated with interpolating statistics. Consider the curve in (X, z)-plane defined by (32) z = z(X) or the curve in the (X,Y )-plane defined by: (33) eY = z(X). They will be called the spectral curves associated with the corresponding interpolating statistics. They are inspired by Eynard-Orantin topological recursion. See Appendix A.8 and Appendix A.10 for some interesting examples. 4. Interpolating Statistics and Umbral Calculus In this Section we relate interpolating statistics to umbral calculus. This makes powerful combinatorial results accessible to the study of interpolating statistics. We will present some examples in Section A. We will also give umbral interpretation of the group structure of the space of interpolating statistics introduced in §3.3. S 4.1. Occupation numbers as deformed binomial coefficients. Recall the generating series of the occupation numbers is given by: ∞ ∞ w W (N)Xk = z(X)N = eNF (X) = exp(N n Xn). k n n=1 Xk=0 X After expanding the right-hand side of the last equality as a series in X, one sees that each Wk(N) is a polynomial of degree k in N. N One can rewrite W (N) as and understand it as a deformed k k binomial coefficient. This is because the following identities clearly hold for n 0 by (17): ≥ x y x + y (34) = . i j n i+j=n X MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 11

This is understood as the deformed Chu-Vandermonde identities.

4.2. Occupation numbers, deformed power functions, and polynomial sequence of binomial type. Let us write the occupation numbers as follows: 1 (35) W (N)= γ (N), k k! k 1 where the factor k! is the Gibbs overcounting correction factor. Then we have by (17):

k k (36) γk(N1 + N2)= γi(N1)γk i(N2). i − i=0 X Recall a sequence of polynomials γn(x) is said to be of binomial type [28] if and only if { }

(37) γ1(x)= x, n n (38) γn(x + y)= γj(x)γn j(y) j − j=0 X for all x, y and n. So we are led to the theory of polynomial sequence of binomial type developed in the setting of umbral calculus [34].

4.3. Examples of classical umbral calculus. Umbral calculus is a symbolic method used by mathematicians in the 19th century to magi- cally “prove” some identities. For example, starting with the binomial formula: n n n k n k (39) (x + y) = x y − , k Xk=0 one gets: n n (40) (x + y)n = (x)k(y)n k, k − Xk=0 where Γ(x + n) (41) (x) = . n Γ(x) Another example involves the Bernoulli polynomials deﬁned by: n n k (42) Bn(x)= Bn kx . k − Xk=0 12 JIANZHOU

Here is a “proof” that d (43) Bn(x)= n Bn 1(x), dx · − n it involves the interchanging of B with Bn: n n d d n k d n n k k Bn(x) = Bn kx B − x dx dx k − ∼ dx k Xk=0 Xk=0 d n n 1 = (B + x) = n (B + x) − dx · n 1 − n 1 n 1 k k = n − B − − x · k Xk=0 n 1 − n 1 k n − Bn 1 kx ∼ · k − − Xk=0 = n Bn 1(x). · − Here is one more example [28, Proposition 4.2.1]. Suppose that a sequence bn n 0is related to an n 0 by the recursion relations: { } ≥ { } ≥ n n b = a n k k Xk=0 for all n, then raising the subscripts one gets:

n n bn = ak =(a + 1)n, k Xk=0 and so n n n n k n k a =(b 1) = ( 1) − b , − − k Xk=0 then, lowering the superscripts, one obtains:

n n k n a = ( 1) − b n − k k Xk=0 for all n. n k n k Now in (39), changing (x + y) to γn(x + y), x to γk(x), and y − to γn k(y), one gets (38). So umbral calculus is the right setting for studying− polynomial sequences of binomial type. MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 13

4.4. Hopf algebra in the modern classical umbral calculus. Since 1960s, Rota and his collaborators have constructed a rigorous mathematical foundation for umbral calculus based on the language of linear functionals, linear operators and their adjoint operators. We will present a brief survey in this and the next Subsections. For more details, see Roman [34]. Note the change from xn to a for n 0 deﬁnes a linear map n ≥ (44) L : C[x] C, xn a . → → n Denote L(p(x)) by L p(x) . Let D = d , then one has h | i dx m n n m (45) D x =(n)mx − , also consider the evaluation map: (46) ev : C[x] C, p(x) p(0). → 7→ Then one has (47) Dk xn = k!δ , h | i k,n and so (48) L p(x) = ev(f (D)(p(x))), h | i L where ∞ a (49) f (t)= n tn L n! n=0 X is a formal power series. So one gets a linear isomorphism:

(50) (C[x])∗ = C[[t]], L f (t), ∼ 7→ L where (C[x])∗ denotes the space of linear functionals on C[x]. The space of formal power series is an algebra with multiplication given by the binomial convolution of the coeﬃcients:

∞ a ∞ b ∞ tn n (51) k tk l tl = a b . k! · l! n! k k l n=0 Xk=0 Xl=0 X kX+l=n It is called the umbral algebra. For f(t) C[[t]], f(D): C[x] C[x] is an operator on C[x]. It can be characterized∈ as follows. For a→ R, let Ea : C[x] C[x] be the translation operator deﬁned by: ∈ → p(x) p(x + a). 7→ An operator P on C[x] is said to be translation invariant if P Ea = EaP for all a R. It is clear that f(D) is translation invariant. Conversely, ∈ 14 JIANZHOU any translation-invariant operator is of this form. Therefore, the umbral algebra is isomorphic to the algebra of translation invariant linear operators on C[x]. Because both C[x] and its dual, C[[t]], have structures of commutative algebras, one gets a bialgebra structure on C[x]. It is called the binomial bialgebra in one variable [26]. Note by (47), for f(t) C[[t]], one has ∈

∞ 1 (52) f(D)= f(D) xk Dk, k!h | i Xk=0 it follows that

∞ 1 n f(D)g(D) = f(D) xk g(D) xl Dn, n! k h | ih | i n=0 X kX+l=n and therefore, n f(D)g(D) xn = f(D) xk g(D) xl . h | i k h | ih | i kX+l=n This means the comultiplication on C[x] is given by:

n n n k n k (53) ∆x = x x − . k ⊗ Xk=0 Deﬁne an isomorphism ϕ : C[x] C[y] C[x, y] by ⊗ → ϕ(xm yn)= xmyn, ⊗ then we have: (54) ϕ(∆(xn))=(x + y)n, and in general, (55) ϕ(∆(p(x))) = p(x + y). The counit is 1 C[[t]]: ∈ (56) 1 p(x) = p(0) = ev(p(x)). h | i In fact, one actually gets a structure of a Hopf algebra on C[x]. The antipode is given by: (57) S(p(x)) = p( x). − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 15

4.5. Polynomial sequences of binomial type, coalgebra isomorphisms, and umbral operators. A coalgebra isomorphism U is a one-to-one onto linear operator on C[x] such that

n n n k n k (58) ∆U(x )= U(x ) U(x − ). k ⊗ Xk=0 n Write U(x )= pn(x). Such an operator is called an umbral operator by Mullin and Rota. Applying the isomorphism ϕ on both sides of (58), one gets:

n n (59) pn(x + y)= pk(x)pn k(y). k − Xk=0 n I.e., the sequence pn(x)= U(x ) is a polynomial sequence of binomial type for an umbral operator U. Converse, given a polynomial sequence n pn(x) of binomial type, setting U(x ) = pn(x) deﬁnes an umbral {operator} (i.e. a coalgebra isomorphism) on C[x].

4.6. Delta series and conjugate sequences. To study the umbral operators on C[x], it is more natural to study in the dual picture, i.e., to study isomorphisms of the dual algebra C[[t]]. Given a linear operator λ : C[x] C[x], its adjoint operator λ∗ : C[[t]] C[[t]] is deﬁned by: → →

(60) λ∗f(D) p(x) = f(D) λp(x) h | i h | i for all p(x) C[x]. The adjoint of an umbral operator is an automorphism of∈ the umbral algebra, and conversely, with due respect to suitable topology on C[[t]]. It follows that U ∗ is determined by U ∗(t) which is of the form

∞ an n (61) U ∗(t)= F (t)= t n! n=1 X with a = 0. Note for p(x) C[x], 1 6 ∈ 1 (62) p(x)= Dk p(x) xk, k!h | i k 0 X≥ 16 JIANZHOU

Therefore,

1 1 p (x) = Dk p (x) xk = Dk U(xn) xk n k!h | n i k!h | i k 0 k 0 X≥ X≥ 1 k n k 1 k n k = U ∗(D ) x x = (U ∗(D)) x x k!h | i k!h | i k 0 k 0 X≥ X≥ 1 k n k = (U ∗(D)) x x k!h | i k 0 X≥ 1 = F (D)k xn xk. k!h | i k 0 X≥ a n ∞ n Given a series of the form F (t)= n=1 n! t with a1 = 0, the sequence of polynomials 6 P 1 (63) p (x)= F (D)k xn xk n k!h | i k 0 X≥ is called the conjugate sequence of the delta series F (t).

Proposition 4.1. The exponential generating series of the conjugate sequence pn(x) n 0 of the delta series F (t) is { } ≥

∞ Xn (64) p (x) = exF (X). n n! n=0 X b k ∞ k Proof. We ﬁrst show that for any series g(t)= k=0 k! t , we have P (65) g(D) eXx = g(X). h | i This is because we have

b b g(D) xn = k Dk xn = k k!δ = b , h | i k!h | i k! · k,n n k 0 k 0 X≥ X≥ and so

Xn Xn g(D) eXx = g(D) xn = b = g(X). h | i h | i· n! n n! n 0 n 0 X≥ X≥ MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 17

Therefore,

∞ Xn ∞ Xn 1 p (x) = F (D)k xn xk n n! n! k!h | i n=0 n=0 k 0 X X X≥ 1 ∞ Xn = F (D)k xn xk k!h | n! i k 0 n=0 X≥ X 1 = F (D)k eXx xk k!h | i k 0 X≥ 1 = F (X)kxk = exF (X). k! k 0 X≥

4.7. Umbral compositions. Let us return to the point of view that the adjoints of umbral operators are the elements in the automorphism group of the umbral algebra C[[t]]. Given two umbral operators U and V , consider the composition U ∗ V ∗ of their adjoints U ∗ and V ∗. One has ◦

(66) U ∗ V ∗(f(D)) p(x) = f(D) V (U(p(x))) h ◦ | i h | i for all f(t) C[[t]] and p(x) C[x]. Suppose that ∈ ∈ n n k n (67) U(x )= pn(x)= an,kx , V (x )= qn(x), Xk=0 an n bn n (68) U ∗(t)= F (t)= t , V ∗(t)= G(t)= t , n! n! n 1 n 1 X≥ X≥ i.e., pn(x) and qn(x) are the conjugate sequences of f(t) and g(t) respectively.{ } Then{ } n n n k (69) rn(x)= V (U(x )) = V (pn(x)) = an,kV (x )= an,kqk(x) Xk=0 Xk=0 is the conjugate sequence of

an n an n U ∗(V ∗(t)) = U ∗(G(t)) = U ∗( t )= (U ∗(t)) n! n! n 1 n 1 a X≥ X≥ = n F (t)n = G(F (t)). n! n 1 X≥ The sequence rn(x) is called the umbral composition of the sequences p (x) and {q (x) }. { n } { n } 18 JIANZHOU

4.8. Group structures on in terms of umbral calculus. Let S γn and ϕn be two polynomial sequences of binomial type, deﬁned {by two} elements{ } G(X), F (X) , respectively. In other words, ∈ S ∞ Xn γ (x) = exp(xG(X)), n n! n=0 X ∞ Xn ϕ (x) = exp(xF (X)). n n! n=0 X Then the sequence θn obtained by the “umbral” substitution of ϕn into γn, i.e.,

(70) θn(x)= γn(x) xk ϕ (x) | 7→ k is a polynomial sequence of polynomial type deﬁned by F (G(X)). 5. Umbral Calculus, Quantum Mechanics, and Coherent States In this Section we relate umbral calculus to quantum mechanics, especially to Stone-von Neumann Theorem, Wick quantization, and coherent states. 5.1. Compositional inverse series and connection coeﬃcients. Suppose that pn(x) and qn(x) are the conjugate sequences of F (t) and G(t) respectively.{ } Then{ } n

(71) qn(x)= cn,kpk(x) Xk=0 for some coefficients cn,k called the connection coefficients. To find these coefficients, rewrite the above equations as n n k (72) V (x )= cn,kU(x ), Xk=0 where U and V are the umbral operators adjoint to F (D) and G(D) 1 respectively. Now apply U − on both sides to get: n 1 n k (73) (U − V )(x )= c x . ◦ n,k Xk=0 This means r (x)= n c xk is the conjugate sequence of { n k=0 n,k } 1 1 (74) (U − V )∗(t)= V ∗((U − )∗(t)) = G(f(t)), ◦ P where f(t) is the compositional inverse series of F (t), i.e., (75) F (f(t)) = f(F (t)) = t. MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 19

5.2. Delta series and associated sequences. Now one can apply the inverse umbral map n 1 (76) x = U − (pn(x)) in many of the formulas above. By (47) we get: k 1 D U − (p (x)) = k!δ . h | n i n,k The left-hand side is equal to 1 k k (U − )∗(D ) p (x) = f(D) p (x) , h | n i h | n i so we get: (77) f(D)k p (x) = k!δ . h | n i k,n n Recall pn(x) is the sequence of deformations of x , in the dual picture,{ Dn }is deformed to f(D)n . In the literature{ } on umbral { } { } calculus, pn(x) is called the associated sequence of the delta series f(t). { } Now we prove the theorems in [34, Section 2.4] in a similar fashion. Proposition 5.1. (The Expansion Theorem [34, Theorem 2.4.1] ) For h(t) C[[t]], ∈ ∞ 1 (78) h(D)= h(D) p (x) f(D)k k!h | k i Xk=0 where f(t) is a delta series and p (x) is its associated sequence. { n } Proof. By (52) we have

∞ 1 1 1 k k h(D) = (U ∗)− (U ∗(h(D))=(U ∗)− U ∗(h(D)) x D k!h | i Xk=0 ∞ 1 k 1 k = h(D) U(x ) (U ∗)− (D ) k!h | i Xk=0 ∞ 1 = h(D) p (x) f(D)k. k!h | k i Xk=0

Proposition 5.2. (The Polynomial Expansion Theorem [34, Theorem 2.4.2]) Let f(t) be a delta series and let pn(x) be the associated sequence of f(t). Then for any p(x) C[x]{, one has} ∈ 1 (79) p(x)= f(D)k p(x) p (x). k!h | i k k 0 X≥ 20 JIANZHOU

Proof. Applying the umbral operator U on both sides of (62), we get 1 U(p(x)) = Dk p(x) U(xk) k!h | i k 0 X≥ 1 k 1 = D U − U(p(x)) p (x) k!h | i k k 0 X≥ 1 1 k = (U − )∗(D ) U(p(x)) p (x) k!h | i k k 0 X≥ 1 = f(D)k U(p(x)) p (x). k!h | i k k 0 X≥ By induction one gets: k (80) f(D) pn(x)=(n)k pn k(x). · − Proposition 5.3. ([34, Theorem 2.4.5]) A sequence pn(x) of polynomials is associated to a delta series f(t) if and only{ if } (81) 1 p (x) = δ , h | n i n,0 (82) f(D)pn(x)= npn 1(x). − Proof. Suppose that pn(x) is associated with f(t). Taking k = 0 in (77), one gets (81). To{ get (82),} one uses (77) again to get: k 1 k 1 f(D) − f(D)p (x) = ev(f(D) − f(D)p (x)) h | n i n = ev(f(D)kp (x)) = f(D)kp (x) = k!δ . n h n i k,n The left-hand side is also equal to 1 k 1 k 1 1 (U − ) (D − ) f(D)p (x) = D − U − (f(D)p (x)) , h ∗ | n i h | n i so we have k 1 1 (83) D − U − (f(D)p (x)) = k!δ . h | n i k,n From this one gets: 1 n 1 (84) U − (f(D)pn(x)) = nx − . It follows that n 1 (85) f(D)pn(x)= nU(x − )= npn 1(x). − Conversely, if (81) and (82) hold, then k k f(D) pn(x) = ev(f(D) pn(x)) = ev((n)kpn k(x)) h | i − = (n)k 1 pn k(x) =(n)k δn k,0 ·h | − i · − = n! δ . · n,k MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 21

So one recovers (77).

5.3. Recursion formula for associated sequences and Heisen- berg commutation relation in umbral calculus.

Proposition 5.4. ( [34, Theorem 2.4.6]) If pn(x) is associated to f(t), then for any h(t) C[[t]], { } ∈ n n (86) h(D)pn(x)= h(D) pk(x) pn k(x). k h | i − Xk=0 Proof. This follows from (78) and (80) as follows:

∞ 1 h(D)p (x) = h(D) p (x) f(D)kp (x) n k!h | k i n Xk=0 ∞ 1 = h(D) pk(x) (n)kpn k(x) k!h | i − Xk=0 n n = h(D) pk(x) pn k(x). k h | i − Xk=0

By taking h(t)= t, one gets the following corollary:

Proposition 5.5. ( [34, Theorem 2.4.9]) If pn(x) is associated with f(t), then { } n n (87) pn′ (x)= D pn k(x) pk(x). k h | − i Xk=0 Proposition 5.6. ([34, Theorem 2.4.8]) If pn(x) is associated to f(t), then n n (88) xpn(x)= f ′(D) pn k(x) pk+1(x). k h | − i Xk=0 b k ∞ k Proof. Write f(t)= k=1 k! t . Then we have: P ∞ b f(D)(xp (x)) = k Dk(xp (x)) n k! n Xk=1 ∞ ∞ bk k bk k 1 = x D p (x)+ kD − p (x) k! n k! · n Xk=1 Xk=1 = xf(D)pn(x)+ f ′(D)pn(x). 22 JIANZHOU

By induction we get: k k k 1 f(D) (xpn(x)) = xf(D) pn(x)+ kf ′(D)f(D) − pn(x) k = xf(D) pn(x)+ k(n)k 1f ′(D)pn k+1(x). − − After evaluating at x = 0, one gets: k f(D) xpn(x) = k(n)k 1 f ′(D) pn k+1(x) . h | i − h | − i By the Polynomial Expansion Theorem,

∞ 1 xp (x) = f(D)k xp (x) p (x) n k!h | n i k Xk=0 ∞ 1 = k(n)k 1 f ′(D) pn k+1(x) pk(x) k! · − h | − i Xk=0 ∞ n = f ′(D) pn k(x) pk(x). k h | − i Xk=0

The above two Theorems concern the following two operators on C[x]: The operator D deﬁned by Dp(x) = p′(x), and the operator X deﬁned by Xp(x) = x p(x). They satisfy the Heisenberg canonical commutation relation: · (89) [D,X]=1. By taking adjoints, we get

(90) [X∗,D∗]=1. These adjoint operators on C[[t]] are given as follows:

(91) D∗g(t)= t g(t), X∗g(t)= g′(t). · 5.4. Umbral shift and more Heisenberg commutation relations in umbral calculus. Since we have pn(0) = 1 pn(x) = 0 for n 1, so one can write h | i ≥ k (92) pn(X)= cn,kx , n =1, 2,.... k 1 X≥ Then by (87),

n 1 − (93) pn′ (x)= cn k,1pk(x). − Xk=0 MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 23

So after integration one gets a recursion relation:

n 1 − x (94) pn(x)= cn k,1 pk(x)dx. − 0 Xk=0 Z

Therefore, the sequence p (x) of polynomials of binomial type is { n } completely determined by the sequence cn,1 . This recovers [36, p. 697, Corollary 1]. { } By (88) we also get a recursion formula:

n 1 1 − n (95) pn+1(x)= xpn(x) f ′(D) pn k(x) pk+1(x) . f ′(0) − k h | − i Xk=0

Recall the generating series of p (x) is { n }

∞ p (x) (96) n tn = exF (t), n! n=0 X

a n ∂ ∞ n where F (t) = n=1 n! t . Taking ∂x and then taking x = 0 on both sides, we get: P

c ∞ a (97) n,1 tn = F (t)= n tn. n! n! n=1 n=1 X X

It follows that

(98) cn,1 = an.

∂ This recovers [36, p. 697, Corollary 2]. In (96), change t to Dy = ∂y , let both sides act on yn and then take y = 0. This yields:

(99) p (x)= exF (Dy) yn . n h | i 24 JIANZHOU

From this one can derive a recursion formula as follows:

xF (Dy) n 1 xF (Dy) n 1 p (x) = e y y − = e Y (y − ) n h | · i h | i

1 = x pn 1(x + y) · f (D ) − ′ x y=0

1 = x pn 1(x). · f ′(Dx) − In the above we have used the following facts: 1 (U − )∗D = f(D ), F (f(t)) = t, F ′(f(t)) f ′(t)=1. y y · So we have recover the following identity in [34, Corollary 3.6.6]: 1 (100) pn(x)= x pn 1(x). · f ′(Dx) − The operator θ on C[x] deﬁned by:

(101) θf pn(x)= pn+1(x) is called the umbral shift. By (100), 1 (102) θf = X . ◦ f ′(D) By taking adjoint: 1 ∂ ∂ (103) θf∗ = = . f ′(t) ∂t ∂f Note:

(104) f(D)∗g(t)= f(t) g(t), · therefore, it is clear that

(105) [θf∗, f(D)∗]=1, MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 25 and so after taking adjoint,

(106) [f(D), θf ]=1. This can be checked directly from (82) and (101). 5.5. Coherent states in umbral calculus. Let us now interpret

∞ p (x) (107) n tn = exF (t) n! n=0 X from the perspective of quantum mechanics. First, inspired by [36, §9], we introduce a Hermitian product on C[x] a n ∞ n as follows. Given a delta series f(t) = n=1 n! t , with associated series pn(x) , let U : C[x] C[x] be the umbral operator deﬁned by: { } → P n (108) U(x )= pn(x). Deﬁne 1 (109) (p(x), q(x)) := [(U − p)(f(D ))q(x) . x¯ |x=0 It is easy to see that

(110) (pn(x),pm(x)) = δm,nn!. Furthermore,

(111) (f(Dx)p(x), q(x))=(p(x), θf q(x)). Also note

(112) f(Dx)1=0, so one can take 1 to be the vacuum vector 0 , a = f(D ) to be the | i x annihilator, and a† = θf to be the creator. Note we have

(113) a 0 =0, [a, a†] =1 | i and

∞ ∞ xF (t) pn(x) n npn 1(x) n xF (t) (114) ae = f(D ) t = − t = t e . x n! n! · n=0 n=0 X X This implies that the family exF (t) is a family of coherent states indexed by t. 5.6. A brief summary of coherent states. For the convenience of the reader who does not have a background in quantum mechanics and coherent states, we summarize some basics in this Subsection. We refer to Berezin-Shubin [4], Combescure-Robert [11], Hall [20] and Takhtajan [38] for proofs and more details. 26 JIANZHOU

5.6.1. Gaussian state as minimum uncertainty state. In Schr¨odinger representation (position representation), the state of a quantum system is described by a function Ψ(x) in x, such that

Ψ 2 = Ψ(x) Ψ(x) = Ψ(x) 2dx =1. || || h | i R | | Z The position and momentum operators are given by the self-adjoint operators: ∂ Q = x , P = i~ , · − ∂x respectively. They satisfy the canonical commutation relation: (115) [Q, P ]= i~. The expectation values of the position and the momentum are defined by x = Ψ(x) Q ψ(x) , p = Ψ(x) P ψ(x) , 0 h | | i 0 h | | i respectively. The variances of the position and the momentum are defined by (∆Q)2 = Ψ(x) (Q x )2 ψ(x) , (∆P )2 = Ψ(x) (P p )2 ψ(x) , h i h | − 0 | i h i h | − 0 | i The uncertainty relation is: ~2 (116) (∆Q)2 (∆P )2 . h i·h i ≥ 4 A state is called a minimum uncertainty state if the equality holds in the above inequality. The Gaussian state defined by 2 1/4 x (117) ϕ (x)=(π~)− exp 0 −2~ is a normalized state with x0 = p0 = 0. For this state we have ~ 2 1 2 x2/~ ϕ0(x) Q ϕ0(x) = x e− dx = , h | | i √π~ R 2 Z 2 ~ 2 1 x2/(2~) 2 d x/(2~) − ~ − ϕ0(x) P ϕ0(x) = e ( ) 2 e dx = , h | | i √π~ R − dx 2 Z and so we have ~2 (118) ϕ (x) (∆Q)2 ϕ (x) ϕ (x) (∆P )2 ϕ (x) = . h 0 | | 0 i·h 0 | | 0 i 4

This means that the Gaussian ϕ0(x) is a minimum uncertainty state. MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 27

5.6.2. Gaussian state as the ground state in the system of quantum harmonic oscillator. The Hamiltonian function of the harmonic oscillator is p2 x2 (119) H = + . 2 2 After the canonical quantization ∂ (120) x xˆ = Q = x , p pˆ = P = i~ , → · → − ∂x one gets the Hamiltonian operator for the quantum harmonic oscillator:

1 1 ∂2 (121) Hˆ = (P 2 + Q2)= ( ~2 + x2). 2 2 − ∂x2

Deﬁne the annihilator a and the creator a† by 1 1 (122) a = (Q + iP ), a† = (Q iP ), √2~ √2~ − then a and a† are adjoint to each other, and ~ (123) Hˆ = ~Nˆ + , 2 where

(124) Nˆ = a†a is the number operator. For any normalized Ψ(x), one has ~ ~ Ψ(x) Hˆ Ψ(x) = ~ aΨ(x) 2 + , h | | i || || 2 ≥ 2 with equality iﬀ

(125) aΨ(x)=0, or equivalently, ∂ (126) ~ Ψ(x)= x Ψ(x). ∂x · It is then clear that the minimum is achieved when Ψ(x) is equal to iθ ϕ0(x) up to a phase, i.e. multiplication by a factor of the form e . 28 JIANZHOU

5.6.3. Coherent state as the generating function of eigenstates in the system of quantum harmonic oscillator. Other eigenstates of quantum harmonic oscillator can be obtained by applying the creator repeatedly:

n (127) ϕn(x)=(a†) ϕ0. It is easy to see that

(128) aϕn(x)= nϕn 1(x), − (129) (ϕm(x),ϕn(x)) = m!δm,n. Physicists prefer to normalize these eigenstates and set:

1 1 n (130) n = ϕn(x)= (a†) ϕ0(x). | i √n! √n! The exponential generating series of the eigenstates is

n ∞ z † (131) ϕ (x) = eza ϕ (x). n n! 0 n=0 X Note

m n † † ∞ z¯ z eza ϕ (x) eza ϕ (x) = ϕ (x) ϕ (x) h 0 | 0 i m!n! h m | n i m,n=0 X ∞ m n z¯ z z 2 = δ = e| | . m!n! m,n m,n=0 X So physicists deﬁne the normalized coherent states by:

∞ n z 2/2 z z 2/2 za† (132) z := e−| | ϕ (x) = e−| | e ϕ (x). | i n n! 0 n=0 X It satisﬁes: (133) a z = z z . | i ·| i 5.6.4. Coherent states as minimum uncertainty states. To understand the properties of z one can do the following computations: | i ~ ~ z Q z = z a + a† z = (z +¯z), h | | i r 2h | | i r 2 ~ ~ z P z = i z a a† z = i (z z¯). h | | i − r 2h | − | i − r 2 − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 29

Similarly, ~ ~ 2 2 2 2 z Q z = z (a + a†) z = z a + a†a + aa† +(a†) z h | | i 2h | | i 2h | | i ~ ~ 2 2 2 = z a +2a†a +(a†) +1 z = [(z +¯z) + 1], 2h | | i 2 ~ ~ 2 2 2 z P z = z (a a†) z = [(z +¯z) 1]. h | | i − 2 h | − | i − 2 − It follows that ~ (∆Q)2 = z Q2 z z Q z 2 = , h i h | | i−h | | i 2 ~ (∆P )2 = z P 2 z z P z 2 = , h i h | | i−h | | i 2 and so ~2 z (∆Q)2 z z (∆P )2 z = . h | | i·h | | i 4 This means that z is a family of minimum uncertainty states, with expectation values| ofi the position and the momentum given by ~ ~ (134) x0 = (z +¯z), p0 = i (z z¯), r 2 − r 2 − respectively. 5.6.5. Coherent states as Gaussian states acted by Weyl-Heisenberg translations. Let us now derive an explicit expression for z . First, | i recall that ϕn(x) can be expressed in terms of Hermite polynomials: (135) 2 2 x x √2x x 1/4 n/2 ~ 1/4 ~ ϕn(x)=(π~)− 2− Hn( )e− 2 =(π~)− Hen( )e− 2 , √h · √~ where Hn(x) and Hen(x) are the physicists’ and probabilists’ Hermite polynomials respectively: n 2 n 2 n x2 d x2 n x d x H (x)=( 1) e e− , He (x)=( 1) e 2 e− 2 . n − dxn n − dxn Then one can use the formula for the exponent generating series of Hermite polynomials:

∞ n ∞ n t 2xt t2 t xt t2/2 (136) H (x) = e − , He (x) = e − , n n! n n! n=0 n=0 X X together with (135) and (132). The result is : 2 1/4 z 2/2 1 2 √2 z (137) z =(π~)− e−| | exp x + xz . | i · −2~ √~ · − 2 30 JIANZHOU

By (134), we have x + ip (138) z = 0 0 . √2~ Plugging this into (137), we get: 2 1/4 i i (x x0) (139) z =(π~)− exp x p exp p x exp − . | i −2~ 0 0 ~ 0 − 2~ If one introduces the Weyl-Heisenberg translation operator: i i i x0p0 p0 Q x0 P (140) Tˆ(z)= e− 2~ e ~ · e− ~ · , then one has (141) z = Tˆ(z)ϕ . | i 0 Using the Heisenberg commutation relation (142) [Q, P ]= i~, one can show that i (143) Tˆ(z) = exp (p Q x P ) . ~ 0 − 0 Another way to obtain (141) and hence also (139) is to rewrite Weyl- Heisenberg translations in terms of creator and annihilator: α 2 (144) Tˆ(z) = exp α a† α¯ a = exp | | exp α a† exp α¯ a . · − · − 2 · · − · Then because

(145) aϕ0 =0, one has α 2 (146) Tˆ(z)ϕ = exp | | exp α a† ϕ = z . 0 − 2 · 0 | i 5.6.6. Coherent states and representation of the Weyl-Heisenberg group. The following two formulas for the multiplication of the operators Tˆ(z) with Tˆ(z′) is well-known: i (147) Tˆ(z)Tˆ(z′) = exp σ(z, z′) Tˆ(z + z′), −2~ i (148) Tˆ(z)Tˆ(z′) = exp σ(z, z′) Tˆ(z′)Tˆ(z), −~ where for z = q + ip, z′ = q′ + ip′, σ(z, z′) is the symplectic product:

(149) σ(z, z′)= q p′ p q′. · − · MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 31

These multiplication formulas lead to the definition of the Heisenberg group H . It is the space of (t, z) R C with the group multiplication 1 ∈ × (150) (t, z)(t′, z′)=(t + t′ + σ(z, z′), z + z′). It is a Lie group of dimension 3. The Schrödinger representation is defined as the following irreducible unitary representation of H1 in L2(R):

it/2 Tˆ(z) (151) ρ(t, z)= e− e .

5.6.7. Stone-von Neumann Theorem. Suppose that Q˜ and P˜ are two self-adjoint operator on a Hilbert space H , such that the following commutation relation is satisﬁed : (152) [Q,˜ P˜]= i, then one can deﬁne a representation of the Heisenberg group H1 as follows:

it/2 i (153)ρ ˜(t, z)= e− exp (p Q˜ x P˜) , ~ 0 − 0 where x0 and p0 are related to z by (134). The Stone-von Neumann Theorem states that when the representationρ ˜ is irreducible, there is a unique unitary map U : H L2(R) such that → 1 (154) Uρ˜(t, z)U − = ρ(t, z).

5.6.8. Weyl quantization. The Lie algebra of the Heisenberg group H1 is the Heisenberg algebra h1, generated by e1, e2, e3, with the following commutation relations:

(155) [e1, e2]= e3, [e3, e1] = [e3, e2]=0.

n3 n2 n1 By the PBW Theorem, elements of the form e3 e2 e1 where each nk is a non-negative integer, span the universal enveloping algebra Uh1 and are linearly independent. The same statement also holds for elements n3 n1 n2 of the form e3 e1 e2 . Another choice for the basis of Uh1 is as follows: For two expressions X and Y , consider the expansion: n n n n k n k (156) (aX + bY ) = Ck,n k(X,Y )a b − , k − Xk=0 then elements of the form

n1!n2! n3 n1+n2 e3 Cn1,n2 (e1, e2) (n1 + n2)! form a basis of Uh1. 32 JIANZHOU

The three diﬀerent bases of Uh1 corresponds to three diﬀerent quan- tizations: the qp-symbol, the pq-symbol and the Weyl symbol. Given an operator of the form: m n (157) f(q,p)= cm,nq p , m,n 0 X≥ its qp-quantization, pq-quantization, and Weyl quantization are given by m n (158) fˆqp = cm,nQ P , m,n 0 X≥ n m (159) fˆpq = cm,nP Q , m,n 0 X≥ m!n! (160) fˆ = c Cm+n(Q, P ) Weyl m,n (m + n)! m,n m,n 0 X≥ respectively. Conversely, f(q,p) is called the qp-symbol of fˆqp, the pq- symbol of fˆpq, and the Weyl symbol of fˆWeyl. n n ~ ∂ n The transformation from p to P = ( i ∂x ) can be achieved by Fourier transform and inverse Fourier transform:−

1 ∞ ixp (161) f˜(p)= f(x) e− dx, √2π Z−∞ 1 ∞ (162) f(x)= f˜(p) eixp dp. √2π Z−∞ First note: ∂ 1 ∞ ( i~ )nψ(x) = (~p)nψ˜(p) eixp dp − ∂x √2π Z−∞ 1 ∞ ∞ n i(x y)p/~ = p e − ψ(y)dydp, 2π~ Z−∞ Z−∞ so we have:

1 ∞ ∞ i(x y)p/~ (163) fˆ ψ(x)= f(x, p)e − ψ(y)dydp. qp 2π~ Z−∞ Z−∞ If we write ∞ (164) fˆqpψ(x)= K(x, y)ψ(y)dy, Z−∞ then the kernel function K(x, y) can be obtained from the symbol f(q,p) by Fourier transform:

1 ∞ ip(y x)/~ (165) K(x, y)= f(x, p)e− − dp. 2π~ Z−∞ MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 33

Conversely, the symbol f(q,p) can be obtained from the kernel function by inverse Fourier transform:

∞ ip(y q)/~ (166) f(q,p)= K(q,y)e − dy. Z−∞ Given two qp-symbols f1 and f2, consider the composition operator ˆ ˆ ˆ ˆ ˆ f1f2. Denote by K, K1, K2 the kernels of the operators f1,qpf2,qp, f1, fˆ2. Then ∞ (167) K(x, y)= K1(x, z)K2(z,y)dz. Z−∞ Using (165) and (166), one can see that fˆ1,qpfˆ2,qp is fˆqp, where 1 i f(q,p)= 2 exp (p(y q) p1(z q) p2(y z)) (2π~) R4 ~ − − − − − Z f (q,p )f (z,p )dp dp dydz. · 1 1 2 2 1 2 Using the formula

1 ∞ i(x q)y/h (168) e − dy = δ(x q), 2π~ − Z−∞ one can perform the integration in y to get: 1 i f(q,p)= exp (p1 p)(q1 q) f1(q,p1)f2(q1,p)dp1dq1. 2π~ R2 −~ − − · Z Now note:

1 i k l exp (p1 p)(q1 q) p1q1dp1dq1 2π~ R2 −~ − − · Z 1 i i l ~ p(q1 q) ~ p1(q1 q) k = q1e − e− − p1dp1 dq1 2π~ R R · Z Z k 1 i i l ~ p(q1 q) (k) = q1e − 2π δ (q1 q)dq1 2π R · ~ − Z k 1 ∂ i k l ~ p(q1 q) ~ k = ( 1) k q1e − 2π(i ) δ(q1 q)dq1 2π R − ∂q · − Z 1 k 1 i i k l j k j ~ p(q1 q) ~ k = ( 1) (l)jq1− ( p) − e − 2π(i ) δ(q1 q)dq1 2π R − ~ · − j=0 Z X k k l j i k j = ( i~) (l) q − ( p) − − j ~ j=0 X ~ ∂ l = (p i )q1 q1=q. − ∂q1 | 34 JIANZHOU

So one gets: ∂ (169) f(q,p)= f (q,p i~ )f (q ,p) . 1 − ∂q 2 1 1 q1=q

From this one gets: n ∞ 1 ~ (170) f(q,p) = ∂nf (q,p) ∂nf (q,p) n! i p 1 · q 2 n=0 X ∂2 (171) = exp i~ [f (q,p )f (q ,p)] . − ∂p ∂q 1 1 2 1 1 1 q1=q,p1=p

Note 1 irQ isP m irq n isp 2 e e ψ(x) q e− dq p e− dp drds (2π) R2 · · Z Z Z 1 irQ isP m (m) n (n) = 2 e e ψ(x)2πi δ (r) 2πi δ (s)drds (2π) R2 · Z 1 irx ~ m (m) n (n) = 2 e ψ(x + s )2πi δ (r) 2πi δ (s)drds (2π) R2 · Z m n m+n ∂ ∂ irx ~ = ( i) m n (e ψ(x + s ))δ(r)δ(s)drds − R2 ∂r ∂s Z ∂ = xm ( ~ )nψ(x) · − ∂x \m n = ((q p )qpψ)(x). So one gets:

1 ir(Q q) is(P p) ˆ − − (172) (fqpψ)(x)= 2 f(p, q)e e ψ(x)dpdqdrds. (2π) R4 Z The Weyl quantization is given by

ˆ 1 ir(Q q)+is(P p) (173) (fWeylψ)(x)= 2 f(p, q)e − − ψ(x)dpdqdrds. (2π) R4 Z By (140) and (143), (174) [ei(rQ+sP )ψ](x)= ei(rx+~rs/2)ψ(x + ~s), whence ˆ 1 ~ (fWeylψ)(x)= 2 exp i(rx + rs/2) (175) (2π) R4 Z i(rq+sp) e− f(q,p)ψ(x + ~s)dpdqdrds. · MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 35

1 By changing s to ~ (y x), integrating over r, then integrating over q, one gets: −

1 i x + y ˆ ~ (x y)p (176) (fWeylψ)(x)= e − f( ,p)ψ(y)dydp. 2π~ R2 2 Z The kernel of this operator is

1 i (x y)p x + y (177) K(x, y)= e ~ − f( ,p)dp. 2π~ R 2 Z Note K(x,y) is obtained as the Fourier transform p x y of the x+y 7→ − function f( 2 ,p). Using the inversion formula, one gets:

(178) f(q,p)= K(q ξ/2, q + ξ/2)eipξ/~dξ. R − Z Using the above two formulas, one can derive the composition formula for Weyl symbols. The result is:

(179) f(q,p)=(f f )(q,p)= ei~L[f (q ,p )f (q ,p )] , 1∗ 2 1 1 1 2 2 2 q1=q2=q,p1=p2=p where

1 ∂2 ∂2 (180) L = . 2 ∂q ∂p − ∂q ∂p 1 2 2 1 This leads to deformation quantization. By comparing this summary with the formalism of the umbral calculus, one can see that the umbral calculus can be embedded in the for- mulism of quantum mechanics. More precisely, one can regard the for- a n a n ∞ n ∞ n mal power series n=0 n! t as the symbol of the operator n=0 n! D , and so the latter is the quantization of the former by any of the above three quantizationP schemes. In the next subsection, we willP see that it is more natural to interpret this quantization as the Wick quantization. 5.7. Segal-Bargmann transform, Wick quantization and umbral calculus. In Section 5.5 we have interpreted the exponential generating series of polynomial sequences of binomial type as coherent states of some quantization. To elaborate this interpretation further, in this Subsection we will relate umbral calculus to further results about coherent states in quantum mechanics. 5.7.1. Overecompleteness of coherent states. Given ψ(x) L2(R), one has a decomposition: ∈

∞ (181) ψ(x)= anϕn(x), n=0 X 36 JIANZHOU where the coeﬃcients an are given by: 1 (182) a = (ψ(x),ϕ (x)). n n! n Because z 2/2 n (183) (ϕ (x),ϕ (x)) = z ϕ (x) = e−| | z¯ , n z h | n i one gets:

∞ z 2/2 n (184) (ψ(x),ϕ (x)) = z ϕ(x) = e−| | a z¯ . z h | i n n=0 X Furthermore, m n n dz dz¯ z 2/2 z¯ z 2/2 z dz dz¯ z ϕm(x) z ∧ = e−| | e−| | ϕn(x) ∧ Ch | i| i 2i C m! · n! 2i n=0 Z − Z X − = ϕm(x), and so dz dz¯ (185) z ψ(x) z ∧ = ψ(x). Ch | i| i 2i Z − Indeed, dz dz¯ ∞ dz dz¯ z ψ(x) z ∧ = z anϕn(x) z ∧ Ch | i| i 2i Ch | i| i 2i n=0 Z − Z X − ∞ dz dz¯ = an z ϕn(x) z ∧ Ch | i| i 2i n=0 X Z − ∞ = anϕn(x)= ψ(x). n=0 X 5.7.2. Segal-Bargmann space and Segal-Bargmann transform. The Fourier- Bargmann transform is deﬁned by the following formula: B ♯ 1/2 (186) F v(z) =: v (z)=(2π~)− u , v . u h z i It is an isometry from L2(R) into L2(R2): the scalar product of two 2 ♯ ♯ states ψ, ψ′ L (R) can be expressed in terms of ψ , (ψ′) : ∈ ♯ ♯ dz dz¯ (187) ψ′, ψ = (ψ′) (z)ψ (z) ∧ . h i C · 2i Z − Using the Fourier-Bargmann transform, the overcompleteness of coherent states (185) can be written as:

♯ dz dz¯ (188) ψ(x)= ψ (z)ϕz(x) ∧ . C 2i Z − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 37

Let x ip (189) ξ = 0 − 0 . √2 The Segal-Bargmann space is the space of holomorphic functions F in ξ such that 1/2 1 2 ξ 2/~ dξ dξ (190) F ~ := F (ξ) e−| | ∧ < + . || || π C | | 2i ∞ Z − This is a Hilbert space in which the space of polynomials in ξ is dense. ♯ ♯ z 2/2 The transformation ψ ψHol(ξ) := ψ (z)e| | is called the Segal- Bargmann transform and7→ is denoted by . One can check that: B 1 (191) (xψ)(ξ)= (~∂ξ + ξ) ψ(ξ), B √2 B 1 (192) (~∂xψ)(ξ)= (~∂ξ ξ)Bψ(ξ), B √2 −

(193) [a†ψ](ξ)= ξ [ψ](ξ), B B (194) [aψ](ξ)= ∂ [ψ](ξ). B ξB The integral kernel of the Bargmann transform is the Bargmann kernel: B 2 2 3/4 1/2 1 x ξ (195) (x, ξ)=(π~)− 2− exp √2x ξ + , B −~ 2 − · 2 where x R, ξ C. The Bargmann kernel is a generating function for the Hermite∈ functions:∈ ξn (196) B(x, ξ)= φ (x). ((2π)nn!)1/2 n 1 X The integral kernel of − is B 2 ¯2 1 31/4 n/2 1 x ξ (197) − (ξ, x)=(π~)− 2− exp √2x ξ¯+ . B −~ 2 − · 2 5.7.3. Wick symbol, Wick quantization and Wick -product. For an operator A in the Wick normal form ∗ k l (198) AˆWick = Ak,l(a†) a , Xk,l deﬁne its Wick symbol by

k l (199) A(ξ, ξ¯)= Ak,lz z¯ . Xk,l 38 JIANZHOU

The operator AˆWickis called the Wick quantization of the Wick symbol A(ξ, ξ¯). The Wick symbol can be obtained by the following formula: (Aˆ ϕ ,ϕ ) (200) A(ξ, ξ¯)= Wick ξ ξ . (ϕξ,ϕξ)

Using the Segal-Bargmann transform, one can identify Wˆ Wick with an operator on the Segal-Bargmann space, also denoted by AˆWick:

1 η¯(η ξ) dη dη¯ (201) (AˆWickf)(ξ)= A(ξ, η¯)f(η)e− − ∧ . π C 2i Z − If A1(ξ, ξ¯) and A2(ξ, ξ¯) are the Wick symbols of operators Aˆ1,Wick and Aˆ2,Wick, then the Wick symbol of the operator AˆWick = Aˆ1,WickAˆ2,Wick is given by

1 (η ξ)(¯η ξ¯) dη dη¯ (202) A(ξ, ξ¯)= A1(ξ, η¯)A2(η, ξ¯)e− − − ∧ . π C 2i Z − Write A(ξ, ξ¯)= A1(ξ, ξ¯) Wick A2(ξ, ξ¯). The following formula holds: (203) ∗ ∂2 A (ξ, ξ¯) A (ξ, ξ¯) = exp ~ (A (ξ, η¯)A (η, ξ¯)) . 1 ∗Wick 2 ∂η∂η¯ 1 2 η=ξ,η¯=ξ¯

5.7.4. Umbral calculus from the point of view of quantum mechanics. Now we can interpret umbral calculus using Stone-von Neumann The- orem. In Section 5.5, given a delta series in the sense of [34], we have constructed an inner product on the space C[ξ] of polynomials in ξ, on which the operator f(Dξ) is interpreted as the annihilator a and its adjoint θf is interpreted as the creator. The umbral operator n U : C[ξ] C[ξ], ξ pn(ξ), can now be interpreted as the intertwin- ing operator→ from the7→ Segal-Bargmann space to the above space, as predicted by the Stone-von Neumann Theorem. As already mentioned in Section 5.5, the exponential generating series of the associated sequence of polynomials can be identiﬁed with the image of the coherent states in the Segal-Bargmann space under the map U. One can use the above interpretation to deﬁne an umbral -product as follows: ∗ 1 1 (204) g(x) h(x) := U(U − (g(x)) U − (h(x))). ∗ · More explicitly,

∞ ∞ ∞ (205) a p (x) b p (x) := a b p (x). k k ∗ l l k l k+l Xk=0 Xl=0 k,lX=0 For a special case of this product, see [?, Appendix]. MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 39

5.8. Segal-Bargmann transform from the point of view of umbral calculus. The Segal-Bargmann transform is essentially a linear map (206) C[x] C[x], He (x) xn. → n 7→ The inverse map is then given by (207) xn He (x). 7→ n From the point of view of umbral calculus, this is a Sheﬀer operator. We will recall in this Subsection the theory of Sheﬀer sequences in umbral calculus. Since it is very similar to the theory of associated sequences, we will be very brief and refer the interested reader to Roman [34] for details.

5.8.1. Sheffer sequences. The reference for this subsection is [34, Sec- tion 2.3]. Let f(t) be a delta series and let g(t) be an invertible series. Then there exists a unique sequence sn(x) n 0 of polynomials satis- fying { } ≥ (208) g(t)f(t)k s (x) = n!δ h | n i n,k for all n, k 0. This sequence is called the Sheffer sequence. for the pair (g(t), f≥(t)). Similar to the associated sequence, one has the Ex- pansion Theorem and the Polynomial Expansion Theorem for Sheffer sequences. There are several ways to characterize the Sheffer sequences: The following statements are equivalent to each other:

(a) sn(x) n 0 is a Sheffer sequence for the pair (g(t),g(t)). { } ≥ (b) sn(x) n 0 satisfies the Sheffer identity: { } ≥ n n (209) sn(x + y)= pk(y)sn k(x), k − Xk=0 where pn(x) n 0 is associated to f(t). { } ≥ (c) The exponential generating series of sn(x) n 0 is given by: { } ≥ ∞ s (x) 1 (210) n tn = exF (t), n! g(F (t)) n=0 X where F (t) is the compositional inverse of f(t). (d) For n 0, ≥ n 1 k n k (211) s (x)= f(F (t))− f(t) x x . n h | i Xk=0 (e) g(D)sn(x) n 0 is the associated sequence for f(t). { } ≥ 40 JIANZHOU

(f) For n 0, ≥ (212) f(D)sn(x)= nsn 1(x). − By comparing (f) with (81) and (82), one is led to the problem finding g(t) from the Sheffer sequence. This can be solved as follows. By (210) one has: ∞ s (x) F (t) n′ tn = exF (t), n! g(F (t)) n=0 X and so ∞ s (0) t (213) n′ f(t)n = . n! g(t) n=0 X 5.8.2. Sheffer operators and umbral compositions. The materials in this subsection is taken from [34, Section 3.5]. Let s (x) be the Sheffer { n } sequence for (g(t), f(t)). Then the linear operator λf,g on C[x] defined by n (214) λg,f x = sn(x) is called the Sheffer operator for sn(x) or for (g(t), f(t)). The adjoint operator of a Sheffer operator can{ be} characterized as follows. Let pn(x) be the sequence of polynomials associated to f(t), and denote by{ λ }: C[x] C[x] the umbral operator defined by λ xn = p (x). f → f n Then the Sheffer operator λg,f is related to λf by: 1 (215) λ = λ . g,f g(D) f For h(t) C[[t]], ∈ 1 1 1 (216) λg,f∗ h(t)= λf∗ h(t)= λf∗ ( h(t)) = h(F (t)). g(D)∗ g(t) g(F (t)) Conversely, an operator on C[[t]] of the form h(t) 1 h(F (t)) is a 7→ g(F (t)) Sheffer operator, where F (t) is a delta series inverse to f(t) and g(t) is an invertible series. From the above characterization of Sheffer operator one can deduce that the set of all Sheffer operators is a group under composition:

(217) λg(t),f(t) λh(t),l(t) = λg(t)h(f(t)),l(f(t)) , 1 ◦ (218) λg−(t),f(t) = λ1/g(F (t)),F (t). As corollaries one has the following results. A Sheffer operator maps Sheffer sequences to Sheffer sequences. If s x) and r (x) are Sheffer { ( } { n } sequences, then the linear operator defined by: sn(x) rn(x) is a Sheffer operator. 7→ MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 41

5.8.3. Sheffer shifts and canonical commutation relation. Let s (x) { n } be the Sheffer sequence for (g(t), f(t)). The linear operator θg,f on C[x] defined by

(219) θg,f sn(x)= sn+l(x) is called the Sheffer shift for s (x) or for (g(t), f(t)). { n } If sn(x) is Sheffer for (g(t), f(t)), then pn(x) = g(D)sn(x) is associated{ to} f(t). and so for n 0, { } ≥ 1 1 1 θg,f g(D)− pn(x)= g(D)− pn+1(x)= g(D)− θf pn(x). Thus, 1 (220) θg,f = g(D)− θf g(D). By a computation using the chain rule, one gets: g (D) 1 (221) θ = x ′ , g,f · − g(D) f (D) ′ and in particular, one gets the following recursion formula for the Shef- fer sequence: g (D) 1 (222) s (x)= x ′ s (x). n+1 · − g(D) f (D) n ′ It is clear that the following commutation relation is satisfied:

(223) [f(Dx), θg,f ]=1. 5.8.4. Exponential generating series of Sheﬀer sequence as coherent states. As in §5.5, we modify the scalar product introduced in [36, §9] to a Hermitian inner product. For p(x), q(x) C[x], deﬁne ∈ 1 (224) (p(x), q(x)) := ((θ− p)(f(D))g(D)q(x) . g,f |x=0 It is easy to see that

(225) (sn(x),sm(x))f,g = δm,nn!. Furthermore,

(226) (f(Dx)p(x), q(x))f,g =(p(x), θg,f q(x))f,g. Also note

(227) f(Dx)1=0, so one can take 1 to be the vacuum vector, a = f(Dx) to be the annihilator, and a† = θg,f to be the creator. Note we have

∞ ∞ xF (t) sn(x) n nsn 1(x) n xF (t) a g(t)e = f(D ) t = − t = t g(t)e . x n! n! · n=0 n=0 X X 42 JIANZHOU

This implies that the family g(t)exF (t) is a family of coherent states indexed by t. It is clear that all the discussion about the relationship between quantum mechanics and associated sequences can be generalized to the Sheﬀer sequences. We omit the details.

5.8.5. Segal-Bargmann transform as a Sheﬀer operator. Recall the exponential generating series of Hen(x) is

∞ xt t2/2 Hen(x) n (228) e − = t , n! n=0 X so by (c) above, Hen(x) n 0 is the Sheﬀer sequence for (f(t),g(t)) = t2/2 { } ≥ n n (t, e− ). Therefore the map He (x) x and its inverse x n 7→ 7→ Hen(x) are both Sheﬀer operators.

6. Interpolating Statistics, Deformed Exponential Functions, and Deformed Entropy Functions As we have discussed above, the associated sequences in umbral calculus can be regarded as deformed power functions, and hence their exponential generating series can be regarded as deformed exponential functions. This leads us to the generalized thermostatics [33] where another kind of deformed exponential functions arise. In this Section we will explain how the interpolating statistics are related to deformed exponential functions and deformed entropy functions. We will show that all interpolating statistics are deformed exponential functions in the sense of [33].

6.1. The φ-logarithm and φ-exponential functions. These were introduced by Naudt [32]. Fix a strictly positive non-decreasing function φ(u), deﬁned on the positive numbers (0, + ). The φ-logarithm function is deﬁned by: ∞ u 1 (229) ln (u)= dv , u> 0. φ φ(v) Z1 The inverse of the function lnφ(x) is called the φ-exponential function and is denoted expφ(x). A typical example used by Naudt [32] is the case of φ(u) = uq. In this case the deformed logarithm function is

u u1−q 1 1 1 q− , if q =1, (230) lnq(u)= dv = − 6 vq log(u), if q =1. Z1 ( MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 43

This was introduced in the context of nonextensive statistical physics in Tsallis [41].

6.2. All interpolating statistics are φT-exponential functions. For our purpose we will take: n φT(p)= p Tn 1p , − − n 2 X≥ and perform the following formal calculation: p du p du lnφT (p): = = n φ(u) u n 2 Tn 1u Z Z − ≥ − m e p ∞ ∞ du = T ui P i u m=0 i=1 Z XX p ∞ m1 + + mn mi n 1 = ··· Ti u − du m1,...,mn · n=0 P m i=n i Z X iXi Y ∞ n m1 + + mn mi p = ln(p)+ ··· Ti . m1,...,mn · n n=1 P m i=n i X iXi Y In other words, pn (231) ln (p) = ln(p)+ a (T) , φT n n n 1 X≥ where each e m + + m (232) a (T)= 1 ··· n T mi n m ,...,m i 1 n i PiXmii=n Y is a weighted homogeneous polynomial in T1,...,Tn of degree n, with deg Tj = j. For example,

a1 = T1, 2 a2 = T2 + T1 , 3 a3 = T3 + T2T1 + T1 , 2 2 4 a4 = T4 +2T3T1 + T2 +3T2T1 + T1 .

Therefore, (231) shows that lnφT is a formal deformation of the ln- function parameterized by T . { n} To ﬁnd an expression for thee formal inverse expφT , introduce the following series in p: ˜ g n (233) XφT (p) := exp(lnφ(p)) = p(1 + bn(T)p ), n 1 X≥ 44 JIANZHOU where each bn(T) is a weighted homogeneous polynomial in a1,...,an, hence it is a weighted homogeneous polynomial in T1,...,Tn:

mi ai (234) bn = m . i i mi! P m i=n i iXi Y For example,

b1 = a1 = T1, a a2 T b = 2 + 1 = 2 + T 2, 2 2 2! 2 1 a a a3 T 7 b = 3 + a 2 + 1 = 3 + T T + T 3. 3 3 1 2 6 3 6 2 1 1 Now applying the Lagrange inversion, one can express p as a formal power series in XΦT :

∞ T n (235) p = XφT + cn( )XφT , n=2 X for example, c = b = a = T , 2 − 1 − 1 − 1 1 3 1 c = b +2b2 = a + a2 = T + T 2, 2 − 2 1 −2 2 2 1 −2 2 1 1 8 1 4 c = b +5b b 5b3 = a +2a a a3 = T + T T T 3. 3 − 3 2 1 − 1 −3 3 2 1 − 3 1 −3 3 3 2 1 − 1 So we get:

∞ u nu (236) expφT (u)= e + cn(T)e . n=2 X It is clearly a deformationg of the ordinary exponential function parameterized by Tn . factor((1/24)*Chi([1,1,1,1],{ } [4])+(1/6)*Chi([1,1,1,1], [1, 3])+(1/4)*Chi([1,1,1,1], [1, 1, 2])+(1/8)*Chi([1,1,1,1], [2, 2])+(1/24)*Chi([1,1,1,1], [1, 1, 1, 1])) Now deﬁne a space Φ by:

n (237) Φ= φT(p)= p T p . { − n } n 2 X≥ Recall the space of interpolating statistics is deﬁned by

(238) = w(X)= X + w Xn . S { n } n 2 X≥ MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 45

n Given any formal series φT(p) = p n 2 Tnp Φ, deﬁne a formal − ≥ ∈ series wφT in as follows: S P n (239) wφT = X + cn(T)X . n 2 X≥ Then φ (p) w (X) deﬁnes a map g : Φ . T 7→ φT → S Proposition 6.1. The map g is a one-to-one map. n Proof. Given a formal power series w(X) = X + n 2 wnX in the space of interpolating series, one can apply Lagrange≥ inversion to get theS inverse series P n X(w)= w + wˆnw . n 2 X≥ Take bn =w ˆn 1. One can solve (234) to express an in terms of bn , and one can solve− (232) to express T in terms{ of} a . Therefore,{ } { n} { n} one can express Tn in terms of bn =w ˆn 1 , and hence also of wn . This gives us an{ inverse} map { Φ. More− explicitly,} { } S → n u + n 2 wˆnu X(u) 1 ≥ (240) φ(u)= n 1 = dX(u) = d . 1+ n 2 nwˆnu − du log X(u) P≥ du P

Proposition 6.2. The following formula for logφ(p) holds when φ is determined by (240): (241) log (p) = log X(u) log X(1). φ − Proof. p du p d X(u) log (p) = = du du φ φ(u) X(u) Z1 Z1 X(p) dX = = log X(u) log X(1). X − ZX(1)

6.3. The deduced logarithm and the φ-entropy. Under the con- dition that 1 v (242) c := dv < + , φ(v) ∞ Z0 Naudt [33, §10.3] introduces a function: 1/u 1 v − (243) χ(u)= dv , φ(v) Z0 46 JIANZHOU and he calls the deformed logarithm lnχ associated with χ the deduced logarithm. Naudts [33, §11.1] deﬁnes the φ-entropy by:

(244) Hφ(p)= pi lnχ(1/pi), i=1 X where p = p ,...,p , p 0, p = 1. After a short calculation: 1 n i ≥ i i n P pi 1 u v (245) H (p)= p dv du. φ − i u2 φ(v) i=1 1 0 X Z Z Proposition 6.3. The partial derivatives of φ-entropy are given by negative the φ-logarithms up to the constant c:

∂Hφ(p) (246) = logφ(pi) c, ∂pi − − where c is a constant deﬁned by: 1 v (247) c = dv. φ(v) Z0 Proof. We have the following computations: ∂H (p) pi 1 u v 1 pi u φ = dv du du ∂p − u2 φ(v) − p φ(u) i Z1 Z0 i Z0 pi u v 1 1 pi u = dv d du φ(v) u − p φ(u) Z1 Z0 i Z0 pi 1 1 v = du dv = ln (p ) c. − φ(u) − φ(v) φ i − Z0 Z0

6.4. All deformed entropy functions are φT-entropy functions. Let ξ(u) be deﬁned by u v (248) ξ (u)(u) := dv. φ φ(v) Z0 It is related to χ(u) as follows: 1 (249) ξ (u)(u)= . φ χ(1/u) MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 47

Denote by ξT(u) the series ξ(u) when φ is the series φT. We have u v T ξ (u) = n dv 0 v n 2 Tn 1v Z − ≥ − u ∞ P m1 + + mn mi n = ··· Ti v dv m1,...,mn · n=0 P m i=n i Z X iXi Y ∞ m + + m un+1 = u + 1 ··· n T mi . m ,...,m i · n +1 n=1 1 n i X PiXmii=n Y w n ∞ n Proposition 6.4. Suppose that F (X) = n=1 n X , so that X = X(w) is obtained by Lagrange inversion from w(X) = X dF (X) , and P dX φ(u) is determined from X(u) by 1 φ(u)= d , du log X(u) then we have

(250) ξφ(u)= F (X(u)). Proof. This is proved by a straightforward computation: u v u d log X(v) ξ (u)(u) = dv = v dv φ φ(v) · dv Z0 Z0 u w X(u) w(X) X(u) dF (X) = dX(w)= dX = dX X(w) X dX Z0 Z0 Z0 = F (X(u)).

We formally deﬁne p ξT(u) (251) H˜ (p) := du. φT u2 Z By a formal calculations, we have:

∞ m + + m pn ˜ 1 n mi H T (p) = p(ln(p)+ ··· T ). φ − m ,...,m i · n(n + 1) n=1 1 n i X PiXmii=n Y Deﬁne a space of deformed entropy by:

∞ pn (252) = Hs(p)= p(ln(p)+ s . H { − n n(n + 1))} n=1 X 48 JIANZHOU

n Then one can deﬁne a map f : Φ by sending φ(T)= n 2 Tn 1p →pn H − ≥ − to Hs(p)= p(ln(p)+ ∞ s ), where − n=1 n n(n+1) P

P m1 + + mn mi (253) sn = ··· Ti , n 1. m1,...,mn ≥ P m i=n i iXi Y Proposition 6.5. The map f : Φ is a one-to-one correspondence. → H Proof. For any deformation of the standard Boltzmann-Gibbs-Shannon pn entropy of the form Hs(p) = p(ln(p) + ∞ s ), one needs − n=1 n n(n+1) to ﬁnd a unique sequence of weighted homogeneous polynomial T = P n Tn(s1,...,sn) of degree n, where deg sn = n, such that Hs(p)= HφT (p). Indeed, one needs to solve the sequence (253) of equations. These are equivalent to e 1 (254) 1+ s pn = . n 1 T pn n 1 n 1 n X≥ − ≥ Its solution is clearly: P 1 (255) 1 T pn = . − n 1+ s pn n 1 n 1 n X≥ ≥ P 6.5. Principle of maximum entropy applied to φ-entropy. Let us formally extremize the modiﬁed φ-entropy functions n pi 1 v u (256) H˜ (p)= p du dv φ − i v2 φ(u) i=1 X Z Z under the constraints (257) p + + p =1, 1 ··· n (258) p E + + p E = E, 1 1 ··· n n by the method of Lagrange multiplier. The solution is given by ∂ H˜φ(p)= a1 + bEi, i =1, . . . , n, ∂pi for some constants a1 and b. By (246), (259) log (p )= a + bE , − φ i i (a+bEi) where a = a1 + c. We let Xi = e− . So the solution is given by (260) p = exp ( (a + bE )), i =1, . . . , n, i φ − i and equivalently,

(261) Xi = exp(lnφ(pi)). MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 49

6.6. Generalized statistics as critical points of φT-entropy. Now we have the following commutative diagrams of one-to-one correspon- dences: g (262) / Φ S ❅❅ ❅❅ ❅❅ f h=fg ❅ H As a corollary, we have: Every generalized statistics can be obtained by applying the Principle of Maximum Entropy to the φT-entropy function for a unique φT Φ, hence also to a unique deformed entropy Hs(p) . ∈ ∈ H 6.7. Generalized Boson-Fermion correspondence on the spaces Φ and . Recall we have deﬁned an involution σ : in §3.2. Now usingH the diagram (262), the involution σ inducesS an→ involutionS τ : Φ Φ and an involution ρ : , so that the following diagrams commute:→ H → H g / Φ ? > ⑦⑦S ⑥⑥ σ ⑦⑦ τ ⑥⑥ ⑦⑦ ⑥⑥ f ⑦⑦ ⑥⑥ / Φ g > S ⑥⑥H f ⑥⑥ ⑥⑥ρ ⑥⑥ H We call these involutions the generalized Boson-Fermion correspondence on the spaces Φ and respectively. H 6.8. Entropy as φ-entropy. In §2.6 we have deﬁned the entropy of the one-particle partition function as the negative of the Legendre transformation of the free energy: dF (X) (263) H(X)= F (X) log X X . − · dX In the above we have also considered the φ-entropy Hφ(p): p 1 u v (264) H (p)= p dv du, φ − u2 φ(v) Z1 Z0 where X dw (265) φ(w)= dX = X , dw dX dF (266) w = X . dX 50 JIANZHOU

The following is the main result of this paper: Theorem 6.1. The entropy function H(X) is related to the φ-entropy function Hφ(p) as follows: (267) H(X)= H (w(X)) + [F (X(1)) log X(1)] w(X), φ − · (268) H (p)= F (X(p)) p log X(p) p[F (X(1)) log X(1)]. φ − · − − Proof. The proof is similar to the proof of (246): p 1 H (p) = p F (X(u))d φ u Z1 F (X(u)) p p 1 = p p dF (X(u)) · u − u 1 Z1 p 1 dF (X(w)) = F (X(p)) pF (X(1)) p dX(w) − − w dX(w) · Z1 X(p) 1 = F (X(p)) pF (X(1)) p dX − − X · ZX(1) = F (X(p)) pF (X(1)) p log X(p)+ p log X(1). − − ·

Corollary 6.1. The following identity holds: d (269) H (p)= log X(p) [F (X(1)) log X(1)]. dp φ − − −

7. Conclusions and Prospects In this paper we first establish a link between interpolating statistics with umbral calculus. This link inspires us to reexamine umbral calculus from the point of view of mathematical physics. An unexpected outcome is that this link also leads to a connection of interpolating statistics with generalized statistical mechanics, this not only opens the door of applications of generalized statistical mechanics to fractional quantum Hall effects, but also enables us to understand the generalized entropy functions as the Legendre transformations of the free energy of interpolating statistics. Originally we are led to this work by spectral curves in Eynard- Orantin topological recursions, which we understand as the genus zero one-point functions in various Gromov-Witten type theory. The theory of interpolating statistics are concerned with one-particle partition MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 51 functions, so this leads to our definition of the spectral curves associated with interpolating statistics. The original motivation to use umbral calculus is to find more examples of interpolating statistics. For- tunately some familiarity with Ramanujan’s second notebook enables us to identify their spectral curves in the setting of Eynard-Orantin topological recursions. We have recorded some examples in Section A.8 and Section A.10. As a result, the role of umbral calculus has changed in several aspects. First, in the beginning it is used as a technical tool to find more examples, but soon it turns out the connections with interpolating statistics and generalized entropy suggest to reexamine the umbral calculus from the point of view of quantum mechanics and statistical physics. We hope this interaction with mathematical physics will lead to new progresses in this important branch of combinatorics. Secondly, because for some examples the spectral curves also arise in topological string theory, so it is natural to expect to find deeper connections between the objects in this paper and string theory. A string theoretical interpretation of either interpolating statistics, or generalized entropy, or umbral calculus, is of course very desirable. We hope to address this in future work. Thirdly, umbral calculus can be regarded as belonging to formal algebraic geometry, e.g. to formal group laws. The latter is related complex cobordism theory and Hirzebruch genera. Hence through this chain of connections, one sees that interpolating statistics are related to deep theories in algebraic topology. This partly reflects that the topological nature of fractional quantum Hall effect. We will elaborate on this point in subsequent work. We remark that it seems to be natural to use generalized statistical physics to study interpolating statistics, in particular the fractional quantum effect. This is because fractional quantum Hall effect is an example of topological orders. Microscopically, topological orders correspond to patterns of long-range quantum entanglement. On the other hand, as Tsallis remarked in [43] nonextensive statistical mechanics seem to be more suitable to describe long-range interactions. In the community of cybernetics and information at least 25 different entropy functions have been advanced (see e.g. Taneja [39]). These functions are introduced for various specific purposes, e.g., image pro- cessing. Except for the Shannon entropy, the other 24 examples in loc. cit. are different from the examples in the following Appendix. It is not clear at present whether the examples in this paper are useful in that area. 52 JIANZHOU

Finally, the examples in the Appendix lead us to connections with operads and cluster algebras. We will report such connections in subsequent work.

Acknowledgements. The author is partly supported by NSFC grants 11661131005 and 11890662. The author thanks Professor Sen Hu and Professor Guowu Meng for introducing him to fractional quantum Hall eﬀects. This is an expanded version of a manuscript that has a limited circulation. Some of the results were reported in a colloquium talk at Peking University in 2018. The author thanks Professor Hunjun Fan for the invitation and the audience for their interest.

Appendix A. Examples of Interpolating Statistics from Umbral Calculus In this Appendix we present some examples from umbral calculus. Most of them can be found in Chapter 3 of Ramanujan’s second notebook [3], or Chapter 4 of Roman’s book on umbral calculus [34]. They are reproduced here because of the following two reasons. On the one hand we want to provide some extra information about their related spectral curves; it is a surprise that some of them reproduce the spectral curves that appear in the literature on Eynard-Orantin topological recursions. On the other hand, we will present the calculations related to φ-logarithm, φ-exponential, and φ-entropy. The motivation is to illustrate the relationship between umbral calculus and generalized entropy studied in §6 by concrete examples. We will also present some examples not in Roman [34]. One of them comes from the dilogarithm, and three others are obtained by averag- ing the Acharya-Swamy statistics. In the former example, a famous identity of Euler for dilogarithm naturally arises, and the φ-entropy is very close to the Rogers dilogarithm. In the latter examples, several integer sequences of rich combinatorial signiﬁcance appear in the computations. We also include the Gentile statistics and an example motivated by it. Some of our examples can be identiﬁed with the examples in the Appendix to Taylor [40] where some some more interesting examples can be found. Let us recall the notations. If γn(x) is a polynomial sequence of binomial type, then { }

∞ Xn (270) γ (x) = exp(xF (X)) n n! n=0 X MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 53 for a formal power series F (X)= X+O(X2). The relationship between the free energy F (X) and γ is given in the following two identities: { n} n k x k (271) γ (x)= n! (F (X) ) n , n · k! |X Xk=1 n where ( ) n means the coeﬃcient of X , and · X ∞ Xn (272) F (X)= γ′ (0) . n n! n=0 X Furthermore, let f(Y ) be the compositional inverse series of F (X), i.e., f(F (X)) = X, F (f(Y )) = Y. Then one has ∞ f(Y )n (273) exY = γ (x) n n! n=0 X and the following recursion relations are satisﬁed:

(274) f(D)γn(x)= n γn 1(x), − · 1 (275) γn(x)= x[f ′(D)]− γn 1(x), n 1, − ≥ where D = d/dx. To find the series φ, we first find: d (276) w(X)= X F (X), dX then we apply Lagrange inversion to get: (277) X = X(w). Next apply (240) to get φ: 1 X(u) (278) φ(u)= d = d . du log X(u) du X(u) The φ-logarithm is computed from φ(u) by the following formula: p du (279) ln (p)= , φ φ(u) Z1 the result is X(p) (280) ln (p) = log . φ X(1) 54 JIANZHOU

The φ-exponential function is the inverse function of lnφ. To get the ϕ-entropy, we ﬁrst compute u vdv (281) ξ (u)= , φ φ(v) Z0 the result is

(282) ξφ(u)= F (X(u)). Next, the φ-entropy can be computed as follows: p ξ(u) p F (X(u)) (283) H (p)= p du = p du. φ − · u2 − · u2 Z1 Z1 The result is (284) H (p)= F (X(p)) p log X(p) p[F (X(1)) log X(1)]. φ − · − − A.1. The Boltzmann-Gibbs statistics. In this case, BG n BG BG BG γk (n)= x , F (X)= X, w (X)= X, f (Y )= Y. The spectral curves are given by: z = eX , X = Y. For this example, wBG(X)= X XBG(w)= w, it follows that φBG(p)= p, ξBG(u)= u. The corresponding entropy is given by: p 1 HBG(p)= p udu = p log p. − u2 · − Z1 A.2. The Fermi-Dirac statistics. In this case we have zx = (1+ X)x, γFD(x)=(x) = x(x 1)(x 2) (x (n 1)), n n − − ··· − − n 1 ∞ ( 1) − F FD(X) = ln(1 + X)= − Xn, n n=1 X f FD(Y )= eY 1. − The spectral curves are given by: z =1+ X, X = eY 1. − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 55

The identity zx zy = zx+y is equivalent to the Chu-Vandermonde identity: · x y x + y (285) = . i j n i+j=n X In this case we have X w wFD(X)= , XFD(w)= . 1+ X 1 w − From this we compute that: 1 φFD(u) = = u(1 u), d FD du log X (u) − u v ξFD(u) = dv = log(1 u)= F FD(XFD(u)). v(1 v) − − Z − The φFD-logarithm is p 1 p FD ln FD (p)= du = log = log X (p). φ u(1 u) 1 p Z − − The φFD-exponential is given by: eX 1 (286) p = X = X . e +1 1+ e− The φFD-entropy is: p 1 HFD(p) = p log(1 u)du u2 − Z1 = p log p (1 p) log(1 p). − − − − A.3. The Bose-Einstein statistics. In this case we have: 1 zx = , (1 X)x − γBE(x)= x(n) = x(x + 1)(x + 2) (x +(n 1)), n ··· − ∞ 1 F BE(X)= ln(1 X)= Xn, − − n n=1 BE Y X f (Y )=1 e− . − The spectral curves are given by:

1 Y z = , X =1 e− . 1 X − − 56 JIANZHOU

The identity zx zy = zx+y is equivalent to the Chu-Vandermonde identities with x,· y changed to x, y respectively: − − x y (x + y) (287) ( 1)i − ( 1)j − =( 1)n − . − i · − j − n i+j=n X In this case we have: X wBE wBE(X)= , XBE(wBE)= . 1 X 1+ wBE − So we have φBE(p)= p(1 + p), ξBE(u) = log(1 + u). The φBE-logarithm is p 1 p BE ln BE (p)= du = log = log X (p). φ u(1 + u) 1+ p Z The φBE-exponential is given by: eX 1 (288) p = = . 1 eX e X 1 − − − The φBE-entropy is: p 1 HBE(p) = p log(1 + v)dv − v2 Z1 = p log p +(1+ p) log(1 + p) 2 log 2. − − A.4. The Acharya-Swamy statistics. In this case one has: γAS(x; ǫ)= x(x + ǫ)(x +2ǫ) (x +(n 1)ǫ), n ··· − ∞ n 1 AS 1 ( 1) − n 1 n F (X; ǫ)= ln(1 + ǫX)= − ǫ − X , ǫ n n=1 X eǫY 1 f AS(Y ; ǫ)= − . ǫ The spectral curves are given by: eǫY 1 z =(1+ ǫX)1/ǫ, X = − . ǫ By taking the logarithmic derivative of F AS we get: X wAS(X; ǫ)= , 1+ ǫX By taking the Lagrange inversion one gets: wAS XAS = . 1 ǫwAS − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 57

A combinatorial interpretation of this by counting trees is given by Parker’s Theorem (see e.g. Gessel [16]). From these we have the following computations:

AS 1 AS 1 φ (p)= d = p(1 ǫp), ξ (u)= log(1 ǫu). dp log XφAS (p) − − ǫ −

p 1 p AS ln AS (p) = du = log = log X (p), φ u(1 ǫu) 1 ǫp Z − − eX 1 p = X = X , 1+ ǫe e− + ǫ p p 1 HAS(p) = log(1 ǫv)dv ǫ v2 − Z1 1 1 = p log p (1 ǫp) log(1 ǫp)+ (1 ǫ) log(1 ǫ). − − ǫ − − ǫ − − A.5. The Gentile statistics. Historically the ﬁrst statistics that was proposed to interpolate between the Bose-Einstein and the Fermi-Dirac statistics is the Gentile statistics. It is a family of one-particle partition functions: (289) z (x)=1+ X + X2 + + Xp, (p) ··· with the following property: FD BE (290) z(1)(X)= z (X), lim z(p)(X)= z (X). p + → ∞ One can compute wn = pn(X) as follows: 1 Xp+1 F = log(1 + X + + Xp) = log − ··· 1 X − ∞ 1 = (Xn X(p+1)n), n − n=1 X and so ∞ X Xp+1 w(X) = (Xn (p + 1)X(p+1)n)= (p + 1) . − 1 X − 1 Xp+1 n=1 X − − A.6. An example motivated by Gentile statistics. We take X (291) F (X)= X + X2 + = . ··· 1 X − Its inverse series is X (292) f(X)= X X2 + = . − ··· 1+ X 58 JIANZHOU

The exponential generating series of the conjugate series of F (X) is

∞ Xn xX (293) γ (x) = exp . n n! 1 X n=0 X − The following are the ﬁrst few examples of γn(X):

γ1(x)= x, 2 γ2(x)= x +2x, 3 2 γ3(x)= x +6x +6x, 4 3 2 γ4(x)= x + 12x + 36x + 24x, for the combinatorial meaning of the coefficients, see A089231, A111596, A066667, A008297 on [37]. They are related to the Lah numbers and the Laguerre polynomials of order 1. See Taylor’s thesis [40, Section 3.4 and p. 97] for its relation to counting− permutations. For this example, we have X w(X)= , (1 X)2 − and so 1+2w √1+4w X = − 2w = w 2w2 +5w3 14w4 + 42w5 132w6 + 429w7 1430w8 + , − − − − ··· where the coefficients of the series expansion are the Catalan numbers up to signs. dw X(1 + X) φ(w) = X = = w√1+4w dX (1 X)3 − = w +2w2 2w3 +4w4 10w5 + 28w6 84w7 + . − − − ··· Up to signs, the coefficients are A002420 on [37]. w w w 1 X w ξφ(w) = dw = dw = dX φ(w) √1+4w X Z0 Z0 Z0 √1+4w 1 = − = F (X). 2 p ξ (w) w √1+4w 1 H (p) = p φ dw = p − dw φ − w2 − 2w2 Z1 Z1 1 2p √1+4p 1+2p √1+4p = − − p log − + cp 2 − 2p for some constant p. MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 59

A.7. The exponential polynomials. For this sequence one has: F (X)= eX 1, − f(Y ) = log(1 + Y ), a combinatorial interpretation by counting trees is given by Drake’s Theorem [13] (see also Gessel [16]). The sequence is given by n k γn(x)= S(n, k)x , Xk=0 where S(n, k) are the Stirling numbers of the second kind: k 1 k j k n (294) S(n, k)= ( 1) − j . k! − j j=0 X The spectral curves are given by: eX 1 X z = e − ,Y = e 1, X = log(1 + Y ). − In this case, dF X w = X dX = Xe , By Lagrange inversion formula: n 1 n 1 ∞ ( 1) − n − X = − wn. n! n=1 X This is essentially Cayley’s famous formula for counting rooted trees. Some other relevant computations are:

1 dw X 2 X φ(w) = d = X = Xe + X e dw log X dX n 1 n 1 ∞ ( 1) − n − = w 1+ − wn . n! n=1 X w w w w X w ξφ(w) = dw = dw = dX φ(w) X dw X Z0 Z0 dX Z0 X = eX dX = eX 1= F (X) − Z0 ∞ n 1 n 1 (n 1) − n = ( 1) − − w . − n! n=1 X log X = log w X − 60 JIANZHOU

p 1 H(p) = p (eX(w) 1)dw − w2 − Z0 n 2 ∞ (n 1) − = p log p + ( 1)n − pn. − − n! n=2 X

A.8. Abel polynomials and Lambert series. In 1826 Abel proved the following deep generalization of the binomial identity (cf. [12, §3.1]):

n n n k 1 n k (295) (x + y) = x(x ka) − (y + ka) − . k − Xk=0 This leads to the sequence of Abel polynomials:

Abel n 1 γ (x)= x(x na) − , n − f Abel(Y )= Y eaY , n 1 ∞ ( an) − F Abel(X)= − Xn, n! n=1 X n 1 ∞ ( an) − wAbel(X)= − Xn. (n 1)! n=1 X − This example is related to Cayley’s formula for rooted trees and Entry 13 of Ramanujan’s second notebook [3]. The spectral curves are:

n 1 ∞ ( an) − z = exp − Xn, X = Y eaY . n! n=1 X Spectral curves of this form coincide exactly with the spectral curves in the studies of Eynard-Orantin topological recursions of Hurwitz numbers [6, 15] and its generalizations [7].

∞ n +1 X = wAbel exp an(wAbel)n · n n=1 X wAbel 1 = exp 1 . 1 awAbel 1 awAbel − − − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 61

Abel p 2 φ (p) = n n = p(1 ap) , 1+ n∞=1(n + 1)a p − u v 1 1 u ξAbel(u) = P dv = = , v(1 av)2 a(1 au) − a 1 au Z0 − − − p 1 p 1 logAbel(p) = dv = log + , v(1 av)2 1 ap 1 ap Z1 − − − p 1 1 1 HAbel(p) = p du − u2 a(1 au) − a Z1 − = p log p + p log(1 ap) p log(1 a). − − − − xf Abel(X) γ∨(X) n ∞ For e = n=0 n! X , see Taylor [40, §3.3 and p. 96]. They are related to A059297 on [37]. P A.9. Gould polynomials. After a change of variables from the original definition in Roman [34], Gould polynomials are defined to be the following sequence of polynomials: n 1 − (296) γGould(x; a, b)= x (x an jb). n − − j=1 Y They form the sequence of polynomials of binomial type associated with ebY 1 (297) f Gould(Y ; a, b)= eaY − , b =0, b 6 whose compositional inverse series is k 1 ∞ − k Gould k 1 X (298) F (X; a, b)= ( 1) − (ak + jb) . − · k! j=1 Xk=1 Y In particular, n 1 ∞ − x (x an bj) Xn − − · n=0 j=1 X Y k 1 ∞ − k k 1 X = exp x ( 1) − (ak + jb) . − · k! j=1 Xk=1 Y The spectral curves are k 1 ∞ − k k 1 X (299) z = exp ( 1) − (ak + jb) − · k! j=1 Xk=1 Y 62 JIANZHOU and ebY 1 (300) X = eaY − . b The latter can be rewritten as follows: aY bY (301) bXe− e +1=0. − One can compare it with Entry 14 in Chapter 3 of Ramanujan’s second notebook [3, (14.4)]. One also has: k 1 ∞ − k Gould k 1 X (302) w (X; a, b)= ( 1) − k (ak + jb) . − · k! j=1 Xk=1 Y It satisfies the following equation: w(1 aw)a/b (303) X = − . (1 (a + b)w)(a+b)/b − This generalizes Wu’s formula (a = α, b = 1): − w X = (1 αw)α[1 + (1 α)w]1 α − − − for Haldane-Wu statistics [19, 49]. We then have: φGould(p) = p(1 ap)(1 (a + b)p), − − u v ξGould(u) = dv v(1 av)(1 (a + b)v) Z0 − − 1 = (log(1 au) log(1 (a + b)u)), b − − − p 1 1 HGould(p) = p (log(1 au) log(1 (a + b)u))du − u2 · b − − − Z 1 = p log p + (1 ap) log(1 ap) − b − − (1 (a + b)p) log(1 (a + b)p) . − − − A.10. Specializations of the spectral curve associated with Gould polynomials. Let us consider the following specializations of the spectral curve (301). Case 1. By taking a = 0, b = ǫ, one gets eǫY 1 (304) X = − , ǫ hence one recovers the spectral curve associated with the Acharya- Swamy statistics. MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 63

Case 2. By taking b 0, one gets → (305) X = Y eaY , and one recovers the Lambert curve and its generalizations mentioned in §A.8. Case 3. When a = g 1, b = 1, one gets the curve: − gY (g 1)Y (306) X = e e − . − This is the spectral curve for the framed topological vertex [6, 9, 51]. Case 4. When b = 2a, one gets − aY aY (307) 2aX = e− e . − − After suitable change of coordinates, one gets the Catalan curve: 1 (308) z = w − w that appears in Eynard-Orantin topological recursions of several geometric problems. A.11. Mittag-Leffler polynomials. For this sequence one has: eX 1 f ML(X) = − , eX +1 1+ X F ML(X) = log , 1 X − n n n 1 γML(x) = − 2k(x) . n k n k k Xk=0 − The exponential generating series of Mittag-Leffler polynomials is x ∞ γ (x) 1+ t (309) n tn = n! 1 t n=0 X − By taking t = X and x = 1, one sees that the spectral curves are: 1+ X 1+ X eY 1 z = ,Y = log , X = − . 1 X 1 X eY +1 − − For simplification of notations we take 1+ X/2 F ML(X) = log . 1 X/2 − Then we have

ML X w = 2 , 1 X − 4 64 JIANZHOU and so

2 ∞ ML 2(√1+ w 1) n 1 (2n 2)! 2n 1 X (w) = − =2 ( 1) − − (w/2) − , w − n!(n 1)! n=1 X − (2n 2)! ML − where the coeﬃcients n!(n 1)! are the Catalan numbers, so X is essentially a generating function− of the Catalan numbers. 1 φML(p) = = p p2 +1 d ML p dp log X (p) p ∞ n 1 (2n 2)! 2n = p 1+2 ( 1) − − (p/2) . − n!(n 1)! n=1 X − u v u dv ξML(u) = dv = = log(u + √1+ u2) φML(v) √v2 +1 Z0 Z0 ∞ n 1 (2n 2)! 2n 1 = 2 ( 1) − − (u/2) − . − (2n 1) (n 1)!2 n=1 X − · −

2 ML 2( 1+ p 1) log ML (p) = log X (p) = log − φ p p 2n ∞ ( 1)n = log p + − n p2n. 22n 2n n=1 X

p ξML(u) p log(u + √1+ u2) HML(p) = p du = p du − u2 − u2 Z Z (1 + 1+ p2) log(p + 1+ p2) = p log p + log + 1 − 2 p − p p ∞ n 1 ( 1) − (2n 2)! 2n 2 = p log p + − − u − . − 22n 2 (2n 1)(2n 2) (n 1)!2 n=2 − X − − · − A.12. Bessel polynomials. For this sequence, f Bessel(t)= t t2/2, − F Bessel(t)=1 √1 2t, −n − n k (2n k 1)! 1 − γBessel(X)= − − xk. n (k 1)!(n k)! 2 Xk=1 − − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 65

The exponential generating function is

∞ γn(x) n x(1 (1 2t)1/2 ) t = e − − , n! n=0 X and so the spectral curves are 2 1 (1 2X)1/2 1/2 Y z = e − − ,Y =1 (1 2X) , X = Y . − − − 2 We have: X wBessel = , √1 2X − XBessel(w) = w(√1+ w2 w), − φBessel(p) = p 1+ p2( 1+ p2 + p). p p v v ξBessel(u) = dv v√1+ v2(√1+ v2 + v) Z1 = u +1 √1+ u2. − One again encounters Catalan numbers in their expansions. We also have:

2 log Bessel (p) = log X(p) = log[p( 1+ p p)] φ − 2n ∞ ( 1)np = log p − n p2n+1. − 22n 2n +1 n=0 X p u +1 √1+ u2 HBessel(p) = p − du − u2 Z = p log p 1+ p2 + p log(p + 1+ p2)+1 − − 2n ∞ ( 1)n = p log p + p − n p p2n+2. − 22n (2n + 1)(2n + 2) n=0 X xf Bessel(X) γ∨(X) n ∞ For e = n=0 n! X , see Taylor [40, §3.2 and p. 96].

A.13. Mott polynomials.P The Mott polynomials (cf. [?, p. 251] and [34, §4.12]) are deﬁned by:

(1 t2)1/2 1 ∞ (310) exp x − − = g (x)tn, · t n n=0 X 66 JIANZHOU where gn(x) are explicitly given by: [n/2] x 2l g (x)=( x/2)n(n 1)! − . n − − l!(n l)!(n 2l 1)! Xl=0 − − − hence they are the associated polynomials with (1 t2)1/2 1 F (t)= − − , t and so 2t f(t)= − . 1+ t2 In (310), change t to 2t, − 1 (1 4t2)1/2 ∞ γ (x) (311) exp x − − = n tn, · 2t n! n=0 X where for n 0, ≥ [n/2] xn 2l γ (x)= n!(n 1)! − n − l!(n l)!(n 2l 1)! Xl=0 − − − form the associated sequence with 1 (1 4t2)1/2 F Mott(t)= − − , 2t t f Mott(t)= . 1+ t2 The spectral curves are: 1 (1 4X2)1/2 1 (1 4X2)1/2 Y z = exp − − ,Y = − − , X = . 2X 2X 1+ Y 2

2 ∞ Mott 1 √1 4X 2n 1 2n 1 (312) w (X)= − − = − X − . √ 2 n 2X 1 4X n=1 − X By Cardano formula, 1 XMott(w) = ( (1+18w2 +33/2w√1+11w2 w4)1/3 6w − 2 1+3w 2 3/2 2 4 1/3 1 + (1+18w +3 w√1+11w w )− . 6w − − 3w The ﬁrst few terms of the series expansion are: XMott(w)= w 3w3 + 17w5 119w7 + 929w9 7755w11 + − − − ··· MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 67

The coefficients 1, 3, 17, 119, 929, 7755,... do not form a sequence on [37]. The φ-function− in− this case − 1 φMott(w)= d Mott dw log X (w) has a complicated explicit expression. The first terms of its series expansion are: φMott(w) = w +6w3 14w5 + 78w7 542w9 + 4214w11 + . − − ··· The sequence of coefficients 1, 6, 14, 78, 542, 4214,... does not appear on [37]. The expression of φMott in terms of X = XMott is given explicitly as follows: XMott dwMott φMott(w) = = X dXMott(w) dX dw 2 1 1 = + 1 (1 4X2)3/2 2X2 − √1 4X2 − − ∞ 2n +1 = (2n + 1) X2n+1 n n=0 X = 1+9X2 + 50X4 + 245X6 + 1134X8 + 5082X10 + . ··· The coefficients 1, 9, 50, 245, 1134,... are the sequence A001818 on OEIS [37]. They are related to the Catalan numbers as follows: 2n +1 1 2n (2n + 1) = (2n + 1)2 . n · n +1 n The ξ-function in this case is: w w w dX ξMott(w) = dw = w dw dX φMott(w) · X Z Z w w u 1 √1 4X2 = dX = − − dX X 2X2√1 4X2 Z Z − 1 √1 4X2 ∞ 1 2n = − − = Xn+1, 2X n +1 n n=0 X 1 2n where the coefficients n+1 n are the Catalan numbers. To expand ξMott(w) to a power series in w, we find: (313) ξMott(w)= w 1 X(w)2. − p 68 JIANZHOU

Write √1 X2 = Y . Then by (312) − 1 Y w2 = − . (1 Y 2)Y 2 − By Cardano formula we find 1 Y = ( w3 +3√ 3w4 + 33w2 +3+18w)1/3 3w − − 2 w 3 3 4 2 1/3 1 + − ( w +3√ 3w + 33w +3+18w)− 3w − − − 3 = 1 2w2 + 10w4 66w6 + 498w8 4066w10 + 34970w12 . − − − −··· The coefficients 1, 2, 10, 66, 498,... are the sequence A027307 on OEIS [37]. These are the 3-Schröder numbers according to Yang-Jiang [50], i.e., the number of paths from (0, 0) to (3n, 0) that stay in first quadrant (but may touch horizontal axis) and where each step is (2, 1), (1, 2) or (1, 1). Therefore, (313) and (312) establish a relation between the Catalan− numbers and the 3-Schröder numbers. See also Drake’s Example 1.6.9 [13]. The entropy function in this case is p 1 HMott(p) = p ξMott(w)dw − w2 Z p 1 1 1 4X2(w) = p − − dw − w2 2X(w) Z p X(p) 1 1 √1 4X2 1 √1 4X2 = p − − d − − 2 2X √ 2 − 1 √1 4X2 2X 1 4X Z − − − 2X√1 4X2 − = p log X(p)+ p 1 4X2(p). − − One can directly check that p dF Mott(X) HMott(p)= F Mott(X) log X X . − · dX

A.14. Conjugate sequence of dilogarithm function. In this Sub- section we present an example not included in Roman [34]. Consider the dilogarithm function:

∞ tn (314) F (t) = Li (t)= . 2 n2 n=1 X MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 69

Although the series is only convergent for t < 1, but it can be rewritten as an integral: | | t log(1 z) t dz z dw (315) − dz = , − z z 1 w Z0 Z0 Z0 − one can use these integral to extend the definition of Li2(t) to a multi- valued function on C 1 . The conjugate sequence−{ } γ (x) for the series F (X) is defined by { n } ∞ tn ∞ tn γ (x) = exp x . n n! n2 n=0 n=1 X X The following are the first few terms of γn:

γ1(x)= x, 1 γ (x)= x2 + x, 2 2 3 2 γ (x)= x3 + x2 + x, 3 2 3 41 3 γ (x)= x4 +3x3 + x2 + x. 4 12 2 To find the series φ, we first find:

d ∞ Xn 1 (316) w(X)= X F (X)= = log , dX n 1 X n=1 X − then we apply Lagrange inversion to get: w (317) X = X(w)=1 e− . − Next apply (240) to get φ:

1 X(u) u (318) φ(u)= d = d = e 1. du log X(u) du X(u) − The φ-logarithm is computed from φ(u) by the following formula: p du p u ep 1 e 1 (319) ln (p)= = du = log − log − , φ φ(u) eu 1 ep − e Z1 Z1 − this is again just log X(p) log X(1). To get the φ-entropy, we ﬁrst compute −

p ∞ np ∞ u p e− 1 ξ (p) = du = p log(1 e− ) + φ eu 1 − − n2 n2 Z0 − n=1 n=1 2 X X p p π = p log(1 e− ) Li (e− )+ . − − 2 6 70 JIANZHOU

On the right-hand side of the last equality we see dilogarithm function again. By Proposition 6.4 we know that p (320) ξ (p)= F (X(p)) = Li (1 e− ). φ 2 − They match each other by the following well-known identity due to Euler [31]: π2 (321) Li (z) + Li (1 z)= log z log(1 z). 2 2 − 6 − − Using

u ∞ B = n un, eu 1 n! n=0 − X one also has ∞ B (322) ξ (p)= n pn+1. φ (n + 1)! n=0 X By taking p = 1 we get:

∞ 1 Bn (323) Li (1 e− )= . 2 − (n + 1)! n=0 X Next, the φ-entropy can be computed as follows: p p u ξφ(u) Li2(1 e− ) Hφ(p) = p 2 du = p −2 du − 1 u − 1 u Z u p Z p Li2(1 e− ) 1 u = p − p d Li (1 e− ) · u − · u 2 − 1 Z1 p p 1 1 udu − − = Li2(1 e ) p Li2(1 e ) p u − − · − − · 1 u e 1 p 1 Z −u p = Li2(1 e− ) p Li2(1 e− ) p log(1 e− ) 1 − p − · − 1 − · − p | 1 = Li (1 e− ) p Li (1 e− ) p log(1 e− )+ p log(1 e− ). 2 − − · 2 − − · − · − This can be expanded as follows.

p ∞ d 1 e− 1 1 Bn n 1 log − = = p − , dp p ep 1 − p n! n=1 − X and so p n 1 e− ∞ B p log − = n . p n! n n=1 X MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 71

By taking p = 1 we get:

∞ 1 Bn 1 log(1 e− )= . − n! n n=1 X Hence we have ∞ B ∞ B pn+1 H (p) = p log p + n pn+1 n φ − (n + 1)! − n! n n=0 n=1 X1 1 X p [Li (1 e− ) log(1 e− )] − · 2 − − − ∞ B pn 1 = p ln p p n . − − n! n(n + 1) − n(n + 1) n=1 X On the other hand, by (322) we have:

p ξ (u) ∞ B pn 1 H (p) = p φ du = p ln p p n . φ − u2 − − n! n(n + 1) − n(n + 1) 1 n=1 Z X It is a match. In this case we have 1 H(X) = F (X) w(X) log X = Li (X) log X log . − 2 − · 1 X − This is very close to Rogers dilogarithm: 1 (324) L(z) = Li (z)+ log z log(1 z). 2 2 − Now p p H(X(p)) = Li (1 e− ) p log(1 e− ). 2 − − − So we get: 1 1 H(X(p)) = H (p)+ p[Li (1 e− ) log(1 e− )]. φ 2 − − − This matches with (268).

A.15. Averaged Acharya-Swamy statistics. In this Subsection we present three more examples not included in Roman [34]. For the ﬁrst example, we take: 1 Y = F (X)= (log(1 + ǫX) log(1 ǫX)), 2ǫ − − 1 when ǫ = 2 this reduces to the case of Mittag-Leﬄer polynomials. We have ǫY ǫY e e− X = f(Y )= ǫY − ǫY , e + e− 72 JIANZHOU

d X 1 1 X w(X) = X F (X)= + = . dX 2 1+ ǫX 1 ǫX 1 ǫ2X2 − − From this we solve for w to get: √1+4ǫ2w2 1 X(w) = − . 2ǫ2w It has the following series expansion:

∞ 1 2n X(w)= ( 1)n ǫ2nw2n+1, − n +1 n n=0 X 1 2n the coeﬃcients n+1 n are the Catalan numbers. Next apply (240) to get φ: 1 X(u) √ 2 2 (325) φ(u)= d = d = u 1+4ǫ u . du log X(u) du X(u) It has the following expansion:

∞ 1 2n φ(u) = u +2u ( 1)n ǫ2n+2u2n+2, − n +1 n n=0 X where we again see the Catalan numbers. The φ-logarithm is computed from φ(u) by the following formula: p du p 1 log (p) = = du φ φ(u) u√1+4ǫ2u Z1 Z1 1+4ǫ2p2 1 √1+4ǫ2 1 = log 2 − log 2 − , p 2ǫ p − 2ǫ this is again just log X(p) log X(1). It has the following expansion: − √1+4ǫ2 1 ∞ 1 2n 1 log (p) = log p log − + ( 1)n − ǫ2np2n, φ − 2ǫ2 − n n n=1 X 2n 1 the coeﬃcients n− are the sequence A001700 on [37]. Note 1 2n 1 1 2n − = , n n 2n n 2n where n are the sequence A000984 on [37]. To get the φ-entropy, we ﬁrst compute p p udu u 1 2 2 ξφ(p) = = du = log( 1+4ǫ p +2ǫp). 2 0 φ(u) 0 u√1+4ǫ u 2ǫ Z Z p MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 73

One can check that it coincides with F (X(p)). It has the following expansion:

∞ ( 1)n 2n ξ (p)= − ǫ2np2n+1, φ 2n +1 n n=0 X

2n the coeﬃcients n are A000984 on [37]. Next, the φ-entropy can be computed as follows: p ξ (u) p 1 H (p) = p φ du = p log(√1+4ǫ2u2 +2ǫu)du. φ − u2 − 2ǫu2 Z1 Z1 We can again check (268) in this case. We have the following expansion:

H(X(p)) = F (X(p)) p log X(p) − ∞ ( 1)n 2n = − ǫ2np2n+1 2n +1 n n=0 X ∞ 1 2n p log p p ( 1)n ǫ2np2n − − · − 2n n n=1 X ∞ 1 2n = p log p ( 1)n ǫ2np2n+1. − − − 2n(2n + 1) n n=1 X Another way to take the average of Acharya-Swamy statistics is as follows:

1 1 1 Y = F (X)= log(1 + ǫX)+ ǫ log(1 + X) , 2 ǫ ǫ

d X 1 1 w(X) = X F (X)= + . dX 2 1+ ǫX 1+ 1 X ǫ

From this we solve for w to get:

1 1 2 2 1 (ǫ + ǫ− )w 1+(ǫ ǫ ) w X(w) = − − − − . 2w (ǫ + ǫ 1) − p − 74 JIANZHOU

It has the following series expansion:

X(w) = w + sw2 + w3 + s( s2 + 2)w4 +( s2 + 2)w5 − − + s(2s4 6s2 + 5)w6 + (2s4 6s2 + 5)w7 − − + s( 5s6 + 20s4 28s2 + 14)w8 − − + ( 5s6 + 20s4 28s2 + 14)w9 − − + (14s9 70s7 + 135s5 120s3 + 42s)w10 − − + (14s8 70s6 + 135s4 120s2 + 42)w11 − − + ( 42s11 + 252s9 616s7 + 770s5 495s3 + 132s)w12 + , − − − ··· where s is deﬁned by ǫ + ǫ 1 s = − . 2 The coeﬃcients appear as A094385, A094385, A157491 or A062991 on [37]. Next apply (240) to get φ: 1 φ(u) = d du log X(u) s = u 1+4(s2 1)u2 + ( 1+4(s2 1)u2 4(s2 1)u2 1). − 2(s2 1) − − − − p − p It has the following expansion:

φ(u) = u + 2[(s2 1)u3 (s2 1)2u5 + 2(s2 1)3u7 5(s2 1)4u9 − − − − − − + 14(s2 1)5u11 + ] − ··· + s[ u2 (s2 1)u4 +2s(s2 1)2u6 5(s2 1)3u8 − − − − − − + +14(s2 1)4u10 + ], − ··· 1 where the coefficients are the Catalan numbers. The function φ(u) has the following expansion: 1 1 = + s +( s2 + 2)u +( 2s3 +3s)u2 + (3s4 8s2 + 6)u3 φ(u) u − − − + (6s5 15s3 + 10s)u4 +( 10s6 + 36s4 45s2 + 20)u5 − − − + ( 20s7 + 70s5 84s3 + 35s)u6 − − + (35s8 160s6 + 280s4 224s2 + 70)u7 − − + (70s9 315s7 + 540s5 420s3 + 126s)u8 − − + ( 126s10 + 700s8 1575s6 + 1800s4 1050s2 + 252)u9 − − − + ( 252s11 + 1386s9 3080s7 + 3465s5 1980s3 + 462s)u10 + , − − − ··· MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 75 to unravel the coefficients, we note s u s − = + s2 1 s(s2 1)u 2(s2 1)2u2 +3s(s2 1)2u3 φ(u) u − − − − − − + 6(s2 1)3u4 10s(s2 1)3u5 20(s2 1)4u6 − − − − − + 35s(s2 1)4u7 + 70(s2 1)5u8 126s(s2 1)5u9 − − − − 252(s2 1)6u10 + , − − ··· the coefficients 1, 1, 2, 3, 6, 10, 20, 35, 70, 126,... are the sequence A001405 n on [37], they are [n/2] for n =0, 1, 2,... . Using also p udu p ξ (u) ξ (p)= , H (p)= p φ du, φ φ(u) φ − u2 Z0 Z after integrations we get the expansion for logφ(p), ξφ and Hφ from the 1 above expansion for φ(u) . Yet another way take the average of Acharya-Swamy statistics is as follows: 1 1 Y = F (X)= log(X + ǫ) + log(X + ) , 2 ǫ

1 1 1 2 X = f(Y )= ǫ + ǫ +4e2Y 2 ǫ − − ǫ s d X 1 1 w(X) = X F (X)= + . dX 2 X + ǫ X + 1 ǫ From this we solve for w to get:

2 ǫ ǫ−1 2 (1 2w) 1 4 − −1 (w w ) X(w) − − s − ǫ+ǫ − = . ǫ + ǫ 1 4(w 1) − − It has the following series expansion:

2X(w) 2 2 2 3 3 1 = (1 t)w + (1 t )w +(1+ t 2t )w ǫ + ǫ− − − − + (1+4t3 5t4)w4 + (1 2t3 + 15t4 14t5)w5 − − − + (1 15t4 + 56t5 42t6)w6 − − + (1+5t4 84t5 + 210t6 132t7)w7 − − + (1 420t6 + 56t5 + 792t7 429t8)w8 − − + (1 14t5 + 420t6 1980t7 + 3003t8 1430t9)w9 + , − − − ··· 76 JIANZHOU where t is deﬁned by

1 2 ǫ ǫ− t = − , ǫ + ǫ 1 − the coeﬃcients 14, 420, 1980, 3003, 1430 etc. are the sequence A117434 on [37]. Next apply (240) to get φ:

1 1 2 2 φ(u)= d = − (1 4t(u u ) 1 4t(u u )). du log X(u) 2t − − − − − p It has the following expansion:

φ(u) = u (t + 1)u2 (2t2 2t)u3 (5t3 6t2 + t)u4 − − − − − (14t4 20t3 +6t2)u5 (42t6 70t5 + 30t3 2t2)u6 − − − − − (132t6 252t5 + 140t4 20t3)u7 − − − (429t7 924t6 + 630t5 140t4 +5t3)u8 − − − (1430t8 3432t7 + 2772t6 840t5 + 70t4)u9 + , − − − ··· where the coeﬃcients 1430, 3432, 2772, 840, 70 etc are the sequence A068763. In this case X(1) = we weil deﬁne the φ-logarithm logφ(p) to be just ∞ 1 log X(p). The function φ(u) has the following expansion:

1 1 = +(t +1)+(3t2 + 1)u + (10t3 3t2 + 1)u2 φ(u) u − + (35t4 20t3 + 1)u3 + (126t5 105t4 + 10t3 + 1)u4 − − + (462t6 504t5 + 105t4 + 1)u5 − + (1716t7 2310t6 + 756t5 35t4 + 1)u6 − − + (6435t8 10296t7 + 4620t6 504t5 + 1)u7 − − + (24310t9 45045t8 + 25740t7 4620t6 + 126t5 + 1)u8 + , − − ··· the coefficients 1, 3, 10, 126,... are the sequence A001700 on [37], they 2n 1 are n− for n =1, 2,... . The coefficients 3, 20, 105, 504, 2310,... are the sequence A000917 on [37]. They are (2n 1)!/(n!(n 2)!) for n = 2, 3,4,.... The sequence of coefficients 10, 105−, 756, 4620,−25740,... are not on [37]. However we note

1 1 4t(u u2)+ (1 4t(u u2) = − − − − , φ(u) 2u(1 u)(1 4t(u u2)) − −p − MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 77 hence 2u(1 u) 1 − = 1+ φ(u) 1 2 t1/2 (2t1/2u)+(2t1/2u)2 − · · 1 1 = 2+ tp(2u)+ ( t +3t2)(2u)2 + ( 3t2 +5t3)(2u)3 · 2 − 2 − 1 1 + (3t2 30t3 + 35t4)(2u)4 + (15t3 70t4 + 63t5)(2u)5 8 − 8 − 1 + ( 5t3 + 105t4 315t5 + 231t6)(2u)6 16 − − 1 + ( 35t4 + 315t5 693t6 + 429t7)(2u)7 16 − − 1 + (35t4 1260t5 + 6930t6 12012t7 + 6435t8)(2u)8 128 − − 1 + ( 9240t6 + 630t5 + 36036t7 51480t8 + 24310t9)(2u)9 + , 128 − − ··· which is essentially the generating function of the Legendre polynomials. Again after integrations we get the expansion for logφ(p), ξφ and 1 Hφ from the above expansion for φ(u) .

A.16. Bell polynomials and universal interpolating statistics. The incomplete Bell polynomials are deﬁned as follows:

Bn,k(t1,...,tn k+1) − m1 mn−k+1 (326) n! t1 tn k+1 = − , m1! mn k+1! 1! ··· (n k + 1)! X ··· − − where the sum is taken over all sequences m1, m2, m3,...,mn k+1 of non-negative integers such that these two conditions are satisﬁe−d:

m1 + m2 + + mn k+1 = k, ··· − m1 +2m2 + +(n k + 1)mn k+1 = n. ··· − − From this deﬁnition it is easy to see that when

∞ Xj (327) F Bell(X)= t , j j! j=1 X one has n ∞ n Bell X k (328) exp(xF (X))=1+ Bn,k(t1,...,tn k+1)x . n! − n=1 X Xk=1 78 JIANZHOU

Therefore, if one takes t1 = 1, then n (n) k (329) γ (x; t)= Bn,k(t1,...,tn k+1)x − Xk=1 defines a polynomial sequence of binomial type. The spectral curves in this case are t t (330) z = exp(X + j Xj),Y = X + j Xj. j! j! j>1 j>1 X X These give the indication that the Bell polynomials give the universal interpolating statistics and the spectral curves can be understood as the universal spectral curve. All other examples are given by choosing tj’s in a suitable way. In this case, the Lagrange inversion formula that finds the coefficient sj in the inverse series ∞ Xj (331) f Bell(X)= s j j! j=1 X from (327) has different combinatorial interpretations involving sum- mations over trees and has appeared in the literature in many times. For example, in Chen’s thesis [10] in 1990, it is interpreted as a summation over all Schröder trees on n vertices. In his 2008 thesis, Drake [13, Example 1.4.7] interpreted it as a summation over phylogenetic trees. This was rediscovered in a 2016 paper by Engbers-Galvin-Smyth [14, Lemma 4.2]. This is a special case of the general formulas for Lagrange inversion in terms of summation over trees published by Wright [48] in 1989. Wright’s formula was rediscovered by Ginzburg and Kapranov [18, Theorem 3.3.9] in the setting of Kozsul duality of operads. This suggests a possible connection of this work with the theory of operads. Indeed, Gingzburg and Kapranov defined the generating series of a k-linear operad (n) by P ∞ xn (332) g (x)= dim (n) . P P · n! n=1 X Denote by s, om, ie the associative operad, the commutative operad, and theA LieC operadL respectively. Then by [18, Example (3.1.12)(a)]:

x x g s(x)= , g om(x)= e 1, g ie(x)= log(1 x). A 1 x C − L − − − BE We recognize that g s(X) = w (X) is the weight function for the Bose-Einstein statistics,A it is also the free energy of the example in MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 79

x §A.6; g om(x)= e 1 as the free energy associated with the exponential C − polynomials (cf. §A.7); g ie(x) as the free energy function of the Bose- Einstein statistics. By [18,L Theorem 3.3.9], if is an admissible dg- operad and is its dual, then P Q (333) g ( g ( x)) = x. Q − P − Note s is self-dual, and om and ie are dual to each other. A C L

References [1] R. Acharya, P. N. Swamy, Statistical mechanics of anyons, J. Phys. A Math. Gen. 27 (1994), 7247 -7263. [2] D. Arovas, J. R. Schrieffer, F. Wilczek, Fractional statistics and the quantum Hall effect, Phys. Rev. Lett. 53 (1984), 722-723. [3] B. C. Berndt, R. J. Evans, B. M. Wilson, Chapter 3 of Ramanujans Second Notebook, Adv. Math. 49(1983), 123-169. [4] F. A. Berezin, M. A. Shubin, The Schrödinger equation. Mathematics and its Applications (Soviet Series), 66. Kluwer Academic Publishers Group, Dor- drecht, 1991. [5] M. V. Berry, Quantal Phase Factors Accompanying Adiabatic Changes. Pro- ceedings of the Royal Society A. 392 (1984): 45-57. [6] V. Bouchard, M. Mariño, Hurwitz numbers, matrix models and enumerative geometry. From Hodge theory to integrability and TQFT tt∗-geometry, 263- 283, Proc. Sympos. Pure Math., 78, Amer. Math. Soc., Providence, RI, 2008. [7] V. Bouchard, D. Hernández Serrano, X. Liu, M. Mulase, Mirror symmetry for orbifold Hurwitz numbers. J. Differential Geom. 98 (2014), no. 3, 375-423. [8] S. Chaturvedi, V. Srinivasan, Interpolations between Bose and Fermi statistics, Physica A 246 (1997) 576-586. [9] L. Chen, Bouchard-Klemm-Marino-Pasquetti Conjecture for C3, arXiv:0910.3739 (2009). [10] W. Y. C. Chen, A general bijective algorithm for trees. Proc. Nat. Acad. Sci. U.S.A. 87 (1990), no. 24, 9635-9639. [11] M. Combescure, D. Robert, Coherent states and applications in mathematical physics. Theoretical and Mathematical Physics. Springer, Dordrecht, 2012. [12] L. Comtet, Advanced combinatorics, The art of finite and infinite expansions. D. Reidel Publishing Company, 1974. [13] B. Drake, An inversion theorem for labeled trees and some limits of areas under lattice paths, PhD Thesis, Brandeis University, 2008. [14] J. Engbers, D. Galvin, C. Smyth, Restricted Stirling and Lah number matrices and their inverses. J. Combin. Theory Ser. A 161 (2019), 271-298. [15] B. Eynard, M. Mulase, B. Safnuk, The Laplace transform of the cut-and-join equation and the Bouchard-Mariño conjecture on Hurwitz numbers. Publ. Res. Inst. Math. Sci. 47 (2011), no. 2, 629-670. [16] I. Gessel, The Carlitz-Scoville-Vaughan Theorem and its generalizations, slides for talk at Joint Mathematics Meeting, San Diego, January 12, 2013 [17] J.W. Gibbs, Elementary principles in statistical mechanics developed with special reference to the rational foundation of thermodynamics. Dover, 1960. 80 JIANZHOU

[18] V. Ginzburg, M. Kapranov, Koszul duality for operads, Duke Math. J. 76 (1994), 203-272. [19] F. D. M. Haldane, “Fractional Statistics” in arbitrary dimensions: A generalization of the Pauli Principle, Phys. Rev. Lett. 67 (1991), 937-940. [20] B. C. Hall, Quantum theory for mathematicians. GTM, 267. Springer, New York, 2013. [21] B.I. Halperin, Statistics of quasiparticles and the hierarchy of fractional quan- tized Hall states, Phys. Rev. Lett 52 (1984), 1583-1586. [22] J. Havrda, F. Charvát Quantification method of classification processes. Con- cept of structural α-entropy, Kybernetika, Vol. 3 (1967), No. 1, 30-35. [23] K. Huang, Statistical mechanics. Second edition. John Wiley & Sons, Inc., New York, 1987. [24] E. T. Jaynes, Information Theory and Statistical Mechanics. Physical Review. Series II (1957) 106 (4): 620-630. [25] E. T. Jaynes, Information theory and statistical mechanics. II. Phys. Rev. (2) 108 (1957), 171-190. [26] S. A. Joni, G.-C. Rota, Coalgebras and bialgebras in combinatorics. Umbral calculus and Hopf algebras (Norman, Okla., 1978), pp. 1-47, Contemp. Math., 6, Amer. Math. Soc., Providence, R.I., 1982. [27] A. Khare, Fractional Statistics and Quantum Theory. World Scientific, 2005. [28] J. P.S. Kung, G.-C. Rota, C. HH Yan, Combinatorics: the Rota way. Cam- bridge University Press, Cambridge, 2009. [29] H. B. Laughlin, Anomalous quantum Hall effect: An incompressible quantum fluid with fractionally charged excitations, Phys. Rev. Lett. 50(1983), 1395- 1398. [30] J. M. Leinaas, J. Myrheim, On the theory of identical particles. Il Nuovo Cimento B. 37 (1) (1977): 1-23. [31] L. Lewin, Polylogarithms and associated functions. With a foreword by A. J. Van der Poorten. North-Holland Publishing Co., New York-Amsterdam, 1981. [32] J. Naudts, Deformed exponentials and logarithms in generalized thermostatistics. Phys. A 316 (2002), no. 1-4, 323-334. [33] J. Naudts, Generalised Thermostatistics. Springer Verlag, 2011. [34] S. Roman, The umbral calculus. Pure and Applied Mathematics, 111. Aca- demic Press, Inc., New York, 1984. [35] S.M. Roman, G.-C. Rota, The umbral calculus. Advances in Math. 27 (1978), no. 2, 95-188. [36] G.-C. Rota, D. Kahaner, A. Odlyzko, On the foundations of combinatorial theory. VIII. Finite operator calculus. J. Math. Anal. Appl. 42 (1973), 684C760. [37] N. J. A. Sloane, The On-Line Encyclopedia of Integer Sequences, http://oeis.org, 2013. [38] L.A. Takhtajan, Quantum mechanics for mathematicians. Graduate Studies in Mathematics, 95. AMS, Providence, RI, 2008. [39] I. J. Taneja, On generalized information measures and their applications, Adv. Electronics and Electron Phys. Vol. 76 (1989), 327-413. [40] J. Taylor, Formal group laws and hypergraph colorings, PhD thesis, Univer- sity of Washington, 2016. MATHEMATICS RELATED TO THE INTERPOLATING STATISTICS 81

[41] C. Tsallis, Possible generalization of Boltzmann-Gibbs statistics. J. Statist. Phys. 52 (1988), no. 1-2, 479-487. [42] C. Tsallis, What are numbers that experiments provide?, Quimica Nova 17 (6) 1994, 468-471. [43] C. Tsallis, Nonextensive statistical mechanics and thermodynamics: histor- ical background and present status. Nonextensive statistical mechanics and its applications (Okazaki, 1999), 3-98, Lecture Notes in Phys., 560, Springer, Berlin, 2001. [44] C. Tsallis, Introduction to nonextensive statistical mechanics. Springer Ver- lag, 2009. [45] D. C., Tsui, H. L.Stormer, A. C. Gossard, Two-dimensional magnetotransport in the extreme quantum limit. Phys. Rev. Lett. 48 (22) (1982): 1559-1562. [46] F. Wilczek, Quantum mechanics of fractional-spin particles. Physical Review Letters. 49 (14) (1982): 957C959. [47] X.-G. Wen, Topological orders in rigid states, Int. J. Mod. Phys. B. 4 (2) (1989): 239. [48] D. Wright, The tree formulas for reversion of power series. J. Pure Appl. Algebra 57 (1989), no. 2, 191-211. [49] Y.-S. Wu, Statistical distribution for generalized ideal gas of fractional- statistics particles, Phys. Rev. Lett. 73 (1994), 922-925. [50] S.-L. Yang, M.-Y. Jiang, The m-Schr¨oder paths and m-Schr¨oder numbers, Disc. Math. (2021) Vol. 344, Issue 2, 112209. [51] J. Zhou, Local mirror symmetry for one-legged topological vertex, arXiv:0910.4320 (2009).

Department of Mathematical Sciences, Tsinghua University, Beijng, 100084, China Email address: [email protected]