Startseite Compatibility of the method of brackets with classical integration rules
Artikel Open Access

Compatibility of the method of brackets with classical integration rules

  • Zachary Bradshaw , Ivan Gonzalez , Lin Jiu , Victor Hugo Moll EMAIL logo und Christophe Vignat
Veröffentlicht/Copyright: 9. Mai 2023

Abstract

The method of brackets is a symbolic approach to the computation of integrals over R n based on a deep result by Ramanujan. Its usefulness to obtain new and difficult integrals has been demonstrated many times in the last few years. This article shows that this method is consistent with most classical rules for the computation of integrals, such as the fundamental theorem of calculus, the Laplace transform, the reduction formula for the integration of functions with spherical symmetry, the Cauchy-Schlömilch transformation, and explicit evaluations for multivariate integrals of product of Bessel functions as obtained by Exton and Srivastava. This work is part of a program dedicated to the derivation of solid theoretical grounds for the use of this attractive integration method.

MSC 2010: 33B99; 33F99

1 Introduction

The problem of evaluating the definite integral

(1.1) I ( f ; a , b ) = a b f ( x ) d x

consists of expressing the value of I in terms of a collection of special functions F , evaluated at a set of special values C . Naturally, the sets F and C depend on the real numbers { a , b } and the integrand f . The typical situation is that the integrand f belongs to a class of functions F 0 and one seeks a larger class of functions F F 0 whose values produce the integral.

Example 1.1

For n N , the evaluation of

(1.2) a b x n d x = b n + 1 a n + 1 n + 1

is elementary. The only special function required is the primitive of x x n . In this case, the class of special functions reduces to F = x x n + 1 n + 1 and the special constants are C = { a , b } . A simple extension of (1.2) shows that the class of polynomials F = f ( x ) = k = 0 m a k x k : a k R , m { 0 , 1 , } is closed under integration.

The aforementioned example is the most elementary case of the Fundamental theorem of Calculus.

Theorem 1.2

Assume g is a differentiable function on the interval [ a , b ] , with derivative g . Then,

(1.3) a b g ( x ) d x = g ( b ) g ( a ) .

This result transforms the problem of evaluating the integral I in (1.1) into the question of finding a primitive; that is, given a function f , find g such that g = f . Then, (1.3) gives the value of the integral I .

Example 1.3

In order to evaluate the integral,

(1.4) I ( a , b ) = a b e x d x ,

observe that g ( x ) = e x is the primitive of f ( x ) = e x and Theorem 1.2 gives

(1.5) I ( a , b ) = a b e x d x = g ( b ) g ( a ) = e b + e a .

A simple generalization of this example shows that every polynomial P in the variable f ( x ) = e t x , with t R a fixed parameter, has a primitive of the same form. Therefore, the integral of P ( e t x ) between a and b can be evaluated using Theorem 1.2.

One of the main difficulties in the evaluation of definite integrals is the existence of elementary functions, such as those encountered in elementary courses, that do not have elementary primitives. In general, this is hard to establish. The notion of elementary function is left, for the purposes of this article, undefined. The reader will find in Ritt [1] more information about this topic.

Example 1.4

The symbolic language Mathematica gives the elementary evaluation

(1.6) 0 1 cos x d x = sin 1

as well as the nonelementary one

(1.7) 0 1 sin x x d x = SinIntegral ( 1 ) ,

where SinIntegral is the Mathematica notation for the sine integral function defined by:

(1.8) Si ( z ) = 0 z sin x x d x .

Therefore, (1.7) simply comes from the introduction of a new function. The question of whether Si may be expressed in terms of more elementary functions is more difficult to decide.

2 The method of brackets

In the study of interaction of elementary particles, R. Feynman introduced his now famous diagrams, from which the physical properties of the reaction can be read. Some parametrizations of these diagrams led to families of quite complicated definite integrals.

A Feynman diagram is a graphical method of representing the interactions of elementary particles (see [2,3] for an introduction). For the present exposition, such a diagram is simply a graph G with E + 1 external lines and N internal lines and L loops. The internal lines are sometimes referred to as propagators. All but one of these external lines are assumed to be independent. From its connection to Physics, the internal and external lines represent the particles that transfer momentum among the vertices of the diagram. Each of these particles carries a mass m i 0 for i = 1 , , N . The vertices represent the interaction of these particles, and conservation of momentum at each vertex assigns the momentum corresponding to the internal lines. A Feynman diagram has an associated integral given by the parametrization of the diagram. The reader will find all the details in [4,5].

Figure 1 
               The massless sunset diagram.
Figure 1

The massless sunset diagram.

Example

The massless sunset diagram is presented in Figure 1. The diagram is parametrized by the integral (in momentum space)

(2.1) G = d D q i π D / 2 d D k i π D / 2 1 q 2 k 2 ( p q k ) 2 ,

where k and q are the D -dimensional vectors. For reasons of convergence, it is traditional to let D = 4 2 ε and then ε 0 + . The Schwinger parametrization is a systematic procedure to express the integral G in terms of some parameters (see [4] for details). It produces the integral

(2.2) G = ( 1 ) D 0 0 0 exp x y z ( x y + x z + y z ) P ( x y + x z + y z ) D / 2 d x d y d z ,

with P being the square of the norm of the 4-vector p .

Several methods have been developed in order to deal with integrals of this form. The method of brackets used here is an effective variation of the method of negative dimension, one of the most used by the physics community, and consists of a small number of heuristic rules, some of which are in the process of being rigorously established. A complete description of the rules and a comparison with the current methods used to evaluate Feynman diagrams appear in [6].

This method of brackets computes the integral

(2.3) 0 f ( x ) d x

from an expansion of the form

(2.4) f ( x ) = n = 0 C ( n ) x α n + β 1

and introduces the bracket of a R by the divergent integral

(2.5) a = 0 x a 1 d x .

Formal replacement of series (2.4) in (2.3) expresses the integral as a bracket series:

(2.6) 0 f ( x ) d x = n = 0 C ( n ) α n + β .

In order to evaluate the integral, the method proposes a rule to assign a value to a bracket series.

Rule 1. For a C , the bracket associated with a is the divergent integral

(2.7) a = 0 x a 1 d x .

A fundamental (open) question is to give a rigorous definition of this bracket, consistent with the operational rules in the following.

Rule 2. The expansion of an arbitrary function. The method of brackets requires the expansion of the integrand in the form:

(2.8) f ( x ) = n ϕ n C ( n ) x β n + α ,

where C ( n ) , α , β C , and ϕ n , called the indicator of n , is defined by:

(2.9) ϕ n = ( 1 ) n Γ ( n + 1 ) .

For functions of several variables, say two, use the expansion

(2.10) f ( x 1 , x 2 ) = n 1 n 2 ϕ 12 C ( n 1 , n 2 ) x 1 β 1 n 1 + α 1 x 2 β 2 n 2 + α 2 ,

with the notation ϕ 12 used for ϕ n 1 ϕ n 2 .

Rule 3. Sum expansion. The expression ( A 1 + + A r ) μ is assigned the bracket expansion

( A 1 + + A r ) μ = n 1 n r ϕ n 1 ϕ n r A 1 n 1 A r n r μ + n 1 + + n r Γ ( μ ) .

This has been established in [7], and it is a direct consequence of the multinomial theorem.

Up to now, the integral is converted into a bracket series. The evaluation of these series is described next.

Rule 4. Evaluation of bracket series. Let a , b R . A one-dimensional bracket series is assigned the value

(2.11) n ϕ n C ( n ) a n + b = 1 a C ( n ) Γ ( n ) ,

where n is obtained from the vanishing of the bracket, that is, n solves a n + b = 0 . The sum reduces to a single value.

Remark 2.1

Observe that this evaluation requires the evaluation of the coefficients C ( n ) outside its original domain, that is, for n N . The conditions on the function f given in Theorem 2.2, guarantee the existence of a unique extension of C . It turns out that Rule 4 is precisely a result of Ramanujan, known these days as Ramanujan’s master theorem. This refers to the formal identity

(2.12) 0 x s 1 λ ( 0 ) x 1 ! λ ( 1 ) + x 2 2 ! λ ( 2 ) d x = Γ ( s ) λ ( s )

stated by S. Ramanujan’s in his Quarterly Reports [8, p. 298]. It was widely used by him as a tool in computing definite integrals and infinite series. In fact, as Hardy puts it in [9], he “was particularly fond of them, and used them as one of his commonest tools.”

Theorem 2.2

(Ramanujan’s Master Theorem). Let φ ( z ) be an analytic (single-valued) function, defined on a half-plane H ( δ ) = { z C : Re z δ } for some 0 < δ < 1 . Suppose that, for some A < π , φ satisfies the growth condition φ ( v + i w ) < C e P v + A w for all z = v + i w H ( δ ) . Then, for all 0 < Re s < δ ,

(2.13) 0 x s 1 { φ ( 0 ) x φ ( 1 ) + x 2 φ ( 2 ) } d x = π sin s π φ ( s ) .

An elementary argument adapts this statement to prove (2.11). See [10] for details.

Note

In the multidimensional evaluation of a bracket series, in the special case when the number of sums and brackets is the same, write

J = n 1 n r ϕ n 1 ϕ n r C ( n 1 , , n r ) a 11 n 1 + + a 1 r n r + c 1 a r 1 n 1 + + a r r n r + c r .

The multiple sum is declared to be

(2.14) J = 1 det ( A ) Γ ( n 1 ) Γ ( n r ) C ( n 1 , , n r ) ,

where A = { a i j } and the values { n i } are the solutions of the linear system obtained from the vanishing of the brackets:

(2.15) a 11 n 1 + + a 1 r n r = c 1 a r 1 n 1 + + a r n r = c r .

An inductive proof follows directly from the one-dimensional case.

3 Consistency with the fundamental theorem of calculus

This section shows that in the case, the integrand in (1.1) is a perfect derivative; then, the method of brackets is consistent with the fundamental theorem of calculus.

Theorem 3.1

The value of a b f ( x ) d x computed by the method of brackets gives a consistent result with the analytic procedure, that is,

(3.1) a b f ( x ) d x = f ( b ) f ( a ) .

Proof

Assume f ( x ) = k = 0 ϕ k C ( k ) x k , then

(3.2) f ( x ) = k = 1 ϕ k C ( k ) k x k 1 = k = 0 ϕ k + 1 C ( k + 1 ) ( k + 1 ) x k

and the last expression can be written as:

(3.3) f ( x ) = k = 0 ϕ k C ( k + 1 ) x k .

The change of variables t = x a b x maps ( a , b ) to ( 0 , ) , and it gives

a b f ( x ) d x = ( b a ) 0 f b t + a t + 1 d t ( t + 1 ) 2 = ( b a ) 0 k ϕ k C ( k + 1 ) ( b t + a ) k ( t + 1 ) k 2 d t = ( a b ) 0 k ϕ k C ( k + 1 ) = 0 k k b t a k m , n ϕ m , n t m m + n + k + 2 Γ ( k + 2 ) d t = ( a b ) k , m , n ϕ k , m , n C ( k + 1 ) Γ ( k + 2 ) = 0 k k b a k m + n + k + 2 m + + 1 .

The system

(3.4) m + n + k + 2 = 0 , m + + 1 = 0

coming from the vanishing of the brackets gives ( m , n ) = ( 1 , k 1 ) . Taking k and as free variables, it gives

