Dichteoperator

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

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

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.

Der Dichteoperator tritt im Kontext der Quantenstatistik auf.

Konstruktion[Bearbeiten]

Der Dichteoperator beschreibt einen gemischten Zustand. In dem betrachteten Ensemble befinden sich mehrere gleichartige Systeme mit den Wahrscheinlichkeiten p_i\, in den orthogonalen Zuständen |\psi_i\rangle. Sind die Zustände nicht orthogonal, so ist das jeweilige Gewichte p_i\, nicht mehr die Wahrscheinlichkeit mit der das Gemisch im jeweiligen Zustand vorliegt. Die Gewichte 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 \left|\psi_i\right\rangle\left\langle\psi_i\right| \quad (1)

Darin ist

\hat {\mathbb{P}}_{\psi_i} = |\psi_i\rangle \; \langle\psi_i|

der Projektionsoperator, der angewandt auf einen beliebigen Zustandsvektor |\Phi\rangle dessen Komponente „parallel“ zum Zustand |\psi_i\rangle herausprojiziert:

\hat {\mathbb{P}}_{\psi_i} |\Phi\rangle = |\psi_i\rangle \; \langle\psi_i|\Phi\rangle \ .

Der Faktor \langle\psi_i|\Phi\rangle darin ist die Wahrscheinlichkeitsamplitude, das im Zustand |\Phi\rangle vorliegende System im Zustand |\psi_i\rangle 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\rangle präpariert sind, ist der Dichteoperator daher einfach der Projektionsoperator selbst:

\hat {\rho} = \hat {\mathbb{P}}_{\psi} = |\psi\rangle \langle\psi|\ .

Der Dichteoperator für das kanonische Ensemble ist:

\hat {\rho} = \frac{e^{-\beta \hat H}}{\rm{Spur}\{e^{-\beta \hat H}\}}[1]

In der Eigenbasis des Hamiltonoperators nimmt \hat {\rho} die Form (1) an. Analoges erhält man für den Dichteoperator des großkanonischen Ensembles \hat \rho = \frac{e^{-\beta (\hat H-\mu \hat N)}}{\rm{Spur}\{e^{-\beta (\hat H-\mu \hat N)}\}}.

Messwerte[Bearbeiten]

Für jeden einzelnen Bestandteil |\psi_i\rangle des Zustandsgemischs ist der Mittelwert der Messergebnisse einer physikalischen Größe A gegeben durch den Erwartungswert \langle A \rangle_{\psi_i} = \langle\psi_i| \hat {A} |\psi_i\rangle \ . 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\rangle ist, ist der Mittelwert aller Messungen an den einzelnen Systemen die gewichtete Summe der einzelnen Erwartungswerte:

\langle A \rangle_{\hat {\rho}} = \sum_i \;p_i \; \langle\psi_i| \hat {A} |\psi_i\rangle \ .

Dies ist gleich der Spur

\langle A \rangle_{\hat {\rho}} = \operatorname{Tr}(\hat {\rho} \hat {A} ) \ ,

wie man mit Hilfe eines vollständigen Systems von orthonormierten Basisvektoren |\varphi_k\rangle sehen kann: Wegen \hat {1} = \sum_k |\varphi_k\rangle \langle\varphi_k| (Einheitsoperator) ist

\langle A \rangle_{\hat {\rho}} = \sum_i \;p_i \; \langle\psi_i| \hat {A}\cdot \hat {1} |\psi_i\rangle = 
\sum_{i,k} \;p_i \; \langle\psi_i|\hat {A}| \varphi_k\rangle \cdot \langle\varphi_k|\psi_i\rangle  =
\sum_{k} \langle\varphi_k| \; \left(\sum_{i} |\psi_i\rangle p_i \langle\psi_i| \hat {A}\right)\;|\varphi_k\rangle=
\sum_{k} \langle\varphi_k| \; \hat {\rho} \hat {A} \;|\varphi_k\rangle
 =\operatorname{Tr}(\hat {\rho} \hat {A} ) \ .

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

\langle A \rangle_{\hat {\rho}} =
\sum_{i,k} \;p_i \; \langle\psi_i| a_k| \varphi_k\rangle \cdot \langle\varphi_k|\psi_i\rangle  =
\sum_{k} a_k \left(\sum_{i} p_i\;\langle\varphi_k| \psi_i\rangle \; \langle\psi_i|\varphi_k\rangle \right)\; =
\sum_{k} a_k \left(\sum_{i} p_i \;|\langle\varphi_k| \psi_i\rangle|^2 \right)\; = \sum_{k} a_k P_k\ .

Darin ist P_k = \sum_{i} p_i\;|\langle\varphi_k| \psi_i\rangle|^2 \ das über das Ensemble gewichtete Mittel für die Wahrscheinlichkeit, ein herausgegriffenes System im Eigenzustand |\varphi_k\rangle 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\rangle unabhängig ist.

Beispiel: Dichteoperator und Dichtematrix für Elektronen-Polarisation[Bearbeiten]

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

\rho_{mn} = \langle \varphi_m | \hat {\rho} | \varphi_n \rangle

Basiszustände[Bearbeiten]

Im Folgenden bezeichnet das Zeichen „\overset{.}{=}“, dass ein Bra, Ket oder ein Operator bezüglich einer Basis dargestellt wird (vergleiche auch Bra-Ket#Darstellung). Die Zustände „Spin auf“ (bezgl. z-Achse) |\!\! \uparrow \rangle \overset{.}{=} \bigl(
 \begin{smallmatrix}
  1\\
  0
 \end{smallmatrix}
\bigr) und „Spin ab“ |\!\! \downarrow \rangle \overset{.}{=} \bigl(
 \begin{smallmatrix}
  0\\
  1
 \end{smallmatrix}
\bigr) werden als ket-Vektoren durch Spalten dargestellt. Die zugehörigen bra-Vektoren sind dann Zeilenvektoren: \langle  \uparrow \!| \overset{.}{=}(1\ 0) bzw. \langle \downarrow \!| \overset{.}{=} (0\ 1). Die Projektionsoperatoren (durch Matrizenmultiplikation):

\hat {P}_{\uparrow} \overset{.}{=} \bigl(
 \begin{smallmatrix}
  1\\
  0
 \end{smallmatrix}
\bigr) \cdot (1\ 0)\  = \bigl(
 \begin{smallmatrix}
  1 & 0\\
  0 & 0
 \end{smallmatrix}
\bigr)\quad, \ 
\hat {P}_{\downarrow} \overset{.}{=} \bigl(
 \begin{smallmatrix}
  0\\
  1
 \end{smallmatrix}
\bigr) \cdot (0\ 1)\ = \bigl(
 \begin{smallmatrix}
  0 & 0\\
  0 & 1
\end{smallmatrix}\bigr)

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

Polarisation in z-Richtung[Bearbeiten]

Die z-Komponente des Spins hat die aus den Eigenwerten gebildete Diagonalmatrix \hat {s_z} \overset{.}{=} \left(
 \begin{smallmatrix}
  1/2 & 0\\
  0 & -1/2
 \end{smallmatrix}
\right)\ . Für das vorausgesagte Messergebnis ergibt sich für das Ensemble \hat {P_{\uparrow}} richtig

\langle \hat {s_z} \rangle = \operatorname{Tr}(\hat P_{_\uparrow}\cdot \hat s_z) \overset{.}{=}
\operatorname{Tr}\left(
\bigl(
 \begin{smallmatrix}
  1 & 0\\
  0 & 0
 \end{smallmatrix}\bigr)
\cdot
\bigl(
 \begin{smallmatrix}
  1/2 & 0\\
  0 & -1/2
 \end{smallmatrix}
\bigr)
\right)\ =
\operatorname{Tr} \bigl(
 \begin{smallmatrix}
  1/2 & 0\\
  0 & 0
 \end{smallmatrix}\bigr)
= \tfrac{1}{2}
.

Für das Ensemble \hat P_{\downarrow} ergibt sich \langle \hat s_z \rangle = \operatorname{Tr}(\hat P_{\downarrow}\cdot \hat s_z) \overset{.}{=}
\operatorname{Tr} \bigl(
 \begin{smallmatrix}
  0 & 0\\
  0 & -1/2
 \end{smallmatrix}\bigr)
= -\tfrac{1}{2}
.

Andere Polarisionsrichtung[Bearbeiten]

Die Zustände von in +x- bzw. -x-Richtung polarisierten Elektronen sind |\!\!\rightarrow \rangle \overset{.}{=} \left(
 \begin{smallmatrix}
  \sqrt{1/2}\\
  \sqrt{1/2}
 \end{smallmatrix}\right)\;,\ 
|\!\!\leftarrow \rangle \overset{.}{=} \left(
\begin{smallmatrix}
 \sqrt{1/2}\\
  -\sqrt{1/2}
\end{smallmatrix}\right)
. Die Projektionsoperatoren dazu haben (in der Basis der s_z-Eigenzustände!) die Matrizen \hat P_{| \rightarrow \rangle} \overset{.}{=} \bigl(
 \begin{smallmatrix}
  1/2 & 1/2\\
  1/2 & 1/2
 \end{smallmatrix}\bigr)\;,\ 
\hat P_{| \leftarrow \rangle} \overset{.}{=}
\bigl(
 \begin{smallmatrix}
  1/2 & -1/2\\
  -1/2 & 1/2
 \end{smallmatrix}
\bigr)\ 
. 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[Bearbeiten]

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

\hat \rho \overset{.}{=} \tfrac{1}{2} \bigl(
 \begin{smallmatrix}
  1 & 0\\
  0 & 0
 \end{smallmatrix}
\bigr)\ +
\tfrac{1}{2} \bigl(
 \begin{smallmatrix}
  0 & 0\\
  0 & 1
 \end{smallmatrix}
\bigr) = \bigl(
 \begin{smallmatrix}
  1/2 & 0\\
  0 & 1/2
 \end{smallmatrix} \bigr)= \tfrac{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[Bearbeiten]

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 \rangle} + p_\leftarrow \;\hat P_{|\leftarrow \rangle} \overset{.}{=} 
p_\uparrow \cdot \bigl(
 \begin{smallmatrix}
  1 & 0\\
  0 & 0
 \end{smallmatrix}
\bigr)\ 
+
p_\leftarrow \cdot \bigl(
 \begin{smallmatrix}
  1/2 & -1/2\\
  -1/2 & 1/2
 \end{smallmatrix}
\bigr) =
\left(
 \begin{smallmatrix}
  p_{_\uparrow} + \tfrac{p_{_\leftarrow}}{2} & -\tfrac{p_{_\leftarrow}}{2}\\
  -\tfrac{p_{_\leftarrow}}{2} & \tfrac{p_{_\leftarrow}}{2}
 \end{smallmatrix}
\right)

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

\langle \hat s_z \rangle = \operatorname{Tr}(\hat \rho_{p_{_\uparrow},p_{_\leftarrow}} \cdot \hat s_z) \overset{.}{=}
\operatorname{Tr} \left(
 \left( \begin{smallmatrix}
  p_{_\uparrow} + \tfrac{p_{_\leftarrow}}{2} & -\tfrac{p_{_\uparrow}}{2}\\
  -\tfrac{p_{_\leftarrow}}{2} & \tfrac{p_{_\leftarrow}}{2}
 \end{smallmatrix} \right)
\cdot
\bigl(
 \begin{smallmatrix}
  1/2 & 0\\
  0 & -1/2
 \end{smallmatrix}\bigr)
\right) = 
\operatorname{Tr}\left(
\left(
 \begin{smallmatrix}
  \tfrac{1}{2} \left(p_{_\uparrow} + \tfrac{p_{_\leftarrow}}{2}\right) & \tfrac{p_{_\uparrow}}{4}\\
  -\tfrac{p_{_\leftarrow}}{4} & -\tfrac{p_{_\leftarrow}}{4}
 \end{smallmatrix} \right)\right)
= \tfrac{1}{2}p_\uparrow 
.

Die in (-x)-Richtung polarisierten Elektronen tragen also erwartungsgemäß nichts zum Erwartungswert \langle \hat s_z \rangle bei.

Formale Definition[Bearbeiten]

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 hermitesch
  2. er ist positiv semidefinit,
  3. 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)= \langle x|\hat \rho|y\rangle

die Dichtematrix in Ortsdarstellung. Sie ist allerdings keine echte Matrix, da die Ortsdarstellung über ein Kontinuum von uneigentlichen Basisvektoren |x\rangle 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}=\langle \mathbf{e}_i|\hat \rho|\mathbf{e}_j\rangle.

