Calibration of spin-light coupling by coherently induced Faraday rotation

: Calibrating the strength of the light-matter interaction is an important experimental task in quantum information and quantum state engineering protocols. The strength of the off-resonant light-matter interaction in multi-atom spin oscillators can be characterized by the readout rate Γ S . Here we introduce the method named Coherently Induced FAraday Rotation (CIFAR) for determining the readout rate. The method is suited for both continuous and pulsed readout of the spin oscillator, relying only on applying a known polarization modulation to the probe laser beam and detecting a known optical polarization component. Importantly, the method does not require changes to the optical and magnetic fields performing the state preparation and probing. The CIFAR signal is also independent of the probe beam photo-detection quantum efficiency, and allows direct extraction of other parameters of the interaction, such as the tensor coupling ζ S , and the damping rate γ S . We verify this method in the continuous wave regime, probing a strongly coupled spin oscillator prepared in a warm cesium atomic vapour.


Introduction
The off-resonant interface of light with atomic ensembles has been widely explored in the last decades [1][2][3][4] in ultra-cold, cold and warm alkali implementations.The spin degree of freedom in the atomic ground state coherences and its coupling to light has been used in protocols ranging from fundamental [5][6][7][8] to technological [9][10][11] applications.Furthermore, the collective spin excitations of the highly polarized atomic ensembles in a static magnetic field can be well approximated by harmonic oscillator-like degrees of freedom -a spin oscillator [12].The oscillator mapping is essential for the interface with nano-mechanical oscillators via back-action evasion [13] and entangling [14] measurements, which promises sensitivity improvements in future gravitational wave detectors [15] and optical quantum control of the hybrid spin-mechanical system [16].
In the interface between atoms and light, characterizing the strength with which the systems couple is paramount for understanding its dynamics.According to the principles of quantum mechanics, the statistical nature of the quantum measurement process leads to fundamental limits in estimation of system's state at a given instant of time [17].Optimizing the role of the measurement performed by the optical probe on the spin oscillator is key.
The optical readout of a highly polarized atomic ensemble prepared in a spin oscillator state contains contributions arising from [18]: optical shot noise, inherent to the quantum nature of light in the process of photo-detection; coherent spin state or ground state noise, from the zero-point energy required to satisfy the Heisenberg uncertainty principle; thermal noise, originating from the extra spin fluctuations of the oscillator having non-zero mean occupancy  S due to imperfect optical pumping; and lastly, quantum back-action noise, originating from the perturbations of the optical readout in the oscillator's dynamics.As we obtain information about the oscillator by performing measurements on the light that has interacted with the system, it is of key importance to faithfully characterize the weight of each of these contributions.
The standard quantum limit (SQL) for a measurement of mechanical displacements, for example, sets the sensitivity to external fields in conventional interferometric measurements [19].At the SQL, the detection shot noise and measurement back-action contribute equally to the measurement imprecision.A highly efficient mapping of the oscillator state to light requires the quantum back-action to dominate over the coupling to the thermal environment [4].Common to these protocols is the importance of the interaction strength parameter between light and the oscillator, here defined as the readout rate Γ S (also commonly known as the measurement rate in the optomechanics community [20]).Another figure of merit that quantifies the efficiency of the coupling is the quantum cooperativity  q , here defined as  q = Γ S 2 S ( S +1/2) .Working in the regime  q 1, in which the coupling is strong, indicates that the measurement significantly influences the oscillator dynamics, allowing for its control and manipulation.
In this paper, we show how the parameter Γ S may be extracted from a measurement based on the interference of the induced Faraday rotation, i.e., the oscillator response to a classical optical polarization modulation, with the modulation itself.We call it Coherently Induced FAraday Rotation, or CIFAR, signal.The method further allows extraction of the effective linewidth,  S , and the tensorial part of the interaction,  S , describing the deviation from the idealized quantum non-demolition (QND) interaction.Crucially this method relies on the same alignment of magnetic and optical fields as well as optical pumping of the atomic ensemble, as used for experiments such as [13,14].The CIFAR method is applicable in all experimental implementations of spin oscillator, from ultra-cold to warm vapors, in ensembles with total angular momentum (per atom) equal or larger than 1  2 .The coherent drive also allows for probing the atomic motion through the laser beam [21] and characterizing the coupling to fast decaying spin modes [22], which give rise to a broadband spin response.We verify the CIFAR in the continuous regime, probing a spin oscillator in the strong coupling regime, prepared in a warm cesium atomic vapour.Lastly, we also present the limits of our linearized oscillator description by driving the system with large classical polarization modulation.
The technique here described is especially suited for continuous wave measurement of (single) spin oscillators, which should be contrasted with the mean value transfer method [10] and thermal noise scaling [23], which rely on time consuming measurements and more dramatic changes of the experimental setup.We nonetheless highlight that this technique can also be employed in the canonical state-preparation-probing experimental cycle.Furthermore, the signal depends only on the interference between the drive and response, it is independent of the overall detection quantum efficiency and thermal noise calibrations, as the calibration used in [13].The CIFAR method is similar to the Optomechanically Induced Transparency (OMIT)-response measurement technique [24], used to characterize optomechanical coupling parameters [25].

