This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Robust photon transmission in nonlinear parity-time-symmetric cavities

Ling-Pu Gong Department of Physics, Chongqing University, Chongqing 401330, China    Xing-Sen Chen Department of Physics, Chongqing University, Chongqing 401330, China    Yin Tan Department of Physics, Chongqing University, Chongqing 401330, China    Rui Zhang Department of Physics, Chongqing University, Chongqing 401330, China    Yu-Yu Zhang yuyuzh@cqu.edu.cn Department of Physics, Chongqing University, Chongqing 401330, China
Abstract

We explore the photon transfer in the nonlinear parity-time-symmetry system of two coupled cavities, which contains nonlinear gain and loss dependent on the intracavity photons. Analytical solution to the steady state gives a saturated gain, which satisfy the parity-time symmetry automatically. The eigen-frequency self-adapts the nonlinear saturated gain to reach the maximum efficiency in the steady state. We find that the saturated gain in the weak coupling regime does not match the loss in the steady state, exhibiting an appearance of a spontaneous symmetry-breaking. The photon transmission efficiency in the parity-time-symmetric regime is robust against the variation of the coupling strength, which improves the results of the conventional methods by tuning the frequency or the coupling strength to maintain optimal efficiency. Our scheme provides an experimental platform for realizing the robust photon transfer in cavities with nonlinear parity-time symmetry.

I Introduction

Considerable progress has been made in parity-time (PT)- symmetric optics to explore intriguing properties makris08 ; klaiman ; kottos ; suchkov ; lv15 , which has been investigated both theoretically bender ; klaiman08 ; jin10 ; hassan15 and experimentally ruter ; peng ; sun14 . PT symmetric system, despite being non-Hermiticity, has been found novel phenomenons, such as PT symmetry breaking ruter ; peng , and nonreciprocal reflectionless transmission feng2013 . Even through the Hermiticity of the quantum observables was never in doubt, a wide class of non-Hermitian Hamiltonians with PT symmetry have motivated discussions on several fronts in physics, including cavity optomechanics jing ; fan ; peng12 , quantum field theories bender04 , PT-symmetric lattice longhi ; zhen15 , and open quantum system rotter ; dembowski . A variety of optical structures provide a alternative platform for testing various theoretical proposals on non-Hermitian PT-symmetric quantum mechanic chang14 .

PT-symmetric systems are invariant under under simultaneous parity-flip and time-flip operations ganainy07 ; chong11 . In optical systems, PT symmetry can be established by incorporating gain and loss in coupled resonators. A linear PT-symmetric system with balanced gain and loss can exhibit a real eigenvalue spectrum and present unusual properties jin17 ; guo09 . If the gain-loss contrast exceeds a certain threshold, the PT symmetry can be spontaneously broken and the spectrum is no longer entirely real, exhibiting an exceptional point. When the gain induces large nonlinear saturation under high pumping and no longer matches the loss, the linear PT symmetry of the coupled cavities breaks down. Theoretically, the nonlinear gain saturation can causes the system to reach a steady state that still contains the PT symmetry characteristics fan17 ; hassan ; ge . Can such nonlinear gain saturation effects be harnessed to enhance the photon transmission in coupled cavities? Photons transmission, such reflection and transmission as a consequence of the nonlinear PT symmetry motivates us to achieve an optimal photon transfer scheme.

We focus on how PT symmetry relates to the nonlinear gain and loss system leads to a robust photon transfer in two coupled cavities. Different from the linear PT-symmetric system with fixed value of gain, we consider a nonlinear PT-symmetric system with a nonlinear gain depending on intracavity photons. An analysis based on quantum Langevin equations is presented to study the steady-state dynamics. It is observed that the eigenfrequency self-adjusts to the nonlinear gain saturation, and the saturated gain in the steady state automatically satisfies PT symmetry. As a consequence of the PT symmetry with nonlinear gain saturation, the photons transfer efficiency obtained is robust independent on the coupling strength, exhibiting an improvement over a conventional transfer scheme, in which the input light frequency is required to adjusted as the coupling strength to maintain optimal efficiency.

II Conventional scheme

To explore the photon transfer scheme in coupled cavities, we consider an input cavity coupling to an output cavity. The Hamiltonian is described by

H=ω1a1a1+ω2a2a2+κ(a1a2+H.c.),H=\omega_{1}a_{1}^{\dagger}a_{1}+\omega_{2}a_{2}^{\dagger}a_{2}+\kappa(a_{1}^{\dagger}a_{2}+H.c.), (1)

where aia_{i}^{\dagger} (aia_{i}) is the creation (annihilation) operator of the cavity ii with the frequency ωi\omega_{i}, κ\kappa is the coupling strength between two cavities.

Refer to caption
Figure 1: (a)Conventional scheme. A driving at a frequency ω\omega interacts with the input cavity at a rate γ1\gamma_{1}. The input cavity couples to the output cavity with the coupling strength κ\kappa, which outputs photons at a rate γ2\gamma_{2}. The intrinsic loss rate of two cavities are γ0\gamma_{0}. (b) PT-symmetric scheme. We consider a nonlinear gain in the input cavity with high pumping theoretically.

Fig. 1(a) shows a conventional scheme of photon transfer in two coupled cavities. A pump light field aina_{in} of a frequency ω\omega transfers photons to the input cavity at a rate γ1\gamma_{1}, and the output cavity out-puts photons at a rate γ2\gamma_{2} due to the leakage of photons. The intrinsic loss rate of each cavity is γ0\gamma_{0}.

The dynamics of the photon transfer process undergoes dissipation due to leakage of photons. It can be described by the quantum Langevin equations for the Heisenberg operators a1a_{1} and a2a_{2},

a˙1(t)\displaystyle\dot{a}_{1}(t) =\displaystyle= (iω1γ0γ1)a1(t)iκa2(t)+f1(t),\displaystyle(-i\omega_{1}-\gamma_{0}-\gamma_{1})a_{1}(t)-i\kappa a_{2}(t)+f_{1}(t), (2)
a˙2(t)\displaystyle\dot{a}_{2}(t) =\displaystyle= (iω2γ0γ2)a2(t)iκa1(t)+f2(t),\displaystyle(-i\omega_{2}-\gamma_{0}-\gamma_{2})a_{2}(t)-i\kappa a_{1}(t)+f_{2}(t), (3)

where f1(t)=2γ1ain(t)f_{1}(t)=\sqrt{2\gamma_{1}}a_{in}(t) is the quantum Langevin force originating from the input light field ain(t)a_{in}(t). For simplicity, the force term f2(t)f_{2}(t) is neglected.

To solve above equations in frequency space, we employ a transforming operator as a(ω)=+eiωta(t)𝑑ta(\omega)=\int_{-\infty}^{+\infty}e^{i\omega t}a(t)dt of a(t)a(t). The input-output relations are given by

a1out(ω)\displaystyle a_{1out}(\omega) =\displaystyle= 2γ1a1(ω)ain(ω)=R(ω)ain(ω),\displaystyle\sqrt{2\gamma_{1}}a_{1}(\omega)-a_{in}(\omega)=R(\omega)a_{in}(\omega), (4)
a2out(ω)\displaystyle a_{2out}(\omega) =\displaystyle= 2γ2a2(ω)=T(ω)ain(ω),\displaystyle\sqrt{2\gamma_{2}}a_{2}(\omega)=T(\omega)a_{in}(\omega), (5)

where a1outa_{1out} and a2outa_{2out} are the reflected and output operators of cavities, respectively. For an input light field at frequency ω\omega, the transmission function T(ω)T(\omega) of the output cavity is obtained as

T(ω)=iκ2γ22γ1κ2+[iΔ1+γ0+γ1][iΔ2+γ0+γ2],T(\omega)=-\frac{i\kappa\sqrt{2\gamma_{2}}\sqrt{2\gamma_{1}}}{\kappa^{2}+[i\Delta_{1}+\gamma_{0}+\gamma_{1}][i\Delta_{2}+\gamma_{0}+\gamma_{2}]}, (6)

and the response function R(ω)R(\omega) of the input cavity for the reflection is

R(ω)=2γ1iΔ1+γ0+γ1+κ2/[iΔ2+γ0+γ2]1,R(\omega)=\frac{2\gamma_{1}}{i\Delta_{1}+\gamma_{0}+\gamma_{1}+\kappa^{2}/[i\Delta_{2}+\gamma_{0}+\gamma_{2}]}-1, (7)

where Δ1(2)=ω+ω1(2)\Delta_{1(2)}=-\omega+\omega_{1(2)}. The corresponding photon transfer efficiency from output cavity is given by

η\displaystyle\eta =\displaystyle= a2outa2outainain=|T(ω)|2,\displaystyle\frac{\langle a_{2out}^{\dagger}a_{2out}\rangle}{\langle a_{in}^{\dagger}a_{in}\rangle}=|T(\omega)|^{2},

which, for a symmetric case (γ1=γ2=γ0\gamma_{1}=\gamma_{2}=\gamma_{0} and ω1=ω2=ω0\omega_{1}=\omega_{2}=\omega_{0}), becomes

η=|2κγκ2+[iΔ+2γ0]2|2.\eta=|\frac{2\kappa\gamma}{\kappa^{2}+[i\Delta+2\gamma_{0}]^{2}}|^{2}. (9)

with the detuning frequency Δ=ω+ω0\Delta=-\omega+\omega_{0}. In the conventional photon transport process, it involves reflections R(ω)R(\omega), which plays a role in the photon transport from the input cavity.

Fig. 2(a) shows the transfer efficiency of photons with a resonance frequency ω=ω0\omega=\omega_{0}. It behaves non-monotonically dependent on the coupling strength κ\kappa. The type of the conventional photon transfer is not robust against coupling strength in the operating conditions. To maximize the transfer efficiency η\eta, it is required ω=ω0±κ24γ02\omega=\omega_{0}\pm\sqrt{\kappa^{2}-4\gamma_{0}^{2}} for κ2γ0\kappa\geqslant 2\gamma_{0} by minimizing the module of the denominator in Eq.(9). One can tune the input frequency according to the coupling strength κ\kappa to achieve the optimal value of η\eta.

Refer to caption
Figure 2: (a)Photon transfer efficiency as a function of the scaled coupling strength κ/κc\kappa/\kappa_{c} for a conventional scheme (red dashed line) with the input frequency ω=ω0\omega=\omega_{0} and a PT-symmetry scheme (blue solid line). (b)Steady-state mode frequency ω\omega and the saturated gain gg as a function of the scaled coupling strength κ/κc\kappa/\kappa_{c}. The parameters are g0=3f0g_{0}=3f_{0}, γ0=0.005ω\gamma_{0}=0.005\omega, and κc=2γ0\kappa_{c}=2\gamma_{0}.

III Nonlinear PT-symmetric scheme

Different from the conventional scheme, the PT-symmetric scheme circumvent photon reflections. The linear PT-symmetric system of two coupled resonators with balanced gain and loss was previously investigated, in which a gain can be induced by pumping erbium ions in the gain resonatorchang14 ; jing . Here, we consider the lasing of the gain medium provides the nonlinear gain theoretically, which depends on the intracavity photons fan17 ; hassan . Fig. 1(b) illustrates a nonlinear saturable gain element into the input cavity for amplification of photons generation and a nonlinear loss in the output cavity. The nonlinear gain rate gg and loss rate γ\gamma depend on the intracavity average photons a1(2)a1(2)=|α1(2)|2\langle a_{1(2)}^{\dagger}a_{1(2)}\rangle=|\alpha_{1(2)}|^{2}, which are expressed as

