CN115577790B - 哈密顿量模拟方法、装置、设备及存储介质 - Google Patents
哈密顿量模拟方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN115577790B CN115577790B CN202211196971.9A CN202211196971A CN115577790B CN 115577790 B CN115577790 B CN 115577790B CN 202211196971 A CN202211196971 A CN 202211196971A CN 115577790 B CN115577790 B CN 115577790B
- Authority
- CN
- China
- Prior art keywords
- target
- quantum
- register
- revolving door
- state
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 84
- 238000004088 simulation Methods 0.000 title claims abstract description 64
- 239000011159 matrix material Substances 0.000 claims abstract description 99
- 238000012549 training Methods 0.000 claims description 209
- 239000002096 quantum dot Substances 0.000 claims description 54
- 230000006870 function Effects 0.000 claims description 52
- 238000012545 processing Methods 0.000 claims description 28
- 230000009471 action Effects 0.000 claims description 26
- 238000004590 computer program Methods 0.000 claims description 15
- 230000021615 conjugation Effects 0.000 claims 4
- 230000000452 restraining effect Effects 0.000 claims 2
- 239000007787 solid Substances 0.000 description 15
- 238000010586 diagram Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010801 machine learning Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 4
- 239000011800 void material Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000003213 activating effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011478 gradient descent method Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000002547 new drug Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 101100001674 Emericella variicolor andI gene Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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
- 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/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
-
- 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/60—Quantum algorithms, e.g. based on quantum optimisation, quantum Fourier or Hadamard transforms
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)
- Logic Circuits (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
Abstract
本公开提供了哈密顿量模拟方法、装置、设备及存储介质,涉及计算机技术领域,尤其涉及量子计算领域。具体实现方案为:确定目标量子电路的子电路中目标可调参数的目标参数值;目标参数值为在目标演化时间t0下所对应的、且满足第一误差条件的参数值;目标量子电路包含有辅助寄存器、块编码寄存器和主寄存器,子电路作用于辅助寄存器;目标量子电路中还包含有目标受控酉门,目标受控酉门用于模拟目标哈密顿量所对应的目标块编码矩阵U;在目标可调参数为目标参数值、辅助寄存器的第一输入态为预设初态、块编码寄存器的第二输入态为预设初态,以及主寄存器的第三输入态为目标哈密顿量所对应的第一量子态的情况下,得到主寄存器的目标输出量子态。
Description
技术领域
本公开涉及计算机技术领域,尤其涉及量子计算、量子模拟领域。
背景技术
随着量子技术的高速发展,量子计算行业正逐步进入规模化和产业化。大量的资金被用于研究容错量子芯片和量子算法,并期望以此来展示量子计算的优越性。量子模拟(Quantum simulation)作为一个关于量子系统动态演化的重要方向,在物理、化学和材料科学等领域有着广泛的应用。
量子模拟是量子计算的一个核心应用,高效且实用的量子模拟方法在量子化学研发新药和电池方面有着很好的应用前景,原因在于量子模拟可以用于模拟微观世界量子系统的演化,从而帮助新材料的研发与化学分子性质的模拟等。同时,量子模拟也是量子机器学习中几个常见量子算法的核心子步骤,比如量子主成分分析(quantum principalcomponent analysis)和量子线性系统求解方案(quantum algorithm for linearsystems of equations)等。
发明内容
本公开提供了一种哈密顿量模拟方法、装置、设备以及存储介质。
根据本公开的一方面,提供了一种哈密顿量模拟方法,包括:
确定目标量子电路的子电路中目标可调参数的目标参数值;其中,所述目标参数值为在目标演化时间t0下所对应的、且满足第一误差条件的参数值;所述目标量子电路包含有辅助寄存器、块编码寄存器和主寄存器,所述子电路作用于所述辅助寄存器;所述目标量子电路中还包含有受控于所述辅助寄存器且作用于所述块编码寄存器和所述主寄存器的目标受控酉门,所述目标受控酉门用于模拟目标哈密顿量所对应的目标块编码矩阵U;所述目标受控酉门包括与目标块编码矩阵U等价的第一受控酉门,和与所述目标块编码矩阵U的共轭转置等价的第二受控酉门;所述目标块编码矩阵U所对应的量子比特的数量为m+n,所述n为所述目标哈密顿量所对应的量子比特的数量,所述m为所述块编码寄存器所包含的量子比特的数量;
在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,得到所述目标量子电路中所述主寄存器的目标输出量子态,其中,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间t0下所对应的目标量子态。
根据本公开的另一方面,提供了一种哈密顿量模拟装置,包括:
参数调整单元,用于确定目标量子电路的子电路中目标可调参数的目标参数值;其中,所述目标参数值为在目标演化时间t0下所对应的、且满足第一误差条件的参数值;所述目标量子电路包含有辅助寄存器、块编码寄存器和主寄存器,所述子电路作用于所述辅助寄存器;所述目标量子电路中还包含有受控于所述辅助寄存器且作用于所述块编码寄存器和所述主寄存器的目标受控酉门,所述目标受控酉门用于模拟目标哈密顿量所对应的目标块编码矩阵U;所述目标受控酉门包括与目标块编码矩阵U等价的第一受控酉门,和与所述目标块编码矩阵U的共轭转置等价的第二受控酉门;所述目标块编码矩阵U所对应的量子比特的数量为m+n,所述n为所述目标哈密顿量所对应的量子比特的数量,所述m为所述块编码寄存器所包含的量子比特的数量;
输出单元,用于在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,得到所述目标量子电路中所述主寄存器的目标输出量子态,其中,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间t0下所对应的目标量子态。
根据本公开的再一方面,提供了一种计算设备,包括:
至少一个量子处理单元;
存储器,耦合到所述至少一个QPU并用于存储可执行指令,
所述指令被所述至少一个量子处理单元执行,以使所述至少一个量子处理单元能够执行以上所述的方法;
或者,包括:
至少一个处理器;以及
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行以上所述的方法。
根据本公开的再一方面,提供了一种存储有计算机指令的非瞬时计算机可读存储介质,当至少一个量子处理单元执行时,所述计算机指令使得所述至少一个量子处理单元执行以上所述的方法;
或者,所述计算机指令用于使所述计算机执行以上所述的方法。
根据本公开的再一方面,提供了一种计算机程序产品,包括计算机程序,所述计算机程序在被至少一个量子处理单元执行时实现以上所述的方法;
或者所述计算机程序在被处理器执行时实现以上所述的方法。
如此,本公开方案在模拟哈密顿量问题上,提供了一种新颖的模拟方案,进一步地,本公开方案可以在近期的量子计算机上实现,因此,实用性强;此外,本公开方案还可以应用于大规模的哈密顿量,因此,还兼具扩展性。
应当理解,本部分所描述的内容并非旨在标识本公开的实施例的关键或重要特征,也不用于限制本公开的范围。本公开的其它特征将通过以下的说明书而变得容易理解。
附图说明
附图用于更好地理解本方案,不构成对本公开的限定。其中:
图1是根据本公开实施例哈密顿量模拟方法的实现流程示意图一;
图2(a)是根据本公开实施例哈密顿量模拟方法的实现流程示意图二;
图2(b)是根据本公开实施例哈密顿量模拟方法的实现流程示意图三;
图2(c)是根据本公开实施例预设参数化量子电路训练方法的实现流程图;
图3(a)至图3(f)是根据本公开实施例预设参数化量子电路的结构示意图;
图4(a)至图4(e)是根据本公开实施例目标量子电路的结构示意图;
图5是根据本公开实施例哈密顿量模拟方法在一具体实施例中的实现流程示意图;
图6是根据本公开实施例哈密顿量模拟装置的结构示意图;
图7是用来实现本公开实施例哈密顿量模拟方法的计算设备的框图。
具体实施方式
以下结合附图对本公开的示范性实施例做出说明,其中包括本公开实施例的各种细节以助于理解,应当将它们认为仅仅是示范性的。因此,本领域普通技术人员应当认识到,可以对这里描述的实施例做出各种改变和修改,而不会背离本公开的范围和精神。同样,为了清楚和简明,以下的描述中省略了对公知功能和结构的描述。
本文中术语“和/或”,仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。本文中术语“至少一种”表示多种中的任意一种或多种中的至少两种的任意组合,例如,包括A、B、C中的至少一种,可以表示包括从A、B和C构成的集合中选择的任意一个或多个元素。本文中术语“第一”、“第二”表示指代多个类似的技术用语并对其进行区分,并不是限定顺序的意思,或者限定只有两个的意思,例如,第一特征和第二特征,是指代有两类/两个特征,第一特征可以为一个或多个,第二特征也可以为一个或多个。
另外,为了更好的说明本公开,在下文的具体实施方式中给出了众多的具体细节。本领域技术人员应当理解,缺少某些具体细节,本公开同样可以实施。在一些实例中,对于本领域技术人员熟知的方法、手段、元件和电路未作详细描述,以便于凸显本公开的主旨。
随着量子技术的高速发展,量子计算行业正逐步进入规模化和产业化。大量的资金被用于研究容错量子芯片和量子算法,并期望以此来展示量子计算的优越性。量子模拟(Quantum simulation)作为一个关于量子系统动态演化的重要方向,在物理、化学和材料科学等领域有着广泛的应用。
量子模拟是量子计算的一个核心应用,高效且实用的量子模拟方法在量子化学研发新药和电池方面有着很好的应用前景,原因在于量子模拟可以用于模拟微观世界量子系统的演化,从而帮助新材料的研发与化学分子性质的模拟等。同时,量子模拟也是量子机器学习中几个常见量子算法的核心子步骤,比如量子主成分分析(quantum principalcomponent analysis)和量子线性系统求解方案(quantum algorithm for linearsystems of equations)等。
一般来说,量子模拟,也可称为哈密顿量模拟,是一项非常困难的任务。经典计算想要完成这样的任务,需要对目标量子系统的哈密顿量进行层析等操作,而且,对于指数增长的量子系统而言,更加困难。目前能够进行哈密顿量模拟的方案在量子电路宽度等方面也都有着较高的要求。因此,迫切需要更加高效与实用的哈密顿量模拟方案,如此,一方面便于研究量子系统,另一方面也可以促进量子计算在化学、机器学习等领域中的发展。
基于此,本公开方案提供了一种哈密顿量模拟方案,而且,能够高效得到演化后的目标量子态。
具体地,图1是根据本公开实施例哈密顿量模拟方法的实现流程示意图一;该方法可选地应用于兼具经典计算能力的量子计算设备中,也可以应用于兼具量子计算能力的经典计算设备中,或者,直接应用于经典计算设备,比如,个人电脑、服务器、服务器集群等具有经典计算能力的电子设备中,或者,直接应用于量子计算机中,本公开方案对此不作限制。
进一步地,该方法包括以下内容的至少部分内容。如图1所示,所述量子计算处理方法,包括:
步骤S101:确定目标量子电路的子电路中目标可调参数的目标参数值。
其中,所述目标参数值为在目标演化时间t0下所对应的、且满足第一误差条件的参数值;所述目标量子电路包含有辅助寄存器、块编码寄存器和主寄存器,所述子电路作用于所述辅助寄存器;所述目标量子电路中还包含有受控于所述辅助寄存器且作用于所述块编码寄存器和所述主寄存器的目标受控酉门,所述目标受控酉门用于模拟目标哈密顿量所对应的目标块编码矩阵U,如此,通过所述目标块编码矩阵U来模拟所述目标哈密顿量的演化。
进一步地,所述目标受控酉门包括与目标块编码矩阵U等价的第一受控酉门,和与所述目标块编码矩阵U的共轭转置等价的第二受控酉门;也即,所述第一受控酉门受控于辅助寄存器,且作用于所述块编码寄存器和所述主寄存器,同理,所述第二受控酉门受控于辅助寄存器,且作用于所述块编码寄存器和所述主寄存器。
这里,所述目标块编码矩阵U所对应的量子比特的数量为m+n,所述n为所述目标哈密顿量所对应的量子比特的数量,也即,所述目标哈密顿量为包含有n个量子比特的目标量子系统的哈密顿量。进一步地,该目标量子系统的量子态为第一量子态,也即所述目标哈密顿量对应所述第一量子态。进一步地,所述m为所述块编码寄存器所包含的量子比特的数量。
这里,所述m为大于等于1的正整数,所述n也为大于等于1的正整数。
可以理解的是,所述子电路中为所述目标量子电路中包含有目标可调参数的至少部分电路,也即该子电路为包含目标可调参数的参数化量子电路。
步骤S102:在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,得到所述目标量子电路中所述主寄存器的目标输出量子态,其中,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间t0下所对应的目标量子态。
这样,本公开方案采用包含辅助寄存器、块编码寄存器和主寄存器的目标量子电路,在目标可调参数为目标参数值的情况下,通过输入第一输入态、第二输入态和第三输入态,来模拟得到所述目标哈密顿量在所述目标演化时间t0下的目标量子态;如此,本公开方案在模拟目标哈密顿量问题上,提供了一种新颖的模拟方案。
进一步地,本公开方案可以在近期的量子计算机上实现,因此,实用性强;此外,本公开方案还可以应用于大规模的哈密顿量,因此,还兼具扩展性。
在一具体示例中,所述辅助寄存器包含有至少一个量子比特,比如,一个、或两个,或两个以上的量子比特。进一步地,所述主寄存器包含有的量子比特的数量与所述目标量子系统所包含的量子比特的数量相关,比如,所述主寄存器所包含的量子比特的数量=所述目标量子系统所包含的量子比特的数量n。
这里,为了便于区分,可将辅助寄存器所包含的量子比特称为辅助量子比特;同理,将所述块编码寄存器所包含的量子比特称为块编码量子比特;将所述主寄存器所包含的量子比特称为主量子比特。
在一具体示例中,所述预设初态可以具体比如|0>,或|1>。本公开方案对此不作具体限制。
图2(a)是根据本公开实施例哈密顿量模拟方法的实现流程示意图二。该方法可选地可以应用于兼具经典计算能力的量子计算设备中,也可以应用于兼具量子计算能力的经典计算设备中,或者,直接应用于经典计算设备,比如,个人电脑、服务器、服务器集群等具有经典计算能力的电子设备中,或者,直接应用于量子计算机中,本公开方案对此不作限制。
可以理解的是,以上图1所示方法的相关内容,也可以应用于该示例中,该示例对相关联内容不再赘述。
进一步地,该方法包括以下内容的至少部分内容。具体地,如图2(a)所示,该方法包括:
步骤S201a:确定目标量子电路的子电路中目标可调参数的目标参数值。
步骤S202a:在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,获取所述目标量子电路中所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息。
步骤S203a:在所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息均满足预设条件的情况下,得到所述主寄存器的目标输出量子态,所述目标输出量子态为模拟所述目标哈密顿量在所述演化时间t下所对应的目标量子态。
比如,在一具体示例中,在所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息均为|0>的情况下,可认为满足所述预设条件。
可以理解的是,实际应用中,当所述第一状态信息和第二状态信息为|0>的概率大于一阈值的情况下,即可认为两者均为|0>,此时,认为满足预设条件;相应地,主寄存器的输出态即为目标量子态。
这样,提供了一种新颖且便于实现的哈密顿量的模拟方案,而且,本公开方案可以在近期的量子计算机上实现,因此,实用性强;此外,本公开方案还可以应用于大规模的哈密顿量,因此,还兼具扩展性。
图2(b)是根据本公开实施例哈密顿量模拟方法的实现流程示意图三。该方法可选地可以应用于兼具经典计算能力的量子计算设备中,也可以应用于兼具量子计算能力的经典计算设备中,或者,直接应用于经典计算设备,比如,个人电脑、服务器、服务器集群等具有经典计算能力的电子设备中,或者,直接应用于量子计算机中,本公开方案对此不作限制。
可以理解的是,以上图1和图2(a)所示方法的相关内容,也可以应用于该示例中,该示例对相关联内容不再赘述。
进一步地,该方法包括以下内容的至少部分内容。具体地,如图2(b)所示,该方法包括:
步骤S201b:将训练完成的预设参数化量子电路中所述目标可调参数的目标参数值,作为所述子电路中目标可调参数的目标参数值。
也就是说,所述预设参数化量子电路中包含有所述目标可调参数,如此,将训练完成的预设参数化量子电路中所述目标可调参数的目标参数值,来作为该子电路中目标可调参数的目标参数值。换言之,该示例中,可以通过训练其他参数化量子电路,来得到该子电路中目标可调参数的目标参数值。
可以理解的是,该示例中,所述子电路和所述目标量子电路的相关说明,可参见以上描述,此处不作赘述。
需要说明的是,所述预设参数化量子电路中还可以包含其他可调参数,本公开方案对此不作具体限制,只要所述预设参数化量子电路包含有子电路所需的目标可调参数即可。
进一步地,所述训练完成的预设参数化量子电路用于模拟目标函数f(x)。
这里,所述目标函数f(x)用于表征演化时间t与自变量x之间关联关系。
进一步地,所述目标量子电路基于如下所得:
将所述预设参数化量子电路中的量子比特作为辅助寄存器,并扩展出块编码寄存器和主寄存器,同时,将所述预设参数化量子电路中作用于所述辅助寄存器的第一目标旋转门替换为所述第一受控酉门,使得所述第一受控酉门受控于所述辅助寄存器,且作用于所述块编码寄存器和所述主寄存器;以及将所述预设参数化量子电路中作用于所述辅助寄存器的第二目标旋转门替换为所述第二受控酉门,使得所述第二受控酉门受控于所述辅助寄存器,且作用于所述块编码寄存器和所述主寄存器。也就是说,所述目标量子电路是在预设参数化量子电路的基础上扩展所得。
这里,所述第一目标旋转门的第一旋转参数和第二目标旋转门的第二旋转参数均为所述目标函数f(x)的自变量x。
进一步地,所述子电路包含有所述预设参数化量子电路中除所述第一目标旋转门和第二目标旋转门之外的至少部分电路;这里,所述第一目标旋转门和第二目标旋转门可统称为目标旋转门,此时,所述子电路包含有所述预设参数化量子电路中除目标旋转门之外的至少部分电路。
可以理解的是,由于所述目标量子电路是在预设参数化量子电路的基础上扩展所得,所以,所述子电路也可以理解为是在预设参数化量子电路的基础上所得,而且,包含有所述预设参数化量子电路中所述目标可调参数对应的部分电路结构,如此,为通过训练该预设参数化量子电路来得到子电路的目标可调参数的目标参数值奠定了基础。
步骤S202b:在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,得到所述目标量子电路中所述主寄存器的目标输出量子态,其中,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间t0下所对应的目标量子态。
可以理解的是,该示例可采用如图2(b)所示的方式得到目标量子态,本公开方案在此不再赘述。
可以理解的是,由于所述预设参数化量子电路,相较于目标量子电路的电路结构简单,因此,通过训练预设参数化量子电路来得到目标可调参数的目标参数值的方式,能够有效减小计算量,为高效模拟哈密顿量奠定了基础。
进一步地,实际应用中,该预设参数化量子电路还可以在经典计算设备中通过模拟的方式得到,相应地,训练得到目标可调参数的目标参数值也可在经典计算设备中实现,所以,本公开方案得到目标可调参数的目标参数值的方式,可以不占用量子计算资源,因此,在为高效模拟哈密顿量奠定基础的同时,还有效降低了计算成本。
而且,本公开方案对第一量子态不作任何限制,换言之,能够模拟任意哈密顿量,通用性强。同时,本公开方案还可以在近期的量子计算机上实现,而且,不需要量子傅里叶变换,实用性强;此外,本公开方案还可以应用于大规模的哈密顿量,因此,还兼具扩展性。综上可知,本公开方案兼具高效性、实用性以及扩展性。
在本公开方案的一具体示例中,所述子电路中目标可调参数的目标参数值通过以下训练方式得到,也即采用如下方式对预设参数化量子电路(如下述第一种方式或第二种方式所构建的)进行训练,并训练得到所述目标可调参数的目标参数值;具体地,如图2(c)所示,所述方法还包括:
步骤S201c:在所述预设参数化量子电路的旋转参数x取值为N个数据点中的任一数据点xj的情况下,获取所述预设参数化量子电路的实际输出态。
这里,所述实际输出态为在预设初态作为输入态的情况下,所述预设参数化量子电路在所述目标可调参数为当前参数值下的输出态。
步骤S202c:基于所述实际输出态与所述预设初态得到实际输出结果yj。
这里,所述实际输出结果yj为所述实际输出态与输入态(也即预设初态)之间的内积;所述N为大于等于1的正整数,所述j为大于等于1小于等于N的正整数;所述旋转参数x包括所述第一旋转参数和所述第二旋转参数。
可以理解的是,在以下所述的图4(a)所示的结构中,不同层中目标旋转门所对应的旋转参数可统称为旋转参数。
步骤S203c:得到N个实际输出结果yj。
也就是说,在j取值为1至N的情况下,即可得到N个实际输出结果yj。
步骤S204c:确定是否满足迭代终止条件;在确定满足迭代终止条件的情况下,执行步骤S205c;否则,执行步骤S206c。
这里,所述迭代终止条件包括以下至少之一:
方式一:基于所述N个实际输出结果yj与N个目标输出结果确定预设损失函数的损失值满足收敛条件;所述目标输出结果
方式二:当前迭代次数到达预设次数。
实际应用中,只要满足以上条件之一,即可为满足迭代终止条件。
步骤S205c:将所述目标可调参数的当前参数值作为训练完成的所述预设参数化量子电路中所述目标可调参数的目标参数值。
步骤S206c:对所述目标可调参数的参数值进行调整,并返回至步骤S201c,以重新获取参数值调整后的实际输出态,以得到实际输出结果yj,进而得到N个实际输出结果yj,并重新确定是否满足迭代终止条件,直至满足迭代终止条件为止。
这里,需要说明的是,实际应用中,演化时间t为指定时间,比如为指定的目标演化时间t0,此时,通过上述训练即可得到该目标演化时间t0所对应的目标参数值。
这样,通过训练其他参数化量子电路,来得到子电路的目标可调参数的目标参数值;这里,由于所述预设参数化量子电路,相较于目标量子电路的电路结构简单,因此,通过训练预设参数化量子电路来得到目标可调参数的目标参数值的方式,能够有效减小计算量,为高效模拟哈密顿量奠定了基础。
在另一具体示例中,还可以采用函数分析方法来得到目标可调参数的目标参数值;具体地,得到目标函数的目标傅里叶级数F(x),其中,该目标傅里叶级数F(x)为目标定义域内近似所述目标函数的傅里叶级数。进一步地,基于所述目标傅里叶级数F(x)得到其他傅里叶级数,比如其他傅里叶级数P(x)和Q(x),其中,
基于预设关系式,即可得到目标可调参数的目标参数值;比如,对于图4(c)所示的目标量子电路而言,该预设关系式可具体为:
这里,所述Q*(x)为Q(x)的复数共轭,P*(x)为P(x)的复数共轭。
这样,能够有效减小计算量,为高效模拟哈密顿量奠定了基础。
可以理解的是,实际应用中,还可以使用任何能以一定精度近似目标函数的三角多项式,来优化得到目标可调参数的最优参数值,本公开方案对此不作具体限制。
以下给出两种方式来构建预设参数化量子电路,具体包括:
第一种方式:
该方式中,所述预设参数化量子电路包括有L个训练层;所述L为大于等于2的偶数,所述L的取值与所述第一误差条件有关;
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数x用于对第一角度进行旋转操作;
用于对第二角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第一旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第一旋转门的旋转角度φ和所述第二旋转门的旋转角度θ为所述目标可调参数。
这里,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;也就是说,将所述预设参数化量子电路中不同训练层的目标旋转门替换成不同的受控酉门,比如,将所述预设参数化量子电路中一训练层中的目标旋转门(为便于描述,可称为第一目标旋转门)替换成第一受控酉门,同时,将所述预设参数化量子电路中另一训练层中的目标旋转门(为便于描述,可称为第二目标旋转门)替换成第二受控酉门,如此,得到目标量子电路。
需要说明的是,实际应用中,所述L个训练层中不同的其他训练层所包含的旋转门种类和数量,可以相同,比如,均包括以上所述的旋转门;或者,也可以不相同,比如,一些其他训练层中包含有以上所述的旋转门中至少一个,另外一些训练层中还包含有其他量子门等,本公开方案对此不作限制,只要至少存在两个训练层包括以上所述量子门即可。
在一具体示例中,所述预设参数化量子电路中包含有一个量子比特,此时,所述目标旋转门、第一旋转门以及第二旋转门均为作用于该量子比特上的单量子比特旋转门。
进一步地,在另一示例中,所述预设参数化量子电路中包含有一个量子比特,而且,所述L个训练层中的每一训练层均包含有目标旋转门、第一旋转门以及所述第二旋转门,也即各训练层的目标旋转门、第一旋转门以及所述第二旋转门均为作用于该量子比特上的单量子比特旋转门。
第二种方式:
该方式中,所述预设参数化量子电路包括有L个训练层;所述L为大于等于2的偶数,所述L的取值与所述第一误差条件有关;
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数x用于对第一角度进行旋转操作;其中,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第二旋转门的旋转角度θ为所述目标可调参数。
也就是说,相较于第一种方式,该第二种方式中,所述至少两个训练层中的训练层中不包含有第一旋转门。可以理解的是,除第一旋转门外,上述第一种方式的有关描述,同样适用于该第二种方式,此处不再赘述。
这样,本公开方案有效提升了预设参数化量子电路的表达能力,同时,所使用的量子门的种类和数量少,为高效模拟哈密顿量奠定了基础。
进一步地,在本公开方案的一具体示例中,各角度满足如下条件之一:
所述第一角度为z轴所对应的角度;
所述第二角度为z轴所对应的角度;
所述第三角度为y轴所对应的角度。
也就是说,在一示例中,所述第一角度为z轴所对应的角度;在另一示例中,所述第二角度为z轴所对应的角度;在再一示例中,所述第三角度为y轴所对应的角度;或者,满足上述条件中的任意两个,如所述第一角度和第二角度均为z轴所对应的角度等。或者,同时满足上述三个条件,即所述第一角度和第二角度均为z轴所对应的角度,第三角度为y轴所对应的角度。
举例来说,在一具体示例中,所述L个训练层中的至少两个训练层包括:
所述目标旋转门,所述旋转参数x用于对z轴所对应的角度进行旋转操作;
用于对z轴所对应的角度进行旋转操作的所述第一旋转门;
用于对y轴所对应的角度进行旋转操作的所述第二旋转门。
或者,在另一示例中,所述L个训练层中的至少两个训练层包括:
所述目标旋转门,所述旋转参数x用于对z轴所对应的角度进行旋转操作;
用于对y轴所对应的角度进行旋转操作的所述第二旋转门。
进一步地,在另一具体示例中,所述预设参数化量子电路中包含有一个量子比特,此时,所述目标旋转门、第一旋转门以及第二旋转门均为作用于该量子比特上的单量子比特旋转门。
进一步地,所述L个训练层中的每一训练层均包含有:
所述目标旋转门,所述旋转参数x用于对z轴所对应的角度进行旋转操作;
用于对z轴所对应的角度进行旋转操作的所述第一旋转门;
用于对y轴所对应的角度进行旋转操作的所述第二旋转门。
或者,所述L个训练层中的每一训练层均包含有:
所述目标旋转门,所述旋转参数x用于对z轴所对应的角度进行旋转操作;
用于对y轴所对应的角度进行旋转操作的所述第二旋转门。
这样,本公开方案有效提升了预设参数化量子电路的表达能力,同时,所使用的量子门的种类和数量少,待训练的目标可调参数的数量少,如此,为高效模拟哈密顿量奠定了基础。
进一步地,在本公开方案的另一具体示例中,在所述L个训练层的任一训练层中包含有所述目标旋转门、所述第一旋转门以及所述第二旋转门的情况下,各旋转门的作用顺序为:
所述第一旋转门、第二旋转门和目标旋转门。
或者,在另一具体示例中,在所述L个训练层的任一训练层中包含有所述目标旋转门和所述第二旋转门的情况下,各旋转门的作用顺序为:第二旋转门和目标旋转门。
也就是说,在一具体示例中,所述L个训练层中的至少两个训练层中的每一训练层所包含的目标旋转门、第一旋转门和第二旋转门,按照旋转门的作用顺序,依次包括:
用于对z轴所对应的角度进行旋转操作的所述第一旋转门;
用于对y轴所对应的角度进行旋转操作的所述第二旋转门;
所述目标旋转门。
或者,在另一具体示例中,所述L个训练层中的至少两个训练层中的每一训练层所包含的目标旋转门和第二旋转门,按照旋转门的作用顺序,依次包括:
用于对y轴所对应的角度进行旋转操作的所述第二旋转门;
所述目标旋转门。
举例来说,以所述预设参数化量子电路中包含有一个量子比特,相对应,所述目标旋转门、第一旋转门以及第二旋转门均为作用于该量子比特上的单量子比特旋转门为例进行说明,如图3(a)所示,所述L个训练层的至少两个训练层中的一个训练层,比如所述L个训练层中的第i个训练层,按照作用顺序,依次包括:
旋转角度φi为z轴所对应的角度的第一旋转门RZ(φi);
旋转角度θi为y轴所对应的角度的第二旋转门RY(θi);
旋转参数xj为z轴所对应的角度的目标旋转门RZ(xj)。
这里,所述第一旋转门RZ(φi)的旋转角度φi和第二旋转门RY(θi)的旋转角度θi为第i个训练层中的目标可调参数,所述i为大于等于1小于等于L的整数。可以理解的是,该示例中,所述L个训练层的至少两个训练层中的另一个训练层,其结构也为如图3(a)所示的结构,这里不再赘述。
进一步地,在另一具体示例中,所述L个训练层中的每一训练层结构均为如图3(a)所示的结构,此处不再赘述。
再举例来说,以所述预设参数化量子电路中包含有一个量子比特,相对应,所述目标旋转门和第二旋转门均为作用于该量子比特上的单量子比特旋转门为例进行说明,如图3(d)所示,所述L个训练层的至少两个训练层中的一个训练层,比如所述L个训练层中的第i个训练层,按照作用顺序,依次包括:
旋转角度θi为y轴所对应的角度的第二旋转门RY(θi);
旋转参数xj为z轴所对应的角度的目标旋转门RZ(xj)。
这里,所述第二旋转门RY(θi)的旋转角度θi为第i个训练层中的目标可调参数,所述i为大于等于1小于等于L的整数。可以理解的是,该示例中,所述L个训练层的至少两个训练层中的另一个训练层,其结构也为如图3(d)所示的结构。这里不再赘述。
进一步地,在另一具体示例中,所述L个训练层中的每一训练层结构均为如图3(d)所示的结构,此处不再赘述。
这样,本公开方案有效提升了预设参数化量子电路的表达能力,同时,所使用的量子门的种类和数量少,待训练的目标可调参数的数量少,如此,为高效模拟哈密顿量奠定了基础。
进一步地,在另一具体示例中,所述预设参数化量子电路的L个训练层之后,还包括其他旋转门。
在一具体示例中,所述预设参数化量子电路的L个训练层之后,还包括:
用于对第四角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第三旋转门;
用于对第五角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第四旋转门;
其中,所述第三旋转门的旋转角度φ0和所述第四旋转门的旋转角度θ0为所述目标可调参数。
在一具体示例中,所述预设参数化量子电路中包含有一个量子比特,此时,所述第三旋转门和第四旋转门均为作用于该量子比特上的单量子比特旋转门。
举例来说,在一示例中,如图3(b)所示,所述预设参数化量子电路在L个训练层之后,还包括:
旋转角度φ0为z轴所对应的角度的第三旋转门RZ(φ0);
旋转角度θ0为y轴所对应的角度的第四旋转门RY(θ0)。
这里,所述旋转角度φ0和旋转角度θ0也为目标可调参数。
基于此,如图3(b)所示的预设参数化量子电路的数学表达式可具体为:
在一具体示例中,所述预设参数化量子电路的L个训练层之后,还包括其他旋转门:
用于对第四角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第三旋转门;
用于对第五角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第四旋转门;
用于对第六角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第五旋转门;
其中,所述第三旋转门的旋转角度φ0和所述第四旋转门的旋转角度θ0为所述目标可调参数;所述第五旋转门的旋转角度α为固定参数,即为不参与训练的参数。或者,所述第三旋转门的旋转角度φ0、所述第四旋转门的旋转角度θ0以及所述第五旋转门的旋转角度α均为所述目标可调参数。
在一具体示例中,所述预设参数化量子电路中包含有一个量子比特,此时,所述第三旋转门、第四旋转门和第五旋转门均为作用于该量子比特上的单量子比特旋转门。
举例来说,在另一示例中,如图3(c)所示,所述预设参数化量子电路在L个训练层之后,还包括:
旋转角度φ0为z轴所对应的角度的第三旋转门RZ(φ0);
旋转角度θ0为y轴所对应的角度的第四旋转门RY(θ0);
以及旋转角度α为z轴所对应的角度的第五旋转门RZ(α)。
这里,所述旋转角度φ0、旋转角度θ0和旋转角度α均为目标可调参数。
基于此,如图3(c)所示的预设参数化量子电路的数学表达式可具体为:
或者,所述旋转角度φ0和旋转角度θ0均为目标可调参数,而旋转角度α为固定参数,不参与训练。
基于此,如图3(c)所示的预设参数化量子电路的数学表达式可具体为:
或者,在另一示例中,所述预设参数化量子电路的L个训练层之后,还包括:
用于对第五角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第四旋转门;
其中,所述第四旋转门的旋转角度θ0为所述目标可调参数。
需要说明的是,所述第四旋转门的相关内容,可参照上述描述,此处不再赘述。也就是说,相较于如图3(b)所示的结构,该示例中,如图3(e)所示,L个训练层之后包括第四旋转门,不包括第三旋转门。
或者,在再一示例中,所述预设参数化量子电路的L个训练层之后,还包括其他旋转门:
用于对第五角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第四旋转门;
用于对第六角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第五旋转门;
其中,所述第四旋转门的旋转角度θ0为所述目标可调参数;所述第五旋转门的旋转角度α为固定参数,即为不参与训练的参数。或者,所述第三旋转门的旋转角度φ0、所述第四旋转门的旋转角度θ0以及所述第五旋转门的旋转角度α均为所述目标可调参数。
需要说明的是,所述第四旋转门和第五旋转门的相关内容,可参照上述描述,此处不再赘述。也就是说,相较于如图3(c)所示的结构,该示例中,如图3(f)所示,L个训练层之后包括第四旋转门和第五旋转门,不包括第三旋转门。
这样,本公开方案有效提升了预设参数化量子电路的表达能力,同时,所使用的量子门的种类和数量少,待训练的目标可调参数的数量少,如此,为高效模拟哈密顿量奠定了基础,同时也为提升结果的精确度奠定了基础。
在本公开方案的一具体示例中,所述目标量子电路包含有M层,所述M为大于等于1小于等于L/2的正整数;
所述M层中的至少一层是基于如下所得:
将两个训练层中第一训练层的第一目标旋转门替换第一受控酉门,将所述两个训练层中第二训练层的第二目标旋转门替换第二受控酉门;其中,所述两个训练层为所述L个训练层中的任意两个训练层。
可以理解的是,该示例适用于上述第一种方式和第二种方式。
这里,由于所述目标量子电路是在预设参数化量子电路的基础上扩展所得,而且,是将预设参数化量子电路中不同层的两个目标旋转门分别替换为第一受控酉门和第二受控酉门所得,所以,所述目标量子电路至多包含有L/2层。
进一步地,在所述预设参数化量子电路中的各训练层均包含有目标旋转门的情况下,比如,各训练层均包含第一种方式的旋转门,即如图3(a)所示的旋转门,或者,各训练层均包含有第二种方式的旋转门,即如图3(d)所示的旋转门,此时,所述目标量子电路均包含有L/2层。
在一具体示例中,所述L个训练层中的至少两个训练层(比如第i个训练层和第i+1(或i+2等,此处仅为示例性说明,还可以为其他层)个训练层)包括:所述目标旋转门,所述第一旋转门,所述第二旋转门,此时,所述目标量子电路中存在一层,比如第层(为向上取整符号),是将第i+1个训练层(可对应以上所述的第一训练层)的目标旋转门(也即第一目标旋转门)替换为所述第一受控酉门,以及将第i个训练层(第二训练层)的目标旋转门(也即第二目标旋转门)替换为所述第二受控酉门后所得。
进一步地,由于所述M层中的至少一层是基于所述预设参数化量子电路中两个训练层所得,所以,在一示例中,所述M层中的至少一层包括:
两个第一旋转门;
两个第二旋转门;
第一受控酉门;
第二受控酉门。
进一步地,在另一示例中,按量子门的作用顺序,所述M层中的至少一层依次包括:
第一旋转门;
第二旋转门;
第一受控酉门;
第一旋转门;
第二旋转门;
第二受控酉门。
或者,在另一示例中,所述M层中的至少一层包括:
两个第二旋转门;
第一受控酉门;
第二受控酉门。
进一步地,在另一示例中,按量子门的作用顺序,所述M层中的至少一层依次包括:
第二旋转门;
第一受控酉门;
第二旋转门;
第二受控酉门。
这里,该示例中量子门的相关介绍可按照以上说明,此处不再赘述。
这样,本公开方案在预设参数化量子电路的基础上构造目标量子电路的过程中,有效提升了目标量子电路的表达能力,同时,所使用的量子门的种类和数量均少,而且,待训练的目标可调参数的数量也少,如此,为高效模拟哈密顿量奠定了基础,同时也为提升结果的精确度奠定了基础。
而且,基于预设参数化量子电路来构造目标量子电路的过程中,可以采用不同的构造方式,因此,本公开方案扩展性强。
在本公开方案的一具体示例中,所述两个训练层为L个训练层中的任意相邻的两个训练层。也就是说,所述M层中的至少一层是基于预设参数化量子电路中相邻的两个训练层所得。
在一具体示例中,所述L个训练层的任意相邻的两个训练层(比如第i个训练层和第i+1个训练层)中的每一训练层均包括:所述目标旋转门,所述第一旋转门,所述第二旋转门,此时,所述目标量子电路中存在一层,比如第层,是将第i+1个训练层(也即第一训练层)的目标旋转门(也即第一目标旋转门)替换为所述第一受控酉门,以及将第i个训练层(第二训练层)的目标旋转门(也即第二目标旋转门)替换为所述第二受控酉门后所得。
进一步地,在一示例中,所述目标量子电路中的各层均是基于预设参数化量子电路中相邻的两个训练层所得,比如,各层均是将所述预设参数化量子电路的相邻的两个训练层中的第一训练层的第一目标旋转门替换为第一受控酉门,以及将所述两个训练层中的第二训练层的第二目标旋转门替换为第二受控酉门后所得。此时,所述目标量子电路中第一受控酉门的数量和第二受控酉门的数量,均为所述预设参数化量子电路中的目标旋转门的数量的一半。
具体来说,在所述预设参数化量子电路中各训练层均包含有所述目标旋转门,第一旋转门和第二旋转门,以及各旋转门的作用顺序如图3(a)所示的情况下,所述目标量子电路的L/2层中的第层是基于如下所得:
将第i+1个训练层中的目标旋转门(也即第一目标旋转门)替换为所述第一受控酉门,以及将第i个训练层中的目标旋转门(也即第二目标旋转门)替换为所述第二受控酉门后所得。
具体地,如图4(a)所示,所述目标量子电路中的第层(i取值为1至L),按量子门的作用顺序包括:
旋转角度φi+1为z轴所对应的角度的第一旋转门RZ(φi+1);
旋转角度θi+1为y轴所对应的角度的第二旋转门RY(θi+1);
第一受控酉门;
旋转角度φi为z轴所对应的角度的第一旋转门RZ(φi);
旋转角度θi为y轴所对应的角度的第二旋转门RY(θi);
第二受控酉门。
或者,在所述预设参数化量子电路中各训练层均包含有所述目标旋转门和第二旋转门,以及各旋转门的作用顺序如图3(d)所示的情况下,所述目标量子电路的L/2层中的第层是基于如下所得:
将第i+1个训练层中的目标旋转门(也即第一目标旋转门)替换为所述第一受控酉门,以及将第i个训练层中的目标旋转门(也即第二目标旋转门)替换为所述第二受控酉门后所得。
具体地,如图4(b)所示,所述目标量子电路中的第层(i取值为1至L),按量子门的作用顺序包括:
旋转角度θi+1为y轴所对应的角度的第二旋转门RY(θi+1);
第一受控酉门;
旋转角度θi为y轴所对应的角度的第二旋转门RY(θi);
第二受控酉门。
需要说明的是,所述目标量子电路中不同层所作用的辅助寄存器、块编码寄存器以及主寄存器均相同。也就是说,实际应用中,可以先将所述预设参数化量子电路中的量子比特作为辅助寄存器,同时扩展出块编码寄存器和主寄存器后,再将预设参数化量子电路中各训练层中的目标旋转门替换为目标受控酉门,进而使得各层共用同一辅助寄存器、块编码寄存器和主寄存器。
可以理解的是,基于图3(a)和图3(c)结合得到的预设参数化量子电路,即可扩展得到如图4(c)所示的目标量子电路。
这样,本公开方案基于预设参数化量子电路来构造目标量子电路,该过程低消耗,而且,通过辅助寄存器即可控制目标受控酉门,并模拟得到目标演化时间t0所对应的目标量子态,相比于现有方案,本公开方案有效降低了所需的量子计算资源,增强了中等规模量子计算设备求解量子特征的可行性。
需要说明的是,本公开方案,如图4(a)或图4(b)所示,当辅助寄存器的量子态为|0>的情况下,激活所述目标量子电路中带空心的受控酉门也即第二受控酉门。当辅助寄存器的量子态为|1>的情况下,激活带实心的受控酉门U,也即第一受控酉门。也就是说,实际应用中,在辅助寄存器的当前的量子态确定的情况下,第一受控酉门工作,或第二受控酉门工作,而非两者共同工作。这样,本公开方案通过辅助寄存器即可控制目标受控酉门,进而模拟得到目标演化时间t0所对应的目标量子态,相比于现有方案,本公开方案有效降低了所需的量子计算资源,增强了中等规模量子计算设备求解量子特征的可行性。而且,本公开方案适用于任意哈密顿量,具有丰富的应用场景。
在本公开方案的一具体示例中,所述目标块编码矩阵可以具体为以下两种形式:
形式一:在所述目标块编码矩阵U能够通过泡利串的组合(比如线性组合)表示的情况下,所述目标量子电路中所述第一受控酉门为所述泡利串的组合(比如线性组合)所对应的等效电路,以及所述目标量子电路中所述第二受控酉门为所述泡利串的组合(比如线性组合)的共轭转置所对应的等效电路。
也就是说,该示例中,所述目标块编码矩阵U可以通过泡利串(比如泡利矩阵和单位矩阵的张量积)的线性组合来表示,此时,由于所述泡利串可以使用简单的预设量子电路来表示,所以,基于实现泡利串的预设量子电路,并通过酉算子的线性组合技术,构造得到第一受控酉门和第二受控酉门的等效电路。
需要说明的是,该形式一中,所述目标量子电路的辅助寄存器的第一输入态为预设初态,所述目标量子电路的块编码寄存器的第二输入态为预设初态,所述目标量子电路的主寄存器的第三输入态为所述第一量子态。
形式二:在所述目标块编码矩阵U不能够通过泡利串的线性组合表示的情况下,所述目标量子电路中所述第一受控酉门为用于实现目标块编码矩阵U的第一等效电路,以及所述目标量子电路中所述第二受控酉门为用于实现所述目标块编码矩阵U的共轭转置的第二等效电路;
其中,所述第一等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第三受控酉门;
受控于所述辅助寄存器和所述第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第四受控酉门;
这里,所述第一组块编码量子比特和第二组块编码量子比特组成所述块编码寄存器;所述第三受控酉门为第一块编码矩阵对应的酉门;所述第四受控酉门为所述第一块编码矩阵的共轭转置对应的酉门。所述第一块编码矩阵为所述目标哈密顿量的块编码矩阵,而且,所述第一块编码矩阵所包含的量子比特的数量小于所述m+n。举例来说,所述第一块编码矩阵所包含的量子比特的数量=m+n-1,此时,所述第一组块编码量子比特包含有一个块编码量子比特;所述第二组块编码量子比特包含m-1个块编码量子比特。
进一步地,所述第二等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的所述第四受控酉门;
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的所述第三受控酉门。
也就是说,该示例中,给出更为一般的量子电路结构,如此,来实现目标块编码矩阵U;具体地,如图4(d)所示,对于给定的任意的第一块编码矩阵UA,所述第一块编码矩阵UA所对应的量子比特的数量为m+n-1,此时,所述目标量子电路中第层的所述目标块编码矩阵U的等效电路,按作用顺序,依次包括:
作用于第一组块编码量子比特上的H门(Hadamard门);
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组量子比特和所述主寄存器的第三受控酉门;
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组量子比特和所述主寄存器的第四受控酉门;
受控于所述辅助寄存器,且作用于第一组量子比特上的泡利X门;
作用于第一组量子比特上的H门;
受控于所述辅助寄存器,且作用于所述块编码寄存器上的反射算子R。
这里,由于所述第三受控酉门为第一块编码矩阵UA的等效电路,所以,为便于描述,所述第三受控酉门可通过字符UA表示;同理,由于所述第四受控酉门为第一块编码矩阵UA的共轭转置的等效电路,所以,为便于描述,所述第四受控酉门可通过字符表示。
这里,所述反射算子R(Reflector)为:
R=2|0(m-1)+1><0(m-1)+1|-I
所述I为单元矩阵。
进一步地,如图4(d)所示,所述目标量子电路中第层的所述目标块编码矩阵U的共轭转置的等效电路,按作用顺序,依次包括:
受控于所述辅助寄存器,且作用于所述块编码寄存器上的反射算子R;
作用于第一组量子比特上的H门;
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组量子比特和所述主寄存器的第四受控酉门;
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组量子比特和所述主寄存器的第三受控酉门;
受控于所述辅助寄存器,且作用于第一组量子比特上的泡利X门;
作用于第一组量子比特上的H门。
可以理解的是,类似于图4(b),该示例中,如图4(e)所示,还可以删除图4(d)中的所有第一旋转门RZ(φi);进一步地,在目标量子电路中包含第三旋转门RZ(φ0)的情况下,还可删除该第三旋转门RZ(φ0),得到基于图3(d)和图3(e)扩展得到的目标量子电路,或者得到基于图3(d)和图3(f)扩展所得到目标量子电路,如此,以模拟一个偶函数,在实现同样效果的同时,可进一步降低一半的电路深度。
需要说明的是,本公开方案,如图4(d)或图4(e)所示,当辅助寄存器的量子态为|0>的情况下,激活所述目标量子电路的第二受控酉门中的带空心的反射算子R、带空心的带空心的UA以及带有空心的泡利X门。当辅助寄存器的量子态为|1>的情况下,激活所述目标量子电路的第一受控酉门中的带实心的带实心的UA、带实心的带有实心的泡利X门,以及带有实心的反射算子R。同理,当块编码寄存器的量子态为|0>的情况下,激活所述目标量子电路的第一受控酉门中的带空心的带空心的UA以及第二受控酉门中的带空心的带空心的UA。当块编码寄存器的量子态为|1>的情况下,激活所述目标量子电路的第一受控酉门中的带实心的带实心的以及所述第二受控酉门中的实心的
需要说明的是,该形式一或形式二中,所述目标函数f(x)=e-it cos x(所述i为虚数)。
基于此,本公开方案具有如下优势:
第一、本公开方案所需目标量子电路的宽度更小。与现有方案所需的辅助量子比特数量相比,本公开方案的目标量子电路中辅助量子比特的数量可以为一个,因此,相比于现有方案,本公开方案所使用的目标量子电路的宽度最小,进而为有效降低计算量、提升处理效率奠定了基础,同时,精度高。
第二、本公开方案更易实现。在量子门的复杂程度和数量上,与现有方案相比,本公开方案目标量子电路中所使用的量子门的数量和种类更少,比如,可以使用单量子比特控制的受控酉门,如第一受控酉门和第二受控酉门,因此,减少了电路深度,降低了所需量子计算资源,同时,增加了在中等规模量子计算设备中执行的可行性。
第三、实用性更强。本公开方案所构造的目标量子电路简洁、且代价更小,更具有实用性。
这样,本公开公开提供了目标块编码矩阵的具体实现形式,极大提高了在中等规模含噪量子设备上的实用性,并具有很强的可扩展性。
以下结合具体示例对本公开方案做进一步详细说明;具体来说,一个目标量子系统的时间演化可以由一个目标哈密顿量H(比如,对于包含有n个量子比特的目标量子系统而言,其目标哈密尔顿量H可表示为作用于n个量子比特上的2n×2n厄米特矩阵)和这个目标量子系统的初始量子态|ψ>决定。具体地说,对于一个由某个时间恒定的目标哈密顿量H决定的目标量子系统,在演化时间t的量子态可以表示为:
|ψt>=e-iHt|ψ>;
这里,是虚数,e-iHt称为演化算子。
这里,量子模拟是指模拟目标量子系统的演化算子e-iHt,比如使用量子设备或参数化量子电路来模拟该目标量子系统的演化算子e-iHt,从而以一定精度在该量子设备或参数化量子电路上近似制备出该目标量子系统在演化时间t下的量子态|ψt>。
进一步地,块编码(Block encoding)是一种量子系统的矩阵编码方式。以包含有n个量子比特的目标量子系统为例,其目标哈密尔顿量H可通过2n×2n厄米特矩阵(为便于描述,可使用字母H表示)表示,此时,存在2n+m×2n+m的矩阵U,且满足如下关系:
此时,矩阵U则称为矩阵H的块编码。这里,矩阵U即为以上所述的目标块编码矩阵;这里,目标块编码矩阵U所对应的量子比特的数量为m+n。
进一步地,若目标块编码矩阵U可以在量子设备上有效制备,那么,若第n个到第m个量子比特上为|0>的情况下,即可获得厄米特矩阵U的相关信息,如此,得到演化算子的相关信息,实现量子演化。
基于此,哈密顿量模拟的任务则可以描述为:
输入:一个目标哈密顿量H的目标块编码矩阵U,以及该目标量子系统的初始量子态(也即第一量子态)|ψ>和目标演化时间t0。
输出:目标量子系统在t时刻的量子态
具体地,本公开方案的目的是给出实用且高效的量子哈密顿量模拟方案。该方案主要分为两个部分,第一部分,基于量子信号处理或者量子神经网络来模拟一个目标函数,比如目标函数f(x)=e-it cos x(t为演化时间,x∈[-π,π])该部分可构建预设参数化量子电路,并训练预设参数化量子电路,使得所述预设参数化量子电路能够模拟该目标函数f(x)。第二部分是使用第一部分中得到的目标参数值来构造用于模拟哈密顿量时间演化的量子电路。
第一部分,即程序一,主要用于计算或优化辅助寄存器上旋转门的目标可调参数,该程序一为会被程序二(也即主程序)调用的子程序。
步骤11:输入目标演化时间误差容忍值∈(也即以上所述的第一误差条件)。
这里,所述误差容忍值∈能够约束用于模拟目标函数f(x)的预设参数化量子电路所输出的实际输出结果与目标输出结果之间的差异程度,如此,来约束模拟得到的目标量子态的精确度。
步骤12:构造待训练的预设参数化量子电路,并根据误差容忍值∈来确定待训练的预设参数化量子电路的训练层数,比如,包含L个训练层;进一步地,还可以基于误差容忍值∈来确定训练数据组的数量N。这里,所述L为大于等于2的偶数;所述N也为大于等于1的正整数。
这里,该示例中,所述预设参数化量子电路为包含有一个量子比特(该示例中可称为辅助量子比特,或称为辅助寄存器)的参数化电路。
需要说明的是,实际应用中,还可以构造包含有两个或两个以上的量子比特的预设参数化量子电路来模拟目标函数f(x),本公开方案对此不作限制,只要能够模拟目标函数,并扩展得到能够求解特征相位的目标参数化量子电路,该预设参数化量子电路就在本公开方案的保护范围内。
该示例中,所述预设参数化量子电路的L个训练层中的每一训练层均包含有量子旋转门序列,且每一训练层中的量子旋转门序列均相同。
可以理解的是,实际应用中,L个训练层中不同训练层所包含的量子旋转门序列可以相同,也可以不相同,或者,部分训练层所包含的量子旋转门序列相同,另外部分训练层所包含的量子旋转门序列不相同等,本公开方案对此不作具体限制。
进一步地,该示例中,以所述L个训练层中第i个训练层所包含的量子旋转门序列为例进行说明。如图3(a)所示,基于所述量子旋转门序列中旋转门的作用顺序可知,第i个训练层所包含的量子旋转门序列依次包括:
旋转角度φi为z轴所对应的角度的第一旋转门RZ(φi);
旋转角度θi为y轴所对应的角度的第二旋转门RY(θi);
旋转参数xj为z轴所对应的角度的目标旋转门RZ(xj)。
这里,所述第一旋转门RZ(φi)的旋转角度φi和第二旋转门RY(θi)的旋转角度θi为第i个训练层中的目标可调参数,所述i为大于等于1小于等于L的整数。
进一步地,该示例中,所述预设参数化量子电路中L个训练层之后,还包括其他旋转门。
具体地,在一示例中,如图3(b)所示,所述预设参数化量子电路在L个训练层之后,还包括:
旋转角度φ0为z轴所对应的角度的第三旋转门RZ(φ0);
旋转角度θ0为y轴所对应的角度的第四旋转门RY(θ0)。
基于此,如图3(b)所示的预设参数化量子电路的数学表达式可具体为:
或者,在另一示例中,如图3(c)所示,所述预设参数化量子电路在L个训练层之后,还包括:
旋转角度φ0为z轴所对应的角度的第三旋转门RZ(φ0);
旋转角度φ0为y轴所对应的角度的第四旋转门RY(θ0);
以及旋转角度α为z轴所对应的角度的第五旋转门RZ(α)。
这里,所述旋转角度φ0、旋转角度θ0和旋转角度α均为目标可调参数。
基于此,如图3(c)所示的预设参数化量子电路的数学表达式可具体为:
或者,所述旋转角度φ0和旋转角度θ0均为目标可调参数,而旋转角度α为固定参数,不参与训练。
基于此,如图3(c)所示的预设参数化量子电路的数学表达式可具体为:
需要说明的是,L个训练层中各训练层的电路结构可以参照图3(a)所示的结构,该图3(b)和图3(c)中未示出。
需要说明的是,由于该预设参数化量子电路包括一个量子比特,所以,可以使用经典计算设备来有效地精确模拟该预设参数化量子电路的运行,即不需要消耗量子计算资源,如此,节约了量子计算资源,同时,也降低了处理成本。
进一步地,可以理解的是,实际应用中,当预设参数化量子电路所包含的量子比特的数量较少(比如,20-30个量子比特)的情况下,均可通过模拟电路的方式,在经典计算设备中计算得到目标可调参数的目标参数值,如此,在计算效率允许范围内,最大程度避免消耗量子计算资源。
步骤13:准备训练数据组;比如,准备N个训练数据点用于训练以上所述的预设参数化量子电路。
该示例以图3(c)所示的预设参数化量子电路为例进行说明,且所述旋转角度α为目标可调参数,参与后续训练流程。相应地,基于图3(c)所示的预设参数化量子电路扩展得到的目标量子电路,如图4(c)所示。
步骤14:随机生成L+1个的参数值θ,以及L+1个的参数值φ,以及1个参数值α。
这里,L+1个的参数值θ,可分别记为θ0和(i为大于等于1小于等于L的正整数)。为便于记录,还可使用向量来表示,即θ={θ0,θ1,…,θi,…,θL}。
同理,L+1个的参数值可分别记为φ0和(i为大于等于1小于等于L的正整数)。为便于记录,还可使用向量φ来表示,即φ={φ0,φ1,…,φi,…φL}。
此时,所述预设参数化量子电路可表示为Ux(α,θ,φ)。
步骤15:对于每个旋转参数xj,1≤j≤N而言,均执行如下操作:
(a)使用经典模拟器(也即在经典计算设备上)模拟以上所述的预设参数化量子电路Ux(α,θ,φ);而且,对于每个xj,可具体得到预设参数化量子电路
(b)输入预设初态,比如|0>,使用经典模拟器模拟获取所述预设参数化量子电路的实际输出态与预设初态(比如|0>)的内积<0|Ux(α,θ,φ)|0>,也即得到实际输出结果,并记为yj。
对于每个xj均执行上述操作后,即完成上述操作后,获得一组实际输出结果共N个。
步骤16:定义函数f(x)=e-itcos(x);将实际输出结果与目标输出结果之间的2-范数作为损失函数,即损失函数L(α,θ,φ)为:
这里,可以理解的是,实际应用中,损失函数还可以是其他任何一种刻画距离的度量函数,例如常用的平均绝对误差函数、均方误差函数和交叉熵函数等。可以根据数据大小、硬件环境、学习精度或者收敛速度等因素选择合适的损失函数,本公开方案对此不作具体限制。
步骤17:基于损失函数L(α,θ,φ)计算得到的损失值,并进行优化,比如通过梯度下降法,来对目标可调参数α,θ和φ进行调整,以最小化L(α,θ,φ);
其中,目标可调参数θ包括θ0和也即θ={θ0,θ1,…,θi,…,θL},目标可调参数φ包括φ0和也即φ={φ0,φ1,…,φi,…φL}。
这里,实际应用中,在经典计算设备上,可使用常见的梯度下降法,也可以使用其他更加科学有效的优化方法,对目标可调参数α,θ0,以及目标可调参数φ0和进行优化,如此,来最小化损失函数的损失值,本公开方案对具体优化方式不作限制。
步骤18:在目标可调参数调整后,重复步骤15-17,直至损失函数L(α,θ,φ)收敛,或达到迭代次数,得到各目标可调参数的最优参数值(也即目标参数值),分别为和
这里,
可以理解的是,重复上述优化过程,使损失函数的损失值达到最小或者达到收敛状态,或者达到迭代次数,此时,则可认为实际输出结果yj趋近于目标输出结果目标可调参数的当前参数值和即为最优参数值。
步骤19:输出最优参数值(也即目标参数值),和共2L+3个。
可以理解的是,实际应用中,在不考虑计算代价的情况下,该程序一可在经典计算设备中运行,也可以在量子计算设备中运行,本公开方案对此不作具体限制。
实际应用中,上述程序一的实现,并非是唯一的,比如,在进行目标可调参数初始化的过程中(比如上述步骤14),可以利用这些目标可调参数的内在性质,或者设置其初始值,提高优化效率;或者,还可以使用函数分析的方法来直接获取目标可调参数的最优参数值。换言之,实际应用中,可以基于具体应用场景和硬件环境等因素来选择合适的实现方式。
举例来说,使用函数分析方法来计算目标可调角度,具体包括:
输入目标函数f(x),可简写为f。计算得到能在目标定义域内近似目标函数f的目标傅里叶级数F(x)。以及,计算得到其他傅里叶级数P(x)和Q(x);其中,
根据下述等式递归计算目标可调参数α,θ和φ的最优参数值:
这里,所述Q*(x)为Q(x)的复数共轭,P*(x)为P(x)的复数共轭。最后,输出最优参数值和
实际应用中,还可以使用任何能以一定精度近似目标函数的三角多项式,来优化得到目标可调参数的最优参数值。
第二部分,即程序二,该程序为主程序,主要用于哈密顿量模拟。
可以理解的是,实际应用中,在不考虑计算代价的情况下,该程序二也可在经典计算设备中运行,也可以在量子计算设备中运行,本公开方案对此不作具体限制。
具体地,如图5所示,该主程序的具体步骤包括:
步骤21:对预设参数化量子电路进行扩展,扩展成具有n+m+1个量子比特的目标量子电路,以使得所述目标量子电路能够输出得到演化后的目标量子态。该示例,以图4(c)所示目标量子电路为例,该预设参数化量子电路中的量子比特为辅助量子比特,可称为辅助寄存器;相应地,新增加的或扩展出的m个量子比特为块编码量子比特,可统称为块编码寄存器,新增加的或扩展出的另外n个量子比特为主量子比特,该n个主量子比特可统称为主寄存器。
这里,所述块编码寄存器所占用的量子位(比如,占用连续的m个量子位),位于所述辅助寄存器(比如,占用第一个量子位)和主寄存器(比如,占用连续的最后n个量子位)之间。
也就是说,所述目标量子电路包括辅助寄存器、块编码寄存器和主寄存器;其中,所述辅助寄存器包括一个辅助量子比特;所述主寄存器包括n个主量子比特。这里,所述n是基于第一量子态所对应的量子比特的数量(也即目标量子系统所包含的量子比特的数量)而确定出的,比如,n为所述目标量子系统中所包含的量子比特的数量。换言之,所述主寄存器所包含的主量子比特的数量,与所述目标量子系统所包含的量子比特的数量相同。
进一步地,所述块编码寄存器包括m个量子比特;所述m是基于所述目标哈密顿量所对应的目标块编码矩阵所对应的量子比特的数量,以及所述目标量子系统中所包含的量子比特的数量而确定的,比如,所述块编码寄存器所对应的量子比特的数量,等于所述目标块编码矩阵所对应的量子比特的数量与目标量子系统中所包含的量子比特的数量的之差。
进一步地,所述目标量子电路基于如下所得:将所述预设参数化量子电路中的量子比特作为辅助寄存器,并扩展出包含有m个块编码量子比特的块编码寄存器和包含有n个主量子比特的主寄存器,同时,将所述预设参数化量子电路中作用于所述辅助寄存器的第一目标旋转门替换为所述第一受控酉门,以及将所述预设参数化量子电路中作用于所述辅助寄存器的第二目标旋转门替换为所述第二受控酉门。
进一步地,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;也就是说,将所述预设参数化量子电路中不同训练层的目标旋转门替换成不同的受控酉门,比如,将所述预设参数化量子电路中一训练层中的目标旋转门(为便于描述,可称为第一目标旋转门)替换成第一受控酉门,同时,将所述预设参数化量子电路中另一训练层中的目标旋转门(为便于描述,可称为第二目标旋转门)替换成第二受控酉门,如此,得到目标量子电路。
可以理解的是,由于所述目标量子电路是在预设参数化量子电路的基础上扩展所得,而且,是将预设参数化量子电路中不同层的两个目标旋转门分别替换为第一受控酉门和第二受控酉门所得,所以,在预设参数化量子电路包含有L层的情况下,所述目标量子电路至多包含有L/2层。
具体地,从所述预设参数化量子电路中扩展出包含有n个主量子比特的主寄存器,以及扩展出包含有m个块编码量子比特的块编码寄存器,同时,将所述预设参数化量子电路的相邻两个训练层中的目标旋转门分别替换为第一受控酉门和第二受控酉门,比如,将第i+1个训练层的目标旋转门替换为所述第一受控酉门,以及将第i个训练层的目标旋转门替换为所述第二受控酉门,得到如图4(a)所示的目标量子电路中的第层的结构示意图,该第层,按各量子门的作用顺序来说,具体包括:
旋转角度φi+1为z轴所对应的角度的第一旋转门RZ(φi+1);
旋转角度θi+1为y轴所对应的角度的第二旋转门RY(θi+1);
第一受控酉门;
旋转角度φi为z轴所对应的角度的第一旋转门RZ(φi);
旋转角度θi为y轴所对应的角度的第二旋转门RY(θi);
第二受控酉门。
这里,为便于描述,可将目标量子电路中作用于辅助量子比特的相关参数化量子电路称为所述目标量子电路的子电路。可以理解的是,该子电路也包含有L/2层。进一步地,如图4(a)所示,该子电路中的每一层中均包含目标可调参数;以该子电路中的第层为例,包含有:
旋转角度φi+1为z轴所对应的角度的第一旋转门RZ(φi+1);
旋转角度θi+1为y轴所对应的角度的第二旋转门RY(θi+1);
旋转角度φi为z轴所对应的角度的第一旋转门RZ(φi);
旋转角度θi为y轴所对应的角度的第二旋转门RY(θi);
这里,所述旋转角度φi+1、旋转角度θi+1以及旋转角度φi和旋转角度θi为当前层的目标可调参数。
可以理解的是,由于所述目标量子电路是在所述预设参数化量子电路的基础上扩展所得,所以,与预设参数化量子电路类似,所述目标量子电路在L/2层之后,还包括其他旋转门。
具体地,在一示例中,所述目标量子电路中L/2层之后,还包括如图3(b)所示的第三旋转门RZ(φ0)和第四旋转门RY(θ0)。这里,所述旋转角度φ0和旋转角度θ0均为目标可调参数。
或者,在另一示例中,所述目标量子电路中L/2层之后,还包括如图3(c)所示的第三旋转门RZ(φ0)和第四旋转门RY(θ0),以及第五旋转门RZ(α)。这里,所述旋转角度φ0和旋转角度θ0均为目标可调参数;而所述旋转角度α为固定值。或者,所述旋转角度φ0、旋转角度θ0以及旋转角度α均为目标可调参数。具体内容可参见以上陈述,此处不再赘述。
步骤22:输入目标哈密顿量H的目标块编码矩阵U,误差容忍值∈>0,目标演化时间t0,将目标量子电路中辅助寄存器的第一输入态设置为预设初态,比如|0>或|1>;将目标量子电路中块编码寄存器的第二输入态设置为预设初态,比如|0>或|1>;将目标量子电路中主寄存器的第三输入态设置为待演化量子态|ψ>(也即以上所述的初始量子态,或第一量子态)。
这里,可以理解的是,所述目标量子电路中的第一受控酉门为目标哈密顿量H的目标块编码矩阵U的等效电路,如图4(a)至图4(c)所示,为便于描述,该第一受控酉门也可使用字符U表示。进一步地,所述目标量子电路中的第二受控酉门为所述目标块编码矩阵U的共轭转置的等效电路,为便于描述,该第二受控酉门可使用字符表示。
进一步地,在一具体示例中,当辅助寄存器的量子态为|0>的情况下,激活所述目标量子电路的第二受控酉门中的带空心的反射算子R、带空心的带空心的UA以及带有空心的泡利X门。当辅助寄存器的量子态为|1>的情况下,激活所述目标量子电路的第一受控酉门中的带实心的带实心的UA、带实心的带有实心的泡利X门,以及带有实心的反射算子R。同理,当块编码寄存器的量子态为|0>的情况下,激活所述目标量子电路的第一受控酉门中的带空心的带空心的UA以及第二受控酉门中的带空心的带空心的UA。当块编码寄存器的量子态为|1>的情况下,激活所述目标量子电路的第一受控酉门中的带实心的带实心的以及所述第二受控酉门中的实心的
步骤23:将目标演化时间t0和误差容忍值∈输入到“程序一”,并运行“程序一”,获取输出的最优参数值(也即目标参数值):和
这里,即
步骤24:如图4(c)所示,输入最优参数值和以及将目标哈密顿量H的目标块编码矩阵U作用在n+m+1个量子比特上的目标量子电路,也即将与目标块编码矩阵U等价的第一受控酉门,以及与目标块编码矩阵的共轭转置等价的第二受控酉门作用到n+m+1个量子比特上的目标量子电路。
步骤25:输出经过目标量子电路演化后,主寄存器上的输出量子态即为目标量子态。
具体地,在所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息均为|0>的情况下,可认为满足所述预设条件。
可以理解的是,实际应用中,当所述第一状态信息和第二状态信息为|0>的概率大于一阈值的情况下,即可认为两者均为|0>,此时,认为满足预设条件;相应地,所述主寄存器的输出态即为目标量子态。
需要说明的是,本公开方案所模拟的目标函数f(x)=e-itcos(x)。
扩展方案
扩展方式一:在上述具体示例中,还定义一个偶函数作为“程序一”中的目标函数。举例来说,可以删除“程序一”和“程序二”中如图3(a)所示的预设参数化量子电路中所有第一旋转门RZ(φi),以及如图3(c)中第三旋转门RZ(φ0),得到如图3(d)和图3(e),或者图3(d)和图3(f)所示的结构,以来模拟一个偶函数,在实现同样效果的同时,可进一步降低一半的电路深度。
扩展方式二:使用上述形式二所示的等效电路(如图4(d)或图4(e)所示)来作为目标块编码矩阵U的等效电路。这里,具体描述可参见上述方式二,此处不再赘述。
案例展示
以下通过具体案例来展示本公开方案。
该案例中,随机选择10个4量子比特的随机哈密顿量(也即目标哈密顿量)作为所要模拟的10个目标量子系统。此外,设置目标演化时间t0=9,并选取随机哈密顿量Hj所对应的目标块编码矩阵U(j),即:
使用图4(c)所示的目标量子电路,通过程序一得到最优参数值,即即α,θ0,φ0,和此时,试验的目标是对于每一个随机哈密顿量Hj,使用目标量子电路模拟实现哈密顿量演化过程。
这里,基于本公开方案数值模拟得到的演化后的目标量子态与理论演化后的量子态进行对比,其平均误差为1.12244×10-13,由此可验证本公开方案的正确性。
综上可知,本公开方案能够适应近期量子计算机,而且,具有如下特征:
第一,本公开方案可以只使用一个辅助量子比特即可得到演化后的目标量子态,降低了所需的量子资源,增强了中等规模量子计算设备求解量子特征的可行性。
第二,本公开方案适用于任意哈密顿量,比如,适用于任何可以有效制备目标哈密顿量所对应的目标块编码矩阵U的应用场景,具有丰富的应用场景。
第三,本公开方案还兼具实用性、高效性、确定性、扩展性以及创新性;具体地,实用性指本公开方案所需电路的宽度低,可以在近期的量子计算机上实现,具有丰富的落地场景;高效性指本公开方案可以低消耗地构造量子电路;确定性指本公开方案可以以极高的概率得到满足精度要求的估计值;扩展性指本公开方案可应用于大规模的哈密顿量演化;创新性指本公开方案提供了新颖且高效的量子电路实现量子哈密顿量模拟。
本公开方案还提供了一种哈密顿量模拟装置,如图6所示,包括:
参数调整单元601,用于确定目标量子电路的子电路中目标可调参数的目标参数值;其中,所述目标参数值为在目标演化时间t0下所对应的、且满足第一误差条件的参数值;所述目标量子电路包含有辅助寄存器、块编码寄存器和主寄存器,所述子电路作用于所述辅助寄存器;所述目标量子电路中还包含有受控于所述辅助寄存器且作用于所述块编码寄存器和所述主寄存器的目标受控酉门,所述目标受控酉门用于模拟目标哈密顿量所对应的目标块编码矩阵U;所述目标受控酉门包括与目标块编码矩阵U等价的第一受控酉门,和与所述目标块编码矩阵U的共轭转置等价的第二受控酉门;所述目标块编码矩阵U所对应的量子比特的数量为m+n,所述n为所述目标哈密顿量所对应的量子比特的数量,所述m为所述块编码寄存器所包含的量子比特的数量;
输出单元602,用于在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,得到所述目标量子电路中所述主寄存器的目标输出量子态,其中,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间t0下所对应的目标量子态。
在本公开方案的一具体示例中,所述输出单元602,还用于:
在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,获取所述目标量子电路中所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息;
在所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息均满足预设条件的情况下,得到所述主寄存器的目标输出量子态,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间t0下所对应的目标量子态。
在本公开方案的一具体示例中,所述参数调整单元601,具体用于:
将训练完成的预设参数化量子电路中所述目标可调参数的目标参数值,作为所述子电路中目标可调参数的目标参数值;所述训练完成的预设参数化量子电路用于模拟目标函数f(x);所述目标函数f(x)用于表征演化时间t与自变量x之间关联关系;
其中,所述目标量子电路是通过如下方式得到:将所述预设参数化量子电路中的量子比特作为辅助寄存器,并扩展出块编码寄存器和主寄存器,将所述预设参数化量子电路中作用于所述辅助寄存器的第一目标旋转门替换为所述第一受控酉门,以及将所述预设参数化量子电路中作用于所述辅助寄存器的第二目标旋转门替换为所述第二受控酉门;
其中,所述第一目标旋转门的第一旋转参数和第二目标旋转门的第二旋转参数均为所述目标函数f(x)的自变量x;所述子电路包含有所述预设参数化量子电路中除所述第一目标旋转门和第二目标旋转门之外的至少部分电路。
在本公开方案的一具体示例中,所述参数调整单元601,还用于:
在所述预设参数化量子电路的旋转参数x取值为N个数据点中的任一数据点xj的情况下,获取所述预设参数化量子电路的实际输出态;其中,所述实际输出态为在预设初态作为输入态的情况下,所述预设参数化量子电路在所述目标可调参数为当前参数值下的输出态;
基于所述实际输出态与所述预设初态得到实际输出结果yj;其中,所述实际输出结果yj为所述实际输出态与预设输入态之间的内积;所述N为大于等于1的正整数,所述j为大于等于1小于等于N的正整数;所述旋转参数x包括所述第一旋转参数和所述第二旋转参数;
得到N个实际输出结果yj;
在确定满足迭代终止条件的情况下,将所述目标可调参数的当前参数值作为训练完成的所述预设参数化量子电路中所述目标可调参数的目标参数值;
其中,所述迭代终止条件包括以下至少之一:
基于所述N个实际输出结果yj与N个目标输出结果确定预设损失函数的损失值满足收敛条件;所述目标输出结果
当前迭代次数到达预设次数。
在本公开方案的一具体示例中,所述参数调整单元601,还用于:
在确定不满足所述迭代终止条件的情况下,对所述目标可调参数的参数值进行调整;
重新在所述预设参数化量子电路的旋转参数x取值为N个数据点中的任一数据点xj的情况下,获取所述预设参数化量子电路的实际输出态,以得到实际输出结果yj;
重新得到N个实际输出结果yj,直至满足所述迭代终止条件为止。
在本公开方案的一具体示例中,所述预设参数化量子电路包括有L个训练层;所述L为大于等于2的偶数,所述L的取值与所述第一误差条件有关;
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数x用于对第一角度进行旋转操作;其中,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;
用于对第二角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第一旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第一旋转门的旋转角度φ和所述第二旋转门的旋转角度θ为所述目标可调参数;
或者,
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数x用于对第一角度进行旋转操作;其中,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第二旋转门的旋转角度θ为所述目标可调参数。
在本公开方案的一具体示例中,还满足以下至少之一:
所述第一角度为z轴所对应的角度;
所述第二角度为z轴所对应的角度;
所述第三角度为y轴所对应的角度。
在本公开方案的一具体示例中,在所述L个训练层的任一训练层中包含有所述目标旋转门、所述第一旋转门以及所述第二旋转门的情况下,各旋转门的作用顺序为:
所述第一旋转门、第二旋转门和目标旋转门;
或者,
在所述L个训练层的任一训练层中包含有所述目标旋转门和所述第二旋转门的情况下,各旋转门的作用顺序为:第二旋转门和目标旋转门。
在本公开方案的一具体示例中,所述预设参数化量子电路的L个训练层之后,还包括其他旋转门。
在本公开方案的一具体示例中,所述目标量子电路包含有M层,所述M为大于等于1小于等于L/2的正整数;
所述M层中的至少一层是基于如下所得:
将两个训练层中第一训练层的第一目标旋转门替换第一受控酉门,将所述两个训练层中第二训练层的第二目标旋转门替换第二受控酉门;其中,所述两个训练层为所述L个训练层中的任意两个训练层。
在本公开方案的一具体示例中,所述两个训练层为L个训练层中的任意相邻的两个训练层。
在本公开方案的一具体示例中,
在所述目标块编码矩阵U能够通过泡利串的组合表示的情况下,所述目标量子电路中所述第一受控酉门为所述泡利串的组合所对应的等效电路,以及所述目标量子电路中所述第二受控酉门为所述泡利串的组合的共轭转置所对应的等效电路;
或者,
在所述目标块编码矩阵U不能够通过泡利串的组合表示的情况下,所述目标量子电路中所述第一受控酉门为用于实现所述目标块编码矩阵U的第一等效电路,以及所述目标量子电路中所述第二受控酉门为用于实现所述目标块编码矩阵U的共轭转置的第二等效电路;
其中,所述第一等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第三受控酉门,第三受控酉门为第一块编码矩阵对应的酉门;
受控于所述辅助寄存器和所述第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第四受控酉门,第四受控酉门为所述第一块编码矩阵的共轭转置对应的酉门;
其中,所述第二等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和所述第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第四受控酉门;
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第三受控酉门;
其中,所述第一组块编码量子比特和第二组块编码量子比特组成所述块编码寄存器;所述第一块编码矩阵为所述目标哈密顿量的块编码矩阵,所述第一块编码矩阵所包含的量子比特的数量小于所述m+n。
本公开实施例的装置的各单元的具体功能和示例的描述,可以参见上述方法实施例中对应步骤的相关描述,在此不再赘述。
本公开方案还提供了一种存储有计算机指令的非瞬时计算机可读存储介质,当至少一个量子处理单元执行时,所述计算机指令使得所述至少一个量子处理单元执行以上应用量子计算设备的所述方法。
本公开方案还提供了一种计算机程序产品,包括计算机程序,所述计算机程序在被至少一个量子处理单元执行时实现应用于量子计算设备所述的方法。
本公开方案还提供了一种计算设备,所述计算设备包括:
至少一个量子处理单元;
存储器,耦合到所述至少一个QPU并用于存储可执行指令,
所述指令被所述至少一个量子处理单元执行,以使所述至少一个量子处理单元能够执行应用于量子计算设备所述的方法。
可以理解的是,本公开方案所述的中使用的量子处理单元(quantum processingunit,QPU),也可称为量子处理器或量子芯片,可以涉及包括多个以特定方式互连的量子比特的物理芯片。
而且,可以理解的是,本公开方案所述的量子比特可以指量子计算设备的基本信息单元。量子比特包含在QPU中,并推广了经典数字比特的概念。
进一步地,根据本公开的实施例,本公开还提供了一种计算设备、一种可读存储介质和一种计算机程序产品。
图7示出了可以用来实施本公开的实施例的示例计算设备700的示意性框图。计算设备旨在表示各种形式的数字计算机,诸如,膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。计算设备还可以表示各种形式的移动装置,诸如,个人数字助理、蜂窝电话、智能电话、可穿戴设备和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例,并且不意在限制本文中描述的和/或者要求的本公开的实现。
如图7所示,设备700包括计算单元701,其可以根据存储在只读存储器(ROM)702中的计算机程序或者从存储单元708加载到随机访问存储器(RAM)703中的计算机程序,来执行各种适当的动作和处理。在RAM 703中,还可存储设备700操作所需的各种程序和数据。计算单元701、ROM 702以及RAM 703通过总线704彼此相连。输入/输出(I/O)接口705也连接至总线704。
设备700中的多个部件连接至I/O接口705,包括:输入单元706,例如键盘、鼠标等;输出单元707,例如各种类型的显示器、扬声器等;存储单元708,例如磁盘、光盘等;以及通信单元709,例如网卡、调制解调器、无线通信收发机等。通信单元709允许设备700通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
计算单元701可以是各种具有处理和计算能力的通用和/或专用处理组件。计算单元701的一些示例包括但不限于中央处理单元(CPU)、图形处理单元(GPU)、各种专用的人工智能(AI)计算芯片、各种运行机器学习模型算法的计算单元、数字信号处理器(DSP)、以及任何适当的处理器、控制器、微控制器等。计算单元701执行上文所描述的各个方法和处理,例如哈密顿量模拟方法。例如,在一些实施例中,哈密顿量模拟方法可被实现为计算机软件程序,其被有形地包含于机器可读介质,例如存储单元708。在一些实施例中,计算机程序的部分或者全部可以经由ROM 702和/或通信单元709而被载入和/或安装到设备700上。当计算机程序加载到RAM 703并由计算单元701执行时,可以执行上文描述的哈密顿量模拟方法的一个或多个步骤。备选地,在其他实施例中,计算单元701可以通过其他任何适当的方式(例如,借助于固件)而被配置为执行哈密顿量模拟方法。
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、现场可编程门阵列(FPGA)、专用集成电路(ASIC)、专用标准产品(ASSP)、芯片上系统的系统(SOC)、复杂可编程逻辑设备(CPLD)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括:实施在一个或者多个计算机程序中,该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释,该可编程处理器可以是专用或者通用可编程处理器,可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令,并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
用于实施本公开的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。这些程序代码可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器,使得程序代码当由处理器或控制器执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在机器上执行、部分地在机器上执行,作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本公开的上下文中,机器可读介质可以是有形的介质,其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的程序。机器可读介质可以是机器可读信号介质或机器可读储存介质。机器可读介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备,或者上述内容的任何合适组合。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦除可编程只读存储器(EPROM或快闪存储器)、光纤、便捷式紧凑盘只读存储器(CD-ROM)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
为了提供与用户的交互,可以在计算机上实施此处描述的系统和技术,该计算机具有:用于向用户显示信息的显示装置(例如,CRT(阴极射线管)或者LCD(液晶显示器)监视器);以及键盘和指向装置(例如,鼠标或者轨迹球),用户可以通过该键盘和该指向装置来将输入提供给计算机。其它种类的装置还可以用于提供与用户的交互;例如,提供给用户的反馈可以是任何形式的传感反馈(例如,视觉反馈、听觉反馈、或者触觉反馈);并且可以用任何形式(包括声输入、语音输入、或者触觉输入)来接收来自用户的输入。
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如,作为数据服务器)、或者包括中间件部件的计算系统(例如,应用服务器)、或者包括前端部件的计算系统(例如,具有图形用户界面或者网络浏览器的用户计算机,用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如,通信网络)来将系统的部件相互连接。通信网络的示例包括:局域网(LAN)、广域网(WAN)和互联网。
计算机系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器,也可以为分布式系统的服务器,或者是结合了区块链的服务器。
应该理解,可以使用上面所示的各种形式的流程,重新排序、增加或删除步骤。例如,本公开中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行,只要能够实现本公开公开的技术方案所期望的结果,本文在此不进行限制。
上述具体实施方式,并不构成对本公开保护范围的限制。本领域技术人员应该明白的是,根据设计要求和其他因素,可以进行各种修改、组合、子组合和替代。任何在本公开的精神和原则之内所作的修改、等同替换和改进等,均应包含在本公开保护范围之内。
Claims (27)
1.一种哈密顿量模拟方法,包括:
目标量子电路的子电路中目标可调参数的目标参数值;其中,所述目标参数值为在目标演化时间下所对应的、且满足第一误差条件的参数值;所述第一误差条件为预设的误差容忍值,用于约束目标参数值的精度,以约束得到的目标输出量子态的精度;所述目标量子电路包含有一个辅助量子比特所组成的辅助寄存器、m个块编码量子比特所组成的块编码寄存器和n个主量子比特所组成的主寄存器;所述子电路作用于所述辅助寄存器;所述目标量子电路中还包含有受控于所述辅助寄存器且作用于所述块编码寄存器和所述主寄存器的目标受控酉门,所述目标受控酉门用于模拟目标哈密顿量所对应的目标块编码矩阵;所述目标受控酉门包括与目标块编码矩阵等价的第一受控酉门,和与所述目标块编码矩阵的共轭转置等价的第二受控酉门;所述目标块编码矩阵所对应的量子比特的数量为m+n,所述n为所述目标哈密顿量所对应的量子比特的数量;
在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,得到所述目标量子电路中所述主寄存器的目标输出量子态,其中,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间下所对应的目标量子态。
2.根据权利要求1所述的方法,还包括:
在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,获取所述目标量子电路中所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息;
其中,所述得到所述目标量子电路中所述主寄存器的目标输出量子态,包括:
在所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息均满足预设条件的情况下,得到所述主寄存器的目标输出量子态,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间下所对应的目标量子态。
3.根据权利要求1所述的方法,其中,确定目标量子电路的子电路中目标可调参数的目标参数值,包括:
将训练完成的预设参数化量子电路中所述目标可调参数的目标参数值,作为所述子电路中目标可调参数的目标参数值;所述训练完成的预设参数化量子电路用于模拟目标函数;所述目标函数用于表征演化时间与自变量之间关联关系;
其中,所述目标量子电路是通过如下方式得到:将所述预设参数化量子电路中的量子比特作为辅助寄存器,并扩展出块编码寄存器和主寄存器,将所述预设参数化量子电路中作用于所述辅助寄存器的第一目标旋转门替换为所述第一受控酉门,以及将所述预设参数化量子电路中作用于所述辅助寄存器的第二目标旋转门替换为所述第二受控酉门;
其中,所述第一目标旋转门的第一旋转参数和第二目标旋转门的第二旋转参数均为所述目标函数的自变量;所述子电路包含有所述预设参数化量子电路中除所述第一目标旋转门和第二目标旋转门之外的至少部分电路。
4.根据权利要求3所述的方法,还包括:
在所述预设参数化量子电路的旋转参数取值为N个数据点中的任一数据点的情况下,获取所述预设参数化量子电路的实际输出态;其中,所述实际输出态为在预设初态作为输入态的情况下,所述预设参数化量子电路在所述目标可调参数为当前参数值下的输出态;
基于所述实际输出态与所述预设初态得到实际输出结果;其中,所述实际输出结果为所述实际输出态与预设初态之间的内积;所述N为大于等于1的正整数,所述j为大于等于1小于等于N的正整数;所述旋转参数包括所述第一旋转参数和所述第二旋转参数;
得到N个实际输出结果;
在确定满足迭代终止条件的情况下,将所述目标可调参数的当前参数值作为训练完成的所述预设参数化量子电路中所述目标可调参数的目标参数值;
其中,所述迭代终止条件包括以下至少之一:
基于所述N个实际输出结果与N个目标输出结果,确定预设损失函数的损失值满足收敛条件;所述目标输出结果;
当前迭代次数到达预设次数。
5.根据权利要求4所述的方法,还包括:
在确定不满足所述迭代终止条件的情况下,对所述目标可调参数的参数值进行调整;
重新在所述预设参数化量子电路的旋转参数取值为N个数据点中的任一数据点的情况下,获取所述预设参数化量子电路的实际输出态,以得到实际输出结果;
重新得到N个实际输出结果,直至满足所述迭代终止条件为止。
6.根据权利要求3-5任一项所述的方法,其中,所述预设参数化量子电路包括有L个训练层;所述L为大于等于2的偶数,所述L的取值与所述第一误差条件有关;
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数用于对第一角度进行旋转操作;其中,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;
用于对第二角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第一旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第一旋转门的旋转角度和所述第二旋转门的旋转角度为所述目标可调参数;
或者,
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数用于对第一角度进行旋转操作;其中,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第二旋转门的旋转角度为所述目标可调参数。
7.根据权利要求6所述的方法,其中,还满足以下至少之一:
所述第一角度为z轴所对应的角度;
所述第二角度为z轴所对应的角度;
所述第三角度为y轴所对应的角度。
8.根据权利要求6所述的方法,其中,
在所述L个训练层的任一训练层中包含有所述目标旋转门、所述第一旋转门以及所述第二旋转门的情况下,各旋转门的作用顺序为:
所述第一旋转门、第二旋转门和目标旋转门;
或者,
在所述L个训练层的任一训练层中包含有所述目标旋转门和所述第二旋转门的情况下,各旋转门的作用顺序为:第二旋转门和目标旋转门。
9.根据权利要求6所述的方法,其中,所述预设参数化量子电路的L个训练层之后,还包括其他旋转门。
10.根据权利要求6所述的方法,其中,所述目标量子电路包含有M层,所述M为大于等于1小于等于的正整数;
所述M层中的至少一层是基于如下所得:
将两个训练层中第一训练层的第一目标旋转门替换第一受控酉门,将所述两个训练层中第二训练层的第二目标旋转门替换第二受控酉门;其中,所述两个训练层为所述L个训练层中的任意两个训练层。
11.根据权利要求10所述的方法,其中,所述两个训练层为个训练层中的任意相邻的两个训练层。
12.根据权利要求10所述的方法,其中,
在所述目标块编码矩阵能够通过泡利串的组合表示的情况下,所述目标量子电路中所述第一受控酉门为所述泡利串的组合所对应的等效电路,以及所述目标量子电路中所述第二受控酉门为所述泡利串的组合的共轭转置所对应的等效电路;
或者,
在所述目标块编码矩阵不能够通过泡利串的组合表示的情况下,所述目标量子电路中所述第一受控酉门为用于实现所述目标块编码矩阵的第一等效电路,以及所述目标量子电路中所述第二受控酉门为用于实现所述目标块编码矩阵的共轭转置的第二等效电路;
其中,所述第一等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第三受控酉门,第三受控酉门为第一块编码矩阵对应的酉门;
受控于所述辅助寄存器和所述第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第四受控酉门,第四受控酉门为所述第一块编码矩阵的共轭转置对应的酉门;
其中,所述第二等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和所述第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第四受控酉门;
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第三受控酉门;
其中,所述第一组块编码量子比特和第二组块编码量子比特组成所述块编码寄存器;所述第一块编码矩阵为所述目标哈密顿量的块编码矩阵,所述第一块编码矩阵所包含的量子比特的数量小于所述m+n。
13.一种哈密顿量模拟装置,包括:
参数调整单元,用于确定目标量子电路的子电路中目标可调参数的目标参数值;其中,所述目标参数值为在目标演化时间下所对应的、且满足第一误差条件的参数值;所述第一误差条件为预设的误差容忍值,用于约束目标参数值的精度,以约束得到的目标输出量子态的精度;所述目标量子电路包含有一个辅助量子比特所组成的辅助寄存器、m个块编码量子比特所组成的块编码寄存器和n个主量子比特所组成的主寄存器;所述子电路作用于所述辅助寄存器;所述目标量子电路中还包含有受控于所述辅助寄存器且作用于所述块编码寄存器和所述主寄存器的目标受控酉门,所述目标受控酉门用于模拟目标哈密顿量所对应的目标块编码矩阵;所述目标受控酉门包括与目标块编码矩阵等价的第一受控酉门,和与所述目标块编码矩阵的共轭转置等价的第二受控酉门;所述目标块编码矩阵所对应的量子比特的数量为m+n,所述n为所述目标哈密顿量所对应的量子比特的数量;
输出单元,用于在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,得到所述目标量子电路中所述主寄存器的目标输出量子态,其中,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间下所对应的目标量子态。
14.根据权利要求13所述的装置,其中,所述输出单元,还用于:
在所述目标可调参数为所述目标参数值、所述辅助寄存器的第一输入态为预设初态、所述块编码寄存器的第二输入态为所述预设初态,以及所述主寄存器的第三输入态为所述目标哈密顿量所对应的第一量子态的情况下,获取所述目标量子电路中所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息;
在所述辅助寄存器的第一状态信息和所述块编码寄存器的第二状态信息均满足预设条件的情况下,得到所述主寄存器的目标输出量子态,所述目标输出量子态为模拟所述目标哈密顿量在所述目标演化时间下所对应的目标量子态。
15.根据权利要求13所述的装置,其中,所述参数调整单元,具体用于:
将训练完成的预设参数化量子电路中所述目标可调参数的目标参数值,作为所述子电路中目标可调参数的目标参数值;所述训练完成的预设参数化量子电路用于模拟目标函数;所述目标函数用于表征演化时间与自变量之间关联关系;
其中,所述目标量子电路是通过如下方式得到:将所述预设参数化量子电路中的量子比特作为辅助寄存器,并扩展出块编码寄存器和主寄存器,将所述预设参数化量子电路中作用于所述辅助寄存器的第一目标旋转门替换为所述第一受控酉门,以及将所述预设参数化量子电路中作用于所述辅助寄存器的第二目标旋转门替换为所述第二受控酉门;
其中,所述第一目标旋转门的第一旋转参数和第二目标旋转门的第二旋转参数均为所述目标函数的自变量;所述子电路包含有所述预设参数化量子电路中除所述第一目标旋转门和第二目标旋转门之外的至少部分电路。
16.根据权利要求15所述的装置,其中,所述参数调整单元,还用于:
在所述预设参数化量子电路的旋转参数取值为N个数据点中的任一数据点的情况下,获取所述预设参数化量子电路的实际输出态;其中,所述实际输出态为在预设初态作为输入态的情况下,所述预设参数化量子电路在所述目标可调参数为当前参数值下的输出态;
基于所述实际输出态与所述预设初态得到实际输出结果;其中,所述实际输出结果为所述实际输出态与预设输入态之间的内积;所述N为大于等于1的正整数,所述j为大于等于1小于等于N的正整数;所述旋转参数包括所述第一旋转参数和所述第二旋转参数;
得到N个实际输出结果;
在确定满足迭代终止条件的情况下,将所述目标可调参数的当前参数值作为训练完成的所述预设参数化量子电路中所述目标可调参数的目标参数值;
其中,所述迭代终止条件包括以下至少之一:
基于所述N个实际输出结果与N个目标输出结果,确定预设损失函数的损失值满足收敛条件;所述目标输出结果;
当前迭代次数到达预设次数。
17.根据权利要求16所述的装置,其中,所述参数调整单元,还用于:
在确定不满足所述迭代终止条件的情况下,对所述目标可调参数的参数值进行调整;
重新在所述预设参数化量子电路的旋转参数取值为N个数据点中的任一数据点的情况下,获取所述预设参数化量子电路的实际输出态,以得到实际输出结果;
重新得到N个实际输出结果,直至满足所述迭代终止条件为止。
18.根据权利要求15-17任一项所述的装置,其中,所述预设参数化量子电路包括有L个训练层;所述L为大于等于2的偶数,所述L的取值与所述第一误差条件有关;
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数用于对第一角度进行旋转操作;其中,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;
用于对第二角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第一旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第一旋转门的旋转角度和所述第二旋转门的旋转角度为所述目标可调参数;
或者,
所述L个训练层中的至少两个训练层包括:
目标旋转门,所述旋转参数用于对第一角度进行旋转操作;其中,所述第一目标旋转门和第二目标旋转门为处于不同训练层中的目标旋转门;
用于对第三角度进行旋转操作的、且作用于所述预设参数化量子电路中量子比特上的第二旋转门;
其中,所述第二旋转门的旋转角度为所述目标可调参数。
19.根据权利要求18所述的装置,其中,还满足以下至少之一:
所述第一角度为z轴所对应的角度;
所述第二角度为z轴所对应的角度;
所述第三角度为y轴所对应的角度。
20.根据权利要求18所述的装置,其中,
在所述L个训练层的任一训练层中包含有所述目标旋转门、所述第一旋转门以及所述第二旋转门的情况下,各旋转门的作用顺序为:
所述第一旋转门、第二旋转门和目标旋转门;
或者,
在所述L个训练层的任一训练层中包含有所述目标旋转门和所述第二旋转门的情况下,各旋转门的作用顺序为:第二旋转门和目标旋转门。
21.根据权利要求18所述的装置,其中,所述预设参数化量子电路的L个训练层之后,还包括其他旋转门。
22.根据权利要求18所述的装置,其中,所述目标量子电路包含有M层,所述M为大于等于1小于等于的正整数;
所述M层中的至少一层是基于如下所得:
将两个训练层中第一训练层的第一目标旋转门替换第一受控酉门,将所述两个训练层中第二训练层的第二目标旋转门替换第二受控酉门;其中,所述两个训练层为所述L个训练层中的任意两个训练层。
23.根据权利要求22所述的装置,其中,所述两个训练层为个训练层中的任意相邻的两个训练层。
24.根据权利要求22所述的装置,其中,
在所述目标块编码矩阵能够通过泡利串的组合表示的情况下,所述目标量子电路中所述第一受控酉门为所述泡利串的组合所对应的等效电路,以及所述目标量子电路中所述第二受控酉门为所述泡利串的组合的共轭转置所对应的等效电路;
或者,
在所述目标块编码矩阵不能够通过泡利串的组合表示的情况下,所述目标量子电路中所述第一受控酉门为用于实现所述目标块编码矩阵的第一等效电路,以及所述目标量子电路中所述第二受控酉门为用于实现所述目标块编码矩阵的共轭转置的第二等效电路;
其中,所述第一等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第三受控酉门,第三受控酉门为第一块编码矩阵对应的酉门;
受控于所述辅助寄存器和所述第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第四受控酉门,第四受控酉门为所述第一块编码矩阵的共轭转置对应的酉门;
其中,所述第二等效电路,按作用顺序至少包括:
受控于所述辅助寄存器和所述第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第四受控酉门;
受控于所述辅助寄存器和第一组块编码量子比特,且作用于第二组块编码量子比特和所述主寄存器的第三受控酉门;
其中,所述第一组块编码量子比特和第二组块编码量子比特组成所述块编码寄存器;所述第一块编码矩阵为所述目标哈密顿量的块编码矩阵,所述第一块编码矩阵所包含的量子比特的数量小于所述m+n。
25.一种计算设备,包括:
至少一个量子处理单元;
存储器,耦合到所述至少一个QPU并用于存储可执行指令,
所述指令被所述至少一个量子处理单元执行,以使所述至少一个量子处理单元能够执行权利要求1-12中任一项所述的方法;
或者,包括:
至少一个处理器;以及
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1-12中任一项所述的方法。
26.一种存储有计算机指令的非瞬时计算机可读存储介质,其特征在于,当至少一个量子处理单元执行时,所述计算机指令使得所述至少一个量子处理单元执行根据权利要求1-12任一项所述的方法;
或者,所述计算机指令用于使所述计算机执行根据权利要求1-12中任一项所述的方法。
27.一种计算机程序产品,包括计算机程序,所述计算机程序在被至少一个量子处理单元执行时实现根据权利要求1-12中任一项所述的方法;
或者所述计算机程序在被处理器执行时实现根据权利要求1-12中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211196971.9A CN115577790B (zh) | 2022-09-28 | 2022-09-28 | 哈密顿量模拟方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211196971.9A CN115577790B (zh) | 2022-09-28 | 2022-09-28 | 哈密顿量模拟方法、装置、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115577790A CN115577790A (zh) | 2023-01-06 |
CN115577790B true CN115577790B (zh) | 2024-08-09 |
Family
ID=84582571
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211196971.9A Active CN115577790B (zh) | 2022-09-28 | 2022-09-28 | 哈密顿量模拟方法、装置、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115577790B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116306950B (zh) * | 2023-02-15 | 2024-08-09 | 北京百度网讯科技有限公司 | 基态特征的确定方法、装置、设备及存储介质 |
CN116468126B (zh) * | 2023-04-06 | 2024-07-09 | 北京邮电大学 | 基于量子梯度下降的求解组合优化问题的迭代量子算法 |
CN116913402B (zh) * | 2023-09-14 | 2024-01-05 | 国开启科量子技术(北京)有限公司 | 基于量子线路模拟生物光电解制氢气的装置、方法、设备及介质 |
CN117951595B (zh) * | 2024-03-27 | 2024-06-18 | 苏州元脑智能科技有限公司 | 生物数据分类方法、装置、电子设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112561068A (zh) * | 2020-12-10 | 2021-03-26 | 北京百度网讯科技有限公司 | 模拟方法、计算设备、经典设备、存储设备及产品 |
CN113379058A (zh) * | 2021-06-08 | 2021-09-10 | 北京百度网讯科技有限公司 | 量子模拟方法及装置、电子设备及存储介质 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11120359B2 (en) * | 2019-03-15 | 2021-09-14 | Microsoft Technology Licensing, Llc | Phase estimation with randomized hamiltonians |
US11562282B2 (en) * | 2020-03-05 | 2023-01-24 | Microsoft Technology Licensing, Llc | Optimized block encoding of low-rank fermion Hamiltonians |
CN114418105B (zh) * | 2020-10-28 | 2023-08-08 | 本源量子计算科技(合肥)股份有限公司 | 一种基于量子线路处理量子应用问题的方法及装置 |
CN114580647B (zh) * | 2022-02-24 | 2023-08-01 | 北京百度网讯科技有限公司 | 量子系统的模拟方法、计算设备、装置及存储介质 |
CN114662694A (zh) * | 2022-03-31 | 2022-06-24 | 北京百度网讯科技有限公司 | 量子系统的特征信息的确定方法、装置、设备及存储介质 |
CN114925840B (zh) * | 2022-05-31 | 2024-07-23 | 北京百度网讯科技有限公司 | 模拟方法、设备及存储介质 |
-
2022
- 2022-09-28 CN CN202211196971.9A patent/CN115577790B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112561068A (zh) * | 2020-12-10 | 2021-03-26 | 北京百度网讯科技有限公司 | 模拟方法、计算设备、经典设备、存储设备及产品 |
CN113379058A (zh) * | 2021-06-08 | 2021-09-10 | 北京百度网讯科技有限公司 | 量子模拟方法及装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN115577790A (zh) | 2023-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115577790B (zh) | 哈密顿量模拟方法、装置、设备及存储介质 | |
CN114580647B (zh) | 量子系统的模拟方法、计算设备、装置及存储介质 | |
CN113011593B (zh) | 消除量子测量噪声的方法及系统、电子设备和介质 | |
US20230054391A1 (en) | Calibration of quantum measurement device | |
CN114418103B (zh) | 基态能量的确定方法、装置、设备及存储介质 | |
CN114418107B (zh) | 酉算子编译方法、计算设备、装置及存储介质 | |
CN115577787B (zh) | 量子振幅估计方法、装置、设备以及存储介质 | |
CN115577776B (zh) | 基态能量的确定方法、装置、设备及存储介质 | |
CN113255922A (zh) | 量子纠缠量化方法和装置、电子设备、计算机可读介质 | |
CN115456189B (zh) | 量子模拟方法、装置、设备及存储介质 | |
CN114580645A (zh) | 随机量子测量的模拟方法、装置、设备及存储介质 | |
CN115577789B (zh) | 量子纠缠程度确定方法、装置、设备及存储介质 | |
CN115577782B (zh) | 量子计算方法、装置、设备及存储介质 | |
CN115577783B (zh) | 量子数据处理方法、装置、设备以及存储介质 | |
CN115577786B (zh) | 量子熵确定方法、装置、设备及存储介质 | |
CN115577781B (zh) | 量子相对熵确定方法、装置、设备及存储介质 | |
CN115456184B (zh) | 量子电路处理方法、量子态制备方法、装置、设备及介质 | |
CN116108926B (zh) | 量子计算方法、装置、设备及存储介质 | |
CN116306950B (zh) | 基态特征的确定方法、装置、设备及存储介质 | |
CN116432766B (zh) | 非局域量子操作的模拟方法、装置、设备以及存储介质 | |
CN115577788B (zh) | 量子熵确定方法、装置、设备及存储介质 | |
CN116090573B (zh) | 非局域量子操作的模拟方法、装置、设备以及存储介质 | |
CN116451794B (zh) | 可蒸馏纠缠的估计方法、装置、设备及存储介质 | |
CN116523065B (zh) | 确定量子设备演化酉矩阵的方法及装置、电子设备和介质 | |
CN112580803B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |