Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Demonstration of optical spring in an un-detuned cavity containing an optical parametric amplifier

Jian Liu ARC Centre of Excellence for Gravitational Wave Discovery, Department of Physics, University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia Juntao Pan ARC Centre of Excellence for Gravitational Wave Discovery, Department of Physics, University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia Carl Blair ARC Centre of Excellence for Gravitational Wave Discovery, Department of Physics, University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia Jue Zhang ARC Centre of Excellence for Gravitational Wave Discovery, Centre of Gravitational Astrophysics, Australian National University, Canberra ACT 2601, Australia Hengxin Sun State Key Laboratory of Quantum Optics and Quantum Optics Devices, Institute of Opto-Electronics, Shanxi University, Taiyuan 030006, China Li Ju ARC Centre of Excellence for Gravitational Wave Discovery, Department of Physics, University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia Chunnong Zhao ARC Centre of Excellence for Gravitational Wave Discovery, Department of Physics, University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia chunnong.zhao@uwa.edu.au
Abstract

Here we demonstrate the capacity to manipulate the optical spring (OS) effect by employing an optical parametric amplifier (OPA) within an optical cavity. We observed more than a factor of 2 increase in the OS frequency shift with the OPA. We also showed for the first time that the OS can be tuned by solely adjusting the OPA phase and showing an un-detuned cavity exhibiting an optical spring. The method can be applied to gravitational wave detectors in the signal recycling configuration to realize narrow bandwidth high sensitivity. The OS can be tuned to align the detector peak sensitivity frequency to known frequency continuous gravitational wave signals, dynamically tuned to track the gravitational wave signal from merging compact binaries or tuned to search for the post-merger signal of known binary coalescence.

journal: opticajournal

1 Introduction

Currently operating second-generation ground-based gravitational wave detectors (GWDs) including Advanced LIGO[1], Advanced Virgo[2], and KAGRA[3] are now being pushed towards the designed sensitivity where quantum noise is dominating the entire band from tens of hertz up to a few kilohertz. Though ’classical’ ways such as increasing the laser power and increasing the test masses weight can shape the quantum noise limited sensitivity in the GWDs, the sensitivity remains limited by the free mass standard quantum limit (SQL). Various ’non-classical’ quantum techniques have been proposed and implemented in GWDs to tackle quantum noise. A sensitivity enhancement of more than 6 dB below the shot noise limit has been demonstrated with external phase squeezed vacuum states injection[4, 5]. Further with an external optical filter cavity or interferometer itself via EPR entanglement[6] to rotate the squeezing angle depending on frequency, a technique called frequency-dependent squeezing[7, 8], the broadband quantum noise reduction has been achieved[9].

Alternatively, the optomechanical interactions between the optical field and suspended test masses inside the GWDs can correlate the phase and amplitude quadratures of the optical field, resulting in pondermotive squeezing that can greatly improve the detector sensitivity to below the SQL[10]. This was analyzed in detail by Buonanno and Chen in the context of GWDs with a blue-detuned signal recycle cavity (SRC)[11] where the detector can be considered an ’optical spring’ detector where the free test mass is coupled with an optical spring. The coupled system can resonate at a pair of frequencies, namely the OS frequency and optical resonance frequency. At these frequencies, the detector quantum noise limited sensitivity is significantly enhanced. The OS frequency is normally around a few to tens of hertz, depending on the SRC detuning and the arm cavity power. Increasing the arm cavity power can shift the OS frequency higher but brings undesired effects[12, 13, 14]. It is also worthwhile to note that along with the extra optical rigidity, the optomechanical interaction also results in optical heating that excites the mechanical motion of the suspended test mass[15], external feedback control or optical damping will be needed[16, 17].

The idea of optical-spring-enhanced detectors was further extended to the dynamic tuning of optical springs by adjusting the SRC parameters, allowing tracking of expected chirps of GW signals[18], but due to complexity in controls, this was not yet employed in GWDs. Somiya etal.𝑒𝑡𝑎𝑙et\leavevmode\nobreak\ al.italic_e italic_t italic_a italic_l . proposed a new scheme with OPA inside the SRC that allows shifting the OS higher without increasing the arm cavity power[19]. Korobko etal.𝑒𝑡𝑎𝑙et\leavevmode\nobreak\ al.italic_e italic_t italic_a italic_l . further found that the OPA phase was crucial for shaping the interferometer’s response[20]. Recently, Somiya etal.𝑒𝑡𝑎𝑙et\leavevmode\nobreak\ al.italic_e italic_t italic_a italic_l . expanded the idea in the context of a dual recycled interferometer without arm cavities, they found that by properly choosing parameters, the OS frequency can be shifted to kilohertz, allowing detections of GW signals from the postmerger remnant of binary neutron star collisions[21]. The experiment carried out by Zhang etal.𝑒𝑡𝑎𝑙et\leavevmode\nobreak\ al.italic_e italic_t italic_a italic_l . confirmed that the OS can be enhanced by OPA in a detuned linear cavity with a strong carrier[22]. In this letter, we demonstrated that both the OPA gain and phase can be used to manipulate the OS effect.

Our model is that of a linear cavity with a moveable end mirror and a nonlinear crystal inside. A 1064 nm laser beam together with a 532 nm beam is injected into the cavity, the infrared beam resonates inside the cavity while the green beam doesn’t. By adopting Huang etal.𝑒𝑡𝑎𝑙et\leavevmode\nobreak\ al.italic_e italic_t italic_a italic_l .’s formalism[23], we can write the system Hamiltonian as

H=𝐻absent\displaystyle H=italic_H = (ωcωL)c^c^ωcLq^c^c^+12(p^2m+mωm2q^2)Planck-constant-over-2-pisubscript𝜔𝑐subscript𝜔𝐿superscript^𝑐^𝑐Planck-constant-over-2-pisubscript𝜔𝑐𝐿^𝑞superscript^𝑐^𝑐12superscript^𝑝2𝑚𝑚superscriptsubscript𝜔𝑚2superscript^𝑞2\displaystyle\hbar\left(\omega_{c}-\omega_{L}\right)\hat{c}^{\dagger}\hat{c}-% \frac{\hbar\omega_{c}}{L}\hat{q}\hat{c}^{\dagger}\hat{c}+\frac{1}{2}\left(% \frac{\hat{p}^{2}}{m}+m\omega_{m}^{2}\hat{q}^{2}\right)roman_ℏ ( italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - italic_ω start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG - divide start_ARG roman_ℏ italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG over^ start_ARG italic_q end_ARG over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( divide start_ARG over^ start_ARG italic_p end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m end_ARG + italic_m italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG italic_q end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) (1)
+iε(c^c^)+iG(eiθc^2eiθc^2),𝑖Planck-constant-over-2-pi𝜀superscript^𝑐^𝑐𝑖Planck-constant-over-2-pi𝐺superscript𝑒𝑖𝜃superscript^𝑐absent2superscript𝑒𝑖𝜃superscript^𝑐2\displaystyle+i\hbar\varepsilon\left(\hat{c}^{\dagger}-\hat{c}\right)+i\hbar G% \left(e^{i\theta}\hat{c}^{\dagger 2}-e^{-i\theta}\hat{c}^{2}\right),+ italic_i roman_ℏ italic_ε ( over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT - over^ start_ARG italic_c end_ARG ) + italic_i roman_ℏ italic_G ( italic_e start_POSTSUPERSCRIPT italic_i italic_θ end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † 2 end_POSTSUPERSCRIPT - italic_e start_POSTSUPERSCRIPT - italic_i italic_θ end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ,

where ωLsubscript𝜔𝐿\omega_{L}italic_ω start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT and ωcsubscript𝜔𝑐\omega_{c}italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT represent the frequency of the laser and the resonant frequency of the cavity, c^^𝑐\hat{c}over^ start_ARG italic_c end_ARG and c^superscript^𝑐\hat{c}^{\dagger}over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT are the annihilation and creation operators of the cavity field; p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG and q^^𝑞\hat{q}over^ start_ARG italic_q end_ARG represent the momentum and position operators for the mechanical resonator; ε=2γP/(ωL)𝜀2𝛾𝑃Planck-constant-over-2-pisubscript𝜔𝐿\varepsilon=\sqrt{2\gamma}\sqrt{P/\left(\hbar\omega_{L}\right)}italic_ε = square-root start_ARG 2 italic_γ end_ARG square-root start_ARG italic_P / ( roman_ℏ italic_ω start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) end_ARG is the cavity field driving strength with input laser power P𝑃Pitalic_P; γ=πc/(2L)𝛾𝜋𝑐2𝐿\gamma=\pi c/(2\mathcal{F}L)italic_γ = italic_π italic_c / ( 2 caligraphic_F italic_L ) is half-linewidth of the cavity with finesse \mathcal{F}caligraphic_F and cavity length L𝐿Litalic_L, G𝐺Gitalic_G is the non-linear gain of the OPA which is proportional to the green pump power, and θ𝜃\thetaitalic_θ is the phase of the pump field.

We can then solve the equations of motion and obtain the steady-state values

ps=0,qs=ωc|cs|2mLωm2,cs=|cs|eiϕ=γiΔ+2Geiθγ2+Δ24G2ε,formulae-sequencesubscript𝑝𝑠0formulae-sequencesubscript𝑞𝑠Planck-constant-over-2-pisubscript𝜔𝑐superscriptsubscript𝑐𝑠2𝑚𝐿superscriptsubscript𝜔𝑚2subscript𝑐𝑠subscript𝑐𝑠superscript𝑒𝑖italic-ϕ𝛾𝑖Δ2𝐺superscript𝑒𝑖𝜃superscript𝛾2superscriptΔ24superscript𝐺2𝜀p_{s}=0,q_{s}=\frac{\hbar\omega_{c}\left|c_{s}\right|^{2}}{mL\omega_{m}^{2}},c% _{s}=\lvert c_{s}\rvert e^{i\phi}=\frac{\gamma-i\Delta+2Ge^{i\theta}}{\gamma^{% 2}+\Delta^{2}-4G^{2}}\varepsilon,italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0 , italic_q start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = divide start_ARG roman_ℏ italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m italic_L italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | italic_e start_POSTSUPERSCRIPT italic_i italic_ϕ end_POSTSUPERSCRIPT = divide start_ARG italic_γ - italic_i roman_Δ + 2 italic_G italic_e start_POSTSUPERSCRIPT italic_i italic_θ end_POSTSUPERSCRIPT end_ARG start_ARG italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_ε , (2)

where

Δ=Δ0(ωc)2|cs|2mL2ωm2,ϕ=arctan2Gsin[θ]Δγ+2Gcos[θ]formulae-sequenceΔsubscriptΔ0Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐2superscriptsubscript𝑐𝑠2𝑚superscript𝐿2superscriptsubscript𝜔𝑚2italic-ϕ2𝐺𝜃Δ𝛾2𝐺𝜃\Delta=\Delta_{0}-\frac{\hbar\left(\omega_{c}\right)^{2}\left|c_{s}\right|^{2}% }{mL^{2}\omega_{m}^{2}},\phi=\arctan{\frac{2G\sin[\theta]-\Delta}{\gamma+2G% \cos[\theta]}}roman_Δ = roman_Δ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - divide start_ARG roman_ℏ ( italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , italic_ϕ = roman_arctan divide start_ARG 2 italic_G roman_sin [ italic_θ ] - roman_Δ end_ARG start_ARG italic_γ + 2 italic_G roman_cos [ italic_θ ] end_ARG (3)

are the effective cavity detuning, and the phase of the cavity field, respectively. We can further solve the equations of motion in the frequency domain and get the modified mechanical susceptibility

χxx1=m(ω2+i[γm+Γos]ω[ωm2+Kos/m]),superscriptsubscript𝜒𝑥𝑥1𝑚superscript𝜔2𝑖delimited-[]subscript𝛾msubscriptΓos𝜔delimited-[]superscriptsubscript𝜔m2subscript𝐾os𝑚\chi_{xx}^{-1}=-m(\omega^{2}+i[\gamma_{\rm{m}}+\Gamma_{\rm{os}}]\omega-[\omega% _{\rm{m}}^{2}+K_{\rm{os}}/m]),italic_χ start_POSTSUBSCRIPT italic_x italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_m ( italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_i [ italic_γ start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT ] italic_ω - [ italic_ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT / italic_m ] ) , (4)

where the optical spring constant Kossubscript𝐾osK_{\rm{os}}italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT and optical damping ΓossubscriptΓos\Gamma_{\rm{os}}roman_Γ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT are defined as:

Kossubscript𝐾os\displaystyle K_{\rm{os}}italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT =Re[χos1]absent𝑅𝑒delimited-[]superscriptsubscript𝜒os1\displaystyle=Re[\chi_{\rm{os}}^{-1}]= italic_R italic_e [ italic_χ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ] (5)
=2(ωcL)2|cs|2(γ2+Δ24G2ω2)(Δ+2Gsin[θ2ϕ])(γ2+Δ24G2ω2)2+4γ2ω2,absent2Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐𝐿2superscriptsubscript𝑐𝑠2superscript𝛾2superscriptΔ24superscript𝐺2superscript𝜔2Δ2𝐺𝜃2italic-ϕsuperscriptsuperscript𝛾2superscriptΔ24superscript𝐺2superscript𝜔224superscript𝛾2superscript𝜔2\displaystyle=\frac{2\hbar(\frac{\omega_{c}}{L})^{2}\lvert c_{s}\rvert^{2}(% \gamma^{2}+\Delta^{2}-4G^{2}-\omega^{2})(\Delta+2G\sin[\theta-2\phi])}{(\gamma% ^{2}+\Delta^{2}-4G^{2}-\omega^{2})^{2}+4\gamma^{2}\omega^{2}},= divide start_ARG 2 roman_ℏ ( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ( roman_Δ + 2 italic_G roman_sin [ italic_θ - 2 italic_ϕ ] ) end_ARG start_ARG ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,
Γos=Im[χos1]meffω=4(ωcL)2|cs|2γ(Δ+2Gsin[θ2ϕ])meff((γ2+Δ24G2ω2)2+4γ2ω2).subscriptΓos𝐼𝑚delimited-[]superscriptsubscript𝜒os1subscript𝑚eff𝜔4Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐𝐿2superscriptsubscript𝑐𝑠2𝛾Δ2𝐺𝜃2italic-ϕsubscript𝑚effsuperscriptsuperscript𝛾2superscriptΔ24superscript𝐺2superscript𝜔224superscript𝛾2superscript𝜔2\Gamma_{\rm{os}}=\frac{Im[\chi_{\rm{os}}^{-1}]}{m_{\rm{eff}}\omega}=\frac{4% \hbar(\frac{\omega_{c}}{L})^{2}\lvert c_{s}\rvert^{2}\ \gamma(\Delta+2G\sin[% \theta-2\phi])}{m_{\rm{eff}}((\gamma^{2}+\Delta^{2}-4G^{2}-\omega^{2})^{2}+4% \gamma^{2}\omega^{2})}.roman_Γ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT = divide start_ARG italic_I italic_m [ italic_χ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ] end_ARG start_ARG italic_m start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT italic_ω end_ARG = divide start_ARG 4 roman_ℏ ( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_γ ( roman_Δ + 2 italic_G roman_sin [ italic_θ - 2 italic_ϕ ] ) end_ARG start_ARG italic_m start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT ( ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG . (6)

It can be seen from Eq.(5,6) that the mechanical frequency ωmsubscript𝜔𝑚\omega_{m}italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT and damping rate γmsubscript𝛾𝑚\gamma_{m}italic_γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT are modified, by both the OPA gain and phase. The new effective mechanical frequency and damping rate are ωeff=ωm2+Kos/msubscript𝜔effsuperscriptsubscript𝜔m2subscript𝐾os𝑚\omega_{\rm{eff}}=\sqrt{\omega_{\rm{m}}^{2}+K_{\rm{os}}/m}italic_ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT = square-root start_ARG italic_ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT / italic_m end_ARG and γeff=γm+Γossubscript𝛾effsubscript𝛾msubscriptΓos\gamma_{\rm{eff}}=\gamma_{\rm{m}}+\Gamma_{\rm{os}}italic_γ start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT = italic_γ start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT. Without OPA, the results reduce to the normal optomechanical cavity case[15]. A more detailed calculation can be found in supplementary 5.1.

2 Experiment setup

The schematic of the experiment is presented in Fig.1. A 2 W Nd:YAG laser produces a linear polarized 1064 nm beam and it is split into two parts. Half of the beam goes to a second harmonic generator (SHG). The SHG is locked via the standard PDH locking technique[24] and can produce a 532 nm green beam with a power of about 500 mW. The other half of the 1064 beam is further split via a polarized beam splitter (PBS) into two paths: the S-polarized locking beam and the P-polarized seed beam.

The locking beam goes through two acoustic-optical modulators (AOM) to shift the frequency relative to the seed beam for cavity detuning purposes. It combines with the seed beam at a second PBS. A half-wave plate is located after this PBS to adjust the beam polarization. The 1064 nm infrared beams are further combined with the S-polarized 532 nm green beam at a dichroic mirror that is HR-coated for 1064 nm and AR-coated for 532 nm. Both the seed beam and locking beam are mode-matched to the cavity and the green beam is focused down to below 50 um in radius to maximize the intensity on the crystal. The whole cavity sits in a high vacuum chamber that has a pressure of about 1e-5 mbar to reduce the air-damping.

Refer to caption
Figure 1: Schematic of the experiment. 1064 nm EOMs are driven at 12 MHz and 532 nm EOM is driven at 300 kHz. All three beams have independent power control allowing arbitrary relative power injection and they are also independently monitored at the transmission port.

The cavity is a Z-shape linear cavity with three half-inch mirrors and a micro-resonator as the end mirror. The total cavity length is 260 mm with a 50-um waist on the input mirror and a second 150-um waist on the resonator. The input mirror is AR coated at 1064 nm and 532 nm at the outside surface while the inner surface is AR coated for 532 nm but HR coated with 99.8%percent\%% power reflectivity for 1064 nm. A 1×\times×2×\times×10 mm PPKTP crystal is placed right after the input mirror, and both the 1×\times×2 mm surfaces are AR coated with power reflectivity less than 0.1%percent\%%. The second cavity mirror is a concave mirror with a radius of curvature of 100 mm and the third one is a flat mirror. Both of them have a power reflectivity of 99.98%percent\%% at 1064 nm and less than 5%percent\%% at 532 nm. By reaching the micro-resonator, the green power would be attenuated by more than three orders of magnitude, so that the thermal effect induced by the green beam absorption on the resonator is negligible. The transmitted beams from the curved mirror are split via a dichroic mirror and then through a PBS depending on wavelength and polarization, then reach photodiodes for power monitoring.

The resonator is made of a 420 μ𝜇\muitalic_μm thick 20 mm wide square silicon frame that is etched with multiple through holes of different sizes leaving windows of a multi-layer AlGaAs/GaAs coating bonded on one side. Each of the coating windows can be used as an individual resonator with different effective masses, resonant frequencies, and quality factors. These windows are located within the central 7 mm ×\times× 7 mm area of the frame. To minimize the loss from clamping, the whole frame is clamped using a 9 mm radius ring shape holder. The resonator used in this experiment is 0.64 mm ×\times× 0.64 mm, the drum shape mode frequency is 85 kHz with an effective mass of 1.75 μ𝜇\muitalic_μg and a quality factor of about 2500 measured in vacuum (see details in supplementary 5.2). The cavity linewidth of the seed beam was measured to be 1.14 MHz, corresponding to a cavity finesse of around =500500\mathcal{F}=500caligraphic_F = 500 for the seed beam.

3 Results

We first measured the optical spring effect without the green pump beam injection. Only a 3 mW locking beam together with a 10 mW seed beam is injected into the cavity. The cavity is stably locked to the locking beam and the seed beam is detuned from the cavity resonance by adjusting the driving frequency of AOM2. The OPA cavity length locking bandwidth is about 3 kHz and the loop gain around the mechanical mode frequency is about -80 dB, so the loop should have little effect on the mechanical mode. The mode frequency and linewidth can be obtained through the OPA cavity length locking error signal via a spectrum analyzer.

Refer to caption
Figure 2: Examples of the 85 kHz drum mode power spectrum density under blue, red, and no detuning and the corresponding Lorenzian fittings.

Examples of the measurement of mechanical mode power spectrum density from a spectrum analyzer are shown in Fig.2. Such curves are measured over a range of cavity detuning of [3γ3𝛾-3\leavevmode\nobreak\ \gamma- 3 italic_γ, 3γ3𝛾-3\leavevmode\nobreak\ \gamma- 3 italic_γ]. Each of the measured curves is then fitted with a Lorentzian shape to obtain the mode central frequency and the linewidth. The OS effect without OPA is summarized in Fig.3. The maximum frequency shift is about ±plus-or-minus\pm±200 Hz and the maximum damping or anti-damping is about ±plus-or-minus\pm±100 Hz at around ±0.25γplus-or-minus0.25𝛾\pm 0.25\leavevmode\nobreak\ \gamma± 0.25 italic_γ detuning. The measurements agree with the theory quite well. In the theory curve, we have taken into account the 20% power loss through mode mismatch, the loss from the vacuum tank window, and 0.65 for the overlap factor between the cavity fundamental mode and the resonator drum mode. Each measurement curve takes about 20 seconds to average the noise. During the whole measurement, the natural frequency of the mode can shift tens of hertz due to the thermal effect and the mode linewidth can increase due to pressure change in the vacuum tank. The detuning has to be set to zero after a few measurements to allow calibration of the mode’s natural frequency and linewidth to reduce errors in the OS results.

Refer to caption
Figure 3: Optical spring effect without OPA. The seed beam power is 10 mW and the transmission PD reading is 2.6 V.

In the following step, we then injected about 50 mW p-polarized green pump laser into the cavity, the pump laser is phase modulated by driving a 532 nm EOM at 300 kHz. The PPKTP crystal temperature is stabilized at 34 degrees to meet the phase-matching condition of the OPA to maximize the down-conversion process. The seed laser 1064 nm beam will be amplified inside the cavity and is monitored through transmission. The transmission signal is then mixed with the 300 kHz local oscillator to generate the phase-locking error signal. Depending on the relative phase between the pump beam and the seed beam, the seed beam power can be either amplified or de-amplified. We define the OPA power gain as

gOPA=|cs(G)|Δ0,θ0cs(G=0)|Δ0,θ0|2=(γγ2G)2.subscript𝑔OPAsuperscriptevaluated-atsubscript𝑐𝑠𝐺formulae-sequenceΔ0𝜃0evaluated-atsubscript𝑐𝑠𝐺0formulae-sequenceΔ0𝜃02superscript𝛾𝛾2𝐺2g_{\rm{OPA}}=\left|\frac{c_{s}(G)|_{\Delta\rightarrow 0,\theta\rightarrow 0}}{% c_{s}(G=0)|_{\Delta\rightarrow 0,\theta\rightarrow 0}}\right|^{2}=\left(\frac{% \gamma}{\gamma-2G}\right)^{2}.italic_g start_POSTSUBSCRIPT roman_OPA end_POSTSUBSCRIPT = | divide start_ARG italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_G ) | start_POSTSUBSCRIPT roman_Δ → 0 , italic_θ → 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_G = 0 ) | start_POSTSUBSCRIPT roman_Δ → 0 , italic_θ → 0 end_POSTSUBSCRIPT end_ARG | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( divide start_ARG italic_γ end_ARG start_ARG italic_γ - 2 italic_G end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (7)

Without the green pump injection, 0.5 mW seed laser registered 0.2 V in the transmission PD, while with 50 mW pump beam and OPA phase locked in amplification, the transmission read 2.0 V, meaning that the gOPAsubscript𝑔OPAg_{\rm OPA}italic_g start_POSTSUBSCRIPT roman_OPA end_POSTSUBSCRIPT = 10, corresponding to a non-linear gain G𝐺Gitalic_G= 0.2 MHz.

We then repeated the OS measurements while keeping the pump beam power stable and the OPA phase locked to maximum amplification. By adjusting the AOM2 driving frequency, we were able to detune the cavity and obtain the OPA-enhanced OS results as shown in Fig.4. The peak frequency shift was enhanced by a factor of 2 and the damping(or anti-damping) was enhanced by a factor of 4.

Refer to caption
Figure 4: OS as a function of the cavity detuning with and without OPA. Note that the x-axis γ𝛾\gammaitalic_γ refers to the original cavity half-linewidth without OPA. The effective cavity linewidth is less with OPA. The w/o OPA traces are from Fig.3 by normalizing the intra-cavity power from 2.6 V to 2.0 V on transmission PD.

As shown in Eq.(5,6), the OPA enhancement is not only achieved by tuning the OPA gain G𝐺Gitalic_G but also the OPA phase θ𝜃\thetaitalic_θ. In our next measurement, we kept the Δ=0Δ0\Delta=0roman_Δ = 0 by adjusting the AOM2. The OPA phase, instead of being locked on θ=0𝜃0\theta=0italic_θ = 0 that has maximum amplification, is locked at an arbitrary value by adding a phase offset to the phase-locking loop error point. By comparing the offset with the full size of the error signal, the phase can be calculated (See supplementary 5.3). The intra-cavity power would change while tuning the OPA phase and ideally, the input power needs to be adjusted to maintain the intra-cavity power the same for a fair comparison. But within a range from 40superscript40-40^{\circ}- 40 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT to 40superscript4040^{\circ}40 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT, the intra-cavity power can be considered stable. The measurements are shown in Fig.5.

Refer to caption
Figure 5: OS as a function of the OPA phase θ𝜃\thetaitalic_θ. The blue trace has the intra-cavity power normalized to 1.8 V while the orange trace assumes the maximum power is 1.8 V and the intra-cavity power decreases while the phase changes. It can be seen that when the phase is small, the OS effects are the same for the two traces.

4 Conclusions

We have experimentally demonstrated that the OS effect can be enhanced by an OPA inside the cavity by more than a factor of two. Both the OPA gain and phase can be tuned independently or together to modify the OS effect. This method can be applied in large-scale gravitational wave detectors to improve detector sensitivity by shifting the OS frequency from a few hertz where low-frequency seismic noise dominates, to a few hundred hertz where the quantum noise dominates. The peak sensitivity frequency can be set to a certain value for monitoring gravitational wave signals that have known frequency or dynamically tuning the peak sensitivity frequency would allow real-time tracking of the gravitational events.

\bmsection

Funding This work was supported by the Australian Research Council

\bmsection

Acknowledgment The authors would like to thank

\bmsection

Disclosures The authors declare no conflicts of interest.

\bmsection

Data availability Data underlying the results presented in this paper are available from the corresponding author upon reasonable request.

\bmsection

Supplemental document See Supplements for supporting content.

References

  • [1] J. Aasi, B. Abbott, R. Abbott, et al., \JournalTitleClassical and quantum gravity 32, 074001 (2015).
  • [2] F. a. Acernese, M. Agathos, K. Agatsuma, et al., \JournalTitleClassical and Quantum Gravity 32, 024001 (2014).
  • [3] Y. Aso, Y. Michimura, K. Somiya, et al., \JournalTitlePhysical Review D 88, 043007 (2013).
  • [4] F. Acernese, M. Agathos, L. Aiello, et al., \JournalTitlePhysical review letters 123, 231108 (2019).
  • [5] M. e. Tse, H. Yu, N. Kijbunchoo, et al., \JournalTitlePhysical Review Letters 123, 231107 (2019).
  • [6] Y. Ma, H. Miao, B. H. Pang, et al., \JournalTitleNature Physics 13, 776 (2017).
  • [7] L. McCuller, C. Whittle, D. Ganapathy, et al., \JournalTitlePhysical review letters 124, 171102 (2020).
  • [8] F. Acernese, M. Agathos, A. Ain, et al., \JournalTitlePhysical review letters 131, 041403 (2023).
  • [9] D. Ganapathy, W. Jia, M. Nakano, et al., \JournalTitlePhysical Review X 13, 041021 (2023).
  • [10] V. Braginsky and F. Y. Khalili, \JournalTitlePhysics Letters A 257, 241 (1999).
  • [11] A. Buonanno and Y. Chen, \JournalTitlePhysical Review D 64, 042006 (2001).
  • [12] C. Zhao, J. Degallaix, L. Ju, et al., \JournalTitlePhysical review letters 96, 231101 (2006).
  • [13] M. Evans, S. Gras, P. Fritschel, et al., \JournalTitlePhysical review letters 114, 161102 (2015).
  • [14] J. Liu, V. Bossilkov, C. Blair, et al., \JournalTitleReview of Scientific Instruments 89 (2018).
  • [15] M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt, \JournalTitleReviews of Modern Physics 86, 1391 (2014).
  • [16] H. Rehbein, H. Müller-Ebhardt, K. Somiya, et al., \JournalTitlePhysical Review D 78, 062003 (2008).
  • [17] R. Singh, G. D. Cole, J. Cripe, and T. Corbitt, \JournalTitlePhysical Review Letters 117, 213604 (2016).
  • [18] B. J. Meers, A. Krolak, and J. A. Lobo, \JournalTitlePhysical Review D 47, 2184 (1993).
  • [19] K. Somiya, Y. Kataoka, J. Kato, et al., \JournalTitlePhysics Letters A 380, 521 (2016).
  • [20] M. Korobko, F. Y. Khalili, and R. Schnabel, \JournalTitlePhysics Letters A 382, 2238 (2018).
  • [21] O. S. Somiya K, Suzuki K and H. K, \JournalTitlePhysical Review D 107, 122005 (2023).
  • [22] J. Zhang, H. Sun, H. Guo, et al., \JournalTitleApplied Physics Letters 122 (2023).
  • [23] S. Huang and G. S. Agarwal, \JournalTitlePhysical Review A 79, 013821 (2009).
  • [24] E. D. Black, \JournalTitleAmerican journal of physics 69, 79 (2001).
\bibliographyfullrefs

sample

5 Supplementary materials

5.1 Theoretical model

The system Hamiltonian can be written as

H=𝐻absent\displaystyle H=italic_H = (ωcωL)c^c^ωcLq^c^c^+12(p^2m+mωm2q^2)Planck-constant-over-2-pisubscript𝜔𝑐subscript𝜔𝐿superscript^𝑐^𝑐Planck-constant-over-2-pisubscript𝜔𝑐𝐿^𝑞superscript^𝑐^𝑐12superscript^𝑝2𝑚𝑚superscriptsubscript𝜔𝑚2superscript^𝑞2\displaystyle\hbar\left(\omega_{c}-\omega_{L}\right)\hat{c}^{\dagger}\hat{c}-% \frac{\hbar\omega_{c}}{L}\hat{q}\hat{c}^{\dagger}\hat{c}+\frac{1}{2}\left(% \frac{\hat{p}^{2}}{m}+m\omega_{m}^{2}\hat{q}^{2}\right)roman_ℏ ( italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - italic_ω start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG - divide start_ARG roman_ℏ italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG over^ start_ARG italic_q end_ARG over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( divide start_ARG over^ start_ARG italic_p end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m end_ARG + italic_m italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG italic_q end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) (8)
+iε(c^c^)+iG(eiθc^2eiθc^2),𝑖Planck-constant-over-2-pi𝜀superscript^𝑐^𝑐𝑖Planck-constant-over-2-pi𝐺superscript𝑒𝑖𝜃superscript^𝑐absent2superscript𝑒𝑖𝜃superscript^𝑐2\displaystyle+i\hbar\varepsilon\left(\hat{c}^{\dagger}-\hat{c}\right)+i\hbar G% \left(e^{i\theta}\hat{c}^{\dagger 2}-e^{-i\theta}\hat{c}^{2}\right),+ italic_i roman_ℏ italic_ε ( over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT - over^ start_ARG italic_c end_ARG ) + italic_i roman_ℏ italic_G ( italic_e start_POSTSUPERSCRIPT italic_i italic_θ end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † 2 end_POSTSUPERSCRIPT - italic_e start_POSTSUPERSCRIPT - italic_i italic_θ end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ,

with each term corresponding to the intra-cavity optical field, the interaction between the intra-cavity field with the mechanical mode, the energy of the mechanical mode, the driving from the injection, and the down-conversion process in the crystal, respectively.

The equation of motion can be expressed as

q^˙=p^/mp^˙=mωm2q^+ωcLc^c^γmp^+ξ^c^˙=i(ωLωc)c^+iωcLq^c^+ε+2Geiθc^γc^+2γc^inc^˙=i(ωcωL)c^iωcLq^c^+ε+2Geiθc^γc^+2γc^in.˙^𝑞^𝑝𝑚˙^𝑝𝑚superscriptsubscript𝜔𝑚2^𝑞Planck-constant-over-2-pisubscript𝜔𝑐𝐿superscript^𝑐^𝑐subscript𝛾𝑚^𝑝^𝜉˙^𝑐𝑖subscript𝜔𝐿subscript𝜔𝑐^𝑐𝑖subscript𝜔𝑐𝐿^𝑞^𝑐𝜀2𝐺superscript𝑒𝑖𝜃superscript^𝑐𝛾^𝑐2𝛾subscript^𝑐insuperscript˙^𝑐𝑖subscript𝜔𝑐subscript𝜔𝐿superscript^𝑐𝑖subscript𝜔𝑐𝐿^𝑞^𝑐𝜀2𝐺superscript𝑒𝑖𝜃^𝑐𝛾superscript^𝑐2𝛾superscriptsubscript^𝑐in\begin{array}[]{l}\dot{\hat{q}}=\hat{p}/m\\ \dot{\hat{p}}=-m\omega_{m}^{2}\hat{q}+\hbar\frac{\omega_{c}}{L}\hat{c}^{% \dagger}\hat{c}-\gamma_{m}\hat{p}+\hat{\xi}\\ \dot{\hat{c}}=i\left(\omega_{L}-\omega_{c}\right)\hat{c}+i\frac{\omega_{c}}{L}% \hat{q}\hat{c}+\varepsilon+2Ge^{i\theta}\hat{c}^{\dagger}-\gamma\hat{c}+\sqrt{% 2\gamma}\hat{c}_{\rm{in}}\\ \dot{\hat{c}}^{\dagger}=i\left(\omega_{c}-\omega_{L}\right)\hat{c}^{\dagger}-i% \frac{\omega_{c}}{L}\hat{q}\hat{c}+\varepsilon+2Ge^{-i\theta}\hat{c}-\gamma% \hat{c}^{\dagger}+\sqrt{2\gamma}\hat{c}_{\rm{in}}^{\dagger}.\end{array}start_ARRAY start_ROW start_CELL over˙ start_ARG over^ start_ARG italic_q end_ARG end_ARG = over^ start_ARG italic_p end_ARG / italic_m end_CELL end_ROW start_ROW start_CELL over˙ start_ARG over^ start_ARG italic_p end_ARG end_ARG = - italic_m italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG italic_q end_ARG + roman_ℏ divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG - italic_γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT over^ start_ARG italic_p end_ARG + over^ start_ARG italic_ξ end_ARG end_CELL end_ROW start_ROW start_CELL over˙ start_ARG over^ start_ARG italic_c end_ARG end_ARG = italic_i ( italic_ω start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT - italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) over^ start_ARG italic_c end_ARG + italic_i divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG over^ start_ARG italic_q end_ARG over^ start_ARG italic_c end_ARG + italic_ε + 2 italic_G italic_e start_POSTSUPERSCRIPT italic_i italic_θ end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT - italic_γ over^ start_ARG italic_c end_ARG + square-root start_ARG 2 italic_γ end_ARG over^ start_ARG italic_c end_ARG start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL over˙ start_ARG over^ start_ARG italic_c end_ARG end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT = italic_i ( italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - italic_ω start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT - italic_i divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG over^ start_ARG italic_q end_ARG over^ start_ARG italic_c end_ARG + italic_ε + 2 italic_G italic_e start_POSTSUPERSCRIPT - italic_i italic_θ end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG - italic_γ over^ start_ARG italic_c end_ARG start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT + square-root start_ARG 2 italic_γ end_ARG over^ start_ARG italic_c end_ARG start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT . end_CELL end_ROW end_ARRAY (9)

Here we included ξ^^𝜉\hat{\xi}over^ start_ARG italic_ξ end_ARG as the Brownian noise operator and c^insubscript^𝑐in\hat{c}_{\rm{in}}over^ start_ARG italic_c end_ARG start_POSTSUBSCRIPT roman_in end_POSTSUBSCRIPT as the input vacuum noise operator both have a mean value of zero. The steady-state solution of the equation would be

ps=0,qs=ωc|cs|2mLωm2,cs=|cs|eiϕ=γiΔ+2Geiθγ2+Δ24G2ε,formulae-sequencesubscript𝑝𝑠0formulae-sequencesubscript𝑞𝑠Planck-constant-over-2-pisubscript𝜔𝑐superscriptsubscript𝑐𝑠2𝑚𝐿superscriptsubscript𝜔𝑚2subscript𝑐𝑠subscript𝑐𝑠superscript𝑒𝑖italic-ϕ𝛾𝑖Δ2𝐺superscript𝑒𝑖𝜃superscript𝛾2superscriptΔ24superscript𝐺2𝜀p_{s}=0,q_{s}=\frac{\hbar\omega_{c}\left|c_{s}\right|^{2}}{mL\omega_{m}^{2}},c% _{s}=\lvert c_{s}\rvert e^{i\phi}=\frac{\gamma-i\Delta+2Ge^{i\theta}}{\gamma^{% 2}+\Delta^{2}-4G^{2}}\varepsilon,italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0 , italic_q start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = divide start_ARG roman_ℏ italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m italic_L italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | italic_e start_POSTSUPERSCRIPT italic_i italic_ϕ end_POSTSUPERSCRIPT = divide start_ARG italic_γ - italic_i roman_Δ + 2 italic_G italic_e start_POSTSUPERSCRIPT italic_i italic_θ end_POSTSUPERSCRIPT end_ARG start_ARG italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_ε , (10)

where

Δ=Δ0(ωc)2|cs|2mL2ωm2,ϕ=arctan2Gsin[θ]Δγ+2Gcos[θ]formulae-sequenceΔsubscriptΔ0Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐2superscriptsubscript𝑐𝑠2𝑚superscript𝐿2superscriptsubscript𝜔𝑚2italic-ϕ2𝐺𝜃Δ𝛾2𝐺𝜃\Delta=\Delta_{0}-\frac{\hbar\left(\omega_{c}\right)^{2}\left|c_{s}\right|^{2}% }{mL^{2}\omega_{m}^{2}},\phi=\arctan{\frac{2G\sin[\theta]-\Delta}{\gamma+2G% \cos[\theta]}}roman_Δ = roman_Δ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - divide start_ARG roman_ℏ ( italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , italic_ϕ = roman_arctan divide start_ARG 2 italic_G roman_sin [ italic_θ ] - roman_Δ end_ARG start_ARG italic_γ + 2 italic_G roman_cos [ italic_θ ] end_ARG (11)

are the effective cavity detuning, and the phase of the cavity field. After linearisation, q^qs+q^^𝑞subscript𝑞𝑠^𝑞\hat{q}\rightarrow q_{s}+\hat{q}over^ start_ARG italic_q end_ARG → italic_q start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + over^ start_ARG italic_q end_ARG, p^ps+p^^𝑝subscript𝑝𝑠^𝑝\hat{p}\rightarrow p_{s}+\hat{p}over^ start_ARG italic_p end_ARG → italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + over^ start_ARG italic_p end_ARG, c^cs+c^^𝑐subscript𝑐𝑠^𝑐\hat{c}\rightarrow c_{s}+\hat{c}over^ start_ARG italic_c end_ARG → italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + over^ start_ARG italic_c end_ARG, and Fourier transformation of Eq. 9, we can solve the equations in frequency domain and derive the position fluctuations of the movable mirror as

q^(ω)=^𝑞𝜔absent\displaystyle\hat{q}(\omega)=over^ start_ARG italic_q end_ARG ( italic_ω ) = 1d(ω){(Δ2+(γiω)24G2)ξ^(ω)\displaystyle-\frac{1}{d(\omega)}\{(\Delta^{2}+(\gamma-i\omega)^{2}-4G^{2})% \hat{\xi}(\omega)- divide start_ARG 1 end_ARG start_ARG italic_d ( italic_ω ) end_ARG { ( roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ - italic_i italic_ω ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) over^ start_ARG italic_ξ end_ARG ( italic_ω ) (12)
\displaystyle-- i2γω0Lf[((ω+iγΔ)c^s+2iGeiθc^s)c^in(ω)+H.c.]},\displaystyle{i\hbar\sqrt{2\gamma}\frac{\omega_{0}}{L_{\rm{f}}}[((\omega+i% \gamma-\Delta)\hat{c}_{s}+2iGe^{i\theta}\hat{c}_{s}^{\dagger})\hat{c}_{in}^{% \dagger}(\omega)+\rm{H.c.}]\}},italic_i roman_ℏ square-root start_ARG 2 italic_γ end_ARG divide start_ARG italic_ω start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_L start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT end_ARG [ ( ( italic_ω + italic_i italic_γ - roman_Δ ) over^ start_ARG italic_c end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + 2 italic_i italic_G italic_e start_POSTSUPERSCRIPT italic_i italic_θ end_POSTSUPERSCRIPT over^ start_ARG italic_c end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ) over^ start_ARG italic_c end_ARG start_POSTSUBSCRIPT italic_i italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT ( italic_ω ) + roman_H . roman_c . ] } ,

with

d(ω)=𝑑𝜔absent\displaystyle d(\omega)=italic_d ( italic_ω ) = 2(ωcL)2|cs|2(Δ+2Gsin[θ2ϕ])2Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐𝐿2superscriptsubscript𝑐𝑠2Δ2𝐺𝜃2italic-ϕ\displaystyle 2\hbar(\frac{\omega_{c}}{L})^{2}\lvert c_{s}\rvert^{2}\left(% \Delta+2G\sin[\theta-2\phi]\right)2 roman_ℏ ( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( roman_Δ + 2 italic_G roman_sin [ italic_θ - 2 italic_ϕ ] ) (13)
+m(ω2ωm2+iωγm)[Δ2+(γiω)24G2].𝑚superscript𝜔2superscriptsubscript𝜔𝑚2𝑖𝜔subscript𝛾𝑚delimited-[]superscriptΔ2superscript𝛾𝑖𝜔24superscript𝐺2\displaystyle+m\left(\omega^{2}-\omega_{m}^{2}+i\omega\gamma_{m}\right)\left[% \Delta^{2}+(\gamma-i\omega)^{2}-4G^{2}\right].+ italic_m ( italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_i italic_ω italic_γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ) [ roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ - italic_i italic_ω ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] .

The susceptibility of the mechanical resonator is modified by the optomechanical interaction as

χxx1=d(ω)Δ2+(γiω)24G2=χ01+χos1superscriptsubscript𝜒𝑥𝑥1𝑑𝜔superscriptΔ2superscript𝛾𝑖𝜔24superscript𝐺2superscriptsubscript𝜒01superscriptsubscript𝜒os1\displaystyle\chi_{xx}^{-1}=-\frac{d(\omega)}{\Delta^{2}+(\gamma-i\omega)^{2}-% 4G^{2}}=\chi_{0}^{-1}+\chi_{\rm{os}}^{-1}italic_χ start_POSTSUBSCRIPT italic_x italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - divide start_ARG italic_d ( italic_ω ) end_ARG start_ARG roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ - italic_i italic_ω ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = italic_χ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT + italic_χ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT (14)
=m(ω2ωm2+iωγm)+2(ωcL)2|cs|2(Δ+2Gsin[θ2ϕ])Δ2+(γiω)24G2,absent𝑚superscript𝜔2superscriptsubscript𝜔m2𝑖𝜔subscript𝛾m2Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐𝐿2superscriptsubscript𝑐𝑠2Δ2𝐺𝜃2italic-ϕsuperscriptΔ2superscript𝛾𝑖𝜔24superscript𝐺2\displaystyle=-m(\omega^{2}-\omega_{\rm{m}}^{2}+i\omega\gamma_{\rm{m}})+\frac{% 2\hbar(\frac{\omega_{c}}{L})^{2}\lvert c_{s}\rvert^{2}\left(\Delta+2G\sin[% \theta-2\phi]\right)}{\Delta^{2}+(\gamma-i\omega)^{2}-4G^{2}},= - italic_m ( italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_i italic_ω italic_γ start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT ) + divide start_ARG 2 roman_ℏ ( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( roman_Δ + 2 italic_G roman_sin [ italic_θ - 2 italic_ϕ ] ) end_ARG start_ARG roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_γ - italic_i italic_ω ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,

where χ0=1m(ω2ωm2+iωγm)subscript𝜒01𝑚superscript𝜔2superscriptsubscript𝜔m2𝑖𝜔subscript𝛾m\chi_{0}=-\frac{1}{m(\omega^{2}-\omega_{\rm{m}}^{2}+i\omega\gamma_{\rm{m}})}italic_χ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - divide start_ARG 1 end_ARG start_ARG italic_m ( italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_i italic_ω italic_γ start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT ) end_ARG is the intrinsic mechanical susceptibility and the second term χos1superscriptsubscript𝜒os1\chi_{\rm{os}}^{-1}italic_χ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT represents the optomechanical modification to the susceptibility with the OPA. To make it clear, we rewrite this equation as

χxx1=m(ω2+i[γmΓos]ω[ωm2+Kos/m]),superscriptsubscript𝜒𝑥𝑥1𝑚superscript𝜔2𝑖delimited-[]subscript𝛾msubscriptΓ𝑜𝑠𝜔delimited-[]superscriptsubscript𝜔m2subscript𝐾os𝑚\chi_{xx}^{-1}=-m(\omega^{2}+i[\gamma_{\rm{m}}-\Gamma_{os}]\omega-[\omega_{\rm% {m}}^{2}+K_{\rm{os}}/m]),italic_χ start_POSTSUBSCRIPT italic_x italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT = - italic_m ( italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_i [ italic_γ start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT - roman_Γ start_POSTSUBSCRIPT italic_o italic_s end_POSTSUBSCRIPT ] italic_ω - [ italic_ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT / italic_m ] ) , (15)

where the optical spring constant Kossubscript𝐾osK_{\rm{os}}italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT and optical damping ΓossubscriptΓos\Gamma_{\rm{os}}roman_Γ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT are defined as:

Kossubscript𝐾os\displaystyle K_{\rm{os}}italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT =Re[χos1]absent𝑅𝑒delimited-[]superscriptsubscript𝜒os1\displaystyle=Re[\chi_{\rm{os}}^{-1}]= italic_R italic_e [ italic_χ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ] (16)
=2(ωcL)2|cs|2(γ2+Δ24G2ω2)(Δ+2Gsin[θ2ϕ])(γ2+Δ24G2ω2)2+4γ2ω2,absent2Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐𝐿2superscriptsubscript𝑐𝑠2superscript𝛾2superscriptΔ24superscript𝐺2superscript𝜔2Δ2𝐺𝜃2italic-ϕsuperscriptsuperscript𝛾2superscriptΔ24superscript𝐺2superscript𝜔224superscript𝛾2superscript𝜔2\displaystyle=\frac{2\hbar(\frac{\omega_{c}}{L})^{2}\lvert c_{s}\rvert^{2}(% \gamma^{2}+\Delta^{2}-4G^{2}-\omega^{2})(\Delta+2G\sin[\theta-2\phi])}{(\gamma% ^{2}+\Delta^{2}-4G^{2}-\omega^{2})^{2}+4\gamma^{2}\omega^{2}},= divide start_ARG 2 roman_ℏ ( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ( roman_Δ + 2 italic_G roman_sin [ italic_θ - 2 italic_ϕ ] ) end_ARG start_ARG ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,
ΓossubscriptΓos\displaystyle\Gamma_{\rm{os}}roman_Γ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT =Im[χos1]meffω=4(ωcL)2|cs|2γ(Δ+2Gsin[θ2ϕ])meff((γ2+Δ24G2ω2)2+4γ2ω2)absent𝐼𝑚delimited-[]superscriptsubscript𝜒os1subscript𝑚eff𝜔4Planck-constant-over-2-pisuperscriptsubscript𝜔𝑐𝐿2superscriptsubscript𝑐𝑠2𝛾Δ2𝐺𝜃2italic-ϕsubscript𝑚𝑒𝑓𝑓superscriptsuperscript𝛾2superscriptΔ24superscript𝐺2superscript𝜔224superscript𝛾2superscript𝜔2\displaystyle=\frac{Im[\chi_{\rm{os}}^{-1}]}{m_{\rm{eff}}\omega}=\frac{4\hbar(% \frac{\omega_{c}}{L})^{2}\lvert c_{s}\rvert^{2}\ \gamma(\Delta+2G\sin[\theta-2% \phi])}{m_{eff}((\gamma^{2}+\Delta^{2}-4G^{2}-\omega^{2})^{2}+4\gamma^{2}% \omega^{2})}= divide start_ARG italic_I italic_m [ italic_χ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ] end_ARG start_ARG italic_m start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT italic_ω end_ARG = divide start_ARG 4 roman_ℏ ( divide start_ARG italic_ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_L end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_γ ( roman_Δ + 2 italic_G roman_sin [ italic_θ - 2 italic_ϕ ] ) end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT ( ( italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_Δ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 4 italic_G start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 4 italic_γ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) end_ARG

The new effective mechanical frequency and damping rate are ωeff=ωm2+Kos/msubscript𝜔effsuperscriptsubscript𝜔m2subscript𝐾os𝑚\omega_{\rm{eff}}=\sqrt{\omega_{\rm{m}}^{2}+K_{\rm{os}}/m}italic_ω start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT = square-root start_ARG italic_ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_K start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT / italic_m end_ARG and γeff=γm+Γossubscript𝛾effsubscript𝛾msubscriptΓos\gamma_{\rm{eff}}=\gamma_{\rm{m}}+\Gamma_{\rm{os}}italic_γ start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT = italic_γ start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT + roman_Γ start_POSTSUBSCRIPT roman_os end_POSTSUBSCRIPT.

5.2 Mechanical resonator

The mechanical resonator used in the experiment is shown in Fig.6. It is a through-etched silicon frame with a multi-layer AlGaAs/GaAs coating bonded at the backside. There are in total 25 square shape windows with the biggest being 1.16 ×\times× 1.16 mm and the smallest being 0.175 ×\times× 0.175 mm. Each window is an independent resonator with different resonance frequencies and quality factors. The one presented in this work is of size 0.64 ×\times× 0.64 mm.

Refer to caption
Figure 6: The mechanical resonator. Left penal: front view of the resonator, showing 25 different size windows. right panel: photo showing the laser beam on the 0.64 ×\times× 0.64 mm window, the beam size is about 150 μ𝜇\muitalic_μm.

The target mechanical mode frequency is about 85 kHz, with a quality factor of around 2500. Multiple modes associated with the whole frame or the frame holder have much lower quality factors, shown in Fig.7.

Refer to caption
Figure 7: Cavity locking error signal PSD, showing the mechanical mode at a frequency around 85 kHz.

5.3 Phase locking loop error signal

The phase-locking loop error signal is obtained by modulating the green laser beam phase through an EOM at 300 kHz while monitoring the intra-cavity power in the seed beam polarization. The intra-cavity power would also be modulated because the OPA amplification is different for different OPA phases. The intra-cavity field is then transmitted through the cavity mirror and detected by the PD outside the vacuum chamber. The PD signal is then mixed with the 300 kHz local oscillator and low-passed to produce the phase-locking error signal. Fig.8 top panel shows the time traces of the transmitted power in blue and the corresponding error signal in orange. When the OPA is locked at the error signal equal to zero, the OPA is either at maximum amplification or maximum de-amplification status. The OPA can also be locked at an arbitrary phase by offsetting the error point of the phase-locking loop. The error signal is a sinusoidal function of the OPA phase, thus by comparing the error signal offset (e.g.formulae-sequence𝑒𝑔e.g.italic_e . italic_g . 40 mV shown in Fig.9) and the full-size error signal (147.5 mV ±plus-or-minus\pm± 7.5 mV), the OPA phase can be calculated (θerrsig=40mV=15.70.8+0.9subscript𝜃errsig40mVsubscriptsuperscriptsuperscript15.7superscript0.9superscript0.8\theta_{\rm{errsig=40mV}}={15.7^{\circ}}^{+0.9^{\circ}}_{-0.8^{\circ}}italic_θ start_POSTSUBSCRIPT roman_errsig = 40 roman_m roman_V end_POSTSUBSCRIPT = 15.7 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT + 0.9 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.8 start_POSTSUPERSCRIPT ∘ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ).

Refer to caption
Figure 8: Time traces of OPA phase scanning, the blue trace in the upper panel represents the power of the laser beam at the seed beam polarisation; the orange trace is the phase-locking loop error signal. The orange horizontal dashed line is the zero level of the phase-locking error signal. The vertical dashed lines indicate that the OPA is either at the maximum amplification (purple) or maximum de-amplification (green). The lower panel is the 2 Hz triangle scanning voltage applied on the PZT mirror where the phase locking loop error signal will be applied to lock the OPA phase.
Refer to caption
Figure 9: Optical spring effect by tuning OPA phase. The left panel from top to bottom shows the power of the beam at the seed beam polarisation, the power of the locking beam, the phase locking loop error signal, and the OPA phase locking loop feedback control signal. The blue traces are when the phase locking loop locked at the error signal equals -40 mV and 40 mV for the orange trace. The right side panel is the corresponding mechanical mode power spectrum measured from the OPA cavity length locking loop error signal, showing that the OS effect can be controlled by adjusting the OPA phase without cavity detuning.