## Large gain quantum-limited qubit measurement using a two-mode nonlinear cavity

S. Khan , R. Vijay 1 2 , I. Siddiqi 3 , and A. A. Clerk 1 1 Department of Physics, McGill University, Montreal, Quebec, Canada H3A 2T8 2 Tata Institute of Fundamental Research, Mumbai 400005, India 3 Department of Physics, University of California, Berkeley, CA 94720

(Dated: December 9, 2014)

We provide a thorough theoretical analysis of qubit state measurement in a setup where a driven, parametrically-coupled cavity system is directly coupled to the qubit, with one of the cavities having a weak Kerr nonlinearity. Such a system could be readily realized using circuit QED architectures. We demonstrate that this setup is capable in the standard linear-response regime of both producing a highly amplified output signal while at the same time achieving near quantum-limited performance: the measurement backaction on the qubit is near the minimal amount required by the uncertainty principle. This setup thus represents a promising route for performing efficient large-gain qubit measurement that is completely on-chip, and that does not rely on the use of circulators or complex non-reciprocal amplifiers.

## I. INTRODUCTION

Among the many attractive aspects of circuit QED approaches to quantum information processing is the ability to use driven microwave cavities for high-fidelity qubit measurements [1-3]. The standard approach is to couple the qubit dispersively to the cavity, such that the resonant frequency of the cavity depends on the state of the qubit. By driving the cavity, the state of the qubit becomes encoded in the phase of the reflected drive tone. For an ideal realization, this initial measurement is quantum limited: the backaction disturbance to the qubit is as small as is permitted by quantum mechanics [1, 4]. In practice though, the qubit signal in the cavity output is too small to be sent directly to a room-temperature amplifier. One thus often first uses a low temperature amplification stage, based on a Josephson parametric amplifier (JPA) [5-13]. This amplifier is yet another driven microwave cavity, one that is made nonlinear by the addition of Josephson junctions. While such setups can yield extremely fast measurement rates, they unfortunately require the use of lossy magnetic circulators. This both adds to experimental complexity and bulkiness, and makes one susceptible to insertion losses.

A promising alternate approach is to combine the measurement and amplification cavities in the standard setup into a single system. One thus has an intrinsically nonlinear cavity which is directly coupled to a qubit, eliminating the need for circulators. Such approaches have been studied in regimes where the nonlinear cavity is driven past a bifurcation point in its dynamics, yielding the Josephson bifurcation amplifier [14-17], and a strong latching-style measurement of the qubit state. To obtain a continuous measurement, one could drive the nonlinear cavity close to (but not past) a point of bifurcation, similar to the driving conditions used to realize a paramp. Setups of this kind have been studied both theoretically [1822] and in experiment [18, 23, 24]. While the nonlinearity in these systems allows for large gain, this comes at the expense of an unavoidable large excess backaction on the qubit [19]. While this problem can be mitigated in more complex regimes having a stronger qubit-cavity coupling [21], for a single nonlinear cavity weakly coupled to a qubit, this excess backaction prevents quantum-limited performance.

FIG. 1. (a) Schematic diagram of the PARI setup. The linear cavity is used to provide tunable coupling between the qubit ˆ σ z and the nonlinear cavity; it can be adjusted to suppress excess backaction, while still allowing the output homodyne current ˆ I to see the amplified quadrature ˆ X e (cf. Eq. (30a)) and therefore the large gain G (cf. Eq. (32)). (b) Possible experimental realisation of the PARI setup in (a) in a superconducting circuit architecture. Both cavities are driven, at frequencies ω dj for j = (1 2) , , and λ 12 is the intercavity coupling strength. All components shown above can be cofabricated on-chip as an integrated qubit-amplifier system.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-b76896b3-4ec6-44da-a90a-f33feb6d84d8.jpg)

In this paper, we consider theoretically a new approach to making continuous qubit measurements. The goal is again to combine the measurement cavity and paramp cavities in a traditional setup into a single, circulator-free device. Unlike the work described above, we now consider directly coupling a qubit to a simple multi-mode structure with intrinsic nonlinearity. Our basic setup is shown in Fig. 1. Similar to a standard dispersive-measurement setup, we have a driven linear cavity dispersively coupled to a qubit, as well as a second driven nonlinear cavity. However, unlike the standard setup, we do not take the output field of the linear cavity and use it to drive (via a circulator) the nonlinear cavity. We instead directly couple these systems, such that the intracavity field of the linear cavity directly modulates the inductance of the nonlinear cavity. We refer to this setup as a 'parametric amplifier with resonant input' (PARI). Note that the amplification properties of a very different kind of two-mode nonlinear cavity were recently studied by Eichler et al. in Ref. 25; unlike our

work, they did not consider a direct coupling to a qubit.

It would appear that the PARI will suffer the same fate encountered when directly coupling a qubit to a single nonlinear cavity: there is nothing to protect the qubit from the large fluctuations associated with the nonlinear cavity, hence there will be backaction on the qubit far in excess of the quantum limit. Protecting the qubit from such fluctuations is precisely the role of the circulator in a standard cavity-circulator-paramp setup; the directional nature of the circulator prevents amplified noisy signals emanating from the amplifier input port from reaching the qubit. Similar protection could be achieved without a circulator, if one had a truly directional amplifier, one that did not produce amplified noise at its input port [26].

In the PARI design, there is no explicitly directional element preventing the qubit from seeing the nonlinear cavity and its large, amplified noise. The only element mitigating this coupling is the linear cavity. As we will see, this is already enough to effectively protect the qubit from excess backaction . There is an emergent directionality here as a result of having a driven system, and having the freedom to control the relative phase between the drive tones applied to each of the two cavities. We stress that the linear cavity does not act as a simple filter, suppressing noise at high/low frequencies. Rather, the driven linear cavity acts to control the phase-sensitive coupling of the qubit to the nonlinear cavity, in a manner that would be impossible if the qubit was coupled directly to the nonlinear cavity (as in, e.g., Ref. [19]).

Of course, simply protecting the qubit from backaction is not enough. We also want the system to truly act as an amplifier, producing a sufficiently large qubit signal as well as a sufficiently short measurement time τ meas = 1 Γ / meas (i.e. how long one must run the measurement before being able to distinguish the qubit ground state from the excited state). Our analysis shows that for experimentally feasible parameters, the PARI system has enough flexibility to allow both these goals to be achieved while at the same time having near-ideal backaction properties. We stress that as it contains no circulators, all components of our detector can be cofabricated on-chip. Hence the PARI setup realizes an integrated qubitamplifier system with a large output signal, fast measurement rate and near-quantum-limited backaction.

The remainder of this paper is organized as follows. In Sec. II, we introduce the model for the PARI system, and outline our basic analysis of the corresponding HeisenbergLangevin equations; of particular note is Sec. II E, where we develop an extremely useful mapping of this two-mode system to a single effective nonlinear cavity. In Sec III, we turn to analyzing and optimizing the properties of the PARI detector, viewed as a qubit state detector. We use linear-response theory to calculate the detector gain and noise properties. Our analysis reveals that optimal performance is possible using realistically weak cavity-cavity couplings, if one designs the system so that the damping rate of the nonlinear cavity is greater than that of the linear cavity, and weak nonlinearities are used. Using realistic experimental parameters, we predict that near quantum-limited performance should be possible with photon number gains ∼ 20 dB (cf. Fig. 3), and with fast measurement times ∼ 100 ns (cf. Fig. 4). The backaction is within a factor of 2 of its quantum limited value:

$$1 / \eta _ { \text{meas} } = \Gamma _ { \varphi } / \Gamma _ { \text{meas} } \simeq 1. 7 5 \text{ \quad \ \ } ( 1 )$$

where Γ ϕ is the rate of backaction-induced qubit dephasing (directly proportional to the backaction noise driving the qubit, cf. Eq. (35)), and Γ meas is the measurement rate (inversely proportional to the detector imprecision noise, cf. Eq. (51)). As discussed explicitly in section IV C, this performance (achieved completely on-chip) is on par with or better than modern circulator-based setups [27-30]. In contrast, if one directly (and weakly) coupled a qubit to a single nonlinear cavity, a gain of 20 dB would imply missing the quantum limit by a large factor of ∼ 50 .

## II. BASIC TWO-CAVITY SYSTEM AND ANALYSIS

## A. Hamiltonian and coupling

The two-cavity detector system we consider consists of a single nonlinear cavity with a Kerr type nonlinearity, coherently coupled to an ordinary linear cavity. Throughout this paper, we label the nonlinear cavity with the index 1 , and the linear cavity with index 2 . The two cavity Hamiltonians are given in terms of the cavity annihilation operators ˆ a j ( j = 1 2 , ) by

$$\underset { \text{pli-} } { n \text{ is} } \quad \hat { H } _ { 1 } = \hat { H } _ { 1 } ^ { s y s } + \hat { H } _ { 1 } ^ { \text{envt} } = \omega _ { c 1 } \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 1 } - \Lambda \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 1 } ^ { \dagger } \hat { a } _ { 1 } \hat { a } _ { 1 } + \hat { H } _ { 1 } ^ { \text{envt} } \underset { \text{(2a)} } { ( 2 a ) }$$

$$\iowarpoonright _ { \alpha _ { i c } } \quad \hat { H } _ { 2 } = \hat { H } _ { 2 } ^ { s y s } + \hat { H } _ { 2 } ^ { \text{env} t } = \omega _ { c 2 } \hat { a } _ { 2 } ^ { \dagger } \hat { a } _ { 2 } + \hat { H } _ { 2 } ^ { \text{env} t } \quad \quad ( 2 b )$$

( glyph[planckover2pi1] = 1 ) where ω cj are the cavity resonance frequencies. For the nonlinear cavity Hamiltonian [cf. Eq. (2a)], the Kerr type nonlinearity has strength Λ glyph[lessmuch] κ 1 , where κ 1 is the nonlinear cavity damping rate; this system could be realised as a SQUID shunted with a capacitor in superconducting circuit architecture [31]. We assume a regime where the phase fluctuations of the SQUID are small enough for only the leading order Kerr nonlinearity to be important; the conditions under which higher order terms in an expansion of the SQUID cosine potential can be dropped are discussed in Sec. IV B. Both cavities are also driven externally, via coupling to external microwave transmission lines. The ˆ H envt j terms describe this driving, and the damping (at rate κ j ) of both cavities in a standard input-output theory approach [4, 32]; we neglect internal cavity losses.

The qubit to be measured is connected to the two-cavity system via a dispersive coupling to the linear cavity [1]; this coupling is described by the Hamiltonian,

$$\hat { H } _ { \sigma _ { z } } = A \hat { \sigma } _ { z } \hat { a } _ { 2 } ^ { \dagger } \hat { a } _ { 2 } \equiv \hat { \sigma } _ { z } \hat { F },$$

where ˆ σ z commutes with the qubit Hamiltonian. The cavity2 resonance frequency is modulated by the qubit state by an amount A , the dispersive coupling strength. We have introduced ˆ F , the effective backaction noise operator. Throughout this paper we will consider values of A small enough that a

linear response analysis of the detector is valid; when considering real experimental parameters, we will explicitly check the validity of this assumption. As is standard with dispersive qubit measurements, we measure the qubit in the energy eigenstate ( ˆ σ z ) basis. Measurements along other axes can be achieved by first performing a qubit rotation, as has been done in countless circuit QED experiments.

Finally, we describe the crucial element of the setup: how the two cavities are coupled to each other to allow the nonlinear cavity to learn about the qubit state. For the twocavity system within superconducting circuit architecture, the qubit state-dependent current in cavity2 generates a flux that threads the SQUID loop in cavity1 , therefore modulating its Josephson energy. This inductive coupling between the cavities takes the nonlinear form [see Appendix A],

$$\hat { H } _ { 1 2 } = \frac { 1 } { 2 } \left [ \lambda _ { 1 2 } ( \hat { a } _ { 2 } + \hat { a } _ { 2 } ^ { \dagger } ) \right ] ( \hat { a } _ { 1 } + \hat { a } _ { 1 } ^ { \dagger } ) ^ { 2 } \quad \quad ( 4 )$$

where λ 12 is the strength of the intercavity coupling. Note that the cavity2 'sum' quadrature effectively modulates the 'spring constant' of cavity1 .

We now make an important choice for the two cavity drive frequencies, ω dj ( j = 1 2) , , taking the linear cavity to be driven at twice the driving frequency of the nonlinear cavity:

$$\omega _ { d 2 } = 2 \omega _ { d 1 } \quad \ \ ( 5 ) \quad \text{an}$$

This choice of frequencies is unusual, in that the nonlinear cavity is effectively being coupled to a higher frequency signal, the field of the linear cavity. In particular, as we will show in Eq. (20), this choice (along with a rotating-wave approximation) leads to an emergence of effective parametric driving and photon hopping terms in the coupling Hamiltonian, as opposed to the usual dispersive ( ˆ ˆ a a † ) terms (which are nonresonant in this case).

The full system Hamiltonian is then given by ˆ H = ˆ H 1 + ˆ H 2 + ˆ H 12 + ˆ H σ z . We emphasise that the results of our analysis hold for a generic Hamiltonian of the form of ˆ H , involving a dispersive coupling between two independent modes. As mentioned earlier, our emphasis will be on weak qubit-detector couplings, allowing the use of quantum linearresponse theory [4]. In what follows, we will thus start by developing a description of the PARI in the absence of any coupling to a qubit (i.e. A → 0 ), and will then use this description combined with linear-response theory to determine its properties as a detector.

## B. Heisenberg-Langevin equations

Treating the cavity driving and dissipation as in standard input-output theory [4, 32], the Heisenberg-Langevin equations for the two cavity modes are ( j = 1 2 , )

$$\frac { d } { d t } \hat { a } _ { j } = - i [ \hat { a } _ { j }, \hat { H } _ { j } ^ { s y s } ] - i [ \hat { a } _ { j }, \hat { H } _ { 1 2 } ] - \frac { \kappa _ { j } } { 2 } \hat { a } _ { j } - \sqrt { \kappa _ { j } } \, \hat { a } _ { j } ^ { \text{in} } ( t ), \, ( 6 )$$

where the driving fields describe both coherent drive tones (amplitudes ¯ a in j ) as well as noise:

$$\hat { a } _ { 1 } ^ { \text{in} } ( t ) & = \bar { a } _ { 1 } ^ { \text{in} } e ^ { - i \omega _ { d 1 } t } + \hat { \xi } _ { 1 } ( t ), \\ \hat { a } _ { 2 } ^ { \text{in} } ( t ) & = \bar { a } _ { 2 } ^ { \text{in} } e ^ { i \delta } e ^ { - i \omega _ { d 2 } t } + \hat { \xi } _ { 2 } ( t ).$$

Here, ω dj = ω cj + ∆ j , with ∆ j denoting the detuning from cavity resonance. As we will see, the relative phase between the cavity drive tones will play an important role in allowing one to tune the system properties. We denote this phase difference δ , and without loss of generality take the ¯ a in j in Eqs. (7) to be real and positive. The operators ˆ ( ) ξ j t describe the classical and quantum vacuum noises that enter each cavity via the drive port; the noise is taken to be delta-correlated vacuum noise in the standard input-output theory approach [4, 32].

## C. Classical Cavity Amplitudes

As a result of the strong coherent drives on both cavities, the average cavity field amplitudes will be non-zero, 〈 ˆ a j 〉 glyph[negationslash] = 0 . Our first task will be to find these amplitudes by solving the classical equations of motion for our system. Consistent with our focus on weak qubit-detector couplings and use of linear response, we will solve for these amplitudes in the absence of any coupling to the qubit, and then use them to transform to a displaced interaction picture (with respect to the cavity drive frequencies):

$$\text{$\mathfrak { a }$} _ { j } \quad \hat { a } _ { j } ( t ) & = e ^ { - i \omega _ { d j } t } \left \{ \alpha _ { j } + e ^ { i \mu _ { j } } \left [ e ^ { i \beta } \hat { d } _ { j } ( t ) \right ] \right \} \quad ( j = 1, 2 ), \\ \text{$\mathfrak { a }$} _ { j } \quad \alpha _ { j } \equiv \langle \hat { a } _ { j } \rangle \equiv e ^ { i \mu _ { j } } \sqrt { \bar { n } _ { j } }.$$

Here ¯ n j are the respective cavity photon numbers due to the classical drives, µ j are the average cavity mode phases, and β is an overall phase that will be chosen later for convenience.

The values of the classical cavity mode amplitudes α j are determined by solving the classical equations of motion for both cavities (in the absence of any qubit coupling). The two coupled nonlinear equations obtained can be reduced to a single equation for the nonlinear cavity mode amplitude [see Appendix B],

$$\inf _ { \mathrm d e _ { \mathrm k } } \alpha _ { 1 } [ - i ( \Delta _ { 1 } + 2 \widetilde { \Lambda } | \alpha _ { 1 } | ^ { 2 } ) + \kappa _ { 1 } / 2 ] - i \eta \alpha _ { 1 } ^ { * } = - \sqrt { \kappa _ { 1 } } \, \bar { a } _ { 1 } ^ { \mathrm i n }. \ ( 9 )$$

We see that the cavity-cavity coupling modifies the classical amplitude equation for α 1 from that of a standard Kerr cavity in two ways. The first is a modification of the bare Kerr constant Λ in Eq. (2a) to the complex valued ˜ Λ ,

$$\widetilde { \Lambda } = \Lambda + \lambda _ { 1 2 } ^ { 2 } \left ( \frac { i \chi _ { 2 } } { 4 } \right ), \text{ \quad \ \ } \text{(10)}$$

where χ 2 ≡ χ 2 [ ω = 0] is the (complex) zero-frequency bare linear cavity susceptibility; for non-zero frequency,

$$\chi _ { 2 } [ \omega ] = \left [ - i ( \omega + \Delta _ { 2 } ) + \kappa _ { 2 } / 2 \right ] ^ { - 1 }. \quad \ \ ( 1 1 )$$

Also new, and more complicated, is the last term on the left hand side in Eq. (9), which involves α ∗ 1 ; this describes an effective parametric driving of cavity1 due to the coupling to cavity2 . The effective parametric drive strength η is given by

$$\eta = \lambda _ { 1 2 } \left ( \sqrt { \kappa _ { 2 } } \, \bar { a } _ { 2 } ^ { \text{in} } \, \chi _ { 2 } \, e ^ { i \delta } \right ). \quad \quad \quad ( 1 2 ) \quad \quad \quad$$

It can be shown [see Appendix B] that Eq. (9) can be written as a quintic polynomial in the nonlinear cavity photon number | α 1 | 2 = ¯ n 1 [33]:

$$\bar { n } _ { 1 } = \left [ \frac { \kappa _ { 1 } ( \bar { a } _ { 1 } ^ { \text{in} } ) ^ { 2 } } { D ^ { 2 } } \right ] & \left \{ \kappa _ { 1 } / 2 + i \eta + i ( \Delta _ { 1 }$$

where

$$D = \left | \kappa _ { 1 } / 2 - i \left ( \Delta _ { 1 } + 2 \widetilde { \Lambda } \bar { n } _ { 1 } \right ) \right | ^ { 2 } - | \eta | ^ { 2 } \quad ( 1 4 ) \quad \text{quadquad } \text{quad }$$

We see that in comparison to a single driven Kerr cavity (where the average photon number is determined from a cubic polynomial), the two-mode structure of the PARI system yields more complex nonlinear physics at even the classical level. In what follows, we will attempt to work in regimes where the cavity-cavity coupling λ 12 is weak enough that classical bifurcation physics described by Eq. (13) is not too different from that of a single driven Kerr cavity. At a point of bifurcation, the driven nonlinear cavity becomes multistable when the cavity photon number ¯ n 1 as a function of drive detuning acquires more than one stable solution. Just before this turning-over point is reached, the slope of ¯ n 1 versus ∆ 1 is nearly vertical, and thus the cavity becomes very sensitive to even small changes in detuning; it is near (but not past) this point of bifurcation that we wish to operate.

## D. Linearization procedure

Once the α j are found, the displacement and interaction picture transformation in Eq. (8) allows obtaining the transformed Heisenberg-Langevin equations of motion from Eqs. (6):

$$\frac { d } { d t } \hat { d } _ { j } = - i [ \hat { d } _ { j }, \hat { H } _ { j } ^ { d } ] - i [ \hat { d } _ { j }, \hat { H } _ { 1 2 } ^ { d } ] - \frac { \kappa _ { j } } { 2 } \hat { d } _ { j } - \sqrt { \kappa _ { j } } \, \hat { \xi } _ { j } ( t ), \, ( 1 5 ) \quad.$$

where ˆ H d j are the respective cavity Hamiltonians in the new frame. We focus on the standard regime of strong driving and weak intrinsic nonlinearity ( Λ glyph[lessmuch] κ 1 ), which allows us to drop terms that are cubic and quartic in ˆ d , d 1 ˆ † 1 ; formally, we are neglecting terms that are suppressed by factors of 1 /n ¯ 1 . The nonlinear cavity Hamiltonian in the displaced interaction picture thus takes the form of a detuned degenerate parametric amplifier:

$$\hat { H } _ { 1 } ^ { d } = - \widetilde { \Delta } _ { 1 } \hat { d } _ { 1 } ^ { \dagger } \hat { d } _ { 1 } - \frac { g _ { 1 } } { 2 } ( e ^ { - i 2 \beta } \hat { d } _ { 1 } ^ { \dagger } \hat { d } _ { 1 } ^ { \dagger } + h. c. ), \quad ( 1 6 ) \quad \text{ term} \quad \text{phase}$$

where the effective detuning ˜ ∆ 1 and parametric interaction strengths g 1 are given by

$$\Delta _ { 1 } = \Delta _ { 1 } + 4 \Lambda \bar { n } _ { 1 } \text{ \quad \ \ } ( 1 7 a )$$

$$g _ { 1 } = 2 \Lambda \bar { n } _ { 1 } \text{ \quad \ \ } ( 1 7 b )$$

The linear cavity Hamiltonian following the displacement transformation takes the form:

$$\hat { H } _ { 2 } ^ { d } = - \Delta _ { 2 } \hat { d } _ { 2 } ^ { \dagger } \hat { d } _ { 2 } + \widetilde { A } \hat { \sigma } _ { z } ( e ^ { i \beta } \hat { d } _ { 2 } + e ^ { - i \beta } \hat { d } _ { 2 } ^ { \dagger } + \frac { 1 } { \sqrt { n } _ { 2 } } \hat { d } _ { 2 } ^ { \dagger } \hat { d } _ { 2 } ) \ ( 1 8 )$$

where

$$\widetilde { A } = \sqrt { \bar { n } _ { 2 } } A \, \quad \, \quad \, \quad \, ( 1 9 )$$

is the dressed dispersive qubit-cavity coupling. The linear term in Eq. (18) is a qubit state-dependent drive on the linear cavity that encodes qubit state information in the linear cavity quadratures. In contrast, the nonlinear term provides a qubitstate dependent shift of the linear cavity drive detuning; this term would be of no consequence without the linear driving terms. Furthermore, we show later that for a strongly detuned linear cavity and a weak dispersive coupling, this shift in linear cavity drive detuning is a relatively small effect. Therefore, in calculating how the linear cavity quadratures learn about the qubit state, the final nonlinear term can be ignored.

However, the backaction due to this nonlinear term can be significant, since we anticipate the linear cavity quadratures to experience some amplification by virtue of the coupling to the nonlinear cavity. Large fluctuations in linear cavity photon number can contribute to a 'quadratic' backaction force noise; we characterize this important effect in detail and show that it can be suppressed using appropriate parameter choices, another key advantage of the PARI.

Finally, consider the form of the cavity-cavity coupling Hamiltonian, ˆ H d 12 in the displaced interaction picture. As already mentioned, we take the drive frequencies to satisfy ω d 2 = 2 ω d 1 , allowing us to make a further rotating wave approximation on this coupling (see discussion after Eq. (5)). Retaining only the leading (quadratic) terms in the displaced, rotating frame as before, we obtain:

$$\hat { H } _ { 1 2 } ^ { d } & = \left [ \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } ( e ^ { i \mu _ { 1 2 } } \hat { d } _ { 1 } ^ { \dagger } \hat { d } _ { 2 } + h. c. ) \\ & \quad + \frac { \widetilde { \lambda } _ { 1 2 } ^ { ( 2 ) } } { 2 } ( e ^ { i \mu _ { 1 2 } } e ^ { - i 2 \beta } ( \hat { d } _ { 1 } ^ { \dagger } ) ^ { 2 } + h. c. ) \right ], \quad ( 2 0 )$$

where the dressed intercavity couplings are defined as:

$$\widetilde { \lambda } _ { 1 2 } ^ { ( j ) } = \sqrt { \bar { n } _ { j } } \lambda _ { 1 2 } \quad ( j = 1, 2 ) \quad \quad ( 2 1 )$$

and

$$\mu _ { 1 2 } = \mu _ { 2 } - 2 \mu _ { 1 }, \text{ \quad \ \ } ( 2 2 )$$

where µ j = arg ( α j ) is the phase of the cavityj classical amplitude (cf. Eq. (8)). Note that this coupling is the only term in our Hamiltonian with an explicit dependence on these phases. Also note that this phase difference µ 12 is sensitive to

the relative phase δ between the two drives, and hence can be controlled in an experiment.

The transformed cavity-cavity coupling Hamiltonian in Eq. (20) has two types of terms. The first describes simple photon hopping between the two cavities; this term will allow the nonlinear cavity to know about the qubit. The second term is not a coupling, but effectively renormalizes the cavity1 Hamiltonian. It describes an effective parametric driving, and thus modifies the strength of the parametric drive term arising from the Kerr nonlinearity in cavity1 . We thus define an effective parametric drive strength g eff as:

$$- \, i g _ { \text{eff} } = ( g _ { 1 } - \widetilde { \lambda } _ { 1 2 } ^ { ( 2 ) } e ^ { i \mu _ { 1 2 } } ) e ^ { - i 2 \beta } \equiv | g _ { \text{eff} } | e ^ { i ( \phi _ { \varepsilon } - 2 \beta ) } \quad ( 2 3 )$$

For simplicity, we now make a choice for the reference phase β , defining β = φ / g 2 + π/ 4 , such that g eff is purely real, g eff = | g eff | [34].

## E. Mapping to a single effective nonlinear cavity

## 1. Eliminating the linear cavity

Our driven two-cavity system has been reduced to a linear system of equations, Eqs. (15). For further analysis, it is useful to exactly eliminate cavity2 from the problem, resulting in a description only involving the nonlinear cavity (albeit with modified properties and driving noises). This will be useful for calculating the gain and output noise of our detector; it also directly mirrors our intuition that the linear cavity acts as a tunable coupler between the qubit and the nonlinear cavity. The elimination procedure is described in Appendix C. Introducing canonical quadratures for the displaced cavity fields via ˆ d j = (ˆ x j + ˆ ) ip j / √ 2 , the Fourier-transformed equations for for the nonlinear cavity take the form:

$$\begin{pmatrix} \hat { x } _ { 1 } [ \omega ] \\ \hat { p } _ { 1 } [ \omega ] \end{pmatrix} = \mathbf M _ { \text{eff} } ^ { - 1 } [ \omega ] \begin{pmatrix} \hat { x } _ { e } ^ { \text{in} } [ \omega ] \\ \hat { p } _ { e } ^ { \text{in} } [ \omega ] \end{pmatrix} \quad \quad ( 2 4 ) \quad \text{and}$$

where the operators ˆ x in e , ˆ p in e describe both the incident noise on the system as well as the signal associated with the qubit (see Eq. (C12) in Appendix C). The matrix M eff [ ω ] is the dressed inverse susceptibility (i.e. inverse Green function) of the nonlinear cavity. It has the form

$$\mathbf M _ { \text{eff} } [ \omega ] \equiv \mathbf M _ { 1 } [ \omega ] - i \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } \mathbf \Sigma [ \omega ]. \quad \ \ ( 2 5 )$$

Here, M 1 is the nonlinear cavity inverse susceptibility in the absence of the photon-hopping term in Eq. (20); it simply corresponds to a detuned degenerate parametric amplifier with parametric interaction strength g eff [cf. Eq. (23)]. In contrast, Σ [ ω ] is an effective self-energy matrix, describing the dynamical modification of the nonlinear cavity's properties due to the photon-hopping coupling in Eq. (20). The full form of these matrices is given in Eqs. (C4), (C8) in Appendix C.

In the analysis to follow, we will focus on the lowfrequency properties of the PARI system. Focusing on ω = 0 in our interaction picture, we find that the dressed cavity1

susceptibility has the general form expected of a parametric amplifier driven by a detuned pump:

$$\mathbf M _ { \text{eff} } [ 0 ] = \begin{bmatrix} g _ { \text{eff} } - \frac { \kappa _ { \text{eff} } [ 0 ] } { 2 } & - \widetilde { \Delta } _ { \text{eff} } [ 0 ] \\ \widetilde { \Delta } _ { \text{eff} } [ 0 ] & - g _ { \text{eff} } - \frac { \kappa _ { \text{eff} } [ 0 ] } { 2 } \end{bmatrix} \quad ( 2 6 )$$

where the effective nonlinear cavity detuning ˜ ∆ eff and damping rate κ eff are given by

$$\widetilde { \Delta } _ { \text{eff} } [ 0 ] & \equiv \widetilde { \Delta } _ { \text{eff} } = \widetilde { \Delta } _ { 1 } - \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } \Delta _ { 2 } | \chi _ { 2 } | ^ { 2 } \\ \kappa _ { \text{eff} } [ 0 ] & \equiv \kappa _ { \text{eff} } = \kappa _ { 1 } + \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } \kappa _ { 2 } | \chi _ { 2 } | ^ { 2 } \quad \quad ( 2 7 )$$

Note that the general structure of M eff remains even if λ 12 = 0 : a single driven Kerr cavity also has linearized dynamics equivalent to a detuned, degenerate parametric amplifier [19]. We have thus mapped our system onto an effective singlecavity system.

It is worth noting that the contributions to M eff from the self-energy (i.e. photon hopping) are ∝ λ 2 12 , whereas the direct parametric driving due to cavity2 yields terms that are first order in λ 12 [cf. Eq. (23)]. As a result, it is this second effect that describes the dominant perturbation of the nonlinear cavity due to the cavity-cavity coupling.

## 2. Amplified and Squeezed Quadratures

We thus see that similar to a single driven Kerr cavity, the PARI system can be mapped onto an effective detuned degenerate parametric amplifier model. To make the amplification generated by such a system clearer, it is useful to introduce two new canonical cavity quadratures [19]. We first introduce the angle θ e to parametrize the relative size of the effective pump detuning to the parametric driving strength [Eqs. (27) and (23) respectively],

$$\sin \theta _ { e } = \Delta _ { \text{eff} } / g _ { \text{eff} }. \text{ \quad \ \ } ( 2 8 )$$

We will see that having amplification necessarily requires | ˜ ∆ eff | &lt; g eff , hence we can take -π/ 2 ≤ θ e ≤ π/ 2 for all regimes of interest. In terms of the displaced nonlinear cavity operators ˆ d , d 1 ˆ † 1 , the rotated quadratures ˆ X ,P e ˆ e are then defined as

$$\hat { X } _ { e } & = \frac { 1 } { \sqrt { 2 } } \left ( e ^ { - i \theta _ { e } / 2 } \hat { d } _ { 1 } + e ^ { i \theta _ { e } / 2 } \hat { d } _ { 1 } ^ { \dagger } \right ) \\ \hat { P } _ { e } & = \frac { - i } { \sqrt { 2 } } \left ( e ^ { - i \theta _ { e } / 2 } \hat { d } _ { 1 } - e ^ { i \theta _ { e } / 2 } \hat { d } _ { 1 } ^ { \dagger } \right )$$

By expressing ( M eff ) -1 in this basis, the solutions of the nonlinear cavity equations of motion (at ω = 0 ) take the form

$$\nu _ { \lambda } ^ { - } \quad \hat { X } _ { e } = - \left [ ( \chi _ { e - } ) \hat { X } _ { e } ^ { \text{in} } - \tan \theta _ { e } ( \chi _ { e - } - \chi _ { e + } ) \hat { P } _ { e } ^ { \text{in} } \right ] \quad ( 3 0 a )$$

$$\hat { P } _ { e } = - ( \chi _ { e + } ) \hat { P } _ { e } ^ { \text{in} }$$

where ˆ X ,P in e ˆ in e are the corresponding quadratures of the effective input fields incident on cavity1 , and where the effective susceptibilities χ e ± are given simply by:

$$\chi _ { e + } [ 0 ] \equiv \chi _ { e + } = \left [ \kappa _ { \text{eff} } / 2 + \sqrt { g _ { \text{eff} } ^ { 2 } - \widetilde { \Delta } _ { \text{eff} } ^ { 2 } } \right ] ^ { - 1 } \quad ( 3 1 a ) \quad \ni \,,$$

$$\chi _ { e - } [ 0 ] \equiv \chi _ { e - } = \left [ \kappa _ { \text{eff} } / 2 - \sqrt { g _ { \text{eff} } ^ { 2 } - \widetilde { \Delta } _ { \text{eff} } ^ { 2 } } \right ] ^ { - 1 } \quad ( 3 1 b ) \quad \stackrel { \times } { | \Xi | } ^ { 2 }$$

Amplification in a detuned DPA model emerges when one tunes parameters so that the susceptibility χ e -diverges. This results both in extremely large fluctuations of the quadrature ˆ X e , and in amplification of signals driving this quadrature. The standard parametric photon number gain G is derived from the system scattering matrix (see e.g. [19]):

$$\mathcal { G } [ \omega ] = \left | 1 - \kappa _ { 1 } \chi _ { e - } [ \omega ] \right | ^ { 2 } \underset { \chi _ { e - } \rightarrow \infty } { \longrightarrow } \kappa _ { 1 } ^ { 2 } \chi _ { e - } ^ { 2 } \quad ( 3 2 ) \underset { \lambda _ { 1 2 } \nu } { \lambda _ { 1 2 } \nu } }$$

We thus see that ˆ X e represents the amplified quadrature of the cavity. In contrast, the orthogonal quadrature ˆ P e is not the squeezed quadrature of the cavity: signals incident in this quadrature can also emerge in ˆ X out e with amplification. As discussed in Ref. 19, the finite effective detuning ˜ ∆ eff leads to the fact that standard amplified and squeezed quadratures are not orthogonal.

For a single driven Kerr cavity, one can achieve a diverging χ e -by driving the system close to a point of bifurcation. In the PARI system, we have seen that the basic classical equations determining the classical amplitudes are modified; we thus first need to understand whether a similar bifurcation (and diverging χ e -) can still be attained, at cavity-cavity coupling strengths λ 12 sufficient to allow efficient qubit measurements.

## F. Bifurcation for weak cavity-cavity couplings

To ensure the PARI also exhibits a simple bifurcation in its classical dynamics (and hence amplification), we will focus on regimes of weak cavity-cavity coupling λ 12 , such that this only weakly perturbs the bifurcation that would occur in the uncoupled nonlinear cavity. Using equation Eq. (9) for the classical cavity1 amplitude, and the modification of the parametric driving strength in Eq. (23), we see that this requires:

$$\lambda _ { 1 2 } \sqrt { \bar { n } _ { 2 } } \ll 2 \Lambda \bar { n } _ { 1 } & & ( 3 3 a ) & & \text{$\text{$this$}$}$$

$$\lambda _ { 1 2 } \ll 2 \sqrt { \Lambda / | \chi _ { 2 } | } \, \text{ \quad \ \ } ( 3 3 b ) \text{ \ \ } t e$$

where ¯ n 2 is the linear cavity photon number. The first condition ensures that the coupling-induced parametric driving term in Eq. (9) ( ∝ η ) is much smaller than the main Kerr term ( ∝ Λ ). Note that this condition needs to be satisfied self-consistently, as ¯ n 1 follows from solutions to the classical equations of motion, and thus depends weakly on λ 12 . The second equation, and the requirement of a strongly detuned linear cavity ∆ 2 glyph[greatermuch] κ 2 , ensures that the imaginary Kerr constant induced by the cavity-cavity coupling (cf. Eq. (10)) is

FIG. 2. Plot of the nonlinear cavity photon number ¯ n 1 , against nonlinear cavity drive detuning ∆ 1 . We take a cavity-cavity coupling λ 12 weak enough (as per Eqs. (33)), to ensure a small perturbation of bifurcation physics; parameter values are listed on the plot. The detuning at which the bifurcation occurs when λ 12 = 0 is indicated by the dashed vertical line. The nonlinear cavity drive is adjusted so that the parametric gain G , is large, thereby pushing the nonlinear cavity close to, but not beyond the bifurcation,. The inset plot indicates G versus ∆ 1 ; the large maximum value indicates the preservation of amplification behaviour of the nonlinear cavity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-ff6de83a-fa18-42d4-a305-825ab77dabe4.jpg)

also much weaker than the main Kerr term Λ . These conditions together define the weak coupling regime in terms of the strength of the cavity-cavity coupling λ 12 , where we expect simple bifurcation physics to be retained.

To demonstrate that these conditions are sufficient, we solve Eq. (9) using an appropriately weak λ 12 , and using realistic values for κ 1 , κ 2 ; for reasons that will become clear in section IV A, we require ∆ 2 ∼ κ 1 glyph[greatermuch] κ 2 . Shown in Fig. 2 are resulting plots of the classical cavity1 intensity ¯ n 1 versus the cavity1 drive detuning ∆ 1 . The steepness of the plotted response curves indicates proximity to a bifurcation, and the possibility of large parametric gain. This is confirmed in the inset, where we plot the corresponding parametric photon number gain G [ ω = 0] for the modified nonlinear cavity.

The upshot of the analysis here is that as long as λ 12 is weak enough to satisfy Eqs. (33), one can still have simple bifurcation in the classical dynamics, and hence large parametric gain. We do see modifications in the position of the bifurcation point and the cavity1 drive strength needed to reach this point; further details of these modifications are given in Appendix B, including a perturbation theory approach to determining how the bifurcation shifts.

While a weak λ 12 ensures simple bifurcation physics, the coupling also needs to be sufficiently strong that the qubit signal in the nonlinear cavity output is suitably large. In Appendix G, we calculate a corresponding minimum useful value of λ 12 . Below this value, the coupling to the nonlinear cavity is so weak that it would be more efficient to use the linear cavity output field to read out the qubit state. We find

$$\lambda _ { 1 2 } ^ { \min } \sim \sqrt { \frac { \kappa _ { 1 } \kappa _ { 2 } } { \bar { n } _ { 1 } \mathcal { G } } }$$

## Conditions to ensure a simple bifurcation

Weak coupling regime

Strongly detuned cavity2

$$\frac { \lim e \sqrt { \frac { \kappa _ { 1 } \kappa _ { 2 } } { n _ { 1 } \mathcal { G } } } \ll \lambda _ { 1 2 } \ll \frac { 2 \Lambda \bar { n } _ { 1 } } { \sqrt { n _ { 2 } } } < 2 \sqrt { \frac { \Lambda } { | x _ { 2 } | } } } } { \Delta _ { 2 } \gg \kappa _ { 2 } }$$

Parameter list

Cavity j : damping rate κ j , drive strength ¯ a in j , drive detuning ∆ j

Nonlinear cavity Kerr constant, Λ

Intercavity coupling, λ 12

Relative drive phase, δ

Dispersive coupling, A

TABLE I. Conditions defining the parameter regimes that must be adhered to for the nonlinear cavity to be modified only weakly by its coupling to the linear cavity - the 'weak coupling' regime. The lower table lists the large number of tunable parameters in the PARI scheme that help optimize performance.

For sufficiently large parametric G , this minimum value is very far from the upper limits established by Eqs. (33). Thus, there is a large range of cavity-cavity coupling strengths where one can retain the simple bifurcation physics of a single nonlinear cavity, while still having enough qubit information transferred to the nonlinear cavity.

## III. ANALYSIS AND OPTIMIZATION OF DETECTOR NOISE AND RESPONSE

In the previous section, we demonstrated that the PARI system exhibits a bifurcation similar to a single, driven Kerr cavity, and that near such a point, it can be mapped on to a degenerate parametric amplifier (driven by a detuned pump). We now reintroduce the coupling to the qubit, and analyze the properties of the system as a continuous qubit detector. The PARI system has a large number of parameters, which we summarize in Table I. Thus, a key issue to address is how to pick these parameters to ensure optimal system performance.

We will start our analysis by addressing what seems to be the largest flaw in the design: the qubit is not sufficiently protected from the large noise associated with the driven nonlinear cavity, and hence will experience backaction noise well in excess of the quantum limit. We will show that this expectation is false: with appropriate parameter tuning, the qubit can be made immune to the large fluctuations in the nonlinear cavity. We will then turn to analyzing the output characteristics of the detector: both the gain (i.e. the size of the qubit signal) and the measurement rate (i.e. the size of this signal relative to the amount of noise in the detector output).

## A. Backaction noise

The basic dispersive coupling between the qubit and cavity2 was given in Eq. (3), where we introduced the backaction force operator ˆ F = Aa a ˆ ˆ † 2 2 . For a weak qubit-detector coupling A , the low frequency fluctuations ˆ F will dephase an ini- tial superposition in the qubit. This dephasing is described by the rate (see, e.g., [1, 4])

$$\Gamma _ { \varphi } = 2 S _ { F F } [ 0 ],$$

where S FF [ ω ] is the noise spectral density of ˆ F , evaluated at zero qubit coupling:

$$S _ { F F } [ \omega ] \equiv \int _ { - \infty } ^ { \infty } d t \, e ^ { i \omega t } \langle \hat { F } ( t ) \hat { F } ( 0 ) \rangle _ { 0 }. \quad \ ( 3 6 )$$

Working in the displaced interaction picture [cf. Eq. (8)], the backaction force operator takes the form:

$$\hat { F } = A \sqrt { \bar { n } _ { 2 } } ( e ^ { i \beta } \hat { d } _ { 2 } + h. c. ) + A \hat { d } _ { 2 } ^ { \dagger } \hat { d } _ { 2 } \equiv \hat { F } _ { \text{L} } + \hat { F } _ { \text{Q} } \ ( 3 7 )$$

We have split the backaction force operator into parts that are linear and quadratic in the ˆ d 2 operators. For a suitably strong cavity2 drive, ¯ n 2 glyph[greatermuch] 1 , the dominant noise contribution will come from ˆ F L ; we thus first study this quantity.

## 1. Linearized backaction and backaction avoidance

Using the linearized Heisenberg-Langevin equations for our system, we can decompose the Fourier-transformed Heisenberg-picture operator ˆ [ F L ω ] as

$$\hat { F } _ { \text{L} } [ \omega ] = \hat { F } _ { \text{L} } ^ { \text{int} } [ \omega ] + \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \left ( f _ { X } [ \omega ] \hat { X } _ { e } [ \omega ] + f _ { P } [ \omega ] \hat { P } _ { e } [ \omega ] \right ) ; \, ( 3 8 )$$

see Eq. (D1) in Appendix D for details. Here, ˆ F int L is the λ 12 -independent contribution to ˆ F L . It describes vacuum noise incident from the cavity2 input port which has not entered cavity1 , and hence has not been amplified by the nonlinear cavity. The remaining terms describe noise contributions emanating from the nonlinear cavity. We have written them in terms of the quadratures of the nonlinear cavity ˆ X e and ˆ P e , as introduced in Sec. II E (cf. Eq. (30a) and (30b)). f X [ ω , f ] P [ ω ] are real coefficients which are independent of the cavity1 susceptibility (cf. Eq. (D1)).

As discussed in Sec. II E, it is only the ˆ X e quadrature which experiences amplification due to the effective parametric interaction; thus, this is the only term in Eq. (38) which describes a diverging, amplified noise contribution near the bifurcation. Eliminating this large noise contribution thus involves tuning parameters such that the coefficient f X [ ω ] = 0 over frequencies of interest. We have in general:

$$\text{lative} \quad f _ { X } [ \omega ] = - \frac { \widetilde { A } } { \sqrt { 2 } } \left [ \chi _ { 2 } ^ { * } [ - \omega ] e ^ { i ( \mu _ { 1 2 } - \nu [ \omega ] ) } + \chi _ { 2 } [ \omega ] e ^ { - i ( \mu _ { 1 2 } - \nu [ \omega ] ) } \right ]$$

Here, the phase µ 12 is the relative phase between the classical cavity amplitudes (cf. Eq. (22)); it determines the effective phase of the linearized cavity-cavity coupling (cf. Eq. (20)). The remaining phase angle ν ω [ ] is

$$\nu [ \omega ] = \phi _ { g } / 2 + \theta _ { e } [ \omega ] / 2 + 3 \pi / 4, \text{ \quad \ \ } ( 4 0 )$$

where φ g is the phase of the effective parametric interaction in Eq. (23), and θ e [ ω ] characterises the rotation which defines the amplified ˆ X e quadrature at each frequency (cf. Eq. (28)). For reasons that will become clear shortly, we focus on the ω = 0 regime. Writing the bare cavity2 susceptibility defined in Eq. (11) as χ 2 [ ω ] = | χ 2 [ ω ] | exp( iφ 2 [ ω ]) , f X [ ω = 0] takes the form

$$f _ { X } [ 0 ] = - \sqrt { 2 } \widetilde { A } \left | \chi _ { 2 } [ 0 ] \right | \, \cos ( \mu _ { 1 2 } - \phi _ { 2 } [ 0 ] - \nu [ 0 ] ). \quad ( 4 1 ) \quad \widetilde { d _ { 2 } ^ { \dagger } } \hat { d _ { 2 } } \right )$$

Therefore, at zero frequency, one can make f X [0] (and hence the largest noise contribution to ˆ F L [ ω ] ) vanish by an appropriate choice of the phase µ 12 ; this phase can in turn be varied by tuning the relative phase δ between the two cavity drives. One needs:

$$\mu _ { 1 2 } - \phi _ { 2 } [ 0 ] - \nu [ 0 ] = ( 2 M + 1 ) \frac { \pi } { 2 }, \ \ M \in \mathbb { Z } \quad ( 4 2 ) \quad \frac { \mu _ { 1 2 } \infty } { 2 }.$$

In this zero-frequency case, when this reduced backaction condition is satisfied, the backaction noise in ˆ F L is

$$S _ { F _ { L } F _ { L } } [ 0 ] \simeq S _ { F _ { L } F _ { L } } ^ { \text{int} } [ 0 ] + S _ { F _ { L } F _ { L } } ^ { N L } [ 0 ], \quad \ \ ( 4 3 ) \quad \text{spec} i \text{v re}$$

$$S _ { F _ { L } F _ { L } } ^ { \text{int} } [ 0 ] = \widetilde { A } ^ { 2 } \kappa _ { 2 } | \chi _ { 2 } | ^ { 2 } & & ( 4 4 ) \quad \text{tribu}$$

$$S _ { F _ { L } F _ { L } } ^ { N L } [ 0 ] = \widetilde { A } ^ { 2 } \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } | \chi _ { 2 } | ^ { 2 } \kappa _ { 1 } \chi _ { e + } ^ { 2 } + \dots, \quad ( 4 5 ) \quad \text{ woul}$$

see Appendix D for an outline of the full calculation. S int F F L L [0] describes the unamplified contribution of vacuum noise incident on cavity2 , while S NL F F L L [0] describes the contribution of all noise processes involving the nonlinear cavity (including correlations between the two terms in Eq. (38)). This latter contribution is solely due to the non-amplified cavity quadrature (it is just a function of non-diverging susceptibility χ e + , which scales as 1 /κ 1 near the bifurcation), and thus the backaction noise and dephasing is independent of the amplified vacuum noise fluctuations in the nonlinear cavity . We see that careful tuning of the drive phases gives the system a kind of directionality: the nonlinear cavity will amplify the qubit signal, but the qubit will not see amplified vacuum noise emanating from the nonlinear cavity. This protection against backaction is impossible in systems where a qubit is directly, dispersively coupled to a nonlinear cavity, as there is no way to tune the cavity quadrature seen by the qubit; this is why these systems have a backaction far in excess of the quantum limit [19].

We stress that with appropriate parameter choices, enforcing Eq. (42) at ω = 0 will be enough to guarantee a noise suppression over a wide bandwidth. To see this, note that all the frequency dependence in this equation is ultimately related to the frequency dependence of χ 2 [ ω ] , the bare cavity2 susceptibility (cf. Eq. (11)). By choosing a cavity2 drive tone that is far detuned from resonance (i.e. ∆ 2 glyph[greatermuch] κ 2 ), χ 2 [ ω ] will effectively become frequency independent over the amplification bandwidth. As such, ensuring Eq. 42 holds at zero frequency ensures that it holds over this entire bandwidth. Fig. 6 in Appendix E indicates this frequency dependent behaviour and confirms that the zero frequency approximation is valid.

## 2. Quadratic backaction

We now turn to the fluctuations of ˆ F Q , the second term in Eq. (37) which is proportional to the photon number ˆ ˆ d d † 2 2 . While this term is smaller by a factor of ¯ n 2 than the leading contribution, the fact that it is quadratic in field operators implies that it can still contribute significant noise near bifurcation. Furthermore, calculating the autocorrelation of ˆ ˆ d d † 2 2 involves a convolution of ˆ d , d 2 ˆ † 2 correlators, which incorporates noise at all frequencies. The calculation of the various noise components of this term is outlined in section D; in the relevant large gain limit, the dominant contributions to the quadratic backaction can be approximately written as:

$$S _ { F _ { Q } F _ { Q } } [ 0 ] \simeq A ^ { 2 } \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) _ { \omega = 0, \ \Delta _ { 2 } } ^ { 4 } \frac { \mathcal { G } [ \omega ] ^ { 2 } } { \kappa _ { 1 } ^ { 2 } } \cdot | \chi _ { 2 } [ \omega ] | ^ { 4 } \cdot \Omega [ \omega ], \ ( 4 6 )$$

where G [ ω ] (parametric photon number gain) and | χ 2 [ ω ] | (cavity2 susceptibility) parametrize the nonlinear and linear cavity responses respectively, and Ω[ ω ] is an effective bandwidth. Near ω = 0 and ω = ∆ 2 , and over a frequency range specified by the bandwidth Ω[ ω ] , the nonlinear and linear cavity responses respectively are strongest, indicating a large contribution to the backaction noise. The G 2 dependence here would also occur in a single nonlinear cavity directly coupled to a qubit. While there is no simple tuning possible in the PARI that lets one completely suppress the amplified noise contribution to ˆ F Q , one can strongly reduce the relative importance of ˆ F Q compared to other terms by increasing the cavity2 photon number: S F Q F Q depends on the bare dispersive coupling A and not the dressed ˜ A (cf. Eq. (19)). In the limit where ¯ n 2 →∞ while the dressed dispersive coupling ˜ A is kept fixed, the quadratic backaction formally goes to zero. The PARI system also has additional tunability compared to a single-cavity system which allows a further suppression of this term; for a full discussion, see Appendix D.

