Quantum Metamaterial for Nondestructive Microwave Photon Counting
Quantum metamaterial for nondestructive microwave photon counting
1,2 1,3 4,5 6,7 6,7 4,5 Arne L. Grimsmo∗ , Baptiste Royer , John Mark Kreikebaum , Yufeng Ye , Kevin O’Brien , Irfan Siddiqi , Alexandre Blais1,8 1 Institut quantique and D´epartment de Physique, Universit´ede Sherbrooke, Sherbrooke, Qu´ebec, Canada J1K 2R1 2 Centre for Engineered Quantum Systems, School of Physics, The University of Sydney, Sydney, NSW 2006, Australia 3 Department of Physics, Yale University, New Haven, CT 06520, USA 4 Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 5 Quantum Nanoelectronics Laboratory, Department of Physics, University of California, Berkeley, California 94720, USA 6 Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 7 Department of Electrical Engineering and Computer Science,Massachusetts Institute of Technology, Cambridge, MA 02139, USA 9 Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 8 Canadian Institute for Advanced Research, Toronto, Canada
Detecting traveling photons is an essential prim- the detection of dark matter axions [20]. itive for many quantum information processing A number of theoretical proposals and experimental tasks. We introduce a single-photon detector de- demonstrations of microwave single-photon detectors have sign operating in the microwave domain, based on emerged recently. These schemes can broadly be divided a weakly nonlinear metamaterial where the nonlin- into two categories: Time-gated schemes where accurate earity is provided by a large number of Josephson information about the photon’s arrival time is needed a junctions. The combination of weak nonlinearity priori [2, 6, 7, 10, 11, 13], and detectors that operate con- and large spatial extent circumvents well-known tinuously in time and attempt to accurately record the obstacles limiting approaches based on a localized photon arrival time [1, 3, 4, 5, 8, 9, 12, 13]. In this work, Kerr medium. Using numerical many-body simu- we are concerned with the last category, which is simul- lations we show that the single-photon detection taneously the most challenging to realize and finds the fidelity increases with the length of the metama- widest range of applications. terial to approach one at experimentally realistic Depending on the intended application, there are sev- lengths. A remarkable feature of the detector is eral metrics characterizing the usefulness of single-photon that the metamaterial approach allows for a large detectors. Not only is high single-photon detection fi- detection bandwidth. In stark contrast to conven- delity required for many quantum information applica- tional photon detectors operating in the optical tions, but large bandwidth, fast detection and short dead domain, the photon is not destroyed by the de- times are also desirable [21]. Moreover, nondestructive tection and the photon wavepacket is minimally photon counting is of fundamental interest and offers new disturbed. The detector design we introduce of- possibilities for quantum measurement and control. In fers new possibilities for quantum information pro- this article, we introduce the Josephson Traveling-Wave cessing, quantum optics and metrology in the mi- Photodetector (JTWPD), a non-destructive single-photon crowave frequency domain. detector which we predict to have remarkably high perfor- mance across the mentioned metrics. In particular, this detector can have detection fidelities approaching unity Introduction without sacrificing detector bandwidth. In contrast to infrared, optical and ultraviolet frequencies The JTWPD exploits a weakly nonlinear, one- dimensional metamaterial, designed to respond to the arXiv:2005.06483v1 [quant-ph] 13 May 2020 where single-photon detectors are a cornerstone of experi- mental quantum optics, the realization of a detector with presence of a single photon. The nonlinearity is provided similar performance at microwave frequencies is far more by a large number of Josephson junctions, inspired by the challenging [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. The Josephson traveling wave parametric amplifier [22]. Be- interest in realizing such a detector is intimately linked cause the detector response does not rely on any resonant to the emergence of engineered quantum systems whose interaction, the detector bandwidth can be designed to natural domain of operations is in the microwaves, includ- range from tens of MHz to the GHz range. The detec- ing superconducting quantum circuits [14], spin ensem- tion and reset times are predicted to be in the range of bles [15], and semiconductor quantum dots [16]. The con- tens of µs for typical parameters. Moreover, the signal- tinuing improvement in coherence and control over these to-noise ratio (SNR) grows linearly with the length of the quantum systems offers a wide range of new applications metamaterial which can be made large, leading to single- for microwave single-photon detection, such as photon- photon detection fidelities approaching unity. By interro- based quantum computing [17], modular quantum com- gating the nonlinear medium with a “giant probe” [23]—a puting architectures [18], high-precision sensing [19], and probe system that couples to the medium over a spatial extent that is large compared to the length of the signal
∗Corresponding author. Email: [email protected] photons—this approach bypasses previous no-go results
1 of coupled LC oscillators, in a configuration known as composite right/left handed (CRLH) metamaterial [28]. The LC oscillators are coupled via an array of nonlinear couplers (inset) to a readout resonator acting as a giant probe (blue). With the metamaterial coupled at x = z/2 to impedance matched linear transmission lines, the inter-± action time between the photon and the giant probe is τ = z/v where v is the speed of light in the metamaterial. Figure 1: a) Sketch of the JTWPD. Standard trans- As an alternative to this transmission mode, the interac- mission lines (black) are coupled to both ends of a one- tion time can be doubled by terminating the metamaterial dimensional metamaterial (orange) of length z and linear at x = +z/2 with an open where the photon wavepacket is dispersion relation, ω = vk. A cross-Kerr interaction χ be- reflected. To simplify the analysis, we consider the trans- tween the metamaterial and the giant probe mode (blue) mission mode in most of the treatment below, but return leads to a phase shift in the strong measurement tone to a discussion of reflection mode when discussing poten- (yellow) while the signal photon (red) travels through the tial experimental implementation and parameters. ˆ metamaterial. b) Phase space picture of the probe mode. The full detector Hamiltonian can be expressed as H = ˆ ˆ ˆ ˆ With respect to the idle coherent state α , the presence of H0 + Hr + Hint, where H0 contains the linear part of the a signal photon displaces the states by |gz/vi , with g = χα. waveguide including the metamaterial as well as the input and output linear waveguides, Hˆr is the probe resonator Hamiltonian and Hˆint describes the nonlinear coupling be- for photon counting based on localized cross-Kerr interac- tween the probe and the metamaterial. As shown in the tions [24, 25, 26, 27]. Supplementary Materials, in the continuum limit where the size a of a unit cell of the metamaterial is small with respect to the extent of the photon wavepacket, Hˆ0 takes Results the form ˆ ˆ ˆ H0 = dω~ωbνω† bνω. (1) Many proposals for itinerant microwave photon detection ν= ZΩ rely on capturing the incoming photon in a localized ab- X± ˆ In this expression, b† ω creates a delocalized right/left- sorber mode that is interrogated using heterodyne de- ± tection [1, 3, 4, 5, 12, 13]. A first challenge associated moving photon with energy ~ω and satisfies the canoni- cal commutation relation [ˆbνω, ˆb† ] = δνµδ(ω ω0). The with this approach is linked to a version of the quantum µω0 − Zeno-effect: continuously and strongly monitoring the ab- subscript Ω in Eq. (1) is used to indicate that we only sorber will prevent the incoming photon from being ab- consider a band of frequencies around which the metama- sorbed [1, 12], limiting the detector’s quantum efficiency. terial’s dispersion relation is approximately linear. The ˆ Another difficulty concerns the tradeoff between efficiency probe resonator Hamiltonian Hr can be written in a dis- and bandwidth. A large detector response to a single pho- placed and rotating frame with respect to the coherent ton requires a sufficiently long interaction time with the drive field as ~K 2 2 photon. In principle, this can be achieved by making the Hˆ 0 = aˆ† aˆ , (2) r 2 absorber mode long-lived. However, as the mode linewidth is inversely proportional to the photon lifetime, this im- where K is a self-Kerr nonlinearity induced by the nonlin- poses a serious constraint on the detector bandwidth. ear couplers (see Methods). Our solution to overcome these obstacles is illustrated The coupling elements also lead to cross-Kerr interac- schematically in Fig. 1: In place of a localized absorber, we tion between the array of oscillators and the probe mode. use a long and weakly nonlinear metamaterial. Backscat- As mentioned above, this coupling is chosen to be locally tering is avoided by using a nonlinearity that is locally weak such that the nonlinearity is only activated by the weak, yet a large response is made possible by having a presence of a strong coherent drivea ˆin(t) on the probe. In ˆ long photon time-of-flight through the metamaterial. The this limit, the nonlinear interaction Hamiltonian Hint is in presence of a photon is recorded using a continuously mon- the same rotating and displaced frame given by itored probe mode that is coupled to the metamaterial z/2 ˆ ˆ ˆ 2 along the full extent of its length. Thanks to a nonlinear Hint0 = ~ dxχ(x)bν† (x)bµ(x) aˆ†aˆ + α z/2 cross-Kerr coupling, in the presence of the measurement νµ Z− X (3) tonea ˆin(t), a single photon in the metamaterial induces a z/2 ˆ ˆ displacement of the output fielda ˆout(t) relative to its idle + ~ dxg(x)bν† (x)bµ(x) aˆ† +a ˆ , z/2 state. While the interaction between the metamaterial νµ Z− X and the probe mode is locally too weak to cause any no- where we have defined the x-dependent photon annihila- ticeable change ina ˆout(t), the displacement accumulates as the photon travels through the metamaterial leading tion operators to a large enough signal to be recorded using homodyne ω¯ dω detection. ˆb (x) = ˆb eνiωx/v, (4) ν 2πv √ω νω r ZΩ JTWPD design and working principle As illus- withω ¯ a nominal center frequency for the incoming photon trated in Fig. 2, the backbone of the metamaterial is a and which is introduced here for later convenience. The waveguide of length z (orange) realized as a linear chain parameter χ(x) is a dispersive shift per unit length given
2 a
Figure 2: Schematic representation of the JTWPD. The probe resonator with ground plane on top and the center conductor below (blue), as well as a readout port on the right, acts as a giant probe. The light blue arrows illustrate the fundamental mode function of a λ/2 resonator. This probe is coupled via a position dependent cross-Kerr interaction χ(x), mediated by an array non-linear couplers (inset), to a metamaterial waveguide (orange). The metamaterial is coupled to impedance matched input/output transmission lines at x = z/2 and x = z/2 (grey). An incoming photon of Gaussian shape ξ(x, t) is illustrated (red). − in Eq. (22), while g(x) = αχ(x) with α the displacement of photon backscattering, is minimized by working with a gi- the probe resonator field under the strong drivea ˆin. The ant probe which, optimally, does not acquire information expression for α, which we take to be real without loss of about the photon’s position. Focusing first on the ideal generality, can be found in Eq. (20) of the Methods. case where the probe mode self-Kerr nonlinearity K and As can be seen from the second term of Eq. (3) which the dispersive shift χ(x) can be neglected compared to dominates for small χ(x) and large α, the combined effect g(x) = αχ(x), we clarify the dominant noise process for of the cross-Kerr coupling and the strong drive results in the probe resonator and the associated backaction on the a longitudinal-like interaction between the metamaterial photon by deriving a perturbative master equation. In and the probe mode [29]. This corresponds to a photon- the subsequent section, we turn to full numerical analysis number dependent displacement of the probe field rela- including the effect of the nonlinearities K and χ(x). tive to the idle state displacement α, which accumulates Considering the ideal case for the moment and ignoring when a photon travels along the metamaterial. By contin- the spatial dependence of g(x), the interaction Hamilto- uously monitoring the output field of the probe mode, a nian takes the simple longitudinal-coupling form photon is registered when the integrated homodyne signal z/2 exceeds a predetermined threshold. This approach shares ˆ ˆ ˆ Hideal = ~g dxbν† (x)bµ(x) aˆ† +a ˆ . (5) similarities with the photodetector design introduced in z/2 νµ Z− Ref. [12], with the important distinction that here the pho- X ton is probed in-flight as it travels through the metamate- We model the incoming photon by an emitter system with rial rather than after interaction with a localized absorber annihilation operatorc ˆ, [ˆc, cˆ†] = 1, located at x0 < z/2 − mode. This distinction is the key to achieving large detec- and initialized in Fock state 1 . The decay rate κc(t) | i tion fidelities without sacrificing bandwidth. of the emitter to the transmission line is chosen such as An important feature of this detector design is that al- to have a Gaussian wavepacket with center frequencyω ¯ though the detection bandwidth is large, the CRLH meta- and full width at half maximum (FWHM) γ propagating material can be engineered such as to have frequency cut- towards the detector [see Eq. (23) of Methods]. Using offs [28]. The low-frequency cutoff avoids the detector Keldysh path integrals, we trace out the waveguide to find from being overwhelmed by low-frequency thermal pho- a perturbative master equation for the joint emitter-probe tons. Decay of the probe mode via the metamaterial to system. As discussed in the Methods, to second order in the input and output waveguides is minimized by choosing the interaction, this master equation takes a remarkably the probe mode resonance frequency to be outside of the simple form metamaterial’s bandwidth. In this situation, the metame- ρˆ˙ = i gn (t)(ˆa +a ˆ†), ρˆ + Γ(t) [ˆa +a ˆ†]ˆρ terial effectively acts as a Purcell filter for the probe mode, − det c D c (6) thereby avoiding degradation of the probe mode quality + κc(t) [ˆc]ˆρ + κa [ˆa]ˆρ. D D factor. Hybridization of the probe resonator and the meta- In this expression, [ˆo] =o ˆ oˆ† 1/2 oˆ†o,ˆ is the usual material is further minimized by using a nonlinear coupler, D • • − { •} illustrated in the inset of Fig. 2. As discussed in Methods, Lindblad-form dissipator and we have definedρ ˆc(t) = cˆρˆ(t)ˆc†/ cˆ†cˆ (t), the coupler is operated at a point where the quadratic cou- h i pling vanishes leaving a quartic potential of strength EQ z/2 1 2 as the dominant contribution. ndet(t) = dx ξ (x, t) , (7) v z/2 | | Z− 2 z/2 4g κa x+ z 2 Backaction and detector noise In the JTWPD, back- Γ(t) = dx 1 e− 2v ( 2 ) ξ (x, t) , (8) action on the incoming photon’s wavevector, and therefore κav z/2 − | | Z− h i
3 with ξ(x, t) = ξ(t x/v) the incoming photon envelop .0050 − and ndet(t) the fraction of the photon that is in the meta- material at time t. A term of order g/ω¯ describing back- .0025 scattering of the photon into the left-moving field has been dropped from Eq. (6). Withω ¯ the carrier frequency of the 0 100 200 300 400 500 600 incoming photon, this contribution is negligible. MPS site In Eq. (6),ρ ˆc is the state of the system conditioned on a photon having been emitted. The first term of Eq. (6) con- sequently has an intuitive interpretation that is consistent with the form of Hˆideal: The probe resonator is condition- ally displaced by a drive equal to the longitudinal coupling amplitude times the photon fraction in the metamaterial, g n (t). Indeed, while the x-quadrature of the probe, × det xˆ = (ˆa† +a ˆ)/√2, is a constant of motion under Eq. (6), Figure 3: The top panel shows snapshots of the photon the y-quadrature,y ˆ = i(ˆa† aˆ)/√2, is displaced. − The second term of Eq. (6), proportional to the rate number population along the MPS sites at three different Γ(t), is the dominant process contributing to noise also times t1 (red) < t2 (green) < t3 (blue). The white re- along the y-quadrature. The origin of the noise term can gion corresponds to the linear waveguide and the orange be understood as follows. When the photon first enters region to the metamaterial with its coupling to the probe the detector and is only partially inside the metamaterial, resonator. The lower three panels show the Wigner func- the probe mode field evolves to a superposition of being tion W (x, y) of the intracavity probe field at the three displaced to different average values ofy ˆ, leading to en- respective times. When the photon is only partially in- hanced fluctuations in this quadrature. This effect can side the metamaterial, the probe is in a superposition be seen clearly in the numerical results of Fig. 3, which of displaced states (middle lower panel). Parameters are are described in more detail below. Finally, the last line κa = χ(x) = K = 0, gτ = 2 and γτ = 2. of Eq. (6) describes the usual decay of the emitter and probe at respective rates κc(t) and κa. state as a quantum trajectory conditioned on the mea- As the increased fluctuations in the y-quadrature arise surement record [32]. With our approach this is simu- due to uncertainty in the photon’s position, a spatially lated using a stochastic MPS algorithm. Further details longer photon is expected to lead to larger fluctuations. on this numerical technique can be found in Methods and A measurement of the probe’s y-quadrature will collapse the Supplementary Materials. the superposition of displaced states and thus lead to a As in the previous section, we consider a Gaussian pho- backaction effect localizing the photon and randomizing ton wavepacket with FWHM γ propagating towards the its wavevector. This effect can be minimized by decreasing detector by an emitter initialized in the state 1 local- the interaction strength g while keeping gz/v constant by ized to the left of the detector. The interaction strength| i is increasing z. In other words, backaction can be minimized quantified by the dimensionless quantity gτ where τ = z/v by increasing the detector length relative to the spatial is the interaction time as before, and the photon width by extent of the photon. This intuitive reasoning is confirmed the dimensionless quantity γτ. Example snapshots of the by numerical results in the next section. photon number distribution along the MPS sites at three different times t1 < t2 < t3 are shown in Fig. 3, along Numerical Matrix Product State simulations We with the corresponding Wigner functions of the probe now turn to numerical simulations of the JTWPD includ- mode field. Because of the impedance match and negli- ing the self- and cross-Kerr nonlinearities K and χ that gible backaction, the photon wavepacket travels without were dropped from the above discussion. To go beyond the deformation along the waveguide. perturbative results of Eq. (6), it is no longer possible to We start by comparing numerical results from MPS integrate out the waveguide degrees of freedom. A brute- simulations to the perturbative master equation obtained force numerical integration of the dynamics is, however, in Eq. (6). To help in directly comparing the simula- intractable, as the JTWPD is an open quantum many- tion results, we first consider the idealized situation where body system with thousands of modes. We overcome this χ(x) = K = 0. In Fig. 4, we show the average probe res- obstacle by using a numerical approach where the systems onator displacement yˆ whose integrated value is linked is represented as a stochastically evolving Matrix Product to the detector signalh andi the noise ∆ˆy2 as a function State (MPS) conditioned on the homodyne measurement of time. To verify the prediction thath fluctuationsi iny ˆ record of the probe output field. increase for spatially longer photons, we compare Gaus- Our approach is based on trotterizing the time evolu- sian wavepackets of different spectral widths γ. Recall tion and discretizing the photon waveguide, including the that a smaller γτ implies a longer photon relative to the nonlinear metamaterial, along the x axis. Building upon detector length. The solid lines in Fig. 4 are obtained us- and extending recent developments of MPS in the con- ing MPS simulations with γτ = 2 (blue), 4 (orange), 6 text of waveguide QED [30, 31], this leads to a picture (green) and 10 (bright purple). The dotted lines are ob- where the waveguide is represented by a “conveyor belt” tained from Eq. (6) for the same parameters. The agree- of harmonic oscillators (referred to as MPS sites below) ment between the approximate analytical results and the interacting with the probe resonator (see Methods). Mea- full non-perturbative MPS results is remarkable. surement backaction under continuous homodyne detec- In panels (c, d) of Fig. 4 we use a spatially varying g(x), tion of the probe resonator is included by representing the and we consequently only show MPS results in these pan-
4 3
1.5 2
2 ˆ y ˆ y
1.0
1 ∆
0 0.5 signal 0 2 0 2 3 0.7 filtere d 2
2 ˆ y ˆ y 0.6
1 ∆ 0 0.5 0 1 2 3 4 5 6 0 2 0 2 τm/τ 1 2
1
ˆ 0.6 y 10− ˆ y ∆
0 0.5 F 2 10− 0.0 2.5 5.0 0.0 2.5 5.0 − 1 t/τ t/τ 3 10− Figure 4: Time evolution of the intra-cavity probe dis- placement yˆ [(a, c, e)] and fluctuations ∆ˆy2 [(b, d, f)], 1.0 1.5 2.0 2.5 3.0 h i h i gτ in the idealized case χ(x) = K = 0. Top row: κa = 0 and gτ = 2. Middle row: κa = 0 and spatially varying g(x) with average valuegτ ¯ = 2. Bottom row: κaτ = 1.0 and Figure 5: (a) 75 filtered homodyne currents (arbitrary 2 units) for gτ = 3, κ τ = 1.0, K /κ = 10− and g/χ = 5. gτ = 2. The solid lines correspond to MPS simulations a | | a with different photon widths γτ = 2 (blue), 4 (orange), 6 Red traces are obtained with an incoming Gaussian pho- (green) and 10 (bright purple), while the dotted lines are ton of unitless width γτ = 6, and blue traces for vacuum. from integrating Eq. (6). The horizontal gray line is the threshold chosen to max- imize the assignment fidelity. (b) Infidelity versus gτ for γτ = 2 (blue), 4 (orange), and 6 (green), found by aver- els. In practice, the probe will be realized from a resonator aging over Ntraj = 2000 trajectories. Other parameters whose vacuum fluctuations vary in space. To confirm the as in (a). The shaded regions indicate the standard error defined as (1 )/N . robustness of the detector to this variation, Fig. 4 (b, c) ± F − F traj 2 shows yˆ and ∆ˆy versus time as obtained from MPS p simulationsh i for hg(x)i = 2¯g cos2(2πx/z) + µ(x). The cosine is larger than a threshold ythr, i.e. maxt J¯hom(t) > ythr. models the dependence on the mode function of a λ/2 res- The filter f(t) yˆ(t) is obtained from averaging over onator while µ(x) is added to take into account potential a large number∝ of h trajectoriesi and is chosen such as to random variations in the coupling strength which we take give more weight to times where the signal is on average here to be as large as 10%. Moreover, to show the effect larger. We maximize t over the time window [ τm, τm] of a non-uniform g(x) more clearly, we use γτ = 10 cor- and chose the threshold to optimize between quantum− ef- responding to spatially shorter photons than in the other ficiency and dark counts. The quantum efficiency η is de- panels. Although additional structures can now be seen, fined as the probability of detecting a photon given that the long-time average displacement remains unchanged one was present. From the above procedure, it can be es- confirming that the detector is robust against spatial vari- timated as η = Nclick 1/Ntraj 1, with Nclick 1 the number ations of the metamaterial-probe coupling. | | | of reported “clicks” and Ntraj 1 the number of simulated | Panels (e, f) of Fig. 4 show results fo κa > 0. In this sit- trajectories with a photon. On the other hand, the dark uation the MPS evolves stochastically with each trajectory count probability is estimated similarly as the fraction of resulting in a measured current Jhom(t) = √κa yˆ traj + reported clicks pD = Nclick 0/Ntraj 0 in a simulation with h i | | ξ(t), where ξ(t) = dWt/dt with dWt a Wiener process no incoming photon. In these simulations, the dark count representing white noise [32]. We compare yˆ and ∆ˆy2 rate is set by the threshold and the vacuum fluctuations h i h i averaged over one thousand stochastic trajectories to the of the probe resonator. A number that incorporates both results obtained by integrating the Keldysh master equa- η and pD, and is thus a good measure of the performance tion Eq. (6). The agreement is excellent for large γτ, but of a photodetector, is the assignment fidelity [4] small deviations are observed when this parameter is de- creased. We attribute this to terms of higher than second 1 = (η + 1 pD) . (10) order in the interaction Hamiltonian, which are neglected F 2 − in Eq. (6). The exponential decay of yˆ at long time ob- h i In practice, if the arrival time of the photon is known to lie served in panel (e) simply results from the finite damping within some time window, one can optimize t in Eq. (9) rate κa. Indeed, the photon-induced displacement stops over this window in a post-processing step [12]. In our once the photon has travelled past the metamaterial at numerical simulations, the arrival time is known such that which point the probe mode relaxes back to its idle state. this optimization is not necessary and we can therefore For a given trajectory, we infer that a photon is detected simply evaluate J¯hom(t) at t = 0. if the homodyne current convolved with a filter [4] Fig. 5 shows 75 typical filtered output records, J¯hom(t = 0), as a function of the measurement window τm. These τm results are obtained from stochastic MPS simulations with J¯hom(t) = dt0Jhom(t0)f(t0 t). (9) − γτ = 6, gτ = 3, κaτ = 1.0, and include self- and cross- Z0
5 2 Kerr couplings with K /κ = 10− and g/χ = 5. The | | a red traces correspond to simulations where a photon was 4000 present, while the blue traces are for incoming vacuum. The horizontal gray line is the threshold chosen to opti- 2000 mize the assignment fidelity. At τm/τ & 3, most traces are correctly identified. Panel (b) shows the assignment fidelity for γτ = 2 (blue), 4 (orange) and 6 (green) as as 0 function of gτ but fixed g/χ = 5. The measurement time 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 τm is chosen sufficiently large to maximize . As expected [2π MHz] from Fig. 4, the fidelity is reduced for smallerF γτ because spatially longer photons (smaller γτ) lead to more noise Figure 6: Number of unit cells needed to reach gτ in the in the measurement. range 1–3 as a function of self-Kerr non-linearity KQ, for A remarkable feature of Fig. 5 is the clear trend of the α = 5,ω/ ¯ (2π) = 5 GHz, Is = EQ/ϕ0 = 1.1 µA and assignment fidelity approaching unity with increasing gτ. Ztml = 50 Ω. The total Kerr non-linearity of the resonator This number can be increased at fixed interaction strength K = KQ + KJ can be tuned close to zero by introducing g by increasing the detector length. In the next section we another non-linearity with KJ < 0. show that values of gτ in the range 1–3 used in Fig. 5 are within reach for experimentally realistic parameters and metamaterial lengths. KJ < 0. Following this approach, we can allow for a de- tector with a larger KQ contributing to reducing Ncells, Towards experimental realization The JTWPD yet still have a total Kerr nonlinearity K 0 to avoid shares similarities with the Josephson Traveling Wave nonlinear response of the probe mode. Similar' ideas have Parametric Amplifier (JTWPA) [22, 33, 34]. State of the recently been used to cancel unwanted cross-Kerr nonlin- art JTWPAs consists of a metamaterial with up to tens of earities [36]. thousands of unit cells, each comprised of a large Joseph- Fig. 6 shows Ncells as a function of the self-Kerr KQ son junction and a shunt capacitance to ground. In addi- to reach gτ in the range 1–3, for a photon center fre- tion, LC oscillators used to engineer the dispersion rela- quency ofω/ ¯ (2π) = 5 GHz. In these plots we use a tion are placed every few unit cells. We envision a JTWPD nonlinearity Is = EQ/ϕ0 = 1.1 µA for the coupling ele- with a similar number of unit cells, albeit with an increase ment, c.f. Methods, and the other parameters are α = 5 in complexity for each unit cell. A significant design dif- and Ztml = 50 Ω. Crucially, it is possible to reach gτ in ference is that in the JTWPD every unit cell is coupled to the range 1–3, as in our numerical simulations above, us- the same probe resonator. In practice, this resonator can ing a few thousand unit cells without an excessively large be a coplanar waveguide resonator or a 3D cavity. KQ. Alternatively, the same value of gτ can be reached As shown in the Supplementary Materials, the number for a smaller KQ by increasing the transmission line char- of unit cells necessary to reach a given value of gτ can be acteristic impedance, Ztml, as is clear from Eq. (11). As approximated by discussed in more detail in the Supplementary Materials, K can be tuned by varying the coupling capacitance be- 1 gτ R 2 ω¯2 Q N K . (11) tween the junctions and the probe resonator, or by tuning cells ' 2 α 8πZ K E / tml Q Q ~ the characteristic impedance of the coupler mode. where we neglect spatial dependence of the parameters for The CRLH metamaterial has a frequency-independent simplicity. In contrast to the simulation results presented characteristic impdeance Ztml = Ln/Cg given that L /C = L /C , referred to as a balanced above, we assume here that the detector is operated in n g g n p reflection mode, effectively halving the number of unit cells CRLH [28]. Close to the center of the CLRH fre- p p needed for a given value of τ. In this expression, α is quency band, the dispersion relation is approximately lin- the displacement of the probe resonator as before, RK = ear, with a speed of light given by v = 1/ 4LnCg. h/e2 is the quantum of resistance, Z the characteristic For typical parameters, discussed in more detail in the tml p impedance of the metamaterial at the center frequency Supplementary Materials, we expect detection times in ω¯, and EQ the nonlinear energy of the coupling elements, the range τ = 1–10 µs. To have κaτ = 1 as in the sim- discussed in more detail in Methods. ulations above, this then suggests a probe decay rate in The parameter K appearing in Eq. (11) is the self- the range κa/2π 0.015–0.15 MHz. Larger values of κaτ Q ' Kerr nonlinearity of the resonator [see Eq. (14)] due to might be preferable in practice, but we found this regime the nonlinear couplers in Fig. 2. An interesting feature too demanding for numerical simulations due to the pro- of the the coupling element we make use of is that the hibitively small time steps needed. A larger κa relaxes the self-Kerr is always positive K > 0, in contrast to a more constraint on reducing the total self-Kerr nonlinearity K . Q | | conventional Josephson junction element [35]. The total Based on the numerical results in the previous section, Kerr non-linearity of the resonator can be adjusted by in- the detection time is of the order τm 3τ, and thus ex- troducing another nonlinear element such as one or more pected to be in the µs to tens of µs range' for the above Josephson junctions galvanically or capacitively coupled value of τ. The detector reset time is naturally of the to the resonator. We can then write the total Kerr non- order 1/κa, but can likely be made faster using active re- linearity as K = KQ + KJ , where KQ > 0 is the contribu- set protocols. To avoid significant backaction effects, the tion from the couplers in Fig. 2, and KJ < 0 comes from photon’s spectral width must not be too small as we have one or more Josephson junctions. The latter elements can shown in the preceding sections. A value for the dimen- moreover be made tunable, allowing an in-situ tuning of sionless photon width of γτ = 2 corresponds to a FWHM
6 of γ/(2π) = 0.25 MHz, for the value τ = 1 µs. We empha- the shape of the photon population wavepacket is min- size that the detection fidelity increases with increasing γ, imally disturbed by the detection. Together with the and from our numerical results we thus expect photons of large bandwith and high detection fidelity, this opens new spectral width in the MHz range or larger to be detectable possibilities for single-photon measurement and control, with very high fidelity. including feedback of photons after measurement, weak The bandwidth of the detector is set by appropriately single-photon measurement, and cascading photon detec- choosing the parameters of the CRLH metamaterial. In tion with other measurement schemes or coherent interac- the Supplementary Materials we show example parameter tions. sets with bandwidths ranging from several GHz to 100s of MHz. For some applications that require very low dark count rates, lowering the bandwidth might be desirable. Methods In principle the CLRH bandwidth can be made arbitrar- Nonlinear coupling element We make use of a circuit ily small, but the circuit parameters required might be- identical to the SNAIL element introduced in Ref. [37], but come challenging to realize. Another option is to replace used at a different operating point. The coupler consists the coupling element shown in the inset of Fig. 2 by a of a loop of n large junctions with Josephson energy E floating coupler, such that the bandwidth is controlled by s J and a single smaller junction with energy βE , leading to a coupling capacitance. All of these various options are J a nonlinear potential discussed in more detail in the Supplementary Materials. ϕˆ UˆQ(ϕ ˆ) = βEJ cos(ϕ ˆ ϕx) nsEJ cos , (12) − − − ns Discussion where ϕx is the dimensionless flux encircled by the loop. The coupler is operated at the point ϕ = π and β = 1/n Previous work have questioned whether cross-Kerr inter- x s where the potential becomes action can be used for high-fidelity single photon count- ing [27], seemingly in contradiction with our results. There EQ Uˆ (ϕ ˆ) = ϕˆ4 + ..., (13) is, however, a fundamental difference between our proposal Q 24 and the approach of Ref. [27]. There, a number of non- and we have introduced E = E (n2 1)/n3. Here we linear absorbers independently couple to a traveling con- Q J s s have expanded the nonlinear potential− aroundϕ ˆ 0, trol field. This is similar to an alternative version of our which is valid based on the fact that each end of the' ele- proposal where each unit cell of the metematerial couples ment is coupled to harmonic modes with small zero-point to an independent probe resonator. More generally, we flux fluctuations. The crucial property of this coupler is can consider a situation where we partition the N unit cells that it provides a purely nonlinear quartic potential while cells of the detector into M blocks, with each block cou- the quadratic contribution cancels out. This minimizes pled to an independent readout probe resonator. With hybridization between the metamaterial and the probe M = N we have a setup similar to Ref. [27], while cells resonator in the JTWPD, and is a very useful tool for M = 1 corresponds to the JTWPD. However, as shown generating non-linear interaction in general [38]. In prac- in Methods, such a setup gives a √M reduction in the tice there will be small deviations from the ideal opera- probe resonator’s displacement. Our proposal thus has an tion point φ = π, β = 1/n , but as we shown in the √N improvement in the SNR scaling. This improve- x s cells Supplementary Materials, the JTWPD is robust to such ment comes from using what we referred to in the intro- imperfections. duction as a giant probe, i.e. a probe resonator that has a The positive quartic potential in Eq. (13) leads significant length compared to the photon. This contrasts to positive self- and cross-Kerr nonlinearities for the with conventional circuit QED-based photodectors relying probe-metamaterial system, in contrast to more con- on point-like probe systems. Such a setup does not have ventional Josephson junction nonlinearities. In the any obvious analog in the optical domain, demonstrating Supplementary Materials we use a black-box quantization the potential of using metamaterials based on supercon- approach [35] to estimate the Kerr nonlinearities. In par- ducting quantum circuits to explore fundamentally new ticular, the self-Kerr nonlinearity of the probe mode in- domains of quantum optics. duced by Ncells coupler elements takes the form In summary, we have introduced the JTWPD, a mi- crowave single-photon detector based on a weakly nonlin- Ncells 4 ~KQ = EQ,n ϕr(xn) , (14) ear metamaterial coupled to a giant probe. This detector | | n=0 is unconditional in the sense that no apriori information X about the photon arrival time or detailed knowledge of with EQ,n the energy of the nth nonlinear coupler and the photon shape is needed for its operation. Detection fi- ϕr(xn) the dimensionless zero-point flux fluctuations of delities approaching unity are predicted for metamaterial the probe mode biasing the nth coupling element. length that are compatible with state-of-the-art experi- ments. Moreover, because the JTWPD does not rely on Dynamics of the probe resonator The probe res- absorption in a resonant mode, large detection bandwidths onator Hamiltonian can be written as are possible. K ˆ ~ 2 2 iωdt A remarkable feature of the JTWPD, which distin- Hr = ~ωraˆ†aˆ + aˆ† aˆ + ~ iεe− aˆ† + H.c. , (15) guishes this detector from photodetectors operating in the 2 optical regime, is the nondestructive nature of the in- witha ˆ the annihilation operator