Theory
The interaction between an atomic ensemble and light has been widely studied in the context of optical pumping [1] and quantum information applications [4].For a far detuned monochromatic optical field with intensity much below saturation, the effective atom-light interaction is given by the coupling of the electronic-ground-state magnetic sublevels to the light polarization.The interaction can be seen as a mutual rotation of the light and spin variables, { Ŝ , Ŝ , Ŝ , Ŝ0 } and { Ĵ , Ĵ , Ĵ , Ĵ0 } respectively, in the form of polarization-dependent ac Stark shifts of the ground state levels and spin-state-dependent index of refraction, according to the atomic polarizability tensor [26,27].Typical ensembles consists of  ≈ 10 8 − 10 12 atoms, with the collective macroscopic spin being represented as Ĵ,, =  =1 F () ,, , where F () ,, are the Cartesian decomposition of the total angular momentum operator of a single atom.
The light-matter interaction, along with the contribution from an external bias magnetic field applied in the -direction, gives the spin Hamiltonian where the first term refers to the magnetic coupling induced bias magnetic field, shifting the ground state magnetic sublevels by ± S , with sign depending on the direction of the magnetic field with respect to the mean spin.The coefficients  0 ,  1 , and  2 as the relative weights of the scalar, vector and tensor contributions of the polarizability tensor [27], and  S is the single-photon coupling rate.The relative weights of the contributions depend on the level structure of the atom and can be controlled by the laser detuning from the atomic resonance.The vector and tensor contributions are related to circular and linear birrefringence of the atomic medium, respectively.The scalar component leads to a polarization independent phase shift.We now focus in the specific case of cesium [26].For a laser beam detuned by Δ from the  = 4 →  = 5 transition in the D 2 cesium line interacting with atoms in the  = 4 ground state manifold, the a  parameters are given by with Δ 35 /2 = 452 MHz and Δ 45 /2 = 251 MHz as the excited state hyperfine splittings between  = 3 and  = 5, and  = 4 and  = 5 [28], respectively.A detuning Δ > 0 (Δ < 0) corresponds to the case with laser frequency above (below) the  = 4 →  = 5 transition.In our experiments, the probe laser is detuned by Δ/2 = 3 GHz, where relative weights are  0 ∼ 3.83,  1 ∼ 1.05, and  2 ∼ 0.004.In practice, the ensemble is not perfectly polarized due to limited optical pumping efficiency and decay due to, e.g., wall collisions, natural lifetime and optical de-pumping.Since the values for   above are calculated for perfect spin polarization, the effective values observed experimentally differ somewhat from those stated here.Further, the imperfect spin polarization gives rise to a thermal, stochastic distribution of the spins in the different magnetic sublevels, which shows as thermal noise in the detection.
The interaction between light and the atomic ensemble in equation ( 1) can be simplified and linearized in the case of large ground state spin polarization and a strong polarized classical laser field.These approximations constitute the mapping of the spin system to an oscillator system.As we will describe in the next section, the ensemble is optically pumped such that   = Ĵ and Ĵ , Ĵ   at any instant of time.The probe is a strong classical field linearly polarized at an angle  to the magnetic field, which is also the quantization axis.The angle  controls the relative contributions of the vector and tensor effects described by the Hamiltonian (3).For simplicity, we change basis of the polarization variables such that the component parallel to the quantization axis, here Ŝ , describes the strong field as a classical variable with mean photon flux Ŝ = Ŝ0 =  , leaving Ŝ⊥ , Ŝ as zero-mean quantum variables.Mathematically, we rotate the polarization variables around the Ŝ components as { Ŝ = Ŝ cos 2 − Ŝ sin 2, Ŝ⊥ = Ŝ sin 2 + Ŝ cos 2, Ŝ , Ŝ0 }.
For a highly polarized ensemble in the  = 4 hyperfine manifold, the Hamiltonian in equation ( 1) can be simplified [29].In the limit of high steady state spin polarization, where only the two extreme magnetic sublevels, i.e. either   = +4, +3 or   = −4, −3, are populated, we perform the Holstein-Primakoff approximation [30] ), PL ( )] = (/2)( −  ) and [ XS , PS ] = .The quantity sgn(  ) refers to the sign of the mean spin, being positive (negative) for the negative (positive) mass oscillator cases.Notice that the sign of  S carries information about the mutual orientation of   and the external bias magnetic field.In the harmonic oscillator language, the mutual orientation defines the effective mass of the spin oscillator, with − S (+ S ) referring to the negative (positive) mass.In the derivation of equation ( 3), we have omitted constant energy terms, as they do not affect the dynamics of the variables of interest.
The parameters are the spin oscillator readout rate and the tensor coupling parameter, respectively.If  S = 0 the light-spins interaction is of the Quantum Non-Demolition (QND) type.In our experimental regime, as  S ≠ 0, the spin-light interactions deviates from the QND interaction, allowing for dynamical cooling/heating of the spin oscillator and changing the total decay rate and effective bath occupation in similar fashion to the effects of light interaction with a mechanical oscillator in the field of optomechanics [31].Experimentally, due to imperfect spin polarization, we will have  S smaller than the value predicted from equation (5).As shown in Section 3, our full model with  S as free parameter describes the measured response sufficiently well.
The spin system evolves coherently due to the Hamiltonian given in equation ( 3), and incoherently due to spin decay and coupling to an external effective spin bath [32].Furthermore, atomic motion leads to a time-dependent light-spin coupling.There is, in principle, an infinite set of collective spin modes that evolve in time accordingly to the ensemble geometry, collisions, dephasing, and diffusion [22].Here, we focus on the so called flat spin mode corresponding to the total spin Ĵ =  =1 F ()  , the mode which is the most resilient to motional dephasing as it is fully symmetric with respect to shuffling atomic positions.In the linearized language introduced above, we assign effective position and momentum variables, XS , PS to this mode.Later on, we will also introduce effective variables to describe the dynamics of the fast decaying spin modes, here denominated as the broadband response, and also introduce a qualitative model that describes its response to light.
The dynamics of the spin variables due to the Hamiltonian presented in equation ( 3), for the case with effective positive mass, is with  S /2 =  S0 /2 +  S Γ S as the dynamical damping rate, including tensor effects [33].Here,  S0 includes the natural (in the dark) decay rate, and laser induced contributions (from the pumping and probing lasers).For notation purposes, we write the light and spin variables in the matrix form as where the superscripts in (out) denote the optical mode before (after) the interaction with the spin oscillator, to be presented below.In the CIFAR experiments, as we will discuss in Section 3, the oscillator is coherently excited, e.g. with a drive  in L ∝ sin  RF .Since the system is driven by a classical driving field, the system response can also be considered as a classically, and we drop the operator description from here on.
Given the linear system of equations ( 6) and a sinusoidal drive input X in L , a solution in the complex plane can be found using the ansatz  S () =  S ( RF ) −  RF  ,  S () =  S ( RF ) −  RF  , where  S ( RF ) and  S ( RF ) are complex numbers.We write this as where the matrices L and Z parametrize the interaction dynamics as with as the spin susceptibility.All spin and light variables are understood to be functions of drive frequency,  RF , the notation of which we suppress from now on.
The output light field, after the interaction with the spin oscillator given in equation ( 8), is where 1 2 is the 2 × 2 identity matrix.The equation above shows that the output light field will have two contributions: one directly from the input field and another from the response of the spin oscillator to the input.Having the light field as a common source, these two contributions can interfere.By inserting ( 9) and ( 10) into (11), we get the expressions for output optical quadratures after the interaction with the spin ensemble In general, we are able to select arbitrary input X in L and the detection X out L components by controlling the ellipticity of the polarization state by the phases  and , respectively.The input light state, without loss of generality, is assumed to be generated from a pure phase modulation  = | |   that, when referenced to a local oscillator (LO) in a Mach-Zehnder interferometer, as we have in Figure 1, can be arbitrarily decomposed into polarization variables and effective input amplitude and phase quadratures.Here, by convention, we have chosen  in L = ,  in L = 0 for  = 0. Path length difference control in the Mach-Zehnder interferometer allows for mixing the drive components via a basis rotation, and polarization homodyning allows for selecting the detection quadrature, as By inserting the equations (13) in equation ( 12), we get to the final form of the Coherently Induced Faraday Rotation (CIFAR) signal.We typically define the measured quadrature as  det L , such that the absolute squared of the detected spin response to an arbitrary input optical modulation is This equation is the main result of this section, being applicable to the description of the flat spin mode response to light.In a broader view, it becomes necessary to consider other spin modes which in contrast to the total spin will have some spatial dependencies.We consider collective operators corresponding to the transverse spin components of mode  given by Ĵ , = √   =1   (x  ) F, where   (x) represents the spatial shape of the spin mode,  is the volume of cell (for the purpose of proper normalization) and x  is the position of -th atom.The coherent evolution of each mode (collective operator) is then governed by equation ( 3) with the readout rate Γ  S now taking into account the overlap between the spin mode and the Gaussian light mode  G  , such that Γ  S ∼ | G  | 2 [22].At the same time, the incoherent part will depend on motion and wall collisions.It has been shown that even in the paraffin-or alkene-coated cells [34][35][36][37] or cells with dilute buffer gas [38] the atomic motion can be effectively described by the diffusion equation with adequate wall boundary condition.The cells with coated walls will feature a slow decay for the flat mode ( 0 (x  ) = 1/ √ ) which depends on intrinsic dynamics and wall decay, and much faster decay for all other modes, which given by   S =   2  , where  is the effective diffusion constant and   is the characteristic wavenumber of -th mode [22].
For the case of quantum noise, it becomes necessary to consider both thermal contributions of each mode (added incoherently), and the coherent interaction of each mode the the Gaussian beam.These broadband spin noises affect the atomic ensembles serving as magnetometers [39] or quantum memories [22].In our case we can consider only the coherent interaction and thus each spatial spin mode responds to the same light modulation.Therefore, we may simply modify equation (11) to the multimode case: with L  contains the susceptibility with the respective decay rate   S .In our case we shall work in a two-mode approximation for which the zeroth mode is the flat, long-lived mode, and the other mode has  S /2 ∼ 1 MHz.This is a justified approach as all broad modes contribute a similar flat background around the resonance which we primarily study here.Following this approach we obtain a two-component CIFAR signal in which the narrow (with linewidth  S and readout rate Γ S ) and broad parts (with linewidth  S,BB and readout rate Γ S,BB ) of the response may interfere according to their phase relation.In the following sections we will nevertheless give simplified formulas for the single mode case to facilitate understanding, and use the two-component model for fitting.
Given the input and detection angles, as well as the spin oscillator parameters and coupling to light, the CIFAR signal exhibits a characteristic frequency response.For developing intuition about the response, let us focus initially in the special case of  = 45°and  = 0°, corresponding to detecting the phase quadrature of light  out L and driving with an equal superposition of amplitude  in L and phase modulation  in L .For the choice of phases described above, the detected signal goes as Notably, the drive, represented by the constant term, and the spin response, proportional to the susceptibility  S ( RF ), are added coherently and interfere.In particular, the readout rate Γ S plays an important role in the interference pattern, modulating its strength.In the high-Q limit ( S  S ) and around resonance ( RF ∼  S ), the spin susceptibility is  S ∼ − S0 / S , for  S0 = 1 2 ((Δ RF +  S /2)) −1 , with Δ RF =  RF −  S as the detuning between the spin resonance and the input modulation tone.In this limit, the equation ( 16) becomes For exemplifying the procedure to extract the readout rate parameter Γ S , let us consider two specific cases.First, we analyze the case of  S = 0, that is, the light-matter interaction is of the QND type.Here,equation ( 17) reduces to The CIFAR 0 signal is a combination of a constant, a Lorentzian, and a dispersive term, representing the drive, the spin response and the interference between the drive and response, respectively.
Importantly, the minimum and maximum of the signal are separated by ∼ √︃ Γ 2 S +  2 S ∼ Γ S in the limit of high coupling, Γ S  S .The readout rate can thus be extracted just by noting this frequency difference, directly from the sweep figure.
For the second specific case, when  S ≠ 0, equation (17) leads to a correction of the separation, as the maximum and minimum are separated by ∼ √︃ (1 +  2 S ) (Γ 2 S (1 +  2 S ) +  2 S − 2Γ S  S  S ).In the high-coupling limit, Γ S  S , this simplifies to ∼ Γ S (1 +  2 S ).Having derived the needed expressions, we now turn to an experimental investigation of the CIFAR signal under different situations.

