# Spin squeezing of high-spin, spatially extended quantum fields

## Abstract

Investigations of spin squeezing in ensembles of quantum particles have been limited primarily to a subspace of spin fluctuations and a single spatial mode in high-spin and spatially extended ensembles. Here, we show that a wider range of spin-squeezing is attainable in ensembles of high-spin atoms, characterized by sub-quantum-limited fluctuations in several independent planes of spin-fluctuation observables. Further, considering the quantum dynamics of an ferromagnetic spinor Bose-Einstein condensate, we demonstrate theoretically that a high degree of spin squeezing is attained in multiple spatial modes of a spatially extended quantum field, and that such squeezing can be extracted from spatially resolved measurements of magnetization and nematicity, i.e. the vector and quadrupole magnetic moments, of the quantum gas. Taking into account several experimental limitations, we predict that the variance of the atomic magnetization and nematicity may be reduced as far as 20 dB below the standard quantum limits.

Spin squeezing describes a form of entanglement in a many-spin system for which the variance in measurements of particular collective spin variables is smaller than the minimum variance achievable in uncorrelated systems [1, 2]. Aside from representing a theoretically tractable form of entanglement in many-body quantum systems, spin squeezing also promises to improve measurement precision in a variety of applications such as chronometry, magnetometry, and atom interferometry. The recent achievement of metrologically beneficial pseudo-spin squeezing in cold atomic gases, on the Zeeman-insensitive hyperfine transitions utilized in atomic clocks, is among the first realizations of this promise [3, 4].

The conventional form of spin squeezing involves the collective vector spin operator, , which is the equal-weight sum of the vector spin operators for all of the ensemble’s constituent particles. Given an average ensemble spin oriented along the axis, the measurement variance of orthogonal spin projections obey an uncertainty relation, . Ensemble states of uncorrelated particles are limited by the standard quantum limited variance, , achieved by partitioning the minimum measurement uncertainties equally between the two vector spin projections.

This form of *vector operator* spin squeezing describes present
experiments not only on pseudo-spin 1/2 atomic gases, i.e. where
each atom is restricted to just two internal
[3, 4] or external [5] states,
but also on gases of higher-spin atoms (e.g. the
hyperfine-state manifold of Cs [6]). However, the
description of higher-spin atoms by their vector spin alone is
incomplete. On the one hand, light-atom interactions that affect
higher-order spin moments may impair efforts to create vector-spin
squeezing in high-spin atomic gases [7, 8]. On
the other hand, the deliberate addressing and measurement of
higher-order spin moments [9] offers benefits to
certain metrological applications, e.g. eliminating directional
dependencies of atomic magnetometers [10].

Here, we consider theoretically the characterization and preparation of spin-squeezed states in spinor Bose gases. We go beyond existing treatments of spin squeezing by evaluating such squeezing in an ensemble that is explicitly higher-dimensional not only in its spin degrees of freedom, but also in its spatial degrees of freedom. Our investigation is motivated by two experimental results: the controlled amplification of spin noise in Rb quantum fluids [11, 12, 13], and the use of a spinor Bose-Einstein condensate as a spatially resolving magnetometer [14]. As we discuss herein, these two results, respectively, provide a complete means for inducing squeezing of both the spin vector and quadrupole moments of the quantum gas, and also a means of applying the subsequent spatially extended non-classical quantum fluid toward spatially resolved magnetometry with sub shot-noise sensitivity.

## 1 Single-mode squeezing for arbitrary spin

We begin by introducing a generalized treatment of spin squeezing that applies equally to ensembles of spin-1/2 and also of higher-spin particles. Such squeezing will be defined with respect to a coherent spin state [2], in which all particles are prepared in an identical single-particle spin- state . Consider an orthogonal basis of the single-particle spin space that includes the state and also the states . We now define many-body observables and where the summation index is taken over all particles in the -particle ensemble. These spin fluctuation operators describe the subspace of spin excitations orthogonal to the coherent spin state. With respect to this uncorrelated state, these observables satisfy the commutation relations . Thus, measurements of each pair independently obey an uncertainty relation . Following the treatment of squeezing for vector spin operators, we identify one mode of spin fluctuations, labeled by , to be spin squeezed when the measurement variance for one quadrature of the - plane is below the standard quantum limit, .

