View source for Physics:Redfield equation

{{Short description|Markovian master equation of a quantum system weakly coupled to its environment}}
In [[Physics:Quantum mechanics|quantum mechanics]], the '''Redfield equation''' is a Markovian [[Quantum master equation|master equation]] that describes the time evolution of the reduced density matrix {{mvar|ρ}} of a strongly coupled quantum system that is weakly coupled to an environment. The equation is named in honor of Alfred G. Redfield, who first applied it, doing so for [[Physics:Nuclear magnetic resonance|nuclear magnetic resonance]] spectroscopy.<ref>{{Cite journal|last=Redfield|first=A.G.|date=1965 |title=The Theory of Relaxation Processes |journal=Advances in Magnetic and Optical Resonance |volume=1 |pages=1–32 |doi=10.1016/B978-1-4832-3114-3.50007-6|issn=1057-2732|isbn=978-1-4832-3114-3}}</ref> It is also known as the '''Redfield relaxation theory'''.<ref name="Poole 2012">{{cite book | last=Poole | first=Charles P. Jr | title=Relaxation in Magnetic Resonance: Dielectric and Mossbauer Applications | publisher=Elsevier Science | year=2012 | isbn=978-0-323-15182-5 |chapter=8.10 Redfield's General Relaxation Theory | chapter-url=https://books.google.com/books?id=kljr1tTZhKcC&pg=PA119 | pages=119–122}}</ref>

There is a close connection to the Lindblad master equation. If a so-called secular approximation is performed, where only certain resonant interactions with the environment are retained, every Redfield equation transforms into a master equation of Lindblad type.

Redfield equations are trace-preserving and correctly produce a thermalized state for asymptotic propagation. However, in contrast to Lindblad equations, Redfield equations do not guarantee a positive time evolution of the density matrix. That is, it is possible to get negative populations during the time evolution. The Redfield equation approaches the correct dynamics for sufficiently weak coupling to the environment.

The general form of the Redfield equation is

<math display=block>\frac{\partial}{\partial t} \rho(t) = -\frac{i}{\hbar} [H, \rho(t)] -\frac{1}{\hbar^2} \sum_m  [ S_m, (\Lambda_m \rho(t) - \rho(t) \Lambda_m^\dagger)]  </math>

where <math>H</math> is the Hermitian Hamiltonian, and the <math>S_m, \Lambda_m</math> are operators that describe the coupling to the environment, and  <math>[A,B] = AB - BA</math>  is the commutation bracket.  The explicit form is given in the derivation below.

==Derivation==

Consider a quantum system coupled to an environment with a total Hamiltonian of <math>H_\text{tot} = H + H_\text{int} + H_\text{env}</math>.  Furthermore, we assume that the interaction Hamiltonian can be written as <math>H_\text{int} = \sum_n S_n E_n</math>, where the <math>S_n</math> act only on the system degrees of freedom, the <math>E_n</math> only on the environment degrees of freedom.

The starting point of Redfield theory is the [[Physics:Nakajima–Zwanzig equation|Nakajima–Zwanzig equation]] with <math>\mathcal{P}</math> projecting on the equilibrium density operator of the environment and <math>\mathcal{Q}</math> treated up to second order.<ref>Volkhard May, Oliver Kuehn: ''Charge and Energy Transfer Dynamics in Molecular Systems.'' Wiley-VCH, 2000 {{isbn|3-527-29608-5}}</ref> An equivalent derivation starts with second-order [[Perturbation theory (quantum mechanics)|perturbation theory]] in the interaction <math>H_\text{int}</math>.<ref>Heinz-Peter Breuer, Francesco Petruccione: ''Theory of Open Quantum Systems.'' Oxford, 2002 {{isbn|978-0-19-852063-4}}</ref> In both cases, the resulting equation of motion for the density operator in the [[Physics:Interaction picture|interaction picture]] (with <math>H_{0,S} = H + H_\text{env}</math>) is

<math display=block>\frac{\partial}{\partial t} \rho_{\rm I}(t) = -\frac{1}{\hbar^2} \sum_{m,n} \int_{t_0}^t dt' \biggl(C_{mn}(t-t') \Bigl[S_{m,\mathrm{I}}(t), S_{n,\mathrm{I}}(t') \rho_{\rm I}(t')\Bigr] - C_{mn}^\ast(t-t') \Bigl[S_{m,\mathrm{I}}(t), \rho_{\rm I}(t') S_{n,\mathrm{I}}(t')\Bigr]\biggr)</math>

Here, <math>t_0</math> is some initial time, where the total state of the system and bath is assumed to be factorized, and we have introduced the bath correlation function <math>C_{mn}(t) = \text{tr}(E_{m,\mathrm{I}}(t) E_{n} \rho_\text{env,eq})</math> in terms of the density operator of the environment in thermal equilibrium, <math>\rho_\text{env,eq}</math>.

This equation is non-local in time: To get the derivative of the reduced density operator at time t, we need its values at all past times. As such, it cannot be easily solved. To construct an approximate solution, note that there are two time scales: a typical relaxation time <math>\tau_r</math> that gives the time scale on which the environment affects the system time evolution, and the coherence time of the environment, <math>\tau_c</math> that gives the typical time scale on which the correlation functions decay. If the relation

<math display=block>\tau_c \ll \tau_r</math>

holds, then the integrand becomes approximately zero before the interaction-picture density operator changes significantly. In this case, the so-called Markov approximation <math>\rho_{\rm I}(t') \approx \rho_{\rm I}(t)</math> holds. If we also move <math>t_0 \to -\infty</math> and change the integration variable <math>t' \to \tau = t - t'</math>, we end up with the Redfield master equation

<math display=block>\frac{\partial}{\partial t} \rho_{\rm I}(t) = -\frac{1}{\hbar^2} \sum_{m,n} \int_0^\infty d\tau \biggl(C_{mn}(\tau) \Bigl[S_{m,\mathrm{I}}(t), S_{n,\mathrm{I}}(t-\tau) \rho_{\rm I}(t)\Bigr] - C_{mn}^\ast(\tau) \Bigl[S_{m,\mathrm{I}}(t), \rho_{\rm I}(t) S_{n,\mathrm{I}}(t-\tau)\Bigr]\biggr)</math>

We can simplify this equation considerably if we use the shortcut <math>\Lambda_m = \sum_n \int_0^\infty d\tau C_{mn}(\tau) S_{n,\mathrm{I}}(t-\tau)</math>. In the Schrödinger picture, the equation then reads

<math display=block>\frac{\partial}{\partial t} \rho(t) = -\frac{i}{\hbar} [H, \rho(t)] -\frac{1}{\hbar^2} \sum_m [ S_m, \Lambda_m \rho(t) - \rho(t) \Lambda_m^\dagger]</math>

==Secular approximation==

Secular ({{lang-la|saeculum|lit=century}}) approximation is an approximation valid for long times <math>t</math>. The time evolution of the Redfield relaxation tensor is neglected as the Redfield equation describes weak coupling to the environment. Therefore, it is assumed that the relaxation tensor changes slowly in time, and it can be assumed constant for the duration of the interaction described by the interaction [[Physics:Hamiltonian (quantum mechanics)|Hamiltonian]]. In general, the time evolution of the reduced density matrix can be written for the element <math>ab</math> as

{{NumBlk|:|<math>\frac{\partial}{\partial t}\rho_{ab}(t) = -i\omega_{ab}\rho_{ab}(t) - \sum_{cd}\mathcal{R_{abcd}}\rho_{cd}(t)</math>|{{EquationRef|1}}}}

where <math>\mathcal{R}</math> is the time-independent Redfield relaxation tensor. 

Given that the actual coupling to the environment is weak (but non-negligible), the Redfield tensor is a small perturbation of the system Hamiltonian and the solution can be written as 

<math display=block>\rho_{ab}(t) = e^{-i\omega_{ab}t}{\rho}_{ab,\mathrm{I}}(t)</math>

where <math>\rho_{\rm I}(t)</math> is not constant but slowly changing amplitude reflecting the weak coupling to the environment. This is also a form of the [[Physics:Interaction picture|interaction picture]], hence the index "I".<ref group="note">The interaction picture describes the evolution of the density matrix in a "frame of reference" where the changes due to Hamiltonian <math>H_0</math> are not manifested. It is essentially the same transformation as entering a rotating frame of reference to solve a problem of combined rotating motion in classical mechanics. The interaction picture then describes only the envelope of the time evolution of the density matrix where only the more subtle effects of the perturbation Hamiltonian manifest. The mathematical formula for a transformation from the Schrödinger picture to the interaction picture is given by <math>\psi_{\rm I}(t)=U^{\dagger}(t)\psi_{\rm S}(t)=e^{iH_0t/\hbar}\psi_{\rm S}(t)</math>, which is the same form as this equation.</ref> 

Taking a derivative of the <math>\rho_{\rm I}(t)</math> and substituting the equation ({{EquationNote|1}}) for <math>\frac{\partial}{\partial t}\rho_{ab}(t)</math>, we are left with only the relaxation part of the equation

<math display=block>\frac{\partial}{\partial t}\rho_{ab,\mathrm{I}}(t) = - \sum_{cd}\mathcal{R_{abcd}}e^{i\omega_{ab}t-i\omega_{cd}t}\rho_{cd,\mathrm{I}}(t)</math>.

We can integrate this equation on condition that the interaction picture of the reduced density matrix <math>\rho_{\rm I}(t)</math> changes slowly in time (which is true if <math>\mathcal{R}</math> is small), then <math>\rho_{ab,\mathrm{I}}(t) \approx \rho_{ab,\mathrm{I}}(0)</math>, getting

<math display=block>\rho_{ab,\mathrm{I}}(t) = \rho_{ab,\mathrm{I}}(0) - \sum_{cd}\int_0^td\tau\mathcal{R_{abcd}}e^{i\omega_{ab}\tau-i\omega_{cd}\tau}\rho_{cd,\mathrm{I}}(t) = \rho_{ab,\mathrm{I}}(0) - \sum_{cd}\mathcal{R_{abcd}} \frac{(e^{i\Delta\omega t} - 1)}{i\Delta\omega} \rho_{cd,\mathrm{I}}(t)</math> 

where <math>\Delta\omega = \omega_{ab} - \omega_{cd}</math>.

In the limit of <math>\Delta\omega</math> approaching zero, the fraction <math>\frac{(e^{i\Delta\omega t} - 1)}{i\Delta\omega}</math> approaches <math>t</math>, therefore the contribution of one element of the reduced density matrix to another element is proportional to time (and therefore dominates for long times <math>t</math>). In case <math>\Delta\omega</math> is not approaching zero, the contribution of one element of the reduced density matrix to another oscillates with an amplitude proportional to <math>\frac{1}{\Delta\omega}</math> (and therefore is negligible for long times <math>t</math>). It is therefore appropriate to neglect any contribution from non-diagonal elements (<math>cd</math>) to other non-diagonal elements (<math>ab</math>) and from a non-diagonal elements (<math>cd</math>) to diagonal elements (<math>aa</math>, <math>a=b</math>), since the only case when frequencies of different modes are equal is the case of [[Physics:Degenerate energy levels|random degeneracy]]. The only elements left in the Redfield tensor to evaluate after the Secular approximation are therefore:

* <math>\mathcal{R}_{aabb}</math>, the transfer of population from one state to another (from <math>b</math> to <math>a</math>);
* <math>\mathcal{R}_{aaaa}</math>, the depopulation constant of state <math>a</math>; and
* <math>\mathcal{R}_{abab}</math>, the pure dephasing of the element <math>\rho_{ab}(t)</math> (dephasing of coherence).

==Notes==
{{reflist|group=note}}

==References==
{{Reflist}}
==External links==
* [http://qutip.org/docs/latest/guide/dynamics/dynamics-bloch-redfield.html brmesolve] Bloch-Redfield master equation solver from [[Software:QuTiP|QuTiP]].

[[Category:Quantum mechanics]]

{{Sourceattribution|Redfield equation}}