Open Access Article
Farshad Shiria,
Morteza Jamshidib,
Saba Hadidi*c,
Robert Stranger
b and
Alireza Ariafard
*b
aDepartment of Chemistry, The Hong Kong University of Science and Technology, Kowloon, Hong Kong 999077, P. R. China
bResearch School of Chemistry, Australian National University, Canberra, Australian Capital Territory 2601, Australia. E-mail: alireza.ariafard@utas.edu.au
cInorganic Chemistry Department, Faculty of Chemistry, Razi University, Kermanshah 6714414971, Iran. E-mail: s.hadidi@razi.ac.ir
First published on 10th December 2025
This work presents a detailed DFT-based mechanistic investigation of copper-photocatalyzed atom transfer radical addition (ATRA) reactions between CF3SO2Cl and alkenes. Depending on the electronic nature of the alkene substrate, these reactions yield either RCl or RSO2Cl products. The unusual divergence in product selectivity has led to the proposal of multiple mechanistic pathways. In this study, we show that all productive pathways proceed exclusively via outer-sphere single-electron transfer and identify two previously unrecognized mechanisms: an S(VI)/S(IV) redox cycling mechanism responsible for RSO2Cl formation, and a 2c–3e Cl-coordination-induced SET mechanism accounting for RCl formation. These two pathways represent the first models to explicitly demonstrate the bifunctional role of the [SO2Cl]− anion in governing divergent product formation. Additionally, we identify a third, cationic mechanism, in which the carbon-centred radical is oxidized to a carbocation by Cu(II), competing with the other pathways and likewise leading to RCl. Taken together, these results provide a useful framework for understanding chemoselectivity in this class of photocatalytic transformations and may help guide the design of future ATRA protocols.
Among photoredox catalysis methods, the copper photocatalyzed atom-transfer radical addition (ATRA) reactions of CF3SO2Cl (triflyl chloride) with alkenes, independently developed by Dolbier et al.39 and Reiser et al.,40 stand out as groundbreaking and original discoveries in the field (Fig. 1a–d). Dolbier et al. reported that the irradiation of electron-deficient alkenes, such as substrate S1, and triflyl chloride in the presence of a [Cu(dap)2]Cl catalyst (Fig. 1e) results in trifluoromethylchlorination with the extrusion of SO2, yielding P1(S1) as the sole product (Fig. 1a).39 In contrast, Reiser et al. observed that the irradiation of electron-neutral alkenes, such as substrate S2, and triflyl chloride, catalyzed by [Cu(dap)2]Cl, results in trifluoromethylchlorosulfonylation without SO2 extrusion, yielding P2(S2) as the major product (Fig. 1b). However, Reiser et al. observed that the product distribution is highly sensitive to the electronic nature of alkenes. For example, replacing the phenyl ring in substrate S2 with a para-aminophenyl group to form substrate S3 led to the selective formation of product P1 (Fig. 1c). Similarly, substrate S4, characterized by a disubstituted alkene at the C1 position, exclusively yielded product P1 (Fig. 1d).40
Building upon the pioneering work of Dolbier et al. and Reiser et al., numerous studies have advanced this specific light-driven transformation, catalysed not only by copper complexes but also activated by alternative systems, leading to significant breakthroughs in synthetic methodology.41–67 However, despite substantial progress on the synthetic side, the mechanistic understanding of this process remains limited.
Various reaction mechanisms, based on both computational68 and experimental39,40,69–72 findings, have been proposed to explain the strikingly divergent product outcomes that arise from subtle changes in alkene structure. During the preparation of this work, a related computational study by Pham et al.68 was published. While their study offers valuable insights into this class of reactions, our independent investigation supports a different mechanistic picture, which is explored in detail herein.
In the following section, we briefly review the mechanism proposed by Pham et al.68 as background to our investigation. In their study, the reaction is proposed to initiate with photoexcitation of the Cu(I) complex [Cu(dap)2]+, followed by single-electron transfer (SET) to triflyl chloride via either an inner-sphere (ISET) or outer-sphere (OSET) pathway, as illustrated in Fig. 1f. Among these, the inner-sphere mechanism was reported to be more favourable, generating a ˙CF3 radical and leaving a [SO2Cl]− anion coordinated to the in situ generated Cu(II) atom (intermediate A in Fig. 1f). The ˙CF3 radical subsequently reacts with the alkene substrate S to form the R˙ radical. They also proposed that a ligand exchange between a Cl− anion and the [SO2Cl]− ligand in intermediate A, yielding the more stable intermediate C. The R˙ radical can subsequently react via two distinct pathways: (a) with the [SO2Cl]− ligand in intermediate A, leading to the formation of product P2 (RSO2Cl), or (b) with the Cl ligand in intermediate C, resulting in the formation of product P1 (RCl). The product distribution has been proposed to depend on the
value (Fig. 1f). A relatively small
favours the formation of product P2, whereas a higher
enables the ligand exchange process to occur, leading to the formation of intermediate C, which subsequently reacts with the R˙ radical to afford product P1. However, this conclusion was reached without calculating the transition structures connecting A to B and B to C.
To evaluate the accuracy of this claim, we successfully located the crucial transition structures in this study and found that A connects to B through a remarkably low activation barrier of only 1.7 kcal mol−1 (vide infra, Fig. 5). Consequently, the conversion of A to B occurs significantly faster, effectively outcompeting the reaction of R˙ with A. It follows that if the reaction were to proceed via the inner-sphere mechanism proposed by Pham et al.,68 product P1 would be formed exclusively, regardless of the alkene substrate. This outcome does not align with the full set of experimental findings, suggesting that additional, previously unrecognised mechanisms may be needed to account for the observed product distributions.
Our efforts to identify new pathways for the formation of P1 and P2 have led to the discovery of two unprecedented mechanisms. We will show in this study that in these two mechanisms, the [SO2Cl]− anion, generated in situ via outer sphere electron transfer from the excited Cu complex to CF3SO2Cl, plays a pivotal role as the key intermediate. This anion can interact with the R˙ radical through two distinct pathways:
Mechanism (i) (S(VI)/S(IV) redox cycling, Fig. 1h): The sulfur atom of [SO2Cl]− binds to R˙, oxidizing sulfur from a formal oxidation state of +4 to +5 and forming intermediate E. This highly reactive intermediate then transfers an electron to the Cu(II) atom via an outer-sphere mechanism, leading to the formation of product P2 and regenerating the Cu(I) catalyst.
Mechanism (ii) (2c–3e Cl-coordination-induced SET, Fig. 1i): The [SO2Cl]− anion binds to R˙ via its Cl atom, forming a 2c–3e bond in intermediate F. This intermediate subsequently transfers an electron to the Cu(II) atom through an outer-sphere mechanism, yielding product P1, accompanied by the release of SO2 and the regeneration of the Cu(I) catalyst.
The formation of a carbocation through outer-sphere electron transfer from R˙ to a Cu(II) complex, followed by its trapping by an available anion, represents another plausible pathway for generating the desired products (mechanism (iii)). Pham et al.68 investigated this pathway using five-coordinate Cu(II) complexes, [CuIICl] and [CuIISO2Cl], as oxidants and concluded that it was unfavourable, ruling out its viability (Fig. 1g). However, in this study, we demonstrate that five-coordinate Cu(II) complexes are significantly weaker oxidants compared to the four-coordinate Cu(II) complex [CuII]. As a result, the formation of carbocations becomes energetically much more favourable when the more potent four-coordinate Cu(II) complex [CuII] serves as the oxidant. In this study, we assert that this carbocationic pathway is indeed the operative mechanism when substrates can generate very stable carbocations.
As mentioned above, it was concluded by Pham et al. that the electron transfer between the photoexcited complex [Cu(dap)2]+* and CF3SO2Cl proceeds via an inner-sphere mechanism.68 In contrast, we demonstrate here that this is not the case. Once CF3SO2Cl and [Cu(dap)2]+* form an outer-sphere complex (the initial step of an inner-sphere mechanism), an electron is immediately transferred from [Cu(dap)2]+* to CF3SO2Cl with an activation energy as low as 0.2 kcal mol−1 (vide infra, Fig. 3c). This result suggests that an inner-sphere mechanism for electron transfer is unlikely. Here, we will demonstrate that all processes, including the electron transfer from [Cu(dap)2]+* to CF3SO2Cl and the formation of products P1 and P2, occur exclusively through outer-sphere mechanisms.
From S1*, we considered two possible deactivation pathways: (i) intersystem crossing (ISC) to the slightly more stable triplet-state complex T1* via the minimum energy crossing point (MECP) MECP1 and (ii) non-radiative relaxation back to the ground-state complex 1 via MECP2. MECP1 is calculated to be 6.1 kcal mol−1 lower in energy than MECP2, indicating that ISC to T1* is energetically favoured and that T1* should be readily populated upon photoexcitation of 1.
From T1*, we then evaluated two subsequent pathways: (i) outer-sphere electron transfer from T1* to CF3SO2Cl via the Marcus crossing point (MCP) MCP1 and (ii) relaxation back to the ground-state complex 1 through MECP3. MECP3 lies 12.7 kcal mol−1 above MECP2, suggesting that, once formed, T1* is much less prone to non-radiative relaxation than S1*.
MCP1 is found to lie 1.6 kcal mol−1 higher in energy than MECP2, indicating a competition between relaxation and electron transfer. Although relaxation from the excited state may also occur through radiative pathways, this small energy difference is reasonably consistent with the experimental quantum yield of 12% reported by Reiser et al.40
It should finally be noted that, although MECPs do not provide quantitative ISC rates, their relative positions nevertheless offer valuable insight into the accessibility of crossing events. In practice, ISC rates are determined not only by the location of MECPs but also by additional factors such as spin–orbit coupling,77 yet MECP energetics remain a useful qualitative tool for identifying feasible crossings in photochemical mechanisms.
Based on our calculations at the SMD/B3LYP-D3/def2-SVP level of theory, all attempts to optimize the outer-sphere complex G (Fig. 3b) consistently resulted in spontaneous electron transfer from [Cu]* to CF3SO2Cl, yielding intermediate 5. This is evidenced by the increase in total spin density on CF3SO2Cl in intermediate 5, reaching 0.653, along with the elongation of the calculated S–Cl and S–CF3 bond distances from 2.112 and 1.910 Å in 2 to 2.727 and 1.912 Å in 5, respectively. The electron transfer renders the formation of the outer-sphere complex 5 highly exergonic, with a calculated free energy release of approximately 23.5 kcal mol−1 (Fig. 3b). This finding supports the conclusion that the SET process should proceed via an outer-sphere mechanism, because in the assumed outer-sphere complex G the two redox partners remain separated, without the orbital overlap required for an inner-sphere SET pathway.
To further validate the conclusion that electron transfer between [Cu]* and CF3SO2Cl occurs essentially immediately once the outer-sphere complex G is formed, we employed a larger basis set (BS3) for its optimization (BS3: def2-TZVP for Cu, Cl and S, and def2-SVP for other atoms). Interestingly, this time we successfully located outer-sphere complex G; however, we found that it readily converts to the more stable complex 5 with an activation energy as low as 0.2 kcal mol−1 (Fig. 3c). This much lower barrier (0.2 kcal mol−1), compared with the ∼8 kcal mol−1 barrier obtained for MCP1 in Fig. 2, can be rationalized by Marcus theory: at longer donor–acceptor distances a substantial barrier exists, whereas at near-contact the reduced distance between the redox partners greatly increases the probability of electron transfer, making the process essentially barrierless. This behaviour is fully consistent with Marcus theory, in which the outer-sphere SET barrier decreases as the distance between the redox partners is reduced.79,80
An analysis of the frontier orbitals of complex G at the SMD/B3LYP/def2-TZVP level of theory provides insight into why this complex is highly reactive toward the SET process through the outer-sphere mechanism (Fig. 3d). The HOMO in this complex is primarily characterized by a π* orbital of the dap ligand that holds the unpaired electron, while the LUMO is predominantly composed of the σ* orbitals of the S–Cl bond (major contribution) and the S–CF3 bond (minor contribution) in CF3SO2Cl. The energy gap between the HOMO and LUMO in complex G is only 0.27 eV, which is sufficiently small to explain why the single electron transfer from the π* orbital of dap to the lowest lying σ* orbital of CF3SO2Cl occurs almost spontaneously.
From complex 5, we investigated two different pathways for the formation of ˙CF3: (i) a Cu(II)-assisted dissociation pathway and (ii) a direct dissociation pathway (Fig. 3b).
In the Cu(II)-assisted pathway, the [CF3SO2Cl]˙− radical anion must first coordinate to the Cu(II) centre, forming complex I. For this coordination to occur, the four-coordinate Cu(II) complex, which adopts a distorted square planar geometry76 in complex 5, must undergo a rearrangement to a trigonal pyramidal geometry to create an open coordination site for [CF3SO2Cl]˙−. This structural rearrangement is energetically demanding, making complex I 5.6 kcal mol−1 less stable than complex 5. The release of the ˙CF3 radical from complex I proceeds through the transition structure TSI–A, with a relative free energy of 22.0 kcal mol−1, leading to the formation of complex A. This radical dissociation process induces a formal oxidation state change of the sulfur atom from +5 to +4.
In the direct dissociation pathway, the outer-sphere complex 5 undergoes fragmentation to yield the Cu(II) complex B and the S(V) species 3 in a thermodynamically favourable process. From intermediate 3, we investigated two possible pathways: (a) dissociation of Cl− through TS3 and (b) dissociation of ˙CF3 through TS3-D. The calculations clearly indicate that 3 is significantly more reactive toward ˙CF3 dissociation, as evidenced by the fact that TS3-D lies 8.0 kcal mol−1 lower in energy than TS3.
From a comparison of the two pathways presented in Fig. 3b, it is evident that TS3-D is approximately 15.6 kcal mol−1 lower in energy than TSI–A, indicating that ˙CF3 release preferentially occurs through the direct dissociation pathway rather than the Cu(II)-assisted pathway. This conclusion is further supported by calculations at the SMD/wB97XD/def2-TZVP//SMD/B3LYP-D3/def2-SVP level of theory, where in this case TS3-D lies about 19.9 kcal mol−1 lower in energy than TSI–A (see Fig. S1).
The results presented in Fig. 3b offer an opportunity to revisit and refine the mechanistic interpretation proposed by Pham et al. in a recent publication.68 In their computational study, they proposed that the transformation I → TSI–A → A + ˙CF3 corresponds to an inner-sphere SET process. However, our detailed analysis reveals that this transformation actually represents the Cu(II)-assisted release of the ˙CF3 radical after the SET has already occurred from [Cu]* to CF3SO2Cl via an outer-sphere mechanism.
As previously proposed,68 once the ˙CF3 radical is formed, it is added to the alkene substrate during the catalytic cycle to form R˙. Here, we begin our discussion by focusing on substrate S2 (Fig. 1b), while the computational analysis of the reactivity of other substrates depicted in Fig. 1a, c, and d will be addressed later. The ˙CF3 radical can selectively add to either the C1 or C2 position of substrate S2. As shown in Fig. 4c, the addition to the C2 atom proceeds with a significantly lower activation barrier compared to C1, a result that aligns well with the experimental observations.68
Pathway A (Fig. 5) begins with the fragmentation of species D into SO2 and Cl−, followed by coordination of the resultant Cl− to the four-coordinate Cu(II) complex B to form [CuIICl] (complex C). Finally, coupling between the R˙ radical and the Cl ligand in complex C on the symmetry-broken open-shell singlet surface facilitates single-electron transfer from the R˙ radical to the Cu(II) centre via transition structure TSC–P1, ultimately affording the product P1(S2) and regenerating the Cu(I) catalyst 1.
Pathway B (Fig. 5) begins with the coordination of species D to the Cu(II) complex B through one of its oxygen atoms, via the transition structure TSB–A, to form complex A. Subsequently, the R˙ radical attacks the sulfur atom in complex A, promoting the single electron on R˙ to transfer to the Cu(II) centre via the transition structure TSA–P2 located on the symmetry-broken open-shell singlet surface, thereby forming the P2(S2) product.
A comparison of pathways A and B reveals that pathway A is energetically more favourable, as the highest energy point along this route (TSC–P1) lies 3.1 kcal mol−1 below the corresponding transition structure in pathway B (TSA–P2). This difference becomes even more pronounced at the SMD/wB97XD/def2-TZVP//SMD/B3LYP-D3/def2-SVP level of theory in acetonitrile, where TSC–P1 lies 8.7 kcal mol−1 lower in energy than TSA–P2 (Fig. S2). These findings suggest that if the final step of the reaction proceeds via a Cu(II)-mediated mechanism, product P1(S2) would be expected as the predominant outcome. However, this prediction does not align with experimental observations, which show that when substrate S2 is employed, product P2 is formed preferentially (Fig. 1b). Together, these results indicate that alternative, more favourable mechanisms beyond those mediated by Cu(II) complexes must be operative to account for the experimentally observed chemoselectivity. In the following sections, we demonstrate how our proposed outer-sphere mechanisms (i)–(iii), as introduced in the Introduction, rationalize the observed selectivity.
Mechanism (i), responsible for the formation of product P2(S2), was introduced earlier in the Introduction (Fig. 1h). In this pathway, the anion D (SO2Cl−) reacts with the radical R˙(S2) via transition structure TSS2, forming intermediate E(S2) (Fig. 6). From this intermediate, an electron is transferred to the Cu(II) complex B through the Marcus crossing point MCP1(S2), ultimately resulting in the formation of product P2(S2) and the regeneration of the Cu(I) catalyst 1.
Mechanism (iii), which leads to the formation of product P1(S2), involves the transfer of an electron from the radical R˙(S2) to the Cu(II) complex B via the Marcus crossing point MCP3(S2), resulting in the formation of the carbocation R+(S2) and the regeneration of the Cu(I) catalyst 1 (Fig. 6). This carbocation gains partial stabilization through interaction with the phenyl substituent, as evidenced by the shortening of the C2–C4 distance from 2.516 Å in R˙(S2) to 1.630 Å in R+(S2); the CASSCF results supported by spin-density analysis suggest that no electronic interaction exists between C2 and C4 in R˙(S2) (Fig. S3). Finally, R+(S2) is trapped by the Cl substituent of SO2Cl− via transition structure TSD–P1, followed by SO2 release, to yield product P1(S2). In contrast, trapping of R+(S2) through the S atom of SO2Cl−, involving transition structure TSD, requires a much higher activation free energy (Fig. 6). This difference reflects the greater nucleophilicity of the Cl substituent compared to the S atom in SO2Cl−, explaining why carbocation R+(S2) preferentially reacts with SO2Cl− to form product P1(S2) rather than product P2(S2).
It follows from Fig. 6 that the highest energy point on mechanism (i), namely MCP1(S2), lies 2.6 kcal mol−1 lower in energy than the corresponding point on mechanism (iii), MCP3(S2). This energy difference indicates that the formation of P2(S2) is more favourable than that of P1(S2), thereby explaining the experimentally observed chemoselectivity for substrate S2, for which P2(S2) is the dominant product (Fig. 1b).
It is worth noting that for substrate S2, mechanism (ii) outlined in Fig. 1i is unlikely to be operative, as we were unable to locate intermediate F featuring the proposed 2c–3e C–Cl bond. This appears to be due to the full localization of the unpaired electron on the C2 atom in R˙(S2), which prevents the orbital overlap required to stabilize a C–Cl 2c–3e bond and thus precludes the formation of intermediate F.
For the electron-deficient substrate S1, we found that mechanism (ii) is indeed feasible. This is because the single electron on the resulting R˙ radical is not exclusively localized on the C2 atom but is delocalized between the C2 and O atoms, as evidenced by the spin density values of 0.852 for C2 and 0.120 for O. This delocalization significantly enhances the accessibility of the C2 atom for interaction with the Cl atom of the SO2Cl− anion (species D), leading to the formation of intermediate F(S1), which features a characteristic C–Cl 2c–3e bond. This intermediate then transfers an electron to the four-coordinate Cu(II) complex B via the Marcus crossing point MCP2(S1) with relative energy of 18.1 kcal mol−1, ultimately yielding product P1(S1).
Although the delocalization of the single electron in radical R˙(S1) facilitates its interaction with the Cl atom of the SO2Cl− anion, it concurrently reduces its ability to oxidize the S(IV) centre in SO2Cl− to S(V). As a result, the formation of intermediate E(S1) is energetically disfavoured, lying 8.4 kcal mol−1 above R˙(S1). This diminished oxidative capability plays a key role in positioning MCP1(S1) 4.4 kcal mol−1 higher in energy than MCP2(S1), thereby explaining why no product P2 is formed when electron-deficient alkenes such as S1 are used in the photocatalytic process.
Mechanism (iii) is highly unfavourable for this substrate due to its electron-deficient nature, which results in MCP3(S1) having a relative energy of 29.9 kcal mol−1, effectively ruling out its feasibility (Fig. 7). Consequently, based on the results presented in Fig. 7, it can be concluded that for the electron-deficient substrate S1, mechanism (ii) is the operative pathway, leading to the exclusive formation of product P1. This finding is in excellent agreement with the experimental observations.68
The CASSCF results supported by spin-density analysis for E(S1) show that the single unpaired electron mainly occupies the σ* orbital of the S–Cl bond, whereas the σ* orbital of the S–C bond is almost completely empty. This finding demonstrates that the formal oxidation state of sulfur in this species is best described as +V. The same analysis for F(S1) reveals that the C–Cl 2c–3e bond is highly polar, with the unpaired electron density residing primarily on the C atom (for details, see Fig. S3).
Fig. 10 compares the free energy profiles for the oxidation of R˙(S3) to R+(S3) by [CuII] and [CuII–Cl]. As shown, the reaction R˙(S3) + [CuII] → R+(S3) + [CuI] is significantly more exergonic (ΔG = −17.0 kcal mol−1) than the corresponding process with [CuII–Cl], which proceeds with a much smaller driving force (ΔG = −4.8 kcal mol−1). This thermodynamic difference is reflected in the position of the Marcus crossing points: MCP3(S3) lies 12.5 kcal mol−1 lower in energy than MCP3′(S3) (Fig. 10). These results clearly demonstrate that [CuII] is a far superior oxidant compared to [CuII–Cl]. The underlying reason is electronic: [CuII] is a 17-electron species that becomes an 18-electron complex upon reduction, an electronically favourable transformation. In contrast, [CuII–Cl] is a 19-electron species; accepting another electron would yield a 20-electron complex, significantly violating the 18-electron rule and making it a much less effective oxidant.
Species 4 then releases the ˙CF3 radical, further reducing the formal oxidation state of sulfur from +5 to +4 in the resulting SO2Cl− anion. The generated ˙CF3 radical subsequently adds to the alkene substrate, forming the radical intermediate R˙, which serves as a key branching point for several plausible mechanistic pathways. We have explored three such pathways, mechanisms (i) to (iii), each capable of explaining the experimentally observed chemoselectivities arising from variations in alkene substrate structure.
Mechanism (i) (SVI/SIV redox cycling – Fig. 11a): The R˙ radical is added to the sulfur atom of SO2Cl− to form intermediate E, oxidizing the sulfur centre from S(IV) to S(V). A subsequent outer-sphere electron transfer from E to the Cu(II) complex B regenerates the Cu(I) catalyst and affords product P2, in which the formal oxidation state of sulfur is restored to +6. As a result, when this mechanism is operative, the sulfur centre undergoes a stepwise redox sequence: S(VI) → S(V) → S(IV) → S(V) → S(VI), thus justifying the designation “S(VI)/S(IV) redox cycling”.
This mechanism is expected to be operative when intermediate E does not lie significantly higher in energy than the R˙ radical, or when the alternative R+ carbocation is not intrinsically very stable. Consequently, tertiary or electron-deficient R˙ radicals are less likely to oxidize S(IV) to S(V), resulting in a significantly less stable intermediate E and making this pathway less likely to operate.
Mechanism (ii) (2c–3e Cl-coordination-induced SET - Fig. 11b): The R˙ radical forms a weak C–Cl 2c–3e bond with the Cl atom of SO2Cl−, generating intermediate F. An electron is then transferred from this intermediate to the Cu(II) complex B via an outer-sphere SET step, affording product P1 along with SO2 extrusion.
This mechanism is found to be operative for electron-deficient alkene substrates, where the radical centre in R˙ is adjacent to a π-acceptor carbonyl group. This electronic arrangement facilitates delocalization of the unpaired electron in R˙, partially unblocking the radical centre for formation of the weak 2c–3e bond and thereby stabilizing intermediate F as a genuine local minimum on the potential energy surface, ultimately enabling mechanism (ii) to proceed.
Mechanism (iii) (Fig. 11c): Alternatively, the R˙ radical directly transfers an electron to the Cu(II) complex B, forming the carbocation R+. This highly reactive species is then trapped by SO2Cl−, yielding product P1 with concurrent SO2 extrusion. We found that this mechanism is operative when the R+ carbocation is intrinsically very stable.
For this pathway to occur, we assumed that Cl− first coordinates to [Cu(dap)2]2+, forming intermediate C (Fig. 12). One dap ligand can then dissociate, yielding structure X1. This dissociation process is endergonic by about 17.6 kcal mol−1. The resulting Cu(II) species X1 can subsequently coordinate with SO2Cl−, either through its Cl atom or through one of the oxygen atoms. Our calculations show that coordination via the Cl atom is much more favourable, giving complex X3, which is 8.1 kcal mol−1 lower in energy than the O-coordinated species X2. However, X3 is highly reactive and readily releases SO2 gas, leading to the formation of X4.
![]() | ||
Fig. 12 Mechanism for the coordination-based pathway leading from [Cu(dap)2]2+ to [Cu(dap)(Cl)(SO2Cl)] and its decomposition products calculated at the SMD/B3LYP-D3/def2-TZVP//SMD/B3LYP-D3/def2-SVP level in acetonitrile. The free-energy barrier for dap dissociation from complex C was estimated following the protocol of Hall and Hartwig,81 in which ΔHdiss for the dissociation process is taken to approximate . Selected bond distances are shown in pink (Å), and relative free energies (kcal mol−1) are given with respect to intermediate B. | ||
From these results, two important conclusions can be drawn. (1) The relative free energy of intermediate X2 (22.7 kcal mol−1, Fig. 12), which would be the species attacked by the R˙ radical to form RSO2Cl, lies much higher than that of MCP1(S2) (15.2 kcal mol−1, Fig. 6). This comparison shows that this coordination-based pathway is considerably less favourable than the S(VI)/S(IV) redox cycling mechanism identified in this study. (2) The SO2Cl− anion shows a clear preference for coordination to Cu(II) through its Cl atom rather than through an oxygen atom. However, when bound through the Cl atom, SO2 is readily released, forming gaseous SO2 and complex X4. Therefore, if the reaction were to proceed via a coordination-based pathway, it would necessarily lead to R–Cl, because in this case the system is less likely to attain the Cu–OSOCl configuration that reacts with the R˙ radical to form RSO2Cl.
The overall solvation free energy values reported in this work were obtained by adding the free energy corrections from the frequency analysis to the single-point solvation free energies calculated using the SMD solvation model. Moreover, the overall free energy values were corrected by adding ΔG1atm→1M = 1.89 kcal mol−1 to account for the free-energy change associated with compressing 1 mol of an ideal gas from 1 atm to the 1 M solution phase standard state.
The activation free energies for outer-sphere single-electron transfer (SET) processes were calculated using the asymmetric Marcus theory, as recommended in the literature for obtaining more accurate barriers when the curvatures of the reactant and product potential energy surfaces differ.99,100 The solvent reorganisation contribution was obtained using the NonEq = write and NonEq = read keywords in Gaussian 16, which store and retrieve the nonequilibrium solvation effects arising from the reorganisation of the solvent cage in response to changes in the solute's electronic structure, as applied in previous studies.101–103
The respective minimum energy crossing points (MECPs), as reported in Fig. 2, were located using the code developed by Harvey and co-workers.104
To more accurately describe the electronic structure of the open-shell intermediates, multiconfigurational CASSCF calculations105 were performed for R˙(S3), F(S1), and E(S1). For F(S1) and E(S1), a (5,4) active space was employed, comprising five electrons in four orbitals corresponding to the σ and σ* orbitals of the C–S and S–Cl bonds. This selection effectively accounts for the bonding and antibonding interactions governing electron delocalization along these bonds. For R˙(S3), a (7, 7) active space was defined to include the six π electrons of the aromatic ring and the singly occupied p orbital on the β-carbon radical center, thereby encompassing the π/π* orbitals of the aromatic system and the benzylic radical orbital. This multireference treatment affords a more rigorous description of the electronic structure of these intermediates, complementing the DFT results and elucidating the extent of spin and charge delocalization. Because spin density is not rigorously defined for multiconfigurational wavefunctions, the distribution of unpaired electrons was examined through the odd-electron density (OED), computed using the Multiwfn 3.8 program.106 The OED provides a physically meaningful representation of the unpaired electron distribution, analogous to the spin density obtained from single-determinant methods.
| This journal is © The Royal Society of Chemistry 2026 |