本申请涉及量子通信技术领域，特别是涉及一种基于超导量子计算的多比特操控系统及方法。
背景技术
超导量子计算是目前量子计算机的主流实现方式之一，它利用约瑟夫森结(josephsonjunction)电路的量子特性进行编码，并按照量子力学规律来操控这些量子比特(qubit)。在这个过程中，对量子比特的操控尤为关键，这是因为对量子比特的操控将直接影响到量子计算机的运行速度与准确率，因此，如何快速而准确的多比特门操控对于实现通用量子计算机具有重要意义。
目前已被证明，量子通用逻辑门组可以仅由一位量子门和经典两位门构成，现有的实现量子多位门的方式主要分为间接实现和直接实现。即实现任意多比特门的间接方式是将其分拆为多个单比特门和双比特门构成的量子线路。然而这将耗费更多的操控时间，由于量子单位门或/及两位门操作存在一定的错误率(保真度小于1)，这种间接实现多比特门的方式将带来更低的保真度，并花费更多的时间。因此，如何实现单次多比特门(single-shotmulti-qubitgate)操作对于量子计算具有十分重要的意义。
技术实现要素：
鉴于以上所述现有技术的相关问题，本申请的目的在于提供一种基于超导量子计算的多比特操控系统及方法，用于解决实现单次多比特门的操作。
为实现上述目的及其他相关目的，本申请的第一方面提供一种基于超导量子计算的多比特操控系统，包括：微波输入模块，超导量子模块，测量模块，以及反馈模块。其中，所述微波输入模块用于将初始量子态及初始微波序列以交频信号输出；所述超导量子模块用于依据接收的所述交频信号操控内置的多比特门并输出待测量子态；所述测量模块用于接收所述待测量子态并进行测量以获得当前多比特门的固有保真度；所述反馈模块用于依据所述固有保真度进行计算生成优化微波序列，并将所述优化微波序列反馈给所述微波输入模块进行再次操控以获得目标多比特门。
在本申请的第一方面的某些实施方式中，所述微波输入模块包括微波装置。
在本申请的第一方面的某些实施方式中，所述微波输入模块输出的交频信号是经拆分不同时间段的信号叠加处理的交频信号。
在本申请的第一方面的某些实施方式中，所述超导量子模块包括共面波导谐振腔，所述共面波导谐振腔包括多个串联结构的超导传输子(transmon)，每一所述超导传输子(transmon)由库珀对盒构成。
在本申请的第一方面的某些实施方式中，所述超导量子模块包括共面波导谐振器。
在本申请的第一方面的某些实施方式中，多个串联结构的超导传输子(transmon)通过总线与所述共面波导谐振腔相耦合。
在本申请的第一方面的某些实施方式中，所述测量模块包括：信号测量单元及保真度计算单元。其中，所述信号测量单元用于接收所述待测量子态；所述保真度计算单元用于通过预设的测量方式对所述待测量子态进行当前多比特门的固有保真度。
在本申请的第一方面的某些实施方式中，所述预设的测量方式为x-y-z方向的投影测量。
在本申请的第一方面的某些实施方式中，所述预设的测量方式为rb(randomizedbenchmarking)测量，断层扫描(tomography)测量，门集合断层扫描(gatesettomography)测量，或者误差估计(errormitigation)测量。
在本申请的第一方面的某些实施方式中，所述反馈模块包括：参数优化单元及反馈单元，其中，所述参数优化单元用于依据所述固有保真度对所述初始微波序列进行优化处理后输出优化微波序列；所述反馈单元用于将所述优化微波序列反馈给所述微波输入模块进行再次操控以获得目标多比特门。
在本申请的第一方面的某些实施方式中，所述反馈模块还包括保真度计算单元，用于通过预设的测量方式对所述待测量子态进行当前多比特门的固有保真度。
在本申请的第一方面的某些实施方式中，所述参数优化单元利用bfgs算法或globalsearch算法对所述初始微波序列进行优化处理。
在本申请的第一方面的某些实施方式中，所述基于超导量子计算的多比特操控系统还包括优化参数生成模块，用于发送预设的初始微波序列给所述微波输入模块，所述预设的初始微波序列是经预优化处理的。
本申请的第二方面提供一种基于超导量子计算的多比特操控方法，包括以下步骤：将初始量子态及初始微波序列以交频信号输出；依据接收的所述交频信号操控内置的多比特门并输出待测量子态；对所述待测量子态并进行测量以获得当前多比特门的固有保真度；以及依据所述固有保真度进行计算生成优化微波序列，并将所述优化微波序列反馈进行再次操控以获得目标多比特门。
在本申请的第二方面的某些实施方式中，所述交频信号为经拆分不同时间段的信号叠加处理的交频信号。
在本申请的第二方面的某些实施方式中，所述多比特门为由多个串联结构的超导传输子(transmon)通过总线与一共面波导谐振腔相耦合组成，每一所述超导传输子(transmon)由库珀对盒构成。
在本申请的第二方面的某些实施方式中，所述预设的测量方式为x-y-z方向的投影测量。
在本申请的第二方面的某些实施方式中，所述预设的测量方式为rb(randomizedbenchmarking)测量，断层扫描(tomography)测量，门集合断层扫描(gatesettomography)测量，或者误差估计(errormitigation)测量。
在本申请的第二方面的某些实施方式中，所述依据所述固有保真度进行计算生成优化微波序列的步骤包括依据所述固有保真度对所述初始微波序列进行优化处理后输出优化微波序列在本申请的第二方面的某些实施方式中，所述参数优化单元利用bfgs算法或globalsearch算法对所述初始微波序列进行优化处理。
在本申请的第二方面的某些实施方式中，所述方法还包括生成预设的初始微波序列的步骤，所述预设的初始微波序列是经预优化处理的。
如上所述，本申请基于超导量子计算的多比特操控系统及方法，使用共面波导谐振模块将多个超导传输子(transmon)耦合起来，并通过简单的微波元器件对这些量子比特进行操控，同时反馈数据并根据设计方案进行优化，获得最优控制(optimalcontrol)参数，实现单次多比特门操作。本申请提供的技术方案具有可拓展性高，保真度高，操作时间短的优点，应用场景普遍，且成本较低，对于超导量子计算具有很强的实用价值。
附图说明
图1显示为本申请的基于超导量子计算的多比特操控系统的原理框图。
图2显示为本申请的基于超导量子计算的多比特操控系统在一实施例中的示意图。
图3显示为本申请的基于超导量子计算的多比特操控系统在另一实施例中的示意图。
图4显示为本申请的基于超导量子计算的多比特操控系统在再一实施例中的示意图。
图5显示为本申请的基于超导量子计算的多比特操控方法在一实施例中的流程示意图。
图6显示为本申请的基于超导量子计算的多比特操控方法在另一实施例中的流程示意图。
具体实施方式
以下由特定的具体实施例说明本申请的实施方式，熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本申请的其他优点及功效。
在下述描述中，参考附图，附图描述了本申请的若干实施例。应当理解，还可使用其他实施例，并且可以在不背离本申请的精神和范围的情况下进行机械组成、结构、电气以及操作上的改变。下面的详细描述不应该被认为是限制性的，并且本申请的实施例的范围仅由公布的专利的权利要求书所限定。
这里使用的术语仅是为了描述特定实施例，而并非旨在限制本申请。空间相关的术语，例如“上”、“下”、“左”、“右”、“下面”、“下方”、“下部”、“上方”、“上部”等，可在文中使用以便于说明图中所示的一个元件或特征与另一元件或特征的关系。另外，虽然在一些实例中术语第一、第二等在本文中用来描述量子态，但是这些量子态不应当被这些术语限制。这些术语仅用来将一个量子态与另一个量子态进行区分。例如，第一测试量子态可以被称作第二测试量子态，并且类似地，第二测试量子态可以被称作第一测试量子态，而不脱离各种所描述的实施例的范围。第一测试量子态和测试量子态均是在描述一个测试量子态，但是除非上下文以其他方式明确指出，否则它们不是同一个测试量子态。相似的情况还包括第三测试量子态及第四测试量子态。
再者，如同在本文中所使用的，单数形式“一”、“一个”和“该”旨在也包括复数形式，除非上下文中有相反的指示。应当进一步理解，术语“包括”、“包括”表明存在所述的特征、步骤、操作、元件、组件、项目、种类、和/或组，但不排除一个或多个其他特征、步骤、操作、元件、组件、项目、种类、和/或组的存在、出现或添加。此处使用的术语“或”和“和/或”被解释为包括性的，或意味着任一个或任何组合。因此，“a、b或c”或者“a、b和/或c”意味着“以下任一个：a；b；c；a和b；a和c；b和c；a、b和c”。仅当元件、功能、步骤或操作的组合在某些方式下内在地互相排斥时，才会出现该定义的例外。
快速而准确的多比特门操控对于实现通用量子计算机具有重要意义，而实现多比特门操控的物理基础是制备出相互耦合的量子单位。在超导量子计算中，使用超导传输子(transmon)作为基本量子单位构建量子比特，超导传输子(transmon)也是目前超导体系中保真度最高、可调控性最高的物理单位之一。此外，如何将这些超导传输子(transmon)相互耦合起来也是一个重要的问题。
需要特别说明的是，由于transmon为在超导量子计算中构成量子比特的物理单位，在本申请中transmon被称之为超导传输子，在相关的技术文献中，transmon还被称之为超导量子比特，或者超导电荷，又或者超导电荷量子位，这些名称均可表示为本申请所述的超导传输子(transmon)。
通用量子计算的实现基于高保真度的量子门构成的量子线路，其中量子多位门具有很重要的作用。量子纠错码的实现需要作用在至少3个量子比特的多位门，例如toffoli门(toffoli，orcontrolled-controlled-notgate，又被称作ccnotgate或“控-控-非”门)，toffoli门在可逆量子计算中也同样有着十分重要的作用。现有的实现量子多位门的方式主要分为间接实现和直接实现。其中，间接实现的方式是将其分拆为多个量子单比特门和双比特门构成的量子线路实现(参阅doi:10.1038/nature10786)，然而这将耗费更多的时间，并带来更高的错误率；直接实现量子多位门也称为单次多比特门，现有的方案有：通过含时操控的微波信号来区分超导传输子(transmon)耦合系统的不同能级，实现多比特门(参阅doi:https://doi.org/10.1103/physreva.87.022309)；将微波信号转化为分段常数信号对多位量子比特进行操控(参阅doi:10.1103/physrevb.85.054504)。
目前已有一些关于多比特门实现方案的研究，然而高保真度的通用的量子多比特门操控尚未被有效实现。如前所述，通过将多比特门拆分为单比特门和双比特门组合的方式将耗费更多的操控时间，并带来更低的保真度，因而迫切需要寻求一种可以直接实现单次多比特门操控的方案。
关于单次多比特门的实现，诚如前述中提到的两种具有代表性的实现方案都有各自的缺陷。对于直接操控系统能级的方案，在transmon耦合体系较为复杂时，构建含时操控的微波信号难度将随之上升，且用这种能级转换的操控方式，在实现任意多比特门时具有较高的难度，可拓展性较差；对于优化分段常数振幅操控信号的方案，未能有效避免泄露误差(leakageerror)，将对多比特门的保真度带来一定影响。基于此，本申请提出一种基于超导量子计算的多比特操控系统，可用于在200ns、150ns或100ns等不同时间范围以内快速地单次对多比特门(single-shotmulti-qubitgate)的操控。
请参阅图1，显示为本申请的基于超导量子计算的多比特操控系统的原理框图，如图所示，所述基于超导量子计算的多比特操控系统包括：微波输入模块11，超导量子模块12，测量模块13，以及反馈模块14。所述微波输入模块11将初始量子态及初始微波序列以交频信号输出给所述超导量子模块12；所述超导量子模块12依据接收的所述交频信号操控内置的多比特门并输出待测量子态给所述测量模块13，所述测量模块13接收所述待测量子态并进行测量以获得当前多比特门的固有保真度；所述反馈模块14依据所述固有保真度进行计算生成优化微波序列，并将所述优化微波序列反馈给所述微波输入模块11进行再次操控以获得目标多比特门，经以上过程将经历多次重复，最后达到最优调制效果，进而实现单次多比特门操作。
请参阅图2，显示为本申请的基于超导量子计算的多比特操控系统在一实施例中的示意图，如图所示，所述微波输入模块11用于将初始量子态及初始微波序列以交频信号输出；在一实施例中，所述微波输入模块11包括微波装置，用于产生交频形式的微波信号，所述微波装置可以通过商用的微波元器件实现，例如为微波发生器。
在本实施例中，所述微波输入模块11输出的交频信号是经拆分不同时间段的信号叠加处理的交频信号，在实施例中，所述经拆分不同时间段的信号叠加处理的交频信号为不同时间段的drag信号的叠加，使得所述超导量子模块12中内置的所有比特门受同一交频信号的调制，避免了对不同比特门进行单独调制时的信道串扰。由于将信号依据不同的时间段分拆为的分段drag信号的组合，也很大程度的减少了泄露误差(leakageerror)。
所述超导量子模块12用于依据接收的所述交频信号操控内置的多比特门，输出待测量子态，所述待测量子态是经过对输入的初始初始量子态通过当前交频信号操控内置的多比特门后输出的另一种形式的量子态，在不同的实施方式中，所述待测量子态亦可被称之为测试量子态或制备量子态等等。所述待测量子态对应调制的目标多比特门ut，其中t表示当前时刻。在实施例中，所述超导量子模块12包括例如为共面波导谐振器(coplanarwaveguideresonator)的共面波导(coplanarwaveguide，cpw)谐振腔121，所述共面波导谐振腔121串联有多个串联结构的超导传输子(transmon)122，每一所述超导传输子(transmon)122由库珀(cooper)对盒构成，所述的多个串联结构的超导传输子(transmon)122通过总线与所述共面波导谐振腔121相耦合，即buscoupling方案。
在本申请提供的实施例中，所述的buscoupling方案是将制备好的超导传输子(transmon)122按列排布耦合到高品质因子的共面波导谐振腔121中，共面波导谐振腔121中的光子可以诱导超导传输子(transmon)122相互之间的远程强耦合，从而构建了超导量子计算机的物理底层模型。
所述测量模块13用于接收所述待测量子态并进行测量，测量的目的是为了获得当前多比特门的固有保真度，所述测量模块13为由多个测量设备和设备构成的信号测量端。在本实施例中，所述测量模块13包括：信号测量单元131及保真度计算单元132。
所述信号测量单元131用于接收所述待测量子态；所述保真度计算单元132用于通过预设的测量方式对所述待测量子态进行当前多比特门的固有保真度，即当前交频信号下的多比特门的固有保真度。在本实施例中，所述预设的测量方式为x方向，y方向，及z方向三个方向上的投影测量。
在某些实施例中，所述预设的测量方式可以为rb(randomizedbenchmarking)测量，断层扫描(tomography)测量(请参阅d'ariano，etal.physicalreviewletters86.19(2001):4195.)，门集合断层扫描(gatesettomography)测量(请参阅greenbaum,etal.arxivpreprintarxiv:1509.02921(2015).)，或者误差估计(errormitigation)测量(请参阅temme，kristanetal.physicalreviewletters119.18(2017):180509.)等，在本实施例中，优先选用rb(randomizedbenchmarking)测量，容后详述。
所述反馈模块14用于依据所述固有保真度进行计算生成优化微波序列，并将所述优化微波序列反馈给所述微波输入模块11进行再次操控以获得目标多比特门。在本实施例中，所述反馈模块14包括：参数优化单元141及反馈单元142。
所述参数优化单元141用于依据所述固有保真度对所述初始微波序列进行优化处理后输出优化微波序列；所述反馈单元142用于将所述优化微波序列反馈给所述微波输入模块11进行再次操控以获得目标多比特门。在一些实施例中，所述参数优化单元141利用bfgs算法(即broyden-fletcher-goldfarb-shanno，简称bfgs)或globalsearch算法对所述初始微波序列进行优化处理，其中，globalsearch算法通过matlab软件系统实现。在本实施例中，优选选用bfgs算法对所述初始微波序列进行优化处理。
在另一种实施方式中，所述反馈模块14还包括保真度计算单元143，即所述保真度计算单元143为反馈模块14的一部分。所述保真度计算单元143用于通过预设的测量方式对所述待测量子态进行当前多比特门的固有保真度，呈如图3所示，图3显示为本申请的基于超导量子计算的多比特操控系统在另一实施例中的示意图。
在具体的实现方式中，所述反馈模块14包含计算机设备中的软件和硬件，可例如装载有本申请所述反馈模块14功能或反馈模块14中各单元功能的计算机程序的计算机来实现，所述计算机包括但不限于服务器、台式电脑、或笔记本电脑等。
在实施例中，本申请将藉由多比特操控系统的哈密顿量来设计系统架构，即通过汉密顿量的演化来实现对系统的设计。所述哈密顿量为系统的能量算符，是一个描述系统总能量的算符，以h表示，则所述多比特操控系统的哈密顿量描述为：h＝h0+hd+hc。
其中，h0可以看作系统中超导传输子(transmon)本身的哈密顿量；
hd可以看作所述交频信号调制下的哈密顿量，
hc可以看做超导传输子(transmon)之间相互耦合的哈密顿量，在本申请所提供的实施例中，所述超导传输子(transmon)之间相互耦合的哈密顿量hc是也固定的，在本申请的一实施例中，采用x-y型耦合，即，对耦合哈密顿量近似到二阶，仅保留x-y之间的耦合，忽略z-z之间的耦合；
则
其中，jij表示为第i个和第j个超导传输子之间的耦合强度。
在本申请中，通过所述交频信号来调节或调制所述哈密顿量hd，即所述交频信号
本申请提供的实施方式中采用drag最优控制法，关于drag最优控制法的描述可以参阅doi:10.1103/physreva.82.040305(optimizeddrivingofsuperconductingartificialatomsforimprovedsingle-qubitgates)。在本申请提供的实施例中，例如在200ns或者150ns内时间范围内实现，通过drag控制法将所述交频信号比如以10ns或20ns为单位拆分成信号的叠加形式，以确定交频信号的频率。
所述drag最优控制法表示为
在具体的实现过程中，所述不同的时间段可以依据具体的实施情况进行设定，比如为在150ns的范围内设定，通过drag控制法将所述交频信号比如以10ns为单位拆分成信号的叠加形式，则所述微波输入模块11将随机生成或者接收的外部优化参数设备生成的15组参数序列{(an,bn)}1,{(an,bn)}2,…,{(an,bn)}15，即n＝1，2，……15。但并不局限于此，在不同的实施环境或实施状态下，采用所述drag控制法对所述交频信号进行叠加处理时可以采用拆分不同时间段以及不同拆分单位的以确定所述交频信号的频率。
所述信号测量单元131用于接收所述待测量子态并通过预设的测量方式对所述待测量子态进行当前多比特门的固有保真度，即当前交频信号下的多比特门的固有保真度。在本实施例中，所述预设的测量方式为x方向，y方向，及z方向三个方向上的投影测量，例如所述测量模块13采用randomizedbenchmarking的方法测定当前控制信号下单次多比特门ut的固有保真度(intrinsicfidelity,i.e.,fidelityneglectingdecoherence)，所述固有保真度由
在给定信号作用时间4σ和频率的情况下，对于确定的参数序列{(an,bn)}，可以获得当前多比特操作的固有保真度
为进一步阐述本申请基于超导量子计算的多比特操控系统的原理及功效，在一个实施例中，藉由本申请所述的系统实现的基于超导量子计算的多比特操控例如为以下步骤：
在物理底层中，将为多个串联结构的超导传输子(transmon)122通过上述buscoupling方案耦合在共面波导谐振器的共面波导谐振腔上121，并将各所述超导传输子(transmon)制备到一种量子态，例如为基态：
请参阅图4，显示为本申请的基于超导量子计算的多比特操控系统在再一实施例中的示意图，如图所示，在另一个实施例中，所述基于超导量子计算的多比特操控系统还包括优化参数生成模块，优化参数生成模块用于生成预设的初始微波序列给所述微波输入模块11，所述生成预设的初始微波序列是经预优化处理的，进而可以控制初始输入的始微波序列样本，而进一步保证操作的准确性和高效性。在图3所示的实施例中，所述优化参数生成模块在初始情况下随机生成m组参数序列{(an，bn)}1，{(an，bn)}2，…，{(an，bn)}m，并将输出给所述微波输入模块11。以令所述微波输入模块11以交频信号
本申请基于超导量子计算的多比特操控系统使用共面波导谐振模块将多个超导传输子(transmon)耦合起来，并通过简单的微波元器件对这些量子比特进行操控，同时反馈数据并根据设计方案进行优化，获得最优控制(optimalcontrol)参数，实现单次多比特门操作。本申请提供的技术方案具有可拓展性高，保真度高，操作时间短的优点，应用场景普遍，且成本较低。
本申请还提供一种基于超导量子计算的多比特操控方法，使用共面波导谐振模块将多个超导传输子(transmon)耦合起来，并通过简单的微波元器件对这些量子比特进行操控，同时反馈数据并根据设计方案进行优化，获得最优控制(optimalcontrol)参数，实现单次多比特门操作，藉此解决如前所述的通过将多比特门拆分为单比特门和双比特门组合的方式将耗费更多的操控时间，并带来更低的保真度。
请参阅图5，显示为本申请的基于超导量子计算的多比特操控方法在一实施例中的流程示意图，如图所示，所述基于超导量子计算的多比特操控方法包括以下步骤：
在步骤s11中，将初始量子态及初始微波序列以交频信号输出；在实施例中，将初始量子态及初始微波序列以交频信号输出的步骤可以通过一微波装置来实现，所述微波装置可以通过商用的微波元器件实现，例如为微波发生器等。
在本实施例中，所述输出的交频信号是经拆分不同时间段的信号叠加处理的交频信号，在实施例中，所述经拆分不同时间段的信号叠加处理的交频信号为不同时间段的drag信号的叠加，使得所有比特门受同一交频信号的调制，避免了对不同比特门进行单独调制时的信道串扰。由于将信号依据不同的时间段分拆为的分段drag信号的组合，也很大程度的减少了泄露误差(leakageerror)。
在步骤s12中，依据接收的所述交频信号操控内置的多比特门并输出待测量子态；在实施例中，所述多比特门被置于一超导量子模块或超导量子装置中，所述超导量子模块或超导量子装置为超导量子计算机的物理底层模型架构，
所述待测量子态是经过对输入的初始初始量子态通过当前交频信号操控内置的多比特门后输出的另一种形式的量子态，在不同的实施方式中，所述待测量子态亦可被称之为测试量子态或制备量子态等等。所述待测量子态对应调制的目标多比特门ut，其中t表示当前时刻。在实施例中，所述超导量子模块或超导量子装置包括例如为共面波导谐振器(coplanarwaveguideresonator)的共面波导(coplanarwaveguide，cpw)谐振腔，所述共面波导谐振腔包括多个串联结构的超导传输子(transmon)，每一所述超导传输子(transmon)由库珀(cooper)对盒构成，所述的多个串联结构的超导传输子(transmon)通过总线与所述共面波导谐振腔相耦合，即buscoupling方案。
所述transmon为在超导量子计算中构成量子比特的物理单位，在本申请中transmon被称之为超导传输子，在相关的技术文献中，transmon还被称之为超导量子比特，或者超导电荷，又或者超导电荷量子位，这些名称均可表示为本申请所述的超导传输子(transmon)。
在本申请提供的实施例中，所述的buscoupling方案是将制备好的超导传输子(transmon)按列排布耦合到高品质因子的共面波导谐振腔中，共面波导谐振腔中的光子可以诱导超导传输子(transmon)相互之间的远程强耦合，从而构建了超导量子计算机的物理底层模型。
在步骤s13中，对所述待测量子态并进行测量以获得当前多比特门的固有保真度；在实施例中，藉由测量装置或测量设备对所述待测量子态并进行测量，接收所述待测量子态并进行测量的目的是为了获得当前多比特门的固有保真度，所述测量模块为由多个测量设备和设备构成的信号测量端。
通过预设的测量方式对所述待测量子态进行当前多比特门的固有保真度，即当前交频信号下的多比特门的固有保真度。在本实施例中，所述预设的测量方式为x方向，y方向，及z方向三个方向上的投影测量。在某些实施例中，所述预设的测量方式可以为rb(randomizedbenchmarking)测量，断层扫描(tomography)测量(请参阅d'ariano，etal.physicalreviewletters86.19(2001):4195.)，门集合断层扫描(gatesettomography)测量(请参阅greenbaum,etal.arxivpreprintarxiv:1509.02921(2015).)，或者误差估计(errormitigation)测量(请参阅temme,kristanetal.physicalreviewletters119.18(2017):180509.)等，在本实施例中，优先选用rb(randomizedbenchmarking)测量，容后详述。
在步骤s14中，依据所述固有保真度进行计算生成优化微波序列，并将所述优化微波序列反馈进行再次操控。在一些实施例中，利用bfgs算法(即broyden-fletcher-goldfarb-shanno，简称bfgs)或globalsearch算法对所述初始微波序列进行优化处理，其中，globalsearch算法通过matlab软件系统实现。在本实施例中，优选选用bfgs算法对所述初始微波序列进行优化处理。
在步骤s15中，判断获得多比特门是否目标多比特门，也就是判断获得的优化微波序列是否最优控制参数序列{(an,bn)}optimal，如果是，则结束所述方法的步骤，若否，则返回步骤s11中，将所述优化微波序列反馈至微波装置，经以上过程将经历多次重复输入、优化、反馈过程若干次，最后给出最优控制参数序列{(an,bn)}optimal，最后达到最优调制效果，进而实现单次多比特门操作。
在具体的实现方式中，所述反馈的操作可以通过包含软件和硬件的计算机设备来实现，可例如装载有本申请所述反馈功能的计算机程序的计算机来实现，所述计算机包括但不限于服务器、台式电脑、或笔记本电脑等。
在实施例中，本申请基于超导量子计算的多比特操控方法通过汉密顿量的演化来实现。所述哈密顿量为系统的能量算符，是一个描述系统总能量的算符，以h表示，则所述多比特操控系统的哈密顿量描述为：h＝h0+hd+hc。
其中，h0可以看作超导传输子(transmon)本身的哈密顿量；
hd可以看作所述交频信号调制下的哈密顿量，
hc可以看做超导传输子(transmon)之间相互耦合的哈密顿量，在本申请所提供的实施例中，所述超导传输子(transmon)之间相互耦合的哈密顿量hc是也固定的，在本申请的一实施例中，采用x-y型耦合，即，对耦合哈密顿量近似到二阶，仅保留x-y之间的耦合，忽略z-z之间的耦合；
则
其中，jij表示为第i个和第j个超导传输子之间的耦合强度。
在本申请中，通过所述交频信号来调节或调制所述哈密顿量hd，即所述交频信号
本申请提供的实施方式中采用drag最优控制法，关于drag最优控制法的描述可以参阅doi:10.1103/physreva.82.040305(optimizeddrivingofsuperconductingartificialatomsforimprovedsingle-qubitgates)。在本申请提供的实施例中，例如在200ns或者150ns内时间范围内实现，通过drag控制法将所述交频信号比如以10ns或20ns为单位拆分成信号的叠加形式，以确定交频信号的频率。
所述drag最优控制法表示为
在具体的实现过程中，所述不同的时间段可以依据具体的实施情况进行设定，比如为在150ns的范围内设定，通过drag控制法将所述交频信号比如以10ns为单位拆分成信号的叠加形式，则所述微波输入装置将随机生成或者接收的外部优化参数设备生成的15组参数序列{(an,bn)}1,{(an,bn)}2,…,{(an,bn)}15，即n＝1，2，……15。但并不局限于此，在不同的实施环境或实施状态下，采用所述drag控制法对所述交频信号进行叠加处理时可以采用拆分不同时间段以及不同拆分单位的以确定所述交频信号的频率。
所述接收所述待测量子态并通过预设的测量方式对所述待测量子态进行当前多比特门的固有保真度，即当前交频信号下的多比特门的固有保真度。在本实施例中，所述预设的测量方式为x方向，y方向，及z方向三个方向上的投影测量，例如所述测量采用randomizedbenchmarking的方法测定当前控制信号下单次多比特门ut的固有保真度(intrinsicfidelity，i.e.，fidelityneglectingdecoherence)，所述固有保真度由
在给定信号作用时间4σ和频率的情况下，对于确定的参数序列{(an，bn)}，可以获得当前多比特操作的固有保真度
在本实施例中，所述反馈在获得m组参数序列{(an,bn)}1,{(an,bn)}2，…，{(an，bn)}m对应的固有保真度
为进一步阐述本申请基于超导量子计算的多比特操控方法的原理及功效，在一个实施例中，藉由本申请所述的系统实现的基于超导量子计算的多比特操控例如为以下步骤：
在物理底层中，将为多个串联结构的超导传输子(transmon)通过上述buscoupling方案耦合在共面波导谐振器的共面波导谐振腔上，并将各所述超导传输子(transmon)制备到一种量子态，例如为基态：
请参阅图6，显示为本申请的基于超导量子计算的多比特操控方法在另一实施例中的流程图，如图所示，在另一个实施例中，所述方法还包括生成预设的初始微波序列的步骤，所述预设的初始微波序列是经预优化处理的。通过预先生成的初始微波序列给所述微波输入端，所述生成预设的初始微波序列是经预优化处理的，进而可以控制初始输入的始微波序列样本，而进一步保证操作的准确性和高效性。在图6所示的实施例中，生成预设的初始微波序列的步骤s10在初始情况下随机生成m组参数序列{(an，bn)}1，{(an,bn)}2,…,{(an,bn)}m，并将输出给所述微波输入端，以令所述微波输入端以交频信号
在本申请的第二方面的某些实施方式中，所述方法还包括生成预设的初始微波序列的步骤，所述预设的初始微波序列是经预优化处理的。
本申请基于超导量子计算的多比特操控方法使用共面波导谐振模块将多个超导传输子(transmon)耦合起来，并通过简单的微波元器件对这些量子比特进行操控，同时反馈数据并根据设计方案进行优化，获得最优控制(optimalcontrol)参数，实现单次多比特门操作。本申请提供的技术方案具有可拓展性高，保真度高，操作时间短的优点，应用场景普遍，且成本较低。
应理解，在本申请的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
上述实施例仅例示性说明本申请的原理及其功效，而非用于限制本申请。任何熟悉此技术的人士皆可在不违背本申请的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本申请所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本申请的权利要求所涵盖。