## A geometric pathway to scalable quantum sensing

Mattias T. Johnsson, Nabomita Roy Mukty, Daniel Burgarth, Thomas Volz and Gavin K. Brennen ∗ Center for Engineered Quantum Systems, Dept. of Physics &amp; Astronomy, Macquarie University, 2109 NSW, Australia

Entangled resources enable quantum sensing that achieves Heisenberg scaling, a quadratic improvement on the standard quantum limit, but preparing large N spin entangled states is challenging in the presence of decoherence. We present a quantum control strategy using highly nonlinear geometric phase gates which can be used for generic state or unitary synthesis on the Dicke subspace with O N ( ) or O N ( 2 ) gates respectively. The method uses a dispersive coupling of the spins to a common bosonic mode and does not require addressability, special detunings, or interactions between the spins. By using amplitude amplification our control sequence for preparing states ideal for metrology can be significantly simplified to O N ( 5 / 4 ) geometric phase gates with parameters that are more robust to mode decay. The geometrically closed path of the control operations ensures the gates are insensitive to the initial state of the mode and the sequence has built-in dynamical decoupling providing resilience to dephasing errors.

Introduction .- Quantum enhanced sensing offers the possibility of using entanglement in an essential way to measure fields with a precision superior to that which can be obtained with unentangled resources [1-4]. Entangled resources allow the measurement sensitivity to scale as 1 /N with respect to the resources applied (socalled Heisenberg scaling), compared to the 1 / √ N obtained otherwise (the standard quantum limit, or shotnoise limit) [4-6].

spins and mode. As such it is adaptable to a variety of architectures at the forefront of quantum control including NV centres in diamond, trapped ion arrays, Rydberg atoms, and superconducting qubits.

Creating large-scale entanglement in multipartite systems for the purposes of metrology is a difficult problem for a number of reasons. There is the difficulty in precisely constructing the required quantum state using realistic quantum operations, the need to protect that quantum state from decoherence and loss [7], and the problem of carrying out a number of quantum measurement operations on the state with precise control.

From a metrology perspective, there is also the issue that many schemes claim to achieve Heisenberg limit by virtue of quadratic scaling of the Fisher information of the system [8]. While this ensures that there is an observable which has an uncertainty that scales as 1 /N with respect to some resource, it does not specify what that observable is, or require it be a convenient experimentally measurable quantity. There have been attempts to address these problems in various ways. For example, to mitigate the decoherence issue, recent work has suggested using quantum error correction assisted metrology (see [9] and references therein) or phase protected metrology [10]. Such workarounds require the ability to perform complex quantum control in the former case or engineered interactions in the latter.

Here we present a state preparation scheme and measurement protocol using geometric phase gates generalizing Ref. [11]. Our protocol addresses the issues of state preparation, decoherence protection, and choice of measurement operator. It is relatively simple to engineer as it involves only the coupling of an ensemble of qubits to a common bosonic mode, e.g. a cavity or mechanical oscillation, as well as simple global control pulses on the

Unlike previous work our scheme does not require special engineering of the physical layout of the spins, special detunings for adiabatic state preparation, addressability, or direct interaction between the spins. Furthermore, it exceeds the performance of spin squeezing protocols because of the highly nonlinear nature of the geometric phase gates used. Another advantage is that due to the geometric nature of the gate, it is completely insensitive to variations or uncertainties in the rate at which the perimeter is traversed.

Furthermore, our protocol has dynamical decoupling built in, providing resilience against dephasing during the state preparation, which is the dominant source of noise in many physical implementations. While dynamical decoupling has been considered [12] in the context of the Møllmer-Sørenson geometric gate [13], our scheme extends this to a highly nonlinear geometric phase gate and a full quantum state preparation algorithm.

Method .-We consider a collection of two-level spin half systems, and define the collective raising and lowering angular momentum operators as J + = ∑ N j =1 σ + j , J -= ( J + ) † , and the components of the total angular momentum vector are J x = ( J + + J -) / 2 , J y = ( J + -J -) / 2 i, J z = ∑ j ( 0 | 〉 j 〈 0 | - | 1 〉 j 〈 1 ) | / 2. Dicke states are simultaneous eigenstates of angular momentum J and J z : | J = N/ ,J 2 z = M 〉 , M = -J, . . . , J .

Consider the measurement of a field which generates a collective spin rotation about an axis perpendicular to ˆ given z by U η ( ) = exp[ iη ( J x sin δ + J y cos δ )]. For a measurement operator O on the system, the single shot estimation of the parameter η has variance

$$( \Delta \eta ) ^ { 2 } = \frac { ( \Delta O ( \eta ) ) ^ { 2 } } { | \partial _ { \eta } \langle O ( \eta ) \rangle | ^ { 2 } }. \quad \quad ( 1 )$$

When the measured observable is O = J z 2 , the parame-

ter variance is [14]

$$( \Delta \eta ) ^ { 2 } & = \left ( ( \Delta J ^ { x 2 } ) ^ { 2 } f ( \eta ) + 4 \langle J ^ { x 2 } \rangle - 3 \langle J ^ { y 2 } \rangle - 2 \langle J ^ { z 2 } \rangle \\ & \times ( 1 + \langle J ^ { x 2 } \rangle ) + 6 \langle J ^ { z } J ^ { x 2 } J ^ { z } \rangle \right ) ( 4 ( \langle J ^ { x 2 } \rangle - \langle J ^ { z 2 } \rangle ) ^ { 2 } ) ^ { - 1 } \quad \ W ( \iota$$

with f ( η ) = (∆ J z 2 ) 2 (∆ J x 2 ) 2 tan 2 ( η ) + tan ( 2 η ). When the initial state is the Dicke state | J, 0 , the estimate angle satisfies 〉 sin 2 ( η ) = 8 〈 J z 2 ( η ) 〉 N N ( +2) and the uncertainty in the measured angle is minimized at η min = 0 such that the quantum Cram´ er-Rao bound is saturated:

$$( \Delta \eta ) ^ { 2 } = \frac { 2 } { N ( N + 2 ) }. \quad \quad ( 2 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Note it is not essential that we know the angle δ of the field direction in the ˆ x -ˆ plane [15]. y

The best known quantum algorithm for deterministically preparing a Dicke state | J, M 〉 requires O (( N/ 2 + M N ) ) gates and has a circuit depth O N ( ) [16]. This works even for a linear nearest neighbour architecture, but requires a universal gate set and full addressability. Other proposals exist [17-20], but they all suffer from drawbacks such as not scaling beyond a few spins, strong adiabaticity or geometry, contraints, requiring large initial Fock states of motional modes, and couplings causing transitions outside the Dicke space.

In our setup we assume N spins with homogeneous energy splittings described by a free Hamiltonian H 0 = ω J 0 z which can be controlled by semi-classical fields performing global rotations generated by J x , J y . Additionally, we assume the ensemble is coupled to a single quantized bosonic mode ˆ. a Our scheme requires a dispersive coupling between the n spins and the bosonic mode of the form V = ga aJ † z .

The geometric phase gate (GPG) makes use of two basic operators [11, 21], the displacement operator D α ( ) = e αa † -α a ∗ and the rotation operator R θJ ( z ) = e iθa † aJ z which perform a closed loop in the mode phase space

$$U _ { G P G } ( \theta, \phi, \chi ) & = D ( - \beta ) R ( \theta J ^ { z } ) D ( - \alpha ) R ( - \theta J ^ { z } ) \quad \text{print} \quad \text{img } t \\ & \quad \times D ( \beta ) R ( \theta J ^ { z } ) D ( \alpha ) R ( - \theta J ^ { z } ) \quad \text{img } t \\ & = e ^ { - i 2 \chi \sin ( \theta J ^ { z } + \phi ) }. \quad \text{(3)} \quad \text{angle}$$

where φ = arg( α ) -arg( β ) and χ = | αβ | .

The system and the mode are decoupled at the end of GPG cycle. Also, if the mode begins in the vacuum state, it ends in the vacuum state and the first operation R ( -θJ z ) in Eq. (3) is not needed. In the GPG it is necessary to evolve by both R θJ ( z ) and R ( -θJ z ). This can be done by conjugating with a global flip of the spins R θJ ( z ) = e -iπJ x R ( -θJ z ) e iπJ x , implying that the GPG can be generated regardless of the sign of the dispersive coupling strength g . Furthermore, because R ( ± θJ z ) commutes with H 0 at all steps, this conjugation will cancel the free evolution accumulated during the GPG.

Assuming the number of spins n is even, we consider N/ 2 sequential applications of the GPG (see also [11]):

$$& \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \\ ) ^ { - 1 } \quad & \quad W ( \ell ) = \prod _ { k = 1 } ^ { N / 2 } U _ { G P G } ( \theta _ { k }, \phi _ { k } ( \ell ), \chi ) \\ \text{ial} \\ \text{ies} \quad & \quad = \sum _ { M = - J } ^ { J } e ^ { - i 2 \chi \sum _ { k = 1 } ^ { N / 2 } \sin ( \theta _ { k } M + \phi _ { k } ( \ell ) ) } | J, M \rangle \langle J, M |, \\ \text{ed} \quad & \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

with glyph[lscript] = 0 , . . . , N .

$$\theta _ { k } = \frac { 2 \pi k } { N + 1 }, \ \phi _ { k } ( \ell ) = \frac { 2 \pi k ( N / 2 - \ell ) } { N + 1 } + \frac { \pi } { 2 }, \ \chi = \frac { \pi } { N + 1 },$$

gives W glyph[lscript] ( ) = e -iπ J,glyph[lscript] | -N/ 2 〉〈 J,glyph[lscript] -N/ 2 | meaning it applies a π phase shift on the symmetric state with glyph[lscript] excitations. For N odd we can use N GPGs with the same angles θ k , φ k ( glyph[lscript] ) as above but with χ = π 2( N +1) .

Given the control toolbox above of global rotations and the GPGs, one can synthesize an arbitrary unitary operator on the Dicke subspace. Writing U = ∑ N k =1 e iλ k | λ k 〉〈 λ k | , where {| λ k 〉} form an orthonormal basis on the Dicke subspace, and λ k ∈ R (note since the global phase is irrelevant we have set λ 2 J +1 = 0). This unitary can be decomposed as U = ∏ N -1 k =1 [ K λ ( k ) e iλ k | J, - 〉〈 J J, - | J K † ( λ k )], where K λ ( k ) is any unitary extension of the state synthesis mapping | J, - 〉 J → | λ k 〉 . The phasing unitary is the same as W (0) but with χ → λ / N k ( + 1). To construct K λ ( k ), we find the decomposition: ˜ K = [ ∏ N -1 s =1 e iβ s J y U GPG ( θ , φ s s , χ s )] e -iJ y π 2 U GPG ( π 2 , 0 , π 4 ) e iJ y π 2 gives very accurate results when optimized over the 4 N -4 free parameters { β , θ s s , φ s , χ s } [22]. The overall complexity in GPG count using this approach is N for state synthesis and 5 N / 2 2 for unitary synthesis.

While the general state synthesis approach above can be used for building Dicke states, the N -1 action angles { χ s } that optimize state fidelity are O (1) and this has implications for noise due to mode decay as described below. We now describe a way, based on amplitude amplification, to improve matters by only using the GPGs that appear in W glyph[lscript] ( ) that have action angles O (1 /N ), and with the added advantage of providing an analytical solution to the Dicke state synthesis problem. Starting with all spins down, i.e. in | J, - 〉 J , let the target state be | w 〉 = | J, 0 〉 . We will make use of the operators U w = e -iπ w | 〉〈 w | = W N/ ( 2) and U s = e -iπ s | 〉〈 s | = e iJ y π/ 2 W (0) e -iJ y π/ 2 . In total the operators U w and U s each use N/ 2 GPGs. The orbit of the initial state | s 〉 under the operators U w and U s , is restricted to a subspace spanned by the orthonormal states | w 〉 and | s ′ 〉 = | s 〉-| w 〉〈 w s | 〉 √ 1 -|〈 w s | 〉| 2 , exactly as in Grover's algorithm. The composite pulse is one Grover step U G = U U s w . Geometrically, relative to the state | s ′ 〉 , the initial state | s 〉 is rotated by an angle δ/ 2 toward | w 〉 , where δ = 2sin -1 ( |〈 w s | 〉| ) , and after each

Grover step is rotated a further angle δ toward the target. The optimal number of Grover iterations to reach the target is # G = ⌊ π 4 |〈 w s | 〉| ⌋ where the relevant overlap is 〈 w s | 〉 = d J -J, 0 ( π/ 2) = 2 -J √ (2 J )! /J !. For J glyph[greatermuch] 1 we have 〈 w s | 〉 ≈ ( πJ ) -1 / 4 . Then the optimal number of Grover steps is

$$\# G = \lfloor \pi ^ { 5 / 4 } N ^ { 1 / 4 } / 2 ^ { 9 / 4 } \rfloor. \text{ \quad \ \ } ( 4 )$$

The fidelity overlap of the output state ρ of the protocol with the target state is F = Tr( | w 〉〈 w ρ | ). For the Grover method it is easily calculated as

$$F = \sin ^ { 2 } ( ( \# G + 1 / 2 ) \delta ) > 1 - \sqrt { 2 / \pi N }. \quad ( 5 )$$

While the fidelity error falls off at least as fast as √ 2 /πN for all N glyph[greatermuch] 1, if N is near a value where the argument in Eq. (4) is a half integer, i.e. glyph[ceilingleft] 32(2 k + 1) 4 /π 5 glyph[ceilingright] , with k ∈ Z , the error will be much lower. For example, at N = (10 70 260 700 1552) the fidelity error is (1 84 , , , , . × 10 -4 , 1 57 . × 10 -5 , 1 68 . × 10 -6 , 3 65 . × 10 -8 , 1 92 . × 10 -8 ). The number of GPGs needed to prepare the Dicke state by the Grover method is c × N 5 / 4 with a constant c &lt; 1.

The effectiveness of our scheme when used for metrology is shown in Figure 1, which shows the precision ∆ η obtainable as a function of N , compared to that obtained from both the standard quantum limit and the ultimate Cram´ er-Rao bound. It also shows the fidelity obtainable as a function of the number of spins N . The achievable fidelity is clearly optimized for specific spin values.

We have focused on preparing the state | J, 0 , 〉 but with simple modifications our protocol can prepare any Dicke state | J, M 〉 . First use the initial state | s 〉 = e iglyph[epsilon1] M J y | J, - 〉 J , and second substitute the operators U w = W M ( + N/ 2) and U s = e iglyph[epsilon1] M J y W (0) e -iglyph[epsilon1] M J y where glyph[epsilon1] M = cos -1 ( M/J ). Now the relevant overlap is |〈 w s | 〉| = | d J M, -J ( -glyph[epsilon1] M ) , | and for J - | M | glyph[greatermuch] 1, | d J M, -J ( -glyph[epsilon1] ) | ≈ ( √ πJ sin glyph[epsilon1] M ) -1 / 2 [23], implying # G = O N ( 1 / 4 ).

Finally, measurement of J z 2 could be done by again coupling the spins to the bosonic mode but now with a linear coupling V m = g m J z ( a † + a ) which generates a mode displacement depending on the collective spin projection. When the mode is in a number diagonal state, e.g. a thermal state, with mean excitation number ¯, the measurement of ˆ after a coupling time n n τ is 〈 ˆ n 〉 = ¯ + ( n g m τ ) 2 〈 J z 2 〉 . If mode number and quadratic spin operator measurements are difficult there are alternatives. One is estimation by a classical average over p experiments: E J [ 〈 z 2 〉 ] = ∑ p k =1 M k ( ) 2 p , where M k ( ) is the outcome of the k th measurement of J z [24]. Another is, after accumulating the signal, one could invert the state preparation and measure | J, - 〉〈 J J, - | J which gives the same precision scaling as Eq. (1), but doesn't require number resolved excitation counts.

There will be errors due to decay of the bosonic mode during the operations, as well as decoherence due to en-

Δ

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01120v5-figure-5a0adddb-e336-49b1-aee9-a9b1dcdc8eee.jpg)

5

10

50

100

5001000

5

10

50

100

5001000

n

N

n

N

n

FIG. 1. Measurement precision ∆ η as a function of number of spins (log-log scale). Shown are the shot noise limit, the quantum Cram´ er-Rao bound Eq. (2), and this protocol (blue). Inset: Fidelity for sets of ensemble sizes using the same number of Grover steps, which grow as N 1 / 4 .

vironmental coupling to the spins, which will degrade the fidelity. We now address these.

Mode damping .- We treat the mode as an open quantum system with decay rate κ . In order to disentangle the spins from the mode, the third and fourth displacement stages of the k -th GPG should be modified to D ( -α k ) → D ( -α e k -κθ k /g ) and D ( -β k ) → D ( -β e k -κθ k /g ). For simplicity we choose | α k | = | β k | . For an input spin state in the symmetric Dicke space ρ = ∑ M,M ′ ρ M,M ′ | J, M 〉〈 J, M ′ | , the process for the k -th GPG with decay on the spins, including the modified displacement operations above, is [43]

$$& \stackrel { = } { = } \quad \mathcal { E } ^ { ( k ) } ( \rho ) \ = \ U _ { G P G } ( \theta _ { k }, \phi _ { k }, \chi _ { k } ) \left [ \sum _ { M, M ^ { \prime } } R _ { M, M ^ { \prime } } ^ { ( k ) } \rho _ { M, M ^ { \prime } }. \\ & \stackrel { | J, M \rangle \langle J, M ^ { \prime } | \right ] \times U _ { G P G } ^ { \dagger } ( \theta _ { k }, \phi _ { k }, \chi _ { k } ) } { = }$$

where χ k = | α k | 2 f ( θ k ), f ( θ k ) = ( e -3 θ k κ/ g 2 + e -θ k κ/ g 2 ) / 2 and Γ M,M ′ and ∆ M,M ′ are given in the Supplementary Material [26].

The full operation is a concatenation of these imperfect processes E ( k ) ( ρ ), and we characterise its accuracy with the process fidelity F pro ( E , U ), which measures how close a quantum operation E is to the ideal operation U [28].

For each GPG we can readily find the lower bound on the process fidelity (see Supplementary Material [26])

$$\text{cr} _ { - } \quad F _ { \text{pro} } ( \mathcal { E }, U _ { \text{GPG} } ) > e ^ { - 6 \pi \chi _ { k } \kappa / f ( \theta _ { k } ) g } \cos ( \chi _ { k } 4 \pi \kappa / f ( \theta _ { k } ) g ).$$

Note the scaling of the exponent with 1 /N since the action angles χ are O (1 /N ). For the composite phasing map, numerically we find the tighter bound for the fidelity F pro ( E , W ( glyph[lscript] )) &gt; e -π κ/g 2 which is notably independent of N .

Dephasing: We next address spin decoherence. We assume that amplitude damping due to spin relaxation is small by the choice of encoding. This can be accommodated by choosing qubit states with very long decay times

either as a result of selection rules, or by being far detuned from fast spin exchange transitions. Hence we will focus on dephasing. Due to the cyclic evolution during each GPG, there is error tolerance to dephasing because if the interaction strength between the system and environment is small compared to g , then the spin flip pulses used between each pair of dispersive gates R θa a ( † ) will echo out this noise to low order.

For a given input density matrix ρ (0), the output after a total time T has off-diagonal matrix elements that decay as ρ M,M ′ ( T ) = ρ M,M ′ (0) e -( M M -′ ) 2 A T ( ) . For the global dephasing map the numbers M,M ′ ∈ -[ N/ ,N/ 2 2] are in the collective Dicke basis, while for local dephasing it is with respect to a local basis M,M ′ ∈ [ -1 / 2 1 , / 2]. Our argument for suppression of dephasing works for both cases. Global dephasing is the most deleterious form of noise when the state has large support over coherences in the Dicke subspace, due to decay rates that scale quadratically in the difference in M number. However, it leaves the total Dicke space, and in particular the target Dicke state, invariant. Local dephasing induces coupling outside the Dicke space, but with a rate that is at most linear in N .

Consider the evolution during the N/ 2 control pulses to realize either of the phasing gates U s or U w . Assuming Gaussian bath statistics, the effective dephasing rate can be written as the overlap of the noise spectrum S ω ( ) and the filter function | f ( ω ) | 2 [47, 48]. As shown in the Supplementary material [26], for our pulse sequence to lowest order in ω/g we find

$$g ^ { 2 } | f ( \omega ) | ^ { 2 } \approx \frac { ( \omega / g ) ^ { 2 } \pi ^ { 4 } N ^ { 2 } ( N + 2 ) ^ { 2 } } { 9 ( N + 1 ) ^ { 2 } }. \quad \ \ ( 6 ) \quad \ 0. c$$

Comparison with the case where no spin flips are applied is plotted in Fig. 2 showing there is substantial decoupling from the dephasing environment when the spectral density has dominant support in the range ω &lt; g/ 2.

The freedom to apply the GPGs in any order allows for further improvement. Consider coupling to a zero temperature Ohmic bath with noise spectrum S ω ( ) = αωe -ω/ω c with cutoff frequency ω /g c = 0 1. For . N = 20, the ratio of the effective decay rate for the linearly ordered sequence of GPGs above to that with no decoupling is A T /A ( ) 0 ( T ) = 0 0085. . However, by sampling over permutations of the ordering of GPGs we find a sequence [49] achieving A T /A ( ) 0 ( T ) = 0 0026. . Examples of the effectiveness of our decoupling protocol on sensitivity for various values of A T ( ) are shown in Figure 3. Effectiveness on the fidelities can be found in the Supplementary Material [26].

Error-tolerant states .- The state preparation method we have described has inherent tolerance to decoherence. However once the state is prepared further errors such as qubit loss or dephasing could accumulate. Strategies to address this by using superpositions of Dicke states were

2

2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01120v5-figure-1db64683-fde3-48bd-899a-cb36df854865.jpg)

FIG. 2. Suppression of dephasing via dynamical decoupling inherent in the sequence of GPGs used for each of the operators U s and U w . Solid curves are filter functions using the GPGs. Dashed curves are plots of Eq. (26), which is a good approximation for ω/g &lt; 1 /πN . Short-dashed curves are the bare case without decoupling. (Red, green, blue) curves correspond to n = (10 100 1000) spins. , ,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01120v5-figure-061838ab-c306-4071-9392-f3fad97a464d.jpg)

