Various Structures of the Roots and Explicit Properties of q -cosine Bernoulli Polynomials and q -sine Bernoulli Polynomials

: In this paper, we deﬁne cosine Bernoulli polynomials and sine Bernoulli polynomials related to the q -number. Furthermore, we intend to ﬁnd the properties of these polynomials and check the structure of the roots. Through numerical experimentation, we look for various assumptions about the polynomials above.


Introduction
In the last three decades, the area of q-calculus has acted as a bridge between applied mathematics and engineering sciences. Many mathematicians studied the various fields of mathematics such as q-differential equations, q-integrals and differentials, q-series, q-trigonometric functions, q-hypergeometric functions, q-gamma and q-beta functions, and many properties, which are based on the discovery of q-numbers by Jackson (see ). This q-calculus plays important roles in many different areas of mathematics such as number theory, combinatorics, special functions and analysis, differential equations, and numerous interesting properties of them (see [1,2]).
The following diagram shows the variations of the different types of degenerate Bernoulli polynomials, Bernoulli polynomials, and q-Bernoulli polynomials. Those polynomials in the first row and the second row of the diagram were studied by Carlitz [5] and Kim and Ryoo [13,14], respectively. The study of these has produced beneficial results in combinatorics and number theory. The motivation of this paper is to investigate some explicit identities for q-cosine Bernoulli polynomials and q-sine Bernoulli polynomials in the third row of the diagram. t e q (t) − 1 e q (xt) = ∑ ∞ n=0 B n,q (x) t n n! (q-Bernoulli polynomials) t e q (t) − 1 e q (xt)COS q (ty) = ∑ ∞ n=0 C B n,q (x, y) t n n! (q-cosine Bernoulli polynomials) t e q (t) − 1 e q (xt)SIN q (ty) = ∑ ∞ n=0 S B n,q (x, y) t n n! (q-sine Bernoulli polynomials) Properties of certain q-orthogonal polynomials are connected to the q-oscillator algebra. The Wall and q-Laguerre polynomials are shown to arise as matrix elements of q-exponentials of the generators in a representation of this algebra (see [21]). Throughout this paper, the symbols N, Z, R and C denote the set of natural numbers, the set of integers, the set of real numbers, and the set of complex numbers, respectively.
We would like to begin by introducing several definitions related to q-numbers used in this paper. Definition 1. The Gaussian binomial coefficients are defined by: where m and r are non-negative integers. For r = 0, the value is one since the numerator and the denominator are both empty products. Like the classical binomial coefficients, the Gaussian binomial coefficients are center-symmetric. There are analogues of the binomial formula, and this definition has a number of properties (see [3,4,[6][7][8][9]11,12,19,22]).

Definition 2.
The q-analogues of (x − a) n and (x + a) n are defined by: Definition 3. Let z be any complex numbers with |z| < 1. Two forms of q-exponential functions can be expressed as: From Definition 3, we note that (1) e q (x)e q (y) = e q (x + y) if yx = qxy, and D q f (0) = f (0).
We can prove that f is differentiable at zero, and it is clear that D q x n = [n] q x n−1 .

Definition 5.
We define the q-integral as: If this function, f (x) is differentiable on the point x, and the q-derivative in Definition 4 goes to the ordinary derivative in the classical analysis when q → 1. Definition 6. The q-trigonometric functions are: where In various mathematics applications, which include number theory, finite difference calculus, combinatorial analysis, p-adic analytic number theory, and other fields, the Bernoulli, Euler, and Genocchi polynomials are widely studied. Acknowledging their significance, many mathematicians are familiar with these numbers and polynomials, and they have been studied for a long time. The previous definitions and theorems are also applied to polynomials, and their properties are studied in various ways in combination with Bernoulli, Euler, and Genocchi polynomials, which are considered important (see [1,5,[13][14][15][16][17][18][20][21][22][23][24][25][26]28]). The definition of q-Bernoulli polynomials is as follows: Definition 7. q-Bernoulli numbers, B n,q , and polynomials, B n,q (z), can be expressed as(see [17]): Recently, in [13], we confirmed the properties of cosine Bernoulli polynomials and sine Bernoulli polynomials. The definitions and representative properties of cosine Bernoulli polynomials and sine Bernoulli polynomials are the following. n (x, y) are defined by means of the generating functions: e tx sin(ty).

(10)
Based on the above, many studies can confirm various polynomials and their properties (see [28]).
The main aim of this paper is to identify the property of q-cosine Bernoulli polynomials and q-sine Bernoulli polynomials. In Section 2, we introduce cosine Bernoulli polynomials and sine Bernoulli polynomials combined with the q-number and confirm various properties and identities of polynomials. Here, we use the properties and exponential functions associated with the q-number. In Section 3, we can affirm the structure of the approximation roots of q-cosine Bernoulli polynomials and q-sine Bernoulli polynomials. By changing the q-numbers, we can find interesting properties and speculations.

Some Properties of q-cosine Bernoulli Polynomials and q-sine Bernoulli Polynomials
In this section, we introduce q-cosine Bernoulli polynomials and q-sine Bernoulli polynomials. By using the q-exponential function and q-trigonometric function, we can find various properties and identities. Lemma 1. Let y ∈ R and i = √ −1 ∈ C. Then, we have: where SI N q (x) = sin q −1 (x), COS q (x) = cos q −1 (x).
Proof. To find a relation between the q-exponential function, E q , and the q-trigonometric function, in particular q-sine and q-cosine functions, we can make the following equation: and: Therefore, we obtain the required results, Lemma 1. Lemma 2. Let x, y ∈ R and i = √ −1 ∈ C. Then, we have: Proof. (i) Using the property of two q-exponential functions, we can find: which is the required result.
(ii) By substituting −ity in (i), we can find the result (ii) in the same manner, so we omit this proof.
Theorem 2. For |q| < 1, we obtain: Proof. (i) By taking (x ⊕ iy) q instead of z in q-Bernoulli polynomials and applying Lemmas 1 and 2, we obtain: In a similar way, we can find the following equation, From (17) and (18), we can find: and: which is the required result shown of Theorem 2.
From Lemmas 1 and 2 and Theorem 2, we need to find some properties of e q (tx)COS q (ty) and e q (tx)SIN q (ty) in order to investigate some identities of special polynomials (see Definition 9). Let: Then, we can find Lemma 3.

Lemma 3.
Let k be a nonnegative integer. Then, we derive: where [x] is the greatest integer not exceeding x.
Proof. (i) We also know that the cosine functions consist of even terms in power series. In the same manner on quantum-calculus, we note that [12]). Multiplying COS q (ty) in the q-exponential function, we have: From Equation (21), it holds: By comparing the coefficients of both sides, we complete the proof of Lemma in quantum calculus (see [12]), and we can derive: By applying (21), we can change Equation (25) as follows: From the equation above, we can find the required result of Lemma 3(ii).
Now, we will introduce the q-cosine Bernoulli polynomials and q-sine Bernoulli polynomials considered in the previous lemmas and theorem. Definition 9. Let x, y be real numbers. Then, q-cosine Bernoulli polynomials and q-sine Bernoulli polynomials are defined by: and: respectively.

Corollary 1.
From Theorem 2 and Definition 9, it holds: Theorem 3. Let x, y ∈ R. Then, we have: where B n,q is the q-Bernoulli numbers.
Proof. (i) Here, we will show a relation between q-cosine Bernoulli polynomials and q-Bernoulli numbers. From the generating function of q-cosine Bernoulli polynomials, we can find the following Equation (31).
From (31), we can clearly obtain the required result of Theorem 3(i).
(ii) We omit the proof of Theorem 3(ii) because we can derive the required result if we use a similar proof method as the proof in Theorem 3(i).
Theorem 4. Let n ∈ Z and x, y ∈ R. Then, we derive: Proof. (i) We can transform the generating function of the q-cosine Bernoulli polynomials as follows when e q (t) = 1.
The left-hand side of Equation (33) can be transformed to: and the right-hand side of Equation (33) is transformed as: From (34) and (35), we complete the proof of Theorem 4(i).
(ii) Since we can find the required result from the same method (i), we omit the proof of Theorem 4(ii).

Theorem 5.
Let a be a real number. Then, we investigate: Proof. (i) In the Introduction, we referred to the two kinds of q-exponential functions and noted that e q (x)E q (−x) = 1. Using this property of q-exponential functions in q-cosine Bernoulli polynomials, we have: Applying Cauchy's product in Equation (37), we obtain: where the required result (i) is completed instantly.
(ii) From the q-sine Bernoulli polynomials, we can consider the following equation in a similar way to (i).
Using the same method (i), we are able to find the required result (ii).
(40) Theorem 6. Let x, y, and r be any real numbers. Then, we investigate: (i) C B n,q ((x ⊕ r) q , y) + S B n,q ((x r) q , y) = n ∑ k=0 n k q q ( n−k 2 ) r n−k C B k,q (x, y) + (−1) n−k S B k,q (x, y) Proof. (i) By substituting (x ⊕ r) q into x in the generating function of q-cosine Bernoulli polynomials, we can see the following equation.
Using Lemma 2(i) in Equation (42), we have: In a similar method, we can find the following equation.
By adding (43) with (44), we can derive the result (i) of Theorem 6.
(ii) We also can find the following equations, = t e q (t) − 1 e q (tx)COS q (ty)E q (−tr).

Corollary 3. From Theorem 6, it holds:
(i) C B n,q ((x ⊕ r) q , y) + C B n,q ((x r) q , y) Corollary 4. Setting r = 1 in Theorem 6, it holds: (47) Theorem 7. For x, y ∈ R, we derive: Proof. (i) Considering the q-partial derivative for x in the q-cosine Bernoulli polynomials, we have: Here, we note that: Therefore, we obtain: [n] qC B n−1,q (x, y) t n [n] q ! .

(51)
For q-exponential function E q (tx), we note that: and the q-derivative of q-trigonometric functions is: To find the required result, we derive: (54) , we can find: Hence, we investigate: and complete the proof of Theorem 7(i).
(ii) We note that: By using Equation (57) and applying a similar proof method as (i), we can find the required result, so we omit the proof of Theorem 7(ii).
Based on the content above, we introduce the new type of q-Bernoulli polynomials, which are polynomials of q-cosine Bernoulli polynomials when x = 0.

Definition 10.
For y ∈ R, we define the new type of q-Bernoulli polynomials as: Theorem 8. Let x, y be real numbers. Then, we have: Proof. (i) From the q-cosine Bernoulli polynomials, we find: and we obtain the required result (i).
(ii) We omit the proof of Theorem 8(ii) because the proof is very similar to (i).

S B n,q (x, y)
In this section, we would like to confirm the specific polynomial of q-cosine Bernoulli polynomials and q-sine Bernoulli polynomials. Mathematica can be used to identify the structure and build-up of the roots of polynomials to think bout a number of assumptions. Example 1. The specific polynomials of C n,q defined in Section 2 are shown below: Let us take a look at the specific shape and structure of the roots for the q-cosine Bernoulli polynomials defined in Definition 9, which are related to C n,q . Example 2. By using the q-cosine Bernoulli polynomials from the theorems obtained in Section 2, we can discover the following: Let us check the structure of the roots of these q-cosine Bernoulli polynomials. First, assume that y = 5. Let us change q and n in these circumstances. Then, the following Figure 1 can be obtained. First, when q is fixed at 0.9, the figure on the left is when n = 10, the center is when n = 20, and the right is when n = 30. Then, the following Figure 1 can be obtained. Based on Figure 1, we can assume that the larger the value of n, the more elliptical the structure becomes, excluding the three real roots. Figure 2 is the structure of approximation roots when n = 50. From the Figure above, we can make the following assumption.

Conjecture 3.
The greater the value of n is in C B n,0.9 (x, 5), the greater the distribution of approximate roots leaving out the three real roots shows an elliptical structure leaving out the three real roots.
Now, we are going to leave y at five and change the value of q. The following Figure 3 shows the distribution of the roots at q = 0.5. The figure on the left is the distribution of approximation roots at n = 10, and the figure on the right is at n = 30. In Figure 3, we can see that as n increases, the roots excluding three real roots become closer to a circle. Here, we can make the following assumption through Mathematica and also by observing Figures 1-3. Conjecture 4. C B ( n, q)(x, 5) always contains only three real roots in the range of 0.5 ≤ q < 1.
Let us look at the structure of the roots that supports the assumption above. Figure 4 is when y is fixed at five. In Figure 4, the left figure shows the structure of approximations at q = 0.9, and the right figure shows when q = 0.5. Let us check Figure 5 by observing the last picture of C B n,0.1 (x, 5). As before, fix the value of y at five. The left figure in Figure 5 shows the location of the approximations at n = 10, while the middle one shows when n = 20, and finally, the right one shows when n = 30. Most importantly, in Figure 5, we can see that as the value of q approaches zero and as n increases, the structure changes from an ellipse to a circular structure. From now on, let us look at the polynomial of S n,p,q associated with the q-sine function and find the q-sine Bernoulli polynomials.

Example 5.
If you obtain a specific polynomial of S n,p,q , this is the following: Here, we will visualize the structure of the roots in C n,p,q and in S n,p,q . Figure 6 shows the build-up of roots for C n,p,q . After fixing y at five, the figure on the far left is shown at q = 0.9 and the figure on the far right is at q = 0.1.
The following Figure 7 shows a stacked structure of roots for S n,p,q . This is also a form obtained under the same conditions shown in Figure 6, and compared to Figure 6, the lower the value of n for q = 0.9, the more likely it will have a slightly different change in location from C n,p,q . While the distribution of approximation roots appears similar, we can figure out that the approximation roots themselves are different. From now on, let us find the q-sine Bernoulli polynomials.
The following Figure 8 shows the build-up of the roots of q-sine Bernoulli polynomials. Similar forms to the structures of q-cosine Bernoulli polynomials, which were obtained earlier, can be found, and the structure in Figure 8 shows the stacked form of approximation roots. When fixed at y = 5, the left figure in Figure 8 is the structure that can be seen when q = 0.9, the middle figure is the structure that can be found when q = 0.5, and the right is when q = 0.1. In Figure 8, blue means the value of n is small, and red appears when the value of n is 30. As a result, we can see that the stacking structure is changing as the q-number changes. From Figure 8, we can make the following assumption. Conjecture 7. S B n,q (x, 5) obtains a value for approximation roots with a diameter of four as n increases and q approaches zero.

Conflicts of Interest:
The authors declare that there are no conflicts of interest regarding the publication of this paper.