All-Atom Photoinduced Charge Transfer Dynamics in Condensed Phase via Multistate Nonlinear-Response Instantaneous Marcus Theory

Photoinduced charge transfer (CT) in the condensed phase is an essential component in solar energy conversion, but it is challenging to simulate such a process on the all-atom level. The traditional Marcus theory has been utilized for obtaining CT rate constants between pairs of electronic states but cannot account for the nonequilibrium effects due to the initial nuclear preparation. The recently proposed instantaneous Marcus theory (IMT) and its nonlinear-response formulation allow for incorporating the nonequilibrium nuclear relaxation to electronic transition between two states after the photoexcitation from the equilibrium ground state and provide the time-dependent rate coefficient. In this work, we extend the nonlinear-response IMT method for treating photoinduced CT among general multiple electronic states and demonstrate it in the organic photovoltaic carotenoid-porphyrin-fullerene triad dissolved in explicit tetrahydrofuran solvent. All-atom molecular dynamics simulations were employed to obtain the time correlation functions of energy gaps, which were used to generate the IMT-required time-dependent averages and variances of the relevant energy gaps. Our calculations show that the multistate IMT could capture the significant nonequilibrium effects due to the initial nuclear state preparation, and this is corroborated by the substantial differences between the population dynamics predicted by the multistate IMT and the Marcus theory, where the Marcus theory underestimates the population transfer. The population dynamics by multistate IMT is also shown to have a better agreement with the all-atom nonadiabatic mapping dynamics than the Marcus theory does. Because the multistate nonlinear-response IMT is straightforward and cost-effective in implementation and accounts for the nonequilibrium nuclear effects, we believe this method offers a practical strategy for studying charge transfer dynamics in complex condensed-phase systems.

Modeling the Effect of Disorder in the Two-Dimensional Electronic Spectroscopy of Poly-3-hexyltiophene in an Organic Photovoltaic Blend: A Combined Quantum/Classical Approach

We introduce a first-principles model of the 12-mer poly-3-hexyltiophene (P3HT) polymer system in the realistic description of an organic photovoltaic blend environment. We combine Molecular Dynamics (MD) simulations of a thin-film blend of P3HT and phenyl-C61-butyric acid methyl ester (PCBM) to model the interactions with a fluctuating environment with Time-Dependent Density Functional Theory (TDDFT) calculations to parametrize the effect of the torsional flexibility in the polymer and construct an exciton-type Hamiltonian that describes the photoexcitation of the polymer. This allows us to reveal the presence of different flexibility patterns governed by the torsional angles along the polymer chain which, in the interacting fluctuating environment, control the broadening of the spectral observables. We identify the origin of the homogeneous and inhomogeneous line shape of the simulated optical signals. This is paramount to decipher the spectroscopic nature of the ultrafast electron-transfer process occurring in organic photovoltaic (OPV) materials.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.