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 Chemical Education

Observing Quantum Dynamics in Real Time: An Excel-Ready Finite-Difference Algorithm for Solving the Time-Dependent Schrödinger Equation

Taylor A. Rossman, Zachary P. Parks, Michael G. Messina

Published Mar 11, 2020• 4 citations
Publisher

Abstract

We present an algorithm for solving the time-dependent Schrödinger equation that is based on the finite-difference expression of the kinetic energy operator. Students who have some knowledge of linear algebra can understand the theory used to derive the algorithm. This is because the finite-difference kinetic energy matrix and the Hückel matrix for linear conjugated hydrocarbons have similar forms. We offer 3 modules that can be used to solve 1-dimensional time-dependent Schrödinger equations, in real time. Module 1 contains two worksheets designed to show quantum dynamics in bound potentials. Module 2 contains 2 worksheets designed to illustrate tunneling. Module 3 contains 2 worksheets set up to illustrate the effect of time-dependent external fields. Further, modules 1 and 3 also include the Verlet algorithm so students can compare the quantum dynamical evolution of the systems to their classical counterpart in real time.

DOI
10.1021/acs.jchemed.9b01138

Related papers

Time Dependent Schrödinger Equation

Pawel Scharoch, Maciej P. Polak, Radosław Szymon +1 more

This project looks into the time evolution of a wave function within a two-dimensional quantum well. We start by solving the time-dependent Schrödinger equation for stationary states in a quantum well. Next, we express any wave function as a linear combination of stationary states, allowing us to understand their time evolution. Two methods are presented: one relies on decomposing the wave function into a basis of stationary states and the other on discretisation of the time-dependent Schrödinger equation, incorporating three-point formulas for derivatives. These approaches necessitate confronting intricate boundary conditions and require maintaining energy conservation for numerical accuracy. We further demonstrate the methods using a wave packet, revealing fundamental phenomena in quantum physics. Our results demonstrate the utility of these methods in understanding quantum systems, despite the challenges in determining stationary states for a given potential. This study enhances our comprehension of the dynamics of quantum states in constrained systems, essential for fields like quantum computing and nanotechnology.

Cambridge University Press eBooks•Jan 31, 2024•0 citations

Quantum Dynamics

Uri Peskin

Explicitly time-dependent Hamiltonians are ubiquitous in applications of quantum theory. It is therefore necessary to solve the time-dependent Schrödinger equation directly. The system’s dynamics is associated with a unitary time-evolution operator (a propagator), formally given as an infinite Dyson series. Time-dependent observables are invariant under unitary time-dependent transformations, where it is sometimes useful to transform the time-evolution from the states into the corresponding operators. This is carried out in part (in full) by transforming to the interaction (Heisenberg) picture. The corresponding equations of motion for the time-dependent operators are introduced. For quadratic potential energy functions, the time evolution of quantum expectation values coincides with the corresponding classical dynamics. This is demonstrated and analyzed in detail for Gaussian wave packets and a coherent state. Finally, we derive exact and approximate expressions for time-dependent transition probabilities and transition rates between quantum states. The validity of time-dependent perturbation theory is analyzed by comparison to exact dynamics.

Cambridge University Press eBooks•May 11, 2023•0 citations

Finite-difference approach to solving Heisenberg’s operator equations of motion: Application to one-dimensional time dependent Hamiltonians

Omar A. Sharafeddin

We reviewed and expanded on the finite difference approach introduced by Moncrief [Phys. Rev. D 28, 2485 (1983)] and Bender [Phys. Rev. Lett. 55, 901 (1985)] for solving Heisenberg’s operator equations of motion. In this approach, finite-difference recurrence relations are used to evolve the matrix representation of the operators in time. The advantages and disadvantages of this approach for the study of quantum processes in real time are discussed. The approach performed very well as illustrated by examples of harmonic and Morse oscillators interacting with continuous wave and pulsed laser fields.

The Journal of Chemical Physics•Jul 8, 1995•2 citations

Solution of Time-Independent Schrodinger Equation (TISE) by Using Finite Difference Approach

I Gusti Agung Widagda, Nengah Artawan, I Gde Antha Kasmawan +2 more

This research investigates the numeric solution of the time-independent Schrödinger equation for the quantum harmonic oscillator by finite difference approach. The harmonic oscillator, described by a quadratic function potential, is a fundamental model in quantum mechanics due to its broad applications, ranging from molecular vibrations to quantum field theory. The time-independent Schrödinger equation is a second-order differential equation that typically poses challenges when solved analytically for complex potentials. The finite difference method become an attractive choice as it transforms the continuous differential equation into a system of linear equations that can be computationally solved through computer programming code. In this study, the spatial domain is discretized, and the second derivative is calculated by using central differences, transforming the TISE into a tridiagonal matrix representing Hamiltonian of system. By finding solutions to this matrix eigenvalue problem, wavefunctions and eigenvalues are obtained. The study results demonstrate that the finite difference approach effectively solves the TISE for the harmonic oscillator. The results obtained by using the finite difference method closely approximate the analytical results. The linear regression results show respectively that the gradient (β1), regression coefficient (β0) and coefficient of determination (R²) approach ideal values of 1, 0, and 1. The z-test results also show that the value of calculated z < critical z, indicating that the wavefunction and probability density, whether estimated by using finite difference approach or analytical methods, are equivalent with confidence level of 95 percent.

BULETIN FISIKA•Feb 28, 2025•0 citations

EXPONENTIAL PROPAGATORS (INTEGRATORS) FOR THE TIME-DEPENDENT SCHRÖDINGER EQUATION

André D. Bandrauk, Huizhong Lu

The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.

Journal of Theoretical and Computational Chemistry•Apr 23, 2013•70 citations

Quantum Dynamics Using The Time-Dependent Schrödinger Equation

Abraham Nitzan

This chapter focuses on the time-dependent Schrödinger equation and its solutions for several prototype systems. It provides the basis for discussing and understanding quantum dynamics in condensed phases, however, a full picture can be obtained only by including also dynamical processes that destroy the quantum mechanical phase. Such a full description of quantum dynamics cannot be handled by the Schrödinger equation alone; a more general approach based on the quantum Liouville equation is needed. This important part of the theory of quantum dynamics is discussed in Chapter 10. Given a system characterized by a Hamiltonian Ĥ , the time-dependent Schrödinger equation is For a closed, isolated system Ĥ is time independent; time dependence in the Hamiltonian enters via effect of time-dependent external forces. Here we focus on the earlier case. Equation (1) is a first-order linear differential equation that can be solved as an initial value problem.

Oxford University Press eBooks•Apr 6, 2006•0 citations
We present an algorithm for solving the time-dependent Schrödinger equation that is based on the finite-difference expression of the kinetic energy operator. Students who have some knowledge of linear algebra can understand the theory used to derive the algorithm. This is because the finite-difference kinetic energy matrix and the Hückel matrix for linear conjugated hydrocarbons have similar forms. We offer 3 modules that can be used to solve 1-dimensional time-dependent Schrödinger equations, in real time. Module 1 contains two worksheets designed to show quantum dynamics in bound potentials. Module 2 contains 2 worksheets designed to illustrate tunneling. Module 3 contains 2 worksheets set up to illustrate the effect of time-dependent external fields. Further, modules 1 and 3 also include the Verlet algorithm so students can compare the quantum dynamical evolution of the systems to their classical counterpart in real time.
This project looks into the time evolution of a wave function within a two-dimensional quantum well. We start by solving the time-dependent Schrödinger equation for stationary states in a quantum well. Next, we express any wave function as a linear combination of stationary states, allowing us to understand their time evolution. Two methods are presented: one relies on decomposing the wave function into a basis of stationary states and the other on discretisation of the time-dependent Schrödinger equation, incorporating three-point formulas for derivatives. These approaches necessitate confronting intricate boundary conditions and require maintaining energy conservation for numerical accuracy. We further demonstrate the methods using a wave packet, revealing fundamental phenomena in quantum physics. Our results demonstrate the utility of these methods in understanding quantum systems, despite the challenges in determining stationary states for a given potential. This study enhances our comprehension of the dynamics of quantum states in constrained systems, essential for fields like quantum computing and nanotechnology.
Explicitly time-dependent Hamiltonians are ubiquitous in applications of quantum theory. It is therefore necessary to solve the time-dependent Schrödinger equation directly. The system’s dynamics is associated with a unitary time-evolution operator (a propagator), formally given as an infinite Dyson series. Time-dependent observables are invariant under unitary time-dependent transformations, where it is sometimes useful to transform the time-evolution from the states into the corresponding operators. This is carried out in part (in full) by transforming to the interaction (Heisenberg) picture. The corresponding equations of motion for the time-dependent operators are introduced. For quadratic potential energy functions, the time evolution of quantum expectation values coincides with the corresponding classical dynamics. This is demonstrated and analyzed in detail for Gaussian wave packets and a coherent state. Finally, we derive exact and approximate expressions for time-dependent transition probabilities and transition rates between quantum states. The validity of time-dependent perturbation theory is analyzed by comparison to exact dynamics.
We reviewed and expanded on the finite difference approach introduced by Moncrief [Phys. Rev. D 28, 2485 (1983)] and Bender [Phys. Rev. Lett. 55, 901 (1985)] for solving Heisenberg’s operator equations of motion. In this approach, finite-difference recurrence relations are used to evolve the matrix representation of the operators in time. The advantages and disadvantages of this approach for the study of quantum processes in real time are discussed. The approach performed very well as illustrated by examples of harmonic and Morse oscillators interacting with continuous wave and pulsed laser fields.
This research investigates the numeric solution of the time-independent Schrödinger equation for the quantum harmonic oscillator by finite difference approach. The harmonic oscillator, described by a quadratic function potential, is a fundamental model in quantum mechanics due to its broad applications, ranging from molecular vibrations to quantum field theory. The time-independent Schrödinger equation is a second-order differential equation that typically poses challenges when solved analytically for complex potentials. The finite difference method become an attractive choice as it transforms the continuous differential equation into a system of linear equations that can be computationally solved through computer programming code. In this study, the spatial domain is discretized, and the second derivative is calculated by using central differences, transforming the TISE into a tridiagonal matrix representing Hamiltonian of system. By finding solutions to this matrix eigenvalue problem, wavefunctions and eigenvalues are obtained. The study results demonstrate that the finite difference approach effectively solves the TISE for the harmonic oscillator. The results obtained by using the finite difference method closely approximate the analytical results. The linear regression results show respectively that the gradient (β1), regression coefficient (β0) and coefficient of determination (R²) approach ideal values of 1, 0, and 1. The z-test results also show that the value of calculated z < critical z, indicating that the wavefunction and probability density, whether estimated by using finite difference approach or analytical methods, are equivalent with confidence level of 95 percent.
The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.
This chapter focuses on the time-dependent Schrödinger equation and its solutions for several prototype systems. It provides the basis for discussing and understanding quantum dynamics in condensed phases, however, a full picture can be obtained only by including also dynamical processes that destroy the quantum mechanical phase. Such a full description of quantum dynamics cannot be handled by the Schrödinger equation alone; a more general approach based on the quantum Liouville equation is needed. This important part of the theory of quantum dynamics is discussed in Chapter 10. Given a system characterized by a Hamiltonian Ĥ , the time-dependent Schrödinger equation is For a closed, isolated system Ĥ is time independent; time dependence in the Hamiltonian enters via effect of time-dependent external forces. Here we focus on the earlier case. Equation (1) is a first-order linear differential equation that can be solved as an initial value problem.