Clausen function

From testwiki
Jump to navigation Jump to search

Template:Short description Template:Redir Template:Redir

File:Mplwp Clausen.svg
Graph of the Clausen function Template:Math

In mathematics, the Clausen function, introduced by Template:Harvs, is a transcendental, special function of a single variable. It can variously be expressed in the form of a definite integral, a trigonometric series, and various other forms. It is intimately connected with the polylogarithm, inverse tangent integral, polygamma function, Riemann zeta function, Dirichlet eta function, and Dirichlet beta function.

The Clausen function of order 2 – often referred to as the Clausen function, despite being but one of a class of many – is given by the integral:

Cl2(φ)=0φlog|2sinx2|dx:

In the range 0<φ<2π the sine function inside the absolute value sign remains strictly positive, so the absolute value signs may be omitted. The Clausen function also has the Fourier series representation:

Cl2(φ)=k=1sinkφk2=sinφ+sin2φ22+sin3φ32+sin4φ42+

The Clausen functions, as a class of functions, feature extensively in many areas of modern mathematical research, particularly in relation to the evaluation of many classes of logarithmic and polylogarithmic integrals, both definite and indefinite. They also have numerous applications with regard to the summation of hypergeometric series, summations involving the inverse of the central binomial coefficient, sums of the polygamma function, and Dirichlet L-series.

Basic properties

The Clausen function (of order 2) has simple zeros at all (integer) multiples of π, since if k is an integer, then sinkπ=0

Cl2(mπ)=0,m=0,±1,±2,±3,

It has maxima at θ=π3+2mπ[m]

Cl2(π3+2mπ)=1.01494160

and minima at θ=π3+2mπ[m]

Cl2(π3+2mπ)=1.01494160

The following properties are immediate consequences of the series definition:

Cl2(θ+2mπ)=Cl2(θ)
Cl2(θ)=Cl2(θ)

See Template:Harvtxt.

General definition

Template:Multiple image

More generally, one defines the two generalized Clausen functions:

Sz(θ)=k=1sinkθkz
Cz(θ)=k=1coskθkz

which are valid for complex z with Re z >1. The definition may be extended to all of the complex plane through analytic continuation.

When z is replaced with a non-negative integer, the standard Clausen functions are defined by the following Fourier series:

Cl2m+2(θ)=k=1sinkθk2m+2
Cl2m+1(θ)=k=1coskθk2m+1
Sl2m+2(θ)=k=1coskθk2m+2
Sl2m+1(θ)=k=1sinkθk2m+1

N.B. The SL-type Clausen functions have the alternative notation Glm(θ) and are sometimes referred to as the Glaisher–Clausen functions (after James Whitbread Lee Glaisher, hence the GL-notation).

Relation to the Bernoulli polynomials

The SL-type Clausen function are polynomials in θ, and are closely related to the Bernoulli polynomials. This connection is apparent from the Fourier series representations of the Bernoulli polynomials:

B2n1(x)=2(1)n(2n1)!(2π)2n1k=1sin2πkxk2n1.
B2n(x)=2(1)n1(2n)!(2π)2nk=1cos2πkxk2n.

Setting x=θ/2π in the above, and then rearranging the terms gives the following closed form (polynomial) expressions:

Sl2m(θ)=(1)m1(2π)2m2(2m)!B2m(θ2π),
Sl2m1(θ)=(1)m(2π)2m12(2m1)!B2m1(θ2π),

where the Bernoulli polynomials Bn(x) are defined in terms of the Bernoulli numbers BnBn(0) by the relation:

Bn(x)=j=0n(nj)Bjxnj.

Explicit evaluations derived from the above include:

Sl1(θ)=π2θ2,
Sl2(θ)=π26πθ2+θ24,
Sl3(θ)=π2θ6πθ24+θ312,
Sl4(θ)=π490π2θ212+πθ312θ448.

Duplication formula

For 0<θ<π, the duplication formula can be proven directly from the integral definition (see also Template:Harvtxt for the result – although no proof is given):

Cl2(2θ)=2Cl2(θ)2Cl2(πθ)

Denoting Catalan's constant by K=Cl2(π2), immediate consequences of the duplication formula include the relations:

Cl2(π4)Cl2(3π4)=K2
2Cl2(π3)=3Cl2(2π3)

For higher order Clausen functions, duplication formulae can be obtained from the one given above; simply replace θ with the dummy variable x, and integrate over the interval [0,θ]. Applying the same process repeatedly yields:

Cl3(2θ)=4Cl3(θ)+4Cl3(πθ)
Cl4(2θ)=8Cl4(θ)8Cl4(πθ)
Cl5(2θ)=16Cl5(θ)+16Cl5(πθ)
Cl6(2θ)=32Cl6(θ)32Cl6(πθ)

And more generally, upon induction on m,m1

Clm+1(2θ)=2m[Clm+1(θ)+(1)mClm+1(πθ)]

Use of the generalized duplication formula allows for an extension of the result for the Clausen function of order 2, involving Catalan's constant. For m1

Cl2m(π2)=22m1[Cl2m(π4)Cl2m(3π4)]=β(2m)

Where β(x) is the Dirichlet beta function.

Proof of the duplication formula

From the integral definition,

Cl2(2θ)=02θlog|2sinx2|dx

Apply the duplication formula for the sine function, sinx=2sinx2cosx2 to obtain

02θlog|(2sinx4)(2cosx4)|dx=02θlog|2sinx4|dx02θlog|2cosx4|dx

Apply the substitution x=2y,dx=2dy on both integrals:

20θlog|2sinx2|dx20θlog|2cosx2|dx=2Cl2(θ)20θlog|2cosx2|dx

On that last integral, set y=πx,x=πy,dx=dy, and use the trigonometric identity cos(xy)=cosxcosysinxsiny to show that:

cos(πy2)=siny2Cl2(2θ)=2Cl2(θ)20θlog|2cosx2|dx=2Cl2(θ)+2ππθlog|2siny2|dy=2Cl2(θ)2Cl2(πθ)+2Cl2(π)
Cl2(π)=0

Therefore,

Cl2(2θ)=2Cl2(θ)2Cl2(πθ).

Derivatives of general-order Clausen functions

Direct differentiation of the Fourier series expansions for the Clausen functions give:

ddθCl2m+2(θ)=ddθk=1sinkθk2m+2=k=1coskθk2m+1=Cl2m+1(θ)
ddθCl2m+1(θ)=ddθk=1coskθk2m+1=k=1sinkθk2m=Cl2m(θ)
ddθSl2m+2(θ)=ddθk=1coskθk2m+2=k=1sinkθk2m+1=Sl2m+1(θ)
ddθSl2m+1(θ)=ddθk=1sinkθk2m+1=k=1coskθk2m=Sl2m(θ)

By appealing to the First Fundamental Theorem Of Calculus, we also have:

ddθCl2(θ)=ddθ[0θlog|2sinx2|dx]=log|2sinθ2|=Cl1(θ)

Relation to the inverse tangent integral

The inverse tangent integral is defined on the interval 0<z<1 by

Ti2(z)=0ztan1xxdx=k=0(1)kz2k+1(2k+1)2

It has the following closed form in terms of the Clausen function:

Ti2(tanθ)=θlog(tanθ)+12Cl2(2θ)+12Cl2(π2θ)

Proof of the inverse tangent integral relation

From the integral definition of the inverse tangent integral, we have

Ti2(tanθ)=0tanθtan1xxdx

Performing an integration by parts

0tanθtan1xxdx=tan1xlogx|0tanθ0tanθlogx1+x2dx=
θlogtanθ0tanθlogx1+x2dx

Apply the substitution x=tany,y=tan1x,dy=dx1+x2 to obtain

θlogtanθ0θlog(tany)dy

For that last integral, apply the transform :y=x/2,dy=dx/2 to get

θlogtanθ1202θlog(tanx2)dx=θlogtanθ1202θlog(sin(x/2)cos(x/2))dx=θlogtanθ1202θlog(2sin(x/2)2cos(x/2))dx=θlogtanθ1202θlog(2sinx2)dx+1202θlog(2cosx2)dx=θlogtanθ+12Cl2(2θ)+1202θlog(2cosx2)dx.

Finally, as with the proof of the Duplication formula, the substitution x=(πy) reduces that last integral to

02θlog(2cosx2)dx=Cl2(π2θ)Cl2(π)=Cl2(π2θ)

Thus

Ti2(tanθ)=θlogtanθ+12Cl2(2θ)+12Cl2(π2θ).

Relation to the Barnes' G-function

For real 0<z<1, the Clausen function of second order can be expressed in terms of the Barnes G-function and (Euler) Gamma function:

Cl2(2πz)=2πlog(G(1z)G(1+z))+2πzlog(πsinπz)

Or equivalently

Cl2(2πz)=2πlog(G(1z)G(z))2πlogΓ(z)+2πzlog(πsinπz)

See Template:Harvtxt.

Relation to the polylogarithm

The Clausen functions represent the real and imaginary parts of the polylogarithm, on the unit circle:

Cl2m(θ)=(Li2m(eiθ)),m1
Cl2m+1(θ)=(Li2m+1(eiθ)),m0

This is easily seen by appealing to the series definition of the polylogarithm.

