Selective cooling and squeezing in a lossy optomechanical closed loop embodying an exceptional surface

A closed-loop, lossy optomechanical system consisting of one optical and two degenerate mechanical resonators is computationally investigated. This system constitutes an elementary synthetic plaquette derived from the loop phase of the intercoupling coefficients. In examining a specific quantum attribute, we delve into the control of quadrature variances within the resonator selected through the plaquette phase. An amplitude modulation is additionally applied to the cavity-pumping laser to incorporate mechanical squeezing. Our numerical analysis relies on the integration-free computation of steady-state covariances for cooling and the Floquet technique for squeezing. We provide physical insights into how non-Hermiticity plays a crucial role in enhancing cooling and squeezing in proximity to exceptional points. This enhancement is associated with the behavior of complex eigenvalue loci as a function of the intermechanical coupling rate. Additionally, we demonstrate that the parameter space embodies an exceptional surface, ensuring the robustness of exceptional point singularities under experimental parameter variations. However, the pump laser detuning breaks away from the exceptional surface unless it resides on the red-sideband by an amount sufficiently close to the mechanical resonance frequency. Finally, we show that this disparate parametric character entitles frequency-dependent cooling and squeezing, which is of technological importance.

A considerable body of these efforts has been devoted to investigating the cooling of mechanical resonators.Various schemes exist for this goal, such as backaction cooling [31,32], and feedback cooling [5,10].Another effective choice utilizes the so-called resolved-sideband regime, which is reached when the cavity decay rate is much lower than the frequency of the mechanical oscillator [6,7,33].Accordingly, performed experiments revealed a sizable amount of cooling of the mechanical vibrations in this regime [31,[34][35][36][37][38][39].Going beyond a single resonator, simultaneous cooling of two mechanical resonators coupled to an optical cavity in the resolvedsideband regime was theoretically suggested [40].It has importance from the standpoint of quantum coherence and also as an experimental realization of cooling of hybridized modes in the mechanical degeneracy regime [41].However, if two mechanical resonators are coupled to an optical cavity, they form a so-called Λ scheme, which generates dark and bright degenerate modes [42].Though * bulutay@fen.bilkent.edu.truseful for storing quantum information, the dark mode creates an obstacle for multiple cooling [32].By means of a phase-dependent phonon-exchange loop-coupling, the dark mode effect is alleviated for specific phase values, and multiple mechanical resonators can be cooled to their ground-state [43].Alternatively, the dark mode effect can be eliminated for multimode optomechanical cooling by introducing an auxiliary cavity or mechanical mode to the system as a substitute cooling channel [44][45][46].
A recently flourishing, yet seemingly disparate topic is artificial gauge potentials [78,79], bestowing features arXiv:2307.09851v4[quant-ph] 26 Jun 2024 such as nonreciprocal photon transport [80][81][82], on-chip optical nonreciprocity, and optical isolation [83,84].Optical nonreciprocity refers to the phenomenon in which the transmission of light is different in one direction compared to the opposite, which yields various applications such as isolators [85], circulators [86], and directional amplifiers [87].In optomechanics, nonreciprocity offers a promising avenue for manipulating thermal fluctuations, such as for cooling phononic resonators [88], thermal-noise cancellation [89], mode routing and thermal management [90].Such nonreciprocity can be created by breaking the time-reversal symmetry in the system, which is one of the prime uses of synthetic gauge fields [91].A frequency difference between mechanical modes or cavity field modulation is a reliable method for creating synthetic fields [92,93].They enable tunable cooling in the optomechanical system by introducing a phase-dependent coupling constant [43,94].Most recently, strong mechanical squeezing has been proposed by breaking the dark-mode effect with the synthetic-gaugefield [95].
Captivating all these fields, a revolutionary wave of non-Hermitian physics is currently underway, offering unique advantages empowered by the so-called exceptional point (EP) singularities [96][97][98].Naturally, non-Hermitian phenomena have become an appealing resource, so far largely capitalized in photonics, with unprecedented consequences such as, loss-induced revival of lasing [99], phonon lasing [100], boosting optomechanical interactions with high-order EPs [101], enhanced quantum sensing [102].As a few other outstanding yet highly relevant studies for our work, mechanical cooling is proposed in a non-Hermitian system [103], also with a synthetic gauge field [94], where the average phonon occupation number is minimized at the EP, and chiral, non-Hermitian tunable mechanical squeezing in nano-optomechanical networks are experimentally realized [104].Another notable development in this front is the concept of an exceptional surface on which the EPs move [105].Its presence assures insensitivity to variations over the parameters that define the exceptional surface.As such, the existence of an EP does not require unrealistic fine-tuning of parameters.Soon after the exceptional surfaces were proposed [106,107], they were experimentally demonstrated for two-mode systems of magnon polariton [108] and optical waveguides [109].
In line with this progress, the aim of this work is to present a unified treatment of both cooling and squeezing of the mechanical quadratures controllable through a closed-loop phase of a photonic cavity coupled to two intercoupled lossy mechanical resonators.Compared to our recent work, where we studied loop phase control of optomechanically-induced transparency in a gain-loss non-Hermitian setting, here we tackle the control of quantum phenomena of quadrature variances and squeezing, in contrast to probe transmission, which is essentially governed by classical mean fields [110].Our numerical framework constitutes the Floquet formalism [111][112][113] for steady-state or modulated evolution.First, mechanical cooling is explored based on an experimentally attainable parameter set.Moreover, mechanical squeezing can be accomplished by imposing an upper-sideband amplitude modulation of the cavity pumping laser while the loop-phase selects the resonator with maximum squeezing.We demonstrate that non-Hermiticity plays an intricate role, as determined by the proximity to two EPs derived from the upper-and lower-frequency polaritonphonon modes.Next, we show that the parameter variations of this three-mode optomechanical system guide the EPs over an exceptional hypersurface, making them robust against fabrication errors and fluctuations in the experimental environment.The exception is the detuning, as it needs to be on the red-sideband by an amount close to the mechanical resonance frequency; otherwise, with further deviation, the system breaks away from the exceptional surface.Favorably, this behavior bestows a frequency-dependent squeezing as a function of laser detuning.
The paper is organized as follows.In Sec.II, we introduce our model and present its theoretical analysis.In Sec.III, we discuss possible experimental realization schemes and our choice of parameters.This is followed by our results in Sec.IV, divided into three subsections: mechanical cooling, exceptional surfaces, and squeezing.Our main conclusions are drawn in Sec.V. Appendix A contains the steady-state values of the cavity and mechanical modes; Appendix B derives the covariance matrix expressions for the Floquet formalism; Appendix C derives integration-free computation of steady-state variances.

II. THEORY A. Optomechanical Network Hamiltonian
We consider the simplest closed optomechanical network consisting of a photonic cavity coupled to two mechanical resonators via coupling rates g 1 and g 2 as shown in Fig (1).The mechanical resonators are bilinearly intercoupled with a rate µ, and both have equal resonance angular frequency ω m , and damping rates γ 1 and γ 2 .The photonic cavity is driven, in the most general case, with an amplitude-modulated laser with carrier frequency ω L and amplitude ε L (t); when applied, the modulation frequency is Ω = 2π/τ so that , with d n = ε n /ε 0 being the modulation amplitude depth associated with the sideband n where amplitude ε 0 is defined in terms of the laser power, ε 0 = 2P ℏω L .The Hamiltonian (ℏ = 1) in the rotating frame of the pump laser at frequency ω L is where ∆ = ω cav − ω L is the frequency detuning between the cavity and input laser, â (â † ), and b1 ( b † 1 ) and b2 ( b † 2 ) are the annihilation (creation) operators of the optical, first and second mechanical resonator modes, respectively, κ is the cavity decay rate, and η is the cavity coupling parameter (escape efficiency).The non-Hermiticity of the system follows from the lossy nature of all the resonators.
FIG. 1. Closed-loop interaction optomechanical system composed of a photonic cavity with the relevant resonance at ωcav, and two mechanical resonators with identical frequencies, ωm.Loss rates are indicated with wavy arrows.The cavity is pumped with a laser (in the case of squeezing also modulated) having a carrier frequency ωL and amplitude εL(t).

B. Closed-loop phase
The three complex coupling coefficients, , and µ = |µ| e iϕµ , in Eq. ( 1) comprise a closed-loop phase, ϕ ℓ := −ϕ 1 +ϕ 2 +ϕ µ [110].However, as we discuss in Sec.III, retaining g 1 , g 2 ∈ R while imparting ϕ µ → ϕ ℓ is the viable implementation from a practical standpoint.Notably, ϕ ℓ , just like a Peierls phase, threads synthetic magnetic flux through this triangle plaquette [114].ϕ ℓ is a gauge-invariant quantity having observable consequences as it controls which mechanical resonator is to be predominantly cooled or squeezed.In the case of identical mechanical resonators, a crucial consequence of this choice is its invariance under simultaneous swapping of the two mechanical resonators, 1 ←→ 2 along with ϕ ℓ → −ϕ ℓ = 2π − ϕ ℓ (see Fig. 2).In other words, one can toggle between the selected mechanical resonators via a phase reversal.We expound on this symmetry in the Results sections.

FIG. 2.
In the closed-loop coupling scheme, swapping the two identical mechanical resonators 1 ←→ 2, is equivalent to

C. Quantum Langevin equations
To account for input noise and losses, we switch to quantum Langevin equations [115], which describe the dynamics of the cavity and mechanical modes as where âin (t), b1,in (t) and b2,in (t) are zero-mean cavity and mechanical input noise operators, respectively.They satisfy the following correlation functions (displaying only the non-zero ones) under the Markovian-reservoir assumption [30], where j = 1, 2, n a , and n m are the mean occupancy of the cavity and mechanical baths, respectively.For cooling and squeezing purposes, our primary focus is on quantum fluctuations.Thus, we linearize the system by writing cavity mode â and mechanical modes b1 , b2 as a sum of the classical mean value and quantum fluctuations operators ‫(א‬t) → ⟨ ‫(א‬t)⟩ + δ ‫(א‬t).The remaining details of the steady-state values are provided in Appendix A.
The quantum fluctuations around the classical mean values, represented by δ ‫(א‬t), have the following equations of motion Next, we switch to experimentally accessible quadrature operators of position and momentum, which are expressed in terms of the fluctuation operators as and the corresponding quadrature noise operators are Finally, the position-momentum quadrature fluctuations of the cavity and mechanical modes can be cast in the form [76] where R(t) = [δ Xa , δ Ŷa , δ Xb1 , δ Ŷb1 , δ Xb2 , δ Ŷb2 ] T , with the superscript T indicating vector or matrix transpose and the dot, Ṙ represents time derivative.M(t) is a 6×6 time-dependent matrix with its elements (not included here) readily extracted from the equations of motion of δ ‫(א‬t).N(t) is the noise operator vector and defined as We will solve Eq. ( 6), which is a first-order inhomogeneous differential equation.

D. Quadrature variances
The formal solution of Eq. ( 6) is where G(t) satisfies Ġ(t) = M(t)G(t) subject to the initial condition G(0) = I, with I being the identity matrix [76].In order to investigate the mechanical squeezing and cooling, we need quadrature fluctuations of mechanical resonators.Hence, we introduce the covariance matrix V(t) to analyze the dynamics of this optomechanical system where its entries are defined as for i, j = 1, 2, . . ., 6. From the definition of the covariance matrix and Eq. ( 8), we obtain where Entries V 33 , V 44 , V 55 and V 66 give the position and momentum variances of first and second mechanical res-onators, respectively.Squeezing beyond the vacuum state occurs when −10 log 10 [V ii /0.5] > 0 dB, since for the vacuum state has δX 2 = δY 2 = 0.5 .

E. Floquet Analysis
The formal solution gives the variance as a function of time.An alternative is the Floquet method, which yields variance when t → ∞, i.e., in the steady state regime.Due to the modulation with frequency Ω, in the absence of noise the dynamical variables, R(t), are also periodic with Ω in the steady-state i.e., R(t) = R(t + 2π/Ω).This allows us to expand these dynamical variables into the Fourier series whereas Fourier coefficients are actually time-dependent due to noise.Introducing the Fourier transformation for these Fourier series coefficients, we obtain We can also expand the so-called drift matrix M(t) into its harmonics as Eq. ( 6) becomes For this equality to hold, terms inside the parenthesis on the left and right-hand sides must be equal, where Rn (ω−Ω) = Rn−1 (ω).In this case, we have the following equality for Eq. ( 6), Here, the advantage is the removal of time dependence in the M matrix.Now, we have a time-independent but infinitely coupled set of algebraic equations among different harmonic contributions.We write Eq. ( 15) as Eq. ( 16) can be written as P(ω) R(ω) = n(ω).Inevitably, we need to limit Fourier components to a maximum size, ±N , which we take as N = 2 after checking its convergence.Forthwith, we can solve R(ω) by matrix inversion as R(ω) = T(ω)n(ω) where T(ω) = P −1 (ω).Importantly, the entries of P in Eq. ( 16) are also matrices, and the overall dimension is 6(2N + 1) × 6(2N + 1).Dimension of R(ω) is 6 × (2N + 1).In Appendix B, we extend this formalism to the covariance matrix V(t).To limit the technical level of the main text, the remaining expressions of this formalism are deferred to Appendix B.Moreover, in Appendix C, we provide an integration-free calculation of steady-state variances in the absence of modulation.

III. POSSIBLE EXPERIMENTAL REALIZATION AND THE PARAMETER SET
In order to substantiate the applicability of our work, we first argue on its possible experimental realization.Here, the main challenge is achieving the closed-loop phase, ϕ ℓ .This being a gauge invariant quantity, it can be imparted to any one of the couplings or divided among them while leading to the same physical consequence [110].We believe that the most viable option is incorporating the plaquette phase solely to the intermechanical coupling µ, where it is mediated by an optical link with its phase being locked to the optical cavity pumping laser up to an offset that is tunable by electrooptic and thermo-optic means [116,117].
We should note that this full parameter set was not pre-optimized for mechanical cooling or squeezing purposes [28].In analyzing any of these parameters' sensitivity, we designate the numerical values above with sub-scripts such as κ c , P c , etc.For this data set, there are two EPs for the intermechanical coupling rate, µ occurring at |µ EP,1 | ≃ 52.5 (γ 1 + γ 2 ) and |µ EP,2 | ≃ 80.45 (γ 1 + γ 2 ), which we utilize in our following analysis.For each |µ EP | there are two distinct ϕ ℓ values corresponding to opposite chirality EPs, such that the one with ϕ ℓ < π is the clockwise and ϕ ℓ > π is the counterclockwise EP.Having κ > γ 1,2 > 0, this makes up a loss-loss system.As shown by Özdemir et al., as long as the mode losses are nonuniform (as in our case) under a gauge transformation, it can be mapped to the prototypical gain-loss non-Hermitian model [98].

A. Mechanical Cooling
We begin our analysis with the optomechanical cooling process.Under red-detuned pumping, the laser photons to be admitted into the cavity need to gain the deficient energy from the phonons of the mechanical resonators, specifically through the â †b 1,2 terms of the Hamiltonian [cf Eq. ( 1)].The steady-state mean phonon number (n(∞)) is found with the Integration-free calculation explained in Appendix C. In Fig. 3, we plot the final mean phonon numbers of each of the mechanical resonators as a function of the closed-loop phase, ϕ ℓ for either of the exceptional-point coupling values, |µ| = |µ EP,1,2 |.Either mechanical resonator can be cooled down to ground state selectively depending on ϕ ℓ , which is indicative of the clockwise versus counterclockwise association of the ϕ ℓ = π/2 and 3π/2 values.That is, for ϕ ℓ = π/2 and 3π/2, the second or first mechanical resonator is respectively favored in cooling [43,94].In the same vein, when ϕ ℓ = 0 or π, resonator discrimination is lost as the timereversal symmetry is restored.Previous studies considered specific synthetic phases of ϕ ℓ = π/2 and 3π/2, and determined peak performance at the EPs [43,94].This fact signals that non-Hermiticity in the system, which is responsible for the EPs, promotes cooling as well.One of these works explained this by relating EPs to field localization in the resonator with less loss [99,103], while the other study emphasized the unidirectionality of phonon transport for this scheme [94].To elaborate on these findings, we examine where the optimal mechanical cooling lies with respect to both intermechanical resonator coupling and global phase.In Fig. 4(a), we plot mean phonon number of the second mechanical resonator at the ϕ ℓ ≃ 0.49π with respect to intermechanical coupling, which harbors the minimum phonon number as extracted from Fig. 4(b) which is marked with a cross.We observe that the minimum final mean phonon number of the second mechanical resonator [n 2,min (∞) ≃ 0.591] appears at |µ| = 1.12|µEP,1 | and ϕ ℓ ≃ 0.49π, which is shifted from ϕ ℓ = π/2, thus not exactly at an EP, unlike [94] that registers optimal cooling at the EP.As a side remark, much higher EPinduced cooling was reported for the open-loop geometry [103]; however, it lacks mechanical resonator selectivity with the plaquette phase, as in here or [94].In order to gain more insight on why the range |µ EP,1 | ≤ |µ| ≤ |µ EP,2 | hosts optimal cooling, we plot in Fig. 5 the complex upper half-plane eigenvalues (z (|µ|, ϕ ℓ ) = α + iω, ω > 0) of the stationary drift matrix, M (0) at six different |µ| values while continuously varying ϕ ℓ ∈ [0, π]; for ϕ ℓ ∈ [π, 2π] the same loci are retraced backwards.The vertical axis designates oscillation angular frequency, and the horizontal axis the damping rate for α < 0. Subplots (b) and (e) in Fig. 5 display the two EPs |µ EP,1 | and |µ EP,2 |, where both real and imaginary parts of two modes as well as their eigenvectors coincide.Markedly, these occur at phase angles away from π/2, unlike in our optomechanically induced transparency study [110]; unfortunately, this hampers an analytical treatment of the EPs.On the big picture, as |µ| is increased from below |µ EP,1 | to above |µ EP,2 | the bare modes undergo successive hybridizations, as we delineate below.Starting from |µ| ≪ |µ EP,1 |, as the intermechanical coupling |µ| is ramped up, the two originally degenerate phonon modes at ω m split, forming lower-and upper-frequency phonon hybridizations.Concomitantly, the radiation pressure couplings between the optical and two mechanical resonators induce phonon-polariton characters.For |µ| ∼ |µ EP,1 |, the upper-frequency polariton branch experiences an EP around 1.04 ω m .As |µ| is increased to about |µ EP,2 |, this time the lower-frequency polariton branch EP is formed around 0.94 ω m .At either of these EPs, the two approaching skewed phononpolariton states become completely parallel, which is the hallmark of non-Hermitian systems [98].The range |µ| > |µ EP,2 | causes strongly intercoupled mechanical modes which are weakly coupled to the optical mode, opposite to the case in |µ| < |µ EP,1 |.Thereby, the two extreme subplots (a) and (f) are the mirror reflections of one another with respect to a midway vertical line.Overall, the most critical aspect of the eigenvalue loci in Fig. 5 for mechanical cooling purposes is the polaritonphonon mode coupling displayed with black arrows in plots (c) and (d), which links to the cold optical cavity mode.Once again, for ϕ ℓ < π/2 (ϕ ℓ > π/2), the second (first) mechanical resonator is preferentially cooled.Final phonon number, n2 (∞) Final phonon number, n2 (∞) FIG. 6. Steady-state mean phonon number of the second mechanical resonator (ϕ ℓ = π/2) with respect to κ and pump power, normalized to κc/(2π) =900 MHz, and Pc = 0.25 mW.
Proceeding with the impact of the other parameters on the final mean phonon number, in Fig. 6, we investigate cavity decay rate and pump power normalized to their values in the original parameter set κ c and P c , respectively.The rapid increase in the final mean phonon number for κ < κ c (upper panel) arises from the fact that the photonic cavity becomes increasingly off-resonant compared to the broadening of the mechanical resonators' linewidths γ 1,2 .On the other hand, the slow rise as κ > κ c is the manifestation of gradual departure from the beneficial resolved-sideband regime.Likewise, in the lower panel, we observe that pump power has an optimum value slightly above P c .This sensitivity is due to the EP dependency on pump laser power.As the latter changes, the proximity to either of the EPs is adversely affected.
In Fig. 7, we plot the steady-state mean phonon number of the second mechanical resonator with respect to the ambient temperature (taken as equal for the cavity and mechanical resonators) from cryogenic to room temperature.In accordance with the Bose-Einstein distribution, the corner thermal energy beyond which n2 (∞) rapidly increases corresponds to mechanical resonance, ω m .The ground-state cooling regime is marked by the shaded area.This is approached at higher temperatures depending on the closed-loop phase, with the optimal value being close to π/2 (see Fig. 3).Final phonon number, n2 (∞) FIG. 7. Steady-state mean phonon number of the second mechanical resonator with respect to reservoir temperature at the exceptional point |µEP,1| ≃ 52.5 (γ1 + γ2) for three different ϕ ℓ values.

B. Exceptional Surfaces
To complement our parameter analysis, we bring in a recently perceived concept of an exceptional surface [105].As we discussed in Figs. 4 and 5, having an EP has important physical consequences.However, if its presence is limited to a very finely chosen parameter set, it remains only a mathematical serendipity.In this regard, the formation of exceptional surfaces is highly welcomed to alleviate such experimental concerns [105].Here, we systematically trace EPs over a range of different parameters to ratify that our system actually possesses exceptional surfaces.As depicted in Fig. 8, we track the two EPs in relation to distinct parameter combinations.For the subplots (a)-(c), we see that they form a hypersurface contingent on the two specified parameters.As such, the system exhibits robustness against the variations of pump power, radiation pressure couplings, and mechanical losses so that it will reside on the exceptional surface.The exception to this is the behavior under detuning (∆) displayed in Fig. 8 (d).The EPs are found only in a very restricted ∆ range around ω m .In other words, a change in ∆ corresponds to a dimension perpendicular to the exceptional surface.Being limited to three-axis visualization, in Fig. 5, we can only plot the magnitudes of µ EP .Because of this, in subplot (d), the two surfaces cross each other, which seemingly suggests the presence of a third-order EP where all three eigenvalues coalesce [122].However, upon further examination, we observed that in Fig. 8 (d), two EPs always have different phases when their magnitudes are the same so that they each remain as second-order EPs.

C. Mechanical Squeezing
To obtain mechanical squeezing beyond the vacuum level, we apply amplitude modulation over the pump laser [74].We opt for the so-called upper sideband mod-ulation, i.e., ϵ 1 ̸ = 0, ϵ −1 = 0.As with the laser carrier component, ϵ 0 causing a shift in the EPs, a similar effect occurs under modulation.In particular, for the modulation depth of d = ε 1 /ε 0 = 0.5, the first EP moves to |µ| ≃ 50.83 (γ 1 + γ 2 ).Mechanical squeezing is achieved when the variance is reduced below the vacuum value of 1/2, as marked by the shaded region in Fig. 9. Here, we explore the effect of modulation depth, which is defined as the ratio d = ε 1 /ε 0 .The same control over the loop phase, ϕ ℓ , also persists for the mechanical squeezing, displayed in Fig. 9(a), and as expected, higher modulation depth enhances squeezing.Specifically, mechanical squeezing is reached for d > 0.3 (corresponding to 0 dB, marked with a dashed line), and towards d → 1 system starts to become numerically unstable.As in cooling, the loop phase angle for maximum squeezing again occurs away from the π/2 value and is modulation depthdependent.Fig. 9(b) indicates that for µ → 0 the mechanical squeezing cannot be attained.As the mechanical coupling constant increases, so does the degree of squeezing.However, squeezing gradually diminishes as the system moves away from the first exceptional point.In Fig. 10(a) and (b), we plot steady-state mean phonon number and the minimum variance of the second mechanical resonator as a function of detuning for different intermechanical coupling constants.Here, the parameter that breaks away from the exceptional surface, i.e., laser detuning, brings quite a beneficial feature, namely frequency-dependent cooling, and squeezing.In that respect, the importance of loop coupling is manifested by the absence of frequency dependence when intermechanical coupling is set to zero (µ = 0).In contrast, under both loop coupling and modulation, frequency-dependent squeezing is attained, which has recently gained importance, distinctly in the context of gravitational wave detection [51].As expected, maximum cooling and squeezing occurs close to µ EP and for the red-detuning of ∆ ≃ ω m and ∆ ≃ 1.1ω m , respectively.On either side, both the mechanical cooling and squeezing fall off as a function of the laser pump frequency.Beyond a certain detuning, cooling below ground state and squeezing below the vacuum level (marked by a shaded region) are lost.

V. THE VALIDITY RANGE OF THE NUMERICAL TECHNIQUES
We would like to end our discussion on which extremes our numerical approaches can fail.Foremost, the com-putation of the steady-state variances hinges on the stability of the optomechanical response.For cooling, this is assured by having the eigenvalues of M (0) on the lefthand complex plane [110], as we displayed in Fig. 5.In the case of squeezing when the modulation depth, d, approaches one, we observed that the numerical convergence of the Floquet method becomes quite demanding even before the system crosses the instability.Seemingly, the integration-free computation also poses a potential limitation.Namely, as explained in Appendix C, it relies on the diagonalization of M (0) , which is granted whenever this matrix is not defective.Therefore, exactly at EPs, it is not diagonalizable, with the closest alternative being the Jordan normal form [102]. Nevertheless, in practice, one can get arbitrarily near to EP while the technique works flawlessly.As a matter of fact, this is the case in Fig. 4 where the horizontal dashed lines mark the first and second EPs, |µ EP,1 | and |µ EP,2 |.Finally, the strong optomechanical-coupling regime though it is experimentally challenging [123] should not pose a numerical difficulty as neither of these techniques is perturbative in nature.

VI. CONCLUSION
This study involves a theoretical exploration of closedloop phase control within an elementary optomechanical plaquette operating in the quantum regime.It highlights the cooling and squeezing of the mechanical resonator selected according to the loop phase.We provide a physical insight into the crucial role of non-Hermiticity in proximity to the EPs, linking it to the behavior of complex eigenvalue loci influenced by the intermechanical resonator coupling.Another advantageous aspect is the presence of a parameter space containing an exceptional surface, ensuring the robustness of EP singularities under parameter variations.The laser detuning corresponds to the parameter causing deviation from the exceptional surface and facilitates frequency-dependent squeezing, a characteristic of practical significance.There are a number of directions in which this research can be further expanded.Foremost, the triangle plaquette can be enriched in the number of resonators and network topology, like being extended to a two-dimensional lattice with various synthetic flux patterns [124].Alternatively, the underlying U (1) Abelian symmetry can be promoted to higher continuous symmetry groups or non-Abelian gauge potentials [78].As opposed to these fundamental aspects, on the practical side, a timely pursuit is to assess the potential of the elementary optomechanical plaquette for quantum sensing and information processing purposes.Both of these require an in-depth analysis of quantum noise, specifically the role of thermal fluctuations, imprecision noise (output shot noise) and quantum back-action noise (radiation-pressure force noise) [30].These can be investigated in both cooling and squeezing settings to comprehend the quantum advantages in real-world sens-ing applications.
We set ℓ max = 1 since |V (2) ii | ≪ V (1) ii , then maximum squeezing corresponding to minimum variance is min(V ii (t)) = V In this appendix, we consider the case where there is no modulation so that Eq. ( 16) reduces to iωI + M (0) R(ω) = − N(ω) .(C1) As M (0) has no ω dependence, for the integration-free computation of steady-state variances, our aim is to extract the explicit frequency dependence of T(ω) when the above equation is inverted to This is achieved by first applying a similarity transformation that diagonalizes M (0) as where U is simply composed of the eigenvectors of M (0) , with the corresponding eigenvalues being {λ i }.After this diagonalization, matrix inverse is trivially performed, which is followed by undoing the similarity transformation, leading to the frequency dependence of the entries of T matrix as With this explicit form, we can perform the integration for the steady-state covariance matrix using the residue calculus, which yields the entries Here, sgn(•) is the sign function, and Re{•} extracts the real part.As an alternative expression, in the literature, the so-called symmetrically ordered covariance matrix is widely used, which has the entries ⟨ Ri (t) Rj (t) + Rj (t) Ri (t)⟩ , (C7) and it satisfies the following Ricatti/Lyapunov matrix differential equation [64,74,75,103,113] By setting to zero the left-hand side of Eq. (C8) for the steady state, and performing the same similarity transformation on M (0) , we obtain an even simpler expression as where

FIG. 4 .
FIG. 4. Steady-state mean phonon number of the second mechanical resonator with respect to (a) |µ| where loop-phase is taken as ϕ ℓ ≃ 0.49π which gives the minimum steady-state phonon number, marked with a cross on the right subplot and (b) with respect to |µ| and ϕ ℓ .The horizontal dashed lines mark the first and second EPs, |µEP,1| and |µEP,2|.
FIG. 10.(a) Steady-state mean phonon number of the second mechanical resonator at |µEP,1| ≃ 52.5 (γ1 + γ2) in absence of modulation (d = 0), (b) minimum variance with respect to the vacuum value (in dB) of the second mechanical resonator at |µEP,1| ≃ 51.36 (γ1 + γ2) with modulation depth d = 0.4 as a function of ∆ normalized to mechanical frequency for different mechanical coupling constants; reservoir temperatures are taken as T = 18.1 K and T = 1.9 K, respectively.