Acta Chimica Sinica ›› 2021, Vol. 79 ›› Issue (7): 948-952.DOI: 10.6023/A21040172 Previous Articles    

Article

Cu催化偶联反应合成烷基芳基醚的DFT机理研究

满清敏a,b, 付尊蕴b, 刘甜甜b, 郑明月b,*(), 蒋华良a,b,*()   

  1. a 中国科学技术大学 化学物理系 合肥 230026
    b 中国科学院上海药物研究所 药物发现与设计中心 上海 201203
  • 投稿日期:2021-04-22 发布日期:2021-06-02
  • 通讯作者: 郑明月, 蒋华良
  • 基金资助:
    国家自然科学基金面上项目(81773634)

DFT Mechanism of Cu Catalyzed Coupling Reaction to Alkyl Aryl Ethers

Qingmin Mana,b, Zunyun Fub, Tiantian Liub, Mingyue Zhengb(), Hualiang Jianga,b()   

  1. a Department of Chemical Physics, University of Science and Technology of China, Hefei 230026, China
    b Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Shanghai 201203, China
  • Received:2021-04-22 Published:2021-06-02
  • Contact: Mingyue Zheng, Hualiang Jiang
  • Supported by:
    National Natural Science Foundation of China(81773634)

The mechanism of Cu catalyzed intramolecular Ullmann C―O coupling reaction of 2-(2-iodophenyl)ethan-1-ol was proposed and verified by density functional theory (DFT). All of the structures have been fully optimized using N,N-dimethylformamide as solvent by the Gaussian16 package with the B3LYP functional. In the process of optimization, the standard 6-31G (d) basis set was used for C, O, N and H atoms, while Stuttgart basis set (SDD) was used for Cu and I atoms. Grimme D3 method was used for dispersion correction of all calculations. The frequency calculation is carried out at the same theoretical level to determine whether the calculated geometry is minimum (zero imaginary frequency) or transition states (one imaginary frequency). The extra single-point calculations, including solvation and dispersion corrections, on the optimized geometries were employed to obtain improved Gibbs energy values with B3LYP/def2-tzvp in the continuum solvation model based on density (SMD). The intrinsic reaction coordinates (IRC) of all transition states were calculated to confirm that the structure connects the reactant and the product. The calculation shows that catalytic cycle starts from the in-situ formation of cuprous complex, and cuprous complex undergo oxidative addition to give Cu(III) complexes. The Cu(III) complex is coordinated with alcohol to form a four-coordinated Cu(III) complex. Then the substrate is activated by O-halogen bond, and the intermediates underwent ligand exchange to form another Cu(III) complexes, which transmits the desired alkyl aryl ethers and regenerates the catalyst cycle. Quantitative calculation shows that the reaction mechanism can be divided into four steps: (1) oxidative addition, (2) coordination effect, (3) ligand exchange, (4) reductive elimination. First, it was found that tert-butoxide was not only a base, but also a catalyst species. Secondly, it is pointed out that ligand exchange is the rate-determining step of the whole reaction. Before the ligand exchange, it is activated by a slowly approaching tert-butoxide to form a halogen bond (XB). The calculated results are in good agreement with the experimental data, which proves the rationality of the calculation mechanism and helps to explain the Ullmann reaction mechanism.

Key words: density functional theory, Cu catalyzed, Ullmann C―O coupling, tert-butoxide