Preprint
Article

This version is not peer-reviewed.

Harmonic Number Series Associated with Certain Generating Functions

Submitted:

14 June 2025

Posted:

19 June 2025

You are already at the latest version

Abstract
We present several series involving central binomial coefficients, Catalan numbers $C_n$, harmonic numbers $H_n$, odd harmonic numbers $O_n$, and certain products of these famous sequences. A highlight of our analysis is the derivation of closed-form expressions for the two classes $$\sum_{n=1}^\infty \binom {2n}{n} H_nO_n x^n\quad \text{and} \quad \sum_{n=1}^\infty C_nH_nO_nx^n.$$ Our results complement the findings of Lehmer, Boyadzhiev, Chen, and Li and Chu, among others. We will also derive new series involving products of harmonic numbers. Some of our series are nontrivial in the sense that the series cannot be expressed in terms of elementary functions but will involve combinations of dilogarithms. Only in exceptional cases, where the golden ratio is involved, can we express them in simple closed forms.
Keywords: 
;  ;  ;  ;  

1. Introduction and Motivation

The central binomial coefficients are defined for nonnegative n by 2 n n = ( 2 n ) ! ( n ! ) 2 . These numbers are closely related to the famous Catalan numbers C n = 1 n + 1 2 n n , which can also be expressed by the recursion
C n = 2 ( 2 n 1 ) n + 1 C n 1 , C 0 = 1 .
Catalan numbers have a long history, possess beautiful properties, and play an outstanding role in combinatorics. Excellent sources on these numbers are the books by Koshy [12], Roman [22] and Stanley [26]. Some examples of recent work involving Catalan numbers include [2,3,4,8,10,20,21,29]. The ordinary generating functions of the numbers 2 n n and C n are well known [14]:
n = 0 2 n n x n = 1 1 4 x and n = 0 C n x n = 2 1 + 1 4 x .
Both series converge for | x | < 1 / 4 .
We proceed with some definitions and basic statements. Let Li n ( x ) and ζ ( s ) denote the polylogarithm functions and the Riemann zeta function, respectively. These functions are defined by [15,25]
Li 0 ( x ) = x 1 x , Li 1 ( x ) = ln ( 1 x ) , and Li n ( x ) = k = 1 x k k n , | x | 1 , n 2 ,
and
ζ ( s ) = k = 1 1 k s = 1 1 2 s k = 1 1 ( 2 k 1 ) s , ( s ) > 1 , 1 1 2 1 s k = 1 ( 1 ) k 1 k s , ( s ) > 0 ; s 1 .
For n 2 , we have the identities Li n ( 1 ) = ζ ( n ) and Li n ( 1 ) = ( 1 2 1 n ) ζ ( n ) . The function Li 2 ( x ) is the familiar dilogarithm, and we have
Li 2 1 2 = π 2 12 1 2 ln 2 2 .
Among the many beautiful properties of the dilogarithm, the following identity will be used frequently in the paper:
Li 2 ( y ) = Li 2 1 y + π 2 3 1 2 ln 2 y i π ln y , y > 1 , i = 1 .
Harmonic numbers of order m, denoted H n ( m ) , and odd harmonic numbers of order m, denoted O n ( m ) , are defined by
H n ( m ) = k = 1 n 1 k m , O n ( m ) = k = 1 n 1 ( 2 k 1 ) m ,
where H 0 ( m ) = O 0 ( m ) = 0 , and H n ( 1 ) = H n and O n ( 1 ) = O n are the ordinary harmonic and odd harmonic numbers, respectively. Harmonic numbers and odd harmonic numbers are related, among other ways, by
H 2 n = 1 2 H n + O n and H 2 n 1 = 1 2 H n 1 + O n .
A fundamental connection between harmonic numbers and the digamma function ψ ( z ) = Γ ( z ) / Γ ( z ) is given by [25]
H z = ψ ( z + 1 ) + γ ,
with γ being the Euler–Mascheroni constant and ψ ( z ) equals
ψ ( z ) = γ + k = 0 1 k + 1 1 k + z .
Similarly, generalized harmonic numbers H n ( m ) are connected to the polygamma functions ψ ( m ) ( z ) of order m, defined by
ψ ( m ) ( z ) = d m d z m ψ ( z ) = ( 1 ) m + 1 m ! j = 0 1 ( j + z ) m + 1 ,
through the relation
H z ( m ) = ζ ( m ) + ( 1 ) m 1 ( m 1 ) ! ψ ( m 1 ) ( z + 1 ) .
These connections imply that generalized harmonic numbers at rational arguments H p / q ( m ) must be interpreted as polygamma functions evaluated at rational arguments. As shown by Kölbig [13], and later by Choi and Cvijović [7], the polygamma function ψ ( m 1 ) ( p / q + 1 ) can be explicitly expressed in terms of polylogarithms or other special functions. For instance, we have
H 3 / 4 = π 2 3 ln 2 and H 1 / 6 ( 3 ) = 2 3 π 3 90 ζ ( 3 ) .
The correspondence between H z ( m ) and ψ ( m ) ( z ) also enables the derivation of further relations between harmonic numbers and odd harmonic numbers. Two such relations relevant to this study are
H n 1 / 2 = 2 O n 2 ln 2 and H n 1 / 2 H 1 / 2 = 2 O n .
This paper is particularly motivated by the papers by Boyadzhiev [5] and Chen [6]. Boyadzhiev studied series (generating functions) involving central binomial coefficients and harmonic numbers H n . His main results include ([5], Theorem 1)
n = 0 2 n n H n x n = 2 1 4 x ln 1 + 1 4 x 2 1 4 x , x [ 1 / 4 ; 1 / 4 ) ,
and the generating function for the product C n H n , valid for | x | 1 / 4 ([5], Corollary 2)
n = 0 C n H n x n = 1 x 1 4 x ln 2 1 4 x 1 + 1 4 x ln 1 + 1 4 x + ln 2 .
Chen added results for the sequences 2 n n ( H 2 n H n ) , C n ( H 2 n H n ) , and 2 n n O n . Among other results, he proved the identities ([6], Theorem 8)
n = 1 2 n n O n x n = ln ( 1 4 x ) 2 1 4 x , x [ 1 / 4 ; 1 / 4 ) ,
as well as ([6], Corollary 9)
n = 1 C n O n x n = 1 2 x 1 1 4 x + 1 2 1 4 x ln ( 1 4 x ) , | x | < 1 / 4 .
See also [3] for alternative proofs of Chen’s results, as well as generalizations and related combinatorial identities.
The goal of the present paper is twofold. First, we derive and evaluate in closed form a range of harmonic number series associated with certain known generating functions. Second, we establish new expressions for two classes of series involving products of harmonic numbers H n and odd harmonic numbers O n , respectively. These classes are
n = 1 2 n n H n O n x n and n = 1 C n H n O n x n .
The method employed to prove these expressions is well known and is based on integration. The evaluation of infinite series via integrals is a classical analytical technique, recently revitalized in works by Xu [28], Sofo and Nimbran [24], Stewart [27], and Li and Chu [16,17,18,19]. We also derive new series involving products of harmonic numbers.
To provide the reader with a sense of the results to come, we present the following example:
n = 1 H n H n + r + s + 1 n + r + s + 1 r + 2 = r + 2 s r + s r ( H r + s H r ) H r + s + 1 r + 1 + H r + s ( 2 ) H r ( 2 ) ,
which holds for 0 r C Z , 0 s C Z . A special case of this identity is
n = 1 H n H n + 1 n ( n + 1 ) = ζ ( 2 ) + 2 ζ ( 3 ) ,
where ζ ( 2 ) = π 2 / 6 . The second example deals with two series involving the product H n O n . In general, these series exhibit a nontrivial structure, in the sense that they typically cannot be expressed in terms of elementary functions, but rather involve combinations of dilogarithms. Two such evaluations are
n = 1 ( 1 ) n 2 n n H n O n 8 n = 2 3 Li 2 2 3 Li 2 2 3 π 2 4 1 2 ln 2 3 ln 12
and
n = 1 ( 1 ) n C n H n O n 4 n = 2 Li 2 ( 2 ) Li 2 ( 2 ) + 2 2 Li 2 1 2 Li 2 1 2 + 2 2 + 2 2 ln 2 ln ( 2 + 1 ) π 2 2 ( 1 + 2 ) + 3 2 ln 2 6 2 4 ln 2 + i π ln 2 .
The latter result is a symmetric expression in dilogarithms and an equivalent form is
n = 1 ( 1 ) n 1 C n H n O n 4 n = 2 2 ( 2 1 ) Li 2 2 2 + Li 2 2 2 2 ( 2 + 2 2 ln 2 ) ln ( 2 + 1 ) + 6 2 7 12 π 2 6 2 1 2 ln 2 2 + 2 ( 3 2 + 2 ) ln 2 .
However, we succeed in evaluating certain series involving the golden ratio α = ( 1 + 5 ) / 2 in closed form. For instance, in Corollary 6 below we show that
n = 1 2 n n H n O n 4 n α n = π 2 12 α + 1 2 ( 5 ln α 4 ln 2 ) α ln α .
The next lemma will be used frequently in this paper.
Lemma 1. 
Let ( a n ) n 1 be a sequence with generating function f ( x ) = n = 1 a n x n . Then we have
n = 1 a n H n x n = 0 1 f ( x ) f ( x t ) 1 t d t
and
n = 1 a n O n x n = 0 1 f ( x ) f ( x t 2 ) 1 t 2 d t .
Proof. 
The results follow directly from the representations
H n = 0 1 1 x n 1 x d x and O n = 0 1 1 x 2 n 1 x 2 d x .
As a simple application of Lemma 1, we obtain the generating function of the harmonic numbers. By choosing a n = 1 for all n and employing the geometric series f ( x ) = x / ( 1 x ) , where | x | < 1 , we have
n = 1 H n x n = 0 1 x 1 x x t 1 x t 1 t d t = x 1 x 0 1 1 1 x t d t = ln ( 1 x ) 1 x , | x | < 1 .

2. Series Associated with Known Generating Functions

In this section, we derive some interesting series identities from known generating functions, hereby focusing on the generating functions (2), (14), (10), and (11).

2.1. Series Associated with Generating Function (2)

First, we state a lemma.
Lemma 2. 
For u , v C Z , we have
0 1 x u ( 1 x ) v d x = 1 u + v + 1 u + 1 ( u + 1 ) ,
0 1 x u ( 1 x ) v ln x d x = H u + v + 1 H u u + v + 1 u + 1 ( u + 1 ) ,
0 1 x u ( 1 x ) v ln 2 x d x = ( H u + v 1 H u ) 2 + H u + v + 1 ( 2 ) H u ( 2 ) u + v + 1 v + 1 ( v + 1 ) ,
and, in particular,
0 1 x u 1 ln ( 1 x ) d x = H u u , 0 u C Z .
Proof. 
The identity (15) is the well-known Beta function integral, while (16) is obtained by differentiating with respect to u. Identity (17) follows from differentiating (16) with respect to u. These identities are also listed in Gradshteyn and Ryzhik’s book [11] and other tables. □
Theorem 1. 
For s 1 / 2 Z , r 0 and r + s 1 / 2 Z , we have
n = 0 2 n n 2 2 n H n + r + s H s n + r + s s + 1 = 2 ( s + 1 ) 2 s + 1 H r + s 1 / 2 H s 1 / 2 r + s 1 / 2 r 1 .
In particular,
n = 0 2 n n 2 2 n H n + r n + r = 2 2 r + 1 2 r r O r r , r 0 ,
and
n = 0 C n H n + 1 2 2 n = 4 .
Proof. 
Write ( 1 x ) / 4 for x into the first identity of (2) and arrange to obtain
n = 0 2 n n 2 2 n ( 1 x ) n + r 1 x s ln x = ( 1 x ) r 1 x s 1 / 2 ln x .
Integrating from 0 to 1 using Lemma 2 results in the main identity (19). Identity (20) follows from identity (19) by setting s = 0 in conjunction with relation (7). □
Corollary 1. 
We have
n = 1 2 n n 2 2 n H n n = π 2 3 ,
Proof. 
Write (20) as
n = 1 2 n n 2 2 n H n + r n + r = 2 2 r + 1 2 r r O r r H r r , r 0 .
Thus,
n = 1 2 n n 2 2 n H n n = 2 lim r 0 O r r lim r 0 H r r .
Now, using L’Hôpital’s rule and (6), we obtain
lim r 0 O r r = lim r 0 H 2 r 1 2 H r r = lim r 0 2 ζ ( 2 ) H 2 r ( 2 ) 1 2 ζ ( 2 ) H r ( 2 ) 1 = 2 ζ ( 2 ) 1 2 ζ ( 2 ) = 3 2 ζ ( 2 )
and
lim r 0 H r r = lim r 0 ζ ( 2 ) H r ( 2 ) 1 = ζ ( 2 ) .
Plug (24) and (25) in (23) to obtain (22). □
Corollary 2. 
We have
n = 0 C n H n 2 2 n = 4 ln 2
and
n = 0 n n + 1 C n H n + 1 2 2 n = 2 π 2 3 4 .
Proof. 
Working with (21) we can calculate
n = 0 C n H n + 1 2 2 n = n = 0 C n H n 2 2 n + n = 0 C n 2 2 n ( n + 1 ) = n = 0 C n H n 2 2 n + 4 4 ln 2 ,
where we used the fact that
n = 0 C n 2 2 n ( n + 1 ) = 4 4 ln 2 .
This evaluation can be found in [1], for instance. The second identity is a combination of (21) with
n = 0 C n H n + 1 2 2 n ( n + 1 ) = 8 2 π 2 3 ,
which is a special case of a more general result also derived in [1]. □
Remark 1. 
Using (7) in conjunction with n 1 / 2 m = 2 2 m 2 n n n m 2 ( n m ) n m 1 , we can restate identity (19) in the form
n = 0 2 n n 2 2 n H n + r + s H s n + r + s s + 1 = s + 1 r 2 2 r + 1 2 s s O r + s O s r + s r 2 ( r + s ) r + s .
Theorem 2. 
We have
n = 0 O n + 1 ( n + 1 ) ( 2 n + 1 ) = π 2 6 ,
n = 0 C n O n + 3 2 2 n ( 2 n + 5 ) = 8 9 + π 32 π ln 2 4 ,
n = 0 C n H n + r + 1 2 2 n + 1 ( n + r + 1 ) = H r r 2 2 r ( O r + 1 1 ) 2 ( r 1 ) r 1 ( 2 r 1 ) ( 2 r + 1 ) , r 0 , r ± 1 2 ,
and more generally, for 0 s , r C Z , r + s 1 / 2 Z ,
n = 0 C n 2 2 n + 1 H n + r + s H s 1 n + r + s s = H r + s 1 H s 1 r + s 1 s 2 s 2 s + 1 H r + s 1 / 2 H s 1 / 2 r + s 1 / 2 r 1 .
Proof. 
Write ( 1 x ) / 4 for x in the second identity in (2) to obtain
n = 0 C n ( 1 x ) n + 1 2 2 n + 1 = 1 x 1 / 2 ,
which can now be arranged as
n = 0 C n ( 1 x ) n + r x s 1 2 2 n + 1 ln x = ( 1 x ) r 1 x s 1 ln x ( 1 x ) r 1 x s 1 / 2 ln x ,
which now produces (27) upon integration from 0 to 1, using Lemma 2. □
Remark 2. 
Identity (26) is also recorded as identity (8) and generalized in a different direction in [1].

2.2. Series Associated with Generating Function (14)

Next, we derive series associated with the generating function of harmonic numbers (14).
Theorem 3. 
We have
n = 1 H n ( n + s ) ( n + s + 1 ) = H s s , s 0 , n = 1 H n n ( n + 1 ) = π 2 6 , n = 1 2 2 n H n ( 2 n + 2 s 1 ) ( 2 n + 2 s + 1 ) ( 2 n + 2 s + 3 ) 2 ( n + s 1 ) n + s 1 = 4 3 s ( 2 s + 1 ) O s + 1 1 2 s s ,
and more generally, 0 r C Z , 0 s C Z ,
n = 1 H n n + r + s + 1 r + 2 = r + 2 s H r + s H r r + s r .
Proof. 
Write 1 x for x in the generating function (14) and arrange as
n = 1 H n x r + 1 ( 1 x ) n + s 1 = x r ( 1 x ) s 1 ln x .
Now integrate from 0 to 1, making use of Lemma 2; this gives (29). The other identities are special cases at r = s = 0 and r = 1 / 2 . □
Remark 3. 
The identity (28) looks very interesting but it is just a restatement of the integral given in (18). This integral is also stated by Xu ([28], Eq. (2.11)), with the restriction u > 0 . The equivalence of the forms follows from a simple calculation using (14):
H u u = 0 1 x u 1 ln ( 1 x ) d x = 0 1 x u 1 ( 1 x ) ln ( 1 x ) 1 x d x = n = 0 H n 0 1 x n + u 1 d x 0 1 x n + u d x .
Indeed, in Theorem 2.2, Xu [28] offers a formula for shifted harmonic series of the form n = 1 H n + u ( n + v ) ( n + w ) with parameters u, v, and w.
Theorem 4. 
For r C Z and m Z + ,
n = 1 H n n + r + 1 m + 1 = H r ζ ( m + 1 ) H r ( m + 1 ) + m + 1 2 ζ ( m + 2 ) H r ( m + 2 ) 1 2 n = 1 m 1 ζ ( n + 1 ) H r ( n + 1 ) ζ ( m n + 1 ) H r ( m n + 1 ) .
In particular,
n = 1 H n n + r + 1 2 = H r ζ ( 2 ) H r ( 2 ) + ζ ( 3 ) H r ( 3 ) .
Proof. 
Sum (28) over s from 1 to r to obtain
n = 1 H n 1 n + 1 1 n + r + 1 = 1 2 H r 2 + H r ( 2 ) ,
which is now valid for every complex number r that is not a negative integer. Identity (30) follows upon differentiating (31) with respect to r, m times, using Leibnitz rule and the facts that
d m d r m H r = ( 1 ) m 1 m ! ζ ( m + 1 ) H r ( m + 1 )
and
d m d r m H r ( 2 ) = ( 1 ) m 1 ( m + 1 ) ! ζ ( m + 2 ) H r ( m + 2 ) .
Remark 4. 
Identity (30) is equivalent to a main result of Sofo and Cvijović [23]. While Sofo’s identity is expressed in terms of the Hurwitz zeta function, our formulation is given directly in terms of the Riemann zeta function.
Theorem 5. 
We have
n = 1 H n H n + 1 n ( n + 1 ) = ζ ( 2 ) + 2 ζ ( 3 ) , n = 1 H n O n ( 2 n 3 ) ( 2 n 1 ) = 1 6 + π 2 36 + 1 2 ln 2 , n = 1 H n H n + s + 1 ( n + s ) ( n + s + 1 ) = 1 s H s + H s 2 + H s ( 2 ) , s 0 ,
and more generally, for 0 r C Z , 0 s C Z ,
n = 1 H n H n + r + s + 1 n + r + s + 1 r + 2 = r + 2 s r + s r ( H r + s H r ) H r + s + 1 r + 1 + H r + s ( 2 ) H r ( 2 ) .
Proof. 
Using (14), write
n = 1 H n ( 1 x ) n + s 1 x r + 1 ln x = ( 1 x ) s 1 x r ln 2 x
and integrate from 0 to 1 using Lemma 2. □

2.3. Series Associated with Generating Function (10)

Next, we establish some identities associated with the generating function (10).
Theorem 6. 
We have
n = 1 2 n n 2 ( n 1 ) n 1 O n 2 4 n n = ln 2 π ,
n = 1 2 n n 2 O n 2 4 n ( n + 1 ) = 4 ( 1 ln 2 ) π ,
and more generally, if r is a non-negative integer, then
n = 1 2 n n 2 ( n + r 1 ) n + r 1 O n 2 4 n ( n + r ) = 2 2 r 2 r 1 O r ln 2 π .
Proof. 
Write x / 4 for x in (10) and arrange as
n = 1 2 n n O n 2 2 n x n + r 1 ( 1 x ) 1 / 2 = 1 2 x r 1 ln ( 1 x ) .
Integrate both sides from 0 to 1 using Lemma 2 and substitute r 1 / 2 for r to obtain
n = 1 2 n n O n 2 2 n n + r 3 / 2 = 3 H r 1 / 2 2 ( 2 r 1 )
and hence deduce (32), since
q 3 2 = q ! ( 3 2 ) ! q 1 1 2 ! = q ! 3 π 4 ! 2 2 ( q 1 ) π 2 ( q 1 ) q 1 ( q 1 ) ! = 2 2 q q 3 π 2 ( q 1 ) q 1
and using the identity H k 1 / 2 = 2 O k 2 ln 2 . □
Theorem 7. 
We have
n = 1 2 n n O n 2 2 n n = n = 1 2 2 n H n n ( n + 1 ) 2 ( n + 1 ) n + 1 , n = 1 C n O n 2 2 n = n = 1 2 2 ( n + 1 ) H n ( n + 1 ) ( n + 2 ) 2 ( n + 2 ) n + 2 ,
and more generally, if r is a non-negative integer, then
n = 1 2 n n O n 2 2 n ( n + r ) = n = 1 2 2 ( n + r ) H n ( n + r ) ( n + r + 1 ) 2 ( n + r + 1 ) n + r + 1 .
Proof. 
We have
n = 1 2 n n 2 2 n O n x n + r 1 = x r 1 2 ln ( 1 x ) 1 x = x r 1 1 x 2 n = 1 x n H n = 1 2 n = 1 x n + r 1 ( 1 x ) 1 / 2 H n ;
from which, integrating from 0 to 1, using Lemma 2, we get
n = 1 2 n n 2 2 n O n n + r = 1 3 n = 1 H n n + r + 1 / 2 3 / 2
and hence (33), since
q + 1 / 2 3 / 2 = q + 1 / 2 q 1 = π 2 2 q + 2 2 q + 2 q + 1 ( q + 1 ) ! 3 π 4 q 1 ! = q ( q + 1 ) 2 2 q 3 2 q + 2 q + 1
and ( q + 1 / 2 ) ! = π 2 2 q + 2 2 q + 2 q + 1 ( q + 1 ) ! .
Remark 5. 
Combining the previous results, we obtain
n = 1 C n O n 2 2 n n = 3 2 n = 1 2 2 n n ( n + 1 ) ( 2 n + 1 ) ( 2 n + 3 ) H n C n .
Next, we present selected series evaluations involving products of Fibonacci (Lucas) numbers and odd harmonic numbers. Recall that the Fibonacci and Lucas numbers both satisfy the recurrence u n = u n 1 + u n 2 , n 2 , but with respective initial conditions F 0 = 0 , F 1 = 1 and L 0 = 2 , L 1 = 1 . For negative subscripts, we have F n = ( 1 ) n 1 F n and L n = ( 1 ) n L n .
The Binet formulas for these numbers state that
F n = α n β n α β , L n = α n + β n , n Z ,
where α = ( 1 + 5 ) / 2 and β = 1 / α .
Theorem 8. 
If r is an even integer and s any integer, then
n = 1 ( 1 ) n 1 2 n n 2 2 n 1 O n F r n + s L r n = 1 5 α s ln 1 + α r / L r 1 + α r / L r β s ln 1 + β r / L r 1 + β r / L r
and
n = 1 ( 1 ) n 1 2 n n 2 2 n 1 O n L r n + s L r n = α s ln 1 + α r / L r 1 + α r / L r + β s ln 1 + β r / L r 1 + β r / L r .
Proof. 
Insert x = α r / L r and x = β r / L r , in turn, into (10), then multiply through by α s (or β s ) and combine the resulting expressions according to the Binet formulas (34). □
Example 1. 
Choosing s = 0 and r = 0 , r = 1 , r = 2 in Theorem 8 gives:
n = 1 ( 1 ) n 1 O n 8 n = 1 6 ln 3 2 , n = 1 ( 1 ) n 1 O n L 2 n 12 n = 57 38 5 α ln α + 4 3 + 4 + α ln 5 α 3 , n = 1 ( 1 ) n 1 O n F 2 n 12 n = 285 190 5 α ln α + 4 3 4 + α ln 5 α 3 , n = 1 ( 1 ) n 1 O n L 4 n 28 n = 231 66 4 α ln 3 α + 9 7 + 3 + α ln 12 3 α 7 , n = 1 ( 1 ) n 1 O n F 4 n 28 n = 1155 330 4 α ln 3 α + 9 7 3 + α ln 12 3 α 7 .

2.4. Series Associated with Generating Function (11)

We now derive identities associated with generating function (11).
Theorem 9. 
We have
n = 1 C n O n H n + 1 2 2 n ( n + 1 ) = 24 14 ζ ( 3 ) 2 π 2 3 ,
n = 1 C n O n H n + 2 2 2 n ( n + 2 ) = 14 27 ,
and more generally, for s 1 / 2 Z , 0 s C Z , r + s + 1 / 2 Z ,
n = 1 C n O n 2 2 n H n + r + s + 1 H s 1 n + r + s + 1 s = 2 H r + s H s 1 r + s s 4 s 2 s + 1 H r + s + 1 / 2 H s 1 / 2 r + s + 1 / 2 r 2 s 2 s + 1 H r + s + 1 / 2 H s 1 / 2 2 + H r + s + 1 / 2 ( 2 ) H s 1 / 2 ( 2 ) r + s + 1 / 2 r .
Proof. 
Write ( 1 x ) / 4 for x in (11) and rearrange to obtain
n = 1 C n O n 2 2 n ( 1 x ) n + r + 1 x s 1 = 2 ( 1 x ) r x s 1 2 ( 1 x ) r x s 1 / 2 + ( 1 x ) r x s 1 / 2 ln x .
Now integrate both sides from 0 to 1 using Lemma 2. We also note that identity (35) was obtained by setting s = 1 in (37) and taking the limit at r = 1 . □

3. Series Associated with the Generating Function of Generalized Harmonic Numbers

In ([28], Lemma 1.2), Xu established the following result.
Lemma 3. 
For 0 s C Z and p Z + ,
0 1 x s 1 Li p ( x ) d x = n = 1 p 1 ( 1 ) n 1 s n ζ ( p + 1 n ) ( 1 ) p H s s p .
Theorem 10. 
For 0 s C Z and p Z + ,
n = 1 H n ( p ) n + s n + s + 1 = n = 1 p 1 ( 1 ) n 1 s n ζ ( p + 1 n ) ( 1 ) p H s s p .
Proof. 
The generating function of the generalized harmonic numbers is given by
n = 1 H n ( p ) x n = Li p ( x ) 1 x ,
which allows the following rearrangement:
n = 1 H n ( p ) ( 1 x ) x n + s 1 = x s 1 Li p ( x ) ,
from which (39) follows by integration, using (38). □
Remark 6. 
Identity (39) generalizes identity (28).
Theorem 11. 
If p is a positive integer and m and r are non-negative integers, then
n = 1 H n + r ( m + 1 ) H n ( m + 1 ) n p = n = 1 p 1 ( 1 ) n 1 m + 1 n m ζ p + 1 n H r ( m + n ) + ( 1 ) p n = 1 m p + m n 1 m n ζ ( n + 1 ) H r ( m + p n ) ( 1 ) p n = 0 m j = 1 r p + m n 1 m n H j ( n + 1 ) j p + m n .
In particular,
n = 1 H n + r n 2 = 2 ζ ( 3 ) + ζ ( 2 ) H r n = 1 r H n n 2 ,
where we used
n = 1 H n n 2 = 2 ζ ( 3 ) .
Proof. 
We begin with the polylogarithm series
n = 1 x n n p = Li p ( x ) ,
which allows us to write
n = 1 x n + s 1 n p = x s 1 Li p ( x ) ,
from which, integrating from 0 to 1, we get, by (38),
n = 1 1 n p n + s = n = 1 p 1 ( 1 ) n 1 s n ζ ( p + 1 n ) ( 1 ) p H s s p .
Differentiating (41) m times with respect to s and simplifying, we find
n = 1 1 n p n + s m + 1 = n = 1 p 1 ( 1 ) n 1 m + 1 n m ζ ( p + 1 n ) s m + n + ( 1 ) p n = 1 m p + m n 1 m n ζ ( n + 1 ) s p + m n ( 1 ) p n = 0 m p + m n 1 m n H s ( n + 1 ) s p + m n .
Summing (42) over s from 1 to r gives (40). □
Corollary 3. 
If p 2 is an integer and r is a non-negative integer, then
n = 1 H n + r n p = p + 2 2 ζ ( p + 1 ) 1 2 n = 2 p 1 ζ ( n ) ζ ( p + 1 n ) + n = 1 p 1 ( 1 ) n 1 ζ p + 1 n H r ( n ) ( 1 ) p n = 1 r H n n p .
Proof. 
Set m = 0 in (40) and use the standard Euler sum
n = 1 H n n p = p + 2 2 ζ ( p + 1 ) 1 2 n = 2 p 1 ζ ( n ) ζ ( p + 1 n ) .
Remark 7. 
An explicit evaluation of each of the following sums is possible:
n = 1 H n + r ( p ) n p , n = 1 H n + r ( m + 1 ) n p ( when m + p is even integer ) , n = 1 H n + r ( 2 ) n 4 , and n = 1 H n + r ( 4 ) n 2 .
This is due to the well-known fact that the sum n = 1 H n ( q ) n p can be expressed in terms of zeta values in the special cases q = 1 , p = q or p + q odd or p + q even, but with ( p , q ) being restricted to ( 2 , 4 ) or ( 4 , 2 ) .

4. Series Associated with the Dilogarithm

