<?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">JAMP</journal-id><journal-title-group><journal-title>Journal of Applied Mathematics and Physics</journal-title></journal-title-group><issn pub-type="epub">2327-4352</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jamp.2015.36080</article-id><article-id pub-id-type="publisher-id">JAMP-57408</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Analytical Approach of the Nonlinear Surface Plasmon at a Left-Handed Material
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ajdi</surname><given-names>S. Hamada</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ali</surname><given-names>H. El-Astal</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mahmoud</surname><given-names>A. Sbiah</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mohammed</surname><given-names>M. Shabat</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Physics, The Islamic University of Gaza, Gaza Strip, Palestinian Authority</addr-line></aff><aff id="aff1"><addr-line>Department of Physics, Al-Aqsa University, Gaza, Gaza Strip, Palestinian Authority</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>shabat@iugaza.edu.ps(MMS)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>17</day><month>06</month><year>2015</year></pub-date><volume>03</volume><issue>06</issue><fpage>673</fpage><lpage>679</lpage><history><date date-type="received"><day>12</day>	<month>May</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>21</month>	<year>June</year>	</date><date date-type="accepted"><day>25</day>	<month>June</month>	<year>2015</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>
 
 
  The surface wave dispersion relations of surface Plasmon at the interface of a left-handed material and a non-linear Kerr medium of arbitrary nonlinearity are derived based on a generalized first integral approach. The normalized power flow is also investigated for various values of frequency. The above study is conducted for both cases: self-focusing (α≺0) and de-focusing (α≻0) nonlinear Kerr coefficient. 
 
</p></abstract><kwd-group><kwd>Left-Handed Materials</kwd><kwd> Nonlinearity</kwd><kwd> Dispersion Relation</kwd><kwd> Power Flow</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Recently there are great interest and investigation of plasmonics. This is due to the increasing of transmission in layered thin films composed of metals and their experimental applications. Surface Plasmons are charges oscillations occurring at the interface between metal and dielectric layers. Plasmonics concerns with the surface Plasmons and the light interaction with metals.</p><p>The interaction of light and surface Plasmon has increased many applications and investigation studies such as developed spectroscopy, high resolution microscopy and sensing, development of light sources and cloaking left-handed materials.</p><p>Surface Plasmon excitation concerns with the free electrons oscillation of the interface between metal and dielectric layers leading to the resultant excitation which depends of the optical properties of the two layers and the interface geometry between the two layers. The surface Plasmon excitation is also investigated at simple geometry interface where the left-handed material and the dielectric layers have an interface which is infinite planar [<xref ref-type="bibr" rid="scirp.57408-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.57408-ref2">2</xref>] .</p><p>In our study, a Kerr-nonlinear type dielectric has been investigated where the dielectric function of the dielectric media depends on the electric field intensity. The surface Plasmons dispersion equation at a planar interface between a metal layer and a linear optical layer (where the wave number is k and the angular frequency is ω) can be expressed as the following:</p><disp-formula id="scirp.57408-formula505"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x8.png"  xlink:type="simple"/></disp-formula><p>where c is the speed of light in vacuum and the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x9.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x10.png" xlink:type="simple"/></inline-formula> are dielectric functions of the dielectric and the metal, respectively, which can depend only on the frequency. Nonlinear dielectric-metal interface cases have also been investigated showing that the surface Plasmon can be excited by both TM and TE waves [<xref ref-type="bibr" rid="scirp.57408-ref3">3</xref>] - [<xref ref-type="bibr" rid="scirp.57408-ref5">5</xref>] . For such cases, the nonlinear surface Plasmon system is leading to the optical bistability, which has been observed recently in the transmission spectrum [<xref ref-type="bibr" rid="scirp.57408-ref6">6</xref>] and in the Goos-Handen shift [<xref ref-type="bibr" rid="scirp.57408-ref7">7</xref>] . It was also found that Equation (1) can be used with the dielectric function <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x11.png" xlink:type="simple"/></inline-formula> of the intensity-dependent medium (nonlinear medium) where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x12.png" xlink:type="simple"/></inline-formula> can be written as:</p><disp-formula id="scirp.57408-formula506"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x13.png"  xlink:type="simple"/></disp-formula><p>where ε is the linear part and α is a nonlinear dielectric coefficients, which frequency-dependent.</p><p>In recent years, there has been an increasing growing interest in new artificial metamaterials. One of the most important reasons is due to the unusual characteristics and behaviors. Some of new interesting application of metamaterial is to use the left handed material or metamaterial in construction optical wave guide sensors [<xref ref-type="bibr" rid="scirp.57408-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.57408-ref9">9</xref>] .</p><p>Instead of a semi-infinite metallic region, we study here the surface Plasmon dispersion relation of a left- handed material (LHM). LHM is a medium with negative permittivity, permeability, and refractive index, which was initially discussed by Veselago [<xref ref-type="bibr" rid="scirp.57408-ref10">10</xref>] in the microwave frequency and could restore the evanescent wave- field components, being thus termed a perfect lens. Nowadays, the efforts have been made to design LHMs for optical frequencies using metallic nanowires [<xref ref-type="bibr" rid="scirp.57408-ref11">11</xref>] . Additionally, it was shown that photonic crystals demonstrate negative refraction under some conditions [<xref ref-type="bibr" rid="scirp.57408-ref12">12</xref>] , thus resembling the most famous feature of LHMs. Many other features of LHM were studied in published papers [<xref ref-type="bibr" rid="scirp.57408-ref13">13</xref>] - [<xref ref-type="bibr" rid="scirp.57408-ref15">15</xref>] . In this work, we also used the correct method (first integral approach) to derive the dispersion relation and present numerical analysis on the limitation of the inexact treatment by comparing its results with the correct theory.</p></sec><sec id="s2"><title>2. Theory</title><p>ollowing the theory and approach of a TM-polarized wave which is considered to be propagated at the interface of a Kerr-type medium and a LHM, and by ignoring any loss in both media [<xref ref-type="bibr" rid="scirp.57408-ref16">16</xref>] . So we can easily find out the dispersive relation as follows; For simplicity, we take an isotropic Kerr medium as described in Equation (2), i.e.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x14.png" xlink:type="simple"/></inline-formula>. For the case of anisotropic media, it is rather mathematically tedious [<xref ref-type="bibr" rid="scirp.57408-ref13">13</xref>] . In TM waves (<xref ref-type="fig" rid="fig1">Figure 1</xref>), we have only y component for H Field, and the electric field in both the nonlinear medium (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x15.png" xlink:type="simple"/></inline-formula>) and the LHM (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x16.png" xlink:type="simple"/></inline-formula>) in general can be written as:</p><disp-formula id="scirp.57408-formula507"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x17.png"  xlink:type="simple"/></disp-formula><p>where c.c. is complex conjugate constant, and the relative phase of the two components <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x18.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x19.png" xlink:type="simple"/></inline-formula> is out of phase<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x20.png" xlink:type="simple"/></inline-formula>. The components of the amplitudes<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x21.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x22.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x23.png" xlink:type="simple"/></inline-formula> in the nonlinear Kerr medium were obtained from the Maxwell equations as following the approach in [<xref ref-type="bibr" rid="scirp.57408-ref16">16</xref>] have the form:</p><disp-formula id="scirp.57408-formula508"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x24.png"  xlink:type="simple"/></disp-formula><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Geometry of the problem</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1720305x25.png"/></fig><p>The solution of Maxwell Equations (4) in the LHM (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x26.png" xlink:type="simple"/></inline-formula>) which has a dielectric function<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x27.png" xlink:type="simple"/></inline-formula>, the x-com-</p><p>ponent electric field amplitude will have the form: <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x28.png" xlink:type="simple"/></inline-formula>where</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x29.png" xlink:type="simple"/></inline-formula>is selected in the form of the commonly used function for plasmon investigation , and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x30.png" xlink:type="simple"/></inline-formula> is constructed in an analogous form, i.e.,</p><disp-formula id="scirp.57408-formula509"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x31.png"  xlink:type="simple"/></disp-formula><p>The losses are neglected. Here <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x32.png" xlink:type="simple"/></inline-formula> is an effective plasma frequency that depends upon the geometry of the system, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x33.png" xlink:type="simple"/></inline-formula>is a resonance frequency and F is a parameter that also depends upon the system structure. For the case of LHM, Equation (4) can be transformed to</p><disp-formula id="scirp.57408-formula510"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x34.png"  xlink:type="simple"/></disp-formula><p>where D<sub>z</sub> is the z component of the electric induction vector D. Equation (6) hold in a linear medium with replaced ε<sub>d</sub> by ε<sub>eff</sub> .</p><p>The continuity of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x35.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x36.png" xlink:type="simple"/></inline-formula> across z = 0 then yields;</p><disp-formula id="scirp.57408-formula511"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x37.png"  xlink:type="simple"/></disp-formula><p>with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x38.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x39.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x40.png" xlink:type="simple"/></inline-formula> is the magnitude of the field.</p><p>A standard treatment of the nonlinear region 1 invokes a “first integral” to get at an equation for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x41.png" xlink:type="simple"/></inline-formula> [<xref ref-type="bibr" rid="scirp.57408-ref17">17</xref>] . Here, the technique is generalized to arbitrary nonlinearity. Differentiating Equation (4) with respect to z and multiplying through by<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x42.png" xlink:type="simple"/></inline-formula>, we get the following equation in region 1, as:</p><disp-formula id="scirp.57408-formula512"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x43.png"  xlink:type="simple"/></disp-formula><p>An integration with respect to z gives</p><disp-formula id="scirp.57408-formula513"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x44.png"  xlink:type="simple"/></disp-formula><p>where C is an integration constant. The key step is to recognize the identity</p><disp-formula id="scirp.57408-formula514"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x45.png"  xlink:type="simple"/></disp-formula><p>The Equation (9) can be rewritten as</p><disp-formula id="scirp.57408-formula515"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x46.png"  xlink:type="simple"/></disp-formula><p>Applying the boundary conditions, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x47.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x48.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x49.png" xlink:type="simple"/></inline-formula> must vanish as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x50.png" xlink:type="simple"/></inline-formula>, then the constant C equals zero. Equation (11) is valid for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x51.png" xlink:type="simple"/></inline-formula>.</p><p>Applying the boundary conditions at Equation (11)<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x52.png" xlink:type="simple"/></inline-formula>, then we have new notations as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x53.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x54.png" xlink:type="simple"/></inline-formula> to give:</p><disp-formula id="scirp.57408-formula516"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x55.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x56.png" xlink:type="simple"/></inline-formula> is the squared magnitude of the electric field at the interface between regions 1 and 2. Requiring this result to be consistent with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x57.png" xlink:type="simple"/></inline-formula> in Equation (4) gives:</p><disp-formula id="scirp.57408-formula517"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x58.png"  xlink:type="simple"/></disp-formula><p>Furthermore, using Equation (7) together with the relation<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x59.png" xlink:type="simple"/></inline-formula>, the following expressions can be derived as:</p><disp-formula id="scirp.57408-formula518"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x60.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x61.png" xlink:type="simple"/></inline-formula> is given from Equation (2) with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x62.png" xlink:type="simple"/></inline-formula>.</p><p>Now let us apply the above results to derive an explicit and exact dispersion relation for the surface Plasmon at a LHM-Kerr dielectric interface. Substituting Equation (14) into Equation (13), the result can finally be reduced to a quadratic equation in the wave number k which, leading to the following relation:</p><disp-formula id="scirp.57408-formula519"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x63.png"  xlink:type="simple"/></disp-formula><p>Putting the limits <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x64.png" xlink:type="simple"/></inline-formula> or<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x65.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x66.png" xlink:type="simple"/></inline-formula>in the Equation (15) leads to the Equation (1) as expected. By taking the derivative <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x67.png" xlink:type="simple"/></inline-formula> of the Equation (15) to be zero, one can find the surface Plasmon resonance frequency<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x66.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x68.png" xlink:type="simple"/></inline-formula>. This is equivalent to find the poles of Equation (15) and by writing Equation (15) in the form:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x69.png" xlink:type="simple"/></inline-formula>, one obtains<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x70.png" xlink:type="simple"/></inline-formula>. However, the poles of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x71.png" xlink:type="simple"/></inline-formula> is to be non-negative. Furthermore, one can show that the additional poles from <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x72.png" xlink:type="simple"/></inline-formula> will be given by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x73.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x74.png" xlink:type="simple"/></inline-formula>,</p><p>which are unacceptable since these will make k = 0 in Equation (15). Thus, we obtain the following implicit expression for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x75.png" xlink:type="simple"/></inline-formula> to satisfy from the poles of Equation (15):</p><disp-formula id="scirp.57408-formula520"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x76.png"  xlink:type="simple"/></disp-formula><p>Then, the surface plasmon frequency can be solved from Equation (16) to give the following relation:</p><disp-formula id="scirp.57408-formula521"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x77.png"  xlink:type="simple"/></disp-formula><p>This result is compared and contrasted with the one from the inexact approach [<xref ref-type="bibr" rid="scirp.57408-ref18">18</xref>] , namely,</p><disp-formula id="scirp.57408-formula522"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x78.png"  xlink:type="simple"/></disp-formula><p>Note that even in the weak field limit, Equation (17) implies that:</p><disp-formula id="scirp.57408-formula523"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x79.png"  xlink:type="simple"/></disp-formula><p>which contains an extra factor of 1/2 compared within Equation(18).</p><p>The power flow in the LHM linear medium is given by:</p><disp-formula id="scirp.57408-formula524"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x80.png"  xlink:type="simple"/></disp-formula><p>And the power flow in the nonlinear medium is written as:</p><disp-formula id="scirp.57408-formula525"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1720305x81.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. Results and Discussions</title><p>It has been noticed that the dispersion curves are obviously changed by the effect of the nonlinearity for both self-focusing (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x82.png" xlink:type="simple"/></inline-formula>) and de-focusing (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x83.png" xlink:type="simple"/></inline-formula>) nonlinear Kerr coefficient, as seen in <xref ref-type="fig" rid="fig2">Figure 2</xref>, related to the linear case where α = 0, i.e., in dielectric medium.</p><p>In <xref ref-type="fig" rid="fig3">Figure 3</xref>, Surface Plasmon frequency <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x84.png" xlink:type="simple"/></inline-formula> is plotted as a function of the nonlinearity (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x85.png" xlink:type="simple"/></inline-formula>). For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x86.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x87.png" xlink:type="simple"/></inline-formula>decreases monotonically with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x88.png" xlink:type="simple"/></inline-formula> while for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x89.png" xlink:type="simple"/></inline-formula> there is a cutoff values of (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x85.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x86.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x89.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x90.png" xlink:type="simple"/></inline-formula>) above which no</p><p>real solution of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x91.png" xlink:type="simple"/></inline-formula> exists. As can be noticed from <xref ref-type="fig" rid="fig3">Figure 3</xref>, the inset shows how the dispersion relation is modified by the field-dependent dielectric constant in a Kerr medium.</p><p>The normalized power <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x92.png" xlink:type="simple"/></inline-formula> for different values of frequency is plotted versus the wave vector k in for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x93.png" xlink:type="simple"/></inline-formula> (in <xref ref-type="fig" rid="fig4">Figure 4</xref>) and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x94.png" xlink:type="simple"/></inline-formula> (in <xref ref-type="fig" rid="fig5">Figure 5</xref>). As can be seen in <xref ref-type="fig" rid="fig4">Figure 4</xref> the curves show the behavior of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x95.png" xlink:type="simple"/></inline-formula> is reciprocal to that of the right-handed material. This behavior is because of the left-handed material layer where the power flow is in the opposite direction. While in <xref ref-type="fig" rid="fig5">Figure 5</xref>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x92.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x93.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x94.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x95.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x96.png" xlink:type="simple"/></inline-formula> the power flow shows two interesting behaviors, i.e., the non-reciprocal behavior and the bistability behavior. The bistability confirms that some of the power flow is in the positive direction while the other in the negative direction.</p></sec><sec id="s4"><title>4. Conclusion</title><p>It is noticed from our above treatment that the dispersion curves are strongly dependent on the intensity of the electric field. The most important conclusions are that for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x97.png" xlink:type="simple"/></inline-formula> the surface plasmon resonance frequency de-</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Dispersion curves for different electric field strength: E<sub>o</sub><sub> </sub><sub>= 0</sub> (curve 1) and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x99.png" xlink:type="simple"/></inline-formula> = 9 &#215; 10<sup>10</sup> V<sup>2</sup>/m<sup>2</sup>, (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x99.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x100.png" xlink:type="simple"/></inline-formula>, curve 2) and (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x99.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x100.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x101.png" xlink:type="simple"/></inline-formula>, curve 3)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1720305x98.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> The surface plasmon frequency versus the field intensity for different value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x103.png" xlink:type="simple"/></inline-formula> curve 1, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x103.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x104.png" xlink:type="simple"/></inline-formula> curve 2</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1720305x102.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> The normalized power P/P<sub>o</sub> for different values of frequency: f = 5.1 GHz (curve 1), f = 5.15 GHz (curve 2), and f = 5.19 GHz (curve 3) for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x106.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1720305x105.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> The normalized power P/P<sub>o</sub> for different values of frequency: f = 5.11 GHz (curve 1), f = 5.15 GHz (curve 2), and f = 5.17 GHz (curve 3) for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x108.png" xlink:type="simple"/></inline-formula></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1720305x107.png"/></fig><p>creases monotonically with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x109.png" xlink:type="simple"/></inline-formula> whereas for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x110.png" xlink:type="simple"/></inline-formula>, there is a cutoff value. It has been shown that the reversed power flow is due to the LHM material which is acting as a mirror for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x111.png" xlink:type="simple"/></inline-formula>, but for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x109.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x110.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x111.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1720305x112.png" xlink:type="simple"/></inline-formula> both re-</p><p>versal and bistability cases have been clearly observed. The two interesting behaviors could lead to new design of future application in Optoelectronic-Microwave technology.</p></sec><sec id="s5"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.57408-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Dionne, J.A., Verhagen, E., Polman, A. and Atwater, H.A. (2008) Are Negative Index Materials Achievable with Surface Plasmon Waveguides? A Case Study of Three Plasmonic Geometries. Optics Express, 16, 19001-19017.  
http://dx.doi.org/10.1364/OE.16.019001</mixed-citation></ref><ref id="scirp.57408-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Tsakmakidis, K.L., Hermann, C., Klaedtke, A., Jamois, C. and Hess O. (2006) Surface Plasmon Polaritons in Generalized Slab Heterostructures with Negative Permittivity and Permeability. Physical Review B, 73, Article ID: 085104.  
http://dx.doi.org/10.1103/PhysRevB.73.085104</mixed-citation></ref><ref id="scirp.57408-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Glass, N.E. and Rogovin, D. (1989) Surface-Polariton and Guided-Wave Excitation in Thin-Film Kerr Media. Physical Review B, 40, 1511. http://dx.doi.org/10.1103/PhysRevB.40.1511</mixed-citation></ref><ref id="scirp.57408-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Smolyaminoy, I.L., et al. (2002) Single-Photon Tunneling via Localized Surface Plasmons. Physical Review Letters, 88, Article ID: 187402. http://dx.doi.org/10.1103/PhysRevLett.88.187402</mixed-citation></ref><ref id="scirp.57408-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Baher, S. and Cottam, M.G. (2004) Theory of Nonlinear s-Polarized Phonon-Polaritons in Multilayered Structures. Journal of Science, 15, 171-177.</mixed-citation></ref><ref id="scirp.57408-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Wurtz, G.A., Potrald, R. and Zayats, A.V. (2006) Optical Bistability in Nonlinear Surface-Plasmon Polaritonic Crystals. Physical Review Letters, 97, Article ID: 057402.  http://dx.doi.org/10.1103/PhysRevLett.97.057402</mixed-citation></ref><ref id="scirp.57408-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Zhou, H., Chen, X., Hou, P. and Li, C.F. (2008) Giant Bistable Lateral Shift Owing to Surface Plasmon Excitation in Kretschmann Configuration with a Kerr Nonlinear Dielectric. Optical Letters, 33, 1249-1251.  
http://dx.doi.org/10.1364/OL.33.001249</mixed-citation></ref><ref id="scirp.57408-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Taya, S.A., Shabat, M.M. and Khalil, H.M. (2009) Enhancement of Sensitivity in Optical Waveguide Sensors Using Left-Handed Materials. Optik, 120, 504-508. http://dx.doi.org/10.1016/j.ijleo.2007.12.001</mixed-citation></ref><ref id="scirp.57408-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Taya, S.A., Shabat, M.M. and Khalil, H.M. (2008) Analysis of TM Nonlinear Asymmetrical Waveguide Optical Sensors. Sensors and Actuators A: Physical, 147, 137-141. http://dx.doi.org/10.1016/j.sna.2008.05.002</mixed-citation></ref><ref id="scirp.57408-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Veselago, V.G. (1968) The Electrodynamics of Substances with Simultaneously Negative Values of ε and μ. Soviet Physics Uspekhi, 10, 509-514. http://dx.doi.org/10.1070/PU1968v010n04ABEH003699</mixed-citation></ref><ref id="scirp.57408-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Podolskiy, V.A., Sarychev, A.K. and Shalaev, V.M. (2002) Plasmon Modes in Metal Nanowires and Left-Handed Materials. Journal of Nonlinear Optical Physics &amp; Materials, 11, 65-74. http://dx.doi.org/10.1142/S0218863502000833</mixed-citation></ref><ref id="scirp.57408-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Notomi, M. (2000) Theory of Light Propagation in Strongly Modulated Photonic Crystals: Refractionlike Behavior in the Vicinity of the Photonic Band Gap. Physical Review B, 62, 10696-10705. 
http://dx.doi.org/10.1103/PhysRevB.62.10696</mixed-citation></ref><ref id="scirp.57408-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Hamada, M.S., El-Astal, A.H. and Shabat, M.M. (2007) Nonlinear TE Surface Waves in a Photosensitive Semiconductor Film Bounded by a Superconductor Cover. International Journal of Modern Physics B, 21, 1817-1825. 
http://dx.doi.org/10.1142/S0217979207037065</mixed-citation></ref><ref id="scirp.57408-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Hamada, M.S., El-Astal, A.H. and Shabat, M.M. (2007) Characteristic of Surface Waves in Nonlinear Left-Handed- Photosensitive-Semiconductor Waveguide Structure. International Journal of Modern Physics B, 21, 5319-5329. 
http://dx.doi.org/10.1142/S0217979207038459</mixed-citation></ref><ref id="scirp.57408-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Hamada, M.S., Ass’ad, A.I., Ashour, H.S. and Shabat, M.M. (2006) Nonlinear Magnetostatic Surface Waves in a Ferrite-Left-Handed Waveguide Structure. Journal of Microwaves and Optoelectronics, 5, 45-54.</mixed-citation></ref><ref id="scirp.57408-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Mihalache, D., Nazmitdinov, R.G. and Fedyanin, V.K. (1989) Nonlinear Optical Waves in Layered Structure. Soviet Journal of Nuclear Physics, 20, 86-107.</mixed-citation></ref><ref id="scirp.57408-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Economou, E.N. (1969) Surface Plasmons in Thin Films. Physical Review, 182, 539-554. 
http://dx.doi.org/10.1103/PhysRev.182.539</mixed-citation></ref><ref id="scirp.57408-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Waysin, G.M., Simon, H.J. and Deck, R.T. (1981) Optical Bistability with Surface Plasmons. Optics Letters, 6, 30-32. 
http://dx.doi.org/10.1364/OL.6.000030</mixed-citation></ref></ref-list></back></article>