表面电极离子阱中五线阱的输运研究.pdf

Modeling and Simulation 建模与仿真建模与仿真,2023,12(5),4707-4715 Published Online September 2023 in Hans.https:/www.hanspub.org/journal/mos https:/doi.org/10.12677/mos.2023.125429 文章引用文章引用:胡天杨,王旭.表面电极离子阱中五线阱的输运研究J.建模与仿真,2023,12(5):4707-4715.DOI:10.12677/mos.2023.125429 表面电极离子阱中五线阱的表面电极离子阱中五线阱的 输运研究输运研究 胡天杨，王胡天杨，王 旭旭*贵州大学大数据与信息工程学院，贵阳 贵州 收稿日期：2023年8月21日；录用日期：2023年9月13日；发布日期：2023年9月20日 摘摘 要要 量子计算机的扩展需要将多个子系统连接合成为可扩展的物理系统，离子作为量子比特载体，在各个子量子计算机的扩展需要将多个子系统连接合成为可扩展的物理系统，离子作为量子比特载体，在各个子系统间传递量子信息，所以离子输运是实现在多个离子囚禁区域或多个子系统间的量子比特扩展方案中系统间传递量子信息，所以离子输运是实现在多个离子囚禁区域或多个子系统间的量子比特扩展方案中必要的操控手段，可见离子输运的重要性极高。因此，我们整理了一种用于计算离子输运过程中分段直必要的操控手段，可见离子输运的重要性极高。因此，我们整理了一种用于计算离子输运过程中分段直流电极的时变电压的方法，在方法的设计中，我们不是从纯理论的角度来研究离子输运，还考虑了电子流电极的时变电压的方法，在方法的设计中，我们不是从纯理论的角度来研究离子输运，还考虑了电子学的实际约束，且让实验方法更加简洁明了。最终的实验结果表明该方法可以让离子按预期的路线输运，学的实际约束，且让实验方法更加简洁明了。最终的实验结果表明该方法可以让离子按预期的路线输运，这说明该方法是可行的，产生的直流电极电压是可靠的。这说明该方法是可行的，产生的直流电极电压是可靠的。关键词关键词 量子计算机量子计算机，离子囚禁离子囚禁，离子输运离子输运，直流电极直流电极 Study on Shuttling of Five-Wire Trap in Surface Electrode Ion Trap Tianyang Hu,Xu Wang*College of Big Data and Information Engineering,Guizhou University,Guiyang Guizhou Received:Aug.21st,2023;accepted:Sep.13th,2023;published:Sep.20th,2023 Abstract The expansion of quantum computers requires the synthesis of multiple subsystems into scalable *通讯作者。胡天杨，王旭 DOI:10.12677/mos.2023.125429 4708 建模与仿真 physical systems,and ions as qubit carriers transfer quantum information between subsystems.Therefore,ion shuttling is a necessary control means to realize the qubit expansion scheme in multiple ion confinement regions or between multiple subsystems,which shows the importance of ion shuttling.Therefore,we have formulated a method for calculating the time-varying voltage of the segmented DC electrode during ion shuttling.In the design of the method,we do not study ion shuttling from a purely theoretical point of view but also take into account the practical con-straints of electronics,and make the experimental method more concise and clear.The experi-mental results show that the method can shuttle ions in the expected route,which shows that the method is feasible and the generated DC electrode voltage is reliable.Keywords Quantum Computer,Ion Confinement,Ion Shuttling,DC Electrode Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 实现可扩展量子计算机的主要候选系统之一是离子阱系统1 2，离子阱是一种通过在电极上施加交变电场实现带电原子囚禁的装置3，广泛应用于质谱分析4 5、精密光谱6、量子模拟7 8和量子计算9 10等领域。在量子计算领域，离子阱也是一种非常重要的物理体系，量子比特具有长相干时间11 12、高保真度13、高保真度量子逻辑门14、比特间全连通15等多个优点，也在少量子比特数的情况下实现了量子纠错实验16 17。当前，制约离子阱系统发展的关键因素是在保持量子态及量子操作高保真度的前提下，如何提高系统内量子比特数目18的方法，即系统的可扩展性19。一种可能的方法是通过设计预实现基于 QCCD 架构的离子阱芯片20，通过表面电极对离子囚禁与不同功能区域进行定义和划分、结合离子在不同功能区间的输运，最终提高全局高保真度量子21比特数量。离子的输运是实现在多个囚禁区域的量子比特扩展方案中必要的操控手段。所谓的多个囚禁区域是指不同区域有着不同的功能，如载入区、操控区、冷却区。载入区载入离子时污染大，则需要运输到洁净的操控区进行操控，操控结束输运到冷却区保证离子在运动的基态。那么离子在不同功能区之间的输运起到一定现实的重要意义。目前，已有很多实验组成功实现了离子的输运，但是他们对 DC 电极电压波形的设计过程过于复杂，不易于理解，在这里，整理出了一种相较于前人方法更简单、更易于理解的一种方法，并在设计过程中考虑了电子学的实际约束，让该方法不仅仅是只适用于理论上。本文的其他安排如下：第 2 节叙述了离子阱结构，以及相应 RF 电压下的赝势情况以及管道高度。第 3 节详细讲解了如何设计 DC 电极电压波形以及其中的约束条件。第 4 节展示了用第 3 节的方法设计出了 DC 电极电压波形以及相应的结果图。最后对全文进行了总结。2.表面电极离子阱结构与表面电极离子阱结构与 RF 电极电极 本文设计了一种 DC 分段电极的表面电极离子阱，也称位五线阱，其结构如图 1 所示：Open AccessOpen Access胡天杨，王旭 DOI:10.12677/mos.2023.125429 4709 建模与仿真 Figure 1.Schematic diagram of surface electrode ion trap structure 图图 1.表面电极离子阱结构示意图 其中红色的电极为 RF 电极，蓝色的电极为 DC 电极，在以下叙述的输运过程中，只用到了图 1 中编好号的十对电极，其余 DC 电极接地或设置为 0 V。在电压施加过程中，电极编号相同的电极应施加相同的直流电压，以达到离子在 y 方向的位置不变或保持在 y=0 处。同理，红色的 RF 电极也应施加相同的RF 电压。理论上，DC 电极的宽度越窄越易于离子的控制，但是由于工程实践的可行性，将 DC 电极宽度设置为 70 um，既能在工程中实现，也能较好的控制离子的输运。DC 电极的长度理论上是越长越好，但是考虑到离子阱芯片的微型化，将长度设置为 500 um，既不影响实验的效果，也考虑到了离子阱芯片的构造。在表面电极离子阱中，离子受到电场束缚被囚禁在势阱中。由于射频电场是时变的，只能观察到某一时刻上的电势分布情况。对于这种时变的电场，赝势近似方法22是一种行之有效的射频电场分析方法。赝势近似可以将离子的电场和电荷分别用电势和电荷来表示，并在计算中将离子的电荷与表面电场相互作用的“真实”势能替换为一个简化的赝势等效势能。这样可以极大地简化计算，并提高计算效率。赝势的计算公式为()()22,4PPRFRFex y zx y zm=(1)其中，e 为离子的电荷量，m 为离子的质量，RF为射频驱动频率，赝势近似将电压值为RFV的直流电压产生的空间静电势(),RFx y z转化为射频电压()cosRFRFVt产生的等效势。当在图一所示的 RF 电极施加 200 V 电压，射频驱动频率为 10 Mhz 时：选取 X=0 处的 yz 面看赝势情况，如图 2 所示：(a)(b)Figure 2.Radial pseudopotential distribution generated by ion traps at surface electrodes 图图 2.表面电极离子阱产生的径向赝势分布 70um60um间隙=7um100um500um1 2 3 4 5 6 7 8 9 101 2 3 4 5 6 7 8109xzy胡天杨，王旭 DOI:10.12677/mos.2023.125429 4710 建模与仿真 因为芯片阱是关于坐标原点对称的，所以在 x=0，y=0 这个位置处看赝势最低点处于的高度，如图2(b)即离子最终被囚禁住的高度。由图 2 可知，在 RF 电极电压为 200 V，射频驱动频率为 10 MHZ 时，在芯片阱上方 72 um 处，y=0 um 处会有一条虚拟的赝势最低点的管道，即后面离子输运时的输运路径。3.DC 电极电压设计方法电极电压设计方法 本节的主要内容是讲在图 1 所示的表面电极离子阱，也称五线阱中离子的直线输运设计方法，本文求解电极电压波形的过程中，主要分为以下几个步骤：1)获得各个 DC 电极的“基函数”。2)确定空间中期望的电势。3)确定求解过程中的约束条件。4)反解出 DC 电压波形。下面就以上四个步骤进行展开讨论。3.1.获得各个获得各个 DC 电极的基函数电极的基函数 所谓“基函数”，是指陷阱阵列中每个 DC 电极的单位电压对空间中离子的贡献。获取 DC 电极基函数的具体方法如下：1)将所有的 RF 电极接地。2)设置第 j 个 DC 电极电压为 1 V，其余 DC 电压均为0 V。3)用物理仿真软件(cpo)计算第 j 个电极在空间中产生的电势 j，j则是第 j 个电极的“基函数”。4)对所有 DC 电极遍历第二、三步，保存各个电极的“基函数”，得出的 10 对电极的“基函数”图如图 3 所示：Figure 3.“Basis function”of 10 pairs of electrodes 图图 3.10 对电极的“基函数”因为要保证离子最终的输运管道处于 y=0 处，所以编号相同的电极必须加相同的电压。将编号相同的电极看作一对电极，图 3 为十对电极的“基函数”。有了 DC 电极的“基函数”，则可以计算出 DC电极在管道上任意一点产生的总电势，即：1nDCjjjU=(2)其中 Uj表示 DC 电极加载的实际电压，j表示该电极在该点处的电势贡献，DC即为该点处的总电势值。胡天杨，王旭 DOI:10.12677/mos.2023.125429 4711 建模与仿真 3.2.确定空间中期望的电势确定空间中期望的电势 在五线阱的模拟中，离子的输运可以假想为谐波电势最小值位置的移动。所以只需要移动谐波势的最小位置即可，在理想情况下，离子总是位于电势最小的位置。期望的电势 Vdes的具体表达式如下：()()()()22min1,2xdesmtVx txxtq=(3)其中 m 和 q 分别代表离子的质量和电荷量，x代表离子输运轴向的久期频率，Vdes表示空间中期望的电势，Xmin为电势最低点的位置。Figure 4.Desired potential at Xmin=0 图图 4.Xmin=0 时的期望电势 图 4 为当离子输运到 Xmin=0 时，期望的 DC 电极在管道上产生的电势如图所示，可以保证在 Xmin处的电势依然是附件电势的最小值，当改变 Xmin的值时，相应电势最小值的位置也会改变，也就意味着离子在输运过程中位置的改变。3.3.确定求解过程中的约束条件确定求解过程中的约束条件 通过 3.1 中提到的“基函数”以及公式(4)，可以计算出加了 DC 电极电压以后在管道上的实际电势情况，此时的总电势即为实际电势，即：()()(),gennnnVt xvtx=(4)Vn即为 n 个 DC 电极的电压，在本文中，n=10。同时，也定义了一个 D 函数来衡量实际产生的电势与理想情况下电势的差距，如公式(5)所示：()()()2,gendesD t xVt xVt x=(5)下面定义了一个高斯项 G，其作用是定义了每一个 Xmin点处，以 Xmin为中心，RX为宽度的感兴趣的区域。因为更在意电势最低点附近的区域，而不是所有空间位置。这样可以确保后面所说的损失函数(8)不受到远离 X=Xmin点的影响。根据公式(7)，可将其他参数化，使其宽度与选择的最小陷阱深度 VR一致。胡天杨，王旭 DOI:10.12677/mos.2023.125429 4712 建模与仿真 ()()()2min2,expxxxtG t xRt=(6)222xRxRV q m=(7)()()()()()0222201211dd,ddddfrlMNtxnnndeftxmnvtvtCG t x D t xxrvtvrrttt=+(8)最后，给出了总的损失函数公式(8)，其中 t0和 tf分别代表输运开始的时间点和结束的时间点，Xl和XR表示输运的起始位置和结束位置，m 表示输运中一共考虑的点数，n 表示电极数，r0，r1，r2分别表示权重系数，在模拟实验中，需要根据自己的需要调整权重系数，Vdef表示设置的默认 DC 电极电压，r0中这一项即表示实际输出的电压与设置的默认电压的差距，r1中的项和 r2中的项是波形的一阶和二阶导数，这用于减少其在高频下的光谱含量，并限制波形的转换速率。其作用是满足在硬件条件下实验的要求，使电压的变化率尽量小。3.4.求解求解 DC 电极电压波形电极电压波形 在求解电极电压波形的实验中，需要用到一种二次规划的数值算法对成本函数(8)进行最小化，该算法在几种方法中被发现是对解决传输问题最有用的23 24，如今，通常是用 Python 里的 MOSEK 求解器来求解。在实际求解电极电压时，需要在把输运管道分为间隔为 1 um 的点的离散网格上进行。4.实验结果与分析实验结果与分析 本文模拟了 Yb+离子从40 um 到 40 um 的输运情况，轴向囚禁频率设置为 500 KHz，其中权重系数r0，r1，r2分别设置为 5*1010，5*105，2*106，权重系数不固定，可根据自己的需要进行调整，默认电压 Vdef设置为 5 V，VR设置为 5。Figure 5.Shows the Vdes diagram during transport 图图 5.输运过程中 Vdes图 图 5 中只展示出了从输运开始到结束以 10 um 为间隔的期望电势结果图。胡天杨，王旭 DOI:10.12677/mos.2023.125429 4713 建模与仿真 Figure 6.Visualization of weight function(7)when Xmin=0 图图 6.Xmin=0 时权重函数(7)可视化 图 6 举例出了输运过程中的一个点(x=0)时的权重函数(7)，可以看到，当离子输运到 x=0 处时，最感兴趣的是 x=0 点，其权重系数为 1，而附件点的权重依次降低，其目的是保证成本函数(8)受到远离Xmin附件的点的影响更小。Figure 7.DC electrode voltage variation 图图 7.DC 电极电压变化 图 7 为输运过程中 10 对 DC 电极的电压变化，右上角的标号分别对应图 1 中的电极编号，可以看到，电极电压变化非常平滑，并且最大值均保持在10 V 以内，完全满足实验过程中的实验要求。在图 7 所示的电压下，对实际产生的空间电势 Vgen与想要空间电势的 Vdes的部分点进行了采样并比较，如图 8 所示：胡天杨，王旭 DOI:10.12677/mos.2023.125429 4714 建模与仿真 Figure 8.Sampling points at x=40,20,0,20,40um in the transport process 图图 8.对输运过程中 x=40、20、0、20、40 um 处的采样点 通过以上图 8 实际电势与理想电势的对比可知，理想电势与实际电势基本重合，尤其在电势等于 0的点附近，这足以说明离子会按照设计的管道路径进行输运。模拟实验证明，这种方法产生的电极电压波形对离子的输运是可靠的。5.结束语结束语 离子输运是构造可扩展的量子计算机和量子模拟器的必要元素。因此本文设计了一套完整的波形产生方法，与前人的设计方法相比25 26，此套设计方法更为简单，最后在本文设计的离子阱芯片上成功模拟实现了离子的输运，足以说明该套方法产生的电极电压波形是可靠的。参考文献参考文献 1 Wineland,D.J.,Monroe,C.,Itano,W.M.,Leibfried,D.,King,B.E.and Meekhof,D.M.(1998)Technology,Experi-mental Issues in Coherent Quantum-State Manipulation of Trapped Atomic Ions.The Journal of Research of the Na-tional Institute of Standards and Technology,103,259-328.https:/doi.org/10.6028/jres.103.019 2 Kaushal,V.,Lekitsch,B.,Stahl,A.,Hilder,J.,Pijn,D.,Schmiegelow,C.,Bermudez,A.,Mller,M.,Schmidt-Kaler,F.and Poschinger,U.(2020)Shuttling-Based Trapped-Ion Quantum Information Processing.AVS Quantum Science,2,Article ID:014101.https:/doi.org/10.1116/1.5126186 3 Siverns,J.D.and Quraishi,Q.(2017)Ion Trap Architectures and New Directions.Quantum Information Processing,16,Article No.314.https:/doi.org/10.1007/s11128-017-1760-2 4 Schwartz,J.C.,Senko,M.W.and Syka,J.(2002)A Two-Dimensional Quadrupole Ion Trap Mass Spectrometer.Jour-nal of the American Society for Mass Spectrometry,13,659-669.https:/doi.org/10.1016/S1044-0305(02)00384-7 5 Decker,T.K.,Zheng,Y.,Ruben,A.J.,Wang,X.,Lammert,S.A.,Austin,D.E.and Hawkins,A.R.(2019)A Microscale Planar Linear Ion Trap Mass Spectrometer.Journal of the American Society for Mass Spectrometry,30,482-488.https:/doi.org/10.1007/s13361-018-2104-x 6 Ozawa,A.,Davila-Rodriguez,J.,Bounds,J.R.,Schuessler,H.A.,Hnsch,T.W.and Udem,T.(2017)Single Ion Fluo-rescence Excited with a Single Mode of an UV Frequency Comb.Nature Communications,8,Article No.44.https:/doi.org/10.1038/s41467-017-00067-9 7 Blatt,R.and Roos,C.F.(2012)Quantum Simulations with Trapped Ions.Nature Physics,8,277-284.https:/doi.org/10.1038/nphys2252 8 Massar,S.,Spindel,P.,Varn,A.F.and Wunderlich,C.J.(2015)Investigating the Emergence of Time in Stationary 胡天杨，王旭 DOI:10.12677/mos.2023.125429 4715 建模与仿真 States with Trapped Ions.Physical Review A,92,Article ID:030102.https:/doi.org/10.1103/PhysRevA.92.030102 9 Nielsen,M.A.and Chuang,I.(2012)Quantum Computation and Quantum Information.Statistical Science,27,373-394.https:/doi.org/10.1214/11-STS378 10 Hffner,H.,Roos,C.F.and Blatt,R.(2008)Quantum Computing with Trapped Ions.Physics Reports,469,155-203.https:/doi.org/10.1016/j.physrep.2008.09.003 11 Wang,P.,Luan,C.-Y.,Qiao,M.,Um,M.,Zhang,J.,Wang,Y.,et al.(2021)Single Ion Qubit with Estimated Cohe-rence Time Exceeding One Hour.Nature Communications,12,Article No.233.https:/doi.org/10.1038/s41467-020-20330-w 12 Barends,R.,Kelly,J.,Megrant,A.,Veitia,A.,Sank,D.,Jeffrey,E.,et al.(2014)Superconducting Quantum Circuits at the Surface Code Threshold for Fault Tolerance.Nature,508,500-503.https:/doi.org/10.1038/nature13171 13 Mehta,K.K.,Bruzewicz,C.D.,McConnell,R.,Ram,R.J.,Sage,J.M.and Chiaverini,J.(2016)Integrated Optical Ad-dressing of an Ion Qubit.Nature Nanotechnology,11,1066-1070.https:/doi.org/10.1038/nnano.2016.139 14 Balance,C.,Harty,T.,Linke,N.,Sepiol,M.and Lucas,D.M.(2016)High-Fidelity Quantum Logic Gates Using Trapped-Ion Hyperfine Qubits.Physical Review Letters,117,Article ID:060504.https:/doi.org/10.1103/PhysRevLett.117.060504 15 Monz,T.,Schindler,P.,Barreiro,J.T.,Chwalla,M.,Nigg,D.,Coish,W.A.,et al.(2011)14-Qubit Entanglement:Cre-ation and Coherence.Physical Review Letters,106,Physical Review Letters,130506.https:/doi.org/10.1103/PhysRevLett.106.130506 16 Ryan-Anderson,C.,Bohnet,J.G.,Lee,K.,Gresh,D.,Hankin,A.,Gaebler,J.,et al.(2021)Realization of Real-Time Fault-Tolerant Quantum Error Correction.Physical Review X,11,Article ID:041058.https:/doi.org/10.1103/PhysRevX.11.041058 17 Schindler,P.,Barreiro,J.T.,Monz,T.,Nebendahl,V.,Nigg,D.,Chwalla,M.,Hennrich,M.and Blatt,R.J.S.(2011)Experimental Repetitive Quantum Error Correction.Science,332,1059-1061.https:/doi.org/10.1126/science.1203329 18 Wu,W.B.,Zhang,T.and Chen,P.(2021)Quantum Computing and Simulation with Trapped Ions:On the Path to the Future.Fundamental Research,1,213-216.19 Decaroli,C.,Matt,R.,Oswald,R.,Axline,C.J.,Ernzer,M.,Flannery,J.,Ragg,S.and Home,J.P.(2021)Technology,Design,Fabrication and Characterization of a Micro-Fabricated Stacked-Wafer Segmented Ion Trap with Two X-Junctions.Quantum Science and Technology,6,Article ID:044001.https:/doi.org/10.1088/2058-9565/ac07ee 20 Kielpinski,D.,Monroe,C.and Wineland,D.J.(2002)Architecture for a Large-Scale Ion-Trap Quantum Computer.Nature,417,709-711.https:/doi.org/10.1038/nature00784 21 Clark,C.R.,Tinkey,H.N.,Sawyer,B.C.,Meier,A.M.,Burkhardt,K.A.,Seck,C.M.,et al.(2021)High-Fidelity Bell-State Preparation with 40Ca+Optical Qubits.Physical Review Letters,127,Article ID:130505.https:/doi.org/10.1103/PhysRevLett.127.130505 22 Dehmelt,H.G.(1968)Radiofrequency Spectroscopy of Stored Ions I:Storage.In:Advances in Atomic and Molecular Physics,Elsevier,Amsterdam,53-72.https:/doi.org/10.1016/S0065-2199(08)60170-0 23 De Clercq,L.E.(2015)Transport Quantum Logic Gates for Trapped Ions.ETH,Zurich.https:/doi.org/10.1103/PhysRevLett.116.080502 24 Oswald,R.(2015)Velocity Control of Trapped Ions for Transport Quantum Logic Gates.M.Sc.Thesis,Swiss Federal Institute of Technology,Zurich.25 刘威.面向多量子比特的表面电极离子阱芯片D:博士学位论文.长沙:国防科学技术大学,2014.26 陈昆.芯片上的离子囚禁与输运D:硕士学位论文.长沙:国防科学技术大学,2013.