F Assemat

D Grosso

A Signoles

A Facon

I Dotsenko

S Haroche

J M M Raimond

Michel Brune

S Gleyzes
email: gleyzes@lkb.ens.fr

Quantum Rabi oscillations in coherent and in mesoscopic "cat" field states

Keywords:

come

The Rabi oscillations of a two-level atom in a resonant, single-mode coherent field state is one of the simplest phenomena in quantum optics. Nevertheless, it exhibits surprisingly complex features at the mesoscopic scale (few tens of photons) [1][2][3][4]. The oscillations, at an angular frequency Ω 0 √ n, collapse and revive (n is the average photon number in the coherent state; Ω 0 is the vacuum Rabi frequency measuring the atom-field coupling). The collapse, occurring on a time scale T c = 2 √ 2/Ω 0 , results from the quantum field amplitude uncertainty and from the corresponding dephasing of the Rabi oscillations. The (first) revival, around T r = 4π √ n/Ω 0 , results from the rephasing of oscillations associated to different photon numbers [5]. This revival provides a landmark illustration of field amplitude quantization [6]. Between collapse and revival, the field evolves into an entangled atom-field state, involving two coherent states with different phases [7,11,[35][36][37]. It is called a "cat state" in memory of Schrödinger's metaphor. Close to t = T r /2, the atomic state factors out of a field "cat", superposition of coherent states with opposite phases [5]. These phenomena can be observed in systems implementing the Jaynes and Cummings model, a spin-1/2 coupled to a one-dimensional harmonic oscillator [START_REF] Jaynes | Proc. IEEE[END_REF]. Ions in traps [START_REF] Meekhof | [END_REF]14], cavity quantum electrodynamics (CQED) [3,6] and circuit quantum electrodynamics (cQED) [16,39] are thus ideal platforms for this observation.

Nevertheless, revival observations have so far been limited to small photon numbers since experiments face formidable challenges. For microwave CQED with superconducting cavities crossed by fast Rydberg atoms, the interaction time is limited to a few vacuum Rabi periods, 2π/Ω 0 . Revivals have been observed only for n 1 [6,17]. Early revivals induced by a time-reversal of the collapse can be observed for larger n values (about 10), but the maximum separation of the cat compo-nents is small [18,19]. Ion traps have similar limitations [START_REF] Meekhof | [END_REF]20]. In cQED, the limited coherence time of tunable superconducting qubits makes it difficult to observe long-term dynamics in the resonant regime [21]. Large cat-state preparation so far relies instead on the dispersive, non-resonant interaction [22,39], in which the atom is simply a transparent dielectric material with a state-dependent index of refraction [3].

In this Letter, we push the quantum revival phenomenon at a much larger scale. Using slow circular Rydberg atoms crossing a high-Q superconducting cavity, we achieve atom-field interaction times up to twenty vacuum Rabi oscillations periods. We observe the complete first revival for n = 13.2. Resetting the atom close to t = T r /2, we leave in the cavity a cat state. We observe the quantum revivals in this initial cat state. We use them to measure the cat state Wigner function and to investigate its fast decoherence.

The experiment is sketched in Fig. 1(a). Additional details are given in [5]. A 2D-MOT and an additional longitudinal velocity selection produce a Rubidium atomic beam propagating upwards along the Ox axis at an average velocity v = 8.4 m/s towards the cavity C.

Inside the cavity, atomic samples are selectively prepared in |e , the circular state 51C with a principal quantum number 51, by laser, radio-frequency and microwave excitation [41]. Each sample contains 0.08 atoms on the average. Events with two atoms simultaneously present in C have thus negligible influence. The cavity is tuned close to resonance with the |e → |g transition at 51.1 GHz (|g is the 50C circular state). An electric field along the cavity axis, Oy, produced by a voltage applied across the mirrors, stabilizes the circular states and makes it possible to Stark-tune the atomic transition frequency in or out of resonance with C. After the interaction with the cavity mode, the atoms drift towards a field-ionization counter D, allowing us to measure the populations in |e or |g (see [5] for details).

The cavity C, cooled down to 1.5 K by a wet 4 He cryostat, sustains a linearly polarized Gaussian standingwave mode with a waist w = 6 mm [3]. Its damping time is T Cav = 8.1 ± 0.3 ms. Its temperature corresponds to an average thermal photon number n th = 0.38. In order to reduce the residual photon number in C, we send absorbing atoms prepared in |g , starting 2.7 ms before the sequence. A microwave source S, coupled to C, performs tunable coherent displacements of its mode. The injected amplitude β is a linear function of the injection time t β . For injections lasting more than ≈ 100 ns, the calibration is |β| = 0.257(t β (µs) -0.05) [5].

We first record the vacuum Rabi oscillations with no injection (β = 0) and an atom initially in |e . The sequence timing is sketched in Fig. 1(b). The initial state |e is prepared in a large electric field, resulting in a δ/2π = 1.4 MHz atom-cavity detuning, for which the interaction is negligible. We abruptly set δ = 0 at t = 0. After a variable interaction time t i , we set back δ to a large value (2π × 4.04 MHz), halting the evolution. We measure the probability, P g , for finding the atom in |g . Figure 2(a) shows the experimental P g (t e ) (dots with statistical error bars, joined by a thin black line) as a function of the effective interaction time t e taking into account the motion of the atom through the Gaussian mode geometry [5]. We observe nearly 20 vacuum Rabi periods, a considerable improvement over previous experiments [6]. The red line in Fig. 2(a) results from a numerical simulation of the experiment taking into ac-count the detection imperfections, resulting in a reduced oscillation contrast, as well as the residual initial thermal field and cavity relaxation towards thermal equilibrium [5]. The Rabi frequency, Ω 0 /2π = 49.8 kHz, the probability p 1 = 0.09 for having initially one photon in C and the precise position of the exciting laser beams inside the cavity mode, defining the relation between the real and effective interaction times are extracted from a fitting procedure described in [5].

We record now the Rabi oscillation in a coherent state |β with |β| 2 = n ≈ 13.2 ± 0.1 [timing in Fig. 1(b), with t β = 14 µs]. Note that β can be assumed to be real without loss of generality. The signal is plotted in Fig. 2(b) (points with statistical error bars joined by a thin black line). It clearly exhibits the first revival around T r = 146 µs. The signal is in excellent agreement with a numerical simulation of the experiment (red solid line) taking into account the same imperfections as those in Fig. 2(a).

The Fourier transform of the revival signal [dots in Fig. 2(c) with a fit to a sum of gaussians -solid red line] exhibits discrete peaks, very close to the expected Rabi frequencies in an n-photon field, Ω 0 √ n + 1/2π (vertical blue lines). They are well-resolved up to about 18 photons. This signal provides a textbook evidence of field quantization. The weight of the peaks directly measures the photon number distribution, p(n) [5]. It is in excellent agreement with the expected Poisson distribution for a mean initial photon number 13.2 (blue solid line).

Close to half the revival time, t e ≈ T r /2 ≈ 73 µs, the atomic state is expected to factor out [5]. The cavity is then in a superposition of two components with nearly opposite phases [35,36]. The photon number parity P = n (-1) n p(n) oscillates with t e [24] and the field Wigner function, W (α), exhibits fringes near the origin of phase space revealing the cat's coherence [38].

In order to probe the cavity state, we make use of its resonant interaction with a 'probe' atom, initially in |e . For a well-chosen [5] effective interaction time t e close to T r /2, the probability P g for finding the atom in |g is P g (t e ) ≈ (1 -P) /2, and P g measures the parity P = πW (0)/2. By displacing the cat state by α before the resonant interaction, we get access to W (-α).

The timing of the experiment is depicted in Fig. 1(c). We set the first resonant effective interaction time to t i = 60 µs. After the atom has been tuned out of resonance, we perform a displacement by a real amplitude α (injection time t α ). Simultaneously, we get rid of the |g part of the atomic state by a resonant radio-frequency pulse transferring |g to low-m states, finally undetected. Thus, any detected atom was in |e at the beginning of a new resonant interaction, starting after a delay t d = 6 µs and lasting for a variable time t i corresponding to the effective time t e .

Figure 3(a) shows P g (t e , α) (dots with statistical error bars) as a function of α for t e = 68.5 µs T r /2. We have chosen this t e value to optimize the contrast of the oscillations, so that P g (t e , α) ≈ [1 -πW (-α)/2]/2. The solid red line results from the simulation of the experiment and is in excellent agreement with the data. From the period of the oscillations in Fig. 3(a), we deduce the size of the cat, D 2 , measured by the square of the distance D of the superposed components in phase space [5]. We get D 2 = 45.1 ± 0.4 photons, in excellent agreement with the value provided by the numerical simulation. This size compares with that of the largest cats ever prepared by dispersive interaction [39]. ment [arrow on Fig. 3(a)] brings one of the extrema of W onto the origin of phase space, resulting in a displaced cat with a maximum parity. We clearly observe here a dual revival, one at the standard revival time T r and one earlier, at T r /2, revealing the displaced cat parity. The agreement of the experimental data with the simulation (red solid line) is excellent. The Fourier transform of the Rabi signal, shown in Fig. 3(c), conspicuously displays the displaced cat parity. The photon number distribution inferred from this spectrum provides a displaced cat parity P = -0.48. The simulation provides P = -0.41 in excellent agreement with the observed value. The parity for an ideal experiment taking into account cavity relaxation, is expected to be -0.49, very close to the observed one. We thus generate large cats with a good fidelity.

The decoherence time scale of the cat in the finite temperature environment of the cavity is expected to be 200 µs [5,40]. In order to investigate this fast decoherence, we monitor the amplitude of the revival at T r /2 as a function of the delay time t d (note that the displacement by α does not change the dynamics of decoherence).

Fig. 4(a-d) shows four Rabi oscillations signals for t e close to T r /2 for α = -0.60 and for four values of t d (dots with statistical error bars and simulation of the experiment as a solid red line). The agreement between simulation and data is excellent, confirming the decoherence time scale. Note that the moderate T c value, much shorter than 130 ms reached in previous experiments [38], is appropriate for an observation of the cat's decoherence in the available interaction time range.

In conclusion, we have observed the resonant Rabi oscillation in a coherent field in an unprecedented range of photon numbers. We have shown the generation of large cat states through the first observation of an early revival for the Rabi oscillation in a cat with a well-defined parity. We have monitored the fast decoherence of these non-classical states.

The resonant interaction generates cat states efficiently and significantly faster than the dispersive method [5]. Using more than one atom, it can lead to the preparation of more complex state superpositions with multiple components [27,28]. It is thus a promising method for fundamental decoherence studies, but also, in the cQED context, for the use of cats in quantum error correction protocols [29,30].

This experiment also opens the way to a new realm for atomic physics CQED, with extremely long interaction times and extremely low damping cavities. Particularly promising is the use of laser-trapped circular Rydberg atoms [31]. They could be combined with a 3Dmicrowave structure sustaining a high-quality resonant mode, but admitting no other modes resonant with the atomic spontaneous emission decay channels. One might then combine the best of the CQED and cQED worlds, with reproducible atoms, well-controlled and understood decoherence channels, nearly infinite interaction times and the slow pace of CQED experiments instrumental for real-time control [32].

We acknowledge funding by the EU under the ERC projet 'DECLIC' (Project ID: 246932) and the RIA project 'RYSQ' (Project ID: 640378).

SUPPLEMENTARY MATERIAL: QUANTUM RABI OSCILLATIONS IN COHERENT AND IN MESOSCOPIC SCHR ÖDINGER "CAT" FIELD STATES

We first give a simple theoretical model of the experiment. We then provide additional details about the experimental apparatus and the data acquisition. We then discuss the relations between the real resonant interaction time and the effective one. Finally, we discuss the data processing and the numerical simulation of the experiment.

Collapse, revival and cat generation: theoretical remarks

We examine here the collapse and revival for an initial coherent or cat-like state. We consider an atom initially in |e . The cavity C contains a quantum state Ψ = n c n |n (|n is the Fock state with n photons), corresponding to the photon number distribution

p(n) = |c n | 2 . For a coherent state, p(n) = p c (n), with p c (n) = e -n n n /n!,
where n is the average photon number. A simple calculation [START_REF] Haroche | Exploring the Quantum: Atoms, Cavities and Photons[END_REF] provides the probability P g (t) for detecting the atom in |g after an effective interaction time t e

P g (t e ) = 1 2 1 - n p(n) cos Ω 0 √ n + 1 t e . (S1)
The Rabi signal is a sum of the Rabi oscillations in the Fock states |n , at a frequency Ω 0 √ n + 1, weighted by the photon number distribution. The complex interplay of the non-rational frequencies appearing in this sum leads to the collapse and revival phenomenon.

Let us first consider the case of a coherent state, |β , with a √ n dispersion of the photon number distribution, p c (n), around n. For short effective interaction times, the dispersion of the oscillation phases in Eq. (S1) is negligible and we find, as expected, a Rabi oscillation at the average angular frequency Ω r = Ω 0 √ n with a close-tounity contrast. The collapse is due to the spread of the Rabi frequencies in Eq. (S1). In order to estimate the collapse time, we can compute the time, T c 2π/Ω 0 , after which the Rabi oscillations for two photon numbers separated by the photon number variance, n = n ± √ n/2, get in phase opposition. A more detailed calculation shows that the envelope of the oscillations collapse is a gaussian with a 1/e time-width T c = 2 √ 2/Ω 0 [START_REF] Haroche | Exploring the Quantum: Atoms, Cavities and Photons[END_REF]. The revivals correspond to the rephasing of the Rabi oscillations associated to different photon numbers. For n n, it is easy to see that the oscillations for n and n+1 photons come back in phase at a time T r = 4π √ n/Ω 0 . Since this time is independent of the chosen n value, all oscillations come back in phase at that time, resulting in the maximum of the revival.

An early quantum revival can be expected near t e = T r /2 if the initial cavity state |Ψ has a well-defined photon number parity P = n (-1) n p(n) = exp(iπN ) (N is the photon number operator). This is the case for an ideal 'phase ' or 'parity cat' (|γ ± | -γ )/ √ 2 (we assume γ 1 for the normalization factor), for which P = ±1. Choosing t e = 2π

√ n 0 /Ω 0 + τ , where n 0 is an even integer close to n, we get from Eq. (S1), expanding √ n + 1 at the first order in 1/ √ n 0

P g (t e ) = 1 2 [1 -P cos Φ] , (S2)
where

Φ ≈ Ω r τ. (S3)
We thus get Rabi oscillations at the frequency Ω r , with an amplitude proportional to the parity P of the initial state. For an initial coherent state, P = exp(-2n) and there is no sizable Rabi oscillation near T r /2. For a parity cat state, there are instead well-contrasted oscillations at half the revival time. This result can be retrieved by a simple argument. For a parity-cat, p(n) cancels for odd or even photon numbers. The least spacing between Fock states contributing to the sum in Eq. (S1) is 2 and the rephasing time is thus reduced by a factor of 2. Choosing a time τ such that cos Φ = 1, the measurement of P g directly provides P, proportional to the field's Wigner function at the origin, W (0) [START_REF] Lutterbach | [END_REF]. Adding a translation of the initial state by an amplitude α [represented by the displacement operator D(α)] before starting the Rabi oscillation thus directly determines the Wigner function, W (-α) at one point in phase space:

W (-α) = 2 π D(-α)PD(α) . (S4)
This simple analysis in terms of a beat between Rabi oscillations leaves aside the evolution of the field state. A much more detailed insight is provided in [START_REF] Haroche | Exploring the Quantum: Atoms, Cavities and Photons[END_REF][35][36][37]. Since the initial field injection in C provides the phase reference, we can assume that β is real and positive without loss of generality. Expanding the full atom-cavity state to the first non-trivial order in (n -n)/n and Ω = Ω 0 /4 √ n. The state |Ψ(t e ) is generally an entangled atom-field state, a cat state in which two atomic states |ψ ± a (t e ) , slowly rotating at frequency ±Ω in the equatorial plane of the Bloch sphere, are correlated to two coherent field components, also rotating at ±Ω in the Fresnel plane from the initial |β state. The two coherent components overlap for t e close to 0 and for t e close to all multiples of T r . The field state then factors out and the highcontrast Rabi oscillations at Ω r result from a quantum beat between the two atomic states. In this simple picture, the revivals are periodic. In fact, the higher order terms in the developments in powers of (n -n)/n result in a progressive phase spread (squeezing) of the coherent components. The long-term revivals are thus longer and less contrasted than the first. Already at the first revival, the phase spread is of the same order as the initial phase uncertainty in |β [START_REF] Haroche | Exploring the Quantum: Atoms, Cavities and Photons[END_REF]. It is thus about twice longer than the collapse and has a maximal ideal contrast close to 50%.

Another interesting situation occurs at half the revival time, T r /2, for which Ωt e = π/2. The two atomic states are thus both equal, within a global phase, to (|e -|g ) √ 2. The atom disentangles from the field. The probability P g for detecting it in g is 1/2 and the field is then left in

|ψ cat = e iπn |iβ -| -iβ , (S8)
a cat state with the parity

P cat = -cos(πn) . ( S9
)
Let us compare this resonant generation of a cat state with the dispersive one, commonly used in Cavity Quantum Electrodynamics [38] and particularly in circuit QED [39]. In the dispersive regime, the atom in state |e or |g causes a shift of the cavity frequency by ±Ω 2 0 /4δ (δ is the atom-field detuning). Preparing the atom in a state in the equatorial plane of the Bloch sphere and detecting it after a π/2 resonant pulse mixing again |e and |g after the interaction with the cavity (and thus erasing any which-path information about the atomic state in the cavity), we prepare a cat state, superposition of the coherent components |β exp(±iϕ) , with ϕ = Ω 2 0 t d /4δ (t d is the effective interaction time for a dispersive interaction, close to the effective interaction time for a resonant interaction, t e ). A parity cat is reached for t d = 2πδ/Ω 2 0 . When δ > Ω 0 √ n, a condition essential to enforce the dispersive regime, the dispersive interaction is slower by a factor δ/Ω 0 √ n than the resonant one to generate a phase cat. A more detailed discussion, taking into account intermediate regimes, shows that the cat preparation time is minimal, all other things equal, when δ = 0 [START_REF] Haroche | Exploring the Quantum: Atoms, Cavities and Photons[END_REF].

The cat states are very sensitive to decoherence. The Lindblad equation formalism leads to an exact expression for the density matrix of a relaxing cat state. In an FIG. S1: Detailed scheme of the experiment with all beams and an indication of distances environment at zero temperature, the coherence between the cats components is found to decrease, for t T Cav , as exp(-2nt/T Cav ) [START_REF] Haroche | Exploring the Quantum: Atoms, Cavities and Photons[END_REF]. The decoherence time scale is thus T D = T Cav /(2n), 2n times shorter than the cavity energy damping time. This fast relaxation is a key feature of the decoherence of mesoscopic state superpositions. It can be qualitatively understood in terms of parity. Starting from an odd cat, a photon loss leads to an even cat. After having lost on the average about one photon, the field is in a mixture of odd and even cats, which is also a statistical mixture of the two coherent components. The coherence is thus lost after about one photon has disappeared, after a time of the order of T Cav /n. At a non-zero temperature, the decoherence time is

T D = T Cav /2[n(1 + 2n th ) + n th ],
where n th is the average thermal photon number [40].

EXPERIMENTAL DETAILS

A detailed scheme of the experiment is shown in Fig. S1. A slow vertical 85 Rubidium atomic beam effuses from a 2D-MOT (SYRTE Laboratory). Its transverse velocity spread is reduced by a transverse optical molasses stage, situated right above the 2D-MOT arrangement. The longitudinal beam velocity distribution, shown in Fig. S2 as a red line, is inferred from a circular-atom timeof-flight measurement from the laser-excitation region, close to the center of the cavity, to the field-ionization detector (note that the velocity-dependent efficiency of the circular state preparation and the finite lifetime of the circular levels slightly distort the measured distribution). The distribution is peaked at about 8.50 m/s, with a width of 2.4 m/s. This dispersion is too large for a good control of the experimental sequence. We therefore perform an additional passive longitudinal velocity-selection stage between the molasses and the cryostat. Two resonant laser beams, separated by 14 cm are used to push the atoms out of the beam axis, interrupting the atomic flux. They are periodically interrupted for 5 ms each 15 ms. The atoms reaching the cavity must have passed the two selection beams while they were off. The resulting velocity distribution is shown as a black line in Fig. S2. It is centered at 8.36 m/s, with a width 1.2 m/s.

The circular atoms are prepared by laser, rf and mw excitation inside the cavity C [41]. We use three lasers, perpendicular to the atomic beam at 780, 776 and 1258 nm resonant on the 5S 1/2 → 5P 3/2 , 5P 3/2 → 5D 5/2 and 5D 5/2 → 52F transitions. The pulsed laser excitation (0.1 µs pulse for the third laser) is performed in a weak electric field F (0.14 V/cm) produced by a voltage applied across the cavity mirrors. The Stark structure being resolved, we only excite levels with a magnetic quantum number m = ±2. We then raise F to 2.4 V/cm. The atomic state adiabatically evolves into the lowest m = ±2 levels in the hydrogenic Stark manifold. We apply a σ +polarized rf field at 230 MHz, slightly below resonance with the Stark frequency in the manifold, while slowly decreasing F . The m = 2 level is efficiently transferred to the circular state 52C by an adiabatic rapid passage process. It leaves the m = -2 state unaffected (transitions to the m = -1 Stark levels are out of resonance due to the finite quantum defects of the low-states). These m = -2 states do not interact with C and are not finally detected. A final 1 µs microwave pulse prepares the initial level 51C (|e ), leaving in the 52 manifold spurious high-m states and about 20 % of the 52C population.

For initial field erasing, we send resonant absorbing atoms in C before the start of the sequence. We use a two-photon mw pulse resonant on the 52C → 50C to prepare directly about 10 atoms in |g . The last fielderasing atoms are prepared 2.7 ms before the the start of the actual sequence so that they do not contaminate the field-ionization detection of the useful samples.

The cavity C is made up of Niobium-coated diamondmachined copper mirrors. We measure its resonant frequency and its lifetime (8.1 ms) by injecting a field into it via its diffraction loss channels. Probe atoms are used to measure this field and its decay versus time [42].

The field injections in the experimental sequence are also performed by S, used with a constant power. The injected amplitude modulus |β|, measured at time t = 0 in the sequence (the time at which the atom is first set at resonance with the mode) is a linear function of the injection time t β . We obtain an approximate calibration of this linear law by recording the first periods of the resonant Rabi oscillation as a function of the injection time and by fitting them to a simulation. A better determination is provided by the Rabi revival signal itself, whose Fourier transform directly measures the photon number distribution. A slight refinement within the error bars of the previous method is provided by adjusting finely the scale of the amplitude axis used for the experimental data in Fig. 3(a) (main text). We finally use the linear law given in the main text. The small offset on the injection time is due to the finite commutation time of the X-band source.

The relative phase, ϕ, of the two injections used for the parity measurements is set to π by adjusting the detuning ∆ of the source with respect to the cavity (in the cavity-field rotating frame, the source field rotates at ∆). The 103.8 µs time separation between the two injections makes ϕ more sensitive to the precise value of the cavity frequency than the Rabi oscillation itself (a cavity shift of only 2.4 kHz, much lower than Ω 0 /2π, shifts ϕ by π). We thus check the cavity resonance frequency every day. We have devoted considerable efforts to the reduction of cavity mechanical vibrations due, in particular, to cryogenic fluids oscillations.

The field-ionization detector D is active between 7.5 and 9 ms after the initial laser pulse, so that we detect the lower-velocity part of the beam distribution, providing an additional velocity-selection stage. During their time of flight towards D, the spontaneous decay of the circular states cannot be neglected (lifetime of the order of 30 ms). A part of the atoms initially in |e is thus transferred into |g and affects the measured population. Similarly, a part of the atoms left in the 52 manifold may decay in the 51 one and be spuriously detected in |e . In order to avoid cross-talks between levels populations, we first transfer, 450 µs after the laser excitation, the atoms in |g into the 48C circular state with a two-photon microwave pulse. We then use, 600 µs after the laser pulse, a resonant rf pulse, transferring the atoms of the 52C manifold into low-m states, which have a ionization threshold very different from that of the circular state (and which, moreover, rapidly decay by optical transitions). A similar procedure is used to get rid of the atoms left in |g after the first Rabi oscillation period in the cat parity measurements.

DATA ANALYZIS

The essential parameters of the experiment are extracted from a direct fit of the raw vacuum Rabi oscillation signal, P g (t i ), recorded as a function of the interaction time t i . This signal is shown in Fig. S3 (note that the signal in Fig. 2(a) of the main text is plotted as a function of the effective interaction time t e ). The fitting function is the expected Rabi oscillation in a mixture of zero and one-photon states describing the initial residual thermal field in the cavity. The ideal P g values, P i g , are then given by

P i g = (1 -p 1 ) cos 2 (Ω 0 t e /2) + p 1 cos 2 ( √ 2Ω 0 t e /2) (S10)
where p 1 is the initial probability for having one photon in C. The time variable in this expression is the effective interaction time t e . It takes into account the motion of the atom along axis Ox through the gaussian transverse electromagnetic TEM mode geometry, characterized by its waist w = 6 mm. It reads, as a function of t i :

t e (t i ) = w v √ π 2 erf (t i -x 0 v) v w -erf -x 0 v v w , (
S11) In this expression, x 0 is the initial (t = 0) position of the atom on the Ox axis, determined by the position of the exciting laser beams. It is not precisely known and must be extracted from the fitting procedure. The atomic velocity v = 8.1 ± 0.1 m/s is well-known from the velocity distribution shown in Fig. S2.

We then have to take into account in the fit the detection errors. They include the 50C -48C probe efficiency, the efficiency of the rf pulse used to get rid of the 52C, the contamination of 51C into 50C due to spontaneous emission between cavity interaction and the pulse transferring 50C into 48C, the contamination of 50C into 51C due to absorption of thermal photons and the detection efficiency difference between the 51 and 48 levels. All these errors can be encapsulated in a simple homographic law transforming the ideal P i g signal into the measured one, P g . For N repetitions of the experiment, the number of atoms actually detected in e or g, N d e,g can be written as a function of the ideal counts, N i e,g , as

N d e,g = A e,g N i e,g + B e,g N i g,e + C e,g .

Using N i g + N i e = N and after simple algebra, we get The Ω 0 value is very close to the expectation based on the cavity geometry (2π × 50 kHz). The slight discrepancy may be due to a slight offset of the laser beam with respect to the center of the cavity in the Oz direction. The p 1 value is in fair agreement with the expectations of an ideal cavity erasure procedure, followed by a return to thermal equilibrium (0.34 photons on the average) for 2.7 ms. The a, b, c, d values are also in qualitative agreement with independently measured detection errors.

We use the same set of parameters for all numerical simulations shown as red lines in the figures of the main text. These simulations are performed with the Qutip library [43], computing the interaction hamiltonian between a 2-level atom and the harmonic oscillator of the mode of the cavity, truncated to the first 50 Fock states. We take into account the finite lifetime of the cavity and of the atom, the original thermal state we start from, and the field relaxation towards the steady thermal state of the experiment.

We use Fourier transform to convert the Rabi oscillation data to the frequency domain and get straight vision of the photon number distribution. In oder to do that and obtain sharp peaks the experimental signal is first linearly interpolated with a 0.1 µs step. We then substract its mean value, symmetrize it with respect to t = 0, and complete the data with 0 until t = ±3000 µs. We finally compute the fast Fourier transform with a rectangle window.

FIG. 1 :

 1 FIG. 1: (a) Sketch of the experimental set-up (only one cavity mirror shown for clarity), with the axis conventions. (b) Timing of the Rabi oscillation sequence. The solid blue line depicts the atom-cavity detuning δ as a function of time (time scale is qualitative). The resonant interaction starts at t = 0 and lasts for a variable interaction time ti. The green bar depicts the mw injection in C starting at t = -71 µs and lasting t β , the red bar the exciting laser pulse at t = -25 µs. The circular state preparation sequence takes place between the laser pulse and t = 0. (c) Timing of the cat revival experiment. Same conventions as for (b). The two resonant interactions last for a fixed time ti = 60 µs and a variable time t i . They are separated by the variable delay t d during which |g is eliminated (purple bar) and an injection lasting tα is performed (green bar).

  FIG. 2: (a) Vacuum Rabi oscillation. Dots with statistical error bars joined by a thin black line: experimental probability Pg(te) as a function of the effective interaction time te. Solid red line: numerical simulation of the experiment. (b) Rabi revival in a 13.2-photon coherent field. Dots with statistical error bars joined by a thin black line: experimental probability Pg(te). Solid red line: numerical simulation of the experiment. The inset shows a zoom on the revival. (c) Fourier transform of the signal in (b) (black dots). The solid red line is a fit to a superposition of Gaussian peaks. The vertical dotted blue lines are at the expected Rabi frequencies for photon numbers, n, given on top. The blue solid line is the envelope of the Poisson photon number distribution for a 13.2 photon coherent state, adjusted in height to fit the abitrary scale of the Fourier Transform.

  FIG. 3: (a) Measurement of the displaced cat's parity proportional to its Wigner function. Pg(t e , α) as a function of α for t e = 68.5 µs (dots with statistical error bars joined by a thin solid black line). The solid red line results from a numerical simulation. (b) Quantum revivals in a displaced cat state. Pg(t e , α) as a function of t e for α = -0.60 (dots with statistical error bars and a connecting thin solid black line). The solid red line results from a numerical simulation. The additional revival around Tr/2 is clearly visible. (c) Fourier transform of the signal in (b) (dots with connecting thin black line). Vertical blue lines depict the expected Rabi frequencies. The absence of even photon numbers in the distribution is conspicuous.

Figure 3 (

 3 FIG. 4: (a)-(d) Rabi oscillations signals around Tr/2 for the probe atom in a displaced cat state for α = -0.60 and four delay times t d = 6, 86, 146 and 206 µs. Dots with statistical error bars and connecting thin line are the experimental data. The solid red line results from the simulation of the experiment. The fast decrease of the oscillation contrast reveals the cat decoherence.

  FIG. S2: Velocity distributions of the atomic beam, with (black dots with connecting thin black line) and without (red dots with connecting line) strobsoscopic selection. The red curve is renormalized to have a comparable scale to the other one. The errors bars represent the statistical dispersion.

  FIG. S3: Vacuum Rabi oscillation. Dots with statistical error bars joined by a thin black line: experimental probability Pg(ti) as a function of the real interaction time ti. Solid red line: fit on the theoretical model we use to get all unknown parameters.

  g -B g , c = A g -B g + A e -B e b = B g + C g , d = B g + C g + B e + C e (S14)We use the homographic law of Eq. (S13) to transform the ideal Rabi oscillation of Eq. (S10) and we fit the resulting signal onto the raw Rabi oscillation. The result of the fit is shown as a red solid line in Fig.S3. We extract from this fit all unknown parameters, namelyΩ 0 = 2π × (49.88 ± 0.03) kHz, x 0 = 1.72 ± 0.02 mm, p 1 = 0.094 ± 0.014 a = 1, b = 0.133 ± 0.02, c = 0.297 ± 0.07, d = 1.136 ± 0.08 (S15)

  1, we find that the atom-field state |Ψ(t e ) at the effective interaction time t e reads

	|Ψ(t e ) =	1 √ 2	[|ψ + a (t e ) ⊗ |ψ + c (t e )
					+ |ψ -a (t e ) ⊗ |ψ -c (t e ) ] ,	(S5)
	where
	|ψ ± a (t e ) =	1 √ 2	e ±iΩrte/2 e ±iΩte |e ∓ |g ,	(S6)
	|ψ ± c (t

e ) = e ∓iΩrte/4 |βe ±iΩte , (S7)