Der Dichteoperator (auch statistischer Operator) ist ein linearer Operator, der den Zustand eines Ensembles von physikalischen Systemen, in der Quantenmechanik auch den Zustand nur eines einzigen Systems beschreibt. Diese Beschreibung ist in physikalischer Hinsicht vollständig. Das heißt, mit Hilfe des Dichteoperators lässt sich für jede am System bzw. Ensemble mögliche Messung der Erwartungswert vorhersagen.[Anm. 1] Befindet sich das System in einem Zustandsgemisch, gibt der Dichteoperator insbesondere an, mit welcher Wahrscheinlichkeit sich ein aus dem Ensemble herausgegriffenes System in einem bestimmten reinen Zustand befindet. Wird der Operator (mit Bezug auf eine Basis) als Matrix dargestellt, so spricht man von der Dichtematrix (bzw. der statistischen Matrix); diese wird in der Quantenstatistik viel verwendet.
Der Dichteoperator wurde ursprünglich im Rahmen der klassischen Physik von George Gabriel Stokes für den Polarisationszustand eines Lichtstrahls entwickelt (Stokes-Parameter). In die Quantenmechanik wurde er 1927 von Lew Landau und John von Neumann[1] eingeführt und dann ausführlich von Paul Dirac in Principles of Quantum Mechanics (1930) und von John von Neumann in Mathematische Grundlagen der Quantenmechanik (1932) dargestellt.
Konstruktion
Dichteoperator für einen reinen quantenmechanischen Zustand
Für einen reinen Zustand mit (normiertem) Zustandsvektor
heißt der Dichteoperator (als dyadisches Produkt in Bra-Ket-Schreibweise)
.
Dieser Operator bleibt ungeändert, wenn man denselben Zustand durch einen Zustandsvektor
beschrieben hätte. Daher besteht, anders als beim Zustandsvektor, eine in beiden Richtungen eindeutige Zuordnung zwischen dem physikalischen Zustand und seinem Dichteoperator.
Dieser Operator ist ein Projektionsoperator
, denn angewendet auf einen beliebigen Zustandsvektor
, projiziert er diesen auf den durch
bestimmten 1-dimensionalen Unterraum des Hilbertraums:
,
wobei der Zahlenfaktor
das Skalarprodukt beider Vektoren ist.
ist hermitesch und idempotent (d. h.
). Seine Eigenwerte sind 1 (für alle Vektoren desselben reinen Zustands) und Null (alle dazu orthogonalen Vektoren).
Für einen kohärenten, also reinen Überlagerungszustand
![{\displaystyle |\Psi \rangle =\alpha |\psi \rangle +\beta |\phi \rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/2fd6e9224139979194782b6526a38696d56e627d)
lässt sich der Dichteoperator
durch die beiden überlagerten Zustände ausdrücken (mit der komplexen Konjugation und
):
.
Wenn
und
orthogonal sind und als Basisvektoren genommen werden, dann ist
durch die Matrix
![{\displaystyle {\hat {\rho }}={\begin{pmatrix}|\alpha |^{2}&\alpha \beta ^{*}\\\alpha ^{*}\beta &|\beta |^{2}\end{pmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/775551d64643b0e5927016daf97da161981fccc0)
dargestellt. Die kohärente Linearkombination drückt sich in den Nichtdiagonalelementen aus. Diese ändern sich, wenn man die Koeffizienten
mit verschiedenen Phasenfaktoren
multipliziert, so wie dann auch der Überlagerungszustand ein anderer wird. Dieselben nichtdiagonalen Matrixelemente treten auch in der Formel für den Erwartungswert eines Operators
auf:
.
Sie bilden dort die Interferenzterme.
Dichteoperator für ein Zustandsgemisch
Ein Ensemble, das aus Subensembles zusammengesetzt ist, in denen sich jeweils die Systeme in demselben reinen Zustand
befinden, ist in einem gemischten Zustand. Hier sind die reinen Zustände inkohärent überlagert. Sind die Zustände
orthogonal, so ist die jeweilige Anzahl
der betreffenden Ensembles die Wahrscheinlichkeit dafür, dass bei einer Messung ein einzelnes System im Zustand
gefunden wird. Die Gewichte sind dann auf 1 normiert:
. Dann ist der Dichteoperator gegeben durch
.
Mit Hilfe der Projektionsoperatoren lässt sich der Dichteoperator auch schreiben als
![{\displaystyle {\hat {\rho }}=\sum _{i}p_{i}\;{\hat {\mathbb {P} }}_{\psi _{i}}\ .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b1ed6ada69c8032ce2841315fd3087fd0eb4fab9)
Der Erwartungswert eines beliebigen Operators
ist dann
![{\displaystyle \left\langle {\hat {O}}\right\rangle =\sum _{i}p_{i}\;\left\langle \psi _{i}\right|{\hat {O}}\left|\psi _{i}\right\rangle \ ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cdde14c4a44524d4a39a6eb464ce928f6c58092a)
also die inkohärente Summe der Erwartungswerte für die einzelnen Subensembles, jeweils gewichtet mit der relativen Anzahl der darin enthaltenen Einzelsysteme. Es gibt keine Interferenz zwischen den Zuständen verschiedener Einzelsysteme.
Wurde zum Beispiel das Ensemble aus zwei Subensembles zusammengesetzt, die jeweils nur Systeme in dem einen oder dem anderen von zwei orthogonalen Zuständen
und
haben, so ist der Dichteoperator
![{\displaystyle {\hat {\rho }}=p_{\psi }\;{\hat {\mathbb {P} }}_{\psi }+p_{\phi }\;{\hat {\mathbb {P} }}_{\phi }\ .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a608a89e705d7b4e09c812b971f899ec5016917)
und
mit
sind dabei die relativen Häufigkeiten.
Mit
und
als Basisvektoren ist die Dichtematrix dieses Zustandsgemischs durch die Diagonalmatrix
![{\displaystyle {\hat {\rho }}={\begin{pmatrix}p_{\psi }&0\\0&p_{\phi }\end{pmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/19d461ba33ce3ef139f8e10928459fce4641cf23)
gegeben. Die inkohärente Überlagerung von Systemen drückt sich im Verschwinden der Nichtdiagonalelemente aus, wenn (wie hier) die Systeme jeweils einen der Basiszustände besetzen.
In einer anderen Basis hat derselbe Dichteoperator im Allgemeinen eine Nichtdiagonalmatrix, ausgenommen der Fall, dass alle Basiszustände mit gleicher Häufigkeit vertreten sind.
Im Fall gleicher Häufigkeit aller
inkohärent überlagerten Basiszustände ist der Dichteoperator das
-fache des Einheitsoperators
𝟙 und hat die Matrix (hier für
:)}
![{\displaystyle {\hat {\rho }}={\begin{pmatrix}{\tfrac {1}{N}}&0\\0&{\tfrac {1}{N}}\end{pmatrix}}\ .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f48867ce4a2dafdf85dc4b97d121e36a6a2f454)
Diese Matrix ist unabhängig davon, ob innerhalb des von den beteiligten Zuständen definierten Unterraums eine andere Basis gewählt wurde. Darin drückt sich die Tatsache aus, dass inkohärente Ensembles physikalisch identisch sind, wenn sie aus orthogonalen Zuständen mit jeweils gleicher Häufigkeit, aber verschieden gewählter Basis des durch die überlagerten Zustände gebildeten Unterraums gebildet sind.
Der Dichteoperator für das kanonische Ensemble ist:
[2]
In der Eigenbasis des Hamiltonoperators nimmt
die Form (1) an. Analog erhält man für den Dichteoperator des großkanonischen Ensembles
.
Zustandsgemisch bei einem einzelnen System
Ein Zustandsgemisch liegt auch bei nur einem einzigen System vor, wenn es vor einer Messung mit einem zweiten System verschränkt war, so dass bestimmte reine Zustände des ersten Systems mit bestimmten reinen Zuständen des zweiten System vollständig korreliert waren. Wenn dann durch diese Messung, die gar nicht auf das erste System einwirkt, der Zustand des zweiten Systems zu einem bestimmten reinen Zustand reduziert wurde, der nicht als solcher zu den korrelierten Zuständen gehört hatte, muss anschließend das erste System als Zustandsgemisch behandelt werden.
Dieser Fall ist häufig, zum Beispiel wenn ein Atom ein anderes stößt, dabei mit gewisser Wahrscheinlichkeit eine Anregung verursacht und dann unter einem bestimmten Ablenkwinkel auf einen Detektor trifft. Das getroffene Atom befindet sich danach in einem Zustandsgemisch in Form einer inkohärenten Überlagerung von angeregtem Zustand und Grundzustand. Wenn man durch eine Messung am gestoßenen Atom die Richtung seines Rückstoßes festgestellt hätte, würde sich umgekehrt das stoßende Atom nun in einem Zustandsgemisch befinden, gebildet aus einer inkohärenten Überlagerung der gestreuten Wellen verschiedener Energie. Zur Beschreibung benutzt man den Reduzierten Dichteoperator, der sich aus dem Dichteoperator des ursprünglichen Gesamtsystems durch partielle Spurbildung ergibt und keine Informationen zu dem Teilsystem, an dem gemessen wurde, mehr enthält. Diese durch Verschränkung vermittelte Veränderung des Zustands eines Systems, ohne dass es Objekt einer physikalischen Einwirkung geworden wäre, stellt einen der für die Anschauung schwierigsten Aspekte der Quantenphysik dar (siehe z. B. Schrödingers Katze, Quantenverschränkung, EPR-Paradoxon, Quantenradierer).
Messwerte
Für jeden einzelnen Bestandteil
des Zustandsgemischs ist der Mittelwert der Messergebnisse einer physikalischen Größe
gegeben durch den Erwartungswert
Darin ist
der zu
gehörige Operator (s. Quantenmechanik, Observable).
Da das Ensemble ein Gemisch von Systemen in den verschiedenen beteiligten Zuständen
ist, ist der Mittelwert aller Messungen an den einzelnen Systemen die gewichtete Summe der einzelnen Erwartungswerte:
![{\displaystyle \langle A\rangle _{\hat {\rho }}=\sum _{i}\;p_{i}\;\langle \psi _{i}|{\hat {A}}|\psi _{i}\rangle \ .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6de6c1fdac73de7ff212b257f21744fb72124bbc)
Dies ist gleich der Spur
![{\displaystyle \langle A\rangle _{\hat {\rho }}=\operatorname {Tr} ({\hat {\rho }}{\hat {A}})\ ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d2b9217a2449cd94e6fd7733aed965475534e849)
wie man mit Hilfe eines vollständigen Systems von orthonormierten Basisvektoren
sehen kann: Wegen
(Einheitsoperator) ist
![{\displaystyle {\begin{aligned}\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}})\ .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c0dc094c67610b99cdb6ab9b77caac9e55e53576)
Sind die
gerade die Eigenzustände zur Observable
(d. h.
mit den Eigenwerten
), dann gilt weiter
![{\displaystyle {\begin{aligned}\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}\ .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f54d9030eb67f338a7e6ef54b65c2f40790c8c0f)
Darin ist
das über das Ensemble gewichtete Mittel für die Wahrscheinlichkeit, ein herausgegriffenes System im Eigenzustand
anzutreffen.
ist also auch die Wahrscheinlichkeit, bei einer einzelnen Messung den Eigenwert
als Ergebnis zu erhalten. Charakteristisch ist, dass
durch eine inkohärente Summe gegeben wird, die von den relativen Phasen der am Ensemble beteiligten Zustände
unabhängig ist.
Umgekehrt lässt sich der Operator
durch die aus seinen Eigenwerten
und den Dichteoperatoren der Eigenzustände
gebildete Summe darstellen:
![{\displaystyle {\hat {A}}=\sum _{k}a_{k}\;{\hat {P}}_{\varphi _{k}}\ .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad0bc9659de0069159ed3ec90d24e7020cae25cd)
Beispiel: Dichteoperator und Dichtematrix für Elektronen-Polarisation
Die Dichtematrix ist die Matrix, mit der der Operator
in Bezug auf eine orthonormierte Basis
dargestellt werden kann:
![{\displaystyle \rho _{mn}=\langle \varphi _{m}|{\hat {\rho }}|\varphi _{n}\rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/2f10812d8add07dcb0796f78acac6268e02d44c7)
Basiszustände
Im Folgenden bezeichnet das Zeichen „
“, 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)
und „Spin ab“
werden als ket-Vektoren durch Spalten dargestellt. Die zugehörigen bra-Vektoren sind dann Zeilenvektoren:
bzw.
. Die Projektionsoperatoren (durch Matrizenmultiplikation):
![{\displaystyle {\hat {P}}_{\uparrow }\doteq {\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 }\doteq {\bigl (}{\begin{smallmatrix}0\\1\end{smallmatrix}}{\bigr )}\cdot (0\ 1)\ ={\bigl (}{\begin{smallmatrix}0&0\\0&1\end{smallmatrix}}{\bigr )}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c7de62c93596d5bbf3c1e42f29f94bb3b59a2c27)
Dies sind auch die Dichtematrizen für vollständig in
- bzw.
-Richtung polarisierte Elektronen.
Polarisation in z-Richtung
Die
-Komponente des Spins hat die aus den Eigenwerten gebildete Diagonalmatrix
Für das vorausgesagte Messergebnis ergibt sich für das Ensemble
richtig
![{\displaystyle \langle {\hat {s}}_{z}\rangle =\operatorname {Tr} ({\hat {P}}_{_{\uparrow }}\cdot {\hat {s}}_{z})=\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}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6a292b65336e8ddeb972e7cd787e51fabc19bf9b)
Für das Ensemble
ergibt sich
Andere Polarisationsrichtung
Die Zustände von in
- bzw.
-Richtung polarisierten Elektronen sind
Die Projektionsoperatoren dazu haben (in der Basis der
-Eigenzustände!) die Matrizen
Charakteristisch ist, dass dies keine Diagonalmatrizen sind und dass sich die verschiedenen Phasen, mit denen die
-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
-Eigenzuständen die
-Eigenzustände gebildet werden.
Unpolarisiertes Ensemble
Sind die Elektronen je zur Hälfte in
-Richtung polarisiert, heißt die Dichtematrix:
![{\displaystyle {\hat {\rho }}\doteq {\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}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cbbb6afaed10d10c20eca8031a5341a201ba1068)
Die gleiche Dichtematrix ergibt sich für ein Gemisch aus Elektronen, die zu je 50 % in
-Richtung (oder in eine beliebige andere Richtung) polarisiert sind. Damit sind auch alle möglichen Messergebnisse identisch zu denen am Ensemble, das aus
-polarisierten Elektronen gebildet wurde. Die ursprünglichen zur Definition des Ensembles benutzten Polarisationsrichtungen sind physikalisch (und damit auch begrifflich) nicht mehr zu unterscheiden: Es ist immer ein und dasselbe Ensemble entstanden.
Gemisch verschiedener Polarisationsrichtungen
Beispielsweise für ein Gemisch aus Elektronen mit Spin in
-Richtung und
-Richtung mit Anteilen
bzw.
heißt die Dichtematrix
![{\displaystyle {\hat {\rho }}_{p_{_{\uparrow }},p_{_{\leftarrow }}}=p_{\uparrow }\;{\hat {P}}_{\left|{\uparrow }\right\rangle }+p_{\leftarrow }\;{\hat {P}}_{\left|{\leftarrow }\right\rangle }\doteq 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)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d16f6d84a66cbf3cb0b53d9ec9252fc0b00f25ef)
Der Erwartungswert des Spins in
-Richtung ist dann
![{\displaystyle \langle {\hat {s}}_{z}\rangle =\operatorname {Tr} ({\hat {\rho }}_{p_{_{\uparrow }},p_{_{\leftarrow }}}\cdot {\hat {s}}_{z})\doteq \operatorname {Tr} \left(\left({\begin{smallmatrix}p_{_{\uparrow }}+{\tfrac {p_{_{\leftarrow }}}{2}}&-{\tfrac {p_{_{\leftarrow }}}{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 }.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/607c00373d666fc6efa8e97d076e4fce38cb6634)
Die in (
)-Richtung polarisierten Elektronen tragen also erwartungsgemäß nichts zum Erwartungswert
bei.
Gegeben sei ein quantenmechanisches System, das auf einem Hilbertraum
modelliert ist. Ein beschränkter linearer Operator
auf
ist ein Dichteoperator, wenn gilt:
- er ist hermitesch
- er ist positiv semidefinit,
- 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
ein Dichteoperator, so bezeichnet
![{\displaystyle \rho (x,y)=\langle x|{\hat {\rho }}|y\rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a10e230fa508d198582bac4ef3965bd40ab08fe)
die Dichtematrix in Ortsdarstellung. Sie ist allerdings keine echte Matrix, da die Ortsdarstellung über ein Kontinuum von uneigentlichen Basisvektoren
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
gewählt wird:
.
Eigenschaften
- 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
an einem System, das durch den Dichteoperator
beschrieben wird, den Messwert
zu erhalten, ist gegeben durch
![{\displaystyle p_{a}=\sum _{i}\left\langle a_{i}\right|{\hat {\rho }}\left|a_{i}\right\rangle =\operatorname {Tr} ({\hat {\mathbb {P} }}_{a}{\hat {\rho }}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/86c5fec25b8c518de2b16df557f1d90812bb0988)
- wobei
die orthonormierten Eigenvektoren zum Eigenwert
sind und
der Projektionsoperator auf den entsprechenden Eigenraum ist. Anschließend befindet sich das System im Zustand ![{\displaystyle {\frac {{\hat {\mathbb {P} }}_{a}{\hat {\rho }}{\hat {\mathbb {P} }}_{a}}{\operatorname {Tr} ({\hat {\mathbb {P} }}_{a}{\hat {\rho }}{\hat {\mathbb {P} }}_{a})}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b1520c9286e77c510e4b8ab9815c6ff8c4f4de8e)
- Der Mittelwert der Messwerte (Erwartungswert) bei Messung einer Observablen
ist
![{\displaystyle \left\langle {\hat {A}}\right\rangle =\operatorname {Tr} ({\hat {A}}{\hat {\rho }}).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f00ac146919f456e783cb4681e15609cc80edad3)
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
.
Für gemischte Zustände gilt stets
.
Dichtematrix für ein gleichverteiltes Ensemble
Ein
-Niveau-System, bei dem alle
Zustände gleich wahrscheinlich sind, hat die Dichtematrix
![{\displaystyle {\hat {\rho }}={\frac {1}{N}}\ \mathbf {1} _{N}\ ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7a15a08422b766e6819b0e5e39f5f34b0860c269)
wobei
die
-dimensionale Einheitsmatrix bezeichnet.
Reduzierter Dichteoperator
Der reduzierte Dichteoperator wurde 1930 durch Paul Dirac eingeführt.[3][4] 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 mitbeobachtet werden.
Sind
und
zwei Systeme mit (normierten) Zuständen
in ihrem jeweiligen Hilbertraum
, dann hat das zusammengesetzte System
den Tensorraum
zum Hilbertraum. Das Gesamtsystem befindet sich in einem separablen Zustand
, wenn feststeht, dass die beiden Teilsysteme sich in den Zuständen
bzw.
befinden. Allgemein befindet sich das Gesamtsystem in einem Zustand
![{\displaystyle |\Psi \rangle =\sum _{ik}\,c_{ik}\,|\psi _{Ai}\rangle \,|\varphi _{Bk}\rangle }](https://wikimedia.org/api/rest_v1/media/math/render/svg/9993c15b57459c49ca026bc8e1eff8374cbb318e)
(mit orthonormierten Basisvektoren
und Konstanten
), der als verschränkt bezeichnet wird, wenn er sich nicht als separabler Zustand darstellen lässt.
Für eine Observable des Teilsystems
ist der Operator
zunächst nur im Hilbertraum
definiert. Für die Messung dieser, nur das System
betreffenden Observablen am Gesamtsystem muss der Operator gemäß
zu einem Operator auf
erweitert werden, wobei
der Einheitsoperator in
ist.
Ist der Zustand des Systems separabel, dann ergibt sich der Erwartungswert
![{\displaystyle \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 \ .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5aad9e7dba0879472a61ce8165be0d6106c18653)
Das stimmt mit dem Ergebnis überein, das man erhält, wenn man das Teilsystem
von vornherein als ein isoliertes System betrachtet.
Im Allgemeinen hingegen folgt für den Erwartungswert:
![{\displaystyle {\begin{aligned}\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 =\operatorname {Tr} ({\hat {\rho }}_{\!A}\,{\hat {O}}_{\!A})\,.\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/79acfaabca6af4750fe871b81d158638f6637d56)
Darin ist mit
![{\displaystyle (\rho _{\!A})_{ii'}=\sum _{k}c_{ik}c_{i'k}^{*}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9958ebe903131c972f7cb8e6d2d4844ebba072eb)
der reduzierte Dichteoperator für das Teilsystem
definiert, wenn das Gesamtsystem im Zustand
ist. Er ist ein Operator im Raum
und entsteht, wenn in der Matrix des Dichteoperators für das Gesamtsystem
![{\displaystyle (\rho _{\!A+B})_{iki'k'}=c_{ik}c_{i'k'}^{*}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/644350a2471e8ea1719eb5bf8552cf0e14b73d89)
durch Summierung über den Index
der Basiszustände des Teilsystems
die partielle Spur gebildet wird.
Eine einfache Interpretation ergibt sich für den Fall, dass es sich bei der Basis
um die Eigenvektoren des Operators
handelt (mit Eigenwerten
). Dann ist der Erwartungswert von
ein inkohärent gewichteter Mittelwert von dessen Eigenwerten:
![{\displaystyle \operatorname {Tr} ({\hat {\rho }}_{\!A}\,{\hat {O}}_{\!A})=\sum _{i}\left(\sum _{k}|c_{ik}|^{2}\right)X_{i}.\ }](https://wikimedia.org/api/rest_v1/media/math/render/svg/45cd9726908f00dc386f14ef87adfc124dd13783)
Für den Fall, dass das Gesamtsystem in einem separablen Zustand vorliegt, z. B.
, ergibt diese Formel das erwartete Ergebnis
denn alle Glieder mit Index
sind Null, und die Summe
ist die Norm von
, also gleich 1.
Einteilchendichteoperator
Der Einteilchendichteoperator[5] 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
. Darin sind
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
![{\displaystyle \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})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d4c82213fdc1363002ff99cfa2e34ec00e03249e)
![{\displaystyle =\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})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c32fca18cb4d8d5a685f5742ae2a51919be92fe3)
Die Einteilchendichtematrix ist dann
![{\displaystyle \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})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2b29f8f2f78c25c2157cb7e9d6eec2a75a31abec)
Die Wahl der (N-1) Integrations- (bzw. Summations-)variablen mit den Nummern 2 bis
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
gibt die Gesamtdichte an, die die
Teilchen am Ort
mit Spinrichtung
bilden.
Da der Einteilchendichteoperator
hermitesch ist, gibt es eine Basis
aus Eigenzuständen:
. Für die Eigenwerte gilt
und
. Die
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
dar, die man im Rahmen eines Einzelteilchenmodells in Bezug auf die gesamte Teilchendichte erreichen kann.
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
![{\displaystyle {\frac {\partial {\hat {\rho }}}{\partial t}}={\frac {i}{\hbar }}\left[{\hat {\rho }},{\hat {H}}\right],}](https://wikimedia.org/api/rest_v1/media/math/render/svg/69e0dd8d814e8df704024c1597d2f07aa9b19c51)
wobei
der Hamilton-Operator des Systems ist. Diese Gleichung ist als von-Neumannsche Bewegungsgleichung bekannt (nicht zu verwechseln mit der Heisenbergschen Bewegungsgleichung).
Diese Differentialgleichung kann man für zeitunabhängige Hamilton-Operatoren lösen und erhält mit dem unitären Zeitentwicklungs-Operator
die Gleichung
.
Diese Lösung kann man durch Einsetzen leicht überprüfen.
Bemerkenswert ist hierbei, dass für den Operator
die übliche Heisenbergsche Bewegungsgleichung nicht gilt, da der Zeitentwicklungsoperator der direkt aus der Schrödingergleichung abgeleiteten Dynamik
gehorcht. Auch die Zeitentwicklung des Operators
durch den Zeitentwicklungsoperator
erfolgt nicht gemäß der üblichen Zeitentwicklungsgleichung für Operatoren (
für eine gewöhnliche Observable A), was jedoch verständlich ist, da
Entropie
Mit Hilfe der Dichtematrix
lässt sich die Von-Neumann-Entropie eines Systems wie folgt definieren:
![{\displaystyle S=-k_{\mathrm {B} }\operatorname {Tr} \left({\hat {\rho }}\ln {\hat {\rho }}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7528903f2461ea1e1aae354590713ffcd3315f4)
wobei
die Boltzmannkonstante ist, und die Spur über dem Raum
genommen ist, in dem
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 (Auf Niederländisch)
Anmerkungen
- ↑ Auch wenn die Einzelmesswerte streuen, wird ihre Streuung durch den Dichteoperator vorhergesagt.
Einzelnachweise
- ↑ J. von Neumann: Wahrscheinlichkeitstheoretischer Aufbau der Quantenmechanik. In: Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse. Band 1927, 1927, S. 245–272 (eudml.org [abgerufen am 14. März 2023]).
- ↑ Anton Amann, Ulrich Müller-Herold: Offene Quantensysteme. Springer, 2011, ISBN 978-3-642-05187-6, S. 80 ff. (eingeschränkte Vorschau in der Google-Buchsuche).
- ↑ P. A. M. Dirac: Note on Exchange Phenomena in the Thomas Atom. In: Mathematical Proceedings of the Cambridge Philosophical Society. 26. Jahrgang, Nr. 3, 1930, S. 376, doi:10.1017/S0305004100016108, bibcode:1930PCPS...26..376D.
- ↑ U. Fano: Description of States in Quantum Mechanics by Density Matrix and Operator Techniques. In: Rev. Mod. Phys. 29. Jahrgang, 1957, S. 74, doi:10.1103/RevModPhys.29.74.
- ↑ Frank L. Pilar: Elementary Quantum Chemistry. McGraw-Hill, NY 1968, S. 354 ff.