Realizing strong photon blockade at exceptional points in the weak coupling regime

1 Introduction

As a significant area of quantum optics, the generation and manipulation of single photons have been making great strides in the past few decades and possess a wide array of applications in the fields of quantum communications [1], quantum cryptography [2], and quantum information processing [3-4]. One of the basic physical mechanisms for generating single photons is the photon blockade (PB) effect. What we called PB is that the first photon within an optical system will block the transmission of the second photon, leading to the phenomenon of photon antibunching in the system. This effect is first produced by Imamoglu et al. in 1997 [5], which plays key roles in exchanging and dealing with photonic quantum information [6–8].

So far, there are two main methods that have been used to generate strong photon blockade effects. One is the conventional photon blockade (CPB), and the other is the unconventional photon blockade (UPB). The CPB schemes require strong non-linear interactions between polaritons, which lead to a quantum anharmonic ladder in the energy spectrum. If a photon is tuned to resonantly excite the system from its ground state to the lowest excited states, the population of the two-photon state will be suppressed and only one photon is allowed in the system. The CPB effect has been achieved in various systems, including atom–cavity QED systems [10-12], cavity optomechanical systems [14], spinning Kerr cavity [16-17], and superconducting qubit systems [18-19].

Different from CPB, the physical mechanism of UPB relies on the quantum destructive interference between two or more quantum transition pathways in weakly non-linear systems. In the experiment, the phenomenon of UPB can be observed in the quantum dot–cavity QED system [21–23] and coupled superconducting resonators [24]. With the development of experiments, theoretical research has also been expanded in different quantum systems, for example, the couple cavities with second-order or third-order non-linearities [25–27], the cavity QED systems based on whispering-gallery-mode resonators [28-29], and the non-reciprocal devices such as spinning optomechanical systems [30].

Although both CPB and UPB can realize photon blockade, each type of PB has its own disadvantages in practice. Specifically, the realization of CPB depends on the strong light–matter interaction of the system, which is a big challenge in a few quantum systems. In particular, a fundamental CPB system typically requires a microcavity with a high Q factor [31], which is difficult to fabricate due to technical limitations. As for UPB, it may be hard to realize strong PB with large average photon numbers, resulting in the difficulty to obtain high-quality single-photons.

To solve these problems existing in the system, researchers try to achieve strong PB at the critical points, especially at exceptional points (EPs). EPs can be treated as critical points of the quantum phase transition from the PT-symmetric phase to the PT-symmetric-broken phase, where two or more eigenvalues and corresponding eigenvectors simultaneously coalesce [32–34]. EPs are one of the peculiar characteristics of the non-Hermitian systems [35], and there are lots of fascinating phenomena around these points such as single-mode lasers [36-37], unidirectional invisibility [38–40], sensitive enhancement [41–45], and topological energy transfer [46-47]. Very recently, EP-tuned purely quantum effects and their applications have been researched like non-reciprocal devices [48–50] and steady Bell-state generation [51].

Additionally, Mostafazadeh defined a new Hamiltonian that exists in the non-Hermitian systems, i.e., pseudo-Hermitian Hamiltonian [52–54]: a Hamiltonian H with a discrete spectrum that satisfies $H^{†}=UH{U}^{-1}$, where U is a linear Hermitian operator. The eigenvalues of this Hamiltonian are either real or complex conjugate pairs. So far, pseudo-Hermiticity plays an important role in the formation of higher-order exceptional points [55] and gives rise to a rich phenomenon in different fields of physics [56–58].

In this work, we theoretically propose a cavity QED system consisting of a gain single-mode cavity and a pair of two-level quantum emitters (QEs). First, we analytically demonstrate that the use of the gain cavity can provide relatively strong PB compared with the loss cavity even if the QE-cavity coupling strength is weak. For this case, we further prove that EP2 can be predicted in parameter space when the cavity and the QEs share the same frequency detuning. At EP2, we can obtain ultrastrong photon blockade effects with large mean photon numbers. Then, we prove that EP3 can be predicted in this system under the pseudo-Hermitian conditions. At this operator regime, the strong PB phenomenon can still be found. Compared with the PB effect that occurs at EP3 and EP2, we find that the PB effect enhanced at EP2 is stronger than that enhanced at EP3. Our proposal provides a new method to realize strong single-photon blockade in the weak coupling regime.

The paper is organized as follows: in Section 2, we give a detailed description of the physical model. By analytically solving a group of dynamics equations, we can obtain the expression of the second-order correlation function and mean photon number. Then, we discuss the origin of the PB effect in the normal loss cavity. In Section 3, we demonstrate that the strong PB effect can be achieved at EP2 in the weak coupling regime, and the physical mechanism can be analyzed in different quantum phase transition regions. In Section 4, we derive the pseudo-Hermitian condition for this considered system; both EP3 and EP2 can be predicted under specific parameter conditions. We study the EP3-enhanced strong PB phenomenon in Section 5. In Section 6, we compared the PB effect enhanced at EP3 and EP2 under different pseudo-Hermitian conditions. Finally, we give the conclusion of the whole work in Section 7.

2 Physical system of the two-level QE-cavity QED system

We construct two two-level QEs (e.g., two-level atoms, molecules, ions, or quantum dots) with the resonant frequency ${\omega }_i\left(i=\mathrm{1,2}\right)$ located in a single-mode cavity with the resonant frequency ${\omega }_c$. The ground (excited) QE state is expressed as ${\left|g\right.〉}_i\left({\left|e\right.〉}_i\right)$, and $g_i$ denotes the coupling strength between the QE and single-mode cavity. This cavity is coherently driven by a classical field with the Rabi frequency $\epsilon$ and pump frequency ${\omega }_d$, as illustrated in Figure 1. Using the rotating-wave approximation, the Hamiltonian of the system can be described as (setting $\hslash =1$)

where ${\Delta }_c={\omega }_c-{\omega }_d$ and ${\Delta }_i={\omega }_i-{\omega }_d$ are the cavity and QEi frequency detunings, respectively. Here, $a\left(a^{†}\right)$ is the annihilation (creation) operator of the cavity mode and ${\sigma }_i={\left|g\right.〉}_i〈\left.e\right|$ is the lowering operator of the ith two-level QE.

FIGURE 1. Schematic illustration of the two-level quantum emitter (QE)–cavity QED system with the cavity-mode frequency ${\omega }_c$ and QEi resonant frequency ${\omega }_i$. A classical field with the intensity $\epsilon$ and the angular frequency ${\omega }_d$ is used to drive the cavity. Here, ${\left|g\right.〉}_i$ and ${\left|e\right.〉}_i$ represent the ground state and excited state of QEi, respectively. $\kappa$ is the effective decay rate of the cavity, and ${\gamma }_i$ is the decay rate of QEi.

The dynamics of this cavity-driven QED system is governed by the quantum master equation:

where $\rho$ is the system density matrix and the Liouvillian operators $L_{\kappa }\left[\rho \right]=\frac{\kappa }{2}\left(2a\rho a^{†}-a^{†}a\rho -\rho a^{†}a\right)$ and $L_{\gamma }\left[\rho \right]={\displaystyle \sum _{i=1}^2}\frac{{\gamma }_i}{2}\left(2{\sigma }_i\rho {\sigma }_i^{†}-{\sigma }_i^{†}{\sigma }_i\rho -\rho {\sigma }_i^{†}{\sigma }_i\right)$ describe the cavity decay rate with $\kappa$ and the QEs with rate ${\gamma }_i$, respectively. In the case of weak driving, we can neglect the quantum jump term $o^{†}\rho o\left(o=a,{\sigma }_i\right)$ to obtain the effective non-Hermitian Hamiltonian:

Here, we take ${\gamma }_1={\gamma }_2=\gamma$ in the following calculation.

In order to give a better understanding of the PB effect from the physical point of view, we need to calculate the zero-delayed second-order correlation function $g^{\left(2\right)}\left(0\right)$. Under the weak driving assumption, i.e., $\epsilon \ll \gamma$ (in this paper, we only analyze the case of weak drive), we assume that the total excitation number of the system is truncated to 2. As a result, the time-dependent wave function can be written as

where $C_{rst}$ is the coefficient of the quantum state $\left|r\right.〉\left|s\right.〉\left|t\right.〉$. r stands for the photon number in the cavity. $\left|s\right.〉,\left|t\right.〉=0$ and 1 represent the two QEs in the ground states and excited states, respectively. First of all, it is necessary to obtain the steady-state solution of $C_{rst}$. We start from solving the Schrodinger equation $i\partial \left|\psi \left(t\right)\right.〉/\partial t=H_{eff}\left|\psi \left(t\right)\right.〉$ and then obtain a set of equations of motion for coefficients:

where ${\Delta }_c^{\prime }={\Delta }_c-i\frac{\kappa }{2},{\Delta }_i^{\prime }={\Delta }_i-i\frac{{\gamma }_i}{2}$. Under the weak driving condition, one can assume that $C_{000}\simeq 1\gg \left\{C_{100},C_{010},C_{001}\right\}\gg \left\{C_{101},C_{110},C_{011},C_{200}\right\}$. By setting $i\partial {\stackrel{.}{C}}_{rst}/\partial t=0$, we can easily obtain the steady-state solution of the aforementioned equations, which are expressed as

where $t_1={\Delta }_1^{\prime }+{\Delta }_c^{\prime }-g_2^2/\left({\Delta }_1^{\prime }+{\Delta }_2^{\prime }\right),t_2=-g_1{g}_2/\left({\Delta }_1^{\prime }+{\Delta }_2^{\prime }\right),t_3={\Delta }_2^{\prime }+{\Delta }_c^{\prime }-g_1^2/\left({\Delta }_1^{\prime }+{\Delta }_2^{\prime }\right)$ and the determinant

According to Eqs 6–8, the second-order correlation function $g^{\left(2\right)}\left(0\right)$ can be approximately yielded by $g^{\left(2\right)}\left(0\right)\approx 2{\left|C_{200}\right|}^2/{\left|C_{100}\right|}^4$ and the mean photon number in the cavity is $I_c=〈a^{†}a〉={\left|C_{100}\right|}^2+{\left|C_{101}\right|}^2+{\left|C_{110}\right|}^2+2{\left|C_{200}\right|}^2$. The expression of $g^{\left(2\right)}\left(0\right)$ can be expanded as

2.1 The PB in the QE-cavity QED system with the gain (loss) cavity

For simplicity, we assume that the two-level QE-cavity coupling strengths are the same (i.e., $g_1=g_2=g$) and that the QE and cavity frequency detunings are also identical (i.e., ${\Delta }_1={\Delta }_2={\Delta }_c=\Delta$) in this section. Next, we prove the validity of our previous calculations by comparing the analytical results with the numerical results given by Eq. 2 under the weak driving assumption shown in Figure 2. The analytical results are in good agreement with the numerical results for the second-order correlation function. In the same strong coupling regime, the use of the gain cavity can show more obvious photon blockade effects (See Figures 2A,C). As for the same weak coupling regime, the choice of the active or passive cavity has a little effect on PB effects (See Figures 2B,D). Therefore, it is worth presenting a new physical model for realizing strong PB effects at a specific area in the weak coupling regime.

FIGURE 2. In the system with the gain (dissipation) cavity, the logarithmic plots of the second-order correlation function $g^{\left(2\right)}\left(0\right)$ as a function of the normalized detuning $\Delta /\gamma$ for two cases: (A, C) the system in the strong coupling regime, while (B, D) in the weak coupling regime. Here, the other parameters are chosen as ${\Delta }_1={\Delta }_2={\Delta }_c=\Delta$. The driving strength takes as $\epsilon =0.1\gamma$ in the following figures.

For a better understanding of the physical mechanism of the PB effect in the system with the loss- or gain-cavity mode, we consider the system by utilizing the dressed-state representation. Specifically, this coupled system has a discrete spectrum consisting of a ladder-type dressed state, $\left|\pm ,n-1\right.〉$ with separated energy levels and other collective states are $\left|gg,n\right.〉$ and $\left|ee,n-2\right.〉$ [59]. Owing to the whole system being under the weak driving assumption, the principal quantum number of the system is truncated to $n\le 2$.

In the case of the loss cavity, when a photon is resonantly excited from the ground state of the system to the states of the lowest doublet, i.e., ${\mathrm{\Psi }}_1^{\left(\pm \right)}$, the absorption of the subsequent photon at the identical pump frequency will be blocked due to the large mismatch energy induced by energy-level anharmonicity (see Figure 3A). This is the blockade mechanism of the well-known CPB scheme.

FIGURE 3. (A) Anharmonic ladder-type energy-level diagram to explain the PB effect in the system with the loss cavity. (B) Quantum transition pathways of the system for different quantum states $\left|r\right.〉\left|s\right.〉\left|t\right.〉$.

Neglecting the dipole–dipole interaction (DDI) between two QEs, there is a direct transition pathway induced by the pump field, i.e., $\left|\mathrm{1,0,0}\right.〉\stackrel{\sqrt{2}\epsilon }{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉,$ and two indirect pathways induced by the QE-cavity coupling strengths, i.e., $\left|\mathrm{1,0,0}\right.〉\stackrel{g_1}{{\displaystyle \to }}\left|\mathrm{0,1,0}\right.〉\stackrel{\epsilon }{{\displaystyle \to }}\left|\mathrm{1,1,0}\right.〉\stackrel{\sqrt{2}{g}_1}{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉$ and $\left|\mathrm{1,0,0}\right.〉\stackrel{g_2}{{\displaystyle \to }}\left|\mathrm{0,0,1}\right.〉\stackrel{\epsilon }{{\displaystyle \to }}\left|\mathrm{1,0,1}\right.〉\stackrel{\sqrt{2}{g}_2}{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉$ (see Figure 3B). The direct transition pathway for the two-photon excited states will be forbidden, owing to the quantum destructive interference with the indirect pathways [32-60]. Consequently, the probability of the two-photon excited states will be reduced, which means that the weak coupling condition can still induce the PB effect. This is the blockade mechanism of the UPB scheme.

2.2 The exceptional point of the system with the gain cavity

In this section, we study the strong PB at a certain characteristic value in a weak coupling limit. This QE-cavity QED system can be described by the Hamiltonian without a driving term in the matrix form as

Here, we choose $\kappa =-\gamma$ and ${\Delta }_1={\Delta }_2={\Delta }_c=\Delta$. The eigenvalues of the system in single-photon space are expressed as

with the corresponding eigenvectors given by

According to Eq. 11, we find that when $g_1^2+g_2^2=\frac{{\gamma }^2}{4}$, the two eigenvalues ${\lambda }_{\pm }$ and the corresponding eigenvectors $\left|{\lambda }_{\pm }\right.〉$ will coalesce simultaneously, which indicates that the second-order exceptional point (EP2) will appear. Considering special circumstances, we think the two-level QE-cavity system has the same coupling strength, i.e., $g_1=g_2=g$. In Figure 4, we plot the real [see Figure 4A] and imaginary parts [see Figure 4B] of eigenvalues as a function of the coupling strength g. For the case of $g>\frac{\gamma }{2\sqrt{2}}$, the two eigenvalues ${\lambda }_{\pm }$ are real and non-degenerate, indicating that the system is in the PT-symmetric phase (PTSP). When $g<\frac{\gamma }{2\sqrt{2}}$, the eigenvalue is a pair of complex conjugates, which is the significant feature of the PT-symmetry-broken phase (PTBP). When $g=g_{EP2}=\frac{\gamma }{2\sqrt{2}}$, both the eigenvalues and corresponding eigenvectors are degenerated. The pink and yellow areas indicate the PT-symmetric phase and PT-symmetry-broken phase, respectively. Subsequently, we will analyze the single-photon blockade effect at and around EP2.

FIGURE 4. (A) Real and (B) imaginary parts of the eigenvalues are obtained by Eq. 12 with $\kappa =-\gamma$ versus the same coupling strength g (i.e.; $g_1=g_2=g$). Green dots show the position of EP2 when $\Delta =0$. The pink and yellow areas represent the PT-symmetric phase (PTSP) and PT-symmetry-broken phase (PTBP) area, respectively.

3 EP2-enhanced strong PB effects in the system

In order to demonstrate the optimal photon blockade at EP2, we plot both the mean photon number $I_c$ and second-order correlation function ${\mathrm{g}}^{\left(2\right)}\left(0\right)$ in two scenarios: i) the system with the gain cavity (i.e., $\kappa <0$) or ii) with the dissipation (i.e., $\kappa >0$) cavity, as shown in Figure 5. It is easy to find that when $\mathrm{g}=\frac{\gamma }{2\sqrt{2}}$, the ideal photon blockade will appear at $\kappa =-\gamma$, where the minimal value of ${\mathrm{g}}^{\left(2\right)}\left(0\right)$ and the maximum value of $I_c$ will be achieved simultaneously. For the case of the system with the gain cavity, there are two dips, one of which is located at EP2. However, for the cases of loss cavity, with the increase in the dissipation rate of the cavity, the PB effect will decrease rapidly. Therefore, compared with the gain cavity and loss cavity, the gain cavity provides a new possibility for imperfect photon blockade.

FIGURE 5. Second-order correlation function $g^{\left(2\right)}\left(0\right)$ and the mean photon number $I_c$ as functions of $\kappa$ in the system with the gain cavity (Panel (A)) and with the loss cavity (Panel (B)). Here, we take the coupling regime at $g=g_{EP2}=\frac{1}{2\sqrt{2}}\gamma$ and $\Delta =0,$ and the other parameters are same as those in the main text.

Furthermore, we can prove the optimal photon blockade at EP2 by calculating the value of ${\mathrm{g}}^{\left(2\right)}\left(0\right)$. By setting $g^{\left(2\right)}\left(0\right)\to 0$ in Eq. 9, we can seek out the positions where the pronounced photon antibunching phenomenon appears. The realization of the minimum value of ${\mathrm{g}}^{\left(2\right)}\left(0\right)$ requires $\left|D_1\right|\to 0$, that is,

Therefore, the strong PB effect can be obtained around $g=\frac{\gamma }{2\sqrt{2}}$ for $\Delta =0$. We noted that the coefficient of $C_{100}\to \infty$ when $\left|D_1\right|\to 0$, which indicates that $I_c$ will reach its maximum value and ${\mathrm{g}}^{\left(2\right)}\left(0\right)$ will reach its maximum value at $g=g_{EP2}$.

For the general case of $g_1\ne g_2$, we plot ${\mathrm{g}}^{\left(2\right)}\left(0\right)$ (see Figure 6A) and $I_c$ (see Figure 6B) in terms of $g_1$ and $g_2$, respectively. The black dashed lines of EP2 denote the optimal condition of the PB effect, i.e., $g_1^2+g_2^2=\frac{{\gamma }^2}{4}$. Under this optimal condition, ${\mathrm{g}}^{\left(2\right)}\left(0\right)$ of the PB effect may reduce to $10^{-4}$ and $I_c$ will increase to 0.32. For better proof that the PB effect enhanced at EP2, as shown in Figure 7A, we provide ${\mathrm{g}}^{\left(2\right)}\left(0\right)$ as functions of normalized detuning $\Delta /\mathrm{\gamma }$ and the same coupling strength $g/\gamma$. It is worth pointing out that the EP2 will be emerged at $g=g_{EP2}=\frac{\gamma }{2\sqrt{2}}<\gamma$, indicating that the strong PB phenomenon has occurred in a weak coupling limit. This result is demonstrated in Figure 7B.

FIGURE 6. When the coupling strengths are different, i.e., $g_1\ne g_2$, (A) the logarithmic plots of $g^{\left(2\right)}\left(0\right)$ and (B) $〈a^{†}a〉$ are as functions of the coupling strengths $g_1/\gamma$ and $g_2/\gamma$. The black dashed lines indicate the area of EP2 in the condition of $g_1^2+g_2^2=\frac{{\gamma }^2}{4}$. The other parameters are similar to those in Figure 4.

FIGURE 7. EP2-enhanced strong PB effect in the condition of $g_1=g_2=g$ and ${\Delta }_c={\Delta }_1={\Delta }_2=\Delta$. (A) Logarithmic plots of $g^{\left(2\right)}\left(0\right)$ as functions of normalized detuning $\Delta /\gamma$ and the same coupling strength $g/\gamma$. The red circle shows the position of EP2, where the extremely small $g^{\left(2\right)}\left(0\right)$ can be achieved. (B) Plots of $g^{\left(2\right)}\left(0\right)$ (see blue curves) and $I_c$ (see red curves) versus the coupling strength $g/\gamma$. The pink dashed line indicates the position of EP2. We choose ${\Delta }_c=0$ in panel (B).

In the following section, we study the reason for the PB effect in different regions. On one hand, to explain the PB effect when $g>\frac{\gamma }{2\sqrt{2}}$ (i.e., in the PT -symmetric region), we draw the ladder-type energy level of the quantum state, as shown in Figure 8A. According to the eigenvalues of the system, this physical mechanism can be understood as follows: when $g>\frac{\gamma }{2\sqrt{2}}$, the single-photon state is a single state. If the driving frequency ${\omega }_d$ equals to $E_1$, a photon is excited from the ground state to the single excited state resonantly, so the single-photon probability increases dramatically. However, the two-photon excitation probability may decrease due to the detuning.

FIGURE 8. Level diagram of the quantum states for the cavity QED with the gain cavity under the weak cavity-driven assumption: (A) dressed states and (B) quantum states.

On the other hand, when $g<\frac{\gamma }{2\sqrt{2}}$, there are two quantum paths that suffer from different efforts, i.e., the photon gain in the path $\left|\mathrm{1,0,0}\right.〉\stackrel{\sqrt{2}\epsilon }{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉$ and the photon loss in the two paths, i.e., $\left|\mathrm{1,0,0}\right.〉\stackrel{g_1}{{\displaystyle \to }}\left|\mathrm{0,1,0}\right.〉\stackrel{\epsilon }{{\displaystyle \to }}\left|\mathrm{1,1,0}\right.〉\stackrel{\sqrt{2}{g}_1}{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉$ and $\left|\mathrm{1,0,0}\right.〉\stackrel{g_2}{{\displaystyle \to }}\left|\mathrm{0,0,1}\right.〉\stackrel{\epsilon }{{\displaystyle \to }}\left|\mathrm{1,0,1}\right.〉\stackrel{\sqrt{2}{g}_2}{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉$. Therefore, the total loss and gain rate are unbalanced in this model. The gain of the photons in the cavity will compensate for the photon loss of the system, which makes the single-photon probability increase. This is the main reason for the PB effect occurring in the PT-symmetric-broken region. Through our calculation of the second-order correlation function by Eq. 9, one can easily adjust the coupling strength to realize the optimal photon blockade. Additionally, in the PT-symmetric-broken region (i.e., $g<g_{EP2}$), the photon blockade will be more obvious with the increase in the coupling strength g. However, in the PT-symmetric region (i.e., $g>g_{EP2}$), the photon antibunching phenomenon will be transformed into photon bunching with the increasing coupling strength.

4 Pseudo-Hermitian conditions for EP3 in the system

In addition to EP2, whether high-order exceptional points (i.e., EP3) will also affect the PB effect is worth studying. In this part, we further show the strong PB at EP3. First, we need to find the pseudo-Hermitian (pH) condition of the system. Following [52–54], the non-Hermitian Hamiltonian $H_r$ without DDI becomes pseudo-Hermitian when its eigenvalues satisfy one of the following conditions: i) all three eigenvalues are real or ii) one of the eigenvalues is real and the others are a pair of complex conjugates. Solving $\mathit{Det}\left(H_r-\mathrm{\Omega }I\right)=0$, i.e.,

where I is an identity matrix. We can obtain three eigenvalues from Eq. 14. Then, in order to meet the pseudo-Hermitian condition, both Eq. 14 and its complex conjugation expression, i.e., $\mathit{Det}\left(H_r^{*}-\mathrm{\Omega }I\right)=0,$ should have the same solutions. Solving these two equations gives rise to the pseudo-Hermitian conditions of the Hamiltonian (14) as

where ${\delta }_{\mathrm{1,2}}={\Delta }_{\mathrm{1,2}}-{\Delta }_c$ is the frequency detuning of the cavity and QEs. In the following calculation, we give the conditions

From the first condition in Eq. 16, it is easy to see that the gain cavity must be introduced to the QED system to keep the gain and loss balanced. From the last equation in Eq. 16, it should satisfy the condition of $\delta \ge 0$. By setting $\delta =0$, the relationship of minimal values of two-level QE-cavity coupling strength is given by ${\left[g_1^2+g_2^2\right]}_{\mathit{min}}=\frac{{\gamma }^2}{2}$, which is a basic condition that should be met in our system. When the system is pseudo-Hermitian, the characteristic equation $\mathit{Det}\left(H_r-\mathrm{\Omega }I\right)=0$ can be specifically written as

where

According to Cardano’s formula and methods [61], the solution of the characteristic equation in Eq. 16 can be determined by the discriminant

with $A=c_2^2-3c_1,B=c_1{c}_2-9c_0,\mathrm{a}\mathrm{n}\mathrm{d}C=c_1^2-3c_0{c}_2$. If $\eta <0$, Eq. 17 has three real solutions, but the solutions are one real root and a pair of complex conjugates if $\eta >0$. In the critical point at $\eta =0$, these three real solutions coalesce to the same value, i.e., $\mathrm{\Omega }={\mathrm{\Omega }}_{EP3}$. In other words, when $A=B=0$, EP3 will appear. If the Hamiltonian in Eq. 11 satisfies the conditions of (16), this non-Hermitian Hamiltonian will transform into a pseudo-Hermitian Hamiltonian.

To prove the aforementioned analysis, we plot the phase transition in Figure 9A, where the green and blank areas represent $\eta <0$ and $\eta >0$, respectively. The black lines, blue dashed lines, and a red curve denote the conditions of $\eta =0$, $B=0$, and $A=0$, respectively. The yellow crossing points produced by the black, blue, and red lines indicate EP3s in math. We can only find one EP3 in our system when two-level QE-cavity coupling strengths are the same, i.e., $g_1=g_2=g_{EP3}$ (see Figure 9B). Furthermore, we analytically prove this critical condition for the existence of EP3.

FIGURE 9. (A) Quantum phase of the discriminant in Eq. 19 under the pseudo-Hermitian conditions in Eq. 16 as a function of coupling strengths $g_1/\gamma$ and $g_2/\gamma$. (B) The yellow dots represent the ranges of $0<g_1/\gamma \le 2$ and $0<g_2/\gamma \le 2$ are plotted to predict EP3. (C, D) Real and imaginary parts of the eigenvalues ${\mathrm{\Omega }}_0$ (see the black lines) and ${\mathrm{\Omega }}_{\pm }$ (see the red and blue lines) versus the same coupling strength $g/\gamma$ in the conditions of Eq. 20.

We noted that when $g_1=g_2$, the pseudo-Hermitian conditions in Eq. 16 reduce to

Moreover, the coefficients in Eq. 18 become

The discriminant in Eq. 19 is $\eta =B^2-4AC=12c_1^3$. We substitute the coefficients in Eq. 21 into Eq. 17 and obtained

Three roots can be obtained by solving the equation

It is obvious that the three real solutions coalesce into one when $c_1=0$, i.e.,

This is EP3 of the proposed QE-cavity QED system. However, when $c_1\ne 0$, two roots of Eq. 17 coalesce to a typical point, $\mathrm{\Omega }={\mathrm{\Omega }}_{\mathrm{E}\mathrm{P}2}$, which means that EP3 transformed into EP2. To verify this result, we plot the real and imaginary parts of the eigenvalues of Eq. 17 in Figures 9C,D, respectively. We noted that the minimum value of the coupling strength $g_{\mathit{min}}=\frac{\gamma }{2}$ is smaller than that of $g_{EP3}=\frac{\gamma }{\sqrt{3}}$. Clearly, when $g_{\mathit{min}}\le g<g_{EP3}$, one eigenvalue ${\mathrm{\Omega }}_0$ is real and the other eigenvalue ${\mathit{\Omega }}_{\pm }$ is a pair of complex conjugate. For $g>g_{EP3}$, all three eigenvalues are real. At a critical point $g=g_{EP3}$, these eigenvalues coalesce to EP3, i.e., ${\mathit{\Omega }}_0={\mathit{\Omega }}_{\pm }={\Delta }_c$.

5 EP3-enhanced strong PB effects in the system

According to the analysis in Section 4, we demonstrate that there is a typical EP3 in the pseudo-Hermiticity condition. In this section, we will study the PB effect at EP3. First, in order to get the minimum value of $g^{\left(2\right)}\left(0\right)$, we should substitute Eq. 20 into Eq. 9, i.e.,

Therefore, we can obtain the conditions for $\left|D_1\right|\to 0$ as

Obviously, when $3g^2-{\gamma }^2=0$ (i.e., $g=g_{EP3}=\frac{\gamma }{\sqrt{3}}$), the corresponding cavity detuning at ${\Delta }_c^0={\Delta }_c^{\left(\pm \right)}=0$. This desired operator regime is found, which result in a strong photon blockade phenomenon, as shown in Figure 10B. In addition, when $g<g_{EP3}$, only detuning ${\Delta }_c^0=0$ is allowed. The maximum value of the mean number photon $I_c$ and the minimum value of $g^{\left(2\right)}\left(0\right)$ can be achieved at this position, as described in Figure 10A. In the case of $g>g_{EP3}$, all of the detunings in Eq. 26 are allowed, as shown in Figure 10C. Therefore, we can observe the PB effect in three regimes with the increase in the coupling strength. The aforementioned results exhibit that under pseudo-Hermiticity conditions, the strong PB phenomenon will be observed at EP3 even if the coupling strengths are weak.

FIGURE 10. Logarithmic plots of $g^{\left(2\right)}\left(0\right)$ and $I_c$ as a function of the cavity frequency detuning ${\Delta }_c/\gamma$ in the cases of (A) $g<g_{EP3}$, (B) $g=g_{EP3}$, and (C) $g>g_{EP3}$. The green dots represent the optimal operator regime for the realization of the strong PB phenomenon at ${\Delta }_c^0$ and ${\Delta }_c^{\left(\pm \right)}$.

We can also find out EP3 in the logarithmic plot of $g^{\left(2\right)}\left(0\right)$ as functions of the normalized detuning ${\Delta }_c/\gamma$ and QE-cavity coupling strength $g/\gamma$, which is shown in Figure 11A. The strong PB phenomenon is obtained at $g=g_{EP3}$, which is shown as the pink area in Figure 11B. In the following section, to explain reasons for the photon blockade in the ${\mathrm{\Omega }}_{\pm }$ region, we give the anharmonic ladder-type energy-level structure in this region (see Figure 12A), where the absorption of the second photon will be blocked, owing to the energy mismatch. This physical mechanism is similar to CPB. On the contrary, in the ${\mathrm{\Omega }}_0$ region, the obvious PB located at the optimal detuning at ${\Delta }_c=0$ comes from the destructive interference between different transition paths (see Figure 12B), which is similar to UPB. At EP3, the PB effect will be significantly enhanced through the coinciding cases of CPB- and UPB-based photon blockade.

FIGURE 11. EP3-enhanced strong PB effect in the pseudo-Hermitian conditions. (A) Plot of $g^{\left(2\right)}\left(0\right)$ as functions of normalized detuning ${\Delta }_c/\gamma$ and the coupling strength $g/\gamma$. EP3 is marked by the red circle. In (B), the pink dashed area indicates the obvious PB effect at $g=g_{EP3}$. We take the parameters as ${\Delta }_c=0,\kappa =-2\gamma$ and ${\delta }_1=-{\delta }_2=\sqrt{g^2-\frac{{\gamma }^2}{4}}$.

FIGURE 12. Explanation of the PB effect of the non-Hermitian system in the pseudo-Hermitian condition for different regions: (A) in the ${\mathrm{\Omega }}_{\pm }$ region and (B) in the ${\mathrm{\Omega }}_0$ region.

6 Comparison with the enhanced PB effect at EP2 and EP3 in the pseudo-Hermitian condition

In practice, the QE-cavity coupling strengths are position-dependent. Hence, the two coupling strengths are usually different for more general cases. For these cases, the transformation between EP3 and EP2 is achieved by changing the coupling strength. As an example, we take $g_2=0.8g_1$, which meet the condition of $\left[g_1^2+g_2^2\right]\ge \frac{{\gamma }^2}{2}$. In Figure 13A, we numerically plot the real parts of eigenvalues of the Hamiltonian (11) as a function of the coupling strength $g_1/\gamma$. It is not difficult to find that there is a typical EP2 at $g_1/\gamma =g_{EP2}/\gamma =0.752$. As expected, the minimum areas in the $g^{\left(2\right)}\left(0\right)$ spectrum (see the blue pattern) are very well in agreement with the real eigenvalues of the effective non-Hermitian Hamiltonian $H_r$ in Eq. 11, as shown by the dashed lines. One of the eigenvalues ${\mathrm{\Omega }}_0$ is always real for arbitrary $g_1/\gamma$ (as the black lines show), and the other two eigenvalues ${\mathrm{\Omega }}_{\pm }$ are a pair of complex conjugates (as the blue and red lines show) when $g_1>g_{EP2}$.

FIGURE 13. In the pseudo-Hermitian conditions when $g_1\ne g_2$, EP3 is transformed into EP2. (A) The minimum value of $g^{\left(2\right)}\left(0\right)$ (see the blue areas) is in good agreement with the real parts of solutions in Eq. 17, where the black dashed line denotes ${\mathrm{\Omega }}_0$ and the red and blue dashed lines denote ${\mathrm{\Omega }}_{\pm }$. (B) Around EP2, the minimum value of $g^{\left(2\right)}\left(0\right)$ and the large mean photon number $I_c$ are achieved simultaneously, indicating the strong PB effect. We take the parameter as ${\Delta }_c/\gamma =0.43,$ and the others are same as those in Figure 11.

Clearly, by comparing the results of Figure 11B and Figure 13B, it is obvious that this new scheme with different coupling strengths exhibits a very strong PB effect at EP2, having approximately two orders of magnitude reduction of $g^{\left(2\right)}\left(0\right)$. However, the mean photon number is almost invariant, which implies that it is independent of the transformation of EPn.

Moreover, considering the two-level QE interaction with the dipole–dipole interaction, we can still find EP2 at $g_{EP2}/\gamma =0.711$ by taking $g_1=g_2=g\mathrm{a}\mathrm{n}\mathrm{d}J=0.2g$ via using the same method explained in Section 4 (Appendix A). By numerically solving Eq. 17, we can obtain the position of EP2, as shown in Figure 14A. Similarly, we can clearly see the obvious PB effect at $g_1=g_{EP2}$ in Figure 14B. By comparing the results in Figures 11B, 14B, the system with DDI shows stronger PB effects at EP2 than the system without DDI.

FIGURE 14. In the pseudo-Hermitian conditions when $g_1=g_2=g$ and $J=0.2g$, there is also an EP2 in the system. (A) Plots of $g^{\left(2\right)}\left(0\right)$ versus the normalized detuning ${\Delta }_c/\gamma$ and the coupling strength $g/\gamma$. In (B), the strong PB effect is shown at $g_{EP2}/\gamma =0.711$. We choose detuning as ${\Delta }_c/\gamma =0.417,$ and other parameters are same as those in Figure 11.

According to aforementioned analysis, we find that the PB effect at EP2 is more obvious than that at EP3 in the condition of the balanced gain–loss rate. When the coupling strength is different, i.e., $g_1\ne g_2$, the photon loss in the paths $\left|\mathrm{1,0,0}\right.〉\stackrel{g_1}{{\displaystyle \to }}\left|\mathrm{0,1,0}\right.〉\stackrel{\epsilon }{{\displaystyle \to }}\left|\mathrm{1,1,0}\right.〉\stackrel{\sqrt{2}{g}_1}{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉$ and $\left|\mathrm{1,0,0}\right.〉\stackrel{g_2}{{\displaystyle \to }}\left|\mathrm{0,0,1}\right.〉\stackrel{\epsilon }{{\displaystyle \to }}\left|\mathrm{1,0,1}\right.〉\stackrel{\sqrt{2}{g}_2}{{\displaystyle \to }}\left|\mathrm{2,0,0}\right.〉$ is asymmetric, which strengthens the quantum interference of the transition paths. In addition, when we consider the influence of DDI between the two emitters, photon loss will emerge in three or more paths. Therefore, the destructive interference between different paths of two-photon excitation will be enhanced, resulting in the more apparent photon blockade effect.

7 Conclusion

In short, we have studied the photon blockade effects in a cavity QED system, where the single-mode cavity is gain and the emitters are loss. Through the analytical solution and numerical results, we, respectively, obtain the equal-time second-order correlation functions to describe the intensity of photon blockade for different cases. We find an interesting phenomenon that there is an EP2 in the system in specific conditions. At this point, the perfect photon antibunching can be observed. Moreover, we find that the physical mechanism of the photon blockade is completely different in PT-symmetric and PT-symmetric-broken regions. For the PT-symmetric region, the anharmonicity of the eigenenergy spectrum occurs, which is similar to CPB. However, in the PT-symmetric-broken region, the interference paths with the photon gain and loss result in UPB. At EP2, the UPB phenomenon is most obvious.

Then, we derive the pseudo-Hermiticity conditions for predicting EP3. The PB effect is also improved at EP3, and we can also explain the photon blockade in different regions. Compared with EP3- and EP2-enhanced PB in different pseudo-Hermiticity conditions, we find that the EP2-enhanced PB may exhibit smaller second-order correlation function. Our work provides a new theoretical foundation for the realization of strong PB effects without strong enough non-linearity of the system under the existing experimental conditions. Our research mainly focused on the theoretical model of photon blockade without experiments. With the development of quantum technologies, we believe that high-quality single-photon sources will be prepared based on EPs in the future.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.

Author contributions

XZ and GZ took the lead on the research work by running the simulations, performing most of the analysis, and producing most of the figures. ZL contributed to the code development. ZL wrote substantial parts of the manuscript. All authors contributed equally to the discussions, read the manuscript, and provided critical feedback.

Acknowledgments

ZL acknowledges the support of the National Natural Science Foundation of China (62075004 and 11804018) and the Beijing Natural Science Foundation (4212051).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References