Improving the Stationary Entanglement of a Laguerre–Gaussian Cavity Mode with a Rotating Mirror via Nonlinear Cross-Kerr Interactions and Parametric Interactions

Quantum entanglement is essential in performing many quantum information tasks. Here, we theoretically investigate the stationary entanglement between a Laguerre–Gaussian (LG) cavity field and a rotating end mirror in an LG-cavity optorotational system with a nonlinear cross-Kerr (CK) interaction and a degenerate optical parametric amplifier (OPA). We calculate the logarithmic negativity of the system to quantify the stationary entanglement. We examine the influence of various system parameters such as the cavity detuning, the strength of the nonlinear CK interaction, the parametric gain and phase of the OPA, the power of the input Gaussian laser, the topological charge of the LG-cavity field, the mass of the rotating end mirror, and the ambient temperature on the stationary entanglement. Under the combined effect of the nonlinear CK interaction and the OPA, we find that the stationary entanglement can be substantially enhanced at lower Gaussian laser powers, smaller topological charges of the LG-cavity field, and larger masses of the rotating end mirror. We show that the combination of the nonlinear CK interaction and the OPA can make the stationary entanglement more robust against the ambient temperature.


Introduction
The experimental achievement of the laser cooling of mechanical modes [1,2] to the quantum ground state has laid a solid foundation for further realizing non-classical quantum effects and their applications.Over the past few decades, significant progress has been made in various aspects of cavity optomechanics [3].Optomechanical systems have potential applications in optical information processing [3].They have been proven to function as all-optical memory elements [4] and have been proposed as a new technology for single-photon detection [5].Moreover, the effective preparation of various quantum states in optomechanical systems has been a subject of extensive research [3].Entangled states, in particular, are valuable quantum resources [6].Notably, optomechanical systems provide a robust platform where scientists can utilize optomechanical interactions to create entanglement between the cavity field and the macroscopic mechanical motion [7]; this enables various studies of quantum entanglement in macroscopic systems [8], thereby bridging the quantum-classical boundary [9].The essence of these systems lies in the linear momentum exchange between photons and mechanical oscillators [7].Additionally, the exchange of photon angular momentum can be utilized to design a Laguerre-Gaussian (LG) cavity optorotational system [10][11][12][13][14], in which a LG-cavityty mode is coupled to a rotating end mirror.The LG-cavityty mode has a well-defined orbital angular momentum equal to lh per photon, where l is the azimuthal mode index or the topological charge [15].It has been demonstrated that the rotation of macroscopic mirrors can be cooled to near their quantum ground state [10].Based on this, more quantum effects have been explored, such as the entanglement between LG-cavity fields and rotating mirrors [11,12] and rotational mirror-mirror entanglement [13,14].
Researchers have continuously striven to achieve stronger optomechanical coupling.One innovative approach [16] involves introducing a quantum two-level system (qubit) into optomechanical systems.Experiments [17] have demonstrated that using a quantum two-level system (qubit or artificial atom) can facilitate nonlinear cross-Kerr (CK) coupling between a mechanical resonator and a microwave cavity mode, enhancing the radiationpressure coupling by six orders of magnitude.This enhancement allows the system to approach a strong coupling regime.Studies have shown that nonlinear CK interaction can improve the entanglement between a cavity field and a moving mirror [18,19].In addition, our previous work has shown that nonlinear CK interaction can enhance the stationary entanglement between an LG-cavity mode and a rotating end mirror [20].The other method is placing a degenerate optical parametric amplifier (OPA) into an optomechanical cavity to increase the optomechanical coupling, which is based on the fact that the squeezed light can be generated from a degenerate OPA [21].It has been reported that a degenerate OPA in an optomechanical system can generate the multipartite optomechanical entanglement among the fundamental optical mode, the second harmonic optical mode, and the mechanical mode [22] and can enhance the stationary entanglement between two charged mechanical oscillators coupled through a Coulomb interaction [23].
In this paper, we investigate the stationary entanglement between an LG-cavity field and a rotating end mirror in an LG-cavity optorotational system with a nonlinear CK interaction and a degenerate OPA.Under the stability conditions of the system, we discuss the influence of the cavity detuning, the strength of the nonlinear CK interaction, the parametric gain and phase of the OPA, the power of the incident Gaussian laser, the topological charge of the LG-cavity field, the mass of the rotating end mirror, and the ambient temperature on the stationary entanglement.Compared to the results in our previous work [20], we find that the addition of the OPA into the optical cavity can substantially enhance the stationary entanglement at lower input powers, smaller topological charges of the LG-cavity field, and larger masses of the rotating end mirror and can make the entanglement more robust against the environmental noise.The entanglement measure, the approximations used, and the computational methods were first introduced by Bhattacharya et al. [11] and used in our previous work [20].
We have organized this paper as follows.In Section 2, we start with an introduction of the theoretical model, give the Hamiltonian of the whole system, derive the Langevin equations of the system, and obtain steady-state solutions.In Section 3, we obtain the linearized quantum Langevin equations and use logarithmic negativity to quantify the stationary entanglement between the LG-cavity field and the rotating end mirror.In Section 4, we discuss the effects of the cavity detuning, the strength of the nonlinear CK interaction, the gain and phase of the OPA, the incident Gaussian laser power, the topological charge of the LG-cavity field, the mass of the rotating end mirror, and the ambient temperature on the stationary optorotional entanglement.A concluding summary is given in Section 5.

Model
The system we consider is an LG-cavity optorotational system formed by two cavity end mirrors with separation L [11], as depicted in Figure 1.The left fixed mirror transmits a small fraction of the incoming beam, and the right rotating mirror fully reflects the incoming beam.Both mirrors are spiral-phase plates [10].The fixed mirror does not modify the topological charge of any laser beam transmitting through it and reduces a topological charge 2l from the laser beam during reflection.The rotating mirror adds a topological charge 2l to the laser beam during reflection.Thus, when a Gaussian laser beam is sent into the cavity, an LG-cavity field is built up inside the cavity.The LG-cavity field provides a torque 2lh 2L/c = hg ϕ to the rotating end mirror, where g ϕ = cl L is the optorotational coupling strength between the LG-cavity field and the rotating end mirror and c is the light speed in vacuum [10].The rotating end mirror with mass m, frequency ω m , and damping rate γ m can be approximated as a harmonic oscillator.It is assumed that a quantum twolevel system is introduced on the rotating end mirror, which generates the nonlinear CK interaction between the LG-cavity field and the rotating end mirror [17].The strength of the nonlinear CK interaction is denoted by g ck .The angular displacement operator ϕ and angular momentum operator L z of the rotating end mirror can be, respectively, written as ϕ = 1 It is assumed that the frequency ω m of the rotating end mirror is much lower than the cavity resonance frequency ω c .Thus, the dynamical Casimir effect can be ignored, which describes the generation of photons from the quantum vacuum due to the rotation of the end mirror [24].Meanwhile, the rotational Doppler effect due to the angular momentum exchange during the optorotational interaction can also be neglected [25].LG-cavity optorotational system with a nonlinear CK interaction and a degenerate OPA.There is a two-level system (red) on the rotating end mirror to produce the nonlinear CK interaction between the LG-cavity field and the rotating end mirror.And a nonlinear optical crystal with a second-order nonlinearity is placed inside the cavity and is pumped by a laser to produce optical parametric amplification.The rotating end mirror is installed on the support S and can rotate about the z axis.The angular displacement of the rotating end mirror from its equilibrium position ϕ 0 = 0 is denoted by ϕ.An external Gaussian laser beam (G) enters the optical cavity through the left fixed mirror.The value of the topological charge of each beam is shown.
In a frame rotating at the incident Gaussian laser frequency ω l , the dynamics of the studied system is described by the following Hamiltonian: where ∆ 0 = ω c − ω l is the detuning of the cavity field with respect to the input laser; a (a † ) is the boson annihilation (creation) operator of the LG-cavity mode; g is the optorotational interaction strength determined by g = g ϕ h 2Iω m , with I = 1 2 mR 2 being the moment of inertia of the rotating end mirror about the z axis and R being the radius of the rotating end mirror; G and θ are the parametric gain and phase of the OPA, respectively; and ε is the amplitude of the incident Gaussian beam determined by ε = 2κ℘ hω l , with ℘ being the power of the incident Gaussian beam and κ being the cavity decay rate.In Equation (1), the fifth term represents the coupling between the LG-cavity field and the OPA.We assume that the degenerate OPA is driven by a pump field at frequency 2ω l , and the three-wave mixing process in the OPA can result in the production of the signal and idler fields with identical frequencies ω l .
Using the Heisenberg equation of motion and considering the corresponding damping terms, we obtain the equations describing the motion of the optical and mechanical modes: where we have included the input vacuum noise operator a in entering the cavity through the left end mirror and the thermal noise operator b in of the rotating end mirror associated with the Brownian motion of the rotating end mirror.Both the noise operators a in and b in have zero mean values.The correlation functions for the noise operators a in and b in in the time domain are given by where − 1} is the thermal phonon number of the rotating end mirror, k B is the Boltzmann constant, and T is the temperature of the surrounding thermal environment.At the steady state, the mean values of the operators a and b are found to be where a s and b s are the amplitudes of the LG-cavity field and the rotating end mirror at steady state.

Quantum Fluctuations
We assume that a strong Gaussian beam drives a single-mode LG-cavity field.In this case, the absolute amplitude |a s | (|b s |) of the cavity (mechanical) mode is much larger than 1; thus, the two operators a and b can be replaced by a = a s + δa and b = b s + δb, respectively, where δa and δb are small quantum fluctuations.We only keep the first order in the quantum fluctuations and obtain the following linearized equations for the fluctuations δa and δb: where g ′ = g + g ck b s .
For convenience, we introduce the quadrature fluctuations of the mechanical and cavity modes as δϕ = 1 (δa + δa † ), and (δa − δa † ), and the quadrature fluctuations of the input noises as The equations of motion for the quadrature fluctuations are given by where U(t) represents the vector of the quadrature fluctuations of the mechanical and cavity modes, N(t) denotes the vector of the input noises, their transposes are and the 4 × 4 matrix A has the form where . The eigenvalues of the matrix A can be used to determine the dynamic stability of the considered system.In order to stablize the system, the real parts of all the eigenvalues of the matrix A must be negative.According to the Routh-Hurwitz stability criterion [26], we find the three stability conditions of the system: ×(σ 2 + ρ 2 ) ωm > 0, In the following numerical simulations, we choose the system parameters to satisfy the above stability conditions.We assume that the incoming noises a in and b in are Gaussian noises, whose mean values are zeros.The time-evolution equations for the fluctuation operators δa and δb have been linearized.Thus, the resulting state of the system is also a Gaussian state and can be fully described by the 4 × 4 covariance matrix V whose elements are given by Based on Equation ( 6), we find the evolution equations for the covariance matrix V: where the diffusion matrix D can be calculated with the aid of the correlation functions of the noise operators The covariance matrix V of the system at the steady state can be obtained by solving the following Lyapunov equation [27]: For the Gaussian continuous-variable bipartite system, we can use the logarithmic negativity E N proposed in Ref. [28] to measure the degree of the stationary entanglement between the LG-cavity field and the rotating end mirror, which is given by where In Equation ( 14), the 2 × 2 matrices V m and V c represent the variances of the rotating end mirror and the LG-cavity field, respectively, and the 2 × 2 matrix V cm describes the correlation between the rotating end mirror and the LG-cavity field.Only when the logarithmic negativity E N is not equal to zero are the LG-cavity field and the rotating end mirror entangled.
For the Gaussian continuous-variable bipartite system, there are some other measures of entanglement, such as the Duan's inseparabiltiy criterion, which is based on the calculation of the total variance of a pair of Einstein-Podolsky-Rosen-type operators [29].

The Stationary Entanglement of the LG-Cavity Field with the Rotating End Mirror in the Presence of the Nonlinear CK Interaction and the OPA
In the following, we analyze the impacts of the cavity detuning ∆ 0 , the parametric gain G and phase θ of the OPA, the strength g ck of the nonlinear CK interaction, the Gaussian laser power ℘, the topological charge l of the cavity field, the mass m of the rotating end mirror, and the ambient temperature T on the stationary entanglement of the cavity field with the rotating end mirror.
The parameters we use in the numerical calculations are similar to those in Ref. [11], which studies the entanglement of an intracavity LG field with a rotating end mirror: the wavelength of the incident Gaussian laser is λ = 810 nm, the cavity length is L = 1 mm, the cavity field decays at a rate of κ = 2π × 1.5 MHz, the rotating end mirror has the radius R = 10 µm, the resonance frequency ω m = 2π × 10 MHz, the mechanical quality factor Q m = 2 × 10 6 , and the damping rate γ m = 2π × 5 Hz.Thus, the studied system is operating in the resolved-sideband limit ω m ≫ κ.
Figure 2 plots the logarithmic negativity E N against the normalized cavity detuning ∆ 0 /ω m for different gains G of the OPA when m = 100 ng, ℘ = 1 mW, θ = 0, 0.2π, l = 15, g ck = 10 −3 g, and T = 0.1 K.When the driving phase θ of the OPA is 0 (Figure 2a), the system stability conditions require ∆ 0 /ω m ≥ 0.117, 0.158, 0.193 for G = 0.1κ, 0.2κ, 0.3κ, respectively, and the corresponding maximum logarithmic negativity E N values just before the unstable regime are 0.250, 0.266, 0.297, respectively.When the driving phase θ of the OPA is set to 0.2π (Figure 2b), the system stability conditions for G = 0.1κ, 0.2κ, 0.3κ restrict the normalized cavity detuning to the range of ∆ 0 /ω m ≥ 0.100, 0.135, 0.168, respectively, the corresponding maximum logarithmic negativity E N just before the unstable regime is 0.252, 0.316, 0.237, respectively.In the absence of the OPA (G = 0), the system stability conditions require ∆ 0 /ω m ≥ 0.023, and we observe that the logarithmic negativity E N increases from 0 to a peak value of about 0.05 and then gradually decreases to 0 as the cavity detuning ∆ 0 increases.Thus, the maximum entanglement in the presence of the OPA far exceeds that in the absence of the OPA (G = 0), indicating that the nonlinear gain G of the OPA significantly enhances the steady-state entanglement.For a fixed phase θ of the OPA, as the gain G of the OPA increases, the regime where the entanglement occurs shifts towards larger cavity detunings ∆ 0 , and the amount of entanglement decreases with increasing cavity detuning ∆ 0 .In Figure 2b, when θ = 0.2π, it is noted that the maximum entanglement does not occur at the highest gain G of the OPA.This shows the influence of the driving phase θ of the OPA on the entanglement, which is discussed further in the following.Figure 3 plots the logarithmic negativity E N against the normalized gain G/κ of the OPA for different strengths g ck of the nonlinear CK interaction when m = 100 ng, ℘ = 1 mW, θ = 0, 0.2π, l = 15, ∆ 0 = 0.12ω m , and T = 0.1 K.When θ = 0 (Figure 3a), for g ck = 0.8 × 10 −3 g, 0.9 × 10 −3 g, 10 −3 g, 1.1 × 10 −3 g, when the gain G of the OPA satisfies G/κ ≤ 0.156, 0.132, 0.108, 0.085, respectively, the system remains stable.For g ck = 0.8 × 10 −3 g, 0.9 × 10 −3 g, 10 −3 g, 1.1 × 10 −3 g, as the gain G of the OPA increases, the logarithmic negativity E N increases and reaches the maximum values of about 0.263, 0.302, 0.293, and 0.300, respectively, just before entering the unstable regime.In Figure 3b (θ = 0.2π), for g ck /g = 0.8 × 10 −3 , 0.9 × 10 −3 , 1.0 × 10 −3 , 1.1 × 10 −3 , the stability condition imposes the limitation G/κ ≤ 0.206, 0.181, 0.155, 0.129, respectively.The maximum logarithmic negativity E N values before approaching the unstable regime are about 0.273, 0.305, 0.279, and 0.292, respectively.This demonstrates that the gain G of the OPA can effectively enhance the steady-state entanglement between the cavity field and the rotating end mirror.This is due to the fact that increasing the gain G of the OPA leads to a larger photon number |a s | 2 in the cavity field and a larger phonon number |b s | 2 in the rotating end mirror, as shown in Figure 4, which plots the number |a s | 2 of the photons in the cavity field and the number |b s | 2 of the phonons in the rotating end mirror at the steady state against the normalized gain G/κ of the OPA when m = 100 ng, ℘ = 1 mW, θ = 0, l = 15, ∆ 0 = 0.12ω m , and g ck = 10 −3 g.A larger photon number |a s | 2 in the cavity field increases the coupling between the cavity field and the rotating end mirror.Meanwhile, a larger photon number |a s | 2 in the cavity field and a larger phonon number |b s | 2 in the rotating end mirror give rise to a stronger CK coupling between the cavity field and the rotating end mirror.Moreover, in Figure 3, for fixed values of θ and G, the larger the CK coupling strength g ck , the larger the entanglement, but, for a fixed value of θ, with increasing the CK coupling strength g ck , the stable regime becomes narrower.Comparing Figure 3a with Figure 3b, it is evident that changing the driving phase θ of the OPA significantly alters the regime of steady-state entanglement.Figure 5 plots the logarithmic negativity E N against the phase θ/π of the OPA for some different gains G of the OPA when m = 100 ng, ℘ = 1 mW, l = 15, ∆ 0 = 0.12ω m , g ck = 10 −3 g, 1.1 × 10 −3 g, and T = 0.1 K. Without the OPA (G = 0), for g ck = 10 −3 g, 1.1 × 10 −3 g, when the phase θ of the OPA is changed from 0 to 2π, the values of the logarithmic negativity E N remain unchanged, which are about 0.011, 0.019, respectively.Next we look at the case with the OPA.When g ck = 10 −3 g (Figure 5a), for G = 0.1κ, 0.2κ, 0.3κ, the system is stable when the phase θ of the OPA satisfies θ ∈ [0, 1.71π] [1.91π, 2π], θ ∈ [0.3π, 1.31π], and θ ∈ [0.45π, 1.16π], respectively.When g ck = 1.1 × 10 −3 g (Figure 5b), for G = 0.1κ, 0.2κ, 0.3κ, the system is stable when the phase θ of the OPA satisfies θ ∈ [0.1π, 1.52π], θ ∈ [0.35π, 1.27π], and θ ∈ [0.47π, 1.14π], respectively.Thus, for a fixed strength g ck of the nonlinear CK interaction, increasing the gain G of the OPA makes the system stable in a narrower range of phases θ.And when the gain G of the OPA is fixed, for g ck = 1.1 × 10 −3 g, the system is stable in a narrower range of phases θ compared to that for g ck = 10 −3 g.When g ck = 10 −3 g (Figure 5a), for G = 0.1κ, 0.2κ, 0.3κ, the logarithmic negativity E N takes its maximum values of about 0.320, 0.254, and 0.171 at θ = 1.91π, 0.3π, 0.45π, respectively.When g ck = 1.1 × 10 −3 g (Figure 5b), for G = 0.1κ, 0.2κ, 0.3κ, the logarithmic negativity E N takes its maximum value of about 0.274, 0.190, and 0.249 at θ = 0.1π, 1.27π, 0.47π, respectively.Hence, for a fixed value of the CK coupling strength g ck , it is noted that the presence of the OPA can enhance the entanglement since the logarithmic negativity E N in the presence of the OPA can be larger than that in the absence of the OPA, but the maximum entanglement does not increase with an increasing gain G of the OPA, and the maximum entanglement happens at a phase θ of the OPA at which the system is close to the unstable regime.Figure 6 depicts the logarithmic negativity E N versus the normalized strength g ck /(10 −3 g) of the nonlinear CK interaction for different gains G of the OPA when m = 100 ng, θ = 0, 0.2π, ℘ = 1 mW, l = 15, ∆ 0 = 0.12ω m , and T = 0.1 K.When θ = 0 (Figure 6a), for different OPA nonlinear gains G = 0, 0.1κ, 0.2κ, 0.3κ, the system stability requires g ck /(10 −3 g) ≤ 1.48, 1.03, 0.63, 0.312.As the strength g ck of the nonlinear CK interaction increases, the logarithmic negativity E N for different G values rapidly rises from 0, reaching a maximum value of about 0.299, 0.265, 0.213, 0.104, respectively, just before the system approaches instability.Similarly, in Figure 6b (θ = 0.2π), for G = 0.1κ, 0.2κ, 0.3κ, the system stability conditions require g ck /(10 −3 g) ≤ 1.48, 1.20, 0.82, 0.44, respectively.As the strength g ck of the nonlinear CK interaction increases, the maximum logarithmic negativity E N values before approaching instability are about 0.299, 0.270, 0.235, 0.106, respectively.Thus, for a fixed phase θ of the OPA, with an increasing gain G of the OPA, the instability is triggered at a weaker strength g ck of the nonlinear CK interaction.For a given nonzero gain G of the OPA, it is noted that the instability for θ = 0.2π is triggered at a stronger strength g ck of the nonlinear CK interaction compared to the case for θ = 0.For a fixed phase θ and gain G of the OPA, it is possible to improve the entanglement by increasing the strength g ck of the nonlinear CK interaction.Additionally, as the gain G of the OPA increases, the entanglement appears at a weaker strength g ck of the nonlinear CK interaction and reaches the maximum value at a weaker strength g ck of the nonlinear CK interaction, but the maximum entanglement decreases.Figure 6.The logarithmic negativity E N against the normalized strength g ck /(10 −3 g) of the nonlinear CK interaction for different gains G of the OPA when m = 100 ng, ℘ = 1 mW, l = 15, ∆ 0 = 0.12ω m , and T = 0.1 K. (a) θ = 0, (b) θ = 0.2π.The black-dotted, red-solid, blue-dashed, and green-dotdashed curves are for G = 0, 0.1κ, 0.2κ, 0.3κ, respectively.
Figure 7 plots the logarithmic negativity E N against the laser power ℘ for different gains G of the OPA when m = 100 ng, θ = 0, 0.2π, ∆ 0 = 0.12ω m , g ck = 10 −3 g, l = 15, and T = 0.1 K.When θ = 0 (Figure 7a), for G = 0, 0.1κ, 0.2κ, 0.3κ, the stability conditions require ℘ ≤ 1.42, 1.03, 0.68, 0.39 mW, and the corresponding driving power thresholds for generating entanglement are 0.36, 0.76, 0.59, 0.37 mW, the maximum logarithmic negativity E N values near the unstable regime are 0.326, 0.298, 0.236, and 0.098, respectively.Similarly, when θ = 0.2π (Figure 7b), the stability conditions require ℘ ≤ 1.42, 1.17, 0.85, 0.52 mW, respectively, and the corresponding driving power thresholds for generating entanglement are 0.36, 0.87, 0.73, 0.48 mW, respectively; then, the maximum logarithmic negativity E N values near the unstable regime are about 0.326, 0.259, 0.293, 0.168, respectively.For a fixed phase θ of the OPA, from Equation (4), it is noted that increasing the laser power ℘ can enhance the photon number |a s | 2 within the cavity; thus, it is possible to improve the steady-state entanglement by increasing the laser power ℘.Additionally, as the gain G of the OPA increases, the entire entanglement regime shifts towards lower laser powers, with the threshold power appearing at lower values.This indicates that it is possible to achieve stronger steady-state entanglement with lower laser powers.Figure 8 plots the logarithmic negativity E N against the topological charge l of the cavity field for different gains G of the OPA when m = 100 ng, θ = 0, 0.2π, ℘ = 1 mW, ∆ 0 /ω m = 0.12, g ck = 10 −3 g, and T = 0.1 K.When θ = 0 (Figure 8a), for G = 0, 0.05κ, 0.1κ, 0.15κ, the entanglement occurs only when the topological charge l is not less than the threshold values l c = 8, 12, 12, 11, respectively, the stable regime is found to be l ≤ 20, 17, 15, 13, and the maximum logarithmic negativity E N close to the unstable regime is about 0.189, 0.122, 0.170, 0.292, respectively.When θ = 0.2π (Figure 8b), for G = 0, 0.05κ, 0.1κ, 0.15κ, the entanglement appears only when the topological charge l is not less than the threshold value l c = 8, 13, 14, 13, respectively, and the stable regime is found to be l ≤ 20, 19, 17, 15, respectively, the maximum logarithmic negativity E N close to the unstable regime is about 0.188, 0.261, 0.195, 0.196, respectively.Notably, when the topological charge l is less than the threshold value l c , the optorotational coupling is not strong enough to generate the entanglement.In addition, the exchange of orbital angular momentum between photons and helical phase elements is fundamental for generating entanglement between the cavity and mechanical modes.Thus, it is possible that an increase in the topological charge l facilitates the enhancement of entanglement.It is found that a higher gain G of the OPA can result in the maximum entanglement occurring at a smaller topological charge l, with a more rapid increase in entanglement and a narrowing stable regime.For a fixed nonzero gain G of the OPA, the system for θ = 0.2π becomes unstable at a larger topological charge l compared to the case for θ = 0. Figure 9 shows the logarithmic negativity E N against the mass m of the rotating end mirror for different gains G of the OPA when ℘ = 1 mW, θ = 0, 0.2π, ∆ 0 /ω m = 0.12, g ck = 10 −3 g, l = 15, and T = 0.1 K.When θ = 0 (Figure 9a), for G = 0, 0.05κ, 0.1κ, 0.15κ, the system's stability conditions require m ≥ 54, 72, 96, 134 ng, respectively, and the logarithmic negativity E N reaches the maximum values of about 0.281, 0.222, 0.251, 0.258 near the unstable regime, respectively.When θ = 0.2π (Figure 9b), for G = 0, 0.05κ, 0.1κ, 0.15κ, the system's stability conditions require m ≥ 54, 62, 76, 98 ng, respectively, and the logarithmic negativity E N reaches the maximum values of about 0.281, 0.281, 0.251, 0.245 near the unstable regime, respectively.For fixed values of θ and G, when the mass m of the rotating end mirror increases, the entanglement decreases.When the phase θ of the OPA is fixed, for larger gain G of the OPA, the system becomes stable at a larger mass m of the rotating end mirror.For a fixed nonzero gain G of the OPA, the system for θ = 0.2π becomes stable at a smaller mass m of the rotating end mirror in comparison with the case for θ = 0.When the phase θ of the OPA is fixed, a larger gain G of the OPA can generate the entanglement between the optical mode and a rotating end mirror with a larger mass m.Even at around 150 ng (Figure 9a), the entanglement can still be achieved, providing valuable insights for experimental realization.Figure 10 plots the logarithmic negativity E N against the ambient temperature T for different gains G of the OPA when m = 100 ng, θ = 0, 0.2π, ℘ = 1 mW, l = 15, ∆ 0 /ω m = 0.12, and g ck = 10 −3 g.For θ = 0, 0.2π, the system is stable when G is less than 0.109κ, 0.156κ, respectively.Thus, the system for θ = 0.2π becomes unstable at a larger gain G of the OPA in contrast to the case for θ = 0.When θ = 0 (Figure 10a), for G = 0, 0.09κ, 0.1κ, 0.108κ, the logarithmic negativity E N takes its maximum values of about 0.012, 0.110, 0.172, and 0.295 at T = 0 K, respectively, and the logarithmic negativity E N is zero when the ambient temperature T is not less than 1.8 K, 8.3 K, 13.8 K, and 27.4 K, respectively.When θ = 0.2π (Figure 10b), for G = 0, 0.1κ, 0.15κ, 0.155κ, the logarithmic negativity E N takes its maximum value of about 0.012, 0.034, 0.198, and 0.281 at T = 0 K, respectively, and the logarithmic negativity E N is zero when the ambient temperature T is not less than 1.8 K, 1.9 K, 13.1 K, and 21.8 K, respectively.When the phase θ and gain G of the OPA are fixed, as the ambient temperature T increases, the entanglement decreases and eventually vanishes, which can be attributed to the detrimental effect of the thermal noise of the environment on the entanglement.For a fixed phase θ of the OPA, an increase in the gain G of the OPA results in a larger entanglement and allows the entanglement to persist at higher temperatures.This indicates that the gain G of the OPA not only enhances the steady-state entanglement, but also improves the system's ability to withstand decoherence in a thermal environment.In our previous work [20], it has been shown that the entanglement between the cavity field and the rotating end mirror with mass 100 ng can be improved by using the nonlinear CK interaction, and the entanglement (E N = 0.459) can be obtained with the nonlinear CK interaction (g ck = 10 −3 g), which is much larger than the entanglement (E N = 0.011) (Figure 5a) between the cavity field and the rotating end mirror with mass 100 ng in the presence of the nonlinear CK interaction (g ck = 10 −3 g) but in the absence of the OPA in this work.The reason is that a large topological charge l = 70 of the cavity field is used in Ref. [20] and a small topological charge l = 15 of the cavity field are used in this work.Furthermore, it has been shown in this work that a larger entanglement (E N = 0.320) (Figure 5a) between the cavity field and the rotating end mirror can be obtained under the combined action of the nonlinear CK interaction (g ck = 10 −3 g) and the OPA (G = 0.1κ, θ = 1.91π).
Finally, we discuss the possibility of the experimental feasibility of this proposal.The optorotational entanglement can be verified by measuring the logarithmic negativity E N , which requires one to measure the 10 independent elements of the correlation matrix V [27].The quadratures of the LG-cavity mode can be measured by homodyne detecting the cavity output field [27].And the quadratures of the rotating end mirror can be measured by homodyning the output field from the adjacent cavity formed by the rotating end mirror and a third fixed spiral phase element [27].On the other hand, the LG-cavity optorotational system has not been realized experimentally yet [12].However, it has been demonstrated experimentally that an LG laser beam with a topological charge of l = 100 can be achieved by using spiral-phase mirrors [30], that the torsional frequency of a nanomechanical resonator can be up to 8.4 MHz [31], and that a micromechanical resonator with mass 25 ng, radius 15 µm, and mechanical quality factor 1.3 × 10 5 can be cooled to 135 mK [32].Thus, with the rapid development of fabrication technologies, it is possible to experimentally realize the optorotational entanglement in the LG-cavity optorotational system with the rotating mirror having mass 100 ng, radius 10 µm, frequency 2π × 10 MHz, and mechanical quality factor 2 × 10 6 at temperature T = 0.1 K.

