Abstract
In the last two decades, considerable advances have been made in the investigation of nano-photonics in photonic crystals. Previous theoretical investigations of photon dynamics were carried out at zero temperature. Here, we investigate micro/nano cavity photonics in photonic crystals at finite temperature. Due to photonic-band-gap-induced localized long-lived photon dynamics, we discover that cavity photons in photonic crystals do not obey Bose-Einstein statistical distribution. Within the photonic band gap and in the vicinity of the band edge, cavity photons combine the long-lived non-Markovain dynamics with thermal fluctuations together to form photon states that memorize the initial cavity state information. As a result, Bose-Einstein distribution is completely broken down in these regimes, even if the thermal energy is larger or much larger than the cavity detuning energy. In this investigation, a crossover phenomenon from equilibrium to nonequilibrium steady states is also revealed.
Similar content being viewed by others
Introduction
Photonic band gap (PBG) structures in photonic crystals (PCs) together with the characteristic dispersion properties have stimulated considerable interest in the study of fundamental photonic science and also in the development of new photonic technology1,2,3. The most significant new features induced by the PBG are the inhibition of atom spontaneous emission and the localization of light4,5,6. As a result, ultrahigh quality-factor cavity has been realized on-chip with PBG structures7. This provides the opportunity to control and manipulate light for photonic information technology8. In the past two decades, quantum optics with a few-level atom placed inside PCs have been extensively explored9 and cavity QED with the features of atomic population trapping and atom-photon bound states in the vicinity of the photonic band edge (PBE) has been examined10,11,12,13,14. These features are obtained mainly at zero temperature by solving the Schrödinger equation in which the PCs contain only one single photon emitted from an atom which is initially in the excited state. On the other hand, when the number of photons increases, light propagating in PCs was understood using the classical Maxwell equations, which is also defined at zero temperature1,3. Indeed, photonic quantum dynamics, even for a pure micro/nano cavity embedded in PCs, has not yet been solved at finite temperature. Practically, understanding photonic quantum dynamics at finite temperature is important for the development of all-optical circuits incorporating cavities and PBG waveguides embedded in PCs in the microwave regime.
All-optical circuits for networks on-chips consist of micro/nano cavities and waveguides8. Micro/nano cavities in PCs are created by point defects and PBG waveguides can be made with coupled defect arrays. Frequencies of cavities and waveguides can be easily tuned by changing the size and/or the shape of defects. In this article, we investigate micro/nano cavity photonics in PCs at finite temperature. Due to PBG-induced localized long-lived non-Markovian photon dynamics, we find that cavity photons in PCs do not obey Bose-Einstein statistical distribution. Within the PBG region and also in the vicinity of the PBE, cavity photons combine the nontrivial non-Markovian dissipations with thermal fluctuations together to form photon states that can memorize the initial cavity state information. As a result, Bose-Einstein statistical distribution for photons is completely broken down in these regimes, even though the photonic thermal energy is larger or much larger than the cavity detuning energy. Also, a crossover phenomenon from equilibrium to nonequilibrium steady states is revealed in low-dimensional photonic crytslas.
Results
Dissipation and dissipationless photon dynamics in photonic crystals
The methodology for the study of cavity photon dynamics in photonic crystals at finite temperature was developed in our previous work15,16,17,18, which is summarized in the Methods Section in the end of the article. To investigate the evolution of cavity photon states in PCs, we shall solve first the dynamics of photon dissipation and fluctuations, determined by Eqs. (19)–(20) given in Methods. Different spectral density J(ω) of the PCs, which requires the knowledge of DOS of the PCs, will provide different photon dissipation and fluctuations. In standard quantum optics, the spectral density J(ω) is usually treated as a constant in the weak-coupling limit, where the Weisskopf-Winger approximation or the Markovian master equation is valid so that the photon damping rate is time-independent (i.e. memoryless). As a result, all cavity photon modes have a finite lifetime and the cavity photon states will ultimately evolve into thermal equilibrium with the environment after decoherence takes place19 and photons inside the cavity must obey Bose-Einstein distribution20,21,22. However, this well-known result is no longer satisfied for micro/nano cavities in PCs, due to the presence of the PBG, as we will show below.
In principle, the DOS for different PCs, denoted by , should be calculated by solving photon eigenfrequencies and eigenfunctions from the Maxwell's equations for different photonic crystal structures2,23,24. In the literature, different photonic DOS have been introduced for different PCs9,10,11,12,14,25,26,27,28,29. For examples, for 1D PCs, the corresponding DOS is given by , where Θ(ω − ωe) is the Heaviside step function and ωe is the frequency at the PBE. This DOS is also often used for isotropic 3D PCs, which could predict qualitatively correct behaviors of the non-Weisskopf-Winger decay and the photon-atom bound state in PCs9,10,11,12,14, but it may overestimate the spontaneous emission rate of atoms due to the absence of the singularity of DOS in reality. Thus, for 3D PCs, the DOS near the PBE may be modeled by an anisotropic DOS9,10,11: , which could be further modified when the vectorial property of EM field is taken into account25,26. For 2D PCs with van Hove saddle point singularity27,28,29, the photonic DOS exhibits a logarithmic divergence near the PBE and can be approximated by , where ω0 is the center of the logarithmic peak. The spectral density J(ω) is microscopically defined as the multiplication of the DOS of the PCs with the photon tunneling amplitude V(ω) between the cavity and the PCs,
Assume that the cavity mode equally couples to all possible modes near the PBE of PCs, i.e. treating V(ω) as a constant, then the corresponding spectral densities J(ω) for different dimensional PCs are fully determined by the corresponding DOS, which are summarized in Table 1 and are also plotted in Fig. 1(a). We take the unit hereafter so that both the frequency ω and the spectral density J(ω) have the same unit as the energy.
Consider the cavity frequency ωc lies not too far away from the PBE. The dissipative photon dynamics is described by the cavity field propagating Green function u(t, t0) in photonic crystals through the relation 〈a(t)〉 = u(t, t0) 〈a(t0)〉 and is determined by the dissipative integro-differential equation (19). The general solution of Eq. (19) is given in Ref. 15,
where and Σ(ω) is the reservoir-induced cavity photon self-energy correction,
The explicit solution of cavity photon self-energy corrections for different spectral densities is also presented in Table 1. The frequency ωb in Eq. (2) is the localized photon mode frequency located inside the PBG (0 < ωb < ωe) but not too far away from the PBE and it is determined by the pole condition: ωb − ωc − Δ(ωb) = 0, where is a principal-value integral.
Equation (2) provides a general solution of the non-Markovian dissipative photon dynamics for micro/nano cavities in various photonic crystal structures. It shows that the cavity photon dynamics in PCs always contains two parts15: a localized photon mode (the first term) plus a non-exponential photon damping (the second term). The localized photon mode is a long-lived non-Markovian effect (dissipationless), induced by the PBG structure in PCs. The corresponding frequency ωb lies always within the PBG. The non-exponential photon damping comes from the non-analyticity of the photon self-energy correction, which is determined by the DOS profile of the PCs. This non-exponential damping is a short-time non-Markovian memory effect and it will become an exponential (Markov) decay in the photonic band (PB) region, as we will show later. The contributions of both the localized photon mode and the non-exponential photon damping strongly rely on the detuning δ = ωc − ωe.
In Fig. 1b, we plot the localized photon mode frequency as a function of the detuning δ for 1D, 2D and 3D PCs, respectively. The values of the localized mode frequencies for three different dimensional PCs are indeed very close to each other. However, a carefully check (see the inset in Fig. 1b) shows that for 3D PCs with an anisotropic DOS, the localized photon mode exists only when , namely, the cavity frequency ωc must be tuned into the PBG or near the PBE. Here χ is a coupling strength between the cavity and the PCs (see Table 1). For 1D and 2D PCs, the localized photon mode exists for any location of the cavity frequency ωc. This is because for 1D and 2D PCs with the DOS given in Table 1, the self-energy corrections are negatively divergent when ω → ωe − 0+ so that the pole condition is always satisfied for the localized photon mode. However, this does not ensure that the localized photon mode in 1D and 2D PCs must have a significant contribution to the photon dynamics for the whole range of the cavity frequency. The importance of the localized photon mode is determined by the localized photon mode amplitude, given in Eq. (2), which is plotted in Fig. 1c as a function of the detuning δ for all the three different dimensional PCs. It shows that for 3D PCs with an anisotropic DOS, the localized photon mode amplitude vanishes for , indicating a critical transition for the occurrence of the localized mode. For 1D and 2D PCs, the localized photon mode amplitude decays to almost zero when ωc is tuned into the PB region, as a crossover phenomenon. However, the overall effect of these localized photon modes on cavity photon dynamics has the very similar behavior for different PCs given in Table 1, as shown in Fig. 1.
The results presented in Fig. 1b–c provide indeed the full steady-state information of the cavity photon dynamics. This is because the non-exponential photon damping [i.e. the second term in Eq. (2)] will decay to zero so that only the localized photon mode contributes to the steady-state cavity photon field. In other words, the steady-state cavity field amplitude quantifies the contribution of the localized photon mode as a dissipationless effect, given by
On the other hand, the non-exponential damping term in Eq. (2) characterizes the time-dependent cavity photon damping rate κ(t) in the intermediate dissipation time region through the relation15,16
The detailed photon dynamics in terms of the cavity field amplitude |u(t, t0)| is provided in Fig. 2(a) for 1D, 2D and 3D PCs, respectively, with several different detuning δ taking from the PBG region to the PB region. The corresponding damping rate κ(t) is given in Fig. 2b, except for the case with δ = 0.1ωe which will be discussed later. The results show that the cavity dynamics changes dramatically when ωc crosses over from the PBG region to the PB region. Since the range of u(t, t0) is given by 1 ≥ |u(t, t0)| ≥ 0, we define the crossover region with the condition which corresponds to , as shown in Fig. 1. The results in Fig. 2b show that the damping rate κ(t) will rapidly approach to zero within the PBG (δ < −0.025ωe) and in the vicinity of the PBE (). In other words, the photon damping rate will approach to zero after some time, due to the existence of the localized photon dynamics. In the PB region (δ > 0.025ωe) where either the localized mode vanishes (for 3D PCs) or becomes negligible (for 1D and 2D PCs), then the cavity photons undergo a full dissipation process and can be approximately characterized by an exponential damping for , as shown by the black curves in Fig. 2a for δ = 0.1ωe. In this case, the damping rate κ(t) will oscillate rapidly in time when u(t, t0) decays to zero, see Eq. (5). This rapidly oscillating damping rate, originated from u(t, t0) approaching zero, has no physical consequence because the photon dissipation is almost completed after this point of time. Thus, alternatively we determine the damping rate κ by fitting the dashed-dot black curves in Fig. 2a with an exponential function ~exp(−κt). The resulting cavity photon damping rate for δ = 0.1ωe is given by , 1.24 × 10−2ωe and 1.72 × 10−2ωe for 1D, 2D and 3D PCs, respectively. Experimentally, the typical photonic band edge frequency ωe ranges from a few GHz to a few tens GHz for most of 1D, 2D and 3D PCs fabricated in the microwave regime30,31, this corresponds to the damping rate κ ranging from a few tens MHz to hundreds MHz when the cavity frequency ωc lies inside the PB region.
Combining the localized photon mode together with the photon dissipation dynamics, as shown in Fig. 1 and Fig. 2, we can see that when the cavity mode is tuned away from the PBE with δ > 0.025ωe, the localized mode has a negligible effect on cavity photon dynamics for 1D and 2D PCs and vanishes completely for 3D PCs. The cavity photon field decays almost exponentially for , which is a Markov decay, as shown in Fig. 2a by the dashed-dot black curves with δ = 0.1ωe. In the crossover region , i.e. the cavity frequency is tuned into the vicinity of the PBE, the contribution of the localized mode to the photon dynamics increases rapidly, while the contribution of the nonexponential photon damping goes in an opposite way, namely the damping rate κ(t) is decreased to zero after some time (see Fig. 2b). When the cavity frequency lies inside the PBG (δ < −0.025ωe), the cavity field has almost no damping (the damping rate κ(t) quickly decays to zero) and the photon dynamics is dominated by the localized photon mode. Thus light can be confined in the defect of the photonic crystal, providing a high-Q micro/nano cavity. In fact, this nontrivial dissipative cavity photon dynamics reproduces the same result with regard to atomic population trapping and atom-photon bound states in the vicinity of the PBG where an atom is placed in photonic crystals9,10,11,12.
The importance of the above results is that the localized photon mode provides dissipationless cavity photon dynamics when the cavity frequency ωc lies inside the PBG or in the vicinity of the PBE. This dissipationless localized-mode contribution to the steady-state cavity photon field is universal for different dimensional PCs with different DOS given in Table 1, because it only relies on the presence of PBG. The localized photon mode decreases rather quickly but smoothly in the crossover region () and becomes negligible for δ > 0.025ωe for 1D and 2D PCs, where the cavity photon dynamics becomes dissipative. For the 3D PCs, the cavity photon dynamics shows the same dissipationless behavior before reaching the critical point , where the localized photon mode dominates the photon dynamics inside the PBG (δ < −0.025ωe) and decays rapidly near the PBE (). When δ > δc, the localized photon mode vanishes for 3D PCs and the cavity photon dynamics becomes fully dissipative. With these results, a general picture for dissipation and dissipationless photon dynamics in photonic crystals is provided.
Thermal photon fluctuations
The above exact solution of dissipation and dissipationless cavity photon dynamics in photonic crystals can be used to describe thermal photon fluctuations through the fluctuated photon correlation function v(t, t), which is fully determined by the generalized non-equilibrium fluctuation-dissipation theorem15,
In Eq. (6), the two-time correlation function depicts photon fluctuations induced by the thermal photonic crystal, where is the initial photon distribution in PCs at temperature T. When the cavity approaches the steady state, photon fluctuations are simply determined by the modified steady-state fluctuation-dissipation theorem, with
where and . The first term is the localized photon mode contribution that modifies the conventional equilibrium fluctuation-dissipation theorem. This additional contribution is negligible when the cavity frequency is tuned into the PB region, i.e. δ > 0.025ωe, where ωb → ωe so that (see Fig. 1b–c). Consequently, the solution (7) is reduced to the standard equilibrium fluctuation-dissipation theorem: . In the high temperature limit, it reproduces Einstein's fluctuation-dissipation relation.
Photon fluctuations are presented in Fig. 3a(i)–(iii) respectively for 1D, 2D and 3D PCs at a given temperature. It shows that photon fluctuations evolve in a similar way for different PCs. When ωc is tuned into the PB region (δ > 0.025ωe), photons continuously flow into the cavity from the photonic crystal until the cavity reaches its steady state. In this case, the cavity thermally equilibrates with the photonic crystal. When ωc lies inside the PBE region (), photons also flow into the cavity at the beginning, but some of photons are then transmitted back into the photonic crystal, due to the environment-induced memory effects and the effect induced by the localized photon mode. When ωc is tuned into the PBG region (δ < −0.025ωe), photons hardly flow from the photonic crystal into the cavity, because the localized photon mode dominates the cavity photon dynamics in the PBG region where the reservoir-induced dissipation is negligible, even though the thermal energy, , is much larger than the detuning energy. Thus, thermal fluctuations are suppressed significantly by the photon localization in the PBG region and also near the PBE.
The physical picture of the above photon fluctuations can be seen clearly by connecting the photon correlations with the measurable average cavity photon number (i.e. field intensity)16,18,32: n(t) = 〈a†(t)a(t)〉 = |u(t, t0)|2 n(t0) + v(t, t), where n(t0) is the initial cavity intensity. It shows that the fluctuated photon correlation function v(t, t) is just the thermal-fluctuation-induced average photon number in the cavity. In Fig. 3b, we present the steady-state photon fluctuations at a given temperature, see the solid-blue curve. The dashed-pink curve is the thermal photon distribution in PCs at the same temperature: . It shows that when δ > 0.025ωe, the thermal-fluctuation-induced steady-state average photon number in the cavity is identical to the thermal photon distribution: . In this regime, the initial cavity photons are totally lost into PCs as u(t, t0) → 0 in the steady-state limit, as a consequence of the Weisskopf-Winger decay. The steady-state cavity photons all come from thermal fluctuations of the reservoir. As a result, the cavity photons approach equilibrium with the reservoir and obey Bose-Einstein distribution.
However, near the PBE (), the thermal-fluctuation-induced photon number deviates significantly from the thermal photon distribution, as shown in Fig. 3b. In other words, Bose-Einstein distribution no longer works for cavity photons in the PBE region, due to the localized photon effect. When the cavity frequency ωc lies inside the PBG (δ < −0.025ωe), the thermal fluctuations approach to zero and thus the cavity is fully determined by the localized photons, even though the cavity mode equally couples to all possible modes in PCs and the thermal energy is much larger than the detuning, . Consequently, Bose-Einstein distribution is broken down. This provides a nontrivial phenomenon of the thermal fluctuation dynamics from equilibrium to nonequilibrium steady states. Here nonequilibrium steady states are defined by these steady states which can memorize the initial state information of the system so that the equilibrium hypothesis cannot be satisfied. Such steady states will be given explicitly in the next subsections.
It is recently shown that non-Markovian dynamics in open systems can induce a critical transition from equilibrium to nonequilibrium steady states, when a localized mode occurs15,22,33. Such a critical transition appears for the 3D PCs with an anisotropic DOS, where there is a critical detuning for the occurrence of localized photon mode, . However, for the 1D and 2D PCs, our result shows that there is no such a critical transition due to the fact that the localized photon mode always exists. The photon fluctuations then undergoes a crossover from equilibrium to nonequilibrium steady states through the change of the detuning δ. No matter it is a crossover or a critical transition, the results in Fig. 3b shows that thermal photon fluctuations cannot make the cavities embedded in PCs to approach equilibrium with the PCs when the cavity frequency lies inside the PBG or near the PBE. This is a general feature for cavity photons in PCs.
Breakdown of Bose-Einstein distribution through the time-evolution of photonic Fock states
Based on the exact cavity photon dissipation and fluctuation dynamics in PCs, we shall explicitly examine the breakdown of Bose-Einstein distribution through the evolution of cavity photonic states and also the crossover from equilibrium to nonequilibrium steady states. This must be done by solving the exact master equation (15) given in the Methods. To be more specific, we consider first the cavity to be initially in a Fock state with an arbitrary photon number n0, i.e. ρ(t0) = |n0〉 〈n0|, which may be prepared experimentally through the real-time quantum feedback control34. By solving the master equation (15), the cavity state at arbitrary time t is given by
where . This result shows that an initial photon Fock state will evolve into a mixed state of different Fock states and the probability in each Fock state |n〉 is .
As we have shown, the breakdown of Bose-Einstein distribution relies on the dissipationless photon dynamics. And different dimensional PCs show almost the same dissipationless photon dynamics, determined by the localized photon mode due to the existence of the PBG. To be specific, we will present in the following the numerical solution of the state evolution for 1D PCs. The 2D and 3D PCs must provide a similar solution, based on the universal dissipation and fluctuation photon dynamics provided by Eqs. (2) and (6) and also the explicit numerical results presented in Fig. 2 and Fig. 3. The time-evolution of the cavity photon distribution for the initial state |n0 = 5〉 is given in Fig. 4. The steady-state limit, , is shown in Fig. 5, where several different initial states have been considered to demonstrate the initial-state dependence of the steady photon states.
The results show that if we tune ωc into the PB region (e.g. δ = 0.1ωe), then u(t, t0) will gradually decay to zero. This indicates that photons in the cavity will gradually be damped into the photonic crystal and photons in the photonic crystal are transferred into the cavity through thermal fluctuations. In this case, the contribution of the localized photon mode is negligible. The cavity state ultimately reaches thermal equilibrium with the photonic crystal and Bose-Einstein statistical distribution is produced,
where the initial state information is completely washed out and the cavity photon steady state is independent of the initial states, as shown in Fig. 4a and Fig. 5a.
Near the PBE (e.g. δ = 0), the cavity still loses photons into the photonic crystal and gains photons from the photonic crystal through thermal fluctuations. At a low temperature (), the photon loss and photon gain make the cavity become a mixed state of several Fock states |n〉 only for n < n0 and mainly distributed around n = n0/2, see Fig. 4b–(i) and also the steady-state limit by Fig. 5b–(i) where the initial state dependence is manifested. The photon distribution deviates significantly from the standard Bose-Einstein distribution. At a relatively high temperature (), the cavity state also becomes a mixed state of several Fock states |n〉, distributed mainly among n ≤ n0/2, but the distribution is broader, see Fig. 4b–(ii) and also the steady-state limit in Fig. 5b–(ii). In this case, the initial-state dependence of the steady state is still shown up and the derivation of the photon distribution from the standard Bose-Einstein distribution is obvious. As the temperature becomes rather high (), the time-evolution of the cavity state behaves quite differently because thermal fluctuations play a more important role now, thus the structure of the initial state is quickly destroyed, as shown in Fig. 4b–(iii). However, the cavity does not really approach thermal equilibrium with the photonic crystal, because is small but not negligible so that the distribution given by Eq. (9) is still different from the thermal state distribution. The cavity steady states for different initial states do not differ from each other as much as these in the low temperature cases, but the initial-state dependence is still observable, as shown in Fig. 5b–(iii).
When the cavity frequency ωc lies inside the PBG, e.g. δ = −0.1ωe given in Fig. 4c and Fig. 5c, the steady-state cavity photon distribution strongly depends on the initial state. In particular, if the photonic crystal temperature is not too high ( and ), the cavity remains in its initial Fock state as a photon localization state. It only has a small chance to decay to the Fock state |n〉 with n < n0 [see Fig. 4c(i) with ] and an even smaller probability to be in the Fock state |n〉 with n > n0 [see Fig. 4c(ii) when kBT increases to ]. The steady-state cavity photons are distributed mainly in the regime n < n0 with the maximum peak at n = n0, see Fig. 5c(i)–(ii). However, at a higher temperature , the thermal fluctuation becomes strong such that some of the initial state information will be lost and the cavity evolves into a mixed state covering several Fock states around the initial one |n0〉 [see Fig. 4c(iii)]. The steady-state cavity photons are distributed more broadly, but are still centered around the initial photon number n0, as shown in Fig. 5c(iii). This is because the localized photon mode dominates the photon dynamics over thermal photon fluctuations. The overall initial-state dependence of the steady states shown in Fig. 5b–c indicates that the equilibrium hypothesis in statistical mechanics is no longer obeyed and Bose-Einstein statistical distribution is completely broken down, even though the initial thermal energy of the photonic crystals is larger or much larger than the detuning energy.
Breakdown of Bose-Einstein distribution through the time-evolution of coherent states
Because Fock sates are highly non-classical and may not be easy to prepare in experiments, here we examine the case the cavity is initially in a coherent state |α0〉. By solving the master equation (15), the cavity state at an arbitrary later time t is given by
where is a displacement operator with α(t) = u(t, t0)α0, and
is a thermal-like state with average particle number v(t, t). Equation (11) shows that the initial cavity state will evolve into a displaced thermal-like state35,36 which is the mixture of displaced number states 37. In the photon number representation, Eq. (11) can be written as
where the probability of the cavity containing n photons (i.e. cavity photon distribution) is given by the diagonal terms of Eq. (13).
The cavity photon distribution in the steady-state limit is presented in Fig. 6. It shows again that when the cavity frequency ωc lies inside the PB region (δ = 0.1ωe), the cavity state (11) will evolve into the thermal state , which indicates that the cavity photon number obeys Bose-Einstein distribution, see Fig. 6a. If the cavity mode is tuned at the PBE (δ = 0), the cavity steady state will depend on its initial state, as shown in Fig. 6b. For low temperature (e.g. ), the thermal fluctuation is suppressed (v(t, t) → 0) and the cavity will evolve into the coherent state |α(t)〉, see Eq. (11). In this case, the cavity photon distribution in the steady state is Poissonian-like, with average photon number 〈n〉 |t→ ∞ = |α(t)|2 ∝ |α0|2, see Fig. 6b–(i). As the temperature goes higher (), the cavity photon distribution is deviated away from the Poisson distribution due to thermal fluctuations, see Fig. 6b–(ii). At high temperature (), since thermal fluctuation plays a more important role, the initial-state dependence is not as significant as that in the low temperature and the cavity photon distribution is not Poissonian-like, but still deviates from Bose-Einstein distribution, see Fig. 6b–(iii). When the cavity frequency lies inside the PBG (δ = −0.1ωe), the cavity steady state strongly depends on its initial state. The photon distribution in the cavity remains almost the same as in its initial state when the temperature is not very high ( and ), see the first two graphs in Fig. 6c. The initial-state dependence of the photon distribution is still strong when the temperature gets higher (), see Fig. 6c–(iii). The photon distribution in this case is slightly broadened from the initial states due to thermal fluctuations at high temperature, but still centered around n ≈ |α0|2, as a manifestation of photon localization. These results obtained for different initial coherent states at different temperature behave indeed in the same way as the results obtained for different initial Fock states under the same temperature.
Discussion
The breakdown of Bose-Einstein distribution in PCs, explored in details through the time-evolution of various cavity photon states for 1D PCs, is also valid for 2D and 3D PCs with different DOS listed in Table 1. This is because the time evolution of the cavity photon states, solved analytically by Eqs. (8) and (11), is fully determined by the photon dissipation and fluctuation dynamics of Eqs. (2) and (6). The steady-state behaviors of Eqs. (2) and (6) for different PCs are almost the same, as shown in Fig. 2 and Fig. 3. There is only a slight difference on the photon dissipation and fluctuation dynamics between 1D and 2D PCs with 3D PCs. For 1D and 2D PCs, varying the detuning leads to a crossover from equilibrium to nonequilibrium steady states for the photon dynamics, while for 3D PCs with an anisotropic DOS, it gives a critical transition, rather than a crossover, due to the existence of a critical detuning for the existence of the localized mode, see explicitly the results presented in Fig. 1 to Fig. 3. Such a critical transition is also found recently in other open systems22,33 but not for the crossover phenomenon. Thus, photon dissipation and fluctuation dynamics investigated in this work reveal a new nontrivial property, i.e. a crossover phenomenon for photon dynamics in low-dimensional photonic crystals.
In conclusion, we show that when the cavity frequency lies inside the PBG or near the PBE in photonic crystals, Bose-Einstein statistical distribution is broken down for cavity photons. This conclusion is generally valid for various photonic band gap structures in PCs. For the 1D and 2D PCs, the breakdown of Bose-Einstein distribution leads to a crossover from equilibrium to nonequilibrium cavity steady states, while for 3D PCs with an anisotropic DOS, the breakdown of Bose-Einstein distribution corresponds to a critical transition rather than a crossover. No matter whether it is a crossover or a critical transition, the breakdown of Bose-Einstein distribution is a consequence of localization photons due to the presence of PBG structures in PCs. Therefore the conclusion is also valid for other nanomaterials with band gap structures. It could provide a hitherto unexplored challenge on photon statistics. Furthermore, this non-trivial photon dynamics can be examined via quantum non-demolition measurement38,39, by sending sequences of circular Rydberg atoms through the photonic crystal micro/nano cavity, which carry information without destroying the cavity photon state. In particular, the cavity photon state can be measured using an experimental setup similar to that given in Ref. 38. Such experiments could be done with microcavities in low-dimensional photonic crystals in the microwave regime.
Methods
To investigate photon dynamics of micro/nano cavities which are coupled each other through waveguides embedded in PCs at finite temperature, we treat both the PCs and waveguides as reservoirs of the cavities. Thus the entire system of the micro/nano cavities (defects) embedded in photonic crystals can be described by the Fano-Anderson model (a model of impurity electrons coupled with continuous states introduced by Anderson40 in solid-state physics and discrete states embedded in a continuum proposed by Fano41 in atomic spectra). The corresponding Fano-Anderson Hamiltonian is given by Refs.16, 40,41,42,43:
where is the annihilation (creation) operator of the micro/nano cavity modes (defects) and the annihilation (creation) operator of the photonic modes of photonic crystals (continuum). The coefficients Vik are tunneling amplitudes of photons between the micro/nano cavities and the photonic crystals. The Hamiltonian (14) has also the same form as the one for cavity photon loss in open space20, by replacing the photonic crystals with an open space.
Consider initially the photonic crystals are in equilibrium state. Integrating out completely the reservoir degrees of freedom of photonic crystals via the influence functional44,45 in the coherent state representation46, arbitrary photon states for micro/nano cavities in PCs are then governed by the following exact master equation15,16,17,18,32,43,47,48
Here ρ (t) is the reduced density matrix for cavity states, is the renormalized Hamiltonian of cavities with the environment-modified cavity frequencies and the environment-induced couplings between different cavities , after the environmental degrees of freedom are completely integrated out. The coefficients κij (t) and characterize photon dissipations and fluctuations in PCs at finite temperature. The renormalized frequency, and the time-dependent dissipation and fluctuation coefficients, κij (t) and , can be exactly and non-perturbatively determined by the dissipation-fluctuation Dyson equation, as shown explicitly in our early works16,17 and will also be given later.
It might be worth pointing out that the first master equation derived from the original Feynman-Vernon influence functional44 was obtained by Caldeira and Leggett for the Brownian motion with a high temperature environment45 which is also called the Caldeira-Leggett master equation in the literature. For the Brownian motion, the system-reservoir coupling Hamiltonian is given by , where x and qk are the positions of the principal harmonic oscillator (as the Brownian particle) and all other harmonic oscillators in the reservoir, respectively. In terms of the second quantization, this system-reservoir coupling Hamiltonian can be rewritten as which shows that the amplitudes for particle tunneling processes between the system and the reservoir (the first two terms) is equal to the amplitudes of particle pair production and annihilation processes (the last two terms), both are given by Vk. However, in quantum optics, it is well-known that the photon tunneling amplitudes cannot be the same as the amplitudes of two-photon pair production and annihilation processes. As a result, the Caldeira-Leggett model is not applicable to photonics.
To make the cavity photon dynamics in thermal photonic crystals more specific, we consider a single-mode micro/nano cavity in photonic crystals. The photon dissipation and fluctuations, characterized by the dissipation and fluctuation coefficients κ (t) and (all the sub-indices (i, j) in Eq. (15) can be dropped now), are determined non-perturbatively and exactly by the nonequilibrium Green functions49,50 through the relations15,16,17,18:
Here u (t, t0) is the cavity photon field propagating Green's function describing the photon field relaxation and v (t, t) characterizes the reservoir-induced photon thermal fluctuations. These two Green functions, u (t, t0) and v (t, t), are determined by the following integrodifferential dissipation equation and the nonequilibrium fluctuation-dissipation theorem, respectively15,16,
where ωc is the original cavity frequency. The integral kernels in Eqs. (19)–(20) characterize all the back-actions between the cavity and photonic crystals and can be determined uniquely by the spectral density J (ε) of photonic crystals through the relations: and , where is the initial photon distribution in PCs at temperature T. The spectral density J(ω) is microscopically defined as a multiplication of the density of states (DOS) of PCs with the photon scattering amplitudes Vk between the cavity and PCs,
In the second equality we have taken the continuous photonic modes of PCs so that Vk → V(ω) and the index i of Vik in Eq. (14) has also been dropped for single-mode cavity. For arbitrary spectral density J(ω), the general solution of Eq. (19) can be obtained exactly15, which is given by Eq. (2). By solving Eqs. (16)–(18) and (15) through Eqs. (19)–(20), the complete solution of photon dissipation and fluctuations can be obtained, as we presented in details in the paper.
References
Joannopoulos, J. D., Johnson, S. G., Winn, J. N. & Meade, R. D. Photonic Crystals: Modeling the Flow of Light (Princeton, New York, 2008).
Busch, K., von Freymann, G., Linden, S., Mingaleev, S. F., Tkeshelashvili, L. & Wegener, M. Periodic nanostructures for photonics. Phys. Rep. 444, 101 (2007).
Notomi, M. Manipulating light with strongly modulated photonic crystals. Rep. Prog. Phys. 73, 096501 (2010).
Yablonovitch, E. Inhibited Spontaneous Emission in Solid-State Physics and Electronics. Phys. Rev. Lett. 58, 2059 (1987).
John, S. Strong localization of photons in certain disordered dielectric superlattices. Phys. Rev. Lett. 58, 2486 (1987).
Yablonovitch, E., Gmitter, T. J. & Leung, K. M. Photonic band structure: The face-centered-cubic case employing nonspherical atoms. Phys. Rev. Lett. 67, 2295 (1991).
Armani, D. K., Kippenberg, T. J., Spillane, S. M. & Vahala, K. J. Ultra-high-Q toroid microcavity on a chip. Nature 421, 925 (2003).
Joannopoulos, J. D., Villeneuve, P. R. & Fan, S. Photnic crystals: putting a new twist on light. Nature 386, 143 (1997).
Lambropoulos, P., Nikolopoulos, G. M., Nielsen, T. R. & Bay, S. Fundamental quantum optics in structured reservoirs. Rep. Prog. Phys. 63, 455 (2000).
John, S. & Wang, J. Quantum Electrodynamics near a Photonic Band Gap: Photon Bound States and Dressed Atoms. Phys. Rev. Lett. 64, 2418 (1990).
John, S. & Wang, J. Quantum optics of localized light in a photonic band gap. Phys. Rev. B 43, 12772 (1991).
Kofman, A. G., Kurizki, G. & Sherman, B. Spontaneous and Induced Atomic Decay in Photonic Band Structures. J. Mod. Opt. 41, 353 (1994).
Zhu, S. Y., Yang, Y., Chen, H., Zheng, H. & Zubairy, M. S. Spontaneous Radiation and Lamb Shift in Three-Dimensional Photonic Crystals. Phys. Rev. Lett. 84, 2136 (2000).
Wang, X. H., Gu, B. Y., Wang, R. & Xu, H. Q. Decay kinetic properties of atoms in photonic crystals with absolute gaps. Phys. Rev. Lett. 91, 113904 (2003).
Zhang, W. M., Lo, P. Y., Xiong, H. N., Tu, M. W. Y. & Nori, F. General Non-Markovian Dynamics of Open Quantum Systems. Phys. Rev. Lett. 109, 170402 (2012).
Lei, C. U. & Zhang, W. M. A quantum photonic dissipative transport theory. Ann. Phys. 327, 1408 (2012).
Tu, M. W. Y. & Zhang, W. M. Non-Markovian decoherence theory for a double-dot charge qubit. Phys. Rev. B 78, 235311 (2008).
Jin, J., Tu, M. W. Y., Zhang, W. M. & Yan, Y. J. Non-equilibrium quantum theory for nanodevices based on the Feynman-Vernon influence functional. New. J. Phys. 12, 083013 (2010).
Brune, M. et al. Process Tomography of Field Damping and Measurement of Fock State Lifetimes by Quantum Nondemolition Photon Counting in a Cavity. Phys. Rev. Lett. 101, 240402 (2008).
Carmichael, H. J. An Open Systems Approach to Quantum Optics, Lecture Notes in Physics Vol. m18, (Springer-Verlag, Berlin, 1993).
Scully, M. O. & Zubairy, M. S. Quantum Optics (Cambridge University press, 1997).
Xiong, H. N., Lo, P. Y., Zhang, W. M., Nori, F. & Feng, D. H. Non-Markovian Dynamics Impact on the Foundations of Statistical Mechanics. arXiv: 1311.1282 (2013) (unpublished).
Noda, S. & Baba, T. Roadmap on photonic crystals (Dordrecht: Kluwer Academic Publishers Groups, 2003).
Sakoda, K. Optical properties of photonic crystals (Berlin, Springer, 2005).
Li, Z. Y., Lin, L. L. & Zhang, Z. Q. Spontaneous Emission from Photonic Crystals: Full Vectorial Calculations. Phys. Rev. Lett. 84, 4341 (2000).
Li, Z. Y. & Xia, Y. Full vectorial model for quantum optics in three-dimensional photonic crystals. Phys. Rev. A 63, 043817 (2001).
McPhedran, R. C., Botten, L. C., McOrist, J., Asatryan, A. A., de Sterke, C. M. & Nicorovici, N. A. Density of states functions for photonic crystals. Phys. Rev. E 69, 016609 (2004).
Birner, A., Wehrspohn, R. B., Gösele, U. M. & Busch, K. Silicon-Based Photonic Crystals. Adv. Mater. 13, 377 (2001).
Ibanescu, M., Reed, E. J. & Joannopoulos, J. D. Enhanced Photonic Band-Gap Confinement via Van Hove Saddle Point Singularities. Phys. Rev. Lett. 96, 033904 (2006).
Smith, D. R., McCall, S. L., Platzman, P. M., Dalichaoich, R., Kroll, N. & Schultz, S. Photonic band structure and defects in one and two dimensions. J. Opt. Soc. Am. B 10, 314 (1993).
Özbay, E. et al. Measurement of a three-dimensional photonic band gap in a crystal structure made of dielectric rods. Phys. Rev. B 50, 1945 (1994).
Xiong, H. N., Zhang, W. M., Wang, X. & Wu, M. H. Exact non-Markovian cavity dynamics strongly coupled to a reservoir. Phys. Rev. A 82, 012105 (2010).
Cai, C. Y., Yang, L. P. & Sun, C. P. Threshold for nonthermal stabilization of open quantum systems. Phys. Rev. A 89, 012128 (2014).
Sayrin, C. et al. Real-time quantum feedback prepares and stabilizes photon number states. Nature 477, 73 (2011).
Lachs, G. Theoretical Aspects of Mixtures of Thermal and Coherent Radiation. Phys. Rev. 138, B1012 (1965).
Marian, P. & Marian, T. A. Squeezed states with thermal noise. I. Photon-number statistics. Phys. Rev. A 47, 4474 (1993).
de Oliveira, F. A. M., kim, M. S., Knight, P. L. & Bužek, V. Properties of displaced number states. Phys. Rev. A 41, 2645 (1990).
Gleyzes, S. et al. Quantum jumps of light recording the birth and death of a photon in a cavity. Nature 446, 297 (2007).
Guerlin, C. et al. Progressive field-state collapse and quantum non-demolition photon counting. Nature 448, 889 (2007).
Anderson, P. W. Localized magnetic states in metals. Phys. Rev. 124, 41 (1961).
Fano, U. Effects of Configuration interaction on intensities and phase shift. Phys. Rev. 124, 1866 (1961).
Mahan, G. D. Many-Body Physics 3rd Ed. (Kluwer Academic/Plenum Publishers, New Yoek, 2000).
Wu, M. H., Lei, C. U., Zhang, W. M. & Xiong, H. N. Non-Markovian dynamics of a microcavity coupled to a waveguide in photonic crystals. Opt. Express 18, 18407 (2010).
Feynman, R. P. & Vernon, F. L. The theory of a general quantum system interacting with a linear dissipative system. Ann. Phys. 24, 118 (1963).
Caldeira, A. O. & Leggett, A. J. Path integral approach to quantum Brownian motion. Physica 121A, 587 (1983).
Zhang, W. M., Feng, D. H. & Gilmore, R. Coherent states: theory and some applications. Rev. Mod. Phys. 62, 867 (1990).
Tan, H. T. & Zhang, W. M. Non-Markovian dynamics of an open quantum system with initial system-reservoir correlations: A nanocavity coupled to a coupled-resonator optical waveguide. Phys. Rev. A 83, 032102 (2011).
Tang, L. P., Cai, C. Y., Xu, D. Z., Zhang, W. M. & Sun, C. P. Master equation and dispersive probing of a non-Markovian process. Phys. Rev. A 87, 012110 (2013).
Schwinger, J. Brownian Motion of a Quantum Oscillator. J. Math. Phys. 2, 407 (1961).
Keldysh, L. V. Diagram Technique for Nonequilibrium Processes. Sov. Phys. JETP 20, 1018 (1965).
Acknowledgements
This research is supported by the National Science Conculs of ROC under Contract No. NSC-102-2112-M-006-016-MY3 and the National Center for Theoretical Science of Taiwan. It is also supported in part by the Headquarters of University Advancement at the National Cheng Kung University, which is sponsored by the Ministry of Education of ROC.
Author information
Authors and Affiliations
Contributions
W.M.Z. proposed the ideas, interpreted the physics and wrote the main manuscript, P.Y.L. performed the theoretical calculations and prepared all the figures and H.N.X. checked the results and participated the discussions. All authors reviewed the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder in order to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Lo, PY., Xiong, HN. & Zhang, WM. Breakdown of Bose-Einstein Distribution in Photonic Crystals. Sci Rep 5, 9423 (2015). https://doi.org/10.1038/srep09423
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep09423
This article is cited by
-
Modeling of global temperature control
Environment, Development and Sustainability (2021)
-
Quantum thermodynamics of single particle systems
Scientific Reports (2020)
-
Spin-polarized magneto-electronic properties in buckled monolayer GaAs
Scientific Reports (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.