Strange India All Strange Things About India and world

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.

Source link

Leave a Reply

Your email address will not be published.