問題
z軸方向に一様な磁場→Bがある。この磁場内で運動する荷電粒子について論ぜよ。但し、荷電粒子の質量m、電荷をqとする。
解答
一様な磁場内→Bで運動する荷電粒子に作用する力はローレンツ力が作用するので運動方程式は
md→vdt=q(→v×→B)
と表される。
成分で表すと
m(dvxdtdvydtdvzdt)=q(vxvyvz)×(00B)=q(vy⋅B−vz⋅0vx⋅0−vx⋅Bvx⋅0−vy⋅0)=q(vyB−vxBo)
となる。
従って、それぞれの加速度は
dvxdt=qBmvydvydt=−qBmvxdvzdt=0
となり、連立微分方程式で表される。
x成分の式をtで微分すると
dvxdt=qBmvyddt(dvxdt)=ddt(qBmvy)d2vxdt2=qBmdvydt=qBm⋅−qBmvx=−(qBm)2vx
となる。ここでω=qBmと置くと
d2vxdt2=−ω2vx
となる。
この微分方程式の一般解は
vx(t)=αcosωt+βcosωt (α,β:任意定数)
となる。
vyについてはx成分の式より
dvxdt=ωvyvy=1ωdvxdt=1ωddt(αcosωt+βcosωt)=1ω[α(−sinωt)⋅ω+βcosωt⋅ω]=−αsinωt+βcosωt
となる。
従って、速度vx(t),vy(t)は
vx(t)=αcosωt+βcosωt vy(t)=−αsinωt+βcosωt
となる。
それぞれからx(t),y(t)を求めると
vx(t)=dx(t)dt=αcosωt+βcosωt ∫dx(t)dtdt=∫(αcosωt+βcosωt)dt∫dx(t)=∫(αcosωt+βcosωt)dtx(t)=1ωαsinωt+1ωβ(−sinωt)+C1(C1:積分定数)=αωsinωt−βωcosωt+C1
続いて
vy(t)=dy(t)dt=−αsinωt+βcosωt∫dy(t)dtdt=∫(−αsinωt+βcosωt)dt∫dy(t)=∫(−αsinωt+βcosωt)dty(t)=−1ωα(−cosωt)+1ωβsinωt+C2(C1:積分定数)=αωcosωt+βωsinωt+C2
となる。
ここでRcosθ=−bω,Rsinθ=αωとなるようにRとωを設定すると
x(t)=αωsinωt−βωcosωt+C1=Rsinθsinωt+Rcosθcosωt+C1=Rcos(θ−ωt)+C1y(t)=αωcosωt+βωsinωt+C2=Rsinθcosωt−Rcosθsinωt+C2=Rsin(θ−ωt)+C2
となる。
z軸については
dvzdt=0∫dvzdtdt=∫0dt∫dvz=∫0dtvz=C3=const.(C3:積分定数)dz(t)dt=vz∫dz(t)dtdt=∫vzdt∫dz(t)=∫vzdtz(t)=vzt+C4(C4:積分定数)
となる。t=0でz(0)=Z0とすると
z(0)=vz⋅0+C4=z0C4=z0
となるので
z(t)=vzt+z0
となる。
以上より、結果をまとめると
x(t)=Rcos(θ−ωt)+C1y(t)=Rsin(θ−ωt)+C2z(t)=vzt+z0
x−y平面では回転運動、z軸方向は等速度運動になり、これらをあわせると螺旋運動になります。
ω<0の時、円運動は反時計回りになり、ω>0の時、円運動は時計回りになります。