Numerical solution of the Schrödinger equation for radiation-matter interaction models using the Runge-Kutta method implemented in Python
DOI:
https://doi.org/10.31349/RevMexFisE.23.020202Keywords:
Schrödinger equation, radiation-matter interaction, Runge-Kutta method, Jaynes-Cummings model, PythonAbstract
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.
Downloads
References
C. Gerry and P. Knight, Introductory Quantum Optics (Cambridge University Press, 2004), https://doi.org/10.1017/CBO9780511791239
A. B. Klimov and S. M. Chumakov, A Group-Theoretical Approach to Quantum Optics (John Wiley & Sons, Ltd, 2009), https://doi.org/10.1002/9783527624003.ch5
G. Agarwal, Quantum Optics, Quantum Optics (Cambridge University Press, 2013), https://books.google.com.mx/books?id=7KKw_XIYaioC
E. Jaynes and F. Cummings, Comparison of quantum and semiclassical radiation theories with application to the beam maser, Proc. IEEE 51 (1963) 89, https://doi.org/10.1109/PROC.1963.1664
B. W. Shore and P. L. Knight, The Jaynes-Cummings Model, J. Mod. Opt. 40 (1993) 1195, https://doi.org/10.1080/09500349314551321
J. Larson and T. K. Mavrogordatos, The Jaynes-Cummings Model and Its Descendants (Institute of Physics Publishing, 2022), https://api.semanticscholar.org/CorpusID:245226335
V. Buzek, N-level Atom Interacting with Single-mode Radiation Field: An Exactly Solvable Model with Multiphoton Transitions and Intensity-dependent Coupling, J. Mod. Opt. 37 (1990) 1033, https://doi.org/10.1080/09500349014551091
M. Nath, S. Sen, and G. Gangopadhyay, Dynamics of cascade three-level system interacting with the classical and quantized field, Pramana-J. Phys. 61 (2003) 1089, https://doi.org/10.1007/BF02704404
C. Sukumar and B. Buck, Multi-phonon generalisation of the Jaynes-Cummings model, Phys. Lett. A 83 (1981) 211, https://doi.org/10.1016/0375-9601(81)90825-2
M. S. Abdalla, M. M. A. Ahmed, and A.-S. F. Obada, Multimode and multiphoton processes in a non-linear JaynesCummings model, Phys. A: Stat. Mech. Appl. 170 (1991) 393, https://doi.org/10.1016/0378-4371(91)90051-D
P. Alsing, D.-S. Guo, and H. J. Carmichael, Dynamic Stark effect for the Jaynes-Cummings system, Phys. Rev. A 45 (1992) 5135, https://doi.org/10.1103/PhysRevA.45.5135
I. A. Bocanegra-Garay et al., Invariant approach to the driven Jaynes-Cummings model, SciPost Phys. 16 (2024) 007, https://doi.org/10.21468/SciPostPhys.16.1.007
L. Hernández-Sánchez et al., Effects of Classical Drivings on the Power Broadening of Atomic Lineshapes, J. Opt. Soc. Am. B 41 (2024) C68,https://doi.org/10.1364/JOSAB.522587
L. Hernández-Sánchez et al., Discriminando Superposiciones de Estados Coherentes Mediante Formas de Línea, Rev. Mex. Fis. 70 (2024) 011302 1-8, https://doi.org/10.31349/RevMexFis.70.011302
L. Hernández-Sánchez et al., Effects of Squeezing on the Power Broadening and Shifts of Micromaser Lineshapes, Photonics 11 (2024), https://doi.org/10.3390/photonics11040371
R. Juarez Amaro, A. Zúñiga-Segundo, and H. Moya-Cessa, Several Ways to Solve the Jaynes-Cummings Model, App. Math. Inf. Sci. 9 (2015) 299, https://doi.org/10.12785/amis/090136
J. C. Butcher, Numerical Methods for Ordinary Differential Equations (John Wiley & Sons, Ltd, 2016), https://onlinelibrary.wiley.com/doi/book/10.1002/9781119121534
W. H. Press, B. P. Flannery, and S. A. Teukolsky, Numerical Recipes: The Art of Scientific Computing (Cambridge University Press, 1986), https://cfatab.harvard.edu/nr/
J. Johansson, P. Nation, and F. Nori, QuTiP 2: A Python framework for the dynamics of open quantum systems, Computer Physics Communications 184 (2013) 1234, https://doi.org/10.1016/j.cpc.2012.11.019
L. Hernández-Sánchez, I. A. Bocanegra-Garay, A. Flores-Rosas, I. Ramos-Prieto, F. Soto-Eguibar, and H. Moya-Cessa, Numerical solution of the Lindblad master equation using the Runge-Kutta method implemented in Python, Rev. Mex. Fis. E 23 (2026) 010214, https://doi.org/10.31349/RevMexFisE.23.010214
L. Hernández-Sánchez et al., Formas de línea atómicas en modelos de tipo Jaynes-Cummings (Editorial Académica Española, 2023), https://www.eae-publishing.com/catalogue/details/es/978-620-2-15886-2/formas-de-l%C3%ADnea-at%C3%B3micas
Downloads
Published
How to Cite
Issue
Section
License
Copyright (c) 2026 L. Hernández-Sánchez, A. Flores-Rosas, I. Ramos-Prieto, F. Soto-Eguibar, H. M. Moya-Cessa

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.
Authors retain copyright and grant the Revista Mexicana de Física E right of first publication with the work simultaneously licensed under a CC BY-NC-ND 4.0 that allows others to share the work with an acknowledgement of the work's authorship and initial publication in this journal.