Probing ultrafast excitation energy transfer of the chlorosome with exciton-phonon variational dynamics

The chlorosome antenna complex is a fascinating structure which due to its immense scale, accurate simulation of excitation energy transfer (EET) dynamics supposes a genuine computational challenge. Resonant vibronic modes have been recently identified in 2D spectra of the chlorosome which motivates our present endeavour of modelling electronic and vibrational degrees of freedom on an equal footing. Following the Dirac-Frenkel time-dependent variational principle, we exploit a general theory of polaron dynamics in two-dimensional lattices based on the Holstein molecular crystal model and investigate a single rod model of pigment aggregates. Unlike reduced formalisms, explicit integration of the degrees of freedom of both the system and the bath requires extensive computational resources. We exploit the architecture of graphic processor units (GPUs) by implementing our simulations on this platform. The simulation of dynamic properties of hundreds or even thousands of pigments is thus achievable in just a few hours. The potential investigation and design of natural or engineered two-dimensional pigment networks can thus be accommodated. Due to the lack of consensus regarding the precise arrangement of chromophores in the chlorosome, helicity and dimerization are investigated independently, extracting their contributions to both optical and EET properties. The presence of dimerization is found to slow down the delocalization process. Exciton delocalization is completed in 100 fs in a single rod aggregate whose dimensions (20 nm) fairly exceed the estimated extent of a coherent domain. Ultrafast energy relaxation in the exciton manifold occurs in 50 fs and the duration of super-diffusive transport is found to last for about 80 fs.