Investigation of Plasma Activated Si-Si Bonded Interface by Infrared Image Based on Combination of Spatial Domain and Morphology

This paper presents a detection method for characterizing the bonded interface of O2 plasma activated silicon wafer direct bonding. The images, obtained by infrared imaging system, were analyzed by the software based on spatial domain and morphology methods. The spatial domain processing methods, including median filtering and Laplace operator, were applied to achieve de-noising and contrast enhancement. With optimized parameters of sharpening operator patterns, disk size, binarization threshold, morphological parameter A and B, the void contours were clear and convenient for segmentation, and the bonding rate was accurately calculated. Furthermore, the void characteristics with different sizes and distributions were also analyzed, and the detailed statistics of the void’s number and size are given. Moreover, the orthogonal experiment was designed and analyzed, indicating that O2 flow has the greatest influence on the bonding rate in comparison with activated time and power. With the optimized process parameters of activated power of 150 W, O2 flow of 100 sccm and time of 120 s, the testing results show that the bonding rate can reach 94.51% and the bonding strength is 12.32 MPa.


Introduction
Wafer direct bonding is a promising technology emerging with the development of integrated circuits (ICs) and microelectromechanical systems (MEMS). It was originally used for micro-cavity sealing of sensor components and vacuum packaging of resonators or detectors [1]. Over the last few decades, wafer direct bonding technology has made significant advances, and played an irreplaceable role in the fabrication of silicon-on-insulator (SOI) and III-V compounds-on-insulator (such as GaSbOI, InPOI and InGaAsOI) substrates [2][3][4][5], fabrication of optoelectronic and photonic devices [6], wafer-level packaging [7], and three-dimensional (3D) integration [8,9]. However, silicon-based wafer direct bonding always needs a high annealing temperature (700-1000 • C) to obtain high bonding strength. Consequently, some severe problems are generated during this process, such as metalized layer and circuits remelting, dopant diffusion, structure fracture caused by mismatch of thermal expansion coefficients between dissimilar materials, and so on.
To reduce the bonding temperature, plasma activated bonding (PAB) and surface activated bonding (SAB) have been reported [10]. O 2 /N 2 /Ar/CF 4 reactive ion etching (RIE) plasma [11][12][13][14], N 2 microwave (MW) radicals [15], ultraviolet radiation [16], and Ar fast atom bombardment (FAB) [17,18] are utilized to activate the bonding surfaces, which make the surfaces hydrophilic and terminate with -OH. Then, these surfaces bond together due to the van der Waals attraction through hydrogen bonds, resulting in a hydrogen-bond interface structure connected by Si-OH-(H 2 O) x -HO-Si. Further thermal annealing at low temperature (200-400 • C) is required to obtain strong Si-O-Si bonds. However, bubbles and voids are formed at the bonded interface due to interfacial trapped hydrogen gas, contaminants, particles, and plasma-induced defects. These voids not only influence the bonding yield, but also decrease the bonding strength. Many researchers have studied the morphology of interfacial voids with different process parameters and attempted to suppress the void formation [19][20][21]. The generally used void detection methods are infrared imaging method [22], ultrasonic method [23], X-ray image method [24], and mirror image method [25]. The resolution of the ultrasonic method is limited by the thickness of the upper surface about the detected material. The X-ray image method is time-consuming, expensive, and the resolution is restricted by the resolution of the detector. The mirror image method requires thinning and polishing the bonding wafer to improve sensitivity. The infrared imaging method is a nondestructive testing technique. Although the resolution is less than with the ultrasonic method and X-ray image method, it is low cost, rapid, and can realize real-time detection. Therefore, the infrared imaging method is widely used for interface detection of silicon wafer bonding. However, the information of the infrared imaging is quite limited, so an accurate image processing technology is needed for further analysis of void characteristics, such as size and distribution. Xiao Tan et al. [26] have developed a precise, nondestructive visual testing method based on wavelet image analysis, which displayed the images with low background noise by using wavelet image denoising, wavelet image enhancement, and contrast enhancement. Nonetheless, wavelet threshold denoising would lead the image to lose some detailed information and produce visual distortion because of the shrinkage of soft threshold and the rough of hard threshold. Haijiao Yun et al. [27] have proposed a novel method combined with platform histogram equalization and fuzzy set theory to solve the low contrast and fuzzy image edge. This method overcame some shortcomings that traditional algorithms showed, such as the phenomenon of the excessive enhancement and weakening of the local gray information. Nevertheless, the decomposition and reconstruction of low or high frequency signal always makes the algorithm more complex. Therefore, a simple but effective image processing method with the functions of denoising and contrast enhancement is needed.
In this paper, a quality detection method of bonded interface for O 2 plasma assisted silicon wafer bonding based on spatial domain and morphology was developed and applied to the analysis of void formation and calculation of bonding rate. The infrared imaging system was established, and the infrared image was obtained. For obtaining the detailed information of interfacial voids, the infrared detection software with Graphical User Interface (GUI) was achieved using Matlab (R2015b) software. The simple and commonly used spatial domain processing methods, such as median filtering and Laplace operator were utilized to achieve denoising and contrast enhancement. Five parameters, sharpening operator patterns, disk size, binarization threshold, morphological parameter A and B, which have a significant influence on image processing quality, were investigate. The void contours in infrared images were clear and convenient for segmentation, and the bonding rate was accurately calculated by using this. In addition, the variation of the numbers and diameters of voids in the bonded interface were also analyzed, and the causes of the void formation and the influence for the bonding rate were explored. Finally, the orthogonal experiment was designed and analyzed, and the optimized process parameters were obtained.

