一种基于光子轨道角动量的高精度波长计的制作方法
【专利摘要】本实用新型公开了一种基于光子轨道角动量的高精度波长计，通过该波长计的特殊光路结构可以将偏振和轨道角动量纠缠的量子态分开，从而实现携带两个极性相反的轨道角动量光子自我干涉，得到对波长很敏感的干涉花瓣图像，再对干涉花瓣图像的变化量进行分析，从而确定波长的变化量。上述方案鲁棒性强，容易集成，而且成本非常的低，同时，测量精度与分辨率也较高。
【专利说明】
一种基于光子轨道角动量的高精度波长计
技术领域
[0001] 本实用新型涉及光量子信息技术、光物理技术领域和计量测试技术领域，尤其涉 及一种基于光子轨道角动量的高精度波长计。
【背景技术】
[0002] 目前，波长计大多数通过经典光学的干涉或者衍射的原理制作的，比如德国 highfinesse公司生产的highfinesse波长计，是目前全球范围内精度最高的波长计，大陆 包括清华大学等高校在内的多个科研院所均采用这种波长计。其 WS6型号测量的分辨率达 到50Mhz，绝对精度200Mhz，光学分辨率(FWHM)接近4000。
[0003] 另一种比较普遍的是SHR-激光波长计，利用反射梯光栅将测量的光谱扩展到整个 波段，整个波长范围是190-1100nm，其近红外波长计SHR-IR波长分辨率在0.02nm，仪器分辨 率(FWHM)接近4000。
[0004] 然而，上述波长计大多采用动态干涉仪和频谱分析的方法。典型的有迈克尔逊干 涉仪、索菲干涉仪等，动态变化的干涉臂容易受外界影响（如电压等），且容易引进振动噪 声，快速傅里叶变换会丢失测量的精度，从而导致分辨率的下降，而科学研究或者技术加工 往往需要一个波长范围窄，精度高的波长计。
[0005] 鉴于此，有必要研究一种可适用于科学研究或者工厂加工的波长相对较窄、精度 相对高且成本较低的波长计。 【实用新型内容】
[0006] 本实用新型的目的是提供一种基于光子轨道角动量的高精度波长计，其测量精度 与分辨率较高，且成本较低。
[0007] 本实用新型的目的是通过以下技术方案实现的：
[0008] 一种基于光子轨道角动量的高精度波长计，包括:激光器、两个半波片HWP、三个偏 振分束器PBS、四个全反镜、两个聚焦透镜、四分之一波片QWP、涡旋光场产生片VPP、设置在 温控系统内的非线性光学晶体KTP、CDD及图像分析器;其中：
[0009] 所述两个半波片HWP分别记为HWP1与HWP2,三个偏振分束器PBS分别记为PBS1、 PBS2与PBS3,四个全反镜分别记为R1、R2、R3与R4,两个聚焦透镜分别记为L1与L2;
[0010] 所述激光器出射的相干光依次经过HWP、QWP、PBS1后射入全反镜R1中；
[0011] 所述全反镜R1将光束反射至PBS2后，一部分光经过PBS2透射后依次经过全反镜 R4、VPP、全反镜R3与全反镜R2后回到TOS2中，再透射至聚焦透镜L1;另一部分光经过PBS2反 射后依次经过全反镜R2、全反镜R3、VPP与全反镜R4后回到PBS2中，再反射至聚焦透镜L1; [0012] 聚焦透镜L1聚焦后的光束再依次通过KTP、聚焦透镜L2、HWP2与TOS3后射入(XD，该 C⑶的输出端还连有图像分析器。
[0013] 进一步的，所述PBS1与PBS3上还分别设有吸光器件。
[0014] 由上述本实用新型提供的技术方案可以看出，通过该波长计的特殊光路结构可以 将偏振和轨道角动量纠缠的量子态分开，从而实现携带两个极性相反的轨道角动量光子自 我干涉，得到对波长很敏感的干涉花瓣图像，再对干涉花瓣图像的变化量进行分析，从而确 定波长的变化量。上述方案鲁棒性强，容易集成，而且成本非常的低，同时，测量精度与分辨 率也较高。
【附图说明】
[0015] 为了更清楚地说明本实用新型实施例的技术方案，下面将对实施例描述中所需要 使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本实用新型的一些实施 例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图 获得其他附图。
[0016] 图1为本实用新型实施例提供的一种基于光子轨道角动量的高精度波长计的光路 结构示意图；
[0017] 图2为本实用新型实施例提供的花瓣重心点与边缘点的分布情况示意图；
[0018] 图3为本实用新型实施例提供的图像重心旋转的示意图
[0019] 图4为本实用新型实施例提供的实际测量波长与理论波长的关系示意图。
【具体实施方式】
[0020] 下面结合本实用新型实施例中的附图，对本实用新型实施例中的技术方案进行清 楚、完整地描述，显然，所描述的实施例仅仅是本实用新型一部分实施例，而不是全部的实 施例。基于本实用新型的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获 得的所有其他实施例，都属于本实用新型的保护范围。
[0021] 本实用新型提出了一种基于光子轨道角动量(0AM)的高精度波长计，其具有精度 高，测量范围大、系统鲁棒性强等优势，相比同类产品，其设备简单，成本低，易操作，易集 成，在科学研究和高精密仪器研制方面具有很大的潜力。
[0022]本实用新型是基于光子轨道角动量的独特性质。光子通过Sagnac(萨格纳克)干涉 仪后会产生自旋角动量和轨道角动量纠缠的光子3
，其中|H >、|V>分别代表水平偏振和垂直偏振的光子态；|+1>与|-1>表示两个极性相反的轨道角动 量，轨道角动量值为1;根据态的叠加原理，I#为偏振和轨道角动量纠缠的纠缠态;e1(a)代表 两个分立态之间的相位差。实验中光子态经过KTP晶体产生相位延迟e lP，其中
? 4为无量纲的比例因子，代表相位差，与KTP晶体的有效长度L，A n(A，T) 有直接的关系。A n(A，T)为水平偏振光和垂直偏振光经过晶体后产生的折射率差，这个折 射率差和光子的波长A、环境温度T有关，一般情况下L保持不变，则KTP晶体的温度T和入射 光的波长A将影响0的变化，而0的变化将影响着光子态|供> 的变化，最终影响光路中干涉图 像的变化。理论和实验均表明，0每变化2JI，对应干涉图像的花瓣将改变|弧度，故波长的 变化将导致图像的变化，反之，通过图像的变化也就可以推得波长的变化。
[0023] 光子既是经典信息也是量子信息的理想载体，单个光子不仅可以携带自旋角动量 (与光波的圆偏振相关），还可以携带轨道角动量(与光波的螺旋相位相关）。携带轨道角动 量的光子由于具有无限维的Hilbert空间，从而在光学通信和传感方面表现出了巨大的优 势。
[0024] 携带不同1值的光束的空间轨迹，以及干涉图像均存在区别，而前文所提到的光子 态是一个偏振（即自旋角动量)和轨道角动量纠缠的纠缠态，无法直接进行干涉，必须创造 出两个极性相反的轨道角动量进行自我干涉。
[0025] 下文中将结合附图1来介绍基于光子轨道角动量的高精度波长计的光路结构。如 图1所示，其主要包括：
[0026]两个半波片HWP、三个偏振分束器PBS、四个全反镜、两个聚焦透镜、四分之一波片 QWP、涡旋光场产生片VPP、设置在温控系统内的非线性光学晶体KTP、⑶D (图像传感器)及图 像分析器；以上器件的工作波长均在780nm附近。实验中晶体的温度用实验室自制的温控系 统，温控误差2mk。
[0027]其中：
[0028] 所述两个半波片HWP分别记为HWP1与HWP2,三个偏振分束器PBS分别记为PBS1、 PBS2与PBS3,四个全反镜分别记为R1、R2、R3与R4,两个聚焦透镜分别记为L1与L2;
[0029] 所述激光器出射的相干光依次经过HWP、QWP、PBS1后射入全反镜R1中；
[0030] 所述全反镜R1将光束反射至PBS2后，一部分光经过PBS2透射后依次经过全反镜 R4、VPP、全反镜R3与全反镜R2后回到TOS2中，再透射至聚焦透镜L1;另一部分光经过PBS2反 射后依次经过全反镜R2、全反镜R3、VPP与全反镜R4后回到PBS2中，再反射至聚焦透镜L1; [0031] 聚焦透镜L1聚焦后的光束再依次通过KTP、聚焦透镜L2、HWP2与TOS3后在CCD中成 像，所述图像分析器从CCD中获取波长变动前后的两幅图像，并根据两幅图像的变化量计算 波长变化。
[0032]如前文所述相位因子0，与KTP晶体的有效长度L、光子的波长A、环境温度T有关;其 中的KTP晶体的有效长度L通常不变，则光子的波长A与环境温度T为变量;本实用新型实施 例中，将KTP晶体置于温控系统中，使得环境温度T恒定，从而使得相位因子0的变化仅与光 子的波长X有关。
[0033] 本实用新型实施例中，PBS1与TOS3上还分别设有吸光器件，用于吸收roSl反射的 光束，以及PBS3透射的光束，以避免这些光束漫反射回光路。
[0034]以上为本实施例所提供的基于光子轨道角动量的高精度波长计的主要结构，下面 针对其工作原理进行介绍。
[0035]本实用新型实施例中，所述激光器出射的相干光依次经过HWP与QWP后聚焦在 Sagnac干涉环的PBS1前，光束的电场g I#表示为：
[0036] Ei](f))-^~En\ll) + ^El |K)； (1)
[0037] 其中，馬j为水平偏振光| H>的态系数，为一个电场强度量，其平方代表的光子数量， 是光强强度量；；爲为垂直偏振光I V〉的态系数，为一个电场强度量，其平方代表的光子数 量，是光强强度量；由于实验条件和器件的原因，所述I与&并非严格相等;e%为水平偏振 光子态和垂直偏振光子态之间的相位差，i为虚数，巾p为相位差因子，0P的数值影响光子的 偏振态，一个周期中（〇-2jt)，<J)P = 0°时，光束为45°线偏振光;么时，光为圆偏振光;其余 情况下为椭圆偏振光。
[0038] 本实用新型实施例中，光束经TOS2透射直至再次返回PBS2的路径方向为逆时针方 向，光束经TOS2反射直至再次返回PBS2的路径方向为顺时针方向；在逆时针方向时，PSB2到 VPP的几何光程记为La;在顺时针方向时，PSB2到VPP的几何光程记为Lb;顺时针方向和逆时 针方向两个路径因散射造成的光子丢失后剩余的光子数比例分别记为n 2和C，严格来讲两 者是不相等的。
[0039] 光束经过PBS2顺时针方向与逆时针方向行径并返回至PBS2后，相位变为如下： [0040]对于水平偏振光| H>，光子态相位的变化情况为：
[0041 ] |~〉=斤綱叫/7%|//〉|+/〉； (：2)
[0042] 对于垂直偏振光| V〉，光子态相位的变化情况为：
[0043] \<pl) = eiKL,+iK"k,)LBJr*l')C 1El |K)|-/); (3>
[0044]按照量子力学态的叠加原理，则有：
[0045] k) = |^;/)+i^, > = j//)| +/)+u)fi\r>|-：/}; (4)
[0046]其中，卜1>与|_1>表示两个极性相反的轨道角动量，轨道角动量值为l;a = kvLA+ (kv+ki) LB_khLB_( kh+ki) La+ 4* p
kindex( index - h、 v、l)分别对应的表不水平、垂直、轨道角动量三种不同光子的传播波矢量Aindex(index = h、 v、l)分别对应的表示水平、垂直、轨道角动量三种不同光子，no为空气中的折射率；当PBS2 所属的Sagnac干涉仪固定后，a仅由巾p决凭
_为水平光I H>与垂直光| V〉的光强相 对比例，当La=Lb，得到a=(i)P，这样不同波长也就不会对测量产生负面影响，即图像变换仅 仅与通过晶体的0变化有关，这大大增加了系统的鲁棒性和准确性。
[0047]本实用新型实施例中，所述聚焦透镜L1聚焦后的光束再依次通过KTP、聚焦透镜 L2、HWP2与PBS3后在(XD中成像，上述过程光子态的变化情况如下：
[0048]聚焦透镜L1聚焦后的光束进入KTP后水平偏振光与垂直偏振光因双折射会产生相 位差0，则叠加的光子态变为：
[0050] 再通过聚焦透镜L2射入HWP2后，叠加的光子态重新分布为：
[0052]最终，经过TOS3将光束的水平或者垂直偏振分开，水平或垂直光路方向将发生轨 道角动量光的干涉，从而实现光子的轨道角动量自我干涉，并在CCD中成像，最终的光子态 为：
[0054] 从上述式（7)可以看出，通过图1所示的光路将偏振和轨道角动量纠缠的态分开， 从而实现了两个光子的轨道角动量自我干涉，并在CCD中形成干涉图像;所述的干涉图像中 的花瓣会随着波长的变化而变化，通过检测花瓣变化情况，可以实现波长变化量的测量。
[0055] 当波长固定时，产生的干涉花瓣图像，是静止不动的，当波长发生一个微小的变化 时(调谐了），干涉花瓣图像会发生一个小角度转动，因为光的速度是3.0*10 8m/s，所以只要 波长发生变化图像会迅速发生转变，通过这种转变可以反映用户调谐的波长量。
[0056] 本实用新型实施例中，由图像分析器从CCD中波长变动前后的两幅图像，并根据两 幅图像的变化量计算波长变化：
[0057] 通过分析波长变动前后的两幅图像中的花瓣变化情况，再根据确定的花瓣变化情 况来确定波长的变化量。这一过程可以通过常规方案来实现，下文中以矢量点乘方法或者 光强互相关方法为例进行介绍。
[0058] 1、矢量点乘方法的步骤如下：
[0059] 1)对于两幅图像均进行二值化与洞处理，再进行边缘检测与重心识别，获得每一 花瓣的重心与边缘的坐标数列。
[0060] 原始的图像Bmp格式的图像，图像几个花瓣的面积不一样，且因为CCD的表面存在 污点，致使出现了细微的干涉条纹，这将大大的影响了光强矩阵的分布。
[0061] 对其二值化并进行洞处理，取一定的阈值实现二值化，然后将因 CCD杂质等造成的 中心暗点填充。
[0062] 之后，对上述对初步处理之后的图像进行边缘检测和重心识别，得到重心以及边 缘的坐标数列。
[0063] 2)将重心与边缘进行分离，获得每一花瓣重心坐标。
[0064] 通过相关的算法对重心与边缘进行分离，通过图2可以发现花瓣重心点与边缘点 的分布情况，计算x坐标数值导数：
，将阶跃点取出，也就是重心点。
[0065] 3)根据花瓣重心坐标来计算两幅图像中对应花瓣重心的旋转角度，从而确定花瓣 变化情况。
[0066] 得到重心坐标后，就可以计算两幅图之间的夹角，其基本思路是利用两幅图中圆 心矢量的点乘。设时间点相邻的两幅图像的重心坐标分别为： xi^i Xt , V,. ^1+15 ^/Vl X- , , J
[0067] 1 ? \ XH 2 5 3；,'+2 -^,12
[x;+3',r/+3'
[0068] 如图3所示，为两幅图像(分别用实线与虚线表示)重心旋转的示意图。
[0069]通过下式可获得重心的旋转角度：
[0071] 其中，？,、弓是CCD接收图像后建立的坐标系；9i与02表示两幅图像中四个花瓣中 的相邻两个花瓣旋转角度，这样做可以减少误差。
[0072] 2、光强互相关方法。
[0073]光强互相关方法是通过选取两幅图像中与花瓣相关的一定区域来进行强度相关 性的计算，计算出的相关性数值即为花瓣变化情况。
[0074]强度的相关函数定义为：
[0076] 上式中f(x，y)和g(x，y)分别代表选定区域的强度矩阵。C是计算后的相关矩阵。
[0077] 该方案为常规方法，旋转圆心的选取对相关系数有较大的影响，为了寻找两幅图 的旋转中心，需要采取一个区域进行搜寻，这样做的优势是非常精确的找到了旋转中心。
[0078] 另一方面，本实用新型实施例中，还需要在使用前进行标定，再以此对测量到的数 值进行修正，并最终呈现在显示界面上。具体如下：
[0079] 1)通过实验对波长计进行标定，获得实验值与理论值的比例。
[0080]标定是得到测量区域内实验的斜率值I，这个数值和理论数值I'有个差距，这个差 距会造成波长的测量值偏小或者偏大，所以需要先得到这两数值的比例
这个比例 和算法的选取以及波长的工作频带有关系，不同的算法显示这个精度在〇.〇1°。
[0081] 2)再利用实验值与理论值的比例对获得的波长的变化量进行修正，获得最终的波 长的变化量。
[0082] 3)最后，通过计算机后续算法处理，将最终的波长的变化量呈现在显示界面上。
[0083] 本实用新型实施例中，可以将获得的波长的变化量乘以或者才除以
来进行 修正。
[0084] 另外，还基于本实用新型上述实施例的方案进行了测试实验。具体如下：
[0085] 通过前述方案，可以知道波长是更容易改变相位因子0的，0随着折射率之差变换， 其中A为入射光的波长，T为晶体的温度，k为压强光学系数。
[0087]保持温度和压强光学系数不变，对相位差0取波长的全微分：
[0089]其中么11 = 112(\，1'，1〇-1^(、1'，1〇为〇?晶体的折射率之差，112、％为101 3晶体的2,7 主轴折射率，不通过的光子经过的折射率不一样。当通过晶体光的波长发生变化时，将导致 花瓣的旋转，其旋转角度为：
[0091 ]下面介绍在此情况下的角度的理论变化。
[0092]对于KTP晶体，在温度为25°时，波长以Mi为单位时，ny的理论值为：
[0094]其中各个系数的值如下所示： 4 =2.09930 A, = 0.922683
[0095] 1 ' A. - 0.0467695 ,4,-0.0138408
[0096] 同样对于n7.有：
[0098]相关的系数为： =2.12725 Z?;, -1.18431
[0099]]尽=5.14852x10 2 ; B4 = 0.6603 B, =：9.68〇56xl〇 ;
[0100] 实验中选取的波段在794 ? 619nm-795 ? 619nm(间隔0 ? 05nm)，这种波段下的波长每 变化0 ? 05个nm，其理论结果约为7 ? 3度。
[0101]在图像分析阶段，采用前文所述的矢量点乘方法与光强互相关方法，其结果如表1 所示：
[0102]
[0103] 表1矢量点乘方法与光强互相关方法的处理结果
[0104] 用线性拟合函数进行拟合，得到拟合方程，具体为：
[0105] Linear model Poly:
[0106] f(x) =pl*x+p2
[0107] Coefficients(with 95 % confidence bounds):
[0108] 理论值：
[0109] pl=7.537(7.535,7.538)
[0110] p2 = 0.05137(0.0326,0.07014)
[0111] Goodness of fit:拟合系数情况
[0112] SSE:0.00666
[0113] R-square:1
[0114] Adjusted R-square:1
[0115] RMSE:0.01924
[0116] 实验值：
[0117] 光强互相关方法：
[0118] pl=7.331(7.241,7.421)
[0119] p2 = 0.4796(-0.5992,1.558)
[0120] Goodness of fit:
[0121] SSE:22
[0122] R-square:0.9994
[0123] Adjusted R-square:0.9994
[0124] RMSE: 1.105 [0125]矢量点乘方法：
[0126] pl=7.374(7.292,7.457)
[0127] p2 = 0.6104(-0.3783,1.599)
[0128] Goodness of fit:
[0129] SSE:18.48
[0130] R-square:0.9995
[0131] Adjusted R-square:0.9995
[0132] RMSE:1.013
[0133] 其中，pi代表拟合后曲线的斜率，P2代表拟合后方程与Y轴的交点坐标。
[0134] 拟合系数中：Goodness of fit为适合度;SSE为拟合误差，这个数值越小代表模型 越好;RMSE root mean square error为均方根误差;R square称为方程的确定系数，0~1 之间，越接近1，表明方程的变量对y的解释能力越强;Adjusted R-square统计学中使用的 术语，跟R-Square类似，越接近1数据预测越好;RMSE是方均根误差，多用于统计学中，这个 数越小越好。
[0135] 同时，实际测量波长与理论波长的关系如图4所示。
[0136] 通过以上测试实验可知，矢量点乘方法的数据表现得更加集中，其均方根误差更趋向于 1。采用A A=7.374degree/0.05歷。这个数值标定后，采用我们的算法，图像检测算法的绝对精度 可以到0 ? 01度，即 A A=7 ? 374degree/0 ? 05nm=0 ? 073degree/0 ? 5pm，数量级上，波长每变化0 ? 1 pm， 角度变化量级为0.01度，根据:心销：-f仏-'喔… 这样得到的分辨率为50Mhz(0. lpm@795nm)，仪器分辨率V A A = 7950000,远远超过现有的 波长计。
[0137] 本实用新型实施例提供的波长计，可应用于780附近的红外波段，实验选定了 lnm 的带宽进行了测量，可以测量出波长在0.1pm量级上的变化，得到的分辨率为50Mhz(0.1pm@ 795nm)，仪器分辨率A/AA = 7950000,远远超过传统的干涉衍射类别的波长计。因此，所提 供的测量波长原理是个高精度的波长计。
[0138] 以上所述，仅为本实用新型较佳的【具体实施方式】，但本实用新型的保护范围并不 局限于此，任何熟悉本技术领域的技术人员在本实用新型披露的技术范围内，可轻易想到 的变化或替换，都应涵盖在本实用新型的保护范围之内。因此，本实用新型的保护范围应该 以权利要求书的保护范围为准。
【主权项】
1. 一种基于光子轨道角动量的高精度波长计，其特征在于，包括:激光器、两个半波片 HWP、三个偏振分束器I3BS、四个全反镜、两个聚焦透镜、四分之一波片QWP、涡旋光场产生片 VPP、设置在温控系统内的非线性光学晶体KTP、⑶D及图像分析器;其中： 所述两个半波片HWP分别记为HWPl与HWP2,三个偏振分束器I3BS分别记为I3BSl、PBS2与 PBS3,四个全反镜分别记为R1、R2、R3与R4,两个聚焦透镜分别记为Ll与L2; 所述激光器出射的相干光依次经过HWP、QWP、PBS1后射入全反镜Rl中； 所述全反镜Rl将光束反射至PBS2后，一部分光经过PBS2透射后依次经过全反镜R4、 VPP、全反镜R3与全反镜R2后回到TOS2中，再透射至聚焦透镜Ll;另一部分光经过TOS2反射 后依次经过全反镜R2、全反镜R3、VPP与全反镜R4后回到PBS2中，再反射至聚焦透镜Ll; 聚焦透镜Ll聚焦后的光束再依次通过KTP、聚焦透镜L2、HWP2与TOS3后射入CCD，该CCD 的输出端还连有图像分析器。2. 根据权利要求1所述的一种基于光子轨道角动量的高精度波长计，其特征在于，所述 PBSl与PBS3上还分别设有吸光器件。
【文档编号】G01J9/02GK205506228SQ201620191192
【公开日】2016年8月24日
【申请日】2016年3月11日
【发明人】刘世隆, 史保森
【申请人】中国科学技术大学