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

Size: px
Start display at page:

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

Transcription

1 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 65, 143 Roma, Italy () Cetro Iterforze Studi per le Appliazioi Militari via della Bigattiera 1, 561 Sa Piero a Grado (Pi), Italy (3) Researh group theory of odesed matter i the Istitute for theoretial physis ad Zerie Istitute for advaed materials Uiversity of Groige Nijeborg 4, NL 9747 AG Groige, Netherlads 1

2 Abstrat I this paper, we use the Quasi Normal Modes (QNM) approah for disussig the trasmissio properties of double-side opeed optial avities: i partiular, this approah is speified for oe dimesioal (1D) Photoi Bad Gap (PBG) strutures. Moreover, we ojeture that the desity of the modes (DOM) is a dyamial variable whih has the flexibility of varyig with respet to the boudary oditios as well as the iitial oditios; i fat, the e.m. field geerated by two moohromati outer-propagatig pump waves leads to iterferee effets iside a quarter-wave (QW) symmetri 1D-PBG struture. Fially, here, for the first time, a large umber of theoretial assumptios o QNM metris for a ope avity, ever disussed i literature, are proved, ad a simple ad diret method to alulate the QNM orm for a 1D-PBG struture is reported. OCIS umbers: 6.11 (Eletromageti optis) (Resoae) 1.7 (Trasmissio) (Photoi rystals)

3 1. Itrodutio. Photoi rystals a be viewed as partiular optial avities havig the properties of presetig allowed ad forbidde bads for the eletromageti radiatio travellig iside, at optial frequeies. For these motivatios these strutures are also amed Photoi Bad Gap (PBG) [1-4]. I these strutures, dispersive properties are usually evaluated assumig ifiite periodi oditios [5]. The fiite dimesios of PBG strutures oeptually modify the alulatio ad the ature of the dispersive properties: this is maily due to the existee of a eergy flow ito ad out of the rystal. A pheomeologial approah to the dispersive properties of oe-dimesioal (1D) PBG strutures has bee preseted [6]. The appliatio of the effetivemedium approah is disussed, ad the aalogy with a simple Fabry-Perot struture is developed by Sipe et al. [7]. The problem of the field desriptio iside a ope avity was disussed by several authors. I partiular, Leug et al. i s. [8-15] itrodued the desriptio of the eletromageti field i a oe side ope miroavities i terms of Quasi Normal Modes (QNMs). Miroavities are mesosopi ope systems. Sie they are ope, ad theore leay, i.e. ooservative, the resoae field eige-futios are QNMs with omplex eige-frequeies. They play importat roles i may optial proesses. The aalogy with ormal modes of oservative systems is emphasized. I partiular, i ertai ases, the QNMs form a omplete set, ad muh of the usual formalism a be arried through. Miroavities are also ope system, i whih there has to be some output ouplig. O aout of the output ouplig, the eletromageti field ad eergy i the miroavity aloe would be otiuously lost to the outside. Thus, i physial terms, the miroavity is a o oservative system, while, from a mathematial stadpoit, the operators whih appear would be o-hermitia (or ot self-adjoit). This the leads to iterestig halleges i attemptig to geeralize the familiar tools of quatum mehais ad mathematial physis to suh a o hermitia settig. The 3

4 issues whih arise, ad the framewor developed for addressig them, are geeri to may other situatios ivolvig ope systems. Theoretial results were obtaied i s. [16-19] whih developed a Hilbert-Shmidt type of theory, leadig to a biliear expasio i terms of the atural modes (QNMs) of the resolvet erel oeted with the itegral equatios of eletromageti ad potetial satterig theory (see erees therei [16-19]). The iitial value problem osidered by s. [16-19] ould be solved usig this Hilbert-Shmidt type of theory taig the temporal Laplae trasform of the Maxwell equatios. Appliative results were obtaied by Bertolotti [] who disussed the liear properties of oe dimesioal resoat avities by usig the matrix ad the ray method (see erees therei []). Fabry-Perot, photoi rystals ad 1D avities i PBG strutures are osidered. QNMs for the desriptio of the eletromageti field i ope avities are itrodued ad some appliatios are give. Theoretial ad appliative results were obtaied by. [1] whih developed a satterig theory for fiite 1D-PBG strutures i terms of the atural modes (QNMs) of the satterer. This theory geeralizes the lassial Hilbert-Shmidt type of biliear expasios of the field propagator to a biliear expasio ito atural modes (QNMs) (see erees therei [1]). It is show that the Sturm-Liouville type of expasios for dispersive media differs osiderably from those for odispersive media, they are e.g. overomplete. Colusive theoretial ad appliative results were obtaied by Masimovi et al. [] who used QNMs to haraterize trasmissio resoaes i 1D optial defet avities ad the related field approximatios. Ref. [] speializes to resoaes iside the badgap of the periodi multilayer mirrors whih elose the defet avities. Usig a field model with the most relevat QNMs, a variatioal priiple permits to represet the field ad the spetral trasmissio lose to resoaes. I s. [3,4], for the first time, the QNM approah was used ad exteded to the desriptio of the salar field behaviour i double-side opeed optial avities, i partiular 1D- 4

5 PBG strutures. The validity of the approah is disussed by provig the QNM ompleteess, disussig the omplex frequeies distributio, as well as the orrespodig field distributios, ad reoverig the behaviour of the desity of modes (DOM). I. [5], the eletromageti field iside a optial ope avity was aalyzed i the framewor of the QNM theory. The role of the omplex frequeies i the trasmissio oeffiiet ad their li with the DOM is larified. A appliatio to a quarter-wave (QW) symmetri 1D-PBG struture is disussed to illustrate the usefuless ad the meaig of the results. I. [6], by usig the QNM formalism i a seod quatizatio sheme, the problem of the outer-propagatio of eletromageti fields iside optial ope avities was studied. The lis betwee QNM operators ad aoial destrutio ad reatio operators desribig the exteral free field, as well as the field orrelatio futios, are foud ad disussed. A appliatio of the theory is performed for ope avities whose rative idex satisfies symmetri properties. I this paper, we use the QNM approah for disussig the trasmissio properties of doubleside opeed optial avities: i partiular, this approah is speified for 1D-PBG strutures. Moreover, we ojeture that the DOM is a dyamial variable whih has the flexibility of varyig with respet to the boudary oditios as well as the iitial oditios; i fat, the e.m. field geerated by two moohromati outer-propagatig pump waves leads to iterferee effets iside a quarter-wave (QW) symmetri 1D-PBG struture. Fially, here, for the first time, a large umber of theoretial assumptios o QNM metris for a ope avity, ever disussed i literature, are proved, ad a simple ad diret method to alulate the QNM orm for a 1D-PBG struture is reported. This paper is orgaized as follows. I setio, the QNM approah is itrodued. I setio 3, a large umber of theoretial assumptios o QNM metris for a ope avity, ever disussed i literature, are proved. I setio 4, the trasmissio oeffiiet for a ope avity is alulated as a superpositio of QNMs. I setio 5, the QNM approah is applied to a 1D-PBG struture ad a trasmissio oeffiiet formula is obtaied for a QW symmetri 1D-PBG struture. I setio 6, 5

6 the trasmissio resoaes are ompared with the QNM frequeies ad the trasmissio modes at the resoaes are alulated as super-positios of the QNM futios. I setio 7, it has bee show that the e.m. field geerated by two moohromati outer-propagatig pump waves leads to iterferee effets iside a QW symmetri 1D-PBG struture, suh that a physial iterpretatio of the DOM is proposed. I setio 8, a fial disussio is reported ad a ompariso is proposed about all the obtaied results ad their assoiated theoretial improvemets with respet to similar topis preseted i literature []. The Appedix desribes a simple ad diret method to alulate the QNM orm for a 1D-PBG struture.. Quasi Normal Mode (QNM) approah. With eree to fig. 1, osider a ope avity as a regio of legth L, filled with a material of a give rative idex ( x ), whih is elosed i a ifiite homogeeous exteral spae. The avity iludes also the termial surfaes, so it is represeted as C = [, L] ad the rest of uiverse as U = (,) ( L, ). The rative idex satisfies [8-15]: the disotiuity oditios, i.e. ( x ) presets a step at x = ad x = L, i this hypothesis a atural demaratio of a fiite regio is provided; the o tail oditios, i.e. ( x) = for x < ad x > L, i this hypothesis outgoig waves are ot sattered ba. The e.m. field E( x, t ) i the ope avity satisfies the equatio [7] x ρ( x) E( x, t) = t, (.1) where ρ ( x) = [ ( x) / ], beig the speed of light i vauum. If there is o exteral pumpig, the e.m. field satisfies suitable outgoig waves oditios [8-15][3,4] 6

7 E( x, t) = ρ E( x, t) for x <, (.) x t E( x, t) = ρ E( x, t) for x > L, (.3) x t where ( / ) ρ =, beig the outside rative idex. I fat, o the left side of the same avity, i.e. x <, the e.m. field is travellig i the egative sese of the x-axis, i.e. E( x, t) = E[ x + ( ) t], so eq. (.) holds as oe a easily prove. O the right side of the avity, i.e. x > L, the e.m. field is travellig i the positive sese of the x-axis, i.e. E( x, t) = E[ x ( ) t], so eq. (.3) holds as well. To tae the avity leaages ito aout, the Laplae trasform of the e.m. field is osidered [8], = Eɶ ( x, ) E( x, t) exp( it) dt, (.4) where is a omplex frequey. The e.m. field has to satisfy the Sommerfeld radiative oditio [7]: lim Eɶ ( x, ) =. (.5) x ± Sie s = i = exp( iπ ) with C [8], eq. (.4) defies a trasformatio whih loos lie a Fourier trasform with a omplex frequey, but it is a π -rotated Laplae trasform. The π -rotated Laplae trasform of the e. m. field overges to a aalytial futio Eɶ ( x, ) oly over the half-plae of overgee Im >. I fat, if the Laplae trasform (.4) is applied to the outgoig waves oditios (.), it follows ɶ(, ) = ρ ɶ (, ) for <, (.6) E x i E x x x ɶ = ρ ɶ >, (.7) E ( x, ) i E x ( x, ) for x L ad, solvig the last equatio (.7): Eɶ ( x, ) exp( i ρ x) = exp( i Re ρ x)exp( Im ρ x) for x > L. (.8) The Sommerfeld radiative oditio (.5) a be satisfied oly if Im >. 7

8 The trasformed Gree futio Gɶ ( x, x, ) a be defied by [7,9] x ɶ ; (.9) + ρ( x) G( x, x, ) = δ ( x x ) it is a e.m. field so, over the half-plae of overgee Im >, it satisfies the Sommerfeld radiative oditios [7,9]: exp( i ρ x) for x G( x, x, ) exp( i ρ x) for x ɶ. (.1) Two auxiliary futios g± ( x, ) a be defied by [8-15] ρ( x) g ( x, ) + ± = x ; (.11) they are ot defied as e.m. fields, beause, over the half-plae of overgee Im >, they satisfy oly the asymptoti oditios [8-15][3,4]: g+ ( x, ) exp( i ρ x) for x. (.1) g ( x, ) exp( i ρ x) for x However, the trasformed Gree futio Gɶ ( x, x, ) a be alulated i terms of the auxiliary futios g± ( x, ). I fat, it a be show that [8-15] the Wrosia W ( x, ) assoiated to the two auxiliary futios g± ( x, ) is x-idepedet, W ( x, ) = g ( x, ) g ( x, ) g ( x, ) g ( x, ) = W ( ), (.13) + + ad for the trasformed Gree futio [8-15]: g ( x, ) g ( x +, ) for x < x W ( ) Gɶ ( x, x, ) =. (.14) g+ ( x, ) g ( x, ) for x < x W ( ) 8

9 I what follows it is proved that, just beause of the asymptoti oditios (.1), the auxiliary futios g± ( x, ) are liearly idepedet over the half-plae of overgee Im >, ad so the Laplae trasformed Gree futio Gɶ ( x, x, ) is aalyti over Im >. The asymptoti oditios establish that, oly over Im >, the auxiliary futio g+ ( x, ) ats as a e.m. field for large x, beause it is expoetially deayig. I fat, from eq. (.1): g x i x i x x + (, ) = exp( ρ ) = exp( Re ρ ) exp( Im ρ ) for x. The, still over Im >, the other auxiliary futio g ( x, ) i geeral does ot at as a e.m. field for large x, so it is expoetially ireasig. I fat, aordig to eq. (.1): g ( x, ) A( ) exp( i ρ x) B( )exp( i ρ x) = + for x, with ( ) B, ad so g x B i ρ x B i ρ x ρ x (, ) ( )exp( ) = ( ) exp( Re )exp(im ), for x. It follows that the auxiliary futios g± ( x, ) are liearly idepedet over Im >, beause the Wrosia W ( ) is ot ull; i fat, from eq. (.13): W ( ) = lim[ g ( x, ) g ( x, ) g ( x, ) g ( x, )] = i ρ B( ) over Im >. Thus, the x + + trasformed Gree futio Gɶ ( x, x, ) is aalyti over Im >, where Gɶ ( x, x, ) does ot diverge; i fat, from eq. (.14): Gɶ ( x, x, ) 1/ W ( ), with W ( ) over Im >. For aalytial otiuatio [8], the trasformed Gree futio Gɶ ( x, x, ) a be exteded also over the lower omplex half-plae Im <. Aordig to. [3], it is always possible to defie a ifiite set of frequeies whih are the poles of the trasformed Gree futio Gɶ ( x, x, ), over the lower omplex half-plae Im <. I other words, there exists a ifiite set of omplex frequeies, = {, ± 1, ±, } whih the Wrosia (.13) is ull [8-15]: Z, with egative imagiary parts Im <, for W ( ) =. (.15) 9

10 The poles of the trasformed Gree futio are erred to as Quasi-Normal-Mode (QNM) eigefrequeies [3]. The defiitio of the QNM eige-frequeies implies that the auxiliary futios g± ( x, ) beome liearly depedet whe they are alulated at the QNM frequeies { }, Z =, ± 1, ±, ; so, the auxiliary futios i the QNM s are suh that [8-15] g ( x, ) = ( ) g ( x, ) = f ( x, ), (.16) + where ( ) is a suitable omplex ostat. The above futios f ( x, ) = f( x) are erred as Quasi-Normal-Mode (QNM) eige-futios [3]. The ouples [, f ( x)] are erred as Quasi- Normal-Modes beause: they are haraterized by omplex frequeies, so they are the ot-statioary modes of a ope avity [8-15] (they are observed i the frequey domai as resoaes of fiite width or i the time domai as damped osillatios); uder the disotiuity ad o tail oditios, the wave futios f ( x ) form a orthogoal basis oly iside the ope avity [8-15] (it is possible to desribe the QNMs i a maer parallel to the ormal modes of a losed avity). Applyig the QNM oditio (.16) to the equatio for the auxiliary futios (.11), it follows that the QNMs [, f ( x)] satisfy the equatio [8-15]: d dx + ( ) ( ) ρ x f x =. (.17) Moreover, applyig the QNM oditio (.16) to the asymptoti oditios for the auxiliary futios (.1), it follows that, oly for log distaes from the ope avity, the QNMs do ot represet e.m. fields beause they satisfy the QNM asymptoti oditios [8-15][3,4], f ( x) = exp( ± i ρ x) for x ±, (.18) 1

11 while, iside the same avity ad ear its termial surfaes from outside, the QNMs represet ot statioary modes, i fat the asymptoti oditios (.18) imply the formal QNM outgoig waves oditios: f ( x ) = i ρ f (), (.19) x x= f ( x) = i ρ f ( L). (.) x x= L The oditios (.19)-(.) are alled as formal beause they are erred to the QNMs whih do ot represet e.m. fields for log distaes from the ope avity, ad as outgoig waves beause they are formally idetial to the real outgoig waves oditios for the e.m. field i proximity to the surfaes of the avity. I fat, eqs. (.19)-(.) for the QNMs a be derived if ad oly if the requiremet of outgoig waves holds for the e.m. field. A ope system is ot oservative beause eergy a esape to the outside. As a result, the time-evolutio operator is ot Hermitia i the usual sese ad the eigefutios (fatorized solutios i spae ad time) are o loger ormal modes but quasi-ormal modes (QNMs) whose frequeies are omplex. QNM aalysis has bee a powerful tool for ivestigatig ope systems. Previous studies have bee mostly system speifi, ad use a few QNMs to provide approximate desriptios. I s. [8-15], the authors review developmets whih lead to a uifyig treatmet. The formulatio leads to a mathematial struture i lose aalogy to that i oservative, Hermitia systems. Hee may of the mathematial tools for the latter a be trasribed. Emphasis is plaed o those ases i whih the QNMs form a omplete set ad thus give a exat desriptio of the dyamis. More expliitly tha s. [8-15], we osider a Laplae trasform of the e.m. field, to tae the avity leaages ito aout, ad remar that, oly i the omplex domai defied by a π - rotated Laplae trasform, the QNMs a be defied as the poles of the trasformed Gree futio, with egative imagiary part. Our iter follows three steps: 11

12 1. the π -rotated Laplae trasform of the e. m. field [see eqs. (.4)-(.5)] overges to a aalytial futio Eɶ ( x, ) oly over the half-plae of overgee Im >.. just beause of the asymptoti oditios [see eq. (.1)], the auxiliary futios g± ( x, ) [see eq. (.11)] are liearly idepedet over the half-plae of overgee Im >, ad so the Laplae trasformed Gree futio Gɶ ( x, x, ) [see eq. (.9)] is aalyti over Im >. 3. with respet s. [8-15], we remar that: the trasformed Gree futio Gɶ ( x, x, ) [see eq. (.14)] a be exteded also over the lower omplex half-plae Im <, for aalytial otiuatio [8]; ad, it is always possible to defie a ifiite set of frequeies whih are the poles of the trasformed Gree futio Gɶ ( x, x, ) [see eqs. (.15)-(.16)], over the lower omplex half-plae Im <, aordig to [8-15]. 3. QNM metris. The QNM orm is defied as [8-15] f f dw =, (3.1) d = ad, with the method proposed for a oe side ope avity [8-15], oe a prove, for a both side ope avity, usig the QNM oditio (.18) ad the eqs. (.11)-(.1), that [8-15][3,4]: L = ρ + ρ + f f ( x) f ( x) dx i [ f () f ( L)]. (3.) Several remars about this geeralized orm are i order: it ivolves f ( x ) rather the f (x) ad it is i geeral omplex; it ivolves the two termial surfae terms If the QNM futio f ( x ) is ormalized, aordig to 1 i ρ f () ad i ρ f L. ( )