Lin(z)=k=1zkknLin(eiθ)=k=1(eiθ)kkn=k=1eikθkn

By Euler's theorem,

eiθ=cosθ+isinθ

and by de Moivre's Theorem (De Moivre's formula)

(cosθ+isinθ)k=coskθ+isinkθLin(eiθ)=k=1coskθkn+ik=1sinkθkn

Hence

Li2m(eiθ)=k=1coskθk2m+ik=1sinkθk2m=Sl2m(θ)+iCl2m(θ)
Li2m+1(eiθ)=k=1coskθk2m+1+ik=1sinkθk2m+1=Cl2m+1(θ)+iSl2m+1(θ)

Relation to the polygamma function

The Clausen functions are intimately connected to the polygamma function. Indeed, it is possible to express Clausen functions as linear combinations of sine functions and polygamma functions. One such relation is shown here, and proven below:

Cl2m(qπp)=1(2p)2m(2m1)!j=1psin(qjπp)[ψ2m1(j2p)+(1)qψ2m1(j+p2p)].

An immediate corollary is this equivalent formula in terms of the Hurwitz zeta function:

Cl2m(qπp)=1(2p)2mj=1psin(qjπp)[ζ(2m,j2p)+(1)qζ(2m,j+p2p)].

Template:Collapse top Let p and q be positive integers, such that q/p is a rational number 0<q/p<1, then, by the series definition for the higher order Clausen function (of even index):

Cl2m(qπp)=k=1sin(kqπ/p)k2m

We split this sum into exactly p-parts, so that the first series contains all, and only, those terms congruent to kp+1, the second series contains all terms congruent to kp+2, etc., up to the final p-th part, that contain all terms congruent to kp+p

Cl2m(qπp)=k=0sin[(kp+1)qπp](kp+1)2m+k=0sin[(kp+2)qπp](kp+2)2m+k=0sin[(kp+3)qπp](kp+3)2m++k=0sin[(kp+p2)qπp](kp+p2)2m+k=0sin[(kp+p1)qπp](kp+p1)2m+k=0sin[(kp+p)qπp](kp+p)2m

We can index these sums to form a double sum:

Cl2m(qπp)=j=1p{k=0sin[(kp+j)qπp](kp+j)2m}=j=1p1p2m{k=0sin[(kp+j)qπp](k+(j/p))2m}

Applying the addition formula for the sine function, sin(x+y)=sinxcosy+cosxsiny, the sine term in the numerator becomes:

sin[(kp+j)qπp]=sin(kqπ+qjπp)=sinkqπcosqjπp+coskqπsinqjπp
sinmπ0,cosmπ(1)mm=0,±1,±2,±3,
sin[(kp+j)qπp]=(1)kqsinqjπp

Consequently,

Cl2m(qπp)=j=1p1p2msin(qjπp){k=0(1)kq(k+(j/p))2m}

To convert the inner sum in the double sum into a non-alternating sum, split in two in parts in exactly the same way as the earlier sum was split into p-parts:

k=0(1)kq(k+(j/p))2m=k=0(1)(2k)q((2k)+(j/p))2m+k=0(1)(2k+1)q((2k+1)+(j/p))2m=k=01(2k+(j/p))2m+(1)qk=01(2k+1+(j/p))2m=12p[k=01(k+(j/2p))2m+(1)qk=01(k+(j+p2p))2m]

For m1, the polygamma function has the series representation

ψm(z)=(1)m+1m!k=01(k+z)m+1

So, in terms of the polygamma function, the previous inner sum becomes:

122m(2m1)![ψ2m1(j2p)+(1)qψ2m1(j+p2p)]

Plugging this back into the double sum gives the desired result:

Cl2m(qπp)=1(2p)2m(2m1)!j=1psin(qjπp)[ψ2m1(j2p)+(1)qψ2m1(j+p2p)]

Template:Collapse bottom

Relation to the generalized logsine integral

The generalized logsine integral is defined by:

snm(θ)=0θxmlognm1|2sinx2|dx

In this generalized notation, the Clausen function can be expressed in the form:

Cl2(θ)=s20(θ)

Kummer's relation

Ernst Kummer and Rogers give the relation

Li2(eiθ)=ζ(2)θ(2πθ)/4+iCl2(θ)

valid for 0θ2π.

Relation to the Lobachevsky function

The Lobachevsky function Λ or Л is essentially the same function with a change of variable:

Λ(θ)=0θlog|2sin(t)|dt=Cl2(2θ)/2

though the name "Lobachevsky function" is not quite historically accurate, as Lobachevsky's formulas for hyperbolic volume used the slightly different function

