Implement also higher order split operator schemes, such as Yoshida-Neri, Kahan-Li, and Blanes-Casa (see also their web site). From the many works of the latter ones, in particular the schemes for time-dependent Hamiltonians should be of interest (because for time-independent Hamiltions we already have the Chebychev polynomial scheme which is considered excellent).
Is the splitting between potential energy couplings and dipole couplings (as currently implemented here and there) really necessary? After all, they should commute with each other.
To be discussed with U.L.
Last edit: Burkhard Schmidt 2023-07-21
Diff: