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

Centered-Difference Applications for Schrödinger's Equation

Matthew Thomas Murachver

Published Jan 1, 2014• 0 citations
OpenAlex

Abstract

This project enumerates methods utilizing discretized centered-difference approximations on the second order differential equation for quantum particles known as Schrodinger’s Equation. An eigenvalue-eigenfunction scheme is developed to sieve for valid solutions to The Time Independent Schrodinger Equation. Additionally the Crank-Nicolson method is applied to the Time Dependent Schrodinger Equation to describe wavefunction (eigenfunction) time evolution. The validity of these methods is discussed with applications to several fundamental pedagogical introductory quantum mechanic systems.

Related papers

Modified Crank-Nicolson method for solving the time-dependent Schrödinger equation

A R Kurniawan, Evi Suaebah, Z. I. Supardi

Abstract The time-dependent Schrödinger equation is fundamental to quantum mechanics, describing the evolution of quantum systems. This study introduces a novel numerical approach that extends and modifies the Crank-Nicholson method to solve the time-dependent Schrödinger equation. This method is designed to enhance both the precision and performance of computational results, especially in the context of quantum tunneling phenomena and applications that require high accuracy over long periods. To assess the precision and performance of the modified Crank-Nicolson method, we benchmark it against two well-known numerical methods: the Runge-Kutta method and the Split-Step Fourier method. Quantum tunneling is selected as the case study due to its relevance in quantum physics applications. Comparative analysis is based on numerical error metrics. The proposed method can be applied to other physical phenomena, expanding its potential application in physics research and technology industries.

Journal of Physics Conference Series•Nov 1, 2024•0 citations

Accurate numerical solutions of the time-dependent Schrödinger equation

W. van Dijk, F. Masafumi Toyama

We present a generalization of the often-used Crank-Nicolson (CN) method of obtaining numerical solutions of the time-dependent Schrödinger equation. The generalization yields numerical solutions accurate to order (Deltax)2r-1 in space and (Deltat)2M in time for any positive integers r and M, while CN employ r=M=1. We note dramatic improvement in the attainable precision (circa ten or greater orders of magnitude) along with several orders of magnitude reduction of computational time. The improved method is shown to lead to feasible studies of coherent-state oscillations with additional short-range interactions, wave-packet scattering, and long-time studies of decaying systems.

arXiv:physics/0701150•Physical Review E•Mar 23, 2007•107 citations

Time—Independent Perturbation Theory

Şaki̇r Erkoç

There are few potentials V (r) for which the Schrödinger equation is exactly solvable. We must therefore develop approximation techniques to obtain the eigenvalues and eigenfunctions for potentials that do not lead to exactly soluble equations. One of these methods is the perturbation theory, which was developed both for time independent and time dependent phenomena. We will discuss the time independent perturbation theory. The formulation of the time independent perturbation theory is different for nondegenerate and degenerate states.

CRC Press eBooks•Jul 28, 2006•0 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

The direct midpoint method as a quantum mechanical integrator III

Ulrich Mutze

This talk is about a method for solving the numerical initial value problem of the time-dependent Schrodinger equation with a Hamiltonian which may depend on time. It is only for time independent Hamiltonians that I was able to describe the accuracy and stability properties in the form of proven theorems. All practical experience shows, that the time-dependent case behaves equally well. In this method, the Hamiltonian comes in only as a hermitian linear operator and not as a mathematical expression that can be reordered or otherwise transformed in order to speed up the necessary calculations. In all cases where it is a feasible computational task to apply the Hamilton operator to a generic state vector, my suggestion is to try this method as a workhorse. It can be considered a ‘Simpson’s rule’ for quantum dynamics. There is no way for me to explain the matter from a technical point of view in 10 minutes, so I’ll try to convey the basic idea of the method. This idea is to treat quantum dynamics in closer analogy to classical dynamics as one usually does. One may argue, that there is a fundamental difference between classical and quantum dynamics since the Schrodinger equation is first order with respect to the time derivative, whereas classical mechanics is inherently second order, since forces determine accelerations and not velocities. This is not so clear a distinction as one might think, since also classical mechanics can be given the first-order form of canonical equations. I propose the opposite of such an approach, namely to give quantum dynamics a second-order form by simply differentiating the time-dependent Schrodinger equation with respect to time. To be sure, this is only a heuristic device for defining a computationally feasible approximation and not a proposal to change quantum mechanics. For the resulting equation for the second time derivative of the wave function, I then carry over and use the integration methods from the field of granular matter dynamics which I have been using in an industrial project concerning simulation of the toning process in electro-photographic copiers/printers over some years [4]. Such methods come in

•Jan 1, 2007•1 citation

Numerical solutions of the Schrödinger equation with source terms or time-dependent potentials

W. van Dijk, F. Masafumi Toyama

We develop an approach to solving numerically the time-dependent Schrödinger equation when it includes source terms and time-dependent potentials. The approach is based on the generalized Crank-Nicolson method supplemented with an Euler-MacLaurin expansion for the time-integrated nonhomogeneous term. By comparing the numerical results with exact solutions of analytically solvable models, we find that the method leads to precision comparable to that of the generalized Crank-Nicolson method applied to homogeneous equations. Furthermore, the systematic increase in precision generally permits making estimates of the error.

arXiv:1412.1802•Physical Review E•Dec 17, 2014•12 citations
This project enumerates methods utilizing discretized centered-difference approximations on the second order differential equation for quantum particles known as Schrodinger’s Equation. An eigenvalue-eigenfunction scheme is developed to sieve for valid solutions to The Time Independent Schrodinger Equation. Additionally the Crank-Nicolson method is applied to the Time Dependent Schrodinger Equation to describe wavefunction (eigenfunction) time evolution. The validity of these methods is discussed with applications to several fundamental pedagogical introductory quantum mechanic systems.
Abstract The time-dependent Schrödinger equation is fundamental to quantum mechanics, describing the evolution of quantum systems. This study introduces a novel numerical approach that extends and modifies the Crank-Nicholson method to solve the time-dependent Schrödinger equation. This method is designed to enhance both the precision and performance of computational results, especially in the context of quantum tunneling phenomena and applications that require high accuracy over long periods. To assess the precision and performance of the modified Crank-Nicolson method, we benchmark it against two well-known numerical methods: the Runge-Kutta method and the Split-Step Fourier method. Quantum tunneling is selected as the case study due to its relevance in quantum physics applications. Comparative analysis is based on numerical error metrics. The proposed method can be applied to other physical phenomena, expanding its potential application in physics research and technology industries.
We present a generalization of the often-used Crank-Nicolson (CN) method of obtaining numerical solutions of the time-dependent Schrödinger equation. The generalization yields numerical solutions accurate to order (Deltax)2r-1 in space and (Deltat)2M in time for any positive integers r and M, while CN employ r=M=1. We note dramatic improvement in the attainable precision (circa ten or greater orders of magnitude) along with several orders of magnitude reduction of computational time. The improved method is shown to lead to feasible studies of coherent-state oscillations with additional short-range interactions, wave-packet scattering, and long-time studies of decaying systems.
There are few potentials V (r) for which the Schrödinger equation is exactly solvable. We must therefore develop approximation techniques to obtain the eigenvalues and eigenfunctions for potentials that do not lead to exactly soluble equations. One of these methods is the perturbation theory, which was developed both for time independent and time dependent phenomena. We will discuss the time independent perturbation theory. The formulation of the time independent perturbation theory is different for nondegenerate and degenerate states.
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.
This talk is about a method for solving the numerical initial value problem of the time-dependent Schrodinger equation with a Hamiltonian which may depend on time. It is only for time independent Hamiltonians that I was able to describe the accuracy and stability properties in the form of proven theorems. All practical experience shows, that the time-dependent case behaves equally well. In this method, the Hamiltonian comes in only as a hermitian linear operator and not as a mathematical expression that can be reordered or otherwise transformed in order to speed up the necessary calculations. In all cases where it is a feasible computational task to apply the Hamilton operator to a generic state vector, my suggestion is to try this method as a workhorse. It can be considered a ‘Simpson’s rule’ for quantum dynamics. There is no way for me to explain the matter from a technical point of view in 10 minutes, so I’ll try to convey the basic idea of the method. This idea is to treat quantum dynamics in closer analogy to classical dynamics as one usually does. One may argue, that there is a fundamental difference between classical and quantum dynamics since the Schrodinger equation is first order with respect to the time derivative, whereas classical mechanics is inherently second order, since forces determine accelerations and not velocities. This is not so clear a distinction as one might think, since also classical mechanics can be given the first-order form of canonical equations. I propose the opposite of such an approach, namely to give quantum dynamics a second-order form by simply differentiating the time-dependent Schrodinger equation with respect to time. To be sure, this is only a heuristic device for defining a computationally feasible approximation and not a proposal to change quantum mechanics. For the resulting equation for the second time derivative of the wave function, I then carry over and use the integration methods from the field of granular matter dynamics which I have been using in an industrial project concerning simulation of the toning process in electro-photographic copiers/printers over some years [4]. Such methods come in
We develop an approach to solving numerically the time-dependent Schrödinger equation when it includes source terms and time-dependent potentials. The approach is based on the generalized Crank-Nicolson method supplemented with an Euler-MacLaurin expansion for the time-integrated nonhomogeneous term. By comparing the numerical results with exact solutions of analytically solvable models, we find that the method leads to precision comparable to that of the generalized Crank-Nicolson method applied to homogeneous equations. Furthermore, the systematic increase in precision generally permits making estimates of the error.