Riemann zeta fractional derivative—functional equation and link with primes

This paper outlines further properties concerning the fractional derivative of the Riemann ζ function. The functional equation, computed by the introduction of the Grünwald–Letnikov fractional derivative, is rewritten in a simplified form that reduces the computational cost. Additionally, a quasisymmetric form of the aforementioned functional equation is derived (symmetric up to one complex multiplicative constant). The second part of the paper examines the link with the distribution of prime numbers. The Dirichlet η function suggests the introduction of a complex strip as a fractional counterpart of the critical strip. Analytic properties are shown, particularly that a Dirichlet series can be linked with this strip and expressed as a sum of the fractional derivatives of ζ . Finally, Theorem 4.3 links the fractional derivative of ζ with the distribution of prime numbers in the left half-plane.


Introduction
In the last years, fractional calculus has played a significant role in many scientific fields and has been utilized for several applications (e.g., electromagnetism, signal processing). It is worth noting that the fractional derivative is not uniquely defined. This lack of uniqueness is undoubtedly the weakest point of the fractional calculus. Nevertheless, remarkable progress has been made toward addressing this and other issues in the last decade with relevant results [1][2][3][4][5][6][7], thus making the problem of independent interest for pure and applied mathematicians.
Fractional calculus has recently been applied to the theory of meromorphic functions (see, e.g., [8][9][10][11]). In particular, the α-order fractional derivative of the Riemann ζ function given by log α n n s , s ∈ C, α ∈ R >0 \ N, (1.1) was computed for the first time in 2015 [12]. By introducing both the Hurwitz ζ function and the Dirichlet series, two coherent generalizations of (1.1) were obtained. These frac-tional derivatives reduce to the derivative in (1.1) for a = 1 and f (n) = 1, which are the cases in which the Hurwitz ζ function and the Dirichlet series are nothing but the Riemann ζ function, respectively [13].
The main problem of fractional derivative in (1.1) is represented by the functional equation. Fortunately, the equation sought can be derived by introducing the forward Grünwald-Letnikov fractional derivative. The aforementioned functional equation was derived since this fractional derivative satisfies the generalized Leibniz rule [14,15]. However, as often happens in mathematics, the solution of one problem gives rise to another more important problem. The right-hand side of the functional equation in (2.7) entails extremely high computational cost due to the presence of three infinite series. Unfortunately, each finite approximation of these series with different numerical simulations produces an overflow. To minimize the computational cost, we follow the approach proposed by Apostol and Spira [16,17] for the integer-order derivatives of ζ . As a result, the aforementioned right-hand side is reduced from three to only one infinite series, with a substantial decrease in terms of computational cost (see Sect. 3). Furthermore, the functional equation of ζ (α) is written as the sum of sines and cosines. In 2017, Guariglia and Silvestrov [14] introduced the problem of the symmetry of (2.7). Accordingly, a quasisymmetric form of the functional equation of ζ (α) is derived and discussed in Sect. 3.
The second part of the paper investigates other properties of ζ (α) . The fractional equivalent of the critical strip [18], that is, (α, 1 + α), is presented and broadly discussed. The infinite series in (2.9) 2 associated with this strip is written in terms of ζ (α) . The Riemann conjecture entails that the critical strip is linked with the nontrivial zeros of ζ , which should all lie on Re s = 1/2 (critical line). Consequently, the distribution of prime numbers along the line Re s = 1/2 + α is considered. Based on these considerations, Sect. 4 discusses the link with prime numbers.
The rest of this paper is organized as follows. The next section provides remarks on ζ (α) . Section 3 reports and describes in detail equivalent forms of the functional equation of ζ (α) . Section 4 outlines the link with prime numbers, together with the analytic properties of the strip (α, 1 + α). Section 5 summarizes the results of this work and draws conclusions.

Notations and preliminary results
For the purposes of this paper α and k are always elements of R and N 0 , respectively, where N 0 = N ∪ {0}. A generic complex variable is denoted by s and the imaginary unit by i. Moreover, for every x ∈ R and n ∈ N 0 , the symbol x n indicates the so-called falling factorial [19].

Fractional calculus and ζ (α)
In recent years, the fractional calculus of holomorphic functions has merited increasing consideration from the international mathematics community. It has several applications in mathematics, theoretical physics, and information engineering [18]. Ortigueira is currently investigating broad applications in information engineering [15] and laying the foundations of fractional linear-time invariant systems. Furthermore, he has provided a suitable generalization of the Caputo derivative to the complex plane [12,20], that is, where f is a complex function of the complex variable s, m -1 < α < m ∈ N and θ = [0, 2π). The aforementioned fractional derivative can be written in terms of distribution theory [21]. The Riemann ζ function is defined as follows: Consequently, ζ converges for all complex numbers s such that Re s > 1. Riemann showed that ζ possesses a unique analytical continuation to the entire complex plane, except for a simple pole at s = 1 with residue 1 [22]. The fractional derivative in (2.1) enables us to compute the α-order fractional derivative of the Riemann ζ function, that is, the derivative in (1.1). The convergence of the complex series associated is shown in [12]. In particular, it pointwise converges in the half-plane Re s > 1 + α, and thus the convergence of ζ (α) depends on the fractional order α. Furthermore, the k-order integer derivative of the Riemann ζ function is given [16,17] by so that the derivative in (1.1) is a straightforward fractional generalization of the derivative in (2.2). As mentioned in the Introduction, (2.1) was applied to two generalizations of the Riemann ζ function (Hurwitz ζ function and Dirichlet series). The results achieved are in accordance with theory of the Riemann ζ function [13].

Functional equation of ζ (α)
The Riemann ζ function, introduced in 1859 [22], is characterized by the following functional equation: The literature calls (2.3) and (2.4) asymmetric equations. In fact, the functional equation of ζ can also be written symmetrically as follows: Introducing a variant of the Riemann ζ function, the Riemann Xi function, defined by (2.5) becomes ξ (s) = ξ (1-s). The function ξ is an entire function of s since π -s/2 Γ (-s/2)ζ (s) has simple poles at s = 0 and s = 1.
Unfortunately, the fractional derivative in (2.1) is unsuitable for deriving a functional equation of ζ (α) . In fact, the integral in (2.1) entails computational problems. Hence, the fractional operator in (2.1) is substituted with the forward Grünwald-Letnikov fractional derivative.
Definition 2.1 (see [15]) Let α ∈ R, h ∈ R ≥0 , and f be a complex-valued function of the complex variable s. The α-order forward Grünwald-Letnikov fractional derivative of f is defined by Unlike the Caputo derivative in (2.1), Definition 2.1 satisfies the generalized Leibniz rule. Thus, the fractional derivative ζ (α) can be recomputed by (2.6). The result achieved is in accordance with (1.1). The proof provided in [14] holds for α ∈ R >0 \ N. Since Definition 2.1 makes sense for any α ∈ R, (1.1) becomes The same reasoning fails for the functional equation of ζ (α) due to the presence of the generalized Leibniz rule (see [15] for more details). Accordingly, for any α ∈ R >0 \ N, it is where A α h,j,n = α h+j+n h!j!n! . The complete proof of (2.7) can be found in [14]. As previously mentioned, this equation suffers from problems of computational complexity. In fact, each numerical approximation of (2.7) produces a buffer overflow. Section 3 discusses such a computational issue.

Critical strip and ζ (α)
Prime numbers have played a central role in mathematics since the age of Euclid. In recent years, the distribution of prime numbers has been widely applied in science as well as engineering. The Riemann ζ function and the distribution of prime numbers are linked through the following representation: where P is the set of prime numbers [22]. ζ (s) = 0 for any negative even integers, that is, for s = -2, -4-, 6, . . . (trivial zeros). As a consequence of (2.8), there are no zeros in the half-plane Re s > 1, so that all the nontrivial zeros must belong to the (critical) strip (0, 1). The critical strip has great importance in the theory of prime numbers. Investigating ζ (α) leads to the introduction of the complex strip (α, 1 + α). Note that (2.1) and (2.6) allow the computation of the derivative η (α) , where η is the Dirichlet η function. Both are given [18] by Moreover, it can be shown [18] that (2.9) ζ (α) and η (α) converge for Re s > 1 + α and Re s > α, respectively [18]. Consequently, for each fixed α ∈ R, the fractional derivative ζ (α) is associated with a unique complex strip (α, 1 + α). Being

Equivalent forms of the functional equation
In the literature, Spira and Apostol were the first to investigate the higher integer derivatives of the Riemann ζ function [16,17]. In particular, Spira determined both the zero-free regions associated with the k-order integer derivative ζ (k) and the functional equation, that is, Note that the terms a hkn and b hkn are constants [23]. Replacing s by 1s yields

Simplified form
The undertaken study aims to repurpose the outcomes contained in [16] to write the functional equation of ζ (α) in a simplified form. In particular, it allows us to reduce the computational cost of (2.7) to only one infinite series (see Theorem 3.3). Furthermore, (3.1) and (3.2) express the functional equation of ζ (k) in terms of sines and cosines, that is, of complex exponentials. The next result, reported below, is due to Apostol [16].
By introducing the function ψ defined by Therefore taking into account the 2πi-periodicity of the complex exponential.

The fractional generalization of the functional equation in (3.3) is based on the following result.
Lemma 3.2 Let α ∈ R and w ∈ C such that Re w < 0. For any s ∈ C, we get Relying on Lemma 3.2, we are able to give the fractional counterpart of (3.3).

Theorem 3.3
Let α ∈ R >0 \ N and w = -log 2πiπ/2. For any s ∈ C, it is The generalized Leibniz rule gives and taking into account Lemma 3.2 it follows ξ is an entire function of s, thus substituting (3.7) into (3.6), the proof follows.
Theorem 3.3 shows that (3.5) has less computational cost than (2.7). Moreover, (3.2) can be generalized as the next result points out.
Theorem 3.4 takes on particular relevance in harmonic analysis. In fact, (3.8) expresses the functional equation of ζ (α) as a sum of sines and cosines. Additionally, (3.5) and (3.8), which are two different forms of (2.7), represent a step forward in our investigation of ζ (α) .

Symmetric form
The problem concerning the symmetric form of (2.7) is dealt with here.
Theorem 3.6 Let α ∈ R >0 \ N. Moreover, we set T α as follows: For any s ∈ C, it is where the symbol . = denotes equality up to a multiplicative constant and B α h,j = α h h!(h-j)! .
Proof Differentiating both members of (2.5) α times, it is Direct calculation leads to Substituting (3.14) and (3.15) into the left-hand side of (3.13) yields The same holds for the right-hand side of (3.13). In fact, we get Thus, Lemma 3.5 gives taking into account the 2πi-periodicity of the complex exponential. Substituting (3.16) and (3.17) into (3.13), the proof follows.
Note that (3.12) is the symmetric form of (2.7) up to the multiplicative constant e -iπα , that is, the fractional equivalent of (2.5). Therefore, Theorem 3.6 provides a quasisymmetric functional equation of ζ (α) by opening up new scenarios in functional fractional analysis.

Distribution of prime numbers and complex strip (α, 1 + α)
This section is devoted to investigation of the strip (α, 1 + α). First and foremost, η (α) is written in terms of ζ (α) . The second part discusses the link between ζ (α) and the distribution of prime numbers.

Representation of η (α) in terms of ζ (α)
Theorem 4.1 Let α ∈ R >0 \ N. For any s ∈ C such that Re s > α, the functions η (α) and ζ (α) are linked by Furthermore, the series in (2.9) 2 can be written in terms of ζ (α) as follows: As mentioned in Sect. 2, the condition Re s > α in Theorem 4.1 assures convergence of η (α) . Additionally, (4.1) is the fractional generalization of (2.9) 1 . In fact, with some caution it is Proof The generalized Leibniz rule reduces to the Leibniz rule for α ∈ N. Accordingly, the proof follows from (4.2).
Corollary 4.2 enables us to write the series in (2.9) 2 in terms of ζ (α) . Note that this result is in accordance with theory. Fix m ∈ N. Therefore

The role of ζ (α) in the distribution of prime numbers
The Euler product in (2.8) cannot be used to investigate the link between ζ (α) and the distribution of prime numbers due to the nonmultiplicativity of log α . In addition, the introduction of the strip (α, 1 + α) raises more delicate problems, such as the fractional generalization of both the Riemann hypothesis and the critical line Re s = 1/2. According to the approach given in [18], the principal candidate for this role is the line Re s = 1/2 + α. Unfortunately, the zeros of ζ (α) represent a major challenge. In fact, even now the zeros of the integer derivative ζ (k) remain an open problem. In 1974, assuming the Riemann hypothesis and k > 1, Levinson and Montgomery [24] proved that ζ has no zeros in the half-plane 0 < Re s < 1/2 and ζ (k) has at most finitely many zeros in Re s < 1/2. Spira showed that most of these zeros lie in 0 ≤ Re s ≤ 1/2 + δ for δ > 0 [23]. Hence, they are located close to the critical line. In particular, many scholars agree that nontrivial zeros of ζ (k) are located randomly to the right of Re s = 1/2. However, the left half-plane Re s < 0 is not zero-free since ζ has one pair of nontrivial zeros near -0.355 ± 3.591i [23]. Thus, the current state of the art does not facilitate any investigation into zeros of ζ (α) , that is, the localization of prime numbers near the line Re s = 1/2 + α. Nevertheless, ζ (α) can be linked with the distribution of prime numbers, as shown below. The aforementioned proof is strongly dependent on the assumption that Re s < 0. In particular, the set of prime numbers and ζ (α) are linked by (4.5) in the left half-plane. As mentioned in Sect. 2, ζ (α) converges for Re s > 1 + α, thus making this theorem true for any α < -1. Accordingly, Theorem 4.3 holds in the strip 1 + α < Re s < 0 with α < -1 (Fig. 2).

Conclusions
This paper analyzed the fractional derivative of the Riemann ζ function. In particular, the functional equation and the link with the distribution of prime numbers have been discussed. Equivalent forms of the functional equation in (2.7) are given. More precisely, such a functional equation is written in a quasisymmetric form. The introduction of the strip (α, 1 + α) has allowed investigation of the link with prime numbers. Analytic properties are given and treated. These results are far from being conclusive. However, Theorem 4.3 represents a step forward towards a final representation of ζ (α) in terms of prime numbers.