Ladder operator

In linear algebra (and its application to quantum mechanics), a raising or lowering operator (collectively known as ladder operators) is an operator that increases or decreases the eigenvalue of another operator. In quantum mechanics, the raising operator is sometimes called the creation operator, and the lowering operator the annihilation operator. Well-known applications of ladder operators in quantum mechanics are in the formalisms of the quantum harmonic oscillator and angular momentum.

Terminology

Main article: Creation and annihilation operators

There is some confusion regarding the relationship between the raising and lowering ladder operators and the creation and annihilation operators commonly used in quantum field theory. The creation operator a_i^† increments the number of particles in state i, while the corresponding annihilation operator a_i decrements the number of particles in state i. This clearly satisfies the requirements of the above definition of a ladder operator: the incrementing or decrementing of the eigenvalue of another operator (in this case the particle number operator).

Confusion arises because the term ladder operator is typically used to describe an operator that acts to increment or decrement a quantum number describing the state of a system. To change the state of a particle with the creation/annihilation operators of QFT requires the use of both annihilation and creation operators. An annihilation operator is used to remove a particle from the initial state and a creation operator is used to add a particle to the final state.

The term "ladder operator" is also sometimes used in mathematics, in the context of the theory of Lie algebras and in particular the affine Lie algebras, to describe the su(2) subalgebras, from which the root system and the highest weight modules can be constructed by means of the ladder operators.^[1] In particular, the highest weight is annihilated by the raising operators; the rest of the positive root space is obtained by repeatedly applying the lowering operators (one set of ladder operators per subalgebra).

General formulation

Suppose that two operators X and N have the commutation relation,

$${\displaystyle [N,X]=cX,}$$

 

for some scalar c. If ${\displaystyle {|n\rangle }}$  is an eigenstate of N with eigenvalue equation,

$${\displaystyle N|n\rangle =n|n\rangle ,}$$

 

then the operator X acts on ${\displaystyle {|n\rangle }}$  in such a way as to shift the eigenvalue by c:

$${\displaystyle {\begin{aligned}NX|n\rangle &=(XN+[N,X])|n\rangle \\&=XN|n\rangle +[N,X]|n\rangle \\&=Xn|n\rangle +cX|n\rangle \\&=(n+c)X|n\rangle .\end{aligned}}}$$

 

In other words, if ${\displaystyle {|n\rangle }}$  is an eigenstate of N with eigenvalue n then ${\displaystyle {X|n\rangle }}$  is an eigenstate of N with eigenvalue n + c or it is zero. The operator X is a raising operator for N if c is real and positive, and a lowering operator for N if c is real and negative.

If N is a Hermitian operator then c must be real and the Hermitian adjoint of X obeys the commutation relation:

$${\displaystyle [N,X^{\dagger }]=-cX^{\dagger }.}$$

 

In particular, if X is a lowering operator for N then X^† is a raising operator for N and vice versa.

Angular momentum

Main article: Angular momentum operator

A particular application of the ladder operator concept is found in the quantum mechanical treatment of angular momentum. For a general angular momentum vector, J, with components, J_x, J_y and J_z one defines the two ladder operators, J₊ and J_–,^[2]

$${\displaystyle J_{+}=J_{x}+iJ_{y},}$$

 

$${\displaystyle J_{-}=J_{x}-iJ_{y},}$$

 

where i is the imaginary unit.

The commutation relation between the cartesian components of any angular momentum operator is given by

$${\displaystyle [J_{i},J_{j}]=i\hbar \epsilon _{ijk}J_{k},}$$

 

where ε_ijk is the Levi-Civita symbol and each of i, j and k can take any of the values x, y and z.

From this, the commutation relations among the ladder operators and J_z are obtained,

$${\displaystyle \left[J_{z},J_{\pm }\right]=\pm \hbar J_{\pm },}$$

 

$${\displaystyle \left[J_{+},J_{-}\right]=2\hbar J_{z}.}$$

 

(Technically, this is the Lie algebra of ${\displaystyle {{\mathfrak {s}}l}(2,\mathbb {R} )}$ ).

The properties of the ladder operators can be determined by observing how they modify the action of the J_z operator on a given state,

$${\displaystyle {\begin{aligned}J_{z}J_{\pm }|j\,m\rangle &=\left(J_{\pm }J_{z}+\left[J_{z},J_{\pm }\right]\right)|j\,m\rangle \\&=\left(J_{\pm }J_{z}\pm \hbar J_{\pm }\right)|j\,m\rangle \\&=\hbar \left(m\pm 1\right)J_{\pm }|j\,m\rangle .\end{aligned}}}$$

 

Compare this result with

$${\displaystyle J_{z}|j\,(m\pm 1)\rangle =\hbar (m\pm 1)|j\,(m\pm 1)\rangle .}$$

 

Thus one concludes that ${\displaystyle {J_{\pm }|j\,m\rangle }}$  is some scalar multiplied by ${\displaystyle {|j\,(m\pm 1)\rangle }}$ ,

$${\displaystyle J_{+}|j\,m\rangle =\alpha |j\,(m+1)\rangle ,}$$

 

$${\displaystyle J_{-}|j\,m\rangle =\beta |j\,(m-1)\rangle .}$$

 

This illustrates the defining feature of ladder operators in quantum mechanics: the incrementing (or decrementing) of a quantum number, thus mapping one quantum state onto another. This is the reason that they are often known as raising and lowering operators.

To obtain the values of α and β first take the norm of each operator, recognizing that J₊ and J₋ are a Hermitian conjugate pair (${\displaystyle {J_{\pm }=J_{\mp }^{\dagger }}}$ ),

$${\displaystyle \langle j\,m|J_{+}^{\dagger }J_{+}|j\,m\rangle =\langle j\,m|J_{-}J_{+}|j\,m\rangle =\langle j\,(m+1)|\alpha ^{*}\alpha |j\,(m+1)\rangle =|\alpha |^{2},}$$

 

$${\displaystyle \langle j\,m|J_{-}^{\dagger }J_{-}|j\,m\rangle =\langle j\,m|J_{+}J_{-}|j\,m\rangle =\langle j\,(m-1)|\beta ^{*}\beta |j\,(m-1)\rangle =|\beta |^{2}.}$$

 

The product of the ladder operators can be expressed in terms of the commuting pair J² and J_z,

$${\displaystyle J_{-}J_{+}=(J_{x}-iJ_{y})(J_{x}+iJ_{y})=J_{x}^{2}+J_{y}^{2}+i[J_{x},J_{y}]=J^{2}-J_{z}^{2}-\hbar J_{z},}$$

 

$${\displaystyle J_{+}J_{-}=(J_{x}+iJ_{y})(J_{x}-iJ_{y})=J_{x}^{2}+J_{y}^{2}-i[J_{x},J_{y}]=J^{2}-J_{z}^{2}+\hbar J_{z}.}$$

 

Thus, one may express the values of |α|² and |β|² in terms of the eigenvalues of J² and J_z,

$${\displaystyle |\alpha |^{2}=\hbar ^{2}j(j+1)-\hbar ^{2}m^{2}-\hbar ^{2}m=\hbar ^{2}(j-m)(j+m+1),}$$

 

$${\displaystyle |\beta |^{2}=\hbar ^{2}j(j+1)-\hbar ^{2}m^{2}+\hbar ^{2}m=\hbar ^{2}(j+m)(j-m+1).}$$

 

The phases of α and β are not physically significant, thus they can be chosen to be positive and real (Condon-Shortley phase convention). We then have:^[3]

$${\displaystyle J_{+}|j,m\rangle =\hbar {\sqrt {(j-m)(j+m+1)}}|j,m+1\rangle =\hbar {\sqrt {j(j+1)-m(m+1)}}|j,m+1\rangle ,}$$

 

$${\displaystyle J_{-}|j,m\rangle =\hbar {\sqrt {(j+m)(j-m+1)}}|j,m-1\rangle =\hbar {\sqrt {j(j+1)-m(m-1)}}|j,m-1\rangle .}$$

 

Confirming that m is bounded by the value of j (${\displaystyle {-j\leq m\leq j}}$ ), one has

$${\displaystyle J_{+}|j\,j\rangle =0,}$$

 

$${\displaystyle J_{-}|j\,(-j)\rangle =0.}$$

 

The above demonstration is effectively the construction of the Clebsch-Gordan coefficients.

Applications in atomic and molecular physics

Many terms in the Hamiltonians of atomic or molecular systems involve the scalar product of angular momentum operators. An example is the magnetic dipole term in the hyperfine Hamiltonian,^[4]

$${\displaystyle {\hat {H}}_{\text{D}}={\hat {A}}\mathbf {I} \cdot \mathbf {J} ,}$$

 

where I is the nuclear spin.

The angular momentum algebra can often be simplified by recasting it in the spherical basis. Using the notation of spherical tensor operators, the "-1", "0" and "+1" components of J⁽¹⁾ ≡ J are given by,^[5]

$${\displaystyle {\begin{aligned}J_{-1}^{(1)}&={\dfrac {1}{\sqrt {2}}}(J_{x}-iJ_{y})={\dfrac {J_{-}}{\sqrt {2}}}\\J_{0}^{(1)}&=J_{z}\\J_{+1}^{(1)}&=-{\frac {1}{\sqrt {2}}}(J_{x}+iJ_{y})=-{\frac {J_{+}}{\sqrt {2}}}.\end{aligned}}}$$

 

From these definitions, it can be shown that the above scalar product can be expanded as

$${\displaystyle \mathbf {I} ^{(1)}\cdot \mathbf {J} ^{(1)}=\sum _{n=-1}^{+1}(-1)^{n}I_{n}^{(1)}J_{-n}^{(1)}=I_{0}^{(1)}J_{0}^{(1)}-I_{-1}^{(1)}J_{+1}^{(1)}-I_{+1}^{(1)}J_{-1}^{(1)},}$$

 

The significance of this expansion is that it clearly indicates which states are coupled by this term in the Hamiltonian, that is those with quantum numbers differing by m_i = ±1 and m_j = ∓1 only.

Harmonic oscillator

Main article: Quantum harmonic oscillator

Another application of the ladder operator concept is found in the quantum mechanical treatment of the harmonic oscillator. We can define the lowering and raising operators as

$${\displaystyle {\begin{aligned}{\hat {a}}&={\sqrt {m\omega \over 2\hbar }}\left({\hat {x}}+{i \over m\omega }{\hat {p}}\right)\\{\hat {a}}^{\dagger }&={\sqrt {m\omega \over 2\hbar }}\left({\hat {x}}-{i \over m\omega }{\hat {p}}\right)\end{aligned}}}$$

 

They provide a convenient means to extract energy eigenvalues without directly solving the system's differential equation.

Hydrogen-like atom

Main article: Hydrogen-like atom

There are two main approaches given in the literature using ladder operators, one using the Laplace–Runge–Lenz vector, another using factorization of the Hamiltonian.

Laplace–Runge–Lenz vector

Another application of the ladder operator concept is found in the quantum mechanical treatment of the electronic energy of hydrogen-like atoms and ions. The Laplace–Runge–Lenz vector commutes with the Hamiltonian for an inverse square spherically symmetric potential and can be used to determine ladder operators for this potential.^[6][7] We can define the lowering and raising operators (based on the classical Laplace–Runge–Lenz vector)

$${\displaystyle {\vec {A}}=\left({\frac {1}{Ze^{2}\mu }}\right)\left\{{\vec {L}}\times {\vec {p}}-{\boldsymbol {i}}\hbar {\vec {p}}\right\}+{\frac {\vec {r}}{r}}}$$

 

where ${\displaystyle {\vec {L}}}$  is the angular momentum, ${\displaystyle {\vec {p}}}$  is the linear momentum, ${\displaystyle \mu }$  is the reduced mass of the system, ${\displaystyle e}$  is the electronic charge, and ${\displaystyle Z}$  is the atomic number of the nucleus. Analogous to the angular momentum ladder operators, one has ${\displaystyle A_{+}=A_{x}+iA_{y}}$  and ${\displaystyle A_{-}=A_{x}-iA_{y}}$ .

The commutators needed to proceed are:

$${\displaystyle [A_{\pm },L_{z}]=\mp {\boldsymbol {i}}\hbar A_{\mp }}$$

 

and

$${\displaystyle [A_{\pm },L^{2}]=\mp 2\hbar ^{2}A_{\pm }-2\hbar A_{\pm }L_{z}\pm 2\hbar A_{z}L_{\pm }}$$

 

. Therefore,

$${\displaystyle A_{+}|?,\ell ,m_{\ell }\rangle \rightarrow |?,\ell ,m_{\ell }+1\rangle }$$

 

and

$${\displaystyle -L^{2}\left(A_{+}|?,\ell ,\ell \rangle \right)=-\hbar ^{2}(\ell +1)((\ell +1)+1)\left(A_{+}|?,\ell ,\ell \rangle \right)}$$

 

so

$${\displaystyle A_{+}|?,\ell ,\ell \rangle \rightarrow |?,\ell +1,\ell +1\rangle }$$

 

where the "?" indicates a nascent quantum number which emerges from the discussion.

Given the Pauli^[8][9] equations Pauli Equation IV:

$${\displaystyle 1-A\cdot A=-\left({\frac {2E}{\mu Z^{2}e^{4}}}\right)(L^{2}+\hbar ^{2})}$$

 

and Pauli Equation III:

$${\displaystyle \left(A\times A\right)_{j}=-\left({\frac {2{\boldsymbol {i}}\hbar E}{\mu Z^{2}e^{4}}}\right)L_{j}}$$

 

and starting with the equation

$${\displaystyle A_{-}A_{+}|\ell ^{*},\ell ^{*}\rangle =0}$$

 

and expanding, one obtains (assuming ${\displaystyle \ell ^{*}}$  is the maximum value of the angular momentum quantum number consonant with all other conditions),

$${\displaystyle \left(1+{\frac {2E}{\mu Z^{2}e^{4}}}(L^{2}+\hbar ^{2})-i{\frac {2i\hbar E}{\mu Z^{2}e^{4}}}L_{z}\right)|?,\ell ^{*},\ell ^{*}\rangle =0}$$

 

which leads to the Rydberg formula:

$${\displaystyle E_{n}=-{\frac {\mu Z^{2}e^{4}}{2\hbar ^{2}(\ell ^{*}+1)^{2}}}}$$

 

implying that ${\displaystyle \ell ^{*}+1=n=?}$ , where ${\displaystyle n}$  is the traditional quantum number.

Factorization of the Hamiltonian

The Hamiltonian for a hydrogen-like potential can be written in spherical coordinates as

$${\displaystyle H=1/(2\mu )[p_{r}^{2}+(1/r^{2})L^{2}]+V(r)}$$

 

where ${\displaystyle V(r)=-Ze^{2}/r}$  and ${\displaystyle p_{r}}$  is the radial momentum

$${\displaystyle p_{r}=(x/r)p_{x}+(y/r)p_{y}+(z/r)p_{z}}$$

 

which is real and self-conjugate.

Suppose ${\displaystyle |nl\rangle }$  is an eigenvector of the Hamiltonian where ${\displaystyle l\rangle }$  is the angular momentum and ${\displaystyle n}$  represents the energy, so ${\displaystyle L^{2}|nl\rangle =l(l+1)\hbar ^{2}|nl\rangle }$  and we may label the Hamiltonian as ${\displaystyle H_{l}}$ 

$${\displaystyle H=1/(2\mu )[p_{r}^{2}+(1/r^{2})l(l+1)\hbar ^{2}]+V(r)}$$

 

The factorization method was developed by Infeld and Hull^[10] for differential equations. Newmarch and Golding^[11] applied it to spherically symmetric potentials using operator notation.

Suppose we can find a factorization of the Hamiltonian by operators ${\displaystyle C_{l}}$  as

${\displaystyle C_{l}^{*}C_{l}=2\mu H_{l}+F_{l}}$

(1)

and

$${\displaystyle C_{l}C_{l}^{*}=2\mu H_{l+1}+G_{l}}$$

 

for scalars ${\displaystyle F_{l}}$  and ${\displaystyle G_{l}}$ . The vector ${\displaystyle C_{l}C_{l}^{*}C_{l}|nl\rangle }$  may be evaluated in two different ways as

$${\displaystyle {\begin{aligned}C_{l}C_{l}^{*}C_{l}|nl\rangle &=(2\mu E_{l}^{n}+F_{l})C_{l}|nl\rangle \\&=(2\mu H_{l+1}+G_{l})C_{l}|nl\rangle \end{aligned}}}$$

 

which can be re-arranged as

$${\displaystyle H_{l+1}(C_{l}|nl\rangle )=[E_{l}^{n}+(F_{l}-G_{l})/(2\mu )](C_{l}|nl\rangle )}$$

 

showing that ${\displaystyle C_{l}|nl\rangle }$  is an eigenstate of ${\displaystyle H_{l+1}}$  with eigenvalue

$${\displaystyle E_{l+1}^{n^{'}}=E_{l}^{n}+(F_{l}-G_{l})/(2\mu )}$$

 

If ${\displaystyle F_{l}=G_{l}}$  then ${\displaystyle n^{'}=n}$  and the states ${\displaystyle |nl\rangle }$  and ${\displaystyle C_{l}|nl\rangle }$  have the same energy.

For the hydrogenic atom, setting

$${\displaystyle V(r)=-{\frac {B\hbar }{\mu r}}}$$

 

with

$${\displaystyle B={\frac {Z\mu e^{2}}{\hbar }}}$$

 

a suitable equation for ${\displaystyle C_{l}}$  is

$${\displaystyle C_{l}=p_{r}+{\frac {i\hbar (l+1)}{r}}-{\frac {iB}{l+1}}}$$

 

with

$${\displaystyle F_{l}=G_{l}={\frac {B^{2}}{(l+1)^{2}}}}$$

 

There is an upper bound to the ladder operator if the energy is negative, (so ${\displaystyle C_{l}|nl_{\max }\rangle =0}$  for some ${\displaystyle l_{\max }}$ ) then from Equation (1)

$${\displaystyle E_{l}^{n}=-F_{l}/{2\mu }=-{\frac {B^{2}}{2\mu (l_{\max }+1)^{2}}}=-{\frac {\mu Z^{2}e^{4}}{2\hbar ^{2}(l_{\max }+1)^{2}}}}$$

 

and ${\displaystyle n}$  can be identified with ${\displaystyle l_{\max }+1}$ 

Relation to group theory

Whenever there is degeneracy in a system, there is usually a related symmetry property and group. The degeneracy of the energy levels for the same value of ${\displaystyle n}$  but different angular momenta has been identified as the SO(4) symmetry of the spherically symmetric Coulomb potential.^[12][13]

3D isotropic harmonic oscillator

The 3D isotropic harmonic oscillator has a potential given by

$${\displaystyle V(r)=\mu \omega ^{2}r^{2}/2}$$

 

It can similarly be managed using the factorization method.

Factorization method

A suitable factorization is given by^[11]

$${\displaystyle C_{l}=p_{r}+{\frac {i\hbar (l+1)}{r}}-i\mu \omega r}$$

 

with

$${\displaystyle F_{l}=-(2l+3)\mu \omega \hbar }$$

 

and

$${\displaystyle G_{l}=-(2l+1)\mu \omega \hbar }$$

 

Then

$${\displaystyle E_{l+1}^{n^{'}}=E_{l}^{n}+{\frac {F_{l}-G_{l}}{2\mu }}=E_{l}^{n}-\omega \hbar }$$

 

and continuing this,

$${\displaystyle {\begin{aligned}E_{l+2}^{n^{'}}&=E_{l}^{n}-2\omega \hbar \\E_{l+3}^{n^{'}}&=E_{l}^{n}-3\omega \hbar \\\dots &\end{aligned}}}$$

 

Now the Hamiltonian only has positive energy levels as can be seen from

$${\displaystyle {\begin{aligned}\langle \psi |2\mu H_{l}|\psi \rangle &=\langle \psi |C_{l}^{*}C_{l}|\psi \rangle +\langle \psi |(2l+3)\mu \omega \hbar |\psi \rangle \\&=\langle C_{l}\psi |C_{l}\psi \rangle +(2l+3)\mu \omega \hbar \langle \psi |\psi \rangle \\&\geq 0\end{aligned}}}$$

 

This means that for some value of ${\displaystyle l}$  the series must terminate with ${\displaystyle C_{l_{\max }}|nl_{\max }\rangle =0}$  and then

$${\displaystyle E_{l_{\max }}^{n}=-F_{l_{\max }}/(2\mu )=(l_{\max }+3/2)\omega \hbar }$$

 

This is decreasing in energy by ${\displaystyle \omega \hbar }$  unless for some value of ${\displaystyle l}$ , ${\displaystyle C_{l}|nl\rangle =0}$ . Identifying this value as ${\displaystyle n}$  gives

$${\displaystyle E_{l}^{n}=-F_{l}=(n+3/2)\omega \hbar }$$

 

It then follows the ${\displaystyle n'=n-1}$  so that

$${\displaystyle C_{l}|nl\rangle =\lambda _{l}^{n}|n-1\,l+1\rangle }$$

 

giving a recursion relation on ${\displaystyle \lambda }$  with solution

$${\displaystyle \lambda _{l}^{n}=-\mu \omega \hbar {\sqrt {2(n-l)}}}$$

 

There is degeneracy caused from angular momentum; there is additional degeneracy caused by the oscillator potential. Consider the states ${\displaystyle |n\,n\rangle ,|n-1\,n-1\rangle ,|n-2\,n-2\rangle ,\dots }$  and apply the lowering operators ${\displaystyle C^{*}}$ : ${\displaystyle C_{n-2}^{*}|n-1\,n-1\rangle ,C_{n-4}^{*}C_{n-3}^{*}|n-2\,n-2\rangle ,\dots }$  giving the sequence ${\displaystyle |nn\rangle ,|n\,n-2\rangle ,|n\,n-4\rangle ,...}$  with the same energy but with ${\displaystyle l}$  decreasing by 2. In addition to the angular momentum degeneracy, this gives a total degeneracy of ${\displaystyle (n+1)(n+2)/2}$ ^[14]

Relation to group theory

The degeneracies of the 3D isotropic harmonic oscillator are related to the special unitary group SU(3)^[14][15]

History

Many sources credit Dirac with the invention of ladder operators.^[16] Dirac's use of the ladder operators shows that the total angular momentum quantum number ${\displaystyle j}$  needs to be a non-negative half integer multiple of ħ.

See also

• Creation and annihilation operators
• Quantum harmonic oscillator
• Chevalley basis

References

1. Fuchs, Jurgen (1992), Affine Lie Algebras and Quantum Groups, Cambridge University Press, ISBN 0-521-48412-X
2. de Lange, O. L.; R. E. Raab (1986). "Ladder operators for orbital angular momentum". American Journal of Physics. 54 (4): 372–375. Bibcode:1986AmJPh..54..372D. doi:10.1119/1.14625.
3. Sakurai, Jun J. (1994). Modern Quantum Mechanics. Delhi, India: Pearson Education, Inc. p. 192. ISBN 81-7808-006-0.
4. Woodgate, Gordon K. (1983-10-06). Elementary Atomic Structure. ISBN 978-0-19-851156-4. Retrieved 2009-03-03.
5. "Angular Momentum Operators". Graduate Quantum Mechanics Notes. University of Virginia. Retrieved 2009-04-06.
6. David, C. W. (1966). "Ladder Operator Solution for the Hydrogen Atom Electronic Energy Levels". American Journal of Physics. 34 (10): 984–985. Bibcode:1966AmJPh..34..984D. doi:10.1119/1.1972354.
7. Burkhardt, C. E.; Levanthal, J. (2004). "Lenz vector operations on spherical hydrogen atom eigenfunctions". American Journal of Physics. 72 (8): 1013–1016. Bibcode:2004AmJPh..72.1013B. doi:10.1119/1.1758225.
8. Pauli, Wolfgang (1926). "Über das Wasserstoffspektrum vom Standpunkt der neuen Quantenmechanik". Z. Phys. 36 (5): 336–363. Bibcode:1926ZPhy...36..336P. doi:10.1007/BF01450175. S2CID 128132824.
9. B. L. Van der Waerden, Sources of Quantum Mechanics, Dover, New York, 1968
10. L., Infeld; Hull, T.E. (1951). "The Factorization Method". Rev. Mod. Phys. 23 (1): 21–68. Bibcode:1951RvMP...23...21I. doi:10.1103/RevModPhys.23.21.
11. Newmarch, J. D.; Golding, R. M. (1978). "Ladder operators for some spherically symmetric potentials in quantum". Am. J. Phys. 46: 658–660. doi:10.1119/1.11225.
12. Weinberg, S. J. (2011). "The SO(4) Symmetry of the Hydrogen Atom" (PDF). {{cite journal}}: Cite journal requires |journal= (help)
13. Lahiri, A.; Roy, P. K.; Bagchi, B. (1989). "Supersymmetry and the Ladder Operator Technique in Quantum Mechanics: The Radial Schrödinger Equation". Int. J. Theor. Phys. 28 (2): 183–189. Bibcode:1989IJTP...28..183L. doi:10.1007/BF00669809. S2CID 123255435.
14. Kirson, M. W. (2013). "Introductory Algebra for Physicists: Isotropic harmonic oscillator" (PDF). Weizmann Institute of Science. Retrieved 28 July 2021.
15. , D. M. "." 33 (3) (1965) 207–211. Fradkin, D. M. (1965). "Three-dimensional isotropic harmonic oscillator and SU3". Am. J. Phys. 33 (3): 207–211. Bibcode:1965AmJPh..33..207F. doi:10.1119/1.1971373.
16. https://www.fisica.net/mecanica-quantica/quantum_harmonic_oscillator_lecture.pdf^{[bare URL PDF]}