Efficient Surface Hopping Approach for Modeling Charge Transport in Organic Semiconductors

Mori generalized master equations offer an efficient route to predict and interpret polaron transport

Predicting how a material's microscopic structure and dynamics determine its transport properties remains a fundamental challenge. To alleviate this task's often prohibitive computational expense, we propose a Mori-based generalized quantum master equation (GQME) to predict the frequency-resolved conductivity of small-polaron forming systems described by the dispersive Holstein model. Unlike previous GQME-based approaches to transport that scale with the system size and only give access to the DC conductivity, our method requires only one calculation and yields both the DC and AC mobilities. We further show how to easily augment our GQME with numerically accessible derivatives of the current to increase computational efficiency, collectively offering computational cost reductions of up to 90%, depending on the transport regime. Finally, we leverage our exact simulations to demonstrate the limited applicability of the celebrated and widely invoked Drude-Smith model in small-polaron forming systems. We instead introduce a cumulant-based analysis of experimentally accessible frequency data to infer the microscopic Hamiltonian parameters. This approach promises to provide valuable insights into material properties and facilitate guided design by linking macroscopic terahertz measurements to the microscopic details of small polaron-forming systems.

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.

Exciton Transport in the Nonfullerene Acceptor O-IDTBR from Nonadiabatic Molecular Dynamics

Theory, computation, and experiment have given strong evidence that charge carriers in organic molecular crystals form partially delocalized quantum objects that diffuse very efficiently via a mechanism termed transient delocalization. It is currently unclear how prevalent this mechanism is for exciton transport. Here we carry out simulation of singlet Frenkel excitons (FE) in a molecular organic semiconductor that belongs to the class of nonfullerene acceptors, O-IDTBR, using the recently introduced FE surface hopping nonadiabatic molecular dynamics method. We find that FE are, on average, localized on a single molecule in the crystal due to sizable reorganization energy and moderate excitonic couplings. Yet, our simulations suggest that the diffusion mechanism is more complex than simple local hopping; in addition to hopping, we observe frequent transient delocalization events where the exciton wave function expands over 10 or more molecules for a short period of time in response to thermal excitations within the excitonic band, followed by de-excitation and contraction onto a single molecule. The transient delocalization events lead to an increase in the diffusion constant by a factor of 3-4, depending on the crystallographic direction as compared to the situation where only local hopping events are considered. Intriguingly, O-IDTBR appears to be a moderately anisotropic 3D "conductor" for excitons but a highly anisotropic 2D conductor for electrons. Taken together with previous simulation results, two trends seem to emerge for molecular organic crystals: excitons tend to be more localized and slower than charge carriers due to higher internal reorganization energy, while exciton transport tends to be more isotropic than charge transport due to the weaker distance dependence of excitonic versus electronic coupling.

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.

Theory for nonlinear conductivity switching in semiconducting organic ferroelectrics

In this work, the ferroelectric and semiconducting properties of the organic semiconducting ferroelectric benzotrithiophene tricarboxamide (BTTTA), and especially their nonlinear coupling, are theoretically investigated. BTTTA is an exponent of a small class of semiconducting organic ferroelectrics for which experiments have established a surprising polarization direction dependence of the bulk conductivity at finite fields. First, molecular dynamics (MD) simulations are used to investigate the occurrence and, under the influence of an external electric field, the inversion of the macroscopic electric dipole that forms along the axis of supramolecular columns of BTTTA. The MD results are consistent with the experimentally observed ferroelectric behavior of the material. Building on the MD results, a QM/MM scheme is used to investigate the charge carrier mobility in the quasi-1D BTTTA stacks in the linear and non-linear regimes. Indeed, at finite electric fields, a clear resistance switching effect was observed in the form of a hole mobility that is a factor ∼2 larger for antiparallel orientations of the polarization and field than for a parallel orientation. This phenomenon can be understood as a microscopic ratchet that is based on the non-equilibrium interaction between the (oriented) dipoles and the (direction of the) charge transport.

Exciton Dissociation in a Model Organic Interface: Excitonic State-Based Surface Hopping versus Multiconfigurational Time-Dependent Hartree

Quantum dynamical simulations are essential for a molecular-level understanding of light-induced processes in optoelectronic materials, but they tend to be computationally demanding. We introduce an efficient mixed quantum-classical nonadiabatic molecular dynamics method termed eXcitonic state-based Surface Hopping (X-SH), which propagates the electronic Schrödinger equation in the space of local excitonic and charge-transfer electronic states, coupled to the thermal motion of the nuclear degrees of freedom. The method is applied to exciton decay in a 1D model of a fullerene-oligothiophene junction, and the results are compared to the ones from a fully quantum dynamical treatment at the level of the Multilayer Multiconfigurational Time-Dependent Hartree (ML-MCTDH) approach. Both methods predict that charge-separated states are formed on the 10-100 fs time scale via multiple "hot-exciton dissociation" pathways. The results demonstrate that X-SH is a promising tool advancing the simulation of photoexcited processes from the molecular to the true nanomaterials scale.