An exciton–phonon (ex–ph) model based on our recently developed block interaction product basis framework is introduced to simulate the distal charge separation (CS) process in aggregated perylenediimide (PDI) trimer incorporating the quantum dynamic method, i.e., the time-dependent density matrix renormalization group. The electronic Hamiltonian in the ex–ph model is represented by nine constructed diabatic states, which include three local excited (LE) states and six charge transfer (CT) states from both the neighboring and distal chromophores. These diabatic states are automatically generated from the direct products of the leading localized neutral or ionic states of each chromophore’s reduced density matrix, which are obtained from ab initio quantum chemical calculation of the subsystem consisting of the targeted chromophore and its nearest neighbors, thus considering the interaction of the adjacent environment. In order to quantum-dynamically simulate the distal CS process with massive coupled vibrational modes in molecular aggregates, we used our recently proposed hierarchical mapping approach to renormalize these modes and truncate those vibrational modes that are not effectively coupled with electronic states accordingly. The simulation result demonstrates that the formation of the distal CS process undergoes an intermediate state of adjacent CT, i.e., starts from the LE states, passes through an adjacent CT state to generate the intermediates (∼200 fs), and then formalizes the targeted distal CS via further charge transference (∼1 ps). This finding agrees well with the results observed in the experiment, indicating that our scheme is capable of quantitatively investigating the CS process in a realistic aggregated PDI trimer and can also be potentially applied to exploring CS and other photoinduced processes in larger systems.
