sarkas.time_evolution.integrators.Integrator.magnetic_pos_verlet
sarkas.time_evolution.integrators.Integrator.magnetic_pos_verlet#
- Integrator.magnetic_pos_verlet(ptcls)[source]#
Update particles’ class based on position verlet method in the case of an arbitrary direction of the constant magnetic field. For more info see eq. (79) of Ref. [Chin, 2008]
- Parameters
ptcls (
sarkas.particles.Particles) – Particles data.- Returns
potential_energy (float) – Total potential energy.
Notes
[Chin, 2008] equations are written for a negative charge. This allows him to write \(\dot{\mathbf v} = \omega_c \hat{B} \times \mathbf v\). In the case of positive charges we will have \(\dot{\mathbf v} = - \omega_c \hat{B} \times \mathbf v\). Hence the reason of the different signs in the formulas below compared to Chin’s.
Warning
This integrator is valid for a magnetic field in an arbitrary direction. However, while the integrator works for an arbitrary direction, methods in
sarkas.tools.observableswork only for a magnetic field in the \(z\) - direction. Hence, if you choose to use this integrator remember to change your physical observables.