Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Materials synthesis at terapascal static pressures


Theoretical modelling predicts very unusual structures and properties of materials at extreme pressure and temperature conditions1,2. Hitherto, their synthesis and investigation above 200 gigapascals have been hindered both by the technical complexity of ultrahigh-pressure experiments and by the absence of relevant in situ methods of materials analysis. Here we report on a methodology developed to enable experiments at static compression in the terapascal regime with laser heating. We apply this method to realize pressures of about 600 and 900 gigapascals in a laser-heated double-stage diamond anvil cell3, producing a rhenium–nitrogen alloy and achieving the synthesis of rhenium nitride Re7N3—which, as our theoretical analysis shows, is only stable under extreme compression. Full chemical and structural characterization of the materials, realized using synchrotron single-crystal X-ray diffraction on microcrystals in situ, demonstrates the capabilities of the methodology to extend high-pressure crystallography to the terapascal regime.


The state of matter is strongly affected by variations in chemical composition and external parameters such as pressure and temperature, enabling tuning of material properties. This gives rise to various phenomena relevant for a broad range of scientific disciplines and technological applications, from fundamental understanding of the Universe to targeted design of advanced materials. Compression is known to facilitate metal-to-insulator transitions4, superconductivity5 and new ‘super’ states of matter6. Recent developments in the diamond anvil cell technique, and, particularly, the invention of double-stage and toroidal diamond anvil cells (dsDACs and tDACs)3,7,8, have enabled breakthroughs in the synthesis of materials and the study of structure–property relationships at high and ultrahigh pressures. Very recent examples are the discovery of a new nitrogen allotrope9, bp-N, which resolved a puzzle in our understanding of the high-pressure behaviour of pnictogen family elements, and the synthesis of a plethora of novel transition metal nitrides and polynitrides10,11,12,13,14,15, including metal–inorganic frameworks11,15, which are a new class of compounds featuring open porous structures at megabar compression. Solving and refining the crystal structures of solids synthesized directly from elements in laser-heated conventional DACs10,11,12,13,14,15 at pressures as high as up to about two megabars12,16 became possible owing to the synergy of our expertise both in generating pressures of several megabars3,17,18 (for details see Supplementary Information section ‘Brief overview of the double-stage DAC (dsDAC) technique’) and in single-crystal X-ray diffraction (XRD) at ultrahigh pressures, which were pioneered a few years ago19,20. As the high-pressure high-temperature synthesis has become a well established technique for materials discovery, extending investigations to the TPa regime has long been desired.

Here we report a methodology for high-pressure high-temperature synthesis experiments that extends the limits of high-pressure crystallography to the terapascal range. To achieve the desired pressures, we combined toroidal7,8 and double-stage3,17,18 anvil designs. A rhenium–nitrogen alloy and rhenium nitride Re7N3 were synthesized in three different experiments in the Re–N system (Supplementary Table 1) in a laser-heated dsDAC. Their full structural and chemical characterization was performed in situ using single-crystal XRD.

The dsDACs were prepared following the procedure outlined below. Conventional Boehler–Almax-type single-bevelled diamond anvils with 40-μm culets were milled by focused ion beam (FIB) in order to produce a toroidal profile on the surface of the culet and to shape a miniature culet of about 10 μm in diameter in its centre (Extended Data Fig. 1). As a gasket we used a strip of a 200-μm-thick Re foil, which was pre-indented in a few steps. The final indentation of 10 μm in diameter (made using anvils with the toroidal profile) had a thickness of about 4 μm (the indentation procedure is described in detail in the legend to Extended Data Fig. 1). A hole of approximately 6 μm in diameter was made in the centre of the indentation using FIB or by tightly focused pulsed near-infrared laser to form a pressure chamber. A schematic of the dsDAC assembly, mounted into a BX-90 DAC21 equipped with toroidal diamond anvils, is shown in Extended Data Fig. 1. To realize a dsDAC design, two transparent nanocrystalline diamond17 hemispheres, FIB-milled from a single ball with a diameter of 12 to 14 μm, were placed over the tip of the 10-μm culet (Extended Data Figs. 1, 2). The hemispheres were small enough to stick on the toroidal anvils, but in one case (dsDAC #2, Supplementary Table 1) paraffin wax was used to affix them. A few grains of a rhenium powder (99.995% purity, Merck) were placed into the pressure chamber, which was then filled with nitrogen (N2) at about 1.4 kbar using the high-pressure gas-loading set-up22 at Bayerisches Geoinstitut (BGI, Bayreuth, Germany), closed, and pressurized.

After closing the cells in the pressure chambers, pressures were about 50 to 80 GPa (Extended Data Fig. 3); pressures on the primary anvils were below 10 GPa, as measured according to refs. 23,24. Our experience suggests that the cell should be pressurized quickly to approximately 40 GPa on the primary anvils to avoid loss of nitrogen. The presence of nitrogen can be monitored on N2 vibrons in the Raman spectra (Extended Data Fig. 3). However, N2 vibrons were not detectable above approximately 150 GPa (Extended Data Fig. 3) in the pressure chamber, because at such compression nitrogen becomes non-transparent and we can no longer detect the Raman signal. In dsDAC #2 we were able to observe the evolution of the Raman signal from the secondary anvil in parallel with that from the primary anvil upon pressurization (Extended Data Fig. 4). Huge stress on the secondary anvil is manifested in the large asymmetry of its corresponding Raman line, the high-frequency edge of which is difficult to determine reliably (Extended Data Fig. 4). Thus, it cannot be used for characterization of pressure in the sample chamber. (We also note that, as a rule, Raman spectra of nanocrystalline diamond are somewhat weak and broad).

In all dsDAC experiments described here, we followed the same protocol. After pressurization of the cells to about 120–140 GPa on the first-stage anvils24, the samples were laser-heated. The dsDACs #2 and #3 were heated by a pulsed laser (1-μs pulse duration, 25-kHz repetition rate, approximately 25 W at each side) at BGI using the set-up specially designed for ultrahigh pressures: the near-infrared (1,070 nm) laser beam is of less than 5 μm full-width at half-maximum (FWHM) in diameter and has an optical magnification of about 300×25,26. The entire pressure chamber of dsDAC #2 was heated at 2,900(200) K for about 3 min, and dsDAC #3 at 3,450(200) K for about 5 min. After laser-heating, the pressures on the primary anvils of dsDAC #2 and dsDAC #3 were about 100 GPa and 120 GPa, respectively.

The dsDAC #1 was heated at 13-IDD at GSECARS (Advanced Photon Source, USA) from both sides using a tightly focused near-infrared laser beam (FWHM of about 8 μm in diameter) in pulsed mode (1-μs pulse duration, 50-kHz repetition rate, approximately 20 W each side) for 5 s at a temperature of 2,200(200) K. Powder diffraction data acquired before laser-heating (Extended Data Fig. 5; at 13-IDD the X-ray beam had a FWHM of approximately 3 × 3 μm2) gave the following lattice parameters for Re: for the gasket, a = 2.5606(5) Å, c = 4.0588(12) Å, V = 23.047(7) Å3, and for the Re sample, a = 2.2214(3) Å, c = 3.5609(8) Å, V = 15.21(1) Å3. These parameters correspond to pressures of 149(3) GPa on the gasket and 930(5) GPa on the sample; the conservative values are given according to the equation of state from ref. 27 (Supplementary Table 1; the uncertainty in pressure corresponds to the statistical error in volume). X-ray powder diffraction patterns collected after laser-heating show that the positions of the diffraction lines of the Re gasket did not change within the accuracy of the measurements, and those from the Re sample changed very slightly (a = 2.2297(2) Å, c = 3.5735(5) Å, V = 15.38(1) Å3) corresponding to a pressure of 895(5) GPa (ref. 27).

After laser-heating for each dsDAC at 13-IDD at GSECARS, numerous diffraction spots were observed (Extended Data Fig. 5), indicating phase transformation(s) and/or chemical reaction(s) in the samples. However, interpreting the powder diffraction data turned out to be impossible, as the patterns were dominated by the diffraction lines from the gasket and untransformed Re, owing to the relatively large X-ray beam and a small sample size. Single-crystal diffraction data were of poor quality that precluded their analysis.

The dsDACs with temperature-quenched material were transported to ID11 at the European Synchrotron Radiation Facility (ESRF, Grenoble, France) and investigated using both powder and single-crystal XRD (see Methods). Despite the nominally small size of the X-ray beam, the reflections from the gasket were present even in the patterns collected from the centre of the sample chamber. Two-dimensional (2D) diffraction maps of still XRD images revealed powder diffraction of the Re gasket and untransformed material that enabled the analysis of the pressure distribution both within and around the sample (Extended Data Fig. 2). In dsDAC #1, for example, pressure at the sample/gasket boundary did not exceed approximately 160 GPa, and pressure at all points within the sample chamber was almost the same, of about 900 GPa (Extended Data Fig. 2). Our observations regarding the pressure distribution (Extended Data Fig. 3) in the sample chamber are consistent with those previously reported for toroidal-type anvils7,8 and give the pressure magnification factor (the ratio of the pressures on the primary and secondary anvils) of about 6, in accordance with previous publications on ds-DACs17,28.

Apart from powder diffraction rings, the diffraction patterns collected at ID11 from certain locations in the sample area show numerous spots (Fig. 1). At these positions we collected single-crystal datasets upon rotation of the DAC around the ω axis from −38° to 38° with an angular step of 0.5° (Methods). For dsDAC #1, particularly, the analysis of single-crystal XRD data revealed the presence of domains of two phases (Supplementary Table 2). The first phase is hexagonal (space group P63/mmc) with lattice parameters a = 2.2269(4) Å, c = 3.5702(15) Å and V = 15.33(1) Å3, as determined using 64 reflections. This was interpreted as Re (Figs. 1, 2) being under a pressure of 905(5) GPa (ref. 27). Within uncertainty, the c/a ratio (1.603(5)) coincides with that reported for pure Re at lower pressures3,27. The structure solution and refinement showed indeed that rhenium recrystallizes upon pulsed laser-heating (Fig. 2 and Supplementary Table 2), but is not contaminated by carbon or nitrogen (at least in the quantities that could be detectable from our XRD data).

Fig. 1: Results of XRD measurements on the sample of Re and N2 pulsed laser-heated in dsDAC #1.
figure 1

a, X-ray 2D map showing the distribution of different phases (recrystallized Re and Re7N3) in the pressure chamber of dsDAC #1. Each pixel on the map corresponds to a 2D XRD pattern collected at the Frelon 4M detector at the ID11 beamline at ESRF (beam size FWHM approximately 0.45 × 0.45 μm2, λ = 0.3099 Å). The map covers the whole pressure chamber (21.5 × 21.5 μm2, steps of 0.5 μm in both directions, 10-s acquisition time per frame). The total collection time was about 8 h. The colour intensity is proportional to the intensity of the following reflections: the (100) reflection of the Re gasket for the dark blue region; the (101) reflection of Re for the light blue region (inside the sample chamber); the inset colour bar corresponds to the sum of intensities of (202) and (420) reflections of Re7N3. b, Example of an as-collected diffraction image with diffraction lines and spots of Re (a = 2.2269(4) Å, c = 3.5702(15) Å) and Re7N3 (a = 6.2788(2) Å, c = 4.000(2) Å). The characteristic diffraction image shown in b is highlighted with a white rectangle in a. c, d, The reconstructed reciprocal lattice planes of Re (c) and Re7N3 (d). In c, d, the reflections of Re and Re7N3 are marked by yellow and green circles, respectively, and the corresponding hkl are given. Powder diffraction lines are due to the Re gasket and untransformed rhenium. In bd, blue circles and the blue rectangle indicate parasitic reflections from diamond anvils.

Fig. 2: Crystal structures of the phases observed in laser-heated dsDACs.
figure 2

a, Hexagonal rhenium at 905(5) GPa in dsDAC #1 (a = 2.2269(4) Å, c = 3.5702(15) Å and V = 15.33(1) Å3). b, Cubic (B1 NaCl-type) rhenium–nitrogen solid solution ReN0.2 at 730(4) GPa (a = 3.3994(7) Å, V = 39.28(2) Å3). c, Hexagonal Re7N3 (a = 6.2788(2) Å, c = 4.000(2) Å and V = 136.53(11) Å3). In Re7N3, the structural units are NRe6 prisms with the nitrogen atom in the centre. Rhenium atoms are grey and nitrogen atoms are blue.

The second phase found in the pressure chamber of the dsDAC #1 after heating is also hexagonal (space group P63mc) and has lattice parameters a = 6.2788(2) Å, c = 4.000(2) Å and V = 136.53(11) Å3. On the basis of 394 independent reflections, the structure of this phase was solved and refined in isotropic approximation of atomic displacement parameters (Fig. 2 and Supplementary Table 2) to R1 = 5.7%. The chemical composition of the phase was refined as Re7N3. Considering the possibility of the reaction between rhenium and carbon from the anvils, we checked if the phase could be interpreted as carbide (Re7C3). In this case, however, the isotropic thermal parameter of carbon becomes negative, supporting the assignment of the atomic positions to nitrogen.

The structure units of Re7N3 are distorted NRe6 trigonal prisms (Fig. 2). Three prisms are connected through shared edges forming triads, which are stacked along the 63 axis. Each triad is rotated by 60° with regard to upper and lower neighbours in the columns (Fig. 2). The columns are connected to each other by the common vertices of the prisms. Crystal structures built of combined triads of prisms are well known among carbides, borides, phosphides and nitrides29. Moreover, there are a number of binary compounds with the A7X3 stoichiometry (A and X are different chemical elements), and especially hexagonal ones with Th7Fe3-type structure (more than 70 entries in the ICSD database)30, the same as that of the Re7N3 compound. We noticed that in Re7N3, the shortest and average distances between the Re–Re nearest neighbours (approximately 2.28 Å and 2.37 Å, respectively) are just slightly longer than the Re–Re distances in metallic rhenium (about 2.23 Å), which is present in the pressure chamber along with the nitride. A comparison of the shortest and average distances between the closest AA neighbours in the Th7Fe3-type structured compounds with the metal–metal distances in corresponding pure metals at the same pressures (Extended Data Fig. 6) indeed shows a clear similarity. (In some cases—for example, in experimentally studied Fe7C3 at 158 GPa (ref. 31), or theoretically predicted Fe7N3 at 150 GPa (ref. 32)—the AA distances are even slightly shorter in compounds than in pure metals). Notably, the average Re–N distance in NRe6 prisms in Re7N3 (Re–N is 1.84 Å) follows the same trend as for other Th7Fe3-type structured compounds when AX is compared with AA (Extended Data Fig. 6). According to our experimental data, the Re–N distances in trigonal prisms in Re7N3 vary from approximately 1.79 Å to 1.94 Å, as expected for pressures of several megabars (the shortest previously reported rhenium–nitrogen distance is approximately 1.96 Å in ReN8·xN2 at 134 GPa)11. We note that in the TPa pressure range, the Re–Re interatomic distances become comparable with those of transition metals of the fourth period (Cr, Mn, Fe, Ni), which are known to form Th7Fe3-type structured (or similar) compounds at ambient (or relatively low) pressure30. It may be an indication that a huge reduction of the Re size promotes formation of Re7N3 at several hundreds of GPa, but the existence of Ru7B3 at ambient pressure30 (in ruthenium the metal–metal distance is approximately 2.68 Å versus approximately 2.75 Å in Re) suggests that the size factor may be important, but not necessarily crucial.

The synthesis of Re7N3 was reproduced in dsDAC #2. Diffraction data collected at ID11 at ESRF shows numerous diffraction spots, and the analysis of the integrated powder diffraction pattern confirmed the presence of the hexagonal phase with the lattice parameters very close to those obtained for Re7N3 in dsDAC #1 (Supplementary Tables 1, 3 and Extended Data Fig. 7). Unfortunately, the quality of the diffraction was insufficient for the single-crystal data analysis; the deterioration of the quality of diffraction data may be due to a pressure drop from around 140 GPa to 100 GPa on primary anvils upon laser-heating. Still, for dsDAC #2 we were able to release pressure to ambient without total destruction of the pressure chamber and found there a particle of almost 2 μm in diameter, which consisted of Re and N in the atomic ratio of about 2:1 (Extended Data Fig. 8). This finding provides additional evidence of the synthesis of rhenium nitride in dsDAC #2.

To elucidate the effect of the extreme compression on the stability of the Re7N3 compound and to characterize its physical properties, we carried out electronic structure calculations in the framework of density functional theory and studied its electronic, thermodynamic and vibrational properties (see Methods and Supplementary Information section ‘Computational details’). The optimized lattice parameters and coordinates of atoms of Re7N3 were found to be in excellent agreement with experiment (Supplementary Table 4). A difference in pressure calculated at experimental volumes for Re7N3 may indicate that the calculated equation of state of Re and/or Re7N3 at ultrahigh compressions is becoming less accurate, which is often the case in generalized gradient approximation calculations. Examination of the electronic band structure (Supplementary Information section ‘Electronic properties of Re7N3’ and Supplementary Fig. 1), electronic density of states (Supplementary Figs. 2, 3), electron localization function (Supplementary Fig. 4), and charge density maps (Supplementary Fig. 5) show that Re7N3 is a metal that has a combination of metallic and ionic bonding with some covalent component.

The dsDAC #3 was laser-heated to a maximum temperature of 3,450(200) K and the lattice parameters of Re measured after heating were found to be a = 2.2803(3) Å, c = 3.622(1) Å and V = 16.31(2) Å3. According to the equation of state27 of Re, the sample was under pressure of 730(4) GPa (Supplementary Table 1 and Supplementary Fig. 6). The analysis of single-crystal XRD data revealed the presence of a cubic phase (space group \({Fm}\bar{3}m\)) with a lattice parameter of approximately 3.40 Å to approximately 3.41 Å depending on the spot from which the XRD pattern was taken. Structural solution suggests that the phase has an NaCl (B1)-type structure (Fig. 2 and Supplementary Fig. 7) with one position occupied by Re and the other by a light element. Attempts to refine the crystal structure assuming that the position of the light element is fully occupied by N or C led to an unreasonably high thermal parameter (approximately 0.1 Å2). For the highly symmetric NaCl-type structure containing heavy Re atoms, simultaneous refinement of the occupancy and the thermal parameter of a lighter element is not reasonable, so we constrained the thermal parameters of all atoms to be equal. In this approximation, the composition of the cubic phase was ReN0.20 (Supplementary Table 2). Of course, on the basis of XRD data alone we could not exclude that the light element might be carbon, but theoretical calculations (see Supplementary Information section ‘Re-based solution phase’) suggest that nitrogen is more plausible. A partial occupation of octahedral voids of the underlying face-centred cubic (fcc) packing of Re atoms by nitrogen predicts negative formation enthalpies of metastable alloys (Supplementary Figs. 8, 9 and Supplementary Table 5), whereas filling them with carbon leads to positive formation enthalpies (Supplementary Fig. 8 and Supplementary Table 6).

Theoretical simulations enabled an insight into the possibility of synthesizing Re7N3 at pressures lower than those achieved in the current study. At 100 GPa the formation enthalpy of metastable Re7N3 is well above the convex hull (Fig. 3, Supplementary Information section ‘Thermodynamic stability of Re7N3’ and Extended Data Fig. 9). Even taking into account the anomalously large (approximately 0.2 eV per atom) metastability range of nitrides33, this compound cannot be considered as synthesizable at 100 GPa. By contrast, at 730 GPa the calculated formation enthalpy of Re7N3, although still above the convex hull, becomes well within the metastability range of nitrides (Fig 3, Supplementary Information section ‘Lattice dynamics of Re7N3’ and Extended Data Fig. 9), and at approximately 900 GPa—the pressure of the realized experimental synthesis—it lies on the convex hull (Fig. 3).

Fig. 3: Formation enthalpy of Re7N3.
figure 3

ac, Data are shown with respect to theoretically predicted34 (black squares) and experimentally known (red squares, Re3N and ReN213 (P21/c), ReN2 (P4/mbm), ReN1011 (Immm)) competing high-pressure phases in the ReNx system, calculated at pressures of 100 GPa (a), 730 GPa (b) and 900 GPa (c). hcp, hexagonal close-packed; CG-type N, cubic gauche nitrogen.

Pressures of more than several megabars have long been thought to have a profound effect on the chemistry and physics of materials1,2 and to lead to formation of phases with exotic crystal structures. In this work we have demonstrated that at pressures as high as those exceeding 600 GPa new compounds can be synthesized in laser-heated dsDACs and their structures can be solved in situ. By extending the experimental field of high-pressure synthesis and structural studies to the terapascal range, our work paves the way towards the discovery of new materials and observations of novel physical phenomena.


Diffraction data were acquired at ID11 beamline at ESRF. Experiments with dsDAC #1 were performed using a Frelon 4M detector, wavelength 0.3099 Å, beam size 0.45 × 0.45 μm2 at FWHM; data for dsDAC #2 and dsDAC #3 were collected with Eiger2 CdTe 4M detector, wavelength 0.2882 Å, beam size 0.5 × 0.5 μm2 at FWHM. 2D mappings of still XRD images (without ω oscillations) were performed with an exposure time up to 10 s; single-crystal datasets were collected via DACs rotation around the ω axis from −38° to 38° with an angular step of 0.5° and an acquisition time of 10 s per step.

For the powder diffraction studies, calibration of instrument model and integration of diffraction patterns were made in the DIOPTAS35 software using CeO2 powder standard (NIST SRM 674b). Integrated patterns from powder XRD experiments were processed using the Le Bail technique implemented in JANA200636 software. XRD imaging of the sample chamber was reconstructed using XDI37 programme and map of still images converted from ‘edf’ to ‘tif’ format. For the single-crystal XRD, integration of the reflection intensities and absorption corrections were performed in CrysAlisPro software38. A single crystal of orthoenstatite (Mg1.93,Fe0.06)(Si1.93,Al0.06)O6 (space group Pbca, a = 8.8117(2) Å, b = 5.18320(10) Å, c = 18.2391(3) Å) was used as calibration standard for refinement of the instrument model of the diffractometer. Diffraction images were converted from ‘edf’ to the native CrysAlisPro format ‘ESPERANTO’ with Freac software38. Detailed information of integration parameters as well as of the data-reduction output files and indicators of the XRD data quality are given in ref. 19. The crystal structures were solved using SHELXT or the superflip method in JANA2006 and Olex236,39,40. Crystal structures were refined by least-squares minimization of adjustable parameters. We performed isotropic refinement of atomic displacement parameters due to limited dataset collected in DAC. Reflections coming from parasite diffraction produced by diamonds and crystallized pressure media were eliminated during the refinement procedure. The software Diamond41 was used for visualization of molecular graphics.

The electronic structure, total energy and forces calculations of the studied rhenium nitrides were carried out in the framework of density functional theory (see Supplementary Information section ‘Computational details’). We used supercells of different sizes with an underlying fcc crystal structure and various amounts of either N or C to simulate the Re–N and Re–C cubic phases with NaCl (B1)-type structure (see Supplementary Information section ‘Re-based solution phase’). To investigate the influence of pressure on the thermodynamic stability of Re7N3, its enthalpy of formation, as well as the enthalpies of formation for various phases of rhenium nitride, known experimentally10,11,13 and predicted theoretically34, were calculated and a thermodynamic convex hull was constructed based on the calculation results (Supplementary Information section ‘Thermodynamic stability of Re7N3’).

Phonon dispersion relations for Re7N3 were calculated in the harmonic approximation at volume 200 Å3 (a = 7.122 Å, c = 4.553 Å) of the unit cell, corresponding to P = 102 GPa, as well as at experimental volume 136.52 Å3 (a = 6.277 Å, c = 4.001 Å) of the unit cell (Supplementary Table 4), which corresponded to calculated pressure 732 GPa (see Extended Data Fig. 9 and Supplementary Information section ‘Computational details’). Because Re7N3 is predicted to be dynamically unstable at the synthesis pressure owing to the presence of imaginary frequencies in this approximation (Extended Data Fig. 9 and Supplementary Information section ‘Lattice dynamics of Re7N3’), we investigated the anharmonic effects of lattice vibrations at finite temperature using the temperature-dependent effective potential (TDEP) method42 with effective second-order and third-order interatomic force constants calculated from first principles43. The calculations are based on modelling the potential energy surface in the vicinity of equilibrium with a Hamiltonian of the form:

$$H={U}_{0}+\sum _{i}\frac{{{\bf{p}}}_{i}^{2}}{2{m}_{i}}+\frac{1}{2!}\sum _{ij\alpha \beta }{\varPhi }_{ij}^{\alpha \beta }{u}_{i}^{\alpha }{u}_{j}^{\beta }+\frac{1}{3!}\sum _{ijk\alpha \beta \gamma }{\varPhi }_{ijk}^{\alpha \beta \gamma }{u}_{i}^{\alpha }{u}_{j}^{\beta }{u}_{k}^{\gamma }+\ldots ,$$

where p and m are the momentum and the mass of ion i, respectively, Ф are interaction parameters (the effective force constants) ofincreasing order, u denotes the displacement of ions (ij or k) from their equilibrium positions, and αβγ are Cartesian components.

We calculated the spectral function S(q, E) at 300 K, which describes the spectrum of lattice excitations with energy E = ħΩ (Ω is the applied frequency) for mode s with harmonic frequency ωqs at wavevector q (refs. 44,45). S(q, E) provides insight into the phonon frequencies as well as strength of three-phonon processes via the broadening in Extended Data Fig. 9. The S(q, E) of Re7N3 is typical of a weakly anharmonic solid with Lorentzian broadening of single peaks. Additionally, the lines are reasonably crisp, without substantial broadening, indicating that the anharmonic interaction strength is well within the range of validity for the perturbation theory. Importantly, Re7N3 is seen to be dynamically stable (there are no imaginary frequencies) at the synthesis pressure (see Supplementary Information section ‘Lattice dynamics of Re7N3’).

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this paper.

Data availability

Data supporting this work are available at Zenodo, Structural data deposit at Cambridge Crystallographic Data Centre (CCDC), CSD-2143754 (

Code availability

The temperature-dependent effective potential method is implemented as a package that deals with finite-temperature lattice dynamics in solids. The package is released under the MIT license, available on GitHub, see Source files for TDEP 1.1 are located at


  1. McMillan, P. F. Chemistry at high pressure. Chem. Soc. Rev. 35, 855–857 (2006).

    CAS  Article  Google Scholar 

  2. Needs, R. J. & Pickard, C. J. Role of structure prediction in materials discovery and design. APL Mater. 4, 053210 (2016).

    ADS  Article  Google Scholar 

  3. Dubrovinsky, L., Dubrovinskaia, N., Prakapenka, V. B. & Abakumov, A. M. Implementation of micro-ball nanodiamond anvils for high-pressure studies above 6 Mbar. Nat. Commun. 3, 1163 (2012).

    ADS  Article  Google Scholar 

  4. Austin, I. G. & Mott, N. F. Metallic and nonmetallic behavior in transition metal oxides: electron correlation effects in narrow d bands and polarons are discussed. Science 168, 71–77 (1970).

    ADS  CAS  Article  Google Scholar 

  5. Shimizu, K. et al. Superconductivity in the non-magnetic state of iron under pressure. Nature 412, 316–318 (2001).

    ADS  CAS  Article  Google Scholar 

  6. Kim, E. & Chan, M. H. W. Probable observation of a supersolid helium phase. Nature 427, 225–227 (2004).

    ADS  CAS  Article  Google Scholar 

  7. Dewaele, A., Loubeyre, P., Occelli, F., Marie, O. & Mezouar, M. Toroidal diamond anvil cell for detailed measurements under extreme static pressures. Nat. Commun. 9, 2913 (2018).

    ADS  Article  Google Scholar 

  8. Jenei, Z. et al. Single crystal toroidal diamond anvils for high pressure experiments beyond 5 megabar. Nat. Commun. 9, 3563 (2018).

    ADS  Article  Google Scholar 

  9. Laniel, D. et al. High-pressure polymeric nitrogen allotrope with the black phosphorus structure. Phys. Rev. Lett. 124, 216001 (2020).

    ADS  CAS  Article  Google Scholar 

  10. Friedrich, A. et al. Novel rhenium nitrides. Phys. Rev. Lett. 105, 085504 (2010).

    ADS  Article  Google Scholar 

  11. Bykov, M. et al. High-pressure synthesis of a nitrogen-rich inclusion compound ReN8·xN2 with conjugated polymeric nitrogen chains. Angew. Chem. Int. Ed. 57, 9048–9053 (2018).

    CAS  Article  Google Scholar 

  12. Bykov, M. et al. Synthesis of FeN4 at 180 GPa and its crystal structure from a submicron-sized grain. Acta Crystallogr. E 74, 1392–1395 (2018).

    CAS  Article  Google Scholar 

  13. Bykov, M. et al. High-pressure synthesis of ultraincompressible hard rhenium nitride pernitride Re2(N2)(N)2 stable at ambient conditions. Nat. Commun. 10, 2994 (2019).

    ADS  Article  Google Scholar 

  14. Bykov, M. et al. Synthesis of arsenopyrite-type rhodium pernitride RhN2 from a single-source azide precursor. Eur. J. Inorg. Chem. 2019, 3667–3671 (2019).

    CAS  Article  Google Scholar 

  15. Bykov, M. et al. High‐pressure synthesis of metal–inorganic frameworks Hf4N20·N2, WN8·N2, and Os5N23·N2 with polymeric nitrogen linkers. Angew. Chem. Int. Ed. 59, 10321–10326 (2020).

    CAS  Article  Google Scholar 

  16. Khandarkhaeva, S. et al. Novel rhenium carbides at 200 GPa. Eur. J. Inorg. Chem. 2020, 2186–2190 (2020).

    CAS  Article  Google Scholar 

  17. Dubrovinskaia, N. et al. Terapascal static pressure generation with ultrahigh yield strength nanodiamond. Sci. Adv. 2, e1600341 (2016).

    ADS  Article  Google Scholar 

  18. Dubrovinsky, L. et al. The most incompressible metal osmium at static pressures above 750 gigapascals. Nature 525, 226–229 (2015).

    ADS  CAS  Article  Google Scholar 

  19. Bykova, E. Single-Crystal X-Ray Diffraction at Extreme Conditions in Mineral Physics and Material Sciences. PhD thesis, Univ. Bayreuth (2015).

  20. Dubrovinskaia, N. & Dubrovinsky, L. Crystallography taken to the extreme. Phys. Scr. 93, 062501 (2018).

    ADS  Article  Google Scholar 

  21. Kantor, I. et al. BX90: a new diamond anvil cell design for X-ray diffraction and optical measurements. Rev. Sci. Instrum. 83, 125102 (2012).

    ADS  CAS  Article  Google Scholar 

  22. Kurnosov, A. et al. A novel gas-loading system for mechanically closing of various types of diamond anvil cells. Rev. Sci. Instrum. 79, 045110 (2008).

    ADS  Article  Google Scholar 

  23. Akahama, Y., Hirao, N., Ohishi, Y. & Singh, A. K. Equation of state of bcc-Mo by static volume compression to 410 GPa. J. Appl. Phys. 116, 223504 (2014).

    ADS  Article  Google Scholar 

  24. Akahama, Y. & Kawamura, H. Pressure calibration of diamond anvil Raman gauge to 410 GPa. J. Phys. Conf. Ser. 215, 012195 (2010).

    Article  Google Scholar 

  25. Aprilis, G. et al. Portable double-sided pulsed laser heating system for time-resolved geoscience and materials science applications. Rev. Sci. Instrum. 88, 084501 (2017).

    ADS  CAS  Article  Google Scholar 

  26. Fedotenko, T. et al. Laser heating setup for diamond anvil cells for in situ synchrotron and in house high and ultra-high pressure studies. Rev. Sci. Instrum. 90, 104501 (2019).

    ADS  Article  Google Scholar 

  27. Anzellini, S., Dewaele, A., Occelli, F., Loubeyre, P. & Mezouar, M. Equation of state of rhenium and application for ultra high pressure calibration. J. Appl. Phys. 115, 043511 (2014).

    ADS  Article  Google Scholar 

  28. Sakai, T. et al. High pressure generation using double-stage diamond anvil technique: problems and equations of state of rhenium. High Press. Res. 38, 107–119 (2018).

    ADS  CAS  Article  Google Scholar 

  29. Li, W.-K., Zhou, G.-D. & Mak, T. Advanced Structural Inorganic Chemistry (Oxford Univ. Press, 2008).

  30. ICSD database;

  31. Prescher, C. et al. High Poisson’s ratio of Earth’s inner core explained by carbon alloying. Nat. Geosci. 8, 220–223 (2015).

    ADS  CAS  Article  Google Scholar 

  32. Sagatov, N., Gavryushkin, P. N., Inerbaev, T. M. & Litasov, K. D. New high-pressure phases of Fe7N3 and Fe 7C3 stable at Earth’s core conditions: evidence for carbon–nitrogen isomorphism in Fe-compounds. RSC Adv. 9, 3577–3581 (2019).

    ADS  CAS  Article  Google Scholar 

  33. Sun, W. et al. The thermodynamic scale of inorganic crystalline metastability. Sci. Adv. 2, e1600225 (2016).

    ADS  Article  Google Scholar 

  34. Zhao, Z. et al. Nitrogen concentration driving the hardness of rhenium nitrides. Sci Rep. 4, 4797 (2015).

    Article  Google Scholar 

  35. Prescher, C. & Prakapenka, V. B. DIOPTAS: a program for reduction of two-dimensional X-ray diffraction data and data exploration. High Press. Res. 35, 223–230 (2015).

    ADS  CAS  Article  Google Scholar 

  36. Petříček, V., Dušek, M. & Palatinus, L. Crystallographic computing system JANA2006: general features. Z. Kristallogr. Cryst. Mater. 229, 345–352 (2014).

    Article  Google Scholar 

  37. Hrubiak, R., Smith, J. S. & Shen, G. Multimode scanning X-ray diffraction microscopy for diamond anvil cell experiments. Rev. Sci. Instrum. 90, 025109 (2019).

    ADS  Article  Google Scholar 

  38. CrysAlis Pro v. 171.40.84 (Rigaku, 2020).

  39. Sheldrick, G. M. SHELXT – integrated space-group and crystal-structure determination. Acta Crystallogr. A 71, 3–8 (2015).

    Article  Google Scholar 

  40. Dolomanov, O. V., Bourhis, L. J., Gildea, R. J., Howard, J. A. K. & Puschmann, H. OLEX2: a complete structure solution, refinement and analysis program. J. Appl. Crystallogr. 42, 339–341 (2009).

    CAS  Article  Google Scholar 

  41. Putz, H. & Brandenburg, K. Diamond – Crystal and Molecular Structure Visualization Version 4 (Crystal Impact, 2021).

  42. Hellman, O., Abrikosov, I. A. & Simak, S. I. Lattice dynamics of anharmonic solids from first principles. Phys. Rev. B 84, 180301 (2011).

    ADS  Article  Google Scholar 

  43. Hellman, O. & Abrikosov, I. A. Temperature-dependent effective third-order interatomic force constants from first principles. Phys. Rev. B 88, 144301 (2013).

    ADS  Article  Google Scholar 

  44. Cowley, R. A. Anharmonic crystals. Rep. Prog. Phys. 31, 123–166 (1968).

    ADS  CAS  Article  Google Scholar 

  45. Maradudin, A. A. & Fein, A. E. Scattering of neutrons by an anharmonic crystal. Phys. Rev. 128, 2589–2608 (1962).

    ADS  CAS  Article  Google Scholar 

Download references


We thank S. Petitgirard (ETH Zürich, Zurich, Switzerland) for help in the FIB preparation of samples. We acknowledge A. Kurnosov, D. Wiesner, S. Übelhack and S. Linhardt for technical assistance. D.L. thanks the Alexander von Humboldt Foundation and the Deutsche Forschungsgemeinschaft (DFG, project LA-4916/1-1) for financial support. N.D. and L.D. thank the Federal Ministry of Education and Research, Germany (BMBF, grant no. 05K19WC1) and the DFG (DFG projects DU 954-11/1, DU 393-9/2 and DU 393-13/2) for financial support. N.D. thanks the Swedish Government Strategic Research Area in Materials Science on Functional Materials at Linköping University (Faculty Grant SFO-Mat-LiU no. 200900971). M.B. acknowledges the support of Deutsche Forschungsgemeinschaft (DFG project BY112/2-1). The European Synchrotron Radiation Facility is acknowledged for beamtime provision at the Material Science beamline ID11. Portions of this work were performed at GeoSoilEnviroCARS (The University of Chicago, Sector 13), Advanced Photon Source (APS), Argonne National Laboratory. GeoSoilEnviroCARS is supported by the National Science Foundation – Earth Sciences (EAR – 1634415). This research used resources of the Advanced Photon Source, a US Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under contract no. DE-AC02-06CH11357. Theoretical analysis of chemical bonding was supported by the Russian Science Foundation (project no. 18-12-00492). Support from the Knut and Alice Wallenberg Foundation (Wallenberg Scholar grant no. KAW-2018.0194), the Swedish Government Strategic Research Areas in Materials Science on Functional Materials at Linköping University (Faculty Grant SFO-Mat-LiU no. 2009 00971), the Swedish e-science Research Center (SeRC), the Swedish Research Council (VR) grant no. 2019-05600, the VINN Excellence Center Functional Nanoscale Materials (FunMat-2) grant 2016–05156 and the Swedish Foundation for Strategic Research (SSF) project no. EM16-0004 is gratefully acknowledged. The computations were carried out at supercomputer cluster at NUST ‘MISIS’ and at resources provided by the Swedish National Infrastructure for Computing (SNIC), partially funded by the Swedish Research Council through grant agreement no. 2016-07213.


Open access funding provided by Universität Bayreuth.

Author information




L.D. and N.D. designed the work. L.D., S.K., D.L. prepared high-pressure experiments. L.D., S.K., T.F., M.B., D.L., C.G., E.L.B., P.S., S.C. and V.P. conducted experiments. L.D., S.K., D.L. and S.C. processed experimental data. A.V.P., E.A.S., M.P.B., F. Tasnádi, N.S., F. Trybel and I.A.A. performed theoretical analysis. The manuscript was written by L.D., N.D. and I.A.A. with contributions from all the authors. All the authors commented on successive drafts and have given approval to the final version of the manuscript.

Corresponding authors

Correspondence to Leonid Dubrovinsky or Igor A. Abrikosov.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature thanks the anonymous reviewers for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Toroidal profile produced by FIB milling on a culet of a conventional anvil.

a, Original Boehler–Almax-type single-bevelled diamond anvils with 40-μm culets (before milling). b, Greyscale bitmap used for milling. c, The milling depth profile. d, Crafted toroidal culet. e, Cross-sectional schematic of the dsDAC assembly (not to scale): a gasket (grey) with the pressure chamber (dark grey; the sample is yellow) squeezed between the two toroidal diamond anvils (light blue) equipped with hemispheres of nanocrystalline diamond (white); the pressure chamber is of 6 μm in diameter and 4 μm in height. The gasket was prepared as follows: We used strips of a 200-μm thick Re foil. To make an indentation with a thickness of about 4 μm, we followed a many-step procedure. First, the Re foil was pre-indented to a thickness of ~20 μm using a pair of single-bevelled diamonds with 80-μm culets. Then, a hole about 30 μm in diameter was laser-drilled in the centre of the indentation, and the indentation was pressurized again between the same bevelled diamonds. This led to closing of the hole and a reduction in the thickness of the indentation. This procedure was repeated a few times until a thickness of about 7 μm was achieved. The gasket was mounted into a BX-90 DAC equipped with toroidal diamond anvil and indented by their miniature 10-μm culets to a thickness of ~4 μm. A hole (~6 μm in diameter) in the centre of the ~10-μm indentation was made using FIB or by tightly focused pulsed near-infrared laser to form a pressure chamber.

Extended Data Fig. 2 A schematic illustration of the experimental set-up.

a, Three-dimensional presentation of a diamond anvil featuring a toroidal profile milled by FIB on the surface of a conventional Boehler–Almax-type single-bevelled diamond anvil with a 40-μm culet. A hemisphere of transparent nanocrystalline diamond (NCD)17, which was FIB-milled from a single ball with a diameter of 12 to 14 μm, was placed over the tip to realize a double-stage DAC design. Two anvils of this kind were forced together as shown in the scheme in Extended Data Fig. 1. b, Pressure profile along the cross-section through the centre of dsDAC #1 after pulsed laser-heating. Diffraction patterns were collected at each point with a step of 0.5 μm at the ID11 beamline at ESRF, and pressure was determined according to the equation of state from ref. 27 using the lattice parameters of Re found from powder diffraction data.

Extended Data Fig. 3 Raman spectra taken upon compression of the dsDAC after nitrogen gas loading featuring nitrogen vibron.

At pressures above ~160 GPa the Raman signal of nitrogen becomes non-detectable. The pressures were determined according to ref. 24. a.u., arbitrary units.

Extended Data Fig. 4 The evolution of Raman spectra of primary and secondary anvils upon compression of dsDAC #2.

The pressure is determined from diamond line Raman shift24. The pressures on primary anvil are shown on the left. The peaks at ~1,600 cm−1 are the pressures on the body of the secondary anvil. The values on the right give estimates of pressures from the broadened Raman line of the secondary anvil (the arrows provide examples of positions found by analysis of the first derivatives of the spectra). Reliable determination of the pressure in the chamber from the Raman shift of the diamond line of the secondary nanocrystalline diamond anvils is not feasible.

Extended Data Fig. 5 Examples of powder diffraction patterns collected from dsDAC #1.

Data collected at 13-IDD (Advanced Photon Source, USA; beam FWHM, 3 × 3 μm2). a, The centre of the pressure chamber of the as-compressed ds-DAC. c, Re LP (low pressure): a = 2.5606(5) Å, c = 4.0588(12) Å and V = 23.047(7) Å3, that is, 149(3) GPa according to the equation of state from ref. 27, or 173(3) GPa according to ref. 3; Re HP (high pressure): a = 2.2214(3) Å, c = 3.5609(8) Å and V = 15.21(1) Å3, that is, 930(5) GPa27 or 1298(10) GPa3. b, After pulsed laser-heating at 2,200(200) K during 5 s. d, Re LP: a = 2.5577(3) Å, c = 4.1095(12) Å and V = 23.282(7) Å3, at 140(3) GPa27 or 162(3) GPa3; Re HP: a = 2.2297(2) Å, c = 3.5735(5) Å and V = 15.38(1) Å3, at 895(5) GPa27 or 1,250(10) GPa3; Re7N3: a = 6.3086(4) Å, c = 4.0048(7) Å and V = 138.04(4) Å3. Structural data for Re7N3 were taken from the results of single-crystal XRD data analysis.

Extended Data Fig. 6 Interatomic distances in Th7Fe3 structured and some other selected A7X3 compounds.

a, b, Comparison of shortest (a) and average between first neighbours (b) AA distances in A7X3 compounds with metal–metal contacts in corresponding pure metals (A) at the same pressures. c, Correlation between average XA distances in XA6 prisms and first neighbours AA distances. The red hexagons correspond to Re7N3 as described in this work; data for orthorhombic Fe7C3 at 158 GPa are from ref. 35, data for predicted Fe7N3 at 150 GPa are from ref. 32, all other data are from ref. 30.

Extended Data Fig. 7 Example of powder diffraction pattern collected from dsDAC #2.

Data were collected at ID11 (ESRF, Grenoble, France, beam FWHM 0.5 × 0.5 μm2) at 646 GPa (see Supplementary Table 1). a, 2D diffraction image shows diffraction spots of Re7N3 and rings of hexagonal close-packed (hcp) Re at different pressures. b, Refinement of powder diffraction pattern with Le Bail fit implemented in JANA2006 software. The values of the lattice parameters are given in Supplementary Table 3.

Extended Data Fig. 8 Characteristics of a sample extracted from the dsDAC #2.

Sample contains Re and N in atomic proportions ~2:1. a, SEM image; b, example of EDX spectra. Images and spectra were collected on a ZEISS SEM, Leo Gemini 1530 with a Schottky field emission gun employing an accelerating voltage of 20 kV.

Extended Data Fig. 9 Phonon dispersion relations for Re7N3 calculated within the harmonic approximation.

a, Theoretical pressure P = 102 GPa. b, Experimental volume 136.52 Å3, which corresponds to theoretical pressure ~730 GPa. c, Vibrational spectral function at P ≈ 730 GPa and T = 300 K, calculated using the temperature-dependent effective potential (TDEP) method. The results show that in the harmonic approximation, Re7N3 is unstable at the synthesis pressure (imaginary frequencies are shown below the zero-frequency line). Including anharmonic effects of lattice vibrations removes the dynamical instability at P = 730 GPa (all the branches are real).

Supplementary information

Supplementary Information

Supplementary Information contains sections describing ‘Computational details’, ‘Electronic properties of Re7N3’, ‘Re-based solution phase’, ‘Thermodynamic stability of Re7N3’, ‘Lattice dynamics of Re7N3’ and ‘Brief overview of the double-stage DAC (dsDAC) technique’; Supplementary Figs. 1–9, Tables 1–7 and references.

Reporting Summary

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Dubrovinsky, L., Khandarkhaeva, S., Fedotenko, T. et al. Materials synthesis at terapascal static pressures. Nature 605, 274–278 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI:


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing