Identification of new pillared-layered carbon nitride materials at high pressure.

The compression of the layered carbon nitride C6N9H3 · HCl was studied experimentally and with density functional theory (DFT) methods. This material has a polytriazine imide structure with Cl(-) ions contained within C12N12 voids in the layers. The data indicate the onset of layer buckling accompanied by movement of the Cl(-) ions out of the planes beginning above 10-20 GPa followed by an abrupt change in the diffraction pattern and c axis spacing associated with formation of a new interlayer bonded phase. The transition pressure is calculated to be 47 GPa for the ideal structures. The new material has mixed sp(2)-sp(3) hybridization among the C and N atoms and it provides the first example of a pillared-layered carbon nitride material that combines the functional properties of the graphitic-like form with improved mechanical strength. Similar behavior is predicted to occur for Cl-free structures at lower pressures.

High pressure investigations of carbon nitride materials have been stimulated by theoretical predictions that dense sp3-bonded C3N4 phases would display low compressibility and high hardness values12. Despite many synthesis attempts including use of high pressure the formation of crystalline C3N4 polymorphs remains an elusive goal. However an sp3-bonded carbon nitride imide C2N3H with a defective wurtzite structure has been produced in laser-heated diamond anvil cell (DAC) experiments and this material was recoverable to ambient conditions34. Low density sp2 bonded polymeric or layered graphitic carbon nitride materials (g-CNMs) are also well known5 and these have attracted attention as metal-free redox catalysts, photocatalysts and electroceramic materials67. They also provide precursors for investigation of pressure-induced transformations into dense sp3 bonded phases. We investigated the room temperature compression behavior of one such well-characterized g-CNM of composition C6N9H3·HCl8 to 70 GPa using synchrotron X-ray diffraction experiments in a diamond anvil cell (DAC) combined with density functional theory (DFT) calculations carried out to 100 GPa. The combined experimental and theoretical results indicate the initial onset of layer buckling and movement of the Cl− ions out of their sites within the planes followed by a phase transformation into a structure containing interlayer C-N bonds between sp3 hybridized atoms. The new material constitutes a new example of a pillared-layered gCNM with mixed sp2–sp3 bonding. Related bonding changes have been recorded in graphitic C and BN materials at high pressure91011121314.

A large class of oligomeric, polymeric and graphitic n class="Chemical">carbon nitride materials are determined to have structures based on heptazine (tri-s-triazine, C6N7) rings linked by -N = or -NH- groups15. A second structural motif based on polymerized triazine (C3N3) units is found among graphitic nanocrystalline C3N4 prepared by chemical vapor deposition16. Triazine units are also the building blocks in bulk crystalline compounds synthesized from precursors including melamine (C3N6H6) and cyanuric chloride (C3N3Cl3) under high-P,T conditions (0.5–3 GPa, 300–500°C). Further studies have obtained crystalline graphitic polytriazine imides (PTI) with related structures by reactions in molten salt mixtures (e.g. KCl-LiCl eutectics)517. The s-triazine rings linked by bridging -N = or -NH- units lead to large voids within the graphitic layers. In the case of C6N9H3·HCl studied here the Cl− ions occupy the centers of the voids at ambient pressure. The layers are stacked to provide an alternating pattern of triazine and C12N12 rings along the c axis (Fig. 1).

Figure 1

Structural features of graphitic C6N9H3·HCl at ambient pressure.

(left) Top view of one plane of the structure showing the triazine (C3N3) rings connected by -NH- groups to form large C12N12 voids that are occupied by Cl− ions. The accompanying extra H+ ion is attached to one of six possible N positions from the triazine units surrounding the large ring616. The H atoms have been omitted for clarity. (right) View down c of two adjacent layers of the structure showing the ABAB stacking sequence that places one triazine ring above the C12N12 void in successive planes. The H and Cl species are not shown for clarity.

Results

Experimental findings

The starting compound g-n class="Chemical">C6N9H3·HCl is descibed by a hexagonal unit cell with P63/m symmetry (Fig. 1)818. The X-ray patterns at low pressure (Fig. 2) are dominated by the (002) interlayer reflection at 2θ ~ 8° that is indicative of a layered graphite-like structure with an ABAB stacking along the c axis. The polytriazine imide layers contain C12N12H3 voids hosting Cl− ions derived from the synthesis reaction between melamine and cyanuric chloride, and additional H+ ions are bound to the N atoms of the triazine rings818. Only one of six available N sites is protonated in this way so that the P63/m space group used to analyze the data represents a spatially averaged solution. The ABAB stacking of the graphitic layers creates a structure in which one half of the triazine rings within each layer is positioned above or below C12N12 voids in adjacent layers, while the others overlap triazine rings in the layers above and below (Fig. 1).

Figure 2

Angle dispersive synchrotron X-ray diffraction data for C6N9H3·HCl obtained up to P = 70 GPa (λo = 0.4441 Å).

Peaks of interest are labelled with hkl designations according to the P6m space group. (a) Selected patterns showing the onset of peak broadening for the layered phase above 10 GPa followed by the rapid emergence of the new pattern above 40 GPa. (b) Variation of d spacings measured for the principal peak maxima assigned to the low pressure graphitic phase as a function of pressure.

The prominent (002) interlayer reflection shifts rapidly to smaller d spacing (2θ ~ 10°) with increasing pressure (Fig. 2). Above 10 GPa tn class="Chemical">he diffraction peaks become significantly broadened but the general features of the graphitic material are still recognizable up to P ~ 36 GPa. However at P = 40 GPa the nature of the pattern has changed to become dominated by a main broad asymmetric peak near 13° 2θ. This result signals a structural change into a new high-density form (Fig. 2). After the transformation has occurred the peak positions vary little with continued compression up to 70 GPa indicating that the high pressure structure is significantly less compressible than the graphitic layered phase. The substantial peak broadening that occurred as a function of pressure did not permit us to carry out detailed refinement of the diffraction data. Instead we investigated the nature of the structural changes at high pressure using density functional theory (DFT) calculations, that enabled us to gain understanding of the pressure induced transformation by comparing the predicted diffraction profiles with experiment (Fig. 3).

Figure 3

Comparison of experimental (black line) with the DFT predicted (ILB0 phase, red line) X-ray patterns for the graphitic phase at various pressures.

At 40 GPa the blue pattern corresponds to dense forms with interlayer bonding (ILB1) as described in the text.

DFT calculations

In a first series of DFT studies we geometry-optimized tn class="Chemical">he structure of g-C6N9H3·HCl between P = 0–100 GPa in 10 GPa steps, starting each optimization from the ambient pressure lattice spacings obtained from experiment (a,b = 8.42 Å; c = 6.34 Å). We generated a proton ordered version of the 40-atom crystallographic unit cell within space group P1, removing all local symmetry constraints to allow the atoms to relax both within and away from the graphitic planes. This differs from our experimental data analysis that assumed a spatially averaged P63/m unit cell818. The proton-ordered structure exhibits a small shear distortion between adjacent planes along (010) to maximize electrostatic interactions between the H+ and Cl− ions within layer voids. This leads to a small separation between a and b cell parameters that increases slightly with compression but that remains insignificant until above 20 GPa (Fig. 4). At P = 0 GPa the optimized lattice parameters were found to be a,b = 8.37 Å; c = 6.23 Å, slightly smaller than those found experimentally (8.42, 6.34Å) as is typical for LDA results. During compression there is a rapid decrease in the c-axis parameter consistent with experiment: between 0–40 GPa there is a 17% reduction in c, compared with 3.7% and 5.8% in the a and b lattice parameters (Fig. 4). Between 60–70 GPa a sudden collapse of the c-axis occurred in these initial compression experiments, and then this parameter remained approximately constant up to 100 GPa (Fig. 4).

Figure 4

Changes in the DFT calculated lattice parameters for a proton ordered version of the g-C6N9H3·HCl structure as a function of pressure.

These initial results showed that a dramatic change had occurred in tn class="Chemical">he structure as the c axis dimension was reduced by compression, implying a separation between the initially planar layers of approximately 2.0 Å, without considering the reduction in distances due to layer buckling. Such a small interlayer separation could lead to the onset of C-N interlayer bonding and we investigated this possibility by analyzing the fully relaxed structures as a function of densification. The results led us to identify the structure optimised at high pressure as a new phase of the material, characterized by the occurrence of bonding between a subset of sp3 hybridized C and N atoms in adjacent layers (Fig. 5). This new phase is termed the ILB1 structure to distinguish it from the densified graphite-like phase with buckled but independent layers obtained at lower compression (ILB0). The C-N bonding occurs between atoms associated with triazine rings that in the ambient pressure phase overlapped each other in adjacent layers; the others located above and below the C12N12 voids remain sp2-hybridised giving rise to a unique new pillared-layered structure that contains a mixture of sp2- and sp3-bonded C and N atoms (Fig. 5).

Figure 5

Left: The ILB1 phase of C6N9H3·HCl obtained at high pressure from DFT calculations.In the real structure studied experimentally the Cl− displacements above and below the layers can occur in a random or locally ordered fashion. Right: picture highlighting the formation of pillars with three-dimensional connectivity between the carbon nitride layers produced by interlayer bonding within a mixed sp2–sp3 bonded state; the pillars are connected within the planes by relatively undistorted triazine rings.

Once we had identified the occurrence and connectivity of tn class="Chemical">he new high-pressure phase, we examined the relative enthalpies of the ILB0 and ILB1 polymorphs over a wide pressure range, starting the studies using the compressed but still graphite-like (ILB0) structure obtained at 60 GPa and then gradually increasing P, compared with results during gradual decompression of the ILB1 phase obtained at 70 GPa (Fig. 6). Each calculation was initiated from the geometry-optimized structure at the immediately adjacent pressure value. The C-N connectivity was maintained throughout in both series of calculations, indicating that the ILB0 and ILB1 phases correspond to separate functions in the potential energy surface of the system and that their interconversion is an activated process that does involve interlayer but not intralayer C-N bond breaking or formation. Comparing the calculated enthalpies we identified 47 GPa as the equilibrium transition pressure between the ILB0 and ILB1 polymorphs at T = 0 K (Fig. 6). Following the compression studies we relaxed the ILB1 structure to ambient pressure, and found that the pillared connectivity was retained. Similar calculations repeated for the hypothetical Cl-free compound showed a similar structural transformation occurring at high pressure, but with the isoenthalpic point between ILB0 and ILB1 phases reduced to ~5 GPa. That result indicates that varying the Cl content of these gCNM phases should cause the transition pressures and recoverability of the pillared-layered material to be adjusted experimentally.

Figure 6

Calculated enthalpy as a function of P for the layered-pillared ILB1 phase relative to the layered graphitic (ILB0) polymorph of C6N9H3·HCl and of the hypothetical Cl-free C6N9H3 system.

In order to fully characterize the structural changes associated with compression of n class="Chemical">C6N9H3·HCl from its ambient pressure graphitic phase to the denser ILB1 polymorph at high pressure, we measured a set of dihedral angles in the minimum energy structures relaxed at each pressure. These angles (Fig. 7a) are chosen to describe the distortion of the initially planar graphitic layers: 1) the internal dihedral angle of a protonated C3N3 ring (ABCD); 2) the dihedral angle formed by the skeletal proton (i.e. the bridging NH group) with respect to the neighbouring triazine ring (ABFG); and 3) the internal dihedral angle of a non-protonated C3N3 ring (HIJK). All these values are identically zero for fully planar layers and any deviations thus record buckling of the gCNM framework. Achieving a value of 60° in any of these parameters indicates a change of hybridization at the central C and N atoms from sp2 to sp3. At each pressure reported in Fig. 7 we employ the thermodynamically stable phase from the DFT calculations, i.e. ILB0 below the transition pressure of 47 GPa, and ILB1 above.

Figure 7

(a) Dihedral angles with reference to the initially planar g-C6N9H3·HCl structure as a function of pressure.Here the dihedral angle changes are shown with reference to the fully relaxed and thermodynamically stable ILB0 or ILB1 states at each pressure. The structural phase transformation occurs at 47 GPa. (b) Displacement of the Cl− ions along the c axis relative to the initial graphitic planes within the calculated ILB0 structure.

At ambient pressure the ILB0 structure is nearly planar with all din class="Chemical">hedral angles <2°, as expected for a graphitic phase with delocalized sp2 bonds. However, this condition is lost as soon as the structure is compressed even slightly. By 10 GPa the dihedral angles HIJK and ABFG differ noticeably from zero but those within the C3N3 rings (ABCD) remain constant. This reveals a buckling of the layers around the bridging -NH- groups (ABFG) beginning at low pressure that is accompanied by a displacement of the Cl− ions out of the graphitic planes (Fig. 7b). Between 40–50 GPa the ABCD dihedral angle jumps from 2° to 65°, indicating the sp2 to sp3 hybridization change that has taken place around the C and N atoms associated with the triazine rings. At this point the reduction in the c axis is such that the interlayer C-N bonding has occurred to yield the ILB1 structure. The structural changes can also be correlated with the behavior of the Cl− ions. While these remain within the layers, they are clamped tightly by strong H-bonding interactions with the NH groups around the C12N12 ring sites. On application of pressure, these voids are diminished in size and they become too small to host the Cl− ions that are progressively squeezed out of the plane. This effect becomes more marked with increasing densification until it reaches nearly 0.6 Å by 70 GPa. These Cl− displacements also lead to structural disorder in the real system, as the Cl− ions can be located either above or below the planes and this effect can contribute to the peak broadening observed in the experimental X-ray diffraction patterns. Once removed from the graphitic layers and into the interlayer space the Cl− ions also become free to be displaced along the (a,b) directions.

