Theoretical Study of the Addition of Cu–carbenes to Acetylenes to

5 days ago - Terminal alkynes have become one of the most versatile building blocks for C–C bond construction in the past few decades, and they are ...
0 downloads 0 Views 1009KB Size
Subscriber access provided by UNIV OF CAMBRIDGE

Article

Theoretical Study of the Addition of Cu– carbenes to Acetylenes to Form Chiral Allenes Kangbao Zhong, Chunhui Shan, Lei Zhu, Song Liu, Tao Zhang, Fenru Liu, Boming Shen, Yu Lan, and Ruopeng Bai J. Am. Chem. Soc., Just Accepted Manuscript • DOI: 10.1021/jacs.8b13055 • Publication Date (Web): 19 Mar 2019 Downloaded from http://pubs.acs.org on March 19, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 9 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Theoretical Study of the Addition of Cu–carbenes to Acetylenes to Form Chiral Allenes Kangbao Zhong, †# Chunhui Shan, †# Lei Zhu, † Song Liu, † Tao Zhang, † Fenru Liu, † Boming Shen, † Yu Lan,* †‡ Ruopeng Bai* † †

School of Chemistry and Chemical Engineering, Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 400030, China ‡

College of Chemistry and Molecular Engineering, Zhengzhou University, Zhengzhou 450001, China

KEYWORDS: acetylene, Cu–carbene, DFT calculations, direct nucleophilic addition. ABSTRACT: Terminal alkynes have become one of the most versatile building blocks for C–C bond construction in the past few decades, and they are usually considered to convert to acetylides before further transformations. In this study, a novel direct nucleophilic addition mode for Cu(I)-catalyzed crosscoupling of terminal alkynes and N-tosylhydrazones to synthesize chiral allenes is proposed, and it was investigated by density functional theory with the M11-L density functional. Three different reaction pathways were considered and investigated. The computational results show that the proposed reaction pathway, which includes direct electrophilic attack of protonated acetylene, deprotonation of the vinyl cation, and catalyst regeneration, is the most favorable pathway. Another possible deprotonation–carbenation–insertion pathway is shown to be unfavorable. The direct electrophilic addition step is the rate- and enantioselectivity-determining step in the catalytic cycle. Noncovalent interaction analysis shows that the steric effect between the methyl group of the carbene moiety and the naphthalyl group of the bisoxazoline ligand is important to control the enantioselectivity. In addition, calculation of a series of chiral bisoxazoline ligands shows that a bulky group on the oxazoline ring is favorable for high enantioselectivity, which agrees with experimental observations. Moreover, copper acetylides are stable and their generation is a favorable pathway in the absence of chiral bisoxazoline ligands.

INTRODUCTION The carbon–carbon triple bond is one of the most basic functional groups, and it widely exists in pharmaceuticals, agrochemicals, and functional materials.1 Among these compounds, terminal alkynes have attracted significant attention because of their unique structural features and reactivities,2-7 and they have been widely investigated in the fields of organic synthesis, natural product chemistry, and materials science.1 The reaction modes of terminal alkynes are summarized in Figure 1. Because of the acidity of the terminal hydrogen atom of terminal alkynes, the C–H bond can be deprotonated to generate the corresponding acetylide, which usually acts as a nucleophile to react with possible electrophiles (type a). 810 The π bonds of alkynes are electron-rich, so electrophilic addition of electrophiles can occur to give the corresponding C–C double bonds (type b).11,12 The C–C triple bond can act as a ligand and coordinate to transition metal centers with moderate π-backbonding, which weakens the C≡C bond and enables subsequent nucleophilic attack to generate vinyl metal intermediates (type c). Moreover, coordination to a transition metal also increases the acidity of the terminal hydrogen atom of the terminal alkyne group, so subsequent deprotonation can also occur to give a metal–acetylide species (type d).13-26 In the past few decades, a series of new synthetic methodologies have enhanced alkyne chemistry and made alkynes one of the most versatile building blocks.27-35 Among these methodologies, transition-metal-catalyzed coupling reactions, such as Glaser–Hay reactions,36,37 Sonogashira reactions,38,39 Glaser–Eglinton reactions,40,41 Cadiot–Chodkiewicz reactions,42 and Castro–Stephens reactions,43,44 have received extensive attention as simple and convenient methods to construct internal C–C triple and double bonds.

Terminal alkynes usually undergo a deprotonation process to generate active acetylides before further conversion, such as transmetalation or nucleophilic addition.45-53 C–C bond formation without C–H bond cleavage of the terminal alkyne has been shown to be an unfavorable pathway. However, further mechanistic validation is required, especially under complex reaction conditions. For instance, in a previous study, we showed that direct nucleophilic addition of phenylacetylene to an imine occurs in the presence of a Ni catalyst, and deprotonation to generate nickel acetylide is unfavorable.23 These results show an

Figure 1. Activation modes of terminal alkynes.

alternative mode for terminal alkyne conversion could be possible in transition metal catalysis. Recently, Wang and co-workers reported a straightforward method to synthesize substituted allenes by Cu(I)-catalyzed cross-coupling of terminal alkynes and N-tosylhydrazones (Scheme 1a).54 Subsequently, they developed a Cu(I)-catalyzed cross-coupling reaction between carbenoids and aryl acetylenes (Scheme 1b).55 The excellent functional group compatibility, high yield, and good enantioselectivity make this reaction an efficient method for enantioselective synthesis of trisubstituted allenes.

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

COMPUTATIONAL METHODS All of the density functional theory (DFT) calculations were performed with the Gaussian 09 series of programs.76 The B3-LYP functional77,78 with the standard 6-31G(d) basis set (SDD basis set for Cu) was used for the geometry optimizations in the gas phase. Harmonic vibrational frequency calculations were performed for all of the stationary points to determine whether they are local minima or transition structures and to derive the thermochemical corrections for the enthalpies and free energies. The M11-L functional79-82 proposed by Truhlar et al. with the 6-311+G(d,p) basis set (SDD basis set for Cu) was used to calculate the single-point energies in chloroform solvent to provide more accurate energy information.83-86 The solvent effect was considered by single-point calculations based on the gas-phase stationary points with the SMD87,88 continuum solvation model. The Gibbs free energies of the stationary points calculated using the M11-L functional are used to discuss the energies. The noncovalent interactions (NCI)89,90, lowest unoccupied molecular orbitals (LUMO), and electrostatic potentials Scheme 2. Possible Reaction Pathways for the Cu(I)-Catalyzed Cross-Coupling Reaction Between the Carbenoid and Acetylene. Figure 2. Possible mechanism for Cu(I)-catalyzed cross-coupling between a diazo group and a terminal alkyne.

In copper catalysis, the copper species generally plays two roles: terminal acetylene activation to generate a copper acetylide complex56-65 and carbenoid activation to give a metal–carbene species66-75 (Figure 2). Wang et al reported that the copper acetylide intermediate can react with a carbenoid to give a copper–carbene acetylide complex. They also speculated that the enantioselectivity is controlled by the intramolecular carbene insertion step. However, in our previous theoretical studies, we found that copper-assisted deprotonation of acetylene is a balanced process, but carbenation of copper is irreversible. Therefore, we infer that carbenation would occur before alkynylation. Indeed, when a copper–carbene species forms, the steric effect of a bulky ligand is unfavorable for further alkynylation. Electrophilic attack of the carbene moiety to acetylene also needs to be considered. Consequently, we became interested in this copper-catalyzed enantioselective trisubstituted allene synthesis reaction. In particular, we wanted to clarify the detailed alkyne conversion mode, that is, whether the copper acetylide is the key intermediate or the terminal alkyne directly adds to the copper–carbene moiety.

Scheme 1. Cu(I)-Catalyzed Cross-Coupling Reaction Between a Diazo group and a Terminal Alkyne.

(ESP) were calculated at the B3LYP/6-31G(d) (SDD for Cu) level.

RESULTS AND DISCUSSION Considering the rule of copper species, we propose three putative pathways for Cu-catalyzed cross-coupling of terminal alkynes and N-tosylhydrazones in Scheme 2. All of the pathways start with ligand-coordinated Cu(I) cationic species I. In the pathway proposed by Wang et al. (path A), after amine-assisted deprotonation of phenylacetylene, copper acetylide intermediate III is generated. Carbenation then occurs to give copper– carbene acetylide species V. An alternative pathway to generate species V involves carbenation of catalyst I followed by deprotonation of phenylacetylene (path B). Migratory insertion of the acetylide group into the carbene moiety then gives propargyl copper intermediate VI. Finally, chiral allene VII is generated by protonation of VI. We propose another path-

ACS Paragon Plus Environment

Page 2 of 9

Page 3 of 9 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society way (path C) where once copper–carbene species VIII is generated, because of the electrophilicity of the carbene moiety, electrophilic addition of the terminal alkyne into the carbene π bond occurs to give vinyl cation

Figure 3. Cu(I)-catalyzed cross-coupling reaction between 1-(p-chlorophenyl)-1-diazoethane and phenylacetylene.

IX, which can be deprotonated by the amine to directly generate the common propargyl copper intermediate VI. With the above proposal in hand, theoretical calculations were performed to reveal the mechanism of the Cu(I)-catalyzed cross-coupling reaction between the carbenoid and acetylene. We initially focused on Wang’s proposed pathway (path A). As in their experimental work, we chose coupling of 1-(p-chlorophenyl)-1-diazoethane and phenylacetylene as the model reaction, in which the bisoxazoline ligand is selected as the chiral ligand. The amine-coordinated cationic copper(I) complex was chosen as the starting point for the free energy profiles (Figure 4). Coordination of phenylacetylene to the copper center generates intermediate 2, which is exergonic by 7.2 kcal/mol. Deprotonation of phenylacetylene in intermediate 2 by triethylamine then generates copper(I) acetylide 5 via transition state

Figure 4. Free energy profile of path A for the Cu(I)-catalyzed cross-coupling reaction between 1-(p-chlorophenyl)-1-diazoethane and phenylacetylene. The energies are in kcal/mol and represent the relative free energies calculated with the DFT/M11-L method in chloroform. The bond distances are in angstroms. Noncovalent interactions analysis of 7-ts (blue, attraction; green, weak interaction; red: steric effect).

4-ts, which is endergonic by 14.6 kcal/mol. The energy barrier for this step is 16.7 kcal/mol. Subsequent carbenation with coordination of carbenoid 6 occurs via transition state 7-ts by dissociation of nitrogen gas with an energy barrier of 26.2 kcal/mol to give copper–carbene acetylide complex 8. In the generated copper–carbene 8, the Cu–acetylide and Cu–carbene bond lengths are 1.93 and 1.88 Å, respectively. Next, migration insertion of the acetylide moiety into the Cu–carbene bond via three-membered ring transition state 9-ts releases 39.8 kcal/mol free energy, and propargyl copper complex 10 is generated. The C-Cu single bond length in 10 is 2.01 Å.

The energy barrier for the migration insertion step is 3.8 kcal/mol, indicating fast conversion of 8. Finally, protonation of 10 releases the final product 12 via transition state 11-ts, which is exergonic by 23.6 kcal/mol, and active catalyst copper species 1 is regenerated. The activation free energy for path A is 40.8 kcal/mol, indicating that it is an unfavorable pathway. The calculated rate-limiting step of path A is denitrogenation to generate copper–carbene acetylide complex 8. NCI analysis was performed to investigate the reactivity of this step. The calculated NCI map of transition state 7-ts is shown in Figure 4. It reveals that the interaction between the

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

