Equilibrium p-T phase diagram of boron: experimental study and thermodynamic analysis.

Solid-state phase transformations and melting of high-purity crystalline boron have been in situ and ex situ studied at pressures to 20 GPa in the 1500-2500 K temperature range where diffusion processes become fast and lead to formation of thermodynamically stable phases. The equilibrium phase diagram of boron has been constructed based on thermodynamic analysis of experimental and literature data. The high-temperature part of the diagram contains p-T domains of thermodynamic stability of rhombohedral β-B106, orthorhombic γ-B28, pseudo-cubic (tetragonal) t'-B52, and liquid boron (L). The positions of two triple points have been experimentally estimated, i.e. β-t'-L at ~ 8.0 GPa and ~ 2490 K; and β-γ-t' at ~ 9.6 GPa and ~ 2230 K. Finally, the proposed phase diagram explains all thermodynamic aspects of boron allotropy and significantly improves our understanding of the fifth element.

The phase diagram of boron and thermodynamic stability of boron allotropes remain of fundamental interest in condensed matter physics and chemistry for a very long time – since the first characterization of pure element's allotropes in 1950s123. Later the existence of multiple boron modifications has been reported (see recent reviews23) but some of them never have been reproduced, e.g. β'- and β″-B4, t-B505, "HP form" of Wentorf6, "HPHT form" of t-B1927, etc. At present time only five allotropes are generally accepted: rhombohedral α-B12 (α-phase)8 and β-B106 (β-phase)9, orthorhombic γ-B28 (γ-phase)10, tetragonal t-B19211 and t-B52 (t-phase in this paper)5121314. t-B52 has been proved to exist only very recently and its crystal structure has not been unambiguously established so far. Two more phases have been predicted using ab initio structural evolution algorithm15, i.e. orthorhombic o-B5216, closely related to t-B52, and metallic boron with α-Ga crystal structure10. The interest to the high-pressure behavior of boron has recently raised due to the discovery of boron superconductivity at high pressure1718; unusual partially ionic character of some B-B bondings in γ-boron10; high-pressure synthesis of novel boron-rich compounds1920212223242526 that are refractory and chemically stable272829, superhard303132, and even could have metallic conductivity33; unusual pressure-induced behavior of boron-containing icosahedral3435 and layered36373839 structures; and prediction of the nonmetal-metal phase transition in boron at a pressure above 89 GPa1040.

The stability of boron allotropes has been intensively investigated during past years using ab initio calculations. It has been predicted that at ambient conditions α-B12 and β-B106 have similar static energies, but disordered β-B106 is more stable at ambient pressure, due to its lower zero-point vibrational energy41. At pressures above 2 GPa, denser α-B12 should be more stable10. The ab initio analysis of stability of boron structures10 showed that at pressures above 20 GPa the α-phase loses its stability and another phase, orthorhombic γ-B28 (confirmed experimentally10), becomes stable. Finally, above 89 GPa transition of semiconductive γ-phase into metallic one should occur10. However, this pressure domain has not been explored experimentally at high temperatures (to overcome a kinetic barrier), and the latter allotrope remains to be discovered. Thus, numerous theoretical predictions require rigorous experimental studies.

The first attempt to analyze the high-pressure phase equilibria in boron was made in 200742 based on ab initio calculations for α-B12 and β-B106 phases and some ambient pressure experimental data. However, the result contradicts the experimental data on boron melting under pressure4344, i.e. dTm/dp is overestimated by a factor of 2 (see Fig. 1c). Moreover, the reported α β equilibrium temperature (Tβ/α) at ambient pressure is lower by ~ 200 K as compared with the maximal temperature of β-to-α recrystallization in the presence of Pt melt45.

Figure 1

(a) Sequence of synchrotron energy-dispersive X-ray diffraction patterns taken in situ in the course of stepwise heating of α-boron at 5.5 GPa (E·d = 78.07 eV Å). Asterisk (*) indicates the position of the escape line of boron nitride (hBN, capsule material). (b) Experimental data on α-to-β transformation in boron. Open and solid circles represent the α- and β-phases observed during our in situ experiments. The up triangle corresponds to the beginning of the α-to-β solid-state transformation4, while the down triangle is the onset of recrystallization of α-B12 from β-B106 in the presence of Pt melt50. The dashed line corresponds to the equilibrium between α-B12 and β-B106 that has been calculated in Ref.42. (c) Experimental data on boron melting. Solid and open squares represent the crystalline and liquid boron observed during our in situ experiments. The crossed squares correspond to the beginning of melting. The open triangles show the literature data: (2350 K, ambient pressure)60, (2370 K, ambient pressure)61 and (2480 K, 7.7 GPa)43. The dashed line represents the melting curve calculated in Ref.42.

The first p-T phase diagram of boron was proposed only in 2009 by Oganov et al.10 and contains 5 allotropes (four experimentally confirmed forms α, β, γ, t-B192, and hypothetical metallic one) and liquid boron. This diagram combined more extended ab initio and experimental data on structural stability and phase relationships, but still some points remained unclear. Though overall correct, this phase diagram contained an uncertainty related to the stability field of the tetragonal boron phase, which at that time was thought to be "HPHT t-B192"7. A second point requiring elaboration is that the equilibrium line between α- and β-boron was an estimate, rather than direct measurement (which would be complicated by kinetics) or calculation (which would be complicated by disorder in β-boron). Recrystallization of α-B12 from β-B106 was experimentally observed at much higher temperatures at both ambient (from melts containing Pt45) and high (during solid-state transformation12) pressures. Finally, t-B52 has been recently obtained at ambient12 and high121314 pressures, and even recovered as a single phase1214. Although this tetragonal allotrope has been interpreted as a metastable one as compared to mysterious "HPHT t-B192", our recent results showed that the latter has a crystal structure related in many aspects to t-B52 phase, rather than to t-B19246.

Very recently pseudo-cubic t'-B52 of the t-B52 structural family has been discovered46. It was recovered after experiments at 20 GPa and 2500 K, the highest temperature reported so far for formation of a tetragonal phase. Contrary to common low-density t-B52 phase(s) and related compounds, pseudo-cubic allotrope is quite dense, very close to γ-B28. This phase seems to be a good candidate for a HPHT allotrope, instead of strongly distorted "t-B192 structure". At such high temperature the diffusion processes are quite intense even at such high pressure as 20 GPa. The observation of reversible transformation in boron at HPHT conditions would be a strong support to equilibrium phase diagram.

In the present work we have studied the high-temperature part of boron phase diagram using in situ and recovery high-pressure experiments (to 20 GPa and 2500 K), as well as thermodynamic analysis. The established equilibria between boron allotropes and liquid phase are self-consistent from the point of view of classical thermodynamics and adequately describe the experimental data obtained by independent research groups.

Results

Phase transformation of α-B12 into β-B106

Our in situ studies of α-boron in the 2–6 GPa pressure range using MAX80 multianvil system and energy-dispersive synchrotron X-ray powder diffraction at HASYLAB-DESY47 have shown that the phase transformation of α-B12 into β-B106 occurs at temperatures above 1600 K (Fig. 1a & b). It passes noticeably higher the stability line predicted in Ref. 42 and shows strongly non-monotone temperature dependence in the 1600–1800 K range. At higher temperatures, 2400–2500 K, all lines of β-B106 disappear due to the melting (Fig. 1a & c).

Two opinions are know in the literature on the thermodynamical stability of α-B12, i.e. (i) at low and moderate pressures it has its own domain in the p-T phase diagram42, (ii) below 8 GPa β-B106 phase is stable allover the p-T region up to the melting temperature4348. Since the α β equilibrium temperature obtained from ab initio calculations4249 is lower by ~ 800 K than the onset temperature of the α-to-β transformation (Fig. 1b), one cannot make a clear conclusion whether the crystallization of α-boron under these conditions is equilibrium or kinetically limited.

