Abstract
Quantum technologies use entanglement to outperform classical technologies, and often employ strong cooling and isolation to protect entangled entities from decoherence by random interactions. Here we show that the opposite strategyâpromoting random interactionsâcan help generate and preserve entanglement. We use optical quantum non-demolition measurement to produce entanglement in a hot alkali vapor, in a regime dominated by random spin-exchange collisions. We use Bayesian statistics and spin-squeezing inequalities to show that at least 1.52(4)âÃâ1013 of the 5.32(12)âÃâ1013 participating atoms enter into singlet-type entangled states, which persist for tens of spin-thermalization times and span thousands of times the nearest-neighbor distance. The results show that high temperatures and strong random interactions need not destroy many-body quantum coherence, that collective measurement can produce very complex entangled states, and that the hot, strongly-interacting media now in use for extreme atomic sensing are well suited for sensing beyond the standard quantum limit.
Similar content being viewed by others
Introduction
Entanglement is an essential resource in quantum computation, simulation, and sensing1, and is also believed to underlie important many-body phenomena such as high-Tc superconductivity2. In many quantum technology implementations, strong cooling and precise controls are required to prevent entropyâwhether from the environment or from noise in classical parametersâfrom destroying quantum coherence. Quantum sensing3 is often pursued using low-entropy methods, for example, with cold atoms in optical lattices4. There are, nonetheless, important sensing technologies that operate in a high-entropy environment, and indeed that employ thermalization to boost coherence and thus sensor performance. Notably, vapor-phase spin-exchange-relaxation-free (SERF) techniques5 are used for magnetometry6,7, rotation sensing8, and searches for physics beyond the standard model9, and give unprecedented sensitivity10. In the SERF regime, strong, frequent, and randomly-timed spin-exchange (SE) collisions dominate the spin dynamics, to produce local spin thermalization. In doing so, these same processes also decouple the spin degrees of freedom from the bath of centre-of-mass degrees of freedom, which increases the spin coherence time5. Whether entanglement can be generated, survive, and be observed in such a high entropy environment is a challenging open question11.
Here, we study the nature of spin entanglement in this hot, strongly-interacting atomic medium, using techniques of direct relevance to extreme sensing. We apply optical quantum non-demolition (QND) measurement12,13âa proven technique for both generation and detection of non-classical states in atomic mediaâto a SERF-regime vapor. We start with a thermalized spin state to guarantee the zero mean of the total spin variable and use a [1, 1, 1] direction magnetic field (see Fig. 1a) to achieve QND measurements on three components of the total spin variable. We track the evolution of the net spin using the Bayesian method of Kalman filtering14, and use spin squeezing inequalities15,16 to quantify entanglement from the observed statistics. We observe that the QND measurement generates a macroscopic singlet state17âa squeezed state containing a macroscopic number of singlet-type entanglement bonds. This shows that QND methods can generate entanglement in hot atomic systems even when the atomic spin dynamics include strong local interactions. The spin squeezing and thus the entanglement persist far longer than the spin-thermalization time of the vapor; any given entanglement bond is passed many times from atom to atom before decohering. We also observe a sensitivity to gradient fields that indicates the typical entanglement bond length is thousands of times the nearest-neighbor distance. This is experimental evidence of long-range singlet-type entanglement bonds. These experimental observations complement recent predictions of coherent inter-species quantum state transfer by spin collision physics18,19.
a Experimental setup. A linearly polarized probe beam, red detuned by 44 GHz from the 87Rb D1 line, passes through a glass cell containing a hot 87Rb vapor and 100 torr of N2 buffer gas, which is housed in a low-noise magnetic enclosure. The transmitted light is detected with a shot-noise-limited polarimeter (a Wollaston prism plus differential detector), which indicates \({\cal{F}}_{z}\), the projection of the collective spin \({\cal{F}}\) on the probe direction, plus optical shot noise. A static magnetic field along the [1, 1, 1] direction causes the spin components to precess as \({\cal{F}}_z \to {\cal{F}}_x \to {\cal{F}}_y\) every one-third of a Larmor cycle. In this way the polarimeter record contains information about all three components17. b Representative sample of the Stokes parameter \(S_y^{\left( {out} \right)}\left( t \right)\) showing raw data (blue dots), optimal estimate for the atomic spin \(g{\cal{F}}_z\left( t \right)S_x\) (red line), and ±4Ï confidence interval (pale red region), as computed by Kalman filter (KF). Signal clearly shows atomic spin coherence over ms time-scales. c Expanded view of early signal, colors as in (b), bars show ±4Ï confidence regions of \(g{\cal{F}}_z\left( t \right)S_x\) for a thermal spin state (TSS) and the standard quantum limit (SQL) for a spin-polarized state. The KF acquires a sub-SQL estimate for \({\cal{F}}_z\) in 20âμs, far less than the coherence time. Probe powerâ=â2âmW, Larmor frequencyâ=â1.3âkHz, cell temperatureâ=â463âK.
Results
Material system
We work with a vapor of 87Rb contained in a glass cell with buffer gas to slow diffusion, and housed in magnetic shielding and field coils to control the magnetic environment, see Fig. 1a. The density is maintained at nRbâ=â3.6âÃâ1014 atoms/cm3, and the magnetic field, applied along the [1, 1, 1] direction, is used to control the Larmor precession frequency ÏL/2Ï. At this density, the spin-exchange collision rate is 325âÃâ103âsâ1. For ÏL below about 2ÏâÃâ5âkHz, the vapor enters the SERF regime, characterized by a large increase in spin coherence time.
Spin thermalization
The spin dynamics of such dense alkali vapors5 is characterized by a competition of several local spin interactions, diffusion, and interaction with external fields, buffer gases, and wall surfaces. While the full complexity of this scenario has not yet been incorporated in a quantum statistical model, in the SERF regime an important simplification allows us to describe the state dynamics in sufficient detail for entanglement detection, as we now show.
If j(l) and i(l) are the lth atomâs electron and nuclear spins, respectively, the spin dynamics, including sudden collisions, can be described by the time-dependent Hamiltonian
where the terms describe the hyperfine interaction, SE collisions, spin-destruction (SD) collisions, and Zeeman interaction, respectively. Ahf is the hyperfine (HF) splitting and \(t_n^{\left( {l,l^{\prime}} \right)}\) is the (random) time of the n-th SE collision between atoms l and \(l^{\prime}\), which causes mutual precession of j(l) and \({\mathbf{j}}^{\left( {l^{\prime}} \right)}\) by the (random) angle θn. We indicate with RSE the rate at which such collisions move angular momentum between atoms. Similarly, the third term describes rotations about the random direction dm by random angle \({\it{\uppsi }}_m\), and causes spin depolarization at a rate RSD. γeâ=â2ÏâÃâ28âGHzâTâ1 is the electron spin gyromagnetic ratio. We neglect the much smaller i·B coupling. We note that short-range effects of the magnetic dipoleâdipole interaction (MDDI) are already included in RSE and RSD, and that long-range MDDI effects are negligible in an unpolarized ensemble, as considered here.
The SERF regime is defined by the hierarchy Ahfââ«âRSEââ«âγe|B|, RSD. Our experiment is in this regime, as we have Ahfâââ109âsâ1, RSEâââ105âsâ1, γe|B|âââ104âsâ1 and RSDâââ102âsâ1. The hierarchy implies the following dynamics: on short times, the combined action of the HF and SE terms rapidly thermalizes the spin state, i.e., generates the maximum entropy consistent with the ensemble total angular momentum F, which is conserved by these interactions (see âMethodsâ, âSpin thermalizationâ section). We indicate this F-parametrized max-entropy state by \(\rho _{\mathbf{F}}^{\left( {{\mathrm{th}}} \right)}\). We note that entanglement can survive the thermalization process; for example, \(\rho _{{\mathbf{F}} = 0}^{\left( {{\mathrm{th}}} \right)}\) is a singlet and thus necessarily describes entangled atoms. On longer time-scales, F experiences precession about B due to the Zeeman term and diffusive relaxation due to the depolarization term.
Non-destructive measurement
We perform a continuous non-destructive readout of the spin polarization using Faraday rotation of off-resonance light. On passing through the cell the optical polarization experiences rotation by an angle \({g{\mathbf{\cal{F}}}_z(t) \ll \pi}\), where z is the propagation axis of the probe, g is a light-atom coupling constant and \({{{\mathbf{\cal{F}}}}}\)ââ¡âFaâââFb, where Fα is the collective spin orientation from atoms in hyperfine state αâââ{1, 2} (see âMethodsâ, âObserved spin signalâ section).
For thermalized spin states ã\({{{\mathbf{\cal{F}}}}}\) ã\(\propto \left\langle {\mathbf{F}} \right\rangle\), so that the observed polarization rotation gives a view into the full spin dynamics. The optical rotation is detected by a balanced polarimeter (BP), which gives a signal proportional to the Stokes parameter
where Sx is the Stokes component along which the input beam is polarized20. \(S_y^{\left( {{\mathrm{in}}} \right)}\left( t \right)\) is a zero-mean Gaussian process, whose variance is dictated by photon shot-noise and is characterized by a power-spectral analysis of the BP signal21.
Spin dynamics and spin tracking
The evolution of \({\boldsymbol{{\cal{F}}}}\left( t \right)\) is described by the Langevin equation (see âMethodsâ, âSpin dynamicsâ section)
where γâ=âγe/q is the SERF-regime gyromagnetic ratio, i.e., that of a bare electron reduced by the nuclear slowing-down factor5, which takes the value qâ=â6 in the SERF regime22. Î is the net relaxation rate including diffusion, spin-destruction collisions, and probe-induced decoherence, Q is the equilibrium variance (see below) and dWh, hâââ{x, y, z} are independent temporal Wiener increments.
Based on Eqs. (3) and (2), we employ the Bayesian estimation technique of Kalman filtering (KF)14 to recover \({\mathbf{\cal{F}}}\left( t \right)\), which is shown as \(g{\mathbf{\cal{F}}}_z\left( t \right)S_x\) to facilitate comparison against the measured \(S_y^{\left( {{\mathrm{out}}} \right)}\left( t \right)\) in Fig. 1b). The KF (see âMethodsâ, âKalman filterâ section) gives both a best estimate and a covariance matrix \({{\Gamma }}_{\boldsymbol{{\cal{F}}}}\left( t \right)\) for the components of \({\boldsymbol{{\cal{F}}}}\left( t \right)\), which gives an upper bound on the variances of the post-measurement state. Figure 1c shows that the \({\mathbf{\cal{F}}}_z\) component of \({{\Gamma }}_{\boldsymbol{{\cal{F}}}}\left( t \right)\) is suppressed rapidly, to reach a steady state value which is below the SQL. The other components are similarly reduced in variance by the measurement, and the total variance \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2 \equiv {\mathrm{Tr}}\left[ {{{\Gamma }}_{\boldsymbol{{\cal{F}}}}} \right]\) can be compared against spin squeezing inequalities15,16 to detect and quantify entanglement: Defining the spin-squeezing parameter \(\xi ^2 \equiv \left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2/{\mathrm{SQL}}\), where SQLââ¡âNA13/8 is the standard quantum limit, ξ2â<â1 detects entanglement, indicating a macroscopic singlet state17. The minimum number of entangled atoms15 is NA(1âââξ2)13/16 (see âMethodsâ, âEntanglement witnessâ section).
Experimental results
The cell temperature was stabilized at 463âK to give an alkali number density of nRbâ=â3.55(6)âÃâ1014 atomsâcmâ3, calibrated as described in âMethodsâ, âDensity calibrationâ section, and thus NAâ=â5.32(12)âÃâ1013 atoms within the 3âcmâÃâ0.0503(8)âcm2 effective volume of the beam. At this density, the SE collision rate is RSEâââ325âÃâ103âsâ1. By varying B we can observe the transition to the SERF regime, and the consequent development of squeezing. Figure 2a shows spin-noise spectra (SNS)21, i.e., the power spectra of detected signal from BP, for different values of B, from which we determine the resonance frequency ÏLâ=âγB, relaxation rate Î and the number density. Using these as parameters in the KF (see âMethodsâ, âKalman filterâ section), we obtain \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2\) as shown in Fig. 2b, including a transition to squeezed/entangled states as the system enters the SERF regime.
a Spin noise spectra with atomic spin signal driven by thermal fluctuations and precessing at the Larmor frequency (νLâ=âÏL/2Ï) rising above shot noise of the Faraday rotation probe. Different spectra correspond to different bias field strengths. Black lines are single Lorentz fits for the spectra. Comparing the red and purple curves, we see a roughly 100-fold improvement in signal to noise ratio (SNR) due to suppression of SE relaxation and consequent line narrowing, which indicated a stronger quantum non-demolition (QND) interaction. b Spin variance versus Larmor frequency. Black solid-line shows the standard quantum limit of total spin (SQLâ=âNA13/8). Dashed horizontal lines indicate standard quantum limit (SQL) ±1Ï statistical uncertainty. Round symbols show \(\left| {{\mathrm{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2\) measured with 0.5âmW probe light, corresponding to the spectra in (a). Diamonds and squares show \(\left| {{\mathrm{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2\) measured with 1âmW and 2âmW probe light respectively. All error bars show ±1Ï uncertainty due to uncertainty in atomic number, including uncertainties in atomic density and effective volume (see âMethodsâ, âDensity calibrationâ section).
At a Larmor frequency of 1.3âkHz, we observe ξ2â=â0.650(2) or 1.88(1)âdB of spin squeezing at optimal probe power 2âmW (see âMethodsâ, âKalman filterâ section), which implies that at least 1.52(4)âÃâ1013 of the 5.32(12)âÃâ1013 participating atoms have become entangled as a result of the measurement. This greatly exceeds the previous entanglement records: 5âÃâ105 cold atoms in singlet states using a similar QND strategy17 and a Dicke state involving 2âÃâ1011 impurities in a solid, made by storing a single photon in a multi-component atomic ensemble23. This is also the largest number of atoms yet involved in a squeezed state; see Bao et al. for a recent record for polarized spin-squeezed states24. We use this power and field condition for the experiments described below, and note that the spin-relaxation time greatly exceeds the spin-thermalization time. In this condition, the entanglement bonds are rapidly distributed amongst the atoms by SE collisions without being lost.
We now study the spatial distribution of the induced entanglement. As concerns the observable \({\boldsymbol{{\cal{F}}}}\), the relevant dynamical processes, including precession, decoherence, and probing, are permutationally-invariant: Eqs. (3) and (2) are unchanged by any permutation of the atomic states. This suggests that any two atoms should be equally likely to become entangled, and entanglement bonds should be generated for atoms separated by Îzâââ[0, L], where Lâ=â3 is the length of the cell. Indeed, such permutational invariance is central to proposals25,26 that use QND measurement to interrogate and manipulate many-body systems. There are other possibilities, however, such as optical pumping into entangled sub-radiant states27, that could produce localized singlets.
We test for long-range singlet-type entanglement by applying a weak gradient \(B^{\prime} \equiv d\left| B \right|/dz\) during the cw probing process. A magnetic field gradient, if present, causes differential Larmor precession that converts low-noise singlets into high-noise triplets, providing evidence of long-range entanglement. For example, singlets with separation Îz will convert into triplets and back at angular frequency28 \({\mathrm{\Omega }} = \gamma B^{\prime}{{\Delta }}z\). The range δÎz of separations then induces a range \(\delta {\mathrm{\Omega }} = \gamma B^{\prime}\delta {{\Delta }}z\) of conversion frequencies, which describes a relaxation rate. In Fig. 3 we show the KF-estimated \(\left| {{{\Delta }}{\mathbf{\cal{F}}}_z} \right|^2\) as a function of \(B^{\prime}\) and of time since the last data point, which clearly shows faster relaxation toward a thermal spin state with increasing \(B^{\prime}\). The observed additional relaxation for \(B^{\prime}\)â=â57.2ânTâmmâ1 (relative to \(B^{\prime}\)â=â0) is δΩâ=â1.54âÃâ103âsâ1, found by an exponential fit. For Îz on the order of a wavelength, as would describe sub-radiant states, we would expect δΩâ~â1âsâ1 at this gradient, which clearly disagrees with observations. The observed r.m.s. separation δÎz is about one millimeter, which is thousands of times the typical nearest-neighbor distance \(n_{{\mathrm{Rb}}}^{ - 1/3} \approx 0.14\)âμm.
Points show the KF-obtained variances \(\left\langle {{\mathrm{\Delta }}{\mathbf{\cal{F}}}_z^2} \right\rangle\) of the rotating-frame spin component \({\mathbf{\cal{F}}}_z\) as a function of delay since \({\mathbf{\cal{F}}}_z\) was last aligned along the laboratory z axis, and thus subject to measurement by Faraday rotation. Error bars show ±4Ï uncertainty, originating in the uncertainty of atom number. As seen in these data, increased gradient \(B^{\prime}\) causes a faster relaxation toward the thermal spin state (TSS) value, as is expected for a gas of singlets with a range of separations (along the z axis) in the mm range. The blue and gray dashed lines are the standard quantum limit (SQL) (2.88(7)âÃâ1013 spins2) and TSS (7.99(18)âÃâ1013âspins2) noise levels, respectively. Probe powerâ=â2âmW.
Discussion
Our observation of complex, long-lived, spatially-extended entanglement in SERF-regime vapors has a number of implications. First, it is a concrete and experimentally tractable example of a system in which entanglement is not only compatible with, but, in fact, stabilized by entropy-generating mechanismsâin this case strong, randomly-timed spin-exchange collisions. It is particularly intriguing that the observed macroscopic singlet state shares several traits with a spin liquid state2, which is conjectured to underlie high-temperature superconductivity, a prime example of quantum coherence surviving in an entropic environment. Second, the results show that optical quantum non-demolition measurement can efficiently produce complex entangled states with long-range entanglement. This confirms a critical assumption of QND-based proposals25,26 for QND-assisted quantum simulation of exotic antiferromagnetic phases. Third, the results show that SERF media are compatible with both spin squeezing and QND techniques, opening the way to quantum enhancement of what is currently the most sensitive approach to low-frequency magnetometry and other extreme sensing tasks.
Methods
Density calibration
In the SERF regime, and in the low spin polarization limit, decoherence introduced by SE collisions between alkali atoms is quantified by 5,29,30
where for 87Rb atomic samples the nuclear spin Iâ=â3/2, and ÏLâ=âγe|B|/q. In Eq. (4) the spin-exchange collision rate \(R_{{\mathrm{SE}}} = \sigma _{{\mathrm{SE}}}n_{{\mathrm{Rb}}}\overline V\) is proportional to the alkali density nRb with proportionality dictated by the SE collision cross-section ÏSE and the relative thermal velocity between two colliding 87Rb atoms \(\overline V\). Using the reported value31 of ÏSEâ=â1.9âÃâ10â14âcm2 and \(\overline V = 4.75 \times 10^4\,{\mathrm{cm}}\;{\mathrm{s}}^{ - 1}\), which is computed for 87Rb atoms at a temperature of 463âK, we then calibrate the alkali density by fitting the measured linewidth Îν as a function of ÏL. The model uses Îνâ=âÎν0â+âÎνSE, where ÎνSE is given by Eq. (4), and Îν0 describes density-independent broadening due to power broadening and transit effects. nRb and Îν0 are free parameters found by fitting, with results shown in Fig. 4.
Graph shows the full-width at half maximum linewidth Îν as a function of ÏL, across the transition into the spin-exchange relaxation free (SERF) regime. Blue dots and error bars show the mean ±1Ï standard error of the mean, obtained from the fitting of 20 spectra. Red line shows Eq. (4) fit to the data with density nRb as a free parameter. Probe powerâ=â1âmW, Tâ=â463âK.
Observed spin signal
For a collection of atoms, we define the collective total atomic spin \({\mathbf{F}} \equiv \mathop {\sum}\nolimits_l {{\mathbf{f}}^{\left( l \right)}}\), where f(l) is the total spin of the lth atom. We identify the contributions of the two hyperfine ground states Faâ=â1 and Fbâ=â2, defined as \({\mathbf{F}}_\alpha \equiv \mathop {\sum}\nolimits_l {{\mathbf{f}}_\alpha ^{\left( l \right)}}\), where \({\mathbf{f}}_\alpha ^{\left( l \right)}\) describes the contribution of atoms in state Fα, such that \({\mathbf{f}}^{\left( l \right)} = {\mathbf{f}}_a^{\left( l \right)} + {\mathbf{f}}_b^{\left( l \right)}\).
The Faraday rotation signal arises from an off-resonance coupling of the probe light to the collective atomic spin. To lowest order in F, as appropriate to the regime of the experiment, the polarization signal Sy is related to the collective spin variables Fa,z, Fb,z through the input-output relation20,32,33,34
where \(S_\alpha \equiv ( {E_ + ^{( -)},E_ - ^{( - )}} )\sigma _\alpha ( {E_ + ^{( + )},E_ - ^{( + )}} )^T/2\) are Stokes operators, Ïα, αâââ{x, y, z} are the Pauli matrices and \(E_\beta ^{\left( \pm \right)}\) is the positive-frequency (negative-frequency) part of the quantized electromagnetic field with polarization βâ=â± for sigma-plus (sigma-minus) polarized light. The factor (gaFa,zâââgbFb,z)ââ¡âÎFR plays the role of a Faraday rotation angle, which in this small-angle regime can be seen to cause a displacement of Sy(t) from its input value. It should be noted that ÎFR is operator-valued, enabling entanglement of the spin and optical polarizations, and that the hyperfine ground states Fa, Fb contribute differentially to it.
The coupling constants are14,33,35
where reâ=â2.82âÃâ10â13âcm is the classical electron radius, foscâ=â0.34 is the oscillator strength of the D1 transition in Rb, c is the speed of light, and νâââνα is the optical detuning of the probe-light. \(\Upsilon = 2.4\)âGHz is the pressure-broadened full-width at half-maximum (FWHM) linewidth of the D1 optical transition for our experimental conditions of 100âTorr of N2 buffer gas. For a far-detuned probe beam, such that |νâââ(νaâ+âνb)/2|ââ«â|νaâââνb| (as in this experiment) one can approximate gââ¡âgaâââgb, such that
Spin thermalization
In a local region containing a mean number of atoms NA, the SE and HF mechanisms will rapidly produce a thermal state Ï. We note that this process conserves F, and thus also conserves the statistical distribution of F, including possible correlations with other regions. Ï is then the maximum-entropy state consistent with a given distribution of F. Partitioning arguments then show that, for weakly polarized states such as those used in this experiment, the mean hyperfine populations are \(\left\langle {N_a} \right\rangle /N_{\mathrm{A}} = 3/8\) and \(\left\langle {N_b} \right\rangle /N_{\mathrm{A}} = 5/8\), and the polarisations are \(\left\langle {{\mathbf{F}}_a} \right\rangle = \left\langle {\mathbf{F}} \right\rangle /6\), \(\left\langle {\mathbf{F}} \right\rangle _b = \left\langle {\mathbf{F}} \right\rangle 5/6\), from which the FR signal is \(\left\langle {{\mathrm{\Theta }}_{{\mathrm{FR}}}} \right\rangle = g\left\langle {{\mathbf{\cal{F}}}_z} \right\rangle = - g\left\langle F \right\rangle _z2/3\). The same relations must hold for spin observables that sum F over larger regions, including the region of the beam, which determines which atoms contribute to the observed signal.
Entanglement witness
We can construct a witness for singlet-type entanglement16 as follows: we define the total variance
Separable states of NA atoms will obey a limit \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2 \ge N_{\mathrm{A}}{\cal{C}}\), where \({\cal{C}}\) is a constant, meaning that \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2 < N_{\mathrm{A}}{\cal{C}}\) witnesses entanglement. To find \({\cal{C}}\), we note that a product state of Na atoms in state Fa and Nb atoms in state Fb has \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2 \ge \mathop {\sum}\nolimits_\alpha {N_\alpha F_\alpha }\). Separable states are mixtures of product states. For such states, due to the concavity of the variance, \(\left| {{{\Delta }}{\mathbf{\cal{F}}}} \right|^2 \ge \mathop {\sum}\nolimits_\alpha {\left\langle {N_\alpha } \right\rangle F_\alpha }\) holds16. In light of the 3:5 ratio resulting from spin thermalization, this gives
or \({\cal{C}} = 13/8\). Therefore, the standard quantum limit (SQL) is NA13/8. We define the degree of squeezing \(\xi ^2 \equiv \left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2/\left( {\mathop {\sum}\nolimits_\alpha {\left\langle {N_\alpha } \right\rangle F_\alpha } } \right)\). Meanwhile the âthermal spin state (TSS),â i.e. the fully-mixed state, has \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2 = N_{\mathrm{A}}9/2\).
Our condition provides also a quantitative measure of the number of entangled atoms. We consider a pure entangled quantum state of the form
where \(| {{\it{\uppsi }}^{\left( l \right)}} \rangle\) are single particle states. Here, Np particles are in a product state, while Neâ=âNAâââNp particles are in an entangled state denoted by \(| {\Phi _{\mathrm{e}}} \rangle\). For the collective variances of \(| \kappa \rangle\) we can write that
for hâ=âx, y, z. Let us try to find a lower bound on (11).
Let us assume that all atoms are in state Fα. Then, we know that \(\left( {{{\Delta }}{\mathbf{\cal{F}}}_h} \right)_{{\it{\uppsi }}^{\left( l \right)}}^2 \ge F_\alpha\) while \(\left( {{{\Delta }}{\mathbf{\cal{F}}}_h} \right)_{\Phi _{\mathrm{e}}}^2\) can even be zero, if the entangled state \(| {\Phi _{\mathrm{e}}} \rangle\) is a perfect singlet. Hence,
Based on these, the number of entangled atoms in this case is bounded from below as Neââ¥â(1âââξ2)NA, where \(\xi ^2 = \left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2/{\mathrm{SQL}}\) and the standard quantum limit SQL in this case is FαNA.
Let us now consider the case when some atoms have F1 others have F2 In particular, let us consider a state of the type (10) such that Nα particles have spin Fα with αâ=â1, 2 such that F1ââ¤âF2. Then, for such a pure state,
holds, where
Note that the bound in Eq. (13) is sharp, since it can be saturated by a quantum state of the type (10). In order to minimize the left-hand side of Eq. (13), the particles corresponding to the product part must have as many spins in F1 as possible, since this way we can obtain a small total variance. In particular, if Npââ¥âN1 then all atoms in the product part must have an F1 spin, otherwise at least N1 atoms of the Np atoms.
It is instructive to rewrite Eq. (10) with a piece-wise linear bound as
The bound in Eq. (15) is plotted in Fig. 5a.
a Lower bound on the sum of the three variances given in Eq. (13) as a function of Np for a quantum state of the type given in Eq. (10). We set N1/NAâ=â3/8, F1â=â1, F2â=â2, corresponding to the experiment. If we had only particles with the same spin, it would just be a straight line. b The lower bound on the number of entangled spins, given in Eq. (16), as a function of the spin-squeezing parameter ξ2.
So far, we have been discussing a bound for a pure state of the form (10). The results can be extended to a mixture of such states straightforwardly, since the bound in Eq. (13) is convex in (N1, Np). Then, in our formulas N1, must be replaced by \(\left\langle {N_1} \right\rangle\). We also have to define the number of entangled particles Ne for the case of a mixed state. A mixed state has Ne entangled particles, if it cannot be constructed as a mixture of pure states, which all have fewer than Ne entangled particles15.
We know that in our experiments F1â=â1, F2â=â2, and \(\left\langle {N_1} \right\rangle = 3/8N_{\mathrm{A}}\). From these, we obtain the minimum number of entangled atoms as
The bound in Eq. (16) is plotted in Fig. 5b. Here, again \(\xi ^2 = \left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2/{\mathrm{SQL}}\) and the standard quantum limit SQL in this case is NAâÃâ13/8. For our experiment, the ξ2âââ¥â3/13 case is relevant.
The macroscopic singlet state gives a metrological advantage in estimating gradient fields28,36 and in detecting displacement of the spin state, e.g. by optical pumping14,37.
Balanced polarimeter signal
The photocurrent I(t) of the balanced polarimeter shown in Fig. 1a is
where the detectorâs responsivity is \(\Re = q_{\mathrm{e}}\eta /E_{{\mathrm{ph}}}\) in terms of the detector quantum efficiency η, charge of the electron qe, and photon energy Eph. To account for its spatial structure in Eq. (17) the integral is carried over the area of the probe. From Eq. (2) and Eq. (17) one obtains the differential photocurrent increment
where \(g' = gq_{\mathrm{e}}\), the stochastic increment dwsn(t), due to photon shot-noise, is given by \(dw_{{\mathrm{sn}}}\left( t \right) = \sqrt {\eta q_{\mathrm{e}}^2\dot N} dW\) with \(\dot N\) being the photon-flux and \(dW\sim {\cal{N}}\left( {0,dt} \right)\) representing a differential Wiener increment. In our experiments the photocurrent I(t) is sampled at a rate Îâ1â=â200âk Samples/s. To formulate the discrete-time version of Eq. (17) we consider the sampling process as a short-term average of the continuous-time measurement. The photocurrent I(tk) recorded at tkâ=âkÎ, with k being an integer, can then be expressed as
where the Langevin noise ξD(tk) obeys \(E\left[ {\xi _D\left( t \right)\xi _D\left( {t^{\prime}} \right)} \right] = \delta \left( {t - t^{\prime}} \right)\eta q_{\mathrm{e}}^2\dot N/{{\Delta }}\), with Îâ1 quantifying the effective noise-bandwidth of each observation.
Spin dynamics
We model the dynamics of the average bulk spin of our hot atomic vapor in the SERF regime,5,29,38 and in the presence of a magnetic field B in the [1,1,1] direction, i.e. \({\mathbf{B}} = {\mathrm{B}}\left( {{\hat{\mathrm{x}}} + {\hat{\mathrm{y}}} + {\hat{\mathrm{z}}}} \right)/\sqrt 3\), as
where the matrix A includes dynamics due to Larmor precession and spin relaxation. It can be expressed as Aijâ=ââγBhεhijâ+âÎij, where h, i, jâ=âx, y, z. The relaxation matrix Î has eigenvalues \(T_1^{ - 1}\) and \(T_2^{ - 1} = T_1^{ - 1} + T_{{\mathrm{SE}}}^{ - 1}\) for spin components parallel and transverse to B, respectively. We note that in the SERF regime the decoherence introduced by SE collisions between alkali atoms is quantified by Eq. (4)5,29.
To account for fluctuations due to spin noise in Eq. (20) we add a stochastic term \(\sqrt {\boldsymbol{\sigma} } d{\mathbf{W}}\) where dWh, hâââ{x, y, z) are independent Wiener increments. Thus the statistical model for spin dynamics reads
where the strength of the noise source Ï, the matrix A, and the covariance matrix in statistical equilibrium \(Q = {\mathrm{E}}\left[ {{\mathbf{\cal{F}}}\left( t \right){\mathbf{\cal{F}}}\left( t \right)^{\rm{T}}} \right]\) are related by the fluctuation-dissipation theorem
from which we obtain Ïâ=â2ÎQ.
Kalman filter
Kalman filtering is a signal recovery method that provides continuously-updated estimates of all physical variables of a stochastic model, along with uncertainties for those estimates. For linear dynamical systems with gaussian noise inputs, e.g. the spin dynamics of Eq. 3 with the readout of Eq. 5, the Kalman filter estimates are optimal in a least-squares sense. The KF estimates, e.g. those shown in Fig. 1b and c, indicate our evolving uncertainty about the values of the physical quantities, e.g. \({\mathbf{\cal{F}}}_z\). As such, they provide an upper bound on the intrinsic uncertainty of these same quantities due to, e.g. quantum noise. As information accumulates, the uncertainty bounds on \({\mathbf{\cal{F}}}_x\), \({\mathbf{\cal{F}}}_y\) and \({\mathbf{\cal{F}}}_z\) contract toward zero, implying the production of squeezing and entanglement. This is measurement-induced, rather than dynamically-generated entanglement. The measured signal, i.e. the optical polarization rotation, indicates a joint atomic observable: the sum of the spin projections of many atoms. For an unpolarized state such as we use here the physical back-actionâwhich consists of small random rotations about the \({\mathbf{\cal{F}}}_z\) axis induced by quantum fluctuations in the ellipticity of the probeâhas a negligible effect.
We construct the estimator \(\widetilde {\mathbf{\cal{F}}}_t\) of the macroscopic spin vector using the continuous-discrete version of Kalman filtering14. This framework relies on a two-step procedure to construct the estimate \({\tilde{\mathbf{x}}}_t\), and its error covariance matrix \({\mathbf{\Sigma }} = {\mathrm{E}}\left[ {\left( {{\mathbf{x}}_t - {\tilde{\mathbf{x}}}_t} \right)\left( {{\mathbf{x}}_t - {\tilde{\mathbf{x}}}_t} \right)^T} \right]\), of the state xt of a continuous-time linear-Gaussian process, in our case \({\boldsymbol{{\cal{F}}}}\left( t \right)\), that is observed at discrete-time intervals Îâ=âtkââââtkâ1. Measurement outcomes are described by the observations vector zk, in our case the scalar Ik, which is assumed to be linearly related to xt via the coupling matrix Hk and to experience independent stochastic Gaussian noise as described previously14.
In the first step of the Kalman filtering framework, also called the prediction step, the values at tâ=âtk, \(\widetilde {\boldsymbol{{\cal{F}}}}_{k|k - 1}\) and \({\mathbf{\Sigma }}_{k|k - 1}\), are predicted conditioned on the process dynamics and the previous instance, \(\widetilde {\boldsymbol{{\cal{F}}}}_{k - 1|k - 1}\) and \({\mathbf{\Sigma }}_{k - 1|k - 1}\), as follows:
where
is the state transition matrix describing the evolution of the dynamical model Eq. (20) within the time interval Î, and
is then the effective covariance matrix of the system noise14.
In the second step, or update step, the information gathered through the fresh photocurrent observation Ik is incorporated into the estimate:
where \({\mathbf{H}}_k = \left[ {\eta g\dot N,0,0} \right]\) and the Kalman gain Kk is defined as
with sensor covariance RÎâ=âR/Î dictated by the power-spectral-density, R, of the photocurrent noise, i.e. due to photon shot-noise, and the sampling period, Î. As dicussed in previous work14 the KF is initialised according to a distribution that represents our prior knowledge about the system at time tâ=ât0 and fixes \(\widetilde {\boldsymbol{{\cal{F}}}}_{0|0}\sim {\cal{N}}\left( {{\boldsymbol{\mu }}_0,{\mathbf{\Sigma }}_0} \right)\), where μ0, and Σ0 are the mean value and total variance of the observed data. After initialization KF estimates for the covariance matrix Σk|k undergo a transient and once this transient has decayed they converge to a steady state value Σss.
In Fig. 6 we observe this behavior for the total variance \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2\, \equiv\, {\mathrm{Tr}}\left( {{{\Gamma }}_{\boldsymbol{{\cal{F}}}}} \right)\) as a function of time tâ=âtk, where \({{\Gamma }}_{\boldsymbol{{\cal{F}}}} = {\mathbf{\Sigma }}_{k|k}\). After about 0.8âms, the total variance reaches steady state value which is used to compare with SQL and indicates squeezing degree. Figure 7 shows squeezing degree at different probe power, and presents the optimal probe power we observed is 2âmW.
Validation
To validate the sensor model we employ three validation techniques sensitive to both the statistics of the optical readout and spin noise. First, we analyse the statistics of the sensor output innovation, i.e., the difference between observations Ik (data) and Kalman estimates (\(\tilde y_k = I_k - {\mathbf{H}}_k\widetilde {\boldsymbol{{\cal{F}}}}_{k|k - 1}\)). In Fig. 8, we show the \(\tilde y_k\) histogram with the sensor output estimation error, which is described by zero-mean Gaussian process with variance equal to \({\mathbf{R}}^{{\Delta }} + {\mathbf{H}}_k{\mathbf{\Sigma }}_{k|k - 1}{\mathbf{H}}_k^T\). We find 94% of \(\tilde y_k\) data lie within a two-sided 95% confidence region of the expected Gaussian distribution, thus indicating a very close agreement of the model and observed statistics. We note that while being a standard technique in the validation of Kalman filtering39, this technique for our experimental conditions is more sensitive to photon shot noise than to spin noise. Therefore, to further validate our estimates we also include two other validation techniques, designed to be sensitive to the atomic statistics on a range of time-scales.
Particularly, we perform Monte Carlo simulations based on the model described by Eqs. (2), (3) and (17) and fed with the operating conditions of our experiments and compare the power spectral density (PSD) of the simulated sensor output (Simulation) to the observed PSD of the measurements (Data), as shown in Fig. 9a. The observed agreement between Data and Simulation suggests the validity of the statistics of the spin dynamics model.
a Spin noise spectroscopy (SNS) of Data (blue dots) and Simulation (green dots). The Lorentz fittings of Data (black line) and Simulation (red line) are totally overlapped. The spin distributions (see text) from Data (b) and Simulation (c) are shown as histograms. Error bars indicate plus/minus one standard deviation of histograms of 20 traces. Probe powerâ=â2âmW, νLâ=â1.3âkHz.
Finally, we employ the Kalman filter to identify the evolution of the atomic state variables based on the Simulation. We can then compare the distribution of Kalman spin-estimates from the Data versus that from Simulation. The results are shown in Fig. 9b and c, respectively. The similarity in the statistics of these two spin estimates validates the spin dynamics model. Together with the above validations, it provides a full validation of both the optical and spin parts of the model.
Gradient field tests
A weak gradient magnetic field is applied along the probe (z) direction by coils implemented inside the magnetic shields. In Fig. 10 we plot the three components of \(\left| {{{\Delta }}{\boldsymbol{{\cal{F}}}}} \right|^2\): \(\left| {{{\Delta }}{\mathbf{\cal{F}}}_z} \right|^2 \equiv {{\Gamma }}_{\boldsymbol{{\cal{F}}}}\left( {1,1} \right)\), \(\left| {{{\Delta }}{\mathbf{\cal{F}}}_x} \right|^2 \equiv {{\Gamma }}_{\boldsymbol{{\cal{F}}}}\left( {2,2} \right)\), and \(\left| {{{\Delta }}{\mathbf{\cal{F}}}_y} \right|^2 \equiv {{\Gamma }}_{\boldsymbol{{\cal{F}}}}\left( {3,3} \right)\), as a function of gradient field. Here \({{\Gamma }}_{\boldsymbol{{\cal{F}}}}\left( {i,i} \right) = {\mathbf{\Sigma }}_{ss}\left( {i,i} \right)\). We observe that the variance of each component increases towards the TSS noise level with gradient field. We note that due to the bias field along the [1,1,1] direction, the current (tâ=âtk) sensor reading indicates \({\mathbf{\cal{F}}}_z\) at that time, while \({\mathbf{\cal{F}}}_x\) and \({\mathbf{\cal{F}}}_y\) describe components that were measured 1/3 and 2/3 Larmor cycles earlier, respectively. The combined variance is used to compute \(\left| {{{\Delta }}{\mathbf{\cal{F}}}_z} \right|^2\), as in Fig. 3. We note that the SternâGerlach (SG) effect, in which a gradient causes wave-functions components to separate in accordance with their magnetic quantum numbers, also contributes to the loss of coherence. The SG contribution is negligible, however, due to the weak gradients used here and the rapid randomization of momentum caused by the buffer gas.
Spin variance components increase with gradient field strength. Red dots, yellow squares, blue diamonds represent \(\left\langle {{\mathrm{\Delta }}{\mathbf{\cal{F}}}_z^2} \right\rangle\), \(\left\langle {{\mathrm{\Delta }}{\mathbf{\cal{F}}}_x^2} \right\rangle\), and \(\left\langle {{\mathrm{\Delta }}{\mathbf{\cal{F}}}_y^2} \right\rangle\), respectively. Error bars show ±4Ï uncertainty in the variance, originating in the uncertainty of the atomic number. The dashed line and solid line show thermal spin state (TSS) (7.99(18)Ã1013spins2) and standard quantum limit (SQL) (2.88(7)âÃâ1013spins2) noise levels, respectively for one spin component \({\mathbf{\cal{F}}}_h,h \in \left\{ {x,y,z} \right\}\).
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request. Open-access datasets from this work available at https://doi.org/10.5281/zenodo.3694692.
References
Giovannetti, V., Lloyd, S. & Maccone, L. Quantum-enhanced measurements: beating the standard quantum limit. Science 306, 1330â1336 (2004).
Anderson, P. W. The resonating valence bond state in La2CuO4 and superconductivity. Science 235, 1196â1198 (1987).
Pezzè, L., Smerzi, A., Oberthaler, M. K., Schmied, R. & Treutlein, P. Quantum metrology with nonclassical states of atomic ensembles. Rev. Mod. Phys. 90, 035005 (2018).
Hosten, O., Engelsen, N. J., Krishnakumar, R. & Kasevich, M. A. Measurement noise 100 times lower than the quantum-projection limit using entangled atoms. Nature 529, 505â508 (2016).
Savukov, I. M. & Romalis, M. V. Effects of spin-exchange collisions in a high-density alkali-metal vapor in low magnetic fields. Phys. Rev. A 71, 023405 (2005).
Dang, H. B., Maloof, A. C. & Romalis, M. V. Ultrahigh sensitivity magnetic field and magnetization measurements with an atomic magnetometer. Applied Physics Letters 97, 151110 (2010).
Griffith, W. C., Knappe, S. & Kitching, J. Femtotesla atomic magnetometry in a microfabricated vapor cell. Optics Express 18, 27167â27172 (2010).
Kornack, T. W., Ghosh, R. K. & Romalis, M. V. Nuclear spin gyroscope based on an atomic comagnetometer. Phys. Rev. Lett. 95, 230801 (2005).
Lee, J., Almasi, A. & Romalis, M. Improved limits on spin-mass interactions. Phys. Rev. Lett. 120, 161801 (2018).
Kominis, I., Kornack, T., Allred, J. & Romalis, M. A subfemtotesla multichannel atomic magnetometer. Nature 422, 596â599 (2003).
Kominis, I. K. Sub-shot-noise magnetometry with a correlated spin-relaxation dominated alkali-metal vapor. Phys. Rev. Lett. 100, 073002â (2008).
Grangier, P., Levenson, J. A. & Poizat, J.-P. Quantum non-demolition measurements in optics. Nature 396, 537â542 (1998).
Sewell, R. J., Napolitano, M., Behbood, N., Colangelo, G. & Mitchell, M. W. Certified quantum non-demolition measurement of a macroscopic material system. Nat Photon 7, 517â520 (2013).
Jiménez-MartÃnez, R. et al. Signal tracking beyond the time resolution of an atomic sensor by kalman filtering. Phys Rev Lett 120, 040503 (2018).
Tóth, G. & Mitchell, M. W. Generation of macroscopic singlet states in atomic ensembles. New Journal of Physics 12, 053007 (2010).
Vitagliano, G., Hyllus, P., Egusquiza, I. L. & Tóth, G. Spin squeezing inequalities for arbitrary spin. Phys. Rev. Lett. 107, 240502 (2011).
Behbood, N. et al. Generation of macroscopic singlet states in a cold atomic ensemble. Phys. Rev. Lett. 113, 093601 (2014).
Katz, O., Shaham, R. & Firstenberg, O. Quantum interface for noble-gas spins. Preprint at http://arxiv.org/1905.12532 (2019).
Dellis, A. T., Loulakis, M. & Kominis, I. K. Spin-noise correlations and spin-noise exchange driven by low-field spin-exchange collisions. Phys. Rev. A 90, 032705 (2014).
Colangelo, G. et al. Quantum atomâlight interfaces in the gaussian description for spin-1 systems. New J. Phys 15, 103007 (2013).
Lucivero, V. G., Jiménez-MartÃnez, R., Kong, J. & Mitchell, M. W. Squeezed-light spin noise spectroscopy. Phys. Rev. A 93, 053802 (2016).
Seltzer, S. J. Developments in Alkali-Metal Atomic Magnetometry,. PhD thesis, (Princeton University, 2008).
Zarkeshian, P. et al. Entanglement between more than two hundred macroscopic atomic ensembles in a solid. Nat. Commun. 8, 906 (2017).
Bao, H. et al. Spin squeezing of 1011 atoms by prediction and retrodiction measurements. https://doi.org/10.1038/s41586-020-2243-7 (2018).
Eckert, K., Zawitkowski, L., Sanpera, A., Lewenstein, M. & Polzik, E. S. Quantum polarization spectroscopy of ultracold spinor gases. Phys. Rev. Lett. 98, 100404 (2007).
Hauke, P., Sewell, R. J., Mitchell, M. W. & Lewenstein, M. Quantum control of spin correlations in ultracold lattice gases. Phys. Rev. A 87, 021601 (2013).
Guerin, W., Araújo, M. O. & Kaiser, R. Subradiance in a large cloud of cold atoms. Phys. Rev. Lett. 116, 083601 (2016).
Urizar-Lanz, I., Hyllus, P., Egusquiza, I. L., Mitchell, M. W. & Tóth, G. Macroscopic singlet states for gradient magnetometry. Phys. Rev. A 88, 013626 (2013).
Allred, J. C., Lyman, R. N., Kornack, T. W. & Romalis, M. V. High-sensitivity atomic magnetometer unaffected by spin-exchange relaxation. Phys. Rev. Lett. 89, 130801 (2002).
Happer, W. & Tang, H. Spin-exchange shift and narrowing of magnetic resonance lines in optically pumped alkali vapors. Phys. Rev. Lett. 31, 273â276 (1973).
Happer, W., Jau, Y. Y. & Walker, T. Optically Pumped Atoms. (Wiley, 2010).
Madsen, L. B. & Mølmer, K. Spin squeezing and precision probing with light and samples of atoms in the gaussian description. Phys. Rev. A 70, 052324 (2004).
Geremia, J. M., Stockton, J. K. & Mabuchi, H. Tensor polarizability and dispersive quantum measurement of multilevel atoms. Phys. Rev. A 73, 042112 (2006).
Koschorreck, M. & Mitchell, M. W. Unified description of inhomogeneities, dissipation and transport in quantum light-atom interfaces. J. Phys. B: At. Mol. Opt. Phys. 42, 195502 (2009). (9pp).
Happer, W. Optical pumping. Rev. Mod. Phys. 44, 169â249 (1972).
Apellaniz, I., Urizar-Lanz, In., Zimborás, Z., Hyllus, P. & Tóth, G. Precision bounds for gradient magnetometry with atomic ensembles. Phys. Rev. A 97, 053603 (2018).
Behbood, N. et al. Feedback cooling of an atomic spin ensemble. Phys. Rev. Lett. 111, 103601 (2013).
Ledbetter, M. P., Savukov, I. M., Acosta, V. M., Budker, D. & Romalis, M. V. Spin-exchange-relaxation-free magnetometry with cs vapor. Phys. Rev. A 77, 033408 (2008).
Bar-Shalom, Y., Li, X. & Kirubarajan, T. Estimation with Applications to Tracking and Navigation: Theory Algorithms and Software (Wiley, 2004).
Acknowledgements
We thank Jan Kolodynski for helpful discussions. This project has received funding from the European Unionâs Horizon2020 research and innovation programme under the Marie SkÅodowska-Curie grant agreements QUTEMAG (no. 654339). The work was also supported by ICFOnest + Marie SkÅodowska-Curie Cofund (FP7-PEOPLE-2013-COFUND), the National Natural Science Foundation of China (NSFC) (grant no. 11935012), and ITN ZULF-NMR (766402); the European Research Council (ERC) projects AQUMET (280169), ERIDIAN (713682); European Union projects QUIC (Grant Agreement no. 641122) and FET Innovation Launchpad UVALITH (800901); Quantum Technology Flagship projects MACQSIMAL (820393) and QRANGE (820405); 17FUN03-USOQS, which has received funding from the EMPIR programme co-financed by the Participating States and from the European Unionâs Horizon 2020 research and innovation programme; the Spanish MINECO projects MAQRO (Ref. FIS2015-68039-P), XPLICA (FIS2014-62181-EXP), OCARINA (Grant Ref. PGC2018-097056-B-I00) and Q-CLOCKS (PCI2018-092973), MCPA (FIS2015-67161-P), the Severo Ochoa programme (SEV-2015-0522); Agència de Gestió dâAjuts Universitaris i de Recerca (AGAUR) project (2017-SGR-1354); Fundació Privada Cellex and Generalitat de Catalunya (CERCA program, QuantumCAT), the EU COST Action CA15220 and QuantERA CEBBEC, the Basque Government (Project No. IT986-16), and the National Research, Development and Innovation Office NKFIH (Contract No. K124351, KH129601). We thank the Humboldt Foundation for a Bessel Research Award.
Author information
Authors and Affiliations
Contributions
J.K. and M.W.M. designed the experiment, analyzed the data and wrote the paper. J.K., C.T., and V.G.L. performed the experiment. J.K., R.J.-M., and M.W.M built the Kalman filter model. G.T. and M.W.M. built the entanglement witness. M.W.M. supervised the project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks Denis Sukachev, Iannis Kominis and the other, anonymous, reviewer for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kong, J., Jiménez-MartÃnez, R., Troullinou, C. et al. Measurement-induced, spatially-extended entanglement in a hot, strongly-interacting atomic system. Nat Commun 11, 2415 (2020). https://doi.org/10.1038/s41467-020-15899-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-020-15899-1
This article is cited by
-
The rise of patient avatars in precision oncology
Nature Biotechnology (2024)
-
A polarized 3He system based on metastability-exchange optical pumping
Science China Physics, Mechanics & Astronomy (2024)
-
Entangling motional atoms and an optical loop at ambient condition
npj Quantum Information (2023)
-
Strong coupling of alkali-metal spins to noble-gas spins with an hour-long coherence time
Nature Physics (2022)
-
Bodies and minds, heaps and syllables
Synthese (2021)