Pnictogen-substituted tetrahedranes

Pnictogen-substituted tetrahedranes are pnictogen-containing analogues of tetrahedranes with the formula R_xC_xPn_4-x (Pn = N, P, As, Sb, Bi). Computational work has indicated that the incorporation of pnictogens to the tetrahedral core alleviates the ring strain of tetrahedrane. Although theoretical work on pnictogen-substituted tetrahedranes has existed for decades, only the phosphorus-containing species have been synthesized. These species exhibit novel reactivities, most often through ring-opening and polymerization pathways. Phosphatetrahedranes are of interest as new retrons for organophosphorus chemistry. Their strain also make them of interest in the development of energy-dense compounds.

History

The first synthetic tetrahedral molecule, tetra-tert-butyltetrahedrane (^tBu₄C₄) was reported in 1978 by Maier and coworkers following the synthesis of other Platonic solid species, like cubane and dodecahedrane. The tert-butyl substituents were used to encumber the tetrahedral core and quell the radical-mediated ring-opening of an otherwise kinetically stable but thermodynamically strained molecule via the corset effect. As of 2023, the unencumbered tetrahedrane (H₄C₄) has yet to be synthesized.

The substitution of carbons in the tetrahedral core of tetrahedrane with pnictogens to stabilize the highly strained molecule has been suggested due to the known existence of elemental tetrahedral pnictogens. Notably, white phosphorus, the most stable allotrope of phosphorus, is tetrahedral with the molecular formula P₄. Arsenic can also exist as a metastable tetrahedral allotrope, As₄, known as yellow arsenic. Furthermore, mixed tetrahedral pnictogen molecules have been synthesized, such as AsP₃ and, more recently, (PbBi₃)^-. Elements on the extreme ends of the pnictogen family have not yet been observed in a tetrahedral Pn₄ configuration, however. Nitrogen's orbitals lack diffusivity, and bismuth’s orbitals undergo minimal hybridization due to relativistic contraction.

Computational studies into mixed pnictogen-tetrel tetrahedranes have suggested that pnictogen-substituted tetrahedranes are more stable than their all tetrel counterparts decades before their first synthesis. In 1990, calculations on azatetrahedranes suggested positive correlation between the number of nitrogens in the tetrahedral core and thermodynamic stability. In the same vein, calculations done in 2010 on pnictacubanes suggested positive correlation between the number of phosphoruses in the cuboidal core and the thermodynamic stability.

All currently synthesized pnictogen-substituted tetrahedranes in the scientific literature as of 2023.

In 2019, Wolf and coworkers synthesized the first pnictogen-substituted tetrahedrane: di-tert-butyldiphosphatetrahedrane (^tBu₂C₂P₂), produced from the reaction of nickel catalyst with phosphaalkynes. Shortly thereafter, in 2020, Cummins and coworkers announced that they had synthesized tri-tert-butylmonophosphatetrahedrane (^tBu₃C₃P) from a phosphorus-containing anthracene derivative. In 2021, Cummins and coworkers published the synthesis of triphosphatetrahedrane (HCP₃), completing the set of tetrahedral molecules with carbon- and phosphorus-containing cores.

Phosphatetrahedrane Synthesis

Despite the presentation of phosphatetrahedranes as a series of incrementally changing tetrahedrane derivatives, their syntheses are vastly different.

Tri-tert-butylmonophosphatetrahedrane

The original synthesis of ^tBu₃C₃P reported by Cummins and coworkers in 2020 begins with the phosphorus derivative of anthracene. The addition of sodium hexamethyldisilazide deprotonates the phosphorus followed by triphenylborane then bonding to the phosphorus. The anionic phosphorus-containing intermediate then forms an ionic interaction with the sodium cation in solution. The addition of tri-tert-butyl cyclopropenium ion produces the thermally stable cyclopropenyl phosphine intermediate. Upon irradiation with 254 nm light in the presence of triflic acid and either tetra-n-butyl ammonium chloride or tetramethylammonium fluoride, the anthracene leaving group is forced out, leaving a halogenated phosphine. Upon addition of lithium tetramethylpiperidide with heating, ^tBu₃C₃P and lithium halide salt is generated.

Synthesis of tri-tert-butylmonophosphatetrahedrane with anthracene leaving group.

An improved version of tri-tert-butylmonophosphatetrahedrane synthesis, where the anthracene is replaced by two trimethylsilyl groups, was reported by Cummins and coworkers a year later. To replace one of the trimethylsilyl groups with a chloride, hexachloroethane is added, generating trimethylsilyl chloride and tetrachloroethylene as byproducts. Tetramethylammonium fluoride is then added to remove the remaining trimethylsilyl group and the chloride, generating the ^tBu₃C₃P. This route has a ^tBu₃C₃P yield of 33% as compared to the original route's 19%.

Di-tert-butyldiphosphatetrahedrane

In 2019, Wolf and coworkers reported the synthesis of ^tBu₂C₂P₂ through the use of a metal catalyst. Ni(IPr)(CO)₃, upon addition of 1 equivalent of tert-butylphosphaacetylene (^tBuCP), loses two carbon monoxide ligands. The addition of a second equivalent of ^tBuCP generates the 1,3-diphosphacyclobutadiene ligand, now binding with η⁴ hapticity. Density functional theory calculations into the catalytic cycle suggest that the 1,3-diphosphacyclobutadiene isomerizes into the desired tetrahedrane. Upon addition of a final ^tBuCP, (^tBu₂C₂P₂) is released and the catalytic cycle can begin again.

Synthesis of di-tert-diphosphatetrahedrane.

Triphosphatetrahedrane

Cummins and coworkers reported the synthesis of HCP₃ in 2021. Due to the similarity of HCP₃ to AsP₃, the [Nb^II(ODipp)₃(P₃)]^- previously shown to be a retron for AsP₃ was used for the synthesis of HCP₃. To add a -CH group to [P₃]^3-, bromodichloromethane undergoes halogen abstraction, leaving a carbon-centered radical. The niobium complex then undergoes P₃ transfer to yield HCP₃. The use of bromodichloromethyl trimethylsilane instead of bromodichloromethane in this process yields trimethylsilyl triphosphatetrahedrane ((Me₃S)CP₃).

Synthesis of triphosphatetrahedrane.

Reactivity

Tri-tert-butylmonophosphatetrahedrane

Lewis Acid-Induced Reactions

Addition of W(CO)₅(THF) to ^tBu₃C₃P generates a phosphorus-containing housene analogue.

The addition of 0.2 equivalents of triphenylborane in benzene can produce several cycloadducts. In the absence of exogenous reagents, ^tBu₃C₃P dimerizes into a ladderane-like compound with a P-P bond. In the presence of excess styrene or an atmosphere of ethylene, [4 + 2] cycloadditions occur to give 1-phosphabicyclo[2.2.0]hexenes.

Lewis acid-induced reactions of triphosphatetrahedranes.

Silylene Reaction

The cage opening of ^tBu₃C₃P can be induced by PhC(N^tBu)₂SiN(SiMe₃)₂ over the course of 24 hours to generate the dark red phosphasilene PhC(N^tBu)₂Si=P(^tBu₃C₃).

Monophosphatetrahedrane reaction with silylene.

Ylide Reaction

Reaction of ^tBu₃C₃P with the ylide Ph₃P=CH₂ over 48 hours and with heat induces cage opening in the same manner as the silylene reaction to generate H₂C=P(^tBuC)₃. Reaction of this product with ^tBu₃C₃P generates the symmetric product (^tBuC)₃P(C)P(^tBuC)₃.

Monophosphatetrahedrane reaction with Ph₃P=CH₂ ylide.

Formation of Phosphirane

^tBu₃C₃P is a retron for phosphirane synthesis. Upon reaction with Ni(COD)₂ (COD = cycloocta-1,5-diene) catalyst in triisopropylphosphine, cage opening occurs. Like the silylene and ylide reactions, the phosphorus bridges the (^tBuC)₃ and the alkene components. The phosphate undergoes cycloaddition with the double bond to form the phosphirane moiety. This reaction pathway has been demonstrated for styrene, ethylene, and neohexene. Furthermore, this reaction pathway is also capable of synthesizing vinyl-substituted phosphirane as evidenced by ^tBu₃C₃P and cyclohexa-1,3-diene.

Formation of phosphiranes from monophosphatetrahedrane.

Ligand Substitution

^tBu₃C₃P can be used to replace the ethylene ligand of (Ph₃P)Pt(C₂H₄) in melting THF.

Monophosphatetrahedrane substitutes an ethylene ligand.

Di-tert-butyldiphosphatetrahedrane

Dimerization Reactions

Above the melting point of ^tBu₂C₂P₂ (–32 °C), ^tBu₂C₂P₂ dimerizes into another ladderane-like structure but it is prone to decomposition. This reaction can be hampered by keeping ^tBu₂C₂P₂ under its melting point and/or by keeping the ^tBu₂C₂P₂ concentration low.

^tBu₂C₂P₂ can also be dimerized using nickel complexes to form a variety of exotic structures. ^tBu₂C₂P₂ reacted with 1 equivalent of Ni(Cp^R)(IPr) (IPr = 1,3-bis(2,6-diisopropylphenyl)imidazolin-2-ylidene, R = H, CH₃, 4-(CH₃CH₂)-C₆H₄) generates 0.5 equivalent of a tetracyclo-compound. Upon addition of another equivalent of the same nickel complex, a butterfly-like geometry is adopted, with two nickel atoms coordinated to opposite phosphorus atoms and two coordinated to adjacent phosphorus atoms on different four membered rings. This butterfly-structured compound is a dark red color. The reaction to the butterfly structure is believed to depend on kinetic access to the middle P-P bond. Bulky substituents on Cp^R kinetically hinder the P-P bond cleavage and transformation into the butterfly-structured product.

Stepwise dimerization of diphosphatetrahedrane with Ni(Cp)(IPr).

^tBu₂C₂P₂ can also be reacted with Ni(IMes)₂ (IMes = 1,3-bis(2,4,6-trimethylphenyl)imidazolin-2-ylidene) in toluene to produce 0.5 equivalents of an asymmetric compound with two hexacoordinate Ni atoms, a Ni-Ni bond, and two weak P-P interactions. This product is an intermediate for further chemistry. Heating the product at 60 °C for 3 hours causes the expulsion of a di-tert-butylacetylene and the reformation of P-P bonds. Another reaction pathway involves the addition of 3 equivalents of CO to the product, leading to the production of Ni(IMes)(CO)₃ and the ladderane-like compound described at the beginning of this section. A third reaction pathway involves the addition of hexachloroethane. This produces a 1,2-diphosphocyclobutadiene ring (vide supra) that is coordinated to both nickel atoms. This third reaction pathway also produces the ladderane analogue.

Di-tert-butyl-diphosphatetrahedrane reaction with Ni(IMes)₂ and further reactions.

Ligand Substitution

In solution with coordination complexes, ^tBu₂C₂P₂ can cause ligand substitution. Most of these reactions cause cage-opening. The reaction of ^tBu₂C₂P₂ with [K([18]crown-6)][Fe(anthracene)₂] in toluene and THF causes the expulsion of one anthracene and the cage opening of ^tBu₂C₂P₂ to form the replacement 1,2-diphosphacyclobutadiene ligand. This cage opening is due to P-C bond cleavage. Mono-ligand substitution is also observed in the reaction of ^tBu₂C₂P₂ with [(^DippBIAN)M(COD)] (Dipp = 2,6-diisopropylphenyl, BIAN = bis(arylimino)acenaphthene, M = Fe, Co). The cobalt product, upon reaction with Cy₂PCl (Cy = cyclohexyl), forms a 1,2,3-triphospholium ligand. ^tBu₂C₂P₂ can also displace the toluene in Ni(toluene)(IPr). Toluene is replaced by the aforementioned ladderane analogue in a η²-fashion.

Diphosphatetrahedrane ligand substitution of anthracene forms 1,2-diphosphacyclobutadiene analogue.

Diphosphatetrahedrane ligand substitution of cycloocta-1,5-diene forms triphosphacyclopentane analogue.

Diphosphatetrahedrane ligand substitution of toluene forms ladderane analogue.

