Journal of Applied Mathematics

Volume 2011 (2011), Article ID 715087, 30 pages

http://dx.doi.org/10.1155/2011/715087

## Asymptotic Analysis of Transverse Magnetic Multiple Scattering by the Diffraction Grating of Penetrable Cylinders at Oblique Incidence

Department of Electrical and Computer Engineering, School of Engineering and Applied Science, The George Washington University, Washington, DC 20052, USA

Received 17 June 2011; Accepted 13 September 2011

Academic Editor: Yongkun Li

Copyright © 2011 Ömer Kavaklıoğlu and Roger Henry Lang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We have presented a derivation of the asymptotic equations for transverse magnetic multiple scattering coefficients of an infinite grating of penetrable circular cylinders for obliquely incident plane electromagnetic waves. We have first deducted an “*Ansatz*” delineating the asymptotic behavior of the transverse magnetic multiple scattering coefficients associated with the most generalized condition of oblique incidence (Kavaklıoğlu, 2000) by exploiting Schlömilch series corresponding to the special circumstance that the grating spacing is much smaller than the wavelength of the incident electromagnetic radiation. The validity of the asymptotic equations for the aforementioned scattering coefficients has been verified by collating them with the Twersky's asymptotic equations at normal incidence. Besides, we have deduced the consequences that the asymptotic forms of the equations at oblique incidence acquired in this paper reduce to Twersky's asymptotic forms at normal incidence by expanding the generalized scattering coefficients at oblique incidence into an asymptotic series as a function of the ratio of the cylinder radius to the grating spacing.

#### 1. Introduction

Rayleigh [1] first treated the problem of the incidence of plane electric waves on an insulating dielectric cylinder as long ago as 1881. He published the classical electromagnetic problem of the diffraction of a plane wave at normal incidence by a homogeneous dielectric cylinder [2]. His solution was later generalized for obliquely incident plane waves when the magnetic vector of the incident wave is transverse to the axis of the cylinder by Wait [3]. Moreover, Rayleigh [4, 5] adduced the first theoretical investigation for the problem of diffraction by gratings. His results have been extended by Wait [6] for the treatment of scattering of plane waves by parallel-wire grids with arbitrary angle of incidence. Wait [6] developed the solution of the problem of the scattering of plane electromagnetic waves incident upon a parallel-wire grid that was backed by a plane-conducting surface. He generalized this result subsequently to a plane wave, incident obliquely with arbitrary polarization on a planar grid [7]. Wait did not treat the scattering of obliquely incident plane electromagnetic waves by the infinite array of thick dielectric cylinders. This configuration has recently been studied by Kavaklıoğlu [8–10], and an analytic expression for the generalized multiple scattering coefficients of the infinite grating at oblique incidence was captured in the form of a convergent infinite series [11].

The formal analytical solution for the scattering of a plane acoustic or electromagnetic wave by an arbitrary configuration of parallel cylinders of different radii and physical parameters in terms of cylindrical wave functions was obtained by Twersky [12] who considered all possible contributions to the excitation of a particular cylinder by the radiation scattered by the remaining cylinders in the grating and extended this solution to expound the case where all the axes of cylinders lie in the same plane [13]. Twersky [14] subsequently introduced the formal multiple scattering solution of a plane wave by an arbitrary configuration of parallel cylinders to the finite grating of cylinders. He later employed Green’s function methods to represent the multiple scattering amplitude of one cylinder within the grating in terms of the functional equation and the single-scattering amplitude of an isolated cylinder [15]. Furthermore, Twersky [16] acquired a set of algebraic equations for the multiple scattering coefficients of the infinite grating in terms of the elementary function representations of Schlömilch series [17] and the well-known scattering coefficients of an isolated cylinder.

In the area of acoustics, Millar [18] studied the problem of scattering of a plane wave by finite number of cylinders equispaced in a row that are associated with scatterers both “soft” and “hard” in the acoustical sense. The solutions in the form of series in powers of a small parameter, essentially the ratio of cylinder dimension to wavelength, were obtained. Besides, Millar [19] investigated the scattering by an infinite grating of identical cylinders. In a more recent investigation, Linton and Thompson [20] formulated the diffracted acoustic field by an infinite periodic array of circles and determined the conditions for resonance by employing the expressions which enable Schlömilch series to be computed accurately and efficiently [21–23].

Previous investigations mentioned above do not include the most general case of oblique incidence although the grating is illuminated by an incident plane -polarized electromagnetic wave at an arbitrary angle to the -axis, whereas in the generalized oblique incidence solution presented in this investigation, the direction of the incident plane wave makes an arbitrary oblique angle with the positive -axis as indicated in Figure 1. As far as can be ascertained by the writers, Sivov [24, 25] first treated the diffraction by an infinite periodic array of perfectly conducting cylindrical columns for the most generalized case of obliquely incident plane-polarized electromagnetic waves in order to determine the reflection and transmission coefficients of the infinite grating of perfectly conducting cylinders in free space under the assumption that the period of the grating spacing was small compared to a wavelength. The configuration of a greater relevance to the problem has recently been investigated by many other researchers. For instance, Lee [26] studied the scattering of an obliquely incident electromagnetic wave by an arbitrary configuration of parallel, nonoverlapping infinite cylinders and acquired the solution for the scattering of an obliquely incident plane wave by a collection of closely spaced radially stratified parallel cylinders that can have an arbitrary number of stratified layers [27]. Moreover, Lee [28] presented a general treatment of scattering of arbitrarily polarized incident light by a collection of radially stratified circular cylinders at oblique incidence, described the solution to the problem of scattering of obliquely incident light by a closely spaced parallel radially stratified cylinders embedded in a semi-infinite dielectric medium [29], and developed a general scattering theory for obliquely incident plane-polarized monochromatic waves on a finite slab containing closely spaced radially stratified circular cylinders [30]. In addition, the formulation for the extinction and scattering cross-sections of closely spaced parallel infinite cylinders in a dielectric medium of finite thickness is presented [31]. In the area of modeling photonic crystal structures, Smith et al. [32] developed a formulation for cylinder gratings in conical incidence using a multipole method and studied scattering matrices and Bloch modes in order to investigate the photonic band gap properties of woodpile structures [33]. This area of research has recently received a lot of attention due to potential applications to microcircuitry, nanotechnology, and optical waveguides.

Three-dimensional generalization of Twersky’s solution [15, 16] for scattering of waves by the infinite grating of dielectric circular cylinders was originally developed by Kavaklıoğlu [8–10] by employing the separation-of-variables method for both TM and TE polarizations, and the reflected and transmitted fields were derived for obliquely incident plane -polarized waves in [34]. Kavaklıoğlu and Schneider [35] presented the asymptotic solution of the multiple scattering coefficients for obliquely incident and vertically polarized plane waves as a function of the ratio of the cylinder radius to grating spacing when the grating spacing, , is small compared to a wavelength.

Furthermore, Kavaklıoğlu and Schneider [11] acquired the exact analytical solution for the multiple scattering coefficients of the infinite grating for obliquely incident plane electromagnetic waves by the application of the direct Neumann iteration technique to two infinite sets of equations describing the exact behavior of the multiple scattering coefficients, which was originally published in [8, 10], in the form of a convergent infinite series and obtained the generalized form of Twersky’s functional equation for the infinite grating in matrix form for obliquely incident waves [11].

The purpose of this paper is to elucidate the derivation of the equations pertaining to the asymptotic behavior of the transverse magnetic multiple scattering coefficients of an infinite array of infinitely long circular dielectric cylinders illuminated by obliquely incident plane electromagnetic waves. The arbitrarily polarized obliquely incident plane wave depicted in Figure 1 can be decomposed into two different modes of polarization. The asymptotic representation associated with the transverse magnetic (TM) mode that is also defined as vertical polarization, for which the incident electric field has a component parallel to the constituent cylinders of the grating, will be treated in this investigation.

#### 2. Problem Formulation

“An infinite number of infinitely long identical dielectric circular cylinders,” which are separated by a distance “,” are placed parallel to each other in the - plane and positioned perpendicularly to the - plane as indicated in Figure 1. For TM mode; is the unit vector associated with the vertical polarization and has a component parallel to the cylinders of the grating. The fact that “the incident -field has a component parallel to all the cylinders of the dielectric grating” does not mean that we deal with the TM mode as it does not exclude the existence of other components of -field. The incident plane wave depicted in Figure 1 makes an angle of obliquity with the positive -axis.

Lemma 2.1 (multiple scattering representation for an infinite grating of dielectric circular cylinders for obliquely incident -polarized plane electromagnetic waves [3, 7, 8]). *A vertically polarized plane electromagnetic wave, which is obliquely incident upon the infinite array of identical insulating dielectric circular cylinders with radius “,” dielectric constant “,” and relative permeability “,” can be expanded in “the individual cylindrical coordinate system of the th cylinder” in terms of the cylindrical waves referred to the axis of th cylinder as
*

The origin of each individual cylindrical coordinate system, namely, , is located at the center of the corresponding cylinder. In the above description of the incident field, is a unit vector that denotes the vertical polarization having a component parallel to all the cylinders, is the angle of incidence in - plane measured from the -axis in such a way that , implying that the wave is obliquely incident in the first quadrant of the coordinate system, and “” stands for a Bessel function of order *n*. In addition, we have the following definitions:
“” time dependence is suppressed throughout the paper, where “” stands for the angular frequency of the incident wave in radians per second, “” is the free-space wave number, “” denotes the speed of light in free space, and “” represents time in seconds. The centers of the cylinders in the infinite grating are located at the positions , and so forth. The exact solution for the -component of the electric field in the exterior of the grating belonging to this configuration can be expressed in terms of the incident electric field in the coordinate system of the th cylinder located at , plus a summation of cylindrical waves outgoing from each individual th cylinder located at , as , that is,

Lemma 2.2 (expressions for the -components of the exterior fields [8]). *Let for all , where “” stands for the set of all integers, denote the set of all multiple scattering coefficients corresponding to the exterior electric and magnetic fields of the infinite grating associated with obliquely incident plane -polarized electromagnetic waves, respectively. Then, the exterior electric and magnetic field intensities associated with vertically polarized obliquely incident plane electromagnetic waves are given as *

In this representation, depicts the set of all undetermined multiple scattering coefficients associated with exterior electric fields defined by the expressions (29) and (34)–(37) in [8], and delineates the set of all undetermined multiple scattering coefficients associated with exterior magnetic fields defined by the expressions (40)–(42) in [8], respectively. In expressions (2.4a) and (2.4b), we have where and “” denotes the th order Hankel function of first kind, for all . The series in expression (2.4b) is the generalization of the “*Schlömilch series* for obliquely incident electromagnetic waves” [10, 17] and converges provided that does not equal integers. The integral values of are known as the “grazing modes” or “Rayleigh values” [17]. The convergence of the series for the scattering coefficients can be found on page 342 in [11]. Moreover, the convergence of the Schlömilch series has been discussed by Twersky [17] in detail, who also gives additional references. The exact expressions corresponding to the radial and angular components of the electric and magnetic field intensities have already been obtained in [8] by employing the -component of the external field in the expressions (2.4a) and (2.4b).

#### 3. Derivation of the Asymptotic Equations for the Multiple Scattering Coefficients of the Infinite Grating at Oblique Incidence

This section is devoted to the formal derivation of the asymptotic equations for the exterior electric and magnetic multiple scattering coefficients of the infinite grating of dielectric cylinders for obliquely incident vertically polarized plane waves. Since the wavelength of the incident radiation is much larger than the grating spacing, the condition is automatically satisfied thereby excluding any special case associated with the grazing modes. In order to demonstrate the procedure of obtaining the asymptotic equations for the TM multiple scattering coefficients of the infinite grating at oblique incidence, we will first introduce the exact equations corresponding to the transverse magnetic multiple scattering coefficients associated with the exterior electric and magnetic fields of the infinite grating of dielectric circular cylinders at oblique incidence by asserting the following lemma.

Lemma 3.1 (exact equations of the transverse magnetic multiple scattering coefficients of the infinite grating of insulating dielectric cylinders at oblique incidence [8]). *Exact equations corresponding to the transverse magnetic multiple scattering coefficients of an infinite grating of insulating dielectric cylinders associated with obliquely incident plane electromagnetic waves are first presented by the equations (85a) and (85b) in [8] as
**
The coefficients arising in this infinite set of linear algebraic equations are defined as
**
Two sets of constants and , in which stands for the relative permittivity and permeability of the dielectric cylinders, respectively, are given as
**
for , and for all ; where is defined as , and
**
for and for all , where in the expression above is a constant and given as
**
In these equations and denote the relative dielectric constant and the relative permeability of the insulating dielectric cylinders; and stand for the permittivity and permeability of the free space, respectively. In addition, , and in expressions (3.2) and (3.3) are defined as
**
which imply the first derivatives of the Bessel and Hankel functions of first kind and of order with respect to their arguments.*

Theorem 3.2 (approximate equations for the scattering coefficients of the infinite grating at oblique incidence when ). *The asymptotic form of the exact equations for the transverse magnetic multiple scattering coefficients of an infinite grating at oblique incidence can be inferred by two different sets, in which the first one contains only the odd coefficients and the second set contains only the even coefficients. Odd multiple scattering coefficients associated with the infinite grating of dielectric circular cylinders at oblique incidence satisfy the following two sets of asymptotic equations:
**
Similarly, the even multiple scattering coefficients satisfy the following two infinite sets of asymptotic equations associated with the transverse magnetic multiple scattering coefficients of the infinite grating of dielectric circular cylinders at oblique incidence as
**
where denotes the set of all natural numbers.*

*Proof. *The exact equations in (3.1) can be solved for , and when the distance between the cylinders of the infinite grating is smaller than the wavelength of the incident wave, that is, for the exact equations take the following form:
where is a matrix defined as
and “” connotes the approximation to the “exact form of the Schlömilch series ” in the limiting case when for . Introducing (3.10)) into (3.9), the approximate set of equations for the scattering coefficients of the infinite grating at oblique incidence can explicitly be written as
In the above, we have
The -dependent constants appearing in (3.10) and (3.11) are defined asThe various constants appearing in the definitions (3.13a)–(3.13d) are expressed as
The elements of the matrix of coefficients in ((3.11) can be calculated using the expressions (3.14), for instance and terms can be written as
In terms of the definitions of (3.13a)–(3.13d), the approximate set of equations for the multiple scattering coefficients of the infinite grating at oblique incidence given in ((3.11) takes the following form:
Statement of Theorem 3.2 follows immediately upon decomposition of (3.16) into its odd and even components as it is designated by (3.7) and (3.8).

The elementary function representations of the *Schlömilch series * in (2.5b) have originally been derived by Twersky [17] for the normal incidence and modified by Kavaklıoğlu [10] for the oblique incidence. We will employ these elementary function representations for the evaluation of the asymptotic forms of the Schlömilch series in the limit of *. *Twersky’s forms [16, 17] are still valid for the case of obliquely incident waves [10] with a slight modification in their arguments.

Lemma 3.3 (approximate expressions for the “Schlömilch series ” in the limit of [10, 17]). *We have obtained for the special case of as
**
where . In (3.17), is defined by the following relationship:
**
The angles are the usual “diffraction angles” of the grating, and (3.18) that provides these discrete angles is called the “grating equation”. “Propagating modes” are determined by , and they correspond to , the ’s being the closest integers to the ’s for which should be satisfied, that is, , such that
**
“Evanescent modes” are determined by , and they correspond to integer values of such that , we have , and are determined by . For this case the “grating equation” takes the form of
**
For the general case, we have , for all as
**
Finally, is the Bernoulli polynomial of argument “” and power “”, in ((3.21).*

*Remark 3.4 (Bessel series , , and ). *The propagating range of the *Schlömilch series *, for all , where in (3.17) and ((3.21), is described by “” *Bessel series, *which can explicitly be written as

*Remark 3.5 (Neumann series , , and ). *The evanescent range of the *Schlömilch series *, for all in (3.17) and (3.21), is described by “” where is known as the *Neumann series. * in (3.17) and ((3.21) can be put into the following form for this limiting case () as
In addition, we can obtain the simplified expressions for and as

*Remark 3.6 (special case when ). *The physical problem under consideration corresponds to the special case for which there is only one propagating mode and the scattering of wavelengths is larger than the grating spacing, that is, . Then the Bessel series for , which implies that the plane wave is incident onto the grating in the first quadrant, for all , reduces to
where stands for the *Kronecker delta function*.

*Remark 3.7 (approximations for Neumann series, , and in the limit of ). *Inserting in (3.23a), (3.23b), and (3.23c), the expression for in (3.23a)) reduces to
The approximation of the *Neumann series **,* for , up to terms of the order can be obtained from (3.25a) as
where , for all , denotes the Riemann zeta function. In the same range, the *Neumann series * reduces to
where ’s in (3.26) are given as

*Remark 3.8 (approximations for Schlömilch series, in the limit of ). *If is small, that is to say if , then there is only one discrete propagating mode. Employing the expansions for the *Bessel and Neumann Series* obtained in the previous sections for ; the *Schlömilch Series* in this range can be expressed as

*Remark 3.9 (leading terms of the Schlömilch series, in the limit of ). *The leading terms of ’s for large “”, for all is given as
where corresponds to the* Bernoulli Polynomial*. From (3.28) and (3.29), we can determine the leading terms of the *Schlömilch series* as
The leading terms of for large “” are given by
where ’s and ’s for large "" are given as
In the above expressions, ’s are the *Bernoulli numbers,* and the relationship between *Bernoulli polynomial* and *Bernoulli numbers *is given as

#### 4. Asymptotic Expansions for the Scattering Coefficients of the Infinite Grating at Oblique Incidence in the Limiting Case of “”

In order to find a solution for the set of equations given in (3.15) and (3.16), we have introduced an *“Ansatz”* [36] for the scattering coefficients of the electric and magnetic fields of the infinite grating assuming , and , as
for all , for the odd multiple coefficients corresponding to the electric and magnetic field intensities of the infinite grating associated with obliquely incident plane electromagnetic waves, andfor all for the even multiple coefficients. In the above expressions, we have delineated the wavelength-independent parts of the multiple scattering coefficients associated with the exterior electric and magnetic field intensities as .

Theorem 4.1 (asymptotic equations for the multiple scattering coefficients corresponding to the exterior electric and magnetic field intensities associated with obliquely incident vertically polarized plane electromagnetic waves). *The multiple scattering coefficients corresponding to the exterior electric and magnetic field intensities associated with obliquely incident vertically polarized plane electromagnetic waves satisfy two infinite sets of asymptotic equations described by**
for the odd multiple scattering coefficients, and
**for the even multiple scattering coefficients.*

*Proof. *We have defined the overall effect of the multiple scattering terms when the wavelength is much larger than the grating spacing, that is, , and as
for the electric field coefficients, and
for the magnetic field coefficients. Employing *the approximations of Schlömilch series given in (3.27)* in the expressions (4.3a) and (4.3b), we can write the overall effect of the multiple scattering terms when the wavelength is much larger than the grating spacing, that is, , and as
for all , for the odd coefficients,
for all , for the even coefficients; and the special case for is given by
Defining the wavelength independent parts of the scattering matrices from (3.10)) as
for all , corresponding to the odd, and
corresponding to the even part. Using the definitions in (4.6a)–(4.6d), and introducing (4.4a)–(4.4d) into (3.7) and (3.8), we have obtained the following set of equations for the approximations of the scattering coefficients:
corresponding to the odd scattering coefficients, and
corresponding to the even scattering coefficients. Splitting the matrices in (4.7a) and (4.7b) into two parts, we have
for the odd scattering coefficients, and
for even scattering coefficients. From (4.4a)–(4.4d), we have established the following terms:
for the multiple interactions corresponding to the scattering coefficients of the electric field, and
for the multiple interactions corresponding to the scattering coefficients of the magnetic field, for all , for the odd scattering coefficients, and
for the multiple interactions corresponding to the scattering coefficients of the electric field,
for the multiple interactions corresponding to the scattering coefficients of the magnetic field, for all , for the even scattering coefficients. Inserting (4.9a)-(4.9b) and (4.10a)-(4.10b) into (4.8a)-(4.8b), we have finally obtained the infinite set of asymptotic equations for the multiple scattering coefficients corresponding to the exterior electric and magnetic field intensities of an infinite grating of dielectric circular cylinders associated with obliquely incident and vertically polarized electromagnetic waves as it is proposed by the statement of Theorem 4.1 introduced in (4.2a) and (4.2b). In addition, we have noticed that the scattering coefficients of the electric and magnetic fields appeared as coupled to each others.

#### 5. Discussion and Comparison of the Generalized Transverse Magnetic Multiple Scattering Coefficients of the Infinite Grating with Twersky’s Normal Incidence Case

*Remark 5.1 (Twersky’s asymptotic solution for the multiple scattering coefficients at normal incidence). *The exact equations for the multiple scattering coefficients of the infinite grating associated with the vertically polarized normally incident waves [16] can be solved by truncation as
where, the numerator terms are given as
and the denominator terms are given as
The ’s, for all , are given by
Finally, ’s appearing in (5.2) represent the asymptotic forms of the single-scattering coefficients associated with an isolated cylinder within the grating at normal incidence and can be approximated for as
for all ,wherefor all . Inserting (5.3), (3.29), and (3.30) into (5.2), ’s can be evaluated asExpressions in (5.5a)–(5.5d) are valid for , and . In general, ’s can be expressed as
Obviously, (5.5a) and (5.6) will asymptotically be written as
for all, where represent -independent parts of ’s for wavelengths larger than the radii, that is, . The numerator terms appearing in (5.1b)–(5.1e) can be approximated as
The denominator terms appearing in (5.1c) can be approximated as
Inserting (5.5a)–(5.5d) to (5.10) into (5.1a)–(5.1g), we have