Local radiative density of optical states (LDOS) offers a tool to control the radiative rate of spontaneous emission from molecules, atoms, and quantum dots, which is proportional to LDOS. This paper presents that LDOS how to make the population of excited-state decay exponentially in time, and how these dynamics can be affected. By adopting the plane-wave expansion method, properties of an inverse-opal photonic crystal are studied with the help of photonic dispersion relations. Results in this paper show that the LDOS is radically modified in photonic crystal, and the rate of spontaneous emission can be described by the functions of position in the crystal and orientation of transition dipole moment.
Peng QIU, Guanglong WANG, Jianglei LU, Hongpei WANG. Local density of states in photonic crystal cavity[J]. Frontiers of Optoelectronics, 2012, 5(3): 341-344. DOI: 10.1007/s12200-012-0271-0
Introduction
The radiative rate of spontaneous emission of elementary light source (atoms, molecules or quantum dots) is proportional to local radiative density of optical states (LDOS), and sources with certain orientations of transition dipole moments couple to LDOS [1]. This projected LDOS represents the number of electromagnetic states at a given frequency, location and orientation of the transition dipole. Therefore, based on the frequency and position dependence of the LDOS, a tool is offered to control the radiative rate of spontaneous emission from molecules, atoms, and quantum dots [2].
Photonic crystal is optical structure, it contains a periodic modulation of their refractive index and form band gaps, which can prevent electromagnetic from propagation in the plane of the crystal. Photonic crystal with defects ruin the translational symmetry of crystal lattice, form a cavity surrounding the defect, and introduce an increased LDOS within the band gap [1]. In this paper, plane-wave expansion method was adopted to calculate the LDOS in photonic crystal cavity, and the effects of photonic crystal cavity on LDOS and radiative rate of spontaneous emission were investigated.
Local density of states
Spontaneous emission of light is called a transition of an emitter (atom, molecule or quantum dot) from its excited state to its lower energy state by emitting one or more photons. In the case of a two-level atom, ϵa and ϵb are the energies of the excited state a and ground state b respectively. And the interaction between the two-level atom in its excited state a with vacuum-field modes with frequency ωk leads to a transition to the atomic ground state b via emission of a photon with energy ћωk = ϵa-ϵb, as shown in Fig. 1.
Fig.1 Scheme of interaction of two-level atom with vacuum-field mode
The Schrodinger equation of the system as follows [3]:
The state is a linear combination of the states and of the unperturbed atom-field system. Here, is the state, in which the atom is under its excited state and the field has n photons. is the state, in which the atom is in its ground state and the field has n photons. Taking at time t= 0 the atom to be in the state and the field in vacuum mode , the state vector can be written aswhere ca(t) and cb,BoldItalic(t) are time dependent probability amplitudes of state, in which initial conditions ca(0) = 1 and cb,BoldItalic(0) = 0. From the Schrodinger equation, the equation of motions for ca(t) and cb,BoldItalic(t) can be obtained, with frequency integration over a Dirac δ function and further approximations, ca(BoldItalic,t) can be expressed as [4]
Equation (3) contains the LDOS, which counts the number of modes per unit volume at a given frequency ω, and the atomic dipole oriented along BoldItalicd and positioned at BoldItalic can couple to this frequency. This projected LDOS is defined as [5]where the summation on the field polarization p has been inserted back. The dielectric function ϵ(BoldItalic) must be real here; otherwise the complex mode density N(BoldItalic, ω, BoldItalicd) is not be defined.
Plane-wave expansion method
The plane-wave method is a widespread technique to calculate eigenvalues and the local density of eigenmodes in perfectly-periodic photonic crystal. According to Bloch theorem, which is the base of the plane-wave method, the Bloch modes of a photonic crystal can be expressed as a product of plane waves and functions describing the periodicity of the crystal lattice. These known periodic functions are expanded into Fourier series and inserted into the wave equation. To compute the eigenfrequencies ω(k) and the expansion coefficients of the eigenmodes , the infinite equation set is truncated: the reciprocal-lattice vectors BoldItalic are restricted to a finite set ζ with NBoldItalic elements[6,7]. This results in a 3NBoldItalic dimensional equation:
The transversality of the field gives an additional condition: , which eliminates one component of . For each BoldItalic + BoldItalic one needs to find two orthogonal unit vectors that form an orthogonal triad with BoldItalic + BoldItalic. By expressing the eigenmode expansion coefficients in the plane normal to BoldItalic + BoldItalic as , one third of the unknown is removed. Equation (5) becomes
The coefficients ηBoldItalic-BoldItalic′ are computed by first Fourier-transforming the dielectric function ϵ(BoldItalic), truncating and inverting the resulting matrix. Solving Eq. (6) gives the frequencies ωn(BoldItalic) and H-field eigenmodes BoldItalicn,BoldItalic(BoldItalic) in the photonic crystal. Then, using the Maxwell equations, the E-fields BoldItalicn,BoldItalic(BoldItalic) are obtained from
To compute the LDOS, the integration over BoldItalic needs to be discredited. The expression for the LDOS is then written aswhere the number of k-points should be as large as possible for much better calculation accuracy. In this paper, the LDOS histograms are plotted versus the reduced frequency a/λ= 2πc, where a is the cubic lattice parameter.
LDOS in photonic crystal cavity
In the photonic crystal cavity, the eigenmodes EnBoldItalic(BoldItalic) = ΛnBoldItalic(BoldItalic)/ are Bloch functions [8], so that the expression for the LDOS in Eq. (4) becomeswhere n is the band index. In this case, the dipole orientation BoldItalicd is averaged over all solid angles, which results in
Further, in order to describe emission dynamics of atoms randomly distributed in photonic crystal cavity, the LDOS can be integrated over the cavity:where ρ(BoldItalic) is a density of atoms at certain cavity in the crystal [9,10].
To examine the effects of photonic crystal on the LDOS and radiative dynamics of internal sources, an example of an inverse-opal photonic crystal shown in Fig. 2 will be discussed. LDOS calculated with the plane-wave method in inverse-opal photonic crystal is shown in Fig. 3. The LDOSs are at three different positions in the crystal: at point (0, 0, 0), (1/4, 1/4, 1/4) and (1/4, 1/4, 0). At the first two points, the LDOS does not depend on the dipole orientation BoldItalicd due to high symmetry, whereas at the third point, the LDOS is strongly orientation-dependent.
Fig.2 View of single unit cell of modeled inverse opal
At frequencies near the lowest-order stopgaps, a/λ = 0.6, the Bloch modes are eliminated for BoldItalic-vectors in a solid angle, which can occupy a considerable part of the whole 4π solid angle. This leads to decreased LDOS and spontaneous-emission rate, especially when the orientation of the transition dipole BoldItalicd and its position in the crystal are not favorable for coupling to allowed Bloch modes. At frequencies in the photonic bandgap, a/λ = 0.9, there are no Bloch modes at all: the LDOS is zero regardless the dipole orientation and position in the crystal. Consequently, spontaneous emission from a dipole emitter vanishes completely. The decay of the excited state in this case can only occur via possible non-radiative channels and via weaker atom-field interaction processes. On the other hand, at frequencies outside the bandgap, the density of Bloch modes is increased leading to enhanced radiative decay rates. Figure 3 reveals sharp peaks in the LDOS just above the band gap, which means that the strong modification of the LDOS in photonic bandgap material can lead to ultimate suppressions as well as enhancements of spontaneous emission.
Fig.3 LDOS calculated with plane-wave method in inverse-opal photonic crystal. LDOS at (a) (0, 0, 0) and (1/4, 1/4, 1/4) are shown by connected circles and triangles, respectively; (b) (1/4, 1/4, 0) projected on (1, 1, 0) and (0, 0, 1) directions are shown by connected squares and triangles, respectively
The calculating of LDOS based on plane-wave expansion method was introduced, and its application in an inverse-opal photonic crystal was discussed. Results indicate that the plane-wave expansion method is valid for the calculating of LDOS, and the LDOS can be radically modified in photonic crystal. In addition, it is also proved that in photonic crystal, the LDOS and the rate of spontaneous emission are the functions of position in the crystal unit-cell and orientation of the transition dipole moment. To conclude, LDOS is important for controlling spontaneous decay rates of elementary light sources.
Acknowledgements
This work was supported by the Innovation Foundation of Ordnance Engineering College (No. YSCX0905).
McPhedran R C, Botten L C, McOrist J, Asatryan A, de Sterke C, Nicorovici N. Density of states functions for photonic crystals. Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, 2004, 1(69): 16609
Fussell D P, McPhedran R C, Sterke C M. Three-dimensional green’s tenser, local density of states, and spontaneous emission in finite two-dimensional photonic crystals composed of clinders. Physical Review E, 2004, 2(70 ): 66608
Englund D, Fattal D, Waks E, Solomon G, Zhang B, Nakaoka T, Arakawa Y, Yamamoto Y, Vucković J. Controlling the spontaneous emission rate of single quantum dots in a two-dimensional photonic crystal. Physical Review Letters, 2005, 95(1): 013904
Lodahl P, Floris Van Driel A, Nikolaev I S, Irman A, Overgaag K, Vanmaekelbergh D, Vos W L. Controlling the dynamics of spontaneous emission from quantum dots by photonic crystals. Nature, 2004, 430(7000): 654–657
Sozuer H S, Haus J W, Inguva R. Photonic bands: convergence problems with the plane-wave method. Physical Review B: Condensed Matter and Materials Physics, 1992, 45(24): 13962–13972
Yamaguchi M, Asano T, Kojima K, Noda S. Quantum electrodynamics of a nanocavity coupled with exciton complexes in a quantum dot. Physical Review Letters, 2009, 15(80): 1–10
8
Zhou W D, Qiang Z X, Chen L. Photonic crystal defect mode cavity modelling: a phenomenological dimensional reduction approach. Applied Physics (Berlin), 2007, 40(9): 2615–2623
9
Wu F M, Meng Y L, Xie G H. Chen P, Zhao Y. Top-emitting white organic light-emitting devices based on microcavity structure. Journal of Optoelectronics·Laser, 2008, 19(10): 1287–1290
10
Wang D H, Huang K Y, Xu Q. Spontaneous emission of a polarized atom in a medium between two parallel mirrors. Chinese Physical Society, 2010, 1(15): 138–144