De hecho, los electrones (al menos los de la capa s) hacer pasan un tiempo no trivial dentro del núcleo.
La razón por la que pasan mucho tiempo fuera del núcleo es esencialmente mecánica cuántica. Para usar una explicación demasiado simple, su momento está restringido a un rango consistente con ser capturado (no libre para volar), y como tal hay una incertidumbre necesaria en su posición.
Un ejemplo de física que surge porque pasan algunos tiempo en el núcleo se llama así "captura beta" desintegración radiactiva en la que $$ e + p \to n + \nu $$ se produce dentro del núcleo. La razón por la que esto no ocurre en la mayoría de los núcleos es también mecánica cuántica y está relacionada con los niveles de energía y la exclusión de Fermi.
Para ampliar un poco esta imagen, recurramos a de Broglie y Bohr. La imagen de Bohr de las órbitas de los electrones restringidas a un conjunto de energías finitas $E_n \propto 1/n^2$ y las frecuencias pueden tener una explicación razonablemente natural en términos de la imagen de de Broglie de que toda la materia está compuesta por ondas de frecuencia $f = E/h$ exigiendo que un número entero de ondas encaje en la órbita circular.
Esto conduce a una imagen del átomo en la que todos los electrones ocupan órbitas circulares ordenadas y alejadas del núcleo, y proporciona una explicación de por qué los electrones no caen simplemente en el núcleo bajo la atracción electrostática.
Pero no es toda la historia por varias razones; para nuestros propósitos la más importante es que el modelo de Bohr predice un momento angular mínimo para los electrones de $\hbar$ cuando el valor experimental es 0.
Siguiendo adelante, podemos resolver la ecuación tridimensional de Schrödinger en tres dimensiones para los átomos de tipo hidrógeno:
$$ \left( i\hbar\frac{\partial}{\partial t} - \hat{H} \right) \Psi = 0 $$
para los electrones en un $1/r^2$ potencial electrostático para determinar la función de onda $\Psi$ . La función de onda está relacionada con la probabilidad $P(\vec{x})$ de encontrar un electrón en un punto $\vec{x}$ en el espacio por
$$ P(\vec{x}) = \left| \Psi(\vec{x}) \right|^2 = \Psi^{*}(\vec{x}) \Psi(\vec{x}) $$
donde $^{*}$ significa el complejo conjugado.
Las soluciones se suelen escribir de la forma
$$ \Psi(\vec{x}) = Y^m_l(\theta,\phi) L^{2l+1}_{n-l-1}(r) e^{-r/2} * \text{normalizing factors} $$
Aquí el $Y$ son los armónicos esféricos y los $L$ son los polinomios de Laguerre generalizados. Pero no nos interesan los detalles. Basta con decir que estas soluciones representan una densidad de probabilidad para los electrones que se extiende por una amplia zona cerca del núcleo. También hay que tener en cuenta que para $l=0$ estados (también conocidos como orbitales s) existe una probabilidad no nula en el centro, es decir en el núcleo (este hecho se debe a que estos orbitales tienen un momento angular nulo, lo cual, como se recordará, no era una característica del átomo de Bohr).