Apéry's constant is defined by
|
(1)
|
(OEIS A002117) where is the Riemann zeta
function. Apéry (1979) proved that
is irrational,
although it is not known if it is transcendental.
Sorokin (1994) and Nesterenko (1996) subsequently constructed independent proofs
for the irrationality of
(Hata 2000). Apéry's proof
involved use of the continued fraction
|
(2)
|
(Raayoni 2021, Elimelech et al. 2023).
Lai (2026) reported that the numbers 1, , and
are linearly independent
over the rationals. Since
, this implies that
is irrational.
Lai (2026) states that ChatGPT GPT-6 Astra generated the proof
and that he verified and rewrote it, taking responsibility for its correctness.
arises naturally in a number of
physical problems, including in the second- and third-order terms of the electron's
gyromagnetic ratio, computed using quantum electrodynamics.
Leibtag et al. (2026) obtained the improved irrationality measure bound
|
(3)
|
This improves the earlier bound 5.513891 of Rhin and Viola (2001). They also gave an independent proof that 1, , and
are linearly independent,
with
for every nonzero pair of rational
numbers
. Leibtag et al. (2026) state that their
main independence construction was developed without an LLM. They report later AI
assistance with technical arithmetic estimates used for the improved bounds and with
formalization.
The following table summarizes progress in computing upper bounds on the irrationality measure
for . Here, the exact values for
is given by
|
(4)
| |||
|
(5)
|
(Hata 2000).
| upper bound | reference | |
| 1 | 5.513891 | Rhin and Viola (2001) |
| 1 | 5.4197 | Leibtag et al. (2026) |
| 2 | 8.830284 | Hata (1990) |
| 3 | 12.74359 | Dvornicich and Viola (1987) |
| 4 | 13.41782 | Apéry (1979), Sorokin (1994), Nesterenko (1996), Prévost (1996) |
Beukers (1979) reproduced Apéry's rational approximation to
using the triple integral of
the form
|
(6)
|
where
is a Legendre polynomial. Beukers's integral
is given by
|
(7)
|
a result that is a special case of what is known as Hadjicostas's formula.
This integral is closely related to using the curious identity
|
(8)
| |||
|
(9)
|
where
is a generalized harmonic number and
is a polygamma function
(Hata 2000).
Sums related to include
|
(10)
| |||
|
(11)
|
(used by Apéry), the related sum
|
(12)
|
as first proved by G. Huvent in 2002 (Gourevitch) and rediscovered by B. Cloitre (pers. comm., Oct. 8, 2004), and
|
(13)
| |||
|
(14)
| |||
|
(15)
| |||
|
(16)
| |||
|
(17)
|
where
is the Dirichlet lambda function. The
above equations are special cases of a general result due to Ramanujan (Berndt 1985).
Apéry's constant is given by an infinite family BBP-type formulas of the form
|
(18)
| |||
|
(19)
| |||
|
(20)
| |||
|
(21)
| |||
|
(22)
| |||
|
(23)
| |||
|
(24)
|
(E. W. Weisstein, Feb. 25, 2006), and the amazing two special sums
|
(25)
| |||
|
(26)
|
Determining a sum of this type is given as an exercise by Bailey et al. (2007, p. 225; typo corrected).
Rapidly converging series include
|
(27)
|
(Berndt 1988, Plouffe 1998) and
|
(28)
|
(Plouffe 2006).
|
(29)
|
(Apostol 1973, given incorrectly in Stark 1974), where tends to infinity through positive integers.
A beautiful double series for is given by
|
(30)
|
where
is a harmonic number (O. Oloa, pers. comm.,
Dec. 30, 2005).
Apéry's proof relied on showing that the sum
|
(31)
|
where
is a binomial coefficient, satisfies the
recurrence relation
|
(32)
|
(van der Poorten 1979, Zeilberger 1991). The characteristic polynomial
has roots
, so
|
(33)
|
is irrational and cannot satisfy a two-term recurrence
relation (Jin and Dickinson 2000).
Apéry's constant is also given by
|
(34)
|
where
is a Stirling number of the first kind.
This can be rewritten as
|
(35)
| |||
|
(36)
|
where
is the
th
harmonic number (Castellanos 1988).
Amdeberhan (1996) used Wilf-Zeilberger pairs with
|
(37)
|
to obtain
|
(38)
|
For ,
|
(39)
|
(Boros and Moll 2004, p. 236; Amdeberhan 1996), and for ,
|
(40)
|
(Amdeberhan 1996). The corresponding for
and 2 are
|
(41)
|
and
|
(42)
|
Amdeberhan and Zeilberger (1997) used a Wilf-Zeilberger pair identity with
|
(43)
|
, and
, to obtain the rapidly converging
series
|
(44)
|
which was used to compute to 1 million decimal digits.
Campbell (2022) used the WZ method to obtain
|
(45)
|
Zuniga (2023, p. 7) gives the rapidly converging series
|
(46)
|
where
|
(47)
|
and
|
(48)
|
The ratio of successive terms tends to , corresponding to about 11.86 decimal
digits per term.
Integrals for include
|
(49)
| |||
|
(50)
|
Gosper (1990) gave
|
(51)
|
A continued fraction involving Apéry's constant is
|
(52)
|
(Apéry 1979, Le Lionnais 1983).
The derivative of the Riemann zeta function satisfies
|
(53)
|
is related to the Glaisher-Kinkelin
constant
and polygamma function
by
|
(54)
|
Gosper (1996) expressed as the matrix product
|
(55)
|
where
|
(56)
|
which gives 12 bits per term. The first few terms are
|
(57)
| |||
|
(58)
| |||
|
(59)
|
which gives
|
(60)
|
Given three integers chosen at random, the probability that no common factor will divide them all is
|
(61)
|