Very recently an attempt to study the vicinity of the α–β–γ triple point has been made by Parakhonskiy et al.50. Although the data on the growth of single-crystal α-B12 in metal systems is interesting by itself, the authors' interpretation of the results in terms of equilibrium thermodynamics is quite contradictory. From one side, at temperatures of 1400–1600 K the numerous metastable phases crystallize1251, especially, on the timescale of Ref.50. Moreover, the extrapolation of the α β equilibrium line reported in Ref.50 down to 1 MPa contradicts the crystallization of α-B12 from Pt melt at higher (at least, 200 K higher than in Ref.50) temperatures45, i.e. in the domain of stability of β-B106 phase. At such low temperatures and for such rigid structures the formation enthalpy of α-B12 from β-B106 (ΔHα/β) is more important for the estimation of thermodynamic stability of a phase than the fact of single crystal growth, which can occur outside the domain of thermodynamic stability, especially, in the case of covalent cage structures (e.g. Si clathrates52). Although Parakhonskiy et al.50 believe that their approach is similar to the classical study of the graphite diamond equilibrium, the close (in contrast to graphite/diamond) structural relationships between boron allotropes410 can result in the metastable crystal growth of boron according to the Ostwald rule of stages. Moreover, the ability of boron to host transition metals (Pt, in particular) in the structure53 was completely neglected by Parakhonskiy et al.50, while it can influence the mutual stability of various structures and the crystallization order (for example, metastable nitrogen-doped t-B52 was observed in the B–BN system already at 5 GPa202229).

To perform thermodynamic analysis of relative stabilities of α-B12 and β-B106, we fitted experimental heat capacities to the Holzapfel equation5455 (Fig. 2a, Tab. 128565758). The estimates for different contributions into total Gibbs energy ΔGα/β (solid curves) are presented in Fig. 2b as dashed curves: (1) thermal due to the difference in heat capacities, and (2) configurational due to non-zero configurational entropy of β-B106149. Taking into account only thermal part would lead to astonishing result: the α-B12 phase should be stable either at high temperature or in the whole temperature range with very narrow (500–750 K) stability domain of β-B106. This fact contradicts to all known experimental data on phase relationships between these two boron allotropes and cannot be overcome by variation of ΔHα/β value. Only including the configurational entropy (the lower estimate based on structural data from Ref.59 is ~ 0.5 R Jmol−1K−1)49 allows us to obtain reasonable ΔGα/β dependencies (Fig. 2b, solid curves). The impact of configurational entropy on thermodynamic properties is remarkable only in the case of β-B106 (standard state): α-B12 and γ-B28 have zero configurational entropy, while that of pseudo-cubic t'-B52 is "hidden" in the corresponding ΔSX/t' (quasi-constant at high temperatures) values describing the X t' equilibria.

Figure 2

(a) Heat capacity of α-B12 and β-B106 at ambient pressure. The triangles show the experimental data565758, while solid lines represent the data fit to the Holzapfel equation5455. (b) Estimate for the formation enthalpy (Hα0K − Hβ0K) of α-B12 at ambient pressure (β-B106 was considered as a standard state). The black dashed line shows the difference between Gibbs energies of α-B12 and β-B106 due to the non-zero configurational entropy of β-phase, i.e. T Sβ0K (Sα0K = 0, Sβ0K ~ 0.5R)49. The green dashed line represents the difference between thermal contributions (due to the difference in heat capacities) into the Gibbs energy of α-B12 and β-B106. Blue solid line indicates the zero energy level chosen as Hβ0K. The highest and lowest temperatures of the synthesis of well-distinguishable α-B12 crystals allow one to suggest that Tα–β should be between 933 K50 and 1500 K51 (red rectangle area). Two red curves show the corresponding Gibbs energies of α-B12 (Gα − Gβ) for two different values of (Hα0K − Hβ0K); while these enthalpy values give the temperature limits (933 K50 and 1500 K51) for the α β equilibrium (Gα(T) − Gβ(T) = 0). The mean value of (Hα0K − Hβ0K) is therefore about −4.5 kJ/mol if crystallization of α = B12 occurs at equilibrium conditions between these two temperatures. (c) Tentative low-temperature part of the boron phase diagram. The α β equilibrium line (dashed) crosses the pressure axis, similar to previous simulations110. The slope (dp/dT)α–β = ΔSα–β/ΔVα–β is positive and has a non-zero value at 0 K due to the configurational enthropy of β-B106. The (dp/dT)α–γ → 0 when T → 0; and at low temperatures the α γ equilibrium line (doted) is parallel to the temperature axis, while at higher temperature it joins the triple point defined by the intersection of the calculated α β (dashed) and experimental β γ (solid) equilibrium lines.

Table 1

The parameters of the Holzapfel equation5455 for heat capacity565758 of α-B12 and β-B106

Parameters*	α-rhombohedral boron	β-rhombohedral boron
θ, K	374(49)	970(10)
C0	2.38613	0.06106
C1	−0.33881	0.12374
C2	1.6714	0.52491
A	0.115	0.125

Now, for the α β "equilibrium" line reported in Ref.50 the ΔHα/β value should be about −4.5 kJ/mol (Fig. 2b) but even this "lower estimate" for ΔHα/β seems to be too high as compared to both the results of ab initio calculations1049 and experimental value of ΔHamorphous/β57. Thus, the experimental estimation of ΔHα/β (as well as Sconfig. of β-B106) seems to be the crucial point for establishing topology of the boron phase diagram in the vicinity of the α–β–γ triple point and at lower temperatures. So far, the α β and α γ equilibrium lines reported in Ref.50 may be only considered as overestimated upper boundaries of the thermodynamic stability domain of α-B12. Moreover, the most recent experimental1 and theoretical10 studies clearly indicate that the α β equilibrium line crosses the pressure axis, in contrast to Ref.50.

Fig. 2c shows the tentative low-temperature part of the boron phase diagram. The α β equilibrium line intersects the pressure axis, similar to previous simulations110. The slope (dp/dT)α–β = ΔSα–β/ΔVα–β is positive and has a non-zero value at 0 K due to the configurational enthropy of β-B106. Our lowest estimate for high-temperature Sconfig. of β-B106 is ~ 0.5 R that is consistent with high-temperature instability of α-B12. According to Monte-Carlo simulations (see Ref.1 and references therein), the low-temperature Sconfig. should be ~ 5 times smaller than that at high temperatures. Thus, at low temperatures the α β equilibrium line should have the slope (dp/dT)α–β ≈ (0.1 R J mol−1 K−1)/(0.153 cm3)1859. As for the (dp/dT)α–γ, it tends to 0 when T → 0; and at very low temperatures the α γ equilibrium line should be parallel to the temperature axis. The low-temperature part of the diagram thus constructed seems to be the most thermodynamically consistent with all available experimental data and ab initio calculations, and is very close to that previously reported by Oganov et al10.

Melting curve of β-B106

The experimental points on the β-B106 melting, both our experimental data and the previously reported values43446061 are presented in Fig. 3. The disappearance of the diffraction lines of crystalline boron is accompanied by appearance of a broad halo typical for a liquid phase unambiguously points to the melting (Fig. 1a). Even slight decrease of the melt temperature leads to crystallization of β-B106, indicative of equilibrium melting point and not metastable melting. At 5.45 GPa the melting temperature was found to be 2440 K. This is an intermediate value between the reported melting temperature of 2350–2370 K6061 at ambient pressure and 2480 K at 7.7 GPa43. The experiments at pressures below 2.5 GPa have always led to blowouts and have not allowed establishing the low-pressure part of the melting curve.

Figure 3

Phase diagram of boron.

The symbols show the experimental data. The solid lines represent the equilibria between different allotropes obtained by fitting the experimental points, as well as the unique melting curve of t'-B52 thermodynamically consistent with other phase equilibria. Down and up triangles show solid and melted β-B106. Squares and circles correspond to γ-B28 and t'-B52, respectively, recovered after quenching down to normal conditions. Small open symbols show literature data13. Color guide for symbols and structures: blue – β-B1069, red – γ-B2810, wine – t'-B5246, black – liquid B (structural units of amorphous boron75).

According to Ref.62, ΔHL/β = 50.2 kJ/mol at ambient pressure. For our calculations, the melting temperature of boron at ambient pressure has been chosen as 2360 K (the mean value of 2350 K60 and 2370 K61). Since ΔVL/β is the value that strongly influences the slope of the melting curve, we have used it as a fitted parameter. The zero-pressure atomic volume of β-boron has been calculated by extrapolation of the data on thermal expansion63. A good fit has been obtained with ΔVL/β of 6.7% of the molar volume of crystalline boron at meting temperature which is close to the ~ 5–10% volume change reported in Refs.61, 64, 65 (Tab. 2) and is noticeably lower than the 20% value taken for the calculation of the melting curve in Ref.42. Thus, the experimental melting curve of β-boron up to 8 GPa may be described by the linear equation

Table 2

Thermodynamic data on phase transformations in boron at pressures up to 25 GPa and temperatures between 1500 and 3000 K

Y/X	ΔHY/X (kJ/mol), ΔSY/X (kJ/mol K) or ΔVY/X (cm3/mol)	Fitting option	Experimental values and/or ab initio predictions
L/β	ΔH = 50	Fixed to experimental value	ΔHmelt = 50.24362, 48.9370
	ΔV = 15.5 ΔH/2360.0 = 0.32	Constrained to experimental melting curve and experimental ΔHL/β	Vmelt = 4.986465 to 5.2261; Vβ = 4.7663 (V300K = 4.5759)ΔV = 0.22 to 0.46
	ΔS = ΔH/2360.0 = 2.1·10−2
γ/β	ΔH = 2.5	Fixed to ab initio prediction	ΔH0K = 2.510
	ΔV = − 707.8 ΔH/4534.0 = − 0.35	Constrained to experimental transformation curve and ΔHγ/β	ΔV300K, 1 MPa = 4.2510 − 4.5759 = − 0.32
	ΔS = − ΔH/4534.0 = − 5.5·10−4
t'/β	ΔH = 8.8	Fitted
	ΔV = − 160.2 ΔH/3755.7 = − 0.38	Constrained to experimental transformation curve and fitted ΔHt'/β	ΔV300K, 1 MPa = 4.2746 − 4.5759 = −0.30
	ΔS = ΔH/3755.7 = 2.3·10−3
t'/γ	ΔH = ΔHt'/β − ΔHγ/β = 6.3	Adjusted to fitted ΔHt'/β and ΔHγ/β
	ΔV = 5.1 ΔH/2176.0 = 0.01	Constrained to experimental transformation curve and adjusted ΔHt'/γ	ΔV300K, 1 MPa = 4.2746 − 4.2510 = 0.02
	ΔS = ΔH/2176.0 = 2.9·10−3
L/t'	ΔH = ΔHL/β − ΔHt'/β = 41.2	Defined by ΔGt'/β and ΔGL/β
	ΔV = ΔVL/β − ΔVt'/β = 0.7
	ΔS = ΔSL/β − ΔSt'/β = 1.9·10−2

Equilibrium line between β-B106 and γ-B28

Our experiments at pressures above 10 GPa have shown that samples quenched from 1600–2000 K contain only pure γ-B28, while the samples quenched from 8 GPa contain only β-B106. These results well agree with experimental data of other groups1350, and the β γ equilibrium line that fits both our and literature data can be described by equation

Equilibrium line between γ-B28 and t'-B52

Above 10 GPa the temperature seems to have more impact on the recovered boron allotrope as compared to the pressure. At 15–20 GPa, the samples quenched from 2500 K contained only mixture of γ-B28 and pseudo-cubic t'-B5246, while below 2200 K t-B52 (as well as α-B12) could be observed as intermediate phase prior to crystallization of γ-B28 in the stability domain of the latter1213. At such high temperatures the kinetic factors do not play a decisive role anymore due to the intense diffusion. Thus, above 10 GPa pseudo-cubic t'-B52 is stable at high temperatures, while γ-B28 – at low and moderate temperatures.

