Non-adiabatic molecular dynamics with divide-and-conquer type large-scale excited-state calculations

Nonadiabatic Exciton Dynamics and Energy Gradients in the Framework of FMO-LC-TDDFTB

We introduce a novel methodology for simulating the excited-state dynamics of extensive molecular aggregates in the framework of the long-range corrected time-dependent density-functional tight-binding fragment molecular orbital method (FMO-LC-TDDFTB) combined with the mean-field Ehrenfest method. The electronic structure of the system is described in a quasi-diabatic basis composed of locally excited and charge-transfer states of all fragments. In order to carry out nonadiabatic molecular dynamics simulations, we derive and implement the excited-state gradients of the locally excited and charge-transfer states. Subsequently, the accuracy of the analytical excited-state gradients is evaluated. The applicability to the simulation of exciton transport in organic semiconductors is illustrated on a large cluster of anthracene molecules. Additionally, nonadiabatic molecular dynamics simulations of a model system of benzothieno-benzothiophene molecules highlight the method's utility in studying charge-transfer dynamics in organic materials. Our new methodology will facilitate the investigation of excitonic transfer in extensive biological systems, nanomaterials, and other complex molecular systems consisting of thousands of atoms.

Implementation of Nonadiabatic Molecular Dynamics for Intersystem Crossing Based on a Time-Dependent Density-Functional Tight-Binding Method

Intersystem crossing (ISC) and internal conversion (IC) are types of nonadiabatic transitions that play important roles in a wide range of fields, including photochemistry, photophysics, and photobiology. The nonadiabatic molecular dynamics (NA-MD) method is a powerful tool for computational simulations of dynamic phenomena involving nonadiabatic transitions. In this study, we implemented the NA-MD method, which treats ISC and IC on an equal footing, where the electronic structure is treated at the level of the time-dependent (TD) density-functional tight-binding (DFTB) method, a low-cost semiempirical analog of TD density functional theory (DFT). In particular, the spin-orbit coupling calculation algorithm was implemented in the TD-DFTB framework, and the results showed trends similar to those obtained using TD-DFT. In addition, the NA-MD method successfully reproduced ultrafast ISC of 2-nitronaphthalene.

Nonadiabatic molecular dynamics with subsystem density functional theory: application to crystalline pentacene

In this work, we report the development and assessment of the nonadiabatic molecular dynamics approach with the electronic structure calculations based on the linearly scaling subsystem density functional method. The approach is implemented in an open-source embedded Quantum Espresso/Libra software specially designed for nonadiabatic dynamics simulations in extended systems. As proof of the applicability of this method to large condensed-matter systems, we examine the dynamics of nonradiative relaxation of excess excitation energy in pentacene crystals with the simulation supercells containing more than 600 atoms. We find that increased structural disorder observed in larger supercell models induces larger nonadiabatic couplings of electronic states and accelerates the relaxation dynamics of excited states. We conduct a comparative analysis of several quantum-classical trajectory surface hopping schemes, including two new methods proposed in this work (revised decoherence-induced surface hopping and instantaneous decoherence at frustrated hops). Most of the tested schemes suggest fast energy relaxation occurring with the timescales in the 0.7-2.0 ps range, but they significantly overestimate the ground state recovery rates. Only the modified simplified decay of mixing approach yields a notably slower relaxation timescales of 8-14 ps, with a significantly inhibited ground state recovery.

Long-range corrected fragment molecular orbital density functional tight-binding method for excited states in large molecular systems

Herein, we present a new method to efficiently calculate electronically excited states in large molecular assemblies, consisting of hundreds of molecules. For this purpose, we combine the long-range corrected tight-binding density functional fragment molecular orbital method (FMO-LC-DFTB) with an excitonic Hamiltonian, which is constructed in the basis of locally excited and charge-transfer configuration state functions calculated for embedded monomers and dimers and accounts explicitly for the electronic coupling between all types of excitons. We first evaluate both the accuracy and efficiency of our fragmentation approach for molecular dimers and aggregates by comparing it with the full LC-TD-DFTB method. The comparison of the calculated spectra of an anthracene cluster shows a very good agreement between our method and the LC-TD-DFTB reference. The effective computational scaling of our method has been explored for anthracene clusters and for perylene bisimide aggregates. We demonstrate the applicability of our method by the calculation of the excited state properties of pentacene crystal models consisting of up to 319 molecules. Furthermore, the participation ratio of the monomer fragments to the excited states is analyzed by the calculation of natural transition orbital participation numbers, which are verified by the hole and particle density for a chosen pentacene cluster. The use of our FMO-LC-TDDFTB method will allow for future studies of excitonic dynamics and charge transport to be performed on complex molecular systems consisting of thousands of atoms.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

How Good Is the Vibronic Hamiltonian Repetition Approach for Long-Time Nonadiabatic Molecular Dynamics?

Multiple applied studies of slow nonadiabatic processes in nanoscale and condensed matter systems have adopted the "repetition" approximation in which long trajectories for such simulations are obtained by concatenating shorter trajectories, directly available from ab initio calculations, many times. Here, we comprehensively assess this approximation using model Hamiltonians with parameters covering a wide range of regimes. We find that state transition time scales may strongly depend on the length of the repeated data, although the convergence is not monotonic and may be slow. The repetition approach may under- or overestimate the time scales by a factor of ≤7-8, does not directly depend on the dispersion of energy gap and nonadiabatic coupling (NAC) frequencies, but may depend on the magnitude of the NACs. We suggest that the repetition-based nonadiabatic dynamics may be inaccurate in simulations with very small NACs, where intrinsic transition times are on the order of ≥100 ps.

Scalable Ehrenfest Molecular Dynamics Exploiting the Locality of Density-Functional Tight-Binding Hamiltonian

To explore the science behind excited-state dynamics in high-complexity chemical systems, a scalable nonadiabatic molecular dynamics (MD) technique is indispensable. In this study, by treating the electronic degrees of freedom at the density-functional tight-binding level, we developed and implemented a reduced scaling and multinode-parallelizable Ehrenfest MD method. To achieve this goal, we introduced a concept called patchwork approximation (PA), where the effective Hamiltonian for real-time propagation of the electronic density matrix is partitioned into a set of local parts. Numerical results for giant icosahedral fullerenes, which comprise up to 6000 atoms, suggest that the scaling of the present PA-based method is less than quadratic, which yields a significant advantage over the conventional cubic scaling method in terms of computational time. The acceleration by the parallelization on multiple nodes was also assessed. Furthermore, the electronic and structural dynamics resulting from the perturbation by the external electric field were accurately reproduced with the PA, even when the electronic excitation was spatially delocalized.

On application of deep learning to simplified quantum-classical dynamics in electronically excited states

Deep learning (DL) is applied to simulate non-adiabatic molecular dynamics of phenanthrene, using the time-dependent density functional based tight binding (TD-DFTB) approach for excited states combined with mixed quantum–classical propagation. Reference calculations rely on Tully’s fewest-switches surface hopping (FSSH) algorithm coupled to TD-DFTB, which provides electronic relaxation dynamics in fair agreement with various available experimental results. Aiming at describing the coupled electron-nuclei dynamics in large molecular systems, we then examine the combination of DL for excited-state potential energy surfaces (PESs) with a simplified trajectory surface hopping propagation based on the Belyaev–Lebedev (BL) scheme. We start to assess the accuracy of the TD-DFTB approach upon comparison of the optical spectrum with experimental and higher-level theoretical results. Using the recently developed SchNetPack (Schütt et al 2019 J. Chem. Theory Comput. 15 448–55) for DL applications, we train several models and evaluate their performance in predicting excited-state energies and forces. Then, the main focus is given to the analysis of the electronic population of low-lying excited states computed with the aforementioned methods. We determine the relaxation timescales and compare them with experimental data. Our results show that DL demonstrates its ability to describe the excited-state PESs. When coupled to the simplified BL scheme considered in this study, it provides reliable description of the electronic relaxation in phenanthrene as compared with either the experimental data or the higher-level FSSH/TD-DFTB theoretical results. Furthermore, the DL performance allows high-throughput analysis at a negligible cost.

Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites

Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.

Trajectory Surface Hopping Approach to Condensed-Phase Nonradiative Relaxation Dynamics Using Divide-and-Conquer Spin-Flip Time-Dependent Density-Functional Tight Binding

Nonradiative relaxation of excited molecules is central to many crucial issues in photochemistry. Condensed phases are typical contexts in which such problems are considered, and the nonradiative relaxation dynamics are expected to be significantly affected by interactions with the environment, for example, a solvent. We developed a nonadiabatic molecular dynamics simulation technique that can treat the nonradiative relaxation and explicitly include the environment in the calculations without a heavy computational burden. Specifically, we combined trajectory surface hopping with Tully's fewest-switches algorithm, a tight-binding approximated version of spin-flip time-dependent density-functional theory, and divide-and-conquer (DC) spatial fragmentation scheme. Numerical results showed that this method can treat systems with thousands of atoms within reasonable computational resources, and the error arising from DC fragmentation is negligibly small. Using this method, we obtained molecular insights into the solvent dependence of the photoexcited-state dynamics of trans-azobenzene, which demonstrate the importance of the environment for condensed-phase nonradiative relaxation.

Fast Nonadiabatic Molecular Dynamics via Spin-Flip Time-Dependent Density-Functional Tight-Binding Approach: Application to Nonradiative Relaxation of Tetraphenylethylene with Locked Aromatic Rings

Nonadiabatic dynamics around conical intersections between ground and excited states are crucial to understand excited-state phenomena in complex chemical systems. With this background in mind, we present an approach combining fewest-switches trajectory surface hopping and spin-flip (SF) time-dependent (TD) density-functional tight binding (DFTB), which is a simplified version of SF-TD density functional theory (DFT) with semiempirical parametrizations, for computationally efficient nonadiabatic molecular dynamics simulations. The estimated computational time of the SF-TD-DFTB approach is several orders of magnitude lower than that of SF-TD-DFT. In addition, the proposed method reproduces the time scales and quantum yields in photoisomerization reactions of azobenzene at a level comparable with conventional ab initio approaches, demonstrating reasonable accuracy. Finally, we report a practical application of the developed technique to explore the nonradiative relaxation processes of tetraphenylethylene and its derivative with torsionally locked aromatic rings and discuss the effect of locking the rings on the excited-state lifetime.