Quantum dynamics of open systems is usually formulated in terms of reduced density operator (diagonal matrix elements referred to as "populations", off-diagonal ones as "coherences"). Its evolution is governed by a quantum master equation i. e. a differential equation for the entire density matrix, both the "densities" (diagonal elements) and the "coherences" (off-diagonal elements), the latter ones being intrinsically quantum mechanical. Note that, due to a formal similarity to the Liouville equation in classical mechanics, the quantum master equation is also known as (quantum) Liouville-von Neumann equation or simply LvNE
The first term of the r.h.s. represents the closed system dynamics in terms of a Liouvillian superoperator for the Hermitian (aka Hamiltonian) part (using atomic units throughout)
with the Hamiltonian introduced on our Wiki page about Schrödinger's equations
for a quantum system interacting (within the semiclassical dipole approximation) with electrical field F(t) and where μ is a Hermitian matrix with dipole moments. In our Matlab codes we make use of an energy representation where the Hamiltonian of the unperturbed system, H_{0}=T+V, is diagonal.
In order to retrieve the expectation value of a quantum-mechanical operator/observable O from the solution of the LvNE, the following trace formula is used
where the time-dependence can arise from time-dependent O or time-dependent ρ or both.
The second term on the r.h.s. of the first equation represents the coupling to the environment accounting for time-directed irreversibility, i. e., dissipation and/or dephasing. The most widely used formalisms in use are the following ones:
In conclusion, the above Ansatz of the quantum-master (Liouville-von Neumann) equation in Born-Markov approximation (weak coupling, no memory) provides a description useful in many cases ranging from NMR over quantum optics to molecular/chemical applications.
In order to bring the above evolution equations into a shape compatible with the standard form of linear/bilinear control systems and also with standard ODE solvers, the density operators (or, rather their matrix representations) have to be mapped onto vectors and, accordingly, the superoperators have to be mapped onto matrices.
Frequently used choices of mapping density matrices onto vectors are the following:
Column-wise: this mimics the memory organization in Fortran and MATLAB, in particular, the (:) operator of the latter. This is equivalent to the mathematical concept of vectorization which converts a matrix into a column vector by stacking the columns of the matrix A on top of one another.
Row-wise: this mimics the memory organization in C and in C++.
Diagonal first: First the populations, then the coherences (where the order of the latter ones is arbitrary but fixed), which leads to block structures of the Lindblad Liouvillian superoperators, see e. g. Appendix A of [1] and/or our Wiki page on linear/bilinear control systems.
Wiki: Numerics.Control.Lindblad
Wiki: Numerics.Control
Wiki: Numerics.LvNE.Lindblad
Wiki: Numerics.LvNE.Redfield
Wiki: Numerics.Main