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
- El módulo es, F=qvB
- La dirección, perpendicular al plano formado por los vectores y
- Sentido, la regla del sacacorchos
Aplicamos la ecuación de la dinámica del movimiento crcular uniforme
Como |B2|≤|B1| entonces r2≥r1
Como vemos en la figura, el vector fuerza es perpendicular al vector velocidad , 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
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)
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
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
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
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
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