Materials and Bonding Process
Commercially available Czochralski (CZ) grown double-side-polished 4 inch n-type (100) oriented Si wafers were used for experiments. The thickness of the Si wafer was 500 ± 20 µm and the resistivity was between 2 and 4 Ω·cm.  flow of 50-150 sccm. Thereafter the wafers were further cleaned with DIW in megasonic cleaner W-357HP for 15 min with the power of 450 W and the frequency of 1 MHz. Then, particles less than 0.3 µm were further removed so more water molecules could be adsorbed by half bonds and hydroxyl groups on the bonding surface. Prebonding was achieved in SUSS SB6e bonder. The tool pressure was set to be 3000 mbar, the chamber pressure was 5 × 10 −4 mbar and the prebonding time was 3 min. Finally, the bonded pairs were annealed at 350 • C for 2 h under nitrogen conditions for promoting dehydration condensation reaction between the Si wafers, and the stable covalent Si-O-Si bonds in the bonding interface were formed.

Infrared Imaging System
An infrared imaging system was used for detecting defects and their distribution of the bonded wafer interface before and after annealing. Figure 2 shows the structure diagram of infrared imaging system. The bonded wafers were placed on the object stage between the light source and the charge coupled device (CCD) camera, and a diaphragm was placed under the bonded wafers for high contrast image. The height position of the CCD camera was adjustable, so that the best shooting angle could be obtained. To get a clear infrared image, the focal length of the lens was adjusted manually. The infrared image was stored in the computer through the image acquisition card for further analysis. As the quality of the captured image is directly influenced by the performance of the CCD camera, the lens and their compatibility, an ultra-low illumination black and white camera WAT-902H2 SUPREME and megapixel prime lens M2514-MP2 were chosen. While testing, the nearinfrared portion of the light-wave passes through the Si wafer. If there are unbonded regions in the bonded interface, the light will be reflected twice, forming the coherent light. After the coherent light is processed by the CCD camera, some interference fringes appear in the image. If the unbonded area is large and the gap between the bonded wafers is small, many large interference fringes will occur in the image. When the unbonded area is small, the small Newton rings are observed. As the gap is large enough (for example particles), the transmittance of the infrared light decreases, and then the black patterns are observed in the image.
However, the infrared imaging is a relatively crude detection method. It can only provide a visual image, from which the information is quite limited. So it must be further processed to obtain detailed information of the bonded interface. Since the main factors that affect the quality of the infrared image are noise interference and sharpness, a combination of spatial domain and morphology method for image processing is used for denoising and edge clarity. As long as the noise is reduced and the sharpness is subsequently improved by image processing, a high-quality infrared image can be obtained which is convenient for investigating the morphology of voids and calculating the bonding rate.

Infrared Imaging System
An infrared imaging system was used for detecting defects and their distribution of the bonded wafer interface before and after annealing. Figure 2 shows the structure diagram of infrared imaging system. The bonded wafers were placed on the object stage between the light source and the charge coupled device (CCD) camera, and a diaphragm was placed under the bonded wafers for high contrast image. The height position of the CCD camera was adjustable, so that the best shooting angle could be obtained. To get a clear infrared image, the focal length of the lens was adjusted manually. The infrared image was stored in the computer through the image acquisition card for further analysis. As the quality of the captured image is directly influenced by the performance of the CCD camera, the lens and their compatibility, an ultra-low illumination black and white camera WAT-902H2 SUPREME and megapixel prime lens M2514-MP2 were chosen. While testing, the near-infrared portion of the light-wave passes through the Si wafer. If there are unbonded regions in the bonded interface, the light will be reflected twice, forming the coherent light. After the coherent light is processed by the CCD camera, some interference fringes appear in the image. If the unbonded area is large and the gap between the bonded wafers is small, many large interference fringes will occur in the image. When the unbonded area is small, the small Newton rings are observed. As the gap is large enough (for example particles), the transmittance of the infrared light decreases, and then the black patterns are observed in the image.
However, the infrared imaging is a relatively crude detection method. It can only provide a visual image, from which the information is quite limited. So it must be further processed to obtain detailed information of the bonded interface. Since the main factors that affect the quality of the infrared image are noise interference and sharpness, a combination of spatial domain and morphology method for image processing is used for denoising and edge clarity. As long as the noise is reduced and the sharpness is subsequently improved by image processing, a high-quality infrared image can be obtained which is convenient for investigating the morphology of voids and calculating the bonding rate.
The Matlab software developed by MathWorks, Inc. (Natick, MA, USA) is utilized to develop corresponding processing modules and design the GUI interface. For multi-image processing, the image effect in real time can be observed, which simplifies the processing operation. So, it is convenient to debug parameters and verify results. The GUI interface can be applied for all subsequent image processing and displayed in the part of software evaluation. The Matlab software developed by MathWorks, Inc. (Natick, MA, USA) is utilized to develop corresponding processing modules and design the GUI interface. For multi-image processing, the image effect in real time can be observed, which simplifies the processing operation. So, it is convenient to debug parameters and verify results. The GUI interface can be applied for all subsequent image processing and displayed in the part of software evaluation.
The image processing process mainly includes the following nine steps: 1. Load an infrared image. A square image is needed, which means that the number of pixel rows in the image must be equal with the number of pixel columns because a square matrix is required to convolve with the sharpened operator during the image enhancement process. 2. Background interference removing by using median filter in the spatial domain, which also called image denoising. As we know, the image acquired by CCD camera is severely affected by discrete impulse noise and salt-pepper noise. These kinds of noises are commonly derived from the electronic circuit and sensor caused by low illumination or high temperature, and can be effectively removed by median filter method. Suppose that {x , (i, j) ∈ I } denotes the gray value of each point of the digital image, and the two-dimensional median filtering window A can be defined as It is a nonlinear smoothing filtering method that removes noise while preserving the sharpness and detail of the image. By this way, most of the noise mixed in the digitized image is eliminated, which lays a good foundation for the subsequent detail information extraction.
3. Image enhancement by Laplacian operator. The Laplacian operator is suitable for refining the image blur caused by diffuse reflection of light. It is an isotropic second-order differential operator. For a continuous binary function f(x, y), the Laplacian at position (x, y) is delineated as The sharpening filter in the spatial domain can be expressed in convolution form as where H(r, s) is the sharpening operator pattern. The edges and contours of the image are generally located at the position where the grayscale is abrupt, that is, the high frequency part of the image. The edges are usually extracted by the grayscale difference, which is often in any directions. So the The image processing process mainly includes the following nine steps: 1.
Load an infrared image. A square image is needed, which means that the number of pixel rows in the image must be equal with the number of pixel columns because a square matrix is required to convolve with the sharpened operator during the image enhancement process.

2.
Background interference removing by using median filter in the spatial domain, which also called image denoising. As we know, the image acquired by CCD camera is severely affected by discrete impulse noise and salt-pepper noise. These kinds of noises are commonly derived from the electronic circuit and sensor caused by low illumination or high temperature, and can be effectively removed by median filter method. Suppose that x ij , (i, j) ∈ I 2 denotes the gray value of each point of the digital image, and the two-dimensional median filtering window A can be defined as It is a nonlinear smoothing filtering method that removes noise while preserving the sharpness and detail of the image. By this way, most of the noise mixed in the digitized image is eliminated, which lays a good foundation for the subsequent detail information extraction.

3.
Image enhancement by Laplacian operator. The Laplacian operator is suitable for refining the image blur caused by diffuse reflection of light. It is an isotropic second-order differential operator. For a continuous binary function f(x, y), the Laplacian at position (x, y) is delineated as The sharpening filter in the spatial domain can be expressed in convolution form as where H(r, s) is the sharpening operator pattern. The edges and contours of the image are generally located at the position where the grayscale is abrupt, that is, the high frequency part of the image. The edges are usually extracted by the grayscale difference, which is often in any directions.
So the differential operations need to be directional. The isotropic edge detection operator has the same detection capability for edge contours regardless of their directions. The method is to achieve image enhancement through two normalized high-pass filter patterns, which weakens or removes the low-frequency part of the image, meanwhile retains the high-frequency part of the image. Their equations are shown as below: The normalization method is to subtract the Laplacian operator from the original image, H 1 is a four-neighbor pattern, H 2 is an eight-neighbor pattern, and the discrete function expressions are respectively: This method appertains to a sort of linear sharpening filtering, which can avoid the brightness deviation of the processed image, and may restore the background data while ensuring the Laplacian sharpening effect, so that highlight the edge information well.

4.
Make use of the combination of Top-Hat and Bottom-Hat filtering method to enhance the contrast. Firstly, the enhanced image is added to the Top-Hat filtered image, and then the Bottom-Hat filtered image is subtracted. Finally, the filtered image is inverted to achieve the negative film effect. In this part, a value, called disk size, is set to adjust the influence of the contrast strength.
The disk size ranges from 1 to 100. In addition, it is also possible to suppress mixed noise composed by Gaussian noise and impulse noise in the infrared image. This method makes the outer boundary of the bonded wafers more obvious, and the texture of the interference fringe is more prominent.

5.
Invert the filtered image to achieve the negative effect. In general, uint8 represents 256-level grayscale in the image file, and the range of data is 0~255. The negative effect of the image is to subtract the original image data by 255. It is convenient to analyze the filtering situation compared with the image in the previous step. 6.
Set the binarization threshold to display the bonding area. This is the vital step for a better result. A large number of image tests show that the extent of threshold is regularly between 0.2 and 0.5. It is requested to conserve the complete information of the image without introducing excessive noise pixels. 7.
Properly perform the morphological processing. The purpose of the first morphological processing is to extract the main body. It requests the effective pixel points that are not bonded and filters out the imported noise pixel points after binarization. To obtain the information as complete as possible of each part of the image in the end. Need to pay attention that the processing times are Inf until the image no longer changes. The reason why the mode can be modified aims to find out how many times it takes to reach a steady state, mostly covering from 1 to 10. 8.
Fill the holes and achieve the bonding region connectivity by several erosion and dilation. In the field of image processing, connectivity means the continuous and adjacent of pixels. There are several cases that will happen when a hole is filled. In the first case, the size of the holes is relatively uniform, and can be fully filled though multiple erosion and dilation. In the second case, the size of the holes is quite different. After the small holes are filled, the large holes are incompletely filled. Furthermore, the large and unfilled holes can be seen that either the outer boundary is closed, meaning the pixels of hole contour are continuous and adjacent, or some values are only left on the boundary, meaning the pixels of hole contour are fragmentary. The number of processing times is as less as possible and the range is primarily between 1 and 10. 9.
Calculate the bonding rate. For the first case above, the parameter is directly gained for calculation. While, the second case demands to invoke the circle fitting function to measure the area of all the holes, including the unsuccessful filled holes. The area refers to the sum of the pixel points. Afterwards, the parameter is secured to calculate the bonding rate.

Orthogonal Experiment Design
The orthogonal experiment is designed in order to determine the optimal activated process parameters for high bonding strength and good interface performance. It is worthwhile to pay attention to the fact that the changes are only in activation parameters. All other parameters and experiment steps are unchanged. For O 2 plasma-activated Si-Si wafer direct bonding, there factors: Activated time, activated power, and O 2 flow are crucial for bonding quality. In addition, a verification factor is needed to examine the accuracy of the final result. The validation factor is a dummy horizontal value and does not refer to a specific parameter. Each factor is set at three levels. The activated time is: 30 s (level 1), 120 s (level 2), and 300 s (level 3). The activated power is: 50 W (level 1), 100 W (level 2), and 150 W (level 3). The oxygen flow is: 50 sccm (level 1), 100 sccm (level 2), and 150 sccm (level 3). The verification factor is: Level 1, level 2, and level 3. Therefore, the four-factor and three-level tests are employed, and the L9 (3 4 ) orthogonal table is used. The designed orthogonal experiments are shown in Table 1.

Software Evaluation
The software is used to inspect the bonding interface and calculate the bonding rate. To understand the impact of the setting parameters of each step on the image processing results, nine infrared transmission annealed images corresponding to nine sets of orthogonal experiments are evaluated, and the original images obtained by CCD camera are shown in Figure 3. Table 2 gives the different adjusting parameters of the software, which are used for further study of the characteristics of the image processing results. To be clear, changing the software parameters only makes the image processing more accurate, but has no impact on the actual bonding results. The accurate evaluation of bonding results requires comparison and verification with other detection methods.  Five parameters have a significant influence on image processing quality. Sharpening operator patterns determine the integrity of the high frequency components of the image and they are employed for image sharpening. The disk size resolves the strength of the image contrast. The larger the disk, the stronger the contrast, and the image details will be clearer. It is required to obtain the high-contrast filtered images with obvious details and no background interference. The binarization threshold can be reasonably set to accurately segment the bonding and non-bonding regions. The purpose of setting the morphological parameter A is to extract the main body, obtain the unbonded effective pixel points, and filter the noise pixel points introduced by binarization. The morphological parameter B is to realize hole filling through repeating corrosion and expansion. Premise of no loss of effective information, the times of corrosion and expansion is requested to be as small as possible, and the connected domain of bonding region is 1, that is, the bonding region is connected and closed without any imperfection.
In order to investigate the effect of these parameters on image quality and bonding rate, the bonding rates are combined with these parameters to optimize the software parameter conditions, the experiments are divided into five groups, as shown below.
No. (1): The four-neighbor is compared with the eight-neighbor sharpening operator pattern. If the different patterns are adopted, the corresponding four parameter values, the disk size, the binarization threshold, the morphological parameter A, and the morphological parameter B, ask to be matched to achieve the optimal effect.
Based on the selection of four neighborhood sharpening operators as pattern, the following research is carried out.
No. (2): The disk size is a variable ranging from 1 to 100, and the binarization threshold, morphological parameter A, and morphological parameter B are unchanged.
No. (3): The binarization threshold is a variable ranging from 0.2 to 0.5, and the disk size, morphological parameter A, and morphological parameter B are unchanged.
No. (4): The morphological parameter A is a variable ranging from 1 to 10, and the disk size, binarization threshold, and morphological parameter B are unchanged.
No. (5): The morphological parameter B is a variable ranging from 1 to 10, and the disk size, binarization threshold, and morphological parameter A are unchanged.

