Dispersión de partículas alfa por un núcleo

La ley de la Gravitación Universal describe la interacción entre cuerpos debido a su masa. La fuerza de atracción entre dos cuerpos es central y conservativa, su módulo es inversamente proporcional al cuadrado de la distancia que separa los centros de ambos cuerpos. Cuando se integra la ecuación diferencial que describe el movimiento de un cuerpo bajo la acción de dicha fuerza, se obtiene una trayectoria que es una cónica. El tipo de cónica depende signo de la energía total del cuerpo.

Trayectoria Energía
Elipse E<0
Parábola E=0
Hipérbola E>0

Los planetas describen elipses estando el Sol en uno de sus focos. El hecho de que la energía sea negativa se debe a que la energía potencial de una fuerza atractiva es negativa, y la energía cinética es menor que la energía potencial (el cuerpo está confinado).

La interacción eléctrica puede ser repulsiva o atractiva según que las cargas sean del mismo o distinto signo. La fuerza que describe la interacción eléctrica es central y conservativa, su módulo, de acuerdo a ley de Coulomb, es inversamente proporcional al cuadrado de la distancia que separa ambas cargas.

Trayectoria de la partícula

La fuerza de repulsión entre dos cargas Q y q del mismo signo es de acuerdo a la ley de Coulomb

F= 1 4π ε 0 qQ r 2 r ^

Esta fuerza es conservativa, y la energía potencial Ep correspondiente es

E p = 1 4π ε 0 qQ r

Supongamos que Q es una carga fija y que una partícula de masa m y carga q se mueve en el campo creado por la carga Q.

Como la fuerza de repulsión es central y conservativa se cumple

  1. La energía total de la partícula cargada es constante
  2. E= 1 2 m v 2 + k r k= qQ 4π ε 0

  3. El momento angular es constante
  4. L=r×mv

    Que L sea constante en dirección y sentido quiere decir que la trayectoria de la partícula está contenida en un plano perpendicular a la dirección del momento angular.

Expresamos la energía y el momento angular en coordenadas polares

E= 1 2 m ( dr dt ) 2 + 1 2 m r 2 ( dθ dt ) 2 + k r L=m r 2 dθ dt

Las ecuaciones de constancia del momento angular y de la energía constituyen un par de ecuaciones diferenciales en las que se puede eliminar el tiempo t. Para obtener la ecuación de la trayectoria r=r(θ) se integra la ecuación diferencial

dθ dr = L m r 2 2 m ( E L 2 2m r 2 k r )

El resultado es una hipérbola

r= d εcosθ1 ε= 1+ 2E L 2 m k 2 d= L 2 mk

Las condiciones iniciales son las de una partícula de masa m y carga q que se mueve a gran distancia del centro de fuerzas con velocidad v0, tal como se muestra en la figura.

Parámetro de impacto

El parámetro de impacto b es la distancia existente entre la dirección de la partícula incidente, cuando se encuentra muy alejada del centro de fuerzas y el centro de fuerzas. En la figura, el parámetro de impacto b es la distancia entre la dirección de la velocidad v0 y la carga fija Q.

El módulo del momento angular es

L=mv0·b

La energía total

E= 1 2 m v 0 2

Los valores de d y la excentricidad ε en la ecuación de la hipérbola expresados en términos de la energía E y del parámetro de impacto b son, respectivamente.

d = 2 E b 2 k ε = 1 + ( 2 E b k ) 2 k = q Q 4 π ε 0

En la figura, se representa la ecuación de la trayectoria, una hipérbola, que tiene dos asíntotas, dos rectas simétricas con respecto al eje X, que forman un ángulo θ0 con éste, cuyo valor es

ε ·cosθ0-1=0 o bien, cos θ 0 = 1 ε

A continuación, giramos la hipérbola en sentido antihorario un ángulo θ0. La ecuación de la trayectoria en coordenadas polares es

r= d εcos(θ θ 0 )1

En la figura, se representa la trayectoria de la partícula cargada y en la que se ha señalado, el parámetro de impacto b y el ángulo de dispersión Φ, que es el ángulo entre las direcciones inicial (asíntota horizontal) y final (asíntota inclinada) de la velocidad de la partícula

Ángulo de dispersión

Cuando la partícula se aleja mucho del centro de fuerzas, sigue una trayectoria que tiende asintóticamente a una línea recta. El ángulo que forma dicha recta con el eje horizontal se denomina ángulo de dispersión.

El ángulo de dispersión vale Φ =180-2θ0, es decir

Φ=π2arccos 1 1+ ( 2Eb k ) 2 tan Φ 2 = k 2Eb

Ejemplo:

Hemos establecido un sistema de unidades en el que k=1.

El ángulo de dispersión para E=0.5 y b=0.75 vale Φ =106.3º, como comprobaremos en el programa interactivo, más abajo.

Representamos el ángulo de dispersión Φ en función de parámetro de impacto, para una energía E determinada

k=1;
E=0.5;
phi=@(b) 2*atan(k./(2*E*b))*180/pi;
fplot(phi,[0,5])
grid on
xlabel('b')
ylabel('\Phi')
title('Relación parámetro de impacto, ángulo de dispersión')

Actividades

Se introduce

Se pulsa el botón titulado Nuevo

Se representa la trayectoria de la partícula y se calcula el ángulo de dispersión Φ


Utilizando el principio de conservación de la energía calculamos la distancia mínima de aproximación de una partícula cargada, que choca de frente contra un núcleo atómico

Para hacer más simple el problema supondremos que la masa del núcleo es mucho mayor que la masa del proyectil.

Si la carga del núcleo es Q y la del proyectil es q. La energía total del proyectil es

