Movimiento de una partícula cargada en un campo no uniforme

La fuerza que ejerce el campo magnético B, paralelo al eje Z, sobre una partícula de masa m, cargada q>0 (positiva), con velocidad v en el plano XY es

F → =q v → × B →

Aplicamos la ecuación de la dinámica del movimiento crcular uniforme

qvB=m v 2 r r= mv qB

Como |B2|≤|B1| entonces r2≥r1

Como vemos en la figura, el vector fuerza F → es perpendicular al vector velocidad v → , el vector fuerza apunta al centro de la circunferencia de radio r.

Esto es lo que necesitamos saber, para trazar la trayectoria de la partícula cargada que entra en la región con campo B1 con velocidad v haciendo un ángulo θ con el eje X

B2 es positivo y B1 es negativo

El caso más sencillo se produce cuando |B2|=|B1|=B

q=1.6e-19; %carga
m=1.67e-27; %masa
B=100e-4; %campo magnético
v=2e5; %velocidad
r=m*v/(q*B); %radio
alfa=pi/3;
x=@(th) r*sin(alfa)+r*cos(th);
y=@(th) r*cos(alfa)+r*sin(th);
hold on
fplot(x,y,[0,2*pi])
plot(r*sin(alfa), r*cos(alfa),'ko','markersize',2,'markerfacecolor','k')
fplot(x,y,[-pi/2-alfa,-pi/2+alfa],'color','r')
x=@(th) 3*r*sin(alfa)+r*cos(th);
y=@(th) -r*cos(alfa)+r*sin(th);
fplot(x,y,[pi/2-alfa,pi/2+alfa],'color','r')
x=@(th) 5*r*sin(alfa)+r*cos(th);
y=@(th) r*cos(alfa)+r*sin(th);
fplot(x,y,[-pi/2-alfa,-pi/2+alfa],'color','r')
hold off
axis equal
axis off

La partícula se desplaza horizontalmente, una longitud L (en color verde), en el tiempo T que tarda en recorrer el arco de cincunferencia de radio r con velocidad constante v

L=2rsinθ T= 2rθ v

Se denomina velocidad de deriva vd al cociente L/T, donde L es el desplazamiento a lo largo del eje X (la línea que separa las dos regiones)

v d = L T =v sinθ θ

Examinamos el caso general |B2|<|B1|

q=1.6e-19; %carga
m=1.67e-27; %masa
B1=400e-4; %campo magnético
B2=100e-4; %campo magnético
v=2e5; %velocidad
r1=m*v/(q*B1); %radio
r2=m*v/(q*B2); %
alfa=pi/3;
x=@(th) r1*sin(alfa)+r1*cos(th);
y=@(th) r1*cos(alfa)+r1*sin(th);
hold on
fplot(x,y,[0,2*pi])
plot(r1*sin(alfa), r1*cos(alfa),'ko','markersize',2,'markerfacecolor','k')
fplot(x,y,[-pi/2-alfa,-pi/2+alfa],'color','r')
x=@(th) (2*r1+r2)*sin(alfa)+r2*cos(th);
y=@(th) -r2*cos(alfa)+r2*sin(th);
fplot(x,y,[pi/2-alfa,pi/2+alfa],'color','r')
x=@(th) (3*r1+2*r2)*sin(alfa)+r1*cos(th);
y=@(th) r1*cos(alfa)+r1*sin(th);
fplot(x,y,[-pi/2-alfa,-pi/2+alfa],'color','r')
plot((2*r1+r2)*sin(alfa), -r2*cos(alfa),'ko','markersize',2,
'markerfacecolor','k')
hold off
axis equal
axis off

La partícula se desplaza horizontalmente, una longitud L (en color verde), en el tiempo T que tarda en recorrer los arcos de circunferencia de radios r1 y r2 con velocidad constante v

L=2 r 1 sinθ+2 r 2 sinθ T= r 1 ·2θ v + r 2 ·2θ v v d = L T =v sinθ θ

La velocidad de deriva vd es hacia la derecha

B2 es negativo y B1 es positivo

q=1.6e-19; %carga
m=1.67e-27; %masa
B1=400e-4; %campo magnético
B2=100e-4; %campo magnético
v=2e5; %velocidad
r1=m*v/(q*B1); %radio
r2=m*v/(q*B2); %
alfa=pi/3;
x=@(th) r1*sin(alfa)+r1*cos(th);
y=@(th) -r1*cos(alfa)+r1*sin(th);
hold on
fplot(x,y,[0,2*pi])
plot(r1*sin(alfa), -r1*cos(alfa),'ko','markersize',2,'markerfacecolor','k')
fplot(x,y,[pi/2+alfa, 5*pi/2-alfa],'color','r')
x=@(th) -r2*sin(alfa)+r2*cos(th);
y=@(th) r2*cos(alfa)+r2*sin(th);
fplot(x,y,[-pi/2+alfa, 3*pi/2-alfa],'color','r')
plot(-r2*sin(alfa), r2*cos(alfa),'ko','markersize',2,'markerfacecolor','k')
x=@(th) (-2*r2-r1)*sin(alfa)+r1*cos(th);
y=@(th) -r1*cos(alfa)+r1*sin(th);
fplot(x,y,[pi/2+alfa, 5*pi/2-alfa],'color','r')
hold off
axis equal
axis off

La partícula se desplaza horizontalmente, una longitud L (en color verde), en el tiempo T que tarda en recorrer los arcos de circunferencia de radios r1 y r2 con velocidad constante v

L=2 r 1 sinθ+2 r 2 sinθ T= r 1 ( 2π−2θ ) v + r 2 ( 2π−2θ ) v v d = L T =v sinθ π−θ

La velocidad de deriva vd es hacia la izquierda

B2 y B1 son positivos

q=1.6e-19; %carga
m=1.67e-27; %masa
B1=400e-4; %campo magnético
B2=100e-4; %campo magnético
v=2e5; %velocidad
r1=m*v/(q*B1); %radio
r2=m*v/(q*B2); %radio
alfa=pi/3;
x=@(th) r1*sin(alfa)+r1*cos(th);
y=@(th) -r1*cos(alfa)+r1*sin(th);
hold on
fplot(x,y,[0,2*pi])
plot(r1*sin(alfa), -r1*cos(alfa),'ko','markersize',2,'markerfacecolor','k')
fplot(x,y,[pi/2+alfa, 5*pi/2-alfa],'color','r')
x=@(th) r2*sin(alfa)+r2*cos(th);
y=@(th) -r2*cos(alfa)+r2*sin(th);
fplot(x,y,[pi/2-alfa, pi/2+alfa],'color','r')
plot(r2*sin(alfa), -r2*cos(alfa),'ko','markersize',2,'markerfacecolor','k')
x=@(th) (2*r2-r1)*sin(alfa)+r1*cos(th);
y=@(th) -r1*cos(alfa)+r1*sin(th);
fplot(x,y,[pi/2+alfa, 5*pi/2-alfa],'color','r')
x=@(th) (3*r2-2*r1)*sin(alfa)+r2*cos(th);
y=@(th) -r2*cos(alfa)+r2*sin(th);
fplot(x,y,[pi/2-alfa, pi/2+alfa],'color','r')
hold off
axis equal
axis off

La partícula se desplaza horizontalmente, una longitud L (en color verde), en el tiempo T que tarda en recorrer los arcos de circunferencia de radios r1 y r2 con velocidad constante v

L=2 r 2 sinθ−2 r 1 sinθ T= r 1 ( 2π−2θ ) v + r 2 ·2θ v v d = L T = 2( r 2 − r 1 )sinθ 2 r 1 ( π−θ ) v + 2 r 2 θ v =v r 2 − r 1 ( r 2 − r 1 )θ+ r 1 π sinθ v d =v 1−k ( 1−k )θ+kπ sinθ, k= r 1 r 2 = B 2 B 1

