William D.
Jobbins
a,
Bono
van IJzendoorn
a,
Inigo J.
Vitorica-Yrezabal
b,
George F. S.
Whitehead
b and
Meera
Mehta
*a
aDepartment of Chemistry, University of Manchester, Oxford Road, Manchester, M13 9PL, UK. E-mail: meera.mehta@manchester.ac.uk
bX-ray Diffraction Facility, University of Manchester, Oxford Road, Manchester, M13 9PL, UK
First published on 31st January 2023
Despite being known for decades, the solution-state molecular chemistry of heptapnictogen ([Pn7]3−; Pn = P, As) clusters is not well established. Here we study heavy element derivatives of tetrel functionalized heptapnictogen clusters towards heteroallene capture, specifically isocyanates, an isothiocyanate and CO2 are probed. Clusters (Me3Ge)3P7 (1), (Et3Ge)3P7 (2), (nBu3Sn)3P7 (3), and (Me3Si)3As7 (4) were all found to capture isocyanates between all three of their tetrel–pnictogen bonds. In the case of phenyl isocyanate insertion, tetrel coordination at the isocyanate nitrogen atoms is preferred, while in the case of p-toluenesulfonyl isocyanate insertion, tetrel coordination at oxygen is preferred. Furthermore, the reaction of (Me3Si)3P7 with CO2 gave NMR spectra consistent with the capture of the greenhouse gas. Heteroallene insertion at these clusters was also studied using density functional theory.
Historically, the molecular chemistry of these clusters has been focused on the preparation of heteroatomic or new cluster morphologies, their coordination chemistry with d- and f-block metals, and salt metathesis with group 14 electrophiles.1,6,7,10,11 The emphasis of these investigations has largely been on understanding the structure, bonding, and physical properties of the material. However, the application of [Pn7] (Pn = P, As) clusters in small molecule activation is in its early stages. (Fig. 1) In 2012 and 2014, Goicoechea and co-workers found that protonated [HPn7]2− clusters hydropnictinated carbodiimides and isocyanates.12–14 The Goicoechea group also found that the reaction of [Pn7]3− with alkynes gave 1,2,3-tripnictolides ([R2C2P3]−),15,16 whereas the reaction with CO gave the phosphaethanyloate anion ([PCO]−).17,18 Recently, we have found that boron-functionalized [P7] cages can be applied as transition-metal free catalysts in hydroboration reactions.19,20
Heteroallene insertion between labile non-cluster tetrel–pnictogen bonds has previously been reported.21–26 In 2022 we found for the first time this reactivity with Zintl-derived clusters.27 Neutral tris-silyl functionalized heptaphosphorus clusters with the general formula (R3Si)3P7 (R = Ph, Me) were found to capture and exchange isocyanates and an isothiocyanate (Fig. 1). Here, we show the broad application of this insertion, and expand the scope of functionalized heptapnictogen clusters investigated with heteroallenes to include (R3E)3Pn7 (E = Si, Ge, Sn; Pn = P, As) derivatives. Isocyanates (RN=CO) are interesting substrates for capture because of their structural relationship with the greenhouse gas carbon dioxide (OCO). In addition, this family of tetrel functionalized heptapnictogen clusters was reacted with CO2.
Efforts were also made to prepare germanium and tin-functionalized arsenic clusters (Et3Ge)3As7 and (nBu3Sn)3As7 by the reaction of the group 14 halide with [As7]3−. When Et3GeCl was reacted with [As7]3−, crystals of the new eleven-atom arsenic cluster (Et3Ge)5As11 (5) were obtained (Fig. 2). Single crystal X-ray diffraction (XRD) studies of 5 revealed an As11 core with 5 exo germanium groups. The arsenic core appears to contain an As7 unit with one basal As–As bond cleaved and further coordination of an [As4] unit to two of the bridging and one of the basal As atoms. This As11 core is distorted compared to previously reported [As11]3− and [P11]3− cores, which appear to contain four fused [Pn5] faces.31–33 An average Ge–As bond length of 2.443(3) Å was observed. Mass spectrometry (MS) studies on the bulk isolated material confirmed the presence of (Et3Ge)5As11 (5) and the targeted product (Et3Ge)3As7. To better understand how compound 5 was formed, the [As7]3− precursor was studied by mass spectrometry, which confirmed the presence of [As11]5−. The extent of this contamination was found to be batch dependent. This [As11]5− contaminant is thought to react with Et3GeCl to form 5. As only resonances from the ethyl signals can be observed by 1H and 13C{1H} NMR spectroscopy, the percentage of 5vs. (Et3Ge)3As7vs. other by-products could not be determined. The reaction of Et3GeCl with a batch of [K3(DME)x][As7] which had no [As11]5− impurity detected by MS still showed minor impurities in the 13C{1H} NMR spectrum. The reaction conditions that lead to contamination of the [K3(DME)x][As7] precursor with [As11]5− are not yet understood. Furthermore, when nBu3SnCl was reacted with [As7]3−, the unreacted tin chloride precursor could be observed by 1H, 119Sn, and 13C{1H} NMR spectroscopy, even after multiple purification efforts including distillation under dynamic vacuum and washing with solvents including toluene, ether, and pentane. Difficulties in clean formation of germanium and tin functionalized arsenic clusters precluded further studies of their reactivity.
Fig. 2 Molecular structure of 5. Anisotropic displacement ellipsoids pictured at 50% probability. Arsenic: plum; germanium: green; and carbon: white. |
Clusters 1–4 were then reacted with heteroallenes phenyl isocyanate, p-toluenesulfonyl isocyanate, and phenyl isothiocyanate. In the case of phenyl isothiocyanate, no reactions in any cases were observed, whereas reactions with phenyl isocyanate and p-toluenesulfonyl isocyanate were much more fruitful.
In the case of phenyl isocyanate, reactions of clusters 1 and 2 both gave complete conversion to the heteroallene inserted products 6 and 7, respectively (Scheme 2). Heteroallene insertion between the tetrel–pnictogen bonds of the cluster can lead to the formation of symmetric and/or asymmetric isomers, with the symmetric isomer being thermodynamically favourable (shown in Fig. 3 and discussed in more detail in ESI section 3†).27 In the case of the symmetric isomer, 3 resonances are expected in the 31P NMR spectrum for the basal, apical and bridging P atoms, whereas for the asymmetric isomer all 7 P atoms are magnetically inequivalent. The NMR spectra of both 6 and 7 showed three new resonances in the 31P NMR spectrum consistent with the exclusive formation of the symmetric isomer. In the case of compound 6, the 13C{1H} NMR resonance for the carbon bound to phosphorus could not be observed, while for compound 7 this resonance appears as a doublet at 176.0 ppm with a 1JCP of 49 Hz. Single crystal XRD studies further elucidated the solid-state structures of 6 and 7 (Fig. 4 and 5). In both structures, it was found that the Ge centre coordinates to the isocyanate nitrogen atoms. Average Ge–N bond lengths of 1.921(3) Å and 1.936(7) Å were observed for 6 and 7, while average CO double bond lengths of 1.227(5) Å and 1.231(7) Å were observed for 6 and 7, respectively. This coordination is in contrast to our previously reported reactivity between phenyl isocyanate and (Me3Si)3P7, where upon tris-insertion the Si centres coordinated to two of the isocyanates via the nitrogen atoms and one via the oxygen atom.27 This difference in coordination is presumably due to the lower oxophilicity of Ge than that of Si.34 Unfortunately, combining phenyl isocyanate with clusters 3 and 4 did not lead to reactions.
Fig. 4 Molecular structure of 6. Anisotropic displacement ellipsoids pictured at 50% probability. Phosphorus: orange; nitrogen: blue; oxygen: red; germanium: green; and carbon: white. |
Fig. 5 Molecular structure of 7. Anisotropic displacement ellipsoids pictured at 50% probability. Phosphorus: orange; nitrogen: blue; oxygen: red; germanium: green; and carbon: white. |
Next, p-toluenesulfonyl isocyanate was allowed to react with clusters 1–4 (Scheme 3). In the case of germanium functionalized phosphorus cluster 1, after 4 days 81% conversion to the inserted product 8 was observed by 31P NMR spectroscopy. The single crystal XRD data of 8 confirmed that the Me3Ge groups were coordinated to the isocyanate oxygen atoms. This coordination is in line with the previously reported reactivity of the (Me3Si)3P7 cluster with p-toluenesulfonyl isocyanate,27 where the electron-withdrawing group of the isocyanate is thought to decrease the basicity of nitrogen and favour oxygen coordination for the tetrel. The XRD data of 8 revealed an average Ge–O bond length of 1.889(2) Å and a CN bond length of 1.294(5) Å (Fig. 6).
Fig. 6 Molecular structure of 8. Anisotropic displacement ellipsoids pictured at 50% probability. Phosphorus: orange; nitrogen: blue; oxygen: red; germanium: green; and carbon: white. |
In a similar fashion, the reaction of cluster 2 with p-toluenesulfonyl isocyanate gave 55% conversion to 9. Single crystals suitable for XRD studies could not be obtained. Similarly, the reaction of 3 with p-toluenesulfonyl isocyanate gave inserted product 10 in 41% conversion after 1 week. However, when the arsenic cluster 4 was reacted with p-toluenesulfonyl isocyanate complete conversion to 11 was observed after the same amount of time. The structure of 11 was crystallographically verified and again the silicon groups were found to be coordinated to the isocyanate oxygen atoms (Fig. 7). The XRD data of 12 showed an average Si–O bond length of 1.727(3) Å and a CN bond length of 1.283(5) Å.
Fig. 7 Molecular structure of 11. Anisotropic displacement ellipsoids pictured at 50% probability. Arsenic: plum; nitrogen: blue; oxygen: red; silicon: purple; and carbon: white. |
The formation of compounds 8–10 gave 31P NMR spectra consistent with the presence of both symmetric and asymmetric isomers (Fig. 3), with 10 resonances observed in the 31P NMR spectrum. The symmetric and asymmetric isomers were identified using 31P COSY NMR experiments and as expected the symmetric isomer was favoured. In the case of compound 11, because the cluster 31P NMR handle is lost, the symmetric and asymmetric isomers cannot be detected.
For compounds 9 and 10, since crystals suitable for XRD studies could not be obtained, NMR spectroscopy, infrared (IR) spectroscopy, and density functional theory (DFT) were used to better understand their tetrel coordination modes (see below).
For compounds 8–10, the carbon bound to phosphorus could not be observed by 13C{1H} NMR spectroscopy, whereas for compound 11 this resonance appears at 185.9 ppm. We have previously observed with our (R3Si)3P7 isocyanate insertions that27 (1) the quaternary carbon bound to phosphorus is not always observed by 13C{1H} NMR spectroscopy; (2) when it is observable and there is O,O,O-coordination of the silyl groups, the imine carbon resonances appear between 176.2 and 181.3 ppm; and (3) when it is observable and there is N,N,N-coordination of the silyl groups, the carbonyl carbon signals appear between 182.5 and 216.3 ppm. Comparison of these values and that of imine carbon of 11 shows that the resonance for the imine carbon of 11 is downfield from this range, while the resonance for the carbonyl carbon of 7 is upfield from where it would be expected. These observations suggest that the electronic properties of the isocyanate and the electrophile bound to either the N or O significantly affect the 13C{1H} NMR resonance of the carbon on phosphorus. Thus, discerning this carbon as either an imine carbon or carbonyl-like from 13C{1H} NMR spectroscopy alone is difficult. However, it is worth noting that the 31P NMR spectrum of 9 is nearly identical to that of compound 8, where O,O,O tetrel coordination was confirmed by XRD studies.
Compounds 6–11 were also investigated by IR spectroscopy, in an effort to observe either the imidate [E–O–C(NR)] CN and C–O stretches or the amide [E–NR–C(O)] CO and C–N stretches. Additionally, the structures of 6–11 with both N,N,N and O,O,O tetrel coordination modes were computed at the PBE1PBE/6-311G(d,p) (the SV(p) basis set was used for Sn atoms) level of theory to obtain predicted IR stretches and thermodynamic data. Predicted and observed IR data are summarized in Table 1 and ESI section 4.† In the case of compounds 6 and 7 both imidate C–N and CO stretches could be observed, consistent with N,N,N-coordination of the tetrel and their XRD structures. However, for compounds 8–11 the observed and predicted IR data were in agreement with the presence of amide CN and C–O stretches, consistent with O,O,O-coordination of the tetrel. Furthermore, the IR data for compounds 8 and 11 are in line with the XRD structures of both compounds. However, it is important to recognize that IR data are not diagnostic in distinguishing between the presence of imidate vs. amide moieties, because first C–N and C–O stretches appear near the fingerprint region and second the expected regions for CO and CN stretches overlap.35–38
Cluster | Label | Computeda (cm−1) | Observed (cm−1) |
---|---|---|---|
a Computed at the PBE1PBE/6-311G(d,p) level of theory, and the SV(p) basis set was used for Sn atoms (10). | |||
6 | CO | 1709 | 1601 |
C–N | 1366 | 1376 | |
7 | CO | 1705 | 1570 |
C–N | 1228 | 1187 | |
8 | CN | 1580 | 1598 |
C–O | 1343 | 1300 | |
9 | CN | 1569 | 1570 |
C–O | 1438 | 1340 | |
10 | CN | 1540 | 1559 |
C–O | 1343 | 1360 | |
11 | CN | 1587 | 1587 |
C–O | 1293 | 1298 |
Next, the free energy difference between the tetrel O,O,O- and N,N,N-coordination modes was calculated for compounds 6–11 (Table 2). Because ΔG is calculated by subtracting the free energy of the N,N,N-coordination mode from that of the O,O,O-coordination mode, a positive value means that the N,N,N-coordination is thermodynamically favoured and a negative value means the opposite. In the case of compounds 6 and 7, the tetrel N,N,N-coordination is thermodynamically favoured by 54.2 kJ mol−1 and 31.4 kJ mol−1, respectively. Meanwhile, the free energy difference between the two tetrel coordination modes of compounds 8 and 10 is not significant. And in the case of compounds 9 and 11, the computed energy differences are consistent with O,O,O-coordination being thermodynamically favourable by −18.9 kJ mol−1 and −48.2 kJ mol−1, respectively.
Cluster | ΔG [O,O,O – N,N,N]a (kJ mol−1) |
---|---|
a Computed at the PBE1PBE/6-311G(d,p) level of theory, and the SV(p) basis set was used for Sn atoms. | |
(Me3Ge-PhNCO)3P7 (6) | 54.2 |
(Et3Ge-PhNCO)3P7 (7) | 31.4 |
(Me3Ge-tosylNCO)3P7 (8) | −1.1 |
(Et3Ge-tosylNCO)3P7 (9) | −18.9 |
(nBu3Sn-tosylNCO)3P7 (10) | 7.2 |
(Me3Si-tosylNCO)3As7 (11) | −48.2 |
Thus, for compound 9, based on the 31P NMR spectrum, the close structural relationship between 8 and 9, the IR data, and the computed energy differences, we postulate that the Ge groups coordinate to the three oxygens of p-toluenesulfonyl isocyanate. For compound 10, based on the IR data and that all the other tetrels coordinate to the oxygens upon p-toluenesulfonyl isocyanate insertion, including in our previous work,27 we propose that the Sn groups also coordinate to the isocyanate oxygen atoms.
Cluster | FIAa (kJ mol−1) | HIAa (kJ mol−1) |
---|---|---|
a Computed at the BP86/SV(p) level of theory. b HIA values could not be calculated as [Me3Si]+ was used as the reference. | ||
(Me3Si)3P7 | 324 | —b |
(Me3Ge)3P7 (1) | 265 | 331 |
(Et3Ge)3P7 (2) | 277 | 340 |
(nBu3Sn)3P7 (3) | 266 | 333 |
(Me3Si)3As7 (4) | 332 | —b |
(Et3Ge)3As7 | 281 | 344 |
(nBu3Sn)3As7 | 269 | 335 |
Encouraged by the insertions observed with isocyanates, the thermodynamic stability of capturing CO2 with clusters (Me3Si)3P7, 1–4, (Et3Ge)3As7 and (nBu3Sn)3As7 was investigated computationally (Table 4). The ΔG for inserting 3 equivalents of CO2 into all three tetrel–pnictogen bonds was computed at the PBE1PBE/6-311G(d,p) level of theory, and the SV(p) basis set was used for Sn atoms. It was found that only the silicon functionalized clusters (Me3Si)3P7 and (Me3Si)3As7 (4) gave CO2 inserted products that are thermodynamically favourable, by −27.1 kJ mol−1 for (Me3Si)3P7 and −32.2 kJ mol−1 for (Me3Si)3As7. However, CO2 insertions with the germanium and tin functionalized phosphorus and arsenic clusters gave products that are significantly unfavourable with ΔG values for the reactions ranging between 91.9 and 136.5 kJ mol−1.
Clusters 1–3 showed no reactivity towards CO2, which is in line with our computational findings. However, CO2 capture with 4 should yield a thermodynamically favourable product, but experimentally no reaction is observed even after several weeks. This is postulated to be because of (Me3Si)3As7 (4) being less reactive towards heteroallene insertion than (Me3Si)3P7. For example, we observed that (Me3Si)3As7 (4) reacted much slower towards p-toluenesulfonyl isocyanate than (Me3Si)3P7 (7 days vs. 1 day). Additionally, we previously reported that (Me3Si)3P7 reacted with phenyl isocyanate; however, under the same reaction conditions (Me3Si)3As7 showed no reaction even after several weeks.
Deposition Numbers 2202018 (for 5), 2202019 (for 6), 2202020 (for 7), 2202022 (for 8), and 2202021 (for 11) contain the supplementary crystallographic data for this paper.†
Footnote |
† Electronic supplementary information (ESI) available: General information, DFT information, XRD information, IR information, synthetic procedures, analytical data, and NMR spectra. CCDC 2202018–2202022. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d2dt04074c |
This journal is © The Royal Society of Chemistry 2023 |