Nonreciprocal Flow of Fluctuations, Populations and Correlations Between Doubly Coupled Bosonic Modes

1. Introduction

It has been demonstrated in recent years that unconventional properties of non-Hermitian parity-time $\left(\mathcal{PT}\right)-$ symmetric systems such as exceptional points and nonreciprocity [1,2,3,4,5,6] are not restricted to the time reversal systems [7,8,9], but can be obtained in quantum systems with a Hermitian Hamiltonian by interfering linear (excitation-preserving) and nonlinear two-mode interactions [10,11,12,13,14,15,16]. Subsequently, similar approaches have been studied including simultaneous application of linear and dissipative interactions [17,18] or creation of two coupling processes by engineering collective atomic systems [19,20]. The simultaneous existence of two coupling processes between quantum systems has been experimentally realized among optomechanical cavities [21,22] and in superconducting circuits [23,24].

To date, most of the treatments of the dynamics of doubly coupled modes have either been concerned with entanglement dynamics [25,26,27], quantum sensing [28,29] and topological phases [30,31,32,33,34]. Here, we consider a system composed of two Gaussian bosonic modes and concentrate on their fluctuation and correlation properties for features indicative of the simultaneous presence of the linear and nonlinear couplings processes. The two mode system provides the simplest example of the effects generated by the double coupling. In particular, we investigate how the coupling influences on the fluctuations, populations and correlation properties of the modes. The treatment includes the dissipation of the modes to local thermal as well as to squeezed reservoirs which, as we will see can have a significant effect on the fluctuation and correlation properties of the modes.

Modes interacting with local thermal or squeezed reservoirs will be sensitive to the number of photons and in the case of squeezed reservoirs also to the correlations present in the reservoirs and in general will evolve in a phase-sensitive fashion to a stationary state that will reflect such correlations. A squeezed reservoir is characterized by the mean photon number n and by the phase dependent two-photon correlations m, which range from zero to $\left|m\right|=\sqrt{n(n+1)}$. Values of the correlation in the range $\sqrt{n(n+1)}\ge \left|m\right|>n$ correspond a quantum squeezed reservoir while values of the correlation in the range $0<\left|m\right|\le n$ correspond to a classically squeezed reservoir, and $m=0,n\ne 0$ correspond to a thermal reservoir [35,36]. Presence of the correlations m leads to a reduction of the fluctuations in one quadrature of the modes.

It is well known that when only linear coupling is applied between two modes, it cannot create any correlations between noisy Gaussian modes [37]. On the other hand when the nonlinear coupling is applied, it can create two-photon correlations between the modes which are necessary for entanglement but the correlated modes are left mutually incoherent [38,39,40,41,42,43]. Moreover the modes are strongly amplified in both populations and fluctuations such that the modes are found in a highly fluctuating thermal state [44,45,46].

In this paper, we demonstrate that when both linear and nonlinear interactions are simultaneously applied the fluctuation and correlation properties of the modes are significantly different. Analytical expressions are obtained for steady state variances of the quadrature components of the mode operators and for correlation functions which show that relative to mutual coupling strengths of the two interactions an exceptional point emerges which is characteristic of non-Hermitian systems. To put it another way, our results demonstrate that in the Hermitian quantum system composed of simultaneously linearly and nonlinearly coupled modes one can construct non-Hermitian dynamics which can lead to nonreciprocal (one directional) influences of the modes on each other. It also creates an exceptional point which separates two distinctive parameter regimes, an exponential amplification regime and an oscillatory regime. In these two regimes the fluctuation properties of the modes are found to be differently altered even to the point of turning the fluctuations from thermal to quantum squeezed fluctuations. Under suitable conditions, the double coupling can establish nonresiprocal transfer of population and correlations between the modes.

The plan of this paper is as follows. In Section 2, we introduce the Hamiltonian of doubly coupled bosonic modes and give a description of independent external reservoirs to which the modes are dumped. Then, we derive the quantum Langevin equations for the quadrature components of the field operators and discuss their nonresiprocal coupling properties. In Section 3, we give general solutions for the time-dependent quadrature operators which exhibit the presence of an exceptional point separating two parameter regimes, exponential amplification and oscillatory regimes. In Section 4, we specialize to the exponential amplification regime and apply the solutions to obtain analytic stationary expressions of physical quantities of interest as the variances of the quadrature operators, populations of the modes, single- and two-mode correlation functions, and to investigate their dependence on the noise properties of thermal and squeezed reservoirs. Analytic expressions of these physical quantities in the oscillatory regime are presented and extensively discussed in Section 5. Differences and similarities of the results obtained in those two regimes are also examined. Finally, a brief discussion of the results and conclusions are given in Section 6.

2. Doubly Coupled Bosonic Modes

The system we study consists of two frequency degenerate radiation modes described by bosonic creation (annihilation) operators, $a^{†}\left(a\right)$ and $b^{†}\left(b\right)$, respectively. The modes are directly coupled to each other through presence of two different types of interaction processes, the linear optical photon exchange (beamsplitter type) and nonlinear two-photon (two-mode squeezing) interaction processes. In addition, both modes are damped with the rate $\kappa$ by coupling to local (independent) reservoirs, which will be considered as thermal or squeezed reservoirs, see Figure 1.

The interaction between the modes is determined by the Hamiltonian $(\hslash =1)$

$$\begin{array}{c}\hfill H=\lambda \left(a^{†}b{e}^{-i({\varphi }_a-{\varphi }_b)}+a{b}^{†}{e}^{i({\varphi }_a-{\varphi }_b)}\right)+g\left(ab{e}^{i({\varphi }_a+{\varphi }_b)}+a^{†}{b}^{†}{e}^{-i({\varphi }_a+{\varphi }_b)}\right),\end{array}$$

(1)

where $\lambda$ and g are real parameters determining the strength, respectively, the linear and nonlinear coupling between the modes, and ${\varphi }_i(i=a,b)$ is the phase of the i-the mode. The simultaneous action of the two interactions will couple the modes and will give the systematic evolution of the modes due to the mutual interaction. Without loss of generality we will assume that the phase of the mode a is fixed at ${\varphi }_a=0$, while the phase ${\varphi }_b\equiv \varphi$ of the mode b can be varied that $\varphi$ will play the role of the relative phase between squeezed reservoirs.

Using the input-output formalism, we readily find that the mode operators obey the quantum Langevin equations

$$\begin{array}{c}\hfill {\displaystyle \frac{du}{dt}}=-i\left[u,H\right]-\kappa u-\sqrt{2\kappa }{u}^{\mathrm{in}},\end{array}$$

(2)

where $u=a,b$. The first term on the right-hand side of Equation (2) gives the evolution due to a coupling between the modes, the second term determines damping of the mode with the rate $\kappa$, which we will assume the same for both modes, resulting from the interaction of the modes with local reservoirs. This interaction also gives rise to fluctuations in the system, these being represented by the input noise operators $u^{\mathrm{in}}$.

It is more convenient to determine dynamics of the modes in terms of Hermitian operators representing the amplitudes of two quadrature phase components of the fields

$$\begin{array}{cc}\hfill X_a& ={\displaystyle \frac{1}{\sqrt{2}}}\left(a+a^{†}\right),Y_a={\displaystyle \frac{1}{\sqrt{2}i}}\left(a-a^{†}\right),\hfill \\ \hfill X_b& ={\displaystyle \frac{1}{\sqrt{2}}}\left(b{e}^{i\varphi }+b^{†}{e}^{-i\varphi }\right),Y_b={\displaystyle \frac{1}{\sqrt{2}i}}\left(b{e}^{i\varphi }-b^{†}{e}^{-i\varphi }\right),\hfill \end{array}$$

(3)

with similar expressions for the input-noise operators. Under the double coupling between the modes the quadrature components obey the following equations of motion

$$\begin{array}{cc}\hfill {\dot{X}}_b& =-\kappa X_b-(g-\lambda )Y_a-\sqrt{2\kappa }{X}_b^{\mathrm{in}},\hfill \\ \hfill {\dot{Y}}_a& =-\kappa Y_a-(g+\lambda )X_b-\sqrt{2\kappa }{Y}_a^{\mathrm{in}},\hfill \\ \hfill {\dot{X}}_a& =-\kappa X_a-(g-\lambda )Y_b-\sqrt{2\kappa }{X}_a^{\mathrm{in}},\hfill \\ \hfill {\dot{Y}}_b& =-\kappa Y_b-(g+\lambda )X_a-\sqrt{2\kappa }{Y}_b^{\mathrm{in}},\hfill \end{array}$$

(4)

We first observe that the equations of motion for the quadratures $X_b$ and $Y_a$ are decoupled from the equations of motion for $X_a$ and $Y_b$. Further, inside each pair of equations, quadratures are coupled to each other in an asymmetric manner that the X quadratures are coupled to Y quadratures with strength $(g-\lambda )$, whereas Y quadratures are coupled to X with strength $(g+\lambda )$. For $g=\lambda$ we have complete coupling asymmetry that the evolution of the X quadratures is decoupled from the evolution of the Y quadratures. In this case there is a unidirectional coupling between the quadratures. Thus, we expect chiral propagation of fluctuations and correlations from Y quadratures to X quadratures of the other mode.

The dynamics of the quadrature operators is influenced by the quantum noise operators $X_i^{\mathrm{in}}$ and $Y_i^{\mathrm{in}}$ of the input modes to which the modes a and b are coupled. They obey Gaussian statistics and are delta correlated in time. We assume that the input noise operators represent two independent reservoirs which are in general in broadband squeezed vacuum states whose the properties are characterised by the correlation functions

$$\begin{array}{cc}\hfill \langle X_a^{\mathrm{in}}\left(t\right)X_a^{\mathrm{in}}\left(t^{\prime }\right)\rangle & =\left({\displaystyle \frac{1}{2}}+n+m\right)\delta (t-t^{\prime }),\hfill \\ \hfill \langle Y_a^{\mathrm{in}}\left(t\right)Y_a^{\mathrm{in}}\left(t^{\prime }\right)\rangle & =\left({\displaystyle \frac{1}{2}}+n-m\right)\delta (t-t^{\prime }),\hfill \\ \hfill \langle X_b^{\mathrm{in}}\left(t\right)X_b^{\mathrm{in}}\left(t^{\prime }\right)\rangle & =\left({\displaystyle \frac{1}{2}}+n+mcos2\varphi \right)\delta (t-t^{\prime }),\hfill \\ \hfill \langle Y_b^{\mathrm{in}}\left(t\right)Y_b^{\mathrm{in}}\left(t^{\prime }\right)\rangle & =\left({\displaystyle \frac{1}{2}}+n-mcos2\varphi \right)\delta (t-t^{\prime }),\hfill \end{array}$$

(5)

where n is the average number of thermal photons and m describes the degree of phase dependent two-photon correlations between photons contained in the reservoirs. The presence of the two-photon correlations leads to an asymmetric distribution of noise between the quadratures. We assume that the phase of the squeezed reservoir coupled to the mode a is fixed at zero, whereas the phase $\varphi$ of the squeezed reservoir coupled to the mode b can be varied. The parameters n and m are related to each other such that $m\le n$ corresponds to a classically squeezed field whereas $m\le \sqrt{n(n+1)}$ corresponds to a quantum squeezed field [35,36].

3. Dynamics of the Doubly Coupled Bosonic Modes

To proceed with the solution of the set of differential Equation (4), it is convenient to introduce the Laplace transforms of the quadratures

$$\begin{array}{c}\hfill X_i\left(p\right)\equiv {\int }_0^{\infty }{X}_i{e}^{-pt}dt,Y_i\left(p\right)\equiv {\int }_0^{\infty }{Y}_i{e}^{-pt}dt,\end{array}$$

(6)

whose the application converts the equations into a set of algebraic equations. The set of the algebraic equations can be readily solved for the transformed quadratures to give

$$\begin{array}{cc}\hfill X_i\left(p\right)& ={\displaystyle \frac{(p+\kappa )L_{ix}\left(p\right)-(g-\lambda )L_{jy}\left(p\right)}{(p-p_1)(p-p_2)}},\hfill \\ \hfill \\ \hfill Y_i\left(p\right)& ={\displaystyle \frac{(p+\kappa )L_{iy}\left(p\right)-(g+\lambda )L_{jx}\left(p\right)}{(p-p_1)(p-p_2)}},i\ne j=a,b,\hfill \end{array}$$

(7)

in which we have introduced the abbreviations

$$\begin{array}{c}\hfill L_{ix}\left(p\right)=\left[X_i\left(0\right)-\sqrt{2\kappa }{X}_i^{\mathrm{in}}\left(p\right)\right],L_{iy}\left(p\right)=\left[Y_i\left(0\right)-\sqrt{2\kappa }{Y}_i^{\mathrm{in}}\left(p\right)\right],\end{array}$$

(8)

and $p_1$ and $p_2$ are the roots of the quadratic equation

$$\begin{array}{c}\hfill p^2+2\kappa p+{\kappa }^2+{\lambda }^2-g^2.\end{array}$$

(9)

The roots can be easily computed, to give

$$\begin{array}{c}\hfill p_{1,2}=-\kappa \pm \sqrt{g^2-{\lambda }^2}.\end{array}$$

(10)

It is clearly seen from Equation (10) that the roots are strongly dependent on the relationship between the coupling constants g and $\lambda$. We distinguish two parameter regimes that depending on whether $\lambda >g$ or $\lambda <g$ the factor $\sqrt{g^2-{\lambda }^2}$ can be either the real or the complex parameter. A threshold at which the roots change character is for $g=\lambda$, which corresponds to the exceptional point. For $g>\lambda$ the roots are real

$$\begin{array}{c}\hfill p_1=-\kappa +\alpha ,p_2=-\kappa -\alpha ,\alpha =\sqrt{g^2-{\lambda }^2},\end{array}$$

(11)

while for $\lambda >g$ the roots are complex

$$\begin{array}{c}\hfill p_3=-\kappa +i\beta ,p_4=-\kappa -i\beta ,\beta =\sqrt{{\lambda }^2-g^2}.\end{array}$$

(12)

Thus there are two distinct regimes in which the solutions have different character. We will call the regime $g>\lambda$ an exponential amplification regime, and $\lambda >g$ an oscillatory regime.

