Description

Apply Hamiltonian separately for both times inside Hamiltonian to capture all features (ions, gauge field, phases, ...). This makes it slower because the full Hamiltonian is applied twice per operation inside the exponential, but at least all terms are treated correctly.