“…Since the introduction of the tilde space doubles the number of nuclear degrees of freedom (DoF), and since a thermal environment can be realistically mimicked only using hundreds of DoFs, especially if the temperature is very high, as in the case of a thermal radiation field, it is essential to use a methodology suitable to treat a large number of dynamical variables. Several techniques have been developed which can, at least in principle, overcome what has been termed the curse of dimensionality . , Several methodologies can nowadays be used to accurately describe the quantum dynamics of a system with a large number of DoFs, among which Multi-Configuration Time-Dependent Hartree (MCTDH), Multiple-Davydov’s Ansätze (MDA), hierarchical equations of motion (HEOM), ,, Quasi-adiabatic path integrals (QUAPI), and the wide family of Tensor-Network representations , are probably the most representative. Here, we employ the so-called Tensor Train (TT) format (or Matrix Product States (MPS) in the physics literature) to efficiently represent the vibronic wave function.…”