Measuring the environment of a Cs qubit with dynamical decoupling sequences

We report the experimental implementation of dynamical decoupling on a small, non-interacting ensemble of up to 25 optically trapped, neutral Cs atoms. The qubit consists of the two magnetic-insensitive Cs clock states $\left| F=3, m_F=0 \right>$ and $\left|F=4, m_F=0\right>$, which are coupled by microwave radiation. We observe a significant enhancement of the coherence time when employing Carr-Purcell-Meiboom-Gill (CPMG) dynamical decoupling. A CPMG sequence with ten refocusing pulses increases the coherence time of 16.2(9) ms by more than one order of magnitude to 178(2) ms. In addition, we make use of the filter function formalism and utilize the CPMG sequence to measure the background noise floor affecting the qubit coherence, finding a power-law noise spectrum $1/\omega^\alpha$ with $\mathit{\alpha} = 0.89(2)$. This finding is in very good agreement with an independent measurement of the noise in the intensity of the trapping laser. Moreover, the measured coherence evolutions also exhibit signatures of low-frequency noise originating at distinct frequencies. Our findings point toward noise spectroscopy of engineered atomic baths through single-atom dynamical decoupling in a system of individual Cs impurities immersed in an ultracold $^{87}$Rb bath.


I. INTRODUCTION
The tremendous progress in manipulating and measuring ultracold quantum gases makes them a versatile platform for quantum technologies. Electric and magnetic field control enable the precise creation of quantum superpositions for targeted applications in quantum computing [1], quantum simulation [2], and quantum sensing [3]. In all these applications, information is stored in the relative phase of such a quantum superposition in so-called quantum bits (qubits). These coherent superpositions, however, are usually very fragile because they are sensitive to external perturbations. In particular, even in highly specialized labs, quantum systems cannot be perfectly isolated from their environment, which will lead to dissipation and eventually to decoherence, i.e., the decay of a pure quantum state into a statistical mixture. In this regard, the reduction of decoherence and the extension of coherence times is a prime challenge in quantum science which will enable the broad applications of quantum technologies.
One of the main sources of decoherence is dephasing caused by time-dependent (classical) fluctuations. It is known that such dephasing can, in principle, be reversed by so-called dynamical decoupling, a coherent controlpulse method that was originally developed in the field of solid-state physics. The primary purpose of dynamical decoupling is the effective reduction of the coupling between the qubit and a mostly generic environment by a sequence of pulses. These single [4] and multi-pulse sequences [5][6][7][8] have already been used to extend coherence in various physical systems, ranging from spin ensembles [9], semiconductor quantum dots [10,11], nitrogenvacancy centers in diamond [12,13], superconducting * widera@rptu.de qubits [14], ensembles of ultracold atoms and ions [15][16][17] to the limit of single atoms and ions [18][19][20]. Consequently, many key characteristics and possible causes of qubit decoherence in ultracold atomic systems have already been discussed and explained in various experimental and theoretical works [21][22][23][24][25][26].
Remarkably, such sequences cannot only be utilized to partially reverse the decoherence dynamics but also to perform noise spectroscopy, i.e., to measure details of the coupled, unknown degrees of freedom of the environment [14,20,[27][28][29][30][31]. This first appears to be counter-intuitive since sensing is usually more effective if the coupling between the probe and the measurement object is strong. The idea of noise spectroscopy through dynamical decoupling, by contrast, is the modification of the spectral properties of the qubit which probes the environmental noise spectrum in different frequency ranges. While one could expect that the measurement of all environmental effects is a hopeless task, the most relevant effects can be directly extracted from the detailed knowledge of the qubit's dynamics. This creates the possibility to recover details of the environment by controlling and measuring the qubit's time evolution. From a fundamental point of view, the dynamical driving of the qubit can be used for the creation of correlated states between the environmental and the qubit's degrees of freedom, closely related to spin-boson [32], polaron [33,34], and Kondo physics [35,36]. From a more practical point of view, the gained knowledge about the environment can be utilized to improve the experiment by applying optimized sequences tailored to the situation in the lab.
In this paper, we exploit the dynamical decoupling method to determine essential properties of the environmental noise spectrum affecting the qubit coherence in a small ensemble of optically trapped, neutral 133 Cs atoms. We use the two magnetic-insensitive Zeeman levels |F = 3, m F = 0 and |F = 4, m F = 0 of the 6S 1/2 ground state as qubit states and couple them by mi-crowave radiation. Such atoms have been immersed into an ultracold atomic gas recently [37] for probing coherence and dephasing dynamics of individual Cs impurities beyond the effect of electromagnetic fields. In this work, we demonstrate how dynamical qubit control can give access to essential properties of a coupled environment. Our results represent a subsequent stepping stone towards utilizing the Cs qubit as a quantum probe for correlated environments consisting of interacting quantum particles.
The paper is structured as follows. Sec. II comprises a description of the experimental setup and a detailed characterization of the qubit. In Sec. III, we introduce the theoretical description of the dynamically driven qubit in a decohering environment. Then, in Sec. IV, we present the experimental realization of the dynamical decoupling sequence and a comparison to theory. After that, we conclude our results in Sec. V and provide a brief outlook.

A. Experimental setup
A schematic overview of the trapping geometry in our experimental setup is given in Fig. 1. The preparation of the Cs atoms is a three-stage process. First, about 50 Cs atoms are captured in a high-gradient magnetooptical trap and loaded into a crossed optical dipole trap formed by the intersection of two 1064 nm laser beams. Subsequently, the Cs atoms are further cooled to a temperature of T Cs = 5.9(3) µK using degenerate Raman sideband cooling (see Refs. [37,38]  to the cooling effect, Raman cooling optically pumps the Cs atoms to the absolute ground state |F = 3, m F = 3 . The Cs atoms are then transferred to the magneticinsensitive state |F = 3, m F = 0 by four successive microwave Landau-Zener sweeps [37]. We remove residual Cs atoms that are not in |F = 3, m F = 0 from the trap by a subsequent cleaning scheme, which is based on the spin selective readout in Ref. [39] and consists of a combination of microwave Landau-Zener sweeps and resonant laser pulses. This preparation process eventually leaves about 15to25 Cs atoms in the desired state |F = 3, m F = 0 and, at most, two Cs atoms in states |F = 3, m F = 0 [40]. In the last step, a one-dimensional optical lattice formed by two counter-propagating 790 nm laser beams is used as a conveyor belt to transport the Cs atoms to a second crossed optical dipole trap located at an axial distance of about 208 µm. Here, a coherent superposition of the Cs clock states is prepared, as explained in the subsequent section below, and the dynamical decoupling sequence is applied. Importantly, the Cs atoms are always trapped in the static optical lattice during the sequence. The lattice creates a repulsive potential for the Cs atoms along the horizontal direction and thereby freezes the atoms' position in this dimension. This results in trapping frequencies of ω ax = 2π × 60.7 kHz and ω rad = 2π × 694 Hz in the axial and the radial direction, respectively. Moreover, it ensures that Cs-Cs interactions are negligible since, at most, one Cs atom can be trapped in each lattice site. At the end of each dynamical decoupling sequence, a resonant laser pulse removes all Cs atoms in |F = 4 from the trap, and the optical lattice is used for spatially-resolved fluorescence imaging of Cs atoms in |F = 3, m F = 0 .

B. Qubit characterization
We start with a detailed characterization of the qubit transition before providing information about the dynamical decoupling sequence in Sec. III. For all measurements shown in the following, the external magnetic field at the atoms' position is calibrated to a value of 198.5 mG (see Appendix V for further details on the magnetic field calibration).
The qubit's transition frequency is measured using microwave spectroscopy, taking care to avoid power broadening. The corresponding resonant Rabi frequency Ω R = 2π × 619(2) Hz is independently measured, and Fig. 2(a) shows the spectrum of a microwave square pulse with τ = 725 µs duration. The signal follows the typical shape of the spectrum of a Fourier-limited square pulse where Ω 2 = 4π 2 (δν 0 −δν 0,res ) 2 +Ω 2 R , and A and C are the amplitude and the offset of the fit, respectively. The peak frequency is shifted from the exact atomic resonance frequency ν 0 = 9.192 631 77 GHz by δν 0,res = −133(12) Hz. The two main contributions to this shift are the differential light shift δ DT of the trapping field and the quadratic Zeeman shift δ B . We find that the measured shift is in good agreement with the expected shift δν 0,theo = −123 Hz calculated from experimental trapping parameters and the Breit-Rabi formula. For dynamical decoupling sequences with multiple (π-)pulses, the pulse durations are important timescales determining the duration of the sequences and must be short compared to all other experimental timescales. We achieve for all dynamical decoupling measurements a resonant Rabi frequency of Ω R = 2π × 21.14(4) kHz [ Fig. 2(b)], which corresponds to a π-pulse duration of τ π = 23.65 µs. The Cs clock transition chosen is electricdipole forbidden so that the excited state is long-lived, and T 1 decay occurs on the timescale of seconds and will therefore be neglected. Another important timescale is the lifetime of the Cs atoms in the optical lattice, which is limited by heating due to the relative phase noise between the two laser beams to about 685(81) ms. Moreover, a sufficient timescale separation between coherence time and π-pulse duration is required for realizing multi-pulse sequences. To measure the coherence time in the absence of dynamical decoupling, we apply a Ramsey sequence consisting of two microwave π/2-pulses with variable sequence length T [ Fig. 2(c)]. The microwave frequency is set to the unperturbed transition frequency ν 0 , i.e., detuned by approximately 133 Hz from the actual qubit resonance so that several fringes can be observed within the decay time. However, the pulses can still be approximated as near-resonant, as the detuning is much smaller than the bare Rabi frequency. For such a sequence, the expected population dynamics can be described by where A and C are the amplitude and offset of the fit. Previously, this functional behavior was a heuristic assumption that describes accurately the dynamics of the experiment. Later in this paper, in Sec. IV, we will provide a justification of this formula which is a consequence of the measured noise spectrum. We fit Eq. (2) to our data in Fig. 2(c) and find a frequency shift of |δν 0,res | = 136.8 (8) Hz, close to the frequency shift measured through microwave spectroscopy, and a coherence time ofτ c = 15.7(7) ms. A detailed study on inhomogeneous and homogeneous dephasing effects causing decoherence in an ensemble of Cs atoms confined by a far-off-resonant standing-wave optical dipole trap can be found in Ref. [22]. Inhomogeneous dephasing is mainly caused by inhomogeneous differential light shifts originating from the thermal motion of the atoms in the trapping potential, which lead to slightly different resonance frequencies among the ensemble. The effect of this dephasing mechanism can be reversed by applying a spin-echo sequence comprising an additional π-pulse between the two Ramsey π/2-pulses. Homogeneous dephasing effects, instead, lead to homogeneous, time-dependent fluctuations of the transition frequency among the ensemble. The main causes of homogeneous dephasing in our experimental setup are intensity fluctuations of the dipole trap laser and heating of the atoms (e.g., due to relative phase noise between the two lattice laser beams). Fluctuating magnetic fields are another source of broadening which is, however, far less pronounced because of the use of magnetic-field-insensitive states. Additional homogeneous dephasing sources, that play a minor role in our experimental setup, are the pointing instability of the dipole trap laser and imperfections of microwave pulses. Dephasing caused by homogeneous broadening cannot be simply reversed by the spin-echo sequence. Here, we need to apply dynamical decoupling sequences, which extend the spin-echo scheme by using multiple π-pulses, to compensate for noise originating at different frequencies. The theoretical background for this is presented in the following section.

III. QUBIT COHERENCE IN A NOISY ENVIRONMENT
We use a semiclassical model to describe the dynamics of the driven qubit coupled to environmental noise [41,42]. Hereby, we study the dynamics of the density matrix ρ of the qubit, which is governed by the von Neumann equation The corresponding Hamiltonian iŝ and reported in the frame rotating with the qubit eigenfrequency. Here, we have introduced the Pauli matriceŝ σ i with i = x, y, z of the two qubit states |F = 3, m F = 0 and |F = 4, m F = 0 . The driving of this qubit is described by the coherent, time-dependent Rabi frequency Ω (t). The influence of the environment is encoded in fluctuations described by the stochastic process β(t). We assume that this noise has a vanishing mean, β(t) = 0, and is determined by the time-correlation function The expectation values are taken over noise realizations. The stochastic fluctuations are characterized by the noise spectrum In our description, we have neglected noise-driven spin flips, which are expected to be relevant only on timescales much longer than we describe in the following. Moreover, in this work, we will consider a spectral density of the form with a positive real number S 0 > 0 and a positive real exponent α > 0, describing an algebraic decay of the spectral density. This spectral density describes a variety of noise processes in various physical scenarios including 1/f -noise. Similar to Ref. [24], we derive a master equation which is valid in the weak coupling approximation. For this, we transform into a frame which is rotating with the driving field The dynamics ofˆ is then governed by the master equation where we introducedĴ(t) =Û (t)σ zÛ † (t). This master equation can be mapped to the Bloch equations for the polarization p = (p x , p y , p z ) T = Tr[σˆ ], with σ = (σ x ,σ y ,σ z ) T . The Bloch equations are given by with the matrix elements of M(t, t ) All other matrix elements are equal to zero. We will consider the case where Ω (t) is a series of δ-shaped π-pulses applied at times t j , described by Ω (t) = N j=1 πδ(t − t j ). In this case, we find that sin(Φ) = 0 for all times. The initial state in this protocol is obtained after applying a π/2-pulse to the ground state. Therefore, the initial state is polarized in the y direction such that our initial conditions are p x (0) = 0, p y (0) = 1, p z (0) = 0. With these initial conditions, we can solve Eq. (8) and find p x (T ) = 0, p y (T ) = V (T ), p z (T ) = 0, where T is the total integration time and where we have introduced the visibility The visibility is determined by the coherence integral and the filter function For the sequence of π-pulses, we obtain Until now, we have neither specified the instances t j for the π-pulses nor the total integration time T . In fact, the result given in Eq. (17) is general and can be used for an arbitrary sequence of δ-shaped π-pulses, i.e., pulses with infinitely short pulse duration. In addition, one can also include a finite duration τ π of the π-pulses, which results in a slight modification of g N (ω, T ) that has been discussed in Refs. [14][15][16]42]. In this paper, we employ the N -Carr-Purcell-Meiboom-Gill (N -CPMG) sequence, whose π-pulses are applied at times t j = (2j − 1)T /(2N ) [5,6], and the 1-CPMG sequence corresponds to the spin-echo sequence. This pulse sequence consists of temporally equidistant π-pulses where the first pulse is applied at time t 1 = T /(2N ) and the last pulse at t N = T − T /(2N ) [see Fig. 3(a)]. We remark at this point that it has been shown that the filter function of the CPMG sequence is hardly modified by the effect of short π-pulses with non-vanishing duration, as used in our experiment (τ π /(T /N ) < 10 −1 ) [42]. This justifies that we can neglect the influence of the finite pulse duration, which we do throughout this paper. The choice of the CPMG sequence enables a very efficient decoupling from noise with a spectrum S(ω) given by Eq. (6), especially if we consider no high-frequency cut-off [15,43]. This results from the filter function g N (ω, T ), which is shown in Fig. 3(b) for different N -CPMG sequences and a fixed sequence length T . The filter functions' peaks shift to higher frequencies as the number of π-pulses is increased, resulting in different frequency ranges of the noise spectrum contributing to decoherence [Eq. (15)]. Consequently, since the spectral density S(ω) [Eq. (6)] is decreasing with ω, we obtain a lower total decoherence for an increasing number N of π-pulses. We mention, however, that we could also have chosen another pulse sequence such as the Uhrig sequence [7] which would show a similar performance in reducing the decoherence.
Besides the decoupling from noise, this band-pass filter property of g N (ω, T ) can also be exploited to perform noise spectroscopy. Here, we vary the number N of π-pulses and the sequence length T such that the filter function samples the environmental noise spectrum S(ω). From the measured decoherence signal, we can then reconstruct the spectrum. To demonstrate this, we first use that the sum in Eq. (17) for the N -CPMG sequences can be written with the help of a geometric series Here, we have used that the expression of g N (ω, T ) for the N -CPMG sequence only depends on u = ωT . The expression (18) exhibits for large N a sharp peak whenever u = (2k + 1)N π, where the denominator in the absolute value vanishes. For large N , we can approximate the filter function as an infinite series with the introduction of sinc(x) = sin(x)/x. Since the slope of the spectral density Eq. (6) where we have defined the coherence time using and the Riemann ζ function ζ(s) = ∞ n=1 1/n s . This finding is in agreement with Ref. [41].
Equation (20) shows that the exponent α can directly be inferred from the visibility V (T ) [27,28,44]. The latter exhibits an exponential behavior where the exponent of the argument is 1 + α, determined by the algebraic behavior of the spectral density. Moreover, the coherence time τ c [Eq. (21)] is algebraically growing with the number N of pulses. This increasing coherence time can be understood by the shift of the filter function's peak towards higher frequencies for an increasing number N of π-pulses [see Eq. (19) and Fig. 3(b)]. The latter leads to a reduction in the net integrated noise and thereby to an enhancement of the coherence time. The exponent of the growth is α/(1 + α) which is monotonically increasing with increasing α. Therefore, the noise reduction is more pronounced if the noise spectrum decays faster. For α → ∞, the coherence time is increasing linearly with N .

IV. N -CPMG MEASUREMENTS
In the experiment, the effect of a N -CPMG sequence on the coherence time is investigated by tracking the visibility of Ramsey fringes for different sequence lengths T . To measure a Ramsey fringe for fixed T , we scan the phase ϕ of the second π/2-pulse from −185°to185°in 10 steps and measure the Cs population in |F = 3, m F = 0 [ Figs. 3(a) and 3(c)]. We fit a sinusoidal function [37] to the data, where a, c and Φ are the amplitude, the offset, and the free phase of the fit, respectively. The fringe visibility is defined as The Ramsey fringes with total Cs atom number [ Fig. 4(a)] need to be normalized to calculate the fringe visibility according to Eq. (24). We rescale the data for a given number N of π-pulses and pulse sequence length T according to where the total number of Cs atoms contributing to the signal N Cs,tot (T ) depends on T due to the finite lifetime of the atoms in the optical lattice. For each data set with a fixed number N of π-pulses, the total atom number N Cs,tot (T ) is deduced from two additional measurements. First, we extract the maximal number of Cs atoms contributing to the signal N Cs,0 from a Ramsey fringe with N = 0 π-pulses and shortest possible pulse sequence length T = 27 µs [ Fig. 4(b)]. Second, we measure the lifetime τ LT of the Cs atoms in the optical lattice by monitoring the Cs atom number as a function of the pulse sequence length T without applying microwave radiation [ Fig. 4(c)]. We find typical lifetimes between τ LT = 539(37) ms and τ LT = 831(72) ms. The total atom number, which sets the scaling factor for the normalization, is then given by Each data point in Fig. 3(c) [4(a)] is an average of typically 12to15 experimental runs resulting in a total measurement time of about 18 hto24 h for each N -CPMG sequence. Figure 5 shows the measured visibility evolution for different numbers of π-pulses ranging from N = 0 to N = 10. The data points predict, in general, a decaying dynamics of the visibility, which is slower for a larger number N of π-pulses. To get a more quantitative access to this data, we fit to the measured visibility. For the calculation of χ N , we calculate the coherence integral in Eq. (15) with the expected form of the spectral density in Eq. (6) and the filter function given by Eq. (17) for every number N of πpulses. This assures that we also take into account small-N effects which are not present in the analytical result shown in Eq. (20). Importantly, each N -CPMG data set owns its specific amplitude V 0 and offset b. In contrast, the parameters S 0 and α are global fitting parameters shared among all N -CPMG data sets. The fits are visible as the solid lines in Fig. 5 and agree well with the data points. From our fit, we find an exponent α = 0.89 (2) and the value S 0 = 1288(122) s 1−α . In addition, we can also extract the coherence time τ c , which is defined as exp[−χ N (τ c )] = 1/e, from these fits. The values τ c are shown in Fig. 6 as data points in a color scheme that matches the one used for the visibility data shown in Fig. 5 (N = 0 black, N = 1 green,  etc.). In the absence of dynamical decoupling (N = 0), we find a coherence time of τ c = 16.2(9) ms. This value is in agreement with the coherence time ofτ c = 15.7(7) ms obtained by fitting the population dynamics in Fig. 2(c) with Eq. (2). We want to underline that we have assumed a heuristic decay with exp −(T /τ c ) 2 in Eq. (2). To validate this assumption, we can now compare this decay to the analytical result in Eq. (20) which predicts a decay with exp −(T /τ c ) 1+α ≈ exp −(T /τ c ) 1.89 . Consequently, the accurate result obtained from Eq. (2) can  be explained by the exponent α = 0.89 (2) which is close to one. For a growing number N of π-pulses, we obtain an increasing coherence time. In particular, we find for the maximum value of N = 10 a coherence time of τ c = 178(2) ms, which corresponds to an enhancement of more than one order of magnitude when compared to the N = 0 case. To compare these results with analytical estimates, we have calculated Eq. (21) from the fitted value of α and S 0 . The result is shown as solid gray line in Fig. 6 and agrees better for larger numbers of π-pulses. The analytical estimate predicts an algebraic growth with N α/(1+α) ≈ N 0.47 . Compared to the theory, the experimental results also suggest that the observed coherence time can be further extended by adding more π-pulses to the CPMG sequence since the Ramsey fringe visibility doesn't drop significantly. This also implies that imperfections in the π-pulses play a minor role in our system. Our analytical model in Eq. (21) predicts an unbounded growth of the coherence time, which is a consequence of neglecting the natural T 1 decay time of the qubit states. Dynamical decoupling sequences can, in principle, enhance the total coherence time only up to the τ c = 2T 1 limit if energy relaxation is taken into account [14]. In our experiment, energy relaxation occurs on the timescale of seconds, but the maximum number of π-pulses is eventually limited by the experimentally achievable minimum pulse duration and maximum sequence length. The former is determined by the available Rabi frequencies, and the latter is constrained by the lifetime of the Cs atoms in the optical lattice of about 685(81) ms. Consequently, we can never reach the physical limit τ c = 2T 1 and expect the coherence time to grow for all experimentally accessible N -CPMG sequences.
In contrast to the strictly monotonous decaying fits visible in Fig. 5, our experimental data also show narrow dips in the visibility for N ≥ 6. These dips are highlighted in the insets of Fig. 5 for N = 6, 8, 10. We expect that these dips originate from several resonances in the noise spectrum. In fact, already a single resonance gives rise to the appearance of multiple dips in the visibility. This can be seen by considering a spectral density which describes phenomenologically a resonance at ω 0 with a finite width of ∆ω. Using Eq. (19) and Eq. (28) in Eq. (15), we expect multiple dips at times T k (N ) = (2k +1)N π/ω 0 for k = 0, 1, 2, . . . . This makes it very difficult to determine from the experimental data whether the dips result from different resonances or from a single resonance. In particular, we could assume that the two dips in Fig. 5 for N = 8 found at T a and T b with T a ≈ 120 ms and T b ≈ 200 ms originate from the same resonance ω 0 . In addition, if they correspond to neighboring dips, we can assume T b /T a = (2k + 3)/(2k + 1) which can be compared to the data, T b /T a ≈ 1.67. This suggests that k = 1 is a good fit with T b = T 2 (N = 8) and T a = T 1 (N = 8). From this, we can estimate the frequency which is approximately ω 0 ≈ 2π × 100 Hz. In order to identify dominant contributions to the experimentally determined noise spectrum S(ω), we investigate common effects, such as intensity fluctuations of the dipole trap laser and magnetic field fluctuations, that have already been seen to cause homogeneous dephasing in ultracold atomic systems [20,22]. The intensity fluctuations are measured by shining the trapping light of the horizontal dipole trap onto a fast photodiode. We record the photodiode voltage as a function of time and calculate the power spectral density (PSD) S DT of our discrete-time signal. In Fig. 7(a), S DT is visible as the red curve. The resulting PSD S DT comprises two noise components where the first component describes power-law noise with a characteristic, positive exponentα and the second component describes white noise with a constant power spectral density S WN . We use Eq. (29) as a fit model for our data in Fig. 7(a) and find an exponentα = 0.8904 (2). This fit is visible as black solid line in Fig. 7(a) and agrees very well with S DT . Importantly, the fitted exponent coincides with α = 0.89 (2)  that some of the dips in the visibility data originate from a resonance of S(ω) at ω 0 ≈ 2π × 100 Hz. Remarkably, a resonance at ω 0 is found in S DT and can therefore explain such features in the visibility.
To quantify the magnetic field fluctuations, we cannot directly measure the fluctuations outside of the vacuum system because they differ from the fluctuations at the atoms' position. Instead, we record the current through the magnetic field coils as a function of time and calculate the PSD of the signal. The resulting PSD, S B , is shown as brown curve in Fig. 7(b). Fundamentally different from S DT , the PSD S B is rather flat with no clear decay. Such a noise spectrum would result in a very different evolution of the visibility. Therefore, we conclude that the background noise floor inferred from the measured coherence evolution is mainly caused by intensity fluctuations of the trapping light and not by magnetic field fluctuations. This finding highlights the magneticinsensitivity of the two Cs clock states which have been used as stable qubit states.

V. CONCLUSION AND OUTLOOK
In conclusion, we have studied the dynamics of a noisecoupled qubit realized in 133 Cs under the effect of dynamical driving. Hereby, we have focused on a dynam-ical decoupling protocol, in particular the CPMG pulse sequence, that has increased the coherence time of the qubit by an order of magnitude. Moreover, from the visibility data, we were able to analyze properties of the environmental noise spectrum including its algebraic decay. From this, we have calculated analytical estimates for the coherence time which are in good agreement with the measured data. To find the origin of the noise affecting the qubit dynamics, we have measured the PSDs determining the noise spectra of the dipole trap laser intensity and the magnetic field. With these, we could show that the main noise source originates from the trapping laser intensity while magnetic field fluctuations play only a minor role. The latter originates from the magneticinsensitivity of the used qubit states. Remarkably, the algebraic decay obtained from the visibility data is in excellent agreement with the algebraic decay found in the PSD of the dipole trap laser intensity. This highlights the ability to measure spectral properties of the environment with a dynamical probe.
In future work, we will investigate the dynamics of the 133 Cs qubit when it is coupled to a bosonic environment of 87 Rb atoms. The qubit dynamics could then be analyzed to potentially measure properties of an interacting environment. This is the next step towards the implementation of a single 133 Cs qubit as a versatile, nondestructive quantum probe for noise spectroscopy of engineered atomic baths. An important prerequisite for this is that the lifetime of the qubit without the bosonic bath exceeds the lifetime in the presence of the 87 Rb atoms. By comparing typical lifetimes obtained for the mixture [37], τ c,Rb ≈ 1 ms, which are one order of magnitude shorter than the shortest lifetime obtained in this paper, τ c = 16.2(9) ms, we conclude that this requirement is fulfilled. This will enable the distinction of dynamical features emerging from, e.g., noise in the dipole trap laser intensity and features emerging from the bosonic 87 Rb bath.
Besides its potential as a quantum sensor, this setup can also be used to explore the entanglement dynamics of the Cs qubit with the surrounding 87 Rb bath under dynamical driving. This creates the possibility to study non-equilibrium phenomena in quantum manybody systems closely connected to spin-boson and polaron physics. Additionally, it will be interesting to investigate the possibility of tailoring the qubit-bath coupling with periodic driving which opens the door to quantumstate engineering in open quantum systems.

DATA AVAILABILITY STATEMENT
The data that support the findings of this study are openly available at the following URL/DOI: https:// doi.org/10.5281/zenodo.7728140.
We employ microwave spectroscopy on the |F = 1, m F = 0 → |F = 2, m F = 1 transition of the 87 Rb ground state to precisely calibrate the magnetic field B exp at the atoms' positions by compensating for the contributions B bg of ambient fields. The Rb atoms are trapped in the second crossed optical dipole trap and are initially prepared in the |F = 1, m F = 0 state. The microwave frequency is fixed to the expected transition frequency given by the Breit-Rabi formula for the desired field magnitude B exp = 198.5 mG, neglecting the small differential light shift δ DT = 39 Hz of the trapping field. We vary the externally applied magnetic field B coil , and thereby also the transition frequency of the Rb atoms, and probe the excited state population by standard absorption imaging using a time-of-flight measurement [ Fig. A1(a)]. The microwave pulse duration τ = 243.9 µs is chosen such that full population transfer occurs as soon as the resonance condition B coil + B bg = B exp is fulfilled. The signal follows the typical shape of the spectrum of a Fourier-limited square pulse where Ω 2 = ∆ 2 + Ω 2 R , ∆ = 2π × (B coil − B coil,res ) × 0.7 MHz G −1 in the limit of small magnetic fields and A and C are the amplitude and the offset of the fit, respectively. Notably, the resonant Rabi frequency Ω R = 2π ×2.180(5) kHz is measured independently [ Fig. A1(b)]