Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am currently reading some material that makes heavy usage of Hypergeometric functions, and there is one particular point about applying Stirling's approximation to various terms consisting of Gamma-Functions that is not very clear to me.

We have the classical Stirling's approximation formula for the Gamma-Function in the form: $$ \Gamma(z)=\sqrt{2\pi}e^{-z}z^{z-1/2}\left(1+O\left(\frac{1}{|z|}\right)\right) $$

for $|\arg(z)|<\pi$ as $|z|\to\infty$. In absolute value: $$ |\Gamma(z)|=\sqrt{2\pi}e^{-\Re(z)}|z|^{\Re(z)-1/2}e^{-\Im(z)\arg(z)}\left(1+O\left(\frac{1}{|z|}\right)\right) $$

Now, there is also the shifted Stirling's approximation for the Gamma-Function due to C. Rowe, if I am not mistaken, that says: $$ \Gamma(z+a)=\sqrt{2\pi}e^{-z}z^{z+a-1/2}\left(1+O\left(\frac{1}{|z|}\right)\right) $$

uniformly for $|\arg(z)|\leq\pi-\varepsilon$, $a$ in a compact subset of $\mathbb{C}$ and some suitable fixed $\varepsilon>0$, as $|z|\to\infty$. In absolute value: $$ |\Gamma(z+a)|=\sqrt{2\pi}e^{-\Re(z)}|z|^{\Re(z+a)-1/2}e^{-\Im(z+a)\arg(z)}\left(1+O\left(\frac{1}{|z|}\right)\right) $$

My first question refers to terms of the type $$ \Gamma(az+b)\Gamma(cz+d) $$ for some complex numbers $a,b,c$ and $d$.

(Q1) Using the classical Stirling's approximation formula (i.e. not the shifted one), how can one obtain a meaningful aggregated asymptotics for the above expression?

I am asking this because there are several places that apply the non-shifted version of Stirling's formula to shifted Gamma factors without mentioning any details, which leaves the impression that this is a fairly standard or even trivial argument. Unfortunately, at this point I am unable to see its triviality. What bothers me here is the term $$ (az+b)^{az+b-1/2}(cz+d)^{cz+d-1/2} $$ as well as $\arg(az+b)$ and $\arg(cz+d)$, since the shifts by $c$ and $d$ break any easy manipulations.

I am naturally assuming that I am missing something very obvious here (as usual).

I have intentionally not specified anything about the parameters $a,b,c$ and $d$ because my question rather aims at the principle of applying the non-shifted Stirling's approximation to Gamma terms like the above one.

(Q2) Are there any other versions or forms of the Stirling's approximation for $\Gamma(z)$ that could be particularly useful for computing such kinds of asymptotics?

I will be extremely thankful if someone could give some insight in (principle of) the application of Stirling's approximatioin formula(s) to terms composed of Gamma factors!

share|cite|improve this question
Check out the book of Paris and Kaminsky – Sasha Sep 9 '11 at 15:29
The Digital Library of Mathematical Functions entry on asymptotics of the gamma function may also be useful. – Mike Spivey Sep 9 '11 at 17:39
@Mike Spivey: thank you for pointing me to this link. It seems to be a very useful reference for various things related to special functions and asymptotics in complex analysis. – ex-falso-quodlibet Sep 12 '11 at 21:59
@ex-falso-quodlibet: You're welcome. I use the DLMF all the time and am happy to introduce it to others! – Mike Spivey Sep 12 '11 at 22:00
up vote 5 down vote accepted

Let me take a stab at (Q1). I find it better working for $\log\Gamma(z)$ rather than $\Gamma(z)$, since $\Gamma(z) = \exp( \log\Gamma(z))$.

Stirling formula reads as follows: $$ \log\Gamma(z) \sim (z-\frac{1}{2}) \log(z) - z + \frac{1}{2} \log(2 \pi) + o(1) $$ for $\vert z \vert \to \infty$ and $ \vert \arg(z) \vert < \pi - \epsilon$.

Notice that shifted formula is a simple consequence of the above:

$$ \begin{eqnarray} \log\Gamma(z+a) &\sim& ( z+a -\frac{1}{2}) \log(z+a) - z - a + \frac{1}{2} \log(2 \pi) + o(1) \\ &\sim& ( z+a -\frac{1}{2}) \log(z) + ( z+a -\frac{1}{2}) \log(1+\frac{a}{z}) - z - a + \frac{1}{2} \log(2 \pi) + o(1) \\ &\sim& ( z+a -\frac{1}{2}) \log(z) + z \log(1+\frac{a}{z}) - z - a + \frac{1}{2} \log(2 \pi) + o(1) \\ &\sim& ( z+a -\frac{1}{2}) \log(z) - z + \frac{1}{2} \log(2 \pi) + o(1) \end{eqnarray} $$

Now $$ \begin{eqnarray} \log\Gamma(a z + b) + \log\Gamma(c z + d) & \sim & ( a z +b -\frac{1}{2}) \log(a z) - a z + \frac{1}{2} \log(2 \pi) + \\ & & ( c z +d -\frac{1}{2}) \log(c z) - c z + \frac{1}{2} \log(2 \pi) + o(1) \end{eqnarray} $$

Then it is a matter of recombining terms as $ ( \mathcal{A} + \mathcal{B} z) + (\kappa z + \rho -\frac{1}{2}) \log (\kappa z) - \kappa z + \frac{1}{2} \log(2 \pi) + o(1) $.

See the book of Paris and Kaminski to fill in the details.

share|cite|improve this answer
Thank you for this clarification, it is much appreciated! The $\log$-notation certainly has its advantages, especially when having to keep track of various coefficients. Also, you appear to prefer equivalence notation instead of $=$-notation, my guess is because you'd like to avoid abuse of notation and because this way you also don't have to write the bounded $arg$-term anyway. Your comment about recombining terms suggests about an useful application of Hadamard's Factorization Theorem (or more generally, Weierstrass Factorization Theorem for entire functions of order $>1$) :-) – ex-falso-quodlibet Sep 12 '11 at 21:49

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.