Implementation and performance of the artificial force induced reaction method in the GRRM17 program

Tetraborylation of p-Benzynes Generated by the Masamune-Bergman Cyclization through Reaction Design Based on the Reaction Path Network

Designing the reactant molecule of an initial reaction, based on quantum chemical pathway exploration, enabled us to access a new reaction, i.e., the tetraborylation reaction of p-benzynes generated from 1,2-diethynylbenzene derivatives, using bis(pinacolato)diborane(4) (B2pin2). Based on the reaction path network generated via the artificial-force-induced reaction (AFIR) method, desired and undesired paths were identified and used to modify the chemical structure of the reactant. After the in silico screening, the optimal structure of the reactant was determined to be a 1,2-diethynylbenzene derivative with a butylene linker. The reaction of the optimized reactant and its derivatives with an excess of B2pin2 gave the tetraborylated products in good yields (up to 58%). It is quite intriguing that the two carbons of p-benzyne behave formally as dicarbenes in this reaction.

Photoinduced dual bond rotation of a nitrogen-containing system realized by chalcogen substitution

Photoinduced concerted multiple-bond rotation has been proposed in some biological systems. However, the observation of such phenomena in synthetic systems, in other words, the synthesis of molecules that undergo photoinduced multiple-bond rotation upon photoirradiation, has been a challenge in the photochemistry field. Here we describe a chalcogen-substituted benzamide system that exhibits photoinduced dual bond rotation in heteroatom-containing bonds. Introduction of the chalcogen substituent into a sterically hindered benzamide system provides sufficient kinetic stability and photosensitivity to enable the photoinduced concerted rotation. The presence of two different substituents on the phenyl ring in the thioamide derivative enables the generation of a pair of enantiomers and E/Z isomers. Using these four stereoisomers as indicators of which bonds are rotated, we monitor the photoinduced C-N/C-C concerted bond rotation in the thioamide derivative depending on external stimuli such as temperature and photoirradiation. Theoretical calculations provide insight on the mechanism of this selective photoinduced C-N/C-C concerted rotation.

An encompassed representation of timescale hierarchies in first-order reaction network

Complex networks are pervasive in various fields such as chemistry, biology, and sociology. In chemistry, first-order reaction networks are represented by a set of first-order differential equations, which can be constructed from the underlying energy landscape. However, as the number of nodes increases, it becomes more challenging to understand complex kinetics across different timescales. Hence, how to construct an interpretable, coarse-graining scheme that preserves the underlying timescales of overall reactions is of crucial importance. Here, we develop a scheme to capture the underlying hierarchical subsets of nodes, and a series of coarse-grained (reduced-dimensional) rate equations between the subsets as a function of time resolution from the original reaction network. Each of the coarse-grained representations guarantees to preserve the underlying slow characteristic timescales in the original network. The crux is the construction of a lumping scheme incorporating a similarity measure in deciphering the underlying timescale hierarchy, which does not rely on the assumption of equilibrium. As an illustrative example, we apply the scheme to four-state Markovian models and Claisen rearrangement of allyl vinyl ether (AVE), and demonstrate that the reduced-dimensional representation accurately reproduces not only the slowest but also the faster timescales of overall reactions although other reduction schemes based on equilibrium assumption well reproduce the slowest timescale but fail to reproduce the second-to-fourth slowest timescales with the same accuracy. Our scheme can be applied not only to the reaction networks but also to networks in other fields, which helps us encompass their hierarchical structures of the complex kinetics over timescales.

Size-dependent reactivity of Rh cationic clusters to reduce NO by CO in the gas phase at high temperatures

The reactivity of the reduction of NO pre-adsorbed on Rh2-9+ clusters by CO was investigated using a combination of an alternate on-off gas injection method and thermal desorption spectrometry. The reduction of RhnNxOy+ clusters by CO was evaluated by varying the CO concentration at T = 903 K. Among the RhnNxOx+ clusters, the Rh3N2O2+ cluster exhibited the highest reduction activity, whereas the other clusters, Rh2,4-9NxOx+, showed lower reactivity. Density functional theory (DFT) calculations for Rh3+ and Rh6+ revealed that the rate-determining step for NO reduction in the presence of CO was NO bond dissociation through the kinetics analysis using the RRKM theory. The reduction of Rh3N2O2+ is kinetically preferable to that of Rh6N2O2+. The DFT results were in qualitative agreement with the experimental results.

Azobenzene as a photoswitchable mechanophore

Azobenzene has been widely explored as a photoresponsive element in materials science. Although some studies have investigated the force-induced isomerization of azobenzene, the effect of force on the rupture of azobenzene has not been explored. Here we show that the light-induced structural change of azobenzene can also alter its rupture forces, making it an ideal light-responsive mechanophore. Using single-molecule force spectroscopy and ultrasonication, we found that cis and trans para-azobenzene isomers possess contrasting mechanical properties. Dynamic force spectroscopy experiments and quantum-chemical calculations in which azobenzene regioisomers were pulled from different directions revealed that the distinct rupture forces of the two isomers are due to the pulling direction rather than the energetic difference between the two isomers. These mechanical features of azobenzene can be used to rationally control the macroscopic fracture behaviours of polymer networks by photoillumination. The use of light-induced conformational changes to alter the mechanical response of mechanophores provides an attractive way to engineer polymer networks of light-regulatable mechanical properties.

Using Mechanochemistry to Activate Commodity Plastics as Initiators for Radical Chain Reactions of Small Organic Molecules

Radical initiators such as azo compounds and organic peroxides have been widely used to facilitate numerous transformations of free radicals, which enable the efficient synthesis of structurally complex molecules, natural products, polymers, and functional materials. However, these high-energy reagents are potentially explosive and thus often require special precautions or delicate operating conditions. We postulated that a more convenient and safer alternative for radical chain initiation could be developed by mechanical activation of thermodynamically stable covalent bonds. Here, we show that commodity plastics such as polyethylene and poly(vinyl acetate) are capable of acting as efficient initiators for radical chain reactions under solvent-free mechanochemical conditions. In this approach, polymeric mechanoradicals, which are generated by homolytic cleavage of the polymer chains in response to the applied mechanical energy provided by ball milling, react with tris(trimethylsilyl)silane to initiate radical chain dehalogenation of organic halides. Preliminary calculations support our proposed force-induced radical chain mechanism.

Ultrafast Raman observation of the perpendicular intermediate phantom state of stilbene photoisomerization

Trans-cis photoisomerization is generally described by a model in which the reaction proceeds via a common intermediate having a perpendicular conformation around the rotating bond, irrespective of from which isomer the reaction starts. Nevertheless, such an intermediate has yet to be identified unambiguously, and it is often called the 'phantom' state. Here we present the structural identification of the common, perpendicular intermediate of stilbene photoisomerization using ultrafast Raman spectroscopy. Our results reveal ultrafast birth and decay of an identical, short-lived transient that exhibits a vibrational signature characteristic of the perpendicular state upon photoexcitation of the trans and cis forms. In combination with ab initio molecular dynamics simulations, it is shown that the photoexcited trans and cis forms are funnelled off to the ground state through the same, perpendicular intermediate.

Effect of Internal Substituents on the Properties of Dibenzo[g,p]chrysene

We investigated the chemical and physical properties of internally functionalized dibenzo[g,p]chrysene (DBC) derivatives. These molecules exhibit chiral double-helicene-like structures that are configurationally stable at ambient temperatures. The internal substituents control the conformational change in the excited state, thereby modulating the emission intensity. Notably, the DBC derivative with a methylenedioxy unit undergoes aromatization through elimination of the internal substituent upon photoexcitation, resulting in the formation of DBC.

Quantum Chemical Calculations of Flexible Tripeptide-Ni(II) Ion-Mediated Supramolecular Fragments and Comparative Analysis of Tripeptide Complexes with Various Metal(II) Ions

An automated conformational search method was employed to efficiently determine the stable conformers and weak hydrogen bonds of a flexible tripeptide coordinated with a solitary metal(II) ion in an aqueous environment. Quantum chemical calculations were performed to investigate the tendency of octahedral coordination formation between different metal(II) ions and various coordination models (ammonia molecule, chelate molecule, and flexible tripeptide). The octahedral coordination was analyzed by decomposing it into tridentate, bidentate, and monodentate coordination model complexes to assess their formation propensities and conformational properties. Additionally, population analysis, including electrostatic potential mapping and natural population analysis, was performed to identify the unique properties of the Ni(II) ion in forming octahedral coordination in crystals and to explore the potential of other metal(II) ions for self-assembling novel coordination configurations in peptide-metal compounds. Two common hydrogen bonding interactions were examined by using artificial forces to facilitate dissociation or reinforcement.

Efficient Acceleration of Reaction Discovery in the Ab Initio Nanoreactor: Phenyl Radical Oxidation Chemistry

Over the years, many computational strategies have been employed to elucidate reaction networks. One of these methods is accelerated molecular dynamics, which can circumvent the expense required in dynamics to find all reactants and products (local minima) and transition states (first-order saddle points) on a potential energy surface (PES) by using fictitious forces that promote reaction events. The ab initio nanoreactor uses these accelerating forces to study large chemical reaction networks from first-principles quantum mechanics. In the initial nanoreactor studies, this acceleration was done through a piston periodic compression potential, which pushes molecules together to induce entropically unfavorable bimolecular reactions. However, the piston is not effective for discovering intramolecular and dissociative reactions, such as those integral to the decomposition channels of phenyl radical oxidation. In fact, the choice of accelerating forces dictates not only the rate of reaction discovery but also the types of reactions discovered; thus, it is critical to understand the biases and efficacies of these forces. In this study, we examine forces using metadynamics, attractive potentials, and local thermostats for accelerating reaction discovery. For each force, we construct a separate phenyl radical combustion reaction network using solely that force in discovery trajectories. We elucidate the enthalpic and entropic trends of each accelerating force and highlight their efficiency in reaction discovery. Comparing the nanoreactor-constructed reaction networks with literature renditions of the phenyl radical combustion PES shows that a combination of accelerating forces is best suited for reaction discovery.

Skeletal rearrangement of a boron-containing annulenic molecule into a macrocycle bridged by an electronically stabilized boron cation

An annulenic molecule containing a three-coordinate chloroborane moiety, which exhibits a borane-olefin proximity effect, undergoes a skeletal rearrangement upon chloride abstraction, to generate a three-dimensional macrocyclic molecule featuring a borocenium (η5-cyclopentadienyl-B+-R) structure.

Quantum chemical calculations for reaction prediction in the development of synthetic methodologies