g

FIG. 3. Precision obtained for 70 spins with a single measurement of J z 2 as a function of mode decay for various global dephasing factors A T ( ): no global dephasing (blue line), A T ( ) = 10 -6 (blue dashed), A T ( ) = 10 -5 (blue dotdashed), A T ( ) = 10 -4 (blue dotted). These dephasings correspond to an underlying decoherence rate of γ gdp = 10 -4 g accumulated over each phasing gate of duration T . For a zero temperature Ohmic bath, the corresponding cuttoff frequencies are: ω /g c = { 0 003 0 022 0 094 . , . , . } . Black line: performance without dynamical decoupling with A 0 ( T ) = 0 0223. .

recently proposed [44]. The states considered were

$$| \varphi _ { u } \rangle = \frac { 1 } { \sqrt { 2 ^ { n } } } \sum _ { j = 0 } ^ { n } \sqrt { n \choose j } | J = \frac { k n u } { 2 }, M = k j - J \rangle.$$

The number of spins N = k × n × u , and the parameters u and n determine the robustness of the states to some number of loss and dephasing errors respectively. A state performing well in the presence of one

erasure error is u = 1 , n = 2 , k = N/ 2, which can be written | ϕ 1 〉 = 1 2 ( | J, - 〉 J + √ 2 | J, 0 〉 + | J, J 〉 ) while u = 2 , n = 1 , k = N/ 2 tolerates one dephasing error and can be written | ϕ 2 〉 = 1 √ 2 ( | J, - 〉 J + | J, 0 ). 〉 Both of these states can be prepared using our protocol. The specific steps are given in the Supplementary Material [26]. Superpositions of Dicke states also arise as code words of so-called permutationally invariant quantum codes [46], and of recently discovered codes which admit Gaussian Clifford operations [34] which our method could also prepare.

Implementations .- The scheme we have presented is amenable to a variety of architectures which allow collective dispersive couplings between spins and an oscillator. These include: trapped Rydberg atoms coupled to an microwave cavity [35, 36], trapped ions coupled to a common motional mode [37] or to an optical cavity mode [38], superconducting qubits coupled to microwave resonators [39], and NV centres in diamond coupled to a microwave mode inside a superconducting transmission line cavity [40].

One contender to test our scheme is Rydberg atoms coupled to microwave cavities. Recently the dispersive detection of small atomic Rydberg ensembles coupled to a high-Q microwave cavity was reported [36]. Their numbers suggest a ratio of κ/g ≈ 0 8 (with . κ = 2 π × 11 8 kHz . and g = 2 π × 14 3 . kHz). The collective coupling rate observed was a few MHz, suggesting an additional pathway to improving γ/g by orders of magnitude by encoding spins through collective subensembles. Consider an encoding where each spin is composed of n physical spins with logical states | 0 〉 = | j = n/ , 2 - 〉 j and | 1 〉 = | j = n/ , 2 -j +1 , i.e. 〉 the permutationally invariant states of zero or one excitation shared among the n spins. If the spins within each logical qubit interact, e.g. via dipole-dipole interactions, there will be a dipoleblockade to larger numbers of excitations. Hence collective rotations frequency-tuned to the transition energy E 1 -E 0 will be collective but only act on this qubit subspace. The dispersive interaction strength is enhanced by g → g √ n , provided the number n is similar for all logical spins. Using this kind of encoding, dispersive coupling with strength g ≈ 2 π × 2 2 MHz was obtained with NV . ensembles in diamond bonded onto a transmission line resonator with quality factor Q ≈ 4300 at the first harmonic frequency ω 1 = 2 π × 2 75 GHz. . Microwave cavities with much higher quality factors, e.g. Q ≈ 10 , have 9 been realized [42] which for the same dispersive coupling could provide κ/g ≈ 10 -6 .

We acknowledge helpful discussions with Jason Twamley and Yingkai Ouyang. This research was funded in part by the Australian Research Council Centre of Excellence for Engineered Quantum Systems (Project number CE170100009).

## SUPPLEMENTAL INFORMATION

## Gate fidelity with cavity decay

Cavity field decay at a rate κ acts as a source of error for the many body interactions which the cavity mode mediates. Consider the the joint evolution of the spins and the mode. The coupling of the mode to its environment is treated as irreversible and thus can be described by the standard master equation in Lindblad form. The equation of motion for the joint state is

$$\mathrm d \text{ is } \quad \dot { \rho } ( t ) & = \mathcal { L } ( \rho ( t ) ) \\ \text{cell} _ { \text{cell} } \\ \text{led} \\ \mathcal { I } \text{ to }$$

The evolution due to decay conserves the quantum number J , and it will be convenient to compute the adjoint action on on a joint state state of the spins and mode with Heisenberg evolved operators e L t A M,M ′ (0) where:

$$A ^ { M, M ^ { \prime } } ( t ) \equiv | J, M \rangle \langle J, M ^ { \prime } | \otimes | \alpha _ { M } \rangle \langle \beta _ { M ^ { \prime } } | ( t ).$$

The solutions are easily verified to be given by

$$\underset { \mathbb { H } z } { \text{l to} } \underset { n = 0 } { A ^ { M, M ^ { \prime } } ( t ) } = \sum _ { n = 0 } ^ { \infty } \frac { b _ { M M ^ { \prime } } ^ { n } ( t ) } { n! } e ^ { - ( i g M + \kappa / 2 ) a ^ { \dagger } a t } \\ \underset { \mathbb { t h } } { \text{late} } \underset { \mathbb { t h } } { \times } a ^ { n } A ^ { M, M ^ { \prime } } ( 0 ) ( a ^ { \dagger } ) ^ { n } e ^ { ( i g M ^ { \prime } - \kappa / 2 ) a ^ { \dagger } a t } \quad ( 8 )$$

where

$$b _ { M M ^ { \prime } } ( t ) = \frac { \kappa ( 1 - e ^ { - [ \kappa + i g ( M - M ^ { \prime } ) ] t } ) } { \kappa + i g ( M - M ^ { \prime } ) }. \quad \ ( 9 )$$

The evolved state is then

$$\rho ( t ) = e ^ { \mathcal { L } t } \rho ( 0 ) = \sum _ { M,, M ^ { \prime } } \rho _ { M, M ^ { \prime } } ( 0 ) A ^ { M, M ^ { \prime } } ( t ).$$

In order to evaluate the effect of cavity decay during the the geometric phase gate, we are particularly interested in the case where initially A M,M ′ (0) = | M M 〉〈 ′ | ⊗ | α M 〉〈 β M ′ | , with | α M 〉 , | β M ′ 〉 coherent states. This kind of factorization is true at any stage of spin coupling to the field. Using Eq. (8), the sum becomes an exponential and the evolved state is

$$\intertext { a } _ { \text{av} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \prod _ { \text{$\alpha$} } \text{$\alpha$} \\ \text{son} \\ \text{ad} \quad.$$

where

$$\Gamma _ { M, M ^ { \prime } } ( t ) = \alpha _ { M } \beta _ { M ^ { \prime } } ^ { * } b _ { M, M ^ { \prime } } ( t ) - ( | \alpha _ { M } | ^ { 2 } + | \beta _ { M ^ { \prime } } | ^ { 2 } ) \frac { 1 - e ^ { - \kappa t } } { 2 }.$$

We ignore decay during the displacement stages of the evolution (i.e. we assume these are done quickly relative to the decay rate), and we assume that the system particles do not interact with the field during these steps. For simplicity we evaluate the performance when the cavity begins in the vacuum state, in which case there are seven time steps to consider:

$$D ( - \beta ) e ^ { - i \tau _ { 5 } V } D ( - \alpha ) e ^ { i \tau _ { 3 } V } D ( \beta ) e ^ { - i \tau _ { 1 } V } D ( \alpha ).$$

Let τ 5 = τ 3 = τ 1 so that the periods of spin field coupling are all equal in duration. In order that the field state return to the vacuum at the end of the sequence, we choose α ′ = α -κτ 1 , β ′ = β -κτ 1 for the parameters of the second two displacement operators. The total sequence then yields the output state:

where we defined R M,M ′ = e d M,M ′ ( t 2 )+ d M,M ′ ( t 4 )+ d M,M ′ ( t 6 ) and χ = | αβ | ( e -3 κτ / 1 2 + e -κτ / 1 2 ) / 2. This can be interpreted as coherent evolution with an evolution operator

$$U = e ^ { - i 2 \chi \sin ( \phi + \theta J ^ {$$

where θ = gτ 1 , followed by further evolution diagonal in the { M } basis and dephasing. Matrix elements diagonal in M are invariant.

For simplicity, we assume | α | = | β | , g &gt; 0 and write θ = gτ 1 . The factor R M,M ′ that dictates the deviation from perfect evolution can be written

$$\rho _ { \text{out} } \ = \ \sum _ { M, M ^$$

$$R _ { M, M ^ { \prime } } = e ^ { - \Gamma _ {$$

where Γ M,M ′ and ∆ M,M ′ are real, and explicitly are

$$\begin{array} {$$

$$\stackrel { \wedge \cdot \cdot \cdot } { \Delta _ {$$

(15)

Notice, Γ M,M ′ = Γ M ,M ′ and Γ M,M = 0 and also ∆ M,M ′ = -∆ M ,M ′ . An expansion up to first order in κ g yields simplified expressions:

$$\Gamma _ { M, M ^ { \prime } } \ = \ \frac { | \alpha | ^ { 2 } \frac { \sigma } { M - M ^ { \prime } } } { M - M ^ { \prime } } \left ( 2 \sin ( \theta M ^ { \prime } + \phi ) - \theta M ^ { \prime } ( \cos ( \theta M ^ { \prime } + \phi ) + \cos ( \theta M + \phi ) + 4 ) + \theta M ( \cos ( \theta M ^ { \prime } + \phi ) ) \\ - 4 ( \sin ( \theta ( M - M ^ { \prime } ) ) ) - 2 ( \sin ( \theta M + \phi ) ) + \theta M ( \cos ( \theta M + \phi ) ) + 4 \theta M \right )$$

$$\Delta _ { M, M ^ { \prime } } \ = \ \begin{array} { c c } | \alpha | ^ { 2 } \theta ^ { \geq } _ { \widehat { g } } ( - \sin ( \theta M ^ { \prime } ) + i ( \cos ( \theta M ^ { \prime } ) ) + \sin ( \theta M ) - i ( \cos ( \theta M ) ) ) ( \cos ( \theta M ^ { \prime } + \theta M + \phi ) \\ - i ( \sin ( \theta M ^ { \prime } + \theta M + \phi ) ) ) ( i ( \sin ( \theta ( M ^ { \prime } + M ) + 2 \phi ) ) + \cos ( \theta ( M ^ { \prime } + M ) + 2 \phi ) + 1 ) \end{array}$$

| α | 2 4 πκ/g .

In order to obtain a more rigorous process fidelity bound, we proceed as follows. For an input spin state in the symmetric Dicke space ρ = ∑ M,M ′ ρ M,M ′ | J, M 〉〈 J, M ′ | , the process for the k -th GPG with decay on the spins, including the modified displacement operations above, is [43]

$$\mathcal { E } ^ { ( k ) } ( \rho ) \ = \ \begin{array} { c } U _ { G P G } ( \theta _ { k }, \phi _ { k }, \chi _ { k } ) \left [ \sum _ { M, M ^ { \prime } } R _ { M, M ^ { \prime } } ^ { ( k ) } \rho _ { M, M ^ { \prime } } \\ | J, M \rangle \langle J, M ^ { \prime } | \right ] \times U _ { G P G } ^ { \dagger } ( \theta _ { k }, \phi _ { k }, \chi _ { k } ) \\ ( 1 8 ) \end{array} \quad \text{Henc}$$

where

$$\chi _ { k } = | \alpha _ { k } | ^ { 2 } ( e ^ { - 3 \theta _ { k } \kappa / 2 g } + e ^ { - \theta _ { k } \kappa / 2 g } ) / 2$$

and

$$R _ { M, M ^ { \prime } } ^ { ( k ) } = e ^ { - \Gamma _ { M, M ^ { \prime } } ( \theta _ { k }, \phi _ { k }, \alpha _ { k } ) } e ^ { i \Delta _ { M, M ^ { \prime } } ( \theta _ { k }, \phi _ { k }, \alpha _ { k } ) }$$

with the factors Γ M,M ′ and ∆ M,M ′ given above. Now if we adjust α such that on the k -th GPG, χ k = π/ N ( +1), then we have

$$| \alpha _ { k } | ^ { 2 } = \frac { 2 \pi } { ( N + 1 ) ( e ^ { - 3 \theta _ { k } \kappa / 2 g } + e ^ { - \theta _ { k } \kappa / 2 g } ) }.$$

Because the coherent and decoherent maps for different GPGs commute, the entire sequence that phases a Dicke state according to W glyph[lscript] ( ) is:

$$\mathcal { E } ( \rho ) \ & = \ \mathcal { E } ^ { ( N / 2 ) } \circ \cdots \circ \mathcal { E } ^ { ( 1 ) } ( \rho ) \quad \text{The } ] \\ & = \ W ( \ell ) \sum _ { M, M ^ { \prime } } \Upsilon _ { M, M ^ { \prime } } ( \ell ) \rho _ { M, M ^ { \prime } } | J, M \rangle \langle J, M ^ { \prime } | W ( \ell ) ^ { \dagger } \quad \text{pressed} \\ \intertext { r u i d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d o d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d e d s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t o d e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n e t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t i n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e o n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n s s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n d e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t o n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e s s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e s s s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e e n s t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t o n t o n t e n t e n t e n t e n t e n t e n t e n t e n t e n t e n t$$

This describes ideal evolution followed by a nonlinear dephasing map, where the decoherence factor is

$$\Upsilon _ { M, M ^ { \prime } } ( \ell ) \ & = \ \prod _ { k = 1 } ^ { N / 2 } R _ { M, M ^ { \prime } } ^ { ( k ) } \\ & = \ \exp \left [ \sum _ { k = 1 } ^ { N / 2 } ( - \Gamma _ { M, M ^ { \prime } } ( \theta _ { k }, \phi _ { k } ( \ell ), \alpha _ { k } ) \\ & \quad + i \Delta _ { M, M ^ { \prime } } ( \theta _ { k }, \phi _ { k } ( \ell ), \alpha _ { k } ) ) \right ].$$

The process fidelity F pro ( E , U ) measures how close a quantum operation E is to the ideal operation U as measured by some suitable metric. The fidelity measure we use is the overlap between the induced Jamioglyph[suppress]kowskil Choi state representations of the operations. The process fidelity is readily computed using the fact that the noise map E ( ρ S (0)) commutes with the target unitary U . Hence, we can compute the fidelity which measures how close the noisy map E ′ ( ρ S (0)) = U † E ( ρ S (0)) U is to the ideal operation, i.e. the identity operation:

$$F _ { \text{pro} } ( \mathcal { E }, U ) = F _ { \text{pro} } ( \mathcal { E } ^ { \prime }, \mathcal { I } ) = _ { S, S ^ { \prime } } \langle \Phi ^ { + } | \rho _ { \mathcal { E } ^ { \prime } } | \Phi ^ { + } \rangle _ { S, S ^ { \prime } }.$$

where

$$| \Phi ^ { + } \rangle _ { S, S ^ { \prime } } = \frac { 1 } { \sqrt { D } } \sum _ { M } | J, M \rangle _ { S } \otimes | J, M \rangle _ { S ^ { \prime } }, \quad ( 2 0 ) \quad \text{                                                                                                                                                                                                        
                                                                                                                                                                                                        where$$

Here we are computing the overlap of the Jamioglyph[suppress] lkowskiChoi representations of the maps as states in the Hilbert space H ⊗H S S ′ containing our system space and a copy each with dimension D :

$$\stackrel { \iota \cdot } { = } \int _ { \substack { \rho \mathcal { E } ^ { \prime } } } & = \mathcal { I } _ { S } \otimes \mathcal { E } ^ { \prime } _ { S ^ { \prime } } ( | \Phi ^ { + } \rangle _ { S, S ^ { \prime } } ) \\ \text{h} \\ \text{d}$$

Hence

$$F _ { \text{pro} } ( \mathcal { E }, W ( \ell ) ) = \frac { 1 } { ( N + 1 ) ^ { 2 } } \sum _ { M, M ^ { \prime } = - J } ^ { J } \Upsilon _ { M, M ^ { \prime } } ( \ell ),$$

From the limits above, we readily find the process fidelity bound

$$F _ { \text{pro} } ( \mathcal { E }, U _ { \text{GPG} } ) > e ^ { - 6 \pi | \alpha | ^ { 2 } \kappa / g } \cos ( | \alpha | ^ { 2 } 4 \pi \kappa / g ). \quad ( 2 1 )$$

Finally we note that numerical simluation suggests a tighter bound given by

$$F _ { \text{pro} } ( \mathcal { E }, W ( \ell ) ) > e ^ { - \pi ^ { 2 } \kappa / g }. \text{ \quad \ \ } ( 2 2 )$$

Notably, this fidelity is independent of N .

## A loose upper bound on precision as function of fidelity

The precision of the estimation parameter η is expressed as

$$( 1 9 ) \text{ar de} ^ { - } ( \Delta \eta ) ^ { 2 } = & ( ( \Delta J ^ { x 2 } ) ^ { 2 } f ( \eta ) + 4 \langle J ^ { x 2 } \rangle - 3 \langle J ^ { y 2 } \rangle - 2 \langle J ^ { z 2 } \rangle \\ & \times ( 1 + \langle J ^ { x 2 } \rangle ) + 6 \langle J ^ { z } J ^ { x 2 } J ^ { z } \rangle ) ( 4 ( \langle J ^ { x 2 } \rangle - \langle J ^ { z 2 } \rangle ) ^ { 2 } ) ^ { - 1 }$$

To check how the precision scales in relation to the fidelity, F , we can calculate the precision assuming an input density matrix, ρ = a J, | 0 〉〈 J, 0 | + b glyph[notforces] , where | J, 0 〉 is our ideal Dicke state, glyph[notforces] is the identity matrix of size ( N +1) × ( N +1) and a and b are related to fidelity as a = (1 + 1 /N F ) -1 /N and b = (1 -F /N ) . We choose this form for the input density matrix so that applying a global dephasing map to the output state of our protocol would make it diagonal in the | J, M 〉 basis but keep the population in | J, 0 〉 constant. Assuming the diagonal matrix elements (except the | J, 0 〉〈 J, 0 | entry) are equally weighted is a maximally unbiased assumption. After calculating the variances and expectation values of the angular momentum operators as they appear in Equation (23), then taking the high fidelity limit F → 1 and assuming large N , the precision is found to be

$$( \Delta \eta ) ^ { 2 } = 2 / ( N ( N + 2 ) ) + \sqrt { ( 1 - F ) / 1 0 }, \quad ( 2 4 )$$

where 2 / N N ( ( +2)) is the Cram´ er-Rao bound. Numerically we find this approximate form is extremely good for 1 -F &lt; 10 -2 .

glyph[negationslash]

For our choice of the input density matrix, there appears to be a lower bound to the precision (as a function of N ) that is set by the fidelity. If we would want the overall expression to fall off as 1 /N 2 , i.e to achieve the Heisenberg scaling, then we would need the error 1 -F to scale as 1 /N 4 . While this requirement is demanding in terms of performance, it should be noted that we have assumed that all the populations in | J, M 〉 , M = 0 are equal when in fact the non-zero terms in the output density matrix are much more concentrated near the target state | J, 0 〉 for our protocol. As the precision involves terms like expectation values of J z 4 , error terms with support on states far away from | J, 0 〉 will give large errors. Thus, we are overestimating the error in this case and this should be viewed as a loose upper bound on the precision.

## Dephasing rates and compensation via dynamical decoupling

Due to the cyclic evolution during each GPG, there is error tolerance to dephasing because if the interaction strength between the system and environment is small compared to g , then the spin flip pulses used between each pair of dispersive gates R θa a ( † ) will echo out this noise to low order.

For a given input density matrix ρ (0), the output after a total time T has off-diagonal matrix elements that decay as ρ M,M ′ ( T ) = ρ M,M ′ (0) e -( M M -′ ) 2 A T ( ) . For the global dephasing map the numbers M,M ′ ∈ -[ N/ ,N/ 2 2] are in the collective Dicke basis, while for local dephasing it is with respect to a local basis M,M ′ ∈ [ -1 / 2 1 , / 2]. Our argument for suppression of dephasing works for both cases. Global dephasing is the most deleterious form of noise when the state has large support over coherences in the Dicke subspace, due to decay rates that scale quadratically in the difference in M number. However, it leaves the total Dicke space, and in particular the target Dicke state, invariant. In contrast, local dephasing induces coupling outside the Dicke space, but with a rate that is at most linear in N .

Consider the evolution during the N/ 2 control pulses to realize either of the phasing gates U s or U w . Assuming Guassian bath statistics, the effective dephasing rate can be written as the overlap of the noise spectrum S ω ( ) and the filter function | f ( ω ) | 2 (see e.g. [47, 48]):

$$A ( T ) = \frac { 1 } { 2 \pi } \int _ { 0 } ^ { \infty } d \omega S ( \omega ) | f ( \omega ) | ^ { 2 }. \quad \ \ ( 2 5 ) \quad \ \ \ s p e c _ { \text{$\bullet$:$d$} }$$

For an initial system-bath state ρ (0) = ρ S (0) ⊗ ρ B (0) with the bath in thermal equilibrium ρ B (0) = ∏ k (1 -e -βω k ) e -βω b b k † k k at inverse temperature β ( k B ≡ 1), the noise spectrum is S w ( ) = 2 ( π n ω ( ) + 1 / 2) I ( ω , ) where I ( ω ) = ∑ k | g k | 2 δ ω ( -ω k ) is the boson spectral density, and n ω ( k ) = ( e βω k -1) -1 is the thermal occupation number in bath mode k . The filter function is obtained from the windowed Fourier transform f ( w ) = ∫ T 0 C t e ( ) iωt , where C t ( ) is the time-dependent control pulse sequence. In the present case C t ( ) is a unit sign function that flips every time a collective spin flip is applied:

$$\underset { \text{$re$} } { \text{$re$} } \\ \text{$n$} _ { \text{$t} } = \begin{cases} \begin{array} { c c } 1 & t \in \cup _ { k = 1 } ^ { N / 2 } \{ [ T _ { k } ^ { ( 0 ) }, T _ { k } ^ { ( 1 ) } ) \cup [ T _ { k } ^ { ( 2 ) }, T _ { k } ^ { ( 3 ) } ) \} \\ - 1 & t \in \cup _ { k = 1 } ^ { N / 2 } \{ [ T _ { k } ^ { ( 1 ) }, T _ { k } ^ { ( 2 ) } ) \cup [ T _ { k } ^ { ( 3 ) }, T _ { k } ^ { ( 4 ) } ) \} \\ 0 & \text{otherwise} \end{array} \end{cases}$$

where T ( m ) k = mθ /g k + 4 ∑ k -1 j =1 θ /g j are the flip times with the duration between pulses growing linearly. The angles are θ k = 2 πk N +1 and the total time is

$$T = T _ { N / 2 } ^ { ( 4 ) } = \frac { \pi N ( N + 2 ) } { g ( N + 1 ) }.$$

The explicit form of the filter function is

$$\begin{array} { c c c } | f ( \omega ) | ^ { 2 } & = & \frac { 1 } { \omega ^ { 2 } } \Big | \sum _ { k = 1 } ^ { N / 2 } ( e ^ { i \omega T _ { k } ^ { ( 0 ) } } - 2 e ^ { i \omega T _ { k } ^ { ( 1 ) } } + 2 e ^ { i \omega T _ { k } ^ { ( 2 ) } } \\ & & - 2 e ^ { i \omega T _ { k } ^ { ( 3 ) } } + e ^ { i \omega T _ { k } ^ { ( 4 ) } } ) \Big | ^ { 2 }. \end{array}$$

In comparison, consider evolution where no spin flips are applied during the sequence, in which case the bare functions are C (0) ( ) t = 1 ∀ t ∈ [0 , T ), and | f (0) ( ω ) | 2 = 4 sin 2 ( Tω/ 2) /ω 2 . Results are plotted in Fig. 3 in the main paper and show there is substantial decoupling from the dephasing environment when the spectral density has dominant support in the range ω &lt; g/ 2. For 2 πkω/g glyph[lessmuch] 1, the summands in f ( ω ) can be expanded in a Taylor series in ω/g and to lowest order we find

$$g ^ { 2 } | f ( \omega ) | ^ { 2 } \approx \frac { ( \ome$$

This approximation is valid for ω/g &lt; 1 /πN , and, as shown in Fig. 3 in the main paper, for 1 /πN &lt; ω/g &lt; 1 / 2 the function is essentially flat with an average value g 2 | f ( ω ) | 2 ≈ 3 independent of N . In the region 1 /πN &lt; ω/g &lt; 1 / 2 the bare filter function is oscillatory and has an average g 2 | f (0) ( ω ) | 2 ≈ 13 63, while . for ω/g &lt; 1 /πN it asymptotes to π N 2 2 ( N +2) 2 ( N +1) 2 . Thus, in the region ω/g &lt; 1 /πN the ratio determining the reduction factor in the dephasing rate is | f ( ω ) | 2 | f (0) ( ω ) | 2 = π ω /g 2 2 2 , while for ω/g ∈ [1 /πN, 1 / 2], the reduction factor can be approximated by | f ( ω ) | 2 | f (0) ( ω ) | 2 ≈ 0 22, provided the noise . spectrum is sufficiently flat there. The effectiveness of this decoupling on precision and state precision fidelity is shown in Figure 4.

Further, the aforementioned freedom to apply the GPGs in any order allows room for further improvement. For example, consider coupling to a zero temperature Ohmic bath with noise spectrum S ω ( ) = αωe -ω/ω c and

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01120v5-figure-16255f42-fdb9-4c1a-a86b-f3aab97d20b7.jpg)

g

FIG. 4. a) Precision obtained for 70 spins with a single measurement of J z 2 as a function of mode decay for various global dephasing factors A T ( ): no global dephasing (blue line), A T ( ) = 10 -6 (blue dashed), A T ( ) = 10 -5 (blue dotdashed), A T ( ) = 10 -4 (blue dotted). These dephasings correspond to an underlying decoherence rate of γ gdp = 10 -4 g accumulated over each phasing gate of duration T . For a zero temperature Ohmic bath, the corresponding cuttoff frequencies are: ω /g c = { 0 003 0 022 0 094 . , . , . } . Black line: performance without dynamical decoupling with A 0 ( T ) = 0 0223. . e.g. if one were to switch the sign of the dispersive coupling during each GPG rather than flipping the spins. (b) Fidelity error for the same environments as above.

having cutoff frequency ω /g c = 0 1. . For N = 20, the ratio of the effective decay rate for the linearly ordered sequence of GPGs above to that with no decoupling is A T /A ( ) 0 ( T ) = 0 0085. . However, by sampling over permutations of the ordering of GPGs we find a sequence [49] achieving A T /A ( ) 0 ( T ) = 0 0026. .

## Preparation of error-tolerant states

The state preparation method we have described so far has some inherent tolerance to decoherence. However, once the state is prepared, further errors could accumulate such as qubit loss or dephasing, while waiting for the accumulation of the measurement signal. Some strategies to address this were recently proposed in Ref. [44] where they suggest using superpositions of Dicke states as probe states. The class of states considered there are

$$| \varphi _ { u } \rangle = \frac { 1 } { \sqrt { 2 ^ { n$$

Here the number of spins N = k × n × u , and the parameters u and n determine the robustness of the states to some number of loss and dephasing errors respectively, while k is a parameter to scale the number of qubits in the superposition (larger k means better performance). The case u = 1 tolerates erasure errors; specifically, the state | ϕ 1 〉 has a large quantum Fisher information obeying Heisenberg scaling when the number of erasure errors is less than n . We will consider the state performing well in the presence of one erasure error: u = 1 , n = 2 , k = N/ 2, which can be written

$$| \varphi _ { 1 } \rangle = \frac { 1 } { 2 } ( | J, - J$$

The case u = 2 tolerates a constant number of dephasing errors. We will focus on the state with u = 2 , n = 1 , k = N/ 2 which tolerates one dephasing error and can be written

$$| \varphi _ { 2 } \rangle = \frac { 1 } { \sqrt { 2 } } ($$

We now describe how to make these states using our protocol. A key ingredient to prepare a superposition of Dicke states is to perform a controlled state preparation. If we introduce an ancilla spin A which can be allowed to couple to the mode when the other spins do not (e.g. by detuning the other spins far away from the dispersive coupling regime), then a controlled displacements of the mode can be done:

$$\underset { \begin{array} {$$

Here R π ( | 1 〉 A 〈 1 ) | = e iπa a † | 1 〉 A 〈 1 | , meaning only the ancilla state | 1 〉 A couples to the mode. Now replacing the displacements D β ( ) and D ( -β ) with the controlled displacements Λ( β ) and Λ( -β ), the effect is a controlled GPG (see also Ref [45]):

$$\Lambda ( U _ { G P G } ) = | 0 \rangle _ { A } \langle 0 | \otimes \mathbf 1 + | 1 \rangle _ { A } \langle 1 | \otimes U _ { G P G }. \quad ( 3 0 )$$

Thus by simply replacing all instances of GPGs with controlled GPS we can achieve a controlled Grover step unitary G . Note the unitary e iJ y π/ 2 conjugating W (0) in U s does not need to be controlled, meaning the entire unitary U # G G can be made into a controlled unitary

$$\cdot \quad \Lambda ( U _ { G } ^ { \# G } ) = ( | 0 \rangle _ { A } \langle 0 | \otimes \mathbf 1 + | 1 \rangle _ { A } \langle 1 | \otimes U _ { G } ) ^ { \# G }. \quad ( 3 1 )$$

This is not quite enough. The state preparation of a Dicke states described above applies U # G G to a particular

initial state, namely the spin coherent state | s 〉 . We will also require a way to perform a controlled rotation on all the spins of the form

$$\Lambda ( e ^ { i J ^ { y } \pi / 2 } ) = | 0 \rangle _ { A } \langle 0 | \otimes \mathbf 1 + | 1 \rangle _ { A } \langle 1 | \otimes e ^ { i J ^ { y } \pi / 2 }. \quad ( 3 2 ) \quad \ 3.$$

Without having direction interactions between the ancilla and the other spins it is not obvious how to do this. However, it is possible to mediate the interaction with the mode by choosing φ = 0 and θ glyph[lessmuch] 1 in one instance of a controlled GPG. This will give Λ( U GPG ( θ, 0 , π/ 4 ) θ ≈ | 0 〉 A 〈 0 | ⊗ 1 + 1 | 〉 A 〈 1 | ⊗ e -iJ z π/ 2 where we have approximated sin( θJ z ) ≈ θJ z . Note, in order for this to be valid we require θ glyph[lessmuch] 1 /N and consequently χ = | α | 2 glyph[greatermuch] N , i.e. the area of the GPG in phase space needs to grow with N , or the gate could be composed into N GPGs each of area of O (1). This will consequently incur a loss of fidelity due to mode decay

$$F _ { \text{pro} } ( \mathcal { E }, U _ { \text{GPG} } ) > e ^ { - 6 \pi | \alpha | ^ { 2 } \kappa / g } \cos ( | \alpha | ^ { 2 } 4 \pi \kappa / g ). \quad ( 3 3 ) \quad \text{To}$$

but no worse than the performance for state preparation without ancilla.

The controlled operation is then

$$\Lambda ( e ^ { i J ^ { \nu } \pi / 2 } ) = e ^ { - i J ^ { \tilde { \tau } } \pi / 2 } \Lambda ( U _ { G P G } ( \theta, 0, \pi / 4 \theta ) e ^ { i J ^ { \tilde { \tau } } \pi / 2 }. \quad \ t$$

We can now write the process to prepare the state | ϕ 2 〉 :

- 1. Prepare the product state 1 √ 2 ( 0 | 〉 + 1 ) | 〉 A ⊗| J, - 〉 J .
- 2. Apply e -iJ x π/ 2 Λ( U GPG ( θ, 0 , π/ 4 ) θ e iJ x π/ 2 .
- 3. Apply Λ( U # G G ). This involves N × # G instances of Λ( U GPG ) for varying parameters.
- 4. Measure the ancilla in the |± 〉 x A basis. The outcomes r = ± 1 each occur with probability 1 / 2. The conditional system state is

$$\frac { 1 } { \sqrt { 2 } } ( | J, - J \rangle \pm | J, 0 \rangle )$$

- 5. Apply the classically controlled product unitary Z r ( ) = Z (1 -r / ) 2 ⊗ N . If we assume N/ 2 is odd then

$$Z ( r ) \frac { 1 } { \sqrt { 2 } } ( | J, - J \rangle \pm | J, 0 \rangle ) = | \varphi _ { 2 } \rangle.$$

To prepare the state | ϕ 1 〉 a similar process can be used. However, rather than the | 0 〉 A state being correlated with the product state | J, - 〉 J we want it correlated with the GHZ state 1 √ 2 ( | J, - 〉 J + | J, J 〉 ). Such a state can be prepared using one additional controlled GPG gate. This follows from the observation that e i J y π 2 U GPG ( π, π/ 2 , π 8 ) e -i J y π 2 | J, - 〉 J = 1 √ 2 ( | J, - 〉 J + | J, J 〉 ). These two processes are summarized by the following circuits:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01120v5-figure-402d5c67-ade0-4afb-9f4b-2f973fa7d7ca.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01120v5-figure-8641f9e3-b6b8-4d78-935c-9bf5b2ce54e5.jpg)

- [7] R. Demkowicz-Dobrza´ski, J. Koglyph[suppress] lody´ski n n and M. Gut ¸˘ a Nat. Comms. 3 , 1063 EP (2012).
- ∗ gavin.brennen@mq.edu.au
- [1] W. Wasilewski, K. Jensen, H. Krauter, J. J. Renema, M. V. Balabas, and E. S. Polzik, Phys. Rev. Lett. 104 , 133601 (2010).
- [2] J. Taylor et al. , Nat. Phys. 4 , 810 (2008).
- [3] W. Wang et al. , Nat. Comm. 10 , 4382 (2019).
- [4] V. Giovannetti, S. Lloyd and L. Maccone, Nat. Photonics 5 , 222 (2011).
- [5] V. Giovannetti, S. Lloyd and L. Maccone, Science 306 , 1330 (2004).
- [6] S. Pirandola, B. R. Bardhan, T. Gehring, C. Weedbrook, and S. Lloyd, Nat. Photonics 12 , 724 (2018).
- [8] M. Zwierz, C. A. P´rez-Delgado and P. Kok, Phys. Rev. e Lett. 105 , 180402 (2010).
- [9] S. Zhou, M. Zhang and L. Jiang, Nat. Comms. 9 , 78 (2018).
- [10] S. D. Bartlett, G. K. Brennen and A. Miyake, Quantum Science and Technology 3 , 014010 (2017).
- [11] X. Wang and P. Zanardi, Phys. Rev. A 65 , 032327 (2002).
- [12] G. Xu and G. Long, Phys. Rev. A 90 , 022323 (2014).
- [13] A. Sørensen and K. Mølmer, Phys. Rev. Lett. 82 , 1971 (1999).
- [14] I. Apellaniz, B. L¨cke, J. Peise, C. Klempt, and G. T´, u o New J. Phys. 17 , 083027 (2015).
- [15] This follows from U η ( ) = exp[ iη ( J x sin δ + J y cos δ )] =

- exp( iδJ z ) exp( iηJ y ) exp( -iδJ z ), the fact our measurement observable is ˆ J z 2 , and our initial state is | J, M = 0 . 〉
- [16] A. B¨ artschi and S. Eidenbenz, Preprint at https://arxiv.org/abs/1904.07358 (2019).
- [17] C. Wu et al. , Phys. Rev. A 95 , 013845 (2017).
- [18] K. D. B. Higgins et al. , Nat. Comms. 5 , 4705 EP (2014).
- [19] T. Keating, C. H. Baldwin, Y.-Y. Jau, J. Lee, G. W. Biedermann and I. H. Deutsch, Phys. Rev. Lett. 117 213601 (2016).
- [20] D. B. Hume, C. W. Chou, T. Rosenband, and D. J. Wineland, Phys. Rev. A 80 052302 (2009).
- [21] L. Jiang et al. , Nat. Phys. 4 , 482 EP (2008).
- [22] This was tried for up to N = 100 spins with the decomposition having fidelity error 1 -|〈 λ K J, | ˜ | - 〉| J 2 &lt; 10 -12 for synthesizing sampled random state vectors | λ 〉 .
- [23] D. J. Rowe, H. de Guise and B. C. Sanders, J. Math. Phys. 42 , 2315 (2001).
- [24] B. L¨cke u et al. , Science 334 , 773 (2011).
- [25] G. K. Brennen, K. Hammerer, L. Jiang, M. D. Lukin, and P. Zoller, Preprint at https://arxiv.org/abs/0901.3920 (2009).
- [26] See Supplemental Material [url] for details, which includes Ref. [45].
- [27] Brennen, G. K., Pupillo, G., Rico, E., Stace, T. M. &amp; Vodola, D. Phys. Rev. Lett. 117 , 240504 (2016).
- [28] M. Nielsen and I. L. Chuang Quantum Computation and Quantum Information (Cambridge University Press, 2010).
- [29] G. S. Agarwal Physica Scripta 82 , 038103 (2010).
- [30] Z.-Y. Wang and R.-B. Liu Quantum Error Correction pp. 351-375 (Cambridge University Press, 2013).
- [49] The ordering { 8 4 5 9 3 7 6 2 10 1 , , , , , , , , , } achieves this. Note an exhaustive search over all 10! permutations was not done.
- [32] Y. Ouyang, N. Shettell and D. Markham, Preprint at
- https://arxiv.org/abs/1908.02378 (2019).
- [33] H. Pollatsek and M. B. Ruskai, Linear Algebra and Applications, 392 , 255 (2004).
- [34] J.A. Gross, arXiv:2005.10910.
- [35] C. Sayrin et al. , Nature 477 , 73 EP (2011).
- [36] S. Garcia, M. Stammeier, J. Deiglmayr, F. Merkt, and A. Wallraff, Phys. Rev. Lett. 123 , 193201 (2019).
- [37] J. S. Pedernales et al. , Sci. Rep. 5 , 15472 EP (2015).
- [38] M. Lee et al. , Phys. Rev. Lett. 122 , 153603 (2019).
- [39] C. Wang et al. , Science 352 , 1087 (2016).
- [40] T. Astner et al. Phys. Rev. Lett. 118 , 140502 (2017).
- [41] S. J. Bosman, V. Singh, A. Bruno, and G. A. Steele, Appl. Phys. Lett. 107 , 192602 (2015).
- [42] A.R. Castellia, L.A. Martinez, J.M. Pate, R.Y. Chiao, and J.E. Sharping, J. Appl. Phys. 126 , 173908 (2019).
- [43] Brennen, G. K., Hammerer, K., Jiang, L., Lukin, M. D. &amp; Zoller, P. Global operations for protected quantum memories in atomic spin lattices. Preprint at https://arxiv.org/abs/0901.3920 (2009).
- [44] Ouyang, Y., Shettell, N. &amp; Markham, D. Robust quantum metrology with explicit symmetric states. Preprint at https://arxiv.org/abs/1908.02378 (2019).
- [45] Brennen, G. K., Pupillo, G., Rico, E., Stace, T. M. &amp; Vodola, D. Loops and strings in a superconducting lattice gauge simulator. Phys. Rev. Lett. 117 , 240504 (2016).
- [46] Pollatsek, H. &amp; Ruskai, M. B. Permutationally invariant codes for quantum error correction. Linear Algebra and Applications , 392 , 255 (2004).
- [47] Agarwal, G. S. Saving entanglement via a nonuniform sequence of π pulses. Physica Scripta 82 , 038103 (2010).
- [48] Wang, Z.-Y. &amp; Liu, R.-B. Quantum Error Correction pp. 351-375 (Cambridge University Press, 2013).
- [49] The ordering { 8 4 5 9 3 7 6 2 10 1 , , , , , , , , , } achieves this. Note an exhaustive search over all 10! permutations was not done.
Mattias T. Johnsson Nabomita Roy Mukty Daniel Burgarth Thomas Volz Gavin K. Brennen
Geometric pathway to scalable quantum sensing
Entangled resources enable quantum sensing that achieves Heisenberg scaling,
a quadratic improvement on the standard quantum limit, but preparing large
scale entangled states is challenging in the presence of decoherence. We
present a quantum control strategy using highly nonlinear geometric phase gates
for preparing entangled states on spin ensembles which can be used for
practical precision metrology. The method uses a dispersive coupling of $N$
spins to a common bosonic mode and does not require addressability, special
detunings, or interactions between the spins. Using a control sequence that
executes Grover's algorithm on a subspace of permutationally symmetric states,
a target entangled resource state can be prepared using $O(N^{5/4})$ geometric
phase gates. The geometrically closed path of the control operations ensures
the gates are insensitive to the initial state of the mode and the sequence has
built-in dynamical decoupling providing resilience to dephasing errors.