Me encontré con esto en el libro Introductory Quantum Optics de Gerry & Knight como un ejercicio destinado al lector.
De todos modos, comencé con las soluciones proporcionadas a las ecuaciones de Maxwell sin una fuente/corriente $$ \begin{cases}E_x(t,z) = \left(\dfrac{2\omega^2}{V\varepsilon_0} \right)^{1/2}q(t)\sin(kz) \\ B_y(t,z) = \left(\dfrac{2\omega^2}{V\varepsilon_0} \right)^{1/2}\dfrac{\mu_0\varepsilon_0}{k}\dot{q}(t)\cos(kz) \end{cases} $$
Luego simplemente verifiqué que realmente satisfacían las ecuaciones de movimiento. Al hacerlo, noté que esto requiere que la función dinámica (también conocida como estado del campo) $q(t)$ satisfaga $\partial_t \dot{q}(t) \propto q(t)$, lo cual tiene sentido aquí ya que el campo $B$ debe satisfacer la ecuación de onda $(\partial_t^2-\partial_z^2)B=0 $ en unidades donde la velocidad $c$ de aquello que está ondeando es $1$. En estas unidades, su solución general tiene la forma $f(x+t)+g(x-t)$, lo cual es equivalente a un producto de un coseno con un argumento temporal y otro con un argumento espacial. Si esto no está claro, puedes consultar la Serie de Conferencias de Feynman Volumen II $\S$ 20.1 (enlace) para una revisión.
Lo siguiente es el Hamiltoniano para un campo de un solo modo que consiste en los campos anteriores $E$ y $B$, $$H = \dfrac{1}{2}\int dV \left[\varepsilon_0 E^2+\dfrac{1}{\mu_0}B^2 \right], $$ donde la integración es sobre una cavidad unidimensional a lo largo del eje $z$ con $z\in[0,L]$ y con cada campo desapareciendo en los límites. Escogiendo $\omega = kc$ (es decir, un solo modo), recordando $c^2 = (\varepsilon_0\mu_0)^{-1}$, y considerando $\dot{q}(t)= p$ como momento canónico, tenemos $$H = \dfrac{1}{L}\left[\omega^2 q^2\int_0^L\sin^2(kz)\hspace{0.1cm}dz+p^2\int_0^L\cos^2(kz)\hspace{0.1cm} dz \right]. $$ Eso reduce el problema a dos integrales. Siguiendo los datos de los límites, uno puede escoger $k = m\pi/L$ con $m\in \mathbb{N}$ de manera que $$\int_0^L \sin^2\dfrac{m\pi z}{L} \hspace{0.1cm} dz = \frac{1}{2}\int_0^L\left[ 1-\cos\frac{2m\pi z}{L}\right]dz = \dfrac{L}{2}$$ ya que solo sobrevive la integración del primer término. De manera similar, $\dfrac{L}{2}$ es todo lo que queda de la otra integral, resultando en \begin{align} H &= \frac{1}{L}\left(\omega^2 q^2 \dfrac{L}{2} + p^2 \dfrac{L}{2}\right)\\ \Rightarrow H &= \dfrac{1}{2}(p^2 + \omega^2 q^2), \end{align} mostrando que el Hamiltoniano para un campo de un solo modo es equivalente al de un oscilador armónico simple con masa unitaria. Aunque este cálculo fue completamente clásico, claramente sitúa a uno en un lugar familiar para empezar a cuantizar el campo electromagnético.