## Demonstration of Entanglement-Enhanced Phase Estimation in Solid

Gang-Qin Liu, 1 Yu-Ran Zhang, 1 Yan-Chun Chang, 1 Jie-Dong Yue, 1 Heng Fan, 1,2, ∗ and Xin-Yu Pan 1,2, †

1 Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 2 Collaborative Innovation Center of Quantum Matter, Beijing 100190, China (Dated: April 9, 2015)

Precise parameter estimation plays a central role in science and technology. The statistical error in estimation can be decreased by repeating measurement, leading to that the resultant uncertainty of the estimated parameter is proportional to the square root of the number of repetitions in accordance with the central limit theorem. Quantum parameter estimation, an emerging field of quantum technology, aims to use quantum resources to yield higher statistical precision than classical approaches. Here, we report the first room-temperature implementation of entanglement-enhanced phase estimation in a solid-state system: the nitrogen-vacancy centre in pure diamond. We demonstrate a super-resolving phase measurement with two entangled qubits of different physical realizations: an nitrogen-vacancy centre electron spin and a proximal 13 Cnuclear spin. The experimental data shows clearly the uncertainty reduction when entanglement resource is used, confirming the theoretical expectation. Our results represent an elemental demonstration of enhancement of quantum metrology against classical procedure.

Information about the world is acquired by observation and measurement, the results of which are subject to error [1]. The classical approach to reduce the statistical error is to increase the number of resources for the measurement in accordance with the central limit theorem, however, this method sometimes seems undesirable and inefficient [2]. Quantum parameter estimation, the emerging field of quantum technology, aims to yield higher statistical precision of unknown parameters by harnessing entanglement and other quantum resources than purely classical approaches [3]. Since this quantum-enhanced measurement will benefit all quantitative science and technology, it has attracted a lot of attention as well as contention. Using N independent particles to estimate a parameter ϕ can achieve at best the standard quantum limit (SQL) or called shot noise limit scaling as δϕ ∝ 1 / √ N while it is believed that using N entangled particles and exotic states such as NOON states in principle is able to achieve the inviolable Heisenberg limit scaling as δϕ ∝ 1 /N [4, 5]. In such circumstances, there are many efforts using non-classical states and quantum strategy for sub-SQL phase estimation in different physical realizations, such as optical interferometry [2, 6-9], atomic systems [10, 11], and Bose-Einstein condensates [12, 13].

In this paper, we report the first room-temperature proofof-principle implementation of entanglement-enhanced phase estimation in a solid-state system: the nitrogen-vacancy (NV) centre in pure diamond single crystal. An individual NV center can be viewed as a basic unit of a quantum computer in which the nuclear spin with a long coherence time performs as the memory and the centre electron spin with a high control speed acts as the probe. This solid-state system is one of the most promising candidates for quantum information processing (QIP), and many coherent control and manipulation processes have been performed with this system [14-29]. Here, we demonstrate a super-resolving phase measurement with two entangled qubits of different physical realizations: a NV centre electron spin and a proximal 13 C nuclear spin. We are able to improve the phase sensitivity by factors close to √ 2 compared with the classical scheme, which conforms to the fundamental Heisenberg limit. As we have entangled two qubits with different physical realizations, our results represent a more generalized and elemental demonstration of enhancement of quantum metrology. Moreover, our system has overcomed the defects of post-selection in the most common optical systems which are fatal due to the fact that the measurement trials abandoned will eliminate the quantum advantage over classical strategy.

## Results

## System description

The phase estimation scheme is implemented by optically detected magnetic resonance (ODMR) [14, 16] technique on

FIG. 1: General scheme and system description. (a) Phase estimation schemes of the independent states and the electron-nuclear entangled state. By harnessing entanglement, quantum metrology yields higher statistical precision than classical approaches. (b) Energy levels and physical encoding of the two-bit system. The electron spin and a nearby 13 C nuclear spin of an NV center are employed to demonstrate the metrology scheme. At excited-state level anti-crossing (ESLAC), both spins can be polarized, manipulated and readout with high fidelity. Two-bit conditional quantum gates are implemented by applying selective microwave (MW) or (RF) pulses.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-f788751d-4d86-46cf-8aa6-52ac53c2c4af.jpg)

a home-built confocal microscope system. The description of the system can be found in [27]. The spin-1 electron spin of NV centre has triplet ground states with a zero-field splitting of ∆ ≈ 2 87 . GHz between the states | 0 〉 and | ± 1 〉 . As an external magnetic field of about 507 Gauss is applied along [111] direction of the diamond crystal, the degeneration of | ± 1 〉 states can be well relieved, and the first qubit is encoded on the | 0 〉 and | -1 〉 subspace ((hereinafter labelled as | 1 〉 )). The electron spin state can be initialized to | 0 〉 state by a short 532 nm laser pulse (3 µ s) and manipulated by resonant microwave (MW) pulses of tunable duration and phase. The electron spin state is readout by collecting the spin-dependent fluorescence. To enhance the fluorescence collection efficiency, a solid immersion lens (SIL) [30] is etched above the selected NV center, typical count rate in this experiment is 250 kps with SIL, see Methods for details.

The second qubit is encoded on the | ↑〉 and | ↓〉 states of a nearby 13 C nuclear spin. See Fig. 1(b) for the energy levels of the two qubit system. The coupling strength between the target nuclear spin and centre electron spin is 12.8 MHz, which indicates the 13 C atom sites on the third shell from the NVcentre [31]. The polarization and readout procedure of the nuclear spin is more complicated than that of a electron spin. The 507 Gauss magnetic field causes excited-state level anticrossing (ESLAC) of centre electron spin, in which the optical spin polarization of centre electron will transfer to nearby nuclear spins [32, 33]. So the host 14 N nuclear spin, the nearby 13 C nuclear spin as well as the center electron spin are polarized by the same laser pulse under this magnetic field. To readout the nuclear spin state, a mapping gate, which transfers nuclear spin state to electron spin, and a following optical readout of electron spin state are employed [17, 33], see Methods for details.

