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

Jurnal Fisika Indonesia

Solution Of The Time Dependent Schrodinger Equation Using The Higher Order Trotter-Suzuki Method

Dhea Chita Octavianty, Pekik Nurwantoro

Published Jul 28, 2022• 0 citations
Publisher

Abstract

We solve the solution of the Time-Dependent Schrodinger Equation (TDSE) using the higher-order Trotter-Suzuki Decomposition. We use The Gaussian function as a stationary state in a linear potential system. The TDSE solution using Baker - Campbell - Hausdorff was used to validate the results and to measure the accuracy of the Trotter - Suzuki decomposition. So that the difference between the Baker - Campbell - Hausdorff result and the Suzuki - Trotter decomposition is considered an error. The error of the TDSE solution by the Trotter - Suzuki second-order decomposition was lower than the first order. Meanwhile, the error of the TDSE solution by second-order hybrid will be lower than the second-order Trotter - Suzuki decomposition when the value of dx=0.1 and 0.05 with dt < 0.0001. The error comparison of these three methods is only valid when time t<1.

DOI
10.22146/jfi.v26i1.66637

Related papers

Time Dependent Schrödinger Equations with Vector Potentials - Numerical Calculations and Visualizations -

Hiroshi Natori, Tomo Munehisa

We study algorithm to calculate the time dependent Schrödinger equation numerically and to visualize the results. We extend the product formula by the Trotter-Suzuki decomposition to cases of Hamiltonians with vector potentials. This method is unconditionally stable and it can be used for a time dependent Hamiltonian so that one can apply it to various problems without any difficulty. After examining numerical errors for time evolutions of one dimensional wave packets, we present the errors under a constant magnetic field. Also we discuss a cyclotron motion in quantum mechanics, where an advantage for the unconditional stability is demonstrated. In order to visualize the results we make animations of the time evolutions, where the phase of the wave function is represented by colors and its absolute value is shown by brightness. Here we suggest a new way to represent a complex number by colors.

Journal of the Physical Society of Japan•Feb 15, 1997•9 citations

Spectral Method for Solving the Time-Dependent Schrödinger Equation on a Non-Uniform Coordinate Grid

М. А. Захаров

Abstract— The fundamental possibility of a numerical solution of the time-dependent Schrödinger equation for a wave function defined on a non-uniform coordinate grid by the spectral method using the fast Fourier transform algorithm is discussed. The method is based on reducing the non-uniform coordinate grid to a uniform one by a non-linear transformation of coordinates and approximating the obtained evolution operator using the Lie–Trotter–Suzuki product formula. Algorithms for the numerical solution of the first and second orders are developed.

Physics of Particles and Nuclei Letters•Apr 1, 2025•0 citations

Efficient Trotter-Suzuki Schemes for Long-time Quantum Dynamics

Marko Maležič, Johann Ostmeyer

Accurately simulating long-time dynamics of many-body systems is a challenge in both classical and quantum computing due to the accumulation of Trotter errors. While low-order Trotter-Suzuki decompositions are straightforward to implement, their rapidly growing error limits access to long-time observables. We present a framework for constructing efficient high-order Trotter-Suzuki schemes by identifying their structure and directly optimizing their parameters over a high-dimensional space. This method enables the discovery of new schemes with significantly improved efficiency compared to traditional constructions, such as those by Suzuki and Yoshida. Based on the theoretical efficiency and practical performance, we recommend two novel highly efficient schemes at $4^{\textrm{th}}$ and $6^{\textrm{th}}$ order. We also demonstrate the effectiveness of these decompositions on the Heisenberg model and the quantum harmonic oscillator, and find that for a fixed final time they perform better across the computational cost. Even when using large time steps, they surpass established low-order schemes like the Leapfrog. Finally, we investigate the in-practice performance of different Trotter schemes and find the decompositions with more uniform coefficients tend to feature improved error accumulation over long times. We have included this observation into our choice of recommended schemes.

quant-pharXiv:2601.18756•Journal of Physics A Mathematical and Theoretical•Jan 27, 2026•1 citation

Quasi Monte-Carlo methods for the TDSE (Time dependent Schrödinger equation)

Gowri Suryanarayana, Dirk Nuyens, Ronald F. A. Cools

We approximate the solution for the time dependent Schrödinger equation (TDSE) in two steps. We first use a pseudo-spectral collocation method that uses samples of the functions on rank-1 lattice points. We then get a system of ordinary differential equations in time, which we solve approximately by stepping in time using the Strang splitting method.

•Nov 23, 2015•0 citations

A METHOD OF GENERATING EXACT SOLUTIONS OF THE TIME-DEPENDENT SCHRÖDINGER EQUATION WITH TIME-DEPENDENT MASS

Axel Schulze‐Halberg

Extending the method presented in our previous paper,12 we map the time-dependent Schrödinger equation (TDSE) with time-dependent mass on a stationary Schrödinger equation for a nonconstant potential. On solving the latter, we can thus generate a large class of exact solutions of the original TDSE. Several examples are given, including potentials of power-law and modified Pöschl–Teller type.

Modern Physics Letters A•Aug 10, 2002•0 citations

A Method to Calculate a Delta Function of the Hamiltonian by the Suzuki-Trotter Decomposition

T. Munehisa, Yasuko Munehisa

We propose a new method to calculate expectation values of a delta function of the Hamiltonian, < Ψ\mid δ(\hat{H} - E)\mid Ψ>. Since the delta function can be replaced with a Gaussian function, we evaluate < Ψ\mid \sqrt{\fracβπ} e^{-β(\hat{H} - E)^2} \mid Ψ> with large βadopting the Suzuki-Trotter decomposition. Errors of the approximate calculations with the finite Trotter number N_t are estimated to be O(1/N_t^K) for the Kth-order decomposition. The distinct advantage of this method is that the convergence is guaranteed even when the state \mid Ψ> contains the eigenstates whose energies spread over the wide range. In this paper we give a full description of our method within the quantum mechanical physics and present the numerical results for the harmonic oscillator problems in one- and three-dimensional space.

arXiv:cond-mat/0001304•Jan 21, 2000•2 citations
We solve the solution of the Time-Dependent Schrodinger Equation (TDSE) using the higher-order Trotter-Suzuki Decomposition. We use The Gaussian function as a stationary state in a linear potential system. The TDSE solution using Baker - Campbell - Hausdorff was used to validate the results and to measure the accuracy of the Trotter - Suzuki decomposition. So that the difference between the Baker - Campbell - Hausdorff result and the Suzuki - Trotter decomposition is considered an error. The error of the TDSE solution by the Trotter - Suzuki second-order decomposition was lower than the first order. Meanwhile, the error of the TDSE solution by second-order hybrid will be lower than the second-order Trotter - Suzuki decomposition when the value of dx=0.1 and 0.05 with dt < 0.0001. The error comparison of these three methods is only valid when time t<1.
We study algorithm to calculate the time dependent Schrödinger equation numerically and to visualize the results. We extend the product formula by the Trotter-Suzuki decomposition to cases of Hamiltonians with vector potentials. This method is unconditionally stable and it can be used for a time dependent Hamiltonian so that one can apply it to various problems without any difficulty. After examining numerical errors for time evolutions of one dimensional wave packets, we present the errors under a constant magnetic field. Also we discuss a cyclotron motion in quantum mechanics, where an advantage for the unconditional stability is demonstrated. In order to visualize the results we make animations of the time evolutions, where the phase of the wave function is represented by colors and its absolute value is shown by brightness. Here we suggest a new way to represent a complex number by colors.
Abstract— The fundamental possibility of a numerical solution of the time-dependent Schrödinger equation for a wave function defined on a non-uniform coordinate grid by the spectral method using the fast Fourier transform algorithm is discussed. The method is based on reducing the non-uniform coordinate grid to a uniform one by a non-linear transformation of coordinates and approximating the obtained evolution operator using the Lie–Trotter–Suzuki product formula. Algorithms for the numerical solution of the first and second orders are developed.
Accurately simulating long-time dynamics of many-body systems is a challenge in both classical and quantum computing due to the accumulation of Trotter errors. While low-order Trotter-Suzuki decompositions are straightforward to implement, their rapidly growing error limits access to long-time observables. We present a framework for constructing efficient high-order Trotter-Suzuki schemes by identifying their structure and directly optimizing their parameters over a high-dimensional space. This method enables the discovery of new schemes with significantly improved efficiency compared to traditional constructions, such as those by Suzuki and Yoshida. Based on the theoretical efficiency and practical performance, we recommend two novel highly efficient schemes at 4th4^{\textrm{th}}4th and 6th6^{\textrm{th}}6th order. We also demonstrate the effectiveness of these decompositions on the Heisenberg model and the quantum harmonic oscillator, and find that for a fixed final time they perform better across the computational cost. Even when using large time steps, they surpass established low-order schemes like the Leapfrog. Finally, we investigate the in-practice performance of different Trotter schemes and find the decompositions with more uniform coefficients tend to feature improved error accumulation over long times. We have included this observation into our choice of recommended schemes.
We propose a new method to calculate expectation values of a delta function of the Hamiltonian, < Ψ\mid δ(\hat{H} - E)\mid Ψ>. Since the delta function can be replaced with a Gaussian function, we evaluate < Ψ\mid \sqrt{\fracβπ} e^{-β(\hat{H} - E)^2} \mid Ψ> with large βadopting the Suzuki-Trotter decomposition. Errors of the approximate calculations with the finite Trotter number N_t are estimated to be O(1/N_t^K) for the Kth-order decomposition. The distinct advantage of this method is that the convergence is guaranteed even when the state \mid Ψ> contains the eigenstates whose energies spread over the wide range. In this paper we give a full description of our method within the quantum mechanical physics and present the numerical results for the harmonic oscillator problems in one- and three-dimensional space.