基于不变量绝热捷径的171Yb中性原子基态-里德伯态 Bell型纠缠态快速制备

上一篇

下一篇

王晓锋. 基于不变量绝热捷径的171Yb中性原子基态-里德伯态 Bell型纠缠态快速制备[J]. 物理学报, 2026, 75(16): 160602. doi: 10.7498/aps.75.20260364
引用本文: 王晓锋. 基于不变量绝热捷径的171Yb中性原子基态-里德伯态 Bell型纠缠态快速制备[J]. 物理学报, 2026, 75(16): 160602. doi: 10.7498/aps.75.20260364
Xiaofeng WANG. Fast preparation of Bell-type entangled states between the ground and Rydberg states of neutral 171Yb atoms based on shortcut to adiabaticity via invariants[J]. Acta Physica Sinica, 2026, 75(16): 160602. doi: 10.7498/aps.75.20260364
Citation: Xiaofeng WANG. Fast preparation of Bell-type entangled states between the ground and Rydberg states of neutral 171Yb atoms based on shortcut to adiabaticity via invariants[J]. Acta Physica Sinica, 2026, 75(16): 160602. doi: 10.7498/aps.75.20260364

基于不变量绝热捷径的171Yb中性原子基态-里德伯态 Bell型纠缠态快速制备

    作者简介: E-mail: andywangxiaofeng@163.com .
    通讯作者: E-mail: andywangxiaofeng@163.com
  • 中图分类号: 03.67.Bg, 32.80.Ee, 42.50.Dv, 03.65.Ta

Fast preparation of Bell-type entangled states between the ground and Rydberg states of neutral 171Yb atoms based on shortcut to adiabaticity via invariants

    Corresponding author: E-mail: andywangxiaofeng@163.com
  • MSC: 03.67.Bg, 32.80.Ee, 42.50.Dv, 03.65.Ta

  • 摘要: 针对171Yb中性原子里德伯阻塞体系中的快速纠缠态制备问题, 提出一种基于Lewis-Riesenfeld不变量逆向工程的绝热捷径方案, 用于实现基态-里德伯态对称单激发Bell型纠缠态的快速高保真制备. 在大失谐绝热消除与强里德伯阻塞近似下, 将原始双原子三能级系统约化为有效两能级模型, 并据此在共振驱动条件下构造解析控制脉冲. 数值结果表明, 在总演化时间T = 1 μs, 里德伯相互作用强度$V=2\text{π}\times 120$ MHz条件下, 系统可实现从初态到基态-里德伯态对称单激发Bell型纠缠态的快速转移, 末态保真度超过0.999, 所需峰值单原子有效两光子Rabi频率仅约为$2\text{π}\times0.66$ MHz. 进一步在完整双原子三能级模型中验证发现, 无耗散时末态保真度可达0.9999以上, 考虑中间态自发辐射和里德伯态衰减后仍保持在97%以上. 与传统绝热方案相比, 该方案可在更短时间内完成目标态制备, 并对驱动误差、失谐偏差、多普勒效应、激光线宽及相互作用波动表现出较好的鲁棒性. 结果表明, 该方案在171Yb里德伯中性原子平台上具有一定实验可行性, 并可为后续纠缠资源态向长寿命逻辑编码子空间的映射提供参考.
  • 加载中
  • 图 1  171Yb双原子系统及里德伯阻塞示意图 (a)两个原子被独立光镊囚禁, 间距为R; (b)双光子里德伯激发路径, 其中$ \left| {g_0} \right\rangle $$ \left| {g_1} \right\rangle $表示核自旋量子比特态, 参与里德伯激发的态$ \left| {g_1} \right\rangle $在后文简记为$ \left| {g} \right\rangle $, 经由中间态$ \left| {e} \right\rangle $耦合至里德伯态$ \left| {r} \right\rangle $, 中间态单光子失谐为$ \varDelta_{\mathrm{m}} $; (c)里德伯阻塞示意图, 范德瓦耳斯相互作用导致双激发态$ \left| {rr} \right\rangle $发生能级平移V, 从而有效抑制第二个原子的激发

    Figure 1.  Schematic of the physical model and Rydberg blockade for two 171Yb atoms: (a) Two atoms are independently trapped in optical tweezers with separation R; (b) two-photon ladder excitation from the nuclear-spin qubit state to the Rydberg state via the intermediate state $ \left| {e} \right\rangle $, with single-photon detuning $ \varDelta_{\mathrm{m}} $; (c) illustration of the Rydberg blockade, where the van der Waals interaction shifts the doubly excited state $ \left| {rr} \right\rangle $ by V, thereby suppressing the second excitation.

    图 2  STA 方案下的控制脉冲与系统布居演化 (a)根据有效两能级STA设计得到的控制脉冲有效两光子Rabi频率$ \varOmega(t)/2\text{π} $和失谐量$ \varDelta(t)/2\text{π} $随时间的变化, 其中失谐量在整个演化过程中保持为零; (b)在中间态绝热消除后的四维有效物理基底模型中得到的主要态布居演化, 蓝色虚线表示初始态$ \left| {G} \right\rangle =\left| {gg} \right\rangle $, 红色实线表示目标 Bell 型纠缠态$ \left| {B} \right\rangle $, 黑色点线表示双激发态$ \left| {rr} \right\rangle $

    Figure 2.  Control pulses and population dynamics for the STA scheme: (a) Time dependence of the effective two-photon Rabi frequency $ \varOmega(t)/2\text{π} $ and the detuning $ \varDelta(t)/2\text{π} $, where the detuning remains zero throughout the process; (b) population evolution of the main states, where the blue dashed, red solid, and black dotted lines denote $ \left| {G} \right\rangle =\left| {gg} \right\rangle $, the target Bell-type entangled state $ \left| {B} \right\rangle $, and the doubly excited state $ \left| {rr} \right\rangle $, respectively.

    图 3  末态保真度随原始双光子失谐$ \delta_{2\mathrm{ph}} $的变化关系. 横坐标为$ \delta_{2\mathrm{ph}}/2\text{π} $(MHz), 纵坐标为目标Bell型纠缠态的末态保真度. 结果表明, 在本文参数范围内, 系统的最优工作点位于$ \delta_{2\mathrm{ph}}/2\text{π} = 0\; \mathrm{MHz} $附近, 因此后续数值模拟均取该共振工作点

    Figure 3.  Final-state fidelity as a function of the initial two-photon detuning $ \delta_{2\mathrm{ph}} $. The horizontal axis denotes $ \delta_{2\mathrm{ph}}/2\text{π} $ (MHz), and the vertical axis denotes the final-state fidelity of the target Bell-type entangled state. The results show that, within the parameter range considered in this work, the optimal operating point is located near $ \delta_{2\mathrm{ph}}/2\text{π} = $$ 0\; \mathrm{MHz} $. Therefore, this resonant working point is adopted in all subsequent numerical simulations.

    图 4  完整模型与耗散效应下目标Bell型纠缠态制备的数值验证 (a)目标Bell型纠缠态保真度随时间的演化, 其中黑色实线表示有效两能级模型, 红色虚线表示无耗散时的完整双原子三能级模型, 蓝色点划线表示考虑中间态衰减与里德伯态有限寿命后的完整开放系统模型; (b)完整模型下主要态布居的时间演化, 其中$ P_{\mathrm{B}} $表示目标Bell型纠缠态布居, $ P_{{\mathrm{rr}}} $表示双里德伯激发态布居, $ P_{\mathrm{e}} $表示中间态总布居. 结果表明, 在完整模型中无耗散时末态保真度可达0.9999以上, 而考虑耗散后仍保持在97%以上, 同时$ P_{{\mathrm{rr}}} $$ P_{\mathrm{e}} $始终维持在很低水平

    Figure 4.  Numerical validation of the target Bell-type entangled state preparation in the full model with dissipation: (a) Time evolution of the target Bell-type entangled state fidelity, where the black solid line denotes the effective two-level model, the red dashed line denotes the full two-atom three-level model without dissipation, and the blue dash-dotted line denotes the full open-system model including the intermediate-state decay and finite Rydberg lifetime; (b) time evolution of the main populations in the full model, where $ P_{\mathrm{B}} $, $ P_{{\mathrm{rr}}} $, and $ P_{\mathrm{e}} $ denote the target Bell-type entangled state population, the doubly excited Rydberg-state population, and the total intermediate-state population, respectively. The results show that the final-state fidelity exceeds 0.9999 in the full model without dissipation and remains above 97% in the presence of dissipation, while both $ P_{{\mathrm{rr}}} $ and $ P_{\mathrm{e}} $ stay at very low levels throughout the process.

    图 5  传统绝热方案与本文STA方案的公平比较 (a)无耗散条件下, 传统绝热方案的末态保真度随总演化时间$ T_{\mathrm{ad}} $的变化. 虚线表示目标保真度阈值$ F=0.99 $, 对应所需总时间约为$ T_{\mathrm{ad}}=8.79\; \text{μs}$; (b)开放系统中两种方案末态保真度的比较, 其中 STA 方案取$ T=1.00$ μs, 传统绝热方案取$ T_{\mathrm{ad}}=8.79\; \text{μs} $. 可以看出, 随着里德伯态衰减率$ \varGamma_{\mathrm{r}} $的增大, STA方案始终保持更高的末态保真度

    Figure 5.  Fair comparison between the conventional adiabatic protocol and the present STA protocol: (a) Final-state fidelity of the conventional adiabatic protocol as a function of the total evolution time $ T_{\mathrm{ad}} $ in the absence of dissipation. The dashed line denotes the target fidelity threshold $ F=0.99 $, corresponding to the required duration $ T_{\mathrm{ad}}=8.79\; \text{μs} $; (b) comparison of the final-state fidelities between the two protocols in the open-system case. Here the STA duration is fixed at $ T=1.00$ μs, while the adiabatic duration is chosen as $ T_{\mathrm{ad}}=8.79\; \text{μs} $. It is seen that the STA protocol consistently maintains a higher final-state fidelity as the Rydberg decay rate $ \varGamma_{\mathrm{r }}$ increases.

    图 6  不同实验误差下的鲁棒性分析 (a)末态保真度随里德伯相互作用相对误差$ \delta_{\mathrm{V }}$的变化; (b)末态保真度随多普勒失谐的变化; (c)有限激光线宽条件下经Monte Carlo平均得到的末态保真度, 误差条表示Monte Carlo平均的标准误差

    Figure 6.  Robustness analysis under different experimental imperfections: (a) Final state fidelity as a function of the relative error of Rydberg interaction $ \delta_{\mathrm{V}} $; (b) final-state fidelity as a function of the Doppler detuning; (c) final state fidelity obtained by Monte Carlo averaging under finite laser linewidth, with error bars representing the standard error of the Monte Carlo average.

    图 7  STA方案在相对Rabi振幅误差$ \delta_{\varOmega} $与归一化失谐误差$ \delta_{\varDelta}/\varOmega_{\max} $下的二维鲁棒性分布, 其中$ \varOmega_{\max} $为本文STA方案的峰值有效两光子Rabi频率. 横坐标以百分比形式给出. 白色虚线分别对应保真度$ F=0.990 $$ F=0.999 $的等值线

    Figure 7.  Two-dimensional robustness map of the STA scheme versus the relative Rabi-amplitude error $ \delta_{\varOmega} $ and the normalized detuning error $ \delta_{\varDelta}/\varOmega_{\max} $, where $ \varOmega_{\max} $ denotes the peak effective two-photon Rabi frequency of the STA protocol. The horizontal axis is given in percentage. The white dashed curves represent the contours of fidelity $ F=0.990 $ and $ F=0.999 $, respectively.

    图 8  在更贴近实验条件的两类非理想因素下, 本文方案的鲁棒性分析结果 (a)在两原子耦合不对称条件下, 目标Bell型纠缠态末态保真度及反对称单激发态泄漏随不对称参数$ \epsilon $的变化, 其中$ \varOmega_A(t)=\varOmega(t)(1+\epsilon) $, $ \varOmega_B(t)=\varOmega(t)(1-\epsilon) $; (b)在时变里德伯相互作用条件下, 末态保真度随调制深度η的变化, 其中采用模型$ V(t)=V_0[1+\eta\sin(2\text{π} \nu t+\phi)] $, 并给出了两个代表性调制频率下的结果

    Figure 8.  Presents the robustness analysis of the proposed scheme under two types of nonideal effects that are closer to realistic experimental conditions: (a) Under asymmetric atom–laser couplings, the final fidelity of the target Bell state and the leakage to the antisymmetric singly excited state are shown as functions of the asymmetry parameter $\epsilon$, where $ \varOmega_A(t)=\varOmega(t)(1+\epsilon) $ and $ \varOmega_B(t)=\varOmega(t)(1-\epsilon) $; (b) under a time-dependent Rydberg interaction, the final-state fidelity is plotted as a function of the modulation depth η, using the model $ V(t)=V_0[1+\eta\sin(2\text{π} \nu t+\phi)] $, with results shown for two representative modulation frequencies.

    表 1  主要符号说明. 为避免不同层次模型中的符号混淆, 本文对主要参数的定义进行统一说明

    Table 1.  Summary of the main symbols used in different levels of the model. To avoid confusion among different model hierarchies, the definitions of the main parameters are listed here.

    符号 所属模型/层次 物理含义
    $ \varOmega_1(t), \, \varOmega_2(t) $ 原始/完整双原子三能级模型 两束单光子驱动场的Rabi频率
    $ \varDelta_{\mathrm{m}} $ 原始/完整双原子三能级模型 中间态单光子失谐
    $ \delta_{2\mathrm{ph}} $ 原始/完整双原子三能级模型 裸双光子失谐
    $ \varOmega(t) $ 中间态消除后的单原子有效模型 单原子有效两光子Rabi频率
    $ \varOmega_{\mathrm{eff}}(t) $ 对称双原子有效两能级模型 $ \left| {G} \right\rangle \leftrightarrow\left| {B} \right\rangle $之间的集体有效Rabi频率
    $ \varDelta(t) $ 有效两能级模型 含裸双光子失谐、AC Stark位移及有限阻塞修正的总有效失谐
    V 各层次模型共用 双里德伯激发态$ \left| {rr} \right\rangle $的相互作用能移
    $ \delta_V, \delta_\varOmega, \delta_\varDelta $ 误差与鲁棒性分析模型 分别表示相互作用强度、驱动幅度和有效失谐的相对偏差或扰动参数
    下载: 导出CSV
  • [1] Preskill J 2018 Quantum 2 79 doi: 10.22331/q-2018-08-06-79
    [2] Kaufman A M, Ni K K 2021 Nat. Phys. 17 1324 doi: 10.1038/s41567-021-01357-2
    [3] Henriet L, Béguin L, Signoles A, Lahaye T, Browaeys A, Reymond G O, Jurczak C 2020 Quantum 4 327 doi: 10.22331/q-2020-09-21-327
    [4] Browaeys A, Lahaye T 2020 Nat. Phys. 16 132 doi: 10.1038/s41567-019-0733-z
    [5] Ebadi S, Wang T T, Levine H, Keesling A, Semeghini G, Omran A, Bluvstein D, Samajdar R, Pichler H, Ho W W, Choi S, Sachdev S, Greiner M, Vuletić V, Lukin M D 2021 Nature 595 227 doi: 10.1038/s41586-021-03582-4
    [6] Endres M, Bernien H, Keesling A, Levine H, Anschuetz E R, Krajenbrink A, Senko C, Vuletić V, Greiner M, Lukin M D 2016 Science 354 1024 doi: 10.1126/science.aah3752
    [7] Barredo D, de Léséleuc S, Lienhard V, Lahaye T, Browaeys A 2016 Science 354 1021 doi: 10.1126/science.aah3778
    [8] Bernien H, Schwartz S, Keesling A, Levine H, Omran A, Pichler H, Choi S, Zibrov A S, Endres M, Greiner M, Vuletić V, Lukin M D 2017 Nature 551 579 doi: 10.1038/nature24622
    [9] Wilk T, Gaëtan A, Evellin C, Wolters J, Miroshnychenko Y, Grangier P, Browaeys A 2010 Phys. Rev. Lett. 104 010502 doi: 10.1103/PhysRevLett.104.010502
    [10] Isenhower L, Urban E, Zhang X L, Gill A T, Henage T, Johnson T A, Walker T G, Saffman M 2010 Phys. Rev. Lett. 104 010503 doi: 10.1103/PhysRevLett.104.010503
    [11] Levine H, Keesling A, Omran A, Bernien H, Schwartz S, Zibrov A S, Endres M, Greiner M, Vuletić V, Lukin M D 2018 Phys. Rev. Lett. 121 123603 doi: 10.1103/PhysRevLett.121.123603
    [12] Graham T M, Song Y, Scott J, Poole C, Phung D, Biswas A, Scholl P, Saffman M 2019 Phys. Rev. Lett. 123 230501 doi: 10.1103/PhysRevLett.123.230501
    [13] Saffman M, Walker T G, Mølmer K 2010 Rev. Mod. Phys. 82 2313 doi: 10.1103/RevModPhys.82.2313
    [14] Jaksch D, Cirac J I, Zoller P, Rolston S L, Côté R, Lukin M D 2000 Phys. Rev. Lett. 85 2208 doi: 10.1103/PhysRevLett.85.2208
    [15] Lukin M D, Fleischhauer M, Côté R, Duan L M, Jaksch D, Cirac J I, Zoller P 2001 Phys. Rev. Lett. 87 037901 doi: 10.1103/PhysRevLett.87.037901
    [16] Saffman M 2016 J. Phys. B: At. Mol. Opt. Phys. 49 202001 doi: 10.1088/0953-4075/49/20/202001
    [17] Adams C S, Pritchard J D, Shaffer J P 2020 J. Phys. B: At. Mol. Opt. Phys. 53 012002 doi: 10.1088/1361-6455/ab52ef
    [18] 张正源, 张天乙, 刘宗凯, 丁冬生, 史保森 2020 物理学报 69 180301 doi: 10.7498/aps.69.20200649 Zhang Z Y, Zhang T Y, Liu Z K, Ding D S, Shi B S 2020 Acta Phys. Sin. 69 180301 doi: 10.7498/aps.69.20200649
    [19] Morgado M, Whitlock S 2021 AVS Quantum Sci. 3 023501 doi: 10.1116/5.0036562
    [20] Maller K M, Lichtman M T, Xia T, Sun Y, Piotrowicz M J, Carr A W, Isenhower L, Saffman M 2015 Phys. Rev. A 92 022336 doi: 10.1103/PhysRevA.92.022336
    [21] Evered S J, Bluvstein D, Kalinowski M, Ebadi S, Manovitz T, Zhou H, Li S H, Geim A A, Wang T T, Maskara N, Levine H, Semeghini G, Greiner M, Vuletić V, Lukin M D 2023 Nature 622 268 doi: 10.1038/s41586-023-06481-y
    [22] 白文杰, 严冬, 韩海燕, 华硕, 谷开慧 2022 物理学报 71 014202 doi: 10.7498/aps.71.20211284 Bai W J, Yan D, Han H Y, Hua S, Gu K H 2022 Acta Phys. Sin. 71 014202 doi: 10.7498/aps.71.20211284
    [23] Su S L 2018 Chin. Phys. B 27 110304 doi: 10.1088/1674-1056/27/11/110304
    [24] Wang C, Liu B B, Yun M R, Zhu X Y, Su S L 2024 Int. J. Theor. Phys. 63 229 doi: 10.1007/s10773-024-05746-0
    [25] Xu J Z, Sun L N, Wei J F, Du Y L, Luo R, Yan L L, Feng M, Su S L 2022 Chin. Phys. Lett. 39 090301 doi: 10.1088/0256-307X/39/9/090301
    [26] Guo F Q, Wu J L, Zhu X Y, Jin Z, Zeng Y, Zhang S, Yan L L, Feng M, Su S L 2020 Phys. Rev. A 102 062410 doi: 10.1103/PhysRevA.102.062410
    [27] Daley A J, Ye J, Zoller P 2008 Phys. Rev. Lett. 101 170504 doi: 10.1103/PhysRevLett.101.170504
    [28] Gorshkov A V, Hermele M, Gurarie V, Xu C, Julienne P S, Ye J, Zoller P, Demler E, Lukin M D, Rey A M 2010 Nat. Phys. 6 289 doi: 10.1038/nphys1535
    [29] Derevianko A, Katori H 2011 Rev. Mod. Phys. 83 331 doi: 10.1103/RevModPhys.83.331
    [30] Ludlow A D, Boyd M M, Ye J, Peik E, Schmidt P O 2015 Rev. Mod. Phys. 87 637 doi: 10.1103/RevModPhys.87.637
    [31] Madjarov I S, Covey J P, Shaw A L, Choi J, Kale A, Cooper A, Pichler H, Schkolnik V, Williams J R, Endres M 2020 Nat. Phys. 16 857 doi: 10.1038/s41567-020-0903-z
    [32] Saskin S, Wilson J T, Grinkemeyer B, Thompson J D 2019 Phys. Rev. Lett. 122 143002 doi: 10.1103/PhysRevLett.122.143002
    [33] Jenkins A, Lis J W, Senoo A, McGrew W F, Kaufman A M 2022 Phys. Rev. X 12 021027 doi: 10.1103/PhysRevX.12.021027
    [34] Cooper A, Covey J P, Madjarov I S, Porsev S G, Safronova M S, Endres M 2018 Phys. Rev. X 8 041055 doi: 10.1103/PhysRevX.8.041055
    [35] Ma S, Burgers A P, Liu G Y, Wilson J, Zhang B, Thompson J D 2022 Phys. Rev. X 12 021028 doi: 10.1103/PhysRevX.12.021028
    [36] Ma S, Liu G, Peng P, Zhang B, Jandura S, Claes J, Burgers A P, Pupillo G, Puri S, Thompson J D 2023 Nature 622 279 doi: 10.1038/s41586-023-06438-1
    [37] Wilson J T, Saskin S, Grinkemeyer B, Meng Y, Thompson J D 2022 Phys. Rev. Lett. 128 033201 doi: 10.1103/PhysRevLett.128.033201
    [38] Scholl P, Williams H J, Bornet G, Wallner F, Barredo D, Lienhard V, Signoles A, Hainaut C, Franz T, Geier S, Tebben A, Salzinger A, Zürn G, Lahaye T, Weidemüller M, Browaeys A 2022 PRX Quantum 3 020303 doi: 10.1103/PRXQuantum.3.020303
    [39] Bergmann K, Theuer H, Shore B W 1998 Rev. Mod. Phys. 70 1003 doi: 10.1103/RevModPhys.70.1003
    [40] Vitanov N V, Rangelov A A, Shore B W, Bergmann K 2017 Rev. Mod. Phys. 89 015006 doi: 10.1103/RevModPhys.89.015006
    [41] Theis L S, Motzoi F, Wilhelm F K, Saffman M 2016 Phys. Rev. A 94 032306 doi: 10.1103/PhysRevA.94.032306
    [42] Guéry-Odelin D, Ruschhaupt A, Kiely A, Torrontegui E, Martínez-Garaot S, Muga J G 2019 Rev. Mod. Phys. 91 045001 doi: 10.1103/RevModPhys.91.045001
    [43] Torrontegui E, Ibáñez S, Martínez-Garaot S, Modugno M, del Campo A, Guéry-Odelin D, Ruschhaupt A, Chen X, Muga J G 2013 Adv. At. Mol. Opt. Phys. 62 117 doi: 10.1016/B978-0-12-408090-4.00002-5
    [44] Berry M V 2009 J. Phys. A: Math. Theor. 42 365303 doi: 10.1088/1751-8113/42/36/365303
    [45] Demirplak M, Rice S A 2003 J. Phys. Chem. A 107 9937 doi: 10.1021/jp030708a
    [46] Demirplak M, Rice S A 2005 J. Phys. Chem. B 109 6838 doi: 10.1021/jp040647w
    [47] Lewis H R, Riesenfeld W B 1969 J. Math. Phys. 10 1458 doi: 10.1063/1.1664991
    [48] Chen X, Ruschhaupt A, Schmidt S, del Campo A, Guéry-Odelin D, Muga J G 2010 Phys. Rev. Lett. 104 063002 doi: 10.1103/PhysRevLett.104.063002
    [49] Chen X, Torrontegui E, Muga J G 2011 Phys. Rev. A 83 062116 doi: 10.1103/PhysRevA.83.062116
    [50] Ruschhaupt A, Chen X, Alonso D, Muga J G 2012 New J. Phys. 14 093040 doi: 10.1088/1367-2630/14/9/093040
    [51] Bason M G, Viteau M, Malossi N, Huillery P, Arimondo E, Ciampini D, Fazio R, Giovannetti V, Mannella R, Morsch O 2012 Nat. Phys. 8 147 doi: 10.1038/nphys2170
    [52] del Campo A 2013 Phys. Rev. Lett. 111 100502 doi: 10.1103/PhysRevLett.111.100502
    [53] Kang Y H, Chen Y H, Wu Q C, Huang B H, Song J, Xia Y 2016 Sci. Rep. 6 36737 doi: 10.1038/srep36737
    [54] 张春玲, 刘文武 2018 物理学报 67 160302 doi: 10.7498/aps.67.20180315 Zhang C L, Liu W W 2018 Acta Phys. Sin. 67 160302 doi: 10.7498/aps.67.20180315
    [55] Kuwamoto T, Honda K, Takahashi Y, Yabuzaki T 1999 Phys. Rev. A 60 R745 doi: 10.1103/PhysRevA.60.R745
    [56] Takasu Y, Maki K, Komori K, Takano T, Honda K, Kumakura M, Yabuzaki T, Takahashi Y 2003 Phys. Rev. Lett. 91 040404 doi: 10.1103/PhysRevLett.91.040404
    [57] Halter C, Miethke A, Sillus C, Hegde A, Görlitz A 2023 J. Phys. B: At. Mol. Opt. Phys. 56 055001 doi: 10.1088/1361-6455/acb74b
    [58] Walker T G, Saffman M 2008 Phys. Rev. A 77 032723 doi: 10.1103/PhysRevA.77.032723
    [59] Beterov I I, Tretyakov D B, Entin V M, Ekers A, Bezuglov N N, Ryabtsev I I, Miculis K, Mahrov B, Ferber R, Auzinsh M, Schmieder H 2009 Phys. Rev. A 79 052504 doi: 10.1103/PhysRevA.79.052504
  • 加载中
图( 8) 表( 1)
计量
  • 文章访问数:  34
  • HTML全文浏览数:  34
  • PDF下载数:  1
  • 施引文献:  0
出版历程
  • 收稿日期:  2026-03-13
  • 刊出日期:  2026-08-20

基于不变量绝热捷径的171Yb中性原子基态-里德伯态 Bell型纠缠态快速制备

    通讯作者: E-mail: andywangxiaofeng@163.com
    作者简介: E-mail: andywangxiaofeng@163.com
  • 太原师范学院物理系, 晋中 030619

摘要: 针对171Yb中性原子里德伯阻塞体系中的快速纠缠态制备问题, 提出一种基于Lewis-Riesenfeld不变量逆向工程的绝热捷径方案, 用于实现基态-里德伯态对称单激发Bell型纠缠态的快速高保真制备. 在大失谐绝热消除与强里德伯阻塞近似下, 将原始双原子三能级系统约化为有效两能级模型, 并据此在共振驱动条件下构造解析控制脉冲. 数值结果表明, 在总演化时间T = 1 μs, 里德伯相互作用强度$V=2\text{π}\times 120$ MHz条件下, 系统可实现从初态到基态-里德伯态对称单激发Bell型纠缠态的快速转移, 末态保真度超过0.999, 所需峰值单原子有效两光子Rabi频率仅约为$2\text{π}\times0.66$ MHz. 进一步在完整双原子三能级模型中验证发现, 无耗散时末态保真度可达0.9999以上, 考虑中间态自发辐射和里德伯态衰减后仍保持在97%以上. 与传统绝热方案相比, 该方案可在更短时间内完成目标态制备, 并对驱动误差、失谐偏差、多普勒效应、激光线宽及相互作用波动表现出较好的鲁棒性. 结果表明, 该方案在171Yb里德伯中性原子平台上具有一定实验可行性, 并可为后续纠缠资源态向长寿命逻辑编码子空间的映射提供参考.

English Abstract

    • 量子计算与量子信息处理是当前物理学与信息科学交叉领域的重要研究方向之一[1]. 在众多物理实现平台中, 基于中性原子的光镊阵列体系因具有良好的可扩展性、灵活可重构的几何结构以及较长的相干时间, 已成为实现量子模拟、量子逻辑门和多体量子态调控的重要候选平台[28]. 近年来, 随着单原子囚禁、单站点操控、高保真读出以及里德伯相互作用调控等关键实验技术的持续发展, 中性原子体系在大规模量子信息处理中的应用前景日益受到关注[912].

      在中性原子量子计算方案中, 里德伯阻塞机制为实现长程、快速且可控的双体相互作用提供了重要手段[1318]. 当一个原子被激发到里德伯态后, 原子间强烈的范德瓦耳斯相互作用会使双激发态能级发生显著平移, 从而有效抑制相邻原子的再次激发, 这一效应可用于实现确定性纠缠态制备和双比特量子逻辑操作. 基于该机制的中性原子体系不仅已广泛用于纠缠态制备, 也被应用于双比特及多比特量子逻辑门、受控相位操作和可扩展量子信息处理任务, 并在理论和实验上取得了重要进展[925].

      值得指出的是, 里德伯中性原子体系的研究意义不仅体现在强相互作用可控平台的构建方面, 更在于其已成为量子信息处理的重要候选体系. 依托里德伯阻塞诱导的可编程相互作用, 该平台已被广泛用于量子逻辑门、几何量子操作、多体纠缠操控以及其他量子信息任务[2326]. 例如, 基于里德伯相互作用的受控相位门和双比特几何门方案展示了该体系在量子逻辑控制中的应用潜力[23,25]; 而围绕多体里德伯纠缠的鲁棒操控与判别, 以及基于几何量子操作的量子信息处理任务等研究, 则进一步拓展了其应用范围[24,26]. 从这一背景出发, 发展面向具体原子平台的快速、高保真且实验简洁的纠缠资源态制备方案, 对于后续逻辑态映射和相关量子逻辑操作具有重要意义.

      与常用的碱金属原子相比, 以镱为代表的类碱土金属原子在量子信息处理中展现出独特优势[2731]. 特别地, 171Yb具有核自旋$ I=1/2 $且电子基态壳层闭合, 量子比特可编码于$ ^{1}{\mathrm{S}}_0 $基态的核自旋自由度中, 因此对磁场噪声具有较强抑制能力, 同时能级结构相对简洁, 便于精密操控[3234]. 此外, 镱原子同时具备适用于激光冷却、窄线操控和里德伯激发的丰富光学跃迁资源, 为其在中性原子量子计算中的应用奠定了基础[2938].

      然而, 在实际纠缠态制备过程中, 系统演化速度与制备保真度之间通常存在明显折中. 传统绝热方法, 如受激拉曼绝热过程(STIRAP)或失谐扫描方案, 虽然能够在一定程度上抑制非目标跃迁, 但往往要求较长演化时间以满足绝热条件[39,40]. 过长的操作时间会导致里德伯态自发辐射、中间态散射、多普勒频移、激光线宽及强度涨落等误差源的累积, 从而限制最终制备保真度[1921,41]. 因此, 如何在有限时间内实现快速、高保真且鲁棒的纠缠态制备, 仍是中性原子量子调控研究中的关键问题.

      为克服传统绝热控制的时间瓶颈, 绝热捷径(shortcuts to adiabaticity, STA)方法近年来受到了广泛关注[4249]. STA的基本思想是在有限时间内构造与绝热演化等效的目标态转移路径, 从而在显著缩短演化时间的同时抑制非绝热泄漏. 其中, 基于Lewis-Riesenfeld不变量的逆向工程方法具有解析结构清晰、边界条件可控以及实验参数设计灵活等优点, 已在原子输运、布居反转和量子态操控等问题中得到广泛应用[4754]. 尽管如此, 针对171Yb里德伯中性原子平台中基态-里德伯态Bell型纠缠态快速制备问题, 尤其是在兼顾强阻塞、低驱动幅度和实验实现简洁性的前提下, 仍有进一步研究空间.

      基于此, 本文研究171Yb中性原子里德伯阻塞体系中的快速纠缠态制备问题, 提出一种基于 Lewis-Riesenfeld不变量逆向工程的绝热捷径控制方案. 需要说明的是, 本文直接制备的目标态为由参与里德伯耦合的基态$ \left| {g} \right\rangle \equiv\left| {g_1} \right\rangle $与里德伯态$ \left| {r} \right\rangle $构成的对称单激发Bell型纠缠态, 而非已经完全存储在长寿命核自旋逻辑基底$ \{\left| {g_0} \right\rangle, \left| {g_1} \right\rangle \} $中的最终逻辑Bell态. 本文首先在大失谐绝热消除与强里德伯阻塞近似下, 将原始三能级双原子系统约化为有效两能级模型; 随后在共振驱动条件下, 基于 Lewis-Riesenfeld不变量逆向工程构造解析控制脉冲. 对于本文采用的对称双光子驱动, 在本文参数范围内所需的有效共振工作点可通过预先谱线标定在裸双光子共振附近选取, 从而避免复杂的实时失谐扫描, 并完成从初态到目标基态–里德伯态对称单激发Bell型纠缠态的快速高保真转移. 与传统绝热方案相比, 该方案可在微秒时间尺度内完成纠缠态制备, 并保持较低驱动幅度, 从而有望简化实验控制要求. 进一步地, 本文系统分析了激光幅度误差、失谐偏差、多普勒效应、激光线宽以及相互作用波动对制备保真度的影响, 表明该方案具有较好的鲁棒性, 在本文讨论的参数范围内呈现出一定的实验可行性.

      与已有相关工作相比, 本文工作的特点主要体现在以下三个方面: 1)面向171Yb里德伯中性原子平台开展基态-里德伯态对称单激发Bell型纠缠态的快速制备研究, 突出该体系在核自旋编码和抗磁噪声方面的优势; 2)在共振驱动条件下基于Lewis-Riesenfeld不变量构造解析控制脉冲, 避免了复杂频率扫描, 提高了实验实现的简洁性; 3)在微秒时间尺度和较弱驱动幅度下实现了高保真且具有较强鲁棒性的纠缠资源态制备, 为后续逻辑态映射和量子门操作提供了可扩展的控制思路.

    • 本文对物理模型中涉及的主要符号及参数进行了统一的定义, 其所属模型层次及具体物理含义如表1所示. 需要指出的是, 原始双光子激发模型及完整双原子三能级模型中的裸双光子失谐统一记为$ \delta_{2\mathrm{ph}} $; 而在有效两能级模型中, 所有由裸双光子失谐、AC Stark位移等修正共同贡献的总失谐统一记为$ \varDelta(t) $, 以避免不同层次模型中的符号混淆.

    • 图1(a)所示, 本文考虑两个被独立光镊囚禁的全同171Yb原子, 原子间距记为R. 在实际实验中, 从基态到里德伯态的直接单光子跃迁通常受选择定则限制, 因此多采用双光子梯形激发方案. 图1(b)给出了相关能级结构: 原子的两个核自旋态$ \left| {g_0} \right\rangle $$ \left| {g_1} \right\rangle $可作为量子比特编码态, 其中参与里德伯激发的态记为$ \left| {g_1} \right\rangle $, 后文简记为$ \left| {g} \right\rangle $; 通过两束单光子Rabi频率分别为$ \varOmega_1(t) $$ \varOmega_2(t) $的激光, 经由中间态$ \left| {e} \right\rangle $激发至里德伯态$ \left| {r} \right\rangle $. 在中间态大失谐条件下, 该双光子梯形激发过程可近似等效为基态$ \left| {g} \right\rangle $与里德伯态$ \left| {r} \right\rangle $之间的有效两光子耦合. 需要特别说明的是, $ \left| {g_0} \right\rangle $$ \left| {g_1} \right\rangle $描述的是171Yb基态流形中的长寿命核自旋逻辑编码态, 而本文动力学中实际参与里德伯耦合的是其中一个编码态$ \left| {g_1} \right\rangle $, 后文简记为$ \left| {g} \right\rangle $. 因此, 本文直接制备的目标纠缠态属于$ \left| {g} \right\rangle $$ \left| {r} \right\rangle $张成的基态–里德伯态子空间中的对称单激发Bell型纠缠态, 其物理性质不同于完全存储在$ \{\left| {g_0} \right\rangle, \left| {g_1} \right\rangle \} $逻辑计算基底中的最终逻辑Bell态.

      图1(c)给出了里德伯阻塞机制的物理图像. 当其中一个原子已被激发到里德伯态时, 由于强烈的范德瓦耳斯相互作用, 双激发态$ \left| {rr} \right\rangle $的能级将发生显著平移, 使得第二个原子的激发偏离共振条件, 从而有效抑制双里德伯激发. 该机制是中性原子体系中实现确定性纠缠态制备和双比特量子门的核心基础[1317].

      在旋波近似下, 双原子系统哈密度量可写为

      其中$ \varDelta_{\mathrm{m}} $为中间态单光子失谐, $ \delta_{2\mathrm{ph}} $为双光子失谐, $ V=C_6/R^6 $为双原子同时处于里德伯态时的相互作用能移, 本文主要关注其幅值满足强阻塞条件时对动力学的影响[26,27].

      本文关注从集体基态

      向对称单激发Bell型纠缠态

      的高保真转移. 需要强调的是, (3)式所定义的$ \left| {B} \right\rangle $是由基态$ \left| {g} \right\rangle $与里德伯态$ \left| {r} \right\rangle $构成的对称单激发 Bell型纠缠态, 而非完全存储在长寿命逻辑编码基底中的最终逻辑Bell态. 为避免冗长表述, 后文在不引起歧义时仍将$ \left| {B} \right\rangle $简称为目标Bell型纠缠态.

    • 由于两个原子受到近似相同的激光驱动, 且体系初始态$ \left| {G} \right\rangle $属于交换对称子空间, 因此系统动力学主要限制在对称子空间内. 引入对称中间态

      则在强里德伯阻塞条件$ V\gg\varOmega_{\mathrm{eff}}(t), |\varDelta(t)| $下, 双激发态$ \left| {rr} \right\rangle $仅以高阶微扰方式参与动力学, 可近似忽略其布居. 与此同时, 在如下大失谐条件下,

      中间态$ \left| {e} \right\rangle $的占据始终保持较小, 因此可以对其进行绝热消除[3739]. 由此, 系统可进一步约化为$ \left| {G} \right\rangle $$ \left| {B} \right\rangle $构成的有效两能级模型. 其中$ \varGamma_e $表示中间态$ \left| {e} \right\rangle $的衰减率. 关于从(1)式到下述有效两能级哈密顿量的详细推导, 包括交换对称子空间中的基底选取与矩阵表示、大失谐条件下对中间态的绝热消除、AC Stark位移的显式近似表达、强里德伯阻塞条件下对双激发通道的进一步约化, 以及相应高阶修正项的量级分析, 见补充材料A (online).

      在忽略整体能量平移并将中间态绝热消除所产生的光移修正并入有效失谐后, 在基矢$ \{\left| {G} \right\rangle, \left| {B} \right\rangle \} $下, 有效哈密顿量可写为

      其中$ \varOmega_{\mathrm{eff}}(t) $为对称子空间中$ \left| {G} \right\rangle $$ \left| {B} \right\rangle $之间的集体有效Rabi频率, $ \varDelta(t) $为有效失谐量. 需要说明的是, 原始双光子激发模型中的双光子失谐记为$ \delta_{2\mathrm{ph}} $, 而在有效两能级模型中, 所有由双光子失谐及光移修正共同贡献的总失谐统一记为$ \varDelta(t) $, 以避免符号混淆.

      对于双光子梯形激发过程, 在大失谐条件下, 单原子基态$ \left| {g} \right\rangle $与里德伯态$ \left| {r} \right\rangle $之间的有效两光子耦合近似为

      其中$ \varOmega(t) $表示实验中实际施加到单原子上的有效两光子Rabi频率. 由于双原子对称激发存在集体增强效应, 从$ \left| {G} \right\rangle $$ \left| {B} \right\rangle $的有效耦合满足

      其来源在于对称驱动算符作用于集体基态时满足

      因此, $ \left| {G} \right\rangle \leftrightarrow\left| {B} \right\rangle $之间的矩阵元相较于单原子情形增强了$ \sqrt{2} $倍. 更详细推导见补充材料A (online).

      在上述近似下, 中间态绝热消除所产生的AC Stark位移及其他小的能级修正均并入有效失谐$ \varDelta(t) $中统一处理. 在保留主导二阶光移修正并将有限阻塞引起的附加修正记为$ \delta_{\rm blk}(t) $后, 总有效失谐可近似写为

      详细推导见补充材料A (online).

    • 为了实现体系从初态到目标态的快速高保真转移, 本文采用Lewis-Riesenfeld (LR)不变量逆向工程方法设计控制脉冲. 对于(6)式所对应的有效两能级哈密顿量, $ \varOmega_{\mathrm{eff}}(t) $$ \varDelta(t) $分别表示集体有效Rabi频率和总有效失谐.

      按照LR不变量理论, 引入满足不变量方程

      的Hermitian不变量算符. 对于两能级体系, 可将其参数写为

      其中μ为具有频率量纲的常数, $ \theta(t) $$ \beta(t) $为待设计的辅助参数.

      将(6)式和(12)式代入不变量条件(11)式, 可得到控制参数应满足的逆向工程关系:

      式中, 字符上方的圆点表示对时间的导数.

      (13) 式和(14) 式给出了有效哈密顿量参数与辅助变量$ \theta(t) $$ \beta(t) $之间的精确对应关系. 因此, 只要合理选取满足边界条件的$ \theta(t) $$ \beta(t) $, 就可以逆向构造出实现目标态转移所需的捷径控制脉冲.

      进一步地, 为保证体系从初始态$ |G\rangle $平滑演化到目标态$ |B\rangle $, 取辅助参数满足如下边界条件:

      其中$ t_{\mathrm{f}} $为总演化时间. 该边界条件使不变量在初末时刻的本征态方向分别与$|G\rangle $$|B\rangle $对应, 从而为沿不变量本征态实现目标态转移提供了条件.

      在本文所采用的共振捷径方案中, 为了使控制脉冲具有平滑、单峰且非负的形式, 选择常数辅助相位

      此时有

      代入(13)式和(14)式可得

      由此可见, 在该辅助相位选择下, 原本的一般逆向工程关系显著简化为共振条件下的STA脉冲设计问题. 这样不仅避免了额外的时变失谐扫描, 而且对于本文采用的对称单光子脉冲映射, 主导的差分AC Stark位移在一阶近似下自动抵消, 因此实验上无需实施复杂的时变失谐扫描, 而只需通过预先标定将裸双光子失谐调至最优共振工作点附近.

    • 本文采用五次多项式插值函数

      显然, 该函数在$ t=0 $$ t=t_{\mathrm{f }}$处满足所要求的边界条件, 并且在整个演化区间内连续光滑.

      对(19) 式求导, 可得

      因此, 根据(18)式, 有效两能级模型中的STA控制脉冲写为

      在前述有效基底$ \{\left| {G} \right\rangle, \left| {B} \right\rangle \} $下, 集体有效耦合与单原子有效两光子耦合满足

      这里$ \varOmega(t) $表示单原子有效两光子Rabi频率.

      由(21)式和(22)式, 可进一步得到实际施加在单原子上的控制脉冲为

      由(23) 式可见, 所设计的STA脉冲在初末时刻均为零, 并在中间演化阶段呈单峰平滑变化, 这种脉冲形状有利于实验实现. 与此同时, 由于整个控制过程在有效两能级模型层面保持$ \varDelta(t)=0 $的共振条件下, 在本文采用的对称双光子驱动条件下, 实验上无需额外引入时变失谐扫描, 只需通过预先标定将系统工作点调至有效共振附近, 从而进一步简化了参数控制复杂度.

      综上, 本文通过选取合适的辅助参数$ \theta(t) $和常数相位$ \beta(t)=-\text{π}/2 $, 构造了从初态$ |gg\rangle $到目标基态-里德伯态对称单激发Bell型纠缠态

      的解析STA控制脉冲. 后续数值模拟将基于(23)式和(24)式对该方案的快速性、高保真度以及鲁棒性进行验证.

    • 需要说明的是, 本文控制脉冲的解析设计基于(6)式所示的有效两能级模型完成. 为了更直观地展示目标Bell型纠缠态的布居转移以及双里德伯激发态$ \left| {rr} \right\rangle $的残余泄漏, 本文在主要动力学结果中进一步采用中间态绝热消除后得到的四维有效模型进行数值求解, 基底取为$ \{\left| {gg} \right\rangle, \left| {gr} \right\rangle, \left| {rg} \right\rangle, \left| {rr} \right\rangle \} $. 该模型并未重新引入中间态动力学, 而是在保留有效两光子耦合和双激发阻塞通道的前提下, 对物理布居演化给出更直观的描述.

      进一步地, 有效两能级模型的适用性来源于两个层次的近似条件: 其一, 在大失谐条件$ |\varDelta_{\mathrm{m}}|\gg \varOmega_1, \varOmega_2, \varGamma_e $下, 中间态$ \left| {e} \right\rangle $的瞬时占据始终较小, 因此可进行绝热消除; 其二, 在强里德伯阻塞条件$ V\gg \varOmega_{\mathrm{eff}}, |\varDelta| $下, 双激发态$ \left| {rr} \right\rangle $远离共振, 仅以高阶微扰方式参与演化. 基于此, 本文首先在低维有效模型中构造解析脉冲, 再在更接近实验观测的有效物理基底中展示主要动力学结果, 最后回到包含中间态与耗散过程的完整双原子三能级模型中进行严格数值检验.

      需要强调的是, 四维物理基底模型并非重新引入独立的中间态动力学, 而是对应于在大失谐条件下已经完成中间态消除后的有效物理描述. 该模型保留了Bell态制备中最关键的两个特征: 一是初态$ \left| {gg} \right\rangle $向对称单激发态$ \left| {B} \right\rangle $的相干转移, 二是有限阻塞强度下双激发态$ \left| {rr} \right\rangle $引起的残余泄漏. 因此, 它适合用于评估所设计控制脉冲在强阻塞实验条件下的实际性能. 为进一步验证绝热消除与有效两能级近似的可靠性, 本文首先在该四维有效物理基底中展示主要动力学与鲁棒性结果, 随后再引入包含中间态与耗散过程的完整双原子三能级模型, 对所设计脉冲进行更严格的数值检验.

    • 基于(19)式和(23)式, 本文对总演化时间$ T=1$ μs的情形进行了数值模拟. 由五次多项式轨迹可知, $ \dot{\theta}(t) $$ t=T/2 $附近达到峰值, 其最大值约为$ 15\text{π}/(8 T) $. 代入$ T=1$ μs后可得有效驱动峰值约为$ 2\text{π} \times0.94$ MHz, 从而相应的峰值有效两光子Rabi频率为

      这一结果说明, 本文方案仅需峰值单原子有效两光子Rabi频率处于亚MHz至MHz量级的较弱驱动, 即可在微秒时间尺度内完成高保真Bell型纠缠态制备, 这对于减小中间态散射、AC Stark频移以及紫外激光功率负担均具有积极意义[13,15].

      图2给出了控制脉冲和主要态布居的时间演化. 图2(a)中, 有效两光子Rabi频率$ \varOmega(t) $呈平滑单峰波形, 并在起止时刻严格回到零; 在有效两能级模型层面, 所设计的总有效失谐$ \varDelta(t) $在整个演化过程中取零, 对应共振驱动. 图2(b)表明, 在强里德伯阻塞条件$ V=2\text{π} \times120\text{ MHz}\gg \varOmega_{\max} $下, 系统布居从初始态$ \left| {G} \right\rangle $平滑转移至目标态$ \left| {B} \right\rangle $, 而双激发态$ \left| {rr} \right\rangle $的布居始终维持在极低水平, 表明有效两能级近似和阻塞条件均得到良好满足. 在无耗散纯态演化近似下, 为定量刻画目标态制备效果, 定义末态保真度为

      最终在$ t = T $时刻, 目标态保真度高于0.999. 需要指出的是, 在理想无耗散的有效两能级模型中, 本文所构造的STA脉冲满足边界条件$ \theta(0) = 0 $$ \theta(T) =\text{π} $, 因此末态接近单位保真度在一定程度上是逆向工程构造所保证的. 基于此, 本节中的高保真结果主要用于验证所设计解析脉冲与目标态转移路径的一致性, 而方案在实际条件下的优势还需结合完整模型及开放系统中的动力学表现进一步评估.

    • 上述结果是在强里德伯阻塞和大失谐绝热消除条件下得到的有效两能级模型基础上获得的. 为了进一步检验所设计解析STA脉冲在更真实物理情形下的可行性, 下面在完整双原子三能级模型中对该方案进行数值验证, 并进一步考虑中间态自发辐射和里德伯态有限寿命所带来的耗散效应.

      需要进一步说明的是, 本文在有效两能级模型中采用的共振条件$ \varDelta(t)=0 $, 并不意味着实验上必须对任意时刻的修正后双光子共振进行复杂的实时频率跟踪. 对于大失谐双光子梯形激发, 中间态绝热消除后, 总有效失谐可写为

      其中$ \delta_{2\mathrm{ph}} $为裸双光子失谐, $ \delta_{\mathrm{AC}}(t) $表示差分AC Stark位移, $ \delta_{\mathrm{h.o.}}(t) $表示更高阶修正项. 在大失谐近似下, 主导的差分AC Stark位移满足

      对于本文完整模型验证中所采用的对称单光子脉冲映射:

      主导的时变差分光移在一阶上自动抵消, 即$ \delta_{\mathrm{AC}}(t)\approx 0 $. 因此, 本文所取的共振条件更准确地说应该是: 在对称驱动和预先频率标定基础上实现的有效共振工作点, 而非要求实验过程中对任意时变修正失谐进行独立实时补偿. 结合图3可以看到, 本文参数范围内, 完整模型中的最优工作点位于$ \delta_{2\mathrm{ph}}/2\text{π}= 0\; \mathrm{MHz} $附近, 这说明修正后的有效共振工作点位于裸双光子共振附近, 因而在本文参数范围内, 该共振STA方案具有较明确的实验实现路径.

      在完整模型中, 每个原子均包含基态$ \left| {g} \right\rangle $、中间态$ \left| {e} \right\rangle $和里德伯态$ \left| {r} \right\rangle $, 系统总哈密顿量可写为

      其中$ \varOmega_{1}(t) $$ \varOmega_{2}(t) $分别为两束驱动场的单光子Rabi频率, $ \varDelta_{\mathrm{m}} $为中间态单光子失谐, $ \delta_{2\mathrm{ph}} $为双光子失谐, V表示双里德伯态之间的相互作用强度. 为与前文有效模型中的解析STA脉冲相对应, 本文采用与有效耦合一致的脉冲映射关系, 使完整模型中单原子两光子有效耦合满足:

      其中$ \theta(t) $仍取前文给出的边界条件多项式形式. 为使完整模型数值计算具有明确可复现性, 本文在数值验证中采用如下对称单光子脉冲映射:

      关于完整模型中单光子脉冲映射的具体形式及其参数估算, 详见补充材料B (online). 由于前文已取

      故有

      从而单原子有效两光子Rabi频率和完整模型中的对称单光子脉冲可分别写为

      以及

      由此可见, 所采用的单光子驱动脉冲在$ t=0 $$ t= T $处均自然回零, 且在整个演化区间内保持光滑. 这样, 在对称单激发子空间中, 由$ \left| {gg} \right\rangle $到目标态$ \left| {B} \right\rangle $

      的集体有效耦合即恢复为前文有效两能级模型中的$ \varOmega_{\mathrm{eff}}(t)=\dot{\theta}(t) $.

      对于本文采用的参数$ T=1$ μs, $ \varDelta_{\mathrm{m}}=2\text{π}\times 300$ MHz以及峰值有效两光子Rabi频率$ \varOmega_{\max}= 2\text{π}\times0.66$MHz, 在满足(32)式的诸多可能实现中, 本文数值验证采用一组对称且便于实验实现的单光子脉冲映射: 由(33)式可得峰值单光子Rabi频率为

      这一结果表明, 完整模型中所需的单光子驱动强度处于当前中性原子实验可达的几十MHz量级.

      进一步地, 在大失谐条件下, 单原子中间态的瞬时占据量级可近似估算为

      对于(33)式所示的对称脉冲映射, (40)式可进一步写为

      因此, 其峰值量级约为

      这一数量级与图4(b)中完整模型数值所得中间态总布居峰值处于$ 10^{-3} $量级的结果相符. 进一步地, 中间态散射概率可估算为

      最后一步代入了本文所取参数$ \varGamma_{\mathrm{e}}=2\text{π}\times0.18$ MHz和$ \varDelta_{\mathrm{m}}=2\text{π}\times300$ MHz. 该估计结果说明, 在本文参数区间内, 中间态散射是受控的. 这里的$ \left| {B} \right\rangle $仍表示由$ \left| {g} \right\rangle $$ \left| {r} \right\rangle $构成的基态-里德伯态对称单激发纠缠态.

      进一步地, 为考察耗散对方案性能的影响, 系统动力学由Lindblad主方程描述,

      其中, 塌缩算符$ L_k $用于描述系统中的耗散通道, 分别对应中间态$ \left| {e} \right\rangle $向基态$ \left| {g} \right\rangle $的自发辐射以及里德伯态$ \left| {r} \right\rangle $的有限寿命衰减. 例如, 可取$ L_{{\mathrm{e}}, j}=\sqrt{\varGamma_{\mathrm{e}}}\times \left| {g} \right\rangle _j\left( e \right)_j $$ L_{{\mathrm{r}}, j}=\sqrt{\varGamma_{\mathrm{r}}}\left| {g} \right\rangle _j\left( r \right)_j $, 其中$ j=1, 2 $. 这里为简化起见, 将里德伯态的有限寿命等效描述为回到基态$ \left| {g} \right\rangle $的有效衰减通道.

      数值计算中采用参数$ V=2\text{π} \times120\text{ MHz},\;\varDelta_{\mathrm{m}}= 2\text{π} \times300\text{ MHz}, $ 并取中间态衰减率和里德伯态衰减率分别为$ \varGamma_{\mathrm{e}}=2\text{π} \times0.18\text{ MHz}, \;\varGamma_{\mathrm{r}}=2\text{π} \times0.01$ MHz.

      在开放系统情形下, 目标态制备保真度定义为

      其中, $ {\boldsymbol{\rho}}(T) $为演化终止时刻 T 的系统密度矩阵.

      为确定完整模型中的最佳工作点, 本文对原始双光子失谐$ \delta_{2\mathrm{ph}} $进行了数值扫描, 结果如图3 所示. 可以看到, 在本文所考虑的参数范围内, 末态保真度在$ \delta_{2\mathrm{ph}}/2\text{π}=0\; \mathrm{MHz} $附近达到最大值, 因此后续完整模型与开放系统模拟中均取$ \delta_{2\mathrm{ph}}=0 $作为工作点. 该结果说明, 在本文所采用的参数条件下, 系统的最优工作区间位于双光子共振附近, 因此实验上可通过预先频率标定将系统调至该工作点, 而无需引入复杂的时变失谐补偿.

      图4给出了有效两能级模型、完整双原子三能级模型以及考虑耗散后的开放系统模型的数值比较结果. 由图4(a)可以看出, 在无耗散条件下, 完整模型与有效两能级模型所得目标Bell型纠缠态保真度演化几乎重合, 末态保真度可达0.9999以上; 在进一步考虑中间态衰减和里德伯态有限寿命后, 末态保真度仍保持在97%以上, 仅出现有限下降. 这表明前文基于有效模型构造的解析STA脉冲在更真实的完整物理模型中仍然有效, 并对耗散具有较好的鲁棒性.

      图4(b)进一步给出了完整模型中主要布居的演化. 可以看到, 目标Bell型纠缠态布居$ P_{\mathrm{B}} $在演化结束时接近于1, 而双里德伯激发态布居$ P_{{\mathrm{rr}}} $在整个过程中始终维持在$ 10^{-5} $量级, 中间态总布居$ P_{\mathrm{e }}$的最大值仅约为$ 10^{-3} $量级. 这说明双激发泄漏和中间态真实占据均受到有效抑制, 从而验证了强阻塞近似和大失谐绝热消除在本文参数区间内的合理性. 综合来看, 本文提出的解析STA方案不仅在有效两能级模型中能够实现Bell型纠缠态的快速高保真制备, 而且在完整开放系统模型中仍保持较高的制备质量, 体现出良好的实验可行性. 需要进一步指出的是, 在理想无耗散条件下, STA方案的末态高保真在一定程度上由逆向工程边界条件所保证, 因此要更全面地评估其实际优势, 还需将其与传统绝热方案在开放系统中进行公平比较.

    • 为更有物理意义地评估本文方案的实际性能, 下面进一步将所提出的STA方案与传统绝热方案进行比较. 需要强调的是, 理想无耗散条件下 STA 方案的高保真转移主要用于验证逆向工程脉冲设计的正确性, 而方案相对于传统控制的实际优势, 更应通过考虑达到给定目标保真度所需的时间代价以及开放系统中的耗散表现来衡量.

      为避免仅在固定短时间条件下比较可能带来的偏向性, 本文首先对传统绝热方案在无耗散条件下实现目标Bell态制备所需的总演化时间$ T_{\mathrm{ad}} $进行了扫描, 结果如图5(a)所示. 比较中, 传统绝热方案仍采用与前文相同形式的失谐扫描脉冲, 即

      其中峰值有效两光子Rabi频率取为

      与本文STA方案的峰值有效两光子Rabi频率保持一致; 同时取

      可以看到, 随着总演化时间$ T_{\mathrm{ad}} $的增加, 传统绝热方案的末态保真度逐渐提高; 当目标保真度阈值取为$ F=0.99 $时, 所需总时间约为$ T_{\mathrm{ad}}=8.79\; \text{μs} $. 这一结果说明, 传统绝热方案并非不能实现较高保真度, 但需要显著长于本文STA方案的演化时间.

      在此基础上, 进一步采用Lindblad主方程将系统耗散纳入动力学模型, 并在相同参数约束下比较STA方案与传统绝热方案在开放系统中的表现. 具体而言, STA方案取总时间$ T=1.00$ μs, 传统绝热方案取其达到$ F=0.99 $所需的代表性时间$ T_{\mathrm{ad}}=8.79\; \text{μs} $, 其余系统参数、里德伯相互作用强度、单光子失谐以及驱动幅度约束均保持一致. 图5(b)给出了不同里德伯态衰减率$ \varGamma_r $下两种方案的末态保真度. 可以看到, 随着$ \varGamma_{\mathrm{r}} $的增大, 两种方案的保真度均有所下降, 但传统绝热方案下降更为迅速. 这是因为传统绝热过程持续时间更长, 系统在演化过程中累积了更多的耗散损失; 相比之下, STA方案能够在更短时间内完成目标态制备, 从而有效抑制退相干对末态保真度的影响.

      上述结果表明, 本文方案的实际优势并不主要体现在理想闭系统中由边界条件保证的“构造性”高保真, 而是体现在存在耗散时仍能保持较高的目标态制备效率. 与传统绝热方案相比, STA方案通过显著缩短演化时间, 有效降低了开放系统中的退相干累积, 因此更具有实际应用意义. 基于此, 下面进一步考察该方案在实验参数波动和噪声扰动下的鲁棒性.

    • 在实际实验中, 原子热运动、光镊位置抖动以及激光强度和相位噪声都会影响控制精度. 为评估本文方案在非理想条件下的稳定性, 定义误差模型为

      其中$ \delta_{\mathrm{V}} $描述由原子间距波动引起的里德伯相互作用相对误差, $ \delta_\varOmega $表示驱动幅度误差, $ \delta_\varDelta $表示附加失谐误差. 多普勒频移可视为附加失谐的一种典型来源, 而激光线宽噪声则通过相位扩散过程进入动力学模型. 具体而言, 将驱动场写为$ \varOmega(t){\mathrm{e}}^{\mathrm{i}\phi(t)} $, 其中随机相位$ \phi(t) $满足相位扩散近似:

      式中的D与激光线宽对应. 这里有限线宽的影响作为独立相干噪声源单独评估, 因此采用随机相位轨迹下的纯态平均描述, 而未与前述主方程耗散过程同时叠加. 这种处理有助于分离激光相位噪声对相干控制性能的独立影响, 避免不同误差源叠加后难以辨别其主导贡献. 对于每一条随机轨迹, 系统仍按照含时薛定谔方程进行纯态演化, 并将末态相对于目标Bell型纠缠态的投影定义为该轨迹对应的保真度. 最终, 通过对所有随机轨迹的保真度取平均, 得到有限线宽条件下的平均末态保真度. 在数值计算中, 对于每一个线宽参数点, 均采样1000条随机相位轨迹. 图中的误差条表示平均值的标准误差. 进一步的轨迹数收敛性检验表明, 本文所采用的采样数已足以保证统计结果收敛, 相关结果见补充材料C (online).

      图6展示了不同实验误差下的数值结果. 图6(a)表明, 在较宽的相互作用误差范围内, 末态保真度始终保持接近于1. 这是因为本文方案工作在深度里德伯阻塞区$ V\gg \varOmega_{\mathrm{eff}} $, 适度的相互作用波动不会显著改变有效两态子空间中的动力学轨迹. 图6(b)给出了多普勒失谐对保真度的影响. 可以看到, 在50 kHz量级的多普勒失谐范围内, 末态保真度仍保持在约0.995附近, 说明所设计脉冲对热运动引起的附加失谐具有较好的容忍度. 图6(c)给出了激光线宽噪声下经Monte Carlo平均得到的末态保真度. 随着激光线宽增大, 平均保真度逐渐下降, 说明相位扩散噪声会削弱有限时间相干控制过程; 尽管如此, 在几十kHz量级的激光线宽范围内, 该方案仍保持较高平均保真度, 表明其对有限激光相位噪声具有一定容忍度.

      为了进一步评估方案对驱动误差与失谐误差的综合容忍度, 图7给出了末态保真度F随相对Rabi振幅误差$ \delta_{\varOmega} $和归一化失谐差$ \delta_{\varDelta}/\varOmega_{\max} $变化的二维分布, 其中$ \varOmega_{\max} $表示本文STA方案的峰值有效两光子Rabi频率. 与采用传统绝热方案参考频率进行归一化的方式相比, 这里采用自身的特征耦合尺度作为归一化基准, 因而更能直接反映所提STA控制协议对静态失谐偏差的固有容忍度. 图中横坐标采用百分比形式表示相对Rabi振幅误差, 纵坐标表示以$ \varOmega_{\max} $归一化后的附加失谐误差. 可以看到, 高保真区域围绕零误差点分布, 并在失谐方向上保持较明显的延展, 表明本文设计的STA脉冲对静态失谐差具有较强鲁棒性; 同时, 在约$ \pm 10{\text{%}} $的Rabi振幅误差范围内, 末态保真度仍可保持在较高水平. 图中的白色虚线分别对应于$ F=0.990 $$ F=0.999 $的等值线, 它们进一步刻画了高保真工作区域的范围.

      除了全局振幅误差和静态失谐偏差外, 实际实验中两个原子所受激光驱动也可能并不完全相同, 例如由于聚焦光斑不均匀、光镊位置偏差或局域光强漂移等因素, 会导致两个原子的有效耦合强度存在不对称性. 为此, 本文进一步引入双原子耦合不对称模型:

      其中$ \epsilon $描述两原子间的相对耦合不对称程度. 图8(a)给出了目标Bell型纠缠态末态保真度以及反对称单激发态泄漏随$ \epsilon $的变化. 可以看到, 在适度耦合不对称范围内, 目标态保真度仍保持较高水平, 但随着不对称程度增大, 系统会逐渐产生反对称态

      的泄漏. 这是由于理想情况下的交换对称驱动被破坏后, 体系不再严格限制在对称单激发子空间内. 尽管如此, 在实验上可接受的不对称范围内, 该方案仍表现出较好的鲁棒性.

      进一步地, 考虑到光镊中原子的残余热运动会引起原子间距离随时间微小变化, 从而导致里德伯相互作用在单次演化过程中呈现弱时变特征, 本文采用如下模型进行描述:

      其中η为相互作用调制深度, ν为调制频率, ϕ为初始相位. 图8(b)给出了不同代表性调制频率下末态保真度随η的变化关系. 结果表明, 在本文考虑的微弱扰动范围内, 只要位置波动所对应的相互作用调制深度保持在较小范围内, 时变里德伯相互作用对目标Bell型纠缠态制备结果保真度的影响仍较有限. 这说明本文方案不仅对静态相互作用误差具有鲁棒性, 对更贴近实验实际的弱时变相互作用扰动同样具有较好的容忍能力.

    • 为了将上述理论控制方案应用于具体实验平台, 下面对171Yb体系中的能级选择与激发方式进行讨论. 本文选取171Yb的$ ^{1}{\mathrm{S}}_0 $基态核自旋态作为量子比特编码, 即

      由于该编码主要依赖核自旋自由度, 因此对外部磁场噪声具有较强抑制作用, 同时避免了碱金属原子中较复杂的超精细结构耦合问题[27,28,3234].

      在里德伯激发方面, 本文考虑将高主量子数三重态$ n^{3}{\mathrm{S}}_1 $作为目标里德伯态, 例如$ n\approx70 $时可提供足够强的范德瓦耳斯相互作用, 从而满足强阻塞条件. 实验上可采用双光子梯形激发路径, 经由窄线宽中间态$ ^{3}{\mathrm{P}}_1 $实现从$ {}^1{\rm S}_0$$ n^{3}{\mathrm{S}}_1 $的耦合. 一种可行的波长组合为: 第一步使用556 nm激光驱动$ ^{1}{\rm S}_0 \leftrightarrow$${}^{3}{\rm P}_1 $跃迁, 第二步使用紫外激光驱动$ ^{3}{\mathrm{P}}_1\leftrightarrow n^{3}{\mathrm{S}}_1 $跃迁[2527,3134]. 这一方案与现有171Yb光镊和窄线宽激光技术的发展方向具有较好的兼容性.

      为避免理论符号化简掩盖具体实验实现路径, 这里给出一条与本文模型相一致的可行子能级示例. 可取参与里德伯激发的有效基态为

      例如, 在弱偏置磁场下, 采用π偏振先耦合到

      上述路径满足电偶极跃迁选择定则$ \Delta m_F=0 $, 同时可借助Zeeman分辨与频率选择, 使另一核自旋编码态在演化过程中作为旁观态保留. 需要强调的是, 本文理论分析并不依赖某一唯一的子能级实现方式, 而只要求实验上存在一条可通过偏振与频率选择实现的有效双光子激发路径. 因而, 文中所记的$ \left| {g} \right\rangle $可理解为在实验中经过态选择后参与里德伯耦合的有效基态.

      需要指出的是, 本文采用的参数区间与现有171Yb中性原子实验平台具有较好的兼容性[3133,3538]. 一方面, 对于主量子数$ n\sim 70 $附近的171Yb里德伯态, 在典型微米量级原子间距下可实现百 MHz 量级的有效相互作用强度, 从而满足强阻塞条件$ V\gg \varOmega_{\mathrm{eff}} $. 另一方面, 本文方案所需峰值有效两光子Rabi频率仅约为$ 2\text{π} \times0.66$ MHz, 处于亚MHz至MHz量级, 显著低于许多直接快速驱动方案对大功率激光的需求, 有利于减小激光功率负担及附加AC Stark频移. 与此同时, 由于控制过程中采用共振驱动而非频率扫描, 实验上只需通过声光调制器或电光调制器对激光振幅进行整形即可实现所设计的脉冲波形, 从而降低了控制复杂度. 综合来看, 本文提出的控制协议在当前实验参数范围内具有一定实现可行性.

    • 在实验实现上, 本文使用的单原子有效两光子 Rabi频率满足

      相应地, 在对称单激发子空间中, 集体有效耦合满足$ \varOmega_{\mathrm{eff}}(t)=\sqrt{2}\varOmega(t) $. 关于完整模型中对称单光子脉冲的映射方式及相应实验参数的量级估算, 详见补充材料B (online).

      因此, 只要适当增大中间态失谐$ \varDelta_{\mathrm{m }}$, 即可在抑制中间态布居和自发辐射散射的同时获得所需的有效耦合. 根据本文数值结果, 所需峰值有效两光子Rabi频率仅约为$ 2\text{π} \times0.66$ MHz, 明显低于许多依赖强驱动的快速量子门方案. 这意味着实验上无需采用极高功率的紫外激光场即可完成目标操作, 从而有助于降低功率噪声和附加AC Stark频移带来的误差.

      对于本文在完整模型数值中采用的对称单光子脉冲映射, 峰值单光子Rabi频率约为

      其中取$ \varDelta_{\mathrm{m}}=2\text{π}\times300$ MHz. 在该参数下, 中间态峰值占据约为

      而中间态散射概率估计约为

      说明大失谐条件下中间态真实占据和散射均受到有效抑制.

      进一步地, $ 556\text{ nm} $窄线跃迁$ ^{1}{\mathrm{S}}_0\leftrightarrow {}^{3}{\mathrm{P}}_1 $所需功率可以作简单量级估算. 该跃迁的自然线宽为$ \varGamma/2\text{π} \approx 182\; \mathrm{kHz} $, 饱和强度$ I_{\rm sat} \approx 0.14\, \rm{mW/cm^2}$[55,56]. 若取典型聚焦束腰$ w_{556}\sim 3\;\text{μm} $, 则由饱和强度估算可知, 达到(58)式所示单光子Rabi频率所需功率处于微瓦量级; 即使考虑Clebsch-Gordan系数、偏振失配及实验裕量, 其需求通常仍可控制在数$ \text{ μW} $量级. 对于第二步紫外激光驱动$ ^{3}{\mathrm{P}}_1\leftrightarrow n^3 {\mathrm{S}}_1 $跃迁, 在典型$ w_{\mathrm{UV}}\sim 3 $$ 5\;\text{μm} $的聚焦条件下, 按典型偶极矩量级估算, 实现$ \varOmega_{2, \max}/2\text{π}\approx 19.9$ MHz所需UV功率可处于数十毫瓦量级. 这里给出的功率需求均应理解为实验量级估算, 而非唯一确定值; 具体数值会随主量子数n、偏振以及跃迁偶极矩而变化[57]. 上述参数量级与现有171Yb光镊平台和类碱土里德伯光镊实验条件相兼容[3538].

      另一方面, 本文取里德伯相互作用强度$ V= 2\text{π}\times120$ MHz. 若选取典型光镊间距

      则由$ V=C_6/R^6 $[13]可对所需相互作用系数量级作出估算. 这里将V以角频率形式表示, 因此$ V/2\text{π} $对应普通频率单位, 从而有

      该结果表明, 本文所取参数对应的$ C_6/h $处于约$ 10^2 $$ 10^3\; \mathrm{GHz}\cdot {\text{μ}}\mathrm{m}^6 $量级, 与文献中高主量子数里德伯态可达到的强相互作用范围相符, 因此本文所采用的V值具有合理的物理基础[58].

      从时间尺度上看, 本文提出的控制过程总演化时间约为1 μs. 对于$ n\approx70 $171Yb里德伯态, 其自然寿命修正后的有效寿命通常仍明显长于该操作时间[13,19,45,59]. 因此, 本文方案有助于在退相干效应充分积累之前完成纠缠态制备. 此外, 由于本文采用对称双光子驱动并在实验前对双光子共振点进行标定, 整个控制过程可在有效共振工作点附近运行, 实验中无需进行复杂的频率扫描, 仅需借助声光调制器或电光调制器对激光振幅进行波形整形即可, 从而有望简化实验控制要求.

      关于本文共振 STA方案对实验设置的基本要求, 还可进一步概括如下: 其一, 需要满足大失谐条件

      以保证中间态可绝热消除, 并抑制真实占据和散射; 其二, 实验上应尽量采用近似对称的双光子驱动脉冲, 以减小差分AC Stark位移对有效双光子共振条件的破坏; 其三, 在正式施加控制脉冲前, 可通过双光子谱线测量或Ramsey标定将裸双光子失谐$ \delta_{2\mathrm{ph}} $调至最优工作点附近; 其四, 只要残余有效失谐保持在小于有效耦合尺度的范围内, 方案仍可维持较高末态保真度. 结合本文对失谐偏差、多普勒频移和激光线宽的数值分析可见, 在典型实验参数范围内, 该共振STA方案对小幅残余失谐仍具有一定容忍度.

    • 除了里德伯态寿命和激光控制精度外, 实验中还会受到原子残余热运动、多普勒频移及光镊位置波动等因素的影响. 本文前述鲁棒性分析表明, 在典型冷原子实验参数范围内, 这些非理想因素不会显著降低目标态保真度. 特别是在深度里德伯阻塞条件下, 适度的原子间距波动仅会引起有限的相互作用涨落, 而不会破坏整体控制路径.

      此外, 需要进一步说明本文所制备纠缠态在量子信息处理中的定位. 本文直接得到的是基态–里德伯态对称单激发Bell型纠缠态, 它可作为后续逻辑纠缠制备的中间资源态. 若实验目标是获得完全存储于长寿命核自旋基态中的逻辑Bell态, 则可在本文控制协议结束后, 利用附加的相干退激发或受控态映射脉冲, 将$ \left| {r} \right\rangle $选择性转移到另一稳定基态(例如另一核自旋编码态或辅助长寿命基态), 从而实现从$ \left(\left| {gr} \right\rangle +\left| {rg} \right\rangle \right)/\sqrt{2} $到逻辑Bell态的转换. 由于本文方案已经在微秒时间尺度内完成主要纠缠资源态的快速制备, 因此可为后续在较低退相干累积条件下实施相干映射留下时间窗口.

    • 本文提出了一种基于Lewis-Riesenfeld不变量的绝热捷径控制方案, 用于171Yb中性原子里德伯阻塞体系中的快速纠缠态制备. 本文直接制备的目标态为基态$ \left| {g} \right\rangle $与里德伯态$ \left| {r} \right\rangle $构成的对称单激发Bell型纠缠态, 而非已经完全存储在长寿命核自旋逻辑基底中的最终逻辑Bell态.

      数值结果表明, 在总演化时间$ T=1$ μs、里德伯相互作用强度$ V=2\text{π}\times120$ MHz的参数下, 本文方案仅需峰值单原子有效两光子Rabi频率约$ 2\text{π}\times0.66$ MHz, 即可实现保真度高于0.999的基态-里德伯态对称单激发Bell型纠缠态制备. 进一步的开放系统比较表明, 与传统绝热方案相比, 所提STA协议可在显著更短的演化时间内完成目标态转移, 从而减少耗散累积, 并保持更高的末态保真度.

      从实验实现角度看, 本文方案具有两个突出特点: 其一, 仅需较低的峰值有效两光子Rabi频率, 处于MHz量级, 有助于降低中间态散射AC Stark频移带来的误差; 其二, 整个控制过程在有效共振工作点附近运行, 无需复杂的时变失谐扫描, 因此在171Yb光镊阵列平台上具有较好的实验可行性. 未来可进一步结合开放系统主方程与原子运动效应, 并将该方案推广到多原子并行纠缠制备情形, 从而拓展其在中性原子量子信息处理中的应用. 若进一步结合实验中的相干态映射步骤, 本文所制备的基态-里德伯态Bell型纠缠态可作为实现长寿命逻辑Bell态制备的前置纠缠资源.

    参考文献 (59)

目录

/

返回文章
返回