Applying this treatment to an ensemble of spin-1/2 particles, we recall that any coherent spin state can be written as the eigenstate of a projection of the dimensionless vector spin. Taking that projection to lie along , we identify the one mode of spin fluctuation operators as the Pauli operators and . More generally, for an ensemble of spin- particles prepared in the eigenstate of the spin operator, the spin fluctuation operators defined by the state act on as times the spin vector operators and . For this one mode of spin-fluctuation operators, the commutation relation matches that of the spin-vector operators. Thus, we recover the conventional description of vector-spin squeezing as pertaining to just one mode of spin fluctuations atop a maximum-spin coherent spin state.

To illustrate further the possibility of squeezing in several spin fluctuation modes, we consider spin squeezing of spin-1 particles prepared initially in the state. Our treatment is facilitated by working in a polar state basis, where is the zero-eigenvalue state of [15]. The pair of fluctuation operators defined according to the state () are identified as () and (), where the former is a component of the quadrupole moment tensor [16]. The coherent spin state is thus regarded as the uncorrelated vacuum state for independent spin fluctuations in the - and the - planes, and correlated states may be spin squeezed with respect to vector-spin or quadrupole-spin components, or linear combinations thereof. Finally, we note that this unconventional form of spin squeezing, about the polar state, may be mapped directly onto the more conventional vector-spin squeezing discussed above by a unitary transformation that rotates the polar basis states onto the proper eigenstates; experimentally, such generic unitary transformations upon alkali-atom hyperfine spin states are produced by rf-pulses and quadratic Zeeman energy shifts [17].

## 2 Spin squeezing in a spatially extended quantum field

We now consider such multi-spin-mode squeezing in a spatially extended, ferromagnetic gaseous Bose-Einstein condensate where, as we show below, such squeezing is produced naturally by spin-dependent interactions [18]. As in recent experiments [11, 12], these condensed atoms (at rest) are prepared initially in the common spin state . Within the volume occupied by the condensate, we define position-space spin-fluctuation measurement densities

(1) | |||||

(2) |

where is now the position-space Bose field operator for the polar state , and the index runs over the transverse spin polarizations. These observables represent components of the local nematicity (rank-2 moments of the atomic spin) and magnetization (rank-1 moments of the atomic spin) of the quantum gas, and obey the commutation relation where , and expectation values are taken with respect to the polar condensate.

These measurement density operators are used to construct operators corresponding to spatially resolved measurements of particular spin fluctuations of the quantum field. Defining the real measurement mode functions , we may define mode measurement operators as follows:

(3) |

where the common polarization index is omitted for clarity. In this definition for (and similarly for ), the local measurement weight is given as and the locally measured quadrature is defined by the angle . The commutation relation

(4) |

is non zero to the extent that the two measurement operators probe non-parallel spin-fluctuation quadratures on a common population of atoms.

We see that our present treatment subsumes the commonly discussed single-mode description of spin squeezing in atomic ensembles. That is, the equal-weight sum of atomic spin operators is measured by spatial-mode operators defined with the uniform functions and . However, our treatment easily allows for a multi-mode description of spin squeezing and collective spin measurements of a spatially extended quantum field. For example, spin squeezing in a spatially resolved measurement of the atomic spin, e.g. in the spinor-gas magnetometer of Ref. [14], may be treated by defining measurement mode functions corresponding to every resolved pixel in the magnetization-sensitive image. Our treatment applies also to optical measurements sensitive to spatially varying linear combinations of magnetization and nematicity, which are achieved by measuring different components of the linear optical susceptibility tensor [19].

## 3 Generation of spin squeezing in a ferromagnetic Bose-Einstein condensate