Quantum chemical calculations have been used in the development of synthetic methodologies to analyze the reaction mechanisms of the developed reactions. Their ability to estimate chemical reaction pathways, including transition state energies and connected equilibria, has led researchers to embrace their use in predicting unknown reactions. This perspective highlights strategies that leverage quantum chemical calculations for the prediction of reactions in the discovery of new methodologies. Selected examples demonstrate how computation has driven the development of unknown reactions, catalyst design, and the exploration of synthetic routes to complex molecules prior to often laborious, costly, and time-consuming experimental investigations.

An energy decomposition and extrapolation scheme for evaluating electron transfer rate constants: a case study on electron self-exchange reactions of transition metal complexes

A simple approach to the analysis of electron transfer (ET) reactions based on energy decomposition and extrapolation schemes is proposed. The present energy decomposition and extrapolation-based electron localization (EDEEL) method represents the diabatic energies for the initial and final states using the adiabatic energies of the donor and acceptor species and their complex. A scheme for the efficient estimation of ET rate constants is also proposed. EDEEL is semi-quantitative by directly evaluating the seam-of-crossing region of two diabatic potentials. In a numerical test, EDEEL successfully provided ET rate constants for electron self-exchange reactions of thirteen transition metal complexes with reasonable accuracy. In addition, its energy decomposition and extrapolation schemes provide all the energy values required for activation-strain model (ASM) analysis. The ASM analysis using EDEEL provided rational interpretations of the variation of the ET rate constants as a function of the transition metal complexes. These results suggest that EDEEL is useful for efficiently evaluating ET rate constants and obtaining a rational understanding of their magnitudes.

A Combined Reaction Path Search and Hybrid Solvation Method for the Systematic Exploration of Elementary Reactions at the Solid-Liquid Interface

We present a combined simulation method of single-component artificial force induced reaction (SC-AFIR) and effective screening medium combined with the reference interaction site model (ESM-RISM), termed SC-AFIR+ESM-RISM. SC-AFIR automatically and systematically explores the chemical reaction pathway, and ESM-RISM directly simulates the precise electronic structure at the solid-liquid interface. Hence, SC-AFIR+ESM-RISM enables us to explore reliable reaction pathways at the solid-liquid interface. We applied it to explore the dissociation pathway of an H2O molecule at the Cu(111)/water interface. The reaction path networks of the whole reaction and the minimum energy paths from H2O to H2 + O depend on the interfacial environment. The qualitative difference in the energy diagrams and the resulting change in the kinematically favored dissociation pathway upon changing the solvation environments are discussed. We believe that SC-AFIR+ESM-RISM will be a powerful tool to reveal the details of chemical reactions in surface catalysis and electrochemistry.

A systematic mechanistic survey on the reactions between OH radical and CH3OH on ice

A systematic mechanistic survey was performed for the CH3OH + OH reaction on ice. ONIOM(ωB97X-D/Def2-TZVP:AMOEBA09) calculations suggested a range of binding energies for the CH2OH radical (0.29-0.69 eV) and CH3OH (0.15-0.72 eV) molecule on hexagonal water ice (Ih) and amorphous solid water (ASW). Computed average binding energies of CH2OH radical (0.49 eV) and CH3OH (0.41 eV) are relatively stronger compared to the CH3O radical binding energies (0.32 eV, Sameera et al., J. Phy. Chem. A, 2021, 125, 387-393). Thus, the CH3OH molecule, CH2OH and CH3O radicals can adsorb on ice, where the binding energies follow the order CH2OH > CH3OH > CH3O. The multi-component artificial force-induced reaction (MC-AFIR) method systematically determined the reaction mechanisms for the CH3OH + OH reaction on ice, where two reaction paths, giving rise to CH2OH and CH3O radicals, were confirmed. A range of reaction barriers, employing the ωB97X-D/Def2-TZVP level of theory, was found for each reaction (0.03-0.11 eV for CH2OH radical formation, and 0.03-0.44 eV for CH3O radical formation). Based on the lowest energy reaction paths, we suspect that both reactions operate on ice. The computed data in this study evidence that the nature of the binding site or the reaction site has a significant effect on the computed binding energies or reaction barriers. Thus, the outcomes of the present study will be very useful for the computational astrochemistry community to determine reliable binding energies and reaction barriers on ice.

Reproducing the Reaction Route Map on the Shape Space from Its Quotient by the Complete Nuclear Permutation-Inversion Group

This study develops an algorithm to reproduce reaction route maps (RRMs) in the shape space from the outputs of potential search algorithms. To demonstrate the algorithm, global reaction route mapping is utilized as a potential search algorithm, but the proposed algorithm should work with other potential search algorithms in principle. The proposed algorithm does not require any encoding of the molecular configurations and is thus applicable to complicated realistic molecules for which efficient encoding is not readily available. We show that subgraphs of an RRM mapped to each other by the action of the symmetry group are isomorphic and also provide an algorithm to compute the set of feasible transformations in the sense of Longuet-Higgins. We demonstrate the proposed algorithm in toy models and in more realistic molecules. Finally, we remark on absolute rate theory from our perspective.

Designing main-group catalysts for low-temperature methane combustion by ozone