Conclusions
To sum up, we have investigated the stationary entanglement between an LG-cavity field and a rotating end mirror in an LG-cavity optorotational system with a nonlinear CK interaction and a degenerate OPA.We show that the combination of the nonlinear CK interaction and the OPA have detrimental effects on the system stability so that the system becomes unstable at a lower Gaussian laser power, a smaller topological charge of the cavity field, and a larger mass of the rotating end mirror.Moreover, we show that the strength of the nonlinear CK interaction and the gain and phase of the OPA have strong impacts on the stationary entanglement.We find that the combination of the nonlinear CK interaction and the OPA can significantly increase the entanglement at lower Gaussian laser powers, the smaller topological charges of the cavity field, and the larger masses of the rotating end mirror and make the entanglement more robust against the ambient temperature compared to the case in the presence of only the nonlinear CK interaction [20].And the maximum entanglement appears near the unstable regime in the presence of the nonlinear CK interaction and the OPA.In contrast, the injection of two-level atoms into the optical cavity of an LG-cavity optorotational system can also significantly improve the optorotational entanglement and enhance the robustness of entanglement to the ambient temperature, but the maximum entanglement does not happen near the unstable regime in the presence of two-level atoms [12].Our findings have important implications for realizing the stationary entanglement experimentally.In the near future, we might study how the combination of the nonlinear CK interaction and the OPA affects the entanglement of multiple rotating mirrors in an LG-cavity optorotational system [33].

Figure 1 .
Figure 1.Schematic diagram of anLG-cavity optorotational system with a nonlinear CK interaction and a degenerate OPA.There is a two-level system (red) on the rotating end mirror to produce the nonlinear CK interaction between the LG-cavity field and the rotating end mirror.And a nonlinear optical crystal with a second-order nonlinearity is placed inside the cavity and is pumped by a laser to produce optical parametric amplification.The rotating end mirror is installed on the support S and can rotate about the z axis.The angular displacement of the rotating end mirror from its equilibrium position ϕ 0 = 0 is denoted by ϕ.An external Gaussian laser beam (G) enters the optical cavity through the left fixed mirror.The value of the topological charge of each beam is shown.

Figure 4 .
Figure 4. (a) The number |a s | 2 of the photons in the cavity field at the steady state against the normalized gain G/κ of the OPA.(b) The number |b s | 2 of the phonons in the rotating end mirror at the steady state against the normalized gain G/κ of the OPA.The parameters are m = 100 ng, ℘ = 1 mW, l = 15, ∆ 0 = 0.12ω m , g ck = 10 −3 g, and θ = 0.