Open Access Article
Chenyang
Hu
ab,
Maren
Pink
b and
Jose M.
Goicoechea
*a
aDepartment of Chemistry, Indiana University, 800 East Kirkwood Ave., Bloomington, Indiana 47405, USA. E-mail: jgoicoec@iu.edu
bDepartment of Chemistry, University of Oxford, Chemistry Research Laboratory, 12 Mansfield Rd., Oxford, OX1 3TA, UK
First published on 10th July 2025
We describe the synthesis of an isolable phosphinidene sulfide and phosphinidene selenide (BnArNP
Ch; Ch = S, Se) by reaction of a sterically protected phosphinidene oxide (BnArNP
O) with Lawesson's or Woollins' reagents, respectively. These compounds reveal short P
E bonds consistent with a high degree of pπ–pπ double bond character. The pendant amido group on the phosphorus atom plays an important role in their stabilization through nitrogen lone pair donation into the P–Ch π* orbital, and evidence for isomerism about the N–P bond was observed for both species. Dimeric intermediates featuring mixed phosphorus(III)/phosphorus(V) centers could be isolated and structurally authenticated during the course of these studies providing some insight into the mechanism of formation of the titular compounds.
Ch; where Ch = group 16 element), are also known to exhibit carbene-like reactivity.14,15 That being said, the chemical properties of phosphinidene chalcogenides remain poorly understood on account of their inherent instability under ambient conditions. The polarity of phosphorus–chalcogen bonds, and their weak π-bonding character, makes them highly reactive, which often leads to deleterious dimerization or oligomerization reactions. Given the difficulties in isolating heavier phosphinidene chalcogenides, they have thus far only been observed at extremely low temperatures (3–15 K) in inert gas matrices (Fig. 1a).16
Two different approaches towards the chemical stabilization of phosphinidene chalcogenides have been previously explored in the chemical literature. The first of these, which we term the intermolecular approach, involves the trapping of in situ generated phosphinidene chalcogenides with chemical reagents (e.g. dienes).14,17–27 Using this approach, the formation of a transient phosphinidene chalcogenide can be inferred by the resulting products (often a heterocycle of some kind), however this does not allow for the isolation of the reactive P
Ch bond. The second strategy involves the intramolecular electronic stabilization of the P
Ch moiety by adjacent metal centres,28–31 or electron-donating R substituents.32–35 This latter strategy was successfully employed by Schmidpeter and co-workers to access the only structurally authenticated R–P
S and R–P
Se compounds reported to date (Fig. 1b). Both compounds are stabilized by electron-delocalization of the ylidyl substituent into the P
Ch π* orbital and were found to be moderately stable in solution for short periods of time.34,35
Recently, we reported the first example of an isolable phosphinidene oxide using a bulky amido supporting ligand for the stabilization of the highly reactive P(III)
O bond (BnArNP
O).36 During these studies, it became apparent that the nitrogen atom lone pair plays an important role in the stabilization of the reactive P
O moiety. We have since gone on to show that this compound can act as an electrophilic phosphinidene precursor for the synthesis of phosphorus-main group element double bonds (Fig. 1c).37 This motivated us to explore the feasibility of employing the same ligand to support heavier phosphinidene chalcogenides. Herein, we expand this strategy and report the transformation of this phosphinidene oxide to its sulfide and selenide analogues using Lawesson's and Woollins' reagents, respectively (Fig. 1d). It is worth noting that recently Tan described the isolation of heavier antimony-containing analogues, the stibinidene chalcogenides (MsFluid*Sb
Ch; Ch = S–Te; MsFluid* is a bulky hydrindacene substituent) in which the reactive Sb
Ch bonds are protected by an extremely bulky hydrindacene ligand.38
Given the structural and electronic similarity between A and carbonyl compounds, we trialled other strategies which have been shown to be effective in the chalcogenation of carbonyl-containing species. Lawesson's reagent (LR) is well established as a versatile thionation reagent. The mechanism which gives rise to its effective reactivity, involves the dissociation of the dimeric structure and the in situ generation of a dithioxo-phosphorane.39 We envisioned that the use of Lawesson's reagent could facilitate the conversion of our phosphinidene oxide, A, to its sulfide analogue, mirroring the thionation reactivity observed towards carbonyl compounds.
Reaction of 0.5 equivalents of LR with 1 equivalent of compound A (Scheme 1) gave rise to multiple products, as evidenced by in situ31P{1H} NMR spectroscopy. Addition of excess LR drove conversion to a major product which features a set of doublets at 42.5 and 47.8 ppm with a 2JP–P coupling constant of 93.2 Hz. The proton-coupled 31P NMR spectrum revealed that the resonance at 42.5 ppm splits into a doublet of triplets, with a 3JP–H coupling constant of 15.2 Hz, indicating that it originates from the LR fragment (i.e. coupling to a p-methoxy-phenyl group is observed). These data suggests that this product, 1, is a dimeric phosphinidene sulfide, incorporating a dithio-phosphorane unit derived from LR dissociation. Similar dimeric structures have been reported previously for germanium and tin, formed via the reaction of germylenes, stannylenes or stannylene sulfides with LR.40 The reaction also affords an intractable mixture of side-products which we believe to be oligomers of (MeO)C6H4P(O)(S) and (MeO)C6H4P(S)2, as previously reported in the literature.41
The structure of compound 1 was confirmed by single crystal X-ray diffraction.† As expected, the structure revealed a four-membered P2S2 ring containing one P(III) centre and one P(V) centre. Given the poor quality of the crystal and the presence of extensive positional disorder, the discussion of bond metric data is not viable, however the connectivity is consistent with a dimeric compound (see ESI, Fig. S27†). It is worth noting that a related dimeric compound DippTerP(μ-S)2P(S)DippTer, was recently reported by Hering-Junghans and co-workers (DippTer 2,6-(2,6-iPr2–C6H3)2–C6H3).42
During the purification of 1, a minor product caught our attention. Upon removal of the excess LR from the reaction mixture by washing with MeCN, two new singlet resonances appeared at 477.0 and 472.1 ppm in the 31P{1H} NMR spectrum in a 4.28
:
1 ratio, the chemical shifts of these resonances comparable favourably with the R–P
S compound reported by Schmidpeter.34 These two peaks are also similar to those observed for the phosphinidene oxide A, which also displays two 31P{1H} NMR singlet resonances arising from cis–trans isomers. Based on this spectral similarity, we assigned these peaks to phosphinidene sulfide 2. A dynamic equilibrium exists between compounds 1, 2 and LR, with the equilibrium favouring the formation of compound 1. Notably, while the removal of excess LR shifts the equilibrium toward compound 2, washing with MeCN alone was insufficient to eliminate all the residual LR.
We hypothesized that introducing an additional reagent capable of consuming LR while leaving compound 2 intact would provide the additional thermodynamic driving force needed to generate 2. LR and its analogues are known to react with Lewis bases, forming four-coordinate monomeric adducts.43,44 For this purpose, 4-dimethylaminopyridine (DMAP) was selected as an external Lewis base to neutralize LR. Satisfactorily, upon the addition of excess DMAP to compound 1, in situ31P{1H} NMR spectroscopy revealed the complete consumption of 1, and increased signal intensity for compound 2, together with a new singlet at 109.4 ppm. Independent reaction of LR with DMAP in THF confirmed that the new 31P{1H} resonance corresponds to the monomeric LR-DMAP adduct. These findings provide strong support for the success of our synthetic strategy. Subsequently, a one-pot synthesis of compound 2 was performed by sequential addition of LR and DMAP to compound A, affording 2 in excellent yields (83%) on a 50 mg scale (Scheme 1). Addition of excess LR to 2 leads to the reformation of 1 again, further confirming the equilibrium nature of this transformation.
The solid-state molecular structure of compound 2 (Fig. 2) was determined by single crystal X-ray diffraction. The crystal structure reveals a trans-arrangement of the BnArN–P
S moiety about the N–P bond, consistent with the solid-state structure of A. The P
S bond length in 2, 1.9166(14) Å, is significantly shorter than the P–S bond in Schmidpeter's ylide-substituted phosphinidene sulfide (1.981(2) Å), which the authors attribute to a significant contribution from the “dipolar resonance formula”.34 Base-stabilized phosphinidene sulfides, such as an intramolecularly stabilized species recently reported by Nikonov, also exhibit longer P–S bonds (e.g. 1.998(2) Å).45 Notably, the P
S bond in 2 is also shorter than the sum of double bond covalent radii (1.96 Å)46 and the P(V)–S bond in SPPh3 (1.9554 Å).47 It is worth noting that the short P–S bond in SPPh3 arises from hyperconjugation of the sulphur atom lone pairs into a doubly degenerate set of P–C σ* antibonding orbitals and that, consequently, such bonds cannot be referred to as double bonds. In contrast, in compound 2, the P
S bond length is attributed to a strong π-interaction between phosphorus and sulphur, reinforcing the P
S double bond character (vide infra). Notably, the nitrogen atom still retains its planar geometry
, consistent with A and previously reported amido phosphaalkenes and phosphaimines.48–50 The N–P bond length (1.674(3) Å) falls between that of single and double bond, closely resembling previously reported N–P bond length in amido phosphaalkenes and phosphaimines. These structural features indicate that the nitrogen atom lone pair is delocalized into P–S π* orbital, which aligns with the restricted N–P bond rotation observed in the aforementioned NMR studies.
To gain deeper insight into the electronic structure of 2, density functional theory (DFT) calculations were performed at the M06-2X/Def2-SVP level of theory. The optimized geometry is in good agreement with that obtained by X-ray crystallography (e.g. P–S: 1.941 Å; N–P: 1.689 Å;
). While the LUMO orbital reveals a highly localized P–S π* orbital, the P–S π orbital is diffuse and low in energy, observed as the HOMO−12. Additionally, one of the lone pairs on sulphur is distinctly characterized in the HOMO−1 (Fig. 3a). Natural bond orbital (NBO) analysis further confirmed the presence of a well-defined σ-bond and π-bond within the P
S moiety (Fig. S29 in ESI†). The π-interaction is predominantly composed of p-orbital overlap between phosphorus (96.97% p-character) and sulphur (98.14% p-character), with a significant polarization toward sulphur (70.47%). Second-order perturbation theory analysis revealed a back-donation energy of 36.05 kcal mol−1 from the nitrogen atom lone pair to the P–S π* orbital, which is comparable to that of A (41.12 kcal mol−1). The back-donation was also evidenced by the low occupancy of nitrogen lone pair (1.699 e−) and the partially filled P–S π* orbital (0.229 e−). Mayer bond indices corroborate these findings, indicating a bond order of 1.09 for N–P and 1.81 for P
S. Since the electronic structure of 2 closely resembles that of A, a similar energy barrier for N–P bond rotation should also be expected. Energy calculations at SMD-M06-2X/def2-TZVP//M06-2X/def2-SVP level revealed a rotation barrier of 20.7 kcal mol−1, closely matching that of A (21.8 kcal mol−1). The trans-isomer of 2 is thermodynamically favoured by 0.6 kcal mol−1, agreeing with the crystallographic observed trans-geometry for the N–P bond (Fig. 3b).
![]() | ||
| Fig. 3 (a) Selected Kohn–Sham molecular orbital depictions (isovalue = 0.03) for compound 2; (b) calculated N–P bond rotation energy profiles for 2 and 4. | ||
The thermodynamic profile for the equilibrium between 1, 2 and LR was further investigated using DFT calculations (see Scheme S8 in ESI† for full details). Reaction of 2 with 0.5 equivalents of LR to generate compound 1 was found to be spontaneous and thermodynamically favoured by 2.9 kcal mol−1, which aligns well with experimental observations, as compound 1 is the predominant species in the reaction mixture. DMAP-induced LR cleavage is determined to be an exothermic reaction (−6.5 kcal mol−1). By coupling these reactions, the overall transformation (formation of 2 from reaction of 1 and DMAP) favours the formation of 2, with an overall free energy change of −3.6 kcal mol−1. These computational results are in strong agreement with experimental observations, where the addition of DMAP effectively shifts the equilibrium between compounds 1, 2 and LR, to afford 2.
The successful thionation of A prompted us to explore the use of Woollins' Reagent (WR) to access a phosphinidene selenide. Woollins' reagent is the selenium analogue of Lawesson's reagent and is widely utilized in organic synthesis for the selenation of carbonyl compounds.51 Building on insights gained from the reaction between compound A and LR, we treated A with excess WR in benzene (Scheme 2). After stirring overnight, a bright-yellow solution was obtained, and in situ31P NMR spectroscopy revealed a doublet at 46.0 ppm with a 2JP–P coupling constant of 80.6 Hz, and a doublet of triplets at −46.0 ppm (2JP–P = 80.6 Hz and 3JP–H = 17.0 Hz). The similarity of this 31P NMR spectrum to that of compound 1 suggests the formation of the dimeric structure 3. Since Woollins' reagent can also be cleaved using a Lewis base,52 we envisioned that addition of DMAP can promote the formation of phosphinidene selenide 4, following an analogous mechanism to that proposed for the formation of 2. Indeed, upon the addition of excess DMAP to the reaction mixture (Scheme 2), the solution rapidly changed colour to light red. In situ31P{1H} NMR spectroscopy revealed two singlet resonances at 535.4 and 528.6 ppm in a 4.81
:
1 ratio, closely resembling the 31P{1H} NMR spectrum of compounds A and 2. Extended 31P{1H} NMR scans further revealed the presence of selenium satellites, with a 31P–77Se coupling constant of 788.1 Hz, confirming the formation of phosphinidene selenide 4. While no signals were observed in the 77Se NMR spectrum at room temperature, a broad doublet appeared at 1081.6 ppm when the NMR spectrum was recorded at −45 °C. These observations are tentatively attributed to the chemical exchange process between the cis–trans isomers at room temperature, leading to broad and weak signals. Upon cooling, the chemical exchange process is effectively frozen, resulting in enhanced spectral resolution and signal intensity in the 77Se NMR spectrum. Upon re-addition of excess Woollins' reagent, 4 can cleanly convert back to 3.
The structures of 3 and 4 were authenticated by single crystal X-ray diffraction and are comparable to those of 1 and 2, demonstrating structural consistency across the series (Fig. 4). As with its sulphur-containing analogue, phosphinidene selenide 4 adopts a trans-geometry in solid state. The P
Se bond length in 4 was determined to be 2.0662(15) Å, which again is notably shorter that the ylidyl-functionalised analogue (2.129(2) Å). This bond is also shorter than the phosphorus(V)–Se bond (2.106(1) Å) in SePPh3,53 due to presence of a strong P–Se π interaction. The planar geometry of the nitrogen atom
and the short N–P bond length (1.669(4) Å) are consistent with those observed for A and 2, indicating the back-donation from the nitrogen lone pair to P–Se σ* antibonding orbital.
DFT calculations were conducted to elucidate the electronic structure of 4. While the LUMO of 4 is predominantly composed of the P–Se π* interaction, the HOMO represents the lone pair of the selenium atom (Fig. 5). The P–Se π bond is highly delocalized and is observed in the HOMO−3 and HOMO−12. Natural bond orbital (NBO) analysis confirmed the presence of both σ- and π-bonding orbital between phosphorus and selenium (Fig. S29 in ESI†). Due to the similar electronegativity of Se (2.55) and S (2.58), the p-type π orbital remains polarized, with Se contributing 64.72% and P contributing 35.28%, which is comparable to the values determined for 2.
Second-Order Perturbation Theory analysis revealed a back-donation energy of 19.04 kcal mol−1 for nitrogen lone pair donating to the P–Se π* bond and 12.49 kcal mol−1 to σ* bond, in agreement with the crystal structure of 4. Notably, while the back-donation to π* bond is weaker than calculated for 2 (36.05 kcal mol−1), the σ* back-donation is much stronger in 4 (12.49 kcal mol−1versus 3.52 kcal mol−1 in 2). We attribute this difference to the lower electronegativity and larger atomic orbital size of Se. The P–Se π-bond in 4 is less polar, which disfavours the π* back-donation, whereas the orbital mismatch between phosphorus and selenium weakens the σ-bond, thereby enhancing the σ* back-donation. Energy calculations determined a cis–trans isomerization barrier of 21.8 kcal mol−1, with the trans-isomer being energetically favoured by 0.2 kcal mol−1, consistent with the crystallographic structure of 4. While the back-donation in 4 is weaker than that in 2, it may seem counterintuitive that the rotation barrier is higher in 4 (21.8 kcal mol−1versus 20.7 kcal mol−1 in 2). This can be rationalized by the longer P–Se bond length in 4 (2.0662(15) Å) compared to the P–S bond in 2 (1.9166(14) Å), which introduces additional steric hindrance in the transition state. This energy barrier also matches the observed 31P{1H} NMR spectrum of 4, which is similar to those observed for A and 2. Mayer bond indices suggest a bond order of 1.11 for the N–P bond and 1.75 for the P
Se bond, reinforcing the proposed bonding model.
Ch moiety and back-donation from the nitrogen atom lone pair to P–Ch π* orbital, which is also observed in A. DFT calculations further support the presence of significant back-donation, and an energetically significant rotation barrier was determined for 2 (20.7 kcal mol−1) and 4 (21.8 kcal mol−1), which is comparable to that observed in A.
:
2-cis = 4.28
:
1. HRMS (m/z): [M + H]+ calcd. For C85H109NPS, 1206.8013; found 1206.8017.
(CH3)2), 5.53 (s, 2H; NC
2Ph), 6.69–6.75 (m, 3H; ArH), 6.97–7.02 (m, 3H; ArH), 7.16 (s, overlapping with C6D6, 8H; ArH), 7.18 (d, 4JH–H = 1.7 Hz, 2H; ArH), 7.20 (d, 4JH–H = 1.7 Hz, 2H; ArH), 7.25 (t, 4JH–H = 1.7 Hz, 2H; ArH), 7.29 (d, 3JH–H = 7.7 Hz, 2H; ArH), 7.53–7.58 (m, 3H; ArH), 8.42 (dd, 3JH–P = 17.0 Hz, 3JH–H = 7.2 Hz, 2H; ArH). (Accurate integration was not possible as, in solution, 3 is in equilibrium with 4 and WR). 13C{1H} NMR (125 MHz, C6D6): δ (ppm) 24.40, 24.44, 24.68, 25.04, 25.32 (s; CH(CH3)2), 30.77, 30.84, 30.90, 31.05, 31.11, 34.84, 34.88 (s; CH(CH3)2), 55.24 (br; NCH2Ph), 120.71, 120.82, 120.96, 128.59, 129.63, 131.46, 131.56, 131.62, 131.82, 132.39, 136.26, 136.95, 137.05, 137.17, 137.38, 138.87 (d, 1JC–P = 25.3 Hz), 140.13, 140.88, 141.40, 141.54, 141.90, 142.40 (d, 2JC–P = 6.2 Hz) 143.84 (br), 146.79, 146.91, 146.96, 147.05, 147.10, 148.56, 148.63 (s; ArC). 31P NMR (202 MHz, C6D6): δ (ppm) −46.0 (dt, 2JP–P = 80.6 Hz, 3JP–H = 17.0 Hz; PSeP(Se)Ph), 46.0 (d, 2JP–P = 80.6 Hz, PSeP(Se)Ph). 31P{1H} NMR (202 MHz, C6D6): δ (ppm) −45.95 (d, 2JP–P = 80.6 Hz; PSeP(Se)Ph), 46.04 (d, 2JP–P = 80.6 Hz; PSeP(Se)Ph). HRMS (m/z): [M + H]+ calcd. For C91H114NP2Se3, 1520.5925; found 1520.5954.
:
4-cis = 4.81
:
1. 77Se NMR (95.4 MHz, d8-tol, −45 °C): δ (ppm) 1081.62 (d, 1JP–Se = 875.0 Hz; 4-trans). HRMS (m/z): [M + H]+ calcd. For C85H109NPSe, 1254.7457; found 1254.7473.
Footnote |
| † Electronic supplementary information (ESI) available: Additional experimental data, spectra, computational information, and crystallographic data can be found in the ESI. Computationally optimized structures (.xyz) are also provided. CCDC 2444805 (for 1), 2444806 (for 2), 2444807 (for 3) and 2444808 (for 4). For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d5sc04352b |
| This journal is © The Royal Society of Chemistry 2025 |