Arxiv:2106.12935V1 [Math.CO]

arXiv:2106.12935v1 [math.CO] 24 Jun 2021 and h elnumbers Bell The n hystsytercrec formula (2) recurrence the satisfy they and h hsclltrtr,wihwsrdsoee yKatriel by operator rediscovered creation was the which literature, physical the (3) formula recursive the satisfy they and where xaso of expansion oevr h trignmeso h eodkind second the of numbers Stirling the Moreover, pvyrlto;Dbnk formula. Dobinski relation; Spivey aa (1) atto e fsize of set a partition reigof ordering † e od n phrases. and words Key 2010 h trignmeso h eodkn,dntdby denoted kind, second the of numbers Stirling The − dx d a X ahmtc ujc Classification. Subject Mathematics † oyoil n trignmes h otqatmanalog post-quantum the numbers, Stirling and polynomials Abstract. sepudd lcdtn ( a ( elucidating the expounded, for is relation connection recurrence The A numbers. established. Stirling and polynomials Touchard stedrvtv prtrwt epc to respect with operator derivative the is a eoe h prtro utpiainwt h aibe t variable, the with multiplication of operator the denotes nteBsnFc pc 1 27]. [1, space Fock Boson the in 1 = ( ,q p, X ) dx AAOUSO H EEAIE TOUCHARD GENERALIZED THE OF -ANALOGUES d nti ae eitouea( a introduce we paper this In n OYOIL N TRIGNUMBERS STIRLING AND POLYNOMIALS B sfollows as n a n r nw ytefloigrelation following the by known are n h niiainoperator annihilation the and ocadplnma;( polynomial; Touchard into S ( k ,k n, X o-mt ust.Te pera offiinsi h normal the in coefficients as appear They subsets. non-empty ,q p, x = ) dx d n -eomto fSie’ relation. Spivey’s of )-deformation = B S 1. n B 17,0A0 16,1S5 05A99. 11S05, 11B65, 05A30, 11B73, n X k ( ACNOUSSI LAHCEN +1 = n =0 n n Introduction − = X k S = ,q p, =0 n 1 ( X k ,q p, ,k n, ,q p, k , X k =0 n -eomdaaouso h eeaie Touchard generalized the of analogues )-deformed S =0 n 1 -eomdgnrlzdTuhr polynomials Touchard generalized )-deformed ( -acls ( )-calculus; − S ,k n, ) ( k j Y )+ 1) ,k n, x n k =0 − h eain()hsde oncinin connection deep has (1) relation The . ) S 1 X ( B ( x ) ,k n, a k , S k − † . S ( n aifigtecmuainrelation commutation the satisfying ( j e fthe of ues ,q p, ,apa loa offiinsi the in coefficients as also appear ), dx ,k n, d ) 2] aey ti once with connected is it Namely, [20]. − , ewe hs eomtosis deformations these between -trignme;( number; )-Stirling 1 ,cuttenme fwy to ways of number the count ), k k , , ) . q dfre generalized -deformed a is, hat Xf ,q p, ( -elnumber; )-Bell x = ) xf ( x ), 2 LAHCEN OUSSI