Now, let us consider the evolution of the paramagnetic condensate under the quantum-coherent spin dynamics produced by the spin-dependent contact interactions between atoms [20, 21]. As in previous experiments, we consider the condensate to be prepared in the uncorrelated polar state and then allowed to evolve freely at a fixed value of the quadratic Zeeman shift , which gives the difference between the average energy of the Zeeman sublevels and that of the state. Within the Bogoliubov approximation, the spin-dependent Hamiltonian may be written as follows:

(5) |

where with being the atomic mass. The spin-mixing interactions are described through where is the density of the condensate and relates to the -wave scattering lengths for binary collisions among atoms with total spin [22, 23]. For a ferromagnetic Bose gas, . Here, we restrict the spin excitations to the volume of the condensate.

We find it convenient to consider separately the real and imaginary parts of , i.e. we consider the fields and , which commute as . Considering the linearized equations of motion for these operators, obtained from given above, one identifies normal dynamical modes by defining a set of functions , spanning the condensate volume, as solutions to the following differential equation:

(6) |

It follows that and are real and that we may normalize the mode functions so that . With the aid of these functions, we define the magnon-excitation mode operators (omitting the polarization index)

(7) | |||||

(8) |

These operators are found to be canonically conjugate, satisfying the equations of motion and and commutation relations . The canonical evolution of magnon excitation mode is akin to that of a harmonic oscillator with frequency , with defining stable magnon excitations and defining unstable modes.

The action of this canonical evolution on the initial fluctuations in each normal mode may be discerned by examining the scaled covariance matrix

(9) |

where the mode operators are scaled by their initial rms fluctuations. It is important to note that this scaling circularizes the initial fluctuations of the normal mode operators, but that these initial fluctuations do not have the minimum uncertainty. One finds, however, that the fluctuations do reach their minimum product for excitations of a homogeneous condensate and approximate the minimum product for inhomogeneous condensates with densities varying only on long length scales.

The Hermitian covariance matrix may be expressed as with being diagonal with eigenvalues

(10) |

with and . We note that is positive for both stable and unstable magnon modes.

For stable magnon modes, the initial spin fluctuations rotate in the - plane and also experience periodic, bounded dilatations and contractions. In contrast, for unstable magnon modes, the spin fluctuations become progressively squeezed along one quadrature axis and amplified along the other, with the squeezed/amplified axes defined by the rotation matrix . The fluctuations and dynamical evolution of these mode operators relate directly to those of normalized spatial-mode measurement observables given as

(11) | |||

(12) |

corresponding to and , respectively, at . Spatially resolved spin measurements tailored to detecting the squeezed and amplified quadratures of the th normal mode are defined as time-dependent linear combinations of and through the rotation matrix .

### 3.1 Evolution of a homogeneous spinor condensate

For a homogeneous condensate, the evolution following a quench of
the initially paramagnetic quantum gas is expressed simply.
Translational symmetry allows one to identify normal modes *a
priori* as possessing Fourier components only with magnitude .
The normal-mode equation (6) is satisfied by the
functions , with being the
condensate volume, and with the magnon-excitation spectrum given as
where [20]. Unstable modes exist for
below a critical value , specifically within the
range .

The initial fluctuations of the normal mode operators defined with respect to the mode functions are given as and , giving . The normalized covariance matrix for modes with wave vector is given as

(13) |

We see that the dynamically unstable momentum-space modes are progressively squeezed over time (figure 1). The maximally unstable mode, with , is accessed for . For this mode, the orientation of the squeezing axis is fixed at an angle in the - plane, and the greatest squeezing is achieved. Away from this maximally squeezed mode, the orientation of the squeezing/amplification axes varies in time. At long times, defined by , spin fluctuations at dynamically unstable wave vectors are squeezed by a factor that decreases exponentially in time and is given as

(14) |

The orientation of the squeezed quadrature axis is given by the relation , with subsequently varying between 0 and .