## B. Detector response and amplifier forward gain

Having characterized the backaction noise of the PARI detector, we now examine how information on the qubit state appears at the detector output. Referring back to Fig. 1, recall that information about the qubit signal is obtained by measuring the field leaving the nonlinear cavity; this signal will thus benefit from amplification in this cavity. We consider a homodyne measurement, where one measures the quadrature of this output field containing the maximal amount of information on the qubit state. Details on the basics of homodyne measurement and their description can be found in, e.g., Refs. [4, 35]

We first want to understand the signal produced by the qubit, i.e. the average homodyne current associated with the two qubit σ z eigenstates. Note from Eq. (18) that in our displaced frame, the qubit both acts as an effective linear driving force on cavity2 , as well as shifts its frequency. We focus initially on the first effect, which is enhanced by the average

cavity2 photon number ¯ n 2 . Using the effective nonlinear cavity picture developed in Sec. II E, this qubit-induced drive will act as a driving force on both the ˆ X e and ˆ P e nonlinear cavity quadratures; it thus yields new terms in the driving fields in Eqs. (30a), (30b):

$$\hat { X } _ { e } ^ { \text{in} } & \to \hat { X } _ { e } ^ { \text{in} } - \left ( \sqrt { \frac { 2 } { \kappa _ { \text{eff} } } } \, \widetilde { A } \, \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \, | \chi _ { 2 } | \text{Im} \left [ e ^ { i ( \mu _ { 1 2 } + \phi _ { 2 } - \nu ) } \right ] \right ) \hat { \sigma } _ { z } } \quad \text{details} \quad \text{interres} \\ \hat { P } _ { e } ^ { \text{in} } & \to \hat { P } _ { e } ^ { \text{in} } + \left ( \sqrt { \frac { 2 } { \kappa _ { \text{eff} } } } \, \widetilde { A } \, \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \, | \chi _ { 2 } | \text{Re} \left [ e ^ { i ( \mu _ { 1 2 } + \phi _ { 2 } - \nu ) } \right ] \right ) \hat { \sigma } _ { z } } \quad \text{where} \quad \hat { \lambda } \, \cdot \,.$$

The measured homodyne current ˆ I 1 is just proportional to a quadrature of the output field leaving cavity1 , i.e. ˆ I 1 ∝ ( e iφ h ˆ d 1 out , + h.c. ) . For simplicity, we assume in what follows an optimal choice of quadrature, by appropriate choice of the homodyne angle φ h . Using standard input-output theory, we can easily calculate its stationary, average value. Considering the case where we are near bifurcation (and hence have a large parametric photon number gain G , cf. Eq. (32)), we find to leading order

$$\langle \hat { I } _ { 1 } \rangle & \equiv \chi _ { I F } [ 0 ] \langle \hat { \sigma } _ { z } \rangle, \\ \chi _ { I F } [ 0 ] & \simeq \sqrt { 2 } \, \widetilde { A } \, \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \, | \chi _ { 2 } | \, \sqrt { \mathcal { G } } \, \rho _ { e }, \quad \quad ( 4 8 ) \quad \text{ where }$$

where χ IF [0] is the zero-frequency detector response coefficient or 'forward gain' [4], and the angle-dependent prefactor ρ e is given by

$$\rho _ { e } = \sec \theta _ { e } \sin ( \mu _ { 1 2 } + \phi _ { 2 } - \nu + \theta _ { e } ). \quad ( 4 9 ) \quad \frac { \min } { \text{$fmd} }$$

All quantities are at zero frequency. For details of the full calculation, see Appendix E.

As expected, the response of the detector is enhanced by the large parametric gain G glyph[greatermuch] 1 generated when the PARI system is operated near bifurcation. Further, note that tuning the system to satisfy the reduced backaction condition [cf. Eq. (42)] does not make the forward gain vanish . The reduced backaction condition ensures that the qubit is coupled to the nonamplified quadrature ˆ P e of the nonlinear cavity; it thus efficiently drives the orthogonal amplified quadrature ˆ X e , leading to the large, amplified signal.

We have thus far ignored the second effect of the qubit on the linear cavity, namely a dispersive shift of the linear cavity drive detuning ∆ 2 [cf. Eq. (20)]. It turns out that for the case of a strongly detuned linear cavity and weak dispersive coupling so that A glyph[lessmuch] ∆ 2 (parameters indicated in Table II), ∆ 2 ± A produces less than a 1 % change in ∆ 2 . The resulting effect on the gain is small, and so the effect of the nonlinear term in the qubit-detector coupling can be neglected.

## C. Measurement imprecision noise and measurement rate

We now turn to analyzing the intrinsic noise in the output homodyne current, ¯ S II [ ω ] , in the PARI setup; along with the forward gain, this will determine the imprecision of our measurement and the measurement rate (the rate at which information on the qubit state is acquired). The calculation of this intrinsic zero-frequency noise S II [0] (i.e. calculated in the absence of any qubit coupling) follows straightforwardly from the linearized Heisenberg-Langevin equations for our system; details are presented in Appendix E. In the large gain limit of interest, we find to leading order

$$S _ { I I } [ 0 ] \simeq \frac { 1 } { 2 } \kappa _ { \text{eff} } \sec ^ { 2 } \theta _ { e } \, \mathcal { G },$$

where κ eff is the effective nonlinear cavity damping rate defined in Eq. (27). Apart from this renormalized value of κ 1 (which is due to the cavity-cavity interaction), this result has the same general form and dependence on G as for a single Kerr cavity driven near bifurcation [19].

We can now calculate the measurement time τ meas associated with our setup: how long will it take to resolve the qubit signal above the noise in the detector output. For the weak qubit-coupling, linear-response regime we consider, τ meas has the standard definition [4]

$$\frac { 1 } { \tau _ { \text{meas} } } = \Gamma _ { \text{meas} } = \frac { | \chi _ { I F } [ 0 ] | ^ { 2 } } { 2 S _ { I I } [ 0 ] } \equiv \frac { 1 } { 2 } S _ { z z } ^ { - 1 } \quad ( 5 1 )$$

where S II [0] is the intrinsic zero-frequency noise in the measured homodyne current (i.e. in the absence of any coupling to the qubit). In the last inequality, we have introduced the imprecision noise spectral density S zz [0] , which is just the output noise referred back to the qubit.

In the large gain limit, making use of Eqs. (48) and (50), we find

$$\colon \text{full} \quad \Gamma _ { \text{meas} } \simeq \left [ \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } | \chi _ { 2 } | ^ { 2 } \right ] \times \frac { 2 \widetilde { A } ^ { 2 } } { \kappa _ { \text{eff} } } \sin ^ { 2 } ( \mu _ { 1 2 } + \phi _ { 2 } - \nu + \theta _ { e } ) \ ( 5 2 )$$

The nonlinear cavity amplifies both the signal and the noise in the cavity output; as a result, the rate of information gain, Γ meas , is independent of G . This is identical to the case where one directly couples a qubit to a nonlinear cavity [19]. However, unlike this direct coupling case, there is additional tunability in the PARI setup. In particular, the coupling prefactors in Eq. (52) (i.e factor in square brackets) can be tuned to be larger than unity, even if λ 12 is so weak that it does not perturb the bifurcation physics of the nonlinear cavity.

In addition to having a large measurement rate, we also want the sheer magnitude of the output signal to be large enough that the noise of any following amplification stages is irrelevant. It is useful to quantify the size of the output noise via an effective number of thermal quanta ¯ n II : if the cavity were linear and in thermal equilibrium, how hot would it have to be to produce the same amount of output noise? This leads to the definition S II [0] = κ 1 2 (1 + 2¯ n II ) . As both the forward gain | χ IF [0] | 2 and output noise S II [0] scale as G , we are in good shape: ¯ n II ∝ G , and hence, in principle the magnitude of the output noise and signal can be made large enough to make the noise added in following amplification stages insignificant.

## IV. PARAMETER OPTIMIZATION FOR LARGE GAIN, QUANTUM-LIMITED PERFORMANCE

The previous section suggested that the PARI design should allow one to both produce a large-magnitude qubit signal, as well as a small backaction disturbance of the qubit (one that is immune to the amplified fluctuations of the nonlinear cavity). In this section, we show this indeed possible with realistic parameter choices, and also explain the rationale governing these choices.

We start by recalling that the fundamental Heisenberg constraint on the noise of any linear detector implies that there is a minimum possible value for the detector backaction: the backaction dephasing cannot be any smaller than the measurement rate [4]. We thus define the efficiency ratio η meas which characterizes the size of the actual detector backaction compared to this ideal minimum:

$$\frac { 1 } { \eta _ { \text{meas} } } = \Gamma _ { \varphi } / \Gamma _ { \text{meas} } = 4 S _ { z z } S _ { F F } \geq 1 \quad \ \ ( 5 3 ) \quad \ \ T A \, \mathfrak { m }$$

The quantum limit corresponds to achieving η meas = 1 .

## A. Procedure for selecting optimal parameters

In choosing parameters, the first step is to take the relative phase of the cavity drives to fulfill the reduced backaction condition of Eq. (42), which eliminates the contribution of amplified noise to the linear backaction operator ˆ F L . Once this is done, the remaining contributions to the backaction noise come from Eq. (43) (the remaining fluctuations in ˆ F L ) and Eq.(46) (fluctuations in the quadratic backaction operator ˆ F Q ); we begin by analysing the former.

The non-amplified fluctuations in ˆ F L arise from the intrinsic linear cavity backaction and the residual nonlinear cavity backaction. For signal extraction from the nonlinear cavity, the linear cavity backaction represents excess noise not directly tied to any information gain. Ideally, we require the residual nonlinear cavity backaction to be larger than this intrinsic linear cavity noise, so that more information is being extracted from the nonlinear cavity than is leaking out the linear cavity drive port; from Eq. (45), this occurs when the dressed coupling ˜ λ (1) 12 satisfies the constraint:

$$\left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } > \kappa _ { 1 } \kappa _ { 2 } \, \quad \ \ \ ( 5 4 ) \quad \ e a$$

As ¯ n 1 is itself dependent on λ 12 , this is a self-consistent equation for λ 12 . As we work near bifurcation, ¯ n 1 will be large and scaling roughly as κ / 1 (2 √ 3 Λ) [19]. For given values of κ , κ 1 2 and Λ , Eq. (54) constrains the minimum value of λ 12 . This then seems to imply that we should couple the cavities as strongly as possible, thus making ˜ λ (1) 12 large. However, an analysis of the quadratic backaction (cf. Eq. (46)) immediately reveals a problem: the dominant contributions to S F Q F Q depend quartically on the same dressed coupling, and therefore are enhanced more strongly than S NL F F L L when ˜ λ (1) 12 is increased.

TABLE II. Table indicating realistic parameter choices used in Figs. 3, 4, which satisfy the constraints in Eqs. (56a), (56b) as indicated by the lower tables. The cavity drive strengths ¯ a in j , relative drive phase δ , and the nonlinear cavity detuning ∆ 1 form a parameter space where optimal operating points are found via a numerical optimization routine.

| Parameter                                             | Value (MHz)                                           |
|-------------------------------------------------------|-------------------------------------------------------|
| Nonlinear cavity damping, κ 1                         | 100                                                   |
| Linear cavity damping, κ 2 ,                          | 30                                                    |
| Linear cavity detuning, ∆ 2                           | 170                                                   |
| Intercavity coupling, λ 12                            | 1 . 10                                                |
| Dispersive coupling, A                                | 1                                                     |
| Kerr constant, Λ                                      | 0 . 005                                               |
| Nonlinear cavity detuning, ∆ 1                        | (see caption)                                         |
| Cavity drive strengths, ¯ a in j                      | (see caption)                                         |
| Relative drive phase, δ                               | (see caption)                                         |
| Quantum limited performance constraint                | Quantum limited performance constraint                |
| κ 2 1 √ G < ( ˜ λ (1) 12 ) 2 < ¯ n 2 (∆ 2 ) 2 G 3 / 2 | κ 2 1 √ G < ( ˜ λ (1) 12 ) 2 < ¯ n 2 (∆ 2 ) 2 G 3 / 2 |
| 0 . 10 < 1 . 14 < 4 . 93                              | 0 . 10 < 1 . 14 < 4 . 93                              |
| Weak coupling constraints                             | Weak coupling constraints                             |
| ˜ λ (2) 12 < 2Λ¯ n 1 , 2 √ Λ¯ n 2 | χ 2 |             | ˜ λ (2) 12 < 2Λ¯ n 1 , 2 √ Λ¯ n 2 | χ 2 |             |
| 0 . 45 < 0 . 80 , 0 . 74                              | 0 . 45 < 0 . 80 , 0 . 74                              |

Fortunately, the PARI system provides enough flexibility of parameter choices to overcome this problem. We enforce selfconsistently that both large terms in the quadratic backaction S F Q F Q , near ω = 0 , ∆ 2 , are smaller than the contribution to S F F L L from the nonlinear cavity. The form of each term is characterised in Appendix D; we find that to suppress the quadratic backaction, the following conditions must hold:

$$\begin{array} { c c c c c c c } \quad \cdot & & \quad \\ \quad \text{ity,} & & \quad S ^ { N L } _ { F _ { L } F _ { L } } > S ^ { ( 1 ) } _ { F _ { Q } F _ { Q } } & \Longrightarrow & \left ( \widetilde { \lambda } ^ { ( 1 ) } _ { 1 2 } \right ) ^ { 2 } < \frac { \bar { n } _ { 2 } ( \Delta _ { 2 } ) ^ { 2 } } { \mathcal { G } ^ { 3 / 2 } } \\ \quad \text{the} \\ \quad \text{ing} & & \quad S ^ { N L } _ { F _ { L } F _ { L } } > S ^ { ( 2 ) } _ { F _ { Q } F _ { Q } } & \Longrightarrow & \left ( \widetilde { \lambda } ^ { ( 1 ) } _ { 1 2 } \right ) ^ { 2 } < \bar { n } _ { 2 } ( \Delta _ { 2 } ) ^ { 2 } \left ( \frac { \kappa _ { 2 } } { \kappa _ { 1 } } \right ) ^ { 3 } & ( 5 5 ) \\ \quad \dots & & \end{array}$$

where S (1 2) , F Q F Q are the contributions to quadratic backaction near ω = 0 , ∆ 2 respectively (cf. Eq. (46)). These constraints are intuitively clear; large ¯ n 2 reduces the effect of the ˆ ˆ d d † 2 2 operator in the qubit-cavity coupling relative to the linearised backaction operator, and strong linear cavity detuning ensures the two cavity resonances are well-separated in frequency space. Both these conditions reduce the effect of the quadratic backaction, thereby enabling a stronger dressed intercavity coupling ˜ λ (1) 12 . This of course is beneficial for having the nonlinear cavity backaction be the dominant noise source, as described in the discussion following Eq. (54).

The final challenge then is to satisfy the above constraints with experimentally feasible parameter values, allowing close to quantum-limited performance while still being within the weak coupling regime and hence having large gain G . Focusing on realistic parameter choices and requiring weak κ 2 to suppress the imaginary modification of the Kerr constant (cf.

Sec. II F), we enforce κ 2 glyph[similarequal] κ / 1 √ G . This choice makes both conditions in Eq. (55) equivalent; combining the conditions for quantum limited performance in Eqs. (54), (55), and for the weak coupling regime in Eqs. (33a), (33b), we have:

$$\frac { \kappa _ { 1 } ^ { 2 } } { \sqrt { \mathcal { G } } } < \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } < \frac { \bar { n } _ { 2 } ( \Delta _ { 2 } ) ^ { 2 } } { \mathcal { G } ^ { 3 / 2 } } \, \begin{array} { c } \quad & \text{where} \\ \text{$56a$} \\ \text{$\text{$the\colon$} \end{array}.$$

$$\widetilde { \lambda } _ { 1 2 } ^ { ( 2 ) } < 2 \Lambda \bar { n } _ { 1 } \,, \, 2 \sqrt { \frac { \Lambda \bar { n } _ { 2 } } { | \chi _ { 2 } | } } \quad \ ( 5 6 b ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

We reiterate that these conditions are to be satisfied selfconsistently. The lower and upper boundaries of inequality (56a) are where S NL F F L L is comparable to the S int F F L L and S F Q F Q respectively; for ideal performance, we require ˜ λ (1) 12 to be not too close to either boundary, so that S NL F F L L is the dominant noise contribution. We will present an example parameter set that satisfies these constraints in the next section.

To summarize, near quantum-limited PARI performance requires (i) tuning the relative drive phase δ to satisfy the reduced backaction condition of Eq. (42), and (ii) choosing a cavity-cavity coupling λ 12 large enough to ensure that the residual nonlinear cavity backaction S F F L L overcomes the intrinsic linear cavity noise, but also not too large , so that the quadratic backaction remains small in comparison, while ensuring the nonlinear cavity is not too strongly modified. The balance required forces ˜ λ (1 2) , 12 , ˜ λ (2) 12 to satisfy Eqs. (56a), (56b) respectively.

These constraints allow us to see the advantage of weak nonlinearities in this setup. As mentioned earlier, large ¯ n 2 values are beneficial for quantum-limited performance, as they suppress the quadratic backaction operator. However, ¯ n 2 appears in the dominant modification of the nonlinear cavity bifurcation (cf. Eq. (56b)), and cannot be increased indefinitely, unless λ 12 is simultaneously decreased so that ˜ λ (2) 12 is held fixed. This, however, raises another potential problem: weakening λ 12 suppresses the gain and the measurement rate, and reduces S NL F F L L relative to S int F F L L . The way out is to ensure the dressed coupling ˜ λ (1) 12 is held fixed while λ 12 is decreased; reducing the Kerr nonlinearity Λ to increase ¯ n 1 achieves precisely this effect.

## B. Optimal performance with realistic parameter values

We now implement the above guidelines for choosing parameters, focusing on experimentally accessible values that yield near quantum-limited performance with parametric photon number gains G (defined in Eq. (32)) of O (100) , or O (20 dB) ; corresponding results for the gain, measurement rate and noise are shown in Figs. 3 and 4.

As discussed, one ideally needs a large cavity1 photon number ¯ n 1 , which requires a relatively weak cavity1 Kerr constant. We thus take a value Λ = 5 × 10 -5 κ 1 . Note that for standard lumped element resonators based on a shunted

SQUID, the strength of the Kerr nonlinearity is given by

$$\Lambda / \kappa _ { 1 } = \frac { \pi } { 2 } \frac { R } { R _ { K } },$$

where R is the effective shunt resistance , and R K = h/e 2 is the quantum of resistance. For typical setups with R ∼ 50 Ω , Kerr nonlinearities on the order of Λ ∼ 10 -2 κ 1 are standard [31]. Our chosen value of Λ /κ thus requires this nonlinearity to be weakened, or 'diluted'. The utility of such weakened nonlinearities, and methods for achieving this, have recently garnered attention [36]. One could for example introduce an additional linear inductance L 0 in series with the SQUID (Josephson inductance L J ). This effectively reduces Λ /κ by a factor of p 3 , where the participation ratio p = L / L J ( 0 + L J ) can be made small by choosing an appropriately large L 0 (see, e.g. Ref. [36]). Our chosen Λ /κ thus requires p glyph[similarequal] 0 1 . . One must also be careful that the weakened nonlinearity (and consequent larger photon numbers) do not cause higher nonlinearities in the SQUID potential to become relevant. This condition results in the constraint Qp glyph[greaterorsimilar] 5 , where Q is the nonlinear cavity quality factor [3740]. Using the parameters in Table II and taking a cavity1 resonance frequency ω c 1 = 6 GHz yields Q = 60 and p glyph[similarequal] 0 1 . . One thus satisfies the requirement on Qp , while at the same time having a not too-large cavity2 resonance frequency ω c 2 ∼ 2 ω c 1 ∼ 12 GHz. Finally, we note that weak nonlinearities can also be achieved by replacing the single SQUID with an array of SQUIDs, or by using transmissionline resonators (see Ref. [36, 41]).

Moving on to the remaining parameters, recall that optimal performance occurs when κ 2 ∼ κ / 1 √ G ; we take κ 2 = 0 3 . κ 1 . We then implement the optimization conditions. The resulting parameter choices are all shown in Table II. We take values that are experimentally realistic; more experimentallychallenging choices of parameters would lead to even more optimal performance. The nonlinear cavity drive strength ¯ a in 1 and detuning ∆ 1 are chosen to be close enough to the bifurcation to achieve the required parametric gain. A first order perturbation approach narrows the region in phase space where this bifurcation exists; a numerical optimization routine in ¯ a in 1 -∆ 1 -δ space is then implemented to tune these parameters to obtain closest-to-desired detector performance.

In Fig. 3, we plot both the parametric gain G (dashed blue), as well as the inverse efficiency ratio, 1 /η meas (solid red), both as a function of the nonlinear cavity detuning ∆ 1 . Wefind that as predicted, one can both achieve a large effective parametric gain of ∼ 20 dB (implying a large output signal), and have a backaction noise which is approximately a factor of 1 75 . more than the minimal quantum limited value. The corresponding efficiency is η meas glyph[similarequal] 0 57 . . Had one directly coupled the qubit to the nonlinear cavity, the backaction would exceed the quantum limit by a large factor ∼ 100 .

In Fig. 4 we plot the measurement rate Γ meas versus ∆ 1 for the same parameter set. We focus on the weak dispersive coupling regime, and take A = 0 01 . κ 1 . It is possible to obtain measurement rates on the order of 10 MHz with a nearly quantum-limited detector.

We note that the ability to reach the quantum limit on mea-

FIG. 3. Plot of the zero frequency parametric gain G [cf. Eq. (32)] (dashed blue, left-hand axis) and the inverse efficiency ratio 1 /η meas (cf. Eq. (53)) (solid red, right hand log axis), as a function of the nonlinear cavity detuning ∆ 1 . For optimal ∆ 1 glyph[similarequal] -0 83 . κ 1 , we obtain 1 /η meas glyph[similarequal] 1 75 . , with G glyph[similarequal] 100 . The drive parameters are all fixed for this plot: (¯ a in 1 ) 2 /κ 1 = (62 2) . 2 , (¯ a in 2 ) 2 /κ 1 = (45 5) . 2 , δ = 0 40 . , ¯ n 1 ∼ 8000 , ¯ n 2 ∼ 1200 . The dashed green curve with intermediate minimal backaction-imprecision product corresponds to the limit where ¯ n 2 → ∞ while ˜ A and ˜ λ (2) 12 are kept fixed; in this limit, the quadratic backaction S F Q F Q (cf. Eq. (46)) does not contribute. The dot-dashed black curve is the inverse efficiency where we also take κ 2 → 0 , thus suppressing the intrinsic linear cavity backaction S int FF (cf. Eq. (45)). In this limit, one can exactly reach the quantum limit.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-a1cd39c8-ae57-4fda-b621-38362ad95322.jpg)

surement is restricted by the non-zero intrinsic linear cavity backaction S int F F L L and the second order backaction S F Q F Q . To see this, we include in Fig. 3 the dashed green curve showing the inverse efficiency ratio 1 /η meas with S F Q F Q taken to zero ( ¯ n 2 → ∞ , ˜ A held finite). Also included is the dot-dashed black curve showing this product with S int F F L L also excluded ( κ 2 → 0 ), indicating precisely quantum-limited measurement if both noise terms can be suppressed. The relative contributions of each term depends on which boundary of Eq. (56a) the dressed coupling ˜ λ (1) 12 is closer to.

As mentioned earlier, weak Kerr nonlinearities are advantageous for optimal PARI performance. Using a stronger value of Λ limits the cavity intensities ¯ n j , hence making the quadratic backaction more important. However, the multiple parameters in the PARI provide a way out; adjusting the linear cavity detuning ∆ 2 allows us to suppress S F Q F Q , albeit at the cost of a reduced measurement rate. In Fig. 8 in Appendix F, we include an analogous plot to Fig. 3, with a stronger nonlinearity, Λ = 10 -4 κ 1 ; we obtain η meas glyph[similarequal] 0 5 . with Γ meas ∼ 1 MHz .

We conclude this section by mentioning the freedom of parameter choices in the PARI scheme, even under the constraints of Eqs. (56a), (56b). The ability to tune the relative drive phase δ means that there exist multiple parameter sets allowing optimal PARI performance. This is explicitly shown in Figs. 7 (a), (b) in Appendix F, where 1 /η meas and G are analyzed in ∆ 1 -δ space; a manifold of optimal operating parameters can clearly be seen to exist.

FIG. 4. Plot of the zero frequency measurement rate Γ meas (dashed blue, left-hand axis) and the inverse efficiency ratio 1 /η meas (solid red, right hand log axis), as a function of the nonlinear cavity detuning ∆ 1 . For optimal ∆ = 1 -0 831 . , we have 1 /η meas glyph[similarequal] 1 75 . and Γ meas glyph[similarequal] 10 MHz . The drive parameters are the same as for Fig. 3: (¯ a in 1 ) 2 /κ 1 = (62 2) . 2 , (¯ a in 2 ) 2 /κ 1 = (45 5) . 2 , δ = 0 40 . , ¯ n 1 ∼ 8000 , ¯ n 2 ∼ 1200 . We take κ 1 = 100 MHz to compute the measurement rate, with the dispersive coupling strength A = 1 MHz; the choice for A is justified in the next section.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-b1344ce0-7302-4d32-b6a7-67543e09294e.jpg)

## C. Comparison to conventional circulator-based setup

The parametric photon number gain G (cf. Eq. (32)) allows a comparison of the amplification performance of the PARI scheme to that of conventional qubit-detector-amplifier setups. It is also useful to develop a meaningful way to compare the noise performance of both schemes in a common framework.

For the conventional setup, the linear cavity output signal is quantum-limited, but the necessary following amplification stage adds imprecision noise that causes deviations from ideal performance. Fortunately, the qubit is protected from additional backaction due to the use of circulators, and so the backaction noise after amplification, S FF , is unchanged from its quantum-limited value. The amplifier boosts the output signal but also makes it noisier; this added noise then modifies the output signal imprecision noise S zz to κ 1 2 (1 + 2¯ n add ) (for no noise added by the following amplifier, ¯ n add = 0 and S zz is unchanged). If, following amplification, the backactionimprecision product deviates from its quantum-limited value by a factor of x , then attributing this deviation solely to an increase in the imprecision noise, we can write S zz = x κ 1 2 ≡ κ 1 2 (1 + 2¯ n add ) , so that ¯ n add = ( x -1) / 2 . For the parameters used in Figs. 3, 4, we find x glyph[similarequal] 1 75 . , and the PARI noise performance is therefore equivalent to that of a standard qubit-cavity-amplifier setup where the paramp stage adds ¯ n add glyph[similarequal] 0 375 . quanta of noise.

FIG. 5. Plot of 〈 ∆ I 1 〉 (cf. Eq. (59)) against the dispersive coupling strength A . We see that for A ≤ 0 01 . κ 1 , the system response is approximately linear in A and the results of our analysis hold. We choose a value of A = 0 01 . κ 1 = 1 MHzfor Fig. 4 to operate within this linear regime while having a fast measurement rate.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-967ca636-88dd-4fac-963e-3c2ad2f733b5.jpg)

## D. Validity of linear response theory

Our analysis assumes operation of the PARI in the limit of weak dispersive coupling, so that the detector can be treated in the linear response regime. It is important to determine whether the parameter choices used in Figs. 3, 4 in fact allow the detector response to remain linear in the dispersive coupling strength A . A standard approach is to include the exact qubit-cavity dispersive coupling at the level of the classical equations; this coupling leads to a qubit state dependent shift of the linear cavity frequency by ± A , for ˆ σ z = {↑ ↓} , = { 1 , - } 1 respectively. We can then compute a classical cavity 1 output homodyne current that knows about A :

$$\langle \hat { I } _ { 1 \sigma _ { z } } \rangle ( A ) = \frac { \kappa _ { 1 } } { 2 } \left ( e ^ { i \phi _ { h } } \alpha _ { 1 \sigma _ { z } } ( A ) + h. c. \right ) \quad ( 5 8 ) \quad \text{systo}$$

The difference between the currents for the two possible states,

$$\langle \Delta I _ { 1 } \rangle = \langle \hat { I } _ { 1 \uparrow } \rangle - \langle \hat { I } _ { 1 \downarrow } \rangle \quad \quad ( 5 9 ) \quad \text{wh}$$

measures how well the detector output discriminates between these states; it is thus a measure of the detector response as a function of A . We plot 〈 ∆ I 1 〉 against A in Fig. 5 for the same parameters as Figs. 3, 4. For this optimal parameter set, choosing dispersive coupling strengths A glyph[lessorsimilar] 1 MHz (for κ 1 = 100 MHz) allows a linear detector response, and our analysis can be safely applied.

## V. CONCLUSIONS

In this paper, we have detailed the theory of a two-cavity amplifier for qubit state measurement in the linear response regime. The setup, where the linear cavity housing the qubit is coupled directly to the nonlinear cavity used for measurement, can be mapped to a single modified nonlinear cavity coupled to the qubit via a tunable coupling. Adhering to a weak intercavity coupling regime ensures retention of the nonlinear cavity's bifurcation, and hence a large amplifier gain. We find further that the linear cavity's response modifies the amplifier's properties. Crucially, the modified qubit-cavity coupling can be tuned to shield the qubit from linear measurement backaction, while the multiple parameter choices allow suppressing the quadratic backaction, all while preserving the large gain. Therefore, quantum-limited qubit state measurement that escapes a single nonlinear cavity detector within weak coupling is possible with the PARI device.

## VI. ACKNOWLEDGEMENTS

We thank Andrew Eddins, David Toyli, and Eli-Levenson Falk for many useful discussions. This work was supported by NSERC and the Army Research Office under Grant #: W911NF-14-1-0078.

## Appendix A: Derivation of intercavity coupling Hamiltonian

For the setup shown in Fig. 1 (b), the Josephson energy of the nonlinear circuit's SQUID can be written in terms of the flux threading the SQUID loop,

$$\hat { \Phi } _ { \text{ext} } = \Phi _ { B } + M \hat { I } _ { 2 } \text{ \quad \ \ } ( A 1 )$$

where Φ B is a static flux bias. The second term contains the contribution to the flux due to the coupling with the linear cavity, and is proportional to the current I 2 flowing in the linear resonator. M is the mutual inductance of the coupled cavity system. The physical Hamiltonian of the nonlinear cavity in Fig. 1 (b) is therefore:

$$\hat { H } _ { 1 } ^ { s y s } = \frac { \hat { Q } _ { 1 } ^ { 2 } } { 2 C _ { 1 } } - E _ { J } [ \hat { \Phi } _ { \text{ext} } ] \, \cos \left ( 2 \pi \frac { \hat { \Phi } _ { 1 } } { \Phi _ { 0 } } \right ) \quad ( A 2 )$$

where ˆ Φ ˆ , Q are the canonical coordinates describing the nonlinear cavity with shunt capacitance C 1 , and Φ 0 = h/ e 2 is the magnetic flux quantum. The first term includes the capacitive energy of the circuit, while the second is the Josephson energy term providing the nonlinearity and the intercavity coupling. This expression can be simplified via Taylor expansion, keeping terms up to and including ˆ Φ 4 1 , ignoring constant terms and the term proportional only to ˆ I 2 (these are accounted for in the total drive on the linear cavity). As a result, we obtain the quadratic- and quartic-inˆ Φ 1 terms that form the nonlinear cavity Hamiltonian with a Kerr type nonlinearity, as is standard, together with a term depending on both Φ 1 and ˆ I 2 ; we focus on this intercavity coupling term:

$$\hat { H } _ { 1 2 } = \left \{ \frac { M } { 2 } \left ( \frac { 2 \pi } { \Phi _ { 0 } } \right ) ^ { 2 } \frac { d E _ { J } } { d \Phi _ { e x t } } \Big | _ { \Phi _ { B } } \right \} \hat { I } _ { 2 } \hat { \Phi } _ { 1 } ^ { 2 } \quad ( A 3 )$$

To generalize and simplify the analysis, we now introduce the usual creation and annihilation operators for both cavities,

$$\hat { a } _ { j } = & \frac { 1 } { \sqrt { 2 \hbar { L } _ { j } \omega _ { c j } } } \hat { \Phi } _ { j } + i \frac { 1 } { \sqrt { 2 \hbar { C } _ { j } \omega _ { c j } } } \hat { Q } _ { j } \\ \hat { a } _ { j } ^ { \dagger } = & \frac { 1 } { \sqrt { 2 \hbar { L } _ { j } \omega _ { c j } } } \hat { \Phi } _ { j } - i \frac { 1 } { \sqrt { 2 \hbar { C } _ { j } \omega _ { c j } } } \hat { Q } _ { j } \quad \ ( A 4 ) \quad \text{The}$$

where ω cj = 1 / √ L C j j is the natural frequency of cavity j ( j = 1 2 , ). Using the fact that ˆ I 2 = ˆ Φ 2 /L 2 , we can rewrite the intercavity coupling term in Eq. (A3) in terms of creation and annihilation operators to get the form written in Eq. (4). From there, the intercavity coupling strength λ 12 can be extracted in terms of physical quantities,

$$\lambda _ { 1 2 } = 2 \omega _ { c 1 } \sqrt { \$$

To obtain a numerical estimate, we use the explicit form for the Josephson energy E J [Φ] = E 0 ∣ ∣ ∣ cos ( π Φ Φ 0 )∣ ∣ ∣ , and rewrite L 2 in terms of the characteristic impedance Z 2 = √ L /C 2 2 and ω c 2 , hence simplifying the expression for λ 12 ,

$$\lambda _ { 1 2 } & = \left [ M \cdot ( \omega _ {$$

i where R K = h/e 2 = 2Φ 0 /e is the quantum of resistance. Taking ω c 1 = (2 π )6 GHz, ω c 2 = (2 π )12 GHz, choosing Z 2 glyph[similarequal] Z 0 , where Z 0 is the impedance of free space, and noting that the angular factor will be O (1) for the static bias fluxes used in such setups, we obtain,

$$\lambda _ { 1 2 } \sim \left ( \frac { M } { p H } \$$

From the considered parameter choices (cf. Table II), we have a value of λ 12 glyph[similarequal] 1 1 . MHz, which would correspond to a mutual inductance M glyph[similarequal] 11 pH .

## Appendix B: Classical equations of motion

We consider the full Heisenberg equations for the two cavities, Eq. (6), and compute the classical equations of motion assuming a steady state coherent solution of the form indicated in Eq. (8), without the quantum operators. Imposing the driving frequency condition [cf. Eq. (5)] followed by a rotating wave approximation allows one to easily obtain the following classical equations of motion for the nonlinear and linear cavities respectively:

$$\begin{array} {$$

The second equation can be used to eliminate α 2 to obtain a complex nonlinear equation for α 1 ,

$$&, \quad \text{comp.} \lambda \text{numeric equation} \, \in$$

The first line in the equation above is just the classical equation of motion for an individual nonlinear cavity, while the second line represents the modifications due to the coupling to the linear cavity. We can rewrite the modifications proportional to λ 2 12 as changing the Kerr constant,

$$\underset { \text{in} } { \text{m} } \quad \alpha _ { 1 } \left [ - i ( \Delta _ { 1 } + 2 \widetilde { \Lambda } | \alpha _ { 1 } | ^ { 2 } ) + \kappa _ { 1 } / 2 \right ] - i \eta \alpha _ { 1 } ^ { * } = - \sqrt { \kappa _ { 1 } } \, \bar { a } _ { 1 } ^ { \text{in} }$$

where the modified Kerr constant ˜ Λ is as defined in Eq. (10), and the emergent parametric driving strength η is as defined in Eq. (12). Eq. (B2) mixes α 1 and α ∗ 1 , so that the phase of α 1 plays an important role. However, we now note that it is possible to write down the above equation and its complex conjugate as a matrix system:

$$\overline { \mathcal { C } _ { 2 } } \quad \begin{bmatrix} \kappa _ { 1 } / 2 - i \zeta & - i \eta \\ i \eta ^ { * } & \kappa _ { 1 } / 2 + i \zeta ^ { * } \end{bmatrix} \begin{pmatrix} \alpha _ { 1 } \\ \alpha _ { 1 } ^ { * } \end{pmatrix} = - \sqrt { \kappa _ { 1 } } \begin{pmatrix} \bar { a } _ { 1 } ^ { \text{in} } \\ \bar { a } _ { 1 } ^ { \text{in} } \end{pmatrix} \quad \text{(B4)}$$

where ζ is given by

$$\zeta = \Delta _ { 1 } + 2 \widetilde { \Lambda } | \alpha _ { 1 } | ^ { 2 }$$

Inverting the 2-by-2 matrix above allows an expression for α 1 that does not include explicit reference to its phase; squaring this expression allows one to obtain Eq. (13). The equation can be solved numerically for the nonlinear cavity photon number profiles indicated in Fig. 2.

The purely real equation is also more tractable for a perturbative expansion in the weak intercavity coupling strength. We expand the detuning, drive strength, and nonlinear cavity occupation around the unmodified nonlinear cavity bifurcation (∆ 1 0 , , a ¯ in 1 0 , ) (where dn /d ¯ 1 ∆ 1 →∞ ) [19],

$$\lambda _ { 1, 0 / } \, \stackrel { \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ \Delta _ { 1 } & = \Delta _ { 1, 0 } + \frac { \lambda _ { 1 2 } } { \kappa _ { 1 } } \cdot \Delta _ { 1, 1 } \\ \bar { a } _ { 1 } ^ { \text{in} } & = \bar { a } _ { 1, 0 } ^ { \text{in} } + \frac { \lambda _ { 1 2 } } { \kappa _ { 1 } } \cdot \bar { a } _ { 1, 1 } ^ { \text{in} } \\ \bar { n } _ { 1 } & = \bar { n } _ { 1, 0 } + \frac { \lambda _ { 1 2 } } { \kappa _ { 1 } } \cdot \bar { n } _ { 1, 1 } \quad \quad ( B 6 ) \\ \ n \, & \text{the influence use income that } \Lambda \,. \, \bar { a } ^ { \text{in} } \quad \bar { n } \,.$$

To retain the bifurcation, we impose that ∆ 1 1 , , ¯ a in 1 1 , , ¯ n 1 1 , be chosen such that the O λ ( 12 ) modifications to both the classical nonlinear equation and the diverging derivative at bifurcation dn /d ¯ 1 ∆ 1 , vanish. This yields a modified set of detuning and driving strengths describing the new position of the bifurcation:

$$& \text{bifurcation} \colon \\ \alpha _ { 2 } \alpha _ { 1 } ^ { * } & \Delta _ { 1 } = \Delta _ { 1, 0 } + \lambda _ { 1 2 } \sqrt { \kappa _ { 2 } } \, \bar { a } _ { 2 } ^ { \text{in} } | \chi _ { 2 } | \left [ \cos ( \delta + \phi _ { 2 } ) + \sqrt { 3 } \sin ( \delta + \phi _ { 2 } ) \right ] \\ \frac { \lambda _ { 1 2 } } { 2 } \alpha _ { 1 } ^ { 2 } & \bar { a } _ { 1 } ^ { \text{in} } = \bar { a } _ { 1, 0 } ^ { \text{in} } + \lambda _ { 1 2 } \sqrt { \kappa _ { 2 } } \, \bar { a } _ { 2 } ^ { \text{in} } \, | \chi _ { 2 } | \left [ \sin ( \delta + \phi _ { 2 } ) \left ( \sqrt { 3 } / 2 \Lambda \right ) ^ { 1 / 2 } \right ] \\ ( B 1 ) &$$

where φ 2 is defined via χ 2 = | χ 2 | exp( iφ 2 ) as before. An important observation - borne out by the numerics - is that the critical detuning and driving strength values are affected by the relative drive phase δ in particular, indicating that tuning δ nontrivially modifies the bifurcation physics.

## Appendix C: Effective nonlinear cavity

The set of four coupled Heisenberg equations of motion in the displaced frame, Eqs. (15), form a four-dimensional system. It is clearest to solve this system by expressing it in the canonical quadrature basis; this transform is carried out using

Eq. (C1) is also used to transform the input noise operators ˆ ξ j to the canonical basis, with the replacements ( ˆ d , d j ˆ ) † j → ( ˆ ξ , ξ j ˆ ) † j , (ˆ x , p j ˆ ) j → (ˆ x in j , p ˆ in j ) . Following this, the full fourdimensional system in the canonical basis can be written in Fourier space as,

$$( i \omega 1 + \mathbf M _ { c } ) \begin{pmatrix} \hat { x } _ { 1 } [ \omega ] \\ \hat { p } _ { 1 } [ \omega ] \\ \hat { x } _ { 2 } [ \omega ] \\ \hat { p } _ { 2 } [ \omega ] \end{pmatrix} = \begin{pmatrix} \sqrt { \kappa _ { 1 } } \, \hat { x } _ { 1 } ^ { \text{in} } [ \omega ] \\ \sqrt { \kappa _ { 1 } } \, \hat { p } _ { 1 } ^ { \text{in} } [ \omega ] \\ \sqrt { \kappa _ { 2 } } \, \hat { x } _ { 2 } ^ { \text{in} } [ \omega ] \\ \sqrt { \kappa _ { 2 } } \, \hat { p } _ { 2 } ^ { \text{in} } [ \omega ] \end{pmatrix} + \sqrt { 2 } \begin{pmatrix} 0 \\ - \cos \left ( \frac { \phi _ { x } } { 2 } + \frac { 3 \pi } { 4 } \right ) \\ \sin \left ( \frac { \phi _ { x } } { 2 } + \frac { 3 \pi } { 4 } \right ) \end{pmatrix} \hat { B } [ \omega ]$$

Here ˆ [ B ω ] = ˜ Aσ ˆ [ z ω ] describes the qubit's frequencydependent driving force on the linear cavity via the linearised dispersive coupling, and the susceptibility matrix M c is:

The off-diagonal matrices - the cross-terms - describe the intercavity coupling; in particular, these appear in the form of an effective rotation defined by the matrix R ( θ ) :

$$\mathbf M _ { c } = \begin{bmatrix} \mathbf M _ { 1 } & - \lambda _ { 1 2 } \sqrt { \bar { n } _ { 1 } } \, \mathbf R ( - \mu _ { 1 2 } - \pi / 2 ) \\ \lambda _ { 1 2 } \sqrt { \bar { n } _ { 1 } } \, \mathbf R ( \mu _ { 1 2 } + \pi / 2 ) & \mathbf M _ { 2 } \\ \quad & ( C 3 ) \\ \text{where the nonlinear cavitv response matrix } \mathbf M. \text{ and linear} \end{bmatrix} \quad \mathbf R ( \theta$$

where the nonlinear cavity response matrix M 1 and linear cavity response matrix M 2 take the forms:

$$M _ { 1 } = \begin{bmatrix} g _ { \text{eff} } - \frac { \kappa _ { 1 } } { 2 } & - \widetilde { \Delta } _ { 1 } \\ \widetilde { \Delta } _ { 1 } & - g _ { \text{eff} } - \frac { \kappa _ { 1 } } { 2 } \end{bmatrix}, \, M _ { 2 } = \begin{bmatrix} - \frac { \kappa _ { 2 } } { 2 } & - \Delta _ { 2 } \\ \Delta _ { 2 } & - \frac { \kappa _ { 2 } } { 2 } \end{bmatrix} \quad \begin{array} { c } \text{The $r$} \\ \text{Eq. (2)} \\ \text{total} \\ \text{to exp} \\ \text{aratel} \vdots \end{array}$$

The rotation is characterised by the angle µ 12 defined in Eq. (22). It is convenient to use the block matrix structure to expand Eq. (C2) into equations for both cavity modes separately:

$$( i \omega 1 + \mathbf M _ { 1 } ) \begin{pmatrix} \hat { x } _ { 1 } [ \omega ] \\ \hat { p } _ { 1 } [ \omega ] \end{pmatrix} - \lambda _ { 1 2 } \sqrt { \bar { n } _ { 1 } } \, \mathbf R \, ( - \mu _ { 1 2 } - \pi / 2 ) \begin{pmatrix} \hat { x } _ { 2 } [ \omega ] \\ \hat { p } _ { 2 } [ \omega ] \end{pmatrix} = \sqrt { \kappa _ { 1 } } \begin{pmatrix} \hat { x } _ { 1 } ^ { \infty } [ \omega ] \\ \hat { p } _ { 1 } ^ { \infty } [ \omega ] \end{pmatrix}$$

$$( i \omega 1 + \mathbf M _ { 2 } ) \begin{pmatrix} \hat { x } _ { 2 } [ \omega ] \\ \hat { p } _ { 2 } [ \omega ] \end{pmatrix} + \lambda _ { 1 2 } \sqrt { \frac { \cdot } { n } } _ { 1 } \mathbf R \left ( \mu _ { 1 2 } + \pi / 2 \right ) \begin{pmatrix} \dot { x } _ { 1 } [ \omega ] \\ \hat { p } _ { 1 } [ \omega ] \end{pmatrix} = \sqrt { \kappa _ { 2 } } \begin{pmatrix} \dot { x } _ { 2 } ^ { \infty } [ \omega ] \\ \hat { p } _ { 2 } ^ { \infty } [ \omega ] \end{pmatrix}$$

Note that we ignore the coupling to the qubit for the moment for clarity. A crucial part of our analysis involves solving Eq. (C6b) for the linear cavity quadratures and substituting the results into Eq. (C6a), to obtain a set of equations for the cavity-1 quadratures only. We thus obtain Eq. (24), an effec- tive two-dimensional system with the linear cavity eliminated:

$$\begin{pmatrix} \hat { x } _ { 1 } [ \omega ] \\ \hat { p } _ { 1 } [ \omega ] \end{pmatrix} = \mathbf M _ { \text{eff} } ^ { - 1 } [ \omega ] \begin{pmatrix} \hat { x } _ { e } ^ { \text{in} } [ \omega ] \\ \hat { p } _ { e } ^ { \text{in} } [ \omega ] \end{pmatrix} \quad \quad ( C 7 )$$

The effective response matrix M eff [ ω ] can be written in terms of the original nonlinear cavity response matrix M 1 together

$$\stackrel { \bullet } { 3 } \right ] _ { \text{or} } \quad \mathbf R ( \theta ) = \begin{bmatrix} \cos \theta & \sin \theta \\ - \sin \theta & \cos \theta \end{bmatrix} ; \quad \mathbf R ^ { - 1 } ( \theta ) = \mathbf R ( - \theta ) \quad ( \mathbf C 5 )$$

the matrix T c ,

$$\begin{smallmatrix} \text{the} \\ \text{by} \\ \hat { g } \, \delta \\ \hat { \hat { p } } _ { 1 } \\ \hat { \hat { x } } _ { 2 } \\ \hat { \hat { p } } _ { 2 } \end{smallmatrix} = \text{T} _ { c } \begin{pmatrix} \hat { d } _ { 1 } \\ \hat { d } _ { 1 } ^ { \dagger } \\ \hat { d } _ { 2 } ^ { \dagger } \end{pmatrix} = \frac { 1 } { \sqrt { 2 } } \begin{bmatrix} 1 & 1 & 0 & 0 \\ - i & i & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & - i & i \end{bmatrix} \begin{pmatrix} \hat { d } _ { 1 } \\ \hat { d } _ { 1 } ^ { \dagger } \\ \hat { d } _ { 2 } ^ { \dagger } \end{bmatrix} \quad ( C 1 )$$

with the self-energy matrix introduced in Eq. (25),

$$\begin{array} { c c c c c } \quad \,.. &. \cdot \cdot \cdot \cdot \\ M _ { \text{eff} } [ \omega ] \equiv \mathbf M _ { 1 } [ \omega ] - i \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } \Sigma [ \omega ] \\ \Sigma [ \omega ] = \frac { i } { 2 } \begin{bmatrix} ( \chi _ { 2 } [ \omega ] + \chi _ { 2 } ^ { * } [ - \omega ] ) & i ( \chi _ { 2 } [ \omega ] - \chi _ { 2 } ^ { * } [ - \omega ] ) \\ - i ( \chi _ { 2 } [ \omega ] - \chi _ { 2 } ^ { * } [ - \omega ] ) & ( \chi _ { 2 } [ \omega ] + \chi _ { 2 } ^ { * } [ - \omega ] ) \end{bmatrix}$$

where the full linear cavity susceptibility takes the form

$$\chi _ { 2 } [ \omega ] = \left ( - i \omega + \kappa _ { 2 } / 2 \pm i \Delta _ { 2 } \right ) ^ { - 1 } \quad \text{(C9)}$$

and M 2 [ ω ] = iω 1 + M M 2 . eff [ ω ] is given explicitly by,

$$\underbrace { \mathbf M _ { \text{eff} } [ \omega ] = \begin{bmatrix} g _ { \text{eff} } - \frac { \kappa _ { \text{eff} } } { 2 } [ \omega ] & - \widetilde { \Delta } _ { \text{eff} } [ \omega ] \\ \widetilde { \Delta } _ { \text{eff} } [ \omega ] & - g _ { \text{eff} } - \frac { \kappa _ { \text{eff} } } { 2 } [ \omega ] \end{bmatrix} \quad ( C 1 0 ) \quad \text{Final} \quad \text{filed,} \quad \text{$\varepsilon$}$$

Finally, the drive vector on the nonlinear cavity is also modified, as describe in Eqs. (24); this can then be expressed as:

$$\binom { \hat { x } _ { e } ^ { \text{in} } [ \omega ] } { \hat { p } _ { e } ^ { \text{in} } [ \omega ] } = \sqrt { \kappa _ { 1 } } \binom { \hat { x } _ { 1 } ^ { \text{in} } [ \omega ] } { \hat { p } _ { 1 } ^ { \text{in} } [ \omega ] } - \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \left \{ \mathbf R \left ( \mu _ { 1 2 } + \pi / 2 \right ] \mathbf M _ { 2 } [ \omega ] \right \} ^ { - 1 } \sqrt { \kappa _ { 2 } } \binom { \hat { x } _ { 2 } ^ { \text{in} } [ \omega ] } { \hat { p } _ { 2 } ^ { \text{in} } [ \omega ] } \right )$$

We now reintroduce the coupling to the qubit in Eqs. (C2); the only change is to the modified drive on the effective system. It can be obtained from Eq. (C12) using the replacement:

where all quadratures are understood to be frequency dependent, the effective two-dimensional system of Eq. (24) can be written in this rotated basis as:

$$\sqrt { \kappa _ { 2 } } \begin{pmatrix} \hat { x } _ { 2 } ^ { \text{in} } \\ \hat { p } _ { 2 } ^ { \text{in} } \end{pmatrix} & \rightarrow \sqrt { \kappa _ { 2 } } \begin{pmatrix} \hat { x } _ { 2 } ^ { \text{in} } \\ \hat { p } _ { 2 } ^ { \text{in} } \end{pmatrix} + \sqrt { 2 } \begin{pmatrix} - \cos \left ( \frac { \phi _ { q } } { 2 } + \frac { 3 \pi } { 4 } \right ) \\ \sin \left ( \frac { \phi _ { q } } { 2 } + \frac { 3 \pi } { 4 } \right ) \end{pmatrix} \hat { B } [ \omega ] \\ \text{(where the frequency labels of the quadratures are $\sin$-\quad where} }$$

(where the frequency labels of the quadratures are suppressed).

The matrix M eff [ ω ] (cf. Eq. (C10)) bears similarity to the susceptibility matrix for a driven, damped nonlinear cavity with no coupling to a linear cavity. The full system is solved easily by first diagonalizing M eff [ ω ] ,

$$M _ { \text{eff} } [ \omega ] = - V _ { e } \begin{bmatrix} ( \chi _ { e - } [ \omega ] ) ^ { - 1 } & 0 \\ 0 & ( \chi _ { e + } [ \omega ] ) ^ { - 1 } \end{bmatrix} V _ { e } ^ { - 1 } \ \ ( C 1 4 ) \quad \text{freque} \quad \text{by ad}$$

where V e , the matrix of eigenvectors of M eff [ ω ] , is:

$$\mathbf V _ { e } = \begin{bmatrix} \cos ( \theta _ { e } / 2 ) & \sin ( \theta _ { e } / 2 ) \\ \sin ( \theta _ { e } / 2 ) & \cos ( \theta _ { e } / 2 ) \end{bmatrix} \quad \text{(C15)}$$

θ e is defined in Eq. (28), and the frequency-dependent effective nonlinear cavity susceptibilities are simply χ e ± [ ω ] -1 = ( -iω + χ -1 e ± ) , with χ e ± defined in Eq. (31a), (31b).

The analysis is most transparent in the rotated basis of amplified and squeezed quadratures, ˆ X ,P e ˆ e respectively, as discussed earlier. Using the transformation:

$$\begin{pmatrix} \hat { X } _ { e } \\ \hat { P } _ { e } \end{pmatrix} = \mathbf R ( \theta _ { e } / 2 ) \begin{pmatrix} \hat { x } _ { 1 } \\ \hat { p } _ { 1 } \end{pmatrix} \ ; \ \begin{pmatrix} \hat { X } _ { e } ^ { \text{in} } \\ \hat { P } _ { e } ^ { \text{in} } \end{pmatrix} = \mathbf R ( \theta _ { e } / 2 ) \begin{pmatrix} \hat { x } _ { e } ^ { \text{in} } \\ \hat { p } _ { e } ^ { \text{in} } \end{pmatrix} \quad \begin{pmatrix} \text{terms} \\ \hat { P } _ { L } ^ { \text{int} } [ \mu \\ ( C 1 6 ) \quad & \text{noise} \, \, \,$$

$$\begin{pmatrix} \hat { X } _ { e } [ \omega ] \\ \hat { P } _ { e } [ \omega ] \end{pmatrix} = \mathbf M _ { \text{eff} } ^ { - 1 } [ \omega ] \begin{pmatrix} \hat { X } _ { e } ^ { \text{in} } [ \omega ] \\ \hat { P } _ { e } ^ { \text{in} } [ \omega ] \end{pmatrix} \quad \text{(C17)}$$

where M -1 eff [ ω ] is given in terms of χ e ± [ ω ] by:

$$\text{the} _ { \text{iv} } \quad M _ { \text{eff} } ^ { - 1 } [ \omega ] = \begin{bmatrix} - \chi _ { e - } [ \omega ] & ( \chi _ { e - } [ \omega ] - \chi _ { e + } [ \omega ] ) \tan \theta _ { e } \\ 0 & - \chi _ { e + } [ \omega ] \end{bmatrix} \ ( C 1 8 )$$

The effective cavity output quadratures are then simply obtained by multiplying out Eqs. (C17), leading to the frequency-dependent versions of Eqs. (30a), (30b) (acquired by adding frequency labels to all quantities).

## Appendix D: Backaction force noise calculations

## a. Linearised backaction calculation

We begin by elaborating on the calculation of the linearised backaction, and start again by express the linearised backaction force operator ˆ F L (cf. Eq. (37)) in the rotated quadrature basis described in Appendix C. This is achieved by eliminating the linear cavity modes from ˆ F L , thereby expressing it in terms of the nonlinear cavity modes only. Upon transforming to the ˆ X , P e ˆ e basis, ˆ F L takes the form given in Eq. (38). The ˆ F int L [ ω ] operator is given in terms of the linear cavity input noise operators ˆ [ ξ 2 ω ] only , and is in particular independent of where the effective nonlinear cavity detuning, damping, and parametric interaction strength are:

$$\widetilde { \Delta } _ { \text{eff} } [ \omega ] & = \widetilde { \Delta } _ { 1 } + \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } \frac { i } { 2 } \left ( \chi _ { 2 } [ \omega ] - \chi _ { 2 } ^ { * } [ - \omega ] \right ) \\ \kappa _ { \text{eff} } [ \omega ] & = \kappa _ { 1 } + \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 2 } ( \chi _ { 2 } [ \omega ] + \chi _ { 2 } ^ { * } [ - \omega ] ) \\ g _ { \text{eff} } & = \left | g _ { 1 } - \widetilde { \lambda } _ { 1 2 } ^ { ( 2 ) } e ^ { i \mu _ { 1 2 } } \right | & ( C 1 1 )$$

λ 12 :

$$\hat { F } _ { L } ^ { \text{int} } [ \omega ] = \sqrt { \kappa _ { 2 } } \, \widetilde { A } & \left ( \chi _ { 2 } [ \omega ] e ^ { i ( \phi _ { g } / 2 - 3 \pi / 4 ) } \hat { \xi } _ { 2 } [ \omega ] & S _ { F _ { i } } \\ & + \chi _ { 2 } ^ { * } [ - \omega ] e ^ { - i ( \phi _ { g } / 2 - 3 \pi / 4 ) } \hat { \xi } _ { 2 } ^ { \dagger } [ \omega ] \right ) \quad ( D 1 )$$

The λ 12 dependent backaction terms involve f X [ ω ] , defined in Eq. (39), and f P [ ω ] , which is defined as:

$$f _ { P } [ \omega ] = i \frac { \widetilde { A } } { \sqrt { 2 } } \, \left [ \chi _ { 2 } ^ { * } [ - \omega ] e ^ { i ( \mu _ { 1 2 } - \nu [ \omega ] ) } - \chi _ { 2 } [ \omega ] e ^ { - i ( \mu _ { 1 2 } - \nu [ \omega ] ) } \right ] _ { \quad \text{Fuctu} \\ ( D 2 ) \quad \text{types} }.$$

We now move on to characterising the noise spectral density of the linearised backaction force. For calculations at nonzero frequencies, we compute the symmetrized noise spectral density, ¯ S FF [ ω ] , where:

$$\bar { S } _ { F F } [ \omega ] = \frac { 1 } { 2 } \left ( S _ { F F } [ \omega ] + S _ { F F } [ - \omega ] \right ) \quad \text{(D3)} \quad \text{if} \quad.$$

The unsymmetrized spectral density S FF [ ω ] , given by Eq. (36), can be computed using an application of the WienerKhinchin theorem, and the form of ˆ [ F L ω ] . It is simple to relate ˆ F L to the vacuum drives on both cavities, ˆ X ,P in j ˆ in j in the rotated basis using the solutions to Eqs. (C17). S FF [ ω ] , and hence ¯ S FF [ ω ] , can then be calculated using the usual correlation functions for these vacuum noise drives at zero temperature,

$$\langle \hat { X } _ { j } ^ { \text{in} } [ \omega ] \hat { X } _ { j } ^ { \text{in} } [ \omega ^ { \prime } ] \rangle & = \pi \delta ( \omega + \omega ^ { \prime } ) & S _ { F _ { c } } ^ { ( 1 } \\ \langle \hat { P } _ { j } ^ { \text{in} } [ \omega ] \hat { P } _ { j } ^ { \text{in} } [ \omega ^ { \prime } ] \rangle & = \pi \delta ( \omega + \omega ^ { \prime } ) & \\ \langle \hat { X } _ { j } ^ { \text{in} } [ \omega ] \hat { P } _ { j } ^ { \text{in} } [ \omega ^ { \prime } ] \rangle & = i \pi \delta ( \omega + \omega ^ { \prime } ) & & \text{(D4)} & & \text{On}$$

## b. Quadratic backaction calculation

The second backaction source is the quadratic backaction operator ˆ F Q (cf. Eq. (37)). The effect of this quadratic backaction can be computed from the relevant autocorrelation function G F Q F Q ( ) t ,

$$& G _ { F _ { Q } F _ { Q } } ( t ) = \langle \langle \hat { F } _ { Q } ( t ) \hat { F } _ { Q } ( 0 ) \rangle \rangle = A ^ { 2 } \langle \langle \hat { d } _ { 2 } ^ { \dagger } ( t ) \hat { d } _ { 2 } ( t ) \hat { d } _ { 2 } ^ { \dagger } ( 0 ) \hat { d } _ { 2 } ( 0 ) \rangle \rangle \quad \text{The                                                                                                                                                                                                        
                                                                                                                                                                                                        }                                                                                                                                                                                                        \end{                                                                                                                                                                                                        0                                                                                                                                                                                                        1                                                                                                                                                                                                        2                                                                                                                                                                                                        3                                                                                                                                                                                                        4                                                                                                                                                                                                        6                                                                                                                                                                                                        7                                                                                                                                                                                                        8                                                                                                                                                                                                        9                                                                                                                                                                                                        10                                                                                                                                                                                                        5                                                                                                                                                                                                        20                                                                                                                                                                                                        19                                                                                                                                                                                                        -                                                                                                                                                                                                        ($$

We retain only connected terms, as indicated by the double angled brackets, and calculate the correlation function in the weak dispersive coupling regime so that the qubit-cavity coupling can be ignored; the remaining system Hamiltonian in terms of the ˆ d j operators - given by the sum ˆ H d 1 + ˆ H d 2 + ˆ H d 12 -is quadratic [cf. Eqs. (16), (18), (20)]. We therefore use Wick's theorem in going from the first line to the second to simplify the average. Fourier transforming the autocorrelation function G F Q F Q ( ) t allows computing the noise spectral density of the quadratic backaction, S F Q F Q [ ω ] , following ap- plication of the Wiener-Khinchin theorem [4],

$$\text{pawait on any number} \text{mm} \text{mm} \text{unused} \{ \tau _ { 1 }, \\ S _ { F _ { \text{Q} } F _ { \text{Q} } } [ \omega ] = A ^ { 2 } \int _ { - \infty } ^ { \infty } \frac { d \omega ^ { \prime } } { 2 \pi } & \left \{ S _ { d _ { 2 } ^ { \dagger } d _ { 2 } } [ \omega ^ { \prime } ] S _ { d _ { 2 } d _ { 2 } ^ { \dagger } } [ - \omega ^ { \prime } - \omega ] \\ + & S _ { d _ { 2 } ^ { \dagger } d _ { 2 } ^ { \dagger } } [ \omega ^ { \prime } ] S _ { d _ { 2 } d _ { 2 } } [ - \omega ^ { \prime } - \omega ] \right \} \\ \text{need} & \quad \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat{ \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat { \ } \hat$$

where the spectral density of two operators ˆ X, Y ˆ is simply:

$$S _ { X Y } [ \omega ] = \langle \hat { X } [ \omega ] \hat { Y } [ - \omega ] \rangle \text{ \quad \ \ } \text{(D7)}$$

Fluctuations in the linear cavity mode ˆ d 2 arise from three types of noise terms; in addition to noise incident directly on the linear cavity (a contribution independent of λ 12 ), ˆ d 2 is driven by noise incident on the nonlinear cavity that enters the linear cavity by virtue of the cavity-cavity coupling, and is hence ∝ λ 12 . Lastly, terms ∝ λ 2 12 involve noise incident on the linear cavity passing to the nonlinear cavity, being amplified, and then returning to drive ˆ d 2 .

For some analytic insight into the quadratic backaction, we focus on the zero frequency noise S F Q F Q [0] . The integrals over frequency space mean that S F Q F Q [0] sees noise at all frequencies (cf. Eq. (D6)), and as a result we expect large contributions near two frequencies, as indicated in Eq. (46). Near ω = 0 , the nonlinear cavity response G [0] is the dominant contribution, which also sets the bandwidth Ω[ ω = 0] = κ / 1 √ G . The linear cavity response here is set by the detuning ∆ 2 since the linear cavity is strongly detuned. The contribution to the quadratic backaction near ω = 0 is then:

$$S _ { F _ { \mathbb { Q } } F _ { \mathbb { Q } } } ^ { ( 1 ) } \simeq A ^ { 2 } \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 4 } \cdot \frac { 1 } { \Delta _ { 2 } ^ { 4 } } \cdot \frac { \mathcal { G } ^ { 2 } } { \kappa _ { 1 } ^ { 2 } } \cdot \frac { \kappa _ { 1 } } { \sqrt { \mathcal { G } } } = A ^ { 2 } \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 4 } \cdot \frac { \mathcal { G } ^ { 3 / 2 } } { \kappa _ { 1 } \Delta _ { 2 } ^ { 4 } } \\ \intertext { a n d }$$

On the other hand, near ω = ∆ 2 , the resonant linear cavity susceptibility dominates the backaction as | χ 2 | glyph[similarequal] 1 /κ 2 , and the bandwidth is therefore Ω[ ω = ∆ ] = 2 κ 2 . The nonlinear cavity response here is suppressed by a factor of the large linear cavity detuning, so that G [ ω = ∆ ] 2 glyph[similarequal] ( κ / 1 ∆ ) 2 2 . The quadratic backaction near ω = ∆ 2 is then

$$\text{lract} \quad S _ { F _ { Q } F _ { Q } } ^ { ( 2 ) } \simeq A ^ { 2 } \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 4 } \cdot \frac { 1 } { \kappa _ { 2 } ^ { 4 } } \cdot \frac { \kappa _ { 1 } ^ { 2 } } { \Delta _ { 2 } ^ { 4 } } \cdot \kappa _ { 2 } = A ^ { 2 } \left ( \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \right ) ^ { 4 } \cdot \frac { \kappa _ { 1 } ^ { 2 } } { \kappa _ { 2 } ^ { 3 } \Delta _ { 2 } ^ { 4 } }$$

The above estimates for the main contributions to the quadratic backaction integrals (cf. Eq. (D6)) are found to approximate the full integrals well. In the optimization procedure detailed in Section (IV A), these terms are required to be small relative to the intrinsic nonlinear cavity backaction, leading to the constraints in Eq. (55).

## Appendix E: Output homodyne current calculations

The full output homodyne current from the nonlinear cavity is the result of interference of the classical homodyne reference beam and the cavity output field; it can be written as [42]:

$$\cdot \quad \hat { I } _ { 1 } [ \omega ] = \sqrt { \kappa _ { 1 } } ( \cos \phi _ { h } \, \hat { x } _ { 1 } ^ { \text{out} } [ \omega ] + \sin \phi _ { h } \, \hat { p } _ { 1 } ^ { \text{out} } [ \omega ] ) \quad ( \text{E1} )$$

FIG. 6. Plot of the amplifier forward gain scaled by the dispersive qubit-cavity coupling, | χ IF [ ω /A ] | (dashed blue, left-hand axis) and the inverse efficiency ratio 1 /η meas (solid red, right-hand axis), both as a function of frequency ω , for the same parameter set as Fig. 3. As frequency increases, 1 /η meas first increases slowly as the reduced backaction condition is no longer satisfied. However, the concurrent reduction in gain G [ ω ] leads eventually to a reduced value of 1 /η meas , as expected. The vertical dashed line indicates the measurement rate ∼ 0 1 . κ 1 for our parameter choices. We see that over this frequency range, | χ IF [ ω ] | decreases somewhat; this is due to Γ meas becoming comparable to the amplification bandwidth, a feature that is controlled by the experimental value of κ 1 . Importantly, 1 /η meas increases relative to the zero frequency value by only a small amount, while the forward gain decreases by less than a factor of 2 ; our interest in the zero frequency performance of the PARI is therefore justified.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-b7ff420e-1396-4407-b51b-6a46e37e8879.jpg)

In standard input-output theory, these cavity output quadratures can be written in terms of the reflected input drives, and the intracavity field leaking out of the cavity [4]. Using the solutions to Eqs. (C17) in the rotated basis, the output homodyne current can be written in terms of the amplified and squeezed quadratures. In the particular case of the average homodyne current, the vacuum noise drive vanishes, leaving:

$$\langle \hat { I } _ { 1 } [ \omega ] \rangle = \kappa _ { 1 } \left [ \cos \phi _ { h } \langle \hat { X } _ { e } [ \omega ] \rangle + \sin \phi _ { h } \langle \hat { P } _ { e } [ \omega ] \rangle \right ] \quad ( \text{E2} ) \quad \underset { \dots \dots } { \overset { \cdots \dots } { \text{ invers} } }$$

where we take φ h → φ h + θ / e 2 to absorb angles associated with the transformation. The cavity drives can be related to the qubit states using the frequency-dependent version of Eq. (47):

$$\hat { \zeta } _ { e } ^ { \text{in} } [ \omega ] & \to \hat { \zeta } _ { e } ^ { \text{in} } [ \omega ] - \\ & \left ( \frac { \widetilde { A } \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } } { \sqrt { 2 } } \left [ i \chi _ { 2 } [ \omega ] e ^ { i ( \mu _ { 1 2 } - \nu ) } - i \chi _ { 2 } ^ { * } [ - \omega ] e ^ { i ( \mu _ { 1 2 } - \nu ) } \right ] \right ) \hat { \sigma } _ { z } } { \quad \text{indica} } \\ \hat { P } _ { e } ^ { \text{in} } [ \omega ] & \to \hat { P } _ { e } ^ { \text{in} } [ \omega ] + \\ & \left ( \frac { \widetilde { A } \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } } { \sqrt { 2 } } \left [ \chi _ { 2 } [ \omega ] e ^ { i ( \mu _ { 1 2 } - \nu ) } + \chi _ { 2 } ^ { * } [ - \omega ] e ^ { i ( \mu _ { 1 2 } - \nu ) } \right ] \right ) \hat { \sigma } _ { z } \quad ( \text{E3} ) \quad \text{drive} \\ & \quad 1 / \eta _ { \text{me} } \\ \text{where} \, we cunnree the frequency label on \, \hat { \zeta } \quad \text{Following this} \quad \text{close}$$

where we suppress the frequency label on ˆ σ z . Following this substitution into Eq. (C17) and taking the ensemble average, the averaged quadratures 〈 ˆ X ω , P e [ ] 〉 〈 ˆ [ e ω ] 〉 can be substituted in Eq. (E2). The amplifier forward gain for arbitrary frequency is given simply by the frequency-dependence version of the linear response relation in Eq. (48),

$$\overline { \frac { 3 } { 2 } } \quad \chi _ { I F } [ \omega ] = \frac { \widetilde { A } \, \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } } { \sqrt { 2 } } \kappa _ { 1 } \sec \theta _ { e } \left ( S _ { 1 } \chi _ { e - } [ \omega ] + S _ { 2 } \chi _ { e + } [ \omega ] \right ) \text{ (E4)}$$

where the coefficients S 1 [ ω , S ] 2 [ ω ] are given by

$$S _ { 1 } & = \cos \phi _ { h } \times \\ & \quad - \left ( i \chi _ { 2 } [ \omega ] e ^ { i ( \mu _ { 1 2 } - \nu + \theta _ { e } ) } - i \chi _ { 2 } ^ { * } [ - \omega ] e ^ { - i ( \mu _ { 1 2 } - \nu + \theta _ { e } ) } \right ) \\ S _ { 2 } & = - \sin \phi _ { h } \left ( \chi _ { 2 } [ \omega ] e ^ { i ( \mu _ { 1 2 } - \nu ) } + \chi _ { 2 } ^ { * } [ - \omega ] e ^ { - i ( \mu _ { 1 2 } - \nu ) } \right ) \text{ (E5)}$$

Near bifurcation, where the S 1 term dominates, the magnitude of the gain reduces to the result in Eq. (48). Furthermore, in the strongly detuned regime where φ 2 → π/ 2 , and once the reduced backaction condition (cf. Eq. (42)) is satisfied, the angular term ρ e in Eq. (48) takes the form:

$$\rho _ { e } & = \sec \theta _ { e } \, \sin \left [ ( 2 M + 1 ) \frac { \pi } { 2 } + \pi + \theta _ { e } \right ] \cos \phi _ { h } \\ & \Longrightarrow \, \sec \theta _ { e } \, \cos \theta _ { e } \, \cos \phi _ { h } = \cos \phi _ { h }$$

Clearly, choosing φ h = 0 sets ρ e → 1 , allowing us to drop this angular factor.

The imprecision noise spectral density of the output signal is defined in terms of the full output homodyne current for zero qubit coupling, analogously to S FF [ ω ] (cf. Eq. (36)),

$$S _ { I I } [ \omega ] \equiv \int _ { - \infty } ^ { \infty } d t \, e ^ { i \omega t } \langle \hat { I } ( t ) \hat { I } ( 0 ) \rangle _ { 0 }. \quad \text{(E7)}$$

S II [ ω ] is straightforward to compute by expressing the full output homodyne current in terms of ˆ X ,P e ˆ e , and making use of the usual correlation functions for the input drives [cf. Eqs (D4)].

While the complete ω -dependent expressions are unwieldy, we can observe the frequency dependence by computing the inverse efficiency ratio 1 /η meas as a function of frequency for parameter choices indicated in Table II; this is plotted in solid red in Fig. 6. We also include the forward gain | χ IF | (dashed blue).

## Appendix F: Freedom of parameter choice

In this section we include supplementary graphs that indicate the freedom of parameter choices afforded by the PARI scheme.

Figs. 7 (a), (b) are 3-D plots, against ∆ 1 and the cavity drive phase difference δ , of the inverse efficiency ratio 1 /η meas and the parametric gain G respectively, for fixed nonlinear cavity drive ¯ a in 1 ; all other parameters are fixed at the values indicated in Table II. The 'valley' (colored magenta) in Fig. 7 (a) indicates a region where the deviation from the quantum limit is at most a factor of 2 ( η meas ≥ 0 5 . ); clearly, a

FIG. 7. 3D plot against nonlinear cavity detuning ∆ 1 and relative phase difference δ between cavity drives, of (a) the inverse efficiency ratio 1 /η meas . The valley (magenta) indicates a region where η meas ≥ 0 5 . ; (b) the parametric gain G , for the same region as in (a). Clearly, in regions where the PARI performance is nearly quantumlimited [the valley in Fig. 7 (a)], G can be large, of O (20 dB) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-95d73da2-65fb-4966-abf6-a8b34231d91f.jpg)

nontrivial set of parameters exists where this is true. Fig. 7 (b) indicates large gain over parameters corresponding to this valley; therefore, there are multiple parameter sets that are good choices for operation of the PARI detector.

Fig. 8 is a plot of the gain G and inverse efficiency ratio 1 /η meas as a function of nonlinear cavity detuning ∆ 1 with a Kerr constant Λ = 10 -4 κ 1 . This is twice as strong as the nonlinearity used for Figs. 3, 4. We find at the optimal detuning an inverse efficiency of glyph[similarequal] 2 01 . , with G still being large, G ∼ 80 . Strengthening Λ reduces the allowed values

- [1] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. A 69 , 062320 (2004).
- [2] D. I. Schuster, A. Wallraff, A. Blais, L. Frunzio, R.-S. Huang, J. Majer, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. Lett. 94 , 123602 (2005).
- [3] J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, A. Wallraff, A. Blais, M. Devoret, S. Girvin, and R. J. Schoelkopf,

of ¯ n j , and hence we pay a cost in terms of the measurement rate; Γ meas now is ∼ 1 MHz , a factor of 10 slower than that in Fig. 4. Here, the drive parameters are given by: (¯ a in 1 ) 2 /κ 1 = (43 8) . 2 , (¯ a in 2 ) 2 /κ 1 = (57 0) . 2 , ¯ n 1 ∼ 3000 , ¯ n 2 ∼ 700 , and δ glyph[similarequal] 0 06 . . Again, an optimization procedure in ¯ a in 1 -∆ 1 -δ space has been performed.

FIG. 8. Plot of the parametric gain G (dashed blue, left-hand axis) and the inverse efficiency ratio 1 /η meas (solid red, right-hand log axis) as a function of nonlinear cavity detuning ∆ 1 , for a stronger Λ than in Figs. 3, 4. For optimal detuning, we find 1 /η meas glyph[similarequal] 2 01 . , with G glyph[similarequal] 80 , and Γ meas glyph[similarequal] 1 MHz .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0226v2-figure-b8063dd8-b7af-46bb-8ad7-c55996f25dc4.jpg)

## Appendix G: Minimum useful value of λ 12

From a description of the output homodyne currents, it is possible to derive the minimum useful value of the intercavity coupling in the PARI scheme. We consider two types of qubit-cavity detector systems: the PARI detector, with output homodyne current ˆ I 1 , and a linear readout cavity capable of quantum limited qubit state measurement with unit gain, having corresponding output current ˆ I 2 . The forms of both currents are as follows:

$$\hat { I } _ { 1 } & \simeq \widetilde { A } \, \hat { \sigma } _ { z } \cdot \sqrt { \kappa _ { 1 } } \cdot | \chi _ { 2 } | \cdot \widetilde { \lambda } _ { 1 2 } ^ { ( 1 ) } \cdot \chi _ { e - } \\ \hat { I } _ { 2 } & \simeq \widetilde { A } \, \hat { \sigma } _ { z } \cdot \sqrt { \kappa _ { 2 } } \cdot | \chi _ { 2 } |$$

If ˆ I 1 and ˆ I 2 are equal, then using the nonlinear cavity serves no advantage over the single cavity system; this condition therefore leads to the minimum useful intercavity coupling strength, given in Eq. (34).

Nature 449 , 443 (2007).

- [4] A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and R. J. Schoelkopf, Rev. Mod. Phys. 82 , 1155 (2010).
- [5] B. Yurke, J. Opt. Soc. Am. B 4 , 1551 (1987).
- [6] B. Yurke, L. R. Corruccini, P. G. Kaminsky, L. W. Rupp, A. D. Smith, A. H. Silver, R. W. Simon, and E. A. Whittaker, Phys. Rev. A 39 , 2519 (1989).

- [7] M. A. Castellanos-Beltran and K. W. Lehnert, Appl. Phys. Lett. 91 , 083509 (2007).
- [8] M. A. Castellanos-Beltran, K. D. Irwin, G. C. Hilton, L. R. Vale, and K. W. Lehnert, Nature Phys. 4 , 929 (2008).
- [9] T. Yamamoto, K. Inomata, M. Watanabe, K. Matsuba, T. Miyazaki, W. D. Oliver, Y. Nakamura, and J. S. Tsai, Appl. Phys. Lett. 93 , 042510 (2008).
- [10] N. Bergeal, F. Schackert, M. Metcalfe, R. Vijay, V. E. Manucharyan, L. Frunzio, D. E. Prober, R. J. Schoelkopf, S. M. Girvin, and M. H. Devoret, Nature 465 , 64 (2010).
- [11] B. Abdo, F. Schackert, M. Hatridge, C. Rigetti, and M. H. Devoret, Appl. Phys. Lett. 99 , 162506 (2011).
- [12] N. Roch, E. Flurin, F. Nguyen, P. Morfin, P. Campagne-Ibarcq, M. H. Devoret, and B. Huard, Phys. Rev. Lett. 108 , 147701 (2012).
- [13] J. Y. Mutus, T. C. White, E. Jeffrey, D. Sank, R. Barends, J. Bochmann, Y. Chen, Z. Chen, B. Chiaro, A. Dunsworth, J. Kelly, A. Megrant, C. Neill, P. J. J. O'Malley, P. Roushan, A. Vainsencher, J. Wenner, I. Siddiqi, R. Vijay, A. N. Cleland, and J. M. Martinis, Appl. Phys. Lett. 103 , 122602 (2013).
- [14] I. Siddiqi, R. Vijay, F. Pierre, C. M. Wilson, M. Metcalfe, C. Rigetti, L. Frunzio, and M. H. Devoret, Phys. Rev. Lett. 93 , 207002 (2004).
- [15] I. Siddiqi, R. Vijay, M. Metcalfe, E. Boaknin, L. Frunzio, R. J. Schoelkopf, and M. H. Devoret, Phys. Rev. B 73 , 054510 (2006).
- [16] A. Lupas ¸cu, E. F. C. Driessen, L. Roschier, C. J. P. M. Harmans, and J. E. Mooij, Phys. Rev. Lett. 96 , 127003 (2006).
- [17] R. Vijay, M. H. Devoret, and I. Siddiqi, Review of Scientific Instruments 80 , 111101 (2009).
- [18] F. R. Ong, M. Boissonneault, F. Mallet, A. Palacios-Laloy, A. Dewes, A. C. Doherty, A. Blais, P. Bertet, D. Vion, and D. Esteve, Phys. Rev. Lett. 106 , 167002 (2011).
- [19] C. Laflamme and A. A. Clerk, Phys. Rev. A 83 , 1155 (2011).
- [20] M. Boissonneault, A. Doherty, F. Ong, P. Bertet, D. Vion, D. Esteve, and A. Blais, Phys. Rev. A 85 , 022305 (2012).
- [21] C. Laflamme and A. A. Clerk, Phys. Rev. Lett. 109 , 123602 (2012).
- [22] M. Boissonneault, A. C. Doherty, F. R. Ong, P. Bertet, D. Vion, D. Esteve, and A. Blais, Phys. Rev. A 89 , 022324 (2014).
- [23] R. Vijay, D. H. Slichter, and I. Siddiqi, Phys. Rev. Lett. 106 , 110502 (2011).
- [24] J. E. Johnson, E. M. Hoskinson, C. Macklin, D. H. Slichter, I. Siddiqi, and J. Clarke, Phys. Rev. B 84 , 220503 (2011).
- [25] C. Eichler, Y. Salathe, J. Mlynek, S. Schmidt, and A. Wallraff, Phys. Rev. Lett. 113 , 110502 (2014).
- [26] B. Abdo, K. Sliwa, S. Shankar, M. Hatridge, L. Frunzio, R. Schoelkopf, and M. Devoret, Phys. Rev. Lett. 112 , 167701 (2014).
- [27] R. Vijay, C. Macklin, D. H. Slichter, S. J. Weber, K. W. Murch, R. Naik, A. N. Korotkov, and I. Siddiqi, Nature 490 , 77 (2012).
- [28] K. W. Murch, S. J. Weber, K. M. Beck, E. Ginossar, and I. Siddiqi, Nature 499 , 62 (2013).
- [29] K. W. Murch, S. J. Weber, C. Macklin, and I. Siddiqi, Nature 502 , 211 (2013).
- [30] G. de Lange, D. Rist` e, M. J. Tiggelman, C. Eichler, L. Tornberg, G. Johansson, A. Wallraff, R. N. Schouten, and L. DiCarlo, Phys. Rev. Lett. 112 , 080501 (2014).
- [31] M. Hatridge, R. Vijay, D. H. Slichter, J. Clarke, and I. Siddiqi, Phys. Rev. B 83 , 134501 (2011).
- [32] C. W. Gardiner and P. Zoller, Quantum Noise (Springer, Berlin, 2000).
- [33] W. Wustmann and V. Shumeiko, Phys. Rev. B 87 , 184501 (2013).
- [34] The input noise operators defined in Eq. (7) also depend on β , but the noise correlators are independent of this overall phase and hence it is not necessary to keep track of this.
- [35] H. Wiseman and G. Milburn, Quantum Measurement and Control (Cambridge University Press, 2010).
- [36] C. Eichler and A. Wallraff, EPJ Quantum Technology 1 , 2 (2014), 10.1140/epjqt2.
- [37] V. E. Manucharyan, E. Boaknin, M. Metcalfe, R. Vijay, I. Siddiqi, and M. Devoret, Phys. Rev. B 76 , 014524 (2007).
- [38] F. D. O. Schackert, A Practical Quantum-Limited Parametric Amplifier Based on the Josephson Ring Modulator , Ph.D. thesis, Yale University (2013).
- [39] B. Abdo, A. Kamal, and M. Devoret, Phys. Rev. B 87 , 014508 (2013).
- [40] A. Narla, K. M. Sliwa, M. Hatridge, S. Shankar, L. Frunzio, R. J. Schoelkopf, and M. H. Devoret, Appl. Phys. Lett. 104 , 232605 (2014).
- [41] X. Zhou, V. Schmitt, P. Bertet, D. Vion, W. Wustmann, V. Shumeiko, and D. Esteve, Phys. Rev. B 89 , 214517 (2014).
- [42] C. W. Gardiner and M. J. Collett, Phys. Rev. A 31 , 3761 (1985).
Saeed Khan R. Vijay I. Siddiqi Aashish A. Clerk
Large gain quantum-limited qubit measurement using a two-mode nonlinear cavity
We provide a thorough theoretical analysis of qubit state measurement in a
setup where a driven, parametrically-coupled cavity system is directly coupled
to the qubit, with one of the cavities having a weak Kerr nonlinearity. Such a
system could be readily realized using circuit QED architectures. We
demonstrate that this setup is capable in the standard linear-response regime
of both producing a highly amplified output signal while at the same time
achieving near quantum-limited performance: the measurement backaction on the
qubit is near the minimal amount required by the uncertainty principle. This
setup thus represents a promising route for performing efficient large-gain
qubit measurement that is completely on-chip, and that does not rely on the use
of circulators or complex non-reciprocal amplifiers.