Bridged Stilbenes: AIEgens Designed via a Simple Strategy to Control the Non-radiative Decay Pathway.

To broaden the application of aggregation-induced emission (AIE) luminogens (AIEgens), the design of novel small-molecular dyes that exhibit high fluorescence quantum yield (Φfl ) in the solid state is required. Considering that the mechanism of AIE can be rationalized based on steric avoidance of non-radiative decay pathways, a series of bridged stilbenes was designed, and their non-radiative decay pathways were investigated theoretically. Bridged stilbenes with short alkyl chains exhibited a strong fluorescence emission in solution and in the solid state, while bridged stilbenes with long alkyl chains exhibited AIE. Based on this theoretical prediction, we developed the bridged stilbenes BPST[7] and DPB[7], which demonstrate excellent AIE behavior.

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.

Computational study on the luminescence quantum yields of terbium complexes with 2,2'-bipyridine derivative ligands.

Terbium complexes are widely used as luminescent materials because of their bright green emission and sharp emission spectra and the independence of their emission wavelengths from the surrounding environment. The luminescence quantum yield (LQY), however, heavily depends on the surroundings, and an appropriate ligand design is indispensable. In this study, we focus on a Tb3+ complex coordinated by a 2,2'-bipyridine derivative ligand (L1), whose LQY is almost zero at room temperature [M. Hasegawa et al., New. J. Chem. 2014, 38, 1225] and compare it with a Tb3+ complex with a bipyridine ligand, which is widely used as a photo-antenna ligand. To discuss the LQYs of the complexes, we computed their energy profiles, i.e. the energetic and structural changes during the emission and quenching processes. The low LQY of the TbL1(NO3)2 complex was explained by the stability of the minimum energy crossing point between the potential energy surfaces of the ligand-centered lowest triplet state and the ground state, which was induced by the out-of-plane bending of the azomethine moiety. The most efficient way to improve the LQY by modification of the ligand is to replace the azomethine moieties by other functional groups, such as ether or reduced azomethine groups, whose minimum energy crossing points are unstable enough to reduce the rate of the quenching processes.

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.

Multistructural microiteration technique for geometry optimization and reaction path calculation in large systems.

We propose a multistructural microiteration (MSM) method for geometry optimization and reaction path calculation in large systems. MSM is a simple extension of the geometrical microiteration technique. In conventional microiteration, the structure of the non-reaction-center (surrounding) part is optimized by fixing atoms in the reaction-center part before displacements of the reaction-center atoms. In this method, the surrounding part is described as the weighted sum of multiple surrounding structures that are independently optimized. Then, geometric displacements of the reaction-center atoms are performed in the mean field generated by the weighted sum of the surrounding parts. MSM was combined with the QM/MM-ONIOM method and applied to chemical reactions in aqueous solution or enzyme. In all three cases, MSM gave lower reaction energy profiles than the QM/MM-ONIOM-microiteration method over the entire reaction paths with comparable computational costs. © 2017 Wiley Periodicals, Inc.

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.

The K-Region in Pyrenes as a Key Position to Activate Aggregation-Induced Emission: Effects of Introducing Highly Twisted N,N-Dimethylamines.

A new design strategy to activate aggregation-induced emission (AIE) in pyrene chromophores is reported. In a previous report, we demonstrated that highly twisted N,N-dialkylamines of anthracene and naphthalene induce drastic AIE when these donors are introduced at appropriate positions to stabilize the S1/S0 minimum energy conical intersection (MECI). In the present study, this design strategy was applied to pyrene: the introduction of N,N-dimethylamine substituents at the 4,5-positions of pyrene, the so-called K-region, are likely to stabilize MECIs. To examine this hypothesis, four novel pyrene derivatives, which contain highly twisted N,N-dimethylamino groups at the 4- (4-Py), 4,5- (4,5-Py), 1- (1-Py), or 1,6-positions (1,6-Py) were tested. The nonradiative transitions of 4,5-Py are highly efficient (knr = 57.1 × 107 s-1), so that its fluorescence quantum yield in acetonitrile decreases to Φfl = 0.04. The solid-state fluorescence of 4,5-Py is efficient (Φfl = 0.49). In contrast, 1,6-Py features strong fluorescence (Φfl = 0.48) with a slow nonradiative transition (knr = 11.0 × 107 s-1) that is subject to severe quenching (Φfl = 0.03) in the solid state. These results underline that the chemistry of the pyrene K-region is intriguing, both from a photophysical perspective and with respect to materials science.

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.

Trihydroborates and Dihydroboranes Bearing a Pentacoordinated Phosphorus Atom: Double Ring Expansion To Balance the Coordination States.

