CN113033812A - 量子操作执行方法、装置及量子操作芯片 - Google Patents

量子操作执行方法、装置及量子操作芯片 Download PDF

Info

Publication number
CN113033812A
CN113033812A CN202110357168.8A CN202110357168A CN113033812A CN 113033812 A CN113033812 A CN 113033812A CN 202110357168 A CN202110357168 A CN 202110357168A CN 113033812 A CN113033812 A CN 113033812A
Authority
CN
China
Prior art keywords
quantum
target
odd
vertex
operations
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.)
Granted
Application number
CN202110357168.8A
Other languages
English (en)
Other versions
CN113033812B (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.)
Tencent Technology Shenzhen Co Ltd
Original Assignee
Tencent Technology Shenzhen Co Ltd
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 Tencent Technology Shenzhen Co Ltd filed Critical Tencent Technology Shenzhen Co Ltd
Priority to CN202110357168.8A priority Critical patent/CN113033812B/zh
Publication of CN113033812A publication Critical patent/CN113033812A/zh
Priority to PCT/CN2021/117519 priority patent/WO2022205781A1/zh
Priority to KR1020227005112A priority patent/KR20220137614A/ko
Priority to EP21839311.4A priority patent/EP4089593A4/en
Priority to JP2022513857A priority patent/JP7381723B2/ja
Application granted granted Critical
Publication of CN113033812B publication Critical patent/CN113033812B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
    • G06N10/80—Quantum programming, e.g. interfaces, languages or software-development kits for creating or handling programs capable of running on quantum computers; Platforms for simulating or accessing quantum computers, e.g. cloud-based quantum computing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
    • G06N10/70—Quantum error correction, detection or prevention, e.g. surface codes or magic state distillation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
    • G06N10/20—Models of quantum computing, e.g. quantum circuits or universal quantum computers
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
    • G06N10/40—Physical realisations or architectures of quantum processors or components for manipulating qubits, e.g. qubit coupling or qubit control
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B82—NANOTECHNOLOGY
    • B82Y—SPECIFIC USES OR APPLICATIONS OF NANOSTRUCTURES; MEASUREMENT OR ANALYSIS OF NANOSTRUCTURES; MANUFACTURE OR TREATMENT OF NANOSTRUCTURES
    • B82Y10/00—Nanotechnology for information processing, storage or transmission, e.g. quantum computing or single electron logic

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Superconductor Devices And Manufacturing Methods Thereof (AREA)

Abstract

本申请公开了一种量子操作执行方法、装置及量子操作芯片，涉及量子技术领域。所述量子操作执行方法包括：获取量子电路对应的n个未调度的量子操作；基于量子计算器件的拓扑结构，以及n个所述未调度的量子操作进行操作调度，获得目标操作组；在所述量子计算器件上并行执行所述目标操作组中包含的各个被调度的量子操作；所述被调度的量子操作的控制波形是基于指定目标生成的；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；上述方案能够适用于任意类型和任意实现原理的量子操作，从而极大的提高了串扰抑制的适用性。

Description

