In mathematics, the Glaisher–Kinkelin constant or Glaisher's constant, typically denoted A, is a mathematical constant, related to special functions like the K-function and the Barnes G-function. The constant also appears in a number of sums and integrals, especially those involving the gamma function and the Riemann zeta function. It is named after mathematicians James Whitbread Lee Glaisher and Hermann Kinkelin. Its approximate value is:
A = 1.28242712910062263687... (sequence A074962 in the OEIS). Glaisher's constant plays a role both in mathematics and in physics. It appears when giving a closed form expression for Porter's constant, when estimating the efficiency of the Euclidean algorithm. It also is connected to solutions of Painlevé differential equations and the Gaudin model.
Definition The Glaisher–Kinkelin constant A can be defined via the following limit:
A = lim n → ∞ H ( n ) n n 2 2 + n 2 + 1 12 e − n 2 4 {\displaystyle A=\lim _{n\rightarrow \infty }{\frac {H(n)}{n^{{\tfrac {n^{2}}{2}}+{\tfrac {n}{2}}+{\tfrac {1}{12}}}\,e^{-{\tfrac {n^{2}}{4}}}}}}
where H ( n ) {\displaystyle H(n)} is the hyperfactorial: H ( n ) = ∏ i = 1 n i i = 1 1 ⋅ 2 2 ⋅ 3 3 ⋅ . . . ⋅ n n {\displaystyle H(n)=\prod _{i=1}^{n}i^{i}=1^{1}\cdot 2^{2}\cdot 3^{3}\cdot {...}\cdot n^{n}} An analogous limit, presenting a similarity between A {\displaystyle A} and 2 π {\displaystyle {\sqrt {2\pi }}} , is given by Stirling's formula as:
2 π = lim n → ∞ n ! n n + 1 2 e − n {\displaystyle {\sqrt {2\pi }}=\lim _{n\to \infty }{\frac {n!}{n^{n+{\frac {1}{2}}}\,e^{-n}}}}
with n ! = ∏ i = 1 n i = 1 ⋅ 2 ⋅ 3 ⋅ . . . ⋅ n {\displaystyle n!=\prod _{i=1}^{n}i=1\cdot 2\cdot 3\cdot {...}\cdot n} which shows that just as π is obtained from approximation of the factorials, A is obtained from the approximation of the hyperfactorials.
Relation to special functions Just as the factorials can be extended to the complex numbers by the gamma function such that Γ ( n ) = ( n − 1 ) ! {\displaystyle \Gamma (n)=(n-1)!} for positive integers n, the hyperfactorials can be extended by the K-function with K ( n ) = H ( n − 1 ) {\displaystyle K(n)=H(n-1)} also for positive integers n, where:
K ( z ) = ( 2 π ) − z − 1 2 exp [ ( z 2 ) + ∫ 0 z − 1 ln Γ ( t + 1 ) d t ] {\displaystyle K(z)=(2\pi )^{-{\frac {z-1}{2}}}\exp \left[{\binom {z}{2}}+\int _{0}^{z-1}\ln \Gamma (t+1)\,dt\right]}
This gives:
A = lim n → ∞ K ( n + 1 ) n n 2 2 + n 2 + 1 12 e − n 2 4 {\displaystyle A=\lim _{n\rightarrow \infty }{\frac {K(n+1)}{n^{{\tfrac {n^{2}}{2}}+{\tfrac {n}{2}}+{\tfrac {1}{12}}}\,e^{-{\tfrac {n^{2}}{4}}}}}} . A related function is the Barnes G-function which is given by
G ( n ) = ( Γ ( n ) ) n − 1 K ( n ) {\displaystyle G(n)={\frac {(\Gamma (n))^{n-1}}{K(n)}}}
and for which a similar limit exists:
1 A = lim n → ∞ G ( n + 1 ) ( 2 π ) n 2 n n 2 2 − 1 12 e − 3 n 2 4 + 1 12 {\displaystyle {\frac {1}{A}}=\lim _{n\rightarrow \infty }{\frac {G(n+1)}{\left(2\pi \right)^{\frac {n}{2}}n^{{\frac {n^{2}}{2}}-{\frac {1}{12}}}e^{-{\frac {3n^{2}}{4}}+{\frac {1}{12}}}}}} . The Glaisher-Kinkelin constant also appears in the evaluation of the K-function and Barnes-G function at half and quarter integer values such as:
K ( 1 / 2 ) = A 3 / 2 2 1 / 24 e 1 / 8 {\displaystyle K(1/2)={\frac {A^{3/2}}{2^{1/24}e^{1/8}}}}
K ( 1 / 4 ) = A 9 / 8 exp ( G 4 π − 3 32 ) {\displaystyle K(1/4)=A^{9/8}\exp \left({\frac {G}{4\pi }}-{\frac {3}{32}}\right)}
G ( 1 / 2 ) = 2 1 / 24 e 1 / 8 A 3 / 2 π 1 / 4 {\displaystyle G(1/2)={\frac {2^{1/24}e^{1/8}}{A^{3/2}\pi ^{1/4}}}}
G ( 1 / 4 ) = 1 2 9 / 16 A 9 / 8 π 3 / 16 ϖ 3 / 8 exp ( 3 32 − G 4 π ) {\displaystyle G(1/4)={\frac {1}{2^{9/16}A^{9/8}\pi ^{3/16}\varpi ^{3/8}}}\exp \left({\frac {3}{32}}-{\frac {G}{4\pi }}\right)}
with G {\displaystyle G} being Catalan's constant and ϖ = Γ ( 1 / 4 ) 2 2 2 π {\displaystyle \varpi ={\frac {\Gamma (1/4)^{2}}{2{\sqrt {2\pi }}}}} being the lemniscate constant. Similar to the gamma function, there exists a multiplication formula for the K-function. It involves Glaisher's constant:
∏ j = 1 n − 1 K ( j n ) = A n 2 − 1 n n − 1 12 n e 1 − n 2 12 n {\displaystyle \prod _{j=1}^{n-1}K\left({\frac {j}{n}}\right)=A^{\frac {n^{2}-1}{n}}n^{-{\frac {1}{12n}}}e^{\frac {1-n^{2}}{12n}}}
The logarithm of G(z + 1) has the following asymptotic expansion, as established by Barnes:
ln G ( z + 1 ) = z 2 2 ln z − 3 z 2 4 + z 2 ln 2 π − 1 12 ln z + ( 1 12 − ln A ) + ∑ k = 1 N B 2 k + 2 4 k ( k + 1 ) z 2 k + O ( 1 z 2 N + 2 ) {\displaystyle \ln G(z+1)={\frac {z^{2}}{2}}\ln z-{\frac {3z^{2}}{4}}+{\frac {z}{2}}\ln 2\pi -{\frac {1}{12}}\ln z+\left({\frac {1}{12}}-\ln A\right)+\sum _{k=1}^{N}{\frac {B_{2k+2}}{4k\left(k+1\right)z^{2k}}}+O\left({\frac {1}{z^{2N+2}}}\right)}
The Glaisher-Kinkelin constant is related to the derivatives of the Euler-constant function:
γ ′ ( − 1 ) = 11 6 ln 2 + 6 ln A − 3 2 ln π − 1 {\displaystyle \gamma '(-1)={\frac {11}{6}}\ln 2+6\ln A-{\frac {3}{2}}\ln \pi -1}
γ ″ ( − 1 ) = 10 3 ln 2 + 24 ln A − 4 ln π − 7 ζ ( 3 ) 2 π 2 − 13 4 {\displaystyle \gamma ''(-1)={\frac {10}{3}}\ln 2+24\ln A-4\ln \pi -{\frac {7\zeta (3)}{2\pi ^{2}}}-{\frac {13}{4}}}
A {\displaystyle A} also is related to the Lerch transcendent:
∂ Φ ∂ s ( − 1 , − 1 , 1 ) = 3 ln A − 1 3 ln 2 − 1 4 {\displaystyle {\frac {\partial \Phi }{\partial s}}(-1,-1,1)=3\ln A-{\frac {1}{3}}\ln 2-{\frac {1}{4}}}
Glaisher's constant may be used to give values of the derivative of the Riemann zeta function as closed form expressions, such as:
ζ ′ ( − 1 ) = 1 12 − ln A {\displaystyle \zeta '(-1)={\frac {1}{12}}-\ln A}
ζ ′ ( 2 ) = π 2 6 ( γ + ln 2 π − 12 ln A ) {\displaystyle \zeta '(2)={\frac {\pi ^{2}}{6}}\left(\gamma +\ln 2\pi -12\ln A\right)}
where γ is the Euler–Mascheroni constant.
Series expressions The above formula for ζ ′ ( 2 ) {\displaystyle \zeta '(2)} gives the following series:
∑ k = 2 ∞ ln k k 2 = π 2 6 ( 12 ln A − γ − ln 2 π ) {\displaystyle \sum _{k=2}^{\infty }{\frac {\ln k}{k^{2}}}={\frac {\pi ^{2}}{6}}\left(12\ln A-\gamma -\ln 2\pi \right)}
which directly leads to the following product found by Glaisher:
∏ k = 1 ∞ k 1 k 2 = ( A 12 2 π e γ ) π 2 6 {\displaystyle \prod _{k=1}^{\infty }k^{\frac {1}{k^{2}}}=\left({\frac {A^{12}}{2\pi e^{\gamma }}}\right)^{\frac {\pi ^{2}}{6}}}
Similarly it is
∑ k ≥ 3 k odd ln k k 2 = π 2 24 ( 36 ln A − 3 γ − ln 16 π 3 ) {\displaystyle \sum _{k\geq 3}^{k{\text{ odd}}}{\frac {\ln k}{k^{2}}}={\frac {\pi ^{2}}{24}}\left(36\ln A-3\gamma -\ln 16\pi ^{3}\right)}
which gives:
∏ k ≥ 3 k odd k 1 k 2 = ( A 36 16 π 3 e 3 γ ) π 2 24 {\displaystyle \prod _{k\geq 3}^{k{\text{ odd}}}k^{\frac {1}{k^{2}}}=\left({\frac {A^{36}}{16\pi ^{3}e^{3\gamma }}}\right)^{\frac {\pi ^{2}}{24}}}
An alternative product formula, defined over the prime numbers, reads:
∏ p prime p 1 p 2 − 1 = A 12 2 π e γ , {\displaystyle \prod _{p{\text{ prime}}}p^{\frac {1}{p^{2}-1}}={\frac {A^{12}}{2\pi e^{\gamma }}},}
Another product is given by:
∏ n = 1 ∞ ( e n n ( n + 1 ) n ) ( − 1 ) n − 1 = 2 1 / 6 e π A 6 {\displaystyle \prod _{n=1}^{\infty }\left({\frac {en^{n}}{(n+1)^{n}}}\right)^{(-1)^{n-1}}={\frac {2^{1/6}e{\sqrt {\pi }}}{A^{6}}}}
A series involving the cosine integral is:
∑ k = 1 ∞ Ci ( 2 k π ) k 2 = π 2 2 ( 4 ln A − 1 ) {\displaystyle \sum _{k=1}^{\infty }{\frac {{\text{Ci}}(2k\pi )}{k^{2}}}={\frac {\pi ^{2}}{2}}(4\ln A-1)}
Helmut Hasse gave another series representation for the logarithm of Glaisher's constant, following from a series for the Riemann zeta function:
ln A = 1 8 − 1 2 ∑ n = 0 ∞ 1 n + 1 ∑ k = 0 n ( − 1 ) k ( n k ) ( k + 1 ) 2 ln ( k + 1 ) {\displaystyle \ln A={\frac {1}{8}}-{\frac {1}{2}}\sum _{n=0}^{\infty }{\frac {1}{n+1}}\sum _{k=0}^{n}(-1)^{k}{\binom {n}{k}}(k+1)^{2}\ln(k+1)}
Integrals The following are some definite integrals involving Glaisher's constant:
∫ 0 ∞ x ln x e 2 π x − 1 d x = 1 24 − 1 2 ln A {\displaystyle \int _{0}^{\infty }{\frac {x\ln x}{e^{2\pi x}-1}}\,dx={\frac {1}{24}}-{\frac {1}{2}}\ln A}
∫ 0 1 2 ln Γ ( x ) d x = 3 2 ln A + 5 24 ln 2 + 1 4 ln π {\displaystyle \int _{0}^{\frac {1}{2}}\ln \Gamma (x)\,dx={\frac {3}{2}}\ln A+{\frac {5}{24}}\ln 2+{\frac {1}{4}}\ln \pi }
the latter being a special case of:
∫ 0 z ln Γ ( x ) d x = z ( 1 − z ) 2 + z 2 ln 2 π + z ln Γ ( z ) − ln G ( 1 + z ) {\displaystyle \int _{0}^{z}\ln \Gamma (x)\,dx={\frac {z(1-z)}{2}}+{\frac {z}{2}}\ln 2\pi +z\ln \Gamma (z)-\ln G(1+z)}
A double integral is given by:
∫ 0 1 ∫ 0 1 − x ( 1 + x y ) 2 ln x y d x d y = 6 ln A − 1 6 ln 2 − 1 2 ln π − 1 2 {\displaystyle \int _{0}^{1}\int _{0}^{1}{\frac {-x}{(1+xy)^{2}\ln xy}}dxdy=6\ln A-{\frac {1}{6}}\ln 2-{\frac {1}{2}}\ln \pi -{\frac {1}{2}}}
Generalizations The Glaisher-Kinkelin constant can be viewed as the first constant in a sequence of infinitely many so-called generalized Glaisher constants or Bendersky constants. They emerge from studying the following product: ∏ m = 1 n m m k = 1 1 k ⋅ 2 2 k ⋅ 3 3 k ⋅ . . . ⋅ n n k {\displaystyle \prod _{m=1}^{n}m^{m^{k}}=1^{1^{k}}\cdot 2^{2^{k}}\cdot 3^{3^{k}}\cdot {...}\cdot n^{n^{k}}} Setting k = 0 {\displaystyle k=0} gives the factorial n ! {\displaystyle n!} , while choosing k = 1 {\displaystyle k=1} gives the hyperfactorial H ( n ) {\displaystyle H(n)} . Defining the following function P k ( n ) = ( n k + 1 k + 1 + n k 2 + B k + 1 k + 1 ) ln n − n k + 1 ( k + 1 ) 2 + k ! ∑ j = 1 k − 1 B j + 1 ( j + 1 ) ! n k − j ( k − j ) ! ( ln n + ∑ i = 1 j 1 k − i + 1 ) {\displaystyle P_{k}(n)=\left({\frac {n^{k+1}}{k+1}