Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex.

A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.

Nonadiabatic Simulation of Exciton Dynamics in Organic Semiconductors Using Neural Network-Based Frenkel Hamiltonian and Gradients.

In this study, we present a multiscale method to simulate the propagation of Frenkel singlet excitons in organic semiconductors (OSCs). The approach uses neural network models to train a Frenkel-type Hamiltonian and its gradient, obtained by the long-range correction version of density functional tight-binding with self-consistent charges. Our models accurately predict site energies, excitonic couplings, and corresponding gradients, essential for the nonadiabatic molecular dynamics simulations. Combined with the fewest switches surface hopping algorithm, the method was applied to four representative OSCs: anthracene, pentacene, perylenediimide, and diindenoperylene. The simulated exciton diffusion constants align well with experimental and reported theoretical values and offer valuable insights into exciton dynamics in OSCs.

Modeling Charge Transport in Organic Semiconductors Using Neural Network Based Hamiltonians and Forces.

The fewest switches surface hopping method has been widely used for the simulation of charge transport in organic semiconductors. In the present study, we perform nonadiabatic molecular dynamics (NAMD) simulations of hole transport in anthracene and pentacene. The simulations employ neural network (NN) based Hamiltonians in two different nuclear relaxation schemes, which utilize either a precalculated reorganization energy or site energy gradients additionally obtained from NN models. The performance of the NN models is evaluated in reproducing hole mobilities and inverse participation ratios in terms of both quality and computational cost. The results show that charge mobilities and inverse participation ratios obtained by models, which were trained on DFTB or DFT training data, are in very good agreement with the respective QM reference method for implicit relaxation and, where available, also for explicit relaxation. Reasonable agreement with experimental hole mobilities is achieved. Utilizing our models in NAMD simulations of charge transfer amounts to a reduction of the computational cost in a range of 1 to 7 orders of magnitude compared to DFTB and DFT. This proves neural networks as promising tools for the improvement of accuracy and efficiency of charge and potentially also exciton transport simulations in complex and large molecular systems.

Charge and Exciton Transfer Simulations Using Machine-Learned Hamiltonians.

Quantum-mechanical simulations of charge and exciton transfer in molecular organic materials are a key method to increase our understanding of organic semiconductors. Our goal is to build an efficient multiscale model to predict charge-transfer mobilities and exciton diffusion constants from nonadiabatic molecular dynamics simulations and Marcus-based Monte Carlo approaches. In this work, we apply machine learning models to simulate charge and exciton propagation in organic semiconductors. We show that kernel ridge regression models can be trained to predict electronic and excitonic couplings from semiempirical density functional tight binding (DFTB) reference data with very good accuracy. In simulations, the models could reproduce hole mobilities along the anthracene crystal axes to within 8.5% of the DFTB reference and 34% of the experimental results with only 1000 training data points. Using these models decreased the cost of exciton transfer simulations by one order of magnitude.