“…7,8,16,126 For high-dimensional nuclear wave functions, the method of choice is the Multiconfigurational Time-Dependent Hartree (MCTDH) algorithm, 127,128 possibly in its multilayer implementation. 129 When potential surfaces can be approximated as harmonic oscillators, tensor network approaches are another powerful way to perform full quantum dynamics 130,131 As a hallmark feature, methods relying on wavepacket propagations guarantee a correct dynamics of the nuclear wavepacket at both electronic and polaritonic avoided crossings, conical intersections, and seams between the PoPESs, including a correct decay of nuclear coherence without needing to resort to artificial corrections as in the semiclassical methods. Secondly, its propagation allows to exactly include decay channels in the dynamics, either through effective non-Hermitian Hamiltonians 11,12,83 that are exact when the dynamics after decay are not of interest, or by direct solution of a Lindblad-style master equation.…”