王雨思,周贤建,李青原,叶钊松,蔡源海,刘近秋,陈 楠,卜轶坤
(1. 厦门大学 电子科学与技术学院(国家示范性微电子学院),厦门 福建 361000;2. 瑞之路(厦门)眼镜科技有限公司,厦门 福建 361000)
随着电脑、智能手机、移动显示设备等使用人群的不断增多,屏幕的使用时长不断增加,数码时代的用眼健康问题日益被关注。相关研究表明,对人眼有害的蓝光波长范围为415~455 nm,高能蓝光不仅会对眼睛产生不可逆转的伤害,同时会引发眼底视网膜启动光氧化机制。长期的蓝光照射会导致视网膜色素上皮细胞衰亡,引起视网膜病变、黄斑变性、白内障、青光眼等眼科疾病[1]。另一方面,460~480 nm左右的蓝光光谱波长对于人体生物钟调节起着非常重要的作用。近年来,对于如何有效防治蓝光损伤的同时,也不影响蓝光衰减对于人体节律调节的影响研究,成为热点课题,在商用光源及防护镜片领域受到了广泛重视[2]。
目前,有效的蓝光防护技术路线主要有两条:一种是采用吸收方式,即直接在镜片材质中加入吸收材料,如褐色素、黄色素等着色剂实现短波蓝光透射衰减,此类技术路线的典型代表为德国哈罗防蓝光护目眼镜,该方案通过在镜片中添加吸收剂,能过滤97%高能短波蓝光,但其局限性在于,该方案由于蓝光吸收过高将导致镜片发黄,严重破坏了色中性,这种产品只适用于专业电子防护,在日光条件下不能满足一般镜片增透镀膜的需要;另一种采用反射方式,即利用光学镀膜方式对不同波段的滤光进行处理[3-4]。这也是目前国内外眼镜企业及激光防护领域所采取的主流技术路线。基于反射方式的防护薄膜种类繁多,通过对特定波长光的滤除,可以起到防护效果[5-7]。例如,用于实现400~500 nm短波截止的滤光片称为干涉截止滤光片,通过高低折射率多层介质薄膜间的干涉原理,可对蓝光波段415 nm~455 nm实现高反射率,对475 nm~760 nm蓝光波段实现高透过率,进而对蓝光实现部分阻隔。现有方案的局限性在于,干涉截止滤光片的截止斜率宽窄仅通过改变高低折射率薄膜的膜对周期数来实现。其光谱目标值的取值方法多采取与工业标准匹配的分段取值方法。目标值的设计往往是线性且随机的,因而优化过程中不仅增加了运算量,还由于解的不唯一性引入额外误差。由于截止斜率的调控精度不高,容易造成445~475 nm的有益蓝光透射性能不达标的问题。
本文提出一种适用于GB/T38120—2019标准的斜率可控蓝光防护薄膜实用解决方案。通过改进薄膜设计过程中目标值的设定方法,利用非线性玻尔兹曼函数拟合长波通薄膜结构,获得非线性目标值;并进一步优化膜系中的膜层厚度使其满足光谱设计要求;改进的非线性目标值结合计算机辅助优化实现了厚度最优解的选取,得到了斜率可高精度调控的全介质薄膜结构。研究发现利用改进的目标值优化膜系结构,光谱通带波段具有明显的波纹减缓振荡效果。采用电子束蒸发离子束辅助沉积方法,成功制备了具有蓝光防护功能的光学多层薄膜。该目标值设定方法能够更好地结合薄膜优化算法简化优化过程,提高最优解选取精度,更好地满足蓝光防护薄膜膜系设计的需要。
由于白光LED的辐射光谱强度主要集中于400~500 nm范围,峰值辐射强度局限在415~455 nm附近。因此,445 nm以下波长的光会对人眼造成伤害,应该适当滤除。根据GB/T 38120—2019《蓝光防护膜的光健康与光安全应用技术要求》,蓝光防护膜的光透射比应满足表1中要求。另一方面,昼夜节律响应函数的峰值在464 nm处,主要位于446~477 nm的蓝光波段。同时根据国家轻工行业标准QB/T 2506—2017 《眼镜镜片光学树脂镜片》可知,445 nm以上波长的光,对人体节律具有有益调节功能,因此保证445~475 nm的高透过率至关重要,这也就引入了透射波段精确调控的问题。
表1 GB/T38120—2019蓝光防护膜的光透射比要求Tab. 1 Light transmittance requirement of blue-blocking protective thin films according to GB/T38120—2019
其中光透射比计算公式为:
式中λ为波长,单位为纳米(nm);E0(λ)为光源初始光谱能量分布;E1(λ)为光路中加入镜片后的光谱能量分布。为方便统一计算,用Es(λ)太阳光谱功率分布代替E0(λ)。为了对上述有害蓝光滤除及有益蓝光增强进行准确评估,利用两个参数精确评价蓝光护目镜片及薄膜的光谱性能:蓝光透射比τs,b可以量化评估由于暴露于蓝光辐射而可能对视网膜造成的损害程度;445~475 nm光透射比τc可用于衡量蓝光阻隔对节律的影响;τs,b与τc计算公式如下:
式(2)中,蓝光危害的波段主要集中于400~500 nm,峰值波长为437 nm。蓝光危害加权函数B(λ)数据来自GB/T 20145—2006。式(3)中,波段445~475 nm的选取可以用于量化镜片抑制蓝光辐射对人体昼夜节律的影响能力。镜片对昼夜节律的影响由τc表示。τs,b的取值范围为0~1,τs,b=0表示镜片完全阻挡蓝光,镜片完全阻挡蓝光对昼夜节律影响较大不利于褪黑素分泌。τs,b=1表示镜片没有蓝光阻隔能力,不会影响正常昼夜节律。
对于在500~800 nm波段具有很高的透过率,在相邻的400~450 nm波段上急速过渡到透射被完全截止状态,这种类型的光谱设计可以通过干涉截止滤光片实现[8]。常用高、低折射率相互交替的透明介质材料TiO2/SiO2来设计膜系,膜系初始结构可表示为(0.5HL0.5H)^s。s为对称膜系的周期;0.5H表示光学厚度为1/8波长的高折射率TiO2材料;L表示光学厚度为1/4波长的低折射率SiO2材料。TiO2/SiO2膜对组合的多层膜系在透射带内利用等效折射率E替换成单层膜,膜系反射率在两个数值之间振荡,分别是:
基片反射率:
上包络反射率:
目前,最普遍的消除波纹的方法为改变周期内的膜层厚度,以实现薄膜和基板以及薄膜和空气之间的相位匹配。结合计算机辅助算法使等效折射率变为更接近的预期值。光学薄膜的优化问题本质是通过构建初始膜系的响应特性R(λ)与所期望的理想特性RD(λ)之间的差别构成评价函数,来评价设计性能的优劣。优化方法主要有:最小二乘法、准牛顿法、共轭梯度法、单纯形法、黄金分割法以及聚类全局优化算法等[9]。对于自动优化设计来说,初始结构的选取对于最终优化结果的优劣尤为重要。初始结构的选取一般要求有经验的人员从平时的积累或者解析分析中得出。例如,滤光片的目标反射率可以取下列特征的理想化参数:
构造不同类型的评价函数,如计算平方和型评价函数最小值,进而得到理想的设计结果。
值得注意的是,目标值的线性理想参数在截止带与通带过渡区域取值不好,尖点与未能重合区域会在评价函数中占一定的比例,给评价函数收敛带来难度与额外的计算量。这就影响了滤光片在通带的光谱性能,导致通带内透射波纹不能很好消除。另外,该方法不能较好地实现斜率调控的目的,收敛过程随机,对于高精度控制斜率及计算机寻找特殊曲线类型的计算结果会令人失望。
本文采用与薄膜类型一致的函数关系获取非线性目标值。图1(彩图见期刊电子版)是常规光学薄膜设计流程图,S1,S2,S3为长波通薄膜结构不同膜对周期数S的仿真透射光谱曲线。可见,曲线斜率随着周期数的增加而变陡峭。以往的设计中,当薄膜材料固定时,周期数S为调节长波通滤光片斜率的主要参数。由图1可知,周期数增多虽然可以使得截止斜率更加陡峭,然而并不能实现对斜率的精准调控。通过OriginPro软件拟合长波通薄膜结构的光谱,得到非线性玻尔兹曼方程。根据国标要求,通过调整非线性玻尔兹曼方程的4个参数,能够实现对斜率及截止位点的参数化调制。
图1 非线性目标值的光学薄膜设计流程图Fig. 1 Flowchart of optical thin film design based on a nonlinear target
利用上述拟合方法得到的函数关系还不能满足蓝光防护的光谱设计要求,需要对函数方程参数进行调整。表2描述了不同参数条件下,玻尔兹曼函数曲线的分布规律。该函数利用参数a1表示初始值;a2表示最终值;(a2−a1)/dx表示曲线斜率;中点位置横坐标用X0表示。通过调节以上参数,使得函数曲线满足国标要求,即可找到满足条件的光谱目标值。A1-D1描述了中心坐标X0在400~460 nm每隔20 nm取值所得的函数曲线,用于分析蓝光衰减波长的最佳位点。A2-D2曲线为调整斜率的曲线,用于分析光谱斜率与蓝光防护性能之间的影响。
表2 玻尔兹曼目标值参数取值Tab. 2 Values of Boltzmann target parameter
将得到的8组曲线A1~D1,A2~D2带入式(2),式(3),得到蓝光阻隔参数τs,b以及昼夜节律影响参数τc,如图2(彩图见期刊电子版)所示。由于滤光片性能需要在实现短波蓝光衰减的同时,保证445~475 nm有益蓝光的透射功能,所以对于430~450 nm波段的斜率边缘陡峭程度的控制就显得尤为关键。图2中蓝色区域表示曲线能够同时兼顾蓝光阻隔及节律调节功能的取值范围。因此,利用B1曲线类型作为目标值用于多层薄膜设计,可以实现符合要求的深截止蓝光防护薄膜。
图2 蓝光阻隔参数τs,b与昼夜节律影响参数τc计算结果Fig. 2 Calculation results of blue hazard factors τs,b andcircadian rhythm factors τc
分别利用玻尔兹曼非线性目标值B1和线性目标值优化得到膜系结构,仿真透射光谱结果如图3(彩图见期刊电子版)所示。以K9为基板(ns=1.52),选用(H/2 L H/2)^S长波通薄膜初始结构。薄膜材料选用较高折射率材料TiO2(nTiO2= 2.25@510 nm)与SiO2(nSiO2= 1.45@510 nm)作为膜对组合,可以在较少的层数下达到光谱要求[10]。SiO2膜层结构为无定型态,在工作波长上色散较小,消光系数低,吸收少,是一种理想的低折射率材料[11-12]。目标值优化过程通过设定与长波通“S”线形一致的光谱目标(Target),利用玻尔兹曼函数得到满足蓝光光谱特性要求的目标函数B1,结合局部优化算法对膜层厚度进行精细修正来完成设计。膜系设计的初始结构保持了相同的14层薄膜结构,带入Macleod中优化,用以说明相同层数和膜系结构条件下,线性及非线性目标值对通带内光谱性能的影响。Simplex Parameters中的迭代次数等相关参数保持一致。目标值数据波长间隔选为1 nm,范围为380~780 nm。基于线性目标值优化得到的膜系Z1为:
膜层总厚度为702.24 nm,总层数为14层。数字代表物理厚度。同理,基于线性目标值优化得到的膜系Z2为:
膜层总厚度为697.37 nm,总层数为14层。数字代表物理厚度。
图3(a)、3(c)分别为基于非线性和线性目标值优化得到的薄膜400~800 nm透射光谱。图3(b)、3(d)为各种透射光谱500~800 nm波段的局部放大图。可见,“S”线形的函数目标值优化得到的透射振荡区间范围明显小于直线段线形目标值优化结果,具有更平坦的光谱性能。为了更好地分析基于线性及非线性目标值优化得到的多层膜通带透射曲线的增透及振荡程度,仿真得到不同层数的蓝光防护膜的导纳终点在500到800 nm处的轨迹图,如图3(e)~3(f)所示。根据反射公式,n0为空气的折射率,Y为膜层的等效导纳。若使薄膜在500~800 nm内高透,由R=1−T可知,n0=Y时可使反射为零。图3(e)中蓝色实线表示的导纳轨迹在500~800 nm内,位于Air(1,0)附近。而图3(f)中,500~800 nm内的导纳轨迹偏离(1,0)程度较高,因而透射光谱具有更高的振荡效果。对比分析上述线性及非线性目标值优化结果可看出,改进的玻尔兹曼目标值优化得到的蓝光防护薄膜结构具有明显的斜率控制和通带波纹压缩特点。分析仿真光谱可知,改进的目标值取值方法具有明显的斜率控制和通带波纹压缩特点。由500~800 nm导纳轨迹可知,通带波纹压缩是由于平滑的目标值曲线,特别是拐点处的目标值与光谱值差异性减少导致的。
图3 基于非线型((a),(b),(e))及线型((c),(d),(f))目标值设计蓝光防护薄膜的仿真透过率和导纳轨迹曲线Fig. 3 Theoretical design of transmittance curve of blue light protective film based on nonlinear ((a), (b), (e)) and linear ((c),(d), (f)) targets
实验采用ZZS850型电子束蒸发离子束辅助沉积系统分别制备蓝光防护单层薄膜和多层薄膜,离子源采用中国科学院空间中心K12考夫曼宽束离子源进行离子辅助沉积[13]。电子枪蒸镀过程中, TiO2容易失氧,会形成不同的结晶态。为了得到完整的TiO2膜,以Ti3O5颗粒材料作为蒸发的源材料最好,因为它可以产生相当恒定的折射率和更好的耐磨性。单层TiO2薄膜沉积参数:离子源偏压为300 V;离子束流为120 mA;沉积速率为0.3 nm/s;采用石英晶振控制,基片温度为65 ℃;Ar气流量为2 cm3/min,O2流量为10 cm3/min;镀膜基片选用厚度为2 mm、直径为20 mm的K9光学玻璃,镀膜前基片采用酒精和乙醚混合液清洗;并采用压强控制仪控制真空内的气体压力,控制真空度为1.5×10−2Pa。单层SiO2薄膜采用环形SiO2作为蒸发材料。沉积速率为1 nm/s,其它沉积参数完全一样。
薄膜软件Opti RE用于拟合测量的透射光谱以获得SiO2层的折射率。确定厚度及光学常数范围时,由于SiO2薄膜的消光系数在可见光波长范围内低于1.0×10−4,对透射的影响可忽略不计。TiO2薄膜短波弱吸收使用柯西定律拟合折射率与波长的关系。拟合得到的单层膜光学常数如表3所示。
表3 TiO2和SiO2单层膜的光学常数Tab. 3 Optical constants of the TiO2 and SiO2 thin films
采用电子束蒸发离子束辅助沉积技术,根据3.3节,膜系Z1结构选择TiO2和SiO2分别作为高低折射率薄膜材料,得到14层蓝光防护薄膜。对制备的单面蓝光防护光学薄膜,采用Lambda 750分光光度计,对可见光400~800 nm内的透射曲线进行测试,实验结果如图4所示。
图4 蓝光防护薄膜的400~800 nm透射曲线Fig. 4 Transmission of the blue light protective thin film at 400~800 nm
由图4可知,蓝光防护薄膜的测试曲线与理论设计曲线基本吻合。薄膜的光谱特性在385~415 nm的平均透过率小于15%,实现高能短波深截止;在415~445 nm的平均透过率小于50%,实现高能短波适当滤除;在445~475 nm的平均透过率大于80%,保障了有益蓝光透射对节律调节的需要;在475~800 nm平均透过率大于95%,满足了剩余可见光波段高透过率需要。
本文根据蓝光防护镜片及薄膜工业标准中500~800 nm波段高透过、400~450 nm波段蓝光深截止的要求,提出了利用非线性玻尔兹曼函数拟合长波通薄膜结构的方法。结果表明,通过改进的非线性光谱目标值优化膜系结构,可以实现更好的光谱性能。对比线性和非线性仿真光谱分析结果可知,改进的目标值取值方法具有明显的斜率控制和通带波纹压缩特点。对比线性和非线性500~800 nm导纳轨迹可知,通带波纹压缩是由于平滑的非线性目标值曲线,特别是拐点处的目标值与光谱值差异性减少导致的。最后,采用电子束蒸发离子束辅助沉积镀膜方法,在K9基片上成功制备了具有短波蓝光深截止及剩余可见光波段高透过率特性的蓝光防护功能薄膜。所制备的蓝光防护薄膜满足GB/T38120—2019要求。