A quantum system for stabilizing a bosonic qubit

EP4720939A1Pending Publication Date: 2026-04-08ALICE & BOB
View PDF 0 Cites 0 Cited by

Patent Information

Authority / Receiving Office
EP · EP
Patent Type
Applications
Current Assignee / Owner
Filing Date
2024-03-05
Publication Date
2026-04-08

AI Technical Summary

Technical Problem

Current methods for stabilizing squeezed cat qubits face challenges due to the need for multiple Hamiltonian terms and high energy injection, which complicates the system's stability and coherence, making it difficult to achieve the desired noise bias for large-scale quantum algorithms.

Method used

A quantum system comprising a non-linear superconducting quantum circuit with an asymmetrically threaded superconducting quantum interference device, utilizing three microwave sources to engineer a Hamiltonian that stabilizes a two-dimensional manifold for a bosonic qubit, reducing the number of required pumps and enhancing noise bias.

Benefits of technology

The system achieves improved noise bias and reduced non-adiabatic errors, enabling more efficient bit-flip suppression and potentially supporting large-scale quantum algorithms with enhanced stability and coherence.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure EP2024055758_28112024_PF_FP_ABST
    Figure EP2024055758_28112024_PF_FP_ABST
Patent Text Reader

Abstract

The invention relates to a quantum system (1) for stabilizing a bosonic qubit comprising: - a command circuit (5) including three microwave sources (11, 13, 15), and - a non-linear superconducting quantum circuit (3) including a four-wave mixing non-linear element (7) and a resonant portion (9), the quantum circuit (3) having a first and a second modes with a first and a second resonant frequencies f a , f b respectively. The quantum circuit (3) is arranged to engineer a Hamiltonian H expressed as H / ħ =g2(a2 + λa†+a – a2)b† + h. c., the Hamiltonian H yielding through dissipation of the second mode an effective dissipator D[a2 + λa†+a – a2] which stabilizes the bosonic qubit. To this end, the first microwave source (11) delivers radiation at a frequency equal to |2f a − f b |, the second microwave source (13) delivers radiation at a frequency equal to f b , and the third microwave source (15) delivers radiation at a frequency equal to f b to drive the second mode.
Need to check novelty before this filing date? Find Prior Art

Description

[0001]A quantum system for stabilizing a bosonic qubit The field of the invention relates to stabilizing a bosonic qubit. Generally speaking, superconducting qubits can be implemented as two-level systems of superconducting electronic circuits. Such qubits can be stored in bosonic modes and thus form a particular class of superconducting qubits known as bosonic qubits. It is known from the state of the art that it is possible to stabilize cat qubits, i.e. bosonic qubits defined by a quantum manifold spanned by the superpositions of two coherent states, which are quasi-classical states of the bosonic mode. To this end, a specific dissipative stabilization mechanism may be performed and consists in engineering a non- linear conversion between two photons of a first mode – also known as memory mode or cat qubit mode – that hosts the stabilized quantum manifold and one photon of a second mode – also known as buffer mode – that is strongly dissipative. The cat qubits offer the advantage of having a bit-flip rate that decreases exponentially while the phase-flip rate increases only linearly. As a result, stabilized cat qubits benefit from a high noise bias, which means that the bit-flip probability is much smaller than the phase-flip probability. Due to their noise structure, cat qubits can be concatenated with a repetition code instead of the usual surface code used for usual superconducting qubits. However, it appears that the noise bias of the cat qubits may not be sufficient to use a simple repetition code for some large-scale algorithms and there is therefore a need to find bosonic qubits with an even higher noise bias. For this purpose, recent studies have focused on a particular class of cat qubits: the squeezed cat qubits. For instance, Q. Xu et al. (2022) have proposed, in the article “Autonomous quantum error correction and fault-tolerant quantum computation with squeezed cat qubits” (arXiv:2210.13406) to implement an autonomous quantum error correction (AQEC) scheme using squeezed cat code against the dominant error source, excitation loss, in continuous-variable systems. This article announces a noise bias of about 10-15with only a mean number of four photons to encode the cat-states. As a further example, T. Hillmann et al. (2022) propose and analyze, in the article “Quantum error correction with dissipatively stabilized squeezed cat qubits” (arXiv:2210.13359, Phys. Rev. A 107, 032423), the error correction performance of a dissipatively stabilized squeezed cat qubit. However, the possibility of stabilizing a squeezed cat qubit faces several obstacles in practice and thus remains purely theoretical as of today. In particular, such a stabilization would require to engineer a multitude of Hamiltonian terms to yield the desired dissipator and thus to use too many pumps to the point of threatening the stability of the system. It can be added that one of these terms, where ^^ is the photon annihilation operator of the first mode and ^^ is the photon annihilation operator of the second mode, requires the injection of a high amount of energy in the system that jeopardizes the coherence of the system. The present invention seeks to improve the situation. To this end, the Applicant proposes a quantum system for stabilizing a bosonic qubit comprising: - a command circuit including a first microwave source, a second microwave source and a third microwave source each arranged for delivering microwave radiation, and - a non-linear superconducting quantum circuit including an asymmetrically threaded superconducting quantum interference device and at least one resonant portion to which said asymmetrically threaded superconducting quantum interference device is connected, said asymmetrically threaded superconducting quantum interference device having flux lines through which radiation can be delivered to modulate a common flux and a differential flux, said non-linear superconducting quantum circuit having a first mode with a first resonant frequency and a second mode with a second resonant frequency, said second resonant frequency being different from twice the first resonant frequency, said second mode being dissipative. The asymmetrically threaded superconducting quantum interference device is arranged such that, when ^ the first microwave source delivers radiation through said flux lines to modulate the common flux at a frequency equal to the absolute value of the difference between twice the first resonant frequency and the second resonant frequency, ^ the second microwave source delivers radiation through said flux lines to modulate the common flux at a frequency equal to the second resonant frequency, and ^ the third microwave source delivers radiation at a frequency equal to the second resonant frequency to the a least one resonant portion to drive the second mode, the non-linear superconducting quantum circuit engineers a Hamiltonian H expressed as + h. c., where g2is proportional to the amplitude of the source, annihilation operator of the first mode, ^^ is a complex number which phase and amplitude result from the amplitude of the second microwave source and are each defined relative to the phase and amplitude of the first microwave source respectively, ^^ is a complex number resulting from the third microwave source and defined relative to the first microwave source, ^^ is the annihilation operator of the second mode, and h. c. is the Hermitian conjugate, said Hamiltonian H yielding through dissipation of the second mode an effective dissipator ^^[^^2^^ ^^ − ^^2]which stabilizes a two-dimensional manifold hosting the bosonic qubit in the first mode. According to one or more embodiments, the first microwave source and the second microwave source are configured such that the respective phases of the first microwave source and the second microwave source are substantially equal to each other. According to one or more embodiments, the second microwave source and the third microwave source are configured such that the respective phases of the second microwave source and the third microwave source are substantially equal to each other. According to one or more embodiments, the first microwave source, the second microwave source and the third microwave source are configured such that the respective phases of the first microwave source, the second microwave source and the third microwave source are substantially equal to each other. According to one or more embodiments, the third microwave source is arranged to deliver radiation through the flux lines to modulate the differential flux at a frequency equal to the second resonant frequency to drive the second mode. According to one or more embodiments, the second microwave source is arranged to deliver radiation whose amplitude causes ^^ to have an amplitude less than or equal to 1. According to one or more embodiments, the second microwave source is arranged to deliver radiation whose amplitude causes ^^ to have an amplitude substantially equal to 1. According to one or more embodiments, the respective phases of the first mode and the second mode across the asymmetrically threaded superconducting quantum interference ^^ device each have a zero-point fluctuation ^^^^, ^^^^, and the ratio^^^^^^is less than 3. Furthermore, the Applicant also proposes a method for stabilizing a bosonic qubit performed by the quantum system described above and comprising the following operations: - delivering, by the first microwave source, radiation through the flux lines to modulate the common flux at a frequency equal to the absolute value of the difference between twice the first resonant frequency and the second resonant frequency, - delivering, by the second microwave source, radiation through the flux lines to modulate the common flux at a frequency equal to the second resonant frequency, - delivering, by the third microwave source, radiation at a frequency equal to the second resonant frequency to the a least one resonant portion to drive the second mode. The implementation of this method causes the non-linear superconducting quantum circuit to engineer a Hamiltonian H expressed as + h. c., where g2is proportional to the amplitude of source, is the annihilation operator of the first mode, ^^ is a complex number which phase and amplitude result from the amplitude of the second microwave source and are each defined relative to the phase and amplitude of the first microwave source respectively, ^^ is a complex number resulting from the third microwave source and defined relative to the first microwave source, ^^ is the annihilation operator of the second mode, and h. c. is the Hermitian conjugate, said Hamiltonian H yielding through dissipation of the second mode an effective dissipator ^^[ ^^2^^ ^^ − ^^2] which stabilizes a two-dimensional manifold hosting the bosonic qubit in the first mode. Other features and advantages of the invention will become apparent from the following description provided for indicative and non-limiting purposes, with reference to the accompanying drawings, wherein: - Figure 1 schematically illustrates a quantum system comprising a non-linear superconducting quantum circuit and a command circuit according to the invention, - Figure 2 illustrates the respective Wigner functions of a standard cat qubit and a bosonic qubit stabilized by the quantum system of figure 1, - Figure 3 illustrates an explicit diagram of the quantum system of figure 1, - Figure 4 illustrates a partial electrical equivalent diagram of a first galvanic example implementation of the quantum system of figure 3, - Figure 5 illustrates a partial electrical equivalent diagram of a second galvanic example implementation of the quantum system of figure 3, - Figure 6 illustrates a partial electrical equivalent diagram of a capacitive example implementation of the quantum system of figure 3, - Figure 7 illustrates the bit-flip rate of the bosonic qubit stabilized by the quantum system of figure 1 as a function of the number of photons, - Figure 8 illustrates the non-adiabatic errors of the quantum system of figure 1 resulting from the implementation of a Z gate, - Figure 9 illustrates a comparison of the non-adiabatic errors between the bosonic qubit stabilized by the quantum system of figure 1 and a theoretical squeezed cat qubit, and - Figure 10 illustrates the bit-flip rate of the bosonic qubit stabilized by the quantum system of figure 1 as a function of the ratio of the respective zero-point fluctuations of the superconducting phase of the bosonic qubit mode and the buffer mode. The drawings and the following description are comprised for the most part of positive and well-defined features. As a result, they are not only useful in understanding the invention, but they can also be used to contribute to its definition, should the need arise. A. Cat qubits So far, the Applicant's work has been generally concerned with the stabilization of cat qubits. The possibility of realizing such a qubit has been demonstrated by R. Lescanne et al. in “Exponential suppression of bit-flips in a qubit encoded in an oscillator” (Nature Physics, 2020) demonstrated that such cat qubits can be stabilized with a non-linear conversion between two photons of a first mode a – memory mode or cat qubit mode – and one photon of a second mode b – buffer mode. A cat qubit is defined as a two-dimensional manifold spanned by the so-called cat states | ^^^which are superpositions of two coherent states | ^^^and |− ^^^: ^^^^^±^^^ with: Stabilized cat qubits are known to benefit from a high noise bias, which means that the bit-flip probability is exponentially smaller than the phase-flip probability. More precisely, an effective error channel (e.g., bit errors or "bit-flips") is suppressed in an exponential way with the "size" – i.e. the average number of photons ^̅^ = | ^^|2– of the Schrödinger cat states of the cat qubits. As previously mentioned, this exponential suppression of bit-flip errors is only at the cost of linear increase of phase-flip errors. According to current knowledge, this suppression should apply to a large class of physical noise processes having a local effect on the phase space of a harmonic oscillator. This includes, but is not limited to, photon loss, thermal excitations, photon dephasing, and various nonlinearities induced by coupling to a Josephson junction. Recent experiments in the context of quantum superconducting circuits have observed this exponential suppression of bit-flip errors with the average number of photons in the cat states. A.1 Stabilization schemes The cat qubits can be stabilized or confined by the following exemplary schemes: a) a parametric dissipative stabilization, with jump operator ^^2=√^^2( ^^2− ^^2), where ^^2is the two-photon dissipation rate, ^^ is the photon annihilation operator of the memory mode a and ^^ is a complex number defining the cat qubit. This jump operator can be realized by coupling a lossy buffer mode b with dissipation rate ^^^^, and a four-wave mixing device – typically a Josephson junction or an ATS – to the cat qubit mode a and by engineering the Hamiltonian = ^^2( ^^2− ^^2) ^^†+ h. c., where ^^ is the photon annihilation operator of the b and ^^2is the two-photon coupling rate, by applying to the four-wave mixing device a pump at frequency |2 ^^^^− ^^^^| and a drive of the buffer mode b at frequency ^^^^provided ^^2< ^^^^. b) a Kerr Hamiltonian^where ^^ is the amplitude of the Kerr Hamiltonian, ^^ is operator, and | ^^|2is the mean photon number. c) a detuned Kerr Hamiltonian^where ^^ is the amplitude of the Kerr operator, ^^ is a complex number defining the cat qubit, and ^^ is the detuning factor. d) a two-photon exchange (TPE) Hamiltonian = ^^2( ^^2− ^^2) ^^++ h. c., where ^^2 is the complex two-photon coupling rate, ^^ photon annihilation operator, ^^ is a complex number defining the cat qubit, and ^^±are the lowering and raising operators of the two-level system. This Hamiltonian can be engineered in the same way as the parametric dissipative stabilization a). e) a DC dissipative stabilization, for which the Applicant filed the European patent application EP 23306839.4, in which a cat qubit is stabilized in the spirit of the previous stabilization scheme a), except that the two-photon pump that engineers the non-linear conversion between two photons of the memory mode a and one photon of the buffer mode b is replaced with a DC voltage source which biases a non-linear element formed exclusively of one or more Josephson junctions; the two-photon coupling rate g2is therefore not limited by the amplitude of the two-photon pump: ^^^^2 =^^ 24 ^^^^^^^^, where ^^^^is the Josephson energy of the one or more Josephson junctions, ^^^^is the zero-point fluctuation of the phase of the memory mode a, and ^^^^is the zero-point fluctuation of the phase of the buffer mode b. f) a resonant dissipative stabilization, with jump operator ^^2=√^^2(^^2− ^^2), where ^^2is the two-photon dissipation rate, ^^ is the photon annihilation operator of the memory mode a and ^^ is a complex number defining the cat qubit. This jump operator can be realized by coupling a lossy buffer mode b with dissipation rate ^^^^, and a three-wave mixing device to the cat qubit mode a which engineers the Hamiltonian = ^^2( ^^2− ^^2) ^^†+ h. c., where ^^ is the photon annihilation operator of the buffer mode b provided the mode frequencies verify substantially 2 ^^^^= ^^^^and ^^2< ^^^^to which a drive of the buffer mode at frequency ^^^^is added. The following subsection focuses on a particular stabilization scheme: the dissipative squeezing stabilization. A.2 The squeezed cat qubit Cat qubits have the property to allow the realization of quantum gates such as the Z or the CNOT gate while maintaining the noise bias. The characteristics which are used to assess the quality of a quantum gate are its execution time, i.e. the time for the quantum gate to operate, and the associated error probability. For cat qubits, the phase-flip error probability at the optimal time of the gate depends on the ratio^^1⁄^^2, where ^^1is the single photon loss rate – or error rate since it leads to phase-flip errors – of the cat qubits used to perform the quantum gate, and ^^2is the two-photon dissipation rate – also called correction or stabilization rate – of these cat qubits. The smaller^^1⁄^^2is, the better the gate fidelity gets. It is generally considered that the use of a single repetition code is sufficient to correct the remaining errors as bit-flip errors are sufficiently rare; and, more particularly, a phase- flip error correction code is sufficient to correct the remaining phase-flip. This can be for instance a repetition code defined in the dual base or any other state-of-the-art error correction code. However, even if the cat qubit has a high noise bias, which is a particularly advantageous feature in the perspective of realizing a quantum error correction and thus designing a reliable quantum computer, it appears that such a noise bias might not be enough for large-scale quantum algorithm with a simple repetition code. In this context, it was recently highlighted that a particular family of cat qubits have an even higher noise bias: the squeezed cat qubits. The squeezed cat states also have less non-adiabatic errors during the period of time of a quantum gate, such as a Z gate or a CNOT gate. Overall, reducing the non-adiabatic errors allows to reduce the requirement on the ratio^^1⁄^^2. The squeezed states | ^^, ^^^are defined as: | ^^, ^^^ = ^^( ^^) ^^( ^^)|0^ with ^^( ^^) the displacement operator: ^^(^^)= ^^^^ ^^†− ^^∗ ^^And ^^(^^)the squeezing operator: where: ^^ is the complex squeezing parameter with norm ^^ = | ^^|. The squeezed cat states are defined as the coherent superposition of two squeezed states with opposite displacement amplitudes and identical squeezing: where: ^^ is a normalization constant. Note that to have proper squeezed cat states the complex phase of ^^ and ^^ must follow: ^^ =|^^|^^^^ ^^^^ = ^^ ^^2 ^^ ^^The dissipative squeezing stabilization is based on the following jump operator: ^^ ^^ 2 where: ^^2is the two-photon dissipation rate, ^^ is the photon annihilation operator, ^^ is a complex number defining the cat qubit, ^^ and ^^ are the modulus and argument of the complex squeezing parameter ^^ = ^^ ^^^^ ^^. This jump operator can be realized by coupling a lossy buffer mode b with dissipation rate ^^^^, and a four-wave mixing device – typically a Josephson junction or an ATS – to the cat qubit mode and by engineering the Hamiltonian ^^ + h. c., where ^^ is the photon annihilation operator of mode b. While the dissipation of the cat qubit is of the form ^^[ ^^2− ^^2], that of the theoreticalsqueezed cat qubit is of the following more complicated form ^^[cosh( ^^)2 ^^2 +sinh( ^^)2 ^^†2+ cosh( ^^) sinh ( ^^)(2 ^^† ^^ + 1) − ^^2], where ξ and ^^ are assumed to be realnumbers without loss of generality. More particularly, the dissipative squeezingstabilization is engineered through the interaction term [cosh( ^^)2 ^^2 +cosh( ^^) sinh( ^^) (2 ^^† ^^ + 1) + sinh( ^^)2 ^^†2− ^^2] ^^† + h. c.Consequently, the theoretical implementation of the stabilization of the squeezed cat qubit relies on generating four particular terms: i) ^^2^^†+ h. c.: this first term corresponds to the non-linear conversion between two photons of a first mode a – the memory – and one photon of a second mode b – the buffer. This term requires a pump at a frequency|2 ^^^^− ^^^^|, where ^^^^is the resonant frequency of the first mode a and ^^^^is the resonant frequency of the second mode b. ii) ^^†^^ ^^†+ h. c.: this second term corresponds to a longitudinal coupling between the first mode a and the second mode b. This term can be obtained with a pump at a frequency ^^^^. iii) ^^†2^^†+ h. c.: this third term requires a pump at a frequency 2 ^^^^+ ^^^^. iv) ^^2^^†+ h. c.: this fourth term requires the drive of the second mode b at a frequency ^^^^. It must be noted that a so-called “standard cat qubit” is stabilized by only generating the first term ^^2^^†+ h. c. and the fourth term ^^2^^†+ h. c.. However, a considerable obstacle to the use of squeezed cat qubits is the difficulty of engineering such a dissipative squeezing stabilization. This obstacle is above all due to the multiplicity of terms to be engineered using pumps, which jeopardizes the stability of the system, increases thermal population or reduces coherence times. In addition, it must be noted that the generation of the third term requires a very high frequency 2 ^^^^+ ^^^^, around 15 gigahertz (GHz) assuming that the frequencies ^^^^and ^^^^are each typically of 5 gigahertz (GHz). The amount of energy to be injected into the system to reach such a frequency has detrimental effects on the coherence of the system. Moreover, the phase and amplitude of the pump at frequency 2 ^^^^+ ^^^^– for generating the third term – must be precisely tuned with respect to the phase and amplitude of the pump at frequency ^^^^– for generating the second term ^^†^^ ^^†– and the drive of the second mode b at frequency ^^^^– for generating the fourth term ^^2^^†. Therefore, it appears that squeezed cat qubits, although theoretically offering a higher noise bias than cat qubits, are currently not possible to stabilize due to the multiple difficulties of implementing the required dissipation. On the other hand, cat qubits are possible to stabilize but impose extreme conditions of use to achieve a sufficient noise bias and thus run a large-scale quantum algorithm. This fact was the subject of a recent publication by E. Gouzien et al. (2023): "Computing 256-bit Elliptic Curve Logarithm in 9 Hours with 126133 Cat Qubits" (arXiv:2502.06639). B. The moon cat qubit Figure 1 illustrates a schematic diagram of a quantum system 1 arranged to stabilize a particular type of bosonic qubit: a moon cat qubit. The quantum system 1 comprises a non-linear superconducting quantum circuit 3 and a command circuit 5. The non-linear superconducting quantum circuit 3 is arranged to make possible four-wave mixing between a first mode a and a second mode b. In the following, the first mode a hosts a moon cat qubit, while the second mode b is dissipative and is used as a buffer in between the moon cat qubit and the external environment. In the remainder of the description, the first mode a is referred to as the memory mode, while the second mode b is referred to as the buffer mode. The memory mode a and the buffer mode b correspond to natural resonant frequencies of the non-linear superconducting quantum circuit 3. Consequently, the memory mode a and the buffer mode b each have a respective resonant frequency. The memory mode a has a resonant frequency ^^^^^=^ ^^ ^^ and the buffer mode b has a resonant frequency ^^^^=2 ^^, where ^^^^and ^^^^are the respective angular frequencies of the memory mode a and the buffer mode b. By “having” a memory mode and a buffer mode, it should be understood here that the non-linear superconducting quantum circuit 3 comprises components operating in a superconducting regime which host the modes independently of each other or concurrently. In other words, the memory mode a and the buffer mode b may be hosted in different subsets of components of the superconducting circuit or on the same subset of components. The memory mode a has a high-quality factor ^^^^while the buffer mode b has a low- quality factor ^^^^. As will be appreciated, the quality factor (Q-factor) can be determined in various ways, for instance by: (a) spectroscopic linewidth measurement, wherein the quality factor is given by ^^ = ^^ / Δ ^^, wherein ^^ is resonant frequency and Δ ^^ is the spectroscopic linewidth; or (b) via a time-domain measurement wherein a tone is sent in, and a return signal is measured after a pre-determined time, wherein ^^ = ^^ ∗ ^^ with ^^ the characteristic decay time. Of course, the skilled person would be aware of various other methods for determining the quality factor of a particular mode. The non-linear superconducting quantum circuit 3 is intended to be subject to microwave radiations delivered by the command circuit 5 in order to engineer various non-linear interaction between the memory mode a and the buffer mode b. The frequency of each microwave radiation is tuned to select specific terms within the rotating-wave approximation. As detailed below, the non-linear superconducting quantum circuit 3 is intended to interact with the command circuit 5 for engineering a particular Hamiltonian H yielding, through dissipation of the buffer mode b, an effective dissipator which stabilizes a two- dimensional manifold hosting the moon cat qubit. Such a Hamiltonian H can be expressed in the following form: ^^2^^ ^^ − Each term of the Hamiltonian H will be explained and detailed below. The non-linear superconducting quantum circuit 3 comprises a four-wave mixing non- linear element 7 and at least one resonant portion 9. The four-wave mixing non-linear element 7 is arranged to be parametrically driven by the command circuit 5 to cause the non-linear superconducting quantum circuit 3 to engineer the required Hamiltonian H for stabilizing the moon cat qubit. More particularly, the four-wave mixing non-linear element 7 is an asymmetrically threaded superconducting quantum interference device – or ATS hereinafter. It is known to the skilled person that an ATS can be used to engineer the 2-to-1 photon conversion, i.e. the first term ^^2^^†+ h. c. of the above-mentioned Hamiltonian H, in order to perform a dissipative stabilization, as successfully shown by R. Lescanne et al. (2020). Contrary to the first implementations of this stabilization scheme proposed by Z. Leghtas et al. in the article “Confining the state of light to a quantum manifold by engineered two- photon loss” (Science, Vol. 347, No. 6224, 2015) and S. Touzard et al. in the article “Coherent Oscillations inside a Quantum Manifold Stabilized by Dissipation” (Physical Review X 8, 023005, 2018) in which the superconducting circuit element used as a four- wave mixer was a transmon with a single Josephson junction, the solution developed by R. Lescanne et al. (2020) exploits the ATS design which has much lower cross-Kerr terms than the transmon and thus allowed to observe the exponential suppression of bit-flips. The ATS has flux lines through which radiation can be delivered to modulate a common flux and / or a differential flux. It is known to the skilled person that, when biased at its flux working point (0 − ^^ or conversely), the Hamiltonian ^^^^ ^^ ^^of the ATS 7 has the following “sin-sin” form: ^^^ where: - ^^ = ^^^^(^^ + ^^†)+ ^^^^(^^ + ^^†)is the total superconducting phase difference across the ATS 7, - ^^^^is the zero-point fluctuation of the phase of the memory mode a across the ATS 7, - ^^^^is the zero-point fluctuation of the phase of the buffer mode b across the ATS 7, - ^^^^is the Josephson energy of the side junctions, - ^^^^is the inductive energy of the central inductance, - ^^ corresponds to the common flux modulation of the two loops of the ATS 7, - corresponds to the differential flux modulation of the two loops of the ATS 7, and - ℏ is the reduced Planck constant. The common flux modulation ^^ can be implemented by delivering microwave radiations through flux lines 7 out of phase, while the differential flux modulation ^^Δ( ^^) can be implemented by delivering microwave radiations through flux lines of the ATS 7 in phase. Parametric pumping of the ATS 7 is typically done by pumping the common flux as pumping the differential flux merely displaces the modes coupled to the ATS 7. As an example, by pumping the common flux at the frequency ^^^^= |2 ^^^^− ^^ ^^ = ^^^^cos (2 ^^ ^^^^^^), the non-linear resonant part of the Hamiltonian writes in the rotating frame: ^^ ^^ ^^ ^^ which is typically the 2-to-1 photon exchange Hamiltonian needed to engineer the two- photon stabilization. To engineer a longitudinal Hamiltonian between the memory mode a and the buffer mode b, the common flux must be pumped at the frequency ^^^^: ^^^^( ^^) = ^^^^cos (2 ^^ ^^^^^^) In the rotating frame, the parametric part of the Hamiltonian writes: 3 This Hamiltonian can be written as follows to highlight the desired dynamics: ^^ with ^^2^^ = ^^^^^^^^^^^^^^^2^ / ℏ. The first term corresponds to the desired longitudinal coupling. The second term is an effective drive of the buffer mode b at the resonant frequency ^^^^. This term can be used to engineer the term ^^2^^†of the Hamiltonian H. However, since this term needs to be tuned independently of ^^, the buffer drive 15 – or direct drive – is required. The third term corresponds to a spurious and possibly detrimental non-linearity. Depending on the noise level, the moon cat qubit stabilization has to deal with, this can lead to a constraint on the ratio ^^^^ / ^^^^to ensure this last term is not too detrimental. The 2-to-1 photon conversion between the memory mode a and the buffer mode b is obtained by parametrically pumping the ATS at the frequency ^^^^=|2 ^^^^− ^^^^|. Advantageously, the pump frequency satisfies ^^^^≫ ^^2to make this parametric pumping work as well as possible, where ^^2is the two-photon coupling rate. The particular situation where 2 ^^^^= ^^^^( ^^^^= 0) is called the resonant case and may be advantageous, as developed in the European patent application EP 21306965.1 filed by the Applicant. However, in such a resonant case, the 2-to-1 photon conversion cannot be activated by parametrically pumping a four-wave mixing non-linear element such as ATS. Instead, a three-wave mixing non-linear element should be used. As we will see later, the longitudinal term of the Hamiltonian H still requires four-wave mixing to be engineered. Hence, the resonant situation, although compatible with the proposed moon cat qubit stabilization requires an extra three-wave mixing non-linear element. In other words, the resonant case requires having both three-wave mixing and four-wave mixing to generate the first term ^^2^^†+ h. c. and the second term ^^†^^ ^^†+ h. c., respectively. The three-wave mixing can arise by operating the ATS at a different flux point than the 0 − ^^ flux point it is usually operated at or by adding another non-linear element as described in the above-mentioned European patent application EP 21306965.1. In the present invention, the Applicant proposes to use an ATS to stabilize the moon cat qubit hence the advantage of verifying|2 ^^^^− ^^^^|≫ ^^2or being outside the resonant regime for the 2-to-1 photon conversion. The resonant portion 9 is arranged to be coupled or connected to the ATS 7 to provide the non-linear superconducting quantum circuit 3 with the memory mode a and the buffer mode b having respective resonant frequencies ^^^^and ^^^^. More particularly, the memory mode a and the buffer mode b “participate” in the ATS 7, which means that a portion or the entirety of the mode magnetic energy is stored in the ATS 7. Such a participation can be quantified by a zero-point fluctuation of the superconducting phase across the ATS, noted ^^^^for the memory mode a and ^^^^for the buffer mode b. In the schematic diagram of the quantum system 1 illustrated in figure 1, the non-linear superconducting quantum system 3 comprises only one resonant portion, i.e. the resonant portion 9. A single resonant portion can be configured to produce both the memory mode a and the buffer mode b. However, the non-linear superconducting quantum system 3 typically comprises two resonant portions to form the memory mode a and the buffer mode b respectively. The command circuit 5 is arranged to deliver microwave radiations. The command circuit 5 is arranged to parametrically drive the ATS 7 to implement the non-linear conversion between two photons of the memory mode a and one photon of the buffer mode b, to parametrically drive the ATS 7 to implement the longitudinal coupling between the memory mode a and the buffer mode b, and to drive the buffer mode b. To this end and as illustrated in figure 1, the command circuit 5 comprises a first microwave source 11, a second microwave source 13 and a third microwave source 15. The first microwave source 11 is arranged to parametrically drive the ATS 7. More particularly, the first microwave source 11 is used to implement a parametric pump delivering radiation at the frequency|2 ^^^^− ^^^^|to the ATS 7. In the remainder of the description, the first microwave source 11 is referred to as the two-photon pump. The second microwave source 13 is also arranged to parametrically drive the ATS 7. More particularly, the second microwave source 13 is used to implement a parametric pump delivering radiation at the frequency ^^^^to the ATS 7. In the remainder of the description, the second microwave source 13 is referred to as the longitudinal pump. Finally, the third microwave source 15 is arranged to drive the buffer mode b by delivering a radiation at the frequency ^^^^to the resonant portion 9. In the remainder of the description, the third microwave source 15 is referred to as the buffer drive. As explained above, the combined effect of the two-photon pump 11 and the buffer drive 15 allows the non-linear superconducting quantum circuit 3, and more exactly the ATS 7, to perform the non-linear conversion between two photons of the memory mode a and one photon of the buffer mode b as well as the drive of the buffer mode b. In view of stabilizing a two-dimensional manifold hosting the moon cat qubit, which is the purpose of the quantum system 1, such a combined effect contribute to the generation of the terms ^^2^^†+ h. c. and ^^2^^†+ h. c. to the Hamiltonian H to be engineered by the non-linear superconducting quantum circuit 3. Moreover, the parametric pump implemented using the longitudinal pump 13 allows a longitudinal coupling between the memory mode a and the buffer mode b.The consequence of such a longitudinal coupling is the generation of the term ^^† ^^ ^^† +h. c. which also contributes to the Hamiltonian H and thus to the stabilization of the two-dimensional manifold hosting the moon cat qubit. In response to the radiations delivered by the command circuit 5, the non-linear superconducting quantum circuit 3 engineers a Hamiltonian H expressed as follows: ^^2^^ ^^ − where: - g2is the two-photon coupling rate – or interaction strength – proportional to the amplitude of the two-photon pump 11, - ^^ is the annihilation operator of the memory mode a, - ^^ is a complex number which phase and amplitude result from the amplitude of the radiation induced by the longitudinal pump 13, - ^^ is a complex number resulting from the buffer drive 15, - ^^ is the annihilation operator of the buffer mode b, - h. c. designates the Hermitian conjugate. In the expression of the Hamiltonian H, it should be noted that the parameters ^^ and ^^ are both defined relative to the two-photon pump 11. In particular, the phase and the amplitude from which the complex number ^^ results are each defined relative to the phase and amplitude of the two-photon pump 11 respectively. In the remainder of the description, the complex number ^^ is referred to as the moon cat parameter. It may be noted that such a Hamiltonian H does not comprise the term As explained previously, such a term requires an additional pump which would increase the complexity and threaten the stability of the system. Moreover, this term implies a very high frequency – corresponding to 2 ^^^^+ ^^^^– and pumps a quantity of energy into the system which affects the coherence of the bosonic qubit. The Hamiltonian H yields through dissipation of the buffer mode b an effective dissipator ^^[ ^^2^^ ^^ − ^^2] which stabilizes the two-dimensional manifold hosting the moon cat qubit. The relative phase of the microwave sources can be adjusted to obtain the desired number of photons in the memory. More particularly, advantageously, the phase difference between the two-photon pump 11 and the longitudinal pump 13 is substantially equal to zero. In order words, the two- photon pump 11 and the longitudinal pump 13 are configured such that their respective phases are substantially equal to each other. Advantageously, the phase difference between the longitudinal pump 13 and the buffer drive 15 is substantially equal to zero. In order words, the longitudinal pump 13 and the buffer drive 15 are configured such that their respective phases are substantially equal to each other. Preferably, the two-photon pump 11, the longitudinal pump 13 and the buffer drive 15 are configured such that their respective phases are substantially equal to each other. By "substantially equal", it should be understood here that, ideally, the respective phases of the concerned microwave source are strictly equal to each other. However, in practice, such equality is difficult to achieve. Typically, the phase of one microwave source deviates by 5% from the phase of the other. Figure 2 illustrates the respective Wigner functions – or Wigner tomographies – of a standard cat qubit and the moon cat qubit stabilized by the quantum system 1. It can be noted that the standard cat qubit and the moon cat qubit stabilized by the quantum system 1 have the same number of photons ^̅^ = 4, but the latter has a bigger separation in phase-space than the standard cat qubit, meaning the two blobs are more separated. This increased separation at the same number of photons explains the increase in the noise bias. It can be seen that each of the two blobs has a crescent moon shape. This is why the bosonic qubit stabilized by the quantum system 1 is called “moon cat qubit”. B.1 Example implementations of the quantum system Figure 3 illustrates more explicitly the schematic diagram of figure 1. The resonant portion 9 takes the form of a linear microwave network connected to the ATS 7 such that, when coupled via a linear coupler 17 to the ATS 7 which acts as an inductive element, the non-linear superconducting quantum circuit 3 has the memory mode a and the buffer mode b at respective resonant frequencies ^^^^and ^^^^which participates in the ATS 7. The two-photon pump 11 and the longitudinal pump 13 are set-up to modulate the common flux in the ATS 7. In the example illustrated in figure 3, the command circuit 5 further comprises a microwave source 19 set-up to modulate the differential flux in the ATS 7. Such a microwave source 19 can be used to drive the buffer mode b by delivering a radiation at the required frequency ^^^^to the resonant portion 9. In the sense of the present invention, the microwave source 19 can thus also be considered as the “third microwave source” or the “buffer drive” instead of the buffer drive 15. It should be noted that the microwave source 19 can also be used to compensate for a parasitic drive of the buffer mode b caused by the longitudinal pump 13. Although the amplitude and phase of such a compensation drive can be computed analytically, it is fine-tuned experimentally. In order to clearly distinguish the roles of the sources in the Hamiltonian, a microwave network 21 which applies the correct phase offset is shown in figure 3. Alternatively, each microwave source can be simply coupled to a single node of the ATS 7 and their relative phase and amplitude can be set so as to achieve the desired flux modulation. In that case, to modulate the common flux, the microwave sources need to address the circuit out of phase and, to modulate the differential flux, the microwave sources need to address the circuit in phase. As explained above, the two-photon pump 11 is arranged to deliver microwave radiation to the ATS 7 at a frequency substantially equal to |2 ^^^^− ^^^^| to cause the non-linear superconducting quantum circuit 3 to perform the 2-to-1 photon conversion between the memory mode a and the buffer mode b, and thus to engineer the term ^^2^^†+ h. c. of the Hamiltonian H. To convert this 2-to-1 photon conversion into two-photon dissipation, the buffer mode b is selectively coupled to a load 23 via a linear coupler 25 and a microwave filter 27 configured as a band pass filter with a frequency ^^^^. Alternatively, the microwave filter 27 may be configured as a band stop filter at a frequency ^^^^, and may be placed in between, on the one hand, the external environment and, on the other hand, the memory mode a and the buffer mode b to isolate the memory mode a and thus prevents the memory mode a from suffering additional losses coming from unwanted coupling to the load 23. Alternatively, it may be configured as a low-pass (respectively high-pass) filter if ^^^^> ^^^^(resp ^^^^> ^^^^). In other embodiments, the microwave filter 27 can be omitted when coupling between the load 23 and substantially only the buffer mode b can be established. As explained above, the memory mode a has a high-quality factor ^^^^while the buffer mode b has a low-quality factor ^^^^. As explained above, the buffer drive 15 is arranged to drive the buffer mode b by delivering microwave radiation to the resonant portion 9 at a frequency substantially equal to the second resonant frequency ^^^^in order to engineer the term ^^2^^†+ h. c. of the Hamiltonian H. Alternatively, the buffer mode b is driven by the microwave source 19 set at frequency ^^^^. In the above, the load 23 can be seen as part of the command circuit 5 of figure 1, while the linear coupler 25 and the microwave filter 27 can be seen as part of the non-linear superconducting quantum circuit 3. As explained above, the longitudinal pump 13 is arranged to deliver microwave radiation to the ATS 7 at a frequency substantially equal to the second resonant frequency ^^^^to implement a longitudinal coupling between the memory mode a and the buffer mode b in order to engineer the term ^^†^^ ^^†+ h. c. of the Hamiltonian H. B.1.1 Galvanic example implementation Figures 4 and 5 illustrate electrical equivalent diagrams of respective embodiments of the non-linear superconducting quantum circuit 3 which takes the form of a galvanic circuit. Such electrical equivalent diagrams are partial since they both only represent the ATS 7 and the resonant portion 9, and thus represent neither the linear coupler 25 nor the microwave filter 27. More particularly, the resonant portion 9 includes a first resonant portion 29 and a second resonant portion 31. The ATS 7 is realized as known in the art, for example in R. Lescanne et al. (2020). The ATS 7 includes a first Josephson junction 33 and a second Josephson junction 35 in parallel, and an inductive element 37 in parallel between them. As a result, the ATS 7 has two connected loops, each loop comprising a Josephson junction 33, 35 in parallel with the shunt inductive element 37. The inductive element 37 can be realized either geometrically or with a chain of junctions. The ATS 7 has both of its loops flux biased in DC and AC. The DC bias sets the working point of the ATS 7. It may be operated near the so-called saddle point, which is a sweet spot in frequency and has small cross-Kerr terms. Both the first resonant portion 29 and the second resonant portion 31 are galvanically coupled to the ATS 7. The first resonant portion 29 confers the memory mode a having the resonant frequency ^^^^to the non-linear superconducting quantum circuit 3, while the second resonant portion 31 confers the buffer mode b having the resonant frequency ^^^^to the non-linear superconducting quantum circuit 3. By "galvanically coupled", it should be understood here that there are short electrically conducting portions which connect the first resonant portion 29 and the second resonant portion 31 to the ATS 7, i.e., a short electrically conducting track or any other mean ensuring a physically continuous conducting junction. The expression "short" means that the electrically conducting track has an impedance which is negligeable compared to the impedance of the ATS 7, the first resonant portion 29 and the second resonant portion 31 at the resonant frequencies ^^^^and ^^^^. These short electrically conducting portions correspond to the linear coupler 17. In the embodiment of figure 4, the first resonant portion 29 includes a capacitive element 39 and an inductive element 41 which are connected in series. Similarly, the second resonant portion 31 includes a capacitive element 43 and an inductive element 45 which are connected in series. In the embodiment of figure 5, the first resonant portion 29 also includes the capacitive element 39 and the inductive element 41. However, in such an embodiment, the capacitive element 39 and the inductive element 41 are connected in parallel. Similarly, the capacitive element 43 and the inductive element 45 of the second resonant portion 31 are connected in parallel. In the respective embodiments of figures 4 and 5, the non-linear superconducting quantum circuit 3 includes two resonant portions. However, as previously explained, the non-linear superconducting quantum circuit 3 may comprise only one resonant portion which is arranged to produce both the memory mode a and the buffer mode b. B.1.2 Capacitive example implementation Figure 6 illustrates an electrical equivalent diagram of an embodiment of the non-linear superconducting quantum circuit 3 which takes the form of a capacitive circuit. Contrary to figures 4 and 5, the external environment is represented. The ATS 7 serves as a central reference point to which the rest of the components are connected. The filtering and external environment of the buffer mode b are also present to this end. The microwave sources have been omitted for simplicity but are arranged in the same way as in figure 3 to run the circuit. In figure 6, the linear microwave network is formed by a capacitive element galvanically coupled to the ATS 7 to form the buffer mode b and a parallel LC resonator strongly capacitively coupled to the ATS 7 to form the memory mode a. The couplings to the ATS 7 correspond to the linear coupler 17. Both the memory mode a and buffer mode b are strongly coupled to the ATS 7 as shown by the zero-point fluctuation of the phase of the two modes in central inductance of the ATS 7. The buffer is coupled to the external environment by a capacitive element which corresponds to the linear coupler 25. B.2 Characterization of the moon cat qubit In the present invention, the non-linear superconducting quantum circuit 3 engineers a Hamiltonian H expressed as = g2(^^2+ ^^ ^^†^^ − ^^2)^^†+ h. c.. The Hamiltonian H yields through dissipation of the buffer mode b the effective dissipator ^^[^^2+ ^^ ^^†^^ − ^^2] which stabilizes the two-dimensional manifold hosting the moon cat qubit. If the respective phases of the second microwave 13 and the buffer drive 15 (or alternatively 19) are properly tuned, ^^ and ^^2possess the same complex phase, and can thus be taken real without loss of generality. Thus, unless otherwise stated, ^^ and ^^ are assumed to be real numbers. In analogy with the even and odd cat states ^^±^ the even, respectively odd, moon cat state of the moon cat qubit hosted by the two- manifold stabilized using the quantum system 1 is a superposition of moon cat states that contain only even, respectively odd, Fock states. The expansion on Fock basis of the even moon cat state | ^^^+^^^^ ^^ ^^and the odd moon cat state | ^^^−^^^^ ^^ ^^is the following (MCS stands for “moon cat state”): with: ^^−1 ^^0= 1 and: with: ^^1= 1 where: ^^+and ^^−are normalization constants. The figure of merit of the implementation of the quantum system 1 is the relative amplitude of the term ^^2^^†+ h. c. and ^^†^^ ^^†+ h. c.. In this regard, if a comparison between the moon cat qubits of the present invention and the theoretical squeezed cat qubits could be established, ^^ could be expressed as a function of the complex squeezing parameter ^^ = ^^ ^^^^ ^^as follows: ^^ = 2tanh ( ^^) Like for the squeezed cat, the bit-flip probability is exponentially decreased by increasing the value of ^^ – or equivalently ^^ for the squeezed cat. B.3 Performances of the moon cat qubit Figure 7 illustrates the bit-flip probability as a function of the number of photons and for different values of the moon cat parameter ^^. More particularly, for each value of the moon cat parameter ^^, there is a solid line which corresponds to the bit-flip rate of the moon cat qubit stabilized by quantum system 1 and a dotted line which corresponds to the bit-flip rate of the comparable squeezed cat qubit. By "comparable", it should be understood here that, for a given value of the moon cat parameter ^^, we consider the squeezed cat qubit obtained with a parameter ^^ = in accordance with the preceding formula. It appears that, for ^^ = 1 and at a number of photons ^̅^ = 4, the bit-flip of the moon cat qubit stabilized by the quantum system 1 is improved by six orders of magnitude compared to the bit-flip of the standard cat qubit at the same number of photons. Figure 7 shows that, in the case in which the value of the moon cat parameter ^^ is less than or equal to 1, the moon cat qubit (solid lines) involves the same improvement in terms of bit-flip as the squeezed cat qubit (dotted lines) or even slightly better. As previously explained, the phase and amplitude of the moon cat parameter ^^ both result from the amplitude of the radiation induced by the longitudinal pump 13. Advantageously, the amplitude of the radiation induced by the longitudinal pump 13 is tuned so that the moon cat parameter ^^ is equal to 1. Beyond this value, the moon cat qubit stabilized by quantum system 1 is affected by a saturation phenomenon and its performance is therefore less good than the theoretical performance of the squeezed cat qubit. The quantum system 1 also shows less non-adiabatic errors during the implementation of the gate. The description below focuses on the Z gate, but it can be directly generalized to the CNOT gate. During the implementation of a Z gate with a known cat qubit, i.e. a cat qubit stabilized by generating only the term ^^2^^†+ h. c. and the term ^^2^^†+ h. c., two kinds of error may occur: - the single photon loss error given by ^̅^ ^^1^^, where ^̅^ is the number of photons in the cat states, ^^1is the single photon loss and ^^ is the duration of the gate, - the non-adiabatic errors due to the fact that the gate is performed in a finite time ^^2two-photon exchange rate. The non-adiabatic errors prevent to gate as fast as possible in order to avoid the increasing of the single photon loss. A compromise has to be found between the single photon loss error and the non-adiabatic errors since the former is proportional to the duration of the gate, whereas the latter is proportional to the inverse of the duration. The quantum system 1 has less non-adiabatic errors than the known cat qubit. The non- adiabatic errors of the quantum system 1 during a Z gate are given by the following formula: Figure 8 shows simulations of the non-adiabatic errors of the quantum system 1 during a Z gate as a function of the gate time in unit of 1 / ^^2and for a fixed number of photons ^̅^ = 4. It can be seen that, when the value of the moon cat parameter ^^ increases, the non- adiabatic errors are reduced. The solid lines show numerical simulations and the dotted lines the formula above-mentioned. Figure 9 illustrates a comparison of the non-adiabatic errors between the moon cat qubit stabilized by the quantum system 1 (solid lines) and the squeezed cat qubit (dotted lines), and shows that the quantum system 1 performs as well as the squeezed cat qubit. The way the moon cat qubit is stabilized by the quantum system 1 implies that a spurious term of the form ^^†^^( ^^†+ ^^) is activated by the longitudinal pump 13, where ^^ is the photon annihilation operator of the buffer mode b, that is dissipative. However, this spurious term is not detrimental to the quantum system 1 if its amplitude is small enough. ^^2 As mentioned above, its amplitude is proportional to ℏ ^^^^2^^ 2 ^^^2^, where ^^^^and ^^^^are the zero-point fluctuations of the superconducting phase of the memory mode a and the buffer mode b respectively across the ATS 7, and where ^^2^^ is the amplitude of the sought-for longitudinal coupling. Finally, figure 10 illustrates the bit-flip rate of the moon cat qubit stabilized by the quantum system 1 as a function of the ratio^^ ^^^^^^with a photon number approximately given by ^̅^^^= 4 in the memory mode a. The different approximate values of ^̅^^^in this figure are due to the fact that, in practice, this value is difficult to fix. This figure clearly shows that, as long as^^ ^^^^^^< 3, the moon cat qubit stabilized by the quantum system 1 preserves an advantage over the standard cat qubit. More particularly, the curve corresponding to the standard cat qubit is that obtained for ^^ = 0. It should be noted that this threshold value for^^ ^^^^^^depends in a non-trivial manner on the noise parameters of the quantum system 1 and on the ratio^^2^^^^. In figure 10, the curves correspond to the following experimentally relevant parameters: - a ratio^^2^^^^= 0.1, - the photon loss rate of the buffer mode b^^ ^^2 ^^ = 60 megahertz (MHz), - the photon loss rate of the memory mode a^^ ^^ = 0.05 megahertz (MHz), - a dephasing rate of the memory mode a^^ ^^2 ^^ = 0.01 megahertz (MHz), and - a thermal population of the memory mode a ^^^^ℎ= 1%.

Claims

Claims

1. A quantum system (1) for stabilizing a bosonic qubit comprising: - a command circuit (5) including a first microwave source (11), a second microwave source (13) and a third microwave source (15, 19) each arranged for delivering microwave radiation, and - a non-linear superconducting quantum circuit (3) including an asymmetrically threaded superconducting quantum interference device (7) and at least one resonant portion (9) to which said asymmetrically threaded superconducting quantum interference device (7) is connected, said asymmetrically threaded superconducting quantum interference device (7) having flux lines through which radiation can be delivered to modulate a common flux and a differential flux, said non-linear superconducting quantum circuit (3) having a first mode (a) with a first resonant frequency and a second mode (b) with a second resonant frequency, said second resonant frequency being different from twice the first resonant frequency, said second mode (b) being dissipative, the asymmetrically threaded superconducting quantum interference device (7) being arranged such that, when ^ the first microwave source (11) delivers radiation through said flux lines to modulate the common flux at a frequency equal to the absolute value of the difference between twice the first resonant frequency and the second resonant frequency, ^ the second microwave source (13) delivers radiation through said flux lines to modulate the common flux at a frequency equal to the second resonant frequency, and ^ the third microwave source (15, 19) delivers radiation at a frequency equal to the second resonant frequency to the a least one resonant portion (9) to drive the second mode (b), the non-linear superconducting quantum circuit (3) engineers a Hamiltonian H expressed as − ^^2)^^†+ h. c., where g2is proportional to the amplitude ofthe source (11), ^^ is the annihilation operator of the first mode (a), ^^ is a complex number which phase and amplitude result from the amplitude of the second microwave source (13) and are each defined relative to the phase and amplitude of thefirst microwave source (11) respectively, ^^ is a complex number resulting from the third microwave source (15, 19) and defined relative to the first microwave source (11), ^^ is the annihilation operator of the second mode (b), and h. c. is the Hermitian conjugate, said Hamiltonian H yielding through dissipation of the second mode (b) an effective dissipator ^^[ ^^2+ ^^ ^^†^^ − ^^2] which stabilizes a two-dimensional manifold hosting the bosonic qubit in the first mode (a).

2. The quantum system (1) of claim 1, wherein the first microwave source (11) and the second microwave source (13) are configured such that the respective phases of the first microwave source (11) and the second microwave source (13) are substantially equal to each other.

3. The quantum system (1) of claim 1 or 2, wherein the second microwave source (13) and the third microwave source (15, 19) are configured such that the respective phases of the second microwave source (13) and the third microwave source (15, 19) are substantially equal to each other.

4. The quantum system (1) of one of the preceding claims, wherein the first microwave source (11), the second microwave source (13) and the third microwave source (15, 19) are configured such that the respective phases of the first microwave source (11), the second microwave source (13) and the third microwave source (15, 19) are substantially equal to each other.

5. The quantum system (1) of one of the preceding claims, wherein the third microwave source (19) is arranged to deliver radiation through the flux lines to modulate the differential flux at a frequency equal to the second resonant frequency to drive the second mode (b).

6. The quantum system (1) of one of the preceding claims, wherein the second microwave source (13) is arranged to deliver radiation whose amplitude causes ^^ to have an amplitude less than or equal to 1.

7. The quantum system (1) of claim 6, wherein the second microwave source (13) is arranged to deliver radiation whose amplitude causes ^^ to have an amplitude substantially equal to 1.

8. The quantum system (1) of one of the preceding claims, wherein the respective phases of the first mode (a) and the second mode (b) across the asymmetrically threaded superconducting quantum interference device (7) each have a zero-point ^^ fluctuation ^^ , ^^ , and where^^^^ ^^in the ratio ^^^^is less than 3.

9. A method for stabilizing a bosonic qubit performed by the quantum system (1) of one of the preceding claims and comprising the following operations: - delivering, by the first microwave source (11), radiation through the flux lines to modulate the common flux at a frequency equal to the absolute value of the difference between twice the first resonant frequency and the second resonant frequency, - delivering, by the second microwave source (13), radiation through the flux lines to modulate the common flux at a frequency equal to the second resonant frequency, - delivering, by the third microwave source (15, 19), radiation at a frequency equal to the second resonant frequency to the a least one resonant portion (9) to drive the second mode (b), the implementation of said method causing the non-linear superconducting quantumcircuit (3) to engineer a Hamiltonian H expressed as +h. c., where g2 is proportional to the amplitude of the source , isthe annihilation operator of the first mode (a), ^^ is a complex number which phase and amplitude result from the amplitude of the second microwave source (13) and are each defined relative to the phase and amplitude of the first microwave source (11) respectively, ^^ is a complex number resulting from the third microwave source (15) and defined relative to the first microwave source (11), ^^ is the annihilation operator of the second mode (b), and h. c. is the Hermitian conjugate, said Hamiltonian H yielding through dissipation of the second mode (b) an effective dissipator ^^[^^2^^ ^^ − ^^2]which stabilizes a two-dimensional manifold hosting the bosonic qubit in the first mode (a).