El espectrómetro de masas incluyendo el peso de los iones

El espectrómetro de masas sin tener en cuenta el peso del ión.

Consideremos un ión de carga +q y masa m que sale del origen con velocidad v0 en la dirección del eje Z tal como se muestra en la figura. El ión se mueve en una región del plano YZ donde hay un campo magnético uniforme B cuya dirección es el eje X.

La fuerza que ejerce el campo magnético sobre la partícula cargada es

Fm=qv×B

La ecuación del movimiento es

m d 2 r d t 2 =q( v y j ^ + v z k ^ )×B i ^ m d 2 r d t 2 =qB( v z j ^ v y k ^ )

o bien,

d 2 y d t 2 =ω dz dt ω= qB m d 2 z d t 2 =ω dy dt

Tenemos que resolver este sistema de dos ecuaciones diferenciales con las siguientes condiciones iniciales: en el instante t=0, el ión parte del origen con velocidad vz=v0

Integramos la primera ecuación diferencial

d dt ( dy dt ωz )=0 dy dt ωz= c 1

Determinamos la constante de integración c1 a partir de las condiciones iniciales: en el origen z=0, vy=dy/dt=0, resultando c1=0.

dy dt ωz=0

Integramos la segunda ecuación diferencial

d dt ( dz dt +ωy )=0 dz dt +ωy= c 2

Determinamos la constante de integración c2 a partir de las condiciones iniciales: en el origen y=0, vz=dz/dt=v0, resultando c2=v0.

dz dt +ωy= v 0

Hemos reducido un sistema de dos ecuaciones diferenciales de segundo orden a un sistema de dos ecuaciones diferenciales de primer orden

dy dt ωz=0 dz dt +ωy= v 0

Derivamos la segunda ecuación respecto de z e introducimos la expresión dy/dt de la primera

d 2 z d t 2 + ω 2 z=0

Esta es la ecuación diferencial de un Movimiento Armónico Simple, la solución es z=Asin(ωt)+Bcos (ωt). Los coeficientes A y B se determinan a partir de las condiciones iniciales: en el instante t=0, z=0, dz/dt=v0. Resultando B=0, y A=v0/ω

z= v 0 ω sin( ωt )

Derivamos respecto de z y despejamos y en la segunda ecuación diferencial

y= v 0 ω v 0 ω cos( ωt )

Estas son las ecuaciones paramétricas de la trayectoria, eliminando el tiempo t entre ambas ecuaciones.

z 2 + ( y v 0 ω ) 2 = ( v 0 ω ) 2

que es la ecuación de una circunferencia centrada en (y, z)=(v0/ω, 0) y de radio

R= v 0 ω = m v 0 qB

que es el resultado que ya obtuvimos en la página anterior.

El espectrómetro de masas teniendo en cuenta el peso del ión.

La ecuación del movimiento de la partícula es ahora

m d 2 r d t 2 =q( v y j ^ + v z k ^ )×B i ^ mg k ^ m d 2 r d t 2 =qB( v z j ^ v y k ^ )mg k ^

o bien,

d 2 y d t 2 =ω dz dt ω= qB m d 2 z d t 2 =ω dy dt g

Tenemos que resolver este sistema de ecuaciones diferenciales con las siguientes condiciones iniciales: en el instante t=0, el ión parte del origen con velocidad vz=v0

Integramos la primera ecuación diferencial

d dt ( dy dt ωz )=0 dy dt ωz= c 1

Determinamos la constante de integración c1 a partir de las condiciones iniciales: en el origen z=0, vy=dy/dt=0, resultando c1=0.

Integramos la segunda ecuación diferencial

d dt ( dz dt +ωy )=g dz dt +ωy=gt+ c 2

Determinamos la constante de integración c2 a partir de las condiciones iniciales: en el instante t=0, y=0, vz=dz/dt=v0,  resultando c2=v0.

dz dt +ωy=gt+ v 0

Hemos reducido un sistema de dos ecuaciones diferenciales de segundo orden a un sistema de dos ecuaciones diferenciales de primer orden

dy dt ωz=0 dz dt +ωy=gt+ v 0

Derivamos la segunda ecuación respecto de z y sustituimos la expresión de dy/dt de la primera ecuación.

d 2 z d t 2 + ω 2 z=g

La solución particular de esta ecuación diferencial es z1=c3. Introduciéndola en la ecuación diferencial obtenemos z1=c3=-g/ω2.

La solución de la ecuación diferencial homogénea es z2=Asin(ωt)+Bcos (ωt).

La solución completa es z=z1+z2

z=Asin( ωt )+Bcos( ωt ) g ω 2

Los coeficientes A y B se determinan a partir de las condiciones iniciales: en el instante t=0, z=0, vz=dz/dt=v0. Resultando B= g/ω2, y A=v0/ω

z= v 0 ω sin( ωt )+ g ω 2 cos( ωt ) g ω 2

Derivamos respecto de z y despejamos y en la segunda ecuación diferencial

y= g ω 2 sin( ωt ) v 0 ω cos( ωt )+ v 0 ω g ω 2 (ωt)

En un espectrómetro real g<<ω2, la trayectoria de un ión es una semicircunfencia, tal como se ha descrito en el apartado anterior.

Intersección del haz de iones con el eje Y

En la primera ecuación paramétrica ponemos z=0, y calculamos las raíces de la ecuación

v 0 ω sin( ωt )+ g ω 2 cos( ωt )= g ω 2 Asin( ωt+φ )= g ω 2 Acosφ= v 0 ω Asinφ= g ω 2 φ=arctan g ω v 0

La ecuación se escribe

Asin(ωt+φ)=Asin(φ)

o bien,

sin(ωt+φ)=sin(φ)

Una solución es ωt=2nπ, n=0, 1, 2, 3…

La solución n=0 corresponde a t=0, cuando el ión sale del origen

La otra solución es

ωt=(2n+1)π-2φ.

Cuando n=0, tenemos ωt=π-2φ que corresponde a la primera intersección del haz de iones con el eje Y, tal como se señala en la figura mediante un círculo de color azul.

La coordenada y del primer punto de impacto es

y=Asinφsin(ωt)Acosφcos(ωt)+ v 0 ω g ω 2 (ωt)= Acos(ωt+φ)+ v 0 ω g ω 2 (ωt)

En el instante t tal que ωt=π-2φ

y 1 =Acos(πφ)+ v 0 ω g ω 2 (π2φ)=2 v 0 ω g ω 2 ( π2arctan g v 0 ω )

Ecuaciones paramétricas de la trayectoria

Escribimos las ecuaciones paramétricas de la trayectoria en términos de dos parámetros

α= g ω 2 β= v 0 ω z=βsin( ωt )+αcos( ωt )α y=αsin( ωt )βcos( ωt )+βα(ωt)

Vamos dando valores a ωt para dibujar la trayectoria.

wt=linspace(0,20,200);
alfa=0.1;
beta=1;
z=beta*sin(wt)+alfa*cos(wt)-alfa;
y=alfa*sin(wt)-beta*cos(wt)+beta-alfa*wt;
plot(y,z)
axis equal
grid on
xlabel('y')
ylabel('z')
title('Espectrómetro de masas')

Los puntos de intersección con el eje Y, es decir, z=0 se calculan del siguiente modo

ωt=π2arctan α β y 1 =2βα( π2arctan α β ) ωt=2π y 2 =2πα ωt=3π2arctan α β y 3 =2βα( 3π2arctan α β ) ωt=4π y 4 =4πα

y así, sucesivamente

Si no se tuvieran en cuenta el peso de los iones g=0, α=0 e y2n+1=2β, y2n=0, con n=0.1,2,3…

Actividades

Se introduce

Se pulsa el botón titulado Nuevo. Utilice los botones pausa || y paso a paso >| para acercarse a la intersección de la trayectoria del ión con el eje Y

Se observa el movimiento del ión “pesado”

Se mide las intersecciones con el eje Y, es decir, cuando z=0, y se comprueban con las calculadas.

Ejemplo.

Sea α=0.1 y β=1.0

ωt yi
2.94 1.71
6.28 -0.63
9.22 1.08
12.57 -1.26

Referencias

Gondar J. L., Costa B. The mass spectrograph revisited. Eur. J. Phys. 20 (1999) 469-475