Lindblad superoperator
Encyclopedia
The Lindblad superoperator is often used to express the quantum master equation for a dissipative system
.
In the canonical formulation of quantum mechanics, a system's time evolution is governed by unitary dynamics. This implies that phase coherence is maintained throughout the process, and is a consequence of the fact that all participating degrees of freedom are considered. However, any real physical system is not absolutely isolated, and will interact with its environment. This interaction with degrees of freedom external to the system results in dissipation of energy into the surroundings, and randomization of phase. This latter effect is the reason quantum mechanics is difficult to observe on a macroscopic scale. Moreso, understanding the interaction of a quantum system with its environment is necessary to understanding many commonly observed phenomena like the spontaneous emission of atoms, or the performance of many quantum technological devices, like the laser.
Certain mathematical techniques have been introduced to treat the interaction of a quantum system with its environment. One of these is the use of the density matrix
, and its associated master equation. While in principle this approach to solving quantum dynamics is equivalent to the Schrödinger picture
or Heisenberg picture
, it allows more easily for the inclusion of incoherent processes, which represent environmental interactions. The density operator has the property that it can represent a classical mixture of quantum states, and is thus vital to accurately describe the dynamics of so-called open quantum systems.
For a collapse operator
, the Lindblad superoperator, acting on the density matrix
, is
Such a term is found regularly in the Lindblad equation
as used in quantum optics
, where it can express absorption or emission of photons from a reservoir. For example, the master equation for a single mode optical resonator (e.g. a Fabry–Perot cavity) coupled to a thermal bath is
where
is the frequency of the optical mode, 
is the mode annihilation operator, 
is the mode linewidth, and 
is the thermal occupation number of the photons in the bath, as given by the Bose–Einstein distribution.
The dynamics of the entire system can be described by the Liouville equation of motion,
. This equation, containing an infinite number of degrees of freedom, is impossible to solve analytically except in very particular cases. What's more, under certain approximations, the bath degrees of freedom need not be considered, and an effective master equation can be derived in terms of the system density matrix, 
. The problem can be analyzed more easily by moving into the interaction picture, defined by the unitary transformation 
, where 
is an arbitrary operator, and 
. It is straightforward to confirm that the Liouville equation becomes
where the Hamiltonian
is explicitly time dependent. This equation can be integrated directly to give
This implicit equation for
can be substituted back into the Liouville equation to obtain an exact differo-integral equation
We proceed with the derivation by assuming the interaction is initiated at
, and at that time there are no correlations between the system and the bath. This implies that the initial condition is factorable as 
, where 
is the density operator of the bath initially. Since the system and bath Hamiltonians act on different Hilbert subspaces, they commute, and thus in the interaction picture we can write
tracing over the bath degrees of freedom,
. Tracing over the bath degrees of freedom of the aforementioned differo-integral equation yields
This equation is exact for the time dynamics of the system density matrix but requires full knowledge of the dynamics of the bath degrees of freedom. A simplifying assumption called the Born approximation rests on the largeness of the bath and the relative weakness of the coupling, which is to say the coupling of the system to the bath should significantly alter the bath eigenstates. In this case the full density matrix is factorable for all times as
. The master equation becomes
The equation is now explicit in the system degrees of freedom, but is very difficult to solve. A final assumption is the Born-Markoff approximation that the time derivative of the density matrix depends only on its current state, and not on its past. This assumption is valid under fast bath dynamics, wherein correlations between bath and system variables are lost extremely quickly, and amounts to replacing
on the right hand side of the equation.
This is the final form of the master equation we need.
for system operators
and bath operators 
, the master equation becomes
which can be expanded as
The expectation values
are with respect to the bath degrees of freedom.
Dissipative system
A dissipative system is a thermodynamically open system which is operating out of, and often far from, thermodynamic equilibrium in an environment with which it exchanges energy and matter....
.
In the canonical formulation of quantum mechanics, a system's time evolution is governed by unitary dynamics. This implies that phase coherence is maintained throughout the process, and is a consequence of the fact that all participating degrees of freedom are considered. However, any real physical system is not absolutely isolated, and will interact with its environment. This interaction with degrees of freedom external to the system results in dissipation of energy into the surroundings, and randomization of phase. This latter effect is the reason quantum mechanics is difficult to observe on a macroscopic scale. Moreso, understanding the interaction of a quantum system with its environment is necessary to understanding many commonly observed phenomena like the spontaneous emission of atoms, or the performance of many quantum technological devices, like the laser.
Certain mathematical techniques have been introduced to treat the interaction of a quantum system with its environment. One of these is the use of the density matrix
Density matrix
In quantum mechanics, a density matrix is a self-adjoint positive-semidefinite matrix of trace one, that describes the statistical state of a quantum system...
, and its associated master equation. While in principle this approach to solving quantum dynamics is equivalent to the Schrödinger picture
Schrödinger picture
In physics, the Schrödinger picture is a formulation of quantum mechanics in which the state vectors evolve in time, but the operators are constant. This differs from the Heisenberg picture which keeps the states constant while the observables evolve in time...
or Heisenberg picture
Heisenberg picture
In physics, the Heisenberg picture is a formulation of quantum mechanics in which the operators incorporate a dependency on time, but the state vectors are time-independent. It stands in contrast to the Schrödinger picture in which the operators are constant and the states evolve in time...
, it allows more easily for the inclusion of incoherent processes, which represent environmental interactions. The density operator has the property that it can represent a classical mixture of quantum states, and is thus vital to accurately describe the dynamics of so-called open quantum systems.
For a collapse operator
, the Lindblad superoperator, acting on the density matrixDensity matrix
In quantum mechanics, a density matrix is a self-adjoint positive-semidefinite matrix of trace one, that describes the statistical state of a quantum system...
, isSuch a term is found regularly in the Lindblad equation
Lindblad equation
In quantum mechanics, the Kossakowski–Lindblad equation or master equation in the Lindblad form is the most general type of markovian and time-homogeneous master equation describing non-unitary evolution of the density matrix \rho that is trace preserving and completely positive for any initial...
as used in quantum optics
Quantum optics
Quantum optics is a field of research in physics, dealing with the application of quantum mechanics to phenomena involving light and its interactions with matter.- History of quantum optics :...
, where it can express absorption or emission of photons from a reservoir. For example, the master equation for a single mode optical resonator (e.g. a Fabry–Perot cavity) coupled to a thermal bath is
where
is the frequency of the optical mode, is the mode annihilation operator, is the mode linewidth, and is the thermal occupation number of the photons in the bath, as given by the Bose–Einstein distribution.Derivation from Liouvillian dynamics
The derivation assumes a quantum system with a finite number of degrees of freedom coupled to a bath containing an infinite number of degrees of freedom. The system and bath each possess a Hamiltonian written in terms of operators acting only on the respective subspace of the total Hilbert space. These Hamiltonians govern the internal dynamics of the uncoupled system and bath. There is a third Hamiltonian that contains products of system and bath operators, thus coupling the system and bath. The most general form of this Hamiltonian isThe dynamics of the entire system can be described by the Liouville equation of motion,
. This equation, containing an infinite number of degrees of freedom, is impossible to solve analytically except in very particular cases. What's more, under certain approximations, the bath degrees of freedom need not be considered, and an effective master equation can be derived in terms of the system density matrix, . The problem can be analyzed more easily by moving into the interaction picture, defined by the unitary transformation , where is an arbitrary operator, and . It is straightforward to confirm that the Liouville equation becomeswhere the Hamiltonian
is explicitly time dependent. This equation can be integrated directly to giveThis implicit equation for
can be substituted back into the Liouville equation to obtain an exact differo-integral equationWe proceed with the derivation by assuming the interaction is initiated at
, and at that time there are no correlations between the system and the bath. This implies that the initial condition is factorable as , where is the density operator of the bath initially. Since the system and bath Hamiltonians act on different Hilbert subspaces, they commute, and thus in the interaction picture we can writetracing over the bath degrees of freedom,
. Tracing over the bath degrees of freedom of the aforementioned differo-integral equation yieldsThis equation is exact for the time dynamics of the system density matrix but requires full knowledge of the dynamics of the bath degrees of freedom. A simplifying assumption called the Born approximation rests on the largeness of the bath and the relative weakness of the coupling, which is to say the coupling of the system to the bath should significantly alter the bath eigenstates. In this case the full density matrix is factorable for all times as
. The master equation becomesThe equation is now explicit in the system degrees of freedom, but is very difficult to solve. A final assumption is the Born-Markoff approximation that the time derivative of the density matrix depends only on its current state, and not on its past. This assumption is valid under fast bath dynamics, wherein correlations between bath and system variables are lost extremely quickly, and amounts to replacing
on the right hand side of the equation.This is the final form of the master equation we need.
Linear Coupling
If the interaction Hamiltonian is assumed to have the formfor system operators
and bath operators , the master equation becomeswhich can be expanded as
The expectation values
are with respect to the bath degrees of freedom.