We require the integration formulas stated in the next lemma.
Lemma 4. 
For s > 0 , we have
0 1 x s 1 Li 2 x d x = π 2 6 s H s s 2 ,
0 1 ( 1 x ) s 1 Li 2 x d x = π 2 6 s H s ( 2 ) s ,
0 1 x s 1 Li 2 ( x ) ln x d x = π 2 3 s 2 + H s ( 2 ) s 2 + 2 H s s 3 ,
and
0 1 1 x s 1 Li 2 ( x ) ln 1 x d x = π 2 6 s 2 2 ζ ( 3 ) s + 2 H s 3 s + H s ( 2 ) s 2 .
Proof. 
Identity (43) is recorded in Lewin ([15], p. 308) (also p = 2 in (38)).
Multiply through the dilogarithm functional equation
Li 2 ( x ) + Li 2 1 x = π 2 6 ln x ln 1 x ,
by x u 1 and integrate from 0 to 1, using (43) and the known integral (e.g., Xu [28], Eq. (2.13)):
0 1 x s 1 ln x ln 1 x d x = π 2 6 s + H s ( 2 ) s + H s s 2 ;
this gives (44). Identities (45) and (46) are obtained by differentiating (43) and (44) with respect to u. □
Theorem 12. 
We have
n = 1 H n + 1 ( 2 ) n n + 1 = 3 π 2 6 .
Proof. 
Write
n = 1 1 x n n Li 2 ( x ) = ln x Li 2 ( x )
and integrate. □
Theorem 13. 
If 0 s C Z , then
n = 1 n + s n + 1 1 H n + s H s 1 ( n + 1 ) n 2 = ζ ( 2 ) s 2 + 2 ζ ( 3 ) s 2 H s ( 3 ) s H s ( 2 ) s 2 .
In particular,
n = 1 H n + 1 n 2 n + 1 = ζ ( 2 ) + 2 ζ ( 3 ) 3 .
Proof. 
The series
n = 1 x n n 2 = Li 2 ( x )
allows us to write
n = 1 1 x s 1 x n n 2 ln 1 x = 1 x s 1 Li 2 ( x ) ln 1 x ,
from which (47) follows upon term-wise integration from 0 to 1, using (16) and (46). □
Corollary 4. 
We have
n = 1 2 2 n 2 n + 1 n + 1 1 O n + 1 n 2 n + 1 = 3 2 ζ ( 2 ) + 7 2 ζ ( 3 ) 6 .
Proof. 
Set x = 1 / 2 in (47) and use n + 1 / 2 n + 1 = 2 2 ( n + 1 ) 2 n + 1 n + 1 , along with the relations H n + 1 / 2 H 1 / 2 = 2 O n + 1 , H 1 / 2 ( 3 ) = 8 6 ζ ( 3 ) , and H 1 / 2 ( 2 ) = 4 2 ζ ( 2 ) . □
The next two results are derived from the following generating function of H n 2 , which is given by (see, for example, [9]):
n = 1 H n 2 x n = ln 2 1 x 1 x + Li 2 ( x ) 1 x .
Theorem 14. 
For 0 s C Z ,
n = 1 H n 2 n + s n + s + 1 = H s H s 1 s + H s ( 2 ) s + ζ ( 2 ) s .
In particular,
n = 1 H n 2 ( n + 1 ) ( n + 2 ) = ζ ( 2 ) + 1 .
Proof. 
Using (48), set up the following identity
n = 1 H n 2 x n + s 1 1 x = x s 1 ln 2 1 x + x s 1 Li 2 ( x )
and integrate, using (17) and Lemma 4 to evaluate the resulting integrals. □
Theorem 15. 
For 0 s C Z ,
n = 1 H n 2 n + s + 1 s + 1 = s + 1 s 2 s 2 + ζ ( 2 ) H s ( 2 ) .
Proof. 
Write
n = 1 H n 2 x n 1 x s = 1 x s 1 ln 2 1 x + 1 x s 1 Li 2 ( x ) ,
and integrate. □

5. A First New Generating Function

