Entrar Página Discusión Historial Go to the site toolbox

Modelo de molécula de hidrógeno

De Laplace

(Diferencias entre revisiones)
(Trabajo total)
(Trabajo total)
Línea 151: Línea 151:
<center><math>W=W_1+W_2=\frac{e^2}{32\pi\varepsilon_0a}</math></center>
<center><math>W=W_1+W_2=\frac{e^2}{32\pi\varepsilon_0a}</math></center>
-
donde las contribuciones debidas a la otra carga puntual se cancelan mutuamente. La razón es sencilla: en este subsistema de dos cargas puntaulaes, primero las separamos y luego las volvemos a aproximar a la misma distancia que estaban, con lo que la distancia relativa no cambia y no se realiza trabajo alguno.
+
donde las contribuciones debidas a la otra carga puntual se cancelan mutuamente. La razón es sencilla: en este subsistema de dos cargas puntuales, primero las separamos y luego las volvemos a aproximar a la misma distancia que estaban, con lo que la distancia relativa no cambia y no se realiza trabajo alguno.
 +
 
 +
Esto era fácil de prever por lo que podíamos haber ido directamente al trabajo para mover dos cargas dentro de la nube esférica.
[[Categoría:Problemas de electrostática en el vacío]]
[[Categoría:Problemas de electrostática en el vacío]]

Revisión de 22:06 28 jun 2010

Contenido

1 Enunciado

Un modelo simple de la molécula de hidrógeno es el siguiente: tenemos dos cargas puntuales (los núcleos) de valor + e inmersas en una nube esférica de radio a, con carga − 2e distribuida uniformemente.

  1. Determine la posición de equilibrio entre las dos cargas puntuales, suponiendo que se encuentran situadas simétricamente respecto al centro de la nube.
  2. Calcule el potencial y el campo eléctrico en todos los puntos del espacio en la situación anterior.
  3. Suponga que las dos cargas positivas se desplazan una cantidad c = a / 4 a lo largo de la recta que las une, manteniendo la distancia entre ellas. En este caso, ¿qué campo se ve en el exterior de la molécula? ¿cuáles son los dos primeros momentos del desarrollo multipolar del potencial eléctrico?
  4. Calcule el trabajo necesario para realizar el desplazamiento del apartado anterior.

2 Posición de equilibrio

Las cargas positivas se repelen mutuamente. Sin embargo no se alejan indefinidamente debido a la fuerza debida a la nube negativa, que las empuja hacia el centro y que es más intensa cuanto más se alejan del centro.

La condición de equilibrio para una de las cargas positivas es que la fuerza total sobre ella se anule

\mathbf{F}_1=e\left(\mathbf{E}_q(\mathbf{r}_1)+\mathbf{E}_\rho(\mathbf{r}_1)\right)=\mathbf{0}

Si suponemos uno de los núcleos, de carga + e situado en \mathbf{r}_1=z_0\mathbf{u}_z y el otro en \mathbf{r}_2=-z_0\mathbf{u}_z, el campo que el segundo produce en la posición del primero es

\mathbf{E}_q(\mathbf{r}_1)=\frac{1}{4\pi\varepsilon_0}\,\frac{e(\mathbf{r}_2-\mathbf{r}_1)}{|\mathbf{r}_2-\mathbf{r}_1|^3}=\frac{1}{4\pi\varepsilon_0}\,\frac{e\mathbf{u}_z}{(2z_0)^2}

El campo debido a la nube de carga puede calcularse aplicando la ley de Gauss y el resultado es

\mathbf{E}=\begin{cases}\displaystyle\frac{(-2e)\mathbf{r}}{4\pi\varepsilon_0 a^3} & r < a \\ & \\ \displaystyle\frac{(-2e)\mathbf{r}}{4\pi\varepsilon_0 r^3} & r > a\end{cases}

En la posición ocupada por el primer núcleo, este campo vale

\mathbf{E}_\rho(\mathbf{r}_1) = -\frac{2ez_0\mathbf{u}_z}{4\pi\varepsilon_0a^3}

La condición de equilibrio es entonces