The catalytic combustion of methane at a low temperature is becoming increasingly key to controlling unburned CH₄ emissions from natural gas vehicles and power plants, although the low activity of benchmark platinum-group-metal catalysts hinders its broad application. Based on automated reaction route mapping, we explore main-group elements catalysts containing Si and Al for low-temperature CH₄ combustion with ozone. Computational screening of the active site predicts that strong Brønsted acid sites are promising for methane combustion. We experimentally demonstrate that catalysts containing strong Bronsted acid sites exhibit improved CH₄ conversion at 250 °C, correlating with the theoretical predictions. The main-group catalyst (proton-type beta zeolite) delivered a reaction rate that is 442 times higher than that of a benchmark catalyst (5 wt% Pd-loaded Al₂O₃) at 190 °C and exhibits higher tolerance to steam and SO₂. Our strategy demonstrates the rational design of earth-abundant catalysts based on automated reaction route mapping.

Reaction mechanism - explored with the unified reaction valley approach

One of the ultimate goals of chemistry is to understand and manipulate chemical reactions, which implies the ability to monitor the reaction and its underlying mechanism at an atomic scale. In this article, we introduce the Unified Reaction Valley Approach (URVA) as a tool for elucidating reaction mechanisms, complementing existing computational procedures. URVA combines the concept of the potential energy surface with vibrational spectroscopy and describes a chemical reaction via the reaction path and the surrounding reaction valley traced out by the reacting species on the potential energy surface on their way from the entrance to the exit channel, where the products are located. The key feature of URVA is the focus on the curving of the reaction path. Moving along the reaction path, any electronic structure change of the reacting species is registered by a change in the normal vibrational modes spanning the reaction valley and their coupling with the path, which recovers the curvature of the reaction path. This leads to a unique curvature profile for each chemical reaction, with curvature minima reflecting minimal change and curvature maxima indicating the location of important chemical events such as bond breaking/formation, charge polarization and transfer, rehybridization, etc. A decomposition of the path curvature into internal coordinate components or other coordinates of relevance for the reaction under consideration, provides comprehensive insight into the origin of the chemical changes taking place. After giving an overview of current experimental and computational efforts to gain insight into the mechanism of a chemical reaction and presenting the theoretical background of URVA, we illustrate how URVA works for three diverse processes, (i) [1,3] hydrogen transfer reactions; (ii) α-keto-amino inhibitor for SARS-CoV-2 Mpro; (iii) Rh-catalyzed cyanation. We hope that this article will inspire our computational colleagues to add URVA to their repertoire and will serve as an incubator for new reaction mechanisms to be studied in collaboration with our experimental experts in the field.

Computational survey of humin formation from 5-(hydroxymethyl)furfural under basic conditions

A comprehensive reaction-path search for the oligomerization of 5-(hydroxymethyl)furfural (HMF) based on quantum chemical calculations was conducted to clarify the mechanism of humin formation in the oxidation of HMF to furan-2,5-dicarboxylic acid (FDCA), in which humin is a typical macromolecular byproduct. The present procedure repeatedly utilizes the multi-component artificial-force-induced reaction (MC-AFIR) method to investigate multistep oligomerization reactions. Although humin formation has been reported even in reagent-grade HMFs with 97-99% purity during their storage at low temperatures, no direct addition path of two HMFs with <185 kJ mol^-1 barrier has been found, suggesting humin formation is caused by a reaction with impurities. Based on the reaction conditions, we considered the reactions of HMF + H₂O, HMF + OH^-, and HMF + O₂ and identified three reaction paths with <65 kJ mol^-1 barrier for the reaction of HMF + OH^-. Further, the suppression of humin formation by the acetal protection of HMF is computationally confirmed.

Real-time observation of the Woodward-Hoffmann rule for 1,3-cyclohexadiene by femtosecond soft X-ray transient absorption

The stereochemistry of pericyclic reactions is explained by orbital symmetry conservation, referred to as the Woodward-Hoffmann (WH) rule. Although this rule has been verified using the structures of reactants and products, the temporal evolution of the orbital symmetry during the reaction has not been clarified. Herein, we used femtosecond soft X-ray transient absorption spectroscopy to elucidate the thermal pericyclic reaction of 1,3-cyclohexadiene (CHD) molecules, i.e., their isomerization to 1,3,5-hexatriene. In the present experimental scheme, the ring-opening reaction is driven by the thermal vibrational energy induced by photoexcitation to the Rydberg states at 6.2 eV and subsequent femtosecond relaxation to the ground state of CHD molecules. The direction of the ring opening, which can be conrotatory or disrotatory, was the primary focus, and the WH rule predicts the disrotatory pathway in the thermal process. We observed the shifts in K-edge absorption of the carbon atom from the 1s orbital to vacant molecular orbitals around 285 eV at a delay between 340 and 600 fs. Furthermore, a theoretical investigation predicts that the shifts depend on the molecular structures along the reaction pathways and the observed shifts in induced absorption are attributed to the structural change in the disrotatory pathway. This confirms that the orbital symmetry is dynamically conserved in the ring-opening reaction of CHD molecules as predicted using the WH rule.

An automated reaction route mapping for the reaction of NO and active species on Ag4 clusters in zeolites

A computational investigation of the catalytic reaction on multinuclear sites is very challenging. Here, using an automated reaction route mapping method, the single-component artificial force induced reaction (SC-AFIR) algorithm, the catalytic reaction of NO and OH/OOH species over the Ag₄²⁺ cluster in a zeolite is investigated. The results of the reaction route mapping for H₂ + O₂ reveal that OH and OOH species are formed over the Ag₄²⁺ cluster via an activation barrier lower than that of OH formation from H₂O dissociation. Then, reaction route mapping is performed to examine the reactivity of the OH and OOH species with NO molecules over the Ag₄²⁺ cluster, resulting in the facile reaction path of HONO formation. With the aid of the automated reaction route mapping, the promotion effect of H₂ addition on the SCR reaction was computationally proposed (boosting the formation of OH and OOH species). In addition, the present study emphasizes that automated reaction route mapping is a powerful tool to elucidate the complicated reaction pathway on multi-nuclear clusters.

Chemical modification of dimethylpolysiloxane for enhancement of CO2 binding enthalpy

The intermittent increase in CO₂ concentration in the atmosphere is a serious problem that contributes to climate change; the combustion of fossil fuels produces the majority of CO₂, and technology is needed to capture it efficiently. Various CO₂ capture materials have been developed so far. Membrane separation, in particular, has an advantage over other capture technologies due to its ease of use. Poly(dimethylsiloxane) (PDMS) has been widely used as a membrane material for CO₂ capture because of its high gas permeability. However, despite their high CO₂ permeance, PDMS membranes are still in their infancy, especially regarding CO₂ selectivity due to the weak interaction between CO₂ and PDMS. Here we evaluated the CO₂ interaction with the PDMS chain at the atomic scale and attempted to improve the CO₂ affinity of the PDMS chain using density functional theory (DFT). Specifically, we substituted elements in the Si-O framework with other elements and substituted the methyl groups with other chemical groups, and incorporated metallic elements such as Mg and Ti. All the chemical modifications by main group elements resulted in physisorption, but chemisorption of CO₂ was observed in PDMS incorporating metallic elements. Since several modes of CO₂ binding were observed in PDMS with incorporated metal elements, the binding enthalpy and binding mode were analyzed. As a result of various chemical modifications, it was found that introducing earth metal elements into PDMS was the most effective way to enhance the interaction between PDMS and CO₂.

Early-Stage Formation of the SIFSIX-3-Zn Metal-Organic Framework: An Automated Computational Study

Metal-organic frameworks (MOFs) have attracted significant attention over the past 2 decades due to their wide applicability as functional materials. However, targeted synthesis of novel MOFs remains problematic as their formation mechanisms are poorly understood, which forces us to rely on serendipity in the synthesis of novel MOFs. Here, we demonstrate a workflow employing the artificial force induced reaction (AFIR) method to investigate the self-assembly process of the node of the SIFSIX-3-Zn MOF, [Zn(pyz)₄(SiF₆)₂]^2- (pyz = pyrazine), in an automated manner. The workflow encompassing AFIR calculations, generation of extensive reaction path networks, propagation simulations of intermediates, and further refinements of identified formation pathways showed that the nodal structure can form through multiple competing pathways involving interconvertible intermediates. This finding provides a plausible rationale for the stochastic multistage processes believed to be key in MOF formation. Furthermore, this work represents the first application of an automated reaction mechanism discovery method to a MOF system using a general workflow that is applicable to study the formation of other MOF motifs as well.

Dinaphthooxepine Bisimide Undergoes Oxygen Extrusion Reaction upon Electron Injection at Room Temperature

We report the synthesis and properties of a dinaphthooxepine bisimide (DNOBI), a nonplanar perylene bisimide (PBI) analogue with an inserted oxygen atom. A DNOBI underwent an oxygen-extrusion reaction smoothly upon electron injection at room temperature, affording PBI in good yield. Studies on the reaction mechanism suggest that the injection of two electrons triggers the isomerization of DNOBI to dinaphthooxanorcaradiene bisimide, which is a key step in inducing the oxygen-extrusion reaction.

Oxidation and Reduction Pathways in the Knowles Hydroamination via a Photoredox-Catalyzed Radical Reaction

Systematic reaction path exploration revealed the entire mechanism of Knowles's light-promoted catalytic intramolecular hydroamination. Bond formation/cleavage competes with single electron transfer (SET) between the catalyst and substrate. These processes are described by adiabatic processes through transition states in an electronic state and non-radiative transitions through the seam of crossings (SX) between different electronic states. This study determined the energetically favorable SET path by introducing a practical computational model representing SET as non-adiabatic transitions via SXs between substrate's potential energy surfaces for different charge states adjusted based on the catalyst's redox potential. Calculations showed that the reduction and proton shuttle process proceeded concertedly. Also, the relative importance of SET paths (giving the product and leading back to the reactant) varies depending on the catalyst's redox potential, affecting the yield.

Multi-state Energy Landscape for Photoreaction of Stilbene and Dimethyl-stilbene

We have recently developed the reaction space projector (ReSPer) method, which constructs a reduced-dimensionality reaction space uniquely determined from reference reaction paths for a polyatomic molecular system and projects classical trajectories into the same reaction space. In this paper, we extend ReSPer to the analysis of photoreaction dynamics and relaxation processes of stilbene and present the concept of a "multi-state energy landscape," incorporating the ground- and excited-state reaction subspaces. The multi-state energy landscape successfully explains the previously established photoreaction processes of cis-stilbene, such as the cis-trans photoisomerization and photocyclization. In addition, we discuss the difference in the excited-state reaction dynamics between stilbene and 1,1'-dimethyl stilbene based on a common reaction subspace determined from the framework part of reference structures with different number of atoms. This approach allows us to target any molecule with a common framework, greatly expanding the applicability of the ReSPer analysis. The multi-state energy landscape provides fruitful insight into photochemical reactions, exploring the excited- and ground-state potential energy surfaces, as well as comprehensive reaction processes with nonradiative transitions between adiabatic states, within the stage of a reduced-dimensionality reaction space.

ReNeGate: A Reaction Network Graph-Theoretical Tool for Automated Mechanistic Studies in Computational Homogeneous Catalysis

Exploration of the chemical reaction space of chemical transformations in multicomponent mixtures is one of the main challenges in contemporary computational chemistry. To remove expert bias from mechanistic studies and to discover new chemistries, an automated graph-theoretical methodology is proposed, which puts forward a network formalism of homogeneous catalysis reactions and utilizes a network analysis tool for mechanistic studies. The method can be used for analyzing trajectories with single and multiple catalytic species and can provide unique conformers of catalysts including multinuclear catalyst clusters along with other catalytic mixture components. The presented three-step approach has the integrated ability to handle multicomponent catalytic systems of arbitrary complexity (mixtures of reactants, catalyst precursors, ligands, additives, and solvents). It is not limited to predefined chemical rules, does not require prealignment of reaction mixture components consistent with a reaction coordinate, and is not agnostic to the chemical nature of transformations. Conformer exploration, reactive event identification, and reaction network analysis are the main steps taken for identifying the pathways in catalytic systems given the starting precatalytic reaction mixture as the input. Such a methodology allows us to efficiently explore catalytic systems in realistic conditions for either previously observed or completely unknown reactive events in the context of a network representing different intermediates. Our workflow for the catalytic reaction space exploration exclusively focuses on the identification of thermodynamically feasible conversion channels, representative of the (secondary) catalyst deactivation or inhibition paths, which are usually most difficult to anticipate based solely on expert chemical knowledge. Thus, the expert bias is sought to be removed at all steps, and the chemical intuition is limited to the choice of the thermodynamic constraint imposed by the applicable experimental conditions in terms of threshold energy values for allowed transformations. The capabilities of the proposed methodology have been tested by exploring the reactivity of Mn complexes relevant for catalytic hydrogenation chemistry to verify previously postulated activation mechanisms and unravel unexpected reaction channels relevant to rare deactivation events.

Guided discovery of chemical reaction pathways with imposed activation

Computational power and quantum chemical methods have improved immensely since computers were first applied to the study of reactivity, but the de novo prediction of chemical reactions has remained challenging. We show that complex reaction pathways can be efficiently predicted in a guided manner using chemical activation imposed by geometrical constraints of specific reactive modes, which we term imposed activation (IACTA). Our approach is demonstrated on realistic and challenging chemistry, such as a triple cyclization cascade involved in the total synthesis of a natural product, a water-mediated Michael addition, and several oxidative addition reactions of complex drug-like molecules. Notably and in contrast with traditional hand-guided computational chemistry calculations, our method requires minimal human involvement and no prior knowledge of the products or the associated mechanisms. We believe that IACTA will be a transformational tool to screen for chemical reactivity and to study both by-product formation and decomposition pathways in a guided way.

Reactive orbital energy theory serving a theoretical foundation for the electronic theory of organic chemistry

It is established that the reactive orbital energy theory (ROET) theoretically reproduces the rule-based electronic theory diagrams of organic chemistry by a comparative study on the charge transfer natures of typical organic carbon-carbon and carbon-heteroatom bond formation reactions: aldol, Mannich, α-aminooxylation, and isogyric reactions. The ROET, which is an expansion of the reaction electronic theories (e.g., the frontier orbital theory) in terms of orbital energies, elucidates the reactive orbitals driving reactions and the charge transferability indices of the reactions. Performing the ROET analyses of these reactions shows that the charge transfer directions given in the rule-based diagrams of the electronic theory are reproduced even for the functional groups of charge transfer destinations in all but only two processes for 38 reaction processes. The ROET analyses also make clear the detailed orbital-based pictures of these bond formation reactions: that is, the use of the out-of-plane antibonding π orbitals in acidic conditions (enol-mode) and in-plane antibonding π orbitals in basic conditions (enolate-mode), which explain the experimentally assumed mechanisms such as the π-bond formations in acidic conditions and σ-bond formations at α-carbons in basic conditions. Furthermore, the ROET analyses explicate that the methyl group initially accepts electrons and then donates them to the bond formations in the target reactions. It is, consequently, suggested that the ROET serves a theoretical foundation for the electronic theory of organic chemistry.

Chemoton 2.0: Autonomous Exploration of Chemical Reaction Networks

Fueled by advances in hardware and algorithm design, large-scale automated explorations of chemical reaction space have become possible. Here, we present our approach to an open-source, extensible framework for explorations of chemical reaction mechanisms based on the first-principles of quantum mechanics. It is intended to facilitate reaction network explorations for diverse chemical problems with a wide range of goals such as mechanism elucidation, reaction path optimization, retrosynthetic path validation, reagent design, and microkinetic modeling. The stringent first-principles basis of all algorithms in our framework is key for the general applicability that avoids any restrictions to specific chemical systems. Such an agile framework requires multiple specialized software components of which we present three modules in this work. The key module, Chemoton, drives the exploration of reaction networks. For the exploration itself, we introduce two new algorithms for elementary-step searches that are based on Newton trajectories. The performance of these algorithms is assessed for a variety of reactions characterized by a broad chemical diversity in terms of bonding patterns and chemical elements. Chemoton successfully recovers the vast majority of these. We provide the resulting data, including large numbers of reactions that were not included in our reference set, to be used as a starting point for further explorations and for future reference.

Effect of alkali metal cations on network rearrangement in polyisoprene ionomers

The effects of cations, Li⁺, Na⁺, and Cs⁺, on the structure of ionic aggregates and network rearrangement in carboxylated polyisoprene (PI) ionomers were studied. We found that network rearrangement via interaggregate hopping of metal carboxylates is improved with a decrease in cation size, even though density functional theory (DFT) calculation indicated the increase in the attractive interaction between metal carboxylates. At the same time, we also found that as the size of the cation decreases, the inclusion of the PI segment in the ionic aggregate increases. The DFT calculation suggested the cation-π interaction between the cation and double bonds in the PI segment as the cause for the inclusion. The inclusion of the PI segment with a low glass transition temperature (T_g) plasticizes the ionic aggregate and would sterically hinder the attractive interaction between metal carboxylates. In fact, the electron spin resonance measurement revealed a decrease in the T_g of the ionic aggregate with a decrease in cation size. Based on our findings, we proposed that the inclusion of PI segments in the ionic aggregate is the possible cause for the enhancement of network rearrangement in the carboxylated PI ionomers with a decrease in the cation size.

Multistructural microiteration combined with QM/MM-ONIOM electrostatic embedding

Multistructural microiteration (MSM) is a method to take account of contributions of multiple surrounding structures in a geometrical optimization or reaction path calculation using the quantum mechanics/molecular mechanics (QM/MM) ONIOM method. In this study, we combined MSM with the electrostatic embedding (EE) scheme of the QM/MM-ONIOM method by extending its original formulation for mechanical embedding (ME). MSM-EE takes account of the polarization in the QM region induced by point charges assigned to atoms in the multiple surrounding structures, where the point charges are scaled by the weight factor of each surrounding structure determined through MSM. The performance of MSM-EE was compared with that of the other methods, i.e., ONIOM-ME, ONIOM-EE, and MSM-ME, by applying them to three chemical processes: (1) chorismate-to-prephenate transformation in aqueous solution, (2) the same transformation as (1) in an enzyme, and (3) hydroxylation in p-hydroxybenzoate hydroxylase. These numerical tests of MSM-EE yielded barriers and reaction energies close to experimental values with computational costs comparable to those of the other three methods.

Leveraging algorithmic search in quantum chemical reaction path finding

Reaction path finding methods construct a graph connecting reactants and products in a quantum chemical energy landscape. They are useful in elucidating various reactions and provide footsteps for designing new reactions. Their enormous computational cost, however, limits their application to relatively simple reactions. This paper engages in accelerating reaction path finding by introducing the principles of algorithmic search. A new method called RRT/SC-AFIR is devised by combining rapidly exploring random tree (RRT) and single component artificial force induced reaction (SC-AFIR). Using 96 cores, our method succeeded in constructing a reaction graph for Fritsch-Buttenberg-Wiechell rearrangement within a time limit of 3 days, while the conventional methods could not. Our results illustrate that the algorithm theory provides refreshing and beneficial viewpoints on quantum chemical methodologies.

Enhancement of the mechanical and thermal transport properties of carbon nanotube yarns by boundary structure modulation

Carbon nanotubes (CNTs) exhibit extremely high nanoscopic thermal/electrical transport and mechanical properties. However, the macroscopic properties of assembled CNTs are significantly lower than those of CNTs because of the boundary structure between the CNTs. Therefore, it is crucial to understand the relationship between the nanoscopic boundary structure in CNTs and the macroscopic properties of the assembled CNTs. Previous studies have shown that the nanoscopic phonon transport and macroscopic thermal transport in CNTs are improved by Joule annealing because of the improved boundary Van-der-Waals interactions between CNTs via the graphitization of amorphous carbon. In this study, we investigate the mechanical strength and thermal/electrical transport properties of CNT yarns with and without Joule annealing at various temperatures, analyzing the phenomena occurring at the boundaries of CNTs. The obtained experimental and theoretical results connect the nanoscopic boundary interaction of CNTs in CNT yarns and the macroscopic mechanical and transport properties of CNT yarns.

Fold-in Synthesis of a Pentabenzopentaaza[10]circulene

A pentabenzopentaaza[10]circulene has been synthesized as the largest fully conjugated hetero[n]circulene via a fold-in oxidative fusion reaction of an ortho-phenylene-bridged cyclic pyrrole pentamer. This circulene takes a saddle-distorted structure with bond lengths of the central ten-membered-ring in the range of 1.455-1.493 Å. Relatively broad absorption and fluorescence spectra were observed, which reflects its flexible structure, in accordance with the low-temperature NMR spectra and theoretical calculations. The energy barrier for saddle-to-saddle interconversion was estimated to be quite small (≈3 kcal mol^-1 ).

Electrochemical Dearomative Dicarboxylation of Heterocycles with Highly Negative Reduction Potentials