The goal of this section is to derive a new generating function involving central binomial coefficients, harmonic numbers, and odd harmonic numbers. Our first result in this section is a new proof of the generating function (10) (see ([6], Theorem 8) or ([3], Theorem 3)).
Theorem 16. 
For all x [ 1 / 4 , 1 / 4 ) , we have
n = 1 2 n n O n x n = ln ( 1 4 x ) 2 1 4 x .
Proof. 
For the sequence a n = 2 n n , we use (2) to see that f ( x ) = 1 1 4 x 1 4 x and by series-integral correspondence (13), we have
n = 1 2 n n O n x n = 0 1 1 1 t 2 1 4 x t 2 1 4 x 1 4 x 1 4 x t 2 d t .
Next, as we have
1 1 t 2 1 4 x t 2 1 4 x 1 4 x 1 4 x t 2 d t = 1 1 4 x arctanh t arctanh 1 4 x t 1 4 x t 2 + c ,
we obtain the identity
n = 1 2 n n O n x n = 1 1 4 x lim t 1 arctanh t arctanh 1 4 x t 1 4 x t 2 ,
where the limit is of the form . Now, using the definition of the inverse hyperbolic tangent function, we get
lim t 1 arctanh t arctanh 1 4 x t 1 4 x t 2 = 1 2 lim t 1 ln 1 1 4 x t 1 4 x t 2 ln ( 1 t ) = 1 2 lim t 1 ln 1 4 x t 2 1 4 x t 1 t ln 1 4 x t 2 = 1 2 ln 1 1 4 x ln 1 4 x = ln 1 4 x ,
where we have used the continuity of the logarithm and
lim t 1 1 4 x t 2 1 4 x t 1 t = 1 1 4 x .
This completes the proof. □
Lemma 5. 
For all x [ 1 / 4 , 1 / 4 ) ,
F ( x ) = lim t 1 ln ( 1 4 x ) ln ( 1 t ) ln ( 1 4 x t ) ln 1 4 x 1 4 x t = ln ( 1 4 x ) ln 1 4 x 2 x , i f x < 0 ; 0 i f x = 0 ; ln ( 1 4 x ) ln 1 4 x 2 x i π , i f x > 0 ,
where i = 1 is the complex unit.
Proof. 
We find
lim t 1 ln ( 1 4 x ) ln ( 1 t ) ln ( 1 4 x t ) ln 1 4 x 1 4 x t = ln ( 1 4 x ) lim t 1 ln ( 1 t ) ln 1 4 x 1 4 x t = ln ( 1 4 x ) lim t 1 ln 1 t 1 4 x 1 4 x t ,
and by applying L’Hopital’s rule, we obtain the result F ( x ) = ln ( 1 4 x ) ln 1 4 x 2 x .
Theorem 17. 
For all x [ 1 / 4 , 1 / 4 ) , we have
n = 1 2 n n H n O n x n = 1 1 4 x Li 2 1 1 4 x Li 2 1 1 4 x π 2 4 + G ( x )
with
G ( x ) = 1 2 ln ( 1 4 x ) ln 1 4 x x , i f x < 0 ; ln 1 4 x x i π , i f x > 0 .
For x = 0 the series on the left-hand side of (53) is obviously zero.
Proof. 
We apply the previous theorem with a n = 2 n n O n and f ( x ) = ln ( 1 4 x ) 2 1 4 x , and combine it with (12). This yields
n = 1 2 n n H n O n x n = I ( x ) 2 1 4 x ,
where
I ( x ) = 0 1 1 4 x ln ( 1 4 x t ) 1 4 x t ln ( 1 4 x ) ( 1 t ) 1 4 x t d t .
Now, one checks that
1 4 x ln ( 1 4 x t ) 1 4 x t ln ( 1 4 x ) ( 1 t ) 1 4 x t d t = 2 Li 2 1 4 x t 1 4 x 2 Li 2 1 4 x t 1 4 x + ln ( 1 4 x ) ln ( 1 t ) + 2 ln ( 1 4 x t ) arctanh 1 4 x t 1 4 x + c .
This produces
I ( x ) = 2 Li 2 1 1 4 x 2 Li 2 1 1 4 x + 2 Li 2 ( 1 ) 2 Li 2 ( 1 ) + lim t 1 ln ( 1 4 x ) ln ( 1 t ) + 2 ln ( 1 4 x t ) arctanh 1 4 x t 1 4 x .
To determine the limit we again use the definition of the arctanh function and simplify. The result is
lim t 1 ln ( 1 4 x ) ln ( 1 t ) + 2 ln ( 1 4 x t ) arctanh 1 4 x t 1 4 x = ln ( 1 4 x ) ln 2 1 4 x + F ( x ) .
The result follows from (52) and further simplification steps. □
Remark 8. 
For all x [ 1 , 0 ) (and analogously for x ( 0 , 1 ) ), identity (53) can also be written as
n = 1 2 n n H n O n 4 n x n = 1 1 x Li 2 1 1 x Li 2 1 1 x π 2 4 + 1 2 ln ( 1 x ) ln 4 ( 1 x ) x .
In addition, we mention that an equivalent form of (53) is
n = 1 2 n n H n O n x n = 1 2 1 4 x × Li 2 ( 1 4 x ) 4 Li 2 1 4 x + π 2 2 + ln ( 1 4 x ) ln 1 4 x | x | 1 sgn x 2 π i ,
so that
n = 1 2 n n H n O n 4 n x n = 1 2 1 x Li 2 ( 1 x ) 4 Li 2 1 x + π 2 2 + ln ( 1 x ) ln 4 ( x 1 ) x i π .
Remark 9. 
To prove Theorem 17, an equivalent approach would be to start with a n = 2 n n H n , use Boyadzhiev’s expression (8) and apply (13). So, what we have shown above is the evaluation
0 1 2 1 4 x ln 1 + 1 4 x 2 1 4 x 2 1 4 x t 2 ln 1 + 1 4 x t 2 2 1 4 x t 2 d t 1 t 2 = 1 1 4 x Li 2 1 1 4 x Li 2 1 1 4 x π 2 4 + G ( x ) ,
where G ( x ) is given in (54).
Corollary 5. 
We have
n = 1 ( 1 ) n 2 n n H n O n 4 n = 2 8 8 Li 2 1 2 7 π 2 6 + 7 ln 2 2 ,
n = 1 2 n n H n O n 8 n = 2 Li 2 1 2 Li 2 ( 2 ) + π 2 12 9 8 ln 2 2 ,
and
n = 1 ( 1 ) n 2 n n H n O n 8 n = 2 3 Li 2 2 3 Li 2 2 3 π 2 4 ln 2 2 + ln 2 ln 3 2 + ln 2 3 2 .
Proof. 
To obtain the first series, set x = 1 / 4 in (53) and (54), and simplify using the identity G ( 1 4 ) = 3 2 ln 2 2 . The second series is obtained by setting x = 1 / 8 in (53) and (54). Note that
G 1 8 = ln 2 2 + i π 2 ln 2 ,
and this yields
n = 1 2 n n H n O n 8 n = 2 Li 2 ( 2 ) Li 2 ( 2 ) π 2 4 ln 2 2 + i π 2 ln 2 .
To remove the imaginary term, we simplify further using the dilogarithm property (4). Finally, insert x = 1 / 8 in (53) and simplify making use of the evaluation
G 1 8 = ln 2 2 + 1 2 ln 2 ln 3 + 1 2 ln 2 3 .
Corollary 6. 
Let α = ( 1 + 5 ) / 2 be the golden ratio. Then
n = 1 2 n n H n O n ( 4 α ) n = π 2 12 α + 5 2 α ln 2 α 2 α ln 2 ln α .
Proof. 
Set x = α / 4 in (53) and (54), and use that 1 1 / α = 1 / α in conjunction with
Li 2 ( α ) = π 2 10 ln 2 α , Li 2 ( α ) = 7 π 2 30 + 1 2 ln 2 α i π ln α ,
where the second dilogarithm identity follows from (4). □
Theorem 18. 
For all x ( 1 / 4 , 1 / 4 ) , we have
n = 1 2 n n n H n O n x n = 2 x ( 1 4 x ) 3 Li 2 1 1 4 x Li 2 1 1 4 x π 2 4 + G ( x ) + ln 1 4 x + 1 1 4 x 1 H ( x ) ,
with G ( x ) given in (54) and H ( x ) given by
H ( x ) = 1 4 x ln ( 1 4 x ) + 4 x ln 1 4 x x , i f x < 0 ; 1 4 x ln ( 1 4 x ) + 4 x ln 1 4 x x + i π , i f x > 0 .
Alternatively, we can write
n = 1 2 n n n H n O n x n = x ( 1 4 x ) 3 Li 2 ( 1 4 x ) 4 Li 2 ( 1 4 x ) + π 2 2 + 2 ln 1 4 x + 1 1 4 x 1 + ln 1 4 x | x | ln ( 1 4 x ) 2 ln ( 1 4 x ) 2 x H ( x ) ,
where
H ( x ) = i π ln ( 1 4 x ) , i f x < 0 ; 2 π i , i f x > 0 .
For x = 0 the series on the left-hand side of (56) is obviously zero.
Proof. 
Differentiate identity (53) with respect to x using the product rule. □
Three particular evaluations are offered in the next corollary.
Corollary 7. 
We have
n = 1 ( 1 ) n 2 n n n H n O n 8 n = 6 18 Li 2 2 3 Li 2 2 3 π 2 4 + 2 ln ( 2 + 3 ) 4 ln 2 + ln 3 + 1 2 ln 2 ln 3 ln 2 2 + 1 2 ln 2 3 ,
n = 1 2 n n n H n O n 8 n = 2 2 Li 2 1 2 Li 2 ( 2 ) + π 2 12 9 8 ln 2 2 + 2 ln ( 1 + 2 ) ,
n = 1 2 n n n H n O n 5 n = 2 5 Li 2 1 5 Li 2 ( 5 ) + π 2 12 + 2 ln α + 5 4 ln 5 1 8 ln 2 5 ,
where α is the golden ratio.
We also present a second remarkable evaluation of the complementary series from Corollary 6.
Corollary 8. 
We have
n = 1 2 n n n H n O n 4 n α n = π 2 α 2 24 + 5 4 α 2 ln 2 α + α 2 ln α 2 2 + α 3 ln α α 2 ln 2 ln α .

6. Series Associated with Generating Function (53) and Dilogarithm Functions

