CN116859300B - 基于金刚石nv色心的量子传感频率跟踪控制方法及系统 - Google Patents
基于金刚石nv色心的量子传感频率跟踪控制方法及系统 Download PDFInfo
- Publication number
- CN116859300B CN116859300B CN202311124042.1A CN202311124042A CN116859300B CN 116859300 B CN116859300 B CN 116859300B CN 202311124042 A CN202311124042 A CN 202311124042A CN 116859300 B CN116859300 B CN 116859300B
- Authority
- CN
- China
- Prior art keywords
- frequency
- magnetic field
- microwave source
- demodulation amplitude
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 229910003460 diamond Inorganic materials 0.000 title claims abstract description 25
- 239000010432 diamond Substances 0.000 title claims abstract description 25
- 238000006243 chemical reaction Methods 0.000 claims abstract description 32
- 238000005259 measurement Methods 0.000 claims abstract description 18
- 230000006698 induction Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 abstract description 5
- 238000010521 absorption reaction Methods 0.000 description 11
- 230000003595 spectral effect Effects 0.000 description 9
- 230000005283 ground state Effects 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000000387 optically detected magnetic resonance Methods 0.000 description 6
- 230000005281 excited state Effects 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 230000010354 integration Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 239000013078 crystal Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000003776 cleavage reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002189 fluorescence spectrum Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000007017 scission Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 125000004432 carbon atom Chemical group C* 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 125000004433 nitrogen atom Chemical group N* 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/032—Measuring direction or magnitude of magnetic fields or magnetic flux using magneto-optic devices, e.g. Faraday or Cotton-Mouton effect
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种基于金刚石NV色心的量子传感频率跟踪控制方法及系统,属于量子传感技术领域。本发明通过在金刚石NV色心传感头处引入辅助磁场测量装置,先根据磁场与NV色心磁共振频率的关系,利用磁场测量装置的测磁结果推算出频率偏差,完成微波源输出信号中心频率的粗调。随后根据解调幅值与频率偏差的转换系数计算解调输出对应的频率偏差,完成微波源输出信号中心频率的细调。通过粗调加细调实现微波源输出信号中心频率对NV色心磁共振频率的宽范围、精准跟踪,提高系统测量动态范围和精度。此外,本发明还设置了保护阈值,当磁场变化较小,频率偏差值处在线性约束范围内时,跳过粗调过程,直接进行频率细调,加快调节速度,提高系统测量带宽。
Description
技术领域
本发明属于量子传感技术领域,更具体地,涉及一种基于金刚石NV色心的量子传感频率跟踪控制方法及系统。
背景技术
量子传感技术是指利用量子系统、量子特性或量子现象来对某一物理量进行测量。现阶段,在量子传感领域,NV色心具有的高灵敏度的量子效应优势以及宽松的工作条件使得它在产业应用方面具有广阔的前景。随着量子科技在各个领域的飞速发展,量子传感也广泛应用于电力能源传输、医疗设备、生物检测、电磁成像等多个领域。
单个NV色心是一个简单的二能级体系,分为基态3A2与激发态3E,均为自旋三重态(包含ms=0、ms=+1、ms=-1三种自旋态)。由于自旋态ms=±1态被激发到激发态后存在非辐射跃迁过程,导致NV色心发出的荧光状态有明暗的区别。施加连续激光,会使得电子自旋态都趋向于ms=0态,此时荧光强度最强。施加微波后,当微波频率满足ms=0态和ms=±1态之间的能级频率差,即达到共振频率时,会让处于ms=0和ms=±1自旋态的电子布居度发生改变,进而使收集的荧光谱线强度发生变化,形成吸收峰。进一步地,由于磁场导致的塞曼分裂会导致自旋态ms=±1态的能级发生退兼并,ms=±1态共振频率会随着外磁场的出现而发生分裂,不考虑微弱的电场和核自旋影响,分裂后的共振频率(分别为ms=0到ms=+1和ms=0到ms=-1)满足下式:
其中,为磁旋比,/>,/>为外加待测磁场与NV轴向的夹角,为外加待测磁场的强度。通常利用微波源扫频后可以根据荧光谱线计算共振频率间的差值推算磁场大小,但该方法每次都需要一次完整的微波频率扫描过程,这个时间通常在秒级以上。因此在待测场快速变化时,此方法难以适用。若固定微波频率,并保证磁场变化时某一吸收峰对应的微波中心频率始终在共振频率附近的一段频率-荧光强度的线性范围内,此时外磁场变化会导致磁共振频率变化,偏离微波源中心频率,荧光强度随之发生变化。利用该范围内荧光强度变化和频率偏移的近似线性关系,结合频率偏移和外磁场强度的关系,即可通过NV色心的荧光强度来表征外磁场强度。这种方法将磁场与荧光强度对应,对信号的响应速度较快,但会受到大量噪声干扰而产生荧光强度波动,造成测量精度较差。此外,由于荧光强度变化的线性区间范围有限,难以实现大量程范围的测量。
为了减轻噪声的干扰,可以利用频率调制方法将荧光信号调制到高频后做锁相解调,利用该方式得到的解调输出为原始荧光信号的一阶微分谱线,这种调制解调方案减弱了低频噪声干扰,并且微分谱线先比于原始信号有了更宽的线性区。但是该线性区仍受限于吸收峰有限的本征半高宽。为了进一步扩大信号的测量范围,可以采用频率跟踪的方式使微波源中心频率实时跟踪磁共振频率的变化,进而使得线性区会随信号变化而移动。通过把锁相放大器解调输出信号当成频率跟踪控制模块的反馈输入信号,频率跟踪控制模块把共振频率附近的解调曲线斜率当作解调输出-频率偏差转换系数,输出微波源中心频率和NV色心共振频率之间的频率偏差信号,并以此调整微波源输出信号的中心频率值。这样可以使每一次磁场测量过程都处于频率-解调幅值线性区约束范围内,具有较高的测量精度。但该方式存在的问题在于系统处于在待测磁场快速、大幅度改变情况下,前后两次测量值之间的差值会超出线性区约束范围而无法获取到有效的调频依据,出现跟踪丢失的问题。
基于这种背景,为避免频率跟踪控制方法超出频率-解调幅值线性约束范围后出现跟踪丢失的问题,保证频率跟踪控制方法在磁场快速、大幅度改变场景下的精准测量,一种用于大动态范围、高精度量子传感的共振频率跟踪控制方法对磁场或其他参量的量子传感系统及其应用具有重要的意义。
发明内容
针对现有技术的缺陷和改进需求,本发明提供了一种基于金刚石NV色心的量子传感频率跟踪控制方法及系统,旨在解决传统频率跟踪控制方法超出线性约束范围后出现跟踪丢失的问题,保证频率跟踪控制方法在磁场快速、大幅度的幅值突变下的精准测量。
为实现上述目的,第一方面,本发明提供了一种基于金刚石NV色心的量子传感频率跟踪控制方法,包括以下步骤:
S1,根据计算得到频率偏差ε1,其中,/>为磁旋比,/>为外加待测磁场与NV轴向的夹角,/>,/>为当前时刻通过辅助磁场测量装置获取的外加待测磁场强度,/>的初始值为零;将微波源的频率调整为f 1,其中,f 1=f 0+ε1,f 0为上一轮循环结束时微波源的频率,f 0的初始值为微波源起始频率;
S2,根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2,其中,f 2=f 1+δ2,所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;
S3,计算频率偏差ε=f 2-f 0,并根据计算得到磁场强度的变化量/>,将磁场强度的变化量/>与/>求和后赋值给/>,再返回S1继续下一轮频率跟踪控制。
进一步地,将S2替换为S2',所述S2'包括:
S21',若f 1频率下的解调幅值V大于误差阈值W,则根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2并执行S3,其中,f 2=f 1+δ2,所述转换系数K通过获取f0前后的多个频率点对应的解调幅值拟合得到;否则,执行S22';
S22',对当前V进行累加,获得累加值Vi,同时计数器值n执行自加1操作;
S23',若累加值Vi大于误差阈值W或者计数器值n大于最大累计次数,根据δ2'=Vi/(n·K)计算频率偏差δ2',再将微波源的频率调整为f 2并执行S3,其中,f 2=f 1+δ2';否则,返回S21'。
为实现上述目的,第二方面,本发明提供了另一种基于金刚石NV色心的量子传感频率跟踪控制方法,包括以下步骤:
S0,若f 0频率下的解调幅值大于保护阈值,则执行S1;否则,执行S3;其中,f 0为上一轮循环结束时微波源的频率,f 0的初始值为微波源起始频率;
S1,根据计算得到频率偏差ε1,其中,/>为磁旋比,/>为外加待测磁场与NV轴向的夹角,/>,/>为当前时刻通过辅助磁场测量装置获取的外加待测磁场强度,/>的初始值为零;将微波源的频率调整为f 1,其中,f 1=f 0+ε1;
S2,根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2并执行S4,其中,f 2=f 1+δ2,所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;
S3,根据解调幅值与频率偏差的转换系数K计算f 0频率下的解调幅值对应的频率偏差δ2',再将微波源的频率调整为f 2并执行S4,其中,f 2=f 0+δ2',所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;
S4,计算频率偏差ε=f 2-f 0,并根据计算得到磁场强度的变化量/>,将磁场强度的变化量/>与/>求和后赋值给/>,再返回S0继续下一轮频率跟踪控制。
进一步地,将S2替换为S2',所述S2'包括:
S21',若f 1频率下的解调幅值V大于误差阈值W,则根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2并执行S4,其中,f 2=f 1+δ2,所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;否则,执行S22';
S22',对当前V进行累加,获得累加值Vi,同时计数器值n执行自加1操作;
S23',若累加值Vi大于误差阈值W或者计数器值n大于最大累计次数,根据δ2''=Vi/(n·K)计算频率偏差δ2'',再将微波源的频率调整为f 2并执行S4,其中,f 2=f 1+δ2'';否则,返回S21'。
进一步地,将S3替换为S3',所述S3'包括:
S31',若f 0频率下的解调幅值V大于误差阈值W,则根据解调幅值与频率偏差的转换系数K计算f 0频率下的解调幅值对应的频率偏差δ2',再将微波源的频率调整为f 2并执行S4,其中,f 2=f 0+δ2',所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;否则,执行S32';
S32',对当前V进行累加,获得累加值Vi,同时计数器值n执行自加1操作;
S33',若累加值Vi大于误差阈值W或者计数器值n大于最大累计次数,根据δ2'''=Vi/(n·K)计算频率偏差δ2''',再将微波源的频率调整为f 2并执行S4,其中,f 2=f 0+δ2''';否则,返回S31'。
为实现上述目的,第三方面,本发明提供了一种基于金刚石NV色心的量子传感频率跟踪控制系统,包括:微波源、锁相放大器、辅助磁场测量装置、频率跟踪控制模块、NV色心传感头;
所述辅助磁场测量装置安装于所述NV色心传感头处,用于测量外加待测磁场的强度并输出至所述频率跟踪控制模块;所述频率跟踪控制模块还用于接收所述锁相放大器输出的解调幅值,以实现第一方面或第二方面所述的基于金刚石NV色心的量子传感频率跟踪控制方法。
进一步地,所述辅助磁场测量装置包括感应线圈和积分器,或者感应线圈、模数转换器和数字积分器。
总体而言,通过本发明所构思的以上技术方案,能够取得以下有益效果:
(1)本发明利用荧光谱线的解调幅值作为频率跟踪控制模块的输入,充分利用了解调幅值-频率曲线中的线性区约束范围,利用该范围内频率偏差和荧光解调幅值的线性关系来使微波源的频率实时跟踪金刚石NV色心的共振频率。通过在金刚石NV色心传感头处引入辅助磁场测量装置,先根据磁场与NV色心磁共振频率的关系,利用磁场测量装置的测磁结果推算出频率偏差,完成微波源输出信号中心频率的粗调。随后根据解调幅值与频率偏差的转换系数计算解调输出对应的频率偏差,完成微波源输出信号中心频率的细调。通过粗调加细调实现微波源输出信号中心频率对NV色心磁共振频率的宽范围、精准跟踪,提高系统测量动态范围和精度。
(2)本发明中引入线性区间保护阈值,以保护阈值为非线性区的判断依据,当超过保护阈值才会读取辅助磁场测量装置输出,避免每个控制周期都需先读取辅助磁场测量装置进行粗调,从而优化系统资源,加快调节速度,同时也可以解决解调幅值-频率曲线在线性区约束范围外出现跟踪丢失的问题,适用于待测信号频率带宽、幅值都较大的情况。
(3)为了避免测量干扰而引起的微波源改变中心频率的问题,通过设置误差阈值的方式,将误差阈值以内的频率偏差进行累加,从而可以减弱部分噪声误差导致的频率跟踪控制误动问题。
附图说明
图1和图2分别为本发明实施例所涉及到的金刚石NV色心晶体结构及其能级结构。
图3为本发明实施例提供的扫频ODMR谱线。
图4为本发明实施例提供的ODMR定频方法示意图。
图5为本发明实施例提供的量子传感系统简化模型。
图6为本发明实施例提供的频率跟踪控制方法示意图。
图7为本发明实施例提供的一种基于金刚石NV色心的量子传感频率跟踪控制方法流程图之一。
图8为本发明实施例提供的一种基于金刚石NV色心的量子传感频率跟踪控制方法流程图之二。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
在本发明中,本发明及附图中的术语“第一”、“第二”等(如果存在)是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。
参见图1和图2,本发明实施例基于图1和图2中所示的金刚石NV色心晶体结构与能级结构。本实施例以金刚石NV色心所构成的量子系统为例,基于其能级结构和量子效应来进行量子传感。图1所示晶体结构中,一个氮原子(Nitrogen)取代其中一个碳原子,并且捕获其临近位的一个空穴,形成了一个基本的NV色心结构。单个NV色心是一个简单的二能级体系,如图2所示,分为基态3A2与激发态3E,均为自旋三重态。在基态与激发态之间存在两个亚稳单态1A1、1E。通过激光泵浦可以实现能级跃迁从基态跃迁到激发态,自旋0态会通过辐射跃迁回到基态,而自旋±1态通过非辐射跃迁(不发光)回到基态,此时自旋态会变到0态,这导致处于基态自旋0态的电子发出的荧光更强,因此可以通过对荧光光谱的探测判断NV色心中的自旋态,进而实现磁场测量。
参见图3,本发明实施例在初次扫频后的ODMR谱线如图。由于磁场会导致自旋态±1发生分离,因此从自旋0态到自旋+1态和从自旋0态到自旋-1态的共振频率不同,从而可以利用ODMR谱中劈裂吸收峰中的共振频率差反应外部磁场。在图3中,NV色心轴向与磁场方向平行,为1,记/>为磁旋比。/>与/>为外加磁场后导致的劈裂吸收峰,与/>的差值为/>,/>为NV色心的感应磁场值。
参见图4,利用原始谱线进行固定微波频率、测量荧光强度推算频率偏差的方法来测量出NV色心的共振频率,从而算出磁场强度和电流幅度。
参见图5,本发明实施例可应用在包括但不限于图5中所示量子传感系统中,以本实施例金刚石NV色心的量子传感系统为例,所述量子传感系统包括光学系统、微波系统、标准生磁系统、NV色心传感头(辅助磁场测量装置可为感应线圈+积分器)、读出系统、频率跟踪控制系统、上位机。所述光学系统包括:激光源、二向色镜、物镜、透镜、滤波片等;所述微波系统包括微波源、功率放大器、环形器、天线等;所述读出系统包括光电探测器、锁相放大器等;所述标准生磁系统包括标准电流源和生磁线圈;微波源、锁相放大器、辅助磁场测量装置、频率跟踪控制模块、NV色心传感头组成了共振频率跟踪控制系统。
其中,微波源用于提供微波信号、进行微波频率扫描以确定起始跟踪频率、对微波信号进行频率调制(对荧光信号进行调制),同时作为频率跟踪控制模块的控制对象,通过改变输出微波的中心频率以执行频率跟踪控制。锁相放大器用于对信号进行解调,输出的解调信号发送至频率跟踪控制模块,作为细调微波源中心频率的依据。辅助磁场测量装置用于同步监测磁场,并将磁场监测值传送至频率跟踪控制模块,用以提供微波源输出信号中心频率偏离NV色心磁共振频率后的频率粗调依据。频率跟踪控制模块接收来自辅助磁场测量装置的输出进行频率跟踪控制,同时也接收锁相放大器的输出信号进行频率跟踪控制。频率跟踪控制模块首先根据来自辅助磁场测量装置的输出信号完成一次频率跟踪调节,然后再根据锁相放大器的解调输出值完成频率的进一步跟踪调节。
参见图6,在本发明实施例进行频率跟踪控制时利用的锁相放大器的解调输出谱线中,实线部分表示当前磁场下锁相放大器输出的解调谱线,虚线部分表示待测磁场发生变化后的锁相放大器解调谱线;当微波源中心频率等于NV色心磁共振频率时,锁相放大器解调输出为零。当磁场变化时,共振频率也会变化,此时锁相放大器解调谱线出现整体偏移,若不改变微波源输出信号的中心频率,此时解调输出为V。此时频率跟踪控制方法会利用解调输出-频率偏差转换系数K将解调输出V转换成频率偏差,同时以此调整微波源中心频率,使中心频率其跟踪磁共振频率。图6中线性约束范围可利用解调谱线的微分曲线依据系统欲实现的精度进行确定;保护阈值可定为线性约束范围端点幅值95%处;误差阈值依靠系统所需达到的灵敏度进行确定,本实施例选为端点幅值0.5%处。
参阅图7,本发明实施例提供了一种适用于上述系统的频率跟踪控制方法,通过粗调加细调实现微波源输出信号中心频率对NV色心磁共振频率的宽范围、精准跟踪,所述方法包括以下步骤:
步骤1,程序上电开始运行。
步骤2,初始化过程,激光源设定发射功率并开启,微波系统设定功率和扫频范围并将微波源FM调制模式设定为外参考模式,设置FM参数,清空控制系统计数器和寄存器等。
步骤3,若检测到荧光信号,则执行步骤4,否则执行步骤2。
步骤4,上位机控制外部标准电流源输出电流,并由线圈转换成静磁场施加到NV色心探头上,微波源开始扫频,读出ODMR谱中磁共振吸收峰的共振频率f a、f b,并计算△f=fb-f a。
步骤5,若A=|1-△f/(2•• B)|<校准阈值,则执行步骤7,否则执行步骤6。
步骤6,上位机校准失败预警,需要进行设备手动校准,执行步骤21。
步骤7,自校准完成,上位机指示灯亮,提升系统工作正常,进入频率跟踪控制部分。
步骤8,上位机设置远端上共振吸收峰f b(或下共振频率f a)作为微波源输出中心频率,即以该吸收峰作为频率跟踪控制方法的起始频率f 0。同时开启FM调制,扫频后计算解调幅值与频率偏差的转换系数K,其中解调幅值与频率偏差的转换系数K选择为共振频率点前后N点范围内的解调幅值与频率曲线的差值斜率。
步骤9,频率跟踪控制模块检测锁相放大器输出的解调幅值,若解调幅值V为0,则不改变中心频率,控制系统向上位机指示频率偏差为0,并继续检测锁相放大器输出的解调幅值直至解调幅值V不为0;否则执行步骤10。
步骤10,控制系统读出辅助磁场测量装置输出,此时输出满足:
其中,为当前时刻通过辅助磁场测量装置获取的外加待测磁场强度,/>是积分线圈的实时输出值,/>是积分线圈的面积。计算此时的频率偏差为:
,/>的初始值为零;
使中心频率f 1=f 0+ε1进行跟踪。
步骤11,若解调幅值V小于误差阈值W,执行步骤12,否则执行步骤16。
步骤12,锁相放大器解调幅值V累加得到累加值Vi,同时计数器的值n执行+1操作。
步骤13,若解调幅值V大于误差阈值W或者n大于最大累加次数Nmax(例如Nmax=5),执行步骤14,否则执行步骤11。
步骤14,计算频率偏差δ2’,计算公式如下:
δ2'=Vi/(n·K)
其中,n为计算器当前值,计算完成后计数器值n清零。
步骤15,使中心频率f2=f1+δ2'进行跟踪,执行步骤18。
步骤16,计算频率偏差,计算频率偏差δ2,计算公式如下:
δ2= V/ K
其中V为锁相放大器解调幅值,K为解调谱线零点处斜率,即线性区斜率。
步骤17,使中心频率f 2=f 1+δ2进行跟踪。
步骤18,计算总频率偏差ε=f 2-f 0,并根据计算得到磁场强度的变化量/>,将磁场强度的变化量/>与/>求和后赋值给/>。
步骤19,若仍检测到荧光信号,则执行步骤9,否则执行步骤20。
步骤20,结束跟踪调节过程。
参见图8,本发明实施例提供了一种适用于上述系统的频率跟踪控制方法,引入线性区间保护阈值,避免每个控制周期都需先读取辅助磁场测量装置进行粗调,从而优化系统资源,所述方法包括以下步骤:
步骤1,程序上电开始运行。
步骤2,初始化过程,激光源设定发射功率并开启,微波系统设定功率和扫频范围并将微波源FM调制模式设定为外参考模式,设置FM参数,清空控制系统计数器和寄存器等。
步骤3,若检测到荧光信号,则执行步骤4,否则执行步骤2。
步骤4,上位机控制外部标准电流源输出电流,并由线圈转换成静磁场施加到NV色心探头上,微波源开始扫频,读出ODMR谱中磁共振吸收峰的共振频率f a、f b,并计算△f=fb-f a。
步骤5,若A=|1-△f/(2•• B)|<校准阈值,则执行步骤7,否则执行步骤6。
步骤6,上位机校准失败预警,需要进行设备手动校准,执行步骤21。
步骤7,自校准完成,上位机指示灯亮,提升系统工作正常,进入频率跟踪控制部分。
步骤8,上位机设置远端上共振吸收峰f b(或下共振频率f a)作为微波源输出中心频率,即以该吸收峰作为频率跟踪控制方法的起始频率f 0。同时开启FM调制,扫频后计算解调幅值与频率偏差的转换系数K,其中解调幅值与频率偏差的转换系数K选择为共振频率点前后N点范围内的解调幅值与频率曲线的差值斜率。
步骤9,频率跟踪控制模块检测锁相放大器输出的解调幅值,若解调幅值V为0,则不改变中心频率,控制系统向上位机指示频率偏差为0,并继续检测锁相放大器输出的解调幅值直至解调幅值V不为0;否则执行步骤10。
步骤10,若解调幅值V大于保护阈值H,即此时解调幅值离开解调的线性区间,执行步骤11,否则执行步骤12。
步骤11,控制系统读出辅助磁场测量装置输出,此时输出满足:
其中,为当前时刻通过辅助磁场测量装置获取的外加待测磁场强度,/>是积分线圈的实时输出值,/>是积分线圈的面积。计算此时的频率偏差为:
,/>的初始值为零;
使中心频率f 1=f 0+ε1进行跟踪。
步骤12,若解调幅值V小于误差阈值W,执行步骤13,否则执行步骤17。
步骤13,锁相放大器解调幅值V累加得到累加值Vi,同时计数器的值n执行+1操作。
步骤14,若解调幅值V大于误差阈值W或者n大于最大累加次数Nmax(例如Nmax=5),执行步骤15,否则执行步骤12。
步骤15,计算频率偏差δ2’,计算公式如下:
δ2'=Vi/(n·K)
其中,n为计算器当前值,计算完成后计数器值n清零。
步骤16,使中心频率f 2=f 1+δ2'进行跟踪,执行步骤19(若执行完步骤10后,不执行步骤11而直接执行步骤12,此时f 2=f 0+δ2')。
步骤17,计算频率偏差,计算频率偏差δ2,计算公式如下:
δ2= V/ K
其中V为锁相放大器解调幅值,K为解调谱线零点处斜率,即线性区斜率。
步骤18,使中心频率f 2=f 1+δ2进行跟踪(若执行完步骤10后,不执行步骤11而直接执行步骤12,此时f 2=f 0+δ2')。
步骤19,计算总频率偏差ε=f 2-f 0,并根据计算得到磁场强度的变化量/>,将磁场强度的变化量/>与/>求和后赋值给/>。
步骤20,若仍检测到荧光信号,则执行步骤9,否则执行步骤21。
步骤21,结束跟踪调节过程。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种基于金刚石NV色心的量子传感频率跟踪控制方法,其特征在于,包括以下步骤:
S1,根据计算得到频率偏差ε1,其中,/>为磁旋比,/>为外加待测磁场与NV轴向的夹角,/>,/>为当前时刻通过辅助磁场测量装置获取的外加待测磁场强度,/>的初始值为零;将微波源的频率调整为f 1,其中,f 1=f 0+ε1,f 0为上一轮循环结束时微波源的频率,f 0的初始值为微波源起始频率;
S2,根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2,其中,f 2=f 1+δ2,所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;
S3,计算频率偏差ε=f 2-f 0,并根据计算得到磁场强度的变化量,将磁场强度的变化量/>与/>求和后赋值给/>,再返回S1继续下一轮频率跟踪控制。
2.根据权利要求1所述的基于金刚石NV色心的量子传感频率跟踪控制方法,其特征在于,将S2替换为S2',所述S2'包括:
S21',若f 1频率下的解调幅值V大于误差阈值W,则根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2并执行S3,其中,f 2=f 1+δ2,所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;否则,执行S22';
S22',对当前V进行累加,获得累加值Vi,同时计数器值n执行自加1操作;
S23',若累加值Vi大于误差阈值W或者计数器值n大于最大累计次数,根据δ2'=Vi/(n·K)计算频率偏差δ2'并将计数器值n清零,再将微波源的频率调整为f 2并执行S3,其中,f 2=f 1+δ2';否则,返回S21'。
3.一种基于金刚石NV色心的量子传感频率跟踪控制方法,其特征在于,包括以下步骤:
S0,若f 0频率下的解调幅值大于保护阈值,则执行S1;否则,执行S3;其中,f 0为上一轮循环结束时微波源的频率,f 0的初始值为微波源起始频率;
S1,根据计算得到频率偏差ε1,其中,/>为磁旋比,/>为外加待测磁场与NV轴向的夹角,/>,/>为当前时刻通过辅助磁场测量装置获取的外加待测磁场强度,/>的初始值为零;将微波源的频率调整为f 1,其中,f 1=f 0+ε1;
S2,根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2并执行S4,其中,f 2=f 1+δ2,所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;
S3,根据解调幅值与频率偏差的转换系数K计算f 0频率下的解调幅值对应的频率偏差δ2',再将微波源的频率调整为f 2并执行S4,其中,f 2=f 0+δ2',所述转换系数K通过获取f0前后的多个频率点对应的解调幅值拟合得到;
S4,计算频率偏差ε=f 2-f 0,并根据计算得到磁场强度的变化量,将磁场强度的变化量/>与/>求和后赋值给/>,再返回S0继续下一轮频率跟踪控制。
4.根据权利要求3所述的基于金刚石NV色心的量子传感频率跟踪控制方法,其特征在于,将S2替换为S2',所述S2'包括:
S21',若f 1频率下的解调幅值V大于误差阈值W,则根据解调幅值与频率偏差的转换系数K计算f 1频率下的解调幅值对应的频率偏差δ2,再将微波源的频率调整为f 2并执行S4,其中,f 2=f 1+δ2,所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;否则,执行S22';
S22',对当前V进行累加,获得累加值Vi,同时计数器值n执行自加1操作;
S23',若累加值Vi大于误差阈值W或者计数器值n大于最大累计次数,根据δ2''=Vi/(n·K)计算频率偏差δ2''并将计数器值n清零,再将微波源的频率调整为f 2并执行S4,其中,f 2=f 1+δ2'';否则,返回S21'。
5.根据权利要求3所述的基于金刚石NV色心的量子传感频率跟踪控制方法,其特征在于,将S3替换为S3',所述S3'包括:
S31',若f 0频率下的解调幅值V大于误差阈值W,则根据解调幅值与频率偏差的转换系数K计算f 0频率下的解调幅值对应的频率偏差δ2',再将微波源的频率调整为f 2并执行S4,其中,f 2=f 0+δ2',所述转换系数K通过获取f 0前后的多个频率点对应的解调幅值拟合得到;否则,执行S32';
S32',对当前V进行累加,获得累加值Vi,同时计数器值n执行自加1操作;
S33',若累加值Vi大于误差阈值W或者计数器值n大于最大累计次数,根据δ2'''=Vi/(n·K)计算频率偏差δ2'''并将计数器值n清零,再将微波源的频率调整为f 2并执行S4,其中,f 2=f 0+δ2''';否则,返回S31'。
6.一种基于金刚石NV色心的量子传感频率跟踪控制系统,其特征在于,包括:微波源、锁相放大器、辅助磁场测量装置、频率跟踪控制模块、NV色心传感头;
所述辅助磁场测量装置安装于所述NV色心传感头处,用于测量外加待测磁场的强度并输出至所述频率跟踪控制模块;所述频率跟踪控制模块还用于接收所述锁相放大器输出的解调幅值,以实现权利要求1至5任一项所述的基于金刚石NV色心的量子传感频率跟踪控制方法。
7.根据权利要求6所述的基于金刚石NV色心的量子传感频率跟踪控制系统,其特征在于,所述辅助磁场测量装置包括感应线圈和积分器,或者感应线圈、模数转换器和数字积分器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311124042.1A CN116859300B (zh) | 2023-09-01 | 2023-09-01 | 基于金刚石nv色心的量子传感频率跟踪控制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311124042.1A CN116859300B (zh) | 2023-09-01 | 2023-09-01 | 基于金刚石nv色心的量子传感频率跟踪控制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116859300A CN116859300A (zh) | 2023-10-10 |
CN116859300B true CN116859300B (zh) | 2023-11-17 |
Family
ID=88232629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311124042.1A Active CN116859300B (zh) | 2023-09-01 | 2023-09-01 | 基于金刚石nv色心的量子传感频率跟踪控制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116859300B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118091507B (zh) * | 2024-04-28 | 2024-07-23 | 国仪量子技术(合肥)股份有限公司 | 磁测量方法、装置及存储介质、控制器 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013055246A2 (en) * | 2011-10-14 | 2013-04-18 | Uniwersytet Jagiellonski | A method and a device for the measurement of changes in magnetic field |
EP2837930A1 (en) * | 2013-08-15 | 2015-02-18 | Latvijas Universitate | Method for detecting the magnitude of a magnetic field gradient |
CN109061295A (zh) * | 2018-06-29 | 2018-12-21 | 北京航空航天大学 | 一种近场微波谐振器谐振频率测量系统及方法 |
CN111521961A (zh) * | 2020-06-09 | 2020-08-11 | 中北大学 | 基于金刚石nv色心的单路三轴磁信息检测方法及系统 |
CN112462305A (zh) * | 2019-09-06 | 2021-03-09 | 北京华航无线电测量研究所 | 基于金刚石nv色心的微波线性调频测量装置及测量方法 |
CN114062982A (zh) * | 2020-07-29 | 2022-02-18 | 国仪量子(合肥)技术有限公司 | 基于上位机的nv色心磁测量装置、磁测量方法、存储介质 |
CN114200363A (zh) * | 2022-02-16 | 2022-03-18 | 国仪量子(合肥)技术有限公司 | 集成电路的磁场检测方法及缺陷检测方法、装置 |
CN114415080A (zh) * | 2022-01-24 | 2022-04-29 | 中国科学技术大学 | 金刚石nv色心磁力计及测量磁场的方法 |
CN114660512A (zh) * | 2022-03-02 | 2022-06-24 | 国仪量子(合肥)技术有限公司 | 基于金刚石nv色心探头的磁异常检测方法、介质、设备 |
WO2022242971A1 (de) * | 2021-05-20 | 2022-11-24 | Endress+Hauser SE+Co. KG | Detektionseinheit für magnetfeldsensor |
WO2022256992A1 (zh) * | 2021-06-07 | 2022-12-15 | 中国科学技术大学 | 基于金刚石nv- 15n耦合自旋体系的原子钟实现方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016118756A1 (en) * | 2015-01-23 | 2016-07-28 | Lockheed Martin Corporation | Apparatus and method for high sensitivity magnetometry measurement and signal processing in a magnetic detection system |
US10359479B2 (en) * | 2017-02-20 | 2019-07-23 | Lockheed Martin Corporation | Efficient thermal drift compensation in DNV vector magnetometry |
US10712408B2 (en) * | 2016-11-08 | 2020-07-14 | Massachusetts Institute Of Technology | Methods and apparatus for optically detecting magnetic resonance |
-
2023
- 2023-09-01 CN CN202311124042.1A patent/CN116859300B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013055246A2 (en) * | 2011-10-14 | 2013-04-18 | Uniwersytet Jagiellonski | A method and a device for the measurement of changes in magnetic field |
EP2837930A1 (en) * | 2013-08-15 | 2015-02-18 | Latvijas Universitate | Method for detecting the magnitude of a magnetic field gradient |
CN109061295A (zh) * | 2018-06-29 | 2018-12-21 | 北京航空航天大学 | 一种近场微波谐振器谐振频率测量系统及方法 |
CN112462305A (zh) * | 2019-09-06 | 2021-03-09 | 北京华航无线电测量研究所 | 基于金刚石nv色心的微波线性调频测量装置及测量方法 |
CN111521961A (zh) * | 2020-06-09 | 2020-08-11 | 中北大学 | 基于金刚石nv色心的单路三轴磁信息检测方法及系统 |
CN114062982A (zh) * | 2020-07-29 | 2022-02-18 | 国仪量子(合肥)技术有限公司 | 基于上位机的nv色心磁测量装置、磁测量方法、存储介质 |
WO2022242971A1 (de) * | 2021-05-20 | 2022-11-24 | Endress+Hauser SE+Co. KG | Detektionseinheit für magnetfeldsensor |
WO2022256992A1 (zh) * | 2021-06-07 | 2022-12-15 | 中国科学技术大学 | 基于金刚石nv- 15n耦合自旋体系的原子钟实现方法及装置 |
CN114415080A (zh) * | 2022-01-24 | 2022-04-29 | 中国科学技术大学 | 金刚石nv色心磁力计及测量磁场的方法 |
CN114200363A (zh) * | 2022-02-16 | 2022-03-18 | 国仪量子(合肥)技术有限公司 | 集成电路的磁场检测方法及缺陷检测方法、装置 |
CN114660512A (zh) * | 2022-03-02 | 2022-06-24 | 国仪量子(合肥)技术有限公司 | 基于金刚石nv色心探头的磁异常检测方法、介质、设备 |
Non-Patent Citations (1)
Title |
---|
金刚石晶体能带计算与光谱分析;李海东等;《量子电子学报》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116859300A (zh) | 2023-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116859300B (zh) | 基于金刚石nv色心的量子传感频率跟踪控制方法及系统 | |
CN104535530B (zh) | 一种高精度气体浓度检测方法及检测装置 | |
CN113376449B (zh) | 基于里德堡原子的低频微波电场传感器及探测方法 | |
CN111609918A (zh) | 基于包络检波电路的光纤分布式振动传感系统 | |
CN108964797A (zh) | 一种基于鉴频法的微波信号源相位噪声自动化测量系统及方法 | |
CN108459040B (zh) | 基于金刚石nv色心的磁悬浮加速度计的差分检测方法 | |
CN105577280A (zh) | 一种光载微波信号动态宽频实时数字解调系统 | |
CN107290753A (zh) | 一种激光测量系统及方法 | |
CN116087604A (zh) | 一种单圆环差分式磁通门电流传感器 | |
CN108931782A (zh) | 一种激光测距装置及其测距方法 | |
CN115561504A (zh) | 反射式全光纤电流传感器最优调制深度的调制信号确定方法 | |
CN113091617B (zh) | 一种多模光纤光程变化测量系统 | |
CN112729544A (zh) | 一种扫描吸收光谱的激光波长锁定系统与方法 | |
CN106375039A (zh) | 一种提高矢量网络分析仪接收机动态范围的方法 | |
CN112097924B (zh) | 基于相位生成载波技术的相位灵敏度标定方法 | |
CN204359684U (zh) | 一种高精度气体浓度检测装置 | |
CN114353685A (zh) | 基于混沌布里渊相位谱的高频动态应变测量装置和方法 | |
JP2008151676A (ja) | Esr装置 | |
CN109870677A (zh) | 一种调频连续波干涉信号幅度归一化方法 | |
RU2324947C1 (ru) | Устройство определения частоты и вида модуляции принимаемых сигналов | |
CN118311486A (zh) | 金刚石nv色心磁强测量方法及系统 | |
CN221123301U (zh) | 一种电涡流传感器 | |
CN219657765U (zh) | 一种高频信号功率和频率同步测量探头电路 | |
CN108233759B (zh) | 一种具有温度补偿系统的质谱仪射频电源 | |
CN116147801A (zh) | 一种基于大调制中心频率的外调制布里渊光相关域反射方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |