Predicting bond dissociation energies of cyclic hypervalent halogen reagents using DFT calculations and graph attention network model

Introduction

Hypervalent iodine reagents are increasingly gaining attention in the fields of organic synthesis and catalysis due to their environmental benefits, accessibility, and cost-efficiency . Over the last three decades, a series of cyclic hypervalent iodine(III) reagents has been developed (Figure 1), including the well-known Zhdankin reagents and Togni reagents . These reagents are popularly used as electrophilic group transfer reagents in a variety of reactions, such as C–H functionalization , unsaturated alkane addition , and cyclization .

[1860-5397-20-127-1]

Figure 1: Examples of cyclic hypervalent halogen reagents.

Despite the rapid development of hypervalent iodine(III) reagents, the exploration of isoelectronic hypervalent bromine(III) and chlorine(III) reagents has been comparatively limited despite their demonstrated potential for unique applications . For example, hypervalent bromine(III) reagents enable C–H amination and alkene aziridination reactions without the need for additional Lewis acid activation . However, challenges in the synthesis and stabilization of cyclic hypervalent bromine and chlorine reagents have impeded their development relative to their iodine(III) analogs . Cyclic hypervalent bromine(III) reagents were pioneered by Miyamoto and have since been developed to a certain extent . Biphenyl hypervalent bromine(III) reagents have been synthesized by Yoshida and Wencel-Delord (Figure 1). Cyclic hypervalent chlorine(III) reagents with similar skeletal structures have not been reported yet, and only biphenyl hypervalent chlorine(III) reagents and cyclic diaryliodonium salts have been synthesized.

Previous investigations have highlighted the critical role of bond dissociation energy (BDE) in understanding the group transfer capabilities and chemical stability of hypervalent iodine(III) reagents. In this context, detailed knowledge of the BDE of hypervalent bromine(III) and chlorine(III) reagents is especially crucial for designing novel reagents. Yet, the BDE values of hypervalent bromine(III) and chlorine(III) reagents remain largely elusive, hampering the design and synthesis of novel reagents.

In recent years, machine learning has emerged as a promising and cost-effective alternative to traditional DFT calculations for predicting key properties of organic molecules such as BDE, nucleophilicity, and electrophilicity . Recently, applications of the Elastic Net model with Avalon fingerprints and the deployment of artificial neural network (ANN) models with the Mordred cheminformatics package have demonstrated considerable success in predicting the BDEs of hypervalent iodine(III) reagents. However, previous studies have been limited to the prediction of hypervalent iodine(III) reagents. Driven by their proven effectiveness and our ongoing interest in hypervalent halogen chemistry , we are motivated to develop a machine learning model for a broader array of cyclic hypervalent halogen reagents, thereby integrating different halogen centers and making it easier to predict the group transfer capacity and chemical stability of different cyclic hypervalent halogen reagents.

Results and Discussion

We selected five different skeletons and twenty common transfer groups for combination (Figure 2) and calculated their BDEs. Referring to the previous computational studies of hypervalent iodine and the computational database of organic species by Paton and co-workers , geometry optimizations and single point energy calculations for homolytic BDEs are both performed using M06-2X/def2-TZVPP in the gas phase at 298.15 K by Gaussian 16 . Frequency calculations confirmed that optimized structures are minima (no imaginary frequency). The accuracy of computational BDEs of halides using M06-2X/def2-TZVPP is also evaluated and compared to experimental BDEs, demonstrating the reliability of the method (see Supporting Information File 1).

[1860-5397-20-127-2]

Figure 2: Common cyclic hypervalent halogen skeletons and transfer groups.

The computational homolytic BDEs are presented in Table 1. From the perspective of halogen centers, hypervalent iodine(III) reagents exhibit the highest homolytic BDEs, followed by hypervalent bromine(III) reagents, while hypervalent chlorine(III) reagents have the lowest. Generally, the homolytic BDEs of cyclic hypervalent iodine(III) reagents are above 30.0 kcal/mol, consistent with their good chemical stability. The homolytic BDEs of some cyclic hypervalent bromine(III) and most cyclic hypervalent chlorine(III) reagents are below 20 kcal/mol, implying these reagents should be too reactive to be isolated. From the perspective of transfer groups, the homolytic BDEs of groups with strong trans effects such as -F, -CCH, -CN, -OCF3, -OTf, -OTs are elevated, while those of -N3, -NH2, -SCF3, etc. are smaller. These results are consistent with our previous studies on the group transfer ability of hypervalent iodine(III) reagents . According to the calculation results, skeleton 5 may be a better candidate for synthesizing cyclic hypervalent bromine(III) and chlorine(III) reagents. The groups with strong trans effects, such as -F, -CCH, -CN, -OTf, can help stabilize cyclic hypervalent bromine(III) and chlorine(III) reagents.

Table 1: Computational homolytic BDEs (kcal/mol) of cyclic hypervalent halogen reagents.

[Graphic 1]   1-X 2-X 3-X 4-X 5-X R X =
I X =
Br X =
Cl X =
I X =
Br X =
Cl X =
I X =
Br X =
Cl X =
I X =
Br X =
Cl X =
I X =
Br X =
Cl a F 78.6 50.5 30.5 81.9 56.6 35.5 80.8 54.3 33.3 80.6 61.4 29.4 84.0 62.0 42.9 b Cl 55.3 29.4 10.1 58.1 35.6 14.0 57.4 33.1 12.1 57.1 40.9 10.7 60.1 40.3 21.3 c Br 44.7 19.8 0.2 47.3 24.9 4.3 46.7 23.0 2.4 46.5 31.2 2.0 49.3 30.2 11.4 d CH3 33.2 13.4 −0.3 42.4 28.7 16.5 42.8 29.6 20.0 49.1 49.3 40.0 39.9 28.9 22.2 e CF3 33.5 13.5 −0.8 39.0 24.5 11.8 38.5 23.6 12.5 41.2 39.0 25.1 37.2 24.5 12.6 f CHCH2 40.4 21.4 8.7 49.6 36.5 25.8 49.6 36.8 28.1 56.0 58.2 47.1 46.8 36.1 27.0 g CCH 66.1 42.0 24.5 73.1 53.9 38.2 72.6 53.0 39.2 76.3 68.5 48.5 71.1 53.8 39.0 h CN 68.8 43.6 24.2 72.6 51.6 33.4 71.6 49.7 32.7 72.8 61.6 37.9 71.9 53.0 34.8 i N3 32.1 7.6 −11.2 35.8 14.6 −4.5 34.8 12.3 −6.3 36.3 23.3 −4.8 36.5 17.9 0.1 j NH2 37.7 12.6 −4.3 45.5 25.1 9.6 45.0 24.1 8.3 49.4 40.2 16.7 44.5 27.5 11.7 k NHAc 47.6 22.0 2.6 53.4 31.6 13.0 52.7 30.0 12.2 55.8 43.4 18.8 53.0 33.7 15.7 l OH 53.0 26.4 6.5 58.7 35.7 15.9 57.7 33.5 14.1 60.1 44.9 16.3 58.7 38.3 19.5 m OCH3 40.7 15.3 −3.5 46.2 24.6 6.1 45.2 22.4 4.3 47.9 34.5 7.7 46.2 26.9 9.1 n OCF3 62.9 37.0 18.4 64.4 40.5 19.8 63.5 38.1 17.4 62.7 45.0 13.3 67.0 46.6 28.6 o OCOCH3 54.6 27.4 7.5 58.1 33.7 12.3 57.5 31.1 9.9 57.6 39.3 8.6 59.1 37.1 17.4 p OCOCF3 62.3 36.1 17.4 63.5 39.0 18.0 62.7 36.7 15.8 61.2 43.4 11.3 65.9 45.1 27.0 q OCOPh 55.8 28.5 8.7 58.9 33.8 12.4 58.4 31.8 10.4 58.4 40.5 9.0 60.1 38.0 19.1 r OTf 66.3 41.9 25.8 65.1 41.3 21.3 64.3 39.0 18.7 62.0 43.6 11.1 69.8 50.8 35.4 s OTs 61.8 36.5 18.4 62.4 38.4 18.3 62.3 36.2 15.5 61.0 43.0 12.9 65.8 45.8 28.9 t SCF3 40.1 16.1 −3.1 44.6 23.3 4.0 43.2 21.6 3.0 44.9 33.7 14.0 44.2 26.2 7.9

In addition, we also calculated the heterolytic BDEs of cyclic hypervalent halogen reagents to comprehensively examine the strength of chemical bonds (Table 2). Geometry optimizations and single point energy calculations for heterolytic BDEs are performed using M06-2X/def2-TZVPP in the SMD (acetonitrile) Implicit solvent model at 298.15 K. Due to the instability of some transfer group cations, such as +OCH3, +OCF3, +OCOCF3, +OCOPh, +OTf and +SCF3, it is difficult for us to investigate their heterolytic BDEs. From Table 2, it can be seen that, except for CF3 and CHCH2, all other transfer groups exhibit high heterolytic BDEs with hypervalent halogen centers.

Table 2: Computational heterolytic BDEs (kcal/mol) of cyclic hypervalent halogen reagents.

[Graphic 2]   1-X 2-X 3-X 4-X 5-X R X =
I X =
Br X =
Cl X =
I X =
Br

留言 (0)

沒有登入
gif