Achieving transparency with plasmonic and metamaterial coatings

Size: px
Start display at page:

Download "Achieving transparency with plasmonic and metamaterial coatings"

Transcription

1 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Ahievig trasparey with plasmoi ad metamaterial oatigs Adrea Alù,2, ad Nader Egheta,* Dept. of Eletrial ad Systems Egieerig, Uiversity of Pesylvaia, Philadelphia, 94 PA, USA 2 Dept. of Applied Eletrois, Uiversity of Roma Tre, Rome, Italy Abstrat: The possibility of usig plasmoi ad metamaterial overs to drastially redue the total satterig ross setio of spherial ad ylidrial objets is disussed. While it is ituitively expeted that ireasig the physial size of a objet may lead to a irease i its overall satterig ross setio, here we see how a proper desig of these lossless metamaterial overs ear their plasma resoae may idue a dramati drop i the satterig ross setio, makig these objet early ivisible or trasparet to a outside observer a pheomeo with obvious appliatios for low-observability ad o-ivasive probe desig. Physial isights ito this pheomeo ad some umerial results are provided. PACS umbers: a, e, 42.5.Gy, 33.2.Fb A. Alù, ad N. Egheta, Physial Review E, 72, 6623, 25. Copyright (25) by the Ameria Physial Soiety. INTRODUCTION Ahievig ivisibility or low observability has bee the subjet of extesive studies i the physis ad egieerig ommuities for deades. The use of absorbig srees (e.g., []) ad ati-refletio oatigs (e.g., [2]) to dimiish the satterig or the refletio from objets, for example, are ommo i several appliatios. If the former require absorptio, ad therefore power dissipatio, the latter are suitable primarily for plaar or early-plaar trasparet objets. I this paper, however, we disuss how low-loss (i the limit eve o-loss) passive overs might be utilized i order to drastially redue the satterig from spherial ad ylidrial objets without requirig high dissipatio, but relyig o a ompletely differet mehaism. For this purpose, we osider here the use of materials with egative or low eletromageti ostitutive parameters, e.g., ertai metals ear their plasma frequey or metamaterials with egative parameters. Several oble metals ahieve this requiremet for their eletrial permittivity at the ifra-red (IR) or visible

2 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE regimes, eve with reasoably low losses [3]-[4]. At lower frequeies, moreover, artifiial materials ad metamaterials may be properly sythesized to meet similar requiremets at the desired frequey (see referees i [5]). At mirowaves, for istae, the tehology of artifiial materials with uusual eletromageti properties has a log traditio i the egieerig ommuity. Materials with low-positive, ear-zero or egative relative permittivity have bee sythesized by embeddig arrays of thi wires i a host material [6]-[8], ad their properties have bee exploited i differet mirowave appliatios over the years, e.g., for highly diretive ateas [9]-[] or more reetly for aomalous tuelig ad trasmissio (e.g., []-[2]). By a similar priiple, arrays of resoatig loops or split-rig resoators affet i a aalogous way the permeability of the bulk material ad have bee studied i the last deade [3]-[4]. The ombiatio of these ilusios i omposite metamaterials with a egative idex of refratio [5] has reetly attrated a great deal of attetio i the sietifi ommuity. The idea of utilizig low-permittivity oatigs to idue ivisibility has bee suggested i the quasi-stati limit for spherial objets (i.e., for the ase of a sphere beig small ompared with the operatig wavelegth) ad for the higher-order terms i power series expasio [6, p. 49-5], [7]-[9]. Here, however, we explore fully the ase of plasmoi ad metamaterial overs with egative or low-permittivity ad/or permeability materials i the dyami full-wave satterig ase (i.e., the sphere does ot eed to be eletrially very small), providig some ew physial isights ad possibility of optimized desigs for reduig the total satterig ross setio of spherial objets whose dimesios are omparable with the wavelegth of operatio. THEORETICAL ANALYSIS Let us osider a spherial satterer of radius a, omposed of a homogeeous material with permittivity ε ad permeability μ, surrouded by free spae (with ostitutive parameters ε ad μ ). (As will be metioed later, i the limit this sphere a also be a perfet eletri odutor.) Would it be possible to over suh a objet with a lossless (or low-loss) spherial shell (with ε, μ ad radius a > a) i suh a way that the total satterig ross setio drastially dimiishes, eve though its overall physial size has * To whom orrespodee should be addressed, egheta@ee.upe.edu 2

3 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE evidetly ireased? Here we explore the oditios uder whih this may be possible. The geometry of the problem is depited i Fig.. Let the illumiatig wave be represeted by a i t e ω moohromati plae wave. All the ostitutive parameters may i geeral be omplex at the frequey of iterest ω, takig ito aout possible losses. Adoptig a stadard Mie expasio, oe a write the sattered field from a plae wave as the well-kow disrete sum of spherial harmois with omplex amplitudes TE ad ( beig a iteger), respetively for the TE ad spherial waves (see e.g. [2]). Their expressios for the problem at had may be writte i the oveiet form [2]: U TE TE = TE TE U + iv, U = U + iv, () where the futios U ad V are real valued i the limit of o ohmi (material) losses, ad i the polarizatio are give by the determiats: U = ( ) ( ) ( ) j ka j k a y k a kaj ka kaj ka kay ka ( ) / ε ( ) / ε ( ) / ε j ( k a ) y ( k a ) j ( k a ) ka j( ka ) / ε ka y( ka) / ε ka j( ka) / ε, (2) V = ( ) ( ) ( ) j ka j k a y k a kaj ka kaj ka kay ka ( ) / ε ( ) / ε ( ) / ε j ( k a ) y ( k a ) y ( k a ) ka j( ka) / ε ka y( ka ) / ε ka y( ka ) / ε. (3) Aalogous expressios for the TE polarizatio are obtaied after substitutig ε with μ i (2) ad (3). I these formulas k ω εμ, k ω εμ ad k ω ε μ are the wave umbers i the three regios ad (). j, (). y are spherial Bessel futios. The satterig oeffiiets are related to the total satterig ross setio through the formula: 2π Q = + +. (4) TE 2 2 ( 2 )( ) s 2 k = Geerally speakig, the satterig ross setio Q s is determied by the satterig oeffiiets up to the order = Nmax, sie the suessive satterig oeffiiets beyod max N will be egligible [6]. The value 3

4 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE of N N usually ireases with the physial/eletrial size of the satterer (a rule of thumb is that max ( k a ) [6, p. 26]), ad that is oe of the reasos that larger objets geerally have a wider max satterig ross setio. Followig this observatio, oe may expet, ituitively, that overig a objet would irease its satterig ross setio, together with its physial/eletrial size, sie more satterig oeffiiets otribute i the summatio. However, as will be show below, this statemet is ot always true, ad with judiiously hoosig the over material ad its size, oe may be able to make a give homogeous (or eve metalli) sphere less detetable to a exteral observer. Let us first oetrate o the ase i whih the satterer is small eough to osider oly the dipolar term ( = ) i the Mie expasio. This ase is osistet with what was doe before i [6]-[9]. I this ase, i.e., whe ka, ka, ka, the expressios i (2) ad (3) are redued to the followig [2]: where U + π ( ka ) ( + / ) ε ( + / ) ε / ε, (5) 4 ( 2+ )( /2 )! γ γ γ / ε γ / ε / ε V ( + ) + ( + ) ( + ) ( + / ) ε ( + / ) ε / ε ( k a ), (6) γ γ γ ε γ ε ε + / / / γ a/ a is the ratio of ore-shell radii. These losed-form expressios reveal iterestig properties for eletrially small satterers. Firstly, as expeted, for small satterers the value of U is small ad teds to zero as a. This is osistet with the fat that usually a small satterer has a very low sattered field. I fat, osiderig stadard dieletris far from their plasma frequey (all permittivities greater or equal to TE ε, all permeabilities equal to μ ) ad the sum i the deomiator of is domiated by V, yieldig the followig approximate expressio for : 2+ ( ) ( ) j k a f γ, (7) where f ( γ ) are positive real futios of γ. This ofirms that for suh tiy dieletri satterers the satterig properties are domiated by the first-order term, orrespodig to the radiatio from a eletri dipole. This also shows that overig this struture with a dieletri shell may result i ireasig its 4

5 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE satterig ross setio together with its physial/eletrial size (i.e., a ). However, there is oe possibility for whih addig a over may lead to a drasti redutio of the total satterig ross setio from a small sphere, ad that is whe the over parameters are hose suh that the determiat U i (5) may approah zero, resultig i, whih is the domiat term i (4), to beome zero. I the limit for whih (5) applies, its value will be idetially zero for a give order, if γ satisfies the followig oditio: γ ( ) ( ) ( ) ( ) ε ε + ε + ε = 2+ T ε ε + ε+ ε. (8) where T is a shorthad for the ( ) 2 + -st root. Clearly, i order for γ to have physial meaig, its value should satisfy γ. This result is osistet with the fidigs i the quasi-stati ase i [6]-[9]. I Fig. 2 the values of permittivities for whih (8) a be fulfilled are reported i a otour plot, with lighter regios idiatig higher values of T. Whe γ = T for a give, the orrespodig satterig oeffiiet beomes idetially zero (i the small-sphere approximatio). Similar results may be obtaied for the TE polarizatio by replaig the permittivities with the orrespodig permeabilities. Oe may otie that if both materials are stadard materials far from their plasma frequeies (i.e., with ostitutive parameters higher tha those of free spae) it is ot possible to ahieve suh a oditio for ay value of γ. If however we are allowed to use overs with permittivity or permeability lower tha the oe i free spae, or with egative values, whih might represet metals lose to their plasma frequey or alteratively low-permittivity, low-permeability, or egative-parameter metamaterials, the followig Fig. 2 a proper hoie of the ratio of ore-shell allows brigig to zero oe partiular. If the ore sphere (the oe to be hidde) is dieletri (ad eletrially small, i order to apply (8)), the the eletri dipole term is the domiat term i satterig, ad all the other terms i (4) are egligible. I this ase, from (8) the oditio for the over radius beomes: a = 3 ( ε ε )[ 2ε+ ε] ( ε ε)[ 2ε + ε ] a, (9) whih is osistet with the formula derived i [6, p. 5] i the stati ase. This applies, of ourse, oly whe the pair of permittivities ε, ε falls i allowable regios i Fig. 2. As a speial ase of iterest, we 5

6 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE may osider the ase i whih the ore sphere is perfetly odutig, for whih get a 3 ( ε ε )/( 2ε ε ) = + a, whih may be satisfied whe ε i. I this ase, we ε is hose i the iterval < ε < ε. It may be oted here that with some other ombiatios of parameters it is possible to derive aother ratio γ for whih the determiat V i (6) beomes zero. As already disussed i our previous work [2]-[22], this ratio would idue a aomalous plasmoi resoae i the ore-shell system, ausig a opposite effet, amely a highly ireased satterig from a eletrially small ore-shell system. Whe the size of the sphere ireases, expressio (9) is o more effetive i determiig the exat radius of the outer ore to ahieve the trasparey oditio. I this ase, however, it is still possible to ahieve aother oditio for whih U =, whih happes at the partiular γ that makes the determiat i (2) zero. This oditio is effetive for bigger dimesios of the sphere for whih the quasi-stati limit does ot apply, eve though this does ot eessarily imply that the satterig from the sphere will beome zero, sie U = is for the dipolar term to be zero, ad for too large spheres higher order terms i (4) may beome more sigifiat tha the first-order term. I this ase, however, it may be possible to exploit a similar approah with a multi-layer over, whih should provide more degrees of freedom, ad thus a proper desig might make simultaeously zero some of the higher-order terms that otribute otieably to the satterig of the partile. This is oe of the subjets i this area we are urretly ivestigatig. It is iterestig to ote that it may be possible to derive aalogous theoretial results for other aoial geometries. I the ylidrial ase, for istae, followig similar steps it may be show that the trasparey oditio for the TE (with respet to the ylider axis, i.e., with the mageti field parallel with the ylider axis) polarizatio beomes: γ γ ( ε ε )( ε+ ε) ( ε ε)( ε + ε ) = 2 μ μ μ μ 2 = = 2 for for, () where γ represets agai the ratio of ore-shell radii, ad the formulas may be derived by duality. Other geometrial shapes ad referee systems would allow obtaiig aalogous oditios, whih are valid i the quasi-stati approximatio, ad may be similarly exteded to the full-wave aalysis. 6

7 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE NUMERICAL RESULTS AND DISCUSSION A importat questio is how large the radius a may beome i the geeral full-wave satterig ase while oe would still be able to ahieve a drasti redutio of satterig ross setio with a sigle over, aelig just the first order term. For this purpose, Fig. 3 shows the full-wave total satterig ross setio of a dieletri partile ( ε = 4ε ) overed by a plasmoi or metamaterial over with ε = 3ε, i terms of the ratio a/ a, for differet sizes of the ore-shell partile. For this ombiatio of parameters the value of T i (8) for = is give by T.6. As a be see from this Figure, whe the outer radius of the overed partile is extremely small ( a = λ /, λ = 2 π /k beig the free-spae wavelegth), its satterig ross setio goes to ear zero very lose to that ratio. Ireasig the partile size ( a = λ / ), the ross setio a still beome very small, makig the objet early ivisible, (sie the higher order satterig oeffiiets are still egligible), eve though the quasi-stati solutio is o more adequate ad the trasparey oditio should be obtaied diretly by equatig Eq. (2) to zero, shiftig dowwards the required value of γ for the trasparey i this ase. Eve with a larger outer radius ( a = λ /5 ) we may get a very low satterig, but of ourse the value of γ for the miimum satterig has bee shifted dowwards agai. I order to give a physial isight ito this trasparey pheomeo, oe may thik of it i the quasi-stati approximatio ad with the goal of simply reduig the dipolar term. As is well kow, the dipolar satterig is the mai respose related to the polarizatio vetor P= ( ) ε ε E idued loally by the loal field as a result of exitig eletri field. A regular dieletri sphere would show a sattered dipolar field due to the total idued eletri dipole momet, whih i the quasi-stati limit orrespods to the itegral of the vetor P over the volume of the sphere. A over with ε < ε, however, would show a loal polarizatio vetor P= ( ε ε ) E i the over, beig ati-parallel to the loal eletri field i the over, whih after beig itegrated over the shell volume may ael the origial dipole momet of the ore itself. With the proper hoie of radius a for the over this aellatio may be omplete, ad this is what formula (9) 7

8 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE heuristially represets. This is summarized i Fig. 4 (a similar ituitive explaatio has bee provided i [23] for the stati seario). We ote that this aellatio pheomeo does ot rely o a resoat pheomeo, ulike the opposite pheomeo of the satterig resoae by suh a ore-shell system (for whih V = ) [2]-[22]. This is a iterestig poit, sie it esures: a) that this pheomeo is ot strogly sesitive to losses or to the spherial symmetry of the partile, ad therefore ohmi losses or imperfetios i the partile shape will ot preipitously alter this pheomeo, ad b) that the pheomeo may ot be strogly sesitive to the ratio oditio (8) (i.e., it has a larger ratio badwidth ), as a also be see i the examples of Fig. 3. Usig a similar argumet, we may apply the trasparey pheomeo to a plasmoi partile ( ε = 3ε, a = λ / ) by oatig it with a dieletri material ( ε = ε ). For this ombiatio T.825 ad therefore the required radius for the over is a.2a i order to ahieve trasparey. This is show i Fig. 5, where the satterig ross setio (a) ad the otributios from the several satterig oeffiiets (b) are depited varyig the over radius. It is iterestig to otie how lose to the ear-zero-satterig ratio we fid a peak i the satterig ross setio (Fig. 5a), whih is learly due to the resoae of the 2 satterig oeffiiet (Fig. 5b). This is expeted whe V 2 = i Eq. (3), whih i the quasi-stati approximatio is expeted for a.5a. It a also be otied that the trasitio i whih approahes zero i Fig. 5b is a relatively smooth oe (very differet form the satterig peak, as disussed i [2]), implyig that the rage of the ratio of radii i Fig. 5a for whih the satterer has a very low satterig ross setio is relatively broad. As was just poited out, this is aother idiatio of the o-resoat behavior of this pheomeo. Fig. 6 reports the same plots as i Fig. 5, but for a larger partile ( a = λ /5 ). The results show that eve for a partile far from the quasi-stati oditio, the satterig ross setio may be drastially redued with the same priiple. (Oe may appreiate how the higher order terms i (4) are sigifiat i Fig. 6b.) The effetiveess of the over is due to the fat that the higher order satterig terms are also i geeral affeted whe suh a omplemetary over is employed (Fig. 6b). It should be uderlied that the ombiatios of parameters hose i Figs. 4, 5 ad 6 also allow the 8

