1.73: Time‐dependent Wigner Function for Harmonic Oscillator Transitions
- Page ID
- 156463
Initial state:
\[
\mathrm{m} :=0 \qquad \mathrm{E}_{\mathrm{m}} :=\mathrm{m}+\frac{1}{2}
\nonumber \]
Final state:
\[
\mathrm{n} :=1 \qquad \mathrm{E}_{\mathrm{n}} :=\mathrm{n}+\frac{1}{2} \qquad t : = \text{FRAME}
\nonumber \]
Define Wigner distribution function for a linear superposition of the initial and final harmonic oscillator state.
\[
W(x,p) : = \frac{1}{\pi^{\frac{3}{2}}}\int_{-\infty}^{\infty}\left[\frac{1}{\sqrt{2^{\mathrm{n}} \cdot \mathrm{n} ! \sqrt{\pi}}} \cdot \operatorname{Her}\left(\mathrm{n}, x+\frac{s}{2}\right) \cdot \exp \left[-\frac{\left(\mathrm{x}+\frac{\mathrm{s}}{2}\right)^{2}}{2}\right] \cdot \exp \left(\mathrm{i} \cdot \mathrm{E}_{\mathrm{n}} \cdot \mathrm{t}\right)+\frac{1}{\sqrt{2^{\mathrm{m}} \cdot \mathrm{m} ! \cdot \sqrt{\pi}}} \cdot \operatorname{Her}\left(\mathrm{m}, \mathrm{x}+\frac{\mathrm{s}}{2}\right) \cdot \exp \left[-\frac{\left(\mathrm{x}+\frac{\mathrm{s}}{2}\right)^{2}}{2}\right] \cdot \exp \left(\mathrm{i} \cdot \mathrm{E}_{\mathrm{m}} \cdot \mathrm{t}\right)\right] \\ \cdot \exp (\mathrm{i} \cdot \mathrm{s} \cdot \mathrm{p}) \cdot \left[\frac{1}{\sqrt{2^{\mathrm{n}} \cdot n ! \cdot \sqrt{\pi}}} \cdot \operatorname{Her}\left(\mathrm{n}, \mathrm{x}-\frac{\mathrm{s}}{2}\right) \cdot \exp \left[-\frac{\left(\mathrm{x}-\frac{\mathrm{s}}{2}\right)^{2}}{2}\right] \cdot \exp \left(-\mathrm{i} \cdot \mathrm{E}_{\mathrm{n}} \cdot \mathrm{t}\right)\\+\frac{1}{\sqrt{2^{\mathrm{m}} \cdot \mathrm{m} ! \cdot \sqrt{\pi}}} \cdot \operatorname{Her}\left(\mathrm{m}, \mathrm{x}-\frac{\mathrm{s}}{2}\right)\cdot\exp\left[-\frac{\left(\mathrm{x}-\frac{\mathrm{s}}{2}\right)^{2}}{2}\right]\cdot \exp\left(\mathrm{-i}\cdot\mathrm{E}_{\mathrm{m}}\cdot\mathrm{t}\right)\right]\mathrm{ds}
\nonumber \]
Display Wigner distribution:
\[
\mathrm{N}=60 \qquad \mathrm{i} :=0 \ldots \mathrm{N} \qquad \mathrm{x}_{\mathrm{i}} :=-2.5+\frac{5 \cdot \mathrm{i}}{\mathrm{N}} \\ \mathrm{j} :=0 \ldots \mathrm{N} \qquad \mathrm{P}_{\mathrm{j}} :=-2.5+\frac{5 \mathrm{j}}{\mathrm{N}} \qquad \text{Wigner}_{i,j}: = W(x_{i},p_{j})
\nonumber \]