Generalized quantum channels

Abstract

According to some aspects, a quantum information system is provided that includes an ancilla qubit; a qudit coupled to the ancilla qubit, a detector configured to generate a detection result based on a quantum state of the ancilla qubit, and a driving source 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 qudit based on the detection result.

Description

FIELD

The technology described herein relates generally to quantum information systems. Specifically, the present application is directed to systems and methods for controlling a quantum mechanical system.

BACKGROUND

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=2. Qudits may be implemented using a physical quantum system with multiple states, such as the multiple energy levels of a quantum oscillator.

BRIEF SUMMARY

The present application relates generally to systems and methods for controlling a quantum mechanical system.

According to some aspects, a quantum information system is provided that includes an ancilla qubit; a qudit coupled to the ancilla qubit; a detector configured to generate a detection result based on a quantum state of the ancilla qubit; and a driving source 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 qudit based on the detection result.

According to some embodiments, the quantum information system further includes a controller coupled to the driving source and the detector, wherein the controller is configured to: receive, from the detector, a detection signal indicative of the detection result; control the driving source to drive the qudit with the at least one qudit driving signal; and control the driving source to drive the qubit with the at least one qubit driving signal.

According to some embodiments, the controller is further configured to: obtain an indication of a desired quantum channel; determine the at least one qudit driving signal based on the indication of the desired quantum channel; and determine the at least one qubit driving signal based on the indication of the desired quantum channel.

According to some embodiments, the indication of the desired quantum channel is received from a user.

According to some embodiments, the indication of the desired quantum channel comprises a plurality of Kraus operators.

According to some embodiments, the plurality of Kraus operators form a minimal Kraus representation of the desired quantum channel.

According to some embodiments, obtaining the indication of the desired quantum channel comprises determining the plurality of Kraus operators that form the minimal Kraus representation of the desired channel from one or more other operators.

According to some embodiments, the controller is further configured to determine a plurality of joint unitary operators, each of the plurality of joint unitary operators configured to act on a joint system comprising the ancilla qubit and the qudit.

According to some embodiments, the controller is further configured to create a binary tree structure associated with the plurality of unitary operators.

According to some embodiments, the controller is further configured to determine each of the plurality of joint unitary operators based on a diagonalization of the sum of a subset of the plurality of Kraus operators that are rooted from an associated node of the binary tree structure.

According to some embodiments, the controller is further configured to determine a first qudit unitary operator, a second qudit unitary operator and a qubit unitary operator associated with and based on each of the plurality of joint unitary operators.

According to some embodiments, the controller is further configured to: determine the at least one qudit driving signal based on the first qudit unitary operator and the second qudit unitary operator; and determine the at least one qubit driving signal based on the qubit unitary operator.

According to some embodiments, the desired quantum channel comprises an initialization channel for the qudit.

According to some embodiments, the desired quantum channel comprises a stabilization channel for the qudit.

According to some embodiments, the desired quantum channel comprises a quantum error correction channel for the qudit.

According to some embodiments, wherein the desired quantum channel comprises a positive operator value measurement of the qudit.

According to some embodiments, the desired quantum channel comprises a quantum instrument channel.

According to some embodiments, the ancilla qubit is a superconducting qubit.

According to some embodiments, the ancilla qubit comprises at least on Josephson junction.

According to some embodiments, the ancilla qubit comprises a transmon qubit.

According to some embodiments, the qudit comprises a quantum oscillator.

According to some embodiments, the quantum oscillator comprises electromagnetic radiation in a cavity.

According to some embodiments, the cavity is a stripline cavity.

According to some embodiments, the cavity comprises a three-dimensional metal cavity.

According to some embodiments, the driving source comprises an electromagnetic radiation source.

According to some embodiments, the electromagnetic radiation source comprises a microwave pulse generator.

According to some embodiments, the detector comprises a read-out cavity. Some aspects are directed to a method of operating a quantum information system that includes a qudit coupled to an ancilla qubit forming a qudit-qubit system, the method comprising: applying a first unitary operation to the qudit-qubit system; generating a detection result based on a quantum state of the ancilla qubit; and applying a second unitary operation to the qudit-qubit system based on the detection result.

According to some embodiments, the method further comprises initializing the ancilla qubit before applying the first unitary operation.

According to some embodiments, the method further comprises: obtaining an indication of a desired quantum channel; determining at least one qudit driving signal based on the indication of the desired quantum channel; and determining at least one qubit driving signal based on the indication of the desired quantum channel.

According to some embodiments, the indication of the desired quantum channel is received from a user.

According to some embodiments, the indication of the desired quantum channel comprises a plurality of Kraus operators.

According to some embodiments, the plurality of Kraus operators form a minimal Kraus representation of the desired quantum channel.

According to some embodiments, obtaining the indication of the desired quantum channel comprises determining the plurality of Kraus operators that form the minimal Kraus representation of the desired channel from one or more other operators.

According to some embodiments, the method further comprises determining a plurality of joint unitary operators, each of the plurality of joint unitary operators configured to act on a joint system comprising the ancilla qubit and the qudit.

According to some embodiments, the method further comprises determining a binary tree structure associated with the plurality of unitary operators.

According to some embodiments, the method further comprises determining each of the plurality of joint unitary operators based on a diagonalization of the sum of a subset of the plurality of Kraus operators that are rooted from an associated node of the binary tree structure.

According to some embodiments, the method further comprises determining a first qudit unitary operator, a second qudit unitary operator and a qubit unitary operator associated with and based on each of the plurality of joint unitary operators. According to some embodiments, the method further comprises: determining the at least one qudit driving signal based on the first qudit unitary operator and the second qudit unitary operator; and determining the at least one qubit driving signal based on the qubit unitary operator.

According to some embodiments, the desired quantum channel comprises an initialization channel for the qudit.

According to some embodiments, the desired quantum channel comprises a stabilization channel for the qudit.

According to some embodiments, the desired quantum channel comprises a quantum error correction channel for the qudit.

According to some embodiments, the desired quantum channel comprises a positive operator value measurement of the qudit.

According to some embodiments, the desired quantum channel comprises a quantum instrument channel.

According to some embodiments, the ancilla qubit is a superconducting qubit.

According to some embodiments, wherein the ancilla qubit comprises at least on Josephson junction.

According to some embodiments, the ancilla qubit comprises a transmon qubit.

According to some embodiments, the qudit comprises a quantum oscillator.

According to some embodiments, the quantum oscillator comprises electromagnetic radiation in a cavity.

According to some embodiments, the cavity is a stripline cavity.

According to some embodiments, the cavity comprises a three-dimensional metal cavity.

According to some embodiments, the driving signal comprises an electromagnetic radiation signal.

According to some embodiments, the electromagnetic radiation signal comprises a microwave radiation signal.

According to some embodiments, the detector comprises a read-out cavity.

Some aspects are directed to 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, wherein the method comprises: obtaining a plurality of Kraus operators associated with a desired quantum channel; generating a plurality of unitary joint operations associated with a joint qudit-qubit system based on the Kraus operators; determining a first unitary joint operation of a plurality of unitary joint operations to perform on the joint qudit-qubit system; and determining two unitary qudit-only operations and one unitary qubit-only operation based on the first unitary joint operation.

According to some embodiments, obtaining the plurality of Kraus operators comprises: receiving an indication of a desired quantum channel; and determining a minimal Kraus representation of the quantum channel based on the indication.

According to some embodiments, the indication of the desired channel comprises a plurality of non-minimal Kraus operators.

According to some embodiments, the indication of the desired channel comprises a super-operator matrix, a Choi matrix, or a Jamiolkowski matrix.

According to some embodiments, the method further comprises generating a binary tree structure associated with the plurality of unitary joint operations.

According to some embodiments, each node of the binary tree structure is associated with a respective unitary joint operation of the plurality of unitary joint operations.

According to some embodiments, each leaf of the binary tree structure is associated with a Kraus operator of the desired quantum channel.

According to some embodiments, the unitary qubit-only operation is a selective number-dependent arbitrary phase operation.

The foregoing is a non-limiting summary of the invention, which is defined by the appended claims.

BRIEF DESCRIPTION OF DRAWINGS

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:

FIG. 1 is block diagram of a quantum information system, according to some embodiments.

FIG. 2 depicts a quantum circuit for constructing an arbitrary quantum channel, according to some embodiments.

FIG. 3 depicts a binary tree structure used to construct an arbitrary quantum channel, according to some embodiments.

FIG. 4 is a flowchart of a method of operating a quantum information system, according to some embodiments.

FIG. 5 is a block diagram of a quantum information system based on cavity quantum electrodynamics, according to some embodiments.

FIG. 6 depicts an example spectrum of a transmon qubit coupled to a quantum oscillator of a storage cavity, according to some embodiments.

FIG. 7 depicts an energy level diagram of a transmon qubit coupled to a quantum oscillator of a storage cavity, according to some embodiments.

FIG. 8 is a block diagram of a computer system, according to some embodiments.

FIG. 9 is a flow chart of a method of operating a quantum information system, according to some embodiments.

FIG. 10 is a block diagram of different types of quantum channels classified by their outputs.

DETAILED DESCRIPTION

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, |0 and |1, e.g., |ψ=α|0+β|1, where α and β are complex numbers representing the probability amplitude of the logical qubit being in state |0 and |1, 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 |0 or |1. 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) 0, |1, |2. . . . , |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 trace-preserving (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

$\rho =\sum _i{p}_i{\psi }_i\rangle \langle {\psi }_i,$ where i labels each pure quantum state that forms the mixed state and the coefficients p_i 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. 1, a quantum information system 101 used to construct a quantum channel includes a qudit 110, a single ancilla qubit 120, a qubit state detector 130 a controller 140, and a driving source 150, according to some embodiments. The ancilla qubit 120 is coupled to the qudit 110 such that the state of the qudit 110 may affect the state of the ancilla qubit 120 and vice versa. For example, the ancilla qubit 120 and the qudit 110 may be dispersively coupled—meaning that a detuning between the ancilla qubit 120 and the qudit 110 (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 120 and the qudit 110, the detuning being the frequency difference between the transition frequency of the ancilla qubit 120 and one or more supported modes of the cavity.

The qudit 110 may be any suitable d-dimensional quantum system with d quantum states. For example, the qudit 110 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 120 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 transmission line shunted plasma oscillation (“transmon”) qubit. In some embodiments, the superconducting ancilla qubit includes at least one Josephson junction.

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

The qubit state detector 130 measures the state of the ancilla qubit 120. In some embodiments, the measurement of the ancilla qubit 120 does not disturb the state of the qudit 110. The qubit state detector 130 transmits a detection result to the controller 140.

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

As will be described in more detail below, the qudit 110 and the ancilla qubit 120 will be operated upon by multiple unitary operations, controlled by the controller 140 and implemented by the driving source 150. After each unitary operation, a measurement of the ancilla qubit 120 is performed by the qubit state detector 130. A detection result from the qubit state detector 130 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 120 is reset to an initial state (for example, the ground state) after each measurement and before the unitary operation is performed on the qudit 110 and the ancilla qubit 120.

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

Obtaining the Minimal Kraus Representation of a Quantum Channel

Quantum channels (i.e., CPTPs) may be represented using the Kraus representation:

$\begin{array}{cc}\left(\rho \right)=\sum _{i=1}^N{K}_i\rho K_i^{†}.& \left(\mathrm{Eqn}.1\right)\end{array}$ In Eqn. 1, (ρ) 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, K_i, where the dagger operator indicates the Hermitian conjugate of a Kraus operator. The Kraus operators are not necessarily unitary, Hermitian, or square matrices. But the CPTP map of Eqn. 1 is trace preserving because

$\sum _{i=1}^N{K}_i^{†}{K}_i=$ 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 according to some embodiments, 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² linearly independent operators for a Hilbert space of dimension d, the Kraus rank is no larger than d². 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²×d² Hermitian matrix) and from there obtain the minimal Kraus representation as described in M. D. Choi, Linear Algebra Appl., 10, 285 (1975), which is incorporated herein by reference in its entirety and at least for its discussion of a technique for determining a minimal Kraus representation. (In a case that any terminology used herein conflicts with the usage of that terminology in Choi, the terminology should be afforded a meaning most consistent with how a person of ordinary skill would understand its usage herein.) A second approach is to calculate an overlap matrix C_ij=Tr(K_iK_j^†) and then diagonalize it, C=V^†DV. The new Kraus operators, {tilde over (K)}_i=Σ_jV_ijK_j, 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 acting on the quantum state represented by the density matrix may be written as:

$\begin{array}{cc}{\tilde{\rho }}_{\mathrm{ij}}=\sum _{m,n}{T}_{\mathrm{ij},\mathrm{mn}}{\rho }_{\mathrm{mn}}& \left(\mathrm{Eqn}.2a\right)\\ \mathrm{or}& \\ \overrightarrow{\tilde{\rho }}=T·\overrightarrow{\rho }& \left(\mathrm{Eqn}.2b\right)\end{array}$ where {tilde over (ρ)}=(ρ), 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-2 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:

$\begin{array}{cc}T=\sum _i^N{K}_i\otimes K_i^{*}& \left(\mathrm{Eqn}.3\right)\end{array}$ where K_i are the Kraus operators (of which there are N different Kraus operators), and K_i* is the complex conjugate of the Kraus operator K_i. 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 for a system with a d-dimensional Hilbert space corresponds (one-to-one) to a state (a density matrix) on two subsystems with a Hilbert space ⊗ as follows:τ=(⊗I)(|Ω)(Ω|)  (Eqn. 4)where

$\Omega \rangle =\frac{1}{\sqrt{d}}\sum _ii\rangle \otimes i\rangle$ 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:T_ij,mn=M_im,.  (Eqn. 5)

Being a density matrix, τ is Hermitian. Moreover, τ is semi-positive definite if and only if is completely positive; is normalized if T is trace preserving. The Choi matrix M may be converted to the Kraus representation using the fact that, if M is diagonalized:

$\begin{array}{cc}M=\sum _i{\lambda }_i{\upsilon }_i{\upsilon }_i^{†},& \left(\mathrm{Eqn}.6\right)\end{array}$ where v_i are d² dimensional eigenvectors of τ. The Kraus operators are then obtained by rearranging √{square root over (λ_i)}v_i 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 10⁻¹⁰ to the value 0.

Theory of Universal Construction of Quantum Channels

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₂N┐. With a single ancilla qubit, the circuit depth of L=┌log₂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. 2, a quantum circuit 200 representing the construction of an general quantum channel includes a sequence of unitary operations 251-253 on the qudit 110 and the ancilla qubit 120 over time. The quantum circuit 200 also includes a measurement 231-233 after each of the unitary operations 251-253. Each of the measurements 231-233 may be made, for example, by a single qubit state detector 130 at different times. While FIG. 2 shows a particular example with a binary tree depth L=3, 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 110 may begin in any arbitrary state represented by a density matrix ρ. In some embodiments, the ancilla qubit 120 is prepared in a predetermined initial state. For example, the initial state of the ancilla qubit 120 may be the ground state, represented by the state |0>. A first unitary operation 251 on the joint qudit-qubit system is performed after initialization of the ancilla qubit 120. After the first unitary operation 251 is complete, a first measurement 231 of the ancilla qubit 120 is performed, resulting in a detection result 241. The first detection result 241 may be output to the controller 140 (not shown in FIG. 2) to determine the second unitary operation 252 to be performed on the joint qudit-qubit system. After the second unitary operation 252 is performed, a second measurement 232 of the ancilla qubit 120 is performed, resulting in a second detection result 242. The second detection result 243 may be output to the controller 140. The first detection result 241 (which may be stored by the controller) and the second detection result 242 are then used to determine the third unitary operation 253 to be performed on the joint qudit-qubit system. After the third unitary operation 253 is performed, a third measurement 233 of the ancilla qubit 120 is performed, resulting in a third detection result 243. The output state of the qudit 120 after all the operations described is {tilde over (ρ)}=(ρ), the state of the qudit after applying the quantum channel . In some embodiments, the classical information contained in the detection results 241-243 are output and may be used for additional processing.

Referring to FIG. 3, a binary tree representation 300 used to create the quantum channel described by the quantum circuit 200 of FIG. 2. The Kraus operators 351-358, represented as K_b(l), are associated with the different leaves of the binary tree 300, represented as b^(L)∈ {0,1}^L. For example the Kraus operators 351-352, K₀₀₀ and K₀₀₁, are associated with a first unitary operation 251 (which is the same for every Kraus operator, and represented as U_Ø), a second unitary operation 320 (represented as U₀) and a third unitary operation 330 represented as U₀₀); the Kraus operators 353-354, K₀₁₀ and K₀₁₁, are associated with the first unitary operation 251, the second unitary operation 320 and a fourth unitary operation 331 (represented as U₀₁); the Kraus operators 355-356, k₁₀₀ and K₁₀₁, are associated with the first unitary operation 251, a fifth unitary operation 321 (represented as U₁) and a sixth unitary operation 332 (represented as U₁₀); and the Kraus operators 357-358, K₁₁₀ and K₁₁₁, are associated with the first unitary operation 251, the fifth unitary operation 321 (represented as U₁) and a seventh unitary operation 333 (represented as U₁₁). Thus, for a binary tree with depth L=3, 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. 2. Which unitary operation is applied in the lth round (represented by U_b(l)) is based on the most recent detection result of the ancilla qubit 120. The binary tree 300 is used to identify the unitary operation to use at each stage based on the detection result of the ancilla qubit 120, the detection result always being one of two results (either a “0” or a “1”). For example, referring to FIGS. 2-3 together, the first unitary operation 251 is not dependent on any measurement results. The second unitary operation 252 depends on the first detection result 241. If the detection result 252 is a “0”, then the controller 140 uses the second unitary operation 320 of FIG. 3 as the second unitary operation 252 of FIG. 2. On the other hand, if the detection result 252 is a “1”, then the controller 140 uses the fifth unitary operation 321 of FIG. 3 as the secondary unitary operation 252 of FIG. 2. This same technique can be applied to subsequent unitary operations of FIG. 2, such as the third unitary operation 252. At each leaf of the binary tree 300, the branch selected by the controller 140 is based on the most recent detection result. If the detection result is “0”, the upper branch is selected; if the detection result is “1”, the lower branch is selected.

Using the quantum circuit 200 and the binary tree 300, 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=1 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_o and K₁. In this situation, the quantum circuit of FIG. 2 simplifies to the following steps: (1) initialize the ancilla qubit 120 to the state |0>, (2) perform a joint unitary operation U∈SU(2d), and (3) discard (“trace over”) the ancilla qubit. No measurement of the ancilla qubit 120 is necessary because there is only one round of operation and therefore no adaptive control or feedback. Thus, in some embodiments, the ancilla qubit 120 is simply ignored after the joint unitary operation is performed. In some embodiments, however, the ancilla qubit 120 may be measured for other reasons other than creating a universal quantum channel of rank 2.

The joint unitary operation U may be represented by a 2d×2d matrix as follows:

$\begin{array}{cc}U=\left(\begin{array}{cc}\langle 0U0\rangle & *\\ \langle 1U0\rangle & *\end{array}\right),& \left(\mathrm{Eqn}.7\right)\end{array}$ where 0|U|0=K₀ and 1|U|0=K₁ 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. 7 is a 2d×d matrix that is an isometry, meaning the following condition is fulfilled:

$\begin{array}{cc}\sum _{b=0,1}{\left(\langle bU0\rangle \right)}^{†}\langle bU0\rangle ={d\times d}_{}.& \left(\mathrm{Eqn}.8\right)\end{array}$ The isometry condition of Eqn. 7 is guaranteed by the trace preserving nature of CPTP maps. When the ancilla qubit 120 is discarded (traced over), the quantum channel _U(ρ)=K₀ρK₀^†+K₁ρK₁^† is realized. Thus, any quantum channel of Kraus rank 2 acting on a qudit 110 can be formed by implementing a single joint unitary operation of Eqn. 7 and a single ancilla qubit 120.

As will be described below, if the ancilla qubit 120 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 2 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 L=┌log₂N┐, 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: (1) initializing the ancilla qubit, (2) 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), (3) detecting the ancilla qubit 120, and (4) 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 2^L−1 possible unitary operations (associated with the 2^L−1 nodes of a binary tree of depth L) and 2^L possible trajectories (associated with the 2^L leaves of the binary tree). For example, in FIG. 3, the depth L=3, resulting in the seven unitary operations (251, 320, 321, 330, 331, 332, 333) and eight possible trajectories leading to the leaves 351-358, represented by the eight different Kraus operators of the quantum channel with Kraus rank eight.

As discussed above in connection with FIG. 2, the lth unitary operation is represented by U_b(l), and is associated with the node of the binary tree, b^(l)=(b₁ b₂ . . . b_L)∈{0,1}^l) with l=0, . . . , L−1. For L=1, there is only one unitary operation for b^(l)=Ø, which is U_b(l)Ø, as determined by Eqn. 7 above. Generally, the unitary gate, U_b(l) is determined in a way similar to Eqn. 7:

$\begin{array}{cc}{U}_{b^{\left(l\right)}}=\left(\begin{array}{cc}\langle 0U_{b^{\left(l\right)}}0\rangle & *\\ \langle 1U_{b^{\left(l\right)}}0\rangle & *\end{array}\right).& \left(\mathrm{Eqn}.9\right)\end{array}$

In embodiments where the ancilla always starts in the ground state 10), it is sufficient to specify the d×d submatrices b_l+1|U_b(l)|0, where |b_l+1 is the projectively measured state of the ancilla qubit for b_l+1=0,1. Each leaf of the binary tree, b^(L)∈{0,1}^L is associated with Kraus operators labeled in binary notation, K_b(l)=K_i, where i=(b₁b₂ . . . b_L)+1 and K_i>N=0, where N is the Kraus rank of the quantum channel, as illustrated by the leaves 351-358 of FIG. 3. Each Kraus operator has an associated singular value decomposition Kraus operator K_b(L)=W_b(L)D_b(L)V_b(L)^†.

In some embodiments, the d×d submatrix <b_i+1|U_b(l)|0> may be constructed as from the known Kraus operators of the minimal Kraus representation as follows. For each node b^(l) with l−1, . . . , L−1, a nonnegative Hermitian matrix is determined and diagonalized as follows:

$\begin{array}{cc}\sum _{b_{l+1},\dots ,b_L}{K}_{b^{\left(L\right)}}^{†}{K}_{b^{\left(L\right)}}=V_{b^{\left(l\right)}}{D}_{b^{\left(l\right)}}^2{V}_{b^{\left(l\right)}}^{†}\equiv M_{b^{\left(l\right)}}^2,& \left(\mathrm{Eqn}.10\right)\end{array}$ where the V_b(l) is a unitary matrix, D_b(l) is a diagonal matrix with nonnegative elements, and M_b(l) is a Hermitian matrix satisfying M_b(l)=V_b(l)D_b(l)V_b(l)^†. For notational convenience, a matrix P_b(l) is introduced and defined as:(P_b(l))_j,k=sgn[(D_b(l))_j,k],  (Eqn. 11)where sgn(0) is defined as zero such that P_b(l)²=P_b(l) and P_b(l)D_b(l)=D_b(l)P_b(l)=D_b(l). An orthogonal projection of matrix P_b(l) is defined as P_b(l)^⊥=−P_b(l), as well as the related projection Q_b(l), which is defined asQ_b(l)≡V_b(l)P_b(l)^⊥V_b(l)^⊥  (Eqn. 12)Further, an inverse matrix is defined as:

$\begin{array}{cc}{\left(D_{b^{\left(l\right)}}^{-1}\right)}_{j,k}=\{\begin{array}{cc}1/{\left(D_{b^{\left(l\right)}}\right)}_{j,k}& \mathrm{if}{D}_{b^{\left(l\right)}}\ne 0,\\ 0& \mathrm{otherwise}.\end{array}& \left(\mathrm{Eqn}.13\right)\end{array}$ Additionally, the Moore-Penrose pseudoinverse of the matrix M_b(l) is defined as:M_b(l)⁺=V_b(l)D_b(l)⁻¹V_b(l)^†.  (Eqn. 14)Finally, for l=0 the following values are fixed: V_b(o)=D_b(o)=D_b(o)⁻¹=P_b(o)= and P_b(o)^⊥=0.

Based on the above definitions and relations, the explicit expression for the relevant submatrices of the unitary matrices is:

$\begin{array}{cc}\langle b_{l+1}U_{b^{\left(l\right)}}0\rangle =M_{b^{\left(l+1\right)}}{M}_{b^{\left(l\right)}}^{+}+\frac{1}{\sqrt{2}}{Q}_{b^{\left(l\right)}},& \left(\mathrm{Eqn}.15\right)\end{array}$ with b⁽ⁱ⁺¹⁾=(b^(l),b_(i+1)) for l=0, . . . ,L−2, and

$\begin{array}{cc}\langle b_{l+1}U_{b^{\left(l\right)}}0\rangle =K_{b^{\left(l+1\right)}}{M}_{b^{\left(k\right)}}^{+}+\frac{1}{\sqrt{2}}{W}_{b^{\left(l+1\right)}}{V}_{b^{\left(l+1\right)}}^{†}{Q}_{b^{\left(l\right)}}& \left(\mathrm{Eqn}.16\right)\end{array}$ for l=L−1. The unitary matrix U_b(l) may therefore be completely determined based on Eqn. 15 and Eqn. 16, using the various aforementioned definitions and W_b(l) is a unitary matrix that ensures that the isometric condition Σ_bl+1=0,1(b_l+1|U_b(l)|0)^†(b_l+1|U_b(l)|0=_d×d is fulfilled. Because each term in Eqn. 15 and Eqn. 16 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 FIGS. 2 and 3 may be determined from the Kraus operators of the minimal Kraus representation. It is noted that for L=1, the above equations simplify to:

$\begin{array}{cc}\begin{array}{c}\langle b_1U_{b^{\left(0\right)}}0\rangle =K_{b^{\left(1\right)}}\\ =\{\begin{array}{cc}{K}_1& \mathrm{for}{b}_1=0,\\ K_2& \mathrm{for}{b}_1=1,\end{array}\end{array}\hspace{1em}& \left(\mathrm{Eqn}.17\right)\end{array}$ which is consistent with the result for a quantum channel of Kraus rank 2 as discussed in connection with Eqn. 7 above.

FIG. 4 depicts a method 400 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 400 may be applied, for example to system 100 shown in FIG. 1 and discussed above, though may also be applied to any suitable quantum system in which an ancilla qubit 120 is coupled to a qudit 110. In some embodiments, the method 400 is used to implement the quantum circuit 200 shown in FIG. 2.

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

At act 404, 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 300 shown in FIG. 3. In some embodiments, the unitary operation is implemented using one or more driving signals that act on the qudit 110 and the ancilla qubit 120 at different times such that the unitary operation is decomposed into multiple simpler unitary operations.

At act 406, a detection result is generated by the qubit state detector 130 based on the quantum state of the ancilla qubit 120. In some embodiments, the qubit state detector 130 may measure whether the ancilla qubit 120 is in the ground state or the excited state. In other embodiments, the qubit state detector 130 may generate a detection result by measuring the ancilla qubit 120 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 140 for later use.

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

Embodiments Based On Circuit QED

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. 5, an example quantum information system 500 based on cQED includes storage cavity 510 and a transmon qubit 520 that are dispersively coupled together. The storage cavity 510 may be a stripline cavity or a three-dimensional cavity. The storage cavity 510 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 510 are used to implement the qudit 110 of FIG. 1. The transmon qubit 520 is used as the ancilla qubit 120.

Operations between the storage cavity 510 and the transmon qubit 520 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 550 controlled by a controller 540.

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

In some embodiments, electromagnetic driving pulses generated by the electromagnetic pulse generator 550 are used to implement unitary operations on the quantum state of the quantum oscillator stored in the storage cavity 510 and the quantum state of the transmon qubit 520. For example, an electromagnetic signal Ω(t) may be applied to the transmon qubit 520 and an electromagnetic signal Ω(t) may be applied to the quantum oscillator within the storage cavity 510. 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 520) of the quantum information system 500 may be described using the Hamiltonian:H/ℏ=ω_câ^†{circumflex over (α)}+ω_q|ee|++χâ^†â|ee|,  (Eqn. 18)where higher order terms are omitted. In Eqn. 18, ω_q is the qubit transition frequency between the ground state |g (sometimes referred to as |0) and the excited state |e (sometimes referred to as |1) of the transmon qubit 520; ω_c is the resonant frequency of the cavity; χ is the dispersive coupling constant between the transmon qubit 520 and the oscillator; â^† and â are the creation and annihilation operators, respectively, for a photon within the storage cavity 510. 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 520 (i.e., by applying an electromagnetic impulse to the transmon qubit 520) 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 520 may have a transition frequency ω_q between 5 GHz and 10 GHz, such as between 7 GHZ and 8 GHZ, or approximately 7.6 GHz; the quantum mechanical oscillator may have a transition frequency ω_c between 6 GHz and 11 GHz, such as between 8 GHz and 9 GHz, or approximately 8.2 GHz; the dispersive shift χ may be between 1 MHz and 10 MHz, such as between 4 MHz and 9 MHz, or such as approximately 8.2 MHz. In some embodiments, the dispersive shift χ may be three orders of magnitude larger than the dissipation of the transmon qubit 520 and the storage cavity 510, which allows for greater unitary control over the joint system.

FIG. 5 depicts an illustrative spectrum 500 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. 6 is an example of a qubit spectrum 600 for a qubit dispersively coupled to an resonant cavity which has an average photon number n≈2. 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 600 of FIG. 5, the different Fock states of the oscillator |0,|1,|2,|3,|4 and |5 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 0 MHz of detuning (and equal to the ostensible qubit transition frequency, which as discussed above may in some embodiments be between 5 GHz and 10 GHz). When the cavity includes a single photon, the transition frequency of the qubit is detuned by approximately 10 MHz; when the cavity includes two photons, the transition frequency of the qubit is detuned by approximately 17 MHz; when the cavity includes three photons, the transition frequency of the qubit is detuned by approximately 26 MHz; when the cavity includes four photons, the transition frequency of the qubit is detuned by approximately 34 MHz; and when the cavity includes five photons, the transition frequency of the qubit is detuned by approximately 43 MHz. 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 8.3 MHz.

FIG. 7 depicts an energy level diagram 700 for the joint system that includes the transmon qubit 520 dispersively coupled to the storage cavity 510, according to some embodiments. Based on this number-dependent detuning of the transition frequency of the transmon qubit 520, 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 10 MHz 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 520 may be driven independently from the storage cavity 510, causing a rotation of the quantum state of the transmon qubit 520. The amount of rotation of the quantum state may be dependent on the quantum state of the storage cavity 510 (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 520, while leaving the state of the transmon unmodified. The different phases θ_i are qualitatively illustrated in FIG. 6, 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 U.S. patent application Ser. No. 15/552,998, entitled “TECHNIQUES OF OSCILLATOR CONTROL FOR QUANTUM INFORMATION PROCESSING AND RELATED SYSTEMS AND METHODS,” and filed Aug. 23, 2017, which is incorporated herein by reference in its entirety and at least for its discussion of implementing SNAP gates. (In a case that any terminology used herein conflicts with the usage of that terminology in U.S. patent application Ser. No. 15/552,998, the terminology should be afforded a meaning most consistent with how a person of ordinary skill would understand its usage herein.)

In some embodiments, SNAP gates may be used to implement the following entangling unitary operation:

$\begin{array}{cc}\begin{array}{c}{U}_{\mathrm{cnt}}\left({\theta }_i\right)=\left(\begin{array}{cc}{S}_0& -S_1\\ S_1& S_0\end{array}\right)\\ =\left(\begin{array}{cccccc}\cos \frac{{\theta }_1}{2}& & & -\sin \frac{{\theta }_1}{2}& & \\ & \ddots & & & \ddots & \\ & & \cos \frac{{\theta }_d}{2}& & & -\sin \frac{{\theta }_d}{2}\\ \sin \frac{{\theta }_1}{2}& & & \cos \frac{{\theta }_1}{2}& & \\ & \ddots & & & \ddots & \\ & & \sin \frac{{\theta }_d}{2}& & & \cos \frac{{\theta }_d}{2}\end{array}\right)\\ =\prod _{n=0}^{d-1}\exp \left(-{\mathrm{iY}}_n{\theta }_n/2\right),\end{array}\hspace{1em}& \left(\mathrm{Eqn}.19\right)\end{array}$ where Y_n≡−i|g,ne,n|+H.c. is the Pauli Y operator for the two-dimensional subspace associated with n photons in the storage cavity 510, H.c. 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 U_ent is associated with a quantum channel described by the Kraus operators {S₀, S₁}. A related entangling operation U′_ent may be formed by first acting on the qudit (e.g., the storage cavity 510) alone with a unitary operation V^† and, after performing U_ent, performing an adaptive unitary operations W₀ or W₁ on the storage cavity 510 alone, with the unitary operations W₀ or W₁ being dependent on the detection result from a previous measurement of the transmon qubit 520. Thus:

$\begin{array}{cc}\begin{array}{c}{U}_{\mathrm{ent}}^{\prime }=\left(\begin{array}{cc}{W}_0& 0\\ 0& W_1\end{array}\right)\left(\begin{array}{cc}{S}_0& -S_1\\ S_1& S_0\end{array}\right)\left(\begin{array}{cc}{V}^{†}& 0\\ 0& V^{†}\end{array}\right)\\ =\left(\begin{array}{cc}{W}_0{S}_0{V}^{†}& *\\ W_1{S}_1{V}^{†}& *\end{array}\right).\end{array}\hspace{1em}& \left(\mathrm{Eqn}.20\right)\end{array}$

The decomposition of Eqn. 20 is referred to as the “cosine-sine decomposition” and matches the relevant two submatrices of the desired unitary operation:

$\begin{array}{cc}U=\left(\begin{array}{cc}\langle 0U0\rangle & *\\ \langle 1U0\rangle & *\end{array}\right),& \left(\mathrm{Eqn}.21\right)\end{array}$ with 0|U|0=W₀S₀V^† and 1|U|0=W₁S₁V^†. Based on this, the quantum circuit similar to that of FIG. 2 for a general quantum channel in cQED can be determined by identifying the matrices W₀, S₀, W₁, S₁, and V for unitary operations at the different rounds U=U_b(l). In this way, potentially complicated joint unitary operations (E.g., U_b(l) of FIG. 2) 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 520.

In some embodiments, the entangling unitary operation U′_ent is determined in the following way. First, the singular value decompositions (0|U|0)=W₀S₀V₀^† and (1|U|0)=W₁S₁V₁^† are determined, with the W matrices and the S matrices set to their desired values based on Eqn. 20. Then, it is ensured that V₀=V₁=V. To uniquely the decomposition, some embodiments may require that the singular values S₀ are arranged in descending order such that (S₀)_j,j≥(S₀)_j+1,j+1, while the singular values in S₁ are arranged in ascending order such that (S_i)_j,j≤(S₁)_j+1,j+1. The isometric condition Σ_b=0,1(b|U|0)†b|U|0=_d×d ensures that V₀^†V₁S₁²V₁^†V₀=_d×d−S₀². Since both S₁² and _d×d−S₀² are diagonal with elements in ascending order, V₁^†V₀ must be the identity, which means V₀=V₁=V. Thus, all the components of U_ent′, which fulfills (0|U|0)=W₀S₀V^† and 1|U|0=W₁S₁V^†, are obtained.

In terms of quantum circuits, the techniques described herein for cQED systems simplifies a complex 2d-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.

Computer and Software Aspects

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 140 and/or the controller 540 may include a computer system that performs such a method. Referring to FIG. 8, an example computer system 800 may include a processor 810, a memory 820, a storage device 830, and input/output device(s) 840. A system bus 801 couples the various components of the computer system 800 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 820 and/or the storage device 830.

The computer system 800 may include a variety of non-transitory computer readable media, including the memory 820 and the storage device 830. 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 800.

The memory 820 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 800, 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 810. By way of example, computer system 800 includes software 822 stored in the memory 820 that is executable by the processor 810.

The computer 800 may also include other removable/non-removable, volatile/nonvolatile computer storage media. By way of example only, FIG. 8 illustrates a storage device 830. The storage device 830 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 800 may include a number of input/output device(s) 840. For example, to facilitate operation in a networked environment, the computer system 800 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) 840 include user interfaces that allow a user of the computer system 800 to input information (such as an indication of a desired quantum channel) and receive feedback. For example, the input/output device(s) 840 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) 840 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 130 and/or the driving source 150. For example, the computer system 800 may receive detection results via the input/output device(s) 840 from the qubit state detector 130, store the detection results in the memory 820 and/or the storage device 830, use the detection result to process data using the processor 810, and then transmit information to control the driving signals implemented by the driving source 150 via the input/output device(s) 840.

FIG. 9 depicts an example method 900 for creating a generalized quantum channel. The acts of method 900 may, by way of example, be performed by the processor 810 of the computer system 800, which is part of the controller 140.

At act 902, the processor 810 obtains an indication of a desired quantum channel. The indication may be obtained via the input/output device(s) 840. 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 904, the processor 810 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 902 includes the Kraus operators of the minimal Kraus representation of the desired quantum channel, then act 902 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 906, the processor 810 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. 15 and Eqn. 16.

At act 908, the processor 908 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. 20.

At act 910, the processor 810 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. 7. In some embodiments, an indication of the determined driving signals is provided to the driving source 150, where the driving signals are generated and directed toward the qudit 110 and the ancilla qubit 120.

Example Applications

The general quantum channels (i.e., CPTP maps) described in the present application may include 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 a, 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 λ_μ≥0 and Σ_μλ_μ=1, one form of the Kraus operators representing the stabilizing quantum channel is {K_i^μ=√{square root over (λ_μ)}|ψ_μi|}, 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 N_s stabilizer generators, the recovery is a CPTP map with Kraus rank 2^Ns. In some embodiments, the ancilla qubit may be used to sequentially measure all N_s 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:U_b(l)=U_l=P₊⊗Ŝ_l+P−⊗I,  (Eqn. 22)with Ŝ_l for the l-th stabilizer and P_±=½(|g±|e)(g|+e|), which is independent of the previous measurement outcomes b^(l-1). Finally, the correction unitary operation U_b(Ns) 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:

$\begin{array}{cc}W_{\uparrow }\rangle \equiv \frac{0\rangle +\sqrt{3}6\rangle }{2},W_{\downarrow }\rangle \equiv \frac{\sqrt{3}3\rangle +9\rangle }{2}.& \left(\mathrm{Eqn}.23\right)\end{array}$

For small loss probability γ for each excitation, this encoding scheme can correct errors up to O(γ²), which includes the following four relevant processes: identity evolution ({circumflex over (l)}), losing one excitation (â), losing two excitations (â²), and back-action induced dephasing ({circumflex over (n)}). Based on the Kraus representation of the QEC recovery (with Kraus rank 4), the following set of unitary operations U_b(l) is obtained for the construction of the QEC recovery channel with an adaptive quantum circuit:

$\begin{array}{cc}{\tilde{U}}_{Ø}=\left(\begin{array}{c}{\hat{P}}_3\\ \hat{I}-{\hat{P}}_3\end{array}\right),{\hat{U}}_0=\left(\begin{array}{c}{\hat{P}}_W\\ \hat{I}-{\hat{P}}_W\end{array}\right),{\hat{U}}_1=\left(\begin{array}{c}{\hat{P}}_1\\ \hat{I}-{\hat{P}}_1\end{array}\right),{\tilde{U}}_{00}=\left(\begin{array}{c}\hat{I}\\ \hat{0}\end{array}\right),{\tilde{U}}_{01}=\left(\begin{array}{c}{U}_{\hat{n}}\\ \hat{0}\end{array}\right),{\tilde{U}}_{10}=\left(\begin{array}{c}{U}_{\hat{a}}\\ \hat{0}\end{array}\right),{\tilde{U}}_{11}=\left(\begin{array}{c}{U}_{{\hat{a}}^2}\\ \hat{0}\end{array}\right),& \left(\mathrm{Eqn}.24\right)\end{array}$ where the projection are defined as {circumflex over (P)}_i≡Σ_k|3k+i3k+i| and {circumflex over (P)}_W≡|W_↑W_↑|+|W_↓W_↓|, and the unitary operators U_Ô (where Ô=â, â², {circumflex over (n)}) transform the error states Ô|W_σ back to |W_σ for σ=↑,↓. In other words:

$\begin{array}{cc}{U}_O=\sum _{\sigma }W_{\sigma }\rangle \frac{\langle W_{\sigma }{\hat{O}}^{†}}{\sqrt{\langle W_{\sigma }O^{†}OW_{\sigma }\rangle }}+U^{\perp },& \left(\mathrm{Eqn}.25\right)\end{array}$ 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⁽²⁾=(0,0), there is no error and the identify operation (Î) is sufficient. If b⁽²⁾=(0,1), there is back-action induced dephasing error, which changes the coefficients of Fock states so we need to correct for that with U_{{circumflex over (n)}}. If b⁽²⁾=(1,1), there is a single excitation loss, which can be corrected with U_â. If b⁽²⁾=(1,0), there are two excitation losses, which can be fully corrected with U_â2. 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:

$\begin{array}{cc}\rho \mapsto {\varepsilon }_{\mathrm{QI}}\left(\rho \right)=\sum _{\mu =1}^M{\varepsilon }_{\mu }\left(\rho \right)\otimes \mu \rangle \langle \mu ,& \left(\mathrm{Eqn}.26\right)\end{array}$ where |μμ| are orthogonal projections of the measurement device with M classical outcomes, ξ_μ are completely positive trace non-increasing maps, and Σ_μ=1^Mξ_μ(ρ) preserves the trace. Note that ξ_μ(ρ) gives the post-measurement state associated with outcome μ.

In some embodiments, the QI is implemented as follows. (1) Find the minimum Kraus representation for ξ_μ(ρ) (each with rank J_μ) with Kraus operators K_μ,j for j=1, 2, . . . , J_μ. (2) Introduce binary labeling of these Kraus operators, K_{{right arrow over (b)}(L)} where the binary label has length L=L₁+L₂, with the first L₁=┌log₂M┐ to encode μ and the remaining first L₂=┌log₂ max_μ(J_μ)┐ bits to encode j (padding with zero operators to make a total of 2^L Kraus operators). (3) Use the quantum circuit with L rounds of adaptive evolution and ancilla measurement. (4) Output the final state of the quantum system as well as b^(L1) that encodes μ associated with the M possible classical outcomes. In this way, the arbitrary QI described in Eq. 25 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, 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:

$\begin{array}{cc}\rho \mapsto {\varepsilon }_{\mathrm{POVM}}\left(\rho \right)=\sum _{\mu =1}^M\mathrm{Tr}\left[{\Pi }_{\mu }\rho \right]\mu \rangle \langle \mu ,& \left(\mathrm{Eqn}.27\right)\end{array}$ which is characterized by a set of Hermitian positive semidefinite operators {Π_μ}_μ=1^M that sum to the identify operator. In some embodiments, the positive semidefinite Π_μ is decomposed as Π_μ=Σ_jK_μ,j^†K_μ,j with a set of Kraus operators

${\left\{K_{\mu ,j}\right\}}_{j=1,\dots ,J_{\mu }}.$ 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. 10: (a) a standard quantum channel 1010 with the quantum system (e.g., qudit 1011) as the output and all detection results from the measurement of the ancilla qubit are discarded; (b) a POVM 1020 with the classical measurement outcomes 1021-1023 as the output and the quantum system discarded; and (c) a QI with both the quantum system 1031 and at least a portion of the classical measurement outcomes 1032-1033 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₁ bits of the ancilla measurement record. The remaining L₂ bits of the measurement record are discarded. In the QI 1030 of FIG. 10, L₁=2 and L₂=1. 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.

Other Considerations

Having thus described several aspects of at least one embodiment of this invention, it is to be appreciated that various alterations, modifications, and improvements will readily occur to those skilled in the art. Such alterations, modifications, and improvements are intended to be part of this disclosure, and are intended to be within the spirit and scope of the invention. Further, though advantages of the present invention are indicated, it should be appreciated that not every embodiment of the invention will include every described advantage. Some embodiments may not implement any features described as advantageous herein and in some instances. Accordingly, the foregoing description and drawings are by way of example only.

Various aspects of the present invention may be used alone, in combination, or in a variety of arrangements not specifically discussed in the embodiments described in the foregoing and is therefore not limited in its application to the details and arrangement of components set forth in the foregoing description or illustrated in the drawings. For example, aspects described in one embodiment may be combined in any manner with aspects described in other embodiments.

Use of ordinal terms such as “first,” “second,” “third,” etc., in the claims to modify a claim element does not by itself connote any priority, precedence, or order of one claim element over another or the temporal order in which acts of a method are performed, but are used merely as labels to distinguish one claim element having a certain name from another element having a same name (but for use of the ordinal term) to distinguish the claim elements.

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

The indefinite articles “a” and “an,” as used herein in the specification and in the claims, unless clearly indicated to the contrary, should be understood to mean “at least one.”

As used herein in the specification and in the claims, the phrase “at least one,” in reference to a list of one or more elements, should be understood to mean at least one element selected from any one or more of the elements in the list of elements, but not necessarily including at least one of each and every element specifically listed within the list of elements and not excluding any combinations of elements in the list of elements. This definition also allows that elements may optionally be present other than the elements specifically identified within the list of elements to which the phrase “at least one” refers, whether related or unrelated to those elements specifically identified.

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 ±5%.

The phrase “and/or,” as used herein in the specification and in the claims, should be understood to mean “either or both” of the elements so conjoined, i.e., elements that are conjunctively present in some cases and disjunctively present in other cases. Multiple elements listed with “and/or” should be construed in the same fashion, i.e., “one or more” of the elements so conjoined. Other elements may optionally be present other than the elements specifically identified by the “and/or” clause, whether related or unrelated to those elements specifically identified. Thus, as a non-limiting example, a reference to “A and/or B”, when used in conjunction with open-ended language such as “comprising” can refer, in one embodiment, to A only (optionally including elements other than B); in another embodiment, to B only (optionally including elements other than A); in yet another embodiment, to both A and B (optionally including other elements); etc.

As used herein in the specification and in the claims, “or” should be understood to have the same meaning as “and/or” as defined above. For example, when separating items in a list, “or” or “and/or” shall be interpreted as being inclusive, i.e., the inclusion of at least one, but also including more than one, of a number or list of elements, and, optionally, additional unlisted items. Only terms clearly indicated to the contrary, such as “only one of” or “exactly one of,” or, when used in the claims, “consisting of,” will refer to the inclusion of exactly one element of a number or list of elements. In general, the term “or” as used herein shall only be interpreted as indicating exclusive alternatives (i.e. “one or the other but not both”) when preceded by terms of exclusivity, such as “either,” “one of,” “only one of,” or “exactly one of” “Consisting essentially of,” when used in the claims, shall have its ordinary meaning as used in the field of patent law.

Also, the phraseology and terminology used herein is for the purpose of description and should not be regarded as limiting. The use of “including,” “comprising,” or “having,” “containing,” “involving,” and variations thereof herein, is meant to encompass the items listed thereafter and equivalents thereof as well as additional items.

Claims

1. A quantum information system comprising: an ancilla qubit;a qudit coupled to the ancilla qubit;a detector configured to generate a detection result based on a quantum state of the ancilla qubit;a driving source 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 qudit based on the detection result; anda controller coupled to the driving source and the detector and configured to: obtain an indication of a desired quantum channel;receive, from the detector, a detection signal indicative of the detection result;determine the at least one qudit driving signal based on the indication of the desired quantum channel;determine the at least one qubit driving signal based on the indication of the desired quantum channel;control the driving source to drive the qudit with the at least one qudit driving signal; andcontrol the driving source to drive the qubit with the at least one qubit driving signal.

2. The quantum information system of claim 1, wherein the indication of the desired quantum channel comprises a plurality of Kraus operators that form a minimal Kraus representation of the desired quantum channel.

3. The quantum information system of claim 2, wherein obtaining the indication of the desired quantum channel comprises determining the plurality of Kraus operators that form the minimal Kraus representation of the desired channel from one or more other operators.

4. The quantum information system of claim 3, wherein the controller is further configured to: determine a plurality of joint unitary operators, each of the plurality of joint unitary operators configured to act on a joint system comprising the ancilla qubit and the qudit; andcreate a binary tree structure associated with the plurality of unitary operators.

5. The quantum information system of claim 4, wherein the controller is further configured to determine each of the plurality of joint unitary operators based on a diagonalization of the sum of a subset of the plurality of Kraus operators that are rooted from an associated node of the binary tree structure.

6. The quantum information system of claim 5, wherein the controller is further configured to: determine a first qudit unitary operator, a second qudit unitary operator and a qubit unitary operator associated with and based on each of the plurality of joint unitary operators;determine the at least one qudit driving signal based on the first qudit unitary operator and the second qudit unitary operator; anddetermine the at least one qubit driving signal based on the qubit unitary operator.

7. The quantum information system of claim 1, wherein the ancilla qubit comprises a transmon qubit.

8. The quantum information system of claim 1, wherein the qudit comprises a quantum oscillator comprising electromagnetic radiation in a cavity.

9. A method of operating a quantum information system that includes a qudit coupled to an ancilla qubit forming a qudit-qubit system, the method comprising: applying a first unitary operation to the qudit-qubit system;generating a detection result based on a quantum state of the ancilla qubit;applying a second unitary operation to the qudit-qubit system based on the detection result;obtaining an indication of a desired quantum channel;determining at least one qudit driving signal based on the indication of the desired quantum channel; anddetermining at least one qubit driving signal based on the indication of the desired quantum channel.

10. The method of claim 9, wherein the indication of the desired quantum channel comprises a plurality of Kraus operators that form a minimal Kraus representation of the desired quantum channel.

11. The method of claim 10, wherein obtaining the indication of the desired quantum channel comprises determining the plurality of Kraus operators that form the minimal Kraus representation of the desired channel from one or more other operators.

12. The method of claim 11, further comprising determining a plurality of joint unitary operators, each of the plurality of joint unitary operators configured to act on a joint system comprising the ancilla qubit and the qudit; and determining a binary tree structure associated with the plurality of unitary operators.

13. The method of claim 12, further comprising determining each of the plurality of joint unitary operators based on a diagonalization of the sum of a subset of the plurality of Kraus operators that are rooted from an associated node of the binary tree structure.

14. The method of claim 13, further comprising: determining a first qudit unitary operator, a second qudit unitary operator and a qubit unitary operator associated with and based on each of the plurality of joint unitary operators;determining the at least one qudit driving signal based on the first qudit unitary operator and the second qudit unitary operator; anddetermining the at least one qubit driving signal based on the qubit unitary operator.

15. 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, wherein the method comprises: obtaining a plurality of Kraus operators associated with a desired quantum channel;generating a plurality of unitary joint operations associated with a joint qudit-qubit system based on the Kraus operators;determining a first unitary joint operation of a plurality of unitary joint operations to perform on the joint qudit-qubit system; anddetermining two unitary qudit-only operations and one unitary qubit-only operation based on the first unitary joint operation.

16. The at least one non-transitory storage medium of claim 15, wherein the method further comprises generating a binary tree structure associated with the plurality of unitary joint operations.

17. The at least one non-transitory storage medium of claim 16, wherein each node of the binary tree structure is associated with a respective unitary joint operation of the plurality of unitary joint operations.

18. The at least one non-transitory storage medium of claim 17, wherein each leaf of the binary tree structure is associated with a Kraus operator of the desired quantum channel.