Modelo de campo atómico
De Laplace
1 Enunciado
El potencial medio temporal de un átomo de hidrógeno neutro viene dado por
![\phi = \frac{q \mathrm{e}^{-\alpha r}}{4\pi\varepsilon_0 r} \left( 1 + \frac{\alpha
r}{2} \right)](/wiki/images/math/3/4/b/34bb91a2b22fcc8ca0cf1c8796290de5.png)
en donde q es la carga electrónica, y α − 1 = a0 / 2. Halle la distribución de carga (continua y discreta) que dará lugar a este potencial e interprete este resultado físicamente.
2 Solución
Tenemos que el potencial posee simetría de rotación, por lo que todas las derivadas e integrales van a ser sobre la coordenada esférica radial r
![\phi = \frac{q}{4\pi\varepsilon_0 r}\mathrm{e}^{-\alpha r}\left(1+\frac{\alpha r}{2}\right)](/wiki/images/math/4/d/7/4d764fa0ed2146f58ad916ae45c1f6ef.png)
Para la densidad de carga de volumen
![\rho = -\varepsilon_0\nabla^2\phi =
-\frac{\varepsilon_0}{r}\frac{\mathrm{d}^2\ }{\mathrm{d}
r^2}\left(r\phi\right)](/wiki/images/math/a/5/8/a584c639f631bc93794539341f659ccd.png)
Calculamos esta derivada. Multiplicando por r
![r\phi = q\mathrm{e}^{-\alpha
r}\left(1+\frac{\alpha}{2}r\right) = \mathrm{e}^{-\alpha
r}f(r)](/wiki/images/math/3/1/6/316bec9cba279507c2fe90b407d675ef.png)
la derivada de este producto es igual a
![\frac{\mathrm{d}^2\ }{\mathrm{d} r^2}\left(r\phi\right) =
\mathrm{e}^{-\alpha r}\left(f''(r)-2\alpha f'(r)+\alpha^2
f(r)\right)](/wiki/images/math/0/0/2/0022e828ba8d2d7d49574ccd3396f728.png)
Sustituyendo f(r)
![\frac{\mathrm{d}^2\ }{\mathrm{d} r^2}\left(r\phi\right) =
q\mathrm{e}^{-\alpha
r}\left(0-2\alpha\left(\frac{\alpha}{2}\right)+\alpha^2\left(1+\frac{\alpha}{2}r\right)\right)
= \frac{q \alpha^3r \mathrm{e}^{-\alpha r}}{2}](/wiki/images/math/9/2/5/925cdab9cdc30f75cd9ad56992183fea.png)
y la densidad de carga de volumen es
![\rho = -\frac{q \alpha^3 \mathrm{e}^{-\alpha
r}}{8\pi}](/wiki/images/math/e/2/9/e296687d6cf741beaedeb20fd68258fe.png)
Queda la carga del núcleo, que es una carga puntual. Esta la sacamos de que el sistema es neutro, (ya que el flujo del campo eléctrico tiende a 0 para ). La carga de volumen podemos escribirla también como
![\rho = -\frac{\varepsilon_0}{r^2}\frac{\mathrm{d}\ }{\mathrm{d} r}\left(r^2\frac{\mathrm{d}\phi}{\mathrm{d}
r}\right)](/wiki/images/math/0/e/c/0eccd073924d10612196e43f1e4994f7.png)
así que la carga total de volumen es
![Q_v = \int \rho\,\mathrm{d}\tau = 4\pi \int_0^\infty\rho\,r^2\,\mathrm{d}r = -4\pi\varepsilon_0 \int_0^\infty\frac{\mathrm{d}\ }{\mathrm{d}r}\left(r^2\frac{\mathrm{d}\phi}{\mathrm{d}r}\right)\,\mathrm{d}r = 4\pi\varepsilon_0 \lim_{r\to 0}r^2\frac{\mathrm{d}\phi}{\mathrm{d}r}](/wiki/images/math/5/1/0/510a1c1877790c337f28a1a0000e5b05.png)
![\phi = \frac{g(r)}{r}](/wiki/images/math/8/f/0/8f0d09f1475e3a024e80487cd88e3509.png)
![r^2 \frac{\mathrm{d}\phi}{\mathrm{d}r}= r^2 \frac{r
g'(r)-g(r)}{r^2} = r g'(r)-g(r)](/wiki/images/math/2/b/4/2b4c94e9ec7fa86754fb6f5570459738.png)
![Q_v = 4\pi\varepsilon_0 \lim_{r\to 0}
r^2\frac{\mathrm{d}\phi}{\mathrm{d}r} =
-4\pi\varepsilon_0\frac{q}{4\pi\varepsilon_0}=-q](/wiki/images/math/d/3/5/d358818bf01a16773a8e6fea817ee167.png)
Por tanto la carga del núcleo es simplemente + q. Así pues, tenemos
- Una carga puntual + q situada en el núcleo
- Una carga − q distribuida de forma exponencial en todo el espacio
![\rho = -\frac{q\alpha^3\mathrm{e}^{-\alpha
r}}{8\pi}](/wiki/images/math/e/2/9/e296687d6cf741beaedeb20fd68258fe.png)
o, juntándolo todo en una sola distribución
![\rho =
q\left(\delta(\mathbf{r})-\frac{\alpha^3\mathrm{e}^{-\alpha
r}}{8\pi}\right)](/wiki/images/math/6/0/d/60d2a99d37353169e11763ba8856a040.png)