The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

Surface hopping modeling of charge and energy transfer in active environments

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems

Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.

Unimolecular decomposition of formamide via direct chemical dynamics simulations

Classical chemical dynamics simulations show that formamide (NH₂CHO) can dissociate via multiple pathways, either by direct dissociations or via intramolecular rearrangements to different isomers followed by dissociation.

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).

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.

Excited state proton transfer dynamics of thioacetamide in S2(ππ*) state: resonance Raman spectroscopic and quantum mechanical calculations study

The photophysics and photochemistry of thioacetamide (CH3CSNH2) after excitation to the S2 electronic state were investigated by using resonance Raman spectroscopy in conjunction with the complete active space self-consistent field (CASSCF) method and density functional theory (DFT) calculations. The A-band resonance Raman spectra in acetonitrile, methanol, and water were obtained at 299.1, 282.4, 266.0, 252.7, and 245.9 nm excitation wavelengths to probe the structural dynamics of thioacetamide in the S2 state. CASSCF calculations were done to determine the transition energies and structures of the lower-lying excited states, the conical intersection points CI(S2/S1) and CI(S1/S0), and intersystem crossing points. The structural dynamics of thioacetamide in the S2 state was revealed to be along eight Franck-Condon active vibrational modes ν15, ν11, ν14, ν10, ν8, ν12, ν18, and ν19, mostly in the CC/CS/CN stretches and the CNH8,9/CCH5,6,7/CCN/CCS in-plane bends as indicated by the corresponding normal mode descriptions. The S2 → S1 decay process via the S2/S1 conical intersection point as the major channel were excluded. The thione-thiol photoisomerization reaction mechanism of thioacetamide via the S2,FC → S'1,min excited state proton transfer (ESPT) reaction channel was proposed.

Photodynamic behavior of electronic coupling in a N-methylformamide dimer

The role of the bridging water molecules has been studied during the excited state photodynamics of a N-methylformamide dimer in complex with water molecules employing the complete active space self-consistent field (CASSCF) and CAS perturbation theory (CASPT2) methods.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Nonadiabatic photodynamics of a retinal model in polar and nonpolar environment

The nonadiabatic photodynamics of the all-trans-2,4-pentadiene-iminium cation (protonated Schiff base 3, PSB3) and the all-trans-3-methyl-2,4-pentadiene-iminium cation (MePSB3) were investigated in the gas phase and in polar (aqueous) and nonpolar (n-hexane) solutions by means of surface hopping using a multireference configuration-interaction (MRCI) quantum mechanical/molecular mechanics (QM/MM) level. Spectra, lifetimes for radiationless deactivation to the ground state, and structural and electronic parameters are compared. A strong influence of the polar solvent on the location of the crossing seam, in particular in the bond length alternation (BLA) coordinate, is found. Additionally, inclusion of the polar solvent changes the orientation of the intersection cone from sloped in the gas phase to peaked, thus enhancing considerably its efficiency for deactivation of the molecular system to the ground state. These factors cause, especially for MePSB3, a substantial decrease in the lifetime of the excited state despite the steric inhibition by the solvent.