The motion of a single charged particle in an ideal Penning trap (with perfect alignment of its magnetic field to the quadrupole potential) is an exactly solvable system in both classical and quantum mechanics. The particle's motion along the trap's axis is simple harmonic motion, and the motion in the trap's xy-plane is a perturbation of cyclotron motion that reduces to cyclotron motion exactly in the zero-electric-field limit.
Classical Treatment We assume for simplicity a positive charge q on the trapped particle of mass m. The electric potential has the form : \phi(x,y,z) = \frac{V_0}{2d^2} \left( z^2 - \frac{x^2+y^2}{2} \right) where V_0 is the potential difference between the "endcap" and "ring" electrodes, and d is a length scale related to the spacing between these electrodes. From the gradient of \phi we obtain the electric field as :\vec E = -\nabla \phi = -\frac{V_0}{d^2}\left(z\hat z -\frac{x\hat x + y\hat y}{2}\right) = -\frac{V_0}{d^2}\left(z\hat z -\frac{1}{2}\vec \rho\right) where \vec \rho is the projection of \vec r into the xy-plane. The magnetic field \vec B = -B\hat z is completely uniform, and by convention we direct it "into the page" along the negative z-axis so that the resulting motion for our positively-charged particle is counter-clockwise. The axial and planar equations of motion can then be obtained directly from the Lorentz Force law \vec F = q(\vec E + \vec v \times \vec B) as :\ddot{z} = -\omega_z^2 z :\ddot\vec\rho = \vec\omega_c \times \dot\vec\rho + \frac{1}{2}\omega_z\vec\rho Where \omega_z and \omega_c are the axial and cyclotron frequencies of the system respectively, given by :\omega_z = \sqrt{\frac{qV_0}{md^2}}, \omega_c = \frac{qB}{m}. These two equations of motion are completely decoupled. Therefore, axially the system is a
simple harmonic oscillator with frequency \omega_z. The planar equation of motion is also exactly solvable. One may use the complex variable u=x+iy to represent planar coordinates, so that the equation of motion takes the form :\ddot u(t) = i \omega_c \dot u(t) +\frac{1}{2} \omega_z^2 u(t) which, upon
Fourier transform to the frequency domain, yields : -\omega^2 \tilde u(\omega) = - \omega_c \omega \tilde u(\omega) + \frac{1}{2} \omega_z^2 \tilde u(\omega) :\rightarrow \omega^2 - \omega_c \omega +\frac{1}{2}\omega_z^2 = 0 The application of the
quadratic formula to this expression shows that the planar motion decomposes into exactly two possible frequency modes: :\omega_\pm = \frac{1}{2}\omega_c \pm \frac{1}{2}\sqrt{\omega_c^2-2\omega_z^2} which are the aforementioned
modified cyclotron mode (\omega_+) and
magnetron mode (\omega_-). The general solution for the planar motion is therefore :u(t)=c_1 e^{i\omega_+ t} + c_2 e^{i\omega_- t} . One can see from the above analysis that \omega_+ and \omega_- are only strictly real numbers in the case that \omega_c \geq \sqrt 2 \omega_z, which may also be taken as the condition for the particle to remain trapped within a finite radius. The substitution of non-real \omega_+ and \omega_- into the general solution produces orbits that escape to infinity (while still satisfying the equation of motion). Substituting the definitions for the axial and cyclotron frequencies, one might also write this condition as :B \geq \sqrt\frac{2m V_0}{qd^2} which shows that a condition for any charged particle to remain confined in a Penning trap is for the magnetic field to be sufficiently strong, while also having a nonzero voltage differential so that axial confinement is possible.
Quantum-Mechanical Treatment The quantum-mechanical treatment of a single charged particle moving in a penning trap requires the construction of a
Hamiltonian for the system, which means the magnetic field \vec B = -B\hat z has to be represented through a choice of
magnetic vector potential. Since the magnetic field is totally uniform, and the system is axially symmetric, the "Landau Gauge" is particularly convenient: :\vec A = \frac{1}{2} \vec B \times \vec r Component-wise, this is A_x = -\frac{1}{2}By, A_y = \frac{1}{2}Bx, A_z = 0. We plug this and the electric potential \phi from the classical treatment into the general Hamiltonian for a spinless charged particle in an electromagnetic field : \hat H = \frac{({\bf \hat p} - q{\bf A}(\hat x,\hat y,\hat z))^2}{2m} + q\phi(\hat x,\hat y,\hat z) to obtain :\hat H = \frac{\hat p_x^2+\hat p_y^2+\hat p_z^2}{2m}+\frac{\omega_c}{2}(\hat x \hat p_y - \hat y \hat p_x)+\frac{1}{2}m\left(\frac{\omega_1}{2}\right)^2(\hat x^2 + \hat y^2) + \frac{1}{2}m\omega_z^2 \hat z^2 where \omega_z and \omega_c are defined the same way as in the classical treatment, and \omega_1 = \sqrt{\omega_c^2 - 2 \omega_z^2} has been introduced as an additional shorthand for the square root of the discriminant that separates \omega_c into \omega_+ and \omega_-. From here, one might attempt to solve the
Schrödinger equation for this Hamiltonian directly in its current form. However, it proves far more analytically bountiful to take advantage of the fact that this system is actually isomorphic to three uncoupled harmonic oscillators. One need only define the ladder operators :\hat a_+ = \frac{1}{2\sqrt{\hbar}}\left( \sqrt\frac{m\omega_1}{2}(\hat x - i \hat y) + \sqrt\frac{2}{m\omega_1}(\hat p_y + i \hat p_x) \right) :\hat a_- = \frac{1}{2\sqrt{\hbar}}\left( \sqrt\frac{m\omega_1}{2}(\hat y - i \hat x) + \sqrt\frac{2}{m\omega_1}(\hat p_x + i \hat p_y) \right) :\hat a_z = \frac{1}{\sqrt{2\hbar}}\left(\sqrt{m\omega_z}\hat z + i\frac{1}{\sqrt{m\omega_z}}\hat p_z\right) Which, along with their
Hermitian conjugate operators \hat a_+^\dagger, \hat a_-^\dagger, and \hat a_z^\dagger, raise and lower the energy level of each of the system's three modes (modified cyclotron, magnetron, and axial respectively), just like the raising and lowering operators from the
quantum simple harmonic oscillator. These obey the commutation relations :[\hat a_+,\hat a_+^\dagger] = [\hat a_-,\hat a_-^\dagger] = [\hat a_z,\hat a_z^\dagger] = 1, with all commutators between operators for different modes vanishing. In terms of these ladder operators, the Hamiltonian can be rewritten in the form :\hat H = \hbar \omega_+ \left(\hat a_+^\dagger \hat a_+ + \frac{1}{2}\right) - \hbar \omega_- \left(\hat a_-^\dagger \hat a_- + \frac{1}{2}\right) + \hbar \omega_z \left(\hat a_z^\dagger \hat a_z + \frac{1}{2}\right) and by isomorphism with the simple harmonic oscillator, we see that all states of this system are bound states, with discrete energy levels :E_{n_+ n_- n_z} = \hbar \omega_+ \left(n_+ + \frac{1}{2}\right) - \hbar \omega_- \left(n_- + \frac{1}{2}\right) + \hbar \omega_z \left(n_z + \frac{1}{2}\right) where n_+, n_- and n_z are quantum numbers describing the excitation of the system's three modes. ==Fourier-transform mass spectrometry==