Experimental implementation
The experimental setup is depicted in Figure 1.We start by describing the atomic spin ensemble.The atomic spin ensemble is a warm gas of cesium atoms, confined to a spin anti-relaxation-coated microcell [40] with a 300 µm × 300 µm cross-section and 10 mm in length.The system is probed with a Gaussian beam that has a waist of  0 ∼ 80 µm in radius (1/e 2 ).The sub-millimeter transverse dimensions of the cell allow for fast motional averaging [21], ensuring an integrated interaction between all atoms and the light.The microcell is positioned in a magnetic shield which contains coils producing a homogeneous magnetic bias field  in the -direction.The strength of the bias field splits the magnetically sensitive Zeeman levels by | S |, i.e., the Larmor frequency.Here, the Larmor frequency is in the range  S /2 ∼ 0.3 MHz to 1.5 MHz.The intrinsic (in the dark) spin damping rate at 59 °C and 1.5 MHz Larmor frequency is  0,dark /2 = 450 Hz.Circularly polarized resonant optical pumping lasers, co-aligned with the magnetic field, spin polarizes the ensemble [41] to ∼ 80%.
The remaining part of the setup in Figure 1 sets the phase sensitive control of the excitation and detection of the optical signal.We effectively generate an arbitrary polarization of the input light by combining two laser beams, here named local oscillator (LO) and modulation, with orthogonal polarizations on a polarization beam splitter (PBS) with a phase delay .The modulation beam is phase modulated at a RF frequency  RF in a fiber electro-optic modulator (EOM) according to  EOM e   e  sin  RF  , for  as the modulation strength.An arbitrary optical polarization state can be set by choosing a given  and relative intensities of the beams.One of the output ports of the PBS is sent to the spin ensemble and the other is used for a polarization detection setup (phase lock).We now describe the input polarization to the spin ensemble.
The phase  controls the relative contributions of circular and diagonal components, represented by  in  and  in  , in the input polarization state.The DC (static) components will lead to a small rotation in the local oscillator's polarization; the AC components will induce CIFAR signal.In the linearized quadrature language, the effective AC input drive is written as (− in L sin  +  in L cos ) sin  RF .In the port that is directed to the phase lock output, the DC frequency interferometric signal, given in equation ( 19) is used for stabilizing , the path length difference between the LO and the EOM arms, and therefore the input modulation to the spin ensemble.Deviations from the ideal phase shift induced by birrefringent elements, e.g., half wave plates in Figure 1, leads to a further mixing between  in  and  in  in equation (19), complicating the calibration of .The same argument applies to setting the detection angle .Experimentally, we first drive the spin oscillator with a RF magnetic field [41] to set .The detection half-wave plate is set to give  the maximum photo-detected response to the applied magnetic field, which happens at  = 0°.Subsequently, we switch to a polarization modulation drive and find the effective  = 90°when the spin response has a Lorentzian line-shape.
The polarization modulation at frequency  RF leads to a phase coherent interaction between the oscillator and light according to equation (12).The signal is recorded by balanced polarimetry photodetection and processed by a lock-in amplifier phase-referenced to the drive, allowing us to extract the slowly varying amplitude  and phase components.In a experimental protocol very similar to the one used in continuous wave Magneto-Optical Resonance Signal measurements [41], scanning  RF around the resonant frequency  S at a rate much smaller than the spin damping rates, we ensure the steady-state performance and extract the signal of interest.To extract the useful parameters from the data, we implemented non-linear optimization and curve fitting routine to a two spin-modes model based in equation (15).
For experimental implementations that operate in the pulsed regime [5,7,10], in which the probing follows a spin-state preparation stage, the CIFAR signal can be extracted in a similar manner to the one prescribed in the continuous readout.For example, during a single repetition, a fixed-frequency polarization modulated probe can map out the time-evolving signal, thus obtaining a single point in the interference signal.Repeating the experiment with different drive frequencies allows for mapping a signal similar to those presented in the Results section.The data analysis therefore borrows the analysis discussed in the Theory section.

