We report the development of an efficient many-body algorithm for simulating open quantum system dynamics that utilizes a time-dependent variational principle for matrix product states to evolve large system-environment states. Capturing all system-environment correlations, we reproduce the non-perturbative, quantum-critical dynamics of the zero temperature spin-boson model, and then exploit the many-body information to visualize the complete time-frequency spectrum of the environmental excitations. Our 'environmental spectra' reveal correlated vibrational motion in polaronic modes which preserve their vibrational coherence during incoherent spin relaxation, demonstrating how environment information could yield valuable insights into complex quantum dissipative processes.

The simulation of open quantum dynamics is a critical tool for understanding how the non-classical properties of matter might be functionalised in future devices. However, unlocking the enormous potential of molecular quantum processes is highly challenging due to the very strong and non-Markovian coupling of ‘environmental’ molecular vibrations to the electronic ‘system’ degrees of freedom. Here, we present an advanced but general computational strategy that allows tensor network methods to effectively compute the non-perturbative, real-time dynamics of exponentially large vibronic wave functions of real molecules. We demonstrate how ab initio modelling, machine learning and entanglement analysis can enable simulations which provide real-time insight and direct visualisation of dissipative photophysics, and illustrate this with an example based on the ultrafast process known as singlet fission.

The process of photosynthesis, the main source of energy in the living world, converts sunlight into chemical energy. The high efficiency of this process is believed to be enabled by an interplay between the quantum nature of molecular structures in photosynthetic complexes and their interaction with the environment. Investigating these effects in biological samples is challenging due to their complex and disordered structure. Here we experimentally demonstrate a technique for studying photosynthetic models based on superconducting quantum circuits, which complements existing experimental, theoretical, and computational approaches. We demonstrate a high degree of freedom in design and experimental control of our approach based on a simplified three-site model of a pigment protein complex with realistic parameters scaled down in energy by a factor of 105. We show that the excitation transport between quantum-coherent sites disordered in energy can be enabled through the interaction with environmental noise. We also show that the efficiency of the process is maximized for structured noise resembling intramolecular phononic environments found in photosynthetic complexes.

The complex dynamics of ultrafast photoinduced reactions are governed by their evolution along vibronically coupled potential energy surfaces. It is now often possible to identify such processes, but a detailed depiction of the crucial nuclear degrees of freedom involved typically remains elusive. Here, combining excited-state time-domain Raman spectroscopy and tree-tensor network state simulations, we construct the full 108-atom molecular movie of ultrafast singlet fission in a pentacene dimer, explicitly treating 252 vibrational modes on 5 electronic states. We assign the tuning and coupling modes, quantifying their relative intensities and contributions, and demonstrate how these modes coherently synchronise to drive the reaction. Our combined experimental and theoretical approach reveals the atomic-scale singlet fission mechanism and can be generalized to other ultrafast photoinduced reactions in complex systems. This will enable mechanistic insight on a detailed structural level, with the ultimate aim to rationally design molecules to maximise the efficiency of photoinduced reactions.

We calculate the exact many-body time dynamics of polaritonic states supported by an optical cavity filled with organic molecules. Optical, vibrational and radiative processes are treated on an equal footing employing the Time-Dependent Variational Matrix Product States algorithm. We demonstrate signatures of non-Markovian vibronic dynamics and its fingerprints in the far-field photon emission spectrum at arbitrary light-matter interaction scales, ranging from the weak to the strong coupling regimes. We analyze both the single and many-molecule cases, showing the crucial role played by the collective motion of molecular nuclei and dark states in determining the polariton dynamics and the subsequent photon emission.