All these results allowed us to define the domain of thermodynamic stability of γ-B28 in the phase diagram (Fig. 3). The triple point between β-B106, γ-B28 and t'-B52 should be located at ~ 9.6 GPa and ~ 2230 K, while the equilibrium line between γ-B28 and t'-B52 is described by the equation Only dense t'-B52 phase can explain such a low pressure slope of the t' γ equilibrium line, contrary to other known tetragonal phases, t-B52 and t-B19246.

Equilibrium line between β-B106 and t'-B52

The position of the β-γ-t' triple point, experimental data on phase stability of β-B106 and the lowest pressure of t'-B28 formation (in some experiments this phase has been recovered at ~ 7.7 GPa) allow one to draw the β t' equilibrium line (Fig. 3) which has negative pressure slope and follows the equation This equilibrium line intersects the β-B106 melting curve at ~ 8.0 GPa and ~ 2490 K that is the triple point between β-B106, t'-B52 and liquid.

Discussion

Among all boron allotropes mentioned above, the experimental thermodynamic values are known only for α-B12 and β-B106, at least at ambient pressure. That is why it seems reasonable to use β-B106 as standard state with known thermodynamic potentials, while for L, t' and γ phases (X or Y) the corresponding Gibbs energy should be corrected by a value of ΔGX/β defined by three parameters, i.e. ΔHX/β, ΔVX/β and ΔSX/β, independent of pressure and temperature in the first approximation. These parameters, if not known from experiment, can be obtained by fitting the experimental p-T data to theoretical isopotential X Y equilibrium lines (i.e. ΔGX/Y(p,T) = ΔGX/β(p,T) − ΔGY/β(p,T) = 0).

Since the majority of experimentally observed equilibria between boron allotropes can be presented by straight lines, we did not complicated our analysis with exact formulae containing integrals. For example, p·ΔVX/β and ∫ΔVX/β dp differ only by ~ 5% (particularly, due to the close compressibilities of all boron allotropes)66676869. Thus, the fitted parameters should be interpreted as average values over p-T domain of experimental data.

In order to describe the high-temperature part of the boron phase diagram, one should know at least 9 parameters e.g. 3 expressions for ΔGX/β each containing 3 parameters, i.e. ΔHX/β, ΔVX/β and ΔSX/β. These dependences follow the equation For example, if one have such parameters for X = L, t' and γ, two other Gibbs energies ΔGt'/γ and ΔGL/t' will be related to them by simple general relationship and will give us thermodynamically consistent equilibrium lines (ΔGX/β (p,T) = 0). In order to establish these 9 parameters, one needs at least 9 independent experimental values characterizing those equilibria. Each experimental equilibrium line gives us two parameters: ΔH/ΔS and ΔV/ΔS according to the equation that is a good approximation at high temperatures.

Four experimental equilibrium lines described above, i.e. L β, β γ, t' β and t' γ give us 8 independent parameters. The ninth one we have chosen as experimental value of enthalpy of boron melting ΔHL/β (see Tab. 2) quite reproducible in independent experiments at ambient pressure436270.

The fitted and fixed values of the thermodynamic parameters are given in Tab. 2 (data from Refs.10, 43, 46, 49, 59, 61, 62, 63, 64, 65, 70). One can observe a reasonable agreement of the fitted parameters with estimates for ΔVL/β and ΔVt'/β based on volume measurements at high temperatures, as well as with ab initio predicted ΔHγ/β value at 0 K. All this additionally indicates self-consistency and uniqueness of thermodynamic equations used, as well as real physical meaning of the fitted parameters.

The thermodynamic parameters of all mentioned above equilibria give the unique melting line of t'-B52, i.e. The parameters that define corresponding Gibbs energy are given in Tab. 2.

Figure 3 shows the phase diagram that represents all equilibria described above. At low pressures, two rhombohedral phases are stable, α-B12 and β-B106. These allotropes have related structures4, rhombohedrally distorted fcc packing, produced by individual B12 icosahedra at low temperatures, and more complicated clusters at high temperatures. At pressures above 10 GPa, two other phases become stable, namely, γ-B28 and t'-B52. Both of them have structures similar to distorted NaCl (orthorhombic and tetragonal distortions, respectively)1046. It should also be noted that the common feature of two high-temperature phases (β-B106 and t'-B52) is a certain degree of intrinsic structural disorder (partially occupied Wickoff positions), which gives to them non-zero configurational entropy already at 0 K and, subsequently, results in their stability at high temperatures.

Finally, the equilibrium phase diagram of boron has been constructed at pressures up to 20 GPa and temperatures up to 2500 K (Fig. 3). It has been experimentally proved that at least four boron phases, i.e. β-B106, γ-B28, t'-B52 and liquid boron, have p-T domains of thermodynamic stability. Two triple points have been established, i.e. the first between β-B106, t'-B52 and liquid at ~ 8.0 GPa and ~ 2490 K; and the second between β-B106, γ-B28 and t'-B52 at ~ 9.6 GPa and ~ 2230 K. Thermodynamic analysis leads to the unique melting curve of t'-B52 that is consistent with other experimentally established phase equilibria.

Methods

Multianvil experiments

Quenching experiments to 20 GPa were performed using 6–8 type large-volume multianvil systems with octahedral pressure assemblies at Laboratoire des Sciences des Procédés et des Matériaux (LSPM–CNRS) and Bayerisches Geoinstitut (BGI). The experimental details are described elsewhere7172. Pressure and temperature have been either directly measured (thermocouples and p-T phase transitions of reference materials) or estimated from previously obtained p-T calibration curves; in all cases, the uncertainties were estimated to be about 1 GPa and 50 K, respectively. Samples were gradually compressed to the desired pressure at ambient temperature, and then the temperature was increased with a rate of about 300 K/min up to the desired value. After heating for 5–10 min, the samples were quenched by switching off the power and then slowly decompressed.

In situ X-ray diffraction

In situ experiments to 7 GPa were carried out using multianvil X-ray system MAX80 at beamline F2.1, DORIS III (HASYLAB-DESY). The experimental setup has been described elsewhere47. Energy-dispersive X-ray diffraction data were collected on a Canberra solid state Ge-detector with fixed Bragg angle θ = 4.555(3) using a white beam collimated down to 100 × 100 μm2.

The sample temperature up to 2200 K was measured by a W3%Re–W25%Re thermocouple. The correction for the pressure effect on the thermocouple emf was made using the data of Li et al73. Above 2200 K the power – temperature calibration curve was linearly extrapolated to the high-temperature region (up to 2600 K). Pressures at different temperatures were found from the p-V-T equation of state of highly ordered (P3 = 0.98 ± 0.02) graphite-like hexagonal boron nitride74.

The samples were gradually compressed to the required pressure at ambient temperature and then diffraction patterns were collected at the stepwise (~50 K) temperature increase. With the storage ring operating at 4.44 GeV and 150 ± 50 mA, diffraction patterns were collected for 1 min in real time. After heating, the samples were quenched by switching off the power, and then the pressure was slowly released down to ambient.

Ex situ X-ray diffraction

The recovered samples were studied by conventional powder X-ray diffraction using G3000 TEXT (Inel) diffractometer (Bragg-Brentano geometry) employing CuKα1 radiation. The goniometer was aligned with high purity silicon (a = 5.431066 Å) and the standard sample of LaB6 (a = 4.15695 Å). Synchrotron X-ray powder diffraction measurements (λ = 1.10347 Å) have been performed at beamline I711, MAX II (MAX-lab); Debye-Scherrer geometry with rotating quartz capillary was used. Unit cell parameters, the size of the blocks of coherent scattering and strains were derived from the LeBail profile analysis performed using the GSAS program.

Author Contributions

V.L.S. and O.O.K. equally participated in experimental work, data analysis and writing the manuscript.