The nearby nuclear spin couples to the centre electron spin through strong dipolar interaction, which provides excellent conditions to implement two-qubit controlled gate. On the one hand, the resonant frequency of | 0 ↑〉 ⇔ | 1 ↑〉 transition and | 0 ↓〉 ⇔ | 1 ↓〉 transition are separated by 12.8 MHz from each other, so we can selectively manipulation one branch of nuclear spin with high fidelity while keep the other branch untouched (using weak MW pulses, see black arrow in Fig. 1(b)). On the other hand, the nuclear spin state evolution is strong affected by the state of electron spin: when electron spin is on the | 0 〉 state (or | -1 〉 state), the dynamics of the nuclear spin is dominated by the external magnetic filed (or the dipolar interaction, respectively), its Zeeman splitting between the | ↑〉 and | ↓〉 states is about 500 kHz, which is far away from the dipolar interaction strength of 12.8 MHz. Therefore, we can selectively manipulate nuclear spin state in one branch of electron spin, as well (using RF pulses, see red arrow in Fig. 1(b) and Supplementary Note 1).

## Phase preparation and measurement

Fig. 2(b) describes the pulse sequence to prepare and mea-

FIG. 2: Phase preparation and measurement (a) For each unknown state sitting on the equatorial plane of the Bloch sphere, we compare the amplitudes of the two Rabi signals from this state (under orthogonal microwave pulses driven) to extract the original phase information: ϕ = arctan( A /A x y ) . (b) Upper pane: pulse sequence to prepare and measure electron-nuclear spin entanglement state 1 √ 2 ( 1 | ↑〉 + e 2 iϕ | 0 ↓〉 ) . Lower pane: pulse sequence to prepare and measure nuclear spin superposition state 1 √ 2 ( 0 | ↑〉 + e iϕ | 0 ↓〉 ) . (c-e) Nuclear Rabi signals for phase measurement, with 0.1 M, 0.4 M and 2 M repetition of pulse sequence in (b) . The input phase ϕ is 30 ◦ . Solid circle with blue fitting line is driven by 0 ◦ RF pulse (X measurement), and square with red dash fitting line is driven by 90 ◦ RF pulse (Y measurement). (f) Dependence of measured phase and its standard deviation on repeat number.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-2c408260-3f73-42e9-bc72-1de9bc3bd262.jpg)

sure the phase of a superposition state. Take nuclear spin for example, the qubit is defined in a rotating frame with frequency equalling to the energy splitting between | 0 ↑〉 and | 0 ↓〉 states. After polarized to | 0 ↑〉 by laser pulse, a resonant RF π 2 pulse brings the system to 1 √ 2 ( 0 | ↑〉 + e iϕ | 0 ↓〉 ) state. The phase of this state is determined by the relative phase of the applied RF pulse, which is tunable in experiment. The phase of electron superposition state is prepared in the same way, with resonant MW pulses.

The phase information of a superposition state is detected by converting it to population information of the spin qubits and a following optical readout. To eliminate the system error in long time measurement, we use a self-calibration measurement scheme as shown in Fig. 2(a). For each unknown state siting on the equatorial plane of the Bloch sphere, we measure the Rabi oscillations driven by two orthogonal microwave pulses (0 ◦ and 90 ◦ ), and compare the amplitudes of the two Rabi signal to extract the original phase information: ϕ = arctan A x A y , where A x and A y are the amplitudes of Rabi signals driven by 0 ◦ and 90 ◦ microwave pulses, respectively. Note that this is a single spin experiment, and we need to repeat the pulse sequence many times to get a reliable signal to noise ratio (SNR). Figs. 2(c-e) presents the nuclear Rabi signal of 0.1 M, 0.4 M and 2 M repeat of the pulse sequence in

FIG. 3: State tomography and phase relation. (a) The state tomography result of a nuclear spin superposition state 1 √ 2 ( | ↑〉 + | ↓〉 ) . (b) State tomography result of the electron-nuclear entangled state 1 √ 2 ( 1 | ↑〉 + 0 | ↓〉 ) . (c) Phase relation of an independent state. (d) Phase relation of an entangled state. Compared with independent state, the phase relation of entangled state has double frequency dependence on input phase, so a more precise phase estimation result can be achieved with entangled state.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-9a1b7731-0f2f-4fce-b1f1-c5239d3b9d78.jpg)

Fig. 2(b), with an input phase of 30 ◦ . The error bar of the data point represent standard deviation (SD) of 10 repeat measurements. It is clear to see from Fig. 2(c) to Fig. 2(e) that the SNR is better as the measurement sequence is repeated more times. In Fig. 2(f), the decrease of phase estimation error can be well described by central limit theorem.

To improve phase estimation accuracy, one can increase the repeat number, which means longer measurement time is needed. An equivalent way is to employ more qubits. As mentioned before, the state of the multi-qubit system, independent or entanglement, determines the accuracy limit of phase estimation. For the investigated two-qubit system, the electron and nuclear spin can be prepared and measured independently. Fig. 3(a) plots the state tomography result of a nuclear spin superposition state. Using such independent state (either nuclear spin or electron spin) will get a phase relation as depicted in Fig. 3(c), the amplitude of Rabi signal has cosine dependence on the phase of input state.

The electron and nuclear spin can be prepared in entangled state by combination of MW and RF pulses. As shown in the upper pane of Fig. 2(b), after the first RF π 2 pulse (with phase ϕ ) brings the system to 1 √ 2 ( 0 | ↑〉 + e iϕ | 0 ↓〉 ) state, a selective MW π pulse of | 0 ↑〉 ⇔ | 1 ↑〉 transition, which has relative phase ϕ to the first RF pulse, brings the system to 1 √ 2 ( 1 | ↑〉 + e 2 iϕ | 0 ↓〉 ) state. The MW and RF channels are synchronized to the same clock reference and relative phase between them is calibrated before each measurement. Typical state tomography result of an electron-nuclear entangled state ( ϕ = 0 ◦ ) is depicted in Fig. 3(b). It is worth noting that

FIG. 4: Dependence of phase error on repeat number and input phases. (a) Standard deviations of the entangled state and independent state against repeat number. Single spin states (electron and nuclear) of the same input phase ϕ = 30 ◦ are prepared and measured independently, the phase extracted from each entangled state 'single measurement' is 2 ϕ = 60 ◦ . (b) The phase uncertainty is fitted by function δϕ = a/ √ ν + c for both single spin state and the entangled state, where ν corresponds to repeat number in unit of million (M). The curves show that the phase uncertainty, phase error represented by δϕ , of entanglement case is apparently lower than the case of single spin state. (c) The phase estimation results of different input phases including 0 ◦ , 5 ◦ , 10 ◦ , 15 ◦ , 20 ◦ , 25 ◦ and 30 ◦ , the repeat number is fixed to 1 M. (d) The phase error of different input phases, the repeat number is fixed to 1 M.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-df10a049-c086-4fe5-879c-3d3b277b7d3c.jpg)

the dephasing time of electron spin (0.7 µ s, see Methods) is very short compared with the typical manipulation time (e.g. 10 µ s for a flip operation) of nuclear spin, which limits the QIP applications of this entangled state [22, 27, 33]. However, in our phase estimation application, the sensitive phase information is converted to population right after its generation and then only limited by T 1 of electron spin, which is about 5 ms. Meanwhile, the coherence of electron is less affected under microwave driving (see Methods), thus the phase of the entangled state is well preserved during the preparation and measurement. As shown in Fig. 3(d), the phase relation of the entangled state has double frequency dependence on the phase of input state, so the phase estimation using the entangled state of two qubit is more precise than that of using two state from independent single qubit.

## Entanglement-enhanced phase estimation

To demonstrate the merits of entangled state over independent state in phase estimation application, we compare their performances on different repeat number and different input phase, the measured results are summarized in Fig. 4. The experimental procedure is: firstly, single spin states (electron and nuclear) of the same input phase ( 30 ◦ ) are prepared and measured independently. Then the output phases extracted

from the same repeat number are counted together, no weight is added for either electron or nuclear spin states. For a fair comparison with entangled state, half of the statistic samples ( ν ) are extracted from electron spin states, and the other half ( ν ) are extracted from nuclear spin states. In the case of entangled state phase estimation, the entangled states are prepared and measured using the same repeat number ( ν ). Note the same MW and RF channels are used to prepare the independent and entangled states.

As shown in Fig. 4(a) and (c), the phase extracted from each entangled state measurement is 2 ϕ , so the phase error of input phase is just half of the standard deviation from ν sample statistic ( δϕ/ 2 ). For the independent-state input, the double sample number ( 2 ν ) only suppresses phase error to δϕ/ √ 2 level, which is larger than entanglement-state input.

Explicitly, we would expect that the phase uncertainty δϕ proportion respectively to 1 / √ νN and 1 / ( √ νN ) for single spin state and the entangled state with N being 2 in our experiment since two-qubit entanglement between electron spin and nuclear spin is used as the quantum resource. In Fig. 4(b), we consider that identical measurement is repeated ν times, then a general formulae, δϕ = a/ √ ν + c , is used to fit experimental data, where c is assumed to be a systematic error depending on specific experimental setup. The parameter a for scheme of entangled state should be smaller than that of the single spin scheme, corresponding to smaller uncertainty about the phase, if we assume that the single spin state and the entangled state are realized by the same physical state. Fig. 4(b) of the experimental data demonstrates clearly that the precision of phase estimation is enhanced by using entanglement which agrees well with theoretical expectation. Here the discrepancy between theory and experiment is possibly due to two related reasons: the electron spin state and the nuclear spin state are not the same, in particular for their decoherence time, while their similarity is assumed theoretically; the readout of NV centre system can only be by intensity of florescence of the electron spin. Besides the data processing method presented here, we have also tried linear fitting in the log-log scale for SD as well as variance (see Supplementary Note 2 and Supplementary Fig. 3 for detail). The obtained results are in good agreement with the results in Fig. 4(b), which confirms the validity of our conclusion.

Fig. 4(c) and (d) show the phase estimation results of different input phases. The phase error of entangled state is smaller than independent state in all input phases, which indicates the enhancement of phase estimation accuracy by entanglement is phase independent.

## Discussion

As summarized in Fig. 4 by different figures of merit quantifying the uncertainty of phase estimation, the entanglementenhanced precision is clearly shown by experimental data. This experiment demonstrates the advantage of the quantum metrology scheme. Practically by using quantum metrology, the measured physical quantity should have the same interaction on the probe system no matter it is prepared as a single qubit or entangled state. In our special designed experiment, the measured phases are artificially encoded to the probe state such that the enhancement of precision can be shown by entangled probe state. However in principle, the confirmation of theoretical expectation by experimental data provides a solid evidence that quantum phase estimation is applicable in this solid state system.

In this experiment, we use repeating measurement to overcome the low photon collection efficiency of NV center. The phase estimation accuracy can be further improved by employing single-shot measurement technique, which is now available in NV system [19, 20, 34]. Although the photon collection efficiency is not perfect ( &lt; 20%, not every measurement is stored and counted), the following two facts guarantee the reliability of the demonstration: (1) we use the same scheme to measure single and entangled states, that is, the phase information is finally converted to fluorescence signal of NV center and detected. (2) The detection efficiency of the system is stable (though not perfect as single-shot readout) for all the measurement, so we can directly compare the measured phase noise of single and entangled states.

As the phase estimation accuracy is determined by the total number of entangled qubits, a straightforward way to improve the phase accuracy is increasing the involved spin number. The large amount of weakly coupled 13 C nuclear spins around NV center are one of the best candidates. With the assistance of dynamical decoupling on center electron spin, up to 6 13 C nuclear spins can be coherent manipulated [3537]. Multi-qubit application such as error correction has been demonstrated in this system [23, 24].

In conclusion, we report the first room-temperature implementation of entanglement-enhanced phase estimation in a solid-state system: the nitrogen-vacancy (NV) centre in pure diamond. We demonstrate a super-resolving phase measurement with two entangled qubits of different physical realizations: a NV centre electron spin and a proximal 13 C nuclear spin. Thus, our results represent a more generalized and elemental demonstration of enhancement of quantum metrology against classical procedure, which fully exploits the quantum nature of the system and probes.

## Methods

## Cram´ er-Rao bound and quantum Fisher information

In the simplest version of the typical quantum parameter estimation problem, we aim to recover the value of a unknown continuous parameter [say phase ϕ in Fig. 1(a)] encoded in a fixed set of states ρ ϕ of a quantum system [3]. We can obtain a single result ξ via performing a measurement on the system and it is useful to express the measurement in terms of set of POVM { ˆ E ξ } . With large number of measurements, it is possible to calculate the estimator ˜( ϕ ξ ) with the observation

conditional probability density function ( pdf ) of result ξ given the true values ϕ : p ξ ϕ ( | ) = Tr ( ˆ E ρ ξ ϕ ) . When the number of measurements ν is sufficiently large and the estimation is unbiased [38, 39], the root-mean square error for the statistical uncertainty can be shown to obey the well-known Cram´ er-Rao bound [40] given by

$$\delta \varphi \equiv \sqrt { \sum _ { \xi } [ \tilde { \varphi } ( \xi ) - \varphi ] ^ { 2 } p ( \xi | \varphi ) } \geq \frac { 1 } { \sqrt { \nu F ( \varphi ) } }, \quad ( 1 ) \quad \begin{array} { c } \frac { 1 } { \sigma } \\ \frac { 1 } { \sigma } \end{array} \right ) \text{where } F ( \varphi ) \equiv \sum _ { \varepsilon } p ( \xi | \varphi ) \left [ \partial _ { \varphi } \ln p ( \xi | \varphi ) \right ] ^ { 2 } \text{ is the Fisher infor-} \quad \begin{array} { c } \frac { 1 } { \sigma } \end{array}$$

where F ϕ ( ) ≡ ∑ ξ p ξ ϕ ( | ) [ ∂ ϕ ln p ξ ϕ ( | )] 2 is the Fisher information corresponding to the selected POVM and the conditional pdf of the result. Eq. (1) provides a lower bound for the achievable lower bound by choosing the optimal measurement expressed by some POVM { ˆ E opt ξ } that maximizes the Fisher information: F Q ( ϕ ) = max { ˆ E opt ξ } F ϕ ( ) which is known as the quantum Fisher information (QFI) [41].

For the classical scheme with separable probe state [( 0 | 〉 + e -iϕ | 1 ) 〉 / √ 2] ⊗ N , the lower bound at best leads to the SQL δϕ se ∝ 1 / √ νN . To implement the quantum counterpart of the Heisenberg limit δϕ en ∝ 1 / √ νN , we can choose the GHZ state ( 0 | 〉 ⊗ N + e -iNϕ | 1 〉 ⊗ N ) / √ 2 as the optimal probe state. Consider that the qubit number N = 2 , the two-qubit maximally entangled state will obtain a √ 2 advantage against the separable state.

## Sample preparation

High purity single crystal diamond (Element Six, N concentration &lt; 5 ppb) is used for this experiment. There is almost no natural NV center in this diamond. NV centres are produced by electron implantation (7.5 Mev) and a following 2 hours vacuum annealing (at 800 ◦ C). Due to the random distribution of 13 C nuclear spins, the spin bath of individual NV center can be very different [31]. We choose NV centers with nearby 13 C nuclear spins, which can be identified by the extra splitting in ODMR signal, to implement the two-bit metrology scheme. Fig. 5(a) presents the physical structure of an NV center and a nearby 13 C nuclear spin. Fig. 5(c) is ODMR signal of this two-bit system. The coupling strength between electron spin and the selected 13 C nuclear is 12.8 MHz. Fig. 5(b) shows two dimensional fluorescence image of the FIB-etched SIL. The cross-cursor marked bright spot (blue) is the one used for this experiment.

## Coherence of electron spin and nearby nuclear spin

In this pure diamond, the coherence of NV electron spin is dominated by the randomly distributed 13 Cnuclear spins (natural abundance, 1.1%). The dephasing time of individual NV centres can be significantly different [42], from less than 1 µ s to nearly 10 µ s. From the free-induction decay signal of this NV center in Fig. 5(d), we extract the dephasing time ( T ∗ 2 ) of this electron spin, which is 0.72 µ s and much larger than the

FIG. 5: Characterization of the two-qubit system. (a) NV centre with a nearby 13 C nuclear spin. (b) 2D fluorescence scan of the FIB-etched solid immersion lens (SIL, 12 µ m diameter). The bright spot is the investigated NV centre. (c) ODMR of NV electron spin under magnetic field of 690 Gauss. The splitting is caused by the nearby 13 C nuclear spin. (d) Rabi oscillation and Free-induction decay (FID) of electron spin ( B = 507 Gauss). Due to the thermal fluctuation of the spin bath, centre electron spin picks up random phase during free precession and loses coherence. With the help of resonant MW pulses, electron spin can be flipped in short time, and is less affected by the bath noise.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-789118e9-1901-4d73-927b-17afa4721841.jpg)

time consumption of single manipulation on it (about 70 ns, see Rabi oscillation of electron spin in the same figure). It is worth noting that the dephasing time is not the direct limitation of electron manipulation duration. The latter is usually named T 1 ρ and can be characterized by the envelop decay time of electron spin Rabi oscillation [43]. The dephasing time of nuclear spin ( T ∗ 2 n = 270 µ s) is much longer than that of electron spin. Meanwhile, the used half π pulse of nuclear spin is only several microseconds, so the dephasing effect of nuclear spin can be ignored, See Supplementary Fig. 1 for the FID signal of the nearby nuclear spin.

## Coherent manipulation of electron spin and nuclear spin at ESLAC

As mentioned in the main text, we work at the excitedstate level anti-crossing (ESLAC) point to achieve fast and high fidelity initialization of the electron-nuclear two-qubit system. Under an external magnetic field of 507 Gauss (along the quantization axis of the selected NV) and laser excitation (532 nm), the electron and nuclear spins are polarized simultaneously. Fig. 6(a) shows ODMR spectrum of centre electron spin at such magnetic field. From the contrast difference of two peaks, which correspond to | ↑〉 and | ↓〉 states of 13 C nuclear spin, we estimate the polarization rate of this nuclear spin is about 85% (in | ↑〉 state). Furthermore, by measuring the pulse-ODMR spectrum of electron spin, we conclude that

the host 14 N nuclear spin is completely polarized under this magnetic field.

Fig. 6(b) shows the pulse sequence of electron and nuclear spin manipulation. After polarization with high fidelity, both spin states can be manipulated with resonant MW (or RF) pulses. For electron spin, the final state is readout by counting the fluorescence intensity of NV center, since | 0 〉 state is brighter than | 1 〉 state. For nuclear spin state, we use a mapping gate, which is composed by a weak pulse of | 0 ↑〉 ⇔ | 1 ↑〉 transition, to transfer the its state to electron spin and then readout optically. For example, an unknown nuclear spin state of | 0 〉 ⊗ ( α | ↑〉 + β | ↓〉 ) is transferred to α | 1 ↑〉 + β | 0 ↓〉 after applying the mapping gate. We carefully tuned the microwave power and pulse duration to maximum the flip efficiency while avoiding the unwanted non-resonant excitation. By comparing the Rabi amplitude of nuclear spin (Fig. 6(d), with mapping gate) and electron spin (Fig. 5(d), without mapping gate), we conclude that the mapping gate has transfer efficiency of more than 92%.

Figs. 6(c) and 6(d) present the pulse-ODMR spectrum and Rabi oscillation of the nearby nuclear spin when electron spin is at | 0 〉 state. The resonant frequency of this nuclear spin is 495 kHz, which is smaller than the Larmor frequency of 13 C nuclear spin under this magnetic field (542 kHz). We attribute this modification to the 'enhance effect' of center electron spin. As the nuclear spin is close to the electron spin, the nonsecular terms of their dipole interaction contribute some electronic character to the nuclear-spin levels and modify its magnetic moment [16]. The Rabi frequency of nuclear spin is about 100 kHz, which reaches 20% of the Zeeman splitting, such fast manipulation also benefits from the electron enhance effect. We discuss the validity of rotating wave approximation (RWA) in Supplementary Note 1.

## Synchronization of pulse generators and phase calibration

Synchronization of the MW and RF generators is one of the main challenges in this experiment. We use the same clock reference for all the generators. For each cable connection and pulse sequence, we measure the phase of prepared state as we scan the phase of input MW pulses. This gives us a phase relation between MW and RF channels, which is used to compensate the difference between the two rotating frames. We check the phase relation before and after each data acquisition. The phase drift of our system is about 2 ◦ in 2 hours measurement.

## Data normalization and state tomography

Since the population information of electron spin is the only directly measurable signal in NV system, we normalize all the data to the fluorescence intensity of electron spin | 0 〉 state. Specifically, we apply two readout pulses (300 ns) at the end of each measurement. See pulse sequences in Fig. 2(b) and

FIG. 6: Coherent manipulation of electron spin and nuclear spin at ESLAC (a) . ODMRspectrum of electron spin at B = 507 Gauss. (b) Pulse sequence to manipulate electron and nuclear spins at ESLAC. (c) ODMRspectrum and (d) Rabi oscillation of the nearby 13 C nuclear spin. At ESLAC, electron spin and nearby nuclear spins (including host 14 N nuclear spin and nearby 13 C nuclear spins) can be polarized by a short laser pulse. Then spin states are manipulated by resonant MW or RF pulses. Electron spin states are readout by counting the fluorescence intensity of NV centre; nuclear spin states are mapped to electron spin and readout in the same way. The resonant frequency of this nuclear spin is 495 kHz, which is slightly modified by centre electron spin.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-8678e6d9-02af-4c6b-ae22-f0481c5aa1a9.jpg)

Fig. 6(b). The first readout pulse gets the instant population information of NV electron spin, and the second readout pulse (1 µ s later) records a reference for the first one, as electron spin is polarized to | 0 〉 again after the 1 µ s laser excitation. The ratio between the first signal and the reference signal is used for further data analysis, such as phase estimation or state tomography.

To carry out state tomography, we adopt the method detailed in Ref. [18, 27]. Total three working transitions, | 0 ↑ 〉 ⇔ | 0 ↑〉 , | 0 ↑〉 ⇔ | 1 ↑〉 and | 0 ↓〉 ⇔ | 1 ↓〉 are selected. The real and imaginary parts of the matrix elements in each working transition are measured by using RF (or MW) pulses of 0 ◦ and 90 ◦ phases, respectively. Other three transitions are measured in the same way, but extra transfer pulses are added before Rabi measurement in the working transitions. The full procedure of state tomography can be found in Supplementary information (Supplementary Fig. 2).

Acknowledgement This work was supported by the National Basic Research Program of China ('973' Program under Grant Nos. 2014CB921402 and 2015CB921103), National Natural Science Foundation of China (under Grant No. 11175248), and the Strategic Priority Research Program of the Chinese Academy of Sciences (under Grant Nos. XDB07010300 and XDB01010000).

Author Contributions X.-Y.P. and H.F. designed the experiment. X.-Y.P. is in charge of experiment, H.F. is in charge of theory. G.-Q.L., Y.-C. C. and X.-Y.P. performed the experiment. Y.-R. Z. and J.-D. Y carried out the theoretical study.

- G.-Q.L. and Y.-R. Z. wrote the paper. All authors analysed the data and commented on the manuscript.
- Competing Interests The authors declare that they have no competing financial interests.

Correspondence Correspondence and requests for materials should be addressed to X.-Y.P or H.F.

- ∗ Electronic address: hfan@iphy.ac.cn
- † Electronic address: xypan@aphy.iphy.ac.cn
- [1] Helstron, C. W. Quantum detection and estimation theory. (Academic Press, New York, 1976).
- [2] Israel, Y ., Rosen, S. &amp; Silberberg, Y . Supersensitive polarization microscopy using NOON states of light. Phys. Rev. Lett. 112 , 103604 (2014).
- [3] Giovannetti, V., Lloyd, S. &amp; Maccone, L. Advances in quanum metrology. Nat. Photon. 5 , 222-229 (2011).
- [4] Giovannetti, V., Lloyd, S. &amp; Maccone, L. Quantum metrology. Phys. Rev. Lett. 96 , 010401 (2006).
- [5] Giovannetti, V. &amp; Maccone, L. Sub-Heisenberg estiamtion strategies are ineffective. Phys. Rev. Lett. 108 , 210404 (2012).
- [6] Nagata, T., Okamoto, R., O'Brien, J. L., Sasaki, K. &amp; Takeuchi, S. Beating the standard quantum limit with four-entangled photons. Science 316 , 726-729 (2007).
- [7] Afek, I., Ambar, O. &amp; Silberberg, Y. High-NOON quantum and classical light states by mixing quantum and classical light. Science 328 , 879-881 (2010).
- [8] Xiang, G. Y., Higgins, B. L., Berry, D. W., Wiseman, H. M. &amp; Pryde G. J. Entanglement-enhanced measurement of a completely unknown optical phase. Nat. Photon. 5 , 43-47 (2011).
- [9] Higgins, B. L., Berry, D. W., Barlett, S. D., Wiseman, H. M. &amp;Pryde, G. J. Entanglement-free Heisenberg-limited phase estimation. Nature 450 , 393-396 (2007).
- [10] Meyer, V., Rowe, M. A., Kielpinski, D., Sackett, C. A., Itano, W. M., Monroe, C., &amp; Wineland, D. J. Experimental demonstration of entanglement-enhanced rotation angle estimation using trapped ions. Phys. Rev. Lett. 86 , 5870 (2001).
- [11] Schleier-Smith, M., Leroux, H. I. D. &amp; Vuleti´ c, V . States of an ensemble of two-level atoms with reduced quantum uncertainty. Phys. Rev. Lett. 104 , 073604 (2010).
- [12] Jo, G. B. et al. Long phase coherence time and number squeezing of two Bose-Einstein condensates on an atom chip. Phys. Rev. Lett. 98 , 030407 (2007).
- [13] Sørensen, A., Duan, L. M., Cirac, J. I. &amp; Zoller, P. Manyparticle entanglement with Bose-Einstein condensates. Nature 409 , 63-66 (2001).
- [14] Gruber, A., Dr¨ abenstedt, A., Tietz, C., Fleury, L., Wr achtrup, J. &amp;von Borczyskowski, C. Scanning confocal optical microscopy and magnetic resonance on single defect centers. Science 276 , 2012-2014 (1997).
- [15] Jelezko, F., Gaebel, T., Popa, I., Domhan, M., Gruber, A. &amp; Wrachtrup, J. Observation of coherent oscillation of a single nuclear spin and realization of a two-qubit conditional quantum gate. Phys. Rev. Lett. 93 , 130501 (2004).
- [16] Childress, L., Dutt, M. V. G., Taylor, J. M., Zibrov, A. S., Jelezko, F., Wrachtrup, J., Hemmer, P. R. &amp; Lukin, M. D. Coherent dynamics of coupled electron and nuclear spin qubits in diamond. Science 314 , 281-285 (2006).
- [17] Gurudev Dutt, M. V., Childress, L., Jiang, L., Togan, E., Maze, J., Jelezko, F., Zibrov, A. S., Hemmer, P. R. &amp; Lukin, M. D. Quantum register based on individual electronic and nuclear spin
- qubits in diamond. Science 316 , 1312-1316 (2007).
- [18] Neumann, P., Mizuochi, N., Rempp, F., Hemmer, P., Watanabe, H., Yamasaki, S., Jacques, V., Gaebel, T., Jelezko, F. &amp; Wrachtrup, J. Multipartite entanglement among single spins in diamond. Science 320 , 1326-1329 (2008).
- [19] Neumann, P., Beck, J., Steiner, M., Rempp, F., Fedder, H., Hemmer, P. R., Wrachtrup, J. &amp; Jelezko, F. Single shot readout of a single nuclear spin. Science 329 , 542-544 (2010).
- [20] Robledo, L., Childress, L., Bernien, H., Hensen, B., Alkemade, P. F. A. &amp; Hanson, R. High-fidelity projective read-out of a solid state spin quantum register. Nature 477 , 574-578 (2011).
- [21] Maurer, P. C. et al. Room-temperature quantum bit memory exceeding one second. Science 336 , 1283-1286 (2012).
- [22] van der Sar, T. et al. Decoherence-protected quantum gates for a hybrid solid-state spin register. Nature 484 , 82-86 (2012).
- [23] Waldherr, G. et al. Quantum error correction in a solid-state hybrid spin register. Nature 506 , 204-207 (2014).
- [24] Taminiau, T. H., Cramer, J., van der Sar, T., Dobrovitski, V. V. &amp; Hanson, R. Universal control and error correction in multi-qubit spin registers in diamond. Nat. Nanotech. 9 , 171-176 (2014).
- [25] Shi, F. Z. et al. Room-temperature implementation of the deutsch-jozsa algorithm with a single electronic spin in diamond. Phys. Rev. Lett. 105 , 040504 (2010).
- [26] Xu, X. K. et al. Coherence-protected quantum gate by continuous dynamical decoupling in diamond. Phys. Rev. Lett. 109 , 070502 (2012).
- [27] Liu, G. Q., Po, H. C., Du, J. F., Liu, R. B. &amp; Pan, X. Y . Noise resilient quantum evolution steered by dynamical decoupling. Nat. Commun. 4 , 2254 (2013).
- [28] Pan, X. Y., Liu, G. Q., Yang, L. L. &amp; Fan H. Solid-state optimal phase-covariant quantum cloning machine. Appl. Phys. Lett. 99 , 051113 (2011).
- [29] Chang, Y. C., Liu, G. Q., Liu, D. Q., Fan H. &amp; Pan, X. Y . Room temperature quantum cloning machine with full coherent phase control in nanodiamond. Sci. Rep. 3 , 1498 (2013).
- [30] Hadden, J. P., Harrison, J. P., Stanley-Clarke, A. C., Marseglia, L., Ho, Y.-L. D., Patton, B. R., O'Brien, J. L. &amp; Rarity, J. G. Strongly enhanced photon collection from diamond defect centers under microfabricated integrated solid immersion lenses. Appl. Phys. Lett. 97 , 241901 (2010).
- [31] Smeltzer, B., Childress, L. &amp; Gali, A. 13 C hyperfine interactions in the nitrogen-vacancy centre in diamond. New J. Phys. 13 , 025021 (2011).
- [32] Jacques, V. et al. Dynamic polarization of single nuclear spins by optical pumping of nitrogen-vacancy color centers in diamond at room temperature. Phys. Rev. Lett. 102 , 057403 (2009).
- [33] Smeltzer, B., McIntyre, J. &amp; Childress, L. Robust control of individual nuclear spins in diamond. Phys. Rev. A 80 , 050302(R) (2009).
- [34] Dr´ eau, A., Spinicelli, P., Maze, J. R., Roch, J.-F. &amp; Jacques V . Single-shot readout of multiple nuclear spin qubits in diamond under ambient conditions. Phys. Rev. Lett. 110 , 060502 (2013).
- [35] Zhao, N. et al. Sensing single remote nuclear spins. Nat. Nanotech. 7 , 657-662 (2012).
- [36] Kolkowitz, S., Unterreithmeier, Q. P., Bennett, S. D. &amp; Lukin, M. D. Sensing distant nuclear spins with a single electron spin. Phys. Rev. Lett. 109 , 137601 (2012).
- [37] Taminiau, T. H., Wagenaar, J. J. T., van der Sar, T., Jelezko, F., Dobrovitski, V. V. &amp; Hanson, R. Detection and control of individual nuclear spins using a weakly coupled electron spin. Phys. Rev. Lett. 109 , 137602 (2012).
- [38] Zhang, Y. L., Zhang, Y. R., Mu, L. Z. &amp; Fan, H. Criterion for remote clock synchronization with Heisenberg-scaling accuracy.

Phys. Rev. A 88 , 052314 (2013).

- [39] Zhang, Y. R., Jin, G. R., Cao, J. P., Liu, W. M. &amp; Fan, H. Unbounded quantum Fisher information in two-path interferometry with finite photon number. J. Phys. A 46 , 035302 (2013).
- [40] Cram´ er, H. Mathematical Methods of Statistics. Ch. 32-34 (Princeton Univ., 1946).
- [41] Braunstein, S. L. &amp; Caves, C. M. Statistical distance and the geometry of quantum states. Phys. Rev. Lett. 72 , 3439-3443 (1994).
- [42] Liu, G. Q., Pan, X. Y., Jiang, Z. F., Zhao, N. &amp; Liu, R. B. Controllable effects of quantum fluctuations on spin free-induction decay at room temperature. Sci. Rep. 2 , 432 (2012).
- [43] Liu, G. Q. et al. Protection of centre spin coherence by dynamic nuclear spin polarization in diamond. Nanoscale 6 , 10134-10139 (2014).

## Supplementary Materials

## Supplementary Note 1: Validity of rotating wave approximation and coherence of nuclear spin

As mentioned in the main text, the achieved Rabi frequency of nuclear spin is not small compared with the energy gap ( | 0 ↑〉 → | 0 ↓〉 ), thus we need to evaluate whether rotating wave approximation (RWA) still works well under this circumstance. We measure nuclear pulse-ODMR spectrum and Rabi oscillation with RF pulses of different driven powers, as shown in Supplementary Fig. 1(a-b). The minimum Rabi fre-

FIG. S1: Coherence of nuclear spin. (a) Pulse-ODMR of 13 C nuclear spin with RF pulses of different driven power. The resonant frequencies are the same for all the measurement (evidence of weak driven). (b) Rabi oscillation of nuclear spin under weak RF driven. The measured Rabi frequency is about 10 kHz, with typical envelop decay time of T 1 ρ = 1 3 . ms. (c) FID of nuclear spin when electron spin is at m S = 0 state, with dephasing time of T ∗ 2 n ( m s = 0) = 270 s. (d) FID of nuclear spin when electron spin is at m S = -1 state, with dephasing time of T ∗ 2 n ( m s = -1) = 212 s. The square, circle and triangle are experiment data with detuning of 0 kHz, 5 kHz and 10 kHz, respectively. Solid lines are fitting to them.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-3190acfe-59da-4267-bef7-ee233f688b24.jpg)

quency is only 10 kHz with a RF power of -24 dBm(at signal generator), which is much less than the energy gap of nuclear spin and RWA works well at this power. We then choose π pulse of this RF power ( 49 µ s) to measure pulse-ODMR spectrum of nuclear spin. The resonant frequency between | 0 ↑〉 and | 0 ↓〉 transition is 495 kHz under this weak driven power. We find that the resonant frequencies are the same for all the measured RF pulses, including the one, which corresponds to a Rabi frequency of about 20% of nuclear spin energy gap. Therefore, we conclude that RWA works well for all the measurements in this experiment. This conclusion conforms with the results in Supplementary Ref. 1 and Supplementary Ref. 2, where RWA works well with a Rabi frequency less than half of the energy gap of two-level system (NV electron spin).

We then consider the coherence of nuclear spin. Supplementary Fig. 1(c) and (d) present the FID signals of nuclear spin under ESLAC, for both m S = 0 and m S = -1 states of electron spin. The dephasing time of nuclear spin ( T ∗ 2 n = 270 µ s for m S = 0 state and 210 µ s for m S = -1 state) is shorter than the result in Supplementary Ref. 3 This may be caused by the complicated spin bath of this NV center. However, similar to the case of electron spin, the dephasing time is not the limitation of nuclear manipulation duration. The Rabi envelope decay time of this nuclear spin ( T 1 ρ ) is more than 1 ms. The half π pulse, which is used to generate the superposition state of nuclear spin, is only 5 µ s and much shorter than the dephasing time. So we ignore the dephasing effect of nuclear spin in the metrology experiment.

FIG. S2: Procedure to carry out state tomography on electronnuclear two-qubit system. (a) The three solid arrows are selected working transitions, which can be driven by MW/RF pulse directly. The real and imaginary parts of the matrix elements in each working transition are measured by RF (or MW) pulses of 0 ◦ and 90 ◦ phases, respectively. (b) For the other three transitions (dash arrow), one or two pulses are applied to transfer the state information to working transition ( | 0 ↑〉 → | 0 ↓〉 ) before Rabi measurement. (c) Pulse duration and phase of state tomography. The diagonal elements of density matrix are measured three times, and the mean of these measurements are used.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-b2cc87ec-4e35-45ad-a53a-030549138b20.jpg)

FIG. S3: Data processing for standard deviation and variance of phase. (a) Linear fitting for standard deviation δϕ 2 against repeat number ν in the log-log scale with adjusted R-squares 0.967 and 0.881 for single state and entangled state. (b) Linear fitting for standard deviation subtracted the system error δϕ se in the log-log scale for single state and entangled state. Adjusted R-squares are 0.962 and 0.904. (c) The variance of phase is fitted by function δϕ 2 = a/ν + c where c represents the squared system error δϕ 2 . (d) Linear fitting for variance subtracted the squared system error in the log-log scale. Adjusted R-squares are 0.916 and 0.893 for single state and entangled state, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0480v2-figure-2c067671-3a34-4a8d-8dc1-85b77b02bee2.jpg)

## Supplementary Note 2: Other data processing methods for standard deviation and variance

Regardless of system error, the variance of phase δϕ 2 with a sufficiently large number of measurements ν will be approximately normally distributed as δϕ 2 ∝ 1 /ν , which is based on the classical central limit theorem. This fact can also be explained by the additive property of Fisher information and is shown in Eq. (1) in the main text. Therefore, in Fig. 4(b) we set the exponent of ν as 0 5 . for standard deviation (SD) and a conclusive result is given. In Supplementary Fig. 3, we try other data processing methods for the presentation of the entanglement-enhanced metrology. In Supplementary Fig. 3(a), we use the linear fitting in the log-log scale to analyze the SD, δϕ 2 , against ν . Ideally the slope should be 0.5 and the intercept gives the value that represents the enhancement. However, there always exists the system error which will insult the linear analysis and give an inconclusive result for large number of measurement ν ∼ 1 M. We thus fix the slopes as 0.5 and give the enhancement by reading the intercepts in Supplementary Fig. 3(a). The adjusted R-square, ¯ R 2 , (ranging from 0 to 1 with larger number indicating better fitting) for single state and entangled state are 0.967 and 0.881, respectively. In Supplementary Fig. 3(b), the system error, δϕ se , read in Supplementary Fig. 3(c) is taken into consideration and subtracted out before analyzing, which leads to better adjusted R-squares: 0.962 and 0.904 for single state and en- tangled state, respectively. Therefore, we conclude that the data processing method we use in the main text for SD provides a better fitting to the experimental data. We also use the function δϕ 2 = a/ν + c to fit experimental data of variance, see results in Supplementary Fig. 3(c). Taking out the effect of squared system error δϕ 2 se , we use the linear fitting in the log-log scale and present the enhancement in Supplementary Fig. 3(d) with ¯ R 2 being 0.916 and 0.893 for single state and entangled state, respectively. It shows that the data processing method we use for SD in the main text is better than that for variance.

- ∗ Electronic address: hfan@iphy.ac.cn
- † Electronic address: xypan@aphy.iphy.ac.cn
- [S1] Fuchs, G. D., Dobrovitski, V. V., Toyli, D. M., Heremans, F. J. &amp; Awschalom, D. D. Gigahertz Dynamics of a Strongly Driven Single Quantum Spin. Science 326 , 1520-1522 (2009).
- [S2] Scheuer, J. et al. Precise qubit control beyond the rotating wave approximation. New J. Phys. 16 , 093022 (2014).
- [S3] Dutt, M. V. G. et al. Quantum register based on individual electronic and nuclear spin qubits in diamond. Science 316 , 13121316 (2007).
Gang-Qin Liu Yu-Ran Zhang Yan-Chun Chang Jie-Dong Yue Heng Fan Xin-Yu Pan
Demonstration of Entanglement-Enhanced Phase Estimation in Solid
Precise parameter estimation plays a central role in science and technology.
The statistical error in estimation can be decreased by repeating measurement,
leading to that the resultant uncertainty of the estimated parameter is
proportional to the square root of the number of repetitions in accordance with
the central limit theorem. Quantum parameter estimation, an emerging field of
quantum technology, aims to use quantum resources to yield higher statistical
precision than classical approaches. Here, we report the first room-temperature
implementation of entanglement-enhanced phase estimation in a solid-state
system: the nitrogen-vacancy centre in pure diamond. We demonstrate a
super-resolving phase measurement with two entangled qubits of different
physical realizations: an nitrogen-vacancy centre electron spin and a proximal
${}^{13}$C nuclear spin. The experimental data shows clearly the uncertainty
reduction when entanglement resource is used, confirming the theoretical
expectation. Our results represent an elemental demonstration of enhancement of
quantum metrology against classical procedure.