Ma · Mathematics

What Are the Main Properties of the Digamma Function?

Published By the 9to5Science editorial team

At a glance

  1. 1. Define the function

    The digamma function is the logarithmic derivative of the gamma function: ψ(z) = Γ′(z)/Γ(z).

  2. 2. Use the recurrence

    Differentiate Γ(z + 1) = zΓ(z) to obtain ψ(z + 1) = ψ(z) + 1/z.

  3. 3. Apply reflection and multiplication

    These identities relate values at reflected arguments and sums of values at evenly spaced arguments.

  4. 4. Check the domain

    The digamma function has poles at zero and the negative integers; avoid these arguments in its formulas.

The digamma function, written ψ(z)\psi(z), is the logarithmic derivative of the gamma function. Its identities are useful for simplifying expressions, evaluating special values, and studying the gamma function in complex analysis. Here are several standard properties, along with a worked evaluation using two of them.

Definition and domain

Wherever Γ(z)\Gamma(z) is finite and nonzero, define

ψ(z)=Γ′(z)Γ(z).\psi(z)=\frac{\Gamma'(z)}{\Gamma(z)}.

The gamma function extends meromorphically to the complex plane and has simple poles at 0,−1,−2,…0,-1,-2,\ldots. The digamma function has simple poles at those same points, each with residue −1-1. It is finite elsewhere. For positive real xx, ψ(x)\psi(x) is real-valued.

The logarithmic derivative is useful because differentiating a product becomes addition: if ff and gg are nonzero, then (fg)′/(fg)=f′/f+g′/g(fg)'/(fg)=f'/f+g'/g. This is why gamma-function identities give corresponding digamma identities after logarithmic differentiation.

Recurrence relation

The gamma recurrence is Γ(z+1)=zΓ(z)\Gamma(z+1)=z\Gamma(z). Taking its logarithmic derivative gives

ψ(z+1)=ψ(z)+1z,\psi(z+1)=\psi(z)+\frac{1}{z},

for zz away from the poles. Equivalently, ψ(z)=ψ(z+1)−1/z\psi(z)=\psi(z+1)-1/z. Repeated application shifts an argument by an integer; for a positive integer mm,

ψ(z+m)=ψ(z)+∑k=0m−11z+k,\psi(z+m)=\psi(z)+\sum_{k=0}^{m-1}\frac{1}{z+k},

provided none of the terms is singular. In particular, since ψ(1)=−γ\psi(1)=-\gamma, where γ\gamma is the Euler–Mascheroni constant,

ψ(n)=Hn−1−γ(n=1,2,3,…),\psi(n)=H_{n-1}-\gamma\qquad(n=1,2,3,\ldots),

with Hj=1+1/2+⋯+1/jH_j=1+1/2+\cdots+1/j and H0=0H_0=0.

Reflection formula

The gamma reflection identity is Γ(z)Γ(1−z)=π/sin⁡(πz)\Gamma(z)\Gamma(1-z)=\pi/\sin(\pi z). Logarithmically differentiating it yields

ψ(1−z)−ψ(z)=πcot⁡(πz),\psi(1-z)-\psi(z)=\pi\cot(\pi z),

for noninteger zz. This formula gives a difference of two digamma values; it does not say that the values are equal. For example, at z=1/3z=1/3, it gives

ψ ⁣(23)−ψ ⁣(13)=π3.\psi\!\left(\frac{2}{3}\right)-\psi\!\left(\frac{1}{3}\right)=\frac{\pi}{\sqrt{3}}.

Multiplication formula

For a positive integer mm, differentiating the gamma multiplication formula gives

∑k=0m−1ψ ⁣(z+km)=mψ(mz)−mlog⁡m.\sum_{k=0}^{m-1}\psi\!\left(z+\frac{k}{m}\right) =m\psi(mz)-m\log m.

This holds when the arguments avoid digamma poles. For m=2m=2, it reads

ψ(z)+ψ ⁣(z+12)=2ψ(2z)−2log⁡2.\psi(z)+\psi\!\left(z+\frac{1}{2}\right)=2\psi(2z)-2\log 2.

It is often called the digamma multiplication formula; the case m=2m=2 is also known as the duplication formula.

Series representation

For zz not equal to a nonpositive integer, a standard convergent series is

ψ(z)=−γ+∑n=0∞(1n+1−1n+z).\psi(z)=-\gamma+\sum_{n=0}^{\infty}\left(\frac{1}{n+1}-\frac{1}{n+z}\right).

The paired terms are essential: each difference is of order 1/n21/n^2 for large nn, so the series converges. This representation follows, for example, by logarithmically differentiating an infinite-product representation of Γ(z)\Gamma(z). At z=1z=1, every summand vanishes, recovering ψ(1)=−γ\psi(1)=-\gamma.

Worked example: evaluate ψ(1/3) and ψ(2/3)

The reflection formula at z=1/3z=1/3 gives the difference:

ψ ⁣(23)−ψ ⁣(13)=π3.\psi\!\left(\frac{2}{3}\right)-\psi\!\left(\frac{1}{3}\right)=\frac{\pi}{\sqrt{3}}.

Use the multiplication formula with m=3m=3 and z=1/3z=1/3 to find the sum. Since ψ(1)=−γ\psi(1)=-\gamma,

ψ ⁣(13)+ψ ⁣(23)+ψ(1)=3ψ(1)−3log⁡3,ψ ⁣(13)+ψ ⁣(23)=−2γ−3log⁡3.\begin{aligned} \psi\!\left(\frac{1}{3}\right)+\psi\!\left(\frac{2}{3}\right)+\psi(1) &=3\psi(1)-3\log 3,\\ \psi\!\left(\frac{1}{3}\right)+\psi\!\left(\frac{2}{3}\right) &=-2\gamma-3\log 3. \end{aligned}

Adding and subtracting these equations gives

ψ ⁣(13)=−γ−32log⁡3−π23,ψ ⁣(23)=−γ−32log⁡3+π23.\psi\!\left(\frac{1}{3}\right)=-\gamma-\frac{3}{2}\log 3-\frac{\pi}{2\sqrt{3}},\qquad \psi\!\left(\frac{2}{3}\right)=-\gamma-\frac{3}{2}\log 3+\frac{\pi}{2\sqrt{3}}.

The sum comes from multiplication, while the difference comes from reflection. Using both identities determines each value individually.

Common mistakes

  • Confusing ψ\psi with Γ\Gamma. The digamma function is Γ′/Γ\Gamma'/\Gamma, not the gamma function itself.
  • Reversing the recurrence sign. The forward shift is ψ(z+1)=ψ(z)+1/z\psi(z+1)=\psi(z)+1/z; shifting backward subtracts 1/z1/z.
  • Reversing the reflection difference. With the convention above, ψ(1−z)−ψ(z)=πcot⁡(πz)\psi(1-z)-\psi(z)=\pi\cot(\pi z).
  • Ignoring poles. Formulas involving ψ(z)\psi(z) cannot be evaluated at z=0,−1,−2,…z=0,-1,-2,\ldots; reflection also excludes integer zz.

These identities provide practical ways to shift arguments, relate reflected values, and evaluate sums of digamma values. For other arguments, the series representation or numerical evaluation can be used once the poles have been ruled out.

References

  1. NIST Digital Library of Mathematical Functions, Chapter 5: Gamma Function
  2. Abramowitz and Stegun, Handbook of Mathematical Functions, Chapter 6: Gamma and Related Functions

Spotted an error? Send a correction. Read how we write and check articles in our editorial policy.