Finally, let us take a look at the case of the pyrazine model. The populations, coherences and total energy deviation are shown in Figure . As in the previous case for the butatriene cation, all methods underestimate the first initial population decay, with SH-EDC outpeforming CTMQC and CTMQC-E with a population behaviour similar to the QD results. In terms of coherences, SH-EDC simulates the decoherence and recoherence closely, and CT-MQC captures the decoherence well but shows a more gradual recoherence than the exact. CTMQC-E, on the other hand, remains overcoherent. The increase in the total energy of the system is notably reduced with the accumulated force correction of CTMQC-E as evidenced in the lower panel.