1 Introduction
Cavity optomechanics has sparked extensive theoretical and experimental research interest in the last decade due to its various applications [
1–
12]: detecting weak-force, mass, displacements, and orbital angular momentum [
13–
16]; creating macroscopic nonclassical states [
17,
18]; achieving mechanical squeezing [
19,
20]; obtaining optical nonreciprocity [
21]. These advancements exhibit potential advantages of cavity optomechanics in quantum metrology [
22], quantum information processing [
23–
35], and fundamental physics questions [
36–
40].
Though optomechanical systems have brought wide attention, the single-photon optomechanical coupling
g0 is small in experiment [
41,
42]. People usually bring in classical driving lasers to improve the effective optomechanical coupling strength [
43–
47]. However, the phase and amplitude of the driving laser have small fluctuation — the so-called laser phase and amplitude noise. Schliesser et al. firstly observed the laser phase noise in experiment [
48]. Generally, laser amplitude noise can be approximately neglected by stabilizing the laser intensity [
49,
50]. Recently, many efforts have been devoted to study the influence of laser phase noise on quantum physics: ground-state cooling [
51–
56]; quantum state estimation [
57]; weak-force sensing [
58,
59]; squeezing mechanical oscillators [
60] and output light of an optical cavity [
61]; quantum memory [
62]; optomechanical entanglement [
63–
65]. Many of these studies demonstrate that phase noise has destructive effects on quantum phenomena. Optomechanical entanglement is a macroscopic quantum phenomenon and is significant for quantum information processing [
66,
67]. Some works have discussed the influence of the Kerr nonlinear medium on the stationary optomechanical entanglement in the presence of laser phase noise [
65,
68,
69]. Although they claim that Kerr nonlinearity can promote entanglement to some extent, they leave an unresolved contradiction: strong optical force coupling accompanied by huge laser phase noise. Specifically, the effect of laser phase noise is described by
where
α and
are the mean value of intracavity field and time derivatives of phase noise, respectively. The effective optomechanical coupling is
where the single-photon coupling
g0 is very weak in experiment. Therefore, people usually improve the effective optomechanical coupling
by raising
α. However, the increase of
α is unavoidable to enlarge the effect of phase noise. Naturally, we raise a novel and interesting idea: Can we simultaneously strengthen the effective optomechanical coupling and inhibit the effect of laser phase noise?
In this paper, we present a scheme to improve the effective optomechanical coupling and inhibit laser phase noise at the same time. The system consists of an optical cavity, a mechanical oscillator, and a Kerr nonlinear medium. We obtain the effective Hamiltonian with optical and mechanical parametric amplification terms after linearizing the systemic Hamiltonian and then transform the system into the squeezing frame. According to the analytical results, we found that the effective optomechanical coupling is enhanced, and the laser phase noise decreases exponentially by adjusting the squeezing parameters. Generally, we can introduce a broadband-squeezed vacuum environment to suppress the increased thermal noise around the squeezed cavity field and mechanical oscillator. However, the enhanced thermal noise around the mechanical oscillator has little influence on the system due to its tiny decay, which means we only need to suppress the enlarged thermal noise around the squeezed cavity field by exploiting the vacuum environment. Our calculation shows that our proposal effectively solves the contradiction: improving effective optomechanical coupling leads to enlarging laser phase noise. To display our scheme, we exploit it to simulate quantum memory and stationary optomechanical entanglement as examples. Our numerical results show our strategy can significantly protect both the quantum memory and the stationary optomechanical entanglement.
This paper is organized as follows: in Section 2, We explain the theoretical model, derive the effective Hamiltonian, and obtain the dynamic equation of the additional mode describing the effect of phase noise. We demonstrate some actual physical phenomena (quantum memory and stationary optomechanical entanglement) in Section 3. Finally, we give a conclusion in Section 4.
2 Physical model
We consider an optomechanical system and its schematic diagram is illustrated in Fig. 1. In this setup, a classical laser drives the optical cavity, and a Kerr medium is located in the cavity. Kumar
et al. [
70] have shown that Kerr medium inside an optomechanical system can effectively inhibit the normal mode splitting. In the rotating frame with frequency
ωL, the Hamiltonian of the system can be written as (
) [
71,
72]
where
and
are the annihilation (creation) operators of the cavity field and the mechanical oscillator, respectively. The optical cavity mode with resonance frequency
ωc contains a Kerr nonlinear medium with Kerr coefficient
u. The detuning ∆
0 satisfies ∆
0 =
ωc−
ωL. The mechanical oscillator with resonance frequency
ωm is accompanied by a Duffing nonlinear term with amplitude
η and coupled to the cavity field with coupling strength
g0. By coupling the mechanical mode to a qubit, a strong Duffing nonlinearity can be obtained, in which the nonlinear amplitude
η can reach 10
−4ωm [
73].
ωL is the frequency of the classical driving laser,
φ(
t) describes the laser phase noise — a zero-mean stationary Gaussian stochastic process, and
EL is the strength of the driving laser. Generally, one can adjust the driving strength
by controlling the input power
Ps where
κ is the decay of the cavity through its input port. The amplitude noise of the driving laser is negligible compared to the phase noise via stabilizing laser source, so we can neglect the amplitude noise of the driving laser in this paper [
62,
63].
The optical Kerr effect has been studied theoretically and experimentally, and the corresponding nonlinear Kerr coefficient has the following form [
74,
75]:
with
where
c represents the speed of light in vacuum,
n0 (
n2) are the linear (nonlinear) refractive index of the material with ranges 2⩽
n0 ⩽4 and 10
−13 ≥
n2 ≥10
−17 cm
2/W [
76].
Veff is defined as the effective mode volume and describes the peak electric field strength within the cavity.
ε(
r) and Φ(
r) are the dielectric constant and the electric field strength, respectively.
Veff is located in between 10
2–10
4 µm
3 when the quality factor of the cavity
Q is limited in (10
6–10
8) [
77]. For a near-infrared wavelength (
λ = 1064 nm), the nonlinear Kerr coefficient
u is estimated on the order between 0.0006–2721 Hz in a silica microsphere by calculating Eq. (3) with experimentally accessible parameters.
We consider the full description of the systemic dynamics including the fluctuation-dissipation processes of the optical and the mechanical modes. After transforming the cavity mode to a randomly rotating frame according to , we derive a set of quantum Langevin equations governing the systemic dynamics
where
γm describes the damping rate of the mechanical mode,
is the thermal noise operator acting on the mechanical oscillator, and
is the squeezed input vacuum noise operator. The squeezed input optical field has been used to enhance sideband cooling and suppress Stokes scattering [
78,
79]. The corresponding noise correlations are written as
where N = sinh2 (re) is the mean photon number of the broadband-squeezed vacuum environment, M = is the strength of the autocorrelation of the squeezed vacuum noise, and Nth = is the equilibrium mean thermal photon number. Here re (Φe) are the squeezing amplitude (angle) of the broadband-squeezed vacuum environment, kB is the Boltzmann constant, and T is the temperature of the mechanical oscillator. The system is linearized by substituting operators and to Eq. (4) where α (β) and describe the mean values and the fluctuations of the optical (mechanical) mode. Therefore, we can simplify the linearized Langevin equation as
where we have ignored the higher-order nonlinear terms, and the effective parameters are listed as: ∆=∆0− ; Ω = 2uα2; ; Ωm = . The mean values of the optical and the mechanical modes can be obtained by solving the steady Langevin equation:
We rewritten the strength of parametric amplification coefficient as Ω = |Ω|e−2iθ (i.e., α = |α|e−iθ) with real angle θ. The phase factors can be absorbed into the operators . Therefore, we can obtain the linearized Hamiltonian
where we have defined the optomechanical coupling g as g = g0|α|. It is obvious that the Kerr nonlinear medium and the Duffing nonlinearity lead to the optical and the mechanical parametric amplification terms with amplitudes |Ω| and Ωm, respectively. We exploit squeezing transformations and acting on the linearized Hamiltonian where the squeezing phase is fixed as π. Here we choose the squeezing strength and with η = |Ω|/∆ and . Therefore, we can obtain the following effective Hamiltonian
where the effective coupling strength is G = g exp (r′) with r′ = rm−r, and the effective detuning of the optical and mechanical modes are and . We adjust the squeezing amplitudes to satisfy rm≥r, and thus the effective G can remain a large value.
Then we discuss the statistical properties of laser phase noise. As shown in Eq. (6), we note that the laser phase noise affects the systemic dynamics by the additional noise term
and the influence of phase noise is mainly depended on the mean value of the cavity field
α. Generally, the single-photon coupling
g0 is very small and one need to improve the effective optomechanical coupling via a large
α. It is a terrible contradiction that the large mean value of cavity field
α will lead to a sizeable effective optomechanical coupling and phase noise when
α is huge. Therefore, it is significant to enhance the effective optomechanical coupling and suppress the influence of phase noise at the same time. If the phase noise correlation satisfies
, the spectrum of the noise is flat and the cut-off frequency
γc→∞(i.e.,
, where Γ
L is the linewidth of the driving laser. However, the spectral density of the phase noise is not a flat spectrum due to the finite non-zero correlation time of phase noise. In other words, it is a finite bandwidth color noise. Generally, the noise spectrum is equivalent to a low pass filtered white noise with the following spectrum and correlation function [
54,
58,
62]
where 1/γc is correlation time of the laser phase noise so that the phase noise is suppressed at frequencies ω > γc. The correlation time decreases and the frequency noise starts reaching the white noise with the increasing of γc. Moreover, the frequency spectrum in Eq. (10) is equivalent to the differential equation where ε(t) is a Gaussian random variable with the noise correlation function
We redefine an additional noise operator where ψ satisfies the following differential equation
Therefore, we can rewrite the Langevin equation with the quadrature fluctuations of the optical field and the mechanical oscillator: ; ; ; . Moreover, the corresponding input noise operators are amended as , , and with . We note that noise and own exponential factors, which means decreasing laser phase noise is accomplished by increasing thermal noise. To suppress the increased thermal noise around the cavity field, we adjust the amplitude and phase of the squeezed vacuum environment. If the squeezing parameters satisfy the conditions r = re and Φe−2θ = π, the effective input noise of the cavity is equivalent to a vacuum noise and we can obtain the following noise correlation function
We derive these noise correlations in detail in Appendix A. By combining Eqs. (9), (12) and (13), we can derive the Langevin equation to describe the dynamic evolution of the system. We rewrite the Langevin equation as a compact matrix form
where we have defined the vector of continuous variable fluctuation operators , and the corresponding input noise vector is n = . Moreover, the drift matrix is the 5×5 matrix
where the element in the drift matrix A describes the coupling between the phase noise operator and the optical momentum operator. It is different from the standard cavity optomechanical system that the effective coupling and the phase noise term multiply exponential factors and e−r, respectively. One can enlarge the squeezing strength r to suppress the influence of phase noise and increase r′ to improve the coupling G. Therefore, the parametric processes, induced by the Kerr medium and Duffing nonlinearity, can simultaneously increase the effective optomechanical coupling and reduce the coupling between the laser phase noise and the cavity field. According to the Eq. (14), we obtain the following dynamical equation of the covariance matrix
where the matrix element of the covariance matrix V can be expressed as
and the corresponding noise matrix is
where we have defined the parameter . Generally, one exploits a squeezed vacuum bath to counteract the influence of the factor . However, the mechanical decay γm is very small so that the factor in Eq. (18) has a little affect on the systemic dynamics. Therefore, we retain this factor in the following calculations.
3 Demonstrating some actual quantum phenomenon
In this section, we take two examples to test the efficiency of our scheme mentioned in the above section. First, we theoretically investigate the performance of our proposal on improving the optomechanical quantum memories against the laser phase noise. Second, we exploit our design to demonstrate the stationary optomechanical entanglement.
3.1 Quantum memory
Quantum memory is indispensable for quantum information processing and has made enormous progress in optics and atoms [
80]. Let us briefly recall the quantum memory in optomechanical system [
81]. As shown in Fig. 2, the state of the optical mode is transferred to the mechanical mode in the writing process with time
. Then the state is stored in the mechanical membrane for a while
τ by decoupling the mechanical and the optical modes. Finally, the optical mode obtains the storied state in the reading process, and the corresponding reading time is
. One of the advantages of quantum memory in optomechanical systems is that the decay rate of the mechanical oscillator is much smaller than the optical cavity.
Here we apply the model to achieve quantum memory. We suppose the effective detuning satisfies the resonate condition ∆e≈∆m≫G which can be achieved by controlling the frequency of the optical driving. This condition limits that η cannot close to one infinitely. In other words, we cannot fully cancel the laser phase noise. In our paper, we limit η in 0–0.9999. is defined as the fidelity between the initial state and final state . In the phase space, one can rewrite the fidelity as , where is the vector of the optical quadratures and Wi (Wf) are the Wigner functions of the optical initial (final) states. To simplify the calculation, we assume the state of the cavity staying in a pure Gaussian state — the Wigner function Wi (Wf) are the Gaussian distribution. Under this assumption, one can obtain the fidelity between the initial and the final states of the optical mode at time
where the parameters , Θ has the following forms
with the initial (final) covariance matrix
Vi (
Vf) and the initial (final) optical mean quadratures
ξi (
ξf). The detailed derivation of Eq. (19) have been proposed by Wang [
82].
To measure the fidelity F, we need to calculate the expectations of optical quadratures (ξi and ξf) and the covariance matrix (Vi and Vf) of the initial and final states. Here we consider the squeezed coherent state (a typical Gaussian state) as the initial state of the optical field in the squeezing frame (i.e., |µ, χ⟩ = D(µ)S(χ)|0⟩), and thus the state in the original frame is S(r)†D(µ)S(χ)|0⟩ where and are displacement and squeezing operators, respectively. Therefore, one can obtain the initial vector u(0) = and the corresponding initial covariance matrix is
where we have assumed .
To numerically analyze the quantum memory effect, we choose the parameters similar to those in Refs. [
63,
64]: length of the cavity
L = 1 mm; wavelength of the cavity field 1064 nm; mass of the mechanical oscillator
m≃10 ng; frequency of the mechanical oscillator
ωm = 2π ×10 MHz; quality factor
Qm = 2 × 10
6; the optical decay rate
κ = 2π × 100 kHz; single-photon coupling strength
g0 = 2π × 100 Hz; the linewidth of driving laser Γ
L in range 1∼50 kHz; the cut-off frequency
γc in range 0.1–50 kHz. Moreover, we assume the mean thermal phonon number
nth = 3 and the parameter
r′ = 0. The storage time is
τ =
= 1.035 µs. We fix the effective coupling
G = 0.05
ωm which can be achieved by controlling the strength of driving laser.
To demonstrate the advantage of our proposal, we simulate the fidelity F as a function of parameter η in Figs. 3(a) and (b) for different cut-off frequency γc and laser linewidth ΓL, respectively, where we have limited the parameter η in interval [0, 0.9999]. It is clear that the fidelity is improving with the increasing of parameter η even the noise spectrum has large cut-off frequency γc and laser linewidth ΓL. In particular, the fidelity is approximately the same with ideal situation (i.e., without laser phase noise γc = 0) for η = 0.9999 though the cutoff frequency γc takes the values 5 kHz, 8 kHz, 10 kHz, and 15 kHz. It indicates that the phase noise is extremely suppressed even can be approximately ignored, which is largely different from the case η = 0 (the standard optomechanical system). We improve the parameter η to increase the squeezing parameter r (i.e., enhancing er) so that the noise term can be effectively suppressed. Moreover, we numerically simulate the variation of |α| with η in Fig. 4. It is obvious that |α| is monotonically decreasing with the increasing of parameter η and the minimum value |α|e−r is 420.45 for η = 0.9999. At this time, the effect of phase noise is reduced about an order of magnitude.
In Fig. 5, we simulate the fidelity F versus to the increasing of the squeezing amplitude of the initial state which is described by the squeezing parameters χ. One can easily find the fidelity decreases by increasing parameter χ, while the variation of the fidelity for η = 0.9 is slower than η = 0.7, η = 0.5, and η = 0. Therefore, our scheme can protect the fidelity and inhibit the phase noise for a small χ. However, the advantage of the scheme slowly disappears and the initial state would be more sensitive to various noises when the initial state becomes more and more non-classical (i.e., with the increasing of χ).To further clarify the promotion effect of large η (i.e., large r), we simulate the fidelity as the function of cut-off frequency γc and laser linewidth ΓL for different η in Figs. 6(a)–(c). The results show that the area of high fidelity shrinks with the decreasing of η, and the destructive effect of the cut-off frequency γc on the fidelity F is more larger than the laser linewidth ΓL. Moreover, the fidelity can arrive at 0.955 for η = 0.9999 even the parameters of the noise spectrum are very huge (γc = 0.01 MHz and ΓL = 0.01 MHz).
3.2 Entanglement
Here we consider the stationary optomechanical entanglement as the second example to demonstrate the advantage of our scheme on suppressing the phase noise. The steady-state of the system is associated with Eq. (16). If and only if all the eigenvalues of the matrix
A have negative real parts, the system can arrive in the steady-state. We derive the stability conditions by using the Routh–Hurwitz criteria [
83–
86]. According to the criteria, we get the following two non-trivial stability conditions:
In the next calculation, we restrict all the parameters to satisfy the stable condition [Eq. (22a)] and [Eq. (22b)]. Here we consider the frequency condition ∆e,∆m > 0. It should be noticed that the critical condition [Eq. (22a)] limits the exponential improvement of coupling and noise suppression (i.e, 0 < r < rmax) where
For our model, the steady-state is a zero-mean Gaussian state because we have linearized the dynamics of the fluctuations, and all noises are Gaussian; as a consequence, it is fully characterized by the 5 × 5 stationary covariance matrix V (∞) with matrix elements,
The stationary covariance matrix V (∞) can be obtained by solving the following Lyapunov equation
We find that the Lyapunov equation (25) is linear for the covariance matrix V, which means the Lyapunov equation can be straightforwardly and analytically solved. One can identify all the quantum properties of the stationary state of the optomechanical system according to the stationary covariance matrix V (∞). Therefore, we can simulate the influence of laser phase noise on achieving quantum entanglement between the mechanical oscillator and optical mode. The stationary optomechanical entanglement relates to the mechanical and optical quadratures, and thus we concentrate on the reduced 4 × 4 covariance matrix of V (∞). This reduced correlation matrix has the following form
where
VA,
VB and
VC are 2×2 matrix. The matrix
VA (
VB) are associated with the optical mode (mechanical oscillator), while
VC describes the optomechanical correlations. The logarithmic negativity is the famous and convenient measure for predicting continuous variable (CV) entanglement [
87], and its definition is
where η− is the symplectic eigenvalue of the bipartite system, and it has the following form
with Σ(V) = det VA + det VB−2 det VC.
Then we study the advantage of our scheme on CV entanglement when the phase noise exists in the system. In Figs. 7(a) and (b), we simulate
EN versus to the optical decay
κ under different
η. Obviously, laser phase noise has a prominent effect on the stationary optomechanical entanglement while the large parameter
η can improve the maximum value of
EN and broaden the parameter region existing entanglement. By comparing with Figs. 7(a) and (b), we find our proposal (
η = 0.9999) has a great advantage than the standard optomechanical coupling model (
η = 0). Therefore, our scheme extremely inhibits the negative effect of laser phase noise on the stationary optomechanical entanglement, which improve the conclusion in Ref. [
63]. Moreover, we simulate
EN as a function of
κ for
η = 0.9999 and
η = 0 in Figs. 8(a) and (b), respectively. According to the numerical results, we summarize the advantages of our scheme versus the standard optomechanical coupling model: the entanglement
EN is more greater for huge
γc and Γ
L;
EN can exist in a wide range of
κ; the entanglement
EN decays more slowly with the increasing of
γc and Γ
L;
EN exists even for larger laser phase noise (
γc, Γ
L) = (20 kHz, 20 kHz). However, the standard optomechanical coupling model (
η = 0) is sensitive to laser phase noise, and the stationary optomechanical entanglement
EN is approximately close to zero for (
γc, Γ
L) = (15 kHz, 15 kHz). In Fig. 9, we simulate
EN versus the cut-off frequency
γc and the laser linewidth Γ
L for (a)
η = 0.9999, (b)
η = 0.999, and (c)
η = 99. The numerical results show that the destructive effect of phase noise is tiny when
η is closer to one. Although the maximum achievable entanglement decreases with the increasing of cut-off frequency
γc and laser linewidth Γ
L, we still obtain a large parameters range to maintain the entanglement for a large
η. Moreover, we also notice that the destructive effect of
γc is more remarkable than Γ
L. Therefore, our method can suppress the laser phase noise, and thus protect the stationary entanglement.
4 Conclusion
In summary, we studied a theoretical proposal to suppress the phase noise and improve the effective optomechanical coupling. The optomechanical system includes a mechanical Duffing nonlinearity and a Kerr medium that can create the mechanical and optical parametric amplification terms. Further calculation shows that we can enhance the effective optomechanical coupling and inhibit the laser phase noise at the same time. In this process, we use the squeezed vacuum environment to inhibit the increased thermal noise. To test the performance of our proposal, we simulate quantum memory and stationary optomechanical entanglement as examples. The numerical results show that our scheme effectively suppresses the destructive influence of laser phase noise on quantum memory and protects the storing fidelity at a high value. Moreover, our proposal can also protect stationary optomechanical entanglement. In particular, the maximal entanglement decreases very slowly with the increasing of the laser phase noise, and it exists in wide ranges of parameters. Our scheme provides a promising way for inhibiting the phase noise of optomechanical systems or other quantum systems driven by lasers and has potential applications for achieving quantum information processes and observing quantum phenomena.