CN103426031B - A kind of optimization method of ellipsometer system parameter - Google Patents
A kind of optimization method of ellipsometer system parameter Download PDFInfo
- Publication number
- CN103426031B CN103426031B CN201310302341.XA CN201310302341A CN103426031B CN 103426031 B CN103426031 B CN 103426031B CN 201310302341 A CN201310302341 A CN 201310302341A CN 103426031 B CN103426031 B CN 103426031B
- Authority
- CN
- China
- Prior art keywords
- ellipsometer
- matrix
- system parameter
- parameter
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000011159 matrix material Substances 0.000 claims abstract description 114
- 238000005259 measurement Methods 0.000 claims description 28
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000009897 systematic effect Effects 0.000 abstract description 2
- 230000010287 polarization Effects 0.000 description 11
- 238000004088 simulation Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000572 ellipsometry Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013433 optimization analysis Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention discloses a kind of ellipsometer system parameter optimization method, comprise the system model first setting up ellipsometer to be optimized, and obtain the matrix W comprising ellipsometer system parameter; Then the conditional number of ellipsometer system matrix W under different ellipsometer system parameter is calculated; Again calculate the difference of the maxima and minima of matrix W conditional number under each ellipsometer system parameter, and by the arrangement of difference descending order; Last putting in order according to difference, optimize each ellipsometer system parameter, after repeating n suboptimization, choose the systematic parameter that arrangement is the most front, difference when trying to achieve minimal condition numerical value and the (n-1)th suboptimization after its n suboptimization between matrix W minimal condition numerical value, when difference size meets the demands, then think that the ellipsometer system parameter after the n-th suboptimization is the optimal system parameter of current ellipsometer.The inventive method process is clear and definite, simple, and the selected standard flexibility and changeability of optimized parameter, goes for the optimization of the ellipsometer of existing different types of structure.
Description
Technical Field
The invention belongs to the technical field of ellipsometers, and particularly relates to an ellipsometer system parameter optimization method.
Background
The ellipsometer is a measuring instrument most widely used in the field of ellipsometry, and is a general optical measuring instrument for acquiring information of a sample to be measured by using polarization characteristics of light.
In recent years, in order to adapt to different measurement conditions and the needs of user groups, ellipsometers of various configuration types have been developed, including a rotating polarizer type, a rotating analyzer type, a single rotation compensator type, a double rotation compensator type, and the like. When an ellipsometer is used to measure a sample to be measured, the measurement result often deviates from its true value to some extent. These deviations are caused by many reasons, including ellipsometer random noise, ellipsometer systematic errors, environmental random noise, and measurement artifacts. Even under the same deviation, the ellipsometer has different measurement results and sensitivity to the deviation according to different configured system parameters. Therefore, in order to reduce the influence of the deviation on the measurement result of the ellipsometer, a certain method must be used to perform optimization analysis on the system parameters of the ellipsometer, so as to obtain the optimal system parameters of the ellipsometer, thereby reducing the influence of the deviation on the measurement result of the ellipsometer as much as possible. During the ellipsometer design process, an optimal system parameter configuration for the current ellipsometer system configuration must be found.
Currently, for parameter optimization of an ellipsometer system, a main concern is to optimize the phase retardation of an ellipsometer compensator, for example, an optimal phase retardation of 127 ° is proposed to reduce the influence of the deviation on the ellipsometer measurement result. However, the ellipsometer has many parameters, such as the azimuth angles of the polarizer and the analyzer, the number of sampling points, the azimuth angle of the compensator, etc., which are also important influencing factors for the relationship between the deviation and the measurement result of the ellipsometer, and must be comprehensively considered for optimization.
Disclosure of Invention
The invention aims to provide an ellipsometer system parameter optimization method, which is simple and feasible and can be applied to ellipsometer measurement systems with various structural types.
In order to achieve the above object, the present invention provides an ellipsometer system parameter optimization method for determining an optimal system parameter of an ellipsometer, which is characterized in that the method includes:
(1) establishing a system model of the ellipsometer to be optimized, and arranging the system model of the ellipsometer into a form of a vector dot product and a vector dot product to obtain a system matrix W containing system parameters of the ellipsometer;
(2) calculating the condition number of the system matrix W under each ellipsometer system parameter;
(3) calculating the difference value between the maximum condition number and the minimum condition number of the system matrix W under each system parameter, and sequencing the corresponding system parameters according to the sequence of the difference values from large to small;
(4) and respectively and sequentially carrying out iterative optimization on each ellipsometer system parameter according to the sequence until the difference value between the minimum condition value of the most-front system parameter and the minimum condition value corresponding to the previous optimization is smaller than a threshold value, and finishing the optimization, wherein each corresponding system parameter is the optimized ellipsometer system parameter.
As a further preferred embodiment of the present invention, the specific process of iteratively optimizing the ellipsometer system parameters is as follows:
(1) for any ellipsometer system parameter, taking the current values of the rest system parameters as fixed values, and obtaining the relation between the system matrix W and the ellipsometer system parameter, so as to obtain the minimum value of the condition number of the system matrix W and the system parameter value corresponding to the minimum value, and take the minimum value and the system parameter value as the updated system parameter after one iteration;
in the same way, the updated values of the rest system parameters after one iteration are sequentially obtained;
(2) and repeating the process for multiple times, wherein each time the process is repeated, a group of updated ellipsometer system parameters are obtained, so that multiple iterative updates of each system parameter are completed until an iteration termination condition is met, and the latest group of ellipsometer system parameters can be obtained.
In the invention, the system model is a linear equation shown in formula (1):
wherein:
Αq=(aq,0aq,1aq,2aq,3)(2)
wq,j,k=aq,jsq,k(5)
wherein, IqIs the light intensity signal of the detector, AqTo analyze the arm vector, SqFor the deflection arm vector, subscript Q =0,1, …, Q-1 indicates the qth measurement component of the instrument, Q is the number of measurements, and M is the Mueller matrix of the sample to be measured, where the elements are represented by Mj,kIs shown as aq,j,sq,kRespectively correspond to AqAnd SqElement in vector, wq,j,kRepresenting the measurement component I of the q-th orderqElement m to be measured for j row and k columnj,kI.e. the transfer characteristic of the ellipsometer measurement system, the subscripts j =0,1,2,3 and k =0,1,2,3 denote the fourth element in the vector, respectively. The system model of the ellipsometer is shown as formula (1).
The form of changing the sample matrix M to be tested of 4 × 4 in the formula (1) into the column vector of 16 × 1After finishing, the method can obtain:
and all of the Q measurements may be represented as a Q × 16 ellipsometer system matrix W, where W isqRepresenting the Q-th row of the matrix W, the Q measurements of the detector can be represented as a measurement vector I, where I isqIs the q-th line of vector I. Equation (6) can be expressed in the form of a vector dot product:
the ellipsometer system matrix W in the formula (7) includes system parameters to be optimized for the ellipsometer, and the system matrix W includes different parameters for ellipsometers with different structure types, and for convenience of the following description, only three variables l, p, and m are used to represent the three parameters of the system matrix W, and actually the types and the number of the parameters vary according to different kinds of the ellipsometers.
In the invention, the condition number of the ellipsometer system matrix W under different ellipsometer system parameters is calculated, the condition number of the ellipsometer system matrix W is used as an evaluation standard for ellipsometer system parameter optimization, and the smaller the condition number is, the smaller the deviation in the measurement vector I is to the solution vector IThe smaller the effect of (c).
The condition number is a measure to mathematically evaluate the sensitivity of the solution of the linear system of equations Ax = b to errors or uncertainties in b. Then for equation (8), if the ellipsometer system matrix W condition number is large, a small deviation in the measurement vector I will cause a solution vectorLarge deviation and poor numerical stability, and if the condition number of the measurement matrix W of the ellipsometer system is small, the solution vector caused by the small deviation in the measurement vector IThe deviation of (2) is small and the numerical stability is good.
Because each ellipsometer system parameter has an influence on the W condition number, only one parameter is needed to be used as a variable quantity in the W condition number calculation process, and the corresponding relation between the condition number of the matrix W and the ellipsometer system parameter can be obtained by fixing the rest parameters.
The condition number of the ellipsometer system matrix W is shown in formula (9)
κp(W)=||W||p·||W-1||p(8)
W-1=(WT·W)-1·WT(10)
Wherein, κp(W) is the condition number of the ellipsometer system matrix W, | | W | | caltroppIs the p-norm of the ellipsometer system matrix W, equation (9) is the calculation method of the p-norm of the matrix W, D (W) is the space domain of the matrix W, x is an arbitrary vector in the space domain of the matrix W, xiDenotes the ith element, W, in the vector x-1The superscript T represents the matrix transpose, which is the generalized inverse of the matrix W.
In the invention, the difference value between the maximum value and the minimum value of the ellipsometer system matrix W condition number under each ellipsometer system parameter is calculated, the larger the difference value is, the larger the influence of the ellipsometer system parameter on the matrix W condition number is, the ellipsometer system parameter is used as a priority optimization object, and each ellipsometer system parameter is sequentially optimized according to the size sequence of the difference value. If the condition number difference of the matrix W corresponding to l in the three parameters of l, m and p is the maximum, and the condition difference of the matrix W corresponding to p is the minimum, the method is arranged as follows: l, m, p.
In the invention, according to the arrangement sequence of the parameters of the ellipsometer system in the third step, the parameters are optimized in sequence. In the optimization process, only one parameter is used as a variable quantity each time, and the numerical values of the rest fixed parameters all use the current optimal value as an initial value.
Step 1, respectively obtaining the values of ellipsometer system parameters l, m and p at the minimum condition number of an ellipsometer system matrix W and the corresponding minimum condition values: l1,κp(W(l1));p1,κp(W(p1));m1,κp(W(m1))。
Step 2, with p1,m1As an initial value, the relation between the matrix W condition number and the ellipsometer system parameter value l can be obtained, and the value of the matrix W condition number at the minimum position and the corresponding minimum condition value l are obtained2,κp(W(l2)). Then use2,m1The ellipsometer system parameter p is optimized as an initial value, and the corresponding p can be obtained2,κp(W(p2)). Also will l2,p2As initial value, corresponding m is obtained2,κp(W(m2)). The parameter value l of the ellipsometer system after the second optimization can be obtained2,p2,m2。
And 3, repeating the optimization process of the step 2 to obtain a new group of ellipsometer system parameters, and obtaining a group of ellipsometer system parameter values and corresponding condition values of the matrix W after nth optimization: ln,κp(W(ln));pn,κp(W(pn));mn,κp(W(mn))。
Now choose the ellipsometer system parameter l with the largest influence on the condition number of the matrix W, let Δ κp(W(l))=κp(W(ln))-κp(W(ln-1) When Δ κ)p(W(l))<The optimization is stopped as soon as possible. The ellipsometer system parameter value obtained after the nth optimization is the optimal system parameter of the current ellipsometer, wherein different values can be selected according to the optimization requirements as the judgment condition for stopping optimization.
Compared with a general instrument optimization method, the ellipsometer system parameter optimization method provided by the invention obtains the matrix W containing the ellipsometer system parameters by establishing an ellipsometer system model and expressing the ellipsometer system model in the form of vector dot product. And calculating the condition number of the matrix W relative to each parameter, and taking the condition number as the evaluation standard of the optimized configuration of the ellipsometer system. And calculating the difference value between the maximum value and the minimum value of the matrix W condition number under each ellipsometer system parameter again, and finally optimizing each ellipsometer system parameter from large to small according to the sequence of the difference values, wherein the optimal system parameter and the repeated calculation times are determined by the design requirements. The method can optimize all system parameters of the ellipsometer, obtain the optimal configuration of the ellipsometer system through repeated calculation for many times, and the optimization calculation times are variable according to different precision requirements, so that the method can be suitable for optimizing the ellipsometers with different structural types.
Drawings
FIG. 1 is a flow chart of a method for ellipsometer system parameter optimization provided in the present invention.
Fig. 2 is a schematic structural diagram of a dual-rotation compensator type muller matrix ellipsometer system.
Fig. 3 is a graph of a simulation result of a measurement of a mueller matrix ellipsometer of a dual rotation compensator type.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings. The description herein is intended to be illustrative only and not limiting with respect to the specific embodiments described.
The ellipsometer system parameter optimization method provided by the invention can be suitable for ellipsometers with different structure types, such as: rotary polarizer type, rotary analyzer type, single rotary compensator type, and dual rotary compensator type. The dual-rotation compensator type muller matrix ellipsometer system is the most complex in structure, and the system parameters to be optimized are the most. The principle and operation of the method of the present invention will be described in detail below with reference to the accompanying drawings by taking a dual-rotation compensator type muller matrix ellipsometer system as an example.
The implementation flow of the ellipsometer system parameter optimization method provided in this embodiment is shown in fig. 1.
(1) And establishing a system model of the ellipsometer to be optimized, and arranging the system model of the ellipsometer into a form of a vector dot product and a vector dot product to obtain a matrix W containing the system parameters of the ellipsometer.
For a dual-rotation compensator type muller matrix ellipsometer system structure, the system structure is shown in fig. 2. The dual-rotation compensator type muller matrix ellipsometer system mainly comprises a light source 1, a polarizer 2, a polarizing arm rotation compensator 3, a polarization analyzing arm rotation compensator 5, a polarization analyzer 6 and a detector 7. The polarizer 2 and the polarizing arm rotating compensator 3 form a polarizing arm 8 of the double-rotating-compensator type muller matrix ellipsometer measuring system, and the polarization analyzing arm rotating compensator 5 and the polarization analyzer form an polarization analyzing arm 9 of the double-rotating-compensator type muller matrix ellipsometer measuring system. The light beam emitted by the light source 1 is unpolarized light, is changed into linearly polarized light through the polarizer 2, is changed into elliptically polarized light through the modulation of the polarizing arm rotating compensator 3, then the polarized light is acted with the sample 4 to be detected, the polarization state of the polarized light is changed, the polarization detecting arm rotating compensator 5 demodulates the polarized light, and finally the light intensity signal of the emergent polarized light is detected through the polarization detector 6 and the detector 7.
The detector 7 detects to obtain a light intensity signal containing information of the sample 4 to be detected, then Fourier analysis is carried out on the light intensity signal detected by the detector 7 to obtain Fourier coefficients containing the Mueller matrix information of the sample 4 to be detected, and finally the Mueller matrix information of the sample 4 to be detected can be solved by the Fourier coefficients.
Wherein the polarizer 2 has a bias angle P relative to the incident surface, and the polarization arm rotation compensator 3 has a phase retardation1And an initial phase angle C with respect to the incident surfaces1Amount of phase delay of the resolver arm rotation compensator 42And an initial phase angle C with respect to the incident surfaces2The analyzer 3 is relative to the incident surface angle a, the number of sampling points N in a single period in two measurements, and the like, which are all ellipsometer system parameters to be optimized in the optimization of the double-rotation compensator type ellipsometer.
And the polarizing arm rotating compensator 3 and the polarization analyzing arm rotating compensator 5 of the double-rotating compensator type muller matrix ellipsometer measuring system continuously rotate in a certain proportion, wherein the proportion is 5: 3, the light intensity signal detected by the detector 7 in this case can be written as an expression shown in equation (11) (R.W.Collins et al., J.Opt.Soc.Am.A, Vol.16, pp.1997-2006,1999),
wherein I (t) represents the emergent light intensity of the ellipsometer measuring system, { I0,I0(α2n,β2n) Expressing the Fourier coefficients of each order of the light intensity signal when the light intensity signal is spread into a Fourier expansion form, wherein the Fourier coefficients contain the Mueller matrix information of the sample 4 to be tested, I0Is a direct current component, (α)2n,β2n) To be a direct current component I0Normalized Fourier coefficients of various orders, n represents the order of the Fourier coefficients, omega is the fundamental frequency of rotation of the two rotation compensators, and t is the rotation time of the rotation compensators.
Hadamard analysis of formula (11) gave the form shown in formula (12) (R.W.Collins et al., J.Opt.Soc.Am.A., Vol.16, pp.1997-2006,1999),
(12)
in equation (12), N represents the total number of integration points in which the number of sampling points in a single period is a Hadamard component, N should be greater than the total number of fourier coefficients, IqIndicating that the detector intensity in the qth measurement is the qth Hadamard integral component.
The relation between the Fourier coefficient and the Mueller matrix of the sample to be detected is shown as a formula (13),
in the formula (13), i represents a row of the mueller matrix of the sample to be measured, j represents a column of the mueller matrix of the sample to be measured, and mijThe element of the ith row and the jth column in the Mueller matrix of the sample to be tested is shown,andrespectively represent α2nAnd β2nAbout mijWhen the dual-rotation compensator type Mueller matrix ellipsometer measuring system is calibrated, the coefficientAndare all known quantities.
Thus, fourier coefficients of respective orders can be obtained by equation (12), and then a mueller matrix of the sample can be obtained by equation (13). Therefore, equations (11) to (13) are the system model of the dual-rotation compensator type muller matrix ellipsometer, and the transmission characteristics of the dual-rotation compensator type muller matrix ellipsometer system expressed by equations (11) to (13) conform to the linear relationship described in equation (1). Finally, after the formula (13) is substituted for the formula (12), the formula is expressed in a vector dot product form:
the formula (14) is a system model representing the dual-rotation compensator type muller matrix ellipsometer in the form of vector dot product, and obtains a matrix W containing parameters of the system to be optimized of the dual-rotation compensator type ellipsometer, wherein P, a,1,2and N is an ellipsometer system parameter to be optimized included in the system matrix W.
(2) And (3) solving the condition number of the matrix W of the dual-rotation compensator type ellipsometer system in the step (1).
The method for solving the p-norm of the ellipsometer system matrix W is defined in equations (9) to (11), and it is preferable to use the 2-norm of the ellipsometer system matrix W to find the condition number in this embodiment, as shown in equation (15):
κ2(W)=||W||2·||W-1||2(15)
wherein, κ2(W) is the corresponding condition number under 2 norms of the ellipsometer system matrix W, | W | | count2Is the 2 norm of the matrix W.
Because each ellipsometer system parameter has an influence on the W condition number, only one parameter is needed to be used as a variable quantity in the W condition number calculation process, and the corresponding relation between the condition number of the matrix W and the ellipsometer system parameter can be obtained by fixing the rest parameters. The parameters P of the dual-rotation compensator ellipsometer system are optimized in turn,1,2a, N. Firstly, taking the bias angle P of the polarizer 2 relative to the incident surface as a variable, taking one of the other 4 parameters to be in accordance with an actual arbitrary fixed value, namely obtaining the relation that the value of the condition number of the matrix W changes along with the value of the P value, and obtaining the P value at the minimum position of the condition number of the matrix W at the moment and the corresponding condition value: p1,κ2(W(P1)). The same method is adopted to obtain the relation between other parameters and the condition number of the matrix W, and the parameter value and the corresponding condition value at the minimum position of the condition number of the matrix W are obtained:1,1,κ2(W(1,1)),2,1,κ2(W(2,1)),A1,κ2(W(A1))。
(3) and (3) calculating the difference between the maximum value and the minimum value of the matrix W condition number of each parameter by using the relation between the parameters of the ellipsometer system of each double-rotation compensator obtained in the step (2) and the matrix W condition number, and arranging the parameters from large to small according to the difference. A larger difference indicates a larger influence of the parameter on the condition number of the matrix W, and priority needs to be given in optimization. The system parameters of the double-rotation compensator type ellipsometer are arranged from large to small according to the difference value of matrix W condition numbers corresponding to the parameters:1,2,P,A,N。
(4) and (4) sequentially optimizing each parameter according to the arrangement sequence of the parameters of the double-rotation compensator type ellipsometer system in the step (3).
Step 1, firstly, optimizing the phase delay amount of the deflection arm rotation compensator 31Will be2,1,P1,A1,N1The system parameters are used as initial values to obtain the condition number and1obtaining the minimum value of the condition number of the matrix W and the corresponding system parameter value:1,2,κ2(W(1,2))。
will be provided with1,2,P1,A1,N1The system parameters are used as initial values to obtain the condition number of the matrix W and the phase delay of the analyzer arm rotation compensator 52The minimum value of the condition number of the matrix W and the corresponding phase delay amount are obtained according to the relation between the matrix W and the condition number:2,2,κ2(W(2,2))。
will be provided with1,2,2,2,A1,N1And (3) taking the system parameters as initial values to obtain the relation between the matrix W condition number and the bias angle P of the polarizer 2 relative to the incident plane, and obtaining the minimum value of the matrix W condition number and the corresponding system parameter values: p2,κ2(W(P2))。
Will be provided with1,2,2,2,P1,N1And (3) taking the system parameters as initial values to obtain the relation between the matrix W condition number and the offset angle A of the analyzer 6 relative to the incident surface, and obtaining the minimum value of the matrix W condition number and the corresponding system parameter values: a. the2,κ2(W(A2))。
Will be provided with1,2,2,2,P2,A2Taking the system parameters as initial values, obtaining the relation between the matrix W condition number and the single-period sampling point number N, and obtaining the minimum value of the matrix W condition number and the corresponding system parameter values: n is a radical of2,κ2(W(N2))。
Thus obtaining a group of optimized parameters of the dual-rotation compensator type Mueller matrix ellipsometer system:1,2,2,2,P2,A2,N2。
and 2, repeating the optimization process of the step 1 for n times to obtain the n-th group of optimized parameters of the double-rotation compensator type ellipsometer system and the minimum condition number of the matrix W corresponding to each parameter:1,n,κ2(W(1,n)),2,n,κ2(W(2,n)),Pn,κ2(W(Pn)),An,κ2(W(An)),Nn,κ2(W(Nn)). Selecting ellipsometer system parameters having the greatest influence on matrix W condition number1The system parameter, Δ κ (W: (W) (W))1))=κ2(W(1,n))-κ2(W(1,n-1) When Δ κ (W: (W) ()1))<The optimization is stopped, wherein the termination conditions for the optimization process are determined according to the actual requirements.
Thus, the parameter values of the optimized double-rotation compensator type muller matrix ellipsometer system in the nth group can be considered as follows:1,n,2,n,Pn,An,Nni.e. the optimal system parameters meeting the current optimization requirements.
Fig. 3 is a simulation result diagram of simulation measurement of a dual-rotation compensator type muller matrix ellipsometer. In simulation, a sample to be tested is ideal air, and the Mueller matrix of the sample is as follows:
in the simulation, the ideal light intensity signal I is measuredqAdding a mean value of 0 and a standard deviation of 0.01IqAfter one gaussian noise is generated, the mueller matrix of the ideal air is solved through equations (12) to (14). FIG. 3 shows the results of the solution under two different ellipsometer system parameters, wherein m is the value00The constant value of 1 for the normalized element is not shown. The system parameter values of the unoptimized dual-rotation compensator type muller matrix ellipsometer are as follows:1=90°,2=90 °, P =10 °, a =10 °, N = 50. The optimized ellipsometer system parameter values are as follows:1=124°,2=124 °, P =45 °, a =45 °, N =70, and it can be seen from the figure that, for the same measurement noise, the simulation result under the optimized ellipsometer system parameter configuration condition fluctuates less than the unoptimized simulation result and is closer to the mueller matrix value of the ideal air. The optimized system parameters of the dual-rotation compensator type Mueller matrix ellipsometer obtained by the simulation are not the optimal system parameters of the structure type, and are the results after two times of optimization, and different optimization results can be obtained in the actual optimization process according to different optimization standards.
The optimization process of the ellipsometer system parameters of other different structure types is similar to the implementation process of the optimization of the dual-rotation compensator type muller matrix ellipsometer, and only a specific system model needs to be changed into a system model suitable for the ellipsometer of a specific structure type.
The above description is an embodiment of the present invention, but the present invention should not be limited to the disclosure of the embodiment and the drawings. Therefore, it is intended that all equivalents and modifications which do not depart from the spirit of the invention disclosed herein are deemed to be within the scope of the invention.
Claims (4)
1. A method for optimizing ellipsometer system parameters for determining optimal system parameters for an ellipsometer, the method comprising:
(1) establishing a system model of the ellipsometer, and arranging the ellipsometer system model into a vector and vector dot product form, so as to obtain a system matrix W containing ellipsometer system parameters;
(2) calculating the condition number of the system matrix W under each ellipsometer system parameter;
(3) calculating the difference value between the maximum condition number and the minimum condition number of the system matrix W under each system parameter, and sequencing the corresponding system parameters according to the sequence of the difference values from large to small;
(4) and respectively and sequentially performing loop iteration optimization on each system parameter of the ellipsometer according to the sequence until the difference value between the minimum condition value of the most-front system parameter and the minimum condition value of the system parameter optimized at the previous time is less than the threshold value, thereby completing the optimization, and at the moment, each corresponding system parameter is the optimized ellipsometer system parameter.
2. The ellipsometer system parameter optimizing method according to claim 1, wherein the specific process of performing the iterative optimization on the ellipsometer system parameter is as follows:
(1) for any ellipsometer system parameter, taking the current values of the rest system parameters as fixed values, obtaining the relation between the system matrix W and the ellipsometer system parameter, and further obtaining the minimum value of the condition number of the system matrix W and the system parameter value corresponding to the minimum value, so as to serve as the updated system parameter after one iteration;
in the same way, the updated values of the rest system parameters after one iteration are sequentially obtained;
(2) repeating the above process for multiple times, wherein each time the process is repeated, a group of updated ellipsometer system parameters is obtained, so that multiple iterative updates of each system parameter are performed, and the latest group of ellipsometer system parameters, namely the optimized ellipsometer system parameters, can be obtained until the iteration termination condition is met.
3. The method of claim 1 or 2, wherein the ellipsometer system matrix W has a condition number of:
κp(W)=||W||p·||W-1||p
wherein, W-1=(WT·W)-1·WT,κp(W) is the condition number of the ellipsometer system matrix W, | | W | | caltroppIs the p-norm of the ellipsometer system matrix W, D (W) is the spatial domain of the matrix W, x is an arbitrary vector in the spatial domain of the matrix W, xiDenotes the ith element, W, in the vector x-1The superscript T represents the matrix transpose, which is the generalized inverse of the matrix W.
4. The ellipsometer system parameter optimization method according to claim 1 or 2, wherein the system model is:
wherein, IqIs the light intensity signal of the detector, AqTo analyze the arm vector, SqFor a polarizing arm vector, subscript Q is 0,1, …, Q-1 represents the qth measurement component of the instrument, Q is the measurement times, M is the Mueller matrix of the sample to be measured, wherein the elements are represented by Mj,kIs shown as aq,j,sq,kRespectively correspond to AqAnd SqElement in vector, wq,j,kRepresenting the measurement component I of the q-th orderqElement m to be measured for j row and k columnj,kThe index j is 0,1,2,3 and k is 0,1,2,3 respectively represent the number of elements in the vector;
any element W in the system matrix WqObtained by the following formula:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310302341.XA CN103426031B (en) | 2013-07-18 | 2013-07-18 | A kind of optimization method of ellipsometer system parameter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310302341.XA CN103426031B (en) | 2013-07-18 | 2013-07-18 | A kind of optimization method of ellipsometer system parameter |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103426031A CN103426031A (en) | 2013-12-04 |
CN103426031B true CN103426031B (en) | 2016-04-13 |
Family
ID=49650737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310302341.XA Active CN103426031B (en) | 2013-07-18 | 2013-07-18 | A kind of optimization method of ellipsometer system parameter |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103426031B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103776445B (en) * | 2014-02-24 | 2017-01-04 | 北京理工大学 | Amplitude-division polarization navigation angle sensing design method and device |
CN112345460B (en) * | 2020-10-14 | 2024-03-22 | 武汉颐光科技有限公司 | Double-rotation and high-yield mode switching method and system for Mueller matrix ellipsometer |
CN113281267B (en) * | 2021-05-14 | 2022-05-20 | 华中科技大学 | Method for calibrating system parameters of dual-rotation compensator type Mueller matrix ellipsometer |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520413A (en) * | 2009-02-18 | 2009-09-02 | 深圳大学 | A heterodyne interference elliptic-deviation measurement nonlinear error compensation method |
CN101701820A (en) * | 2009-11-02 | 2010-05-05 | 北京航空航天大学 | Method for extracting optical fiber gyro random error characteristics based on wavelet variance |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7277172B2 (en) * | 2005-06-06 | 2007-10-02 | Kla-Tencor Technologies, Corporation | Measuring overlay and profile asymmetry using symmetric and anti-symmetric scatterometry signals |
-
2013
- 2013-07-18 CN CN201310302341.XA patent/CN103426031B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520413A (en) * | 2009-02-18 | 2009-09-02 | 深圳大学 | A heterodyne interference elliptic-deviation measurement nonlinear error compensation method |
CN101701820A (en) * | 2009-11-02 | 2010-05-05 | 北京航空航天大学 | Method for extracting optical fiber gyro random error characteristics based on wavelet variance |
Non-Patent Citations (2)
Title |
---|
基于修正等效介质理论的微纳深沟槽结构反射率快速算法研究;刘世元 等;《物理学报》;20080930;第57卷(第9期);5996-6001 * |
基于非线性最小二乘拟合法的Muller矩阵椭偏仪;侯俊峰 等;《中国激光》;20130430;第40卷(第4期);1-8 * |
Also Published As
Publication number | Publication date |
---|---|
CN103426031A (en) | 2013-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103163077B (en) | Rotary device type spectroscopic ellipsometers systematic parameter calibration steps | |
CN111122460B (en) | Single-rotation compensator type spectroscopic ellipsometer parameter calibration method and device | |
CN107490397B (en) | High-accuracy self-adaptation filters the quick Peak Search Method of FBG spectrum | |
KR101888722B1 (en) | Numerical aperture integration in raleigh wavelengths for optical critical dimension(ocd) metrology | |
CN110793556B (en) | Distributed three-dimensional shape sensing demodulation method based on optical frequency domain reflection parameter optimization | |
CN102080988A (en) | Device and method for detecting single photon polarization quantum state in real time | |
CN103426031B (en) | A kind of optimization method of ellipsometer system parameter | |
CN102980748A (en) | Method for measuring optical constants of thin film with non-uniform refractive index | |
CN106643835B (en) | Optical fiber Fabry-Perot cavity demodulation method and device and optical fiber Fabry-Perot interferometer | |
CN115391726B (en) | Fitting optimization method and related device for ellipsometry system | |
CN104568248A (en) | Measuring method of microzone stress of amorphous optical thin film | |
Harasim et al. | Improvement of FBG peak wavelength demodulation using digital signal processing algorithms | |
CN111121653A (en) | Single-layer film critical thickness estimation value calculation method | |
CN102519597A (en) | Phase correction apodization method for Fourier transform spectrometer | |
CN103217385B (en) | Method for assessing random error of ellipsometer measurement system | |
CA2927376A1 (en) | Method and system for tracking the centre of a peak from a plurality of sample points in an optical system | |
CN103411890A (en) | Rotating compensator type ellipsometer system error assessment and elimination method thereof | |
CN104677837A (en) | Calibration method of full-Mueller matrix ellipsometer | |
CN105004707A (en) | CCD array on-line Raman spectrometer spectra standardization method | |
CN110567883B (en) | Method for measuring Mueller matrix spectrum in real time | |
EP2306166B1 (en) | Measurement of fourier coefficients using integrating photometric detector | |
JP2015105909A (en) | Ofdr optical fiber measuring method using group delay calculation and device implementing the same | |
CN107643135B (en) | Stress analysis method based on time domain phase unwrapping | |
Harasim et al. | Accuracy improvement of Fiber Bragg Grating peak wavelength demodulation using wavelet transform and various center wavelength detection algorithms | |
Varentsova et al. | Restoration of a signal and its instantaneous spectral characteristics by the sliding window method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |