Lindblad-Gleichung

In der Quantenmechanik bezeichnet die Kossakowski-Lindblad-Gleichung (benannt nach Andrzej Kossakowski und Göran Lindblad) oder Mastergleichung in Lindblad-Form den allgemeinsten Typ einer zeit-homogenen Mastergleichung. Sie beschreibt eine nicht-unitäre Evolution des Dichteoperators <math>\rho</math>, welche spurerhaltend und komplett positiv für jede Anfangsbedingung ist.

Hintergrund

Die Lindblad-Gleichung für eine auf das <math>N</math>-dimensionale (Teil-)System reduzierte Dichtematrix <math>\rho</math> kann geschrieben werden als:

<math>\dot\rho = -\frac \mathrm i {\hbar}[H, \rho] + \sum_{n,m = 1}^{N^2 - 1} h_{n,m} \left( L_n \, \rho \, L_m^\dagger - \frac 1 2 \left( \rho \, L_m^\dagger \, L_n + L_m^\dagger \, L_n \, \rho \right) \right)</math>

Dabei bezeichnet

• der erste Summand den reversiblen Teil der Zeitentwicklung mit
  • der imaginären Einheit <math>\mathrm i</math>
  • der reduzierten Planck-Konstante <math>\hbar</math>
  • einem (hermiteschen) Hamilton-Operator <math>H</math>; <math>H</math> ist jedoch nicht notwendigerweise gleich dem Hamilton-Operator des Systems, sondern beinhaltet zusätzlich die effektive unitäre Dynamik der Wechselwirkung zwischen System und Umgebung.
• die Summe <math>\sum</math> den irreversiblen Teil mit
  • den Konstanten <math>h_{n,m}</math>, die die Dynamik festlegen. Sie bilden eine Koeffizientenmatrix <math>h = (h_{n,m})</math>, die positiv semidefinit sein muss, um sicherzustellen, dass die Gleichung spurerhaltend und komplett positiv ist.
  • den Operatoren <math>L_m</math>, die eine beliebige lineare Basis im Hilbertraum des Systems bilden.

Die Summation läuft nur über <math>N^2 - 1</math>, weil wir <math>L_{N^2}</math> proportional zum Identitätsoperator genommen haben, wodurch der Summand verschwindet. Unsere Konvention impliziert, dass die <math>L_m</math> für <math>m < N^2</math> spurlos sind.

Die Terme in der Summation, bei denen <math>m = n</math> gilt, können mit Lindblad-Superoperatoren beschrieben werden:

<math>L(C) \, \rho = C \, \rho \, C^\dagger - \frac 1 2 \left( C^\dagger \, C \, \rho + \rho \, C^\dagger \, C \right).</math>

Falls die Terme <math>h_{m,n}</math> alle Null sind, reduziert sich die Lindblad-Gleichung auf die Von-Neumann-Gleichung, das Quanten-Analogon der klassischen Liouville-Gleichung. Eine verwandte Gleichung, das Ehrenfest-Theorem, beschreibt die zeitliche Entwicklung der Erwartungswerte der Observablen.

Auch die folgenden Gleichungen für Quantenobservablen <math>A</math> werden Lindblad-Gleichungen genannt:

<math>\dot A = -\frac 1{\mathrm i\hbar} [H, A] + \frac 1{2\hbar} \sum^\infty_{k=1} \big(V^\dagger_k [A, V_k] + [V^\dagger_k, A] V_k \big)</math>

Diagonalisierung

Da die Matrix <math>h = (h_{n,m})</math> positiv semidefinit ist, kann sie mit einer unitären Transformation <math>u</math> diagonalisiert werden:

<math>u^\dagger h u = \begin{bmatrix}

\gamma_1 & 0 & \cdots & 0 \\ 0 & \gamma_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \gamma_{N^2-1} \end{bmatrix}</math>

wobei die Eigenwerte <math>\gamma_i</math> nicht negativ sind.

Wenn wir eine andere orthonormale Operator-Basis <math>A</math> definieren:

<math>A_i = \sum_{j = 1}^{N^2-1} u_{j,i} \, L_j</math>

können wir die Lindblad-Gleichung in diagonaler Form umschreiben:

<math>\dot\rho = -\frac \mathrm i {\hbar}[H,\rho] +\sum_{i = 1}^{N^2-1} \gamma_i \left( A_i \, \rho \, A_i^\dagger - \frac 1 2 \left( \rho \, A_i^\dagger \, A_i + A_i^\dagger \, A_i \, \rho \right) \right).</math>

Diese Gleichung ist invariant unter unitärer Transformation der Lindblad-Operatoren und -Konstanten,

