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 Applied Mathematics and Physics

Wavelet Basis and Fourier Transform Approach to Path Integral Formulation of Electron Dynamics

Tommi Höynälänmaa

Published Jan 1, 2025• 0 citations
Publisher

Abstract

We present a method for the path integral formulation of electronic structure simulation. The time evolution operator is represented as matrix in a basis consisting of Deslauriers-Dubuc or Daubechies wavelets. We present an approximation of the path integral kernel and a method for calculating wavefunctions. The kernel is tested by finding wavefunctions and eigenenergies of one-dimensional and three-dimensional harmonic oscillators and hydrogen atom in one and three dimensions.

DOI
10.4236/jamp.2025.1310186

Related papers

A path-integration calculation method based on the real-space finite-difference scheme

Hidekazu Goto, Tomoya Ono, Kikuji Hirose

We propose a new path-integration calculation method to treat the time evolution of a wavefunction within the framework of the real-space finite-difference formalism, and also develop an effective scheme to compute the scattering wavefunction for an incident electron with arbitrary energy, in which an impulse wavefunction is adopted as an initial state of the time evolution. In this method, once the time evolution of the initial impulse wavefunction is calculated, all of the solutions in the scattering problem can be derived by means of Fourier analysis of the time-evolved wavefunction, which leads to a reduction of the calculation time. In order to test the applicability of our newly developed simulation procedures, we implemented simulations for the one-dimensional scattering problem. Each simulation showed the usefulness of the present scheme by yielding the steady scattering states in agreement with exact ones.

Journal of Physics Condensed Matter•Aug 24, 2007•2 citations

Path integral formulation of dissipative quantum dynamics

Alexey Novikov

In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.

•Jan 20, 2005•0 citations

Path integral density matrix dynamics: A method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems

Scott Habershon

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.

The Journal of Chemical Physics•Sep 12, 2013•8 citations

Numerical Path Integral Approach to Quantum Dynamics and Stationary Quantum States

Ilkka Ruokosenmäki, Tapio T. Rantala

Applicability of Feynman path integral approach to numerical simulations of quantum dynamics of an electron in real time domain is examined. Coherent quantum dynamics is demonstrated with one dimensional test cases (quantum dot models) and performance of the Trotter kernel as compared with the exact kernels is tested. Also, a novel approach for finding the ground state and other stationary sates is presented. This is based on the incoherent propagation in real time. For both approaches the Monte Carlo grid and sampling are tested and compared with regular grids and sampling. We asses the numerical prerequisites for all of the above.

arXiv:1403.1448•Communications in Computational Physics•Jan 1, 2015•0 citations

The Feynman path integral, from the harmonic oscillator to quantum electrodynamics

Rafael Carlos Silva Tonhon

The Feynman path-integral formulation of quantum mechanics is developed. The math- ematical framework is worked out from Schrödinger's equation, illustrating the compati- bility of the two formulations. We introduce the sum over all paths, known as the kernel K, and discuss its implications for events occurring in succession. We compute the kernel for quadratic Lagrangians. Then, making use of imaginary time, we consider the discrete version of the kernel, known as the Euclidean path integral, and use it to solve the (quan- tum) harmonic and anharmonic oscillators through numerical Monte Carlo simulations. Lastly, we use the framework to quantize the electromagnetic field, as a motivation to discuss perturbation theory via Feynman diagrams.

•Jan 1, 2022•0 citations

Mathematical aspects of the Feynman path integral

Daniel W. Boutros

Path integration methods are of crucial importance to quantum mechanics and quantum field theory. There are multiple ways the path integral can be constructed, one method uses the link between the Fokker-Planck and the Langevin equations, as is covered in this thesis. This is related to the standard derivation of the path integral in physics, in which the time is discretised. We study the mathematical problems related to the Feynman path integral, in particular the impossibility of a Lebesgue-type measure on the space of paths. It is discussed how oscillatory integrals can be used to have a well-defined `integral' on an infinite dimensional space. This formalism is subsequently applied to both the harmonic and the anharmonic oscillator. We prove an infinite dimensional oscillatory integral exists and obtain a convergent expression for both these cases, under conditions. Examples of these conditions are the initial wavefunctions being Schwartz functions as well as conditions on the endtime, angular frequency and coupling constant. Finally, a possible approach to establish an oscillatory integral for the hydrogen atom is discussed. It is proven that the result is no longer independent of the sequence of projection operators, which is a key step towards a rigorous path integral for this system.

•Jan 1, 2020•0 citations
We present a method for the path integral formulation of electronic structure simulation. The time evolution operator is represented as matrix in a basis consisting of Deslauriers-Dubuc or Daubechies wavelets. We present an approximation of the path integral kernel and a method for calculating wavefunctions. The kernel is tested by finding wavefunctions and eigenenergies of one-dimensional and three-dimensional harmonic oscillators and hydrogen atom in one and three dimensions.
We propose a new path-integration calculation method to treat the time evolution of a wavefunction within the framework of the real-space finite-difference formalism, and also develop an effective scheme to compute the scattering wavefunction for an incident electron with arbitrary energy, in which an impulse wavefunction is adopted as an initial state of the time evolution. In this method, once the time evolution of the initial impulse wavefunction is calculated, all of the solutions in the scattering problem can be derived by means of Fourier analysis of the time-evolved wavefunction, which leads to a reduction of the calculation time. In order to test the applicability of our newly developed simulation procedures, we implemented simulations for the one-dimensional scattering problem. Each simulation showed the usefulness of the present scheme by yielding the steady scattering states in agreement with exact ones.
In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.
We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.
Applicability of Feynman path integral approach to numerical simulations of quantum dynamics of an electron in real time domain is examined. Coherent quantum dynamics is demonstrated with one dimensional test cases (quantum dot models) and performance of the Trotter kernel as compared with the exact kernels is tested. Also, a novel approach for finding the ground state and other stationary sates is presented. This is based on the incoherent propagation in real time. For both approaches the Monte Carlo grid and sampling are tested and compared with regular grids and sampling. We asses the numerical prerequisites for all of the above.
The Feynman path-integral formulation of quantum mechanics is developed. The math- ematical framework is worked out from Schrödinger's equation, illustrating the compati- bility of the two formulations. We introduce the sum over all paths, known as the kernel K, and discuss its implications for events occurring in succession. We compute the kernel for quadratic Lagrangians. Then, making use of imaginary time, we consider the discrete version of the kernel, known as the Euclidean path integral, and use it to solve the (quan- tum) harmonic and anharmonic oscillators through numerical Monte Carlo simulations. Lastly, we use the framework to quantize the electromagnetic field, as a motivation to discuss perturbation theory via Feynman diagrams.
Path integration methods are of crucial importance to quantum mechanics and quantum field theory. There are multiple ways the path integral can be constructed, one method uses the link between the Fokker-Planck and the Langevin equations, as is covered in this thesis. This is related to the standard derivation of the path integral in physics, in which the time is discretised. We study the mathematical problems related to the Feynman path integral, in particular the impossibility of a Lebesgue-type measure on the space of paths. It is discussed how oscillatory integrals can be used to have a well-defined `integral' on an infinite dimensional space. This formalism is subsequently applied to both the harmonic and the anharmonic oscillator. We prove an infinite dimensional oscillatory integral exists and obtain a convergent expression for both these cases, under conditions. Examples of these conditions are the initial wavefunctions being Schwartz functions as well as conditions on the endtime, angular frequency and coupling constant. Finally, a possible approach to establish an oscillatory integral for the hydrogen atom is discussed. It is proven that the result is no longer independent of the sequence of projection operators, which is a key step towards a rigorous path integral for this system.