Artificial Force Induced Reaction Research Articles

An Unexpected Ireland-Claisen Rearrangement Cascade During the Synthesis of the Tricyclic Core of Curcusone C: Mechanistic Elucidation by Trial-and-Error and Automatic Artificial Force-Induced Reaction (AFIR) Computations.

In the course of a total synthesis effort directed toward the natural product curcusone C, the Stoltz group discovered an unexpected thermal rearrangement of a divinylcyclopropane to the product of a formal Cope/1,3-sigmatropic shift sequence. Since the involvement of a thermally forbidden 1,3-shift seemed unlikely, theoretical studies involving two approaches, the "trial-and-error" testing of various conceivable mechanisms (Houk group) and an "automatic" approach using the Maeda-Morokuma AFIR method (Morokuma group) were applied to explore the mechanism. Eventually, both approaches converged on a cascade mechanism shown to have some partial literature precedent: Cope rearrangement/1,5-sigmatropic silyl shift/Claisen rearrangement/retro-Claisen rearrangement/1,5-sigmatropic silyl shift, comprising a quintet of five sequential thermally allowed pericyclic rearrangements.

Combined Automated Reaction Pathway Searches and Sparse Modeling Analysis for Catalytic Properties of Lowest Energy Twins of Cu13.

In nanocatalysis, growing attention has recently been given to investigation of energetically low-lying structural isomers of atomic clusters, because some isomers can demonstrate better catalytic activity than the most stable structures. In this study, we present a comparative investigation of catalytic activity for NO dissociation of a pair of the energetically degenerated isomers of Cu13 cluster having C2 and C s symmetries. It is shown that although these isomers have similar structural, electronic, and optical properties, they can possess very different catalytic activities. The effect of isomerization between cluster isomers is considered using state-of-the-art automated reaction pathway search techniques such as an artificial force induced reaction (AFIR) method as a part of a global reaction route mapping (GRRM) strategy. This method allows effectively to locate a large number of possible reaction pathways and transition states (TSs). In total, 12 TSs for NO dissociation were obtained for Cu13, of C2, C s, as well as I h isomers. Sparse modeling analysis shows that LUMO is strongly negatively correlated with total energy of TSs. For most TSs, LUMO has the antibonding character of NO, consisting of the interaction between π* of NO and SOMO of Cu13. Therefore, an increase in the strength of interaction between NO molecule and Cu13 cluster causes the rise in energy of the LUMO, resulting in lowering of the TS energy for NO dissociation. The combination of the automated reaction pathway search technique and sparse modeling represents a powerful tool for analysis and prediction of the physicochemical properties of atomic clusters, especially in the regime of structural fluxionality, where traditional methods based on random geometry search analyses are difficult.

UV-Vis, Fluorescence, and Resonance Raman Spectroscopic and Density Functional Theoretical Studies on 3-Amino-1,2,4-triazole: Microsolvation and Solvent-Dependent Nonadiabatic Excited State Decay in Solution.

The microsolvation and photophysics of 3-amino-1,2,4-triazole (3AT) after excitation to the light-absorbing S2(nπ*) state were studied by using resonance Raman spectroscopy and single component artificial force-induced reaction (SC-AFIR) in a global reaction route mapping (GRRM) strategy. The vibrational spectra were assigned on the basis of experimental data and density functional theory (DFT) calculations. The resonance Raman spectra of 3AT were measured to probe the excited state structural dynamics in the Franck-Condon region. The conformations of 3AT(CH3CN)1, 3AT(CH3OH)2, and 3AT(H2O)2 clusters were determined by combining vibrational spectrum experiments and B3LYP/6-311++G(d,p) computations. DFT calculations were carried out to obtain the minimal excitation energies of the lower-lying singlet excited states, and the curve-crossing points. It was revealed that the short-time structural dynamics of 3AT were dominated by the N-N stretching coordinates. An excited state decay mechanism is proposed: 3AT is initially excited to the S2(nπ*) state, then the conical intersection (CI) of the S2(nπ*)/S1(ππ*) potential energy surfaces is crossed, and 3AT then decays to the lower solvent-dependent excited state S1(ππ*). It subsequently returns to the S0 state, accompanied by a large Stokes fluorescence shift, which was interpreted as the stabilized S1(ππ*) excited state bonding to several water molecules via intermolecular hydrogen bonding.

Excited-State Reactivity of [Mn(im)(CO)3 (phen)]+ : A Structural Exploration.

