Movatterモバイル変換


[0]ホーム

URL:


Zum Inhalt springen
WikipediaDie freie Enzyklopädie
Suche

Dichteoperator

aus Wikipedia, der freien Enzyklopädie

DerDichteoperator (auchstatistischer Operator) ist ein linearer Operator, der denZustand einesEnsembles vonphysikalischen Systemen, in derQuantenmechanik 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 derErwartungswert vorhersagen.[Anm. 1] Befindet sich das System in einemZustandsgemisch, gibt der Dichteoperator insbesondere an, mit welcherWahrscheinlichkeit sich ein aus dem Ensemble herausgegriffenes System in einem bestimmtenreinen Zustand befindet. Wird der Operator (mit Bezug auf eine Basis) als Matrix dargestellt, so spricht man von derDichtematrix (bzw. derstatistischen Matrix); diese wird in derQuantenstatistik viel verwendet.

Der Dichteoperator wurde ursprünglich im Rahmen derklassischen Physik vonGeorge Gabriel Stokes für denPolarisationszustand einesLichtstrahls entwickelt (Stokes-Parameter). In dieQuantenmechanik wurde er 1927 vonLew Landau undJohn von Neumann[1] eingeführt und dann ausführlich vonPaul Dirac inPrinciples of Quantum Mechanics (1930) und von John von Neumann inMathematische Grundlagen der Quantenmechanik (1932) dargestellt.

Konstruktion

[Bearbeiten |Quelltext bearbeiten]

Dichteoperator für einen reinen quantenmechanischen Zustand

[Bearbeiten |Quelltext bearbeiten]

Für einen reinen Zustand mit (normiertem)Zustandsvektor|ψ{\displaystyle |\psi \rangle } heißt der Dichteoperator (alsdyadisches Produkt inBra-Ket-Schreibweise)

ρ^=|ψψ|{\displaystyle {\hat {\rho }}=\left|\psi \right\rangle \left\langle \psi \right|}.

Dieser Operator bleibt ungeändert, wenn man denselben Zustand durch einen Zustandsvektoreiφ|ψ{\displaystyle \mathrm {e} ^{\mathrm {i} \varphi }|\psi \rangle } 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 einProjektionsoperatorρ^=P^ψ{\displaystyle {\hat {\rho }}={\hat {\mathbb {P} }}_{\psi }}, denn angewendet auf einen beliebigen Zustandsvektor|ϕ{\displaystyle |\phi \rangle }, projiziert er diesen auf den durch|ψ{\displaystyle |\psi \rangle } bestimmten 1-dimensionalen Unterraum des Hilbertraums:

ρ^|ϕ=|ψψ|ϕ{\displaystyle {\hat {\rho }}|\phi \rangle =\left|\psi \right\rangle \left\langle \psi \right|\phi \rangle },

wobei der Zahlenfaktorψ|ϕ{\displaystyle \langle \psi |\phi \rangle } dasSkalarprodukt beider Vektoren ist.ρ^{\displaystyle {\hat {\rho }}} ist hermitesch und idempotent (d. h.ρ^2=ρ^{\displaystyle {\hat {\rho }}^{2}={\hat {\rho }}}). Seine Eigenwerte sind 1 (für den ket-Vektor des reinen Zustands und alle seine Vielfachen) und Null (für alle dazu orthogonalen Vektoren).

Für einen kohärenten, also reinen Überlagerungszustand

|Ψ=α|ψ+β|ϕ{\displaystyle |\Psi \rangle =\alpha |\psi \rangle +\beta |\phi \rangle }

lässt sich der Dichteoperatorρ^=|ΨΨ|{\displaystyle {\hat {\rho }}=\left|\Psi \right\rangle \left\langle \Psi \right|} durch die beiden überlagerten Zustände ausdrücken (mit derkomplexen Konjugation undxx=|x|2{\displaystyle x^{*}x=|x|^{2}}):

ρ^=|ΨΨ|=(α|ψ+β|ϕ)(αψ|+βϕ|)=|α|2|ψψ|+αβ|ψϕ|+αβ|ϕψ|+|β|2|ϕϕ|{\displaystyle {\begin{aligned}{\hat {\rho }}&=\left|\Psi \right\rangle \left\langle \Psi \right|\\&=\left(\alpha |\psi \rangle +\beta |\phi \rangle \right)\left(\alpha ^{*}\langle \psi |+\beta ^{*}\langle \phi |\right)\\&=|\alpha |^{2}|\psi \rangle \langle \psi |+\alpha \beta ^{*}|\psi \rangle \langle \phi |+\alpha ^{*}\beta |\phi \rangle \langle \psi |+|\beta |^{2}|\phi \rangle \langle \phi |\end{aligned}}}.

Wenn|ψ{\displaystyle |\psi \rangle } und|ϕ{\displaystyle |\phi \rangle } orthogonal sind und alsBasisvektoren genommen werden, dann istρ^{\displaystyle {\hat {\rho }}} durch die Matrix

ρ^=(|α|2αβαβ|β|2){\displaystyle {\hat {\rho }}={\begin{pmatrix}|\alpha |^{2}&\alpha \beta ^{*}\\\alpha ^{*}\beta &|\beta |^{2}\end{pmatrix}}}

dargestellt. Dass die Überlagerung in Form einer Linearkombination eine kohärente Überlagerung ist, drückt sich in den Nichtdiagonalelementen aus. Diese ändern sich, wenn man die Koeffizientenα, β{\displaystyle \alpha ,\ \beta } mitverschiedenen Phasenfaktoreneiφ1,2{\displaystyle \mathrm {e} ^{\mathrm {i} \varphi _{1,2}}} multipliziert, so dass sich ihre relative Phase ändert. Auch der Überlagerungszustand wird dann ein anderer, obwohl die Anteile|α|2, |β|2{\displaystyle |\alpha |^{2},\ |\beta |^{2}} der beiden Basiszustände gleich bleiben.

Dieselben nichtdiagonalen Matrixelemente machen auch den Erwartungswert eines beliebigen OperatorsO^{\displaystyle {\hat {O}}} und damit alle Messergebnisse von der relativen Phase abhängig:

Ψ|O^|Ψ=|α|2ψ|O^|ψ+|β|2ϕ|O^|ϕ+αβψ|O^|ϕ+αβϕ|O^|ψ{\displaystyle \left\langle \Psi \right|{\hat {O}}\left|\Psi \right\rangle =|\alpha |^{2}\left\langle \psi \right|{\hat {O}}\left|\psi \right\rangle +|\beta |^{2}\left\langle \phi \right|{\hat {O}}\left|\phi \right\rangle +\alpha ^{*}\beta \left\langle \psi \right|{\hat {O}}\left|\phi \right\rangle +\alpha \beta ^{*}\left\langle \phi \right|{\hat {O}}\left|\psi \right\rangle }.

Die nichtdiagonalen Matrixelemente bilden dort dieInterferenzterme.

Dichteoperator für ein Zustandsgemisch

[Bearbeiten |Quelltext bearbeiten]

Ein Ensemble gleicher Systeme, die sich in verschiedenen reinen Zuständen|ψi  (i=1,2,){\displaystyle |\psi _{i}\rangle \ \ (i=1,2,\ldots )} befinden, ist einZustandsgemisch. Es wird durch einen Dichteoperator

ρ^=ipiP^ψi{\displaystyle {\hat {\rho }}=\sum _{i}p_{i}\;{\hat {\mathbb {P} }}_{\psi _{i}}}

beschrieben. InBra-Ket-Schreibweise:

ρ^=ipi|ψiψi|(1){\displaystyle {\hat {\rho }}=\sum _{i}p_{i}\left|\psi _{i}\right\rangle \left\langle \psi _{i}\right|\quad (1)}.

Die Gewichtepi{\displaystyle p_{i}} sind reelle Zahlen mit0pi1{\displaystyle 0\leq p_{i}\leq 1} undipi=1{\displaystyle \sum _{i}p_{i}=1}. Als reelle Linearkombination von hermiteschen Projektionsoperatoren ist auch der Dichteoperator hermitesch.

Die komplexen Phasen der Zustände|ψi{\displaystyle |\psi _{i}\rangle } haben keinen Einfluss auf den Dichteoperator. Dementsprechend ist für einen beliebigen OperatorO^{\displaystyle {\hat {O}}} der Erwartungswertdie Summe der Erwartungswerte für die einzelnen Zustände, jeweils gewichtet mit der relativen Häufigkeit bei der Präparation des Gemischs.

O^=ipiψi|O^|ψi .{\displaystyle \left\langle {\hat {O}}\right\rangle =\sum _{i}p_{i}\;\left\langle \psi _{i}\right|{\hat {O}}\left|\psi _{i}\right\rangle \ .}

Die reinen Zustände sindinkohärent überlagert, es gibt keine Interferenzterme zwischen den verschiedenen überlagerten Zuständen.

Für das folgende wird angenommen, dass die Zustände|ψi{\displaystyle |\psi _{i}\rangle } orthogonal sind. Dann ist Gleichung (1) dieSpektraldarstellung des Dichteoperators, die|ψi{\displaystyle |\psi _{i}\rangle } sind seine Eigenzustände und diepi{\displaystyle p_{i}} die Eigenwerte dazu.

