Harmonic Number
A harmonic number is a number of the form
arising from truncation of the harmonic series. A harmonic number can be expressed analytically as
where π gamma
is the Euler-Mascheroni constant and
π Psi(x)=psi_0(x)
is the digamma
function.
The first few harmonic numbers π H_n
are 1, π 3/2
, π 11/6
, π 25/12
, π 137/60
, ... (OEIS A001008
and A002805). The numbers of digits in the
numerator of π H_(10^n)
for π n=0
,
1, ... are 1, 4, 41, 434, 4346, 43451, 434111, 4342303, 43428680, ... (OEIS A114467),
with the corresponding number of digits in the denominator given by 1, 4, 40, 433,
4345, 43450, 434110, 4342302, 43428678, ... (OEIS A114468).
These digits converge to what appears to be the decimal digits of π log_(10)e=0.43429448...
(OEIS A002285).
The first few indices π n
such that the numerator of π H_n
is prime are given by 2, 3, 5, 8, 9, 21, 26, 41, 56, 62,
69, ... (OEIS A056903). The search for prime
numerators has been completed up to π 81780
by E. W. Weisstein (May 13, 2009), and the following
table summarizes the largest known values.
| π n | decimal digits | discoverer |
| 63942 | 27795 | E. W. Weisstein (Feb. 14, 2007) |
| 69294 | 30067 | E. W. Weisstein (Feb. 1, 2008) |
| 69927 | 30301 | E. W. Weisstein (Mar. 11, 2008) |
| 77449 | 33616 | E. W. Weisstein (Apr. 4, 2009) |
| 78128 | 33928 | E. W. Weisstein (Apr. 9, 2009) |
| 78993 | 34296 | E. W. Weisstein (Apr. 17, 2009) |
| 81658 | 35479 | E. W. Weisstein (May. 12, 2009) |
The denominators of π H_n
appear never to be prime except for the case π H_2=3/2
. Furthermore, the denominator is never a prime
power (except for this case) since the denominator is always divisible by the
largest power of 2 less than or equal to π n
, and also by any prime π p
with π n/2<p<=n
.
The harmonic numbers are implemented as [n].
The values of π n
such that π H_n
equals or exceeds 1, 2, 3, ... are given by 1, 4, 11, 31, 83, 227, 616, 1674, ...
(OEIS A004080). Another interesting sequence
is the number of terms in the simple continued
fraction of π H_(10^n)
for π n=0
,
1, 2, ..., given by 1, 8, 68, 834, 8356, 84548, 841817, 8425934, 84277586, ... (OEIS
A091590), which is conjectured to approach
π 12ln2/pi^2=0.8427659...
(OEIS A089729).
The definition of harmonic numbers can also be extended to the complex plane, as illustrated above.
Based on their definition, harmonic numbers satisfy the obvious recurrence equation
with π H_1=1
.
The number formed by taking alternate signs in the sum also has an explicit analytic form
| π H_n^' | π = | π sum_(k=1)^(n)((-1)^(k+1))/k |
(4)
|
| π Image | π = | π ln2+1/2(-1)^n[psi_0(1/2n+1/2)-psi_0(1/2n+1)] |
(5)
|
| π Image | π = | π ln2+1/2(-1)^n[H_((n-1)/2)-H_(n/2)]. |
(6)
|
π H_(2n)^'
has the particularly beautiful
form
Alternating harmonic numbers may be supported in a future version of the Wolfram Language.
The harmonic number π H_n
is never an integer except for π H_1
, which can be proved by using the strong triangle inequality
to show that the 2-adic value of π H_n
is greater than 1 for π n>1
. This result was proved in 1915 by Taeisinger, and the
more general results that any number of consecutive terms not necessarily starting
with 1 never sum to an integer was proved by KΕ±rschΓ‘k in 1918 (Hoffman
1998, p. 157).
The harmonic numbers have odd numerators and even denominators.
The π n
th
harmonic number is given asymptotically by
where π gamma
is the Euler-Mascheroni constant (Conway
and Guy 1996; Havil 2003, pp. 79 and 89), where the general π (2n)
th term is π zeta(1-2n)
, giving π -12
, 120, π -252
, 240, ... for π n=1
, 2, ... (OEIS A006953).
This formula is a special case of an Euler-Maclaurin
integration formulas (Havil 2003, p. 79).
Inequalities bounding π H_n
include
(Young 1991; Havil 2003, pp. 73-75) and
(DeTemple 1991; Havil 2003, pp. 76-78).
An interesting analytic sum is given by
(Coffman 1987). Borwein and Borwein (1995) show that
where π zeta(z)
is the Riemann zeta function. The first of
these had been previously derived by de Doelder (1991), and the third by Goldbach
in a 1742 letter to Euler (Borwein and Bailey 2003, pp. 99-100; Bailey et
al. 2007, p. 256). These identities are corollaries of the identity
(Borwein and Borwein 1995). Additional identities due to Euler are
| π sum_(n=1)^(infty)(H_n)/(n^2) | π = | π 2zeta(3) |
(23)
|
| π 2sum_(n=1)^(infty)(H_n)/(n^m) | π = | π (m+2)zeta(m+1)-sum_(n=1)^(m-2)zeta(m-n)zeta(n+1) |
(24)
|
for π m=2
,
3, ... (Borwein and Borwein 1995), where π zeta(3)
is ApΓ©ry's constant.
These sums are related to so-called Euler sums.
A general identity due to B. Cloitre (pers. comm., Jan. 7, 2006) is
where π (x)_n
is a Pochhammer symbol.
Gosper gave the interesting identity
| π sum_(i=0)^(infty)(z^iH_i)/(i!) | π = | π -e^zsum_(k=1)^(infty)((-z)^k)/(kk!) |
(26)
|
| π Image | π = | π e^z[lnz+Gamma(0,z)+gamma], |
(27)
|
where π Gamma(0,z)
is the incomplete gamma function and
π gamma
is the Euler-Mascheroni
constant.
G. Huvent (2002) found the beautiful formula
A beautiful double series is given by
(Bailey et al. 2007, pp. 273-274). Another double sum is
| π sum_(i=0)^(k-1)sum_(j=k)^n((-1)^(i+j-1))/(j-i)(n; i)(n; j)=sum_(i=1)^(k-1)(n; i)^2(H_(n-i)-H_i) |
(30)
|
for π 1<=k<=n
(Sondow 2003, 2005).
There is an unexpected connection between the harmonic numbers and the Riemann hypothesis.
Generalized harmonic numbers in power π r
can be defined by the relationship
where
| π H_(n,1)=H_n. |
(32)
|
These number are implemented as [n, r].
The numerators of the special case π H_(n,2)
are known as Wolstenholme
numbers. B. Cloitre (pers. comm., ) gave the surprising identity
which relates π H_(n,2)
to an indefinite version of a famous series for π zeta(2)
. π H_(n,2)
also satisfies
where π zeta(2)
is the Riemann zeta function. This follows
from the identity
where π gamma_1(z)
is the trigamma function since
For odd π r>=3
,
the generalized harmonic numbers have the explicit form
where π psi_r(n)
is the polygamma function, π Gamma(r)
is the gamma function,
and π zeta(r)
is the Riemann zeta function.
The 2-index harmonic numbers satisfy the identity
(P. Simon, pers. comm., Aug. 30, 2004).
Sums of the generalized harmonic numbers π H_(n,r)
include
for π |z|<1
,
where π Li_r(z)
is a polylogarithm,
where equations (40), (41), (42), and (44) are due to B. Cloitre (pers. comm., Oct. 4,
2004) and π Li_2(z)
is a dilogarithm. In general,
(P. Simone, pers. comm. June 2, 2003). The power harmonic numbers also obey the unexpected identity
(M. Trott, pers. comm.).
P. Simone (pers. comm., Aug. 30, 2004) showed that
where
This gives the special results
for π t=0,pi/2,pi
,
respectively.
Conway and Guy (1996) define the second-order harmonic number by
| π H_n^((2)) | π = | π sum_(i=1)^(n)H_i |
(54)
|
| π Image | π = | π (n+1)(H_(n+1)-1) |
(55)
|
| π Image | π = | π (n+1)(H_(n+1)-H_1), |
(56)
|
the third-order harmonic number by
and the π k
th-order
harmonic number by
A slightly different definition of a two-index harmonic number π c_n^((j))
is given by Roman (1992) in connection with the harmonic logarithm. Roman (1992) defines this
by
| π c_n^((0)) | π = | π {1 for n>=0; 0 for n<0 |
(59)
|
| π c_0^((j)) | π = | π {1 for j=0; 0 for j!=0 |
(60)
|
plus the recurrence relation
For general π n>0
and π j>0
,
this is equivalent to
and for π n>0
,
it simplifies to
For π n<0
,
the harmonic number can be written
where π |_n]!
is the Roman factorial and π s
is a Stirling
number of the first kind.
A separate type of number sometimes also called a "harmonic number" is a harmonic divisor number (or Ore number).
See also
ApΓ©ry's Constant, Book Stacking Problem, Egyptian Fraction, Euler Sum, Faulhaber's Formula, Harmonic Divisor Number, Harmonic Logarithm, Harmonic Series, Unit Fraction, Wolstenholme NumberRelated Wolfram sites
http://functions.wolfram.com/GammaBetaErf/HarmonicNumber/, http://functions.wolfram.com/GammaBetaErf/HarmonicNumber2/Portions of this entry contributed by Jonathan Sondow (author's link)
Explore with Wolfram|Alpha
More things to try:
References
Bailey, D. H.; Borwein, J. M.; Calkin, N. J.; Girgensohn, R.; Luke, D. R.; and Moll, V. H. Experimental Mathematics in Action. Wellesley, MA: A K Peters, 2007.Borwein, J. and Bailey, D. Mathematics by Experiment: Plausible Reasoning in the 21st Century. Wellesley, MA: A K Peters, 2003.Borwein, D. and Borwein, J. M. "On an Intriguing Integral and Some Series Related to π zeta(4)." Proc. Amer. Math. Soc. 123, 1191-1198, 1995.Coffman, S. W. "Problem 1240 and Solution: An Infinite Series with Harmonic Numbers." Math. Mag. 60, pp. 118-119, 1987.Conway, J. H. and Guy, R. K. The Book of Numbers. New York: Springer-Verlag, pp. 143 and 258-259, 1996.de Doelder, P. J. "On Some Series Containing π Psi(x)-Psi(y)
and π (Psi(x)-Psi(y))^2
for Certain Values of π x
and π y
." J. Comp. Appl. Math. 37, 125-141, 1991.DeTemple, D. W. "The Non-Integer Property of Sums of Reciprocals of Consecutive Integers." Math. Gaz. 75, 193-194, 1991.Flajolet, P. and Salvy, B. "Euler Sums and Contour Integral Representation." Experim. Math. 7, 15-35, 1998.Graham, R. L.; Knuth, D. E.; and Patashnik, O. "Harmonic Numbers" and "Harmonic Summation." Β§6.3 and 6.4 in Concrete Mathematics: A Foundation for Computer Science, 2nd ed. Reading, MA: Addison-Wesley, pp. 272-282, 1994.Gosper, R. W. "harmonic Summation and exponential gfs." math-fun@cs.arizona.edu posting, Aug. 2, 1996.Havil, J. Gamma: Exploring Euler's Constant. Princeton, NJ: Princeton University Press, 2003.Hoffman, P. The Man Who Loved Only Numbers: The Story of Paul ErdΕs and the Search for Mathematical Truth. New York: Hyperion, 1998.Huvent, G. "Autour de la primitive de π t^pcoth(alphat/2)
." Feb. 3, 2002. http://perso.orange.fr/gery.huvent/articlespdf/Autour_primitive.pdf.Roman, S. "The Logarithmic Binomial Formula." Amer. Math. Monthly 99, 641-648, 1992.Roman, S. The Umbral Calculus. New York: Academic Press, p. 99, 1984.Savio, D. Y.; Lamagna, E. A.; and Liu, S.-M. "Summation of Harmonic Numbers." In Computers and Mathematics (Ed. E. Kaltofen and S. M. Watt). New York: Springer-Verlag, pp. 12-20, 1989.Sloane, N. J. A. Sequences A001008/M2885, A002285/M3210, A002805/M1589, A004080, A006953/M2039, A056903, A082912, A089729, A091590, A096618, A114467, and A114468 in "The On-Line Encyclopedia of Integer Sequences."Sondow, J. "Criteria for Irrationality of Euler's Constant." Proc. Amer. Math. Soc. 131, 3335-3344, 2003.Sondow, J. "Problem 11026: An Identity Involving Harmonic Numbers." Amer. Math. Monthly 112, 367-369, 2005.Trott, M. "The Mathematica Guidebooks Additional Material: Harmonic Numbers Inversion." http://www.mathematicaguidebooks.org/additions.shtml#S_3_06.Young, R. M. "Euler's Constant." Math. Gaz. 75, 187-190, 1991.
Referenced on Wolfram|Alpha
Harmonic NumberCite this as:
Sondow, Jonathan and Weisstein, Eric W. "Harmonic Number." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/HarmonicNumber.html