量子操作执行方法、装置及量子操作芯片
技术领域
本申请实施例涉及量子技术领域，特别涉及一种量子操作执行方法、装置及量子操作芯片。
背景技术
在量子计算器件中，量子比特之间的耦合会导致量子比特之间存在串扰，而如何抑制量子比特之间的串扰，是量子计算中急需解决的问题。
在相关技术中，量子计算中的串扰抑制通常有硬件方式和软件方式，其中，软件方式主要是在利用交叉共振效应实现的两比特操作的场景中，引入额外的波形来纠正两比特操作的失真，从而减轻两比特操作下的串扰影响。
然后，上述基于软件的串扰抑制方式只适用于利用交叉共振效应实现的两比特操作的场景，导致串扰抑制的适用性较差。
发明内容
本申请实施例提供了一种量子操作执行方法、装置及量子操作芯片，能够通过软件方式实现对任意类型量子操作的串扰抑制，提高了串扰抑制的适用性。所述技术方案如下：
根据本申请实施例的一个方面，提供了一种量子操作执行方法，所述方法包括：
获取量子电路对应的n个未调度的量子操作；n大于等于1，且n为整数；
基于量子计算器件的拓扑结构，以及n个所述未调度的量子操作进行操作调度，获得目标操作组；所述目标操作组包含各个被调度的量子操作；
在所述量子计算器件上并行执行所述目标操作组中包含的各个被调度的量子操作；
其中，所述被调度的量子操作的控制波形是基于指定目标生成的；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是以是否执行量子操作，对所述量子计算器件中的各个量子比特进行划分得到的。
根据本申请实施例的一个方面，提供了一种量子操作的控制波形生成方法，所述方法包括：
获取被调度的量子操作；
基于指定目标，获取所述被调度的量子操作的控制波形函数的函数参数；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是在量子计算器件中，以包含的量子比特是否执行量子操作进行划分的；
基于所述函数参数生成所述被调度的量子操作的控制波形。
根据本申请实施例的一个方面，提供了一种量子操作执行装置，所述装置包括：
量子操作获取模块，用于获取量子电路对应的n个未调度的量子操作；n大于等于1，且n为整数；
量子操作调度模块，用于基于量子计算器件的拓扑结构，以及n个所述未调度的量子操作进行操作调度，获得目标操作组；所述目标操作组包含各个被调度的量子操作；
量子操作执行模块，用于在所述量子计算器件上并行执行所述目标操作组中包含的各个被调度的量子操作；
其中，所述被调度的量子操作的控制波形是基于指定目标生成的；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是所述量子计算器件中的各个量子比特以是否执行量子操作进行区分得到的。
根据本申请实施例的一个方面，提供了一种量子操作的控制波形生成装置，所述装置包括：
量子操作获取模块，用于获取被调度的量子操作；
参数获取模块，用于基于指定目标，获取所述被调度的量子操作的控制波形函数的函数参数；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是以是否执行量子操作，对所述量子计算器件中的各个量子比特进行划分得到的；
波形生成模块，用于基于所述函数参数生成所述被调度的量子操作的控制波形。
根据本申请实施例的一个方面，提供了一种量子操作芯片，所述量子操作芯片用于实现如上所述的量子操作执行方法，或者，所述量子操作芯片用于实现如上所述的量子操作的控制波形生成方法。
本申请实施例提供的技术方案至少包括如下有益效果：
基于量子电路，以及量子计算器件的拓扑结构，进行量子操作的调度，得到被调度的量子操作，并在量子计算器件上执行被调度的量子操作，其中，执行的量子操作的波形，是在量子计算器件上保证量子操作不失真的情况下，使得区域间串扰最小化的控制波形；上述方案基于量子电路的拓扑结构对量子操作的执行顺序和控制波形进行优化即可以实现串扰抑制，不受量子操作的操作类型以及实现原理的限制，能够适用于任意类型和任意实现原理的量子操作，从而极大的提高了串扰抑制的适用性。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一个实施例提供的方案应用场景的示意图；
图2是本申请一个实施例提供的量子操作执行方法的流程图；
图3是本申请一个实施例提供的量子操作执行方法的流程图；
图4是图3所示实施例涉及的一种ZZ串扰抑制示意图；
图5是图3所示实施例涉及的一种最大区域的量子比特数和不可被抑制的ZZ串扰总个数相互制约的示意图；
图6是图3所示实施例涉及的拓扑图及其对偶图；
图7是图3所示实施例涉及的最优抑制求解示意图；
图8是图3所示实施例涉及的最优抑制问题求解示例图；
图9是图3所示实施例涉及的量子电路示意图；
图10是图3所示实施例涉及的抑制方案示意图；
图11是图3所示实施例涉及的抑制效果差异示意图；
图12是图3所示实施例涉及的一种单比特操作示意图；
图13是图3所示实施例涉及的单比特波形示意图；
图14是图3所示实施例涉及的单比特波形示意图；
图15是图3所示实施例涉及的单比特操作和两比特操作的示意图；
图16是图3所示实施例涉及的一种单比特波形拓扑图；
图17是图3所示实施例涉及的串扰抑制情况示意图；
图18是图3所示实施例涉及的另一种串扰抑制情况示意图；
图19是图3所示实施例涉及的一种两比特波形拓扑图；
图20是图3所示实施例涉及的实现Rxx(π/2)操作时对ZZ串扰的抑制效果的对比示意图；
图21是图3所示实施例涉及的网格拓扑；
图22是图3所示实施例涉及的两种调度方式的对比示意图；
图23是图3所示实施例涉及的执行量子电路的示意图；
图24是图3所示实施例涉及的量子内存的示意图；
图25是图3所示实施例涉及的恒等量子电路示意图；
图26本申请一个实施例提供的量子操作执行装置的框图；
图27本申请一个实施例提供的量子操作的控制波形生成装置的框图；
图28是本申请一个实施例提供的计算机设备的结构框图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。
在对本申请实施例进行介绍说明之前，首先对本申请中涉及的一些名词进行解释说明。
1、量子计算(Quantum Computation)：利用量子态的叠加和纠缠等性质快速完成计算任务的一种计算方式。
2、量子比特(Qubit)：量子信息的承载形式，下文中部分位置简写为比特。
3、量子操作(QuantumOperation)：对量子比特进行操纵，从而对量子比特承载的量子信息进行处理。常见的量子操作有泡利X、Y、Z变换(或写作σx、σy、σz)，哈达玛变换(H)，受控泡利X变换(CNOT)等。只使用单比特操作和两比特操作，可完成任意的量子计算，下文中部分位置简写为操作。
4、量子电路(QuantumCircuit)：量子计算的一种描述模型，由量子比特以及在量子比特上的量子操作组成。
5、量子计算器件(Quantum Computing Device)：执行量子计算的物理装置。
6、量子电路执行并行度(ExecutionParallelism)：在量子计算器件上执行量子电路时，同时执行的量子操作个数。
7、哈密顿量(Hamiltonian)：描述量子系统(例如量子比特)演化的数学工具。有效哈密顿量指的是对原哈密顿量进行坐标变换并忽略无关因素后得到的简化哈密顿量。
8、超导量子计算器件(SuperconductingQuantumComputingDevice)：指量子比特采用超导技术实现的量子计算器件。在超导量子计算器件中，施加控制波形(ControlPulse)影响哈密顿量中特定项，可影响量子比特演化，从而实现量子操作。
9、ZZ串扰(ZZCrosstalk)：由于量子比特高能级之间的耦合而产生的
Figure BDA0003003841120000051
相互作用。在超导量子计算器件中，具有物理耦合的两个比特间存在ZZ串扰，这是一种永久存在(always-on)的相互作用。
10、量子内存(QuantumMemory)：保存量子信息的装置，与量子计算器件的不同之处在于，其设计目的不在于计算，而在于长时间保存信息，但量子内存上也可以进行某些量子操作(或称刷新(Refresh))以延长信息保存时间。
11、量子纠错码(Quantum Error Correction Code)：真实的量子计算器件存在各种错误，量子纠错码用于编码量子比特，提高量子计算的容错率。
12、保真度(Fidelity)：衡量实际值与真值相似度的指标，具体场景中有不同定义形式。保真度越高，实际值与真值越相似。
13、图(Graph)：数学上，图由顶点与连接顶点的边组成，记为G＝(V,E)，其中V是顶点集，E是边集。量子计算器件的拓扑结构可用图表示，其中顶点代表量子比特，边代表量子比特间的耦合。
14、平面图(PlanarGraph)、对偶图(DualGraph)：能够以某种方式画在平面上并保证不同边不相交的图称为平面图。平面图存在对偶图，对偶图的顶点对应平面图的面，当平面图中两个面相邻时，对偶图中对应顶点之间连边，这样对偶图的边就与平面图的边一一对应。平面图G的对偶图记为G*，平面图的边(u,v)对应对偶图的边(u,v)*。
15、最短路(ShortestPath)：连接两个顶点的所有路径中最短的路径。
16、连通分量(ConnectedComponent)：连通分量是图的子图，满足各顶点之间均存在路径相连的条件(连通性条件)。当增加图中的任何其他顶点都会违背连通性条件时，此时的连通分量称为最大连通分量。
17、割(Cut)、割集(Cut-Set)：图G＝(V,E)的割C＝(S,T)是将图的顶点集V划分成两个不相交的集合S和T的一种划分，其对应的割集EC是图的边集E中两端点分属S和T的边组成的集合，即EC＝{(u,v)|u∈S,V∈T}。
18、奇度顶点配对(Odd-vertexPairing)：当移除一个边集D后剩下的图中没有奇数度的顶点，则称该边集为奇度顶点配对。一个顶点的度数为其所连边的个数。任何奇度顶点配对均包含图中所有奇度顶点。
19、最大匹配(Maximum Matching)：图的一个匹配是图的一个子图，其中每两条边都没有公共顶点。图的最大匹配指的是边数最多的匹配。
20、顶点二着色(Vertex 2-Coloring)：为顶点赋予两种颜色中的一种，使得有边相连的顶点颜色不同。如果一个图可以顶点二着色，则称该图是二分图(BipartiteGraph)。
21、搜索(Search)：在问题的众多可能解中，高效地求得较优解的方法。常见搜索方法有穷尽搜索(Exhaustive Search)：穷尽所有可能选取最优解，集束搜索(BeamSearch)：问题的解由若干个子问题的解组成，在求解过程中对于每个子问题仅保留有限个数解的搜索方法。
在相关技术中，抑制ZZ串扰的技术方案可分为四类：
1)采用可调节的耦合。ZZ串扰需要通过量子比特间的耦合介导，在不需要耦合时，可通过调低耦合强度，等效地“关闭”耦合，从而抑制串扰。
2)采用异构量子比特。ZZ串扰的强度与耦合相连的两个量子比特的非谐性有关，当它们的非谐性等大反号(绝对值相同，符号相反)时，ZZ串扰强度为0。通过使用异构的量子比特，可以实现非谐性的等大反号，从而实现对ZZ串扰的抑制。
3)引入多种耦合方式。不同耦合方式介导的ZZ串扰具有不同的强度，特别是它们的符号可以相反。通过引入多种耦合方式，使得不同方式介导的串扰强度相互抵消，可抑制ZZ串扰。
4)施加额外波形。在执行两比特操作时，ZZ串扰会使两比特操作失真，通过施加额外波形，可以一定程度上纠正这种失真，减轻ZZ串扰的影响。
上述技术方案中前三类为硬件实现方案，其缺点在于需要复杂的硬件结构与控制逻辑。例如，可调节耦合方案需要引入额外的耦合器和控制耦合器的控制线，异构量子比特对芯片的制造工艺提出了更高的要求，多种耦合方式则需要精确地控制不同耦合的参数，这些复杂的硬件结构增加了控制复杂性并引入了额外的退相干因素。与这些方案相比，本申请后续实施例提供的方案可用在具有单个固定耦合的同构量子比特芯片上；此外，本申请所示的的方案也可用于这三类硬件实现方案的硬件上，具有良好的通用性。
而对于施加额外波形这一方案，其缺点在于适用性差。这一方案针对利用交叉共振效应实现的两比特操作设计，目前未有工作表明其可用于其他方式实现的操作或单比特操作。与之相比，本申请后续实施例提供的方案既适用于单比特操作也适用于双比特操作，且对操作的实现技术没有特殊要求。
总而言之，本申请提出了一种更加简单、通用的串扰抑制方法，无需复杂硬件结构支持，且具有高适用性。
请参考图1，其示出了本申请一个实施例提供的方案应用场景的示意图。如图1所示，该应用场景可以是超导量子计算平台，该应用场景包括：量子计算器件11、稀释制冷机12、控制设备13和计算机14。
量子计算器件11是一种作用在物理量子比特上的电路，量子计算器件11可以实现成为量子芯片，如处于绝对零度附近的超导量子芯片。稀释制冷机12用于为超导量子芯片提供绝对零度的环境。
控制设备13用于对量子计算器件11进行控制，计算机14用于对控制设备13进行控制。例如，编写好的量子程序经过计算机14中的软件编译成指令发送给控制设备13(如电子/微波控制系统)，控制设备13将上述指令转换为电子/微波控制信号输入到稀释制冷机12，控制处于小于10mK温度的超导量子比特。读取的过程则与之相反，读取波形被输送到量子计算器件11。
在对本申请方法实施例进行介绍说明之前，先对该方法的运行环境进行介绍说明。本申请实施例提供的量子操作执行方法，其可以由经典计算机(如PC(PersonalComputer，个人计算机))执行实现，例如通过经典计算机执行相应的计算机程序以实现该方法；也可以在经典计算机和量子计算机的混合设备环境下执行，例如由经典计算机执行控制波形生成等步骤，而由量子计算机实现基于控制波形的量子操作执行等步骤。
在下述方法实施例中，为了便于说明，仅以各步骤的执行主体为计算机设备进行介绍说明。应当理解的是，该计算机设备可以是经典计算机，也可以包括经典计算机和量子计算机的混合执行环境，本申请实施例对此不作限定。
请参考图2，其示出了本申请一个实施例提供的量子操作执行方法的流程图。该方法各步骤的执行主体可以是计算机设备。该方法可以包括如下几个步骤：
步骤21，获取量子电路对应的n个未调度的量子操作；n大于等于1，且n为整数。
上述量子电路对应的n个未调度的量子操作，是量子电路中未被执行的部分或者全部量子操作。
步骤22，基于量子计算器件的拓扑结构，以及n个该未调度的量子操作进行操作调度，获得目标操作组；该目标操作组包含各个被调度的量子操作。
在本申请实施例中，量子操作的操作调度是指：确定在某个时间点上，在量子计算器件中并行执行的一个或者多个量子操作。
步骤23，在该量子计算器件上并行执行该目标操作组中包含的各个被调度的量子操作；该被调度的量子操作的控制波形是基于指定目标生成的；该指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；该区域间串扰是分属于不同的连通区域的量子比特之间的串扰；该连通区域是以是否执行量子操作，对该量子计算器件中的各个量子比特进行划分得到的。
在本申请实施例中，量子计算器件中的各个量子比特在某一时刻A可以分为两类，一类是在该时刻A执行量子操作的量子比特(可以多个量子操作并行执行)，另一类是时刻A未执行量子操作的量子比特；这两类量子比特在量子电路的拓扑结构中形成至少两个连通区域。
对于一个连通区域，该连通区域中全部量子比特的类型相同，并且，当该连通区域包含两个或者两个以上量子比特时，该连通区域中的任意一个量子比特与该连通区域中的至少一个其它量子比特相邻(两个量子比特相邻，是指两个量子比特之间存在边，或者说存在耦合)。
并且，对于量子计算器件中的相邻的两个量子比特，若在时刻A，这两个量子比特的类型不同，则这两个量子比特在时刻A分属于不同的连通区域。
其中，上述区域间串扰可以是上述连通区域间的ZZ串扰，或者，也可以是上述连通区域间的其他类型串扰。
综上所述，本申请实施例所示的方案，基于量子电路，以及量子计算器件的拓扑结构，进行量子操作的调度，得到被调度的量子操作，并在量子计算器件上执行被调度的量子操作，其中，执行的量子操作的波形，是在量子计算器件上保证量子操作不失真的情况下，使得区域间串扰最小化的控制波形；上述方案基于量子电路的拓扑结构对量子操作的执行顺序和控制波形进行优化即可以实现串扰抑制，不受量子操作的操作类型以及实现原理的限制，能够适用于任意类型和任意实现原理的量子操作，从而极大的提高了串扰抑制的适用性。
在本申请的一个实施例中，还提供一种量子操作的操作波形生成方法，该方法包括：获取被调度的量子操作；基于指定目标，获取该被调度的量子操作的控制波形函数的函数参数；该指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；该区域间串扰是分属于不同的连通区域的量子比特之间的串扰；该连通区域是在量子计算器件中，以包含的量子比特是否执行量子操作进行划分的；基于该函数参数生成该被调度的量子操作的控制波形。
本申请上述图2所示实施例在波形控制层面提出了一种量子操作的控制波形的优化方案，通过该方案得到的控制波形能够在实现量子操作的同时有效地抑制区域间串扰；在此基础上，还进一步在量子电路层面提出一种基于量子电路的操作调度方案(也可以称为电路调度/量子电路调度)。该调度方案能够在实现最优抑制的同时保证量子电路执行的高并行度，从而最大程度地提高量子电路执行的保真度。
本申请所示的方案可以工作在软件层面，基于通用的量子比特控制模型，能够避免现有方案复杂的硬件组件与控制逻辑，从而更为简单、通用。本方案可应用于超导量子计算器件以提供抑制串扰的高保真度量子计算，亦可应用于改进量子内存，提高量子数据的保存时间。
请参考图3，其示出了本申请一个实施例提供的量子操作执行方法的流程图。该方法各步骤的执行主体可以是计算机设备。该方法可以包括如下几个步骤：
步骤301，获取量子电路对应的n个未调度的量子操作；n大于等于1，且n为整数。
在本申请实施例中，当量子计算系统中存在待执行的量子计算任务时，计算机设备可以获取到该量子计算任务对应的量子电路，进而获取到量子电路对应的n个未调度的量子操作。
其中，上述未调度的量子操作，是指尚未被计算机设备调度在量子计算器件上执行的量子操作。
步骤302，基于n个未调度的量子操作获取可调度集合；该可调度集合中包含n个未调度的量子操作中可调度的量子操作。
在一个量子电路中，两个或者两个以上的量子操作之间可能会存在依赖关系，或者说是执行顺序关系；也就是说，一个量子操作A，需要在前一个量子操作B已经执行完成的情况下，才可以被调度执行，上述的量子操作B可以称为量子操作A的前置量子操作。
在一种可能的实现方式中，上述n个未调度的量子操作中可调度的量子操作，可以是量子电路未被调度执行的量子操作中，当前时刻可执行的量子操作；其中，一个量子操作可执行，可以是指量子电路中不存在该量子操作对应的，未执行的前置量子操作。
步骤303，从该可调度集合中提取出候选操作组，并获取量子比特划分方案。
其中，该量子比特划分方案是基于该拓扑结构，按照指定划分方式得到的；该量子比特划分方案指示将该量子计算器件中的各个量子比特划分为两个量子比特集合，以在该量子计算器件上形成至少两个该连通区域。
在本申请实施例中，计算机设备可以基于可调度集合中的量子操作是否包含两比特量子操作，选择使用不同的方案来提取候选操作组，以及使用不同的方案来获取量子比特划分方案。下面分两个方面，对可调度集合中的量子操作是否包含两比特量子操作的量子比特划分方案进行介绍。
一)可调度集合中的各个量子操作为单比特量子操作。
在一种可能的实现方式中，响应于上述可调度集合中的各个量子操作为单比特量子操作，计算机设备将可调度集合获取为上述候选操作组。
当可调度集合中的各个量子操作均为单比特量子操作时，计算机设备不需要考虑是否有同一个量子操作对应的两个比特被同时划分到两个集合的情况，只需要将可调度集合中的全部量子操作构成的操作组，获取为候选操作组。
其中，若候选操作组是可调度集合中包含的量子操作时单比特量子操作时提取的操作组，则计算机设备可以基于该拓扑结构，按照第二指定划分方式获得目标奇度顶点配对方案；然后，在拓扑图中去除目标奇度顶点配对方案对应的边，获得去边后的拓扑图；再对去边后的拓扑图进行顶点二着色，获得上述量子比特划分方案。
在一种可能的实现方式中，上述基于该拓扑结构，按照第二指定划分方式获得目标奇度顶点配对方案，包括：
获取完全图，该完全图中的顶点为目标对偶图中的奇度顶点；该目标对偶图是拓扑图的对偶图，该拓扑图是该量子计算器件的拓扑结构的平面图；
在该完全图中确定配对的各个奇度顶点对；
获取各个该奇度顶点对的至少一种奇度顶点配对方案；该奇度顶点配对方案中包含各个该奇度顶点对分别在的该目标对偶图中的一条路径；
获取至少一种该奇度顶点配对方案各自的最大区域量子比特数，以及，至少一种该奇度顶点配对方案各自的不可抑制的串扰总数；该不可抑制的串扰总数由该奇度顶点配对方案中的路径对应在该拓扑图中的边数进行表征；该最大区域量子比特数由该拓扑图除去割集之后的最大连通区域中的顶点数进行表征，该割集是该拓扑图中除去该奇度顶点配对方案中的路径对应在该拓扑图中的边之后剩余的边组成的集合；
基于至少一种该奇度顶点配对方案各自的最大区域量子比特数，以及，至少一种该奇度顶点配对方案各自的不可抑制的串扰总数，获取至少一种该奇度顶点配对方案的串扰抑制指数；
将对应该串扰抑制指数最小的该奇度顶点配对方案，获取为目标奇度顶点配对方案。
在一种可能的实现方式中，上述获取各个该奇度顶点对的至少一种奇度顶点配对方案，包括：
获取各个该奇度顶点对各自在该完全图中的前k短路径；k≥1，且k为整数；
对各个该奇度顶点对各自在该完全图中的前k短路径进行组合，获得至少一种该奇度顶点配对方案。
在一种可能的实现方式中，上述串扰抑制指数为：αNQ+NC；
其中，NQ表示该最大区域量子比特数，NC表示该不可抑制的串扰总数，α表示NQ和NC的重要程度。
结合上述内容的介绍，当可调度集合中的各个量子操作为单比特量子操作时，计算机设备不需要考虑是否有同一个量子操作对应的两个比特被同时划分到两个集合的情况，只需要将可调度集合中的全部量子操作构成的操作组，获取为候选操作组，而后按照第二指定划分方式，得到最终的目标奇度顶点配对方案，再基于目标奇度顶点配对方案得到量子比特划分方案。
其中，对于第二指定划分方式，我们旨在通过选择性地在一些量子比特上施加恒等操作来最大程度抑制ZZ串扰。请参考图4，其示出了本申请实施例涉及的一种ZZ串扰抑制示意图。其中，图4展示了一个在9量子计算器件上抑制ZZ串扰的例子。在这个例子中，通过在黑色顶点上施加恒等量子操作，量子计算器件上所有的ZZ串扰均被抑制，本申请中，可以称这种情况为完全抑制。
其中，通过在量子比特上施加量子操作，将每个量子比特划分在单独的区域内，使得量子计算器件上所有的ZZ串扰均被抑制，则称该量子计算器件存在完全抑制方案。
其中，若一个量子计算器件的拓扑图是二分图，则该量子计算器件存在完全抑制方案。
对于不存在完全抑制的拓扑结构，至少存在一个区域，其包含两个或者两个以上的量子比特，这两个或者两个以上的量子比特组成的连通区域内的串扰不被抑制。对于这种情况，需要关注两个指标：最大区域的量子比特数和不可被抑制的ZZ串扰总个数；这两个指标相互制约。请参考图5，其示出了本申请实施例涉及的一种最大区域的量子比特数和不可被抑制的ZZ串扰总个数相互制约的示意图。
为了表征这种制约，本申请提出了“最优抑制”概念：
通过在量子比特上施加量子操作，将量子计算器件中的量子比特划分在若干连通区域内，记最大连通区域的量子比特数为NQ，不可被抑制的ZZ串扰总个数为NC，给定表征二者重要程度的系数α，若某种抑制方案使得αNQ+NC最小，则称该量子计算器件存在最优抑制方案。
例如，在图5所示的拓扑结构中，图5中的(a)部分所示的方案为最优抑制方案时，对应α＝0，图5中的(b)部分所示的方案为最优抑制方案时，对应α＝1.5，图5中的(c)部分所示的方案为最优抑制方案时，对应α＝2。
为了求解最优抑制方案，可以将其转化为图的割问题：
图G＝(V,E)的割C＝(S,T)，是将图的顶点集V划分成两个不相交的集合S和T的一种划分方式，其对应的割集EC是图的边集E中两端点分属S和T的边组成的集合，即EC＝{(u,v)|u∈S,V∈T}。
寻找最优抑制方案的等价问题：给定图G＝(V,E)和重要度系数α，寻找一个割C＝(S,T)，使得αNQ+NC最小，其中NQ为G移除割集EC后的最大连通分量的顶点数，NC为G移除割集EC后剩下的边的总数。
通过割导出的两个顶点集S和T分别对应“施加波形的顶点集”和“不施加波形的顶点集”，割集对应的耦合跨越不同区域，因此割集对应的耦合介导的串扰可被抑制，与之相反，对于割集的补对应的耦合，由于它们相连的两个量子比特属于同一连通区域内，因此割集的补对应的耦合介导的串扰不可被抑制。
在给出上述图问题的解后，只需要对“施加波形的顶点集”中的量子比特施加量子操作，即可实现最优抑制。例如图4中，黑白两色顶点分别构成了顶点集S和T，虚线边含于割集内，故这些耦合介导的串扰可被抑制，而实线边不属于割集，其介导的串扰不可被抑制。
解决上述图的割问题是困难的：考虑α＝0的情况，需要找到一个割，使得除去割集后剩下的边最少，即割集包含的边最多，这是最大割问题，属于非确定性多项式时间完全(NP-Complete)问题，不存在普适的多项式算法。因此，本申请考虑采用启发式方法进行求解。此外，考虑到大多数拓扑量子纠错码的拓扑结构均为平面图结构，本申请针对平面图设计本申请的方案。本申请所示的方案基于以下概念：
奇度顶点配对：当移除一个边集D后剩下的图中没有奇数度的顶点，则称该边集为奇度顶点配对。任何奇度顶点配对均包含图中所有奇度顶点，因为任何图中奇度顶点的个数均为偶数，因而称边集D为配对。
对偶图：平面图存在对偶图，对偶图的顶点对应平面图的面，当平面图中两个面相邻时，对偶图中对应顶点之间连边，这样对偶图的边就与平面图的边一一对应。平面图G的对偶图记为G*，平面图的边(u,v)对应对偶图的边(u,v)*。
定理：仅当图G的一个边集D的补
Figure BDA0003003841120000146
在对偶图G*中对应的边集
Figure BDA0003003841120000147
是对偶图G*的一个奇度顶点配对时，图G的一个边集D含于一个割集EC。
例如，请参考图6，其示出了一种拓扑图及其对偶图。其中，图6给出了图5所示例子的拓扑图及其对偶图，图中边界上只有一个顶点的边均与s相连对偶图中的奇度顶点为i和j。图5中的(a)部分中割集的补为{(10，11)，(4，11)，(11，12)}，对应对偶图中{(c，i)，(c，d)，(d，j)}是对偶图的一个奇度顶点配对；类似地，图5中的(b)部分割集的补以及对应对偶图的奇度顶点配对为：{(10，11)，(11，12)，(3，4)，(4，5)→(c，i)，(d，j)，(c，s)，(d，s)}。
引理：图G的一个边集D含于一个割集EC，则移除割集EC剩下的边
Figure BDA0003003841120000141
根据上述定理及引理，在拓扑图中寻找割的问题，可转化为在对偶图中寻找奇度顶点配对
Figure BDA0003003841120000142
的问题。考虑到奇度顶点配对需要包含所有奇度顶点，则最简单的奇度顶点配对仅仅包含若干简单路径，每条路径连接两个奇度顶点，且不同路径连接的奇度顶点各异(如图6中{(c，i)，(c，d)，(d，j)}为一条连接奇度顶点i和j的简单路径)。此外，由于拓扑图中移除割集剩下的边
Figure BDA0003003841120000143
因而可以为上述每条路径选取最短路，这样试图最小化
Figure BDA0003003841120000144
也就是最小化
Figure BDA0003003841120000145
这可以作为求解最优抑制问题的起点。为了综合考虑最大连通分量的顶点数NQ，本方案可以对最短路进行松弛：通过考虑前k短路，寻找最优抑制问题的解。请参考图7，其示出了本申请实施例涉及的最优抑制求解示意图，本申请所示的最优抑制求解算法分为三个步骤。
步骤一：顶点配对。
构造以对偶图中奇度顶点为顶点的完全图，边(u,v)的边权为M-d(u,v)，d(u,v)是u,v间最短路长度，M＝1+max{d(u,v)|u,v是对偶图中的奇度顶点}。在完全图上采用最大匹配算法确定配对方案，这样，在相互匹配的顶点间采用最短路得到的奇度顶点匹配所含路径最少，可作为求解问题的起点。例如，针对图11，构造以顶点i和j为顶点的完全图，(i,j)边权为最短路径长度(1+3)-3＝1。由于只有两个点，因而最大匹配即为{(i,j)}。
步骤二：路径松弛。
对于每个匹配的顶点对，依次尝试其前k短路，对于每次尝试，求得NQ和NC，计算αNQ+NC，求得问题的解。例如，i、j间最短路为(1)i→c→d→j，对应的NQ＝4,NC＝3，考虑其前3短路，有(2)i→c→s→d→j，(3)i→c→s→l→k→j，根据α的不同，会得到不同的解。
步骤三：导出割。
根据步骤二得到的奇度顶点配对，在拓扑图中去除奇度顶点配对对应的对偶边，剩下的图是二分图，故通过顶点二着色可导出割C＝(S,T)。例如，步骤二得到的三种方案(1)、(2)、(3)依次对应图11中的(a)、(b)、(c)。
二)可调度集合中的各个量子操作包含两比特量子操作。
在一种可能的实现方式中，响应于该可调度集合中包含两比特量子操作，获取该可调度集合中的各个两比特量子操作，得到两比特操作集合；响应于该两比特操作集合中包含3个或者3个以上操作，基于该两比特操作集合中，距离最近的两个量子操作构建两个操作组；基于该两比特操作集合中，除了该两个操作之外的其它量子操作与该两个操作组之间的距离，执行i次操作添加，i大于等于1，且i为整数；基于该i次操作添加的结果，获取该候选操作组。
在本申请实施例中，当可调度集合中包含两比特量子操作时，计算机设备可以先对其中的两比特量子操作进行调度，也就是以两比特操作为主设计量子比特划分方案。
为了尽可能的提高量子操作的并行度，计算机设备可以选取同时执行不会降低抑制效果的操作。在量子计算器件中，距离较近的两个量子操作同时执行时，通常会降低串扰抑制效果，而距离较远的两个量子操作同时执行，则不会对串扰抑制效果产生较大的影响。因此，在本申请实施例中，当两比特操作集合中包含3个或者3个以上操作时，计算机设备可以首先选取距离最近的两个两比特量子操作，将其分在不同的组，以确保两者不会同时执行，而后基于组外的操作到对应的操作组之间的距离，逐步扩充这两个组的规模，以在保证串扰抑制效果的前提下，获取到具有尽可能大的并行度的候选操作组。
在一种可能的实现方式中，基于该两比特操作集合中，除了该两个操作之外的其它量子操作与该两个操作组之间的距离，执行i次操作添加，包括：
在第j次操作添加过程中，获取该两比特操作集合中的该其它量子操作分别与该两个操作组之间的距离；1≤j≤i，且j为整数；将最大距离对应的指定量子操作，添加入该最大距离对应的指定操作组；该最大距离是该其它量子操作分别与该两个操作组之间的距离中的最大值；基于该拓扑结构以及该指定操作组，按照第一指定划分方式获得目标奇度顶点配对方案；该目标奇度顶点配对方案包括目标对偶图中的各个奇度顶点对之间的路径，且该目标奇度顶点配对方案用于指示对应的操作组的该量子比特划分方案；该目标对偶图是拓扑图的对偶图，该拓扑图是该拓扑结构的平面图；响应于该目标奇度顶点配对方案满足串扰抑制要求，将该指定量子操作从该两比特操作集合中删除；响应于删除该指定量子操作之后的该其它量子操作的数量不为0，执行下一次操作添加过程。
在一种可能的实现方式中，该基于该至少一次操作添加的结果，获取该候选操作组，包括：
响应于该第i次操作添加后的该其它量子操作的数量为0，或者，响应于该第i次操作添加过程中得到的该目标奇度顶点配对方案不满足该串扰抑制要求，将该第i次操作添加后的该两个操作组中，包含操作数量最多的操作组获取为该候选操作组。
在一种可能的实现方式中，该在第j次操作添加过程中，获取该两比特操作集合中的该其它量子操作分别与该两个操作组之间的距离，包括：
在第j次操作添加过程中，获取第一量子操作与第一操作组中已有的各个量子操作之间的操作距离；该操作距离用于指示两个两比特操作对应的量子比特在该拓扑结构中的最短路径长度之和；该第一量子操作是该其它量子操作中的任意一个，该第一操作组是该两个操作组中的任意一个；
基于该第一量子操作与该第一操作组中已有的各个量子操作之间的操作距离，获取该第一量子操作与该第一操作组之间的距离。
在本申请实施例中，对于两个两比特量子操作a＝(a1,a2)和b＝(b1,b2)，两者间的距离定义为D(a,b)＝∑i∑jd(ai,bj)，其中，d(ai,bj)表示量子比特ai和bj之间的最短路长度。操作a＝(a1,a2)和一组操作G的距离定义为：
ming∈GD(a,g)。
由于量子计算器件的拓扑结构可以通过拓扑图来表示，而操作组中的操作对应的量子比特，对应的就是图中的顶点，量子比特之间的耦合关系，对应的就是图中的边；因此，在本申请实施例中，计算机设备可以通过在拓扑图的对偶图中寻找奇度顶点配对方案的方式，来搜索满足串扰抑制要求的量子比特划分方案。而在逐步扩充两个操作组，以获取到具有尽可能大的并行度的候选操作组的过程中，每次向一个操作组中添加一个两比特量子操作之后，即基于添加操作后的操作组中的两比特量子操作，结合量子计算器件的拓扑结构获取对应的目标奇度顶点配对方案，并判断获取到的目标奇度顶点配对方案是否满足串扰抑制要求，如果满足，则继续添加，如果不满足，则需要停止添加；当不满足串扰抑制要求时触发停止添加，或者，不存在需要添加的两比特量子操作时，计算机设备可以将当前两个操作组中包含量子操作最多的一个作为候选操作组；可选的，如果当前两个操作组中包含的量子操作数量相同，则可以随机选择一个作为候选操作组。
在一种可能的实现方式中，该基于该拓扑结构以及该指定操作组，按照第一指定划分方式获得目标奇度顶点配对方案，包括：
在该目标对偶图中删除该指定操作组中的各个量子操作对应的边；获取以该目标对偶图中的奇度顶点为顶点的完全图；在该完全图中确定各个该奇度顶点对；获取各个该奇度顶点对的至少一种奇度顶点配对方案；该奇度顶点配对方案中包含各个该奇度顶点对分别在的该目标对偶图中的一条路径；在至少一种该奇度顶点配对方案中添加该指定操作组中的各个量子操作对应在该目标对偶图中的边；从至少一种该奇度顶点配对方案中，筛除该两比特操作对应的量子比特不属于同一个连通区域的奇度顶点配对方案，获得筛除后的各个该奇度顶点配对方案；获取筛除后的各个该奇度顶点配对方案各自的最大区域量子比特数，以及，筛除后的各个该奇度顶点配对方案各自的不可抑制的串扰总数；该不可抑制的串扰总数由该奇度顶点配对方案中的路径对应在该拓扑图中的边数进行表征；该最大区域量子比特数由该拓扑图除去割集之后的最大连通区域中的顶点数进行表征，该割集是该拓扑图中除去该奇度顶点配对方案中的路径对应在该拓扑图中的边之后剩余的边组成的集合；基于筛除后的各个该奇度顶点配对方案各自的最大区域量子比特数，以及，筛除后的各个该奇度顶点配对方案各自的不可抑制的串扰总数，获取筛除后的各个该奇度顶点配对方案的串扰抑制指数；将对应该串扰抑制指数最小的该奇度顶点配对方案，获取为目标奇度顶点配对方案。
其中，在可调度集合中包含两比特量子操作的情况下，上述基于目标奇度顶点配对方案，获取量子比特划分方案的过程包括：在该拓扑图中去除目标奇度顶点配对方案对应的边，获得去边后的该拓扑图；对去边后的该拓扑图进行顶点二着色，获得量子比特划分方案(类似上述可调度集合中的各个量子操作为单比特量子操作的情况下，最优抑制算法的步骤三中的导出割的步骤)。
对于可调度集合中包含两比特量子操作的情况，由于两比特量子操作对应的量子比特为两个，因此，如果直接按照上述第二指定划分方式进行量子比特的划分，则可能会将一个两比特量子操作对应的两个量子比特划分到不同的集合中。进一步地，若同时执行多个量子操作，原则上这些操作涉及的量子比特需要属于同一个集合，但第二指定划分方式亦不能保证这一点。对此，本申请实施例在上述第二指定划分方式的基础上进行改进，以适用于在指定量子比特上进行指定量子操作时的量子比特划分场景。该改进后的划分方式即为上述第一指定划分方式。
上述第一指定划分方式对应的情况可视为带有约束的最优抑制问题：要求指定的量子操作所涉及的量子比特在割中属于同一集合。本方案可以采用“先删边，后加边，额外检查”的策略将带有约束的问题转为无约束问题，采用上述第二指定划分方式给出的方法求解。记操作涉及的比特集合为Q，要求Q中的量子比特在割中属于同一集合。
先删边：在对偶图中删去边{(u,v)*|(u,v)∈E且u,v∈Q}。
后加边：在步骤二给出删边后的图的奇度顶点配对后，将方才删去的边加入到奇度顶点配对中。
额外检查：上述两个步骤在某些情况下无法确保Q中的量子比特均在同一集合内，因而在步骤二中，每给出一个割候选时，均需要额外检查该候选是否满足同一集合的条件，如果不满足则放弃该候选。
例如，请参考图8，其示出了本申请实施例涉及的最优抑制问题求解示例图。假设需要在图8中的量子比特11、28上执行两比特操作，先在对偶图中删去(11,28)*＝(i,j)，使用第二指定划分方式给出的最优抑制算法得到奇度顶点配对i→c→d→j，而后加入(i,j)，根据加入(i,j)后的奇度顶点配对导出割，可以检查这个方案满足11、28同在一个集合的条件，因而可作为一个解。
给定拓扑图G＝(V,E)及顶点集
Figure BDA0003003841120000181
边集{(u,v)|(u,v)∈E且u,v∈Q}构成了拓扑图G的若干连通分量C1,C2,…,Cn，上述“先删边，再加边”策略能够确保同一连通分量内的顶点在导出的割中处于同一集合。
根据上述结论，单个两比特量子操作涉及的边集只构成一个连通分量(例如图8中的例子)，本方案的策略能给出合法解。对于多个两比特量子操作，其涉及的边集可能构成多个连通分量，此情形需要“额外检查”来确保不同连通分量的量子比特同属一个集合。
本方案将完整的算法总结在表1中。其中遍历由前k短路组成的奇度顶点配对，可使用包括穷尽搜索(Exhaustive Search)、集束搜索(BeamSearch)在内的搜索算法。
表1
Figure BDA0003003841120000191
在一种可能的实现方式中，该串扰抑制要求包括：
该目标奇度顶点配对方案的最大区域量子比特数小于比特数阈值；
以及，
该目标奇度顶点配对方案的不可抑制的串扰总数小于串扰总数阈值。
在本申请实施例中，在调度这些两比特量子操作时，需要在满足串扰抑制要求的情况下尽可能地提高执行并行度。在本申请中，串扰抑制要求可以是对最大区域的量子比特数NQ和不可被抑制的ZZ串扰总个数NC的综合限制，例如要求二者各自小于某些阈值。
其中，上述算法1对应的是确定一个包含两比特量子操作的操作组，且得到的目标奇度顶点配对方案满足串扰抑制要求的情况下，输出该操作组对应的量子比特划分方案的情况。而在本申请实施例所示的方案中，对于计算机设备首先需要确定出满足串扰抑制要求的候选操作组，也就是说，在上述逐步扩充两个操作组的过程中，每次向一个操作组中添加一个两比特量子操作之后，即基于添加操作后的操作组中的两比特量子操作，运行上述算法1中的顶点配对和路径松弛步骤，得到对应的目标奇度顶点配对方案，然后通过该目标奇度顶点配对方案，确定当前的操作组是否满足串扰抑制要求。
在一个示例性的方案中，当两比特操作集合中只包含1个量子操作时，计算机设备可以将该量子操作获取为候选操作组。
在另一个示例性的方案中，当两比特操作集合中包含2个量子操作时，计算机设备可以先将这2个量子操作加入同一个操作组，并确定收满足抑制要求，若是，则将这2个量子操作组成的操作组作为候选操作组，否则，将这2个量子操作中任意一个操作组成的操作组，作为候选操作组。
在一种可能的实现方式中，若候选操作组是可调度集合中包含两比特量子操作时提取的操作组，则计算机设备可以基于上述拓扑结构，按照第一指定划分方式获得该目标奇度顶点配对方案，然后再基于目标奇度顶点配对方案，获取量子比特划分方案(即上述基于候选操作组重新执行上述算法1)。
在另一种可能的实现方式中，若候选操作组是可调度集合中包含两比特量子操作时提取的操作组，则计算机设备也可以直接获取上述候选操作组提取过程中，按照上述第一指定划分方式获取的目标奇度顶点配对方案，然后再基于目标奇度顶点配对方案，获取量子比特划分方案。
本申请实施例所示的方案，对于执行量子电路来说，需要选择何时执行哪些操作，即对量子电路进行调度，不同的调度方案会导致串扰抑制效果和执行并行度的不同。本申请涉及的调度方法综合考虑这两种因素，在实现最优抑制的同时尽可能地提高电路执行的并行度。本申请涉及的调度算法流程可以如表2所示。
表2
Figure BDA0003003841120000211
Figure BDA0003003841120000221
请参考图9，其示出了本申请实施例涉及的量子电路示意图。其中，图9中的(a)部分示出了量子计算器件的拓扑结构，图9中的(b)部分示出了待执行的量子电路。本申请以执行图9所示的量子电路为例解释上述表2所示的算法。给定量子电路，构造可调度操作集合(对应表2算法中的行2)，当且仅当根据数据依赖，其前序依赖的操作均已被执行时，可以认为一个操作是可以被调度执行的。例如，图9中，未执行任何操作时，可调度操作集合为{H1,H2,…,H8}，当所有H操作执行后，可调度操作集合为{CNOT1,4,CNOT2,5CNOT3,6,X8}。然后分类讨论可调度操作集合中操作的类型(行3、6)。如果只有单比特操作，依据以下结论设计执行方案：
对于一组可调度的单比特操作，可以以最优抑制方式在至多两个时间步内执行这组操作。也就是，采用第二指定划分方式对应的最优抑制算法，可将量子比特分为S,T两个集合，相应地这些单比特操作也被分在这两个集合中，则可以在第一个时间步内在S上施加波形，第二个时间步内在T上施加波形，从而在至多两个时间步内执行完毕。
对于本申请所示的电路调度算法，执行含有操作最多的集合(对应表2算法中的行4、5)，而对于另一个集合中的操作，可以留待后续调度。通过这种方式，在最好的抑制效果下实现了最高的并行度。
例如：请参考图10，其示出了本申请实施例涉及的抑制方案示意图。对图15中的拓扑而言，本申请涉及的最优抑制算法给出图10的量子比特划分方案，黑色顶点含有的H操作个数为4，白色顶点含有的操作只有3个，因而选取黑色顶点集合执行。
另一种情况，如果可执行操作集合含有两比特操作，那么通过本申请所示的方案中，对应表2算法中的行7，以两比特量子操作为主涉及量子比特划分方案。
例如，对于图9中的所有CNOT操作，D(CNOT1,4，CNOT2,5)＝D(CNOT2,5,CNOT3,6)＝6，D(CNOT1,4,CNOT3,6)＝10。针对图9中的CNOT操作，请参考图11，其示出了本申请实施例涉及的不同执行方案的距离和抑制效果差异示意图。
对于两比特量子操作，本申请首先选取距离最近的两个操作，将其分在不同的组，以确保两者不会同时执行；而后逐步扩充这两个组的规模：通过依次选取距离最远的操作，直到不满足抑制要求为止。例如，对于图9中的CNOT操作，选取距离最近的CNOT1,4,CNOT2,5并分为两组，A＝{CNOT1,4}，B＝{CNOT2,5}，对于剩下的CNOT3,6，其距A组比B组远，因而将其加入A组，最终得到A＝{CNOT1,4,CNOT3,6}，B＝{CNOT2,5}。
对于一组可执行的两比特量子操作，多次调用算法3，直到该组操作均被调度，如果最终这组操作被分在了K个时间步执行，则这组操作中距离前K小的操作一定分属不同时间步。也就是说，在距离衡量抑制效果的情形下，执行时间的增加(并行度的减少，退相干错误的增加)能带来抑制效果的提升，因为每增加一个时间步，就会有两个距离最近的操作被分开执行。
在给出上述A,B组后，本方案可以选择操作最多的组执行。例如，对于图9的CNOT操作，选取组A＝{CNOT1,4,CNOT3,6}，此时的量子比特划分方案为图11中的(b)部分所示，根据可执行操作集合{CNOT1,4,CNOT2,5,CNOT3,6,X8}，最终执行CNOT1,4,CNOT3,6,X8。
步骤304，基于候选操作组以及量子比特划分方案，获取目标操作组。
在一种可能的实现方式中，基于候选操作组以及该量子比特划分方案，获取该目标操作组，包括：
基于该量子比特划分方案，将该量子计算器件中的各个量子比特划分为两个量子比特集合；将该候选操作组划分为该两个量子比特集合分别对应的子操作组；该子操作组中的量子操作对应的量子比特，处于对应的该量子比特集合内；将目标子操作组中的量子操作以及填充操作，合并为该目标操作组；该目标子操作组是该两个量子比特集合分别对应的该子操作组中，包含量子操作数量最多的子操作组；该填充操作是填充量子比特对应的恒等量子操作，该填充量子比特是该目标子操作组对应的量子比特集合中，除了该目标子操作组对应的量子比特之外的其它量子比特。
在本申请实施例中，对于获取到的量子比特划分方案对应的两个量子比特集合，计算机设备可以将两个量子比特集合中，对应在量子电路中的操作较多的一个量子比特集合确定为执行操作的目标量子比特集合，对于该目标量子比特集合中除了量子电路中的操作对应的量子比特之外的其它量子比特，可以在这些其它量子比特上填充恒等量子操作，以实现最优的串扰抑制效果。
在一种可能的实现方式中，基于候选操作组以及量子比特划分方案，获取该目标操作组之后，还包括：
响应于在该目标操作组之外，还存在m个未调度的量子操作，基于该拓扑结构，以及m个未调度的量子操作进行操作调度，获得新的该目标操作组；m大于等于1，且m为整数；
在该量子计算器件上并行执行新的该目标操作组中包含的各个被调度的量子操作。
在本申请实施例中，计算机设备每获取到一个时间步上的目标操作组之后，可以判断量子电路中是否还存在剩余的未被调度的量子操作，若是，则基于剩余的未被调度的量子操作，重新执行上述步骤301至步骤304，得到下一个时间步上的目标操作组，计算机设备迭代执行上述步骤301至步骤304，直至量子电路中的量子操作全部被调度完成。
步骤305，在该量子计算器件上并行执行该目标操作组中包含的各个被调度的量子操作。
其中，该被调度的量子操作的控制波形是基于指定目标生成的；该指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；该区域间串扰是分属于不同的连通区域的量子比特之间的串扰；该连通区域是以是否执行量子操作，对该量子计算器件中的各个量子比特进行划分得到的。
在本申请实施例中，计算机设备对于要执行的量子操作，在控制对应的量子操作不失真的情况下，优化量子操作的控制波形，使得区域间串扰被最小化。
在一种可能的实现方式中，该指定目标是基于对应的量子操作在控制波形作用下的平均操作保真度，以及该区域间串扰获得的单标量目标。
在一种可能的实现方式中，该单标量目标的公式为：
Figure BDA0003003841120000251
其中，L为该单标量目标，U(T)表示对应的量子操作在控制波形作用下的的酉变换，UTarget(T)表示对应的量子操作的目标酉变换，F(U(T),UTarget(T))表示U(T)与UTarget(T)之间的平均操作保真度；
Figure BDA0003003841120000252
表示在该区域间串扰影响下的封闭量子系统的演化。
在一种可能的实现方式中，计算机设备基于该单标量目标，对目标量子操作的控制波形函数的函数参数进行R次迭代更新；R≥2且R为整数；该目标量子操作是各个该可调度量子操作以及该恒等量子操作中的任意一个；
基于迭代更新后的该函数参数生成该目标量子操作的控制波形。
在本申请实施例中，计算机设备可以基于上述单标量目标，对待执行的量子操作的控制波形函数的参数进行多轮迭代，每一轮迭代过程中，在上一轮迭代的基础上，对控制波形函数的参数进行更新，直至达到迭代次数，或者参数的变化小于某一阈值。
在一种可能的实现方式中，该基于该单标量目标，对目标量子操作的控制波形函数的函数参数进行R次迭代更新，包括：
获取第r次迭代过程中的第一系统演化，该第一系统演化是在第r个控制波形影响下的封闭量子系统的演化；2≤r≤R，且r为整数；该第r个控制波形是第r-1次迭代更新后的该函数参数对应的控制波形；
基于该第一系统演化获取第二系统演化，该第二系统演化是第r次迭代过程中的该区域间串扰影响下的封闭量子系统的演化；
基于该第二系统演化，以及该目标量子操作在该第r个控制波形作用下的平均操作保真度，获取第r次迭代过程中的该单标量目标；
基于第r次迭代过程中的该单标量目标对第r-1次迭代更新后的该函数参数进行更新，获得第r次迭代更新后的该函数参数。
在本申请实施例中，上述区域间串扰可以是区域间的ZZ串扰，或者，也可以是其他类型串扰。
以上述区域间串扰为ZZ串扰为例，本申请实施例中的波形优化方案旨在给定量子操作，设计能够实现量子操作同时抑制ZZ串扰的波形。本方案首先考虑只有一个单比特操作的区域，然后考虑只有一个两比特操作的区域，最后考虑由多个这两种操作组成的区域。
1)只有一个单比特操作的区域。
请参考图12，其示出了本申请实施例涉及的一种单比特操作示意图。如图12所示，本方案需要在量子比特q上实现单比特操作的同时，抑制其与周围量子比特间的ZZ串扰。
该区域的有效哈密顿量可以写作：
Figure BDA0003003841120000261
其中，Ωx,y(t)表示控制波形，
Figure BDA0003003841120000262
表示量子比特q上的泡利X/Y/Z操作，恒等操作省略不写，
Figure BDA0003003841120000263
表示量子比特q与i之间的ZZ串扰，λi表示其强度。
取λ＝∑iλi，上面的哈密顿量可整理为：
H(t)＝HCtrl(t)+λHXtalk；
其中，
Figure BDA0003003841120000264
Figure BDA0003003841120000265
本方案可以为Ωx(t)与Ωy(t)选取适当的波形，以抑制ZZ串扰HXtalk的影响。
2)只有一个两比特操作的区域。
请参考图13，其示出了本申请实施例涉及的两比特操作示意图。如图13所示，本方案需要在量子比特p,q上实现两比特操作的同时抑制其与周围量子比特间的ZZ串扰。
该区域的有效哈密顿量可以写作：
Figure BDA0003003841120000266
其中，λ＝∑iλi+∑jλj，λ′＝λpq。HCoupling为量子比特的耦合形式，例如，
Figure BDA0003003841120000271
可实现CNOT操作，
Figure BDA0003003841120000272
可实现iSWAP操作。通过为Ω(t)选取适当的波形，以抑制跨区域的ZZ串扰HXtalk的影响。然而这种方法要求硬件能够对两比特耦合HCoupling具有较强的调节能力，请参考图14，其示出了本申请实施例涉及的单比特波形示意图，如果硬件无法提供这种能力，可采用如图14所示的单比特波形进行辅助。
如图14所示，波形A、C是施加在量子比特p上的单比特控制波形；波形B、D是施加在量子比特q上的单比特控制波形；波形E为施加在耦合项HCoupling上的波形。
描述上述控制波形的哈密顿量为：
Figure BDA0003003841120000273
在优化时，只优化波形A～D，而对E直接采用原本用于实现两比特操作的波形。通过这种方式，保证波形E是硬件可实现的，同时引入波形A～D来抑制ZZ串扰。
3)波形优化。
由上面的介绍可以看出，无论是单比特操作还是两比特操作情形，哈密顿量部分可写作：
H(t)＝HCtrl(t)+λHXtalk；
其中，HCtrl(t)是可控制部分，用于实现单/两比特操作，HXtalk是要抑制的跨区域串扰，λ是串扰强度。首先给出以下结论：
封闭量子系统按照H(t)＝HCtrl(t)+λHXtalk进行演化，演化对应的酉变换可写作：
U(t)＝UCtrl(t)UXtalk(t)；
其中，UCtrl(t)由
Figure BDA0003003841120000274
定义，表示控制波形下系统的演化。
UXtalk(t)表示串扰影响下系统的演化。则UXtalk(t)可写作串扰强度λ的表达式：
Figure BDA0003003841120000275
记
Figure BDA0003003841120000281
则
Figure BDA0003003841120000282
由上述结论可知，当
Figure BDA0003003841120000283
均为0时，
Figure BDA0003003841120000284
Figure BDA0003003841120000285
即彻底消除了串扰的影响。在本申请的一种可能的实现方案中，可以选择消除λ前n阶的影响，此时，可以将n阶串扰抑制问题总结为以下多目标优化问题：
Figure BDA0003003841120000286
subject to U(T)＝UTarget(T)；
其中，T表示操作持续的时间，UTarget(T)表示要实现的目标操作。
本方案进一步将其转化为单标量目标多参数优化问题，比如，可以采用梯度下降方法求解。该过程如下：
首先，为波形Ω(t)选取一种带参数的函数形式。例如，选用N个不同频率的余弦函数之和的形式：
Figure BDA0003003841120000287
其中，A＝(A1,A2,…,AN)为待优化参数，T为波形总时间。
其次，构造如下的单标量目标：
Figure BDA0003003841120000288
其中，F(U,V)表示操作U和V之间的平均操作保真度(AverageGateFidelity)，
Figure BDA0003003841120000289
为矩阵X的F-范数(Frobenius-范数)，wi为权重系数，通常取值在10-3～10-5。
在求解时，本方案可以通过求解薛定谔方程获得每个时刻的UCtrl(t)，进而采用数值积分的方法求得
Figure BDA00030038411200002810
最终求得L。通过数值求导的方法求得
Figure BDA00030038411200002811
采用梯度下降方法更新参数A。
4)由多个单/两比特操作组成的区域。
当区域由多个单/两比特操作组成时，类似上面的情形，区域的有效哈密顿量可以写作：
H(t)＝HCtrl(t)+λHXtalk+λ′HInner-Xtalk；
其中，λ表示跨区域的ZZ串扰的等效强度，λ′表示区域内ZZ串扰的等效强度。
在这种情况下，为了抑制跨区域的串扰，无需重新优化波形，只需要将按照上面方法优化出的波形应用在对应的操作上即可。例如，请参考图15，其示出了一种单比特操作和两比特操作的示意图。其中，图15展示了一个由一个单比特操作(3)和一个两比特操作(1、2)组成的区域及区域周围的耦合。该区域的有效哈密顿量可写作：
Figure BDA0003003841120000291
对比之前情形的哈密顿量，可以得到，在量子比特3上应用优化的单比特操作波形
Figure BDA0003003841120000292
可以抑制
Figure BDA0003003841120000293
串扰，在1、2上应用优化的两比特操作波形，可以抑制1-a、1-b、2-c、2-e串扰，从而抑制所有跨区域的串扰。也就是说，独立地应用优化过的波形即可抑制所有跨区域的串扰，这一特点保证了本申请所示方案的可扩展性：只需对只有一个单/两比特操作的区域进行优化，即可扩展到任意形状的区域上。
本申请后续内容通过计算机模拟仿真实验来说明本申请所示方案带来的改进。首先对比分析相关技术中的波形和使用本方案涉及的优化后的波形对ZZ串扰的抑制效果，然后对比分析我们的方案在量子计算任务上对相关技术方案的提升。
一、波形对ZZ串扰的抑制效果。
1)单比特波形。
相关技术所示的方案通常采用高斯波形实现单比特操作，本方案采用N项余弦波形叠加实现单比特操作，取N＝5。
为了刻画单比特波形对ZZ串扰的抑制效果，请参考图16，其示出了本申请实施例涉及的一种单比特波形拓扑图。在图16所示的拓扑中，在量子比特1上施加波形实现特定操作，如果ZZ串扰被抑制，则量子比特2的状态在施加波形前后应该保持一致。
本申请中，可采用施加波形后量子比特2的状态与施加波形前状态的差异刻画波形对ZZ串扰的抑制效果。这里采用不保真度来刻画状态差异，两个状态的不保真度定义为：
Figure BDA0003003841120000301
其中，<α|β>为|α>和|β>的内积。
请参考图17，其示出了本申请实施例涉及的串扰抑制情况示意图，请参考图18，其示出了本申请实施例涉及的另一种串扰抑制情况示意图。其中，图17展示了在量子比特1上实现Rx(π/2)操作时的串扰抑制情况，图18展示了实现恒等I操作时的串扰抑制情况。横轴表示不同的串扰强度，纵轴表示不保真度衡量的串扰抑制效果，不保真度越低，串扰抑制效果越好。
通过图17和图18的对比分析可得到如下结果：
相比于高斯波形，本方案涉及的优化后的波形对ZZ串扰的抑制效果提升了4～5个数量级。
本申请所示的方案涉及的优化后的波形对不同强度的串扰均存在强抑制效果。例如，以不保真度10-8为标准，按照本方案优化的Rx(π/2)能够抑制6MHz以内的串扰，按照本方案优化的I波形能够抑制8MHz以内的串扰。考虑到实验上观测到的ZZ串扰通常小于1MHz，可以说本方案实现了对ZZ串扰的强抑制。
2)两比特波形。
相关技术中的方案通常采用平顶高斯波形实现两比特操作，本方案考虑
Figure BDA0003003841120000302
类型的耦合，并采用直接在耦合项上施加波形的方式实现两比特操作。本方案采用N项余弦波形叠加，取N＝5。
为了刻画两比特波形对ZZ串扰的抑制效果，请参考图19，其示出了本申请实施例涉及的一种两比特波形拓扑图。在图19所示的拓扑中，在量子比特2、3上施加波形实现特定操作，如果ZZ串扰被抑制，则量子比特1、4的状态在施加波形前后应该保持一致。因而，采用施加波形前后量子比特1、4的状态不保真度来刻画抑制效果。
请参考图20，其示出了实现Rxx(π/2)操作时对ZZ串扰的抑制效果的对比示意图。通过图20的分析可得到与单比特波形相似的结果：
相比于平顶高斯波形，本方案涉及的优化后的波形对ZZ串扰的抑制效果提升了4～5个数量级。也就是说，本方案涉及的优化后的波形对不同强度的串扰均存在强抑制效果。
综合而言，本方案以一种简单、统一的方式，在单比特操作和两比特操作均实现了对ZZ串扰的强抑制。
二、量子计算任务上的提升。
本申请选择Hidden-Shift(HS)算法，Bernstein–Vazirani(BV)算法，量子傅立叶变换(QuantumFourierTransform，QFT)，Ising模型模拟，量子容量估计(QuantumVolume，QV)作为量子计算任务。请参考图21，其示出了本申请实施例涉及的网格拓扑。在图21所示的网格拓扑上进行计算机仿真模拟，拓扑图中相邻的两比特之间存在ZZ串扰，根据实验结果设定其强度为200kHz。
本申请提供量子操作调度的方案与最大化并行方式调度(ParallelSchedule)进行对比，本申请涉及的调度方式记为串扰抑制调度(X-suppressedSchedule)。请参考图22，其示出了两种调度方式的对比示意图。
在图22中，不同量子计算任务上，本申请提供的调度方式(X-suppressedSchedule)相比最大化并行方式调度(ParallelSchedule)带来明显的提升。
其中，图22中的横轴说明：HS-4表示4比特的HiddenShift算法，其他以此类推。
从图22中可得出如下结果：
本申请提供的方法相比最大化并行方式调度，大幅降低了ZZ串扰对量子计算任务的影响。在本申请评测的量子计算任务上，本申请提供的方法相比原先方法取得了10～110倍的提升。
本申请提供的方法具有良好扩展性：随着量子计算任务所用量子比特数目的增加，本申请提供的方法带来的提升也增加。本申请提供的方法不仅能适用于当下量子比特较少的量子计算器件，也能适用于未来大规模量子计算器件。
本申请涉及的波形优化方案对串扰HXtalk的形式没有要求，除了应用于ZZ串扰外，如果有其他类型串扰的话，只需要将HXtalk替换为相应的形式，本申请的方法亦可适用。类似地，上述的“最优抑制”算法和量子操作调度算法同样可适用于其他串扰。
本申请涉及的波形优化方法最终归结为求解一个多目标优化问题。上述实施例中给出了一种将其转为单标量目标多参数优化的求解方法，该方法中的梯度下降法可替换为其它优化方法，例如Nelder-Mead(NM)算法，Broyden-Fletcher-Goldfarb-Shanno(BFGS)算法等。除此之外，其他多目标优化算法也可用于求解该问题，例如遗传算法，模拟退火等。
本申请所示的方案可以用于在超导量子计算器件上执行量子电路时，抑制器件中存在的ZZ串扰对计算的影响，提高计算的保真度。请参考图23，其示出了本申请实施例涉及的执行量子电路的示意图。如图23所示，本申请涉及的量子操作调度方法可以用于处理量子电路，生成量子电路执行方案，本申请涉及的波形优化方法用于为执行方案中所涉及到的量子操作生成对应的控制波形；依据量子电路执行方案，使用控制波形控制量子计算器件，可实现在量子计算器件上执行量子电路的同时抑制ZZ串扰对计算的影响。
本申请涉及也可以应用于改进量子内存，抑制量子内存器件中存在的ZZ串扰，提高量子内存的数据保存时间。请参考图24和图25，其中图24示出了本申请实施例涉及的量子内存的示意图，图25示出了本申请实施例涉及的恒等量子电路示意图。如图24和25所示，量子内存的刷新操作等效于在量子内存上执行恒等量子电路。
综上所述，本申请实施例所示的方案，对于量子计算器件对应的量子操作，在保证量子操作不失真的情况下，获取使得量子操作在所述量子计算器件上执行时产生的区域间串扰最小化的控制波形，并基于该控制波形在量子计算器件上执行该量子操作，从而通过对量子操作的控制波形进行优化的方式来抑制该量子操作引起的区域间串扰，由于任意类型的量子操作都可以进行控制波形的优化，因此，该方案能够实现对任意类型量子操作的串扰抑制，提高了串扰抑制的适用性，从而提高串扰抑制的效果。
请参考图26，其示出了本申请一个实施例提供的量子操作执行装置的框图。该装置具有实现上述方法示例的功能，所述功能可以由硬件实现，也可以由硬件执行相应的软件实现。该装置可以是上文介绍的计算机设备，也可以设置在计算机设备中。如图26所示，该装置可以包括：
量子操作获取模块2601，用于获取量子电路对应的n个未调度的量子操作；n大于等于1，且n为整数；
量子操作调度模块2602，用于基于量子计算器件的拓扑结构，以及n个所述未调度的量子操作进行操作调度，获得目标操作组；所述目标操作组包含各个被调度的量子操作；
量子操作执行模块2603，用于在所述量子计算器件上并行执行所述目标操作组中包含的各个被调度的量子操作；
其中，所述被调度的量子操作的控制波形是基于指定目标生成的；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是所述量子计算器件中的各个量子比特以是否执行量子操作进行区分得到的。
在一种可能的实现方式中，所述量子操作调度模块2602，包括：
可调度集合获取子模块，用于基于n个所述未调度的量子操作获取可调度集合；所述可调度集合中包含n个所述未调度的量子操作中可调度的量子操作；
候选操作组提取子模块，用于从所述可调度集合中提取出候选操作组；
划分方案获取子模块，用于获取量子比特划分方案；所述量子比特划分方案是基于所述拓扑结构，按照指定划分方式得到的；所述量子比特划分方案指示将所述量子计算器件中的各个量子比特划分为两个量子比特集合，以在所述量子计算器件上形成至少两个所述连通区域；
目标操作组获取子模块，用于基于所述候选操作组以及所述量子比特划分方案，获取所述目标操作组。
在一种可能的实现方式中，所述装置还包括：
所述量子操作调度模块2602，还用于响应于在所述目标操作组之外，还存在m个所述未调度的量子操作，基于所述拓扑结构，以及m个所述未调度的量子操作进行操作调度，获得新的所述目标操作组；m大于等于1，且m为整数；
所述量子操作执行模块2603，还用于在所述量子计算器件上并行执行新的所述目标操作组中包含的各个被调度的量子操作。
在一种可能的实现方式中，所述候选操作组提取子模块，包括：
两比特操作集合获取单元，用于响应于所述可调度集合中包含两比特量子操作，获取所述可调度集合中的各个两比特量子操作，得到两比特操作集合；
操作组构建单元，用于响应于所述两比特操作集合中包含3个或者3个以上操作，基于所述两比特操作集合中，距离最近的两个量子操作构建两个操作组；
操作添加单元，用于基于所述两比特操作集合中，除了所述两个操作之外的其它量子操作与所述两个操作组之间的距离，执行i次操作添加，i大于等于1，且i为整数；
候选操作组获取单元，用于基于所述i次操作添加的结果，获取所述候选操作组。
在一种可能的实现方式中，所述操作添加单元，用于，
在第j次操作添加过程中，获取所述两比特操作集合中的所述其它量子操作分别与所述两个操作组之间的距离；1≤j≤i，且j为整数；
将最大距离对应的指定量子操作，添加入所述最大距离对应的指定操作组；所述最大距离是所述其它量子操作分别与所述两个操作组之间的距离中的最大值；
基于所述拓扑结构以及所述指定操作组，按照第一指定划分方式获得目标奇度顶点配对方案；所述目标奇度顶点配对方案包括目标对偶图中的各个奇度顶点对之间的路径，且所述目标奇度顶点配对方案用于指示对应的操作组的所述量子比特划分方案；所述目标对偶图是拓扑图的对偶图，所述拓扑图是所述拓扑结构的平面图；
响应于所述目标奇度顶点配对方案满足串扰抑制要求，将所述指定量子操作从所述两比特操作集合中删除；
响应于删除所述指定量子操作之后的所述其它量子操作的数量不为0，执行下一次操作添加过程。
在一种可能的实现方式中，所述候选操作组获取单元，用于响应于所述第i次操作添加后的所述其它量子操作的数量为0，或者，响应于所述第i次操作添加过程中得到的所述目标奇度顶点配对方案不满足所述串扰抑制要求，将所述第i次操作添加后的所述两个操作组中，包含操作数量最多的操作组获取为所述候选操作组。
在一种可能的实现方式中，在基于所述拓扑结构以及所述指定操作组，按照第一指定划分方式获得目标奇度顶点配对方案时，所述操作添加单元，用于，
在所述目标对偶图中删除所述指定操作组中的各个量子操作对应的边；
获取以所述目标对偶图中的奇度顶点为顶点的完全图；
在所述完全图中确定各个所述奇度顶点对；
获取各个所述奇度顶点对的至少一种奇度顶点配对方案；所述奇度顶点配对方案中包含各个所述奇度顶点对分别在的所述目标对偶图中的一条路径；
在至少一种所述奇度顶点配对方案中添加所述指定操作组中的各个量子操作对应在所述目标对偶图中的边；
从至少一种所述奇度顶点配对方案中，筛除所述两比特操作对应的量子比特不属于同一个连通区域的奇度顶点配对方案，获得筛除后的各个所述奇度顶点配对方案；
获取筛除后的各个所述奇度顶点配对方案各自的最大区域量子比特数，以及，筛除后的各个所述奇度顶点配对方案各自的不可抑制的串扰总数；所述不可抑制的串扰总数由所述奇度顶点配对方案中的路径对应在所述拓扑图中的边数进行表征；所述最大区域量子比特数由所述拓扑图除去割集之后的最大连通区域中的顶点数进行表征，所述割集是所述拓扑图中除去所述奇度顶点配对方案中的路径对应在所述拓扑图中的边之后剩余的边组成的集合；
基于筛除后的各个所述奇度顶点配对方案各自的最大区域量子比特数，以及，筛除后的各个所述奇度顶点配对方案各自的不可抑制的串扰总数，获取筛除后的各个所述奇度顶点配对方案的串扰抑制指数；
将对应所述串扰抑制指数最小的所述奇度顶点配对方案，获取为所述目标奇度顶点配对方案。
在一种可能的实现方式中，所述串扰抑制要求包括：
所述目标奇度顶点配对方案的最大区域量子比特数小于比特数阈值；
以及，
所述目标奇度顶点配对方案的不可抑制的串扰总数小于串扰总数阈值。
在一种可能的实现方式中，在第j次操作添加过程中，获取所述两比特操作集合中的所述其它量子操作分别与所述两个操作组之间的距离时，所述操作添加单元，用于，
在第j次操作添加过程中，获取第一量子操作与第一操作组中已有的各个量子操作之间的操作距离；所述操作距离用于指示两个两比特操作对应的量子比特在所述拓扑结构中的最短路径长度之和；所述第一量子操作是所述其它量子操作中的任意一个，所述第一操作组是所述两个操作组中的任意一个；
基于所述第一量子操作与所述第一操作组中已有的各个量子操作之间的操作距离，获取所述第一量子操作与所述第一操作组之间的距离。
在一种可能的实现方式中，所述候选操作组提取子模块，用于响应于所述可调度集合中的各个量子操作为单比特量子操作，将所述可调度集合获取为所述候选操作组；
所述划分方案获取子模块，包括：
配对方案获取单元，用于基于所述拓扑结构，按照第二指定划分方式获得目标奇度顶点配对方案；
去边单元，用于在所述拓扑图中去除所述目标奇度顶点配对方案对应的边，获得去边后的所述拓扑图；
着色单元，用于对去边后的所述拓扑图进行顶点二着色，获得量子比特划分方案。
在一种可能的实现方式中，所述配对方案获取单元，用于，
获取完全图，所述完全图中的顶点为目标对偶图中的奇度顶点；所述目标对偶图是拓扑图的对偶图，所述拓扑图是所述量子计算器件的拓扑结构的平面图；
在所述完全图中确定配对的各个奇度顶点对；
获取各个所述奇度顶点对的至少一种奇度顶点配对方案；所述奇度顶点配对方案中包含各个所述奇度顶点对分别在的所述目标对偶图中的一条路径；
获取至少一种所述奇度顶点配对方案各自的最大区域量子比特数，以及，至少一种所述奇度顶点配对方案各自的不可抑制的串扰总数；所述不可抑制的串扰总数由所述奇度顶点配对方案中的路径对应在所述拓扑图中的边数进行表征；所述最大区域量子比特数由所述拓扑图除去割集之后的最大连通区域中的顶点数进行表征，所述割集是所述拓扑图中除去所述奇度顶点配对方案中的路径对应在所述拓扑图中的边之后剩余的边组成的集合；
基于至少一种所述奇度顶点配对方案各自的最大区域量子比特数，以及，至少一种所述奇度顶点配对方案各自的不可抑制的串扰总数，获取至少一种所述奇度顶点配对方案的串扰抑制指数；
将对应所述串扰抑制指数最小的所述奇度顶点配对方案，获取为目标奇度顶点配对方案。
在一种可能的实现方式中，在获取各个所述奇度顶点对的至少一种奇度顶点配对方案时，所述操作添加单元或者所述配对方案获取单元，用于，
获取各个所述奇度顶点对各自在所述完全图中的前k短路径；k≥1，且k为整数；
对各个所述奇度顶点对各自在所述完全图中的前k短路径进行组合，获得至少一种所述奇度顶点配对方案。
在一种可能的实现方式中，所述串扰抑制指数为：αNQ+NC；
其中，NQ表示所述最大区域量子比特数，NC表示所述不可抑制的串扰总数，α表示NQ和NC的重要程度。
在一种可能的实现方式中，所述目标操作组获取子模块，用于，
基于所述量子比特划分方案，将所述量子计算器件中的各个量子比特划分为两个量子比特集合；
将所述候选操作组划分为所述两个量子比特集合分别对应的子操作组；所述子操作组中的量子操作对应的量子比特，处于对应的所述量子比特集合内；
将目标子操作组中的量子操作以及填充操作，合并为所述目标操作组；所述目标子操作组是所述两个量子比特集合分别对应的所述子操作组中，包含量子操作数量最多的子操作组；所述填充操作是填充量子比特对应的恒等量子操作，所述填充量子比特是所述目标子操作组对应的量子比特集合中，除了所述目标子操作组对应的量子比特之外的其它量子比特。
在一种可能的实现方式中，所述指定目标是基于对应的量子操作在控制波形作用下的平均操作保真度，以及所述区域间串扰获得的单标量目标。
在一种可能的实现方式中，所述单标量目标的公式为：
Figure BDA0003003841120000371
其中，L为所述单标量目标，U(T)表示对应的量子操作在控制波形作用下的的酉变换，UTarget(T)表示对应的量子操作的目标酉变换，F(U(T),UTarget(T))表示U(T)与UTarget(T)之间的平均操作保真度；
Figure BDA0003003841120000372
表示在所述区域间串扰影响下的封闭量子系统的演化。
在一种可能的实现方式中，所述装置还包括：
迭代更新模块，用于基于所述单标量目标，对目标量子操作的控制波形函数的函数参数进行R次迭代更新；R≥2且R为整数；所述目标量子操作是各个所述可调度量子操作以及所述恒等量子操作中的任意一个；
波形生成模块，用于基于迭代更新后的所述函数参数生成所述目标量子操作的控制波形。
在一种可能的实现方式中，所述迭代更新模块，用于，
获取第r次迭代过程中的第一系统演化，所述第一系统演化是在第r个控制波形影响下的封闭量子系统的演化；2≤r≤R，且r为整数；所述第r个控制波形是第r-1次迭代更新后的所述函数参数对应的控制波形；
基于所述第一系统演化获取第二系统演化，所述第二系统演化是第r次迭代过程中的所述区域间串扰影响下的封闭量子系统的演化；
基于所述第二系统演化，以及所述目标量子操作在所述第r个控制波形作用下的平均操作保真度，获取第r次迭代过程中的所述单标量目标；
基于第r次迭代过程中的所述单标量目标对第r-1次迭代更新后的所述函数参数进行更新，获得第r次迭代更新后的所述函数参数。
在一种可能的实现方式中，所述区域间串扰包括区域间ZZ串扰。
综上所述，本申请实施例所示的方案，对于量子计算器件对应的量子操作，在保证量子操作不失真的情况下，获取使得量子操作在所述量子计算器件上执行时产生的区域间串扰最小化的控制波形，并基于该控制波形在量子计算器件上执行该量子操作，从而通过对量子操作的控制波形进行优化的方式来抑制该量子操作引起的区域间串扰，由于任意类型的量子操作都可以进行控制波形的优化，因此，该方案能够实现对任意类型量子操作的串扰抑制，提高了串扰抑制的适用性，从而提高串扰抑制的效果。
请参考图27，其示出了本申请一个实施例提供的量子操作的控制波形生成装置的框图。如图27所示，该装置可以包括：
量子操作获取模块2701，用于获取被调度的量子操作；
参数获取模块2702，用于基于指定目标，获取所述被调度的量子操作的控制波形函数的函数参数；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是以是否执行量子操作，对所述量子计算器件中的各个量子比特进行划分得到的；
波形生成模块2703，用于基于所述函数参数生成所述被调度的量子操作的控制波形。
其中，上述量子操作的控制波形生成装置生成被调度的量子操作的控制波形的过程，可以参考上述图3所示实施例中的相关介绍，此处不再赘述。
综上所述，本申请实施例所示的方案，对于量子计算器件对应的量子操作，在保证量子操作不失真的情况下，获取使得量子操作在所述量子计算器件上执行时产生的区域间串扰最小化的控制波形，从而通过对量子操作的控制波形进行优化的方式来抑制该量子操作引起的区域间串扰，由于任意类型的量子操作都可以进行控制波形的优化，因此，该方案能够实现对任意类型量子操作的串扰抑制，提高了串扰抑制的适用性，从而提高串扰抑制的效果。
需要说明的是，上述实施例提供的装置，在实现其功能时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的装置与方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
请参考图28，其示出了本申请一个实施例提供的计算机设备的结构框图。该计算机设备可以用于实施上述实施例中提供的方法。以该计算机设备为经典计算机为例：
该计算机设备2800包括处理单元(如CPU(Central Processing Unit，中央处理器)、GPU(Graphics Processing Unit，图形处理器)和FPGA(Field Programmable GateArray，现场可编程逻辑门阵列)等)2801、包括RAM(Random-Access Memory，随机存储器)2802和ROM(Read-Only Memory，只读存储器)2803的系统存储器2804，以及连接系统存储器2804和中央处理单元2801的系统总线2805。该计算机设备2800还包括帮助服务器内的各个器件之间传输信息的基本输入/输出系统(Input Output System，I/O系统)2806，和用于存储操作系统2813、应用程序2814和其他程序模块2815的大容量存储设备2807。
该基本输入/输出系统2806包括有用于显示信息的显示器2808和用于用户输入信息的诸如鼠标、键盘之类的输入设备2809。其中，该显示器2808和输入设备2809都通过连接到系统总线2805的输入输出控制器2810连接到中央处理单元2801。该基本输入/输出系统2806还可以包括输入输出控制器2810以用于接收和处理来自键盘、鼠标、或电子触控笔等多个其他设备的输入。类似地，输入输出控制器2810还提供输出到显示屏、打印机或其他类型的输出设备。
该大容量存储设备2807通过连接到系统总线2805的大容量存储控制器(未示出)连接到中央处理单元2801。该大容量存储设备2807及其相关联的计算机可读介质为计算机设备2800提供非易失性存储。也就是说，该大容量存储设备2807可以包括诸如硬盘或者CD-ROM(Compact Disc Read-Only Memory，只读光盘)驱动器之类的计算机可读介质(未示出)。
不失一般性，该计算机可读介质可以包括计算机存储介质和通信介质。计算机存储介质包括以用于存储诸如计算机可读指令、数据结构、程序模块或其他数据等信息的任何方法或技术实现的易失性和非易失性、可移动和不可移动介质。计算机存储介质包括RAM、ROM、EPROM(Erasable Programmable Read-Only Memory，可擦写可编程只读存储器)、EEPROM(Electrically Erasable Programmable Read-Only Memory，电可擦写可编程只读存储器)、闪存或其他固态存储其技术，CD-ROM、DVD(Digital Video Disc，高密度数字视频光盘)或其他光学存储、磁带盒、磁带、磁盘存储或其他磁性存储设备。当然，本领域技术人员可知该计算机存储介质不局限于上述几种。上述的系统存储器2804和大容量存储设备2807可以统称为存储器。
根据本申请实施例，该计算机设备2800还可以通过诸如因特网等网络连接到网络上的远程计算机运行。也即计算机设备2800可以通过连接在该系统总线2805上的网络接口单元2811连接到网络2812，或者说，也可以使用网络接口单元2811来连接到其他类型的网络或远程计算机系统(未示出)。
所述存储器还包括至少一条指令、至少一段程序、代码集或指令集，该至少一条指令、至少一段程序、代码集或指令集存储于存储器中，且经配置以由一个或者一个以上处理器执行，以实现上述量子操作执行方法。
在一个示例性实施例中，还提供了一种计算机可读存储介质，所述存储介质中存储有至少一条指令、至少一段程序、代码集或指令集，所述至少一条指令、所述至少一段程序、所述代码集或指令集在被处理器执行时以实现上述量子操作执行方法。
可选地，该计算机可读存储介质可以包括：ROM(Read-Only Memory，只读存储器)、RAM(Random-Access Memory，随机存储器)、SSD(Solid State Drives，固态硬盘)或光盘等。其中，随机存取记忆体可以包括ReRAM(Resistance Random Access Memory，电阻式随机存取记忆体)和DRAM(Dynamic Random Access Memory，动态随机存取存储器)。
在一个示例性实施例中，还提供了一种计算机程序产品或计算机程序，所述计算机程序产品或计算机程序包括计算机指令，所述计算机指令存储在计算机可读存储介质中。计算机设备的处理器从所述计算机可读存储介质中读取所述计算机指令，所述处理器执行所述计算机指令，使得所述计算机设备执行上述实施例所示的方法。
在一个示例性实施例中，还提供了一种量子操作芯片，该量子操作芯片用于实现如上所述实施例所示的方法。
应当理解的是，在本文中提及的“多个”是指两个或两个以上。“和/或”，描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。另外，本文中描述的步骤编号，仅示例性示出了步骤间的一种可能的执行先后顺序，在一些其它实施例中，上述步骤也可以不按照编号顺序来执行，如两个不同编号的步骤同时执行，或者两个不同编号的步骤按照与图示相反的顺序执行，本申请实施例对此不作限定。
以上所述仅为本申请的示例性实施例，并不用以限制本申请，凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (15)

1.一种量子操作执行方法，其特征在于，所述方法包括：
获取量子电路对应的n个未调度的量子操作；n大于等于1，且n为整数；
基于量子计算器件的拓扑结构，以及n个所述未调度的量子操作进行操作调度，获得目标操作组；所述目标操作组包含各个被调度的量子操作；
在所述量子计算器件上并行执行所述目标操作组中包含的各个被调度的量子操作；
其中，所述被调度的量子操作的控制波形是基于指定目标生成的；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是以是否执行量子操作，对所述量子计算器件中的各个量子比特进行划分得到的。
2.根据权利要求1所述的方法，其特征在于，所述基于量子计算器件的拓扑结构，以及n个所述未调度的量子操作进行操作调度，获得目标操作组，包括：
基于n个所述未调度的量子操作获取可调度集合；所述可调度集合中包含n个所述未调度的量子操作中可调度的量子操作；
从所述可调度集合中提取出候选操作组；
获取量子比特划分方案；所述量子比特划分方案是基于所述拓扑结构，按照指定划分方式得到的；所述量子比特划分方案指示将所述量子计算器件中的各个量子比特划分为两个量子比特集合，以在所述量子计算器件上形成至少两个所述连通区域；
基于所述候选操作组以及所述量子比特划分方案，获取所述目标操作组。
3.根据权利要求2所述的方法，其特征在于，所述从所述可调度集合中提取出候选操作组，包括：
响应于所述可调度集合中包含两比特量子操作，获取所述可调度集合中的各个两比特量子操作，得到两比特操作集合；
响应于所述两比特操作集合中包含3个或者3个以上操作，基于所述两比特操作集合中，距离最近的两个量子操作构建两个操作组；
基于所述两比特操作集合中，除了所述两个操作之外的其它量子操作与所述两个操作组之间的距离，执行i次操作添加，i大于等于1，且i为整数；
基于所述i次操作添加的结果，获取所述候选操作组。
4.根据权利要求3所述的方法，其特征在于，所述基于所述两比特操作集合中，除了所述两个操作之外的其它量子操作与所述两个操作组之间的距离，执行i次操作添加，包括：
在第j次操作添加过程中，获取所述两比特操作集合中的所述其它量子操作分别与所述两个操作组之间的距离；1≤j≤i，且j为整数；
将最大距离对应的指定量子操作，添加入所述最大距离对应的指定操作组；所述最大距离是所述其它量子操作分别与所述两个操作组之间的距离中的最大值；
基于所述拓扑结构以及所述指定操作组，按照第一指定划分方式获得目标奇度顶点配对方案；所述目标奇度顶点配对方案包括目标对偶图中的各个奇度顶点对之间的路径，且所述目标奇度顶点配对方案用于指示对应的操作组的所述量子比特划分方案；所述目标对偶图是拓扑图的对偶图，所述拓扑图是所述拓扑结构的平面图；
响应于所述目标奇度顶点配对方案满足串扰抑制要求，将所述指定量子操作从所述两比特操作集合中删除；
响应于删除所述指定量子操作之后的所述其它量子操作的数量不为0，执行下一次操作添加过程。
5.根据权利要求4所述的方法，其特征在于，所述基于所述拓扑结构以及所述指定操作组，按照第一指定划分方式获得目标奇度顶点配对方案，包括：
在所述目标对偶图中删除所述指定操作组中的各个量子操作对应的边；
获取以所述目标对偶图中的奇度顶点为顶点的完全图；
在所述完全图中确定各个所述奇度顶点对；
获取各个所述奇度顶点对的至少一种奇度顶点配对方案；所述奇度顶点配对方案中包含各个所述奇度顶点对分别在的所述目标对偶图中的一条路径；
在至少一种所述奇度顶点配对方案中添加所述指定操作组中的各个量子操作对应在所述目标对偶图中的边；
从至少一种所述奇度顶点配对方案中，筛除所述两比特操作对应的量子比特不属于同一个连通区域的奇度顶点配对方案，获得筛除后的各个所述奇度顶点配对方案；
获取筛除后的各个所述奇度顶点配对方案各自的最大区域量子比特数，以及，筛除后的各个所述奇度顶点配对方案各自的不可抑制的串扰总数；所述不可抑制的串扰总数由所述奇度顶点配对方案中的路径对应在所述拓扑图中的边数进行表征；所述最大区域量子比特数由所述拓扑图除去割集之后的最大连通区域中的顶点数进行表征，所述割集是所述拓扑图中除去所述奇度顶点配对方案中的路径对应在所述拓扑图中的边之后剩余的边组成的集合；
基于筛除后的各个所述奇度顶点配对方案各自的最大区域量子比特数，以及，筛除后的各个所述奇度顶点配对方案各自的不可抑制的串扰总数，获取筛除后的各个所述奇度顶点配对方案的串扰抑制指数；
将对应所述串扰抑制指数最小的所述奇度顶点配对方案，获取为所述目标奇度顶点配对方案。
6.根据权利要求2所述的方法，其特征在于，所述从所述可调度集合中提取出候选操作组，包括：
响应于所述可调度集合中的各个量子操作为单比特量子操作，将所述可调度集合获取为所述候选操作组；
所述获取量子比特划分方案，包括：
基于所述拓扑结构，按照第二指定划分方式获得目标奇度顶点配对方案；
在所述拓扑图中去除所述目标奇度顶点配对方案对应的边，获得去边后的所述拓扑图；
对去边后的所述拓扑图进行顶点二着色，获得所述量子比特划分方案。
7.根据权利要求6所述的方法，其特征在于，所述基于所述拓扑结构，按照第二指定划分方式获得目标奇度顶点配对方案，包括：
获取完全图，所述完全图中的顶点为目标对偶图中的奇度顶点；所述目标对偶图是拓扑图的对偶图，所述拓扑图是所述量子计算器件的拓扑结构的平面图；
在所述完全图中确定配对的各个奇度顶点对；
获取各个所述奇度顶点对的至少一种奇度顶点配对方案；所述奇度顶点配对方案中包含各个所述奇度顶点对分别在的所述目标对偶图中的一条路径；
获取至少一种所述奇度顶点配对方案各自的最大区域量子比特数，以及，至少一种所述奇度顶点配对方案各自的不可抑制的串扰总数；所述不可抑制的串扰总数由所述奇度顶点配对方案中的路径对应在所述拓扑图中的边数进行表征；所述最大区域量子比特数由所述拓扑图除去割集之后的最大连通区域中的顶点数进行表征，所述割集是所述拓扑图中除去所述奇度顶点配对方案中的路径对应在所述拓扑图中的边之后剩余的边组成的集合；
基于至少一种所述奇度顶点配对方案各自的最大区域量子比特数，以及，至少一种所述奇度顶点配对方案各自的不可抑制的串扰总数，获取至少一种所述奇度顶点配对方案的串扰抑制指数；
将对应所述串扰抑制指数最小的所述奇度顶点配对方案，获取为所述目标奇度顶点配对方案。
8.根据权利要求5或7所述的方法，其特征在于，所述串扰抑制指数为：αNQ+NC；
其中，NQ表示所述最大区域量子比特数，NC表示所述不可抑制的串扰总数，α表示NO和NC的重要程度。
9.根据权利要求2所述的方法，其特征在于，所述基于所述候选操作组以及所述量子比特划分方案，获取所述目标操作组，包括：
基于所述量子比特划分方案，将所述量子计算器件中的各个量子比特划分为两个量子比特集合；
将所述候选操作组划分为所述两个量子比特集合分别对应的子操作组；所述子操作组中的量子操作对应的量子比特，处于对应的所述量子比特集合内；
将目标子操作组中的量子操作以及填充操作，合并为所述目标操作组；所述目标子操作组是所述两个量子比特集合分别对应的所述子操作组中，包含量子操作数量最多的子操作组；所述填充操作是填充量子比特对应的恒等量子操作，所述填充量子比特是所述目标子操作组对应的量子比特集合中，除了所述目标子操作组对应的量子比特之外的其它量子比特。
10.根据权利要求1所述的方法，其特征在于，所述指定目标是基于对应的量子操作在控制波形作用下的平均操作保真度，以及所述区域间串扰获得的单标量目标。
11.根据权利要求10所述的方法，其特征在于，所述方法还包括：
基于所述单标量目标，对目标量子操作的控制波形函数的函数参数进行R次迭代更新；R≥2且R为整数；所述目标量子操作是各个所述可调度量子操作以及所述恒等量子操作中的任意一个；
基于迭代更新后的所述函数参数生成所述目标量子操作的控制波形。
12.根据权利要求11所述的方法，其特征在于，所述基于所述单标量目标，对目标量子操作的控制波形函数的函数参数进行R次迭代更新，包括：
获取第r次迭代过程中的第一系统演化，所述第一系统演化是在第r个控制波形影响下的封闭量子系统的演化；2≤r≤R，且r为整数；所述第r个控制波形是第r-1次迭代更新后的所述函数参数对应的控制波形；
基于所述第一系统演化获取第二系统演化，所述第二系统演化是第r次迭代过程中的所述区域间串扰影响下的封闭量子系统的演化；
基于所述第二系统演化，以及所述目标量子操作在所述第r个控制波形作用下的平均操作保真度，获取第r次迭代过程中的所述单标量目标；
基于第r次迭代过程中的所述单标量目标对第r-1次迭代更新后的所述函数参数进行更新，获得第r次迭代更新后的所述函数参数。
13.一种量子操作的控制波形生成方法，其特征在于，所述方法包括：
获取被调度的量子操作；
基于指定目标，获取所述被调度的量子操作的控制波形函数的函数参数；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是在量子计算器件中，以包含的量子比特是否执行量子操作进行划分的；
基于所述函数参数生成所述被调度的量子操作的控制波形。
14.一种量子操作执行装置，其特征在于，所述装置包括：
量子操作获取模块，用于获取量子电路对应的n个未调度的量子操作；n大于等于1，且n为整数；
量子操作调度模块，用于基于量子计算器件的拓扑结构，以及n个所述未调度的量子操作进行操作调度，获得目标操作组；所述目标操作组包含各个被调度的量子操作；
量子操作执行模块，用于在所述量子计算器件上并行执行所述目标操作组中包含的各个被调度的量子操作；
其中，所述被调度的量子操作的控制波形是基于指定目标生成的；所述指定目标用于在控制对应的量子操作不失真的情况下使区域间串扰最小化；所述区域间串扰是分属于不同的连通区域的量子比特之间的串扰；所述连通区域是所述量子计算器件中的各个量子比特以是否执行量子操作进行区分得到的。
15.一种量子操作芯片，其特征在于，所述量子操作芯片用于实现如权利要求1至12任一所述的量子操作执行方法，或者，所述量子操作芯片用于实现如权利要求13所述的量子操作的控制波形生成方法。
CN202110357168.8A 2021-04-01 2021-04-01 量子操作执行方法、装置及量子操作芯片 Active CN113033812B (zh)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN202110357168.8A CN113033812B (zh) 2021-04-01 2021-04-01 量子操作执行方法、装置及量子操作芯片
PCT/CN2021/117519 WO2022205781A1 (zh) 2021-04-01 2021-09-09 量子操作执行方法、装置、芯片、设备及存储介质
KR1020227005112A KR20220137614A (ko) 2021-04-01 2021-09-09 양자 연산 실행 방법 및 장치, 칩, 디바이스, 및 저장 매체
EP21839311.4A EP4089593A4 (en) 2021-04-01 2021-09-09 METHOD AND APPARATUS FOR PERFORMING QUANTUM OPERATION, AND CHIP, DEVICE AND STORAGE MEDIA
JP2022513857A JP7381723B2 (ja) 2021-04-01 2021-09-09 量子操作実行方法及び装置、量子操作の制御波形の生成方法及び装置、量子操作チップ、並びにコンピュータ装置及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110357168.8A CN113033812B (zh) 2021-04-01 2021-04-01 量子操作执行方法、装置及量子操作芯片

Publications (2)

Publication Number Publication Date
CN113033812A true CN113033812A (zh) 2021-06-25
CN113033812B CN113033812B (zh) 2022-03-18

Family

ID=76454105

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110357168.8A Active CN113033812B (zh) 2021-04-01 2021-04-01 量子操作执行方法、装置及量子操作芯片

Country Status (5)

Country Link
EP (1) EP4089593A4 (zh)
JP (1) JP7381723B2 (zh)
KR (1) KR20220137614A (zh)
CN (1) CN113033812B (zh)
WO (1) WO2022205781A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114418107A (zh) * 2022-01-17 2022-04-29 北京百度网讯科技有限公司 酉算子编译方法、计算设备、装置及存储介质
WO2022205781A1 (zh) * 2021-04-01 2022-10-06 腾讯科技（深圳）有限公司 量子操作执行方法、装置、芯片、设备及存储介质

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105912070A (zh) * 2016-04-08 2016-08-31 中国科学院物理研究所 用于量子比特操控的数字波形调整方法
US20160314407A1 (en) * 2008-03-24 2016-10-27 D-Wave Systems Inc. Systems, devices, and methods for analog processing
US20180114138A1 (en) * 2013-08-02 2018-04-26 University Of Maryland Fault-tolerant scalable modular quantum computer architecture with an enhanced control of multi-mode couplings between trapped ion qubits
CN109213603A (zh) * 2018-05-31 2019-01-15 合肥本源量子计算科技有限责任公司 一种用于对接量子计算机与用户的云平台及平台运行方法
US20190044044A1 (en) * 2018-02-15 2019-02-07 Intel Corporation Reducing crosstalk from flux bias lines in qubit devices
CN109800883A (zh) * 2019-01-25 2019-05-24 合肥本源量子计算科技有限责任公司 量子机器学习框架构建方法、装置及量子计算机
CN110942152A (zh) * 2019-10-29 2020-03-31 合肥本源量子计算科技有限责任公司 量子芯片系统、量子计算处理系统及电子设备
CN111091195A (zh) * 2019-12-24 2020-05-01 北京百度网讯科技有限公司 一种超导电路结构及超导量子芯片、超导量子计算机
CN111183432A (zh) * 2017-09-29 2020-05-19 国际商业机器公司 减少和/或减轻量子位门中的串扰
CN111461334A (zh) * 2020-03-30 2020-07-28 北京百度网讯科技有限公司 量子电路的处理方法、装置及设备
CN111967603A (zh) * 2020-09-01 2020-11-20 腾讯科技（深圳）有限公司 量子芯片、量子处理器及量子计算机
CN112529193A (zh) * 2020-12-04 2021-03-19 北京百度网讯科技有限公司 基于量子系统的数据处理方法及量子设备

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11194642B2 (en) 2018-11-29 2021-12-07 International Business Machines Corporation Noise and calibration adaptive compilation of quantum programs
US10657304B1 (en) * 2019-01-03 2020-05-19 International Business Machines Corporation Mapping logical qubits on a quantum circuit
CN109685216B (zh) * 2019-01-11 2020-12-29 清华大学 一种量子计算机
CN110085094A (zh) * 2019-04-16 2019-08-02 国仪量子(合肥)技术有限公司 量子计算和量子精密测量的教学设备
US11615333B2 (en) * 2019-06-24 2023-03-28 International Business Machines Corporation Quantum circuit topology selection based on frequency collisions between qubits
CN113033812B (zh) * 2021-04-01 2022-03-18 腾讯科技（深圳）有限公司 量子操作执行方法、装置及量子操作芯片

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160314407A1 (en) * 2008-03-24 2016-10-27 D-Wave Systems Inc. Systems, devices, and methods for analog processing
US20180114138A1 (en) * 2013-08-02 2018-04-26 University Of Maryland Fault-tolerant scalable modular quantum computer architecture with an enhanced control of multi-mode couplings between trapped ion qubits
CN105912070A (zh) * 2016-04-08 2016-08-31 中国科学院物理研究所 用于量子比特操控的数字波形调整方法
CN111183432A (zh) * 2017-09-29 2020-05-19 国际商业机器公司 减少和/或减轻量子位门中的串扰
US20190044044A1 (en) * 2018-02-15 2019-02-07 Intel Corporation Reducing crosstalk from flux bias lines in qubit devices
CN109213603A (zh) * 2018-05-31 2019-01-15 合肥本源量子计算科技有限责任公司 一种用于对接量子计算机与用户的云平台及平台运行方法
CN109800883A (zh) * 2019-01-25 2019-05-24 合肥本源量子计算科技有限责任公司 量子机器学习框架构建方法、装置及量子计算机
CN110942152A (zh) * 2019-10-29 2020-03-31 合肥本源量子计算科技有限责任公司 量子芯片系统、量子计算处理系统及电子设备
CN111091195A (zh) * 2019-12-24 2020-05-01 北京百度网讯科技有限公司 一种超导电路结构及超导量子芯片、超导量子计算机
CN111461334A (zh) * 2020-03-30 2020-07-28 北京百度网讯科技有限公司 量子电路的处理方法、装置及设备
CN111967603A (zh) * 2020-09-01 2020-11-20 腾讯科技（深圳）有限公司 量子芯片、量子处理器及量子计算机
CN112529193A (zh) * 2020-12-04 2021-03-19 北京百度网讯科技有限公司 基于量子系统的数据处理方法及量子设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A.D. PATTERSON 等: "Calibration of a Cross-Resonance Two-Qubit Gate Between Directly Coupled Transmons", 《PHYSICAL REVIEW APPLIED》 *
SHENGYU ZHANG 等: "Quantum algorithms for graph problems with cut queries", 《ARXIV》 *
储继 等: "超导量子线路中的两比特绝热门", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
罗均文 等: "基于波长-模式双复用的量子保密通信系统", 《光学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022205781A1 (zh) * 2021-04-01 2022-10-06 腾讯科技（深圳）有限公司 量子操作执行方法、装置、芯片、设备及存储介质
CN114418107A (zh) * 2022-01-17 2022-04-29 北京百度网讯科技有限公司 酉算子编译方法、计算设备、装置及存储介质

Also Published As

Publication number Publication date
KR20220137614A (ko) 2022-10-12
WO2022205781A1 (zh) 2022-10-06
JP7381723B2 (ja) 2023-11-15
EP4089593A1 (en) 2022-11-16
CN113033812B (zh) 2022-03-18
JP2023523109A (ja) 2023-06-02
EP4089593A4 (en) 2022-11-16

Similar Documents

Publication Publication Date Title
JP7186797B2 (ja) 量子計算のための方法及びシステム
US12039407B2 (en) Systems and methods for analog processing of problem graphs having arbitrary size and/or connectivity
Harding et al. Meta-parametric design
Cortese et al. Loading classical data into a quantum computer
Lobato et al. Multi-objective genetic algorithm for missing data imputation
US20230267358A1 (en) Distributed Quantum Computing Simulation Method and Apparatus
CN113033812B (zh) 量子操作执行方法、装置及量子操作芯片
CN114492782B (zh) 基于强化学习的神经网络的片上核心编译映射方法及装置
Van Dijk et al. Tagged BDDs: Combining reduction rules from different decision diagram types
CN111915011B (zh) 一种单振幅量子计算模拟方法
US9594863B2 (en) Method for determining by optimization a multi-core architecture
Sundaram et al. Distribution of quantum circuits over general quantum networks
Dellnitz Big data efficiency analysis: Improved algorithms for data envelopment analysis involving large datasets
CN111931939A (zh) 一种单振幅量子计算模拟方法
CA2947578C (en) Method and system for generating an embedding pattern used for solving a quadratic binary optimization problem
WO2020204741A1 (en) Device and methods for a quantum circuit simulator
Verhetsel et al. A 44-element mesh of Schneiders’ pyramid: Bounding the difficulty of hex-meshing problems
Singh et al. Optimizing graph processing on gpus using approximate computing: Poster
EP3903242B1 (en) Device and methods for a quantum circuit simulator
Sumer et al. Fast parallel and adaptive updates for dual-decomposition solvers
CN116167447B (zh) 量子电路处理方法、装置及电子设备
WO2024116608A1 (ja) 計算機システム及び情報処理方法
Barkalov et al. High performance computing for global optimization problems
Huang et al. Qubit Mapping: The Adaptive Divide-and-Conquer Approach
CN117391206A (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
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 40046375

Country of ref document: HK

GR01 Patent grant
GR01 Patent grant