Струм ймовірності

В квантовій механіці, струм ймовірності (або потік ймовірності) описує зміни функції щільності ймовірності.

Визначення

Струм ймовірності ${\vec {j}}$ визначається як

{\vec {j}}={\frac {\hbar }{2mi}}\left(\Psi ^{*}{\vec {\nabla }}\Psi -\Psi {\vec {\nabla }}\Psi ^{*}\right)={\frac {\hbar }{m}}{\mbox{Im}}(\Psi ^{*}{\vec {\nabla }}\Psi )

та задовліьняє квантово-механічне рівняння неперервності

{\frac {\partial \rho }{\partial t}}+{\vec {\nabla }}\cdot {\vec {j}}=0

зі щільністю ймовірності $\rho$ , заданою

\rho =|\Psi |^{2}

.

Рівняння неперевності є еквівалетним наступному інтегральному рівнянню:

{\frac {\partial }{\partial t}}\int \limits _{V}|\Psi |^{2}dV+\int \limits _{S}{\vec {j}}\cdot {\vec {dS}}=0

де $V$ — объём и $S$ − межа об'єму $V$ . Це закон збуреження для щільності ймовірності в квантовій механіці.

Зокрема, якщо $\Psi$ — хвильова функція окремої частинки, інтеграл в першому доданку попереднього рівняння (без похідної по часу) — ймовірність отримання значення в межах $V$ , коли стан частинки виміряно. Другий доданок — швидкість, з якою ймовірність «витікає» з об'єму $V$ .

Загалом рівняння свідсить, що похідна по часу ймовірності знаходження частинки в $V$ дорівнює швидкості, по якій ймовірність «витікає» з $V$ .

Приклади

Плоска хвиля

Струм ймовірності, який можна зіставити плоскій хвилі

\Psi =Ae^{i{\vec {k}}\cdot {\vec {r}}}e^{i\omega t}

запишеться у вигляді

{\vec {j}}={\frac {\hbar }{2mi}}|A|^{2}\left(e^{-i{\vec {k}}\cdot {\vec {r}}}{\vec {\nabla }}e^{i{\vec {k}}\cdot {\vec {r}}}-e^{i{\vec {k}}\cdot {\vec {r}}}{\vec {\nabla }}e^{-i{\vec {k}}\cdot {\vec {r}}}\right)=|A|^{2}{\frac {\hbar {\vec {k}}}{m}}.

Це частка квадрата амплітуди на швидкість частинки:

{\vec {v}}={\frac {\vec {p}}{m}}={\frac {\hbar {\vec {k}}}{m}}

.

Зауважте, що струм ймовірності є відмінним від нуля не дивлячись, на те, що плоскі хвилі це стаціонарний стан і отже

{\frac {d|\Psi |^{2}}{dt}}=0

всюди. Це демонструє, що частинка може рухатись, навіть якщо її площинна щільність ймовірності не має ніякої явної залежності від часу.

Частинка в ящику

Для одновимірного ящика з нескінченним стінками довжиною $L$ ( $0<x<L$ ), хвильові функції запишуться у вигляді

\Psi _{n}={\sqrt {\frac {2}{L}}}\sin \left({\frac {n\pi }{L}}x\right)

та нуль справа і зліва від ями. Тоді струм запишеться у вигляді

j_{n}={\frac {\hbar }{2mi}}\left(\Psi _{n}^{*}{\frac {\partial \Psi _{n}}{\partial x}}-\Psi _{n}{\frac {\partial \Psi _{n}^{*}}{\partial x}}\right)=0

оскільки $\Psi _{n}=\Psi _{n}^{*}.$

Виведення рівняння неперевності

В цьому розділі рівняння неперевності виводиться із визначення струму ймовірності та основних принципів квантової механіки

Припустимо, що $\Psi$ — хвильова функція, яка залежить від трьох змінних $x$ , $y$ , та $z$ ). Тоді

P=\int \limits _{V}|\Psi |^{2}dV

визначає ймовірність виміряти позицію частинки в об'ємі V. Похідна по часу запишеться у вигляді

{\frac {dP}{dt}}={\frac {\partial }{\partial t}}\int \limits _{V}|\Psi |^{2}dV=\int \limits _{V}\left({\frac {\partial \Psi }{\partial t}}\Psi ^{*}+\Psi {\frac {\partial \Psi ^{*}}{\partial t}}\right)dV

де останнє рівння припускає, що часткову похідну по часу можна внести під інтеграл (форма об'єму $V$ не залежить від часу). Для подальшого спрощення роглянемо нестаціонарне Рівняння Шредінгера

i\hbar {\frac {\partial \Psi }{\partial t}}={\frac {-\hbar ^{2}}{2m}}\nabla ^{2}\Psi +V\Psi

і використаємо його, щоб виділити похідну по часу від $\Psi$ :

{\frac {\partial \Psi }{\partial t}}={\frac {i\hbar }{2m}}\nabla ^{2}\Psi -{\frac {i}{\hbar }}V\Psi

Реузльат підстановки в попереднє рівняння для ${\frac {dP}{dt}}$ дає

{\frac {dP}{dt}}=-\int \limits _{V}{\frac {\hbar }{2mi}}\left(\Psi ^{*}\nabla ^{2}\Psi -\Psi \nabla ^{2}\Psi ^{*}\right)dV

.

Тепер після переходу до дивергенції

\nabla \cdot \left(\Psi ^{*}{\vec {\nabla }}\Psi -\Psi {\vec {\nabla }}\Psi ^{*}\right)={\vec {\nabla }}\Psi ^{*}\cdot {\vec {\nabla }}\Psi +\Psi ^{*}\nabla ^{2}\Psi -{\vec {\nabla }}\Psi \cdot {\vec {\nabla }}\Psi ^{*}-\Psi {\vec {\nabla }}^{2}\Psi ^{*}

і, оскільки, перший та третій доданки скорочуються:

{\frac {dP}{dt}}=-\int \limits _{V}{\vec {\nabla }}\cdot {\frac {\hbar }{2mi}}\left(\Psi ^{*}{\vec {\nabla }}\Psi -\Psi {\vec {\nabla }}\Psi ^{*}\right)dV

Якщо згадаємо вираз для $P$ і зауважимо, що вираз на який діє оператор набла є ${\vec {j}}$ тоді запишем вираз

\int \limits _{V}\left({\frac {\partial |\Psi |^{2}}{\partial t}}+{\vec {\nabla }}\cdot {\vec {j}}\right)dV=0

який є інтегральною формою рівняння неперевності. Диференціальна форма випливає з того факту, що попередні рівняння виконана для всіх об'ємів $V$ , і інтегралом можна знехтувати

{\frac {\partial |\Psi |^{2}}{\partial t}}+{\vec {\nabla }}\cdot {\vec {j}}=0.

Посилання

Савула Я. Метод скінченних елементів (окремі сторінки посібника 1993 р.) http://old.ami.lnu.edu.ua/books/AMI/savula.pdf

Шинкаренко Г. Чисельні методи математичної фізики (окремі сторінки чорновика посібника)http://old.ami.lnu.edu.ua/books/AMI/nmmf.pdf

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.