Implementation of surface hopping molecular dynamics using semiempirical methods

An Efficient, Augmented Surface Hopping Algorithm That Includes Decoherence for Use in Large-Scale Simulations

We propose and implement a highly efficient augmented surface hopping algorithm that (i) can be used for large simulations (with many nuclei and many electronic states) and (ii) includes the effects of decoherence without parametrization. Our protocol is based on three key modifications of the surface hopping methodology: (a) a novel separation of classical and quantum degrees of freedom that treats avoided and trivial crossings efficiently, (b) a multidimensional approximation of the time derivative matrix that avoids explicit construction of the derivative coupling at most time steps, and (c) an efficient approximation for the augmented fewest-switches surface hopping decoherence rate. We will show that this protocol can be several orders of magnitude more efficient than the traditional protocol for large multidimensional problems. Furthermore, the marginal cost for including decoherence effects is now negligible.

An On-the-Fly Surface-Hopping Program JADE for Nonadiabatic Molecular Dynamics of Polyatomic Systems: Implementation and Applications

Nonadiabatic dynamics simulations have rapidly become an indispensable tool for understanding ultrafast photochemical processes in complex systems. Here, we present our recently developed on-the-fly nonadiabatic dynamics package, JADE, which allows researchers to perform nonadiabatic excited-state dynamics simulations of polyatomic systems at an all-atomic level. The nonadiabatic dynamics is based on Tully's surface-hopping approach. Currently, several electronic structure methods (CIS, TDHF, TDDFT(RPA/TDA), and ADC(2)) are supported, especially TDDFT, aiming at performing nonadiabatic dynamics on medium- to large-sized molecules. The JADE package has been interfaced with several quantum chemistry codes, including Turbomole, Gaussian, and Gamess (US). To consider environmental effects, the Langevin dynamics was introduced as an easy-to-use scheme into the standard surface-hopping dynamics. The JADE package is mainly written in Fortran for greater numerical performance and Python for flexible interface construction, with the intent of providing open-source, easy-to-use, well-modularized, and intuitive software in the field of simulations of photochemical and photophysical processes. To illustrate the possible applications of the JADE package, we present a few applications of excited-state dynamics for various polyatomic systems, such as the methaniminium cation, fullerene (C20), p-dimethylaminobenzonitrile (DMABN) and its primary amino derivative aminobenzonitrile (ABN), and 10-hydroxybenzo[h]quinoline (10-HBQ).

Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Theory, Implementation, and Parameters

Semiempirical orthogonalization-corrected methods (OM1, OM2, and OM3) go beyond the standard MNDO model by explicitly including additional interactions into the Fock matrix in an approximate manner (Pauli repulsion, penetration effects, and core–valence interactions), which yields systematic improvements both for ground-state and excited-state properties. In this Article, we describe the underlying theoretical formalism of the OMx methods and their implementation in full detail, and we report all relevant OMx parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine. For a standard set of mostly organic molecules commonly used in semiempirical method development, the OMx results are found to be superior to those from standard MNDO-type methods. Parametrized Grimme-type dispersion corrections can be added to OM2 and OM3 energies to provide a realistic treatment of noncovalent interaction energies, as demonstrated for the complexes in the S22 and S66×8 test sets.

Excited-State Proton-Transfer-Induced Trapping Enhances the Fluorescence Emission of a Locked GFP Chromophore

The chemical locking of the central single bond in core chromophores of green fluorescent proteins (GFPs) influences their excited-state behavior in a distinct manner. Experimentally, it significantly enhances the fluorescence quantum yield of GFP chromophores with an ortho-hydroxyl group, while it has almost no effect on the photophysics of GFP chromophores with a para-hydroxyl group. To unravel the underlying physical reasons for this different behavior, we report static electronic structure calculations and nonadiabatic dynamics simulations on excited-state intramolecular proton transfer, cis–trans isomerization, and excited-state deactivation in a locked ortho-substituted GFP model chromophore (o-LHBI). On the basis of our previous and present results, we find that the S₁ keto species is responsible for the fluorescence emission of the unlocked o-HBI and the locked o-LHBI species. Chemical locking does not change the parts of the S₁ and S₀ potential energy surfaces relevant to enol–keto tautomerization; hence, in both chromophores, there is an ultrafast excited-state intramolecular proton transfer that takes only 35 fs on average. However, the locking effectively hinders the S₁ keto species from approaching the keto S₁/S₀ conical intersections so that most of trajectories are trapped in the S₁ keto region for the entire 2 ps simulation time. Therefore, the fluorescence quantum yield of o-LHBI is enhanced compared with that of unlocked o-HBI, in which the S₁ excited-state decay is efficient and ultrafast. In the case of the para-substituted GFP model chromophores p-HBI and p-LHBI, chemical locking hardly affects their efficient excited-state deactivation via cis–trans isomerization; thus, the fluorescence quantum yields in these chromophores remain very low. The insights gained from the present work may help to guide the design of new GFP chromophores with improved fluorescence emission and brightness.

Photophysics of Auramine-O: electronic structure calculations and nonadiabatic dynamics simulations

Sequential vs. concerted S₁ relaxation pathways.

Diabatic Population Matrix Formalism for Performing Molecular Mechanics Style Simulations with Multiple Electronic States

An accurate description of nonbonded interactions is important in investigating dynamics of molecular systems. In many situations, fixed point charge models are successfully applied to explaining various chemical phenomena. However, these models with conventional formulations will not be appropriate in elucidating the detailed dynamics during nonadiabatic events. This is mainly because the chemical properties of any molecule, especially its electronic populations, significantly change with respect to molecular distortions in the vicinity of the surface crossing. To overcome this issue in molecular simulations yet within the framework of the fixed point charge model, we define a diabatic electronic population matrix and substitute it for the conventional adiabatic partial charges. We show that this matrix can be readily utilized toward attaining more reliable descriptions of Coulombic interactions, in combination with the interpolation formalism for obtaining the intramolecular interaction potential. We demonstrate how the mixed formalism with the diabatic charges and the interpolation can be applied to molecular simulations by conducting adiabatic and nonadiabatic molecular dynamics trajectory calculations of the green fluorescent protein chromophore anion in aqueous environment.

Simulation of femtosecond two-dimensional electronic spectra of conical intersections

We have simulated femtosecond two-dimensional (2D) electronic spectra for an excited-state conical intersection using the wave-function version of the equation-of-motion phase-matching approach. We show that 2D spectra at fixed values of the waiting time provide information on the structure of the vibronic eigenstates of the conical intersection, while the evolution of the spectra with the waiting time reveals predominantly ground-state wave-packet dynamics. The results show that 2D spectra of conical intersection systems differ significantly from those obtained for chromophores with well separated excited-state potential-energy surfaces. The spectral signatures which can be attributed to conical intersections are discussed.

Decoherence Allows Model Reduction in Nonadiabatic Dynamics Simulations

A nonadiabatic (NA) molecular dynamics (MD) simulation requires calculation of NA coupling matrix elements, the number of which scales as a square of the number of basis states. The basis size can be huge in studies of nanoscale materials, and calculation of the NA couplings can present a significant bottleneck. A quantum-classical approximation, NAMD overestimates coherence in the quantum, electronic subsystem, requiring decoherence correction. Generally, decoherence times decrease with increasing energy separation between pairs of states forming coherent superpositions. Since rapid decoherence stops quantum dynamics, one expects that decoherence-corrected NAMD can eliminate the need for calculation of NA couplings between energetically distant states, notably reducing the computational cost. Considering several types of dynamics in a semiconductor quantum dot, we demonstrate that indeed, decoherence allows one to reduce the number of needed NA coupling matrix elements. If the energy levels are spaced closer than 0.1 eV, one obtains good results while including only three nearest-neighbor couplings, and in some cases even with just the first nearest-neighbor coupling scheme. If the energy levels are spaced by about 0.4 eV, the nearest-neighbor model fails, while three or more nearest-neighbor schemes also provide good results. In comparison, the results of NAMD simulation without decoherence vary continuously with changes in the number of NA couplings. Thus, decoherence effects induced by coupling to a quantum-mechanical environment not only provide the physical mechanism for NAMD trajectory branding and improve the accuracy of NAMD simulations, but also afford significant computational savings.

Excited States and Photochemistry of Chromophores in the Photoactive Proteins Explored by the Combined Quantum Mechanical and Molecular Mechanical Calculations

A photoactive protein usually contains a unique chromophore that is responsible for the initial photoresponse and functions of the photoactive protein are determined by the interaction between the chromophore and its protein surroundings. The combined quantum mechanical and molecular mechanical (QM/MM) approach is demonstrated to be a very useful tool for exploring structures and functions of a photoactive protein with the chromophore and its protein surroundings treated by the QM and MM methods, respectively. In this review, we summarize the basic formulas of the QM/MM approach and emphasize its applications to excited states and photoreactions of chromophores in rhodopsin protein, photoactive yellow protein, and green fluorescent protein.

Generalized trajectory surface-hopping method for internal conversion and intersystem crossing

Trajectory-based fewest-switches surface-hopping (FSSH) dynamics simulations have become a popular and reliable theoretical tool to simulate nonadiabatic photophysical and photochemical processes. Most available FSSH methods model internal conversion. We present a generalized trajectory surface-hopping (GTSH) method for simulating both internal conversion and intersystem crossing processes on an equal footing. We consider hops between adiabatic eigenstates of the non-relativistic electronic Hamiltonian (pure spin states), which is appropriate for sufficiently small spin-orbit coupling. This choice allows us to make maximum use of existing electronic structure programs and to minimize the changes to available implementations of the traditional FSSH method. The GTSH method is formulated within the quantum mechanics (QM)/molecular mechanics framework, but can of course also be applied at the pure QM level. The algorithm implemented in the GTSH code is specified step by step. As an initial GTSH application, we report simulations of the nonadiabatic processes in the lowest four electronic states (S0, S1, T1, and T2) of acrolein both in vacuo and in acetonitrile solution, in which the acrolein molecule is treated at the ab initio complete-active-space self-consistent-field level. These dynamics simulations provide detailed mechanistic insight by identifying and characterizing two nonadiabatic routes to the lowest triplet state, namely, direct S1 → T1 hopping as major pathway and sequential S1 → T2 → T1 hopping as minor pathway, with the T2 state acting as a relay state. They illustrate the potential of the GTSH approach to explore photoinduced processes in complex systems, in which intersystem crossing plays an important role.

Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.

Local Control Theory in Trajectory Surface Hopping Dynamics Applied to the Excited-State Proton Transfer of 4-Hydroxyacridine

The application of local control theory combined with nonadiabatic ab initio molecular dynamics to study the photoinduced intramolecular proton transfer reaction in 4-hydroxyacridine was investigated. All calculations were performed within the framework of linear-response time-dependent density functional theory. The computed pulses revealed important information about the underlying excited-state nuclear dynamics highlighting the involvement of collective vibrational modes that would normally be neglected in a study performed on model systems constrained to a subset of the full configuration space. This study emphasizes the strengths of local control theory for the design of pulses that can trigger chemical reactions associated with the population of a given molecular excited state. In addition, analysis of the generated pulses can help to shed new light on the photophysics and photochemistry of complex molecular systems.

Constructing polyatomic potential energy surfaces by interpolating diabatic Hamiltonian matrices with demonstration on green fluorescent protein chromophore

Simulating molecular dynamics directly on quantum chemically obtained potential energy surfaces is generally time consuming. The cost becomes overwhelming especially when excited state dynamics is aimed with multiple electronic states. The interpolated potential has been suggested as a remedy for the cost issue in various simulation settings ranging from fast gas phase reactions of small molecules to relatively slow condensed phase dynamics with complex surrounding. Here, we present a scheme for interpolating multiple electronic surfaces of a relatively large molecule, with an intention of applying it to studying nonadiabatic behaviors. The scheme starts with adiabatic potential information and its diabatic transformation, both of which can be readily obtained, in principle, with quantum chemical calculations. The adiabatic energies and their derivatives on each interpolation center are combined with the derivative coupling vectors to generate the corresponding diabatic Hamiltonian and its derivatives, and they are subsequently adopted in producing a globally defined diabatic Hamiltonian function. As a demonstration, we employ the scheme to build an interpolated Hamiltonian of a relatively large chromophore, para-hydroxybenzylidene imidazolinone, in reference to its all-atom analytical surface model. We show that the interpolation is indeed reliable enough to reproduce important features of the reference surface model, such as its adiabatic energies and derivative couplings. In addition, nonadiabatic surface hopping simulations with interpolation yield population transfer dynamics that is well in accord with the result generated with the reference analytic surface. With these, we conclude by suggesting that the interpolation of diabatic Hamiltonians will be applicable for studying nonadiabatic behaviors of sizeable molecules.

Modeling ultrafast exciton deactivation in oligothiophenes via nonadiabatic dynamics

Nonadiabatic excited-state dynamics reveal the exciton relaxation processes in oligothiophenes. Ultrafast deactivation and exciton localization are predicted to occur within 200 fs, involving bond stretching, ring puckering, and torsional oscillations.