Double-ligand substitution is seen in the reaction of ^tBu₂C₂P₂ with [Co(COD)₂][K(THF)_0.2]. The major product has its ligands doubly substituted by 1,2-diphosphacyclobutadiene. A minor product with double substitution by 1,3-diphosphacyclobutadiene ligand is also observed. However, the cobalt complex with one 1,2-diphosphacyclobutadiene ligand and one 1,3-diphosphacyclobutadiene ligand is not observed; this is likely due to steric clash between the tert-butyl substituents. The preference for 1,2-diphosphacyclobutadiene makes ^tBu₂C₂P₂ a potentially valuable retron as phosphalkynes are known to produce 1,3-diphosphacyclobutadiene.

Diphosphatetrahedrane ligand substitution of cycloocta-1,5-diene forms diphosphacyclobutadiene analogues.

Only two of ^tBu₂C₂P₂'s ligand substitution reactions are known to preserve the tetrahedral cage. Reacting (pftb)[Ag(CH₂Cl₂)₂] (pftb = Al[PFTB]^- = Al[OC(CF₃)₃]₄^-) with ^tBu₂C₂P₂ in lightless conditions leads to the generation of a disilver complex wherein each of the two ^tBu₂C₂P₂ ligates to one silver atom and each of the two ladderane analogues (vida infra) ligates to both silver atoms. By reacting Ni(CO)₄ in THF at –80 °C and in the absence of light, three ^tBu₂C₂P₂ molecules coordinate to Ni in an η²-fashion. Calculations by Intrinsic Bond Orbital (IBO) theory suggest that the coordination occurs through a 3-center-2-electron bond.

Diphosphatetrahedrane ligand substitution with [Ag(C₂H₂Cl₂)₂](pftb).

Reaction of diphosphatetrahedrane with Ni(CO)₄ yields intact diphosphatetrahedrane ligands.

Reactions with N-Heterocyclic Carbenes

^tBu₂C₂P₂ can be used a retron to form phosphirenes or phosphaalkenes with the addition of 1 equivalent or 2 equivalents of N-heterocyclic carbenes (NHC), respectively. Upon the addition of 1 equivalent of IPr, IMes, or ^MesDAC (1,3-bis(2,4,6-trimethylphenyl)diamidocarbene), ^tBu₂C₂P₂ undergoes ring opening at one phosphorus atom's P-C bonds, creating structures with a bridging P-P bond between the NHC and the phosphirene. IBO calculations and crystallographic evidence support the assignment of double bonding to the P=C bond to resultant molecule. This reaction is very slow, taking several weeks to reach completion.

Upon the addition of 2 equivalents of TMC (2,3,4,5-tetramethylimidazolin-2-ylidene) in benzene, both phosphorus atoms bond to THC. The P-P bond is broken. A double bond also forms between the two carbons of ^tBu₂C₂P₂, generating a phosphaalkene. This reaction happens significantly faster, with a reported speed of 1 hour.

Selectivity between the two reactions is suggested to be achieved by changing the steric bulk of the NHC used. A bulky NHC should prefer generating a phosphirene, whereas a smaller NHC should prefer generating a phosphaalkene.

Diphosphatetrahedrane reactivity with NHCs is controlled by sterics.

Triphosphatetrahedrane

Reaction with (dppe)Fe(Cp*)Cl

HCP₃, upon addition of [(dppe)Fe(Cp*)]Cl (dppe = 1,2-bis-(diphenylphosphino)ethane) in sodium tetraphenylborate and THF, undergoes salt metathesis and produces [(dppe)Fe(Cp*)(HCP₃)][BPh₄]. This product is crystallizable, producing purple crystals. This product is prone to decomposition back to HCP₃.

Triphosphatetrahedrane reaction with (dppe)Fe(Cp*)Cl.

Theoretical Work

Azatetrahedranes

Bonding Parameters

Ring and cage strain results in poor angular overlap of orbitals, leading to non-linear bonding. Due to the interest in tetrahedrane and their azatetrahedral analogues as highly strained molecules, Politzer and Seminario introduced the "bond deviation index" to determine the deviation of the bond path — defined as the path following maxima between nuclei — and the linear bond between the nuclei.

