VYSOKÉ UENÍ TECHNICKÉ V BRN BRNO UNIVERSITY OF TECHNOLOGY
FAKULTA STROJNÍHO INENÝRSTVÍ ÚSTAV FYZIKÁLNÍHO INENÝRSTVÍ
FACULTY OF MECHANICAL ENGINEERING INSTITUTE OF PHYSICAL ENGINEERING
OPTIMALIZACE PROCEDURY M ENÍ OPTICKÉHO ZÁENÍ ROZPTÝLENÉHO PEVNÝMI T LESY PROVÁD NÉHO SCATTERMETEREM OPTIMIZING MEASUREMENT PROCEDURE OF THE OPTICAL RADIATION SCATTERED BY SOLID SURFACE PERFORMED BY THE SCATTERMETER
DIPLOMOVÁ PRÁCE MASTER'S THESIS
AUTOR PRÁCE Bc. JAKUB KLUS AUTHOR
VEDOUCÍ PRÁCE prof. RNDr. MILOSLAV OHLÍDAL, CSc. SUPERVISOR
BRNO 2013
Abstrakt Tato práce se zabývá optimalizací procedury měření optického záření rozptýleného pevnými tělesy prováděného pomocí scatterometru. V práci je popsáno několik metod použitelných k nalezení souboru pozic nepřekrývajících se kruhových detektorů na povrchu polokoule (i koule). Výsledky práce jsou aplikovány na samotný měřicí přístroj a je ukázána shoda experimentálních výsledků získaných před i po optimalizaci. Výhodu nového uspořádání lze nalézt především v převedení dřívějšího způsobu měření na soubor nezávislých měření, což má význam pro matematické zpracování výsledků.
Summary This work concerns with optimizing the measurement procedure of the optical radi- ation scattered from a solid surface performed by the scattermeter. The major part is devoted to the description of methods of finding set of positions of non-overlapping circular detectors on a surface of a hemisphere (or a sphere). Results of the optimiza- tion are applied to the measurement device. The agreement of the previous and new measurement procedures is shown in the end of this work. The crucial advantage of the new procedure is that the new set of point allows us to treat the measured values independently. This is useful for the mathematical processing of the results.
Klíčová slova Rozptyl elektromagnetického záření, scattermeter, drsné povrchy, solární články
Keywords Scattering of the electromagnetic radiation, scattermeter, rought surfaces, solar cells
KLUS, J.Optimalizace procedury měření optického záření rozptýleného pevnými tělesy prováděného scattermeterem. Brno: Vysoké učení technické v Brně, Fakulta strojního inženýrství, 2013. 60 s. Vedoucí prof. RNDr. Miloslav Ohlídal, CSc.
Hereby I declare that I have created this work autonomously under a scientific supervi- sion of Prof. Miloslav Ohlídal. All sources, references and literature used or excerpted during the elaboration of this work are properly cited and listed in the complete refer- ence.
Bc. Jakub Klus
I would like to thank to my supervisor Prof. Miloslav Ohlídal for his suggestions and advices. Furthermore I would like to thank to Ing. Pavel Nádaský for fruitful discus- sions and help in the laboratory. I also thank my parents for the financial support dur- ing studies.
Bc. Jakub Klus
CONTENTS Contents
1 Introduction3
2 Scattermeter II5 2.1 Scanning procedure...... 5 2.2 Angular distribution of the scattered light intensity...... 6 2.3 New type of experiment...... 6
3 Mathematical description of the problem9 3.1 Tammes‘s problem...... 9 3.2 Thomson’s problem...... 9 3.3 Riesz’s s-energy...... 10
4 Optimizing positions 11 4.1 Implementation of the steepest descent method...... 11 4.2 Parallelization of the procedure...... 13 4.3 Other approaches...... 16
5 Initial point configurations 17 5.1 Spiral points...... 17 5.2 Polyhedral breakdown systems...... 18 5.3 Parallel rings...... 23
6 Results 27 6.1 Polyhedral breakdown systems...... 27 6.2 Spiral points...... 37 6.3 Parallel rings...... 38 6.4 Repulsion forces...... 39
7 Applications 43 7.1 Visualizing spherical points...... 43 7.2 Generating the measurement sequence...... 45
8 Experiment 51 8.1 Time consumption...... 51 8.2 Angular distribution...... 52 8.3 Comparison with integrating sphere...... 53
9 Conclusions 55
10 Bibliography 57
11 List of symbols and abbreviations 59
1
1. INTRODUCTION 1. Introduction Since 2006 a device measuring an angular distribution of the light intensity scattered from a surface of solids had been developed in The Laboratory of Coherent Optics, IPE FME BUT. The basic purpose of the device was to verify simulations of the light scattering from studied surfaces, which were done within the framework of geomet- rical optics. This led to the quite extraordinary name of the device; Scattermeter. An original design from Gründling [1], nowadays called Scattermeter I (SM I), was com- plicated. Because of too many degrees of freedom the device was unable to repeat the same measurement conditions automatically. The construction itself had some disad- vantages, namely a large shadowed area in the detector space and a high sensitivity to the ambient light. An improved design came from the cooperation of Gründling and Brilla in 2009. New construction of the device is described in Brilla’s master thesis [2], the device is fa- miliarly called Scattermeter II (SM II). The ambient light was eliminated by housing the device into a lightproof box. Degrees of freedom were reduced to only two computer controlled coordinates. The light source was dismounted from the device and placed on an optical bench, to achieve exchangeability of various sources, which brought more measurement setups. The nearly completed device was ignored for more than a year, during this time all measurements were done on SM I. In 2011 a new master thesis was assigned to Nádaský focusing on measurements of solar cells by means of SM II. Since Brilla left SM II unfinished, Nádaský has struggled with getting experimental results comparable with SM I. The whole process of finalizing SM II is described in his master thesis [3]. The main results are changes in the construction of the detector holder including the detector arm and the detailed description of setting up the experiment. Nevertheless the aim of [3] was not only to get the device running but also to classify certain solar cells provided by the company Solartec s.r.o. With samples provided, SM II was improved to determine the angular distribution of the light intensity scattered from the surface of solids as precisely as we believe is possible. The real capabilities are to be summarized in the next chapter of this work. During the work of Nádaský on his thesis, there was an unanswered question: “Can we somehow compare our experimental results with other types of measurements?” Fortunately there are some measurements done directly on solar cells in Solartec s.r.o. using an integrating sphere enabling us to determine the energy scattered from the solar cell surface. When Nádaský finally received all the data he needed to compare performance of SM II with the integrating sphere, the results were completely differ- ent. After discussions in our laboratory we came to the conclusion that this difference can be caused by two main factors: a non-zero background signal on the detector and measuring in overlapping detector areas. The aim of this work is to find a setup of non-overlapping detector area positions, so that the results acquired using SM II are directly comparable with the data mea- sured by means of the integrating sphere. This should allow us to quantify the energy scattered in different directions. Another improvement is a reduction of the detector measurement positions which should result in a shorter acquisition time. The next chapter is devoted to the current state of SM II, its capabilities, measure- ment process and problems which should be solved further in this thesis. Chapter 3 in- troduces the mathematical background of this thesis, basic definitions and approaches.
3 Chapter 4 describes a simple numerical method of optimizing detector positions. The following chapter shows three methods of generating non-overlapping detector setups, which are used in this work as initial setups. Further then in the sixth chapter the final tested configurations of detector positions are studied and their features are character- ized. Chapter 8 is devoted to the application of the results on the measurement process of SM II. In the last two chapters we compare old and new measurement setups and then summarize the results of this work.
4 2. SCATTERMETER II 2. Scattermeter II Scattermeter II (SM II) is a second generation of a unique device developed in The Lab- oratory of coherent optics, Institute of Physical Engineering, Faculty of Mechanical En- gineering, Brno University of Technology. The main goal of it is to measure the angular distribution of the electromagnetic radiation scattered from the rough surfaces. Fig- ure 2.1 illustrates the basic parts of SM II. The centre of the base plate determines the z-rotation axis (red dashed line). By rotating around this axis one can setup different in- cident angles of the source beam. The main motor unit, which moves the azimuth ring, is hidden inside the main stand, which is attached to the base plate with four screws. The centre of the azimuth ring defines the azimuth rotation axis (blue dashed line). A sample is mounted in the holder, which can be extruded and tilted to align the normal of the sample with the azimuth rotation axis. The declination pivot consists of a holder and a declination motor which is placed on the azimuth ring. To ease the main motor unit load there is a counterweight on the opposite side of the ring. The declination motor shaft defines the declination axis (green dashed line).
Figure 2.1: Schematic illustration of SM II without the source of radiation. There are three axes of rotation in the figure: z-rotation axis (red dashed line), azimuth rotation axis (blue dashed line) and declination axis (green dashed line). Note: The centre of the sample is aligned with the rotation centre before each measurement.
2.1. Scanning procedure
SM II is capable of measuring the scattered light intensity in discrete positions on a surface of a hemisphere. The positions are achieved by rotation around the azimuth and the declination axis with a constant radius. The radius of the hemisphere is de- termined by the detector arm length. Angular coordinates are not defined in the same way as it is usual for spherical coordinates in mathematics. Figure 2.2 shows that we
5 2.2. ANGULAR DISTRIBUTION OF THE SCATTERED LIGHT INTENSITY define the azimuth angle φ from 0◦ to 180◦ and the declination angle θ from 90◦ to 90◦, − the declination 0◦ is on the pole of the hemisphere. Scanning starts by setting up the azimuth and the declination angle to the initial values. Then the declination angle in- creases by steps1 towards the end position. In the last position of the declination angle the azimuth angle increases once and then the declination decreases towards the initial position. This repeats until the azimuth angle reaches its last position.
0◦ θ
detector incident beam
sample
0◦ φ
Figure 2.2: Modified hemispherical coordinates.
2.2. Angular distribution of the scattered light intensity
So far the best investigated measurement according to [3] is the angular distribution of the scattered light intensity. Using a red (λ = 635 nm) laser source, SM II is capable of measuring the intensity of light reflected from a rough sample. From these results one can determine the main scattering directions for a certain sample. Rotation along the z-axis changes the incident angle of the laser beam. SM II is also equipped with an infrared (λ = 1550 nm) laser source. Current re- search is done on silicon solar cells, which are transparent for this laser. This allows us to also measure the angular distribution of the radiation transmitted through a rough sample. Unfortunately the construction of SM II does not allow us to change the angle of incidence for the transmission regime. Nevertheless with both the total transmission and the total reflection we can try to determine the energy trapped inside the solar cell as the difference between them. This new type of the measurement requires some data processing which is described in the next section.
2.3. New type of experiment
For angular measurements the actual values acquired from the device are not impor- tant. Because it was validated in [3] that the detector response is linear, we can nor-
1This step also determines the resolution of the measurement, the best resolution is limited by the controlling software and the smallest step is 1◦.
6 2. SCATTERMETER II malize the distribution without influencing the results. Since we want to measure the total energy reflected and transmitted we need to convert the signal from photodiode to power and then integrate it over the whole hemisphere. As the first approximation we have summed the calibrated directional values. This led to the total power reflected from the sample greater than the power of the incident beam. We have found out that the problem was in overlapping of the detector positions. According to [2] the hemi- sphere radius is chosen in the way that detector’s positions on the equator are touching each other. With a decreasing declination angle the detector positions become more and more overlapping. This can be compensated by multiplying the value measured with the detector by the ratio of the detector’s active area and the area determined by the neighbouring detectors positions [3]. In this work we deal with a completely differ- ent approach, we will find a non-overlapping detector positions configuration with the maximal surface coverage. The next chapter deals with the mathematical definition of such a problem.
7
3. MATHEMATICAL DESCRIPTION OF THE PROBLEM 3. Mathematical description of the problem Although a grid of equidistant points in angular coordinates of the hemispherical sys- tem seems to be useful for the visualization, it is not applicable for integrating. Instead of a scanning detector let us consider an array of detectors1. Setting up this array in the grid is not possible. So even if it is not the goal of this work, its results are also valid for constructing the Scattermeter as a static array of equivalent detectors. To achieve better mathematical description we can consider not the hemispherical distribution, but the spherical one. Then we are searching for distribution of non- overlapping spherical caps with equal viewing angle on a surface of hemisphere. These problems are also referred as spherical packing problems.
3.1. Tammes‘s problem
The spherical packing problem is also referred as Tammes’s problem [4]. Tammes was a Dutch botanist, who studied exit places on the surface of spherical pollen grains. The exit places on the surface of the pollen grain can be assumed to be circular caps of the same area. Tammes found out that the number of exits is related to the size of pollen grain and an area needed for one exit place [5]. The problem can be formulated as follows: "How many spherical caps of the same radius can be placed on a unit sphere without overlapping?" Another formulation is to find an arrangement of N points on a unit sphere that maximizes the minimum distance between any two of them. This is equivalent to the original Tammes’s problem and also to our problem by thinking of the N points as the centres of spherical caps of an equal maximum radius.
3.2. Thomson’s problem
A different approach for the similar problem was done by Sir Joseph John Thomson. In 1904 he proposed the model of the atom structure as follows: The atom is composed of electrons surrounded by a soup of positive charge to balance electron’s negative charge. Thomson’s problem concerns with position of electrons in a spherical atom shell and can be formulated in this manner: Find a minimum energy configuration of N electrons that repel each other with a force given by Coulomb’s law on the surface of a sphere. The more mathematical way to express this problem is: Let r1, r2,..., rN be a col- lection of N distinct points on the unit sphere centred at the origin of the coordinate system. The energy of this configuration of points is defined to be: X 1 . (3.1) ri rj i 1 This layout was rejected because of the total cost and extremely complicated calibration of the de- tector array. 9 3.3. RIESZ’S S-ENERGY 3.3. Riesz’s s-energy The bridge between Tammes’s problem and Thomson’s problem is Riesz’s s-energy. The Riesz’s energy is called after M. Riesz, who derived a number of important proper- ties of Riesz potentials [6]. For a collection of N 2 distinct points XN = x1,..., xN ≥ { } and s > 0, the discrete Riesz’s s-energy of XN is defined by: X 1 Es (XN ) = s . (3.2) xj xk j6=k | − | Searching for the minimal Es, when s = 1 and N is fixed we are solving Thomson’s problem. On the other hand when s and N is fixed, we get the best packing prob- → ∞ lem, which is the problem of finding N-point configuration with the largest separation radius [7]. From the solution of best packing problems in the plane, one can expect that the best solution for the spherical surface will consist of hexagons [8]. But according to Euler’s characteristic2: F E + V = 2, where F is the number of faces, E is the number of − edges and V is the number of vertices, sphere cannot be tiled only using hexagons. Let us suppose that sphere will be divided using only hexagons and pentagons. Assuming that exactly 3 edges originate from each vertex it follows from Euler’s characteristic that there has to be exactly 12 pentagons in order to fit on the sphere [8]. A common way to solve the best packing problem is based on numerical optimiza- tion using a computer, but with growing number of points the problem starts to be very challenging. Largest list of spherical codes is provided on http://neilsloane.com/. The author of the web is also co-author of [9] where the packing problem is described by means of mathematics. Computations uncovered that there are many local min- ima for the s-energy problem, which are not global minima, but their energies are very close to the global minimum [8]. It is estimated that the number of distinct local min- ima grows exponentially with N. This is one of the reasons why the packing problem can be used for testing the global optimizing methods [8]. 2Euler’s characteristic is an topological invariant which is for the case of a convex polyhedron (sphere) equal to 2 [8]. 10 4. OPTIMIZING POSITIONS 4. Optimizing positions In this chapter we will introduce the numerical optimizing method for finding a mini- mum of the s-energy. There are many optimization methods, but as the first approach most commonly used is the steepest descent method [4, 10] as the basic and simple method [11]. The steepest descent (gradient) method is a first-order optimization algorithm. Let us have a function of n variables f (x1, . . . , xn) which is differentiable on R (for con- venience we will write x = (x1, . . . , xn) and f(x) in vector notation). The idea of this method is to find a minimum of f(x) by repeatedly computing minima of a function g(t) of a single variable t [12]. Suppose that f(x) has a minimum at X0 and we are starting at a point x1. Then in the first step we search for a minimum of f(x) along the straight line in the direction of x f(x), which is closest to X0. We change the x −∇| 1 starting from the point x1 by the step t1 x f(x): − · ∇| 1 z (t) = x t1 x f(x). (4.1) − ∇| 1 To find t1 we minimize the function: g(t) = f (z (t)) . (4.2) At the minimum of g(t), z(t) is closest to X0, so we take z(t1) in the next step as a start point x2. The optimization process is stopped when changes of subsequent ti are small enough. 4.1. Implementation of the steepest descent method In our case the steepest descent method comes from the physical background of the problem. We assume a set of N points and their position vectors r1,..., rN. Let us define scalar potential of the i-th point ϕi(r) as follows: 1 ϕi (r) = Q s , (4.3) · r ri | − | where r is a position in R3, s defines the potential (see Chapter 3) and Q is an arbitrary constant. The direction of a steepest descent is: Φ(r, r1,..., rN) = F(r, r1,..., rN), (4.4) − ∇ where F(r, r1,..., rN) is the net force and Φ(r, r1,..., rN) is the total potential in the position r given by the equation: N N X X 1 Φ(r, r1,..., rN) = ϕi (r) = Q s (4.5) · r ri i=1 i=1 | − | The value of Φ(r, r1,..., rN) can be determined numerically by evaluating all the −∇ partial derivatives of Φ(r, r1,..., rN) for all components using for example the finite difference approximation: ∂Φ Φ(x, x1, . . . , xi + h, . . . , xN ) Φ(x, x1, . . . , xi, . . . , xN ) = − , (4.6) ∂xi h 11 4.1. IMPLEMENTATION OF THE STEEPEST DESCENT METHOD ∂Φ Φ(y, y1, . . . , yi + h, . . . , yN ) Φ(y, y1, . . . , xi, . . . , yN ) = − , (4.7) ∂yi h ∂Φ Φ(z, z1, . . . , zi + h, . . . , zN ) Φ(z, z1, . . . , zi, . . . , zN ) = − , (4.8) ∂zi h where h 0 is the step of differentiation. This process involves 4N calculations of → Φ(r, r1,..., rN) and every calculation of Φ(r, r1,..., rN) requires N calculations of ϕi(r). This can lead to a very intensive usage of computing power. However, we can avoid calculating Φ(r, r1,..., rN) numerically by just expressing the explicit equation for −∇ net force F(r, r1,..., rN). From the equation (4.3) we can derive the formula for the force caused by the i-th point Fi(r) analytically: r ri Fi(r) = Q1 − s+2 , (4.9) · r ri | − | where Q1 is arbitrary constant. Then the net force F(r, r1,..., rN) of the set r1,..., rN in the position r can be estimated as a vector sum of all N forces caused by points of our set: N N X X r ri F(r, r1,..., rN) = Fi(r) = Q1 − s+2 = Φ(r, r1,..., rN). (4.10) · r ri −∇ i=1 i=1 | − | According to equation (3.2) we define the energy of our system of N particles XN = r1,..., rN as follows: { } N N X 1 X X Es (XN)) = s = ϕk (rj) . (4.11) rj rk j6=k | − | j>0 k>0,k6=j The direction of the steepest descent for i-th point is then equal to the net force Fi,net acting on the i-th point: N X rj ri Fi,net = Q1 − s+2 . (4.12) rj ri j=1,j6=i | − | 0 0 0 After determining the forces Fi,net we search for a new system of particles X = r ,..., r N { 1 N} which fulfil: 0 r = ri + t Fi,net, (4.13) i · 0 where t is changed until Es(XN ) > Es(XN ). The algorithm works with the Cartesian coordinate system in Euclidean space R3. It can be described in eight steps: 1. Let us have an initial set of N points XN = r1,..., rN on a unit sphere and an { } initial value of t. Calculate the s-energy Es(XN ). 2. Calculate the net force Fi,net acting on each point. 3. Subtract the normal component from the net force. 4. Calculate set of points X0 = r0 ,..., r0 , according to (4.13). N { 1 N} 0 0 5. Normalize vectors r1,..., rN. 12 4. OPTIMIZING POSITIONS 0 6. Calculate s energy Es(XN ). 0 7. If Es(XN ) > Es(XN ) decrease the value of t and continue with step 4. 0 0 8. If Es(XN ) < Es(XN ) use the XN as a new value of XN and continue with step 2. Steps 3 and 5 are important for keeping the points on the unit sphere. Normalization in step 5 is straightforward. All vectors with the same length are on the same spherical surface. As it is depictured in Figure 4.1, subtracting the normal component Fi,n from the net force Fi,net gives us a direction in the tangential plane Fi,t which has better confinement to the unit sphere surface. Fi,net ri Fi,t Fi,n Figure 4.1: Illustration of the subtraction of the normal component Fi,n from the net force Fi,net. The resulting vector Fi,t is lying in the blue tangential plane of the sphere. 4.2. Parallelization of the procedure The optimization procedure described above has quadratic complexity (n2), the num- O ber of calculations can be reduced to half. In the case of calculating the energy Es(XN ), we can use the fact that the potential is symmetric and is a function only of the mutual distance: 1 1 ϕi (rj) = C s = C s = ϕj (ri) . (4.14) · rj ri · ri rj | − | | − | Then the s-energy of the system Es(XN ) can be evaluated by only half the additions: N N X X Es (XN ) = 2 ϕk (rj) . (4.15) · j>0 k>j To achieve the same performance for calculating the net force we can use Newton’s third law: "The forces of two bodies on each other are always equal and are directed in 13 4.2. PARALLELIZATION OF THE PROCEDURE opposite orientation." There is no simple formula to prove the impact on the procedure, but when we calculate the force that the j-th point is acting on the i-th point we can store it also as an opposite force of the i-th point acting on the j-th and do not evaluate it again: N N X ri rj X ri rj Fj(ri) = C1 − s+2 = Fi(rj) = C1 − s+2 . (4.16) ri rj − − ri rj j=1,j6=i | − | j=1,j6=i | − | Although we have halved the total number of calculations we still have to deal with an enormous number of operations. This is the reason why we have implemented the algorithm in programming language C for parallel computers1. We will not discuss the implementation here, but we will illustrate how the work is distributed between workers to shorten the time needed for the one iteration of the optimization. Because the load distribution process is the same for both the energy and the net force computation we will show it only for the energy. Assume that we have w workers and a set of N points (we treat the point as a task which has to be solved by a worker). To calculate the i-th point contribution we have to evaluate (4.3) (N i)times, the total N(N−1) − number of calculations is then n = 2 , please note that last point has no other to interact with, so we actually have to distribute only N-1 tasks. The first approach is to assign the same number of tasks wi to each worker. In principle it is not exactly possible, so we have assigned N−1 tasks to the first w 1 workers and the rest to the last worker. w − Let us demonstrate the distribution on the example of 30 points and 6 workers. We have assigned four tasks to each of first five workers and then the remaining ten tasks to the last one. As we can see in Figure 4.2, the first worker is doing more than twice