Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics

In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.

Accurate Computation of Nonadiabatic Coupling with Projector Augmented-Wave Pseudopotentials

Synergy of nonadiabatic molecular dynamics with real-time time-dependent density functional theory has led to significant progress in modeling excited-state dynamics in nanoscale and condensed matter systems over the past decade. Nonadiabatic coupling (NAC) is the central quantity in such simulations, and its accurate and efficient evaluation is an enduring challenge in time-dependent Kohn-Sham theory, particularly in conjunction with planewave basis sets and projector augmented-wave (PAW) pseudopotentials because of the complexity of the PAW "all-electron" wave function. We report a method for rigorous evaluation of the NAC with PAW wave functions and demonstrate an efficient approximation to the rigorous NAC that gives comparable accuracy. As a validation, we intensely examine the NAC matrix elements calculated using both pseudo- and all-electron wave functions under the PAW formalism in six representative systems. The approximate NAC obtained with pseudowave functions is close to the exact all-electron NAC, with the largest deviations observed when subshell d-electrons are involved in the transitions. The developed approach provides a rigorous and convenient methodology for the numerical computation of NAC in the Kohn-Sham theory framework.

Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications

This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.

Dependence of electron transfer dynamics on the number of graphene layers in π-stacked 2D materials: insights from ab initio nonadiabatic molecular dynamics

Recent time-resolved transient absorption studies demonstrated that the rate of photoinduced interfacial charge transfer (CT) from Zn-phthalocyanine (ZnPc) to single-layer graphene (SLG) is faster than to double-layer graphene (DLG), in contrast to the expectation from Fermi's golden rule. We present the first time-domain non-adiabatic molecular dynamics (NA-MD) study of the electron injection process from photoexcited ZnPc molecules into SLG and DLG substrates. Our calculations suggest that CT occurs faster in the ZnPc/SLG system than in the ZnPc/DLG system, with 580 fs and 810 fs being the fastest components of the observed CT timescales, respectively. The computed timescales are in close agreement with those reported in the experiment. The computed CT timescales are determined largely by the magnitudes of the non-adiabatic couplings (NAC), which we find to be 4 meV and 2 meV, for the ZnPc/SLG and ZnPc/DLG systems, respectively. The transitions are driven mainly by the ZnPc out-of-plane bending mode at 1100 cm-1 and an overtone of fundamental modes in graphene at 2450 cm-1. We find that dephasing occurs on the timescale of 20 fs and is similar in both systems, so decoherence does not notably change the qualitative trends in the CT timescales. We highlight the importance of proper energy level alignment for capturing the qualitative trends in the CT dynamics observed in experiment. In addition, we illustrate several methodological points that are important for accurately modeling nonadiabatic dynamics in the ZnPc/FLG systems, such as the choice of surface hopping methodology, the use of phase corrections, NAC scaling, and the inclusion of Hubbard terms in the density functional and molecular dynamics calculations.

Ultrafast charge transfer coupled with lattice phonons in two-dimensional covalent organic frameworks

Covalent organic frameworks (COFs) have emerged as a promising light-harvesting module for artificial photosynthesis and photovoltaics. For efficient generation of free charge carriers, the donor–acceptor (D-A) conjugation has been adopted for two-dimensional (2D) COFs recently. In the 2D D-A COFs, photoexcitation would generate a polaron pair, which is a precursor to free charge carriers and has lower binding energy than an exciton. Although the character of the primary excitation species is a key factor in determining optoelectronic properties of a material, excited-state dynamics leading to the creation of a polaron pair have not been investigated yet. Here, we investigate the dynamics of photogenerated charge carriers in 2D D-A COFs by combining femtosecond optical spectroscopy and non-adiabatic molecular dynamics simulation. From this investigation, we elucidate that the polaron pair is formed through ultrafast intra-layer hole transfer coupled with coherent vibrations of the 2D lattice, suggesting a mechanism of phonon-assisted charge transfer.

The donor–acceptor (D-A) conjugation has been adopted for two-dimensional (2D) covalent organic frameworks (COFs) for efficient generation of free charge carriers. Here, the authors investigate the dynamics of photogenerated charge carriers in 2D D-A COFs by combining femtosecond optical spectroscopy and non-adiabatic molecular dynamics simulation.

A Simple Phase Correction Makes a Big Difference in Nonadiabatic Molecular Dynamics

The outcomes of nonadiabatic molecular dynamics (NA-MD) calculations are modulated by the parameters entering the time-dependent Schrödinger equation (TD-SE). The adiabatic states are commonly used as the basis in which the TD-SE is integrated. However, the phase inconsistencies of such states along the nuclear trajectories obtained in NA-MD simulations may render the wave function and other relevant properties ill-behaving, adversely affecting the dynamics. This work illustrates the consequence of adiabatic state phase inconsistencies in nonadiabatic Ehrenfest dynamics. A simple phase-correction approach is proposed and is demonstrated to alter the dynamics to make it consistent with the reference calculations done in the phase-consistent diabatic representation.

Stochastic and Quasi-Stochastic Hamiltonians for Long-Time Nonadiabatic Molecular Dynamics

In the condensed-matter environments, the vibronic Hamiltonian that describes nonadiabatic dynamics often appears as an erratic entity, and one may assume it can be generated stochastically. This property is utilized to formulate novel stochastic and quasi-stochastic vibronic Hamiltonian methodologies, which open a new route to long-time excited state dynamics in atomistic solid-state systems at negligible computational cost. Using a model mimicking a typical solid-state material in noisy environment, general conclusions regarding the simulation of nonadiabatic dynamics are obtained: (1) including bath is critical to complete excited state relaxation; (2) a totally stochastic modulation of energies and couplings has a net effect of no bath and inhibits relaxation; (3) including a single or several dominant electron-phonon modes may be insufficient to complete the excited state relaxation; (4) only the multiple modes, even those that have negligible weights, can represent both the deterministic modulation of system's Hamiltonian and stochastic effects of bath.