La versión rápida y sucia es que se modelan todas las partículas de un tipo determinado como excitaciones de una serie de osciladores armónicos cuánticos:
$$ H = \int\frac{\mathrm{d}^{3}\vec{p}}{(2\pi)^3} E_{\vec{p}} \left(a_{\vec{p}}^{\dagger}a_{\vec{p}} + \frac{1}{2}\right) $$
por lo que una partícula de momento $\vec{p}$ sería el $|1\rangle$ estado del oscilador armónico de momento $\vec{p}$ . Nota $E_{\vec{p}}^2 - \vec{p}^2 = m^2$ en unidades naturales y $E_{\vec{p}}$ es una frecuencia angular por la relación de Broglie. Para simplificar esto se define una cosa llamada "operador de campo" que permite trabajar en el espacio de la posición en lugar del momento:
$$ \phi = \int\frac{\mathrm{d}^{3}\vec{p}}{(2\pi)^3}\frac{1}{2E_{\vec{p}}}\left(a_{\vec{p}} \mathrm{e}^{-ipx} + a^{\dagger}_{\vec{p}} \mathrm{e}^{ipx}\right) $$
donde $p$ y $x$ sin flechas indica cuatro vectores y cuatro posiciones. Si se introduce esto y se hace el álgebra, se obtiene el hamiltoniano teórico de campo estándar para un campo libre (escalar):
$$ H = \int \mathrm{d}^{3}\vec{x}\left(\left(\frac{\partial\phi}{\partial t}\right)^2 + \left(\nabla\phi\right)^2 - m^2\phi^2\right) $$
El espacio de Hilbert para este hamiltoniano es justo lo que se espera de un conjunto de osciladores armónicos:
$$ \mathcal{H} = \bigotimes_{\vec{p}}\mathcal{H}_{\vec{p}} $$
donde $\mathcal{H}_{\vec{p}}$ es el espacio de Hilbert para un solo oscilador armónico, y en la expresión para el Hamiltoniano hemos suprimido realmente una serie incontable de $\otimes \mathbb{I} \otimes$ antes y después de cada operador de escalera. A veces la gente llama a esto un espacio de Fock, pero no es realmente un espacio de Fock. Tiene propiedades similares, pero su construcción es muy diferente [1].
Para la dinámica, se utiliza la imagen de Heisenberg y, en particular, la ecuación de Heisenberg ( no la ecuación de Schrodinger):
$$ \frac{\mathrm{d}\phi}{\mathrm{d}t} = i\left[H,\phi\right] \\ \frac{\mathrm{d}\pi}{\mathrm{d}t} = i\left[H,\pi\right] $$
donde $\pi = \frac{\partial\phi}{\partial t}$ es el momento conjugado del campo definido de la forma habitual a partir de la Lagrangiana. Una vez más, al revisar el álgebra se encuentra que el campo obedece a la ecuación de Klein-Gordon:
$$ \left(\Box + m^2\right)\phi = 0 $$
Naturalmente, esta es una afirmación bastante extraña para hacer sobre el universo. ¿Por qué todas las partículas son excitaciones de un oscilador armónico? ¿Es sólo una aproximación, como tantas cosas en la física que se modelan mediante osciladores armónicos, o hay algo más fundamental?
Obviamente, la respuesta es que hay algo más fundamental. Para verlo, hay que fijarse en la estructura geométrica diferencial de la variedad del espaciotiempo, y en particular en las diferentes representaciones de su grupo de isotropía (el grupo de Lorentz). Al hacer esto, se ve que la imagen del espacio de posición es el punto de partida natural y, sorprendentemente, se convierte en osciladores armónicos cuando se hace una transformación de Fourier. En esencia, éste es el verdadero formalismo matemático de la cuantificación canónica.
Estoy encantado de entrar en los detalles técnicos de esta construcción si quieres (también explica los campos vectoriales y los campos espinores, cosa que el enfoque anterior no hace), pero es sobre todo de interés matemático y filosófico más que algo práctico con los cálculos. (También es útil si quieres ver la unificación y esas cosas, supongo).
[1]: En particular, $\mathcal{H}$ viene ya equipado con la idea de indistinguibilidad incorporada, porque si vas a llamar al estado $|2\rangle$ de un oscilador armónico un estado de 2 partículas (donde ambas tienen el mismo momento), ya no existe el concepto de "qué partícula es 1 y cuál es 2".