This momentum-space treatment describes how the dynamical instabilities of a ferromagnetic spinor Bose-Einstein condensate may be used as a mode-by-mode parametric amplifier, both for the amplification of small spin fluctuations atop the polar condensate, as explored in Refs. [12, 13, 24], and also for the generation of a spin-squeezed quantum field. This field may be manipulated further to suit a particular metrological application. For example, spin squeezing may first be prepared at a particular wave vector by allowing the system to evolve for a given time at for which . Thereafter, the quadratic shift may be adjusted to at which the spin modes at wave vector are dynamically stable. Under such stable evolution, the spin fluctuations undergo rotation in the - plane, so that squeezing can be generated along any desired quadrature axis. More general control of the -dependent rotation angle could be realized by applying a -dependent optically induced quadratic Zeeman shift, e.g. by taking advantage of the band structure of an appropriately tuned and polarized optical lattice.

## 4 Locality in spin squeezing and measurement

It is illuminating to consider the evolution of a spatially extended quantum field in position space. The spin-dependent contact interactions to which we have ascribed the generation of a correlated quantum field act locally; thus, we would expect spatial correlations to propagate at a finite rate across the spinor condensate. Yet, the above description for the homogeneous system (section 3.1) dealt with measurement observables that incorporate correlations between particles at infinite range. In contrast, given the local origin of spin correlations, we expect that spin squeezing should be effectively observed by local measurements of the collective atomic spin on a finite number of atoms.

One approach to tracing the spatial evolution of spin correlations is to consider correlations among the position-space spin-fluctuation measurement densities, via the covariance matrix defined as

(15) |

For the uncorrelated condensate prepared at time , we have , i.e. correlations are purely local. At later times, for the initially homogeneous condensate,

(16) |

Now, let us inquire as to the spatial range of spin-squeezing produced by spin-dependent contact interactions by asking whether finite-volume spatial-mode operators may be defined that will capture the spin squeezing produced in a homogeneous spinor condensate. We consider two possibilities. First, we define “quadrature-tuned” measurement mode functions through their Fourier transform as

(17) |

where is real. Now, following the evolution time , evaluating the covariance matrix with respect to these spin-fluctuation observables, we obtain simply

(18) |

We consider also a second, simpler, “uniform quadrature” mode function with Fourier transform

(19) |

where is the wave vector for which maximal squeezing is achieved.

Now we are in a position to consider forms of the measurement operators appropriate for capturing squeezing achieved by the spinor Bose-Einstein condensate amplifier. We examine the long-time limit (equation 14), and consider amplification under two conditions.

### 4.1 Shallow quench

Under the condition , the maximally unstable mode occurs at zero wave vector. Expanding about , we approximate

(20) |

We note that , proportional to , has the units of a diffusion constant. Let us now define the measurement operator through

(21) |

so that the quantity is a 3D Gaussian with peak value of unity at the origin, and with a volume of .

Evaluating the correlation matrix, we now find fluctuations of the operator are squeezed with respect to their variance at . An estimate for the squeezing achieved using the quadrature-tuned mode function may be obtained analytically. If we assume that , then the reduction in the variance is given approximately as

(22) |

That is, the spin correlations produced in the quenched spinor Bose-Einstein condensate are indeed local, contained within a volume that grows diffusively in time.

These approximations are confirmed by numerical calculations, results of which are shown in figure 2. The quadrature-tuned spin-fluctuation measurement is indeed effective in capturing the spin squeezing produced in this system within a measurement volume that grows only slowly with evolution time. In contrast, a uniform-quadrature measurement requires larger volumes to capture this squeezing. This requirement is explained by the finite differences between the quadrature axes along which squeezing is achieved at different wave vectors. The higher the degree of squeezing one wishes to detect, the more precisely must the measurement quadrature axes be tuned to the spin-squeezed axis by restriction to a narrow range of wave vectors.

### 4.2 Deep quench

For a deep quench, with the maximum instability occurs on a sphere in -space, with radius defined by . About this shell, the gain is given approximately as

(23) |

To capture this squeezing, we define a measurement mode through the following relation:

(24) |

with chosen so that at . For large , the volume defined by this function tends toward , so that the effective radius of the measurement volume is .

Now we determine the degree of spin squeezing captured by this measurement operator. Repeating calculations as above, with the requirement , we have

(25) |

Again, we see that local measurements will capture the spin squeezing produced by the spinor Bose-Einstein condensate amplifier. Here, the correlation volume grows slower than diffusion, with radius scaling only as . Numerical evaluations of are presented in figure 3, confirming the rapid convergence onto the single-mode maximal degree of squeezing for finite-volume, quadrature-tuned measurement modes.

## 5 Experimental considerations

We conclude with a discussion of experimental prospects for
realizing and detecting multi-mode spin squeezing in a quenched
ferromagnetic condensate. We focus on the case of Rb,
which has been shown experimentally
[25, 26, 27] and theoretically
[28, 29] to be ferromagnetic, with with being the Bohr radius. The behaviour of
condensates prepared in the state and quenched to
the regime of dynamic instabilities by a rapid change of have
been studied using *in situ* magnetization-sensitive imaging.
Pertinent to the present discussion, the amplification of spin
fluctuations following quenches to different values of was
measured precisely in Ref. [12], and such amplification
was found to be roughly consistent with quantum-limited
amplification of quantum fluctuations in the initial state, with a
gain, measured by the increase in the variance of spin fluctuations,
as high as +30 dB.

In contrast to previous studies on the parametric amplification of spin fluctuations [12, 24], here we have considered the squeezing that may accompany such amplification. We assume experimental conditions similar to those of Refs. [11, 12], i.e. an optically trapped Rb condensate of atoms with peak density . At this density, the characteristic time scale for spin mixing in the spinor condensate is and the characteristic length scale is . Let us estimate the maximum observable degree of spin squeezing, taking into account non-linearities, atom loss, and the spatial inhomogeneity of the condensate.

### 5.1 Single-mode nonlinearities and atom loss

Thus far, we have treated spin dynamics in a quenched spinor condensate exclusively through the Bogoliubov approximation (equation 5), i.e. neglecting high-order terms in the spin fluctuation operators and also the depletion of the condensate due to spin mixing. Within this linear, undepleted-pump approximation, we obtain the unphysical result, expressed in equation 10, that the minimum spin-quadrature variance in each unstable spin mode tends exponentially to zero. As in quantum-optical systems, sensible limits to squeezing are recovered by going beyond the linearized equations of motion and the undepleted-pump approximation.

These limits are partly elucidated in the full dynamics of a single spin mode. We consider therefore the dynamics of a ferromagnetic spinor condensate in the single-mode approximation. Such a treatment would be appropriate for a spinor condensate spanning a volume , realized experimentally with an atom number of . The full Hamiltonian is now given as

(26) |

where the Bose operators annihilate particles in polar state and in the single spatial mode of the condensate.

In this expression, the single-polarization term matches the two-axis counter-twisting Hamiltonian described in Ref. [2]. On its own, this number-conserving Hamiltonian produces squeezing that saturates with time due to depletion of the pump (the population ), reaching the Heisenberg limit with [30]. However, the Hamiltonian contains additional terms: the externally imposed and pump-dependent quadratic Zeeman shifts and inter-polarization coupling. These nonlinear terms restrict spin squeezing to sub-Heisenberg scaling, with . Nevertheless, the expected single-mode squeezing is significant. We have performed exact single-mode calculations for a 1000-atom sample, and find a reduction of spin fluctuations by as much as 20 dB.

The single-mode limit allows us to discuss an additional limitation arising from the decay of atoms from the trap. Atom loss can be modeled by adding a noise source and a density-dependent loss rate in the coherent evolution of the Bose fields [30]. We have performed numerical calculations including such terms, representing quantum noise in the evolution of a classical field as in the truncated Wigner approximation (TWA) [31]. In accord with experimental observations, we assume an initial per-atom loss rate of 1 . We find squeezing in a atom sample can still be the level of -20 dB, a result that is not surprising given the very short spin-mixing time as compared to the lifetime of trapped atoms.

### 5.2 Effects of inhomogeneity and multi-mode nonlinearities

