Hurwitz Zeta Function

DOWNLOAD Mathematica Notebook

The Hurwitz zeta function zeta(s,a) is a generalization of the Riemann zeta function zeta(s) that is also known as the generalized zeta function. It is classically defined by the formula

 zeta(s,a)=sum_(k=0)^infty1/((k+a)^s)
(1)

for R[s]>1 and by analytic continuation to other s!=1, where any term with k+a=0 is excluded. It is implemented in this form in the Wolfram Language as HurwitzZeta[s, a].

The slightly different form

 zeta^*(s,a)=sum_(k=0)^infty1/([(a+k)^2]^(s/2))
(2)

is implemented in the Wolfram Language as Zeta[s, a]. Note that the two are identical only for R[a]>0.

HurwitzZetaFunction

The plot above shows zeta(s,a) for real s and a, with the zero contour indicated in black.

For a>-1, a globally convergent series for zeta(s,a) (which, for fixed a, gives an analytic continuation of zeta(s,a) to the entire complex s-plane except the point s=1) is given by

 zeta(s,a)=1/(s-1)sum_(n=0)^infty1/(n+1)sum_(k=0)^n(-1)^k(n; k)(a+k)^(1-s)
(3)

(Hasse 1930).

The Hurwitz zeta function is implemented in the Wolfram Language as Zeta[s, a].

For a=1, zeta(s,a) reduces to the Riemann zeta function zeta(s),

 zeta(s,1)=zeta(s).
(4)

If the singular term is excluded from the sum definition of zeta(s,a), then zeta(s,0)=zeta(s) as well.

The Hurwitz zeta function is given by the integral

 zeta(s,a)=1/(Gamma(s))int_0^infty(t^(s-1)dt)/(e^(at)(1-e^(-t)))
(5)

for R[s]>1 and R[a]>0.

HurwitzZetaZeros

The plot above illustrates the complex zeros of zeta(s,a) (Trott 1999), where s=x+iy. Here, the complex s-plane is horizontal and the real a-line is vertical and runs from a=1/2 at the bottom to a=1 at the top. The upper line is the critical line R[s]=1/2, which contains zeros of zeta(s)=zeta(s,1). The lower two lines are R[s]=0 and R[s]=1/2 (again), which contain zeros of 2^s-1 and zeta(s), respectively, since zeta(s,1/2)=(2^s-1)zeta(s); cf. equation (9) below.

This plot also appeared on the cover of the March 2004 issue of FOCUS, the Mathematical Association of America's news magazine.

The Hurwitz zeta function can also be given by the functional equation

 zeta(s,p/q)=2Gamma(1-s)(2piq)^(s-1)sum_(n=1)^qsin((pis)/2+(2pinp)/q)zeta(1-s,n/q)
(6)

(Apostol 1995, Miller and Adamchik 1999), or the integral

 zeta(s,a)=1/2a^(-s)+(a^(1-s))/(s-1)+2int_0^infty(a^2+y^2)^(-s/2){sin[stan^(-1)(y/a)]}(dy)/(e^(2piy)-1).
(7)

If R[z]<0 and 0<a<=1, then

 zeta(z,a)=(2Gamma(1-z))/((2pi)^(1-z))[sin((piz)/2)sum_(n=1)^infty(cos(2pian))/(n^(1-z))+cos((piz)/2)sum_(n=1)^infty(sin(2pian))/(n^(1-z))]
(8)

(Hurwitz 1882; Whittaker and Watson 1990, pp. 268-269).

The Hurwitz zeta function satisfies

 zeta(-n,a)=-(B_(n+1)(a))/(n+1)
(9)

for n>=0 (Apostol 1995, p. 264), where B_k(a) is a Bernoulli polynomial, giving the special case

 zeta(0,a)=1/2-a.
(10)

In addition,

zeta(s,1/2)=sum_(k=0)^(infty)(k+1/2)^(-s)
(11)
=2^ssum_(k=0)^(infty)(2k+1)^(-s)
(12)
=2^s[zeta(s)-sum_(k=1)^(infty)(2k)^(-s)]
(13)
=2^s(1-2^(-s))zeta(s)
(14)
=(2^s-1)zeta(s).
(15)

Derivative identities include

d/(ds)zeta(0,a)=ln[Gamma(a)]-1/2ln(2pi)
(16)
d/(ds)zeta(0,0)=-1/2ln(2pi),
(17)

where Gamma(z) is the gamma function (Bailey et al. 2006, p. 179). The definition (1) implies that

 d/(da)zeta(s,a)=-szeta(s+1,a)
(18)

for s!=0,1.

In the limit,

 lim_(s->1)[zeta(s,a)-1/(s-1)]=-psi_0(a)
(19)

(Whittaker and Watson 1990, p. 271; Allouche 1992), where psi_0(z) is the digamma function.

The polygamma function psi_m(z) can be expressed in terms of the Hurwitz zeta function by

 psi_m(z)=(-1)^(m+1)m!zeta(1+m,z).
(20)

For positive integers k, p, and q>p,

 zeta^'(-2k+1,p/q)=([psi(2k)-ln(2piq)]B_(2k)(p/q))/(2k)-([psi(2k)-ln(2pi)]B_(2k))/(q^(2k)2k)+((-1)^(k+1)pi)/((2piq)^(2k))sum_(n=1)^(q-1)sin((2pipn)/q)psi_((2k-1))(n/q)+((-1)^(k+1)2(2k-1)!)/((2piq)^(2k))sum_(n=1)^(q-1)cos((2pipn)/q)zeta^'(2k,n/q)+(zeta^'(-2k+1))/(q^(2k)),
(21)

where B_n is a Bernoulli number, B_n(x) a Bernoulli polynomial, psi_n(z) is a polygamma function, and zeta(z) is the Riemann zeta function (Miller and Adamchik 1999). Miller and Adamchik (1999) also give the closed-form expressions (where a large number of typos have been corrected in the expressions below)

zeta^'(1-2k,1/2)=-(B_(2k)ln2)/(4^kk)-((2^(2k-1)-1)zeta^'(-2k+1))/(2^(2k-1))
(22)
zeta^'(1-2k,1/3; 2/3)=∓(sqrt(3)(9^k-1)B_(2k)pi)/(8k·9^k)-(3B_(2k)ln3)/(4k·9^k)∓((-1)^kpsi_(2k-1)(1/3))/(2sqrt(3)(6pi)^(2k-1))-((9^k-1)zeta^'(1-2k))/(2·9^k)
(23)
zeta^'(1-2k,1/4; 3/4)=∓((4^k-1)B_(2k)pi)/(4^(k+1)k)+((4^(k-1)-1)B_(2k)ln2)/(k·2^(4k-1))∓((-1)^kpsi_(2k-1)(1/4))/(4(8pi)^(2k-1))-((4^k-2)zeta^'(1-2k))/(2^(4k))
(24)
zeta^'(1-2k,1/6; 5/6)=∓((9^k-1)(2^(2k-1)+1)B_(2k)pi)/(8sqrt(3)k·6^(2k-1))+(B_(2k)(3^(2k-1)-1)ln2)/(4k·6^(2k-1))+(B_(2k)(2^(2k-1)-1)ln3)/(4k·6^(2k-1))∓((-1)^k(2^(2k-1)+1)psi_(2k-1)(1/3))/(2sqrt(3)(12pi)^(2k-1))+((2^(2k-1)-1)(3^(2k-1)-1)zeta^'(1-2k))/(2·6^(2k-1)),
(25)

where zeta^'(z_0,a) means dzeta(z,a)/dz|_(z=z_0), zeta^'(z_0) means dzeta(z)/dz|_(z=z_0), and the upper and lower fractions on the left side of the equations correspond to the plus and minus signs, respectively, on the right side.

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.