\mathbf{0}=e\left(\mathbf{E}_q(\mathbf{r}_1)+\mathbf{E}_\rho(\mathbf{r}_1)\right)=e\left(\frac{1}{4\pi\varepsilon_0}\,\frac{e\mathbf{u}_z}{(2z_0)^2}-\frac{2ez_0\mathbf{u}_z}{4\pi\varepsilon_0a^3}\right) = \frac{e^2\mathbf{u}_z}{4\pi\varepsilon_0}\left(\frac{1}{4z_0^2}-\frac{2z_0}{a^3}\right)

que tiene por solución

\frac{1}{4z_0^2}-\frac{2z_0}{a^3} = 0   \Rightarrow    z_0 = \frac{a}{2}

Por tanto, la molécula está formada por dos núcleos situados simétricamente a ambos lados del centro de la nube, a una distancia del centro igual a la mitad del radio de la nube.

3 Potencial y campo eléctrico

3.1 Campo eléctrico

El campo eléctrico en cualquier punto del espacio será la suma del de las dos cargas puntuales, más el de la nube de carga, que hemos enunciado más arriba. Si distinguimos entre el interior y el exterior de ésta tenemos:

Exterior de la molécula
Desde fuerza de la molécula, las dos cargas puntuales se ven como tales, mientras que la nube de carga se ve como una tercera carga puntual situada en el centro de la nube, y de magnitud − 2e. El campo total es entonces
\mathbf{E}(\mathbf{r})=\mathbf{E}_1(\mathbf{r})+\mathbf{E}_2(\mathbf{r})+\mathbf{E}_\rho(\mathbf{r})=\frac{e}{4\pi\varepsilon_0}\left(\frac{\mathbf{r}-(a/2)\mathbf{u}_z}{|\mathbf{r}-(a/2)\mathbf{u}_z|^3}+\frac{\mathbf{r}+(a/2)\mathbf{u}_z}{|\mathbf{r}+(a/2)\mathbf{u}_z|^3}-\frac{2\mathbf{r}}{r^3}\right)
Nótese que aunque la molécula es neutra (la nube tiene la misma carga que los dos núcleos juntos y con signo opuesto) el campo exterior no es nulo.
Interior de la molécula
Dentro de la molécula, las dos cargas puntuales se siguen viendo como puntuales, mientras que la nube de carga produce un campo que varía linealmente con la distancia. El campo total es entonces
\mathbf{E}(\mathbf{r})=\mathbf{E}_1(\mathbf{r})+\mathbf{E}_2(\mathbf{r})+\mathbf{E}_\rho(\mathbf{r})=\frac{e}{4\pi\varepsilon_0}
\left(\frac{\mathbf{r}-(a/2)\mathbf{u}_z}{|\mathbf{r}-(a/2)\mathbf{u}_z|^3}+\frac{\mathbf{r}+(a/2)\mathbf{u}_z}{|\mathbf{r}+(a/2)\mathbf{u}_z|^3}-\frac{2\mathbf{r}}{a^3}\right)

3.2 Potencial eléctrico

Para el potencial eléctrico aplicamos de nuevo el principio de superposición. El potencial total será el de las dos cargas más el de la nube. El potencial de cada carga puntual es

\phi_1(\mathbf{r})=\frac{e}{4\pi\varepsilon_0|\mathbf{r}-(a/2)\mathbf{u}_z|}        \phi_2(\mathbf{r})=\frac{e}{4\pi\varepsilon_0|\mathbf{r}+(a/2)\mathbf{u}_z|}

Para el potencial debido a la nube integramos el campo desde el infinito hasta una cierta distancia del centro.

Exterior de la nube
El campo es en todos los puntos el de una carga puntual, por lo que el potencial también lo es
\phi_\rho(r>a) = -\int_{\infty}^\mathbf{\mathbf{u}_r}\mathbf{E}\cdot\mathrm{d}\mathbf{r}= \int_\infty^r \frac{(-2e)\mathrm{d}r}{4\pi\varepsilon_0 r^2} = \frac{(-2e)}{4\pi \varepsilon_0 r}
Interior de la nube
Para llegar al interior, la integral se compone de dos tramos, uno desde el infinito hasta el borde de la nube, y otro desde ahí hasta el punto donde queremos hallar el potencial
\phi_\rho(r<a) = -\int_{\infty}^\mathbf{\mathbf{r}}\mathbf{E}\cdot\mathrm{d}\mathbf{r}= \int_\infty^a \frac{(-2e)\mathrm{d}r}{4\pi\varepsilon_0 r^2}+\int_a^r \frac{(-2e)r}{4\pi\varepsilon_0 a^3}\,\mathrm{d}r = \frac{(-2e)}{4\pi \varepsilon_0 a}+\frac{(-2e)}{4\pi \varepsilon_0 a^3}\left(\frac{r^2}{2}-\frac{a^2}{2}\right)=\frac{(-2e)(3a^2-r^2)}{8\pi\varepsilon_0 a^3}

Combinando ambos resultados:

\phi_\rho(r)=\begin{cases}\displaystyle -\frac{2e(3a^2-r^2)}{8\pi\varepsilon_0 a^3} & r < a \\ & \\ \displaystyle -\frac{2e}{4\pi\varepsilon_0 r} & r > a\end{cases}

Sumando ahora este potencial al de los dos núcleos obtenemos el potencial total:

Exterior de la molécula
Fuera de la molécula, el potencial es el de tres cargas puntuales
\phi(\mathbf{r})=\frac{e}{4\pi\varepsilon_0}\left(\frac{1}{|\mathbf{r}-(a/2)\mathbf{u}_z|}+\frac{1}{|\mathbf{r}+(a/2)\mathbf{u}_z|}-\frac{2}{r}\right)
Interior de la molécula
En el interior de la molécula, el potencial será el de las dos cargas puntuales más el debido a la nube, que varía cuadráticamente con r:
\phi(\mathbf{r})=\frac{e}{4\pi\varepsilon_0}\left(\frac{1}{|\mathbf{r}-(a/2)\mathbf{u}_z|}+\frac{1}{|\mathbf{r}+(a/2)\mathbf{u}_z|}-\frac{3a^2-r^2}{a^3}\right)

4 Desplazamiento de las cargas

4.1 Campo exterior

Tras el desplazamiento de las cargas, en el exterior de la molécula seguimos viendo tres cargas puntuales, sólo que ahora una de ellas está a una distancia a / 2 − a / 4 = a / 4 del centro y la otra a a / 2 + a / 4 = 3a / 4, por lo que el campo exterior es

\mathbf{E}(\mathbf{r})=\frac{e}{4\pi\varepsilon_0}\left(\frac{\mathbf{r}-(3a/4)\mathbf{u}_z}{|\mathbf{r}-(3a/4)\mathbf{u}_z|^3}+\frac{\mathbf{r}+(a/4)\mathbf{u}_z}{|\mathbf{r}+(a/4)\mathbf{u}_z|^3}-\frac{2\mathbf{r}}{r^3}\right)

4.2 Primeros momentos

Debemos hallar los momentos monopolar (carga) y dipolar de la molécula, formada por dos cargas puntuales, y la nube esférica. Sin embargo, el que la tercera carga sea en realidad una esfera es irrelevante desde el punto de vista del cálculo de los momentos, ya que el desarrollo multipolar sólo se aplica en puntos alejados de la distribución y para cualquier punto exterior a una esfera de carga uniforme, ésta se ve como una carga puntual situada en su centro. Por ello, para el desarrollo multipolar, el sistema equivale a tres cargas puntuales:

  • Una carga puntual q1 = + e situada en \mathbf{r}_1=3a\mathbf{u}_z/4.
  • Una carga puntual q2 = + e situada en \mathbf{r}_2=-a\mathbf{u}_z/4.
  • Una carga puntual q3 = − 2e situada en \mathbf{r}_3=\mathbf{0}.

Así tenemos

Momento monopolar
Es la carga neta de la distribución
Q=e+e-2e = 0\,
Es nula, como era de esperar para una molécula formada por dos átomos con sus respectivos electrones.
Momento dipolar
\mathbf{p}=\sum_i q_i\mathbf{r}_i = (+e)\left(\frac{3a}{4}\mathbf{u}_z\right)+(+e)\left(-\frac{a}{4}\mathbf{u}_z\right)+(-2e)\left(\mathbf{0}\right)=\frac{ea}{2}\mathbf{u}_z

5 Trabajo eléctrico

