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

arXiv:2506.19227

Numerical study of the radiation-matter interaction quantum systems through the time-dependent Schrödinger dynamics

Leonardi Hernández-Sánchez, I. Ramos-Prieto, Ariel Flores-Rosas, Francisco Soto, H. M. Moya-Cessa

Published Jun 24, 2025• 0 citations
arXivHTMLPDF

Abstract

Obtaining exact solutions to the Schrödinger equation in complex quantum systems poses significant challenges. In this context, numerical methods emerge as valuable tools for analyzing such systems. This article proposes a numerical approach using the fourth-order Runge-Kutta method, implemented in Python, to tackle radiation-matter interaction systems. This methodology is applicable to various Hamiltonians, including that of the Jaynes-Cummings model. The accuracy of the numerical results is validated by comparing them with analytical solutions in simplified cases, demonstrating its effectiveness in studying quantum systems where exact solutions cannot be derived.

DOI
10.48550/arxiv.2506.19227

Related papers

Numerical solution of the Schrödinger equation for radiation-matter interaction models using the Runge-Kutta method implemented in Python

Leonardi Hernández-Sánchez, Ariel Flores-Rosas, I. Ramos-Prieto +2 more

Obtaining exact solutions to the time-dependent Schrodinger equation in complex quantum systems presents significant challenges. In this ¨ context, numerical methods offer powerful alternatives for exploring their dynamics. In this pedagogical article, we present a numerical approach based on the fourth-order Runge-Kutta method, implemented in Python, to simulate radiation-matter interaction models. The methodology is illustrated using the well-known Jaynes-Cummings model, and the numerical results are compared with its exact analytical solution for illustrative purposes. Although only this model is explicitly solved, the numerical framework is general and can be readily extended to more complex Hamiltonians, including time-dependent cases. This makes the approach a practical and accessible tool for exploring a wide range of quantum systems.

Revista Mexicana de Física E•Jul 1, 2026•0 citations

Numerical solution of the Lindblad master equation using the Runge-Kutta method implemented in Python

Leonardi Hernández-Sánchez, Ivan A. Bocanegra-Garay, Ariel Flores-Rosas +3 more

The dynamics of open quantum systems is governed by the Lindblad master equation, which provides a consistent framework for incorporating environmental effects into the evolution of the system. Since exact solutions are rarely available, numerical methods become essential tools for analyzing such systems. This article presents a step-by-step implementation of the fourth-order Runge-Kutta method in Python to solve the Lindblad equation for a single quantized field mode subject to decay. A coherent state is used as the initial condition, and the time evolution of the average photon number is investigated. The proposed methodology enables transparent and customizable simulations of dissipative quantum dynamics, emphasizing a pedagogical approach that helps readers understand the numerical structure without relying on external libraries such as QuTiP. This standalone implementation offers full control over each integration step, making it particularly suitable for educational contexts and for exploring non-standard dynamics or introducing custom modifications to the Liouvillian.

arXiv:2506.19238•Revista Mexicana de Física E•Jan 1, 2026•2 citations

Numerical Case Study of an Atom-Photon Interaction in a Cavity Exploring Quantum Control

Javier Jalandoni

We study Magnus expansion (ME) approximation scheme for the interaction between an atom and a single quantized cavity mode (Jaynes-Cumming model) in a closed quantum system in resonance or near resonance for a time-dependent coupling coefficient g(t) in both the interaction and rotating picture by implementing a novel numerical method called MG4 and compare out results to the Runge-Kutta 4th (RK4) order solution to demonstrate the conservation of unitary evolution of the ME. A cursory study of open quantum system is given to encourage the study of ME for dissipative systems. Furthermore, we assume that our time-dependent coupling coefficient g(t) can take on two forms, Gaussian and sinusoidal, which are introduced as pulses to study the behavior and response of the cavity. Our results show that ME is a sufficient approximation scheme in our study of closed quantum systems which may have applications in quantum control.

•Jan 1, 2016•1 citation

Numerical Method for a Confined Atomic System

Clark, Deleonne, A. Sarsa, Ehme, Jeffry

Quantum confinement is a relatively new subject matter in quantum mechanics. Quantum confinement traps the atom in a cavity whose dimensions are small enough to alter its properties. Some physicists have used the time-independent Schrodinger equation to obtain exact results for confinement problems. But this paper focuses on the dynamic of a hydrogen atom when it is moved from captivity. The goal of this project is to study and understand different methods used to predict solutions that cannot be obtained analytically which was done by relying on methods that approximates solutions. To do this, a differential equation with boundary conditions was manipulated with the computer programs C++ and Eclipse using Runge-Kutta, B-Splines, and the Schrodinger equation. While the Runge-Kutta method was utilized as a basis for the boundary problem, its simplicity was not viable for this research which contained so many complex variables. Instead, B-Splines were used to solve the time-dependent Schrodinger equation. B-Spline functions are designed to generalize polynomials for the purpose of approximating arbitrary functions. By taking the summation of the timeindependent Schrodinger equation, we obtained our desired results for the timedependent Schrodinger equation.

•Jul 1, 2015•0 citations

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

A generalized unitary transformation for solving time-independent and time-dependent Jaynes–Cummings models under resonant and selected off-resonance conditions

ANTONINO MESSINA, Hiromichi Nakazato, Antonio S. M. de Castro +4 more

Jaynes–Cummings (JC) models are by far the most widely used to study the dynamics of two-level systems, such as atoms or molecules under suitable conditions, coupled to quantized electromagnetic fields. Nowadays, analytical solutions of JC models are more useful than ever for understanding and controlling radiation-matter interactions relevant to quantum information processing, quantum computing, and other quantum technologies. This study presents the analysis and generalization of a method for the solution of the dynamical problem in the time-independent JC model (Carbonaro et al., 1979) that enables an exact description of the atom-field dynamics in the resonant time-dependent JC model and in a nontrivial class of off-resonance time-dependent JC models. Our analysis produces the ready-to-use forms of the evolution operator in the resonant and non-resonant cases, which are then used to derive the evolution of the atomic population inversion for some types of atom-field coupling. The central idea behind the proposed approach offers opportunities for its extension to other time-dependent cases and its generalization to JC-related models.

Annals of Physics•Jul 7, 2026•0 citations
Obtaining exact solutions to the Schrödinger equation in complex quantum systems poses significant challenges. In this context, numerical methods emerge as valuable tools for analyzing such systems. This article proposes a numerical approach using the fourth-order Runge-Kutta method, implemented in Python, to tackle radiation-matter interaction systems. This methodology is applicable to various Hamiltonians, including that of the Jaynes-Cummings model. The accuracy of the numerical results is validated by comparing them with analytical solutions in simplified cases, demonstrating its effectiveness in studying quantum systems where exact solutions cannot be derived.
Obtaining exact solutions to the time-dependent Schrodinger equation in complex quantum systems presents significant challenges. In this ¨ context, numerical methods offer powerful alternatives for exploring their dynamics. In this pedagogical article, we present a numerical approach based on the fourth-order Runge-Kutta method, implemented in Python, to simulate radiation-matter interaction models. The methodology is illustrated using the well-known Jaynes-Cummings model, and the numerical results are compared with its exact analytical solution for illustrative purposes. Although only this model is explicitly solved, the numerical framework is general and can be readily extended to more complex Hamiltonians, including time-dependent cases. This makes the approach a practical and accessible tool for exploring a wide range of quantum systems.
The dynamics of open quantum systems is governed by the Lindblad master equation, which provides a consistent framework for incorporating environmental effects into the evolution of the system. Since exact solutions are rarely available, numerical methods become essential tools for analyzing such systems. This article presents a step-by-step implementation of the fourth-order Runge-Kutta method in Python to solve the Lindblad equation for a single quantized field mode subject to decay. A coherent state is used as the initial condition, and the time evolution of the average photon number is investigated. The proposed methodology enables transparent and customizable simulations of dissipative quantum dynamics, emphasizing a pedagogical approach that helps readers understand the numerical structure without relying on external libraries such as QuTiP. This standalone implementation offers full control over each integration step, making it particularly suitable for educational contexts and for exploring non-standard dynamics or introducing custom modifications to the Liouvillian.
We study Magnus expansion (ME) approximation scheme for the interaction between an atom and a single quantized cavity mode (Jaynes-Cumming model) in a closed quantum system in resonance or near resonance for a time-dependent coupling coefficient g(t) in both the interaction and rotating picture by implementing a novel numerical method called MG4 and compare out results to the Runge-Kutta 4th (RK4) order solution to demonstrate the conservation of unitary evolution of the ME. A cursory study of open quantum system is given to encourage the study of ME for dissipative systems. Furthermore, we assume that our time-dependent coupling coefficient g(t) can take on two forms, Gaussian and sinusoidal, which are introduced as pulses to study the behavior and response of the cavity. Our results show that ME is a sufficient approximation scheme in our study of closed quantum systems which may have applications in quantum control.
Quantum confinement is a relatively new subject matter in quantum mechanics. Quantum confinement traps the atom in a cavity whose dimensions are small enough to alter its properties. Some physicists have used the time-independent Schrodinger equation to obtain exact results for confinement problems. But this paper focuses on the dynamic of a hydrogen atom when it is moved from captivity. The goal of this project is to study and understand different methods used to predict solutions that cannot be obtained analytically which was done by relying on methods that approximates solutions. To do this, a differential equation with boundary conditions was manipulated with the computer programs C++ and Eclipse using Runge-Kutta, B-Splines, and the Schrodinger equation. While the Runge-Kutta method was utilized as a basis for the boundary problem, its simplicity was not viable for this research which contained so many complex variables. Instead, B-Splines were used to solve the time-dependent Schrodinger equation. B-Spline functions are designed to generalize polynomials for the purpose of approximating arbitrary functions. By taking the summation of the timeindependent Schrodinger equation, we obtained our desired results for the timedependent Schrodinger equation.
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.
Jaynes–Cummings (JC) models are by far the most widely used to study the dynamics of two-level systems, such as atoms or molecules under suitable conditions, coupled to quantized electromagnetic fields. Nowadays, analytical solutions of JC models are more useful than ever for understanding and controlling radiation-matter interactions relevant to quantum information processing, quantum computing, and other quantum technologies. This study presents the analysis and generalization of a method for the solution of the dynamical problem in the time-independent JC model (Carbonaro et al., 1979) that enables an exact description of the atom-field dynamics in the resonant time-dependent JC model and in a nontrivial class of off-resonance time-dependent JC models. Our analysis produces the ready-to-use forms of the evolution operator in the resonant and non-resonant cases, which are then used to derive the evolution of the atomic population inversion for some types of atom-field coupling. The central idea behind the proposed approach offers opportunities for its extension to other time-dependent cases and its generalization to JC-related models.