“…Nonadiabatic effects due to the interactions of multiple electronic excited states are ubiquitous in many molecular systems, ranging from biological light-harvesting complexes to quantum dots and organic optoelectronic materials. − They arise when nuclear motion mixes adiabatic electronic states, the true eigenstates of the electronic Hamiltonian, representing a breakdown of the Born–Oppenheimer approximation. , Signatures of such nonadiabatic effects , can be commonly observed in linear and nonlinear optical spectra of molecules and molecular assemblies, with prominent examples such as dark (dipole-forbidden) transitions contributing to spectral line shapes through intensity borrowing effects, − anomalous emission properties like those observed for molecules undergoing twisted intramolecular charge transfer (TICT), − or molecules violating Kasha’s rule. ,− First-principles modeling of these effects is highly challenging, as it must account for explicit couplings between the electronic and nuclear degrees of freedom. , Additionally, most systems of interest are embedded in condensed-phase environments such as solvents or proteins, and understanding how environmental degrees of freedom influence energy relaxation processes in multiple coupled electronic states poses a major theoretical challenge. − In this work, we introduce a first-principles approach, based on combining molecular dynamics (MD) sampling of vibronic couplings with tensor network methods, ,− that is capable of capturing environmental interactions on several coupled excited states in absorption, fluorescence, and time-resolved spectroscopy experiments. We showcase the strengths of the approach by uncovering the origin of the recently reported dual fluorescence in the proflavine molecule …”