1. Introduction
The detection and identification of surface defects is a key aspect of the quality evaluation of optical elements. With the wide application of precision optical instruments such as ultra-smooth optical elements [
1], there is an increasing demand for the detection and identification of surface defects at the nanometer scale [
2]. For example, in the inertial confinement fusion system (ICF), the aperture of optical elements can reach about 1 m, while the measurement accuracy of surface defects is required to reach the micro-nanometer level. Surface defect detection includes not only the characterization of surface morphology, surface shape errors, and fluctuations in roughness but also the identification of defects such as surface pockmarks and subsurface damage. Surface defects can be categorized as pockmarks, SSBDs, and dust particles, among others. The different types of damage and the corresponding repair techniques vary greatly; so, it is important to detect and classify various types of component defects [
3].
It has been common to use manual visual inspection to evaluate the quality of ultra-precision optical components [
4]. However, this method relies on manual subjective judgment, and the evaluation results are unstable and inefficient. Therefore, more and more academics are devoted to the research of inspection methods that can realize the quality evaluation of optical components, such as photothermography [
5], optical interferometry [
6], optical coherence tomography [
7], and microscopic dark-field scattering imaging [
8]. The microscopic dark-field scattering imaging method is currently the main surface defect detection method and uses a specific device to collect and image the scattered light induced by surface defects. However, the method only makes use of the radiation intensity and geometric properties of electromagnetic waves and is based on limited information. It is therefore difficult to distinguish dirty particles and pockmarks that appear as dots on dark-field intensity images [
1,
9]. Since different types of microscopic defects on the surface of optical elements affect the polarization state of light to different degrees, Thomas et al. [
10,
11] added polarization information to the original scattered light intensity information to analyze the defects existing on the surface of a sample. Bickel and Wang [
12] used ray tracing to find the scattered field distribution based on the theory of polarization and put forward the idea of the Mueller’s matrix. Letnes et al. [
13] used p- and s-polarized light to establish a model of the polarization effect and used all the elements in the Mueller matrix to characterize the polarization properties of the scattered light from a damaged metal surface with roughness in the order of micrometers. However, the accurate and fast detection and classification of various types of nanoscale defects on the surface of optical elements is still a great challenge.
The BRDF [
10] is a commonly used function for modeling and investigating the surface scattering properties of optical elements and is defined as the ratio of the irradiance exiting through the target surface to the irradiance incident on the target surface. The polarized bidirectional reflection distribution function (pBRDF) is generalized to the polarization case by means of the Mueller matrix based on the BRDF [
14]. In the field of defect detection, pBRDF can consider the spatial scattering distribution of various defects on the surface of an optical element more comprehensively than the BRDF. It is able not only to completely describe the light scattering distribution of defects on the surface of the element, but also to characterize the light scattering polarization properties of the target defects. Therefore, the pBRDF has become an important tool for obtaining polarization detection information of target defects [
15]. When the incident and reflected light are described using Jones vectors, the pBRDF becomes a Mueller matrix characterizing the polarization transmission of the optical element, according to the transformation relationship between the Jones matrix and the Mueller matrix [
3].
Based on the principle that different physical properties of the surface defects of optical elements lead to a different modulation of the polarization state of the reflected light, this paper proposes a new method for the optical identification of surface roughness [
14,
16], pockmarks, dirty particles above the surface, and SSBDs [
17] in optical elements using RMMS. In this method, firstly, the combination of the BRDF and the Rayleigh–Rice vector scattering theory is innovatively proposed to construct the pBRDF model for various types of defects on the surface of optical elements, and then the BRDF is numerically analyzed to obtain the optimal detection conditions for defects on the surface of optical elements. On this basis, surface roughness and pockmarks on the surface of optical elements made of borosilicate at the nanometer level were simulated and analyzed via RMMS, and the accuracy of the simulation results was verified through spectroscopic ellipsometry experiments. Then, dirty particles and pockmarks, as well as SSBD, on the surface of the optical element were simulated and analyzed via RMMS [
13]. Based on the RMMS of various types of defects on the surface of optical elements, the detection and classification of pockmarks, roughness, dirty particles, and SSBD on the surface of optical elements were realized [
10,
18,
19]. Finally, the analysis of the one-to-one correspondence between the RMMS measurements and the size of various types of defects allowed the calculation of the dimensions of various types of defects utilizing Muller’s matrix elements [
17].
2. Theory
When a defect on the surface of a light element is illuminated, it displays various reflection characteristics based on the angles of incidence, zenith, and azimuth. This characteristic can be described using the BRDF [
20], which is defined as the ratio of the irradiance
exiting the defect surface in direction
to the irradiance
incident upon the defect surface in direction
, i.e.,
where the incident angle is represented by
, the scattering angle by
, the incident azimuth by
, and the scattering azimuth by
. A schematic of the BRDF scattering geometry is shown in
Figure 1 [
14].
The Rayleigh–Rice theory [
21] generalized the first-order vector perturbation theory. It is a scattering model based on physical optics and applies to surface scattering cases where the magnitude of the variation in the surface height profile is less than the order of the wavelength of the visible light. The Rayleigh–Rice theory deconstructs the surface profile microrough topography into numerous sinusoidal gratings with different magnitudes, spatial periods, and orientations. It then considers the scattered field of the surface as the height of each of these sinusoidal grating components. To determine the intensity of the scattered field, it is necessary to know the height of each sinusoidal grating component of the surface height profile. The power spectral density
is used, which represents the square of the average height of the surface sinusoidal grating components with spatial frequency
, and is defined as follows:
where
represents the area of the scatterer, while
is the surface height profile function of the scatterer, which represents the deviation of the heights of the points on its surface from the average surface height. Thus,
contains the complete dependency of the scattered field on the surface topography.
2.1. Scattering Model for Dig Defects
Digs refer to pits and blemishes on the surface of an optical element, as shown in
Figure 2, which gives a composite scattering schematic of digs. For tiny-sized dig defects with depths smaller than the wavelength of visible light, often between tens to hundreds of nanometers, the Rayleigh–Rice theory effectively describes the scattering characteristics of such surface features.
According to the electromagnetic theory of physical optics, the relationship between the corresponding
s- and
p-components
and
of the incident field and the
s- and
p-components
and
of the scattered field can be established through the Jones matrix [
3]:
The following transformation relationship exists between the Mueller and Jones matrices [
22]:
The expression for the scattering characteristics of the Mueller matrix (pBRDF) due to the modulation of the incident optical field by the surface digs of the optical element is given by the following equation:
the pBRDF comprises 16 different elements, each with a distinct meaning. The value of
M11 reflects the target’s ability to transmit, scatter, and reflect incident light.
M12,
M13, and
M14 reflect the target’s ability to attenuate incident light horizontally, vertically, and in a circular bidirectional manner, respectively.
M21,
M31, and
M41 reflect the target’s ability to polarize incident unpolarized light, and the remaining nine elements reflect the target’s ability for incident light depolarization and phase delay [
15].
2.2. Scattering Model for Microrough Surfaces
The microroughness of the surface of a precision optical element causes the surface profile to be smaller than the visible wavelength by an order of magnitude. The scattering properties of this surface topography can be described using the Rayleigh–Rice theory. In this theory, the calculation of the pBRDF depends on the power spectral density (PSD) of the surface. As shown in Equation (2), the PSD is a function of the height profile
of the scatterer surface, and since the microscopic morphology structure of the microrough surface is determined through stochastic processes with a certain average nature, the height profiles of the different regions of the surface are not the same, and therefore the scattering characteristics will fluctuate with a certain statistical law. Therefore, in this study, the K-correlation model function [
23] is considered to fit the overall average PSD of the microrough surface as follows:
This equation contains three adjustable parameters
A,
B, and
C, which are related to the statistical properties of the surface roughness
, as shown in the following equation:
where
denotes the correlation length of the microrough surface and the value of the parameter
C is 3. After determining the PSD of the surface roughness from the K-correlation model function, the pBRDF of the surface microroughness of the optical element can be calculated.
2.3. Scattering Model of Dirty Particles
Dirty particles above the surface of an optical element generally have equivalent diameters between tens of nanometers and tens of micrometers. The sample sizes of the dirty particles in the present study are smaller than the visible wavelengths. Therefore, these particles can be utilized to simplify the problem into a free-space particle scattering scenario using the two-interaction model [
24] as proposed by Germer [
10] and Nahm [
25], as shown in
Figure 3.
When the incident light irradiates the target surface at wavelength
and angle of incidence
, the pBRDF expression for a particle of radius
at a distance
above the smooth surface in the direction of the scattering zenith angle
and the relative azimuth angle
is as follows:
where
N/
A represents the density of particle scatterers in the irradiated region. The structure factor
F depends on the correlation between scatterers and equals 1 for random and uncorrelated examples.
2.4. Scattering Model for Subsurface Bubble Defects
The bubbles in the optical components can be regarded as spherical particles with a nanometer radius, while the surrounding medium can be considered a uniform background. Due to the distribution of bubbles within the components being random, the light scattered by the bubbles can be regarded as a random distribution of spherical particles interacting with light. If the size of the bubble particles is much smaller than the wavelength of the visible light, the Rayleigh scattering theory is still applicable for studying the characteristics of bubble light scattering, even when dealing with random particle distributions.
Figure 4 gives a schematic diagram of the composite scattering of bubble particles below the surface of the optical element.
When the incident light irradiates the target surface at wavelength
and angle of incidence
, the pBRDF expression for a particle of radius
at a distance
below the smooth surface in the direction of the scattering zenith angle
and the relative azimuth angle
is as follows:
where
is the refractive index of the surface element material,
represents the depth of light penetration in the material,
is the phase delay produced by the refractive action of the surface of the optical element on the incident field, and
is the phase delay produced by the refraction of the surface of the element on the scattered field.
In summary, by analyzing and calculating the pBRDF of the above four defects, the Mueller matrices of all types of microscopic defects on the surface of the optical element can be expressed as Equation (10). Among the 16 characteristic polarization detection elements,
M11 is normalized and has no significant value. The detect and classify pockmarks, roughness, dirty particles, and SSBD on the surface of the optical element on the nanoscale can be achieved using the elements of
M12 (
M21),
M33 (
M44), and
M34 (
M43).
3. Calculation of Optimal Testing Conditions
It is necessary to determine the optimal detection conditions before performing RMMS to detect and classify various types of defects on the surface of optical elements. This study uses dirty particles on the surface of optical elements and SSBDs as an example. For ultra-smooth optical elements, most of the dirty particles of interest to the surface defect detection system have a diameter of about tens of nanometers, the size of which is much smaller than the wavelength of visible light. The Rayleigh–Rice scattering model can be used to approximate the interaction of a single particulate dust with the incident light. The BRDF is used instead of the pBRDF to calculate the optimal experimental conditions, as in Equations (11) and (13), which effectively simplifies the tedious calculation process.
where
where
In this study, we take the dirty particles and SSBDs with a particle size of
a = 5 nm on the surface of the optical element as an example, and draw the pseudo-color maps with radial coordinates
, azimuthal angle
, and incidence angle of 0° and 70°, respectively, as shown in
Figure 5 and
Figure 6.
The different combinations of polarization states are indicated by the subscripts of BRDF, with the first one representing the polarization state of the incident light and the second one representing the polarization state of the outgoing light. The distributions of BRDFss, BRDFsp, and BRDFps (corresponding to (i), (ii), and (iii) in
Figure 5 and
Figure 6, respectively) for dirty particles above the surface and SSBD defects do not change significantly for different incidence angles, and all of them show symmetric distributions. On the other hand, BRDFpp (iv) exhibits an asymmetric distribution.
The reason for this, as can be seen from Equations (12) and (14), is that the three polarization factors Jss, Jsp, and Jps in the Jones matrix have the same function term , so the scattering of the two kinds of particles will be superimposed in the scattering component measurement, resulting in the pseudo-color maps of the three polarization factors, Jss, Jsp, and Jps, which are almost the same. In contrast, the polarization factor Jpp lacks the common functional term as described above, and its respective expression is significantly different. Therefore, it is theoretically possible to characterize the optimal detection conditions for dirty particles and SSBD above the surface of an optical element by using the relationship between the BRDFpp and the refractive index, the angle of incidence, the scattering angle, and the azimuthal angle.
Following this, the BRDFpp term was selected to solve for the optimal detection conditions [
18]. The simulation results of the optimal detection conditions are shown in
Figure 7 when the size of the dirty particles
a = 5 nm, the refractive index of the particles
= 1.457, and the refractive index of the material
= 1.5151, and
= 633 nm.
As
increases, the maximum BRDFpp corresponds to
, which gradually increases until the maximum value is obtained when
is about 70° and
, as shown in
Figure 7a. Based on
Figure 7b, the amplitude of BRDFpp at
is changed with an increase in particle size, but the trend of the scattering distribution remains unaffected, and the maximum BRDFpp value is reached at
. Similarly, analyzing surface roughness, pockmarks, and SSBD determined the optimal detection conditions for all types of defects on the optical element surface as
.
5. Conclusions
In this paper, based on the distinct reflection characteristics of pBRDF caused by the modulation of the incident light field by different kinds of nanoscale defects on the surface of the optical element, we proposed a method for detecting and categorizing nanoscale roughness, digs, dirty particles, and SSBDs on the surface of the optical elements via RMMS. Firstly, a numerical method was used to simulate and analyze the BRDFpp in order to obtain the optimal detection conditions for detecting and classifying various types of defects accurately on the surface of the optical element being tested. In addition, the pBRDF scattering models of surface roughness, digs, dirty particles, and SSBDs of optical elements were established, and the experimental results showed that the Mueller matrix elements M12 = M21, M33 = M44, and M34 = −M43 can significantly distinguish the various types of defects. The dimensions of the various types of defects on the surface can be calculated using the Mueller matrix elements. This study provides a theoretical and experimental foundation for the detection and classification of surface defects on ultra-precision optical components. The Mueller matrix contains rich polarization information, and in future work the polarization parameters can be separated and extracted by decomposing the transmission Mueller matrix to further analyze the internal structure of the optical element under test.