In this section, we derive a second remarkable evaluation of the companion series from Corollary 6, along with related series identities rising from the generating function (53) and properties of the dilogarithm function.
Writing ( 1 x ) / 4 for x and substituting x 2 for x, identity (53) can be written for x ( 0 , 1 ) as
x n = 1 2 n n H n O n 2 2 n ( 1 x 2 ) n = 2 χ 2 1 x π 2 4 + 2 ln 2 ln x + 2 ln 2 x ln x ln ( 1 x 2 ) ,
where
χ 2 y = 1 2 Li 2 ( y ) 1 2 Li 2 ( y )
is Legendre’s chi function ([15], p. 19) with the known values:
χ 2 tan π 8 = χ 2 2 1 = π 2 16 1 4 ln 2 2 1 , χ 2 ( β ) = π 2 12 3 4 ln 2 α , χ 2 ( β 3 ) = π 2 24 3 4 ln 2 α .
Lemma 6. 
If 0 u / 2 Z , then
0 1 x u 1 ln x d x = 1 u 2 , 0 1 x u 1 ln 2 x d x = 2 u 3 ,
0 1 x u 1 ln x ln ( 1 x 2 ) d x = ζ ( 2 ) H u / 2 ( 2 ) 2 u + H u / 2 u 2 ,
0 1 x u 1 χ 2 1 x d x = 2 H u H u / 2 2 u 2 + π 2 8 u .
Proof. 
We prove (59) and (60). A change of variable in (18) yields
0 1 x u 1 ln ( 1 x 2 ) d x = H u / 2 u , 0 u C Z ,
which, upon differentiation with respect to u, gives (59).
A simple change of variable in the integration formula (43) gives
0 1 x u 1 Li 2 ( x 2 ) d x = π 2 6 u 2 H u / 2 u 2 .
Identities (43) and (61) together with the following functional relations for the dilogarithm:
Li 2 x + Li 2 x = 1 2 Li 2 ( x 2 ) , Li 2 1 x + Li 2 x = π 2 3 1 2 ln 2 x + i π ln x , x < 1 , Li 2 1 x + Li 2 x = π 2 6 1 2 ln 2 x , x > 0 ,
lead directly to the determination of (60). □
Theorem 19. 
We have
n = 1 C n H n O n 2 2 n = π 2 2 + 4 ln 2 ,
n = 1 C n H n O n 2 2 n ( n + 2 ) = 14 9 + π 2 6 + 4 9 ln 2 ,
n = 1 H n O n ( 2 n + 1 ) ( 2 n + 3 ) = 1 2 + π 2 24 1 2 ln 2 ,
and more generally, for 0 r C Z ,
n = 1 2 n n n + r + 1 / 2 n H n O n 2 2 n = r + 1 r 2 2 H r 2 H r / 2 r 2 H r / 2 ( 2 ) + π 2 12 r 2 ln 2 + 4 r .
Proof. 
Multiply through (58) by x r 1 and integrate from 0 to 1 using Lemma 6 and the identity
0 1 ( 1 x 2 ) n x r d x = Γ ( n + 1 ) Γ ( r + 1 ) / 2 2 Γ n + 1 + ( r + 1 ) / 2 = 1 ( r + 1 ) n + ( r + 1 ) / 2 n .
Identities (62), (63), and (64) are obtained by evaluating (65) at r = 1 , r = 3 , and r = 2 , respectively. □
Corollary 9. 
If 0 r C Z , 2 r Z , then
n = 1 n + r n 2 n + 2 r + 1 2 n H n O n = 2 r + 1 4 r 2 2 ( H 2 r H r ) r H r ( 2 ) + π 2 6 r 2 ln 2 + 2 r .
Proof. 
Write 2 r for r in (65) and use r + 1 / 2 s = 2 r + 1 2 s 2 s s r s 1 2 2 s .
Theorem 20. 
We have
n = 1 C n H n H n + 1 O n 2 2 n = 16 ln 2 + 5 π 2 3 + 14 ζ ( 3 ) , n = 1 C n H n H n + 2 O n 2 2 n ( n + 2 ) = 236 27 + 28 27 ln 2 + 19 π 2 54 + 14 3 ζ ( 3 ) ,
and more generally, for 0 r C Z ,
n = 1 2 n n 2 2 n H n H n + ( r + 1 ) / 2 O n n + ( r + 1 ) / 2 n = 4 ( r + 2 ) r 3 H r H r / 2 + 4 ( r + 1 ) r 2 H r ( 2 ) 2 r + 3 r 2 H r / 2 ( 2 ) r + 1 r H r / 2 ( 3 ) ( 2 r + 1 ) r 2 π 2 6 + ( r + 1 ) r ζ ( 3 ) 4 ( r + 2 ) r 3 ln 2 + 8 ( 2 r + 3 ) r 4 + r + 1 r 2 H ( r + 1 ) / 2 2 H r 2 H r / 2 r 2 H r / 2 ( 2 ) + π 2 12 r 2 ln 2 + 4 r .
Proof. 
Differentiate (65) with respect to r to obtain
n = 1 2 n n 2 2 n H n O n H n + r + 1 / 2 H r + 1 / 2 n + r + 1 / 2 n = 4 ( r + 2 ) r 3 H r H r / 2 + 4 ( r + 1 ) r 2 H r ( 2 ) 2 r + 3 r 2 H r / 2 ( 2 ) r + 1 r H r / 2 ( 3 ) ( 2 r + 1 ) r 2 π 2 6 + ( r + 1 ) r ζ 3 4 ( r + 2 ) r 3 ln 2 + 8 ( 2 r + 3 ) r 4 ,
from which (66) follows upon second use of (65). □

7. A Second New Generating Function

This section deals with the generating function involving the product C n H n O n .
Theorem 21. 
For all x [ 1 / 4 , 1 / 4 ) , we have
n = 1 C n H n O n x n = 1 4 x 2 Li 2 ( 1 4 x ) 2 Li 2 ( 1 4 x ) + π 2 2 1 4 x + ln ( 1 4 x ) 4 ln ( 1 + 1 4 x ) ln ( 1 4 x ) ln ( 1 1 4 x ) + 1 4 x ln ( 1 4 x ) 2 K ( x ) + 4 ln 2 + π 2 2 + 2 1 4 x Li 2 1 1 4 x Li 2 1 1 4 x ln 1 4 x + 1 1 4 x 1
with K ( x ) given by
K ( x ) = ln x 1 4 x , i f x < 0 ; ln x 1 4 x + i π , i f x > 0 .
For x = 0 the series on the left-hand side of (67) is obviously zero.
Proof. 
Combining (11) with (12) yields
n = 1 C n H n O n x n = 1 2 x 0 1 J ( x , t ) t ( 1 t ) d t
with
J ( x , t ) = t 1 + 1 4 x t 1 4 x t 1 4 x t 2 ln ( 1 4 x t ) + 1 4 x 2 ln ( 1 4 x ) t .
We have
I C ( x , t ) = J ( x , t ) t ( 1 t ) d t = 1 2 ( 2 Li 2 1 4 x t 2 Li 2 1 4 x t 2 1 4 x Li 2 1 4 x t 1 4 x + 2 1 4 x Li 2 1 4 x t 1 4 x 1 4 x ( ln ( 1 4 x ) 2 ) ln ( 1 t ) ln ( 1 4 x t ) ln ( 1 1 4 x t ) + ( ln ( 1 4 x t ) 4 ) ln ( 1 + 1 4 x t ) + 1 4 x ln ( 1 4 x t ) ln 1 1 4 x t 1 4 x 1 4 x ln ( 1 4 x t ) ln 1 + 1 4 x t 1 4 x + 4 1 4 x arctanh 1 4 x t 1 4 x ) + c .
From here, we calculate
2 I C ( x , 1 ) = 2 Li 2 1 4 x 2 Li 2 1 4 x + π 2 2 1 4 x ln ( 1 4 x ) ln ( 1 1 4 x ) + ln ( 1 4 x ) 4 ln 1 + 1 4 x ln 2 1 4 x ln ( 1 4 x ) + 2 ln 2 1 4 x + 1 4 x lim t 1 ln ( 1 4 x ) 2 ln ( 1 t ) + ln 1 1 4 x t 1 4 x ln ( 1 4 x t ) 2
with
lim t 1 ln ( 1 4 x ) 2 ln ( 1 t ) + ln 1 1 4 x t 1 4 x ln ( 1 4 x t ) 2 = ln ( 1 4 x ) 2 lim t 1 ln 1 4 x 1 4 x t 1 t 1 2 ln ( 1 4 x ) = ln ( 1 4 x ) 2 ln 2 x 1 4 x .
By a similar calculation, we get
2 I C ( x , 0 ) = π 2 2 4 ln 2 + 2 1 4 x ln 1 4 x + 1 1 4 x 1 2 1 4 x Li 2 1 1 4 x Li 2 1 1 4 x ,
where it was used that
lim t 0 ln ( 1 4 x t ) ln 1 1 4 x t = 0 .
Gathering terms and simplifying proves the theorem. □
In the following examples, to retain and acknowledge the symmetrical structures, it is even better to keep the purely imaginary terms in the evaluations (which can be canceled out in view of (4)).
Corollary 10. 
We have
n = 1 ( 1 ) n 4 n C n H n O n = 2 Li 2 ( 2 ) Li 2 ( 2 ) 2 2 Li 2 1 2 Li 2 1 2 + 4 2 2 ln 2 + 4 ln ( 2 + 1 ) π 2 2 ( 1 + 2 ) + ( 3 2 ln 2 6 2 4 ) ln 2 + i π ln 2
and
n = 1 ( 1 ) n C n H n O n 8 n = 4 Li 2 3 2 Li 2 3 2 2 6 Li 2 2 3 Li 2 2 3 π 2 3 2 + 1 + 4 6 ln 3 + 2 8 ln 2 2 ln 3 2 4 ln 1 + 3 2 + 2 ln 3 2 ln 3 2 1 + i π + 6 ln 3 2 2 2 ln 2 + ln 3 .
Remark 10. 
An alternative version of (67) is
n = 1 C n H n O n x n + 1 = π 2 8 + ln 2 + 1 4 ( 1 1 4 x Li 2 ( 1 4 x ) 4 Li 2 ( 1 4 x ) π 2 2 2 1 4 x ln ( 1 4 x ) 2 1 4 x ln 1 1 4 x 1 + 1 4 x 4 ln ( 1 + 1 4 x ) + 1 4 x ln ( 1 4 x ) 2 ln ( | x | ) ( ln ( 1 4 x ) 1 ) 2 ) + 1 sgn x 8 i π 1 4 x ln ( 1 4 x ) 2 .

8. Conclusions

In this work, we have explored a wide range of infinite series involving central binomial coefficients, Catalan numbers, and harmonic-type sequences such as H n and O n . Through analytic techniques, including integral representations and generating functions, we derived several new identities, many of which are not expressible in terms of elementary functions and instead involve special functions like the dilogarithm.
A central contribution of this paper is the evaluation of two nontrivial classes of series
n = 1 2 n n H n O n x n and n = 1 C n H n O n x n ,
for which we obtained closed-form expressions or representations involving polylogarithmic terms. These results complement and extend earlier findings by many mathematicians and provide a deeper understanding of the interplay between harmonic numbers and combinatorial coefficients.
While many of the series considered resist expression in closed form, particularly those involving harmonic number products, we have identified special cases–such as those involving the golden ratio–where elegant simplifications are possible. These findings illustrate the rich structure and complexity inherent in such series and open avenues for further investigation, particularly in the context of special values of polylogarithms and their connections to combinatorics.

References

  1. K. Adegoke and R. Frontczak, Series associated with a forgotten identity of Nörlund, Integers 25 (2025), Paper No. A18.
  2. K. Adegoke, R. Frontczak and T. Goy, Fibonacci–Catalan series, Integers 22 (2022), Paper No. A110.
  3. K. Adegoke, R. Frontczak and T. Goy, Combinatorial sums, series and integrals involving odd harmonic numbers, Afr. Mat. (2025), to appear.
  4. H. Alzer and G. V. Nagy, Some identities involving central binomial coefficients and Catalan numbers, Integers 20 (2020), Paper No. A59.
  5. K. N. Boyadzhiev, Series with central binomial coefficients, Catalan numbers, and harmonic numbers, J. Integer Seq. 15 (2012), Article 12.1.7.
  6. H. Chen, Interesting series associated with central binomial coefficients, Catalan numbers and harmonic numbers, J. Integer Seq. 19 (2016), Article 16.1.5.
  7. J. Choi and D. Cvijović, Values of the polygamma functions at rational arguments. J. Phys. A Math. Theor. 40 (2007), no. 50, 15019–15028; Erratum in J. Phys. A Math. Theor. 43 (2010), 239801.
  8. W. Chu, Alternating convolutions of Catalan numbers, Bull. Braz. Math. Soc. (N.S.) 53 (2021), no. 1, 95–105.
  9. M. W. Coffey, On some log-cosine integrals related to ζ(2), ζ(4) and ζ(6), J. Comput. Appl. Math. 159 (2003), 205–215.
  10. T. Goy and M. Shattuck, Determinant formulas of some Toeplitz-Hessenberg matrices with Catalan entries, Proc. Indian Acad. Sci. Math. Sci. 129 (2019), no. 4, Paper No. 46.
  11. I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products, Elsevier/Academic Press, Amsterdam, 2015.
  12. T. Koshy, Catalan Numbers with Applications, Oxford University Press, Oxford, 2009.
  13. K. Kölbig, The polygamma function ψ ( k ) ( x ) for x = 1 4 and x = 3 4 , J. Comput. Appl. Math. 75 (1996), 43–46.
  14. D. H. Lehmer, Interesting series involving the central binomial coefficient, Amer. Math. Monthly 92 (1985), no. 7, 449–457.
  15. L. Lewin, Polylogarithms and Associated Functions, North-Holland Publ., New York-Amsterdam, 1981.
  16. C. Li and W. Chu, Evaluation of infinite series by integrals, Mathematics 10 (2022), no. 14, Paper No. 2444.
  17. C. Li and W. Chu, Improper integrals involving powers of inverse trigonometric and hyperbolic functions, Mathematics 10 (2022), no. 16, Paper No. 2980.
  18. C. Li and W. Chu, Infinite triple series by integrals, Publ. Inst. Math. (Beograd) (N.S.) 112(126) (2022), 131–142.
  19. C. Li and W. Chu, Evaluating infinite series involving harmonic numbers by integration, Mathematics 12 (2024), no. 4, Paper No. 589.
  20. J. Mikić, Two new identities involving the Catalan numbers and sign-reversing involutions, J. Integer Seq. 19 (2019) Article 19.7.7.
  21. F. Qi and B.-N. Guo, Integral representations of the Catalan numbers and their applications, Mathematics 5 (2017), no. 3, Paper No 40.
  22. S. Roman, An Introduction to Catalan Numbers, Birkhäuser/Springer, Cham, 2015.
  23. A. Sofo, D. Cvijović, Extensions of Euler harmonic sums, Appl. Anal. Discrete Math. 6 (2012), 317–328.
  24. A. Sofo and A. S. Nimbran, Euler sums and integral connections, Mathematics 7 (2019), no. 9, Paper No. 833.
  25. H. M. Srivastava and J. Choi, Series Associated with the Zeta and Related Functions, Kluwer Academic Publishers, Dordrecht, 2001.
  26. R. P. Stanley, Catalan Numbers, Cambridge University Press, Cambridge, 2015.
  27. S. M. Stewart, Explicit evaluation of some quadratic Euler-type sums containing double-index harmonic numbers, Tatra Mt. Math. Publ. 77 (2020), no. 1, 73–98.
  28. C. Xu, Identities for the shifted harmonic numbers and binomial coefficients, Filomat 31 (2017), no. 19, 6071–6086.
  29. W. Zhang and L. Chen, On the Catalan numbers and some of their identities, Symmetry 11 (2019), no. 1, Paper No. 62.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated