Patent Description:
Quantum information processing uses quantum mechanical properties to extend the capabilities of information processing. For example, security of information transfer may be enhanced, the amount of information encoded in a communication channel may be increased, and the number of operations required to perform certain computations may be reduced. Just as in conventional information processing where information is stored in one or more bits, quantum information may be stored in one or more quantum bits, known as "qubits. " A qubit may be implemented physically in any two-state quantum mechanical system, such as photon polarization, electron spin, nuclear spin, or various properties of a superconducting Josephson junction, such as charge, energy, or the direction of a current.

Additionally, quantum information may be stored and processed using "qudits," which are quantum systems with "d" number of discrete quantum states. Qubits are a specific example of a qudit with d=<NUM>. Qudits may be implemented using a physical quantum system with multiple states, such as the multiple energy levels of a quantum oscillator. The article entitled "<NPL>), and the article entitled "<NPL>), disclose examples of quantum information processing systems and methods.

The dependent claims set out preferred embodiments.

Various aspects and embodiments are described with reference to the following drawings. The drawings are not necessarily drawn to scale. For the purposes of clarity, not every component may be labeled in every drawing. In the drawings:.

Conventional quantum information processing schemes encode information in one or more two-level quantum systems (i.e., "qubits"). The state of a single qubit may be represented by the quantum state |ψ〉, which may be in any arbitrary superposition of the two quantum states, |<NUM>〉 and |<NUM>〉, e.g., |ψ〉 = α|<NUM>〉 + β|<NUM>〉, where α and β are complex numbers representing the probability amplitude of the logical qubit being in state |<NUM>〉 and |<NUM>〉, respectively.

To perform a useful quantum information process, conventional quantum information systems initialize a set of qubits to a particular quantum state, implement a set of quantum gates on the qubits, and measure the final quantum state of the qubits after performing the quantum gates. A first type of conventional quantum gate is a single-qubit gate, which transforms the quantum state of a single qubit from a first quantum state to a second quantum state. Examples of single-qubit quantum gates include the set of rotations of the qubit on a Bloch sphere. A second type of conventional quantum gate is a two-qubit gate, which transforms the quantum state of a first qubit based on the quantum state of a second qubit. Examples of two-qubit gates include the controlled NOT (CNOT) gate and the controlled phase gate. Conventional single-qubit gates and two-qubit gates unitarily evolve the quantum state of the qubits from a first quantum state to a second quantum state.

Conventional quantum information systems typically perform detections of the qubits by measuring which quantum state of a set of possible quantum states each qubit is in. This type of measurement is referred to as a projective measurement (sometimes called a projection). Examples of projective measurements include measuring the quantum state of a qubit in a particular basis to determine a detection result of either |<NUM>〉 or |<NUM>〉. Another example of a projective measurement is a measurement of the Fock basis (i.e., the photon-number eigenbasis) of a quantum oscillator. In this example, a detection result for the quantum harmonic oscillator indicates which of the infinite number of Fock states (i.e., photon-number eigenstates) |<NUM>〉|<NUM>〉 | <NUM>〉 , | <NUM>〉 , | <NUM>〉 ,. , | n〉 the quantum oscillator is in.

The inventors have recognized and appreciated that unique and powerful quantum information processes can be implemented using a more general type of quantum operation referred to as generalized quantum channels and also known as completely positive and tracepreserving (CPTP) maps. CPTP maps include not only the unitary quantum gates and projective measurements described above, but also include nonunitary quantum state evolution and generalized quantum measurements known as positive-operator valued measures (POVMs). Additionally, CPTP maps can transform a pure quantum state |ψ〉 to a mixed state, represented by the density matrix <MAT>, where i labels each pure quantum state that forms the mixed state and the coefficients pi are nonnegative and sum to one.

Previous theoretical proposals have suggested that implementing arbitrary CPTP maps for a d-dimensional system require more than one ancilla qubit and/or a circuit depth that scales poorly with the dimension, d. The inventors have, however, recognized and appreciated that any arbitrary CPTP map may be implemented for a qudit with the additional of only a single ancilla qubit and a circuit depth that is logarithmic with the dimension of the qudit. Additionally, the inventors have recognized and appreciated a technique for implementing this efficient CPTP scheme using cavity quantum electrodynamics (cQED).

Referring to <FIG>, a quantum information system <NUM> used to construct a quantum channel includes a qudit <NUM>, a single ancilla qubit <NUM>, a qubit state detector <NUM> a controller <NUM>, and a driving source <NUM>, according to some embodiments. The ancilla qubit <NUM> is coupled to the qudit <NUM> such that the state of the qudit <NUM> may affect the state of the ancilla qubit <NUM> and vice versa. For example, the ancilla qubit <NUM> and the qudit <NUM> may be dispersively coupled - meaning that a detuning between the ancilla qubit <NUM> and the qudit <NUM> (e.g., a quantum oscillator in a cavity) is much larger (e.g., an order of magnitude larger) than the coupling strength between the ancilla qubit <NUM> and the qudit <NUM>, the detuning being the frequency difference between the transition frequency of the ancilla qubit <NUM> and one or more supported modes of the cavity.

The qudit <NUM> may be any suitable d-dimensional quantum system with d quantum states. For example, the qudit <NUM> may include a quantum oscillator (e.g., harmonic or anharmonic). For example, the qudit may be physically implemented using a cavity that supports electromagnetic radiation, such as a stripline cavity or a three-dimensional conductive cavity (e.g., made from a metal such as aluminum). A subset of the photon number states of the oscillator may be the d quantum states of the qudit.

The ancilla qubit <NUM> may be any suitable two-dimensional quantum system with two quantum states. For example, the ancilla qubit may be a superconducting qubit. Examples of superconducting qubits include a superconducting charge qubit where the two quantum states relate to the charge of a superconductor, a superconducting flux qubit where the two quantum states are the direction of a current, and a superconduting phase qubit where the two quantum states are two energy eigenstates. A specific implementation of a superconducting charge qubit is a a transmission line shunted plasma oscillation ("transmon") qubit. In some embodiments, the superconducting ancilla qubit includes at least one Josephson junction.

The driving source <NUM> is coupled to the qudit <NUM> and the ancilla qubit <NUM> to enable the driving source <NUM> to control the quantum state of the qudit <NUM> and the ancilla qubit <NUM>. For example, the driving source <NUM> can implement unitary operations on the qudit <NUM> and the ancilla qubit <NUM>. For example, in embodiments where the qudit <NUM> is a quantum oscillator and the ancilla qubit <NUM> is a transmon qubit, the driving source <NUM> may create electromagnetic signals for driving the qudit <NUM> and the ancilla qubit <NUM>.

The qubit state detector <NUM> measures the state of the ancilla qubit <NUM>. In some embodiments, the measurement of the ancilla qubit <NUM> does not disturb the state of the qudit <NUM>. The qubit state detector <NUM> transmits a detection result to the controller <NUM>.

The controller <NUM> is configured to control the driving signals generated by the driving source <NUM>. In some embodiments, the controller <NUM> may determine the driving signals based on the detection result received from the qubit state detector <NUM>.

As will be described in more detail below, the qudit <NUM> and the ancilla qubit <NUM> will be operated upon by multiple unitary operations, controlled by the controller <NUM> and implemented by the driving source <NUM>. After each unitary operation, a measurement of the ancilla qubit <NUM> is performed by the qubit state detector <NUM>. A detection result from the qubit state detector <NUM> is used by the controller to determine a subsequent unitary operation to perform on the joint qudit-qubit system. In some embodiments, the quantum state of the ancilla qubit <NUM> is reset to an initial state (for example, the ground state) after each measurement and before the unitary oepration is performed on the qudit <NUM> and the ancilla qubit <NUM>.

The above-procedure is discussed in theoretical detail below, followed by an example implementation based on a particular embodiment that uses cQED devices.

Quantum channels (i.e., CPTPs) may be represented using the Kraus representation: <MAT> In Eqn. <NUM>, <IMG>(ρ) represents the CPTP map acting on a density matrix ρ, which represents the quantum state being operated upon. The CPTP map is represented using N different Kraus operators, Ki, where the dagger operator indicates the Hermitian conjugate of a Kraus operator. The Kraus operators are not necessarily unitary, HermitianD(ε), or square matrices. But the CPTP map of Eqn. <NUM> is trace preserving because <MAT>. The Kraus representation is not unique because an arbitrary different set of Kraus operators can be formed using a unitary matrix and the resulting Kraus operators represent the same CPTP map as the original Kraus operators.

To efficiently construct a CPTP map, it is convenient to work with the Kraus representation with the minimum number of Kraus operators, which is referred to as the Kraus rank of the CPTP map. Since there are at most d<NUM> linearly independent operators for a Hilbert space of dimension d, the Kraus rank is no larger than d<NUM>. A non-minimal representation of a CPTP map may be converted to the minimal Kraus representation using efficient computational techniques known in the art. For example, the Kraus representation can be converted into the Choi matrix (a d<NUM> × d<NUM> Hermitian matrix) and from there obtain the minimal Kraus representation as described in <NPL>). A second approach is to calculate an overlap matrix Cij =Tr(KiKj†) and then diagonalize it, C = V†DV. The new Kraus operators, K̃i = Σj VijKj, are the most economic representation with some of them being zero matrices if the original representation is redundant.

In some embodiments, a quantum channel may be provided in a representation other than a Kraus representation (e.g., superoperator matrix representation, Jamiolkowski/Choi matrix representation). Such alternative representations may also be converted into a minimal Kraus representation. In embodiments where the quantum channel is provided in superoperator matrix representation, the quantum channel is first converted to the Choi matrix representation. The Choi matrix is then converted to a minimal Kraus representation.

Since CPTP maps are linear in the density matrix ρ, ρ may be treated as a vector and the matrix form of the superoperator <IMG> acting on the quantum state represented by the density matrix may be written as: <MAT> or <MAT> where ρ̃ = <IMG>(ρ), is the state of the quantum system after the applying the quantum channel.

The above matrix form is a convenient representation of the quantum channel when considering the concatenation of multiple quantum channels because applying a first channel followed by a second quantum channel results in an overall channel represented by the matrix multiplication of the two superoperators representations of the two quantum channels. The matrix form of the quantum channel also allows the quantum channel to be characterized using the determinant of the matrix. For example, for Markovian channels or Kraus rank-<NUM> channels, the determinant of the matrix representation is always positive. The matrix representation, however, makes it difficult to determine whether a given matrix representing a possible quantum channel qualifies as a CPTP map. To make such a determination, the superoperator is converted to the Jamiolkowski/Choi matrix representation or Kraus representation.

Obtaining the superoperator matrix representation of a quantum channel from a Kraus representation is relatively straightforward as compared to obtaining the Kraus representation from the superoperator matrix representation. Given a channel in Kraus form, the superoperator matrix T can be obtained as follows: <MAT> where Ki are the Kraus operators (of which there are N different Kraus operators), and Ki* is the complex conjugate of the Kraus operator Ki. Obtaining the Kraus representation from the superoperator matrix T, however, uses the channel-state duality (i.e., Jamiolkowski-Choi isomorphism), from which it is known that each channel <IMG> for a system with a d-dimensional Hilbert space <IMG> corresponds (one-to-one) to a state (a density matrix) on two subsystems with a Hilbert space <IMG>⊗<IMG> as follows: <MAT> where <MAT> is the maximally entangled state of the two subsystems and τ is the Jamiolkowski matrix representation of the quantum channel. The Choi matrix M is simply a constant multiple of the Jamiolkowski matrix τ by a constant d, the dimension of the Hilbert space. The Choi matrix M and the super-operator matrix T are related as follows: <MAT>.

Being a density matrix, τ is Hermitian. Moreover, τ is semi-positive definite if and only if <IMG> is completely positive; τ is normalized if <IMG> is trace preserving. The Choi matrix M may be converted to the Kraus representation using the fact that, if M is diagonalized: <MAT> where vi are d<NUM> dimensional eigenvectors of τ. The Kraus operators are then obtained by rearranging <MAT> as d × d matrices. The number of non-zero eigenvalues λi is the Kraus rank of the corresponding quantum channel. In some embodiments, numerical calculations are performed where the eigenvalues may be truncated by, for example, setting all eigenvalues with a value less than <NUM>-<NUM> to the value <NUM>.

Having described above techniques for obtaining the minimal Kraus representation of a particular quantum channel, according to some embodiments, techniques for physically constructing a desired quantum channel are described. In some embodiments, a binary-tree scheme is used to construct any arbitrary CPTP map. The procedure to construct a CPTP map with Kraus rank N is associated with a binary tree of depth L = [log<NUM>N]. With a single ancilla qubit, the circuit depth of L = [log<NUM>N] is the lowest possible, which is what is meant herein when a construction of a quantum channel is referred to as "efficient.

Referring to <FIG>, a quantum circuit <NUM> representing the construction of an general quantum channel <IMG> includes a sequence of unitary operations <NUM>-<NUM> on the qudit <NUM> and the ancilla qubit <NUM> over time. The quantum circuit <NUM> also includes a measurement <NUM>-<NUM> after each of the unitary operations <NUM>-<NUM>. Each of the measurements <NUM>-<NUM> may be made, for example, by a single qubit state detector <NUM> at different times. While <FIG> shows a particular example with a binary tree depth L=<NUM>, it should be understood that this technique may be extended to any binary tree depth, the binary tree depth being a function of the Kraus rank of the desired quantum channel.

The quantum state of the qudit <NUM> may begin in any arbitrary state represented by a density matrix ρ. In some embodiments, the ancilla qubit <NUM> is prepared in a predetermined initial state. For example, the initial state of the ancilla qubit <NUM> may be the ground state, represented by the state |<NUM>〉. A first unitary operation <NUM> on the joint qudit-qubit system is performed after initialization of the ancilla qubit <NUM>. After the first unitary operation <NUM> is complete, a first measurement <NUM> of the ancilla qubit <NUM> is performed, resulting in a detection result <NUM>. The first detection result <NUM> may be output to the controller <NUM> (not shown in <FIG>) to determine the second unitary operation <NUM> to be performed on the joint qudit-qubit system. After the second unitary operation <NUM> is performed, a second measurement <NUM> of the ancilla qubit <NUM> is performed, resulting in a second detection result <NUM>. The second detection result <NUM> may be output to the controller <NUM>. The first detection result <NUM> (which may be stored by the controller) and the second detection result <NUM> are then used to determine the third unitary operation <NUM> to be performed on the joint qudit-qubit system. After the third unitary operation <NUM> is performed, a third measurement <NUM> of the ancilla qubit <NUM> is performed, resulting in a third detection result <NUM>. The output state of the qudit <NUM> after all the operations described is ρ̃ = <IMG>(ρ), the state of the qudit after applying the quantum channel <IMG>. In some embodiments, the classical information contained in the detection results <NUM>-<NUM> are output and may be used for additional processing.

Referring to <FIG>, a binary tree representation <NUM> used to create the quantum channel described by the quantum circuit <NUM> of <FIG>. The Kraus operators <NUM>-<NUM>, represented as Kb(l), are associated with the different leaves of the binary tree <NUM>, represented as b(L) ∈ {<NUM>,<NUM>}L. For example the Kraus operators <NUM>-<NUM>, K<NUM> and K<NUM>, are associated with a first unitary operation <NUM> (which is the same for every Kraus operator, and represented as Uø), a second unitary operation <NUM> (represented as U<NUM>) and a third unitary operation <NUM> represented as U<NUM>); the Kraus operators <NUM>-<NUM>, K<NUM> and K<NUM>, are associated with the first unitary operation <NUM>, the second unitary operation <NUM> and a fourth unitary operation <NUM> (represented as U<NUM>); the Kraus operators <NUM>-<NUM>, K<NUM> and K<NUM>, are associated with the first unitary operation <NUM>, a fifth unitary operation <NUM> (represented as U<NUM>) and a sixth unitary operation <NUM> (represented as U<NUM>); and the Kraus operators <NUM>-<NUM>, K<NUM> and K<NUM>, are associated with the first unitary operation <NUM>, the fifth unitary operation <NUM> (represented as U<NUM>) and a seventh unitary operation <NUM> (represented as U<NUM>). Thus, for a binary tree with depth L=<NUM>, there are seven different unitary operations that are determined to implement the desired quantum channel.

The series of joint qudit-qubit unitary operations is applied, as described in <FIG>. Which unitary operation is applied in the lth round (represented by Ub(l)) is based on the most recent detection result of the ancilla qubit <NUM>. The binary tree <NUM> is used to identify the unitary operation to use at each stage based on the detection result of the ancilla qubit <NUM>, the detection result always being one of two results (either a "<NUM>" or a "<NUM>"). For example, referring to <FIG> together, the first unitary operation <NUM> is not dependent on any measurement results. The second unitary operation <NUM> depends on the first detection result <NUM>. If the detection result <NUM> is a "<NUM>", then the controller <NUM> uses the second unitary operation <NUM> of <FIG> as the second unitary operation <NUM> of <FIG>. On the other hand, if the detection result <NUM> is a "<NUM>", then the controller <NUM> uses the fifth unitary operation <NUM> of <FIG> as the secondary unitary operation <NUM> of <FIG>. This same technique can be applied to subsequent unitary operations of <FIG>, such as the third unitary operation <NUM>. At each leaf of the binary tree <NUM>, the branch selected by the controller <NUM> is based on the most recent detection result. If the detection result is "<NUM>", the upper branch is selected; if the detection result is "<NUM>", the lower branch is selected.

Using the quantum circuit <NUM> and the binary tree <NUM>, any arbitrary quantum channel may be constructed and efficiently implemented.

Before describing the details of how to generate a general quantum channel according to some embodiments, a simplified example of the simple case where the tree depth L=<NUM> is described, which correspond to quantum channels with a Kraus rank of less than two. In such a situation, the quantum channel is characterized by two Kraus operators: K<NUM> and K<NUM>. In this situation, the quantum circuit of <FIG> simplifies to the following steps: (<NUM>) initialize the ancilla qubit <NUM> to the state |<NUM>〉, (<NUM>) perform a joint unitary operation U ∈ SU( <NUM>d ) , and (<NUM>) discard ("trace over") the ancilla qubit. No measurement of the ancilla qubit <NUM> is necessary because there is only one round of operation and therefore no adaptive control or feedback. Thus, in some embodiments, the ancilla qubit <NUM> is simply ignored after the joint unitary operation is performed. In some embodiments, however, the ancilla qubit <NUM> may be measured for other reasons other than creating a universal quantum channel of rank <NUM>.

The joint unitary operation U may be represented by a 2d × 2d matrix as follows: <MAT> where 〈<NUM>|U|<NUM>〉 = K<NUM> and 〈<NUM>|U|<NUM>〉 = K<NUM> are both d × d submatrices and the asterisks (*) denote other submatrices that are irrelevant in the case where U is unitary. Thus the left column of the matrix U in Eqn. <NUM> is a <NUM>d × d matrix that is an isometry, meaning the following condition is fulfilled: <MAT> The isometry condition of Eqn. <NUM> is guaranteed by the trace preserving nature of CPTP maps. When the ancilla qubit <NUM> is discarded (traced over), the quantum channel <IMG>(ρ) = K<NUM>ρK<NUM>†+K<NUM>ρK<NUM>† is realized. Thus, any quantum channel of Kraus rank <NUM> acting on a qudit <NUM> can be formed by implementing a single joint unitary operation of Eqn. <NUM> and a single ancilla qubit <NUM>.

As will be described below, if the ancilla qubit <NUM> is measured rather than ignored, "which trajectory" information is obtained in the form of a detection result. This information can be used to determine additional operations to perform when constructing a quantum channel with a Kraus rank greater than two.

Having thus described the simplified embodiment of constructing an arbitrary Kraus rank <NUM> quantum channel, a similar but more complicated technique can be used to construct an arbitrarily complex quantum channel with any Kraus rank N. In some embodiments, a quantum channel of Kraus rank N is implemented using a quantum circuit with a circuit depth of <MAT>, which is the number of joint unitary operations the quantum circuit performs in series to achieve the desired result. The quantum circuit repeats a number of "rounds" of operations, each round including: (<NUM>) initializing the ancilla qubit, (<NUM>) performing a unitary operation over the joint qudit-qubit system, the unitary operation being based on the detection result from the previous rounds (except the very first round, where the unitary is not based on a measurement result), (<NUM>) detecting the ancilla qubit <NUM>, and (<NUM>) storing the classical detection result information for use in a subsequent round. For a quantum circuit consisting of L rounds of operations with adaptive control (based on the binary detection results), there are <NUM>L -<NUM> possible unitary operations (associated with the <NUM>L -<NUM> nodes of a binary tree of depth L) and <NUM>L possible trajectories (associated with the <NUM>L leaves of the binary tree). For example, in <FIG>, the depth L=<NUM>, resulting in the seven unitary operations (<NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>, <NUM>) and eight possible trajectories leading to the leaves <NUM>-<NUM>, represented by the eight different Kraus operators of the quantum channel with Kraus rank eight.

As discussed above in connection with <FIG>, the lth unitary operation is represented by Ub(l), and is associated with the node of the binary tree, b(l) = (b<NUM>b<NUM> ··· bL) ∈ {<NUM>,<NUM>}l) with l = <NUM>,. For L=<NUM> , there is only one unitary operation for b(l) = Ø, which is Ub(l)=Ø , as determined by Eqn. <NUM> above. Generally, the unitary gate, Ub(l) is determined in a way similar to Eqn. <NUM>: <MAT>.

In embodiments where the ancilla always starts in the ground state |<NUM>〉, it is sufficient to specify the d × d submatrices 〈bl+<NUM>|Ub(l)|<NUM>〉, where |bl+<NUM>〉 is the projectively measured state of the ancilla qubit for bl+<NUM> = <NUM>,<NUM>. Each leaf of the binary tree, b(L) ∈ {<NUM>,<NUM>}Lis associated with Kraus operators labeled in binary notation, Kb(l) = Ki, where i = (b<NUM>b<NUM> ··· bL) + <NUM> and Ki>N = <NUM>, where N is the Kraus rank of the quantum channel, as illustrated by the leaves <NUM>-<NUM> of <FIG>. Each Kraus operator has an associated singular value decomposition Kraus operator Kb(L) = Wb(L)Db(L)Vb(L)†.

In some embodiments, the d × d submatrix 〈bl+<NUM>|Ub(l)|<NUM>〉 may be constructed as from the known Kraus operators of the minimal Kraus representation as follows. For each node b(l) with l - <NUM>,. , L - <NUM>, a nonnegative Hermitian matrix is determined and diagonalized as follows: <MAT> where the Vb(l) is a unitary matrix, Db(l) is a diagonal matrix with nonnegative elements, and Mb(l)is a Hermitian matrix satisfying Mb(l) = Vb(l)Db(l)Vb(l)†. For notational convenience, a matrix Pb(l) is introduced and defined as: <MAT> where sgn(<NUM>) is defined as zero such that <MAT> and Pb(l)Db(l) = Db(l)Pb(l) = Db(l). An orthogonal projection of matrix Pb(l) is defined as <MAT>, as well as the related projection Qb(l), which is defined as <MAT>. Further, an inverse matrix is defined as: <MAT> Additionally, the Moore-Penrose pseudoinverse of the matrix Mb(l) is defined as:
<MAT>
Finally, for l = <NUM> the following values are fixed: <MAT> and <MAT>.

Based on the above definitions and relations, the explicit expression for the relevant submatrices of the unitary matrices is: <MAT> with b(l+<NUM>) = (b(l),b(l+<NUM>)) for l = <NUM>,. , L - <NUM>, and
<MAT>
for l = L - <NUM>. The unitary matrix Ub(l) may therefore be completely determined based on Eqn. <NUM> and Eqn. <NUM>, using the various aforementioned definitions and Wb(l) is a unitary matrix that ensures that the isometric condition
<MAT>
is fulfilled. Because each term in Eqn. <NUM> and Eqn. <NUM> may be determined from the Kraus operators of the minimal Kraus representation, each unitary operation needed to construct a quantum circuit and binary tree similar to the examples shown in <FIG> and <FIG> may be determined from the Kraus operators of the minimal Kraus representation. It is noted that for L=<NUM>, the above equations simplify to: <MAT> which is consistent with the result for a quantum channel of Kraus rank <NUM> as discussed in connection with Eqn. <NUM> above.

<FIG> depicts a method <NUM> of operating a quantum information system that includes a qudit coupled to an ancilla qubit forming a qudit-qubit system, according to some embodiments. Method <NUM> may be applied, for example to system <NUM> shown in <FIG> and discussed above, though may also be applied to any suitable quantum system in which an ancilla qubit <NUM> is coupled to a qudit <NUM>. In some embodiments, the method <NUM> is used to implement the quantum circuit <NUM> shown in <FIG>.

At act <NUM>, the ancilla qubit <NUM> is initialized to a predetermined quantum state. In some embodiments, the ancilla qubit <NUM> is initialized to the ground state of the ancilla qubit <NUM>. The initialization of the quantum state of the ancilla qubit <NUM> may be perfomed by driving the state of the ancilla qubit <NUM> with a driving signal from the driving source <NUM> and/or performing a measurement of the ancilla <NUM> using the qubit state detector to project the ancilla qubit <NUM> into a particular state.

At act <NUM>, the driving source applies a unitary operation to the qudit-qubit system based on a previous detection result, if available. For the first round of operations, there is no previous detection result, so the unitary operation is independent of measurement results. For all other rounds of operations, one or more previous measurement results may be used to determine the unitary operation to apply to the joint qudit-qubit system. The choice of unitary operator may be determined using a binary tree structure, such as the binary tree structure <NUM> shown in <FIG>. In some embodiments, the unitary operation is implemented using one or more driving signals that act on the qudit <NUM> and the ancilla qubit <NUM> at different times such that the unitary operation is decomposed into multiple simpler unitary operations.

At act <NUM>, a detection result is generated by the qubit state detector <NUM> based on the quantum state of the ancilla qubit <NUM>. In some embodiments, the qubit state detector <NUM> may measure whether the ancilla qubit <NUM> is in the ground state or the excited state. In other embodiments, the qubit state detector <NUM> may generate a detection result by measuring the ancilla qubit <NUM> in a basis that includes superpositions of the ground state and the excited state. In some embodiments, the detection result is stored in a storage medium associated with the controller <NUM> for later use.

At act <NUM>, the controller <NUM> determines whether there are additional rounds of operations to be performed. If yes, then the method <NUM> returns to act <NUM>. If no, then the method <NUM> ends. In some embodiments, the number of rounds of operations is determined by the Kraus rank of the desired quantum channel.

The previous section describes how any arbitrary quantum channel (i.e., CPTP map) can be created using a series of unitary operations and measurements with adaptive control of the unitary operations used based on detection results from the measurements. Now, an embodiment based on a physical implementation of the qudit and ancilla qubit using cQED is described.

Referring to <FIG>, an example quantum information system <NUM> based on cQED includes storage cavity <NUM> and a transmon qubit <NUM> that are dispersively coupled together. The storage cavity <NUM> may be a stripline cavity or a three-dimensional cavity. The storage cavity <NUM> supports electromagnetic radiation, such as microwave radiation, to create a quantum oscillator. A predetermined number d of the photon number states of the quantum oscillator stored within the storage cavity <NUM> are used to implement the qudit <NUM> of <FIG>. The transmon qubit <NUM> is used as the ancilla qubit <NUM>.

Operations between the storage cavity <NUM> and the transmon qubit <NUM> may be used to perform entangling operations between the two quantum systems. These operations may implemented using driving signals generated by an electromagnetic pulse generator <NUM> controlled by a controller <NUM>.

The quantum information system includes a read-out cavity <NUM> that is also coupled to the transmon qubit <NUM>. Operations between the read-out cavity <NUM> and the transmon qubit <NUM> may map the quantum state of the transmon qubit to the state of a quantum oscillator within the read-out cavity <NUM>. These operations may be controlled by a controller <NUM>, which controls driving signals that control the operations performed on the read-out cavity and the transmon qubit <NUM>. In operation, the read-out cavity <NUM> may be operated as a fast "readout" oscillator whereas the storage cavity <NUM> may be operated as a "storage" oscillator. In some embodiments, the read-out cavity <NUM> may have a shorter decoherence time (and a lower quality factor) than the storage cavity <NUM>. When the state of the readout oscillator is detected using a cavity state detector <NUM>, the state of the storage cavity <NUM> remains undisturbed by the measurement. By transferring quantum state information from the transmon qubit <NUM> to the read-out cavity <NUM> and then detecting the quantum state of the read-out cavity using the cavity state detector <NUM>, the state of the transmon qubit <NUM> may be determined without disturbing the state of the quantum oscillator within the storage cavity <NUM>. This may be referred to as a quantum non-demolition measurement. In some embodiments, detection results from the cavity state detector <NUM> may be stored by the controller <NUM> for use in determining subsequent driving signals for controlling the transmon qubit <NUM> and the storage cavity <NUM>.

In some embodiments, electromagnetic driving pulses generated by the electromagnetic pulse generator <NUM> are used to implement unitary operations on the quantum state of the quantum oscillator stored in the storage cavity <NUM> and the quantum state of the transmon qubit <NUM>. For example, an electromagnetic signal Ω(t) may be applied to the transmon qubit <NUM> and an electromagnetic signal ε(t) may be applied to the quantum oscillator within the storage cavity <NUM>. Generally in the discussion below, application of such an electromagnetic signal or pulse may also be referred to as "driving" of the qubit or oscillator.

According to some embodiments, the joint qudit-qubit system (e.g., the joint system of the quantum oscillator and the transmon qubit <NUM>) of the quantum information system <NUM> may be described using the Hamiltonian:
<MAT>
where higher order terms are omitted. <NUM>, ωq is the qubit transition frequency between the ground state |g〉 (sometimes referred to as |<NUM>〉) and the excited state |e〉 (sometimes referred to as |<NUM>〉) of the transmon qubit <NUM>; ωc is the resonant frequency of the cavity; χ is the dispersive coupling constant between the transmon qubit <NUM> and the oscillator; â† and â are the creation and annihilation operators, respectively, for a photon within the storage cavity <NUM>. As a result of the dispersive coupling, when a photon is added to the cavity the qubit transition frequency changes by χ. Driving signals may thereby modify a particular Fock state |n〉 of the oscillator by driving the transmon qubit <NUM> (i.e., by applying an electromagnetic impulse to the transmon qubit <NUM>) at a frequency ωq + nχ. According to some embodiments, such a driving signal may modify the Fock state |n〉 by altering the phase of the state.

As illustrative yet non-limiting examples, the transmon qubit <NUM> may have a transition frequency ωq between <NUM> and <NUM>, such as between <NUM> and <NUM>, or approximately <NUM>; the quantum mechanical oscillator may have a transition frequency ωc between <NUM> and <NUM>, such as between <NUM> and <NUM>, or approximately <NUM>; the dispersive shift χ may be between <NUM> and <NUM>, such as between <NUM> and <NUM>, or such as approximately <NUM>. In some embodiments, the dispersive shift χ may be three orders of magnitude larger than the dissipation of the transmon qubit <NUM> and the storage cavity <NUM>, which allows for greater unitary control over the joint system.

<FIG> depicts an illustrative spectrum <NUM> of a transmon qubit coupled to a quantum oscillator, according to some embodiments. As discussed above, dispersive coupling between a physical qubit and a quantum mechanical oscillator causes the number states of the oscillator |n〉 to resolve to different frequencies of the transmon qubit. This configuration is sometimes referred to the "number-split regime.

<FIG> is an example of a qubit spectrum <NUM> for a qubit dispersively coupled to an resonant cavity which has an average photon number n̅ ≈ <NUM>. The horizontal axis of the figure represents the shift in the qubit transition frequency for excitations of different Fock states of the coupled resonant cavity. Put another way, the figure illustrates that the transition frequency of the transmon qubit depends on the number of photons in the cavity.

In the example spectrum <NUM> of <FIG>, the different Fock states of the oscillator |<NUM>〉, |<NUM>〉, |<NUM>, |<NUM>〉, |<NUM>〉 and |<NUM>〉 are each associated with different transition frequencies of the transmon qubit. For example, the transition frequency of the qubit where there are no photons in the cavity is defined as <NUM> of detuning (and equal to the ostensible qubit transition frequency, which as discussed above may in some embodiments be between <NUM> and <NUM>). When the cavity includes a single photon, the transition frequency of the qubit is detuned by approximately <NUM>; when the cavity includes two photons, the transition frequency of the qubit is detuned by approximately <NUM>; when the cavity includes three photons, the transition frequency of the qubit is detuned by approximately <NUM>; when the cavity includes four photons, the transition frequency of the qubit is detuned by approximately <NUM>; and when the cavity includes five photons, the transition frequency of the qubit is detuned by approximately <NUM>. This number-dependent detuning of the transition frequency can be approximated as an nχ detuning, where n is the excitation number of the cavity and χ is a detuning per photon number. For example, χ may be approximately <NUM>.

<FIG> depicts an energy level diagram <NUM> for the joint system that includes the transmon qubit <NUM> dispersively coupled to the storage cavity <NUM>, according to some embodiments. Based on this number-dependent detuning of the transition frequency of the transmon qubit <NUM>, the qubit may be addressed selectively using driving pulses with narrow spectral widths and central frequencies tuned to match the detuned transition frequencies for a particular excitation number. For example, driving the qubit at a frequency with <NUM> detuning will cause the quantum state of the cavity to change only if there is a single photon in the cavity. Thus, a driving pulse may be applied to adjust the quantum phase of a particular Fock state of the oscillator by selecting the appropriate frequency to match the targeted state. A driving pulse to implement a unitary operation may also include multiple pulses each targeting different Fock states within the same signal or separate signal. Since individual pulses may be of a different frequency, the multiple frequency components can be combined into a single pulse.

In some embodiments, the transmon qubit <NUM> may be driven independently from the storage cavity <NUM>, causing a rotation of the quantum state of the transmon qubit <NUM>. The amount of rotation of the quantum state may be dependent on the quantum state of the storage cavity <NUM> (e.g., the rotation can be photon number dependent). Such rotations induce a photon number dependent Berry phase to the quantum state of the transmon qubit <NUM>, while leaving the state of the transmon unmodified. The different phases θi are qualitatively illustrated in <FIG>, with the phase decreasing as a function of the photon number state of the quantum oscillator. This type of operation is referred to as a Selective Number-dependent Arbitrary Phase (SNAP) operation and is described in detail in <CIT>.

In some embodiments, SNAP gates may be used to implement the following entangling unitary operation: <MAT> where Yn ≡ -i|g, n〉〈e, n| + H. is the Pauli Y operator for the two-dimensional subspace associated with n photons in the storage cavity <NUM>, H. represents the Hermitian conjugate, and d is the dimensionality of the qudit, physically implemented using d energy levels of the quantum oscillator. The entangling operation Uent is associated with a quantum channel described by the Kraus operators {S<NUM>, S<NUM>}. A related entangling operation U'ent may be formed by first acting on the qudit (e.g., the storage cavity <NUM>) alone with a unitary operation V† and, after performing Uent, performing an adaptive unitary operations W<NUM> or W<NUM> on the storage cavity <NUM> alone, with the unitary operations W<NUM> or W<NUM> being dependent on the detection result from a previous measurement of the transmon qubit <NUM>. Thus: <MAT>.

The decomposition of Eqn. <NUM> is referred to as the "cosine-sine decomposition" and matches the relevant two submatrices of the desired unitary operation: <MAT> with 〈<NUM>|U|<NUM>〉 = W<NUM>S<NUM>V† and 〈<NUM>|U|<NUM>〉 = W<NUM>S<NUM>V†. Based on this, the quantum circuit similar to that of <FIG> for a general quantum channel in cQED can be determined by identifying the matrices W<NUM>, S<NUM>, W<NUM>, S<NUM>, and V for unitary operations at the different rounds U = Ub(l). In this way, potentially complicated joint unitary operations (E. , Ub(l) of <FIG>) into three simpler unitary operations: two unitary operations that act on only the oscillator and one unitary operation that acts on only the transmon qubit <NUM>.

In some embodiments, the entangling unitary operation <MAT> is determined in the following way. First, the singular value decompositions 〈<NUM>|U|<NUM>〉 = W<NUM>S<NUM>V<NUM>† and 〈<NUM>|U|<NUM>〉 = W<NUM>S<NUM>V<NUM>† are determined, with the W matrices and the S matrices set to their desired values based on Eqn. Then, it is ensured that V<NUM> = V<NUM> = V. To uniquely the decomposition, some embodiments may require that the singular values S<NUM> are arranged in descending order such that (S<NUM>)j,j ≥ (S<NUM>)j+<NUM>,J+<NUM>, while the singular values in S<NUM> are arranged in ascending order such that (S<NUM>)j,j ≤ (S<NUM>)j+<NUM>,j+<NUM>. The isometric condition
<MAT>
ensures that
<MAT>
Since both <MAT> and <MAT> are diagonal with elements in ascending order, V<NUM>†V<NUM> must be the identity, which means V<NUM> = V<NUM> = V. Thus, all the components of <MAT>, which fulfills 〈<NUM>|U|<NUM>〉 = W<NUM>S<NUM>V† and 〈<NUM>|U|<NUM>〉 = W<NUM>S<NUM>V†, are obtained.

In terms of quantum circuits, the techniques described herein for cQED systems simplifies a complex <NUM>d-dimensional unitary operation to two unitary operations acting on the qudit (e.g., quantum oscillator) alone, an entangling operation and a measurement, where the unitary operations used may be based on the detection results from the measurement.

In some embodiments, at least one non-transitory storage medium is encoded with executable instructions that, when executed by at least one processor, cause the at least one processor to carry out a method of creating a generalized quantum channel. In some embodiments, the controller <NUM> and/or the controller <NUM> may include a computer system that performs such a method. Referring to <FIG>, an example computer system <NUM> may include a processor <NUM>, a memory <NUM>, a storage device <NUM>, and input/output device(s) <NUM>. A system bus <NUM> couples the various components of the computer system <NUM> to allow the exchange of information between components. In some embodiments, the at least one non-transitory storage medium encoded with executable instructions that, when executed by at least one processor, cause the at least one processor to carry out a method of creating a generalized quantum channel may include the memory <NUM> and/or the storage device <NUM>.

The computer system <NUM> may include a variety of non-transitory computer readable media, including the memory <NUM> and the storage device <NUM>. Computer readable media can be any available media including both volatile and nonvolatile media, removable and non-removable media. Examples of computer readable media includes storage media such as RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to store the desired information and which can accessed by computer system <NUM>.

The memory <NUM> may include computer storage media in the form of volatile and/or nonvolatile memory such as read only memory (ROM) and random access memory (RAM). A basic input/output system BIOS), containing the basic routines that help to transfer information between elements within computer system <NUM>, such as during start-up, is typically stored in ROM. RAM typically contains data and/or program modules that are immediately accessible to and/or presently being operated on by the processor <NUM>. By way of example, computer system <NUM> includes software <NUM> stored in the memory <NUM> that is executable by the processor <NUM>.

The computer <NUM> may also include other removable/non-removable, volatile/nonvolatile computer storage media. By way of example only, <FIG> illustrates a storage device <NUM>. The storage device <NUM> may be a hard disk drive that reads from or writes to non-removable, nonvolatile magnetic media, a magnetic disk drive that reads from or writes to a removable, nonvolatile magnetic disk, and an optical disk drive that reads from or writes to a removable, nonvolatile optical disk such as a CD ROM or other optical media. Other removable/non-removable, volatile/nonvolatile computer storage media that can be used in the exemplary operating environment include, but are not limited to, magnetic tape cassettes, flash memory cards, digital versatile disks, digital video tape, solid state RAM, solid state ROM, and the like.

The computer system <NUM> may include a number of input/output device(s) <NUM>. For example, to facilitate operation in a networked environment, the computer system <NUM> may use the network interface to store information in network storage drives or receive information from external computer systems connected to the network. Another example of the input/output device(s) <NUM> include user interfaces that allow a user of the computer system <NUM> to input information (such as an indication of a desired quantum channel) and receive feedback. For example, the input/output device(s) <NUM> may include a keyboard, a touchscreen interface, a mouse, a microphone, a speaker, and/or a display. Another example of the input/output device(s) <NUM> is a communication interface that allows the computer system to send and receive data from to and from other devices, such as the qubit state detector <NUM> and/or the driving source <NUM>. For example, the computer system <NUM> may receive detection results via the input/output device(s) <NUM> from the qubit state detector <NUM>, store the detection results in the memory <NUM> and/or the storage device <NUM>, use the detection result to process data using the processor <NUM>, and then transmit information to control the driving signals implemented by the driving source <NUM> via the input/output device(s) <NUM>.

<FIG> depicts an example method <NUM> for creating a generalized quantum channel. The acts of method <NUM> may, by way of example, be performed by the processor <NUM> of the computer system <NUM>, which is part of the controller <NUM>.

At act <NUM>, the processor <NUM> obtains an indication of a desired quantum channel. The indication may be obtained via the input/output device(s) <NUM>. For example, a user may enter the indication using a user interface. Alternatively, the indication may be received from a different computer system via the network interface. In some embodiments, the indication may be a set of Kraus operators associated with the desired quantum channel. In another embodiment, the indication may be a super-operator matrix, a Choi matrix, or a Jamiolkowski matrix.

At act <NUM>, the processor <NUM> determines a minimal Kraus representation of the desired quantum channel. The exact procedure for achieving this is dependent on the form of the indication of the desired quantum channel. For example, if the indication of the desired quantum channel obtained in act <NUM> includes the Kraus operators of the minimal Kraus representation of the desired quantum channel, then act <NUM> may simply include verifying that the obtained Kraus operators are indeed the minimal Kraus representation. Alternatively, if the indication of the desired quantum channel is a super-operator matrix, a Choi matrix, or a Jamiolkowski matrix, the indication is converted to the minimal Kraus representation using the techniques described above.

At act <NUM>, the processor <NUM> creates a binary tree structure and generates associated unitary joint operations to be performed on the qudit-qubit system. As described above, in some embodiments, each node of the binary tree structure is associated with a respective one of the multiple unitary joint operations and each leaf of the binary tree structure is associated with a Kraus operator of the minimal Kraus representation of the desired quantum channel. In some embodiments, the unitary joint operations are generated using Eqn. <NUM> and Eqn.

At act <NUM>, the processor <NUM> converts each unitary joint operation into two unitary oscillator-only operations and one unitary qubit-only operation. In some embodiments, the unitary qubit-only operation is a SNAP operation. In some embodiments, the unitary oscillator-only operations and the unitary qubit-only operation are determined using Eqn.

At act <NUM>, the processor <NUM> determines driving signals associated with the two oscillator-only operations and the one unitary qubit-only operation. In some embodiments, the driving signals for the unitary qubit-only operation have spectral properties based on the photon number dependent transition frequencies described in <FIG>. In some embodiments, an indication of the determined driving signals is provided to the driving source <NUM>, where the driving signals are generated and directed toward the qudit <NUM> and the ancilla qubit <NUM>.

The general quantum channels (i.e., CPTP maps) described in the present application may inlcude multiple physical operations including cooling, quantum gates, measurements, and dissipative dynamics. The capability to construct an arbitrary CPTP map offers a unified approach to many aspects of quantum technology. To illustrate the wide range of impact of quantum channel construction, several example applications are described below. Embodiments are not limited to any of these applications.

A first application of constructing a generalized quantum channel is the initialization and/or stabilization of the quantum state of a qudit. Many quantum information processing tasks require working with a well-defined (often pure) initial state. One common approach is to sympathetically cool the system to the ground state by coupling to a cold bath, or optically pumping to a specific dark state, and then performing unitary operations to bring the system to a desired initial state. This can be slow if the system has a large relaxation time scale. The techniques described above, however, can actively cool the system by measurement and adaptive control. The above channel construction technique can be applied to discretely pump the qudit from an arbitrary state into the target state σ, which can be pure or mixed. The pumping time depends on the quantum gate and measurement speed, instead of the natural relaxation rate.

The quantum channel ρ ↦ εinit(ρ) =Tr(ρ)σ stabilizes the quantum state of the qudit to the target state σ. If the target state has diagonal representation σ = ∑µ λµ|ψµ〉 〈ψµ|, where λµ ≥ <NUM> and ∑µ λµ = <NUM>, one form of the Kraus operators representing the stabilizing quantum channel is <MAT>, where |i〉 are basis vectors of the Hilbert space of the qudit. Contrary to the conventional approaches discussed in the previous paragraph, this dissipative map bundles the cooling and state preparation steps and pumps an arbitrary state into the target state σ. In the case where the target state is pure, this quantum channel reduces to the "measure and rotate" procedure. Depending on the purity of the target state, entropy can be extracted from or injected into the system using the ancilla qubit. If a quantum circuit for this quantum channel is constructed using the techniques of the present application is implemented repeatedly, state stabilization can be achieved. In some embodiments, this allows a nonclassical resource state to be kept alive in a noisy quantum memory.

A second application of the generalized quantum channel construction technique described herein is in quantum error correction (QEC). In this application, multiple steady quantum states or even a subspace of steady states may be stabilized. The multiple quantum states which may be used to encode useful classical or quantum information. In some embodiments, using subspaces of steady states for QEC may include implementing a recovery map of QEC. Due to ubiquitous coupling between the qudit and the environment of the qudit, the quantum information initially stored in the qudit will unavoidably decohere as the qudit becomes entangled with the environment. Conventional QEC schemes encode quantum information in some carefully chosen logical subspaces and use syndrome measurement and conditional recovery operations to actively decouple the qudit from the environment. Despite the variety of QEC codes and recovery schemes, the operation of any QEC recovery can always be identified as a quantum channel.

For qubit-based stabilizer codes with Ns stabilizer generators, the recovery is a CPTP map with Kraus rank <NUM>Ns. In some embodiments, the ancilla qubit may be used to sequentially measure all Ns stabilizer generators to extract the syndrome, and finally perform a correction unitary operation conditioned on the syndrome pattern. Since the stabilizer generators commute with each other, their ordering does not change the syndrome. Moreover, the stabilizer measurement does not require conditioning on previous measurement outcomes, because the unitary operation at the l-th round is simply: <MAT> with Šl for the l-th stabilizer and <MAT>, which is independent of the previous measurement outcomes b(l-<NUM>). Finally, the correction unitary operation <MAT>is performed, conditioned on the syndrome b(Ns),.

In some embodiments, QEC codes that fulfill the quantum error-correction conditions associated with a set of error operations may be used. For these QEC codes, the Kraus representation of the QEC recovery map may be obtained and efficiently implemented with the construction of general quantum channels described herein. In a particular non-limiting example, a QEC code known as the binomial code uses the larger Hilbert space of higher excitations to correct excitation loss errors in bosonic systems. In order to correct up to two excitation losses, the binomial code encodes the two logical basis states as: <MAT>.

For small loss probability γ for each excitation, this encoding scheme can correct errors up to O(γ<NUM>), which includes the following four relevant processes: identity evolution (Î), losing one excitation (â), losing two excitations (â<NUM>), and back-action induced dephasing (n̂). Based on the Kraus representation of the QEC recovery (with Kraus rank <NUM>), the following set of unitary operations Ub(l) is obtained for the construction of the QEC recovery channel with an adaptive quantum circuit: <MAT> where the projection are defined as P̂i = Σk |<NUM>k + i〉〈<NUM>k + i| and P̂W = |W↑〉〈W↑| + |W↓〉〈W↓|, and the unitary operators UÔ (where Ô=â, â<NUM>, n̂) transform the error states Ô|Wσ〉back to |Wσ〉 for σ =↑, ↓. In other words: <MAT> where U⊥ is an isometry that takes the complement of the syndrome subspace to the complement of the logical subspace. In some embodiments, the first two rounds of operations, projective measurements are performed to extract the error syndrome. In the last round, a correction unitary operation is applied to restore the logical states. For example, if the measurement outcome b(<NUM>) = (<NUM>,<NUM>), there is no error and the identify operation (Î) is sufficient. If b(<NUM>) = (<NUM>,<NUM>), there is back-action induced dephasing error, which changes the coefficients of Fock states so we need to correct for that with Un̂. If b(<NUM>) = (<NUM>,<NUM>), there is a single excitation loss, which can be corrected with Uâ. If b(<NUM>) = (<NUM>,<NUM>), there are two excitation losses, which can be fully corrected with Uâ<NUM>. Repetitive application of the above QEC recovery channel can stabilize the system in the code space spanned by |W↑〉 and |W↓〉. Note that for more complicated QEC codes (e.g. GKP code [GKP_PRA_2001]) and the QEC.

