This review provides the fundamental theoretical tools for the development of a complete wave‐function formalism for the study of time‐evolution of chemico‐physical systems at finite temperature. The methodology is based on the non‐equilibrium thermo‐field dynamics (NE‐TFD) representation of quantum mechanics, which is alternative to the commonly used density matrix representation. TFD concepts are extended and integrated with the tensor‐train (TT) numerical tools leading to a novel and powerful theoretical and computational framework for the study of complex quantum dynamical problems. In addition, NE‐TFD techniques are extended to enable the study of dissipative open systems via a new formulation of the hierarchical equations of motion (HEOM) fully integrated with TT methodologies. We demonstrate that the combination of the TFD machinery with computational advantages of TTs results in a powerful theoretical and computational framework for scrutinizing dynamics of complex multidimensional electron‐vibrational systems. We illustrate the validity and the computational advantages of the developed methodologies by applying them to the study of quantum coherence effects in the energy‐transfer processes in antenna systems, to the analysis of fingerprints of vibrational modes in electron‐transfer and charge‐transfer processes in various model and realistic multidimensional molecular systems, as well as to simulation of other fundamental models of physical chemistry.This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and Kinetics Theoretical and Physical Chemistry > Statistical Mechanics

By employing the time-dependent variational principle and the versatile multi-D₂ Davydov trial states, in combination with the Green’s function method, we study the dynamics of the Tavis–Cummings model and the Holstein–Tavis–Cummings model in the presence of diagonal disorder and cavity–qubit coupling disorder. For the Tavis–Cummings model, time evolution of the photon population, the optical absorption spectra, and the hetero-entanglement between the qubits and the cavity mode are calculated by using the Green’s function method to corroborate numerically exact results of Davydov’s Ansätze. For the Holstein–Tavis–Cummings model, only the latter is utilized to simulate effects of disorder on the photon population dynamics and the absorption spectra. We have demonstrated that the complementary employment of analytical and numerical methods permits uncovering a fairly comprehensive picture of a variety of complex behaviors in disordered multidimensional polaritonic cavity quantum electrodynamics systems.