For each of the regimes we will investigate properties of the steady state $(t\to \infty )$ variances of the quadrature operators, populations of the modes, $\langle a^{†}a\rangle$ and $\langle b^{†}b\rangle$, single-mode two-photon correlations, $\langle a^2\rangle$ and $\langle b^2\rangle$, which we will call as "local two-photon correlations", and two-mode two-photon correlations $\langle ab\rangle$,47,48] which we will call as "global two-photon correlations". Apart from the global two-photon correlations we will also consider the single-photon two-mode correlations $\langle a^{†}b\rangle$, which are known to carry information about coherence properties the modes [37]. This restricted class of the correlation functions considered here results from the fact that the dynamics of the modes considered are Gaussian. Note that the local two-photon correlations are responsible for squeezing of the single-mode fluctuations, whereas the global two-photon correlations are necessary for entanglement [49,50,51,52,53].

By inverse Laplace transformation and application of the Cauchy residue theorem we then have from Equation (7), with the help of Equation (11) that in the case of $g>\lambda$ the time evolution of the quadrature operators is given by

$$\begin{array}{cc}\hfill X_i\left(t\right)& ={\displaystyle \frac{1}{2}}\left\{[L_{ix}\left(p_1\right)-u{L}_{jy}\left(p_1\right)]e^{p_{1}t}+[L_{ix}\left(p_2\right)+u{L}_{jy}\left(p_2\right)]e^{p_{2}t}\right\},\hfill \\ \hfill Y_i\left(t\right)& ={\displaystyle \frac{1}{2u}}\left\{-[L_{jx}\left(p_1\right)-u{L}_{iy}\left(p_1\right)]e^{p_{1}t}+[L_{jx}\left(p_2\right)+u{L}_{iy}\left(p_2\right)]e^{p_{2}t}\right\},i\ne j=a,b.\hfill \end{array}$$

(13)

where $u=\sqrt{(g-\lambda )/(g+\lambda )}$.

Similarly, the inverse Laplace transformation of Equation (7) with the complex roots (12) yields to

$$\begin{array}{cc}\hfill X_i\left(t\right)& ={\displaystyle \frac{1}{2}}\left\{[L_{ix}\left(p_3\right)-iw{L}_{jy}\left(p_3\right)]e^{p_{3}t}+[L_{ix}\left(p_4\right)+iw{L}_{jy}\left(p_4\right)]e^{p_{4}t}\right\},\hfill \\ \hfill Y_i\left(t\right)& ={\displaystyle \frac{i}{2w}}\left\{[L_{jx}\left(p_3\right)-iw{L}_{iy}\left(p_3\right)]e^{p_{3}t}-[L_{jx}\left(p_4\right)+iw{L}_{iy}\left(p_4\right)]e^{p_{4}t}\right\},i\ne j=a,b.\hfill \end{array}$$

(14)

where $w=\sqrt{(\lambda -g)/(\lambda +g)}$. In this case the quadrature evolve in an oscillatory manner.

The dynamics of the quadrature operators is influenced by the quadratures $X_i^{\mathrm{in}}\left(p\right)$ and $Y_i^{\mathrm{in}}\left(p\right)$ of the input modes, whose the statistics are given by the statistics of the reservoirs to which the modes are coupled. The effect of the statistics of the reservoirs will be evident in the fluctuation and correlation properties of the modes.

4. Exponential Amplification Regime, $g>\lambda$.

We begin our discussion of the correlation and fluctuation properties of the modes by considering first the exponential regime $g>\lambda$.

Let us first consider the average number of photons of the modes, i.e., populations of the modes. Since

$$\begin{array}{cc}\hfill a& ={\displaystyle \frac{1}{\sqrt{2}}}\left(X_a+i{Y}_a\right),a^{†}={\displaystyle \frac{1}{\sqrt{2}}}\left(X_a-i{Y}_a\right),\hfill \\ \hfill b& ={\displaystyle \frac{1}{\sqrt{2}}}\left(X_b+i{Y}_b\right),b^{†}={\displaystyle \frac{1}{\sqrt{2}}}\left(X_b-i{Y}_b\right),\hfill \end{array}$$

(15)

and we have solutions for the time evolution of the quadrature operators, the populations can be evaluated from

$$\begin{array}{cc}\hfill \langle a^{†}a\rangle & ={\displaystyle \frac{1}{2}}\left(\langle X_a^2\rangle +\langle Y_a^2\rangle -1\right),\langle b^{†}b\rangle ={\displaystyle \frac{1}{2}}\left(\langle X_b^2\rangle +\langle Y_b^2\rangle -1\right),\hfill \end{array}$$

(16)

where $\langle X_i^2\rangle$ and $\langle Y_i^2\rangle$ are variances of the quadrature operators of the modes. Thus, evaluation of the populations of the modes requires calculations of the variances of the quadrature operators of the modes.

The Gaussian character of the system enables the quadrature variances to be readily calculated. With the use of Equations (13) and (5), we find that the variances of the quadrature components of the mutually doubly coupled modes are

$$\begin{array}{cc}\hfill \langle X_a^2\rangle & =\left({\displaystyle \frac{1}{2}}+n+m\right)\left(1+{\displaystyle \frac{1}{2}}{sinh}^2\psi \right)+{\displaystyle \frac{1}{2}}{\left(\overline{g}-\overline{\lambda }\right)}^2\left({\displaystyle \frac{1}{2}}+n-mcos2\varphi \right){cosh}^2\psi ,\hfill \\ \hfill \langle Y_a^2\rangle & =\left({\displaystyle \frac{1}{2}}+n-m\right)\left(1+{\displaystyle \frac{1}{2}}{sinh}^2\psi \right)+{\displaystyle \frac{1}{2}}{\left(\overline{g}+\overline{\lambda }\right)}^2\left({\displaystyle \frac{1}{2}}+n+mcos2\varphi \right){cosh}^2\psi ,\hfill \\ \hfill \langle X_b^2\rangle & =\left({\displaystyle \frac{1}{2}}+n+mcos2\varphi \right)\left(1+{\displaystyle \frac{1}{2}}{sinh}^2\psi \right)+{\displaystyle \frac{1}{2}}{\left(\overline{g}-\overline{\lambda }\right)}^2\left({\displaystyle \frac{1}{2}}+n-m\right){cosh}^2\psi ,\hfill \\ \hfill \langle Y_b^2\rangle & =\left({\displaystyle \frac{1}{2}}+n-mcos2\varphi \right)\left(1+{\displaystyle \frac{1}{2}}{sinh}^2\psi \right)+{\displaystyle \frac{1}{2}}{\left(\overline{g}+\overline{\lambda }\right)}^2\left({\displaystyle \frac{1}{2}}+n+m\right){cosh}^2\psi .\hfill \end{array}$$

(17)

where all system parameters have been normalised by the damping rate $\kappa$: $\overline{g}=g/\kappa$, $\overline{\lambda }=\lambda /\kappa$, $\overline{\alpha }=\alpha /\kappa$, and the angle $\psi$ is defined through the relation $\psi =arctanh\left(\overline{\alpha }\right)$. In what follows we will drop the bar on $\overline{g},\overline{\lambda }$, and $\overline{\alpha }$ with the understanding that we are dealing with dimensionless rescaled quantities.

Equation (17) shows that the double coupling between the modes enhances the variances indicating enhancement of the fluctuations of the modes. However, the variances are unequally enhanced that the variances of the X quadratures are less enhanced than the variances of the Y quadratures. This is easily understood since, according to Equation (4) there is an asymmetry in the coupling between the X and Y quadratures leading to unbalanced transfer of fluctuations between the quadratures.

Once the variances have been determined, it is only a matter of substitution of Equation (17) in Equation (16) to derive explicit expressions for the steady state populations. Thus, we obtain the following expressions for the populations

$$\begin{array}{cc}\hfill \langle a^{†}a\rangle & =n+\left[\left({\displaystyle \frac{1}{2}}+n\right)g^2+g\lambda mcos2\varphi \right]{cosh}^2\psi ,\hfill \end{array}$$

(18)

and

$$\begin{array}{cc}\hfill \langle b^{†}b\rangle & =n+\left[\left({\displaystyle \frac{1}{2}}+n\right)g^2+g\lambda m\right]{cosh}^2\psi .\hfill \end{array}$$

(19)

The first terms on the right-hand sides of Equations (18) and (19) represent the population of the modes in the absence of the direct coupling between the modes. The second terms on the right-hand sides are due to the coupling between the modes. Of particular interest to us is the dependence of the populations on the term $g\lambda$, which accounts for double coupling between the modes.

It is quite evident from Equations (18) and (19) that the populations are amplified by the mutual couplings and the major role in the amplification plays the nonlinear coupling process. Since $tanh\psi$ is limited to one, it means that the amplifier is operating below threshold for the stable steady-state solutions [44,45,54].

Note that the manner in which the populations are amplified by the couplings is different for $\langle a^{†}a\rangle$ and $\langle b^{†}b\rangle$, and the populations are very dependent on the nature of the reservoirs. Perhaps the most significant is the insensitivity of the population of the mode b to the phase of the squeezed reservoir to which the mode is coupled. As it is seen from Equation (18) the phase sensitivity is transferred to the mode a. Thus, a variation of the phase $\varphi$ of the squeezed reservoir coupled to the mode b will modify the population of the other mode.

The squeezed reservoirs provide the mechanism for controlled amplification of the mode a that the population of the mode can vary with the phase in such a way that it could be possible to cease the amplification process. Specifically, at the exceptional point, when $\lambda =g$ with $m=\sqrt{n(n+1)}$ and in the strong squeezing limit $n\gg 1$, Equation (18) takes a simple form

$$\begin{array}{cc}\hfill \langle a^{†}a\rangle & =n+\left({\displaystyle \frac{1}{2}}+n\right)\left(1+cos2\varphi \right)g^2.\hfill \end{array}$$

(20)

This result directly shows that for strongly squeezed reservoirs the population is extremely phase sensitive at the exceptional point, and the choice of phase $\varphi =\pi /2$ leads to ceasing of the amplification process.

It is worth noting that the phenomenon of ceasing of the amplification process of mode a, derived using $m=\sqrt{n(n+1)}$, is unique to quantum squeezing that it is not possible to cease the the amplification process when $m=n$, i.e., the modes interact with classically squeezed fields of the reservoirs. In addition, this phase-dependent amplification of the mode gives us a new control over nonlinear effects in the system of coupled modes.

Figure 2 shows the variation of the population of the mode a with the strength of the linear coupling $\lambda$ for two limiting cases of the phase $\varphi =0$ and $\varphi =\pi /2$ and two different values of the nonlinear coupling strength g. Note that the limit $g\to 1$ corresponds to approaching the instability threshold in the nonlinear amplifying interaction [44,45,54]. For g away from the instability threshold the population varies with the phase and the linear coupling such that for $\varphi =0$ it is further amplified and de-amplified for $\varphi =\pi /2$. For g close to the instability threshold, the population exhibits a weak dependence on the phase, rapidly degrades with increasing strength of the linear interactions and returns to its initial value as $\lambda$ approaches the exceptional point $\lambda =g$.

4.1. Two-Photon Correlations Inside the Modes

We now turn to the calculation of the single mode two-photon correlation functions, which are known to give rise to squeezing of the fluctuations of the quadrature operators. To evaluate the correlation functions we use Equation (15) and find

$$\begin{array}{cc}\hfill 〈aa〉& ={\displaystyle \frac{1}{2}}\left\{〈X_a^2〉-〈Y_a^2〉+2i〈X_a{Y}_a〉+1\right\},\hfill \\ \hfill 〈bb〉& ={\displaystyle \frac{1}{2}}\left\{〈X_b^2〉-〈Y_b^2〉+2i〈X_b{Y}_b〉+1\right\}.\hfill \end{array}$$

(21)

It is seen that the two-photon correlation functions are not only related to the variances but also to the correlation functions $\langle X_i{Y}_i\rangle$, which are readily evaluated with the use of Equation (13) and are given by

$$\begin{array}{cc}\hfill 〈X_a\left(t\right)Y_a\left(t\right)〉& ={\displaystyle \frac{1}{2}}i+{\displaystyle \frac{1}{2}}msin2\varphi {sinh}^2\psi ,\hfill \end{array}$$

(22)

and

$$\begin{array}{cc}\hfill 〈X_b\left(t\right)Y_b\left(t\right)〉& ={\displaystyle \frac{1}{2}}i+msin2\varphi \left(1+{\displaystyle \frac{1}{2}}{sinh}^2\psi \right).\hfill \end{array}$$

(23)

Then, the two-photon correlation functions obtained from Equation (21) are

$$\begin{array}{cc}\hfill 〈aa〉& =m-\left[\left({\displaystyle \frac{1}{2}}+n\right)g\lambda +msin\varphi e^{-i(\varphi +{\displaystyle \frac{\pi }{2}})}{g}^2+mcos\varphi e^{i\varphi }{\lambda }^2\right]{cosh}^2\psi ,\hfill \\ \hfill 〈bb〉& =m{e}^{2i\varphi }-\left[\left({\displaystyle \frac{1}{2}}+n\right)g\lambda +msin\varphi e^{i(\varphi -{\displaystyle \frac{\pi }{2}})}{g}^2+mcos\varphi e^{i\varphi }{\lambda }^2\right]{cosh}^2\psi .\hfill \end{array}$$

(24)

The contribution of the first terms on the right-hand side of Equation (24) obviously correspond to correlations transferred to the modes in the absence of the inter-mode couplings. The second terms correspond to creation of the correlations inside the modes by the inter-mode couplings. It is seen that the effect of the presence of the double coupling between the modes is to create phase-insensitive and phase sensitive contributions to the correlations.

The phase insensitive contributions corresponds to those created when the reservoirs are in thermal states $(m=0)$, and then the correlations are

$$\begin{array}{cc}\hfill 〈aa〉& =〈bb〉=-\left({\displaystyle \frac{1}{2}}+n\right)g\lambda {cosh}^2\psi ,\hfill \end{array}$$

(25)

This shows that nonzero two-photon correlations arise solely from the presence of the double coupling between the modes. Consequently, uncoupled modes being in thermal states are turned by the double coupling to squeezed states. This is clearly seen when one considers variances of the quadrature components, Equation (17), which in the case of $m=0$ reduce to

