In mathematics, the n-th harmonic number is the sum of the reciprocals of the first n natural numbers:
Starting from n = 1, the sequence of harmonic numbers begins:
Harmonic numbers are related to the harmonic mean in that the n-th harmonic number is also n times the reciprocal of the harmonic mean of the first n positive integers.
Harmonic numbers have been studied since antiquity and are important in various branches of number theory. They are sometimes loosely termed harmonic series, are closely related to the Riemann zeta function, and appear in the expressions of various special functions.
The harmonic numbers roughly approximate the natural logarithm function [1] : 143 and thus the associated harmonic series grows without limit, albeit slowly. In 1737, Leonhard Euler used the divergence of the harmonic series to provide a new proof of the infinity of prime numbers. His work was extended into the complex plane by Bernhard Riemann in 1859, leading directly to the celebrated Riemann hypothesis about the distribution of prime numbers.
When the value of a large quantity of items has a Zipf's law distribution, the total value of the n most-valuable items is proportional to the n-th harmonic number. This leads to a variety of surprising conclusions regarding the long tail and the theory of network value.
The Bertrand-Chebyshev theorem implies that, except for the case n = 1, the harmonic numbers are never integers. [2]
n | Harmonic number, Hn | |||
---|---|---|---|---|
expressed as a fraction | decimal | relative size | ||
1 | 1 | 1 | ||
2 | 3 | /2 | 1.5 | |
3 | 11 | /6 | ~1.83333 | |
4 | 25 | /12 | ~2.08333 | |
5 | 137 | /60 | ~2.28333 | |
6 | 49 | /20 | 2.45 | |
7 | 363 | /140 | ~2.59286 | |
8 | 761 | /280 | ~2.71786 | |
9 | 7 129 | /2 520 | ~2.82897 | |
10 | 7 381 | /2 520 | ~2.92897 | |
11 | 83 711 | /27 720 | ~3.01988 | |
12 | 86 021 | /27 720 | ~3.10321 | |
13 | 1 145 993 | /360 360 | ~3.18013 | |
14 | 1 171 733 | /360 360 | ~3.25156 | |
15 | 1 195 757 | /360 360 | ~3.31823 | |
16 | 2 436 559 | /720 720 | ~3.38073 | |
17 | 42 142 223 | /12 252 240 | ~3.43955 | |
18 | 14 274 301 | /4 084 080 | ~3.49511 | |
19 | 275 295 799 | /77 597 520 | ~3.54774 | |
20 | 55 835 135 | /15 519 504 | ~3.59774 | |
21 | 18 858 053 | /5 173 168 | ~3.64536 | |
22 | 19 093 197 | /5 173 168 | ~3.69081 | |
23 | 444 316 699 | /118 982 864 | ~3.73429 | |
24 | 1 347 822 955 | /356 948 592 | ~3.77596 | |
25 | 34 052 522 467 | /8 923 714 800 | ~3.81596 | |
26 | 34 395 742 267 | /8 923 714 800 | ~3.85442 | |
27 | 312 536 252 003 | /80 313 433 200 | ~3.89146 | |
28 | 315 404 588 903 | /80 313 433 200 | ~3.92717 | |
29 | 9 227 046 511 387 | /2 329 089 562 800 | ~3.96165 | |
30 | 9 304 682 830 147 | /2 329 089 562 800 | ~3.99499 | |
31 | 290 774 257 297 357 | /72 201 776 446 800 | ~4.02725 | |
32 | 586 061 125 622 639 | /144 403 552 893 600 | ~4.05850 | |
33 | 53 676 090 078 349 | /13 127 595 717 600 | ~4.08880 | |
34 | 54 062 195 834 749 | /13 127 595 717 600 | ~4.11821 | |
35 | 54 437 269 998 109 | /13 127 595 717 600 | ~4.14678 | |
36 | 54 801 925 434 709 | /13 127 595 717 600 | ~4.17456 | |
37 | 2 040 798 836 801 833 | /485 721 041 551 200 | ~4.20159 | |
38 | 2 053 580 969 474 233 | /485 721 041 551 200 | ~4.22790 | |
39 | 2 066 035 355 155 033 | /485 721 041 551 200 | ~4.25354 | |
40 | 2 078 178 381 193 813 | /485 721 041 551 200 | ~4.27854 |
By definition, the harmonic numbers satisfy the recurrence relation
The harmonic numbers are connected to the Stirling numbers of the first kind by the relation
The harmonic numbers satisfy the series identities
and
These two results are closely analogous to the corresponding integral results
and
There are several infinite summations involving harmonic numbers and powers of π: [3] [ better source needed ]
An integral representation given by Euler [4] is
The equality above is straightforward by the simple algebraic identity
Using the substitution x = 1 −u, another expression for Hn is
The nth harmonic number is about as large as the natural logarithm of n. The reason is that the sum is approximated by the integral
whose value is ln n.
The values of the sequence Hn − ln n decrease monotonically towards the limit
where γ ≈ 0.5772156649 is the Euler–Mascheroni constant. The corresponding asymptotic expansion is
where Bk are the Bernoulli numbers.
A generating function for the harmonic numbers is
where ln(z) is the natural logarithm. An exponential generating function is
where Ein(z) is the entire exponential integral. The exponential integral may also be expressed as
where Γ(0, z) is the incomplete gamma function.
The harmonic numbers have several interesting arithmetic properties. It is well-known that is an integer if and only if , a result often attributed to Taeisinger. [5] Indeed, using 2-adic valuation, it is not difficult to prove that for the numerator of is an odd number while the denominator of is an even number. More precisely,
with some odd integers and .
As a consequence of Wolstenholme's theorem, for any prime number the numerator of is divisible by . Furthermore, Eisenstein [6] proved that for all odd prime number it holds
where is a Fermat quotient, with the consequence that divides the numerator of if and only if is a Wieferich prime.
In 1991, Eswarathasan and Levine [7] defined as the set of all positive integers such that the numerator of is divisible by a prime number They proved that
for all prime numbers and they defined harmonic primes to be the primes such that has exactly 3 elements.
Eswarathasan and Levine also conjectured that is a finite set for all primes and that there are infinitely many harmonic primes. Boyd [8] verified that is finite for all prime numbers up to except 83, 127, and 397; and he gave a heuristic suggesting that the density of the harmonic primes in the set of all primes should be . Sanna [9] showed that has zero asymptotic density, while Bing-Ling Wu and Yong-Gao Chen [10] proved that the number of elements of not exceeding is at most , for all .
The harmonic numbers appear in several calculation formulas, such as the digamma function
This relation is also frequently used to define the extension of the harmonic numbers to non-integer n. The harmonic numbers are also frequently used to define γ using the limit introduced earlier:
although
converges more quickly.
In 2002, Jeffrey Lagarias proved [11] that the Riemann hypothesis is equivalent to the statement that
is true for every integer n ≥ 1 with strict inequality if n > 1; here σ(n) denotes the sum of the divisors of n.
The eigenvalues of the nonlocal problem
are given by , where by convention , and the corresponding eigenfunctions are given by the Legendre polynomials . [12]
The nth generalized harmonic number of order m is given by
(In some sources, this may also be denoted by or )
The special case m = 0 gives The special case m = 1 reduces to the usual harmonic number:
The limit of as n → ∞ is finite if m > 1, with the generalized harmonic number bounded by and converging to the Riemann zeta function
The smallest natural number k such that kn does not divide the denominator of generalized harmonic number H(k, n) nor the denominator of alternating generalized harmonic number H′(k, n) is, for n=1, 2, ... :
The related sum occurs in the study of Bernoulli numbers; the harmonic numbers also appear in the study of Stirling numbers.
Some integrals of generalized harmonic numbers are
and
where A is Apéry's constant ζ(3),
and
Every generalized harmonic number of order m can be written as a function of harmonic numbers of order using
for example:
A generating function for the generalized harmonic numbers is
where is the polylogarithm, and |z| < 1. The generating function given above for m = 1 is a special case of this formula.
A fractional argument for generalized harmonic numbers can be introduced as follows:
For every integer, and integer or not, we have from polygamma functions:
where is the Riemann zeta function. The relevant recurrence relation is
Some special values are
where G is Catalan's constant. In the special case that , we get
where is the Hurwitz zeta function. This relationship is used to calculate harmonic numbers numerically.
The multiplication theorem applies to harmonic numbers. Using polygamma functions, we obtain
or, more generally,
For generalized harmonic numbers, we have
where is the Riemann zeta function.
The next generalization was discussed by J. H. Conway and R. K. Guy in their 1995 book The Book of Numbers . [1] : 258 Let
Then the nth hyperharmonic number of order r (r>0) is defined recursively as
In particular, is the ordinary harmonic number .
The Roman Harmonic numbers [13] , named after Steven Roman, were introduced by Daniel Loeb and Gian-Carlo Rota in the context of a generalization of umbral calculus with logarithms [14] . There are many possible definitions, but one of them, for , is
and
Of course,
If , they satisfy
Closed form formulas are
where is Stirling numbers of the first kind generalized to negative first argument, and
which was found by Donald Knuth.
In fact, these numbers were defined in a more general manner using Roman numbers and Roman factorials, that include negative values for . This generalization was useful in their study to define Harmonic logarithms.
The formulae given above,
are an integral and a series representation for a function that interpolates the harmonic numbers and, via analytic continuation, extends the definition to the complex plane other than the negative integers x. The interpolating function is in fact closely related to the digamma function
where ψ(x) is the digamma function, and γ is the Euler–Mascheroni constant. The integration process may be repeated to obtain
The Taylor series for the harmonic numbers is
which comes from the Taylor series for the digamma function ( is the Riemann zeta function).
When seeking to approximate Hx for a complex number x, it is effective to first compute Hm for some large integer m. Use that as an approximation for the value of Hm+x. Then use the recursion relation Hn = Hn−1 + 1/n backwards m times, to unwind it to an approximation for Hx. Furthermore, this approximation is exact in the limit as m goes to infinity.
Specifically, for a fixed integer n, it is the case that
If n is not an integer then it is not possible to say whether this equation is true because we have not yet (in this section) defined harmonic numbers for non-integers. However, we do get a unique extension of the harmonic numbers to the non-integers by insisting that this equation continue to hold when the arbitrary integer n is replaced by an arbitrary complex number x,
Swapping the order of the two sides of this equation and then subtracting them from Hx gives
This infinite series converges for all complex numbers x except the negative integers, which fail because trying to use the recursion relation Hn = Hn−1 + 1/n backwards through the value n = 0 involves a division by zero. By this construction, the function that defines the harmonic number for complex values is the unique function that simultaneously satisfies (1) H0 = 0, (2) Hx = Hx−1 + 1/x for all complex numbers x except the non-positive integers, and (3) limm→+∞ (Hm+x − Hm) = 0 for all complex values x.
This last formula can be used to show that
where γ is the Euler–Mascheroni constant or, more generally, for every n we have:
There are the following special analytic values for fractional arguments between 0 and 1, given by the integral
More values may be generated from the recurrence relation
or from the reflection relation
For example:
Which are computed via Gauss's digamma theorem, which essentially states that for positive integers p and q with p < q
Some derivatives of fractional harmonic numbers are given by
And using Maclaurin series, we have for x < 1 that
For fractional arguments between 0 and 1 and for a > 1,
In mathematics, the gamma function is one commonly used extension of the factorial function to complex numbers. The gamma function is defined for all complex numbers except the non-positive integers. For every positive integer n,
The Riemann zeta function or Euler–Riemann zeta function, denoted by the Greek letter ζ (zeta), is a mathematical function of a complex variable defined as
Euler's constant is a mathematical constant, usually denoted by the lowercase Greek letter gamma, defined as the limiting difference between the harmonic series and the natural logarithm, denoted here by log:
In mathematics, the polygamma function of order m is a meromorphic function on the complex numbers defined as the (m + 1)th derivative of the logarithm of the gamma function:
In mathematics, the digamma function is defined as the logarithmic derivative of the gamma function:
In mathematics, the Hurwitz zeta function is one of the many zeta functions. It is formally defined for complex variables s with Re(s) > 1 and a ≠ 0, −1, −2, … by
In mathematics, in the area of analytic number theory, the Dirichlet eta function is defined by the following Dirichlet series, which converges for any complex number having real part > 0:
In mathematics, the polylogarithm (also known as Jonquière's function, for Alfred Jonquière) is a special function Lis(z) of order s and argument z. Only for special values of s does the polylogarithm reduce to an elementary function such as the natural logarithm or a rational function. In quantum statistics, the polylogarithm function appears as the closed form of integrals of the Fermi–Dirac distribution and the Bose–Einstein distribution, and is also known as the Fermi–Dirac integral or the Bose–Einstein integral. In quantum electrodynamics, polylogarithms of positive integer order arise in the calculation of processes represented by higher-order Feynman diagrams.
The Basel problem is a problem in mathematical analysis with relevance to number theory, concerning an infinite sum of inverse squares. It was first posed by Pietro Mengoli in 1650 and solved by Leonhard Euler in 1734, and read on 5 December 1735 in The Saint Petersburg Academy of Sciences. Since the problem had withstood the attacks of the leading mathematicians of the day, Euler's solution brought him immediate fame when he was twenty-eight. Euler generalised the problem considerably, and his ideas were taken up more than a century later by Bernhard Riemann in his seminal 1859 paper "On the Number of Primes Less Than a Given Magnitude", in which he defined his zeta function and proved its basic properties. The problem is named after Basel, hometown of Euler as well as of the Bernoulli family who unsuccessfully attacked the problem.
In mathematics, the Glaisher–Kinkelin constant or Glaisher's constant, typically denoted A, is a mathematical constant, related to the K-function and the Barnes G-function. The constant appears in a number of sums and integrals, especially those involving gamma functions and zeta functions. It is named after mathematicians James Whitbread Lee Glaisher and Hermann Kinkelin.
In mathematics, the explicit formulae for L-functions are relations between sums over the complex number zeroes of an L-function and sums over prime powers, introduced by Riemann (1859) for the Riemann zeta function. Such explicit formulae have been applied also to questions on bounding the discriminant of an algebraic number field, and the conductor of a number field.
In mathematics, the Dirichlet beta function is a special function, closely related to the Riemann zeta function. It is a particular Dirichlet L-function, the L-function for the alternating character of period four.
In mathematics, the Riemann zeta function is a function in complex analysis, which is also important in number theory. It is often denoted and is named after the mathematician Bernhard Riemann. When the argument is a real number greater than one, the zeta function satisfies the equation
In mathematics, the multiplication theorem is a certain type of identity obeyed by many special functions related to the gamma function. For the explicit case of the gamma function, the identity is a product of values; thus the name. The various relations all stem from the same underlying principle; that is, the relation for one special function can be derived from that for the others, and is simply a manifestation of the same identity in different guises.
In probability theory and directional statistics, a wrapped Cauchy distribution is a wrapped probability distribution that results from the "wrapping" of the Cauchy distribution around the unit circle. The Cauchy distribution is sometimes known as a Lorentzian distribution, and the wrapped Cauchy distribution may sometimes be referred to as a wrapped Lorentzian distribution.
In mathematics, Ramanujan's Master Theorem, named after Srinivasa Ramanujan, is a technique that provides an analytic expression for the Mellin transform of an analytic function.
In mathematical physics, the Wu–Sprung potential, named after Hua Wu and Donald Sprung, is a potential function in one dimension inside a Hamiltonian with the potential defined by solving a non-linear integral equation defined by the Bohr–Sommerfeld quantization conditions involving the spectral staircase, the energies and the potential .
Gregory coefficientsGn, also known as reciprocal logarithmic numbers, Bernoulli numbers of the second kind, or Cauchy numbers of the first kind, are the rational numbers that occur in the Maclaurin series expansion of the reciprocal logarithm
The Bernoulli polynomials of the second kindψn(x), also known as the Fontana-Bessel polynomials, are the polynomials defined by the following generating function:
This article incorporates material from Harmonic number on PlanetMath, which is licensed under the Creative Commons Attribution/Share-Alike License.