Integrated conversion and photodetection of virtual photons in an ultrastrongly coupled superconducting quantum circuit

The ground-state of an artificial atom ultrastrongly coupled to quantized modes is entangled and contains an arbitrary number of virtual photons. The problem of their detection has been raised since the very birth of the field but despite the theoretical efforts still awaits experimental demonstration. Recently experimental problems have been addressed in detail showing that they can be overcome by combining an unconventional design of the artificial atom with advanced coherent control. In this work we study a simple scheme of control-integrated continuous measurement which makes remarkably favourable the tradeoff between measurement efficiency and backaction showing that the unambiguous detection of virtual photons can be achieved within state-of-the art quantum technologies.

Integrated conversion and photodetection of virtual photons in an ultrastrongly coupled superconducting quantum circuit.

Introduction
Advancement of fabrication technologies has allowed producing solid-state systems which exhibit the physics of atoms ultrastrongly coupled (USC) to quantized modes [1][2][3][4].In these systems, the coupling constant g between the artificial atom (AA) and the mode is comparable or even larger than the bare excitation frequencies ϵ and ω c of the AA and of the mode.This regime has been achieved on several different platforms [2,3], semiconductors [5][6][7] and superconductors [8][9][10] being the most promising for applications.
In the USC regime, nonperturbative physics is predicted to emerge [2,3] which is undetectable in the standard strong-coupling regime of quantum optics [11,12].Higher-order antiresonant terms in the Hamiltonian break the conservation of the number of excitations this occurrence being at the heart of most of the phenomenology USC is expected to exhibit.A striking feature is that eigenstates are highly entangled, and in particular the ground-state |Φ⟩ contains virtual photons (VP).The simplest instance is the two-level quantum Rabi model [13,14] where {|g⟩, |e⟩} are the eigenstates of a two-level atom and a (a † ) are the annihilation (creation) operators of a quantized harmonic mode.The ground state has the form {|n⟩} being the number eigenstates of the mode [2,3].It is seen that |Φ⟩ contains an even number N = n + |e⟩⟨e| of excitations while in the absence of USC the ground state is |0g⟩ thus it is factorized and it does not contain VPs.The question of how to detect ground-state VPs in USC systems has been posed since the birth of the field [1].They cannot be probed by standard photodetection because the USC vacuum |Φ⟩ cannot radiate [15], thus VPs must be converted to real excitations which are then detectable.Early theoretical proposals of VP detection leverage time-dependent coupling constants as in the dynamical Casimir effect [16][17][18] but require modulation of quantum hardware at subnanosecond times which is still unavailable.Another class of proposals formulated in the last decade [19][20][21][22][23] introduces an additional lower energy AA level |u⟩ not coupled to the mode (see Fig. 1a) making |Φ⟩ a false vacuum which can undergo radiative decay.This option also poses several experimental challenges, from the low yield of detectable photons to the fact that conventional quantum hardware does not ensure that the conversion is faithful, i.e. that output photons are really produced due solely to USC [23].For these reasons, despite the huge theoretical effort, detection of VPs still awaits demonstration.In a recent work [24], it has been shown that the above experimental problems can be overcome.Efficient, faithful and selective conversion of VPs to real ones can be achieved by combining an unconventional superconducting multilevel AA design [25][26][27], with coherent amplification of the conversion of ground-state VPs by advanced control [21] and a tailored measurement protocol.
In this work, we discuss a toy model of integrated protocol combining VP conversion by STIRAP [24] with photodetection by continuous measurement of the mode [28].In particular, we consider decay into a transmission line coupled to the mode during the whole protocol, which is the simplest experimental option.A detectable signal is obtained if the decay rate κ of the mode is large enough which however, determines a backaction of the continuous measurement inducing decoherence which may affect the efficiency of the coherently amplified conversion.Our results show that STIRAP is resilient to this backaction making favourable the tradeoff of the integrated protocol with a continuous measurement.

Model
We illustrate VPs conversion/detection considering the Hamiltonian of a three-level AA coupled to the mode (see Fig. 1a) (3) Hamiltonian (3) describes a three-level AA with the two excited states |g⟩ and |e⟩ ultrastrongly coupled to the mode, as described in Eq. (1), and the ground state |u⟩ uncoupled.It is a three-level approximation of the multilevel Hamiltonian described in ref. [24], which is implemented by a fluxonium-like superconducting AA galvanically coupled to a mode also implemented by a superconducting LC resonator.Parameters for the lumped elements of this quantum circuit have been found such that coupling of |u⟩ to the mode is very small thus guaranteeing that only VPs are converted in real photons.For instance, the first level splitting ϵ ′ > 0 of the "uncoupled" state |u⟩ is much larger than ω c = ϵ, the latter being the second atomic splitting.For a detailed discussion we refer to ref. [24].The eigenstates of H are classified in two sets (see Fig. 1a), namely the factorized states {|n⟩ ⊗ |u⟩} with energies −ϵ ′ + nω c , and the entangled Rabi-like eigenstates {|Φ l ⟩} of H R , with eigenvalues E l .

Conversion of VPs
Conversion of VPs employs a STIRAP protocol [29] where the system is driven by a two-tone field W (t) = W s (t) cos ω s t + W p (t) cos ω p t mainly coupled with the u − g transition of the AA.We take the field resonant with the two relevant transitions Standard approximations yield the Λ driving configuration [30] of Fig. 1a described in a rotating frame by the control Hamiltonian [21] where the Rabi amplitudes Ω p (t) = W p (t) γ ug ⟨0g|Φ⟩ and Ω s (t) = W s (t) γ ug ⟨2g|Φ⟩ depend on the matrix element γ ug := ⟨u|γ|g⟩ of the AA "dipole" operator [24].Operating the "counterintuitive" pulse sequence [29] 1b, population transfer by STIRAP occurring in the first part of the protocol.As shown in Ref. [24] the nearly ideal scenario described in this section can be implemented also in the USC regime by state-of-the-art superconducting quantum technologies in an unconventional design of the quantum circuit, and using superinductors [25][26][27] and advanced control at microwave frequencies.Control based on STIRAP has been proposed [33][34][35][36] and demonstrated [37,38] in standard superconducting quantum devices.

Toy-model for photodetection
Ideally, once the population has been transferred in |2u⟩ the converted VP pair can be detected.For oscillators with quality factor Q ≳ 10 4 the population of the mode remains large enough [39] to allow photons to be detected (and even counted) by single-shot nondemolition measurements performed by a quantum probe coupled dispersively to the mode [40,41].A much simpler procedure is a continuous measurement [28] which uses radiative decay of converted VPs with a rate κ into a transmission line.In this case, a key advantage is that the initial state is faithfully prepared by simply letting the system relax [42] due to photodetection.Thus the protocol can be repeated over and over yielding a detectable signal if κ is large enough.Since the oscillator selection rules prevent direct |2u⟩ → |0u⟩ decay, photodetection involves the sequential decay of the mode |2u⟩ → |1u⟩ → |0u⟩.Therefore we formulate a minimal model restricting the analysis to the four-dimensional Hilbert space spanned by {|0u⟩, |1u⟩, |2u⟩, |Φ⟩}.In Fig. 1b we show the population histories with an oscillator decay rate κ(t) turned on after the completion of STIRAP.Photons in |2u⟩ first decay to |1u⟩ and then to |0u⟩ at each step a photon being emitted into the transmission line.A minimal model of decay is described by a Lindblad equation with Lindbladian given by where the dissipator is defined as Â † Â ρ + ρ Â † Â and ρ is the density matrix of the system.The first term describes emission with decay rate κ and the second describes absorption whose rate has been written making the phenomenological assumption that detailed balance at thermal equilibrium, β = 1/(k B T ), can be used also for the driven system.
We point out that some care is required to guarantee a physically consistent picture of photon decay.We now briefly describe how to interpret Eq.( 5) to obtain the minimal model of photodetection.The key point is that the operator â must be defined such to avoid photon annihilation in the Rabi ground state |Φ⟩ otherwise we could have photon emission even in the absence of the level |u⟩, which is unphysical.The complete theory requires using "dressed" field operators, say â → Â [3].In our  1.Parameters used in the simulations and corresponding physical values for the superconducting flux-based architecture including a superconducting junction with Josephson energy EJ /(2π) = 10 GHz considered in Ref. [24], where an extended Rabi model was studied with coupling constant g/ωc = 0.5 yielding squared matrix elements |⟨0g|Φ⟩| 2 = 0.42 and |⟨2g|Φ⟩| 2 = 0.05 for the probability of finding respectively n = 0, 2 virtual photons in the ground state.
case, this brings a simplification since the new operators are such that Â|Φ⟩ = 0 and they reduce to â when acting on |nu⟩.Since truncation to the four-level space also implies that Â † |Φ⟩ = 0 we simply have to use a projected version of the bare jump operators acting only on states |nu⟩.This provides the correct minimal description of both decoherence and photodetection.
The measured quantity is the extra current of photons emitted into the transmission line defined as where the first term is the average total current for the driven system where t M is the duration of the whole integrated conversion and measurement protocol and the second term is the equilibrium thermal current in the undriven system.Fig. 1b also shows the number of emitted photon pairs (black dashed line) in a cycle at zero temperature when the thermal current is also zero.The parameters used are such that the VPs conversion is complete and the system is reset to the initial state |0u⟩.

Atomic decay
In principle, atomic decay is not relevant in the ideal protocol since in STIRAP the intermediate state |Φ⟩ is expected not to be populated (see Fig. 1b).However, we will see that AA decay plays a role in the integrated conversion/measurement protocol (see Fig. 2a).Thus we take into account it in the Lindblad formalism by introducing two dissipators with jump operators |0, u⟩⟨Φ| and |2u⟩⟨Φ|.Nonradiative decay rates could be explicitly calculated by the Fermi Golden Rule if the atomic environment were specified.Even in the absence of detailed information, we know that the rates |Φ⟩ → |nu⟩ are proportional to the square of the matrix elements ⟨nu| 1 osc ⊗ |u⟩⟨g| |Φ⟩ = ⟨ng|Φ⟩, thus once again we express them in terms of a single parameter γ as The excitation rates at equilibrium are written assuming again that detailed balance holds and finally, we obtain the Lindbladian 3

Results and conclusions
The main result of this work is shown in Fig. 2a where we plot population histories at finite T = 50 mK for an always-on detector, i.e. taking κ constant in our equations.
Having in mind repeating the cycle over and over in order to collect a large enough signal for detection, we show a cycle starting from the thermal state of the system (horizontal dashed lines) that for the parameters chosen also ends in the same state, being thus a limiting cycle.Thermal effects reduce the net population transferred, thus the system converts pairs of VPs with a smaller probability.The backaction of the always-on detector is expected to further reduce population transfer and VP conversion due to dephasing during adiabatic passage in STIRAP and to decay of the mode when |2u⟩ starts to be populated.The former effect leads to a reduction of the final population of |2u⟩ after the completion of STIRAP estimated by [43] where in our case κ ϕ = 3κ/2 [24], which for the value of κ in Table1 turns out to be small.On the contrary, decay of the mode after the adiabatic passage phases has a significant impact on ideal STIRAP since it determines in a strong reduction of P 2 (see Fig. 2a).However, this population loss results in the detection of converted VPs photon pairs which progressively populate |2u⟩ during STIRAP thus the probability of detecting a photon pair per cycle remains large.Notice that the total number of photons decaying in the transmission line (black dotted line in Fig. 2a) is larger than two and increases linearly at very short and large times.This linear component is due to the constant current of thermal photons κ⟨a † a⟩ th which has nothing to do with VPs.By subtracting the thermal part, we obtain the number of detected VP pairs, which turns out to be almost equal to the thermal population of the initial state |0u⟩ (gray dotted line in Fig. 2a).
Comparing the continuous measurement of Fig. 2a with the switchable probe protocol of Fig. 1b we notice that photodetection during the protocol strongly modifies population histories.However, coherent amplification of the VP is preserved since the population of |2u⟩ decaying before the completion of STIRAP is also due to converted VPs which are being detected.Besides being much simpler to implement, the continuous measurement scheme is faster (notice the time scales in Fig. 2a and Fig. 1b) lowering the relative contribution of the stray thermal current.This is apparent in Fig. 2b where we plot for each cycle after thermalization the instantaneous photon current (blue dotted) and the averaged total (blue) and thermal (orange) currents which are proportional to the power respectively emitted into the transmission line.Their difference is the signal due to the converted VPs, which is related to the grey dotted line in Fig. 2a.Summing up, for the continuous measurement scheme the tradeoff between efficient measurement and decoherence is positive, yielding a sufficiently large extra output power from converted VPs.
We briefly comment on the role of atomic decay.Fig. 2a shows that in the integrated protocol, some population appears in |Φ⟩ before the completion of STIRAP since the pump pulse repumps to |Φ⟩ population which just decayed in |0u⟩ because of the always-on coupling to the meter.In the absence of atomic decay, this population would remain trapped in |Φ⟩ after the completion of STIRAP.If γ ̸ = 0 this population relaxes non-radiatively to |0u⟩ resetting efficiently the system for the next cycle.
In Ref. [24] it has been estimated that a signal corresponding to the case study of Table 1a could be amplified by standard HEMT circuitry and discriminated from thermal noise, this task requiring hundreds of repetitions of the conversion/measurement cycle, which is a reasonable figure.In this work, we have analysed in detail the dynamics of the continuous measurement, showing that STIRAP is resilient to measurement backaction.It is likely that combining optimal control theory and advanced computational methods of data analysis [31,44] yields even better figures.For instance, Fig. 2b suggests that shortening the after-STIRAP part of the protocol yields a larger average output power for the converted VPs, while the thermal floor is unchanged.However, the steady-state population of |0u⟩ may be smaller and that population of the intermediate state may trigger leakage from the four-level subspace, requiring an analysis which takes into account the multilevel nature of the setup proposed in Ref. [24] and the subtleties of the physics of an open system in the USC regime.

Fig. 1 .
Fig. 1.(a) Spectrum of the multilevel AA extended Rabi model with the additional uncoupled state |u⟩, Eq.(3), as a function of g and the scheme of the Lambda configuration used to drive the system, Eq.(4).(b) Top panel: Gaussian pulses Ω s/p (t) in the counterintuitive sequence for the coherent amplification protocol and the sigmoid function κ(t−tsm) mimicking a switchable meter.Bottom panel: population histories; for t < tsm population transfer |0u⟩ → |2u⟩ by STIRAP is completed; for t > tsm the mode decays, emitting two photons and resetting the system to the initial state.Parameters are given in Table1.The protocol starts at Ω0ti = −75 and ends at Ω0t f = 450, the meter is switched on at Ω0tsm = 90.We used the shorthand notation for the populations Pn := ⟨n, u|ρ(t)|n, u⟩ for n = 0, 1, 2 and PΦ := ⟨Φ|ρ(t)|Φ⟩, where ρ(t) is the density matrix of the system.

Fig. 2 .
Fig. 2. (a) A cycle of integrated conversion/photodetection protocol at finite temperatures.The meter is always on, i.e. κ ̸ = 0 (not represented) is independent on time, allowing for an overall time tM shorter than in Fig.1(here Ω0t f = 350).Parameters are given in Table1and for the case study considered tM = 425/Ω0 = 1.36 µs.The black dotted line is the total number of photons decaying in the transmission line: it becomes larger than two because of the thermal contribution which is linear in t.By subtracting it, we obtain the number of converted VPs (grey dotted line).The four horizontal dashed lines are the thermal populations.(b) The instantaneous emitted current (dotted line) and its average per cycle (blue solid) compared to the thermal current (orange).