The first trihydroborate bearing a pentacoordinated phosphorus atom was synthesized as a new P-B bonded compound. Hydride abstraction of the trihydroborate gave an intermediary dihydroborane, which showed hydroboration reactivity and was trapped with pyridine whilst maintaining the P-B bond. The dihydroborane underwent a rearrangement, which involved a double ring expansion to compensate for the unbalanced coordination states of the phosphorus and boron atoms, to give a new fused bicyclic phosphine-boronate.

Highly Twisted N,N-Dialkylamines as a Design Strategy to Tune Simple Aromatic Hydrocarbons as Steric Environment-Sensitive Fluorophores.

The steric-environment sensitivity of fluorescence of 9,10-bis(N,N-dialkylamino)anthracenes (BDAAs) was studied experimentally and theoretically. A new design strategy to tune simple aromatic hydrocarbons as efficient aggregation-induced emission (AIE) luminogens and molecular rotors is proposed. For a variety of BDAAs, prominent Stokes shifts and efficient solid-state fluorescence were observed. Calculations on BDAA-methyl suggested that in the ground state (S0) conformations, the pyramidal amine groups are highly twisted, so that their lone-pair orbitals cannot conjugate with the anthracene π orbitals. Fluorescence takes place from the S1 minima, in which one or both amine groups are planarized. The stability of the S1 excited state minima as well as destabilization of the S0 state is the origin of large Stokes shift. Experimental measurement of the nonadiabatic transition rate suggests that para disubstitution by dialkylamino (or strongly electron-donating) groups is a key for fast internal conversion. Minimum energy conical intersection (MECI) between S1 and S0 states was found to have a Dewar-benzene like structure. Although this can be reached efficiently in liquid phase for fast internal conversion, a large amplitude motion is required to reach this MECI, which is prohibited in the solid state and caused efficient AIE. This strategy is used to find experimentally that naphthalene analogues are also efficient AIE luminogens. The flexibility of alkyl chains on amino groups is also found to be important for allowed charge-transfer transition. Thus, three points [(1) highly twisted N,N-dialkylamines, (2) substitution at the para positions, (3) with flexible alkyl groups] were proposed for activation of small aromatic hydrocarbons.

Trimeric cluster of lithium amidoborane-the smallest unit for the modeling of hydrogen release mechanism.

A detailed first-principle DFT M06/6-311++G(d.p) study of dehydrogenation mechanism of trimeric cluster of lithium amidoborane is presented. The first step of the reaction is association of two LiNH2 BH3 molecules in the cluster. The dominant feature of the subsequent reaction pathway is activation of H atom of BH3 group by three Li atoms with formation of unique Li3 H moiety. This Li3 H moiety is destroyed prior to dehydrogenation in favor of formation of a triangular Li2 H moiety, which interacts with protic H atom of NH2 group. As a result of this interaction, Li2 H2 moiety is produced. It features N(-) H(+) H(-) group suited near the middle plane between two Li(+) in the transition state that leads to H2 release. The transition states of association and hydrogen release steps are similar in energy. It is concluded that the trimer, (LiNH2 BH3 )3 , is the smallest cluster that captures the essence of the hydrogen release reaction.

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.

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.

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.

Theoretical insight into the wavelength-dependent photodissociation mechanism of nitric acid.

Photodissociation pathways of HNO3 involving the four lowest electronic singlet states (S0, S1, S2 and S3) were studied by the MS-CAS(12e,8o)PT2/6-31+G* method. All critical points, i.e. minima, transition states and minimum energy conical intersections, were explored systematically by the global reaction route mapping (GRRM) strategy utilizing the anharmonic downward distortion following (ADDF) method. Some key structures were also optimized at the MS-CAS(16e,12o)PT2/aug-cc-pVDZ level. Based on structures and relative energies of these critical points, we discussed the wavelength-dependent photodissociation mechanism of HNO3 in detail. The OH(X(2)Π) dissociation was found to occur through four pathways on the S1, S2, and S3 surfaces depending on the excitation energy, and one of the four pathways was shown to undergo the excited state roaming mechanism. The experimental appearance energy of the O((1)D) dissociation channel at 538.8 kJ mol(-1) was attributed to a conical intersection between the S2 and S3 surfaces. It has been shown experimentally that the O((1)D) channel was dominant at the higher excitation energy. In this study, its mechanistic reason was explained by the shape of the S3 surface and locations of important critical points.

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.

Different conical intersections control nonadiabatic photochemistry of fluorene light-driven molecular rotary motor: A CASSCF and spin-flip DFT study.

We report the light-driven isomerization mechanism of a fluorene-based light-driven rotary motor (corresponding to Feringa's 2nd generation rotary motor, [M. M. Pollard et al., Org. Biomol. Chem. 6, 507-512 (2008)]) at the complete active space self-consistent field (CASSCF) and spin-flip time-dependent density functional theory (TDDFT) (SFDFT) levels, combined with the complete active space second-order perturbation theory (CASPT2) single-point energy corrections. The good consistence between the SFDFT and CASSCF results confirms the capability of SFDFT in investigating the photoisomerization step of the light-driven molecular rotary motor, and proposes the CASPT2//SFDFT as a promising and effective approach in exploring photochemical processes. At the mechanistic aspect, for the fluorene-based motor, the S1/S0 minimum-energy conical intersection (MECIs) caused by pyramidalization of a fluorene carbon have relatively low energies and are easily accessible by the reactive molecule evolution along the rotary reaction path; therefore, the fluorene-type MECIs play the dominant role in nonadiabatic decay, as supported by previous experimental and theoretical works. Comparably, the other type of MECIs that results from pyramidalization of an indene carbon, which has been acting as the dominant nonadiabatic decay channel in the stilbene motor, is energetically inaccessible, thus the indene-type MECIs are "missing" in previous mechanistic studies including molecular dynamic simulations. A correlation between the geometric and electronic factors of MECIs and that of the S1 energy profile along the C═C rotary coordinate was found. The findings in current study are expected to deepen the understanding of nonadiabatic transition in the light-driven molecular rotary motor and provide insights into mechanistic tuning of their performance.

Nearly exclusive growth of small diameter semiconducting single-wall carbon nanotubes from organic chemistry synthetic end-cap molecules.

The inability to synthesize single-wall carbon nanotubes (SWCNTs) possessing uniform electronic properties and chirality represents the major impediment to their widespread applications. Recently, there is growing interest to explore and synthesize well-defined carbon nanostructures, including fullerenes, short nanotubes, and sidewalls of nanotubes, aiming for controlled synthesis of SWCNTs. One noticeable advantage of such processes is that no metal catalysts are used, and the produced nanotubes will be free of metal contamination. Many of these methods, however, suffer shortcomings of either low yield or poor controllability of nanotube uniformity. Here, we report a brand new approach to achieve high-efficiency metal-free growth of nearly pure SWCNT semiconductors, as supported by extensive spectroscopic characterization, electrical transport measurements, and density functional theory calculations. Our strategy combines bottom-up organic chemistry synthesis with vapor phase epitaxy elongation. We identify a strong correlation between the electronic properties of SWCNTs and their diameters in nanotube growth. This study not only provides material platforms for electronic applications of semiconducting SWCNTs but also contributes to fundamental understanding of the growth mechanism and controlled synthesis of SWCNTs.

From roaming atoms to hopping surfaces: mapping out global reaction routes in photochemistry.

The photodissociation of small molecules occurs upon irradiation by ultraviolet or visible light, and it is a very important chemical process in Earth's atmosphere, in the atmospheres of other planets, and in interstellar media. Photodissociation is an important method used to thoroughly investigate the fundamental issues of chemical reactivity. Photodissociation involves molecules and reaction fragments moving over ground- and excited-state potential surfaces (PESs). Molecules can move on a single PES (adiabatic pathway) or can cross over from one PES to another (nonadiabatic pathways). For a full theoretical understanding of a photodissociation mechanism, all of the important nonadiabatic and adiabatic pathways must be determined. This is not an easy task. We have developed an efficient computational method, called the global reaction route mapping (GRRM) strategy, that allows a theoretical exploration of ground- and excited-state PESs and their crossing seams in an automatic manner. In this Perspective, we summarize our approaches and present examples of their application together with newly determined chemical insights. These include the complex photodissociation mechanism of the formaldehyde molecule, the exclusive excited-state roaming dynamics of the nitrate radical, and all product channels and conformational memory in the photodissociation of the formic acid molecule. Finally, perspectives for the theoretical design of photofunctional molecules are discussed.

The Mechanism of Iron(II)-Catalyzed Asymmetric Mukaiyama Aldol Reaction in Aqueous Media: Density Functional Theory and Artificial Force-Induced Reaction Study.

Density functional theory (DFT), combined with the artificial force-induced reaction (AFIR) method, is used to establish the mechanism of the aqueous Mukaiyama aldol reactions catalyzed by a chiral Fe(II) complex. On the bases of the calculations, we identified several thermodynamically stable six- or seven-coordinate complexes in the solution, where the high-spin quintet state is the ground state. Among them, the active intermediates for the selectivity-determining outer-sphere carbon-carbon bond formation are proposed. The multicomponent artificial force-induced reaction (MC-AFIR) method found key transition states for the carbon-carbon bond formation, and explained the enantioselectivity and diastereoselectivity. The overall mechanism consists of the coordination of the aldehyde, carbon-carbon bond formation, the rate-determining proton transfer from water to aldehyde, and dissociation of trimethylsilyl group. The calculated full catalytic cycle is consistent with the experiments. This study provides important mechanistic insights for the transition metal catalyzed Mukaiyama aldol reaction in aqueous media.