E= 1 2 m v 2 + Qq 4π ε 0 r = 1 2 m v 0 2

Cuando el proyectil está a mucha distancia del núcleo, su velocidad es v0 y toda la energía es cinética. En el punto C de máximo acercamiento (véase la figura), la velocidad v es transversal (perpendicular a la dirección radial, dr/dt=0) de modo que el momento angular es L=mRv. La ecuación de la conservación de la energía en dicho punto de máximo acercamiento se escribe

L 2 2m R 2 + Qq 4π ε 0 R = 1 2 m v 0 2

Ecuación de segundo grado en 1/R que permite obtener R en función de la energía y del momento angular de la partícula.

Para una colisión de frente, L=0 y se despeja R

R= qQ 4π ε 0 ( 1 2 m v 0 2 )

En una colisión frontal, la velocidad de la partícula en el punto de máximo acercamiento es cero, v=0

Relación entre parámetro de impacto y ángulo de dispersión

En el apartado anterior, hemos obtenido la relación cuantitativa entre el parámetro de impacto y el ángulo de dispersión para una energía E de la partícula cargada. Mediante el programa interactivo que viene a continuación, se puede establecer una relación cualitativa, observando las trayectorias de un número de partículas de la misma energía que inciden sobre el centro fijo de fuerzas con parámetros de impacto espaciados regularmente.

Se introduce

Se pulsa el botón titulado Empieza

Se observa las trayectorias de las partículas que inciden sobre el centro fijo de fuerzas

La envolvente

Escribimos la ecuación de la trayectoria de la forma

1 r = ε d cos(θ θ 0 ) 1 d 1 r = ε d ( cosθ·cos θ 0 +sinθ·sin θ 0 ) 1 d 1 r = ε d cos θ 0 ( cosθ+sinθ·tan θ 0 ) 1 d

Teniendo en cuenta las relaciones

cos θ 0 = 1 ε ε= 1+ ( 2Eb k ) 2 d= 2E b 2 k cos θ 0 = 1 1+ tan 2 θ 0 tan θ 0 = 2Eb k

La ecuación de la trayectoria se escribe en términos del parámetro b

1 r = k 2E b 2 (cosθ1)+ sinθ b

La ecuación de la trayectoria depende del parámetro de impacto b,
f(r, θ, b)=0

La ecuación de la envolvente de las trayectorias se obtiene derivando con respecto a b e igualando a cero.

f b =0 k E b 3 (cosθ1)+ sinθ b 2 =0 b= (1cosθ) sinθ k E

y combinando ésta con la ecuación de la trayectoria para eliminar el parámetro b. Es decir, se introduce la expresión de b en la ecuación de la trayectoria.

1 r = E 2k (1+cosθ)

Esta es la ecuación de la envolvente en coordenadas polares. Para escribirla en coordenadas rectangulares, ponemos
x=r·cosθ, y=r·sinθ

2k E = x 2 + y 2 +x y 2 = 4k E ( k E x )

En este último programa interactivo, se dibuja la envolvente (en color azul) de las trayectorias hiperbólicas que describen las partículas

Trayectorias

Las ecuaciones del movimiento de la partícula de masa m y carga q, en el campo eléctrico producido por la carga fija Q son:

m d 2 x d t 2 = 1 4π ε 0 qQ r 2 x r m d 2 y d t 2 = 1 4π ε 0 qQ r 2 y r

Resolvemos el sistema de ecuaciones diferenciales por procedimientos numéricos. Superponemos en líneas a trazos, la curva de energía potencial Ep(r)

E p (r)= 1 4π ε 0 qQ r

Para una fuerza repulsiva, k=1, el resultado es el que se ha visto en los programas interactivos, más arriba. Mejores resultados se obtienen haciendo que la partícula cargada parta muy lejos del origen (teóricamente de -∞) por ejemplo x=-100, pero no hay mucha diferencia desde el punto de vista cualitativo.

k=1;
energia=1;
N=10; %partículas
tspan=[0,10];
hold on
potencial=@(r) k./r;
for i=1:N
    x0=[-10,sqrt(2*energia),2.5*i/N,0];
    fg=@(t,x)[x(2);k*x(1)/(x(1)^2+x(3)^2)^(3/2); x(4);
k*x(3)/(x(1)^2+x(3)^2)^(3/2)];
    [t,x]=ode45(fg,tspan,x0);
    plot(x(:,1),x(:,3)) 
end
fp=fplot(potencial,[0.05,5],'--k');
plot(-fp.XData,fp.YData,'--k')
ylim([-3,3]);
xlim([-4,4])
hold off
grid on
xlabel('x')
ylabel('y');
title('dispersión')

Para una fuerza atractiva, k=-1, el resultado es el siguiente

k=-1;
energia=1;
N=10; %partículas
tspan=[0,10];
hold on
potencial=@(r) k./r;
for i=1:N
    x0=[-10,sqrt(2*energia),2.5*i/N,0];
    fg=@(t,x)[x(2);k*x(1)/(x(1)^2+x(3)^2)^(3/2); x(4); 
k*x(3)/(x(1)^2+x(3)^2)^(3/2)];
    [t,x]=ode45(fg,tspan,x0);
    plot(x(:,1),x(:,3)) 
end
fp=fplot(potencial,[0.05,5],'--k');
plot(-fp.XData,fp.YData,'--k')
ylim([-3,3]);
xlim([-4,4])
hold off
grid on
xlabel('x')
ylabel('y');
title('dispersión')

Referencias

French A. P. The envelopes of some families of fixed-energy trajectories. Am. J. Phys. 61 (9) September 1993