Probing ultrastrong light–matter coupling in open quantum systems

Dynamically probing systems of ultrastrongly coupled light and matter by advanced coherent control have been recently proposed as a unique tool for detecting peculiar quantum features of this regime. Coherence allows in principle on-demand conversion of virtual photons dressing the entangled eigenstates of the system to real ones, with unitary efficiency and remarkable robustness. Here, we study this effect in the presence of decoherence, showing that also in far from ideal regimes is it possible to probe such peculiar features.


Introduction
Light-matter interaction is a fundamental building block of Nature leading to countless applications, whose scope broadened in recent years with the advent of quantum technologies [1]. Circuit-QED solid-state systems are one of the forefront platforms for quantum hardware [2] where besides applications, fundamental physics from measurement theory [3] to quantum thermodynamics [4] and quantum communication [5] can be studied. The basic physics is described by the Rabi model [6], where a two-level atom interacts with a single quantized oscillating mode. In the strong coupling regime (SC), the system exhibits exquisitely quantum dynamics, conserving the number of excitations to an enormous degree of accuracy; this approximate symmetry breaks down in the ultrastrong coupling (USC) regime [7]. Systems in this regime have been demonstrated in the last few years in circuit-QED architectures of solid-state artificial atoms (AA), THz metamaterials, intersubband polaritons and other physical systems [8,9]. USC systems have been subject of extensive theoretical investigation predicting novel and outstanding non-perturbative effects [8][9][10][11]. Amongst them, it has been pointed out that non-conservation of the excitation number implies that the ground state |Φ 0 of the Rabi model in the USC regime contains virtual photons. Their detectability has been the subject of intense theoretical investigation [12][13][14][15][16][17][18] but experimental demonstrations are still lacking. To understand and fill this gap, the very recent work [17] addressed key issues on hardware design and control, showing how to demonstrate efficient and unambiguous photon pair convera e-mail: giuseppe.falci@unict.it (corresponding author) sion, amplified using an advanced coherent control protocol analogous to stimulated Raman adiabatic passage (STIRAP) [19], which has been proposed and demonstrated in AAs [20][21][22][23][24]. In this work, we focus on the evaluating the impact of decoherence processes in AAs [25] on the dynamics of the ideally quantum-controlled systems.

Coherent amplification of photon pair production
The simple version of the coherent amplification protocol technique was proposed in Ref. [15]. A three-level atom (basis {|u , |g , |e }) is ultrastrongly coupled to a quantized mode with oscillation frequency ω c , resonant with the e.g transition energy ε. The Hamiltonian is ( = 1) where ε is the atomic u − g splitting, a and a † are the photon annihilation and creation operators, the eigenstates of a † a forming the basis of the Fock states {|n }.
Here λ is the coupling between the mode and the e−g transition, whereas |u is uncoupled (as for − ω c λ). The so-called diamagnetic term [8,9] is implicitly accounted for by renormalized parameters, but it plays a negligible role in the regime of couplings considered in this work. Eigenstates of H are partitioned in two sets: they are strongly entangled atom-mode states, dressed by virtual photon pairs; (2) factorized states |nu , with energy ε + n ω c , which will be used as ancillas. Extra terms may trigger transitions between these subspaces. In particular atomic u − g couplings may determine the conversion of virtual photons of the Rabi subspace in real photons of the ancillary subspace, where they can be detected. At resonance, ω c = virtual photons in the (false) Rabi vacuum |Φ 0 are directly converted to real photon pairs, the main channel involving the target state {|u 2 }, with corresponding amplitude c 0 2 (λ/ω c ). Since this latter is large enough only in the USC regime, the detection of photon pairs, while the AA is found in |u , provides a unique dynamical signature of USC. The above photon pair production channel could be amplified, and thus made detectable, using coherent control. The atom is driven by two-tone control field in Λ configuration described by the Hamiltonian (see Fig. 1a) We take with width T , delay τ > 0 and Stokes/pump peak field amplitude W s/p [15,19]. This protocol yields non-vanishing coherent population transfer |0u → |2u only if Φ 0 |H c |nu = 0, i.e., when |Φ 0 contains pairs of virtual photons. In the regime 0.2 λ/ω c 0.5, which is very interesting for experiments, the dynamics can be conveniently visualized as a Λ-STIRAP population transfer in the "STI-RAP subspace" {|u 0 , |u 2 , |Φ 0 } and control may be tailored to yield ∼ 100% probability to find eventually two photons [26], the AA always remaining in the state |u . We remark that the implementation of dynamical detection schemes in real systems poses severe constraints and challenges. This has been exhaustively discussed in the works [17,26], where detection of virtual photon pairs was studied in superconducting systems using the flux-qubit design and advanced control [27]. Here, we focus on the single issue of evaluating the impact of decoherence on the dynamics of the ideally quantum-controlled systems.

Partially coherent photon pair conversion
We stress that for the detection of USC, it would be sufficient to monitor an increase of the population of the Fock states |n ≥ 2 during part of the protocol. Some transient population of the intermediate state is also tolerable, softening the adiabaticity requirement. Decoherence times T φ ∼ T can also be tolerated [28] since at worst efficiency of the USC-selective channel would be larger than 30%. This opens perspectives also for semiconducting structures, where USC-selective Λ-STIRAP could be observed with some progress in techniques for detecting excess THz photons. Therefore, we address the problem of observability of photon pair conversion when coherence is partially lost.

