Movimiento de una molécula en un campo eléctrico

Vamos a ver como es posible acelerar una molécula empleando un campo eléctrico. En la figura, se muestra un par de esferas con cargas +Q y –Q respectivamente, separadas una distancia 2D, y una molécula de momento dipolar p=q·d, siendo 2d la separación entre los centros de distribución de las cargas positivas y negativas, respectivamente.

La molécula está sometida a cuatro fuerzas: dos de repulsión entre cargas del mismo signo y otras dos de atracción entre cargas de signo contrario. Como las cargas de signo contrario están más cerca, predomina la fuerza de atracción sobre la de repulsión.

Si en el instante en el que la molécula llega al centro de las esferas, se conectasen a potencial cero, la molécula continuaría con velocidad constante.

Un dispositivo de este tipo formado por varias etapas constituye un acelerador de moléculas que permite estudiar su estructura a partir de los choques inelásticos que experimentan.

Un acelerador de moléculas tiene unos electrodos de forma esférica de 0.25 mm de radio, separados 2D=1 mm, que se mantienen a un potencial de ±40 kV. Una molécula de momento dipolar p=2·10^-29 C·m adquiere una energía de aproximadamente 0.01 eV. Como esta energía es muy pequeña, el acelerador consta de 700 etapas separadas 1.4 cm, los electrodos se alimentan por un potencial alterno de 500 kHz. (véase Lorrain P. Corson D., Campos y Ondas Electromagnéticas, Editorial Selecciones Científicas (1972), págs. 144-145).

En la experiencia simulada mantendremos la carga de los dos electrodos esféricos constante e igual a ±Q. Observaremos como la molécula describe un movimiento oscilatorio, que no es Armónico Simple (MAS).

Ecuación del movimiento

Por simetría, las componentes a lo largo del eje Y de las fuerzas que actúan sobre la molécula se anulan de dos en dos. La resultante está dirigida a lo largo del eje X y vale

F=-2·F_a·cosθ_a+2·F_r·cosθ_r

Siendo F_a el módulo de la fuerza atractiva entre cargas de distinto signo, y F_r el módulo de la fuerza repulsiva entre cargas del mismo signo.

F a = 1 4π ε 0 qQ ( x 2 + (D−d) 2 ) F r = 1 4π ε 0 qQ ( x 2 + (D+d) 2 )

Los ángulos que forman los vectores fuerza con el eje X son, respectivamente

cos⁡ θ a = x x 2 + (D−d) 2 cos⁡ θ r = x x 2 + (D+d) 2

La expresión de la fuerza resultante es

F= qQx 2π ε 0 ( 1 ( x 2 + (D+d) 2 ) 3/2 − 1 ( x 2 + (D−d) 2 ) 3/2 )

Si x>0, la resultante de las fuerzas F<0
Si x<0, la resultante de las fuerzas F>0

La fuerza F es de signo contrario al desplazamiento x, pero no es proporcional al desplazamiento, que es la característica distintiva de un MAS.

La ecuación del movimiento es

m d 2 x d t 2 = qQx 2π ε 0 ( 1 ( x 2 + (D+d) 2 ) 3/2 − 1 ( x 2 + (D−d) 2 ) 3/2 )

Esta ecuación diferencial se resuelve por procedimientos numéricos con las condiciones iniciales t=0, x=x₀ dx/dt=0. La molécula parte de la posición x₀ con velocidad inicial nula.

Balance energético

En la posición x, la energía de la molécula es la suma de la energía cinética y potencial

E= 1 2 m v 2 + 1 4π ε 0 ( 2qQ x 2 + (D+d) 2 − 2qQ x 2 + (D−d) 2 )

Cuando pasa por el origen, la energía de la molécula es

E= 1 2 m v 2 + 1 4π ε 0 ( 2qQ D+d − 2qQ D−d )

Cuando parte de la posición inicial x₀ o llega al extremo opuesto de la trayectoria –x₀, la velocidad de la molécula v=0.

E= 1 4π ε 0 ( 2qQ x 0 2 + (D+d) 2 − 2qQ x 0 2 + (D−d) 2 )

Aplicando el principio de conservación de la energía, calculamos la velocidad v de la molécula en cualquier posición x.

Actividades

En la simulación se fijado los siguientes parámetros:

1/(4πє₀)=1,
La masa del la molécula es m=1, y sus cargas ±q=±1. La separación d es una fracción de la separación 2D entre las esferas.
La carga Q de las esferas es múltiplo de q, y su separación es 2D=1 unidad de longitud.
La posición de la molécula en el instante t=0, es x₀=3 unidades de longitud

Se introduce

La carga Q de las dos esferas, seleccionado un número en el control de selección titulado Carga
La separación entre las cargas ±q de la molécula, seleccionado una fracción en el control de selección titulado Separación

Se pulsa el botón titulado Empieza

Se observa el movimiento de la molécula entre las posiciones extremas ±x₀, y se puede medir el periodo P de la oscilación. La fuerza sobre la molécula es de corto alcance ya que solamente tiene un valor apreciable cuando pasa por las proximidades del origen.

En la parte superior del applet, se representa la energía potencial E_p(x), que como podemos observar presenta un mínimo en x=0.

E p (x)= 1 4π ε 0 ( 2qQ x 2 + (D+d) 2 − 2qQ x 2 + (D−d) 2 )

Se representa mediante una recta horizontal la energía total E. Un segmento vertical de color rojo, representa la energía cinética E_k, diferencia entre la energía total E y la energía potencial E_p(x). Una flecha representa la fuerza sobre la partícula, la pendiente cambiada de signo de la curva de la energía potencial en la posición x.

F=− d E p (x) dx

La fuerza F es nula en el origen x=0 (posición de equilibrio), donde la energía potencial E_p(x) presenta un mínimo

Ejemplo:

Datos

Separación entre los electrodos esféricos, 2D=1
Separación entre los centros de distribución de carga de la molécula polar 2d=1/2
Cargas de la molécula polar q=1
Carga de los electrodos Q=6
Valor de la constante 1/(4πє₀)=1

La fuerza sobre la partícula en la posición x=2 es

F= qQx 2π ε 0 ( 1 ( x 2 + (D+d) 2 ) 3/2 − 1 ( x 2 + (D−d) 2 ) 3/2 ) F=2 1·6·2 1 ( 1 ( 2 2 + (0.5+0.25) 2 ) 3/2 − 1 ( 2 2 + (0.5−0.25) 2 ) 3/2 )=−0.468

La energía de la partícula en la posición inicial x₀=3 es

E= 1 4π ε 0 ( 2qQ x 0 2 + (D+d) 2 − 2qQ x 0 2 + (D−d) 2 ) E=1·( 2·1·6 3 2 + (0.5+0.25) 2 − 2·1·6 3 2 + (0.5−0.25) 2 )=−0.1056

La energía cuando pasa por la posición x=2 es

E= 1 2 m v 2 + 1 4π ε 0 ( 2qQ x 2 + (D+d) 2 − 2qQ x 2 + (D−d) 2 ) -0.1056= 1 2 1· v 2 +1·( 2·1·6 2 2 + (0.5+0.125) 2 − 2·1·6 2 2 + (0.5−0.125) 2 )

Se despeja la velocidad de la molécula v=0.678

LineasApplet aparecerá en un explorador compatible JDK 1.1