Ramanujan's Master Theorem

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.

The result is stated as follows:

If a complex-valued function $f(x)$ has an expansion of the form $f(x)=\sum_^\infty \frac(-x)^k$

then the Mellin transform of $f(x)$ is given by

$\int_0^\infty x^ f(x) \, dx = \Gamma(s)\,\varphi(-s)$

where $\Gamma(s)$ is the gamma function.

It was widely used by Ramanujan to calculate definite integrals and infinite series.

Higher-dimensional versions of this theorem also appear in quantum physics through Feynman diagrams.

A similar result was also obtained by Glaisher.

Alternative formalism

An alternative formulation of Ramanujan's master theorem is as follows:

$\int_0^\infty x^\left(\,\lambda(0) - x\,\lambda(1) + x^2\,\lambda(2) -\,\cdots\,\right) dx = \frac\,\lambda(-s)$

which gets converted to the above form after substituting $\lambda(n) \equiv \frac$ and using the functional equation for the gamma function.

The integral above is convergent for $0 < \operatorname(s) < 1$ subject to growth conditions on $\varphi$.

Proof

A proof subject to "natural" assumptions (though not the weakest necessary conditions) to Ramanujan's master theorem was provided by G. H. Hardy(chapter XI) employing the residue theorem and the well-known Mellin inversion theorem.

Recall Euler's representation of the Gamma function $\int_0^\infty x^e^dx=m^\Gamma(n)$

choosing $m=r^k$, multiplying both sides by $\frac$ and then summing over $0\le k$ to obtain:
$\sum_^\infty\frac\int_0^\infty x^e^dx=\Gamma(n)\sum_^\infty\frac$

Observing the sum in the RHS is a Taylor Series and writing $e^$ in its series expansion:

$\sum_^\infty\frac\int_0^\infty x^\sum_^\infty \fracdx=\Gamma(n)f(hr^+a)$

Rewriting the LHS:

$\sum_^\infty \frac\int_0^\infty x^ \sum_^\infty \fracdx=\Gamma(n)f(hr^+a)$

Then once again observing the sum over $k$ is a Taylor series:
$\int_0^\infty x^\sum_^\infty \fracdx=\Gamma(n)f(hr^+a)$

Finally defining $f(hr^+a)= \phi(N)$ and letting $F(x)=\sum_^\infty \frac$ we gain the master theorem:
$\int_0^\infty x^F(x)dx=\Gamma(n)\phi(-n)$

Application to Bernoulli polynomials

The generating function of the Bernoulli polynomials $B_k(x)$ is given by:

$\frac=\sum_^\infty B_k(x)\,\frac$

These polynomials are given in terms of the Hurwitz zeta function:

$\zeta(s,a) = \sum_^\infty \frac$

by $\zeta(1-n,a) = -\frac$ for $~ n \geq 1$.
Using the Ramanujan master theorem and the generating function of Bernoulli polynomials one has the following integral representation:

$\int_0^\infty x^\left(\frac-\frac\right) dx = \Gamma(s)\,\zeta(s,a) \!$

which is valid for $0 < \operatorname(s) < 1$.

Application to the gamma function

Weierstrass's definition of the gamma function

$\Gamma(x) = \frac\,\prod_^\infty \left(\,1 + \frac\,\right)^ e^ \!$

is equivalent to expression

$\log\Gamma(1+x) = -\gamma\,x + \sum_^\infty \frac\,(-x)^k$

where $\zeta(k)$ is the Riemann zeta function.

Then applying Ramanujan master theorem we have:

$\int_0^\infty x^ \frac \mathrm d x = \frac\frac \!$

valid for $0 < \operatorname(s) < 1$.

Special cases of $s = \frac$ and $s = \frac$ are

$\int_0^\infty \frac \, \mathrm d x = \frac\,\zeta\left( \frac \right)$

$\int_0^\infty \frac \, \mathrm d x = \sqrt \frac \zeta\left(\frac 5 4\right)$

Application to Bessel functions

The Bessel function of the first kind has the power series
$J_\nu(z)=\sum_^\infty \frac\bigg(\frac\bigg)^$

By Ramanujan's master theorem, together with some identities for the gamma function and rearranging, we can evaluate the integral

$\frac \int_0^\infty z^J_\nu(\sqrt)\,dz = \Gamma(s)\Gamma(s-\nu)$

valid for $0 < 2\operatorname(s) < \operatorname(\nu)+\tfrac$.

Equivalently, if the spherical Bessel function $j_\nu(z)$ is preferred, the formula becomes

$\frac \int_0^\infty z^j_\nu(\sqrt)\,dz = \Gamma(s)\Gamma\bigg(\frac+s-\nu\bigg)$

valid for $0 < 2\operatorname(s) < \operatorname(\nu)+2$.

The solution is remarkable in that it is able to interpolate across the major identities for the gamma function. In particular, the choice of $J_(\sqrt z)$ gives the square of the gamma function, $j_0(\sqrt)$ gives the duplication formula, $z^J_(\sqrt z)$ gives the reflection formula, and fixing to the evaluable $s=\frac$ or $s=1$ gives the gamma function by itself, up to reflection and scaling.

Bracket integration method

The bracket integration method (method of brackets) applies Ramanujan's master theorem to a broad range of integrals. The bracket integration method generates the integrand's series expansion, creates a bracket series, identifies the series coefficient and formula parameters and computes the integral.

Integration formulas

This section identifies the integration formulas for integrand's with and without consecutive integer exponents and for single and double integrals. The integration formula for double integrals may be generalized to any multiple integral. In all cases, there is a parameter value $n^$ or array of parameter values $N^$ that solves one or more linear equations derived from the exponent terms of the integrand's series expansion.

Consecutive integer exponents, 1 variable

This is the function series expansion, integral and integration formula for an integral whose integrand's series expansion contains consecutive integer exponents.
$\begin &f(y)=\sum_^\infty \frac \ \varphi(n) \ y^ \\ &\int_0^\infty y^ f(y) \, dy \\ &= \int_0^\infty \sum_^\infty \frac \ \varphi(n) \ y^ dy \\ &=\Gamma(-n^)\,\varphi(n^). \end$
The parameter $n^$ is a solution to this linear equation.
$n^+c=0, \ n^= -c$

General exponents, 1 variable

Applying the substitution $y=x^$ generates the function series expansion, integral and integration formula for an integral whose integrand's series expansion may not contain consecutive integer exponents.
$\begin &f(x)=\sum_^\infty \frac \ \varphi(n) \ x^ \\ &\int_0^\infty x^ f(x) \, dx \\ &=\int_0^\infty \sum_^\infty \frac \ \varphi(n) \ x^ dx \\ &=a^ \ \Gamma(-n^)\,\varphi(n^). \\ \end$
The parameter $n^$ is a solution to this linear equation.
$a \ n^+c=0, \ n^=-a^ c$

Consecutive integer exponents, double integral

This is the function series expansion, integral and integration formula for a double integral whose integrand's series expansion contains consecutive integer exponents.
$\begin &f(y_,y_)=\sum_^\infty \frac \frac \ \varphi(n_,n_) \ y_^ \ y_^ \\ &\int_0^\infty y^_ y^_ \ f(y_,y_) \ dy_ \ dy_ \\ &=\int_0^\infty \int_0^\infty \sum^\infty_ \sum^\infty_ \frac \frac \ \varphi(n_,n_) \ y_^ \ y_^ \ dy_ \ dy_ \\ &=\Gamma(-n^_) \ \Gamma(-n^_) \ \varphi(n^_,n^_). \\ \end$
The parameters $n^_$ and $n^_$ are solutions to these linear equations.
$n^_+c_=0, \ n^_+c_=0, \ n^_=-c_, \ n^_=-c_$

General exponents, double integral

This section describes the integration formula for a double integral whose integrand's series expansion may not contain consecutive integer exponents. Matrices contain the parameters needed to express the exponents in a series expansion of the integrand, and the determinant of invertible matrix $A$ is $\det, A,$.
$A= \begin a_ & a_ \\ a_ & a_ \end, \ C= \begin c_ \\ c_ \end, \ \ N^= \begin n^_ \\ n^_ \end$
Applying the substitution $y_=x^_ x^_, \quad y_=x^_ x^_$ generates the function series expansion, integral and integration formula for a double integral whose integrand's series expansion may not contain consecutive integer exponents. The integral and integration formula are
$\begin &\int_0^\infty \int_0^\infty \sum^\infty_ \sum^\infty_ \frac \frac \ \varphi(n_,n_) \ x_^ \ x_^ \ dx_ \ dx_ \\ &=\det, A, ^ \ \Gamma(-n^_) \ \Gamma(-n^_) \ \varphi(n^_,n^_). \end$ The parameter matrix $N^$ is a solution to this linear equation.
$AN^+C=0, \ N^=-A^ C$.

Positive complexity index

In some cases, there may be more sums then variables. For example, if the integrand is a product of 3 functions of a common single variable, and each function is converted to a series expansion sum, the integrand is now a product of 3 sums, each sum corresponding to a distinct series expansion.
*The number of brackets is the number of linear equations associated with an integral. This term reflects the common practice of bracketing each linear equation.
* The complexity index is the number of integrand sums minus the number of brackets (linear equations). Each series expansion of the integrand contributes one sum.
* The summation indices (variables) are the indices that index terms in a series expansion. In the example, there are 3 summation indices $n_, n_$ and $n_$ because the integrand is a product of 3 series expansions.
* The free summation indices (variables) are the summation indices that remain after completing all integrations. Integration reduces the number of sums in the integrand by replacing the series expansions (sums) with an integration formula. Therefore, there are fewer summation indices after integration. The number of chosen free summation indices equals the complexity index.

Integrals with a positive complexity index

The free summation indices $\bar_,\ldots, \bar_$ are elements of set $F$. The matrix of free summation indices is $\bar$ and the coefficients of the free summation indices is matrix $\bar$.
$\bar= \begin \bar_ & \ldots &\bar_ \\ \vdots & & \vdots \\ \bar_ & \ldots & \bar_ \end, \ \bar= \begin \bar_ \\ \vdots \\ \bar_ \end$
The remaining indices are set $B$ containing indices $n_,\ldots, n_$.
Matrices $A, C$ and $N^$ contain matrix elements that multiply or sum with the non-summation indices. The selected free summation indices must leave matrix $A$ non-singular.
$A= \begin a_ & \ldots &a_ \\ \vdots & & \vdots \\ a_ & \ldots & a_ \end, \ C= \begin c_ \\ \vdots \\ c_ \end, \ \ N^= \begin n^_ \\ \vdots \\ n^_ \end$. This is the function's series expansion, integral and integration formula.
$\begin &f(x_, \ldots, x_)\\ &=\sum^\infty_ \sum^\infty_ \frac \frac \ldots \frac \frac \varphi(n_, \ldots, n_, \bar_, \dots, \bar_)\prod_ x_^ \\ &\int^_ \ldots \int^_ x^ \dots x^ f(x_, \ldots, x_) \ dx_ \ldots dx_ \\ &=\det, A, ^ \sum^\infty_ \frac \ldots \frac \ \Gamma(-n^_) \ldots \Gamma(-n^_) \ \varphi(n^_, \ldots, n^_, \bar_, \dots, \bar_). \end$
The parameters $n^_, \ldots, n^_$ are linear functions of the parameters $\bar^_, \ldots, \bar^_$.
$A \ N^+\bar \ \bar+C=0, \ N^= -A^(\bar \ \bar+C)$.

Bracket series

Bracket series notations are notations that substitute for common power series notations (Table 1). Replacing power series notations with bracket series notations transforms the power series to a bracket series. A bracket series facilitates identifying the formula parameters needed for integration. It is also recommended to replace a sum raised to a power: $\frac$
with this bracket series expression:$\sum^_ \ldots \sum^_ \ \phi_ \ x^_ \dots x^_\frac.$

Algorithm

This algorithm describes how to apply the integral formulas.

:Input Integral expression
:Output Integral value or integral cannot be assigned a value
:# Express the integrand as a power series.
:# Transform the integrand's power series to a bracket series.
:# Obtain the complexity index, formula parameters and series coefficient function.
:## Complexity index is the number of integrand sums minus number of brackets.
:## Parameters $n^$ or array $N^$ are solutions to linear equations $an^+c=0$ (zero complexity index, single integral), $AN^+C=0$ (zero complexity index, single integral) or $AN^+\bar \bar+C=0$ (positive complexity index).
:## Identify parameter $a$ or (zero complexity index, single integral) or compute $\det , A,$ (all other cases) from the associated linear equations.
:## Identify the series coefficient function $\varphi()$ of the bracket series.
:# If the complexity index is negative, return integral cannot be assigned a value.
:# If the complexity index is zero, select the formula from table 2 for zero complexity index, single or multiple integral, compute the integral value with this formula, and return this integral value.
:# If the complexity index is positive, select the formula from table 2 for positive complexity index, and compute the integral value as a series expansion with this formula for all possible choices of the free summation indices. Select the lowest complexity index, convergent series expansion, adding series that converge in the same region.
:## If all series expansions are divergent series or null series (all series terms zero), then return integral cannot be assigned a value.
:## If the series expansion is non-null and non-divergent, return this series expansion as the integral value.

Examples

Zero complexity index

The bracket method will integrate this integral.
$\int_0^\infty x^ \ e^ \ dx$

1. Express the integrand as a power series.$\int_0^\infty \sum_^\infty 2^ \ \frac \ x^ \ dx$


2. Transform the power series to a bracket series.
   $\sum_^\infty 2^ \ \phi(n) \cdot \left \langle 3 \ n+ \frac \right \rangle$


3. Obtain the complexity index, formula parameters and series coefficient function.

: Complexity index is zero.
:$3 \ n^+ 5/2=0$
:$n^= -5/6, \ a=3$
: $\varphi(n)=2^$.

4. Use table 2 to compute the integral.

$\int_0^\infty x^ \cdot e^ \ dx$
$=a^ \ \Gamma(-n^) \ \varphi(n^)$
$= \frac$

Positive complexity index

The bracket method will integrate this integral.
$\int^_ \frac \ dx$
1. Express the integrand as a power series. Use the sum raised to a power formula.
$\int^_ \sum_ \ \frac \phi_ 1^ x^ \langle n_+n_+n_+1/2 \rangle \ dx$
2. Transform the power series to a bracket series.
$\int^_ \sum_ \frac \phi_ \langle 5 \ n_+ 3 \ n_+1 \rangle \langle n_+n_+n_+1/2 \rangle$
3. Obtain the complexity index, formula parameters and series coefficient function.
: Complexity index is 1 as 3 sums and 2 brackets.
: Select $n_$ as the free index, $\bar_$. The linear equations, solutions, determinant and series coefficient are
$5n^_+ 3 \bar_+1=0, \ n^_+n^_+\bar_+1/2=0$
$\begin 1 & 1 \\ 0 & 5 \end \begin n^_ \\ n^_ \end + \begin 1 \\ 3 \end \begin \bar_\end+ \begin 1/2 \\ 1 \end =0$
$AN^ + \bar \bar + C =0$
$\det , A, =5$
$n^_= -\frac \bar_-\frac, \ n^_= -\frac \bar_-\frac.$
$\varphi(n^_,n^_, \bar_)= \frac=\frac$
4. Use table 2 to compute the integral
$\begin &\int^_ \frac \ dx \\ &=\sum^_ \frac \det , A, ^\Gamma(-n^_) \Gamma(-n^_) \varphi(n^_,n^_,\bar_) \\ &=\sum^_\frac \frac \end$

Citations

References

*
*
*
*
*
*
*
*
*
*
*
*

External links

*{{MathWorld, id=RamanujansMasterTheorem, title=Ramanujan's Master Theorem

Srinivasa Ramanujan
Theorems in analytic number theory