13 f N ( x) = f( x), (3.3) f f the: f N f N =. (3.4) The QNM ier produt is defied as [8-15][3,4] L N N N ˆ N ˆ N N N N N N m = m + m + ρ m + m + f f i f ( x) f ( x) f ( x) f ( x) dx i f () f () f ( L) f ( L), (3.5) if the QNM ojugate mometum f ˆ ( x ) is itrodued, aordig to: ˆ N N f ( x) = i ρ( x) f ( x). (3.6) Oe a prove that: the ier produt (3.5) is i aordae with the QNM orm (3.) [8-15]; the QNMs form a orthogoal basis iside the ope avity [8-14], i.e. [3,4] N N f fm = δ, m( + m), (3.7) beig δ, m the Kroeer delta, i.e. δ, m 1, =, = m m. Let us itrodue the overlappig itegral of the th QNM L N = ( ) ρ( ), (3.8) I f x x dx whih is lied to the statistial weight i the desity of modes (DOM) of the If the ope avity is haraterized by very slight leaages, th QNM [3,4]. Im << Re, (3.9) the the overlappig itegral of the th QNM overges to [3,4]: I 1. (3.1) 13

14 Proof: L eq.(3.3) N I = f ( x) ρ( x) dx = = f f L eq.(3.) f ( x) ρ( x) dx = L L f ( x) ρ( x) dx f ( x) ρ( x) dx = 1 L L ρ f ( ) ( ) [ () ( )] ( ) ( ) x ρ x dx + i f f x x dx + f L ρ ; (3.11) if Im << Re, the f ( x) f ( x), so eq. (3.1) holds. Istead, if the ope avity is haraterized by a some leaage, the the overlappig itegral of the th QNM a be alulated as [3,4]: ρ I f f L. (3.1) N N = () ( ) Im + Proof: L N N N ˆ N ˆ N N N N N N m m m ρ m m f f = i f ( x) f ( x) + f ( x) f ( x) dx + i f () f () + f ( L) f ( L) = L N N N N N N m m m m = ( + ) f ( x) ρ( x) f ( x) dx + i ρ f () f () + f ( L) f ( L) = = δ ( + ), m m eq.(3.7) eq.(3.6), (3.13) so L ρ f x x f x dx i f f f L f L N N N N N N ( ) ρ( ) ( ) = m δ, m [ () m () ( ) m ( )] + + m ; (3.14) if m =, the =, f ( x) = f ( x) ad more δ, =, so eq. (3.14) a be redued to eqs. (3.8) ad (3.1). Fially, if the ope avity is haraterized by very slight leaages [eq. (3.9)] ad its rative idex ( x ) satisfies the symmetry property ( L / x) = ( L / + x), (3.15) 14

15 the the QNM orm a be approximated i modulus: f f ρ. (3.16) Im Proof of eq. (3.16): ρ ρ eq.(3.1) I = f () f ( L) f () f ( L) + = + 1, (3.17) eq.(3.3) N N Im f f Im so f f f f L f L Hip: f () 1 = ρ () + ( ) = ρ 1 + ( ) Im Im ; (3.18) if the ope avity is symmetri [eq. (3.15) holds], the f ( L) = ( 1) f () = ( 1), so eq. (3.18) a be redued to eq. (3.16). Physial iterpretatio of eq. (3.16) : The modulus of QNM orm is expressed oly i terms of the QNM frequeies; The modulus of the QNM orm is high ( f f >> ρ ) whe the leaages of the ope avity are very slight ( Im << Re ); The QNM theory a be applied to ope avities ad is based o the outgoig waves oditios [eqs. (.)-(.3)] whih formalize some leaages ( Im < ), so the QNM theory a ot ilude the oservative ase, whe the avities are losed ad are ot haraterized by ay leaages ( Im = ). 4. Calulatio of trasmissio oeffiiet. With eree to fig. 1.b., let us ow osider a ope avity of legth L, filled with a rative idex ( x ), i the presee of a pump iomig from the left side. The avity iludes 15

16 also the surfaes, so it is represeted as C = [, L] ad the rest of uiverse as U = (,) ( L, ). The rative idex satisfies the disotiuity ad o tail oditios [8-15], as speified above. Uder these oditios, the QNMs form a omplete basis oly iside the avity, ad the e.m. field a be alulated as a superpositio of QNMs [8-15] N E( x, t) = a( t) f ( x), for x L, (4.1) N where f ( x ) are the ormalized QNM futios [see eq. (3.3)]. The superpositio oeffiiets a( t ) satisfy the dyami equatio [8-15] i aɺ ( t) + i a ( t) = f () b( t), (4.) N ρ where b( t ) is the drivig fore b( t) = ρ xep ( x, t). x= The left-pump EP ( x, t ) satisfies the iomig wave oditio [8-15]: b ( t ) = ρ E ( x, t ) = ρ E (, t ). (4.3) x P x= t P Eah QNM is drive by the drivig fore b( t ) ad at the same time deays beause of Im <. N The ouplig to the fore is determied by the surfae value of the QNM wave futio f (). The avity is i a steady state, so the Fourier trasform with a real frequey = Eɶ ( x, ) E( x, t)exp( it) dt a be applied to equatios (4.1)-(4.3), ad it follows N Eɶ ( x, ) = aɶ ( ) f ( x) N f () bɶ ( ) aɶ ( ) =, (4.4) ρ bɶ ( ) = iρeɶ P(, ) so N N f () f ( x) Eɶ ( x, ) = Eɶ P(, ) i ρ. (4.5) ( ) 16

17 With eree to fig. 1.b., the e.m. field is otiuous at avity surfaes x = ad x = L, so E ɶ + (, ) = E ɶ (, ) ad E ɶ( L, ) E ( + = ɶ L, ), ad the e.m. field Eɶ (, ) at surfae x = is the superpositio of the iomig pump E ɶ (, ) ad the leted field Eɶ (, ), so P R P E ɶ (, ) = E ɶ (, ) + E ɶ (, ), while the e.m. field Eɶ ( L, ) at the surfae x = L is oly the trasmitted field E ɶ ( L, ), so E ɶ( L, ) = E ɶ ( L, ). T T It follows that the trasmissio oeffiiet t( ) for a ope avity of legth L a be defied as the ratio betwee the trasmitted field Eɶ ( L, ) at the surfae x surfae x = [7]: = L ad the iomig pump Eɶ (, ) at E( L, ) t( ) = ɶ E (, ). (4.6) ɶ P The trasmissio oeffiiet is obtaied as superpositio of QNMs isertig eq. (4.5) i eq. (4.6): R P N N f () f ( L) t( ) = i ρ ( ). (4.7) Applyig the QNM ompleteess oditio [8-15] N N f ( x) f ( x ) = for x, x L, the trasmissio oeffiiet simplifies as: t( ) = i N N f () f ( L) ρ. (4.8) N Isertig the ormalized QNM futios f ( x ) [see eq.(3.3)], with f () = 1, equatio (4.8) beomes: f( L) t( ) = i ρ f f. (4.9) For a symmetri avity, suh that f ( L) = ( 1) f () = ( 1), fially [5] ( 1) t( ) =, (4.1) γ 17

18 where γ = f f / i ρ. I aordae with. [5], the trasmissio oeffiiet t( ) of a ope avity a be alulated as a superpositio of suitable futios, with QNM orms f f as weightig oeffiiets ad QNM frequeies as parameters. 5. Oe dimesioal Photoi Bad Gap (1D-PBG) strutures. With eree to fig., let us ow osider a symmetri oe dimesioal (1D) Photoi Bad Gap (PBG) struture whih osists of N periods plus oe layer; every period is omposed of two layers respetively with legths h ad l ad with rative idies h ad l, while the added layer is with parameters h ad h. The symmetri 1D-PBG struture osists of N + 1 layers with a total legth L = N( h + l) + h. If the two layers exteral to the symmetri 1D-PBG struture are osidered, the 1D-spae x a be divided ito N + 3 layers; they are L = [ x, x + 1], =,1,,N + 1, N +, with x =, x =, x = L, x = +. 1 N + N + 3 The rative idex ( x ) taes a ostat value i every layer L, =,1,, N + 1, N +, i.e. for x L, L N + ( x) = h for x L, = 1,3,,N 1, N + 1. (5.1) l for x L, =,4,,N For the symmetri 1D-PBG struture with rative idex (5.1), the QNM orm f f a be obtaied i terms of the frequey ad of the A ( ) N +, B ( ) N + oeffiiets for the g ( x, ) auxiliary futio i the L N + layer o the right side of the 1D-PBG struture (see Appedix): 18

19 db N + f f = i A N + ( ). (5.) d = As proved i. [3], for a quarter wave (QW) symmetri 1D-PBG struture with N periods ad as eree frequey, there are N + 1 families of QNMs, i.e. F QNM, [, N] ; the F family of QNMs osists of ifiite QNM frequeies, i.e. m = { ± ± } QNM,, Z, 1,,, m whih have the same imagiary part, i.e. Im, = Im,, Z, ad are lied with a step m m =, i.e. Re, m = Re, + m, m Z. It follows that, if the omplex plae is divided ito S = m < Re < ( m + 1), m Z, eah of the QNM family some rages, i.e. { } m QNM F drops oly oe QNM frequey over the rage S, i.e., = (Re, + m, Im,) ; there are N + 1 QNM m m frequeies over the rage S m ad they a be erred as, m =, + m, [, N ]. Thus, there N + QNM frequeies over the basi rage, i.e. S { } are 1, [, ], N. The QNM orm (5.) beomes: = < Re < ; they orrespod to, m db N + f, m f, m = i A N + (, m) d =, m. (5.3) The oeffiiets A ( ) N +, B ( ) N + are of the type exp( iδ ) where δ = π ( ) [3]. If oeffiiets A ( ) N +, B ( ) N + are alulated at QNM frequeies, m =, + m, where =, it follows A N + (, m) = A N + (,) ad ( db / d) = ( db / d), so N + =, m N + =, from eq. (5.3): γ f f γ = = ( + m ). (5.4), m, m,, m, i( / ), Thus, for the QNM F QNM family, the, m γ orm is periodi with a step Γ = γ, (,). The trasmissio oeffiiet (4.1) beomes 19

20 t( ) = ( 1) + m N + m,, (5.5) = m= γ, m, m ad, taig ito aout eq. (5.4), N, m 1 t( ) = ( 1) ( 1) γ m, (5.6) =, m=, it overges to: t π π. (5.7) N, ( ) = ( 1) s[ (,)] = γ, 6. QNMs ad trasmissio peas of a symmetri 1D-PBG struture. 6.a. QNM frequeies ad trasmissio resoaes. The trasmissio oeffiiet of a QW symmetri 1D-PBG struture, with N periods ad as eree frequey, is alulated as a superpositio of N + 1 QNM-oseats, whih orrespod to the umber of QNMs i [, ) rage. I fig. 3.a., the trasmissio oeffiiet predited by QNM theory ad the by umeri methods existig i the literature [31-33] are plotted for a QW symmetri 1D-PBG struture, where the eree wavelegth is λ = 1µ m, the umber of periods is N = 6 ad the values of rative idies are =, = 1.5. I figure 3.b., the two trasmissio oeffiiets are plotted for the same h l struture of fig. 3.a. but with = 3, =. I figure 3.., the two trasmissio oeffiiets are h l plotted for the same struture as i figure 3.b., with a ireased umber of periods N = 7. I aordae with. [5], we ote the good agreemet betwee the trasmissio oeffiiet predited by the QNM theory ad the oe obtaied by some umeri methods [31-33].

21 I the tables 1, the same examples of QW symmetri 1D-PBG strutures are osidered: (a) λ = 1µ m, N = 6, =, = 1,5, (b) λ = 1µ m, N = 6, = 3, =, () λ = 1µ m, N = 7, h l = 3, =. For eah of the three examples, the low ad high frequey bad-edges are h l desribed i terms of their resoaes ( / ) ad phases t(.. / ) ; besides, there is Bad Edge oe QNM lose to every sigle bad-edge: the real part of the QNM-frequey Re( / ) is reported together with the relative shift from the bad-edge resoae (Re..) /... h l B E QNM QNM B E B E From tables 1.a. ad 1.b. ( N = 6 ), the low ad high frequey bad-edges are haraterized by egative phases, whih sum is early ( π ), while, from table 1.. ( N = 7 ), their respetive phases are positive with a sum π approximately. I fat, a QW symmetri 1D-PBG struture, with N periods ad as eree frequey, presets a trasmissio spetrum with N+1 trasmissio pea pea resoaes i the [, ) rage, i.e. t( ) = exp( iϕ ), =,1,,N. Oly by umerial alulatios (MATLAB) o eq. (5.7), oe a sueed i provig that: ϕ ϕ + ϕ π = + 1 (N + 1) ( 1), 1,, N. (6.1) pea So, the lowest resoae i the rage [, ) has always a phase ϕ. The phases ϕ N ad ϕ N + 1 of the low ad high frequey bad-edges pea N ad pea N + 1 are suh that 1 1 ( 1) N + ϕn ϕn + π + ; so, for a QW symmetri 1D-PBG struture with N periods, the phases of the two bad-edges are suh that their sum is early π whe N is a odd umber or ( π ) approximately if N is a eve umber. Agai from tables 1, there is oe QNM lose to every sigle bad-edge, with a relative shift from the bad-edge [5]. I fat, a QW symmetri 1D-PBG struture, with N periods ad as eree frequey, presets a trasmissio oeffiiet whih is a superpositio of N + 1 QNMoseats (5.7), etred to the N + 1 QNM frequeies i [, ) rage, i.e. 1

22 QNM, =,1,,N. The QNM-oseats are ot sharp futios, so, whe they are superposed, a aliasig ours; the pea of the th QNM oseat a feel the tails of the ( 1) th ad ( + 1) th QNM-oseats, so the th resoae of the trasmissio spetrum results effetively shifted from the th QNM frequey i.e., pea = Re, =,1,,N. Comparig values i the tables 1.a. ( = =.5 ) ad table 1.b. ( = 1), the relative shift h l (Re ) / dereases with the rative idex step, while, omparig values i QNM B. E. B. E. tables 1.b. ( N = 6 ) ad 1.. ( N = 7 ), the same shift ireases with the umber of periods N. I aordae with. [5], some umerial simulatios prove that, the more the 1D-PBG struture presets a large umber of periods ( L >> λ ) with a high rative idex step ( >> ), the more the th QNM desribes the th trasmissio pea i the sese that Re, omes ear to pea. 6.b. QNM futios ad trasmissio modes. Let us osider oe moohromati pump Eɶ ( x, ), whih is iomig from the left side with uit amplitude ad zero iitial phase beig is the rative idex of the uiverse. p Eɶ (, ) ( ) exp( ) p x pɶ x = i x, (6.) The e.m. field Eɶ ( x, ), iside the QW symmetri 1D-PBG struture, a be alulated as a superpositio of QNMs [see setios 3 ad 4]:

23 Eɶ ( x, ) Eɶ ( x) = f () f ( x) f () f ( x) f ( x) f ( x) = pɶ () i ρ = i ρ = i ρ = = QNMs N N N N ( ) f f γ f ( x) + m f ( x) QW symmetri 1D-PBG struture N N, m,,, m = = = m= γ, m, m γ, m, m = = (6.3) If the pump (6.), iomig from the left side, is tued at the trasmissio pea whih is lose to the th QNM of the [, ) pea rage, i.e. Re, with =,1, N, the the trasmissio mode Eɶ ( x) a be approximated as a superpositio oly of the domiat terms N pea ɶ, m, m = m= Eɶ ( x) = a ( ) f ( x) = pea pea ɶ, m, m ɶ, m, m m= m= = a ( ) f ( x) + a ( ) f ( x), (6.4) m= aɶ ( ) f ( x) pea, m, m where / γ aɶ = = N + m Z = ± ±. (6.5) { } pea,,, m( ),,1, 1,, 1,, pea, m pea I fat, if the superpositio (6.4)-(6.5) is alulated at the frequey Re,, the eah pea domiat term is haraterized by a oeffiiet aɶ, ( ) with a deomiator iludig the m pea pea addedum, ( Re, ) i Im, i Im, =, so, ( pea pea aɶ m ) 1, 1 Im,, pea beig =,1, N + 1, while eah egleted term is haraterized by a oeffiiet aɶ, ( ) with pea pea a deomiator iludig the addedum, >>,, so m aɶ ( ) < 1 << 1 aɶ ( ), beig. pea pea pea pea, m,,, m Thus, the e.m. mode Eɶ ( x) iside the QW symmetri 1D-PBG struture, tued at the pea trasmissio pea lose to the th QNM of the [, ) 3 rage, beig =,1, N, a be approximated as the superpositio of the QNM futios f, ( x ) whih belog to the th QNM- m

