OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 11, Iss. 6 — Mar. 24, 2003
  • pp: 541–551
« Show journal navigation

Moving least-square method for the band-structure calculation of 2D photonic crystals

Sukky Jun, Young-Sam Cho, and Seyoung Im  »View Author Affiliations


Optics Express, Vol. 11, Issue 6, pp. 541-551 (2003)
http://dx.doi.org/10.1364/OE.11.000541


View Full Text Article

Acrobat PDF (645 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

The moving least-square (MLS) basis is implemented for the real-space band-structure calculation of 2D photonic crystals. A value-periodic MLS shape function is thus proposed in order to represent the periodicity of crystal lattice. Through numerical examples, this MLS method is proved to be a promising scheme for predicting band gaps of photonic crystals.

© 2003 Optical Society of America

1. Introduction

The photonic band gap material is nowadays one of central issues in photophysics and optoelectronics communities [1

1. J.D. Joannopoulos, R.D. Meade, and J.N. Winn, Photonic Crystals: Molding the Flow of Light (Princeton University Press, Princeton, 1995).

, 2

2. Y. Xia, “Photonic crystals,” Adv. Mater. 13, 369 (2001) and papers in this special issue. [CrossRef]

]. Accurate computation of the band structure is thus indispensable for the development of various optoelectronic devices and optical fibers based on these materials [3

3. K. Busch, “Photonic band structure theory: assesment and perspectives,” C. R. Physique 3, 53–66 (2002). [CrossRef]

, 4

4. D. Cassagne, “Photonic band gap materials,” Ann. Phys. Fr. 23(4), 1–91 (1998). [CrossRef]

, 5

5. J.B. Pendry, “Calculating photonic band structure,” J. Phys.: Condens. Matter 8, 1085–1108 (1996). [CrossRef]

]. It is well known that the conventional plane wave method has convergence problem arising from the abrupt change in the value of dielectric function across the interface between matrix and inclusion [4

4. D. Cassagne, “Photonic band gap materials,” Ann. Phys. Fr. 23(4), 1–91 (1998). [CrossRef]

, 6

6. H.S. Sözüer, J.W. Haus, and R. Inguva, “Photonic bands: Convergence problems with the plane-wave method,” Phys. Rev. B 45, 13962–13972 (1992). [CrossRef]

]. Special techniques are therefore required for the plane wave method to achieve faster convergence, such as interpolating dielectric function [7

7. R.D. Meade, A.M. Rappe, K.D. Brommer, J.D. Joannopoulos, and O.L. Alerhand, “Accurate theoretical analysis of photonic band-gap materials,” Phys. Rev. B 48, 8434–8437 (1993). [CrossRef]

] for example. Other numerical methods have also been tried for more efficient calculations of photonic band structures. Among them are the finite-difference time-domain (FDTD) method [8

8. C.T. Chan, Q.L. Yu, and K.M. Ho, “Order-N spectral method for electromagnetic waves,” Phys. Rev. B 51, 16635–16642 (1995). [CrossRef]

, 9

9. A.J. Ward and J.B. Pendry, “Calculating photonic Green’s functions using a nonorthogonal finite-difference time-domain method,” Phys. Rev. B 58, 7252–7259 (1998). [CrossRef]

, 10

10. M. Qiu and S. He, “A nonorthogonal finite-difference time-domain method for computing the band structure of a two-dimensional photonic crystal with dielectric and metallic inclusions,” J. Appl. Phys. 87, 8268–8275 (2000). [CrossRef]

], the multiple-scattering theory [11

11. K.M. Leung and Y. Qiu, “Multiple-scattering calculation of the two-dimensional photonic band structure,” Phys. Rev. B 48, 7767–7771 (1993). [CrossRef]

, 12

12. X. Wang, X.G. Zhang, Q. Yu, and B.N. Harmon, “Multiple-scattering theory for electromagnetic waves,” Phys. Rev. B 47, 4161–4167 (1993). [CrossRef]

], the transfer matrix method [13

13. J.B. Pendry and A. MacKinnon, “Calculation of photon dispersion relations,” Phys. Rev. Lett. 69, 2772–2775 (1992). [CrossRef] [PubMed]

], the finite difference method with effective medium theory [14

14. L. Shen, S. He, and S. Xiao, “A finite-diference eigenvalue algorithm for calculating the band structure of a photonic crystal,” Comput. Phys. Comm. 143, 213–221 (2002). [CrossRef]

], the finite element method [15

15. W. Axmann and P. Kuchment, “An efficient finite element method for computing spectra of photonic and acoustic band-gap materials: I. Scalar case,” J. Comput. Phys. 150, 468–481 (1999). [CrossRef]

, 16

16. D.C. Dobson, “An efficient band structure calculations in 2D photonic crystals,” J. Comput. Phys. 149, 363–376 (1999). [CrossRef]

, 17

17. C. Mias, J.P. Webb, and R.L. Ferrari, “Finite element modelling of electromagnetic waves in doubly and triply periodic structures,” IEE Proc.-Optoelectron. , 146(2), 111–118 (1999). [CrossRef]

], and most recently, the cell method [18

18. M. Marrone, V.F. Rodriguez-Esquerre, and H.E. Hernandez-Figueroa, “Novel numerical method for the analysis of 2D photonic crystals: the cell method,” Opt. Express 10, 1299–1304 (2002), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-10-22-1299 [CrossRef] [PubMed]

].

In this paper, we apply the novel moving least-square (MLS) method to the band-structure computation of 2D photonic crystals. The MLS method is also known as the meshfree method in the communities of solid and fluid mechanics. It is a real-space scheme, and thus the discontinuity of dielectric function can naturally be modelled without requiring the Fourier transform. In addition, any irregular shape of interface between materials is easily recognized in MLS modelling by simply putting nodes along the interface, which is similar to finite element method. On the other hand, it is a point-based method that does not need grid or mesh in MLS approximation procedure, being different from other real-space techniques. Owing to its distinctive features, the MLS method has widely been applied to various kinds of engineering problems (for the most recent review, refer to [19

19. S. Li and W.K. Liu, “Meshfree and particle methods and their applications,” Applied Mechanics Review , 55, 1–34 (2002). [CrossRef]

]). Specifically the method may be well-suited to the adaptive computation, which motivates us to investigate here the feasibility of this method for the band-structure calculation of photonic crystals. While executing the photonic band calculation, the most important issue is how to realize the crystal’s periodicity in the MLS modelling procedure. To this end, we have developed an algorithm introducing a value-periodic MLS shape function. Therefore, the goal of this paper is to demonstrate that this periodic MLS shape function can be applied to highly accurate computations of photonic band structure. We show that its performance is outstanding well beyond the conventional plane wave method throughout numerical examples.

2. The moving least-square (MLS) approximation

The MLS representation of an arbitrary function begins with discretizing a domain by a set of nodes. Nodal value u J is assigned to each node. The approximation u h(x) of a function u(x) at any point x in the domain of computation is then expressed by

uh(x)=JNPNJ(x)uJ
(1)

where NP is the total number of nodes. By minimizing the localized error residual functional that is expressed by the moving least-square procedure (see Appendix or [19

19. S. Li and W.K. Liu, “Meshfree and particle methods and their applications,” Applied Mechanics Review , 55, 1–34 (2002). [CrossRef]

, 20

20. D.W. Kim and Y. Kim, “Point collocation method using the fast moving least-square reproducing kernel approximation,” Int. J. Numer. Methods Engrg. 56, 1445–1464 (2003). [CrossRef]

] for details), the MLS shape function NJ (x) is defined as

NJ(x)=N(xJx)pT(xJx)b(x)W(xJx)
(2)

where p(xJ-x) is the polynomial basis vector and W(xJ-x) the window function. They are here abbreviated as pJ(x)=p(xJ-x) and W J (x)=W(xJ-x). The vector b(x) in Eq. (2) is determined by solving the matrix equation of

M(x)b(x)=p(0)
(3)

where p(0)=[1, 0, …, 0]T and the moment matrix M(x) is defined by

M(x)J=1NPWJ(x)pJ(x)pJT(x)
(4)

The MLS shape function is highly localized. That is, it vanishes in the most of problem domain, except a small region (so called compact support) of which the center is at the position x. The comprehensive description of MLS-based meshfree methods is found in a recent review paper [19

19. S. Li and W.K. Liu, “Meshfree and particle methods and their applications,” Applied Mechanics Review , 55, 1–34 (2002). [CrossRef]

] and references therein.

3. Galerkin formulation and matrix eigen-equations

For the computation of band structures in 2D photonic crystals, we begin with the two scalar equations for solving the unknown field intensity ψ(x) as below,

·1ε(x)ψ(x)=λψ(x)(TE modes)
(5)
1ε(x)2ψ(x)=λψ(x)(TM modes)
(6)

where λ=(ω/c)2 and ε(x) is the dielectric function. ω is the frequency of the monochromatic electromagnetic wave, c the speed of light, and x the position vector in x-y plane. Applying the Floquet-Bloch theorem, the problems are reduced to solving for u(x) as

(+ik)·1ε(x)(+ik)u(x)=λu(x)(TE modes)
(7)
(+ik)·(+ik)u(x)=λε(x)u(x)(TM modes)
(8)

where k is the wave vector and u(x) is the function fulfilling the periodicity of lattice structure of which the lattice vector is L, i.e., u(x)=u(x+L).

The Galerkin formulations for Eqs. (7) and (8) result in matrix eigen-equations as,

Au=λBu
(9)

where, for TE modes,

AIJ=Ω1ε(x)(+ik)NIP(x)·(+ik)NJP(x)¯dV,BIJ=ΩNIP(x)NJP(x)¯dV
(10)
Fig. 1. The translation-and-searching algorithm for constructing a periodic meshfree shape function. Supports of the shape function are in gray within solid circle. The dashed circle illustrates the rod of photonic crystal.

and, for TM modes,

AIJ=Ω(+ik)NIP(x)·(+ik)NJP(x)¯dV,BIJ=Ωε(x)NIP(x)NJP(x)¯dV
(11)

When we discretize the unit cell by a set of nodes (or points), it is geometrically natural to locate nodes exactly on the interface between matrix and inclusion. Unfortunately in this situation, it is not clear which dielectric constant has to be assigned to the interface nodes, because the matrix and inclusion have different dielectric constants from each other. In contrast, the MLS-based meshfree method usually employs a set of integration points, apart from the set of nodes, to calculate the integrals in matrix Eqs. [19]. And the integration points are located in either inclusion or matrix, not on the interface. Therefore, the given value of dielectric constant is straightforwardly assigned to integration points. It makes our implementation very simple, not requiring additional techniques such as, for instance, averaging or smoothing the dielectric function.

4. Value-periodic MLS shape functions

The description is now generalized to an arbitrary point x in the unit cell. Using the lattice vector L=n 1a1+n 2a2, the MLS approximation can be written in general form as

uh(x)=LJ=1NPNJ(x+L)uJ=J=1NP[LNJ(x+L)]uJ=J=1NPNJP(x)uJ
(12)

where N J (x+L) is defined by

NJ(x+L)=pJT(x+L)b(x)WJ(x+L)
(13)

Note that b is a function of x, not of x+L. The vector b(x) is obtained by solving

M(x)b(x)=p(0),whereM(x)LJ=1NPWJ(x+L)pJ(x+L)pJT(x+L)
(14)

The summation with respect to L is performed over the current unit cell itself and the neighboring cells only, because, for any farther L, the corresponding support does not overlap the unit cell of computation. In conclusion, the 2D value-periodic MLS shape function for a parallelogram unit cell whose primitive vectors are a1 and a2, can be written as

NJP(x)LNJ(x+L)=n1,n2NJ(x+n1a1+n2a2)
(15)

Therefore the summation over n i’s results in 9 cells in 2D (or 27 cells for 3D case) at most. That is, n i=-1, 0, +1 in Eq. (15) is enough for any x in the unit cell. It should be carefully noted that the opposing boundaries are physically identical, i.e., Γ13 and Γ24 as illustrated in Fig. 1. Therefore, for the nodal summation over J=1, …, NP, we must involve one boundary only among each identical pair in order to avoid over-summing. In other word, the nodes on Γ3 and Γ4 have to make no contribution to the summation with respect to J, once the corresponding nodes along Γ1 and Γ2 are summed up. An example of 2D periodic MLS shape function is plotted in Fig. 2 for a parallelogram unit cell.

5. Numerical examples

Fig. 2. The periodic MLS shape function. The position arrowed is the center of the shape function.

5.1 Electromagnetic Kronig-Penney problem

The electromagnetic Kronig-Penney problem is now considered. The unit cell is depicted in Fig. 3, which is composed of two different dielectric materials. The analytic solution of this case can be obtained elsewhere [17

17. C. Mias, J.P. Webb, and R.L. Ferrari, “Finite element modelling of electromagnetic waves in doubly and triply periodic structures,” IEE Proc.-Optoelectron. , 146(2), 111–118 (1999). [CrossRef]

]. This example is taken into consideration in order to investigate how accurately the MLS method handles the discontinuity of dielectric function. Dielectric constants of ε1=1.0 and ε2=9.0 are used for the simulation. Frequency band structures of both TM and TE modes are given in Fig. 4. Lines are for the analytic solution and open circles for the MLS results using 11×11 nodes for TM modes and 41×41 nodes for TE modes. Convergence rates of the lowest five eigenvalues at point of k=(π/a)[0.5, 0.5] are given in Fig. 5 in which En ’s denote the normalized frequency eigenvalues (ωn a/2πc) and Ee is the analytic solution. The average slope of TM and TE modes are 3.12 and 2.77 respectively. The results for TM modes show better performance over the TE modes case, not only in convergence but also in accuracy. This is due to the difference of smoothness properties of eigenvectors [16

16. D.C. Dobson, “An efficient band structure calculations in 2D photonic crystals,” J. Comput. Phys. 149, 363–376 (1999). [CrossRef]

]. There may be several ways to enhance the accuracy of MLS results on TE modes. The meshfree treatment of discontinuity using Lagrange multiplier has been applied to solid mechanics by [21

21. L.W. Cordes and B. Moran, “Treatment of material discontinuity in the Element-Free Galerkin method,” Comput. Meth. Appl. Mech. Eng. 139, 75–89 (1996). [CrossRef]

]. Its application to photonic crystals remains for our future study.

5.2 Square lattice of circular rods

Photonic band structures of square lattice comprised of circular rods are computed in this section. The square lattice structure and the square unit cell are illustrated in Fig. 6. The unit cell in the figure is discretized by 1697 nodes. The radius of circular cross section is r=0.2a for this example. Dielectric constants are εm=1.0 and εr=8.9 for matrix and rods respectively. In this case, there exists a wide band gap in TM modes as widely known in literatures. Our results on band structures also confirm the wide band gap in TM modes as shown in Fig. 7. In this example, the width of band gap is determined by the 1st TM mode at point M and the 2nd TM mode at point X. The numerical results at these points obtained by both the MLS method and plane wave method are compared in Fig. 8, where the faster converging behavior of the method over the plane wave method are demonstrated. In the case of MLS method, the x-axis indicates the size of discretization such as h=a/(NP -1), where NP denotes the number of nodes in one direction. For the plane wave method, it similarly implies h=a/(NPW -1), where NPW means the number of plane wave used for the computation. This notation is arbitrary because the MLS method solves generalized eigenvalue problems, while the plane wave method does standard eigenvalue problems. The MLS method thus requires more computing time than the conventional plane wave method in general. Nonetheless, it is convenient for the purpose of direct comparison as in Fig. 8, because both (NP)2 and (NPW )2 determine the size of system matrix to be solved.

Fig. 3. Unit cell of the electromagnetic Kronig-Penney problem.
Fig. 4. Band structures of the electromagnetic Kronig-Penney problem. MLS results are denoted by open circles and analytic results by solid lines.
Fig. 5. Convergence rates of the lowest five eigenvalues for the electromagnetic Kronig-Penney problem at k=(π/a)(0.5, 0.5).
Fig. 6. A square lattice of circular rod and its unit cell discretized by 1697 nodes.
Fig. 7. Results on band structures of square lattice composed of circular rods.
Fig. 8. Comparisons of MLS method and plane wave method (TM modes) for the square lattice of circular rods.

6. Concluding remarks

In this paper, the value-periodicMLS shape function is implemented in order to calculate the frequency band-structure of 2D photonic crystals. Through numerical examples, it has been verified that MLS results are more accurate than those of the conventional plane wave methods. Specifically TM modes results of meshfree method are outstanding, compared with TE mode results. This property has also been found in other real-space techniques such as finite element method and finite difference method [14

14. L. Shen, S. He, and S. Xiao, “A finite-diference eigenvalue algorithm for calculating the band structure of a photonic crystal,” Comput. Phys. Comm. 143, 213–221 (2002). [CrossRef]

, 16

16. D.C. Dobson, “An efficient band structure calculations in 2D photonic crystals,” J. Comput. Phys. 149, 363–376 (1999). [CrossRef]

]. In order to enhance the accuracy and convergence of TE modes solutions by the method, special techniques may be applied in our future study. Among them is the meshfree treatment of material discontinuity using Lagrange multipliers, which has been developed for solid mechanics by [21

21. L.W. Cordes and B. Moran, “Treatment of material discontinuity in the Element-Free Galerkin method,” Comput. Meth. Appl. Mech. Eng. 139, 75–89 (1996). [CrossRef]

]. Extending MLS band-structure computation to 3D photonic crystal is an example to which our periodic MLS shape function can immediately be applied. The effect of mechanical deformation on the change of photonic band is also an interesting subject.

Appendix

Here we describe a typical procedure how to define the shape function of the moving least-square (MLS) method. Although there are other routes to defining it, we mainly follow a recent Ref. [20

20. D.W. Kim and Y. Kim, “Point collocation method using the fast moving least-square reproducing kernel approximation,” Int. J. Numer. Methods Engrg. 56, 1445–1464 (2003). [CrossRef]

]. The MLS method begins with the local approximation of an arbitrary function u(x) as

uh(x,x¯)=p(xx¯)·a(x¯)=pT(xx¯)a(x¯)
(16)

in a small region around x̄. The column vector p(x) must contain at least 1 and monomials in order to fulfill the linear reproducing conditions. For example, p(x)=[1, x, y, z]T for 3D. By reproducing conditions, we mean that the MLS approximation can reproduce exactly the functions which appears in p(x). This is the well-known feature of MLS-based meshfree methods. Any polynomials of higher order or arbitrary types of functions can be added in the component of p(x) for the specific purpose of enrichment in MLS approximation.

After discretizing a problem domain by a set of points (I=1, …, NP), the idea of moving least-square interpolant is used in determining the coefficient vector a(x̄). That is, a(¯x) is determined by minimizing the local error residual functional that is expressed as

J(a(x¯))=I=1NPu(xI)p(xIx¯)·a(x¯)2W(xIx¯)
(17)

where the window function W(x) is a compactly-supported continuous function. By solving ∂J/∂a=0 for a(x̄), we have

a(x¯)=M1(x¯)I=1NPp(xIx¯)W(xIx¯)u(xI)
(18)

where the moment matrix M(x̄) is defined by

M(x¯)I=1NPp(xIx¯)pT(xIx¯)W(xIx¯)
(19)

uh(x,x¯)=p(xx¯)·a(x¯)=pT(xx¯)M1(x¯)I=1NPp(xIx¯)W(xIx¯)u(xI)
(20)

Finally, the global MLS approximation of u(x) is obtained by taking the limit of Eq. (20) as x̄→x,

uh(x)=I=1NPpT(0)M1(x)p(xIx)W(xIx)u(xI)
(21)

where p(0)=[1, 0, …, 0]T. The moment matrix M(x) becomes

M(x)=I=1NPp(xIx)pT(xIx)W(xIx)
(22)

Equation (21) can now be rewritten by given sampling values u I=u(xI) and the MLS shape function N I (x)=N(xI-x) as

uh(x)=I=1NPNI(x)uI
(23)

where

NI(x)=pT(0)M1(x)p(xIx)W(xIx)
(24)

Due to the symmetry of M(x), it is straightforward to verify that the shape function can also be written as

NI(x)=pT(xIx)M1(x)p(0)W(xIx)=pT(xIx)b(x)W(xIx)
(25)

Acknowledgements

This research was supported by a grant (02-K14-01-012-1-0) from Center for Nanoscale Mechatronics & Manufacturing of 21st Century Frontier Research Program of Korea.

References and links

1.

J.D. Joannopoulos, R.D. Meade, and J.N. Winn, Photonic Crystals: Molding the Flow of Light (Princeton University Press, Princeton, 1995).

2.

Y. Xia, “Photonic crystals,” Adv. Mater. 13, 369 (2001) and papers in this special issue. [CrossRef]

3.

K. Busch, “Photonic band structure theory: assesment and perspectives,” C. R. Physique 3, 53–66 (2002). [CrossRef]

4.

D. Cassagne, “Photonic band gap materials,” Ann. Phys. Fr. 23(4), 1–91 (1998). [CrossRef]

5.

J.B. Pendry, “Calculating photonic band structure,” J. Phys.: Condens. Matter 8, 1085–1108 (1996). [CrossRef]

6.

H.S. Sözüer, J.W. Haus, and R. Inguva, “Photonic bands: Convergence problems with the plane-wave method,” Phys. Rev. B 45, 13962–13972 (1992). [CrossRef]

7.

R.D. Meade, A.M. Rappe, K.D. Brommer, J.D. Joannopoulos, and O.L. Alerhand, “Accurate theoretical analysis of photonic band-gap materials,” Phys. Rev. B 48, 8434–8437 (1993). [CrossRef]

8.

C.T. Chan, Q.L. Yu, and K.M. Ho, “Order-N spectral method for electromagnetic waves,” Phys. Rev. B 51, 16635–16642 (1995). [CrossRef]

9.

A.J. Ward and J.B. Pendry, “Calculating photonic Green’s functions using a nonorthogonal finite-difference time-domain method,” Phys. Rev. B 58, 7252–7259 (1998). [CrossRef]

10.

M. Qiu and S. He, “A nonorthogonal finite-difference time-domain method for computing the band structure of a two-dimensional photonic crystal with dielectric and metallic inclusions,” J. Appl. Phys. 87, 8268–8275 (2000). [CrossRef]

11.

K.M. Leung and Y. Qiu, “Multiple-scattering calculation of the two-dimensional photonic band structure,” Phys. Rev. B 48, 7767–7771 (1993). [CrossRef]

12.

X. Wang, X.G. Zhang, Q. Yu, and B.N. Harmon, “Multiple-scattering theory for electromagnetic waves,” Phys. Rev. B 47, 4161–4167 (1993). [CrossRef]

13.

J.B. Pendry and A. MacKinnon, “Calculation of photon dispersion relations,” Phys. Rev. Lett. 69, 2772–2775 (1992). [CrossRef] [PubMed]

14.

L. Shen, S. He, and S. Xiao, “A finite-diference eigenvalue algorithm for calculating the band structure of a photonic crystal,” Comput. Phys. Comm. 143, 213–221 (2002). [CrossRef]

15.

W. Axmann and P. Kuchment, “An efficient finite element method for computing spectra of photonic and acoustic band-gap materials: I. Scalar case,” J. Comput. Phys. 150, 468–481 (1999). [CrossRef]

16.

D.C. Dobson, “An efficient band structure calculations in 2D photonic crystals,” J. Comput. Phys. 149, 363–376 (1999). [CrossRef]

17.

C. Mias, J.P. Webb, and R.L. Ferrari, “Finite element modelling of electromagnetic waves in doubly and triply periodic structures,” IEE Proc.-Optoelectron. , 146(2), 111–118 (1999). [CrossRef]

18.

M. Marrone, V.F. Rodriguez-Esquerre, and H.E. Hernandez-Figueroa, “Novel numerical method for the analysis of 2D photonic crystals: the cell method,” Opt. Express 10, 1299–1304 (2002), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-10-22-1299 [CrossRef] [PubMed]

19.

S. Li and W.K. Liu, “Meshfree and particle methods and their applications,” Applied Mechanics Review , 55, 1–34 (2002). [CrossRef]

20.

D.W. Kim and Y. Kim, “Point collocation method using the fast moving least-square reproducing kernel approximation,” Int. J. Numer. Methods Engrg. 56, 1445–1464 (2003). [CrossRef]

21.

L.W. Cordes and B. Moran, “Treatment of material discontinuity in the Element-Free Galerkin method,” Comput. Meth. Appl. Mech. Eng. 139, 75–89 (1996). [CrossRef]

OCIS Codes
(000.4430) General : Numerical approximation and analysis
(260.2110) Physical optics : Electromagnetic optics

ToC Category:
Research Papers

History
Original Manuscript: February 6, 2003
Revised Manuscript: March 12, 2003
Published: March 24, 2003

Citation
Sukky Jun, Young-Sam Cho, and Seyoung Im, "Moving least-square method for the band-structure calculation of 2D photonic crystals," Opt. Express 11, 541-551 (2003)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-11-6-541


Sort:  Journal  |  Reset  

References

  1. J.D. Joannopoulos, R.D. Meade, and J.N. Winn, Photonic Crystals: Molding the Flow of Light (Princeton University Press, Princeton, 1995).
  2. Y. Xia, �??Photonic crystals,�?? Adv. Mater. 13, 369 (2001) and papers in this special issue. [CrossRef]
  3. K. Busch, "Photonic band structure theory: assesment and perspectives,�?? C. R. Physique 3, 53-66 (2002). [CrossRef]
  4. D. Cassagne, "Photonic band gap materials,�?? Ann. Phys. Fr. 23(4), 1-91 (1998). [CrossRef]
  5. J.B. Pendry, �??Calculating photonic band structure,�?? J. Phys.: Condens. Matter 8, 1085-1108 (1996). [CrossRef]
  6. H.S. Sozuer, J.W. Haus, and R. Inguva, �??Photonic bands: Convergence problems with the planewave method,�?? Phys. Rev. B 45, 13962-13972 (1992). [CrossRef]
  7. R.D. Meade, A.M. Rappe, K.D. Brommer, J.D. Joannopoulos, O.L. Alerhand, �??Accurate theoretical analysis of photonic band-gap materials,�?? Phys. Rev. B 48, 8434-8437 (1993). [CrossRef]
  8. C.T. Chan, Q.L. Yu, and K.M. Ho, �??Order-N spectral method for electromagnetic waves,�?? Phys. Rev. B 51, 16635-16642 (1995). [CrossRef]
  9. A.J. Ward and J.B. Pendry, �??Calculating photonic Green�??s functions using a nonorthogonal finite difference time-domain method,�?? Phys. Rev. B 58, 7252-7259 (1998). [CrossRef]
  10. M. Qiu and S. He, �??A nonorthogonal finite-difference time-domain method for computing the band structure of a two-dimensional photonic crystal with dielectric and metallic inclusions,�?? J. Appl. Phys. 87, 8268-8275 (2000). [CrossRef]
  11. K.M. Leung and Y. Qiu, �??Multiple-scattering calculation of the two-dimensional photonic band structure,�?? Phys. Rev. B 48, 7767-7771 (1993). [CrossRef]
  12. X. Wang, X.G. Zhang, Q. Yu, and B.N. Harmon, �??Multiple-scattering theory for electromagnetic waves,�?? Phys. Rev. B 47, 4161-4167 (1993). [CrossRef]
  13. J.B. Pendry and A. MacKinnon, �??Calculation of photon dispersion relations,�?? Phys. Rev. Lett. 69, 2772-2775 (1992). [CrossRef] [PubMed]
  14. L. Shen, S. He, and S. Xiao, �??A finite-difference eigenvalue algorithm for calculating the band structure of a photonic crystal,�?? Comput. Phys. Comm. 143, 213-221 (2002). [CrossRef]
  15. W. Axmann and P. Kuchment, �??An efficient finite element method for computing spectra of photonic and acoustic band-gap materials: I. Scalar case,�?? J. Comput. Phys. 150, 468-481 (1999). [CrossRef]
  16. D.C. Dobson, �??An efficient band structure calculations in 2D photonic crystals,�?? J. Comput. Phys. 149, 363-376 (1999). [CrossRef]
  17. C. Mias, J.P. Webb and R.L. Ferrari, �??Finite element modelling of electromagnetic waves in doubly and triply periodic structures,�?? IEE Proc. Optoelectron. 146(2), 111-118 (1999). [CrossRef]
  18. M. Marrone, V.F. Rodriguez-Esquerre, and H.E. Hernandez-Figueroa, �??Novel numerical method for the analysis of 2D photonic crystals: the cell method,�?? Opt. Express 10, 1299-1304 (2002), <a href="http://www.opticsexpress.org/abstract.cfm?URI=OPEX-10-22-1299">http://www.opticsexpress.org/abstract.cfm?URI=OPEX-10-22-1299</a> [CrossRef] [PubMed]
  19. S. Li and W.K. Liu, �??Meshfree and particle methods and their applications,�?? Appl. Mechanics Rev. 55, 1-34 (2002). [CrossRef]
  20. D.W. Kim and Y. Kim, �??Point collocation method using the fast moving least-square reproducing kernel approximation,�?? Int. J. Numer. Methods Engrg. 56, 1445 - 1464 (2003). [CrossRef]
  21. L.W. Cordes and B. Moran, �??Treatment of material discontinuity in the Element-Free Galerkin method,�?? Comput. Meth. Appl. Mech. Eng. 139, 75-89 (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