Vibración de una molécula diatómica

prev.gif (1231 bytes)home.gif (1232 bytes)next.gif (1211 bytes)

Oscilaciones

Osciladores
Oscilaciones libres
y amortiguadas
Oscilaciones forzadas
Oscilador no lineal
Oscilaciones
anarmónicas
marca.gif (847 bytes)Molécula diatómica
El oscilador caótico

Osciladores acoplados
Modos normales
de vibración
De las oscilaciones
a las ondas
java.gif (886 bytes)Molécula diatómica

java.gif (886 bytes)El potencial de Lennard-Jones

 

Molécula diatómica

Dos átomos ligados formando una molécula estable tendrán una energía potencial Ep(x) que es una función cuadrática de la diferencia entre su separación x y la separación de equilibrio x0.

molecula2.gif (1216 bytes) La fuerza F que describe la interacción entre los dos átomos está dirigida a lo largo de la línea que conecta los dos átomos y viene dada por

El movimiento de un sistema aislado de dos partículas lo podemos reducir al movimiento de una partícula de masa m que se mueve alrededor de un centro fijo de fuerzas bajo la acción de la fuerza de interacción mutua.

m ·a12=F12

donde m es la masa reducida, a12 es la aceleración relativa de la partícula de masa m1 respecto de la partícula de masa m2, y F12 es la fuerza F(x) que describe la interacción mutua. En forma de ecuación diferencial

Esta es la ecuación diferencial de un Movimiento Armónico Simple de frecuencia angular

Ejemplo:

Determinar la constante k de fuerza de interacción entre los átomos de la molécula HCl sabiendo que su frecuencia de vibración f0=9·1013 Hz.

Datos: la masa del hidrógeno m1=1 y la masa del cloro m2=35 u.m.a. Una unidad de masa atómica vale 1.67·10-27 kg.

Teniendo en cuenta que w0=2p f0, obtenemos k=519.2 N/m

Movimiento de los átomos

La solución de la ecuación diferencial es

x=x0+A·sen(w0t+j )

Donde x0 es la separación de equilibrio, A es la amplitud y j la fase inicial que se determinan a partir de las condiciones iniciales, es decir, la posición x y la velocidad v en el instante t=0.

molecula1.gif (1728 bytes)

Conocido el movimiento relativo de la partícula de masa m1 respecto de la partícula de masa m2, y suponiendo que el centro de masas de este sistema aislado está en reposo, podemos determinar fácilmente el movimiento de cada una de las dos partículas.

 

Actividades

En el programa interactivo se ha fijado el valor de k=2p2 y se permite variar las masas de los átomos m1 y m2 de la molécula diatómica.

Por ejemplo, cuando introducimos las masas de las moléculas m1=m2=1, la frecuencia angular de vibración vale w 0=2p , y el periodo P=1 unidad de tiempo.

ForzadasApplet aparecerá en un explorador compatible con JDK 1.1.
                 
 

El potencial de Lennard-Jones

La energía potencial Ep(r) de una molécula debida a las fuerzas intermoleculares es una función de r, la distancia entre los centros de las dos moléculas interactuantes. Esta función se puede expresar mediante una función conocida como potencial de Lennard-Jones.

El potencial de Lennard- Jones tiene una forma similar al potencial de Morse, una fórmula empírica que describe bastante bien la energía potencial del estado ligado de una molécula diatómica para una configuración electrónica dada.

Los parámetros r0 y E0 están determinados por la estructura de las moléculas individuales. Para r>r0 la pendiente de Ep(r) es positiva, la fuerza es atractiva, para r<r0, la pendiente es negativa y la fuerza es fuertemente repulsiva, mientras que para r=r0 la fuerza es nula (mínimo de la energía potencial).

Algunos valores típicos de r0 y E0 son los siguientes

Gas E0 (10-23 J) r0 (angstroms)
Hidrógeno (H2) 43 3.3
Nitrógeno (N2) 131 4.2
Oxígeno (O2) 162 3.9

Fuente: Physics, Mechanics, Waves and Thermodynamics. Roller, Blum Edt. Holden-day (1981), pág. 655.

En la simulación vamos usar un potencial e(x) que es independiente de tipo de molécula, tomando x=r/r0.

La función e(x) se hace muy grande para valores de x próximos a cero, y tiende hacia cero para valores grandes de x. El eje X es su asíntota horizontal. Este potencial nos suministra un ejemplo más de oscilaciones anarmónicas.

Energía

La energía total de la partícula en una posición x es la suma de la energía cinética y potencial

E=Ek+Ep(x)

anarmonico1.gif (2652 bytes) En la figura el segmento AB representa la energía total E, el segmento AC la energía potencial y el segmento BC la energía cinética.

La energía potencial es máxima en los dos puntos más alejados del equilibrio x=1, y es mínima en x=1, (posición de equilibrio). Lo contrario le ocurre a la energía cinética que es mínima en los dos puntos más alejados del origen y es máxima en el mínimo x=1.

Fuerza

La fuerza sobre la partícula será

como vemos la posición de equilibrio f(x)=0 es x=1, que corresponde a un mínimo cuyo valor es -1.

A la izquierda del origen x<0, la pendiente de la función Ep es negativa y su valor es muy grande (el ángulo es un poco más de 90º), la fuerza es positiva (hacia la derecha). A la derecha del origen x>0 la pendiente de la curva es positiva y su valor es pequeño (próximo a cero para valores grandes de x), la fuerza es negativa (hacia la izquierda).

Ecuación del movimiento

La ecuación del movimiento de una partícula de masa m=1 es

Resolvemos esta ecuación por el método de Runge-Kutta para las siguientes condiciones iniciales.

En el instante t=0, la posición inicial de la partícula es x0 y su velocidad inicial es v0=0. En la posición inicial la energía cinética de la partícula es cero y por tanto, la energía potencial es igual a la energía total e.

Dada la energía total e de la partícula podemos hallar la posición inicial x0, para la cual la partícula está en reposo. Para ello, hay que resolver la ecuación.

 

Actividades

En el programa interactivo, se introduce un valor para la energía total de la partícula en el intervalo (-1, 0), en el control de edición titulado Energía

En el applet se puede observar tres zonas:

  • A la derecha, la representación de la función energía potencial e(x). Los valores de la energía cinética y potencial de la partícula representados por barras de color rojo y azul, respectivamente.
  • En la parte inferior izquierda, el movimiento de la partícula unida a un imaginario muelle elástico no ideal. Se representa mediante una flecha la fuerza sobre la partícula. Desde el punto de vista del concepto de impulso, vemos que una fuerza intensa actuando durante un corto intervalo de tiempo produce el mismo cambio de velocidad que un fuerza relativamete pequeña durante un tiempo más grande.
  • En la parte superior izquierda se representa la trayectoria de la partícula en el espacio de las fases.

En este caso, es mucho más evidente que cuando la energía total es próxima al mínimo (-1) , la partícula describe aproximadamente un M.A.S. La trayectoria en el espacio de las fases es aproximadamente una elipse. Si la energía total de la partícula es próxima a cero, se observa una trayectoria en el espacio de las fases que se asemeja a elipse fuertemente deformada. La partícula pasa la mayor parte del tiempo que emplea en completar la trayectoria en la región x>1.

ForzadasApplet aparecerá en un explorador compatible con JDK 1.1.