A novel first-order nonparaxial scalar theory for calculating the angular scattering that is caused by the interface roughness in an optical multilayer was proposed. As in the case that the interface roughness is moderate, the analytic expressions of angular-resolved scattering for a typical p-layer design were derived. Notably, these formulas are general because they do not depend on the prior restrictive hypothesis for the correlation degree of the various interfaces in a stack. In order to verify the theory, the formulas in the case of single-surface are presented and are exactly identical to those of the generalized Harvey-Shack theory. Also, their smooth-surface approximations are the same in form as those given by the typical first-order vector perturbation theories and are validated by numerically comparing with the typical vector theory for three representative multilayer design types with slightly rough interfaces. In addition, the usability of the novel theory in the case of moderate roughness is discussed by comparing this theory to the typical theories for optical coatings at different roughness levels.
© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
Light scattering in optical thin-films remains a crucial issue and is still limiting the development of advanced optical applications [1–3]. The scattering usually originates from the structural irregularities at the interfaces and in the bulk of the materials . However, in most cases, the roughness-induced light scattering is predominant for scattering levels of stacks. To understand the origins and to further control the magnitude of off-specular optical losses in thin film technologies [5,6], a lot of theoretical tools in vector and scalar have been developed.
Based on the first-order vector perturbation approach, several vector theories have been presented by Bousquet et al. , Elson et al. [8,9], Amra et al. , etc. These theories have been widely used to predict the angular behavior of surface scattering in optical multilayers. It is worth noting that all these theories lead to identical results , even though the derivations are quite different. These theories exhibit some common advantages : firstly, they are general without a prior restrictive hypothesis for the interfacial correlation properties. Secondly, they include the polarization properties of the incident and scattered light. On the basis of these theories, a lot of works have been carried out by Amra et al. for the sake of solving the inverse scattering problem [13–15], that is, retrieving interface irregularity information from light scatter data. Recently, there are numerous literatures on dealing with ellipsometry measurement , localized defects , multiscale roughness , multiple scattering , and light scattering reduction , which further promotes the understanding of light scattering in thin film systems.
However, there is by far no obvious progress of scalar methods. In general, scalar theories are derived from the classical Beckmann-Kirchhoff (BK) theory . Inevitably, their paraxial assumption limits their abilities to handle wide-angle scattering and large angles of incidence accurately. For example, the most prominent scalar theory, presented by Eastman [22,23] and extended by Carniglia , only permits the calculation of the total amount of scattered light. In recent years, an empirically modified BK theory  that is nonparaxial has been developed by Krywonos et al. based on the typical generalized Harvey-Shack (GHS) theory , and its prediction is in good agreement with experimental wide-angle scatter measurement for surfaces with moderate roughness at arbitrary incident. This provides the motivation for developing a nonparaxial scalar scattering (NSS) theory to predict the angular surface scattering from optical multilayers, as reported in this paper. It is worth noting that the above studies are performed all for quite low-loss high-quality coatings. However, not all applications of interest satisfy the smooth-interface limitation especially at short wavelengths. To our knowledge, both existing vector and scalar theories adequately describe the behavior of scattering in multilayer coatings with moderately rough interfaces, which is also the interest of this paper.
This paper is organized as follows. In Sec. 2, the detail derivation of the NSS theory that is nonparaxial and able to handle wide-angle scattering at arbitrary incident angles for an optical multilayer with rough interfaces is given. The single surface and smooth approximations for the NSS theory are presented in Sec. 3, as well as the predictions for a narrow-band filter, an antireflective coating, and a 24-layer mirror. Modeling results for a single-layer coating and the 24-layer mirror at different roughness levels are presented and discussed in Sec. 5, followed by a conclusion in Sec. 6.
2. First-order NSS theory for surface scattering from typical p-layer design
2.1 Amplitudes of fields scattered by moderately rough interface within multilayer
A typical multilayer stack is shown in Fig. 1. There are s interfaces in this stack, forming the boundaries of p layers, where s = p + 1. Assume that all media i are homogeneous in their bulk, each bulk with a complex refractive index ni = vi – iχi and a thickness ei. The surrounding medium such as air has an index of n0 and the index of the substrate is ns. The ith interface separates media i and i + 1. Expect for interface i, all the interfaces are considered to be perfectly smooth. The ith interface is rough and presents a random height irregularity hi(x, y) that represents the displacement of the real interface from its mean position.
In our model, several assumptions for the interfacial irregularity are given as follows. Firstly, the rms value of hi is not larger than the wavelength λ of the incident light. Secondly, the slope of hi is much less than one. Thirdly, the presence of hi does not affect the indices and thicknesses of media i and j where j = i + 1. Besides, scattering induced by local defects (dust, pits, scratches, etc.) as well as multiple scattering are not taken into consideration, whereas multiple reflections in multilayered media are taken into account.
Clearly, there would be scattered light produced in whole space when a monochromatic plane wave strikes the multilayer. To illustrate the definition of quantities for this phenomenon, a diagram is drawn in Fig. 2. Our objective in this section is to derive the formulas of a scalar method that are used for calculating the amplitudes of the fields, denoted by Ei(d) ± , which are scattered in the entrance and exit media by the microstructural irregularity at the ith interface of the multilayer.
2.1.1 Theoretical method
To derive the expressions for the amplitudes of the scattered fields Ei(d) ± , we hereby follow a similar method that was used for Raman scattering by multilayers . The foundation of the derivation is the initially scattered fields (Fig. 3) generated at interface i.
As shown in Fig. 3, there are four initially scattered fields Eik(s) (where k = 1,2,3,4 are the indices of the fields) at interface i that are induced by two main incident fields Ei1(0) and Ei4(0) on the interface. Based on the modified BK theory, the amplitudes of the initially scattered fields are calculated in Sec. 2.1.2, which is the key to the NSS theory. Regarding the main field, its calculation is carried out in Sec. 2.1.3. It is worth noting that when the rms value of hi are far less than λ, the main field is determined as the stationary field inside a multilayer that is free of irregularities. However, for large roughness values, the adjustment of the main field for the scattering loss is extremely essential.
Four initially scattered fields are firstly added at the two sides of the ith interface:
Next, the fields Ei(s) ± leave the ith interface, continue to propagate through the rest of the film, and finally flow into the two outer semi-infinite media. In this process, the scattered fields obey the usual rules of field propagation in a multilayer. Fortunately, the rules have been formulized by Eastman in terms of a standard matrix solution . Besides, multiple reflections, introduced by multilayered media, are also considered that cause portions of each field to pass out of the multilayer to the entrance and exit media (as shown in Fig. 4), i.e.,
Note that the functions Pik do not depend on the interfacial irregularities but on the ideal multilayer characteristics (thickness and index of each layer), as well as on the observation conditions (wavelength, incidence and scattering angles). After Eik(s), Eik(0), and Pik are obtained, the expressions for the amplitudes of the scattered fields Eik(d) are easily given in Sec. 2.1.5.
In what follows we sum up every step of the calculation.
2.1.2 Amplitudes of initially scattered fields Eik(s)
Recalling the scattering coefficient, it is defined as the ratio of the scattered filed to the field reflected in the specular direction by a perfectly smooth surface . Utilizing this definition, the amplitudes of the initially scattered fields Eik(s) at interface i can be expressed as:
It should be pointed out that Eq. (5) is for reflection scattering, while sinθ0 of Eq. (5) must be replaced by n0sinθ0/ns for transmission scattering. Notice that here we are more concerned with the modification conditions in the modified BK theory rather than the improved closed-form solutions that only exist for surfaces with Gaussian autocovariance (ACV) functions. For this reason, the NSS theory does not depend on the kind of function (Gaussian, exponential, etc.) that describes the correlation function at interfaces.
In addition, aik in Eq. (3) refers to the reflected or transmitted field in the specular direction, and ηikhi in Eq. (4) labels the optical path difference between the incident and scattering rays. The respective terms aik and ηik are written as:Eqs. (7)-(9), i.e.,
2.1.3 Amplitudes of the main fields Eik(0)
To make it possible for our model to satisfy the law of energy conservation, it is quite necessary to adjust the main field amplitude for light scattering loss. Based on the well-established method of matrix calculation (see Eqs. (39)-(45) of ), the amplitudes of the main fields Eik(0) can be calculated by recurrence from the values Es(0), i.e.,26], the terms dmk that cause the reduction of the main field can be expressed as:
Note that the functions PSD and ACV are Fourier transforms with each other.
2.1.4 Field propagation functions Pik
By utilizing Eq. (10), the propagation function Pik for the scattered field Ei(s) ± can be given by:
Additionally, the terms Wi, Ai, Bi, and Ci satisfies:
2.1.5 Amplitudes of scattered fields Eik(d)
Equations (1)-(2) and the knowledge of Eik(s), Eik(0), and Pik lead to the following relationship for the amplitudes of the scattered fields Ei(d) ± , i.e.,
By solving the above-mentioned problems, we can calculate the Ei(d) ± field amplitudes and give a simple analytic expression of Eq. (23). Notice that Eq. (3) should be valid for moderately rough interfaces because it is directly gleaned from the modified BK theory, and the adjustment for the main field amplitude has been made in Eq. (12). Thus Eq. (23) should be suitable for a multilayer where only the ith interface is moderately rough.
2.2 Expressions for angle-resolved scattering (ARS)
In the previous section we dealt with the amplitudes of the fields that are scattered by interface i in each of the two media that surround the multilayer. We now concentrate on the case of all interface i consisting of the roughness hi(x, y). Based on the assumption that multiple scattering is neglected, the total scattered fields are simply the sum of the fields scattered by all interfaces in the p-layer stack, which are given by:
The energy scattered in a solid angle element around the direction (θs, ϕs) is proportional to:Eq. (4) leads to:
Compared with Eq. (26), Eq. (27) directly describes the correlation properties of the various initially scattered fields. In order to remove random variables hi, we take the expected value for Eq. (27) and obtain:Eq. (28), the expectation term in Eq. (29) can be written as:
Further, we denote by Hik;jl the expectation term in the right side of Eq. (30). Supposed that both hi′ and hj″ are normal random variables, Hik;jl can be derived by:Eq. (31) is in the same form as the transfer function in the GHS theory (see Eqs. (32)-(35) in ). Similar to the transfer function, Eq. (31) can be written in the from:
According to the modified BK theory, the square amplitude of the total scattered field equals to the scattered radiance defined as radiant power per unit solid angle per unit projected source area . A quantity describing the angular scattered intensity is the angle-resolved scattering (ARS) defined as the power ΔPs scattered into a small solid angle ΔΩs that is normalized to the solid angle and the incident power Pi . By multiplying the radiance by cosθs and integrating over the illumination area As, the expressions for the ARS in reflection and transmission for the multilayer can be given by:
The Bidirectional scattering distribution function (BSDF) is related to the ARS by BSDF = ARS/cosθs. From Eq. (37), we can see that the roughnesses hi of the different interfaces enter into the ARS expressions, partly through their autocorrelation and cross-correlation functions and partly through their total and band-limited rms roughness value. The cik ± terms only rely on σi,rel and the ACVij terms only exist in Gik;jl. In addition, Eq. (37) is more general compared to those of the typical scalar theories because it is not necessary to make a prior hypothesis on the ACVij terms. However, in order to deal with the calculation of Eq. (37), it is necessary to know ACVij in all i and j. The interface roughness of designs would be assumed to be fully correlated in the following sections for simplicity. This assumption means that all interfaces have the same shape as the substrate profile, yielding ACVij = ACVii = ACV for all i and j.
By comparing Eq. (37) with Eqs. (22)-(23) of  (for the popular Bousquet’s vector scattering (BVS) theory that are only applied to optical multilayers with smooth interfaces), one can notice that the scalar ARS distinctly differs from the vector ARS, owing to the fact that the scalar ARS, i.e., Eq. (37) is derived based on four initially scattered fields Eik(s) (k = 1, 2, 3, 4) at each of the interfaces within a multilayer. For comparison, the vector ARS is derived based on two initially scattered fields on the two sides of each interface. It is also noticed that in the case of the fully correlated interface roughness, the scalar ARS does not meet the linear relationship with the interface PSD, which differs from the vector ARS. Besides, the polarization properties of the incident and scattering light are not included in the scalar ARS.
3. Single-surface and smooth-interface approximations to the first-order NSS theory
In previous works, researchers have derived analytic expressions for the ARS of a moderately rough surface (see Eqs. (32)-(36) of ) and that of an optical coating (see Eqs. (22)-(23) of ). However, their results are confined with the assumption that the number of interfaces is one or that the interfaces are smooth enough. It was not existent for such a theory that can be used to predict scattering behavior for multilayered moderately rough interfaces. Due to the fact that multiple scattering is very remarkable in such a multilayer, it is almost impossible to sufficiently confirm the first-order NSS theory by experimental measurement. In this section, Eq. (37) in the single-surface and smooth-surface approximations was discussed and analyzed.
3.1 Case of single surfaces
In this subsection, we are concerned with an uncoated random rough surface. The roughness is described by function h(x, y), whose rms value can be the same order of the incident wavelength λ. The media of indices n0 and ns are separated by this surface. The angular scattered intensity is still calculated by Eq. (37). For the case of s = 1, Eq. (37) in refection is approximated as:
Correspondingly, B01;01 and G01;01 are still calculated by Eqs. (34) and (35), i.e.,
It is noticed that Eqs. (39)-(41) are exactly the same as Eqs. (32)-(36) of  given by the GHS theory. Consequently, the NSS theory can be applied for uncoated surfaces even if the roughness is moderate. Meanwhile, it is also validated that the initially scattered field amplitude is given by Eq. (3) for moderately rough interfaces.
3.2 Case of smooth-surfaces
By assuming the interfacial roughness to be small compared to the incident wavelength (i.e., hi << λ for all i), the scattering function Eq. (37) can be approximated as:Eqs. (13)-(15) are omitted by setting dik = 1 at all i and k. In addition, σi,rel is set to be equal to σi,tot in Eq. (42) for simplicity.
By comparing Eq. (42) with Eqs. (22)-(23) of , it is found that the approximated scalar ARSs are in the same form as the vector ARSs given in the BVS theory. Especially in the case of the fully correlated interface roughness (PSDij = PSDii = PSD), Eq. (42) is simplified to:Eq. (42) and Eqs. (22)-(23) of , Figs. 5(a)-5(c) are plotted for three representative design types (where the interfaces are isotropic and smooth, of which the rms roughnesses are in the range of validity of the BVS theory) to illustrate how close the predictions of the NSS and BVS theories are to each other. For the prediction of the BVS theory, both the incident and scattered light fields are assumed to be s polarized and remain in the plane of incidence.
Figure 5(a) is for a narrow-band filter SH(LH)2L6(HL)2H (cited from ); Fig. 5(b) is for an antireflective coating LH2SH2L; Fig. 5(c) is for a dielectric reflector S(HL)11HL2 (cited from ), where S represents the substrate, H and L are for high- and low-index quarterwave optical thickness layers, respectively. The parameters used for the calculation of Fig. 5(a) are as follows: nH = 2.3 and nL = 1.3. The interfaces are described by the sum of a Gaussian and an exponential , i.e.,Fig. 5(b) are identical to those of Fig. 5(a), except that nH = 1.90 and nL = 1.38. In Fig. 5(c), the parameters nH and nL are set to be nH = 2.3 and nL = 1.45. The autocovariance function of interface is :
It can be seen from Fig. 5(a) that the NSS prediction agrees well with the prediction of the BVS, although there is a slight disagreement at near specular directions. The BVS curve here is consistent with the theoretical curve in Fig. 17 of . The excellent agreement between the NSS and BVS curves is also found in Fig. 5(b) especially at the specular direction. In Fig. 5(c), the NSS method still predicts the curve that is almost exactly identical to the BVS curve. The difference can be seen in the regions of 60°~150° and −60°~−150° (that is, the angular range away from the specular direction), but their weight is negligible for the total scattering. Besides, it is noticed that the predicted curves in Fig. 5(c) agree well with the theoretical curve in Fig. 11 of . However, compared to the measured curve in Fig. 11 of , the NSS curve of Fig. 5(c) is decreased at the near specular direction. The decreased scattering may be the result of neglecting scattering of surface defects or dust particles in the NSS model. In the following, we further give Fig. 6 for the reflector as the example to seek how the validity of the first-order NSS theory varies with (a) the incident angle θ0 and (b) the correlation length lc2.
It is clear from Fig. 6 that the NSS prediction agree well with the prediction of the BVS in all cases. However, there are slight differences mainly in the angular range that is away from the specular direction. These differences are acceptable owing to their negligible proportion for the amount of scattered light. And the differences are increased with either θ0 or lc2 for different reasons. The reason why the differences increase with the increasing of θ0 is that large values of θ0 cause the lateral displacement of various reflected waves inside the mirror, which makes the interfacial correlation properties changed. Small values of lc2 result in the violation of the small roughness slope assumption, which makes the differences increased.
To summarize, the expression of the NSS theory in the case of smooth-interfaces, i.e., Eq. (42), is not only in the same form as those given by the BVS theory (that is, Eqs. (22)-(23) of ), but also their predictions for different stacks are in good agreement especially at small incident angles and large correlation lengths. This indicates the validity of the NSS theory for optical coatings with smooth interfaces. It is worth pointing out that the BVS theory is also first-order limited just like the NSS theory. The reason why the first-order theories are valid here is that high-quality coatings have quite low scattering levels.
4. Interface roughness dependence of the validity of the first-order NSS theory
In this section, we provided examples to show the dependence of the validity of the NSS theory on the interface roughness. To this effort, the results of ARS calculations using the NSS, BVS, and GHS theories in the following discussion were given for optical coatings at four different roughness levels. It is worth pointing out that the NSS theory is only used to calculate the scattering from bare substrate or top surface of a stack. The purpose is to show the regularity and magnitude of scattering from moderately rough surfaces as a reference to judge whether the prediction of the NSS theory is reasonable in the case of moderate roughness. Two typical designs were analyzed: a single low index layer and a 24-layer mirror.
4.1 Single low index layer coating
Consider a single-layer coating deposited on an opaque substrate. In this case, the impact of multiple scattering is relatively low because only two rough interfaces are involved. Besides, the parameters used for simulations are as follows: n0 = 1.0, n1 = 1.3, and ns = 1.52. e1 = λ/4 at i0 = 0° and λ0 = λ = 633 nm. For simplifying the calculation, we suppose that the interfaces are described by a one-dimension Gaussian autocovariance function:
The results of ARS calculations using the NSS theory for the single layer coating at different roughness levels (σtot/λ = 0.01, 0.05, 0.1, and 0.5) are given in Fig. 7. It is very intuitive that the intensity of scattering at the near specular direction increases firstly with the increasing of the roughness and then decreases slightly with the further increasing of the roughness, while the large angle scattering keeps increased with the increasing roughness. This redistribution of power is in conformity with the law of energy conservation.
Moreover, the curves of ARS calculated by the NSS, BVS, and GHS theories for the single layer coating are plotted in Fig. 8.
In Fig. 8(a1), for the smoothest interface case under investigations (σtot/λ = 0.01), i.e. the condition that the interfaces are smooth enough is satisfied, the prediction of the NSS agrees well with the BVS prediction. In theories of single surface scattering, the criterion to evaluate whether or not an optical surface is smooth is : σtotcosθi/λ≤ 0.02. We suppose that this criterion is valid likewise for interfaces inside a multilayer. The GHS theory predicts the higher ARS curve than the NSS and BVS theories, which implies that it is possible to reduce scattering in multilayers by taking a proper design. The reason why the scattering is reduced may be that the optical factors, cik± , at different interfaces are in opposite phase, which leads to a destructive interference effect .
Figure 8(a2) for σtot/λ = 0.05 shows that the BVS predicts slightly stronger scattering than the NSS over all scattering angles. This is reasonable because the smooth interface lamination is slightly violated. Except this, the prediction in Fig. 8(a2) is quite similar to that in Fig. 8(a1) in terms of the curve trend. With further increasing of the roughness in Fig. 8(a3) for σtot/λ = 0.1 and Fig. 8(a4) for σtot/λ = 0.5, the differences between the BVS and NSS increase quickly. The reason is that the ARS predicted by the BVS is proportional to the PSD of interface and thereby to the square roughness while the predictions of the NSS at low scattering angles are decreased especially at large values of roughness. As a result, the BVS model predicts the highest but relatively reasonable curve in Fig. 8(a3), while it is notably failed for the BVS prediction in Fig. 8(a4) because of its TS value greater than one.
Although the smooth interface approximation is violated seriously in Fig. 8(a4), the reason why we consider that the NSS theory show reasonable accuracy is based on the fact that the NSS results should be close to those of the GHS. The reflection scattering of a multilayer is the superposition of the two components: the first consisting of reflection scattering of the top-surface and the second consisting of portion of transmission scattering of the top-surface that is then reflected back to the incident medium. Due to low reflectance in the single low-index layer coating, the scattering level is predominated by the reflection scattering at the top-surface, while the contribution of the other interfaces is relatively low. Moreover, it is noticed that the results in Figs. 8(a1)-8(a3) are consistent with the conclusion of  that a scattering-reduction effect can be obtained by coating a rough surface with an antireflection layer. However, the case of quite rough surfaces was not discussed in , and according to the results in Fig. 8(a4), we think that this method of reducing scattering is not suitable for very rough surfaces.
To seek how the ARS varies with the roughness, the results of normalized ARS are shown in Fig. 8(b). In Fig. 8(b1), similar to Fig. 8(b3), the width of the calculated curve increases with the increasing roughness. However, in Fig. 8(b2), the BVS result does not rely on the roughness. This means that the variation of roughness only causes the change in the magnitude of the ARS predicted by the BVS without affecting its profile. Besides, Fig. 8(b4) for the case of θ0 = 45° and lc = 2λ shows that different roughness levels lead to slightly shift in the location of ARS curve peaks corresponding to the specular direction.
4.2 24-layer dielectric mirror
To reveal the interface roughness dependency of the validity of the first-order NSS theory, it is necessary to discuss the influence of the roughness on the scattering properties of a system with several layers. The results of ARS and normalized-ARS for the 24-layer mirror are shown in Figs. 8(a) and 8(b), respectively. Note that a 1-D Gaussian ACV function was used to describe the interfaces within the mirror for simplicity purpose.
Figure 9(a1) shows that the NSS and BVS results agree well with each other, as expected. Since the validity of the BVS theory is largely restricted to smooth interface criteria, the differences between the BVS and NSS predictions are dramatically increased with the increasing roughness until extremely remarkable, as shown in Figs. 9(a2)-9(a4). Regarding the GHS, its prediction in Fig. 9(a) is much smaller than the predictions of the NSS and GHS, which differs from that in Fig. 8(a). This can be explained by the fact that the reflectance of the mirror is rather higher than that of the single layer coating, leading to scattering from interface i ≠ 0 that contributes mainly to scattering level instead of scattering from interface i = 0. Hence, we conclude that the NSS theory predicts ARS curves in reasonable values in Fig. 9. Also, it can be summarized some features for scattering of a mirror from Fig. 9(b). In Fig. 9(b2), similar to Fig. 8(b2), the predicted curve is not related to the roughness. In Fig. 9(b1), the angular width of the predicted curve decreases with the increasing of the roughness, which is opposite to that in Fig. 9(b3). In general, rougher surfaces contain larger slopes of the surfaces leading to wider angle ranges of scattering, which is for the case of single surface scattering. In Figs. 8(b1), for a surface having a single low index layer, the presence of this effect may be due to the predomination of the top-surface scattering for the reflected scattering. However, in Fig. 9(b1), for a surface with a high-reflection film, the opposite effect may be caused by the fact that the scattered field of which the propagation is along the specular direction is easier to be reflected back into the incident medium. Besides, in the case of the oblique incidence shown in Fig. 9(b4), there are still shift in the location of the normalized ARS curve peak.
In conclusion, a novel first-order nonparaxial scalar scattering theory for multilayer coatings in which the interfaces are assumed to be rough was proposed. In this theory, neglecting multiple scattering and based on the modified BK model, analytic expressions for the ARS of a typical p-layer design were derived. The theory, as a result of scalar treatment, first gives access to the angular scattering for multilayer coatings and does not contain a prior restrictive hypothesis for the correlation degree of the various interfaces. Moreover, we show that the formulas in the single surface approximations are exactly identical to those given by the GHS theory and they in the smooth-surface approximation are in the same form as those given in the BVS theory. Numerical graphs for three representative design types are plotted to confirm the validity of the NSS theory for multilayer coatings assuming slightly rough interfaces. Results indicate that the NSS theory is valid especially for small incident angles and large correlation lengths. Also, the simulations for two multilayer design types at four different roughness levels reveal that the first-order NSS theory enable to predict the magnitude or regularity of scattering in multilayer coatings with moderately rough interfaces, which demonstrates its unique advantage compared with typical vector and scalar theories. Besides, the prediction of the NSS theory fits measured data roughly for a 24-layer mirror with smooth interfaces. However, the NSS theory is first-order limited since multiple scattering is neglected, which makes it quite difficult to verify the usability of the NSS theory in the case of moderate roughness by experiment, which is the aim of our future work.
Sichuan Science and Technology Program (Grant No. 2018JY0550) and Anhui Province Key Laboratory of Non-Destructive Evaluation, Hefei ZC Optoelectronic Technologies Ltd (Grant No. CGHBMWSJC06).
1. V. K. Sakharov, “Model of lock-in in a ring laser and a semiconductor laser gyro,” Tech. Phys. 56(8), 1135–1141 (2011). [CrossRef]
2. H. Miao, H. Yang, R. X. Adhikari, and Y. Chen, “Quantum limits of interferometer topologies for gravitational radiation detection,” Class. Quantum Gravity 31(16), 165010 (2014). [CrossRef]
3. S. Zeidler, T. Akutsu, Y. Torii, E. Hirose, Y. Aso, and R. Flaminio, “Calculation method for light scattering caused by multilayer coated mirrors in gravitational wave detectors,” Opt. Express 25(5), 4741–4760 (2017). [CrossRef] [PubMed]
5. J. A. Thornton, “Structure and topography of sputtered coatings,” Annu. Rev. Mater. Sci. 7, 239–260 (1977). [CrossRef]
6. W. M. Tong and R. S. Williams, “Kinetics of surface growth: phenomenology, scaling, and mechanisms of smoothening and roughening,” Annu. Rev. Chem. 45(1), 401–438 (1994). [CrossRef]
7. P. Bousquet, F. Flory, and P. Roche, “Scattering from multilayer thin films: theory and experiment,” J. Opt. Am. 71(9), 1115–1123 (1981). [CrossRef]
8. J. M. Elson, “Diffraction and diffuse scattering from dielectric multilayers,” J. Opt. Soc. Am. 69(1), 48–54 (1979). [CrossRef]
9. J. M. Elson, “Surface scattering in optical interference coatings,” J. Opt. Soc. Am. 66, 230–234 (1974).
10. C. Amra, “First-order vector theory of bulk scattering in optical multilayers,” J. Opt. Soc. Am. A 10(2), 365–374 (1993). [CrossRef]
11. S. Schröder, T. Herffurth, H. Blaschke, and A. Duparré, “Angle-resolved scattering: an effective method for characterizing thin-film coatings,” Appl. Opt. 50(9), C164–C171 (2011). [CrossRef] [PubMed]
15. C. Amra, “Light scattering from multilayer optics. I. Tools of investigation,” J. Opt. Soc. Am. A 11(1), 197–210 (1994). [CrossRef]
19. G. Xu, P. G. Stegmann, S. D. Brooks, and P. Yang, “Modeling the single and multiple scattering properties of soot-laden mineral dust aerosols,” Opt. Express 25(24), A990–A1008 (2017). [CrossRef] [PubMed]
20. J. Zhang, H. Wu, H. Jiao, S. Schröder, M. Trost, Z. Wang, and X. Cheng, “Reducing light scattering in high-reflection coatings through destructive interference at fully correlated interfaces,” Opt. Lett. 42(23), 5046–5049 (2017). [CrossRef] [PubMed]
21. P. Beckmann and A. Spizzichino, “The Scattering of Electromagnetic Waves from Rough Surfaces,” (Pergamon Press, 1963).
22. J. M. Eastman, “Scattering by All- Dielectric Multilayer Bandpass Filters and Mirrors for Lasers,” in Physics of Thin Films. Advances in Research and Development 10, 167–226 (1978).
23. J. M. Eastman, “Surface Scattering in Optical Interference Coatings,” J. Opt. Soc. Am. 66, 230–234 (1974).
24. K. Carniglia, “Scalar scattering theory for multilayer optical coatings,” Opt. Eng. 18(2), 104–115 (1979). [CrossRef]
25. J. E. Harvey, A. Krywonos, and C. L. Vernold, “Modified Beckmann-Kirchhoff scattering model for rough surfaces with large incident and scattering angles,” Opt. Eng. 46(7), 078002 (2007). [CrossRef]
26. A. Krywonos, J. E. Harvey, and N. Choi, “Linear systems formulation of scattering theory for rough surfaces with arbitrary incident and scattering angles,” J. Opt. Soc. Am. A 28(6), 1121–1138 (2011). [CrossRef] [PubMed]
27. R. A. Craig, G. J. Exarhos, W. T. Pawlewicz, and R. E. Williford, “Interference-enhanced Raman scattering from Ti02/Si02 multilayers: measurement and theory,” Appl. Opt. 26(19), 4193–4197 (1987). [CrossRef] [PubMed]
29. S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, and J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Express 19(10), 9820–9835 (2011). [CrossRef] [PubMed]