La forma más sencilla de calcular el trabajo para realizar el desplazamiento es imaginar un proceso que lleve desde el estado inicial al final. Este proceso se compone de dos pasos:

  1. La carga q1 se lleva de \mathbf{r}_1=(a/2)\mathbf{u}_z a \mathbf{r}'_1=(3a/4)\mathbf{u}_z.
  2. La carga q2 se lleva de \mathbf{r}_2=-(a/2)\mathbf{u}_z a \mathbf{r}'_2=-(a/4)\mathbf{u}_z.

Hay que hacerlo en dos pasos porque si movemos dos cargas a la vez podemos tener problemas al suponer un potencial dependiente del tiempo (el de la otra carga que se está moviendo a su vez).

5.1 Desplazamiento de la primera carga

El trabajo para mover la primera carga es

W_1 = q_1\left(\phi'_1(\mathbf{r}'_1)-\phi'_1(\mathbf{r}_1)\right)

donde φ' representa el potencial creado por el resto de cargas del sistema (en este caso, el otro núcleo y la nube)

\phi'_1(\mathbf{r})=\frac{e}{4\pi\varepsilon_0}\left(\frac{1}{|\mathbf{r}+(a/2)\mathbf{u}_z|}-\frac{3a^2-r^2}{a^3}\right)

Calculamos la diferencia de potencial entre la posición inicial y la final

\phi'_1\left(\frac{3a}{4}\mathbf{u}_z\right)-\phi'_1\left(\frac{a}{2}\mathbf{u}_z\right)=\frac{e}{4\pi\varepsilon_0}\left(\left(\frac{4}{5a}-\frac{39}{16a}\right)-\left(\frac{1}{a}-\frac{11}{4a}\right)\right)=\left(-\frac{1}{5}+\frac{5}{16}\right)\frac{e}{4\pi\varepsilon_0a}

Mantenemos separadas las dos fracciones para poner de manifiesto cuánto procede de la interacción entre las dos cargas puntuales y cuánto de la nube.

El primer trabajo es entonces

W_1= \left(-\frac{1}{5}+\frac{5}{16}\right)\frac{e^2}{4\pi\varepsilon_0a}

5.2 Desplazamiento de la segunda carga

Del mismo modo, el trabajo para mover la segunda carga es

W_2 = q_2\left(\phi'_2(\mathbf{r}'_2)-\phi'_2(\mathbf{r}_2)\right)

donde φ'2 representa el potencial creado la nube y el otro núcleo que ahora se encuentra en z = 3a / 4

\phi'_2(\mathbf{r})=\frac{e}{4\pi\varepsilon_0}\left(\frac{1}{|\mathbf{r}-(3a/4)\mathbf{u}_z|}-\frac{3a^2-r^2}{a^3}\right)

Calculamos la diferencia de potencial entre la posición inicial y la final

\phi'_2\left(-\frac{a}{4}\mathbf{u}_z\right)-\phi'_2\left(-\frac{a}{2}\mathbf{u}_z\right)=\frac{e}{4\pi\varepsilon_0}\left(\left(\frac{1}{a}-\frac{47}{16a}\right)-\left(\frac{4}{5a}-\frac{11}{4a}\right)\right)=\left(\frac{1}{5}-\frac{3}{16}\right)\frac{e}{4\pi\varepsilon_0a}

con lo que el segundo trabajo es

W_2= \left(\frac{1}{5}-\frac{3}{16}\right)\frac{e^2}{4\pi\varepsilon_0a}

5.3 Trabajo total

Sumando las dos contribuciones resulta un trabajo total

W=W_1+W_2=\frac{e^2}{32\pi\varepsilon_0a}

donde las contribuciones debidas a la otra carga puntual se cancelan mutuamente. La razón es sencilla: en este subsistema de dos cargas puntuales, primero las separamos y luego las volvemos a aproximar a la misma distancia que estaban, con lo que la distancia relativa no cambia y no se realiza trabajo alguno.

Esto era fácil de prever por lo que podíamos haber ido directamente al trabajo para mover dos cargas dentro de la nube esférica.

Herramientas:

Herramientas personales
TOOLBOX
LANGUAGES
licencia de Creative Commons
Aviso legal - Acerca de Laplace