Quadpack computation of Feynman loop integrals

The paper addresses a numerical computation of Feynman loop integrals, which are computed by an extrapolation to the limit as a parameter in the integrand tends to zero. An important objective is to achieve an automatic computation which is effective for a wide range of instances. Singular or near singular integrand behavior is handled via an adaptive partitioning of the domain, implemented in an iterated/repeated multivariate integration method. Integrand singularities possibly introduced via infrared (IR) divergence at the boundaries of the integration domain are addressed using a version of the Dqags algorithm from the integration package Quadpack, which uses an adaptive strategy combined with extrapolation. The latter is justified for a large class of problems by the underlying asymptotic expansions of the integration error. For IR divergent problems, an extrapolation scheme is presented based on dimensional regularization.     

Computation of Feynman integrals of functionals that are functions of quadratic functionals

If the arguments of a function G: ^{d 1} are taken as quadratic functional defined on a space C of continuous functions, we obtain a functional G: C ¹.We give a formula for computing analytic Feynman integrals of such functionals. We also propose a method of approximate computation of sequential Feynman integrals based on replacing the kernel of an integral operator by a degenerate kernel.Translated fromMatematicheskie Metody i Fiziko-Mekhanicheskie Polya, Issue 32, 1990, pp. 58–61.     

Transformations of Feynman path integrals and generalized densities of Feynman pseudomeasures

Applications of transformations of Feynman path integrals and Feynman pseudomeasures to explain arising quantum anomalies are considered. A contradiction in the literature is also explained.     

Operator-valued Feynman integral via conditional Feynman integrals on a function space

Let C ₀ ^r [0; t] denote the analogue of the r-dimensional Wiener space, define X _t : C ^r [0; t] → ?^2r by X _t (x) = (x(0); x(t)). In this paper, we introduce a simple formula for the conditional expectations with the conditioning function X _t . Using this formula, we evaluate the conditional analytic Feynman integral for the functional $$ \Gamma _t \left( x \right) = exp \left\{ {\int_0^t {\theta \left( {s,x\left( s \right)} \right)d\eta \left( s \right)} } \right\}\varphi \left( {x\left( t \right)} \right) x \in C^r \left[ {0,t} \right] $$ , where η is a complex Borel measure on [0, t], and θ(s, ·) and φ are the Fourier-Stieltjes transforms of the complex Borel measures on ? ^r . We then introduce an integral transform as an analytic operator-valued Feynman integral over C ^r [0, t], and evaluate the integral transform for the function Γ _t via the conditional analytic Feynman integral as a kernel.     

A method of calculating massive Feynman integrals

Scientific-Research Institute of Nuclear Physics at the Moscow State University. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 89, No. 1, pp. 56–72, October, 1991.     

Feynman path integrals as infinite-dimensional oscillatory integrals: Some new developments

We discuss a basic mathematical approach to Feynman path integrals as infinite-dimensional oscillatory integrals. We present new results on asymptotics of such integrals which exploit recently developed approximation techniques via finite dimensional oscillatory integrals. Applications are also given, namely to the study of the trace of the time evolution operator in quantum mechanics and to the interpretation of Gutzwiller's trace formula as a leading term in an asymptotic expansion around classical periodic orbits.The second named author is an Alexander von Humboldt Stiftung fellow.     

Feynman path integrals for polynomially growing potentials

A general class of infinite dimensional oscillatory integrals with polynomially growing phase functions is studied. A representation formula of the Parseval type is proved, as well as a formula giving the integrals in terms of analytically continued absolutely convergent integrals. These results are applied to provide a rigorous Feynman path integral representation for the solution of the time-dependent Schrödinger equation with a quartic anharmonic potential. The Borel summability of the asymptotic expansion of the solution in power series of the coupling constant is also proved.     

Computation of exponential integrals

Let P ^d be a convex polyhedron and f: ^d a linear function. One studies the computational complexity of the integral _pexp f(xdx. It is shown that these integrals satisfy nontrivial algebraic relations, which makes possible the construction of polynomial algorithms for certain polyhedra. Examples are given of the application of exponential integrals to the calculation of volume and nonlinear programming.Translated from Zapiski Nauchnykh Seminarov Leningradksogo Otdeleniya Matematicheskogo Instituta im. V. A. Skeklova AN SSSR, Vol. 192, pp. 149–162, 1991.     

Computation of oscillating integrals

An adaptive quadrature method for the automatic computation of integrals with strongly oscillating integrand is presented. The integration method is based on a truncated Chebyshev series approximation. The algorithm uses a global subinterval division strategy. There is a protection against the influence of round-off errors. A Fortran implementation of the algorithm is given.     

Relationships between generalized Wiener integrals and conditional analytic Feynman integrals over continuous paths

Let C[0, t] denote a generalized Wiener space, the space of real-valued continuous functions on the interval [0, t], and define a random vector Z _n: C[0, t] → R ⁿ⁺¹ by \({Z_n}\left( x \right) = \left( {x\left( 0 \right) + a\left( 0 \right),\int_o^{{t_1}} {h\left( s \right)dx\left( s \right) + x\left( 0 \right) + a\left( {{t_1}} \right),...,\int_0^{{t_n}} {h\left( s \right)dx\left( s \right) + x\left( 0 \right) + a\left( {{t_n}} \right)} } } \right)\), where a ∈ C[0, t], h ∈ L ₂[0, t], and 0 < t ₁ <... < t _n ≤ t is a partition of [0, t]. Using simple formulas for generalized conditional Wiener integrals, given Z _n we will evaluate the generalized analytic conditional Wiener and Feynman integrals of the functions F in a Banach algebra which corresponds to Cameron-Storvick’s Banach algebra S. Finally, we express the generalized analytic conditional Feynman integral of F as a limit of the non-conditional generalized Wiener integral of a polygonal function using a change of scale transformation for which a normal density is the kernel. This result extends the existing change of scale formulas on the classical Wiener space, abstract Wiener space and the analogue of the Wiener space C[0, t].     

Computation of poisson type integrals

We consider problems of computing the Poisson integral when the point at which the integral is evaluated approaches the ball surface. Techniques are proposed enabling one to improve the computation efficiency.