Effects of static load and residual stress on fused silica direct bonding interface properties

Defect free direct bonding of rigid and large area glass samples, such as prisms, becomes increasingly important for the manufacturing of modern optical and optomechanical components. Typically, in order to apply a static load during the annealing step, specialized heat-resistant pressure mountings are required. This makes manufacturing effortful and cost-intensive. In this paper, we present plasma activated bonding experiments conducted on fused silica plates where residual stress has been introduced prior to the contacting step and where annealing is performed with and without a static load. We find that in case of a sufficiently smooth surface, bonding strength is insensitive towards residual stress or static load, or more precisely, towards the interface stress. Furthermore, the residual Fresnel reflection losses of the realized bonding interface were optically measured and they amount to only 10-6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{-6}$$\end{document}. We propose that a consideration of the change in Gibbs free energy, dG, allows qualitatively predicting the resulting bonding strength and its spatial distribution, where dG is determined by surface energy and interface stress. At the end of this article, conceivable applications are discussed.


Introduction
Direct bonding is a bonding technique which allows permanent and high strength joining of a variety of similar and dissimilar materials (including optical glasses, silicon, germanium, metals and laser crystals) without introducing any intermediate layer [1][2][3][4][5]. Adhesion is provided solely by intermolecular forces between the surfaces which, however, have to be planar, smooth and clean [6,7]. It has been developed in the 1980s for electronic applications based on silicon wafers [8], which is still today the most important field of application [1,6,9]. Generally, direct bonding employs preceding hydrophilic surface activation steps and a successive annealing step during which a static load is applied onto the workpieces to be bonded [10][11][12]. In case of wafer direct bonding, the static load during the annealing step is typically provided by two "chucks", i. e., well planarized pressure plates made of highly stiff materials, e. g., technical ceramics [3], because it is proposed that the application of a static load enhances bonding strength and reduces defects by overcoming the wafers' local waviness [10].
Considering the surface physics mechanisms, it has been widely accepted that after the contacting and prior to the annealing step adhesion is mainly provided by surface molecules which exhibit van-der-Waals bonding and hydrogen bonding. Then, during the annealing step, these initially weakly bonded surface molecules react to become covalently bonded [6,8,13,14]. This temperature activated process strongly increases the adhesion [7,[15][16][17]. The underlying chemical reaction is the so-called silanol condensation, sometimes also referred to as siloxane condensation, with equation which summarizes the hydrophilic mechanism for silicon or silica direct bonding [15,16,18,19]. Considering the removal of interface H 2 O, a two phase mechanism has been (1) Si-OH + HO-Si → SiO 2 + H 2 O proposed [20]. In the first phase, the contacted surface area enlarges and the interface width increases, where most of the water stays trapped at the not yet contacted zones. In the second phase, water diffuses out of the interface and reacts with the SiO 2 network, cf. Ref. [21], leading to a reversible formation of immobile Si-OH groups, cf. Ref. [22], with a resulting interface width of 3 to 4 nm [20].
Due to more temperature sensitive parts in the fields of microelectronics, three-dimensional device integration, micro-electro-mechanical systems [23,24], micro-fluidics [25] and optical, optomechanical and photonics applications [19,26,27], increased interest has emerged on direct bonding techniques at lower annealing temperatures in order to avoid thermal stress but still obtain permanent and high bonding strength. This is typically made possible by employing a plasma activation step prior to the contacting step, allowing annealing temperatures about or below 250 • C [28][29][30]. The corresponding bonding technique is then referred to as plasma activated bonding (PAB). In addition, many authors employ a de-ionized water rinse step immediately after the plasma-activation step and prior to the contacting step [17,18,[31][32][33][34][35][36] which is proposed to increase the hydrogen bonding interaction distance and thus to improve the spontaneity of adhesion during the contacting step [37,38].
More recently, PAB of optical materials has especially become interesting for the manufacturing of large-scale optical and optomechanical components for space-based, lithographic and high power laser applications [4,19,[39][40][41][42], driven by properties like homogeneously high bonding strength, superior transmittance, chemical resistance and long-term stability [3,17,19,34,39,43]. Finally, the ability to directly bond large optical surfaces is mandatory for specific applications where throughput and imaging resolution scales with the achievable size of the optical aperture, cf. Ref. [44].
In order to obtain a homogeneous bonding interface, it is important that a careful preceding surface planarization is performed. This is in particular true when the workpieces to be bonded are stiff and rigid (e. g., consider a large, rigid prism in contrast to a thin, elastic wafer). In addition, as a state of the art, during the annealing step a static load of up to several MPa is employed [9,40,45]. However, considering the static load, from a theoretical point of view there is no reason to do so: Once a full-area contact between two planarized surfaces is established, the initial "weak" adhesion should maintain this contact throughout the annealing step. Taking this hypothesis further, one could suggest that even initially slightly convex surfaces which are brought into full-area contact and maintain this contact, so that the workpieces are now mechanically stressed, should be bondable without any static load during the annealing step. This suggests the investigation of the effects of a static load during the annealing step and of residual stress on direct bonding interface properties.
As to knowledge of the authors, there is no published research on glass direct bonding concerning the influence of static load applied during the annealing step on direct bonding success. Note that there is some literature on metal direct bonding concerning that topic [9,[46][47][48], but it is not transferable to surface physics mechanisms of glass direct bonding because of very different material ductilities [6]. For glass direct bonding, this research topic is interesting for two reasons: Firstly, from a practical point of view, it is worthwhile noting that specialized (i. e., cost-intensive) heat-resistant pressure mountings are typically being employed today when manufacturing large-scale optical and optomechanical components via direct bonding. These mountings would become no longer necessary which implies that manufacturability of such components would be greatly improved. Secondly, from an academic point of view, when accessing new aspects of direct bonding the results might shed light on undiscovered underlying surface physics mechanisms.
Internal preliminary PAB experiments (performed on stiff fused silica round plates with 1" diameter) suggest that the employment of a static load during the annealing step does not affect bonding strength: We measured values of 56.7 ± 12.5MPa and 55.7 ± 10.3MPa (given as average bonding strength value ± 99 % confidence interval) in case of a static load and no static load, where sample size was 33 and 26, respectively. However, this data is not sufficiently clear, because samples where from different manufacturing charges. Also, we did not gain information about spatial distribution of bonding strength. Furthermore, a small bonding area is not fully representative to a large one, because, on the one hand, there may be edge effects [49,50] and on the other hand, sample residual stress may be different. In this paper, our motivation is the investigation of direct bonding of large area and stiff fused silica samples, because of its important role for state-of-the-art large-aperture-optics applications, as described above. The key properties for such applications are bonding strength, interface reflectivity, as well as its spatial distribution. Therefore, for our findings to be relevant for the application, we examined these properties via performing PAB experiments on large area and stiff fused silica samples with focus on the effects of static load and residual stress.

Experimental
Samples are 6.35 mm thick rectangular plates with 110 × 66 mm 2 bonding surface area. They have been sawn out of uncoated fused silica blanks with a planarity verified to be smaller than 127 nm as peak-to-valley (PV) value. Thus, the samples can be considered as stiff and highly planar. A smooth surface is a key factor for successful direct bonding, too [51]. This is why for all samples but one, chemical-mechanical polishing (CMP) is performed. From surface profiles that have been measured via atomic-force-microscopy (AFM) on both a 1 × 1 μm 2 and a 10 × 10 μm 2 measurement field, the root-mean-square (RMS) surface roughness, S q , is calculated.
Although a highly planar surface is also a key factor for successful direct bonding [52], we proceed in the following way. Via lapping, we prepared upper the vacuum chuck (which is holding the upper of the bonding partners during the contacting step) in such a way that its surface is convex with a PV value of 5 μm . The bottom vacuum chuck is concave with a PV value of 3 μm . This causes the corresponding bonding surface to become convex when chucked, with a PV value of approximately 4 μm , see Fig. 1. As the top chuck has a higher PV value, this method is proposed to enhance contact front propagation control and thus mitigate the risk of creating interface voids. Also, it introduces radial symmetric residual stresses. Planarity measurements have been conducted via Fourier transform phase shifting interferometry (using the device "Verifire MST" from ZYGO).
Direct bonding is performed using typical PAB process steps, i. e., surfaces to be bonded are treated with a mixture of ammonia water and hydrogen peroxide, highly diluted in de-ionized water (DIW), based on Ref. [53] -often referred to as "RCA1"-cleaning -, low-pressure plasma activation with the process gases O 2 and N 2 , consecutively (at a substrate-anode distance of about 100 mm in a static setup, using a commercial plasma activation tool by EV Group, Austria), a mega-sonic assisted DIW rinse, and dried on a spinning stage for 120 s before brought into contact. Contacting is performed in atmospheric clean room conditions (class ISO 2) to avoid particle contamination. To obtain a full area contact, which is verified by completely vanishing Newton rings on the full area, contact front propagation is monitored, and controlled by applying a temporary, locally adjustable static load, cf. Ref. [40]. Finally, in a vacuum furnace with pressure below 0.1 mbar, each sample pair is annealed at 200 • C for 20 h. During the annealing step, some sample pairs are subject to a homogeneous static load, P, of several MPa, a value that is typical for wafer direct bonding, and some are subject to no static load. An overview of the experiments is given in Table 1.
After annealing, sample pairs are characterized with respect to their interface's residual reflectivity (also known as Fresnel reflectivity) using angle resolved optical light measurement at wavelength 633 nm and incident angle 25 • . Finally, sample pairs are characterized with respect to their bonding strength. Typically, in wafer direct bonding, the blade method (also referred to as dual-cantilever bending (DCB) method, cf. Ref. [54]) is used. In our case, however, the blade method cannot be applied because of the high thickness and corresponding rigidity of the bonded samples, cf. Ref. [19]. We therefore performed a three-point flexural test. It is a common and standardized method in fracture strength quantification of brittle materials [55,56]. To this end, sample pairs are cut by diamond sawing into numerous pieces à 5 × 5 × 12.7 mm 3 . These pieces are then put into setup for destructively obtaining the flexural strength, f , which we refer to as bonding strength. Here, support span width was 11 mm. Prior to sawing, pieces are labeled using an automatic laser writing tool in order to later allocate their original position. Thus, we are able to investigate the spatial distribution and homogeneity of the interface's bonding strength. The corresponding process steps are shown in Fig. 2. In addition, flexural strength measurements have also been performed on bulk fused silica samples in order to compare the bonding strength with the bulk material str ength.

Results
From the surface profiles, measured via AFM prior to direct bonding, we have calculated the RMS surface roughness values. We find S q = 0.5 nm on the original surfaces. On the CMP treated surfaces we find S q = 0.2 nm which can be considered as highly smooth. Exemplary surface profiles are shown in Fig. 4 which are representative for the other samples. In case of all five direct bonding experiments but one, the controlling of the contact front propagation went very well, i. e., we were able to establish full-area and completely defect free contact. A photograph of an exemplary, representative successful directly bonded sample pair is shown in Fig. 5. In one case, namely sample pair #1, unfortunately we captured a small particle between the bonding surfaces. That particle evoked a large defect at the interface's edge; its position and extent is schematically shown in Fig. 6, a photograph is shown in Fig. 7.
After performing the direct bonding experiments, all sample pairs have been measured with regards to their bonding interface reflectivity. The results are shown in Table 2. It becomes evident that the application of a static load, P, during the annealing step does not affect the reflectivity. However, an increase in surface roughness does, see sample pair #5. In case of sample pair #5, the reflectivity has been measured at two positions, namely at the center and approximately 10 mm next to one of edges. The underlying       sample pairs exhibit very similar average values while the last one, i. e., sample pair #5 where one bonding partner has not been CMP treated, has a smaller average value. However, there is nearly an overlap between all confidence intervals so that no fully convincing statement can be made about a significant difference. Therefore, in order to access the data in a more meaningful way, a more sophisticated statistical analysis is necessary. For this, let us consider the survival distribution, S( f ) , i. e., the function which gives the probability that a piece breaks beyond a specific stress value. Figure 10 depicts the measured survival distributions. It becomes evident that the first four sample pairs' curves lie relatively closely together, but that of sample pair #5 is clearly shifted towards smaller bonding strength values. In order to model the measured data, let us assume that the survival distribution follows a Weibull distribution which is the most established and convenient method for reliability analysis of brittle materials [57]. According to that model, the survival function is given by where the fitting parameters, f ,c and m, are referred to as critical strength and Weibull modulus, respectively.
The following statistical analysis is performed by using the software package reliability by reid [58]. Firstly, for each sample pair's data set, we fit the parameters from Eq. 2. Here, we find that the resulting 99 % confidence intervals largely overlap for sample pairs #1 to #4 but do not overlap for sample pair #5 (not depicted). This is why we may hypothesize that sample pairs #1 to #4 stem from the same underlying distribution, i. e., do not differ in average bonding strength. This hypothesis is referred to as H 0 .
As such, we may combine their data and calculate the fitting parameters of their joint survival distribution. We find f ,c = 39 MPa ± 1 MPa and m = 4.8 ± 0.5 (with uncertainty given as standard error) and covariance 0.14MPa. In Fig. 10, the error bands corresponding to the regions of 1 , 2 , 3 and 4 are depicted in different shades of gray, where is the calculated estimated standard error of S( f , f ,c , m) , cf. Ref. [59]. In the following, we refer to that distribution as S * . Eventually, by performing a KolmoGoroW-smirnoW test, we find that we must reject the hypothesis that curve #5 stems from S * , where type I error is = 0.01 . In contrast, we find that we can not reject H 0 , i. e., the hypothesis that the other four curves stem from S * .
By further regarding Fig. 10, note that curves #1 and #2 look slightly shifted towards higher bonding strength values as compared to curves #3 and #4. This suggests that the application of a static load might increase bonding strength. Now, a question concerning statistical sensitivity arises: If we conclude that the application of static load does not change bonding strength, what is the probability that this conclusion is false, i. e., what it the type II error of accepting H 0 ? Vice versa, if we assume a type II error of = 0.05 , what is the minimal bonding strength difference, | f ,c ("with static load") − f ,c ("no static load")| , that would keep us from mistakenly accept H 0 ? As result, we find that this difference is 6.2MPa. Thus, we have demonstrated, firstly, that curve #5 does significantly differ from the other four curves, which implies that a larger surface roughness significantly decreases bonding strength. Secondly, we have demonstrated that if the application of a static load during the annealing step (in case of highly smooth surfaces) has a significant effect on bonding strength, this effect would be  Note that at the bottom right corner, no data could be retrieved because of a defect due to a particle, cf. Figs. 6 and 7 smaller than 6.2MPa. (Note that 6.2MPa is a relatively small value as compared to bonding strength variation from, e. g., temperature modification, using another optical glass material, cf. Ref. [17,39], or using a different surface preparation as shown in this article). Now, the spatial distribution of bonding strength is considered. It has been measured for sample pairs #1, #3 and #5. The results are depicted as heat maps in Figs. 11, 12 and 13, respectively. In case of sample pairs #1 and #3, we find an isotropic and homogeneous distribution of bonding strength values (disregarding the area of locally smaller values in the particle position's vicinity in case of sample pair #1). In case of sample pair #5, we find a rather inhomogeneous distribution of bonding strength values, where smaller values occur more frequently in the center area. (Or, more precisely, when fitting a centered paraboloid, f (x, y) = a 0 + a 2,x x 2 + a 2,y y 2 , to the data, the coefficients a 2,x and a 2,y significantly differ from zero only in case of sample pair #5, where again = 0.01 has been used.) This indicates that a CMP treatment results in a much more homogeneous bonding strength distribution. Note that eight data points from sample pair #1, namely those taken in the particle position's vicinity, have been neglected for the statistical analysis above, because they are considered to have emerged from a systematic rather than a statistical bias. The raw data supporting the findings described above is available from the corresponding author upon reasonable request.

Discussion
The results presented above indicate that, in case of highly smooth surfaces, the application of a static load during the annealing step neither affects bonding strength (measured within a reasonable statistical sensitivity) nor reflectivity, as shown for sample pairs #1 to #4. An increase in surface roughness however does affect bonding strength and interface reflectivity, as shown for sample pair #5, where bonding strength is reduced while its spatial distribution is less homogeneous, and interface reflectivity is increased while its spatial distribution is still homogeneous.
Considering sample pair #1, where an interface defect due to an enclosed particle occurred, it is worthwhile noting that the particle position (see Figs. 6 and 7) clearly correlates with the bonding strength distribution (see Fig. 11). Thus, it is a fully deterministic effect rather than a statistical or sawing preparation artifact. In particular, this result clearly shows that the bonding strength distribution is spatially homogeneous except for the area around the particle.
Considering bonding strength measurements in general, our results demonstrate that a reliability analysis under the assumption of a Weibull distribution is a powerful tool for comparing different strength data sets. This was expected because the Weibull distribution is known to be a suitable model for describing reliability of brittle materials in general [57,60], especially if the number of strength measurements is large [61,62] as it is the case in our experiments.
Typical bonding strength average values (received from earlier fused silica direct bonding experiments) are ranging from one half to two thirds of the bulk material strength [3,63]. However, it is worthwhile noting that those experiments have been conducted on smaller bonding areas and also, which is proposed to be more decisive, at an annealing temperature of 250 • C rather than 200 • C. The positive correlation between bonding strength and annealing temperature for PAB in that temperature range is well documented [17,34]. Still, as our results show bonding strength values of about one third of the bulk material strength. For intended optics applications where a large bonding area is needed, this value can be considered as sufficiently high, cf. Ref. [39]. Note that the local occurrence of much higher bonding  Note that in this case, smaller values tend to be found in the center (except for two outliers on the right edge). The colors are subject to the same colorbar as used in Fig. 11 strength values (e. g., 65 MPa in Fig. 12 top right) is a consequence of statistical scatter which is typical for fracture measurements.
An important observation has been made in case of sample pairs #3, #4 and #5 which, as a reminder, have been annealed without the application of a static load during the annealing step. Here, prior to the annealing step, no spontaneous delamination has been observed although the convex vacuum chuck exposed the sample pairs to residual stresses. In the following, as originally proposed by proposed by Turner and spearinG [64,65] (based on earlier works by TonG, Gösele and Yu [52,66]) we assume that this observation implies that the initially "weak" adhesion, or more precisely, the surface energy, , which keeps the contacted surfaces adhered together is larger than the change of strain energy, d U E , which would have been released with change of contacted surface area, dA, due to delamination. Thus, is the condition for preventing spontaneous delamination. Note that the factor 2 is used to accommodate adhesion from both surfaces. The key to use Eq. 3 is to develop an expression for U E as a function of A [64].
For estimating the minimum required surface energy to prevent spontaneous delamination, min = dU E ∕dA ∕ 2 , we consider that the PV value of the chucked bonding surface is s = 4 μm and that the surface itself is spherically convex. As the PV value is much smaller than the sample's thickness, t, it is sufficient to only regard stresses resulting from pure bending [67]. Using classical plate theory, it can be shown that dU E ∕dA = E t 3 2 ∕ [12 (1 − )] , where E and are the material's YounG's modulus and poisson's ratio, respectively, and is the plate's surface curvature [65]. Eventually, by applying the small sagitta approximation, s = L 2 ∕ 8 , we obtain from which, when assuming L = 110 mm, t = 6.35 mm, E = 72 GPa and = 0.18 , we find min = 6.6 mJ m −2 .
For highly smooth surfaces, i. e., S q = 0.2nm, as it is the case for sample pairs #1 to #4, it is known from direct measurements [7,16,17,51,[68][69][70] that the initial surface energy approximately ranges from = 80 to 100mJ m −2 . Evidently, this is much larger than min which explains the observation that delamination did not occur. In the following, this relation is denoted as min ≪ . For a less smooth surface, i. e., S q = 0.5nm, as it is the case for sample pair #5, we can estimate the surface energy based on measurements by larreY et al. [51] to be approximately = 20mJ m −2 . This value is still larger than min (but in a similar order of magnitude) which again explains that delamination did not occur. In the following, this relation is denoted as min ≲ . Now let us again turn our attention to the observation that for sample pair #5 the spatial distribution of bonding strength is not homogeneous anymore but rather resembles a radial symmetric distribution which correlates to the residual stress distribution that has been introduced by the vacuum chuck, see Fig. 13. In contrast, for sample pair #3 the spatial distribution is completely homogeneous, see Fig. 12, although residual stress has been introduced, as well. Evidently, the underlying residual stress affects the bonding strength distribution in case of min ≲ but not in case of min ≪ . Considering Fig. 13 again, note that the reason that smaller bonding strength values occur more frequently in the center area is that in case of slightly convex samples min has its maximum value directly at the point of initial contact, i. e., in the center area, cf. Ref. [65].
Bonding strength is determined during the annealing step by near-interface properties, for example, by the elastic strain energy stored near the bonding interface [7,17] which can also be described as interface stress. Here, it is irrelevant if the interface stress is determined by internal or external conditions, i. e., by residual stresses or a static load, respectively. Note that residual stresses can both increase and decrease the interface stress [7], while a static load typically only decreases it. Here, the latter case is proposed to occur in such a way that nano-asperities are flattened, so that less energy has to be applied to overcome the intermolecular gaps between the bonding surfaces, thus reducing interface stress. As it has been suggested by plach et al. [17], on a microscopic scale formation of covalent bonds occurs site specific such that the Gibbs free energy, G, that is in this case the elastic strain energy minus the surface energy times newly contacted area, is minimized. It can be shown that the condition ≥ min from above is equivalent to the condition d G ≤ 0 , which is known as criterion of spontaneous change [71]. Furthermore, note that in case of d G ≤ 0 a decrease of dG increases the chemical reaction rate.
Taking all the observations from above into account, we conclude that in case of a not so smooth surface, bonding strength is affected by interface stresses, and in case of a highly smooth surface, bonding strength is not affected by interface stresses anymore. This could be regarded as a saturation effect where the reaction rates of the underlying chemical reactions, as they are summarized in Eq. 1, are in case of min ≪ determined solely by chemical gradients, i. e., differences in molecular concentration, rather than interface stress which, according to plach et al., is correlated with contacted area. We therefore propose the following categorization.
(i) If d G > 0 , spontaneous delamination will occur right after the contacting step. Defect free direct bonding is possible only if there is a permanent static load during the annealing step. This is the case if surface roughness is rather high. (ii) If d G ≲ 0 , spontaneous delamination will not occur after the contacting step. The resulting bonding strength will be interface stress sensitive. Therefore, it may vary as a function of residual stress or static load. (iii) If d G ≪ 0 , spontaneous delamination will not occur after the contacting step, either. However, the resulting bonding strength will be spatially homogeneous, independent of interface stresses. Therefore, a static load does not affect bonding strength. This is the case if surfaces are highly smooth.
A schematic overview of these three cases is shown in Fig. 14. It is worthwhile noting that this model is a twodimensional simplification, because it does not consider influences of, e. g., the plasma activation type or the annealing temperature. Still, this model allows understanding the emergence of two fundamental transitions, namely those between spontaneous delamination, interface stress sensitive bonding strength and interface stress insensitive bonding strength. While the first phase transition is known from literature [10][11][12]65], the second one has not been described before. The schematic overview, see Fig. 14, suggests that only if there is an extraordinarily high interface stress, the bonding strength may spatially vary. Indeed, this has been observed as well, namely, in case of sample pair #1 where a particle has been included between the bonding interfaces. Here, we indeed observe a localized decrease in bonding strength, see Fig. 11. A few mm further away from the particle inclusion, no change in bonding strength in observable anymore. This suggests that the transition from " G ≲ 0 " to " G ≪ 0 ", i. e., the transition from an interface stress sensitive bonding strength to an insensitive one, respectively, may rather be abrupt than continuous.
Finally, considering the measurements of interface reflectivity, it is suggested that this quantity is determined by surface roughness rather than interface stress. This is proposed to be due to an increased emergence of stray light at rough surfaces which known to be strongly correlated with increased reflectivity [72]. However, results may vary for shorter wavelengths which may be sensitive towards remaining bonding interface's nano-asperities.

Conclusion and outlook
In this paper, our motivation was the investigation of direct bonding of large area and stiff fused silica samples, because of its important role for state-of-the-art large-aperture-optics applications. We investigated the effects of firstly, residual stresses introduced by using a convex vacuum chuck during the contacting step and secondly, the application of a static load during the annealing step on interface properties. Our key finding is that in case of a sufficiently smooth surface neither residual stress nor static load, which both are conditions that determine the interface stress, affect bonding strength in a significant way. We propose that a consideration of change in Gibbs free energy, dG, is sufficient for estimating the resulting bonding strength and its spatial distribution, where dG is determined by surface energy and interface stress. Note that interface stress is correlated with contacted area. In doing so, we introduce three cases: spontaneous delamination prior to the annealing step, interface stress sensitive bonding strength, and interface stress insensitive bonding strength, see Fig. 14. Another finding is that, in case of visible light, interface reflectivity is determined by surface roughness rather than interface stress.
When applying direct bonding, or more precisely, PAB, for the manufacturing of modern optical and optomechanical applications, typically the objective is to obtain a high bonding strength with large area homogeneity as well as a very low reflectivity. In practical terms, as demonstrated, this desired insensitive bonding strength (of about one third of the bulk material strength) and very low reflectivity (of about 10 −6 ) can be achieved even at annealing temperatures as low as 200 • C by applying a CMP treatment on both surfaces in order to gain an RMS surface roughness of S q = 0.2 nm. Furthermore, as under these conditions the application of a static load during the annealing step does not affect bonding strength, there is no need using specialized (and therefore Fig. 14 Phase diagram which gives the relation between change of Gibbs free energy, dG, on a microscopic scale, and the resulting spatial distribution of the bonding strength, f , on a macroscopic scale. It is irrelevant if the interface stress (here given as relative value of the bulk strength) is determined by internal or external conditions, i. e., by residual stresses or a static load during the annealing step, respectively cost-intensive) heat-resistant pressure mountings anymore. This greatly improves manufacturability.
Although we quantitatively derived a condition to prevent undesired spontaneous delamination prior to the annealing step, see Eq. 4, further research is needed to predict the exact conditions of the transition from interface stress sensitive to insensitive bonding strength. Note that this transition has not been described in the literature before. Therefore, it will be interesting for the authors to see if this transition can be replicated by other research groups, furthermore, as soon as more data is available, to see if it can be characterized within ehrenFesT phase transition theory [73] or if it is just a saturation effect. Since this paper shows results for large area and stiff fused silica samples, it would be interesting to see if also small or thin samples show this effect.
The authors expect that further research would eventually allow the prediction of the bondability of substrates where residual stresses play a much more pronounced role than it is the case in today's applications and where the application of a static load during the annealing step is likewise impracticable. One conceivable application is direct bonding of thin, initially planar wafers on substrates which are intentionally non-planar, such as optical lenses or mirrors. This is particularly of technological interesting if the wafer's surface exhibits optically relevant structures (for example those which, as a state of the art, can only be inscribed on planar substrates) because it would allow realizing components which combine several optical properties. This and related ideas are proposed to even further enhance the manufacturing possibilities for novel optical applications.