Dann istpi{\displaystyle p_{i}} die Wahrscheinlichkeit dafür, dass bei einer Messung ein einzelnes System im Zustand|ψi{\displaystyle |\psi _{i}\rangle } gefunden wird. Denn der Operator für diese Messung istO^i=|ψiψi|{\displaystyle {\hat {O}}_{i}=\left|\psi _{i}\right\rangle \left\langle \psi _{i}\right|} und der ErwartungswertO^i=ipiψi|ψiψi|ψi =ipiδii=pi{\displaystyle \left\langle {\hat {O}}_{i}\right\rangle =\sum _{i'}p_{i'}\;\langle \psi _{i'}|\psi _{i}\rangle \;\langle \psi _{i}|\psi _{i'}\rangle \ =\sum _{i'}p_{i'}\delta _{ii'}=p_{i}}.

In der Basis seiner Eigenzustände hat der Dichteoperator eine Diagonalmatrix mit den Gewichtenpi{\displaystyle p_{i}} auf der Hauptdiagonale. Weil für allepi2pi{\displaystyle p_{i}^{2}\leq p_{i}} gilt, istρ^2ρ^{\displaystyle {\hat {\rho }}^{2}\leq {\hat {\rho }}}, wobei das Gleichheitszeichen nur dann gilt, wenn genau einpi=1{\displaystyle p_{i}=1} und alle anderenpi=0{\displaystyle p_{i'}=0}, wenn also ein reiner Zustand vorliegt.

Sind allepi{\displaystyle p_{i}} einander gleich, gilt wegen der Normierungpi=1N{\displaystyle p_{i}={\tfrac {1}{N}}} mitN{\displaystyle N} als der Anzahl verschiedener Zustände, die zusammengemischt wurden und eine Basis des betreffenden Teilraums des Zustandsraums des Systems bilden. In diesem Fall der vollständigenEntartung ist die Dichtematrix das1N{\displaystyle {\tfrac {1}{N}}}-Fache der Einheitsmatrix, so dass ein Zustandsgemisch aus anderen Basiszuständen dieses Teilraums die gleiche Dichtematrix haben und daher physikalisch dasselbe Zustandsgemisch beschreiben würde. Es ist daher in seinem solchen Gemisch unmöglich zu ermitteln, aus welchen Basiszuständen das Gemisch zusammengesetzt ist.

Wurde zum Beispiel ein Gemisch nur aus zwei Basiszuständen|ψ1{\displaystyle |\psi _{1}\rangle } und|ψ2{\displaystyle |\psi _{2}\rangle } zusammengesetzt, so ist der Dichteoperator im Allgemeinen

ρ^=p1P^ψ1+p2P^ψ2 .{\displaystyle {\hat {\rho }}=p_{1}\;{\hat {\mathbb {P} }}_{\psi _{1}}+p_{2}\;{\hat {\mathbb {P} }}_{\psi _{2}}\ .}

p1,2{\displaystyle p_{1,2}} mitp1+p2=1{\displaystyle p_{1}+p_{2}=1} sind dabei die Gewichte oder relativen Häufigkeiten.

Die Dichtematrix dieses Zustandsgemischs ist durch die Diagonalmatrix

ρ^=(p100p2){\displaystyle {\hat {\rho }}={\begin{pmatrix}p_{1}&0\\0&p_{2}\end{pmatrix}}}

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, deren Zustände dann kohärente Überlagerungen der ursprünglichen Basiszustände sind, hat derselbe Dichteoperator im Allgemeinen auch Elemente außerhalb der Hauptdiagonale. Ausgenommen ist nur der Fall vollständiger Entartung, bei dem alle Basiszustände mit gleicher Häufigkeit vertreten sind.

Im entgegengesetzten Fall sind alle Gewichtepi{\displaystyle p_{i}} voneinander verschieden. Dann liegen die Eigenzustände eindeutig fest, und es lässt sich (prinzipiell) durch Messungen eindeutig erkennen, welche Basiszustände nach Gleichung (1) das Zustandsgemisch bilden. Dann ist auch die aus der klassischen Physik geläufige Vorstellung zulässig, dass ein aus dem Ensemble herausgegriffenes Objekt mit Sicherheit genau einen dieser Basiszustände einnimmt und nicht etwa in einem Zustand ihrer kohärenten Überlagerung ist.[2] Insofern verhält sich ein solches Zustandsgemisch wie ein Gemisch vieler gleicher Systeme in der alltäglichen Anschauung oder der klassischen Physik. Allerdings muss man in der Quantenmechanik oft auch ein einzelnes System als Zustandsgemisch beschreiben. Dies ist z. B. notwendig, wenn ein System isoliert betrachtet werden soll, nachdem es mit einem anderen System in Wechselwirkung war und das aus beiden Systemen zusammengesetzte Gesamtsystem dabei in einen reinen, aber verschränkten Zustand übergegangen ist. Dieser Fall tritt bei jederquantenmechanischen Messung ein.

Dichteoperator für ein kanonisches Ensemble

[Bearbeiten |Quelltext bearbeiten]

Der Dichteoperator für daskanonische Ensemble ist:

ρ^=eβH^Spur{eβH^}.{\displaystyle {\hat {\rho }}={\frac {\mathrm {e} ^{-\beta {\hat {H}}}}{\rm {{Spur}\{e^{-\beta {\hat {H}}}\}}}}.}[3]

In der Eigenbasis des Hamiltonoperators nimmtρ^{\displaystyle {\hat {\rho }}} die Form (1) an. Analog erhält man für den Dichteoperator des großkanonischen Ensembles

ρ^=eβ(H^μN^)Spur{eβ(H^μN^)}.{\displaystyle {\hat {\rho }}={\frac {\mathrm {e} ^{-\beta ({\hat {H}}-\mu {\hat {N}})}}{\rm {{Spur}\{\mathrm {e} ^{-\beta ({\hat {H}}-\mu {\hat {N}})}\}}}}.}

Zustandsgemisch bei einem einzelnen System

[Bearbeiten |Quelltext bearbeiten]

Ein Zustandsgemisch liegt auch bei nur einem einzigen System vor, wenn es vor einer Messung mit einem zweiten System zu einem Gesamtsystem verschränkt war, so dass bestimmte reine Zustände des ersten Systems mit bestimmten reinen Zuständen des zweiten Systems 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 reinenZustand 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 getroffenen 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 weiteren Beschreibung des einen Teilsystems benutzt man denReduzierten Dichteoperator, der sich aus dem vollen Dichteoperator des ursprünglichen Gesamtsystems durch partielle Spurbildung über das andere Teilsystem ergibt und danach keine Informationen über dieses 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 Aspekte der Quantenphysik dar, die für die Anschauung am schwierigsten sind (siehe z. B.Schrödingers Katze,Quantenverschränkung,EPR-Paradoxon,Quantenradierer).

Messwerte

[Bearbeiten |Quelltext bearbeiten]

Für jeden einzelnen Bestandteil|ψi{\displaystyle |\psi _{i}\rangle } des Zustandsgemischs ist der Mittelwert der Messergebnisse einer physikalischen GrößeA{\displaystyle A} gegeben durch den ErwartungswertAψi=ψi|A^|ψi .{\displaystyle \langle A\rangle _{\psi _{i}}=\langle \psi _{i}|{\hat {A}}|\psi _{i}\rangle \ .}Darin istA^{\displaystyle {\hat {A}}} der zuA{\displaystyle A} gehörige Operator (s.Quantenmechanik,Observable).

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

Aρ^=ipiψi|A^|ψi .{\displaystyle \langle A\rangle _{\hat {\rho }}=\sum _{i}\;p_{i}\;\langle \psi _{i}|{\hat {A}}|\psi _{i}\rangle \ .}

Dies ist gleich der Spur

Aρ^=Tr(ρ^A^) ,{\displaystyle \langle A\rangle _{\hat {\rho }}=\operatorname {Tr} ({\hat {\rho }}{\hat {A}})\ ,}

wie man mit Hilfe eines vollständigen Systems von orthonormierten Basisvektoren|φk{\displaystyle |\varphi _{k}\rangle } sehen kann: Wegen1^=k|φkφk|{\displaystyle {\hat {1}}=\sum _{k}|\varphi _{k}\rangle \langle \varphi _{k}|} (Einheitsoperator) ist

Aρ^=ipiψi|A^1^|ψi=i,kpiψi|A^|φkφk|ψi=kφk|(i|ψipiψi|A^)|φk=kφk|ρ^A^|φk=Tr(ρ^A^) .{\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}}}

Sind die|φk{\displaystyle |\varphi _{k}\rangle } gerade die Eigenzustände zur ObservableA{\displaystyle A} (d. h.A^|φk=ak|φk{\displaystyle {\hat {A}}|\varphi _{k}\rangle =a_{k}|\varphi _{k}\rangle } mit den Eigenwertenak{\displaystyle a_{k}}), dann gilt weiter

Aρ^=i,kpiψi|ak|φkφk|ψi=kak(ipiφk|ψiψi|φk)=kak(ipi|φk|ψi|2)=kakPk .{\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}}}

Darin istPk=ipi|φk|ψi|2 {\displaystyle 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|φk{\displaystyle |\varphi _{k}\rangle } anzutreffen.Pk{\displaystyle P_{k}} ist also auch die Wahrscheinlichkeit, bei einer einzelnen Messung den Eigenwertak{\displaystyle a_{k}} als Ergebnis zu erhalten. Charakteristisch ist, dassPk{\displaystyle P_{k}} durch eineinkohärente Summe gegeben wird, die von den relativen Phasen der am Ensemble beteiligten Zustände|ψi{\displaystyle |\psi _{i}\rangle } unabhängig ist.

Umgekehrt lässt sich der OperatorA^{\displaystyle {\hat {A}}} durch die aus seinen Eigenwertenak{\displaystyle a_{k}} und den Dichteoperatoren der EigenzuständeP^φk=|φkφk|{\displaystyle {\hat {P}}_{\varphi _{k}}=|\varphi _{k}\rangle \langle \varphi _{k}|} gebildete Summe darstellen:

A^=kakP^φk .{\displaystyle {\hat {A}}=\sum _{k}a_{k}\;{\hat {P}}_{\varphi _{k}}\ .}

Beispiel: Dichteoperator und Dichtematrix für Elektronen-Polarisation

[Bearbeiten |Quelltext bearbeiten]

DieDichtematrix ist die Matrix, mit der der Operatorρ^{\displaystyle {\hat {\rho }}} in Bezug auf eine orthonormierte Basis|φk{\displaystyle |\varphi _{k}\rangle } dargestellt werden kann:

ρmn=φm|ρ^|φn.{\displaystyle \rho _{mn}=\langle \varphi _{m}|{\hat {\rho }}|\varphi _{n}\rangle .}

Basiszustände

[Bearbeiten |Quelltext bearbeiten]

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

P^(10)(1 0) =(1000), P^(01)(0 1) =(0001).{\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 )}.}

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

Polarisation in z-Richtung

[Bearbeiten |Quelltext bearbeiten]

Diez{\displaystyle z}-Komponente des Spins hat die aus den Eigenwerten gebildete Diagonalmatrixs^z(1/2001/2) .{\displaystyle {\hat {s}}_{z}\doteq \left({\begin{smallmatrix}1/2&0\\0&-1/2\end{smallmatrix}}\right)\ .} Für das vorausgesagte Messergebnis ergibt sich für das EnsembleP^{\displaystyle {\hat {P_{\uparrow }}}} richtig

s^z=Tr(P^s^z)=Tr((1000)(1/2001/2)) =Tr(1/2000)=12.{\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}}.}

Für das EnsembleP^{\displaystyle {\hat {P}}_{\downarrow }} ergibt sichs^z=Tr(P^s^z)=Tr(0001/2)=12.{\displaystyle \langle {\hat {s}}_{z}\rangle =\operatorname {Tr} ({\hat {P}}_{\downarrow }\cdot {\hat {s}}_{z})=\operatorname {Tr} {\bigl (}{\begin{smallmatrix}0&0\\0&-1/2\end{smallmatrix}}{\bigr )}=-{\tfrac {1}{2}}.}

Andere Polarisationsrichtung

[Bearbeiten |Quelltext bearbeiten]

Die Zustände von in+x{\displaystyle +x}- bzw.x{\displaystyle -x}-Richtung polarisierten Elektronen sind|(1/21/2), |(1/21/2).{\displaystyle \left|{\rightarrow }\right\rangle \doteq \left({\begin{smallmatrix}{\sqrt {1/2}}\\{\sqrt {1/2}}\end{smallmatrix}}\right)\;,\ \left|{\leftarrow }\right\rangle \doteq \left({\begin{smallmatrix}{\sqrt {1/2}}\\-{\sqrt {1/2}}\end{smallmatrix}}\right).}Die Projektionsoperatoren dazu haben (in der Basis dersz{\displaystyle s_{z}}-Eigenzustände!) die MatrizenP^|(1/21/21/21/2), P^|(1/21/21/21/2) .{\displaystyle {\hat {P}}_{\left|{\rightarrow }\right\rangle }\doteq {\bigl (}{\begin{smallmatrix}1/2&1/2\\1/2&1/2\end{smallmatrix}}{\bigr )}\;,\ {\hat {P}}_{\left|{\leftarrow }\right\rangle }\doteq {\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 diesz{\displaystyle s_{z}}-Eigenzustände als ket-Vektoren hier überlagert wurden, in den Matrixelementen außerhalb derHauptdiagonale wiederfinden. Das ist Ausdruck der kohärenten Überlagerung, durch die aussz{\displaystyle s_{z}}-Eigenzuständen diesx{\displaystyle s_{x}}-Eigenzustände gebildet werden.

Unpolarisiertes Ensemble

[Bearbeiten |Quelltext bearbeiten]

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

ρ^12(1000) +12(0001)=(1/2001/2)=121^.{\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}}.}

Die gleiche Dichtematrix ergibt sich für ein Gemisch aus Elektronen, die zu je 50 % in±x{\displaystyle \pm x}-Richtung (oder in eine beliebige andere Richtung) polarisiert sind. Damit sind auch alle möglichen Messergebnisse identisch zu denen am Ensemble, das aus±z{\displaystyle \pm z}-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 immerein und dasselbe Ensemble entstanden.

Gemisch verschiedener Polarisationsrichtungen

[Bearbeiten |Quelltext bearbeiten]

Beispielsweise für ein Gemisch aus Elektronen mit Spin in(+z){\displaystyle (+z)}-Richtung und(x){\displaystyle (-x)}-Richtung mit Anteilenp{\displaystyle p_{\uparrow }} bzw.p{\displaystyle p_{\leftarrow }} heißt die Dichtematrix

ρ^p,p=pP^|+pP^|p(1000) +p(1/21/21/21/2)=(p+p2p2p2p2).{\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).}

Der Erwartungswert des Spins in±z{\displaystyle \pm z}-Richtung ist dann

s^z=Tr(ρ^p,ps^z)Tr((p+p2p2p2p2)(1/2001/2))=Tr((12(p+p2)p4p4p4))=12p.{\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 }.}

Die in (x{\displaystyle -x})-Richtung polarisierten Elektronen tragen also erwartungsgemäß nichts zum Erwartungswerts^z{\displaystyle \langle {\hat {s}}_{z}\rangle } bei.

Formale Definition

[Bearbeiten |Quelltext bearbeiten]

Gegeben sei einquantenmechanisches System, das auf einemHilbertraum H{\displaystyle \mathbf {H} } modelliert ist.Ein beschränkterlinearer Operatorρ^{\displaystyle {\hat {\rho }}} aufH{\displaystyle \mathbf {H} } ist ein Dichteoperator, wenn gilt:

  1. er isthermitesch,
  2. er istpositiv semidefinit,
  3. er istSpurklasse mit Spur gleich 1.

Obwohl die BegriffeDichtematrix undDichteoperator oft synonym gebraucht werden, besteht ein mathematischer Unterschied. Genau wie in derlinearen Algebra eineMatrix die Basisdarstellung eines linearen Operators ist, kann in der Quantenmechanik zwischen abstraktemDichteoperator und einer konkretenDichtematrix in einer bestimmten Darstellung unterschieden werden. Istρ^{\displaystyle {\hat {\rho }}} ein Dichteoperator, so bezeichnet

ρ(x,y)=x|ρ^|y{\displaystyle \rho (x,y)=\langle x|{\hat {\rho }}|y\rangle }

die Dichtematrix inOrtsdarstellung. Sie ist allerdings keine echte Matrix, da die Ortsdarstellung über ein Kontinuum vonuneigentlichenBasisvektoren|x{\displaystyle |x\rangle } definiert ist, sondern ein so genannterIntegralkern.

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

ρij=ei|ρ^|ej{\displaystyle \rho _{ij}=\langle \mathbf {e} _{i}|{\hat {\rho }}|\mathbf {e} _{j}\rangle }.

Eigenschaften

[Bearbeiten |Quelltext bearbeiten]
  • Die Menge aller Dichteoperatoren ist einekonvexe Menge, derenRand die Menge der reinen (quantenmechanischen) Zustände ist. Die Menge ist im Gegensatz zu klassischen Theorien keinSimplex, d. h. ein Dichteoperator ist im Allgemeinen nicht eindeutig alsKonvexkombination von reinen Zuständen darstellbar.
pa=iai|ρ^|ai=Tr(P^aρ^),{\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 }}),}
wobei|ai{\displaystyle \left|a_{i}\right\rangle } die orthonormiertenEigenvektoren zumEigenwert a{\displaystyle a} sind undP^a{\displaystyle {\hat {\mathbb {P} }}_{a}} derProjektionsoperator auf den entsprechendenEigenraum ist. Anschließend befindet sich das System im ZustandP^aρ^P^aTr(P^aρ^P^a).{\displaystyle {\frac {{\hat {\mathbb {P} }}_{a}{\hat {\rho }}{\hat {\mathbb {P} }}_{a}}{\operatorname {Tr} ({\hat {\mathbb {P} }}_{a}{\hat {\rho }}{\hat {\mathbb {P} }}_{a})}}.}
A^=Tr(A^ρ^).{\displaystyle \left\langle {\hat {A}}\right\rangle =\operatorname {Tr} ({\hat {A}}{\hat {\rho }}).}

Dichtematrix für reine Zustände

[Bearbeiten |Quelltext bearbeiten]

Besteht dasEnsemble nur aus Systemen im selben reinen Zustand, so gilt für die DichtematrixTr(ρ^2)=Tr(ρ^)=1{\displaystyle \operatorname {Tr} \,({\hat {\rho }}^{2})=\operatorname {Tr} \,({\hat {\rho }})=1}.

Für Zustandsgemische gilt stetsTr(ρ^2)<1{\displaystyle \operatorname {Tr} \,({\hat {\rho }}^{2})<1}.

Dichtematrix für ein gleichverteiltes Ensemble

[Bearbeiten |Quelltext bearbeiten]

EinN{\displaystyle N}-Niveau-System, bei dem alleN{\displaystyle N} Zustände gleich wahrscheinlich sind, hat die Dichtematrix

ρ^=1N 1N ,{\displaystyle {\hat {\rho }}={\frac {1}{N}}\ \mathbf {1} _{N}\ ,}

wobei1N{\displaystyle \mathbf {1} _{N}} dieN{\displaystyle N}-dimensionaleEinheitsmatrix bezeichnet.

Reduzierter Dichteoperator

[Bearbeiten |Quelltext bearbeiten]

Der reduzierte Dichteoperator 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.[4] Der reduzierte Dichteoperator wurde 1930 durchPaul Dirac eingeführt.[5]

SindA{\displaystyle A} undB{\displaystyle B} zwei Systeme mit (normierten) Zuständen|ψA,|φB{\displaystyle |\psi _{A}\rangle \,,|\varphi _{B}\rangle } in ihrem jeweiligen HilbertraumHA, HB{\displaystyle \mathbb {H} _{A},\ \mathbb {H} _{B}}, dann hat das zusammengesetzte SystemA+B{\displaystyle A+B} denTensorraumHAHB{\displaystyle \mathbb {H} _{A}\otimes \mathbb {H} _{B}} zum Hilbertraum. Das Gesamtsystem befindet sich in einemseparablen Zustand|ψA|φBHAHB{\displaystyle |\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|ψA{\displaystyle |\psi _{A}\rangle } bzw.|φB{\displaystyle |\varphi _{B}\rangle } befinden. Allgemein befindet sich das Gesamtsystem in einem Zustand

|Ψ=ikcik|ψAi|φBk{\displaystyle |\Psi \rangle =\sum _{ik}\,c_{ik}\,|\psi _{Ai}\rangle \,|\varphi _{Bk}\rangle }

(mit orthonormierten Basisvektoren|ψAi,|φBk{\displaystyle |\psi _{Ai}\rangle \,,\,|\varphi _{Bk}\rangle } und Konstantencik{\displaystyle c_{ik}} ), der alsverschränkt bezeichnet wird, wenn er sichnicht als separabler Zustand darstellen lässt.

Für eine Observable des TeilsystemsA{\displaystyle A} ist der OperatorO^A{\displaystyle {\hat {O}}_{\!A}} zunächst nur im HilbertraumHA{\displaystyle \mathbb {H} _{A}} definiert. Für die Messung dieser, nur das SystemA{\displaystyle A} betreffenden Observablen am Gesamtsystem muss der Operator gemäßO^A1^B{\displaystyle {\hat {O}}_{\!A}\otimes {\hat {\mathbf {1} }}_{B}} zu einem Operator aufHAHB{\displaystyle \mathbb {H} _{A}\otimes \mathbb {H} _{B}} erweitert werden, wobei1^B{\displaystyle {\hat {\mathbf {1} }}_{B}} der Einheitsoperator inHB{\displaystyle \mathbb {H} _{B}} ist.

Ist der Zustand des Systems separabel, dann ergibt sich der Erwartungswert

ψA|φB|(O^A1^B)|ψA|φB=ψA|O^A|ψAφB|1^B|φB=ψA|O^A|ψA .{\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 \ .}

Das stimmt mit dem Ergebnis überein, das man erhält, wenn man das TeilsystemA{\displaystyle A} von vornherein als ein isoliertes System betrachtet.

Im Allgemeinen hingegen folgt für den Erwartungswert:

Ψ|(O^A1^B)|Ψ=ikikcikcikψAi|O^A|ψAiφBk|1^B|φBk=ii(kcikcik)ψAi|O^A|ψAi=Tr(ρ^AO^A).{\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}}}

Darin ist mit

(ρA)ii=kcikcik{\displaystyle (\rho _{\!A})_{ii'}=\sum _{k}c_{ik}c_{i'k}^{*}}

der reduzierte Dichteoperator für das TeilsystemA{\displaystyle A} definiert, wenn das Gesamtsystem im ZustandΨ{\displaystyle \Psi } ist. Er ist ein Operator im RaumHA{\displaystyle \mathbb {H} _{A}} und entsteht aus der Matrix des Dichteoperators für das Gesamtsystem

(ρA+B)ikik=cikcik{\displaystyle (\rho _{\!A+B})_{iki'k'}=c_{ik}c_{i'k'}^{*}},

wenn nur die Glieder mitk=k{\displaystyle k=k'} betrachtet werden und durch Summierung über den Indexk=k{\displaystyle k=k'} der Basiszustände des TeilsystemsB{\displaystyle B} einepartielle Spur gebildet wird.

Der reduzierte Dichteoperator hängt nicht von der Wahl des OperatorsO^A{\displaystyle {\hat {O}}_{\!A}} ab. Daher gilt die Formel

Ψ|(O^A1^B)|Ψ=Tr(ρ^AO^A){\displaystyle \langle \Psi |\,\left({\hat {O}}_{\!A}\otimes {\hat {\mathbf {1} }}_{B}\right)\,|\Psi \rangle =\operatorname {Tr} ({\hat {\rho }}_{\!A}\,{\hat {O}}_{\!A})}

mit demselbenρ^A{\displaystyle {\hat {\rho }}_{\!A}} für jeden OperatorO^A{\displaystyle {\hat {O}}_{\!A}}. Mithin ermöglicht der reduzierte Dichteoperator die Berechnung sämtlicher Erwartungswerte von Observablen, die nur das TeilsystemA{\displaystyle A} betreffen, und stellt damit die vollständige Beschreibung des Zustands dieses Teilsystems dar.

Eine einfache Interpretation des reduzierten Dichteoperators ergibt sich, wenn man bei gegebenemO^A{\displaystyle {\hat {O}}_{\!A}} (mit EigenwertenXi{\displaystyle X_{i}}) für die Basiszustände|ψAi{\displaystyle |\psi _{Ai}\rangle } die Eigenvektoren dieses Operators wählt. Dann ist der Erwartungswert vonO^A{\displaystyle {\hat {O}}_{\!A}} ein inkohärent gewichteter Mittelwert von dessen Eigenwerten:

Tr(ρ^AO^A)=i(k|cik|2)Xi. {\displaystyle \operatorname {Tr} ({\hat {\rho }}_{\!A}\,{\hat {O}}_{\!A})=\sum _{i}\left(\sum _{k}|c_{ik}|^{2}\right)X_{i}.\ }

Für den Fall, dass das TeilsystemA{\displaystyle A} sich in einem dieser Eigenzustände|ψAi0{\displaystyle |\psi _{A{i_{0}}}\rangle } befindet, so dass das Gesamtsystem in einem separablen Zustand vorliegt, z. B.|ψAi0|φB{\displaystyle |\psi _{A{i_{0}}}\rangle |\varphi _{B}\rangle }, ergibt diese Formel das erwartete ErgebnisTr(ρ^AO^A)=Xi0,{\displaystyle \operatorname {Tr} ({\hat {\rho }}_{\!A}\,{\hat {O}}_{\!A})=X_{i_{0}},} denn alle Glieder mit Indexii0{\displaystyle i\neq i_{0}} sind Null, und die Summe(k|ci0k|2){\displaystyle \left(\sum _{k}|c_{i_{0}k}|^{2}\right)} ist die Norm von|φB{\displaystyle |\varphi _{B}\rangle }, also gleich 1.

Einteilchendichteoperator

[Bearbeiten |Quelltext bearbeiten]

Der Einteilchendichteoperator[6] 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 derN-Teilchen-WellenfunktionΨ(r1,ms1,r2,ms2,,rN,msN,){\displaystyle \Psi ({\vec {r}}_{1},m_{s1},\,{\vec {r}}_{2},m_{s2},\ldots ,\,{\vec {r}}_{N},m_{sN},)}. Darin sindri,msi,{\displaystyle {\vec {r}}_{i},m_{si},} die Orts- und Spinkoordinate desi-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

ρ(r1,ms1,r2,ms2,,rN,msN, r1,ms1,r2,ms2,,rN,msN){\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})}
=Ψ(r1,ms1,r2,ms2,,rN,msN)Ψ(r1,ms1,r2,ms2,,rN,msN).{\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}).}

Die Einteilchendichtematrix ist dann

ρ1(r,s, r,s)=ms2,msNdV2dVNΨ(r,s,r2,ms2,,rN,msN)Ψ(r,s,r2,ms2,,rN,msN).{\displaystyle \rho _{1}({\vec {r}}',s',\ {\vec {r}},s)=\sum _{m_{s2},\ldots m_{sN}}\int _{\mathrm {d} V_{2}\ldots \mathrm {d} V_{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 bisN{\displaystyle 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ρ1(r,s,r,s){\displaystyle \rho _{1}({\vec {r}},s,\,{\vec {r}},s)} gibt die Gesamtdichte an, die dieN{\displaystyle N} Teilchen am Ortr{\displaystyle {\vec {r}}} mit Spinrichtungms{\displaystyle m_{s}} bilden.

Da der Einteilchendichteoperatorρ^1{\displaystyle {\hat {\rho }}_{1}} hermitesch ist, gibt es eine Basis{|χn,n=1,2,}{\displaystyle \{|\chi _{n}\rangle \,,n=1,2,\ldots \}} aus Eigenzuständen:ρ^1|χn=λn|χn{\displaystyle {\hat {\rho }}_{1}|\chi _{n}\rangle =\lambda _{n}|\chi _{n}\rangle }. Für die Eigenwerte gilt0λn1{\displaystyle 0\leq \lambda _{n}\leq 1} undnλn=N{\displaystyle \sum _{n}\lambda _{n}=N}. DieN{\displaystyle N} Eigenzustände mit den größten Eigenwerten heißennatürliche Orbitale. Wenn man jedes natürliche Orbital mit einem Teilchen besetzt, also einen Zustand in Form derSlater-Determinante bildet, stellt diese die beste Annäherung an die ursprünglicheN-Teilchen-WellenfunktionΨ{\displaystyle \Psi } dar, die man im Rahmen eines Einzelteilchenmodells in Bezug auf die gesamte Teilchendichte erreichen kann.

Zeitentwicklung

[Bearbeiten |Quelltext bearbeiten]
Hauptartikel:Von-Neumann-Gleichung

Aus derSchrödingergleichung, die die Zeitentwicklung (Dynamik) reiner Quantenzustände beschreibt, kann man unmittelbar die Zeitentwicklung einesZustandsgemischs 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 Zustandsgemischs zu

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

wobeiH^{\displaystyle {\hat {H}}} derHamilton-Operator des Systems ist. Diese Gleichung ist alsvon-Neumannsche Bewegungsgleichung bekannt (nicht zu verwechseln mit derHeisenbergschen Bewegungsgleichung).

DieseDifferentialgleichung kann man für zeitunabhängige Hamilton-Operatoren lösen und erhält mit demunitärenZeitentwicklungs-OperatorU^(t)=eiHt/{\displaystyle {\hat {U}}(t)=\mathrm {e} ^{-\mathrm {i} Ht/\hbar }} die Gleichung

ρ^(t)=U^(t)ρ^(0)U^(t){\displaystyle {\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 OperatorU^(t){\displaystyle {\hat {U}}(t)} die üblicheHeisenbergsche Bewegungsgleichungnicht gilt, da der Zeitentwicklungsoperator der direkt aus der Schrödingergleichung abgeleiteten DynamikitU(t)=H(t)U(t){\displaystyle \mathrm {i} \hbar \partial _{t}U(t)=H(t)U(t)} gehorcht. Auch die Zeitentwicklung des Operatorsρ{\displaystyle \rho } durch den ZeitentwicklungsoperatorU^(t){\displaystyle {\hat {U}}(t)} erfolgt nicht gemäß der üblichen Zeitentwicklungsgleichung für Operatoren (U(t)AU(t){\displaystyle U(t)^{\dagger }AU(t)} für eine gewöhnliche Observable A), was jedoch verständlich ist, daρ^(t)=U^(t)ρ^(0)U^(t)=ipiU(t)|ψ(0)ψ(0)|U(t)=ipi|ψ(t)ψ(t)|.{\displaystyle {\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 |Quelltext bearbeiten]

Mit Hilfe des Dichteoperatorsρ^{\displaystyle {\hat {\rho }}} lässt sich dieVon-Neumann-Entropie eines Systems wie folgt definieren:

S=kBTr(ρ^lnρ^),{\displaystyle S=-k_{\mathrm {B} }\operatorname {Tr} \left({\hat {\rho }}\ln {\hat {\rho }}\right),}

wobeikB{\displaystyle k_{\mathrm {B} }} dieBoltzmannkonstante ist, und die Spur über dem RaumH{\displaystyle \mathbf {H} } genommen ist, in demρ^{\displaystyle {\hat {\rho }}} operiert.

Die Entropie jedes reinen Zustands ist Null, da als Eigenwerte des Dichteoperators nur 0 und 1 vorkommen. Dies stimmt mit derheuristischen Argumentation überein, dass keine Unsicherheit über die Präparation des Zustandes herrscht. Bei allen echten Zustandsgemischen liegen mindestens zwei Eigenwerte des Dichteoperators zwischen 0 und 1, Gemische haben daher positive Entropie.

Die zeitliche Entwicklung des Systems gemäß einer Schrödingergleichung lässt die Entropie konstant, weil der Dichteoperator dabei immer eine unitäre Transformation durchläuft. Das stellt eine Verbindung zwischenReversibilität eines Prozesses und seiner eventuellen Entropieänderung her – ein fundamentales Ergebnis, das die Quantenmechanik mit derInformationstheorie und derThermodynamik verbindet. Ein echtes Zustandsgemisch kann demnach auf diesem Weg nicht in einen reinen Zustand übergehen.[7]

Weblinks

[Bearbeiten |Quelltext bearbeiten]

Anmerkungen

[Bearbeiten |Quelltext bearbeiten]
  1. Auch wenn die Einzelmesswerte streuen, wird ihre Streubreite ebenso wie alle weiteren Charakteristika ihrer Verteilung durch den Dichteoperator vorhergesagt.

Einzelnachweise

[Bearbeiten |Quelltext bearbeiten]
  1. 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]). 
  2. Siehe Ref. Fano 1957, Kap. 3(e). Es gilt auch bei jedem hermiteschen Operator mit nicht-entarteten Eigenwerten, dass die Basis aus Eigenvektoren eindeutig bestimmt ist. S.Claude Cohen-Tannoudji:Quantenmechanik. de Gruyter, 1999,ISBN 3-11-016458-2., Kap. 2.4.1
  3. 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). 
  4. 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. 
  5. 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. 
  6. Frank L. Pilar:Elementary Quantum Chemistry. McGraw-Hill, NY 1968,S. 354 ff. 
  7. J. v. Neumann:Mathematische Grundlagen der Quantenmechanik. Springer (1932, 1968, 1996), Kap. V.3 .
Abgerufen von „https://de.wikipedia.org/w/index.php?title=Dichteoperator&oldid=252062880
Kategorien:
Versteckte Kategorie:

[8]ページ先頭

©2009-2025 Movatter.jp