Experimental realizations of spinor Bose-Einstein condensates also differ from our idealized treatment in that the condensates are inhomogeneous, and, thus, the adequacy of a linearized normal-mode description to long-time dynamical evolution is less clear than in the homogeneous case. To assess the influence of mixing between normal modes and also atom losses for this situation, we have performed numerical simulations of the evolution of a quasi-one-dimensional quenched spinor Bose-Einstein condensate. Such a condensate is assumed sufficiently tightly confined in two radial dimensions that the spinor wavefunction is radially uniform, while being more weakly confined in the third dimension, along which the wavefunction is allowed to vary. Specifically, we consider a condensate with Thomas-Fermi radii of 3 and 200 along the radial and axial dimensions, respectively. The central, radially averaged condensate density determines the spin mixing interaction strength to be .

The quantum evolution of this spinor condensate is calculated using the TWA to simulate the effects of quantum noise and squeezing [24]. This approximation involves adding random noise to the initial mean-field condensate wavefunction, with the noise variance determined, in this case, by the one-dimensional condensate density. Temporal evolution is thereafter considered according to classical wave equations, which include the full non-linear dynamics given by spin-mixing Hamiltonian (the non-linear Schrödinger equation) and also the effects of atom loss (through additional noise and loss terms, described above). Expectation values, variances, and correlations for various measurement operators are quantified by repeating the calculations for many instances of the initial noise. These calculations were found convergent at a temporal resolution of s and spatial resolution of m.

We present results of a numerical simulation of an initially polar, trapped, Rb spinor condensate quenched at to a quadratic Zeeman shift of . The results for several numerical runs, i.e. for three random instances of initial noise as evaluated according to the TWA, are presented in figure 4a. Dynamical instabilities following the quench to low result in amplified fluctuations in the transverse components of the magnetization and the related components of the nematicity; only one polarization of the observables, and , is shown in the figure. The strong local correlation between these observables is clearly evident, reflecting the angle of the amplified quadrature axis for the maximally unstable normal modes.

Less evident in these position-space plots is the sharp reduction in the initial fluctuations of the normal-mode measurement observables. To reveal this reduction, we identify the normal dynamical modes based on the initial condensate distribution. The temporal gain determined for the first 41 unstable modes is shown in the inset of figure 5, sorted in decreasing order of temporal gain () and labeled by an integer mode index. Normalized spatial-mode measurement operators are then formed as defined in equations 11 and 12, and the values of these operators at different evolution times are evaluated for an ensemble of 600 numerical runs. From these values, we evaluate the scaled covariance matrices numerically and determine therefrom the quadrature axes that are amplified and squeezed, and also the measurement variance along those axes.

The squares of the measurement outcomes along the amplified and squeezed quadrature axes for the three simulation runs discussed above are shown in figure 4b. Aside from pronounced variations according to the initial random values assigned to the quantum noise in these simulations, these data indicate a suppression of spin fluctuations between 20 and 30 dB.

The ensemble averaged reduction in the initial spin fluctuations, determined for each mode and at different evolution times, is shown in figure 5. We find the temporal evolution of the dynamically unstable modes causes the measurement variance of the squeezed quadrature operators to diminish exponentially in time for evolution times up to around 60 ms, after which the variance saturates and then begins growing. The minimum variance of around or -24 dB is achieved for the unstable modes with highest gain (lowest mode index) at . This maximum degree of squeezing matches well with the estimated effects of nonlinearities and atom loss presented above.

## 6 Conclusion

We have presented a formalism for the evaluation of spin squeezing in ensembles of high-spin atoms. This formalism identifies several independent planes of spin fluctuations atop an arbitrary coherent spin state, revealing a greater resource for applications of spin squeezing in quantum information and metrology than for ensembles of spin-1/2 particles. Moreover, we developed an understanding of spin squeezing in spatially extended quantum fields, where squeezing in many spatial modes can be separately generated and measured. This formalism is applicable to many contemporary experiments on spin squeezing of atomic gases or of other distributed quantum objects.

