Quantum memory and non-demolition measurement of single phonon state with nitrogen-vacancy centers ensemble

TL;DR

This paper proposes a method for quantum non-demolition measurement and memory of single phonons in diamond using NV centers, achieving high fidelity and overlap in phonon state transfer.

Contribution

It introduces a novel scheme utilizing strain-induced coupling in diamond NV centers for efficient phonon quantum memory and QND measurement.

Findings

High fidelity (up to 98.57%) in phonon state transfer.

Effective phonon absorption and emission by NV centers ensemble.

Numerical simulations confirm scheme feasibility under realistic conditions.

Abstract

In diamond, the mechanical vibration induced strain can lead to interaction between the mechanical mode and the nitrogen-vecancy (NV) centers. In this work, we propose to utilize the strain induced coupling for the quantum non-demolition (QND) single phonon measurement and memory in diamond. The single phonon in a diamond mechanical resonator can be perfectly absorbed and emitted by the NV centers ensemble (NVE) with adiabatically tuning the microwave driving. An optical laser drives the NVE to the excited states, which have much larger coupling strength to the mechanical mode. By adiabatically eliminating the excited states under large detuning limit, the effective coupling between the mechanical mode and the NVE can be used for QND measurement of the single phonon state. Under realistic experimental conditions, we numerically simulate the scheme. It is found that the fidelity of the absorbing and emitting process can reach a much high value. The overlap between the input and the output phonon shapes can reach $98.57\%$.

Equations

$$H=$$

$$H^{\prime }=$$

$$H^{\prime }=g_m(a^{†}{a}_m+a{a}_m^{†})+\frac{\Omega (t)}{2}(d^{†}+d).$$

$$\displaystyle H^{\prime}=\left[\begin{array}[]{ccc}0&g_{m}&0\\ g_{m}&0&\frac{\Omega(t)}{2}\\ 0&\frac{\Omega(t)}{2}&0\\ \end{array}\right].$$

$$H$$

$$|\Psi ⟩=c_d|D⟩\otimes |{\varphi }_0⟩+|N{⟩}_{-1}|0{⟩}_{+1}|0{⟩}_0|0{⟩}_m\otimes |{\varphi }_1⟩,$$

$$|{\varphi }_1⟩={\int }_{-\Delta {\omega }_e}^{+\Delta {\omega }_e}d\omega c_{\omega }{e}^{†}\left(\omega \right)|{\varphi }_0⟩,$$

$$\dot{c_d}$$

$$\dot{c_{\omega }}$$

$$c_{\omega }\left(t\right)=\sqrt{\kappa /2\pi }{\int }_0^t{e}^{-i\omega \left(t-\tau \right)}{c}_d\left(\tau \right)\cos \theta \left(\tau \right)d\tau ,$$

$$\dot{c_d}=-\frac{\kappa \cos \theta }{\pi }{\int }_0^t\frac{\sin \left[\delta \omega \left(t-\tau \right)\right]}{\left(t-\tau \right)}{c}_d\left(\tau \right)\cos \theta \left(\tau \right)d\tau .$$

$$\delta \left(x\right)=\underset{k\to \infty }{\lim }\frac{1}{\pi }\frac{\sin kx}{x},$$

$$\dot{c_d}=-\frac{\kappa }{2}{c}_d\left(t\right){\cos }^2\theta ,$$

$$c_d=e^{-\frac{\kappa }{2}{\int }_0^t{\cos }^2\theta \left(\tau \right)d\tau },$$

$$c_{\omega }\left(t\right)$$

$$f\left(t\right)$$

$$=$$

$$f\left(t\right)=\frac{\sqrt{\kappa }\sqrt{1+\exp \left(-\frac{\kappa T}{2}\right)}}{\exp \left[-\kappa \left(t-\frac{T}{2}\right)/2\right]+\exp \left[\kappa \left(t-\frac{T}{2}\right)/2\right]},$$

$$H$$

$$U=\exp [\frac{g}{{\omega }_m}\left(a_m^{†}-a_m\right)|E⟩⟨E|]$$

$$\tilde{H}$$

$${\tilde{H}}_r$$

$${\tilde{H}}_r^{{}^{\prime }}=\frac{g{\Omega }_o}{2{\omega }_m}{e}^{i\delta t}{a}_m^{†}|E⟩⟨0|+\frac{g{\Omega }_o}{2{\omega }_m}{e}^{-i\delta t}{a}_m|0⟩⟨E|$$

$${\tilde{H}}_{eff}^{{}^{\prime }}$$

$$\Delta f_s=\frac{g^2{\Omega }_o^2}{2{\omega }_m^2\delta }.$$

$$H_c$$

$$|\Psi >$$

$${\dot{c}}_1=-\frac{{\gamma }_m}{2}{c}_1-i{g}_m{c}_2+{\kappa }_e\sum _{j=1}^n{\tilde{b}}_j,$$

$${\dot{c}}_2=-\frac{{\gamma }_1}{2}{c}_2-i{g}_m{c}_1-i\frac{\Omega \left(t\right)}{2}{c}_3,$$

$${\dot{c}}_3=-\frac{{\gamma }_0}{2}{c}_3-i\frac{\Omega \left(t\right)}{2}{c}_2,$$

Full text

I Introduction

In quantum information processing, one of the key challenges is to realize the effective coupling, or communication, between distant qubits. Usually, the photon is used as flying qubit for quantum communications DLCZ , or inducing the effective coupling between the distant qubits Yin2007 . For example, in superconducting quantum circuit systems, the microwave photon confined in transmission lines resonator is used as a quantum bus Song2017 . Recently, with the advances in the fabrication and manipulation of the mechanical systems QizheIn1 ; WangIn4 ; VladimirIn6 ; ConnellIn7 , an alternative way of using phonon to couple distant qubits was proposed SchuetzIn2 . The surface acoustic wave (SAW) has been successfully coupled with superconducting qubits Patricio2016 ; Martin2014 . The advantages of the phonon quantum bus compared with the photon quantum bus are the high quality factor and small effective size Patricio2016 ; SchuetzIn2 ; painter2015 . The speed of SAW is $5$ orders of magnitude slower than the speed of light Martin2014 ; ManentiIn3 . Therefore, the wavelength of SAW at GHz is also $5$ orders of magnitude smaller than the microwave light. In this way, the individual superconducting qubit addressing by SAW is possible. There are many schemes concerning single phonons, for example, the phonon states preparation Galland3 ; Matsuda5 ; Matsuda6 , detection Galland3 ; Matsuda5 ; Matsuda6 ; Yanay4 ; Chai7 ; Bai8 ; Woolley11 ; Ohm9 , and the phonon mediated interface Neto10 .

As we know, in quantum optics experiments, the single photon detectors are widely used. For microwave photons, the quantum non-demolish (QND) measurement for photon number states has been realized in both cavity QED and superconducting circuit QED systems Haroche2007 . The QND measurement in circuit QED systems can be used for quantum error corrections. In order to further develop the phonon based quantum information processing, the efficient single phonon detector is needed. The ultimate goal is to realize QND measurement on single phonon state. Inspired by the QND measurement for photon number states, we use the strain induced phonon and the NVE coupling in the diamond to get the strong nonlinearity for the QND measurement of the single phonon state. The geometry of the mechanical crystal structure contains a diamond crystal with rectangular holes arranged periodically which can precisely manipulate mechanical vibrations painter2009 .

The nitrogen-vacancy (NV) center, which consists a substitutional nitrogen atom and adjacent vacancy in diamond, is one of the most promising system for solid-state quantum information processing. NV centers can be controlled by microwave with long coherence time even at room temperature. The diamond lattice vibration, or phonon, can couple with the NV centers ensemble (NVE) Ovartchaiyapong2014 ; Golter2016 and be cooled by the NVE MacQuarrie2017 . There are two different kind of mechanism for coupling the NV center and the phonons. The first on is magnetic field induced coupling Rabl2008 ; Rabl2010 ; du2009 ; zhou2009 ; yin2013 ; ma2016 ; cai2017 , and the second one is the strain induced coupling Bennett2013 ; wang2016x ; wang2016l . Here we focus on the second method. The strain induced coupling between NV centers electron spins in the ground state and the phonons is usually quite small, e.g. less than $1$ kHz Bennett2013 . Recently, the excited-state electron-phonon coupling has been reported in experiments wang2016l . It is about $6$ orders of magnitude stronger than the ground-state electron-phonon coupling maze2011 ; doherty2011 . By taking the advantage of strong coupling, the quantum control of the internal states of a NV center has been realized by using optomechanical sideband transitionwang2016l . We can use this strong coupling mechanism to induce the strong effective nonlinear coupling between the phonon and the NVE, and realized QND phonon number measurement.

In this work, the strain-induced spin-phonon coupling is employed in designing the scheme for the single phonon absorption, emission and QND measurement. We consider the NVE situated a few $\mu m$ below the diamond surface which is shown in Fig. 1(a). The coupling strength between the NVE and the single phonon can be enhanced by a factor $\sqrt{N}$ through the collective excitation, which can reach the strong coupling regime. The resonant frequency of the phonons propagating in the diamond can be controlled by the rectangular holes shown in the upper diagram of Fig. 1(a), which are periodically arranged on the diamond chip. The single phonon with the resonant frequency can be absorbed by the NVE, and the absorption will induce the resonant frequency shift which can be revealed in the phonon absorption spectrum, then we can detect the single phonon state through probing the resonance frequency shift of the phonons in the diamond.

The inverse process of the single phonon detection is a emission process. The emitted phonon shape can be controlled by regulating the driving pulse acting on the NVE. The similarity between the emitted and the absorbed single phonon reaches $98.57\%$ if we inverse the driving pulse when emitting compared with the absorbing process. If the single phonon state is in quantum state before absorbing, the size and shape of the phonon will not be changed when emitting, this process is a quantum non-demolition (QND) measurement process, which realizes an ideal projective measurement that leaves the system in an eigenstate of the phonon number Haroche2007 ; Braginsky14r1 ; Grangier14r2 ; Nogues14r3 .

II Model

As schematically shown in Fig. 1(a), we consider a diamond chip with periodically arranged rectangle holes to adjust the refractive index for the surface acoustic wave (SAW). The NVE are located near the diamond surface coupling to the laser, microwave and SAW.

With a zero magnetic field, the spin-triplet ground state of the NV center splits into two energy levels, $M_{s}=0$ and the nearly degenerate sublevels $M_{s}=\pm 1$ Manson1r49 . We apply an external magnetic field $\overrightarrow{B}_{ext}$ along the crystalline direction $[100]$ of the NV center Saito12 to split the degenerate sublevels $M_{s}=\pm 1$, which results in a three-level system denoted by $|0\rangle=|^{3}A,M_{s}=0\rangle$, $|-1\rangle=|^{3}A,M_{s}=-1\rangle$, $|+1\rangle=|^{3}A,M_{s}=+1\rangle$, respectively. $|E\rangle$ is an excited state with an energy gap of about 1.189eV to the ground state.Rogers2008 The schematic diagram of the NV center energy structure is shown in Fig. 1(b).

The NV centers exist in a diamond crystal and they can sense the lattice vibrations. If there is a phonon produced in the lattice, the NV center can absorb it when the mechanical frequency and the energy splitting of the NV centers are matched Bennett2013 ; Yin2015a . For this system, the phonon mode $a_{m}$ with the frequency $\omega_{m}$ couples to the NV centers with transition of $|-1\rangle$ to $|+1\rangle$. A classical microwave field $\varepsilon\left(t\right)=\widetilde{\varepsilon}\left(t\right)e^{-i\omega_{01}t}$ drives the transition $|+1\rangle$ to $|0\rangle$ with a Rabi oscillation frequency $\Omega(t)$ ($\omega_{01}$ is the frequency splitting between the levels |0\rangle\and$\ |+1\rangle$). The Rabi frequency $\Omega(t)$ can be written as $\Omega(t)=\Omega_{0}\alpha\left(t\right)$. Within a good approximation, we assume that $\widetilde{\varepsilon}\left(t\right)$ increase gradually from zero with $\widetilde{\varepsilon}\left(0\right)\simeq 0$ to a finite strength. Because $\alpha(t)$ is proportional to $\widetilde{\varepsilon}(t)$ ($\alpha(t)\propto\widetilde{\varepsilon}(t)$), then we can get $\Omega\left(t\right)=\Omega_{0}^{\prime}\widetilde{\varepsilon}\left(t\right)$, which discribes that the shape of the amplitude of the classical driving microwave field totally decides the Rabi frequency.

Neglecting the dissipation of the NVE system, the Hamiltonian of the NVE system including the coupling between the mode $a_{m}$ and the output can be written as (setting $\hbar=1$) Duan2003 ; Yin2015

[TABLE]

Where $G_{m}$ is the coupling rate between single NV center and the phonon, which is typically small. $e\left(\omega\right)$ denotes the one-dimentional phonon modes which are free in the diamond and couple to the phonon mode $a_{m}$. We need to consider the free propagating modes within a finite bandwidth $\left[\omega_{e}-\Delta\omega_{e},\omega_{e}+\Delta\omega_{e}\right]$ that couple to the NV mode with the carrier frequency $\omega_{e}$. Within this bandwidth, the coupling between $e(\omega)$ and the NV is a constant approximately which is denoted by $\sqrt{\kappa/2\pi}$ for convenience. $\kappa$ is the effective decay rate.

To obtain the relationship between the driving pulse shape and the output phonon pulse shape, a simple picture by neglecting the dissipation of the NVE system and the coupling of the mode $a_{m}$ to the output is studied. Initially, the NVE are cooled to the ground state $|0\rangle$, and a $\pi$ pulse is driven between the state $|0\rangle$ and $|-1\rangle$, and then the NVE is in $|-1\rangle^{\otimes N}$. We assume that the phonon mode frequency $\omega_{m}$ equals to $\omega_{e}$ ($\omega_{e}$ is the frequency splitting between the levels $|-1\rangle$ and $|+1\rangle$). In the rotating frame, the Hamiltonian of the NVE system is

[TABLE]

We map the operator $|+1\rangle_{i}\langle-1|$ to the bosonic operators. $\sqrt{N}a^{{\dagger}}=\sum_{i=1}^{N}|+1\rangle_{i}\langle-1|$, $a^{{\dagger}}a|n\rangle_{+1}=n|n\rangle_{+1}$ means there are $n$ NVs that are in the state $|+1\rangle$. $d^{{\dagger}}=|0\rangle_{i}\langle+1|$ and $d=|+1\rangle_{i}\langle 0|$ are the creation and annihilation operators for the ith NV center, and then the Hamiltonian can be written as follow, where $g_{m}=\sqrt{N}G_{m}$.

[TABLE]

In the bases of $|N\rangle_{-1}|0\rangle_{+1}|0\rangle_{0}|1\rangle_{m}$, $|N-1\rangle_{-1}|1\rangle_{+1}|0\rangle_{0}|0\rangle_{m}$, $|N-1\rangle_{-1}|0\rangle_{+1}|1\rangle_{0}|0\rangle_{m}$, where the state $|p\rangle_{-1}|q\rangle_{+1}|r\rangle_{0}|s\rangle_{m}$ represents that the numbers of the NVs which are in the states $|-1\rangle$, $|+1\rangle$, and $|0\rangle$ are $p$, $q$ and $r$ respectively, and the number of phonons in the diamond is $s$. The Hamiltonian can be written as a matrix as follow in the three bases,

[TABLE]

This Hamiltonian has the well-known dark state $|D\rangle$ with the form $|D\rangle=\frac{-\Omega(t)/2}{\sqrt{g_{m}^{2}+\left[\Omega(t)/2\right]^{2}}}|N\rangle_{-1}|0\rangle_{+1}|0\rangle_{0}|1\rangle_{m}+\frac{g_{m}}{\sqrt{g_{m}^{2}+\left[\Omega(t)/2\right]^{2}}}|N-1\rangle_{-1}|0\rangle_{+1}|1\rangle_{0}|0\rangle_{m}$. The dark state means, the state will remain on the state $|D\rangle$, if the driving pulse $\Omega(t)$ and the coupling strength satisfy the relationship $\cos\theta=\frac{\Omega(t)/2}{\sqrt{g_{m}^{2}+\left[\Omega(t)/2\right]^{2}}}$ and $\sin\theta=\frac{g_{m}}{\sqrt{g_{m}^{2}+\left[\Omega(t)/2\right]^{2}}}$.

The effects of non-zero dissipation will be analyzed in section IV where the phonon detecting efficiency is calculated as well as the overlap of the absorbing and emitting a single phonon with the real lossy system.

III Pulse shape

To obtain the relationship of the pulse shape between the driving and the phonon, we first consider the emitting process from state $|N-1\rangle_{-1}|0\rangle_{+1}|1\rangle_{0}|0\rangle_{m}$ to $|N\rangle_{-1}|0\rangle_{+1}|0\rangle_{0}|1\rangle_{m}$ under ideal conditions. The Hamiltonian of the system can be written as

[TABLE]

Assume that, at the time $t=0$, the Raby frequency $\Omega\left(t\right)=0$, the NVE system are in the state $|N-1\rangle_{-1}|0\rangle_{+1}|1\rangle_{0}|0\rangle_{m}$, after applying a classical driving pulse $\varepsilon\left(t\right)$, the Rabi frequency, which is proportional to the intensity of the driving pulse, changes slowly and is within the adiabatic approximation. Then we can expend the state $|\Psi\rangle$ of the whole system into the following form Duan2003

[TABLE]

where $|\phi_{0}\rangle$ denotes the vacuum state of mode $e\left(\omega\right)$ in the diamond, and

[TABLE]

represents the single phonon output state. Initially, $c_{d}=1$, $c_{\omega}=0$ and $\Omega\left(0\right)=0$. Within the adiabatic approximation, we would like to calculate the time evolution of the whole NVE system state $|\Psi\rangle$ by substituting it into the schr$\ddot{o}$dinger equation $i\partial_{t}|\Psi\rangle=H|\Psi\rangle$. The coefficients $c_{d}$ and $c_{\omega}$ can be got, which satisfy the following evolution equations:

[TABLE]

We can get the solution of $c_{\omega}$ as follow

[TABLE]

substituting the solution into the equation of $c_{d}$, leads to

[TABLE]

The bandwidth $\delta\omega$ satisfies $\delta\omega T\gg 1$, where the operation time $T$ characterizes the time scale for a significant change of $c_{d}$ and $\sin\theta$, so the above function satisfies a $\delta$ function

[TABLE]

then we can obtain

[TABLE]

The solution of $c_{d}$ and $c_{\omega}$ are

[TABLE]

[TABLE]

The single-phonon pulse shape $f\left(t\right)$ can be obtained by the Fourier transformation

[TABLE]

If the Rabi frequency $\Omega\left(t\right)$, which is proportional to the strength of the driving pulse, satisfies the function as $\Omega\left(t\right)=2g_{m}\exp\left[\kappa\left(t-\frac{T}{2}\right)/2\right]$, the output phonon shape will be symmetric. The output phonon shape function can be figured out as

[TABLE]

the results is shown in Fig. 3 for $\kappa/2\pi=1.6\times 10^{5}$, $T=\frac{20}{\kappa}$. We assume that, in the absorption process, we know the phonon shape in advance, such as the shape in Fig. 2. As the absorption process is the time reversal of the emission process, and the phonon shape is symmetrical in the time domain, if we reverse the temporal driving pulse shape on NV centers in the absorbing process, the phonon can be completely absorbed by the NVE Cirac13 . If the phonon shape is not symmetrical in the time domain, it can also be completely absorbed as long as it is reversed simultaneously with the driving pulse.

IV QND measurement

As schematically shown in Fig. 1, we consider the energy levels $|0\rangle$ and $|E\rangle$. After the absorption process, we stop the driving pulse $\Omega$ between the state $|0\rangle$ and $|+1\rangle$ and start the driving pulse $\Omega_{o}$ between the state $|0\rangle$ and $|E\rangle$ to detect the single phonon state.

The Hamiltonian of the detecting system is

[TABLE]

Applying the Schrieffer-Wolff transformation

[TABLE]

to the Hamiltonian gives

[TABLE]

Apply the RWA, and set $H_{0}=\omega_{m}a_{m}^{{\dagger}}a_{m}+(\omega_{0E}-\frac{g^{2}}{\omega_{m}}-\delta)|E\rangle\langle E|$, the Hamiltonian becomes

[TABLE]

Where $\omega_{0E}=\delta+\omega_{m}+\frac{g^{2}}{\omega_{m}}+\omega_{o}$.

Then we apply another RWA, and set $H_{0}^{\prime}=\delta|E\rangle\langle E|$, the Hamiltonian is calculated as

[TABLE]

When the interaction strength of $\frac{\Omega_{o}}{2}$ is far small compared with the detuning $\delta$, the time-average effective Hamiltonian is obtained asJames2

[TABLE]

From the Hamiltonian we can see that the resonant frequency of the phonon crystal will be affected by the state of NVE, which means the resonant frequency shift is proportional to the excitation number in state $|0\rangle$. When there is one NV center absorbing a single phonon, the frequency shift of the phonon crystal is

[TABLE]

The effective dissipation of the excited state $|E\rangle$ is $\gamma_{e}=(\frac{\overline{g}}{\delta})^{2}\gamma_{E}$, where $\overline{g}=\frac{g\Omega_{o}}{2\omega_{m}}$ is the equivalent coupling strength, $\gamma_{E}$ is the decay rate of the energy level $|E\rangle$ with the population of 1. If we set $g=2\pi\times 5MHz$, $\Omega_{o}=2\pi\times 290MHz$, $\omega_{m}=2\pi\times 900MHz$ wang2016l , $\gamma_{E}=2\pi\times 3MHz$ Robledo , and $\delta=2\pi\times 30MHz$, the effective dissipation and the frequency shift can be calculated as $\gamma_{e}=2\pi\times 2.16kHz$ and $\Delta f_{s}=2\pi\times 43.26kHz$ respectively, which implies that the single phonon absorption induced frequency shift can be detected distinctively. In fact, as long as we have $\delta\gg\gamma_{E}$, the demand for the resolution can be satisfied, and the scheme can be realized under the realistic experimental conditions.

V Dissipation effects

In the above sections, the dynamical evolution has been analysed in ideal condition and the relationship between the driving pulse shape and the emitting phonon shape has been obtained, however, there are various losses in the real systems. If we consider the dynamical evolution of the NVE system with the conditional Hamiltonian which include the possible losses, the whole conditional Hamiltonian has the following form:

[TABLE]

Where $\gamma_{0}$ is the dephasing rate in state $|0\rangle$, $\gamma_{1}$ denotes the spontaneous emission of state $|+1\rangle$, and $\gamma_{m}$ denotes the phonon dissipation in the single NV evolution.

For numerical simulations, we need to discretize the field $e(\omega)$ by introducing a finite but small frequency interval $\delta\omega$ between two adjacent mode frequency. Then the number of the total modes we have is $n=\frac{2\Delta\omega_{e}}{\delta\omega}+1$, the frequency of the $jth$ mode $\omega_{j}$ which is denoted by $e_{j}$ is given by $\omega_{j}=\left(j-\frac{n+1}{2}\right)\delta\omega$. The Hamiltonian can be transformed to the form:

[TABLE]

where $\kappa_{e}=\sqrt{\kappa\delta\omega/2\pi}$. Similarly, the state $|\Psi\rangle$ of the NVE system can be expanded by the discretized phonon pulse state with the form $|\phi_{1}\rangle=\sum_{j=1}^{n}\delta\omega b_{j}e_{j}^{\dagger}|vac\rangle$, which is

[TABLE]

Substituting $|\Psi\rangle$ into the schr$\ddot{o}$dinger equation $i\partial_{t}|\Psi\rangle=H|\Psi\rangle$, we can get

[TABLE]

[TABLE]

[TABLE]

[TABLE]

where $\widetilde{b}_{j}=\sqrt{\delta\omega}b_{j}$. We assume that, initially, there is a single phonon state in the diamond, which shape is shown as Fig. 2, and all of the NV centers are in the state $|-1\rangle$, which means $c_{1}=c_{2}=c_{3}=0$, $\sum\left|b_{j}\right|^{2}\neq 0$. If we would like to absorb the phonon effectively, the driving pulse $\Omega^{{}^{\prime}}\left(t\right)$ should be the time reversal of $\Omega\left(t\right)$ in Fig. 2, that is $\Omega^{{}^{\prime}}\left(t\right)=g_{m}\exp\left[\kappa\left(-t+\frac{T}{2}\right)/2\right]$, where $T=20/\kappa$. The shape control of the driving microwave pulse can be easily achieved by modulating an arbitrary wave generator. We set $\gamma_{0}/2\pi=0.16kHz$, $\gamma_{1}/2\pi=0.16kHz$ Gill2013 , $\gamma_{m}/2\pi=0.16kHz$ Tao2014 ; Ovartchaiyapong2014 , $g_{m}/2\pi=0.96MHz$, and $\kappa/2\pi=0.32MHz$. The numerical simulation results of absorbing process is shown in Fig. 3a with the black lines. Under the condition of strong coupling between the NVE and the phonon, the process is quasi-adiabatic, the value of $c_{2}$ is distinctly smaller than $c_{1}$ and $c_{3}$. The first half is the process that the free single phonon entering the strong coupling area with the NVE plays a leading role, and the later half is predominantly the absorbing process. at time $t=T$, the fidelity between the actual state and the ideal state is $99.38\%$. As the dissipations increase, the fidelity decreases obviously which is also shown in Fig. 3a with the blue lines. The red lines in Fig. 3a implies that, it will also reduce the fidelity of the state if the coupling strength and the effective decay rate $\kappa$ do not satisfy the adiabatic condition.

After absorbing and detecting the single phonon, we can apply a driving pulse $\Omega(t)$ on the NVE, the phonon we have just absorbed will then be emitted. The comparison of the input and output phonon pulse shape is shown in Fig. 3b with the $\kappa=0.32MHz$. The black line is the input phonon shape, the red, blue and green ones are the output phonon shape with $g_{m}/2\pi=0.96MHz$, $0.64MHz$ and $0.32MHz$ respectively. we can see that, the more stronger the coupling is, the larger overlap it has between the input and output single phonon shape. When the coupling strength $g_{m}/2\pi=0.96MHz$, the overlap can reach $98.57$. In the first step, the NVE absorbs one phonon from the diamond or not. And then, we can detecting wether the single phonon state was absorbed through the frequency shift of the phononic crystal. Finally, the NVE can emit the phonon with nearly the same shape compared with the absorbed one. In this process, the single phonon state has been measured without changing the phonon shape, which is a true QND measurement. The process of emitting can act as a single phonon source. The whole courses including absorbing, detecting and emitting a single phonon can also serve as a single phonon memory.

VI Conclusion

We have proposed a scheme to realize the QND single phonon state detecting and emitting based on the strain mediated interaction between the NVE and the single phonon. By analyze the dynamical evolution of the real system, we are able to calculate the fidelity of the absorbing process and the overlap between the input and output phonon number state, both of which can reach a very high value. The emitting process can act as a single phonon source and the whole courses including absorbing, detecting and emitting a single phonon can also serve as a single phonon memory. In future, the similar method may also be used to realize the QND measurement for arbitrary Fock states of phonons.

Acknowledgements.

R.X.W. is supported by China Postdoctoral Science Foundation (2016M600999). Z.Q.Y. is supported by National Natural Science Foundation of China (NSFC) (61435007 and 61771278), and Joint Fund of the Ministry of Education of the People’s Republic of China (MOE) (6141A02011604). G.L.L is supported by National Natural Science Foundation of China (NSFC) (20141300566).