CN114492815B - 一种基于量子化学计算目标体系能量的方法、装置及介质 - Google Patents
一种基于量子化学计算目标体系能量的方法、装置及介质 Download PDFInfo
- Publication number
- CN114492815B CN114492815B CN202210103684.2A CN202210103684A CN114492815B CN 114492815 B CN114492815 B CN 114492815B CN 202210103684 A CN202210103684 A CN 202210103684A CN 114492815 B CN114492815 B CN 114492815B
- Authority
- CN
- China
- Prior art keywords
- target system
- quantum
- average energy
- energy
- expectation
- 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 60
- 238000012360 testing method Methods 0.000 claims abstract description 82
- 238000004364 calculation method Methods 0.000 claims abstract description 52
- 238000005259 measurement Methods 0.000 claims abstract description 28
- 238000013507 mapping Methods 0.000 claims description 18
- 238000013461 design Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 15
- 230000005284 excitation Effects 0.000 claims description 13
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 description 20
- 238000004088 simulation Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 13
- 229910052739 hydrogen Inorganic materials 0.000 description 13
- 239000001257 hydrogen Substances 0.000 description 13
- 239000002096 quantum dot Substances 0.000 description 12
- 150000002431 hydrogen Chemical class 0.000 description 10
- 238000011161 development Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000005283 ground state Effects 0.000 description 7
- 239000002245 particle Substances 0.000 description 7
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 230000001131 transforming effect Effects 0.000 description 6
- 230000005428 wave function Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000005610 quantum mechanics Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000005381 potential energy Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000005366 Ising model Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000004776 molecular orbital Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000005359 quantum Heisenberg model Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000005303 weighing Methods 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于量子化学计算目标体系能量的方法、装置及介质，方法包括：获取待求解目标体系的试验态，并测量试验态的平均能量期望，判断平均能量期望是否满足目标体系能量的计算终止条件，其中，计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度，若是，则将当前平均能量期望作为待求解目标体系的能量，否则，更新试验态，测量更新后的当前试验态的平均能量期望，继续执行判断平均能量期望是否满足目标体系能量的计算终止条件的步骤，直至获取满足终止条件的待求解目标体系的能量，它能够为量子化学模拟计算目标体系能量的实现提供支持，提高计算速度和计算精度，促进量子化学模拟应用的进一步发展。
Description
技术领域
本发明属于量子计算技术领域，特别是一种基于量子化学计算目标体系能量的方法、装置及介质。
背景技术
量子计算机是一类遵循量子力学规律进行高速数学和逻辑运算、存储及处理量子信息的物理装置。当某个装置处理和计算的是量子信息，运行的是量子算法时，它就是量子计算机。量子计算机因其具有相对普通计算机更高效的处理数学问题的能力，例如，能将破解RSA密钥的时间从数百年加速到数小时，故成为一种正在研究中的关键技术。
量子计算模拟是一个借助数值计算和计算机科学来仿真遵循量子力学规律的模拟计算，作为一个仿真程序，它依据量子力学的量子比特的基本定律，利用计算机的高速计算能力，刻画量子态的时空演化。
随着量子化学理论的不断完善，计算化学已经成了化学工作者解释实验现象、预测实验结果、指导实验设计的重要工具，在药物的合成、催化剂的制备等方面有着广泛的应用。但是，面对计算化学所涉及的巨大计算量，经典计算机在计算精度、计算尺寸等方面显得能力有限，这就在一定程度上限制了计算化学的发展，由此导致用户对化学体系进行模拟计算的应用不强，影响量子化学模拟应用的进一步发展。
发明内容
本发明的目的是提供一种基于量子化学计算目标体系能量的方法、装置及介质，以解决现有技术中的不足，它能够为量子化学模拟计算目标体系能量的实现提供支持，提高了计算速度和计算精度，促进量子化学模拟应用的进一步发展。
本申请的一个实施例提供了一种基于量子化学计算目标体系能量的方法，所述方法包括：
获取待求解目标体系的试验态，并测量所述试验态的平均能量期望；
判断所述平均能量期望是否满足所述目标体系能量的计算终止条件，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度；
若是，则将当前所述平均能量期望作为所述待求解目标体系的能量，否则，更新所述试验态，测量更新后的当前试验态的平均能量期望，继续执行所述判断所述平均能量期望是否满足所述目标体系能量的计算终止条件的步骤，直至获取满足所述终止条件的所述待求解目标体系的能量。
可选的，所述获取待求解目标体系的试验态，包括：
根据待求解目标体系的电子数和轨道信息，获取所述目标体系的Hartree Fock态；
根据所述目标体系的Hartree Fock态，获取待求解目标体系的试验态。
可选的，所述根据所述目标体系的Hartree Fock态，获取待求解目标体系的试验态，包括：
根据预先选择的拟设方式，对所述目标体系的Hartree Fock态进行演化，得到演化后的量子态作为待求解目标体系的试验态。
可选的，所述根据预先选择的拟设方式，对所述目标体系的Hartree Fock态进行演化，得到演化后的量子态作为待求解目标体系的试验态，包括：
根据预先选择的拟设方式，计算所述目标体系费米子形式的簇算符；
选择映射方式并将所述目标体系费米子形式的簇算符变换为泡利算子形式的簇算符；
将所述泡利算子形式的簇算符分解为对应酉算子形式并进行演化，得到演化后的量子态作为待求解目标体系的试验态。
可选的，所述拟设方式包括：单激发耦合簇或单双激发耦合簇。
可选的，所述映射方式为Jordan－Wigner变换、Parity变换、Bravyi－Kitaev变换和SegmentParity变换中的一种。
可选的，所述测量所述试验态的平均能量期望，包括：
所述测量所述试验态的平均能量期望，包括：
获取所述目标体系对应的费米子哈密顿量，并将所述目标体系对应的费米子哈密顿量变换为所述目标体系的泡利哈密顿量；
根据所述目标体系的泡利哈密顿量分解的各个子项，构建所述目标体系的泡利哈密顿量各个子项对应的量子线路；
利用所述目标体系的泡利哈密顿量各个子项对应的量子线路，测量所述试验态的平均能量期望。
本申请的又一实施例提供了一种基于量子化学计算目标体系能量的装置，所述装置包括：
获取模块，用于获取待求解目标体系的试验态，并测量所述试验态的平均能量期望；
判断模块，用于判断所述平均能量期望是否满足所述目标体系能量的计算终止条件，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度；
更新模块，用于若是，则将当前所述平均能量期望作为所述待求解目标体系的能量，否则，更新所述试验态，测量更新后的当前试验态的平均能量期望，继续执行所述判断所述平均能量期望是否满足所述目标体系能量的计算终止条件的步骤，直至获取满足所述终止条件的所述待求解目标体系的能量。
可选的，所述获取模块，包括：
第一获取单元，用于根据待求解目标体系的电子数和轨道信息，获取所述目标体系的Hartree Fock态；
第二获取单元，用于根据所述目标体系的Hartree Fock态，获取待求解目标体系的试验态。
可选的，所述第二获取单元，包括：
演化单元，用于根据预先选择的拟设方式，对所述目标体系的Hartree Fock态进行演化，得到演化后的量子态作为待求解目标体系的试验态。
可选的，所述演化单元，包括：
计算单元，用于根据预先选择的拟设方式，计算所述目标体系费米子形式的簇算符；
第一变换单元，用于选择映射方式并将所述目标体系费米子形式的簇算符变换为泡利算子形式的簇算符；
分解单元，用于将所述泡利算子形式的簇算符分解为对应酉算子形式并进行演化，得到演化后的量子态作为待求解目标体系的试验态。
可选的，所述获取模块，包括：
第二变换单元，获取所述目标体系对应的费米子哈密顿量，并将所述目标体系对应的费米子哈密顿量变换为所述目标体系的泡利哈密顿量；
构建单元，用于根据所述目标体系的泡利哈密顿量分解的各个子项，构建所述目标体系的泡利哈密顿量各个子项对应的量子线路；
测量单元，用于利用所述目标体系的泡利哈密顿量各个子项对应的量子线路，测量所述试验态的平均能量期望。
本申请的一个实施例提供了一种存储介质，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行上述任一项所述的方法。
本申请的一个实施例提供了一种电子装置，包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行上述任一项所述的方法。
本申请的又一实施例提供了一种量子计算机操作系统，所述量子计算机操作系统根据上述任一项中所述方法实现基于量子化学计算目标体系的能量。
本申请的又一实施例提供了一种量子计算机，所述量子计算机包括所述的量子计算机操作系统。
与现有技术相比，本发明首先获取待求解目标体系的试验态，并测量试验态的平均能量期望，判断平均能量期望是否满足目标体系能量的计算终止条件，若是，则将当前平均能量期望作为待求解目标体系的能量，否则，更新试验态，测量更新后的当前试验态的平均能量期望，继续执行判断平均能量期望是否满足目标体系能量的计算终止条件的步骤，直至获取满足终止条件的待求解目标体系的能量，它能够为量子化学模拟计算目标体系能量的实现提供支持，提高计算速度和计算精度，促进量子化学模拟应用的进一步发展。
附图说明
图1为本发明实施例提供的一种基于量子化学计算目标体系能量的方法的计算机终端的硬件结构框图；
图2为本发明实施例提供的一种基于量子化学计算目标体系能量方法的流程示意图；
图3为本发明实施例提供的一种拟设方式对应的量子线路结构示意图；
图4为本发明实施例提供的一种根据泡利算子形式簇算符构建对应的量子线路示意图；
图5为本发明实施例提供的一种构建氢分子泡利哈密顿量各个子项对应的量子线路示意图；
图6为本发明实施例提供的一种氢分子泡利哈密顿量各个子项期望展开形式的测量线路示意图；
图7为本发明实施例提供的一种基于量子化学计算目标体系能量装置的结构示意图。
具体实施方式
下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
本发明实施例首先提供了一种基于量子化学计算目标体系能量的方法，该方法可以应用于电子设备，如计算机终端，具体如普通电脑、量子计算机等。
下面以运行在计算机终端上为例对其进行详细说明。图1为本发明实施例提供的一种基于量子化学计算目标体系能量的方法的计算机终端的硬件结构框图。如图1所示，计算机终端可以包括一个或多个(图1中仅示出一个)处理器102(处理器102可以包括但不限于微处理器MCU或可编程逻辑器件FPGA等的处理装置)和用于存储数据的存储器104，可选地，上述计算机终端还可以包括用于通信功能的传输装置106以及输入输出设备108。本领域普通技术人员可以理解，图1所示的结构仅为示意，其并不对上述计算机终端的结构造成限定。例如，计算机终端还可包括比图1中所示更多或者更少的组件，或者具有与图1所示不同的配置。
存储器104可用于存储应用软件的软件程序以及模块，如本申请实施例中的基于量子化学计算目标体系能量的方法对应的程序指令/模块，处理器102通过运行存储在存储器104内的软件程序以及模块，从而执行各种功能应用以及数据处理，即实现上述的方法。存储器104可包括高速随机存储器，还可包括非易失性存储器，如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中，存储器104可进一步包括相对于处理器102远程设置的存储器，这些远程存储器可以通过网络连接至计算机终端。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
传输装置106用于经由一个网络接收或者发送数据。上述的网络具体实例可包括计算机终端的通信供应商提供的无线网络。在一个实例中，传输装置106包括一个网络适配器(Network Interface Controller，NIC)，其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中，传输装置106可以为射频(Radio Frequency，RF)模块，其用于通过无线方式与互联网进行通讯。
需要说明的是，真正的量子计算机是混合结构的，它包含两大部分：一部分是经典计算机，负责执行经典计算与控制；另一部分是量子设备，负责运行量子程序进而实现量子计算。而量子程序是由量子语言如QRunes语言编写的一串能够在量子计算机上运行的指令序列，实现了对量子逻辑门操作的支持，并最终实现量子计算。具体的说，量子程序就是一系列按照一定时序操作量子逻辑门的指令序列。
在实际应用中，因受限于量子设备硬件的发展，通常需要进行量子计算模拟以验证量子算法、量子应用等等。量子计算模拟即借助普通计算机的资源搭建的虚拟架构(即量子虚拟机)实现特定问题对应的量子程序的模拟运行的过程。通常，需要构建特定问题对应的量子程序。本发明实施例所指量子程序，即是经典语言编写的表征量子比特及其演化的程序，其中与量子计算相关的量子比特、量子逻辑门等等均有相应的经典代码表示。
量子线路作为量子程序的一种体现方式，也称量子逻辑电路，是最常用的通用量子计算模型，表示在抽象概念下对于量子比特进行操作的线路，其组成包括量子比特、线路(时间线)，以及各种量子逻辑门，最后常需要通过量子测量操作将结果读取出来。
不同于传统电路是用金属线所连接以传递电压信号或电流信号，在量子线路中，线路可看成是由时间所连接，亦即量子比特的状态随着时间自然演化，在这过程中按照哈密顿运算符的指示，一直到遇上逻辑门而被操作。
一个量子程序整体上对应有一条总的量子线路，本发明所述量子程序即指该条总的量子线路，其中，该总的量子线路中的量子比特总数与量子程序的量子比特总数相同。可以理解为：一个量子程序可以由量子线路、针对量子线路中量子比特的测量操作、保存测量结果的寄存器及控制流节点(跳转指令)组成，一条量子线路可以包含几十上百个甚至千上万个量子逻辑门操作。量子程序的执行过程，就是对所有的量子逻辑门按照一定时序执行的过程。需要说明的是，时序即单个量子逻辑门被执行的时间顺序。
需要说明的是，经典计算中，最基本的单元是比特，而最基本的控制模式是逻辑门，可以通过逻辑门的组合来达到控制电路的目的。类似地，处理量子比特的方式就是量子逻辑门。使用量子逻辑门，能够使量子态发生演化，量子逻辑门是构成量子线路的基础，量子逻辑门包括单比特量子逻辑门，如Hadamard门(H门，阿达马门)、泡利－X门(X门)、泡利－Y门(Y门)、泡利－Z门(Z门)、RX门、RY门、RZ门等等；多比特量子逻辑门，如CNOT门、CR门、iSWAP门、Toffoli门等等。量子逻辑门一般使用酉矩阵表示，而酉矩阵不仅是矩阵形式，也是一种操作和变换。一般量子逻辑门在量子态上的作用是通过酉矩阵左乘以量子态右矢对应的矩阵进行计算的。
量子态，即量子比特的逻辑状态，在量子算法(或称量子程序)中用二进制表示，例如，一组量子比特为q0、q1、q2，表示第0位、第1位、第2位量子比特，从高位到低位排序为q2q1q0，该组量子比特对应的量子态是该组量子比特对应的本征态的叠加，该组量子比特对应的本征态共有2的量子比特总数次方个，即8个本征态(确定的状态)：|000>、|001>、|010>、|011>、|100>、|101>、|110>、|111>，每个本征态的位与量子比特对应一致，如|000>态，000从高位到低位对应q2q1q0，|>为狄拉克符号。
以单个量子比特说明，单个量子比特的逻辑状态可能处于|0>态、|1>态、|0>态和|1>态的叠加态(不确定状态)，具体可以表示为/>其中，c和d为表示量子态振幅(概率幅)的复数,振幅模的平方|c|2和|d|2分别表示|0>态、|1>态的概率，|c|2+|d|2＝1。简言之，量子态是各本征态组成的叠加态，当其它本征态的概率为0时，即处于唯一确定的本征态。
长期以来，在原子水平上对分子和材料的能量和性质的理论解释一直被认为是量子计算最直接的应用之一，其作为一种新的计算范式受到广泛关注。与经典计算相比，量子计算的计算能力随量子比特数的增加呈指数增长。随着不断的发展，在许多领域都出现突破性的进展，包括制药、光伏、航空、电子和能源发电等。量子计算机最有可能的应用之一是模拟量子系统，其中分子是自然界常见的量子系统，计算分子系统的能量是量子化学的主要目标之一。
参见图2，图2为本发明实施例提供的一种基于量子化学计算目标体系能量的方法的流程示意图。
本实施例提供一种基于量子化学计算目标体系能量的方法的一实施例，所述基于量子化学计算目标体系能量的方法包括：
S201：获取待求解目标体系的试验态，并测量所述试验态的平均能量期望。
具体的，获取待求解目标体系的试验态，可以包括：
1.根据待求解目标体系的电子数和轨道信息，获取所述目标体系的Hartree Fock态。
首先对于待求解的目标体系，电子数，就是目标体系包含的电子数量，电子是一种基本粒子，一般是指目标体系的核外电子的数目；轨道信息，是以数学的方法描述目标体系在原子核外的特定空间中，找到电子的几率，并指出电子在三维空间中的可能位置。
示例性的，对于待求解氢分子目标体系，其含有四个单电子自旋分子轨道和两个电子，根据氢分子电子数和轨道信息，若用一个量子比特代表一个自旋分子轨道，即0表示空轨道，1表示占据轨道，因此氢分子目标体系的Hartree Fock(哈特里-福克)态可以用量子态|0011>来表示。
对于待求解氢分子目标体系，只需要分别在两个量子比特上添加一个NOT门，就可以在量子线路中将|0000>初始化成|0011>。因此对于任意一个含有M个自旋分子轨道的N电子体系，它的Hartree Fock态都可以这样简单的表示，只要在量子线路中给定M个量子比特，然后在前N个量子线路上添加NOT门即可得到所需要的N电子体系的Hartree Fock态。
需要说明的是，在量子计算中，波函数的选取需要一个参考波函数作为基矢，比如量子化学中一般使用Hartree Fock态矢作为参考波函数，以满足：
ψ(θ)＝U(θ)|ψ>Hartree-Fock
其中，ψ(θ)表示在一组参数集合θ下对应的波函数，U(θ)表示在一组参数集合θ下对应的矩阵算子，参考波函数|ψ>Hartree-Fock对应化学中Hartree Fock基态，表示分子的电子都处在最低轨道处。
2.根据所述目标体系的Hartree Fock态，获取待求解目标体系的试验态。
具体的，根据预先选择的拟设方式，对所述目标体系的Hartree Fock态进行演化，得到演化后的量子态作为待求解目标体系的试验态，可包括如下步骤：
步骤1：根据预先选择的拟设方式，计算所述目标体系费米子形式的簇算符。
具体的，簇算符可以理解为是一种人为定义的算符种类，用于表示轨道上电子的跳跃。拟设是一种将制备好的分子初态，例如|ψ>Hartree-Fock演化到量子线路上的方法，预先选择的拟设方式可为耦合簇法(Coupled Cluster，CC)，它是一种从Hartree Fock分子轨道出发，通过拟设得到试验态|ψ＞的方法。这里的拟设为指数耦合簇算符eT，满足：|ψ＞＝eT|ψ＞Hartree-Fock，拟设中的T就是N电子簇算符，其定义式为若干激发算符之和，即：
T＝T1+T2+...+TN
其中，T1是单粒子激发算符，T2是双粒子激发算符，其余项以此类推。由于在一个多电子体系中，三激发、四激发发生的概率很小，所以通常在双激发处进行“截断”，最终只剩T1和T2两项，即：
T＝T1+T2
其中，为创造算符，ar、as为湮灭算符，p、q、r、s代表轨道，这里的待定系数tpq、tpqrs需要通过优化器来寻找的参数/>满足/>
需要注意的是，通过拟设方式将目标体系的初态转化为费米子形式的簇算符之后，由于eT指数耦合簇算符不是酉算子，因此无法直接将eT指数耦合簇算符通过预设映射方式映射到量子比特上，无法构建出对应的量子线路，所以，需要构造出酉算子版本的指数耦合簇算符，即酉耦合簇算符(Unitary Coupled Cluster，UCC)。
示例性的，首先可定义一个等效的厄米哈密顿量令/>然后，以/>为生成元生成UCC算符：/>其中，若UCC中的簇算符T只含有T1这一项，则称这一拟设为单激发耦合簇(UCCS)；若UCC中的簇算符T含有T1和T2两项，则称这个拟设为单双激发耦合簇(UCCSD)。
对应的，对于UCCS和UCCSD，拟设对应的量子线路相同，例如为图3所示，图3为一种拟设方式对应的量子线路结构示意图，具体为UCC方法对应的四比特的量子线路示意图，图示为4量子比特q0、q1、q2、q3的量子线路示意图，其中，X-π/2、Xπ/2表示参数分别为-π/2、π/2的X门，Y门，同理，图标及其与实心连线表示CNOT门，Zθ表示参数为θ的Z门。显示的拟设原理可以包括：拟设公式，例如可以为量子线路对应的矩阵算子U(θ)。对于UCC，对应拟设公式为：
其中，即为拟设，Pi为生成元。
可选的，对于描述氢分子的|0011>Hartree-Fock态，此时的簇算符T正是费米子哈密顿量，即
当T＝T1时，即由前四项单激发所构造成的哈密顿量；当T＝T1+T2时，即由单激发和双激发所共同构造成的哈密顿量。
步骤2：选择映射方式并将所述目标体系费米子形式的簇算符变换为泡利算子形式的簇算符。
具体的，映射方式可为Jordan-Wigner变换、Parity变换、Bravyi-Kitaev变换和SegmentParity变换中的一种。
本领域技术人员能够理解的是，每种映射方式对应的映射原理可以包括：态映射原理和算符映射原理，例如，对于Jordan－Wigner变换，显示的态映射为：
其中，表示量子比特的计算态，/>表示变换矩阵，/>表示费米子体系的占据态。显示的算符映射为：
其中，表示升降算符，j表示量子比特序号，P表示宇称集，ZP(j)表示作用属于宇称集P的量子比特上的一组泡利Z矩阵，X表示泡利X矩阵，Y表示泡利Y矩阵。
同等的，算符映射也可显示为：
其中，表示产生算符，aj表示湮灭算符，/>与aj统称费米子体系的升降算符，表示量子比特上的产生算符/湮灭算符，/>表示宇称算子，n表示量子比特数。
其他变换的态映射和算符映射显示方式与Jordan－Wigner变换的原理相同，在此不再赘述。
在一种可选的方式中，若将费米子形式的簇算符通过Jordan－Wigner变换变成泡利算子形式时，它是若干子项的和，表达式为：
其中，σ为泡利算子，α、β∈(X，Y，Z，I)，而i、j则表示簇算符子项所作用的子空间，h为实数。
步骤3：将所述泡利算子形式的簇算符分解为对应酉算子形式并进行演化，得到演化后的量子态作为待求解目标体系的试验态。
具体的，接上述示例，将费米子形式的簇算符通过Jordan-Wigner变换变成泡利算子形式时，它是若干子项的和，表达式为：
但是，如果对这些子项进行求和，最后得到的泡利算子形式簇算符想要对角化以生成酉算子，是比较困难的。因此为了能够以每个子项Hk为生成元将UCC算符分解成有限个酉算子来进行模拟，有必要引进渐进近似定理，即托特公式(Trotter Foluma)，该定理是量子仿真算法的核心：limn→∞(eiAt/neiBt/n)n＝ei(A+B)t，其中，A、B均为厄米算符，t为实数，n为正整数。
需要说明的是，通过托特公式，指数函数就可以分解为数个子指数函数项的近似形式。托特分解强调的是当n取的越大，其越接近于原式的趋势，而非具体考虑n取何值。
示例性的，假设经Jordan-Wigner变换后的泡利算子形式的簇算符T的表达式如下：
根据托特公式，构造泡利算子型哈密顿量H对应的量子线路，可通过逐项模拟哈密顿量，即先对H1项进行模拟：
通过推导，得出在q0量子比特上直接添加RZ门即可模拟H1项。
对于H2、H3、H4、H5项进行模拟，可以参照H1项，得到：
U2(H2，θ2)＝CNOT(0，1)RZ(1，2θ2)CNOT(0，1)
U3(H3，θ3)＝CNOT(0，2)CNOT(1，2)RZ(2，2θ3)CNOT(1，2)CNOT(0，2)
U4(H4，θ4)＝H(0)CNOT(0，1)RZ(1，2θ4)CNOT(0，1)H(0)
那么，模拟T的量子线路构造如图4所示，最终可根据如图4所示的一种根据泡利算子形式簇算符构建对应的量子线路示意
图测量所述试验态的平均能量期望。
其中，所述测量所述试验态的平均能量期望，包括：
步骤a：获取所述目标体系对应的费米子哈密顿量，并将所述目标体系对应的费米子哈密顿量变换为所述目标体系的泡利哈密顿量。
具体的，哈密顿量是所有粒子的动能的总和加上与系统相关的粒子的势能。对于不同的情况或数量的粒子，哈密顿量是不同的，因为它包括粒子的动能之和以及对应于这种情况的势能函数，一般用H表示。在量子力学中，经典力学的物理量变为相应的算符，哈密顿量对应的正是哈密顿算符。
具体的，基于目标系统的力学分析，就可以获得此系统的哈密顿量，获取目标体系对应的费米子哈密顿量需要借助创造算符和湮灭算符aq来实现，它们满足反对易关系。
示例性的，对于氢分子体系，其对应的费米子哈密顿量为：
需要说明的是，在量子计算中，费米子形式的哈密顿量是无法直接在线路上演化的，所以需要存在将积分形式的期望值求解转换成量子线路可读的过程，这个过程称作映射。需要注意的是，映射仅仅是将哈密顿量换个形式表达，而各类型哈密顿量所表示的系统能量信息是等价的。此外，对于一个量子模拟线路或者真实量子芯片来说，泡利算符是更加容易操作和生成，所以可以将目标体系对应的费米子哈密顿量变换为目标体系的泡利哈密顿量，便于后续的模拟操作。
接上述示例，对于氢分子体系，将其对应的费米子哈密顿量变换为泡利哈密顿量，具体为：
步骤b：根据所述目标体系的泡利哈密顿量分解的各个子项，构建所述目标体系的泡利哈密顿量各个子项对应的量子线路。
具体的，获取待求解目标体系的试验态|ψn>后，需要开始利用量子期望估计算法来计算试验态|ψn>在分子哈密顿量上的期望。所谓量子期望估计，是指对于多电子体系、Heisenberg模型(海森堡模型)、量子Ising模型(易辛模型)等体系的哈密顿量H可以展开成多个子项的和，即：
其中，h为实数，σ为泡利算子，α、β和γ∈(X，Y，Z，I)，而i、j、k表示哈密顿量子项所作用的子空间。
由于可观测量是线性的，因此在利用下式计算体系的平均能量时：
E＝<ψ*|H|ψ>
其中，ψ*与ψ是正交归一的，等式右边也可以展开成这种形式：
由此可知，只须先对每个子项求期望，然后对各个期望求和，就能得到体系的平均能量E。需要说明的是，每个子项期望的测量可以在量子处理器上进行，利用经典处理器可负责对各个期望进行求和。
示例性的，假设某一个体系的哈密顿量为H，它最终可以展开成这种形式：
在该式中，所有子项系数h均是1，并假设所获取的试验态为如下形式：
|ψ>＝a|00>+b|01>+c|10>+d|11>
其中，a2、b2、c2、d2分别是指测量试验态时，坍塌到|00>、|01>、|10>、|11>的概率PS，将哈密顿量的各个子项H1、H2、H3分别作用于试验态上，可以依次得到期望E1、E2、E3，具体的：
E1＝<ψ*|H1|ψ>
E2＝<ψ*|H2|ψ>
E3＝<ψ*|H3|ψ>
以E1、E2、E3为例，对于期望E1，系数h就是期望，无须构造线路测量，即对于期望E2，其哈密顿量为/>由于测量操作是在σZ上(以σZ的特征向量为基向量所构成的子空间)进行的，所以只需要在量子比特上加上测量门即可，然后将测量结果传递给经典处理器求和。
接上述示例，以氢分子为例，/> 其泡利哈密顿量可以分解为15个子项，分别构造整个氢分子哈密顿量Hp的15个子项的测量线路，得到如图5所示的一种构建氢分子泡利哈密顿量各个子项对应的量子线路示意图，即可获得各个子项的期望E(i)。
步骤c：利用所述目标体系的泡利哈密顿量各个子项对应的量子线路，测量所述试验态的平均能量期望。
具体的，将目标体系的泡利哈密顿量各个子项期望测量线路展开，即可得到各个子项期望E(i)的测量线路，随后量子处理器将E(i)依次传给经典处理器求和，就得到了目标体系在该试验态下的平均能量期望。
需要说明的是，由于测量操作是在σZ上进行的，对于包含σx、σy的哈密顿量，此时不能直接测量，需要对σx和σy进行换基操作，也就是让试验态再演化一次，由于σx＝H×σZ×H，即对于σx和σy，在测量前，需要分别在相应的量子比特上添上Hadamard门和/>门，然后将测量结果传递给经典处理器求和即可。
示例性的，接上述氢分子的示例，将氢分子体系的泡利哈密顿量各个子项期望测量线路展开，即可得到如图6所示的一种氢分子泡利哈密顿量各个子项期望展开形式的测量线路示意图，接着，量子处理器将E(i)依次传给经典处理器求和，就得到了氢分子在该试验态下的平均能量E(n)，且
S202：判断所述平均能量期望是否满足所述目标体系能量的计算终止条件，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度。
对于描述某一目标体系(如多电子体系)的哈密顿量的特征值E1E2...En，进而求得目标体系的能量E0，通过利用该目标体系的哈密顿量作用于试验态时，可以得到该体系在这一状态下的平均能量E，它将大于或接近于体系的基态能量E0，即：
从上述表达式中可以看出，如果所获取的试验态|ψ>正好就是体系的基态|ψ0>，那么不等式中的等号成立，直接得到了目标体系的基态能量E0；但往往更多的情况是，获取的试验态|ψ>与目标体系的基态相比有一定差距，导致计算得到的E大于E0很多，这时就需要引入一组参数通过不断调节/>来更新试验态，使其最终非常接近目标体系的基态能量。
具体的，通过获取待求解目标体系的试验态并测量试验态/>的平均能量En，判断当前平均能量期望与前次测量后的平均能量期望的差值符合精度，其中精度可以由用户根据计算需求自行设定。
S203：若是，则将当前所述平均能量期望作为所述待求解目标体系的能量，否则，更新所述试验态，测量更新后的当前试验态的平均能量期望，继续执行所述判断所述平均能量期望是否满足所述目标体系能量的计算终止条件的步骤，直至获取满足所述终止条件的所述待求解目标体系的能量。
具体的，若待求解目标体系的试验态对应的平均能量期望满足计算终止条件，则所获取的试验态正好就是体系的基态，直接得到了目标体系的能量E0；否则优化器会利用梯度无关算法，如Nelder－Mead算法或梯度相关算法，如梯度下降法等等优化参数然后传给量子处理器，继续演化和测量，通过不断迭代参数/>来更新试验态，使其最终获取满足终止条件的待求解目标体系的能量。
例如，当获取待求解目标体系的试验态为第一次演化时，则前次演化不存在，可以默认前次演化后测量的试验态的平均能量为0，直接进入下一次迭代；在当前平均能量期望与前次测量后的平均能量期望的差值不符合精度时，利用优化方法，对拟设中量子线路的可调参数进行优化，更新优化后试验态，测量更新后的当前试验态的平均能量期望，返回执行S202的步骤，直到某次演化后的能量差值符合精度，并将该次演化后的试验态下的平均能量期望确定为目标体系对应能量。
需要强调的是，以上提出的拟设方式、映射方式及优化方法等仅仅作为示例，并不构成对本发明的限定，例如，拟设方式还包括HE(Hardware Efficient，硬件高效)、SP(Symmetry Preserved，对称保持)等方法。
可见，本发明首先获取待求解目标体系的试验态，并测量试验态的平均能量期望，判断平均能量期望是否满足目标体系能量的计算终止条件，若是，则将当前平均能量期望作为待求解目标体系的能量，否则，更新试验态，测量更新后的当前试验态的平均能量期望，继续执行判断平均能量期望是否满足目标体系能量的计算终止条件的步骤，直至获取满足终止条件的待求解目标体系的能量，它能够为量子化学模拟计算目标体系能量的实现提供支持，提高计算速度和计算精度，促进量子化学模拟应用的进一步发展。
参见图7，图7为本发明实施例提供的一种基于量子化学计算目标体系能量装置的结构示意图，与图2所示的流程相对应，所述装置包括：
获取模块701，用于获取待求解目标体系的试验态，并测量所述试验态的平均能量期望；
判断模块702，用于判断所述平均能量期望是否满足所述目标体系能量的计算终止条件，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度；
更新模块703，用于若是，则将当前所述平均能量期望作为所述待求解目标体系的能量，否则，更新所述试验态，测量更新后的当前试验态的平均能量期望，继续执行所述判断所述平均能量期望是否满足所述目标体系能量的计算终止条件的步骤，直至获取满足所述终止条件的所述待求解目标体系的能量。
具体的，所述获取模块，包括：
第一获取单元，用于根据待求解目标体系的电子数和轨道信息，获取所述目标体系的Hartree Fock态；
第二获取单元，用于根据所述目标体系的Hartree Fock态，获取待求解目标体系的试验态。
具体的，所述第二获取单元，包括：
演化单元，用于根据预先选择的拟设方式，对所述目标体系的Hartree Fock态进行演化，得到演化后的量子态作为待求解目标体系的试验态。
具体的，所述演化单元，包括：
计算单元，用于根据预先选择的拟设方式，计算所述目标体系费米子形式的簇算符；
第一变换单元，用于选择映射方式并将所述目标体系费米子形式的簇算符变换为泡利算子形式的簇算符；
分解单元，用于将所述泡利算子形式的簇算符分解为对应酉算子形式并进行演化，得到演化后的量子态作为待求解目标体系的试验态。
具体的，所述获取模块，包括：
第二变换单元，获取所述目标体系对应的费米子哈密顿量，并将所述目标体系对应的费米子哈密顿量变换为所述目标体系的泡利哈密顿量；
构建单元，用于根据所述目标体系的泡利哈密顿量分解的各个子项，构建所述目标体系的泡利哈密顿量各个子项对应的量子线路；
测量单元，用于利用所述目标体系的泡利哈密顿量各个子项对应的量子线路，测量所述试验态的平均能量期望。
与现有技术相比，本发明首先获取待求解目标体系的试验态，并测量试验态的平均能量期望，判断平均能量期望是否满足目标体系能量的计算终止条件，若是，则将当前平均能量期望作为待求解目标体系的能量，否则，更新试验态，测量更新后的当前试验态的平均能量期望，继续执行判断平均能量期望是否满足目标体系能量的计算终止条件的步骤，直至获取满足终止条件的待求解目标体系的能量，它能够为量子化学模拟计算目标体系能量的实现提供支持，提高计算速度和计算精度，促进量子化学模拟应用的进一步发展。
本发明实施例还提供了一种存储介质，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行上述任一项方法实施例中的步骤。
具体的，在本实施例中，上述存储介质可以被设置为存储用于执行以下步骤的计算机程序：
S201：获取待求解目标体系的试验态，并测量所述试验态的平均能量期望；
S202：判断所述平均能量期望是否满足所述目标体系能量的计算终止条件，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度；
S203：若是，则将当前所述平均能量期望作为所述待求解目标体系的能量，否则，更新所述试验态，测量更新后的当前试验态的平均能量期望，继续执行所述判断所述平均能量期望是否满足所述目标体系能量的计算终止条件的步骤，直至获取满足所述终止条件的所述待求解目标体系的能量。
具体的，在本实施例中，上述存储介质可以包括但不限于：U盘、只读存储器(Read－Only Memory，简称为ROM)、随机存取存储器(Random Access Memory，简称为RAM)、移动硬盘、磁碟或者光盘等各种可以存储计算机程序的介质。
本发明实施例还提供了一种电子装置，包括存储器和处理器，其特征在于，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行上述任一项方法实施例中的步骤。
具体的，上述电子装置还可以包括传输设备以及输入输出设备，其中，该传输设备和上述处理器连接，该输入输出设备和上述处理器连接。
具体的，在本实施例中，上述处理器可以被设置为通过计算机程序执行以下步骤：
S201：获取待求解目标体系的试验态，并测量所述试验态的平均能量期望；
S202：判断所述平均能量期望是否满足所述目标体系能量的计算终止条件，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度；
S203：若是，则将当前所述平均能量期望作为所述待求解目标体系的能量，否则，更新所述试验态，测量更新后的当前试验态的平均能量期望，继续执行所述判断所述平均能量期望是否满足所述目标体系能量的计算终止条件的步骤，直至获取满足所述终止条件的所述待求解目标体系的能量。
本发明实施例还提供了一种量子计算机操作系统，所述量子计算机操作系统根据本发明实施例中提供的上述任一方法实施例实现基于量子化学计算目标体系的能量。
本申请的实施例还提供了一种量子计算机，所述量子计算机包括所述的量子计算机操作系统。
以上依据图式所示的实施例详细说明了本发明的构造、特征及作用效果，以上所述仅为本发明的较佳实施例，但本发明不以图面所示限定实施范围，凡是依照本发明的构想所作的改变，或修改为等同变化的等效实施例，仍未超出说明书与图示所涵盖的精神时，均应在本发明的保护范围内。
Claims (10)
1.一种基于量子化学计算目标体系能量的方法，其特征在于，所述方法包括：
根据预先选择的拟设方式，对待求解目标体系的态进行演化，得到演化后的量子态作为所述目标体系的试验态，并测量所述试验态的平均能量期望，其中，所述目标体系的/>态根据目标体系的电子数和轨道信息确定，所述演化包括：计算所述目标体系费米子形式的簇算符，选择映射方式并将所述目标体系费米子形式的簇算符变换为泡利算子形式的簇算符，将所述泡利算子形式的簇算符分解为对应酉算子形式并进行演化；
响应于当前所述平均能量期望满足所述目标体系能量的计算终止条件，以将当前所述平均能量期望作为所述目标体系的能量，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度。
2.根据权利要求1所述的方法，其特征在于，所述方法还包括：
若当前所述平均能量期望不满足所述目标体系能量的计算终止条件，更新所述试验态，测量更新后的当前试验态的平均能量期望，直至获取满足所述终止条件的所述目标体系的能量。
3.根据权利要求1所述的方法，其特征在于，所述拟设方式包括：单激发耦合簇或单双激发耦合簇。
4.根据权利要求3所述的方法，其特征在于，所述映射方式为Jordan-Wigner变换、Parity变换、Bravyi-Kitaev变换和SegmentParity变换中的一种。
5.根据权利要求4所述的方法，其特征在于，所述测量所述试验态的平均能量期望，包括：
获取所述目标体系对应的费米子哈密顿量，并将所述目标体系对应的费米子哈密顿量变换为所述目标体系的泡利哈密顿量；
根据所述目标体系的泡利哈密顿量分解的各个子项，构建所述目标体系的泡利哈密顿量各个子项对应的量子线路；
利用所述目标体系的泡利哈密顿量各个子项对应的量子线路，测量所述试验态的平均能量期望。
6.一种基于量子化学计算目标体系能量的装置，其特征在于，所述装置包括：
获取模块，用于根据预先选择的拟设方式，对待求解目标体系的态进行演化，得到演化后的量子态作为所述目标体系的试验态，并测量所述试验态的平均能量期望，其中，所述目标体系的/>态根据目标体系的电子数和轨道信息确定，所述演化包括：计算所述目标体系费米子形式的簇算符，选择映射方式并将所述目标体系费米子形式的簇算符变换为泡利算子形式的簇算符，将所述泡利算子形式的簇算符分解为对应酉算子形式并进行演化；
判断模块，用于响应于当前所述平均能量期望满足所述目标体系能量的计算终止条件，以将当前所述平均能量期望作为所述目标体系的能量，其中，所述计算终止条件为当前平均能量期望与前次测量后的平均能量期望的差值符合精度。
7.一种存储介质，其特征在于，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行所述权利要求1至5任一项所述的方法。
8.一种电子装置，包括存储器和处理器，其特征在于，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行所述权利要求1至5任一项所述的方法。
9.一种量子计算机操作系统，其特征在于，所述量子计算机操作系统根据权利要求1至5任一项所述的方法实现基于量子化学计算目标体系的能量。
10.一种量子计算机，其特征在于，所述量子计算机包括权利要求9所述的量子计算机操作系统。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210103684.2A CN114492815B (zh) | 2022-01-27 | 2022-01-27 | 一种基于量子化学计算目标体系能量的方法、装置及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210103684.2A CN114492815B (zh) | 2022-01-27 | 2022-01-27 | 一种基于量子化学计算目标体系能量的方法、装置及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114492815A CN114492815A (zh) | 2022-05-13 |
| CN114492815B true CN114492815B (zh) | 2023-08-08 |
Family
ID=81476352
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210103684.2A Active CN114492815B (zh) | 2022-01-27 | 2022-01-27 | 一种基于量子化学计算目标体系能量的方法、装置及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114492815B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115527626B (zh) * | 2022-08-16 | 2023-04-25 | 腾讯科技（深圳）有限公司 | 分子处理方法、装置、电子设备、存储介质及程序产品 |
| CN116340714B (zh) * | 2023-03-10 | 2024-07-16 | 本源量子计算科技（合肥）股份有限公司 | 一种张量数据处理方法及相关装置 |
| CN116612823B (zh) * | 2023-05-12 | 2024-09-03 | 合肥本源量子计算科技有限责任公司 | 一种基态能量的确定方法和系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109379187A (zh) * | 2018-12-03 | 2019-02-22 | 三峡大学 | 一种基于氢原子的完全量子通信方法 |
| CN110637310A (zh) * | 2017-05-15 | 2019-12-31 | 谷歌有限责任公司 | 量子计算系统内的算子平均 |
| CN111599414A (zh) * | 2020-03-25 | 2020-08-28 | 清华大学 | 一种基于量子计算机的全量子分子模拟方法 |
| CN111712839A (zh) * | 2018-01-30 | 2020-09-25 | 河道研究有限公司 | 确定态能的方法 |
| EP3718060A1 (en) * | 2017-11-30 | 2020-10-07 | 1QB Information Technologies Inc. | Methods and systems for quantum computing enabled molecular ab initio simulations using quantum-classical computing hardware |
| CN113496285A (zh) * | 2021-07-07 | 2021-10-12 | 北京百度网讯科技有限公司 | 基于量子电路的数据处理方法及装置、电子设备和介质 |
-
2022
- 2022-01-27 CN CN202210103684.2A patent/CN114492815B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110637310A (zh) * | 2017-05-15 | 2019-12-31 | 谷歌有限责任公司 | 量子计算系统内的算子平均 |
| EP3718060A1 (en) * | 2017-11-30 | 2020-10-07 | 1QB Information Technologies Inc. | Methods and systems for quantum computing enabled molecular ab initio simulations using quantum-classical computing hardware |
| CN111712839A (zh) * | 2018-01-30 | 2020-09-25 | 河道研究有限公司 | 确定态能的方法 |
| CN109379187A (zh) * | 2018-12-03 | 2019-02-22 | 三峡大学 | 一种基于氢原子的完全量子通信方法 |
| CN111599414A (zh) * | 2020-03-25 | 2020-08-28 | 清华大学 | 一种基于量子计算机的全量子分子模拟方法 |
| CN113496285A (zh) * | 2021-07-07 | 2021-10-12 | 北京百度网讯科技有限公司 | 基于量子电路的数据处理方法及装置、电子设备和介质 |
Non-Patent Citations (1)
| Title |
|---|
| 孔祥宇 ; 朱垣晔 ; 闻经纬 ; 辛涛 ; 李可仁 ; 龙桂鲁 ; .核磁共振量子信息处理研究的新进展.物理学报.2018,(第22期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114492815A (zh) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114492814B (zh) | 基于量子计算模拟目标体系能量的方法、装置及介质 | |
| JP7471736B2 (ja) | 量子系の基底状態エネルギーの推定方法、およびシステム | |
| CN114492815B (zh) | 一种基于量子化学计算目标体系能量的方法、装置及介质 | |
| CN114528996B (zh) | 一种目标体系试验态初始参数的确定方法、装置及介质 | |
| CN113379058B (zh) | 量子模拟方法及装置、电子设备及存储介质 | |
| CN114512194B (zh) | 用于量子化学模拟中获取目标体系试验态的方法及装置 | |
| CN115116552B (zh) | 一种量子化学模拟的分子能量计算方法及装置 | |
| CN114519429B (zh) | 获取目标体系的可观测量的方法、装置及介质 | |
| CN115618663B (zh) | 一种网格方程与物理方程耦合的量子求解方法及装置 | |
| CN117313873A (zh) | 基于量子计算的化学反应路径的生成方法、装置及介质 | |
| CN117634620A (zh) | 一种静电势计算方法、装置、存储介质及电子装置 | |
| CN114512193B (zh) | 基于自旋对称性和等同粒子特性制备体系试验态的方法 | |
| CN117151231A (zh) | 利用变分量子线路求解线性系统的方法、装置及介质 | |
| CN114511090B (zh) | 基于自旋对称性制备体系试验态的方法、装置及介质 | |
| CN114511091B (zh) | 基于等同粒子特性制备体系试验态的方法、装置及介质 | |
| CN117114119A (zh) | 一种计算目标体系能量的方法、装置及介质 | |
| CN116842704A (zh) | 一种组合物最优结构的确定方法、装置及介质 | |
| CN116862004A (zh) | 利用量子化学计算目标体系能量的方法、装置及介质 | |
| CN116976447A (zh) | 一种试验波函数的确定方法、装置及介质 | |
| CN117010516A (zh) | 一种目标体系哈密顿量的确定方法及装置 | |
| CN117669757A (zh) | 一种哈密顿量的构造方法及装置 | |
| CN118782161A (zh) | 一种电子数守恒约束的基态能量计算方法及相关装置 | |
| CN117951429A (zh) | 一种偏微分方程组的求解方法及装置 | |
| CN118227928A (zh) | 一种基于量子计算的线性方程组的求解方法及装置 | |
| CN117669758A (zh) | 基于分布式vqe的量子模拟方法、装置及存储介质 |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 230088 6th floor, E2 building, phase II, innovation industrial park, 2800 innovation Avenue, Hefei high tech Zone, Hefei City, Anhui Province Applicant after: Benyuan Quantum Computing Technology (Hefei) Co.,Ltd. Address before: 230088 6th floor, E2 building, phase II, innovation industrial park, 2800 innovation Avenue, Hefei high tech Zone, Hefei City, Anhui Province Applicant before: ORIGIN QUANTUM COMPUTING COMPANY, LIMITED, HEFEI |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |