CN114511093B - Simulation method of boson subsystem - Google Patents

Simulation method of boson subsystem Download PDF

Info

Publication number
CN114511093B
CN114511093B CN202011281177.5A CN202011281177A CN114511093B CN 114511093 B CN114511093 B CN 114511093B CN 202011281177 A CN202011281177 A CN 202011281177A CN 114511093 B CN114511093 B CN 114511093B
Authority
CN
China
Prior art keywords
hilbert
pauli
space
boson
operators
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
Application number
CN202011281177.5A
Other languages
Chinese (zh)
Other versions
CN114511093A (en
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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202011281177.5A priority Critical patent/CN114511093B/en
Publication of CN114511093A publication Critical patent/CN114511093A/en
Application granted granted Critical
Publication of CN114511093B publication Critical patent/CN114511093B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N10/00Quantum computing, i.e. information processing based on quantum-mechanical phenomena
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Artificial Intelligence (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Evolutionary Computation (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a simulation method of a boson subsystem, which assumes that the maximum occupation number cutoff of the boson subsystem is 2 n -1, the method comprising: establishing n-bit qubit basis vectors and 2 n A code mapping relation of a space base vector of the Hilbert is maintained; will 2 n The generated operators after the space interception of the Hilbert are expressed by Pauli matrixes; pauli matrix representation to generate operators from 2 n Vehicrt space recurrence to 2 n+1 The Hilbert space is maintained. The invention uses 2 brought by the quantum bit entanglement and state superposition principle n The space complexity of the simulation of the boson subsystem is reduced to a logarithmic relationship by the aid of the space characteristics of the Hilbert, so that the simulation method can be used for a quantum computer, saves quantum bit computing resources, reduces computing difficulty, can be used for environments such as a quantum virtual machine of a classical computer, and solves the problem of memory index growth of the simulation boson subsystem of the classical computer.

Description

Simulation method of boson subsystem
Technical Field
The invention belongs to the field of quantum computation, and particularly relates to a boson subsystem simulation method.
Background
Quantum simulation algorithms are an important class of quantum algorithms, with fermi and boson quantum simulation algorithms being particularly important. For the fermi simulation algorithm, there is a well-known Jordan-Wigner transform and a modified Bravyi-Kitaev transform. For quantum simulation of fermi or bosons, the Pauli matrix expression generating operators can be used for constructing quantum circuits for simulating related problems by using the existing method so as to further calculate the whole set of quantum simulation algorithm of the fermi or bosons. The Jordan-Wigner transform, the Bravyi-Kitaev transform, is Pauli matrix expression that gives two different Fermi-generating operators. In the prior art, a boson simulation algorithm generally adopts an algorithm based on one-to-one correspondence between quantum bits and occupied basis vectors. When the algorithm is executed on a quantum computer, the required quantum bit number is O (N) (N is the maximum occupancy number of bosons), and on a quantum virtual machine of a classical computer, the required calculation resource grows exponentially along with the increase of the occupancy number, so that the classical computer cannot complete the calculation when the occupancy number is large.
Disclosure of Invention
Based on the above, the invention aims to solve the problem of large consumption of computing resources in the boson simulation algorithm in the prior art by establishing a new space mapping relation between the quantum bits and Hilbert and defining a Pauli matrix expression of a new boson generating operator on the basis of the space mapping relation.
The invention aims at realizing the following technical scheme: a simulation method of a boson subsystem, wherein the maximum occupation number cutoff of the boson subsystem is 2 n -1, the method comprising:
establishing n-bit qubit basis vectors and 2 n A code mapping relation of a space base vector of the Hilbert is maintained;
the 2 n The generated operators after the space interception of the Hilbert are expressed by Pauli matrixes;
expressing the Pauli matrix of the generated operator from 2 n Vehicrt space recurrence to 2 n+1 The Hilbert space is maintained.
Specifically, the n-bit qubit basis vector is combined with 2 n The coding mapping relation of the space base vector of the Hilbert is as follows
0 >=|0 1 ,0 2 ,...,0 n-1 ,0 n >,
1 >=|0 1 ,0 2 ,...,0 n-1 ,1 n >,
2 >=|0 1 ,0 2 ,...,1 n-1 ,0 n >,
3 >=|0 1 ,0 2 ,...,1 n-1 ,1 n >,
...
Figure BDA0002780835320000029
wherein ,|φi >(i=0,1,2,...,2 n -1) is represented at 2 n Basic vector with the number of bosons in Hilbert space being i, i 0>And |1>Respectively representing two bases of the qubit.
Specifically, said 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000021
wherein ,
Figure BDA0002780835320000022
is said 2 n Generating operators after space truncation of the Hilbert>
Figure BDA0002780835320000023
Is the effect operator of the ith item on the 1 st to nth quantum bits, and the effect operator is a tensor product formed by certain arrangement and combination of elements in the Pauli combination matrix;
the Pauli combined matrix is Pauli matrix sigma xy ,σ z The Pauli combination matrix includes at least four elements:
Figure BDA0002780835320000024
specifically, the generation algorithmPauli matrix expression of symbols from 2 n Vehicrt space recurrence to 2 n+1 The method for maintaining Hilbert space is that
Will be
Figure BDA0002780835320000025
The qubit position acted by each quantum gate is added by 1
Figure BDA0002780835320000026
1 st to 2 nd n -1 is written as
Figure BDA0002780835320000027
Will be 2 n The term is written as
Figure BDA0002780835320000028
Will be 2 n +1 to 2 n+1 -1 is written as
Figure BDA0002780835320000031
Said 2 n+1 The generation operator after the space truncation of the Hilbert is the addition of the above items
Figure BDA0002780835320000032
Specifically, when n=1, the ratio of 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000033
specifically, when n=2, the 2 n Dimension HilThe resulting operator after the bert space truncation is expressed as a Pauli matrix:
Figure BDA0002780835320000034
similarly, the technical scheme of the invention can write Pauli matrix expression of any truncated quantum bit number to generate an operator.
The technical scheme disclosed by the invention has the following beneficial effects:
1) The method reduces the number of qubits required for the simulation calculation of the boson subsystem with the maximum occupation number truncated to N from O (N) to O (log) 2 N);
2) The method can be simultaneously applied to quantum simulation and classical numerical simulation based on a quantum virtual machine.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the following detailed description of the present invention will be made with reference to examples. It should be understood that the examples described herein are for illustrative purposes only and are not intended to limit the scope of the present invention.
The invention discloses a simulation method of a boson subsystem, which assumes that the maximum occupation number cutoff of the boson subsystem is 2 n -1, the method comprising:
s001, establishing n-bit quantum bit basic vector and 2 n A code mapping relation of a space base vector of the Hilbert is maintained;
s002. will 2 n The generated operators after the space interception of the Hilbert are expressed by Pauli matrixes;
s003, expressing Pauli matrix for generating operator from 2 n Vehicrt space recurrence to 2 n+1 The Hilbert space is maintained.
Specifically, in embodiments of the present invention, an n-bit qubit basis vector is combined with 2 n The coding mapping relation of the space base vector of the Hilbert is as follows
0 >=|0 1 ,0 2 ,...,0 n-1 ,0 n >,
1 >=|0 1 ,0 2 ,...,0 n-1 ,1 n >,
2 >=|0 1 ,0 2 ,...,1 n-1 ,0 n >,
3 >=|0 1 ,0 2 ,...,1 n-1 ,1 n >,
...
Figure BDA0002780835320000045
wherein ,|φi >(i=0,1,2,...,2 n -1) is represented at 2 n Basic vector with the number of bosons in Hilbert space being i, i 0>And |1>Respectively representing two bases of the qubit.
In particular, 2 n The matrix form of the generated operators after the dimensional Hilbert spatial truncation can be written as:
Figure BDA0002780835320000041
wherein N=2n -1。
Specifically, in the embodiment of the present invention, 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000042
wherein ,
Figure BDA0002780835320000043
is said 2 n Generating operators after space truncation of the Hilbert>
Figure BDA0002780835320000044
Is the effect operator of the ith item on the 1 st to nth quantum bits, and the effect operator is a tensor product formed by certain arrangement and combination of elements in the Pauli combination matrix;
pauli combined matrix is Pauli matrix sigma xy ,σ z The Pauli combination matrix includes at least four elements:
Figure BDA0002780835320000051
three Pauli matrices sigma xy ,σ z The following relationship is satisfied:
σ x σ y =iσ zy σ z =iσ xz σ x =iσ yx σ y =-σ y σ xy σ z =-σ z σ yz σ x =-σ x σ zx σ x =σ y σ y =σ z σ z =I;
wherein I is a 2×2 identity matrix
Figure BDA0002780835320000052
Pauli matrix and identity matrix I form a closed algebraic system, and can be used as the basis of hermite 2X 2 matrix space.
Specifically, in embodiments of the present invention, pauli matrix expression to generate operators is from 2 n Vehicrt space recurrence to 2 n+1 The method for maintaining Hilbert space is that
Will be
Figure BDA0002780835320000053
The qubit position acted by each quantum gate is added by 1
Figure BDA0002780835320000054
1 st to 2 nd n -1 is written as
Figure BDA0002780835320000055
Will be 2 n The term is written as
Figure BDA0002780835320000056
Will be 2 n +1 to 2 n+1 -1 is written as
Figure BDA0002780835320000057
2 n+1 The generation operator after the space truncation of the Hilbert is the addition of the above items
Figure BDA0002780835320000058
Specifically, in the embodiment of the present invention, when n=1, 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000061
when n=2, 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000062
when n=3, 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000063
when n=4, said 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000064
when n=5, said 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000065
/>
Figure BDA0002780835320000071
when n=6, said 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure BDA0002780835320000072
/>
Figure BDA0002780835320000081
the truncated generating operators given by the Pauli combination matrix can be simplified and the same class items can be obtained to obtain Pauli matrix expression of the generating operators. For example when n=3,
Figure BDA0002780835320000082
/>
Figure BDA0002780835320000091
in particular, according to the recurrence relation disclosed in the present invention and the Pauli expression of the generating operator after n=1 truncation, in principle, any Pauli matrix expression of the generating operator after n-value truncation and the corresponding annihilation operator can be written, and will not be described in detail herein. Specifically, in quantum simulation calculations, any expression of an operator combined by producing annihilation operators can be deduced from Pauli expressions that produce annihilation operators. In an embodiment of the present invention, some common operators combined by generating annihilation operators when n=3 are given below.
Truncated particle arithmetic symbols are expressed as
Figure BDA0002780835320000092
The square expression of truncated particle arithmetic is
Figure BDA0002780835320000093
Other common operator members
Figure BDA0002780835320000094
Finally, it should be noted that the above-mentioned embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made to the technical solution of the present invention without departing from the spirit and scope of the technical solution of the present invention, and all such modifications and equivalents are intended to be encompassed in the scope of the claims of the present invention.

Claims (8)

1. A simulation method of a boson subsystem, wherein the maximum occupation number cutoff of the boson subsystem is 2 n -1, the simulation method comprising, on a quantum computer or a quantum virtual machine of a classical computer, the following steps:
generating an operator using Pauli matrix construction;
building a quantum circuit with the generated operators; and
applying the quantum circuit to a boson initial quantum state to simulate evolution of the boson system; wherein the method for generating operators by Pauli matrix construction comprises the following steps:
establishing n-bit qubit basis vectors and 2 n Constructing 2 by using the coding mapping relation of the space basic vectors of the Hilbert n A dimension Hilbert space;
the 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure FDA0004181132960000011
expressing the Pauli matrix of the generated operator from 2 n Vehicrt space recurrence to 2 n+1 A dimension Hilbert space;
wherein ,
Figure FDA0004181132960000012
is said 2 n Generating operators after space truncation of the Hilbert>
Figure FDA0004181132960000013
Is the effect operator of the ith item on the 1 st to nth quantum bits, and the effect operator is a tensor product formed by certain arrangement and combination of elements in the Pauli combination matrix;
the Pauli combined matrix is Pauli matrix sigma xy ,σ z The Pauli combination matrix includes at least four elements:
Figure FDA0004181132960000014
the Pauli matrix expression of the generated operator is from 2 n Vehicrt space recurrence to 2 n+1 The method for maintaining the Hilbert space comprises the following steps:
will be
Figure FDA0004181132960000015
The qubit position acted by each quantum gate is added by 1
Figure FDA0004181132960000016
1 st to 2 nd n -1 is written as
Figure FDA0004181132960000017
Will be 2 n The term is written as
Figure FDA0004181132960000018
Will be 2 n +1 to 2 n+1 -1 is written as
Figure FDA0004181132960000021
Said 2 n+1 The generation operator after the space truncation of the Hilbert is the addition of the above items
Figure FDA0004181132960000022
2. The boson system simulation method of claim 1, wherein the n-bit qubit basis vector is 2 n The coding mapping relation of the space base vector of the Hilbert is as follows
0 >=|0 1 ,0 2 ,...,0 n-1 ,0 n >,
1 >=|0 1 ,0 2 ,...,0 n-1 ,1 n >,
2 >=|0 1 ,0 2 ,...,1 n-1 ,0 n >,
3 >=|0 1 ,0 2 ,...,1 n-1 ,1 n >,
...
Figure FDA0004181132960000023
wherein ,|φi >(i=0,1,2,...,2 n -1) is represented at 2 n Basic vector with the number of bosons in Hilbert space being i, i 0 j> and |1j >Each represents two bases of a j-th bit (j=1, 2,., n) qubit.
3. The boson system simulation method according to any one of claims 1-2, wherein when n=1, the ratio of 2 is n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure FDA0004181132960000024
4. a boson system simulation method according to any of claims 1-2, characterized in that when n = 2, the 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure FDA0004181132960000025
5. a boson system simulation method according to any of claims 1-2, characterized in that when n=3, the ratio of 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure FDA0004181132960000026
6. the boson system simulation method according to any one of claims 1-2, wherein when n=4, the ratio of 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure FDA0004181132960000031
7. the boson system simulation method according to any one of claims 1-2, wherein when n=5, the ratio of 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure FDA0004181132960000032
Figure FDA0004181132960000041
8. the boson system simulation method according to any one of claims 1-2, wherein when n=6, the ratio of 2 n The resulting operators after the dimensional Hilbert spatial truncation are expressed as Pauli matrices:
Figure FDA0004181132960000042
/>
Figure FDA0004181132960000051
/>
CN202011281177.5A 2020-11-16 2020-11-16 Simulation method of boson subsystem Active CN114511093B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011281177.5A CN114511093B (en) 2020-11-16 2020-11-16 Simulation method of boson subsystem

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011281177.5A CN114511093B (en) 2020-11-16 2020-11-16 Simulation method of boson subsystem

Publications (2)

Publication Number Publication Date
CN114511093A CN114511093A (en) 2022-05-17
CN114511093B true CN114511093B (en) 2023-06-09

Family

ID=81546892

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011281177.5A Active CN114511093B (en) 2020-11-16 2020-11-16 Simulation method of boson subsystem

Country Status (1)

Country Link
CN (1) CN114511093B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018033823A1 (en) * 2016-08-17 2018-02-22 International Business Machines Corporation Efficient reduction of resources for the simulation of fermionic hamiltonians on quantum hardware
CN108290733A (en) * 2015-12-04 2018-07-17 耶鲁大学 Use the technology and related system and method for Bose subpattern progress quantum error correction
WO2018172629A1 (en) * 2017-03-24 2018-09-27 Bull Sas Method for simulating, on a conventional computer, a quantum circuit
WO2019078907A1 (en) * 2017-10-18 2019-04-25 Google Llc Simulation of quantum circuits
WO2019236712A2 (en) * 2018-06-06 2019-12-12 Microsoft Technology Licensing, Llc Layouts for fault-tolerant quantum computers
CN111160560A (en) * 2019-12-31 2020-05-15 合肥本源量子计算科技有限责任公司 Method and system for predicting resources required by analog quantum computation
WO2020151129A1 (en) * 2019-01-25 2020-07-30 合肥本源量子计算科技有限责任公司 Quantum machine learning framework construction method and apparatus, and quantum computer and computer storage medium
CN111738448A (en) * 2020-06-23 2020-10-02 北京百度网讯科技有限公司 Quantum line simulation method, device, equipment and storage medium

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108290733A (en) * 2015-12-04 2018-07-17 耶鲁大学 Use the technology and related system and method for Bose subpattern progress quantum error correction
WO2018033823A1 (en) * 2016-08-17 2018-02-22 International Business Machines Corporation Efficient reduction of resources for the simulation of fermionic hamiltonians on quantum hardware
WO2018172629A1 (en) * 2017-03-24 2018-09-27 Bull Sas Method for simulating, on a conventional computer, a quantum circuit
WO2019078907A1 (en) * 2017-10-18 2019-04-25 Google Llc Simulation of quantum circuits
WO2019236712A2 (en) * 2018-06-06 2019-12-12 Microsoft Technology Licensing, Llc Layouts for fault-tolerant quantum computers
WO2020151129A1 (en) * 2019-01-25 2020-07-30 合肥本源量子计算科技有限责任公司 Quantum machine learning framework construction method and apparatus, and quantum computer and computer storage medium
CN111160560A (en) * 2019-12-31 2020-05-15 合肥本源量子计算科技有限责任公司 Method and system for predicting resources required by analog quantum computation
CN111738448A (en) * 2020-06-23 2020-10-02 北京百度网讯科技有限公司 Quantum line simulation method, device, equipment and storage medium

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Entropy Exchange and Thermodynamic Properties of the Single Ion Cooling Process;Ping-Xing Chen;《entropy》;第1-8页 *
用经典计算机模拟量子计算机;范洪强;胡滨;袁征;;密码学报(03);第23-35页 *
量子纠缠与经典关联在双光子关联成像中的作用;陈平形;《量子光学学报》;第174-179页 *

Also Published As

Publication number Publication date
CN114511093A (en) 2022-05-17

Similar Documents

Publication Publication Date Title
CN113935491B (en) Method, device, equipment, medium and product for obtaining eigenstates of quantum system
Li et al. Chaotifying linear Elman networks
US20230186138A1 (en) Training of quantum neural network
CN111914378B (en) Single-amplitude quantum computing simulation method and device
JP2023541741A (en) Quantum state preparation circuit generation method, device, chip, equipment and program
CN112633508A (en) Quantum line generation method and device, storage medium and electronic device
JP7372996B2 (en) Node grouping method, device and electronic equipment
JP2023545595A (en) Method, apparatus, device, and storage medium for acquiring eigenstates of quantum systems
Zhou et al. Strongly universal Hamiltonian simulators
JP2023521804A (en) quantum data loader
CN114511093B (en) Simulation method of boson subsystem
CN111414961A (en) Task parallel-based fine-grained distributed deep forest training method
Araujo et al. Low-rank quantum state preparation
Tchórzewski et al. Quantum inspired evolutionary algorithm to improve the accuracy of a neuronal model of the electric power exchange
US20030084013A1 (en) Encoding of data and extraneous information into synthetic gene sequences and the retrieval of same
CN114550849A (en) Method for solving chemical molecular property prediction based on quantum graph neural network
CN108615078B (en) Secret data communication method
Qi et al. Qubit neural tree network with applications in nonlinear system modeling
Sarkar et al. Scalable quantum circuits for n-qubit unitary matrices
Zhong et al. Quantum Competition Network Model Based On Quantum Entanglement.
Chakravarty et al. On reductions of self-dual Yang-Mills equations
Gerdt et al. A mathematica package for simulation of quantum computation
CN117291271A (en) Efficient quantum circuit simulation method based on distributed system
D'Ariano Physics as quantum information processing: quantum fields as quantum automata
Ponnambalam et al. Modified ART1 neural networks for cell formation using production data

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Li Desheng

Inventor after: Wu Wei

Inventor after: Zhong Ming

Inventor after: Wu Chunwang

Inventor after: Chen Pingxing

Inventor before: Li Desheng

Inventor before: Wu Wei

Inventor before: Zhong Ming

Inventor before: Wu Chunwang

Inventor before: Chen Pingxing

GR01 Patent grant
GR01 Patent grant