Open access

Features of Electromagnetic Waves Scattering by Surface Fractal Structures

Written By

O. Yu. Semchuk and M. Willander

Submitted: 07 October 2010 Published: 05 July 2011

DOI: 10.5772/16281

From the Edited Volume

Behaviour of Electromagnetic Waves in Different Media and Structures

Edited by Ali Akdagli

Chapter metrics overview

2,947 Chapter Downloads

View Full Metrics

1. Introduction

Accurate measurement of surface roughness of machined workpieces is of fundamental importance particularly in the precision engineering and manufacturing industry. This is mainly due to the more stringent demand on material quality as well as the miniaturization of product components in these industries [1-3]. For instance, in the disk drive industry, to maintain the quality of the electrical components mounted on an optical disk, the surface roughness of the disk must be accurately measured and controlled. Hence, the surface finish, normally expressed in terms of surface roughness, is a critical parameter used for the acceptance or rejection of a product.

Surface roughness is usually determined by a mechanical stylus profilometer. However, the stylus technique has certain limitations: the mechanical contact between the stylus and the object can cause deformations or damage on the specimen surface and it is a pointwise measurement method and is time consuming. Hence a noncontact and more speed optical method would be attractive. Different optical noncontact methods for measuring surface roughness have been developed mainly based on reflected light detection, focus error detection, laser scattering, speckle and the interference method [4-10]. Some of these have good resolutions and are being applied in some sectors where mechanical measuring methods previously enjoyed clear predominance. Among these methods, the light scattering method [11] which is a noncontact area-averaging technique, is potentially more speedy for surface inspection than other profiling techniques particularly the traditional stylus technique. Other commercially available products such as the scanning tunneling microscope (STM), the atomic force microscope (AFM) and subwavelength photoresist gratings [12-15], which are pointwise techniques, are used mainly for optically smooth surfaces with roughnesses in the nanometer range.

In this chapter in the frame of the Kirchhoff method (scalar model) the average coefficient of light scattering by surface fractal structures was calculated. A normalized band-limited Weierstrass function is presented for modeling 2D fractal rough surfaces. On the basis of numerical calculation of average scattering coefficient the scattering indicatrises diagrams for various surfaces and falling angles were calculated. The analysis of the diagrams results in the following conclusions: the scattering is symmetrically concerning a plane of fall; with increase the degree of a surface calibration the picture becomes complicated; the greatest intensity of a scattering wave is observed in a mirror direction; there are other direction in which the bursts of intensity are observed.

Advertisement

2. Fractal model for two-dimensional rough surfaces

At theoretical research of processes of electromagnetic waves scattering selfsimilar heterogeneous objects (by rough surfaces) is a necessity to use the mathematical models of dispersive objects. As a basic dispersive object we will choose a rough surface. As is generally known, she is described by the function z ( x , y ) of rejections z of points of M of surface from a supporting plane (x,y) (fig.1) and requires the direct task of relief to the surface.

Figure 1.

Schematic image of rough surface

There are different modifications of Weierstrass–Mandelbrot function in the modern models of rough surface are used. For a design a rough surface we is used the Weierstrass limited to the stripe function [3,4]

z ( x , y ) = c w n = 0 N 1 m = 1 M q ( D 3 ) n sin { K q n [ x cos 2 π m M + y sin 2 π m M ] + φ n m } , E1

where c w is a constant which ensures that W(x, y) has a unit perturbation amplitude; q(q> 1) is the fundamental spatial frequency; D (2 < D< 3) is the fractal dimension; K is the fundamental wave number; N and M are number of tones, and φ n m is a phase term that has a uniform distribution over the interval [ π , π ] .

The above function is a combination of both deterministic periodic and random structures. This function is anisotropic in the two directions if M and N are not too large. It has a large derivative and is self similar. It is a multi-scale surface that has same roughness down to some fine scales. Since natural surfaces are generally neither purely random nor purely periodic and often anisotropic, the above proposed function is a good candidate for modeling natural surfaces.

The phases φ n m can be chosen determinedly or casually, receiving accordingly determine or stochastic function z ( x , y ) . We further shall consider φ n m as casual values, which in regular distributed on a piece [ π ; π ] . With each particular choice of numerical meanings all N × M phases φ n m (for example, with the help of the generator of random numbers) we receive particular (with the beforehand chosen meanings of parameters c w , q , K , D , N , M ) realization of function z ( x , y ) . The every possible realizations of function z ( x , y ) form ensemble of surfaces.

A deviation of points of a rough surface from a basic plane proportional c w , therefore this parameter is connected to height of inequalities of a structure of a surface. Further it is found to set a rough surface, specifying root-mean-square height of its structure σ , which is determined by such grade:

σ ( h 2 ) , E2

where h = z ( x , y ) , ( ... ) = n = 0 N 1 m = 1 M π π d φ n m 2 π ( ... ) - averaging on ensemble of surfaces.

The connection between c w and σ can be established, directly calculating integrals:

σ = [ n = 0 N 1 m = 1 M π π d φ n m 2 π z 2 ( x , y ) ] 1 2 = c w [ M ( 1 q 2 N ( D 3 ) ) 2 ( 1 q 2 ( D 3 ) ) ] 1 2 . E3

So, the rough surface in our model is described by function from six parameters: c w (or ), q , K , D , N , M . The influence of different parameters on a kind of a surface can be investigated analytically, and also studying structures of surfaces constructed by results of numerical accounts of Weierstrass function. Analysis of the surface profiles built by us on results of numeral calculations (fig. 2) due to the next conclusions:

  • the wave number K sets length of a wave of the basic harmonic of a surface;

  • the numbers N , M , D and q determine a degree of a surface calibration at the expense of imposing on the basic wave from additional harmonics, and N and M determine the number of harmonics, which are imposed;

  • D
  • q

Let's notice that with increase N , M , D and q the spatial uniformity of a surface on a large scale is increased also.

Figure 2.

Examples of rough surface by the Weierstrass function K = 2 π ; N = M = 5 ; σ = 1 . D = 2 , 1 ; D = 2 , 5 ; D = 2 , 9 (from above to the bottom) q = 1 , 1 ; q = 3 ; q = 7 (from left to right)

By means of the original program worked out by us in the environment of Mathematika 5.1 there was the created base of these various types of fractal dispersive surfaces on the basis of Weierstrass function.

Influence each of parameters q , K , D , N , M on character of profile of surface it appears difficult enough and determined by values all other parameters. So, for example, at a value D = 2 , 1 , what near to minimum ( D = 2 ), the increase of size q does not almost change the type of surface (see the first column on fig.2). With the increase of size D the profile of surface becomes more sensible to the value q (see the second and third columns on fig.2). Will notice that with an increase N , M , D and q increases and spatial homogeneity of surface on grand the scale: large-scale "hills" disappear, and finely scale heterogeneities remind a more mesentery on a flat surface.

Advertisement

3. Electromagnetic wave scattering on surface fractal structures

At falling of electromagnetic wave there is her dispersion on the area of rough surface - the removed wave scattering not only in direction of floppy, and, in general speaking, in different directions. Intensity of the radiation dissipated in that or other direction is determined by both parameters actually surfaces (by a reflectivity, in high, by a form and character of location of inequalities) and parameters of falling wave (frequency, polarization) and parameters of geometry of experiment (corner of falling). The task of this subdivision is establishing a connection between intensity of the light dissipated by a fractal surface in that or other direction, and parameters of surface.

Figure 3.

The scheme of experiment on light scattering by fractal surface: S is a scattering surface; D-detector, θ 1 is a falling angle; θ 2 is a polar angle; θ 3 is an azimuthally angle

The initial light wave falls on a rough surface S under a angle θ 1 and scattering in all directions. The scattering wave is observed by means of the detector D in a direction which is characterized by a polar angle θ 2 and an azimuthally angle θ 3 . The measured size is intensity of light I s scattered at a direction ( θ 2 , θ 3 ) . Our purpose is construction scattering indicatrise of an electromagnetic wave by a fractal surface (1).

As I s = E s E s * (where E s is an electric field of the scattering wave in complex representation) that the problem of a finding I s is reduced to a finding of the scattered field E s .

The scattered field we shall find behind Kirchhoff method [16], and considering complexity of a problem, we shall take advantage of more simple scalar variant of the theory according to which the electromagnetic field is described by scalar size. Thus we lose an opportunity to analyze polarizing effects

The base formula of a Kirchhoff method allows to find the scattered field under such conditions:

  • the falling wave is monochromatic and plane;

  • a scattered surface rough inside of some rectangular (-X <x0<X, -Y <y0<Y) and corpulent outside of its borders;

  • the size of a rough site much greater for length of a falling wave;

  • all points of a surface have the ended gradient;

  • the reflection coefficient identical to all points of a surface;

  • the scattering field is observed in a wave zone, i.e. is far enough from a scattering surface.

Under these conditions the scattered field is given by

E s ( r ) = i k r F ( θ 1 , θ 2 , θ 3 ) exp ( i k r ) 2 π r S 0 exp [ i k φ ( x 0 , y 0 ) ] d x 0 d y 0 + E e ( r ) E4

Where k is the wave number of falling wave; F ( θ 1 , θ 2 , θ 3 ) = R 2 C ( A 2 + B 2 + C 2 ) is a angle factor; R - scattering coefficient; φ ( x 0 , y 0 ) = A x 0 + B y 0 + C h ( x 0 , y 0 ) is the phase function; h ( x 0 , y 0 ) = z ( x 0 , y 0 ) ; A = sin θ 1 sin θ 2 cos θ 3 ; B = sin θ 2 sin θ 3 ; C = cos θ 1 cos θ 2 ; E e ( r ) = R C exp ( i k r ) 4 π r ( A I 1 + B I 2 ) ,

I 1 = Y Y [ e i k φ ( X , y 0 ) e i k φ ( X , y 0 ) ] d y 0 , I 2 = X X [ e i k φ ( x 0 , Y ) e i k φ ( x 0 , Y ) ] d x 0 . E5

After calculation of integrals (4) and (5) by means of the formula

e i z sin ϕ = l = I l ( z ) e i l ϕ , E6

(where I l ( z ) is the Bessel function of the whole order), we receive

E s ( r ) = 2 i k F X Y exp ( i k r ) π r l { r s } { [ u v I l u v ( ξ u ) ] e i n m l n m φ n m } sin c ( k c X ) sin c ( k s Y ) + E e ( r ) E7

where

F = F ( θ 1 , θ 2 , θ 3 ) , l { r s } l 0 , 1 = l 0 , 2 = ... l ( N 1 ) , M = , u v u = 1 N 1 v = 0 M , n m n = 1 N 1 m = 0 M , ξ u k c w C q ( D 3 ) u , sin c , x sin x x k c k A + K n m q n l n m cos 2 π m M , k s k B + K n m q n l n m sin 2 π m M , E e ( r ) = i k X Y R C ( A 2 + B 2 ) e i k r π r sin c ( k A X ) sin c ( k B Y ) E8

Thus, expression (6) gives the decision of a problem about finding a field scattering by a fractal surface, within the limits of Kirchhoff method.

Now under the formula (4) it is possible to calculate intensity of scattered waves if to set parameters of a disseminating surface c w (or) σ D , q , K , N , M , X , Y , ϕ n m , parameter k (or λ = 2 π k ) a falling wave and parameters θ 1 , θ 2 , θ 3 of geometry of experiment. This intensity will be to characterize scattering on concrete realization of a surface z ( x , y ) (with a concrete set of casual phases ϕ n m ). For comparison of calculations with experimental data it is necessary to operate with average on ensemble of surfaces intensity ( I s ) = ( E s E s ) . Such intensity has appeared proportional intensity I 0 = ( 2 k X Y cos Θ 1 π r ) 2 of the wave reflected from the corresponding smooth basic surface, therefore for the theoretical analysis of results it is more convenient to use average scattering coefficient

( ρ s ) = ( I s ) I 0 . E9

After calculation ( I s ) and leaving from (6), we shall receive exact expression

( ρ s ) = [ F ( θ 1 , θ 2 , θ 3 ) cos Θ 1 ] 2 l { r s } { u v I l u v 2 ( ξ u ) sin c 2 ( k c X ) sin c 2 ( k s Y ) } + [ R ( A 2 + B 2 ) 2 C cos θ 1 ] 2 sin c 2 ( k A X ) sin c 2 ( k B Y ) E10

As expression (7) consist the infinite sum to use it for numerical calculations inconveniently. Essential simplification is reached in case of ξ n 1 . Using thus decomposition function in a line

I ν ( z ) = ( 3 2 ) ν k = 0 ( z 2 / 4 ) k k ! Γ ( ν + k + 1 ) E11

that rejecting members of orders, greater than ξ n 2 . We shall receive the approached expression for average scattering coefficient

( ρ s ) [ F ( θ 1 , θ 2 , θ 3 ) cos θ 1 ] 2 { [ 1 ( k σ C ) 2 ] sin c 2 ( k A X ) sin c 2 ( k B Y ) + 1 2 c f 2 n m q 2 ( D 3 ) n sin c 2 [ ( k A + K q n cos 2 π m M ) X ] sin c 2 [ ( k B + K q n sin 2 π m M ) Y ] } E12

where

c f k c w C = k σ C [ 2 M 1 q 2 ( D 3 ) 1 q 2 N ( D 3 ) ] 1 2 E13
Advertisement

4. Results of numerical calculations

On the basis of numerical calculations of average factor of dispersion under the formula (8) we had been constructed the average scattering coefficient ( ρ s ) from θ 2 and θ 3 . (scattering indicatrix diagrams) for different types of scattering surfaces. At the calculations we have supposed R = 1 , and consequently did not consider real dependence of reflection coefficient R from the length of a falling wave λ and a falling angle θ 1 . The received results are presented on fig. 4.

Figure 4.

Dependencies of the log ( ρ s ) from the angles θ2 and θ3 for the various type of fractal surfaces: a, a’, a’’ – the samples of rough surfaces, which the calculation of dispersion indexes was produced; from top to bottom the change of scattering index is rotined for three angles of incidence θ 1 = 30 , 40 , 60 0 (a-d, a’-d’, a’’-d’’) at N=5, M=10, D=2.9, q=1.1; n=2, M=3, D=2.5, q=3; N=5, M=10, D.2.5, q=3 accordingly

The analysis of schedules leads to such results:

  • Scattering is symmetric concerning of a falling plane;

  • The greatest intensity of the scattering wave is observed in a direction of mirror reflection;

  • There are other directions in which splashes in intensity are observed;

  • With increase in a calibration degree of surfaces (or with growth of its large-scale heterogeneity) the picture of scattering becomes complicated. Independence of the type of scattering surface there is dependence of the scattering coefficient from the incidence angle of light wave. As far as an increase of the incidence angle from 300 to 600 amounts of additional peaks diminishes. Is their most number observed at θ 1 = 30 0 . It is related to influence on the scattering process of the height of heterogeneity of the surface. At the increase of the angle of incidence of the falling light begins as though not to “notice” the height of non heterogeneity and deposit from them diminishes.

The noted features of dispersion are investigation of combination of chaoticness and self-similarity relief of scattering surface.

Advertisement

5. Conclusion

In this chapter in the frame of the Kirchhoff method the average coefficient of light scattering by surface fractal structures was calculated. A normalized band-limited Weierstrass function is presented for modeling 2D fractal rough surfaces. On the basis of numerical calculation of average scattering coefficient the scattering indicatrises diagrams for various surfaces and falling angles were calculated. The analysis of the diagrams results in the following conclusions: the scattering is symmetrically concerning a plane of fall; with increase the degree of a surface calibration the picture becomes complicated; the greatest intensity of a scattering wave is observed in a mirror direction; there are other direction in which the bursts of intensity are observed.

References

  1. 1. Bifano T. G. Fawcett H. E. Bierden P. A. 1997 Precision manufacture of optical disc master stampers. Precis. Eng. 20(1), 53 EOF 62 EOF .
  2. 2. Wilkinson P. et al. 1997 Surface finish parameters as diagnostics of tool wear in face milling, Wear 25(1), 47 EOF 54 EOF .
  3. 3. Sherrington I.’. Smith E. H. 1986 The significance of surface topography in engineering. Precis. Eng. 8(2). 79 87 .
  4. 4. Kaneami J. Hatazawa T. 1989 Measurement of surface profiles bv the focus method. Wear 134 221 229 .
  5. 5. Mitsui K. 1986 In-process sensors for surface roughness and their applications. Precis. Eng. 8(40), 212-220.
  6. 6. Baumgart J. W. Truckenbrodt H. 1998 Scatterometrv of honed surfaces. Opt. Eng. 37(5), 1435-1441.
  7. 7. Tay C. J. Toh L. S. Shang H. M. Zhang J. B. 1995 Whole-field determination of surface roughness bv speckle correlation. Appl. Opt. 34(13). 2324 2335 .
  8. 8. Peiponen K. E. Tsuboi T. 1990 Metal surface roughness and optical reflectance Opt. Laser Technol. 22(2). 127 130 .
  9. 9. Whitley J. Q. Kusy R. P. Mayhew M. J. Buckthat J. E. 1987 Surface roughness of stainless steel and electroformed nickle standards using HeNe laser. Opt. Laser Technol 19(4), 189-196.
  10. 10. Mitsui M. Sakai A. Kizuka O. 1988 Development of a high resolution sensor for surface roughness, Opt. Eng. 27(^6). 498 502 .
  11. 11. Vorburger T. V. Marx E. Lettieri T. R. 1993 Regimes of surface roughness measurable with light scattering. Appl Opt. 32(19!. 3401 3408 .
  12. 12. Raymond C. J. Murnane M. R. Naqvi H. S. S. Mcneil J. R. 1995 Metrology of subwavelength photoresist gratings usine optical scat-terometry. J. Vac. Sci. Technol B 13(4), 1484-1495.
  13. 13. Whitehouse D. J. 1991 Nanotechnologv instrumentation. Meas. Control 24(3), 37-46.
  14. 14. Madsen L. L J.Srgensen J. F. Carneiro K. Nielsen H. S. (1993-1994). Roughness of smooth surfaces: STM versus profilometers. Metro-logia 30 513-516.
  15. 15. Stedman M. 1992 The performance and limits of scanning probe microscopes. In Proc. Int. Congr. X-ray Optics and Microanalysis, 347 352 , Manchester. IOP Publishing Ltd.
  16. 16. Berry M. V. Levis Z. V. 1980 On the Weirstrass-Mandelbrot fractal function. Proc.Royal Soc. London A, 370 459

Written By

O. Yu. Semchuk and M. Willander

Submitted: 07 October 2010 Published: 05 July 2011