Controlling quasi-parametric amplifications: From multiple 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions to non-Hermitian sensing

Xiaoxiong Wu State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China    Kai Bai Key Laboratory of Artificial Micro- and Nano-structures of Ministry of Education and School of Physics and Technology, Wuhan University, Wuhan 430072, China    Penghong Yu State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China    Zhaohui Dong State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China    Yanyan He State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China    Jingui Ma Key Laboratory for Laser Plasmas (MOE), Collaborative Innovation Center of IFSA (CICIFSA), School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China    Vladislav V. Yakovlev Texas A&M University, College Station, TX 77843, USA    Meng Xiao phmxiao@whu.edu.cn Key Laboratory of Artificial Micro- and Nano-structures of Ministry of Education and School of Physics and Technology, Wuhan University, Wuhan 430072, China Wuhan Institute of Quantum Technology, Wuhan 430206, China    Xianfeng Chen State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China Jinan Institute of Quantum Technology, Jinan 250101, China Collaborative Innovation Center of Light Manipulation and Applications, Shandong Normal University, Jinan 250358, China    Luqi Yuan yuanluqi@sjtu.edu.cn State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
(July 4, 2024)
Abstract

Quasi-parametric amplification (QPA) is a nonlinear interaction in which the idler wave is depleted through some loss mechanism. QPA plays an important role in signal amplification in ultrafast photonics and quantum light generation. The QPA process has a number of features characterized by the non-Hermitian parity-time (𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T) symmetry. In this report, we explore new interaction regimes and uncover multiple 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions in such QPA process where transitions are particularly sensitive to external parameters. In particular, we demonstrate the feasibility of detection of 10βˆ’11superscript101110^{-11}10 start_POSTSUPERSCRIPT - 11 end_POSTSUPERSCRIPT inhomogeneities of the doped absorber, which is order of magnitude more sensitive than similar measurements performed in a linear absorption regime. In doing so, we reveal a family of 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions appearing in the QPA process and provide a novel nonlinear optical sensing mechanism for precise optical measurements.

I Introduction

In the past decade, non-Hermiticity introduced by the nonreciprocal coupling or exchange with surrounding environment has ignited tremendous research interest, and explorations of non-Hermitian physics using light or waves have attracted flourishing attention where various exotic phenomena and many unique features different from those in conventional Hermitian systems are explored [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. One of the key findings is the exceptional point (EP) surpporting parity-time (𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T) symmetry with balanced loss and gain distributions in space, i.e., [𝒫⁒𝒯,HP⁒T]=0𝒫𝒯subscript𝐻𝑃𝑇0[\mathcal{PT},H_{PT}]=0[ caligraphic_P caligraphic_T , italic_H start_POSTSUBSCRIPT italic_P italic_T end_POSTSUBSCRIPT ] = 0, where multiple bands coalesce with each other at EP, and afterwards fruitful applications have been proposed including unidirectional invisibility [11, 12, 13, 14] and single-mode lasing [15, 16, 17]. Due to the extreme responsivity of the EP from the surrounding change of environment [18, 19, 20, 21, 22], it has been proposed that one can use EP to perform sensing in photonics [23, 24]. Although there are some debates on the efficiency [18, 25, 4, 26, 27, 28], the idea of using its nonlinear dependence of parameters near EP is fundamentally intriguing.

Most of previous studies on the non-Hermitian EP sensing fall in photonic platforms. Recently it has been noticed that the nonlinear wave-mixing process in the ultrafast science may also support the non-Hermitian engineering, where the property of EP has been discussed [29, 30, 31, 32]. In particular, the connection between non-Hermiticity and the nonlinear wave-mixing processes such as optical parametric amplifications has been unveiled [33, 34, 35, 36, 29, 37, 30], which provides deeper understanding on the nonlinear photon conversion and proposes the quasi-parametric amplification (QPA) with high nonlinear conversion efficiency [38, 33, 34, 35, 39, 36]. Yet, so far, the behaviors of phase transitions and EP in this QPA system still lack further elaborate studies, especially for the dynamical evolution and its potential application with the use of the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transition such as aforementioned sensing. Especially, high-precision sensing has remarkable boost for scientific researches [40, 41], while QPA system may serve a unique candidate to detect the sensing information with very weak fluctuations.

In this report, we investigate a broad class of novel QPA systems including the parametric down-conversion (PDC) process that generates signal photons and linear depletion that absorbs idler photons. This hybrid process can be associated with the non-Hermitian 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry system dependent on the signal photon flux density along the propagation distance. Under different parametric conditions, multiple 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions exist, which results from the competition between the absorption of the idler photon and the signal amplification from PDC. We find that multiple transitions around EPs provides a dramatic variation of signal photons with sharp dips. Such unique feature can be utilized for a sensing detection to measure the homogeneousness of the doped absorber, i.e., very small perturbations can be measured with ultra-sensitivity, which can be verified in our simulations. Thus, our work opens a new perspective on the relationship between the nonlinear wave-mixing process and non-Hermiticity and unpacked an important avenue for the ultra-sensitive sensing at EPs in nonlinear optics.

II Model and theory

Refer to caption
Figure 1: Principle of QPA based on nonlinear PDC and linear absorption. (a) Schematic of a typical PDC process where the red, green and blue arrows represent pump (Ο‰asubscriptπœ”π‘Ž\omega_{a}italic_Ο‰ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT), signal (Ο‰bsubscriptπœ”π‘\omega_{b}italic_Ο‰ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT) and idler beams (Ο‰csubscriptπœ”π‘\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT), respectively. The pump, signal, and/or idler beams are injected into the nonlinear crystal which includes both the second-order nonlinearity Ο‡n⁒r(2)subscriptsuperscriptπœ’2π‘›π‘Ÿ\chi^{(2)}_{nr}italic_Ο‡ start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n italic_r end_POSTSUBSCRIPT and the linear absorber with χα(1)superscriptsubscriptπœ’π›Ό1\chi_{\alpha}^{(1)}italic_Ο‡ start_POSTSUBSCRIPT italic_Ξ± end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT. The panel below the schematic represents measuring intensities for crystal either with the uniform doping or with the doping with minor fluctuations, which reveals the sensing concept using EPs. (b) Illustration of PDC including a conversion between one pump photon Ο‰asubscriptπœ”π‘Ž\omega_{a}italic_Ο‰ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT, one signal photon Ο‰bsubscriptπœ”π‘\omega_{b}italic_Ο‰ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT and one idler photon Ο‰csubscriptπœ”π‘\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, as well as the linear absorption of the idler photon.

We consider a PDC process where light is injected into a nonlinear optical crystal with the second-order nonlinearity Ο‡n⁒r(2)subscriptsuperscriptπœ’2π‘›π‘Ÿ\chi^{(2)}_{nr}italic_Ο‡ start_POSTSUPERSCRIPT ( 2 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_n italic_r end_POSTSUBSCRIPT, doped with linear absorbers with Ο‡a(1)superscriptsubscriptπœ’π‘Ž1\chi_{a}^{(1)}italic_Ο‡ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 1 ) end_POSTSUPERSCRIPT. As depicted in Fig. 1, such second-order nonlinear process can support the photon conversion from a pump photon at the frequency Ο‰asubscriptπœ”π‘Ž\omega_{a}italic_Ο‰ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT to two photons, i.e., the signal photon at Ο‰bsubscriptπœ”π‘\omega_{b}italic_Ο‰ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT and the idler photon at Ο‰csubscriptπœ”π‘\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT (Ο‰a=Ο‰b+Ο‰csubscriptπœ”π‘Žsubscriptπœ”π‘subscriptπœ”π‘\omega_{a}=\omega_{b}+\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = italic_Ο‰ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT), and vice versa. The doped absorber is chosen to have the energy transition close to Ο‰csubscriptπœ”π‘\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, so it plays the role of loss for photons at Ο‰csubscriptπœ”π‘\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT. Such linear depletion of idler photons can lead to QPA if one injects strong pump beam together with a weak signal probe, where the nonlinear process can be operated only inside the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-broken phase to reach the high-efficient generation of the signal [38, 33, 35, 34, 36, 30]. In this work, different from Ref. [38, 33, 35, 34, 36, 30], we stick ourselves to this system but explore the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions based on different combination of inject beams and then point out the potential for sensing in nonlinear optical processes.

For the light propagation along the z𝑧zitalic_z direction in the nonlinear crystal, the electric field can be described by [35]

βˆ‡2E⁒(x,z,t)βˆ’nr2c02β’βˆ‚2βˆ‚t2⁒E⁒(x,z,t)=1Ο΅0⁒c02β’βˆ‚2PN⁒Lβˆ‚t2,superscriptβˆ‡2𝐸π‘₯𝑧𝑑superscriptsubscriptπ‘›π‘Ÿ2superscriptsubscript𝑐02superscript2superscript𝑑2𝐸π‘₯𝑧𝑑1subscriptitalic-Ο΅0superscriptsubscript𝑐02superscript2subscript𝑃𝑁𝐿superscript𝑑2\displaystyle\nabla^{2}E(x,z,t)-\frac{{n_{r}}^{2}}{c_{0}^{2}}\frac{\partial^{2% }}{\partial t^{2}}E(x,z,t)=\frac{1}{\epsilon_{0}c_{0}^{2}}\frac{\partial^{2}P_% {NL}}{\partial t^{2}},βˆ‡ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_E ( italic_x , italic_z , italic_t ) - divide start_ARG italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG βˆ‚ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG βˆ‚ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_E ( italic_x , italic_z , italic_t ) = divide start_ARG 1 end_ARG start_ARG italic_Ο΅ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG βˆ‚ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_N italic_L end_POSTSUBSCRIPT end_ARG start_ARG βˆ‚ italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (1)

where c0subscript𝑐0c_{0}italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the speed of light in vacuum, nrsubscriptπ‘›π‘Ÿn_{r}italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT describes the refractive index, and PN⁒L=Ο΅0⁒de⁒f⁒f⁒|E|2subscript𝑃𝑁𝐿subscriptitalic-Ο΅0subscript𝑑𝑒𝑓𝑓superscript𝐸2P_{NL}=\epsilon_{0}d_{eff}|E|^{2}italic_P start_POSTSUBSCRIPT italic_N italic_L end_POSTSUBSCRIPT = italic_Ο΅ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_d start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT | italic_E | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is the nonlinear polarization with the second-order nonlinear coefficient de⁒f⁒fsubscript𝑑𝑒𝑓𝑓d_{eff}italic_d start_POSTSUBSCRIPT italic_e italic_f italic_f end_POSTSUBSCRIPT and dielectric constant in vacuum Ο΅0subscriptitalic-Ο΅0\epsilon_{0}italic_Ο΅ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The electric field can be further expanded as E=12β’βˆ‘jAj⁒(z)⁒Bj⁒(x)⁒ei⁒(Ο‰j⁒tβˆ’Ξ²j⁒z)+c.c.formulae-sequence𝐸12subscript𝑗subscript𝐴𝑗𝑧subscript𝐡𝑗π‘₯superscript𝑒𝑖subscriptπœ”π‘—π‘‘subscript𝛽𝑗𝑧ccE=\frac{1}{2}\sum\limits_{j}A_{j}(z)B_{j}(x)e^{i(\omega_{j}t-\beta_{j}z)}+\rm{% c.c.}italic_E = divide start_ARG 1 end_ARG start_ARG 2 end_ARG βˆ‘ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) italic_B start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x ) italic_e start_POSTSUPERSCRIPT italic_i ( italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_t - italic_Ξ² start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_z ) end_POSTSUPERSCRIPT + roman_c . roman_c . where j=a,b,cπ‘—π‘Žπ‘π‘j=a,b,citalic_j = italic_a , italic_b , italic_c, Ξ²j=nr⁒(Ο‰j)⁒β0subscript𝛽𝑗subscriptπ‘›π‘Ÿsubscriptπœ”π‘—subscript𝛽0\beta_{j}=n_{r}(\omega_{j})\beta_{0}italic_Ξ² start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT = italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) italic_Ξ² start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT giving the wave vector at the frequency Ο‰jsubscriptπœ”π‘—\omega_{j}italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT, Bj⁒(x)subscript𝐡𝑗π‘₯B_{j}(x)italic_B start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_x ) is the transverse modal profile, and Aj⁒(z)subscript𝐴𝑗𝑧A_{j}(z)italic_A start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) is the field amplitude for the pump, signal, or idler field respectively [34, 35]. Then, together with the slowly varying envelope approximation and the rotating-wave approximation, Eq. (1) becomes [33, 39, 30]

{d⁒ua⁒(z)/d⁒z=i⁒Γ⁒ub⁒(z)⁒uc⁒(z)⁒eβˆ’i⁒Δ⁒β⁒z,d⁒ub⁒(z)/d⁒z=i⁒Γ⁒ua⁒(z)⁒ucβˆ—β’(z)⁒ei⁒Δ⁒β⁒z,d⁒uc⁒(z)/d⁒z=i⁒Γ⁒ua⁒(z)⁒ubβˆ—β’(z)⁒ei⁒Δ⁒β⁒zβˆ’Ξ±β’uc⁒(z),cases𝑑subscriptπ‘’π‘Žπ‘§π‘‘π‘§π‘–Ξ“subscript𝑒𝑏𝑧subscript𝑒𝑐𝑧superscript𝑒𝑖Δ𝛽𝑧𝑑subscript𝑒𝑏𝑧𝑑𝑧𝑖Γsubscriptπ‘’π‘Žπ‘§superscriptsubscript𝑒𝑐𝑧superscript𝑒𝑖Δ𝛽𝑧𝑑subscript𝑒𝑐𝑧𝑑𝑧𝑖Γsubscriptπ‘’π‘Žπ‘§superscriptsubscript𝑒𝑏𝑧superscript𝑒𝑖Δ𝛽𝑧𝛼subscript𝑒𝑐𝑧\left\{\begin{array}[]{l}{du_{a}(z)}/{dz}=i{\Gamma}{u_{b}}(z){u_{c}}(z)e^{-i% \Delta\beta z},\\ {du_{b}(z)}/{dz}=i{\Gamma}{u_{a}}(z)u_{c}^{*}(z)e^{i\Delta\beta z},\\ {du_{c}(z)}/{dz}=i{\Gamma}{u_{a}}(z)u_{b}^{*}(z)e^{i\Delta\beta z}-\alpha{u_{c% }}(z),\end{array}\right.{ start_ARRAY start_ROW start_CELL italic_d italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) / italic_d italic_z = italic_i roman_Ξ“ italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z ) italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) italic_e start_POSTSUPERSCRIPT - italic_i roman_Ξ” italic_Ξ² italic_z end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_d italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z ) / italic_d italic_z = italic_i roman_Ξ“ italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT βˆ— end_POSTSUPERSCRIPT ( italic_z ) italic_e start_POSTSUPERSCRIPT italic_i roman_Ξ” italic_Ξ² italic_z end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_d italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) / italic_d italic_z = italic_i roman_Ξ“ italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT βˆ— end_POSTSUPERSCRIPT ( italic_z ) italic_e start_POSTSUPERSCRIPT italic_i roman_Ξ” italic_Ξ² italic_z end_POSTSUPERSCRIPT - italic_Ξ± italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) , end_CELL end_ROW end_ARRAY (2)

where ΓΓ\Gammaroman_Ξ“ represents the nonlinear coupling coefficient of the PDC process, Ξ±β‰₯0𝛼0\alpha\geq 0italic_Ξ± β‰₯ 0 is the linear absorption coefficient, Δ⁒β=Ξ²aβˆ’(Ξ²b+Ξ²c)Δ𝛽subscriptπ›½π‘Žsubscript𝛽𝑏subscript𝛽𝑐\Delta\beta=\beta_{a}-(\beta_{b}+\beta_{c})roman_Ξ” italic_Ξ² = italic_Ξ² start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - ( italic_Ξ² start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_Ξ² start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ), and uj⁒(z)≑2⁒nr⁒(Ο‰j)⁒ϡ0⁒c0NT⁒(0)⁒ℏ⁒ωj⁒Aj⁒(z)subscript𝑒𝑗𝑧2subscriptπ‘›π‘Ÿsubscriptπœ”π‘—subscriptitalic-Ο΅0subscript𝑐0subscript𝑁𝑇0Planck-constant-over-2-pisubscriptπœ”π‘—subscript𝐴𝑗𝑧u_{j}(z)\equiv\sqrt{\frac{2n_{r}(\omega_{j})\epsilon_{0}c_{0}}{N_{T}(0)\hbar% \omega_{j}}}A_{j}(z)italic_u start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) ≑ square-root start_ARG divide start_ARG 2 italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) italic_Ο΅ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_N start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( 0 ) roman_ℏ italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG end_ARG italic_A start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) are the corresponding dimensionless field amplitude with NT⁒(0)subscript𝑁𝑇0N_{T}(0)italic_N start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( 0 ) being the total photon flux at z=0𝑧0z=0italic_z = 0. Here, the intensity of field at Ο‰jsubscriptπœ”π‘—\omega_{j}italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT is Ij⁒(z)=2⁒nr⁒(Ο‰j)⁒ϡ0⁒c0⁒Aj⁒(z)⁒Ajβˆ—β’(z)subscript𝐼𝑗𝑧2subscriptπ‘›π‘Ÿsubscriptπœ”π‘—subscriptitalic-Ο΅0subscript𝑐0subscript𝐴𝑗𝑧superscriptsubscript𝐴𝑗𝑧I_{j}(z)=2n_{r}(\omega_{j})\epsilon_{0}c_{0}A_{j}(z)A_{j}^{*}(z)italic_I start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) = 2 italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) italic_Ο΅ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_A start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) italic_A start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT βˆ— end_POSTSUPERSCRIPT ( italic_z ) so the total intensity of the light in the crystal gives IT⁒(z)=βˆ‘jIj⁒(z)subscript𝐼𝑇𝑧subscript𝑗subscript𝐼𝑗𝑧I_{T}(z)=\sum\limits_{j}I_{j}(z)italic_I start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_z ) = βˆ‘ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_I start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) and the photon flux is then defined as NT⁒(z)=βˆ‘jNj⁒(z)=βˆ‘jIj⁒(z)/ℏ⁒ωjsubscript𝑁𝑇𝑧subscript𝑗subscript𝑁𝑗𝑧subscript𝑗subscript𝐼𝑗𝑧Planck-constant-over-2-pisubscriptπœ”π‘—N_{T}(z)=\sum\limits_{j}N_{j}(z)=\sum\limits_{j}I_{j}(z)/{\hbar\omega_{j}}italic_N start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ( italic_z ) = βˆ‘ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_N start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) = βˆ‘ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT italic_I start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ( italic_z ) / roman_ℏ italic_Ο‰ start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT.

In the following, we take Δ⁒β=0Δ𝛽0\Delta\beta=0roman_Ξ” italic_Ξ² = 0 for the simplicity to explore the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions under the phase-matching condition. We define g⁒(z)≑Γ⁒ub⁒(z)𝑔𝑧Γsubscript𝑒𝑏𝑧{g}(z)\equiv\Gamma{u_{b}}(z)italic_g ( italic_z ) ≑ roman_Ξ“ italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z ) and {u~a⁒(z),u~c⁒(z)}≑eα⁒z/2⁒{ua⁒(z),uc⁒(z)}subscript~π‘’π‘Žπ‘§subscript~𝑒𝑐𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsuperscript𝑒𝛼𝑧2subscriptπ‘’π‘Žπ‘§subscript𝑒𝑐𝑧\{{\begin{array}[]{*{20}{c}}{{{\tilde{u}}_{a}}(z),}&{{{\tilde{u}}_{c}}(z)}\end% {array}}\}\equiv e^{\alpha z/2}\{u_{a}(z),u_{c}(z)\}{ start_ARRAY start_ROW start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) , end_CELL start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY } ≑ italic_e start_POSTSUPERSCRIPT italic_Ξ± italic_z / 2 end_POSTSUPERSCRIPT { italic_u start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) , italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) }, so Eq. (2) gets converted to

βˆ’i⁒dd⁒z⁒[u~a⁒(z)ub⁒(z)u~c⁒(z)]=[βˆ’i⁒α/20g⁒(z)Γ⁒ucβˆ—β’(z)00gβˆ—β’(z)0i⁒α/2]⁒[u~a⁒(z)ub⁒(z)u~c⁒(z)]≑H~Q⁒P⁒A⁒(z)⁒[u~a⁒(z)ub⁒(z)u~c⁒(z)].𝑖𝑑𝑑𝑧delimited-[]subscript~π‘’π‘Žπ‘§missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsubscript𝑒𝑏𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsubscript~𝑒𝑐𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressiondelimited-[]𝑖𝛼20𝑔𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionΞ“superscriptsubscript𝑒𝑐𝑧00missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsuperscript𝑔𝑧0𝑖𝛼2missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressiondelimited-[]subscript~π‘’π‘Žπ‘§missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsubscript𝑒𝑏𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsubscript~𝑒𝑐𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsubscript~𝐻𝑄𝑃𝐴𝑧delimited-[]subscript~π‘’π‘Žπ‘§missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsubscript𝑒𝑏𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionsubscript~𝑒𝑐𝑧missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpression-i\frac{d}{{dz}}\left[{\begin{array}[]{*{20}{c}}{{{\tilde{u}}_{a}}(z)}\\ {{u_{b}}(z)}\\ {{{\tilde{u}}_{c}}(z)}\end{array}}\right]=\left[{\begin{array}[]{*{20}{c}}{-i% \alpha/2}&0&{g(z)}\\ {\Gamma u_{c}^{*}(z)}&0&0\\ {{g^{*}}(z)}&0&{i\alpha/2}\end{array}}\right]\left[{\begin{array}[]{*{20}{c}}{% {{\tilde{u}}_{a}}(z)}\\ {{u_{b}}(z)}\\ {{{\tilde{u}}_{c}}(z)}\end{array}}\right]\equiv\tilde{H}_{QPA}(z)\left[{\begin% {array}[]{*{20}{c}}{{{\tilde{u}}_{a}}(z)}\\ {{u_{b}}(z)}\\ {{{\tilde{u}}_{c}}(z)}\end{array}}\right].- italic_i divide start_ARG italic_d end_ARG start_ARG italic_d italic_z end_ARG [ start_ARRAY start_ROW start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY ] = [ start_ARRAY start_ROW start_CELL - italic_i italic_Ξ± / 2 end_CELL start_CELL 0 end_CELL start_CELL italic_g ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL roman_Ξ“ italic_u start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT βˆ— end_POSTSUPERSCRIPT ( italic_z ) end_CELL start_CELL 0 end_CELL start_CELL 0 end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_g start_POSTSUPERSCRIPT βˆ— end_POSTSUPERSCRIPT ( italic_z ) end_CELL start_CELL 0 end_CELL start_CELL italic_i italic_Ξ± / 2 end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY ] [ start_ARRAY start_ROW start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY ] ≑ over~ start_ARG italic_H end_ARG start_POSTSUBSCRIPT italic_Q italic_P italic_A end_POSTSUBSCRIPT ( italic_z ) [ start_ARRAY start_ROW start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( italic_z ) end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY ] . (3)

H~Q⁒P⁒A⁒(z)subscript~𝐻𝑄𝑃𝐴𝑧\tilde{H}_{QPA}(z)over~ start_ARG italic_H end_ARG start_POSTSUBSCRIPT italic_Q italic_P italic_A end_POSTSUBSCRIPT ( italic_z ) represents the dynamical Hamiltonian of such nonlinear QPA system with three eigenvalues Ξ»~0,Β±=0subscript~πœ†0plus-or-minus0\tilde{\lambda}_{0,\pm}=0over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT 0 , Β± end_POSTSUBSCRIPT = 0 and Β±|g⁒(z)|2βˆ’(Ξ±2)2plus-or-minussuperscript𝑔𝑧2superscript𝛼22\pm\sqrt{{{\left|{{g}(z)}\right|}^{2}}-{{\left({\frac{\alpha}{2}}\right)}^{2}}}Β± square-root start_ARG | italic_g ( italic_z ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( divide start_ARG italic_Ξ± end_ARG start_ARG 2 end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG, together with corresponding eigenvectors v~0,±⁒(z)=[0,1,0]Tsubscript~𝑣0plus-or-minus𝑧superscriptdelimited-[]010missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpression𝑇\tilde{v}_{0,\pm}(z)={\left[{\begin{array}[]{*{20}{c}}{\begin{array}[]{*{20}{c% }}{0,}&1,&{0}\end{array}}\end{array}}\right]^{T}}over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT 0 , Β± end_POSTSUBSCRIPT ( italic_z ) = [ start_ARRAY start_ROW start_CELL start_ARRAY start_ROW start_CELL 0 , end_CELL start_CELL 1 , end_CELL start_CELL 0 end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT and 12⁒[1,0βˆ’i⁒|Ξ±|Β±βˆ’Ξ±2+|g|2gβˆ—]T12superscriptdelimited-[]10plus-or-minus𝑖𝛼superscript𝛼2superscript𝑔2superscript𝑔missing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpressionmissing-subexpression𝑇\frac{1}{{\sqrt{2}}}{\left[{\begin{array}[]{*{20}{c}}{1,}&0&{\frac{{-i\left|% \alpha\right|\pm\sqrt{-{\alpha^{2}}+{{\left|g\right|}^{2}}}}}{{g^{*}}}}\end{% array}}\right]^{T}}divide start_ARG 1 end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG [ start_ARRAY start_ROW start_CELL 1 , end_CELL start_CELL 0 end_CELL start_CELL divide start_ARG - italic_i | italic_Ξ± | Β± square-root start_ARG - italic_Ξ± start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + | italic_g | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG italic_g start_POSTSUPERSCRIPT βˆ— end_POSTSUPERSCRIPT end_ARG end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL start_CELL end_CELL end_ROW end_ARRAY ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT, respectively. One can find that there exists a pair of EPs at

Ξ±=Β±2⁒|g⁒(z)|,𝛼plus-or-minus2𝑔𝑧\alpha=\pm 2|g(z)|,italic_Ξ± = Β± 2 | italic_g ( italic_z ) | , (4)

labeled as zE⁒Psubscript𝑧𝐸𝑃z_{{}_{EP}}italic_z start_POSTSUBSCRIPT start_FLOATSUBSCRIPT italic_E italic_P end_FLOATSUBSCRIPT end_POSTSUBSCRIPT. The 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase of this system can be analyzed based on the basis {u~a,u~c}subscript~π‘’π‘Žsubscript~𝑒𝑐\{\tilde{u}_{a},\tilde{u}_{c}\}{ over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT }, where the amplitude of ub⁒(z)subscript𝑒𝑏𝑧u_{b}(z)italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z ) (or g⁒(z)𝑔𝑧g(z)italic_g ( italic_z )) determines the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions while the field propagates along z𝑧zitalic_z, i.e., oscillations at the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase and amplifications at the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-broken phase [33, 30].

III Analysis and results

III.1 Features of 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions under different parametric conditions

We now show the appearance of multiple 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions under different parametric conditions in the QPA system in Fig. 1. We first analyze the density ratio η⁒(Ξ±,z)≑n~a/(n~a+n~c)=na/(na+nc)πœ‚π›Όπ‘§subscript~π‘›π‘Žsubscript~π‘›π‘Žsubscript~𝑛𝑐subscriptπ‘›π‘Žsubscriptπ‘›π‘Žsubscript𝑛𝑐\eta(\alpha,z)\equiv\tilde{n}_{a}/(\tilde{n}_{a}+\tilde{n}_{c})=n_{a}/(n_{a}+n% _{c})italic_Ξ· ( italic_Ξ± , italic_z ) ≑ over~ start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT / ( over~ start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + over~ start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) = italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT / ( italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT + italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) on the basis {u~a,u~c}subscript~π‘’π‘Žsubscript~𝑒𝑐\{\tilde{u}_{a},\tilde{u}_{c}\}{ over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT , over~ start_ARG italic_u end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT }, and also the imaginary part of the eigenvalue Ξ»~Β±subscript~πœ†plus-or-minus\tilde{\lambda}_{\pm}over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT, where ni⁒(z)≑|ui⁒(z)|2subscript𝑛𝑖𝑧superscriptsubscript𝑒𝑖𝑧2n_{i}(z)\equiv|u_{i}(z)|^{2}italic_n start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_z ) ≑ | italic_u start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( italic_z ) | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT gives the dimensionless photon flux density for each field. Here, η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) labels the flux density ratio between the pump density and idler density that represents signatures of 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase. To plot η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) and ||||ImΞ»~Β±|\tilde{\lambda}_{\pm}|over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT | versus the propagation distance z𝑧zitalic_z, we numerically solve Eq. (2) with different choices of the absorption coefficient α𝛼\alphaitalic_Ξ±, and the initial combination of densities, i.e., {na⁒(0),nb⁒(0),nc⁒(0)}subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐0\left\{n_{a}(0),n_{b}(0),n_{c}(0)\right\}{ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) } (na⁒(0)+nb⁒(0)+nc⁒(0)=1subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐01n_{a}(0)+n_{b}(0)+n_{c}(0)=1italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) + italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) + italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) = 1). Hence ||||ImΞ»~Β±|\tilde{\lambda}_{\pm}|over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT | (dependent on ub⁒(z)subscript𝑒𝑏𝑧u_{b}(z)italic_u start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z )) and η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) (dependent on nasubscriptπ‘›π‘Žn_{a}italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT and ncsubscript𝑛𝑐n_{c}italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT) can be obtained.

Refer to caption
Figure 2: 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions of QPA system under different initial input conditions {na⁒(0),nb⁒(0),nc⁒(0)}subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐0\{n_{a}(0),n_{b}(0),n_{c}(0)\}{ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) }. (a) The evolution of density ratio η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ). In each subfigure, the horizontal axis is dimensionless propagation distance z𝑧zitalic_z and the vertical axis is linear absorption coefficient α𝛼\alphaitalic_Ξ±. (a1) and (a2) Zoom-in plots for η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) under {na⁒(0),nb⁒(0),nc⁒(0)}={0.2,0,0.8}subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐00.200.8\{n_{a}(0),n_{b}(0),n_{c}(0)\}=\{0.2,0,0.8\}{ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) } = { 0.2 , 0 , 0.8 } and {na⁒(0),nb⁒(0),nc⁒(0)}={0,0.2,0.8}subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐000.20.8\{n_{a}(0),n_{b}(0),n_{c}(0)\}=\{0,0.2,0.8\}{ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) } = { 0 , 0.2 , 0.8 }, respectively. (b) Imaginary eigenvalues ||||ImΞ»~Β±|\tilde{\lambda}_{\pm}|over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT | where white (purple) regions corresponds to the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry(-broken) phase, respectively. (b1) and (b2) Zoom-in plots for ||||ImΞ»~Β±|\tilde{\lambda}_{\pm}|over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT | under {na⁒(0),nb⁒(0),nc⁒(0)}={0.2,0,0.8}subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐00.200.8\{n_{a}(0),n_{b}(0),n_{c}(0)\}=\{0.2,0,0.8\}{ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) } = { 0.2 , 0 , 0.8 } and {na⁒(0),nb⁒(0),nc⁒(0)}={0,0.2,0.8}subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐000.20.8\{n_{a}(0),n_{b}(0),n_{c}(0)\}=\{0,0.2,0.8\}{ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) } = { 0 , 0.2 , 0.8 }, respectively.

Figs. 2(a)-(b) plot respectively, the variations of density ratios η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) and imaginary eigenvalues ||||Im(Ξ»~Β±)|(\tilde{\lambda}_{\pm})|( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) | on different nc⁒(0)subscript𝑛𝑐0n_{c}(0)italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) (from 0 to 1 horizontally) and nb⁒(0)subscript𝑛𝑏0n_{b}(0)italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) (from 0 to 1 vertically), where each subfigure shows the result at different z𝑧zitalic_z under various α𝛼\alphaitalic_Ξ±. One can see that Fig. 2 exhibits different characteristics of 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions during the light propagation when photon conversion under PDC (Ο‰aβ†’Ο‰b+Ο‰cβ†’subscriptπœ”π‘Žsubscriptπœ”π‘subscriptπœ”π‘\omega_{a}\to\omega_{b}+\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT β†’ italic_Ο‰ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT) and its inverse conversion (Ο‰b+Ο‰cβ†’Ο‰aβ†’subscriptπœ”π‘subscriptπœ”π‘subscriptπœ”π‘Ž\omega_{b}+\omega_{c}\to\omega_{a}italic_Ο‰ start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT β†’ italic_Ο‰ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT) occur simultaneously, i.e., η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) oscillates as a function of z𝑧zitalic_z at a fixed α𝛼\alphaitalic_Ξ±. Most subfigures in Fig. 2(a) show such oscillation features for small α𝛼\alphaitalic_Ξ±. For example, subfigures in the second row at nb⁒(0)=0.2subscript𝑛𝑏00.2n_{b}(0)=0.2italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) = 0.2 exhibit the features that photon conversions continuously happen as the absorption on photons at Ο‰csubscriptπœ”π‘\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is small. Subfigures in Fig. 2(b) also reflect such feature. The aforementioned regions in which η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) oscillates in Fig. 2(a) corresponds to the white regions in Fig. 2(b) where ||||Im(Ξ»~Β±)|β†’0(\tilde{\lambda}_{\pm})|\to 0( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) | β†’ 0, meaning the system is under the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase and supports two different ||||Re(Ξ»~Β±)|β‰ 0(\tilde{\lambda}_{\pm})|\neq 0( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) | β‰  0 resulting in oscillating dynamics. As the obvious distinction, one notes that for the choice of larger α𝛼\alphaitalic_Ξ±, η⁒(Ξ±,z)πœ‚π›Όπ‘§\eta(\alpha,z)italic_Ξ· ( italic_Ξ± , italic_z ) saturates towards unity in Fig. 2(a), indicating ncβ†’0β†’subscript𝑛𝑐0n_{c}\to 0italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT β†’ 0. This result corresponds to regions where ||||Im(Ξ»~Β±)|β‰ 0(\tilde{\lambda}_{\pm})|\neq 0( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) | β‰  0 but ||||Re(Ξ»~Β±)|β†’0(\tilde{\lambda}_{\pm})|\to 0( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) | β†’ 0 in Fig. 2(b) as the system goes into the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-broken phase, where the population of all three types of photons becomes stable. The critical point of α𝛼\alphaitalic_Ξ± where the phase transition occurs at large distance z𝑧zitalic_z gradually shifts to smaller value with nc⁒(0)subscript𝑛𝑐0n_{c}(0)italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) increasing. For example, for the case of nb⁒(0)=0.2subscript𝑛𝑏00.2n_{b}(0)=0.2italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) = 0.2, the increase of nc⁒(0)subscript𝑛𝑐0n_{c}(0)italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) leads to a larger portion of the input idler photons at Ο‰csubscriptπœ”π‘\omega_{c}italic_Ο‰ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT that promotes the photon inversion, which results in the fast inverse conversion to pump photons at Ο‰asubscriptπœ”π‘Ž\omega_{a}italic_Ο‰ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT under the absorption of idler photons. Hence pump photons become dominant and the region of η⁒(Ξ±,z)β†’1β†’πœ‚π›Όπ‘§1\eta(\alpha,z)\to 1italic_Ξ· ( italic_Ξ± , italic_z ) β†’ 1 occupies larger region in the plots in Fig. 2(a), which corresponds to the larger region of non-zero ||||Im(Ξ»~Β±)|(\tilde{\lambda}_{\pm})|( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) | in Fig. 2(b).