Applying this formalism, we consider a source of spin squeezing in an extended spinor Bose-Einstein condensate. Referring to recent experiments studying the evolution of Rb condensates quenched to a regime of dynamical instabilities [11, 12, 13], we show that the strong amplification of spin fluctuations demonstrated in those experiments is accompanied also by spin squeezing, i.e. that the dynamically unstable normal modes act as quantum parametric amplifiers of spin fluctuations. Treating experimental realities such as atom losses and effects of confinement, we apply analytic and numerical tools to determine that measurable spin squeezing better than 20 dB may be achieved in such systems.

This treatment is relevant also to the question of symmetry breaking in a quantum many-body system. An initially polar-state spinor Bose-Einstein condensate preserves the rotational symmetry about the alignment axis. The quantum quench considered in this work, a rapid reduction in the quadratic Zeeman shift to induce dynamical instabilities and parametric amplification of spin fluctuations, traverses the atomic system across a symmetry-breaking quantum phase transition, where the spin-dependent contact interactions favor broken symmetry states of maximum magnetization [11]. In a closed noiseless quantum system, bona fide symmetry breaking does not occur; rather, the system evolves coherently to a symmetric superposition of broken-symmetry states. The spin-squeezed states described in this work indeed represent such a coherent superposition.

The metrological significance of spin squeezed quantum fields, and also the question of symmetry breaking in quantum systems, provide strong motivation for the detection and characterization of squeezing in quenched spinor Bose gases. For this purpose, spatially resolved measurements of the magnetization of optically trapped condensates have been demonstrated [14], and optical means of measuring components of both the magnetization and the nematicity have been suggested [19]. However, sensitivity below the atomic shot noise level, required for the direct detection of squeezing, is experimentally challenging and was not achieved in those measurements.

Alternately, we propose that one could use the condensate itself as a pre-measurement amplifier. Following the generation of inhomogeneous spin squeezing with the quadratic shift at , one may quickly increase the quadratic Zeeman shift. The previously unstable dynamical modes now evolve stably as the squeezed and amplified spin-fluctuation quadratures rotate in their quadrature spaces. After their rotation, the spin fluctuations may be again parametrically amplified by returning the quadratic shift to . Spin squeezing produced by the first stage of amplification would result in a reduction of spin fluctuations observed after the second amplification stage, as compared to a single-stage amplification of vacuum spin noise. We have performed TWA-based numerical simulations which confirm that this technique is applicable to trapped spinor Bose-Einstein condensates at experimentally accessible settings.

This work was supported by the NSF and the Army Research Office with funding from the DARPA OLE program. Partial personnel and equipment support was provided by the LDRD Program of LBNL under the Dept. of Energy Contract No. DE-AC02-05CH11231. D.M.S.-K. acknowledges support of the Miller Institute for Basic Research in Science, and J.S. acknowledges the JQI-NSF-PFC for support.

## References

### References