benzyl group of the bisoxazoline ligand and the phenyl group of the reacting diazo compound (highlighted by the red circle) leads to significant steric repulsion, which is the reason for the high activation energy. The calculated high energy barrier of carbenation after alkynylation suggests that if the reaction order is changed, the extra steric repulsion will be released. Therefore, another pathway (path B) was also considered. The free energy profile of path B is shown in Figure 5. Ligand exchange with diazo 6 generates complex 13, which is endergonic by 7.6 kcal/mol. Subsequent denitrogenation of intermediate 13 occurs via transition state 14-ts by cleavage of the C–N bond with an energy barrier of only 5.2 kcal/mol to give copper–carbene intermediate 15, which releases 15.0 kcal/mol free energy. Phenylacetylene coordinated to intermediate 15 is deprotonated by the amine via transition state 16-ts with generation of the corresponding copper–carbene acetylide 8, which is endergonic by 29.3 kcal/mol. The energy barrier of deprotonation is 35.1 kcal/mol. Subsequent carbene insertion into the copper–acetylide bond and protonation are the same as path A, and they are not the rate-determining step. The total activation energy of path B is 35.1 kcal/mol. Because of the high energy barrier of the alkynylation step in path B, an alternative pathway is taken into account (path C), in which the terminal

Page 4 of 9

alkyne is directly attacked by the copper–carbene species (Figure 5). The energy barrier of direct intermolecular electrophilic addition of copper– carbene to phenylacetylene via transition state 17-ts is only 26.9 kcal/mol, which is 8.2 kcal/mol lower than that of the previously mentioned deprotonation step. In transition state 17-ts, the carbene carbon atom acts as electrophile to attack the C–C triple bond of phenylacetylene Vinyl cation species 18 is generated, in which the C–Cu bond length is 2.10 Å, and the length of the newly generated C–C single bond is 1.42 Å. The distances of Cu-C2 and Cu-C3 is 2.28 and 3.24 Å, and the bond angle of Cu-C1-C2 is 78.1o, which is indicates that there is no strong bonding between C-C double bond and copper center. To further validate our propose, we calculated the molecular orbitals of 18. The computational result shows that there is a strong orbital overlap between Cu and C1, and no obvious bonding was found between Cu and C2/C3 in occupied orbitals. Therefore, there is no coordination between copper and C-C double bond in 18. Direct nucleophilic addition is endergonic by only 7.1 kcal/mol. Vinyl cation 18 is deprotonated by the amine via transition state 19-ts, generating propargyl copper complex 10.

Figure 5. Free energy profiles of paths B and C. The energy values are in kcal/mol and represent the relative free energies calculated by the DFT/M11-L method in chloroform. The bond distances are in angstroms.

ACS Paragon Plus Environment

Page 5 of 9 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure 6. Optimized structures of transition states 16-ts, 17-ts, 17′-ts, and 19-ts, intermediates 18 and 18′, and the molecular orbitals of intermediate 18. The bond lengths are in angstroms. AC2–C3–C4 is the C2–C3–C4 bond angle and ACu–C1–C2 is the Cu–C1–C2 bond angle.

The deprotonation step is exergonic by 17.6 kcal/mol and the energy barrier for this step is only 11.3 kcal/mol. Finally, similar to paths A and B, allene product 12 is released by protonation of 10 and active catalytic species 1 is regenerated by coordination of amine. The rate-determining step in path C is nucleophilic addition of acetylene and the overall activation energy is 26.9 kcal/mol, which agrees with experimental observation. To further investigate the possibility of the nucleophilic addition step with the carbene species, the electrostatic potential (ESP) surface was calculated and molecular orbital (MO) analysis was performed for copper– carbene species 15 with the M11-L density functional. The ESP surface shows that the positive charge in the carbene moiety is mainly located on the carbene carbon atom (Figure 7a). Moreover, the LUMO of 15 mainly consists of the p orbital of the carbene carbon atom (Figure 7b), which clearly reveal typical Fischer-type metal–carbene character. These results indicate that the carbene carbon atom is electron deficient and can be easily attacked by weak nucleophiles, such as terminal alkynes.

this step. The (R)-enantiomer 12 is generated by Si-face attack of the carbene via transition state 17-ts. Generation of the corresponding (S)-enantiomer 12′ is accomplished by the corresponding Re-face attack. The energy barrier of (S)-enantiomer generation via 17′-ts is 28.7 kcal/mol, which is 1.8 kcal/mol higher than that of 17-ts, indicating that generation of 12′ is unfavorable. Therefore, (R)-12 is considered to be the main product. The calculated enantiomeric excess (ee) is 91%, which is in good agreement with the experimental value (88% ee).

Figure 8. NCI analysis of 17-ts and 17′-ts (blue, attraction; green, weak interaction; red: steric effect).

Figure 7. (a) Electrostatic potential map of carbenoid copper intermediate 15. (b) LUMO energy of carbenoid copper intermediate 15 and the ratio of C1.

In the theoretical calculations, we found that nucleophilic addition is the rate-determining step and the enantioselectivity is also controlled in

To investigate the origin of the enantioselectivity, NCI analysis of the nucleophilic addition transition states 17-ts and 17′-ts was performed (Figure 8). There is a clear steric effect between the methyl group of the carbene moiety and the naphthalyl group of the bisoxazoline ligand in transition state 17′-ts, which reveals repulsion in this case. Furthermore, there is weak π-stacking between the phenyl group of the reacting acetylene and the naphthalyl group of transition state 17-ts. Therefore, the enantioselectivity is mainly controlled by the steric effect of the electrophilic addition transition state, which means that chiral bisoxazoline with large substituent groups is important for the enantioselectivity. From the above results, the substituent effect of the chiral bisoxazoline ligand is important for the enantioselectivity of copper-catalyzed crosscoupling of diazo compounds and terminal

Table 1. Theoretically Predicated and Experimentally Observed Enantioselectivities for a Series of Chiral Bisoxazoline Ligands.

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

R1

R2

ΔΔG‡ (kcal/mol)

