I am trying to find out relativistic description of a quantum harmonic oscillator. For a classical relativistic oscillator mass is a function of co-ordinates(http://arxiv.org/abs/1209.2876). $$m(x)=m-\frac{m\omega^{2}x^{2}}{c^{2}}.$$ I want to get a quantum mechanical model for this particular hamiltonian $$H=\sqrt{p^{2}c^{2}+m(x)^{2}c^{4}}.$$ First, I thought of re-writing this hamiltonian with the help of $\alpha$ and $\beta$ matrices. I got $$H=\vec{\alpha}.\vec{p}+\beta m(x).$$ But by adopting this prescription implies wave function should be $4$-component spinors. How can I get $4$-component spinors which would reduce to $1$-component wave function in non-relativistic quantum mechanics? If someone can suggest some other methods to solve for the hamiltonian that would also help.

- 201,751

- 399
-
Just a terminological comment, not just questioning your wording but the wording in the aforementioned paper as well. "Harmonic" is something that is composed of sines and cosines (the solutions to the motion of the oscillator in this case). If the solutions aren't sines and cosines, e.g. because the energy isn't simply quadratic, one shouldn't call it "harmonic" oscillator. It's a relativistic generalization of the harmonic oscillator. – Luboš Motl Apr 22 '13 at 16:53
-
You are right.Changed it. – WInterfell Apr 22 '13 at 16:57
-
Your hamiltonian should have $m^2 c^4$. Also I would suggest that you can try to find (and I believe you will succeed) the eigenvalues and eigenstates of the operator under the square root. On the first thought, $H=\sqrt{Q}$ and $Q$ should have the same eigenvectors and correspondingly related eigenvalues. Note the reasoning behind the Klein-Gordon equation. You should also keep in mind that you should be carefull with physical interpretation of your results -- relativistic QM is not totally physically consitent and should be replaced with QFT. – Peter Kravchuk Apr 22 '13 at 19:33
-
Moreover, I think that your hamiltonian in $\alpha$ and $\beta$ matrices is not obviously the correct square root -- $m$ and $p$ do not commute. Finally, these matrices are some non-commuting operators, not initially present in the theory -- you are introducing new degrees of freedom, specifically spin. Just thought that it can be not exactly what you want. – Peter Kravchuk Apr 22 '13 at 19:42
2 Answers
One could include harmonic potential in Dirac equation in the usual way (see, e.g., here): $$ \left(\gamma^\mu\left(i\hbar\partial_\mu - eA_\mu\right)-mc\right)\psi=0 $$ or in more mundane notation (and stationary in time): $$ \begin{pmatrix} mc^2-E+e\phi & c\mathbf{\sigma}\cdot(\mathbf{p}-e\mathbf{A})\\ -c\mathbf{\sigma}\cdot(\mathbf{p}-e\mathbf{A}) & mc^2+E-e\phi \end{pmatrix} \begin{pmatrix} \psi_+\\ \psi_- \end{pmatrix} = \begin{pmatrix} 0\\ 0 \end{pmatrix} $$ where $$e\phi = \frac{kx^2}{2}$$
Remark
Although relativistically one might be justified interpreting the potential energy as mass, one should not be misled by the usual notation for the spring constant $k=m\Omega^2$, which is just a matter of convenience.

- 58,522
The answer is simple, you need to use the Hamiltonian $$H=mc^2(u_0^2-\sum_{k=1}^3u_k^2)+\frac{k(x_0^2-\sum_{k=1}^3 x_k^2)}{2}$$
Then we get 4 independent standard problems about the harmonic oscillator

- 1
-
The formula $m(x)=m-\frac{m \omega^2 x^2}{c^2}$ is not relativistically invariant and therefore not correct, the correct formula is the dependence of energy on velocity, which is done in my formulas. – Evgeniy Jan 15 '23 at 13:39