24 family, beig m Z = {, ± 1, ±, } ; moreover, the weigh-oeffiiets of the superpositio pea pea aɶ, m( ) are alulated i the trasmissio resoae ad deped from the th QNMfamily. By umerial alulatios (MATLAB) oe o eqs. (6.4)-(6.5) oe a prove that, whe the pea trasmissio pea is loser to the QNM Re,, it is suffiiet to alulate the trasmissio mode (6.4) by the I order approximatio [] Eɶ ( x) aɶ ( ) f ( x), (6.6) pea,, where [] / γ aɶ N, (6.7) pea,,,( ) =,,1, 1 pea = +, otherwise, by the II order approximatio, more ied, pea Eɶ ( x) aɶ,( ) f,( x) +, (6.8) + aɶ ( ) f ( x) + aɶ ( ) f ( x) pea pea, 1, 1, + 1, + 1 where / γ aɶ N. (6.9) pea,,, 1( ) =,,1, 1 ± pea (, ± ) = + pea I aordae with. [], i the hypothesis that Re,, it follows pea pea pea, << (, ± ) < < (, + m ) <, so pea pea pea aɶ,( ) >> aɶ, ± 1 ( ) > > aɶ, ( ) >, m beig m : i fat, the e.m. mode (6.4)-(6.5) tued at the trasmissio pea pea, for the I order, a be approximated by eq. (6.6), orrespodig to the th QNM of the [, ) rage (the pea resoae is lose eough to the QNM frequey Re, ) []; the, for the II order, by eq. (6.8), iludig the two adjaet ( ± 1) th pea QNMs (if the resoae is far from the QNM frequeies Re =, 1 ±, but the ( ) pea ±, aɶ, futio tued at feels the tails of the aɶ ( ), ± 1 4

25 futios tued i Re, ± 1 ); fially, the ext orders a be egleted (beause the resoae is too far from the QNM frequeies Re, m =, + m, so the aɶ ( ), futio tued at does ot feel the tails of the aɶ, m( ) futios tued i Re, m, beig m ). pea pea By umerial alulatios (MATLAB) o eqs. (6.6) ad (6.8) oe a prove that, if the QNM is pea eve loser to the trasmissio pea of a QW symmetri 1D-PBG struture, i.e. Re,, whih is haraterized by a large umber of periods, i.e. L >> λ, ad a high rative idex step, i.e. >>, the the I order approximatio (6.6)-(6.7) a be redued to the QNM approximatio: Eɶ ( x) f ( x). (6.1), pea I fat, if Re,, the eq. (6.7) a be approximated as / γ / γ aɶ = = N + ; (6.11) pea,,,,,( ),,1, 1 pea ( Re,) i Im, i Im, if L >> λ, >>, the the 1D-PBG struture is haraterized by very slight leaages [see eq. (3.9)], ad, sie the rative idex satisfies suitable symmetry properties [see eq. (3.15)], the approximated QNM orm (3.16) a be used γ, = f, f,,,,1, N 1 ρ Im = + (6.1), ad, i modulus, the weigh-oeffiiet (6.11) overges to: 1 aɶ N. (6.13) pea,,( ) 1,,1, 1 γ, Im = +, Thus, i aordae with. [5], some umerial simulatios prove that the more the 1D-PBG struture presets a large umber of periods with a high rative idex step, the more the QNMs desribe the trasmissio peas i the sese that the QNM futios approximate the e.m. modes i the trasmissio resoaes. 5

26 I fig. 4, with eree to a QW symmetri 1D-PBG struture ( λ = 1µ m, N = 6, = 3, = ), the e.m. mode itesities (a) I I bad edge at the low frequey bad edge ( I bad edge / =.8 ) h l ad (b) I at the high frequey bad edge ( / = ), i uits of the II bad edge II bad edge itesity for a iomig pump I pump, are plotted as futios of the dimesioless spae x / L, where L is the legth of the 1D-PBG struture. It is lear the shiftig betwee the QNM approximatio (6.1) (the e.m. mode itesity I I bad edge is due to QNM lose to the bad-edge) ad the first order approximatio (6.6)-(6.7) [more ied tha (6.1), Ibad edge is alulated as the produt of the QNM, lose to the frequey bad-edge, for a weigh oeffiiet, whih taes ito aout the shift betwee the bad-edge resoae ad the QNM frequey] or the seod order approximatio (6.8)-(6.9) [eve more ied, but almost superimposed to (6.6)-(6.7), Ibad edge is alulated as the sum of the first order otributio due to the QNM, lose to the frequey badedge, ad the seod order otributios of the two adjaet QNMs, with the same imagiary part, so belogig to the same QNM family]. 7. Quarter wave (QW) symmetri 1D-PBG strutures exited by two outer-propagatig field pumps. Let us osider a QW symmetri 1D-PBG struture with N periods ad as eree frequey. If the uiverse iludes the termial surfaes of the 1D-PBG, the avity is represeted as C = [, L] ad the rest of uiverse as U = (,) ( L, ). The 1D-PBG presets a rative idex ( x ) whih satisfies the symmetry properties ( L / x) = ( L / + x). The avity is haraterized by a trasmissio spetrum with N + 1 resoaes i the [, ) rage [3,4], i.e. t( ) = exp[ iϕ ], r ( ) =, =,1, N [34]. 6

27 Let us osider oe pump, omig from the left side, whih is tued at the the 1D-PBG struture, with uit amplitude ad zero iitial phase th resoae of ( ) pɶ ( x) = exp( i x), (7.1) beig the rative idex of the uiverse. It exites i the avity a e.m. mode satisfies the equatio [7] ( E ɶ ) ( x) whih d Eɶ dx ad the boudary oditios are [7]: ( ) ( ) + ( x) E ( ) x = ɶ, (7.) ( ) ( ) ( ) ( ) ( ) E ɶ () = pɶ () + r ( ) pɶ () = pɶ () = 1. (7.3) ( ) ( ) Eɶ ( L) = t( ) pɶ () = t( ) Let us osider aother pump, omig from the right side, tued at the 1D-PBG struture with uit amplitude ad ostat phase-differee ϕ, i.e. th resoae of the pɶ ( x) = exp[ i ( )( x L)]exp[ i ϕ]. (7.4) ( ) If aloe, this pump exites i the avity a e.m. mode ( E ɶ ) ( x) whih satisfies a equatio similar to (7.) [7] d Eɶ dx ( ) ( ) + ( x) E ( ) x = ɶ, (7.5) beause of the symmetry properties of the rative idex, with boudary oditios whih are differet from (7.3) [7]: ( ) ( ) ( ) ( ) ( ) E ɶ ( L) = pɶ ( L) + r ( ) pɶ ( L) = pɶ ( L) = exp[ i ϕ]. (7.6) ( ) ( ) Eɶ () = t( ) pɶ ( L) = t( )exp[ i ϕ] It is easy to verify that the e.m. mode ( E ɶ ) ( x), whih is exited by the pump (7.4) ad solves the equatio (7.5) with the boudary oditios (7.6), is related to the e.m. mode ( E ɶ ) ( x), whih is exited by the pump (7.1) ad solves equatio (7.) with the boudary oditios (7.3), by the li: 7

28 Eɶ ( x) = [ Eɶ ( x)] t( )exp[ i ϕ]. (7.7) ( ) ( ) The two outer-propagatig pumps, whih are tued at the th resoae of the 1D-PBG struture, exite a e.m. mode Eɶ ( x) iside the avity, whih is the liear super-positio of ( E ɶ ) ( x) ad ( E ɶ ) ( x), i.e. [7] Eɶ x Eɶ x Eɶ x Eɶ x Eɶ x t i. (7.8) ( ) ( ) ( ) ( ) ( ) = ( ) + ( ) = ( ) + [ ( )] ( ) exp[ ϕ] It is easy to alulate from (7.8) the e.m. mode itesity iside the avity [7] Iɶ ( x) = Eɶ ( x)[ Eɶ ( x)] = ( ) = Iɶ ( x) + Re{[ Eɶ ( x)] t ( ) exp[ i ϕ]} = = Iɶ x + x + =, (7.9) ( ) ( ) ( ){1 os[ φ ( ) ϕ ϕ ]} ϕ ϕ = Iɶ x x + ( ) ( ) 4 ( )os [ φ ( ) ] if the e.m. mode exited by the pump (7.1) is represeted by Eɶ ( ) ( ) ( x I ) ( x) exp[ i ( = ɶ φ ) ( x)]. Eq. (7.9) shows that the e.m. field geerated by two moohromati outer-propagatig pump waves leads to iterferee effets iside a QW symmetri 1D-PBG struture Desity of modes (DOM). I fig. 5, with eree to a symmetri QW 1D-PBG struture ( λ = 1µ m, N = 6, = 3, l = ), the e.m. mode itesity exited iside the ope avity by oe field pump (- - -) omig from the left side [I order approximatio (6.6)] ad by two outer-propagatig field pumps i phase ( ) [see eq. (7.9)] are ompared whe eah of the two field pumps are tued at (a) the low frequey bad-edge ( I bad edge / =.8 ) or (b) the high frequey bad-edge ( II bad edge / = ). The e.m. mode itesities I I bad edge ad I II bad edge, i uits of the itesity for the two pumps I pump, are plotted as futios of the dimesioless spae x / L, where L is the legth of the 1D-PBG struture. If the two outer-propagatig field pumps i phase are 8 h

29 tued at the low frequey bad-edge, there is a ostrutive iterferee ad the field distributio i that bad-edge is remiiset of the distributio exited by oe field pump i the same trasmissio pea; while, if the two pumps are tued at the high frequey bad-edge, there is a destrutive iterferee ad almost o e.m. field peetratio ours i the struture. Some umerial simulatios o eq. (7.9) prove that, whe a QW symmetri 1D-PBG struture with N periods is exited by two outer-propagatig field pumps with a phase-differee ϕ, if N is a eve umber, the e.m. mode itesities i the low (high) frequey bad-edge ireases i stregth whe the two field pumps are i phase ϕ = (out of phase ϕ = π ) ad almost flags whe the two pumps are out phase ϕ = π (i phase ϕ = ). Moreover, if N is a odd umber, the two field distributios i the low ad high frequey bad-edges exhage their physial respose with respet to the phase-differee of the two pumps. The duality betwee the field distributios i the two bad-edges, for whih the oe ireases i stregth whe the other almost flags, a be explaied reallig that the trasmissio phases of the two bad edges are suh that their sum is π ( π ) whe N is a odd (eve) umber. Theore, if a 1D-PBG struture is exited by two outer-propagatig field pumps whih are tued at oe trasmissio resoae, o e.m. mode a be produed beause it depeds ot oly o the boudary oditios (the two field pumps) but also o the iitial oditios (the phasedifferee of the two pumps) [6]. We ojeture that the desity of modes is a dyami variable whih has the flexibility of varyig with respet to the boudary oditios as well as the iitial oditios, i aordae with. [6]. 8. A fial disussio. The Quasi Normal Mode (QNM) theory osiders the realisti situatio i whih a optial avity is ope from both sides ad is elosed i a ifiite homogeeous exteral spae; the la 9