1

Bn

2-Nap

1.8

96 : 4

94 : 6

2

Bn

Ph

0.7

77 : 23

87 :13

3

Me

Ph

0.1

55 : 45

67 : 33

Entry

er Calc.

Exp.

alkynes. Therefore, further theoretical calculations were performed to investigate how the chiral ligand affects the enantioselectivity of the allene. When the R group of bisoxazoline is changed from the 2-naphthalyl group to the phenyl group, the energy gap between Re-face and Si-face attack decreases to 0.7 kcal/mol, which is 1.1 kcal/mol lower than that with the naphthalyl group (Table 1). The corresponding calculated enantiomeric ratio (er) is 77:23, indicating an obvious enantioselectivity decrease. When the R group is changed from the benzyl to methyl (Table 1, entry 3), the energy gap between Re-face and Si-face attack decreases to only 0.1 kcal/mol, indicating poor enantioselectivity. In fact, the calculated er based on the energy profile is 55:45, which reasonably well agrees with the experimental result (67:33). These computational results are in agreement with our previous suggestion that the enantioselectivity is mainly controlled by the steric effect of the electrophilic addition transition state. Based on the calculated free energy profiles of the possible reaction pathways and the prediction of the origin of enantioselectivity, we can conclude that the catalytic cycle of copper-catalyzed cross-coupling of diazo compounds and terminal alkynes includes carbenation of the copper(I) catalyst, electrophilic attack of the terminal alkyne by the carbene moiety, and proton transfer. Both the rate- and enantioselectivity-determining steps are considered to be the electrophilic attack step. The theoretical calculations predict that a bulky ligand is unfavorable for alkynylation after carbenation of copper. Moreover, the Fischer-type cationic copper–carbene complex shows significant electrophilicity, so subsequent electrophilic attack of the terminal alkyne can easily occur. Following this idea, we propose that if a small ligand or ligand-free catalytic system is applied with a cuprate, an alternative pathway would occur. Therefore, we also investigated the mechanism of another copper-catalyzed cross-coupling of acetylene and N-tosylhydrazones pathway for achiral synthesis of allenes proposed by Wang et al. Differing from the above model reaction, a bidentate ligand-free copper-catalysis system, in which LiOtBu was used as base and ligand, and gives racemic allenes, was also reported experimentally. To further validate the effects of LiOtBu, the free energy profiles of paths A–C for Cu(I)-catalyzed cross-coupling between N-tosylhydrazones and phenylacetylene are also systematically taken into account (for details, see Supporting Information (SI)). The computational result shows that the deprotonation with LiOtBu is drastic exothermic, and the generated copper(I) acetylene intermediate is thermodynamic stable, which is different from the above results. Therefore, the deprotonation–carbenation–insertion pathway is the most favorable under strongly basic conditions.

CONCLUSIONS A novel electrophilic addition–deprotonation pathway for cationic Cu(I)-catalyzed cross-coupling of terminal alkynes and N-tosylhydrazones to synthesize chiral allenes has been investigated by DFT with the M11-L density functional. This novel pathway consists of four main steps:

(1) carbenation of the chiral bisoxazoline-coordinated copper species to generate a copper–carbene intermediate, (2) direct electrophilic addition of the carbene carbon moiety to the terminal alkyne to generate a vinyl cation species, (3) deprotonation of the vinyl cation to form a propargyl copper intermediate, and (4) and protonation to release the chiral substituted allene with regeneration of the copper catalyst. Differing from a previously reported mechanism, copper acetylide is not present in this new pathway. The calculated total activation energy is 26.9 kcal/mol. The direct electrophilic addition step is the rate-determining step in the catalytic cycle. DFT calculations reveal that the carbenation/migration insertion pathway proposed by Wang et al. is unfavorable because the overall activation energy is 40.8 kcal/mol. In addition, the enantioselectivity-determining step is predicted to be electrophilic attack of copper–carbene to the terminal alkyne. NCI analysis shows that the steric effect between the methyl group of the carbene moiety and the naphthalyl group of the bisoxazoline ligand is important to control the enantioselectivity. A series of chiral bisoxazoline ligands were considered. We found that a bulky group on the oxazoline is favorable for high enantioselectivity, which agrees with experimental observations. The theoretical calculations predict that a bulky ligand is unfavorable for alkynylation after carbenation of copper. Moreover, the Fischer-type cationic copper–carbene complex shows significant electrophilicity, so the subsequent electrophilic attack of acetylene can easily occur. Following this idea, we propose that if a small ligand or ligand-free catalytic system is applied with a cuprate, an alternative pathway will occur. A contrast study of the ligand-free cuprate-catalyzed cross-coupling reaction of carbenoids and terminal alkynes was also performed. The computational results show that the regular pathway involving alkynylation, carbenation, carbene insertion, and protonation is favorable in this case. We believe that this study will aid in the understanding of copper-catalyzed cross-coupling of carbenoids and terminal alkynes, and it will provide a practical theoretical guide for further experimental investigations.

ASSOCIATED CONTENT Supporting Information This material is available free of charge via the Internet at http://pubs.acs.org. Computational details, Cartesian coordinates and energies of all reported structures.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] *E-mail: [email protected] ORCID Yu Lan: 0000-0002-2328-0020 Ruopeng Bai: 0000-0002-1097-8526

Author Contributions #

These authors contributed equally.

Notes The authors declare no competing financial interests.

ACKNOWLEDGMENT This project was supported by the National Science Foundation of China (Grants 21822303 and 21772020), Fundamental Research Funds for the Central Universities (Chongqing University) (No. 2018CDXZ0002; 2018CDPTCG0001/4; 2018CDYJSY0055) and Chongqing Science and Technology Committee. We are also grateful for support by Chongqing Postdoctoral Science Special Foundation (XmT2018085) and graduate research and innovation foundation of Chongqing, China (Grant No. CYB18043).

ACS Paragon Plus Environment

Page 6 of 9

Page 7 of 9 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

REFERENCES 1 Acetylene Chemistry: Chemistry, Biology and Material Science; Diederich, F., Stang, P. J., Tykwinski, R. R., Eds.; Wiley-VCH: Weinheim, Germany, 2005. 2 Chinchilla, R.; Najera, C., Chemicals from alkynes with palladium catalysts. Chem. Rev. 2014, 114, 1783-1826. 3 Shi, W.; Lei, A., 1,3-Diyne chemistry: synthesis and derivations. Tetrahedron Lett. 2014, 55, 2763-2772. 4 Liu, C.; Yuan, J.; Gao, M.; Tang, S.; Li, W.; Shi, R.; Lei, A. Oxidative Coupling between Two Hydrocarbons: An Update of Recent C-H Functionalizations. Chem. Rev. 2015, 115, 12138-12204. 5 Boyarskiy, V. P.; Ryabukhin, D. S.; Bokach, N. A.; Vasilyev, A. V., Alkenylation of Arenes and Heteroarenes with Alkynes. Chem. Rev. 2016, 116, 58945986. 6 Trost, B. M.; Masters, J. T., Transition metal-catalyzed couplings of alkynes to 1,3-enynes: modern methods and synthetic applications. Chem. Soc. Rev. 2016, 45, 2212-2238. 7 de Orbe, M. E.; Amenos, L.; Kirillova, M. S.; Wang, Y.; Lopez-Carrillo, V.; Maseras, F.; Echavarren, A. M., Cyclobutene vs 1,3-Diene Formation in the GoldCatalyzed Reaction of Alkynes with Alkenes: The Complete Mechanistic Picture. J. Am. Chem. Soc. 2017, 139, 10302-10311. 8 Dureen, M. A.; Stephan, D. W., Terminal alkyne activation by frustrated and classical Lewis acid/phosphine pairs. J. Am. Chem. Soc. 2009, 131, 8396-8397. 9 Li, C. J., The development of catalytic nucleophilic additions of terminal alkynes in water. Acc. Chem. Res. 2010, 43, 581-590. 10 Kim, S. H.; Yoon, J.; Chang, S., Palladium-catalyzed oxidative alkynylation of heterocycles with terminal alkynes under air conditions. Org. Lett. 2011, 13, 14741477. 11 Ananikov, V. P.; Gayduk, K. A.; Beletskaya, I. P.; Khrustalev, V. N.; Antipin, M. Y., Remarkable ligand effect in Ni- and Pd-catalyzed bisthiolation and bisselenation of terminal alkynes: solving the problem of stereoselective dialkyldichalcogenide addition to the C triple chemical bond C Bond. Chem.-Eur. J. 2008, 14, 2420-2434. 12 Ohmura, T.; Oshima, K.; Taniguchi, H.; Suginome, M., Switch of regioselectivity in palladium-catalyzed silaboration of terminal alkynes by liganddependent control of reductive elimination. J. Am. Chem. Soc. 2010, 132, 1219412196. 13 Le Lagadec, R.; Roman, E.; Toupet, L.; Mueller, U.; Dixneuf, P. H., (C5Me5)Ru-vinylidene complexes from terminal alkynes and propargyl alcohol derivatives. Organometallics 1994, 13, 5030-5039. 14 Alonso, F.; Beletskaya, I. P.; Yus, M., Transition-metal-catalyzed addition of heteroatom-hydrogen bonds to alkynes. Chem. Rev. 2004, 104, 3079-3159. 15 Rodionov, V. O.; Presolski, S. I.; Diaz, D. D.; Fokin, V. V.; Finn, M. G., Ligand-accelerated Cu-catalyzed azide-alkyne cycloaddition: a mechanistic report. J. Am. Chem. Soc. 2007, 129, 12705-12712. 16 Shen, R.; Chen, T.; Zhao, Y.; Qiu, R.; Zhou, Y.; Yin, S.; Wang, X.; Goto, M.; Han, L. B., Facile regio- and stereoselective hydrometalation of alkynes with a combination of carboxylic acids and group 10 transition metal complexes: selective hydrogenation of alkynes with formic acid. J. Am. Chem. Soc. 2011, 133, 17037-17044. 17 Han, J.; Paton, R. S.; Xu, B.; Hammond, G. B., Synthesis of Cyclic alphaAminophosphonates through Copper Catalyzed Enamine Activation. Synthesis 2013, 45, 463-470. 18 Bai, R.; Zhang, G.; Yi, H.; Huang, Z.; Qi, X.; Liu, C.; Miller, J. T.; Kropf, A. J.; Bunel, E. E.; Lan, Y.; Lei, A., Cu(II)-Cu(I) synergistic cooperation to lead the alkyne CH activation. J. Am. Chem. Soc. 2014, 136, 16760-16763. 19 Qi, X.; Li, Y.; Zhang, G.; Li, Y.; Lei, A.; Liu, C.; Lan, Y., Dinuclear versus mononuclear pathways in zinc mediated nucleophilic addition: a combined experimental and DFT study. Dalton Trans. 2015, 44, 11165-11171. 20 Qi, X.; Zhang, H.; Shao, A.; Zhu, L.; Xu, T.; Gao, M.; Liu, C.; Lan, Y., Silver Migration Facilitates Isocyanide-Alkyne [3 + 2] Cycloaddition Reactions: Combined Experimental and Theoretical Study. ACS Catal. 2015, 5, 6640-6647. 21 Qi, X.; Bai, R.; Zhu, L.; Jin, R.; Lei, A.; Lan, Y., Mechanism of Synergistic Cu(II)/Cu(I)-Mediated Alkyne Coupling: Dinuclear 1,2-Reductive Elimination after Minimum Energy Crossing Point. J. Org. Chem. 2016, 81, 1654-1660. 22 Straker, R. N.; Peng, Q.; Mekareeya, A.; Paton, R. S.; Anderson, E. A., Computational ligand design in enantio- and diastereoselective ynamide [5+2] cycloisomerization. Nat. Commun. 2016, 7, 10109. 23 Liu, R. R.; Zhu, L.; Hu, J. P.; Lu, C. J.; Gao, J. R.; Lan, Y.; Jia, Y. X., Enantioselective alkynylation of N-sulfonyl alpha-ketiminoesters via a Friedel-Crafts alkylation strategy. Chem. Commun. 2017, 53, 5890-5893.

24 Qi, X.; Li, Y.; Bai, R.; Lan, Y., Mechanism of Rhodium-Catalyzed C-H Functionalization: Advances in Theoretical Investigation. Acc. Chem. Res. 2017, 50, 2799-2808. 25 Yue, X.; Qi, X.; Bai, R.; Lei, A.; Lan, Y., Mononuclear or Dinuclear? Mechanistic Study of the Zinc-Catalyzed Oxidative Coupling of Aldehydes and Acetylenes. Chem.-Eur. J. 2017, 23, 6419-6425. 26 Shan, C.; Zhu, L.; Qu, L. B.; Bai, R.; Lan, Y., Mechanistic view of Rucatalyzed C-H bond activation and functionalization: computational advances. Chem. Soc. Rev. 2018, 47, 7552-7576. 27 Thorand, S.; Krause, N., Improved procedures for the palladium-catalyzed coupling of terminal Alkynes with aryl bromides (Sonogashira coupling). J. Org. Chem. 1998, 63, 8551-8553. 28 Langille, N. F.; Dakin, L. A.; Panek, J. S., Sonogashira Coupling of Functionalized Trifloyl Oxazoles and Thiazoles with Terminal Alkynes:  Synthesis of Disubstituted Heterocycles. Org. Lett. 2002, 4, 2485-2488. 29 Wei, C.; Li, C.-J., Enantioselective Direct-Addition of Terminal Alkynes to Imines Catalyzed by Copper(I)pybox Complex in Water and in Toluene. J. Am. Chem. Soc. 2002, 124, 5638-5639. 30 Shintani, R.; Fu, G. C., A new copper-catalyzed [3 + 2] cycloaddition: enantioselective coupling of terminal alkynes with azomethine imines to generate fivemembered nitrogen heterocycles. J. Am. Chem. Soc. 2003, 125, 10778-10779. 31 Wei, C.; Mague, J. T.; Li, C. J., Cu(I)-catalyzed direct addition and asymmetric addition of terminal alkynes to imines. Proc. Natl. Acad. Sci. U. S. A. 2004, 101, 5749-5754. 32 Shang, M.; Wang, H. L.; Sun, S. Z.; Dai, H. X.; Yu, J. Q., Cu(II)-mediated ortho C-H alkynylation of (hetero)arenes with terminal alkynes. J. Am. Chem. Soc. 2014, 136, 11590-11593. 33 Tang, S.; Liu, K.; Long, Y.; Qi, X.; Lan, Y.; Lei, A., Tuning radical reactivity using iodine in oxidative C(sp(3))-H/C(sp)-H cross-coupling: an easy way toward the synthesis of furans and indolizines. Chem. Commun. 2015, 51, 8769-8772. 34 Tang, S.; Wang, P.; Li, H.; Lei, A., Multimetallic catalysed radical oxidative C(sp(3))-H/C(sp)-H cross-coupling between unactivated alkanes and terminal alkynes. Nat. commun. 2016, 7, 11676. 35 Till, N. A.; Smith, R. T.; MacMillan, D. W. C., Decarboxylative Hydroalkylation of Alkynes. J. Am. Chem. Soc. 2018, 140, 5701-5705. 36 Hay, A., Communications- Oxidative Coupling of Acetylenes. J. Org. Chem. 1960, 25, 1275-1276. 37 Hay, A. S., Oxidative Coupling of Acetylenes. II1. J. Org. Chem. 1962, 27, 3320-3321. 38 Sonogashira, K.; Tohda, Y.; Hagihara, N., A convenient synthesis of acetylenes: catalytic substitutions of acetylenic hydrogen with bromoalkenes, iodoarenes and bromopyridines. Tetrahedron Lett. 1975, 16, 4467-4470. 39 Negishi, E.; Anastasia, L., Palladium-catalyzed alkynylation. Chem. Rev. 2003, 103, 1979-2017. 40 Glaser, C., Beiträge zur Kenntniss des Acetenylbenzols. Ber. Dtsch. Chem. Ges. 1869, 2, 422-424. 41 Glaser, C., Untersuchungen über einige Derivate der Zimmtsäure. Annalen der Chemie und Pharmacie 1870, 154, 137-171. 42 Chodkiewicz, W., Synthesis of acetylenic compounds. Ann. Chim. Paris. 1957, 2, 819. 43 Castro, C. E.; Stephens, R. D., Substitutions by ligands of low valent transition metals. A preparation of tolans and heterocyclics from aryl iodides and cuprous acetylides. J. Org. Chem. 1963, 28, 2163. 44 Stephens, R. D.; Castro, C. E. The Substitution of Aryl Iodides with Cuprous Acetylides. A Synthesis of Tolanes and Heterocyclics1. J. Org. Chem. 1963, 28, 3313-3315. 45 Trost, B. M.; Sorum, M. T.; Chan, C.; Rühter, G., Palladium-Catalyzed Additions of Terminal Alkynes to Acceptor Alkynes. J. Am. Chem. Soc. 1997, 119, 698708. 46 Chen, Z.; Li, J.; Jiang, H.; Zhu, S.; Li, Y.; Qi, C., Silver-catalyzed difunctionalization of terminal alkynes: highly regio- and stereoselective synthesis of (Z)-beta-haloenol acetates. Org. Lett. 2010, 12, 3262-3265. 47 Yazaki, R.; Kumagai, N.; Shibasaki, M., Direct catalytic asymmetric conjugate addition of terminal alkynes to alpha,beta-unsaturated thioamides. J. Am. Chem. Soc. 2010, 132, 10275-10277. 48 Xiao, Q.; Xia, Y.; Li, H.; Zhang, Y.; Wang, J., Coupling of NTosylhydrazones with Terminal Alkynes Catalyzed by Copper(I): Synthesis of Trisubstituted Allenes. Angew. Chem. Int. Ed. 2011, 50, 1114-1117. 49 Ye, F.; Ma, X.; Xiao, Q.; Li, H.; Zhang, Y.; Wang, J., C(sp)-C(sp3) bond formation through Cu-catalyzed cross-coupling of N-tosylhydrazones and trialkylsilylethynes. J. Am. Chem. Soc. 2012, 134, 5742-5745.

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

