The intrinsic reaction coordinate (IRC) approach has been used extensively in quantum chemical analysis and prediction of the mechanism of chemical reactions. The IRC gives a unique connection from a given transition structure to local minima of the reactant and product sides. This allows for easy understanding of complicated multistep mechanisms as a set of simple elementary reaction steps. In this article, three topics concerning the IRC approach are discussed. In the first topic, the first ab initio study of the IRC and a recent development of an IRC calculation algorithm for enzyme reactions are introduced. In the second topic, cases are presented in which dynamical trajectories bifurcate and corresponding IRC connections can be inaccurate. In the third topic, a recent development of an automated reaction path search method and its application to systematic construction of IRC networks are described. Finally, combining these three topics, future perspectives are discussed.

This article reports implementation and performance of the artificial force induced reaction (AFIR) method in the upcoming 2017 version of GRRM program (GRRM17). The AFIR method, which is one of automated reaction path search methods, induces geometrical deformations in a system by pushing or pulling fragments defined in the system by an artificial force. In GRRM17, three different algorithms, that is, multicomponent algorithm (MC‐AFIR), single‐component algorithm (SC‐AFIR), and double‐sphere algorithm (DS‐AFIR), are available, where the MC‐AFIR was the only algorithm which has been available in the previous 2014 version. The MC‐AFIR does automated sampling of reaction pathways between two or more reactant molecules. The SC‐AFIR performs automated generation of global or semiglobal reaction path network. The DS‐AFIR finds a single path between given two structures. Exploration of minimum energy structures within the hypersurface in which two different electronic states degenerate, and an interface with the quantum mechanics/molecular mechanics method, are also described. A code termed SAFIRE will also be available, as a visualization software for complicated reaction path networks. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

In this account, a technical overview of the artificial force induced reaction (AFIR) method is presented. The AFIR method is one of automated reaction path search methods developed by the authors, and has been applied extensively to a variety of chemical reactions such as organocatalysis, organometallic catalysis, and photoreactions.There are two modes in the AFIR method, i.e., multi-component mode (MC-AFIR) and

A comprehensive picture of the ultrafast nonradiative decay mechanisms of three cytosine tautomers (amino-keto, imino-keto, and amino-enol forms) is revealed by high-level ab initio potential energy calculations using the multistate (MS) CASPT2 method and also by on-the-fly excited-state molecular dynamics simulations employing the CASSCF method. To obtain a reliable potential energy profile along the deactivation pathways, the MS-CASPT2 method is employed even for the optimization of minimum energy structures in the excited state and conical intersection (CI) structures between the ground and excited states. In the imino (imino-keto) form, we locate a new CI structure involving the twisting of the imino group, and the decay pathway leading to this CI is found to be barrierless, suggesting a remarkably efficient deactivation of imino cytosine. In the keto (amino-keto) form, the MS-CASPT2 calculations exhibit an efficient decay path to the ethylene-like CI involving the twisting of C-C double bond in the six-membered ring, with a barrier of ~0.08 eV from the minimum of the 1 ππ* state. In the enol (amino-enol) form, three types of CIs are identified for the first time. Among them, the ethylene-like CI with a similar molecular structure to the keto form provides the most preferred deactivation pathway of enol cytosine. This pathway exhibits a higher barrier of ~0.22 eV and a higher energy of CI than those of keto cytosine. Nonadiabatic molecular dynamics simulations provide a time-dependent picture of the deactivation processes, including the excited-state lifetime of each tautomer. In particular, the decay time of the imino tautomer is predicted to be only ~100 fs. Our computational results are in remarkably good agreement with the experimental findings in recent femtosecond pump-probe photoionization spectroscopy [J. Am. Chem. Soc. 131, 16939 (2009); J. Phys. Chem. A 115, 8406 (2011)], supporting the coexistence of more than one tautomer in the photophysics of isolated cytosine and that each tautomer exhibits a different excited-state lifetime.1

Hund’s multiplicity rule states that a higher spin state has a lower energy for a given electronic configuration1. Rephrasing this rule for molecular excited states predicts a positive energy gap between spin-singlet and spin-triplet excited states, as has been consistent with numerous experimental observations over almost a century. Here we report a fluorescent molecule that disobeys Hund’s rule and has a negative singlet–triplet energy gap of −11 ± 2 meV. The energy inversion of the singlet and triplet excited states results in delayed fluorescence with short time constants of 0.2 μs, which anomalously decrease with decreasing temperature owing to the emissive singlet character of the lowest-energy excited state. Organic light-emitting diodes (OLEDs) using this molecule exhibited a fast transient electroluminescence decay with a peak external quantum efficiency of 17%, demonstrating its potential implications for optoelectronic devices, including displays, lighting and lasers.

On-the-fly dynamics simulations were carried out using spin-flip time dependent density functional theory (SF-TDDFT) to examine the photoisomerization and photocyclization mechanisms of cis-stilbene following excitation to the ππ* state. A state tracking method was devised to follow the target state among nearly degenerate electronic states during the dynamics simulations. The steepest descent path from the Franck-Condon structure of cis-stilbene in the ππ* state is shown to reach the S1-minimum of 4,4-dihydrophenanthrene (DHP) via a cis-stilbene-like structure (referred to as (S1)cis-min) on a very flat region of the S1-potential energy surface. From the dynamics simulations, the branching ratio of the photoisomerization is calculated as trans:DHP = 35:13, in very good agreement with the experimental data, trans:DHP = 35:10. The discrepancy between the steepest descent pathway and the significant trans-stilbene presence in the branching ratio observed experimentally and herein computationally is clarified from an analysis of geometrical features along the reaction pathway, as well as the low barrier of 0.1 eV for the pathway from (S1)cis-min to the twisted pyramidal structure on the S1-potential energy surface. It is concluded that ππ*-excited cis-stilbene propagates primarily toward the twisted structural region due to dynamic effects, with partial branching to the DHP structural region via the flat-surface region around (S1)cis-min.

The mechanism of ultrafast intersystem crossing in rhenium(I) carbonyl bipyridine halide complexes Re(X)(CO) 3 (bpy) (X = Cl, Br, I) is studied by exploring the structural deformations when going from Franck−Condon (FC) to critical geometries in the low-lying singlet and triplet excited states and by selecting the key vibrational modes. The luminescent decay observed in [Re(Br)(CO) 3 (bpy)] is investigated by means of wavepacket propagations based on the multiconfiguration time-dependent Hartree (MCTDH) method. The dominant coordinates underlying the nonradiative decay process are extracted from minima, minimum energy seam of crossing (MESX) and minimum energy conical intersection (MECI) geometries obtained by the seam model function (SMF)/single-component artificial force induced reaction (SC-AFIR) approach. By choosing the normal modes used in MCTDH from the MECI and MESX geometries, not only the degenerate energy points but also the low-energy-gap regions are included. For this purpose a careful vibrational analysis is performed at each critical geometry and analyzed under the light of the pertinent nonadiabatic coupling terms obtained from the linear vibronic coupling (LVC) model augmented by spin−orbit coupling (SOC) in the electronic diabatic representation.

This article provides an overview on an automated reaction path search method called artificial force induced reaction (AFIR). The AFIR method induces various chemical transformations by applying force between pairs of fragments in a system. By pushing fragments from their various mutual orientations or by applying force between various fragment pairs using the AFIR method, many reaction paths can be explored systematically. In this article, the basic ideas and several different implementations are introduced first. Then, its thoroughness in the automated reaction path search is discussed with its applications to two small molecules. In the later part, its versatility is shown with discussing some previous application examples to organic reaction, organometallic catalysis, photoreaction, surface reaction, phase transition, and enzyme reaction. In addition, an attempt of predicting an idea of new synthesis method from scratch on the basis of the concept quantum chemistry‐aided retrosynthetic analysis (QCaRA) is presented, where the AFIR method was used as a reaction path search engine in QCaRA. Finally, future outlook and a comment on the GRRM program in which the AFIR method is available are given. This article is categorized under: Structure and Mechanism > Reaction Mechanisms and Catalysis Theoretical and Physical Chemistry > Reaction Dynamics and Kinetics