Results
We will now present experimental support to the model described in Section 2. We fit the CIFAR model given in equation ( 14) to the recorded data, present its performance on different Fig. 3. Scaled fit residuals from Figure 2. Residuals between the model and data, both in CIFAR amplitude (left column) and phase (right column), for the various drive voltages shown in Figure 2. In the right-most column we show the histogram of the residuals along with a unity width Gaussian curve (dashed lines) to guide the eye.We also print the reduced  2 .Some outliers are not shown.
experimental conditions and discuss the overall validity of the model.We start by the studying the response of spin oscillator to increasing modulation amplitudes .In the data present in figures 2 to 4, we fix  = 45°, the probe power at 500 µW, and cancel the non-linear quadratic Zeeman shift with tensor Stark shifts [43] by setting  ∼ 60°.In Figure 2, for each modulation amplitude, we record 3 scans and show the average as points, and the statistical 1 uncertainty error bars (vertical bars).We double the amplitude starting from 31 mV (in blue) until 500 mV (in purple), showing |CIFAR| (top panel) and the phase response (bottom panel), the amplitude and phase of the CIFAR signal, respectively.The grey line is the response of the spin oscillator to a shot noise-limited drive, in which the coherent polarization modulation is turned off.We see that the amplitude of the CIFAR signal follows the drive increase, doubling as the amplitude doubles.As the drive amplitude increases, the coherent response dominates the signal and the spread around the mean values decreases.
The average traces for each drive amplitude in Figure 2 is fit by the CIFAR model.The fits are displayed as solid lines, showing the good agreement to the measured amplitude and phase data.In the Figure 2 inset, we show the readout rate Γ S /2 returned by the fitting routine as a function of the drive amplitude.For the 31 mV drive amplitude, the value for the readout rate is Γ S /2 = 10.685+0.008 −0.18 kHz.For increasing drive, nonetheless, the trend is that the fitting routine returns smaller values.
The asymmetric parameter errors are obtained with the conf_invertal function of the lmfit Python package.The function returns the parameter values for which  2 =  2 min + 1, i.e., the interval containing the usual 68.27% probability, which for a Gaussian parameter error corresponds to the 1 uncertainty.Similar results was obtained by Markov Chain Monte Carlo optimization (not shown).The asymmetry arises due to a strong correlations between Γ S and other fit parameters, mainly the parameter describing the overall response ∝ .
We further note that successfully fitting the model to the data relies on reliably ascribing errors to the individual data points; the individual traces spans 2 orders of magnitude, and failing to account for this in the optimization routine leads to discrepancies in either the peak or valley of the traces.Curiously, the data errors largely inherit the shape of the undriven atomic ensemble, i.e., a Lorentzian centered on the spin frequency (not shown).This places the condition that to obtain good fit values, all measurements must be repeated a number of times, to obtain trustworthy statistics.
Looking at the fitting residuals for the different traces, presented in Figure 3, we see that structured deviations between model and data appear as the spin oscillator is driven with larger amplitudes.The residuals to the traces with drive voltages above 250 mV, shown in red and purple points, present significant deviations from our model.We believe that the deviations from our model appear as we start to drive the spin system significantly away from the small oscillation amplitude limited approximation that leaves the system beyond the linearized oscillator model.We have, therefore, experimentally found the limit on the drive strength that our model can describe and the driving amplitude regime that must be used to return results free of systematic effects.
In Figure 4 we present the dependence of the readout rate on   , the mean spin length.According to equation (4), Γ  ∝   .The spin length is controlled by the temperature of the vapor cell, which sets number of atoms.When heated, the cesium vapor pressure increases [28].For the data on Figure 4, we record CIFAR scans while the cell is heated from ∼34 °C (blue points) to ∼59 °C (purple points).The solid lines are fits to the data, with the frequency axis shifted according to Δ RF and re-scaled to the returned spin linewidth  S .The extracted readout rate increases from 1.1 kHz to 10.0 kHz.As the temperature, and consequently   , is increased, the peak signal increases and the minimum is shifted up in frequency.Importantly, the frequency detuning Δ RF for which the CIFAR signal is minimal follows the readout rate Γ S , as shown in the inset.There is, approximately, a one-to-one correspondence between Γ S / S and the frequency of the CIFAR signal minimum value, as shown by the line with slope 1 (solid line).Therefore, by choosing the input modulation type  = 45°, an approximate readout rate can be easily extracted from the CIFAR signal as the frequency difference between the maximum and minimum of the trace.We also note that at the highest temperature setting (∼59 °C), with linewidth  S0 /2 = 1.3 kHz and estimated spin thermal occupation  S ∼ 0.75 [44], we have Γ S / S ∼ 7 and estimate  q ∼ 3, indicating that the spin oscillator is strongly coupled to light.
In Figure 5, we present the CIFAR signal for different strengths of the tensor coupling parameter  S .For a given detuning from the atomic resonance Δ, it is modified by selecting , the angle between the LO linear polarization and the magnetic field  direction.According to equation ( 5), the angle  = 45°turns off the tensor coupling.For this experiment, we reduced the spin resonance frequency to  S ∼ 400 kHz to avoid non-linear Zeeman splitting [41], the probe power was set to 250 µW to reduce probe-induced power broadening, and the temperature to  = 55°.In Figure 5, we show the amplitude of the CIFAR signal for  = 45°(top panel) and  = 90°(bottom panel).The data for  = {0°, 45°, 90°} are shown in blue, orange and green dots, respectively.The choice of  = 45°gives responses similar to those presented in Figure 2.For this data set, we have Γ S = 4.9 kHz.The setting  = 90°, nonetheless, gives a rather different picture.According to equation (14), the detected signal goes as where in the last passage we used the high-Q ( S  S and  RF ∼  S ), and the small tensor coupling ( S 1) limit.For this configuration, the CIFAR is dominated by the constant term, since we mostly detect the input modulation.Near the spin resonance, the oscillator will add ( S < 0) or remove ( S > 0) signal according to the tensor coupling sign.The obtained tensor coupling parameters are  S = −0.045± 0.002 and  S = 0.040 ± 0.003 for  = 0°and  = 90°, respectively.For reference, the expected tensor parameter for a perfectly spin polarized ensemble is | th S | = 0.053.For  = 45°the spin contribution is, according to our theory, null; the returned value is  S = 0.000 ± 0.001.
In our last study we present the broadband spin contributions to the CIFAR signal.The measurements presented in Figure 6 are taken in the same experimental conditions as the data in Figure 2, but now scanning the drive tone in a ∼ 600 kHz band around  S .The CIFAR signal amplitude (top panels) and phase (bottom panels), including the model fits, are shown for  ∈ {−45°, 0°, 45°} in blue, orange and green, respectively.Apart from the symmetric changes in the response as  is changed from −45°to 45°, the  ∼ 0°amplitude and phase responses display characteristic features of a broadband spin response.The broadband spin response can be clearly seen by setting Γ S = 0 (dashed orange line) or in our full model fit (dark orange line).The light orange line corresponds to the predicted response of the spin oscillator in the case Γ S,BB = 0.The broadband spin response, having a damping rate  S,BB /2 = 0.93 MHz, couples to the drive with rate Γ S,BB /2 = 33.4kHz, distorting the phase response and adding a pedestal to the detected amplitude.Remarkably, although having a complex origin [22], the broadband response is qualitatively well described by a single effective mode.

Conclusion
In summary, we have presented a novel approach for calibrating the light-matter interaction between off-resonant optical beams and collective spin systems, the CIFAR technique.Experimentally, the calibration method relies only on applying a known input modulation and detecting a known optical quadrature, the variables parametrized by  and , respectively.Supplied with the input-output relations of the spin-light interaction, a simple procedure for determining the interaction parameters, among those most importantly the readout rate Γ S , is described.Fitting the recorded signal to the full model provides a full characterization of the system parameters.The technique does not rely on knowing the photo-detection efficiency or the ensemble spin polarization.We have verified the good agreement between data and the CIFAR method by continuously probing a strong coupled spin oscillator prepared in a warm cesium atomic vapour.
Theoretical refinement of our model can be envisioned with the consideration of the full Zeeman structure of the ground state manifold [45] ( = 4, in the present example).With that, the model will be able to account for the non-unity spin polarization and return consistent values for the tensor parameter.Furthermore, in the case of large excursions by the transverse angular momentum variables induced by the drive light, going beyond the linearized regime can account for the mismatch between theory and data shown in our residuals analysis.
The technique here presented is also a powerful method for studying the coupling of light to the spin modes under diffusion and spatial averaging [22].Our works provides an evidence for the coherent coupling and classical back-action of the short-lived spin modes with light, as opposed to previous observations of just broadband spin noise [36].The CIFAR technique also paves the way for probing and engineering the optical coupling of higher order spin modes to light, a source of inefficiencies and unwanted noise in quantum limited measurements.By preparing the optical field in a suitable spatial mode, the multimode capabilities of the spin-light platform can be utilized.

Fig. 1 . 2 :
Fig. 1.Experimental setup.A strong linearly polarized LO (red) is mode matched to a weaker, phase modulated beam from an EOM (blue) on a PBS.Part of the light is sent to a polarization sensitive detection setup, which is used to stabilize the relative phase between the two beams.The atomic input is sent through an optically polarized room temperature alkali spin ensemble (orange circles), situated in a homogeneous magnetic field  (purple).The collective total spin vector (orange) modulates the input light polarization due to the Faraday interaction, while the light also drives the spin.The field at the output of the ensemble is detected in a polarization self-homodyning setup.EOM: Electro-optic modulator.LO: Local Oscillator.PBS: Polarizing beamsplitter./2: Half wave plate.

Fig. 2 .
Fig. 2. CIFAR as a function of modulation amplitude.CIFAR response amplitude (top) and phase (bottom) for different electrical EOM drive voltage .The average of 3 scans (dots) is presented along with their statistical 1 uncertainty error bars (vertical bars).The solid lines are the model fits to the individual curves.The grey line in the top panel is the measured response without any modulation at the input.Inset: fitted readout rate and error bars as function of the drive voltage.For a discussion of the error bars, see the main text.

Fig. 4 .
Fig. 4. CIFAR scans for different  S / / /      .We vary the readout rate by changing the temperature of the cell from ∼34 °C to ∼59 °C.Inset: The location of the minimum of the CIFAR response in units of Δ RF / S as a function of the normalized readout rate Γ S / S .Solid line: line with slope 1.

Fig. 5 .
Fig. 5. CIFAR signal different tensor coupling parameters  S .The overall response of the spin oscillator to light depends on  S , here controlled by the angle  between the input linear polarization LO and the direction of the magnetic bias field .The CIFAR signals for input modulation with  = 45°(top panel, logarithmic scale) and with  = 90°(bottom panel, linear scale), and  = {0°, 45°, 90°} are shown in blue, orange and green, respectively.

Fig. 6 .
Fig.6.Coherent interference between the responses of the narrow and broadband spin modes.CIFAR response amplitude (top row) and phase (bottom row) as a function frequency detuning for three different modulation phases,  ∈ {−45°, 0°, 45°}.The data was taken under the same experimental conditions as the 62 mV drive trace (orange curve) in Figure2.For  = 0°we plot the fit result evaluated with the broadband readout rate set to Γ S,BB = 0 (solid light orange curves, top and bottom panels) and narrowband readout rate Γ S = 0 (dashed light orange curves).
The electric field of the laser light after the PBS is  LO () ê +  EOM () ê ∼ | LO |ê  + | EOM |e   (1 +  sin  RF ) ê , to the first order in .Assuming | EOM | | LO |, the equivalent input Stokes parameters [42] can be written, to the first order in | EOM |, as