- D.J. Wineland et al. Spin squeezing and reduced quantum noise in spectroscopy. Phys. Rev. A, 46:R6797, 1992.
- M. Kitagawa and M. Ueda. Squeezed spin states. Phys. Rev. A, 47(6):5138–5143, 1993.
- M. H. Schleier-Smith, I. D. Leroux, and V. Vuletic. States of an ensemble of two-level atoms with reduced quantum uncertainty. Phys. Rev. Lett., 104(7):073604, 2010.
- J. Appel et al. Mesoscopic atomic entanglement for precision measurements beyond the standard quantum limit. Proc. Natl. Acad. Sci. USA, 106(27):10960–10965, 2009.
- J. Esteve et al. Squeezing and entanglement in a Bose-Einstein condensate. Nature, 455(7217):1216–1219, 2008.
- J. Hald, J.L. Sorensen, C. Schori, and E.S. Polzik. Spin squeezed atoms: a macroscopic entangled ensemble created by light. Phys. Rev. Lett., 83(7):1319, 2000.
- D. V. Kupriyanov et al. Multimode entanglement of light and atomic ensembles via off-resonant coherent forward scattering. Phys. Rev. A, 71(3):032348, 2005.
- J.M. Geremia, J.K. Stockton, and H. Mabuchi. Tensor polarizability and dispersive quantum measurement of multilevel atoms. Phys. Rev. A, 73(042112), 2006.
- V. V. Yashchuk et al. Selective addressing of high-rank atomic polarization moments. Phys. Rev. Lett., 90(25):253001, 2003.
- V. M. Acosta et al. Production and detection of atomic hexadecapole at Earth’s magnetic field. Optics Express, 16(15):11423, 2008.
- L.E. Sadler et al. Spontaneous symmetry breaking in a quenched ferromagnetic spinor Bose condensate. Nature, 443:312, 2006.
- S.R. Leslie et al. Amplification of fluctuations in a spinor Bose Einstein condensate. Phys. Rev. A, 79:043631, 2009.
- C. Klempt et al. Multiresonant spinor dynamics in a Bose-Einstein condensate. Phys. Rev. Lett., 103(19):195302, 2009.
- M. Vengalattore et al. High-resolution magnetometry with a spinor Bose-Einstein condensate. Phys. Rev. Lett., 98(20):200801, 2007.
- Specifically, we use the basis , , and . With this definition, the vector spin operator takes the form where the indices run over the Cartesian coordinates.
- We adopt the definition where and denotes an angular momentum component operator in matrix form.
- Paolo Giorda, Paolo Zanardi, and Seth Lloyd. Universal quantum control in irreducible state-space sectors: Application to bosonic and spin-boson systems. Phys. Rev. A, 68(6):062320, 2003.
- We emphasize that Bose-Einstein condensation, or even Bose-Einstein statistics, are not necessary for the existence of spin squeezing in an extended, multi-mode quantum field. Much of the notation developed in this work may be extended also to squeezing with more generic constituents. Here, we focus on spinor Bose-Einstein condensates not only for notational convenience, but also to describe the role of spin-dependent interactions in creating such squeezing.
- Iacopo Carusotto and Erich J. Mueller. Imaging of spinor gases. J. Phys. B, 37:S115, 2004.
- A. Lamacraft. Quantum quenches in a spinor condensate. Phys. Rev. Lett., 98(16):160404, 2007.
- G. I. Mias, N. R. Cooper, and S. M. Girvin. Quantum noise, scaling, and domain formation in a spinor Bose-Einstein condensate. Phys. Rev. A, 77(2):023616, 2008.
- T.-L. Ho. Spinor Bose condensates in optical traps. Phys. Rev. Lett., 81:742, 1998.
- T. Ohmi and K. Machida. Bose-Einstein condensation with internal degrees of freedom in alkali atom gases. J. Phys. Soc. Jpn., 67:1822, 1998.
- J.D. Sau, S.R. Leslie, M.L. Cohen, and D.M. Stamper-Kurn. Theory of domain formation in inhomogeneous ferromagnetic dipolar condensates within the truncated Wigner approximation. Phys. Rev. A, 80:023622, 2009.
- H. Schmaljohann et al. Dynamics of F = 2 spinor Bose-Einstein condensates. Phys. Rev. Lett., 92:040402, 2004.
- M.-S. Chang et al. Observation of spinor dynamics in optically trapped Rb Bose-Einstein condensates. Phys. Rev. Lett., 92:140403, 2004.
- A. Widera et al. Precision measurement of spin-dependent interaction strengths for spin-1 and spin-2 Rb-87 atoms. New Journal of Physics, 8:152, 2006.
- N.N. Klausen, J.L. Bohn, and C.H. Greene. Nature of spinor Bose-Einstein condensates in rubidium. Phys. Rev. A, 64:053602, 2001.
- E.G.M. van Kempen, S.J.J.M.F. Kokkelmans, D.J. Heinzen, and B.J. Verhaar. Interisotope Determination of Ultracold Rubidium Interactions from Three High-Precision Experiments. Phys. Rev. Lett., 88:093201, 2002.
- A. André and M.D. Lukin. Atom correlations and spin squeezing near the Heisenberg limit: Finite-size effect and decoherence. Phys. Rev. A, 65:053819, 2002.
- C.W. Gardiner and P. Zoller. Quantum Noise. Springer-Verlag, Berlin, Heidelberg, 2000.