The acceleration is given by

dE = - e*(VOLT+V1 x+V20 x^2/2+V11 x y+V02 y^2/2)
*Sin[2 Pi FREQ (t-ts) + PHI +DPHI],

とありますが，磁場の影響はSAD中に入っているのでしょうか？

上式は進行方向電場に関する方程式のようなので，ベクトルポテンシャルは
Az=(VOLT+V1 x+V20 x^2/2+V11 x y+V02 y^2/2)
*Cos[2 Pi FREQ (t-ts) + PHI +DPHI]/2 Pi FREQ,
になると思います．
この場合，
md^2x/dt^2=-e*v_z*(V1+V20 x+V11 y)*Cos[2 Pi FREQ (t-ts) + PHI +DPHI]/2 Pi FREQ,
md^2y/dt^2=-e*v_z*(V02 y+V11 y)*Cos[2 Pi FREQ (t-ts) + PHI +DPHI]/2 Pi FREQ,
といった効果です．

ご指摘の効果ははいっています。そうでなければsymplecticになりません。


> The acceleration is given by
>
> dE = - e*(VOLT+V1 x+V20 x^2/2+V11 x y+V02 y^2/2)
> *Sin[2 Pi FREQ (t-ts) + PHI +DPHI],
>
> とありますが，磁場の影響はSAD中に入っているのでしょうか？
>
> 上式は進行方向電場に関する方程式のようなので，ベクトルポテンシャルは
> Az=(VOLT+V1 x+V20 x^2/2+V11 x y+V02 y^2/2)
> *Cos[2 Pi FREQ (t-ts) + PHI +DPHI]/2 Pi FREQ,
> になると思います．
> この場合，
> md^2x/dt^2=-e*v_z*(V1+V20 x+V11 y)*Cos[2 Pi FREQ (t-ts) + PHI +DPHI]/2 Pi FREQ,
> md^2y/dt^2=-e*v_z*(V02 y+V11 y)*Cos[2 Pi FREQ (t-ts) + PHI +DPHI]/2 Pi FREQ,
> といった効果です．