La velocidad de deriva vd es hacia la derecha, depende del ángulo θ y del cociente k entre los dos campos mágnéticos

Representamos vd/v en función del ángulo θ para tres valores de k=0.1, 0.3 y 0.5

hold on
for k=0.1:0.2:0.5
    f=@(th) sin(th)*(1-k)./(th*(1-k)+k*pi);
    fplot(f,[0,pi/2],'displayName',num2str(k))
end
hold off
set(gca,'XTick',0:pi/12:pi/2)
set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'})
grid on
legend('-DynamicLegend','location','best')
xlabel('\theta')
ylabel('v/v_d')
title('Velocidad deriva')

A medida que k se incrementa, vd alcanza un máximo para cierto valor del ángulo θ

B2 y B1 son negativos

q=1.6e-19; %carga
m=1.67e-27; %masa
B1=200e-4; %campo magnético
B2=100e-4; %campo magnético
v=2e5; %velocidad
r1=m*v/(q*B1); %radio
r2=m*v/(q*B2); %radio
alfa=pi/3;
x=@(th) r1*sin(alfa)+r1*cos(th);
y=@(th) r1*cos(alfa)+r1*sin(th);
hold on
fplot(x,y,[0,2*pi])
plot(r1*sin(alfa), r1*cos(alfa),'ko','markersize',2,'markerfacecolor','k')
fplot(x,y,[-pi/2-alfa, -pi/2+alfa],'color','r')
x=@(th) (2*r1-r2)*sin(alfa)+r2*cos(th);
y=@(th) r2*cos(alfa)+r2*sin(th);
fplot(x,y,[-pi/2+alfa, 3*pi/2-alfa],'color','r')
plot((2*r1-r2)*sin(alfa), r2*cos(alfa),'ko','markersize',2,
'markerfacecolor','k')
x=@(th) (3*r1-2*r2)*sin(alfa)+r1*cos(th);
y=@(th) r1*cos(alfa)+r1*sin(th);
fplot(x,y,[-pi/2-alfa, -pi/2+alfa],'color','r')
x=@(th) (4*r1-3*r2)*sin(alfa)+r2*cos(th);
y=@(th) r2*cos(alfa)+r2*sin(th);
fplot(x,y,[-pi/2+alfa, 3*pi/2-alfa],'color','r')
hold off
axis equal
axis off

La partícula se desplaza horizontalmente, una longitud L (en color verde), en el tiempo T que tarda en recorrer los arcos de circunferencia de radios r1 y r2 con velocidad constante v

L=2 r 2 sinθ−2 r 1 sinθ T= r 1 ·2θ v + r 2 ( 2π−2θ ) v v d = L T = 2( r 2 − r 1 )sinθ 2 r 1 θ v + 2 r 2 ( π−θ ) v =v r 2 − r 1 r 2 π−( r 2 − r 1 )θ sinθ v d =v 1−k π−( 1−k )θ sinθ, k= r 1 r 2 = B 2 B 1

La velocidad de deriva vd es hacia la izquierda

Representamos vd/v en función del ángulo θ para tres valores de k=0.1, 0.3 y 0.5

hold on
for k=0.1:0.2:0.5
    f=@(th) sin(th)*(1-k)./(pi-(1-k)*th);
    fplot(f,[0,pi/2],'displayName',num2str(k))
end
hold off
set(gca,'XTick',0:pi/12:pi/2)
set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'})
grid on
legend('-DynamicLegend','location','best')
xlabel('\theta')
ylabel('v/v_d')
title('Velocidad deriva')

vd crece con el ángulo θ

Referencias

Vladimir Ivchenko. Motion of charged particles in a stepped magnetic field. Phys. Educ. 61 (2026) 065007