Dichtematrix

In der Physik beschreibt die Dichtematrix bzw. der Dichteoperator (auch statistische Matrix bzw. statistischer Operator) für ein Ensemble gleichartiger Systeme, mit welcher Wahrscheinlichkeit sich ein herausgegriffenes System in einem bestimmten Zustand befindet.

Der Dichteoperator wurde ursprünglich im Rahmen der klassischen Physik von Stokes für den Polarisationszustand eines Lichtstrahls entwickelt (Stokes-Parameter). Im folgenden wird der Dichteoperator im Zusammenhang mit Systemen in quantenmechanisch definierten Zuständen dargestellt.

Konstruktion

In dem betrachteten Ensemble befinden sich mehrere gleichartige Systeme mit den Wahrscheinlichkeiten $p_i$ in den Zuständen $|{\psi }_i⟩$. Diese Zustände müssen nicht notwendigerweise orthogonal sein. Die Wahrscheinlichkeiten sind auf 1 normiert: $\sum _i{p}_i=1.$ Dann ist (in bra-ket-Schreibweise) der Dichteoperator gegeben durch

$\hat{\rho }=\sum _i{p}_i|{\psi }_i⟩⟨{\psi }_i|.$

Darin ist

${\hat{\mathbb{P}}}_{{\psi }_i}=|{\psi }_i⟩⟨{\psi }_i|$

der Projektionsoperator, der angewandt auf einen beliebigen Zustandsvektor $|\mathrm{\Phi }⟩$ dessen Komponente „parallel“ zum Zustand $|{\psi }_i⟩$ herausprojiziert:

${\hat{\mathbb{P}}}_{{\psi }_i}|\mathrm{\Phi }⟩=|{\psi }_i⟩⟨{\psi }_i|\mathrm{\Phi }⟩.$

Der Faktor $⟨{\psi }_i|\mathrm{\Phi }⟩$ darin ist die Wahrscheinlichkeitsamplitude, das im Zustand $|\mathrm{\Phi }⟩$ vorliegende System im Zustand $|{\psi }_i⟩$ vorzufinden.

Mit Hilfe der Projektionsoperatoren lässt sich der Dichteoperator auch schreiben als

$\hat{\rho }=\sum _i{p}_i{\hat{\mathbb{P}}}_{{\psi }_i}.$

Für ein Ensemble, in dem alle Systeme im selben Zustand $|\psi ⟩$ präpariert sind, ist der Dichteoperator daher einfach der Projektionsoperator selbst:

$\hat{\rho }={\hat{\mathbb{P}}}_{\psi }=|\psi ⟩⟨\psi |.$

Messwerte

Für jeden einzelnen Bestandteil $|{\psi }_i⟩$ des Zustandsgemischs ist der Mittelwert der Messergebnisse einer physikalischen Größe $A$ gegeben durch den Erwartungswert $⟨A{⟩}_{{\psi }_i}=⟨{\psi }_i|\hat{A}|{\psi }_i⟩.$ Darin ist $\hat{A}$ der zu $A$ gehörige Operator (s. Quantenmechanik, Observable).

Da das Ensemble ein Gemisch von Systemen in den verschiedenen beteiligten Zuständen $|{\psi }_i⟩$ ist, ist der Mittelwert aller Messungen an den einzelnen Systemen die gewichtete Summe der einzelnen Erwartungswerte:

$⟨A{⟩}_{\hat{\rho }}=\sum _i{p}_i⟨{\psi }_i|\hat{A}|{\psi }_i⟩.$

Dies ist gleich der Spur

$⟨A{⟩}_{\hat{\rho }}=\mathrm{Tr}(\hat{\rho }\hat{A}),$

wie man mit Hilfe eines vollständigen Systems von orthonormierten Basisvektoren $|{\phi }_k⟩$ sehen kann: Wegen $\hat{1}=\sum _k|{\phi }_k⟩⟨{\phi }_k|$ (Einheitsoperator) ist

$⟨A{⟩}_{\hat{\rho }}=\sum _i{p}_i⟨{\psi }_i|\hat{A}\cdot \hat{1}|{\psi }_i⟩=\sum _{i,k}{p}_i⟨{\psi }_i|\hat{A}|{\phi }_k⟩\cdot ⟨{\phi }_k|{\psi }_i⟩=\sum _k⟨{\phi }_k|(\sum _i|{\psi }_i⟩p_i⟨{\psi }_i|\hat{A})|{\phi }_k⟩=\sum _k⟨{\phi }_k|\hat{\rho }\hat{A}|{\phi }_k⟩=\mathrm{Tr}(\hat{\rho }\hat{A}).$

Sind die $|{\phi }_k⟩$ gerade die Eigenzustände zur Observable $A$ (d. h. $\hat{A}|{\phi }_k⟩=a_k|{\phi }_k⟩$ mit den Eigenwerten $a_k$), dann gilt weiter

$⟨A{⟩}_{\hat{\rho }}=\sum _{i,k}{p}_i⟨{\psi }_i|a_k|{\phi }_k⟩\cdot ⟨{\phi }_k|{\psi }_i⟩=\sum _k{a}_k(\sum _i{p}_i⟨{\phi }_k|{\psi }_i⟩⟨{\psi }_i|{\phi }_k⟩)=\sum _k{a}_k(\sum _i{p}_i|⟨{\phi }_k|{\psi }_i⟩{|}^2)=\sum _k{a}_k{P}_k.$

Darin ist $P_k=\sum _i{p}_i|⟨{\phi }_k|{\psi }_i⟩{|}^2$ das über das Ensemble gewichtete Mittel für die Wahrscheinlichkeit, ein herausgegriffenes System im Eigenzustand $|{\phi }_k⟩$ anzutreffen. $P_k$ ist also auch die Wahrscheinlichkeit, bei einer einzelnen Messung den Eigenwert $a_k$ als Ergebnis zu erhalten. Charakteristisch ist, dass $P_k$ durch eine inkohärente Summe gegeben wird, die von den relativen Phasen der am Ensemble beteiligten Zustände $|{\psi }_i⟩$ unabhängig ist.

Beispiel: Dichteoperator und Dichtematrix für Elektronen-Polarisation

Die Dichtematrix ist die Matrix, mit der der Operator $\hat{\rho }$ in Bezug auf eine orthonormierte Basis $|{\phi }_k⟩$ dargestellt werden kann:

${\rho }_{mn}=⟨{\phi }_m|\hat{\rho }|{\phi }_n⟩$

Basiszustände

Die Zustände „Spin auf“ (bezgl. z-Achse) $|\uparrow ⟩=(\begin{array}{c}1\\ 0\end{array})$ und „Spin ab“ $|\downarrow ⟩=(\begin{array}{c}0\\ 1\end{array})$ werden als ket-Vektoren durch Spalten dargestellt. Die zugehörigen bra-Vektoren sind dann Zeilenvektoren: $⟨\uparrow |=(10)$ bzw. $⟨\downarrow |=(01)$. Die Projektionsoperatoren (durch Matrizenmultiplikation):

${\hat{P}}_{\uparrow }=(\begin{array}{c}1\\ 0\end{array})\cdot (10)=(\begin{array}{cc}1& 0\\ 0& 0\end{array}),{\hat{P}}_{\downarrow }=(\begin{array}{c}0\\ 1\end{array})\cdot (01)=(\begin{array}{cc}0& 0\\ 0& 1\end{array})$

Dies sind auch die Dichtematrizen für vollständig in $+z$- bzw. $-z$-Richtung polarisierte Elektronen.

Polarisation in z-Richtung

Die z-Komponente des Spins hat die aus den Eigenwerten gebildete Diagonalmatrix ${\hat{s}}_z=\left(\begin{array}{cc}1/2& 0\\ 0& -1/2\end{array}\right).$ Für das vorausgesagte Messergebnis ergibt sich für das Ensemble ${\hat{P}}_{\uparrow }$ richtig

$⟨{\hat{s}}_z⟩=\mathrm{Tr}({\hat{P}}_{{}_{\uparrow }}\cdot {\hat{s}}_z)=\mathrm{Tr}((\begin{array}{cc}1/2& 0\\ 0& -1/2\end{array})\cdot (\begin{array}{cc}1& 0\\ 0& 0\end{array}))=\mathrm{Tr}(\begin{array}{cc}1/2& 0\\ 0& 0\end{array})=\frac{1}{2}.$

Für das Ensemble ${\hat{P}}_{\downarrow }$ ergibt sich $⟨{\hat{s}}_z⟩=\mathrm{Tr}({\hat{P}}_{\downarrow }\cdot {\hat{s}}_z)=\mathrm{Tr}(\begin{array}{cc}0& 0\\ 0& -1/2\end{array})=-\frac{1}{2}.$

Andere Polarisionsrichtung

Die Zustände von in $+x$- bzw. $-x$-Richtung polarisierten Elektronen sind $|\to ⟩=\left(\begin{array}{c}\sqrt{1/2}\\ \sqrt{1/2}\end{array}\right),|\leftarrow ⟩=\left(\begin{array}{c}\sqrt{1/2}\\ -\sqrt{1/2}\end{array}\right).$ Die Projektionsoperatoren dazu haben (in der Basis der $s_z$-Eigenzustände!) die Matrizen ${\hat{P}}_{|\to ⟩}=(\begin{array}{cc}1/2& 1/2\\ 1/2& 1/2\end{array}),{\hat{P}}_{|\leftarrow ⟩}=(\begin{array}{cc}1/2& -1/2\\ -1/2& 1/2\end{array}).$ Charakteristisch ist, dass dies keine Diagonalmatrizen sind und dass sich die verschiedenen Phasen, mit denen die $s_z$-Eigenzustände als ket-Vektoren hier überlagert wurden, in den Matrixelementen außerhalb der Hauptdiagonale wiederfinden. Das ist Ausdruck der kohärenten Überlagerung, durch die aus $s_z$-Eigenzuständen die $s_x$-Eigenzustände gebildet werden.

Unpolarisiertes Ensemble

Sind die Elektronen je zur Hälfte in $\pm z$-Richtung polarisiert, heißt die Dichtematrix:

$\hat{\rho }=\frac{1}{2}(\begin{array}{cc}1& 0\\ 0& 0\end{array})+\frac{1}{2}(\begin{array}{cc}0& 0\\ 0& 1\end{array})=(\begin{array}{cc}1/2& 0\\ 0& 1/2\end{array})=\frac{1}{2}\cdot \hat{1}$

Die gleiche Dichtematrix ergibt sich für ein Gemisch aus Elektronen, die zu je 50 % in $\pm x$-Richtung polarisiert sind. Damit sind auch alle möglichen Messergebnisse identisch zu denen am Ensemble, das aus $\pm z$-polarisierten Elektronen gebildet wurde. Die ursprünglichen Polarisationsrichtungen sind physikalisch (und damit auch begrifflich) nicht mehr zu unterscheiden: Es ist beide Male ein und dasselbe Ensemble geworden.

Gemisch verschiedener Polarisationsrichtungen

Für ein Gemisch aus Elektronen mit Spin in $+z$-Richtung und $-x$-Richtung (Anteile $p_{\uparrow }$ bzw. $p_{\leftarrow }$), heißt die Dichtematrix

${\hat{\rho }}_{p_{{}_{\uparrow }},p_{{}_{\leftarrow }}}=p_{\uparrow }{\hat{P}}_{|\uparrow ⟩}+p_{\leftarrow }{\hat{P}}_{|\leftarrow ⟩}=p_{\uparrow }\cdot (\begin{array}{cc}1& 0\\ 0& 0\end{array})+p_{\leftarrow }\cdot (\begin{array}{cc}1/2& -1/2\\ -1/2& 1/2\end{array})=\left(\begin{array}{cc}{p}_{{}_{\uparrow }}+\frac{p_{{}_{\leftarrow }}}{2}& -\frac{p_{{}_{\leftarrow }}}{2}\\ -\frac{p_{{}_{\leftarrow }}}{2}& \frac{p_{{}_{\leftarrow }}}{2}\end{array}\right)$

Der Erwartungswert des Spins in $\pm z$-Richtung ist dann

$⟨{\hat{s}}_z⟩=\mathrm{Tr}({\hat{\rho }}_{p_{{}_{\uparrow }},p_{{}_{\leftarrow }}}\cdot {\hat{s}}_z)=\mathrm{Tr}(\left(\begin{array}{cc}{p}_{{}_{\uparrow }}+\frac{p_{{}_{\leftarrow }}}{2}& -\frac{p_{{}_{\uparrow }}}{2}\\ -\frac{p_{{}_{\leftarrow }}}{2}& \frac{p_{{}_{\leftarrow }}}{2}\end{array}\right)\cdot (\begin{array}{cc}1/2& 0\\ 0& -1/2\end{array}))=\mathrm{Tr}\left(\left(\begin{array}{cc}\frac{1}{2}(p_{{}_{\uparrow }}+\frac{p_{{}_{\leftarrow }}}{2})& \frac{p_{{}_{\uparrow }}}{4}\\ -\frac{p_{{}_{\leftarrow }}}{4}& -\frac{p_{{}_{\leftarrow }}}{4}\end{array}\right)\right)=\frac{1}{2}{p}_{\uparrow }.$

Die in ($-x$)-Richtung polarisierten Elektronen tragen also erwartungsgemäß nichts zum Erwartungswert $⟨{\hat{s}}_z⟩$ bei.

Formale Definition

Gegeben sei ein quantenmechanisches System, das auf einem Hilbertraum $\mathbf{H}$ modelliert ist. Ein beschränkter linearer Operator $\hat{\rho }$ auf $\mathbf{H}$ ist ein Dichteoperator, wenn gilt:

1. er ist positiv semidefinit,
2. er ist Spurklasse mit Spur gleich 1.

Obwohl die Begriffe Dichtematrix und Dichteoperator oft synonym gebraucht werden, besteht ein mathematischer Unterschied. Genau wie in der linearen Algebra eine Matrix die Basisdarstellung eines linearen Operators ist, kann in der Quantenmechanik zwischen abstraktem Dichteoperator und einer konkreten Dichtematrix in einer bestimmten Darstellung unterschieden werden. Ist $\hat{\rho }$ ein Dichteoperator, so bezeichnet

$\rho (x,y)=⟨x|\hat{\rho }|y⟩$

die Dichtematrix in Ortsdarstellung. Sie ist allerdings keine echte Matrix, da die Ortsdarstellung über ein Kontinuum von uneigentlichen Basisvektoren $|x⟩$ definiert ist, sondern ein so genannter Integralkern.

In endlichdimensionalen Hilberträumen (z. B. bei Spinsystemen) ergibt sich dagegen dann eine positiv semidefinite Matrix mit Spur 1, also eine echte Dichtematrix, wenn eine Orthonormalbasis ${\mathbf{e}}_i$ gewählt wird:

${\rho }_{ij}=⟨{\mathbf{e}}_i|\hat{\rho }|{\mathbf{e}}_j⟩$

.

Eigenschaften

• Jeder Dichteoperator ist selbstadjungiert (oder hermitesch), da positive Operatoren immer selbstadjungiert sind.

• Die Menge aller Dichteoperatoren ist eine konvexe Menge, deren Rand die Menge der reinen (quantenmechanischen) Zustände ist. Die Menge ist im Gegensatz zu klassischen Theorien kein Simplex, d. h. ein Dichteoperator ist im Allgemeinen nicht eindeutig als Konvexkombination von reinen Zuständen darstellbar.

• Die Wahrscheinlichkeit, bei der Messung einer Observablen $A$ an einem System, das durch den Dichteoperator $\hat{\rho }$ beschrieben wird, den Messwert $a$ zu erhalten, ist gegeben durch

$p_a=\sum _i⟨a_i|\hat{\rho }|a_i⟩=\text{Tr}({\hat{\mathbb{P}}}_a\hat{\rho }),$

wobei

$|a_i⟩$

die orthonormierten Eigenvektoren zum Eigenwert

$a$

sind und

${\hat{\mathbb{P}}}_a$

die Projektion auf den entsprechenden Eigenraum ist.

• Der Mittelwert der Messwerte (Erwartungswert) bei Messung einer Observablen $A$ ist

$⟨\hat{A}⟩=\text{Tr}(\hat{A}\hat{\rho }).$

Dichtematrix für reine Zustände

Ist das betrachtete Ensemble ein reines Ensemble, besteht das System also nur aus einem reinen Zustand, so gilt für die Dichtematrix $\mathrm{Tr}({\hat{\rho }}^2)=\mathrm{Tr}(\hat{\rho })=1$.

Für gemischte Zustände gilt stets $\mathrm{Tr}({\hat{\rho }}^2)<1$.

Dichtematrix für ein gleichverteiltes Ensemble

Ein $N$-Niveau-System, bei dem alle $N$ Zustände gleich wahrscheinlich sind, hat die Dichtematrix

$\hat{\rho }=\frac{1}{N}{\mathbf{1}}_N,$

wobei ${\mathbf{1}}_N$ die $N$-dimensionale Einheitsmatrix bezeichnet.

Zeitentwicklung

Aus der Schrödingergleichung, die die Zeitentwicklung (Dynamik) reiner Quantenzustände beschreibt, kann man unmittelbar die Zeitentwicklung gemischter Zustände ableiten. Dazu benutzt man eine beliebige Zerlegung der Dichtematrix in reine Zustände, deren Dynamik der Schrödinger-Gleichung genügt, und berechnet daraus die Dynamik des gemischten Zustandes zu

$\dot{\hat{\rho }}=\frac{i}{\hslash }[\hat{\rho },\hat{\mathbb{H}}],$

wobei $\hat{\mathbb{H}}$ der Hamilton-Operator des Systems ist. Diese Gleichung ist als von-Neumann'sche Bewegungsgleichung bekannt (nicht zu verwechseln mit der Heisenberg'schen Bewegungsgleichung).

Diese Differentialgleichung kann man für zeitunabhängige Hamilton-Operatoren integrieren und erhält mit dem unitären Zeitentwicklungs-Operator $\hat{U}(t)=e^{-iHt/\hslash }$ die Gleichung

$\hat{\rho }(t)=\hat{U}(t)\hat{\rho }(0){\hat{U}}^{†}(t).$

Bemerkenswert ist hierbei, dass für den Operator $\hat{U}(t)$ die übliche Heisenberg'sche Bewegungsgleichung nicht gilt, da $\rho$ keine Observable ist. Auch die Transformation mit dem Operator $\hat{U}(t)$ ist nicht gemäß der üblichen Zeitentwicklungsgleichung für Operatoren.

Entropie

Mit Hilfe der Dichtematrix $\hat{\rho }$ lässt sich die Entropie eines Systems wie folgt definieren (siehe auch Von-Neumann-Entropie):

$S=-k_B\text{Tr}(\hat{\rho }\ln \hat{\rho }),$

wobei $k_B$ die Boltzmannkonstante ist, und die Spur über dem Raum $\mathbf{H}$ genommen ist, in dem $\hat{\rho }$ operiert.

Die Entropie jedes reinen Zustands ist Null, da die Eigenwerte der Dichtematrix Null und Eins sind. Dies stimmt mit der heuristischen Argumentation überein, dass keine Unsicherheit über die Präparation des Zustandes herrscht.

Man kann zeigen, dass auf einen Zustand angewendete unitäre Operatoren (wie der aus der Schrödinger-Gleichung gewonnene Zeitentwicklungs-Operator) die Entropie des Systems nicht ändern. Das verbindet die Reversibilität eines Prozesses mit seiner Entropieänderung - ein fundamentales Ergebnis, das die Quantenmechanik mit der Informationstheorie und der Thermodynamik verbindet.

Weblinks

• Artikel von Lieven Smits aus Antwerpen über De dichtheidsmatrix in de statistische mechanica (in niederländisch)