近年来, 波长探测作为科学技术中的一个基本工具, 已被广泛应用于分析化学^{[1, 2]}、 生物传感^{[3, 4]}、 材料分析^[5]以及光学通信^{[6, 7, 8]}等众多领域, 光谱分辨率、 工作带宽和器件尺寸是衡量探测性能的三个重要指标。 传统的波长计是通过色散组件分离到对应探测器位置^{[9, 10]}或不同谐振腔选择对应谐振分量^{[11, 12, 13]}来进行波长标定, 直接从空间通道即可得到对应的频谱分量。 随着光学器件的高性能和集成化发展, 基于模式干涉效应的高性能光谱计浮出水面^{[14, 15, 16, 17, 18]}。 该模式干涉来源于同一波长的多个传输模式间耦合效应, 总强度图案是每个模式电磁场分布的矢量叠加和, 传播常数差会随着传输距离增加逐渐积累, 并产生以2π 为周期的强弱分布。 通过采用探测器收集波长依赖的干涉图案, 得到每个波长唯一对应的指纹特征, 最后将待测信号图案进行光谱重建实现光谱探测。 这种机制通过引入更多模式(增加基模与高阶模传播常数差)来解决积累可探测相位差的问题, 在提升探测性能的同时可以极大缩小器件的物理尺寸。 Redding等^[16]使用无序光子晶体中依赖于波长的散射干涉图案, 在半径为25 μm尺寸实现了0.75 nm的分辨率和25 nm的工作带宽, 探测指标不够理想。 Paudel等^[18]提出了基于片上集成马赫-曾德干涉仪的傅里叶变换光谱计, 在达到140 MHz光谱分辨率和12 nm工作带宽的同时, 需要繁琐的制备工艺。 因此, 在微纳器件中同时实现分辨率、 工作带宽和探测效率优异的波长探测始终是人们追求的目标。

本文设计并实现了基于锥形多模光纤的波长计, 拉制多模光纤形成的锥形区域与共焦显微系统和可见光相机相集成, 后者可收集得到强烈依赖波长变化的模式干涉图案。 在光信号波长精细扫描得到的包含丰富细节图像中, 内积相关性运算联合主成分分析法可以高效处理频谱通道和空间通道构成的庞大数据量。 在400~700 nm的工作带宽内, 计算粗略校准矩阵的内积相关性可以初步实现精度为1 nm的波长探测; 在此基础上, 采用主成分分析法可将分辨率提高到20 pm, 波长探测准确率达到96.7%。 该算法不仅具有较高的波长探测精度, 运算效率也比传统重建算法提升近50倍。 该器件尺寸为π×(20 μm)²×0.5 mm, 具有高性能、 便携式和低成本等优势, 为基于光纤传输系统的波长实时探测提供了可能途径。

如图3所示, 当输入波长为635 nm时, 在不同锥形区域(光纤直径分别为40, 50和60 μm)沿光纤轴线提取强度干涉图像, 明显观察到类似周期性分布。 与高品质因数的谐振腔不同, 该结构可对其工作带宽内的波长分量进行相位差积累, 并形成细节丰富的频谱特征, 包括传导谐振产生的尖锐峰值、 法布里-珀罗谐振产生的宽谱背景变化和法诺干涉产生的不规则形状等。 这种频谱特征与尺寸参数和信号波长密切相关, 可以视为单色光波长信号的独特空间编码。 随着直径的增大, 一方面参与干涉的模式数急剧增多(增长率约20个· μm^-1), 另一方面反射路径变长, 两者都会导致空间周期变短。 这里通过拍长Δ z=2π /Δ β 进行描述: 由于Δ β /β 远小于1, 可以用准直波导和线性增长的传播常数进行模型近似, 得到Δ z为43 μm。 在实验中, 根据显微物镜放大倍数和相机像素尺寸进行计算, 拍长为41 μm, 与估计值基本一致。 虽然拍长在微调信号波长值时变化不明显, 但局部强度分布仍然有所改变(闪烁现象), 说明充分利用图像细节可以得到更高分辨率, 势必引入更多空间信息带来计算冗余度。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 不同锥形区域即光纤直径(a) 40 μm, (b) 50 μm和(c) 60 μm下, 沿轴线提取的强度分布图Fig.3 The intensity distribution patterns are extracted along the axis of the fiber with different taper diameters (a) 40 μm, (b) 50 μm and (c) 60 μm

由于模式干涉图案具有强烈波长依赖性, 与整个工作带宽对应的校准矩阵可以表征该器件的内在特性。 内积相关性反映某个波长与其他波长产生的强度图像间相似程度: 归一化后值为0说明两张图像完全不相同; 值为1说明图像完全相同; 而值在0~1之间说明有一定相似性。 这从另一个角度衡量了器件的波长探测能力。 采用超连续谱光源+单色器组合的激光光源(带宽1 nm)获取粗略校准矩阵, 进行内积运算后可以粗略实现纳米级别的光谱分辨率。 这里, 对不同结构参数组合的锥形光纤进行测量计算: 如图4(a)所示, 光纤直径越大, 激发出的干涉模式越多, 对波长分辨能力就更强(对应半高全宽越小); 如图4(b)所示, 探测长度越长, 会包含进更多的空间通道作为采样点, 波长分辨能力也越强, 从实验上印证了理论模型的准确性。 在特定参数器件中, 扫描信号波长形成的模式干涉图像“ 闪烁” 现象说明空间与波长通道一一映射, 内积曲线应当具有较高信噪比。 如果选用窄带信号光源、 克服工艺不完美性和提升可见光相机性能(如降低像素尺寸、 提高量子效率等), 通过内积运算即可达到最优分辨率。 综合考虑波长计各项性能指标, 选择测量直径为40 μm和长度为500 μm的光纤区域。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 测试波长选用635 nm时, 内积曲线与锥形光纤(a)直径和(b)长度的关系Fig.4 The relationship between the inner product and (a) diameter and (b) length of the tapered fiber at the wavelength of 635 nm

采用窄带可调谐激光器634.80~639.40 nm进行步长为10 pm的波长扫描获得精细校准矩阵。 随机选取其中两个波长并在不同功率下收集强度图像, 经过主成分分析算法可以准确探测出分立的单色光信号, 如图5(a)所示, 探测结果与校准波长基本吻合。 光谱分辨率是衡量光谱器件辨识邻近波长能力的重要指标, 定义为重建单色光信号的半高全宽(FWHM)。 在此基础上, 选取两个相邻波长的强度图像并逐渐增加波长间隔, 图像叠加后与精细校样矩阵进行内积运算得到如图5(b), 根据瑞利判据得到分辨率为20 pm, 证明该器件具有皮米量级的分辨能力。 将一组待测单色光信号同时输入该器件和商用光谱仪(Princeton Instrument SP2500, 分辨率0.26 nm)中, 前者通过两步法(内积相关性运算和主成分分析算法)探测波长后, 与后者测量结果(校准波长)进行比较, 如图5(c)所示: 在实验测量过程中, 设置激光器632.0~638.0 nm范围、 以0.1 nm为步长进行波长扫描, 将波长计测量结果作为探测波长并与校准波长进行比较, 两者吻合并保持理想的线性关系(为了清楚起见, 62个波长值中只选择奇数序号波长值进行展示)。 对于波长结果的认定, 定义差值超过分辨率一半即10 pm即为探测误差, 由测量结果得到波长探测准确率为96.7%, 如图5(d)所示, 误差主要源于激光器中心波长漂移和器件工艺缺陷。

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 (a)两个分立单色光信号的探测; (b)波长计分辨率为20 pm; (c)探测波长与校准波长的关系曲线; (d)波长探测结果统计分析图Fig.5 (a) Detection of two separate monochromatic light signals; (b) spectral resolution of 20 pm; (c) plot of the detected wavelength versus calibrated wavelength; (d) statistical histogram of wavelength detection results

基于锥形多模光纤波长计的工作带宽由多模光纤和CCD探测范围共同决定。 受限于激光光源种类和参数, 实验中仅选择630~640 nm波长范围进行波长探测。 考虑到多模光纤传输窗口和硅材料探测曲线, 采用超连续谱光源与单色器组合对该范围进行步长为1 nm的粗略扫描, 每个波长值均得到单峰内积相关性曲线, 说明在该范围内模式干涉图案与波长一一映射, 工作带宽达到400~1 100 nm范围。 得益于内积相关性和主成分分析联合算法, 该锥形光纤波长计不必扫描感兴趣区域以外的精细波长, 减少了近50倍数据存储量和重建时间(=1 μm/20 nm), 探测效率(与波长探测时间成反比)较常规重建算法提高近50倍^{[15, 16, 17, 18]}。 与宽谱连续信号相比, 分立窄带信号仅有有限频谱通道, 重建时不易遗漏信号分量, 该器件更加适合波长探测的应用。 在实际应用中, 器件的机械性能、 温度特性和可重复性至关重要。 前两者可以通过固定封装和反馈控制进行克服, 而后者主要来自校准光路的漂移和散焦。 将1, 5和12 h后扫描得到的校准矩阵与初始矩阵进行内积运算, 得到如图6所示的内积带宽曲线。 较长时间后, 光路未经校准的半高全宽增加0.1 nm(变化率为13.8%), 而经校准仅增加0.01 nm(2.5%), 说明可重复性能够通过自动校准系统或将多模光纤与探测器阵列相互集成得到有效改善。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 内积带宽随时间变化规律Fig.6 Variation of inner product bandwidth with time

在基于锥形多模光纤的波长计中, 证实了内积相关性运算可以初步实现粗略波长探测, 而进一步采用主成分分析法可以实现20 pm的精细探测。 利用搭建的显微共焦成像系统, 在微纳尺度器件中不仅具有较高波长探测精度, 运算效率也较传统算法大大提升。

强度分布图像中的噪声主要包括散粒噪声和暗噪声两种成分。 对于散粒噪声, 假设探测信号在某波长处光子数为Np, 则在每个频谱通道中光子平均数为Np/N, 总散射噪声功率为两者乘积平方根, 这与光栅光谱仪的散射功率相同; 而对于暗噪声, 需要在所有空间通道进行积累, 对光谱重建作用与信号强弱、 探测深度等因素密切相关。 将基于锥形多模光纤波长计的性能向极限推进, 可行的途径有两个: 通过模式干涉图案去相关、 降低实验噪声和误差后处理等手段进一步优化探测算法; 改善可见光相机性能, 如降低像素尺寸和探测噪声、 提高量子探测效率等, 这将是下一步研究方向。 由于锥形多模光纤侧面泄露的信号能量仅占总入射能量的20%~30%, 在基于光纤连接和传输的通信系统中, 可通过剥皮、 拉锥和成像三个步骤进行波长实时探测。