Na2CO3的高压拉曼光谱研究
徐超文1,2,3, 高静4,*, 李营1,*, 秦霏5, 刘红1, 易丽1, 崔月菊1, 孙凤霞1, 房雷鸣6
1.中国地震局地震预测研究所, 北京 100036
2.中国工程物理研究院流体物理研究所, 四川 绵阳 621900
3.防灾科技学院, 河北 三河 065201
4.中国科学院地质与地球物理研究所岩石圈演化国家重点实验室, 北京 100029
5.School of Earth Sciences, University of Bristol, Bristol BS81RJ, United Kingdom
6.中国工程物理研究院核物理与化学研究所, 四川 绵阳 621999
*通讯作者 e-mail: gaojing@mail.iggcas.ac.cn; subduction6@hotmail.com

作者简介: 徐超文, 1990年生, 中国地震局地震预测研究所博士后 e-mail: dkchaowen@126.com

摘要

碳酸盐是碳在地球内部的重要载体之一, 其在地幔高温高压条件下的晶体化学是理解地球深部碳的赋存状态和循环过程的关键, 而结构稳定性和相变是晶体化学最基本的研究内容。 碳酸钠(Na2CO3)是一种常见的碱性碳酸盐矿物, 在产自地幔过渡带-下地幔的金刚石中已发现含钠的碳酸盐矿物包裹体, 这成为碳酸钠能够俯冲进入地幔深部的直接矿物学证据。 前人利用拉曼光谱技术研究了Na2CO3在常温常压下的晶格振动模式, 但其在高压下的稳定性和结构变化却鲜有报道。 利用金刚石压腔装置结合先进的共聚焦拉曼光谱技术, 以硅油作为传压介质, 在准静水压力条件下, 在0.001~27.53 GPa压力区间对Na2CO3粉末在600~1 200 cm-1波段的振动特征进行了细致地分析。 本次实验重点分析了[CO3]2-基团振动模式在升压和卸压过程中的行为。 结果表明, 在0.001~11.88 GPa压力范围内, [CO3]2-基团对称伸缩振动 γ1(1 088.06和1 070.76 cm-1)、 反对称伸缩振动 γ3(865.10和797.50 cm-1)和面内弯曲振动 γ4(720.10和696.71 cm-1)都出现了振动峰的分裂。 随着压力增加, 所有振动峰都向高频率漂移, 半高宽也逐渐增加。 在13.40 GPa时, Na2CO3发生结构相变, 具体表现为690.08 cm-1处出现1条新的拉曼峰, 并且随着压力升高该峰的强度逐渐增大。 同时反对称伸缩振动峰 γ3以及面内弯曲振动峰 γ4的强度持续减弱, 半高宽也继续变大。 这些现象表明Na2CO3结构相变源于[CO3]2-内部晶格变化。 当压力卸载到4.18 GPa时, [CO3]2-的振动模式与常温常压下的完全吻合, 相变出现的新峰也已经消失, 表明该相变是由[CO3]2-基团畸变引起的并且具有可逆性。 继续升压至27.53 GPa, 拉曼光谱继续蓝移, Na2CO3的拉曼谱线再没有变化, 说明高压相在这一压强范围内保持稳定。 在整个加压过程中, 反对称伸缩振动 γ3和面内弯曲振动 γ4处的拉曼峰出现强度减弱现象。 同时也计算了各个峰频率对压力的依赖系数d γ/d P, 结果显示[CO3]2-基团内各个振动模式对压力的响应是不同的, 这很可能与C—O键的键长有关。 最后, 对比发现, 对称伸缩振动 γ1峰的强度比反对称伸缩振动 γ3和面内弯曲振动 γ4峰的强度大, 并且[CO3]2-基团对称伸缩振动 γ1受压力影响相对较小, 可以用来区别不同种类的碳酸盐矿物。

关键词: 碳酸钠; 拉曼散射; 高压; 相变
中图分类号:O521+.2 文献标志码:A
High Pressure Raman Spectrum Study of Na2CO3
XU Chao-wen1,2,3, GAO Jing4,*, LI Ying1,*, QIN Fei5, LIU Hong1, YI Li1, CUI Yue-ju1, SUN Feng-xia1, FANG Lei-ming6
1. Institute of Earthquake Forecasting, China Earthquake Administration (CEA), Beijing 100036, China
2. Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China
3. Institute of Disaster Prevention, Sanhe 065201, China
4. State Key Laboratory of Lithospheric Evolution, and Institutions of Earth Science, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
5. School of Earth Sciences, University of Bristol, Bristol BS81RJ, United Kingdom
6. Institute of Nuclear Physics and Chemistry, China Academy of Engineering Physics, Mianyang 621999, China
*Corresponding authors
Abstract

Carbonate is one of the important carriers of carbon in the earth's interior. Therefore, its crystal chemistry under the condition of temperature and pressure corresponding to the mantle is the key to understand the carbon occurrence state and cycle process of deep earth, but structural stability and phase transition are the basic research contents of crystal chemistry. Na2CO3 is a common alkaline carbonate, which enters the earth's interior by subduction of oceanic crust. The existence of sodium carbonate in the subducted slab can significantly reduce the melting temperature of peridotite, promote partial melting and induce mantle heterogeneity. The inclusions of sodium carbonates have been found in the diamonds from the mantle transition zone and the lower mantle, providing direct mineralogical evidence that sodium carbonate can deeply subductin to the deep mantle. The lattice vibration modes of Na2CO3 at ambient condition were reported previously by Raman spectroscopy, but its stability and structural changes under high pressure are poorly reported. In this study, we used silicone oil as pressure transmitting medium and the Raman spectrum of Na2CO3 powder have been carefully ascertained in the pressure range of 0.001~27.53 GPa and the wavelength range of 600~1 200 cm-1 using diamond anvil cell combined with advanced confocal Raman spectroscopy. This experiment focused on the analysis of the behavior of [CO3]2- group vibration mode in the process of compression and decompression. The results showed that splitting of vibration peaks respectively appeared in symmetric stretching vibration γ1, antisymmetric stretching vibration γ3 and the in-plane bending vibration γ4 of the [CO3]2- at the pressure range of 0.001~11.88 GPa. With the increase of pressure, all peaks shift to high frequency, and the full width at half maximum (FWHM) increases gradually. The phase transition occurred at 13.40 GPa, accompanied by a new Raman peak at 690.08 cm-1 and the intensity of the peak increases with the increase of pressure. At the same time, the intensity of antisymmetric stretching vibration and in-plane bending vibration continued to weaken, and the FWHM of Na2CO3 also continued to increase, indicating that the phase transition of Na2CO3 originates from the internal lattice vibration of [CO3]2-. When the pressure is decompressed to 4.18 GPa, it is found that the vibration mode of [CO3]2- is identical with that at ambient condition, and the new peak has disappeared, indicating that the phase transition is caused by the distortion of [CO3]2- group and is recoverable. The Raman peaks continued shifting to high frequencies when the pressure increased to 27.53 GPa, suggesting this new phase can remain stable in this pressure range. The intensity of Raman peaks at the antisymmetric stretching vibration γ3 and in-plane bending vibration γ4 decreased during compression. Meanwhile, the calculated dependence coefficient of relative pressure-shift of each Raman peak showed that the response of each vibration mode to pressure is different in [CO3]2-. This is probably related to the length of the C—O bond. Finally, by comparison, the intensity of symmetric stretching vibration γ1 peak is higher than that of antisymmetric stretching vibration γ3 and in-plane bending vibration γ4. The pressure also has little effect on the typical Raman peak γ1 of [CO3]2-, and therefore can be used to distinguish different kinds of carbonates.

Keyword: Sodium carbonate; Raman scattering; High pressure; Phase transition
引言

高温高压实验和理论模拟计算的结果显示, 洋壳中的碳酸盐矿物可以随着板块俯冲作用进入地球深部[1, 2, 3]。 来自过渡带-下地幔的金刚石中发现碳酸盐的包裹体, 为碳酸盐能被带入到下地幔提供最直接的矿物学证据[4, 5, 6], 此外, 研究表明, 地球深部碳酸盐与超深金刚石的成因密切相关[7]。 因此, 研究碳酸盐矿物在高温高压条件下的稳定性, 对于理解碳在地球深部的分布和赋存以及深部碳循环动力学过程都具有重要的指示意义。

碳酸钠Na2CO3是一种常见的碳酸盐矿物, 在室温条件下为单斜结构(γ相), 空间群为C2/m, Z=4, 晶胞参数为a=8.920 Å, b=5.245 Å, c=6.050 Å, V=277.50 Å3, β =101.35°。 常明等[8]利用红外光谱对Na2CO3的结构进行研究, 发现[CO3]2-基团存在4种红外吸收模式。 目前, 还没有通过拉曼光谱研究Na2CO3在高压下的稳定性和结构的实验报道。 为此, 本文利用金刚石压腔(diamond anvil cell, DAC)装置结合原位激光拉曼光谱技术, 研究了Na2CO3在室温高压下的拉曼光谱特征。 该研究结果有助于更好的理解碳酸盐在地幔的分布状况。

1 实验部分

本次高压实验采用Mao-Bell型金刚石压腔装置, 金刚石砧面直径为400 μm, 高压密封垫片为铼片, 首先将其预压到厚度~40 μm, 然后在凹坑中心钻一个直径为200 μm的圆孔作为样品腔, 接着将带孔的铼片按照预压时的方位放置在金刚石顶砧上, 将Na2CO3粉末(Alfa Aesar)压成薄片后放入样品腔内, 红宝石作为标压物质, 硅油作为传压介质。 拉曼光谱分析在中国科学院地质与地球物理研究所岩石圈演化国家重点实验室进行, 所用仪器为共聚焦的WITec 300R激光拉曼光谱仪, 实验装置如图1所示。 采用波长为532 nm激光器激发样品, 激光发射功率为50 mW。 实验重点测定了600~1 200 cm-1波段, 扫描时间为60 s, 扫描次数为1次。

图1 金刚石压腔结合拉曼散射示意图Fig.1 Sketch of diamond anvil cell combined with Raman scattering

2 结果与讨论

拉曼效应是基于光与物质内化学键或晶格之间的相互作用, 每种化学成分的拉曼光谱都是独一无二的, 通过拉曼光谱可对材料进行定性和定量研究分析。 本次实验重点关注Na2CO3中的碳酸根离子[CO3]2-的内振动模式。 [CO3]2-为平面构型, C原子在3个O原子围成的三角形中心, [CO3]2-基团有3个拉曼活性的振动模: (1)对称伸缩振动γ1, (2)反对称伸缩振动γ3和(3)面内弯曲振动γ4。 本次实验的最高压力为27.53 GPa, 共采集了22条升压拉曼谱和1条降压谱。 图2(a)是常温常压下Na2CO3的拉曼光谱, 图2(b)是常温条件下, Na2CO3的拉曼光谱随着压力变化的情况, 图2(c)是降压过程中拉曼光谱。

图2 不同压力下Na2CO3的原位拉曼光谱Fig.2 In-situ Raman spectra of Na2CO3 at different pressure

2.1 常温常压下Na2CO3拉曼光谱特征

从图2(a)可以看出: 由于本次实验所用的WITec共聚焦拉曼光谱技术具有高的分辨率, 发现Na2CO3原始晶胞内, [CO3]2-基团在同一振动模内出现了波数相近的双峰: 对称伸缩振动γ1, 1 088.06和1 070.76 cm-1; 反对称伸缩振动γ3, 865.10和797.50 cm-1; 面内弯曲振动γ4, 720.10和696.71 cm-1。 相同振动模式区域拉曼峰的出现与消失可能是由于晶体的不同取向导致的, 并非全部的振动模式能都出现, 而Na2CO3振动模式γ1, γ3γ4区域内相近波数的拉曼双峰出现, 是由于晶体结构中存在八配位的Na+与O2-组成相互共面的八面体与七配位的Na+与O2-组成多面体, Na—O键的细微变化, 带动局域C—O键长的改变, 最终产生波数相近的拉曼双峰。

2.2 高压下Na2CO3拉曼光谱特征

图2(b)是常温高压条件下Na2CO3的拉曼光谱图随着压力变化的情况。 从图中可以明显的看出, 随着压力增大, 6条拉曼峰均向高波数方向移动, 这是因为高压下, C—O键受压缩变短, 使振动频率增加。 但由于Na2CO3晶体结构具有无序性, 在移动过程中, [CO3]2-反对称伸缩振动γ3和面内弯曲振动γ4产的的拉曼峰强度逐渐变弱, 半高宽也变大。 但是, 在0~11.88 GPa压力范围并没有观察到新峰的出现与旧峰的消失, 说明Na2CO3没有发生结构变化。

当压力达到13.40 GPa时, 拉曼光谱发生细微的变化, 我们发现在690.08 cm-1出现一条新拉曼峰(红色箭头处), 表明Na2CO3结构相变源于[CO3]2-内部晶格振动。 虽然观察到新峰出现, 但是原来[CO3]2振动产生的6条拉曼峰均稳定存在, 暗示着可能是由于Na在晶格中的移动导致八面体或者多面体发生畸变, 使[CO3]2-振动改变, 从而发生压致相变, 同时也观察到反对称伸缩振动γ3(865.10和797.50 cm-1)和面内弯曲振动γ4(720.10和696.71 cm-1)强度持续减弱, 半高宽也继续变大。 研究发现[9], 同为碱金属碳酸盐的BaCO3在10.6 GPa附近发生相变, 而本次实验观察到Na2CO3相变压力高于BaCO3。 因为Na+离子半径(Ⅶ -1.12 Å, Ⅷ -1.18 Å)小于Ba2+(Ⅵ -1.35 Å), 所以Na2CO3晶体排列更加致密, 所需的相变压力也更高。 截至目前, 还没有高压下Na2CO3晶体的X射线衍射、 拉曼光谱、 红外光谱等其他相关光谱学分析数据, 因此, 其在相变后具体结构还有待进一步研究。

相变之后, 当压力继续升高, 位于反对称伸缩振动γ3区域内右侧的峰强度已经非常微弱, 原面内弯曲振动γ4区域内的两条峰强度持续变弱, 半高宽也变大, 导致720.10和696.71 cm-1出现重叠现象。 但总体的看, 所有峰持续发生不同程度的蓝移。 值得注意的是, 在压力达到20.36 GPa以后, 位于面内弯曲振动γ4的新峰继续向高波数移动, 而原来的两条拉曼峰开始红移, 表明随着压力增大面内弯曲振动出现软化现象。

当压力卸载到4.18 GPa时, 我们采集了拉曼光谱[图2(c)], 发现此压力下[CO3]2-振动模式与常温常压下的完全吻合, 相变出现的新峰也已经消失, 所有拉曼峰的频率随着压力的降低而降低, 表明Na2CO3晶格振动导致的结构相变在低压条件下是可逆的。 利用PeakFit软件, 我们拟合出了不同压力下的拉曼峰峰位列于表1, 由于高压下部分拉曼峰强度减弱、 宽化、 重合, 表1并没有全部列出高压下所有拉曼峰的峰位。

表1 室温Na2CO3在不同压力下的拉曼光谱图 Table 1 Raman spectra of Na2CO3 at room temperature under different pressure