The dearomative dicarboxylation of stable heteroaromatics using CO₂ is highly challenging but represents a very powerful method for producing synthetically useful dicarboxylic acids, which can potentially be employed as intermediates of biologically active molecules such as natural products and drug leads. However, these types of transformations are still underdeveloped, and concise methodologies with high efficiency (e.g., high yield and high selectivity for dicarboxylations) have not been reported. We herein describe a new electrochemical protocol using the CO₂ radical anion (E_1/2 of CO₂ = -2.2 V in DMF and -2.3 V in CH₃CN vs SCE) that produces unprecedented trans-oriented 2,3-dicarboxylic acids from N-Ac-, Boc-, and Ph-protected indoles that exhibit highly negative reduction potentials (-2.50 to -2.94 V). On the basis of the calculated reduction potentials, N-protected indoles with reduction potentials up to -3 V smoothly undergo the desired dicarboxylation. Other heteroaromatics, including benzofuran, benzothiophene, electron-deficient furans, thiophenes, 1,3-diphenylisobenzofuran, and N-Boc-pyrazole, also exhibit reduction potentials more positive than -3 V and served as effective substrates for such dicarboxylations. The dicarboxylated products thus obtained can be derivatized into useful synthetic intermediates for biologically active compounds in few steps. We also show how the dearomative monocarboxylation can be achieved selectively by choice of the electrolyte, solvent, and protic additive; this strategy was then applied to the synthesis of an octahydroindole-2-carboxylic acid (Oic) derivative, which is a useful proline analogue.

New Tools for Taming Complex Reaction Networks: The Unimolecular Decomposition of Indole Revisited

The level of detail attained in the computational description of reaction mechanisms can be vastly improved through tools for automated chemical space exploration, particularly for systems of small to medium size. Under this approach, the unimolecular decomposition landscape for indole was explored through the automated reaction mechanism discovery program AutoMeKin. Nevertheless, the sheer complexity of the obtained mechanisms might be a hindrance regarding their chemical interpretation. In this spirit, the new Python library amk-tools has been designed to read and manipulate complex reaction networks, greatly simplifying their overall analysis. The package provides interactive dashboards featuring visualizations of the network, the three-dimensional (3D) molecular structures and vibrational normal modes of all chemical species, and the corresponding energy profiles for selected pathways. The combination of the joined mechanism generation and postprocessing workflow with the rich chemistry of indole decomposition enabled us to find new details of the reaction (obtained at the CCSD(T)/aug-cc-pVTZ//M06-2X/MG3S level of theory) that were not reported before: (i) 16 pathways leading to the formation of HCN and NH₃ (via amino radical); (ii) a barrierless reaction between methylene radical and phenyl isocyanide, which might be an operative mechanism under the conditions of the interstellar medium; and (iii) reaction channels leading to both hydrogen cyanide and hydrogen isocyanide, of potential astrochemical interest as the computed HNC/HCN ratios greatly exceed the calculated equilibrium value at very low temperatures. The reported reaction networks can be very valuable to supplement databases of kinetic data, which is of remarkable interest for pyrolysis and astrochemical studies.

Substituent and Solvent Effects on the Photoisomerization of Cinnamate Derivatives: An XMS-CASPT2 Study

Cinnamate derivatives show a variety of photo-induced reactions. Among them is trans-cis photoisomerization, which may involve the nonradiative decay (NRD) process. The extended multistate complete active space second-order perturbation (XMS-CASPT2) method was used in this study as a suitable theory for treating multireference electronic nature, which was frequently manifested in the photoisomerization process. The minimum energy paths of the trans-cis photoisomerization process of cinnamate derivatives were determined, and the activation energies were estimated using the resulting intrinsic reaction coordinate (IRC) paths. Natural orbital analysis revealed that the transition state's (TS) electronic structure is zwitterionic-like, elucidating the solvent and substituent effect on the energy barrier of photoisomerization paths. Furthermore, it was found that the charge on the pyramidalized carbon atom at the TS structure was strongly correlated with the activation energy barrier for the cinnamate derivatives. Thus, it seemingly provided a physical picture of the photoisomerization of cinnamates and was a good descriptor potentially applicable to molecular design for controlling the rate constant of the photoisomerization reaction.

Bringing out the potential of organoboron compounds by designing the chemical bonds and spaces around boron

Since the structures, reactivity and properties of organoboron compounds stem from the electron deficiency and low electronegativity of boron, the design of the chemical bonds attached to boron as well...

Radical Difunctionalization of Gaseous Ethylene Guided by Quantum Chemical Calculations: Selective Incorporation of Two Molecules of Ethylene

Ethylene, of which about 170 million tons are produced annually worldwide, is a fundamental C₂ feedstock that is widely used on an industrial scale for the synthesis of polyethylenes and polyvinylchlorides. Compared to other alkenes, however, the direct use of ethylene for the synthesis of fine chemicals such as pharmaceuticals and agrochemicals is limited, probably due to its small and gaseous character. We, herein, report a new radical difunctionalization strategy of ethylene, aided by quantum chemical calculations. Computationally proposed imidyl and sulfonyl radicals can be introduced into ethylene in the presence of an Ir photocatalyst under irradiation with blue light-emitting diodes (LEDs) (λ_max = 440 nm). The present reaction systems led to the selective incorporation of two molecules of ethylene into the substrate, which could be rationally explained by computational analysis.