The implicit midpoint method described by Chambliss and Franklin [Am. J. Phys. **88**, 1075 (2020)] is known not to exactly conserve energy for a spatially varying magnetic field, nor does it conserve angular momentum in a constant magnetic field. It is, therefore, not competitive with the widely used Boris solver in plasma physics, which conserves both.

The implicit midpoint algorithm

where $r\u2032=r(t+\Delta t),\u2009v\u2032=v(t+\Delta t)$, and $r=r(t),\u2009\u200av=v(t)$, corresponding to Chambliss and Franklin's^{1} Eqs. (5) and (6), takes the form of

for a charged particle *q* in a general magnetic field $B(r)=B(r)B\u0302(r)$, with cyclotron frequency $\omega (r)=qB(r)/m$. This algorithm is known *not* to conserve energy. As noted by Boris,^{2} dotting both sides of Eq. (3) with $(v\u2032+v)$ gives

In a spatially varying magnetic field, there is no reason why the triple product on the RHS of Eq. (4) should vanish, and therefore, the kinetic energy is generally not conserved. However, for a constant magnetic field where $\omega \u2032=\omega $ and $B\u0302\u2032=B\u0302$, it would vanish,

which is why Chambliss and Franklin's Fig. 1 showed zero energy error.

This energy conserving feature can be generalized to the case of a nonconstant magnetic field, since as long as the same magnetic field is used to evaluate the acceleration at **r′** and **r**, the RHS of Eq. (4) will still vanish as Eq. (5). This suggests that one should use the common magnetic field at the midpoint, resulting in the Boris solver^{2}

widely used in plasma physics.^{3}

Also, it should be further noted that, in a constant magnetic field, both velocity vectors (3) and (6) are rotated the same way. However, the Boris position (7) will remain on the gyrocircle orbit of radius $rg=v/\omega $ regardless of the value of $\Delta t$. As a consequence, it also exactly conserves angular momentum, $r\u2032\xd7v\u2032=r\xd7v$, at all $\Delta t$. This is because for motion perpendicular to $B\u0302$, since

the constancy of *r*^{2} implies the conservation of angular momentum. The implicit midpoint position (2), on the other hand, will yield positions on a *larger* circle tangent to the gyrocircle at the initial position. As a consequence, its angular momentum $r(t)\xd7v(t)$ will oscillate sinusoidally *above* the initial value as a function of time. This is clearly shown in Chambliss and Franklin's Fig. 2. In contrast, Boris's algorithm would have produced zero error to machine precision.