$$\begin{array}{cc}\hfill \langle X_a^2\rangle & =\langle X_b^2\rangle =\left({\displaystyle \frac{1}{2}}+n\right)\left[1+g(g-\lambda ){cosh}^2\psi \right],\hfill \\ \hfill \langle Y_a^2\rangle & =\langle Y_b^2\rangle =\left({\displaystyle \frac{1}{2}}+n\right)\left[1+g(g+\lambda ){cosh}^2\psi \right].\hfill \end{array}$$

(26)

Evidently, $\langle X_i^2\rangle <\langle Y_i^2\rangle ,(i=a,b)$ when $\lambda \ne 0$ indicating that the modes are in squeezed states with the variances of X quadratures reduced at the expense of increased variances of Y quadratures. Further inside into the variances reveals that $\langle X_i^2\rangle$ decreases as $\lambda$ is increasing and a maximum squeezing, i.e., the minimum value of $\langle X_i^2\rangle$ is reached for $\lambda =g$, i.e., at the exceptional point, in which case $\langle X_i^2\rangle =(1/2+n)$. Since $\langle X_i^2\rangle \ge 1/2$, it follows that the maximum reduction of the variances (fluctuations) corresponds to maximally classically squeezed fields.

Conditions for squeezed fluctuations of the modes can be more conveniently examined by investigating the so-called degrees of squeezing defined as

$$\begin{array}{cc}\hfill {\eta }_{aa}& ={\displaystyle \frac{|〈aa〉|-〈a^{†}a〉}{〈a^{†}a〉}},{\eta }_{bb}={\displaystyle \frac{|〈bb〉|-〈b^{†}b〉}{〈b^{†}b〉}}.\hfill \end{array}$$

(27)

Negative values of ${\eta }_{ii}$, with the limiting value $-1$, indicate classically squeezed field, whereas positive values of ${\eta }_{ii}$ indicate a quantum squeezed field, and the largest the positive value the greater quantum squeezing.

When we specify to thermal reservoirs interacting with the modes, we get

$$\begin{array}{cc}\hfill {\eta }_{aa}& ={\eta }_{bb}=-{\displaystyle \frac{n+\left(\frac{1}{2}+n\right)g\left(g-\lambda \right){cosh}^2\psi }{n+\left(\frac{1}{2}+n\right)g^2{cosh}^2\psi }},\hfill \end{array}$$

(28)

which is always negative. Thus, fluctuations of both modes cannot be reduced below the vacuum limit of ${\displaystyle \frac{1}{2}}$. Therefore, both modes display classically squeezed fluctuations.

Further, we would like to point out that Equation (25) provides an interesting example of an unexpected result that interactions applied between the modes create correlations inside the modes. Such effect is absent when only one type of the interaction, linear or nonlinear, is applied between the modes.

The phase sensitive contributions to the correlations (24) are due to correlations m existing in squeezed reservoirs. A close look at the correlations, Equation (24), reveals that the two-photon correlations inside the modes are created by the linear and nonlinear coupling processes from the correlations existing in the reservoirs under completely different conditions. It can be seen in Equation (24) that for $\varphi =0$ the phase sensitive contributions are due to the linear coupling process $\left(\lambda \right)$, whereas for $\varphi =\pi /2$ the contribution to the correlations is done by the nonlinear process $\left(g\right)$.

Let us examine how the correlations behave at the exceptional point, when $\lambda =g$. In this case, we obtain from Equation (24) that

$$\begin{array}{cc}\hfill 〈aa〉& =m-\left({\displaystyle \frac{1}{2}}+n+mcos2\varphi \right)g^2,\hfill \\ \hfill 〈bb〉& =m{e}^{2i\varphi }-\left({\displaystyle \frac{1}{2}}+n+m\right)g^2,\hfill \end{array}$$

(29)

from which it is apparent that the correlations contained in the mode a are sensitive to phase of the reservoir to which mode b is coupled. Remembering that in the large quantum squeezing limit, $n-m\to -{\displaystyle \frac{1}{2}}$, we see that the proces of enhancement of the correlations inside the mode a can be significantly reduced and ultimately ceased if one choses $\varphi =\pi /2$. Thus, for strongly squeezed reservoirs the same dramatic phase dependence exhibited by the population of the mode a also appears in the single-mode two-photon correlations.

4.2. Correlations Between the Modes

Now we turn our attention to the two-mode correlations which could exist between the modes. There could be one-photon and two-photon correlations existing between the two modes which are determined by the correlation functions $\langle a^{†}b\rangle$ and $\langle ab\rangle$, respectively. From Equation (15), we find that the correlations can be evaluated by calculating the correlations between quadratures of different modes

$$\begin{array}{cc}\hfill 〈a^{†}b〉& ={\displaystyle \frac{1}{2}}\left[〈X_a{X}_b〉+〈Y_a{Y}_b〉+i\left(〈X_a{Y}_b〉-〈Y_a{X}_b〉\right)\right],\hfill \end{array}$$

(30)

and

$$\begin{array}{cc}\hfill 〈ab〉& ={\displaystyle \frac{1}{2}}\left[〈X_a{X}_b〉-〈Y_a{Y}_b〉+i\left(〈X_a{Y}_b〉+〈Y_a{X}_b〉\right)\right].\hfill \end{array}$$

(31)

The correlation functions can be evaluated using Equation (13) with the help of Equation (5), and their explicit steady-state expressions are given in the Appendix A. Hence, by inserting Equations (A1) and (A2) into Equations (30) and (31), we obtain

$$\begin{array}{cc}\hfill \langle a^{†}b\rangle & =-mgsin\varphi e^{i\varphi }{cosh}^2\psi ,\hfill \end{array}$$

(32)

and

$$\begin{array}{cc}\hfill 〈ab〉& =-i\left[\left({\displaystyle \frac{1}{2}}+n\right)g+m\lambda cos\varphi e^{i\varphi }\right]{cosh}^2\psi .\hfill \end{array}$$

(33)

From Equation (32), we see that the interaction of the modes with squeezed reservoirs $(m\ne 0)$ is essential for generating the first-order correlations between the modes that it is not possible to have $\langle a^{†}b\rangle \ne 0$ when the modes interact with thermal reservoirs. In addition, the correlations are specifically dependent on the relative phase of the squeezed reservoirs. A comparison of Equations (32) and (33) shows that the difference that the phase $\varphi$ has on the one-photon and two-photon correlations. Namely, the choice of phase $\varphi =0$ results in $\langle a^{†}b\rangle =0$ with simultaneous enhancement of the two-photon correlations $\langle ab\rangle$. Conversely, the choice of phase $\varphi =\pi /2$ results in $\langle a^{†}b\rangle$ maximal and $\langle ab\rangle$ minimal. Thus, the first-order coherence between the modes can be destroyed or preserved depending on the relative phase of the reservoirs.

In addition, Equations (32) and (33), reveal that the one-photon correlations between the modes are created solely by the nonlinear coupling process from the two-photon correlations existing in the squeezed reservoirs. However, the two-photon correlations between the modes are created by both linear and nonlinear coupling processes. The effect of the nonlinear coupling is to create a phase insensitive contribution from the reservoirs thermal noise $({\displaystyle \frac{1}{2}}+n)$, whereas the phase sensitive contribution is created by the linear coupling process from the two-photon correlations existing in the squeezed reservoirs. Note that $\langle ab\rangle$ becomes phase insensitive when phase $\varphi =\pi /2$.

