CN112003592B - 一种实现高分辨率量子传感的脉冲整形算法 - Google Patents

一种实现高分辨率量子传感的脉冲整形算法 Download PDF

Info

Publication number
CN112003592B
CN112003592B CN202010817380.3A CN202010817380A CN112003592B CN 112003592 B CN112003592 B CN 112003592B CN 202010817380 A CN202010817380 A CN 202010817380A CN 112003592 B CN112003592 B CN 112003592B
Authority
CN
China
Prior art keywords
pulse
time
resolution
data
sequence
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
Application number
CN202010817380.3A
Other languages
English (en)
Other versions
CN112003592A (zh
Inventor
刘俊
唐军
李中豪
郭浩
刘文耀
温焕飞
马宗敏
周彦汝
崔凌霄
张健
潘鼎文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North University of China
Original Assignee
North University of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North University of China filed Critical North University of China
Priority to CN202010817380.3A priority Critical patent/CN112003592B/zh
Publication of CN112003592A publication Critical patent/CN112003592A/zh
Application granted granted Critical
Publication of CN112003592B publication Critical patent/CN112003592B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03K—PULSE TECHNIQUE
    • H03K5/00—Manipulating of pulses not covered by one of the other main groups of this subclass
    • H03K5/01—Shaping pulses
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Physics & Mathematics (AREA)
  • Nonlinear Science (AREA)
  • Measurement Of Radiation (AREA)
  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)

Abstract

本发明涉及量子检测的微波脉冲领域，具体涉及一种实现高分辨率量子传感的脉冲整形算法，可应用于实现电子自旋的高光谱分辨、高频率分辨率交流磁场检测以及特定核磁共振信号的检测，本发明利用NV色心对磁信号和微波信号变化的高敏感度，通过算法改变脉冲电压幅值，能够突破现有信号发生器的自身硬件最高时间分辨率的限制，在高速量子操控中，实现更精确的量子相位翻转及控制。该方法实现简单，提升谱信噪比的同时还可以增强其应用于量子传感器的响应灵敏度，特别适用于需要多个π脉冲量子动态解耦技术以及高保真度量子调控领域当中。

Description

