The magnetic quantum number, designated by the letter \(m_l\), is the third quantum numbers which describe the unique quantum state of an electron. The magnetic quantum number distinguishes the orbitals available within a subshell, and is used to calculate the azimuthal component of the orientation of the orbital in space. As with our discussion of rigid rotors, the quantum number \(m_l\) refers to the projection of the angular momentum in this arbitrarily chosen direction, conventionally called the \(z\) direction or quantization axis. \(L_z\), the magnitude of the angular momentum in the z direction, is given by the formula
The quantum number \(m_l\) refers, loosely, to the direction of the angular momentum vector. The magnetic quantum number \(m_l\) only affects the electron's energy if it is in a magnetic field because in the absence of one, all spherical harmonics corresponding to the different arbitrary values of \(m_l\) are equivalent. The magnetic quantum number determines the energy shift of an atomic orbital due to an external magnetic field (this is called the Zeeman effect) - hence the name magnetic quantum number. However, the actual magnetic dipole moment of an electron in an atomic orbital arrives not only from the electron angular momentum, but also from the electron spin, expressed in the spin quantum number, which is the fourth quantum number. \(m_s\) and discussed in the next chapter.
The answer is complicated; while \(m_l=0\) corresponds to the \(p_z\), the orbitals for \(m_l=+1\) and \(m_l=−1\) lie in the xy-plane (see Spherical Harmonics), but not on the axes. The reason for this outcome is that the wavefunctions are usually formulated in spherical coordinates to make the math easier, but graphs in the Cartesian coordinates make more intuitive sense for humans. The \(p_x\) and \(p_y\) orbitals are constructed via a linear combination approach from radial and angular wavefunctions and converted into \(xy\) (this was discussed previously). Thus, it is not possible to directly correlate the values of \(m_l=±1\) with specific orbitals. The notion that we can do so is sometimes presented in introductory courses to make a complex mathematical model just a little bit simpler and more intuitive, but it is incorrect.
The three wavefunctions for \(n=2\) and \(l=1\) are as follows.
\[ \begin{align} |\psi_{2,1,0} \rangle &=r \cos θR(r) \\[4pt] |\psi_{2,1,+1} \rangle &=−\dfrac{r}{2} \sinθ e^{iϕ} R(r) \\[4pt] |\psi_{2,1,-1} \rangle &=+\dfrac{r}{2} \sinθ e^{-iϕ} R(r) \end{align} \nonumber \]
The notation is \(|\psi_{n,l,m_l} \rangle\) with \(R(r)\) is the radial component of this wavefuction, \(θ\) is the angle with respect to the z-axis and \(ϕ\) is the angle with respect to the \(xz\)-plane.
\[R(r)=\sqrt{\dfrac{Z^5}{32\pi a_0^5}}\mathrm{e}^{-Zr/2a_0} \nonumber \]
in which \(Z\) is the atomic number (or probably better nuclear charge) and \(a_0\) is the Bohr radius.
In switching from spherical to Cartesian coordinates, we make the substitution \(z=r \cosθ\), so:
\[|\psi_{2,1,0} \rangle =z R(r) \nonumber \]
This is \(\psi_{2p_z}\) since the value of \(\psi \) is dependent on \(z\): when \(z=0\); \(\psi =0\), which is expected since \(z=0\) describes the \(xy\)-plane.
The other two wavefunctions are degenerate in the \(xy\)-plane. An equivalent statement is that these two orbitals do not lie on the x- and y-axes, but rather bisect them. Thus it is typical to take linear combinations of them to make the equation look prettier. If any set of wavefunctions is a solution to the Schrödinger equation, then any set of linear combinations of these wavefunctions must also be a solution (Section 2.4). We can do this because of the linearity of the Schrödinger equation.
In the equations below, we're going to make use of some trigonometry, notably Euler's formula:
\[ \begin{align} \mathrm{e}^{\mathrm{i}\phi} &=\cos{\phi}+\mathrm{i}\sin{\phi}\\[4pt] \sin{\phi} &= \dfrac{\mathrm{e}^{\mathrm{i}\phi}-\mathrm{e}^{-\mathrm{i}\phi}}{2\mathrm{i}}\\[4pt] \cos{\phi} &= \dfrac{\mathrm{e}^{\mathrm{i}\phi}+\mathrm{e}^{-\mathrm{i}\phi}}{2} \end{align} \nonumber \]
We're also going to use \(x=\sin θ\cos ϕ\) and \(y=\sin θ \sinϕ \).
\begin{align*} \psi_{2p_x} &=\dfrac{1}{\sqrt{2}}\left(\psi_{2,1,+1}-\psi_{2,1,-1}\right) \\[4pt] &=\dfrac{1}{2}\left(\mathrm{e}^{\mathrm{i}\phi}+\mathrm{e}^{-\mathrm{i}\phi} \right)r\sin{\theta} f(r) \\[4pt] &=r\sin{\theta}\cos{\phi}f(r)=xf(r) \\[4pt] \psi_{2p_y} &=\dfrac{\mathrm{i}}{\sqrt{2}}\left(\psi_{2,1,+1}+\psi_{2,1,-1}\right)\\[4pt] &=\dfrac{1}{2\mathrm{i}}\left(\mathrm{e}^{\mathrm{i}\phi}-\mathrm{e}^{-\mathrm{i}\phi} \right)r\sin{\theta}f(r)\\[4pt] &=r\sin{\theta}\sin{\phi}f(r)=yf(r)\\ \end{align*}
So, while \(m_l=0\) corresponds to \(|p_z \rangle\), \(m_l=+1\) and \(m_l=−1\) cannot be directly assigned to either \(|p_x \rangle\) or \(|p_y \rangle\), but rather a combination of \(|p_x \rangle\) and \(|p_y \rangle\). An alternative description is that \(m_l=+1\) might correspond to \((|p_x \rangle\ + |p_y \rangle )\) and \(m_l=−1\) might correspond to \((|p_x \rangle\ - |p_y \rangle)\).