Next Article in Journal
Thermal Properties of New Insulating Juncus Maritimus Fibrous Mortar Composites/Experimental Results and Analytical Laws
Next Article in Special Issue
Superfluids, Fluctuations and Disorder
Previous Article in Journal
A Comparative Study of TiO2 Paste Preparation Methods Using Solvothermally Synthesised Anatase Nanoparticles in Dye-Sensitised Solar Cells
Previous Article in Special Issue
Topological Qubits as Carriers of Quantum Information in Optics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Photon Blockade in a Coupled Cavity System Mediated by an Atom

1
College of Information School, East China Jiaotong University, Nanchang 330013, China
2
School of Science, Zhejiang Sci-Tech University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(5), 980; https://doi.org/10.3390/app9050980
Submission received: 15 January 2019 / Revised: 2 March 2019 / Accepted: 4 March 2019 / Published: 8 March 2019
(This article belongs to the Special Issue Quantum Optics for Fundamental Quantum Mechanics)

Abstract

:
We investigate theoretically the photon statistics in a coupled cavity system mediated by a two-level atom. The system consists of a linear cavity weakly driven by a continuous laser, and a nonlinear cavity containing an atom inside. We find that there exists a photon blockade in the linear cavity for both parameter regimes where the coupling strength between the atom and the nonlinear cavity is greater (or less) than the dissipation rate of the linear cavity. We also extend our model by pumping the two cavities simultaneously and find that the conventional photon blockade is apparent in the linear cavity, whereas the unconventional photon blockade appears in the nonlinear cavity. These results show that our work has potential applications for a single photon source in a weakly nonlinear system.

1. Introduction

Single photon source and single photon detectors have great potential applications in precision measurement instruments, quantum information computing, and quantum communications [1,2]. For instances, in quantum information science, single photons can be promising quantum bits on which information can be encoded using interferometers or polarization, and can be used in quantum key distribution and other applications. In other domains, single photon detectors are used in a wide range of applications, such as bioluminescence detection, Deoxyribonucleic acid (DNA) sequencing, and picosecond imaging circuit analysis. Single photon emitters deliver photons one at a time, which is also called photon antibunching phenomenon [3]. This single photon is quantized by its energy ν k , where is Planck’s constant, and ν k is the frequency of the k th mode in the quantized electromagnetic field. A single photon source can be generated by different mechanisms. Clauser [4] first generated single photons in 1974 based on a cascade transition of calcium atoms. Diedrich and Walther investigated the antibunching and sub-Poissonian properties of the resonance fluorescence of a single atomic ion stored in a radio-frequency trap [5]. Nano-objects in condensed matter are also possible sources of single photons [6], such as semiconductor nanocrystals [7], color centers in diamonds [8], and quantum dots [9,10,11].
In cavity quantum electrodynamics (QED), a photon blockade is typically achieved by the coupling of nonlinear freedoms that give rise to the anharmonicity of the Jaynes–Cummings ladder of the system eigenstates. A conventional photon blockade is the result of the strong nonlinearity of the system, which requires the single-photon nonlinearity to be at least larger than the mode linewidth κ [12,13,14,15]. A photon blockade created in weakly nonlinear systems, which are far more natural in many areas of photonics and are more feasible to integration and scalability, is referred to as an unconventional photon blockade (UPB) [16,17,18]. The underlying mechanism for the UPB is the destructive quantum interference between different excitation pathways, which usually requires that the system involves several degrees of freedom to establish multi-transition paths. The UPB was first put forward by Liew and Savona [16] in a pair of coupled quantum modes. The UPB has been, since then, studied in various cavity systems, such as two coupled Kerr nonlinear cavity structures with one laser-driven [19], one linear cavity coupled to a Kerr nonlinear cavity with one driving laser [20], two nonlinear cavities with two driving laser [21], two coupled Kerr cavities in a semiconductor [22] and in all-silicon [23], two overlapped cavities corresponding to low frequency and high frequency modes [24], and one cavity coupled to a quantum dot with two driving laser [25].
We know that physical mechanisms of conventional and unconventional photon blockades are different. For conventional photon blockades, the system needs to have a strong nonlinarity effect, whereas UPBs need a weak nonlinearity effect. For obtaining the unconventional photon blockade, we introduce a two-level atom system in Cavity B. We use a pump laser to couple the cavity field, and the interaction between the cavity and the atom can then be effectively manipulated, which provide a method for realization of the unconventional photon blockade. In this paper, we go back to the basic Jaynes–Cummings model and study the conventional and unconventional photon blockade effect in a linear cavity coupled to a nonlinear cavity. The nonlinearity of the second cavity is mediated by a two-level atom. An atom itself is a potential emitter of light, and the cavity can enhance the spontaneous emission rate and the photon production rate. The second cavity provides different transition paths for the photons in the first cavity. We find that there are conventional and unconventional photon blockade effects when the first cavity, the second cavity, and the atom are resonant with the driving laser. So far, the coexistence of two different photon blockades has not been reported. Laser manipulation and atomic coupling are relatively mature technologies, so our scheme is feasible in experiment. A photon blockade is an effective method for realization of a single photon source, which plays an important role in quantum communication and quantum information processing.
The rest of this paper is organized as follows. In Section 2, we introduce the model, the photon transition paths, and the ladder of the dressed states. In Section 3, we demonstrate the conventional and unconventional single photon blockade by numerically solving the master equation and calculating the equal-time second correlation function. In Section 4, we study the optimal parameters for equal-time second correlation function using the non-Hermitian effective Hamiltonian method. In Section 5, we provide some extension discussions of our proposal by adding a simultaneous pump laser to the nonlinear cavity. Conclusions are given in Section 6.

2. Model

As shown in Figure 1, we consider a coupled cavity system that consists of two cascaded cavities. The first cavity (Cavity A) is driven by a classical light field with frequency ω L and amplitude ε . The Hamiltonian of Cavity A with the driven leaser is given by
H A = ω a a a + ε ( a e i ω L t + a e i ω L t )
where a is the annihilation operator of Cavity A whose resonant frequency is ω a , and a is the creation operator. We apply the rotating wave approximation. The Hamiltonian of the second cavity (Cavity B) reads as
H B = ω b b b + ω 1 σ σ + g ( σ b + σ b )
Here, b and b are the annihilation operator and creation operator of the mode in Cavity B, whose resonant frequency is ω b . The lowering (up) operator of the two-level atom in Cavity B is denoted by σ ( σ ), and the transition frequency of the atom is ω 1 . g describes the coupling strength between the atom and Cavity B. These two cavities are coupled through photon hopping interactions with strength of J .
In a rotating frame with respect to H 0 = ω L a a + ω L b b + ω L σ σ , the Hamiltonian of the whole system is given by ( = 1 ):
H = Δ a a a + Δ b b b + Δ 1 σ σ + g ( σ b + σ b ) + J ( a b + a b ) + ε ( a + a ) .
Δ a , Δ b and Δ 1 are, respectively, the frequency detuning of Cavity A, Cavity B, and the atom with respect to the driving field.
In the weak-driving regime, where the driving amplitude ε is far less than the dissipation rate of Cavity A, we consider the total excitation number N 2 in our system. Then the bare states of the system are | 0 , 0 , g , | 0 , 0 , e , | 0 , 1 , g , | 1 , 0 , g , | 0 , 1 , e , | 1 , 0 , e , | 1 , 1 , g , | 0 , 2 , g , and | 2 , 0 , g , with the first number denoting the photon number in Cavity A, the second number denoting the photon number in Cavity B, and the third part denoting the status of the atom in Cavity B. For simplification, we denote these nine bare states by | 0 , | 1 , , | 8 , respectively. The Hamiltonian in Equation (3) without the driving field can be expressed in the following matrix:
[ 0 0 0 0 0 0 0 0 0 0 Δ 1 g 0 0 0 0 0 0 0 g Δ b J 0 0 0 0 0 0 0 J Δ a 0 0 0 0 0 0 0 0 0 Δ b + Δ 1 J 0 g 0 0 0 0 0 J Δ a + Δ 1 g 0 0 0 0 0 0 0 g Δ a + Δ b J J 0 0 0 0 g 0 J 2 Δ b 0 0 0 0 0 0 0 J 0 2 Δ a ]
which is defined based on the bare states mentioned above. Assuming Δ a = Δ b = Δ 1 = Δ , the eigenvalues and eigenstates of the system within the two excitation subspace are given by
E 0 = 0 , | ψ 0 = | 0 , 0 , g
E 1 0 = Δ , | ψ 1 0 = J g 2 + J 2 | 0 , 0 , e + g g 2 + J 2 | 1 , 0 , g
E 1 = g 2 + J 2 + Δ , | ψ 1 = g 2 J 2 + 2 g 2 | 0 , 0 , e 1 2 | 0 , 1 , g + J 2 J 2 + 2 g 2 | 1 , 0 , g
E 1 + = g 2 + J 2 + Δ , | ψ 1 + = g 2 J 2 + 2 g 2 | 0 , 0 , e + 1 2 | 0 , 1 , g + J 2 J 2 + 2 g 2 | 1 , 0 , g
E 2 0 = 2 Δ
E 2 = 1 2 ( 2 2 g 2 + 3 J 2 J 2 ( 16 g 2 + J 2 ) + 4 Δ )
E 2 = 1 2 ( 2 2 g 2 + 3 J 2 + J 2 ( 16 g 2 + J 2 ) + 4 Δ )
E 2 + = 1 2 ( 2 2 g 2 + 3 J 2 J 2 ( 16 g 2 + J 2 ) + 4 Δ )
E 2 + + = 1 2 ( 2 2 g 2 + 3 J 2 + J 2 ( 16 g 2 + J 2 ) + 4 Δ ) .
The eigenstates corresponding to the eigenvalues E 2 0 , E 2 , E 2 , E 2 + , and E 2 + + are too cumbersome. We will not list them here. The transition channels are demonstrated in Figure 2a. The multiple excitation paths from state | 0 , 0 , g to state | 2 , 0 , g are the fundamental and essential factors to the single photon blockade in our system. A single photon blockade is apparent when one can find the proper quantum destructive interference conditions among different transition paths from | 1 , 0 , g to | 2 , 0 , g . The eigenenergy levels are shown in Figure 2b. The blue arrows and red arrows illustrate the anharmonicity of the Jaynes–Cummings ladder. When the resonant absorption frequency ω is detuned to reach the state | ψ 1 (corresponding to the eigenvalues E 1 ) or | ψ 1 + , it will block the absorption of a second photon at frequency ω because the states corresponding to eigenenergy E 2 , E 2 , E 2 0 , E 2 + , and E 2 + + are all detuned from ω .

3. Numerical Computation Conventional and Unconventional Photon Blockade

In this section, we study the conventional and unconventional photon blockade effects by numerically calculating the equal-time second-order correlation function of Cavity A.

3.1. Quantum Master Equation