(3.5) a b f ( x ) d x = ( a b ) k ϕ k C ( k + 1 ) = 0 k k b a k Γ ( + 1 ) Γ ( k + 1 ) Γ ( k + 2 ) .

Simplifying the previous summand gives

a b f ( x ) d x = ( a b ) k ϕ k C ( k + 1 ) = 0 k b a k k + 1 = k ϕ k C ( k + 1 ) b k + 1 a k + 1 k + 1 = k ϕ k + 1 C ( k + 1 ) ( b k + 1 a k + 1 ) = ( f ( b ) f ( 0 ) ) ( f ( a ) f ( 0 ) ) = f ( b ) f ( a ) .

The proof is complete. The solutions corresponding to other combination of free variables yield the same result.□

4 The Laplace transform

A common variety of definite integrals appear as the Laplace transform of a function. This is defined by:

(4.1) ( f ) ( s ) = 0 e x s f ( x ) d x .

If the function f admits the power series

(4.2) f ( x ) = n = 0 a n x α n + β 1 ,

then

(4.3) ( f ) ( s ) = n = 0 a n 0 e x s x α n + β 1 d x = n = 0 a n Γ ( n α + β ) s n α + β .

The next statement shows that (4.3) can be obtained directly by the method of brackets.

Theorem 4.1

The evaluation of the Laplace transform (4.1) by the method of brackets also produces (4.3).

Proof

The method of brackets gives

( f ) ( s ) = 0 e s x f ( x ) d x = 0 m , n ϕ m x m s m a n x α n β 1 d x = m , n ϕ m a n m + α n + β s m .

The vanishing of the bracket gives m = α n β , and therefore,

(4.4) ( f ) ( s ) = n a n s m Γ ( m ) = n = 0 a n Γ ( n α + β ) s n α + β ,

confirming (4.3).□

5 Radially symmetric multidimensional integrals

This section confirms that an integral of the form

(5.1) I = R m f ( x 1 2 + + x m 2 ) d x 1 d x 2 d x m

is evaluated correctly by the method of brackets.

The classical evaluation is performed using the n -dimensional spherical coordinates: r 2 = x 1 2 + + x m 2 and

x 1 = r cos φ 1 , x 2 = r sin φ 1 cos φ 2 , x m 2 = r sin φ 1 sin φ m 3 cos φ m 2 , x m 1 = r sin φ 1 sin φ m 2 cos φ m 1 , x m = r sin φ 1 sin φ m 2 sin φ m 1 ,

with 0 φ 1 π , 0 φ j 2 π for 2 j m and 0 r < and

(5.2) d x 1 d x m = r m 1 sin m 2 φ 1 sin φ m 2 d r d φ 1 d φ m 1 .

This gives

(5.3) I = 0 r m 1 f ( r 2 ) d r j = 1 m 2 0 π sin j ( φ m 1 j ) d φ m 1 j 0 2 π d φ m 1 .

The last integral is 2 π , and the integrals in the product are evaluated using the classical formula [11, 3.261.5]

(5.4) 0 π / 2 sin 2 a 1 φ cos 2 b 1 φ d φ = 1 2 B ( a , b ) = Γ ( a ) Γ ( b ) 2 Γ ( a + b )

to produce

(5.5) I = 2 π 0 r m 1 f ( r 2 ) d r j = 1 m 2 π Γ j 2 Γ ( j 2 + 1 )

and using Γ ( a + 1 ) = a Γ ( a ) , this reduces to

(5.6) I = 2 π m / 2 Γ ( m 2 ) 0 r m 1 f ( r 2 ) d r .

The change of variables t = r 2 gives, with the notation n = m / 2 ,

(5.7) I = π n Γ ( n ) 0 t n 1 f ( t ) d t .

The computation of I by the method of brackets begins with the expansion:

(5.8) f ( t ) = = 0 ϕ a ( ) t .

Then,

I = R m = 0 ϕ a ( ) ( x 1 2 + + x m 2 ) d x 1 d x m = 2 m R + m = 0 ϕ a ( ) ( x 1 2 + + x m 2 ) d x 1 d x m = 2 m R + m = 0 ϕ a ( ) n 1 , , n m ϕ n 1 , , n m x 1 2 n 1 x m 2 n m n 1 + + n m Γ ( ) d x 1 d x m = 2 m , n 1 , , n m ϕ , n 1 , , n m a ( ) Γ ( ) n 1 + + n m j = 1 m 2 n j + 1 .

The associated linear system coming from the vanishing of the brackets is

(5.9) A n 1 n m = 0 1 1 , where A = 1 1 1 1 0 2 0 0 0 0 0 2 .

The solution is ( , n 1 , , n m ) = m 2 , 1 2 , , 1 2 and det A = 2 m . Therefore,

(5.10) I = 2 m 1 2 m a m 2 Γ m 2 j = 1 m Γ 1 2 = π m / 2 a m 2 = π n a ( n ) .

Using the method of brackets in the opposite order as usual yields

a ( n ) = 1 Γ ( n ) k = 0 a k k + n = 1 Γ ( n ) 0 k = 0 ϕ k a k t k + n 1 d t = 1 Γ ( n ) 0 t n 1 f ( t ) d t .

Then, (5.10) yields

(5.11) I = π n Γ ( n ) 0 t n 1 f ( t ) d t ,

as in (5.7).

6 Cauchy-Schlömilch transformation

In this section, we show that the method of brackets is consistent with the Cauchy-Schlömilch transformation.

Theorem 6.1

The method of brackets is consistent with the Cauchy-Schlömilch transformation. That is, using the method of brackets to evaluate the integrals, we have

(6.1) f ( x 1 / x ) d x = f ( x ) d x .

Proof