As $\langle a^{†}b\rangle$ can be nonzero, it follows that the modes can be coherent to the first-order, and we can determine the degree of the first-order coherence by considering the quantity ${\gamma }_{ab}=\left|〈a^{†}b〉\right|/\sqrt{〈a^{†}a〉〈b^{†}b〉}$, which lies between 0 for mutually incoherent and 1, for mutually perfectly coherent fields. For the case considered here, the explicit analytical form of ${\gamma }_{ab}$ can be found using Equations (32), (18) and (19). With the choice of phase $\varphi =\pi /2$, at which $|\langle a^{†}b\rangle |$ is maximal, we obtain

$$\begin{array}{cc}\hfill {\gamma }_{ab}& ={\displaystyle \frac{mg{cosh}^2\psi }{\sqrt{{\left[n+\left(\frac{1}{2}+n\right)g^2{cosh}^2\psi \right]}^2-{\left(mg\lambda {cosh}^2\psi \right)}^2}}}.\hfill \end{array}$$

(34)

Figure 3 shows the coherence function ${\gamma }_{ab}$ as a function of n and as a function of $\lambda /\kappa$ for $g=0.99\kappa$. It is seen that the degree of coherence increases with n such that for $n>1$ this is close to unity, indicating that the modes are almost perfectly coherent.

In order to study two-photon correlations between the modes it is useful to define the degree of two-photon correlations between the modes by introducing the quantity ${\eta }_{ab}=\left|〈ab〉\right|/\sqrt{〈a^{†}a〉〈b^{†}b〉}$, which lies between 0 and ∞. However, there is a threshold value of ${\eta }_{ab}$ which distinguishes between classically and quantum two-photon correlated fields. Namely, for classically correlated fields ${\eta }_{ab}$ can be no larger than 1, and ${\eta }_{ab}>1$ can be achieved only for quantum correlated fields. It should be pointed out that a value of ${\eta }_{ab}$ above one indicates that the modes may be entangled [43,49,50,51,52,53].

When the modes interact with squeezed reservoirs, we find that the choices of phase $\varphi =0$ and $\varphi =\pi /2$ require significantly different conditions for ${\eta }_{ab}$ to be larger than one. With the help of Equations (33), (18) and (19), we find that

$$\begin{array}{cc}\hfill {\eta }_{ab}& =1+{\displaystyle \frac{(1-g)\left[\left(\frac{1}{2}+n\right)g+m\lambda \right]{cosh}^2\psi -n}{n+g\left[\left(\frac{1}{2}+n\right)g+m\lambda \right]{cosh}^2\psi }},\hfill \end{array}$$

(35)

when $\varphi =0$, and

$$\begin{array}{cc}\hfill {\eta }_{ab}& ={\displaystyle \frac{\left(\frac{1}{2}+n\right)g{cosh}^2\psi }{\sqrt{{\left[n+\left(\frac{1}{2}+n\right)g^2{cosh}^2\psi \right]}^2-{\left(mg\lambda {cosh}^2\psi \right)}^2}}}.\hfill \end{array}$$

(36)

when $\varphi =\pi /2$. Since in the exponential amplification regime, the condition for stable steady-state solutions [37,44,45,54] requires $g<1$, we see from Equation (35) that in the case of the ordinary vacuum $(m=n=0)$ it ensures that always ${\eta }_{ab}$ is larger than one. However, for the thermal or squeezed vacuum, where $n\ne 0$ the first term in the numerator of Equation (35) may not be large enough to enforce ${\eta }_{ab}>1$.

To search for limits imposed by n on ${\eta }_{ab}$ being larger than one, we plot in Figure 4 the degree of two-photon correlations ${\eta }_{ab}$ as a function of $\lambda$ and as a function of n for $g=.99$ and $m=\sqrt{n(n+1)}$. The largest value of ${\eta }_{ab}$ occurs for $n=0$ at which ${\eta }_{ab}>1$, degrades with an increasing n and turns into values smaller than one for $n\approx 0.5$. The linear interaction has the effect to decrease further the correlations.

A comparison of Figure 3 and Figure 4 immediately shows the difference that the number of thermal photons has on the degree of coherence ${\gamma }_{ab}$ and on the degree of two-photon correlations ${\eta }_{ab}$. As it is clearly seen, ${\eta }_{ab}$ is large in the region of $n<0.5$ in which ${\gamma }_{ab}$ is small, while for the values of n where ${\gamma }_{ab}$ is large, ${\eta }_{ab}$ is smaller than one and even tends to zero. For $n\gg 1$ the degree of coherence can, in principle, approach one $({\gamma }_{ab}=1)$, the perfect coherence. Thus, unlike the degree of the first-order correlation which approaches maximal values for large n, the degree of the two-photon correlations has maximal values for small n. In addition, the perfect coherence ${\gamma }_{ab}=1$ can be achieved only if the reservoirs are in in quantum squeezed states, $m=\sqrt{n(n+1)}$. For maximally classically squeezed reservoirs, $m=n$, the degree of coherence is significantly less than one.

5. Oscillatory Regime $\lambda >g$

We now turn to consider the case of $\lambda >g$, in which the quadrature components exhibit oscillatory behaviour and, as in the previous section, we study the fluctuation properties of the modes, their populations, single and two-mode correlations. The starting point of our calculations this time is Equation (14), the solutions for the time evolution of the quadrature components for $\lambda >g$ and we will first consider the steady-state properties of the variances of the quadrature components.

5.1. Fluctuation Properties of the Modes

Following the same procedure as described in the preceding section, if we evaluate variances of the quadrature components according to Equation (14), we arrive to the following expressions

$$\begin{array}{cc}\hfill \langle X_a^2\rangle & =\left({\displaystyle \frac{1}{2}}+n+m-m{cos}^2\varphi {sin}^2\chi \right)-\left({\displaystyle \frac{1}{2}}+n-mcos2\varphi \right)g(\lambda -g){cos}^2\chi ,\hfill \\ \hfill \langle Y_a^2\rangle & =\left({\displaystyle \frac{1}{2}}+n-m+m{cos}^2\varphi {sin}^2\chi \right)+\left({\displaystyle \frac{1}{2}}+n+mcos2\varphi \right)g(\lambda +g){cos}^2\chi ,\hfill \\ \hfill \langle X_b^2\rangle & =\left({\displaystyle \frac{1}{2}}+n+mcos2\varphi -m{cos}^2\varphi {sin}^2\chi \right)-\left({\displaystyle \frac{1}{2}}+n-m\right)g(\lambda -g){cos}^2\chi ,\hfill \\ \hfill \langle Y_b^2\rangle & =\left({\displaystyle \frac{1}{2}}+n-mcos2\varphi +m{cos}^2\varphi {sin}^2\chi \right)+\left({\displaystyle \frac{1}{2}}+n+m\right)g(\lambda +g){cos}^2\chi ,\hfill \end{array}$$

(37)

where $\chi =arctan\left(\beta \right)$. These results are in marked contrast to those found in the exponential region $g>\lambda$ given by Equation (17). In particular, we can see immediately that outside the exceptional point $(g\ne \lambda )$ the variances of the X quadratures are reduced while the variances of the Y quadratures are enhanced by the double coupling.

The first important fact we can derive from Equation (37) is that even when the modes interact with thermal reservoirs, the variances of the X quadratures can be reduced below their vacuum value of $1/2$, which indicates that the modes can be found in quantum squeezed states. To show it more explicitly, we set $m=0$ in Equation (37) and find that

$$\begin{array}{cc}\hfill \langle X_a^2\rangle & =\langle X_b^2\rangle =\left({\displaystyle \frac{1}{2}}+n\right)\left[1-g(\lambda -g){cos}^2\chi \right],\hfill \\ \hfill \langle Y_a^2\rangle & =\langle Y_b^2\rangle =\left({\displaystyle \frac{1}{2}}+n\right)\left[1+g(\lambda +g){cos}^2\chi \right].\hfill \end{array}$$

(38)

Clearly, the fluctuations in $X_i$ quadratures can be reduced below $1/2$ if $g\ne \lambda$. This implies that outside the exceptional point the fields of both modes can display quantum squeezed fluctuations.

In Figure 5 we plot the variances of the X quadratures as a function of the scaled nonlinear coupling strength $g/\kappa$ for $n=0$ and several different values of $\lambda /\kappa$. The figure shows clearly that the variances are reduced below $1/2$ for all values of $g<\lambda$ indicating that the modes are in quantum squeezed states and the maximum squeezing occurs for g in the vicinity of the exceptional point. It is seen that approaching the exceptional point cause cessation of quantum squeezing. The largest reduction of the variances is achieved for $\lambda \gg \kappa$, in which case ${\Delta }^2{X}_i\approx 0.25$, so that we may speak of $50\%$ reduction of the fluctuations below the quantum level. Again, we point out that this result is in sharp contrast with the exponential case of $g>\lambda$, Equation (26), where the variances were reduced to only the vacuum level.

We may readily adapt the expressions (37) for the variances to calculate populations of the modes. On making use of variances (37) in Equation (16), we obtain

$$\begin{array}{cc}\hfill 〈a^{†}a〉& =n+\left[\left({\displaystyle \frac{1}{2}}+n\right)g^2+mcos2\varphi g\lambda \right]{cos}^2\chi ,\hfill \end{array}$$

(39)

and

$$\begin{array}{cc}\hfill 〈b^{†}b〉& =n+\left[\left({\displaystyle \frac{1}{2}}+n\right)g^2+mg\lambda \right]{cos}^2\chi .\hfill \end{array}$$

(40)

Apart from the appearance of $cos\chi$ in place of $cosh\psi$, Equations (39) and (40) are formally identical with the corresponding results found in Section 4 for the exponential case $g>\lambda$. As for the case $g>\lambda$, only the population of the mode a depends on the phase such that for $\varphi =0$ the modes are equally populated and by varying the phase, it is possible to vary the population such that effect of amplification of the population ultimately ceases to occur.

5.2. Two-Photon Correlations Inside the Modes

We now consider the single-mode two-photon correlations. Using the variances (37) in Equation (21), we find

$$\begin{array}{cc}\hfill 〈aa〉& =m-mcos\varphi e^{i\varphi }{sin}^2\chi -\left[\left({\displaystyle \frac{1}{2}}+n\right)g\lambda +m{g}^{2}cos2\varphi \right]{cos}^2\chi ,\hfill \\ \hfill 〈bb〉& =m{e}^{2i\varphi }-mcos\varphi e^{i\varphi }{sin}^2\chi -\left[\left({\displaystyle \frac{1}{2}}+n\right)g\lambda +m{g}^2\right]{cos}^2\chi .\hfill \end{array}$$

(41)

Further inside into the correlation properties of the modes is gained by considering the degrees of the correlations. If we restrict ourselves to the case where the doubly coupled modes interact with local thermal reservoirs, we get that the degrees of the correlations are

$$\begin{array}{cc}\hfill {\eta }_{aa}& ={\eta }_{bb}={\displaystyle \frac{\left(\frac{1}{2}+n\right)g\left(\lambda -g\right){cos}^2\chi -n}{n+\left(\frac{1}{2}+n\right)g^2{cos}^2\chi }}.\hfill \end{array}$$

(42)

Apparently, the degrees of correlations can have positive values, which is quite different from the exponential amplification case of $g\ge \lambda$, where ${\eta }_{aa}$ and ${\eta }_{bb}$ were always negative. It is particularly interesting to note that in the case of the ordinary vacuum state $(n=0)$ of the reservoirs, ${\eta }_{aa}$ and ${\eta }_{bb}$ are always positive. Therefore, the fluctuations of the modes can be reduced below the vacuum limit so that the modes can display quantum squeezed fluctuations. The effect of thermal photons n is to diminish the correlations and eventually to turn ${\eta }_{aa}$ and ${\eta }_{bb}$ into negative values. In other words, thermal fluctuations can turn the fluctuations of the modes from quantum to classically squeezed fluctuations.

The feature of the correlations just described are easily seen in Figure 6, which shows ${\eta }_{aa}$ and ${\eta }_{bb}$, as given in Equation (42) as a function of $g/\kappa$ and n for $\lambda =5\kappa$. The largest value of ${\eta }_{aa}$ and ${\eta }_{bb}$ occurs for small g and n. It is also seen that with an increasing g the threshold for ${\eta }_{aa}$ and ${\eta }_{bb}$ to be positive shifts towards larger n.

To continue the analysis of the two-photon correlations, we now assume that the modes interact with squeezed reservoirs. In this case, the degrees of the two-photon correlations depend on the phase $\varphi$ such that the choice of the phase $\varphi =0$ results in the degrees equal for both modes

$$\begin{array}{cc}\hfill {\eta }_{aa}& ={\eta }_{bb}={\displaystyle \frac{\left\{\left|m-g\left[\left(\frac{1}{2}+n\right)\lambda +mg\right]\right|-g\left[\left(\frac{1}{2}+n\right)g+m\lambda \right]\right\}{cos}^2\chi -n}{n+g\left[\left(\frac{1}{2}+n\right)g+m\lambda \right]{cos}^2\chi }}.\hfill \end{array}$$

(43)

The absolute values appearing in the numerator of this equation involve a difference between the two-photon correlations existing in the squeezed reservoirs and correlations generated by the double coupling between the modes.

In Figure 7, we plot ${\eta }_{aa}$ and ${\eta }_{bb}$ as a function of $g/\kappa$ and n for $\lambda =5\kappa$ and quantum squeezed reservoirs with $m=\sqrt{n(n+1)}$. We can distinguish two separate ranges of the parameters for which the degrees are positive. For the first, occurring for $g<0.2\kappa$ the quantum correlations are those transferred to the modes due to their interaction with the reservoirs. The correlations decrease with an increasing g and ultimately cease at $g=.2\kappa$. As g increases further quantum correlations appear again showing that the double coupling can generate quantum correlations inside the modes independent of the state of the reservoirs.

The choice of phase $\varphi =\pi /2$ leads to widely different behaviour of the degrees of the correlations of the modes that

$$\begin{array}{cc}\hfill {\eta }_{aa}& ={\displaystyle \frac{\left|m-g\left[\left(\frac{1}{2}+n\right)\lambda -mg\right]{cos}^2\chi \right|-g\left[\left(\frac{1}{2}+n\right)g-m\lambda \right]{cos}^2\chi -n}{n+g\left[\left(\frac{1}{2}+n\right)g-m\lambda \right]{cos}^2\chi }},\hfill \end{array}$$

(44)

and

$$\begin{array}{cc}\hfill {\eta }_{bb}& ={\displaystyle \frac{m-n+g\left(\frac{1}{2}+n-m\right)(\lambda -g){cos}^2\chi }{n+g\left[\left(\frac{1}{2}+n\right)g+m\lambda \right]{cos}^2\chi }}.\hfill \end{array}$$

(45)

Since $\lambda >g$, it is clearly seen from the structure of the numerator in Equation (45) that ${\eta }_{bb}$ is always positive. Hence, the mode b will exhibit quantum squeezed fluctuations even if the reservoirs are in classically squeezed states. This is illustrated in Figure 8, where we plot ${\eta }_{bb}$ as a function of g for $\lambda =5,n=1$ and several different values of $m\le n$. For $m=0.5n$, ${\eta }_{bb}$ is negative over the entire range of g indicating classically squeezed correlations present in the mode. For $m=0.75n$ quantum correlations appear in the restricted range of large g. For $m=0.9n$ quantum correlations appear in less restricted range of g and for $m=n$ quantum correlations occur over the entire range of g.

5.3. Correlations Between the Modes

Finally, we consider properties of the two-mode correlations which in the case of $\lambda >g$ follow from Equations (30) and (31) with the correlations functions involving quadrature operators of different modes given by Equations (A3) and (A4). Thus, we find that

$$\begin{array}{cc}\hfill \langle a^{†}b\rangle & =-mgsin\varphi e^{i\varphi }{cos}^2\chi ,\hfill \end{array}$$

(46)

and

$$\begin{array}{cc}\hfill 〈ab〉& =-i\left[\left({\displaystyle \frac{1}{2}}+n\right)g+m\lambda cos\varphi e^{i\varphi }\right]{cos}^2\chi .\hfill \end{array}$$

(47)

From the form of the correlations we can see that the one-photon correlation function is different from zero only for choices of phase between 0 and $\pi /2$ and attains its maximal values for phase $\varphi =\pi /4$, at which the linear and nonlinear processes equally contribute to the two-photon correlation function. This result is in marked contrast to that found in the case of $g>\lambda$, where $\langle a^{†}b\rangle$ attained a maximal value at $\varphi =\pi /2$. At that phase the two-photon correlation function was independent of the correlations m.

A better inside into the correlations is obtained by considering the degrees of the correlations. Consider first the degree of the first-order coherence ${\gamma }_{ab}$. Using Equation (46) with the populations given in Equations (39) and (40), we readily find that for $\varphi =\pi /2$ the degree of coherence ${\gamma }_{ab}$ is

$$\begin{array}{cc}\hfill {\gamma }_{ab}& ={\displaystyle \frac{mg{cos}^2\chi }{\sqrt{{\left[n+\left(\frac{1}{2}+n\right)g^2{cos}^2\chi \right]}^2-{\left(mg\lambda {cos}^2\chi \right)}^2}}}.\hfill \end{array}$$

(48)

The expression is plotted in Figure 9 for $\lambda =\kappa$ and $m=\sqrt{n(n+1)}$. The degree of the first-order coherence increases with g and attains maximal values which is no larger than $0.6$. This result is significantly different from that obtained in the exponential amplification regime $g>\lambda$, where the degree of the first-order coherence can be as large as 1.

To examine properties of the two-mode two-photon correlations, we use Equation (47), which together with the populations given by Equations (39) and (40) gives

$$\begin{array}{cc}\hfill {\eta }_{ab}& =1+{\displaystyle \frac{(1-g)\left[\left(\frac{1}{2}+n\right)g+m\lambda \right]{cos}^2\chi -n}{n+g\left[\left(\frac{1}{2}+n\right)g+m\lambda \right]{cos}^2\chi }},\hfill \end{array}$$

(49)

for the choice of phase $\varphi =0$, and

$$\begin{array}{cc}\hfill {\eta }_{ab}& ={\displaystyle \frac{\left(\frac{1}{2}+n\right)g{cos}^2\chi }{\sqrt{{\left[n+\left(\frac{1}{2}+n\right)g^2{cos}^2\chi \right]}^2-{\left(mg\lambda {cos}^2\chi \right)}^2}}},\hfill \end{array}$$

(50)

for the choice of phase $\varphi =\pi /2$.

It is clear from Equation (49) that for $\varphi =0$ the minimum requirement for ${\eta }_{ab}$ to be larger than one is that g should be smaller than one. Since in the oscillatory regime there are no limits imposed on g, we have that in the case of the ordinary vacuum $(m=n=0)$ the requirement $g<1$ is necessary and sufficient. However, for the thermal or squeezed vacuum it is necessary but not sufficient. For $n\ne 0$ the first term in the numerator of Equation (49) may not be large enough to enforce ${\eta }_{ab}>1$. To search for limits imposed by n on ${\eta }_{ab}$ being larger than one, we plot in Figure 10 the degree of two-photon correlations ${\eta }_{ab}$ as a function of g and as a function of n for $\lambda =0.8$ and $m=n$. The degree of the correlations is initially for $g=0$ equal to one, but increases to values greater than 1 with an increasing g. For $g>0$ and small n the correlations attain maximal values indicating that at the parameter regime strong quantum two-photon correlations exist between the modes.

6. Discussion and Conclusions

In this paper we have studied the fluctuation and correlation properties of two bosonic modes mutually coupled through a linear photon exchange and nonlinear parametric type processes. This Hermitian quantum system is known to exhibit non-Hermitian dynamics, which can lead to nonreciprocal (one directional) influences of the modes on each other [16,26]. In addition, it results in the appearance of an exceptional point separating two parameter regimes, exponential amplification and oscillatory regimes. The fluctuation and correlation properties of the modes have been investigated by evaluating the stationary expressions for the variances of the quadrature components of the field operators, populations of the modes, and single- and two-mode correlations. We have assumed that apart from the presence of the double coupling, the modes are in contact (interact) with local thermal or squeezed reservoirs, and have shown that the creation of the correlations by the nonreciprocal coupling is strongly influenced by the noise properties of the reservoirs. The differences in the fluctuations and correlations for these two regimes were studied in details. In particular, for the exponential amplification regime the nonreciprocal coupling tends to convert thermal fluctuations of independent modes being in contact with thermal reservoirs into classically squeezed fluctuations. In the oscillatory regime, however, thermal fluctuations of the modes can be turned by the nonreciprocal coupling into quantum squeezed fluctuations. Apart from the fluctuations, the nonreciprocal coupling can have significant effect on the amplification of the population of the modes. We have shown that when the modes interact with strongly squeezed reservoirs the amplification of the population of one of the modes can be controlled by varying phase of the reservoir interacting with the other mode, and even can be completely ceased. We have also discuss conditions under which the modes could be coherent and simultaneously entangled. We have found that the coherence and entanglement exclude each other that strongly coherent modes exhibit classical rather than quantum two-photon correlations, and vice versa, strongly two-photon correlated modes are behaving as being mutually incoherent.