## Abstract

We investigate the dynamics of a microwave-driven Josephson junction capacitively coupled to a lumped-element LC oscillator. In the regime of driving where the Josephson junction can be approximated as a Kerr oscillator, this minimal nonlinear system has been previously shown to exhibit a bistability in phase and amplitude. In the present study, we characterize the full phase diagram and show that besides a parameter regime exhibiting bistability, there is also a regime of self-oscillations characterized by a frequency comb in its spectrum. We discuss the mechanism of comb generation which appears to be different from those studied in microcavity frequency combs and mode-locked lasers. We then address the fate of the comblike spectrum in the regime of strong quantum fluctuations, reached when nonlinearity becomes the dominant scale with respect to dissipation. We find that the nonlinearity responsible for the emergence of the frequency combs also leads to its dephasing, leading to broadening and ultimate disappearance of sharp spectral peaks. Our study explores the fundamental question of the impact of quantum fluctuations for quantum systems which do not possess a stable fixed point in the classical limit.

Original language | English (US) |
---|---|

Article number | 153601 |

Journal | Physical review letters |

Volume | 120 |

Issue number | 15 |

DOIs | |

State | Published - Apr 9 2018 |

## All Science Journal Classification (ASJC) codes

- Physics and Astronomy(all)

## Access to Document

## Other files and links

## Fingerprint Dive into the research topics of 'Frequency Combs in a Lumped-Element Josephson-Junction Circuit'. Together they form a unique fingerprint.

## Cite this

- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver

*Physical review letters*,

*120*(15), [153601]. https://doi.org/10.1103/PhysRevLett.120.153601

**Frequency Combs in a Lumped-Element Josephson-Junction Circuit**. In: Physical review letters. 2018 ; Vol. 120, No. 15.

}

*Physical review letters*, vol. 120, no. 15, 153601. https://doi.org/10.1103/PhysRevLett.120.153601

**Frequency Combs in a Lumped-Element Josephson-Junction Circuit.** / Khan, Saeed; Türeci, Hakan E.

Research output: Contribution to journal › Article › peer-review

TY - JOUR

T1 - Frequency Combs in a Lumped-Element Josephson-Junction Circuit

AU - Khan, Saeed

AU - Türeci, Hakan E.

N1 - Funding Information: Khan Saeed Türeci Hakan E. Department of Electrical Engineering, Princeton University , Princeton, New Jersey 08544, USA 9 April 2018 13 April 2018 120 15 153601 3 November 2017 © 2018 American Physical Society 2018 American Physical Society We investigate the dynamics of a microwave-driven Josephson junction capacitively coupled to a lumped-element L C oscillator. In the regime of driving where the Josephson junction can be approximated as a Kerr oscillator, this minimal nonlinear system has been previously shown to exhibit a bistability in phase and amplitude. In the present study, we characterize the full phase diagram and show that besides a parameter regime exhibiting bistability, there is also a regime of self-oscillations characterized by a frequency comb in its spectrum. We discuss the mechanism of comb generation which appears to be different from those studied in microcavity frequency combs and mode-locked lasers. We then address the fate of the comblike spectrum in the regime of strong quantum fluctuations, reached when nonlinearity becomes the dominant scale with respect to dissipation. We find that the nonlinearity responsible for the emergence of the frequency combs also leads to its dephasing, leading to broadening and ultimate disappearance of sharp spectral peaks. Our study explores the fundamental question of the impact of quantum fluctuations for quantum systems which do not possess a stable fixed point in the classical limit. U.S. Department of Energy 10.13039/100000015 DE-SC0016011 Basic Energy Sciences 10.13039/100006151 In superconducting quantum circuits, the Josephson junction (JJ) is a lossless nonlinear element that provides critical functionality for various quantum information processing tasks [1–3] , from gate operations to readout and amplification, made possible by controlling JJ dynamics via its embedding circuit and effective drive. For instance, under strong coupling and weak excitation (relative to the intrinsic nonlinearity), JJ-based artificial atoms have enabled cavity QED implementations [4–9] that have been extensively discussed using open Jaynes-Cummings or Rabi models in single [10] and multimode regimes [11] . However, applications employing JJs under strong excitation conditions, for readout [12–14] and quantum-limited amplification [15–19] , require an understanding of dynamical instabilities that sensitively depend on the model of nonlinearity employed [20,21] . Here we investigate the dynamics of a shunted JJ when capacitively coupled to a microwave-driven linear resonator. The dynamics of such a system under strong driving have been theoretically [10,22] and experimentally [12] studied in the context of a high-power read-out scheme, and are found to exhibit a bistability between two states with distinct phase and amplitude. In the adiabatic regime where mode coupling is weaker than losses, the coupled system maps to a single coherently driven Kerr oscillator with renormalized parameters, exhibiting precisely this bistability [13,23–25] . However, we find that in the strong-coupling regime the nonlinear mode acquires a retarded self-interaction mediated by the linear mode, which changes the classical phase diagram dramatically: for certain drive and detuning ranges the system may have no stable fixed points, a phase not exhibited by the single coherently driven Kerr oscillator [26] . In this dynamical regime, nonzero frequency instabilities emerge as limit cycles, yielding discrete, equally spaced comblike spectra in the frequency domain. Such comb formation in coherently driven microresonators [27–31] and incoherently pumped mode-locked lasers [32,33] is often understood as an instability towards symmetric sideband growth via four-wave mixing, in an underlying resonator geometry supporting multiple spatial modes [34,35] and a distributed nonlinearity (while exceptions have been discussed as well [36–38] ). Our results indicate that the minimal manifestation of Kerr-mediated comb formation is embodied in a Kerr oscillator coupled to a linear oscillator. While limit cycles [39,40] and their modification under classical noise [38,41] have been extensively studied in classical systems, they are far less explored in the deep quantum regime [42,43] accessible to the lumped-element JJ circuit discussed here. Using master equation and phase-space simulations, we investigate the fate of comblike spectra as the nonlinearity is tuned from weak to strong (equivalently, high to low mode occupation at the instability threshold), so that the system moves from an expected semiclassical regime towards a well-defined quantum regime where a single photon can in principle trigger the comb instability. We find that while the nonlinearity, together with strong enough coupling to the linear mode, is necessary for limit cycles to emerge, this very nonlinearity introduces quantum noise that dephases the limit cycle; for weak noise, the dephasing time typically scales inversely with the strength of the nonlinearity. Model .—The model we study is realizable in lumped-element setups [Fig. 1(c) ], as well as JJ-embedded transmission-line resonators [44] . We assume that the nonlinearity of the junction can be approximated by its lowest order Kerr nonlinearity. The resulting model described by the Hamiltonian H ^ = H ^ a + H ^ b + H ^ g + H ^ d (See Supplemental Material [45] ) then is generic, consisting of a driven linear oscillator a (frequency ω a ) coupled to a nonlinear oscillator b (frequency ω b ) [Fig. 1(a) ]. H ^ d = η ( a ^ + a ^ † ) describes the drive (strength η ) in the frame rotating at the drive frequency ω d . The corresponding drive frame Hamiltonians of the modes are H ^ a = - Δ d a a ^ † a ^ and H ^ b = - Δ d b b ^ † b ^ - ( Λ / 2 ) b ^ † b ^ † b ^ b ^ , respectively, with Λ > 0 being the strength of the nonlinearity, and the frequency detunings defined as Δ d a / d b = ω d - ω a / b . The two oscillators are coupled linearly, with Hamiltonian H ^ g = g ( a ^ † b ^ + b ^ † a ^ ) . The system dynamics including damping for both the linear (rate κ ) and nonlinear mode (rate γ ) are then governed by the master equation ρ ^ ˙ = - i [ H , ρ ^ ] + κ D [ a ^ ] ρ ^ + γ D [ b ^ ] ρ ^ , where D [ o ^ ] is the standard dissipative superoperator D [ o ^ ] ρ = o ^ ρ o ^ † - 1 2 { o ^ † o ^ , ρ } . 1 10.1103/PhysRevLett.120.153601.f1 FIG. 1. (a) Schematic representation of the two-mode system. (b) Mode and drive frequencies. (c) Lumped-element circuit QED implementation of (a). We begin with the classical dynamics of the two-mode system, obtaining operator equations of motion and making the replacement ( ( ⟨ a ^ ⟩ , ⟨ b ^ ⟩ ) → ( α , β ) . This yields α ˙ = i Δ d a α - κ 2 α - i g β - i η β ˙ = i Δ d b β - γ 2 β + i Λ | β | 2 β - i g α . (1) Note that Eqs. (1) are invariant if Λ → Λ / c and ( α , β , η ) → c ( α , β , η ) , for c ∈ R + [45] . Physically, scaling Λ → Λ / c and η → c η yields the same classical dynamics, except with mode amplitudes scaled by c . This simple Λ -dependence is not true of the quantum dynamics, as we shall see later. Next, the linearity of mode a ^ and the coupling allows it to be integrated out exactly, leading to an effective dynamical equation for the nonlinear mode: β ˙ = i Δ d b β - γ 2 β + i Λ | β | 2 β + g χ a η [ e ( i Δ d a - κ 2 ) t - 1 ] - g 2 ∫ 0 t d τ F ( τ ) β ( t - τ ) , (2) where the linear mode susceptibility χ a - 1 = - i Δ d a + κ / 2 . The first line is the classical equation of motion for a coherently driven Kerr oscillator; the drive term is ∝ g χ a since the linear mode is driven. More interesting is the term in the second line, which describes the delayed self-interaction of the nonlinear mode—mediated by the linear mode—with a memory kernel F ( τ ) = e ( i Δ d a - ( κ / 2 ) ) τ . When F ( τ ) decays rapidly relative to the timescale of system dynamics ( κ ≫ g ), we may set β ( t - τ ) ≈ β ( t ) within a Markov approximation; this is also equivalent to adiabatically eliminating the linear mode ( α ˙ ≈ 0 ). We then obtain an effective Markov regime equation for the (long-time) dynamics of the nonlinear mode, β ˙ = i Δ ˜ d b β - γ ˜ 2 β + i Λ | β | 2 β - η ˜ . (3) This is the classical dynamical equation for a renormalized Kerr oscillator, with modified detuning Δ ˜ d b = ω d - ( ω b + g 2 Δ d a | χ a | 2 ) , damping γ ˜ = γ + g 2 κ | χ a | 2 and drive η ˜ = g χ a η . Therefore, when the linear mode can be adiabatically eliminated, the two-mode system behaves like an effective Kerr oscillator [47] . From here, the classical fixed points ( α ¯ , β ¯ ) of the two-mode system are found by setting α ˙ = β ˙ = 0 , or, equivalently, setting β ˙ = 0 in the Markov regime equation, Eq. (3) . The equation relating the fixed points | β ¯ | 2 to the drive strength | η | 2 is found to be the standard cubic polynomial for a Kerr oscillator, with the modified parameters defined earlier [45] . The relationship is single valued for Δ d b > Δ d b MP but becomes multivalued for Δ d b < Δ d b MP , defining a region of multiple fixed points; here the critical detuning Δ d b MP = - ( 3 / 2 ) ( γ + g 2 κ | χ a | 2 ) + g 2 Δ d a | χ a | 2 [45] . Dropping terms ∝ g 2 arising from the linear mode yields the standard result for a single driven Kerr oscillator. Stability analysis .—To treat the memory term exactly we perform a Laplace domain linear stability analysis around the above fixed points; details can be found in Ref. [45] . Instability is determined by the dominant pole (pole with largest real part) of the linearized dynamical matrix. For a resonantly driven linear mode, Δ d a = 0 , an analysis of the real and imaginary parts of the poles separately allows the phase diagram to be mapped out entirely analytically; we focus on this case from here on (for nonzero Δ d a , see Supplemental Material [45] ). Two distinct parameter regimes are obtained, determined by the relative strength of g and κ . For g < κ / 2 , the typical phase diagram in η - Δ d b space is shown in Fig. 2(a) . For Δ d b > Δ d b MP , the system has only one fixed point (FP), as discussed earlier; the stability analysis indicates that this FP is always stable. For Δ d b < Δ d b MP , the orange hatched region emerges where one of the system’s FPs is unstable, and the unstable pole s has Im s = 0 . In this region, the typical curve relating | β ¯ | 2 to η ( S curve) is shown in the inset, with green (purple) segments showing unstable (stable) FPs. The unstable FPs coincide exactly with the region of 3 total FPs; dynamically, instabilities from the unstable branch settle into one of the two stable fixed points at the same drive strength. This is precisely the stability diagram of the modified Kerr oscillator defined by Eq. (3) . 2 10.1103/PhysRevLett.120.153601.f2 FIG. 2. Phase diagram in η - Δ d b space. The possible phases are listed in the top-left table. (a) Phase diagram for g < ( κ / 2 ) (here, g = ( κ / 4 ) ). (b) Phase diagram for g > ( κ / 2 ) (here, g = 2 κ ). (c) through (f) Plots of | β ¯ | 2 against η showing the change in the S curve as Δ d b is swept across the four dynamical regions in (b). Green (purple) segments of the S curve depict unstable (stable) | β ¯ | 2 values. Much more interesting is the case g > κ / 2 , for which the phase diagram is shown in Fig. 2(b) . We first consider Δ d b > Δ d b MP , where the classical equations admit only 1 FP. We find that for Δ d b above a minimum critical detuning Δ d b L C = - ( 3 / 2 ) ( γ + κ ) , the single fixed point that exists is never unstable (region 1). For Δ d b MP < Δ d b < Δ d b L C (region 2), this is no longer the case. A typical | β ¯ | 2 - η plot in region 2, Fig. 2(c) , shows emergent unstable | β ¯ | 2 values in green, where the system has only one, unstable FP, hinting at the emergence of limit cycle solutions; this regime is not possible for the single coherently driven Kerr oscillator. The minimum and maximum unstable | β ¯ | 2 values occur at drive strengths η - (open square) and η + (solid circle), respectively. At these drives, the dominant pole reaches the threshold of instability, now with nonzero Im s = ± g 2 - κ 2 / 4 ≡ ± Ω . Note that for η > η + and η < η - , the system always has at least one stable fixed point (lying on one of the upward pointing purple segments). Also, in region 2, the unstable green segments lie entirely in the drive range η - < η < η + . However, with more negative detuning, the latter does not remain so. For Δ d b < Δ d b MP (regions 3, 4), the S curve can be multivalued, as seen in Fig. 2(d) , and eventually further deforms to Fig. 2(e) , where unstable | β ¯ | 2 are not all contained in the range η - < η < η + . Here, for η > η + and η < η - the system now admits 3 FPs, of which only one is stable (checkered purple). On the other hand, within the range η - < η < η + (shaded green as before), all three fixed points are unstable . This is clearly seen in Fig. 2(e) : only green segments of the S curve lie in the green shaded region. Furthermore, the range η - < η < η + is detuning dependent; as Δ d b becomes more negative, this region shrinks, and vanishes when η - = η + . We mark this as the terminal boundary of region 3, which occurs at a critical detuning Δ d b crit = - ( Δ d b MP ) 2 + 1 2 [ ( Δ d b MP ) 2 - ( Δ d b L C ) 2 ] [dashed orange line in Fig. 2(b) ]. For Δ d b < Δ d b crit , region 4 begins where η - > η + . The S curve typically looks like Fig. 2(f) . Since at least one stable fixed point always exists for η > η + and η < η - , and since η - > η + in region 4, we easily conclude that at least one stable fixed point (SFP) now exists for all driving strengths. The green shaded region with 0 SFPs thus gives way to the orange hatched region with 3 FPs, 2 SFPs. We note that unlike Δ d b L C , which is a strict minimal detuning for instability, Δ d b crit is not a strict maximal detuning. Beyond Δ d b crit , limit cycle solutions can ostensibly still emerge, since unstable fixed points with nonzero Im s still exist. However, if excursions from these unstable fixed points are large enough, the system can always find a stable fixed point to settle into in this region. Finally, as g → κ / 2 , Δ d b crit , Δ d b MP , Δ d b L C all become equal; both green and purple regions shrink and eventually vanish, such that for g < κ / 2 , only the orange hatched region persists, and we recover the phase diagram in Fig. 2(a) . Dynamics in the unstable region can be studied numerically by simulating Eqs. (1) . We calculate the steady state power spectrum of the linear mode, S α [ ω ] = | F { α ( t ) } | 2 , where F { α ( τ ) } = ∫ - ∞ ∞ d τ e - i ω τ α ( τ ) is the Fourier transform. This quantity is of particular relevance for circuit QED realizations of our model, where S α [ ω ] is the power spectrum of the resonator mode, which can be directly monitored in experiments [48] . For each spectrum S α [ ω ] , the frequency spacing Δ ω is plotted in η - Δ d b space in Fig. 3 , scaled by Ω . We find multifrequency limit cycles in a region that has excellent agreement with the (green shaded) analytic region of 0 SFPs. The spacing Δ ω is close to Ω for Δ d b ∼ Δ d b L C , but decreases as Δ d b becomes more negative; this reduction is observed for general parameters in this system (see additional phase diagrams included in Ref. [45] ). 3 10.1103/PhysRevLett.120.153601.f3 FIG. 3. Numerically simulated phase diagram in η - Δ d b space for g = 2 κ > ( κ / 2 ) . Surface plot shows spectral spacing Δ ω obtained from S α [ ω ] ; the “island” of multifrequency solutions overlaps perfectly with the analytically predicted unstable region below (shaded green). Plots (a) through (c) show S α [ ω ] at the correspondingly labeled points on the phase diagram. Quantum regime .—To study the modification of limit cycle dynamics in the quantum regime, we employ both master equation simulations and a stochastic approach based on the positive- P representation of the density matrix ρ [49] . The latter allows access to normal-ordered operator averages and correlation functions via a set of stochastic differential equations (SDEs) for the independent complex variables α → ≡ ( α , α † , β , β † ) : d α → = A → d t + Λ B d W . → (4) The drift vector A → describes deterministic classical evolution, equivalent to Eqs. (1) . Then, by construction, any quantum effects must appear as stochastic “noise” terms, involving the vector of independent Wiener increments d W → . The scale and nature of this noise is set by the matrix B . In the absence of thermal noise (which we neglect), B = i diag ( 0 , 0 , β , i β † ) is purely quantum in origin, and depends on β , β † (as opposed to being constant for thermal noise). Equations (4) are thus driven by multiplicative noise. Crucially, the classical and quantum contributions depend differently on the nonlinearity. Scaling Λ → Λ / c and ( α → , η ) → c ( α → , η ) in Eqs. (4) leaves the classical drift term unchanged as discussed earlier, but scales the stochastic term by a ( 1 / c ) factor [45] . Hence, decreasing Λ ( c > 1 ; equivalently, increasing mode occupation) suppresses the relative impact of quantum dynamics, moving the system closer to effectively classical evolution. By varying Λ and η according to the above scaling, we can stay on fixed positions on the classical phase diagram (a) just within the limit cycle region, and (b) just outside [see Fig. 4 ], while modifying the stochastic dynamics. For regimes of weak nonlinearity relative to the damping rates, Λ ∼ [ 0.001 , 0.1 ] γ , κ , we find that Eqs. (4) may be simulated controllably; here, modal occupations of O ( [ 10 2 , 10 3 ] ) make master equation and even Monte Carlo simulations unfeasible. For stronger nonlinearities Λ ≳ γ , P representation SDE simulations run into notorious difficulties [50,51] ; however, weaker excitation numbers then render density matrix simulations tractable again [52] . Combining the two methods yields a complete picture of dynamics as the nonlinearity is increased to the quantum regime. The precise scaling of stochastic terms with nonlinearity does depend on the nonlinear model employed [43,53,54] ; more generally, our analysis may be regarded as a study of dynamics under transition from high to low modal occupations. We analyze again the linear mode power spectrum, using the Wiener-Khinchin theorem: S a [ ω ] = ∫ - ∞ ∞ d τ e - i ω τ ⟨ a ^ † ( τ ) a ^ ( 0 ) ⟩ . For SDE simulations, the required correlation function is determined via averaging, ⟨ a ^ † ( τ ) a ^ ( 0 ) ⟩ = lim t → ∞ ( 1 / N s ) ∑ i = 1 N s α i † ( t + τ ) α i ( t ) , summing over at least N s = 10 5 trajectories for each calculation. 4 10.1103/PhysRevLett.120.153601.f4 FIG. 4. | S a [ ω ] | , log scale, (a) within, and (b) outside the limit cycle region, as a function of increasing nonlinearity Λ from left to right. Master equation (ME), SDE simulations (PP), and stable regime linearized spectra (Lin.) are shown. Within the limit cycle region [Fig. 4(a) ], for the weakest nonlinearity Λ = 0.001 γ , the spectrum appears close to the classical result [Fig. 3(c) ]. However, as the nonlinearity becomes stronger, the peaks in the discrete spectrum broaden. A weak-noise phase dynamics analysis [45] indicates a phase diffusion time ∝ 1 / Λ (equivalently, comb peak linewidths ∝ Λ ), with a proportionality coefficient of order one determined by local properties of the limit cycle attractor. Outside the unstable region [Fig. 4(b) ], the classical FP with mode occupations ( | α ¯ | 2 , | β ¯ | 2 ) is stable. Fluctuations around this FP yield a quantum spectrum that we determine analytically [45] using a linearized analysis [47,49,55] ; the result agrees well with S a [ ω ] for weak nonlinearities, but deviates as Λ increases and the fluctuations are no longer small relative to ( | α ¯ | 2 , | β ¯ | 2 ) . For intermediate Λ = 0.1 γ , we are able to compare SDE and master equation simulations in both regions, finding very good agreement. Conclusion .—The driven, strongly coupled nonlinear system of a Kerr oscillator and a linear mode admits a phase with no classical SFPs. Here, classical dynamics features limit cycles with sharp peaks in the mode spectra; however, the quantum dynamics introduce dephasing due to the very nonlinearity that gives birth to the limit cycles, broadening and ultimately washing out these spectral peaks as Λ increases, even if all other noise sources are absent. Our study is relevant for on-chip microwave domain frequency comb generation using quantum circuits with weak nonlinearities (realized in recent circuit QED experiments [18,56,57] ), and for further understanding stable operating regimes of JJ-based nonlinear multimode circuit QED systems. We would like to thank Howard J. Carmichael and Mike Metcalfe for useful discussions. This work was supported by the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Award No. DE-SC0016011. [1] 1 J. Clarke and F. K. Wilhelm , Nature (London) 453 , 1031 ( 2008 ). NATUAS 0028-0836 10.1038/nature07128 [2] 2 P. Bertet , F. R. Ong , M. Boissonneault , A. Bolduc , F. Mallet , A. C. Doherty , A. Blais , D. Vion , and D. Esteve , arXiv:1111.0501 . [3] 3 M. H. Devoret and R. J. Schoelkopf , Science 339 , 1169 ( 2013 ). SCIEAS 0036-8075 10.1126/science.1231930 [4] 4 A. Blais , R.-S. Huang , A. Wallraff , S. M. Girvin , and R. J. Schoelkopf , Phys. Rev. A 69 , 062320 ( 2004 ). PLRAAN 1050-2947 10.1103/PhysRevA.69.062320 [5] 5 A. Wallraff , D. I. Schuster , A. Blais , L. Frunzio , R.-S. Huang , J. Majer , S. Kumar , S. M. Girvin , and R. J. Schoelkopf , Nature (London) 431 , 162 ( 2004 ). NATUAS 0028-0836 10.1038/nature02851 [6] 6 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. H. Devoret , S. M. Girvin , and R. J. Schoelkopf , Nature (London) 449 , 443 ( 2007 ). NATUAS 0028-0836 10.1038/nature06184 [7] 7 L. DiCarlo , J. M. Chow , J. M. Gambetta , L. S. Bishop , B. R. Johnson , D. I. Schuster , J. Majer , A. Blais , L. Frunzio , S. M. Girvin , and R. J. Schoelkopf , Nature (London) 460 , 240 ( 2009 ). NATUAS 0028-0836 10.1038/nature08121 [8] 8 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 ). PRLTAO 0031-9007 10.1103/PhysRevLett.94.123602 [9] 9 J. Gambetta , A. Blais , D. I. Schuster , A. Wallraff , L. Frunzio , J. Majer , M. H. Devoret , S. M. Girvin , and R. J. Schoelkopf , Phys. Rev. A 74 , 042318 ( 2006 ). PLRAAN 1050-2947 10.1103/PhysRevA.74.042318 [10] 10 L. S. Bishop , E. Ginossar , and S. M. Girvin , Phys. Rev. Lett. 105 , 100505 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.105.100505 [11] 11 N. M. Sundaresan , Y. Liu , D. Sadri , L. J. Szőcs , D. L. Underwood , M. Malekakhlagh , H. E. Türeci , and A. A. Houck , Phys. Rev. X 5 , 021035 ( 2015 ). PRXHAE 2160-3308 10.1103/PhysRevX.5.021035 [12] 12 M. D. Reed , L. DiCarlo , B. R. Johnson , L. Sun , D. I. Schuster , L. Frunzio , and R. J. Schoelkopf , Phys. Rev. Lett. 105 , 173601 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.105.173601 [13] 13 R. Vijay , M. H. Devoret , and I. Siddiqi , Rev. Sci. Instrum. 80 , 111101 ( 2009 ). RSINAK 0034-6748 10.1063/1.3224703 [14] 14 F. Mallet , F. R. Ong , A. Palacios-Laloy , F. Nguyen , P. Bertet , D. Vion , and D. Esteve , Nat. Phys. 5 , 791 ( 2009 ). NPAHAX 1745-2473 10.1038/nphys1400 [15] 15 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 ). PLRAAN 0556-2791 10.1103/PhysRevA.39.2519 [16] 16 M. A. Castellanos-Beltran and K. W. Lehnert , Appl. Phys. Lett. 91 , 083509 ( 2007 ). APPLAB 0003-6951 10.1063/1.2773988 [17] 17 B. Abdo , K. Sliwa , L. Frunzio , and M. Devoret , Phys. Rev. X 3 , 031001 ( 2013 ). PRXHAE 2160-3308 10.1103/PhysRevX.3.031001 [18] 18 C. Eichler , Y. Salathe , J. Mlynek , S. Schmidt , and A. Wallraff , Phys. Rev. Lett. 113 , 110502 ( 2014 ). PRLTAO 0031-9007 10.1103/PhysRevLett.113.110502 [19] 19 C. Macklin , K. O’Brien , D. Hover , M. E. Schwartz , V. Bolkhovsky , X. Zhang , W. D. Oliver , and I. Siddiqi , Science 350 , 307 ( 2015 ). SCIEAS 0036-8075 10.1126/science.aaa8525 [20] 20 J. M. Fink , A. Dombi , A. Vukics , A. Wallraff , and P. Domokos , Phys. Rev. X 7 , 011012 ( 2017 ). PRXHAE 2160-3308 10.1103/PhysRevX.7.011012 [21] 21 M. Fitzpatrick , N. M. Sundaresan , A. C. Y. Li , J. Koch , and A. A. Houck , Phys. Rev. X 7 , 011016 ( 2017 ). PRXHAE 2160-3308 10.1103/PhysRevX.7.011016 [22] 22 M. Boissonneault , J. M. Gambetta , and A. Blais , Phys. Rev. Lett. 105 , 100504 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.105.100504 [23] 23 M. I. Dykman and M. A. Krivoglaz , Sov. J. Exper. Theor. Phys. 50 , 30 ( 1979 ). [24] 24 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 ). PRLTAO 0031-9007 10.1103/PhysRevLett.93.207002 [25] 25 I. Siddiqi , R. Vijay , F. Pierre , C. M. Wilson , L. Frunzio , M. Metcalfe , C. Rigetti , R. J. Schoelkopf , M. H. Devoret , D. Vion , and D. Esteve , Phys. Rev. Lett. 94 , 027005 ( 2005 ). PRLTAO 0031-9007 10.1103/PhysRevLett.94.027005 [26] 26 M. Dykman , Fluctuating Nonlinear Oscillators: From Nanomechanics to Quantum Superconducting Circuits ( Oxford University Press , Oxford, 2012 ). [27] 27 P. Del’Haye , A. Schliesser , O. Arcizet , T. Wilken , R. Holzwarth , and T. J. Kippenberg , Nature (London) 450 , 1214 ( 2007 ). NATUAS 0028-0836 10.1038/nature06401 [28] 28 P. Del’Haye , O. Arcizet , A. Schliesser , R. Holzwarth , and T. J. Kippenberg , Phys. Rev. Lett. 101 , 053903 ( 2008 ). PRLTAO 0031-9007 10.1103/PhysRevLett.101.053903 [29] 29 J. S. Levy , A. Gondarenko , M. A. Foster , A. C. Turner-Foster , A. L. Gaeta , and M. Lipson , Nat. Photonics 4 , 37 ( 2010 ). NPAHBY 1749-4885 10.1038/nphoton.2009.259 [30] 30 T. Herr , V. Brasch , J. D. Jost , C. Y. Wang , N. M. Kondratiev , M. L. Gorodetsky , and T. J. Kippenberg , Nat. Photonics 8 , 145 ( 2014 ). NPAHBY 1749-4885 10.1038/nphoton.2013.343 [31] 31 R. P. Erickson , M. R. Vissers , M. Sandberg , S. R. Jefferts , and D. P. Pappas , Phys. Rev. Lett. 113 , 187002 ( 2014 ). PRLTAO 0031-9007 10.1103/PhysRevLett.113.187002 [32] 32 H. A. Haus , IEEE J. Selected Topics Quant. Electron. 6 , 1173 ( 2000 ). 10.1109/2944.902165 [33] 33 J. Faist , G. Villares , G. Scalari , M. Rösch , C. Bonzon , A. Hugi , and M. Beck , Nanophotonics 5 , 272 ( 2016 ). 10.1515/nanoph-2016-0015 [34] 34 Y. K. Chembo and N. Yu , Phys. Rev. A 82 , 033801 ( 2010 ). PLRAAN 1050-2947 10.1103/PhysRevA.82.033801 [35] 35 T. S. Mansuripur , C. Vernet , P. Chevalier , G. Aoust , B. Schwarz , F. Xie , C. Caneau , K. Lascola , C.-e. Zah , D. P. Caffey , T. Day , L. J. Missaggia , M. K. Connors , C. A. Wang , A. Belyanin , and F. Capasso , Phys. Rev. A 94 , 063807 ( 2016 ). PLRAAN 2469-9926 10.1103/PhysRevA.94.063807 [36] 36 A. Ganesan , C. Do , and A. Seshia , Appl. Phys. Lett. 111 , 064101 ( 2017 ). APPLAB 0003-6951 10.1063/1.4985266 [37] 37 A. Ganesan , C. Do , and A. Seshia , Phys. Rev. Lett. 118 , 033903 ( 2017 ). PRLTAO 0031-9007 10.1103/PhysRevLett.118.033903 [38] 38 K. Rayanov , B. L. Altshuler , Y. G. Rubo , and S. Flach , Phys. Rev. Lett. 114 , 193901 ( 2015 ). PRLTAO 0031-9007 10.1103/PhysRevLett.114.193901 [39] 39 B. van der Pol , The London, Edinburgh, and Dublin Philos. Mag. J. Sci. 2 , 978 ( 1926 ). 10.1080/14786442608564127 [40] 40 S. H. Strogatz , Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering ( Avalon Publishing , New York, 1994 ). [41] 41 S. Louca , arXiv:1506.00756 . [42] 42 N. Lörch , J. Qian , A. Clerk , F. Marquardt , and K. Hammerer , Phys. Rev. X 4 , 011015 ( 2014 ). PRXHAE 2160-3308 10.1103/PhysRevX.4.011015 [43] 43 C. Navarrete-Benlloch , T. Weiss , S. Walter , and G. J. de Valcárcel , Phys. Rev. Lett. 119 , 133601 ( 2017 ). PRLTAO 0031-9007 10.1103/PhysRevLett.119.133601 [44] 44 J. Bourassa , F. Beaudoin , J. M. Gambetta , and A. Blais , Phys. Rev. A 86 , 013814 ( 2012 ). PLRAAN 1050-2947 10.1103/PhysRevA.86.013814 [45] 45 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.120.153601 for details of linear stability analysis, weak noise approximations, nonzero Δ d a results, and additional numerical simulations, which includes Ref. [46]. [46] 46 S. M. Girvin , Quantum Machines: Measurement and Control of Engineered Quantum Systems ( Oxford University Press , New York, 2014 ), p. 113 . [47] 47 P. D. Drummond and D. F. Walls , J. Phys. A 13 , 725 ( 1980 ). JPHAC5 0305-4470 10.1088/0305-4470/13/2/034 [48] 48 M. P. da Silva , D. Bozyigit , A. Wallraff , and A. Blais , Phys. Rev. A 82 , 043804 ( 2010 ). PLRAAN 1050-2947 10.1103/PhysRevA.82.043804 [49] 49 H. J. Carmichael , Statistical Methods in Quantum Optics 1—Master Equations and Fokker-Planck Equations ( Springer , New York, 2002 ). [50] 50 A. Gilchrist , C. W. Gardiner , and P. D. Drummond , Phys. Rev. A 55 , 3014 ( 1997 ). PLRAAN 1050-2947 10.1103/PhysRevA.55.3014 [51] 51 C. Gardiner , Quantum Noise—A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods with Applications to Quantum Optics ( Springer , New York, 2004 ). [52] 52 J. R. Johansson , P. D. Nation , and F. Nori , Comput. Phys. Commun. 184 , 1234 ( 2013 ). CPHCBZ 0010-4655 10.1016/j.cpc.2012.11.019 [53] 53 A. Kamal , A. Marblestone , and M. Devoret , Phys. Rev. B 79 , 184301 ( 2009 ). PRBMDO 1098-0121 10.1103/PhysRevB.79.184301 [54] 54 W. Casteels and C. Ciuti , Phys. Rev. A 95 , 013812 ( 2017 ). PLRAAN 2469-9926 10.1103/PhysRevA.95.013812 [55] 55 S. Chaturvedi , C. W. Gardiner , I. S. Matheson , and D. F. Walls , J. Stat. Phys. 17 , 469 ( 1977 ). JSTPBS 0022-4715 10.1007/BF01014350 [56] 56 C. Eichler and A. Wallraff , Euro. Phys. J. Quant. Technol. 1 , 2 ( 2014 ). 10.1140/epjqt2 [57] 57 X. Zhou , V. Schmitt , P. Bertet , D. Vion , W. Wustmann , V. Shumeiko , and D. Esteve , Phys. Rev. B 89 , 214517 ( 2014 ). PRBMDO 1098-0121 10.1103/PhysRevB.89.214517 Publisher Copyright: © 2018 American Physical Society.

PY - 2018/4/9

Y1 - 2018/4/9

N2 - We investigate the dynamics of a microwave-driven Josephson junction capacitively coupled to a lumped-element LC oscillator. In the regime of driving where the Josephson junction can be approximated as a Kerr oscillator, this minimal nonlinear system has been previously shown to exhibit a bistability in phase and amplitude. In the present study, we characterize the full phase diagram and show that besides a parameter regime exhibiting bistability, there is also a regime of self-oscillations characterized by a frequency comb in its spectrum. We discuss the mechanism of comb generation which appears to be different from those studied in microcavity frequency combs and mode-locked lasers. We then address the fate of the comblike spectrum in the regime of strong quantum fluctuations, reached when nonlinearity becomes the dominant scale with respect to dissipation. We find that the nonlinearity responsible for the emergence of the frequency combs also leads to its dephasing, leading to broadening and ultimate disappearance of sharp spectral peaks. Our study explores the fundamental question of the impact of quantum fluctuations for quantum systems which do not possess a stable fixed point in the classical limit.

AB - We investigate the dynamics of a microwave-driven Josephson junction capacitively coupled to a lumped-element LC oscillator. In the regime of driving where the Josephson junction can be approximated as a Kerr oscillator, this minimal nonlinear system has been previously shown to exhibit a bistability in phase and amplitude. In the present study, we characterize the full phase diagram and show that besides a parameter regime exhibiting bistability, there is also a regime of self-oscillations characterized by a frequency comb in its spectrum. We discuss the mechanism of comb generation which appears to be different from those studied in microcavity frequency combs and mode-locked lasers. We then address the fate of the comblike spectrum in the regime of strong quantum fluctuations, reached when nonlinearity becomes the dominant scale with respect to dissipation. We find that the nonlinearity responsible for the emergence of the frequency combs also leads to its dephasing, leading to broadening and ultimate disappearance of sharp spectral peaks. Our study explores the fundamental question of the impact of quantum fluctuations for quantum systems which do not possess a stable fixed point in the classical limit.

UR - http://www.scopus.com/inward/record.url?scp=85045261250&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85045261250&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.120.153601

DO - 10.1103/PhysRevLett.120.153601

M3 - Article

C2 - 29756875

AN - SCOPUS:85045261250

VL - 120

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 15

M1 - 153601

ER -