9 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE presee of the resoae peaks desribed i [2] (at whih the differet V may go to zero) for other ratios of radii (they are learly see i Figs. 5b ad 6b). For larger partiles the presee of this peak worses the trasparey performae of these plasmoi ad metamaterial overs: the miimum dip we see i the plot of satterig ross setio is losely surrouded by the regios where the resoat maxima are preset, ad as a result a good trasparey oditio may ot be effetively ahieved for larger partiles (at least with a sigle-layer homogeeous over with egative-parameter metarials). Therefore, i order to ahieve better trasparey performae for partiles that do ot satisfy the quasi-stati approximatio, the trasparey ratio T ad the high-satterig resoat ratios should be as far apart as possible. There is ideed a rage of material parameters i Fig. 2 where this oditio is possible: whe oe of the two materials has a positive permittivity, but lower tha the permittivity of the outside regio (e.g., free spae). I this rage of material parameters, it a ideed be show that o resoae peak may arise i the quasi-stati approximatio (the determiat i (6) aot be zero for material parameters i this rage [2]). Fig. 7, as a example, shows the results for a sphere with ε = ε ad a = λ / overed with a material with ε =.5ε, as a futio of a / a. I this ase, as a be see, the redutio of total satterig ross setio a also be ahieved, ad it is learly evidet from Fig. 7b that o high-satterig resoat peak is aroud this trasparey poit. The total satterig ross setio of larger objets may be suessfully redued by further reduig the positive permittivity of the over. Fig. 8 shows some results for the same ier material as i Fig. 7, but twie its radius (ad so eight times its volume). The over utilized ow has ε =.ε. We ote how the satterig ross setio has bee drastially redued at the optimum radius for whih =. The residual satterig ross setio at the miimum dip i Fig. 8a is maily due to the otributio of the TE mageti dipole otributio, sie this is ot diretly affeted by the low permittivity of the over, as a be see from Fig. 8b. Provided that we may fid (or sythesize, as disussed i the itrodutio) a material whose permeability may be brought dow towards zero at the same frequey, however, oe might be able to further redue the satterig from the objet usig a similar priiple. Choosig a over TE permeability as μ =.25μ, for whih U = at the same a =.2a where we get the miimum i Fig. 8a, we are able to make the objet eve less visible, as show i Fig. 9. 9

10 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Fig. shows the otour plots of the magitude of the radial ompoet of the ear-zoe sattered eletri field for the ase of a dieletri sphere with ε = ε, μ = μ ad a = λ /5 i four differet ases: o over (Fig. a), a dieletri over with ε = 5ε ad a =.2a (Fig. b), the same over as i Fig. 8 (Fig. ) ad the same over as i Fig. 9 (Fig. d). I all these ases the system is assumed to be exited by a plae wave propagatig alog the ẑ diretio with a uit-amplitude eletri field polarized alog ˆx, ad the plots are i the x-y plae, where the eletri multipole otributio is maximum. We ote that the otour plots are all ormalized to the same value, so we a ompare the field amplitudes i the differet figures. As a be easily see, i the first ase without over (Fig. a) the eletri field idues a strog eletri dipole momet iside the sphere, whih osequetly geerates a dipolar sattered field i the outside regio. This is the usual situatio. Addig a over made of a regular dieletri with arbitrarily seleted ε = 5ε the sattered field is further ehaed (Fig. b), as expeted. The sattered radiatio beomes eve less symmetri due to the otributio of the quadrupole ad higher order multipoles. However, whe the over used i Fig. 8 with ε =.ε, μ = μ ad a =.2a is used (Fig. ), a sattered field omparable with the previous ases (a ad b) is observed iside the objet, but the over almost aels the sattered radiatio outside. The residual sattered radiatio (whih is very weak ompared with the sattered radiatio i Fig. a or b) is learly due to the quadrupolar ad otupolar terms, but it has bee drastially redued i magitude, osistet with what was show i Fig. 8a. Whe the over of Fig. 9 (with ε =.ε, μ =.25μ ad a =.2a) is employed, the sattered field is eve further redued, eve though the effet of this over is maily see i the mageti field plots (Fig. ), sie it has bee desiged to redue the mageti dipole radiatio. Fig. shows the orrespodig plots similar to Fig., but for the magitude of the radial ompoet of the ear-zoe sattered mageti field, i the y-z plae where the mageti multipoles otribute maximally to the radial ompoets. Fig. a shows the radial ompoets of the mageti field sattered by a sigle dieletri sphere, ad it is domiated by the mageti dipole otributio. Ireasig the size by addig aother dieletri over, as show i Fig. b, produes a irease i the sattered field, ad ow the mageti quadrupole beomes also otieable. Whe the over of Fig. is used here, istead, the mageti ear field is similar to the oe obtaied without ay over, as expeted from the results of Fig. 8. However,

11 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE addig the low-permeability oditio to the over, as show i Fig. d, the mageti dipole otributio to the sattered field is essetially aeled. Now the objet is almost ivisible i both the E ad H plaes, i.e., i the whole spae. Fig. 2 shows the total time-averaged Poytig vetor distributio i the x-z plae for the same four ases. Whe the trasparey overs are used (Figs. 2, 2d), oe otes how the iomig plae wave seems to pass through the sphere without beig affeted otieably by the sattered field. I other words, a observer seated outside the overed sphere, eve i the ear field, would ot be able to essetially sese the presee of the sphere (sie there is almost o sattered field outside), although the fields are otieably haged iside the sphere ad its over. If a ovetioal satterer is plaed behid suh a trasparet sphere, a observer would iterestigly otie the presee of that satterer without essetially seeig the trasparet sphere betwee him/her ad the satterer itself. Fig. 3 reports the similar results as i Fig., but for a smaller sphere ( a = λ /2 ). I this ase, the sattered fields from the sphere aloe ad from the sphere overed with a stadard dieletri (Figs. 3a ad 3b) are domiated by the eletri dipole term, ad the mageti dipole otributio is egligible. Coverig it with the proper over (sie we are i this quasi-stati ase we may simply use the approximate expressio (9), whih yields for the radius a =.9a), the trasparey is effetively ahieved. Due to its small dimesios, there is o eed to employ low- μ materials, sie the H-plae field otributio is already egligible. The related plots therefore are ot reported here. We eed to poit out the sesitivity/robustess of the results preseted here to some realisti parameters. As disussed previously, the trasparey pheomeo is ot due to a resoae, but just relies o the overall aellatio of the multipolar satterig fields (see Fig. 4). This is why this trasparey effet has a relatively broader rage for seletio of ratio of radii ad over sizes (we are ot restrited to oe sigle value for the over, but there is a relatively broader rage of values for whih a over with egative polarizability may hide a dieletri sphere). By the same toke, this pheomeo may be robust with the variatio of the material losses. Fig. 4, as a example, shows how the results of Fig. 7a ad Fig. 8a are modified by the presee of losses i the over material. As a be see, the trasparey effet is ot affeted muh by the presee of reasoable

12 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE ohmi (material) losses. For the dispersio with the Drude or Loretzia models, i the frequey regios ear the plasma frequeies where the ostitutive parameters have values lose to zero (whih we have show to be more effetive i iduig this trasparey effet) there are o sharp resoaes i the medium dispersio, ad therefore, followig the Kramers-Kroig relatios, we expet less ohmi losses tha ear the resoat frequeies [4]. We speulate that this setup is ot too sesitive to possible imperfetios i the ostrutio ad geometry of the over. Beig a o-resoat pheomeo, we do ot expet substatial hages i the results reported here for quasi-spherial or quasi-ylidrial objets with small imperfetios o their surfae. The fat that a egative polarizability material would effetively ael oe or more multipolar otributios of the sattered field is ideed ot a property of the spherial or the ylidrial geometry, ad it is preditable that formulas aalogous to (9) ad () may be derived as well for several other geometries ad for more omplex objets. As a fial remark, we should poit out that this trasparey pheomeo may somehow be related to the issues of o-radiatig soures i the iverse satterig ad soure problems i eletromagetis ad aoustis [24]-[25]. We are urretly ivestigatig this poit. SUMMARY I this otributio we have studied how it is possible to drastially redue the satterig ross setio of spherial ad ylidrial objets usig lossless plasmoi or metamaterial overs. We have provided physial isights ad umerial examples of how a proper desig of these metamaterial overs ear their plasma resoae may idue a dramati drop of the satterig ross setio (eve i the absee of material loss) makig the objet early trasparet to a exteral observer. We have also disussed how this effet is isesitive to the possible losses or other imperfetios i strutures. This pheomeo a provide umerous potetial appliatios i the desig of low-observable targets, low-ouplig i desely-paked devies, ad o-ivasive field ao-probes. These issues are urretly uder ivestigatio. ACKNOWLEDGEMENTS This work is supported i part by the U.S. Defese Advaed Researh Projets Agey (DARPA) Grat Number HR-4-P-42. Adrea Alù has bee partially supported by the 24 SUMMA Graduate 2

13 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Fellowship i Advaed Eletromagetis. We thak Dr. W. Ross Stoe for poitig us to the referees [24] ad [25] ad the topi of iverse satterig problems ad Dr. Alexey Viorgradov for brigig to our attetio Ref [7]. REFERENCES [] R. L. Fate, ad M. T. MCormak, IEEE Tras. Ateas Propag., 3, 443 (968). [2] J. Ward, Vauum, 22, 369 (972). [3] P. B. Johso, ad R. W. Christy, Phys. Rev. B, 6, 437 (972). [4] L. Ladau, ad E. M. Lifshitz, Eletrodyamis of otiuous media (Elsevier, 984). [5] R. W. Ziolkowski, ad N. Egheta (guest editors), IEEE Tras. Ateas Propag., Speial Issue o Metamaterials, 58 (23). [6] W. Rotma, IRE Tras. Ateas Propag.,, 82 (962). [7] J. B. Pedry, A. J. Holde, D. J. Robbis, ad W. J. Stewart, J. Phys.: Cod. Matter,, 4785 (998). [8] A. N. Lagarkov, S. M. Matytsi, K. N. Rozaov, ad A. K. Saryhev, J. Appl. Phys., 84, 386 (998). [9] J. Brow, Pro. IEE, IV, 5, (953). [] I. J. Bahl, ad K. C. Gupta, IEEE Tras. Ateas Propag., 22, 9 (974). [] A. Alù, ad N. Egheta, IEEE Tras. Ateas Propag., 5, 2558 (23). [2] A. Alù, F. Bilotti, N. Egheta, ad L. Vegi, Pro. of IEE Semiar o Metamaterials for Mirowave ad (Sub) millimetre Wave Appliatios, Lodo, UK, pp. /-/6, November 24, 23. [3] J. B. Pedry, A. J. Holde, D. J. Robbis, ad W. J. Stewart, IEEE Tras. Mirowave Theory Teh, 47, 275 (999). [4] P. Gay-Balmaz ad O. J. F. Marti, Appl. Phys. Lett., 8, 939 (22). [5] R. A. Shelby, D. R. Smith, ad S. Shultz, Siee, 292, 77 (2). [6] C. F. Bohre, ad D. R. Huffma, Absorptio ad Satterig of Light by Small Partiles (Wiley, New York, 983). [7] Z. Hashi, ad S. Shtrikma, J. Meh. Phys. Solids,, 27 (963). [8] M. Kerker, J. Opt. So. Am., 65, 376 (975). [9] H. Chew, ad M. Kerker, J. Opt. So. Am., 66, 445 (976). 3

14 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE [2] J. A. Stratto, Eletromageti Theory (MGraw-Hill Comp., 94). [2] A. Alù, ad N. Egheta, J. Appl. Phys., 97, 943 (2 pages) (25). [22] A. Alù, ad N. Egheta, Pro. of 24 IEEE Radar Coferee, Philadelphia, USA, April 26-29, 24, pp [23] M. Kerker, Aerosol Siee ad Tehology,, 275 (982). [24] N. Bleistei ad J. Cohe, J. Math. Phys. 8, 94 (977) [25] B. J. Hoeders, i Iverse Soure Problems i Optis, H. P. Baltes (ed.), Vol. 9 of Topis i Curret Physis, Spriger, Berli, 4 (978) FIGURES Figure Cross setio of a spherial satterer omposed of two oetri layers of differet isotropi materials i a suitable spherial referee system. Figure 2 Rages of material permittivities for whih the oditio (8) admits physial solutios ad otour plot for the values of the ratio T ( polarizatio). 4

15 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Figure 3 Normalized satterig ross setio for a spherial partile with ε = 4ε, μ = μ, ad the over with ε = 3ε, μ = μ, versus the ratio a/ a for three differet sizes of the outer radius of the over. Figure 4 Heuristi iterpretatio of the trasparey pheomeo: aellatio of the overall dipole 5

16 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE momet through a idued egative polarizatio vetor. Figure 5 (a) Normalized total satterig ross setio ad (b) otributios of the several (o-zero) satterig oeffiiets, i terms of the over radius for a fixed ier ore ( ε = 3ε, ε = ε, μ = μ = μ, a = λ / ). 6

17 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Figure 6 The same as Figure 5, but a = λ /5 ad higher order terms beome sigifiatly higher. But still, a sigifiat redutio of the satterig ross setio may be ahieved. 7

18 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Figure 7 (a) Normalized total satterig ross setio ad (b) otributios of the several (o-zero) satterig oeffiiets, i terms of the over radius for a fixed ier ore ( ε = ε, ε =.5ε, μ = μ = μ, a = λ / ). 8

19 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Figure 8 (a) Normalzied total satterig ross setio ad (b) otributios of the several (o-zero) satterig oeffiiets, i terms of the over radius for a fixed ier ore ( ε = ε, ε =.ε, μ = μ = μ, a = λ /5 ). 9

20 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Figure 9 (a) Normalized total satterig ross setio ad (b) otributios of the several (o-zero) satterig oeffiiets, i terms of the over radius for a fixed ier ore ( ε = ε, ε =.ε, μ = μ, μ =.25μ, a = λ /5 ). Figure Cotour plots of the distributio of the magitude of radial ompoet of the sattered eletri field i the x-z plae, idued by a plae wave travelig alog the z diretio with a eletri field polarized alog the x axis: (a) for a sphere with ε = ε, μ = μ, a = λ /5 ; (b) for the same 2

21 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE sphere, but overed with ε 5ε same as (), but with μ =.25μ. =, μ = μ, a =.2a; () the same as (b), but with ε =.ε ; (d) the Figure Cotour plots of the distributio of the radial ompoet of the sattered mageti field i the y-z plae, aalogous to Fig.. Figure 2 Total time-averaged Poytig vetor distributios i the x-z plae for the ases of Figs. -. 2

22 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE Figure 3 Cotour plots of the distributio of the radial ompoets of the sattered eletri field i 22

23 Phys. Rev. E, 72, 6623 (9 pages), 25 DOI:.3/PhysRevE the x-z plae, idued by a plae wave travelig alog the z axis with a eletri field polarized alog the x axis: (a) for a sphere with ε = ε, μ = μ, a = λ /2 ; (b) for the same sphere, but overed with ε = 5ε, μ = μ, a =.9a; () the same as (b), but with ε =.ε. Figure 4 (a) The same as Fig. 7a; ad (b) the same as Fig. 8a, osiderig ohmi losses i the over materials. 23

After the completion of this section the student. V.4.2. Power Series Solution. V.4.3. The Method of Frobenius. V.4.4. Taylor Series Solution

After the completion of this section the student. V.4.2. Power Series Solution. V.4.3. The Method of Frobenius. V.4.4. Taylor Series Solution Chapter V ODE V.4 Power Series Solutio Otober, 8 385 V.4 Power Series Solutio Objetives: After the ompletio of this setio the studet - should reall the power series solutio of a liear ODE with variable

More information

Basic Waves and Optics

Basic Waves and Optics Lasers ad appliatios APPENDIX Basi Waves ad Optis. Eletromageti Waves The eletromageti wave osists of osillatig eletri ( E ) ad mageti ( B ) fields. The eletromageti spetrum is formed by the various possible

More information

Fluids Lecture 2 Notes

Fluids Lecture 2 Notes Fluids Leture Notes. Airfoil orte Sheet Models. Thi-Airfoil Aalysis Problem Readig: Aderso.,.7 Airfoil orte Sheet Models Surfae orte Sheet Model A aurate meas of represetig the flow about a airfoil i a

More information

4. Optical Resonators

4. Optical Resonators S. Blair September 3, 2003 47 4. Optial Resoators Optial resoators are used to build up large itesities with moderate iput. Iput Iteral Resoators are typially haraterized by their quality fator: Q w stored

More information

Observer Design with Reduced Measurement Information

Observer Design with Reduced Measurement Information Observer Desig with Redued Measuremet Iformatio I pratie all the states aot be measured so that SVF aot be used Istead oly a redued set of measuremets give by y = x + Du p is available where y( R We assume

More information

λ = 0.4 c 2nf max = n = 3orɛ R = 9

λ = 0.4 c 2nf max = n = 3orɛ R = 9 CHAPTER 14 14.1. A parallel-plate waveguide is kow to have a utoff wavelegth for the m 1 TE ad TM modes of λ 1 0.4 m. The guide is operated at wavelegth λ 1 mm. How may modes propagate? The utoff wavelegth

More information

ANOTHER PROOF FOR FERMAT S LAST THEOREM 1. INTRODUCTION

ANOTHER PROOF FOR FERMAT S LAST THEOREM 1. INTRODUCTION ANOTHER PROOF FOR FERMAT S LAST THEOREM Mugur B. RĂUŢ Correspodig author: Mugur B. RĂUŢ, E-mail: m_b_raut@yahoo.om Abstrat I this paper we propose aother proof for Fermat s Last Theorem (FLT). We foud

More information

Sx [ ] = x must yield a

Sx [ ] = x must yield a Math -b Leture #5 Notes This wee we start with a remider about oordiates of a vetor relative to a basis for a subspae ad the importat speial ase where the subspae is all of R. This freedom to desribe vetors

More information

Physics 3 (PHYF144) Chap 8: The Nature of Light and the Laws of Geometric Optics - 1

Physics 3 (PHYF144) Chap 8: The Nature of Light and the Laws of Geometric Optics - 1 Physis 3 (PHYF44) Chap 8: The Nature of Light ad the Laws of Geometri Optis - 8. The ature of light Before 0 th etury, there were two theories light was osidered to be a stream of partiles emitted by a

More information

Cloaked Sensors. Romain Fleury, Jason Soric and Andrea Alù * Dept. of Electrical and Computer Engineering, The University of Texas at Austin

Cloaked Sensors. Romain Fleury, Jason Soric and Andrea Alù * Dept. of Electrical and Computer Engineering, The University of Texas at Austin Physial Bouds o Absorptio ad Satterig for Cloaked Sesors Romai Fleury, Jaso Sori ad Adrea Alù * Dept. of Eletrial ad Computer Egieerig, The Uiversity of Texas at Austi 1 Uiversity Statio C83, Austi, TX

More information

= 47.5 ;! R. = 34.0 ; n air =

= 47.5 ;! R. = 34.0 ; n air = Setio 9: Refratio ad Total Iteral Refletio Tutorial Pratie, page 449 The agle of iidee is 65 The fat that the experimet takes plae i water does ot hage the agle of iidee Give:! i = 475 ;! R = 340 ; air

More information

Bernoulli Numbers. n(n+1) = n(n+1)(2n+1) = n(n 1) 2

Bernoulli Numbers. n(n+1) = n(n+1)(2n+1) = n(n 1) 2 Beroulli Numbers Beroulli umbers are amed after the great Swiss mathematiia Jaob Beroulli5-705 who used these umbers i the power-sum problem. The power-sum problem is to fid a formula for the sum of the

More information

Analog Filter Synthesis

Analog Filter Synthesis 6 Aalog Filter Sythesis Nam Pham Aubur Uiversity Bogda M. Wilamowsi Aubur Uiversity 6. Itrodutio...6-6. Methods to Sythesize Low-Pass Filter...6- Butterworth Low-Pass Filter Chebyshev Low-Pass Filter Iverse

More information

I. Existence of photon

I. Existence of photon I. Existee of photo MUX DEMUX 1 ight is a eletromageti wave of a high frequey. Maxwell s equatio H t E 0 E H 0 t E 0 H 0 1 E E E Aos( kzt ) t propagatig eletrial field while osillatig light frequey (Hz)

More information

Certain inclusion properties of subclass of starlike and convex functions of positive order involving Hohlov operator

Certain inclusion properties of subclass of starlike and convex functions of positive order involving Hohlov operator Iteratioal Joural of Pure ad Applied Mathematial Siees. ISSN 0972-9828 Volume 0, Number (207), pp. 85-97 Researh Idia Publiatios http://www.ripubliatio.om Certai ilusio properties of sublass of starlike

More information

Chapter 8 Hypothesis Testing

Chapter 8 Hypothesis Testing Chapter 8 for BST 695: Speial Topis i Statistial Theory Kui Zhag, Chapter 8 Hypothesis Testig Setio 8 Itrodutio Defiitio 8 A hypothesis is a statemet about a populatio parameter Defiitio 8 The two omplemetary

More information

Nonstandard Lorentz-Einstein transformations

Nonstandard Lorentz-Einstein transformations Nostadard Loretz-istei trasformatios Berhard Rothestei 1 ad Stefa Popesu 1) Politehia Uiversity of Timisoara, Physis Departmet, Timisoara, Romaia brothestei@gmail.om ) Siemes AG, rlage, Germay stefa.popesu@siemes.om

More information

Achieving transparency with plasmonic and metamaterial coatings

Achieving transparency with plasmonic and metamaterial coatings University of Pennsylvania ScholarlyCommons Departmental Papers (ESE) Department of Electrical & Systems Engineering 7-1-2005 Achieving transparency with plasmonic and metamaterial coatings Andrea Alù

More information

16th International Symposium on Ballistics San Francisco, CA, September 1996

16th International Symposium on Ballistics San Francisco, CA, September 1996 16th Iteratioal Symposium o Ballistis Sa Fraiso, CA, 3-8 September 1996 GURNEY FORULAS FOR EXPLOSIVE CHARGES SURROUNDING RIGID CORES William J. Flis, Dya East Corporatio, 36 Horizo Drive, Kig of Prussia,

More information

Epsilons Near Zero limits in the Mie scattering theory

Epsilons Near Zero limits in the Mie scattering theory Epsilos Near Zero limits i the Mie scatterig theory M. Tagviashvili * Adroikashvili Istitute of Physics, 6 Tamarashvili st, 77, Tbilsi, Georgia * madoa.tagviashvili@grt.ge Abstract: The classical Mie theory

More information

Chapter 4: Angle Modulation

Chapter 4: Angle Modulation 57 Chapter 4: Agle Modulatio 4.1 Itrodutio to Agle Modulatio This hapter desribes frequey odulatio (FM) ad phase odulatio (PM), whih are both fors of agle odulatio. Agle odulatio has several advatages

More information

Effect of Magnetic Field on Marangoni Convection in Relatively Hotter or Cooler Liquid Layer

Effect of Magnetic Field on Marangoni Convection in Relatively Hotter or Cooler Liquid Layer Iteratioal Joural of Advaed Researh i Physial Siee (IJARPS) Volume, Issue, Jauary 05, PP 7-3 ISSN 349-7874 (Prit) & ISSN 349-788 (Olie) www.arjourals.org ffet of Mageti Field o Maragoi Covetio i Relatively

More information

Summation Method for Some Special Series Exactly

Summation Method for Some Special Series Exactly The Iteratioal Joural of Mathematis, Siee, Tehology ad Maagemet (ISSN : 39-85) Vol. Issue Summatio Method for Some Speial Series Eatly D.A.Gismalla Deptt. Of Mathematis & omputer Studies Faulty of Siee

More information

Lecture 8. Dirac and Weierstrass

Lecture 8. Dirac and Weierstrass Leture 8. Dira ad Weierstrass Audrey Terras May 5, 9 A New Kid of Produt of Futios You are familiar with the poitwise produt of futios de ed by f g(x) f(x) g(x): You just tae the produt of the real umbers

More information

ME260W Mid-Term Exam Instructor: Xinyu Huang Date: Mar

ME260W Mid-Term Exam Instructor: Xinyu Huang Date: Mar ME60W Mid-Term Exam Istrutor: Xiyu Huag Date: Mar-03-005 Name: Grade: /00 Problem. A atilever beam is to be used as a sale. The bedig momet M at the gage loatio is P*L ad the strais o the top ad the bottom

More information

SYNTHESIS OF SIGNAL USING THE EXPONENTIAL FOURIER SERIES