<math> \sqrt{\gamma_i} A_i \to \sqrt{\gamma_i'} A_i' = \sum_{j = 1}^{N^2-1} v_{j,i} \sqrt{\delta_i} A_j ,</math>

und auch unter inhomogener Transformation

<math> A_i \to A_i' = A_i + a_i,</math>

<math> H \to H' = H + \frac{1}{2\mathrm i} \sum_{j=1}^{N^2-1} \gamma_j \left (a_j^* A_j - a_j A_j^\dagger \right ).</math>

Allerdings zerstört die erste Transformation die Orthonormalität der Operatoren <math>A_i</math> (solange nicht alle <math>\gamma_i</math> identisch sind) und die zweite die Spurlosigkeit. Folglich, bis auf Entartung der <math>\gamma_i</math>, sind die <math>A_i</math> der Diagonalform der Lindblad-Gleichung eindeutig bestimmt durch die Dynamik, solange wir von ihnen fordern orthonormal und spurlos zu sein.

Beispiel Harmonischer Oszillator

Ein häufiges Beispiel ist die Beschreibung der Dämpfung eines quantenmechanischen harmonischen Oszillators. Für diesen gilt

<math>\begin{align}

L_1 &= a \\ L_2 &= a^{\dagger} \\ h_{n,m} &= \begin{cases} \tfrac{\gamma}{2} \left (\bar{n}+1 \right ) & n=m=1 \\ \tfrac{\gamma}{2} \bar{n} & n=m=2 \\ 0 & \text{sonst} \end{cases} \end{align}</math>

Hier ist

• <math>\bar{n}</math> die mittlere Anzahl von Anregungen im Reservoir, die den Oszillator dämpfen, und
• <math>\gamma</math> die Zerfallsrate.

Zusätzliche Lindblad-Operatoren können hinzugefügt werden, um diverse Formen von Dephasierung und Vibrationsdämpfung (vibrational relaxation) zu modellieren. Diese Methoden sind in gitterbasierte Dichteoperator-Propagationsmethoden zur Beschreibung offener Quantensysteme aufgenommen.

Literatur

• A. Kossakowski: On quantum statistical mechanics of non-Hamiltonian systems. In: Reports on Mathematical Physics. Band 3, Nr. 4, 1972, doi:10.1016/0034-4877(72)90010-9, bibcode:1972RpMP....3..247K. 
• G. Lindblad: On the generators of quantum dynamical semigroups. In: Communications in Mathematical Physics. Band 48, Nr. 2, 1. Juni 1976, ISSN 0010-3616, S. 119–130, doi:10.1007/BF01608499, bibcode:1976CMaPh..48..119L. 
• Vittorio Gorini, Andrzej Kossakowski, E. C. G. Sudarshan: Completely positive dynamical semigroups of N‐level systems. In: Journal of Mathematical Physics. Band 17, Nr. 5, 1. Mai 1976, ISSN 0022-2488, S. 821–825, doi:10.1063/1.522979 (aip.org). 
• G. Lindblad: Non-Equilibrium Entropy and Irreversibility. Springer Verlag, 1983, ISBN 1-4020-0320-X (books.google.com). 
• Thomas Banks, Leonard Susskind, Michael E. Peskin: Difficulties for the evolution of pure states into mixed states. In: Nuclear Physics B. Band 244, Nr. 1, 1984, doi:10.1016/0550-3213(84)90184-6, bibcode:1984NuPhB.244..125B. 
• Quantum dynamical semigroups and applications. Springer Verlag, Berlin 1987, ISBN 0-387-18276-4. 
• Roman S. Ingarden, A. Kossakowski, M. Ohya: Information dynamics and open systems. Classical and quantum approach. Springer Verlag, Berlin 1997, ISBN 0-7923-4473-1. 
• Luigi Accardi, Yun Gang Lu, Igor V. Volovič: Quantum theory and its stochastic limit. Springer Verlag, Berlin / Heidelberg / New York / Barcelona / Hong Kong / London / Mailand / Paris / Tokyo 2002, ISBN 3-540-41928-4. 
• The Theory of Open Quantum Systems. Oxford University Press, New York 2002, ISBN 0-19-852063-8. 
• Open quantum systems. 2. The Markovian approach. Springer Verlag, Berlin / Heidelberg / New York 2006, ISBN 3-540-30992-6. 
• Quantum mechanics of non-Hamiltonian and dissipative systems. Elsevier Science, Amsterdam / Boston / London / New York 2008, ISBN 978-0-08-055971-1. 
• C.W. Gardiner, Peter Zoller: Quantum noise. A handbook of Markovian and non-Markovian quantum stochastic methods with applications to quantum optics (= Springer Series in Synergetics). 3. Auflage. Springer Verlag, Berlin / Heidelberg 2010, ISBN 978-3-642-06094-6. 

Weblinks

• The Lindblad master equation (Memento vom 23. Februar 2015 im Internet Archive) (englisch)