4.1 Solving Maxwell s equations in linear, isotropic, homogeneous, circularly symmetric media
|
|
- Valerie White
- 5 years ago
- Views:
Transcription
1 Maxwell s Equations in Cylindrical Coordinates Masud Mansuripur College of Optical Sciences, The University of Arizona, Tucson, Arizona 8571 Introduction In many systems of practical interest, electromagnetic waves propagate within linear, isotropic, homogeneous, circularly symmetric media. Examples include hollow-tube microwave waveguides, optical fibers, fiber lasers, nano-wires, and nano-rods. In this chapter we discuss the solution of Maxwell s equations for systems that exhibit circular symmetry around a given axis. We then proceed to examine the fundamental characteristics of these solutions, which are generally referred to as the modes (or eigenmodes) of the system. 4.1 Solving Maxwell s equations in linear, isotropic, homogeneous, circularly symmetric media An electromagnetic wave in a system with cylindrical symmetry has eigenfunctions of the form E(r,φ,z,t) = [E r (r)r + E φ (r)φ + E z (r)z ]exp(imφ)exp(ik o σ z z)exp( iωt), with a similar expression for the magnetic field H(r,φ,z,t). Here k o =ω/c= /λ o, where λ o is the vacuum wavelength; the integer m is the azimuthal mode number; and the complex-valued σ z is the propagation constant along the z-axis. In general, the beam resides within an environment having complex permittivity and permeability ε o ε(ω) and µ o µ(ω). As usual, the speed of light in vacuum c = 1/ µ o ε o and the impedance of free space Z o = µ o /ε o. In the absence of free charges and free currents, Maxwell s equations are written E r + r E r / r + ime φ + ik o rσ z E z = 0, (1) εe r = σ z Z o H φ (mz o /k o r)h z, εe φ = σ z Z o H r (iz o /k o ) H z / r, εe z = (mz o /k o r)h r + (iz o /k o r)h φ + (iz o /k o ) H φ / r, (a) (b) (c) µz o H r = (m/k o r)e z σ z E φ, (3a) µz o H φ = σ z E r + (i/k o ) E z / r, (3b) µz o H z = (m/k o r)e r (i/k o r)e φ (i/k o ) E φ / r, (3c) H r + r H r / r + imh φ + ik o rσ z H z = 0. (4) Substituting for E r and E φ from Eqs.(a) and (b) into Eqs.(3a) and (3b) yields (µε σ z )H r = (mε /Z o k o r)e z + i(σ z /k o ) H z / r, (5) (µε σ z )H φ = (mσ z /k o r)h z + i(ε/k o Z o ) E z / r. (6) Placing the above expressions for H r and H φ in Eq.(c) results in the following nd order differential equation for E z :
2 E z / r + (1/r) E z / r +[k o (µε σ z ) (m/r) ]E z = 0. (7) We define the radial propagation constant σ r = µε σ z, so that the radial coordinate r may be written in normalized form as ρ = k o σ r r. [For given values of µ, ε, and σ z, we shall always choose the value of σ r such that Im(σ r ) 0.] Equation (7) may now be written as follows: E z / ρ + (1/ρ) E z / ρ + (1 m /ρ )E z = 0. (8) In similar fashion, we substitute for H r and H φ from Eqs.(3a) and (3b) into Eqs.(a) and (b) to obtain (µε σ z )E r = i(σ z /k o ) E z / r (mµ Z o /k o r)h z, (9) (µε σ z )E φ = (mσ z /k o r)e z i(µz o /k o ) H z / r. (10) Placing the above expressions for E r and E φ in Eq.(3c) results in the following nd order differential equation for H z : H z / r + (1/r) H z / r +[k o (µε σ z ) (m/r) ]H z = 0. (11) Once again, use of the normalized radial coordinate ρ = k o σ r r leads to the standard Bessel equation for H z (ρ), that is, H z / ρ + (1/ρ) H z / ρ +(1 m /ρ )H z = 0. (1) We mention in passing that Eq.(1) could also have been obtained by substituting H r and H φ of Eqs.(5) and (6) into Eq.(4). Similarly, Eq.(8) could have been obtained by substituting E r and E φ of Eqs.(9) and (10) into Eq.(1). This should not be surprising, considering that Maxwell s equations are interdependent. The solutions of Eqs.(8) and (1) may now be expressed in terms of Bessel functions of the radial distance r from the z-axis. In general, any linear combination of Bessel functions of the first and second kind, order m, namely, J m (ρ) and Y m (ρ), satisfies the above equations. Thus the general form of the cylinder function appearing in the solution to Eqs.(8) and (1) will be C m (ρ) = A J m (ρ)+by m (ρ), where A and B are arbitrary complex constants. The derivative of the cylinder function with respect to its complex argument ρ, denoted by C m (ρ), will appear in the expressions for H r, H φ, E r, and E φ, as can be readily observed by examining Eqs.(5), (6), (9), and (10). The following identities will be useful: ρ C m (ρ) + ρc m (ρ) + (ρ m )C m (ρ) = 0, C m (ρ) = (ρ/m)[c m 1 (ρ) + C m+1 (ρ)], C m (ρ) = C m 1 (ρ) (m/ρ)c m (ρ), C m (ρ) = (m/ρ)c m (ρ) C m+1 (ρ), C m (ρ) = ( 1) m C m (ρ), C m (ρ) = ( 1) m C m (ρ). (13a) (13b) (13c) (13d) (13e) (13f) We divide the solutions to Eqs.(1-4) into two sets; for one set, we let E z = 0, in which case the solution will be referred to as Transverse Electric (TE); for the other, we let H z = 0, and call the solution Transverse Magnetic (TM). Occasionally, it will be possible to express the complete solution to a given problem in terms of either TE or TM modes; this occurs, for instance, when
3 σ z = 0 or when m = 0. Often, however, the boundary conditions cannot be matched with just one set of solutions; in such cases both TE and TM modes appear in the solution to Maxwell s equations, and the resulting mode will be referred to as EH or HE, depending on whether E z or H z dominates the longitudinal field component. The TE mode wavefunctions are written below: TE: E m (r) = (mµ /k o rσ r )C m (k o σ r r)r (iµ /σ r )C m (k o σ r r)φ, (14a) Z o H m (r) = (iσ z /σ r )C m (k o σ r r)r (mσ z /k o rσ r )C m (k o σ r r)φ + C m (k o σ r r)z. (14b) Similarly, the TM mode wavefunctions are given by: TM: E m (r) = (iσ z /σ r )C m (k o σ r r)r (mσ z /k o rσ r )C m (k o σ r r)φ + C m (k o σ r r)z, (15a) Z o H m (r) = (mε /k o rσ r )C m (k o σ r r)r + (iε /σ r )C m (k o σ r r)φ. (15b) The case of m = 0 may be simplified as follows: TE (m = 0): E 0 (r,φ,z) = (iµ/σ r )C 1 (k o σ r r)exp(ik o σ z z)φ, (16a) Z o H 0 (r,φ,z) = [ (iσ z /σ r )C 1 (k o σ r r)r + C 0 (k o σ r r)z ]exp(ik o σ z z). (16b) TM (m = 0): E 0 (r,φ,z) = [ (iσ z /σ r )C 1 (k o σ r r)r + C 0 (k o σ r r)z ]exp(ik o σ z z), (17a) Z o H 0 (r,φ,z) = (iε/σ r )C 1 (k o σ r r)exp(ik o σ z z)φ. (17b) In general, the field components in the xy-plane, E m and H m, may be expressed in terms of a superposition of right- and left-circularly polarized states, as follows: TE: E m (r,φ,z) = (µ /σ r ){C m 1 (k o σ r r)exp[i(m 1)φ](x +iy )+C m+1 (k o σ r r)exp[i(m+1)φ](x iy )} exp(ik o σ z z), (18a) Z o H m (r,φ,z)={(iσ z /σ r )[C m 1 (k o σ r r)exp[i(m 1)φ](x +iy ) C m+1 (k o σ r r)exp[i(m+1)φ](x iy )] + C m (k o σ r r)exp(imφ)z }exp(ik o σ z z). (18b) TM: E m (r,φ,z)= {(iσ z /σ r )[C m 1 (k o σ r r)exp[i(m 1)φ](x +iy ) C m+1 (k o σ r r)exp[i(m+1)φ](x iy )] + C m (k o σ r r)exp(imφ)z }exp(ik o σ z z), (19a) Z o H m (r,φ,z) = (ε/σ r ){C m 1 (k o σ r r)exp[i(m 1)φ](x +iy )+C m+1 (k o σ r r)exp[i(m+1)φ](x iy )} exp(ik o σ z z). (19b) When the region of interest includes the z-axis, the only acceptable cylinder function in the above equations will be C m (ρ) = J m (ρ); this is because Y m (ρ) as ρ 0. However, if the region of interest excludes the z-axis, the cylinder function will be a superposition of H m (1) (ρ)= 3
4 J m (ρ)+ iy m (ρ), which is an outgoing wave (i.e., one that propagates away from the cylinder axis), and H m () (ρ)= J m (ρ) iy m (ρ), which is an incoming wave. Bessel function J m ( ) as a superposition of plane-waves Consider the following p-polarized (TM) plane-wave within a homogeneous, isotropic medium specified by its (ε, µ) parameters: E(r,φ,z,t) = (E x x + E y y + E z z )exp[ik o (σ x x + σ y y + σ z z)]exp( iωt), H(r,φ,z,t) = (H x x + H y y )exp[ik o (σ x x + σ y y + σ z z)]exp( iωt). (0a) (0b) Here (x, y) = r (cosφ x + sinφ y ), and σ r = σ x x + σ y y = (σ x + iσ x )x + (σ y + iσ y )y = (σ x x + σ y y ) + i(σ x x +σ y y ) = σ r + iσ r. If σ r and σ r happen to be oriented in the same direction within the xyplane, say, at an angle θ relative to the x-axis, then σ r = (σ r + iσ r )(cosθ x + sinθ y ). Maxwell s equations then demand that σ r +σ z =µε; moreover, the field amplitudes must be related as follows: E(r,φ,z,t) = [ (σ z /σ r )(cosθ x + sinθ y ) + z ]exp{ik o [σ r r cos(θ φ)+σ z z ct]}, H(r,φ,z,t) = (ε /Z o σ r )(sinθ x cosθ y )exp{ik o [σ r r cos(θ φ)+σ z z ct]}. (1a) (1b) In these equations the field amplitudes have been normalized to yield a magnitude of unity for E z. Next, we fix σ z and σ r, and consider a superposition of all such plane-waves covering the range of angles θ from 0 to. For the amplitude distribution of these plane-waves we choose f (θ ) = ( i m ) 1 exp(imθ ), with m being an arbitrary integer. The z-component of the resulting E- field will then be E z (r,φ) = ( i m ) 1 0 exp(imθ )exp[iko σ r r cos(θ φ)]dθ = () 1 exp(imφ) 0 exp(imθ iko σ r r sinθ )dθ = () 1 [ 0 exp(imθ iko σ r r sinθ )dθ + exp(imθ iko σ r r sinθ )dθ ]exp(imφ). () A change of variable from θ to θ reveals the second integral to be the complex conjugate of the first; we thus have E z (r,φ) = [(1/) 0 cos(mθ ko σ r r sinθ )dθ]exp(imφ) = J m (k o σ r r)exp(imφ). (3a) For the component of the E-field in the xy-plane we then have E (r,φ) = (σ z /σ r )( i m ) 1 0 exp(imθ )(cosθ x + sinθ y)exp[iko σ r r cos(θ φ)]dθ = (σ z /σ r )( i m ) 1 {exp[i(m 1)θ ](x + iy )+exp[i(m+1)θ ](x iy )}exp[ik o σ r r cos(θ φ)]dθ = (iσ z /σ r ){J m 1 (k o σ r r)exp[i(m 1)φ ](x + iy ) J m+1 (k o σ r r)exp[i(m+1)φ ](x iy )}. (3b) 4
5 The above expressions for the various E-field components are in agreement with Eq.(19a), corresponding to a TM-polarized cylindrical wave. The corresponding expression for the H-field may be similarly derived from Eq.(1b). Derivation of the expressions for a TE-polarized cylindrical wave is straightforward. Hankel functions H m (1,) (ρ) as superpositions of plane-waves In general, for a plane-wave such as that described by Eqs. (0), the radial propagation constant may be written as σ r = (σ r + iσ r )(cosθ x + sinθ y ), where σ r = σ r σ r = σ r + iσ r = µε σ z is the (complex) length of the vector σ r. Thus far we have considered only real-valued angles θ, but θ can be complex-valued as well. For example, Fig.1 shows a path in the complex θ-plane that, starting from ½ + i, goes down along a vertical leg to the real axis, covers the range ( ½, ½) on this axis, then descends along a second vertical leg to ½ i. (The complex magnitude of the vector σ r does not vary as θ traces an arbitrary trajectory in the complex θ- plane, because sin θ + cos θ is always equal to unity, irrespective of whether θ is real or complex.) Since σ r = σ x x + σ y y = (σ x x + σ y y ) + i(σ x x + σ y y ), it is clear that Real(σ r ) = σ x x + σ y y, a real-valued vector in the xy-plane, defines the direction of phase propagation, whereas Imag(σ r ) = σ x x + σ y y, another real-valued vector in the xy-plane, specifies the direction of (1) attenuation. The Hankel function H m is obtained from a superposition of the plane-waves located on a trajectory such as T 1 of Fig.1(a), when the complex-amplitude distribution for these plane-waves is ( i m ) 1 exp(imθ). In other words, H m (1) (k o σ r r)exp(imφ) = ( i m ) 1 T1 exp(imθ )exp[ik o σ r r cos(θ φ)]dθ. (4) A horizontal shift of T 1 by φ does not alter the value of the above integral, provided that ½ <φ <½, which is needed to ensure that the integrand goes to zero at the initial and final points of T 1. In fact, any deformation of T 1 will not affect the resulting integral, so long as the trajectory s start and finish points remain within (,0)+ i and (0,) i, respectively. We conclude that H m (1) (k o σ r r) = ( i m ) 1 T1 exp(imθ )exp(ik o σ r r cosθ )dθ. (5) Therefore, in the region x > 0 of the xy-plane, where ½ <φ <½, Eq. (5) describes H m (1) (k o σ r r) as a superposition of outgoing plane-waves. Similar considerations apply to H m () (k o σ r r) with the exception that the θ-plane trajectory in this case is T of Fig.1(b). The plane-waves that constitute type Hankel functions are thus seen to be plane-waves that propagate toward the center of the coordinate system. For a generalized plane-wave, given the z-component k o σ z of the propagation vector, the radial phase-factor may be written as follows: ik o σ r r = ik o (σ r + iσ r )(cosθ x + sinθ y ) r(cosφ x + sinφ y ) = ik o r(σ r + iσ r )cos(θ + iθ φ) = k o r[σ r cos(θ φ)chθ σ r sin(θ φ)shθ ]+ ik o r[σ r cos(θ φ)chθ +σ r sin(θ φ)shθ ]. (6) 5
6 When θ +, shθ and chθ approach each other on their way to + ; it is thus required for the real part of the above expression to be negative, that is, σ r cos(θ φ) σ r sin(θ φ)>0. When θ, however, we must have σ r cos(θ φ) +σ r sin(θ φ)>0. Let us define the real-valued angle θ o = tan 1 (σ r /σ r ); clearly, since σ r 0, we have 0 θ o ½ when σ r 0, and ½ θ o 0 when σ r 0. The starting point of the trajectory must, therefore, lie in the region (θ o,θ o )+i when θ o 0 and in (θ o,θ o +)+i when θ o 0. Similarly, the trajectory s terminal point must be in ( θ o, θ o ) i when θ o 0, and in ( θ o, θ o ) i when θ o 0; see Fig.. θ θ T 1 Complex θ -plane (θ =θ +iθ ) T ½ ½ θ ½ θ Fig. 1. (a) Trajectory T 1 in the complex θ-plane for type 1 Hankel function H m (1) (ρ). (b) Trajectory T for type Hankel function H m () (ρ). Any deformation of these trajectories is allowed, so long as the start and finish points remain in the corresponding shaded regions. Note that either trajectory may be translated horizontally by an amount φ, where ½ <φ <½, without affecting the resultant integral. θ θ T 1 θ o = tan 1 (σ r /σ r ) 0 θ o ½ ½ θ o 0 T 1 θ o θ o θ o θ o θ θ o θ o θ o θ o + θ Fig.. Trajectory T 1 in the complex θ-plane for the defining integral of H m (1) (ρ) in the general case when σ r = σ r + iσ r. In general σ r 0. In (a) σ r 0, whereas in (b) σ r 0. 6
7 Figure 3 shows, for a given value of σ r = σ r + iσ r, the trajectories of Real(σ r ) and Imag(σ r ), corresponding to the θ-trajectory T 1 depicted in Fig.(a). While the phase-fronts of different plane-waves could propagate in various directions, their amplitudes essentially decay along the horizontal axis. Within the x > 0 half-space, the superposition of all these plane-waves produces, in accordance with Eq.(4), the function H m (1) (k o σ r r)exp(imφ)exp(ik o σ z z). Note that, since θ is real-valued on the circular sections of the trajectories, Real(σ r ) and Imag(σ r ) are parallel to each other on the circular trajectories. In contrast, on the vertical legs of T 1, the trajectories of Real(σ r ) and Imag(σ r ) are nearly orthogonal to each other. σ x x + σ y y σ y σ x x + σ y y σ x Fig. 3. σ x σ y -plane trajectories of Real(σ r ) (red) and Imag(σ r ) (blue) for the case of µε =.5+0.7i and σ z = 0.5, corresponding to σ r = i and θ o 8. The θ-plane integration path is depicted in Fig.(a). 4. Guided modes and surface-plasmon-polaritons in systems of cylindrical symmetry The analysis presented in this section is applicable whenever the medium hosting the electromagnetic wave consists of two regions: inside a cylinder of radius r o the material parameters are (ε 1,µ 1 ), while outside the cylinder, the parameters are (ε,µ ). Examples include step-index optical fibers supporting guided modes, and metallic nano-wires or nano-rods (with or without a dielectric coating) that host surface-plasmon-polariton excited waves. In general, a superposition of TE and TM waves in both regions (i.e., inside and outside the cylinder of radius r o ) is needed to satisfy the continuity of E z, H z, E φ, H φ at the r = r o boundary. Inside the cylinder, the appropriate cylinder function is J m (k o σ r1 r), and the (initially unknown) coefficients for TM and TE modes are A 1 and B 1, respectively. Outside the cylinder, the appropriate function is H m (1) (k o σ r r), while the corresponding coefficients are A and B. The continuity equations are thus written Continuity of E z : A 1 J m (k o σ r1 r o ) = A H m (1) (k o σ r r o ), (7a) Continuity of H z : B 1 J m (k o σ r1 r o ) = B H m (1) (k o σ r r o ), (7b) 7
8 Continuity of E φ : A 1 (mσ z /k o r o σ r1 ) J m (k o σ r1 r o ) + B 1 (iµ 1 /σ r1 ) J m (k o σ r1 r o ) = A (mσ z /k o r o σ r )H m (1) (k o σ r r o ) + B (iµ /σ r )H m (1) (k o σ r r o ), (7c) Continuity of H φ : A 1 (iε 1 /σ r1 ) J m (k o σ r1 r o ) B 1 (mσ z /k o r o σ r1 ) J m (k o σ r1 r o ) = A (iε /σ r )H m (1) (k o σ r r o ) B (mσ z /k o r o σ r )H m (1) (k o σ r r o ). (7d) The continuity of D r = ε o εe r is guaranteed by the continuity of H z and H φ ; similarly, the continuity of B r = µ o µh r is guaranteed by the continuity of E z and E φ. The above equations can be solved for A 1, B 1, A, B provided that the following characteristic equation is satisfied: µ 1 µ 1 ε 1 σ z J m (k o σ r1 r o ) J m (k o σ r1 r o ) µ H m (1) (k o σ r r o ) µ ε σ z H (1) m (k o σ r r o ) (8) ε 1 µ 1 ε 1 σ z J m (k o σ r1 r o ) J m (k o σ r1 r o ) ε H m (1) (k o σ r r o ) µ ε σ z H m (1) (k o σ r r o ) = m σ z 1 k o r o µ 1 ε 1 σ z 1 µ ε σ z Given the values of m, k o, r o, µ 1, ε 1, µ, ε, one may solve Eq. (8) numerically for acceptable values of σ z. In lossless media (i.e., real-valued µ 1, ε 1, µ, ε ), if a real-valued σ z is found, it will correspond to a guided mode (e.g., in glass fibers, where µ ε <σ z < µ 1 ε 1 ). In general, however, σ z will be complex, corresponding to an attenuated mode along the z-axis. A special case would arise when σ z = 0; here TE and TM modes are decoupled, and the characteristic equation for each may be written independently of the other, as follows: TE: J m (k o µ 1 ε 1 r o ) µ 1 /ε 1 Jm (k o µ 1 ε 1 r o ) = µ /ε H m (1) (k o µ ε r o ) H m (1) (k o µ ε r o ) (9) TM: J m (k o µ 1 ε 1 r o ) ε 1 /µ 1 Jm (k o µ 1 ε 1 r o ) = ε /µ H m (1) (k o µ ε r o ) H m (1) (k o µ ε r o ) This case probably never arises in practice, because either Eq.(9) or Eq.(30) can be satisfied only by coincidence, when the cylinder radius r o happens to be just right for the two media inside and outside the cylinder. Even then it is hard to see how energy can be conserved, because the entire wave is a standing wave in cylindrical coordinates. One should investigate the properties of the Poynting vector in order to gain further insight into this type of problem. Another special case is when m = 0; here once again TE and TM modes are decoupled, and the characteristic equation for each mode becomes: (30) TE: µ 1 J 1 (k o µ 1 ε 1 σ z r o ) µ H (1) 1 (k o µ ε σ z r o ) = µ 1 ε 1 σ z J 0 (k o µ 1 ε 1 σ z r o ) µ ε σ z H (1) 0 (k o µ ε σ z r o ) (31) 8
9 TM: ε 1 J 1 (k o µ 1 ε 1 σ z r o ) ε H (1) 1 (k o µ ε σ z r o ) = µ 1 ε 1 σ z J 0 (k o µ 1 ε 1 σ z r o ) µ ε σ z H (1) 0 (k o µ ε σ z r o ) (3) As before, given the values of k o, r o, µ 1, ε 1, µ, ε, one can solve Eq.(31) or Eq.(3) numerically to determine the acceptable values of σ z. 9
10 4.3 Energy flux and the Poynting vector For the TE and TM modes described by Eqs.(14) and (15), the time-averaged Poynting vector < S > = ½Real(E H*) is given by TE: < S r > = exp[ k o Im(σ z )z] m Im(µ * σ Im[µ * σ r C m (k o σ r r)c * r ) m+1(k o σ r r)] Z o µε σ z k o r µε σ z C m(k o σ r r), (33a) < S φ > = m Re[µ /(µε σ z )]exp[ k o Im(σ z )z] Z o k o r C m (k o σ r r), (33b) < S z > = Re(µ * σ z )exp[ k o Im(σ z )z] C m (k o σ r r) + Z o µε σ z k o r µε σ z C m(k o σ r r) m (33c) TM: < S r > = exp[ k o Im(σ z )z] Z o µε σ z Im[ε * σ r C m (k o σ r r)c * m+1(k o σ r r)] m Im(ε * σ r ) k o r µε σ z C m(k o σ r r), (34a) < S φ > = m Re[ε /(µε σ z )]exp[ k o Im(σ z )z] Z o k o r C m (k o σ r r), (34b) < S z > = Re(ε * σ z )exp[ k o Im(σ z )z] C m (k o σ r r) + Z o µε σ z k o r µε σ z C m(k o σ r r) m (34c) When a TE mode of amplitude A is superimposed on a TM mode of amplitude B, the following cross-terms will appear in the expression of the Poynting vector: Cross terms: m Im(A * B) < S r > = Z o k o r exp[ k o Im(σ z )z]im[σ z /(µε σ z )] C m (k o σ r r), (35a) < S φ > = Im(A * B) Z o exp[ k o Im(σ z )z] Re[(σ z /σ r )C * m(k o σ r r)c m (k o σ r r)], (35b) < S z > = m Im[A * B(µ * ε + σ z )] exp[ k o Im(σ z )z] Re[C m (k o σ r r)c * m (k o σ r r)/σ r ]. Z o k o r µε σ z (35c) The above equations will be useful when analyzing the behavior of optical fibers, hollowtube waveguides, nano-wires, nano-rods, and many similar systems of cylindrical symmetry. The analysis of propagating or guided modes in such systems is highly complex and specialized, 10
11 requiring a textbook of its own. We do hope, however, that the present chapter has provided the reader with the necessary background to facilitate his/her future forays into the subject. 11
12 References 1. For an extensive listing of the various properties of Bessel functions see I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, seventh edition, Academic Press, For applications of modal analysis to step-index optical fibers see G. Keiser, Optical Fiber Communications, 3 rd edition, McGraw-Hill, New York, For applications of modal analysis to cylindrical waveguides see S. Ramo, J. A. Whinnery, and T. Van Duzer, Fields and Waves in Communications Electronics, Wiley, New York, For a discussion of surface-plasmon-polaritons see C. A. Pfeiffer, E. N. Economou, and K. L. Ngai, Surface polaritons in a circularly cylindrical interface: Surface plasmons, Physical Review B 10, (1974). 1
13 Appendix Bessel Functions and Their Properties Bessel functions are solutions of the second-order differential equation Gradshteyn & Ryzhik d yx ( ) 1 d yx ( ) ν + + (1 ) yx ( ) = 0. (A1) dx x dx x Although the function y( ), the variable x, and the constant parameter ν are generally complex-valued, we shall primarily be interested in cases where y( ), x, and ν are real. For a given value of ν, the nd order equation (A1) has two independent solutions: Bessel functions of the first kind, Jν ( x), and Bessel functions of the second kind, Yν ( x); the latter are sometimes referred to as Neumann functions. Bessel functions of the third kind, types 1 and, also known as Hankel functions, are constructed from J ν ( ) and Y ν ( ) as follows: H ν (1) (x)=jν (x)+iy ν (x); H ν () (x)=jν (x) iy ν (x). (Aa) (Ab) In what follows, we list some of the important properties of the Bessel functions without attempting to prove them; such proofs may be found in any standard treatise on Bessel functions. For ν = n 0 (i.e., ν a non-negative integer), the solutions of Eq.(A1) are given by G&R 8.40 k k n ( 1) ( x/) Jn( x) = ( x/) n= 0,1,, 3... (A3) k!( n+ k)! k = 0 G&R n n 1 ( x/) ( n k 1)! c n k = 0 k! Yn( x) = [ + ln( x/)] J ( x) ( x/) k x x x + ( n!) k= 1 k k= 1 k!( k+ n)! m= 1m m= 1m n ( /) n n k k 1 ( /) ( 1) ( /) n + k 1 k 1. (A4) In the above equations, 0! is defined as 1. Also, when the upper limit of a sum over k is below its lower limit, the term containing the sum must be set to zero; this occurs for n = 0 in the nd and 3 rd terms on the right-hand-side of Eq.(A4), thus reducing the expression for Y 0 (x) to one containing the 1 st and 4 th terms only. Finally, the Euler constant c appearing in Eq.(A4) is defined as n 1 c = lim ln n n = k= 1 k (A5) Given below are expanded forms of J( x), Y( x), J( x), and Y ( ); 1 x plots of these functions appear in Fig. A ( x/) ( x/) ( x/) J... 0( x ) = 1 + +, (A6) (1!) (!) (3!) 13
14 ( x/) ( x/) ( x/) J... 1( x) = x + +, (A7) 1!!!3! 3!4! 4 6 ( x/) 1 ( x/) 1 1 ( x/)... Y0( x) = [ c+ ln ( x/)] J0( x) + (1 ) (1 ), (1!) (!) (3!) (A8) k k k+ 1 k x ( 1) ( x/) 1 1 Y1( x) = [ c+ ln( x/)] J1( x) x k= 1 k! ( k+ 1)! m= 1m m= 1m (A9) Fig.A1. (a) Plots of J 0 (x) and J 1 (x). (b) Plots of Y 0 (x) and Y 1 (x). Changing the sign of a Bessel function s argument modifies the function in the following way: J n ( x) = ( 1) n J n (x), Y n ( x)= ( 1) n [Y n (x)+i J n (x)]. When n is fixed and x 0, the small-argument limiting forms of the Bessel functions are (A10) (A11) n ( x/) Jn( x) ~ ; n 0, (A1) n! ( n 1)! n Yn ( x) ~ ( x/) ; n 1, (A13) Y0( x) ~ ln x. Similarly, when x, the large-argument limiting forms become (A14) G&R J ( x) /( x) cos[ x ( n/) ( /4)], (A15) n G&R Y( x) /( x) sin[ x ( n/) ( /4)]. (A16) n The first-derivatives of the Bessel functions are given by G&R G&R J' ( x) = J ( x), (A17) 0 1 Y' ( x) = Y( x), (A18)
15 G&R , G&R , J' ( x) = J ( x) ( n/ x) J ( x) = J ( x) + ( n/ x) J ( x), (A19) n n 1 n n+ 1 n Y' ( x) = Y ( x) ( nxy / ) ( x) = Y ( x) + ( n/ xy ) ( x). (A0) n n 1 n n+ 1 n The modified Bessel functions of imaginary argument are the real-valued functions I n ( ) and K n ( ) defined as follows: I n (x) = i n J n (ix), K n (x)=i n+1 ( /) H n (1) (ix). Figure A shows plots of the functions I 0 (x), I 1 (x), K 0 (x), and K 1 (x). G&R G&R (A1) (A) Fig.A. Plots of I 0 (x), I 1 (x), K 0 (x), and K 1 (x). Two useful identities involving Bessel functions of different kinds and different orders are Yn( x) Jn+ 1( x) Jn( xy ) n+ 1( x) =, x G&R (A3) 1 In( xk ) n+ 1( x) + Kn( xi ) n+ 1( x) =. x G&R (A4) Finally, a frequently-encountered integral representation of the Bessel function of the first kind, n th order is given below. 0 cos(nx)exp(iβ cosx)dx = i n J n (β ). G&R (A5) For an extensive listing of the various properties of the Bessel functions see I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, seventh edition, Academic Press,
1 The formation and analysis of optical waveguides
1 The formation and analysis of optical waveguides 1.1 Introduction to optical waveguides Optical waveguides are made from material structures that have a core region which has a higher index of refraction
More informationWaves in Linear Optical Media
1/53 Waves in Linear Optical Media Sergey A. Ponomarenko Dalhousie University c 2009 S. A. Ponomarenko Outline Plane waves in free space. Polarization. Plane waves in linear lossy media. Dispersion relations
More informationLecture 3 Fiber Optical Communication Lecture 3, Slide 1
Lecture 3 Optical fibers as waveguides Maxwell s equations The wave equation Fiber modes Phase velocity, group velocity Dispersion Fiber Optical Communication Lecture 3, Slide 1 Maxwell s equations in
More informationFundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation
Fundamentals on light scattering, absorption and thermal radiation, and its relation to the vector radiative transfer equation Klaus Jockers November 11, 2014 Max-Planck-Institut für Sonnensystemforschung
More informationMultilayer Reflectivity
Multilayer Reflectivity John E. Davis jed@jedsoft.org January 5, 2014 1 Introduction The purpose of this document is to present an ab initio derivation of the reflectivity for a plane electromagnetic wave
More informationEigenfunction Expansions of Source-excited Electromagnetic Fields on Open Cylindrical Guiding Structures in Unbounded Gyrotropic Media
128 Eigenfunction Expansions of Source-excited Electromagnetic Fields on Open Cylindrical Guiding Structures in Unbounded Gyrotropic Media A. V. Kudrin 1, E. Yu. Petrov 1, and T. M. Zaboronkova 2 1 University
More informationElectromagnetic Waves
Electromagnetic Waves Maxwell s equations predict the propagation of electromagnetic energy away from time-varying sources (current and charge) in the form of waves. Consider a linear, homogeneous, isotropic
More informationElectromagnetic Theory for Microwaves and Optoelectronics
Keqian Zhang Dejie Li Electromagnetic Theory for Microwaves and Optoelectronics Second Edition With 280 Figures and 13 Tables 4u Springer Basic Electromagnetic Theory 1 1.1 Maxwell's Equations 1 1.1.1
More informationElectromagnetic Waves Across Interfaces
Lecture 1: Foundations of Optics Outline 1 Electromagnetic Waves 2 Material Properties 3 Electromagnetic Waves Across Interfaces 4 Fresnel Equations 5 Brewster Angle 6 Total Internal Reflection Christoph
More informationElectromagnetic Theory for Microwaves and Optoelectronics
Keqian Zhang Dejie Li Electromagnetic Theory for Microwaves and Optoelectronics Translated by authors With 259 Figures Springer Contents 1 Basic Electromagnetic Theory 1 1.1 Maxwell's Equations 1 1.1.1
More informationSurface Plasmon Polaritons on Metallic Surfaces
Surface Plasmon Polaritons on Metallic Surfaces Masud Mansuripur, Armis R. Zakharian and Jerome V. Moloney Recent advances in nano-fabrication have enabled a host of nano-photonic experiments involving
More informationLecture Notes on Wave Optics (03/05/14) 2.71/2.710 Introduction to Optics Nick Fang
Outline: A. Electromagnetism B. Frequency Domain (Fourier transform) C. EM waves in Cartesian coordinates D. Energy Flow and Poynting Vector E. Connection to geometrical optics F. Eikonal Equations: Path
More informationElectromagnetic fields and waves
Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell
More informationCylindrical Dielectric Waveguides
03/02/2017 Cylindrical Dielectric Waveguides Integrated Optics Prof. Elias N. Glytsis School of Electrical & Computer Engineering National Technical University of Athens Geometry of a Single Core Layer
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationTypical anisotropies introduced by geometry (not everything is spherically symmetric) temperature gradients magnetic fields electrical fields
Lecture 6: Polarimetry 1 Outline 1 Polarized Light in the Universe 2 Fundamentals of Polarized Light 3 Descriptions of Polarized Light Polarized Light in the Universe Polarization indicates anisotropy
More informationII Theory Of Surface Plasmon Resonance (SPR)
II Theory Of Surface Plasmon Resonance (SPR) II.1 Maxwell equations and dielectric constant of metals Surface Plasmons Polaritons (SPP) exist at the interface of a dielectric and a metal whose electrons
More information9 The conservation theorems: Lecture 23
9 The conservation theorems: Lecture 23 9.1 Energy Conservation (a) For energy to be conserved we expect that the total energy density (energy per volume ) u tot to obey a conservation law t u tot + i
More informationElectromagnetic optics!
1 EM theory Electromagnetic optics! EM waves Monochromatic light 2 Electromagnetic optics! Electromagnetic theory of light Electromagnetic waves in dielectric media Monochromatic light References: Fundamentals
More informationMAT389 Fall 2016, Problem Set 4
MAT389 Fall 2016, Problem Set 4 Harmonic conjugates 4.1 Check that each of the functions u(x, y) below is harmonic at every (x, y) R 2, and find the unique harmonic conjugate, v(x, y), satisfying v(0,
More informationElectromagnetic (EM) Waves
Electromagnetic (EM) Waves Short review on calculus vector Outline A. Various formulations of the Maxwell equation: 1. In a vacuum 2. In a vacuum without source charge 3. In a medium 4. In a dielectric
More informationIntroduction to Polarization
Phone: Ext 659, E-mail: hcchui@mail.ncku.edu.tw Fall/007 Introduction to Polarization Text Book: A Yariv and P Yeh, Photonics, Oxford (007) 1.6 Polarization States and Representations (Stokes Parameters
More informationPropagation of EM Waves in material media
Propagation of EM Waves in material media S.M.Lea 09 Wave propagation As usual, we start with Maxwell s equations with no free charges: D =0 B =0 E = B t H = D t + j If we now assume that each field has
More informationAzimuthally polarized surface plasmons as effective terahertz waveguides
Azimuthally polarized surface plasmons as effective terahertz waveguides Qing Cao and Jürgen Jahns Optische Nachrichtentechnik, FernUniversität Hagen, Universitätsstrasse 27/PRG, D-58084 Hagen, Germany
More informationHomework 1. Nano Optics, Fall Semester 2017 Photonics Laboratory, ETH Zürich
Homework 1 Contact: mfrimmer@ethz.ch Due date: Friday 13.10.2017; 10:00 a.m. Nano Optics, Fall Semester 2017 Photonics Laboratory, ETH Zürich www.photonics.ethz.ch The goal of this homework is to establish
More informationElectromagnetic Waves
Physics 8 Electromagnetic Waves Overview. The most remarkable conclusion of Maxwell s work on electromagnetism in the 860 s was that waves could exist in the fields themselves, traveling with the speed
More informationSummary of Beam Optics
Summary of Beam Optics Gaussian beams, waves with limited spatial extension perpendicular to propagation direction, Gaussian beam is solution of paraxial Helmholtz equation, Gaussian beam has parabolic
More information4. Integrated Photonics. (or optoelectronics on a flatland)
4. Integrated Photonics (or optoelectronics on a flatland) 1 x Benefits of integration in Electronics: Are we experiencing a similar transformation in Photonics? Mach-Zehnder modulator made from Indium
More informationGeneral review: - a) Dot Product
General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,
More informationElectromagnetic Scattering from an Anisotropic Uniaxial-coated Conducting Sphere
Progress In Electromagnetics Research Symposium 25, Hangzhou, China, August 22-26 43 Electromagnetic Scattering from an Anisotropic Uniaxial-coated Conducting Sphere You-Lin Geng 1,2, Xin-Bao Wu 3, and
More informationOverview in Images. S. Lin et al, Nature, vol. 394, p , (1998) T.Thio et al., Optics Letters 26, (2001).
Overview in Images 5 nm K.S. Min et al. PhD Thesis K.V. Vahala et al, Phys. Rev. Lett, 85, p.74 (000) J. D. Joannopoulos, et al, Nature, vol.386, p.143-9 (1997) T.Thio et al., Optics Letters 6, 197-1974
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More information1. Consider the biconvex thick lens shown in the figure below, made from transparent material with index n and thickness L.
Optical Science and Engineering 2013 Advanced Optics Exam Answer all questions. Begin each question on a new blank page. Put your banner ID at the top of each page. Please staple all pages for each individual
More informationECE 604, Lecture 17. October 30, In this lecture, we will cover the following topics: Reflection and Transmission Single Interface Case
ECE 604, Lecture 17 October 30, 2018 In this lecture, we will cover the following topics: Duality Principle Reflection and Transmission Single Interface Case Interesting Physical Phenomena: Total Internal
More informationFocal shift in vector beams
Focal shift in vector beams Pamela L. Greene The Institute of Optics, University of Rochester, Rochester, New York 1467-186 pgreene@optics.rochester.edu Dennis G. Hall The Institute of Optics and The Rochester
More informationElectromagnetic waves in free space
Waveguide notes 018 Electromagnetic waves in free space We start with Maxwell s equations for an LIH medum in the case that the source terms are both zero. = =0 =0 = = Take the curl of Faraday s law, then
More informationChapter 2 Basic Optics
Chapter Basic Optics.1 Introduction In this chapter we will discuss the basic concepts associated with polarization, diffraction, and interference of a light wave. The concepts developed in this chapter
More informationPhysics 506 Winter 2004
Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationPhysics 504, Lecture 9 Feb. 21, 2011
Last Latexed: February 17, 011 at 15:8 1 1 Ionosphere, Redux Physics 504, Lecture 9 Feb. 1, 011 Let us return to the resonant cavity formed by the surface of the Earth a spherical shell of radius r = R
More informationOn Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material Surfaces
Vol. 114 2008) ACTA PHYSICA POLONICA A No. 6 A Optical and Acoustical Methods in Science and Technology On Electromagnetic-Acoustic Analogies in Energetic Relations for Waves Interacting with Material
More informationChapter 33. Electromagnetic Waves
Chapter 33 Electromagnetic Waves Today s information age is based almost entirely on the physics of electromagnetic waves. The connection between electric and magnetic fields to produce light is own of
More informationElectromagnetic Waves. Chapter 33 (Halliday/Resnick/Walker, Fundamentals of Physics 8 th edition)
PH 222-3A Spring 2007 Electromagnetic Waves Lecture 22 Chapter 33 (Halliday/Resnick/Walker, Fundamentals of Physics 8 th edition) 1 Chapter 33 Electromagnetic Waves Today s information age is based almost
More information9. Dispersion relation of metal nanorods and nanotips
9. Dispersion relation of metal nanorods and nanotips D. E. Chang, A. S. Sørensen, P. R. Hemmer, and M. D. Lukin, Strong coupling of single emitters to surface plasmons, PR B 76,035420 (2007) M. I. Stockman,
More informationPROPAGATION IN A FERRITE CIRCULAR WAVEGUIDE MAGNETIZED THROUGH A ROTARY FOUR-POLE MAGNETIC FIELD
Progress In Electromagnetics Research, PIER 68, 1 13, 2007 PROPAGATION IN A FERRITE CIRCULAR WAVEGUIDE MAGNETIZED THROUGH A ROTARY FOUR-POLE MAGNETIC FIELD M. Mazur Analog Techniques Department Telecommunication
More informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationin Electromagnetics Numerical Method Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD
2141418 Numerical Method in Electromagnetics Introduction to Electromagnetics I Lecturer: Charusluk Viphavakit, PhD ISE, Chulalongkorn University, 2 nd /2018 Email: charusluk.v@chula.ac.th Website: Light
More informationarxiv: v1 [physics.optics] 30 Mar 2010
Analytical vectorial structure of non-paraxial four-petal Gaussian beams in the far field Xuewen Long a,b, Keqing Lu a, Yuhong Zhang a,b, Jianbang Guo a,b, and Kehao Li a,b a State Key Laboratory of Transient
More informationHomework 1. Nano Optics, Fall Semester 2018 Photonics Laboratory, ETH Zürich
Homework 1 Contact: mfrimmer@ethz.ch Due date: Friday 12 October 2018; 10:00 a.m. Nano Optics, Fall Semester 2018 Photonics Laboratory, ETH Zürich www.photonics.ethz.ch The goal of this homework is to
More informationScattering cross-section (µm 2 )
Supplementary Figures Scattering cross-section (µm 2 ).16.14.12.1.8.6.4.2 Total scattering Electric dipole, a E (1,1) Magnetic dipole, a M (1,1) Magnetic quardupole, a M (2,1). 44 48 52 56 Wavelength (nm)
More informationPhysics 3323, Fall 2014 Problem Set 13 due Friday, Dec 5, 2014
Physics 333, Fall 014 Problem Set 13 due Friday, Dec 5, 014 Reading: Finish Griffiths Ch. 9, and 10..1, 10.3, and 11.1.1-1. Reflecting on polarizations Griffiths 9.15 (3rd ed.: 9.14). In writing (9.76)
More informationREFLECTION AND REFRACTION AT A SINGLE INTERFACE
REFLECTION AND REFRACTION AT A SINGLE INTERFACE 5.1 THE BEHAVIOUR OF LIGHT AT A DIELECTRIC INTERFACE The previous Chapters have been concerned with the propagation of waves in empty space or in uniform,
More informationGreen s Tensors for the Electromagnetic Field in Cylindrical Coordinates
Chapter 25 Green s Tensors for the Electromagnetic Field in Cylindrical Coordinates The solutions of the vector Helmholtz equation in three dimensions can be expressed by a complete set of vector fields
More informationLow Losses Left Handed Materials Using Metallic Magnetic Cylinders.
Low Losses Left Handed Materials Using Metallic Magnetic Cylinders. N. García and E.V. Ponizovskaia Laboratorio de Física de Sistemas Pequeños y Nanotecnología, Consejo Superior de Investigaciones Científicas,
More informationLecture 19 Optical MEMS (1)
EEL6935 Advanced MEMS (Spring 5) Instructor: Dr. Huikai Xie Lecture 19 Optical MEMS (1) Agenda: Optics Review EEL6935 Advanced MEMS 5 H. Xie 3/8/5 1 Optics Review Nature of Light Reflection and Refraction
More informationAntennas and Propagation
Antennas and Propagation Ranga Rodrigo University of Moratuwa October 20, 2008 Compiled based on Lectures of Prof. (Mrs.) Indra Dayawansa. Ranga Rodrigo (University of Moratuwa) Antennas and Propagation
More informationElectromagnetic Waves
Electromagnetic Waves Our discussion on dynamic electromagnetic field is incomplete. I H E An AC current induces a magnetic field, which is also AC and thus induces an AC electric field. H dl Edl J ds
More informationElectromagnetism. Christopher R Prior. ASTeC Intense Beams Group Rutherford Appleton Laboratory
lectromagnetism Christopher R Prior Fellow and Tutor in Mathematics Trinity College, Oxford ASTeC Intense Beams Group Rutherford Appleton Laboratory Contents Review of Maxwell s equations and Lorentz Force
More informationE E D E=0 2 E 2 E (3.1)
Chapter 3 Constitutive Relations Maxwell s equations define the fields that are generated by currents and charges. However, they do not describe how these currents and charges are generated. Thus, to find
More informationMicrowave Phase Shift Using Ferrite Filled Waveguide Below Cutoff
Microwave Phase Shift Using Ferrite Filled Waveguide Below Cutoff CHARLES R. BOYD, JR. Microwave Applications Group, Santa Maria, California, U. S. A. ABSTRACT Unlike conventional waveguides, lossless
More informationHomogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different Galilean Reference Systems
Progress In Electromagnetics Research Symposium Proceedings, Cambridge, USA, July 5 8, 2010 489 Homogeneous Bianisotropic Medium, Dissipation and the Non-constancy of Speed of Light in Vacuum for Different
More information1. Reflection and Refraction of Spherical Waves
1. Reflection and Refraction of Spherical Waves Our previous book [1.1] was completely focused on the problem of plane and quasi-plane waves in layered media. In the theory of acoustic wave propagation,
More informationLecture 36 Date:
Lecture 36 Date: 5.04.04 Reflection of Plane Wave at Oblique Incidence (Snells Law, Brewster s Angle, Parallel Polarization, Perpendicular Polarization etc.) Introduction to RF/Microwave Introduction One
More informationarxiv: v2 [physics.optics] 10 Oct 2018
Dissipation Effect on Optical Force and Torque near Interfaces Daigo Oue 1, 1 Division of Frontier Materials Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka, Japan 560-8531 arxiv:1809.00445v2
More informationUniform Plane Waves Page 1. Uniform Plane Waves. 1 The Helmholtz Wave Equation
Uniform Plane Waves Page 1 Uniform Plane Waves 1 The Helmholtz Wave Equation Let s rewrite Maxwell s equations in terms of E and H exclusively. Let s assume the medium is lossless (σ = 0). Let s also assume
More informationFiber Optics. Equivalently θ < θ max = cos 1 (n 0 /n 1 ). This is geometrical optics. Needs λ a. Two kinds of fibers:
Waves can be guided not only by conductors, but by dielectrics. Fiber optics cable of silica has nr varying with radius. Simplest: core radius a with n = n 1, surrounded radius b with n = n 0 < n 1. Total
More informationScalar electromagnetic integral equations
Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how
More informationElectrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Electro Dynamic
Electrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Name Electro Dynamic Instructions: Use SI units. Short answers! No derivations here, just state your responses clearly. 1. (2) Write an
More informationWave Phenomena Physics 15c. Lecture 15 Reflection and Refraction
Wave Phenomena Physics 15c Lecture 15 Reflection and Refraction What We (OK, Brian) Did Last Time Discussed EM waves in vacuum and in matter Maxwell s equations Wave equation Plane waves E t = c E B t
More informationAnalysis of second-harmonic generation microscopy under refractive index mismatch
Vol 16 No 11, November 27 c 27 Chin. Phys. Soc. 19-1963/27/16(11/3285-5 Chinese Physics and IOP Publishing Ltd Analysis of second-harmonic generation microscopy under refractive index mismatch Wang Xiang-Hui(
More informationElectromagnetic Waves
May 7, 2008 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 7 Maxwell Equations In a region of space where there are no free sources (ρ = 0, J = 0), Maxwell s equations reduce to a simple
More informationCHAPTER 9 ELECTROMAGNETIC WAVES
CHAPTER 9 ELECTROMAGNETIC WAVES Outlines 1. Waves in one dimension 2. Electromagnetic Waves in Vacuum 3. Electromagnetic waves in Matter 4. Absorption and Dispersion 5. Guided Waves 2 Skip 9.1.1 and 9.1.2
More information(a) Show that the amplitudes of the reflected and transmitted waves, corrrect to first order
Problem 1. A conducting slab A plane polarized electromagnetic wave E = E I e ikz ωt is incident normally on a flat uniform sheet of an excellent conductor (σ ω) having thickness D. Assume that in space
More informationPHY3128 / PHYM203 (Electronics / Instrumentation) Transmission Lines
Transmission Lines Introduction A transmission line guides energy from one place to another. Optical fibres, waveguides, telephone lines and power cables are all electromagnetic transmission lines. are
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationUSPAS Accelerator Physics 2017 University of California, Davis
USPAS Accelerator Physics 2017 University of California, Davis Chapter 9: RF Cavities and RF Linear Accelerators Todd Satogata (Jefferson Lab) / satogata@jlab.org Randika Gamage (ODU) / bgama002@odu.edu
More informationSupplementary Information: Non-collinear interaction of photons with orbital angular momentum
Supplementary Information: Non-collinear interaction of photons with orbital angular momentum T. oger, J. F. Heitz, D. Faccio and E. M. Wright November 1, 013 This supplementary information provides the
More informationScattering of ECRF waves by edge density fluctuations and blobs
PSFC/JA-14-7 Scattering of ECRF waves by edge density fluctuations and blobs A. K. Ram and K. Hizanidis a June 2014 Plasma Science and Fusion Center, Massachusetts Institute of Technology Cambridge, MA
More informationANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. E = jωb. H = J + jωd. D = ρ (M3) B = 0 (M4) D = εe
ANTENNAS Vector and Scalar Potentials Maxwell's Equations E = jωb H = J + jωd D = ρ B = (M) (M) (M3) (M4) D = εe B= µh For a linear, homogeneous, isotropic medium µ and ε are contant. Since B =, there
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationSURFACE PLASMONS AND THEIR APPLICATIONS IN ELECTRO-OPTICAL DEVICES
SURFACE PLASMONS AND THEIR APPLICATIONS IN ELECTRO-OPTICAL DEVICES Igor Zozouleno Solid State Electronics Department of Science and Technology Linöping University Sweden igozo@itn.liu.se http://www.itn.liu.se/meso-phot
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationarxiv: v2 [physics.acc-ph] 27 Oct 2014
Maxwell s equations for magnets A. Wolski University of Liverpool, Liverpool, UK and the Cockcroft Institute, Daresbury, UK arxiv:1103.0713v2 [physics.acc-ph] 27 Oct 2014 Abstract Magnetostatic fields
More informationGUIDED MICROWAVES AND OPTICAL WAVES
Chapter 1 GUIDED MICROWAVES AND OPTICAL WAVES 1.1 Introduction In communication engineering, the carrier frequency has been steadily increasing for the obvious reason that a carrier wave with a higher
More informationIN conventional optical fibers, light confinement is achieved
428 JOURNAL OF LIGHTWAVE TECHNOLOGY, VOL. 20, NO. 3, MARCH 2002 Asymptotic Matrix Theory of Bragg Fibers Yong Xu, George X. Ouyang, Reginald K. Lee, Member, IEEE, and Amnon Yariv, Life Fellow, IEEE Abstract
More informationSpherical Waves. Daniel S. Weile. Department of Electrical and Computer Engineering University of Delaware. ELEG 648 Spherical Coordinates
Spherical Waves Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Spherical Coordinates Outline Wave Functions 1 Wave Functions Outline Wave Functions 1
More informationSurface Plasmon-polaritons on thin metal films - IMI (insulator-metal-insulator) structure -
Surface Plasmon-polaritons on thin metal films - IMI (insulator-metal-insulator) structure - Dielectric 3 Metal 2 Dielectric 1 References Surface plasmons in thin films, E.N. Economou, Phy. Rev. Vol.182,
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationLecture 1 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture 1 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell 1. Overview of the Course - Last semester we covered electrostatics, magnetostatics,
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationGeneral Appendix A Transmission Line Resonance due to Reflections (1-D Cavity Resonances)
A 1 General Appendix A Transmission Line Resonance due to Reflections (1-D Cavity Resonances) 1. Waves Propagating on a Transmission Line General A transmission line is a 1-dimensional medium which can
More informationMIDSUMMER EXAMINATIONS 2001
No. of Pages: 7 No. of Questions: 10 MIDSUMMER EXAMINATIONS 2001 Subject PHYSICS, PHYSICS WITH ASTROPHYSICS, PHYSICS WITH SPACE SCIENCE & TECHNOLOGY, PHYSICS WITH MEDICAL PHYSICS Title of Paper MODULE
More informationPhysics Letters A 374 (2010) Contents lists available at ScienceDirect. Physics Letters A.
Physics Letters A 374 (2010) 1063 1067 Contents lists available at ScienceDirect Physics Letters A www.elsevier.com/locate/pla Macroscopic far-field observation of the sub-wavelength near-field dipole
More informationTheory of Electromagnetic Fields
Theory of Electromagnetic Fields Andrzej Wolski University of Liverpool, and the Cockcroft Institute, UK Abstract We discuss the theory of electromagnetic fields, with an emphasis on aspects relevant to
More informationFRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES
Progress In Electromagnetics Research, PIER 5, 3 38, 000 FRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES Q. A. Naqvi and A. A. Rizvi Communications Lab. Department of Electronics Quaidi-i-Azam University
More informationIntroduction to optical waveguide modes
Chap. Introduction to optical waveguide modes PHILIPPE LALANNE (IOGS nd année) Chapter Introduction to optical waveguide modes The optical waveguide is the fundamental element that interconnects the various
More informationDepartment of Physics Preliminary Exam January 2 5, 2013
Department of Physics Preliminary Exam January 2 5, 2013 Day 2: Electricity, Magnetism and Optics Thursday, January 3, 2013 9:00 a.m. 12:00 p.m. Instructions: 1. Write the answer to each question on a
More informationElectromagnetic Scattering from a PEC Wedge Capped with Cylindrical Layers with Dielectric and Conductive Properties
0. OZTURK, ET AL., ELECTROMAGNETIC SCATTERING FROM A PEC WEDGE CAPPED WIT CYLINDRICAL LAYERS... Electromagnetic Scattering from a PEC Wedge Capped with Cylindrical Layers with Dielectric and Conductive
More informationNumerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material
Chapter 30 Numerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material Hiroshi Maeda Additional information is available at the end of the chapter http://dx.doi.org/10.5772/50555
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 9
WiSe 202 20.2.202 Prof. Dr. A-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg
More informationLECTURE 18: Horn Antennas (Rectangular horn antennas. Circular apertures.)
LCTUR 18: Horn Antennas (Rectangular horn antennas. Circular apertures.) 1 Rectangular Horn Antennas Horn antennas are popular in the microwave bands (above 1 GHz). Horns provide high gain, low VSWR (with
More information