is a somewhat typical representative of even zeta values — the values of Riemann’s zeta function at positive even integers. It is shadowed by the far more famous , which was a main subject of Euler’s resolution of Basel’s problem, and — an objet de l’étude of Apéry’s iconic proof of the irrationality of the latter (and also of ) [4, 30]. Though known to be irrational (and transcendental!), serves as a natural guinea pig for extending Apéry’s machinery to other zeta values. Apéry-type approximations to the number were discovered and rediscovered on several occasions [8, 29, 33], however they were not good enough to draw conclusions about its irrationality. An unexpected difficulty to control the ‘true’ arithmetic of those rational approximations to generated further research [14, 34], which eventually led to producing sufficient approximations and establishing a new world record for the irrationality measure of .
In this note we turn our attention to a rational side of the coin and prove the following two-parametric identity.
For integers , define two rational functions
The fact that both sides of (1) are linear forms in and with rational coefficients is verifiable by standard techniques [14, 33, 34] which employ the partial-fraction decomposition of the rational functions. A remarkable outcome of this identity is the coincidence of two different-looking rational approximations to the zeta value. Such coincidences are often a source of deep algorithmic and analytical developments — check  for another exploration of this theme (see also ).
The main difficulty in establishing equality (1) (in contrast to tackling, for example, Apéry’s sums in  for ) is that its both sides are not hypergeometric functions but rather derivatives of hypergeometric functions. Another issue is that the summation range on the left-hand side is somewhat unnatural.
2. Symbolic summation
(A) We compute the linear recurrence
which holds simultaneously for and for all with .
(B) We show that the following initial values hold:
Combined with (A) this proves that holds true for all .
In order to carry out the steps (A) and (B), advanced symbolic summation techniques in the setting of difference rings are utilized. Among them the following three summation paradigms play a decisive role, that are available within the summation package Sigma .
(i) Creative telescoping. Given a sum and , one searches for constants , free of , and such that
holds for all . Thus summing (6) over one obtains the recurrence
By specializing further — e.g., to and , or sending to if the limit exists — one obtains recurrence relations for more specific sums. The computed creative telescoping solution is also called a proof certificate for the recurrence (7) found: usually it allows one to verify that is a solution of (7) by simple polynomial arithmetic, without analyzing the usually complicated computation steps of the underlying summation algorithm. The algorithmic version of creative telescoping has been introduced in [32, 18] for hypergeometric sums. In order to prove (1), we will employ a generalized machinery for creative telescoping  where the summand can be composed not only in terms of hypergeometric products, but over indefinite nested sums defined over hypergeometric products. We emphasize that all recurrences produced below (using the Sigma-command GenerateRecurrence) are accompanied by such proof certificates which guarantee the correctness of all the calculations. Since the output is rather large and can be easily reproduced with Sigma, any explicit printout of the proof certificates is skipped.
(ii) Recurrence solving. Given a linear recurrence of the form (7), one can search for solutions that are expressible within certain classes function spaces. Using the Sigma-command SolveRecurrence one can search for hypergeometric solutions [17, 18] and, more generally, for all solutions that are expressible in terms of indefinite nested sums defined over hypergeometric products. Such solutions are also called d’Alembertian solutions [2, 19] a subclass of Liouvillian solutions .
(iii) Simplification of expressions. Within Sigma the expressions in terms of indefinite nested sums defined over hypergeometric products are represented in the setting of difference rings and fields [12, 23, 27]. Utilizing this difference ring machinery [24, 28] (compare also ) one can apply, e.g., the Sigma-command SigmaReduce to an expression in terms of indefinite nested sums. Then the output is a simplified expression where the arising sums and products (except products such as ) are independent among each other as functions of their external parameter. In particular, the input expression evaluates to zero (from a certain point on) if and only if Sigma reduces the expression to the zero-expression.
These summation paradigms can be used to transform a definite (multi-)sum to an expression in terms of indefinite nested sums by deriving a linear recurrence, solving the recurrence found in terms of indefinite nested sums, and, in case that sufficiently many solutions are found, combining them to an expression that evaluates to the same sequence as the input sum. Recently this machinery has been used for large scale problems coming from particle physics (see, e.g.,  and references therein). In this regard, also the package EvaluateMultiSum, which automatizes this summation mechanism, has been utilized non-trivially in the sections below.
3. A linear recurrence in for the left-hand side
By definition all Pochhammer symbols in the former expression are of the form for some and . Thus, we can apply the well-known formula
for with , where
denotes the harmonic number of order . Employing this formula we get for all the following representation:
In addition, we will prepare the summand of . Here we note that the rule (9) cannot be applied to the arising factor , however we can easily overcome this issue by using the following: For with and any differentiable function , we have
Therefore, for all with we get
Because of the factor , we have for all with . In a nutshell, and can be written as
where the summands and are given in terms of hypergeometric products and linear combinations of harmonic numbers. Since these sums fit to the input class of Sigma, we can apply the command GenerateRecurrence to both sums and compute the recurrences
where the coefficients are given in (3) and where is too large to be reproduced here (verification of the latter equality required an extra simplification step with Sigma). To compute the recurrence for the hypergeometric sum one can alternatively use the Mathematica package fastZeil  based on .
4. A linear recurrence in for the right-hand side
In order to calculate a linear recurrence for we will follow the same strategy as for in Section 3 by utilizing more advanced summation tools of Sigma. Define
As before, we apply the formula (9) and its relatives to get the following monster summand of in terms of the harmonic numbers of order and :
In order to tackle the summand of , we have to differentiate twice. With and
we conclude that for all we have
the last equality follows since for all . In addition, we can use the following calculation: For and , we have
For , we use and the rule
valid for any differentiable function , in place of (9). It follows that
Denote by , , , and the four resulting sums in (11) and use Sigma to compute four linear recurrences of with . A routine calculation demonstrates that each of the recurrences found can be brought to the form
where the coefficients are given in (3) and where only the right-hand sides for differ. As an illustration, we provide with details about how we treat