<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JEMAA</journal-id><journal-title-group><journal-title>Journal of Electromagnetic Analysis and Applications</journal-title></journal-title-group><issn pub-type="epub">1942-0730</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jemaa.2011.311074</article-id><article-id pub-id-type="publisher-id">JEMAA-8339</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Studies on Electromagnetic Waves Properties in Metallic-Dielectric Photonic Crystal
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>rafa</surname><given-names>H. Aly</given-names></name><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ehab</surname><given-names>Abdel-Rahman</given-names></name></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hassan</surname><given-names>S. Hanafey</given-names></name></contrib></contrib-group><author-notes><corresp id="cor1">* E-mail:<email>arafa16@yahoo.com(RHA)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>04</day><month>11</month><year>2011</year></pub-date><volume>03</volume><issue>11</issue><fpage>465</fpage><lpage>470</lpage><history><date date-type="received"><day>July</day>	<month>15th,</month>	<year>2011</year></date><date date-type="rev-recd"><day>August</day>	<month>12th,</month>	<year>2011</year>	</date><date date-type="accepted"><day>September</day>	<month>2nd,</month>	<year>2011.</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  We have performed numerical analyses to investigate the propagation characteristics of a simple-one dimensional metallic-dielectric structure. Several different metals have been studied (aluminum, gold, copper, and silver). Copper gives the smallest absorption and aluminum is more absorptive. We have compared between the transmittance of these structures with different metals and the same dielectric material. Also we used these different metals with the same layer thickness, incidence angels, number of layers and the same electromagnetic waves range.
 
</p></abstract><kwd-group><kwd>Transmittance</kwd><kwd> PC</kwd><kwd> Pasmonic</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Photonic crystals (PCs) with a periodic arrangement of dielectric or metallic elements in one, two or three dimensions affect the properties of photons in the same way as a periodic potential affects the properties of electrons in semiconductor crystals. (PCs) are structures with periodically modulated dielectric constants whose distribution follows a periodicity of the order of a fraction of the optical wavelength. Photonic crystals or photonic band gap materials (PBG) are a new class of optical materials which have the properties to manipulate the flow of electromagnetic waves within it. Researchers have proposed many new applications of photonic devices which may revolutionize the field of photonics in much the same way as semiconductors revolutionized electronics. They can generate spectral regions named photonic band gaps (PBGs) where light cannot propagate in a manner analogous to the formation of electronic band gaps in semiconductors [1,2]. There are several studies of metallic [3-7] and superconducting photonic crystals [8-11] which are mostly concentrated at microwave, millimeterwave, and far-infrared frequencies. In those frequencies, metals act like nearly perfect reflectors with no significant absorption problems.</p><p>In the recent years there are many authors using a composite structure which is composed of (metal/dielectric) structure instead of (dielectric/dielectric) structure, since the metallic-dielectric photonic crystals (MDPCs) can be useful for a large number of applications and also will open the door to new optical devices that depend mainly on the optical properties of the metal.</p><p>There are certain advantages to introducing 1D MDPCs. These include reduced size and weight, easier fabrication methods, and lower costs. Here, we use the transfer matrix method to study the effect of absorption of metal at the near-infrared and visible frequencies. In particular, we study a simple-one dimensional metallic-dielectric structure. Silver, Gold, Copper, and Aluminum have been used in order to study the effect of different metals on optical properties.</p></sec><sec id="s2"><title>2. Theoretical Analysis</title><p>In order to study the electromagnetic waves propagating through photonic crystal structures and determine the optical properties of these periodic structures. We will explain in brief a mathematical treatment with a simple one dimensional photonic crystal structure (1DPC) (see <xref ref-type="fig" rid="fig1">Figure 1</xref>) which is composed of two materials with thicknesses (<img src="6-9801214\4ca59aaa-5bc2-4230-b8d7-ed893b5178a4.jpg" />and<img src="6-9801214\6f7b40ae-9013-4b6a-ab14-50199022f9d5.jpg" />) and refractive indices (<img src="6-9801214\5dabe4fd-fee3-483d-9e77-cba30aa014fd.jpg" />and<img src="6-9801214\1969bc4a-45e5-4a68-9d92-9812a876686d.jpg" />) respectively. The analysis of the incident electromagnetic radiation on this structure will be performed using the transfer matrix method (TMM).</p><p>We will produce TMM in the form of the dynamical and the propagating matrices which had been used to describe the wave interaction through each layer and the wave response at the interfaces between these layers [<xref ref-type="bibr" rid="scirp.8339-ref12">12</xref>]. Our analysis will be start considering that, the electro-</p><p>magnetic waves are propagating through the multilayer structure along x-direction, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. The used periodic structure can be defined along this direction as:</p><disp-formula id="scirp.8339-formula129376"><label>(1)</label><graphic position="anchor" xlink:href="6-9801214\60f45ee6-132e-43b5-bddf-fde0640ae067.jpg"  xlink:type="simple"/></disp-formula><p>The electric field of a general plane-wave solution can be written as<img src="6-9801214\4ee70e16-28ef-4924-8c1d-0c7b91074603.jpg" />, where <img src="6-9801214\12152aed-8872-4016-b326-739c723bb2e2.jpg" /> is the z-component of the wave vector; <img src="6-9801214\d2f14cbf-84c1-4a23-ba1b-e701dcd76c7c.jpg" /> <img src="6-9801214\a37762b9-3fb9-46ff-9bdf-0932bdaecace.jpg" />, <img src="6-9801214\9d108bc6-474a-4f27-9b5e-1facfb6a292e.jpg" />is the angle of the incident ray in each layer and c is the speed of light then the electric field distribution <img src="6-9801214\564a0fbb-c9f2-4625-8534-296c96de7a79.jpg" />can be written as:</p><disp-formula id="scirp.8339-formula129377"><label>(2)</label><graphic position="anchor" xlink:href="6-9801214\22aa4d2b-aebb-4f49-879c-c5ceadea49d4.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="6-9801214\85e8b0d9-10e6-4b63-b0aa-9cac9cbcc394.jpg" />is the x-component of the wave vector, <img src="6-9801214\ad417620-df85-4bc8-8993-59b23639246f.jpg" />, <img src="6-9801214\7889387b-d223-49dd-8554-671f6af90525.jpg" />and <img src="6-9801214\82da9364-9a84-43f2-a0f7-dbba089bdd34.jpg" />represent the amplitudes of the plane wave at each interface<img src="6-9801214\d883f0b2-3e6c-4d96-9fc2-77c1d8fbff40.jpg" />. It’s wellknown that this distribution of the electric field was taken for<img src="6-9801214\32ff1283-f3f7-40c7-85d3-bd7c10ba51de.jpg" />, and we will now obtain the general case for N period. The amplitudes of the plane waves at different layers can be related by</p><disp-formula id="scirp.8339-formula129378"><label>(3)</label><graphic position="anchor" xlink:href="6-9801214\31d78028-1fc3-4da5-8beb-0eb3a253591e.jpg"  xlink:type="simple"/></disp-formula><p>with<img src="6-9801214\58cb9192-e2cf-4651-a407-2bb18db69d9b.jpg" />.</p><p>where matrices D (dynamical matrix) and P(propagation matrix) can be written as:</p><disp-formula id="scirp.8339-formula129379"><label>(4a)</label><graphic position="anchor" xlink:href="6-9801214\4f08444b-3255-4df6-bfbc-65df8e7d8f28.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.8339-formula129380"><label>(4b)</label><graphic position="anchor" xlink:href="6-9801214\50ff6c58-1b7f-4179-b674-c60398e74e43.jpg"  xlink:type="simple"/></disp-formula><p>Since the propagation matrix can be written in the form of sine and cosine functions instead of the exponential function for simplifying</p><p><img src="6-9801214\b8f83521-00d7-4f35-87e5-5ac005b89944.jpg" />(5)</p><p>where</p><p><img src="6-9801214\b333576f-3855-4184-ac91-35f31eadc6f0.jpg" />&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160;&#160; <img src="6-9801214\e6fb932d-ac69-470f-97e5-1de3daaa18a4.jpg" /></p><p>The relation between the amplitudes of the plane wave can be written as</p><disp-formula id="scirp.8339-formula129381"><label>(6)</label><graphic position="anchor" xlink:href="6-9801214\aed899c7-9a69-4e87-92c5-19cf02cc1304.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="6-9801214\bd90979c-e778-473b-8a16-a2be996276c2.jpg" /> is the matrix for one period without the effect of the vacuum and substrate, whose elements m<sub>11</sub>, m<sub>12</sub>, m<sub>21</sub> and m<sub>22</sub> was computed for TE-waves, which can be obtained also for TM-waves using the same analysis, and<img src="6-9801214\e410162a-0b94-4881-937d-3d4d5ede3f0b.jpg" />, is the lattice constant. This matrix can take the following form [<xref ref-type="bibr" rid="scirp.8339-ref12">12</xref>].</p><disp-formula id="scirp.8339-formula129382"><label>(7)</label><graphic position="anchor" xlink:href="6-9801214\012e240b-031a-4fa8-b2a5-e2df361b0c56.jpg"  xlink:type="simple"/></disp-formula><p>with</p><disp-formula id="scirp.8339-formula129383"><label>(8)</label><graphic position="anchor" xlink:href="6-9801214\138270f9-2479-4f70-8665-fdd610785233.jpg"  xlink:type="simple"/></disp-formula><p>We can write the matrix of N periods in the form [<xref ref-type="bibr" rid="scirp.8339-ref12">12</xref>]:</p><disp-formula id="scirp.8339-formula129384"><label>(9a)</label><graphic position="anchor" xlink:href="6-9801214\b9427968-f64e-4196-842d-d9e8bd6a1a7b.jpg"  xlink:type="simple"/></disp-formula><p>Since,</p><disp-formula id="scirp.8339-formula129385"><label>(9b)</label><graphic position="anchor" xlink:href="6-9801214\49ba68d7-4be6-42fd-b8e9-9bf33a120648.jpg"  xlink:type="simple"/></disp-formula><p>where M(Na) is the matrix for N period, and Q<sub>11</sub>, Q<sub>12</sub>, Q<sub>21</sub> and Q<sub>22</sub> are parts of their elements that can be related to the elements of the single period matrix(m<sub>11</sub>, m<sub>12</sub>, m<sub>21</sub>, and m<sub>22</sub>) by the relations[<xref ref-type="bibr" rid="scirp.8339-ref13">13</xref>]:</p><disp-formula id="scirp.8339-formula129386"><label>(10)</label><graphic position="anchor" xlink:href="6-9801214\1a88c99c-22f8-4c48-a85c-21a26023ae82.jpg"  xlink:type="simple"/></disp-formula><p>with</p><disp-formula id="scirp.8339-formula129387"><label>(11)</label><graphic position="anchor" xlink:href="6-9801214\b0a761f4-65f0-44c5-8138-e3bb41fd48bf.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="6-9801214\8e4cc123-5c40-411e-8dc1-f2395d618a46.jpg" />, is Chebyshev polynomials of the second kind. Then by using the above expressions we can obtain the reflection and the transmission coefficients as:</p><disp-formula id="scirp.8339-formula129388"><label>(12)</label><graphic position="anchor" xlink:href="6-9801214\39d4d040-1cd1-4319-88fb-2d1ec080833a.jpg"  xlink:type="simple"/></disp-formula><p>Finally we can calculate the transmittance and the reflectance using the following expressions:</p><disp-formula id="scirp.8339-formula129389"><label>(13)</label><graphic position="anchor" xlink:href="6-9801214\ef64d866-5981-4888-bbc4-191799efe2c8.jpg"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.8339-formula129390"><label>(14)</label><graphic position="anchor" xlink:href="6-9801214\6e27f985-f5e0-411d-b140-94ff7a5598c3.jpg"  xlink:type="simple"/></disp-formula><p>This gives us the final form of the transmission and reflection coefficients for the incident electromagnetic waves inside this 1D-binary periodic structure as a function of the dynamical and the propagating matrices depending on the numerical solution of Maxwell’s using the TMM.</p></sec><sec id="s3"><title>3. Results and Discussions</title><p>The periodicity of the permittivity plays the same role for the photons that propagate inside the structure than the atomic potential for the electrons. Leading further this analogy, the thicknesses and the index contrast of the photonic crystal determinate many of its optical properties as it does for conduction properties of semiconductors. We have used for our simulation Essential Macleod software. The Essential Macleod is a comprehensive software package for the design, analysis, manufacture and trouble shooting of thin film optical coatings. The program has a true multiple document interfaces. It can handle coatings from rugates to ultrafast, from singlelayer antireflection coatings to demanding color separation beamsplitters, providing all that is necessary for their design, analysis, production planning and even reverse engineering of a failed production attempt. A wide range of performance parameters, from straightforward transmittance and reflectance to color coordinates in different color spaces, is built into the software. It can synthesize designs from scratch or refine existing ones, investigate errors and extract optical constants of film materials for use in designs.</p><p>We can obtain the range of frequencies for which the kind of light cannot propagate inside the materials and the others can propagate. By reducing the size of the elementary cell of the periodic lattice, the frequencies can shift to higher values. The consequence of this property is the possibility to transpose a photonic crystal design from the microwave domain to infrared or visible range. In our results we have studied one dimensional metallic (different metals) dielectric (SiO<sub>2</sub>) photonic crystals (MDPC’s). For our calculations we have used the thickness of the dielectric (SiO<sub>2</sub>) layer is 344 nm, the thickness layer of metal is 5nm, the incidence of angle is 35˚ and the number of periods 10 for all our results.</p><p>In the case of Ag-SiO<sub>2</sub> (see <xref ref-type="fig" rid="fig2">Figure 2</xref>) we can see the magnitude of transmittance near the unity with visible light and there is a small photonic band gap (PBG) appeared within NIR (750 - 930 nm). The magnitude of transmittance in the case of Au-SiO<sub>2</sub> (<xref ref-type="fig" rid="fig3">Figure 3</xref>) is less than 10% around 400 nm and increased to ~70 % within 550 - 700 nm and then we obtained the second and large PBG in the range 775 - 925 nm. In <xref ref-type="fig" rid="fig4">Figure 4</xref>, the magnitude of transmittance in the case of Cu-SiO<sub>2</sub> is go downing around 60% within the range of visible light (400 - 700) and we obtained the same PBG between 775 - 925 nm.</p><p>In <xref ref-type="fig" rid="fig5">Figure 5</xref>, we have examined Al-SiO<sub>2</sub> with the same conditions and we found Aluminum that gives the lowest possible transmittance within the visible and near IR frequencies with the complete PBG between 600 – 900 nm. The significance of the magnitude of incident angles also is obtained in <xref ref-type="fig" rid="fig6">Figure 6</xref> with θ = 25˚, 30˚, 35˚ and 40˚. The behavior of incident angle has different influences on the forbidden bands of TE and TM modes. This work will be a useful reference when designing and producing some apparatus such as, optical filter, eyes-protecting to laser and heat-absorbing window.</p></sec><sec id="s4"><title>4. Conclusions</title><p>We have studied metallic photonic crystals at visible and</p><p>near-infrared wavelengths with the transfer matrix method. We focused on the transmittance of these structures, thus we studied the simple 1D structure. We expect that our conclusions regarding the transmittance will hold to any other metallic structures. By comparing the results for different metals, we found that copper gives the lowest possible transmittance within the visible frequencies. Silver gives slightly higher transmittance within the visible with small PBG near IR. Gold gives low possible transmittance within short visible and going to near 70% between 550 - 700 nm with small PBG within 750 - 850 nm. Aluminum is very lossy and is not recommended for optical photonic crystals. Silver and gold are acceptable although slightly lower in performance. We also found our results are very sensitive to the thickness of layers as well as PBG is sensitive to the kind of polarizations. This work will be a useful reference when designing and producing some apparatus such as, eyes-protecting to laser and heat-absorbing window.</p></sec><sec id="s5"><title>5. Acknowledgements</title><p>The corresponding author Arafa H. Aly thankful Prof Macleod for his offers me the Essential Macleod program to achieve this work.</p></sec><sec id="s6"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.8339-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">E. Yablonovitch, “Inhibited Spontaneous Emission in Solid-State Physics and Electronics,” Physical Review Letters, Vol. 58, No. 20, 1987, pp. 2059-2062.  
doi:10.1103/PhysRevLett.58.2059</mixed-citation></ref><ref id="scirp.8339-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">S. John, “Strong Localization of Photons in Certain Disordered Dielectric Superlattices,” Physical Review Letters, Vol. 58, No. 23, 1987, pp. 2486-2489.  
doi:10.1103/PhysRevLett.58.2486</mixed-citation></ref><ref id="scirp.8339-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">J. A. Sánchez-Gil and A. A. Maradudin, J. Q. Lu and V. D. Freilikher, “Transmission of Electromagnetic Waves through Thin Metal Films with Randomly Rough Surfaces,” Physical Review B, Vol. 51, No. 23, 1995, pp. 17100-17115. doi:10.1103/PhysRevB.51.17100</mixed-citation></ref><ref id="scirp.8339-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">M. M. Sigalas, C. T. Chan, K. M. Ho and C. M. Soukoulis, “Metallic Photonic Band-Gap Materials,” Physical Review B, Vol. 52, No. 16, 1995, pp. 11744-11751.  
doi:10.1103/PhysRevB.52.11744</mixed-citation></ref><ref id="scirp.8339-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">J. S. McCalmont, M. M. Sigalas, G. Tuttle, K. M. Ho and C. M. Soukoulis, “A Layer-by-Layer Metallic Photonic Band- Gap Structure,” Applied Physics Letters, Vol. 68, No. 5, 1996, pp. 2759-2761. doi:10.1063/1.115589</mixed-citation></ref><ref id="scirp.8339-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">A. H. Aly and S.-W. Ryu, “Study of Optical Transmission in Nanometallic Photonic Crystal,” Journal of Computational and Theoretical Nanoscience, Vol. 5, No. 4, 2008, pp. 597-601. doi:10.1166/jctn.2008.022</mixed-citation></ref><ref id="scirp.8339-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">D. Nau, A. Sch?nhardt, C. Bauer, A. Christ, T. Zentgraf, J. Kuhl and H. Giessen, “Disorder Issues in Metallic Photonic Crystals,” Physica Status Solidi (B), Vol. 243, No. 10, 2006, pp. 2331-2343. doi:10.1002/pssb.200668054</mixed-citation></ref><ref id="scirp.8339-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">A. H. Aly, “The Transmittance of Two Types of One-Di- mensional Periodic Structures,” Materials Chemistry and Physics, Vol. 115, No. 1, 2009, pp. 391-394.  
doi:10.1016/j.matchemphys.2008.12.009</mixed-citation></ref><ref id="scirp.8339-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">A. H. Aly, S.-W. Ryu, H.-T. Hsu and C.-J. Wu, “THz Transmittance in One-Dimensional Superconducting Nano- material-Dielectric Superlattice,” Materials Chemistry and Physics, Vol. 113, No. 1, 2009, pp. 382-384. 
doi:10.1016/j.matchemphys.2008.07.123</mixed-citation></ref><ref id="scirp.8339-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">A. H. Aly, H.-T. Hsu, T.-J. Yang, C.-J. Wu and C.-K. Hwangbo, “Extraordinary Optical Properties of a Super- conducting Periodic Multilayer in near-Zero-Permittivity Operation Range,” Journal of Applied Physics, Vol. 105, No. 8, 2009, pp. 083917-083923. </mixed-citation></ref><ref id="scirp.8339-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">J. D. Joannopoulos, S. G. Johnson, J. N. Winn and R. D. Meade, “Photonic Crystals Moleding the Flow of Light,” Princeton University Press, Princeton, 2008.</mixed-citation></ref><ref id="scirp.8339-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">P. Yeh, “Optical Waves in Layered Media,” John Wiley and Sons, New York, 1988.</mixed-citation></ref><ref id="scirp.8339-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">M. Born and E. Wofl, “Principles of Optics,” Cambridge, London, 1999.</mixed-citation></ref></ref-list></back></article>