Search for Articles:

Contents

Laguerre-type 2nd kind hypergeometric Bernoulli polynomials and numbers

Paolo Emilio Ricci1, Pierpaolo Natalini2
1Sezione di Matematica ”Luciano Modica”, International Telematic University UniNettuno, 39 Corso Vittorio Emanuele II, I-00186 Rome, Italy
2Department of Industrial, Electronic and Mechanical Engineering, Roma Tre University, Via Vito Volterra, 62, I-00146 Rome, Italy
Copyright © Paolo Emilio Ricci, Pierpaolo Natalini. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We introduce the Laguerre-type 2nd kind hypergeometric Bernoulli polynomials and numbers. After showing their recursive computation, we exploit the Laguerre-type Blissard problem to derive a representaion formula of the relevant numbers in terms of Bell’s polynomials.

Keywords: Laguerre-type exponentials, 2nd kind hypergeometric Bernoulli polynomials and numbers, Bell’s polynomials, Blissard’s problem

1. Introduction

In a preceding article [1], we have considered the 2nd kind hypergeometric Bernoulli polynomials, and the relevant numbers, according to what follows.

Let \[\begin{aligned} \label{1} T_r(x):=\sum\limits_{\ell =0}^{r-1} \frac{x^\ell}{\ell!} \ , \end{aligned} \tag{1}\] that is the \(r\)th partial sum of the exponential series, which is polynomial of degree \(r-1\). Recalling the Stirling numbers of the second kind which can be computed as \[\begin{aligned} \label{2} S(n,k) =\frac{1}{k!} \sum\limits_{m=0}^k (-1)^{k-m} {k \choose m} m^n , \end{aligned} \tag{2}\] and the \(r\)-associate Stirling numbers of the second kind \(S(n,k;r)\), defined by \[ \left(\sum\limits_{\ell =r}^\infty \frac{x^\ell}{\ell!}\right)^{k} = \left(e^x – T_r(x)\right)^k = k! \sum\limits_{n=rk}^\infty S(n,k;r) \, \frac{x^n}{n!} , \tag{3}\] we find as a particular case: \(S(n,k;1)=S(n,k)\).

Then the 2nd kind hypergeometric Bernoulli polynomials are defined by the generating function: \[\begin{aligned} \label{4} \frac{ \frac{x^r}{r!} \ e^{tx}}{e^x – T_r(x)} = \frac{e^{tx}}{_1F_1(1,r+1,x)} =\sum\limits_{n=0}^\infty B^{[r-1,\, 1]}_n (t) \, \frac{x^n}{n!} \,, \end{aligned} \tag{4}\] and, by definition, the relevant numbers are their values at the origin: \[\begin{aligned} \label{5} B^{[r-1,\, 1]}_n := B^{[r-1,\, 1]}_n (0) \,. \end{aligned} \tag{5}\]

We recall that the study of Bernoulli numbers, together with Stirling numbers, is an important topic in Number theory because they intervene in many fields and are closely related to Riemann’s zeta-function and with Fermat’s last theorem.

Several extensions of these mathematical entities also appeared in the literature. The interested reader could look at the articles [28], where the main generalizations are recalled.

All these extensions were obtained by replacing the denominator \(e^x-1\), in the generating function of Bernoulli polynomials, with the remainder term of the truncated exponential \(e^x – T_r(x)\), where \(T_r(x)\) is defined in (1).

In [1] we have determined the exponential generating function of the polynomials \(B^{[r-1,\, 1]}_n (t)\), providing a recursive computation for the \(B^{[r-1,\, 1]}_n\) numbers and their representation in terms of Bell polynomials.

Further information can be found in the above quoted article and in the references therein, where the connection with the preceding literature by several authors [317], is reported.

In this article we exploit the Laguerre-type exponential of the first order \(e_1(x)\) defined in [18], (see also [19], where further information is given), which is used instead of the ordinary exponential for defining new sets of 2nd kind hypergeometric Bernoulli polynomials and numbers, which are called 2nd kind hypergeometric Laguerre-type Bernoulli polynomials and numbers.

We first give a recursive method for computing the 2nd kind hypergeometric Laguerre-type Bernoulli numbers, denoted as \({\cal B}^{[(r-1)^2,\, 1]}_n\). This method avoids the definition of these numbers by using partition of integers (as it is done in [17]), because the use of partitions is a computationally more expensive method. Then, by using the Laguerre-type Blissard problem, we derive a connection of these 2nd kind hypergeometric Laguerre-type Bernoulli numbers with Bell’s polynomials.

A Table of the \({\cal B}^{[(r-1)^2,\, 1]}_n\) numbers, for rational values of \(r\), derived by the second author by using the Mathematica\(^\copyright\) computer algebra program, is reported in the Appendix.

Of course a further extension is possible using the Laguerre-type exponential of order \(n\), \(e_n(x)\), also defined in [18], and a formal extension of this type is recalled in the last Remark, but the relevant numbers are not shown here, as the method to be used is essentially the same.

2. The Laguerre derivative and the relevant exponentials

The Laguerre derivative (shortly L-derivative) of order 1, is defined by [18, 19] \[\begin{aligned} \label{6} D_L:= D x D = D + x D^2, \end{aligned} \tag{6}\] where \(D=D_x = d/dx\).

An eigenfunction of the operator (6) is given by the Laguerre-type exponential (shortly L-exponential), of order 1 \[\begin{aligned} \label{7} e_1(x) = \sum\limits_{k=0}^\infty \frac{x^k}{{(k!)}^2} \,, \end{aligned} \tag{7}\] since it results \[\begin{aligned} \label{8} D_L \; e_1(ax) = a e_1(ax)\,. \end{aligned} \tag{8}\]

Furthermore, for every positive integer \(n\), the function \[\begin{aligned} \label{9} e_n(x):= \displaystyle \sum\limits_{k=0}^\infty \frac{x^k}{{(k!)}^{n+1}}, \end{aligned} \tag{9}\] is an eigenfunction of the operator \[\begin{aligned} \label{10} \begin{array}{c} D_{(n-1)L}:= D x \cdots D x D x D = D \left( xD + x^2 D^2 + \dots + x^{n-1}D^{n-1} \right) = \\ = S(n,1) D + S(n,2) x D^2 + \dots + S(n,n) x^{n-1} D^n \,, \end{array} \end{aligned} \tag{10}\] where \(S(n,k)\), \(k=1,2,\dots, n\) denote Stirling numbers of the second kind. In fact, for every constant \(a\), it results \[\begin{aligned} \label{11} D_{nL} \, e_n(ax) = a e_n(ax)\,. \end{aligned} \tag{11}\]

Note that \(e_n(x)\) reduces to the exponential function when \(n=0\), so that we put by definition: \[\begin{aligned} e_0(x):= e^x, \qquad D_{0L} := D, \qquad ({ obviously \ it \ results} \;\; D_{1L} = D_L). \end{aligned}\]

2.1. The Laguerre-type Blissard problem

Consider the formal exponential \[\begin{aligned} \label{12} e_1(a t) = \sum\limits_{k=0}^\infty \frac{a^k t^k}{[k!]^2} = \sum\limits_{k=0}^\infty \frac{a_k t^k}{[k!]^2} \,, \end{aligned} \tag{12}\] associated with the umbral sequence \(a = \{ a_k \}\), where \[\begin{aligned} \label{13} a^k := a_k, \qquad \forall k \geq 0, \qquad a_0:= 1. \end{aligned} \tag{13}\]

According to Blissard [20], the solution \(b = \{ b_n \}\), of the umbral equation \[\begin{aligned} \label{14} e_1(at)\, e_1(bt) = 1, \end{aligned} \tag{14}\] is expressed by the Bell’s polynomials \(Y_n(f_1,g_1;f_2,g_2;\dots;f_n,g_n)\)   [20], as follows \[\begin{aligned} \label{15} \left\{ \begin{array}{l} b_0 := 1 , \, \\ \phantom{\rule{1 pt}{16 pt}} b_n = Y_n ( -1!, a_1; [2!]^2, a_2; -[3!]^2, a_3; \dots; (-1)^n [n!]^2, a_n), \quad (\forall \, n>0)\,. \end{array} \right. \end{aligned} \tag{15}\]

At present, the Blissard symbolic method is called the umbral calculus, a term coined by J.J. Sylvester.

The modern version of the umbral calculus is due to Rota and Roman [21, 22].

3. Basic definitions

In what follows, dealing with hypergeometric functions, and for typographical convenience, we use for the rising factorial the Pochhammer symbol according to the notation: \[\begin{aligned} \label{16} (x)_n = \left\{ \begin{array}{lr} x(x+1)\cdots (x+n-1) =\frac{\Gamma(x+n)}{\Gamma(x)} \,, \quad & n \geq 1 \,,\\ 1\,, \quad & n=0 \,. \end{array} \right. \end{aligned} \tag{16}\]

We put, for shortness: \[\begin{aligned} \label{17} {\cal T}_r(x):=\sum\limits_{\ell =0}^{r-1} \frac{x^\ell}{[\ell!]^2} , \end{aligned} \tag{17}\] that is the \(r\)th partial sum of the exponential series, which is polynomial of degree \(r-1\). Then, the \(r\)-associate Laguerre-type Stirling numbers of the second kind \({\cal S}(n,k;r)\) are defined by \[\begin{aligned} \label{18} \left(\sum\limits_{\ell =r}^\infty \frac{x^\ell}{[\ell!]^2}\right)^{k} = \left(e_1(x) – {\cal T}_r(x)\right)^k = [k!]^2 \sum\limits_{n=rk}^\infty {\cal S}(n,k;r) \, \frac{x^n}{[n!]^2} \ . \end{aligned} \tag{18}\]

Of course, when \(k=1\), we put by definition \({\cal S}(n,k;1)={\cal S}(n,k)\), and the \({\cal S}(n,k)\) are the Laguerre-type version (of order 1) of the above recalled \(S(n,k)\).

Note that \[\begin{aligned} \label{19} e_1(x) – {\cal T}_r(x) =\frac{x^r}{[r!]^2} \ _1F_2(1;r+1;r+1;x) =\frac{x^r}{[r!]^2} \,\sum\limits_{n=0}^\infty \frac{(1)_n}{[(r+1)_n]^2} \,\frac{x^n}{n!} \,. \end{aligned} \tag{19}\]

4. Hypergeometric Laguerre-type Bernoulli polynomials

In what follows we use the definition by Booth and Hassen [9] for the generalized Bernoulli polynomials, which is different from that introduced by Kurt [9] and in a previous article by Natalini and Bernardini [5].

The hypergeometric Laguerre-type Bernoulli polynomials are defined by the generating function: \[\begin{aligned} \label{20} \frac{\frac{x^r}{[r!]^2} \ e_1(t\, x)}{e_1(x) – {\cal T}_r(x)} = \frac{e_1(t\, x)}{_1F_2(1;r+1; r+1;x)} =\sum\limits_{n=0}^\infty {\cal B}^{[(r-1)^2,\, 1]}_n (t) \, \frac{x^n}{[n!]^2} \,. \end{aligned} \tag{20}\]

4.1. The exponential generating function

Since \[\begin{aligned} e_1(x) – {\cal T}_r(x) =\sum\limits_{n=r}^\infty\frac{x^n}{[n!]^2} = x^r \,\sum\limits_{n=0}^\infty \frac{1}{[(n+r)!]^2}\ x^n \,, \end{aligned}\] and \[\begin{aligned} \label{21} \frac{\frac{x^r}{[r!]^2} \ e_1(t\, x)}{e_1(x) – {\cal T}_r(x)} = \frac{e_1(t\, x)}{[r!]^2 \,\sum\limits_{n=0}^\infty \frac{1}{[(n+r)!]^2}\ x^n} \,, \end{aligned} \tag{21}\] we have the exponential generating function of the hypergeometric Bernoulli polynomials \({\cal B}^{[(r-1)^2,\, 1]}_n (t)\) \[\begin{aligned} \label{22} \frac{e_1(t\, x)}{[\Gamma(r+1)]^2 \,\sum\limits_{n=0}^\infty \frac{1}{[\Gamma(n+r+1)]^2} \, x^n } = \sum\limits_{n=0}^\infty {\cal B}^{[(r-1)^2,\, 1]}_n (t) \, \frac{x^n}{[n!]^2} \ . \end{aligned} \tag{22}\]

4.2. Recursive computation of the 2nd kind hypergeometric Laguerre-type Bernoulli numbers

In (19), putting \(t=0\), gives the exponential generating function [23] of the generalized hypergeometric Bernoulli numbers \({\cal B}^{[(r-1)^2,\, 1]}_n:= {\cal B}^{[(r-1)^2,\, 1]}_n(0)\): \[\begin{aligned} \label{23} \sum\limits_{n=0}^\infty {\cal B}^{[(r-1)^2,\, 1]}_n \, \frac{x^n}{[n!]^2} = \frac{1}{[\Gamma(r+1)]^2 \,\sum\limits_{n=0}^\infty \frac{1}{[\Gamma(n+r+1)]^2} \, x^n }\,, \end{aligned} \tag{23}\] which is valid even for non integer (in particular for fractional) values of the parameter \(r\).

From Eq. (20) we find: \[\begin{aligned} \sum\limits_{n=0}^\infty \sum\limits_{h=0}^n {n \choose h} {\cal B}^{[(r-1)^2,\, 1]}_{n-h} \ \frac{h!\, [\Gamma(r+1)]^2}{[\Gamma(h+r+1)]^2} \,\frac{x^n}{[n!]^2} = 1, \end{aligned}\] and therefore \[\begin{aligned} {\cal B}^{[(r-1)^2,\, 1]}_{0} = 1 \,, \end{aligned}\] and for \(n=1,2,3, \dots\) we find the \({\cal B}^{[(r-1)^2,\, 1]}_n\) numbers solving by recursion the triangular system: \[\begin{aligned} \sum\limits_{h=0}^n {n \choose h} {\cal B}^{[(r-1)^2,\, 1]}_{n-h} \ \frac{h!\, [\Gamma(r+1)]^2}{[\Gamma(h+r+1)]^2} = 0 \,. \end{aligned}\]

5. Representation in terms of Bell polynomials

We use a Laguerre-type version of the Blissard problem [20] showing that the reciprocal of a Laguerre-type Taylor series can be expressed in terms of Bell polynomials.
In fact, consider the sequences \(a := \{ a_k \} = (1, a_1, a_2, a_3, \dots )\), and \(b := \{ b_k \} = (b_0, b_1, b_2, b_3, \dots )\), and the function: \[\begin{aligned} \label{24} \begin{array}{c} \frac{1}{\sum\limits_{h=0}^\infty a_h \frac{x^h}{[h!]^2}} \, \quad (t\geq 0)\,. \end{array} \end{aligned} \tag{24}\]

Using the umbral formalism (that is, letting \(a_k \equiv a^k\) and \(b_k \equiv {\cal B}^k\)), the solution of the equation, \[\begin{aligned} \label{25} \begin{array}{c} \frac{1}{\sum\limits_{n=0}^\infty \frac{a^n t^n}{[n!]^2}} = \sum\limits_{n=0}^\infty \frac{{\cal B}^n t^n}{[n!]^2}\,, \qquad { i.e.} \qquad e_1(a\, t)\, e_1(b\, t) = 1\,, \end{array} \end{aligned} \tag{25}\] is given by \[\begin{aligned} \label{26} \left\{ \begin{array}{l} b_0 := 1 , \, \\ \phantom{\rule{1 pt}{26 pt}} b_n = Y_n ( -1!, a_1; [2!]^2, a_2; -[3!]^2, a_3; \dots; (-1)^n [n!]^2, a_n), \quad (\forall \, n>0), \end{array} \right. \end{aligned} \tag{26}\] where \(Y_n\) is the \(n\)th Bell polynomial [20].

The Bell polynomials are usually written in the form \[\begin{aligned} \quad Y_n(f_1,g_1;f_2,g_2;\dots;f_n,g_n)=\sum\limits_{k=1}^{n}B_{n,k}(g_1,g_2,\dots ,g_{n-k+1})\, f_{k}\,. \end{aligned}\] where the \(B_{n,h}\) are called Bell polynomials of the second kind [24].

By using this equation, the function (21) can be rewritten as \[\begin{aligned} \label{27} \begin{array}{c} \frac{1}{\sum\limits_{n=0}^\infty \frac{a_n t^n}{[n!]^2}} = 1+\sum\limits_{n=1}^\infty \sum\limits_{h=1}^n (-1)^h [h!]^2 \, B_{n,h}(a_1, a_2, \dots, a_{n-h+1}) \,\frac{t^n}{[n!]^2} \,. \end{array} \end{aligned} \tag{27}\]

It is convenient to introduce the definition \[\begin{aligned} \label{28} \begin{array}{c} C_n(a):=\sum\limits_{h=1}^n (-1)^h [h!]^2 \, B_{n,h}(a_1, a_2, \dots, a_{n-h+1}) \,, \quad C_0(a):= 1 \,, \end{array} \end{aligned} \tag{28}\] so that Eq. (24) becomes \[\begin{aligned} \label{29} \begin{array}{c} \frac{1}{\sum\limits_{n=0}^\infty \frac{a_n t^n}{[n!]^2}} =\sum\limits_{n=0}^\infty C_n(a) \,\frac{t^n}{[n!]^2} \,. \end{array} \end{aligned} \tag{29}\]

Applying this result to Eq. (20), with \[\begin{aligned} \begin{array}{c} a_n =\frac{[n!]^2 \, [\Gamma(r+1)]^2}{[\Gamma(n+r+1)]^2} \ , \end{array} \end{aligned}\] we find our result

Theorem 1. The 2nd kind hypergeometric Bernoulli numbers \({\cal B}^{[(r-1)^2,\, 1]}_n\) are expressed in terms of Bell polynomials by the equation \[\begin{aligned} \begin{array}{c} {\cal B}^{[(r-1)^2,\, 1]}_n = C_n \left(\frac{[n!]^2 \, [\Gamma(r+1)]^2}{[\Gamma(n+r+1)]^2} \, \right) \,, \end{array} \end{aligned}\] that is \[\begin{aligned} \label{30} {\cal B}^{[(r-1)^2,\, 1]}_n =&\sum\limits_{h=1}^n (-1)^h [h!]^2 \, B_{n,h}\left(\frac{[\Gamma(r+1)]^2}{[\Gamma(r+2)]^2}, \frac{2! \, [\Gamma(r+1)]^2}{[\Gamma(r+3)]^2}, \dots, \frac{(n-h+1)! \, [\Gamma(r+1)]^2}{[\Gamma(n-h+r+2)]^2}\right) \phantom{\rule{1 pt}{24 pt}} \notag\\ =&\sum\limits_{h=1}^n (-1)^h [h!]^2 \, B_{n,h}\left(\frac{1}{[(r+1)]^2}, \frac{2!}{[(r+1)_2]^2}, \dots, \frac{(n-h+1)!}{[(r+1)_{n-h+1}]^2}\right) . \end{aligned} \tag{30}\]

Many representation formulas in terms of Bell’s polynomials are reported in [25].

Remark 1. Since Laguerre-type exponentials exist for any integer \(n\), the above formulas can be generalized in a straightful way. It is sufficient to substitute the exponent 2 in the preceding equations (17)–(30).

For example, Eq. (23) becomes \[\begin{aligned} \sum\limits_{n=0}^\infty {\cal B}^{[(r-1)^{n+1},\, 1]}_n \, \frac{x^n}{[n!]^{n+1}} = \frac{1}{[\Gamma(r+1)]^{n+1} \,\sum\limits_{n=0}^\infty \frac{1}{[\Gamma(n+r+1)]^{n+1}} \, x^n }\,, \end{aligned}\] and Eq. (30) writes \[\begin{aligned} {\cal B}^{[(r-1)^{n+1},\, 1]}_n =\sum\limits_{h=1}^n (-1)^h [h!]^{n+1} \, B_{n,h}\left(\frac{1}{[(r+1)]^{n+1}}, \frac{2!}{[(r+1)_2]^{n+1}}, \dots, \frac{(n-h+1)!}{[(r+1)_{n-h+1}]^{n+1}}\right) \ . \end{aligned}\]

The proofs are obtained in the same way, so it is not necessary to further insist on this subject.

6. Conclusion

Starting from some recent generalizations of Bernoulli polynomials we have introduced in [17] a generalized hypergeometric versions of these polynomials and of the corresponding \(r\)-associated Stirling numbers.

In the same article we have introduced values of these mathematical entities for rational values of \(r\), thus expanding the original definitions.

In this article we have computed the Laguerre-type hypergeometric Bernoulli numbers of the second kind \({\cal B}^{[(r-1)^2,\, 1]}_n\) (of order 1), by using the same methodology of the preceding article. As this technique can be generalized to any integer order \(n\), we have provided a method to introduce infinite many sequences of numbers, which extend the classical hypergeometric Bernoulli polynomials and their \(r\)-associated Stirling numbers.

Appendix

Examples of Laguerre-type hypergeometric Bernoulli numbers of the second kind (of order 1) for fractional values of \(r\) are reported in the following table.

\(B_{n}^{[(r-1)^2,1]}\)

\(r = \dfrac{1}{2}\) \(r = \dfrac{3}{2}\) \(r = \dfrac{5}{2}\)
\(n=0\) \(1\) \(1\) \(1\)
\(n=1\) \(-\dfrac{4}{9}\) \(-\dfrac{4}{25}\) \(-\dfrac{4}{49}\)
\(n=2\) \(\dfrac{512}{2025}\) \(\dfrac{768}{30625}\) \(\dfrac{1024}{194481}\)
\(n=3\) \(-\dfrac{54272}{297675}\) \(-\dfrac{69632}{20671875}\) \(-\dfrac{44032}{384359283}\)
\(n=4\) \(\dfrac{10944512}{66976875}\) \(\dfrac{583221248}{3064088671875}\) \(-\dfrac{8227422208}{257813217024363}\)
\(n=5\) \(\dfrac{2592800768}{14587563375}\) \(\dfrac{10162405376}{95894626953125}\) \(\dfrac{273529372672}{63164238170968935}\)
\(n=6\) \(\dfrac{2076128595410944}{9059970923128125}\) \(-\dfrac{4831515683323904}{85636299234814453125}\) \(\dfrac{960183899411972096}{2922229493606610586243845}\)
\(n=7\) \(-\dfrac{3994590858182656}{116485344044021875}\) \(\dfrac{1265768498802458624}{88388894567362060546875}\) \(-\dfrac{1420764981357694681088}{7384473930343904951438196315}\)
\(n=8\) \(\dfrac{443310617419493408768}{757445926122522421875}\) \(\dfrac{72766396544266019012608}{92737671507965767730712890625}\) \(\dfrac{6249695099236333124373088}{49532117117991408029131635998445715}\)
\(n=9\) \(\dfrac{15043421255400186580041728}{13398460988773309120546875}\) \(-\dfrac{3409251665395799255146299392}{1064164780553907184709930419921875}\) \(\dfrac{16424267841858310022449512448}{172616369453000976699737962595675200035}\)
\(n=10\) \(\dfrac{174535890043789072451453845504}{72954620083870668161377734375}\) \(\dfrac{45170646204261335263651924364404864}{2330866190884346473961859920501708984375}\) \(\dfrac{23073711488077198082354872998011015584}{85639269294870119960158372035884358617364375}\)

References

  1. Ricci, P. E., & Natalini, P. (2022). Bell Polynomials and 2nd Kind Hypergeometric Bernoulli Numbers. In Seminar of I. Vekua Institute of Applied Mathematics (Vol. 48, pp. 36-45).

  2. Bretti, G., Natalini, P., & Ricci, P. E. (2004). Generalizations of the Bernoulli and Appell polynomials. In Abstract and Applied Analysis (Vol. 2004, No. 7, pp. 613-623). Hindawi Publishing Corporation.

  3. Carlitz, L. (1960). Note on Nörlund’s \({\cal B}^{(z)}_n\) Polynomial. Proceedings of the AMS, 11, 452–455.

  4. Elezović, N. (2016). Generalized Bernoulli polynomials and numbers, revisited. Mediterranean Journal of Mathematics, 13(1), 141-151.

  5. Natalini, P., & Bernardini, A. (2003). A generalization of the Bernoulli polynomials. Journal of Applied Mathematics, 2003(3), 155-163.

  6. Srivastava, H. M., & Todorov, P. G. (1988). An explicit formula for the generalized Bernoulli polynomials. Journal of Mathematical Analysis and Applications, 130(2), 509-513.

  7. Srivastava, H. M., Garg, M., & Choudhary, S. (2010). A new generalization of the Bernoulli and related polynomials. Russian Journal of Mathematical Physics, 17(2), 251-261.

  8. Ricci, P. E., Srivastava, R., & Natalini, P. (2021). A family of the r-associated Stirling numbers of the second kind and generalized Bernoulli polynomials. Axioms, 10(3), 219.

  9. Booth, R., & Hassen, A. (2011). Hypergeometric Bernoulli polynomials. Journal for Algebra and Number Theory Academia, 2(1), 1-7.

  10. Broder, A. Z. (1984). The \(r-\)Stirling numbers. Discrete Mathematics, 49(3), 241-259.

  11. Geleta, H. L., & Hassen, A. (2016). Fractional hypergeometric zeta functions. The Ramanujan Journal, 41(1), 421-436.

  12. Hassen, A., & Nguyen, H. D. (2008). Hypergeometric Bernoulli polynomials and Appell sequences. International Journal of Number Theory, 4(05), 767-774.

  13. Hu, S., & Kim, M. S. (2018). On hypergeometric Bernoulli numbers and polynomials. Acta Mathematica Hungarica, 154(1), 134-146.

  14. Kurt, B. (2010). A further generalization of the Bernoulli polynomials and on the 2D-Bernoulli polynomials \(B^2_n (x, y)\). Applied Mathematical Sciences, 4(47), 2315-2322.

  15. Kurt, B. (2013). Some relationships between the generalized Apostol-Bernoulli and Apostol-Euler polynomials. Turkish Journal of Analysis and Number Theory, 1(1), 54-58.

  16. Mihoubi, M., & Tiachachat, M. (2020). A new class of the t-stirling numbers and the generalized bernoulli polynomials. Integers: Electronic Journal of Combinatorial Number Theory, 20, 1-8.

  17. Ricci, P. E., & Natalini, P. (2022). Hypergeometric Bernoulli polynomials and r-associated stirling numbers of the second kind. Integers: Electronic Journal of Combinatorial Number Theory, 22, A56.

  18. Dattoli, G., & Ricci, P. E. (2003). Laguerre-type exponentials, and the relevant \(L\)-circular and \(L\)-hyperbolic functions. Georgian Mathematics Journal, 10, 481–494.

  19. Ricci, P. E. (2020). Laguerre-type exponentials, Laguerre derivatives and applications. A survey. Mathematics, 8(11), 2054.

  20. Riordan, J. (1958). An Introduction to Combinatorial Analysis. J. Wiley & Sons, Chichester.

  21. Roman, S. M., & Rota, G. C. (1978). The umbral calculus. Advances in Mathematics, 27(2), 95-188.

  22. Roman, S. M. (1984). The Umbral Calculus. Academic Press, New York.

  23. Manocha, H. L., & Srivastava, H. M. (1984). A Treatise on Generating Functions. Ellis Horwood, Chichester.

  24. Comtet, L. (2012). Advanced Combinatorics: The Art of Finite and Infinite Expansions. Springer Science & Business Media.

  25. Qi, F., Niu, D. W., Lim, D., & Yao, Y. H. (2020). Special values of the Bell polynomials of the second kind for some sequences and functions. Journal of Mathematical Analysis and Applications, 491(2), 124382.