一种实现高分辨率量子传感的脉冲整形算法
技术领域
本发明涉及量子检测的微波脉冲领域，具体涉及一种实现高分辨率量子传感的脉冲整形算法。
背景技术
近年来，脉冲整形应用于量子调控和量子传感的动态解耦光谱学、核自旋的纳米光谱等新兴领域，具有提高信号灵敏度的作用；最简单的序列由一系列的周期性方波脉冲组成，可以在非常宽的频率范围内激发自旋，但不具备选择性激发的能力，故常采用脉冲整形，应用于高分辨率量子传感，令重复时间为τ，对于众多的脉冲数N，序列的光谱响应类似于窄带滤波器的频谱响应，其带宽为(1/Nτ)，中心频率为1/(2τ)，能抑制与重复时间相对应的频率以外的所有频率上的噪声，若调整τ使其在特定频率下与信号共振，其频率为Fac≈1/2τ，在提高噪声抑制作用的同时，可以提高信号的灵敏度。
当使用多个控制脉冲时，滤波器带宽会变窄，因此重复时间τ必须精确设置为信号频率Fac≈1/2τ，但对于任意信号发生器硬件，τ只能以采样率对应的时间t的增量进行调整，这在实验过程中会限制频率的分辨率；具体地说，当检测到具有频率Fac的信号时，最小频率增量由下式给出：
在我们的实验过程中，用于控制自旋量子位和超导的任意波形发生器具有典型的2.5GS/s的采样率，对应于时间分辨率为ts＝0.4ns，在高频下操作时，这种固定的时间分辨率无法满足实验需求，由于电子瓶颈的限制，如果直接在频域内产生短脉冲信号，将会受到带宽的限制而带来许多技术问题，从而提高大量成本；虽然存在采样率较高的脉冲发生器，但价格昂贵且基本只能勉强达到足够的时间分辨率，硬件的时间采样率是量子调控实验的一个严重限制。
发明内容
为解决上述技术问题，本发明提供一种实现高分辨率量子传感的脉冲整形算法，具体步骤如下：
S1：以列序为主的数组数据使用前，对所有行从小到大进行编号，同时对其编号另起一列存入数组中；
S2：以列序为主的数组的行数编号完成后，将该行编号与时间分辨率相乘的结果作为时间变量；
S3：在转换为行数的过程中，当需计算的时间不为时间分辨率整数倍时，对行数以四舍五入为原则进行处理；
S4：在转换后所对应的行数处输入的幅值是以转换前的时间所计算而得；
S5：使用的高斯脉冲为其中τ为单个高斯脉冲中的左右幅值降至的时间长度，t0为单个高斯脉冲的中心对称轴，Δt为时间分辨率的增量；
进一步地，t0＝c*τ，其中系数c可根据实验所需条件调整单个周期的高斯脉冲长度；脉冲有效宽度为其单个完整脉冲的半高全宽，令单次脉冲时长为tpw，半高全宽为tFWHM，则tpw＝2*t0，
通过绘制低浓度样品状态的金刚石NV色心的Rabi频率曲线来测试改良后脉冲的效果，具体方法如下：
a.首先调试系统，设置激光强度为15毫瓦，设置某一型号任意波形发生器的微波电压幅值为1.5v，基频设置为2.73GHz，使用的微波功率放大器增益为20dBm，将任意波形发生器设置为IQ调制模式，采样率设置为任意波形发生器的最大值2.5Gs/s，采用低浓度金刚石，进行实验；
b.设置同步系统的TTL信号的时序，每一周期初始化激光时长为50us，关断时间为6us，初始化激光的上升沿信号为Readout信号，序列周期为400个脉冲周期，通过改变幅值生成的高时间分辨率高斯脉冲脉宽从0.2ns步进到80ns，将制作好的序列加载至任意波形发生器的一个通道并执行运行，信号发生器将微波脉冲信号发送至射频功率放大器，微波脉冲信号经射频功率放大器放大后通过共振天线施加到低浓度的金刚石NV色心上；
c.接着进行Rabi振荡实验，使用数据采集卡采集光电探测器输出端数据；
d.通过编写算法，求解各周期微波频率下的荧光强度值，对于任一周期脉冲数求和为将每一脉冲周期的数据段首部去掉6个无用的上升沿数据点，即/>为保证信噪比，再求均值/>得出结果为这个周期微波频率下的荧光强度值；
e.处理后的全部脉冲周期的荧光强度数据绘图即为Rabi振荡谱；
f.重复步骤b至步骤e，比较普通高斯脉冲和高分辨率高斯脉冲测得的Rabi谱数据。
一种实现高分辨率量子传感的脉冲整形算法，适用于垂直分辨率高的任意波形发生器，可使用于任意编程环境中，脉冲幅值数据存储采用以列序为主的数组存储方式，按照行号从小到大的顺序，依次存储每一时刻的数据，当以列序为主的数组中无所选择的脉冲频率数据，将行值为1的列中的存储数据替换为新数据。
与现有技术相比，本发明取得的有益效果，通过特定的算法改变电压幅值，并结合脉冲整形优化技术可以使时间采样率超出硬件限制，能够实现更好的优化脉冲之间延迟处的自旋演化。
附图说明
图1微波脉冲调制示意图
图2原脉冲与幅值修改后脉冲的形状示意图
图3 Rabi振荡实验流程图
图4原脉冲与幅值修改后脉冲应用于Rabi振荡实验对比图
具体实施方案
下面结合附图对本发明做进一步的说明，以令本领域技术人员参照说明书文字能够实施。
一种实现高分辨率量子传感的脉冲整形算法，具体包括以下步骤：
S1：在以列序为主的数组数据使用前，对所有行从小到大进行编号，同时对其编号另起一列存入数组中；
S2：列序为主的数组的行数编号完成后，将该行编号与时间分辨率相乘的结果作为时间变量；
S3：在转换为行数的过程中，当需计算的时间不为时间分辨率整数倍时，对行数以四舍五入为原则进行处理；
动态解耦光谱学是基于这种有精确定时脉冲重复时间τ的量子比特控制的周期性脉冲调制的领域，如图1；
S4：在转换后所对应的行数处输入的幅值是以转换前的时间所计算而得；
S5：使用的高斯脉冲为其中τ为单个高斯脉冲中的左右幅值降至的时间长度，t0为单个高斯脉冲的中心对称轴，Δt为时间分辨率的增量；t0＝c*τ，其中系数c可根据实验所需条件调整单个周期的高斯脉冲长度；脉冲有效宽度为其单个完整脉冲的半高全宽，令单次脉冲时长为tpw，半高全宽为tFWHM，则/>如图2，用整形脉冲调制微波信号，将时间分辨率限制为脉冲发生器硬件采样时间ts的倍数，其中，正方体和圆形的轮廓显示原始和时移脉冲；
一种实现高分辨率量子传感的脉冲整形算法应用到脉冲ODMR光谱及Rabi振荡的检测中，如图3为进行相关量子调控实验的流程图，其具体步骤如下：
a.首先调试量子调控系统，使用功率为15毫瓦的激光器，选用某一型号任意波形发生器，设置任意波形发生器的微波电压幅值为1.5v，设置为IQ调制模式，采样率设置为任意波形发生器的最大值2.5Gs/s，将任意波形发生器的基频设置为所记录共振峰对应频率。使用的微波功率放大器增益为20dBm；
b.设置同步系统的TTL信号的时序，每周期初始化激光时长为50us，关断时间为6us，初始化激光的上升沿信号为Readout信号，序列周期为400个脉冲周期，通过改变幅值生成的高时间分辨率高斯脉冲脉宽从0.2ns步进到80ns。将制作好的序列加载至任意波形发生器的一个通道并执行运行，信号发生器将微波脉冲信号发送至射频功率放大器，微波脉冲信号经射频功率放大器放大后通过共振天线施加到低浓度的金刚石NV色心上；
c.接着进行Rabi振荡实验，使用数据采集卡采集光电探测器输出端数据；
d.通过编写算法，求解各周期微波频率下的荧光强度值，对于任一周期脉冲数求和为将每一脉冲周期的数据段首部去掉6个无用的上升沿数据点，即/>为保证信噪比，再求均值/>得出结果为这个周期微波频率下的荧光强度值；
e.处理后的全部脉冲周期的荧光强度数据绘图即为Rabi振荡谱；
f.重复步骤b至步骤e，比较普通高斯脉冲和高分辨率高斯脉冲测得的Rabi谱数据，如图4。
一种实现高分辨率量子传感的脉冲整形算法，适用于垂直分辨率高的任意波形发生器，可使用于任意编程环境中，其特征在于脉冲幅值数据存储采用以列序为主的数组存储方式，按照行号从小到大的顺序，依次存储每一时刻的数据，当以列序为主的数组中无所选择的脉冲频率数据，将行值为1的列中的存储数据替换为新数据。
以上实施例仅为本发明的优选技术方案，其并不仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的图例。

Claims (3)

1.一种实现高分辨率量子传感的脉冲整形算法，其特征在于，包括以下步骤：
S1：以列序为主的数组的数据使用前，对所述数组的所有行从小到大进行编号，同时将其编号另起一列存入数组中；
S2：当所述数组的行数编号完成后，将所述行编号与时间分辨率相乘的结果作为时间变量；
S3：当需计算的时间不为时间分辨率整数倍时，我们在转换为行数时以四舍五入为原则进行处理；
S4：在转换后所对应的行数处输入的幅值是以转换前的时间所计算而得；
S5：使用的高斯脉冲为其中τ为单个高斯脉冲中的左右幅值降至/>的时间长度,t0为单个高斯脉冲的中心对称轴，Δt为时间分辨率的增量；t0＝c*τ，其中系数c可根据实验所需条件调整单个周期的高斯脉冲长度；脉冲有效宽度为其单个完整脉冲的半高全宽，令单次脉冲时长为tpw，半高全宽为tFWHM，则/>
2.根据权利要求1所述的一种实现高分辨率量子传感的脉冲整形算法，所述脉冲整形算法应用到脉冲ODMR光谱及Rabi振荡的检测中，其特征在于，绘制金刚石NV色心Rabi频率曲线的具体方法如下：
a.首先调试量子调控系统，我们使用低浓度的金刚石NV色心，使用功率为15毫瓦的激光，选用某一型号任意波形发生器，设置任意波形发生器的微波电压幅值为1.5v，设置为IQ调制模式，采样率设置为任意波形发生器的最大值2.5Gs/s，将任意波形发生器的基频设置为所记录共振峰对应频率，使用的微波功率放大器增益为20dBm；
b.设置同步系统的TTL信号的时序，每周期初始化激光时长为50us，关断时间为6us，初始化激光的上升沿信号为Readout信号，序列周期为400个脉冲周期，通过改变幅值生成的高时间分辨率高斯脉冲脉宽从0.2ns步进到80ns；将制作好的序列加载至任意波形发生器的一个通道并执行运行；
c.接着进行Rabi振荡实验，使用数据采集卡采集光电探测器输出端数据；
d.通过编写算法，求解各周期微波频率下的荧光强度值；对于任一周期脉冲数的荧光强度，求和为其中N参与计算的总样点数，x(Ni)为第i个样点的模拟信号幅值强度，而当t≠Ni时，单位冲击信号δ(t-Ni)在整个区间内对时间t-Ni的积分为1，为将每一脉冲周期的数据段首部去掉6个无用的上升沿数据点，即为保证信噪比，再求均值/>得出结果为这个周期微波频率下的荧光强度值；
e.处理后的全部脉冲周期的荧光强度数据绘图即为Rabi振荡谱；
f.重复步骤b至步骤e，比较普通高斯脉冲和高分辨率高斯脉冲测得的Rabi谱数据。
3.根据权利要求1所述的一种实现高分辨率量子传感的脉冲整形算法，其特征在于，适用于垂直分辨率高的任意波形发生器，可使用于任意编程环境中，脉冲幅值数据存储采用以列序为主的数组存储方式，按照行号从小到大的顺序，依次存储每一时刻的数据，当以列序为主的数组中无所选择的脉冲频率数据，将行值为1的列中的存储数据替换为新数据。
CN202010817380.3A 2020-08-14 2020-08-14 一种实现高分辨率量子传感的脉冲整形算法 Active CN112003592B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010817380.3A CN112003592B (zh) 2020-08-14 2020-08-14 一种实现高分辨率量子传感的脉冲整形算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010817380.3A CN112003592B (zh) 2020-08-14 2020-08-14 一种实现高分辨率量子传感的脉冲整形算法

Publications (2)

Publication Number Publication Date
CN112003592A CN112003592A (zh) 2020-11-27
CN112003592B true CN112003592B (zh) 2023-09-12

Family

ID=73473291

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010817380.3A Active CN112003592B (zh) 2020-08-14 2020-08-14 一种实现高分辨率量子传感的脉冲整形算法

Country Status (1)

Country Link
CN (1) CN112003592B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113447863B (zh) * 2021-06-04 2022-06-03 电子科技大学 面向高频交变磁场的金刚石nv色心磁力仪频率测量方法
CN116593949B (zh) * 2023-01-03 2024-05-14 中国电力科学研究院有限公司 量子高速调控磁测量方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3787760A (en) * 1971-06-09 1974-01-22 Spectrospin Ag Method and apparatus for recording spin resonance spectra using two sequences of rf exciting pulses
JP2008108927A (ja) * 2006-10-26 2008-05-08 Nec Corp 量子ビット素子用読み出し方法、及びそれを適用した量子ビット素子用読み出し回路
JP2010099539A (ja) * 2010-02-10 2010-05-06 Toshiba Corp 磁気共鳴装置
CN104054266A (zh) * 2011-10-25 2014-09-17 中国科学院空间科学与应用研究中心 一种时间分辨单光子或极弱光多维成像光谱系统及方法
CN107271075A (zh) * 2017-06-27 2017-10-20 广州市禹成消防科技有限公司 基于光量子反演的光纤传感方法及系统
CN107356820A (zh) * 2017-06-07 2017-11-17 南京邮电大学 一种基于脉冲光探测磁共振的电磁场近场成像系统及方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005217837A (ja) * 2004-01-30 2005-08-11 Sony Corp サンプリングレート変換装置およびその方法、並びに、オーディオ装置
US7375802B2 (en) * 2005-08-04 2008-05-20 Lockheed Martin Corporation Radar systems and methods using entangled quantum particles
US7989775B2 (en) * 2005-08-04 2011-08-02 Lockheed Martin Corporation Sensor systems and methods using entangled quanta
US7583222B2 (en) * 2005-08-26 2009-09-01 Vaisala Oyj Method for using pulse compression in weather radar

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3787760A (en) * 1971-06-09 1974-01-22 Spectrospin Ag Method and apparatus for recording spin resonance spectra using two sequences of rf exciting pulses
JP2008108927A (ja) * 2006-10-26 2008-05-08 Nec Corp 量子ビット素子用読み出し方法、及びそれを適用した量子ビット素子用読み出し回路
JP2010099539A (ja) * 2010-02-10 2010-05-06 Toshiba Corp 磁気共鳴装置
CN104054266A (zh) * 2011-10-25 2014-09-17 中国科学院空间科学与应用研究中心 一种时间分辨单光子或极弱光多维成像光谱系统及方法
CN107356820A (zh) * 2017-06-07 2017-11-17 南京邮电大学 一种基于脉冲光探测磁共振的电磁场近场成像系统及方法
CN107271075A (zh) * 2017-06-27 2017-10-20 广州市禹成消防科技有限公司 基于光量子反演的光纤传感方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
闫树斌.SOI纳米光波导激光熔凝光整加工温度场的数值模拟及分析.《强激光与粒子束》.2017,98-105. *

Also Published As

Publication number Publication date
CN112003592A (zh) 2020-11-27

Similar Documents

Publication Publication Date Title
CN112003592B (zh) 一种实现高分辨率量子传感的脉冲整形算法
US4703267A (en) High dynamic range in NMR data acquisition
Mentink-Vigier et al. Increasing sensitivity of pulse EPR experiments using echo train detection schemes
RU2522179C2 (ru) Уменьшение sar в параллельной передаче посредством зависимого от каонного пространства выбора rf-импульсов
CN103763051A (zh) 实现瞬态信号捕获和频谱分析的系统
US3824452A (en) Fourier transform nuclear magnetic resonance spectrometer employing means for generating random pulse intervals to avoid anomalies due to net transverse magnetization
CN106772648A (zh) 提高质子旋进类传感器调谐精度和信噪比的方法及系统
JP2006522924A5 (zh)
CN101828911B (zh) 一种基于奇异谱熵的神经元动作电位特征提取方法
US3681680A (en) Rf spectrometer employing modulation of a dc magnetic field to excite resonance
CN116562379B (zh) 量子比特高能级跃迁频率的测量方法、装置和量子计算机
JP6994903B2 (ja) 電子スピン共鳴測定装置及び方法
Freeman The concept of spatially‐encoded single‐scan NMR
US4614907A (en) Method of obtaining pseudofiltering effect in process of data accumulation and nuclear magnetic resonance spectrometry utilizing same
CA2105677C (en) Method and apparatus for obtaining an nmr signal having a preselected frequency domain
JPS62103555A (ja) Ｎｍｒイメ−ジング装置
US7652475B2 (en) Cartesian continuous sampling with unequal gradients
CN109959852A (zh) 一种高速量测器件本征随机涨落特性的系统及方法
Souloumiac A Blind Reconstruction Algorithm for Level-Crossing Analog-to-Digital Conversion
CN115561271B (zh) 一种复合绝缘子老化核磁质子密度和弛豫谱联合检测方法
CN217305873U (zh) 一种在不完美系统中实施任意几何量子逻辑门的系统
HOCH et al. Maximum entropy reconstruction and nonuniform sampling in multidimensional NMR
CN118133982A (zh) 一种二维量子比特阵列耦合强度的相位校准方法
CN100365425C (zh) 一种具有选择性激发的正负相位组合脉冲的实现方法
CN117805075A (zh) 一种基于傅里叶微波脉冲激发的nv色心光探测磁共振系统

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