The electronic excited state reactivity of [Mn(im)(CO)3 (phen)]+ (phen = 1,10-phenanthroline; im = imidazole) ranging between 420 and 330 nm have been analyzed by means of relativistic spin-orbit time-dependent density functional theory and wavefunction approaches (state-average-complete-active-space self-consistent-field/multistate CAS second-order perturbation theory). Minimum energy conical intersection (MECI) structures and connecting pathways were explored using the artificial force induced reaction (AFIR) method. MECIs between the first and second singlet excited states (S1 /S2 -MECIs) were searched by the single-component AFIR (SC-AFIR) algorithm combined with the gradient projection type optimizer. The structural, electronic, and excited states properties of [Mn(im)(CO)3 (phen)]+ are compared to those of the Re(I) analogue [Re(im)(CO)3 (phen)]+ . The high density of excited states and the presence of low-lying metal-centered states that characterize the Mn complex add complexity to the photophysics and open various dissociative channels for both the CO and imidazole ligands. © 2018 Wiley Periodicals, Inc.

Designing the Backbone of Hexasilabenzene Derivatives with a High Unimolecular Kinetic Stability.

It is an important subject to theoretically predict the kinetic stability of transient species. In this study, we have studied the kinetic stability of hexasilabenzene Si6 H6 and its derivatives, that is, decasilanaphthalene Si10 H8 and Li-substituted hexasilabenzene Si6 Li6 , theoretically by the artificial force induced reaction (AFIR) method combined with the rate constant matrix contraction (RCMC) method. Molecular design was further conducted to extend the unimolecular lifetime of hexasilabenzene derivatives. Although both Si10 H8 and Si6 Li6 were shown to possess shorter lifetimes than Si6 H6 , we found that the lifetimes of Si6 Li6 changed depending on arrangements of Li atoms around the monocyclic Si6 backbone. Based on this knowledge, we found that a compound of an atomic composition Si6 H4 Li2 with a planar, monocyclic Si6 backbone has a relatively long unimolecular lifetime. Moreover, substitution of the two Li atoms by Na atoms further increased the lifetime.

DFT Study on Fe(IV)-Peroxo Formation and H Atom Transfer Triggered O2 Activation by NiFe Complex

The mechanism for dioxygen activation using the biomimetic model complex of [NiFe]-hydrogenase, [NiLFe(η5-C5Me5)]+ [L = N,N′-diethyl-3,7-diazanonane-1,9-dithiolato] was established using density functional theory (DFT) and artificial force-induced reaction (AFIR) methods. Our computational results suggest that O2 binds to the FeII center in an end-on fashion and forms a high-valent iron complex, NiFe–peroxo (NiIIFeIV(η2-O2)), which has been experimentally observed. The O–O bond cleavage occurs in the presence of borohydride (BH4–) through hydrogen atom transfer (HAT). Once the HAT occurs, the generated BH3 radical anion (BH3•–) binds to the terminal oxygen of NiFe–OOH, giving rise to BH3OH– and NiIIFeIV═O. The second HAT from BH4– to the oxygen of NiIIFeIV═O leads to BH3OH– and Fe-reduced NiIIFeII complex. Importantly, the dioxygen activation is triggered by HAT, not by proton transfer or hydride transfer. The O2 is activated by the Fe center, and the oxidation state of Fe varies during the process, while...

Global Reaction Route Mapping for Surface Adsorbed Molecules: A Case Study for H2O on Cu(111) Surface

In this study, artificial force induced reaction (AFIR) method, which has been utilized extensively in molecular systems, was applied to surface adsorbed species. A systematic strategy for obtaining path networks for surface adsorbed species is described for future applications to various surface reactions. As a case study, global reaction route mapping on the potential energy surface (PES) of H2O adsorbed on Cu(111) surface was performed.

DFT and AFIR Study on the Mechanism and the Origin of Enantioselectivity in Iron-Catalyzed Cross-Coupling Reactions.

The mechanism of the full catalytic cycle for Fe-chiral-bisphosphine-catalyzed cross-coupling reaction between alkyl halides and Grignard reagents (Nakamura and co-workers, J. Am. Chem. Soc. 2015, 137, 7128) was rationalized by using density functional theory (DFT) and multicomponent artificial force-induced reaction (MC-AFIR) methods. The computed mechanism consists of (a) C-Cl activation, (b) transmetalation, (c) C-Fe bond formation, and (d) C-C bond formation through reductive elimination. Our survey on the prereactant complexes suggested that formation of FeII(BenzP*)Ph2 and FeI(BenzP*)Ph complexes are thermodynamically feasible. FeI(BenzP*)Cl complex is the active intermediate for C-Cl activation. FeII(BenzP*)Ph2 complex can be formed if the concentration of Grignard reagent is high. However, it leads to biphenyl (byproduct) instead of the cross-coupling product. This explains why slow addition of Grignard reagent is critical for the cross-coupling reaction. The MC-AFIR method was used for systematic determination of transition states for C-Fe bond formation and C-C bond formation starting from the key intermediate FeII(BenzP*)PhCl. According to our detailed analysis, C-C bond formation is the selectivity-determining step. The computed enantiomeric ratio of 95:5 is in good agreement with the experimental ratio (90:10). Energy decomposition analysis suggested that the origin of the enantioselectivity is the deformation of Ph-ligand in Fe-complex, which is induced by the bulky tert-butyl group of BenzP* ligand. Our study provides important mechanistic insights for the cross-coupling reaction between alkyl halides and Grignard reagents and guides the design of efficient Fe-based catalysts for cross-coupling reactions.

Copper-Catalyzed Enantioselective Boron Conjugate Addition: DFT and AFIR Study on Different Selectivities of Cu(I) and Cu(II) Catalysts

We present a mechanistic survey on the LCu-catalyzed (L = chiral 2,2′-bipyridine ligand) enantioselective boron conjugate addition reaction, carried out using density functional theory (DFT) and artificial force induced reaction (AFIR) methods. The computed catalytic cycle for Cu(I)- and Cu(II)-based catalysts consists of three steps: (a) boron–boron bond cleavage of B2(pin)2, (b) boron conjugate addition on the β carbon of chalcone, and (c) protonation. The enantioselectivity of the reaction with LCuI or LCuII catalysts is solely governed at the boron conjugate addition step. The multicomponent (MC)-AFIR search and the subsequent DFT calculations for the LCuI catalyst determined transition states (TSs), which lead to CuI-O-enolate and CuI-C-enolate, and both equally contribute to the C–B bond formation with no enantioselectivity. On the other hand, a MC-AFIR search and the subsequent DFT calculations for the analogous LCuII catalyst showed that only the transition state (TS) leading to CuII-O-enolate con...

Global search for low-lying crystal structures using the artificial force induced reaction method: A case study on carbon

We propose an approach to perform the global search for low-lying crystal structures from first principles, by combining the artificial force induced reaction (AFIR) method and the periodic boundary conditions (PBCs). The AFIR method has been applied extensively to molecular systems to elucidate the mechanism of chemical reactions such as homogeneous catalysis. The present PBC/AFIR approach found 274 local minima for carbon crystals in the ${\mathrm{C}}_{8}$ unit cell described by the generalized gradient approximation--Perdew-Burke-Ernzerhof functional. Among many newly predicted structures, three low-lying structures, which exhibit somewhat higher energy compared with those previously predicted, such as $\mathrm{Cco}\ensuremath{-}{\mathrm{C}}_{8}$ ($Z$-carbon) and $M$-carbon, are further discussed with calculations of phonon and band dispersion curves. Furthermore, approaches to systematically explore two- or one-dimensional periodic structures are proposed and applied to the ${\mathrm{C}}_{8}$ unit cell with the slab model. These results suggest that the present approach is highly promising for predicting crystal structures.

Transition-Metal-Free Boryl Substitution Using Silylboranes and Alkoxy Bases

Silylboranes are used as borylation reagents for organohalides in the presence of alkoxy bases without transition-metal catalysts. PhMe2Si–B(pin) reacts with a variety of aryl, alkenyl, and alkyl halides, including sterically hindered examples, to provide the corresponding organoboronates in good yields with high borylation/silylation ratios, showing good functional group compatibility. Halogenophilic attack of a silyl nucleophile on organohalides, and subsequent nucleophilic attack on the boron electrophile are identified to be crucial, based on the results of extensive theoretical and experimental studies. This boryl­ation reaction is further applied to the first direct dimesitylboryl (BMes2) substitution of aryl halides using Ph2MeSi–BMes2 and Na(O-t-Bu), affording aryldimesitylboranes, which are regarded as an important class of compounds for organic materials.1 Introduction2 Boryl Substitution of Organohalides with PhMe2Si–B(pin)/Alkoxy Bases3 Mechanistic Investigations4 DFT Mechanistic Studies Using an Artificial Force Induced Reaction (AFIR) Method5 Dimesitylboryl Substitution of Aryl Halides with Ph2MeSi–BMes2/Na(O-t-Bu)6 Conclusion

Theoretical Study of Addition Reactions of L4M(M = Rh, Ir) and L2M(M = Pd, Pt) to Li+@C60

The addition reaction of M(Cl)(CO)(PPh3)2 (M = Rh, Ir) and M(PPh3)2 (M = Pd, Pt) fragments with X@C60 (X = 0, Li+) were characterized by density functional theory (DFT) and the artificial force-induced reaction (AFIR) method. The calculated free energy profiles suggested that the η2[6:6]-addition is the most favorable reaction, which is consistent with the experimental observations. In the presence of Li+ ion, the reaction is highly exothermic, leading to η2[6:6] product of L4IrLi+@C60. In contrast, an endothermic reaction was observed in the absence of a Li+ ion. The encapsulated Li+ ion can enhance the thermodynamic stability of the η2[6:6] product. The energy decomposition analysis showed that the interaction between metal fragment and X@C60 fragment is the key for the thermodynamic stability. Among the group IA and IIA metal cations, Be2+ encapsulation is the best candidate for the development of new fullerene-transition metal complexes, which will be useful for future potential applications such as solar cells, catalysts, and electronic devices.

Isomerization in Gold Clusters upon O2 Adsorption

A systematic investigation is performed on structural transformations in small neutral gold clusters (Au3–Au12) induced by O2 adsorption, with the use of the fully automated reaction path search techniques, i.e., anharmonic downward distortion following (ADDF) and artificial force induced reaction (AFIR) methods, implemented in the global reaction route mapping (GRRM) program. For each size of Au cluster, the most stable structure, low-energy isomers, and a network of isomerization pathways are determined. The located Aun–O2 adsorption forms can be classified into two groups: η1-AunO2, where only one oxygen atom is adsorbed on Aun, and η2-AunO2, where both oxygen atoms are adsorbed on Aun in a bridged manner. These two adsorption forms can be transformed to each other with a low barrier. The isomerization pathways of gold clusters upon O2 adsorption are compared with those obtained for the pure gold clusters without O2. It is demonstrated that O2 adsorption promotes structural transformations in gold clus...

Artificial Force Induced Reaction Method for Systematic Determination of Complex Reaction Mechanisms.

Nowadays, computational studies are very important for the elucidation of reaction mechanisms and selectivity of complex reactions. However, traditional computational methods usually require an estimated reaction path, mainly driven by limited experimental implications, intuition, and assumptions of stationary points. However, the artificial force induced reaction (AFIR) method in the global reaction route mapping (GRRM) strategy can be used for unbiased and automatic reaction path searches for complex reactions. In this account, we highlight applications of the AFIR method to a variety of reactions (organic, organometallic, enzymatic, and photochemical) of complex molecular systems. In addition, the AFIR method has been successfully used to rationalise the origin of stereo- and regioselectivity. The AFIR method can be applied from small to large molecular systems, and will be a very useful tool for the study of complex molecular problems in many areas of chemistry, biology, and material sciences.

Artificial Force Induced Reaction (AFIR) Method for Exploring Quantum Chemical Potential Energy Surfaces.

In this account, a technical overview of the artificial force induced reaction (AFIR) method is presented. The AFIR method is one of the 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., a multicomponent mode and a single-component mode. The former has been applied to bimolecular and multicomponent reactions and the latter to unimolecular isomerization and dissociation reactions. Five numerical examples are presented for an Aldol reaction, a Claisen rearrangement, a Co-catalyzed hydroformylation, a fullerene structure search, and a nonradiative decay path search in an electronically excited naphthalene molecule. Finally, possible applications of the AFIR method are discussed.

Computational SN 2-Type Mechanism for the Difluoromethylation of Lithium Enolate with Fluoroform through Bimetallic C-F Bond Dual Activation.

The reaction mechanism for difluoromethylation of lithium enolates with fluoroform was analyzed computationally (DFT calculations with the artificial force induced reaction (AFIR) method and solvation model based on density (SMD) solvation model (THF)), showing an SN 2-type carbon-carbon bond formation; the "bimetallic" lithium enolate and lithium trifluoromethyl carbenoid exert the C-F bond "dual" activation, in contrast to the monometallic butterfly-shaped carbenoid in the Simmons-Smith reaction. Lithium enolates, generated by the reaction of 2 equiv. of lithium hexamethyldisilazide (rather than 1 or 3 equiv.) with the cheap difluoromethylating species fluoroform, are the most useful alkali metal intermediates for the synthesis of pharmaceutically important α-difluoromethylated carbonyl products.

ChemInform Abstract: Computational Catalysis Using the Artificial Force Induced Reaction Method

AbstractReview: 33 refs.

Computational Catalysis Using the Artificial Force Induced Reaction Method.

The artificial force induced reaction (AFIR) method in the global reaction route mapping (GRRM) strategy is an automatic approach to explore all important reaction paths of complex reactions. Most traditional methods in computational catalysis require guess reaction paths. On the other hand, the AFIR approach locates local minima (LMs) and transition states (TSs) of reaction paths without a guess, and therefore finds unanticipated as well as anticipated reaction paths. The AFIR method has been applied for multicomponent organic reactions, such as the aldol reaction, Passerini reaction, Biginelli reaction, and phase-transfer catalysis. In the presence of several reactants, many equilibrium structures are possible, leading to a number of reaction pathways. The AFIR method in the GRRM strategy determines all of the important equilibrium structures and subsequent reaction paths systematically. As the AFIR search is fully automatic, exhaustive trial-and-error and guess-and-check processes by the user can be eliminated. At the same time, the AFIR search is systematic, and therefore a more accurate and comprehensive description of the reaction mechanism can be determined. The AFIR method has been used for the study of full catalytic cycles and reaction steps in transition metal catalysis, such as cobalt-catalyzed hydroformylation and iron-catalyzed carbon-carbon bond formation reactions in aqueous media. Some AFIR applications have targeted the selectivity-determining step of transition-metal-catalyzed asymmetric reactions, including stereoselective water-tolerant lanthanide Lewis acid-catalyzed Mukaiyama aldol reactions. In terms of establishing the selectivity of a reaction, systematic sampling of the transition states is critical. In this direction, AFIR is very useful for performing a systematic and automatic determination of TSs. In the presence of a comprehensive description of the transition states, the selectivity of the reaction can be calculated more accurately. For relatively large molecular systems, the computational cost of AFIR searches can be reduced by using the ONIOM(QM:QM) or ONIOM(QM:MM) methods. In common practice, density functional theory (DFT) with a relatively small basis set is used for the high-level calculation, while a semiempirical approach or a force field description is used for the low-level calculation. After approximate LMs and TSs are determined, standard computational methods (e.g., DFT with a large basis set) are used for the full molecular system to determine the true LMs and TSs and to rationalize the reaction mechanism and selectivity of the catalytic reaction. The examples in this Account evidence that the AFIR method is a powerful approach for accurate prediction of the reaction mechanisms and selectivities of complex catalytic reactions. Therefore, the AFIR approach in the GRRM strategy is very useful for computational catalysis.

Palladium-catalyzed regioselective and stereo-invertive ring-opening borylation of 2-arylaziridines with bis(pinacolato)diboron: experimental and computational studies.

A palladium catalyzed regioselective borylative ring opening reaction of 2-arylaziridines to give β-amino-β-arylethylborates was developed. The reaction reported herein represents the first example of ring-opening borylation of non-vinylic aziridines and direct borylative C(sp3)-N bond cleavage of neutral organic substrates. NMR studies and density functional theory (DFT) calculations suggested that the active intermediate for the reaction is a PdL2 complex [L = P(t-Bu)2Me]. The multi-component artificial force-induced reaction method (MC-AFIR) located the transition states for the regioselectivity-determining aziridine ring opening that proceeds in an SN2 fashion, and explained the selectivity of the reaction. The full catalytic cycle consists of a selectivity-determining aziridine ring opening (oxidative addition), a proton transfer, phosphine ligand dissociation from the catalyst, boron-boron bond cleavage, and reductive elimination. Water is important to the drive the transmetalation step. The calculated overall mechanism and selectivity are consistent with the experimental results.

Mechanisms of Hydrogen Generation from Tetrameric Clusters of Lithium Amidoborane.

The first-principles study of dehydrogenation mechanism of tetrameric clusters of lithium amidoborane LiNH2BH3, (LiAB)4, is presented. The choice of tetramer is based on the suspicion that dimeric cluster models used in previous theoretical studies are too small to capture the essence of the reaction. Dehydrogenation pathways starting from three isomers of (LiAB)4 tetramers were explored by applying the artificial force induced reaction (AFIR) method at the M06 level of theory. All obtained reaction pathways feature initial dimerization of two LiAB molecules in the tetramer. Formation of intermediates containing the Li3H moiety is a very characteristic feature of all pathways. In the succeeding rate-limiting step of the release of H2 molecule, a hydridic H atom of the Li3H moiety activates a protic H atom of the NH2 group with formation of the Li2H2 moiety in transition state. The most kinetically favorable pathway has the activation enthalpy of 26.6 kcal mol(-1), substantially lower than that found for dimeric cluster. The obtained results suggest that only three LiAB molecules directly participate in the elementary reactions.