Phy5645/Energy conservation

(1) The energy operator in three dimensions is: $H=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V$ so the average energy in state $\Psi$ is: $\left\langle E\right\rangle =\iiint \Psi ^{\ast }H\Psi \,d^{3}{\textbf {r}}=\iiint \Psi ^{\ast }\left(-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\Psi +V\Psi \right)\,d^{3}{\textbf {r}}$

Using the identity, $\Psi ^{*}\nabla ^{2}\Psi =\nabla \cdot \left(\Psi ^{*}\nabla \Psi \right)-\nabla \Psi ^{\ast }\cdot \nabla \Psi ,$ we obtain

$\left\langle E\right\rangle =-{\frac {\hbar ^{2}}{2m}}\iiint \left(\nabla \cdot \left(\Psi ^{\ast }\nabla \Psi \right)-\nabla \Psi ^{\ast }\cdot \nabla \Psi \right)\,d^{3}{\textbf {r}}+\iiint \Psi ^{\ast }V\Psi \,d^{3}{\textbf {r}}$ $=-{\frac {\hbar ^{2}}{2m}}\iiint \nabla \cdot \left(\Psi ^{\ast }\nabla \Psi \right)\,d^{3}{\textbf {r}}+{\frac {\hbar ^{2}}{2m}}\iiint \nabla \Psi ^{\ast }\cdot \nabla \Psi \,d^{3}{\textbf {r}}+\iiint \Psi ^{\ast }V\Psi \,d^{3}{\textbf {r}}$

If we apply Gauss' Theorem to the first term,

$-{\frac {\hbar ^{2}}{2m}}\iiint \nabla \left(\Psi ^{\ast }\nabla \Psi \right)\,d^{3}{\textbf {r}}=\iint \Psi ^{\ast }\nabla \Psi \cdot d{\textbf {S}},$

as well as the condition, $\lim _{r\to \infty }\Psi ^{\ast }\nabla \Psi =0,$ we obtain

$\left\langle E\right\rangle =\int W\,d^{3}{\textbf {r}}=\int \left[{\frac {\hbar ^{2}}{2m}}\nabla \Psi ^{\ast }\cdot \nabla \Psi +\Psi ^{\ast }V\Psi \right]d^{3}{\textbf {r}}$

(2):first we find the time derivative of energy density:

${\frac {\partial W}{\partial t}}={\frac {\partial }{\partial t}}\left(\nabla \psi ^{*}\nabla \psi +\psi ^{*}\nabla \psi \right)={\frac {\hbar ^{2}}{2m}}\left(\nabla \psi ^{*}\nabla {\frac {\partial \psi }{\partial t}}+\nabla {\frac {\partial \psi ^{*}}{\partial t}}\nabla \psi \right)+{\frac {\partial \psi ^{*}}{\partial t}}\nabla \psi +\psi ^{*}\nabla {\frac {\partial \psi }{\partial t}}$, $={\frac {\hbar ^{2}}{2m}}\left(\nabla \cdot \left(\nabla \psi ^{*}\cdot {\frac {\partial \psi }{\partial t}}+{\frac {\partial \psi ^{*}}{\partial t}}\cdot \nabla \psi \right)-\left({\frac {\partial \psi }{\partial t}}\nabla ^{2}\psi ^{*}+{\frac {\partial \psi ^{*}}{\partial t}}\nabla ^{2}\psi \right)\right)+{\frac {\partial \psi ^{*}}{\partial t}}\nabla \psi +\psi ^{*}\nabla {\frac {\partial \psi }{\partial t}}$ $={\frac {\hbar ^{2}}{2m}}\nabla \cdot \left(\nabla \psi ^{*}\cdot {\frac {\partial \psi }{\partial t}}+{\frac {\partial \psi ^{*}}{\partial t}}\cdot \nabla \psi \right)+{\frac {\partial \psi ^{*}}{\partial t}}\left(-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi +\nabla \psi \right)+{\frac {\partial \psi }{\partial t}}\left(-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi ^{*}+\nabla \psi ^{*}\right)$,

Using Schrodinger Equations: $i\hbar {\frac {\partial \psi }{\partial t}}=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi +\nabla \psi$, and, $-i\hbar {\frac {\partial \psi ^{*}}{\partial t}}=-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi ^{*}+\nabla \psi ^{*}$,

Also the energy flux density is: ${\textbf {S}}=-{\frac {\hbar ^{2}}{2m}}\left({\frac {\partial \psi ^{*}}{\partial t}}\nabla \psi +{\frac {\partial \psi }{\partial t}}\nabla \psi ^{*}\right)$,

So:${\frac {\partial W}{\partial t}}=-\nabla \cdot {\textbf {S}}+{\frac {\partial \psi ^{*}}{\partial t}}{\frac {\partial \psi }{\partial t}}-{\frac {\partial \psi }{\partial t}}{\frac {\partial \psi ^{*}}{\partial t}}=-\nabla \cdot {\textbf {S}}$, Hence: ${\frac {\partial W}{\partial t}}+\nabla \cdot {\textbf {S}}=0$

Back to Relation Between the Wave Function and Probability Density