Let f e ( x ) = ( f ( x ) + f ( x ) ) / 2 and f o ( x ) = ( f ( x ) f ( x ) ) / 2 be the even and odd parts of f , respectively. Then, f ( x ) = f e ( x ) + f o ( x ) , and we have

(6.2) f ( x 1 / x ) d x = f e ( x 1 / x ) d x + f o ( x 1 / x ) d x

(6.3) = 2 0 f e ( x 1 / x ) d x .

Similarly,

(6.4) f ( x ) d x = 2 0 f e ( x ) d x .

Thus, we have reduced (6.1) to the identity

(6.5) 0 f e ( x 1 / x ) d x = 0 f e ( x ) d x .

Put another way, it suffices to consider even functions, for the odd part vanishes. Without loss of generality, let us assume that f is even. Suppose now that f ( x ) = n = 0 ϕ n g ( n ) x α n + β 1 . Let ε > 0 . Then, we have

(6.6) f ( x ε / x ) = n = 0 ϕ n g ( n ) ( x ε / x ) α n + β 1

and

(6.7) f ( x + ε / x ) = n = 0 ϕ n g ( n ) ( x + ε / x ) α n + β 1 .

Define

(6.8) I ( ε ) 0 f ( x ε / x ) d x ,

so that we have

(6.9) d d ε I ( ε ) 0 f ( x ε / x ) x d x .

Now applying the method of brackets, we find

(6.10) d d ε I ( ε ) = 0 1 x n = 0 ϕ n g ( n ) ( α n + β 1 ) ( x ε / x ) α n + β 2 d x

(6.11) = 0 n = 0 k = 0 l = 0 ϕ n ϕ k ϕ l g ( n ) ( α n + β 1 ) ( 1 ) l ε l Γ ( 2 α n β ) x k l 1 2 α n β + k + l d x

(6.12) = n = 0 k = 0 l = 0 ϕ n ϕ k ϕ l g ( n ) ( α n + β 1 ) ( 1 ) l ε l Γ ( 2 α n β ) k l 2 α n β + k + l .

The vanishing of the brackets produces three series, one for each choice of free index. These are

(6.13) I l = 1 α l = 0 ϕ l g 2 l + 2 β α ( 2 l + 1 ) ( 1 ) l ε l Γ ( 2 l ) Γ ( l ) Γ β 2 l 2 α .

(6.14) I k = 1 α k = 0 ϕ k g 2 k + 2 β α ( 2 k + 1 ) ( 1 ) k ε k Γ ( 2 k ) Γ ( k ) Γ β 2 k 2 α .

(6.15) I n = 1 2 n = 0 ϕ n g ( n ) ( α n + β 1 ) ( 1 ) ( α n + β ) / 2 1 ε ( α n + β ) / 2 1 Γ ( 2 α n β ) Γ 1 α n + β 2 .

Since f is even, f is odd. Therefore, we have

(6.16) d d ε I ( ε ) 0 f ( ε / x x ) x d x .

Applying the method of brackets to this form of d d ε I ( ε ) , we have

(6.17) d d ε I ( ε ) = 0 1 x n = 0 ϕ n g ( n ) ( α n + β 1 ) ( ε / x x ) α n + β 2 d x

(6.18) = 0 n = 0 k = 0 l = 0 ϕ n ϕ k ϕ l g ( n ) ( α n + β 1 ) ( 1 ) l ε k Γ ( 2 α n β ) x l k 1 2 α n β + k + l

(6.19) = n = 0 k = 0 l = 0 ϕ n ϕ k ϕ l g ( n ) ( α n + β 1 ) ( 1 ) l ε k Γ ( 2 α n β ) l k 2 α n β + k + l .

The vanishing of the brackets produces three series, one for each choice of free index. These are

(6.20) I l + = 1 α l = 0 ϕ l g 2 l + 2 β α ( 2 l + 1 ) ( 1 ) l ε l Γ ( 2 l ) Γ ( l ) Γ β 2 l 2 α ,

(6.21) I k + = 1 α k = 0 ϕ k g 2 k + 2 β α ( 2 k + 1 ) ( 1 ) k ε k Γ ( 2 k ) Γ ( k ) Γ β 2 k 2 α ,

(6.22) I n + = 1 2 n = 0 ϕ n g ( n ) ( α n + β 1 ) ( 1 ) ( α n + β ) / 2 1 ε ( α n + β ) / 2 1 Γ ( 2 α n β ) Γ 1 α n + β 2 .

But these series exactly cancel the series arising from the first form of d d ε I ( ε ) . Thus, summing the solutions, we find that d d ε I ( ε ) = 0 . Then, I ( 1 ) = I ( 0 ) , and this completes the proof.□

Remark 6.2

An alternative form of the Schlömilch transformation on the half-line is usually written as:

(6.23) 0 f ( ( a x b x 1 ) 2 ) d x = 1 a 0 f ( x 2 ) d x .

Starting with the expansion

(6.24) f ( x ) = k 0 ϕ k F ( k ) x k ,

the usual procedure of the method of brackets evaluated the right-hand side as

(6.25) 1 a 0 f ( x 2 ) d x = 1 a k 0 ϕ k F ( k ) 2 k + 1 = Γ 1 2 2 a F 1 2 ,

and this is equivalent to

(6.26) 1 2 0 f ( x 2 ) d x = π 2 a F 1 2 .

In order to evaluate the left-hand side of (6.23), start with

(6.27) 0 f ( ( a x b x 1 ) 2 ) d x = k 0 ϕ k F ( k ) 0 ( a x b x 1 ) 2 k d x .

Rule 3 (sum expansion) of the method of brackets is now used to produce two expressions for the term ( a x b x 1 ) 2 k , namely,

(6.28) ( a x b x 1 ) 2 k = m = 0 Γ ( 2 k + m ) m ! Γ ( 2 k ) a m b 2 k m x 2 m 2 k

