CuS2的高压拉曼光谱和X射线衍射

姜峰 赵慧芳 谢亚飞 姜昌国 谭大勇 肖万生

引用本文:
Citation:

CuS2的高压拉曼光谱和X射线衍射

    作者简介: 姜 峰(1995-),男,硕士研究生,主要从事高压矿物学研究. E-mail:fengjiang@gig.ac.cn;
    通讯作者: 肖万生, wsxiao@gig.ac.cn
  • 中图分类号: O521.2;P311.9

High Pressure Raman Spectroscopy and X-ray Diffraction of CuS2

    Corresponding author: XIAO Wansheng, wsxiao@gig.ac.cn
  • CLC number: O521.2;P311.9

  • 摘要: 利用金刚石压腔在高温高压条件下合成了黄铁矿结构CuS2,结合显微拉曼光谱和同步辐射X射线衍射实验,发现黄铁矿结构CuS2在低于30 GPa的压力范围内保持稳定,无结构相变发生。高压拉曼光谱研究结果显示,CuS2的所有拉曼模频率随压力升高呈连续单调线性增加。对X射线衍射实验获得的CuS2体积随压力变化关系进行Birch-Murnaghan状态方程拟合,得到零压晶胞体积V0 = 193.8(5)Å3,体积模量K0 = 99(2)GPa,K0'=4(固定)。运用第一性原理理论计算,得到CuS2的拉曼模频率和晶胞体积随压力变化关系均与实验观测结果保持一致。与其他黄铁矿结构过渡金属二硫化物MS2(M = Mn,Fe,Co,Ni)对比,发现MS6配位八面体大小(M—S键长)主导晶胞的体积大小和晶体的压缩性,并推测CuS2中Cu离子可能以+2价形式存在。研究结果弥补了黄铁矿结构CuS2高压拉曼光谱和X射线衍射研究的缺失,证实其在高温高压下的结构稳定性,对全面认识CuS2的物理化学性质及黄铁矿结构物质的统一规律具有重要价值,对于探索Cu在地球深部的价态和赋存形式也具有指示意义。
  • 图 1  黄铁矿(Pa3)MX2晶体结构示意图(蓝色大球和黄色小球分别表示M和X原子)

    Figure 1.  Pyrite-type(Pa3) MX2 crystal structure (M and X atoms are shown as big blue and small yellow spheres, respectively.)

    图 2  不同条件下4组实验反应前后的特征拉曼光谱(LH和RH分别代表激光加温和电阻丝加温,4组实验中加温前所测拉曼信号均为单质硫的高压相,加温后产物为黄铁矿结构CuS2

    Figure 2.  Raman spectra of the four experiments at different conditions before and after reaction (LH and RH represent laser heating and resistance heating, respectively. Before heating, all the Raman peaks in four experiments belong to the high-pressure phase of elemental sulfur. After heating, the reaction products are pyrite structure CuS2.)

    图 3  CuS2常温加压(a)、卸压(b)过程的代表性拉曼光谱

    Figure 3.  Representative Raman spectra of CuS2 in the process of increasing pressure (a) and releasing pressure (b)

    图 4  实验和理论计算CuS2的拉曼频率随压力变化关系(黑色圆形和菱形分别对应Exp4合成的CuS2在常温加压和卸压过程的结果,三角形表示Exp2合成CuS2卸压过程的数据,星形为理论计算值,实线和虚线分别为实验和计算数据线性拟合结果)

    Figure 4.  Experimental and theory calculated pressure dependence of Raman vibrational modes of CuS2(Black spheres and diamonds correspond to the results of compression process and decompression process of CuS2 synthesised at Exp4,respectively. Black triangles are data of decompression process of CuS2 synthesised at Exp2. The stars represent theory calculated points. Solid line and dotted line are linnear fitted with all experimental and calculated data, respectively.)

    图 5  6.2 GPa压力X射线粉末衍射谱Rietveld结构精修结果

    Figure 5.  Rietveld refinement analyses of the powder XRD pattern at 6.2 GPa

    图 6  在加压(a)和卸压(b)过程中CuS2在不同压力下的X射线衍射图谱(下三角指示产物中的杂质峰,黄铁矿结构CuS2和传压介质KCl的衍射峰位在图中用短棒标出)

    Figure 6.  X-ray powder diffraction patterns at different pressures on compression (a) and decompression (b)(The inverted triangles indicate impurity peaks. The diffraction peaks of pyrite structure CuS2 and pressure transmitting medium KCl are represented with vertical bars.)

    图 7  理论计算和衍射实验的CuS2及其它过渡金属二硫化物MS2(M = Mn,Fe,Co,Ni,Cu)的体积-压力变化关系(黑色圆形和黑色三角形分别对应X射线衍射实验和第一性原理理论计算结果。黑色实线和黑色虚线是对实验和计算数据进行Birch-Murnaghan 状态方程拟合结果。蓝色、绿色、紫色、棕色和红色实线分别为文献给出的黄铁矿结构MnS2、FeS2、CoS2、NiS2和CuS2状态方程结果。)

    Figure 7.  Pressure-volume relationship of experimental and calculated CuS2 and other transition metal disulfides MS2(M = Mn,Fe,Co,Ni,Cu) (The black circles and triangles represent the experimental and observed data of the CuS2, respectively. The solid black line and dotted black line are the Birch-Murnaghan equation of state fit respectively with listed parameters. The blue, green, purple, yellow, and red line shows the equation of state of the pyrite structure MnS2, FeS2, CoS2, NiS2, and CuS2 in references.)

    表 1  黄铁矿结构CuS2的拉曼频率随压力变化及格临爱森参数(γ)

    Table 1.  Pressure dependences of Raman modes and the Grüneisen parameters (γ) of pyrite-type CuS2

    Symmetry classificationω0/cm−1(dω/dp)/(cm−1·GPa−1)γMethod
    This workRef.[38]
    Eg + Tg(1)2132072.451.14Exp.
    Tg(2)2662643.061.14
    Ag + Tg(3)5125121.750.34
    Tg(1)2062.330.97Calc.
    Eg2072.310.96
    Tg(2)2573.111.04
    Tg(3)4951.800.31
    Ag4991.410.24
    Note: The K0(Exp.) = 99 and K0(Calc.) = 85.6 are used respectively to calculate the Grüneisen parameters. They are obtained from the fitted Birch-Murnaghan EOS in this work.
    下载: 导出CSV

    表 2  黄铁矿结构CuS2与其他黄铁矿结构过渡金属二硫化物的零压体积模量和晶胞体积对比

    Table 2.  Comparison of zero-pressure bulk modulus (K0) and unit formula volume (V0) of CuS2 with that of other pyrite structure transition-metal disulfides

    CompositionsPressure range/GPaV03K0/GPaK0Reference
    MnS20−11.7225.74(0)65.9(3)5.1(2)Ref.[14]
    FeS20−80159.00(7)140.2(15)5.52(19)Ref.[41]
    CoS20−20169.68(1)94(2)6.9(5)Ref.[42]
    NiS20−150178.32102.14.6Ref.[43]
    CuS20−7194.10(10)98.8(6)4(fix)Ref.[27]
    CuS20−29.6193.8(5)99(2)4(fix)This study (Exp.)
    CuS20−30196.5(2)85.6(7)4(fix)This study (Calc.)
    Note: All MS2(M = Mn,Fe,Co,Ni,Cu) are pyrite structure (Pa3) in p-range. All the V0 are uniformly transformed to same units for comparison purposes.
    下载: 导出CSV
  • [1] BROSTIGEN G, KJEKSHUS A. Redetermined crystal structure of FeS2 (Pyrite) [J]. Acta Chemica Scandinavica, 1969, 23(6): 2186–2188. doi: 10.3891/acta.chem.scand.23-2186
    [2] BITHER T A, BOUCHARD R, CLOUD W, et al. Transition metal pyrite dichalcogenides. High-pressure synthesis and correlation of properties [J]. Inorganic Chemistry, 1968, 7(11): 2208–2220.
    [3] NOWACK E, SCHWARZENBACH D, HAHN T. Charge densities in CoS2 and NiS2 (pyrite structure) [J]. Acta Crystallographica Section B: Structural Science, 1991, 47(5): 650–659. doi: 10.1107/S0108768191004871
    [4] MAKOVICKY E. Crystal structures of sulfides and other chalcogenides [J]. Reviews in Mineralogy and Geochemistry, 2006, 61(1): 7–125. doi: 10.2138/rmg.2006.61.2
    [5] TEMPLETON D H, DAUBEN C H. The crystal structure of sodium superoxide [J]. Journal of the American Chemical Society, 1950, 72(5): 2251–2254. doi: 10.1021/ja01161a103
    [6] KJEKSHUS A, RAKKE T. Preparation and properties of magnesium, copper, zinc and cadmium dichalcogenides [J]. Acta Chemica Scandinavica A, 1979, 33(8): 617–620. doi: 10.3891/acta.chem.scand.33a-0617
    [7] KUWAYAMA Y, HIROSE K, SATA N, et al. The pyrite-type high-pressure form of silica [J]. Science, 2005, 309(5736): 923–925. doi: 10.1126/science.1114879
    [8] SHIRAKO Y, WANG X, TSUJIMOTO Y, et al. Synthesis, crystal structure, and electronic properties of high-pressure PdF2-type oxides MO2(M= Ru, Rh, Os, Ir, Pt) [J]. Inorganic Chemistry, 2014, 53(21): 11616–11625. doi: 10.1021/ic501770g
    [9] YU R, ZHAN Q, DE JONGHE L C. Crystal structures of and displacive transitions in OsN2, IrN2, RuN2, and RhN2 [J]. Angewandte Chemie International Edition, 2007, 46(7): 1136–1140. doi: 10.1002/anie.200604151
    [10] HU Q, KIM D Y, YANG W, et al. FeO2 and FeOOH under deep lower-mantle conditions and Earth’s oxygen-hydrogen cycles [J]. Nature, 2016, 534(7606): 241–244. doi: 10.1038/nature18018
    [11] LIU J, HU Q Y, BI W L, et al. Altered chemistry of oxygen and iron under deep Earth conditions [J]. Nature Communications, 2019, 10(1): 153. doi: 10.1038/s41467-018-08071-3
    [12] KLEPPE A K, JEPHCOAT A P. High-pressure Raman spectroscopic studies of FeS2 pyrite [J]. Mineralogical Magazine, 2004, 68(3): 433–441. doi: 10.1180/0026461046830196
    [13] HARRAN I, CHEN Y Z, WANG H Y, et al. High-pressure induced phase transition of FeS2: electronic, mechanical and thermoelectric properties [J]. Journal of Alloys and Compounds, 2017, 710: 267–273. doi: 10.1016/j.jallcom.2017.03.256
    [14] KIMBER S A J, SALAMAT A, EVANS S R, et al. Giant pressure-induced volume collapse in the pyrite mineral MnS2 [J]. Proceedings of the National Academy of Sciences of the United States of America, 2014, 111(14): 5106–5110. doi: 10.1073/pnas.1318543111
    [15] FUJII T, TANAKA K, MARUMO F, et al. Structural behaviour of NiS2 up to 54 kbar [J]. Mineralogical Journal, 1987, 13(7): 448–454. doi: 10.2465/minerj.13.448
    [16] ELGHAZALI M A, NAUMOV P G, MU Q, et al. Pressure-induced metallization, transition to the pyrite-type structure, and superconductivity in palladium disulfide PdS2 [J]. Physical Review B, 2019, 100(1): 014507. doi: 10.1103/PhysRevB.100.014507
    [17] HUANG S X, WU X, QIN S. Ultrahigh-pressure phase transitions in FeS2 and FeO2: implications for super-earths' deep interior [J]. Journal of Geophysical Research, 2018, 123(1): 277–284. doi: 10.1002/2017JB014766
    [18] BITHER T A, PREWITT C T, GILLSON J L, et al. New transition metal dichalcogenides formed at high pressure [J]. Solid State Communications, 1966, 4(10): 533–535. doi: 10.1016/0038-1098(66)90419-4
    [19] MUNSON R A. The synthesis of copper disulfide [J]. Inorganic Chemistry, 1966, 5(7): 1296–1297. doi: 10.1021/ic50041a055
    [20] BAYLISS P. Crystal chemistry and crystallography of some minerals within the pyrite group [J]. American Mineralogist, 1989, 74(9/10): 1168–1176.
    [21] UEDA H, NOHARA M, KITAZAWA K, et al. Copper pyrites CuS2 and CuSe2 as anion conductors [J]. Physical Review B, 2002, 65(15): 155104. doi: 10.1103/PhysRevB.65.155104
    [22] KAKIHANA M, MATSUDA T D, HIGASHINAKA R, et al. Superconducting and fermi surface properties of pyrite-type compounds CuS2 and CuSe2 [J]. Journal of the Physical Society of Japan, 2019, 88(1): 014702. doi: 10.7566/JPSJ.88.014702
    [23] KING H E, PREWITT C T. Structure and symmetry of CuS2 (pyrite structure) [J]. American Mineralogist, 1979, 64(11/12): 1265–1271.
    [24] MOSSELMANS J F W, PATTRICK R A D, VAN DER LAAN G, et al. X-ray absorption near-edge spectra of transition metal disulfides FeS2 (pyrite and marcasite), CoS2, NiS2 and CuS2, and their isomorphs FeAsS and CoAsS [J]. Physics and Chemistry of Minerals, 1995, 22(5): 311–317. doi: 10.1007/BF00202771
    [25] FOLMER J C W, JELLINEK F, CALIS G H M. The electronic structure of pyrites, particularly CuS2 and Fe1- xCuxSe2: an XPS and Mössbauer study [J]. Journal of Solid State Chemistry, 1988, 72(1): 137–144. doi: 10.1016/0022-4596(88)90017-5
    [26] TOSSELL J A, VAUGHAN D J, BURDETT J K. Pyrite, marcasite, and arsenopyrite type minerals: crystal chemical and structural principles [J]. Physics and Chemistry of Minerals, 1981, 7(4): 177–184. doi: 10.1007/BF00307263
    [27] HÜPEN H, WILL G, HÖFFNER C, et al. X-ray diffraction of CuS2 under high pressure [J]. Materials Science Forum, 1991, 79/80/81/82: 697–702. doi: 10.4028/www.scientific.net/MSF.79-82.697
    [28] NIWA K, TERABE T, SUZUKI K, et al. High-pressure stability and ambient metastability of marcasite-type rhodium pernitride [J]. Journal of Applied Physics, 2016, 119(6): 065901. doi: 10.1063/1.4941436
    [29] TSE J S, KLUG D D, UEHARA K, et al. Elastic properties of potential superhard phases of RuO2 [J]. Physical Review B, 2000, 61(15): 10029–10034. doi: 10.1103/PhysRevB.61.10029
    [30] MAO H K, XU J, BELL P M. Calibration of the ruby pressure gauge to 800 kbar under quasi-hydrostatic conditions [J]. Journal of Geophysical Research, 1986, 91(B5): 4673–4676. doi: 10.1029/JB091iB05p04673
    [31] DATCHI F, LETOULLEC R, LOUBEYRE P. Improved calibration of the SrB4O7: Sm2+ optical pressure gauge: advantages at very high pressures and high temperatures [J]. Journal of Applied Physics, 1997, 81(8): 3333–3339. doi: 10.1063/1.365025
    [32] HAMMERSLEY A P, SVENSSON S O, HANFLAND M, et al. Two-dimensional detector software: from real detector to idealised image or two-theta scan [J]. High Pressure Research, 1996, 14(4/5/6): 235–248. doi: 10.1080/08957959608201408
    [33] GONZALEZ-PLATAS J, ALVARO M, NESTOLA F, et al. EosFit7-GUI: a new graphical user interface for equation of state calculations, analyses and teaching [J]. Journal of Applied Crystallography, 2016, 49(4): 1377–1382. doi: 10.1107/S1600576716008050
    [34] TOBY B H, VON DREELE R B. GSAS-II: the genesis of a modern open-source all purpose crystallography software package [J]. Journal of Applied Crystallography, 2013, 46(2): 544–549. doi: 10.1107/S0021889813003531
    [35] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letters, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865
    [36] ECKERT B, SCHUMACHER R, JODL H J, et al. Pressure and photo-induced phase transitions in sulphur investigated by Raman spectroscopy [J]. High Pressure Research, 2000, 17(2): 113–146. doi: 10.1080/08957950008200934
    [37] PEIRIS S M, SWEENEY J S, CAMPBELL A J, et al. Pressure-induced amorphization of covellite, CuS [J]. The Journal of Chemical Physics, 1996, 104(1): 11–16. doi: 10.1063/1.470870
    [38] ANASTASSAKIS E, PERRY C H. Light scattering and ir measurements in XS2 pryite-type compounds [J]. The Journal of Chemical Physics, 1976, 64(9): 3604–3609. doi: 10.1063/1.432711
    [39] VOGT H, CHATTOPADHYAY T, STOLZ H J. Complete first-order Raman spectra of the pyrite structure compounds FeS2, MnS2 and SiP2 [J]. Journal of Physics and Chemistry of Solids, 1983, 44(9): 869–873. doi: 10.1016/0022-3697(83)90124-5
    [40] SOURISSEAU C, CAVAGNAT R, FOUASSIER M. The vibrational properties and valence force fields of FeS2, RuS2 pyrites and FeS2 marcasite [J]. Journal of Physics and Chemistry of Solids, 1991, 52(3): 537–544. doi: 10.1016/0022-3697(91)90188-6
    [41] THOMPSON E C, CHIDESTER B A, FISCHER R A, et al. Equation of state of pyrite to 80 GPa and 2 400 K [J]. American Mineralogist, 2016, 101(5): 1046–1051. doi: 10.2138/am-2016-5527
    [42] BRAZHKIN V V, DZHAVADOV L N, EL'KIN F S. Study of the compressibility of FeSi, MnSi, and CoS2 transition-metal compounds at high pressures [J]. JETP Letters, 2016, 104(2): 99–104. doi: 10.1134/S0021364016140083
    [43] YUY G, ROSS N L. Prediction of high-pressure polymorphism in NiS2 at megabar pressures [J]. Journal of Physics: Condensed Matter, 2010, 22(23): 235401. doi: 10.1088/0953-8984/22/23/235401
  • [1] 李东飞张可为里佐威刘承志郭瑞孙成林李海波 . 高压下Td-WTe2单晶体材料的拉曼光谱研究. 高压物理学报, 2016, 30(5): 369-374. doi: 10.11858/gywlxb.2016.05.004
    [2] 田雨刘雪廷何运鸿赵慧芳姜峰谭大勇肖万生 . NaCl-O2体系高温高压化学反应的拉曼光谱证据. 高压物理学报, 2017, 31(6): 692-697. doi: 10.11858/gywlxb.2017.06.003
    [3] 高玲玲马艳梅刘丹郝健金云霞王峰王秋实邹广田崔启良 . 环庚烷的高压拉曼光谱研究. 高压物理学报, 2008, 22(2): 192-196 . doi: 10.11858/gywlxb.2008.02.013
    [4] 韩茜吴也黄海军 . BiFeO3高压拉曼光谱研究. 高压物理学报, 2018, 32(5): 051202-1-051202-5. doi: 10.11858/gywlxb.20170698
    [5] 袁真张少鹏靳常青王晓慧 . PbZr0.52Ti0.48O3陶瓷的高压拉曼光谱研究. 高压物理学报, 2015, 29(2): 95-98. doi: 10.11858/gywlxb.2015.02.002
    [6] 何运鸿田雨赵慧芳姜峰谭大勇肖万生 . 高氯酸钠高压相变的拉曼光谱证据. 高压物理学报, 2018, 32(4): 041201-1-041201-9. doi: 10.11858/gywlxb.20180543
    [7] 刘振先崔启良赵永年邹广田 . 传压介质对晶格振动和相变压力的影响-Bi2O3的高压拉曼光谱研究. 高压物理学报, 1990, 4(2): 81-86 . doi: 10.11858/gywlxb.1990.02.001
    [8] 赵慧芳谭大勇姜峰谢亚飞姜昌国罗兴丽肖万生 . Re-H2O体系高温高压化学反应的拉曼证据. 高压物理学报, 2020, 34(4): 040102-1-040102-10. doi: 10.11858/gywlxb.20200518
    [9] 刘莉君陈良臣陈红王汝菊张芝婷车荣钲周镭王积方李玉良姚幼新 . C60高压X射线衍射研究. 高压物理学报, 1993, 7(3): 202-207 . doi: 10.11858/gywlxb.1993.03.006
    [10] 赵金郑海飞 . 0.1~800 MPa压力下方解石拉曼光谱的实验研究. 高压物理学报, 2003, 17(3): 226-229 . doi: 10.11858/gywlxb.2003.03.012
    [11] 秦善刘景李海建祝向平李晓东 . 绿柱石的原位高压X射线研究. 高压物理学报, 2008, 22(1): 1-5 . doi: 10.11858/gywlxb.2008.01.001
    [12] 秦振兴张建波 . 四甲基硅烷的高压拉曼散射研究. 高压物理学报, 2016, 30(5): 375-379. doi: 10.11858/gywlxb.2016.05.005
    [13] 瞿清明郑海飞 . 0.1~3 000 MPa下碳化硅顶砧拉曼光谱作为压力计的研究. 高压物理学报, 2007, 21(3): 332-336 . doi: 10.11858/gywlxb.2007.03.020
    [14] 王一峰苏文辉千正男马贤锋阎学伟 . 稀土焦绿石化合物R2Fe4/3W2/3O7的高温高压稳定性研究和X射线相分析. 高压物理学报, 1988, 2(4): 296-304 . doi: 10.11858/gywlxb.1988.04.002
    [15] 宋海鹏刘云贵李想靳树宇王欣宇巫翔 . 羟碳铈矿的高压拉曼光谱研究. 高压物理学报, 2019, 33(6): 060105-1-060105-8. doi: 10.11858/gywlxb.20190847
    [16] 范大伟许金贵魏舒怡陈志强谢鸿森 . 天然绿帘石原位高压同步辐射X射线衍射研究. 高压物理学报, 2014, 28(3): 257-261. doi: 10.11858/gywlxb.2014.03.001
    [17] 乔二伟郑海飞孙樯 . 甲醇作为压标的拉曼研究. 高压物理学报, 2004, 18(4): 368-373 . doi: 10.11858/gywlxb.2004.04.014
    [18] 张红肖万生谭大勇罗崇举李延春刘景 . 斜锆石(ZrO2)高温高压相变的Raman光谱研究. 高压物理学报, 2007, 21(3): 264-268 . doi: 10.11858/gywlxb.2007.03.008
    [19] 顾惠成陈良辰鲍忠兴李凤英 . Pb0.8Sn0.2Te晶体高压相变的X射线衍射研究. 高压物理学报, 1996, 10(3): 227-230 . doi: 10.11858/gywlxb.1996.03.010
    [20] 刘景赵菁车荣钲杨洋 . 高压下的同步辐射能量色散粉末衍射. 高压物理学报, 2000, 14(4): 247-252 . doi: 10.11858/gywlxb.2000.04.002
  • 加载中
图(7)表(2)
计量
  • 文章访问数:  285
  • 阅读全文浏览量:  236
  • PDF下载量:  5
出版历程
  • 收稿日期:  2020-02-17
  • 录用日期:  2020-02-26
  • 网络出版日期:  2020-04-25
  • 刊出日期:  2020-08-05

CuS2的高压拉曼光谱和X射线衍射

    作者简介:姜 峰(1995-),男,硕士研究生,主要从事高压矿物学研究. E-mail:fengjiang@gig.ac.cn
    通讯作者: 肖万生, wsxiao@gig.ac.cn
  • 1. 中国科学院广州地球化学研究所矿物学与成矿学重点实验室,广东 广州 510640
  • 2. 广东省矿物物理与材料研究开发重点实验室,广东 广州 510640
  • 3. 中国科学院大学,北京 100049

摘要: 利用金刚石压腔在高温高压条件下合成了黄铁矿结构CuS2,结合显微拉曼光谱和同步辐射X射线衍射实验,发现黄铁矿结构CuS2在低于30 GPa的压力范围内保持稳定,无结构相变发生。高压拉曼光谱研究结果显示,CuS2的所有拉曼模频率随压力升高呈连续单调线性增加。对X射线衍射实验获得的CuS2体积随压力变化关系进行Birch-Murnaghan状态方程拟合,得到零压晶胞体积V0 = 193.8(5)Å3,体积模量K0 = 99(2)GPa,K0'=4(固定)。运用第一性原理理论计算,得到CuS2的拉曼模频率和晶胞体积随压力变化关系均与实验观测结果保持一致。与其他黄铁矿结构过渡金属二硫化物MS2(M = Mn,Fe,Co,Ni)对比,发现MS6配位八面体大小(M—S键长)主导晶胞的体积大小和晶体的压缩性,并推测CuS2中Cu离子可能以+2价形式存在。研究结果弥补了黄铁矿结构CuS2高压拉曼光谱和X射线衍射研究的缺失,证实其在高温高压下的结构稳定性,对全面认识CuS2的物理化学性质及黄铁矿结构物质的统一规律具有重要价值,对于探索Cu在地球深部的价态和赋存形式也具有指示意义。

English Abstract

  • 黄铁矿结构(Pa3,Z = 4)是MX2型化合物的一种典型结构,原型为天然矿物黄铁矿(FeS2)。在该立方结构中,M原子占据4a(0, 0, 0)位置,X原子占据8c(u, u, u)位置,形成共角顶的AB6配位八面体[1],如图1所示。一些二硫化物MS2(M = Mn,Fe,Co,Ni,Cu,Zn,Ru,Os等)[2-4]和二氧化物MO2(M = Be,Mg,Zn,Cd,Na等)[5-6]均能够结晶成黄铁矿结构。另外,黄铁矿结构也是许多二氧化物和二氮化物在高压下的稳定多形,如SiO2、GeO2、SnO2、PbO2、FeO2、RuO2、OsO2、IrO2、PtO2、PtN2、PdN2[7-9]。由于Fe和SiO2是地球内部的重要组成元素和化合物,深入研究高压下黄铁矿型结构在地球科学和材料科学的发展过程中具有重要意义。例如,近期黄铁矿型结构FeO2和FeO2Hx的高温高压实验和理论研究进展[10-11]为理解下地幔环境中各种元素的化学性质和原子间的相互作用提供了新视角。

    图  1  黄铁矿(Pa3)MX2晶体结构示意图(蓝色大球和黄色小球分别表示M和X原子)

    Figure 1.  Pyrite-type(Pa3) MX2 crystal structure (M and X atoms are shown as big blue and small yellow spheres, respectively.)

    作为黄铁矿结构的典型代表,过渡金属二硫化物MS2(M = Mn,Fe,Co,Ni,Cu等)由于具有良好的光电性质和催化性能被广泛地研究和应用,其高压行为(如高压稳定性、结构相变和金属-绝缘体转变等)也受到高度关注[12-16]。同时,过渡金属二硫化物还可以成为地球深部和气态巨行星内部SiO2等重要矿物的类比研究对象[17],有利于发现SiO2等黄铁矿结构的高压多形,并了解这些高压多形的物理化学性质。

    与其他黄铁矿结构二硫化物不同,CuS2需要在高温高压条件下才能合成[18-19],导致相关研究相对较少。以CuS2为主要组成成分的天然矿物黑硫铜镍矿(Cu,Ni,Fe)S2需要Ni、Fe等类质同像替代进入晶格才能保持稳定[20]。CuS2是一种深紫色晶体,具有金属导电性[21],并在低温下显示超导电性[22],King等[23]对高温高压合成的CuS2单晶进行了详细的晶体结构研究,并与MS2(M = Mn,Fe,Co,Ni)等晶体结构参数和MS6八面体畸变程度进行比较,认为CuS2中Cu的价态为+2价。然而,CuS2的多种光谱和理论研究认为Cu为+1价态[21, 24-26]。晶体化学与光谱学测量结果之间的矛盾或许隐含有价值的科学信息。

    据相关研究报道,迄今只有Hupen等[27]对CuS2开展过高压原位研究工作。Hupen等利用同步辐射能量色散X射线衍射技术在0~6 GPa压力范围对CuS2开展原位测量,认为黄铁矿型CuS2在约1 GPa压力下会发生向白铁矿结构(Pnnm)的可逆相变。白铁矿结构与文献中常见的CaCl2结构等型[4, 28],而对于大多数MX2化合物的高压结构演化来说,CaCl2结构是黄铁矿结构的低压相[7, 29]。因此,CuS2相变还需要更高精度的高压X射线衍射及其他高压测量方法证实。

    本研究利用金刚石压腔高压装置开展了CuS2的高温高压合成及其高压拉曼光谱和X射线衍射测试。黄铁矿结构CuS2可在6~33 GPa压力范围,通过电阻外加热方式(800 K)和激光加热方式(1 500~2 000 K)合成,并对其进行高压拉曼光谱和X射线衍射测试,探索黄铁矿结构CuS2在0~30 GPa压力范围的结构稳定性。同时,对黄铁矿结构CuS2开展第一性原理理论计算,并将计算结果与实验结果进行对比验证。研究结果对于认识 Cu 在地球深部环境中的价态和赋存形式具有重要意义。

    • 将分析纯化学试剂Cu、S以及CuS、S分别按化学计量比1∶2和1∶1称量,并在玛瑙研钵中研磨1.5 h充分混合,作为合成CuS2的初始反应物。高温高压合成实验采用两种类型的金刚石压腔(DAC)装置:一种为可与显微激光加热系统连接的对称型DAC装置(LH-DAC),光谱仪接收样品热辐射谱,并用黑体辐射公式拟合实验温度;另一种为电阻丝加热的DAC装置(RH-DAC),用NiCr-NiAl热电偶测量实验温度。两种DAC装置的金刚石对顶砧砧面直径均为300 μm。实验采用T301不锈钢片封垫,预压厚度约35 μm,用电火花打孔器在中心钻一个直径105 μm的小孔作为样品腔。共进行了4组CuS2的合成实验。激光加热实验(Exp1和Exp2)初始样品均为Cu+2S混合物,隔热层和压力介质为KCl,采用“三明治”方式装样,并加入红宝石微粒用于测压[30],实验压力分别为7.0和34.1 GPa,实验温度1 500~2 000 K。电阻丝加热实验的装样方法与激光加热实验相同,实验样品分别为Cu+2S(Exp3)和CuS+S(Exp4),采用SrB4O7∶Sm2+作压标[31],实验压力分别为12.4和16.6 GPa,实验温度为800 K,恒温约1 h。通过拉曼光谱鉴定4组实验合成产物。

      拉曼光谱测试采用Renishaw 2000 型显微共聚焦拉曼光谱仪。选用波长为532 nm的激光作为激发光源,用20 × 物镜聚焦激光光斑,样品处光斑大小为3 μm × 3 μm。采用背反射方式收集样品信号,通过1 800 线每毫米光栅分光,由热电冷却的CCD探测器采集信号。采谱波段范围100~800 cm−1,光谱分辨率1 cm−1,采谱时间40 s。高压原位角色散X射线衍射实验在中国科学院高能物理研究所同步辐射装置4W2线站开展。入射X射线波长λ为0.619 9 Å,光束被聚焦到约60 μm × 20 μm大小。Pilatus探测器收集衍射图谱,样品到探测器的距离由CeO2标定。利用Fit2D软件[32]对采集的二维衍射图进行积分处理。Birch-Murnagham状态方程参数由EosFit7程序[33]拟合得到。运用GSAS-II软件[34]对6.2 GPa压力点的粉末衍射谱进行Rietveld结构精修。

      CuS2的第一性原理理论计算采用Material Studio软件的CASTEP计算模块。运用Perdew-Burke-Ernzerhof泛函的广义梯度近似(GGA)方法估算交换相关能[35]。赝势选择OTFG超软赝势,截断能设置为700 eV,第一布里渊区k点取样为6 × 6 × 6网格。为确保较好的收敛,能量、最大力场、最大应力和最大位移分别为5.0 × 10−6 eV/atom、0.01 eV/Å、0.02 GPa和5.0 × 10−4 Å。

    • 图2为4组合成实验高温处理前后的拉曼光谱。加热前,由于CuS高压非晶化,在约15 GPa实验压力下基本难以观测到拉曼峰,而Cu单质无拉曼活性峰,因此4组实验在反应前测得的拉曼光谱均属于单质S的高压相[36-37]。加热后,4组实验都得到由3个拉曼峰组成的全新谱图,表明实验样品在不同压力经不同温度处理后得到了相同的新相。新相在常温常压下的3个拉曼峰频率(ω0)在表1列出,与Anastassakis等[38]给出的黄铁矿结构CuS2的拉曼峰位和峰形特征基本一致,可初步确认实验合成的新相为黄铁矿结构CuS2,后述同步辐射X射线衍射实验结果将进一步证实。已有的CuS2合成实验条件基本均在7.0 GPa、1 500 K温压范围以内[2, 18-19, 21],本次4组实验在更宽范围的温压条件下均合成出黄铁矿结构CuS2,拓展了CuS2的合成温压范围,同时说明了CuS2在34.1 GPa、2 000 K温压范围内是热力学稳定态。

      Symmetry classificationω0/cm−1(dω/dp)/(cm−1·GPa−1)γMethod
      This workRef.[38]
      Eg + Tg(1)2132072.451.14Exp.
      Tg(2)2662643.061.14
      Ag + Tg(3)5125121.750.34
      Tg(1)2062.330.97Calc.
      Eg2072.310.96
      Tg(2)2573.111.04
      Tg(3)4951.800.31
      Ag4991.410.24
      Note: The K0(Exp.) = 99 and K0(Calc.) = 85.6 are used respectively to calculate the Grüneisen parameters. They are obtained from the fitted Birch-Murnaghan EOS in this work.

      表 1  黄铁矿结构CuS2的拉曼频率随压力变化及格临爱森参数(γ)

      Table 1.  Pressure dependences of Raman modes and the Grüneisen parameters (γ) of pyrite-type CuS2

      图  2  不同条件下4组实验反应前后的特征拉曼光谱(LH和RH分别代表激光加温和电阻丝加温,4组实验中加温前所测拉曼信号均为单质硫的高压相,加温后产物为黄铁矿结构CuS2

      Figure 2.  Raman spectra of the four experiments at different conditions before and after reaction (LH and RH represent laser heating and resistance heating, respectively. Before heating, all the Raman peaks in four experiments belong to the high-pressure phase of elemental sulfur. After heating, the reaction products are pyrite structure CuS2.)

      群论分析结果表明,黄铁矿结构CuS2振动模式的不可约表示为Γ = Ag + Eg + 3Tg + 2Au + 2Eu + 6Tu。其中Ag + Eg + 3Tg为5个拉曼活性模,Ag为完全对称模,Eg为二重简并模,Tg为三重简并模。AgTg(2)分别对应S—S原子对同向和反向伸缩振动,Eg摆动模反映S原子在垂直于S—S键轴方向的位移,Tg(1)和Tg(3)则对应伸缩模和摆动模的多种组合[39-40]。利用第一性原理计算黄铁矿结构CuS2的声子谱,得到5个拉曼峰,它们的常压振动频率也在表1中列出。可以清楚地发现,Tg(1)与Eg以及Tg(3)与Ag之间的波数差在5 cm−1以内,使得在非偏振拉曼光谱中难以分辨[39-40],导致本次实验和Anastassakis等[38]的拉曼光谱实验中仅观测到Tg(1)+ EgTg(2)和Tg(3)+ Ag共3组拉曼峰。此外,本研究的实验值和计算值以及Anastassakis等[38]的实验结果三者高度吻合,准确地描述了黄铁矿结构CuS2的常压拉曼振动特征。

    • 黄铁矿结构CuS2由常压加压到33.0 GPa,并随后卸压至常压过程的拉曼光谱变化特征如图3所示。图3(a)中,随着压力增加,3个拉曼峰均向高频方向连续移动,至33.0 GPa没有拉曼峰增加或消失;卸压过程(图3(b))则保持与加压过程可逆的变化特征。该现象表明,高压仅使CuS2晶体结构致密化,S—S键和Cu—S键压缩引起的振动频率增大,没有导致结构相变,反映黄铁矿结构CuS2在0~33.0 GPa压力范围内具有结构稳定性。

      图  3  CuS2常温加压(a)、卸压(b)过程的代表性拉曼光谱

      Figure 3.  Representative Raman spectra of CuS2 in the process of increasing pressure (a) and releasing pressure (b)

      图4展示了拉曼光谱实验测得的3组CuS2拉曼模Tg(1)+ EgTg(2)、Tg(3)+ Ag和第一性原理理论计算得到的5个拉曼模EgTg(1)、Tg(2)、Tg(3)、Ag的频率随压力的变化关系。可以看出,拉曼实验值与理论计算值在对应压力下非常接近,两者的最大频率相差小于4%,反映理论计算能很好地描述实验结果。拉曼实验和理论计算数据显示,所有拉曼模的振动频率随压力升高呈线性增加,对拉曼模分别进行线性拟合,拟合给出的频率随压力变化(dω/dp)在表1中列出,实验和拟合计算结果一致。理论计算数据的拟合结果显示,EgTg(1)模的频率变化速率(dω/dp)非常接近,且常压频率相近,对应拉曼光谱中始终难以分离的情况。而AgTg(3)的频率变化速率略有差异,但Ag模的强度远大于Tg(3)模,在拉曼光谱中基本表现为Ag模信号。

      图  4  实验和理论计算CuS2的拉曼频率随压力变化关系(黑色圆形和菱形分别对应Exp4合成的CuS2在常温加压和卸压过程的结果,三角形表示Exp2合成CuS2卸压过程的数据,星形为理论计算值,实线和虚线分别为实验和计算数据线性拟合结果)

      Figure 4.  Experimental and theory calculated pressure dependence of Raman vibrational modes of CuS2(Black spheres and diamonds correspond to the results of compression process and decompression process of CuS2 synthesised at Exp4,respectively. Black triangles are data of decompression process of CuS2 synthesised at Exp2. The stars represent theory calculated points. Solid line and dotted line are linnear fitted with all experimental and calculated data, respectively.)

    • 图5所示为6.2 GPa压力下(Exp1温度处理后)产物的同步辐射X射线衍射谱,黑色圆形、红色实线和黑色实线分别对应观测值、计算值和误差曲线,在图中以相同比例显示,衍射谱下方黑色垂直短棒为计算CuS2Pa3)和KCl(Pm3m)的衍射峰位。除去传压介质KCl(B2)相的衍射峰外,剩余的衍射峰可以指标化为Pa3结构。随后选择KCl(Pm3m)和CuS2Pa3)作为初始结构模型,利用GSAS-II软件对衍射谱进行结构精修。精修结果如图5所示,可见实验和计算结果的X射线衍射谱吻合较好,误差因子Rp = 1.65%、Rwp = 2.85%。因此,确定本次高温高压实验中合成了黄铁矿结构CuS2Pa3),与拉曼光谱测试结果一致。通过结构精修,得到黄铁矿结构CuS2在6.2 GPa压力下的晶胞参数a = 5.657 3(6)Å,其中原子占位情况为Cu(0,0,0)、S(0.397 4,0.397 4,0.397 4)。

      图  5  6.2 GPa压力X射线粉末衍射谱Rietveld结构精修结果

      Figure 5.  Rietveld refinement analyses of the powder XRD pattern at 6.2 GPa

      图6展示了CuS2由6.2 GPa加压至29.6 GPa,随后卸至常压过程中的代表性X射线衍射谱。随着压力增加,传压介质KCl和产物CuS2的衍射峰均逐渐宽化,反映了非静水压力效应及其导致的样品结晶度降低。然而,直到达到最高压力29.6 GPa,仍可以辨识出黄铁矿结构CuS2的(111)、(200)、(210)、(211)、(220)、(222)、(320)、(321)等衍射峰,表明CuS2依然保持黄铁矿结构。卸压至常压过程中,CuS2的衍射峰仍保持宽状态直到常压,没有再次显现出类似6.2 GPa激光加热后的尖锐峰形,所以压力对样品结晶度变差的影响是不可恢复的。此外,16.7 GPa时开始出现两个杂峰(图6中下三角指示),但它们与其他衍射峰之间无明确几何关系,不能指认为同一种结构。因此,两个杂峰可能对应于激光高温下可能产生的其他铜硫化合物(如CuS等)杂质。黄铁矿结构CuS2在低于29.6 GPa压力范围内X射线衍射谱的变化特征说明其在低于30.0 GPa压力范围内具有结构稳定性,不发生结构相变,与上述拉曼实验结果保持一致,与Hupen等[27]发现的结构相变结论不同。

      图  6  在加压(a)和卸压(b)过程中CuS2在不同压力下的X射线衍射图谱(下三角指示产物中的杂质峰,黄铁矿结构CuS2和传压介质KCl的衍射峰位在图中用短棒标出)

      Figure 6.  X-ray powder diffraction patterns at different pressures on compression (a) and decompression (b)(The inverted triangles indicate impurity peaks. The diffraction peaks of pyrite structure CuS2 and pressure transmitting medium KCl are represented with vertical bars.)

    • 同步辐射X射线衍射实验和第一性原理理论计算得到的黄铁矿结构CuS2的晶胞体积随压力变化关系如图7所示。实验和计算得到的晶胞体积在对应压力下非常接近,对应压力最大体积差小于2%。运用Birch-Murnaghan 状态方程分别对衍射实验数据和理论计算数据进行拟合。将压力导数K0′均固定为4,拟合得到黄铁矿结构CuS2零压晶胞体积V0和体积模量K0分别为:V0(Exp.)= 193.8(5) Å3K0(Exp.)= 99(2) GPa,V0(Calc.)= 196.5(2) Å3K0(Calc.)= 85.6(7) GPa。而Hupen等[27]通过CuS2高压衍射实验给出黄铁矿结构CuS2的状态方程参数为V0(Hupen)= 194.1(2) Å3K0(Hupen)= 98.8(6) GPa且K0' = 4(固定)。所有状态方程参数在表2中列出。可以看出,本次实验所得零压体积和体积模量值与Hupen等[27]给出的结果基本一致,而第一性原理理论计算值与两个衍射实验值之间有少许差异。此外,依据格临爱森参数计算公式

      CompositionsPressure range/GPaV03K0/GPaK0Reference
      MnS20−11.7225.74(0)65.9(3)5.1(2)Ref.[14]
      FeS20−80159.00(7)140.2(15)5.52(19)Ref.[41]
      CoS20−20169.68(1)94(2)6.9(5)Ref.[42]
      NiS20−150178.32102.14.6Ref.[43]
      CuS20−7194.10(10)98.8(6)4(fix)Ref.[27]
      CuS20−29.6193.8(5)99(2)4(fix)This study (Exp.)
      CuS20−30196.5(2)85.6(7)4(fix)This study (Calc.)
      Note: All MS2(M = Mn,Fe,Co,Ni,Cu) are pyrite structure (Pa3) in p-range. All the V0 are uniformly transformed to same units for comparison purposes.

      表 2  黄铁矿结构CuS2与其他黄铁矿结构过渡金属二硫化物的零压体积模量和晶胞体积对比

      Table 2.  Comparison of zero-pressure bulk modulus (K0) and unit formula volume (V0) of CuS2 with that of other pyrite structure transition-metal disulfides

      图  7  理论计算和衍射实验的CuS2及其它过渡金属二硫化物MS2(M = Mn,Fe,Co,Ni,Cu)的体积-压力变化关系(黑色圆形和黑色三角形分别对应X射线衍射实验和第一性原理理论计算结果。黑色实线和黑色虚线是对实验和计算数据进行Birch-Murnaghan 状态方程拟合结果。蓝色、绿色、紫色、棕色和红色实线分别为文献给出的黄铁矿结构MnS2、FeS2、CoS2、NiS2和CuS2状态方程结果。)

      Figure 7.  Pressure-volume relationship of experimental and calculated CuS2 and other transition metal disulfides MS2(M = Mn,Fe,Co,Ni,Cu) (The black circles and triangles represent the experimental and observed data of the CuS2, respectively. The solid black line and dotted black line are the Birch-Murnaghan equation of state fit respectively with listed parameters. The blue, green, purple, yellow, and red line shows the equation of state of the pyrite structure MnS2, FeS2, CoS2, NiS2, and CuS2 in references.)

      $ \gamma = - \frac{{{\rm{d}}\ln \omega}}{{{\rm{d}}\ln p }} = \frac{{{K_0}}}{{{\omega _0}}}\frac{{{\rm{d}}\omega }}{{{\rm{d}}p}} $

      式中:ωω0分别为各拉曼模在高压和常压的波数;p为压力;dω/dp为各拉曼模频率随压力的变化;K0选择本研究中状态方程拟合结果K0(Exp.)= 99(2) GPa和K0(Calc.)= 85.6(7) GPa。计算得到的各振动模的格临爱森参数值(γ)在表1中给出。

      作为对比,一些黄铁矿结构二硫化物MS2(M = Mn,Fe,Co,Ni)的体积随压力变化曲线也在图7中给出,它们的状态方程参数[14, 27, 41-43]列于表2中。对应压力下,5种黄铁矿结构二硫化物的体积大小关系为VFe < VCo < VNi < VCu < VMn,与King等[23]给出的M—S键长序列LFe—S < LCo—S < LNi—S < LCu—S < LMn—S一致。5种物质的体积模量中,FeS2的体积模量最大,MnS2的体积模量最小,而CoS2、NiS2和CuS2的体积模量之间相差不大。体积模量大小同样与M—S键长关系密切,往往M—S键长越长体积模量越小,则压缩性越强。说明在黄铁矿结构二硫化物MS2(M = Mn,Fe,Co,Ni)中,MS6配位八面体(M—S键长)主导着晶胞体积大小和晶体压缩性,S—S键长的影响较小,而S—S键对拉曼振动频率影响更大。根据CuS2与其他黄铁矿结构二硫化物之间的相似性推测,Cu与其他二硫化物中过渡金属阳离子(Mn2+,Fe2+等)类似,以+2价形式存在。

    • 利用金刚石压腔在不同温压条件下合成了黄铁矿结构CuS2,拓展了已知的CuS2热力学稳定区间。通过显微拉曼光谱和高压同步辐射X射线衍射实验并结合第一性原理理论计算,发现了CuS2在低于30.0 GPa压力范围内保持结构稳定,不发生结构相变,与已知的Hupen等[27]所给结论不同。拉曼光谱研究表明,CuS2所有拉曼模均随压力呈连续单调线性变化,并给出了各个模的线性关系。利用Birch-Murnaghan状态方程对同步辐射X射线衍射实验和第一性原理计算得到的晶胞体积随压力变化数据进行拟合得到,V0(Exp.)= 193.8(5) Å3K0(Exp.)= 99(2) GPa,V0(Calc.)= 196.5(2) Å3K0(Calc.)= 85.6(7) GPa,K0′均固定为4,与Hupen等[27]研究给出的黄铁矿结构CuS2状态方程参数基本吻合。通过和其他黄铁矿结构二硫化物MS2(M = Mn,Fe,Co,Ni)对比发现,MS6配位八面体大小(M—S键长)主导着晶胞体积大小和晶体压缩性。它们之间的一致性变化规律似乎指示黄铁矿结构CuS2中的Cu为+2价。黄铁矿结构CuS2的高温高压稳定性对于认识Cu在地球深部环境中的价态和赋存形式具有重要指示意义。

参考文献 (43)

目录

    /

    返回文章
    返回