$\lambda ={\frac {1}{R}}{\bigg [}{\frac {1}{N}}\sum _{i}^{N}{r_{i}^{2}}{\bigg ]}^{\frac {1}{2}}$

The strain of experienced by H₄C₄ is calculated to be partially alleviated upon the substitution of carbon atoms by nitrogen atoms. The bond deviation index decreases with the number of nitrogens in azatetrahedranes from 0.114 to 0.087 from 0 to 2 nitrogen atoms. The engendered stability is countered by the propensity of the N-N bond in the strained system to escape as dinitrogen. This is evidenced by the calculated bond length: 1.59 Å, which is much higher than that of aromatic N-N bonds: 1.21-1.36 Å.

On the basis of the higher electronegativity of nitrogen than carbon, azatetrahedranes have less negative electrostatic potentials at their C-C bonds than H₄C₄, leading to greater stability against electrophilic attacks. C-C and N-C bonds both contract in azatetrahedranes. This contraction leads to smaller bond deviations and stronger bonding, although the latter is dependent on the distance from the nitrogen atom(s). Azatetrahedranes also seem to distort from an ideal tetrahedral symmetry. These trends directly correlates with the number of nitrogen atoms. Nitration also contracts the tetrahedral core. However, essentially thermodynamically neutral nitro group rotation leads to small amounts of C-N lengthening, weaking the interaction and localizing the electrons.

Ionization of all azatetrahedranes cores in the series led to cage-opening at the G4MP2 and G4 levels of theory.

Thermodynamics

Early work on azatetrahedranes has also utilized isodesmic reactions — aphysical reactions where compounds are changing but bond types are not — to understand molecular stability. The isodesmic reaction energy is thus a metric of the stabilization/destabilization relative to the starting reagents. Overall, as the number of nitrogens increase, the stability of the system increases. The isodesmic reaction energy goes from 151.8 kcal/mol to 81.8 kcal/mol for 0 to 4 nitrogen atoms. Nitration destabilizes the energy dramatically, as exemplified by 2,4-dinitro-1,3-diazatetrahedrane ((O₂NC)₂N₂) having an isodesmic reaction energy of 152.9 kcal/mol.

Due to the instability of many azatetrahedranes, isodesmic comparisons to azacyclobutadiene analogues have been used to determine which core structures are the most synthetically feasible. Alkorta, Elguero, and Rozas reported that every member of the azatetrahedrane core series is always slightly more unstable than their azacyclobutadiene analogue(s). Jursic's calculations suggest that the energetic differential between the azatetrahedrane and the azacyclobutadiene starts off large and decreases as the number of nitrogen atoms increase in the cage. Furthermore, Jursic's calculations suggest that tetraazatetrahedrane may be slightly more stable (difference of 4.4 kcal/mol at 0 K and the CBSQ level of theory) than its azacyclobutadiene analogue. Despite the instability, non-sterically hindered azatetrahedranes may still be detectable in gaseous matrices.

Phosphatetrahedranes

Intrinsic bond orbital analysis of monophosphatetrahedrane cage closing. Note that the intrinsic bond orbital corresponding P-Cl σ bond is localized on the chloride anion upon cage closing.

Bonding Parameters

Molecular graph of tri-tert-butylmonophosphatetrahedrane optimized at the B3LYP-D3/6-31G** level of theory shows significant bond deviations.

Much like azatetrahedranes, phosphatetrahedranes show bond deviation. The molecular graph of ^tBu₃C₃P shows significantly warped bond paths; the bond critical points lie far from the linear representation.

Investigations into the bonds of phosphatetrahedranes used quasi-atomic orbital analysis with isodesmic reactions. In contrast to azatetrahedranes, the lower electronegativity of phosphorus relative to carbon leads to a localization of negative charge in the carbons and the localization of positive charge in the phosphorus(es). This leads to greater occupancy in the C-C-C ring from 4.17 to 4.28 electrons with the addition of one phosphorus atom. Overall, the extent of charge transfer increases with the number of phosphorus atoms in the tetrahedral core.

Of the phosphatetrahedranes, only the triphosphatetrahedrane core did not show evidence of cage-opening upon ionization.

