1. Introduction
Leonhard Euler in 1755 applied the Euler–Maclaurin sum formula to find the famous standard Euler asymptotic expansion for
as
:
where
is the Bernoulli number defined by
, and
is the Euler–Mascheroni constant.
Ramanujan [
1] proposed the following asymptotic expansion for
:
where
is the
n-th triangular number. However, Ramanujan did not give any formulas for the general terms and also without any proof. Rewrite the above formula as the following notation:
In 2008, Villarino [
2] established an explicit expression for the coefficient sequence
:
where
are the Bernoulli polynomials defined by
. In 2015, Chen and Cheng [
3] reconsidered Ramanujan’s formula and gave the following recurrence relation for
:
In 2019, Chen [
4] improved the recurrence relation as
Another Ramanujan-type harmonic number expansion was given by Wang [
5] in 2018,
where
h is a parameter and
is a coefficient sequence
In this paper, we rewrite Ramanujan’s harmonic number expansion into a similar form of Euler’s asymptotic expansion:
where
is the
nth pronic number, twice the
nth triangular number. In fact, we prove that the number
is the negative of the median Bernoulli number. The median Bernoulli number is studied by the author [
6] in 2005. Then, we have for
,
Then, we could rewrite Wang’s expansion Equation (
7) as follows:
We give simpler asymptotic expansion representations for
using Equations (
10) and (
11), which effectively integrate the results of Villarino, Chen and Cheng, Chen, and Wang (see Equations (
4)–(
7)) and make their representations more meaningful. We discuss some properties of the numbers
and the polynomials
in
Section 3 and
Section 4, respectively. For example, the Hankel determinant of
for any
x can be evaluated as
Furthermore, Chen [
7] gave a new asymptotic expansion. For any nonzero real number
r, the
n-th harmonic number
may have an asymptotic expansion as
n approaches infinity:
where the parameters
satisfy the following recurrence relation
Inspired by this, we give a more general asymptotic expansion in
Section 5 using Equation (
12). Given
real numbers with
,
, we get
We know that the formula with
is Equation (
14) (see ([
7], Theorem 2.3)). Since
,
will remove the
term. This will improve the approximation. Thus, it can be seen that there are a lot of investigations for the
case, see [
4,
8,
9,
10].
If
, then the asymptotic expansion will become
The parameters
and
in Equations (
15) and (
16) are determined by some recurrence relations, which will be illustrated in Theorems 2 and 3, respectively. At the end of this paper, we will compare how close these asymptotic formulas are to
.
2. Median Bernoulli Numbers and
Set
, for
. And for
,
or equivalently,
The corresponding matrix is represented as follows.
This matrix is called the “
-matrix” in [
6], which is a special Euler–Seidel matrix. Let
be the lower diagonal sequence of the
-matrix. The number
is the negative of the median Bernoulli number
, which is the upper diagonal sequence of the
-matrix, i.e.,
(ref. [
6]). Therefore, by [
6], and Equations (8), (15) and (16), we have
Let the ordinary generating function of
as follows.
Let
be the formal Laplace transform of
. Then, the following relation was obtained ([
6], Theorem 4.2, Equation (29))
Since
, we have that for
([
6], Equation (32)),
Using Villarino’s explicit formula for
, Equation (
4), we have for
,
On the other hand, we substitute
as Equation (
20) in Wang’s formula for
, (see Equation (
8)), we have for
,
Therefore, we conclude our result in the following.
Theorem 1. For , we havewhere is the nth pronic number, 3. Some Properties of
Let the ordinary generating function of
as follows.
Using the relation between the ordinary generating functions of
,
, and
of the
-matrix, we have the following relation ([
6], Theorem 4.2, Equation (29))
Then, the following identity is obtained ([
6], Equation (27)).
In the above formula, the formula obtained by substituting
appears in the recurrence relation of
given by Chen and Cheng [
3] in 2015 (see Equation (
5)). Furthermore, if we substitute
into the above identity, we obtain Equation (
6) given by Chen [
4] in 2019.
There are a lot of properties of
obtained from [
6]. For example, let the denominators and the numerators of the rational number
be
, respectively. We have the following properties ([
6], Theorem 1.1):
The denominator is a square-free integer.
The set of the all odd prime divisors of is
.
The denominator is an odd integer, for .
The largest power of 2 that divides the numerator is .
The ordinary generating function
has the following continued fraction representation ([
6], Theorem 5.5)
where for
,
Using this representation, we have the Hankel determinant of
(ref. [
6], Theorem 5.5)
Since the finite product in Equation (
28) is the Hankel determinant of
(see [
6], Equation (41)), we have
By Equation (
19) and an integral representation of
([
11], Equation (28))
we have an integral representation of
, for
,
4. Some Properties of
We first list
for
.
Differentiating Equation (
11) with respect to
x we obtain
On the other hand, we use Equation (
30) to get an integral representation of
:
Let us consider the function
. We express
as its binomial expansion.
We interchange the order of summation and the inner sum becomes
:
Using the inversion binomial theorem to Equation (
19) we have
Applying the above identity and Equation (
29), we know that the Hankel determinant of
is the same as the Hankel determinant of
. However, we use ([
12], Proposition 1), indeed that for any value of
x,
5. New Asymptotic Expansions
To derive our new asymptotic expansions are inspired by ([
7], Theorem 2.3). We need the following lemma.
Lemma 1 ([
7], Lemma 1).
If is an asymptotic expansion for as x approaches infinity. Given any real number r, the parameters are defined by and for ,Then is an asymptotic expansion for . Our new asymptotic expansions are derived from Equation (
23). It is note that
. Therefore, we divide into two cases depending on whether
h is
or not.
Theorem 2. Let r and h be any given real numbers with and . The n-th harmonic number has the following asymptotic expansion as n approaches infinity:where the parameters given by the recurrence relation Proof. Rewrite Equation (
37) as the following representation:
In view of Equation (
23), we have
Comparing the above two expressions, we know that
We apply Lemma 1 and get the result we want. □
Using a similar approach, we can easily derive the following theorem for the situation .
Theorem 3. Given a real number r with . The n-th harmonic number has the asymptotic expansion as n approaches infinity:where the parameters are defined by the following relation Chen [
7] discussed many properties of the
case. Therefore, we mainly deal with the case of
here.
The first few parameters
are:
For
in Equation (
39), the resulting asymptotic expansion is as follows ([
10], Equation (3.24)):
as
.
For
in Equation (
39), we obtain a new asymptotic expansion:
as
.
For
in Equation (
39), we obtain a new asymptotic expansion:
as
.
From a computational point of view, the formulas Equations (
42) and (
43) are better than Equation (
41).
It follows from Equations (
41)–(
43) that for
,
From
Table 1, we observe that, among approximation formulas Equations (
44)–(
46), for
, the formula Equation (
46) would be the best one. There seems to be an optimal real number
r in Equation (
39), and when we substitute it in this formula, the resulting approximation should be optimal. We guess that this real number
r should be close to