In this paper, the timedomain representation for a plane wave is derived based on a inverse fourier transform b laplace domain representation and c plane wave expansion, all of which are mathematically equivalent in the continuous time domain. From those we will derive a set of equations, whose solutions can be approximated by the plane wave expansion method. From the geometrical point of view, the wave diffraction can be described as fourier decomposition of light beam into plane waves which, in propagationacquire, phase shifts depending on their propagation direction. Numerical investigations of twodimensional photonic. Osa planewave expansion method for calculating band.
Let us examine this set of examples in further detail. Plane wave expansion method pwe refers to a computational technique in electromagnetics to solve the maxwells equations by formulating an eigenvalue problem out of the equation. Plane wave expansion pwe method the maxwell equation in a transparent, timeinvariant, source free, and nonmagnetic medium can be written in the following form. In actual calculations, we have to limit the plane wave expansion at some cutoffs are always given in energy units such as rydberg or ev corresponding to the kinetic. The plane wave expansion of a spherical wave is surprisingly hard to find on the web. Wave functions waveguides and cavities scattering separation of variables the special functions vector potentials the spherical bessel equation each function has the same properties as the corresponding cylindrical function. Fourier series and fourier transforms the fourier transform is one of the most important tools for analyzing functions. Introduction to planewave basis sets and pseudopotential. Plane wave expansion method used to engineer photonic. Note that emtl was initially developed mostly for fdtd simulations. Planewave expansion method for calculating band structure of photonic crystal slabs with perfectly matched layers. The plane wave expansion method for the mpcs is briefly considered in sec. The plane wave expansion method uses the set of planes waves expik r as the complete basis.
It expands the field into a set of plane waves and converts maxwells equations to matrix formbyassigning each plane wave a complex amplitude. This algorithm combines the plane wave expansion method with perfectly matched layers for the termination of the computational region in the direction out of the plane. This method is popular among the photonic crystal community as a method of solving for the band structure dispersion relation of specific photonic crystal geometries. This solution is a wave \traveling in the direction of k in the sense that a point of constant phase, meaning k. In particular, a plane to plane method for the determination of the aperture. Such method has been used for treating the twodimensional scattering by a set of perfectly conducting parallel objects. Spherical wave expansion of vector plane wave in discussing the scattering or absorption of electromagnetic radiation by localized systems, it is useful to be able to express a plane electromagnetic wave as a superposition of spherical waves. In, the expansion into plane waves of a cylindrical wave propagating in a lossless medium has been derived. The plane wave expansion method pwem provides a highly efficient numerical solution to maxwells equations for devices with low to moderate dielectric contrast.
Outofplane dispersion and homogenization in photonic. That is, when the average interaction energy between the incident particle and the scattering potential is much smaller than the particles incident kinetic energy, the scattered wave can be considered to be a plane wave. This file discussess the following topics related to the main manuscript. By doing so we will derive the dispersion relation and the density.
This manual is based on my implementation of pwm plane wave method using matlab. The photonic energy band structure of periodic composite materials, consisting of regular spherical arrays or scaffoldings of one kind of dielectric material in a different background material, are studied using the plane wave expansion method in order. The extended plane wave expansion method for 3d anisotropic. The guidedmode expansion method for photonic crystal slabs. Pw expansion of nic has to be converged with density cuto. The method described using the expansion creates spurious modes with zero frequency, which were removed. Planewave expansion method for calculating band structure. One reason for choosing the plane wave basis is that it is the. Can we decompose a spherical wave as a sum of plane waves. The plane wave expansion method pwem provides a highly efficient numerical solution to maxwell us s equations for devices with low to moderate dielectric contrast. The feed lies outside the closed surface, so that the volume v. In particular, this allows for the possibility that the projected characteristics may cross each other.
Rami mehrem 2009, the plane wave expansion, infinite integrals and identities involving spherical bessel functions, arxiv. Planewave expansion of cylindrical functions in lossy media. Oct 04, 20 this lecture steps the student through the formulation and implementation of the plane wave expansion method. Planewave expansion method for calculating band structure of. The plane wave expansion method was carried out to construct dispersion curves for along the brillouin zone shown in figure 4. It expands the field into a set of plane waves and converts maxwells equations to matrix form by assigning each plane wave a complex amplitude. Plane wave expansion of cylindrical functions in lossy media. Plane wave expansion method for photonic band gap calculation using matlab manual version 1. The theoretical investigation is performed using the plane wave expansion method in both the optical and acoustical regimes.
The expansion can be used to solve many wave propagation problems in complicated structures, including laterally inhomogeneous media with curved interfaces. Matlab implementation of the rigorous coupled wave analysis and plane wave expansion method chrisfaddenrcwa. Plane wave expansion dielectric slab waveguide waveguideembedded photonic crystals 2. For the most part, however, it is intended to be applied to structures with some degree of periodicity in the refractive index distribution. Digital library of mathematical functions, equation 10. Plane wave approximation since many articles are devoted to the mathematical aspects of the fmm, we will skip the mathematical derivation of the algorithm. Arpes band structure and photoemission intensity simulations using the electron plane wave expansion epwe method. Research article an enhanced plane wave expansion method. The plane wave pseudopotential method 1 bornoppenheimer approximation 2 apply densityfunctional theory dft to calculate the electronic structure. We show through a numerical convergence study and rigorous analysis on a simplified problem that the convergence of the planewave expansion method is. The convergence criterion for the electronic structure iteration was set. Numerical computation of band gaps in photonic crystal fibres. The projector augmented wave paw method was used to describe the ionic cores.
Timedomain representations of a plane wave with spatial. This method is popular among the photonic crystal community as a method of solving for the band structure dispersion relation of specific photonic crystal. The calculated eigenfrequencies and field patterns for a defect in a square array of dielectric rods are consistent with those obtained by the plane wave expansion method. The plane wave expansion method allows the calculation of dispersion curves i.
Growth details and zigzag chain alignment on a ag111 curved vicinal crystal. The extended plane wave expansion method for 3d anisotropic photonic crystals1 y. This thesis presents the results of numerical studies of scalar and electromagnetic wave propagation in periodic and random dielectric media. Fourier factorization in the plane wave expansion method in modeling photonic crystals 3 and 7 will propose how to factorize anisotropic and threedimensional 3d photonic crystals, respectively. Simultaneous photonic and phononic band gap, defect. Br effect on the electronic structure of zigzag chains studied by. Research article an enhanced plane wave expansion method to solve piezoelectric phononic crystal with resonant shunting circuits ziyanglian, 1,2 shanjiang, 1,2 hongpinghu, 1,2 longxiangdai, 1,2 xuedongchen, 3 andweijiang 3 department of mechanics, huazhong university of science and technology, wuhan, china. Photonic band structure and eigenmodes of magnetophotonic. We implemented plane wave expansion pwe method as a part of electromagnetic template library emtl. The plane wave expansion method is one of the numerical methods that has been used. Planewave expansion methods for photonic crystal fibres.
There is a further truncation of plane wave expansion in plane wave calculations. Pdf band gap calculation using the plane wave expansion. In this document we will apply the plane wave expansion method on a speci c photonic crystal. Lecture 18 cem plane wave expansion method youtube. Starting from maxwells equations the propagation of electromagnetic waves can be described by maxwells equations. A general complex eigenvalue problem is then obtained. Pdf fourier factorization in the plane wave expansion. We present a new algorithm for calculation of the band structure of photonic crystal slabs. The basic underlying idea is that a function fx can be expressed as a linear combination of elementary functions speci cally, sinusoidal waves. Pdf we present a new algorithm for calculation of the band structure of photonic crystal slabs.
That means that one must be able to expand one in terms of the other. For that we refer the reader to previous articles 21, for example, or 55, which give the details of the method. Plane wave approximation of homogeneous helmholtz solutions a. The first is energy cutoffs, which is the cutoff for the wave the input files page 4 of 20. Also, an alternative method to obtain the expansion coefficients of the plane.
The plane wave expansion pwe method have been formulated in a general form for arbitrary 2d photonic crystal and applied for detailed analysis of the. In the plane wave appproximation, a particle incident on a nucleus has a wavefunction approximated by a plane wave of the form eik. Plane wave expansion method for photonic band gap calculation. For example, a 4 x 4 supercell with a central defect can give reasonable accuracy because the missing holes are spaced four lattice units apart. Band structure of honeycomb photonic crystal slabs taii weng and g. This method is popular among the photonic crystal community as a method of. Summary the plane wave expansion method allows the calculation of dispersion curves i. The band structures are calculated by the plane wave expansion method and the. Pdf planewave expansion method for calculating band. Plane wave expansion method used to engineer photonic crystal sensors with high efficiency.
It describes how to construct electromagnetic band diagrams and isofrequency contours. Computational methods for the analysis and design of photonic. Plane waves can be expanded in terms of free spherical waves by. The larger the rc i, the heavier the terms in direct space. Expansion of a plane wave into gaussian beams springerlink.
The dephasing of the plane wave components results in the diffractive broadening of the beam, as showed in fig. Photonic crystals are composed of periodic, dielectric structures that affect the propaga tion of electromagnetic waves. Plane wave expansion method in computational electromagnetism. Consider, first of all, the expansion of a scalar plane wave as a set of scalar spherical waves. Plane wave expansion method for calculating band structure of photonic crystal slabs with perfectly matched layers shouyuan shi, caihua chen, and dennis w. There are some tricks in the derivation process, so we show the process in detail. Fourier factorization in the plane wave expansion method in modeling photonic crystals 18 willbesetbyintech the numerical ef. Application of the plane wave expansion method to a two.
We combine cartesian coordinates and polar coordinates wave number eigenvalue equations based on the plane wave expansion pwe method to calculate and optimize the band structures of the twodimensional 2d metal photonic crystals phcs. Guidedmode expansion method photonic mode dispersion intrinsic diffraction losses convergence tests, comparison with other methods examples of applications books yariv a. Introduction to two dimensional scatering 3 let the total wave be the sum of the incident and scattered waves. Plane waves and free spherical waves both form an onshell complete orthnormal set on with or without the origin. The plane wave spectrum is used to determine the far. Course paperwork syllabus homework course topics resources. Plane wave expansion method 1d photonic crystal considering electromagnetic wave in materials with.
Fourier factorization in the plane wave expansion method. Plane wave approximation of homogeneous helmholtz solutions. Acknowledgements this work would not have been successful without the guidance and the help of several individuals. For the planewave expansion, a 450 ev kinetic energy cutoff was used. Yes, a spherical wave can be expanded as a sum of plane waves, and vice versa. Namely, only the reciprocal lattice vectors whose kinetic energy lower than a predefined maximum cutoff energy, are kept in the expansion, while the rest of the coefficients are set to zero. An integral expansion which expresses a plane monochromatic wave as a superposition of gaussian beams is found. Reformulation of the plane wave method to model photonic crystals. Fullpotential linearized augmented planewave method. In addition, the effectivemedium tensor is applied to improve convergence. Furthermore, we have a plane wave, by which we mean that a surface of constant phase is a plane.
In this paper, we study the approximation of solutions of the homogeneous helmholtz equation. Improved planewave expansion method for band structure. The term is also used, even more specifically, to mean a monochromatic or sinusoidal plane wave. There are several methods to calculate the dispersion relation and therewith the range of the band gaps. The spatial frequency spectrum of the cws can obviously be obtained through numerical fourier transform techniques, however, the knowledge of the analytical closedform expression of the spectral representation gives substantial benefits. To study defect modes, the same plane wave expansion method can be used. Plane wave 5 represented in the first illustration toward the right is a linearly polarized, electromagnetic wave. Gu, institute of physics, academia sinica china conventional plane wave expansion pwe methods, good for calculating such properties as photonic band gaps for. I show how i deal with the 1d, 2d and 3d cases, so you can deal with other problems in. Plane wave approximation of homogeneous helmholtz solutions 3. The matrix method has been implemented, and its potentialities are presented. The basic idea is to replace the unit cell by a more complicated unit cell and preserve the periodicity. The planewave pseudopotential method 1 bornoppenheimer approximation 2 apply densityfunctional theory dft to calculate the electronic structure. Compared with the traditional pwe methods for metal phcs, the band structures can.
In fact, an incoming plane wave along the direction can be expanded in bessel functions. Partial wave analysis of scattering we can take a quick look at scattering from a potential in 3d we assume that far from the origin so the incoming and outgoing waves can be written in terms of our solutions for a constant potential. I am not going to list all source files here, they are pretty short and can be easily modified for your own purpose if you know a little bit knowledge about matlab programming. The condition for solving fors and t in terms ofx and y requires that the jacobian matrix be nonsingular. Pdf planewave expansion method for calculating band structure. The maxwell equation in a transparent, timeinvariant, source free, and nonmagnetic medium can be written in the following form.
Classical wave propagation in periodic and random media by. Unfortunately, the function that describes the material in the. By means of it one can find as many defect modes as possible, including those that are inactive to the incident plane wave in the fdtd transmission spectra method. We apply a generalization of the plane wave expansion method and calculate the pbs and eigenmodes of mpcs in different magnetization geometries and analyze their polarization state. The plane wave expansion, or the rayleigh equation, is commonly used in nuclear physics. Second choice is a core of plane wave expansion pwe method that is very suitable for calculation of photonic crystals band structures. The timeindependent schrodinger equation for the energy eigenstates in the coordinate representation is given by. Because this is a plane wave, each blue vector, indicating the perpendicular displacement from a point on the axis out to the sine.
417 700 62 1305 1397 1148 1275 594 648 938 242 888 754 390 733 262 340 25 63 425 842 484 469 33 231 433 1500 111 1217 506 603 1055 638 504 440 1130 761 553