As the simplest alkyne, acetylene (C
2H
2) is an excellent model for the research of PIP. The addition polymerization of acetylene at AP was widely investigated in the last century and Nobel Prize was awarded to the research on conductive polyacetylene (PA). The PIP of acetylene has completely different reaction routes to those under AP. Under HP, acetylene solidifies into crystal and starts to polymerize at 3.5 GPa at room temperature, with the main product confirmed as
trans-PA
[10][11][12][13], whereas, at 77 K, solid acetylene polymerizes at 12 GPa, consequently forming
cis-PA
[14]. Since it is well known that the
trans-PA is more stable, authors need to clarify the reaction mechanism of this selective PIP of acetylene. Authors used Paris-Edinburgh Press (PE press) as a HP generation apparatus to avoid the photoactivation and facilitates the in situ neutron diffraction
[15][16][17]. The crystal structure of deuterated acetylene at the critical reaction pressure of 5.7 GPa was determined by in situ neutron diffraction, from which authors found the possible reaction route to produce
cis-PA
[18]. As shown
Figure 1a (left), the previous proposed bonding route along a+b or a−b direction will result in
trans-PA (grey line route)
[12][13], while there is another bonding route along the a+c and a−c direction with D
C of 3.1 Å (yellow line), which will promote the formation of
cis-polymer. The fact that anisotropic thermal factors of carbon and hydrogen (deuterium) reach maximum in the a−c plane also supports this bonding route. On the other hand, in the metadynamics simulation, it also shows that the reaction also takes place in the a+c direction, and the formation of
cis-PA is reproduced (
Figure 1a (right)). Additionally, the simulation predicted the production of graphane in a further step of PIP, which was demonstrated by authors' subsequential experiment. This is a typical example that shows the sp-sp
2-sp
3 sequence in the PIP, and also reveals the neighbor-first topochemical rule in understanding the PIP process (
Figure 1b).
Figure 1. (
a) (left) Crystal structure of acetylene at 5.7 GPa. The brown and light blue ellipsoids are carbon and deuterium atoms, respectively. The yellow lines indicate the forming of
cis-PA while the grey for
trans-PA. (right) Metadynamics simulation results of acetylene under 15 GPa. Reprinted with permission from ref.
[18]. Copyright 2017 Angewandte Chemie International Edition; (
b) PIP of acetylene, via
cis-PA to graphane. Hydrogen is omitted for clarity.
The acetylide anions such as C≡C
2− and HC≡C
− were also predicted to polymerize into conjugated extended carbon framework with excellent electrical conductivity under HP
[19][20][21][22][23][24]. As the simplest acetylides, CaC
2 and Li
2C
2 were investigated under HP
[22][25][26][27][28], and authors firstly evidenced their PIP experimentally, which resulted in a remarkable conductivity enhancement,10
7-fold for CaC
2, and 10
9-fold for Li
2C
2 [26][27]. After releasing to the atmospheric pressure, due to the extended covalent conjugated C-C bonding, the high conductivity of the recovered products maintains.
Authors also monitored the crystal structure variation under HP by in situ diffraction and complementary theoretical optimization. CaC
2 crystallizes in a tetragonal phase at ambient and low pressure, and distorts into a monoclinic phase above ~10 GPa. This accelerates the approaching of neighbored C
22−, and authors estimate that the D
C is ~2.9 Å at around 20 GPa
[26]. Above 20 GPa at room temperature, diffraction cannot provide enough structural details since CaC
2 polymerizes into amorphous phase. High resolution gas chromatography mass spectrometry (HR GC-MS) analysis of the hydrated products of the recovered sample provided solid evidence of oligomerization of acetylide. Advanced hydrocarbons were observed besides acetylene. Most of components have the C:H ratio of around 1:1 (
Figure 2a), which indicates no obvious disproportionation in the PIP. The slight deviation from 1:1 is likely due to the termination of the free radicals, which was discussed in detail in the work of NaC
2H
[29]. More importantly, benzene (C
6H
6) was identified accordingly to HR-GC-MS data. This consists with the reaction sequence of forming chain-belt-sheet as predicted by metadynamics simulation (
Figure 2b).
Figure 2. (
a) Total ion chromatograms (TIC) of the product of CaC
2 (synthesized by PE press) recovered from 26 GPa and CaC
2 raw material. (The hydrolysis of the polymerized C
mx− anions in the recovered sample produced C
2nH
2n according to the chemical reaction equation: Ca
nC
2n + 2nH
2O = nCa(OH)
2 + C
2nH
2n); (
b) Simulated chain and ribbon model of CaC
2 at 30 GPa by metadynamics
[26].
By comparing the experimental result of Li
2C
2 under high pressure with the theoretical calculated IR spectra, the carbon-ribbon structure in the Li
2C
2 polymeric product was identified, and many advanced hydrocarbons in the hydrated products of the recovered sample were detected by GC-MS, which also supports the addition PIP of C
22−, and the ribbon structures (
Figure 3a)
[27]. However, different from CaC
2, there were a portion of products with the C:H ratio significantly deviating from 1:1, which revealed the occurrence of the disproportionation reaction. Compared with the Ca
2+ in CaC
2, the Li
+ cation is more active to diffuse. This leads to disproportion and composition segregation under compression, promoting the formation of Li-rich and C-rich phases. In HPHT experiments, crystalline phase Li
3C
4 with zigzag carbon-ribbon structure (lithium polyacenide) and LiC
2 with intercalation structure (lithium graphenide) were indeed obtained (
Figure 3b)
[28]. Authors also predicted a series of lithium polycarbides with a stoichiometric ratio of Li
n+1C
2n (
Figure 3c), which share the common feature of zigzag nano-ribbon structure with various width, and LiC
2 (n = ∞) and Li
3C
4 (n = 2) are the two terminal compounds (
Figure 3b)
[28].
Figure 3. (
a) The IR spectra of the Li
2C
2 (experiment) and the predicted Li-C phases (simulation) under external pressure (above); The relative area of the deconvoluted peaks of the hydrocarbons in the hydrolyzed product (below) in the GC-MS results. The components with a score lower than 90 are counted in “others”. Reprinted with permission from ref.
[27]. Copyright 2017 American Chemical Society; (
b) The crystal structures of Li
3C
4 (left) and LiC
2 (right); (
c) The predicted crystal structures of Li
n+1C
2n (n = 1, 2, …). The black spheres: C. The purple spheres: Li. Reprinted with permission from ref.
[28]. Copyright 2018 American Chemical Society.
As mentioned, acetylene normally reacts at around 5.0 GPa
[18], while the polymerization pressures of disubstituted metal acetylides are mostly over 20 GPa (~22 GPa for CaC
2, and 35 GPa for Li
2C
2 [26][27]) due to the repulsion of the electrostatic interaction and spatial blocking of the cations. Authors also studied the monosubstituted acetylene sodium (NaC
2H) with only one charge on the carbon group
[29]. By the method of in situ spectroscopy and X-ray diffraction (XRD) technology, PIP of NaC
2H was proved beginning from 14 GPa, which is between that of acetylene and acetylides. According to the crystal structure of NaC
2H at 11 GPa (the maximum pressure at which authors can still have usable diffraction data,
Figure 4a (left), the D
C is around 3.0 Å (d
3), and by extrapolation of the d-P curve (
Figure 4a (right)) to 14 GPa, the D
C was around 2.9 Å. This is similar to that of CaC
2 [26]. The specific reaction is shown in
Figure 4b. This is a free radical addition process. Due to the produced species needing to absorb or lose at most two hydrogen atoms for pairing their two radicals, the hydrolysis products always have a C:H ratio within C
xH
x±2 (
Figure 4b).
Figure 4. (
a) Crystal structure of NaC
2H and the evolution of d
1, d
2, and d
3 as a function of pressure. d
1, d
2, and d
3 represent different intermolecular C···C distance; (
b) the reaction routes of the NaC
2H dimerization and the following H-transfer process. Reprinted with permission from ref.
[29]. Copyright 2019 American Chemical Society.
The PIP of C≡C groups should be the simplest PIP reactions. The reactions are likely to follow a free radical route. Though the HP environment prevents us to trap or detect the radicals, the hydrated products of the recovered polyacetylide sample CnHn+2, and CnHn−2 etc., provides indirect evidence. The reactions seem to follow the neighbor-first rule, which is a feature of topochemical route. Regardless of the reaction pressure, the bonding share approximately the same DC, ~3.0 Å, which is a good reference for other unsaturated groups. The hybridization of carbon in PIP follows a sp-sp2-sp3 sequence. The PIP can proceed to the end, or stop at half-way, depending on the reactivity of the unsaturated groups and the pressure applied, which provides more controlling factors for synthesizing advanced materials.
2.2. Aromatic Compounds
As important building blocks for PIP, aromatic compounds also attract continuously attentions. Since the 1980s, benzene (C
6H
6), as the proto model of aromatic systems, has been systematically explored in the high-pressure research
[30][31][32][33][34][35][36][37][38][39][40]. In recent years, through optimizing the synthetic conditions (a slow compression/decompression rate), one-dimensional sp
3 carbon nanothreads were identified from the sample recovered from compressing C
6H
6 at room temperature
[41][42]. Compared with the traditional sp
2 nanotube, higher strength and stiffness were expected for the diamond-type nanothreads. Similar results were observed in other aromatics cases, such as pyridine (C
5H
5N)
[43][44], thiophene (C
4H
4S)
[45], furan (C
4H
4O)
[46], and aniline (C
6H
5NH
2)
[47][48]. Several reaction routes, including [4+2] addition and 1-1′ addition/coupling, were investigated for benzene, and many models were proposed for the nanothreads
[49][50][51]. However, probably due to the diversity of the local structures in the solid sample before reaction, there is still existing considerable confusion about the reaction pathway and mechanism.
To investigate a model compound with well-defined ordered structure (“monodispersed” local structure), authors explored the 1:1 benzene (C
6H
6) and hexafluorobenzene (C
6F
6) cocrystal
[52][53]. It has a strong electrostatic attraction between the two components, which leads to the C
6H
6-C
6F
6 rings being alternately stacked in column regularly
[52]. After releasing from 20 GPa, PIP consequently resulted in a layered structure of H/F substituted graphane. Before reaction occurrence, the structure maintained the column stacking, the column was remarkably titled, and, in fact, the molecules formed closed packed layers (
Figure 5a). The D
C between C
6H
6 and C
6F
6 at 20 GPa were decreased to approximately 2.8 Å at room temperature (without thermal correction), reaching the critical reaction distance of previous benzene molecules (2.8~3.0 Å
[54]). However, since there are several intermolecular C-C distances closed to D
C, a [4+2] reaction is probably preferred according to the neighbor-first rule, but is not exclusive (
Figure 5b).
Figure 5. (
a) Crystal structure of C
6D
6-C
6F
6 before reaction (20 GPa): (left) viewed perpendicular to (102) and (right) approximately along. Reprinted with permission from ref.
[55]. Copyright 2020 Nature; (
b) Selected elemental reactions of the CHCF cocrystal under high pressure; (
c) A histogram of the comparison between the relative abundance of major components in different high-pressure samples from GC-MS results. Reprinted with permission from ref.
[53]. Copyright 2019 Angewandte Chemie International Edition.
Then, authors used HR GC-MS to investigate the reaction intermediates to find the primary reaction models in PIP. Taking the standard mass library and the retention time of commercial samples as reference, authors concluded that the elemental reactions include Diels-Alder, retro-Diels Alder, and 1-1′ coupling reactions (Figure 5b). By comparing the compositional contents of products recovered from different pressures, authors found that product content from the [4+2] Diels-Alder addition route was increased under higher pressure (Figure 5c), rather than the coupling reaction, which is also supported by spectroscopic analysis. Unexpectedly, authors did not observe the 1-D nanothread reported in the PIP of benzene, etc., but found an ordered 2-D H/F substituted graphane structure.
After authors' work, more similar aromatic cocrystal systems were recently reported, including two polycylic 1:1 arene-perfluoroarene cocrystals, naphthalene/octafluorona-phthalene (NOFN), and anthrancene/octafluoronaphthalene (AOFN)
[56][57]. Owing to the strong intermolecular reactions between the larger conjugation, their π-π stacked column structure was more stable under compression, and the molecules only tilted slightly with respect to the stacking direction. The D
C was along the interplanar centroid-centroid line of polycylic rings, determined as around 2.8 Å for NOFN and 2.7 Å for AOFN, and similar to that of reported aromatics. In contrast, the PIP products are one-dimensional polymers for both NOFN and AOFN.
Besides benzene-based derivates, heterocyclic compounds, especially nitrogen-containing cyclic molecules, also attract great scientific attentions
[58][59][60]. Authors investigated the structural evolution and chemical reaction processes of 1H-tetrazole (H
2CN
4)
[60], a nitrogen-rich compound with a high nitrogen content of 80%, up to 100 GPa for the purpose of creating a novel polynitrogen material with high energy density. Through the in situ synchrotron XRD, an irreversible reaction was demonstrated in the range of 60 to 100 GPa, which was also evidenced by the strong fluorescence signal detected in the recovered sample in Raman spectroscopy. Intermolecular bonding between the C and N atoms from neighbored molecules, respectively, was proposed, according to theoretical investigation (
Figure 6a). However, the bonding of nitrogen seems much harder than that of carbon. No chemical reaction proceeds even when the closest intermolecular C…N reaches 2.52 Å at 48.5 GPa, which was much smaller than the D
C of unsaturated carbon groups (
Figure 6b). The consequence obviously shows that the D
C values are highly elemental dependent, probably due to the different electronegativity and/or bonding formation energy.
Figure 6. (
a) Polymeric structure of 1H-tetrazole under 100 GPa, optimized from the structure obtained from the XRD result at 100.3 GPa; (
b) Closet intermolecular C…N and N…N distance as a function of pressure. Reprinted with permission from ref.
[60]. Copyright 2021 Royal Society of Chemistry (RSC) Publishing.
From authors' research, authors can demonstrate that the PIP of aromatics has various elemental reactions, which can coexist and compete mutually. The intermolecular distances are not far from each other, so the neighbor-first rule may not dominate. By modifying the synthetic conditions, the product may change, which brings great challenges, but also great opportunities.
2.3. Alkynlphenyl Systems
Generally, the PIP of alkynes proceeds at a relatively lower pressure, while the products often lack long-range ordering at room temperature owing to the thermal effect. In contrast, for the aromatics, the products often have better ordering, but the pressures required are much higher. To synthesize ordered polymeric materials via topochemical PIP at relatively low pressure, authors investigated the PIP of 1,4-diphenylbutadiyne (DPB), aiming to take the advantage of both aryl and ethynyl
[61]. Authors obtained crystalline armchair graphitic nanoribbons, and, importantly, revealed a new PIP route, dehydro-Diels−Alder (DDA) reaction, which conforms exactly to the topochemical strategy.
Figure 7a displays the in situ XRD patterns under HP. Before chemical reactions, DPB maintained the low-pressure structure (monoclinic phase). The crystal structure at the reaction threshold pressure of 10 GPa showed that the intermolecular distance 3.22 Å (between the alkynyl and phenyl group) and 3.24 Å (between phenyl groups) were both much smaller than those between diynes, where the shortest distance is 4.29 Å, and is still greater than the commonly required distance (4 Å) for 1,4-addition polymerization of diynes at ambient conditions (
Figure 7b). It is the aryl groups that provide a large steric hindrance to prevent the direct addition reaction of triple bonds from approaching and bonding, and hence avoid the possible 1,4-addition product. In such a local structure, the new reaction path, [4+2] DDA reaction between phenyl and alkyne groups, is activated (
Figure 7c), and results in the formation of a graphitic ribbon structure.
Figure 7. (
a) Selected in situ XRD of DPB under high pressure. The peaks from the PIP-DPB are marked by the asterisks; (
b) Crystal structure of DPB optimized by density functional theory (DFT) calculation based on the experimental lattice under 10 GPa. (
c) Proposed DDA reaction. Reprinted with permission from ref.
[61]. Copyright 2020 American Chemical Society.
Subsequently, a similar DDA reaction was also observed for 1,3,5-triethynylbenzene (TEB), where three ethynyl groups were located in the meta-position of a phenyl ring
[62]. Under compression, a one-dimensional nanoribbon structure containing sp, sp
2, and sp
3 carbons was obtained (
Figure 8a). Through in situ spectroscopic analysis, authors found that the phenyl moiety and the ethynyl groups reacted synergistically at a low pressure of 4 GPa. As displayed in
Figure 8b, the crystal structure at 3.6 GPa determined by in situ XRD revealed that the D
C for PIP was around 3.3 Å, longer than that of aromatics (2.8 Å) and even alkynes (3.1 Å), but close to the distance of the forementioned DPB case (3.24 Å). Accompanied by the structure by Rietveld refinement and NMR experimental data, the theoretical simulation of PIP process proposed reasonable product models, most of which were formed via the intermolecular [4+2] DDA reaction route between the ethynylphenyl groups and the phenyl ring (
Figure 8c). Although the original location of the involved ethynyl groups is relatively far from the benzene ring (d
2 = 3.89 Å and d
7 = 4.10 Å), the tendency of forming six-member rings is still favored, which is supposedly due to the flexibility of the terminal ethynyl, in contrast to the restrained alkynyl moieties in DPB. Moreover, it also could reasonably explain the remarkable decrease in the reaction pressure.
Figure 8. (
a) High-resolution transmission electron microscope (TEM) of the recovered sample PE5, which were synthesized at 5 GPa using PE Press; (
b) Crystal structure of TEB optimized at 3.6 GPa. The close intermolecular C⋯C distances are d
1 (3.39 Å), d
2 (3.89 Å), d
3,4,5 (3.30 Å), d
6 (3.33 Å) and d
7 (4.10 Å); (
c) Models of ribbons (3-1, 3-1-H, 3-2, 3-2-H, 3-4, 3-4-H, 3-5 and 3-5-H). New bonds are in blue, and double bonds are in red. Reprinted with permission from ref.
[62]. Copyright 2021 American Chemical Society.
In summary of the addition polymerization via PIP, authors can conclude the critical reaction distance rule for the PIP of specific unsaturated functional groups. When approaching a certain distance (~3.0 Å), unsaturated atoms generally tend to bond. The critical distance rule greatly helps us to design suitable precursors for synthesizing novel carbon materials via PIP. As presented Table 1, there are still several factors influencing the critical distance, including the functional groups, reaction element types, and electron charge numbers. The DC of pure aromatics (2.8 Å) is relatively shorter than that of alkynes (3.1 Å), which means that more compression is required to initiate the reaction. Aromatics always have extra stability but worse reactivity. When the ethynyl groups are introduced into a benzene ring, the delocalized conjugation is expanded spatially, and the benzene ring is “activated” by the ethynyl. When the local structure allows, a synergistical reaction [4+2] is preferred with longer DC. It seems that the synergistical bonding naturally has longer DC, since the two bonding can promote each other. In addition, comparing DPB with TEB, the latter molecule has a terminal ethynyl group, and the DC is accordingly increased, probably since the terminal ethynyl has more dynamic freedom, while for the metal carbide, as the charge of C22− caused strong electrostatic repulsion, it is harder to be activated, and a smaller DC is required for triggering reactions.
Table 1. Intermolecular distance threshold for PIP of unsaturated compounds.
Compounds |
DC (Å) |
Reaction at Initiation |
Pressure (GPa) |
Acetylene (C2H2) |
3.1 [18] |
free radical reaction |
5.7 |
Sodium Monoacetylene (NaHC2) |
2.9 [29] |
11 |
Calcium Acetylene (CaC2) |
2.9 [26] |
20 |
Lithium Acetylene (Li2C2) |
- |
33 |
Benzene (C6H6) |
2.8 [54] |
Diels-Alder/1-1′ coupling, etc. reactions |
18 |
Benzene-hexafluorobenzene (C6H6-C6F6) |
2.8 [53] |
20 |
1,4-diphenylbutadiyne (DPB) |
3.2 [61] |
Diels-Alder reaction |
10 |
1,3,5-Triethynylbenzene (TEB) |
3.3 [62] |
4 |