Proposal to test quantum wave-particle superposition on massive mechanical resonators

We present and analyze a proposal for a macroscopic quantum delayed-choice experiment with massive mechanical resonators. In our approach, the electronic spin of a single nitrogen-vacancy impurity is employed to control the coherent coupling between the mechanical modes of two carbon nanotubes. We demonstrate that a mechanical phonon can be in a coherent superposition of wave and particle, thus exhibiting both behaviors at the same time. We also discuss the mechanical noise tolerable in our proposal and predict a critical temperature below which the morphing between wave and particle states can be effectively observed in the presence of environment-induced fluctuations. Furthermore, we describe how to amplify single-phonon excitations of the mechanical-resonator superposition states to a macroscopic level, via squeezing the mechanical modes. This approach corresponds to the phase-covariant cloning. Therefore, our proposal can serve as a test of macroscopic quantum superpositions of massive objects even with large excitations. This work, which describes a fundamental test of the limits of quantum mechanics at the macroscopic scale, would have implications for quantum metrology and quantum information processing.


INTRODUCTION
Wave-particle duality lies at the heart of quantum physics. According to Bohr's complementarity principle, 1 a quantum system may behave either as a wave or as a particle depending on the measurement apparatus, and both behaviors are never observed simultaneously. This can be well demonstrated via a single photon Mach-Zehnder interferometer, as depicted in Fig. 1a. An incident photon is split, at an input beam-splitter BS 1 , into an equal superposition of being in the upper and lower paths. This is followed by a phase shift ϕ in the upper path. At the output beamsplitter BS 2 , the paths are recombined and the detection probability in the detector D 1 or D 2 depends on the phase ϕ, heralding the wave nature of a single photon. If, however, BS 2 is absent, the photon is detected with probability 1/2 in each detector, and thus, shows its particle nature. In Wheeler's delayedchoice experiment, 2,3 the decision of whether or not to insert BS 2 is randomly made after a photon is already inside the interferometer. The arrangement rules out a hidden-variable theory, which suggests that the photon may determine, in advance, which behavior, wave or particle, to exhibit through a hidden variable. [4][5][6][7][8][9][10][11] Recently, a quantum delayed-choice experiment, where BS 2 is engineered to be in a quantum superposition of being present and absent, has been proposed. 12 Such a version allows a single system to be in a quantum superposition of a wave and a particle, so that both behaviors can be observed in a single measurement apparatus at the same time. 13,14 This extends the conventional boundary of Bohr's complementarity principle. The quantum delayed-choice experiment has already been implemented in nuclear magnetic resonance, [15][16][17] optics, [18][19][20][21][22][23] and superconducting circuits. 24,25 However, all these experiments were performed essentially at the microscopic scale.
Here, as a step in the macroscopic test for a coherent waveparticle superposition on massive objects, we propose and analyze an approach for a mechanical quantum delayed-choice experiment. Mechanical systems are not only being explored now for potential quantum technologies, 26,27 but they also have been considered as a promising candidate to test fundamental principles in quantum theory. 28 In this manuscript, we demonstrate that, similar to a single photon, the mechanical phonon can be prepared in a quantum superposition of both a wave and a particle. The basic idea is to use a single nitrogen-vacancy (NV) center in diamond to control the coherent coupling between two separated carbon nanotubes (CNTs). 29, 30 We focus on the electronic ground state of the NV center, which is a spin S = 1 triplet with a zero-field splitting D ≃ 2π × 2.87 GHz between spin states |0〉 and |±1〉 [see Fig. 1b]. If the spin is in |0〉, the mechanical modes are decoupled, and otherwise are coupled. Moreover, the mechanical noise tolerated by our proposal is evaluated and we show a critical temperature, below which the coherent signal is resolved.

Physical model
We consider a hybrid system 31,32 consisting of two (labelled as k = 1, 2) parallel CNTs and an NV electronic spin, as illustrated in Fig. 1c. The CNTs, both suspended along thex-direction, carry dc currents I 1 and I 2 , respectively, while the spin is placed between them, at a distance d 1 from the first CNT and at a distance d 2 from the second CNT. When vibrating along theŷ-direction, the CNTs can parametrically modulate the Zeeman splitting of the intermediate spin through the magnetic field, yielding a magnetic coupling to the spin. [33][34][35][36][37] For simplicity, below we assume that the CNTs are identical such that they have the same vibrational frequency ω m and the same vibrational mass m. The mechanical vibrations are modelled by quantized harmonic oscillators with a Hamiltonian where b k (b y k ) denotes the phonon annihilation (creation) operator. The Hamiltonian characterizing the coupling of the mechanical modes to the spin is where S z = |+1〉〈+1| − |−1〉〈−1| is the z-component of the spin, q k ¼ b k þ b y k represents the canonical phonon position operator, and g k = μ B g s y zp G k /ħ refers to the Zeeman shift corresponding to the zero-point motion y zp = [ħ/(2mω m )] 1/2 . Here, μ B is the Bohr magneton, g s ≃ 2 is the Landé factor, and G k ¼ μ 0 I k = 2πd 2 k À Á is the magnetic-field gradient, where μ 0 is the vacuum permeability. In order to mediate the coherent coupling of the CNT mechanical modes through the spin, we apply a time-dependent magnetic field B x t ð Þ ¼ B 0 cos ω 0 t ð Þ; with amplitude B 0 and frequency ω 0 , along thex-direction, to drive the |0〉 → |±1〉 transitions with Rabi frequency We apply a static magnetic field along theẑ-direction to eliminate the Zeeman splitting between the spin states |±1〉. 36 This causes the same Zeeman shift, where Δ ± = D ± ω 0 , to be imprinted on |±1〉, and a coherent coupling, of strength Ω 2 /Δ + , between them, as shown in Fig. 1b. We can, thus, introduce a dark state and a bright state with an energy splitting ≃2Ω 2 /Δ + . In this case, the spin state |0〉 is decoupled from the dark state, and is dressed by the bright state. Under the assumption of Ω=Δ ( 1, the dressing will only increase the energy splitting between the dark and bright states to This yields a spin qubit with |D〉 as the ground state and |B〉 as the exited state. The spin-CNT coupling Hamiltonian is accordingly transformed to where σ x = σ + + σ − , with σ − = |D〉〈B| and σ þ ¼ σ y À . When we further restrict our discussion to a dispersive regime ω q ± ω m ) jg k j, the spin qubit becomes a quantum data bus, allowing for mechanical excitations to be exchanged between the CNTs. By using a time-averaging treatment, 38,39 the unitary dynamics of the system is then described by an effective Hamiltonian (see Supplementary Section 1 for a detailed derivation), H eff = H cnt ⊗ σ z , where and σ z = |B〉〈B| − |D〉〈D|. The Hamiltonian H cnt includes a coherent spin-mediated CNT-CNT coupling in the beam-splitter form, which is conditioned on the spin state. Here, we neglect the direct CNT-CNT coupling much smaller than the spin-mediated coupling, as is described in Supplementary Section 1. Furthermore, we find that the decoupling of one CNT from the spin gives rise to a spin-induced shift of the vibrational resonance of the other CNT. Hence, the dynamics described by H eff can be used to implement controlled Hadamard and phase gates.
Quantum delayed-choice experiment with mechanical resonators Let us first discuss the Hadamard gate. Having I k = I and d k = d gives a symmetric coupling g k = g, and a mechanical beam-splitter coupling of strength Fig. 1 a Demonstration of the wave-particle duality using a Mach-Zehnder interferometer. A single photon is first split at the input beam-splitter BS 1 , then undergoes a phase shift ϕ and finally is observed at detectors D 1 and D 2 . The photon behaves as a wave if the output beam-splitter BS 2 is inserted, or as a particle if BS 2 is removed. In quantum delayed-choice experiments, BS 2 is set in a quantum superposition of being present and absent, and consequently, the photon can simultaneously exhibit its wave and particle nature. b Level structure of the driven NV spin in the electronic ground state. Here we have assumed that the Zeeman splitting between the spin states |±1〉 is eliminated by applying an external field. c Schematic representation of a mechanical quantum delayedchoice experiment with an NV electronic spin and two CNTs. The mechanical vibrations of the CNTs are completely decoupled or coherently coupled, depending, respectively, on whether or not the intermediate spin is in the spin state |0〉, with the dc current I k through the kth CNT, and the distance d k between the spin and the kth CNT Unitary evolution for a time τ 0 = π/(4J) then leads to For the phase gate, we can turn off the current, for example, of the second CNT, so that g 1 = g and g 2 = 0. In this case, a dispersive shift of ≃J is imprinted into the vibrational resonance of the first CNT, which in turn introduces a relative phase ϕ ≃ Jτ 1 after a time τ 1 under unitary evolution. Note that, here, both Hadamard and phase gates are controlled operations conditional on the spin state, as mentioned before. The two gates and their timing errors are analyzed in detail in the Supplementary Section 2.
We now turn to the quantum delayed-choice experiment with the macroscopic CNTs. We assume that the hybrid system is initially prepared in the state where |vac〉 refers to the phonon vacuum and I k is the identity operator for the kth CNT. After the initialization, the currents are tuned to be I k = I, to drive the system for a time τ 0 , and the resulting Hadamard operation splits the single phonon into an equal superposition across both CNTs. Then, we turn off I 2 for a time τ 1 to accumulate a relative phase between the CNTs. While achieving the desired phase ϕ, we turn on I 2 following a spin single-qubit rotation |D〉 → cos(φ)|0〉 + sin(φ)|D〉 40-42 with φ a rotation angle, and hold for another τ 0 for a Hadamard operation. Therefore, this Hadamard gate is in a quantum superposition of being both present and absent. The three steps correspond, respectively, to the input beam-splitter, the phase shifter and the quantum output beam-splitter acting in sequence on a single photon in the Mach-Zehnder interferometer, as shown in Fig. 1a.
The final state of the system therefore becomes where describe the particle and wave behaviors, respectively. The coherent evolution of the system is given in more detail in Supplementary Section 2. We find from Eq. (16) that the mechanical phonon is in a quantum superposition of both a wave and a particle, and thus can exhibit both characteristics simultaneously. By applying microwave pulse sequences to tune the rotation angle φ, an arbitrary wave-particle superposition state can be prepared on demand. In the case of φ = 0, the single phonon behaves completely as a particle, but as a wave for φ = π/ 2. The morphing between them can also be observed by tuning the rotation angle φ. The probability, P k , of finding a phonon in the kth CNT is given by which includes two physical contributions, one from the particle nature and the other from the wave nature. Note that the spin in a mixed state cos 2 φ ð Þj0ih0j þ sin 2 φ ð ÞjDihDj is capable of reproducing the same measured statistics as in Eq. (19). 11 Thus, in order to exclude the classical interpretation and prove the existence of the coherent wave-particle superposition, the quantum coherence between the states |0〉 and |D〉 should be verified. 19,20,24,25 Experimentally, such a verification can be implemented by performing quantum state tomography to show all elements of the density matrix of the spin. 42 Next, we consider how to initialize and measure the mechanical system. Initially, the NV spin needs to be in the state |D〉 (i.e., the ground state of the spin qubit), one CNT, e.g., the first CNT, needs to be in its single-phonon state, and the other CNT, e.g., the second CNT, needs to be in its vacuum state. To prepare such an initial state, we can begin with an arbitrary state ρ ini = ρ 1 ⊗ ρ 2 ⊗ ρ spin , where ρ k (k = 1, 2) and ρ spin are the density matrices of the kth CNT resonator and the spin, respectively. One can apply a 532 nm laser pulse to initialize the spin qubit in the state |0〉, and then apply a microwave π/2-pulse to it, to obtain the superposition state 1 ffiffi 2 p j0i þ j À 1i ð Þ , which is followed by a microwave π-pulse to obtain the spin-qubit excited state |B〉. By using the sideband-cooling technique, 43-47 the CNT resonators can be cooled down to their quantum ground state, i.e., the acoustic vacuum |vac〉. For example, one can couple an auxiliary qubit with a large spontaneous-emission rate to the CNT resonators. 48 Once the mechanical ground state is achieved, one can tune the spinqubit transition frequency ω q to be close to the CNT resonance frequency ω m , such that the spin-CNT coupling is then approximately given by a Jaynes-Cummings-type Hamiltonian When acting for a time equal to π/(2g), such a Hamiltonian can, with the spin qubit in the excited state |B〉, transfer a mechanical excitation to the left CNT. 49 Meanwhile, the spin qubit goes to its ground state |D〉. The desired initial state jDi is then obtained. For the phonon number measurement, we still need ω q ≃ ω m as in the initialization, but the spin qubit is required to be in the ground state |D〉. In this situation, the Rabi frequency between the spin and the mechanical resonator depends on the number of phonons in the resonator. 49-53 Thus by directly measuring the occupation probability of |B〉, the phonon number in each CNT can be obtained. The measurement of the spin state is enabled by the different fluorescence of the states |0〉 and |±1〉. 54 To measure the state of the spin qubit, one can first apply a microwave π pulse to map jDi ! 1 ffiffi 2 p j0i À j À 1i ð Þ and jBi ! 1 ffiffi 2 p j0i þ j À 1i ð Þ , and then apply a microwave π/2 pulse to map 1 ffiffi 2 p j0i À j À 1i ð Þ!j 0i and 1 ffiffi 2 p j0i þ j À 1i ð Þ!j À 1i. By measuring the Rabi oscillations between the states |0〉 and |−1〉 according to spin-state-dependent fluorescence, 55 one can readout the spin-qubit state. If one employs the repetitive-readout technique with auxiliary nuclear spins, the readout fidelity can be further improved. 56 Mechanical noise Before discussing the mechanical noise, we need to analyze the total operation time, τ T = 2τ 0 + τ 1 , required for our quantum delayed-choice experiment. Note that during τ T , we have neglected the spin single-qubit operation time due to the driving pulse length~ns. 57,58 Since 0 ≤ τ 1 ≤ 2π/J, we focus on the maximum τ T : τ max T ¼ 5π= 2J ð Þ. A modest spin-CNT coupling g/2π = 100 kHz, which can be obtained by tuning the current I and the distance d (see Supplementary Section 1), is able to mediate an effective CNT-CNT coupling J/2π ≃ 12 kHz, thus giving τ max T ' 0:1 ms. The relaxation time T 1 of a single NV spin at low temperatures can reach up to a few minutes. Moreover, with spin echo techniques, a single spin in an ultra-pure diamond example typically has a dephasing time T 2 ≃ 2 ms even at room temperature, 59 corresponding to a dephasing rate γ s /2π ≃ 80 Hz. When dynamical decoupling pulse sequences are employed, the dephasing time can be made even close to one second at low temperatures. 60 These justify neglecting the spin decoherence. In this case, the mechanical noise dominates the dissipative processes. The dynamics of the system is therefore governed by the following master equation, where ρ(t) is the density operator of the system, γ m is the mechanical decay rate, n th = [exp(ħω m /k B T) − 1] −1 is the equilibrium phonon occupation at temperature T, and L o ð Þρ t ð Þ ¼ o y oρ t ð Þ À 2oρ t ð Þo y þ ρ t ð Þo y o is the Lindblad superoperator. Here, H(t) is a binary Hamiltonian of the form, with and H 1 ¼ Jb y 1 b 1 σ z . In Eq. (22), we did not include the spin singlequbit operation before the third time interval because the length of the driving pulse is very short, as mentioned above. The master equation in Eq. (21) drives the phonon occupation of the kth CNT to be at time t = τ T . For a realistic CNT, we can set the mechanical linewidth to be γ m /2π = 0.4 Hz, 61 leading to a single-phonon lifetime of τ m = 1/γ m ≃ 400 ms. In this situation, τ m is much longer than the total operation time τ T , γ m τ T ( 1 and, thus, we obtain This shows that, in addition to the coherent signal P k , the final occupation has a thermal contribution n th γ m τ T . In Fig. 2, we demonstrate the morphing behavior between particle and wave at T ≃ 10 mK, according to Eq. (25). To confirm this, we also plot numerical simulations, which are in exact agreement with our analytical expression. The thermal occupation, n th γ m τ T , increases as the phase ϕ, because such a phase arises from the dynamical accumulation as discussed above. However, an extremely long phonon lifetime causes it to become negligible even at finite temperatures, as shown in Fig. 2. We now consider the fluctuation noise. In the limit γ m τ T ( 1, the fluctuation noise δn noise k in the phonon occupation n k is expressed, according to the analysis in the Supplementary Section 4, as where the first term is the vacuum fluctuation, which can be neglected, and the second term is the thermal fluctuation, which increases with temperature. To quantitatively describe the ability to resolve the coherent signal from the fluctuation noise, we typically employ the signal-to-noise ratio defined as The signal-resolved regime often requires R k >1 for any P k . However, the probability P k in the range zero to unity indicates that there always exist some P k such that R k <1, in particular, at finite temperatures. Nevertheless, we find that the total fluctuation noise is kept below an upper bound and further that assuming B 2 <1=2 can make either or both of R 1 and R 2 greater than 1. In this case, at least one CNT signal is resolved for each measurement. The conservation of the coherent phonon number equal to 1 ensures that the unresolved signal can be inferred from the resolved one, which allows the morphing between wave and particle to be effectively observed from the fluctuation noise. To quantify this, we define a signal visibility as, in analogy to the signal-to-noise ratio R k . The ratio R describes the visibility of the total signal rather than the single CNT signals. At zero temperature (n th = 0), the noise originates only from the vacuum fluctuation, and this yields R ) 1. However, at finite temperatures, n th increases as T, causing a decrease in R, as shown in Fig. 3. Therefore, the requirement of R>1 sets an upper bound on the temperature, and as a result, leads to a critical Fig. 2 Morphing between particle and wave characteristics of a CNT mechanical phonon. Phonon occupation a n 1 and b n 2 as a function of the relative phase ϕ and the rotation angle φ. The analytical results (colored surfaces) are in excellent agreement with the numerical simulations (black symbols). Here, in addition to γ s /2π = 200γ m /2π = 80 Hz, we assume that g/2π = 100 kHz, ω m /2π = 2 MHz, Ω = 10ω m , and Δ − = 142ω m , resulting in ω q ≃ 1.5ω m and then J/2π ≃ × 12 kHz, and that n th = 100, corresponding to an environmental temperature of ≃10 mK temperature, The critical temperature linearly increases with J/γ m , as plotted in the inset of Fig. 3. To increase J, we can increase the current I through the CNTs, decrease the distance d between the CNTs, or decrease the spin-qubit transition frequency ω q . Furthermore, the increase in the CNT resonance frequency ω m or the decrease in the CNT loss rate γ m can also lead to an increase in the critical temperature. For modest parameters of J/2π = 12 kHz and γ m /2π = 0.4 Hz, a critical temperature T c of ≃47 mK, which is routinely accessible in current experiments, can be achieved.
Test of macroscopicity We have described the implementation of a quantum paradox with massive mechanical objects with experimentally distinguishable single-phonon excitations. The question arises whether this proposal can be considered as a test of macroscopicity. 62,63 Typical proposals of such tests (as cited below) have been based on implementing superpositions of macroscopically distinguishable states of classical-like systems, which are often referred to as Schrödinger's cat states (see, e.g., ref. 64 ). Sometimes, the meaning of Schrödinger's cat states is limited to "superposition states of macroscopic systems, where the amplitude of their excitations is large". 65 Note, however, that the term "large amplitude" can be understood in various ways. These include the cases (criteria) when (i) the amplitudes of the constituent states of a given superposition are large as in classical systems, or (ii) when these amplitudes are large enough concerning their experimental distinguishability (i.e., compared to the resolution of detectors). Strictly speaking, a state satisfying one of these conditions, does not necessarily satisfy the other. For example, a superposition of coherent states, jψi ¼ N ðjαi þ jβiÞ with N being a normalization constant, is a cat state according to criterion (i) if jαj; jβj ) 1, but cannot be considered as a cat state according to criterion, (ii) if ϵ jα À βj ( 1 is beyond the resolution of detectors. Conversely, |ψ〉 is a cat state according to criterion (ii) if ϵ can be resolved experimentally even if |α|, |β| ≈ 1, i.e., when criterion (i) is not satisfied. In the latter case, when the amplitude of such excitations is not large in classical terms, but still macroscopically distinguishable, the states are sometimes referred to as Schrödinger's kitten states, as, e.g., those generated and measured in ref. 66 . In this sense, the single-phonon wave-particle superposition, given in Eq. (16), can be referred to as a Schrödinger kitten state, since the excitations of the macroscopic mechanical systems are small, i.e., at the single-phonon level. Indeed, the amplitudes of singlephonon excitations are not large enough to satisfy criterion (i). However, such superpositions of single phonons are large enough that the constituent states of the superposition, given in Eqs. (17) and (18), are experimentally distinguishable, thus satisfying criterion (ii). Therefore, such a test of a quantum principle at the low-excitation level of massive mechanical objects can also be viewed as a test at the macroscopic scale, as claimed, e.g., in refs. 67-69 and references therein. We note that a collective degree of freedom of many atoms does not necessarily imply that the system is in a macroscopic quantum state. However, we showed that the studied system of macroscopic resonators can be in a maximally entangled twomode state. This state is described by a non-positive Glauber-Sudarshan P function. This implies that the system itself is quantum. Below we describe the method to amplify the smallexcitation kitten states, given in Eqs. (17) and (18), to a cat state with large excitation.

Amplification of the Schrödinger kitten states
Here we apply the idea and method of ref. 70 to show how to amplify the phonon numbers of the single-phonon superposition states |particle〉 and |wave〉, given in Eqs. (17) and (18), by squeezing the mechanical modes b 1 and b 2 . Thus, these states can become Schrödinger's cat-like states. For simplicity, but without loss of generality, here we consider a squeezing operator acting on the mode b k (k = 1, 2), with r being a squeezing parameter. This squeezing leads to where we have defined the phonon squeezed Fock states |S 0 〉 k = U k |0〉 k and jS 1 i k ¼ U k b y k j0i k , with |0〉 k being the vacuum state of the mechanical-mode b k . As a result, the states |particle〉 and |wave〉 become respectively. The final state |Ψ〉 f becomes The modes b k for k = 1, 2 are transformed, via squeezing, to the Bogoliubov modes described by By using this unitary transformation, one obtains the average phonon numbers of |S 0 〉 k and |S 1 〉 k equal to We note that by applying this unconditional amplification method, one can exponentially increase the distinguishability of the states |S 10 〉 and |S 01 〉. Although, a single-shot distinguishability of the mechanical-mode states jP r i and jW r i is not increased, a tomographic distinguishability of these states in the phase space is increased with the amplified amplitudes of the mechanicalmode excitations. Indeed, the distinguishability of jPi r and jW r i, as measured by the infidelity, IF ¼ 1 À jhW r jP r ij 2 ¼ 1 À jhwavejparticleij 2 , is independent of the squeezing parameter r for a given ϕ. For any ϕ ≠ ±π/2, the states are distinguishable, and the highest distinguishability is for ϕ = 0,π, for which the infidelity is IF = 1/2. Thus, even for such optimal values of ϕ, it is impossible to deterministically distinguish the states jP r i and jW r i from each other in a single-shot experiment. We refer to this property as a single-shot distinguishability. Anyway, these mechanical states can be macroscopically distinguished by performing, e.g., Wigner-function tomography on a number of their copies. Such tomographic distinguishability in phase space indeed increases with the squeezing parameter r, as shown in Fig. 4. Finally, we note that the famous optical prototypes of the Schrödinger's cat states, which are given by the odd and even coherent states, |ψ ± 〉 = N (|α〉 ± |−α〉), cannot be distinguished deterministically in a single-shot experiment either. This is because the coherent states |α〉 and |−α〉 are not orthogonal for finite values of α. Their overlap decreases exponentially with increasing α, so |α〉 and |−α〉 become orthogonal in the limit of large |α|. However, this amplification of α cannot be done deterministically, because this process is prohibited by the nocloning and no-signalling theorems. Indeed, non-orthogonal states cannot be deterministically transformed to orthogonal (thus, completely distinguishable) states. Note that popular methods of amplifying small-amplitude states are based on either (i) probabilistic but accurate amplification or (ii) deterministic but inaccurate cloning. For example, the method described, e.g., in refs. 66,71 is probabilistic, because it is based on conditional measurements performed on two copies of |ψ ± 〉. In contrast to this, the amplification method in ref. 70 , as applied here, corresponds to approximate quantum cloning, i.e., phasecovariant cloning by stimulated emission.

DISCUSSION
We have presented a proposal for a quantum delayed-choice experiment with nanomechanical resonators, which enables a macroscopic test of an arbitrary quantum wave-particle superposition. The ability to tolerate the mechanical noise has also been given here, demonstrating that our proposal can be implemented with current experimental techniques. While we have chosen to focus on a spin-nanomechanical setup, the present method could be directly extended to other hybrid systems, for example, mechanical devices coupled to a superconducting atom. 32,49,72 Recently, an experimental work reported that photons can be entangled in their wave-particle degree of freedom. 22 This indicates that the wave-particle nature of photons may be used to encode flying qubits for long-distance quantum communication. Photons are ideal quantum information carriers, but they are difficult to store. In contrast to photons, long-lived phonons could be used for optical information storage. 73 Our study shows that phonons can also be prepared in a wave-particle superposition state, and that the wave-particle nature of phonons is not more special than their other degrees of freedom. Thus, the waveparticle degree of freedom of phonons may be exploited for storing quantum information encoded in the wave-particle degree of freedom of photons. In addition, optomechanical interactions can couple a mechanical mode to optical modes at different frequencies. 74 Thus, the mechanical wave-particle degree of freedom may be employed to map quantum information encoded in the wave-particle degree of freedom from photons at a given frequency to photons at any desired frequency. The mechanical wave-particle nature, as a new degree of freedom, may find various applications in quantum information.
We believe that the macroscopicity of our single-phonon waveparticle superposition is highly counter-intuitive, as based on a refined version of the quantum paradox, even if the mechanical resonators are in the single-phonon-excitation regime. Indeed, we analyzed a "nested" kitten state, as given in Eq. (16), where the particle and wave states, given in Eqs. (17) and (18), are purely mechanical kitten states for ϕ ≠ ±π/2. Moreover, we have described a method, based on mechanical-mode squeezing, which enables the amplification of small-excitation Schrödinger kitten states, given in Eqs. (17) and (18), to large-excitation Schrödinger cat states of the massive mechanical resonators. For these reasons, an experimental realization of our proposal can be a fundamental test of a coherent wave-particle superposition of massive objects with phonon excitations, which can be increased exponentially by squeezing. Hence, this proposed quantum delayed-choice experiment of massive mechanical resonators not only leads to a better understanding of quantum theory at the macroscopic scale, but also indicates that, like the vertical and horizontal polarizations of photons, the mechanical wave-particle nature, as an additional degree of freedom of phonons, may be widely exploited for quantum information applications.

DATA AVAILABILITY
The data that supports the findings of this study are available in the Supplementary Information file. Additional data are also available from the corresponding authors upon reasonable request.
For both plots, we have assumed Re(α) = 0. It is seen that the region marked in yellow, corresponding to the negative Wigner function for |S 1 〉 k , increases with increasing squeezing parameter r  Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.
Here, we, first, in Sec. I present more details of how to obtain the spin-controlled coherent coupling between separated mechanical resonators. Second, in Sec. II, we show the detailed implementation of the controlled Hadamard gate, the phase gate, and the mechanical quantum delayed-choice experiment. Next, in Sec. III, we derive in detail the phonon occupation of each CNT at finite temperatures. Then, Sec. IV describes the detailed derivation of the fluctuation noise and the detailed analysis of the requirement of resolving the coherent signal from the environment-induced fluctuation. Finally, in Sec. V we show the method of the numerical simulation used in this work. The effective Hamiltonian H eff in the article describes a spin-mediated CNT-CNT coupling conditioned on the NV spin state. This is the basic element underlying our proposal. To understand more explicitly the spin-controlled coupling between the CNTs, in this section we derive in detail the effective Hamiltonian. We consider a hybrid quantum system consisting of two parallel CNTs and an NV electronic spin (a qutrit), as depicted in Fig. 1(a). Here, for convenience, illustrations in Figs. 1(b) and 1(c) in the article are reproduced in Figs. 1(b) and 1(a), respectively. The CNTs, respectively, carry dc currents I 1 and I 2 , both along the +x-direction. A spin is placed between them, at a distance d 1 (d 2 ) from the first (second) CNT. According to the Biot-Savart law, the CNTs can, at the position of the spin, generate a magnetic field B

I. Spin-controlled coherent coupling between separated mechanical resonators
( = x, y, z) is a unit vector in theˆ -direction, µ 0 is the vacuum permeability, and the subscript "cnt" refers to the CNTs. When the CNTs vibrate along theŷ-direction, the magnetic field is parametrically modulated by their mechanical displacements y 1 and y 2 , and then is reexpressed, up to first order, as B cnt = B cntẑ is a first-order modification, and where B (1) cnt = k=1,2 G k y k , with a magnetic-field gradient, (2) Note that, here, y 1 > 0 (y 2 < 0) indicates a decrease in d 1 (d 2 ). Therefore, the sign, (−1) k−1 , in Eq. (1) does not appear in Eq. (2). Furthermore, an external magnetic field, B ext = B x (t)x + B zẑ , is applied to the NV spin. We have assumed, as required below, that B x (t) is a time-dependent component but B z is a dc component. The Hamiltonian governing the NV spin is therefore given by where g s 2 is the Landé factor, µ B the Bohr magneton, D 2π × 2.87 GHz the zero-field splitting, and S the -component of the spin operator S ( = x, y, z). In terms of the eigenstates, {|m s , m s = 0, ±1}, of S z , the operator S is expanded as and accordingly, the Hamiltonian H NV is transformed to We find that the magnetic field along theẑ-direction causes different Zeeman shifts to be imposed, respectively, on the spin states | ± 1 , and also that the magnetic field along thex-direction drives the transition between the spin states |0 and| ± 1 . The quantum treatment of the mechanical motion demonstrates that the mechanical vibrations of the CNTs can be modelled by two single-mode harmonic oscillators with a Hamiltonian where ω k is the phonon frequency and b k (b † k ) is the phonon annihilation (creation) operator. Here, we have subtracted the constant zero-point energy ω k /2. The mechanical displacement y k is accordingly expressed as where q k is the canonical phonon position operator, and y 1/2 , with m k being the effective mass, describes the zero-point (zp) motion. Combining Eqs. (5), (6), and (7) gives the full Hamiltonian of the hybrid system, The last line in Eq. (8) describes a magnetic coupling between the spin and the mechanical modes. In order to realize a tunable detuning between them, B x (t) is chosen to be B x (t) = B 0 cos (ω 0 t) with amplitude B 0 and frequency ω 0 . In a frame rotating at H rot = ω 0 (| − 1 −1| + | + 1 +1|), the full Hamiltonian can be divided into two parts, account for the low-and high-frequency components, respectively. Here, we have defined Roughly, having δ ± = δ ± + 2ω 0 Ω allows one to make the rotating-wave approximation (RWA), and to straightforwardly remove H high . However, as demonstrated in Sec. V, the accumulated error increases during the evolution, causing the dynamics driven by H low to deviate largely from that driven by H F . Thus, we are not using the RWA here. In order to suppress the error accumulation, we need to analyze the effects of H high in the limit δ ± Ω. In such a limit, we can employ a time-averaging treatment for the high-frequency component H high [1,2], and as a result, its effective behavior is described by the following time-averaged Hamiltonian, where the first line corresponds to the energy shifts of the spin states | ± 1 , and the second line describes a coherent coupling between these. Accordingly, the full Hamiltonian H F is approximated to be a time-independent form, As seen in Sec. V, the error accumulation is strongly suppressed when H high is included.
Furthermore, from Eq. (22) it is found that the sequential actions of the terms σ + b 1 and σ − b † 2 , as well as of the counter-rotating terms σ − b 1 and σ + b † 2 , can transfer a mechanical phonon from the left to the right CNT, and the reverse process is caused by their Hermitian conjugates. When restricting our discussion to a dispersive regime, this phonon transfer becomes dominant. Hence, in the dispersive regime the dynamics described by H F in Eq. (22) enables a spin quantum bus for the mechanical phonons and can be used to realize a coherent CNT-CNT coupling.
In order to show more explicitly, we rewrite H F in the interaction picture as The condition in Eq. (23) justifies to use a time-averaging treatment of the Hamiltonian H F [1,2]. In the timeaveraging treatment, all terms in Eq. (24) are considered as high-frequency components and exhibit time-averaged behaviors. Based on this, the dynamics of the system can be determined by an effective Hamiltonian Here, we have assumed that ω k = ω m . As expected, Eq. (25) shows a coherent spin-mediated CNT-CNT coupling, corresponding to the standard linear coupler transformation, which can give rise to a direct phonon exchange. Thus in this case, the spin qubit works as a quantum bus. At the same time, it also shows that the CNT-CNT coupling can be turned off if the intermediate spin is in the state |0 . This is because the NV spin in the state |0 is decoupled from the CNTs, and the mechanical phonons can no longer be transferred from one CNT to another. Specifically, if the spin is in the state |D or |B , the CNTs are coupled; however, if the spin is instead in the state |0 , they are decoupled. Note that in Eq. (25) ac Stark shifts caused to be imposed on the qubit have been excluded because we focus only on the quantum states of the CNTs.
In the last part of this section, we evaluate the direct coupling between the CNTs. For simplicity, we assume that I k = I, d k = d, and that the CNTs have the same length L. The attractive force acting on the kth CNT is where is the force size. The work done by the force is given straightforwardly by After applying a perturbation expansion and then a quantization, this direct CNT-CNT coupling is found to be where For a modest setup [3][4][5][6][7][8][9], m = 1.0 × 10 −22 kg, ω m = 2π × 2 MHz, L = 10 nm, d = 2 nm, and I = 380 nA, we have which is much smaller than the mechanical resonance frequency ω m , and also have which is much smaller than the spin-mediated CNT-CNT coupling, for example, 2π × 12 kHz, as shown in the section below. Therefore, the direct CNT-CNT coupling can be neglected in our setup.

II. Controlled Hadamard gate, phase gate, and mechanical quantum delayed-choice experiment
In order to implement a quantum delayed-choice experiment with macroscopic CNT mechanical resonators, we need a controlled Hadamard gate and a phase gate to act on the CNT mechanical modes. Below, we demonstrate how the effective Hamiltonian in Eq. (25) can be used to make all required gates. Let us first consider the controlled Hadamard gate. Tuning the currents to be I k = I and, at the same time, the distances to be d k = d results in a symmetric coupling g k = g. The effective Hamiltonian H eff is accordingly reduced to H eff = H cnt ⊗ σ z , where is a beam-splitter-type interaction, and where is an effective CNT-CNT coupling strength. In our discussion, the NV spin is restricted to a subspace spanned by {|0 , |D }, where the spin is a control qubit of a Hadamard gate. The spin in the state |D mediates the coherent coupling between the separated CNTs, and causes them to evolve under the Hamiltonian H cnt in Eq. (34). According to the Heisenberg equation of motion, b k (t) = exp (iH cnt t/ ) b k exp (−iH cnt t/ ), the unitary evolution for a time t = τ 0 ≡ π/ (4J) corresponds to a Hadamard-like gate, However, when the spin state is |0 , the two CNTs decouple from each other. In this case, their quantum states remain unchanged under the unitary evolution, yielding We have therefore achieved a spin-controlled Hadamard gate between the CNTs. That is, if the NV spin is in the state |D , then the Hadamard operation is applied to the CNTs, and if the NV spin is in the state |0 , then the states of the CNTs are unchanged. We next consider the phase gate. For the phase gate, we tune the currents to be I 1 = 0 and I 2 = 0, such that g 1 = g and g 2 = 0, causing the effective Hamiltonian in Eq. (25) to become We find from Eq. (40) that there exists a spin-induced shift, J, of the mechanical resonance. This dispersive shift can, in turn, introduce a dynamical phase, φ (t) = Jt, onto the first CNT. With the spin being in the state |D , we solve the Heisenberg equations of motion for the CNTs, and then obtain a phase gate, In fact, similar to the controlled Hadamard gate discussed above, the phase gate can also be controlled by the spin according to Eq. (40).
Having achieved all required gates, we now turn to the detailed description of the macroscopic quantum delayedchoice experiment with CNT resonators. The hybrid system is initially prepared in the state |Ψ i ≡ |Ψ (0) = b † 1 ⊗ I 2 |vac ⊗ |D , where |vac refers to the phonon vacuum of the CNTs and I 2 is the identity operator on the second CNT. First, we turn on the currents of the CNTs and ensure I k = I. After a time τ 0 , a Hadamard operation is applied to the CNTs and accordingly, |Ψ i becomes Then, we turn off the current of the second CNT for a phase accumulation for a time τ 1 . As a consequence, the system further evolves to While achieving the desired phase φ, we make a spin single-qubit rotation |D → cos (ϕ) |0 + sin (ϕ) |D , and have Here, note that, we have ignored the length of the driving pulse of the spin rotation as being of the order of ns, and thus assumed that the state of the CNTs remains unchanged. At the end of the driving pulse, we turn on the current of the second CNT again and hold for another τ 0 to perform a Hadamard gate. This gate is in a quantum superposition of being present and absent. The three operations on the mechanical phonon correspond to the actions, on a single photon, of the input beam splitter, the phase shifter, and the output beam splitter, respectively, in quantum delayed-choice experiments with a Mach-Zehnder interferometer. The final state is therefore given by P2 as a function of the rotation angle ϕ and the relative phase φ. This represents a continuous transition between a particle-type behavior (ϕ = 0) and a wave-type behavior (ϕ = π/2). where describe particle and wave behaviors, respectively. This reveals that the CNT mechanical phonon is in a quantum superposition of both a particle and a wave. The probability of finding a single phonon in the kth CNT is expressed as according to Eq. (46). In Fig. 2, we have plotted the probabilities P 1 and P 2 versus the rotation angle ϕ and the relative phase φ. In this figure we find that the mechanical phonon shows a morphing behavior between particle (ϕ = 0) and wave (ϕ = π/2).
We now consider the timing errors of the Hadamard and phase gates. We first consider the Hadamard gate. We assume that the error of the time required for performing the Hadamard gate is δ 0 , such that the actual evolution time for the gate becomes τ 0 = τ 0 + δ 0 . In order to estimate the effect of this timing error on the gate performance, we introduce a gate fidelity, defined as where |Ψ target,0 is the target state given by the ideal Hadamard gate, and ρ actual,0 (τ 0 ) is the actual state obtained by integrating the exact master equation, given by Eq. (122). In Fig. 3, we plot the gate fidelity, F 0 , versus the timing error δ 0 (red curve). In this figure, the initial state for the Hadamard gate is assumed to be b † 1 |vac |D , so that the target state is |Ψ target,0 = 1 Here, |vac represents the acoustic vacuum state of the CNT resonators. From this figure, we find that for −0.32τ 0 δ 0 0.34τ 0 , the gate fidelity F 0 can be kept above 0.9.
For the phase gate, we assume, as above, that the timing error is δ 1 . Thus, the actual evolution time for the phase gate is τ 1 = τ 1 + δ 1 . We also introduce a gate fidelity, defined as where |Ψ target,1 is the target state given by the ideal phase gate, and ρ actual,1 (τ 1 ) is the actual state obtained from the exact master equation, given in Eq. (122). The gate fidelity F 1 is plotted as a function of the timing error δ 1 in Fig. (3) (blue curve). There, we assumed that the initial state for the phase gate is 1 and that the phase accumulated is equal to π. The target state |Ψ target,1 is, therefore, given by 1 It is seen from this figure that, as long as −0.15τ 1 δ 1 0.12τ 1 , we can obtain the gate fidelity of F 1 > 0.9. We have assumed that δ0/τ0 = δ1/τ1 = δ, and that the initial state is: (i) b † 1 |vac |D for the Hadamard gate, and (ii) 1 |vac |D for the phase gate that accumulates a relative phase π. Here, in addition to γs/2π = 200γm/2π = 80 Hz, we have assumed that g/2π = 100 kHz, ωm/2π = 2 MHz, Ω = 10ωm, and ∆− = 142ωm, resulting in ωq 1.5ωm and then J/2π 12 kHz. We have also assumed that n th = 100, which corresponds to the environment temperature of 10 mK. Here, in addition to γs/2π = 200γm/2π = 80 Hz, we have assumed that g/2π = 100 kHz, ωm/2π = 2 MHz, Ω = 10ωm, and ∆− = 142ωm, resulting in ωq 1.5ωm and then J/2π 12 kHz. We have also assumed that n th = 100, which corresponds to an environmental temperature of 10 mK.
Note that the spin, in a classical mixed state of the form cos 2 (ϕ) |0 0| + sin 2 (ϕ) |D D|, would lead to the same measured statistics in Eq. (49), that is, a local hidden variable model is capable of reproducing the quantum predictions. This is a loophole [10][11][12][13]. However, as discussed in Refs. [14][15][16][17], this loophole can be avoided as long as the second Hadamard operation is ensured to be in a truly quantum superposition of being present and absent. In our proposal, the second Hadamard operation is conditioned on the spin state. If the spin is in the |0 state, then the Hadamard operation is absent; if the spin is in the |D state, then the Hadamard operation is present; if the spin is in a quantum superposition of the |0 and |D states, then the Hadamard operation is in a quantum superposition of being present and absent. To confirm such a quantum superposition, in Fig. (4) we numerically calculate the fidelity, F = f Ψ|ρ actual (τ T ) |Ψ f , between the desired state |Ψ f in Eq. (46) and the actual state ρ actual (τ T ) obtained from the exact master equation in Eq. (122). From this figure, we find that the fidelity is very close to unity even for the finite temperature of T 10 mK. Furthermore, in experiments, in order to exclude the classical interpretation and prove the existence of the coherent wave-particle superposition, the quantum coherence between the states |0 and |D should be verified. Experimentally, this coherence can be prepared by a spin single-qubit operation [18][19][20], and can be verified by performing quantum state tomography to show all the elements of the density matrix of the spin [20].

III. Phonon occupation at finite temperatures
We begin by considering the total operation time, which is given by τ T = 2τ 0 + τ 1 , as discussed in Sec. II. Here, τ 0 = π/ (4J) is the time for the Hadamard gate and τ 1 ∈ [0, 2π/J] is the time for the phase gate. In a realistic setup, we can assume ω m /2π 2 MHz, ω q /2π 3 MHz, and g/2π = 100 kHz, such that J/2π 12 kHz, yielding a maximum total time τ max T = 2τ 0 + τ max 1 0.1 ms, where τ max 1 = 2π/J is the maximum phase gate time. Note that, the operation time τ T depends inversely on the CNT-CNT coupling strength J, but the enhancement in J is limited by the validity of the effective Hamiltonian H eff .
The total decoherence in our setup can be divided into two parts, one from the spin and the other from the CNTs. The spin decoherence in general includes the relaxation and the dephasing. For an NV electronic spin, the relaxation time T 1 can reach up to several minutes at low temperatures and the dephasing time can be T 2 2 ms even at room temperature [21,22]. These justify neglecting the spin decoherence. For the mechanical decoherence, despite a long phonon life, the low mechanical frequency makes the CNT mechanical modes very sensitive to the environmental temperature. In this section and in Sec. IV, we discuss the effects of the mechanical noise on our quantum delayedchoice experiment, and demonstrate that the morphing between wave and particle can still be effectively observed even at finite temperatures.
As a result, the dissipative processes, in the hybrid system considered here, are induced only by the mechanical decoherence, which arises from the vacuum fluctuation and thermal noise. The full dynamics of the system can then be governed by the following master equatioṅ where ρ is the density operator of the system, γ m is the mechanical decay rate, n th = [exp ( ω m /k B T ) − 1] −1 is the equilibrium phonon occupation at temperature T , and Here, H (t) is a binary Hamiltonian of the form, with The three time intervals in Eq. (53) correspond to the first Hadamard gate, the phase gate and the second Hadamard gate, respectively. Note that in Eq. (53), we did not include the spin single-qubit rotation before the third interval because the length of the driving pulse is of the order of ns. We can derive the system evolution step by step.
Let us now consider the first evolution interval 0 < t ≤ τ 0 . During this interval, the coupling of the CNT mechanical modes introduces two delocalized phononic modes, such that H 0 is diagonalized to be R k = P k δn noise k . (113) Note that, here, we did not use δn k to define R k because δn signal k in δn k results from quantum fluctuations of the desired signal, as mentioned previously; and therefore this is not the environmental noise. In order to resolve a signal from the fluctuation noise, the ratio R k is required to be R k > 1. However, Eq. (113) demonstrates that this criterion is not always met for all values of P k , in particular, at finite temperatures. For example, P k = 0 leads directly to R k = 0. To address this problem, we now consider the total fluctuation noise, We further assume that S 2 < P 2 1 + P 2 2 .
Under this assumption, if R k < 1, then R 3−k > 1 for k = 1, 2; otherwise R 1 > 1, R 2 > 1. This means that at least one of the signals, P 1 or P 2 , is resolved for each measurement. Because the coherent phonon number equal to 1 is conserved, and therefore the signals in the two CNTs are complementary, the unresolved signal can be completely deduced from the resolved one. Thus, the criterion in Eq. (115) ensures that the morphing behavior between wave and particle can be observed from the environment-induced fluctuation noise. In fact, for any value of P k , the total noise S is limited by an upper bound, which is independent of P k . Meanwhile, P 2 1 + P 2 2 is also limited by a lower bound √ 2/2. Thereby, in order to meet the criterion given in Eq. (115), it is required that Based on this condition, we can define a signal visibility in analogy to R k . When R > 1, the morphing between wave and particle can be observed, and cannot otherwise. This, in turn, leads to an upper bound on the equilibrium phonon occupation, and therefore an upper bound on the temperature, .
Because τ max T 5π/2J, the critical temperature is In Fig. 6 we plot the signal-to-noise ratios R 1 and R 2 at the temperature T 10 mK. We find that almost all signals can be resolved, and also, as expected, find that when the signal in one CNT is unresolved, the signal in the other CNT is resolved. In fact, the upper bound B is the fluctuation noise in the total phonon occupation b † 1 b 1 + b † 2 b 2 at time τ T . The criterion R > 1 heralds that to resolve the morphing behavior, the fluctuation noise in b † 1 b 1 + b † 2 b 2 (τ T ) is required to be smaller than √ 2/2.

V. Numerical simulations
In order to confirm our analytical results, we need to numerically simulate the dynamics with the full master equation given byρ where σ z = |D D| − |0 0|, and H F is the full Hamiltonian of Eq. (8). Here, we use the Python framework QuTiP [24,25] to set up this problem. However, the full Hamiltonian is time-dependent, and it takes a long time to integrate the corresponding Schrödinger equation or the master equation, in particular, for our case, where all quantum gates result from the deterministic time evolution of the system. Thus, in our numerical simulations, we replace H F with H low + H high , as in Eq. (13). This is a reasonable replacement because in our proposal Ω (tens of MHz) is required to be much smaller than ∆ (up to ∼ GHz). In Fig. 7, we plot the unitary evolution of the phonon occupations, b † 1 b 1 and b † 2 b 2 , of the CNTs. Symbols are the exact results from the full Hamiltonian H F and solid curves are given by the approximate Hamiltonian H low + H high . We find an excellent agreement for a very long evolution time, and thus H F can be very well approximated by H low + H high . For additional comparison, we also plot the phonon occupation evolution driven only by the low-frequency component H low , corresponding to dotted curves. As seen in Fig. 7, owing to the error accumulation, the dynamics of H low deviates largely from the full dynamics of H F , even within one oscillation cycle. With the above replacement, we obtain the numerical simulations plotted in