In 2008, Spivey [45] generalized the relation (3) as follows m n m (4) B = km−jS(n, k) B . n+m j j j k X=0 X=0 The above formula is known in the literature as “Spivey’s Bell number formula”. Subsequently, in the same year, Katriel [21] derived a q-deformed formula of (4): m n m (5) B (q)= [k]m−jqjkS (n, k) B (q), n+m q q j j j X=0 Xk=0 where k−1 (6) Sq(n, k)= q S(n 1, k 1) + [k]qSq(n 1, k), − − − are the q-Stirling numbers of the second kind with the initial value Sq(0, 0) = 1, and [k]q := 1−qk 1−q denotes the q-number [17, 18, 19]. Consequently, the associated q-Bell numbers are given by n Bn(q)= Sq(n, k). Xk=0 Recently, the present author [38] derived a (p,q)-deformed variant of Eq. (4). The resulting relation can be expressed as follows n m k m−j jk m n+m−j (7) B (p,q;1) = p(2)[k] q S (n, k) B (p,q; p ), n+m p,q p,q j j j Xk=0 X=0 where n−k k−1 (8) Sp,q(n, k)= p q S(n 1, k 1) + [k]p,qSp,q(n 1, k), − − − pk−qk are the (p,q)-Stirling numbers of the second kind, with Sp,q(0, 0) = 1, [k]p,q := p−q denotes n k (2) the (p,q)-number, which reduces to [k]q for p = 1, and Bn(p,q;1) = k=0 p Sp,q(n, k) are the corresponding (p,q)-Bell numbers in analogy to the classical and deformed cases. Obviously, for p = 1, Eqs. (7) and (8) reduce to (5) and (6), respectively.P If further q = 1, they reduce to (4) and (2), respectively. It is worthwhile to mention that various generalizations of Stirling and Bell numbers have been extended in several ways by several authors, see for instance [16, 3, 6, 11, 26, 25, 37, 44, 33, 34]. In particular, Mansour et al. considered in [33] the noncommutative operators U and V satisfying the commutation relation UV = VU + hV s (for s R and h C 0 ), and ∈ ∈ \ { } introduced a new family of generalized Stirling numbers Ss;h(n, k) by n n s(n−k)+k k (9) (VU) = Ss;h(n, k)V U . k X=1 Moreover, they introduced in [34] the q-version of the above generalized Stirling numbers, which appear in the normal ordering of (VU)n of the noncommutative operators U and V satisfy the q-commutation relation UV = qV U + hV s. Namely n n s(n−k)+k k (10) (VU) = Ss h(n, k q)V U , ; | Xk=1 (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 3 and s(n+1−k)+k−1 (11) Ss h(n + 1, k q)= q Ss h(n, k 1 q)+ h[s(n k)+ k]qSs h(n, k q) ; | ; − | − ; | for all n 0, k 1, with the initial values Ss;h(n, 0 q) = δn,0 and Ss;h(0, k q) = δ0,k for all n, k N ≥:= 0, 1≥, 2,... , with the corresponding generalized| Bell numbers | ∈ 0 { } n (12) Bs h(n q) := Ss h(n, k q). ; | ; | k X=0 In 2015, Mansour and Schork [31] derived the following q-deformed generalization of Spivey’s relation (4) (13) n l m n B m−1 n l q ;[m]q ( + )= m | k m−1 j q ,hm,j+(m−1)l(q) Xk=0 X=0 n k m k j m l m − S − ( +( −1) )B − [j + (m 1)l]q m 1 ;[m] (l, j q )q m 1 ;[m] (k q ) . × − m q | m q | s [m−1]q n n k−1 where hm,s(q) := q and := (1 + h[i]q) is the (q, h)-binomial coeffi- [s]q k k i=0 (q,h) q n Q cients [2], with := [n]q! denotes the q-Gauss binomial coefficients. k [k]q![n−k]q! q In the present work (Section 3), we derive the (p,q)-analogues of the generalized Stirling numbers and Spivey’s relation above. The Touchard polynomials (also known as exponential polynomials or Bell polynomials) have many applications and play important rule in areas of probability and statistics, see e.g., [4, 7, 23, 24, 35, 32]. In particular, the Touchard polynomial Tn(λ) is the n-th moment of a random variable Y which has Poisson distribution and expected value λ. In an operational framework, the Touchard polynomials can be defined [4, 39, 40, 5, 9] by the Rodrigues-like formula: n −x d x (14) Tn(x)= e X e for n N . dx ∈ 0 In 2012, Dattoli et al. [9] introduced a generalization of the Touchard polynomials (14) for higher order m N as follows: ∈ 0 d n (15) T (m)(x)= e−x Xm ex, for n N , n dx ∈ 0 which reduce to the classical Touchard polynomials (14) for m = 1. Moreover, they established the following recursive formula: d (16) T (m)(x)= xm 1+ T (m)(x). n+1 dx n Mansour et al. extended the generalized Touchard polynomials for negative integer order m Z in [34], and furthermore, for arbitrary real order m R in [30]. On the other hand, Mansour∈ and Schork [31] introduced a q-deformation of the generalized∈ Touchard polynomials for m R by the following ∈ (m) m n (17) T (x)= Eq( x)(X Dq) eq(x), n;q − 4 LAHCEN OUSSI

d where Dq denotes the analogue of the derivative operator dx (also called Jackson derivative), 2 (n) n n given by D f(x) = f(x)−f(qx) , E (x) = q x and e (x) = x are the q- q 1−q q n≥0 [n]q! q n≥0 [n]q! exponential functions, and [n] != n [j] is the q-factorial function of the q-number [n] . q j=1 qP P q The q-deformed generalized Touchard polynomials satisfy the recursion formula [31]: Q (m) m (m) (18) T (x)= x 1+ Eq( x)eq(qx)Dq T (x), n;q − n;q which reduces to (16) for q=1. Using the normal ordering (1) and the fact that (see [31])

n m n n m m k k ( −1) S − (X Dq) = X m 1 ;[m] (n, k q )X Dq , m q | Xk=0 one obtains directly from (14) and (17) the following relationships between the Stirling numbers S(n, k) (resp. q-generalized Stirling numbers Ss;h(n, k q)) and the Touchard polynomials (m) | Tn(x) (resp. q-generalized Touchard polynomials Tn;q (x)) :

n k (19) Tn(x)= S(n, k)x =: Bn(x), k X=0 and n m n m m k n m ( ) ( −1) S − ( −1)B − (20) Tn;q (x)= x m 1 ;[m] (n, k q )x = x m 1 ;[m] |n;qm (x), m q | m q k X=0 B n S k where Bn(x) and s;h|n;q(x) := k=0 s;h(n, k q)x are the associated Bell polynomials and q-generalized Bell polynomials, respectively. | P (m) Noting that for q = 1, the undeformed generalized Touchard polynomials Tn (x) satisfy

n (m) n(m−1) k Tn (x)= x Sm,1(n.k)x , Xk=0

S − where, Sm,1(n, k) = m 1 ;m(n, k 1) are another version of the generalized Stirling numbers m | of the second kind considered for instance by Lang [26]. The structure of the paper is outlined as follows: In Section 2 we recall and ﬁx some elementary notions and terminology about (p,q)-calculus, which are useful for the rest of the paper. Section 3, contains the main result and is divided into two subsections, in Subsection 3.1, the (p,q)-deformed generalized Stirling numbers of the second kind are introduced, the case p = 1, corresponds to the q-version. In Subsection 3.2, we discuss the (p,q)-deformed generalized Touchard polynomials for arbitrary real order, with its relation to the (p,q)-deformed generalized Stirling numbers of the second kind. A recurrence relation for these polynomials and some other properties are derived as well as the (p,q)-deformed generalized Spivey formula and the (p,q)-Dobi´nski formula. Finally, in Section 4 we outline some remarks and mention some possible open questions related to this study for future research. (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 5

2. Preliminaries In this section we review some basic tools and terminology about (p,q)-calculus which we will use in the subsequent sections. For more details about (p,q)-calculus and related aspects, we refer the reader to [43, 8, 13, 14]. For x R with p = q, the (p,q)-number (or twin-basic number) [x]p,q is deﬁned as ∈ 6 px qx (21) [x] := − . p,q p q − If n is non-negative integer, then

n n n p −q = pn−kqk−1, if n N := 1, 2, 3,... ; [n] := p−q ∈ { } p,q  k=1 0, P if n = 0 . The (p,q)-factorial is given by  n [k]p,q, if n N; [n] ! := ∈ p,q k=1 1Q, if n = 0 , and the (p,q)-Gaussian binomial coefficient is defined by [n] ! n p,q , if 0 k n; := [n−k]p,q![k]p,q! ≤ ≤ k p,q (0, if k>n 0 or k < 0 . ≥ Moreover, we define the following operators:

(1) The (p,q)-derivative operator Dp,q: f(px)−f(qx) x(p−q) , if x = 0; Dp,qf(x)= 6 (f ′(0), if x = 0, for functions f which are diﬀerentiable at x = 0. (2) The operator of multiplication X with the variable: Xf(x)= xf(x)

(3) The Fibonacci operator Np [36]:

Npf(x)= f(px). k k−1 It is clear that Dp,qx = [k]p,qx , and the (p,q)-Leibniz formula reads n n n n−k k k n−k (22) Dp,q(f(x)g(x)) = Dp,q (f)(p x)Dp,q(g)(q x), k p,q Xk=0 n n−1 0 where, Dp,qf(x) := Dp,q(Dp,q f(x)) for n = 1, 2, 3,..., and Dp,q denotes the identity operator. The operators Dp,q, X and Np satisfy the following (p,q)-commutation relations:

(23) Dp,qX = qXDp,q + Np,

(24) NpX = pXNp, and

(25) Dp,qNp = pNpDp,q. 6 LAHCEN OUSSI

The (p,q)-deformed exponential functions Ep,q(x) and ep,q(x) can be deﬁned, respectively, as follows: n q(2) (26) E (x)= xn, p,q [n] ! n p,q X≥0 and (n) p 2 n (27) ep,q(x)= x , [n] ! n p,q X≥0 which satisfy the basic identity

(28) ep,q(x)Ep,q( x) = 1. − The nth-(p,q)-derivative of the exponential function ep,q(x) satisﬁes the diﬀerential equation [42]: n (n) n (29) Dp,qep,q(x)= p 2 ep,q(p x). The q-settings (which correspond to p = 1) of the above notions are the ones considered in [12, 38, 31, 17, 19]. We conclude this section by recalling the following proposition from [38] which will be useful in the subsequent computations. Proposition 2.1 ([38]). Let n be a nonnegative integer. Then the following relation holds true: n n k n−k k (30) (XDp,q) = Sp,q(n, k)X Np Dp,q, Xk=0 where n−k k−1 (31) Sp,q(n, k)= p q Sp,q(n 1, k 1) + [k]p,qSp,q(n 1, k), − − − are the (p,q)-Stirling numbers of the second kind, with the initial value Sp,q(0, 0) = 1, and consequently the numbers n (32) Bn(p,q)= Sp,q(n, k), Xk=0 can be considered as the associated (p,q)-Bell numbers.

3. Main result 3.1. The (p,q)-deformed generalized Stirling numbers of the second kind. In this subsection, we introduce a (p,q)-analogue of the generalized Stirling numbers of the second kind which were introduced by Mansour et al. in [33]. For h C 0 and s R, let U, V and Wp be noncommutative operators satisfy by the following∈ (p,q\ {)-commutation} ∈ relations: s (33) UV qV U = hV Wp, − (34) WpV = pV Wp,

(35) UWp = pWpU,

(36) W1 = I ( identity operator). (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 7

The case p = 1 and s = 0, reduces to the q-commutation relation UV qV U = h of the q-deformed Heisenberg-Weyl algebra. Furthermore, q = 1 yields the usual− Heisenberg-Weyl algebra ([U, V ] := UV VU = h). On the other hand, choosing p = q = s = 1, one obtains the commutation relation− UV VU = hV of the two-dimensional non-Abelian Lie algebra [41, 46, 47], and p = q = 1 and−s = 2 yield the Jordan quantum plane [10, 28]. Proposition 3.1. For all k N, we have the following identity: ∈ k k k s+k−1 (37) UV q V U = h[k]p,qV Wp, − where U, V and Wp are noncommutative operators satisfying Eqs. (33)-(36). Proof. The proof will be by induction on k. For k = 1, Eq. (37) holds, since Eq. (33). Assume that Eq. (37) holds for k and let us prove it for k + 1. Using induction hypothesis and Eqs. (33) and (34), we obtain k+1 k k s+k−1 UV = q V UV + h[k]p,qV WpV k k s s+k−1 = q V (qV U + hV Wp)+ h[k]p,qV WpV k+1 k+1 k k+s s+k−1 = q V U + hq V Wp + h[k]p,qV WpV k+1 k+1 k k+s s+k = q V U + hq V Wp + p[k]p,qV Wp k+1 k+1 k s+k = q V U + h(q + p[k]p,q)V Wp k+1 k+1 s+k = q V U + h[k + 1]p,qV Wp, which completes the proof.

Corollary 3.2. One can reformulate Proposition 3.1 in terms of the operators Dp,q, X and Np as follows s s s s−1 (38) Dp,qX = q X Dp,q + [s]p,qX Np. Definition 3.3. Let h C 0 and s R. We define the (p,q)-deformed generalized Stirling ∈ \{ } ∈ n numbers of the second kind Ss;h(n, k p,q), as normal ordering coefficients of (VU) by the following: | n n s(n−k)+k n−k k (39) (VU) = Ss h(n, k p,q)V W U , ; | p Xk=0 with U, V and Wp satisfy Eqs. (33)-(36). Moreover, the corresponding (p,q)-deformed generalized Bell polynomials and Bell numbers are defined, respectively, by n k (40) B (x) := Ss h(n, k p,q)x , s;h|n;p,q ; | Xk=0 and n (41) Bs h(n p,q) := Ss h(n, k p,q). ; | ; | Xk=0 Examples 3.4. From Definition 3.3 and using Proposition 3.1 with the commutation relation (35) we can determine directly the first few instances of the (p,q)-deformed generalized Stirling numbers as follows: 8 LAHCEN OUSSI n = 1: Ss h(1, 1 p,q) = 1 and Bs h(1 p,q) = 1. ; | ; | n = 2: Ss;h(2, 1 p,q)= h, Ss;h(2, 2 p,q)= q, and consequently, Bs;h(2 p,q)= h + q. | 2 | s | 3 n = 3: Ss;h(3, 1 p,q) = h [s + 1]p,q, Ss;h(3, 2 p,q) = hq ([2]p,q + pq ) , Ss;h(3, 3 p,q) = q , and consequently,| the associated Bell numbers| | 3 s 2 Bs h(3 p,q)= q + hq([2]p,q + q p)+ h [s + 1]p,q. ; | Theorem 3.5. The (p,q)-deformed generalized Stirling numbers Ss;h(n, k p,q) satisfy the following recursive formula: | (42) n−k+1 s(n−k+1)+k−1 Ss h(n + 1, k p,q)= p q Ss h(n, k 1 p,q)+ h[s(n k)+ k]p,qSs h(n, k p,q), ; | ; − | − ; | for all n 0 and k 1, with the initial conditions Ss h(n, 0 p,q) = δn, for all n 0, and ≥ ≥ ; | 0 ≥ Ss h(0, k p,q)= δ ,k for all k 0. ; | 0 ≥ Proof. By Definition 3.3, one has n+1 n+1 s(n−k+1)+k n−k+1 k (43) (VU) = Ss h(n + 1, k p,q)V W U . ; | p Xk=0 On the other hand, (VU)n+1 = (VU)(VU)n n s(n−k)+k n−k k = VU Ss h(n, k p,q)V W U ; | p k=0 n X s(n−k)+k n−k k = Ss h(n, k p,q)V (UV )W U ; | p k=0 Xn s(n−k)+k s(n−k)+k = Ss h(n, k p,q)V q V U + h[s(n k)+ k]p,q ; | − k X=0 s(n−k)+k+s−1 n−k k V Wp W U × × p n s(n−k)+k s(n−k)+k+1 n−k k = Ss h(n, k p,q) q V UW U + h[s(n k)+ k]p,q ; | p − k X=0 V s(n−k)+k+sW n−k+1U k × p n n−k s(n−k)+k s(n−k)+k+1 n−k k+1 = Ss h(n, k p,q) p q V W U + h[s(n k)+ k]p,q ; | p − k X=0 V s(n−k+1)+kW n−k+1U k , × p where in the fourth equality we use Eq. (37) and in the last equality we use Eq. (35). Comparing the coefficients with Eq. (43) finishes the proof.

Corollary 3.6. The (p,q)-deformed generalized Stirling numbers Ss;h(n, k p,q) satisfy the following identity: | n−k Ss h(n, k p,q)= h Ss (n, k p,q). ; | ;1 | Proof. The proof follows directly from (42). (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 9

Remark 3.7. For s = 0 and h = 1, the recurrence formula (42) becomes n−k+1 k−1 S (n + 1, k p,q)= p q S (n, k 1 p,q) + [k]p,qS (n, k p,q), 0;1 | 0;1 − | 0;1 | which is the recurrence formula for the (p,q)-deformed Stirling numbers of the second kind Sp,q(n, k), introduced recently by the present author in [38]. Therefore,

S (n, k p,q)= Sp,q(n, k). 0;1 | Moreover, for p = 1 we obtain the well known recursive formula for the q-Stirling numbers of the second kind (see [34]). In addition, for p = q = 1, the recurrence formula (42) reduces to the one obtained by Mansour et al. for the generalized Stirling numbers of the second kind [33, 29]. 3.2. The (p,q)-deformed generalized Touchard polynomials. In this subsection we establish a (p,q)-version of the generalized Touchard polynomials introduced in [30], and some related properties.

Definition 3.8. For m R (n N0), we define the nth (p,q)-deformed generalized Touchard polynomials for order m∈by the∈ formula (m) n m n (44) T (x)= Ep,q( p x)(X Dp,q) ep,q(x). n;p,q − Remark 3.9. Obviously, if p = 1, Eq. (44) reduces to the q-generalized Touchard polynomials of order m (17). If further q = p = 1, we obtain the classical generalized Touchard polynomials (14). Examples 3.10. For m = 0: (0) n (n) T (x)= Ep,q( x)D ep,q(x)= p 2 for all n N . n,p;q − p,q ∈ 0 0 N Note that this is different from q-deformed case, where Tn;q(x) = 1, for all n 0. For m = 1: ∈ (1) n n T (x)= Ep,q( p x)(XDp,q) ep,q(x) n;p,q − n n k n−k k = Ep,q( p x) Sp,q(n, k)X N D ep,q(x) − p p,q k=0 Xn n (k) k n−k k = Ep,q( p x) p 2 Sp,q(n, k)X N ep,q(p x) − p k=0 Xn n k n = Ep,q( p x) S˜p,k(n, k)x ep,q(p x) − k=0 n X k = S˜p,q(n, k)x Xk=1 = Bñ(p,q; x), (k) where in the second equality we use Proposition 2.1, and S˜p,q(n, k) = p 2 Sp,q(n, k) are the other version of (p,q)-deformed Stirling numbers of the second kind introduced in [38], and Bñ(p,q; x) are the related (p,q)-Bell polynomials. 10 LAHCEN OUSSI

Proposition 3.11. Let m R 0 . Then, ∈ \ { } m m m m−1 (45) Dp,qX = q X Dp,q + [m]p,qX Np. Moreover, n m n n m m k n k k ( −1) S − − (46) (X Dp,q) = X m 1 ;[m] (n, k p,q )X Np Dp,q m p,q | Xk=0 Proof. To proof (45), we use the (p,q)-Leibniz formula (22): m m Dp,q(X f(x)) = Dp,q(x f(x)) m−1 m = [m]p,qx f(px) + (qx) Dp,qf(x) m−1 m = [m]p,qx Npf(x) + (qx) Dp,qf(x) m−1 m m = [m]p,qX Npf(x)+ q X Dp,qf(x). Hence, m m m m−1 (Dp,qX q X Dp,q) f(x) = [m]p,qX Npf(x) − m The Eq. (46) follows from Definition 3.3 of Ss h(n, k p,q), with the association U X ,V ; | 7→ 7→ Dp,q, Wp Np, and using Eq. (45). 7→ Remark 3.12. Since S (n, k p,q)= Sp,q(n, k), the formula (46) reduces to (30) for m = 1. 0;1 | Theorem 3.13. For m R 0 , the (p,q)-deformed generalized Touchard polynomials of order m are given by the∈ following\ { } formula: n m n m m k n m ( ) ( −1) S˜ − ( −1)B˜ − (47) Tn;p,q(x)= x m 1 ;[m] (n, k p,q )x = x m 1 ;[m] |n;p,qm (x), m p,q | m p,q k X=0 (k) where S˜ s;h(n, k p,q) = p 2 Ss;h(n, k p,q) is one possible definition of the (p,q)-deformed generalized Stirling| numbers with the| associated (p,q)-deformed generalized Bell polynomi- n ˜ ˜ k als Bs;h|n;p,q(x) := Ss;h(n, k p,q)x . k=0 | P Proof. The proof follows from Definition 3.8 and Eq. (46). In the next proposition, we establish a (p,q)-analogue of the recurrence relations (16) and (18) satisfied by the classical Touchard polynomials and by the q-deformed generalized Touchard polynomials, respectively. (m) Proposition 3.14. The (p,q)-deformed generalized Touchard polynomials Tn;pq(x) satisfy the recursive formula:

(m) m n+1 n (m) (48) T (x)= x Np + Ep,q( p x)ep,q(p qx)Dp,q T (x). n+1;p,q − n;p,q Proof. According to Deﬁnition 3.8, we have

(m) n+1 m m n T (x)= Ep,q( p x)(X Dp,q)(X Dp,q) ep,q(x). n+1;p,q − On the other hand, m n n (m) (X Dp,q) ep,q(x)= ep,q(p x)Tn;p,q(x). (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 11

Then,

(m) n+1 m n (m) T (x)= Ep,q( p x)X Dp,q ep,q(p x)T (x) n+1;p,q − n;p,q n+1 m n+1 (m) n (m) = Ep,q(p x)X ep,q(p x)Tn;p,q(px)+ ep,q(p qx)Dp,qTn;p,q(x) m (m) m n+1 n (m) = x T (px)+ x Ep,q( p x)ep,q(p qx)Dp,qT (x) n;p,q − n;p,q m n+1 n (m) = x Np + Ep,q( p x)ep,q(p qx)Dp,q T (x), − n;p,q where in the second equality, we use the (p,q)-Leibniz rule (22).

Remark 3.15. Clearly, for p = 1, the recurrence (48) reduces to (18), and for p = q = 1, it reduces to (16).

Now we are in the position to derive a (p,q)-analogue of the generalized Spivey’s formula.

Theorem 3.16. Let m R 0 . Then, we have the following (p,q)-deformed generalized Spivey’s relation: ∈ \ { }

(49) n l n (m−1)((n−k)(1+k)+kl) B˜ m−1 m p ;[m]p,q|n+l;p,q (1) = m n k m−1 m−1 j (p ,q ,hm,j+(m−1)l(p,q)) Xk=0 X=0 − k(j+l(m−1)) m n+l−k q S˜ m−1 l, j p,q B˜ m−1 m p ;[m]p,q ( ) ;[m]p,q|k;p,q ( ) × m | ! m Let us ﬁrst prove the following two lemmas which will be the main ingredients for the proof of Theorem 3.16.

Lemma 3.17. For any positive integer n,s R and m R 0 , we have ∈ ∈ \ { } m n s s m−1 s m n (50) (X Dp,q) X = X [s]p,qX Np + q (X Dp,q) n n n−k s m−1 n−k ks m k (51) = [s]p,q X (X Np) q (X Dp,q) , n k m− m− (p 1,q 1,hm,s(p,q)) Xk=0 − s q [m−1]p,q where h (p,q) := m− . m,s [s]p,qp 1 Proof. By virtue of Corollary 3.2, one obtains

m s s m−1 s m (X Dp,q)X = X [s]p,qX Np + q (X Dp,q) ,

and a simple induction gives the ﬁrst Eq. (50). For the second Eq. (51) we can use the (q, h)-binomial formula due to Benaoum [2]. Let us recall the (q, h)-binomial formula from [2]; for two variables R and S which satisfy SR qRS = hR2 , one has − n n (R + S)n = RkSn−k, n k (q,h) Xk=0 − 12 LAHCEN OUSSI

n n where := k−1(1 + h[j] ) is the (q, h)-binomial coeﬃcients. n k n k j=0 q − (q,h) − q m−1 s m Consider S = [s]p,qX Np andQR = q (X Dp,q), then

s m m−1 (RS)f(x)= q [s]p,qX Dp,qX Npf(x) s m m−1 = q [s]p,qX Dp,qX f(px) s m m−2 2 s m−1 m m−1 = q [s]p,q[m 1]p,qX X f(p x)+ q [s]p,qq X X Dp,qf(px) − s 2m−2 2 s m−1 m m−1 = q [s]p,q[m 1]p,qX N f(x)+ q [s]p,qq pX X NpDp,qf(x) − p s m−1 q [m 1]p,q q = − S2f(x)+ (SR)f(x). [s] pm−1 p p,q Therefore,

q m−1 RS = h (p,q)S2 + SR. m,s p n Let us introduce the (p,q,h)-binomial coeﬃcients as follows: n k − (p,q,h)

k−1 n n k(k−1) n 1−n := = p (1 + hp [j]p,q). n k n k q n k (p,q,h) ( ,h) p,q j − − p − Y=0

n n Clearly, for p = 1, = . n k n k − (1,q,h) − (q,h) Thus,

n n m−1 s m n n−k m−1 n−k [s]p,qX Np + q (X Dp,q) = [s]p,q (X Np) n k m− m− (p 1,q 1,hm,s(p,q)) Xk=0 − sk m k q (X Dp,q) , × and this proves Eq. (51).

Lemma 3.18. The (p,q)-generalized Touchard polynomials of order m R 0 satisfy the recursive formula ∈ \ { }

(52) n l n−k (m) n k(j+(m−1)l) (m−k)(m−1) Tn+l,p,q(x)= q p j + l(m 1) n k m−1 m−1 − p,q j (p ,q ,hm,s(p,q)) Xk=0 X=0 − h i m j m l n k (m) n k l S˜ − +( −1)( + − ) − + m 1 ;[m] (l, j p,q )x Tk,p,q(p x). × m p,q | (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 13

Proof. We have

m n+l m n m l (X Dp,q) = (X Dp,q) (X Dp,q) l m n l m m j l j j ( −1) S − − = (X Dp,q) X m 1 ;[m] (l, j p,q )X Np Dp,q m p,q | j X=0 n l n k(j+l(m−1)) = [j + l(m 1)]p,qq n k m−1 m−1 − j (p ,q ,hm,s(p,q)) Xk=0 X=0 − m j l m m n k m k l j j S − + ( −1) −1 − − m 1 ;[m] (l, j p,q )X (X Np) (X Dp,q) Np Dp,q × m p,q | n l n (n−k)(m−1) k(j+l(m−1)) = p [j + l(m 1)]p,qq n k m−1 m−1 − k j (p ,q ,hm,s(p,q)) X=0 X=0 − m j l n k m n k m k l j j S − +( + − )( −1) − − m 1 ;[m] (l, j p,q )X Np (X Dp,q) Np Dp,q, × m p,q | where in the second equality, we use Proposition 3.11, in the third equality, we use Lemma 3.17 with s = j + l(m 1), and in the last equality we use (24). − Applying the above identity to the (p,q)-exponential function ep,q(x), yields

n l m l+n n (n−k)(m−1)+(j) m X D e x p 2 S m−1 l, j p,q ( p,q) p,q( )= ;[m]p,q ( ) n k m−1 m−1 m | j (p ,q ,hm,s(p,q)) Xk=0 X=0 − k(j+l(m−1)) j+(l+n−k)(m−1) n−k m k l [j + l(m 1)]p,qq X N (X Dp,q) ep,q(p x) × − p n l n (n−k)(m−1)+(j) m p 2 S m−1 l, j p,q = ;[m]p,q ( ) n k m−1 m−1 m | k j (p ,q ,hm,s(p,q)) X=0 X=0 − T (m) (plx) j l m qk(j+l(m−1))Xj+(l+n−k)(m−1)N n−k k,p,q [ + ( 1)]p,q p k+l × − Ep,q( p x) − n l n (n−k)(m−1)+(j) m p 2 S m−1 l, j p,q = ;[m]p,q ( ) n k m−1 m−1 m | k j (p ,q ,hm,s(p,q)) X=0 X=0 − T (m) (pn−k+lx) j l m qk(j+l(m−1))xj+(l+n−k)(m−1) k,p,q . [ + ( 1)]p,q n+l × − Ep,q( p x) − n+l Multiplying both sides of the above identity by Ep,q( p x) gives −

n l (m) n (n−k)(m−1) k(j+l(m−1)) Tn+l,p,q(x)= p [j + l(m 1)]p,qq n k m−1 m−1 − j (p ,q ,hm,s(p,q)) Xk=0 X=0 − m j l n k m (m) n k l S˜ − +( + − )( −1) − + m 1 ;[m] (l, j p,q )x Tk,p,q(p x). × m p,q | 14 LAHCEN OUSSI

Proof of Theorem 3.16. By virtue of Lemma 3.18 and Theorem 3.13, one obtains the recurrence relation for (p,q)- deformed generalized Bell polynomials as follows

n l n (m−1)((n−k)(1+k)+kl) B˜ m−1 m x p (53) ;[m]p,q|n+l;p,q ( )= m n k m− m− k j (p 1,q 1,hm,s(p,q)) X=0 X=0 − k j l m m j n l k ( + ( −1))S˜ − B˜ − + − [j + l(m 1)]p,qq m 1 ;[m] (l, j p,q )x m 1 ;[m] |k;p,qm(p x). × − m p,q | m p,q Plugging x = 1 into this identity gives the assertion.

Remark 3.19. If m = 1, the recurrence relation (49) reduces to the (p,q)-deformed Spivey’s formula obtained by the present author [38], and if p = 1, one obtains the recurrence relation for q-deformed generalized Spivey’s formula given by Mansour and Schork [31].

In the next result we present an analogue of Dobi´nski’s formula for (p,q)-deformed gener- ˜ alized Bell polynomials Bs;h|n;p,q(x).

Proposition 3.20 (Dobi´nski’s formula). The (p,q)-deformed generalized Bell polynomials

B˜ m−1 m (x) satisfy the following relation: m ;[m]p,q|n;p,q

n−1 [k + j(m 1)]p,q − 1 j=0 ! (k) k (54) B˜ m−1 m (x)= p 2 x , ;[m]p,q|n;p,q n Q m ep,q(p x) [k]p,q! k X≥0

B˜ − m and consequently, the (p,q)-deformed generalized Bell numbers m 1 ;[m] (n p,q ) satisfy m p,q |

n−1 [k + j(m 1)]p,q − m 1 j=0 ! (k) (55) B˜ m−1 (n p,q )= p 2 . ;[m]p,q n Q m | ep,q(p ) [k]p,q! Xk≥0 n−1 m n k k+n(m−1) Proof. Since (X Dp,q) x = [k + j(m 1)]p,q x , we have j=0 − ! Q n−1 [k + j(m 1)]p,q − n m n n j=0 ! (k) k+n(m−1) Ep,q( p x) (X Dp,q) ep,q(x)= Ep,q( p x) Q p 2 x . − − [k]p,q! Xk≥0 Then, by Deﬁnition 3.8

n−1 [k + j(m 1)]p,q − (m) n j=0 ! (k) k+n(m−1) Tn,p,q(x)= Ep,q( p x) Q p 2 x . − [k]p,q! k X≥0 (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 15

Hence, Theorem 3.13 yields

n−1 [k + j(m 1)]p,q − 1 j=0 ! (k) k B˜ m−1 m (x)= p 2 x . ;[m]p,q|n;p,q n Q m ep,q(p x) [k]p,q! Xk≥0 Setting x = 1, we arrive at (55). Remark 3.21. Note that, for m = 1, Eq. (54) reduces to the (p,q)-deformed Dobiński formula given by the present author [38], and further, for p = q = m = 1, one obtains the ordinary Dobiński formula for Bell polynomials 1 k B (x)= xk. n e k! Xk≥0 Thus, Eq. (54) can be considered as the (p,q)-deformed generalized Dobiński formula.

4. Concluding remarks (1) Let us mention that for m = 1, the (p,q)-deformed Touchard polynomials Tn;p,q(x) considered in this paper are slightly diﬀerent from those considered by Mansour et al. in [15, 22]. It would be very interesting to ﬁnd new connections of the (p,q)-deformed generalized Touchard (m) polynomials Tn;p,q(x) to well known combinatorial numbers and polynomials. In particular, 1 it was shown in [30] that the undeformed generalized Touchard polynomials of order m = 2 can be expressed in terms of Hermite polynomials

1 n ( 2 ) i Tn (x)= Hn( i√x). 2 − On the other hand, Mansour and Schork [31] showed that the q-deformed generalized Touchard (m) polynomials Tn;q (x) are connected to q-deformed Laguerre polynomials and q-deformed Bessel polynomials, respectively, for m = 2 and m = 1. Therefore, it would be interesting to determine whether− the above relations hold true for the n;p,q (p,q)-deformed generalized Touchard polynomials T(m) (x) considered in this paper. However, the (p,q)-deformed generalized Stirling numbers Ss;h(n p,q) should be established deeply, as well as several special choices for the parameters s and| h should be exhibited explicitly. These and other questions warrant further study.

5. Acknowledgment Research partially supported by the Austrian FWF Project P29355-N35 and by the Polish National Science Center (NCN) grant 2016/21/B/ST1/00628. The author would like to thank Franz Lehner for careful reading and helpful comments.

References [1] M. Arik and D. Coon, Hilbert spaces of analytic functions and generalized coherent states, J. Math. Phys. 17 (4), 524–527 (2008). [2] H. B. Benaoum, (q, h)-analogues of Newton’s binomial formula, J. Phys. A: Math. Gen.32 (10) 2037–2040 (1999). [3] P. Blasiak, K.A. Penson and A.I. Solomon, The boson normal ordering problem and generalized Bell numbers, Ann. Comb. 7 127–139 (2003). 16 LAHCEN OUSSI

[4] K. N. Boyadzhiev, A series transformation formula and related polynomials, Int. J. Math. Math. Sci, 23 3849–3866 (2005). [5] K. N. Boyadzhiev, Exponential polynomials, Stirling numbers, and evaluation of some gamma integrals, Abstr. Appl. Anal. Article ID 168672 (2009) [6] L. Carlitz, On arrays of numbers, Amer. J. Math. 54 739–752 (1932). [7] L. Carlitz, Degenerate Stirling, Bernoulli and Eulerian numbers, Utilitas Mathematica 15 51–88 (1979). [8] R. B. Corcino, On P,Q-Binomial coefficients, Electron. J. Combin. Number Theory, 8 #A29 (2008). [9] G. Dattoli, B. Germano, M. R. Martinelli, P. E. Ricci, Touchard like polynomials and generalized Stirling numbers, Appl. Math. Comput. 218 (12) 6661–6665 (2012). [10] E. E. Demidov, Y. I. Manin, E. E. Mukhin and D. V. Zhdanovich Non-standard quantum deformations of GL(n) and constant solutions of the Yang-Baxter equation, Prog. Theor. Phys. 102 (Suppl.), 203 (1990). [11] B.S. El-Desouky, N.P. Cakićand T. Mansour, Modified approach to generalized Stirling numbers via differential operators, Appl. Math. Lett. 23 115–120 (2010). [12] T. Ernst, A method for q-calculus, J. Nonlinear Math. Phys. 10 (4) 487–525 (2003). [13] V. Gupta, Bernstein Durrmeyer operators based on two parameters, Facta Universitatis,Series: Mathe- matics and Informatics, 31(1), 79–95. [14] V. Gupta, T. M. Rassias, P. N. Agrawal, A. M. Acu, Basics of Post-quantum Calculus. Recent Advances in Constructive Approximation Theory. Springer Optimization and Its Applications, vol.138. Springer, Cham, (2018). [15] O. Herscovici and T. Mansour, Tsallis p, q-deformed Touchard polynomials and Stirling numbers, Russ. J. Math. Phys. 24 37–50 (2017). [16] F.T. Howard, Degenerate weighted Stirling numbers, Discrete Math. 57 45–58 (1985). [17] F. H. Jackson, On q-functions and a certain difference operator, Trans. Royal Soc. Edinburgh 46 253–281 (1909). [18] F. H. Jackson, On q-definite integrals, Pure Appl. Math. Q. 41 193–203 (1910) . [19] V. Kac and P. Cheung, Quantum calculus, New York: Springer, (2002). [20] J. Katriel, Combinatorial aspects of Boson algebra, Lett. Nuovo Cimento 10 565–567 (1947). [21] J. Katriel, On a generalized recurrence for Bell numbers, J. Integer Seq. 11 Article 08.3.8. (2008). [22] T. Kim, O. Herscovici, T. Mansour and S.-H. Rim, Differential equations for p, q-Touchard polynomials, Open Math. 14 908–912 (2016). [23] D. S. Kim and T. Kim, Some identities for Bernoulli numbers of the second kind arising from a non-linear differential equation, Bull. Korean Math. Soc. 52(6) 2001–2010 (2015) . [24] D. S. Kim, T. Kim, T. Mansour and J. -J. Seo, Linear differential equations for families of polynomials, J. Ineq. Appl. 2016, Article 95 (2016) . [25] W. Lang, Combinatorial interpretation of generalized Stirling numbers, J. Integer Seq. 12 Article 09.3.3. (2009). [26] W. Lang, On generalization of the Stirling number triangles, J. Integer Seq. 3 Article 00.2.4. (2000). [27] M. M. Mangontarum, Spivey’s Bell number formula revisited, J. Integer seq. 21 Article 18.1.1. (2018) [28] Y. I. Manin, Topics in Noncommutative Geometry, Princeton University Press, Princeton (1991). [29] T. Mansour and M. Schork, The commutation relation xy = qyx + hf(y) and Newton’s binomial formula, Ramanujan J. 25 405–445 (2011). [30] T. Mansour, M. Schork The generalized Touchard polynomials revisited, Apll. Math. Comput. 219 9978– 9991 (2013). [31] T. Mansour, M. Schork, On q-deformed generalized Touchard polynomials, Indag. Math. 26 634–645 (2015). [32] T. Mansour and M. Schork, Commutation relations, normal ordering, and Stirling numbers, CRC Press, (2015). [33] T. Mansour, M. Schork and M. Shattuck On a new family of generalized Stirling and Bell numbers, Electron. J. Combin. 18 #P77 (2011) . [34] T. Mansour, M. Schork and M. Shattuck, The Generalized Stirling and Bell Numbers Revisited, J. Integer Seq. 15, Article 12.8.3.(2012). [35] T. Mansour, S. Mulay and M. Shattuck, A general two-term recurrence and its solution, Europ. J. Com- bin.33( 20–26 2012). [36] J. C. Mason and D. C. Handscomb, Chebyshev Polynomials, Chapman and Hall/CRC, (2003). (p,q)-TOUCHARDPOLYNOMIALSANDSTIRLINGNUMBERS 17

[37] M.A. Mendez, P. Blasiak and K.A. Penson, Combinatorial approach to generalized Bell and Stirling numbers and boson normal ordering problem, J. Math. Phys. 46 Article 083511 (2005). [38] L. Oussi, A (p, q)-Deformed Recurrence for the Bell numbers, J. Integer Seq. 23, Article 20.5.2. (2020). [39] S. Roman, The Umbral Calculus, Dover (2005). [40] G.-C. Rota, Finite Operator Calculus, Academic Press (1975). [41] R. A. Sack, Taylor’s theorem for shift operators, Philos. Mag. VIII, 497–503 (1958). [42] P. N. Sadjang, On two (p, q)-analogues of the Laplace transform, J. Diﬀerence Eqn. Appl. 23 1562–1583 (2017). [43] P. N. Sadjang, On the fundamental theorem of (p, q)-calculus and some (p, q)-Taylor theorems, Results Math. 73:39 (2018). doi.org/10.1007/s00025-018-0783-z. [44] M. Schork, On the combinatorics of normal ordering bosonic operators and deformations of it, J. Phys. A: Math. Gen. 36 4651–4665 (2003). [45] M. Z. Spivey, A generalization recurrence for Bell numbers, J. Integer Seq. 11, Article 08.2.5. (2008) [46] R. M. Wilcox, Exponential operators and parameter diﬀerentiation in quantum physics, J. Math. Phys. 8, 962–982 (1967). [47] W. Witschel, Ordered operator expansions by comparison, J. Phys. A, Math. Gen. 8, 143–154 (1974)

Institute of Mathematics, University of Wroclaw, Poland. Email address: [email protected]