and

(6.29) ( a x b x 1 ) 2 k = n = 0 Γ ( 2 k + n ) n ! Γ ( 2 k ) a 2 k n b n x 2 k 2 n .

Replacing in (6.27) and using the standard rule for the method of brackets lead to three expressions for the integral I on the left-hand side.

Index k is free. This yields

(6.30) I ( 1 ) = i 2 b a k 0 ϕ k F ( k ) Γ ( 1 2 k ) Γ ( 1 2 k ) Γ ( 2 k ) ( a b ) k .

This is discarded since every term vanishes.

Index m is free. This yields

(6.31) I ( 2 ) = b 2 m = 0 F m + 1 2 Γ ( m + 1 ) Γ m 1 2 Γ ( 2 m 1 ) ( a b ) m m ! .

Index n is free. This yields

(6.32) I ( 3 ) = 1 2 a n = 0 F n 1 2 Γ ( n + 1 ) Γ n + 1 2 Γ ( 2 n + 1 ) ( a b ) n n ! .

The value of the desired integral is I ( 2 ) + I ( 3 ) . A simple manipulation of these series confirms the Schlömilch transformation (6.23).

7 Borwein integrals and an extension to Bessel functions

7.1 Borwein integrals

A Borwein integral is a definite integral of the form:

(7.1) 0 k = 1 N sin ( b k x ) b k x d x .

The integrals were introduced by David and Borwein in [12] and are somewhat famous for exhibiting consistency patterns that eventually break down. As an example, the following evaluations hold:

(7.2) 0 sin ( x ) x d x = π 2 ,

(7.3) 0 sin ( x ) x sin ( x / 3 ) x / 3 d x = π 2 ,

(7.4) 0 sin ( x ) x sin ( x / 3 ) x / 3 sin ( x / 5 ) x / 5 d x = π 2 ,

(7.5) 0 sin ( x ) x sin ( x / 3 ) x / 3 sin ( x / 5 ) x / 5 sin ( x / 13 ) x / 13 d x = π 2 .

However, the pattern fails at the next step. Indeed, we have

(7.6) 0 sin ( x ) x sin ( x / 3 ) x / 3 sin ( x / 5 ) x / 5 sin ( x / 13 ) x / 13 sin ( x / 15 ) x / 15 d x π 2 .

The correct evaluation of the integral is off from π / 2 by roughly 2.31 × 1 0 11 . It turns out that whenever b 1 = 1 and k = 2 N b k < 1 , the integral will evaluate to π / 2 . We will now show this with the method of brackets. We start from (7.1) and expand the sin ( b k x ) terms. Denoting the integral by I , we have

(7.7) I = 0 k = 1 N n k = 0 ϕ n k Γ ( n k + 1 ) Γ ( 2 n k + 2 ) b k 2 n k x 2 n k d x

(7.8) = n 1 , , n N = 0 ϕ n 1 , , n N Γ ( n 1 + 1 ) Γ ( n N + 1 ) Γ ( 2 n 1 + 2 ) Γ ( 2 n N + 2 ) b 1 2 n 1 b N 2 n N 2 n 1 + 2 n N + 1 .

Now, applying the Legendre duplication formula, this becomes

(7.9) I = n 1 , , n N = 0 ϕ n 1 , , n N 2 1 2 n 1 2 1 2 n N π N / 2 Γ ( n 1 + 3 / 2 ) Γ ( n N + 3 / 2 ) b 1 2 n 1 b N 2 n N 2 n 1 + + 2 n N + 1

(7.10) = π 2 N n 1 , , n N = 0 ϕ n 1 , , n N b 1 2 n 1 b N 2 n N 4 n 1 + + n N Γ ( n 1 + 3 / 2 ) Γ ( n N + 3 / 2 ) 2 n 1 + + 2 n N + 1 .

Now, the vanishing of the bracket yields N independent solutions, which we will label with a subscript j indicating that we are applying the method of brackets to the index n j = 1 / 2 i j n i . Then, the N series produced by the method of brackets are given by:

(7.11) I j = 1 2 π 2 N n i i j 2 i j ϕ n i b i 2 n i Γ ( n i + 3 / 2 ) b j 1 2 i j n i Γ 1 i j n i Γ 1 / 2 + i j n i

(7.12) = π 2 N 1 b j n i i j Γ 1 / 2 + i j n i ( i j Γ ( n i + 3 / 2 ) ) Γ 1 i j n i i j ( b i 2 / b j 2 ) n i n i !

(7.13) = π 2 N 1 b j n i i j Γ ( 1 / 2 ) ( 1 / 2 ) i j n i Γ ( 3 / 2 ) N 1 ( 1 ) i j n i ( i j ( 3 / 2 ) n i ) i j ( b i 2 / b j 2 ) n i n i ! ,

where ( x ) n denotes the Pochhammer symbol. According to [13], any presence of a Pochhammer symbol with a negative index must be treated, in the application of the method of brackets, using the rule:

(7.14) ( k m ) m = k k + 1 ( 1 ) m ( k m ) ! ( ( k + 1 ) m ) ! .

Replacing in (7.13)

(7.15) I j = π 2 b j n i i j ( 1 / 2 ) i j n i ( 0 ) i j n i ( i j ( 3 / 2 ) n i ) i j ( b i 2 / b j 2 ) n i n i !

(7.16) = π 2 b j F C ( N 1 ) [ 1 / 2 , 0 ; 3 / 2 , , 3 / 2 ; ( b 1 / b j ) 2 , , ( b j 1 / b j ) 2 , ( b j + 1 / b j ) 2 , , ( b N / b j ) 2 ] ,

where F C ( N 1 ) denotes the N 1 variable Lauricella hypergeometric function of type C. But this converges when i j b i / b j < 1 . That is, whenever i j b i < b j . Moreover, due to the zero Pochhammer in the numerator, only the n i = 0 for all i j term will survive. Thus, the evaluation in this region is given by π / ( 2 b j ) . Returning to our example, this explains why the pattern terminates. Indeed, we have b 1 = 1 and b i < 1 for all i 1 . So the only I j that has a possibility of surviving is the first one. But I 1 only survives when i = 2 N b i < 1 . This inequality fails at the eighth integral, in which a factor of sin ( x / 15 ) / ( x / 15 ) appears. Note that the method of brackets does not produce an evaluation for the cases where this inequality fails.

7.2 Extension to Bessel functions

The result on Borwein integrals can be extended to the integral of a product of Bessel functions. Indeed, recall the Bessel function of the first kind is given by:

(7.17) J ν ( x ) = 2 ν n = 0 ϕ n x 2 n + ν 4 n Γ ( n + ν + 1 ) .

Let us now compute

(7.18) I = 0 k = 1 N J ν k ( b k x ) ( b k x ) ν k d x .

We have

(7.19) I = 2 i = 1 N ν i 0 n 1 , , n N = 0 ϕ n 1 , , n N b 1 2 n 1 b N 2 n N x 2 n 1 + + 2 n N 4 n 1 + + n N Γ ( n 1 + ν 1 + 1 ) Γ ( n N + ν N + 1 ) d x

(7.20) 2 i = 1 N ν i n 1 , , n N = 0 ϕ n 1 , , n N b 1 2 n 1 b N 2 n N 4 n 1 + + n N Γ ( n 1 + ν 1 + 1 ) Γ ( n N + ν N + 1 ) 2 n 1 + + 2 n N + 1 .

Now, there are N choices for the set of free indices. The vanishing of the bracket corresponding to each choice is given by n j = 1 / 2 i j n i . Let M = i j n i . Then, the series representations produced by the method of brackets for each choice are

(7.21) I j = 2 1 i = 1 N ν i n i i j 2 i j ϕ n i b i 2 n i Γ ( n i + ν i + 1 ) b 1 2 M Γ ( 1 / 2 + M ) Γ ( 1 / 2 + ν j M )

(7.22) = 1 b j 2 i = 1 N ν i n i i j ( 1 ) M Γ ( 1 / 2 ) ( 1 / 2 ) M i j ( b i 2 / b j 2 ) n i n i ! Γ ( 1 / 2 + ν j ) ( 1 / 2 + ν j ) M i j Γ ( ν i + 1 ) ( ν i + 1 ) n i

(7.23) = π b j 2 i = 1 N ν i Γ ( 1 / 2 + ν j ) i j Γ ( ν i + 1 ) n i i j ( 1 / 2 ) M ( 1 / 2 ν j ) M i j ( ν i + 1 ) n i i j ( b i 2 / b j 2 ) n i n i !

(7.24) = π b j 2 i = 1 N ν i Γ ( 1 / 2 + ν j ) i j Γ ( ν i + 1 ) F C ( N 1 ) 1 2 , 1 2 ν j ; ν j ; b j ,

where ν j is the vector of length N 1 with components ν i + 1 for i j , b j is the vector of length N 1 with components ( b i / b j ) 2 for i j , ( x ) n denotes the Pochhammer symbol, and F C ( N 1 ) is the N 1 dimensional Lauricella hypergeometric function of type C, which converges in the region i j b i < b j .

Acknowledgements

Some of the results presented here are part of the doctoral dissertation of the first author. The last author wishes to thank the hospitality of the Mathematics Department of Tulane University while this work was being conducted.

  1. Conflict of interest: The authors state no conflict of interest.

References

[1] J. F. Ritt, Integration in Finite Terms. Liouville’s Theory of Elementary Functions, Columbia University Press, New York, 1948. 10.7312/ritt91596Suche in Google Scholar

[2] M. Polyak, Feynman diagrams for pedestrians and mathematicians, In: M Lyubich, L. Takhtajan (Eds.), Graphs and Patterns in Mathematics and Theoretical Physics, Vol. 73, Proceedings of Symposia in Pure Mathematics, American Mathematical Society, USA, 2005. 10.1090/pspum/073/2131010Suche in Google Scholar

[3] M. Veltman, Diagrammatica. The Path to Feynman Diagrams, Vol. 4, Cambridge Lecture Notes in Physics, Cambridge University Press, Cambridge, 1994. 10.1017/CBO9780511564079Suche in Google Scholar

[4] M. E. Peskin and D. V. Schroder, An Introduction to Quantum Field Theory, CRC Press, Boca Raton, 2019. 10.1201/9780429503559Suche in Google Scholar

[5] L. A. Takhtadzhyan, Quantum Mechanics for Mathematicians, American Mathematical Society, Providence, 2008. Suche in Google Scholar

[6] K. VanDusen, A Comparison of Negative-Dimensional Integration Techniques, PhD thesis, Tulane University, Louisiana, 2021. Suche in Google Scholar

[7] I. Gonzalez and V. H. Moll, Definite integrals by the method of brackets. Part 1, Adv. Appl. Math. 45 (2010), 50–73. 10.1016/j.aam.2009.11.003Suche in Google Scholar

[8] B. C. Berndt, Ramanujan’s Notebooks, Part I, Springer-Verlag, New York, 1985. 10.1007/978-1-4612-1088-7Suche in Google Scholar

[9] G. H. Hardy, Ramanujan. Twelve Lectures on Subjects Suggested by His Life and Work, 3rd Edition, Chelsea Publishing Company, New York, 1978. Suche in Google Scholar

