DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex
Sayan Maity, Beatrix M. Bold, Jigneshkumar Dahyabhai Prajapati, M.V. Sokolov, Tomáš Kubař, Marcus Elstner, Ulrich Kleinekathöfer
Abstract
Because of the size of light-harvesting complexes and the involvement of electronic degrees of freedom, computationally these systems need to be treated with a combined quantum–classical description. To this end, Born−Oppenheimer molecular dynamics simulations have been employed in a quantum mechanics/molecular mechanics (QM/MM) fashion for the ground state followed by excitation energy calculations again in a QM/MM scheme for the Fenna−Matthews−Olson (FMO) complex. The self-consistent-charge density functional tight-binding (DFTB) method electrostatically coupled to a classical description of the environment was applied to perform the ground-state dynamics. Subsequently, long-range-corrected time-dependent DFTB calculations were performed to determine the excitation energy fluctuations of the individual bacteriochlorophyll a molecules. The spectral densities obtained using this approach show an excellent agreement with experimental findings. In addition, the fluctuating site energies and couplings were used to estimate the exciton transfer dynamics.