We assume our system satisfies the Born approximation and is in a “short-memory environment” such that it also satisfies the Markov approximation. The time-evolution of the system density matrix is then governed by the Lindblad master equation [26,27]:
ρ = i [ ρ , H ] + κ a 2 ( n ¯ t h a + 1 ) [ a ] + κ a 2 n ¯ t h a [ a ] + κ b 2 ( n ¯ t h b + 1 ) [ b ] + κ b 2 n ¯ t h b [ b ] + γ 2 ( n ¯ t h 1 + 1 ) [ σ ] + γ 2 n ¯ t h 1 [ σ ]
The super operator [ o ] = 2 o ρ o ρ o o o o ρ is the Lindblad term [28,29] accounting for the losses to the environment. n t h x = [ exp ( ω x / K B T x ) 1 ] 1 , x = { a , b , 1 } are the average thermal excitation numbers of the bath at temperature T x , with K B as the Boltzmann constant. κ a and κ b are the dissipation rates of Cavity A and Cavity B, respectively. The transition rate from | e to | g of the atom is described by the term proportional to γ 2 ( n ¯ t h 1 + 1 ) , which contains a rate for spontaneous transitions and a rate for stimulated transitions induced by thermal photons. The last term in Equation (14) describes the transition rate from | g to | e , which is obtained by absorbing thermal photons from the cavity field. At optical frequencies and laboratory temperatures, the thermal photon number n ¯ t h x can be completely negligible [26].
The time evolution of the mean photon numbers in Cavity A can be obtained by solving the master Equation (14). The result is shown in Figure 3 in which Cavity A, Cavity B, and the atom have the same detuning with respect to the driving laser.
The photon-number statistics in Cavity A can be evaluated by the result of the photon correlation function, which is defined in the steady-time limit as
g ( 2 ) ( τ ) = lim t > a ( t ) a ( t + τ ) a ( t + τ ) a ( t ) a ( t ) a ( t ) 2 .
We assume that there is zero delay between photons, i.e., τ = 0 , and Equation (15) can then be simplified to
g ( 2 ) ( 0 ) = a a a a s s a a s s 2
where ‘ss’ means steady state. In the weak excitation limit, g ( 2 ) ( 0 ) can be roughly estimated by the system density matrix element expanded in the bare state subspace. That is,
a a a a s s = T r ( a a a a ρ s s ) k = 0 8 k | a a a a m , n = 0 8 ρ m n | m n | k = 2 ρ 88
a a s s = T r ( a a ρ s s ) k = 0 8 k | a a m , n = 0 8 ρ m n | m n | k = ρ 33 + ρ 55 + ρ 66 + 2 ρ 88 ρ 33
where k in | k is the index corresponding to the bare states listed in Section 2.
The detuning with respect to the driving light is an important factor affecting the photon blockade effect. We plot the equal-time second-order correlation function g ( 2 ) ( 0 ) as a function of the scaled detuning Δ / κ a in Figure 4, where Δ = Δ a = Δ b = Δ 1 . The red dashed curve in Figure 4 is calculated according to Equations (16)–(18), and the blue curve is obtained by solving Equation (14) numerically. These two curves match well in the parameter domain. Thus, we can analyze the peaks and dips of the second correlation curve based on the profiles of the ρ 33 and ρ 88 curves. As shown in Figure 4, ρ 33 ρ 88 with ρ 33 = 1 , 0 , g | ρ s s | 1 , 0 , g and ρ 88 = 2 , 0 , g | ρ s s | 2 , 0 , g . There are three peaks in the ρ 33 curve that are caused by the resonance transitions of the system. As shown in Figure 2b, there are three resonant transition channels for a single photon, i.e., | ψ 0 | ψ 1 , | ψ 0 | ψ 1 0 , and | ψ 0 | ψ 1 + . The positions of these three peaks can be calculated from the eigenvalues of these states. The positions of the left peak, the middle one, and the right one are obtained based on the resonant transition conditions in the rotating frame E 1 + E 0 = 0 , E 1 0 E 0 = 0 , and E 1 E 0 = 0 . According to Equations (5)–(8) and based on parameters used in Figure 4, one can obtain the positions of these three peaks: −9.1782, 0, and 9.1782. The five peaks in the ρ 88 curve are caused by the resonant transition of two photons, whose position can be calculated in the same way as those in the ρ 33 curve, and the values are −6.6446, −4.2683, 0, 4.2683, and 6.6446, respectively. The other two peaks in the ρ 88 curve, whose positions are the same as the left one and right one in the ρ 33 curve, are induced by a single photon resonant transition. The dips in the ρ 33 and ρ 88 curves are induced by the destructive quantum interference effect between different transition paths, as shown in Figure 2b. The locations of the two peaks in the correlation function g ( 2 ) ( 0 ) correspond to the two dips in the ρ 33 curve, where the excitation of the single photon resonant transition is low. There is a single photon resonant transition and a two-photon resonant transition at Δ / κ a = 0 in the ρ 33 and ρ 88 curves. The dip at Δ / κ a = 0 in the correlation function g ( 2 ) ( 0 ) curve means the single photon resonant transition is far stronger than the two-photon resonant transition at that position. The locations of the other dips in the g ( 2 ) ( 0 ) curve correspond to those in the ρ 88 curve.
To investigate the influence of the coupling strength g and the hopping interaction strength J on the correlation function, in Figure 5, we plot the equal-time second correlation function as the two scaled parameters, g / κ a and J / κ a , with the fixed detuning Δ = 0 , which is a detuning of the single photon resonant transition. From Figure 5, we can see that, when the coupling strength between Cavity B and the atom is larger than κ a , i.e., g / κ a > 1 , one can obtain a minimum value of g ( 2 ) ( 0 ) at the ratio between the hopping interaction strength and the coupling strength J / g 1.25 . At this ratio, the two photon transition is suppressed and the single photon transition is far stronger. As illustrated in Figure 5, two-photon antibunching can also be obtained when the coupling strength between Cavity B and the atom is weaker than the decay rate of Cavity A. That is, the single photon blockade effect can be obtained in our system under weak nonlinearity, which is induced by the atom in Cavity B in our system with g < κ a . Around the other region at g 0.5 κ a , one can obtain a low value of g ( 2 ) ( 0 ) as described in Figure 5. However, if we continue to decrease the coupling strength below 0.3 κ a , the photon blockade effect will disappear, and a two-photon bunching effect will become obvious.

3.2. Unconventional Photon Blockade with Atom-Cavity Detuning

In this section, we focus on the unconventional photon blockade when Cavity A, Cavity B, and the atom are out-of resonance based on weak nonlinearities. By solving Equation (14) numerically, we plot the equal-time second order correlation function log 10 g ( 2 ) ( 0 ) as a function of the scaled cavity detuning Δ / κ a and the scaled atom detuning Δ 1 / κ a in Figure 6, where Cavity B is resonant with Cavity A. Based on the parameters of Figure 6, when the atom is weakly coupled to Cavity B by g = 0.67 κ a and the interaction strength between the two cavities is J = 2 κ a , one can obtain an irregular single photon blockade region at the center of the figure while the detuning of Cavity A, Cavity B, and the atom are all within [ 2 κ a , + 2 κ a ] . By increasing the coupling strength and the hopping strength to g = 0.8 κ a and J = 4 κ a , we can obtain the other two symmetric photon block detuning regimes located at about [ κ a , 3 κ a ] and [ + κ a , + 3 κ a ] , as shown in Figure 6b. With the fixed coupling strength g = 0.67 κ a and based on the condition that the atom in Cavity B is resonant with the driving laser, i.e., Δ 1 = 0 , Figure 7 shows the correlation function log 10 g ( 2 ) ( 0 ) as a function of different detunings of Cavities A and B for photon hopping strength between the two cavities J = 3 κ a in Figure 7a and J = 6 κ a in Figure 7b. One can see from Figure 7 that the single photon blocked region ( log 10 g ( 2 ) ( 0 ) < 0) is symmetrical and centered on Δ a = 0 and Δ b = 0 . As the moderate increase of the hopping strength J , the single photon blocked region stretches to both the red sideband and blue sideband of Cavity A.

4. Optimal Parameters for Sub-Poissonian Characters Derived by the Non-Hermitian Effective Hamiltonian Method

In this section, we focus on the analytical optimal conditions for the system parameters that maximize the sub-Poissonian character of Cavity A. In the weak driving limit, the system states can be expressed as an expansion on the bare states mentioned in Section 2:
| ψ ( t ) = c 000 | 00 g + c 001 | 00 e + c 010 | 01 g + c 100 | 10 g + c 011 | 01 e + c 101 | 10 e + c 110 | 11 g + c 020 | 02 g + c 200 | 20 g
The system is governed by a stochastic Schrodinger equation as (assuming = 1 )
i d d t | ψ ( t ) = H e f f | ψ ( t )
where the non-Hermitian effective Hamiltonian [17,18,26] can be written as
H e f f = ( Δ a i κ a 2 ) a a + ( Δ b i κ b 2 ) b b + ( Δ 1 i γ 2 ) σ σ + g ( σ b + σ b ) + J ( a b + a b ) + ε ( a + a )
Here, we consider the system evolution under the condition that the three components Cavity A, Cavity B, and the atom are resonant, i.e., Δ a = Δ b = Δ 1 = Δ , and the bandwidth of Cavity A equals that of Cavity B, i.e., κ a = κ b = κ . The time-dependent coefficients c m n p ( t ) can be obtained from Equations (20) and (21). According to the idea of performing a consistent expansion of the elements to dominant order in powers of ε / κ a [30], we neglect the subleading order in the driving laser amplitude and obtain the equations of c m n p ( t ) :
i c 001 = ( Δ i γ 2 ) c 001 + g c 010
i c 010 = ( Δ i κ 2 ) c 010 + g c 001 + J c 100
i c 100 = ( Δ i κ 2 ) c 100 + J c 010 + ε c 000
i c 011 = ( Δ i κ 2 ) c 011 + ( Δ i γ 2 ) c 011 + 2 g c 020 + J c 101
i c 101 = ( Δ i κ 2 ) c 101 + ( Δ i γ 2 ) c 101 + g c 110 + J c 011 + ε c 001
i c 110 = ( 2 Δ i κ ) c 110 + g c 101 + 2 J c 020 + 2 J c 200 + ε c 010
i c 020 = ( 2 Δ i κ ) c 020 + 2 g c 011 + 2 J c 110
i c 200 = ( 2 Δ i κ ) c 200 + 2 J c 110 + 2 ε c 100 .
The relationship between different c m n p ( t ) in Equations (22)–(29) once again illustrates the multi-transition channels between different system states, which is consistent with that in Figure 2a. The amplitudes of the coefficient c m n p ( t ) in Equation (19) can be divided into three levels with respect to ε / κ a . That is, the amplitudes of c 000 ( t ) is of the zero order of ε / κ a , { c 001 ( t ) , c 010 ( t ) , and c 100 ( t ) } of the first order and { c 011 ( t ) , c 101 ( t ) , c 110 ( t ) , c 020 ( t ) , c 200 ( t ) } of the second order, i.e., c 000 ( t ) c 001 ( t ) , c 010 ( t ) , c 100 ( t ) c 011 ( t ) , c 101 ( t ) , c 110 ( t ) , c 020 ( t ) , c 200 ( t ) . The average photon occupations in Cavity A and the equal-time second-order correlations then approximate to
n a = a a = | c 100 | 2 + | c 101 | 2 + | c 110 | 2 + 2 | c 200 | 2 | c 100 | 2
g ( 2 ) ( 0 ) = a a a a n a 2 2 | c 200 | 2 | c 100 | 4 .
The expressions of c 100 and c 200 can be obtained by solving Equations (22)–(29) for the steady state | ψ s s with c m n p = 0 and c 000 = 1 .
c 100 = 2 [ 4 i g 2 + ( γ + 2 i Δ ) ( i κ + 2 Δ ) ] ε 4 g 2 ( κ + 2 i Δ ) + [ 4 J 2 + ( κ + 2 i Δ ) 2 ( γ + 2 i Δ ) ]
c 200 = 2 2 { 32 g 6 + 2 g 2 X 1 2 ( 16 J 2 + X 2 ) 8 g 4 ( 4 J 2 X 1 X 3 ) + X 1 3 X 5 } ε 2 M 1
M = 32 g 6 X 1 2 + 2 g 2 ( 4 J 2 + X 1 2 ) ( X 1 2 X 2 + 4 J 2 X 4 ) + 8 g 4 X 1 [ X 1 2 X 3 + J 2 ( 8 γ 4 κ + 8 i Δ ) ] + X 5 [ 4 J 2 + X 1 2 ] 2 X 1
where, X 1 = κ + 2 i γ Δ , X 2 = 2 κ 2 + 7 κ γ + 5 γ 2 + 22 i κ Δ + 34 i γ Δ 56 Δ 2 , X 3 = 3 κ + 4 γ + 14 i Δ , X 4 = 2 κ 2 κ γ + γ 2 + 6 i κ Δ + 2 i γ Δ 8 Δ 2 and X 5 = [ 4 J 2 + ( κ + γ + 4 i Δ ) 2 ] ( γ + 2 i Δ ) .
The optimal parameter for single photon blockade in Cavity A can be obtained by set c 200 = 0 under the condition c 100 0 . Based on Equation (33), one can obtain the optimal parameter pairs of the coupling strength g and the detuning Δ with fixed hopping strength J and atom decay rate γ . When J is fixed at 2.5 κ a and γ = 0.6 κ a , we can obtain the real solution of the optimal parameter pairs {{ g 7.63415 κ a , Δ 6.89979 κ a }, { g 7.63415 κ a , Δ 6.89979 κ a }, { g 0.924508 κ a , Δ 1.32075 κ a }, { g 0.924508 κ a , Δ 1.32075 κ a }, { g 0.924508 κ a , Δ 1.32075 κ a }, { g 0.924508 κ a , Δ 1.32075 κ a }, { g 7.63415 κ a , Δ 6.89979 κ a }, { g 7.63415 κ a , Δ 6.89979 κ a }}.
We verify the solution by plotting the equal-time second-order correlation function g ( 2 ) ( 0 ) as a function of the scaled detuning for different coupling strength g in Figure 8. As expected, the values of g ( 2 ) ( 0 ) are much smaller when we use the { g , Δ } pairs listed above than others, i.e., it shows two strong photon antibunching effects at Δ 1.3 κ a with g = 0.924508 κ a in Figure 8a and Δ 6.9 κ a with g = 7.63415 κ a in Figure 8b. We also see the symmetric property from optimal { g , Δ } solutions, that is, when the coupling strength g is fixed at one of the optimal values, there are two symmetric optimal detunings for Cavity A, Cavity B, and the atom with respect to the driving laser, and vice versa. The underlying physical mechanism of the optimal solution for a single photon blockade lies in the destructive interference between different paths of two photon excitation in Cavity A. As illustrated in Figure 2a and Equation (29), the occupation of | 2 , 0 , g can be minimized by reducing the contribution from | 1 , 0 , g and | 1 , 1 , g . When destructive interference between the direct excitation path | 0 , 0 , g | 1 , 0 , g | 2 , 0 , g and the indirect excitation paths | 0 , 0 , g | 1 , 0 , g | 0 , 1 , g | 0 , 0 , e ( | 1 , 0 , e | 0 , 1 , e | 0 , 2 , g | 1 , 1 , g ) | 2 , 0 , g occurs, the occupation of | 2 , 0 , 0 will be reduced; hence, the value of g ( 2 ) ( 0 ) is decreased.
Figure 9 shows several optimal { g , Δ } pairs for different atom decay rate γ . When γ = 0.25 κ a , one can get { { g ± 0.782092 κ a ,   Δ ± 1.19804 κ a } , { g ± 7.89333 κ a ,   Δ ± 7.17599 κ a } }, and { { g ± 0.817538 κ a ,   Δ ± 1.22938 κ a } , { g ± 7.81893 κ a ,   Δ ± 7.09409 κ a } } with γ = 0.35 κ a , { { g ± 0.857111 κ a ,   Δ ± 1.26361 κ a } , { g ± 7.74537 κ a ,   Δ ± 7.01522 κ a } } with γ = 0.45 κ a . As shown in Figure 9, with the increase in the atom lifetime, the optimal coupling strength g decreases in the UPB regime where g < κ a and increases in the conventional photon blockade regime where g > κ a .

5. Supplementary Discussion

In the above sections, we find that Cavity A is weakly driven by a continuous laser field and Cavity B is not. If the cavities in our system are both driven by a laser field, the Hamiltonian of the system reads
H = Δ a a a + Δ b b b + Δ 1 σ σ + g ( σ b + σ b ) + J ( a b + a b ) + ε a ( a + a ) + ε b ( b + b ) .
There will then also be direct and indirect excitation paths from state | 0 , 0 , g to | 0 , 2 , g , which provides a fundamental condition for the photon blockade in Cavity B.
We plot the equal-time second-order correlation function in Figure 10 as a function of the two scaled parameters g / κ a and J / κ a under the condition that the two driving lasers are of the same frequency. With zero detuning, i.e., Cavity A, Cavity B, and the atom are resonant with the pump laser; we can see from Figure 10 that there is a single photon blockade in both Cavity A and Cavity B. However, the parameter regimes for the photon blockade are different in Cavity A and Cavity B. As shown in Figure 10a, the two-photon antibunching effect is apparent in Cavity A when g κ a and the ratio between J and g is approximately equal to one for a minimum g ( 2 ) ( 0 ) , whereas the photon blockade effect in Cavity B appears when g < κ a ( κ b ) with J > g . We obtain two different types of photon blockades at the same time: a conventional photon blockade is apparent in the linear cavity A and an unconventional photon blockade is dominant in the nonlinear cavity B, here Cavity A, Cavity B, and the atom in Cavity B are all resonant with the pump laser.

6. Conclusions

In summary, we have studied the photon statistics of a coupled cavity QED system that consists of a linear cavity and a nonlinear cavity that consists of a two-level atom. The linear cavity is weakly driven by a continuous laser field, and the two cavities are coupled by hopping strength J . The single photon blockade effect can be observed in the linear cavity, both conventional and unconventional. Compared with existing schemes, we can achieve a conversion between a conventional photon blockade and an unconventional photon blockade in the same system; more importantly, we can also achieve the coexistence of two different photon blockades. In addition, we obtain the coupling between two cavities via photon hopping interactions, which can be realized by current experimental techniques. The photon blockade provides an effective method for obtaining a single photon source, which is the basis of realizing optical quantum information technology. Thus, our work has potential applications in quantum information processing and quantum communication. Being a theoretical work, we provide detailed theoretical derivation and numerical analysis. In our model, these photon blockade effects are mainly achieved by the quantum destructive interference between different excitation paths. We have analyzed the values of the equal-time second-order correlation function in the truncated Hilbert space by solving the master equation of the coupled cavities. In the same parameter regime, we also derived the optimal atom-cavity coupling strength and cavity detuning parameter pairs, which once again verify the symmetric photon blockade regime for the linear cavity. We also extend our model by adding a pump laser to Cavity B. We find that, under the condition of driving these two cavities simultaneously, the conventional photon blockade in Cavity A is dominant, while in Cavity B the unconventional photon blockade is apparent. In our scheme, the introduction of a two-level atom provides an efficient method for the realization of an unconventional photon blockade in Cavity B. When Cavity B is manipulated by the pump laser, the interaction between the cavity and atom can be effectively controlled, which may affect the nonlinearity of the total system. Thus, we obtain the unconventional photon blockade in the case of weak nonlinearity. Here the pump laser and the two-level atom are well controlled experimentally, and the unconventional photon blockade is not obviously dependent on the decay rate of cavity. Therefore, we provide a feasible scheme for the realization of an unconventional photon blockade.

Author Contributions

Conceptualization, M.-C.L. and A.-X.C.; methodology, A.-X.C.; software, M.-C.L.; formal analysis, M.-C.L.; writing—original draft preparation, M.-C.L.; writing—review and editing, A.-X.C.; supervision, A.-X.C.

Funding

This work was supported by the project of National Natural Science Foundation of China [grant number 11775190] and by Science and Technology of Jiangxi province [grant number 20171BAB212006]; the project of science and technology awarded by the Education Bureau of Jiangxi Province [grant number JJ160503] and by the Science Foundation of Zhejiang SCI-TECH University (17062071-Y).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hadfield, R.H. Single-photon detectors for optical quantum information applications. Nat. Photonics 2009, 3, 696–705. [Google Scholar] [CrossRef]
  2. Eisaman, M.D.; Fan, J.; Migdall, A.; Polyakov, S.V. Invited review article: Single-photon sources and detectors. Rev. Sci. Instrum. 2011, 82, 071101. [Google Scholar] [CrossRef] [PubMed]
  3. Lounis, B.; Orrit, M. Single-photon sources. Rep. Prog. Phys. 2005, 68, 1129–1179. [Google Scholar] [CrossRef]
  4. Clauser, J.F. Experimental distinction between the quantum and classical field-theoretic predictions for the photoelectric effect. Phys. Rev. D 1974, 9, 853. [Google Scholar] [CrossRef]
  5. Diedrich, F.; Walther, H. Nonclassical radiation of a single stored ion. Phys. Rev. Lett. 1987, 58, 203–206. [Google Scholar] [CrossRef] [PubMed]
  6. De Martini, F.; Di Giuseppe, G.; Marrocco, M. Single-mode generation of quantum photon states by excited single molecules in a microcavity trap. Phys. Rev. Lett. 1996, 76, 900–903. [Google Scholar] [CrossRef] [PubMed]
  7. Lounis, B.; Moerner, W.E. Single photons on demand from a single molecule at room temperature. Nature 2000, 407, 491–493. [Google Scholar] [CrossRef] [Green Version]
  8. Brouri, R.; Beveratos, A.; Poizat, J.P.; Grangier, P. Photon antibunching in the fluorescence of individual color centers in diamond. Opt. Lett. 2000, 25, 1294–1296. [Google Scholar] [CrossRef]
  9. Michler, P.; Imamoğlu, A.; Mason, M.D.; Carson, P.J.; Strouse, G.F.; Buratto, S.K. Quantum correlation among photons from a single quantum dot at room temperature. Nature 2000, 406, 968–970. [Google Scholar] [CrossRef]
  10. He, Y.M.; He, Y.; Wei, Y.J.; Wu, D.; Atatüre, M.; Schneider, C.; Höfling, S.; Kamp, M.; Lu, C.Y.; Pan, J.W. On-demand semiconductor single-photon source with near-unity indistinguishability. Nat. Nanotechnol. 2013, 8, 213–217. [Google Scholar] [CrossRef] [Green Version]
  11. Gschrey, M.; Thoma, A.; Schnauber, P.; Seifried, M.; Schmidt, R.; Wohlfeil, B.; Krüger, L.; Schulze, J.H.; Heindel, T.; Burger, S.; et al. Highly indistinguishable photons from deterministic quantum-dot microlenses utilizing three-dimensional in situ electron-beam lithography. Nat. Commun. 2015, 6, 7662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Imamoḡlu, A.; Schmidt, H.; Woods, G.; Deutsch, M. Strongly interacting photons in a nonlinear cavity. Phys. Rev. Lett. 1997, 79, 1467. [Google Scholar] [CrossRef]
  13. Leoński, W.; Tanaś, R. Possibility of producing the one-photon state in a kicked cavity with a nonlinear Kerr medium. Phys. Rev. A 1994, 49, R20–R23. [Google Scholar] [CrossRef] [PubMed]
  14. Birnbaum, K.M.; Boca, A.; Miller, R.; Boozer, A.D.; Northup, T.E.; Kimble, H.J. Photon blockade in an optical cavity with one trapped atom. Nature 2005, 436, 87–90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Peyronel, T.; Firstenberg, O.; Liang, Q.Y.; Hofferberth, S.; Gorshkov, A.V.; Pohl, T.; Lukin, M.D.; Vuletić, V. Quantum nonlinear optics with single photons enabled by strongly interacting atoms. Nature 2012, 488, 57–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Liew, T.C.H.; Savona, V. Single photons from coupled quantum modes. Phys. Rev. Lett. 2010, 104, 183601. [Google Scholar] [CrossRef] [PubMed]
  17. Flayac, H.; Savona, V. Unconventional photon blockade. Phys. Rev. A 2017, 96, 053810. [Google Scholar] [CrossRef]
  18. Bin, Q.; Lü, X.Y.; Bin, S.W.; Wu, Y. Two-photon blockade in a cascaded cavity-quantum-electrodynamics system. Phys. Rev. A 2018, 98, 043858. [Google Scholar] [CrossRef]
  19. Bamba, M.; Imamoğlu, A.; Carusotto, I.; Ciuti, C. Origin of strong photon antibunching in weakly nonlinear photonic molecules. Phys. Rev. A 2011, 83, 021802. [Google Scholar] [CrossRef]
  20. Zou, F.; Lai, D.G.; Liao, J.Q. Photon blockade effect in a coupled cavity system. arXiv, 2018; arXiv:1803.06642. [Google Scholar]
  21. Xu, X.W.; Li, Y. Tunable photon statistics in weakly nonlinear photonic molecules. Phys. Rev. A 2014, 90, 043822. [Google Scholar] [CrossRef]
  22. Shen, H.Z.; Zhou, Y.H.; Yi, X.X. Tunable photon blockade in coupled semiconductor cavities. Phys. Rev. A 2015, 91, 063808. [Google Scholar] [CrossRef]
  23. Flayac, H.; Gerace, D.; Savona, V. An all-silicon single-photon source by unconventional photon blockade. Sci. Rep. 2015, 5, 11223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Gerace, D.; Savona, V. Unconventional photon blockade in doubly resonant microcavities with second-order nonlinearity. Phys. Rev. A 2014, 89, 031803. [Google Scholar] [CrossRef]
  25. Tang, J.; Geng, W.; Xu, X. Quantum interference induced photon blockade in a coupled single quantum dot-cavity system. Sci. Rep. 2015, 5, 9252. [Google Scholar] [CrossRef] [PubMed]
  26. Carmichael, H. Master equations and source. In An Open Systems Approach to Quantum Optics; Springer: Berlin/Heidelberg, Germany, 1993; pp. 22–35. [Google Scholar]
  27. Breuer, H.P.; Petruccione, F. Quantum master equations. In The Theory of Open Quantum Systems; Oxford University Press: Oxford, UK, 2002; pp. 109–216. [Google Scholar]
  28. Carmichael, H.J. Quantum trajectory theory for cascaded open systems. Phys. Rev. Lett. 1993, 70, 2273–2276. [Google Scholar] [CrossRef] [PubMed]
  29. Johansson, J.R.; Nation, P.D.; Nori, F. QuTiP: An open-source Python framework for the dynamics of open quantum systems. Comput. Phys. Commun. 2012, 183, 1760–1772. [Google Scholar] [CrossRef] [Green Version]
  30. Carmichael, H.J.; Brecha, R.J.; Rice, P.R. Quantum interference and collapse of the wavefunction in cavity QED. Opt. Commun. 1991, 82, 73–79. [Google Scholar] [CrossRef]
Figure 1. Schematic of the studied system. The system consists of two cavities, which are coupled to each other through a photon-hopping interaction. The first cavity is a linear cavity driven by a continuous laser. The second cavity is nonlinear with a two-level atom.
Figure 1. Schematic of the studied system. The system consists of two cavities, which are coupled to each other through a photon-hopping interaction. The first cavity is a linear cavity driven by a continuous laser. The second cavity is nonlinear with a two-level atom.
Applsci 09 00980 g001
Figure 2. (a) Energy levels within two excitation subspaces and the corresponding excitation paths; (b) the eigenenergy spectrum of the dressed states within two excitation subspaces.
Figure 2. (a) Energy levels within two excitation subspaces and the corresponding excitation paths; (b) the eigenenergy spectrum of the dressed states within two excitation subspaces.
Applsci 09 00980 g002
Figure 3. Time evolution of the average photon numbers in Cavity A in unit of 2 π / κ a . The parameters are κ b = κ a , γ = 0.6 κ a , Δ a = 7 κ a , Δ a = Δ b = Δ 1 , ε = 0.5 κ a , J = 9 κ a , g = 2 κ a , and n ¯ t h a = n ¯ t h b = n ¯ t h 1 = 0 .
Figure 3. Time evolution of the average photon numbers in Cavity A in unit of 2 π / κ a . The parameters are κ b = κ a , γ = 0.6 κ a , Δ a = 7 κ a , Δ a = Δ b = Δ 1 , ε = 0.5 κ a , J = 9 κ a , g = 2 κ a , and n ¯ t h a = n ¯ t h b = n ¯ t h 1 = 0 .
Applsci 09 00980 g003
Figure 4. Plots of the density matrix element ρ 33 , ρ 88 and the second-order correlation function g ( 2 ) ( 0 ) as functions of the scaled detuning of Cavity A for ε = 0.1 κ a and g = 1.8 κ a . Other parameters are the same as in Figure 3.
Figure 4. Plots of the density matrix element ρ 33 , ρ 88 and the second-order correlation function g ( 2 ) ( 0 ) as functions of the scaled detuning of Cavity A for ε = 0.1 κ a and g = 1.8 κ a . Other parameters are the same as in Figure 3.
Applsci 09 00980 g004
Figure 5. The equal-time second-order correlation function g ( 2 ) ( 0 ) as a function of the scaled coupling strength g / κ a and the scaled hopping interaction strength J / κ a under the resonant condition Δ = Δ a = Δ b = Δ 1 = 0 . Other parameters are the same as in Figure 4.
Figure 5. The equal-time second-order correlation function g ( 2 ) ( 0 ) as a function of the scaled coupling strength g / κ a and the scaled hopping interaction strength J / κ a under the resonant condition Δ = Δ a = Δ b = Δ 1 = 0 . Other parameters are the same as in Figure 4.
Applsci 09 00980 g005
Figure 6. The equal-time second-order correlation function log 10 g ( 2 ) ( 0 ) as a function of the scaled detuning Δ / κ a = Δ a / κ a = Δ b / κ a and Δ 1 / κ a for (a) g = 0.67 κ a and J = 2 κ a ; (b) g = 0.8 κ a and J = 4 κ a . Other parameters are the same as in Figure 4.
Figure 6. The equal-time second-order correlation function log 10 g ( 2 ) ( 0 ) as a function of the scaled detuning Δ / κ a = Δ a / κ a = Δ b / κ a and Δ 1 / κ a for (a) g = 0.67 κ a and J = 2 κ a ; (b) g = 0.8 κ a and J = 4 κ a . Other parameters are the same as in Figure 4.
Applsci 09 00980 g006
Figure 7. The equal-time second-order correlation function log 10 g ( 2 ) ( 0 ) as a function of scaled detuning Δ a / κ a and Δ b / κ b under the condition of the atom resonant with the driving laser Δ 1 = 0 for (a) J = 3 κ a and (b) J = 6 κ a . Other parameters are the same as in Figure 6a.
Figure 7. The equal-time second-order correlation function log 10 g ( 2 ) ( 0 ) as a function of scaled detuning Δ a / κ a and Δ b / κ b under the condition of the atom resonant with the driving laser Δ 1 = 0 for (a) J = 3 κ a and (b) J = 6 κ a . Other parameters are the same as in Figure 6a.
Applsci 09 00980 g007
Figure 8. Plots of equal-time second-order correlation function g ( 2 ) ( 0 ) as a function of the scaled detuning for (a) g = 0.6 κ a , g = 0.8 κ a , and g = 0.924508 κ a ; (b) g = 7 κ a , g = 7.63415 κ a , and g = 8 κ a with fixed J = 2.5 κ a . Other parameters are the same as in Figure 4.
Figure 8. Plots of equal-time second-order correlation function g ( 2 ) ( 0 ) as a function of the scaled detuning for (a) g = 0.6 κ a , g = 0.8 κ a , and g = 0.924508 κ a ; (b) g = 7 κ a , g = 7.63415 κ a , and g = 8 κ a with fixed J = 2.5 κ a . Other parameters are the same as in Figure 4.
Applsci 09 00980 g008
Figure 9. Plots of g ( 2 ) ( 0 ) versus the scaled detuning Δ / κ a for (a) g < κ a ; (b) g > κ a with γ = 0.25 κ a , γ = 0.35 κ a , and γ = 0.45 κ a . Other parameters are the same as in Figure 8.
Figure 9. Plots of g ( 2 ) ( 0 ) versus the scaled detuning Δ / κ a for (a) g < κ a ; (b) g > κ a with γ = 0.25 κ a , γ = 0.35 κ a , and γ = 0.45 κ a . Other parameters are the same as in Figure 8.
Applsci 09 00980 g009
Figure 10. Plots of equal-time second-order correlation function g ( 2 ) ( 0 ) of (a) Cavity A and (b) Cavity B as a function of scaled atom-cavity coupling strength g / κ a and cavity interaction strength J / κ a for κ b = κ a , γ = 0.6 κ a , ε a = ε b = 0.1 κ a , and Δ a = Δ b = Δ 1 = 0 .
Figure 10. Plots of equal-time second-order correlation function g ( 2 ) ( 0 ) of (a) Cavity A and (b) Cavity B as a function of scaled atom-cavity coupling strength g / κ a and cavity interaction strength J / κ a for κ b = κ a , γ = 0.6 κ a , ε a = ε b = 0.1 κ a , and Δ a = Δ b = Δ 1 = 0 .
Applsci 09 00980 g010

Share and Cite

MDPI and ACS Style

Li, M.-C.; Chen, A.-X. A Photon Blockade in a Coupled Cavity System Mediated by an Atom. Appl. Sci. 2019, 9, 980. https://doi.org/10.3390/app9050980

AMA Style

Li M-C, Chen A-X. A Photon Blockade in a Coupled Cavity System Mediated by an Atom. Applied Sciences. 2019; 9(5):980. https://doi.org/10.3390/app9050980

Chicago/Turabian Style

Li, Ming-Cui, and Ai-Xi Chen. 2019. "A Photon Blockade in a Coupled Cavity System Mediated by an Atom" Applied Sciences 9, no. 5: 980. https://doi.org/10.3390/app9050980

APA Style

Li, M. -C., & Chen, A. -X. (2019). A Photon Blockade in a Coupled Cavity System Mediated by an Atom. Applied Sciences, 9(5), 980. https://doi.org/10.3390/app9050980

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop