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.

We develop a systematic variational coherent-state expansion for the many-body ground state of the spin-boson model, in which a quantum two-level system is coupled to a continuum of harmonic oscillators. Energetic constraints at the heart of this technique are rationalized in terms of polarons (displacements of the bath states in agreement with classical expectations) and antipolarons (counterdisplacements due to quantum tunneling effects). We present a comprehensive study of the ground-state two-level system population and coherence as a function of tunneling amplitude, dissipation strength, and bias (akin to asymmetry of the double-well potential defining the two-state system). The entanglement among the different environmental modes is investigated by looking at spectroscopic signatures of the bipartite entanglement entropy between a given environmental mode and all the other modes. We observe a drastic change in behavior of this entropy for increasing dissipation, indicative of the entangled nature of the environmental states. In addition, the entropy spreads over a large energy range at strong dissipation, a testimony to the wide entanglement window characterizing the underlying Kondo state. Finally, comparisons to accurate numerical renormalization-group calculations and to the exact Bethe ansatz solution of the model demonstrate the rapid convergence of our variationally optimized multipolaron expansion, suggesting that it should also be a useful tool for dissipative models of greater complexity, as relevant for numerous systems of interest in quantum physics and chemistry.

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 analyze the role of coherent, non-perturbative system-bath interactions in a photosynthetic heat engine. Using the reaction-coordinate formalism to describe the vibrational phonon-environment in the engine, we analyze the efficiency around an optimal parameter regime predicted in earlier studies. We show that, in the limit of high-temperature photon irradiation, the phonon-assisted population transfer between bright and dark states is suppressed due to dephasing from the photon environment, even in the Markov limit where we expect the influence of each bath to have an independent and additive effect on the dynamics. Manipulating the phonon bath properties via its spectral density enables us to identify both optimal low- and high-frequency regimes where the suppression can be removed. This suppression of transfer and its removal suggests that it is important to consider carefully the non-perturbative and cooperative effects of system-bath environments in designing artificial photosynthetic systems and also that manipulating inter-environmental interactions could provide a new multidimensional “lever” by which photocells and other types of quantum devices can be optimized.

Recent tensor network techniques for simulating system-environment wave functions have provided profound insights into non-Markovian dissipation and decoherence in open quantum systems. Here, we propose a dynamically adaptive one-site time-dependent-variational-principle (A1TDVP) method for matrix product states in which local bond dimensions grow to capture developing system-bath entanglement. This avoids the need for multiple convergence runs with respect to bond dimensions and the unfavorable local Hilbert space scaling of two-site methods. A1TDVP is thus ideally suited for open quantum dynamics in finite-temperature bosonic environments, as the initial states typically have low bond dimension but require very large local physical dimensions. We demonstrate this with simulations of nonequilibrium heat flows through a qubit spin, finding a $30\times$ and $10\times$ speed-up over 2TDVP and 1TDVP, respectively.

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S₁ and the dark S₂ states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.