As the most striking feature, for cases at nc⁒(0)=0.8subscript𝑛𝑐00.8n_{c}(0)=0.8italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) = 0.8 with either na⁒(0)=0.2subscriptπ‘›π‘Ž00.2n_{a}(0)=0.2italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) = 0.2, nb⁒(0)=0subscript𝑛𝑏00n_{b}(0)=0italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) = 0 or na⁒(0)=0subscriptπ‘›π‘Ž00n_{a}(0)=0italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) = 0, nb⁒(0)=0.2subscript𝑛𝑏00.2n_{b}(0)=0.2italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) = 0.2 (zoomed-in in Figs. 2(a)-(b), respectively), there are multiple 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions in the small α𝛼\alphaitalic_Ξ± regime originated from the competition between PDC and its inverse conversion. As the immediate consequence at large z𝑧zitalic_z, with the continual absorption and the inverse conversion, this competition makes signal photons decrease and so Figs. 2(b1)-(b2) depict sharp peaks, which is also reflected as multiple staggered stripes in Figs. 2(a1)-(a2). Other details can be seen in Supplementary Material.

III.2 Relationships between 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions and photon flux densities

Refer to caption
Figure 3: (a) Imaginary eigenvalues ||||ImΞ»~Β±|\tilde{\lambda}_{\pm}|over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT | under the initial condition na⁒(0)=0.15subscriptπ‘›π‘Ž00.15n_{a}(0)=0.15italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) = 0.15 and nc⁒(0)=0.85subscript𝑛𝑐00.85n_{c}(0)=0.85italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) = 0.85. The dashed line corresponds to z=25⁒z0𝑧25subscript𝑧0z=25z_{0}italic_z = 25 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, with I and II in the vicinities of two transition points. (b) Photon flux density nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT as a function of α𝛼\alphaitalic_Ξ± (the vertical axis) plotted at z=25⁒z0𝑧25subscript𝑧0z=25z_{0}italic_z = 25 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Two sharp dips of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT can be seen at similar regimes I and II of α𝛼\alphaitalic_Ξ± respectively. (c1) and (d1) Phase rigidity |r|π‘Ÿ|r|| italic_r |, (c2) and (d2) Re(Ξ»~Β±)subscript~πœ†plus-or-minus(\tilde{\lambda}_{\pm})( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) of two dips versus the zoomed-in α𝛼\alphaitalic_Ξ± where the horizontal axis corresponds to j=[Ξ±βˆ’Ξ±I0⁒(II0)]/Δ⁒α𝑗delimited-[]𝛼subscript𝛼subscriptI0subscriptII0Δ𝛼j=\left[\alpha-\alpha_{\rm{I_{0}(II_{0})}}\right]/\Delta\alphaitalic_j = [ italic_Ξ± - italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT ] / roman_Ξ” italic_Ξ± with Ξ±I0⁒(II0)=0.235718⁒(0.113317)subscript𝛼subscriptI0subscriptII00.2357180.113317\alpha_{\rm{I_{0}(II_{0})}}=0.235718(0.113317)italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT = 0.235718 ( 0.113317 ) and Δ⁒α=10βˆ’6Δ𝛼superscript106\Delta\alpha=10^{-6}roman_Ξ” italic_Ξ± = 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT. In (c2) and (d2) black (red) lines denote to the eigenvalues Re(Ξ»~Β±)subscript~πœ†plus-or-minus(\tilde{\lambda}_{\pm})( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ).

To investigate the relationship between 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions and photon flux densities, we consider in detail the case of {na⁒(0),nb⁒(0),nc⁒(0)}={0.15,0,0.85}subscriptπ‘›π‘Ž0subscript𝑛𝑏0subscript𝑛𝑐00.1500.85\{n_{a}(0),n_{b}(0),n_{c}(0)\}=\{0.15,0,0.85\}{ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( 0 ) , italic_n start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 0 ) } = { 0.15 , 0 , 0.85 } and plot the imaginary eigenvalues ||||Im(Ξ»~Β±)|(\tilde{\lambda}_{\pm})|( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) | in Fig. 3(a). One can see similar multiple phase transitions as those in Figs. 2(b1)-(b2). At the distance z=25⁒z0𝑧25subscript𝑧0z=25z_{0}italic_z = 25 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, we plot the photon flux density of the signal density nb⁒(z=25⁒z0)subscript𝑛𝑏𝑧25subscript𝑧0n_{b}(z=25z_{0})italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( italic_z = 25 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) in Fig. 3(b). One can see depletion dips of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT at each value of α𝛼\alphaitalic_Ξ± that the phase transitions occur at Ξ±I0subscript𝛼subscriptI0\alpha_{\rm{I_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT=0.235718 and Ξ±II0subscript𝛼subscriptII0\alpha_{\rm{II_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT=0.113317 (labeled as the dips I and II) respectively in Fig. 3(a). Such dips are more sharper for larger α𝛼\alphaitalic_Ξ±, i.e., Ξ±>0.1𝛼0.1\alpha>0.1italic_Ξ± > 0.1, while for Ξ±<0.1𝛼0.1\alpha<0.1italic_Ξ± < 0.1, the photon flux of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT exhibits oscillation feature. Note that the number of significant digits in α𝛼\alphaitalic_Ξ± may affect the precision in the non-Hermitian sensing proposed in the next section, so we leave the number of the significant digits in α𝛼\alphaitalic_Ξ± upto six here to see its sensing capability in the theoretical limit. Decreasing the number of the significant digits in α𝛼\alphaitalic_Ξ± will not affect the sensing mechanism with EPs but might decrease the sensing precision.

Other than the imaginary eigenvalue plotted in Fig. 2(b), we also calculate the phase rigidity |r|=1/⟨v~Β±|v~Β±βŸ©π‘Ÿ1inner-productsubscript~𝑣plus-or-minussubscript~𝑣plus-or-minus|r|=1/\left\langle{\tilde{v}_{\pm}|\tilde{v}_{\pm}}\right\rangle| italic_r | = 1 / ⟨ over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT | over~ start_ARG italic_v end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ⟩ [42, 27, 43] and include the real eigenvalue Re(Ξ»~)~πœ†(\tilde{\lambda})( over~ start_ARG italic_Ξ» end_ARG ), both of which are plotted in Figs. 3(c1)-(c2) and (d1)-(d2) with the zoom-in α𝛼\alphaitalic_Ξ± in the vicinities of two phase transition dips (I & II). We can see the phase rigidity vanishes at the EPs, which can also been seen in the coalescence of the real eigenvalues in Figs. 3(c2)-(d2). It is worthy to note that the difference Δ⁒αE⁒PΞ”subscript𝛼𝐸𝑃\Delta\alpha_{EP}roman_Ξ” italic_Ξ± start_POSTSUBSCRIPT italic_E italic_P end_POSTSUBSCRIPT between two transition points in Figs. 3(c1)-(d2) is not equivalent to the full width at half maximum of dips of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT in Fig. 3(b) owing to non-immediate reaction to convert the signal density in such small scopes of 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-broken phase. Nevertheless, the sharp dip of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT between phase transitions holds the potential capability for the sensing under the small variation of surrounding conditions. In other words, such multiple phase transitions can provide a new perspective to explore non-Hermitian sensing in a nonlinear system.

III.3 Feasibility of non-Hermitian sensing from multiple phase transitions

As an example, we consider the problem of sensing the nearly homogeneousness of the doped absorber along z𝑧zitalic_z direction in the nonlinear system. Small spatial perturbation on α𝛼\alphaitalic_Ξ± may be not easy to be captured from a simple absorption measurement. In this case, the sharp dips in the photon flux density nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT provide a key candidate for sensing quantitatively. The idea is to tune the system operated at the regime between phase transitions (i.e., the dip I or II in Fig. 3(a)) and then measure nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT. For an absolute homogeneous distribution of α𝛼\alphaitalic_Ξ± in the z𝑧zitalic_z direction, one expects a measurement of the variation of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT to be zero. Nevertheless, due to the extreme sensitivity from 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions, small perturbation in α𝛼\alphaitalic_Ξ± can result in different measurements of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT.

Refer to caption
Figure 4: Photon flux densities of the signal photon (nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT) calculated with different amplitudes of perturbation (΢𝜁\zetaitalic_ΞΆ) with the average α𝛼\alphaitalic_Ξ± chosen as (a1)-(a3) Ξ±I0subscript𝛼subscriptI0\alpha_{\rm{I_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT for the dip I and (b1)-(b3) Ξ±II0subscript𝛼subscriptII0\alpha_{\rm{II_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT for the dip II in Fig. 3(a). Different distances of the nonlinear crystal (zm⁒a⁒xsubscriptπ‘§π‘šπ‘Žπ‘₯z_{max}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT) are taken in simulations, i.e., (a1) and (b1) zm⁒a⁒x=10⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯10subscript𝑧0z_{max}=10z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 10 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, (a2) and (b2) zm⁒a⁒x=20⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯20subscript𝑧0z_{max}=20z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 20 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, (a3) and (b3) zm⁒a⁒x=30⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯30subscript𝑧0z_{max}=30z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 30 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The grey shadow denotes the standard deviations from twenty sets of simulation results under the same parameters but different randomness. The red dashed line is the result of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT in the homogeneous crystal (ΞΆ=0𝜁0\zeta=0italic_ΞΆ = 0). The parameters are set to the condition of exact phase-matching at the center of the dip I (Ξ±I0subscript𝛼subscriptI0\alpha_{\rm{I_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT).

To calibrate the perturbation in α𝛼\alphaitalic_Ξ±, we define the amplitude of perturbation ΢𝜁\zetaitalic_ΞΆ and randomly set α∈[(1βˆ’ΞΆ)⁒α0,(1+ΞΆ)⁒α0]𝛼1𝜁subscript𝛼01𝜁subscript𝛼0\alpha\in[(1-\zeta)\alpha_{0},(1+\zeta)\alpha_{0}]italic_Ξ± ∈ [ ( 1 - italic_ΞΆ ) italic_Ξ± start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , ( 1 + italic_ΞΆ ) italic_Ξ± start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ] continuously in the nonlinear system, where Ξ±0subscript𝛼0\alpha_{0}italic_Ξ± start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the value of the homogeneous absorption (see more details in Supplementary Material). Figs. 4(a1)-(a3) and (b1)-(b3) show the statistics of the output nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT for the choice of different ΢𝜁\zetaitalic_ΞΆ with the choice of Ξ±I0⁒(II0)subscript𝛼subscriptI0subscriptII0\alpha_{\rm{I_{0}(II_{0})}}italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT, i.e., the dips I and II in Fig. 3(b). Three different distances of the nonlinear crystal (z=10⁒z0,20⁒z0,30⁒z0𝑧10subscript𝑧020subscript𝑧030subscript𝑧0z=10z_{0},20z_{0},30z_{0}italic_z = 10 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , 20 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , 30 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT) are taken corresponding to the same transition regime in Fig. 3(a), where larger distance gives sharper transition on α𝛼\alphaitalic_Ξ± but less sensitivity on ||||Im(Ξ»~Β±)|(\tilde{\lambda}_{\pm})|( over~ start_ARG italic_Ξ» end_ARG start_POSTSUBSCRIPT Β± end_POSTSUBSCRIPT ) |. For each choice of parameters, we count twenty sets of simulation data with the random set of α𝛼\alphaitalic_Ξ±. In Figs. 4(a1)-(a3), one can see that for small perturbation ΢𝜁\zetaitalic_ΞΆ, three choices of the crystal distance all give a stable measurement of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT. The threshold for ΢𝜁\zetaitalic_ΞΆ (labeled as ΞΆ0subscript𝜁0\zeta_{0}italic_ΞΆ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT) that the measurement of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT becomes unstable occurs at different values for three cases. Namely, the threshold ΞΆ0β‰ˆ10βˆ’1.5,10βˆ’4,10βˆ’6.5subscript𝜁0superscript101.5superscript104superscript106.5\zeta_{0}\approx 10^{-1.5},10^{-4},10^{-6.5}italic_ΞΆ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT β‰ˆ 10 start_POSTSUPERSCRIPT - 1.5 end_POSTSUPERSCRIPT , 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT , 10 start_POSTSUPERSCRIPT - 6.5 end_POSTSUPERSCRIPT for z=10⁒z0𝑧10subscript𝑧0z=10z_{0}italic_z = 10 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, 20⁒z020subscript𝑧020z_{0}20 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, 30⁒z030subscript𝑧030z_{0}30 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, respectively. It is clear to see that, in this choice of Ξ±0subscript𝛼0\alpha_{0}italic_Ξ± start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the longer distance the crystal is, the better sensitivity for the spacial perturbations on ΢𝜁\zetaitalic_ΞΆ we obtain. For z=30⁒z0𝑧30subscript𝑧0z=30z_{0}italic_z = 30 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT one can expect to measure the variation of the absorber upto 10βˆ’7superscript10710^{-7}10 start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT in Fig. 4(a3). One shall note that here we adjust parameters to satisfy the exact phase-matching condition at the center of the dip I and hence for other α𝛼\alphaitalic_Ξ± there are small phase-mismatching on wave vectors under the spacial perturbations. Other choices of parameters for the phase-matching condition at other values of α𝛼\alphaitalic_Ξ± can be found in Supplementary Material. As a comparison, under this phase-matching condition for the dip I, we also study cases with the choice of Ξ±II0subscript𝛼subscriptII0\alpha_{\rm{II_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT at the dip II. In Figs. 4(b1)-(b3), the threshold for ΞΆ0subscript𝜁0\zeta_{0}italic_ΞΆ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT occurs at larger perturbation for all three crystal distances, meaning the less sensitivity for this case.

Refer to caption
Figure 5: Photon flux densities of the signal photon (nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT) calculated with different amplitudes of perturbation (΢𝜁\zetaitalic_ΞΆ) with the average α𝛼\alphaitalic_Ξ± chosen as (a1)-(a3) Ξ±I0subscript𝛼subscriptI0\alpha_{\rm{I_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT for the dip I and (b1)-(b3) Ξ±II0subscript𝛼subscriptII0\alpha_{\rm{II_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT for the dip II in Fig. 3(a). Different distances of the nonlinear crystal (zm⁒a⁒xsubscriptπ‘§π‘šπ‘Žπ‘₯z_{max}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT) are taken in simulations, i.e., (a1) and (b1) zm⁒a⁒x=10⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯10subscript𝑧0z_{max}=10z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 10 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, (a2) and (b2) zm⁒a⁒x=20⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯20subscript𝑧0z_{max}=20z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 20 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, (a3) and (b3) zm⁒a⁒x=30⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯30subscript𝑧0z_{max}=30z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 30 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The grey shadow denotes the standard deviations from twenty sets of simulation results under the same parameters but different randomness. The red dashed line is the result of nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT in the homogeneous crystal (ΞΆ=0𝜁0\zeta=0italic_ΞΆ = 0), The parameters are set to the condition of exact phase-matching at the center of the dip II (Ξ±II0subscript𝛼subscriptII0\alpha_{\rm{II_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT).

One can also switch the phase-matching condition, i.e., to set it at the center of the dip II, and perform the simulations for different perturbations on ΢𝜁\zetaitalic_ΞΆ. The corresponding results are plotted in Fig. 5. In this case, for the choice of Ξ±II0subscript𝛼subscriptII0\alpha_{\rm{II_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT at the dip II, we can see ΞΆ0=10βˆ’1.5,10βˆ’4,10βˆ’7subscript𝜁0superscript101.5superscript104superscript107\zeta_{0}=10^{-1.5},10^{-4},10^{-7}italic_ΞΆ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 10 start_POSTSUPERSCRIPT - 1.5 end_POSTSUPERSCRIPT , 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT , 10 start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT for z=10⁒z0,20⁒z0,30⁒z0𝑧10subscript𝑧020subscript𝑧030subscript𝑧0z=10z_{0},20z_{0},30z_{0}italic_z = 10 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , 20 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , 30 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, respectively, while for the choice at the dip I with the phase-mismatching, ΞΆ0subscript𝜁0\zeta_{0}italic_ΞΆ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is roughly near 10βˆ’1.5superscript101.510^{-1.5}10 start_POSTSUPERSCRIPT - 1.5 end_POSTSUPERSCRIPT for all three crystal distances.

Although one can tell the threshold value of ΞΆ0subscript𝜁0\zeta_{0}italic_ΞΆ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT when it varies away from its stable value in Figs. 4-5, it is useful to give a quantitative measurement. Thus, we define the contrast ratio Ξ›=1βˆ’nΒ―bΞΆ=ΞΆT⁒(zm⁒a⁒x)na⁒(z=0)βˆ’nbΞΆ=0⁒(zm⁒a⁒x)Ξ›1superscriptsubscriptΒ―π‘›π‘πœsubscriptπœπ‘‡subscriptπ‘§π‘šπ‘Žπ‘₯subscriptπ‘›π‘Žπ‘§0superscriptsubscriptπ‘›π‘πœ0subscriptπ‘§π‘šπ‘Žπ‘₯\Lambda=1-\frac{{\bar{n}_{b}^{\zeta={\zeta_{T}}}({z_{max}})}}{{{n_{a}}(z=0)-n_% {b}^{\zeta=0}({z_{max}})}}roman_Ξ› = 1 - divide start_ARG overΒ― start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ΞΆ = italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) end_ARG start_ARG italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z = 0 ) - italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ΞΆ = 0 end_POSTSUPERSCRIPT ( italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) end_ARG, where na⁒(z=0)subscriptπ‘›π‘Žπ‘§0n_{a}(z=0)italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z = 0 ) is the input photon flux density for the pump and nbΞΆ=0⁒(zm⁒a⁒x)superscriptsubscriptπ‘›π‘πœ0subscriptπ‘§π‘šπ‘Žπ‘₯n_{b}^{\zeta=0}(z_{max})italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ΞΆ = 0 end_POSTSUPERSCRIPT ( italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) give the output signal photon flux density in the case of no perturbation on the absorption coefficient, and also ΞΆTsubscriptπœπ‘‡\zeta_{T}italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, i.e., the perceived precision, nΒ―bΞΆ=ΞΆT⁒(zm⁒a⁒x)=nbΞΆ=0⁒(zm⁒a⁒x)+[na⁒(z=0)βˆ’nbΞΆ=0⁒(zm⁒a⁒x)]β‹…5%superscriptsubscriptΒ―π‘›π‘πœsubscriptπœπ‘‡subscriptπ‘§π‘šπ‘Žπ‘₯superscriptsubscriptπ‘›π‘πœ0subscriptπ‘§π‘šπ‘Žπ‘₯β‹…delimited-[]subscriptπ‘›π‘Žπ‘§0superscriptsubscriptπ‘›π‘πœ0subscriptπ‘§π‘šπ‘Žπ‘₯percent5\bar{n}_{b}^{\zeta={\zeta_{T}}}({z_{max}})=n_{b}^{\zeta=0}({z_{max}})+\left[{{% n_{a}}(z=0)-n_{b}^{\zeta=0}({z_{max}})}\right]\cdot 5\%overΒ― start_ARG italic_n end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ΞΆ = italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) = italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ΞΆ = 0 end_POSTSUPERSCRIPT ( italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) + [ italic_n start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( italic_z = 0 ) - italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ΞΆ = 0 end_POSTSUPERSCRIPT ( italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT ) ] β‹… 5 %. If the contrast ratio satisfies Ξ›β†’1β†’Ξ›1\Lambda\to 1roman_Ξ› β†’ 1, the system holds a sharp dip for the relation between nbsubscript𝑛𝑏n_{b}italic_n start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT and α𝛼\alphaitalic_Ξ± and the sensing is supported to be distinguishable.

Refer to caption
Figure 6: The contrast ratio ΛΛ\Lambdaroman_Ξ› and the transitive perturbation rate ΞΆTsubscriptπœπ‘‡\zeta_{T}italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT versus different nonlinear crystal distances (zm⁒a⁒xsubscriptπ‘§π‘šπ‘Žπ‘₯z_{max}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT) from taking parameters meeting phase-matching conditions (a1), (b1) at the centers of the dip I and (a2), (b2) at the center of the dip II. Red and blue dots denote the results using Ξ±I0subscript𝛼subscriptI0\alpha_{\rm{I_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_I start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT at the dip I and Ξ±II0subscript𝛼subscriptII0\alpha_{\rm{II_{0}}}italic_Ξ± start_POSTSUBSCRIPT roman_II start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT at dip II respectively. The red shadows in (a1) and (b1) (the blue shadows in (a1) and (b2)) refers to the effective perceiving scopes for the dips I and II under perturbations, respectively.

We plot ΛΛ\Lambdaroman_Ξ› and ΞΆTsubscriptπœπ‘‡\zeta_{T}italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT with the propagation distance zm⁒a⁒xsubscriptπ‘§π‘šπ‘Žπ‘₯z_{max}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT in Fig. 6. In Figs. 6(a1)-(b1) for the phase-matching condition at the dip I, one can see that the contrast ratio for the dip I is close to unity from zm⁒a⁒x=5⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯5subscript𝑧0z_{max}=5z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 5 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT to 47.5⁒z047.5subscript𝑧047.5z_{0}47.5 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT because of a sharp peak in the 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transition. In this regime, the system is suitable for sensing. The corresponding quantity βˆ’--log(ΞΆT)subscriptπœπ‘‡(\zeta_{T})( italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ) increases linearly and we can obtain more sensitive precision at a larger distance. In particular, for zm⁒a⁒x=47.5⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯47.5subscript𝑧0z_{max}=47.5z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 47.5 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, we have Ξ›βˆΌ0.6similar-toΞ›0.6\Lambda\sim 0.6roman_Ξ› ∼ 0.6 and ΞΆT∼10βˆ’11similar-tosubscriptπœπ‘‡superscript1011\zeta_{T}\sim 10^{-11}italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT ∼ 10 start_POSTSUPERSCRIPT - 11 end_POSTSUPERSCRIPT, indicating ultra-sensitive measurement for the small perturbation on the absorption coefficient. As the comparison, when we choose α𝛼\alphaitalic_Ξ± near the dip II, ΛΛ\Lambdaroman_Ξ› drops quickly after zm⁒a⁒x>17.5⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯17.5subscript𝑧0z_{max}>17.5z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT > 17.5 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, with ΞΆTsubscriptπœπ‘‡\zeta_{T}italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT upto ∼10βˆ’1similar-toabsentsuperscript101\sim 10^{-1}∼ 10 start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT in Fig. 6(b1). Similar results can be found in Figs. 6(a2)-(b2), where we choose the phase-matching at the dip II. In this case, we obtain a larger contrast when α𝛼\alphaitalic_Ξ± is chosen at the dip II, with ΞΆTsubscriptπœπ‘‡\zeta_{T}italic_ΞΆ start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT reaching 10βˆ’5superscript10510^{-5}10 start_POSTSUPERSCRIPT - 5 end_POSTSUPERSCRIPT at zm⁒a⁒x=27.5⁒z0subscriptπ‘§π‘šπ‘Žπ‘₯27.5subscript𝑧0z_{max}=27.5z_{0}italic_z start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT = 27.5 italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. One can see that the sensing can be more sensitive when we tune the phase-matching at the dip I because of the larger transition regime shown in Fig. 3(a). More discussions on sensing results of cases when one sets the phase-mismatching at centers of dips are given in Supplementary Material.

IV Conclusion

In summary, we have shown that the dynamical evolution of the QPA is controlled through absorption parameter of the system and its spatial distribution. The complexity of nonlinear optical interactions and internal competition for photon conversions supports the abundance of 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions. We reveal that a small variation of surrounding conditions can directly cause evident changes in output photon flux densities at the phase transition around EPs, which enlightens towards sensing for the homogeneousness of the absorber. We showcase the ultra-sensitive measurement of small perturbations on the distribution of absorbers utilizing the EP and transition points. The proposed novel scheme therefore holds important promise for investigating multiple 𝒫⁒𝒯𝒫𝒯\mathcal{PT}caligraphic_P caligraphic_T-symmetry phase transitions as well as demonstrating EP sensing in nonlinear optics, and these sensitive dips can be potentially extended to evaluate thermal fluctuation due to rapidly changing light field distributions [44, 45, 46].

V Acknowledgments

This research was supported by the National Natural Science Foundation of China (12122407 and 12192252) and the National Key Research and Development Program of China (No. 2023YFA1407200 and No. 2021YFA1400900). V.Y. acknowledges support of the NIH (R21CA269099, R21GM142107, and R01GM127696) and AFOSR (FA9550-20-1-0366). L.Y. thanks the sponsorship from Yangyang Development Fund.

References

  • Feng et al. [2017] L. Feng, R. El-Ganainy, and L. Ge, Non-Hermitian photonics based on parity–time symmetry, Nature Photonics 11, 752 (2017).
  • Zhao and Feng [2018] H. Zhao and L. Feng, Parity–time symmetric photonics, National Science Review 5, 183 (2018).
  • El-Ganainy et al. [2018] R. El-Ganainy, K. G. Makris, M. Khajavikhan, Z. H. Musslimani, S. Rotter, and D. N. Christodoulides, Non-Hermitian physics and PT symmetry, Nature Physics 14, 11 (2018).
  • Γ–zdemir et al. [2019] Ş. K. Γ–zdemir, S. Rotter, F. Nori, and L. Yang, Parity–time symmetry and exceptional points in photonics, Nature Materials 18, 783 (2019).
  • Assawaworrarit and Fan [2020] S. Assawaworrarit and S. Fan, Robust and efficient wireless power transfer using a switch-mode implementation of a nonlinear parity–time symmetric circuit, Nature Electronics 3, 273 (2020).
  • Parto et al. [2021] M. Parto, Y. G. N. Liu, B. Bahari, M. Khajavikhan, and D. N. Christodoulides, Non-Hermitian and topological photonics: optics at an exceptional point, Nanophotonics 10, 403 (2021).
  • Ding et al. [2022] K. Ding, C. Fang, and G. Ma, Non-Hermitian topology and exceptional-point geometries, Nature Reviews Physics 4, 745 (2022).
  • Yan et al. [2023] Q. Yan, B. Zhao, R. Zhou, R. Ma, Q. Lyu, S. Chu, X. Hu, and Q. Gong, Advances and applications on non-Hermitian topological photonics, Nanophotonics 12, 2247 (2023).
  • Cheng et al. [2022] D. Cheng, B. Peng, M. Xiao, X. Chen, L. Yuan, and S. Fan, Truncation-dependent PT phase transition for the edge states of a two-dimensional non-Hermitian system, Physical Review B 105, L201105 (2022).
  • Wu et al. [2023] X. Wu, L. Wang, S. Chen, X. Chen, and L. Yuan, Transition characteristics of non-Hermitian skin effects in a Zigzag lattice without chiral symmetry, Advanced Physics Research 2, 2300007 (2023).
  • Lin et al. [2011] Z. Lin, H. Ramezani, T. Eichelkraut, T. Kottos, H. Cao, and D. N. Christodoulides, Unidirectional invisibility induced by PT-symmetric periodic structures, Physical Review Letters 106, 213901 (2011).
  • Yin and Zhang [2013] X. Yin and X. Zhang, Unidirectional light propagation at exceptional points, Nature Materials 12, 175 (2013).
  • Huang et al. [2017] Y. Huang, Y. Shen, C. Min, S. Fan, and G. Veronis, Unidirectional reflectionless light propagation at exceptional points, Nanophotonics 6, 977 (2017).
  • Zhou et al. [2023] Z. Zhou, B. Jia, N. Wang, X. Wang, and Y. Li, Observation of perfectly-chiral exceptional point via bound state in the continuum, Physical Review Letters 130, 116101 (2023).
  • Zhu et al. [2022] B. Zhu, Q. Wang, D. Leykam, H. Xue, Q. J. Wang, and Y. D. Chong, Anomalous single-mode lasing induced by nonlinearity and the non-Hermitian skin effect, Physical Review Letters 129, 013903 (2022).
  • Schumer et al. [2022] A. Schumer, Y. G. N. Liu, J. Leshin, L. Ding, Y. Alahmadi, A. U. Hassan, H. Nasari, S. Rotter, D. N. Christodoulides, P. LiKamWa, and M. Khajavikhan, Topological modes in a laser cavity through exceptional state transfer, Science 375, 884 (2022).
  • Li et al. [2023] A. Li, H. Wei, M. Cotrufo, W. Chen, S. Mann, X. Ni, B. Xu, J. Chen, J. Wang, S. Fan, et al., Exceptional points and non-Hermitian photonics at the nanoscale, Nature Nanotechnology 18, 706 (2023).
  • Hodaei et al. [2017] H. Hodaei, A. U. Hassan, S. Wittek, H. Garcia-Gracia, R. El-Ganainy, D. N. Christodoulides, and M. Khajavikhan, Enhanced sensitivity at higher-order exceptional points, Nature 548, 187 (2017).
  • Chen et al. [2017] W. Chen, Ş. Kaya Γ–zdemir, G. Zhao, J. Wiersig, and L. Yang, Exceptional points enhance sensing in an optical microcavity, Nature 548, 192 (2017).
  • Miri and Alu [2019] M.-A. Miri and A. Alu, Exceptional points in optics and photonics, Science 363, eaar7709 (2019).
  • Wiersig [2020] J. Wiersig, Review of exceptional point-based sensors, Photoics Research 8, 1457 (2020).
  • Peters and Rodriguez [2022] K. J. H. Peters and S. R. K. Rodriguez, Exceptional precision of a nonlinear optical sensor at a square-root singularity, Physical Review Letters 129, 013901 (2022).
  • Feng and Sun [2022] Z. Feng and X. Sun, Harnessing dynamical encircling of an exceptional point in anti-𝒫⁒𝒯𝒫𝒯\mathcal{P}\mathcal{T}caligraphic_P caligraphic_T-symmetric integrated photonic systems, Physical Review Letters 129, 273601 (2022).
  • Xu et al. [2023] Y. Xu, L. Li, H. Jeong, S. Kim, I. Kim, J. Rho, and Y. Liu, Subwavelength control of light transport at the exceptional point by non-Hermitian metagratings, Science advances 9, eadf3510 (2023).
  • Langbein [2018] W. Langbein, No exceptional precision of exceptional-point sensors, Physical Review A 98, 023805 (2018).
  • Chen et al. [2019] C. Chen, L. Jin, and R.-B. Liu, Sensitivity of parameter estimation near the exceptional point of a non-Hermitian system, New Journal of Physics 21, 083002 (2019).
  • Bai et al. [2022] K. Bai, L. Fang, T.-R. Liu, J.-Z. Li, D. Wan, and M. Xiao, Nonlinearity-enabled higher-order exceptional singularities with ultra-enhanced signal-to-noise ratio, National Science Review , nwac259 (2022).
  • Bai et al. [2024] K. Bai, T.-R. Liu, L. Fang, J.-Z. Li, C. Lin, D. Wan, and M. Xiao, Observation of nonlinear exceptional points with a complete basis in dynamics, Physical Review Letters 132, 073802 (2024).
  • Flemens et al. [2021] N. Flemens, N. Swenson, and J. Moses, Efficient parametric amplification via simultaneous second harmonic generation, Optics Express 29, 30590 (2021).
  • Flemens and Moses [2022] N. Flemens and J. Moses, Hermitian nonlinear wave mixing controlled by a PT-symmetric phase transition, Physical Review Letters 129, 153901 (2022).
  • Chen et al. [2023] Y. Chen, Z. Fu, B. Li, S. Peng, B. Zhu, G. Fan, Y. Liu, C. Ding, C. Jin, and Z. Tao, Phase-matched high-harmonic generation under nonadiabatic conditions: Model and experiment, Ultrafast Science 3, 0045 (2023).
  • Ni et al. [2023] F. Ni, H. Liu, Y. Zheng, and X. Chen, Nonlinear harmonic wave manipulation in nonlinear scattering medium via scattering-matrix method, Advanced Photonics 5, 046010 (2023).
  • El-Ganainy et al. [2015] R. El-Ganainy, J. I. Dadap, and R. M. Osgood, Optical parametric amplification via non-Hermitian phase matching, Optics letters 40, 5086 (2015).
  • Zhong et al. [2016] Q. Zhong, A. Ahmed, J. Dadap, R. Osgood Jr, and R. El-Ganainy, Parametric amplification in quasi-PT symmetric coupled waveguide structures, New Journal of Physics 18, 125006 (2016).
  • Miri and AlΓΉ [2016] M.-A. Miri and A. AlΓΉ, Nonlinearity-induced PT-symmetry without material gain, New Journal of Physics 18, 065001 (2016).
  • El-Ganainy et al. [2019] R. El-Ganainy, M. Khajavikhan, D. N. Christodoulides, and S. K. Ozdemir, The dawn of non-Hermitian optics, Communications Physics 2, 37 (2019).
  • Chen and Segev [2021] Z. Chen and M. Segev, Highlighting photonics: looking into the next decade, ELight 1, 2 (2021).
  • Ma et al. [2015] J. Ma, J. Wang, P. Yuan, G. Xie, K. Xiong, Y. Tu, X. Tu, E. Shi, Y. Zheng, and L. Qian, Quasi-parametric amplification of chirped pulses based on a sm3+-doped yttrium calcium oxyborate crystal, Optica 2, 1006 (2015).
  • Ma et al. [2017] J. Ma, J. Wang, B. Zhou, P. Yuan, G. Xie, K. Xiong, Y. Zheng, H. Zhu, and L. Qian, Broadband, efficient, and robust quasi-parametric chirped-pulse amplification, Optics Express 25, 25149 (2017).
  • Tam [1986] A. C. Tam, Applications of photoacoustic sensing techniques, Reviews of Modern Physics 58, 381 (1986).
  • Degen et al. [2017] C. L. Degen, F. Reinhard, and P. Cappellaro, Quantum sensing, Reviews of Modern Physics 89, 035002 (2017).
  • Ding et al. [2018] K. Ding, G. Ma, Z. Zhang, and C. T. Chan, Experimental demonstration of an anisotropic exceptional point, Physical Review letters 121, 085702 (2018).
  • Bai et al. [2023] K. Bai, J.-Z. Li, T.-R. Liu, L. Fang, D. Wan, and M. Xiao, Nonlinear exceptional points with a complete basis in dynamics, Physical Review Letters 130, 266901 (2023).
  • Gavartin et al. [2012] E. Gavartin, P. Verlot, and T. J. Kippenberg, A hybrid on-chip optomechanical transducer for ultrasensitive force measurements, Nature nanotechnology 7, 509 (2012).
  • Monzel and Sengupta [2016] C. Monzel and K. Sengupta, Measuring shape fluctuations in biological membranes, Journal of Physics D: Applied Physics 49, 243002 (2016).
  • Lu et al. [2019] P. Lu, N. Lalam, M. Badar, B. Liu, B. T. Chorpening, M. P. Buric, and P. R. Ohodnicki, Distributed optical fiber sensing: Review and perspective, Applied Physics Reviews 6, 041302 (2019).