Optically levitated nanoparticle as a model system for stochastic bistable dynamics F. Ricci,1 R. A. Rica,1 M. Spasenovi´c,1 J. Gieseler,1, 2 L. Rondin,3 L. Novotny,3 and R. Quidant1, 4 1ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, 08860 Castelldefels (Barcelona), Spain 2Physics Department, Harvard University, Cambridge, Massachusetts 02138, US. 3ETH Z¨urich,Photonics Laboratory, 8093 Z¨urich,Switzerland 4ICREA-Instituci´oCatalana de Recerca i Estudis Avancats, 08010 Barcelona, Spain Nano-mechanical resonators have gained an up to ∼ 50 dB amplification of non-resonant harmonic increasing importance in nanotechnology owing force at the atto-newton scale. In addition, the unprece- to their contributions to both fundamental and dented level of control achieved will enable the use of applied science. Yet, their small dimensions and levitated nanoparticles as a model system for stochastic mass raises some challenges as their dynamics bistable dynamics with applications to a wide range gets dominated by nonlinearities that degrade of fields including biophysics16,17, chemistry18,19 and their performance, for instance in sensing appli- nanotechnology20. cations. Here, we report on the precise control of the nonlinear and stochastic bistable dynamics Results of a levitated nanoparticle in high vacuum. Our experiment is sketched in Fig. 1a. A single sil- We demonstrate how it can lead to efficient ica nanoparticle (d ∼ 177 nm in diameter) is optically signal amplification schemes, including stochastic trapped in vacuum by a tightly focused laser beam, and resonance. This work paves the way for the use its 3D trajectory monitored with a balanced split detec- of levitated nanoparticles as a model system for tion scheme. The detector signal controls an electro-optic stochastic bistable dynamics, with applications modulator (EOM) that feedback-cools all three transla- to a wide variety of fields. tional degrees of freedom by modulating the laser in- tensity. In high-vacuum conditions (P ∼ 10−6 mbar) Introduction the mechanical Q-factor, given by the gas damping γ0, Nano-mechanical resonators need to meet criteria of is measured to be Q = 1:2 × 108. Under feedback, light mass and high-Q factor in order to maximize their the nanoparticle oscillates at non-degenerate angular fre- performances when operated as linear force transducers. (x;y;z) quencies Ω0 with very small amplitudes, correspond- However, these features lead to intrinsically nonlinear 11 behaviour1, with consequent vanishing dynamical range. ing to sub-Kelvin effective temperatures . By optimiz- To overcome this limitation, modern nanotechnology ing the feedback settings (Supplementary Note 2) and requires new sensing schemes that take nonlinearities carefully screening important sources of noise such as into account and even benefit from them2. Many of the mechanical vibrations and air turbulences, we obtain a proposed solutions operate inside an instability region3,4 highly stable resonator with frequency fluctuations im- 5,6 proved by one to two orders of magnitude compared to or close to a bifurcation point , where the system 1,21 ideally becomes infinitely sensitive. Others exploit previous works (See Supplementary Note 3). When one of the spatial modes is parametrically driven at res- fluctuations of noisy environments to trigger stochastic 21 resonances7 that amplify weak harmonic signals8{10. onance , the particle explores the anharmonic part of the optical potential, which can be modelled as a Duffing For all these sensing applications, a single nanopartticle 1 optimally decoupled from the environment represents a nonlinearity . As a result, the equation of motion for the particularly interesting system in the family of high-Q driven coordinate, which we choose to be x, reads resonators. In fact, since its very first realization11, optical levitation of a nano-scale object in vacuum has 2 2 2 F x¨ + γ0 + Ω0ηx x_ + Ω0 1 + ξx + cos (Ωmt) x = enabled several ground-breaking experiments, including m (1) demonstration of zeptonewton force sensitivity1,12, tests arXiv:1705.04061v1 [physics.optics] 11 May 2017 where η is the non-linear feedback-induced damping co- of fluctuation theorems and stochastic thermodynam- efficient, ξ is the Duffing term prefactor, and Ω are, ics13,14, as well as the observation of photon recoil m respectively, the modulation depth and modulation fre- heating15. quency of the driving signal. Finally, the fluctuating In this work, we demonstrate full control on the force F = F + F has two contributions. The first linear, non-linear and bistable dynamics of a levitated th noise one, F , represents the stochastic force arising from ran- nanoparticle in high vacuum and under the effect of th dom collisions with residual air molecules in the chamber, external noise. The potential of our platform is validated while the second, F , represents artificial parametric by the implementation of two nonlinear amplification noise noise that we add through the EOM. See Supplementary schemes, including stochastic resonance7. Remarkably, Note 1 for further experimental set-up details. we demonstrate, in excellent agreement with theory, The nonlinear response of the resonator is fully char- 2 a d Down-sweep Z ) Noise Driving Feedback Up-sweep 2 ) 2 Fitted amplitude 8000 Σ Y 8000 (nm 2 (nm y t 6000 i 2 l Adder Σ i X b 4000 a 4000 t mode s Vacuum i 2000 Chamber B Amplitude x 0 Amplitude x 0 Split y x Detection z b c Amplitude Effective Potential Uopt)H = 85 nm A Δ Uopt)H = 70 nm 0.08 A Δ = 55 nm OR Δ x x (nm) A A Modulation Depthϵ = 40 nm high 0.06 Δ High amplitude A Alow state 50 0.04 0.04 t 0.02 x 0 0 Low amplitude 0.02 -0.02 -50 m state -0.04 0 -0.06 Detuningδ Figure 1. Experimental set-up and nonlinear response. (a) A microscope objective (OBJ) focuses a laser beam inside a vacuum chamber, where a single silica nanoparticle is trapped. Its motion is measured with a split detection scheme and parametric feedback is applied via an electro-optical modulator (EOM) in order to cool its center of mass motion. (b) A linear low amplitude and a nonlinear high amplitude oscillation states coexist when bistability is induced by the coherent driving. (c) The bistable dynamics can be modelled by an effective amplitude double-well potential. Injected optical noise activates stochastic switching between the two states. (d) (top) Measured 2D false color map representing the particle amplitude response in the driving parameter phase space (δm; ). The red and blue cross-cuts, correspond respectively to down- and up- sweeps, and emphasize the typical hysteretic non-linear response of the resonator. Yellow lines are fits to the high amplitude solution A(δm) and allow to retrieve the nonlinear coefficients η and ξ. The bistable regime (bottom map), can be isolated subtracting up-sweep from down-sweep maps. Solid lines correspond to iso-amplitude lines (2): subsets of the parameter space where Ahigh, and consequently A∆ = Ahigh − Alow, are kept constant. acterized by measuring the oscillation amplitude in the provided in Supplementary Note 4. driving parameters phase space (δm; ), where δm = Inside the bistable region, natural thermal noise can 21 Ωm=Ω0 −2 is the normalized detuning . When the mod- activate spontaneous transitions between low and high ulation frequency is decreased (down-sweep) across the amplitude states. However, in high vacuum these events resonant condition δm = 0, the amplitude response spans are extremely rare, which makes our system a promis- a triangular region corresponding to the first instability ing candidate for nonlinear sensing. In order to observe tongue22. This is shown in the top 2D false color map stochastic activation within a reasonable time we there- of Fig. 1d. In the opposite case of a frequency up-sweep, fore add to the EOM a Gaussian noise ζ(t) of amplitude however, the system displays hysteresis, and the tongue NV , that in turns produces a position dependent para- results in a narrower width. The red and blue cross metric stochastic force of the form Fnoise / NV ζ(t)x. cuts in Fig. 1d clearly show this behaviour, identifying Noise calibration can be performed by looking at the par- a region where two stable oscillation states coexist (see ticle's oscillation amplitude A in the absence of paramet- Fig. 1b,c). The bistability region, shown at the bottom ric driving23. For low noise levels, the particle's dynamics of Fig. 1d, is obtained by subtracting up- from down- remains unaffected, being still dominated by the thermo- sweeps maps. By fitting our experimental data with the mechanical noise of the environment. However, for suffi- analytical solution A(δm; ) (See Supplementary Equa- ciently high noise, the particle's amplitude starts to in- 21 2 2 tion 1 and Ref. ) we retrieve the nonlinear coefficients mΩ0A crease. Assuming an effective temperature Teff = , η = (18:1±0:4) µm−2 and ξ = (−9:68±0:15) µm−2. Such kB the noise amplitudes NV are calibrated into meaning- precision measurements are made possible by the highly ful temperature units TN (See Supplementary Note 5). stable oscillation frequencies that lead to minimal drifts Interestingly, the measured trend differs from what re- < −3 of the detuning parameter (∆δm ≈ 2∆Ω=Ω0 ∼ 10 , over ported in other systems8, where a quadratic dependence the over the entire measurement time), therefore allow- between TN and NV was found. This is to be ascribed ing to predict with excellent agreement the boundaries of to the parametric (i.e. multiplicative) nature of the noise the instability