Comparison with experiment: ILB0 vs ILB1 formation at high pressure

Following our theoretical prediction of formation of a novel ILB1 structured phase at high pressure, we compared our calculated X-ray diffraction patterns for this material with experimental data (Fig. 8). Tn class="Chemical">he diffraction pattern recorded experimentally at 70 GPa showed broad features extending throughout the 5–18° 2θ range, that could not be refined to provide structural details. However comparison with the calculated results for the predicted ILB1 and ILB0 structures at the same pressure is revealing. The predicted X-ray pattern for the optimized ILB1 phase at 70 GPa matches the observed diffraction profile, with the main diffraction intensity concentrated near 13° 2θ, along with other features occurring near 10° and 6.5° 2θ. The pattern is significantly different from that expected for the ILB0 structure at the same pressure, leading us to conclude that the high pressure structural transformation observed experimentally corresponds to a transition between ILB0 and ILB1 phases.

Figure 8

Comparison of the experimental XRD pattern obtained at 70 GPa (top) vs simulated patterns generated from DFT calculations (bottom), derived from geometry optimised ILB1 and ILB0 structures of C6N9H3·HCl, with or without interlayer bonding occurring at the same pressure.

Low pressure V(P) results

The computational results enabled us not only to identify tn class="Chemical">he structural features of the new high-pressure ILB1 polymorph of the carbon nitride material under study, but also to further characterize the evolution of the graphitic ILB0 phase at intermediate pressures, highlighting the occurrence of buckling distortions and displacements of the Cl− ions. We carried out a detailed analysis of our angle dispersive X-ray diffraction data below 36 GPa to study the compression behavior of the ILB0 structure and validate the computational results in relation to the experimental measurements. Refinements were performed within P63/m using PowderCell software by first assigning the position of the main diffraction peak to the (002) reflection of the ABAB layer-ordered structure, using the DFT results to guide the structural analysis in the highest pressure ranges where substantial peak broadening was present (Fig. 2). The variations in refined a,b,c cell parameters led to the V(P) plot shown in Figure 9.

Figure 9

V(P) plot obtained from our experimental data (circles) for g-C6N9H3.HCl analyzed within the spatially averaged space group P63/m.

The data were fitted up to 21 GPa using a third order Birch-Murnaghan equation of state (K = 32.8 ± 7.3 GPa and K′ = 4.8 ± 0.5) and extrapolated to 36 GPa (line).

Between 20–30 GPa we note a slight deviation of the V(P) relation upwards from tn class="Chemical">he expected equation of state (EoS) curve obtained by analyzing the data using the spatially averaged space group P63/m. This can occur due to layer buckling along with displacements of the Cl− ions away from the graphitic planes, that become significant above 20 GPa as shown by the DFT results (Fig. 7b). We attempted to analyze our data in this range using lower symmetry space group solutions P6 P6 and P1 but these all gave similar Rwp and χ2 values despite the increased structural degrees of freedom.

We fitted the V(P) data up to 21 GPa using a third order Birch-Murnaghan (BM-3) equation of state: with K tn class="Chemical">he zero pressure bulk modulus and K′ its pressure derivative, and the volume strain variable expressed by

This fit yielded K = 32.8 ± 7.3 GPa and K′ = 4.8 ± 0.5, in good agreement with our DFT results for the ILB0 polymorph (K = 29.1 GPa; K′ = 6.7). Tn class="Chemical">he compressibility is comparable with crystalline graphite (K = 34 GPa), although polycrystalline samples have been reported with K values extending down 10 GPa19.

Fitting our DFT V(P) results for the high-pressure ILB1 n class="Chemical">C6N9H3·HCl structure to a BM-3 EoS yields a bulk modulus K0 = 130 ± 5 GPa and K0′ = 3.57 ± 0.2, clearly showing that the ILB1 structure has much lower compressibility than the sp2 bonded graphite-like ILB0 phase (K0 = 29.1 ± 4.1 GPa). The new high pressure ILB1 phase, despite partly retaining the void intralayer structure of the low pressure ILB0 polymorph, has therefore substantially higher structural strength, a feature of potential interest for practical applications.

Discussion

We conducted a high pressure study of graphitic n class="Chemical">C6N9H3·HCl using complementary experimental and theoretical approaches. The graphitic layered compound C6N9H3·HCl was subjected to pressures of up to 70 GPa in a diamond anvil cell, and its structural behaviour was examined using synchrotron X-ray diffraction. The structural model proposed by Zhang et al. with a layered hexagonal (P63/m) structure was employed for the experimental data analysis, for which it provides a useful description up to a pressure of ~40 GPa. Close examination of the V(P) relation indicates an anomaly near 20 GPa that can be explained by a deviation of the structure from P63/m symmetry that can be interpreted from the computational results. On increasing pressure from ambient conditions, we observe three distinct structural rearrangements: at ambient pressure the carbon nitride layers are atomically flat and host the Cl− ions within the layers; between 10–20 GPa the layers buckle and the Cl− ions displace vertically but not horizontally with respect to the layers; above 20 GPa the Cl− ions can displace significantly both horizontally and vertically, making the P63/m model less and less appropriate for a detailed structural description.

The calculations identify tn class="Chemical">he formation of a new pillared structure (referred to as ILB1), that is thermodynamically stable compared to the graphitic phase above 47 GPa. Simulated XRD patterns using graphitic and pillared phases from the computational work were compared with the experimental XRD data, and these indicate that the predicted pillared structure provides the most plausible match to the observed data above 40 GPa. This pillared structure is formed by buckled sp2-hybridized platelets joined by sp3-hybridised "tubes" running along the (001) direction (Fig. 5). The close match between the experimental data and that of the ILB1 structure suggests that the initially graphite-like system undergoes a change in hybridization and experiences interlayer C-N bonding as predicted by the calculations.

In conclusion, our experimental and DFT results indicate that compression of the layered compound g-n class="Chemical">C6N9H3·HCl at ambient temperature results in a polymorphic transformation leading to a new interlayer bonded phase. This has mixed sp2–sp3 bonding between C and N atoms in adjacent layers and it gives rise to a new family of pillared-layered compounds that can combine the essential characteristics of the 2D material with buckled CxNy(Hz) planes and the mechanical properties of a 3D-bonded compound.

Methods

The starting compound g-n class="Chemical">C6N9H3·HCl was formed from a 1:2 mixture of melamine (C3N6H6) and cyanuric chloride (C3N3Cl3) reacted under high-P,T conditions (1.5 GPa, 550°C) in a piston-cylinder device. The synthesis details and structural characterization of this material were described previously818.

Synchrotron X-ray diffraction studies were carried out at up to 70 GPa using angle dispersive techniques at ESRF ID27 (Grenoble, France: λo = 0.3738 Å: Figure 2). Additional data were obtained at Diamond Light Source I15, UK (λo = 0.4441 Å). Samples were loaded into mechanically (four post Mao-type symmetrical cells) or gas-driven He membrane diamond anvil cells (DACs). Diamonds with 200–300 μm culet sizes were used to collect data up to 70 GPa. Re gaskets pre-indented to ~30 μm with gasket holes ~80 μm in diameter were drilled by electro-erosion. Samples were packed into tn class="Chemical">he cell chamber inside a glove box without a pressure-transmitting medium (PTM). Ruby fluorescence spectra indicated that the pressurisation environment remained reasonably "hydrostatic" throughout most of the runs: pressure gradients of only ~1–2 GPa across the cell were detected and the ruby signals showed little broadening. The samples were loaded without a PTM to avoid interference with the weak and often broadened diffraction signals from the low-Z carbon nitride material, and to avoid potential effects associated with intercalation of fluid species such as He, N2 or H2 into the layered structure. X-ray diffraction data were recorded using MAR image plate or CCD detectors and analysed using Fit2D20, FullProf21 and GSAS2223 software.

Electronic structure calculations were carried out under periodic boundary conditions within the local density approximation (n class="Chemical">LDA) of density functional theory (DFT) using the planewave code CASTEP24. We used ultra-soft pseudopotentials and chose the cut-off energy (450 eV) and k-point grid size (3 × 3 × 4) to obtain best converged results25. Calculations were performed under hydrostatic pressure conditions: full geometry optimizations were carried out at each pressure within space group P1 using a 40-atom unit crystallographic cell. Prior to each geometry optimization, the atoms were subjected to small random displacements to produce symmetry-free initial structures.

Author Contributions

A.S. carried out experiments and analysis assisted by R.Q.C. Calculations were carried out by F.C. and M.D., P.F.M., A.S. and F.C. conceived and wrote the paper.