arXiv is now an independent nonprofit! Learn more
License: arXiv.org perpetual non-exclusive license
arXiv:2403.00847v1 [quant-ph] 29 Feb 2024

Negative refraction without absorption via both coherent and incoherent fields in a four-level left-handed atomic systemSupported by National Natural Science Foundation of China (Grant No.60768001 and No.10464002).

Shun-Cai ZhaoID Corresponding author: zscnum1@126.com. Affiliation: School of Materials Science and Engineering, Nanchang University, Nanchang 330031, PR China Affiliation: Engineering Research Center for Nanotechnology, Nanchang University, Nanchang 330047, PR China Affiliation: Institute of Modern Physics,Nanchang University, Nanchang 330031, PR China    Zheng-Dong Liu Corresponding author: lzdgroup@ncu.edu.cn. Affiliation: School of Materials Science and Engineering, Nanchang University, Nanchang 330031, PR China Affiliation: Engineering Research Center for Nanotechnology, Nanchang University, Nanchang 330047, PR China Affiliation: Institute of Modern Physics,Nanchang University,Nanchang 330031, PR China    Qi-Xuan Wu Affiliation: College English department,Hainan University, Danzhou 571737, PR China
Abstract

This paper attempts a probe into negative refraction without absorption by means of an incoherent pump field and a strong coherent field coupling the dense four-level atomic system.With the application of the incoherent pump field to manipulate the populations in atomic levels and the variable strong coherent field to create quantum coherence, the constraint condition of two equal transition frequencies responding to the probe field in the atomic system isn’t required.And these lead to the propagation transparency and strong magnetic response of the probe field,left-handedness with vanishing absorption in the atomic system.However,an excessive coherent field intensity would increase the absorption.

PACs

42.50.Gy

Keywords

negative refraction; without absorption; left-handed; coherent field; incoherent field.

I Introduction

Negative refraction of electromagnetic radiation[1]has recently attracted considerable attention because of its surprising and counterintuitive electromagnetical and optical effects,such as the reversals of both Doppler shift and Cherenkov effect,negative refraction[2],amplification of evanescent waves[3]and subwavelength focusing [3-5], negative Goos-Ha¨\ddot{a}nchen shift [6] and quenching spontaneous emission[7,8]and so on[9].And material with negative refraction index intrigues the researchers because of its many significant potential applications.The“ perfect lens ”is one of them.Since a slab of such materials has an ability to focus all frequency components of a two-dimensional image, it may become possible to make a “ perfect lens ”in which imaging resolution does not limited by the diffraction limit[3].Up to now, there have been several approaches to the realization of negative refractive index materials, including artificial composite metamaterials [10-11], photonic crystal structures [12],transmission line simulation[13] and chiral media[14-15] as well as photonic resonant materials(coherent atomic vapour)[16-19].In such a type of negative refractive materials, negative refraction is usually accompanied by a strong absorption especially towards higher frequencies. Thus,the realization of negative refraction material without absorption is of great significance.And some effort[20-23]has been made to realize negative refraction without absorption.Ka¨\ddot{a}stel et al[20] realized negative refraction with minimal absorption in a dense atomic gas via electromagnetically induced chirality. The key ingredient of the scheme is the electromagnetic chirality that results from coherently coupling a magnetic dipole transition with an electric dipole transition via atomic coherence induced by the two-photon resonant Raman transitions. In Ref.[21],Ka¨\ddot{a}stel et al also discussed negative refraction with reduced absorption due to destructive quantum interference in coherently driven atomic media.The negative refraction with deeply depressing absorption and without simultaneously requiring both negative electric permittivity and magnetic permeability (left-handedness)[2]was obtained by F.L. Li[23].Ref.[23] shows this at the ideal situation that the two chirality coefficients have the same amplitude but the opposite phase.

In this paper we propose an indirect coupling way of the atoms responding to the probe field via an incoherent pump field and a coherent coupling field.With the application of the incoherent pump field to manipulate the population distribution of each level and the variable coherent field to create quantum coherence,the magnetic response is amplified and the probe field propagates transparently. The atomic system displays negative refraction without absorption and left-handedness.The constraint condition of two equal transition frequencies responding to the probe field does not require,which is different from Ref.[18,24].

II Theoretical model

The level configuration of atoms under consideration is shown in Figure 1. The properties of the four atomic states are as following: levels |1|1\rangle,|3|3\rangle, and |4|4\rangle have same parity, and the parity of level|2|2\rangle is opposite with theirs. The two lower levels |1|1\rangle and |2|2\rangle have opposite parity and so 2|\langle 2|d^\hat{\vec{d}}|1|1\rangle0\neq 0 where d^\hat{\vec{d}} is the electric dipole operator. The two upper levels,|3|3\rangle and |4|4\rangle have the same parity with 4|\langle 4|μ^\hat{\vec{\mu}}|3|3\rangle0\neq 0 where μ^\hat{\vec{\mu}} is the magnetic-dipole operator. As observed in Fig.1,three electromagnetic fields are introduced to couple the four states:The electric and magnetic components of the probe light (corresponding Rabi frequency Ωp\Omega_{p}=EPd21=\frac{\vec{E_{P}}\vec{d_{21}}}{\hbar} ,ΩB\Omega_{B}=BPμ43=\frac{\vec{B_{P}}\vec{\mu_{43}}}{\hbar}) interact with the transitions |2|2\rangle and|1|1\rangle as well as|4|4\rangle and|3|3\rangle, respectively.The incoherent pump field with pumping rate denoted by Γ\Gamma pumps atoms in level |1|1\rangle into upper level |3|3\rangle, and then the atoms decay into metastable level |2|2\rangle via rapid nonradiative transitions, whose decay rate is denoted as Γ32\Gamma_{32}.The strong coherent field Ωc\Omega_{c} takes two effects on the system:(i) populations in level |2|2\rangle are being pumped into level |4|4\rangle and causes population reversion gain for magnetic-dipole transition |3|3\rangle\leftrightarrow|4|4\rangle. (ii)quantum coherence is induced by it and causes energy levels |2|2\rangle and |4|4\rangle split into two dressed sublevels and results in enough response of the electric part of probe field at certain frequency extent.

Refer to caption
Figure 1: Four-level atomic medium interacting with a coherent Ωc\Omega_{c} , an incoherent pumps Γ\Gamma and a probe field Ωp\Omega_{p}.

Using the density-matrix approach, the time-evolution of the system is described as

dρdt=i[H,ρ]+Λρ,\frac{d\rho}{dt}=-\frac{i}{\hbar}[H,\rho]+\Lambda\rho, (1)

Where Λρ\Lambda\rho represents the irreversible decay part in the system.Under the dipole approximation and the rotating wave approximation the density matrix equations described the system are written as follows:

ρ11˙=Γ(ρ33ρ11)+Γ21ρ22+Γ31ρ33+Γ41ρ44+iΩp(ρ21ρ12),\dot{\rho_{11}}=\Gamma(\rho_{33}-\rho_{11})+\Gamma_{21}\rho_{22}+\Gamma_{31}\rho_{33}+\Gamma_{41}\rho_{44}+i\Omega_{p}(\rho_{21}-\rho_{12}), (2)
ρ21˙=(γ21+iΔp)ρ21iΩp(ρ22ρ11)+iΩcρ41,\dot{\rho_{21}}=-(\gamma_{21}+i\Delta_{p})\rho_{21}-i\Omega_{p}(\rho_{22}-\rho_{11})+i\Omega_{c}\rho_{41}, (3)
ρ22˙=Γ21ρ22+Γ32ρ32+Γ42ρ44+iΩp(ρ12ρ21)+iΩc(ρ42ρ24),\dot{\rho_{22}}=-\Gamma_{21}\rho_{22}+\Gamma_{32}\rho_{32}+\Gamma_{42}\rho_{44}+i\Omega_{p}(\rho_{12}-\rho_{21})+i\Omega_{c}(\rho_{42}-\rho_{24}), (4)
ρ31˙=[γ31+i(Δc+δ)]ρ31+iΩpρ23+iΩBρ41,\dot{\rho_{31}}=-[\gamma_{31}+i(\Delta_{c}+\delta)]\rho_{31}+i\Omega_{p}\rho_{23}+i\Omega_{B}\rho_{41}, (5)
ρ32˙=[γ32+i(ΔcΔp+δ)]ρ32iΩcρ34iΩpρ31+iΩBρ42,\dot{\rho_{32}}=-[\gamma_{32}+i(\Delta_{c}-\Delta_{p}+\delta)]\rho_{32}-i\Omega_{c}\rho_{34}-i\Omega_{p}\rho_{31}+i\Omega_{B}\rho_{42}, (6)
ρ33˙=Γ(ρ33ρ11)Γ31ρ33Γ32ρ33+Γ43ρ44+iΩB(ρ43ρ34),\dot{\rho_{33}}=-\Gamma(\rho_{33}-\rho_{11})-\Gamma_{31}\rho_{33}-\Gamma_{32}\rho_{33}+\Gamma_{43}\rho_{44}+i\Omega_{B}(\rho_{43}-\rho_{34}), (7)
ρ41˙=[γ41+i(Δp+Δc)]ρ41+iΩcρ21iΩpρ42+iΩBρ31,\dot{\rho_{41}}=-[\gamma_{41}+i(\Delta_{p}+\Delta_{c})]\rho_{41}+i\Omega_{c}\rho_{21}-i\Omega_{p}\rho_{42}+i\Omega_{B}\rho_{31}, (8)
ρ42˙=(γ42+iΔc)ρ42+iΩc(ρ22ρ44)iΩpρ41+iΩBρ32,\dot{\rho_{42}}=-(\gamma_{42}+i\Delta_{c})\rho_{42}+i\Omega_{c}(\rho_{22}-\rho_{44})-i\Omega_{p}\rho_{41}+i\Omega_{B}\rho_{32}, (9)
ρ43˙=[γ43+i(Δpδ)]ρ43+iΩcρ23+iΩB(ρ33ρ44),\dot{\rho_{43}}=-[\gamma_{43}+i(\Delta_{p}-\delta)]\rho_{43}+i\Omega_{c}\rho_{23}+i\Omega_{B}(\rho_{33}-\rho_{44}), (10)

where the above density matrix elements obey the conditions:ρ11+ρ22+ρ33+ρ44\rho_{11}+\rho_{22}+\rho_{33}+\rho_{44}=1 and ρij\rho_{ij}=ρji\rho_{ji}^{\ast}.And Γij\Gamma_{ij}(i,j=1,2,3,4)is the spontaneous emission decay rate from level |i|i\rangle to level|j|j\rangle,ignoring the collision broaden effect.γ21\gamma_{21}=Γ21\Gamma_{21}/2,γ31\gamma_{31}=(Γ31+Γ32\Gamma_{31}+\Gamma_{32})/2, γ41\gamma_{41}=(Γ43+Γ42\Gamma_{43}+\Gamma_{42})/2, γ42\gamma_{42}=(Γ43+Γ31+Γ21\Gamma_{43}+\Gamma_{31}+\Gamma_{21})/2, γ43\gamma_{43}=(Γ43+Γ42+Γ31+Γ32\Gamma_{43}+\Gamma_{42}+\Gamma_{31}+\Gamma_{32})/2,γ32\gamma_{32}= (Γ32+Γ31+Γ21\Gamma_{32}+\Gamma_{31}+\Gamma_{21})/2 are the decay rates to the corresponding transitions.The detuning of the fields defined as Δp\Delta_{p}=ω21\omega_{21}-ωp\omega_{p},Δc\Delta_{c}=ω42\omega_{42}-ωc\omega_{c}, ΔB\Delta_{B}=ω43\omega_{43}-ωP\omega_{P},respectively.ωij\omega_{ij}=ωi\omega_{i}-ωj\omega_{j} is the transition frequency of level|i|i\rangle and |j|j\rangle (i, j=1,2,3,4) and we have δ\delta =Δp\Delta_{p}-ΔB\Delta_{B}.

According to the classical electromagnetic theory,the electric polarizability is a rank 2 tensor defined by its Fourier transform Pe(ωP)\vec{P}_{e}(\omega_{P}) =ϵ0=\epsilon_{0}αe(ωP)\alpha_{e}(\omega_{P})E(ωP)\vec{E}(\omega_{P}),which is calculated as the mean value of the atomic electric-dipole moment operators by the definition Pe\vec{P}_{e} =Tr{\{ρ^d{\hat{\rho}\vec{d}}}\}=ρ12d21\rho_{12}d_{21}+c.c.where Tr stands for trace.In the following, we only consider the polarizability at the frequency ωP\omega_{P} of the incoming field Ep\vec{E}_{p}.Therefore we drop the explicit ωP\omega_{P} dependence αe(ωP)αe\alpha_{e}(\omega_{P})\equiv\alpha_{e}.Moreover,we choose Ep\vec{E}_{p} parallel to the atomic dipole d21\vec{d}_{21} so that αe\alpha_{e}is a scalar,and its expression is as follows:

αe=d21ρ12ϵ0Ep=d212ρ12ϵ0Ωp,\alpha_{e}=\frac{\vec{d}_{21}\rho_{12}}{\epsilon_{0}\vec{E}_{p}}=\frac{\mid{d_{21}}\mid^{2}\rho_{12}}{\epsilon_{0}\hbar\Omega_{p}}, (11)

In the same way, the classical magnetic polarizations of the medium Pm(ωP)\vec{P}_{m}(\omega_{P})=μ0αmE(ωP)\mu_{0}\alpha_{m}\vec{E}(\omega_{P}),which is related to the mean value of the atomic dipole moment operator through Pm\vec{P}_{m} =Tr{\{ρ^μ{\hat{\rho}\vec{\mu}}}\}=ρ34\rho_{34}μ43+c.c\mu_{43}+c.c.According to the classical Maxwell’s electromagnetic wavevector relation,we choose magnetic dipole is perpendicular to the induced electric dipole so that the magnetizability αm\alpha_{m} is scalar, and its expression is as follows:

αm=μ0μ43ρ34Bp=μ0μ432ρ34ΩB.\alpha_{m}=\frac{\mu_{0}\vec{\mu}_{43}\rho_{34}}{\vec{B}_{p}}=\frac{\mu_{0}\mid\mu_{43}\mid^{2}\rho_{34}}{\hbar\Omega_{B}}. (12)

According to the Clausius-Mossotti relations considering the local effect in dense medium[25], the relative permittivity and relative permeability are expressed as [18,26]

ϵr=1+23Nαe113Nαe,\displaystyle\epsilon_{r}=\frac{1+\frac{2}{3}N\alpha_{e}}{1-\frac{1}{3}N\alpha_{e}}, (13)
μr=1+23Nγm113Nγm.\displaystyle\mu_{r}=\frac{1+\frac{2}{3}N\gamma_{m}}{1-\frac{1}{3}N\gamma_{m}}. (14)

In the above, we obtained the expressions for the electric permittivity and magnetic permeability of the atomic media.In the section that follows, we will demonstrate that the simultaneously negative both permittivity and permeability,negative refraction without absorption can be observed in the four-level atomic system.

III Results and discussion

In the following,with the stationary solutions to the density-matrix equations(2)-(10), we explore the property of simultaneously negative both electric permittivity and magnetic permeability through the numerical calculations.And several typical parameters should be selected before the calculation.The parameters for the electric and magnetic polarizabilities of atoms can be chosen as: electric and magnetic transition dipole moments d21d_{21}=2.5×\times102910^{-29}C\cdotm and μ34\mu_{34}=7.0 ×\times102310^{-23}Cm2s1\cdot m^{2}s^{-1}[16], respectively.The density of atoms N was choosed to be 6.5×1025m36.5\times 10^{25}m^{-3}[27,28].And the other parameters are scaled by γ=106s1\gamma=10^{6}s^{-1}: Γ31\Gamma_{31}=Γ41\Gamma_{41}=Γ42\Gamma_{42}=Γ43\Gamma_{43}= 0.01γ\gamma,Γ21\Gamma_{21}=0.5γ\gamma, Γ32\Gamma_{32}=1γ\gamma. The Rabi frequency of the probe field is Ωp=0.5γ\Omega_{p}=0.5\gamma ,and δ\delta = - 20γ\gamma without the condition of the two transition frequencies ω43\omega_{43}=ω21\omega_{21}. The incoherent pumping rate is Γ=1.0γ\Gamma=1.0\gamma. The strong coherent optical field couples the atomic system with variational Rabi frequencies Ωc=15γ\Omega_{c}=15\gamma, 18γ\gamma,21γ\gamma, 24γ\gamma and fixed frequency detuning Δc\Delta_{c}=20γ\gamma.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 2: The relative electric permittivity ϵr\epsilon_{r} and magnetic permeability μr\mu_{r} of the system as a function of probe frequency detuning Δp/γ\Delta_{p}/\gamma for different values of ΩC\Omega_{C}.

Figure 2 shows the calculated electric permittivity ϵr\epsilon_{r} and magnetic permeability μr\mu_{r} as a function of probe field detuning.Under the cooperation of the strong coherent field and incoherent pumping,the dense four-level atomic system exhibits left-handedness with simultaneously negative permittivity and permeability at some frequency extents of the probe field.From the profiles of the real part of the relative electric permittivity ϵr\epsilon_{r},we observe their values being proximately equal -2 in the probe frequency range[-26γ\gamma,26γ\gamma],although the phase of the figures are proximately opposite at Δp\Delta_{p}=12.5γ\gamma when the Rabi frequency ΩC\Omega_{C} varies by 15γ\gamma, 18γ\gamma, 21γ\gamma and 24γ\gamma.The influence of the strong coherent field on the real part of the relative magnetic permeability μr\mu_{r} is different.When ΩC\Omega_{C}=15γ\gamma,the value is all negative in the probe frequency extent.When ΩC\Omega_{C}=18γ\gamma,21γ\gamma, 24γ\gamma,the positive values emerge in the ranges of [-3γ\gamma, -1.1γ\gamma], [-3.4γ\gamma, -0.4γ\gamma], [-3.85γ\gamma,0.36γ\gamma],respectively.And their amplitude values are increasing with the variation of ΩC\Omega_{C}.However,the ranges for simultaneously negative electric permittivity ϵr\epsilon_{r} and magnetic permeability μr\mu_{r} still exist.In Figure 2,the transparency is shown by the imaginary part of both the relative electric permittivity ϵr\epsilon_{r} and the relative magnetic permeability μr\mu_{r}in some frequency extents.This is a very significant result for us.As is well known, the photon absorption of atom can be greatly depressed via electromagnetically induced transparency (EIT)[29].And the zero absorption phenomena may occur in the EIT extents.With observation of the profile of the imaginary part of ϵr\epsilon_{r}, the excessive strong coherent field intensity causes the shrinking gain and increasing absorption at Δp\Delta_{p}=12.5γ\gamma.However,the imaginary part of μr\mu_{r} displays increasing gain,and its amplitude is gradually amplified and gradually approaching to the resonant point when the coherent field is varied the Rabi frequencies by 15γ\gamma, 18γ\gamma,21γ\gamma and 24γ\gamma. Comparing the images of the relative electric permittivity ϵr\epsilon_{r} with the relative magnetic permeability μr\mu_{r} in Figure 2,we notice that the magnetic response of the probe field is stronger than the electric response.The reason may come from that: the magnetic component of the probe field couples levels |4|4\rangle and |3|3\rangle,and the strong coherent field drives levels |2|2\rangle and |4|4\rangle.The effect of quantum coherence caused by the strong coherent field on the magnetic component is stronger than on the electric component which couples levels |2|2\rangle and |1|1\rangle.

Refer to caption
Refer to caption
Figure 3: The refractive index as a function of the probe field detuning Δp/γ\Delta_{p}/\gammafor different values of ΩC\Omega_{C}.
Refer to caption
Figure 4: The figure of merit (FOM: |real(n)/imag(n)||real(n)/imag(n)|) as a function of the detuning parameter detuning Δp/γ\Delta_{p}/\gammafor different values of ΩC\Omega_{C}.

In Figure 3,the refraction index according to definition of the left-handed material ( n(ω)=ϵr(ω)μr(ω)n(\omega)=-\sqrt{\epsilon_{r}(\omega)\mu_{r}(\omega)} ) [2] is plotted for different values of parameter ΩC\Omega_{C}.As observed in Figure 3,the real part of the refractive index shows negative values and their amplitudes enlarge gradually with the increasing of ΩC\Omega_{C}. The imaginary part of the refractive index displays absorption enhancing near the resonant point and absorption depression on the both sides,when the coherent field is varied the Rabi frequencies by 15γ\gamma, 18γ\gamma,21γ\gamma and 24γ\gamma.The figure of merit(FOM)|Re(n)/Im(n)||Re(n)/Im(n)| shows how much the absorption is suppressed .When the FOM is much larger than unity, it means that there is almost no absorption in this area. As shown in Figure 4, the FOM is far less than unity in the area near the resonant point.This illustrates that the increasing absorption occurs in this area ,which is consistent with the peaks of Im[n] near the resonance in Figure 3. And the FOM has the largest value at the identical re-scaled detuning parameter Δp/γ\Delta_{p}/\gamma when ΩC\Omega_{C}= 15γ\gamma. It means that the absorption is depressed deeply when the coherent field varies its Rabi frequency to ΩC\Omega_{C}= 15γ\gamma,and the excessive coherent field intensity wouldn’t help to realize zero absorption.The markedly feature of present scheme is the absorption depression almost to zero in the left-handed atomic system. One may get the reason for this from the results shown in Figure 2:Because of the EIT effect, the electric and magnetic non-absorption extents appear in the re-scaled detuning parameter Δp/γ\Delta_{p}/\gamma extent,which brings about the zero value of the imaginary part of refraction index. As a result, the negative refraction without absorption occurs in the left-handed atomic system .

In experimental investigation, the solid sample of Er3+E_{r}^{3+} doped into calcium fluorophosphate at room temperature might be a good candidate because it has abundant energy levels, various electric magnetic transitions and high density in the order 1025m310^{25}m^{-3}[27,28].And the four levels |1|1\rangle,|2|2\rangle, |3|3\rangle,and |4|4\rangle in Figure 1 correspond to the energy level configuration I15/24I_{15/2}^{4},I13/24I_{13/2}^{4}, I9/24I_{9/2}^{4} and I11/24I_{11/2}^{4} of the rare earth ion Er3+E_{r}^{3+} doped in calcium fluorophosphate, respectively.

IV Conclusion

In conclusion,we have demonstrated a scheme for realizing negative refraction without absorption via an incoherent pump field and a strong coherent field in the dense four-level atomic system.With the application of the incoherent pump field to manipulate the population distribution of each level and the variable strong coherent field to create quantum coherence,the magnetic response is amplified and the probe field propagates transparently in some frequency extents.Without the constraint condition of two equal transition frequencies responding to the probe field,the atomic system displays negative refraction with vanishing absorption and left-handedness. Therefore,our aim for searching the low-loss negative refraction is possible by choosing appropriate parameters in the scheme,in respect that the main applied limitation of the negative refractive materials is the large amount of dissipation and absorption.However,an excess coherent field intensity would increase the absorption near the resonance.This can also be obviously observed here.

Acknowledgments

The work is supported by the National Natural Science Foundation of China ( Grant No.60768001 and No.10464002 ).

References

  • (1) M.O.Scully,M.S.Zubairy, Quantum Optics(Cambridge University Press, Cambridge, 1997).
  • (2) V.G.Veselago,Sov.Phys.Usp. 10 (1968) 509-514 .
  • (3) J.B.Pendry,Phys.Rev.Lett. 85 (2000) 3966-3969.
  • (4) L.Chen,S.He,L.Shen,Phys.Rev.Lett. 92 (2004)107404.
  • (5) K.Aydin,I.Bulu,E.Ozbay,Appl.Phys.Lett 90 (2007) 254102
  • (6) P.R.Berman, Phys.Rev.E 66 (2002) 067603
  • (7) Y.P.Yang,J.P.Xu,H.Chen,and S.Y.Zhu,Phys.Rev.Lett. 100 (2008) 043601.
  • (8) V.Yannopapas,E.Paspalakis,and N.V.Vitanov,Phys.Rev.Lett. 103 (2008) 063602.
  • (9) Z.M.Zhang,C.J.Fu, Appl. Phys. Lett. 80 (2002) 1097-1099.
  • (10) R.A.Shelby,D.R.Smith,S.Schultz,Science 292 (2001) 77-79.
  • (11) J.Pendry,Nature 423 (2003) 22-23.
  • (12) E.Cubukcu,Nature 423 (2003) 604-605.
  • (13) G.V.Eleftheriades,A.K.Iyer,P.C.Kremer,IEEE Trans.Microwave Theory Tech. 50(2002)2702-2712.
  • (14) J.B.Pendry,Science 306 (2004) 1353-1355.
  • (15) V.Yannopapas,J. Phys.: Condens. Matter 18(2006) 6883-6890
  • (16) Q.Thommen,P.Mandel,Phys.Rev.Lett. 96(2006) 053601.
  • (17) M.o¨\ddot{o}.Oktel,o¨\ddot{o}.E.Mu¨\ddot{u}tecaplgˇ\check{g}u,Phys.Rev.A 70(2004)053806.
  • (18) C.S.Zhao,D.Z.Liu, Int.J.Quant.Inf. 7 (2009) 747-754.
  • (19) J.Q.Shen,Phys.Lett.A 357(2006) 54-60.
  • (20) J.Ka¨\ddot{a}stel,M.Fleischhauer,S.F.Yelin,R.L.Walsworth,Phys.Rev.Lett. 99(2007) 073602.
  • (21) J.Ka¨\ddot{a}stel,M.Fleischhauer,S.F.Yelin,R.L.Walsworth,Phys.Rev. A 79(2009) 063818.
  • (22) P.Tassin,L.Zhang,Th.Koschny,E.N.Economou,and C.M.Soukoulis,Phys.Rev.Lett. 102(2009) 053901.
  • (23) F.L. Li,A.P. Fang,M.Wang,J. Phys.B42(2009) 199505.
  • (24) H. J. Zhang,Y.P Niu,S. Q. Gong,Phys. Lett. A 363(2007) 497-501.
  • (25) G. S. Agarwal,R.W.Boyd,Phys.Rev.A 60(1999) R2681.
  • (26) J. Q. Shen,J.Mod.Opt.53 (2006) 2195-2205.
  • (27) A. A. Kaminskii, V. Mironov, S. A. Kornienko, S. N. Bagaev,G. Boulon, A. Brenier,B. Di Bartolo, Phys.Status Solidi A 151(1995) 231-255.
  • (28) D.K.Sardar,C.H.Coeckelenbergh,R.M.Yow,J.B.Gruber,T.H.Allik,J.Appl.Phys.98(2005) 033535.
  • (29) N.Papasimakis,V.A.Fedotov,N.I.Zheludev,S.L.ProsvirninPhys.Rev.Lett. 101(2008) 253903.