© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

Journal of Physics Conference Series

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

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

Published Nov 1, 2024• 0 citations

Abstract

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.

DOI: 10.1088/1742-6596/2900/1/012009

Related papers

Numerical Simulation of the Time-Dependent Schrodinger Equation Using the Crank-Nicolson Method

Kabir, Adib

This study presents a numerical simulation of a quantum electron confined in a 10 nm potential well, using the Crank-Nicolson numerical technique to solve the time-dependent Schrodinger equation. The results capture the evolution of the electron's wave function at the 2000th time step, illustrating distinct standing wave patterns and probability densities that align with quantum mechanical predictions. Additionally, both 2D and 3D simulations across multiple time steps reveal the dynamic nature of quantum superposition and interference within the well. These findings highlight the method's stability and accuracy, offering a valuable tool for exploring quantum phenomena in constrained quantum systems.

arXiv:2410.10060•Oct 14, 2024•0 citations

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

Centered-Difference Applications for Schrödinger's Equation

Matthew Thomas Murachver

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.

•Jan 1, 2014•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

Numerical solution of the time-dependent Schrödinger equation in spherical coordinates by Fourier-transform methods

Christopher E. Dateo, Horia I. Metiu

We present a numerical Fourier-transform method for solving the time-dependent Schrödinger equation in spherical coordinates. The time propagation is accomplished with a novel implementation of a Crank–Nicholson scheme combined with a Feit–Fleck split-operator method. Calculations on model problems give results that are in excellent agreement with those obtained by other methods. The procedure avoids the numerical instabilities associated with the division by the sine of the polar angle. The propagator is unitary, does not require the calculation of any matrix elements, and uses no matrix multiplications other than fast Fourier transforms.

The Journal of Chemical Physics•Nov 15, 1991•44 citations

Numerical Solution of the Time Dependent 3D Schrödinger Equation Describing Tunneling of Atoms from Anharmonic Traps

I. S. Ishmukhamedov, Altay Ishmukhamedov, Vladimir S. Melezhik

We present an efficient numerical method for the integration of the 3D Schrödinger equation. A tunneling problem of two interacting bosonic atoms confined in a 1D anharmonic trap has been successfully solved by means of this method. We demonstrate fast convergence of the final results with respect to spatial and temporal grid steps. The computational scheme is based on the operator-splitting technique with the implicit Crank-Nicolson algorithm on spatial sixth-order finite-differences. The computational time is proportional to the number of spatial grid points.

EPJ Web of Conferences•Jan 1, 2018•3 citations

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.

This study presents a numerical simulation of a quantum electron confined in a 10 nm potential well, using the Crank-Nicolson numerical technique to solve the time-dependent Schrodinger equation. The results capture the evolution of the electron's wave function at the 2000th time step, illustrating distinct standing wave patterns and probability densities that align with quantum mechanical predictions. Additionally, both 2D and 3D simulations across multiple time steps reveal the dynamic nature of quantum superposition and interference within the well. These findings highlight the method's stability and accuracy, offering a valuable tool for exploring quantum phenomena in constrained quantum systems.

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.

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.

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.

We present a numerical Fourier-transform method for solving the time-dependent Schrödinger equation in spherical coordinates. The time propagation is accomplished with a novel implementation of a Crank–Nicholson scheme combined with a Feit–Fleck split-operator method. Calculations on model problems give results that are in excellent agreement with those obtained by other methods. The procedure avoids the numerical instabilities associated with the division by the sine of the polar angle. The propagator is unitary, does not require the calculation of any matrix elements, and uses no matrix multiplications other than fast Fourier transforms.

We present an efficient numerical method for the integration of the 3D Schrödinger equation. A tunneling problem of two interacting bosonic atoms confined in a 1D anharmonic trap has been successfully solved by means of this method. We demonstrate fast convergence of the final results with respect to spatial and temporal grid steps. The computational scheme is based on the operator-splitting technique with the implicit Crank-Nicolson algorithm on spatial sixth-order finite-differences. The computational time is proportional to the number of spatial grid points.