Macierz gęstości

Szablon:Dopracować Macierz gęstości (ang. density matrix) lub operator gęstości (ang. density operator) to matematyczna reprezentacja stanu układu kwantowego. Jest ogólniejsza od reprezentacji wektorowej, gdyż oprócz stanów czystych (reprezentowanych przez wektor) obejmuje również półklasyczne stany mieszane.

Formalizm operatorów gęstości został wprowadzony przez Johna von Neumanna w 1927Szablon:Fakt.

Tworzenie macierzy gęstości

Dla stanu czystego reprezentowanego przez wektor $|\psi ⟩$ odpowiadający mu operator to

$\rho =|\psi ⟩⟨\psi |,$

czyli operator rzutowy rzutujący na jednowymiarową podprzestrzeń $|\psi ⟩$ przestrzeni Hilberta $\mathscr{H}.$

Z kolei dla stanu mieszanego składającego się z (nieinterferujących ze sobą) składników $|{\psi }_i⟩$ odpowiadający mu operator gęstości to

$\rho =\sum _i{\lambda }_i|{\psi }_i⟩⟨{\psi }_i|,$

gdzie ${\lambda }_i$ to prawdopodobieństwa znalezienia poszczególnego składnika. Muszą one spełniać $0⩽{\lambda }_i<1$ dla każdego i oraz $\sum _i{\lambda }_i=1.$

Własności

Dla układu kwantowego opisywanego w przestrzeni Hilberta $\mathscr{H},$ operator gęstości to dowolny operator liniowy ciągły $\rho :\mathscr{H}\to \mathscr{H}$ spełniający poniższe warunki

• $\rho$ jest samosprzężony,
• $\rho$ jest dodatnio określony,
• $\mathrm{Tr}(\rho )=1$
• $\mathrm{Tr}({\rho }^2)⩽1,$ przy czym równość zachodzi wyłącznie dla stanu czystego.

Gdy układ jest opisywany w skończenie wymiarowej przestrzeni Hilberta, macierz gęstości jest rzeczywiście reprezentowana przez macierz operatora liniowego.

Zbiór wszystkich macierzy gęstości przestrzeni Hilberta $\mathscr{H}$ jest oznaczany jako $S(\mathscr{H}).$ Jest to zbiór wypukły, co oznacza, że każdy operator gęstości może być przedstawiony jako kombinacja wypukła:

$\rho =\sum _i{a}_i{\sigma }_i$

gdzie ${\sigma }_i\in S(\mathscr{H}),$ $a_i⩾0,$ dla każdego $i$ oraz $\sum _i{a}_i=1.$

Przedstawienie to jest niejednoznaczne, co oznacza że stan mieszany układu kwantowego może być zrealizowany jako próbka stanów czystych na wiele sposobów.

Stany czyste są punktami ekstremalnymi zbioru macierzy gęstości i jako takie mają jednoznaczne przedstawienie.

Równanie von Neumanna dla macierzy gęstości

Tak jak dla funkcji falowych istnieje równanie Schrödingera, również dla macierzy gęstości istnieje odpowiednie równanie zwane równaniem von Neumanna (lub Liouville’a-von Neumanna)

$i\hslash \frac{\partial \rho }{\partial t}=[\hat{H},\rho ],$

gdzie $[\hat{H},\rho ]$ to komutator hamiltonianu z macierzą gęstości.

Dzięki temu, że powyższe równanie jest liniowe, w wyprowadzeniu można ograniczyć się do stanów czystych. Istotna jest także samosprzężoność hamiltonianu.

$\begin{array}{cc}\frac{\partial |\psi ⟩⟨\psi |}{\partial t}& =\frac{\partial |\psi ⟩}{\partial t}⟨\psi |+|\psi ⟩\frac{\partial ⟨\psi |}{\partial t}\\ & =(\frac{1}{i\hslash }\hat{H}|\psi ⟩)⟨\psi |+|\psi ⟩{(\frac{1}{i\hslash }\hat{H}|\psi ⟩)}^{*}\\ & =\frac{1}{i\hslash }(\hat{H}|\psi ⟩⟨\psi |-|\psi ⟩⟨\psi |\hat{H})\\ & =\frac{1}{i\hslash }[\hat{H},|\psi ⟩⟨\psi |]\end{array}$

Obliczanie wartości oczekiwanej

Dla operatora obserwabli $\hat{A}$ wartość średnia na wektorze $|\psi ⟩$ to

$⟨\hat{A}⟩=⟨\psi |\hat{A}|\psi ⟩.$

W przypadku mieszania stanów wartość średnią operatora należy uśrednić po wszystkich stanach podlegających mieszaniu wagowaną przez prawdopodobieństwa ich wystąpienia

$⟨\hat{A}⟩=\sum _i{\lambda }_i⟨{\psi }_i|\hat{A}|{\psi }_i⟩.$

Do wnętrza powyższego wyrażenia możemy wstawić operator jednostkowy:

$⟨\hat{A}⟩=\sum _i{\lambda }_i⟨{\psi }_i|\hat{A}\sum _j|{\psi }_j⟩⟨{\psi }_j|{\psi }_i⟩.$

Możemy przestawić ${\lambda }_i$ pod znak sumy oraz zmienić indeks sumowania w (), dzięki czemu otrzymujemy:

$⟨\hat{A}⟩=\sum _i⟨{\psi }_i|\hat{A}{\lambda }_i\sum _j|{\psi }_j⟩⟨{\psi }_j|{\psi }_i⟩,$

$⟨\hat{A}⟩=\sum _i⟨{\psi }_i|\hat{A}\hat{\rho }|{\psi }_i⟩,$

$⟨\hat{A}⟩=\mathrm{tr}(\hat{A}\hat{\rho }).$

Formuła Borna-von Neumana

W wyniku pomiaru obserwabli $X$ na układzie opisanym przez operator gęstości $\rho ,$ otrzymujemy rozkład prawdopodobieństwa na przestrzeni możliwych wyników opisany wzorem

${\mu }_{\rho }^X(\omega )=\mathrm{Tr}(f_X(\omega )\rho ),$

gdzie $f_X:ℬ(R)\to X$ to rozkład spektralny obserwabli $X.$

Przykład

W odniesieniu do spektroskopii NMR operator macierzy gęstości $\hat{\rho }$ opisuje średnią statystyczną układu spinów po wszystkich stanach $\psi ,$ w których się one znajdują: $\hat{\rho }=\overline{|\psi ⟩⟨\psi |}.$ Elementy diagonalne macierzy ${\rho }_{\alpha \alpha }$ oraz ${\rho }_{\beta \beta },$ odpowiadające stanom własnym energii Zeemana $|\alpha ⟩,|\beta ⟩$ odpowiadają spinom będącym na głównych poziomach energetycznych spinów w polu magnetycznym, są to tzw. populacje. Elementy poza diagonalne ${\rho }_{\beta \alpha },$ ${\rho }_{\alpha \beta }$ w macierzy nazwane są koherencjami, odpowiadają one superpozycjom stanów.