为进一步理解高压下Na2CO3晶格中[CO3]2-振动模式, 我们将不同振动模式对应的频率与压强关系绘制在图3中。 同一种符号代表相同的振动模。 整体来看, 随着压力增加, 几乎所有的峰斜率都是正的, 暗示着各个模式振动频率增加。 当压力升高到13.40 GPa(图中长虚线), 新出现红色菱形的拉曼峰, 表明发生结构相变; 当压力升高到20.36 GPa(图中短虚线), γ4振动模式区域绿色菱形的频移降低, 斜率为负值, 发生软化现象。 线性拟合结果显示, 不同振动模式对压力的依赖系数(即dγ/dP值)分别为: dγ1/dP=3.08 cm-1·GPa-1(实心方块), dγ1/dP=3.14 cm-1·GPa-1(空心方块), dγ3/dP=2.18 cm-1·GPa-1(实心圆形), dγ3/dP=2.73 cm-1·GPa-1(空心圆形), dγ4/dP=3.78 cm-1·GPa-1(实心菱形, 小于20.36 GPa), dγ4/dP=-1.33 cm-1·GPa-1(实心菱形, 大于20.36 GPa), dγ4/dP=2.49 cm-1·GPa-1(绿色空心菱形, 小于20.36 GPa), dγ4/dP=1.02 cm-1·GPa-1(绿色空心菱形, 大于20.36 GPa), dγ4/dP=1.74 cm-1·GPa-1(红色菱形), 可见[CO3]2-基团内各个振动模式对压力的响应是不同的, 这很可能与C—O键的键长有关。

图3 不同压力下Na2CO3的振动模式Fig.3 Pressure dependence of all vibrational models of Na2CO3

整体来看, 对称伸缩振动γ1峰的强度比反对称伸缩振动γ3和面内弯曲振动γ4峰的强度大, 并且对称伸缩振动γ1的半高宽只是随压力升高略微增加, 本次实验还发现相转变对[CO3]2-对称伸缩振动γ1影响相对较小。 因此, γ1处的特征峰可以用来区别不同种类的碳酸盐矿物。

3 结论

(1)在室温条件下, 由常压增加到13.40 GPa时, Na2CO3发生相变, 在此之后到最高27.53 GPa再未出现相变, 该相变是可逆的。

(2)在本次实验压力范围内, 随压力升高, 除了面内弯曲振动γ4区域的两条峰从20.36 GPa处开始发生软化之外, 其余峰均发生蓝移并且压力对不同振动模式产生的影响不同。

(3)[CO3]2-对称伸缩振动γ1可以用来区别不同种类的碳酸盐矿物。

参考文献
[1] Kang N, Schmidt M W, Poli S, et al. Contributions to Mineralogy and Petrology, 2016, 171(8-9): 74. [本文引用:1]
[2] Ghosh D B, Bajgain S K, Mookherjee M, et al. Sentific Reports, 2017, 7(1): 848. [本文引用:1]
[3] Tappe S, Romer R L, Stracke A, et al. Earth and Planetary Science Letters, 2017, 466: 152. [本文引用:1]
[4] Kaminsky F V, Ryabchikov I D, Wirth R. Mineralogy and Petrology, 2016, 110: 387. [本文引用:1]
[5] Drewitt J W E, Walter M J, Zhang H, et al. Earth and Planetary Science Letters, 2019, 511: 213. [本文引用:1]
[6] Cerantola V, Bykova E, Kupenko I, et al. Nature Communications, 2017, 8: 15960. [本文引用:1]
[7] Martirosyan N S, Litasov K D, Lobanov S S, et al. Geoscience Frontiers, 2019, 10(4): 1449. [本文引用:1]
[8] CHANG Ming, YUN Hai-li, DONG Si-si, et al(常明, 郧海丽, 董思思, ). Infrared Technology(红外技术), 2016, 38(9): 803. [本文引用:1]
[9] MUBARAK Molutjan, ANWAR Hushur, WANG Jing, et al(穆巴拉克·木里提江, 艾尼瓦尔·吾术尔, 王静, ). Materials Reports(材料导报), 2019, 33(24): 4062. [本文引用:1]