In some embodiments, the QEC application may implement approximate QEC codes, which can also efficiently correct errors but only approximately fulfill the QEC criterion. For approximate QEC codes, it is challenging to analytically obtain the optimal QEC recovery map, but one can use semi-definite programming to numerically optimize the entanglement fidelity and obtain the optimal QEC recovery map. Alternatively one can use the transpose channel or quadratic recovery channels which are known to be near-optimal. All these recovery channels can be efficiently implemented using the general construction of CPTP maps described herein.

In another application of the techniques described here, the construction of generalized quantum channels can be further extended if the intermediate measurement outcomes are part of the output together with the state of the quantum system, which leads to an interesting class of quantum channels called a quantum instrument (QI). QIs enable the tracking of both the classical measurement outcome and the post-measurement state of the qudit. In some embodiments, the quantum instrument has the following CPTP map: <MAT> where |µ〉〈µ| are orthogonal projections of the measurement device with M classical outcomes, εµ are completely positive trace non-increasing maps, and <MAT> preserves the trace. Note that εµ(ρ) gives the post-measurement state associated with outcome µ.

In some embodiments, the QI is implemented as follows. (<NUM>) Find the minimum Kraus representation for εµ(ρ) (each with rank Jµ) with Kraus operators Kµ,j for j = <NUM>, <NUM>,. (<NUM>) Introduce binary labeling of these Kraus operators, Kb(L) where the binary label has length L = L<NUM> + L<NUM>, with the first <MAT> to encode µ and the remaining first L<NUM> = <MAT> bits to encode j (padding with zero operators to make a total of <NUM>L Kraus operators). (<NUM>) Use the quantum circuit with L rounds of adaptive evolution and ancilla measurement. (<NUM>) Output the final state of the quantum system as well as b(L<NUM>) that encodes µ associated with the M possible classical outcomes. In this way, the arbitrary QI described in Eq. <NUM> is constructed. In some embodiments, the QI is a used to implement complicated conditional evolution of the system. In some embodiments, the QI is a used for quantum information processing tasks that require measurement and adaptive control.

In some embodiments, which are not encompassed by the wording of the claims, the qudit is not included in the QU output. In such embodiment, the quantum channel that is constructed is effectively a positive operator valued measure (POVM), which is also referred to as a generalized quantum measurement. A POVM is a CPTP map from the quantum state of the system to the classical state of the measurement device, as represented by: <MAT> which is characterized by a set of Hermitian positive semidefinite operators <MAT> that sum to the identify operator. In some embodiments, the positive semi definite Πµ is decomposed as
<MAT>
with a set of Kraus operators {Kµ,j}j=<NUM>. Thus, in some embodiments, the quantum circuit for the quantum instrument also implements the POVM if the qudit state is removed from the QI output. In some embodiments, this reduces the binary tree construction scheme of a POVM.

In some embodiments, a POVM is used for quantum state discrimination. It is impossible for any detector to perfectly discriminate a set of non-orthogonal quantum states. An optimal detector can achieve the so-called Hellstrom bound, however, by properly designing a POVM to optimize the discrimination between the non-orthogonal state. For example, in optical communication, quadrature phase shift keying uses four coherent states with different phases to send two classical bits of information. Using the techniques for quantum channel construction described herein, an optimized POVM may be constructed.

Summarizing the above application, there are three different classifications of CPTP maps based on the output of the map, which are illustrated in <FIG> a standard quantum channel <NUM> with the quantum system (e.g., qudit <NUM>) as the output and all detection results from the measurement of the ancilla qubit are discarded; (b) a POVM <NUM> with the classical measurement outcomes <NUM>-<NUM> as the output and the quantum system discarded, with (b) not being encompassed by the wording of the claims, and (c) a QI with both the quantum system <NUM> and at least a portion of the classical measurement outcomes <NUM>-<NUM> for the output. In some embodiments, the QI keeps both the post-measurement state of the system and the outcome µ, encoded by the first L<NUM> bits of the ancilla measurement record. The remaining L<NUM> bits of the measurement record are discarded. In the QI <NUM> of <FIG>, L<NUM> = <NUM> and L<NUM> =<NUM>. In principle, all three situations can be reduced to the standard quantum channel with an expanded quantum system that includes an additional measurement device to keep track of the classical measurement outcomes. In some embodiments, however, it is more resource efficient to use a classical memory for classical measurement outcomes, so that the quantum system does not expand unnecessarily and become overly complex.

All definitions, as defined and used herein, should be understood to control over dictionary definitions, definitions, and/or ordinary meanings of the defined terms.

As used herein in the specification and in the claims, the phrase "equal" or "the same" in reference to two values (e.g., distances, widths, etc.) means that two values are the same within manufacturing tolerances. Thus, two values being equal, or the same, may mean that the two values are different from one another by ±<NUM>%.

Claim 1:
A quantum information system (<NUM>) comprising:
a physical ancilla qubit (<NUM>);
a physical qudit (<NUM>) coupled to the ancilla qubit;
a detector (<NUM>) configured to generate a detection result based on a quantum state of the ancilla qubit;
a driving source (<NUM>) coupled to the qudit and the ancilla qubit and configured to apply at least one qudit driving signal to the qudit based on the detection result and at least one qubit driving signal to the ancilla qubit based on the detection result; and
a controller (<NUM>) coupled to the driving source and the detector, wherein the controller is configured to:
obtain an indication of a desired quantum instrument channel, the desired quantum instrument channel being a completely positive and trace preserving (CPTP) map;
determine a plurality of Kraus operators corresponding to a minimal Kraus representation of the desired quantum instrument channel;
determine the at least one qudit driving signal and the at least one qubit driving signal based on the determined plurality of Kraus operators;
control the driving source to drive the qudit with the at least one qudit driving signal;
control the driving source to drive the ancilla qubit with the at least one qubit driving signal; and
subsequent to driving the qudit with the at least one qudit driving signal and driving the ancilla qubit with the at least one qubit driving signal:
- output a final state of the qudit; and
- output one or more detection results generated by the detector in response to driving the qudit with the at least one qudit driving signal and driving the ancilla qubit with the at least one qubit driving signal.