CROSS REFERENCE

[0001] This International PCT Patent Application is based on and claims priority to U.S. Provisional patent Application No. 63/144,230, filed on February 1, 2021, the entirety of which is herein incorporated by reference.

GOVERNMENT FUNDING

[0002] This invention was made with government support under grant numbers 1640959 and 1936118 awarded by the National Science Foundation, grant number DE-SC0019406 awarded by the U.S. Department of Energy, grant numbers FA9550- 15-1-0015 and FA9550-19-1-0399 awarded by the U.S. Air Force Office of Scientific Research, and grant numbers W911NF-15-2-0067, W911NF-16-1-0349, W911NF-18- 1-0020 and W911NF-18-1-0212 awarded by the U.S. Army Research Office. The U.S. government has certain rights in the invention.

FIELD OF THE INVENTION

[0003] This disclosure relates to a method and system for controlling properties of electromagnetic resonators coupled to ancilla qubits for quantum information processing and quantum simulation. BACKGROUND OF THE INVENTION

[0004] Electromagnetic resonators (alternatively referred to as cavities) support electromagnetic resonance modes and may be used as a platform for applications such as quantum information processing and simulation of many-body quantum mechanical systems. These applications generally rely on quantum operations across a plurality of photon-number states of the cavity. Such quantum operations may be facilitated by highly nonlinear interactions among these photon-number states. However, native nonlinear interactions among photons in a cavity may be weak and untunable.

SUMMARY OF THE INVENTION

[0005] A toolbox is implemented and developed for photon-number dependent Hamiltonian engineering by off-resonantly driving ancilla qubit (s) of electromagnetic cavities. A general approach to design and optimize the drive fields for engineering arbitrary single-cavity target Hamiltonian and performing quantum operations is provided, with examples including three-photon interaction, parity-dependent energy, error-transparent Z-rotation for rotation-symmetric bosonic qubits, and cavity self-Kerr cancellation. This scheme is also generalized to implement error-transparent controlled-rotation between two cavities. The flexible and thus highly nonlinear engineered Hamiltonian for photons enables versatile applications for quantum simulation and quantum information processing. As examples, these schemes can be implemented with dispersively coupled microwave cavities and transmon qubits in the cQED platform.

BRIEF DESCRIPTION OF THE DRAWINGS

[0006] For a more complete understanding of the invention, reference is made to the following description and accompanying drawings, in which:

[0007] FIG. 1 illustrates an example system for engineering a photon-number dependent Hamiltonian in a cavity;

[0008] FIG. 2 illustrates an example microwave cavity coupled to an ancilla qubit;

[0009] FIG. 3 illustrates two coupled microwave cavities with three ancilla qubits;

[0010] FIG. 4 illustrates example logic flow for generating a photon-number dependent target Hamiltonians in a cavity system;

[0011] FIG. 5 illustrates example simulated energy shift as a result of engineering photon-number dependent Hamiltonians in a cavity system;

[0012] FIG. 6 illustrates simulated evolution of an example cat state of a cavity as a superposition of coherent states with opposite phases;

[0013] FIG. 7 illustrates simulated engineered energy shifts for error-transparent controlled-rotation in a cavity as a function of photon numbers in coupled cavities and a fidelity of a photon- number controlled π /8-gate as a function of a total cavity relaxation rate;

[0014] FIG. 8 illustrates example fidelity between states evolved by photon number dependent drive and states evolved by a target Hamiltonian as a function of time with or without ancilla qubit relaxation during a π /8-gate operation by abruptly turning on a photon number dependent drive with and smoothly turning on and off the PND drive with the drive ramping function being shown in dashed line.;

[0015] FIG. 9 illustrates example qubit-induced infidility of the photon- number dependent rotation gate as a function of rotation angle with fixed photon number dependent drive parameters and with a θ-dependent gate time TG = as a function of the gate time TG, with fixed parameters X0, and with Tc-dcpcndcnt drive amplitudes assuming that that the initial state is a ⁿ d ^a smooth ramping function for the drive is applied.;

[0016] FIG. 10 illustrates example Wigner function snapshots and cat fidelitiies with self-Kerr nonlinearity and with canceled self-Kerr nonlinearity;

[0017] FIG. 11 illustrates example cat state fidelity during a Kerr cancellation operation;

[0018] FIG. 12 illustrates example micromotion of the cat infidelity under photonnumber dependent Kerr cancellation; and

[0019] FIG. 13 illustrates cat infidelity induced by ancilla qubit errors during photon-number dependent Kerr cancellation. DETAILED DESCRIPTION

[0020] Quantum information processing or many-body quantum simulation applications may be implemented in a platform containing one or more electromagnetic cavities (alternatively referred to as electromagnetic resonators, or resonators for brevity throughout this disclosure). These resonators may support photon- number states in their cavity electromagnetic modes with controllable energies (or resonant frequencies). Logical quantum bits (qubits) may be coded as coherent superpositions of these photon-number states. Quantum operations on these logical qubits may be implemented for quantum information processing and quantum simulation purposes. Each resonator may be associated with an infinite-dimensional Hilbert space. Such an infinite-dimensional Hilbert space of a single resonator, for example, may enable flexible design of quantum error correction codes and result in extension of logical qubit lifetime and coherence time. As such, the rapidly growing Hilbert space of the resonators may be controlled to perform quantum operations for quantum information processing and for emulating the dynamics of the classically intractable many-body quantum systems. These quantum operations may be facilitated by some types of highly nonlinear interactions among the photon-number states of the resonators. However, native nonlinear interactions among photons in a resonator may be too weak and untunable for these quantum operations.

[0021] In the implementations of this disclosure described in more detail below, ancilla qubit (s) coupled to one or more resonators and driven by engineered external electromagnetic fields (alternatively referred to as controlled drives) are used to control the nonlinear interactions of photons in the one or more resonators to facilitate the various quantum operations as needed for a particular quantum information processing or quantum simulation task. Specifically, the frequency components and magnitude, phase, and timing of each frequency component of the external electromagnetic fields may be controlled to drive the ancilla qubit (s) to achieve an engineered interaction represented by a particular target photon-number dependent (PND) Hamiltonian. As such, a toolbox is implemented and developed for engineering a PND Hamiltonian for a resonator-ancilla qubit system. In some example implementations, a target PND Hamiltonian is engineered by off-resonantly driving the ancilla qubit (s) to minimize adverse effect from excitation of excited states of the ancilla qubit (s). A general approach to design and optimize the drive fields for engineering arbitrary a single-cavity target Hamiltonian and performing quantum operations is provided, with examples including target Hamiltonians for emulating three-photon interaction, parity-dependent energy levels, and error-transparent Z- rotation for rotation-symmetric bosonic qubits, and for achieving cancellation of cavity self-Kerr effect. This scheme is also generalized to implement error-transparent controlled-rotation between two coupled cavities. The flexible and thus engineered highly nonlinear Hamiltonian for photon interactions enables versatile applications for quantum simulation and quantum information processing.

[0022] FIG. 1 shows an example system 1 for implementing Hamiltonian engineering in an electromagnetic resonator platform. The system 1 includes a resonator- ancilla qubit system 2, a controller 3, one or more electromagnetic wave sources 4, one or more D.C. electric field sources 5, one or more magnetic field sources 6, and one or more detectors 7. In some implementations, the system 1 may include at least the resonator-ancilla qubit system 2 and the electromagnetic wave sources 4, and may further optionally include one or more of the other components above. [0023] In FIG. 1, the controller 3 may be of any form of dedicated circuits, general- purpose instruments, and mobile or fixed electronic devices including but not limited to desktop computer, laptop computers, tablets, mobile phones, personal digital assistants, and the like. The controller 3 may be configured to provide a user interface for controlling the Hamiltonian engineering in the resonator-ancilla qubit system 2. The controller 3 may include communication interfaces, processors, input/output (I/O) interfaces, storages, and display circuitry. The controller 3 may communicate with one or more of the other components of the system 1 to trigger operations and to send/receive commands, data, or the like, as shown by the various arrows in FIG. 1 connecting the controller 3 to other components of the system 1.

[0024] The electromagnetic wave sources 4 of FIG. 1 may be of any type of sources for generating electromagnetic radiation in various suitable and controllable spectral ranges. The electromagnetic wave sources 4 may include one or more signal generators, electromagnetic transducers/amplifiers, inductive/capacitive devices, terahertz wave generator, and/or optical emitters (e.g., arc lamps, lasers, light emitting diodes). In some implementations, frequency components of the output electromagnetic fields including frequencies, magnitudes, and phases may be adjusted/controlled via the electromagnetic wave sources 4. For example, the output electromagnetic fields may include radio-frequency waves and microwaves in, e.g., 10 KHz-300 GHz, inclusive. For another example, the output electromagnetic fields may include optical waves in the visible, infrared, and/or ultraviolet spectral ranges. For yet another example, the output electromagnetic fields may include terahertz waves. The output electromagnetic fields (alternatively referred to as control drives) may be directed to the resonator-ancilla system 2, via free space, transmission lines, microwave guides, optical fibers, or any other suitable means. The output electromagnetic fields may be further manipulated (e.g., focused, collimated, and the like using optical/microwave components) before being directed to the resonator-ancilla qubit system 2. The timing of the electromagnetic fields may be further controlled. For example, the electromagnetic field may be turned on and off abruptly or ramped up and down smoothly or adiabatically. The electromagnetic fields may be focused and may be spatially scanned to select different regions of the resonator-ancilla qubit system 2.

[0025] The system 1 may optionally include the D.C. electric or magnetic sources 5 and/or 6. The D.C. sources 5 and/or 6 may be used to generate and apply D.C. electric and/or magnetic fields to the ancilla qubits to modify the energy levels of the quantum states and/or to mix quantum states in the ancilla qubits. For example, a D.C. electric field may be applied for device characterization, e.g., for characterizing D.C. Stark effect on the various quantum states in the ancilla qubits. For another example, a magnetic field may be applied to induce energy level splitting. The D.C. electric and or magnetic fields for any of the purposes above may be applied to the entire or a particular region of the resonator-ancilla qubit system 2.

[0026] The resonator-ancilla qubit system 2 of FIG. 1 may include one or more resonators and one or more ancilla qubits. The resonators, for example, may be implemented in the form of microwave resonators, or optical resonators, or the like. A microwave cavity, for example, may be implemented as a planar or three-dimensional superconducting microwave resonator, as described in more detail in the examples given in relation to FIG. 2 and FIG. 3. The ancilla qubits may each include a device having at least two quantum levels. Each ancilla qubits may be coupled to one or more of the resonators in a particular manner. Specific examples are given in the disclosure below in relation to FIG. 2 and FIG. 3.

[0027] The detectors 7 of FIG. 1 may be configured to perform readout and/or other measurements of the resonator-ancilla qubit system 2. In the microwave spectral range, for example, the readout detector may include additional resonator circuits controllably coupled to the resonator-ancilla qubit system for measuring output microwave photons. The readout detector may be associated with other microwave signal filtering and mixing components. The readout detectors may be controlled via the controller 3. In optical spectral range, optical detectors may be implemented for measuring optical emission from the resonator-ancilla qubit system in the form of one or more devices including but not limited to photo diodes, photomultiplier tubes (PMT), microchannel plates, charge-coupled devices (CCDs), and complementary metal-oxide semiconductor (CMOS) sensors. The optical detectors may be configured to detect single or multiple optical photons. The optical detectors may include additional and or auxiliary components including power supplies, electric bias circuit, active or passive quenching circuit, and the like. The optical detector 140 may further include other optical components for collection and processing of the optical signal from the resonator-ancilla qubit system 2, including but not limited to spectral filters (such as spectrometers and dielectric/metal-film filters with suitable optical pass bands) and microscope optics (such as confocal optical microscope and multiple-photon microscope.

[0028] As an example for the resonator-ancilla qubit system 2 of FIG. 1, FIG. 2 shows a single microwave resonator 8 capacitively coupled to an ancilla qubit 11 via a coupling capacitor 14. The singlie microwave resonator 8, for example, may be implemented as an LC circuit including a capacitor 9 and an inductor 10. Such an LC resonator, for example, may be implemented as a planar or three-dimensional superconducting resonator.

[0029] The ancilla qubit 11 may be implemented as, for example, a superconducting phase, charge, or magnetic flux qubit, or their hybridization, such as a fluxonium, a transmon, a xmon, or a quantronium. These superconducting ancilla qubits may rely on one or more Josephson junctions for providing anharmonic electron energy levels. FIG. 2 shows an example transmon qubit that includes a capacitor 12 and an electron cooper-pair box in a form of capacitively shunted Josephson junction 13.

[0030] As another example for the resonator-ancilla qubit system 2 of FIG. 1, FIG. 3 shows two microwave resonators 15 and 16, denoted as resonator a and resonantor b, respectively. The microwae resonators 15 and 16 may be capacitively coupled to an ancilla qubit 17 (labeled as ancilla qubit a) and an ancilla qubit 19 (labeled as ancilla qubit b) via a coupling capacitors 18 and 20, respectively. The two microwave resonators 15 and 16 may further be coupled via a third ancilla qubit 21 (labeled as ancilla qubit c) and capacitors 22 and 23, as illustrated in FIG. 3. Each of the microwave resonators 15 and 16, for example, may be implemented as a superconducting LC circuit including a capacitor and an inductor, similar to the single microwave resonator 8 of FIG. 2. Each of the ancilla qubits 17 and 19 likewise may be implemented as, for example, a superconducting phase, charge, or magnetic flux qubit, or their hybridization, such as a fluxonium, a transmon, a xmon, or a quantronium. These superconducting ancilla qubits may rely on one or more Josephson junctions for providing anharmonic electron energy levels.

[0031] Other components may be coupled to each of resonator-ancilla qubit systems shown in FIG. 2 and FIG. 3. For example, additional superconducting transmission line resonators may be included for coupling controlled microwave signals for initialization, operation, driving, and readout of the resonators 8, 15, and 16, and ancilla qubits 11,17, 19, and 21. Metal traces in the entire resonator-ancilla qubit systems 2 of FIG. 2 and FIG. 3 may include one or more superconducting materials, e.g., alu- minum or niobium. The resonator-ancilla qubit systems 2 of FIG. 2 and FIG. 3 may be placed at a temperature lower than a threshold for achieving superconductivity in the one or more superconducting materials.

[0032] The resonator-ancilla qubit system shown in FIG. 2 and FIG. 3 may be implemented as a planar or three-dimensional circuit Quantum Electrodynsmics (cQED) system. The various superconducting components therein may be fabricated by various means. For example, these components may be lithographically patterned on a planar substrate.

[0033] Each of the microwave resonators may be associated with a cavity mode having a resonance frequency, as determined, for example by the inductance and capacitance of the resonator. In some implementations, the inductance and or capacitance of a resonator may be adjustable. Quantum mechanically, such a microwave resonator may be characterized by various photon-number states, each representing a number of microwave photons in the cavity mode. These photon-number states may be denoted by |0), |1), ..., |n), wherein n represents the number of photons in a particular photon-number state. Quantum information may be based on a logic qubit in such a resonator which may be coded as a coherent superposition of the photon-number states.

[0034] Each of the ancilla qubit above may include multiple quantum levels. Energy spacing between adjacent levels may be sufficiently different such that a transition between two of the multiple quantum levels may be selectively addressed by an external drive field without being significantly affected by the presence of other quantum levels. For example, as described above, the ancilla qubit may be implemented as a transmon, in which a normally harmonic cooper-pair box (having energy levels with equal spacing) is modified to include capacitive shunting of the Josephson junction therein to create anharmonicity in the energy levels. The two levels among the multiple levels may be selected as the basis states of the ancilla qubit. A corresponding microwave transition between the two levels is characterized by a microwave frequency corresponding to the energy spacing between these two quantum levels. The lower energy level of the two levels may be referred to as the ground state, denoted by \g), whereas the higher energy level of the two levels may be referred to as the excited state, denoted by |e). The resonant frequency between the ground state \g) and the excited state |e) of the uncoupled ancilla quibt may be represented by ω _q .

[0035] The coupling strength between a microwave resonator and an ancilla qubit as described above may be configured at various levels. For example, the coupling may be in the dispersive regime, as described in more detail below. In other words, the resonance frequency of the microwave resonator and the two-level transition frequency are sufficiently detuned such that the coupling between the resonator and the ancilla qubit is in the dispersive regime. Under such dispersive coupling, as described in more detail below, the resonant frequencies of the ancilla qubit between the ground state \g) and the excited state |e) may become a series of resonances ω _q,n each shifted from the its uncoupled resonant frequency by a amount x ⁿ - where x represents the dispersive coupling strength between the resonator and the ancilla qubit and n represents the photon-number in the resonator, i.e. , ω _q,n = ω _q — Xn. Quantum evolution of such a dispersively coupled resonator-ancilla qubit system may be described by a Hamiltonian including a dispersive interaction linear in photon-number. For another example, the coupling between the resonator and the ancilla qubit may be stronger and thus nonlinear interactions (e.g. quadratic to photon-number) of higher order than dispersive coupling may need to be included in the Hamiltonian. A cavity selfKerr term which is also quadratic to photon-number would also be included in the Hamiltonian.

[0036] The ancilla qubit may be further driven by external electromagnetic fields. Such external electromagnetic fields may be alternatively referred to as control fields or drive fields. Field parameters or knobs that may be controlled in the drive fields include, for example, the number of frequency components and the frequencies, magnitude, and phase of these frequency components, Further, timing of the drive field may be controlled. For example, the drive fields may be abruptly applied to the ancilla qubit. For another example, the drive fields may be gradually or adiabatically turned on and off. The presence of the drive fields modifies the Hamiltonian of the coupled resonator-ancilla qubit system by including additional interactions with the drive fields. By controlling the various field parameters of the drive fields, a desired target Hamiltonian due to the driven coupled ancilla qubit having a PND nonlinear photon-photon interaction may be engineered and the quantum state of the coupled resonator-ancilla system under the drive fields would evolve as dictated by the target Hamiltonian.

[0037] In some implementations, a specific target Hamiltonian may be generated to depend nonlinearly on photon-numbers in the resonator by adjusting the field parameters of the drive fields in the dispersive regime. Such a photon-number dependent (PND) Hamiltonian may include engineered nonlinear photon-photon interactions (that is otherwise not attainable via native nonlinear interactions in the coupled resonator-ancilla qubit system) underlying quantum evolution of various quantum operation in quantum information processing or quantum simulation.

[0038] In some implementations, the field parameters of the drive fields may be optimally determined for a particular target Hamiltonian such that the excitation and effect from excited population in the excited state |e) of the ancilla qubit may be minimized. By minimizing the excited state population in the ancilla qubit, dephasing of the resonator photon-number states caused by the decoherence of the ancilla qubit may also be minimized.

[0039] In some example implementations, the field parameters of the drive fields may be determined for a particular target Hamiltonian in which the self-Kerr nonlinearity of the resonator described above is canceled as a result of driving the ancilla qubit by the drive fields. The self-Kerr nonlinearity is usually undesired as it reduces fidelity of quantum operations based on the resonator photon-number states. As such, the cancellation of the self-Kerr nonlinearity in the resonator via the Hamiltonian engineering above help enhance the fidelity of the quantum operations.

[0040] The Hamiltonian engineering by adjusting the field parameters of the drive fields may be implemented in various manners. An example for determining field parameters of drive fields is shown in FIG. 4. In step 24, a target Hamiltonian due to coupling to the driven ancilla qubit and suitable for a desired quantum information processing operation or a desired quantum simulation may be determined. The target Hamiltonian may be a predetermined function of the photon-number states of the resonator.

[0041] In step 25, a finite set of possible frequency detunings for driving each of the series of transitions of the ancilla qubit coupled to the resonanters at ω _q,n = ω _q - x ⁿ may be determined. For example, a set of possible frequency detunings may include {±x/2, ±x/4}. The set of possible frequency detunings may include, for another example, {±y/2}. Another set of detunings that commensurate with x may be chosen such that an overall micromotion vanishes periodically, as described in further detail below. [0042] In step 26, choices of frequency detunings from the set of possible frequency detunings are assigned to each frequency component u _q.n = uj _q — x ⁿ of the drive fields and the field magnitudes are adjusted such that the target Hamiltonian is generated while a summation of the average ancilla qubit excited state population due to micromtion is minimized. Such an optimization process may also minimize the decoherence induced by ancilla qubit relaxation.

[0043] The optimization procedure described in FIG. 4 is merely illustrative. Other optimization procedures for obtaining the target Hamiltonian may be used with a general goal of minimizing the micromotion.

[0044] The implementations above may be generalized to coupled resonator-ancilla qubit system such as the example illustrated in FIG. 3 and other systems involving more than two coupled resonators and corresponding ancilla qubits. Specifically, a target Hamiltonian as a result of driving the ancilla qubits by the external drive fields that is dependent on photon-numbers in two or more resonators maybe be achieved while minimizing the average ancilla qubit excited state populations due to micromotion. An example implementation is provided below in more detail.

[0045] While the description above and details provided below use microwave resonators) coupled with transmon(s) as examples, the underlying principles are applicable to other types of resonators and ancilla qubits operating in the microwave spectral range, and other resonators and ancilla qubits operating in other electromagnetic spectral range, e.g., optical range.

[0046] The various sections below further provide more detailed implementations of the principles described above.

Dispersive model with off-resonant drives for Hamiltonian engineering [0047] The engineering of the Hamiltonian may be achieved in a single resonator dispersively coupled with an ancilla qubit. Using a dispersive model, the qubit-cavity system of FIG. 1 and FIG. 2 without external electromagnetic driving fields may be described by the Hamiltonian where is the frequency of the cavity mode a, u _q is the qubit transition frequency between qubit states and and x is the dispersive coupling strength. The effective qubit transition frequency is dependent upon the number state of the cavity, , with resonant frequencies Applying a time-dependent drive Ω(t) supplied from the electromagnetic wave sources 4 of FIG. 1 to the qubit with an interaction of the form and operating in the number-split regime where x is larger than the transition linewidth of both the qubit and the cavity, the qubit may be driven near selective number-dependent transition frequencies to address individual number states of the cavity. Compared to a scheme of imparting selective number-dependent arbitrary phases (SNAP) to photon Fock-states by directly exciting qubit transitions, a driving regime with large detuning to the effective number state dependent qubit transition resonance frequencies may be utilized to engineer a continuous photon-number dependent target Hamiltonian.

[0048] Specifically, the electromagnetic drive of the form may be used to drive the qubit. Moving the total Hamiltonian to an interaction picture with the unitary transformation the interaction becomes

Assuming Vm, a time-dependent perturbation theory may be used to derive an effective Hamiltonian, which governs the long time dynamics of the cavity-qubit system up to the initial and final kicks, as described in further detail below. Since the qubit is only driven off-resonantly with Vm, it can be assumed that it stays in its ground state. Moving back to the original frame, the effective Hamiltonian seen by the photon while the qubit stays in its ground state is

The off-resonant control drives on the ancilla qubit thus effectively generate a photonnumber dependent Hamiltonian for the cavity .

[0049] Rapidly-oscillating micromotion may be dictated by the kick operator further described below. The leading order kick operator is

[0050] To the first order in V, an initial state |n, g) under the kick operator above may evolve to at time t, showing an oscillating small population of the qubit excited state component with a time-averaged probability ^{where the scc} “ ond term include a contribution from the initial kick at t = 0. This excited state component can be viewed as coherent oscillations assuming a closed qubit-cavity system. If detunings commensurate with the dispersive coupling strength x is chosen for the driving electromagnetic felds, the overall micromotion vanishes at a period denotes the greatest common divi¬ sor among all the detunings and the dispersive shift, and averages to zero at long-time. For quantum gates implemented by PND Hamiltonian, the electromagnetic drive may be designed such that TQ = CTM for some c G N in order to achieve maximum gate fidelity. Alternatively, this constraint on TG may be relaxed by smoothly turning on and off the drive to remove the effect of the initial and the final kicks assuming that the drive is abruptly turned on at t = 0. Alternative situation where the drive is smoothly ramped up is described in more detail below.

PND Hamiltonian engineering

[0051] A desired PND Hamiltonian may be chosen as an target and an electromagnetic drive may be designed to generate the interactions that achieves the target Hamiltonian. For example, given a target Hamiltonian of the the form appropriate values of Q _m and 6 _m may be identified such that HE = H _T . The solution for Q _m and 8 _m for a given target Hamiltonian (with reasonable strengths ) may not unique and may be identified in various manners.

[0052] In some example implementations, various steps in FIG. 4 may be followed for determining the values of Q _m and 8 _m for the electromagnetic drive to achieve the target Hamiiltonian with minimal excitation of the ancilla qubit excited state. In the first step (labeled as step 25) of FIG. 4, a finite set of possible detunings, e.g., may be considered and the detunings of the various frequency component of the electromagnetic drive may be selected from this set of possible detunings. By selecting detunings commensurate with x, surprising near-resonant higher-order contributions to the ancilla qubit excited state population may be minimized, and it would also become easier to determine the periodicity at which the micromotion vanishes, e.g., TM = 8TT/X, for the chosen set of (or alternatively, is chosen as the set of possible de tunings). Those detunings are comparable to x which allows the largest possible engineered Hamiltonian strength.

[0053] In the second step (labeled as step 26) of FIG. 4, random choices of drive detunings from the possible set of detunings 8 _m 's for each number state may be made, and optimized parameters (e.g., field magnitude of various frequency components of the electromagnetic drive) that generate the target Hamiltonian may be determined according to Eq. (4) plus additional fourth-order perturbation theory terms while minimizing 52 _n Pn,e, the summation of the average qubit excited state probability due to micromotion. As described in further detail below, the optimized choice of detunings also minimizes the decoherence induced by ancilla qubit relaxation.

[0054] Generation of several example PND target Hamiltonians are provided below with numerical simulation results shown in FIG. 5. As shown inFIG. 5, a dispersive shift may be assumed as an example. The Hamiltonian strengths may be engineered up to with high precision. Even larger strengths z may be achievable even though they are subject to imperfections due to sixth- and higher-order terms in the perturbation theory. This energy scale of the engineered Hamiltonian, is much larger than typical cavity decoherence and thus may be favorable to achieve high fidelity quantum gates or to perform quantum simulation. (i) Example Target Hamiltonian: Photon-Photon Interaction

[0055] In one example application for implementing a PND Hamiltonian engineering, a tunable photon-photon nonlinear interaction may be generated to emulate dynamics of quantum many-body systems with cavity photons. Such nonlinearities are typically weak in native interactions. For example, a purely three-photon interaction for cavity photons may be engineered by setting

[0056] A simulation of the PND Hamiltonian engineering frequency for the three- photon interaction above is shown in 27 of FIG. 5, with the target Hamiltonian shown in dotted line with open markers and the engineered Hamiltonian shown by solid markers. Markers of different shape represent different engineered Hamiltonian strengths.

(ii) Example Target Hamiltonian: Parity-Dependent Energy

[0057] In another example implementation, a photon-number parity dependent Hamiltonian may be generated. Photon-number parity, for example, may serve as an error-syndrome in various bosonic quantum error correction codes such as cat codes and binomial codes. By engineering a Hamiltonian of the form the cavity can distinguish photon-number parity by energy, which may allow for a design of error-detection code or dynamical stabilization of the code states for bosonic quantum error correction.

[0058] A simulation of the energy of the photon-number parity dependent Hamiltonian above is shown in 28 of FIG. 5, with the target Hamiltonian shown in dotted line with open markers and the engineered Hamiltonian shown by solid markers. Markers of different shape represent different engineered Hamiltonian strengths.

(Hi) Example Target Hamiltonian: Error- Transparent Z-Rotation

[0059] In a resonantor, continuous rotation of an encoded logical qubit around the Z-axis may generate the whole family of phase shift gates Ro, including ^-gatc and Z-gate, which are common elements of single-qubit gates for universal quantum computing. For quantum information encoded in rotational-symmetric bosonic code that can correct up to d _n — 1 photon loss errors, the logical states may include with code-dependent coefficients f _n ’s. Phase shift gates at an angle 0 for logical states can be implemented via the cavity Kerr effect the or by four-photon interaction

[0060] To achieve fault-tolerant quantum computation, an error-transparent Hamiltonian that commutes with and is thus uninterrupted by the photon-loss error may be designed to perform continuous logical Z-rotations. By engineering the same positive energy shift hgR for and all of its recoverable error states while engineering an equal but opposite energy shift and all of its recoverable error states, the resulting Z-rotation may be ‘transparent’ to photon-loss-errors. Specifically, for cat codes or binomial codes with d _n = 2, the target Hamiltonian may be of the form [0061] For example, consider the This rotation can be implemented by applying Hz for a time t = g^, by imparting phase and phase with a SNAP gate, or by applying for a time The gate performance may be characterized in the presence of photon loss by performing the rotation gate on over the same gate time, followed by instantaneous single- photon-loss-error recovery in 30 of FIG. 5. Comparing the final fidelities, the PND gate and the SNAP gate show much higher resilience to photon-loss-error than due to their error-transparent structure.

[0062] A simulation of the energy of the error-transparent Z-rotation Hamiltonian above is shown in 29 of FIG. 5, with the target Hamiltonian shown in dotted line with open markers and the engineered Hamiltonian shown by solid markers. Markers of different shape represent different engineered Hamiltonian strengths.

Ancilla qubit-induced decoherence

[0063] In practice, the decoherence of the qubit may induce cavity dephasing during the PND process. Specifically, the qubit relaxation jump operator at a rate may cause dephasing for off-diagonal density matrix elements of the cavity number states at a rate while the qubit dephasing jump operator at a rate causes cavity dephasing at a rate

A choice of the optimized parameters for minimizing the micromotion also minimizes the decoherence induced by qubit relaxation, which is the dominant source of imperfection in typical cQED devices with a kHz-order Γ _q .

[0064] Smoothly turning on the PND drive may remove the contribution to p _n,e from the initial kick and further reduce the cavity dephasing. In FIG. 8, the 7r/8-gate operation via the abrupt PND drive and the smooth PND drive is compared. At the end of the gate operation, the simulated final gate fidelity is 99.929% for the abrupt drive and 99.934% for the smooth drive. The additional infidelity induced by ancilla relaxation is 0.075% for the abrupt drive and 0.055% for the smooth drive, showing a reduction in the qubit-induced cavity dephasing by using smooth ramping.

[0065] In contrast to the resonantly-driven SNAP gate which has an averaged qubit excited-state probability 1/2 during the operation, our scheme has a suppressed qubit excitation and thus has a much smaller decoherence rate during the operation. At the end of the gate operation, the overall qubit-induced decoherence for the PND gate scales as is the phase imparted on the number state while the qubit-induced overall decoherence for the SNAP gate scales as regardless of the phase (limited by ). In FIG. 9, the qubit-induced infidelity for RQ gate implemented by smooth PND drive is studied. The qubit-induced gate infidelity is proportional to the rotation angle G (and thus the total phase) while relatively independent of the gate time while 0 is fixed, as predicted.

[0066] The SNAP and PND schemes complement each other for photon-number- dependent operations. The SNAP gate is ideal for one-shot operation to impart large phases. On the other hand, the PND Hamiltonian engineering scheme is better suited for quantum simulation, continuous operation, and quantum gate with small phases. In 30 of FIG. 5, the 7r/8-gate fidelity modified by a lossy qubit in dashed lines is shown. The off-resonantly driven PND gate accumulates much less decoherence (qubit-induced infidelity=0.055%) than the SNAP gate (qubit-induced infidelity=0.91%), assuming no cavity relaxation. Since the qubit-induced decoher- ence for the PND gate is proportional to the imparted phase, the maximal qubit- induced PND Rg gate infidelity is 0.44% for 0 = 2-zr, which suggests that the PND scheme shall outperform the SNAP scheme (with the given gate time) for arbitrary error-transparent Re gate.

PND Hamiltonian engineering with Kerr

[0067] The description and implementations above are based on a dispersive model of the qubit-cavity coupling. In reality, the underlying microscopic model of coupled qubit-cavity system may also include higher-order coupling terms. For example, for a generalized model with photon self-Kerr K and second-order dispersive shift y', the Hamiltonian may be represented in the form of

[0068] Adding control electromagnetic drives Ω(t) = and assuming the time-dependent perturbation theory can again be used to determine an effective Hamiltonian similar to eq. (4) but with every nX replaced by due to the second-order dispersive shift.

[0069] The effective Hamiltonian seen by the photon while the qubit stays in its ground state is

[0070] The self-Kerr effect may be the leading order correction to cavity resonators that can cause unwanted rotations and (in the presence of photon loss can) introduce extra decoherence. The Hamiltonian engineering scheme described above may be implemented to cancel the cavity self-Kerr by choosing or to engineer a target Hamiltonian while canceling Kerr. Numerical simulation of the benifit of PND Kerr cancellation is presented in FIG. 6. Specifically, 31-36 of FIG. 6 show simulated evolution of the cavity cat state which rep resents an even superposition of coherent states with opposite phases. In particular, 31-33 of FIG. 6 shows simulated Wigner function snapshots of the cat state evolving under the cavity self-Kerr nonlinearity with different fidelity F, whereas 34-36 of FIG. 6 shows simulated Wigner function snapshots of the cat state with PND Kerr cancellation implemented as describe above corresponding to 31-33 of FIG. 6. Further, 37 of FIG. 6 further shows cat state fidelity as a function of time under PND Kerr cancellation and a signature of the micromotion. In the simulations above, it is assumed that cutoff photon-number and there is no photon loss. The cat state can be preserved with close to unit fidelity for t = 20/μs and 99.2 % fidelity for t = 100/zs with PND Kerr cancellation, as shown in 37 of FIG. 6.

PND Hamiltonian engineering for coupled cavities

[0071] The PND scheme may be further generalized to the case of coupled cavities, an example of which is shown in FIG. 3. Specifically, two coupled cavities may be used and the two cavity modes a and b may be dispersively coupled to their own ancilla qubits σ ^a and σ ^b . and to another joint qubit σ ^c with a dispersive shift y _c (assumed to be equal for both cavity modes), yielding a Hamiltonian (without an external electromagnetic drive) in the form of where w _a / _b are the frequencies of the cavities, ^ _q.a / _b / _c are the qubit transition frequencies between and are the dispersive coupling strengths. The coupled qubit <r ^c can be driven to control cavity states dependent on The qubits associated with the two cavities <J ^a and a ^b may likewise be driven to control cavity states dependent on n _a and n _b , respectively. Altogether, a two-cavity Hamiltonian may be engineered.

[0072] The generalized PND scheme described above may be applied to implement controlled-Z-rotations for realizing controlled-phase gates CPHASE(0), which are one class of essential two-qubit entangling gates for universal quantum computing. Using the PND scheme, an error-transparent operation of controlled-Z-rotation which is tolerant against photon loss in the cavities may be implemented.

[0073] For example, an error-transparent Hamiltonian H _C R for CPHASE(0) may be designed such that within a total number distance and its error states have the same negative energy shift hg _c ^ up to residual energy shifts on error states with total photon loss number exceeding [0074] FIG. 7 shows simulated engineered energy shift for error-transparent controlled- rotation as a function of the photon-numbers of coupled cavities. For example, the targeted energy shifts to implement H _C R for d _n = d _na = d _nb = 2 and the numerically simulated engineered energy shifts by the generalized PND as a function of the photon-numbers of coupled cavities n _a and n _b are shown in 38 of FIG. 7. The simulated fidelity of a PND controlled-π /8-gate as a function of the total cavity relaxation rate is shown in 39 of FIG. 7. The simmulated fidelity in 39 of FIG. 7 starting from the kitten-code encoded state followed by instantaneous single-photon-loss recovery in both cavities, is larger than 99.8% even in the presence of the relaxation of all three ancilla qubits. In 39 of FIG. 7, the dashed line represents the modifications due to ancilla qubit relaxation at a rate F _qa /2π = Γ _qb /2π = Fq _C /2π = 3 kHz. It is assumed above that Xa/^ = Xb/^ = X _c /2π = 2.56 MHz, the controlled-π /8-gate time TG = 16-π /x = 2TM , and that the two cavities have the same relaxation rate.

More details for the dispersive model

(i) Unitoary Evalution

[0075] A perturbative expansion of the unitary evolution for an off-resonantly driven, dispersively coupled qubit-oscillator system may be described by the Hamiltonian assumed that a tri-partition ansatz for the evolution operator Us(tf, ti) such that for an initial state and a final state where the evolution is separated into a time-independent effective Hamiltonian H _c g governing the long-time dynamics, as well as initial and final kicks, and The subscript S denotes evolution in the Schrodinger picture.

[0076] Move to the interaction picture with the unitary transformation U = exp(iH _o t/h the interaction term is where the subscript I denotes operators in the interaction picture, 0

The evolution operator in the interaction picture is connected to the Schrodinger picture one by

[0077] Assuming the time-dependent perturbation theory may be used to calculate Ui(tf, ti) in powers of VI and determine the perturbative expansion

[0078] Specifically,

For the zeroth order, It may be determined that and

Thus, it may be determined that

[0079] For

It may be obtained that

[0080] The third- and forth-order terms in the effective Hamiltonian are

[0081] Consider special cases where δ _m are commensurate with x, this problem reduces to a Floquet Hamiltonian with a single periodicity, and the Floquet effective Hamiltonian and the kick operator may be calculated to obtain identical results.

(ii) Unitary Evolution with Smooth Ramping

[0082] It has been assumed above that the drive is abruptly turned on at an initial time U and lasts till a final time t _f . Alternatively, a ramping function λ(t) such that H(t) = H0 + λ(t)V(t) may be applied to smoothly turn on (and off) the drive, which will remove the effect associated with the initial (and the final) kick operator if the ramping time scale is much longer than 1/x- The choice of the ramping function A(t) is not unique.

[0083] For mathematical simplicity, the case of applying a sinusoidal envelope A(t) = sin(γt) to a short-time gate operation from t = 0 to t = TG = π /γ may be considered. Using the time-dependent perturbation theory, it can be found that and

[0084] In the limit the resulting time evolution with this smooth sinusoidal envelope is thus equivalent to having an effective Hamiltonian generated by but without any initial or final kick effects. To compensate for the factor, the same gate (by accumulating the same phase) may be implemented as the abrupt version by rescaling the ramping function to or by letting the system evolve twice as long. In the abrupt version a gate time is chosen at which the micromotion vanishes, while with the sinusoidal envelope there is no such requirement because the micromotion has already been removed by the smooth ramping.

[0085] For long-time operation of the PND Hamiltonian engineering scheme, a ramp-up function A _up (t) = 0 — > 1 and a ramp-down function λ _down (t) = 1 — > 0 may be designed at the beginning and the end of the drive. An example ramp-up and ramp-down functions may be and otherwise. Here may be implemented as a special chosen value to generate the same accumulated phase as the abrupt case during the ramp-up and ramp-down periods. More Details of Evolution with Qubit-Induced Dephasing

[0086] In considering how errors in the ancilla qubit propagate to the cavity mode under off-resonant drives, the ancilla errors may be described by the qubit relaxation jump operator and the qubit dephasing jump operator time-dependent Lindblad master equation where p _to t is the total density matrix of the coupled qubit-oscillator system, and

Here k _a represents the relaxation rate of the cavity.

[0087] Moving to the interaction picture, d_ becomes while stays the same. Under a rotating wave approximation, when such that qubit decay releases a photon-number-dependent energy the relaxation jump operator may be treated as a set of independent jump operators in the cavity number state manifold,

[0088] Assuming again a tri-partition ansatz for the evolution superoperator such that there is a time-independent Liouvillian and a kick superoperator that absorbs the time dependence. For and may be expanded in perturbative orders of

[0089] The time-independent evolution superoperator may then be determined as where and the kick superoperator as with

[0090] Choosing δ _m commensurate with X such that all the time-dependent terms have a common periodicity T _M , then for tf = ti + cT _M for some integer for an Floquet generator

Taking ti = 0 and tracing over the ancilla qubit degree of freedom assuming p _gg = 1 and ρ _ee = p _ge = p _eg = 0, then the cavity density matrix in the interaction picture p _c,I follows a Floquet effective master equation with jump operators [0091] The jump operators cause dephasing for off-diagonal density matrix elements of the cavity number states p _n1n2 at a rate where is the time-averaged probability of the qubit excited state component The second term in p _n.e , is the contribution from the initial kick. Smoothly ramping up the drive can thus reduce the qubit-induced dephasing by removing the effect of the kick.

(iii) Optimized PND Parameters and additional simulations

[0092] In some inpiementations, optimized PND Hamiltonian engineering parameters that minimizes may be determined. All the engineered frequency shifts may be subject to a Fourier transformation precision of, e.g., ±0.5 kHz. The optimized parameters with real Ω _m 's may be obtained. The condition requiring real Ω _m ’s may be relaxed and complex values of Ω _m ’s may be solved, which may give rise to similar performance.

[0093] Additional simulations of fidelity for PND error-transparent Z-rotation is shown in FIG. 8 and FIG. 9. In 40 and 41of FIG. 8, the 7π/8-gate operation via abruptly turning on the PND drive with A(t) = 1 versus smoothly turning on and off the PND drive with are compared, showing that the smooth ramping scheme reduces the qubit-induced dephasing. More specifically in FIG. 8, the simulated fidelity between the states evolved by the PND drive and the states evolved by the target Hamiltonian Hz as a function of time is shown in solid curve (without ancilla qubit relaxation) and dashed curve overlayed by the solid curve (with ancilla qubit relaxation). The ramping function A(t) in presented as the filled orange curve in FIG. 8. The parameters in table 5. is used and it is assumed that kHz and TG = 16π /X = 3.125 /zs. The initial state is At the end of the gate operation, the final gate fidelity is 99.929% for the abrupt drive and 99.934% for the smooth drive. The additional infidelity induced by ancilla relaxation is 0.075% for the abrupt drive and 0.055% for the smooth drive.

[0094] In FIG. 9, the qubit-induced infidelity for the PDN Rg gate implemented by smooth PND drives is shown as being proportional to the rotation angle 0 (and thus the total phase) while relatively independent of the gate time while 0 is fixed, as predicted (values of the Rg gate infidelity with qubit relaxation or dephasing subtracted by the gate infidelity without qubit relaxation or dephasing). More specifically, 42 of FIG. 9 shows qubit-induced infidelity as a function of the rotation angle 0, with fixed PND drive parameters X ₀ , {Ω0}, {δ0} and with a θ-dependent gate time In 42 of FIG. 9, qubit-induced infidelity as a function of the gate time TG, with fixed parameters 0 = π /4, X0, {δ0}, and with TG-dependent drive amplitudes is shown. Here X0, {^o}> and {δo} are the parameters used in table 5.. It is assumed that F _Q /27F = = 3 kHz, and The initial state is ^(|0) _fc + |l) _fe ) ® \g), and a smooth ramping function is used for the drive.

[0095] The various tables below contain data used to generate the FIG. 5.

Table 1. Three-Photon Interaction: X/2TF = 2.56 MHz, = 0.5 kHz (27 of FIG. 5: square points) Table 2. Three-Photon Interaction: x/^ = 2.56 MHz, K3/2 π = 1 kHz (27 of FIG. 5: circular points)

Table 3. Parity-Dependent Energy: x/^ ⁼ 2-56 MHz, P = 20 kHz (28 of FIG. 5: square points)

Table 4. Parity-Dependent Energy: x/^ ⁼ 2-56 MHz, P = 40 kHz (28 of FIG. 5: circular points)

Table 5. Error-Transparent Z-rotation: X/2 π = 2.56 MHz, gR^ir = 20 kHz, d _n =2 (29 of FIG. 5: square points and 30 of FIG. 5) Table 6. Error-Transparent Z-rotation: d _n =2 (29 of FIG. 5: circular points)

More Details of Microscopic model and Kerr Corrections

[0096] The microscopic model of a resonator mode a coupled to another bosonic mode d with anharmonicity a is revisited below. Specifically,

[0097] For a small coupling g, perturbation theory may be used to estimate the frequency shifts as a function of photon-number in the resonator n _a and the anhar- monic mode n^ Expanding up to the order of g^ and keeping only na = 0, 1 (states I#), |e)), the generic Hamiltonian of the coupled system reads

Here

[0098] Consider an off-resonantly driven coupled system with the photon self-Kerr K and the second-order dispersive shift x' , here

[0099] A tri-partition ansatz for the time-evolution operator is again assumed, and move to the interaction picture with the unitary transformation the remaining drive term reads

[0100] The time-dependent perturbation theory may again be used to calculate in powers of to obtain the perturbative expansions with additional contributions from the Kerr term

[0101] For the zeroth order, it may be determined that

[0102] For It may be found that [0103] The third- and forth-order terms in the effective Hamiltonian are

AS = o, (57) where the last term satisfies the condition

(i) Optimized PND Parameters with Kerr and additional simulations

[0104] In this subsection tables of optimized PND Hamiltonian engineering parameters with the additional Kerr term are described. All the engineered frequency shifts are subject to a Fourier transformation precision of ±0.5 kHz.

[0105] Additional simulations of the cat state evolution under PND Kerr cancellation is shown in FIG. 10 and FIG. 11. It is assumed that the system starts with an even cavity cat state and the qubit in its ground state then simulate the state evolution in the rotating frame with With

PND Kerr cancellation, the cat state is preserved at a high fidelity ~ 99.2% even after a long time t = 100 ps. More specifically, Wigner function snap shots and cat fidelities under cavity self-Kerr effect, under PND Kerr cancellation, and and under PND Kerr cancellation with a lossy ancilla qubit and loss cavity are shown in 44, 45, 46, and 46-2 of FIG. 10, respectively, with different fidelities. An abrupt PND drive with parameters shown in table 7. is used, and it is assume as an example a qubit relaxation rate of and a cat size of The snap shots are chosen at multiples of the micromotion period

[0106] Cat state fidelity during the Kerr cancellation operation by abruptly turning on the PND drive and smoothly turning on and off the PND drive is described in Eq. (34) and (35) and are shown in 47 and 48 of FIG. 11, respectively. The cat fidelity as a function of time is shown as solid curves (without ancilla qubit relaxation) and dashed curves (with ancilla qubit relaxation). The ramping function A(t) in shown as the shaded dash curve. Here, parameters shown in table 7. are used as an example, and it is assumed that and At the end of the operation, the final cat fidelity is 99.180% for the abrupt drive and 99.184% for the smooth drive. The additional infidelity induced by ancilla relaxation is 2.568% for the abrupt drive and 2.276% for the smooth drive. The smooth ramping scheme reduces the qubit-induced dephasing as predicted. In the middle of the operation (after the ramp-up period and before the ramp-down period), the cat state fidelity is higher for the case of abrupt drive at each micromotion period because the final kick is canceled out by the initial kick.

[0107] In FIG. 12 and FIG. 13, the micromotion and the qubit-induced infidelity for PND Kerr cancellation are shown. It is shown that the amplitude of the micromotion and the qubit-induced infidelity both scale as as predicted. The Kerr cancellation with y' = 2K is used as a special example case such that both the cavity self-Kerr and the second-order dispersive shift are cancelled out by the engineered Hamiltonian and thus there is a perfect micromotion periodicity TM unperturbed by the total energy.

[0108] More specifically, micromotion of the cat state infidelity under PND Kerr cancellation is shown in FIG. 12. In this table of plots, the infidelity (x-axis, %) of a cat state as a function of time (y-axis, /zs) with varying sets of parameters and is shown. The insets indicate the corresponding cavity self-Kerr value relative to K with Here, it is assumed as an example that abrupt PND drive is used, and and are the parameters used in table 7.. It may be inferred that the micromotion amplitude scales as and that the micromotion period is as predicted.

[0109] Further, cat infidelity induced by ancilla qubit errors during PND Kerr cancellation are shown in FIG. 13. In particular, 49 of FIG. 13 shows values of cat infidelity at time t = 12 μs with qubit relaxation at a rate _Q subtracted by the cat infidelity without qubit relaxation. The left panel has a fixed Ω = Ω0 with varying x, and the right panel has a fixed X = Xo with varying Q. 50 of FIG. 13 shows values of cat infidelity at time t = 12 /zs with qubit dephasing at a rate r^,/27r = 3kHz subtracted by the cat infidelity without qubit dephasing. The left panel has a fixed Q = Qo with varying x, and the right panel has a fixed X ⁼ Xo with varying Q. Here it is assumed a abrupt PND drive of are the parameters used in table 7., and a cat size

Table 7. Kerr Cancellation: (FIG. 6 (b-c)) Table 8. Error- Transparent Z-rotation with Kerr Cancellation: X/2π = 2 MHz, g^ = 20 kHz, K/2π = 3 kHz, 2π = 6 kHz

More Details for Hamiltonian Engineering for Two Coupled Cavities

[0110] Consider two cavity modes d and b dispersively coupled to two ancilla qubits σ ^a and σ ^b respectively, and to another qubit σ ^c jointly with a dispersive shift y _c , assumed as an example to be equal for both modes,

[0111] Drive term can be added as [0112] The fourth order terms can be found in the description above. The engi- neered Hamiltonian is

(i) Error-transparent controlled- Z-rotation

[0113] The form of the engineered Hamiltonian does not have the full degree of freedom to create arbitrary structure of E _nanb but is enough to design error-transparent controlled-rotation along the Z axis for implementing CPHASE gate on rotational- symmetric bosonic codes. Specifically, by solving for target frequency shifts which takes { for all Within a total number dis tance the overall Hamiltonian should accumulate the same phase for and its error states while keeping and their error states unchanged.

[0114] Further shown below are tables of optimized parameters for implementing such error-transparent controlled-rotation with = 2. All the engineered frequency shifts are subject to a Fourier transformation precision of ±0.5 kHz.

Table 9. Error- Transparent Controlled- Rotation (drives on qubit σ ^a ): Xα/2π = 2.56 MHz

Table 10. Error- Transparent Controlled-Rotation (drives on qubit σ ^b ): X _b /2π = 2.56 MHz

Table 11. Error- Transparent Controlled-Rotation (drives on qubit σ ^c ): X _c /2π = 2.56 MHz [0115] In summary, a toolbox is implemented and developed for photon- number dependent Hamiltonian engineering by off-resonantly driving ancilla qubit (s). A general formalism to design and optimize the control drives for engineering arbitrary singlecavity target Hamiltonian and performing quantum gates is provided, with examples including three-photon interaction, parity-dependent energy, error-transparent Z-rotation for rotation-symmetric bosonic qubits, and cavity self-Kerr cancellation. This scheme is also generalized to implement error-transparent controlled-rotation between two cavities. The flexible and thus highly nonlinear engineered Hamiltonian for photons admits versatile applications for quantum simulation and quantum information processing. These schemes can be implemented with dispersively coupled microwave cavities and transmon qubits in the cQED platform.

[0116] The description and accompanying drawings above provide specific example embodiments and implementations. Drawings containing device structure and composition, for example, are not necessarily drawn to scale unless specifically indicated. Subject matter may, however, be embodied in a variety of different forms and, therefore, covered or claimed subject matter is intended to be construed as not being limited to any example embodiments set forth herein. A reasonably broad scope for claimed or covered subject matter is intended. Among other things, for example, subject matter may be embodied as methods, devices, components, or systems. Accordingly, embodiments may, for example, take the form of hardware, software, firmware or any combination thereof.

[0117] Throughout the specification and claims, terms may have nuanced meanings suggested or implied in context beyond an explicitly stated meaning. Likewise, the phrase “in one embodiment/implementation” as used herein does not necessarily refer to the same embodiment and the phrase “in another embodiment/implementation” as used herein does not necessarily refer to a different embodiment. It is intended, for example, that claimed subject matter includes combinations of example embodiments in whole or in part.

[0118] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of skill in the art to which the invention pertains. Although any methods and materials similar to or equivalent to those described herein can be used in the practice or testing of the present invention, the preferred methods and materials are described herein.

[0119] In general, terminology may be understood at least in part from usage in context. For example, terms, such as “and”, “or”, or “and/or,” as used herein may include a variety of meanings that may depend at least in part on the context in which such terms are used. Typically, “or” if used to associate a list, such as A, B or C, is intended to mean A, B, and C, here used in the inclusive sense, as well as A, B or C, here used in the exclusive sense. In addition, the term “one or more” as used herein, depending at least in part upon context, may be used to describe any feature, structure, or characteristic in a singular sense or may be used to describe combinations of features, structures or characteristics in a plural sense. Similarly, terms, such as “a,” “an,” or “the,” may be understood to convey a singular usage or to convey a plural usage, depending at least in part upon context. In addition, the term “based on” may be understood as not necessarily intended to convey an exclusive set of factors and may, instead, allow for existence of additional factors not necessarily expressly described, again, depending at least in part on context.

[0120] Reference throughout this specification to features, advantages, or similar language does not imply that all of the features and advantages that may be realized with the present solution should be or are included in any single implementation thereof. Rather, language referring to the features and advantages is understood to mean that a specific feature, advantage, or characteristic described in connection with an embodiment is included in at least one embodiment of the present solution. Thus, discussions of the features and advantages, and similar language, throughout the specification may, but do not necessarily, refer to the same embodiment.

[0121] Furthermore, the described features, advantages and characteristics of the present solution may be combined in any suitable manner in one or more embodiments. One of ordinary skill in the relevant art will recognize, in light of the description herein, that the present solution can be practiced without one or more of the specific features or advantages of a particular embodiment. In other instances, additional features and advantages may be recognized in certain embodiments that may not be present in all embodiments of the present solution.