许思维 杨晓宁 杨大鑫 王训四 沈祥
1) (湖南文理学院数理学院, 常德 415000)
2) (宁波大学高等技术研究院, 浙江省光电探测材料及器件重点实验室, 宁波 315211)
本文制备了硫系玻璃Ge11.5As24Se64.5-xSx (x = 0, 16.125%, 32.25%, 48.375%和64.5%)并研究了其光学性质, 目的在于筛选可用于光学器件的最佳组分.通过测试该系列玻璃的激光损伤阈值、折射率、三阶非线性折射率以及吸收光谱, 结果发现, 玻璃中的Se被S原子逐渐替代后, 玻璃的线性和三阶非线性折射率逐渐降低, 玻璃光学带隙和激光损伤阈值不断升高.我们进一步利用拉曼散射光谱和高分辨率X射线光电子能谱研究导致这些物理性能变化的结构起源, 通过分析玻璃中不同结构单元的演变过程, 发现在这些玻璃网络结构中均以异极键(Ge-Se/S, As-Se/S)为主, 且相对于Se而言, Ge和As优先与S结合成键.随着玻璃结构中S/Se比例的增加, 与Se相关的化学键(Ge-Se, As-Se和Se-Se)数量逐渐减少, S相关化学键(Ge-S,As-S和S-S)数量逐渐增加, 但这对玻璃的拓扑结构几乎没有影响.由此可以断定引起玻璃物理性质变化的主要原因是玻璃结构体系中各个化学键强度之间的差异.
硫系玻璃因其自身具有极短的非线性光学响应时间、低的声子能、高的非线性系数和良好的红外透过区域而受到众多学者的广泛关注[1-6], 是在光学存储、非线性光学、红外激光和生物传感等中红外光学领域具有广阔应用前景的基质材料[7-13].例如, 在其高折射率(n)和低声子能方面, 近年来被用作中红外光纤的核心材料, 并成为拉曼增益介质的理想候选材料.此外, 硫系玻璃一般具有宽的玻璃形成区域, 因此可以通过调整玻璃化学组分来获取其某些特定的光学性能.
为了更好地满足新型非线性波导器件应用的需求, 硫系玻璃必须具备三个重要条件: 1) 高的玻璃 转 变 温 度(glass transition temperature, Tg);2) 高的三阶非线性折射率(n2); 3) 较大的激光损伤阈值(laser damage threshold, Ith).与二元硫系玻璃As-S和As-Se相比, 三元硫系玻璃Ge-As-S和Ge-As-Se具有更高的Tg.一般而言, 随着四配位键的Ge原子引入, 可使得二元体系As-S和As-Se玻璃网络结构的三维立体程度更高, 从而提高其Tg, 并且可以预期在更大范围内调节三元玻璃的物理性能[1,2].对于三阶非线性而言, Ge-As-Se和As2S3玻璃是硅基玻璃的成百上千倍[14].最近,对三元体系Ge-As-Se玻璃的研究表明, 当玻璃的平均配位数(mean coordination number, MCN)在2.45-2.5之间时, 玻璃具有最小的玻璃转变活化能和脆性指数, 因此得以证明这类玻璃具有最小的结构弛豫动力[15].但是, 在Ith方面, Se基玻璃的表现却有点差强人意.当激光功率一旦超过其损伤阈值, Se基玻璃表面便会出现烧伤、蒸发和破斑等灾难性损伤, 而在同样激光功率的情况下, S基玻璃却很难观察到[16].这主要是归因于S基玻璃相比于Se基玻璃具有更大的平均化学键强.因此与Se基玻璃相比, S基玻璃表现出更高的激光耐受程度[12].由于玻璃Ge11.5As24Se64.5-xSx的MCN均为2.47, 因此它的热稳定性可以很好地应用于波导器件.为了筛选最佳的光学应用组分, 有必要了解玻璃的Ith和n2等物理参数随玻璃组分配比中Se/S比例变化的演变情况.再者, 从化学元素替代的角度出发了解S取代Se的过程中玻璃的结构是如何改变同样也是很有趣的事情.
本文选用Ge11.5As24Se64.5作为起始材料, 逐步将玻璃结构体系中的Se原子用S原子取代, 对其演变过程中光学性质(Ith, n, n2和吸收光谱)进行了对比研究.同时借助拉曼散射光谱(Raman scattering spectra)和X射线光电子能谱(X-ray photoelectron spectroscopy, XPS)技术系统而深入地研究其微观结构中结构单元和化学键(同极键和异极键)的变化情况, 进一步分析了物理性能变化与玻璃结构参数之间的关系.
采用熔融-淬冷技术制备了Ge11.5As24S64.5-xSex样品玻璃.关于玻璃制备的流程见图1所示.
图1 Ge11.5As24Se64.5-xSx玻璃制备的流程图Fig.1.Flow chart of Ge11.5As24Se64.5-xSx glasses fabrication.
值得注意的是在制备过程中, 按照玻璃不同的化学成分称重Ge, As, S和Se四种高纯度(5N)的化学单质, 总计20 g; 而对于石英管需要在已配置好的王水中浸泡2-3 h, 随后取出用去离子水和无水乙醇对其反复冲洗, 晾干之后放置于温度为180 ℃的干燥箱中干燥; 将称量准确的原材料置于石英管中, 用真空泵将其抽至10-4Pa以下,然后用氢氧焰封管后放入摇摆炉中, 升温至850 ℃,使其均匀熔化, 熔制时间控制在10 h以上; 采用水冷法对其进行淬冷后退火, 退火温度控制在Tg以下20 ℃以下, 时间为3 h左右.
首先, 采用X射线衍射(X-ray diffraction, XRD)和红外(Infra-red)显微镜对玻璃样品的质量进行了检测评估.XRD谱图没有显示任何尖锐的衍射峰, 表明玻璃结构为非晶态, 红外显微镜显示所有样品中没有任何可见气泡.用X射线能谱仪成分分析(energy dispersive spectrometer, EDS)检测了玻璃样品的组份(考察了五个不同的点位).总的来说, 实际制备的玻璃样品成分与称重的成分之间的差别均小于0.5%.为了获得玻璃样品的Ith, 将其放置于波长为830 nm激光光束之下照射(保证激光聚焦光斑的直径为1 µm), 照射时间为2 min,之后检测玻璃表面是否存在由于激光辐射引起的表面损伤.由于激光功率计在功率为0.1 mW时表现非常不稳定, 所以对于富Se玻璃的不确定的程度影响较大.玻璃的n采用Metricon Moder 2010棱镜耦合器测量, 激光光源为1550 nm, 误差小于0.001.玻璃的三阶非线性采用是Z-扫描技术, 样品的厚度为1-2 mm, 所用的激光器为钛:蓝宝石飞秒激光器(Clark-MXR CPA 2001型), 泵浦源为Quantronix Palitra系列的光参数放大器, 激发光源的波长为1.55 µm, 脉冲宽度为260 fs, 重复频率为1 kHz.通过紫外-可见-近红外分光光度计(Varian Cary 5000型)测试了玻璃样品的吸收光谱, 波长范围为400-1500 nm, 采用经典的Tauc方程对所测得的吸收光谱进行计算, 得出玻璃样品各自的光学带隙(Eg).
拉曼散射光谱是由配备液氮冷却CCD探测器的显微拉曼光谱仪(Jobin-Yvon-Horiba T64000型)记录的.采用830 nm激光线作为激励源, 并尽可能减小激光功率以避免光诱导效应.光谱仪的分辨率约为0.5 cm-1.为了消除热效应对抽取不同结构单元振动信息的影响, 去除原始数据中的玻色-爱因斯坦(Bose-Einstein)热因子, 其原理如下[17]:
式中, ω和ω0分别为拉曼频率(cm-1)和激发光源的绝对频率(cm-1); Imea(ω)和Ire(ω)分别是测试的拉曼强度和扣除热效应后的拉曼强度; n是玻色-爱因斯坦(Bose-Einstein)因子, 其公式如下:
其中ħ, kB和T分别为普朗克常数、波尔兹曼常数和拉曼散射光谱测量时的温度.利用PeakFit软件将记录的拉曼散射光谱进一步分解为不同的结构单元, 并通过迭代法模拟得出每个拉曼振动峰的相关信息, 即位置(position)、半峰宽度(width at half maximum)和强度(intensity)等.
采用X射线光电子能谱仪(Scienta ESCA-300型)记录XPS光谱, 由单色Al KαX射线(1486.6 eV)作为激发源, 其激发功率为250 W.对于所有的测量, 玻璃样品表面和探测器之间的角度都是90º.为了获得具有代表性的数据需要进行以下两各方面的操作.第一, 在分析腔内的玻璃试样在真空度为1.0 × 10-8Pa的压力下破损表面一层, 这样可以消除由于氧或其他污染物对玻璃样品表面的影响; 第二, 在玻璃样品不同的破损区域重复性测量,并进行检查比对筛选.在重复测量过程中需注意采用低能电子枪(< 10 eV)对由于光电子发射产生的玻璃样品表面电荷进行中和, 至于这个电势的大小取决于玻璃样品表面的几何形状和电导率.Ge3d, As 3d, Se 3d和S 2p光电子能谱的数据分析采用标准的Casa-XPS软件包进行.为了能更准确地描述真实的原子核内层电子的光电能谱, 扣除掉Shirley背景, 并采用卷积Voigt线形峰进行拟合.每个Ge 3d, As 3d, Se 3d和S 2p核级谱都包含一个或多个自旋轨道双峰, 由Ge 3d5/2-3d3/2,As 3d5/2-3d3/2, Se 3d5/2-3d3/2和2p3/2-2p1/2组分组成, 通过迭代曲线拟合确定XPS特征峰内部双峰的个数.
Ge, As和S(Se)原子在硫系玻璃中的配位数分别是4, 3和2.因此, 在理想情况下, 一个Ge原子与2个S (Se)原子结合, 一个As原子与1.5个S (Se)原子结合, 那么本文中Ge和As所需要的S (Se)原子数量为59个.然而对于所制备的Ge11.5As24Se64.5-xSx玻璃成分而言, S (Se)原子数量超过了此数值, 故所有玻璃样品均存在轻微的富S(或富Se)[18].表1列出了所有玻璃样品的组分、n, Ith,Eg和n2.
表1 Ge11.5As24Se64.5-xSx的 组 分 与 光 学 参 数(n, Ith,Eg和n2)Table 1.Compositions and optical parameters (n, Ith, Eg and n2) of Ge11.5As24Se64.5-xSx glasses.
结果发现, Ith和Eg随着玻璃组分中S含量的增加而增大, 同时玻璃的n和n2随之降低.对于Ge11.5As24S64.5玻璃, 在激光器能够达到的最大输出功率均未有损伤的出现, 但所测得的Ith与之前同领域的研究结果的数据相当[2,12].
众所周知, 材料的结构决定其性质.因此, 了解玻璃的微观结构, 并研究其组分的变化对结构的影响具有重要意义.对于缺乏长程有序的非晶硫系玻璃材料而言, 拉曼散射光谱测试技术被认为是了解玻璃结构的局部化学环境最有利的手段之一.
Ge11.5As24Se64.5-xSx玻璃的拉曼散射光谱如图2所示.由于玻璃成分中S或Se原子数量略微大于化学配比的组分, Ge和As之间所形成的阳离子同极键(Ge-Ge和As-As)的数量则可以忽略.由此可以认为在Ge11.5As24Se64.5-xSx玻璃网络结构中, 主要由GeSe4/2和GeS4/2四面体结构(共边和共角)、AsSe3/2和AsS3/2三角锥结构、以及S(Se)原子与自身结合形成的S(Se)链状或环状结构组成.
图2 Ge11.5As24Se64.5-xSx玻璃的拉曼散射光谱分峰拟合图Fig.2.Raman scattering spectra of Ge11.5As24Se64.5-xSx glasses and their decompositions.
对于三元体系的Ge-As-Se硫系玻璃而言, 其内部结构中共顶点(corner-sharing, CS)GeSe4/2和共边(edge-sharing, ES)GeSe4/2四面体结构单元所对应的拉曼振动峰位置分别为198 cm-1和215 cm-1[19,20].AsSe3/2三角锥结构单元对应的拉曼振动峰为230 cm-1, As-Se键振动峰为265 cm-1, Se链状或环状结构单元(Se chains or rings)内阴离子同极键Se-Se的振动峰为255 cm-1左右[19,21-23].另一方面, 在硫系玻璃Ge-As-S中, 344和370 cm-1附近的拉曼特征峰分别归因于GeS4/2(CS)和GeS4/2(ES)四面体结构单元[24-26].三角锥结构单元AsS3/2中的As-S键振动出现在338 cm-1附近, 而异极键As-S引起了380 cm-1附近的拉曼振动峰的出现[26,27].S链状或环状单元结构中的阴离子同极键S-S的拉曼振动峰在490 cm-1附近出现[26,28,29].
为了对Ge11.5As24Se64.5-xSx玻璃微观结构中各个单元结构的演变趋势进一步量化, 将这些玻璃的拉曼散射光谱进行分峰拟合, 通过分解所得到的每条谱线, 分别对应不同的结构单元, 如图2所示.各个分解峰的积分面积与整个拉曼散射光谱积分面积的相对比例如表2所列.结果表明, 研究对象Ge11.5As24Se64.5-xSx玻璃结构体系中均以异极键Ge-Se/S和As-Se/S为主导地位.总体情况是随着玻璃组分中S原子含量的增加, 与S原子结合形成的化学键(Ge-S, As-S和S-S)的数量也逐渐增加; 与此同时, 与Se原子结合形成的化学键(Ge-Se, As-Se和Se-Se)的数量却逐渐减少.但是, 值得关注的是在玻璃样品Ge11.5As24Se48.375S16.125结构中的同极键S-S被完全抑制, 该情况表明在玻璃Ge11.5As24Se64.5-xSx结构体系中, 相比Se和S原子, Ge和As原子更加优先与S原子结合成键,形成结构单元; 这一点从原子电负性的角度也可以得到解释, 原子之间电负性的差异越大, 越容易形成化学键.由于原子Ge, As, Se和S的电负性分别为2.01, 2.18, 2.48和2.58, 所以Ge或As原子与S原子之间的结合能要高于其二者与Se原子之间的结合能[30-33].
表2 拉曼散射光谱分峰拟合中各个结构单元的相对比例Table 2.Relative ratio of the different structural units derived from the decomposed Raman scattering spectra.
利用X射线光电子能谱仪对Ge11.5As24Se64.5-xSx玻璃表面的结构信息(化学状态和化学键)进行了测量, 并分解了Ge 3d, As 3d, Se 3d和S 2p的XPS光谱.由于Ge, As和S(Se)的配位数分别为4,3和2, 过配位和欠配位原子的数目很少.因此,Ge和As原子在电负性方面的差异可以忽略, 将S 2p光谱分解为三种不同的近邻构型, 按照结合能依次递减的顺序分别对应S-S-S, S-S-Ge/As和As/Ge-S-Ge/As; 与S 2p光谱类似, Se 3d光谱按照结合能从高到低分解为Se-Se-Se, Se-Se-Ge/As和 As/Ge-Se-Ge/As三种可能的成键结构,分别如图3(a)和图3(b)所示.
图3 (a) Ge11.5As24Se64.5-xSx 玻璃的S2p 的XPS分解; (b) Ge11.5As24Se64.5-xSx 玻璃的Se3d 的XPS分解Fig.3.(a) S2p spectra of Ge11.5As24Se64.5-xSx glasses and their decompositions; (b) Se3d spectra of Ge11.5As24Se64.5-xSx glasses and their decompositions.
对于前者S 2p光谱, 在大约162.3和162.0 eV左右的位置分别可以获得一个较高结合能的S-SGe/As结构单元和一个低结合能的As/Ge-SGe/As结构单元[30,34].在Ge11.5As24Se64.5-xSx玻璃中, 所有的S 2p光谱均以As/Ge-S-Ge/As结构为主, 并没有发现S-S-S结构单元的存在.随着Se/S比例的增加比例的增加, As/Ge-S-Ge/As结构单元的积分面积随之增大, 而S-S-Ge/As结构单元的积分面积逐渐减小, 最终消失于Ge11.5As24Se48.375S16.125玻璃结构中.而关于后者Se 3d光谱, 高结合能的Se-Se-Ge/As结构单元和低结合能的As/Ge-Se-Ge/As结构单元分别出现在54.9 eV和54.5 eV位置附近, 在所有光谱中, 主要是以As/Ge-Se-Ge/As结构单元的形式存在, 只有少量的Se-Se-Ge/As结构单元.在任何玻璃样品中, 均未发现有Se-Se-Se结构单元的出现.随着Se/S比例的减小, As/Ge-Se-Ge/As结构单元的积分面积逐渐增加, Se-Se-Ge/As结构单元的积分面积逐渐减小.
关于Ge 3d和As 3d光谱, Ge原子一般存在四种最近邻配置构型, 分别是Ge-4Se/S, Ge/As-Ge-3Se/S, 2Ge/As-Ge-2Se/S和3Ge/As-Ge-Se/S成键结构.再次忽略Ge和As的电负性的差异.根据第一原理分子动力学模拟可知, 1个Ge原子与其他四个Ge原子结合成键的情况是几乎不可能[35].为了便于分析结构单元的变化情况, 故将所有的Ge3d光谱分解成完美的GeSe4/2四面体结构和Ge-Ge-相关结构这两种结构单元.另一方面,As原子有三种不同的最近邻构型, 即As-3Se/S,As/Ge-As-2Se/S和2As/Ge-As-Se/S.与Ge 3d谱类似, 同样将所有的As 3d谱分为AsSe3/2三角锥结构和As-As-相关结构.Ge 3d和As 3d谱及其分解峰分别对应图4(a)和图4(b).其中位于31.3和42.9 eV附近分别对应GeSe/S4/2四面体和AsSe/S3/2三角锥结构单元[31,32].在拟合过程中加入任何其他结构单元都无法改善其拟合的结果, 由此说明阳离子之间所形成的同极键数量可以忽略不计, 这与拉曼散射光谱分析的结果相一致.
图4 (a) Ge11.5As24Se64.5-xSx 玻璃的Ge3d 的XPS分解; (b) Ge11.5As24Se64.5-xSx 玻璃的As3d 的XPS分解Fig.4.(a) Ge3d spectra of Ge11.5As24Se64.5-xSx glasses and their decompositions; (b) As3d spectra of Ge11.5As24Se64.5-xSx glasses and their decompositions.
XPS谱中的各个不同的结构单元及其对应的结合能(binding energy, BE)、半高全宽(full width at half maximum, FWHM)和各结构单元的积分面积占XPS峰总面积的百分比(content), 如表3所列.
从S和Se基结构单元的逐渐演化情况表明,由于S和Se都是2配位, 原子替代对玻璃拓扑结构的影响微乎其微.这与表2中不同结构单元的比例随着S含量的增加而逐渐变化的结果相吻合.分析了表2和表3中S-S和Se-Se键的演变过程,可以发现S原子与Ge和As原子的结合成键数量不同于Se原子和二者(Ge和As)结合成键的数量.虽然在Ge11.5As24Se64.5-xSx玻璃体系内S和Se原子的总量相同, 但是从拉曼散射光谱和XPS光谱分析结果来看, 同极键S-S的数量相对比值要小于Se-Se键.此外, 在Ge11.5As24Se48.375S16.125玻璃样品中没有发现S-S键的存在, 而Se-Se键仍然出现在Ge11.5As24Se16.125S48.375玻璃样品结构体系中.同极键S-S在玻璃Ge11.5As24Se48.375S16.125结构中被完全抑制, 这一情况可以说明在Ge11.5As24Se64.5-xSx玻璃结构中, S要比Se原子更容易与Ge和As原子结合成键, 这与它们之间电负性差异有关(在之前分析拉曼散射光谱时已作过类似的解释).
表3 Ge11.5As24Se64.5-xSx 玻璃的Ge3d, As3d, Se3d 和S2p 的XPS的拟合参数Table 3.The fitting parameters for the decomposed Ge3d, As3d, Se3d and S2p spectra of Ge11.5As24Se64.5-xSx glasses.
这里需要指出, 无论表2中各个拉曼分解峰的积分面积所表示的结构单元的相对数值, 还是表3中XPS各分解峰对应的结构单元的相对比例均不等于玻璃结构中各个结构单元的实际数量.其主要因素是由于不同实验方法的探测范围和各个结构单元的敏感因子不同, 从而导致拉曼散射光谱和XPS的测试结果中每个结构单元的相对数量不完全相同.但是, 我们可以从这两种测试手段的分析结果中发现, 在玻璃结构中各个结构单元随成分变化的演变趋势是一致的.由于元素的置换导致玻璃拓扑结构的变化可以忽略不计, 这更容易解释n和Ith等物理性质的变化.
硫系玻璃n主要取决于玻璃内部的离子极化率和玻璃密度[36].首先, 在Ge11.5As24Se64.5-xSx玻璃体系原子总数量不变的情况下, 离子极化率较高的Se(12.57 Å3)逐渐被离子极化率较低的S(9.72 Å3)所取代, 导致玻璃内部体系中各个离子的极化程度总体降低.其次, 在玻璃结构保持一致的情况下,随着Se原子逐渐被S原子所取代, 相对原子质量较小的S原子的增加会导致玻璃密度减小.综上所述, 随着玻璃体系中S含量的增加, n几乎呈线性下降.根据米勒的预测, 高n的玻璃通常表现出高的三阶非线性特性, 换而言之, 玻璃n的高低决定了其n2的大小[37].随着S原子含量的增加, Ith的增加是由于玻璃中形成了更多更强的化学键.化学键Ge-S(133.0 kcal/mol), As-S(90.7 kcal/mol)和S-S(101.6 kcal/mol)的键强度均高于化学键Ge-Se(115.8 kcal/mol), As-Se(22.9 kcal/mol)和Se-Se(79.0 kcal/mol)[33].
本文所研究的Ge11.5As24Se64.5-xSx玻璃具有拓扑结构相似而物理性质不同的特点.这有利于更好地使用该类材料制备平面波导或光纤器件.随着Se含量的增加, n几乎呈线性增加, 可以帮助我们设计玻璃组分, 使其具有理想n.另一方面, 通过改变玻璃成分可以调控其n2和Ith等物理特性.例如, Se基玻璃无法很好的抵御激光辐射; 而S基玻璃则相反, 但是n2相对较低.所以, 在玻璃中混合配比Se和S原子, 可以取长补短, 在降低Ith的前提下增加其n2.
本文以高非线性Ge11.5As24Se64.5玻璃为起始材料, 并将S原子逐步替代Ge11.5As24Se64.5玻璃中的Se原子, 以提高玻璃的Ith.制备了一系列Ge11.5As24Se64.5-xSx玻璃, 并对其物理性能进行测试与分析.从中发现, 玻璃结构体系中Se被S取代不仅改变了玻璃的Eg, n和n2, 而且玻璃的Ith得以显著改善.进一步利用拉曼散射光谱和高分辨率XPS来探测玻璃的结构.结果表明, 异极键(Ge-Se/S, As-Se/S)在玻璃中占主导地位, S原子取代Se原子只会使得与Se相关的化学键(Ge-Se, As-Se和Se-Se)逐渐演变少, 与S相关的化学键(Ge-S, As-S和S-S)逐渐变多;与Se基化学键相比, S基化学键的强度更大, 从而造成Ge或As原子先与S原子结合成键.因为Se和S都同属二配位, 所以在玻璃内, S与Se的相互替换对于玻璃拓扑结构的变化几乎没有影响.由此表明, 玻璃结构体系中各个化学键之间的强度差异是引起玻璃物理性质变化的主要因素.