The programs described in this article and distributed with it aim (1) at integrating the optical Bloch equations governing the time evolution of the density matrix representing the quantum state of an atomic system driven by laser or microwave fields, and (2) at integrating the 1D Maxwell-Bloch equations for one or two laser fields co-propagating in an atomic vapour. The rotating wave approximation is assumed. These programs can also be used for more general quantum dynamical systems governed by the Lindblad master equation . They are written in Fortran 90; however, their use does not require any knowledge of Fortran programming. Methods for solving the optical Bloch equations in the rate equations limit, for calculating the steady-state density matrix and for formulating the optical Bloch equations in the weak probe approximation are also described. Program summary Program Title: CoOMBE CPC Library link to program files: https://doi.org/10.17632/5wsg9d52dk.1 Developers' repository link: https://github.com/durham-qlm/CoOMBE Licensing provisions: GPLv3 Programming language: Fortran 90 Nature of problem: The present programs can be used for the following operations: (1) Integrating the optical-Bloch equations within the rotating wave approximation for a multi-state atomic system. At the choice of the user, the calculation will return either the time-dependent density matrix at given times or the density matrix in the long time limit if the system evolves into a steady state in that limit. The calculation can be done with or without averaging over the thermal velocity distribution of the atoms. The number of atomic states which can be included in the calculation is limited only by the CPU time available and possibly by memory requirements. An arbitrarily large number of laser or microwave fields can be included in the calculation if these fields are all CW. This number is currently limited to one or two for fields that are not all CW. The calculation can be done in the weak probe approximation, or in the rate equations approximation, or without assuming either of these two approximations. Calculating refractive indexes , absorption coefficients and complex susceptibilities is also possible. (2) Integrating the 1D Maxwell-Bloch equations in the slowly varying envelope approximation for one or two fields co-propagating in a single-species atomic vapour. Although geared towards the case of atoms interacting with laser fields , this code can also be used for more general quantum systems with similar equations of motion (e.g., molecular systems, spin systems, etc.). Solution method: The Lindblad master equation is expressed as a system of homogeneous first order linear differential equations , which are transformed as required and solved to obtain the density matrix representing the state of the atomic system. A variety of methods are offered to this end. The same approach is also used in the calculation of the polarisation of the medium when integrating the Maxwell-Bloch equations. The latter are integrated over space using predictor-corrector methods. The library includes a general driving program making it possible to use these codes without additional program development. The distribution also includes examples of the use of a container for running these programs without a pre-installed Fortran compiler.
The programs described in this article and distributed with it aim (1) at integrating the optical Bloch equations governing the time evolution of the density matrix representing the quantum state of an atomic system driven by laser or microwave fields, and (2) at integrating the 1D Maxwell-Bloch equations for one or two laser fields co-propagating in an atomic vapour. The rotating wave approximation is assumed. These programs can also be used for more general quantum dynamical systems governed by the Lindblad master equation . They are written in Fortran 90; however, their use does not require any knowledge of Fortran programming. Methods for solving the optical Bloch equations in the rate equations limit, for calculating the steady-state density matrix and for formulating the optical Bloch equations in the weak probe approximation are also described. Program summary Program Title: CoOMBE CPC Library link to program files: https://doi.org/10.17632/5wsg9d52dk.1 Developers' repository link: https://github.com/durham-qlm/CoOMBE Licensing provisions: GPLv3 Programming language: Fortran 90 Nature of problem: The present programs can be used for the following operations: (1) Integrating the optical-Bloch equations within the rotating wave approximation for a multi-state atomic system. At the choice of the user, the calculation will return either the time-dependent density matrix at given times or the density matrix in the long time limit if the system evolves into a steady state in that limit. The calculation can be done with or without averaging over the thermal velocity distribution of the atoms. The number of atomic states which can be included in the calculation is limited only by the CPU time available and possibly by memory requirements. An arbitrarily large number of laser or microwave fields can be included in the calculation if these fields are all CW. This number is currently limited to one or two for fields that are not all CW. The calculation can be done in the weak probe approximation, or in the rate equations approximation, or without assuming either of these two approximations. Calculating refractive indexes , absorption coefficients and complex susceptibilities is also possible. (2) Integrating the 1D Maxwell-Bloch equations in the slowly varying envelope approximation for one or two fields co-propagating in a single-species atomic vapour. Although geared towards the case of atoms interacting with laser fields , this code can also be used for more general quantum systems with similar equations of motion (e.g., molecular systems, spin systems, etc.). Solution method: The Lindblad master equation is expressed as a system of homogeneous first order linear differential equations , which are transformed as required and solved to obtain the density matrix representing the state of the atomic system. A variety of methods are offered to this end. The same approach is also used in the calculation of the polarisation of the medium when integrating the Maxwell-Bloch equations. The latter are integrated over space using predictor-corrector methods. The library includes a general driving program making it possible to use these codes without additional program development. The distribution also includes examples of the use of a container for running these programs without a pre-installed Fortran compiler.
We present a general-purpose, open-source software suite, SQUIRREL (Streamlined Quantum Unified Interface for Researching Real-time Excitations with Light), for propagating the time-dependent Schrödinger equation on complex geometries in the presence of time-dependent electric and/or magnetic fields. To handle large systems that can be executed on a conventional desktop computer, the SQUIRREL software suite uses a suite of efficient propagation methods for various quantum dynamics applications, including a new perturbation-based element-dropping algorithm that improves computational performance with minimal loss of accuracy. We analyze the efficacy of these optimizations for Crank-Nicolson, scaled Taylor series approximation, and split-operator propagation methods and discuss the range of their applicability to a variety of quantum dynamics problems. In addition, we provide several examples of time-dependent dynamics calculations and extensive documentation for generating custom geometries, potentials, and time-propagation approaches. Our numerical benchmarks and results demonstrate the versatility of the SQUIRREL software suite for efficiently calculating quantum dynamics in complex nanoscale geometries, particularly in the presence of time-dependent magnetic fields, which have received less attention in previous quantum dynamics studies. Program summary Program Title: SQUIRREL CPC Library link to program files: https://doi.org/10.17632/kfvs5s88sj.1 Licensing provisions: GNU General Public License 3 Programming language: MATLAB Supplementary material: Additional details on calculations with an anisotropic effective mass, determination of timesteps, propagation timings, and element-dropping timings. Nature of problem: The SQUIRREL software suite solves the time-dependent Schrödinger equation for electronic states in the presence of time-dependent electric and/or magnetic fields. The code is well-suited for calculating electron dynamics of complex nanostructures in the effective mass approximation, using various propagation and sparse-matrix techniques to reduce the computational complexity of these problems. The program provides a user-friendly interface for testing these optimizations and generating custom geometries and potentials for systems of interest. Solution method: The SQUIRREL software suite uses Crank-Nicolson, split-operator, Padé approximant, and scaled-Taylor-series-based propagation methods to calculate electron dynamics on a finite-element basis. The code includes a novel perturbation-based element-dropping method, which enables sparse representations of the typically dense finite element Hamiltonian matrix. This capability enables a highly efficient and flexible approach for calculating driven dynamics in quantum systems with complex geometries and time-dependent electric/magnetic fields.
We present QDT (“quantum dynamics toolbox”), an open-source Matlab software package that enables users to simulate coupled quantum systems in the subsystem energy eigenbasis using modular functions. QDT requires no user knowledge of operator matrix assembly and automatically performs all necessary operator constructions and Hilbert space expansions. Density matrix propagation is performed by numerically solving the Liouville–von-Neumann equation. In order to simulate dissipation and decoherence effects, the Lindblad formalism is implemented. Furthermore, QDT supplies practical analysis and plotting functions, such as visualization of density matrix and expectation value dynamics, that facilitate the evaluation of simulation results. QDT further provides a module for the simulation of coherent multidimensional spectroscopy. Program Title: QDT — A Quantum Dynamics Toolbox CPC Library link to program files: https://doi.org/10.17632/786vb4rgy3.1 Code Ocean capsule: https://codeocean.com/capsule/4608215 Licensing provisions: MIT Programming language: Matlab R2018a (and above) with parallel computing toolbox (optional) [1] Supplementary material: Available on the QDT website [2] Nature of problem: Coupled quantum systems can be easily set up and discussed as level schemes in the subsystem energy eigenbasis. The simulation of the time evolution of these systems by numerically solving differential equations is not straightforward, however. Simulation code that is written from scratch is time-consuming to create and prone to errors during implementation. Furthermore, if code bases deviate drastically between research groups, reproducibility of results becomes challenging, and communication is hindered. Solution method: QDT addresses the above issues by providing a modular toolbox that contains, among others, elements for the simulation of multi-level quantum systems, quantum harmonic oscillators, their interaction with external electric fields, and dissipation. These elements can be combined by the user to build a specific system of interest. All operators are created and expanded automatically based on a single script file provided by the user. Simulation tasks are automated by a single function call, and parallelized in case of extensive tasks.
A generalized Bloch-Maxwell formulation for laser-field-coupled semiconductors is derived from a two-band model which includes direct Coulomb interactions. The momentum-dependent, microscopic, electron-hole equations of motion in the time-dependent Hartree-Fock approximation and neglecting interband exchange interactions form the starting point for the formulation. A self-consistent set of coupled equations in four dynamical variables for the medium, together with the electric-field amplitude coupled through the Maxwell wave equation in a semiclassical approximation, are obtained to lowest order in the coherent Coulomb exchange interaction and the density-of-states distribution. Intrinsic optical bistability is predicted in a steady state due to a carrier density-dependent redshift of the band edge due explicitly to coherent Coulomb exchange. Integration of the dynamical equations for conditions which correspond to the ultrafast time regime exhibit intrinsic adiabatic inversion, adiabatic following, anomalous Rabi cycling, and unique, as well as fast, optical switching, all of which depend upon the coherent Coulomb exchange interaction.
We introduce a non-Hermitian approximation of Bloch optical equations. This approximation provides a complete description of the excitation, relaxation, and decoherence dynamics of ensembles of coupled quantum systems in weak laser fields, taking into account collective effects and dephasing. In the proposed method, one propagates the wave function of the system instead of a complete density matrix. Relaxation and dephasing are taken into account via automatically adjusted time-dependent gain and decay rates. As an application, we compute the numerical wave packet solution of a time-dependent non-Hermitian Schrödinger equation describing the interaction of electromagnetic radiation with a quantum nano-structure, and compare the calculated transmission, reflection, and absorption spectra with those obtained from the numerical solution of the Liouville-von Neumann equation. It is shown that the proposed wave packet scheme is significantly faster than the propagation of the full density matrix while maintaining small error. We provide the key ingredients for easy-to-use implementation of the proposed scheme and identify the limits and error scaling of this approximation.
In this paper we report on the study of two and multi-level atoms interacting with multiple laser beams. The semi-classical approach is used to describe the system in which the atoms are treated quantum mechanically via the density matrix operator, while the laser beams are treated classically using Maxwells equations. We present results of a two level atom interacting with single and multiple laser beams and demonstrate Rabi oscillations between the levels. The effects of laser modulation on the dynamics of the atom (atomic populations and coherences) are examined by solving the optical Bloch equations. Plots of the density matrix elements as a function of time are presented for various parameters such as laser intensity, detuning, modulation etc. In addition, phase-space plots and Fourier analysis of the density matrix elements are provided. The atomic polarization, estimated from the coherence terms of the density matrix elements, is used in the numerical solution of Maxwells equations to determine the behaviour of the laser beams as they propagate through the atomic ensemble. The effects of saturation and hole-burning are demonstrated in the case of two counter propagating beams with one being a strong beam and the other being very weak. The above work is extended to include four-wave mixing in four level atoms in a diamond configuration. Two co-propagating beams of different wavelengths drive the atoms from a ground state |1〉 to an excited state |3〉 via an intermediate state |2〉. The atoms then move back to the ground state via another intermediate state |4〉, resulting in the generation of two additional correlated photon beams. The characteristics of these additional photons are studied.
wavepacket is a program to simulate the dynamics of a wave packet interacting with a time-dependent potential, with the time-dependent Schrödinger equation being solved on a one-, two-, or three-dimensional spatial grid using the split operator method. This new version fixes bugs present in the original program. Program Title: wavepacket CPC Library link to program files: https://doi.org/10.17632/gnxmkygxrw.1 Licensing provisions: GPLv3 Programming language: C (iso C99) Journal reference of previous version: Comput. Phys. Commun. 185 (2014) 407–414. Does the new version supersede the previous version?: Yes Reasons for the new version: Bug fixes. Summary of revisions: This new version fixes bugs in the program wavepacket [1]. In the original version, the function print_wf_bin, used (optionally) to write out the wave function as a binary file, contained a bug that would result in only part of the full wave function being written out when the program was run in parallel. Additionally, negative values of the parameters nt and/or nprint could lead to the program being stuck in an infinite loop. Additional checks have been added to ensure that both these parameters are positive (or zero). Nature of problem: Solves the time-dependent Schrödinger equation for a single particle interacting with a time-dependent potential. Solution method: The wave function is described by its value on a spatial grid and the evolution operator is approximated using the split-operator method [2,3], with the kinetic energy operator calculated using a Fast Fourier Transform. Additional comments including restrictions and unusual features: Simulation can be in one, two, or three dimensions. The fftw [4] implementation of the fast Fourier transform is used, so the corresponding library must be installed. Serial and parallel versions are compiled from the same source files. The parallel version is based on the Message-Passing Interface (mpi) [5]. C.M. Dion, A. Hashemloo, G. Rahali, Comput. Phys. Commun. 185 (2014) 407–414. M.D. Feit, J.A. Fleck, Jr., A. Steiger, Solution of the Schrödinger equation by a spectral method, J. Comput. Phys. 47 (1982) 412–433. M.D. Feit, J.A. Fleck, Jr., Solution of the Schrödinger equation by a spectral method II: Vibrational energy levels of triatomic molecules, J. Chem. Phys. 78 (1) (1983) 301–308. http://www.fftw.org http://www.mpi-forum.org