0θlog|sec(t)|dt=Λ(θ+π/2)+θlog2.

Relation to Dirichlet L-functions

For rational values of θ/π (that is, for θ/π=p/q for some integers p and q), the function sin(nθ) can be understood to represent a periodic orbit of an element in the cyclic group, and thus Cls(θ) can be expressed as a simple sum involving the Hurwitz zeta function.Template:Citation needed This allows relations between certain Dirichlet L-functions to be easily computed.

Series acceleration

A series acceleration for the Clausen function is given by

Cl2(θ)θ=1log|θ|+n=1ζ(2n)n(2n+1)(θ2π)2n

which holds for |θ|<2π. Here, ζ(s) is the Riemann zeta function. A more rapidly convergent form is given by

Cl2(θ)θ=3log[|θ|(1θ24π2)]2πθlog(2π+θ2πθ)+n=1ζ(2n)1n(2n+1)(θ2π)2n.

Convergence is aided by the fact that ζ(n)1 approaches zero rapidly for large values of n. Both forms are obtainable through the types of resummation techniques used to obtain rational zeta series Template:Harv.

Special values

Recall the Barnes G-function, the Catalan's constant K and the Gieseking constant V. Some special values include

Cl2(π2)=K
Cl2(π3)=V
Cl2(π3)=3πlog(G(23)G(13))3πlogΓ(13)+πlog(2π3)
Cl2(2π3)=2πlog(G(23)G(13))2πlogΓ(13)+2π3log(2π3)
Cl2(π4)=2πlog(G(78)G(18))2πlogΓ(18)+π4log(2π22)
Cl2(3π4)=2πlog(G(58)G(38))2πlogΓ(38)+3π4log(2π2+2)
Cl2(π6)=2πlog(G(1112)G(112))2πlogΓ(112)+π6log(2π231)
Cl2(5π6)=2πlog(G(712)G(512))2πlogΓ(512)+5π6log(2π23+1)

In general, from the Barnes G-function reflection formula,

Cl2(2πz)=2πlog(G(1z)G(z))2πlogΓ(z)+2πzlog(πsinπz)

Equivalently, using Euler's reflection formula for the gamma function, then,

Cl2(2πz)=2πlog(G(1z)G(z))2πlogΓ(z)+2πzlog(Γ(z)Γ(1z))

Generalized special values

Some special values for higher order Clausen functions include

Cl2m(0)=Cl2m(π)=Cl2m(2π)=0
Cl2m(π2)=β(2m)
Cl2m+1(0)=Cl2m+1(2π)=ζ(2m+1)
Cl2m+1(π)=η(2m+1)=(22m122m)ζ(2m+1)
Cl2m+1(π2)=122m+1η(2m+1)=(22m124m+1)ζ(2m+1)

where β(x) is the Dirichlet beta function, η(x) is the Dirichlet eta function (also called the alternating zeta function), and ζ(x) is the Riemann zeta function.

Integrals of the direct function

The following integrals are easily proven from the series representations of the Clausen function:

0θCl2m(x)dx=ζ(2m+1)Cl2m+1(θ)
0θCl2m+1(x)dx=Cl2m+2(θ)
0θSl2m(x)dx=Sl2m+1(θ)
0θSl2m+1(x)dx=ζ(2m+2)Cl2m+2(θ)

Fourier-analytic methods can be used to find the first moments of the square of the function Cl2(x) on the interval [0,π]:[1]

0πCl22(x)dx=ζ(4),
0πtCl22(x)dx=22190720π64ζ(5,1)2ζ(4,2),
0πt2Cl22(x)dx=23π[12ζ(5,1)+6ζ(4,2)2310080π6].

Here ζ denotes the multiple zeta function.

Integral evaluations involving the direct function

A large number of trigonometric and logarithmo-trigonometric integrals can be evaluated in terms of the Clausen function, and various common mathematical constants like K (Catalan's constant), log2, and the special cases of the zeta function, ζ(2) and ζ(3).

The examples listed below follow directly from the integral representation of the Clausen function, and the proofs require little more than basic trigonometry, integration by parts, and occasional term-by-term integration of the Fourier series definitions of the Clausen functions.

0θlog(sinx)dx=12Cl2(2θ)θlog2
0θlog(cosx)dx=12Cl2(π2θ)θlog2
0θlog(tanx)dx=12Cl2(2θ)12Cl2(π2θ)
0θlog(1+cosx)dx=2Cl2(πθ)θlog2
0θlog(1cosx)dx=2Cl2(θ)θlog2
0θlog(1+sinx)dx=2K2Cl2(π2+θ)θlog2
0θlog(1sinx)dx=2K+2Cl2(π2θ)θlog2

References

Template:Reflist