Four-Neighbor and Eight-Neighbor Sharpening Operator Patterns
The analysis results using different sharpening operator patterns are shown in Figures 4 and 5 respectively. It can be seen that more salt-pepper noise is introduced by using the eight-neighbor pattern compared to four-neighbor pattern. Although most of noise can be filtered after the first Five parameters have a significant influence on image processing quality. Sharpening operator patterns determine the integrity of the high frequency components of the image and they are employed for image sharpening. The disk size resolves the strength of the image contrast. The larger the disk, the stronger the contrast, and the image details will be clearer. It is required to obtain the high-contrast filtered images with obvious details and no background interference. The binarization threshold can be reasonably set to accurately segment the bonding and non-bonding regions. The purpose of setting the morphological parameter A is to extract the main body, obtain the unbonded effective pixel points, and filter the noise pixel points introduced by binarization. The morphological parameter B is to realize hole filling through repeating corrosion and expansion. Premise of no loss of effective information, the times of corrosion and expansion is requested to be as small as possible, and the connected domain of bonding region is 1, that is, the bonding region is connected and closed without any imperfection.
In order to investigate the effect of these parameters on image quality and bonding rate, the bonding rates are combined with these parameters to optimize the software parameter conditions, the experiments are divided into five groups, as shown below.
No. (1): The four-neighbor is compared with the eight-neighbor sharpening operator pattern. If the different patterns are adopted, the corresponding four parameter values, the disk size, the binarization threshold, the morphological parameter A, and the morphological parameter B, ask to be matched to achieve the optimal effect.
Based on the selection of four neighborhood sharpening operators as pattern, the following research is carried out.
No. (2): The disk size is a variable ranging from 1 to 100, and the binarization threshold, morphological parameter A, and morphological parameter B are unchanged.
No. (3): The binarization threshold is a variable ranging from 0.2 to 0.5, and the disk size, morphological parameter A, and morphological parameter B are unchanged.
No. (4): The morphological parameter A is a variable ranging from 1 to 10, and the disk size, binarization threshold, and morphological parameter B are unchanged.
No. (5): The morphological parameter B is a variable ranging from 1 to 10, and the disk size, binarization threshold, and morphological parameter A are unchanged.

Four-Neighbor and Eight-Neighbor Sharpening Operator Patterns
The analysis results using different sharpening operator patterns are shown in Figures 4 and 5 respectively. It can be seen that more salt-pepper noise is introduced by using the eight-neighbor pattern compared to four-neighbor pattern. Although most of noise can be filtered after the first morphological processing (step 7), the retention effect on the original information is worse than that of the four-neighbor pattern.
Micromachines 2019, 10, x FOR PEER REVIEW 8 of 15 morphological processing (step 7), the retention effect on the original information is worse than that of the four-neighbor pattern.   Figure 6 compares the bonding rates of the nine different images using four-neighbor and eightneighbor sharpening operator patterns. It can be clarified that, the bonding rates show good consistency for both methods, which indicates that each method can be used for the numerical calculation of the bonding rate. However, the bonding rates for some images are slightly different, and the difference is between 1% and 3%. Because the four-neighbor sharpening pattern has lower image noise, it is thought to be more accurate, and then applied for image enhancement in subsequent experiments. morphological processing (step 7), the retention effect on the original information is worse than that of the four-neighbor pattern.   Figure 6 compares the bonding rates of the nine different images using four-neighbor and eightneighbor sharpening operator patterns. It can be clarified that, the bonding rates show good consistency for both methods, which indicates that each method can be used for the numerical calculation of the bonding rate. However, the bonding rates for some images are slightly different, and the difference is between 1% and 3%. Because the four-neighbor sharpening pattern has lower image noise, it is thought to be more accurate, and then applied for image enhancement in subsequent experiments.  Figure 6 compares the bonding rates of the nine different images using four-neighbor and eight-neighbor sharpening operator patterns. It can be clarified that, the bonding rates show good consistency for both methods, which indicates that each method can be used for the numerical calculation of the bonding rate. However, the bonding rates for some images are slightly different, and the difference is between 1% and 3%. Because the four-neighbor sharpening pattern has lower image noise, it is thought to be more accurate, and then applied for image enhancement in subsequent experiments.

Disk Size and Binarization Threshold
The relationship between the bonding rates and the disk size is shown in Figure 7a. It can be seen that, the disk size has little effect on bonding rate. As the disk size increases from 10 to 100, the bonding rates only fluctuate within 2%. Especially when the disk size is set to be 20-50, the bonding rates are almost steady. However, a larger disk size may lead to the background interference. It was better to set it to 20 in the experiment. Figure 7b illustrates the relationship of bonding rate and binarization threshold. As the threshold rises up, the bonding rate tends to fall down. When the binarization threshold value is greater than 0.35, the bonding rate drops sharply. For the stable bonding rates, the binarization threshold is preferably 0.25-0.35, and the variation range is between 2% and 5%. Provided that don't import too many noise pixels, the complete information of the image will be kept as far as possible in the experiments. After repeated image tests of multiple groups, the binarization threshold was put to 0.3374.

Morphological Parameter A and B
The effect of the morphological parameter A on the bonding rate is shown in Figure 7c. When the parameter A changed, the bonding rates vibrate little, which are within 1%. Whereas, the bonding rate suddenly increases from 91.2% to 92.8% for sample No. 3 when the morphological parameter A changed from 7 to 8. It may be considered that the unbonded part is added to the connected domain due to the special morphology of the individual holes. As a consequence, the morphological parameter A universally takes between 1 and 7, and it was set to 1 in the experiments.
The effect of the morphological parameter B on the bonding rate is shown in Figure 7d. As the morphological parameter B grows up, the bonding rates tend to decrease slowly, which is because the voids or unbonded areas are excessively corroded. The sample No. 5 is most affected by the parameter B. The bonding rate decreases from 90.2% to 85.3% as the morphological parameter B varies from 1 to 10. When the morphological parameter B changes from 1 to 2, there is also a sudden change for the sample No. 1. Comparing the sample No. 1 and No. 5, it is found that the obvious change of bonding rate is related to the small holes distributing in the infrared image. After the repetitive corrosion and expansion treatment, the small holes will continue to expand. When the number of holes is large enough, the small holes join together, resulting in a huge change of bonding rate. In order to obtain an accurate bonding rate for all the samples, it is recommended that the morphological parameter B should be set to 4 and 6, and the value in the experiment was 4.

Disk Size and Binarization Threshold
The relationship between the bonding rates and the disk size is shown in Figure 7a. It can be seen that, the disk size has little effect on bonding rate. As the disk size increases from 10 to 100, the bonding rates only fluctuate within 2%. Especially when the disk size is set to be 20-50, the bonding rates are almost steady. However, a larger disk size may lead to the background interference. It was better to set it to 20 in the experiment.

Void Formation and Characteristics
The infrared image processing software is used for studying the void formation and characteristics. Five representative infrared images were selected, and the process parameters are shown in Table 3. The histogram statistics were performed on the number and diameter of voids, as shown in the Figure 8. It illustrates the size and distribution of voids at different process parameters.  Figure 7b illustrates the relationship of bonding rate and binarization threshold. As the threshold rises up, the bonding rate tends to fall down. When the binarization threshold value is greater than 0.35, the bonding rate drops sharply. For the stable bonding rates, the binarization threshold is preferably 0.25-0.35, and the variation range is between 2% and 5%. Provided that don't import too many noise pixels, the complete information of the image will be kept as far as possible in the experiments. After repeated image tests of multiple groups, the binarization threshold was put to 0.3374.

Morphological Parameter A and B
The effect of the morphological parameter A on the bonding rate is shown in Figure 7c. When the parameter A changed, the bonding rates vibrate little, which are within 1%. Whereas, the bonding rate suddenly increases from 91.2% to 92.8% for sample No. 3 when the morphological parameter A changed from 7 to 8. It may be considered that the unbonded part is added to the connected domain due to the special morphology of the individual holes. As a consequence, the morphological parameter A universally takes between 1 and 7, and it was set to 1 in the experiments.
The effect of the morphological parameter B on the bonding rate is shown in Figure 7d. As the morphological parameter B grows up, the bonding rates tend to decrease slowly, which is because the voids or unbonded areas are excessively corroded. The sample No. 5 is most affected by the parameter B. The bonding rate decreases from 90.2% to 85.3% as the morphological parameter B varies from 1 to 10. When the morphological parameter B changes from 1 to 2, there is also a sudden change for the sample No. 1. Comparing the sample No. 1 and No. 5, it is found that the obvious change of bonding rate is related to the small holes distributing in the infrared image. After the repetitive corrosion and expansion treatment, the small holes will continue to expand. When the number of holes is large enough, the small holes join together, resulting in a huge change of bonding rate. In order to obtain an accurate bonding rate for all the samples, it is recommended that the morphological parameter B should be set to 4 and 6, and the value in the experiment was 4.

