Journal of Alloys and Compounds 823 (2020) 153793

Contents lists available at ScienceDirect

Journal of Alloys and Compounds

journal homepage: http://www.elsevier.com/locate/jalcom

High pressure X-ray diffraction study of sodium (Na2O): Observations of amorphization and equation of state measurements to 15.9 GPa

Xiaoxin Wu a, b, Yan Zhang a, b, Junkai Zhang a, b, Ran Liu c, Jinghai Yang a, b, Bin Yang d, * Hongxin Xu a, b, Yanz Ma a, b, a Key Laboratory of Functional Materials Physics and Chemistry of the Ministry of Education, Jilin Normal University, Changchun, 130103, China b College of Physics, Jilin Normal University, Siping, 136000, China c State Key Laboratory of Superhard Materials, College of Physics, Jilin University, Changchun, 130012, China d Center for High Pressure Science and Technology Advanced Research (HPSTAR), Changchun, 130012, China article info abstract

Article history: We have investigated the high-pressure behavior of (Na2O) up to 30 GPa by synchrotron Received 31 August 2019 angle-dispersive powder X-ray diffraction in a diamond anvil cell at room temperature. Between 15.9 and Received in revised form 17.3 GPa crystalline Na2O transforms to an amorphous state, and on decompression the amorphous 8 January 2020 structure recrystallizes back between 4.9 and 10.3 GPa. The zero-pressure bulk modulus of Na Ois Accepted 9 January 2020 2 obtained in experimental observations for the first time. The pressure-volume data is fitted according to Available online 11 January 2020 the third-order Birch-Murnaghan equation of state and yields four bulk modulus values by setting the V0 or B ’ either as a constant or variable. The bulk modulus values of Na O are smaller than a-Li O, which is Keywords: 0 2 2 High-pressure in good agreement with the previous calculation results. © Alkali-metal oxide 2020 Elsevier B.V. All rights reserved. Bulk modulus Amorphous state

1. Introduction Alkali-metal have already been studied extensively by a variety of theoretical methods to investigate the lattice constants, Alkali-metal oxides are the materials with potential techno- elastic properties, band structures, superionic behavior, structural logical utility due to their remarkable and interesting physical and optical properties [9e13]. Nevertheless, the study based on properties, which could be used in catalytic reactions [1,2], solid- high pressure experiment is limited and always focused on Li2O. state batteries [3,4], gas storage [5], gas detectors [6], and fuel Kunc et al. identified a high pressure phase transition of Li2O in the cells [7]. The X2O(X¼ Li, Na, K and Rb) compounds crystallize in the range of 45e50 GPa by powder X-ray diffraction [14]. Lazicki et al. cubic anti-fluorite structure (space group Fm3m) [8]. The crystal exhibited the concrete data of the anti-fluorite a-Li2O structure and structure and atomic fractional coordinates are shown in Fig. 1 and the anti-cotunnite b-Li2O structure by high pressure powder X-ray Table 1. In oxides, the cation and anion sublattices have diffraction and Raman scattering [15]. However, there is hardly any different symmetries. The negative oxygen ions form a face- high pressure experimental study on other alkali-metal oxides. centered (FCC) lattice and occupy the corner positions, while the Sodium oxide (Na2O) plays an important role among the alkali- positive alkali-metal ions form a simple cubic lattice and are situ- metal oxides as the mantle’s sixth abundant oxide [16]. Na2Ois ated in oxygen tetrahedral. Compared with fluorite type-structure not stable in nature, and it is difficult to prepare and preserve at in CaF2, the interchange of cation and anion positions lead to the ambient conditions. Despite the role it plays in geosciences and its formation of large cavities in the crystal structure. potential technological utility like other alkali-metal oxides, little experimental work has been reported on Na2O. Zintl et al. deter- mined that the lattice parameter of Na2O crystal is 5.55 Å using powder diffraction [17]. Barrie et al. reported measurements on Na * Corresponding author. Key Laboratory of Functional Materials Physics and Chemistry of the Ministry of Education, Jilin Normal University, Changchun, 130103, metal and on Na2O using XPS, X-ray excited Auger spectroscopy and China. electron excited Auger spectroscopy [18]. More concrete data on E-mail address: [email protected] (Y. Ma). https://doi.org/10.1016/j.jallcom.2020.153793 0925-8388/© 2020 Elsevier B.V. All rights reserved. 2 X. Wu et al. / Journal of Alloys and Compounds 823 (2020) 153793

6e13.7 GPa, and 200s above 15.9 GPa, respectively. All measure- ments were made at room temperature. The software FIT2D was used to integrate the diffraction pattern image into a 2-theta rela- tive intensity. Material Studio software was used to build the crystal structure and calculate the XRD patterns of Na2O and NaOH-I. Selected high-pressure synchrotron XRD patterns were fitted by Rietveld method [27], using the GSAS software [28]. The Rietveld refinement allows crystal structural information to be obtained from powder diffraction data. The peaks of the gasket have been removed before refnements using peakfit software.

3. Results and discussion

Representative high-pressure XRD patterns of the sample at 0.5 GPa are shown in Fig. 2. The diffraction peaks of the Na2O with cubic Fm3m symmetry and the NaOH-I (marked with I in Fig. 2) with orthorhombic Cmcm symmetry are denoted with Miller

Fig. 1. FCC crystal structure of X2O(X¼ Li, Na, K and Rb). indices. A small portion of NaOH is contained in the sample and several peaks are merged into the XRD patterns of gasket, only a part of diffraction peaks of NaOH-I can be observed. The corre- Table 1 sponding lattice constants are a ¼ 5.544 ± 0.002 Å for Na O and ¼ 2 Atomic Wyckoff positions and fractional coordinates of X2O(X Li, Na, K and Rb) at ¼ ± ¼ ± ¼ ± ambient conditions. a 3.304 0.003 Å, b 3.304 0.003 Å, c 11.32 0.02 Å for NaOH-I, in consistence with the previously reported values under Atom Wyckoff x y z ambient conditions [17,29].

X2OO 4a 000 The sample was compressed up to 30.1 GPa at room tempera- X8h 0.25 0.25 0.25 ture. Fig. 3 displays the representative XRD patterns at selected pressures. The diffraction peaks of Na2O (marked with Miller indices) are observable up to 15.9 GPa but all disappear at 17.3 GPa. cohesive energy, equilibrium geometry, lattice constant, bulk No new diffraction peak is observed above 15.9 GPa, except the modulus, and G-point phonon, electronic structure, ground-state remaining diffraction peaks of NaOH-II and the gasket, which in- properties and quasiparticle band structures of Na2Ohavebeen dicates Na2O transforms from a crystalline structure to an amor- obtained from theoretical calculations by various methods including phous state between 15.9 and 17.3 GPa. The transformation e LCAO, HF, DFT, TB-LMTO, LDA, GGA, etc [9,10,19 26]. These theo- pressure (15.9 GPae17.3 GPa) of Na2O is lower than that of Li2O retical calculation results still need the experimental support. (50 ± 5 GPa) [15], and analogous to alkali-metal sulfides, of which In this work, we present the high pressure structural behavior of transformation pressures decrease with increasing cation size (Li2S Na2O by angle-dispersive X-ray diffraction at room temperature up (12.5 GPa) [30], Na2S(7GPa)[31], K2S (2.7 GPa) [32], and Rb2S to 30.1 GPa. An amorphous transformation of Na2O has been (~2 GPa) [32]). However, the above-mentioned alkali-metal sulfides fl observed. The bulk modulus of cubic anti- uorite Na2O has been and Li2O have been shown or predicted to undergo a trans- fi investigated by experiments for the rst time and compared with formation from antifluorite to anticotunnite structure (Cs2S has an previous calculations. anticotunnite structure at ambient conditions). But in our work, no new diffraction peak or appreciable drop in signal can be observed 2. Experimental (the intensity of NaOH diffraction peaks show no sign of abating). A similar trend among two compounds with cations from the same The sample is acquired commercially from Alfa Aesar. It consists main group and the same anions, which one changes into a of 80% Na2O and 20% Na2O2 in the powdered form. The Na2O2 amorphous state but another dose not can also be observed in transformed to NaOH during the sample loading process due to the Ca(OH)2 and Mg(OH)2 [33]. The crystalline Ca(OH)2 transforms into presence of moisture (2Na2O2 þ 2H2O ¼ 4NaOH þ O2[). The a glassy state between 10.7 and 15.4 GPa, but Mg(OH)2 does not vaporized oxygen has been detected in the golvebox used for amorphize up to the highest pressures [34]. Therefore, unlike Li2O, loading; and NaOH has been confirmed by its diffraction patterns as it could be possible for Na2O to transforms into an amorphous state shown in the later section of this paper. High-pressure experiments at the pressure between 15.9 and 17.3 GPa. were performed in a symmetric diamond anvil cell (DAC) with flat The diffraction peaks marked I and II represent NaOH-I and anvil of 500 mm in diameter. A T301 steel sheet was served as gasket NaOH-II, respectively. Beck [29] and Loveday [35] reported that with a drilled hole of 200 mm in diameter and 50 mm in thickness. NaOH exhibits a structural phase transition from orthorhombic The Na2O samples and a ruby ball were placed in the sample (Cmcm, NaOH-I) to orthorhombic (Pbcm, NaOH-II) structure at chamber. The ruby ball was used to determine pressure by using ~0.9 GPa, which is observed in our study as well. The diffraction the standard ruby fluorescent technique. Na2O is extremely sensi- peaks of NaOH-I can be clearly observed at 0.5 GPa but disappear at tive to moisture and readily hydrolyzes with common pressure- 2.1 GPa. Meanwhile, the diffraction peaks of NaOH-II start to transmitting media. In order to restrain hydrolyzation, no pres- emerge at 2.1 GPa and retain up to 30.1 GPa. The peaks of NaOH-II sure transmission medium was added in the sample chamber and seem almost disappear at 30.1 GPa. But the peak at ~24 deg. (the the sample loading was processed in a glove box. In-situ synchro- miller index is (1 2 2)) also can be observed, which indicate the tron X-ray diffraction experiments were performed at BL14B1 sta- NaOH-II structure is stable up to at least 30.1 GPa. This pressure is tion of Shanghai Synchrotron Radiation Facility (SSRF) by using a higher than previous reports of stability of NaOH-II structure up to beam wavelength of 0.6887 Å. Diffraction images were collected 9.7 GPa by XRD spectra [36] and 25 GPa by Raman spectra, using Mar 345 image plate detector. The collecting times of the respectively [37]. The small amount of NaOH does not affect the diffraction images were 100s below 5.2 GPa, 150s between structural phase transition of Na2O between 2.1 GPa and 30.1 GPa X. Wu et al. / Journal of Alloys and Compounds 823 (2020) 153793 3

Fig. 2. The XRD pattern of the sample and pure gasket at 0.5 GPa. The diffraction peaks of Na2O and NaOH-I are marked with Miller indices (I is used to substitute NaOH-I). The diffraction peaks of gasket are marked with *. The broad peaks denoted by ‘*’ are attributed to those of the gasket by comparing with the pattern of blank gasket measured in the experiment. The X-ray beam had inevitably passed through the gasket due to the large beam size (~200 mm). The calculated XRD patterns of Na2O and NaOH-I based on literatures [17,29] at ambient conditions are also given for comparison.

Fig. 3. (a) Representative XRD patterns of the sample at different increasing pressures. The peaks of Na2O, NaOH-I (marked with I) and NaOH-II (marked with I) are denoted by miller indices. The peaks of the gasket are denoted by asterisks (*). (b) Rietveld refnements of the XRD patterns at 0.5 GPa and 7.0 GPa. The peaks of the gasket have been removed before refnements. because NaOH-II is stable in this range. The high-pressure changes phases of NaOH coexist at the pressure. NaOH-II completely returns of the gasket peaks marked with * are depicted by the solid lines. to NaOH-I at 0.5 GPa. It is also a reversible phase transition between These peaks have no effect on the observation of high pressure NaOH-I and NaOH-II. behaviors of Na2O and NaOH, because they could be easily distin- The lattice parameters of FCC-phase Na2O as a function of guished by the wide shapes. pressure are shown in Fig. 5a. The cell parameter a shrink by 5.47% As shown in Fig. 4, during decompression, several diffraction at 13.7 GPa. Fig. 5b shows that the experimental pressure-volume peaks of the FCC-phase Na2O can be observed below 4.9 GPa, which data were fitted by third-order BirchMurnaghan (BM) equation indicates Na2O returns to the ambient-pressure phase in the of state (EOS). The solid points of lattice parameters a and volume V pressure region between 4.9 GPa and 10.3 GPa. There is a hysteresis are the data from literature at ambient conditions [17] and the in this reversible phase transition. The diffraction peaks of NaOH-I hollow points with error bars are the experimental data from our structure start to appear at 0.9 GPa during decompression. Two study. We obtain four sets of B0 data of Na2O. If the V0 is not 4 X. Wu et al. / Journal of Alloys and Compounds 823 (2020) 153793

Fig. 4. Representative XRD patterns of the sample at different pressures during the process of decreasing pressure. The peaks of Na2O, NaOH-I (marked with I) and NaOH-II (marked with II) are denoted by miller indices.

Fig. 5. (a) Pressure dependence of the cell parameter a of FCC structure Na2O. (b) Unit cell volume dependence on pressure. The solid symbols are the data at ambient conditions from literature and the hollow symbols with error bars are the experiment data form our study.

restricted, the fit of the pressure-volume (P-V) data to the third- B0 ¼ 53.4 ± 0.5 GPa and the ambient pressure volume is 3 3 order BM-EOS yields a bulk modulus value of B0 ¼ 44 ± 2 GPa V0 ¼ 171.7 ± 0.2 Å . If the V0 is fixed as 170.95 Å [17], we obtain that with B0’ ¼ 6.0 ± 0.4 and an ambient pressure volume of B0 ¼ 55.6 ± 0.8 GPa with B0’ ¼ 4.0 ± 0.3. When B0’ is fixed as 4, the 3 V0 ¼ 172.7 ± 0.3 Å . When B0’ is fixed as 4, the bulk modulus is bulk modulus is B0 ¼ 55.7 ± 0.2 GPa. In Table 2, the lattice X. Wu et al. / Journal of Alloys and Compounds 823 (2020) 153793 5

Table 2 Acknowledgments Comparison of lattice parameter a0 and bulk modulus B0 between our results and previous calculations. This work was supported by the National Natural Science ’ Method a0 (Å) B0 (GPa) B0 Foundation of China (Grant No. 11847126, 11904128, 21801092, and

This work (Expt.) fixed V0 5.55 [17] 55.6(8) 4.0(3) 21878119) and the Thirteenth Five-Year Program for Science and unfixed B0’ Technology of Education Department of Jilin Province (Item No. fi xed V0 5.55 [17] 55.7(2) 4 JJKH20180769KJ and No. 20180201079GX). In-situ high-pressure fi ’ xed B0 XRD experiments were performed at the BL14B1 station of unfixed V0 5.568(1) 44.1(7) 6.0(3) Shanghai Synchrotron Radiation Facility (SSRF).. unfixed B0’ unfixed V0 5.558(3) 53.4(5) 4 fixed B0’ Dovesi [19] (Cal.) LCAO 5.484 57.5 e References Shukla [20] (Cal.) HF 5.481 61.1 e Cancarevic[9] (Cal.) HF 5.45 62.180 e [1] C.T. Campbell, The selective epoxidation of ethylene catalyzed by Ag(111): a e DFT 5.47 58.631 e comparison with Ag(110), J. Catal. 94 (1985) 436 444. Eithiraj [10] (Cal.) TB-LMTO 5.465 59.00 e [2] H. Tsuji, F. Yagi, H. Hattori, Basic sites on alkali ion-added zeolite, Chem. Lett. 20 (1991) 1881e1884. Thompson [22] (Cal.) LDA 5.398 56 e [3] M. Jamal, G. Venugopal, M. Shareefuddin, M. Narasimha Chary, Sodium ion GGA 5.583 54 e conducting glasses with mixed glass formers NaIeNa2OeV2O5eB2O3: appli- cation to solid state battery, Mater. Lett. 39 (1999) 28e32. [4] K.I. Cho, S.H. Lee, K.H. Cho, D.W. Shin, Y.K. Sun, Li2OeB2O3eP2O5 solid elec- trolyte for thin film batteries, J. Power Sources 163 (2006) 223e228. parameters a0 and the bulk modulus B0 are compared between our [5] S. Wang, S. Yan, X. Ma, J. Gong, Recent advances in capture of experimental results and previous calculated results using different using alkali-metal-based oxides, Energy Environ. Sci. 4 (2011) 3805e3819. [6] K.B. Lee, M.G. Beaver, H.S. Caram, S. Sircar, Performance of Na O promoted methods. The B0 in the literatures is calculated from the second 2 alumina as CO chemisorbent in sorption-enhanced reaction process for derivative of total energy versus volume. It could be seen that our 2 simultaneous production of fuel-cell grade H2 and compressed CO2 from results are similar with the calculations and agree well with the synthesis gas, J. Power Sources 176 (2008) 312e319. results of Thompson et al. [22]. The agreement between experi- [7] H.-K. Lee, J.-P. Shim, M.-J. Shim, S.-W. Kim, J.-S. Lee, The characteristics of mental and calculated EOS of Na O suggests that pressure is synthesized potassium hexatitanate and the manufacturing process of the 2 matrix, Mater. Chem. Phys. 45 (1996) 243e247. reasonably hydrostatic in this phase. The experimental bulk [8] F.A. Shunk, Constitution of Binary Alloys, 2nd Suppl., McGraw-Hill, New York, modulus value of Na2O is smaller than a-Li2O(B0 ¼ 90 ± 1 GPa) [15], 1969. € suggesting that the compressibility of Na O is larger than a-Li O. [9] Z. Cancarevic, J.C. Schon, M. Jansen, Stability of alkali-metal oxides as a 2 2 function of pressure: theoretical calculations, Phys. Rev. B 73 (2006) 224114. The different sizes of the cation atoms could account for this. [10] R.D. Eithiraj, G. Jaiganesh, G. Kalpana, Electronic structure and ground-state properties of alkali-metal oxideseLi2O, Na2O, K2O and Rb2O: a first- e 4. Conclusion principles study, Phys. B Condens. Matter 396 (2007) 124 131. [11] M. Moakafi, R. Khenata, A. Bouhemadou, H. Khachai, B. Amrani, D. Rached, M. Rerat, Electronic and optical properties under pressure effect of alkali e Na2O crystallizes in the cubic-antifluoride structure at ambient metal oxides, Eur. Phys. J. B 64 (2008) 35 42. conditions. This substance is sensitive to moisture. The commer- [12] M.V. Aleinikova, Y.N. Zhuravlev, D.V. Korabelnikov, Ab initio investigation of the effect of pressure on the structure and electronic properties of alkali metal cially available Na2O powder contains 20% Na2O2, and the Na2O2 oxides and peroxides, Russ. Phys. J. 55 (2012) 495e500. transforms into NaOH because of the deliquescence. In this work, [13] B. Baumeier, P. Krüger, J. Pollmann, G.V. Vajenine, Electronic structure of alkali-metal fluorides, oxides, and nitrides: density-functional calculations the mixture of Na2O and NaOH sample was compressed up to including self-interaction corrections, Phys. Rev. B 78 (2008) 125111. 30.1 GPa for the first time at room temperature. Na O shows a 2 [14] K. Kunc, I. Loa, A. Grzechnik, K. Syassen, Li2O at high pressures: structural reversible structural phase transition at between 15.9 GPa and properties, phase-transition, and phonons, Phys. Status Solidi (B) 242 (2005) 17.3 GPa from antifluorite to an amorphous state. A hysteresis is 1857e1863. [15] A. Lazicki, C.S. Yoo, W.J. Evans, W.E. Pickett, Pressure-induced antifluorite-to- found during decompression. The zero-pressure bulk modulus of anticotunnite phase transition in , Phys. Rev. B 73 (2006) Na2O is the first experiment based reports. Fitting the pressure- 184120. volume data to the third-order Birch-Murnaghan equation of [16] Gunter Faure, Principles and Applications of Geochemistry : a Comprehensive fi ’ Textbook for Geology Students, Prentice Hall, Upper Saddle River, N.J., 1998. state we obtain four sets of B0 data of Na2Oby xing the V0 or B0 [17] E. Zintl, A. Harder, B. Dauth, Gitterstruktur der Oxyde, Sulfide, Selenide und either as a constant or variable. NaOH also exhibits a reversible Telluride des Lithiums, Natriums und Kaliums, Zeitschrift für Elektrochemie structural phase transition from orthorhombic (Cmcm, NaOH-I) to und angewandte physikalische Chemie 40 (1934) 588e593. orthorhombic (Pbcm, NaOH-II) structure at the range of [18] A. Barrie, F.J. Street, An Auger and X-ray photoelectron spectroscopic study of sodium metal and sodium oxide, J. Electron. Spectrosc. Relat. Phenom. 7 0.5e2.1 GPa, which is in line with previous study. The NaOH-II (1975) 1e31. structure is stable up to at least 30.1 GPa, which is higher than [19] R. Dovesi, C. Roetti, C. Freyria-Fava, M. Prencipe, V.R. Saunders, On the elastic previous reports that NaOH-II structure is found to be stable up to properties of lithium, sodium and . An ab initio study, Chem. Phys. 156 (1991) 11e19. 9.7 GPa by XRD spectra and 25 GPa by Raman spectra, respectively. [20] A. Shukla, M. Dolg, P. Fulde, H. Stoll, Towards a quantum-chemical description of crystalline insulators: a Wannier-function-based Hartree-Fock study of Li2O e Declaration of competing interest and Na2O, J. Chem. Phys. 108 (1998) 8521 8527. [21] E.A. Mikajlo, K.L. Nixon, M.J. Ford, Electron momentum spectroscopy and linear combination of atomic orbitals calculation of bulk Na2O, J. Phys. Con- The authors declare that they have no known competing dens. Matter 15 (2003) 2155. financial interests or personal relationships that could have [22] M. Thompson, X. Shen, P.B. Allen, Density functional calculation of electronic fl structure and phonon spectra of Na2O, Phys. Rev. B 79 (2009) 113108. appeared to in uence the work reported in this paper. [23] M. Thompson, Density Functional Calculations of Physical Properties of So- dium Oxide (Master Thesis), Stony Brook University of Geosciences, New CRediT authorship contribution statement York, 2011. [24] C. Sommer, P. Krüger, J. Pollmann, Quasiparticle band structure of alkali-metal fluorides, oxides, and nitrides, Phys. Rev. B 85 (2012) 165119. Xiaoxin Wu: Writing - original draft. Yan Zhang: Data curation. [25] M.M. Sinha, A. Sharma, Phonon spectra of superionic sodium and potassium Junkai Zhang: Formal analysis. Ran Liu: Resources. Jinghai Yang: oxides, Solid State Ion. 225 (2012) 211e213. [26] M.C. Masedi, P.E. Ngoepe, H.M. Sithole, Beyond lithium-ion batteries: a Project administration. Bin Yang: Resources. Hongxin Xu: Super- computational study on Na-S and Na-O batteries, IOP Conf. Ser. Mater. Sci. vision. Yanz Ma: Writing - review & editing. Eng. 169 (2017), 012001. 6 X. Wu et al. / Journal of Alloys and Compounds 823 (2020) 153793

[27] H.M. Rietveld, A profile refinement method for nuclear and magnetic struc- phases for alkali metal sulfides, J. Chem. Phys. 121 (2004) 2289e2304. tures, J. Appl. Crystallogr. 2 (1969) 65e71. [33] C. Meade, R. Jeanloz, Static compression of Ca(OH)2 at room temperature: [28] B.H. Toby, EXPGUI, a graphical user interface for GSAS, J. Appl. Crystallogr. 34 observations of amorphization and equation of state measurements to 10.7 (2001) 210e213. GPa, Geophys. Res. Lett. 17 (1990) 1157e1160. [29] H.P. Beck, G. Lederer, High pressure transformations of NaOH, J. Chem. Phys. [34] S.M. Sharma, S.K. Sikka, Pressure induced amorphization of materials, Prog. 98 (1993) 7289e7294. Mater. Sci. 40 (1996) 1e77. [30] A. Grzechnik, A. Vegas, K. Syassen, I. Loa, M. Hanfland, M. Jansen, Reversible [35] J.S. Loveday, W.G. Marshall, R.J. Nelmes, S. Klotz, G. Hamel, J.M. Besson, The antifluorite to anticotunnite phase transition in Li2S at high pressures, J. Solid structure and structural pressure dependence of sodium deuteroxide-V by State Chem. 154 (2000) 603e611. neutron powder diffraction, J. Phys. Condens. Matter 8 (1996) L597. [31] A. Vegas, A. Grzechnik, K. Syassen, I. Loa, M. Hanfland, M. Jansen, Reversible [36] J.W. Otto, W.B. Holzapfel, High-pressure structural studies of MOH layered phase transitions in Na2S under pressure: a comparison with the cation array compounds (M¼Na, K, Rb, Cs), J. Phys. Condens. Matter 7 (1995) 5461. in Na2SO4, Acta Crystallogr. B 57 (2001) 151e156. [37] M.P. Krobok, P.G. Johannsen, W.B. Holzapfel, Raman and FTIR study of NaOH [32] J.C. Schon,€ Z. Cancarevic, M. Jansen, Structure prediction of high-pressure and NaOD under pressure, J. Phys. Condens. Matter 4 (1992) 8141.