Iterative solution for particle in magnetic field

Question

I have some questions relating movement of a particle calculating with mathematica.

A particle (mass: m, charge: q) moves on x-axis from left with speed $v_0$ to a magnetic field $\vec{B} = (0,0,B)$. Entry is at $\vec{r}(t=0)=(0,0,0)$.

1. How would equation of movement look like for this problem ?

I tried this:

In[23]:= r[t_] := {x[t], y[t], z[t]}

In[24]:= electricField = 0, Efield, 0;

In[24]:= magneticField = 0 {0, 0, Bfield};

In[26]:= force = q (electricField + Cross[r'[t], magneticField]);

In[27]:= DSolve[{m r''[t] == force, r[0] == {0, 0, 0}, 
 r[0] == {v1, v2, v3}}, r[t], t] //. {Bfield -> m \[Omega] / q, 
Efield -> V Bfield} // ExpandAll // Simplify

During evaluation of In[27]:= DSolve::bvnul: For some branches of the general solution, the given boundary conditions lead to an empty solution.

Out[27]= {}

In[16]:= Clear[eqns]

In[29]:= eqns = Map[Thread, {m r''[t] == force, r[0] == {0, 0, 0}, 
  r'[0] == {v1, v2, v3}}] //. {Bfield -> m \[Omega] /q, 
 Efield -> V Bfield} // ExpandAll // Simplify

Out[29]= {{electricField q == m (x^\[Prime]\[Prime])[t], 
electricField q == m (y^\[Prime]\[Prime])[t], 
electricField q == m (z^\[Prime]\[Prime])[t]}, {x[0] == 0, 
y[0] == 0, z[0] == 0}, {v1 == Derivative[1][x][0], 
v2 == Derivative[1][y][0], v3 == Derivative[1][z][0]}}

In[31]:= solution1 = 
Dsolve[eqns, {x[t], y[t], z[t]}, t] // ExpandAll // Simplify

Out[31]= Dsolve[{{electricField q == m (x^\[Prime]\[Prime])[t], 
electricField q == m (y^\[Prime]\[Prime])[t], 
electricField q == m (z^\[Prime]\[Prime])[t]}, {x[0] == 0, 
y[0] == 0, z[0] == 0}, {v1 == Derivative[1][x][0], 
v2 == Derivative[1][y][0], v3 == Derivative[1][z][0]}}, {x[t], 
y[t], z[t]}, t]

2. When we start with $\vec{r}^{(0)}(t)=q(v_0t,0,0)$ as zeroth approximation and insert it in Lorentz-force, how would equation of movement look like ?

How can I write a routine that automatises this iteration till 6 iterations ? I read about "Nest" command but how to use it ?

3. Charged particles describe circuits in homogenous magnetic field. How can I find exact solution with mathematica and see what angular frequency and radius they have ? There should be a possibility to develop exact solutions in power of $B$, how tovisualise it with mathematica ?

Answer 1

Change In[27] to

DSolve[{m r''[t] == force, r[0] == {0, 0, 0},r'[0] == {v0, 0, 0}} //. {Bfield -> m \[Omega]/q,Efield -> V Bfield}, r[t], t]
(*{{x[t] -> (t V \[Omega] - V Sin[t \[Omega]] + v0 Sin[t \[Omega]])/\[Omega], 
y[t] -> (1/\[Omega])(-v0 - V Cos[t \[Omega]] + v0 Cos[t \[Omega]] +V Cos[t \[Omega]]^2 + V Sin[t \[Omega]]^2), z[t] -> 0}}*)

and MMA finds a unique solution.