Void Formation and Characteristics
The infrared image processing software is used for studying the void formation and characteristics. Five representative infrared images were selected, and the process parameters are shown in Table 3. The histogram statistics were performed on the number and diameter of voids, as shown in the Figure 8. It illustrates the size and distribution of voids at different process parameters. For the Sample No. 1 , there were many interference fringes distributed in the infrared image. The void diameters were mainly in the range of 14 to 22 mm. The interference fringes were multilayer and irregular circular, which were caused by organic pollutants. As the annealing temperature increased to 350 • C, the organic contaminant decomposed into gases which were trapped in the bonded interface and the large voids were formed. In addition, half of the void diameters of 8~14 mm were mostly attracted by small particles. It was obvious that the interference fringes are clear. All unbonded areas were successfully divided and the connected domain was captured. However, the bonding rate was only 69.12%, so the bonding wafers must be further cleaned by chemical cleaning or plasma activation process.
For the sample No. 2 , a large number of voids (total number >100) was detected. The voids were mainly divided into three groups: Smaller than 2 mm, 2-4 mm, and larger than 8 mm. The smaller voids were mainly caused by H 2 molecules produced during the annealing process. The voids with the diameter of 2-4 mm may have been the result of the aggregation of water molecules. Moreover, large regular circular voids with the diameter of 8~14 mm were induced by particle contamination. The large number of voids with the diameter of 2-4 mm means that, there was surplus water between the bonding interface, and the bonding surface was over treated. With a suitable binarization threshold which can keep all the details of the image, the correct segmentation can be achieved from the resultant image. The calculated bonding rate was 86.08%. For the sample No. ②, a large number of voids (total number >100) was detected. The voids were mainly divided into three groups: Smaller than 2 mm, 2-4 mm, and larger than 8 mm. The smaller voids were mainly caused by H2 molecules produced during the annealing process. The voids with the diameter of 2-4 mm may have been the result of the aggregation of water molecules. Moreover, large regular circular voids with the diameter of 8~14 mm were induced by particle contamination. The large number of voids with the diameter of 2-4 mm means that, there was surplus water between the bonding interface, and the bonding surface was over treated. With a suitable binarization threshold which can keep all the details of the image, the correct segmentation can be achieved from the resultant image. The calculated bonding rate was 86.08%.
The interface of the sample No. ③ was evenly occupied by a great quantity of voids with a diameter smaller than 2 mm, which accounted for more than 95%. These small voids were caused by the byproduct H2 formed by the condensation reaction between Si-OH. Meanwhile, a small amount of bigger voids (2-4 mm) were also observed. It indicated that there was still a slight excess of water molecules at the interface, and the surplus reaction also produced excessive H2, which could not be absorbed by the Si substrate. By re-filtering and enhancement, the voids can be successfully captured to undertake more precise results. Although the tiny voids (mainly <1 mm in diameter) could not be successfully extracted, it had only a one in ten thousand (voids to wafer area ratio) impact on the final bonding result.
For the sample No. ④, the infrared image had fewer voids. The small H2 molecules were absorbed by the amorphous oxide layer between the bonded interfaces, so the small voids disappeared. The voids with the diameter ranging from 2 to 8 mm were mostly produced by small particles which could be eliminate by proper cleaning. There were also multiple unbonded areas caused by tweezers at the edge of the bonded wafer. Via sharpening, they can be easily obtained and completely removed to ensure the accuracy of result. Worthy of mention is that the bonding rate approximated that of No. ③, which was ~94%. In the light of Table 3, it can be seen that the number of voids is much less than sample No. ③. It is beneficial to have sufficient contact of the bonded wafers so that higher bonding strength can be achieved. The interface of the sample No. 3 was evenly occupied by a great quantity of voids with a diameter smaller than 2 mm, which accounted for more than 95%. These small voids were caused by the byproduct H 2 formed by the condensation reaction between Si-OH. Meanwhile, a small amount of bigger voids (2-4 mm) were also observed. It indicated that there was still a slight excess of water molecules at the interface, and the surplus reaction also produced excessive H 2 , which could not be absorbed by the Si substrate. By re-filtering and enhancement, the voids can be successfully captured to undertake more precise results. Although the tiny voids (mainly <1 mm in diameter) could not be successfully extracted, it had only a one in ten thousand (voids to wafer area ratio) impact on the final bonding result.
For the sample No. 4 , the infrared image had fewer voids. The small H 2 molecules were absorbed by the amorphous oxide layer between the bonded interfaces, so the small voids disappeared. The voids with the diameter ranging from 2 to 8 mm were mostly produced by small particles which could be eliminate by proper cleaning. There were also multiple unbonded areas caused by tweezers at the edge of the bonded wafer. Via sharpening, they can be easily obtained and completely removed to ensure the accuracy of result. Worthy of mention is that the bonding rate approximated that of No. 3 , which was~94%. In the light of Table 3, it can be seen that the number of voids is much less than sample No. 3 . It is beneficial to have sufficient contact of the bonded wafers so that higher bonding strength can be achieved.
There was nearly no defect in the sample No. 5 . The only two tiny voids (diameter~2 mm) may have been caused by incompletely diffused H 2 molecules. The boundary of silicon wafer was intact and distinct, the bonding rate approached 100%.
To sum up, the software can accurately analyze the various voids with different shapes and sizes. The cause of void formation and the bonding quality can be identified through the number and diameter of voids.

Orthogonal Experiment Analysis
Nine groups of bonding experiments with different activation conditions were performed based on Table 1. Each group of experiments was implemented, and the infrared transmission image was collected as well as the bonding rate calculated for pre-bonding and annealing. The orthogonal test table was applied to solve the investigation index, which not only reduces the number of trials, but also achieves the best results.
Through computing the mean and extremum value, the impact of different factors can be resolved, as shown in Table 4  According to the orthogonal experiment, the O 2 plasma activated process is optimized with the power of 150 W, O 2 flow of 100 sccm, and time of 60 s. The infrared image is shown in Figure 9a.
There are a few voids in the bonded interface after annealing. It shows relatively higher bonding rate (~94.51%). Moreover, the bonding strength is performed by the universal tensile testing machine SUNS UTM4104X. The bonded wafer is diced into small pieces with the size of 5 × 5 mm 2 , and then glued on the testing bar. The pulling force is applied on both sides of the testing bar until the samples are cracked. The result shows that the fracture force is 308 N, and the fracture occurs at the silicon bulk, as shown in Figure 9b. This indicates that the strong bonding strength is 12.32 MPa.

Conclusions
In this study, we developed a quality detection method of bonded interface for O2 plasma assisted silicon wafer bonding based on spatial domain and morphology. This method needed simpler decomposition and reorganization process compared to wavelet denoising. The normalized high-pass filtering avoided image brightness offset, protected sharpening effect, and restored background details. Besides, the infrared imaging system and the image analysis software were constructed for inspecting the bonded interface. The calculation of wafer bonding rate was realized by image analysis software. The influence of software parameters including sharpening operator patterns, disk size, binarization threshold, morphological parameter A and B on bonding rate was studied as well. At the same time, the voids and their distributions of the bonded wafer interface before and after annealing were investigated. The various void formation and characteristics were accurately analyzed by the software. Ultimately, the orthogonal experiments were designed to determine the optimal activated process parameters for high bonding strength and good interface performance. Based on the comparison of different methods and parameters what are employed in the software evaluation and void characteristics, the following conclusions are finally drawn: (1) Comparing the effects of four-neighbor and eight-neighbor sharpening operator pattern, it is well known that more salt and pepper noise will be imported using the eight-neighbor patter. After the first morphological processing, although most of the noise could be filtered out, the retention effect on the original information was worse than that of the four-neighbor pattern. As a result, the four-neighbor pattern was used for image enhancement processing in subsequent experiments. (2) Among the four parameters of disk size, binarization threshold, morphological parameter A and morphological parameter B, the binarization threshold had the most significant effect on the result, and the variation range of result was primarily maintained at 2%~5%. The other parameters all had an effect on the results within 2%. In addition, while processing images, it was also necessary to pay attention to the holes of individual features to bring about the effect of jump. (3) The number of voids in the bonding interface can identify the bonding quality, while the diameter can infer the cause of void formation. For the unclear outline of voids (with <1 mm in diameters) caused by the byproduct H2 formed by the condensation reaction between Si-OH, they cannot be captured during software processing. Whereas, the final bonding result has only a one in ten thousand effect. The voids with a diameter ranging from 2 to 14 mm caused by the aggregation of water molecules and small particles can be easily acquired and completely removed via software processing, ensuring the accuracy of the results. Caused by large dust or organic contaminants, the voids were successfully divided and the connected domain was gained between 14 and 22 mm in diameter, which had a serious impact on the bonding rate. (4) The orthogonal experimental data show that the O2 flow has the greatest influence on bonding rate in comparison with activated time and power. The subsequent experiments particularly focused on the parameter selection of O2 flow. With the optimized process parameters of

Conclusions
In this study, we developed a quality detection method of bonded interface for O 2 plasma assisted silicon wafer bonding based on spatial domain and morphology. This method needed simpler decomposition and reorganization process compared to wavelet denoising. The normalized high-pass filtering avoided image brightness offset, protected sharpening effect, and restored background details. Besides, the infrared imaging system and the image analysis software were constructed for inspecting the bonded interface. The calculation of wafer bonding rate was realized by image analysis software. The influence of software parameters including sharpening operator patterns, disk size, binarization threshold, morphological parameter A and B on bonding rate was studied as well. At the same time, the voids and their distributions of the bonded wafer interface before and after annealing were investigated. The various void formation and characteristics were accurately analyzed by the software. Ultimately, the orthogonal experiments were designed to determine the optimal activated process parameters for high bonding strength and good interface performance. Based on the comparison of different methods and parameters what are employed in the software evaluation and void characteristics, the following conclusions are finally drawn: (1) Comparing the effects of four-neighbor and eight-neighbor sharpening operator pattern, it is well known that more salt and pepper noise will be imported using the eight-neighbor patter. After the first morphological processing, although most of the noise could be filtered out, the retention effect on the original information was worse than that of the four-neighbor pattern. As a result, the four-neighbor pattern was used for image enhancement processing in subsequent experiments. (2) Among the four parameters of disk size, binarization threshold, morphological parameter A and morphological parameter B, the binarization threshold had the most significant effect on the result, and the variation range of result was primarily maintained at 2%~5%. The other parameters all had an effect on the results within 2%. In addition, while processing images, it was also necessary to pay attention to the holes of individual features to bring about the effect of jump. (3) The number of voids in the bonding interface can identify the bonding quality, while the diameter can infer the cause of void formation. For the unclear outline of voids (with <1 mm in diameters) caused by the byproduct H 2 formed by the condensation reaction between Si-OH, they cannot be captured during software processing. Whereas, the final bonding result has only a one in ten thousand effect. The voids with a diameter ranging from 2 to 14 mm caused by the aggregation of water molecules and small particles can be easily acquired and completely removed via software processing, ensuring the accuracy of the results. Caused by large dust or organic contaminants, the voids were successfully divided and the connected domain was gained between 14 and 22 mm in diameter, which had a serious impact on the bonding rate.