50 Jie, X.; Shang, Y.; Hu, P.; Su, W., Palladium-Catalyzed Oxidative CrossCoupling between Heterocycles and Terminal Alkynes with Low Catalyst Loading. Angew. Chem. Int. Ed. 2013, 52, 3630-3633. 51 Xiao, Q.; Zhang, Y.; Wang, J., Diazo compounds and N-tosylhydrazones: novel cross-coupling partners in transition-metal-catalyzed reactions. Acc. Chem. Res. 2013, 46, 236-247. 52 Ye, F.; Wang, C.; Ma, X.; Hossain, M. L.; Xia, Y.; Zhang, Y.; Wang, J., Synthesis of Terminal Allenes through Copper-Mediated Cross-Coupling of Ethyne with N-Tosylhydrazones or alpha-Diazoesters. J. Org. Chem. 2015, 80, 647-652. 53 Zhang, J.; Chen, H.; Lin, C.; Liu, Z.; Wang, C.; Zhang, Y., Cobalt-Catalyzed Cyclization of Aliphatic Amides and Terminal Alkynes with Silver-Cocatalyst. J. Am. Chem. Soc. 2015, 137, 12990-12996. 54 Hossain, M. L.; Ye, F.; Zhang, Y.; Wang, J., CuI-catalyzed cross-coupling of N-tosylhydrazones with terminal alkynes: synthesis of 1,3-disubstituted allenes. J. Org. Chem. 2013, 78, 1236-1241. 55 Chu, W. D.; Zhang, L.; Zhang, Z.; Zhou, Q.; Mo, F.; Zhang, Y.; Wang, J., Enantioselective Synthesis of Trisubstituted Allenes via Cu(I)-Catalyzed Coupling of Diazoalkanes with Terminal Alkynes. J. Am. Chem. Soc. 2016, 138, 14558-14561. 56 Black, D. A.; Arndtsen, B. A., Copper-catalyzed coupling of imines, acid chlorides, and alkynes: a multicomponent route to propargylamides. Org. Lett. 2004, 6, 1107-1110. 57 Park, S. B.; Alper, H., An efficient synthesis of propargylamines via C-H activation catalyzed by copper(I) in ionic liquids. Chem. Commun. 2005, 1315-1317. 58 Gao, Y.; Wang, G.; Chen, L.; Xu, P.; Zhao, Y.; Zhou, Y.; Han, L. B., Coppercatalyzed aerobic oxidative coupling of terminal alkynes with H-phosphonates leading to alkynylphosphonates. J. Am. Chem. Soc. 2009, 131, 7956-7957. 59 Hattori, G.; Sakata, K.; Matsuzawa, H.; Tanabe, Y.; Miyake, Y.; Nishibayashi, Y., Copper-catalyzed enantioselective propargylic amination of propargylic esters with amines: copper-allenylidene complexes as key intermediates. J. Am. Chem. Soc. 2010, 132, 10592-10608. 60 Hein, J. E.; Fokin, V. V., Copper-catalyzed azide-alkyne cycloaddition (CuAAC) and beyond: new reactivity of copper(I) acetylides. Chem. Soc. Rev. 2010, 39, 1302-1315. 61 Makarem, A.; Berg, R.; Rominger, F.; Straub, B. F., A Fluxional Copper Acetylide Cluster in CuAAC Catalysis. Angew. Chem. Int. Ed. 2015, 54, 7431-7435. 62 Cheung, K. P. S.; Tsui, G. C., Copper(I)-Catalyzed Interrupted Click Reaction with TMSCF3: Synthesis of 5-Trifluoromethyl 1,2,3-Triazoles. Org. Lett. 2017, 19, 2881-2884. 63 Jin, L.; Hao, W.; Xu, J.; Sun, N.; Hu, B.; Shen, Z.; Mo, W.; Hu, X., NHeterocyclic carbene copper-catalyzed direct alkylation of terminal alkynes with nonactivated alkyl triflates. Chem. Commun. 2017, 53, 4124-4127. 64 Nallagangula, M.; Namitharan, K., Copper-Catalyzed Sulfonyl AzideAlkyne Cycloaddition Reactions: Simultaneous Generation and Trapping of CopperTriazoles and -Ketenimines for the Synthesis of Triazolopyrimidines. Org. Lett. 2017, 19, 3536-3539. 65 Sakamoto, R.; Kato, T.; Sakurai, S.; Maruoka, K., Copper-Catalyzed C(sp)C(sp(3)) Coupling of Terminal Alkynes with Alkylsilyl Peroxides via a Radical Mechanism. Org. Lett. 2018, 20, 1400-1403. 66 Davies, H. M.; Beckwith, R. E., Catalytic enantioselective C-H activation by means of metal-carbenoid-induced C-H insertion. Chem. Rev. 2003, 103, 2861-2904. 67 Davies, H. M.; Manning, J. R., Catalytic C-H functionalization by metal carbenoid and nitrenoid insertion. Nature 2008, 451, 417-424. 68 Barluenga, J.; Lonzi, G.; Riesgo, L.; Lopez, L. A.; Tomas, M., Pyridine activation via copper(I)-catalyzed annulation toward indolizines. J. Am. Chem. Soc. 2010, 132, 13200-13202. 69 Gillingham, D.; Fei, N., Catalytic X-H insertion reactions based on carbenoids. Chem. Soc. Rev. 2013, 42, 4918-4931. 70 Besora, M.; Braga, A. C.; Sameera, W. M. C.; Urbano, J.; Fructos, M.; Perez, P. J.; Maseras, F., A computational view on the reactions of hydrocarbons with coinage metal complexes. J. Organomet. Chem. 2015, 784, 2-12. 71 Liu, Y.; Luo, Z.; Zhang, J. Z.; Xia, F., DFT Calculations on the Mechanism of Transition-Metal-Catalyzed Reaction of Diazo Compounds with Phenols: O-H Insertion versus C-H Insertion. J. Phys. Chem. A 2016, 120, 6485-6492.

72 Yu, J.; Zhou, Y.; Lin, Z.; Tong, R., Regioselective and Stereospecific CopperCatalyzed Deoxygenation of Epoxides to Alkenes. Org. Lett. 2016, 18, 4734-4737. 73 Fructos, M. R.; Besora, M.; Braga, A. A. C.; Díaz-Requejo, M. M.; Maseras, F.; Perez, P. J., Mechanistic Studies on Gold-Catalyzed Direct Arene C–H Bond Functionalization by Carbene Insertion: The Coinage-Metal Effect. Organometallics 2017, 36, 172-179. 74 Xia, Y.; Qiu, D.; Wang, J., Transition-Metal-Catalyzed Cross-Couplings through Carbene Migratory Insertion. Chem. Rev. 2017, 117, 13810-13889. 75 Lv, X.; Kang, Z.; Xing, D.; Hu, W., Cu(I)-Catalyzed Three-Component Reaction of Diazo Compound with Terminal Alkyne and Nitrosobenzene for the Synthesis of Trifluoromethyl Dihydroisoxazoles. Org. Lett. 2018, 20, 4843-4847. 76 Frisch, M. J. T., G. W.; Schlegel, H. B.; Scuseria, G. E.;; Robb, M. A. C., J. R.; Scalmani, G.; Barone, V.; Mennucci,; B.; Petersson, G. A. N., H.; Caricato, M.; Li, X.; Hratchian, H.; P.; Izmaylov, A. F. B., J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.;; Ehara, M. T., K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima,; T.; Honda, Y. K., O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.;; Peralta, J. E. O., F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin,; K. N.; Staroverov, V. N. K., R.; Normand, J.; Raghavachari, K.;; Rendell, A. B., J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega,; N.; Millam, J. M. K., M.; Knox, J. E.; Cross, J. B.; Bakken, V.;; Adamo, C. J., J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.;; Austin, A. J. C., R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.;; Morokuma, K. Z., V. G.; Voth, G. A.; Salvador, P.;; Dannenberg, J. J. D., S.; Daniels, A. D.; Farkas, O.;; Foresman, J. B. O., J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09,; revision D.01; Gaussian, I. W., CT, 2013. 77 Becke, A. D., Density‐functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648-5652. 78 Lee, C.; Yang, W.; Parr, R. G., Development of the Colle-Salvetti correlationenergy formula into a functional of the electron density. Phys. Rev. B. 1988, 37, 785789. 79 Hopmann, K. H., How Accurate is DFT for Iridium-Mediated Chemistry. Organometallics 2016, 35, 3795-3807. 80 Peverati, R.; Truhlar, D. G., Improving the Accuracy of Hybrid Meta-GGA Density Functionals by Range Separation. J. Phys. Chem. Lett. 2011, 2, 2810-2817. 81 Peverati, R.; Truhlar, D. G., M11-L: A Local Density Functional That Provides Improved Accuracy for Electronic Structure Calculations in Chemistry and Physics. J. Phys. Chem. Lett. 2012, 3, 117-124. 82 Marenich, A. V.; Ho, J.; Coote, M. L.; Cramer, C. J.; Truhlar, D. G., Computational electrochemistry: prediction of liquid-phase reduction potentials. Phys. Chem. Chem. Phys. 2014, 16, 15068-15106. 83 Lin, Y. S.; Tsai, C. W.; Li, G. D.; Chai, J. D., Long-range corrected hybrid metageneralized-gradient approximations with dispersion corrections. J. Chem. Phys. 2012, 136, 154109. 84 Peverati, R.; Truhlar, D. G., Performance of the M11 and M11-L density functionals for calculations of electronic excitation energies by adiabatic time-dependent density functional theory. Phys. Chem. Chem. Phys. 2012, 14, 11363-11370. 85 Steckel, J. A., Ab Initio Calculations of the Interaction between CO2 and the Acetate Ion. J. Phys. Chem. A 2012, 116, 11643-11650. 86 Zhao, Y.; Ng, H. T.; Peverati, R.; Truhlar, D. G., Benchmark Database for Ylidic Bond Dissociation Energies and Its Use for Assessments of Electronic Structure Methods. J. Chem. Theory Comput. 2012, 8, 2824-2834. 87 Cancès, E.; Mennucci, B.; Tomasi, J., A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics. J. Chem. Phys. 1997, 107, 3032-3041. 88 Cossi, M.; Barone, V.; Cammi, R.; Tomasi, J., Ab initio study of solvated molecules: a new implementation of the polarizable continuum model. Chem. Phys. Lett. 1996, 255, 327-335. 89 Johnson, E. R.; Keinan, S.; Mori-Sanchez, P.; Contreras-Garcia, J.; Cohen, A. J.; Yang, W., Revealing Noncovalent Interactions. J. Am. Chem. Soc. 2010, 132, 64986506. 90 Lu, T.; Chen, F., Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580-592.

ACS Paragon Plus Environment

Page 8 of 9

Page 9 of 9 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

TOC

ACS Paragon Plus Environment