OSA's Digital Library

Optics Express

Optics Express

  • Editor: J. H. Eberly
  • Vol. 3, Iss. 1 — Jul. 6, 1998
  • pp: 12–18
« Show journal navigation

Novel approach to photonic bands with frequency-dependent dielectric constants

Kazuaki Sakoda and Jun Kawamata  »View Author Affiliations


Optics Express, Vol. 3, Issue 1, pp. 12-18 (1998)
http://dx.doi.org/10.1364/OE.3.000012


View Full Text Article

Acrobat PDF (362 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We formulated a novel method to calculate the dispersion relations of arbitrary photonic crystals with frequency-dependent dielectric constants based on the numerical simulation of dipole radiation. As an example, we applied this method to a two-dimensional square lattice of metallic cylinders and obtained a good agreement with the previous result by means of the plane-wave expansion method by Kuzmiak et al. [Phys. Rev. B 50, 16 835 (1994)]. In addition to the dispersion relations, we could obtain the symmetries and the wave functions of the eigenmodes.

© Optical Society of America

1. Introduction

Photonic crystals, which have spatially periodic dielectric structures, can control the nature of radiation field.[1–3

1. J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals, (Princeton University Press, Princeton, 1995).

] They can create photonic band gaps, or the frequency ranges in which no eigenmodes of radiation field exist. Then the fundamental optical properties of atoms and molecules such as the rate of spontaneous emission may change drastically when they are embedded in the photonic crystals. In addition, localized eigenmodes of radiation field may appear in the photonic band gaps when appropriate defects are introduced to the regular dielectric structures of the photonic crystals.[4

4. S. L. McCall, P. M. Platzman, R. Dalichaouch, D. Smith, and S. Schultz, “Microwave propagation in two-dimensional dielectric lattices,” Phys. Rev. Lett. 67, 2017–2020 (1991). [CrossRef] [PubMed]

, 5

5. E. Yablonovitch, T. J. Gmitter, R. D. Meade, A. M. Rappe, K. D. Brommer, and J. D. Joannopoulos, “Donor and acceptor modes in photonic band structure,” Phys. Rev. Lett. 67, 3380–3383 (1991). [CrossRef] [PubMed]

] Because the localized modes have high quality factors and show extremely narrow transmission spectra, they can be used as good resonators and optical filters. Moreover, as was pointed out by Yablonovitch,[6

6. E. Yablonovitch, “Photonic band-gap structures,” J. Opt. Soc. Am. B 10, 283–295 (1993). [CrossRef]

] the single-mode light-emitting diodes that utilize the emission through the localized modes may have such properties as good temporal and spatial coherence, low noise, high efficiency, and high modulation rate. Then the micro-fabrication of the photonic crystals with semiconducting materials have been tried energetically.

When we analyze and design the photonic crystals made from semiconductors that work in visible and near infrared frequency ranges, we have to take into account the large frequency dependence of their dielectric constants. However, most of the theoretical analyses of the photonic band structures have been made for the frequency-independent case. An exception is the work reported by Kuzmiak et al.[7

7. V. Kuzmiak, A. A. Maradudin, and F. Pincemin, “Photonic band structures of two-dimensional systems containing metallic components,” Phys. Rev. B 50, 16 835–16 844 (1994). [CrossRef]

] They calculated the photonic band structure, or the dispersion relation, of a two-dimensional (2D) array of metallic cylinders with a frequency-dependent dielectric constant of Drude type by means of the plane-wave expansion method. They showed that the wave equation derived from the Maxwell’s equations leads to a well-defined eigenvalue problem for this particular case when the electric field of the eigenmode is parallel to the cylinder axis (E polarization) and the plane-wave expansion method gives a good convergence. But we have to devise another method when the wave equation does not lead to a well-defined eigenvalue problem, and then the good convergence of the plane-wave expansion is not necessarily guaranteed.

Then, in this paper, we would like to present a new and general method to calculate the dispersion relations of the photonic crystals with frequency-dependent dielectric constants. Our method is based on the numerical simulation of the radiation process of a virtual oscillating dipole moment embedded in the photonic crystal. Our algorithm is suitable to vector processing and the CPU time can be reduced by a considerable amount by using a vector processor. This method was first formulated and applied to the problem of mid-gap eigenmodes localized at point [8

8. K. Sakoda and H. Shiroma, “Numerical method for localized defect modes in photonic lattices,” Phys. Rev. B 56, 4830–4835 (1997). [CrossRef]

] and line defects.[9

9. K. Sakoda, T. Ueta, and K. Ohtaka, “Numerical analysis of eigenmodes localized at line defects in photonic lattices,” Phys. Rev. B 56, 14 905–14 908 (1997). [CrossRef]

] In Sec. 2, we will extend this method in order to treat frequency-dependent dielectric constants and to calculate the dispersion relations. Then it will be applied to the particular problem of Kuzmiak et al, [7

7. V. Kuzmiak, A. A. Maradudin, and F. Pincemin, “Photonic band structures of two-dimensional systems containing metallic components,” Phys. Rev. B 50, 16 835–16 844 (1994). [CrossRef]

] i.e., the photonic band structure of the 2D square array of metallic cylinders for E polarization. The numerical results will be presented and compared with those of Kuzmiak et al in Sec. 3.

2. Theory

First, let us summarize our method briefly. Originally we considered the radiation process of the virtual oscillating dipole moment P d(r, t) that was embedded in the photonic crystal with a structural defect:

Pdrt=μδ(rr0)exp(iωt),
(1)

where μ, and r 0 are the magnitude and the position of the dipole moment, respectively, ω is the angular frequency of the oscillation, and δ is Dirac’s delta function. Then we showed by means of a Green’s-function method[10

10. K. Sakoda and K. Ohtaka, “Optical response of three-dimensional photonic lattices: Solutions of inhomogeneous Maxwell’s equations and their applications,” Phys. Rev. B 54, 5732–5741 (1996). [CrossRef]

] that the induced electric field, E(r, t), and the electromagnetic energy radiated in a unit time, U, are given as follows.[8

8. K. Sakoda and H. Shiroma, “Numerical method for localized defect modes in photonic lattices,” Phys. Rev. B 56, 4830–4835 (1997). [CrossRef]

]

Ert2πωd{μ·Ed*(r0)}Ed(r)exp(iωt)V(ωωd+),
(2)
Uπωd2γμ·Ed(r0)2V{(ωωd)2+γ2},
(3)

where ωd and E d(r) denote the eigen-angular frequency and the eigenfunction of the relevant defect mode, γ is its decay rate, V is the volume of the photonic crystal. E d(r) is normalized as follows.

Vϵ(r)Ed(r)2dr=V.
(4)

Here ϵ(r) is the position-dependent dielectric constant of the photonic crystal. When we derived Eqs. (2) and (3), we assumed that ω was close to ωd that was isolated in the photonic band gap and neglected the contribution from all other eigenmodes. If we can evaluate the frequency dependence of U, we can obtain the eigen-angular frequency as the resonance frequency [see Eq. (3)], and then the induced electric field is proportional to the eigenfunction [see Eq. (2)]. Then we discretized the wave equation derived from the Maxwell’s equations to obtain a difference equation and solved it numerically. We could obtain ωd as the resonance frequency as we had expected, and the calculated values were very close to observed ones. [8

8. K. Sakoda and H. Shiroma, “Numerical method for localized defect modes in photonic lattices,” Phys. Rev. B 56, 4830–4835 (1997). [CrossRef]

, 9

9. K. Sakoda, T. Ueta, and K. Ohtaka, “Numerical analysis of eigenmodes localized at line defects in photonic lattices,” Phys. Rev. B 56, 14 905–14 908 (1997). [CrossRef]

]

Now, we extend this method in order to treat the frequency dependence of the dielectric constants and to calculate the dispersion relations. We begin with the following two equations.

×Ert=1ctHrt,
(5)
×Hrt=1ct{D0rt+4πPdrt},
(6)

where c is the light velocity in vacuum. In Eq. (6), P d(r,t) stands for the virtual oscillating dipole moment as before and D 0(r,t) denotes the electric displacement due to the regular dielectric structure of the photonic crystal. The latter is generally given by the convolution integral of the electric field E(r, t) and the dielectric response function Φ(r,t):

D0rt=dtΦrttErt.
(7)

Φ(r, t) is given by the Fourier transform of the dielectric constant ϵ(r, ω), which is now a function of the frequency as well as the spatial coordinates:

Φrt=12πdωϵrωexp(iωt).
(8)

In order to show the detail of our numerical method, let us treat the particular problem discussed by Kuzmiak et al.[7

7. V. Kuzmiak, A. A. Maradudin, and F. Pincemin, “Photonic band structures of two-dimensional systems containing metallic components,” Phys. Rev. B 50, 16 835–16 844 (1994). [CrossRef]

] in what follows. So, we consider a photonic crystal that contains metallic components and assume a dielectric constant of Drude type in the metal:

ϵrω=ϵ[1ωp2τω(ωτ+i)],
(9)

where ωp is the plasma frequency, τ is the relaxation time, and ϵ is the dielectric constant at sufficiently high frequencies. In Eq. (9), we took into account the imaginary part of the dielectric constant in order to fulfill the Kramers-Kronig relation, and hence, the causality. Then Eq. (8) leads to

Φrt=ϵδ(t)+ϵωp2τ[1exp(tτ)]θ(t),
(10)

where θ(t) is a step function. From Eqs. (5), (6), (7), and (10), we obtain the following wave equation for the metallic region.

×{×Ert}
=1c2[ϵ2t2Ert+ϵωp2Ertϵωp2τ0dtexp(tτ)Ertt
4πω2μδ(rr0)exp(iωt)].
(11)

On the other hand, we assume that ϵ(r, ω) = 1 (air) outside the metal. Then for this region,

×{×Ert}=1c2[2t2Ert4πω2μδ(rr0)exp(iωt)].
(12)

We discretized Eqs. (11) and (12) to obtain difference equations and solved them numerically with an initial condition E(r, 0) = 0 and a boundary condition

E(r+a,t)=exp(ik·a)Ert,
(13)

where k is a wave vector in the first Brillouin zone and a is the elementary lattice vector. The latter condition extracts the contribution to the radiated electromagnetic field from particular eigenmodes with the specified wave vector. Therefore, we can calculate the resonance frequency as a function of k, i.e., we can obtain the dispersion relations. According to Kuzmiak et al.,[7

7. V. Kuzmiak, A. A. Maradudin, and F. Pincemin, “Photonic band structures of two-dimensional systems containing metallic components,” Phys. Rev. B 50, 16 835–16 844 (1994). [CrossRef]

] we considered a 2D photonic crystal composed of a square array of metallic cylinders with a radius R for the numerical calculation. The following parameters were assumed: R/a = 0.472, ϵ = 1.0, and ωp a/27πc = 1.0, where a denotes the lattice constant. Because of the boundary condition, Eq. (13), it was enough to deal with only one unit cell, and therefore, the CPU time necessary for the numerical calculation was small. In the actual calculation, the unit cell was divided into 40 × 40 parts, and one period of the oscillation was divided into 160 parts in order to discretize the wave equation. The further decrease of the size of the spatial and temporal meshes did not bring about an apparent change in the resonance frequencies.

3. Results and discussion

The numerical calculation was done for the wave vectors in the (1, 0) direction in the first Brillouin zone of the 2D square lattice. As an example, Fig. 1 shows the electromagnetic energy radiated by the oscillating dipole moment as a function of the oscillation frequency for k = 0. The abscissa represents the normalized frequency. □, ○, ◇, and • denote the accumulated electromagnetic energy after 10, 20, 50, and 100 cycles of the oscillation, respectively. A resonance at ωa/2πc = 0.745 is clearly observed in this figure.

Figure 1. The electromagnetic energy radiated by an oscillating dipole moment located at r 0/a = (0.3,0) as a function of the oscillation frequency for k = 0. The electric field parallel to the cylinder axis (E polarization) was calculated. □, ○, ◇, and • denote the accumulated electromagnetic energy after 10, 20, 50, and 100 cycles of the oscillation, respectively. The abscissa represents the normalized frequency. According to the previous calculation by Kuzmiak et al. (Ref. [7]), the following parameters were used: R/a = 0.472, ω = 1.0, ϵ ωp a/2πc = 1.0. We also assumed that τ-1 = 0.01 × ωp. A resonance at ωa/2πc = 0.745 is clearly observed.
Figure 2. The dispersion relations in the (1, 0) direction (or the Γ-X direction) of the 2D photonic crystal composed of the metallic cylinders. The ordinate is the normalized frequency and the abscissa is the normalized wave vactor. The same parameters as for Fig. 1 were used.

The dispersion relations thus obtained are presented in Fig. 2 where the symmetry of each eigenmode is also shown. In this figure, the number in parentheses is given in order of the ascending frequency when there are more than one modes of the same symmetry in the analysed spectral region. Note that there is no eigenmode for ωa/2πc < 0.745. We can show that this cut-off frequency is consistent with the result of the long-wavelength approximation of the Maxwell’s equations. When we compare this figure with Fig. 1(b) of Ref. [7

7. V. Kuzmiak, A. A. Maradudin, and F. Pincemin, “Photonic band structures of two-dimensional systems containing metallic components,” Phys. Rev. B 50, 16 835–16 844 (1994). [CrossRef]

], it is found that the overall feature is quite common to both calculations except the degeneracy of the fourth and the fifth bands at the Γ point. As will be shown later, the symmetry of the eigenfunctions of these two modes is really that of the E representation of the C 4v point group that is doubly degenerate. Therefore, we believe that our calculation is more accurate concerned with this matter.

The distribution of the electric fields of the eigenmodes at the Γ point is presented in Fig. 3. The maximum of each electric field is normalized to unity. It is very clear that each eigenmode shows its peculiar symmetry. Especially, Fig. 3(d) is a replica of Fig. 3(e) given by a 90° rotation, which implies that they are degenerate and attributed to the E representation. Finally, Fig. 4 shows the field distribution of the eigenmodes at the X point. They are consistent with the symmetry assignments given in Fig. 2.

Figure 3. The distribution of the electric fields of the (a) A 1(1), (b) A 1(2), (c) B 1, and (d) (e) E modes at the Γ point that has the C 4v symmetry. The maximum of each electric field is normalized to unity. For all eigenmodes, the eigenfunctions show their peculiar symmetries. For doubly degenerate E mode, one eigenfunction is a replica of the other given by a 90° rotation.

We would like to conclude this section by making three remarks. First, we have treated the frequency-dependent dielectric constant of a particular type in this paper. But it is easy to deal with other types of frequency dependence which can be understood from the formulation given in Sec. 2. Second, we can treat three-dimensional (3D) photonic crystals as well as 2D ones. In that case, the spatial mesh for the discretization is 3D and the number of the data that should be dealt with in the numerical calculation increases by a considerable amount. But we believe that the CPU time is still allowable. Thirdly, we did not calculate the decay rate γ of the eigenmodes in this paper because we were mainly interested in the evaluation of the accuracy and efficiency of the present method as a tool for the photonic band calculation. γ is inevitably non-zero when the dielectric constant depends on frequency and has the non-zero imaginary part as a result of the Kramers-Kronig relation. γ can be evaluated by analysing the width of the resonance curve, or more practically by calculating the temporal evolution of the accumulated electromagnetic energy after switching off the oscillation of the dipole.

Figure 4. The distribution of the electric fields of the (a) B 1(1), (b) A 1(1), (c) A 2, (d) B 1(2), (e) A 1(2), and (f) B 2 modes at the X point that has the C 2v symmetry. The maximum of each electric field is normalized to unity as before.

4. Conclusion

We formulated a new and general method to calculate the dispersion relations of arbitrary photonic crystals with frequency-dependent dielectric constants based on the numerical simulation of the radiation process of a virtual oscillating dipole moment embedded in the crystals. As an example, we apllied the present method to the 2D square lattice composed of metallic cylinders with the dielectric constant of Drude type. In addition to the dispersion relations that were consistent with the symmetry of the crystal, we could calculate the wave functions of the eigenmodes.

Acknowledgments

This work was financially supported by the Sumitomo Foundation and the Research Foundation for Opto-Science and Technology.

References

1.

J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals, (Princeton University Press, Princeton, 1995).

2.

Photonic Band Gaps and Localization, edited by C. M. Soukoulis (Plenum, New York, 1993).

3.

Photonic Band Gap Materials, edited by C. M. Soukoulis (Kluwer, Dordrecht, 1996).

4.

S. L. McCall, P. M. Platzman, R. Dalichaouch, D. Smith, and S. Schultz, “Microwave propagation in two-dimensional dielectric lattices,” Phys. Rev. Lett. 67, 2017–2020 (1991). [CrossRef] [PubMed]

5.

E. Yablonovitch, T. J. Gmitter, R. D. Meade, A. M. Rappe, K. D. Brommer, and J. D. Joannopoulos, “Donor and acceptor modes in photonic band structure,” Phys. Rev. Lett. 67, 3380–3383 (1991). [CrossRef] [PubMed]

6.

E. Yablonovitch, “Photonic band-gap structures,” J. Opt. Soc. Am. B 10, 283–295 (1993). [CrossRef]

7.

V. Kuzmiak, A. A. Maradudin, and F. Pincemin, “Photonic band structures of two-dimensional systems containing metallic components,” Phys. Rev. B 50, 16 835–16 844 (1994). [CrossRef]

8.

K. Sakoda and H. Shiroma, “Numerical method for localized defect modes in photonic lattices,” Phys. Rev. B 56, 4830–4835 (1997). [CrossRef]

9.

K. Sakoda, T. Ueta, and K. Ohtaka, “Numerical analysis of eigenmodes localized at line defects in photonic lattices,” Phys. Rev. B 56, 14 905–14 908 (1997). [CrossRef]

10.

K. Sakoda and K. Ohtaka, “Optical response of three-dimensional photonic lattices: Solutions of inhomogeneous Maxwell’s equations and their applications,” Phys. Rev. B 54, 5732–5741 (1996). [CrossRef]

OCIS Codes
(260.2110) Physical optics : Electromagnetic optics
(350.3950) Other areas of optics : Micro-optics

ToC Category:
Focus Issue: Photonic crystals

History
Original Manuscript: March 23, 1998
Revised Manuscript: February 16, 1998
Published: July 6, 1998

Citation
Kazuaki Sakoda and Jun Kawamata, "Novel approach to photonic bands with frequency-dependent dielectric constants," Opt. Express 3, 12-18 (1998)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-3-1-12


Sort:  Journal  |  Reset  

References

  1. J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals, (Princeton University Press, Princeton, 1995).
  2. Photonic Band Gaps and Localization, edited by C. M. Soukoulis (Plenum, New York, 1993).
  3. Photonic Band Gap Materials, edited by C. M. Soukoulis (Kluwer, Dordrecht, 1996).
  4. S. L. McCall, P. M. Platzman, R. Dalichaouch, D. Smith, and S. Schultz, "Microwave propagation in two-dimensional dielectric lattices," Phys. Rev. Lett. 67, 2017-2020 (1991). [CrossRef] [PubMed]
  5. E. Yablonovitch, T. J. Gmitter, R. D. Meade, A. M. Rappe, K. D. Brommer, and J. D. Joannopoulos, "Donor and acceptor modes in photonic band structure," Phys. Rev. Lett. 67, 3380-3383 (1991). [CrossRef] [PubMed]
  6. E. Yablonovitch, "Photonic band-gap structures," J. Opt. Soc. Am. B 10, 283-295 (1993). [CrossRef]
  7. V. Kuzmiak, A. A. Maradudin, and F. Pincemin, "Photonic band structures of two-dimensional systems containing metallic components," Phys. Rev. B 50, 16 835-16 844 (1994). [CrossRef]
  8. K. Sakoda and H. Shiroma, "Numerical method for localized defect modes in photonic lattices," Phys. Rev. B 56, 4830-4835 (1997). [CrossRef]
  9. K. Sakoda, T. Ueta, and K. Ohtaka, "Numerical analysis of eigenmodes localized at line defects in photonic lattices," Phys. Rev. B 56, 14 905-14 908 (1997). [CrossRef]
  10. K. Sakoda and K. Ohtaka, "Optical response of three-dimensional photonic lattices: Solutions of inhomogeneous Maxwell's equations and their applications," Phys. Rev. B 54, 5732-5741 (1996). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited