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

Molecular Physics

Exploring a spectral filtering approach to electronic structure calculations

Alexis J. Kiessling, Jeffrey A. Cina

Published Oct 5, 2020• 0 citations
Publisher

Abstract

We make an initial exploration of a spectral filtering approach to numerically approximating solutions to the time-independent Schrödinger equation for a model system of two spinless fermions with harmonic interactions. Fourier transformation of the time evolution of an arbitrary antisymmetrized wave packet results in resonance at frequencies proportional to eigenenergies of the system Hamiltonian, while the Fourier coefficient approximates the corresponding eigenfunction. Spectral filtering is explored using the two-particle model through both direct numerical integration of the Schrödinger equation and a semiclassical parametrisation. In the former, a discrete position basis describes the time-evolution operator, which is then repeatedly applied to an antisymmetrized wave packet to obtain the packet's time evolution. In the semiclassical parametrisation, equations of motion for the parameters are obtained from the Dirac-Frenkel-McLachlan functional. Integration of the parameter equations of motion describes the packet's time evolution. The semiclassical approach requires less memory and accurately obtains the eigenspectrum. In a step towards application of spectral filtering to realistic systems, a possible ansatz for use with Coulomb potentials is explored, using a simple model Hamiltonian.

DOI
10.1080/00268976.2020.1827178

Related papers

Representation of electrons on symmetric electron-wave stub-filters by waves and particles

Shota Hayakawa, Nobuyuki Hirami, Hisato Fujisaka

This paper presents a method of computing sample trajectories of particle models of electrons on symmetric electron-wave stub-filters. The time-independent Schrödinger equation describing the electron-wave filters is transformed to a system of linear algebraic equations. The Fourier series expansion transforms the eigenvalues of the system to continuous and differentiable eigenfunctions. A electron wave packet is constructed from the obtained smooth eigenfunctions and a nonlinear stochastic ordinary differential equation (NSODE) which governs the behavior of the electron particle model is established by applying the stochastic quantization to the wave packet. Sample trajectories of the particles are computed by integrating the NSODE. Numerical experiments have shown that the statistical behavior of the particles is almost equal to that of the wave packets.

Nonlinear Theory and Its Applications IEICE•Jan 1, 2019•3 citations

Energy eigenvalue spectroscopy

S. M. Blinder

Abstract A wave function which is other than an exact eigenfunction, if it obeys appropriate analytical conditions, can be considered to represent the initial configuration of a nonstationary state. In the course of its subsequent time development the quantum system exhibits implicitly its entire eigenvalue spectrum. A method based, in principle, on Fourier analysis of the evolving quantum system is applied to the direct calculation of the energy eigenvalue spectrum. The spectral function is expressed as a moment expansion, in terms of expectation values of powers of the Hamiltonian. When the expansion is truncated, as it must be in any practical application, the corresponding spectral function represents a smeared‐out eigenvalue spectrum. An alternative approximation leads to the quantum‐mechanical method of moments. As the number of terms is increased, the computed spectrum becomes sharper and more accurate. In certain cases the moment expansion can be circumvented, if the action of the evolution operator can be evaluated in closed form. This is equivalent to finding some solution of the time‐dependent Schrödinger equation. The various methods of eigenvalue spectroscopy are applied to the harmonic oscillator.

International Journal of Quantum Chemistry•May 1, 1967•9 citations

A strategy for time dependent quantum mechanical calculations using a Gaussian wave packet representation of the wave function

Shin‐ichi Sawada, Robert Heather, Bret Jackson +1 more

We develop methodology for performing time dependent quantum mechanical calculations by representing the wave function as a sum of Gaussian wave packets (GWP) each characterized by a set of parameters such as width, position, momentum, and phase. The problem of computing the time evolution of the wave function is thus reduced to that of finding the time evolution of the parameters in the Gaussians. This parameter motion is determined by minimizing the error made by replacing the exact wave function in the time dependent Schrödinger equation with its Gaussian representation approximant. This leads to first order differential equations for the time dependence of the parameters, and those describing the packet position and the momentum of each packet have some resemblance with the classical equations of motion. The paper studies numerically the strategy needed to achieve the best GWP representation of time dependent processes. The issues discussed are the representation of the initial wave function, the numerical stability and the solution of the differential equations giving the evolution of the parameters, and the analysis of the final wave function. Extensive comparisons are made with an approximate method which assumes that the Gaussians are independent and their width is smaller than the length scale over which the potential changes. This approximation greatly simplifies the calculations and has the advantage of a greater resemblance to classical mechanics, thus being more intuitive. We find, however, that its range of applications is limited to problems involving localized degrees of freedom that participate in the dynamic process for a very short time. Finally we give particular attention to the notion that the GWP representation of the wave function reduces the dynamics of one quantum degree of freedom to that of a set of pseudoparticles (each represented by one packet) moving according to a pseudoclassical (i.e., classical-like) mechanics whose ‘‘phase space’’ is described by a position and momentum as well as a complex phase and width.

The Journal of Chemical Physics•Sep 15, 1985•205 citations

Variational Gaussian approximation for the magnetic Schrödinger equation

Selina Burkhard, Benjamin Dörich, Marlis Hochbruck +1 more

In the present paper we consider the semiclassical magnetic Schrödinger equation, which describes the dynamics of particles under the influence of a magnetic field. The solution of the Schrödinger equation is approximated by Gaussian wave packets via the time-dependent variational formulation by Dirac and Frenkel. For the numerical approximation we derive ordinary differential equations for the parameters of the variational solution. Moreover, we prove $L^2$-error bounds and observable error bounds for the approximating Gaussian wave packet.

arXiv:2310.16734•Jan 1, 2023•0 citations

Time-integration of Gaussian variational approximation for the magnetic Schrödinger equation

Malik Scheifinger, Kurt Busch, Marlis Hochbruck +1 more

In the present paper we consider the semiclassical magnetic Schrödinger equation, which describes the dynamics of charged particles under the influence of a electro-magnetic field. The solution of the time-dependent Schrödinger equation is approximated by a single Gaussian wave packet via the time-dependent Dirac--Frenkel variational principle. For the approximation we use ordinary differential equations of motion for the parameters of the variational solution and extend the second-order Boris algorithm for classical mechanics to the quantum mechanical case. In addition, we propose a modified version of the classical fourth order Runge--Kutta method. Numerical experiments explore parameter convergence and geometric properties. Moreover, we benchmark against the analytical solution of the Penning trap.

arXiv:2504.03407•Apr 4, 2025•0 citations

Time-integration of Gaussian variational approximation for the magnetic Schrödinger equation

Malik Scheifinger, Kurt Busch, Marlis Hochbruck +1 more

In the present paper we consider the semiclassical magnetic Schrödinger equation, which describes the dynamics of charged particles under the influence of an electro-magnetic field. The solution of the time-dependent Schrödinger equation is approximated by a single Gaussian wave packet via the time-dependent Dirac–Frenkel variational principle. For the approximation we use ordinary differential equations of motion for the parameters of the variational solution and extend the second-order Boris algorithm for classical mechanics to the quantum mechanical case. In addition, we propose a modified version of the classical fourth-order Runge–Kutta method. Numerical experiments explore parameter convergence and geometric properties. Moreover, we benchmark against the analytical solution of the Penning trap.

Journal of Computational Physics•Sep 5, 2025•0 citations
We make an initial exploration of a spectral filtering approach to numerically approximating solutions to the time-independent Schrödinger equation for a model system of two spinless fermions with harmonic interactions. Fourier transformation of the time evolution of an arbitrary antisymmetrized wave packet results in resonance at frequencies proportional to eigenenergies of the system Hamiltonian, while the Fourier coefficient approximates the corresponding eigenfunction. Spectral filtering is explored using the two-particle model through both direct numerical integration of the Schrödinger equation and a semiclassical parametrisation. In the former, a discrete position basis describes the time-evolution operator, which is then repeatedly applied to an antisymmetrized wave packet to obtain the packet's time evolution. In the semiclassical parametrisation, equations of motion for the parameters are obtained from the Dirac-Frenkel-McLachlan functional. Integration of the parameter equations of motion describes the packet's time evolution. The semiclassical approach requires less memory and accurately obtains the eigenspectrum. In a step towards application of spectral filtering to realistic systems, a possible ansatz for use with Coulomb potentials is explored, using a simple model Hamiltonian.
This paper presents a method of computing sample trajectories of particle models of electrons on symmetric electron-wave stub-filters. The time-independent Schrödinger equation describing the electron-wave filters is transformed to a system of linear algebraic equations. The Fourier series expansion transforms the eigenvalues of the system to continuous and differentiable eigenfunctions. A electron wave packet is constructed from the obtained smooth eigenfunctions and a nonlinear stochastic ordinary differential equation (NSODE) which governs the behavior of the electron particle model is established by applying the stochastic quantization to the wave packet. Sample trajectories of the particles are computed by integrating the NSODE. Numerical experiments have shown that the statistical behavior of the particles is almost equal to that of the wave packets.
Abstract A wave function which is other than an exact eigenfunction, if it obeys appropriate analytical conditions, can be considered to represent the initial configuration of a nonstationary state. In the course of its subsequent time development the quantum system exhibits implicitly its entire eigenvalue spectrum. A method based, in principle, on Fourier analysis of the evolving quantum system is applied to the direct calculation of the energy eigenvalue spectrum. The spectral function is expressed as a moment expansion, in terms of expectation values of powers of the Hamiltonian. When the expansion is truncated, as it must be in any practical application, the corresponding spectral function represents a smeared‐out eigenvalue spectrum. An alternative approximation leads to the quantum‐mechanical method of moments. As the number of terms is increased, the computed spectrum becomes sharper and more accurate. In certain cases the moment expansion can be circumvented, if the action of the evolution operator can be evaluated in closed form. This is equivalent to finding some solution of the time‐dependent Schrödinger equation. The various methods of eigenvalue spectroscopy are applied to the harmonic oscillator.
We develop methodology for performing time dependent quantum mechanical calculations by representing the wave function as a sum of Gaussian wave packets (GWP) each characterized by a set of parameters such as width, position, momentum, and phase. The problem of computing the time evolution of the wave function is thus reduced to that of finding the time evolution of the parameters in the Gaussians. This parameter motion is determined by minimizing the error made by replacing the exact wave function in the time dependent Schrödinger equation with its Gaussian representation approximant. This leads to first order differential equations for the time dependence of the parameters, and those describing the packet position and the momentum of each packet have some resemblance with the classical equations of motion. The paper studies numerically the strategy needed to achieve the best GWP representation of time dependent processes. The issues discussed are the representation of the initial wave function, the numerical stability and the solution of the differential equations giving the evolution of the parameters, and the analysis of the final wave function. Extensive comparisons are made with an approximate method which assumes that the Gaussians are independent and their width is smaller than the length scale over which the potential changes. This approximation greatly simplifies the calculations and has the advantage of a greater resemblance to classical mechanics, thus being more intuitive. We find, however, that its range of applications is limited to problems involving localized degrees of freedom that participate in the dynamic process for a very short time. Finally we give particular attention to the notion that the GWP representation of the wave function reduces the dynamics of one quantum degree of freedom to that of a set of pseudoparticles (each represented by one packet) moving according to a pseudoclassical (i.e., classical-like) mechanics whose ‘‘phase space’’ is described by a position and momentum as well as a complex phase and width.
In the present paper we consider the semiclassical magnetic Schrödinger equation, which describes the dynamics of particles under the influence of a magnetic field. The solution of the Schrödinger equation is approximated by Gaussian wave packets via the time-dependent variational formulation by Dirac and Frenkel. For the numerical approximation we derive ordinary differential equations for the parameters of the variational solution. Moreover, we prove L2L^2L2-error bounds and observable error bounds for the approximating Gaussian wave packet.
In the present paper we consider the semiclassical magnetic Schrödinger equation, which describes the dynamics of charged particles under the influence of a electro-magnetic field. The solution of the time-dependent Schrödinger equation is approximated by a single Gaussian wave packet via the time-dependent Dirac--Frenkel variational principle. For the approximation we use ordinary differential equations of motion for the parameters of the variational solution and extend the second-order Boris algorithm for classical mechanics to the quantum mechanical case. In addition, we propose a modified version of the classical fourth order Runge--Kutta method. Numerical experiments explore parameter convergence and geometric properties. Moreover, we benchmark against the analytical solution of the Penning trap.
In the present paper we consider the semiclassical magnetic Schrödinger equation, which describes the dynamics of charged particles under the influence of an electro-magnetic field. The solution of the time-dependent Schrödinger equation is approximated by a single Gaussian wave packet via the time-dependent Dirac–Frenkel variational principle. For the approximation we use ordinary differential equations of motion for the parameters of the variational solution and extend the second-order Boris algorithm for classical mechanics to the quantum mechanical case. In addition, we propose a modified version of the classical fourth-order Runge–Kutta method. Numerical experiments explore parameter convergence and geometric properties. Moreover, we benchmark against the analytical solution of the Penning trap.