g(α1)=γ0+g01+|α1|2,γ(α2)=γ0+f01+|α2|2,g(\alpha_{1})=-\gamma_{0}+\frac{g_{0}}{1+|\alpha_{1}|^{2}},\quad\gamma(\alpha_{2})=\gamma_{0}+\frac{f_{0}}{1+|\alpha_{2}|^{2}}, (10)

with the unsaturated gain and loss rates g0g_{0} and f0f_{0}, respectively.

The Heisenberg operators a1a_{1} and a2a_{2} obey

a˙1(t)=[iω0+g(α1)]a1(t)iκa2(t),\dot{a}_{1}(t)=[-i\omega_{0}+g(\alpha_{1})]a_{1}(t)-i\kappa a_{2}(t), (11)
a˙2(t)=[iω0γ(α2)]a2(t)iκa1(t).\dot{a}_{2}(t)=[-i\omega_{0}-\gamma(\alpha_{2})]a_{2}(t)-i\kappa a_{1}(t). (12)

Since the nonlinear gain gg depends on the mean photons in the input cavity, we explore the stable state with a saturable gain by solving the above equations analytically.

PT-symmetric regime.–A enhancement of the nonlinear gain and loss generates large steady-state amplitudes in two cavities. We use the mean-field description for the operators αi=ai=αi0eiωt\alpha_{i}=\langle a_{i}\rangle=\alpha_{i0}e^{-i\omega t} lv16 ; lv , where αi0\alpha_{i0} is the steady-state amplitude of the cavity mode. According to Eqs.(11) and (12), the steady-state amplitudes is given by the following equations

i(ωω0)α10=(γ0+gs)α10iκα20,-i(\omega-\omega_{0})\alpha_{10}=(-\gamma_{0}+g_{s})\alpha_{10}-i\kappa\alpha_{20}, (13)
i(ωω0)α20=(γ0fs)α20iκα10,-i(\omega-\omega_{0})\alpha_{20}=(-\gamma_{0}-f_{s})\alpha_{20}-i\kappa\alpha_{10}, (14)

where the saturated gain and loss are gs=g0/(1+|α10|2)g_{s}=g_{0}/(1+|\alpha_{10}|^{2}), and fs=f0/(1+|α20|2)f_{s}=f_{0}/(1+|\alpha_{20}|^{2}), respectively. The above equations suggest the relation α02=ρeiϕα01\alpha_{02}=\rho e^{i\phi}\alpha_{01} with a phase shift ϕ\phi and the modal ratio ρRe+\rho\in{Re}^{+}. It leads to the equation for the eigen-frequencies

(ωω0)2+i(ωω0)(2γ0fsgs)\displaystyle(\omega-\omega_{0})^{2}+i(\omega-\omega_{0})(2\gamma_{0}-f_{s}-g_{s})
+(γ0+gs)(γ0fs)κ2=0.\displaystyle+(-\gamma_{0}+g_{s})(-\gamma_{0}-f_{s})-\kappa^{2}=0. (15)

To obtain a real ω\omega, it is required gs/2γ0fs/2γ0=1g_{s}/2\gamma_{0}-f_{s}/2\gamma_{0}=1. Then it is reasonable to give gs=2γ0cosh2ηg_{s}=2\gamma_{0}\cosh^{2}\eta and fs=2γ0sinh2ηf_{s}=2\gamma_{0}\sinh^{2}\eta with a positive real quantity η\eta. It leads to the relations for the mean photons

|α10|2\displaystyle|\alpha_{10}|^{2} =\displaystyle= g02γ0cosh2η1,\displaystyle\frac{g_{0}}{2\gamma_{0}\cosh^{2}\eta}-1, (16)
|α20|2\displaystyle|\alpha_{20}|^{2} =\displaystyle= f02γ0sinh2η1.\displaystyle\frac{f_{0}}{2\gamma_{0}\sinh^{2}\eta}-1. (17)

Then the eigenfrequencies in Eq. (III) reduces to (ωω0)2=κ2γ02cosh2(2η)(\omega-\omega_{0})^{2}=\kappa^{2}-\gamma_{0}^{2}\cosh^{2}(2\eta). By substituting into Eqs. (13) and (14), one obtains ρ=±1\rho=\pm 1 and tanhη=f0/g0\tanh\eta=\sqrt{f_{0}/g_{0}}.

The saturated gain and loss are obtained as

gs=2γ0g0g0f0,fs=2γ0f0g0f0.\displaystyle g_{s}=2\gamma_{0}\frac{g_{0}}{g_{0}-f_{0}},f_{s}=2\gamma_{0}\frac{f_{0}}{g_{0}-f_{0}}. (18)

Therefore, the stable state is in the PT-symmetric phase with the balanced gain and loss

g𝙿𝚃(α1)=γ𝙿𝚃(α2)=γ0g0+f0g0f0.g_{\mathtt{PT}}(\alpha_{1})=\gamma_{\mathtt{PT}}(\alpha_{2})=\gamma_{0}\frac{g_{0}+f_{0}}{g_{0}-f_{0}}. (19)

The eigenfrequencies are obtained as

(ωω0)2=κ2γ02(g0+f0g0f0)2,(\omega-\omega_{0})^{2}=\kappa^{2}-\gamma_{0}^{2}(\frac{g_{0}+f_{0}}{g_{0}-f_{0}})^{2}, (20)

which is real for κ>kc=γ0(g0+f0)/(g0f0)\kappa>k_{c}=\gamma_{0}(g_{0}+f_{0})/(g_{0}-f_{0}). By setting gs=3γ0g_{s}=3\gamma_{0} and fs=γ0f_{s}=\gamma_{0} with g0=3f0g_{0}=3f_{0}, the eigenfrequencies reduce into

ω=ω0±κ24γ02,\omega=\omega_{0}\pm\sqrt{\kappa^{2}-4\gamma_{0}^{2}}, (21)

with the PT-symmetry g𝙿𝚃=γ𝙿𝚃=2γ0g_{\mathtt{PT}}=\gamma_{\mathtt{PT}}=2\gamma_{0}. The eigenfrequencies are consistent with the optimal frequency to maximize the transfer efficiency in Eq. (9) in the conventional scheme. In contrast to the conventional scheme with frequency tuning, the advantage of the PT-symmetric scheme lies in the self-selected eigenfrequencies without any active tunings.

Refer to caption
Figure 3: Mean photons for each cavity in the PT-symmetry regime a1a1\langle a_{1}^{\dagger}a_{1}\rangle (blue dashed line) and a2a2\langle a_{2}^{\dagger}a_{2}\rangle (red solid line) obtained numerically as a function of time tt for different coupling strength (a)κ/κc=2\kappa/\kappa_{c}=2 and (b) κ/κc=4\kappa/\kappa_{c}=4. The initial state is choose as α1(0)=1\alpha_{1}(0)=1 and α2(0)=0.1eiπ/2\alpha_{2}(0)=0.1e^{i\pi/2}. The parameters are κc=1.1γ0\kappa_{c}=1.1\gamma_{0}, g0=0.1ωg_{0}=0.1\omega, f0=0.005ωf_{0}=0.005\omega, γ0=0.0113ω\gamma_{0}=0.0113\omega and ω0=1\omega_{0}=1.

Fig. 2 (b) show the saturated gain and the self-adjusting eigen-frequency ω\omega in the steady state. We observe a bifurcation in the real part of the frequency at the critical coupling strength κc\kappa_{c}, exhibiting a PT symmetry-breaking transition. In the strong coupling regime κ>kc\kappa>k_{c}, the system have two modes with eigenfrequencies ω=ω0±κ24γ02\omega=\omega_{0}\pm\sqrt{\kappa^{2}-4\gamma_{0}^{2}}. The corresponding saturated gain in Eq. (19) is a constant independent on the coupling strength, and can balance out the loss, g𝙿𝚃=γ𝙿𝚃g_{\mathtt{PT}}=\gamma_{\mathtt{PT}}. It demonstrate that the system posses the PT symmetry in the strong coupling regime κ>kc\kappa>k_{c}.

The steady state in the PT-symmetric phase is given by

(α10α20)=g0f02γ2γ(1eiϕ)\left(\begin{array}[]{c}\alpha_{10}\\ \alpha_{20}\end{array}\right)=\sqrt{\frac{g_{0}-f_{0}-2\gamma}{2\gamma}}\left(\begin{array}[]{c}1\\ e^{i\phi}\end{array}\right) (22)

with sinϕ=±γ0/κ\sin\phi=\pm\gamma_{0}/\kappa. One can easily obtain the mean photons in the steady state for the input and output cavities, which have the same value,

a2a2a1a1=|α20|2|α10|2=1.\frac{\langle a_{2}^{\dagger}a_{2}\rangle}{\langle a_{1}^{\dagger}a_{1}\rangle}=\frac{|\alpha_{20}|^{2}}{|\alpha_{10}|^{2}}=1. (23)

Fig. 3 shows the time evolution of the mean photons obtained by solving Eqs.(11) and (12) numerically for different coupling strength κ\kappa. It is observed that in the steady state a2a2\langle a_{2}^{\dagger}a_{2}\rangle equals to a1a1\langle a_{1}^{\dagger}a_{1}\rangle, which is consistent with the analytical results in Eq. (23). In the steady state, the mean photons a1(2)a1(2)\langle a_{1(2)}^{\dagger}a_{1(2)}\rangle for κ/κc=2\kappa/\kappa_{c}=2 is the same as that for κ/κc=4\kappa/\kappa_{c}=4, which is independent on the coupling strength.

With the input and output relations a^in=2gsa^1,a^2out=2fsa^2\hat{a}_{in}=\sqrt{2g_{s}}\hat{a}_{1},\hat{a}_{2out}=\sqrt{2f_{s}}\hat{a}_{2}, the photon transfer efficiency is obtained as

η=a^2outa^2outa^ina^in=fsgs,\eta=\frac{\langle\hat{a}_{2out}^{\dagger}\hat{a}_{2out}\rangle}{\langle\hat{a}_{in}^{\dagger}\hat{a}_{in}\rangle}=\frac{f_{s}}{g_{s}}, (24)

where gsg_{s} and fsf_{s} are the saturated gain and loss in Eq. (18). So the efficiency is independent on the coupling strength κ\kappa in the PT-symmetric regime.

Fig. 2(a) shows the transfer efficiency in the strong coupling regime κ>kc\kappa>k_{c}. The efficiency η\eta is robust against the variation of the coupling strength κ\kappa. And η\eta is larger than that obtained by the conventional transfer scheme, which varies dependent on the coupling strength for a fixed input frequency. Since the eigen-frequency in the non-linear PT-symmetric scheme can self-adjust to the variation of the the coupling strength for an optimal transfer efficiency. It exhibits an improvement over the conventional scheme due to the robust efficiency induced by the saturated gain and self-adjusting frequency in the steady state.

PT-broken regime.– In the PT-broken regime for a weak coupling strength κ<κc\kappa<\kappa_{c}, the real eigenfrequency is given by ω=ω0\omega=\omega_{0} from Eq. (21). The mean-field approximation for the stationary solutions is assumed to be αi=ai=αi0\alpha_{i}=\langle a_{i}\rangle=\alpha_{i0}. The equations of motion in the steady state can be given in the interacting representation

0=(γ0+gs)α10iκα20,0=(-\gamma_{0}+g_{s})\alpha_{10}-i\kappa\alpha_{20}, (25)
0=(γ0fs)α20iκα10.0=(-\gamma_{0}-f_{s})\alpha_{20}-i\kappa\alpha_{10}. (26)

By solving the above equations, one obtains the saturated gain

gs=γ0+κ2γ0+fs.g_{s}=\gamma_{0}+\frac{\kappa^{2}}{\gamma_{0}+f_{s}}. (27)

The corresponding gain and loss in the stable state are obtained by setting fs=γ0f_{s}=\gamma_{0}

g𝚋=γ0+gs=κ22γ0,γ𝚋=γ0+fs=2γ0,g_{\mathtt{b}}=-\gamma_{0}+g_{s}=\frac{\kappa^{2}}{2\gamma_{0}},\gamma_{\mathtt{b}}=\gamma_{0}+f_{s}=2\gamma_{0}, (28)

Obviously, it exhibits the PT symmetry broken, because the saturated gain does not match the loss g𝚋γ𝚋g_{\mathtt{b}}\neq\gamma_{\mathtt{b}}. Our results are consistent with the saturated gain in wireless power transfer system with magnetic resonators fan17 .

The ratio of mean photons in the input and output cavity is obtained analytically as

a2a2a1a1=|α20|2|α10|2=κ24γ02,\frac{\langle a_{2}^{\dagger}a_{2}\rangle}{\langle a_{1}^{\dagger}a_{1}\rangle}=\frac{|\alpha_{20}|^{2}}{|\alpha_{10}|^{2}}=\frac{\kappa^{2}}{4\gamma_{0}^{2}}, (29)

where α20/α10=iκ/(2γ0)\alpha_{20}/\alpha_{10}=-i\kappa/(2\gamma_{0}) is obtained from Eq.(25). Obviously, the ratio depends on the coupling strength κ\kappa, which is different from that in Eq.(23) in the PT-symmetry regime. Since the ration is smaller than 11, giving κ2/4γ021\kappa^{2}/4\gamma_{0}^{2}\leq 1, it leads to the PT-broken regime κ<κc\kappa<\kappa_{c} with the critical value κc=2γ0\kappa_{c}=2\gamma_{0}.

One obtains the photon transfer efficiency

η=fsa2a2gsa1a1=κ22(κ2+2γ02).\eta=\frac{f_{s}\langle a_{2}^{\dagger}a_{2}\rangle}{g_{s}\langle a_{1}^{\dagger}a_{1}\rangle}=\frac{\kappa^{2}}{2(\kappa^{2}+2\gamma_{0}^{2})}. (30)

Fig. 2(a) shows the transfer efficiency in the PT-broken regime for κ<κc\kappa<\kappa_{c}. η\eta increases as the coupling strength κ\kappa increases. Fig. 2(b) shows that the saturated gain in Eq. (28) is proportional to κ2\kappa^{2}. It is different the photon transfer efficiency and the saturated gain in the PT-symmetric regime for κ>κc\kappa>\kappa_{c}, which is robust independent on the coupling strength.

IV Conclusion

We have studied theoretically a nonlinear PT-symmetric scheme for photon transfer in two coupled cavities by introducing a nonlinear gain and loss. The gain induces nonlinearity dependent on the intracavity photons, which is different from the linear PT-symmetric system with the balanced gain and loss. We obtain the analytical solution to the steady state, which agree well with the numerical ones. The nonlinear saturated gain and loss in the steady state satisfy the PT symmetry automatically. And the eigen-frequency self-adjusts to the nonlinear saturated gain, and exhibits a bifurcation at the critical coupling strength. Consequence, we obtain the robust photon transmission efficiency independent on the coupling strength, exhibiting an improvement over the conventional scheme. In contrast to the conventional scheme of the photon transfer involving the photon reflections, in which the transfer efficiency depends on the coupling strength and input frequency tuning, the nonlinear saturated gain guarantees the stable state with the PT-symmetry and a self-selected frequency, which do not require active tuning to maintain optimal efficiency. The proposal of the nonlinear PT-symmetric system with nonlinear gain and loss provides a powerful platform for investigating intriguing properties prior to those in linear PT-symmetric systems, especially for compound photonic system.

Acknowledgements.
We acknowledge useful discussions with Xin-You Lü. This work was supported by Fundamental Research Funds for the Central Universities Grant No. 2020CQJQY-Z003.

References

  • (1) K. G. Makris, R. EI-Ganainy, D. N. Christodoulides, and Z. H . Musslimani, Phys. Rev. Lett. 100, 103904 (2008)
  • (2) S. Klaiman, U. Günther, and N. Moiseyev, Phys. Rev. Lett. 101, 080402 (2008).
  • (3) T. Kottos, Nat. Phys. 6, 166 (2010).
  • (4) X. Y. Lü, H. Jing, J. Y. Ma, and Y. Wu, Phys. Rev. Lett. 114, 253601 (2015).
  • (5) S. V. Suchkov, S. V. Dmitriev, B. A. Malomed, and Y. S. Kivshar, Phys. Rev. A 85, 033825 (2012).
  • (6) C. M. Bender, and S. Boettcher, Phys. Rev. Lett. 80, 5243 (1998).
  • (7) S. Klaiman, U. Güuther, and N. Moiseyev, Phys. Rev. Lett. 101, 080402 (2008).
  • (8) L. Jin and Z. Song, Phys. Rev. A 81, 032109 (2010).
  • (9) A. U. Hassan, H. Hodaei, M. A. Miri, M. Khajavikhan, and D. N. Christodoulides, Phys. Rev. A 92, 063807 (3015).
  • (10) Y. Sun, W. Tan, H. Q. Li, J. Li, and H. Chen, Phys. Rev. Lett. 112, 143903 (2014).
  • (11) C. E. Rüter, et al., Nat. Phys. 6, 192 (2010).
  • (12) B. Peng, et al., Nature Phys. 10, 394 (2014).
  • (13) L. Feng, et al., Nat. Mater. 12, 108 (2013).
  • (14) H. Jing, S. K. Özdemir, X. Y. Lü, J. Zhang, L. Yang, and F. Nori, Phys. Rev. Lett. 113, 053604 (2014).
  • (15) J. Fan, and L. Zhu, Opt. Express 20,20790 (2012).
  • (16) B. Peng, S. K. Özdemir, J. Zhu, and L. Yang, Opt. Lett. 37, 3435 (2012).
  • (17) C. M. Bender, D. C. Brody, and H. F. Jone, Phys. Rev. D 70, 025001 (2004).
  • (18) S. Longhi, Phys. Rev. Lett. 103, 123601 (2009).
  • (19) B. Zhen, et al., Nature (London) 525, 354 (2015).
  • (20) I. Rotter, J. Phys. A 42, 1-51 (2009).
  • (21) C. Dembowski, et al., Phys. Rev. Lett. 90, 034101 (2003).
  • (22) L. Chang, et al., Nat. Photon. 8, 524 (2014).
  • (23) R. EI-Ganainy, K. G. Makris, D. N. Christodoulides, and Z. H. Mussilimani, Opt. Lett. 32, 2632 (2007).
  • (24) Y. D. Chong, L. Ge, and A. D. Stone, Phys. Rev. Lett. 196, 093902 (2011).
  • (25) L. Jin, and F. Xin, Phys. Rev. A 96, 043821 (2017).
  • (26) A. Guo, et al., Phys. Rev. Lett. 103,093902 (2009).
  • (27) S. Assawaworrarit, X. F. Yu, and S. H. Fan, Nature 546, 387 (2017).
  • (28) A. U. Hassan, H. Hodaei, M. A. Miri, M. Khajavikhan, and D. N. Christodoulides, Phys. Rev. A 92, 63807 (2015).
  • (29) L. Ge, R. EI-Ganainy, Sci. Rep. 6, 24889 (2016).
  • (30) G. S. Agarwal, and K. Qu, Phys. Rev. A 85, 031802 (2012).
  • (31) X. Y. Lü, J. Q. Liao, L. Tian, and F. Nori, Phys. Rev. A 91, 013834 (2015).
  • (32) M. Wang, X. Y. Lü, Y. D. Wang, J. Q. You, and Y. Wu, Phys. Rev. A 94, 053807 (2016).