Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Physics Conference Series

Energy levels and expectation values via accelerated path integral Monte Carlo

Danica Stojiljković, Aleksandar Bogojević, Antun Balaž

Published Aug 1, 2008• 1 citation
Publisher

Abstract

A recently developed method systematically improved the convergence of generic path integrals for transition amplitudes, partition functions, expectation values and energy spectra. This was achieved by analytically constructing a hierarchy of discretized effective actions indexed by a level number p and converging to the continuum limit as 1/N p . Here we apply the above general method to numerical calculations using Metropolis Monte Carlo simulations of energy expectation values and energy spectra. We analyze and compare the ensuing increase in efficiency of several orders of magnitude.

DOI
10.1088/1742-6596/128/1/012062

Related papers

Accelerated path integral calculations for many-body systems

Antun Balaž, Ivana Vidanović, Aleksandar Bogojević

We present a generalization of a recently developed method for accelerated Monte Carlo calculation of path integrals of generic many particle systems. Using this method and the derived hierarchy of efiective actions we calculate the energy spectra of a two particle model with quartic interaction for several values of coupling and demonstrate agreement with analytical results governing the increase in efficiency of the new method.

Journal of Physics Conference Series•Aug 1, 2008•0 citations

Recursive Schrödinger equation approach to faster converging path integrals

Antun Balaž, Aleksandar Bogojević, Ivana Vidanović +1 more

By recursively solving the underlying Schrödinger equation, we set up an efficient systematic approach for deriving analytic expressions for discretized effective actions. With this, we obtain discrete short-time propagators for both one and many particles in arbitrary dimension to orders that have not been accessible before. They can be used to substantially speed up numerical Monte Carlo calculations of path integrals, as well as for setting up an alternative analytical approximation scheme for energy spectra, density of states, and other statistical properties of quantum systems.

arXiv:0806.4774•Physical Review E•Mar 9, 2009•25 citations

SPEEDUP Code for Calculation of Transition Amplitudes via the Effective Action Approach

Antun Balaž, Ivana Vasić, Danica Stojiljković +3 more

We present Path Integral Monte Carlo C code for calculation of quantum mechanical transition amplitudes for 1D models. The SPEEDUP C code is based on the use of higher-order short-time effective actions and implemented to the maximal order $p$=18 in the time of propagation (Monte Carlo time step), which substantially improves the convergence of discretized amplitudes to their exact continuum values. Symbolic derivation of higher-order effective actions is implemented in SPEEDUP Mathematica codes, using the recursive Schrödinger equation approach. In addition to the general 1D quantum theory, developed Mathematica codes are capable of calculating effective actions for specific models, for general 2D and 3D potentials, as well as for a general many-body theory in arbitrary number of spatial dimensions.

arXiv:1105.0542•Communications in Computational Physics•May 22, 2012•6 citations

The Feynman path integral

Mark E. Tuckerman

Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations

•Aug 1, 2023•0 citations

Fast Converging Path Integrals for Time-Dependent Potentials I: Recursive Calculation of Short-Time Expansion of the Propagator

Antun Balaž, Ivana Vidanović, Aleksandar Bogojević +2 more

In this and subsequent paper [1] we develop a recursive approach for calculating the short-time expansion of the propagator for a general quantum system in a time-dependent potential to orders that have not yet been accessible before. To this end the propagator is expressed in terms of a discretized effective potential, for which we derive and analytically solve a set of efficient recursion relations. Such a discretized effective potential can be used to substantially speed up numerical Monte Carlo simulations for path integrals, or to set up various analytic approximation techniques to study properties of quantum systems in time-dependent potentials. The analytically derived results are numerically verified by treating several simple models.

•Jan 1, 2011•12 citations

Fast converging path integrals for time-dependent potentials: I. Recursive calculation of short-time expansion of the propagator

Antun Balaž, Ivana Vidanović, Aleksandar Bogojević +2 more

In this and the subsequent paper (Balaž et al 2011 J. Stat. Mech. P03005) we develop a recursive approach for calculating the short-time expansion of the propagator for a general quantum system in a time-dependent potential to orders that have not been accessible before. To this end the propagator is expressed in terms of a discretized effective potential, for which we derive and analytically solve a set of efficient recursion relations. Such a discretized effective potential can be used to substantially speed up numerical Monte Carlo simulations for path integrals, or to set up various analytic approximation techniques to study properties of quantum systems in time-dependent potentials. The analytically derived results are numerically verified by treating several simple models.

arXiv:0912.2743•Journal of Statistical Mechanics Theory and Experiment•Mar 4, 2011•9 citations
A recently developed method systematically improved the convergence of generic path integrals for transition amplitudes, partition functions, expectation values and energy spectra. This was achieved by analytically constructing a hierarchy of discretized effective actions indexed by a level number p and converging to the continuum limit as 1/N p . Here we apply the above general method to numerical calculations using Metropolis Monte Carlo simulations of energy expectation values and energy spectra. We analyze and compare the ensuing increase in efficiency of several orders of magnitude.
By recursively solving the underlying Schrödinger equation, we set up an efficient systematic approach for deriving analytic expressions for discretized effective actions. With this, we obtain discrete short-time propagators for both one and many particles in arbitrary dimension to orders that have not been accessible before. They can be used to substantially speed up numerical Monte Carlo calculations of path integrals, as well as for setting up an alternative analytical approximation scheme for energy spectra, density of states, and other statistical properties of quantum systems.
We present Path Integral Monte Carlo C code for calculation of quantum mechanical transition amplitudes for 1D models. The SPEEDUP C code is based on the use of higher-order short-time effective actions and implemented to the maximal order ppp=18 in the time of propagation (Monte Carlo time step), which substantially improves the convergence of discretized amplitudes to their exact continuum values. Symbolic derivation of higher-order effective actions is implemented in SPEEDUP Mathematica codes, using the recursive Schrödinger equation approach. In addition to the general 1D quantum theory, developed Mathematica codes are capable of calculating effective actions for specific models, for general 2D and 3D potentials, as well as for a general many-body theory in arbitrary number of spatial dimensions.
Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations
In this and subsequent paper [1] we develop a recursive approach for calculating the short-time expansion of the propagator for a general quantum system in a time-dependent potential to orders that have not yet been accessible before. To this end the propagator is expressed in terms of a discretized effective potential, for which we derive and analytically solve a set of efficient recursion relations. Such a discretized effective potential can be used to substantially speed up numerical Monte Carlo simulations for path integrals, or to set up various analytic approximation techniques to study properties of quantum systems in time-dependent potentials. The analytically derived results are numerically verified by treating several simple models.
In this and the subsequent paper (Balaž et al 2011 J. Stat. Mech. P03005) we develop a recursive approach for calculating the short-time expansion of the propagator for a general quantum system in a time-dependent potential to orders that have not been accessible before. To this end the propagator is expressed in terms of a discretized effective potential, for which we derive and analytically solve a set of efficient recursion relations. Such a discretized effective potential can be used to substantially speed up numerical Monte Carlo simulations for path integrals, or to set up various analytic approximation techniques to study properties of quantum systems in time-dependent potentials. The analytically derived results are numerically verified by treating several simple models.