Metoda podziału operatora – w mechanice kwantowej podstawowa metoda numeryczna rozwiązywania równania Schrödingera zależnego od czasu. Polega na zastąpieniu potencjału w równaniu Schrödingera jego postacią spróbkowaną[1].
Niech równanie Schrödingera zależne od czasu będzie dane w ogólnej postaci
![{\displaystyle {\hat {H}}{\big |}\psi (t)\rangle =i{\tfrac {\partial }{\partial t}}{\big |}\psi (t)\rangle ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/66b127522912d55f3f3841980d43104f5aa06820)
gdzie:
![{\displaystyle {\hat {H}}={\hat {T}}+{\hat {V}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/caf5bc48dfd61ffde8c2a3c0d2d49ac58b2712e4)
z ogólnym rozwiązaniem
![{\displaystyle |\psi (t)\rangle =e^{-i{{\int \limits _{0}^{t}H(\tau )}d\tau }}|\psi (0)\rangle .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce519692f357f48db60924e1c5656f2df6e4ac3e)
Wtedy zastępujemy potencjał potencjałem spróbkowanym
![{\displaystyle V(t)\approx \sum _{n=-\infty }^{+\infty }{\delta (t-n\Delta t)}V(t)\Delta t.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47dc2b6a0b43be0476254bef7110f106097b799d)
Umożliwia to scałkowanie równania dokładnie w każdym przedziale pomiędzy funkcjami
tzn.
![{\displaystyle |\psi (t)\rangle =e^{-i{\hat {V}}(0)\Delta t/2}e^{-i{\hat {T}}\Delta t}e^{-i{\hat {V}}(\Delta t)\Delta t}...e^{-i{\hat {T}}\Delta t}e^{-i{\hat {V}}(n\Delta t)\Delta t/2}|\psi (0)\rangle .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/096fa6220334e3b399fe7e7e942ac84a183a3e45)
Każdy z kroków w tym wyrażeniu jest łatwy do wykonania numerycznie. Aby jeszcze obniżyć błąd, pomija się energię kinetyczną dla członu zerowego, tzn. w chwili
oraz w zerowym i ostatnim członie używa się połowy potencjału.
Przypisy
- ↑ M.D. Feit, J.A. Jr. Fleck, A. Steiger. Solution of the Schrödinger equation by a spectral method. „J. Comput. Phys.”. 47, s. 412, 1982.