Khinchin's Constant

DOWNLOAD Mathematica Notebook

Let

 x=[a_0;a_1,...]=a_0+1/(a_1+1/(a_2+1/(a_3+...)))
(1)

be the simple continued fraction of a "generic" real number x, where the numbers a_i are the partial quotients. Khinchin (1934) considered the limit of the geometric mean

 G_n(x)=(a_1a_2...a_n)^(1/n)
(2)

as n->infty. Amazingly, except for a set of measure 0, this limit is a constant independent of x given by

 K=2.685452001...
(3)

(OEIS A002210), as proved in Kac (1959).

The constant is known as Khinchin's constant, and is commonly also spelled "Khintchine's constant" (Shanks and Wrench 1959, Bailey et al. 1997).

It is implemented as Khinchin, where its value is cached to 1100-digit precision. However, the numerical value of K is notoriously difficult to calculate to high precision, so computation of more digits get increasingly slower.

It is not known if K is irrational, let alone transcendental.

While it is known that almost all numbers x have limits G_n(x) that approach K, this fact has not been proven for any explicit real number x, e.g., a real number cast in terms of fundamental constants (Bailey et al. 1997).

KhinchinsConstant

The values (a_1,a_2,...,a_n)^(1/n) are plotted above for n=1 to 500 and x=pi, sin1, the Euler-Mascheroni constant gamma, and the Copeland-Erdős constant C. Interestingly, the shape of the curves is almost identical to the corresponding curves for the Lévy constant.

If p_n/q_n is the nth convergent of the continued fraction of x, then

lim_(n->infty)q_n^(1/n)=lim_(n->infty)((p_n)/x)^(1/n)
(4)
=e^(pi^2/(12ln2))
(5)
=3.27582...
(6)

(OEIS A086702) for almost all real x (Lévy 1936, Finch 2003). This number is sometimes called the Lévy constant.

Product and sum expressions for K include

 K=product_(n=1)^infty[1+1/(n(n+2))]^(lnn/ln2)
(7)

(Shanks and Wrench 1959; Khinchin 1997, p. 93; Borwein and Bailey 2003, p. 25; Havil 2003, p. 161) and

 K=exp[1/(ln2)sum_(m=1)^infty(H_(2m-1)^'[zeta(2m)-1])/m],
(8)

where zeta(z) is the Riemann zeta function and H_n^' is a an alternating harmonic number (Bailey et al. 1997). Gosper (pers. comm., Jun. 25, 1996) gave

 K=exp[1/(ln2)sum_(k=2)^infty((-1)^k(2-2^k)zeta^'(k))/k],
(9)

where zeta^'(z) is the derivative of the Riemann zeta function. An extremely rapidly converging sum originally due to Gosper (pers. comm., Jun. 25, 1996) and streamlined by O. Pavlyk (pers. comm., Apr. 24, 2006) is given by

 K=exp{-zeta^'(2,2)+1/(ln2)[sum_(k=2)^infty2(-1)^kf(k)]},
(10)

where

 f(k)=(lnk)/((k+2)k^(k+2))[2^(k+1)_2F_1(1,k+2;k+3;-2/k)-_2F_1(1,k+2;k+3;-1/k)]+((2^k-1)zeta^'(k+1,k))/(k+1),
(11)

zeta(s,a) is a Hurwitz zeta function, zeta^'(s,a)=partialzeta/partials, and _2F_1(a,b;c;z) is a hypergeometric function.

Khinchin's constant is also given by the integrals

K=2exp{1/(ln2)int_0^11/(x(1+x))ln[(pix(1-x^2))/(sin(pix))]dx}
(12)
=2exp[1/(ln2)int_0^11/(x(1+x))ln[Gamma(2-x)Gamma(2+x)]dx]
(13)

(Shanks and Wrench 1959) and

 K=exp[(pi^2)/(12ln2)+1/2ln2+1/(ln2)int_0^pi(ln(theta|cottheta|)dtheta)/theta].
(14)

Corless (1992) showed that

 lnK=int_0^1(ln|_x^(-1)_|)/((x+1)ln2)dx,
(15)

with an analogous formula for the Lévy constant.

KhinchinsConstant2

Real numbers x for which lim_(n->infty)G_n(x)!=K include x=e, sqrt(2), sqrt(3), and the golden ratio phi, plotted above.

Amazingly, the constant K is simply the limiting case K=K_0 of a class of means defined by

 K_p=lim_(n->infty)((a_1^p+a_2^p+...+a_n^p)/n)^(1/p)
(16)

for real p<1 whose values are given by

 K_p={sum_(k=1)^infty-k^plg[1-1/((k+1)^2)]}^(1/p)
(17)

(Ryll-Nardzewski 1951; Bailey et al. 1997; Khinchin 1997). An integral representation for K_p is given by

K_p=[1/(ln2)int_0^1(|_1/t_|^p)/(t+1)dt]^(1/p)
(18)
=[1/(ln2)sum_(k=1)^(infty)k^pln(1+1/(k(k+2)))]^(1/p)
(19)

for p=-1, -2, ... (Iosifescu and Kraaikamp 2002, p. 231).

The constant

 K_(-1)=lim_(n->infty)n/(a_1^(-1)+a_2^(-1)+...+a_n^(-1))
(20)

is sometimes known as the Khinchin harmonic mean. The following table summarizes the values for the first few nonpositive integers (Bailey et al. 1997, Plouffe).

pSloanevalue
0A0022102.685452001065306445309714835481795693820382293994462
-1A0874911.745405662407346863494596309683661067294936618777984
-2A0874921.450340328495630406052983076680697881408299979605904
-3A0874931.313507078687985766717339447072786828158129861484792
-4A0874941.236961809423730052626227244453422567420241131548937
-5A0874951.189003926465513154062363732771403397386092512639671
-6A0874961.156552374421514423152605998743410046840213070718761
-7A0874971.133323363950865794910289694908868363599098282411797
-8A0874981.115964408978716690619156419345349695769491182230400
-9A0874991.102543136670728013836093402522568351022221284149318
-10A0875001.091877041209612678276110979477638256493272651429656

Define the following quantity in terms of the kth partial quotient q_k,

 M(s,n,x)=(1/nsum_(k=1)^nq_k^s)^(1/s).
(21)

Then

 lim_(n->infty)M(1,n,x)=infty
(22)

for almost all real x (Khintchine 1934, 1936, Knuth 1981, Finch 2003), and

 M(1,n,x)∼O(lnn).
(23)

Furthermore, for s<1, the limiting value

 lim_(n->infty)M(s,n,x)=K(s)
(24)

exists and is a constant K(s) with probability 1 (Rockett and Szüsz 1992, Khinchin 1997).

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.