This application claims priority to and the benefit of Korean Patent Application No. 10-2021-0083317 filed in the Korean Intellectual Property Office on Jun. 25, 2021, the entire contents of which are incorporated herein by reference.
The present disclosure relates to a quantum computing device and a method of mitigating a detection crosstalk.
Quantum computing performs quantum algorithms or general quantum information processing and involves a measurement process to read outcomes. When the quantum computing attempts to measure several qubits independently, independent measurements may not be performed due to interactions among measurement devices. In particular, the current quantum technology contains a lot of noise, so independent measurements may not be performed. This is called a detection crosstalk.
Some embodiments may provide a quantum computing device and a method of mitigating a detecting crosstalk, capable of mitigating the detection crosstalk.
According to an embodiment, a quantum computing device including a quantum device including a qubit system, a detector, and a mitigation module may be provided. The detector may perform measurements on a plurality of qubits of the qubit system. The mitigation module may perform quantum pre-processing on the plurality of qubits before performing the measurements on the plurality of qubits in the detector, and perform classical post-processing on a measurement outcome of the plurality of qubits in the detector to mitigate a detection crosstalk included in the measurement outcome.
In some embodiments, the mitigation module may perform the quantum pre-processing by applying a single-qubit gate to each of the plurality of qubits.
In some embodiments, the single-qubit gate may represent a local unitary.
In some embodiments, the mitigation module may perform the classical post-processing through analysis of a quantum measurement operator describing the detector.
In some embodiments, the mitigation module may perform the quantum pre-processing and the classical post-processing, by performing decomposition on a quantum measurement operator describing the detector based on a parameter and an operator to optimize the parameter and the operator.
In some embodiments, the mitigation module may normalize the quantum measurement operator, and decompose the normalized quantum measurement operator as a following equation: (1−∈) V{right arrow over (a)}|{right arrow over (a)}{right arrow over (a)}|V{right arrow over (a)}†+∈P. Here, the ∈ denotes the parameter, the V{right arrow over (a)} denotes a transformation of the quantum pre-processing, the {right arrow over (a)} denotes the measurement outcome of the plurality of qubits, and the P denotes the operator that is determined by the ∈ and the V{right arrow over (a)}.
In some embodiments, the operator may be an operator of unit trace.
In some embodiments, the mitigation module may mitigate the detection crosstalk included in the measurement outcome based on the parameter and the operator.
In some embodiments, the mitigation module may mitigate the detection crosstalk included in the measurement outcome as in a following equation:
Here, the ∈ denotes the parameter, the η denotes a constant determined by the operator, the {right arrow over (a)} denotes the measurement outcome of the plurality of qubits, the Π{right arrow over (a)} denotes an operator describing a quantum measurement of the detector, the tr denotes a trace function, and pe({right arrow over (a)}) denotes a probability of obtaining the {right arrow over (a)}.
In some embodiments, the mitigation module may minimize an upper bound of an error between a measurement outcome without the detection crosstalk and an outcome obtained by approximating the measurement outcome without the detection crosstalk on an assumption that an expected value of a Hermitian operator in the decomposition is a constant. In this case, the Hermitian operator may be defined by a local unitary of the quantum pre-processing and the operator P in the decomposition.
In some embodiments, the mitigation module may optimize the parameter and operator P to minimize the upper bound of the error.
In some embodiments, the upper bound of the error may be defined based on the parameter, a maximum eigenvalue of the operator, and a minimum eigenvalue of the operator.
In some embodiments, the mitigation module may optimize the parameter and the operator P while constraining the operator to a predetermined condition.
In some embodiments, the predetermined condition may be a condition in which the operator is equal to or greater than zero.
According to another embodiment, a quantum computing device including a quantum device including a qubit system, a detector, a plurality of single-qubit gates, and a mitigation module may be provided. The detector may perform measurements on a plurality of qubits of the qubit system. The plurality of single-qubit gates may be positioned in front of the detector, and may be applied to the plurality of qubits, respectively. The mitigation module may perform classical post-processing on a measurement outcome of the plurality of qubits.
In some embodiments, the mitigation module may mitigate a detection crosstalk included in the measurement outcome by performing the classical post-processing.
In some embodiments, the mitigation module may decompose a quantum measurement operator describing the detector based on a parameter and an operator, and optimize the parameter and the operator to mitigate the detection crosstalk.
According to yet another embodiment, a method of mitigating a detection crosstalk in a quantum computing device may be provided. The method may include performing quantum pre-processing on a plurality of qubits, performing measurements on the plurality of qubits on which the quantum pre-processing is performed, and performing classical post-processing on a measurement outcome of the plurality of qubits to mitigate a detection crosstalk included in the measurement outcome.
In some embodiments, performing the quantum pre-processing may include applying a single-qubit gate to each of the plurality of qubits.
In some embodiments, performing the classical post-processing may include decomposing, based on a parameter and an operator, a quantum measurement operator describing a detector that perform the measurements, and optimizing the parameter and the operator to mitigate the detection crosstalk.
In the following detailed description, only certain example embodiments of the present invention have been shown and described, simply by way of illustration. As those skilled in the art would realize, the described embodiments may be modified in various different ways, all without departing from the spirit or scope of the present invention. Accordingly, the drawings and description are to be regarded as illustrative in nature and not restrictive. Like reference numerals designate like elements throughout the specification.
As used herein, the singular forms “a” and “an” are intended to include the plural forms as well, unless the context clearly indicates otherwise.
The sequence of operations or steps is not limited to the order presented in the claims or figures unless specifically indicated otherwise. The order of operations or steps may be changed, several operations or steps may be merged, a certain operation or step may be divided, and a specific operation or step may not be performed.
Referring to
The quantum device 110 may be a device that includes a physical quantum system (e.g., a qubit system) which is used to perform quantum computations or simulations.
The measurement device 120 may communicate with the quantum device 110, and may include a detector 121 that performs measurements on the physical quantum system included in the quantum device 110. For example, the detector 121 may measure each of qubits included in the quantum device 110. In some embodiments, an operator describing the detector 121 may be defined by using quantum detector tomography (QDT). In some embodiments, the detector 121 may include a quantum circuit for measurement of a quantum state. In some embodiments, the detector 121 may be provided for each qubit to be measured.
The measurement device 120 may further include a mitigation module 122 that mitigates a detection crosstalk included in a measurement outcome of the detector 121. In some embodiments, the detector 121 may store the measurement outcome in the storage device, and the mitigation module 122 may retrieve the measurement outcome from the storage device. In some embodiments, the detector 121 may provide the measurement outcome to the mitigation module 122 directly. In some embodiments, the mitigation module 122 may be implemented as a physical processor.
Next, a detection crosstalk in a quantum computing device is described.
In quantum computing, a measurement process for reading out a quantum state, i.e., a measurement outcome, is performed independently for a plurality of qubits. If interactions among detectors exist, the independent measurement may be not performed. An error caused by noise among the detectors when obtaining the measurement outcome by measuring the plurality of qubits is called “detection crosstalk”.
The existence of the detection crosstalk may be investigated after obtaining an operator describing the quantum measurement (positive operator value measurement, POVM), for example, through quantum detector tomography. The POVM for the n-qubit measurements {right arrow over (a)}=a1 a2 . . . an, may be described as in Equation 1. Here, ai is a measurement outcome of qubit qi (i=1, . . . , n).
Π{right arrow over (a)}=Πa
In Equation 1, Πa
When the POVM for the measurement outcome of n qubits cannot be expressed as in Equation 1 (that is, when the POVM for the measurement outcomes of n qubits cannot be factorized in a product of the POVM elements of individual qubits as in Equation 2), it may be said that detection crosstalk exists.
Π{right arrow over (a)}≠Πa
When the POVM obtained through the quantum detector tomography is normalized to {tilde over (Π)}{right arrow over (a)}, the normalized POVM may be expressed as in Equation 3.
In Equation 3, tr[ ] denotes a trace function.
When an optimization process as expressed in Equation 4 is considered, the existence of the detection crosstalk and a value of DC being a positive number in Equation 4 are equivalent. That is, the existence of detection crosstalk may be determined through the value of DC.
Next, the quantum detector tomography is described.
The quantum detector tomography is a method of investigating an operator that describes measurements of a quantum device. In some embodiments, a quantum detector tomography method using a maximum likelihood estimation (MLE) technique may be used. Referring to
{|00|,|11|,|++|,|−−|,|+i+i|,|−i−i|} Equation 5
As shown in
A measurement outcome of the detector is given as a probability. The measurement is performed by using the Z-measurement, i.e., |00|, |11|, and outputs when detected by measurement operators are 0 and 1, respectively. The probability that an arbitrary quantum state ρ of a qubit qi is detected by a detector described by the measurement operator Πa
First, the measurement device may measure the MUB to obtain a probability fn,i that the n-th state ρn is detected by the i-th detector, and calculate a log-likelihood function as expressed in Equation 6 based on the probability fn,i.
In Equation 6, {tilde over (Π)}a
Next, the measurement device finds a maximum value of the log-likelihood function expressed in Equation 6 that satisfies a condition of Equation 7. In this case, an optimization equation that finds the maximum value of the log-likelihood function of Equation 6 satisfying the condition of Equation 7 using a Lagrange multiplier method may be expressed as in Equation 8. In Equation 8, a Lagrange multiplier Ri may be given as in Equation 9.
In Equation 7, I denotes an identity matrix.
{tilde over (Π)}a
In Equation 8, Ri† is a conjugate transpose of Ri.
While repeating the calculation expressed in Equation 8, the measurement device ends the optimization when a difference between {tilde over (Π)}a
Next, a method of mitigating a detection crosstalk according to an embodiment is described.
A measurement device obtains a quantum measurement operator describing a quantum detector through quantum detector tomography and analyzes the quantum measurement operator to perform mitigation of a detection crosstalk. Referring to
In some embodiments, the mitigation module of the measurement device may decompose the quantum measurement operator of the detector to analyze the quantum measurement operator. The mitigation module may decompose the quantum measurement operator, for example, as in Equation 11, based on a parameter ∈ and an operator P. The decomposition may vary depending on the choice of the parameter ∈ and operator P. In some embodiments, the mitigation module may perform the classical post-processing by optimizing the parameter ∈ and the operator P. In some embodiments, the operator P may be an operator of unit trace.
In Equation 11, V{right arrow over (a)} is a local unitary transformation and represents a plurality of gates respectively acting on a plurality of qubits. That is, V{right arrow over (a)}=Va
In Equation 11, V{right arrow over (a)} may be obtained through an optimization process of Equation 12.
The quantum pre-processing uses V{right arrow over (a)} obtained in Equation 12. In some embodiments, V{right arrow over (a)} may be obtained in advance by using the optimization process of Equation 12 through repeated experiments, and V{right arrow over (a)} obtained in advance may be used in the quantum pre-processing.
In some embodiments, as shown in
Equation 13 may be derived by applying V{right arrow over (a)} to both sides of Equation 11.
{tilde over (Π)}{right arrow over (a)}V{right arrow over (a)}†{tilde over (Π)}{right arrow over (a)}V{right arrow over (a)}=(1−∈)|{right arrow over (a)}{right arrow over (a)}|+∈V{right arrow over (a)}†PV{right arrow over (a)} Equation 13
In some embodiments, ∈ in Equation 13 may be a minimum ∈ that satisfies a predetermined constraint, for example, P≥0. Here, P≥0 means that all eigenvalues of P are greater than zero.
Next, the measurement device performs the classical post-processing through analysis of the measurement operator. A probability that the measurement device obtains an outcome {right arrow over (a)} by measuring an arbitrary state ρ through measurements described in Equation 13 (i.e., measurements on n-qubit state) may be given as in Equation 14. In the decomposition of Equation 11, if an expected value of a Hermitian operator V{right arrow over (a)}†PV{right arrow over (a)}, defined by the local unitary V{right arrow over (a)} and the operator P is expressed as q({right arrow over (a)})=tr[V{right arrow over (a)}†PV{right arrow over (a)}], Equation 14 may be expressed as Equation 15.
p
e({right arrow over (a)})=tr[ρV{right arrow over (a)}†Π{right arrow over (a)}V{right arrow over (a)}] Equation 14
(tr[Π{right arrow over (a)}])−1pe({right arrow over (a)})=(1−∈)p0({right arrow over (a)})+∈q({right arrow over (a)}) Equation 15
In Equation 15, p0({right arrow over (a)})=tr[ρ|{right arrow over (a)}{right arrow over (a)}|] is a noise-free measurement outcome wherein the noise may be, for example, a detection crosstalk.
Since a value of the probability distribution q({right arrow over (a)}) is not known for an arbitrary state ρ, in some embodiments, for convenience, q({right arrow over (a)}) may be assumed to be a constant η∈[0,1]. When rewriting Equations 14 and 15 using the constant η, p0({right arrow over (a)}) may be derived as in Equation 16. That is, in the decomposition of Equation 15, the noise-free measurement outcome approximated on the assumption that the expected value of the Hermitian operator is the constant may be given as in Equation 16.
The measurement device performs a restoration process of Equation 16 by choosing η for the probability distribution pe({right arrow over (a)}) of the measurement outcome obtained from the actual measurement. Equation 16 includes an error according to the choice of η, and a degree of error according to η may be calculated as in Equation 17.
This error depends on the parameter ∈ and the choice of the constant η, both of which are of order 10−1. Thus, the error may be O(∈|q({right arrow over (a)})−η|) which is found in a percent level. Here, O( ) denotes a big O function.
The error of Equation 17 depends on a difference between a maximum eigenvalue and a minimum eigenvalue of the operator P. In some embodiments, the value of η may be chosen as in Equation 18.
η:=qc=½(b+(P)+b−(P)) Equation 18
In Equation 18, the maximum eigenvalue of the operator P is
and the minimum eigenvalue of the operator P is
When η of Equation 18 is used, |q({right arrow over (a)})−η| in Equation 17 is changed to δ(P), and δ(P)=(b+(P)−b−(P))/2. Therefore, an upper bound of the error in Equation 17 may be given as in Equation 19.
In some embodiments, the measurement device may perform the optimization to minimize the upper bound of Equation 19 while constraining the operator P to a predetermined condition. In some embodiments, the predetermined condition may be P≥0. Referring to Equation 11, the measurement device may find the parameter ∈ and the operator P by, for example, performing the optimization of Equation 20.
The measurement device may mitigate the detection crosstalk based on the parameter ∈ and the operator P found through the optimization (e.g., Equation 20). In some embodiments, the measurement device may output the measurement outcome in which the detection crosstalk is mitigated through the process of Equation 16 based on the parameter ∈ and the operator P.
Currently, although the error rate for the single-qubit measurements in IBMQ and Rigetti with the mitigation protocol, through the cloud-based quantum computing services is in a level of 1-10%, the error rate for two-qubit measurements rate increases to 10-30%. However, according to the above-described embodiments, the error rate is about 10% on average. As such, according to an embodiment, the error rate due to the detection crosstalk can be maintained within 10%.
While this invention has been described in connection with what is presently considered to be practical embodiments, it is to be understood that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.
Number | Date | Country | Kind |
---|---|---|---|
10-2021-0083317 | Jun 2021 | KR | national |