Photon-phonon-photon transfer in optomechanics

We consider transfer of a highly nonclassical quantum state through an optomechanical system. That is we investigate a protocol consisting of sequential upload, storage and reading out of the quantum state from a mechanical mode of an optomechanical system. We show that provided the input state is in a test-bed single-photon Fock state, the Wigner function of the recovered state can have negative values at the origin, which is a manifest of nonclassicality of the quantum state of the macroscopic mechanical mode and the overall transfer protocol itself. Moreover, we prove that the recovered state is quantum non-Gaussian for wide range of setup parameters. We verify that current electromechanical and optomechanical experiments can test this complete transfer of single photon.

Quantum optomechanics 1,2 can reversibly interconnect radiation to mechanical oscillators at the level of individual quanta [3][4][5][6] . This connection has both fundamental and also practical consequences. At fundamental side, it controllably converts a single photon to single phonon and then, back to individual photons. It proves fundamental aspect that the information carried by the physical state of single quantum can be transferred between totally different physical systems. Also, it says message to quantum thermodynamics 7 , that minimal portion of energy can substantially change its form. At practical side, such the high-fidelity photon-phonon-photon transducers can be used to sense mechanical motion at levels of single phonon level. It can stimulate quantum metrology of mechanical motion 8 beyond the current state of the art. Gaussian quantum entanglement 9 and squeezed states 10 between microwave radiation and mechanical oscillator have been experimentally generated and read out using resolved sideband regime. However, Gaussian states are not sufficient resources for advanced quantum technology. It requires non-Gaussian quantum states 11 to open full space of the applications. Single photon state is basic non-Gaussian resource for optical technology, therefore the single phonon states are as well basic for non-Gaussian quantum optomechanics. A mutual interconversion between a photon and a phonon is thus an important proof of this fundamental correspondence. A single photon is simultaneously a very good test-bed to verify quantum nature of the physical interface between light and mechanics. Recently, nonclassical correlations of phonon with photon at optical frequencies have been unconditionally generated 12 . This opens future possibilities to exploit optomechanics operated at single phonon level. The transducers between optical, mechanical and microwave modes have been considered [13][14][15][16][17] and realized in continuous-wave regime [18][19][20][21] . However, the continuous wave regime with steady states is not suitable when quantum processing and metrology runs in time slots. In this case, all operations are pulsed with optimized temporal shapes. Theory proposals of photon-to-phonon transfer in the pulsed [22][23][24][25] and stroboscopic regime 26 have been recently analyzed. Also squeezed light improved read-out of mechanical state has been separately investigated 27 .
In this paper, we analyze complete photon-phonon-photon coherent transfer for microwave and optical experiments in the pulsed and red side-band resolved regime. We investigate in detail possibility to transfer the negativity of Wigner function and more general quantum non-Gaussianity through an optomechanical system in contrast to studies of engineering of a nonclassical state of a mechanical mode 28,29 . Transfer of both negativity and non-Gaussianity is much more demanding than the transfer of Gaussian entanglement. We verify that the electromechanical platforms 9,30 with microwave radiation are already feasible for such transduction preserving negativity of Wigner function. The optomechanical experiments require coupling higher than critical to reach this possibility. To preserve quantum non-Gaussianity of single photon, optical experiment 12 is however already sufficient. Our analysis can be extended to other systems 31 . We numerically verify all the predictions using analysis beyond adiabatic elimination of intracavity field. We use single photon as a basic test bed for photon-phonon-photon channel, however, it verifies basic quality of coherent transfer for other non-Gaussian states. Our results therefore open path to use photon-phonon-photon transducer for many future experiments in quantum metrology, quantum nonlinear optomechanics and quantum thermodynamics.

Summary of the main results
It has been shown previously 32 that in presence of a strong red-detuned coherent pulse, the quantum states of the mechanical mode of an optomechanical cavity and the resonant with the cavity pulse can in principle be swapped. That is, the quantum state of an incident pulse is transferred to the mechanical mode and the state of the latter is mapped to the reflected light. The scheme analyzed here uses this state swap twice: first time to upload a highly nonclassical (single-photon) quantum state to the mechanical mode and second time to read the evolved state out (see Fig. 1). We were inspired by the quality of state read out in experiment 9 .
First, the pulse in single-photon state in the optical mode described by annihilation operator A in is shone upon the cavity at its resonance frequency ω cav . A strong classical coherent pulse with amplitude α 1 at frequency ω ω − m cav enables the state swap, so that the mechanical state after the pulse reads 27,32 is the transmittivity associated with the partial state swap. Here κ is the cavity decay rate and g 0 is the single-photon optomechanical coupling strength enhanced by the mean intracavity photon number α = N 1 1 2 . If the interaction time τ 1 or optomechanical coupling g N 0 1 is sufficient, the optical state is perfectly mapped to the mechanical mode. In practice, high optomechanical coupling may be inaccessible, or cause unwanted heating due to absorption in mirrors. The temporal duration τ 1 of the pulse is as well limited by mechanical decoherence.
After the upload, the mechanical state is left to evolve for time τ s . During this time the mechanical environment admixes thermal noise to the quantum state of the mechanical mode. This effect can be described by an effective transmittance δ = γτ − e s , where γ is the mechanical damping rate. This admixture happens as well during the upload stage (and consequent readout), but we omit it for the state swap stages as this effect can be reduced by shortening the pulses durations.
Finally, another state swap is performed to read the mechanical state out. In a full analogy with Eq. (1) the state of the output optical mode reads where A vac describes a vacuum mode incident to the cavity, 2 and τ 2 is the temporal duration of the readout pulse.
The complete channel from the input optical state with annihilation operator A in to the output one with A out will be degraded by optical losses caused by e.g. imperfect coupling to the optomechanical cavity. The optical loss performs admixture of vacuum to the signal mode and is characterized by transmittivity η κ κ = / e , where κ and κ e are respectively total cavity decay rate and the cavity decay rate due to the coupling to the waveguide.
All the steps constituting the channel thus can be described by known transmittivities and therefore the quadratures of the output optical state can be expressed in terms of quadratures Q in of input state and quadratures of added noise Q N Here the transmittivity equals  η δ = T T 1 2 2 . The noisy mode Q N in the equation above comprises optical noises coming from imperfect coupling, thermal mechanical noise and the initial mechanical state. We assume all these modes to be zero-mean Gaussian noises with known variances. From the latter we construct the covariance matrix of Q N which is a diagonal 2 × 2 matrix with the diagonal elements equal Here we assumed that the mechanical mode and its environment are initially in thermal states with occupation n 0 and n th respectively and all the auxiliary optical modes are in vacuum.
In order to transfer negativity of the Wigner function of a single-photon state, the variance Eq.
Our goal is thus to increase the transmittivity  preserving moderate added noise variance V N . From the expression for  it is evident that it can be increased by increasing either of T 1 and T 2 . However, increasing T 1 not only increases  but decreases V N by diminishing the ratio of the initial thermal mechanical state in the output. Increase of T 2 is not so crucially important and it causes increase of both  and V N . These simple considerations can be proven by calculating partial derivatives  The principal dependence of V N on  is shown for different transmittances T 1,2 in Fig. 2(a). The mechanical mode is considered to be at equilibrium with the environment at = = n n 1 0 t h . Even this low occupation creates a demand for a very efficient state upload, so that increasing T 1 with constant T 2 helps to suppress the added noise dramatically (dashed lines). It is worth noting that for certain values of T 1 the successful transfer of the negativity is impossible regardless of T 2 . Note, it is evident that is much easier to reach entanglement preserving channel than to transfer negative Wigner function. From the Fig. 2(b) that shows the analysis following the full solution (see Methods) we conclude that the current electromechanical setups 9,30 are capable of the transfer of the negativity of the Wigner function provided there is a single-photon Fock state at input. The correspondence between Fig. 2(a) and (b) allows us to conclude that the simple model based on the effective transmittances provides qualitatively correct understanding of the physics of the transfer.
A major impediment to the successful photon-phonon-photon transfer preserving the negativity of Wigner function is the initial occupation of the mechanical mode. The influence of the mechanical bath during the state swap stages can in principle be suppressed by taking shorter pulses, however, the nonzero initial occupation has to be compensated for either by precooling of the mechanical mode or by a very efficient upload. In spite of our considered experimental scheme the precooling happens to be an additional resource that can relax the requirements to the upload of the quantum state stage. The considerations without the additional precooling therefore happen to be more conservative as such regime is more challenging experimentally and it requires the upload efficiencies close to unity. These efficiencies, however, are realistic at least in the domain of state-of-the-art electromechanics 9 , where the transmittances can reach level as high as > − − T 1 10 1 5 . There is therefore an open way to experimentally test photon-phonon-photon channel preserving negativity of Wigner function.
We should note that in order to transmit negativity of the Wigner function, one should have  > 1/2 which places restrictions on each of individual transmittances, particularly on the efficiency of the coupling, η. That is, with r ∈ (0, + ∞ ). To evaluate the non-Gaussian properties of the restored state of the field we therefore need to evaluate the probabilities p 0 and p 1 corresponding to this state and certify if they belong to the region accessible to Gaussian states. The latter is denoted as  and shaded at Fig. 3. Experimentally both probabilities can be determined from homodyne detection scheme used to verify negativity of Wigner function 33 . Alternatively, single photon counters can be used 35 .
In order to evaluate p 0 and p 1 , we use our estimates for the transmittance  and added noise V N to calculate the Wigner function of the recovered state, from which we find the needed probabilities. Results of this calculation are visualized for recent electro-and optomechanical experiments at Fig. 3. Example Wigner functions are demonstrated in Fig. 4. The demonstration of quantum non-Gaussianity transferred by a photon-phonon-photon channel is feasible with both platforms. The required efficiencies, particularly ones of the readout, are significantly lower for the electromechanics than that needed for the transfer of negativity of Wigner function. Note that for the weakly coupled optomechanical system, the transfer of negativity is impossible regardless of the efficiencies of the state swap. Due to the initial occupation of the mechanical modes a very efficient state upload is required.

Methods
For a rigorous analysis that accounts for the imperfections properly, we start with the linearised Heisenberg-Langevin equations for the optomechanical system with a strong resonantly red-detuned pump. In terms of annihilation operators a c and a m of optical and mechanical modes respectively the equations read 27 : Here ω m is the mechanical mode frequency, γ is mechanical viscous damping coefficient, a in and a vac are the input and auxiliary vacuum optical modes respectively, and a th describes the Markovian noises from the mechanical environment. The coupling equals = ≡ g g g N k k 0 for k-th pulse and = g 0 during the storage between the pulses. Usually the condition of resolved sideband κ ω  m holds and therefore the rotating wave approximation (RWA) can be applied thus eliminating the rapid terms proportional to ω i t exp (2 ) m in the equations above. To prove validity of our results obtained with help of RWA, on Fig. 2 we present a solution that takes into account the rapid terms by the perturbation method. The fact that this complete solution coincides with the solution taking advantage of RWA justifies use of the latter for the rest of the paper.
With help of the input-output relation  The pulse shapes f t ( ) in and f t ( ) out are determined from the Eqs (9) and (10). Generally expressions for these functions are quite complicated hence we do not write those explicitly here. In the limit of high κ γ  g, , where the cavity mode can be approximately eliminated 32 ,  Fig. 2(b). The negativity-and entanglement-preserving thresholds V N neg and V N ent can be easily derived as follows. The value of V N neg is naturally obtained from the value of the Wigner function in the origin of the single-photon that passed the beamsplitter. The V N ent is derived from taking the maximum over the two-mode squeezed vacuum in a thermal channel. These boundaries are derived carefully in ref. 37.
The effective transmittivities T 1 and T 2 associated with respectively upload and readout stages present a convenient way to parametrize this dependence. Each of these transmittivities can be increased via increase of either corresponding pump power (intracavity photon number) or pulse duration. The latter way is accompanied by increase of the mechanical environment impact. In the Fig. 2(b) we assume the coupling close to the highest reported and change τ 1,2 to produce change of T 1,2 . We use the following dimensionless parameters to produce the plots: for the green lines (as reported in ref.  Evaluation of the transfer of quantum nongaussianity requires computation of the probability of vacuum and of single photon. If the restored optical quantum state is ρ f , these probabilities equal ρ = 〈 | | 〉 p 0 0 f 0 and ρ = 〈 | | 〉 p 1 1 f 1 and can be determined by homodyne tomography 33 . The quantum state ρ f is found as a result of a single photon |1〉 passing a virtual beamsplitter with transmittance  that adds noise with variance V N . The curves in Fig. 3a are produced for same numerical parameters as Fig. 2(b). To produce Fig. 3 , η = . 0 5 and η = . 0 25 respectively for yellow and purple lines.