SYNTHESIS OF SIGNAL USING THE EXPONENTIAL FOURIER SERIES SYNTHESIS OF SIGNAL USING THE EXPONENTIAL FOURIER SERIES Sadro Adriao Fasolo ad Luiao Leoel Medes Abstrat I 748, i Itrodutio i Aalysi Ifiitorum, Leohard Euler (707-783) stated the formula exp( jω = os(

More information

Principal Component Analysis

Principal Component Analysis Priipal Compoet Aalysis Nuo Vasoelos (Ke Kreutz-Delgado) UCSD Curse of dimesioality Typial observatio i Bayes deisio theory: Error ireases whe umber of features is large Eve for simple models (e.g. Gaussia)

More information

Digital Signal Processing. Homework 2 Solution. Due Monday 4 October Following the method on page 38, the difference equation

Digital Signal Processing. Homework 2 Solution. Due Monday 4 October Following the method on page 38, the difference equation Digital Sigal Proessig Homework Solutio Due Moda 4 Otober 00. Problem.4 Followig the method o page, the differee equatio [] (/4[-] + (/[-] x[-] has oeffiiets a0, a -/4, a /, ad b. For these oeffiiets A(z

More information

The Use of L-Moments in the Peak Over Threshold Approach for Estimating Extreme Quantiles of Wind Velocity

The Use of L-Moments in the Peak Over Threshold Approach for Estimating Extreme Quantiles of Wind Velocity The Use of L-Momets i the Pea Over Threshold Approah for Estimatig Extreme Quatiles of Wid Veloity M.D. Padey Uiversity of Waterloo, Otario, Caada P.H.A.J.M. va Gelder & J.K. Vrijlig Delft Uiversity of

More information

Principal Component Analysis. Nuno Vasconcelos ECE Department, UCSD

Principal Component Analysis. Nuno Vasconcelos ECE Department, UCSD Priipal Compoet Aalysis Nuo Vasoelos ECE Departmet, UCSD Curse of dimesioality typial observatio i Bayes deisio theory: error ireases whe umber of features is large problem: eve for simple models (e.g.

More information

THE MEASUREMENT OF THE SPEED OF THE LIGHT

THE MEASUREMENT OF THE SPEED OF THE LIGHT THE MEASUREMENT OF THE SPEED OF THE LIGHT Nyamjav, Dorjderem Abstrat The oe of the physis fudametal issues is a ature of the light. I this experimet we measured the speed of the light usig MihelsoÕs lassial

More information

Effects of Air Humidity on the Performance of a Polymer Insulator under Lightning Induced Voltage Conditions

Effects of Air Humidity on the Performance of a Polymer Insulator under Lightning Induced Voltage Conditions Effets of Air Humidity o the Performae of a Polymer Isulator uder Lightig Idued Voltage Coditios Mahdi Izadi *, Mohd Zaial Abidi Ab Kadir 2, Chadima Gomes 3, Mohd Syahmi 4, Maryam Hajihai 5,2,3,4,5 Cetre

More information

Class #25 Wednesday, April 19, 2018

Class #25 Wednesday, April 19, 2018 Cla # Wedesday, April 9, 8 PDE: More Heat Equatio with Derivative Boudary Coditios Let s do aother heat equatio problem similar to the previous oe. For this oe, I ll use a square plate (N = ), but I m

More information

Nonparametric Goodness-of-Fit Tests for Discrete, Grouped or Censored Data 1

Nonparametric Goodness-of-Fit Tests for Discrete, Grouped or Censored Data 1 Noparametri Goodess-of-Fit Tests for Disrete, Grouped or Cesored Data Boris Yu. Lemeshko, Ekateria V. Chimitova ad Stepa S. Kolesikov Novosibirsk State Tehial Uiversity Departmet of Applied Mathematis

More information

ε > 0 N N n N a n < ε. Now notice that a n = a n.

ε > 0 N N n N a n < ε. Now notice that a n = a n. 4 Sequees.5. Null sequees..5.. Defiitio. A ull sequee is a sequee (a ) N that overges to 0. Hee, by defiitio of (a ) N overges to 0, a sequee (a ) N is a ull sequee if ad oly if ( ) ε > 0 N N N a < ε..5..

More information

Lesson 8 Refraction of Light

Lesson 8 Refraction of Light Physis 30 Lesso 8 Refratio of Light Refer to Pearso pages 666 to 674. I. Refletio ad Refratio of Light At ay iterfae betwee two differet mediums, some light will be refleted ad some will be refrated, exept

More information

Physics 30 Lesson 8 Refraction of Light

Physics 30 Lesson 8 Refraction of Light Physis 30 Lesso 8 Refratio of Light Refer to Pearso pages 666 to 674. I. Refletio ad refratio of light At ay iterfae betwee two differet mediums, some light will be refleted ad some will be refrated, exept

More information

Lesson 10: Limits and Continuity

Lesson 10: Limits and Continuity www.scimsacademy.com Lesso 10: Limits ad Cotiuity SCIMS Academy 1 Limit of a fuctio The cocept of limit of a fuctio is cetral to all other cocepts i calculus (like cotiuity, derivative, defiite itegrals

More information

Construction of Control Chart for Random Queue Length for (M / M / c): ( / FCFS) Queueing Model Using Skewness

Construction of Control Chart for Random Queue Length for (M / M / c): ( / FCFS) Queueing Model Using Skewness Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5 Costrutio of Cotrol Chart for Radom Queue Legth for (M / M / ): ( / FCFS) Queueig Model Usig Skewess Dr.(Mrs.) A.R. Sudamai

More information

Metasurface Cloak Performance Near-by Multiple Line Sources and PEC Cylindrical Objects

Metasurface Cloak Performance Near-by Multiple Line Sources and PEC Cylindrical Objects Metaurfae Cloa Performae Near-by Multiple Lie Soure ad PEC Cylidrial Objet S. Arlaagić, W. Y. amilto, S. Pehro, ad A. B. Yaovlev 2 Departmet of Eletrial Egieerig Eletromageti Sytem Tehial Uiverity of Demar

More information

Chemistry 2. Assumed knowledge. Learning outcomes. The particle on a ring j = 3. Lecture 4. Cyclic π Systems

Chemistry 2. Assumed knowledge. Learning outcomes. The particle on a ring j = 3. Lecture 4. Cyclic π Systems Chemistry Leture QuatitativeMO Theoryfor Begiers: Cyli Systems Assumed kowledge Be able to predit the umber of eletros ad the presee of ougatio i a rig otaiig arbo ad/or heteroatoms suh as itroge ad oxyge.

More information

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c. 5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =

More information

2. SCHWARZSCHILD GEOMETRY REVISITED The metric for Schwarzschild Geometry is given by, ) (1) For constant values of time we have, c r

2. SCHWARZSCHILD GEOMETRY REVISITED The metric for Schwarzschild Geometry is given by, ) (1) For constant values of time we have, c r urved Spae-Tie ad the Speed of Light aitra Palit uthor/teaher, P-54 Motijheel veue, Motijheel Housig ooperative soiety, Flat- 4, Kolkata-700074, Idia, Eail: palit.aaitra@gail.o Keywords: Shwarzshild Geoetry,

More information

COMP26120: Introducing Complexity Analysis (2018/19) Lucas Cordeiro

COMP26120: Introducing Complexity Analysis (2018/19) Lucas Cordeiro COMP60: Itroduig Complexity Aalysis (08/9) Luas Cordeiro luas.ordeiro@mahester.a.uk Itroduig Complexity Aalysis Textbook: Algorithm Desig ad Appliatios, Goodrih, Mihael T. ad Roberto Tamassia (hapter )

More information

Quasi Normal Modes description. of transmission properties. for Photonic Band Gap structures.

Quasi Normal Modes description. of transmission properties. for Photonic Band Gap structures. Quasi Normal Modes desriptio of trasmissio properties for Photoi Bad Gap strutures. A. Settimi (1), S. Severii (), B. J. Hoeders (3) (1) INGV (Istituto Nazioale di Geofisia e Vulaologia) via di Viga Murata

More information

Principles of Communications Lecture 12: Noise in Modulation Systems. Chih-Wei Liu 劉志尉 National Chiao Tung University

Principles of Communications Lecture 12: Noise in Modulation Systems. Chih-Wei Liu 劉志尉 National Chiao Tung University Priiples of Commuiatios Leture 1: Noise i Modulatio Systems Chih-Wei Liu 劉志尉 Natioal Chiao ug Uiversity wliu@twis.ee.tu.edu.tw Outlies Sigal-to-Noise Ratio Noise ad Phase Errors i Coheret Systems Noise

More information

Statistics 511 Additional Materials

Statistics 511 Additional Materials Cofidece Itervals o mu Statistics 511 Additioal Materials This topic officially moves us from probability to statistics. We begi to discuss makig ifereces about the populatio. Oe way to differetiate probability

More information

Solutions 3.2-Page 215

Solutions 3.2-Page 215 Solutios.-Page Problem Fid the geeral solutios i powers of of the differetial equatios. State the reurree relatios ad the guarateed radius of overgee i eah ase. ) Substitutig,, ad ito the differetial equatio

More information

Quasi Normal Modes description of transmission properties for Photonic Band Gap structures.

Quasi Normal Modes description of transmission properties for Photonic Band Gap structures. Quasi ormal Modes desriptio of trasmissio properties for Photoi Bad Gap strutures. A. Settimi (1-), S. Severii (3), B. J. Hoeders (4) (1) FILAS (Fiaziaria Laziale di Sviluppo) via A. Farese 3, 19 Roma,

More information

High-Frequency Fluctuations of a Modulated, Helical Electron Beam

High-Frequency Fluctuations of a Modulated, Helical Electron Beam Publiatios 1996 High-Frequey Flutuatios of a Modulated, Helial Eletro Beam Mark Athoy Reyolds Embry-Riddle Aeroautial Uiversity, reyoldb@erau.edu Follow this ad additioal works at: https://ommos.erau.edu/publiatio

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

ME203 Section 4.1 Forced Vibration Response of Linear System Nov 4, 2002 (1) kx c x& m mg

ME203 Section 4.1 Forced Vibration Response of Linear System Nov 4, 2002 (1) kx c x& m mg ME3 Setio 4.1 Fored Vibratio Respose of Liear Syste Nov 4, Whe a liear ehaial syste is exited by a exteral fore, its respose will deped o the for of the exitatio fore F(t) ad the aout of dapig whih is

More information

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model.

Finally, we show how to determine the moments of an impulse response based on the example of the dispersion model. 5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).

More information

(8) 1f = f. can be viewed as a real vector space where addition is defined by ( a1+ bi

(8) 1f = f. can be viewed as a real vector space where addition is defined by ( a1+ bi Geeral Liear Spaes (Vetor Spaes) ad Solutios o ODEs Deiitio: A vetor spae V is a set, with additio ad salig o elemet deied or all elemets o the set, that is losed uder additio ad salig, otais a zero elemet

More information

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray

More information

SNAP Centre Workshop. Basic Algebraic Manipulation

SNAP Centre Workshop. Basic Algebraic Manipulation SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)

More information

POWER SERIES METHODS CHAPTER 8 SECTION 8.1 INTRODUCTION AND REVIEW OF POWER SERIES

POWER SERIES METHODS CHAPTER 8 SECTION 8.1 INTRODUCTION AND REVIEW OF POWER SERIES CHAPTER 8 POWER SERIES METHODS SECTION 8. INTRODUCTION AND REVIEW OF POWER SERIES The power series method osists of substitutig a series y = ito a give differetial equatio i order to determie what the

More information

Math 2784 (or 2794W) University of Connecticut

Math 2784 (or 2794W) University of Connecticut ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

Analysis of MOS Capacitor Loaded Annular Ring MICROSTRIP Antenna

Analysis of MOS Capacitor Loaded Annular Ring MICROSTRIP Antenna Iteratioal OPEN AESS Joural Of Moder Egieerig Research (IJMER Aalysis of MOS apacitor Loaded Aular Rig MIROSTRIP Atea Mohit Kumar, Suredra Kumar, Devedra Kumar 3, Ravi Kumar 4,, 3, 4 (Assistat Professor,

More information

Frequentist Inference

Frequentist Inference Frequetist Iferece The topics of the ext three sectios are useful applicatios of the Cetral Limit Theorem. Without kowig aythig about the uderlyig distributio of a sequece of radom variables {X i }, for

More information

Polynomials with Rational Roots that Differ by a Non-zero Constant. Generalities

Polynomials with Rational Roots that Differ by a Non-zero Constant. Generalities Polyomials with Ratioal Roots that Differ by a No-zero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by

More information

Ray-triangle intersection

Ray-triangle intersection Ray-triagle itersectio ria urless October 2006 I this hadout, we explore the steps eeded to compute the itersectio of a ray with a triagle, ad the to compute the barycetric coordiates of that itersectio.

More information

Monte Carlo Integration

Monte Carlo Integration Mote Carlo Itegratio I these otes we first review basic umerical itegratio methods (usig Riema approximatio ad the trapezoidal rule) ad their limitatios for evaluatig multidimesioal itegrals. Next we itroduce

More information

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka)

17 Phonons and conduction electrons in solids (Hiroshi Matsuoka) 7 Phoos ad coductio electros i solids Hiroshi Matsuoa I this chapter we will discuss a miimal microscopic model for phoos i a solid ad a miimal microscopic model for coductio electros i a simple metal.

More information

Economics 250 Assignment 1 Suggested Answers. 1. We have the following data set on the lengths (in minutes) of a sample of long-distance phone calls

Economics 250 Assignment 1 Suggested Answers. 1. We have the following data set on the lengths (in minutes) of a sample of long-distance phone calls Ecoomics 250 Assigmet 1 Suggested Aswers 1. We have the followig data set o the legths (i miutes) of a sample of log-distace phoe calls 1 20 10 20 13 23 3 7 18 7 4 5 15 7 29 10 18 10 10 23 4 12 8 6 (1)

More information

Lesson 4. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Lesson 4. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER) Lesso 4 Thermomehaial Measuremets for Eergy Systems (MENR) Measuremets for Mehaial Systems ad Produtio (MMER) A.Y. 15-16 Zaaria (Rio ) Del Prete RAPIDITY (Dyami Respose) So far the measurad (the physial

More information

B. Maddah ENMG 622 ENMG /20/09

B. Maddah ENMG 622 ENMG /20/09 B. Maddah ENMG 6 ENMG 5 5//9 Queueig Theory () Distributio of waitig time i M/M/ Let T q be the waitig time i queue of a ustomer. The it a be show that, ( ) t { q > } =. T t e Let T be the total time of

More information

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER / Statistics

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER / Statistics ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER 1 018/019 DR. ANTHONY BROWN 8. Statistics 8.1. Measures of Cetre: Mea, Media ad Mode. If we have a series of umbers the

More information

Topic 1 2: Sequences and Series. A sequence is an ordered list of numbers, e.g. 1, 2, 4, 8, 16, or

Topic 1 2: Sequences and Series. A sequence is an ordered list of numbers, e.g. 1, 2, 4, 8, 16, or Topic : Sequeces ad Series A sequece is a ordered list of umbers, e.g.,,, 8, 6, or,,,.... A series is a sum of the terms of a sequece, e.g. + + + 8 + 6 + or... Sigma Notatio b The otatio f ( k) is shorthad

More information

Recurrences: Methods and Examples

Recurrences: Methods and Examples Reurrees: Methods ad Examples CSE 30 Algorithms ad Data Strutures Alexadra Stefa Uiversity of exas at Arligto Updated: 308 Summatios Review Review slides o Summatios Reurrees Reursive algorithms It may

More information

Confidence intervals summary Conservative and approximate confidence intervals for a binomial p Examples. MATH1005 Statistics. Lecture 24. M.

Confidence intervals summary Conservative and approximate confidence intervals for a binomial p Examples. MATH1005 Statistics. Lecture 24. M. MATH1005 Statistics Lecture 24 M. Stewart School of Mathematics ad Statistics Uiversity of Sydey Outlie Cofidece itervals summary Coservative ad approximate cofidece itervals for a biomial p The aïve iterval

More information

(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?

(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m? MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle

More information

Estimation of a population proportion March 23,

Estimation of a population proportion March 23, 1 Social Studies 201 Notes for March 23, 2005 Estimatio of a populatio proportio Sectio 8.5, p. 521. For the most part, we have dealt with meas ad stadard deviatios this semester. This sectio of the otes

More information

RADICAL EXPRESSION. If a and x are real numbers and n is a positive integer, then x is an. n th root theorems: Example 1 Simplify

RADICAL EXPRESSION. If a and x are real numbers and n is a positive integer, then x is an. n th root theorems: Example 1 Simplify Example 1 Simplify 1.2A Radical Operatios a) 4 2 b) 16 1 2 c) 16 d) 2 e) 8 1 f) 8 What is the relatioship betwee a, b, c? What is the relatioship betwee d, e, f? If x = a, the x = = th root theorems: RADICAL

More information

Growing Evanescent Envelopes and Anomalous Tunneling in Cascaded Sets of Frequency-Selective Surfaces in Their Stop Bands

Growing Evanescent Envelopes and Anomalous Tunneling in Cascaded Sets of Frequency-Selective Surfaces in Their Stop Bands Growing Evanesent Envelopes and Anomalous Tunneling in Casaded Sets of Frequeny-Seletive Surfaes in Their Stop ands Andrea Alù Dept. of Applied Eletronis, University of Roma Tre, Rome, Italy. Nader Engheta

More information

Riemann Sums y = f (x)

Riemann Sums y = f (x) Riema Sums Recall that we have previously discussed the area problem I its simplest form we ca state it this way: The Area Problem Let f be a cotiuous, o-egative fuctio o the closed iterval [a, b] Fid

More information

Probability & Statistics Chapter 8

Probability & Statistics Chapter 8 I. Estimatig with Large Samples Probability & Statistis Poit Estimate of a parameter is a estimate of a populatio parameter give by a sigle umber. Use x (the sample mea) as a poit estimate for (the populatio

More information

Probability, Expectation Value and Uncertainty

Probability, Expectation Value and Uncertainty Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such

More information

Absorption and Emission of Radiation: Time Dependent Perturbation Theory Treatment

Absorption and Emission of Radiation: Time Dependent Perturbation Theory Treatment Absorptio ad Eissio of Radiatio: Tie Depedet Perturbatio Theory Treatet Wat Hailtoia for Charged Partile i E & M Field Need the potetial U. Fore o Charged Partile: 1 F e E V B Fore (geeralized for i Lagragia

More information

Dr R Tiwari, Associate Professor, Dept. of Mechanical Engg., IIT Guwahati,

Dr R Tiwari, Associate Professor, Dept. of Mechanical Engg., IIT Guwahati, Dr R Tiwari, Assoiate Professor, Dept. of Mehaial Egg., IIT Guwahati, (rtiwari@iitg.eret.i).3 Measuremet ad Sigal Proessig Whe we ivestigate the auses of vibratio, we first ivestigate the relatioship betwee

More information

PROPERTIES OF AN EULER SQUARE

PROPERTIES OF AN EULER SQUARE PROPERTIES OF N EULER SQURE bout 0 the mathematicia Leoard Euler discussed the properties a x array of letters or itegers ow kow as a Euler or Graeco-Lati Square Such squares have the property that every

More information

Michelson's Repetition of the Fizeau Experiment:

Michelson's Repetition of the Fizeau Experiment: Mihelso's Repetitio of the Fizeau Experimet: A Review of the Derivatio ad Cofirmatio of Fresel's Drag Coeffiiet A. A. Faraj a_a_faraj@hotmail.om Abstrat: I this ivestigatio, Mihelso's 1886 repetitio of

More information

The Growth of Functions. Theoretical Supplement

The Growth of Functions. Theoretical Supplement The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that

More information

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014 UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.-5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved

More information

LINEAR STABILITY ANALYSIS OF A PLANE-POISEUILLE HYDROMAGNETIC FLOW USING ADOMIAN DECOMPOSITION METHOD

LINEAR STABILITY ANALYSIS OF A PLANE-POISEUILLE HYDROMAGNETIC FLOW USING ADOMIAN DECOMPOSITION METHOD .P.B. Si. Bull., Series A, Vol. 75, Iss., 13 ISSN 13-77 LINEAR STABILITY ANALYSIS OF A PLANE-POISEILLE HYDROMAGNETIC FLOW SING ADOMIAN DECOMPOSITION METHOD Samuel O. ADESANYA 1 I this paper, the small-disturbaes

More information

Access to the published version may require journal subscription. Published with permission from: Elsevier.

Access to the published version may require journal subscription. Published with permission from: Elsevier. This is a author produced versio of a paper published i Statistics ad Probability Letters. This paper has bee peer-reviewed, it does ot iclude the joural pagiatio. Citatio for the published paper: Forkma,

More information

Infinite Sequences and Series

Infinite Sequences and Series Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet

More information

6.003 Homework #3 Solutions

6.003 Homework #3 Solutions 6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the

More information

Effect of Uniform Magnetic Field on Rayleigh Bènard Marangoni Convection in a Relatively Hotter or Cooler Layer of Liquid

Effect of Uniform Magnetic Field on Rayleigh Bènard Marangoni Convection in a Relatively Hotter or Cooler Layer of Liquid Iteratioal Joural of Advaed esearh i Physial Siee (IJAPS) Volume Issue 5 ay 5 PP 4-3 ISSN 349-7874 (Prit) & ISSN 349-788 (Olie) www.arjourals.org Effet of Uiform ageti Field o ayleigh Bèard aragoi Covetio

More information

Homework 6: Forced Vibrations Due Friday April 6, 2018

Homework 6: Forced Vibrations Due Friday April 6, 2018 EN40: Dyais ad Vibratios Hoework 6: Fored Vibratios Due Friday April 6, 018 Shool of Egieerig Brow Uiversity 1. The vibratio isolatio syste show i the figure has 0kg, k 19.8 kn / 1.59 kns / If the base

More information

x a x a Lecture 2 Series (See Chapter 1 in Boas)

x a x a Lecture 2 Series (See Chapter 1 in Boas) Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio

More information

1. (a) From Fig we find the smaller wavelength in question to be about 515 nm.

1. (a) From Fig we find the smaller wavelength in question to be about 515 nm. Chapter (a) From Fig - we fid the smaller wavelegth i questio to be about 55 m (b) Similarly, the larger wavelegth is approximately 6 m () From Fig - the wavelegth at whih the eye is most sesitive is about

More information

Chap.4 Ray Theory. The Ray theory equations. Plane wave of homogeneous medium

Chap.4 Ray Theory. The Ray theory equations. Plane wave of homogeneous medium The Ra theor equatio Plae wave of homogeeou medium Chap.4 Ra Theor A plae wave ha the dititive propert that it tregth ad diretio of propagatio do ot var a it propagate through a homogeeou medium p vae

More information

Anna Janicka Mathematical Statistics 2018/2019 Lecture 1, Parts 1 & 2

Anna Janicka Mathematical Statistics 2018/2019 Lecture 1, Parts 1 & 2 Aa Jaicka Mathematical Statistics 18/19 Lecture 1, Parts 1 & 1. Descriptive Statistics By the term descriptive statistics we will mea the tools used for quatitative descriptio of the properties of a sample

More information

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.

6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer. 6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio

More information

Section 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations

Section 1.1. Calculus: Areas And Tangents. Difference Equations to Differential Equations Differece Equatios to Differetial Equatios Sectio. Calculus: Areas Ad Tagets The study of calculus begis with questios about chage. What happes to the velocity of a swigig pedulum as its positio chages?

More information

exist for the work of spherical aggregate formation.

exist for the work of spherical aggregate formation. ISSN 06-933X Colloid oural 0 Vol 73 No 3 pp 406 47 Pleiades Publishig Ltd 0 Origial Russia Text AK Shhei S Kshevetsiy OS Pelevia 0 published i Kolloidyi Zhural 0 Vol 73 No 3 pp 404 46 iellizatio Kietis

More information

(Dependent or paired samples) Step (1): State the null and alternate hypotheses: Case1: One-tailed test (Right)

(Dependent or paired samples) Step (1): State the null and alternate hypotheses: Case1: One-tailed test (Right) (epedet or paired samples) Step (1): State the ull ad alterate hypotheses: Case1: Oe-tailed test (Right) Upper tail ritial (where u1> u or u1 -u> 0) H0: 0 H1: > 0 Case: Oe-tailed test (Left) Lower tail

More information

On Capacity and Delay of Multi-channel Wireless Networks with Infrastructure Support

On Capacity and Delay of Multi-channel Wireless Networks with Infrastructure Support O Capaity ad Delay of Multi-hael ireless Networks with Ifrastruture Support Hog-Nig Dai, Member, IEEE, Raymod Chi-ig og, Seior Member, IEEE, Hao ag, Member, IEEE Abstrat I this paper, we propose a ovel

More information