30 of eergy oservatio for optial ope avities gives resoae field eige-futios with omplex eige-frequeies. The time-spae evolutio operator for the avity is ot hermitia: the modes of the field are quasi-ormal, i.e. they form a orthogoal basis oly iside the avity. The behavior of the eletromageti field i the optial domai, iside oe-dimesioal photoi rystals, is aalyzed by usig a extesio of the QNM theory. Oe-dimesioal (1D) Photoi Bad Gap (PBG) strutures are partiular optial avities, with both sides ope to the exteral eviromet, with a stratified material iside. These 1D-PBG strutures are fiite i spae ad, whe we wor with eletromageti pulses of a spatial extesio loger tha the legth of the avity, the 1D-PBG aot be studied as a ifiite struture: rather we have to osider the boudary oditios at the two eds of the avity. I the preset paper, the QNM approah has bee applied to a double-side opeed optial avity, ad speified for a 1D-PBG struture. All the osideratios made are exhaustively demostrated ad all the subsequetly results agree with the oes preseted i literature [3-5]. I this paper (se. ), suitable outgoig wave oditios for a ope avity without exteral pumpig are formalized [see eqs. (.)-(.3)]. More expliitly tha s. [8-15], here we have osidered a Laplae trasform of the e.m. field, to tae the avity leaages ito aout, ad we have remared that, oly i the omplex domai defied by a π -rotated Laplae trasform, the QNMs a be defied as the poles of the trasformed Gree futio, with egative imagiary part. Our iter has followed three steps: 1. the π -rotated Laplae trasform of the e. m. field [see eqs. (.4)-(.5)] overges to a aalytial futio Eɶ ( x, ) oly over the half-plae of overgee Im >.. just beause of the asymptoti oditios [see eq. (.1)], the auxiliary futios g± ( x, ) [see eq. (.11)] are liearly idepedet over the half-plae of overgee Im >, ad so the Laplae trasformed Gree futio Gɶ ( x, x, ) [see eq. (.9)] is aalyti over Im >. 3

31 3. Differetly from s. [8-15], here we have remared that: the trasformed Gree futio Gɶ ( x, x, ) [see eq. (.14)] a be exteded also over the lower omplex half-plae Im <, for aalytial otiuatio [8]; ad, it is always possible to defie a ifiite set of frequeies whih are the poles of the trasformed Gree futio Gɶ ( x, x, ) [see eqs. (.15)-(.16)], over the lower omplex half-plae Im <, aordig to. [3]. Moreover (se. 3), if the ope avity is haraterized by very slight leaages ad its rative idex satisfies suitable symmetry properties, more tha s. [8-15], here we have proved that the modulus of the QNM orm i modulus a be expressed oly i terms of the QNM frequeies [see eq. (3.16)] ad, more tha s. [3-5], here we have proposed its physial iterpretatio: the QNM orm is high whe the leaages of the ope avity are very slight; the QNM theory a be applied to ope avities ad is based o the outgoig waves oditios whih formalize some leaages ( Im < ), so the QNM theory a ot ilude the oservative ase, whe the avities are losed ad are ot haraterized by ay leaages ( Im = ). The paper (se. 4) has give a proof that the trasmissio oeffiiet t( ) a be alulated as a superpositio of suitable futios, with QNM orms f f as weightig oeffiiets ad QNM frequeies as parameters [see eq.(4.1)], thesis oly postulated i. [5]. The QNM approah has bee applied to a 1D-PBG struture (se. 5): more tha. [5], here the trasmissio oeffiiet of a quarter wave (QW) symmetri 1D-PBG struture, with N periods ad as eree frequey, has bee alulated as a superpositio of N + 1 QNM-oseats, whih orrespod to the umber of QNMs i [, ) rage [see eq. (5.7)]; differetly from the assertio of. [1] for whih the frequey depedee of t( ) ould suggest at first glae that the positios of maxima of T ( ) = t( ) are realized whe the real frequey oiides with the real part of the QNM frequeies, here we have dedued that this fat is true oly i the further hypothesis of real parts of QNM frequeies well separated o the real frequey axis, whereas, 31

32 o the otrary, overlappig of the tails of two ear QNM-oseat otributios to T ( ) lead to a sigifiatly displayed positio of the maxima of T ( ), a sort of aliasig proess (see figs. 4.a-. ad tab. 1.a-.). I se. 6., the trasmissio resoaes have bee ompared with the QNM frequeies (se. 6.a.) ad the trasmissio modes i the resoaes have bee alulated as super-positios of the QNM futios (se. 6.b.): differetly from. [], here we affirm expliitly that a proper field represetatio for the trasmittae problem a be established i terms of QNMs, ad the boudary oditios are ot violated [see eqs. (6.)-(6.3)]; eve if the fields are etirely differet whe osidered o the whole real axis, the e.m. field a be represeted as a superpositio of QNMs iside the ope avity ad o the limitig surfaes too. I the preset paper (se. 7), it has bee show that the e.m. field geerated by two moohromati outer-propagatig pump waves leads to iterferee effets iside a QW symmetri 1D-PBG struture [see eq. (7.9)]. Eve if this paper is hroologially subsequet to. [5], ayway it is oeived as a trait d uio betwee s. [5] ad [6], i fat the paper forestalls oeptually, aordig to lassi eletrodyamis, what. [6] develops i terms of quatum mehais. We have show that if a 1D-PBG struture is exited by two outer-propagatig pumps whih are tued at oe trasmissio resoae, o e.m. mode a be produed beause it depeds ot oly o the boudary oditios (the two pumps) but also o the iitial oditios (the phase-differee of the two pumps) [see fig. 5.a-b.]; i aordae with. [6], we have ojetured (se 7.1) that the desity of the modes (DOM) is a dyamial variable whih has the flexibility to adjust with respet to the boudary oditios as well the iitial oditios. I Appedix, a diret method has bee desribed to obtai the QNM orm for a QW symmetri 1D-PBG struture: the QNM orm f f a be obtaied i terms of the 3

33 frequey ad of the A ( ) N +, B ( ) N + oeffiiets for the g ( x, ) auxiliary futio i the L N + layer o the right side of the 1D-PBG struture A ompariso with Masimovi et al.. []. Let us propose a ritial ompariso betwee Masimovi ad our methods: Masimovi speializes. [] to fiite periodi strutures whih possess trasmissio properties with a badgap, i.e. with a regio of frequeies of very low trasmissio. Ref. [] hooses a field model for the trasmittae problem ad taes the relevat QNMs as those with the real part of their omplex frequey i the give frequey rage. To determie the deompositio oeffiiets i the field model, a variatioal form of the trasmittae problem is employed. The trasmittae problem orrespods to the equatio ad atural boudary oditios, arisig form the oditio of statioarity of a futioal F[ E( x, t )]. If F[ E( x, t )] beomes statioary, i.e. if the first variatio of F[ E( x, t )] vaishes for arbitrary variatios of the e.m. field E( x, t ), the E( x, t ) satisfies the atural boudary oditios. The optimal oeffiiets a be obtaied as solutios of a system of liear equatios. By solvig the system for eah value of the real frequey, the deompositio oeffiiets i the field represetatio for trasmittae problem are obtaied. This eables approximatio of the spetral trasmittae ad letae ad the related field profile. Masimovi applies a method whih is based o a mathematial priiple, ivolvig arduous alulatios. Istead, we have proposed, for a ope avity, a physial method whih starts from the boudary oditios at the surfaes of the avity (se. 4). Moreover, we have provided, for a symmetri 1D-PBG struture (se. 5), a ew algorithm to alulate the QNM orm by a derivative [see eq. (5.)], more diret tha the itegral of. []; ad, we have proved that a QW symmetri 1D-PBG struture with N periods is haraterized by 33

34 N + 1 families of QNM frequeies, periodially distributed o the omplex plae, suh that, for eah family, the QNM orm is periodial with a proper step [see eq. (5.4)]: as a result, the alulatio of the trasmissio oeffiiet is eve more simplified with respet to. []. A ommo assumptio made i the literature [8-15][][3,4][35] is that the spetral trasmissio for sigle resoae situatio, as desribed, is of a Loretia lie shape. We have falle ito lie of this assumptio (se. 6.b): i fat, the e.m. field, iside a QW symmetri 1D-PBG struture, has bee alulated as a superpositio of QNMs; so, eah sigle resoae is desribed by a Loretia lie shape [see eqs. (6.3), (6.4)-(6.5) ad ommets below]. The method of. [] justifies aalytially this ommo assumptio. Masimovi osiders the otributio of a sigle QNM i the field model. Here (se. 6.b.), we have osidered the otributios of the adjaet QNMs, by more ad more ied approximatios iludig suessive orders: i fat, some umerial simulatios have proved that (see figs. 4.a-b.), whe the trasmissio pea is lose to the QNM, it is suffiiet to alulate the trasmissio mode by the I order approximatio [see eqs. (6.6)-(6.7)], otherwise, by the II order approximatio, more ied [see eqs. (6.8)-(6.9) ad ommets below]; so, i our opiio, the otributio of a sigle QNM is suffiiet for the field model oly i the limit of arrow resoaes [see eqs. (6.1), (6.11)-(6.13) ad ommets below]. Appedix: QNM orm for 1D-PBG strutures. With eree to fig., let us ow osider a symmetri oe dimesioal photoi bad gap (1D-PBG) struture whih osists of N period plus oe layer; every period is omposed of two 34

35 layers respetively with legths h ad l ad with rative idies h ad l, while the added layer is with parameters h ad h. The symmetri 1D-PBG struture osists of N + 1 layers with a total legth L = N( h + l) + h. If the two layers exteral to the symmetri 1D-PBG struture are osidered, the 1D-spae x a be divided ito N + 3 layers; the two exteral layers, L = (, ) ad L N + = ( L, + ), ad the 1D- PBG layers, the odd oes L = 1 [( 1)( h + l), ( 1)( h + l) + h], = 1,,, N + 1 ad the eve oes L = [( 1)( h + l) + h, ( h + l)], = 1,,, N. The rative idex ( x ) taes a ostat value i every layer L, =,1,, N + 1, N +, i.e. for x L, L N + ( x) = h for x L, = 1,3,,N 1, N + 1. (A.1) l for x L, =,4,,N At first, the alulatio of quasi ormal modes (QNMs) is stated for the 1D-PBG struture with rative idex (A.1). The homogeeous equatio is solved for the auxiliary futios (.11), with the "asymptoti oditios" (.1). The auxiliary futio g ( x, ) is i x i x g ( x, ) = Ae + Be ϑ( x) + N + 1 ih x ih x + A 1e + B 1e ϑ [ x ( 1)( h + l)] ϑ[( 1)( h + l) + h x] + = 1, (A.) N il x il x + A e + Be ϑ [ x ( 1)( h + l) h] ϑ[ ( h + l) x] + = 1 i x i x + A N + e + B N + e ϑ( x L) where ϑ ( x) is the uit step futio; while the auxiliary futio g+ ( x, ) has a similar expressio to (A.), where some oeffiiets C ( ) ad D ( ) replae A ( ) ad B ( ) with =,1,, N + 1, N +. The auxiliary futio g ( x, ) is a right to left wave for x <, so 35

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

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

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

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

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

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

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

Quasi-normal-modes description of transmission properties for photonic bandgap structures Settimi, A.; Severini, S.; Hoenders, B. J.

Quasi-normal-modes description of transmission properties for photonic bandgap structures Settimi, A.; Severini, S.; Hoenders, B. J. Uiversity of Groige Quasi-ormal-modes descriptio of trasmissio properties for photoic badgap structures Settimi, A.; Severii, S.; Hoeders, B. J. Published i: Joural of the Optical Society of America B-Optical

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

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

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

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

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

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

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

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

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

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

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

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

(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

The beta density, Bayes, Laplace, and Pólya

The beta density, Bayes, Laplace, and Pólya The beta desity, Bayes, Laplae, ad Pólya Saad Meimeh The beta desity as a ojugate form Suppose that is a biomial radom variable with idex ad parameter p, i.e. ( ) P ( p) p ( p) Applyig Bayes s rule, we

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

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

= 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

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

λ = 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

ε > 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

MULTILEVEL ANALYSIS OF DELAMINATION INITIATED NEAR THE EDGES OF COMPOSITE STRUCTURES

MULTILEVEL ANALYSIS OF DELAMINATION INITIATED NEAR THE EDGES OF COMPOSITE STRUCTURES MULTILEVEL ANALYSIS OF DELAMINATION INITIATED NEAR THE EDGES OF COMPOSITE STRUCTURES N. Carrere 1, T. Vadellos 1, E. Marti 1 ONERA, 9 av. de la Divisio Leler, 930 Châtillo, Frae LCTS, 3 Allée de la Boétie,

More information

Appendix: The Laplace Transform

Appendix: The Laplace Transform Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio

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

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

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

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

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

Explicit and closed formed solution of a differential equation. Closed form: since finite algebraic combination of. converges for x x0

Explicit and closed formed solution of a differential equation. Closed form: since finite algebraic combination of. converges for x x0 Chapter 4 Series Solutios Epliit ad losed formed solutio of a differetial equatio y' y ; y() 3 ( ) ( 5 e ) y Closed form: sie fiite algebrai ombiatio of elemetary futios Series solutio: givig y ( ) as

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

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001. Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;

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

Chapter 7: The z-transform. Chih-Wei Liu

Chapter 7: The z-transform. Chih-Wei Liu Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability

More information

END CONDITIONS OF PIANO STRINGS Palaiseau Cedex,

END CONDITIONS OF PIANO STRINGS Palaiseau Cedex, END CONDITIONS OF PIANO STRINGS Kerem Ege, Atoie Chaige aboratory for Solid Mehais, Eole Polytehique, UMR7649, 98 Palaiseau Cedex, erem.ege@lms.polytehique.fr Uité de Méaique, Eole Natioale Supérieure

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

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

Lecture 1: Semiconductor Physics I. Fermi surface of a cubic semiconductor

Lecture 1: Semiconductor Physics I. Fermi surface of a cubic semiconductor Leture 1: Semiodutor Physis I Fermi surfae of a ubi semiodutor 1 Leture 1: Semiodutor Physis I Cotet: Eergy bads, Fermi-Dira distributio, Desity of States, Dopig Readig guide: 1.1 1.5 Ludstrom 3D Eergy

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

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

Thermodynamics of the Primary Eigen Gas and the Postulates of Quantum Mechanics

Thermodynamics of the Primary Eigen Gas and the Postulates of Quantum Mechanics Thermodyamis of the Primary Eige Gas ad the Postulates of Quatum Mehais V.A.I. Meo, Gujarat Uiversity Campus, Ahmedabad-380009, Idia. Abstrat The author shows that that for eah quatum mehaial property

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

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

Matsubara-Green s Functions

Matsubara-Green s Functions Matsubara-Gree s Fuctios Time Orderig : Cosider the followig operator If H = H the we ca trivially factorise this as, E(s = e s(h+ E(s = e sh e s I geeral this is ot true. However for practical applicatio

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

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

Office: JILA A709; Phone ;

Office: JILA A709; Phone ; Office: JILA A709; Phoe 303-49-7841; email: weberjm@jila.colorado.edu Problem Set 5 To be retured before the ed of class o Wedesday, September 3, 015 (give to me i perso or slide uder office door). 1.

More information

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc.

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc. Quatum Mechaics I 1 April, 14 Assigmet 5: Solutio 1 For a particle icidet o a potetial step with E < V, show that the magitudes of the amplitudes of the icidet ad reflected waves fuctios are the same Fid

More information

Mass Transfer Chapter 3. Diffusion in Concentrated Solutions

Mass Transfer Chapter 3. Diffusion in Concentrated Solutions Mass Trasfer Chapter 3 Diffusio i Coetrated Solutios. Otober 07 3. DIFFUSION IN CONCENTRATED SOLUTIONS 3. Theor Diffusio auses ovetio i fluids Covetive flow ours beause of pressure gradiets (most ommo)

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

Société de Calcul Mathématique SA Mathematical Modelling Company, Corp.

Société de Calcul Mathématique SA Mathematical Modelling Company, Corp. oiété de Calul Mathéatique A Matheatial Modellig Copay, Corp. Deisio-aig tools, sie 995 iple Rado Wals Part V Khihi's Law of the Iterated Logarith: Quatitative versios by Berard Beauzay August 8 I this

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

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

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

The Use of Filters in Topology

The Use of Filters in Topology The Use of Filters i Topology By ABDELLATF DASSER B.S. Uiversity of Cetral Florida, 2002 A thesis submitted i partial fulfillmet of the requiremets for the degree of Master of Siee i the Departmet of Mathematis

More information

arxiv:hep-th/ Mar 97

arxiv:hep-th/ Mar 97 Priiples of Disrete ime Mehais: II. Classial Field heory George Jaroszkiewiz ad Keith Norto Departmet of Mathematis, Uiversity of Nottigham Uiversity Park, Nottigham NG7 2RD, UK 19 th Deember 199 arxiv:hep-th/9703080

More information

Basic Probability/Statistical Theory I

Basic Probability/Statistical Theory I Basi Probability/Statistial Theory I Epetatio The epetatio or epeted values of a disrete radom variable X is the arithmeti mea of the radom variable s distributio. E[ X ] p( X ) all Epetatio by oditioig

More information

Local Estimates for the Koornwinder Jacobi-Type Polynomials

Local Estimates for the Koornwinder Jacobi-Type Polynomials Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 93-9466 Vol. 6 Issue (Jue 0) pp. 6 70 (reviously Vol. 6 Issue pp. 90 90) Appliatios ad Applied Mathematis: A Iteratioal Joural (AAM) Loal Estimates

More information

Fuzzy Dynamic Characteristic of Concrete. Material under Impact Loads

Fuzzy Dynamic Characteristic of Concrete. Material under Impact Loads Proeedigs of the 2d WSEAS It. Coferee o Applied ad Theoretial Mehais, Veie, Italy, November 2-22, 26 222 Fuzzy Dyami Charateristi of Corete Material uder Impat Loa GAO SHIQIAO LIU HAIPENG JIN LEI Shool

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

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

Lesson 03 Heat Equation with Different BCs

Lesson 03 Heat Equation with Different BCs PDE & Complex Variables P3- esso 3 Heat Equatio with Differet BCs ( ) Physical meaig (SJF ) et u(x, represet the temperature of a thi rod govered by the (coductio) heat equatio: u t =α u xx (3.) where

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

A first-order phase-transition, a super-cooled fluid, and a glass in a twodimensional. and A. Baram (2)

A first-order phase-transition, a super-cooled fluid, and a glass in a twodimensional. and A. Baram (2) A first-order phase-trasitio, a super-ooled fluid, ad a glass i a twodimesioal lattie gas model E. Eiseberg ( ) ad A. Baram (2). Shool of Physis ad Astroomy, Raymod ad Beverly Sakler Faulty of Exat Siees,

More information

Analysis of spatial temporal converters for all-optical communication links

Analysis of spatial temporal converters for all-optical communication links Aalysis of spatial temporal overters for all-optial ommuiatio liks Da M. Marom, Pag-Che Su, ad Yeshaiahu Faima We aalyze parallel-to-serial trasmitters ad serial-to-parallel reeivers that use ultrashort

More information

Math 213b (Spring 2005) Yum-Tong Siu 1. Explicit Formula for Logarithmic Derivative of Riemann Zeta Function

Math 213b (Spring 2005) Yum-Tong Siu 1. Explicit Formula for Logarithmic Derivative of Riemann Zeta Function Math 3b Sprig 005 Yum-og Siu Expliit Formula for Logarithmi Derivative of Riema Zeta Futio he expliit formula for the logarithmi derivative of the Riema zeta futio i the appliatio to it of the Perro formula

More information

First, note that the LS residuals are orthogonal to the regressors. X Xb X y = 0 ( normal equations ; (k 1) ) So,

First, note that the LS residuals are orthogonal to the regressors. X Xb X y = 0 ( normal equations ; (k 1) ) So, 0 2. OLS Part II The OLS residuals are orthogoal to the regressors. If the model icludes a itercept, the orthogoality of the residuals ad regressors gives rise to three results, which have limited practical

More information

The Riemann Zeta Function

The Riemann Zeta Function Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we

More information

REVERSIBLE NON-FLOW PROCESS CONSTANT VOLUME PROCESS (ISOCHORIC PROCESS) In a constant volume process, he working substance is contained in a rigid

REVERSIBLE NON-FLOW PROCESS CONSTANT VOLUME PROCESS (ISOCHORIC PROCESS) In a constant volume process, he working substance is contained in a rigid REVERSIBLE NON-FLOW PROCESS CONSTANT VOLUME PROCESS (ISOCHORIC PROCESS) I a ostat olume roess, he workig substae is otaied i a rigid essel, hee the boudaries of the system are immoable, so work aot be

More information

Abstract. Fermat's Last Theorem Proved on a Single Page. "The simplest solution is usually the best solution"---albert Einstein

Abstract. Fermat's Last Theorem Proved on a Single Page. The simplest solution is usually the best solution---albert Einstein Copyright A. A. Frempog Fermat's Last Theorem Proved o a Sigle Page "5% of the people thik; 0% of the people thik that they thik; ad the other 85% would rather die tha thik."----thomas Ediso "The simplest

More information

An improved car-following model considering variable safety headway distance. Yu-han Jia a,*, Jian-ping Wu a

An improved car-following model considering variable safety headway distance. Yu-han Jia a,*, Jian-ping Wu a A improved ar-followig model osiderig variable safety headway distae Yu-ha Jia a,*, Jia-pig Wu a ( a Departmet of Civil Egieerig, Tsighua Uiversity, Beijig 00084, Chia) Abstrat: Cosiderig high speed followig

More information

Chapter 2 Solutions. Prob. 2.1 (a&b) Sketch a vacuum tube device. Graph photocurrent I versus retarding voltage V for several light intensities.

Chapter 2 Solutions. Prob. 2.1 (a&b) Sketch a vacuum tube device. Graph photocurrent I versus retarding voltage V for several light intensities. Chapter Solutios Prob..1 (a&b) Sketh a vauum tube devie. Graph photourret I versus retardig voltage V for several light itesities. I light itesity V o V Note that V o remais same for all itesities. ()

More information

Non Linear Dynamics of Ishikawa Iteration

Non Linear Dynamics of Ishikawa Iteration Iteratioal Joural of Computer Appliatios (975 8887) Volume 7 No. Otober No Liear Dyamis of Ishiawa Iteratio Rajeshri Raa Asst. Professor Applied Siee ad Humaities Departmet G. B. Pat Egg. College Pauri

More information

Production Test of Rotary Compressors Using Wavelet Analysis

Production Test of Rotary Compressors Using Wavelet Analysis Purdue Uiversity Purdue e-pubs Iteratioal Compressor Egieerig Coferee Shool of Mehaial Egieerig 2006 Produtio Test of Rotary Compressors Usig Wavelet Aalysis Haishui Ji Shaghai Hitahi Eletrial Appliatio

More information

Organic Electronic Devices

Organic Electronic Devices Orgai letroi Devies Week 3: Charge rasport Leture 3.3: Multiple rap ad Release MR Model Brya W. Boudouris Chemial gieerig Purdue Uiversity Leture Overview ad Learig Objetives Coepts to be Covered i this

More information

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig

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

One way Analysis of Variance (ANOVA)

One way Analysis of Variance (ANOVA) Oe way Aalysis of Variae (ANOVA) ANOVA Geeral ANOVA Settig"Slide 43-45) Ivestigator otrols oe or more fators of iterest Eah fator otais two or more levels Levels a be umerial or ategorial ifferet levels

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

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

An Introduction to the Theory of Imprecise Soft Sets

An Introduction to the Theory of Imprecise Soft Sets IJ Itelliget Systems ad Appliatios 22 75-83 Published Olie Otober 22 i MECS (http://wwwmes-pressorg/) DOI: 585/isa229 A Itrodutio to the Theory of Impreise Soft Sets Tridiv Jyoti eog Dept of Mathematis

More information

Optimal Management of the Spare Parts Stock at Their Regular Distribution

Optimal Management of the Spare Parts Stock at Their Regular Distribution Joural of Evirometal Siee ad Egieerig 7 (018) 55-60 doi:10.1765/16-598/018.06.005 D DVID PUBLISHING Optimal Maagemet of the Spare Parts Stok at Their Regular Distributio Svetozar Madzhov Forest Researh

More information

Equivalence of the empirical shifted Deng-Fan oscillator potential for diatomic molecules

Equivalence of the empirical shifted Deng-Fan oscillator potential for diatomic molecules Equivalee of the empirial shifted Deg-Fa osillator potetial for diatomi moleules M. Hamzavi *, S. M. Ikhdair, K.-E. Thylwe 3 Departmet of Basi Siees, Shahrood Brah, Islami Azad Uiversity, Shahrood, Ira

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

Chaoyang University of Technology -- General KT Theory --

Chaoyang University of Technology -- General KT Theory -- Departmet of Costrutio Egieerig Advaed Soil Mehais Chaoyag iversity of Tehology -- Geeral KT Theory -- NIT 3 APPLICATION OF TERZAGHI S THEORY OF ONE DIMENSIONAL CONSOLIDATION TO PROBLEMS INVOLVING VARIOS

More information

Kinetics of Complex Reactions

Kinetics of Complex Reactions Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet

More information

Study on Solution of Non-homogeneous Linear Equation based on Ordinary Differential Equation Driving Jing Zhang

Study on Solution of Non-homogeneous Linear Equation based on Ordinary Differential Equation Driving Jing Zhang Iteratioal Coeree o Automatio Meaial Cotrol ad Computatioal Egieerig AMCCE 05 Stud o Solutio o No-omogeeous Liear Equatio based o Ordiar Dieretial Equatio Drivig Jig Zag Meaial ad Eletrial Egieerig College

More information

RENEWAL AND AVAILABILITY FUNCTIONS

RENEWAL AND AVAILABILITY FUNCTIONS Advaes ad Appliatios i Statistis 14 Pushpa Publishig House, Allahabad, Idia Available olie at http://pphmj.om/jourals/adas.htm Volume 43, Number 1, 14, Pages 65-89 RENEWAL AND AVAILABILITY FUNCTIONS Dilu

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

STK4011 and STK9011 Autumn 2016

STK4011 and STK9011 Autumn 2016 STK4 ad STK9 Autum 6 ypothesis testig Covers (most of) the followig material from hapter 8: Setio 8. Setios 8.. ad 8..3 Setio 8.3. Setio 8.3. (util defiitio 8.3.6) Ørulf Borga Departmet of Mathematis Uiversity

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

Supplement S1: RNA secondary structure. structure + sequence format

Supplement S1: RNA secondary structure. structure + sequence format Supplemet S1: RN seodary struture RN struture is ofte expressed shematially y its ase pairig: the Watso-rik (W) ase pairs (deie) with (rail), ad G (Guaie) with (ytosie) ad also the o-watso-rik (o-w) ase

More information

Optimal design of N-Policy batch arrival queueing system with server's single vacation, setup time, second optional service and break down

Optimal design of N-Policy batch arrival queueing system with server's single vacation, setup time, second optional service and break down Ameria. Jr. of Mathematis ad Siees Vol., o.,(jauary Copyright Mid Reader Publiatios www.jouralshub.om Optimal desig of -Poliy bath arrival queueig system with server's sigle vaatio, setup time, seod optioal

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