Tracking the primary photoconversion events in rhodopsins by ultrafast optical spectroscopy

We review the most recent experimental and computational efforts aimed at exposing the very early phases of the ultrafast isomerization in visual Rhodopsins and we discuss future advanced experiments and calculations.

Photoswitching of salicylidene methylamine: a theoretical photodynamics study

Photoswitching of simple photochromic molecules attracts substantial attention because of its possible role in future photon-driven molecular electronics. Here we model the full photoswitching cycle of a minimal photochromic Schiff base–salicylidene methylamine (SMA). We perform semiempirical nonadiabatic on-the-fly photodynamics simulations at the OM2/MRCI level and thoroughly analyze the structural time evolution and switching efficiency of the system. We also identify and examine in detail the crucial steps in the SMA photochemistry ruled by excited-state intramolecular proton transfer. The results place the investigated model aromatic Schiff base among the promising candidates for novel photoswitching molecular materials. Our study also shows the potential of the semiempirical multireference photodynamics simulations as a tool for early stage molecular photodevice design.

A multireference configuration interaction study of the photodynamics of nitroethylene

Extended multireference configuration interaction with singles and doubles (MR-CISD) calculations of nitroethylene (H₂C=CHNO₂) were carried out to investigate the photodynamical deactivation paths to the ground state. The ground (S₀) and the first five valence excited electronic states (S₁–S₅) were investigated. In the first step, vertical excitations and potential energy curves for CH₂ and NO₂ torsions and CH₂ out-of-plane bending starting from the ground state geometry were computed. Afterward, five conical intersections, one between each pair of adjacent states, were located. The vertical calculations mostly confirm the previous assignment of experimental spectrum and theoretical results using lower-level calculations. The conical intersections have as main features the torsion of the CH₂ moiety, different distortions of the NO₂ group and CC, CN, and NO bond stretchings. In these conical intersections, the NO₂ group plays an important role, also seen in excited state investigations of other nitro molecules. Based on the conical intersections found, a photochemical nonradiative deactivation process after a π–π* excitation to the bright S₅ state is proposed. In particular, the possibility of NO₂ release in the ground state, an important property in nitro explosives, was found to be possible.

Evaluation of the Time-Derivative Coupling for Accurate Electronic State Transition Probabilities from Numerical Simulations

Spikes in the time-derivative coupling (TDC) near surface crossings make the accurate integration of the time-dependent Schrödinger equation in nonadiabatic molecular dynamics simulations a challenge. To address this issue, we present an approximation to the TDC based on a norm-preserving interpolation (NPI) of the adiabatic electronic wave functions within each time step. We apply NPI and two other schemes for computing the TDC in numerical simulations of the Landau-Zener model, comparing the simulated transfer probabilities to the exact solution. Though NPI does not require the analytical calculation of nonadiabatic coupling matrix elements, it consistently yields unsigned population transfer probability errors of ∼0.001, whereas analytical calculation of the TDC yields errors of 0.0-1.0 depending on the time step, the offset of the maximum in the TDC from the beginning of the time step, and the coupling strength. The approximation of Hammes-Schiffer and Tully yields errors intermediate between NPI and the analytical scheme.

A Simple Solution to the Trivial Crossing Problem in Surface Hopping

Surface hopping studies on supramolecular and nanoscale systems suffer severely from the trivial crossing problem, arising due to high density of adiabatic potential energy surfaces. We present a straightforward solution to the problem by introducing a self-consistency test to the well-known fewest switches surface hopping (FSSH) procedure. If the test is failed, the hopping probabilities are corrected with a simple procedure. The novel self-consistent fewest switches surface hopping (SC-FSSH) approach is applied to the Holstein Hamiltonian to study the time-dependence of the electron population. Already in the five-state system, SC-FSSH allows us to reduce the simulation time 10(4)-fold to achieve the FSSH accuracy. The reliable performance and simple formulation of SC-FSSH greatly expands the applicability range of the surface hopping method.

Photochemistry of aldehyde clusters: cross-molecular versus unimolecular reaction dynamics

Semi-empirical on-the-fly dynamics simulations reveal the importance of cross-molecular reactions in the photochemistry of aldehyde clusters.

Ab initio and semi-empirical Molecular Dynamics simulations of chemical reactions in isolated molecules and in clusters

Recent progress in “on-the-fly” trajectory simulations of molecular reactions, using different electronic structure methods is discussed, with analysis of the insights that such calculations can provide and of the strengths and limitations of the algorithms available.