Aberrant electronic and structural alterations in pressure tuned perovskite NaOsO3

The perovskite NaOsO3 has a metal–insulator transition at temperature 410 K, which is delicate, intriguing, and provokes a lot of debate on its nature. Our combined electrical resistance, Raman, and synchrotron x-ray diffraction experiments show that the insulating ground state in this osmate endures under high pressure up to at least 35 GPa. In this pressure range, compression reveals hidden hysteretic resistance properties with a transient metallic state near 200 K, manifested three electronic character anomalies (at 1.7, 9.0, and 25.5 GPa), and a structural transition to the singular polar phase (at ~18 GPa). We distinguish NaOsO3 from the regular crystallographic behavior of perovskites, though the electrical specificities resemble iridates and nickelates. The theoretical first-principle band structure and lattice dynamics calculations demonstrate that the magnetically itinerant Lifshitz-type mechanism with spin–orbit and spin–phonon interactions is responsible for these pressure-induced changes. Our findings provide another new playground for the emergence of new states in 5d materials by using high-pressure methods.


INTRODUCTION
The composition of NaOsO 3 in the primitive cubic form has been known since 1974 1 . The orthorhombic perovskite NaOsO 3 was only obtained in 2008 2,3 . The new orthorhombic phase attracted a lot of attention, because of its rare metal-insulator transition (MIT) which coincides with a long-range commensurate threedimensional antiferromagnetic (AFM) ordering at a rather high Néel temperature T N = T MIT = 410 K 3-6 . Initially, a Slater-like mechanism was proposed to explain this unusual MIT behavior 4 . Subsequent insights have revealed that this transition is better interpreted by a relativistic magnetic Lifshitz mechanism where MIT is driven by magnetic fluctuations and spin-orbit-induced renormalization of the electron correlation 7,8 .
At high temperatures (T > 410 K), NaOsO 3 shows paramagnetic Fermi liquid behavior, but with decreasing temperature Fermi surface undergoes continuous reconstruction prompted by the changes in the amplitude and direction of the spins which eventually causes a narrow band gap of~0.1 eV at ambient conditions. On the other hand, a huge interaction between spins and lattice vibrations was also observed in NaOsO 3 where pressure or strain directly influences the OsO 6 octahedra 9 . The anomalous spin-phonon behavior was linked to the extended 5d orbitals which relate to the magnetic structure through the Os-O-Os superexchange interactions. In this regard, pressure becomes a favorable tool to induce many new features to emerge, because it can substantially affect the Fermi surface reconstruction 10,11 and, in turn, unpredictably vary the band gap. For this reason, we have initiated a combined experimental and theoretical study to reveal the electronic and structural properties of NaOsO 3 under high pressure.
In general, for perovskites, the relative stability in different modifications can be estimated with a wide range of methods where the tendencies for plausible symmetry transformations can be predicted by considering possible octahedral distortions [12][13][14][15] . Hence, the orthorhombic distortion in NaOsO 3 arises from the concerted rotation of the OsO 6 octahedra in the a-c plane, coupled with a tilting about the b-axis. According to the mentioned theoretical sources, an orthorhombic configuration should be lowered to a monoclinic one if more distortion occurs. The sequence of structural transitions is predicted in theory as: cubic (Pm3m) → orthorhombic (Pnma) → monoclinic (P2 1 =c). Another indicator for possible high-pressure trends is the ratio of cation coordination polyhedra volume (V A /V B in ABO 3 ) that combines both the the tilt angle and the deformation of the A and B-site polyhedra. An investigation on existing perovskite structures (~1300 perovskites) shows that the orthorhombic (Pnma) can reach a polyhedral volume ratio as low as 3.8 16 . Among perovskites, there are very few materials reaching values of V A / V B < 4, which is a distortion factor and also an indicator for a perovskite to undergo the so-called 'post-perovskite' transition 17 .
In the orthorhombic Pnma arrangement NaOsO 3 has a value of 4.6 for the V A /V B ratio 15,16 . Increasing pressure is known to stimulate the rotation of the BO 6 octahedra in ABO 3 perovskites, and therefore decreases the polyhedral volume ratio. Hence, the perovskite NaOsO 3 at 16 GPa and 1135 K was found to stabilize in a post-perovskite form 18 . Strikingly, we discovered that, without increasing the temperature, the orthorhombic phase of NaOsO 3 could transform into a polar perovskite phase with a space group of Pna2 1 , which violates the general trend of the perovskite oxides and has only been detected in PbRuO 3 (Pbnm → Pbn2 1 ) 19 and LiOsO 3 (R3c → R3c) 20 so far. The latter two materials are nonmagnetic, while NaOsO 3 has obvious magnetic properties 2-9 and preserves the band gap during the transition with a slow reduction of the MIT towards lower temperatures (T < 410 K). In addition, compression of NaOsO 3 manifests much more complicated emergent structural and electronic phenomena, which could be only characterized by using multiple techniques, such as electrical transport, Raman, synchrotron x-ray diffraction techniques (XRD), and first-principles calculations that are applied in this study.

Electrical anomalies
As already mentioned, the G-type long-range commensurate AFM order occurs at the same temperature as the continuous MIT. The neutron diffraction measurements confirmed that there is no symmetry change through the MIT at ambient pressure 4,21 . Hence, a purely electronic band gap opening mechanism was considered. The most fundamental property, resistivity, measured at ambient pressure indicated a very clear distinction between metallic behavior at T > 410 K and insulating behavior at T < 410 K 3 (see the inset of Fig. 1b), while more detailed analysis indicates that there are actually three regions in the resistivity ρ(T) curve: metallic (T > 410 K), bad-metal (BM) (T = 30-410 K), and insulating (T < 30 K) 7 .
The BM state was distinguished from the change of ρ taking place around 30 K, where the insulating indirect band gap closes and NaOsO 3 starts to develop the Fermi surface characterized by several hole and electron pockets with a direct gap (pseudogap). The pseudogap continues to decrease with increasing temperature and eventually closes at T MIT . These changes are also observed in our high-pressure electrical transport measurements (see Fig. 1). Furthermore, the warming and cooling curves slightly deviate, forming a narrow thermal hysteresis loop below MIT (Fig. 1a). The hysteresis is progressively attenuated upon pressure but eventually disappears at about 18 GPa. Intriguingly, there is no hysteresis below 30 K, where the deep insulating state is developed. One more anomaly revealed under pressure emerges at around 200 K where the R(T) curve bends down. It is likely that this feature exists even at ambient pressure, because the resistivity in the BM regime does not show straight linear dependence on temperature ( Supplementary Fig. 1) and, therefore, is naturally enhanced by pressure. As highlighted in the Arrhenius plot, this feature at 0.63 GPa manifests a kink which develops into the metallic state with increased pressure (see Fig. 2). The measurements at higher pressures clearly show that the metallic state unfolds most at 10.60 GPa and reenters the original electronic state at around 18.70 GPa, which is on the border of the structural transition. We refer to this anomaly as the BM-M-BM transition because of the following change: BM → M → BM. The origin of this transition is discussed in the theoretical modeling section.
In terms of the MIT, the pressure-dependent electrical transport measurements indicate that the metallic state extends to the lower temperatures very slowly (Fig. 1b). As one can observe, the T MIT scales almost linearly upon pressure. At around 32 GPa, the MIT becomes much broader, but can still be identified. Importantly, up to this pressure, NaOsO 3 preserves the insulating ground state. It is known that nickelates such as PrNiO 3 and NdNiO 3 exhibit similar MITs from the paramagnetic metal to the AFM insulator that are electronically driven as well 22 . The suppression of the first-order MIT in PrNiO 3 is very fast, within 0-1.3 GPa range, but ρ(T) curves manifest a highly broadened thermal hysteresis loop under pressure. In comparison, pressurestimulated hysteretic properties in NaOsO 3 are relatively small and likely to be hard to distinguish at ambient pressure. However, this questions if MIT is really the second-order type that was initially assigned 3 . On the other hand, this feature may emerge for various b a Fig. 1 Suppression of the temperature-dependent MIT in NaOsO 3 . a Temperature dependence of the resistance at various pressures. The inset emphasizes the hysteresis evolution over pressure at T = 295 K. b T MIT vs. pressure. The inset for MIT observed at ambient pressure is taken from ref. 3 . The inset demonstrating MIT at 32.5 GPa represents data from this work. The resistance at high pressure becomes low and T MIT is washed. The error bars represent the uncertainty in determining T MIT and P. c The electrical resistivity measurement during sample compression at fixed temperature of T = 295 K. The inset shows the actual set-up for the resistivity measurement between a pair of diamond anvils: four thin platinum probes attached to the sample, and the ruby sphere loaded for pressure determination.
reasons such as the spin canting effect, inhomogeneity, or defects, and, therefore, could be addressed in detail in the future. Meanwhile, precise compression of NaOsO 3 at a fixed temperature of T = 295 K revealed four anomalies of which none are attributable to the MIT (Fig. 1c). This underlines the uniqueness of this sample and discerns it from nickelates, which have a short track to the metallic state. The first anomaly in compression appears at a sufficiently low pressure~1.7 GPa; the resistivity drops most sharply here. The second anomaly manifests at 9 GPa as deceleration in resistivity with increased pressure. The third, at 18.8 GPa, is related to the structural change, as explained in the next section. The last and fourth anomaly at 25.5 GPa indicates the limit where the resistivity value no longer decreases.
Structural transformation under high pressure Although orthorhombic perovskite NaOsO 3 (Space group Pnma, no. 62. Lattice parameters a = 5.38420(1) Å, b = 7.58038(1) Å, and c = 5.32817(1) Å) exhibits pronounced changes in its magnetic and electronic properties, the structure is reliably stable over the wide range of temperatures from 2 to 480 K 2-5 and up to 16 GPa pressure investigated so far 18  The mentioned study on NaOsO 3 also showed that pentavalent post-perovskite structure may be obtained even if all the criteria normally used to anticipate such transitions are not fulfilled 18 . Our pressure-dependent Raman measurements support the fact that the crystal symmetry does not change up to 16 GPa at room temperature and indicates that further pressure increase causes structural transformation to a different symmetry. Figure 3 shows the Raman spectra observed in the low-frequency range at different temperature and pressure conditions. Here, seven optical peaks are clearly visible at ambient conditions. The magnon excitations identified using resonant inelastic x-ray scattering (RIXS) can be observed in the Raman spectra in many cases as well. A peak at 400 cm −1 seen in Fig. 3a could represent a singlemagnon excitation observed in the NaOsO 3 RIXS experiments 8,23 . The Raman spectra of NaOsO 3 at different pressure-temperature conditions. a Raman spectrum at ambient pressure, and room temperature acquired for 120 s and five accumulations using a 488 nm wavelength. b Selected Raman spectra for the different high-pressure states of NaOsO 3 . The measuring time of 200 s were kept constant for all spectra. c Raman spectra at ambient pressure and different temperatures revealing three components α, β, and γ hidden within a broad peak at nearly 400 cm −1 . Data acquired for 120 s and five accumulations using a 532 nm wavelength. d Pressure evolution of the phonon frequencies. Here the dashed lines indicate notable changes associated with the pressure-induced phenomena. The Raman mode marked as no. 8 is the weak mode, marked with an asterisk in b, and is possibly the reappearance of mode no. 3. Fig. 2 Arrhenius plots of the electrical resistance R at fixed pressures. A small kink observed at 0.63 GPa evolves when pressure reaches a maximum at 10.60 GPa, and disappears for pressures more than 18.70 GPa. Here, the metallic regime is considered in the pressure range from 0.9 to~18 GPa.
By reducing the temperature, we show that this peak exhibits three components α, β, and γ (see Fig. 3c). Here, β shows the most pronounced changes, but the intensity does not change much and, therefore, it can be assumed that these three components are regular Raman modes rather than the manifestation of a magnon. Overall, the optical data is in good agreement with the electrical transport measurements confirming four anomalies at similar pressures. Ambient crystal symmetry should be maintained up to 18 GPa, even though some of the weak modes disappear in the background during compression (see Fig. 3b, d). The phonons substantially rearrange and increase in number at 18 GPa, signaling the onset of a structural change. Then, at about 26 GPa, the continuous large-scale reduction in intensity is observed as the pressure increases. Finally, the Raman modes almost vanish at 35 GPa, indicating that sample is approaching a metallic state, that is the MIT. It should be noted that these observed changes are reversible when the sample is released to the ambient pressure.
To determine the crystal symmetry in the new structure above 18 GPa, we performed synchrotron powder X-ray diffraction (see Fig. 4 and Supplementary Figs. 2-4). The new structure seems to fit well to the space group of Pna2 1 (orthorhombic and noncentrosymmetric) which agrees with the Raman results. We notice the analogy of this rare transition with the one observed in PbRuO 3 samples, whose high-pressure phase is claimed to have the most distorted octahedra in a perovskite structure 19 . The PbRuO 3 demonstrated that applied pressure can substantially shorten A-B range and distort the octahedra thereby compacting ABO 3 perovskite structure. The perovskite structure proves flexible through the A-B bonding and bypasses the collapse into the postperovskite structure. Our discovered high-pressure structural transition for NaOsO 3 also avoids a possible collapse into the post-perovskite structure and manifests more distorted octahedra in the polar perovskite phase. Besides, the polar phase exhibits bad metal characteristics for both PbRuO 3 and NaOsO 3 . In addition, a non-polar to polar transition is also observed in metallic LiOsO 3 20 . The recent study shows that the transition temperature to the polar phase of metallic LiOsO 3 can reach to room temperature at 16 GPa 24 . This structural transition is closely linked to the changes in the local bonding environment where Li ions form bonds with out-of-plane oxygen atoms forming the LiO 6 octahedra. This reaffirms that pressure-induced octahedral deformations play a key role in the formation of the polar phase.
A remark concerning the symmetry of the new structure above 18 GPa at room temperature is in place here. From our Raman results (Fig. 3), we are certain that NaOsO 3 experiences a structural change above 18 GPa. When the pressure is increased, the Raman spectra in particular demonstrate the enhancement of the number of phonons and the pressure-induced-splitting of phonon mode (Fig. 3b). This is consistent with a shift in symmetry toward lower space groups. The reduction of the space group in the original orthorhombic crystal system is the most appropriate assertion, given that the changes in the diffraction patterns are only minor alterations from the new high-pressure phase. It should be noted that the space group change in the powder XRD measurements is subtle (Fig. 4a), because the patterns between Pnma and Pna2 1 symmetries are in many cases not remarkably different 19 . The transition is very smooth giving a tiny volume contraction (<0.1%) as well as slight fracture in the lattice parameters (Fig. 4b, c). A Birch-Murnaghan third-order equation of state fitted to the volume versus pressure shows a bulk modulus increase from K 0 ¼ 152:0 GPa to K 0 ¼ 200:0 GPa for the Pnma and Pna2 1 phases, respectively (see Fig. 4b). This increase is typical for such samples. For instance, the change of bulk modulus in PbRuO 3 is similar:  45 . c Pressure dependence of a/a 0 , b/b 0 , and c/c 0 ratios. Here, a 0 and c 0 are the lattice constants at ambient pressure. Inset is the zoomed area for lattice parameter change at the smaller pressure. Figure 5 shows the phonon dispersion curves for NaOsO 3 in the Pnma phase at ambient pressure. No imaginary frequencies are obtained, which indicates the stability of the structure at the given condition. It is worth noting that the inclusion of the G-AFM ordering is needed to obtain stable phonon data, showing the importance of the spin ordering in stabilizing the structure. The phonon bands are clearly divided into two groups: the lower frequency (0-52 meV) and higher frequency (75-97 meV). The results of partial phonon density of states (DOS) indicate that all atoms contribute to the lattice vibrations at lower frequency range, while vibrations of O atoms are dominant at higher frequency range (Fig. 5). Symmetry-mode analysis of the atomic displacements predicts 24 Raman active modes in NaOsO 3 at the Brillouin zone center 25 : 7A g + 5B 1g + 7B 2g + 5B 3g . All these modes were observed in our DFT calculations as well. Table 1 represents a comparison of the frequencies of the calculated and observed Raman-active modes, showing overall good agreement between experiment and simulation, in particular at low temperature.

DISCUSSION
The structural evolution under pressure is more problematic. The pressure evolution of the simulated XRD spectra agrees very well with the measured XRD spectra [compare Supplementary Fig.  2a, d], but, in contrast to the Raman results, our DFT calculations do not discern any marked symmetry change even up to high pressure (90 GPa) as evinced by the pressure-volume curves shown in Supplementary Fig. 5. The phase with initial Pna2 1 symmetry is found to be unstable and finally reduces to the Pnma symmetry after structural relaxations (see Supplementary Fig. 5). The apparent disagreement might originate from thermal effects, not considered in our DFT analysis. Another possible reason for the experimentally observed structure change above 18 GPa could be the existence of a certain degree of anisotropy in the applied pressure. However, the nonpolar-polar phase transition is normally disfavored by temperature. On the other hand, the LiOsO 3 case unambiguously demonstrates the enhanced stability of the polar phase under high-pressure 20 . Thus, in contrast to conventional ferroelectrics, pressure can give rise to a polar phase or enhance the existing transition temperature in these systems.
Moving to the electronic properties, Fig. 6 shows our DFT calculated band gap and magnetic moment evolution as a function of pressure. Upon pressure, the indirect gap is closed at around 20 GPa, while the direct gap survives even up to very high pressure, indicating the insulator to BM transition upon pressure. As expected, the magnetic moment is suppressed by the pressure, but the suppression rate is relatively slow and the system remains magnetically ordered even at~40 GPa. The detailed evolution of the band structure as a function of pressure is shown in Fig. 7. Figure 8 is a combined temperature-pressure electronic and structural phase diagram of NaOsO 3 based on our experimental and theoretical results, which shows the different types of transitions occurring in the system when pressure and temperature are applied: insulator-to-bad metal, BM-to-M, the anomalous metal island in the BM region, and the subtle non-polar to polar structural transition. At low temperature the system remains insulating up to a certain critical pressure (~20 GPa in DFT) and   then transforms into a bad metal due to the closing of the indirect gap. In this pressure range the valence and conduction bands are still separated by a direct gap (Figs. 6 and 7). This direct gap closes at very large pressure (>90 GPa in DFT), indicating that the evolution of the electronic properties upon pressure share similarities with the temperature-induced band gap closing process through a magnetically itinerant Lifshitz mechanism 7 . As the temperature increases, a metal island emerges in the BM region under increasing pressure. Considering that the pressureenhanced structural distortion in NaOsO 3 favors an insulating state 26 , we speculate that depending on the distortion track, reentrance transitions are possible and we believe that the unique existence of the BM state enables such manifestations to take place. On the other hand, the resistance might contain specific anisotropic properties that are not considered in this work. Finally, above a critical temperature T MIT NaOsO 3 becomes metallic. T MIT is less affected by pressure, and the experimentally traced rate of T MIT decrease is quite slow~0.36 GPa/K. We note that such sluggish process of metallization during compression among 5d systems (especially iridates) is common with the presence of strong SOC, which usually sustains a rigid band gap [27][28][29][30][31] . In NaOsO 3 with the t 3 2g electronic configuration, the SOC formally should not play a dominant role because the resulting L eff = 0 state yields a vanishing orbital moment, yet it is found to be responsible for the weakening of the electron-electron correlation 7 , magnetic anisotropy, and large spin wave gap formation 32 . As the MIT is considered to be the SOC-driven for the most part, its evolution during compression is likely governed by the same source. In this aspect, NaOsO 3 is akin to iridates. However, a high degree of spin-phonon-electronic coupling was proclaimed as a dominant force in this compound 9 , and, therefore, some features manifest differently. Our combined experimental-theoretical study proves a close spin-phonon-electronic interaction in that all anomalies, including the MIT, are captured in both the electronic and optical spectra and exhibit relatively good accuracy. We found that pressure and temperature both induce a MIT accompanied by the formation of an intermediate BM regime and by the progressive decrease of the local magnetic moment. This process can be assimilated to an itinerant Lifshitz MIT associated  with magnetic fluctuations. This scenario is apparently different from the MIT induced by electron doping where NaOsO 3 undergoes a straightforward transition to a metallic state even at a moderately low doping level 33 , hereby highlighting the complex nature of this compound.
In summary, a comprehensive temperature-dependent electrical transport, Raman scattering, synchrotron XRD, and DFT study has been carried out to understand the effect of external pressure on perovskite NaOsO 3 . A sluggish temperature-dependent MIT was followed up to 35 GPa where T MIT = 410 K was suppressed only to a near room temperature, reiterating the common behavior of iridium-based insulators that are known to have strong SOC and own narrow band gaps in Mott, Mott-Hubbard, or Slater paradigms. However, NaOsO 3 manifests a magnetically itinerant Lifshitz-type transition regardless of external conditions and possesses a high degree of spin-phonon-electronic coupling. The optical and electrical spectra revealed clear anomalies at 1.7, 9, 18, and 25.5 GPa. Furthermore, during compression, the BM regime manifested hysteretic resistance anomalies and a BM-M-BM transition at~200 K during compression. The later anomalies are canceled by the structural change that takes place at 18 GPa from a centrosymmetric Pnma to a noncentrosymmetric (polar) Pna2 1 phase. The high-pressure phase has more distorted octahedra and is unexpected according to commonly used theoretical prediction models. Thus, NaOsO 3 is in line with PbRuO 3 and LiOsO 3 , where pressure favors the polar phase 19,20,24 . However, magnetism is absent in the latter examples and, therefore, NaOsO 3 might be an example that can acquire the unique shape of a polar-magnetic metal.

Sample synthesis
The starting materials Na 2 O 2 (97%, Sigma-Aldrich) and OsO 2 (laboratorymade from 99.95% Os, Nanjing Dongrui Platinum Co., Ltd. by heating at 400°C in a few % oxygen in Argon balance) were sealed in a platinum capsule at a molar ratio of Na:Os = 2:1 with 0.5 mol of NaCl (99.99%, Rare Metallic Co.) per f.u. The procedure was conducted in an Ar-filled glove box. The capsule was heated at 1500°C for 60 min under a pressure of 6 GPa in a multi-anvil-type apparatus (CTF-MA1500P, C&T Factory Co., Ltd, Japan). Then, the pressure was gradually released. The product was rinsed with water in an ultrasonic bath for a few minutes several times, separated mechanically, and then tested to confirm there was no damage by singlecrystal XRD and powder XRD. Small, ≤0.5 mm-sized NaOsO 3 crystals with a silver-purple tint were obtained together with black crystals of the hexavalent Na 2 OsO 4 2,34 in the same batch.
Optical Raman measurements A Mao-type symmetric diamond anvil cell (DAC) with 400 μm culet-sized anvils was used for the Raman experiments. A stainless-steel gasket was precompressed to 35 μm thickness and a hole of 150 μm was drilled to load the sample, a ruby sphere for pressure determination 35 , and neon gas to serve as a pressure-transmitting medium 36 . Initially, we made several attempts to use a methanol-ethanol mixture as a pressure-transmitting medium, but the results were unreliable because the pressure medium reacted with the sample when pressure was applied. The high-pressure Raman spectra were measured on a Renishaw inVia spectrometer with a 488 nm laser wavelength, which gave the best signal. The data collection time was 200 s and a laser power of 8 mW was maintained for each spectrum. The sample was placed in the cryostat under vacuum conditions for the temperature-dependent data at ambient pressure. In this case, Raman spectra were measured on a MonoVista CRS+ system with a 532 nm laser.

Synchrotron x-ray diffraction measurements
The in situ high-pressure powder XRD measurements were carried out in an angle-dispersive mode at the beamline 16-BM-D of the Advanced Photon Source (APS), Argonne National Laboratory. The same type of DAC with the loaded ruby and neon pressure medium was used as the Raman experiments. The experiments were performed two times where the incident monochromatic x-ray beam had the wavelengths of λ = 0.4133 Å and λ = 0.3100 Å. Each time, the wavelength of the x-ray was calibrated using a CeO 2 standard. Diffraction patterns were recorded on a MAR345 image plate and then integrated by using DIOPTAS software 37 . Indexing and Rietveld refinements were carried out in EXPO2014 38 .

Electrical measurements
The electronic transport properties under high pressure and low temperature were investigated via the four-probe electrical conductivity method in a DAC made of CuBe alloy. Pressure was generated by a pair of diamonds with a 400 μm diameter culet. A gasket made of stainless steel was pressed and a hole in the center of the gasket was drilled to the diameter of the diamond culet. A cubic boron nitride (c-BN)-epoxy insulation layer was prepared to protect the electrode leads from the metallic gasket. Finally, a hole of 280 μm diameter was drilled in the center of the c-BN gasket. NaCl powder was used as the pressure-transmitting medium. Four platinum strips were arranged to contact the sample in the chamber. The pressure-dependent resistance measurements were performed using the laboratory designed electrical transport system where a Keithley 6221 current source, 2182A nanovoltmeter, and 7001 switch device were used as the current supply, voltmeter, and voltage/current switcher, respectively. The experiments at different pressure-temperature conditions were performed using the Quantum Design Physical Properties Measurement System. The pressure was measured via the ruby fluorescence method at room temperature, before and after cooling.

First-principles calculations
The computational study was based on the density functional theory (DFT) as implemented in the VASP software 39 . We used the projector-augmented wave method (PAW) 40 and the revised Perdew-Burke-Ernzerhof (PBEsol) functional 41 , which is known to give a better description of densely packed structures than the standard PBE functional 42 . NaOsO 3 with Pnma symmetry with a long b-axis was employed 3 . A comparison with Pna2 1 symmetry with a long c-axis was made to inspect the structural transitions. A fully relativistic setup was used with the inclusion of spin-orbit coupling (SOC) in the framework of the DFT + U 43 (U eff = 1.0 eV) and G-type AFM order. For the energetics and electronic structure calculations, the planewave cutoff for the orbitals was set to 600 eV and Brillouin zone sampling of 10 × 6 × 10 (10 × 10 × 6) k-point grid according to Monkhorst-Pack scheme was used for the Pnma (Pna2 1 ) phase. For the phonon calculation, we used a 2 × 2 × 2 supercell with 5 × 3 × 5 k-point mesh considering AFM ordering (SOC not included) using PHONOPY software 44 .