Eigenschaften[Bearbeiten]

  • 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\left\langle a_i\right|\hat  {\rho} \left|a_i\right\rangle=\mbox{Tr}(\hat {\mathbb{P}}_a \hat {\rho}),
wobei \left|a_i\right\rangle die orthonormierten Eigenvektoren zum Eigenwert a \!\, sind und \hat {\mathbb{P}}_a der Projektionsoperator auf den entsprechenden Eigenraum ist. Anschließend befindet sich das System im Zustand \frac{\hat{\mathbb{P}}_a \hat{\rho} \hat{\mathbb{P}}_a}{\operatorname{Tr}(\hat{\mathbb{P}}_a \hat{\rho} \hat{\mathbb{P}}_a)}
  • Der Mittelwert der Messwerte (Erwartungswert) bei Messung einer Observablen A \!\, ist
\left\langle \hat A \right\rangle=\mbox{Tr}(\hat A \hat \rho).

Dichtematrix für reine Zustände[Bearbeiten]

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

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

Dichtematrix für ein gleichverteiltes Ensemble[Bearbeiten]

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.

Reduzierter Dichteoperator[Bearbeiten]

Der reduzierte Dichteoperator wurde 1930 durch Paul Dirac eingeführt.[2][3] Er bezieht sich auf ein herausgegriffenes Teilsystem eines zusammengesetzten Systems und dient dazu, die Ergebnisse von Messungen an dem Teilsystem vorherzusagen, wenn die übrigen Teile des Systems gar nicht mit beobachtet werden.

Sind A und B zwei Systeme mit (normierten) Zuständen  |\psi_A\rangle\, , |\varphi_B\rangle in ihrem jeweiligen Hilbertraum \mathbb{H}_A,\ \mathbb{H}_B, dann hat das zusammengesetzte System A+B den Tensorraum  \mathbb{H}_A \otimes \mathbb{H}_B zum Hilbertraum. Das Gesamtsystem befindet sich in einem separablen Zustand  |\psi_A\rangle\, |\varphi_B\rangle \in  \mathbb{H}_A \otimes \mathbb{H}_B, wenn feststeht, dass die beiden Teilsysteme sich in den Zuständen  |\psi_A\rangle bzw. |\varphi_B\rangle befinden. Allgemein befindet sich das Gesamtsystem in einem Zustand

 |\Psi \rangle = \sum_{ik}\, c_{ik}\, |\psi_{Ai}\rangle\, |\varphi_{Bk}\rangle

(mit orthonormierten Basisvektoren  |\psi_{Ai}\rangle\,,\, |\varphi_{Bk}\rangle und Konstanten c_{ik} ), der als verschränkt bezeichnet wird, wenn er sich nicht als separabler Zustand darstellen lässt.

Für eine Observable des Teilsystems A ist der Operator \hat O_{\!A} zunächst nur im Hilbertraum  \mathbb{H}_A definiert. Für die Messung dieser, nur das System A betreffenden Observablen am Gesamtsystem muss der Operator gemäß \hat O_{\!A} \otimes \hat \mathbf{1}_B zu einem Operator auf  \mathbb{H}_A \otimes \mathbb{H}_B erweitert werden, wobei \hat \mathbf{1}_B der Einheitsoperator in  \mathbb{H}_B ist. Für einen separablen Zustand ergibt sich der Erwartungswert

\langle\psi_A|\, \langle\varphi_B| \,\left( \hat O_{\!A} \otimes \hat \mathbf{1}_B \right) \,|\psi_A\rangle\, |\varphi_B\rangle = \langle\psi_A|\, \hat O_{\!A} |\psi_A\rangle \cdot \langle\varphi_B|\hat \mathbf{1}_B \, |\varphi_B\rangle =  \langle\psi_A|\, \hat O_{\!A} |\psi_A\rangle \ .

Das stimmt mit dem Ergebnis überein, das man erhält, wenn man das Teilsystem A von vornherein als ein isoliertes System betrachtet. Im Allgemeinen hingegen folgt für den Erwartungswert:

\langle\Psi| \,\left( \hat O_{\!A} \otimes \hat \mathbf{1}_B \right) \,|\Psi\rangle 
= \sum_{ik\,i'k'} c_{ik} c^*_{i'k'}   \langle\psi_{Ai'}|\, \hat O_{\!A} |\psi_{Ai}\rangle \cdot \langle\varphi_{Bk'}|\hat \mathbf{1}_B \,  |\varphi_{Bk}\rangle
= \sum_{ii'} \left( \sum_{k} c_{ik} c^*_{i'k} \right)  \langle\psi_{Ai'}|\, \hat O_{\!A} |\psi_{Ai}\rangle  
= \hbox{Tr}(\hat \rho_{\!A} \, \hat O_{\!A})

Darin ist \hat \rho_{\!A} der reduzierte Dichteoperator für das System A, wenn das Gesamtsystem im Zustand \Psi ist. Er ist ein Operator im Raum \mathbb{H}_A und hat die Matrixelemente (in der Basis  |\psi_{Ai}\rangle )

 (\rho_{\!A})_{ii'} =   \sum_{k} c_{ik} c^*_{i'k}

Allgemein ausgedrückt geht der reduzierte Dichteoperator für das Teilsystem A aus dem Dichteoperator  \hat \rho_{\!A+B} für das Gesamtsystem, der die Matrixelemente  (\rho_{\!A+B})_{iki'k'} = c_{ik} c^*_{i'k'}  hat, durch Bildung der partiellen Spur über den Raum \mathbb{H}_B des Teilsystems B hervor.

Eine einfache Interpretation ergibt sich für den Fall, dass es sich bei der Basis  |\psi_{Ai}\rangle um die Eigenvektoren des Operators \hat O_{\!A} handelt (mit Eigenwerten  X_i). Dann ist der Erwartungswert von \hat O_{\!A} ein inkohärent gewichteter Mittelwert von dessen Eigenwerten:

\hbox{Tr}(\hat \rho_{\!A} \, \hat O_{\!A}) = \sum_{i} \left( \sum_{k} |c_{ik}|^2 \right)  X_i .\

Für den Fall, dass das Gesamtsystem in einem separablen Zustand vorliegt, z. B.  |\psi_{A{i_0}}\rangle |\varphi_B\rangle , ergibt diese Formel das erwartete Ergebnis \hbox{Tr}(\hat \rho_{\!A} \, \hat O_{\!A}) = X_{i_0} , denn alle Glieder mit Index i \ne i_0 sind Null, und die Summe \left( \sum_{k} |c_{i_0 k}|^2 \right)  ist die Norm von |\varphi_B\rangle, also gleich 1.

Einteilchendichteoperator[Bearbeiten]

Der Einteilchendichteoperator[4] ist bei einem Vielteilchensystem der auf den Hilbertraum eines Teilchens reduzierte Dichteoperator. Bei Systemen identischer Teilchen genügt die Kenntnis des Einteilchendichteoperators, um Erwartungswerte und Übergangsmatrixelemente jedes Operators auszurechnen, der die Summe von Einteilchenoperatoren ist. Das betrifft z. B. die kinetische Energie und die potenzielle Energie in einem äußeren Feld und ist daher ein wichtiges Hilfsmittel bei der Modellierung der Elektronenhülle von Atomen und Molekülen. Die Berechnungen werden häufig in Ortsdarstellung durchgeführt, also basierend auf der N-Teilchen-Wellenfunktion \Psi(\vec r_1, m_{s1},\,\vec r_2, m_{s2},\ldots,\,\vec r_N, m_{sN},). Darin sind \vec r_i, m_{si}, die Orts- und Spinkoordinate des i-ten Teilchens. In der Matrixdarstellung treten sie hier als z.T. kontinuierliche Indizes auf und werden deshalb nicht als unterer Index sondern wie das Argument einer Funktion geschrieben. Die Dichtematrix des Gesamtsystems heißt

 \rho(\vec r_1', m_{s1}',\,\vec r_2', m_{s2}',\ldots,\,\vec r_N', m_{sN}',\ \vec r_1, m_{s1},\,\vec r_2, m_{s2},\ldots,\,\vec r_N, m_{sN})
=\Psi^*(\vec r_1', m_{s1}',\,\vec r_2', m_{s2}',\ldots,\,\vec r_N', m_{sN}')\cdot \Psi(\vec r_1, m_{s1},\,\vec r_2, m_{s2},\ldots,\,\vec r_N, m_{sN})

Die Einteilchendichtematrix ist dann

 \rho_1(\vec r', s',\ \vec r, s) = \sum_{m_{s2},\ldots m_{sN}}\int_{dV_2 \ldots dV_N} \Psi^*(\vec r', s',\,\vec r_2, m_{s2},\ldots,\,\vec r_N, m_{sN})\cdot \Psi(\vec r, s,\,\vec r_2, m_{s2},\ldots,\,\vec r_N, m_{sN})

Die Wahl der (N-1) Integrations- (bzw. Summations-)variablen mit den Nummern 2 bis N ist beliebig, da die Wellenfunktion bei identischen Teilchen gegenüber Umnummerierung höchstens das Vorzeichen wechselt und daher für die Einteilchendichtematrix immer dasselbe Ergebnis herauskommt.

Das Diagonalelement  \rho_1(\vec r, s,\, \vec r, s) gibt die Gesamtdichte an, die die N Teilchen am Ort \vec r mit Spinrichtung m_s bilden.

Da der Einteilchendichteoperator \hat \rho_1 hermitesch ist, gibt es eine Basis  \{|\chi_n \rangle\, , n=1,2,\ldots\} aus Eigenzuständen: \hat \rho_1 |\chi_n\rangle = \lambda_n | \chi_n\rangle. Für die Eigenwerte gilt 0 \le \lambda_n \le 1 und \sum_n \lambda_n = N. Die N Eigenzustände mit den größten Eigenwerten heißen natürliche Orbitale. Wenn man jedes natürliche Orbital mit einem Teilchen besetzt, also einen Zustand in Form der Slater-Determinante bildet, stellt diese die beste Annäherung an die ursprüngliche N-Teilchen-Wellenfunktion \Psi dar, die man im Rahmen eines Einzelteilchenmodells in Bezug auf die gesamte Teilchendichte erreichen kann.

Zeitentwicklung[Bearbeiten]

Hauptartikel: Von-Neumann-Gleichung

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

\frac{\partial\hat{\rho}}{\partial t}=\frac{i}{\hbar}\left[\hat{\rho},\hat{H}\right],

wobei \hat {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 lösen und erhält mit dem unitären Zeitentwicklungs-Operator \hat U(t) = e^{-i H t/\hbar} die Gleichung

\hat  \rho(t)=\hat U(t)\;\hat \rho(0)\;\hat U^\dagger(t).

Diese Lösung kann man durch Einsetzen leicht überprüfen.

Bemerkenswert ist hierbei, dass für den Operator \hat U(t) die übliche Heisenberg'sche Bewegungsgleichung nicht gilt, da der Zeitentwicklungsoperator der direkt aus der Schrödingergleichung abgeleiteten Dynamik i\hbar \partial_t U(t)=H(t)U(t) gehorcht. Auch die Zeitentwicklung des Operators \rho durch den Zeitentwicklungsoperator \hat U(t) erfolgt nicht gemäß der üblichen Zeitentwicklungsgleichung für Operatoren (U(t)^\dagger A U(t) für eine gewöhnliche Observable A), was jedoch verständlich ist, da \hat  \rho(t)=\hat U(t)\;\hat \rho(0)\;\hat U^\dagger(t)=\sum_i p_i U(t) |\psi(0)\rangle \langle \psi(0)|U(t)^\dagger=\sum_i p_i |\psi(t)\rangle \langle \psi(t)|

Entropie[Bearbeiten]

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

S = - k_B \hbox{Tr} \left( \hat \rho \ln{\hat \rho} \right),

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[Bearbeiten]

Einzelnachweise[Bearbeiten]

  1.  Anton Amann, Ulrich Müller-Herold: Offene Quantensysteme. Springer, 2011, ISBN 9783642051876, S. 80 ff (eingeschränkte Vorschau in der Google-Buchsuche).
  2. P. A. M. Dirac: Note on Exchange Phenomena in the Thomas Atom. In: Mathematical Proceedings of the Cambridge Philosophical Society. 26, Nr. 3, 2008, S. 376. Bibcode: 1930PCPS...26..376D. doi:10.1017/S0305004100016108.
  3. U. Fano: Description of States in Quantum Mechanics by Density Matrix and Operator Techniques. In: Rev. Mod. Phys.. 29, 1957, S. 74. doi:10.1103/RevModPhys.29.74.
  4.  Frank L. Pilar: Elementary Quantum Chemistry. McGraw-Hill, NY 1968, S. 354ff..