Diffraction gratings - davidar/scholarpedia GitHub Wiki
Born in the eighteenth century (Rittenhouse, 1786), the diffraction grating has been one of the most valuable instruments in the history of science and technology. At one extreme, it has enabled the study of celestial bodies. On the other hand, it has been a crucial tool in the study of atomic and molecular structures. These outstanding capabilities are the consequence of one basic property: it can disperse light, i.e. separate the frequencies contained in light radiation, thus allowing the measurement of relative intensities.
Diffraction grating, notations
Figure 1 represents a diffraction grating. The periodic profile <math>\mathcal{P}</math> of period d along the x-axis separates air from a grating material which is generally a metal or a dielectric. The y-axis is perpendicular to the average profile plane and the z-axis is the axis of invariance of the structure. We suppose that the incident light can be described by a sum of monochromatic radiations of different frequencies. Each of these can in turn be described in a time-harmonic regime, which allows us to use the complex notation (with an <math>\exp(-i\omega t)</math> time-dependence). The electromagnetic properties of the grating material (assumed to be non-magnetic) are represented by its complex refractive index <math>\nu</math> which depends on the wavelength <math>\lambda=2\pi c/\omega</math> in vacuum (<math>c=1/\sqrt{\varepsilon_0\mu_0}</math> being the speed of light, with <math>\varepsilon_0</math> and <math>\mu_0</math> the permittivity and the permeability of vacuum). This complex index respectively includes the conductivity (for metals) and/or the losses (for lossy dielectrics). It becomes a real number for lossless dielectrics.
In the air region, the grating is illuminated by an incident plane wave, assumed to be s-polarized (electric field parallel to the z-axis) for simplicity, but the analysis is similar in the other case (p-polarization). With these conventions, the incident electric field <math>\overrightarrow{E^i} =E^i\hat z</math> (<math>\hat z</math> unit vector of the z-axis) is given by :
<math>\label{eq1} E^i=\exp\bigl( ik_0 x\sin(\theta)-ik_0 y\cos(\theta)\bigr) , </math>
with <math>\theta</math> being the angle of incidence, from the y-axis to the incident direction, measured in the counterclockwise sense, and <math>k_0 = 2\pi/\lambda</math> (we take an index equal to unity for air).
Let us now introduce, without demonstration (Petit, 1980), three intuitive properties of the total electric field. Like the incident electric field <math>\overrightarrow{E^i}</math> , the total field <math>\overrightarrow{E}</math> is:
- invariant along the z-axis,
- parallel to the z-axis, thus <math>\overrightarrow{E}=E(x,y)\hat z</math> .
- pseudo-periodic , i.e. by definition :
Notice that in normal incidence (<math>\theta=0</math> ), pseudo-periodicity becomes ordinary periodicity, which in that case is a straightforward property since both grating and incident wave are periodic.
Using these three properties and Maxwell's equations, let us show that the field above the grating is a sum of plane waves. With this aim, we use the first two Maxwell equations:
<math>\label{eq3} \nabla\times\overrightarrow{E}=i\omega\mu_0\overrightarrow{H},\qquad \nabla\times\overrightarrow{H}=-i\omega\varepsilon\overrightarrow{E}, </math>
with <math>\overrightarrow{H}</math> being the magnetic field and <math>\varepsilon</math> being the permittivity, equal to <math>\varepsilon_0</math> in air region and to <math>\varepsilon_0 \nu^2</math> in the grating material. In the following, equations placed at the left and right-hand sides of \eqref{eq3} will be called first and second Maxwell equations respectively. Combining the two Maxwell equations the electric field satisfies in the entire space a Helmholtz equation:
<math>\label{eq4} \nabla^2E+k^2E=0, </math>
with: <math> k=\begin{cases}k_0 &\text{in the air} \\k_0 \nu &\text{in the material} \end{cases} </math>
Moreover, it can be deduced from equation \eqref{eq2} that <math>E(x,y)\exp\bigl( -ik_0 x\sin(\theta)\bigr) </math> has a period d along the x-axis and can thus be expanded in a Fourier series:
<math>\label{eq5} E(x,y)\exp\bigl( -ik_0 x\sin(\theta)\bigr) =\sum_{n=-\infty}^{+\infty} {E_n (y)\exp(2i\pi nx/d)}. </math>
Multiplying both members of equation \eqref{eq5} by <math>\exp\bigl(ik_0 x\sin(\theta)\bigr)</math> yields :
<math>\label{eq6} E(x,y)=\sum_{n=-\infty}^{+\infty} {E_n (y)\exp(i\alpha_n x)},\qquad \alpha_n=k_0 \sin(\theta)+2\pi n/d. </math>
Introducing this expression of <math>E(x,y)</math> in equation \eqref{eq4}, we find :
<math>\label{eq7} \sum_{n=-\infty}^{+\infty}{\bigl( \mathrm{d^2}E_n (y)/\mathrm{d}y^2+\beta_n^2 E_n (y)\bigr) \exp(i\alpha_n x ) }=0, \qquad \beta_n=\sqrt{(k_0^2-\alpha_n^2)}. </math>
Multiplying both members of the above equation by <math>\exp\bigl (-ik_0 x\sin(\theta)\bigr) </math>, we obtain a Fourier series in the left-hand side. It is vital to bear in mind that the nullity of this Fourier series entails the nullity of its coefficients, provided that the Fourier series vanishes for any value of x. Thus, equation \eqref{eq7} must be satisfied for any value of x. Since this equation is satisfied in the air region only, <math>y</math> must be greater than the ordinate <math>y_M</math> of the top of the grooves. In other words, we must restrict the property to the region located above the dotted blue line of figure 1. Finally, in this region:
<math>\label{eq8} \text{if } y>y_M,\qquad E_n(y)=I_n\exp(-i\beta_n y) + R_n\exp(+i\beta_n y), </math>
and therefore, using equation \eqref{eq6}, <math>\label{eq9} \text{if } y>y_M,\qquad E(x,y)=\sum_{n=-\infty}^{+\infty}{\bigl( I_n\exp(i\alpha_n x-i\beta_n y) + R_n\exp(i\alpha_n x+i\beta_n y)\bigr)}. </math>
Let us remark that equation \eqref{eq7} does not assign to <math>\beta_n</math> a unique value. However, equation \eqref{eq9} shows that its determination can be chosen arbitrarily since a sign change does not modify the value of the field, provided that <math>I_n</math> and <math>R_n</math> are permuted. The constants <math>\beta_n=\sqrt{(k_0^2-\alpha_n^2)}</math> being either real or purely imaginary, we impose the real or imaginary part to be positive. Equation \eqref{eq9} shows that the field above the grooves can be represented by a plane wave expansion. The propagation constants of the plane waves along the x- and y-axes are respectively equal to <math>\alpha_n</math> and <math> \pm \beta_n</math>. Some of these plane waves must be rejected since they do not obey the radiation condition (Petit, 1980). This rather intuitive condition states that the scattered field <math>E^s</math> defined in air region by <math>E^s=E-E^i </math> cannot include plane waves propagating downwards or exponentially increasing at infinity. Assuming that the incident electric field has a unit amplitude, this condition entails that <math> I_n=\delta_{n,0} </math> with <math> \delta_{n,0} </math> denoting the Kronecker symbol. Finally, equation \eqref{eq9} becomes:
<math>\label{eq10} \text{if } y>y_M,\qquad E(x,y)=\exp(i\alpha_0 x-i\beta_0 y) +\sum_{n=-\infty}^{+\infty}{ R_n\exp(i\alpha_n x+i\beta_n y)}, </math>
the sum being the expression of the scattered field in air region. The unknown complex coefficients <math> R_n </math> are the amplitudes of the reflected waves.
The conclusion of the last section is that above the grooves, the field scattered by the grating takes the form of a sum of plane waves, each of them being characterized by its order n. According to equation \eqref{eq7}, almost all these waves (an infinite number) are evanescent and propagate along the x-axis at the vicinity of the grating profile. They correspond to the orders n such that <math> \alpha_n^2 \ge k_0^2 </math>, thus rendering <math> \beta_n=i\sqrt{(\alpha_n^2-k_0^2)} </math> a purely imaginary number. Only a finite number of them, called y-propagative orders, propagate towards <math> y=+\infty </math> , with <math> \alpha_n^2 \le k_0^2 </math> , thus <math> \beta_n=\sqrt{(k_0^2-\alpha_n^2)}</math> being real. Let us notice that among these orders, the <math>0^{th}</math> order always exists. It propagates in the direction specularly reflected by the mean plane of the profile, whatever the wavelength may be. In contrast, the other y-propagative orders are dispersive. Indeed, their propagation constants along the x- and y-axes are equal to <math> \alpha_n </math> and <math> \beta_n </math>, in such a way that the diffraction angle <math> \theta_n </math> of one of these waves, measured clockwise from the y-axis, can be deduced from <math> \alpha_n = k_0 \sin(\theta_n) </math>. Using the expression of <math> \alpha_n </math> given by equation \eqref{eq6}, the angle of diffraction is given by :
<math>\label{eq11} \sin(\theta_n) = \sin(\theta)+n\displaystyle\frac{2\pi }{k_0 d}= \sin(\theta)+n\displaystyle\frac {\lambda}{d} . </math>
This is the famous grating formula, often deduced from heuristic arguments of physical optics.
The field below the grooves has the same characteristics: as the field above the grooves, it can be represented by a sum of plane waves. However, it does not contain incident waves and, due to the radiation condition, the plane waves must propagate towards <math> y=-\infty </math> or be attenuated in that direction. The wavenumber <math> k_0 </math> is replaced by <math> k_0 \nu </math> , which yields:
<math>\label{eq12} \text{if } y<0,\qquad E(x,y)=\sum_{n=-\infty}^{+\infty}{T_n\exp(i\alpha_n x-i\gamma_n y)} \qquad \gamma_n=\sqrt{(k_0^2 \nu^2 -\alpha_n^2)}. </math>
It is important to notice that this equation is obtained by assuming that the constant propagations <math> \alpha_n </math> inside the material are the same as those in the air. This property is the direct consequence of pseudo-periodicity (equation \eqref{eq2}) which is valid in the entire space.
If the grating material is a lossless dielectric, the directions of propagation of the transmitted field obey a grating formula as well. This formula is similar to equation \eqref{eq11} , but since the constants of propagation <math>\alpha_n</math> remain the same as in the air and since the wavenumber <math> k_0 </math> is multiplied by <math> \nu </math>, the angles of transmission <math>\theta'_n</math> can be deduced from <math>\alpha_n = k_0 \nu \sin(\theta'_n) </math>, which yields:
<math>\label{eq13} \nu \sin(\theta'_n) = \sin(\theta)+n\displaystyle\frac{2\pi }{k_0 d}= \sin(\theta)+n\displaystyle\frac {\lambda}{d} . </math>
When the grating material is metallic, the transmitted plane waves are absorbed by the metal and the y-propagating orders below the grooves no longer exist anymore.
In conclusion of this section, the reflected and transmitted fields include, outside the grooves, a finite number of plane waves propagating to infinity with scattering angles given by the grating formulae. All the orders are dispersive, except the <math>0^{th}</math> orders. The reflected <math>0^{th}</math> order takes the specular direction while for a lossless material, the transmitted <math>0^{th}</math> order takes the direction transmitted by a plane interface. Consequently, a polychromatic incident plane wave generates in a given order a sum of plane waves scattered in different directions, i.e. a spectrum. The measurement of the intensity along this spectrum allows one to determine the spectral power of the incident wave, provided that the grating has been calibrated. The calibration, which can be achieved experimentally (using a known polychromatic incident plane wave) or theoretically (from the grating shape and the refractive index) provides the intensity ratio between the scattered wave in the spectral order and the incident wave, versus wavelength. This ratio is called efficiency of the grating in the given order and is smaller than unity. Of course, precise measurements of the power spectrum require large efficiencies inside the spectrum range. This requirement, which is crucial in spectroscopy, is not easy to obtain a priori. Indeed, we have to bear in mind that the efficiency of a dispersive order is always challenged by the <math>0^{th}</math> order at least. This remark reveals the vital importance of the electromagnetic theory of gratings, which allows one to calculate the efficiency curve in a given order, and thus providing vital information about the spectroscopic quality of the grating before fabrication. Of course, this grating must be constructed using available technologies, a fact which points out that the progress of grating fabrication is vital.
For a long time, diffraction gratings have been considered as dispersive components of optics and used as tools for spectroscopic purposes only. The dispersive properties of gratings play a key role in many optical instruments, the most important being spectrometers and monochromators. The constant deviation mounting. The entrance and exit slits are fixed and the grating rotates around an axis perpendicular to the figure In these instruments, the dispersion of light is not obtained using a fixed angle of incidence mounting as in figure 1. Indeed, for this configuration, the light source and the grating are fixed, therefore power measurements of the spectrum require a rotation of the detector. It is much easier to keep the source and the detector fixed and to rotate the grating, as illustrated in figure 2. The polychromatic light source is placed on the entrance slit, which is located on the focus of the lens L1. the grating is illuminated under the variable incidence <math>\theta</math> by a parallel beam and generates spectra corresponding to the various propagative orders with <math>n \neq 0</math> . The light which goes to the exit slit contains a set of wavelengths corresponding to different orders n, each of these wavelengths satisfying the grating formula (equation \eqref{eq11}) with angle of diffraction <math>\theta_n=-D-\theta.</math> With our notation, D is the fixed angle between the axes of lenses <math>L_1</math> and <math>L_2</math>, measured in the counterclockwise sense. It can be easily shown that, using the spectrum in the <math>-1^{st}</math> order, a single wavelength is selected, provided that the range of wavelengths in the polychromatic light is not too large (Hutley, 1982). This wavelength is given by:
<math>\label{eq14} \sin(\theta_{-1}) = \sin(\theta)-\displaystyle\frac {\lambda}{d} . </math>
Thus, the wavelength is linked to the angle of incidence by the equation :
<math>\label{eq15} \lambda = d\bigl (\sin(\theta)+\sin(D+\theta) \bigr). </math>
In general, the deviation D is small and this constant deviation mount is identified as the Littrow mount, in which the deviation D vanishes. In this mounting, the incident wave and the <math>-1^{st}</math> order have opposite directions and the wavelength is given from the incidence angle by :
<math>\label{eq16} \lambda = 2d\sin(\theta). </math>
It can be verified from the grating formula that the <math>0^{th}</math> and <math>-1^{st}</math> orders are the only y-propagative orders in Littrow mount as long as <math>\lambda>2d/3</math> (i.e. equivalently <math>\theta>19.5°</math> ). This characteristic holds until the maximum wavelength <math> \lambda=2d </math> (or equivalently <math>\theta=90°</math> ). This range of wavelength is especially interesting for a spectroscopist using a metallic grating since the scattered power is shared between the <math>-1^{st}</math> and <math>0^{th}</math> orders only. As a consequence, the so-called blaze effect may be expected, at least for a few wavelengths in the spectrum. Blaze effect for an echelette grating The blaze effect is obtained when the <math>-1^{st}</math> order is the only reflected order, the intensity in the <math>0^{th}</math> order being null. A simple means to get the blaze effect is to use an echelette grating, as shown in figure 3 (Maréchal and Stroke, 1959). This kind of grating was for a long time the only manufactured grating. Its profile is made of two plane orthogonal facets, the large one making an angle b (blaze angle) with the mean plane of the grating. In Littrow mounting, when the angle of incidence is equal to the blaze angle, we can expect a very strong intensity in the <math>-1^{st}</math> order since its direction of propagation is that of the specular reflection on the large facet. This heuristic blaze property can be demonstrated rigorously if the metal is considered as perfectly conducting (i.e. with an infinite conductivity) for p-polarization (magnetic field parallel to the grooves). In that case, there is no loss in the metal and thus, when the blaze effect arises, the totality of the incident power goes to the <math>-1^{st}</math> order. It should be noticed that this blaze property holds, even when the number of reflected orders is greater than 2. For s-polarization (electric field parallel to the grooves), the efficiency is generally close (but not equal) to unity. It must be emphasized that the blaze effect can be obtained with other kinds of gratings, for example holographic gratings (with profiles close to a sinusoid) or lamellar gratings (with rectangular grooves), provided that the number of y-propagative orders is equal to 2 (Breidne and Maystre, 1980). When this is not the case, the echelette grating is, by far, the best one for blaze properties in Littrow mount.
In a spectrometer, the mount of figure 2 is used in order to know the spectral density of an unknown source, which can belong to a range extending from the ultraviolet to the infrared. For example, spectrographs are associated with terrestrial or embarked telescopes. They allow the astronomers to obtain crucial information on celestial objects like composition or speed. They also are of vital importance in the study of atomic or molecular structures via the light they emit. The Littrow mount is also used at the end of die laser cavities: a rotation of the grating permits the variation of the wavelength emitted by the laser (Hänsch, 1972).
X-ray spectroscopy is used in the study of synchrotron radiation. In that case, the period of the grating is much larger than the wavelength and, in general, a spectrum from large or very large orders must be used. Thus, in order to obtain significative efficiencies, it is necessary to employ mountings with specular reflection on the facet, i.e. mountings close to that of figure 3, with the direction of the employed order deduced from the incident direction by specular reflection on a facet, but with the incidence on the facet different from 0. Furthermore, due to the small reflectivity of metals under normal incidence in this domain, a large incidence angle on the large facet is used in general. Another means for enhancing the efficiency is to cover the large facet with alternate thin films of metallic and dielectric materials. More generally, the use of gratings with specular reflection on a facet are used in the so-called echelle gratings, viz. echelette gratings with large periods (Loewen and Popov, 1997). These gratings are often used in a large order, which increases the resolution. The resolution is the ability of a grating to separate two wavelengths, taking into account that a grating always has a finite size, thus that a monochromatic incident plane wave does not generate plane waves in the different orders but rather scattering lobes (Hutley, 1980).
Among the non-spectroscopic applications of diffraction gratings, one of the oldest one is the use as a light polarizer in the infrared. A metallic grid composed of a periodic set of circular metallic rods is a very good reflector for s-polarized light above a cut-off wavelength, the sharpness of the cut-off increasing with rod thickness . One can have a heuristic explanation of this phenomenon by considering that the space located between two rods behaves like a truncated two-dimensional waveguide with variable section. Such a waveguide has a cut-off wavelength for s-polarization, above which it cannot transmit light, except by tunneling effect, which explains why the grid is a polarizer. The same high-pass filtering property holds for unpolarized light with inductive grids, i.e. biperiodic set of holes perforated in a metallic screen. This property may be used for example for solar absorption (McPhedran and Maystre, 1977). The inductive grid is placed above an absorbing body. Such a device can absorb solar energy if the cut-off wavelength is of the order of 1500 nm. Indeed, nearly all the solar power (mainly confined in the range 350-1500 nm) is transmitted to the absorbing body and, if this body is cooled by adequate liquid pipes, the thermic radiation it emits in the infrared above 1500 nm cannot escape from the absorber, thanks to the inductive grid.
A large part of the more recent applications of diffraction gratings is based on their ability to provoke strong resonance phenomena. One of the most striking examples of this property is the grating coupler, which can transmit the power of a light beam inside a dielectric waveguide (Neviere, Petit and Cadilhac, 1973). Grating coupler The waveguide (figure 4) is formed by a dielectric film of index <math> \nu_g </math>(blue) placed above a dielectric substrate (orange) of index <math> \nu_s < \nu_g </math> and covered by a very thin dielectric film (green) on a part of which has been made the grating profile. If the grating is removed, this device is a planar waveguide and it can guide light by successive total reflections on both sides of the central dielectric film (blue). There is no loss in the propagation, but consequently the power of an incident light beam cannot be transmitted inside the waveguide. Indeed, the propagation constant <math> k_g </math> of a guided wave is greater than the wavenumber <math> k_0 </math> of the light in air region: as a consequence, the propagation constant along the x-axis of an incident plane wave or a beam, which it is smaller than <math> k_0 </math>, cannot match <math> k_g </math>. We have seen that a grating can scatter evanescent waves having propagation constants <math> \alpha_n </math> greater than <math> k_0 </math> (see equation \eqref{eq6}). Let us suppose that the entire upper profile of the waveguide is modulated. If, in this modulated waveguide, it exists an evanescent order n such that <math> \alpha_n = k_g </math>, the incident light can enter the waveguide through this evanescent order by resonant excitation. Conversely, the guided light can escape from the modulated waveguide and is thus exponentially attenuated in the propagation. However, if the modulated region has a finite width like in figure 4, the part of the incident power which is transmitted in the non-modulated region of the waveguide can propagate without any loss. This property is crucial in the technology of optical communications. Furthermore, due to the resonant process, the orders reflected and transmitted by the modulated waveguide present violent variations of efficiency and phase. Thus, if the planar waveguide is chosen in such a way that it transmits a large part of a polychromatic light, the modulated waveguide may present an efficiency very close to zero in the <math>0^{th}</math> reflected order (like the planar waveguide), except at the vicinity of the resonant wavelength where it can reach unity (Popov et al., 1986). This phenomenon has been used for DWDM (Dense Wavelength Division Multiplexing) in optical communications (Boyko et al., 2009).
Another vital phenomenon of resonance, very close to that described in figure 4, is obtained by replacing the dielectric waveguide by a metal. Indeed, a plane interface between air and metal regions can guide a surface plasmon propagating along the interface with a propagation constant <math> k_p > k_0 </math> (Raether, 1988). In contrast with the case of the dielectric waveguide, the surface plasmon is attenuated during propagation, due to the losses inside the metal. If the surface is modulated in order to obtain a metallic grating, an incident beam may generate the surface plasmon by resonant excitation, provided that the propagation constant of an evanescent wave matches the propagation constant of the surface plasmon. As a consequence, the amplitude of the surface plasmon becomes very large and violent phenomena happen. First, a large part of the incident power may be absorbed by Joule effect. Secondly, the efficiencies vary very rapidly at the vicinity of the resonance when the wavelength or the incidence angle are modified. This phenomenon, classified as Wood’s anomaly, was discovered at the beginning of the 20th century (Wood, 1902). This is a major defect for spectroscopy but its selective nature constitutes the basis of Plasmonics. It has been shown theoretically and verified experimentally that the absorption can reach the totality of the incident power (Hutley and Maystre, 1976). Surprisingly, this total absorption is generally obtained with very shallow gratings.
Plasmonics has a crucial importance in nanophotonics, for example for its potential applications to optical communications. Indeed, the surface plasmon is nearly confined at the vicinity of the air-metal interface, thus it is possible to generate several uncoupled surface plasmons on closely separated regions of this air-metal interface. The propagation of surface plasmons in the visible, near infrared and near ultraviolet regions cannot be envisaged over long distances, due to the high losses, but distances of several dozen microns can be reached. This is sufficient for envisaging the use of surface plasmon propagation in photonic circuits (Ozbay, 2006).
Another kind of absorption can be achieved from very deep metallic gratings. Since the modulated region of the grating behaves like an impedance adaptor between air and metal regions, the grating acts like a non-selective light absorber. It can be used for example for increasing the efficiency of photovoltaic cells (Teperik et al., 2008). Finally, let us notice that excitation of surface plasmons by gratings is one of the keys that explain the phenomenon of extraordinary transmission by hole arrays like metallic inductive grids (Ebbesen et al., 1998).
It is not possible to give an exhaustive description of all the applications of diffraction gratings. Some of the most important have to be cited here, like their use as instruments for distance and shape measurements (Hutley, 1980), beam samplers for high power lasers, photolithography masks(see next section), light filters for recording color images as surface-relief structures (Knop, 1978), superlenses in nanophotonics, light pulse compression (Treacy, 1969). The diffraction grating is also the basic instrument of Diffractive Optics (Turunen and Wyrovsky, 1998).
A rigorous grating theory must be deduced rigorously from the basic laws of electromagnetics in the form of a mathematical problem which can be solved on a computer. One may ask if a rigorous theory is necessary to investigate the properties of gratings? Many approximate theories are available. The most famous one uses the Kirchhoff approximation, which considers that any point of the grating profile behaves like an infinite plane interface tangent to the profile at that point (Beckmann and Spizzichino, 1987). These approximate theories are very simple to handle and to implement on computers. They may provide accurate results, under certain drastic conditions. Among these, the wavelength of light must be much smaller than the grating period, and the grating profile must be shallow. However, we have noticed in the preceding section that for spectroscopic purpose, gratings are generally used in conditions where the number of y-propagative orders is very small (often equal to 2). From the grating formula, this property entails that the wavelength of light has the same order of magnitude as the grating period and in that case such approximations fail.
First, we describe briefly one of the most popular rigorous theories of gratings: the Rigorous Coupled-Wave Analysis (RCWA), which is also the most simple one (Moharam and Gaylord, 1986). We have expressed the field above and below the intermediate region of the grating (<math></math> 0<y<y_m></y_m>) y="y_M" z)="\nabla" =" *" H_z="0. &lt;/math&gt; Finally," \sum_{n="&#45;\infty&#125;^&#123;+\infty&#125;&#123;" \sum_{m="1&#125;^&#123;4N+2&#125;&#123;" d\mathbf{F}/dy="\mathbf&#123;A&#125;(y)\mathbf&#123;F&#125;&lt;/math&gt;" \label{eq22} d\mathbf{F}/dy="\mathbf&#123;A&#125;(y)\mathbf&#123;F&#125;,\qquad" \mathbf{F}="\begin&#123;bmatrix&#125;" <math>\mathbf{S}_q="\mathbf&#123;V&#125;_q" \label{eq23} \mathbf{F}="1&#125;^&#123;4N+2&#125;" <math>y="y_M&lt;/math&gt;"
="" electromagnetics.
="" used.
="" 1949). Holographic grating: a photoresist is exposed to a fringe system
="" grooves).
="" micro-lithography.
="" 396.
=""></y<y_m></y_m>)>3-540-10193-4.