Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Open Mathematics

formerly Central European Journal of Mathematics

Editor-in-Chief: Gianazza, Ugo / Vespri, Vincenzo

1 Issue per year

IMPACT FACTOR 2016 (Open Mathematics): 0.682
IMPACT FACTOR 2016 (Central European Journal of Mathematics): 0.489

CiteScore 2016: 0.62

SCImago Journal Rank (SJR) 2016: 0.454
Source Normalized Impact per Paper (SNIP) 2016: 0.850

Mathematical Citation Quotient (MCQ) 2016: 0.23

Open Access
See all formats and pricing
More options …
Volume 14, Issue 1


Volume 13 (2015)

Pochhammer symbol with negative indices. A new rule for the method of brackets

Ivan Gonzalez / Lin Jiu / Victor H Moll
Published Online: 2016-09-30 | DOI: https://doi.org/10.1515/math-2016-0063


The method of brackets is a method of integration based upon a small number of heuristic rules. Some of these have been made rigorous. An example of an integral involving the Bessel function is used to motivate a new evaluation rule.

Key Words: Method of brackets; Pochhammer symbol; Bessel integral

MSC 2010: 33C99; 33C10

1 Introduction

The evaluation of definite integrals is connected, in a surprising manner, to many topics in Mathematics. The last author has described in [15] and [17] some of these connections. Many of these evaluations appear in the classical table [13] and proofs of these entries have appeared in a series of papers starting with [16] and the latest one is [1]. An interesting new method of integration, developed in [10] in the context of integrals coming from Feynman diagrams, is illustrated here in the evaluation of and entry from the classical table mentioned above. Each of these entries have an assigned number and the one discussed here is 6:671:7 in [13]: I:=0J0(ax)sin(bx)dx={0if0<b<a,1/b2a2,if0<a<b.(1)

Here J0 is the Bessel function. The reader is referred to [3] for an introduction to these classical functions.

The new integration procedure has been named the method of brackets. It is based upon a small number of heuristic rules. This method is close in spirit to the classical umbral method, developed in the Combinatorial world by G. C. Rota. The reader will find in [18] an elementary introduction to this topic. The application of this method to Special Functions may be found, for instance, in a series of papers by G. Dattoli and collaborators, such as [4], [5] and [6]. The connections between the umbral method and the method of brackets is the subject of current work and results will be reported elsewhere.

The rules for the method of brackets are described in the next section. Sections 3 and 4 present the evaluation of (1). It turns out that a naive limiting process leads to a conflict between these two evaluations. The issue comes from a limit computation involving the Pochhammer symbol (a)k. It is the goal of this paper to clarify this conflict. This is presented in Section 5 leading to the proposal of a new rule dealing with the extension of the Pochhammer symbol to negative integer indices.

2 The method of brackets

A method to evaluate integrals over the half line [0, ∞), based on a small number of rules has been developed in [1012]. This method of brackets is described next. The heuristic rules are currently being made rigorous in [2] and [14]. The reader will find in [79] a large collection of evaluations of definite integrals that illustrate the power and flexibility of this method.

For a ∈ ℂ, the symbol a0xa1dx(2) is the bracket associated to the (divergent) integral on the right. The symbol ϕn:=(1)nΓ(n+1)(3) is called the indicator associated to the index n. The notation ϕi1i2...ir, or simply ϕ12...r, denotes the product ϕi1 ϕi2 ... ϕir.

Rules for the production of bracket series

Rule P1. If the function f is given by the formal power series f(x)=n=0anxαn+β1,(4) then the improper integral of f over the positive real line is formally written as the bracket series 0f(x)dx=nanαn+β.(5)

Rule P2. For α ∈ ℂ, the multinomial power (a1 + a2 + . . . + ar)α is assigned the r-dimension bracket series n1n2nrϕn1n2n1a1n1a1nrα+n1++nrΓ(α).(6)

Rules for the evaluation of a bracket series

Rule E1 The one-dimensional bracket series is assigned the value nϕnf(n)an+b=1|a|f(n)Γ(n),(7) where n* is obtained from the vanishing of the bracket; that is, n* solves an + b = 0.

The next rule provides a value for multi-dimensional bracket series where the number of sums is equal to the number of brackets.

Rule E2 Assume the matrix A = (aij) is non-singular, then the assignment is n1nrϕn1nrf(n1,,nr)a11n1++a1rnr+c1ar1n1++arrnr+cr=1|det(A)|f(n1,nr)Γ(n1)Γ(nr) where {ni} is the (unique) solution of the linear system obtained from the vanishing of the brackets. There is no assignment if A is singular.

Rule E3 Each representation of an integral by a bracket series has associated an index of the representation via

index = number of sums – number of brackets. (8)

It is important to observe that the index is attached to a specific representation of the integral and not just to integral itself. The experience obtained by the authors using this method suggests that, among all representations of an integral as a bracket series, the one with minimal index should be chosen.

The value of a multi-dimensional bracket series of positive index is obtained by computing all the contributions of maximal rank by Rule E2. These contributions to the integral appear as series in the free parameters. Series converging in a common region are added and divergent series are discarded.

Note 2.1The rule E1is a version of the Ramanujan’s Master Theorem. A complete justification of this rule is provided in [2]. The justification of the remaining rules is the subject of current work.

Note 2.2A systematic procedure in the simplification of the series obtained by this procedure has been used throughout the literature: express factorials in terms of the gamma function and the transform quotients of gamma terms into Pochhammer symbol, defined by (a)k=Γ(a+k)Γ(a).(9)

Any presence of a Pochhammer with a negative index k is transformed by the rule (a)k=(1)k(1a)k.(10)

The example discussed in the next two section provides motivation for an additional evaluation rule for the method of brackets.

3 A first evaluation of entry 6:671:7 in Gradshteyn and Ryzhik

The evaluation of (1) uses the series J0(ax)=m=0ϕma2mΓ(m+1)22mx2m(11) and sin(bx)=n=0ϕnΓ(n+1)Γ(2n+2)b2n+1x2n+1.(12)

Therefore the integral in (1) is given by I=m,nϕm,na2mb2n+1Γ(n+1)22mΓ(m+1)Γ(2n+2)2m+2n+2.(13)

The duplication formula for the gamma function transforms this expression to I=π2m,nϕm,na2mb2n+122m+2nΓ(m+1)Γ(n+3/2)2m+2n+2.(14)

Eliminating the parameter n using Rule E1 gives n* = –m – 1 and produces I=π2m=0ϕm1Γ(m+12)(ab)2m=1bm=0(1)mm!(a2b2)m(12)m=1bm=0(12)m1m!(a2b2)m=1b2a2.

The condition |b| > |a| is imposed to guarantee the convergence of the series on the third line of the previous argument.

The series obtained by eliminating the parameter m by m* = –n – 1 vanishes because of the factor Γ(m + 1) in the denominator. The formula has been established.

4 An alternative evaluation

A second evaluation of (1) begins with 0J0(ax)sin(bx)dx=m,nϕm,na2mb2n+122mm!(2n+1)!2m+2n+2.(15)

The evaluation of the bracket series is described next.

Case 1. Choose n as the free parameter. Then m* = –n – 1 and the contribution to the integral is I1=12n=0ϕn(a2)2n2b2n+11Γ(n)n!(2n+1)!Γ(n+1).(16)

Each term in the sum vanishes because the gamma function has a pole at the negative integers.

Case 2. Choose m as a free parameter. Then 2m + 2n + 2 = 0 gives n* = –m – 1. The contribution to the integral is I2=12mϕm(a2)2mb2m11m!n!(2n+1)!|n=m1Γ(m+1).(17)

Now write n!(2n+1)!=1(n+1)n+1(18) and (17) becomes I2=12bm=0(1)mm!(a2b)2m1(m)m.(19)

Transforming the term (–m)–m by using (x)n=(1)n(1x)n.(20) gives (m)m=(1)m(1+m)m.(21)

Replacing in (19) produces I2=12bk=0(12)k1k!(a2b2)k=121b2a2.

The method of brackets produces half of the expected answer. Naturally, it is possible that the entry in [13] is erroneous (this happens once in a while). Some numerical computations and the evaluation in the previous section, should convince the reader that this is not the case. The source of the error is the use of (20) for the evaluation of the term (–m)–m. A discussion is presented in Section 5.

5 Extensions of the Pochhammer symbol

Rule E1 of the method of brackets requires the evaluation of f(n*). In many instances, this involves the evaluation of the Pochhammer symbol (x)m for m ∉ ℕ. In particular, the question of the value (m)mformN(22) is at the core of the missing factor of 2 in Section 4.

The first extension of (x)m to negative values of n comes from the identity (x)m=(1)m(1x)m.(23)

This is obtained from (x)m=Γ(xm)Γ(x)=Γ(xm)(x1)(x2)(xm)Γ(xm)(24) and then changing the signs of each of the factors. This is valid as long as x is not a negative integer. The limiting value of the right-hand side in (23) as x → –km, with k ∈ ℕ, is (km)m=(1)m(km)!((k+1)m)!.(25)

On the other hand, the limiting value of the left-hand side is limε0(k(m+ε))(m+ε)=limε0Γ((k+1)m(k+1)ε)Γ(km+kε)=limε0Γ((k+1)ε)((k+1)ε)(k+1)mΓ(kε)(kε)km=limε0Γ((k+1)ε)Γ(kε)(1)(k+1)m(1+(k+1)ε)(k+1)m(1+kε)km(1)km=(1)m(km)!((k+1)m)!kk+1.

Therefore the function (x)–m is discontinuous at x = –km, with Direct(km)mLimiting(km)m=k+1k.(26)

For k = 1, this ratio becomes 2. This explains the missing ½ in the calculation in Section 4. Therefore it is the discontinuity of (10) at negative integer values of the variables, what is responsible for the error in the evaluation of the integral (1).

This example suggest that the rules of the method of brackets should be supplemented with an additional one:

Rule E4 Let k ∈ ℕ be fixed. In the evaluation of series, the rule (km)m=kk+1(1)m(km)!((k+1)m)!(27) must be used to eliminate Pochhammer symbols with negative index and negative integer base.

A variety of other examples confirm that this heuristic rule leads to correct evaluations.


The last author acknowledges the partial support of NSF-DMS 1112656. The second author is a graduate student, partially supported by the same grant.


  • [1]

    Amdeberhan T., Dixit A., Guan X., Jiu L., Kuznetsov A., Moll V., Vignat C., The integrals in Gradshteyn and Ryzhik. Part 30: Trigonometric integrals. Scientia, 27:47-74, 2016 Google Scholar

  • [2]

    Amdeberhan T., Espinosa O., Gonzalez, Harrison M., Moll V., Straub A., Ramanujan Master Theorem. The Ramanujan Journal, 29:103-120, 2012 Google Scholar

  • [3]

    Andrews G.E., Askey R., Roy R., Special Functions, volume 71 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, New York, 1999 Google Scholar

  • [4]

    Dattoli G., Cesarano C., Saccetti D., Miscellaneous results on the generating functions of special functions. Integral Transforms and Special Functions, 12:315-322, 2001 Google Scholar

  • [5]

    Dattoli G., Ricci P., Cesarano C., The Lagrange polynomials, the associated generalizations, and the Umbral Calculus. Integral Transforms and Special Functions, 14:181-186, 2003 Google Scholar

  • [6]

    Dattoli G., Ricci P., Cesarano C., Beyond the monomiality: the Monumbrality Principle. Jour. Comp. Anal. Appl., 6:77-83, 2004 Google Scholar

  • [7]

    Gonzalez I., Kohl K., Moll V., Evaluation of entries in Gradshteyn and Ryzhik employing the method of brackets. Scientia, 25:65-84, 2014 Google Scholar

  • [8]

    Gonzalez I., Moll V., Definite integrals by the method of brackets. Part 1. Adv. Appl. Math., 45:50-73, 2010Google Scholar

  • [9]

    Gonzalez I., Moll V., Straub A., The method of brackets. Part 2: Examples and applications. In T. Amdeberhan, L. Medina, and Victor H. Moll (Eds.), Gems in Experimental Mathematics, volume 517 of Contemporary Mathematics, pages 157-172. American Mathematical Society, 2010 Google Scholar

  • [10]

    Gonzalez I., Schmidt I., Optimized negative dimensional integration method (NDIM) and multiloop Feynman diagram calculation. Nuclear Physics B, 769:124-173, 2007 Google Scholar

  • [11]

    Gonzalez I., Schmidt I., Modular application of an integration by fractional expansion (IBFE) method to multiloop Feynman diagrams. Phys. Rev. D, 78:086003, 2008Google Scholar

  • [12]

    Gonzalez I., Schmidt I., Modular application of an integration by fractional expansion (IBFE) method to multiloop Feynman diagrams II. Phys. Rev. D, 79:126014, 2009 Google Scholar

  • [13]

    Gradshteyn I.S., Ryzhik I.M., Table of Integrals, Series, and Products. Edited by D. Zwillinger and V. Moll. Academic Press, New York, 8th edition, 2015.POCHHAMMER SYMBOL WITH NEGATIVE INDICES 7 Google Scholar

  • [14]

    Jiu L., On the method of brackets: rules, examples, interpretations and modifications. Submitted for publication, 2016 Google Scholar

  • [15]

    Moll V., The evaluation of integrals: a personal story. Notices of the AMS, 49:311-317, 2002 Google Scholar

  • [16]

    Moll V., The integrals in Gradshteyn and Ryzhik. Part 1: A family of logarithmic integrals. Scientia, 14:1-6, 2007 Google Scholar

  • [17]

    Moll V., Seized opportunities. Notices of the AMS, pages 476-484, 2010 Google Scholar

  • [18]

    Roman S., The Umbral Calculus. Dover, New York, 1984 Google Scholar

About the article

Received: 2015-07-11

Accepted: 2016-09-08

Published Online: 2016-09-30

Published in Print: 2016-01-01

Citation Information: Open Mathematics, Volume 14, Issue 1, Pages 681–686, ISSN (Online) 2391-5455, DOI: https://doi.org/10.1515/math-2016-0063.

Export Citation

© 2016 Gonzalez et al., published by De Gruyter Open. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License. BY-NC-ND 3.0

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

Karen T. Kohl
Advances in Applied Mathematics, 2017, Volume 89, Page 102

Comments (0)

Please log in or register to comment.
Log in