[10] T. Amdeberhan, O. Espinosa, I. Gonzalez, M. Harrison, V. H. Moll, and A. Straub, Ramanujan’s Master Theorem, The Ramanujan J. 29 (2012), 103–120. 10.1007/s11139-011-9333-ySuche in Google Scholar

[11] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, In: D. Zwillinger, V. H. Moll (Eds.), 8th Edition, Academic Press, New York, 2015. Suche in Google Scholar

[12] D. Borwein and J. M. Borwein, Some remarkable properties of Sinc and related integrals, The Ramanujan J. 5 (2001), 73–89. 10.1023/A:1011497229317Suche in Google Scholar

[13] I. Gonzalez, L. Jiu, and V. H. Moll, Pochhammer symbols with negative indices. A new rule for the method of brackets, Open Math. 14 (2016), 681–686. 10.1515/math-2016-0063Suche in Google Scholar

Received: 2022-11-03
Accepted: 2023-04-01
Published Online: 2023-05-09

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Artikel in diesem Heft

  1. Special Issue on Future Directions of Further Developments in Mathematics
  2. What will the mathematics of tomorrow look like?
  3. On H 2-solutions for a Camassa-Holm type equation
  4. Classical solutions to Cauchy problems for parabolic–elliptic systems of Keller-Segel type
  5. Control of multi-agent systems: Results, open problems, and applications
  6. Logical perspectives on the foundations of probability
  7. Subharmonic solutions for a class of predator-prey models with degenerate weights in periodic environments
  8. A non-smooth Brezis-Oswald uniqueness result
  9. Luenberger compensator theory for heat-Kelvin-Voigt-damped-structure interaction models with interface/boundary feedback controls
  10. Special Issue on Fractional Problems with Variable-Order or Variable Exponents (Part II)
  11. Positive solution for a nonlocal problem with strong singular nonlinearity
  12. Analysis of solutions for the fractional differential equation with Hadamard-type
  13. Hilfer proportional nonlocal fractional integro-multipoint boundary value problems
  14. A comprehensive review on fractional-order optimal control problem and its solution
  15. The θ-derivative as unifying framework of a class of derivatives
  16. Review Articles
  17. On the use of L-functionals in regression models
  18. Minimal-time problems for linear control systems on homogeneous spaces of low-dimensional solvable nonnilpotent Lie groups
  19. Regular Articles
  20. Existence and multiplicity of solutions for a new p(x)-Kirchhoff problem with variable exponents
  21. An extension of the Hermite-Hadamard inequality for a power of a convex function
  22. Existence and multiplicity of solutions for a fourth-order differential system with instantaneous and non-instantaneous impulses
  23. Relay fusion frames in Banach spaces
  24. Refined ratio monotonicity of the coordinator polynomials of the root lattice of type Bn
  25. On the uniqueness of limit cycles for generalized Liénard systems
  26. A derivative-Hilbert operator acting on Dirichlet spaces
  27. Scheduling equal-length jobs with arbitrary sizes on uniform parallel batch machines
  28. Solutions to a modified gauged Schrödinger equation with Choquard type nonlinearity
  29. A symbolic approach to multiple Hurwitz zeta values at non-positive integers
  30. Some results on the value distribution of differential polynomials
  31. Lucas non-Wieferich primes in arithmetic progressions and the abc conjecture
  32. Scattering properties of Sturm-Liouville equations with sign-alternating weight and transmission condition at turning point
  33. Some results for a p(x)-Kirchhoff type variation-inequality problems in non-divergence form
  34. Homotopy cartesian squares in extriangulated categories
  35. A unified perspective on some autocorrelation measures in different fields: A note
  36. Total Roman domination on the digraphs
  37. Well-posedness for bilevel vector equilibrium problems with variable domination structures
  38. Binet's second formula, Hermite's generalization, and two related identities
  39. Non-solid cone b-metric spaces over Banach algebras and fixed point results of contractions with vector-valued coefficients
  40. Multidimensional sampling-Kantorovich operators in BV-spaces
  41. A self-adaptive inertial extragradient method for a class of split pseudomonotone variational inequality problems
  42. Convergence properties for coordinatewise asymptotically negatively associated random vectors in Hilbert space
  43. Relating the super domination and 2-domination numbers in cactus graphs
  44. Compatibility of the method of brackets with classical integration rules
  45. On the inverse Collatz-Sinogowitz irregularity problem
  46. Positive solutions for boundary value problems of a class of second-order differential equation system
  47. Global analysis and control for a vector-borne epidemic model with multi-edge infection on complex networks
  48. Nonexistence of global solutions to Klein-Gordon equations with variable coefficients power-type nonlinearities
  49. On 2r-ideals in commutative rings with zero-divisors
  50. A comparison of some confidence intervals for a binomial proportion based on a shrinkage estimator
  51. The construction of nuclei for normal constituents of Bπ-characters
  52. Weak solution of non-Newtonian polytropic variational inequality in fresh agricultural product supply chain problem
  53. Mean square exponential stability of stochastic function differential equations in the G-framework
  54. Commutators of Hardy-Littlewood operators on p-adic function spaces with variable exponents
  55. Solitons for the coupled matrix nonlinear Schrödinger-type equations and the related Schrödinger flow
  56. The dual index and dual core generalized inverse
  57. Study on Birkhoff orthogonality and symmetry of matrix operators
  58. Uniqueness theorems of the Hahn difference operator of entire function with a Picard exceptional value
  59. Estimates for certain class of rough generalized Marcinkiewicz functions along submanifolds
  60. On semigroups of transformations that preserve a double direction equivalence
  61. Positive solutions for discrete Minkowski curvature systems of the Lane-Emden type
  62. A multigrid discretization scheme based on the shifted inverse iteration for the Steklov eigenvalue problem in inverse scattering
  63. Existence and nonexistence of solutions for elliptic problems with multiple critical exponents
  64. Interpolation inequalities in generalized Orlicz-Sobolev spaces and applications
  65. General Randić indices of a graph and its line graph
  66. On functional reproducing kernels
  67. On the Waring-Goldbach problem for two squares and four cubes
  68. Singular moduli of rth Roots of modular functions
  69. Classification of self-adjoint domains of odd-order differential operators with matrix theory
  70. On the convergence, stability and data dependence results of the JK iteration process in Banach spaces
  71. Hardy spaces associated with some anisotropic mixed-norm Herz spaces and their applications
  72. Remarks on hyponormal Toeplitz operators with nonharmonic symbols
  73. Complete decomposition of the generalized quaternion groups
  74. Injective and coherent endomorphism rings relative to some matrices
  75. Finite spectrum of fourth-order boundary value problems with boundary and transmission conditions dependent on the spectral parameter
  76. Continued fractions related to a group of linear fractional transformations
  77. Multiplicity of solutions for a class of critical Schrödinger-Poisson systems on the Heisenberg group
  78. Approximate controllability for a stochastic elastic system with structural damping and infinite delay
  79. On extremal cacti with respect to the first degree-based entropy
  80. Compression with wildcards: All exact or all minimal hitting sets
  81. Existence and multiplicity of solutions for a class of p-Kirchhoff-type equation RN
  82. Geometric classifications of k-almost Ricci solitons admitting paracontact metrices
  83. Positive periodic solutions for discrete time-delay hematopoiesis model with impulses
  84. On Hermite-Hadamard-type inequalities for systems of partial differential inequalities in the plane
  85. Existence of solutions for semilinear retarded equations with non-instantaneous impulses, non-local conditions, and infinite delay
  86. On the quadratic residues and their distribution properties
  87. On average theta functions of certain quadratic forms as sums of Eisenstein series
  88. Connected component of positive solutions for one-dimensional p-Laplacian problem with a singular weight
  89. Some identities of degenerate harmonic and degenerate hyperharmonic numbers arising from umbral calculus
  90. Mean ergodic theorems for a sequence of nonexpansive mappings in complete CAT(0) spaces and its applications
  91. On some spaces via topological ideals
  92. Linear maps preserving equivalence or asymptotic equivalence on Banach space
  93. Well-posedness and stability analysis for Timoshenko beam system with Coleman-Gurtin's and Gurtin-Pipkin's thermal laws
  94. On a class of stochastic differential equations driven by the generalized stochastic mixed variational inequalities
  95. Entire solutions of two certain Fermat-type ordinary differential equations
  96. Generalized Lie n-derivations on arbitrary triangular algebras
  97. Markov decision processes approximation with coupled dynamics via Markov deterministic control systems
  98. Notes on pseudodifferential operators commutators and Lipschitz functions
  99. On Graham partitions twisted by the Legendre symbol
  100. Strong limit of processes constructed from a renewal process
  101. Construction of analytical solutions to systems of two stochastic differential equations
  102. Two-distance vertex-distinguishing index of sparse graphs
  103. Regularity and abundance on semigroups of partial transformations with invariant set
  104. Liouville theorems for Kirchhoff-type parabolic equations and system on the Heisenberg group
  105. Spin(8,C)-Higgs pairs over a compact Riemann surface
  106. Properties of locally semi-compact Ir-topological groups
  107. Transcendental entire solutions of several complex product-type nonlinear partial differential equations in ℂ2
  108. Ordering stability of Nash equilibria for a class of differential games
  109. A new reverse half-discrete Hilbert-type inequality with one partial sum involving one derivative function of higher order
  110. About a dubious proof of a correct result about closed Newton Cotes error formulas
  111. Ricci ϕ-invariance on almost cosymplectic three-manifolds
  112. Schur-power convexity of integral mean for convex functions on the coordinates
  113. A characterization of a ∼ admissible congruence on a weakly type B semigroup
  114. On Bohr's inequality for special subclasses of stable starlike harmonic mappings
  115. Properties of meromorphic solutions of first-order differential-difference equations
  116. A double-phase eigenvalue problem with large exponents
  117. On the number of perfect matchings in random polygonal chains
  118. Evolutoids and pedaloids of frontals on timelike surfaces
  119. A series expansion of a logarithmic expression and a decreasing property of the ratio of two logarithmic expressions containing cosine
  120. The 𝔪-WG° inverse in the Minkowski space
  121. Stability result for Lord Shulman swelling porous thermo-elastic soils with distributed delay term
  122. Approximate solvability method for nonlocal impulsive evolution equation
  123. Construction of a functional by a given second-order Ito stochastic equation
  124. Global well-posedness of initial-boundary value problem of fifth-order KdV equation posed on finite interval
  125. On pomonoid of partial transformations of a poset
  126. New fractional integral inequalities via Euler's beta function
  127. An efficient Legendre-Galerkin approximation for the fourth-order equation with singular potential and SSP boundary condition
  128. Eigenfunctions in Finsler Gaussian solitons
  129. On a blow-up criterion for solution of 3D fractional Navier-Stokes-Coriolis equations in Lei-Lin-Gevrey spaces
  130. Some estimates for commutators of sharp maximal function on the p-adic Lebesgue spaces
  131. A preconditioned iterative method for coupled fractional partial differential equation in European option pricing
  132. A digital Jordan surface theorem with respect to a graph connectedness
  133. A quasi-boundary value regularization method for the spherically symmetric backward heat conduction problem
  134. The structure fault tolerance of burnt pancake networks
  135. Average value of the divisor class numbers of real cubic function fields
  136. Uniqueness of exponential polynomials
  137. An application of Hayashi's inequality in numerical integration
Heruntergeladen am 21.10.2025 von https://www.degruyterbrill.com/document/doi/10.1515/math-2022-0581/html
Button zum nach oben scrollen