Thermodynamics

Depiction of the corset effect in tri-tert-butylmonophosphatetrahedrane where the attractive interaction between hydrogen atoms of neighboring tert-butyl groups are highlighted in red.

Ivanov, Bozhenko, and Boldyrev studied the energetic landscape of the phosphatetrahedrane series ((HC)_xP_4-x). Their calculations suggest that substitution of phosphorus for carbon increasingly favors the tetrahedral structure. The tetrahedral structure is the absolute minima starting for triphosphatetrahedrane, but diphosphatetrahedrane is only 2.3 kcal/mol higher in energy than the absolute minima. They attribute the stabilization of the tetrahedral structure to phosphorus' amiability towards more acute bond angles. The more diffuse orbitals of phosphorus versus carbon also favor the tetrahedral structure's σ-interactions over the planar phosphacyclobutadiene's π-interactions.

Riu, Ye, and Cummins report similar computational findings. Their calculations show decreasing strain energy with the number of phosphorus atoms in the tetrahedral cage. They also attribute the stabilization to the diffusitivity of phosphorus orbitals. They also note that the accumulation of p-character on the bond orbitals leads to greater s-character on the lone pairs.

The isolability of ^tBu₃C₃P was attributed to the controversial hydrogen-hydrogen bonds (HHB), which some chemists have argued may not exist. Each HHB of the tert-butyl network were calculated (in absence of steric repulsion) to contribute 0.7 kcal/mol of stabilization. Calculations with one of the tert-butyl substituents with a methyl, ethyl, or isopropyl group result in net repulsion due to the loss of HHBs. In total, this forms the basis of the corset effect. Non-Lewis donation of electron density from the tetrahedral core to the tert-butyl substituents also stabilizes ^tBu₃C₃P according to natural bond orbital theory. This effect was also demonstrated in silico for unsubstituted monophosphatetrahedrane.

Substitution by Heavier Congeners

Schaefer and coworkers, in light of the synthesis of ^tBu₃C₃P, ran calculations on the mono-pnictogen-substituted tetrahedrane series, represented by R₃C₃Pn (R = H, ^tBu, Pn = N, P, As, Sb, Bi). These calculations yielded a series of well correlated trends.

Consistent with a perturbation of the pnictogen residing above a C-C-C ring, the C-Pn bonds elongate from 1.493 Å to 2.289 Å, and C-Pn-C angle decreases from 58.0° to 37.1° as heavier congeners are used. This is due to the larger atomic radius of the heavier pnictogens. The H-[C-C-C plane] angle increases from 9.1° to 31.1°, which is also attributed to the diffusivity of the heavier congener's orbitals.

As noted above with the aza- and phosphatetrahedranes, the change in pnictogen electronegativity changes the interaction between the Pn atom and the C-C-C ring. The C-C-C ring becomes increasingly more negatively charged with the heavier pnictogens.

Tabulated bond angles and lengths of the mono-pnictogen-substituted tetrahedrane series calculated at the CCSD(T)/aug-cc-pwCVTZ(-PP) level of theory.

Isodesmic reactions show greater stabilization of the cage structure due to the diffusivity of the pnictogen's orbitals, although even with bismuth, the mono-pnictogen-substituted tetrahedrane is unstable. Delocalization plays a large part in the stabilization of the heavier analogues. For example, electron density is increasingly transferred from the Pn-C bonds into the Pn lone pair in the heavier congeners. These lone pairs are also noted to follow Bent's rule.

As noted above with ^tBu₃C₃P, non-Lewis interactions stabilize the tetrahedral core. These effects also become more pronounced with the heavier pnictogens. Second order perturbations suggest that the key non-Lewis interactions are C-C to C-R* and C-Pn to C-H* (i.e., cage-opening), as well as interactions to Pn-C*. The former set of interactions stabilize the tetrahedral core most when the substituent is an electron-withdrawing group (e.g., fluoride), although decreased electron density in C-C and C-Pn can facilitate cage-opening as well.

See also

• Tetrahedrane
• White Phosphorus
• Yellow Arsenic
• Prismane
• Cubane
• Cyclobutadiene