Modeling decoherence
The open quantum system dynamics for Λ STIRAP is studied considering Markovian decay processes of both the cavity and the AA, accounted for by the following Lindblad equatioṅ The dissipator L cav [ρ] = κ 2 (2aρa † − a † aρ − ρa † a) describes losses of the cavity mode with rate κ. Spontaneous decay of the three-level atom is modeled by where σ ug = |u g|, σ gg = |g g|, σ ge = |g e| and σ ee = |e e|. This is a minimal model aiming to understand the impact of leading processes, assuming suitable external conditions. Temperature much smaller than ω c is assumed as well as and not so strong driving, in order to neglect atomic absorption and stimulated emission. We discard e − u transitions, having in mind a three-level AA implemented by a devices biased at a symmetry point, which enforces parity selection rules. Symmetries also cancel at lowest order pure dephasing due to low-frequency noise both of the AA [25,[29][30][31] and the cavity. We moreover assume that the environments has power spectrum slowly varying in a frequency interval δω ∼ λ around ω = ω c , so that rates are weakly dependent on λ.

Results for the Λ configuration
The dynamics for Λ scheme is carried out by numerically solving Eq. (3). An example of population histories for κ = 0 and γ = 0 is shown in Fig. 1b, colored lines, which exhibit the increase of P |2u due to USC while P |0u and P |Φ0 remain small. Notice that population eventually leaks outside the "STIRAP subspace" due to photon decay |2u → |1u in the cavity whose effect is studied in Fig. 1c. The same figure shows that on the contrary, atomic decay γ is almost irrelevant (gray dotted and dashed lines) as long as the population in P |Φ0 remains small. Actually e − g processes are even less relevant since these states are practically not populated and neglecting it would not make any appreciable difference.
In the same figure, we show that population transfer may also occur due to a stray corotant coupling of the cavity to the u − g transition, H ug = λ |u g| a † + |g u| a , while λ = 0. As shown in Ref. [17], in the presence of the stray coupling, detecting two photons is not anymore a smoking-gun for USC in available AA-based quantum device, spoiling the experimental significance of all the proposals based on transitions in the Λ scheme [13][14][15][16]. On the other hand, Fig. 1c shows that if stray processes could be somehow suppressed, there is room

The Vee configuration
To circumvent the problem of stray couplings, an alternative control scheme operating in the Vee configuration has been introduced in Ref. [17]. This work shows that remarkably stray couplings in this case play no role. Moreover control fields may be coupled more easily to the relevant transitions allowing faster operations. On the other hand, noise acts in a different way; thus, also in this case, the regime of partial coherence has to be investigated in detail. In the Vee scheme the lower doublet of the artificial atom is coupled to the mode, while the ancillary level has higher energy. The Hamiltonian is the same Eq.
The results in Fig. 2 show that in this case, the atomic transition rate γ has in general an impact even larger than the cavity losses κ. This is evident both from the population histories, Fig. 2a, and the maximal efficiency Fig. 2b (dashed and dotted lines vs the red line). On the other hand, the trade-off of this drawback with the possibility of faster processes (notice the time scale in Fig. 2a) is favorable, showing that Vee-STIRAP is also robust against decoherence. Indeed in Fig. 2b, where the efficiency of the protocol is reported at the maximum value obtained during the evolution, figures of ∼ 50% up to a decay rate κ/ω c 2×10 −3 are displayed. Also in this case, decay of level e → g is negligible since |e does not get populated during the dynamics owing to the high efficiency of Vee-STIRAP. It is finally shown that the effect of stray corotating couplings is irrelevant also in the presence of new channels opened by dissipation, thus making Vee-STIRAP a unique tool for the dynamical detection of virtual photons in the USC regime.

Conclusions
We studied a theoretical proposal of dynamical detection of virtual photons in the highly entangled eigenstates of an AA coupled with a quantized oscillating mode in the USC regime. Our proposal leverages on the coherent amplification of photon pair production obtained in the ideal case of a closed system by operating a control scheme similar to STIRAP. We show that detection of photon pairs is also possible with partial coherence only. In particular by driving the system in Λ configuration, the most relevant decoherence mechanism is the photon losses of the cavity. The fundamental difficulty for this protocol is that in the so far available hardware, the unavoidable presence of stray couplings between the cavity and unwanted transitions in the atom may yield pair production also in the absence of USC. Our results show that if this problem is mitigated by new design or control solutions, then decoher- On the other hand, the problem of stray couplings can be circumvented by resorting to STIRAP in the Vee configuration. We have shown that in this case, the stray channels are suppressed also in the presence of decay processes. Moreover, while the contribution of atomic spontaneous decay has also a direct impact in lowering the efficiency, operations are faster the tradeoff between the two effects being positive. This opens realistic perspectives to the detection of virtual USC photons in state of the art solid-state architectures.
Acknowledgements This work was supported by the QuantERA grant SiUCs by CNR GA nr. 731473 Quan-tERA, and by University of Catania, Piano per la Ricerca 2016-18-linea di intervento "Chance", Piano di Incentivi per la Ricerca di Ateneo 2020/2022, proposal Q-ICT.
Funding Information Open access funding provided by Università degli Studi di Catania within the CRUI-CARE Agreement.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.