## Abstract

Fresnel reflectance and transmittance coefficients of a thin film system consisting of an arbitrary number of layers are expressed explicitly in the form of a power series of Fresnel coefficients for individual boundaries and phase terms for the individual films. The series is based on the evaluation of all possible paths light can pass through the system. However, the series is written as consolidated, i. e. all paths corresponding to the same powers are represented using a single term in the series, with multiplicity which is a simple product of binomial coefficients. This result is used to express the normal reflectance of a thin film system with arbitrarily correlated randomly rough boundaries and it is shown that such approach can be computationally efficient in practice.

© 2014 Optical Society of America

## 1. Introduction

Interaction of light with a layered medium, such as a system of thin films, can be viewed as a series of interactions (reflections and refractions) at layer boundaries and propagation through the individual films. However, this view is mostly relegated to introductory optics textbooks where the enumeration of all the interactions may be used for the derivation of reflectance and transmittance of a single thin film. In most practical calculations in optics of thin films a matrix formalism is used instead. Several such formalisms exists, e. g. a 2 × 2 or 4 × 4 formalisms [1] based on the formalism developed by Jacobson [2], or an admittance-based approach [3], nevertheless, they all express the relation between incoming and outgoing light wave amplitudes at one end and the other end of the thin film system using a matrix. Another frequently used approach to multi-layer systems is the recursive or inductive application of formulae for a single layer [36].

Both the matrix formalism and recursive formulae are excellent for the calculation of reflectances and transmittances of multi-layer systems with smooth boundaries. However, in the case of systems with rough boundaries the situation is different. The only theoretical approach that fits well into the matrix formalism is the effective medium approximation (EMA) [79]. This approach is very popular, also due to its simplicity, however, it is valid only for very fine roughness with high spatial frequencies [1]. More advanced approaches used in optical characterisation of rough thin films and their systems include in particular the scalar diffraction theory (SDT) [1017] and perturbation theories [1826] such as the Rayleigh–Rice theory (RRT) [2732]. An overview of the theoretical approaches can be found in [23, 33, 34]. New methods of calculation of light scattered from rough surfaces and thin films are still being actively developed. However, the SDT and second-order RRT remain the state of the art in the field of ellipsometric and spectrophotometric characterisation of rough thin film systems using specularly reflected light. Both SDT and second-order RRT are able to describe roughness outside of the validity range of EMA, each for a different range of magnitudes and spatial frequencies of the roughness. But both also express the reflectances and transmittances globally for the entire thin film system and are not factorable to matrices. This means that they cannot be included within the matrix formalism. Note that a simple matrix approximation of the RRT, applicable in certain special cases, has been formulated [1, 30]. In the case of SDT, an approximate approach was published, based on averaging of individual elements of the matrices entering the matrix formalism [14, 15]. However, this approximation is too rough and often leads to erroneous results. An improved but still approximate method was based on the approximation 〈a/b〉 ≈ 〈a〉/〈b〉 in the recursive formulae for reflectance of a multi-layer thin film [35, 36] (〈...〉 denotes averaging).

In the SDT, it is advantageous to express the reflection coefficients of a layered system using a power series in the Fresnel coefficients and phase terms for individual boundaries and layers instead [11,12,37,38]. However, neither the matrix formalism nor the recursive formulae express the Fresnel coefficients of multi-layer systems as power series. They express the coefficients as relatively complicated rational functions of Fresnel coefficients of individual boundaries and phase terms describing the propagation inside the layers. The sequence-of-interactions view, on the other hand, leads to a power series suitable for the SDT. As noted above, writing the series for a single layer is trivial. Utilisation of the series within the SDT leads to the well-known expressions for the reflectance of rough thin films [12]. Multiple layers, however, represent a challenge.

All possible sequences of reflections and refractions that start with light wave falling onto the uppermost boundary from the ambient and ending with the light wave leaving the uppermost or lowermost boundary towards the ambient were enumerated by Butler in a mathematical journal probably not read by many physicists [39]. The connection to Catalan paths was also noted. Catalan paths for systems from two to six layers are enumerated for instance in the on-line encyclopaedia of integer sequences [4044]; they converge for an infinite number of layers to Catalan numbers [45].

For any multi-layer system the total number of possible ‘paths’ grows exponentially with the number of interactions. However, many of the paths lead to the same combinations of Fresnel coefficient and phase term powers. An illustration of two simple sequences of reflections and refractions that produce equivalent terms in the series is shown in parts (a) and (b) of Fig. 1. Therefore, it would be useful to write a consolidated power series which contains each such term only once, with the corresponding multiplicity as a mere combinatorial factor. And it would be useful to group these terms further by the powers of phase terms as this is also useful in SDT calculations. A path that corresponds to different powers of Fresnel coefficients but the same phase term powers is illustrated in Fig. 1(c).

Fig. 1 Illustration of equivalent light paths through a two-layer system. Paths (a) and (b) lead to identical powers of both the Fresnel coefficients and phase terms, whereas (c) leads only to the same phase terms powers as (a) and (b); the Fresnel coefficients differ.

Both consolidations are possible and the construction of such series will be presented in this work. Furthermore, it will be shown how this series expansion can be applied within SDT to calculation of normal-incidence reflectance of a thin film system with rough boundaries.

Note that the path summation approach was used for the expression of reflection and transmission coefficients of individual inhomogeneous layers [46], although only the lowest order terms were expressed explicitly. In this work the full infinite series containing paths of arbitrary lengths will be formulated. The complete series is equivalent to the expressions obtained using the matrix formalism or recursive formulae.

## 2. Construction of the series

Since any practical calculation employing summation of the series may include only a finite number of terms, the series will be constructed from the beginning as truncated to a certain number of interaction at boundaries. Letting the number of interactions grow to infinity in order to obtain the infinite series is then trivial, whereas the limiting the infinite series to a specific number of interactions could be considerably less straightforward. It should be noted the ordering by the number of interactions is used because it is convenient for the construction of the infinite series. In practical calculations, the series can be truncated to a finite number of terms differently as will be discussed at the end of paragraph 2.2.

#### 2.1. Parametrisation

The media in the layered system will be indexed j = 0, 1,...,L, with 0 denoting the ambient, L denoting the last medium which is presumably the substrate. Index j will be referred to as medium depth for brevity. The reflection and transmission coefficients and phase terms will be then denoted j, r̂′j j, t̂′j and exp(ij) as displayed in Fig. 2, where

$X^j=2πλN^jhjcosϑ^j.$
The symbols j, hj and ϑ̂j denote the complex refractive index, thickness, and complex angle of incidence corresponding to medium j; and λ is the wavelength. No specific relations among the coefficients will be assumed because it would not bring any significant simplification while it would prevent application of the results in other situations, namely when the layer boundaries are not simple ideal boundaries.

Fig. 2 Reflection and transmission coefficient and phase term notation.

A light path which is continuous and starts and ends in medium 0 contains an odd number m of interactions at the boundaries. It is more useful to introduce path length p using the formula

$m=2p+1.$
The path length is non-negative and p = 0 corresponds to a single reflection at the 0–1 boundary. The total power of all Fresnel coefficients in the term corresponding to this path is 2p + 1 and the total power of all phase terms is 2p. Beside path length we will also define path depth d which is equal to the depth of the deepest medium the path passes through (so dL). Any two paths that differ in either p or d are obviously different.

A suitable set of parameters describing a path of depth d is then formed by non-negative integers m0, m1, . . . , md called medium counts and v0, v1, . . . vd called visit counts.

The medium count mj expresses how many times the light:

• passes through medium j downwards,
• passes through medium j upwards,
• is incident on the j–(j + 1) boundary from the side of medium j,
• is incident on the (j − 1)– j boundary from the side of medium j.
It can be seen that all the counts are equal. The count mj is also equal to the power of exp(2im) in the term corresponding to paths described by this value of mj. The medium count can also be interpreted as the number of times the light is reflected off the entire subsystem from medium j + 1 onward, which will be useful below. The path length is related to mj as follows
$p=∑j=1dmj$
and for all paths and all layer systems it holds m0 = 1.

The visit count vj expresses how many times the light:

• visits the subsystem from medium j downward,
• passes through the (j − 1)– j boundary from the side of medium j − 1,
• passes through the (j − 1)– j boundary from the side of medium j.
The first point should be understood as follows: The entire path can be divided into segments contained in media 0 to j − 1 and segments contained in media from j downward. Visit count vj is the number of the latter segments. For all paths and all layer systems it holds v0 = 1.

#### 2.2. Parameter ranges

Medium and visit counts are bound by a set of simple inequalities that follow from continuity of the path (i. e. that it cannot skip any layers). First, for all jd it holds

$vj>0, mj>0,$
while for j > d it holds vj = mj = 0. Once the light enters medium j it must pass through it at least once which leads to the inequality vjmj. Finally, since visits of medium j + 1 must be separated by at least one pass through medium j it must hold vj+1mj. The preceding two relations can be written summarily for vj or mj in a more compact form
$mj≥max(vj,vj+1) or vj≤min(mj,mj−1).$
It is evident that conditions (4)(5) are necessary, i. e. all paths must satisfy them. It may not be immediately evident that they are also sufficient, i. e. that each set of positive m1, m2, . . . , md and accompanying set of v1, v2, . . . , vd bound by these conditions corresponds to at least one possible path.

As can be seen with the help of Fig. 3, if the path

• enters medium j from medium j − 1 the same number of times as it leaves medium j to medium j − 1,
• enters medium j + 1 from medium j the same number of times as it leaves medium j + 1 to medium j,
• leaves the (j − 1)– j boundary downwards the same number of times as it is incident on the j–(j + 1) boundary from the side of medium j, and
• leaves the j–(j + 1) boundary upwards the same number of times as it is incident on the (j − 1)– j boundary from the side of medium j then it is always possible to ‘match the arrows’ in medium j close to the (j − 1)– j boundary with those close to the j–(j + 1) boundary. Hence a path without any loose ends can be always formed provided that mjvj and mjvj+1 are non-negative – which corresponds exactly to conditions (5) (if the matching is done in order left to right no cycles are formed).

Fig. 3 Relations between medium and visit counts connected with medium j and the numbers of reflections within this medium and refractions to or from this medium.

Using the second inequality in Eq. (5), ranges of visit counts can be easily expressed if the medium counts are given. The overall summation scheme for paths of positive depth can thus be constructed as follows:

$∑p=1P∑d=1min(L,p)∑m∑v$
where P denotes the maximum path length which can be either specified explicitly or determined implicitly by a termination condition, e. g. by terminating the series once the absolute value of the entire length-p term falls below a certain small value. The vectors m and v were introduced as shorthand for m1, m2,..., md−1 and v2, v3,..., vd. The m-summation is then
$∑m=∑m1=1p−(d−1)−M0∑m2=1p−(d−2)−M1⋯∑md−1=1p−1−Md−2$
while the v summation is
$∑v=∑v2=1min(m1,m2)∑v3=1min(m2,m3)⋯∑vd=1min(md−1,md).$
The symbol Mj denotes the sum of medium counts for media from 1 to j
$Mj=∑l=1jml, M0≡0$
and md = pMd−1 which follows from Eq. (3). No summation is thus carried out over md because it is uniquely determined by the path length and other mj. Similarly, no summation is carried out over v1 because v1 ≤ min(m0, m1) = 1, hence v1 = 1.

It should be noted that the ordering by path length in Eq. (6) was the only reason why any relations between ranges of individual mj appeared in Eq. (7). If this requirement is removed the upper limits of mj summations can be chosen differently, e. g. each can be chosen independently based on physical considerations. Therefore, a criterion based on the behaviour of individual reflection and refraction coefficients can be used, taking into account more reflections in specific selected layers.

The above scheme does not include the singular term corresponding to a single reflection off the uppermost boundary (p = 0, d = 0). It could be formally included in the scheme but only at the cost of impairing its clarity. Since it would still represent a special case, only made more obscure, the result would not translate well into a computer implementation. Hence it was excluded and must be added separately.

#### 2.3. Coefficient and phase term powers

The powers of Fresnel coefficients and phase terms can be expressed trivially with the help of Fig. 3. Layer j contributes the factor

$(t^jt^j′)vj(r^j′)mj−vjr^j+1mj−vj+1exp(2imjX^j)$
to the sum term. For reference, we will write out explicit forms of the summation terms for depths from 0 to 3 that occur in the series for systems consisting of up to three layers.

Depth 0 corresponds to the single term $P0=r1m0=r^1$ representing the reflection off the uppermost boundary.

Depth 1 represents reflections within the uppermost layer, i. e. the classic series

$P1(m1)=(t^1t^1′)v1(r^1′)m1−v1r^2m1exp(2im1X^1)=t^1t^1′(r^1′)m1−1r^2m1exp(2im1X^1),$
where v1 = 1 was used to obtain the latter expression from the former which corresponds formally to Eq. (10). The following formulae will already be written with this simplification.

Adding the contribution of layer 2 results in the expression of depth-2 terms:

$P2(m1,m2,v2)=t^1t^1′(r^1′)m1−1r^2m1−v2(t^2t^2′)v2(r^2′)m2−v2r^3m2exp(2im2X^2)exp(2im1X^1).$
The depth 3 terms can be obtained analogously:
$P3(m1,m2,m3,v2,v3)=t^1t^1′(r^1′)m1−1r^2m1−v2(t^2t^2′)v2(r^2′)m2−v2r^3m2−v3(t^3t^3′)v3(r^3′)m3−v3r^4m3×exp(2im3X^3)exp(2im2X^2)exp(2im1X^1).$

#### 2.4. Multiplicities

The only missing piece is the multiplicity factor for each term in the series, i. e. the number of possible paths that give raise to each unique combination of Fresnel coefficient and phase term powers. Again, we will refer to Fig. 3 from which it can be seen that in each layer two independent choices are made:

• distribution of the vj visits of medium j among the mj−1 incidences of light on (j − 1)– j boundary from the side of medium j − 1,
• distribution of the mj downward passes of light through medium j among the vj visits of this medium.

The first distribution can be done

$(mj−1vj)$
ways while the second can be done
$(mj−1vj−1)$
ways. Therefore, layer j contributes to the term multiplicity with an independent factor that will be denoted F(vj, mj, mj−1), i. e.
$F(vj,mj,mj−1)=(mj−1vj)(mj−1vj−1)$
and the entire multiplicity Cd for a depth-d summation term can be obtained simply by multiplying these factors for all d layers:
$Cd(m1,m2,…,md−1;v2,v3,…,vd)=Cd(m,v)=∏j=1dF(vj,mj,mj−1).$
For reference, the multiplicities will be again listed for depths from 0 to 3:
$C0≡1,C1(m1)=C0F(v1,m1,m0)=(11)(m1−10)=1,C2(m1,m2;v2)=C1(m1)F(v2,m2,m1)=(m1v2)(m2−1v2−1),C3(m1,m2,m3;v2,v3)=C2(m1,m2;v2)F(v3,m3,m2)=(m1v2)(m2−1v2−1)(m2v3)(m3−1v3−1).$

#### 2.5. Second consolidation and summation algorithm

The order of summation in Eq. (6) permits grouping all identical phase terms together as follows:

$r^=r^1+∑p=1P∑d=1min(L,p)∑mexp(2i∑j=1dmjX^j) Q^(m),$
where (m) is a homogeneous polynomial of degree 2p+1 in Fresnel coefficients. The ranges of mj, which are omitted here for brevity, are the same as in Eq. (6) where they were written out explicitly. The polynomial (m) is constructed using the tail part of Eq. (6), expression (10) with phase terms excluded, and formula (17):
$Q^(m)=t^1t^1′(r^1′)m1−1r^d+1md∏j=2d∑vj=1min(mj,mj−1)F(vj,mj,mj−1)r^jmj−1−vj(t^jt^j′)vj(r^j′)mj−vj.$
In this expression the term $r^d+1md$ was factored out of the innermost summation in front of the entire sum as it does not depend on any non-zero vj. Conversely, terms $r^j+1mj−vj+1$ in Eq. (10) were pulled into one level deeper summations since they belong to layer j conceptually but depend on vj+1, not on vj. The product from j = 2 to d should be understood as empty (equal to 1) for d = 1.

The full factorisation of Eq. (20) means that if the Fresnel coefficients can be treated independently, the number of arithmetic operations necessary to calculate the sum (20) grows only as p instead of the pd/d! growth which corresponds to the fully expanded sum. If the further mathematical operations can be only carried out with complete terms, the expansion is straightforward: the summations over vj are performed sequentially, in a nested manner indicated in Eq. (6), instead of factoring the sums. In paragraph 2.6 where the number of terms is discussed, the latter, more pessimistic, variant will be assumed. It should be noted that both variants coincide for a two-layer system because only one vj summation (over v2) is done.

#### 2.6. Number of terms

As noted in the introductory remarks, the number of possible paths grows exponentially with path length. Specifically for paths of depth 2, 3 and 4 it grows as 2p−1, F2p+1 and (3p − 1)/2+1, respectively, where Fn denotes the n-th Fibonacci number [4042].

In contrast to the number of possible paths, the number of terms in the consolidated series grows polynomially with p. For paths of depth 2 and p ≥ 2 it immediately follows from parameter ranges given by Eqs. (4) and (5) that the number of terms is

$⌊p2⌋(p−⌊p2⌋)≈p24.$
Here ⌊...⌋ denotes rounding down to the nearest integer (the ‘floor’ function). For paths of depth 3 the exact formula is convoluted, however, the number of terms grows as p4/72. Therefore, the number of terms in the series including all paths with length up to p, inclusively, grows as p3/12 for a two-layer system and as p5/360 for a three-layer system. Generally, for an L-layer system it grows as p2L−1. This means that while the combinatorial explosion may be much slower asymptotically in the consolidated series than in the expanded series, it is still fast if the system consists of many layers.

Furthermore, it is not possible to terminate the series after a small number of terms (small p) if L is large and if any contribution of the deeper layers is to be included because the j-th medium contributes only to paths of length j or more.

Both properties limit the usability of the series expansion, even consolidated, in the case of many-layer systems. However, in the case of systems consisting only of several layers it has a certain potential.

The reduction of the number of terms by consolidation can be observed in Table 1 which compares the number of terms corresponding to all paths of length p in systems consisting of two to five layers (paths of any depth are counted, not just those that of maximum possible depth). It can be seen that twofold reduction occurs if paths of length at least 5 are considered so the presented scheme does not lead to a significant simplification if only a handful of terms is taken into account. For paths longer than 5, however, the gain from consolidation is considerable and it increases rapidly with increasing p.

Table 1. Number of terms in the series for individual path lengths p and two to five-layer systems. Columns labelled ‘all’ represent all possible paths and correspond to Table 2 in [39]. Columns labelled ‘cons.’ represent the consolidated number of terms corresponding to the series presented in this work.

#### 2.7. Transmission

The summation scheme formulated for reflection can be adapted for transmission in a fairly straightforward manner. All paths starting with light wave falling onto the uppermost boundary from the ambient and ending with the light leaving the lowermost boundary to the ambient have depth d = L so summation over d is avoided in formulae such as (6) and d is replaced with L. The number of terms in the consolidated series grows somewhat less rapidly, though still as p2L−1.

The meaning of mj and vj changes slightly: they still denote the numbers of downward passes and refractions, however, the numbers of upwards passes and refractions are now mj − 1 and vj − 1, respectively. Figure 3 remains largely unchanged, only with upward arrow counts replaced with vj − 1 and vj+1 − 1. Definition (3) of path length p can be retained, however, the total number of interactions is 2p + 1 − L now and the total power of phase terms is similarly 2pL. Condition (4) holds for all 0 ≤ jL + 1 and condition (5) remains unchanged.

The contribution of one layer to powers of Fresnel coefficients (10) changes to

$t^jvj(t^j′)vj−1(r^j′)mj−vjr^j+1mj−vj+1exp(i(2mj−1)X^j)$
and an additional coefficient $t^L+1vL+1=t^L+1$ appears in all terms. Of course, the singular term 1 disappears.

Concerning multiplicities, the reasoning leading to expression (15) remains valid, whereas Eq. (14) needs to be modified. To obtain a path leaving the layer system through the lowermost boundary, the last incidence of light must result in a refraction, not reflection. Therefore, only the remaining vj − 1 visits can be freely chosen among the remaining mj−1 − 1 incidences and the factor F(vj, mj, mj−1), expressed by Eq. (16) for reflection, becomes

$F(vj,mj,mj−1)=(mj−1−1vj−1)(mj−1vj−1).$
With these modification it is possible to utilise the same summation approach as for reflection.

## 3. Application to scalar diffraction theory for normal incidence

The summation scheme presented above can be applied in a relatively straightforward manner to SDT for the normal reflectance of a thin film system with rough boundaries exhibiting arbitrary cross-correlations. It was shown [35,47,48] that the normal reflectance R is expressed

$R=〈r^(u)exp(ivt1)〉〈r^*(u)exp(−ivt1)〉=|〈r^(u)exp(ivt1)〉|2,$
where u = u(x,y) is the vector with L + 1 components representing the heights of irregularities of boundaries 1, 2,...,L+1 at the Cartesian coordinates (x,y) in the sample plane; u1 is the first component of u. The symbol (u) denotes the Fresnel reflection coefficient of the corresponding system of smooth films with thicknesses hj replaced with hj + ujuj+1 (putting u0uL+1 ≡ 0 to simplify the notation). The thickness hj must be now understood as the mean thickness of the j-th film. The factor v is equal to
$v=−4πλN0,$
where N0 is refractive index of the ambient. The averaging over u, denoted with angle brackets 〈...〉, is carried out within the illuminated area. Under the assumption that the roughness corresponds to a stationary random process (in the wide sense), the sample-plane averaging can be replaced by statistical averaging with the multi-dimensional probability density function of u, which will be denoted ρ(u):
$〈r^(u)exp(ivt1)〉=∫∫⋯∫r^(u)exp(ivt1)ρ(u) du.$
For Gaussian roughness, the probability density of u is
$ρ(u)=1detS1(2π)(L+1)/2exp(−12uTS−1u),$
where subscripts T and −1 denote transposition and matrix inversion, respectively. The matrix S is a symmetrical positive semi-definite matrix with components
$Si,j=〈uiuj〉=σiσjCi,j,$
where
$σj2=〈uj2〉$
is the rms value of height irregularities of the j-th boundary and the cross-correlation coefficients between i-th and j-th boundary are equal to [35]
$Ci,j=〈uiuj〉σiσj.$

The integral (26) can be evaluated for each term in the consolidated series (19) using the well known formula for multi-dimensional Laplace transform of a general Gaussian function

$∫∫⋯∫1detS1(2π)(L+1)/2exp(ifTu)exp(−12uTS−1u) du=exp(−12fTSf).$
(f is an arbitrary vector) and the following series is obtained:
$〈r^(u)exp(ivt1)〉=r^1exp(−v2σ122)+∑p=1P∑d=1min(L,p)∑mexp(2i∑j=1dmjX^j) Q^(m)H^(m).$
For the normal incidence cosϑ̂j = 1 and thus formula (1) reduces to
$X^j=2πλN^jhj.$
The factor Ĥ(m) is given by
$H^(m)=exp(−12∑q=1d+1∑s=1d+1D^qD^sSq,s),$
where
$D^j=4πλ(mjN^j−mj−1N^j−1).$
Note this expression makes use of the relations m0 ≡ 1 and mj = 0 for j > d to include the edge terms 1 and d+1 in a single expression. For a single layer-system the sum (32) coincides with the well-known sum [12].

Formula (32) is applicable to thin film systems exhibiting arbitrary correlations between the Gaussian boundary roughness of individual layers. Several specific types of rough thin film systems have been discussed in the literature, for some of them simplifications can be achieved due to a special form of S [38, 49].

## 4. Numerical example

Alternative methods of evaluation of the SDT for normal reflectance of rough thin film systems exist. One possibility is the utilisation of the approximate formulae [35]. For a small number of layers, one can also try to evaluate the integral in Eq. (24) by a direct integration. Furthermore, it is possible to attempt an explicit tracing of the largest contributions to the result instead of the indiscriminate inclusion of all contributions, whether large or small, in the series.

To illustrate the trade-offs between precision and calculation efficiency of the various approaches, their performance was compared for a simple example of a three-layer thin film system. In specific, the following procedures were compared:

• Direct evaluation of Eq. (24) by means of a numeric quadrature. The change of variables u = Uv, where UUT = S, transforms the probability density function to a spherically symmetrical multi-dimensional Gaussian with unit dispersion in v. The resulting integral was evaluated using a Cartesian product of one-dimensional Hermite–Gauss quadratures [50], which represents a fairly efficient integration method. The reflectance coefficient (u) in the integrand was expressed using the 2 × 2 matrix formalism. Note that using the recursive formulae would lead to practically identical results because these two approaches are very similar from the point of view of computational efficiency.
• Utilisation of consolidated series (32) presented in this work.
• Ray tracing. This means the light paths were constructed recursively. All possible reflections and transmissions were raced, stopping the recursion when the estimated contribution to the reflection coefficient of the entire system dropped below a prescribed threshold. This stopping criterion permits the exclusion of combinations of reflections and transmissions that contribute little to the result, which is not easy using the series. Formula (34) was then used for each generated path leaving the upper boundary to include the effect of roughness. The exclusion of many small terms that contribute little to the result might be expected to lead to a better computation efficiency than the utilisation of the series. On the other hand, the effect of the combinatorial explosion can be much stronger in this approach.
• Application of the approximate formulae [35].
• Disregarding the roughness, i. e. using the formulae for a system of thin films with smooth boundaries.
This first three procedures represent essentially all reasonable approaches to the exact calculation the reflectance of a multilayer system with arbitrarily correlated randomly rough boundaries within the SDT.

The simulated thin film system is schematically depicted in Fig. 4. It consists of three layers in the following order (from top): Si3N4 (150 nm), SiO2 (130 nm) and Si3N4 (100 nm). The substrate is formed by crystalline silicon. Optical constants of all materials were taken from standard tables [5153].

Fig. 4 Schema of the system of three thin films studied in the numerical example.

The boundary roughness of a three-film system is described by ten parameters, namely the ten independent elements of matrix S or, equivalently, by four σj and six Cjk. For the simulated thin film system, they were constructed based on the notion of a gradual built-up of roughness during deposition. It was assumed that the k-th layer adds independent roughness contribution εk(x, y) of the same rms value σ on top of the k + 1-th layer. The heights of irregularities were thus related:

$uk(x,y)=uk+1(x,y)+εk(x,y) for 1≤k≤L,$
where
$〈uL+1(x,y)uL+1(x,y)〉=σL+12, 〈εk(x,y)εj(x,y)〉=δjkσ2, 〈εk(x,y)uL+1(x,y)〉=0.$
The symbols σL+1 and σ denote the rms of heights of the lowermost boundary and the roughness added at each subsequent film boundary, respectively. It can be seen that this results in the following values of S matrix elements:
$Sjk=σL+12+(L+1−max(j,k))σ2.$
Values σL+1 = 10nm and σ = 5nm were used in the calculation.

The calculated spectral dependence of reflectance in the range from 200 to 800 nm is plotted in Fig. 5 for all procedures. Of course, the series, numeric quadrature and ray tracing represent different numerical methods for the calculation of the same result. Therefore, they all permit, in principle, calculation with an arbitrarily high precision (within the SDT). Hence one curve is used to represent all in the figure as they are indistinguishable when calculated to a sufficient precision. The approximate formulae provide reflectance values that are relatively close to the values obtained by precise SDT. Nevertheless, deviations up to several percents are noticeable, namely in the short-wavelength region. As can be also seen from the figure, the formulae for a system of thin films with smooth boundaries would be inadequate for a thin film system similar to the simulated system.

Fig. 5 Calculated spectral dependence of normal reflectance for the three-layer system. The curve labelled ‘series, integration or ray tracing’ represents the result of all three methods when calculated to a sufficient precision.

It was already noted above that the series, numeric quadrature and ray tracing are not fixed approximations (unlike the approximate formulae) and allow a trade-off between the precision and speed of computation. It is, therefore, interesting to compare them also from this point of view. Choosing a specific precision is relatively straightforward for the series because it suffices to terminate the series at path length p for which the entire p-length term falls below a prescribed value.

In the case of Hermite–Gauss quadratures, the precision control is more difficult as it is not obvious how the number of quadrature points translates to precision of the result. An adaptive multi-dimensional quadrature could provide an error estimate. For the purpose of this illustration, a simpler approach was used: The reflectance was evaluated using quadratures with different numbers of points from 3 to 25, the computation time was measured and the precision was determined ex post by finding the maximum difference from a reference result obtained by evaluation of the series to full double precision. For consistency, the computation time and precision of the other methods was measured in the same manner.

The result of this benchmark is shown in Fig. 6. The time plotted in this figure is the mean computation time per one spectrum point for the three-layer system discussed above, measured on author’s Athlon 64 PC. The approximate formulae occupy the lower right corner of the graph, corresponding to short computation time but low precision. The performance of numeric quadrature may be perhaps acceptable in applications requiring only about three significant digits, even for this precision range the consolidated series represents a considerably better choice. A notable property of the series is the slow increase of computation time required for higher numbers of significant digits. This is illustrated by extending the abscissa to precision 10−12 that would be, of course, absurd for characterisation of rough thin film systems in practice. The slow computation time growth is in a stark contrast with the behaviour of ray tracing. Although ray tracing offers a reasonable speed–precision trade-off at lower precisions (10−3), the computation time grows rapidly with increasing precision due to the combinatorial explosion discussed in paragraph 2.6. For the number of significant digits that would be required for instance for numerical differentiation in data fitting, the computation time is already unacceptable. Furthermore, precision control is more difficult for ray tracing since the recursion stopping threshold does not translate directly to precision of the result. Therefore, it can be concluded that of the five approaches tested, the consolidated series preformed best. In fact, it turned out to be the only method capable of calculation of the reflectance to the full double precision in a reasonable time.

Fig. 6 Trade-off between computation speed and precision for individual calculation procedures. The approximate formulae are represented with a single point since their precision is not tunable.

Finally, a peculiar feature of the series should be noted. It becomes more efficient with increasing roughness, whereas the opposite holds for the application of the matrix formalism and integration. This is because the exp(−x2)-like factors in the terms decrease more rapidly with increasing roughness. For integration, however, increased roughness means larger ranges of u and this implies more abscissa points are required to achieve the same precision of the numeric quadrature. This effect is illustrated in Fig. 7, where the computation efficiencies of both methods are compared in dependence of the roughness magnitude. The computations were performed for rough thin film systems that were identical as in the preceding example, except for the overall roughness magnitude. The overall roughness magnitude was varied by scaling all σ values by the same factor. The value of the lowermost roughness σL+1, which is plotted on the abscissa, was varied in the interval [1.25, 15]nm. The mean computation times plotted in the figure correspond to times needed for obtaining the entire reflectance spectra with precision at least 10−6. Again, the precision was determined as the difference to the values calculated using the series to the full double precision. From Fig. 7 it can be seen that although the efficiency of the consolidated series is decreasing with decreasing roughness, the series is more efficient for all roughness values above approximately 2 nm for this thin film system. Of course, this threshold value depends on the structure and optical constants of the thin film system and would differ for other systems. For zero roughness, i. e. a smooth system, the series is not very efficient. However, for such system, one would use the matrix formalism or recursive formulae in practice.

Fig. 7 Dependence of computation efficiencies of the series and direct integration on the magnitude of roughness. The mean times correspond to computations achieving the precision at least 10−6.

## 5. Conclusion

Fresnel reflectance and transmittance coefficients of the system of an arbitrary number of thin films were expressed in the form of an explicit power series in Fresnel coefficients for individual boundaries and phase terms for the individual films. The series was constructed by considering all possible paths of light through the multiple layers and it is ordered by the length of this path. It was, however, constructed as consolidated, i. e. all paths leading to the same powers are represented using a single term in the series and a numeric multiplicity. Notably, the multiplicities were found to be simple products of binomial coefficients and their evaluation does not require any advanced number-theoretic methods. This fundamental result is interesting from a mathematical point of view, however, the series expression is not a mere mathematical curiosity. Since the expression of Fresnel coefficients of a thin film system as a power series is advantageous for the calculation of optical quantities of rough thin film systems using the scalar diffraction theory, the series was utilised to obtain formulae for the normal reflectance of a thin film system with arbitrarily correlated Gaussian rough boundaries. Using a numerical example, it was illustrated for a three-layer system that such approach is viable and can be used for an efficient calculation of optical properties of rough thin film systems using the scalar diffraction theory.

## Acknowledgments

This work was supported by the projects “CEITEC – Central European Institute of Technology” (CZ.1.05/1.1.00/02.0068) and “R&D center for low-cost plasma and nanotechnology surface modifications” (CZ.1.05/2.1.00/03.0086) from European Regional Development Fund.

1. I. Ohlídal and D. Franta, “Ellipsometry of thin film systems,” in Progress in Optics, E. Wolf, ed. (Elsevier, 2000), vol. 41, pp. 181–282. [CrossRef]

2. R. Jacobson, “Light reflection from films of continuously varying refractive index,” in Progress in Optics, E. Wolf, ed. (North-Holland, 1966), Vol. 5, pp. 249–286.

3. Z. Knittl, Optics of Thin Films (John Wiley, 1976).

4. P. Rouard, “Etudes des propriétés optiques des lames métalliques très minces,” Ann. Phys. 7, 291–384 (1937).

5. A. Vašíček, “The reflection of light from glass with double and multiple films,” J. Opt. Soc. Am. 37, 623–634 (1947). [CrossRef]

6. A. W. Crook, “The reflection and transmission of light by any system of parallel isotropic films,” J. Opt. Soc. Am. 38, 954–963 (1948). [CrossRef]   [PubMed]

7. D. E. Aspnes, J. B. Theeten, and F. Hottier, “Investigation of effective-medium models of microscopic surface roughness by spectroscopic ellipsometry,” Phys. Rev. B 20, 3292–3302 (1979). [CrossRef]

8. J. Szczyrbowski, K. Schmalzbauer, and H. Hoffmann, “Optical properties of rough thin films,” Thin Solid Films 130, 57–73 (1985). [CrossRef]

9. L. Névot, B. Pardo, and J. Corno, “Characterization of X-UV multilayers by grazing incidence X-ray reflectometry,” Rev. Phys. Appl. 23, 1675–1686 (1988). [CrossRef]

10. H. E. Bennett, “Specular reflectance of aluminized ground glass and the height distribution of surface irregularities,” J. Opt. Soc. Am. 53, 1389–1394 (1963). [CrossRef]

11. I. Ohlídal, K. Navrátil, and F. Lukeš, “Reflection of light on a system of non-absorbing isotropic film–non-absorbing isotropic substrate with rough boundaries,” Opt. Commun. 3, 40–44 (1971). [CrossRef]

12. I. Ohlídal, K. Navrátil, and F. Lukeš, “Reflection of light by a system of nonabsorbing isotropic film–nonabsorbing isotropic substrate with randomly rough boundaries,” J. Opt. Soc. Am. 61, 1630–1639 (1971). [CrossRef]

13. J. M. Eastman and P. W. Baumeist, “Measurement of microtopography of optical surfaces using a scanning Fizeau interferometer,” J. Opt. Soc. Am. 64, 1369 (1974).

14. J. M. Eastman, “Scattering in all-dielectric multilayer bandpass filters and mirrors for lasers,” in Physics of Thin Films, G. Hass and H. M. Francombe, eds. (Academic, 1978), Vol. 10, p. 167.

15. C. K. Carniglia, “Scalar scattering theory for multilayer optical coatings,” Opt. Eng. 18, 104–115 (1979). [CrossRef]

16. J. Bauer, L. Biste, and D. Bolze, “Optical-properties of aluminum nitride prepared by chemical and plasmachemical vapor-deposition,” Phys. Status Solidi A 39, 173–181 (1977). [CrossRef]

17. J. M. Zavislan, “Angular scattering from optical interference coatings: scalar scattering predictions and measurements,” Appl. Opt. 30, 2224–2244 (1991). [CrossRef]   [PubMed]

18. C. Amra, J. H. Apfel, and E. Pelltier, “Role of interface correlation in light scattering by a multilayer,” Appl. Opt. 31, 3134–3151 (1992). [CrossRef]   [PubMed]

19. C. Amra, C. Grèzes-Besset, and L. Bruel, “Comparison of surface and bulk scattering in optical multilayers,” Appl. Opt. 32, 5492–5503 (1993). [CrossRef]   [PubMed]

20. C. Amra, “Light scattering from multilayer optics I. tools of investigation,” J. Opt. Soc. Am. A 11, 197–210 (1994). [CrossRef]

21. A. Duparré, J. Ferre-Borrull, S. Gliech, G. Notni, J. Steinert, and J. M. Bennett, “Surface characterization techniques for determining the root-mean-square roughness and power spectral densities of optical components,” Appl. Opt. 41, 154–171 (2002). [CrossRef]   [PubMed]

22. J. Meunier, “Optical reflectivity of thin rough films: Application to ellipsometric measurements of liquid films,” Phys. Rev. E 75, 061601 (2007). [CrossRef]

23. A. A. Maradudin, ed., Light Scattering and Nanoscale Surface Roughness (Springer, 2010).

24. S. Schröder, T. Herffurth, H. Blaschke, and A. Duparré, “Angle-resolved scattering: an effective method for characterizing thin-film coatings,” Appl. Opt. 50, C164–C171 (2011). [CrossRef]   [PubMed]

25. J. E. Harvey, S. Schröder, N. Choi, and A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002). [CrossRef]

26. T. Herffurth, S. Schröder, M. Trost, and A. D. A. Tünnermann, “Comprehensive nanostructure and defect analysis using a simple 3D light-scatter sensor,” Appl. Opt. 52, 3279–3287 (2013). [CrossRef]   [PubMed]

27. J. W. S. Rayleigh, Theory of Sound (Macmillan, 1877), Vol. 2.

28. S. O. Rice, “Reflection of electromagnetic waves from slightly rough surfaces,” Commun. Pure Appl. Math. 4, 351–378 (1951). [CrossRef]

29. G. R. Valenzuela, “Depolarization of EM waves by slightly rough surfaces,” IEEE Trans. Antennas Propag. 15, 552–557 (1967). [CrossRef]

30. R. Schiffer, “Reflectivity of a slightly rough surface,” Appl. Opt. 26, 704–712 (1987). [CrossRef]   [PubMed]

31. V. Freilikher, M. Pustilnik, I. Yurkevich, and V. I. Tatarskii, “Polarization of light scattered from slightly rough dielectric film,” Opt. Lett. 19, 1382–1384 (1994). [CrossRef]   [PubMed]

32. D. Franta and I. Ohlídal, “Ellipsometric parameters and reflectances of thin films with slightly rough boundaries,” J. Mod. Opt. 45, 903–934 (1998). [CrossRef]

33. A. Krywonos, “Predicting surface scatter using a linear systems formulation of non-paraxial scalar diffraction,” Ph.D. thesis, University of Central Florida, Orlando (2006).

34. 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. Eng. 19, 9820–9835 (2011).

35. I. Ohlídal and F. Vižd’a, “Optical quantities of multilayer systems with correlated randomly rough boundaries,” J. Mod. Opt. 46, 2043–2062 (1999). [CrossRef]

36. M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, and H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005). [CrossRef]

37. I. Ohlídal and F. Lukeš, “Ellipsometric parameters of rough surfaces and of a system substrate-thin film with rough boundaries,” Opt. Acta 19, 817–843 (1972). [CrossRef]

38. I. Ohlídal, F. Lukeš, and K. Navrátil, “Rough silicon surfaces studied by optical methods,” Surf. Sci. 45, 91–116 (1974). [CrossRef]

39. J. T. Butler, “On the number of propagation paths in multilayer media,” Fibonacci Q. 28, 334–339 (1990).

40. L. D. Killough, “The on-line encyclopedia of integer sequences,” http://oeis.org/A011782. Sequence A011782.

41. N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A001519. Sequence A001519.

42. C. Mallows, N. J. A. Sloane, S. Plouffe, and R. G. Wilson, “The on-line encyclopedia of integer sequences,” http://oeis.org/A007051. Sequence A007051.

43. H. Bottomley, “The on-line encyclopedia of integer sequences,” http://oeis.org/A080937. Sequence A080937.

44. N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A024175. Sequence A024175.

45. N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A000108. Sequence A000108.

46. M. Kildemo, O. Hunderi, and B. Drévillon, “Approximation of reflection coefficients for rapid real-time calculation of inhomogeneous films,” J. Opt. Soc. Am. A 14, 931–939 (1997). [CrossRef]

47. I. Ohlídal, “Reflectance of multilayer systems with randomly rough boundaries,” Opt. Commun. 71, 323–326 (1989). [CrossRef]

48. I. Ohlídal, “Approximate formulas for the reflectances, transmittances, and scattering losses of nonabsorbing multilayers systems with randomly rough boundaries,” J. Opt. Soc. Am. A 10, 158–170 (1993). [CrossRef]

49. I. Ohlídal, K. Navrátil, and M. Ohlídal, “Scattering of light from multilayer systems with rough boundaries,” in Progress in Optics, E. Wolf, ed. (Elsevier, 1995), Vol. 34, pp. 249–331. [CrossRef]

50. A. H. Stroud and D. Secrest, Gaussian Quadrature Formulas (Prentice-Hall, 1966).

51. H. R. Philipp, “Silicon dioxide (SiO2) (glass),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 749–763. [CrossRef]

52. H. R. Philipp, “Silicon nitride (Si3N4) (noncrystalline),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 771–774. [CrossRef]

53. C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, and W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998). [CrossRef]

### References

• View by:
• |
• |
• |

1. I. Ohlídal, D. Franta, “Ellipsometry of thin film systems,” in Progress in Optics, E. Wolf, ed. (Elsevier, 2000), vol. 41, pp. 181–282.
[CrossRef]
2. R. Jacobson, “Light reflection from films of continuously varying refractive index,” in Progress in Optics, E. Wolf, ed. (North-Holland, 1966), Vol. 5, pp. 249–286.
3. Z. Knittl, Optics of Thin Films (John Wiley, 1976).
4. P. Rouard, “Etudes des propriétés optiques des lames métalliques très minces,” Ann. Phys. 7, 291–384 (1937).
5. A. Vašíček, “The reflection of light from glass with double and multiple films,” J. Opt. Soc. Am. 37, 623–634 (1947).
[CrossRef]
6. A. W. Crook, “The reflection and transmission of light by any system of parallel isotropic films,” J. Opt. Soc. Am. 38, 954–963 (1948).
[CrossRef] [PubMed]
7. D. E. Aspnes, J. B. Theeten, F. Hottier, “Investigation of effective-medium models of microscopic surface roughness by spectroscopic ellipsometry,” Phys. Rev. B 20, 3292–3302 (1979).
[CrossRef]
8. J. Szczyrbowski, K. Schmalzbauer, H. Hoffmann, “Optical properties of rough thin films,” Thin Solid Films 130, 57–73 (1985).
[CrossRef]
9. L. Névot, B. Pardo, J. Corno, “Characterization of X-UV multilayers by grazing incidence X-ray reflectometry,” Rev. Phys. Appl. 23, 1675–1686 (1988).
[CrossRef]
10. H. E. Bennett, “Specular reflectance of aluminized ground glass and the height distribution of surface irregularities,” J. Opt. Soc. Am. 53, 1389–1394 (1963).
[CrossRef]
11. I. Ohlídal, K. Navrátil, F. Lukeš, “Reflection of light on a system of non-absorbing isotropic film–non-absorbing isotropic substrate with rough boundaries,” Opt. Commun. 3, 40–44 (1971).
[CrossRef]
12. I. Ohlídal, K. Navrátil, F. Lukeš, “Reflection of light by a system of nonabsorbing isotropic film–nonabsorbing isotropic substrate with randomly rough boundaries,” J. Opt. Soc. Am. 61, 1630–1639 (1971).
[CrossRef]
13. J. M. Eastman, P. W. Baumeist, “Measurement of microtopography of optical surfaces using a scanning Fizeau interferometer,” J. Opt. Soc. Am. 64, 1369 (1974).
14. J. M. Eastman, “Scattering in all-dielectric multilayer bandpass filters and mirrors for lasers,” in Physics of Thin Films, G. Hass, H. M. Francombe, eds. (Academic, 1978), Vol. 10, p. 167.
15. C. K. Carniglia, “Scalar scattering theory for multilayer optical coatings,” Opt. Eng. 18, 104–115 (1979).
[CrossRef]
16. J. Bauer, L. Biste, D. Bolze, “Optical-properties of aluminum nitride prepared by chemical and plasmachemical vapor-deposition,” Phys. Status Solidi A 39, 173–181 (1977).
[CrossRef]
17. J. M. Zavislan, “Angular scattering from optical interference coatings: scalar scattering predictions and measurements,” Appl. Opt. 30, 2224–2244 (1991).
[CrossRef] [PubMed]
18. C. Amra, J. H. Apfel, E. Pelltier, “Role of interface correlation in light scattering by a multilayer,” Appl. Opt. 31, 3134–3151 (1992).
[CrossRef] [PubMed]
19. C. Amra, C. Grèzes-Besset, L. Bruel, “Comparison of surface and bulk scattering in optical multilayers,” Appl. Opt. 32, 5492–5503 (1993).
[CrossRef] [PubMed]
20. C. Amra, “Light scattering from multilayer optics I. tools of investigation,” J. Opt. Soc. Am. A 11, 197–210 (1994).
[CrossRef]
21. A. Duparré, J. Ferre-Borrull, S. Gliech, G. Notni, J. Steinert, J. M. Bennett, “Surface characterization techniques for determining the root-mean-square roughness and power spectral densities of optical components,” Appl. Opt. 41, 154–171 (2002).
[CrossRef] [PubMed]
22. J. Meunier, “Optical reflectivity of thin rough films: Application to ellipsometric measurements of liquid films,” Phys. Rev. E 75, 061601 (2007).
[CrossRef]
23. A. A. Maradudin, ed., Light Scattering and Nanoscale Surface Roughness (Springer, 2010).
24. S. Schröder, T. Herffurth, H. Blaschke, A. Duparré, “Angle-resolved scattering: an effective method for characterizing thin-film coatings,” Appl. Opt. 50, C164–C171 (2011).
[CrossRef] [PubMed]
25. J. E. Harvey, S. Schröder, N. Choi, A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002).
[CrossRef]
26. T. Herffurth, S. Schröder, M. Trost, A. D. A. Tünnermann, “Comprehensive nanostructure and defect analysis using a simple 3D light-scatter sensor,” Appl. Opt. 52, 3279–3287 (2013).
[CrossRef] [PubMed]
27. J. W. S. Rayleigh, Theory of Sound (Macmillan, 1877), Vol. 2.
28. S. O. Rice, “Reflection of electromagnetic waves from slightly rough surfaces,” Commun. Pure Appl. Math. 4, 351–378 (1951).
[CrossRef]
29. G. R. Valenzuela, “Depolarization of EM waves by slightly rough surfaces,” IEEE Trans. Antennas Propag. 15, 552–557 (1967).
[CrossRef]
30. R. Schiffer, “Reflectivity of a slightly rough surface,” Appl. Opt. 26, 704–712 (1987).
[CrossRef] [PubMed]
31. V. Freilikher, M. Pustilnik, I. Yurkevich, V. I. Tatarskii, “Polarization of light scattered from slightly rough dielectric film,” Opt. Lett. 19, 1382–1384 (1994).
[CrossRef] [PubMed]
32. D. Franta, I. Ohlídal, “Ellipsometric parameters and reflectances of thin films with slightly rough boundaries,” J. Mod. Opt. 45, 903–934 (1998).
[CrossRef]
33. A. Krywonos, “Predicting surface scatter using a linear systems formulation of non-paraxial scalar diffraction,” Ph.D. thesis, University of Central Florida, Orlando (2006).
34. S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).
35. I. Ohlídal, F. Vižd’a, “Optical quantities of multilayer systems with correlated randomly rough boundaries,” J. Mod. Opt. 46, 2043–2062 (1999).
[CrossRef]
36. M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]
37. I. Ohlídal, F. Lukeš, “Ellipsometric parameters of rough surfaces and of a system substrate-thin film with rough boundaries,” Opt. Acta 19, 817–843 (1972).
[CrossRef]
38. I. Ohlídal, F. Lukeš, K. Navrátil, “Rough silicon surfaces studied by optical methods,” Surf. Sci. 45, 91–116 (1974).
[CrossRef]
39. J. T. Butler, “On the number of propagation paths in multilayer media,” Fibonacci Q. 28, 334–339 (1990).
40. L. D. Killough, “The on-line encyclopedia of integer sequences,” http://oeis.org/A011782 . Sequence A011782.
41. N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A001519 . Sequence A001519.
42. C. Mallows, N. J. A. Sloane, S. Plouffe, R. G. Wilson, “The on-line encyclopedia of integer sequences,” http://oeis.org/A007051 . Sequence A007051.
43. H. Bottomley, “The on-line encyclopedia of integer sequences,” http://oeis.org/A080937 . Sequence A080937.
44. N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A024175 . Sequence A024175.
45. N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A000108 . Sequence A000108.
46. M. Kildemo, O. Hunderi, B. Drévillon, “Approximation of reflection coefficients for rapid real-time calculation of inhomogeneous films,” J. Opt. Soc. Am. A 14, 931–939 (1997).
[CrossRef]
47. I. Ohlídal, “Reflectance of multilayer systems with randomly rough boundaries,” Opt. Commun. 71, 323–326 (1989).
[CrossRef]
48. I. Ohlídal, “Approximate formulas for the reflectances, transmittances, and scattering losses of nonabsorbing multilayers systems with randomly rough boundaries,” J. Opt. Soc. Am. A 10, 158–170 (1993).
[CrossRef]
49. I. Ohlídal, K. Navrátil, M. Ohlídal, “Scattering of light from multilayer systems with rough boundaries,” in Progress in Optics, E. Wolf, ed. (Elsevier, 1995), Vol. 34, pp. 249–331.
[CrossRef]
50. A. H. Stroud, D. Secrest, Gaussian Quadrature Formulas (Prentice-Hall, 1966).
51. H. R. Philipp, “Silicon dioxide (SiO2) (glass),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 749–763.
[CrossRef]
52. H. R. Philipp, “Silicon nitride (Si3N4) (noncrystalline),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 771–774.
[CrossRef]
53. C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### 2011 (2)

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

#### 2007 (1)

J. Meunier, “Optical reflectivity of thin rough films: Application to ellipsometric measurements of liquid films,” Phys. Rev. E 75, 061601 (2007).
[CrossRef]

#### 2005 (1)

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

#### 2002 (2)

J. E. Harvey, S. Schröder, N. Choi, A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002).
[CrossRef]

#### 1999 (1)

I. Ohlídal, F. Vižd’a, “Optical quantities of multilayer systems with correlated randomly rough boundaries,” J. Mod. Opt. 46, 2043–2062 (1999).
[CrossRef]

#### 1998 (2)

D. Franta, I. Ohlídal, “Ellipsometric parameters and reflectances of thin films with slightly rough boundaries,” J. Mod. Opt. 45, 903–934 (1998).
[CrossRef]

C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### 1990 (1)

J. T. Butler, “On the number of propagation paths in multilayer media,” Fibonacci Q. 28, 334–339 (1990).

#### 1989 (1)

I. Ohlídal, “Reflectance of multilayer systems with randomly rough boundaries,” Opt. Commun. 71, 323–326 (1989).
[CrossRef]

#### 1988 (1)

L. Névot, B. Pardo, J. Corno, “Characterization of X-UV multilayers by grazing incidence X-ray reflectometry,” Rev. Phys. Appl. 23, 1675–1686 (1988).
[CrossRef]

#### 1985 (1)

J. Szczyrbowski, K. Schmalzbauer, H. Hoffmann, “Optical properties of rough thin films,” Thin Solid Films 130, 57–73 (1985).
[CrossRef]

#### 1979 (2)

D. E. Aspnes, J. B. Theeten, F. Hottier, “Investigation of effective-medium models of microscopic surface roughness by spectroscopic ellipsometry,” Phys. Rev. B 20, 3292–3302 (1979).
[CrossRef]

C. K. Carniglia, “Scalar scattering theory for multilayer optical coatings,” Opt. Eng. 18, 104–115 (1979).
[CrossRef]

#### 1977 (1)

J. Bauer, L. Biste, D. Bolze, “Optical-properties of aluminum nitride prepared by chemical and plasmachemical vapor-deposition,” Phys. Status Solidi A 39, 173–181 (1977).
[CrossRef]

#### 1974 (2)

J. M. Eastman, P. W. Baumeist, “Measurement of microtopography of optical surfaces using a scanning Fizeau interferometer,” J. Opt. Soc. Am. 64, 1369 (1974).

I. Ohlídal, F. Lukeš, K. Navrátil, “Rough silicon surfaces studied by optical methods,” Surf. Sci. 45, 91–116 (1974).
[CrossRef]

#### 1972 (1)

I. Ohlídal, F. Lukeš, “Ellipsometric parameters of rough surfaces and of a system substrate-thin film with rough boundaries,” Opt. Acta 19, 817–843 (1972).
[CrossRef]

#### 1971 (2)

I. Ohlídal, K. Navrátil, F. Lukeš, “Reflection of light on a system of non-absorbing isotropic film–non-absorbing isotropic substrate with rough boundaries,” Opt. Commun. 3, 40–44 (1971).
[CrossRef]

#### 1967 (1)

G. R. Valenzuela, “Depolarization of EM waves by slightly rough surfaces,” IEEE Trans. Antennas Propag. 15, 552–557 (1967).
[CrossRef]

#### 1951 (1)

S. O. Rice, “Reflection of electromagnetic waves from slightly rough surfaces,” Commun. Pure Appl. Math. 4, 351–378 (1951).
[CrossRef]

#### 1937 (1)

P. Rouard, “Etudes des propriétés optiques des lames métalliques très minces,” Ann. Phys. 7, 291–384 (1937).

#### Aspnes, D. E.

D. E. Aspnes, J. B. Theeten, F. Hottier, “Investigation of effective-medium models of microscopic surface roughness by spectroscopic ellipsometry,” Phys. Rev. B 20, 3292–3302 (1979).
[CrossRef]

#### Bauer, J.

J. Bauer, L. Biste, D. Bolze, “Optical-properties of aluminum nitride prepared by chemical and plasmachemical vapor-deposition,” Phys. Status Solidi A 39, 173–181 (1977).
[CrossRef]

#### Baumeist, P. W.

J. M. Eastman, P. W. Baumeist, “Measurement of microtopography of optical surfaces using a scanning Fizeau interferometer,” J. Opt. Soc. Am. 64, 1369 (1974).

#### Biste, L.

J. Bauer, L. Biste, D. Bolze, “Optical-properties of aluminum nitride prepared by chemical and plasmachemical vapor-deposition,” Phys. Status Solidi A 39, 173–181 (1977).
[CrossRef]

#### Bolze, D.

J. Bauer, L. Biste, D. Bolze, “Optical-properties of aluminum nitride prepared by chemical and plasmachemical vapor-deposition,” Phys. Status Solidi A 39, 173–181 (1977).
[CrossRef]

#### Bonanni, A.

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

#### Bottomley, H.

H. Bottomley, “The on-line encyclopedia of integer sequences,” http://oeis.org/A080937 . Sequence A080937.

#### Butler, J. T.

J. T. Butler, “On the number of propagation paths in multilayer media,” Fibonacci Q. 28, 334–339 (1990).

#### Carniglia, C. K.

C. K. Carniglia, “Scalar scattering theory for multilayer optical coatings,” Opt. Eng. 18, 104–115 (1979).
[CrossRef]

#### Choi, N.

J. E. Harvey, S. Schröder, N. Choi, A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002).
[CrossRef]

#### Coriand, L.

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

#### Corno, J.

L. Névot, B. Pardo, J. Corno, “Characterization of X-UV multilayers by grazing incidence X-ray reflectometry,” Rev. Phys. Appl. 23, 1675–1686 (1988).
[CrossRef]

#### Duparré, A.

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

J. E. Harvey, S. Schröder, N. Choi, A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002).
[CrossRef]

#### Eastman, J. M.

J. M. Eastman, P. W. Baumeist, “Measurement of microtopography of optical surfaces using a scanning Fizeau interferometer,” J. Opt. Soc. Am. 64, 1369 (1974).

J. M. Eastman, “Scattering in all-dielectric multilayer bandpass filters and mirrors for lasers,” in Physics of Thin Films, G. Hass, H. M. Francombe, eds. (Academic, 1978), Vol. 10, p. 167.

#### Franta, D.

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

D. Franta, I. Ohlídal, “Ellipsometric parameters and reflectances of thin films with slightly rough boundaries,” J. Mod. Opt. 45, 903–934 (1998).
[CrossRef]

I. Ohlídal, D. Franta, “Ellipsometry of thin film systems,” in Progress in Optics, E. Wolf, ed. (Elsevier, 2000), vol. 41, pp. 181–282.
[CrossRef]

#### Harvey, J. E.

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

J. E. Harvey, S. Schröder, N. Choi, A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002).
[CrossRef]

#### Herzinger, C. M.

C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### Hoffmann, H.

J. Szczyrbowski, K. Schmalzbauer, H. Hoffmann, “Optical properties of rough thin films,” Thin Solid Films 130, 57–73 (1985).
[CrossRef]

#### Hottier, F.

D. E. Aspnes, J. B. Theeten, F. Hottier, “Investigation of effective-medium models of microscopic surface roughness by spectroscopic ellipsometry,” Phys. Rev. B 20, 3292–3302 (1979).
[CrossRef]

#### Jacobson, R.

R. Jacobson, “Light reflection from films of continuously varying refractive index,” in Progress in Optics, E. Wolf, ed. (North-Holland, 1966), Vol. 5, pp. 249–286.

#### Johs, B.

C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### Killough, L. D.

L. D. Killough, “The on-line encyclopedia of integer sequences,” http://oeis.org/A011782 . Sequence A011782.

#### Knittl, Z.

Z. Knittl, Optics of Thin Films (John Wiley, 1976).

#### Krywonos, A.

A. Krywonos, “Predicting surface scatter using a linear systems formulation of non-paraxial scalar diffraction,” Ph.D. thesis, University of Central Florida, Orlando (2006).

#### Lukeš, F.

I. Ohlídal, F. Lukeš, K. Navrátil, “Rough silicon surfaces studied by optical methods,” Surf. Sci. 45, 91–116 (1974).
[CrossRef]

I. Ohlídal, F. Lukeš, “Ellipsometric parameters of rough surfaces and of a system substrate-thin film with rough boundaries,” Opt. Acta 19, 817–843 (1972).
[CrossRef]

I. Ohlídal, K. Navrátil, F. Lukeš, “Reflection of light on a system of non-absorbing isotropic film–non-absorbing isotropic substrate with rough boundaries,” Opt. Commun. 3, 40–44 (1971).
[CrossRef]

#### Mallows, C.

C. Mallows, N. J. A. Sloane, S. Plouffe, R. G. Wilson, “The on-line encyclopedia of integer sequences,” http://oeis.org/A007051 . Sequence A007051.

#### McGahan, W. A.

C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### Meunier, J.

J. Meunier, “Optical reflectivity of thin rough films: Application to ellipsometric measurements of liquid films,” Phys. Rev. E 75, 061601 (2007).
[CrossRef]

#### Montaigne-Ramil, A.

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

#### Navrátil, K.

I. Ohlídal, F. Lukeš, K. Navrátil, “Rough silicon surfaces studied by optical methods,” Surf. Sci. 45, 91–116 (1974).
[CrossRef]

I. Ohlídal, K. Navrátil, F. Lukeš, “Reflection of light on a system of non-absorbing isotropic film–non-absorbing isotropic substrate with rough boundaries,” Opt. Commun. 3, 40–44 (1971).
[CrossRef]

I. Ohlídal, K. Navrátil, M. Ohlídal, “Scattering of light from multilayer systems with rough boundaries,” in Progress in Optics, E. Wolf, ed. (Elsevier, 1995), Vol. 34, pp. 249–331.
[CrossRef]

#### Névot, L.

L. Névot, B. Pardo, J. Corno, “Characterization of X-UV multilayers by grazing incidence X-ray reflectometry,” Rev. Phys. Appl. 23, 1675–1686 (1988).
[CrossRef]

#### Ohlídal, I.

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

I. Ohlídal, F. Vižd’a, “Optical quantities of multilayer systems with correlated randomly rough boundaries,” J. Mod. Opt. 46, 2043–2062 (1999).
[CrossRef]

D. Franta, I. Ohlídal, “Ellipsometric parameters and reflectances of thin films with slightly rough boundaries,” J. Mod. Opt. 45, 903–934 (1998).
[CrossRef]

I. Ohlídal, “Reflectance of multilayer systems with randomly rough boundaries,” Opt. Commun. 71, 323–326 (1989).
[CrossRef]

I. Ohlídal, F. Lukeš, K. Navrátil, “Rough silicon surfaces studied by optical methods,” Surf. Sci. 45, 91–116 (1974).
[CrossRef]

I. Ohlídal, F. Lukeš, “Ellipsometric parameters of rough surfaces and of a system substrate-thin film with rough boundaries,” Opt. Acta 19, 817–843 (1972).
[CrossRef]

I. Ohlídal, K. Navrátil, F. Lukeš, “Reflection of light on a system of non-absorbing isotropic film–non-absorbing isotropic substrate with rough boundaries,” Opt. Commun. 3, 40–44 (1971).
[CrossRef]

I. Ohlídal, K. Navrátil, M. Ohlídal, “Scattering of light from multilayer systems with rough boundaries,” in Progress in Optics, E. Wolf, ed. (Elsevier, 1995), Vol. 34, pp. 249–331.
[CrossRef]

I. Ohlídal, D. Franta, “Ellipsometry of thin film systems,” in Progress in Optics, E. Wolf, ed. (Elsevier, 2000), vol. 41, pp. 181–282.
[CrossRef]

#### Ohlídal, M.

I. Ohlídal, K. Navrátil, M. Ohlídal, “Scattering of light from multilayer systems with rough boundaries,” in Progress in Optics, E. Wolf, ed. (Elsevier, 1995), Vol. 34, pp. 249–331.
[CrossRef]

#### Pardo, B.

L. Névot, B. Pardo, J. Corno, “Characterization of X-UV multilayers by grazing incidence X-ray reflectometry,” Rev. Phys. Appl. 23, 1675–1686 (1988).
[CrossRef]

#### Paulson, W.

C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### Penalver, D. H.

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

#### Philipp, H. R.

H. R. Philipp, “Silicon nitride (Si3N4) (noncrystalline),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 771–774.
[CrossRef]

H. R. Philipp, “Silicon dioxide (SiO2) (glass),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 749–763.
[CrossRef]

#### Plouffe, S.

C. Mallows, N. J. A. Sloane, S. Plouffe, R. G. Wilson, “The on-line encyclopedia of integer sequences,” http://oeis.org/A007051 . Sequence A007051.

#### Rayleigh, J. W. S.

J. W. S. Rayleigh, Theory of Sound (Macmillan, 1877), Vol. 2.

#### Rice, S. O.

S. O. Rice, “Reflection of electromagnetic waves from slightly rough surfaces,” Commun. Pure Appl. Math. 4, 351–378 (1951).
[CrossRef]

#### Rouard, P.

P. Rouard, “Etudes des propriétés optiques des lames métalliques très minces,” Ann. Phys. 7, 291–384 (1937).

#### Schmalzbauer, K.

J. Szczyrbowski, K. Schmalzbauer, H. Hoffmann, “Optical properties of rough thin films,” Thin Solid Films 130, 57–73 (1985).
[CrossRef]

#### Schröder, S.

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

J. E. Harvey, S. Schröder, N. Choi, A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002).
[CrossRef]

#### Secrest, D.

A. H. Stroud, D. Secrest, Gaussian Quadrature Formulas (Prentice-Hall, 1966).

#### Šiler, M.

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

#### Sitter, H.

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

#### Sloane, N. J. A.

N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A024175 . Sequence A024175.

C. Mallows, N. J. A. Sloane, S. Plouffe, R. G. Wilson, “The on-line encyclopedia of integer sequences,” http://oeis.org/A007051 . Sequence A007051.

N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A000108 . Sequence A000108.

N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A001519 . Sequence A001519.

#### Stifter, D.

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

#### Stroud, A. H.

A. H. Stroud, D. Secrest, Gaussian Quadrature Formulas (Prentice-Hall, 1966).

#### Szczyrbowski, J.

J. Szczyrbowski, K. Schmalzbauer, H. Hoffmann, “Optical properties of rough thin films,” Thin Solid Films 130, 57–73 (1985).
[CrossRef]

#### Theeten, J. B.

D. E. Aspnes, J. B. Theeten, F. Hottier, “Investigation of effective-medium models of microscopic surface roughness by spectroscopic ellipsometry,” Phys. Rev. B 20, 3292–3302 (1979).
[CrossRef]

#### Tünnermann, A.

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

#### Valenzuela, G. R.

G. R. Valenzuela, “Depolarization of EM waves by slightly rough surfaces,” IEEE Trans. Antennas Propag. 15, 552–557 (1967).
[CrossRef]

#### Vižd’a, F.

I. Ohlídal, F. Vižd’a, “Optical quantities of multilayer systems with correlated randomly rough boundaries,” J. Mod. Opt. 46, 2043–2062 (1999).
[CrossRef]

#### Wilson, R. G.

C. Mallows, N. J. A. Sloane, S. Plouffe, R. G. Wilson, “The on-line encyclopedia of integer sequences,” http://oeis.org/A007051 . Sequence A007051.

#### Woollam, J. A.

C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### Ann. Phys. (1)

P. Rouard, “Etudes des propriétés optiques des lames métalliques très minces,” Ann. Phys. 7, 291–384 (1937).

#### Commun. Pure Appl. Math. (1)

S. O. Rice, “Reflection of electromagnetic waves from slightly rough surfaces,” Commun. Pure Appl. Math. 4, 351–378 (1951).
[CrossRef]

#### Fibonacci Q. (1)

J. T. Butler, “On the number of propagation paths in multilayer media,” Fibonacci Q. 28, 334–339 (1990).

#### IEEE Trans. Antennas Propag. (1)

G. R. Valenzuela, “Depolarization of EM waves by slightly rough surfaces,” IEEE Trans. Antennas Propag. 15, 552–557 (1967).
[CrossRef]

#### J. Appl. Phys. (1)

C. M. Herzinger, B. Johs, W. A. McGahan, J. A. Woollam, W. Paulson, “Ellipsometric determination of optical constants for silicon and thermally grown silicon dioxide via a multi-sample, multi-wavelength, multi-angle investigation,” J. Appl. Phys. 83, 3323–3336 (1998).
[CrossRef]

#### J. Mod. Opt. (3)

I. Ohlídal, F. Vižd’a, “Optical quantities of multilayer systems with correlated randomly rough boundaries,” J. Mod. Opt. 46, 2043–2062 (1999).
[CrossRef]

M. Šiler, I. Ohlídal, D. Franta, A. Montaigne-Ramil, A. Bonanni, D. Stifter, H. Sitter, “Optical characterization of double layers containing epitaxial ZnSe and ZnTe films,” J. Mod. Opt. 52, 583–602 (2005).
[CrossRef]

D. Franta, I. Ohlídal, “Ellipsometric parameters and reflectances of thin films with slightly rough boundaries,” J. Mod. Opt. 45, 903–934 (1998).
[CrossRef]

#### J. Opt. Soc. Am. (5)

J. M. Eastman, P. W. Baumeist, “Measurement of microtopography of optical surfaces using a scanning Fizeau interferometer,” J. Opt. Soc. Am. 64, 1369 (1974).

#### Opt. Acta (1)

I. Ohlídal, F. Lukeš, “Ellipsometric parameters of rough surfaces and of a system substrate-thin film with rough boundaries,” Opt. Acta 19, 817–843 (1972).
[CrossRef]

#### Opt. Commun. (2)

I. Ohlídal, “Reflectance of multilayer systems with randomly rough boundaries,” Opt. Commun. 71, 323–326 (1989).
[CrossRef]

I. Ohlídal, K. Navrátil, F. Lukeš, “Reflection of light on a system of non-absorbing isotropic film–non-absorbing isotropic substrate with rough boundaries,” Opt. Commun. 3, 40–44 (1971).
[CrossRef]

#### Opt. Eng. (3)

J. E. Harvey, S. Schröder, N. Choi, A. Duparré, “Total integrated scatter from surfaces with arbitrary roughness, correlation widths, and incident angles,” Opt. Eng. 51, 013402 (2002).
[CrossRef]

C. K. Carniglia, “Scalar scattering theory for multilayer optical coatings,” Opt. Eng. 18, 104–115 (1979).
[CrossRef]

S. Schröder, A. Duparré, L. Coriand, A. Tünnermann, D. H. Penalver, J. E. Harvey, “Modeling of light scattering in different regimes of surface roughness,” Opt. Eng. 19, 9820–9835 (2011).

#### Phys. Rev. B (1)

D. E. Aspnes, J. B. Theeten, F. Hottier, “Investigation of effective-medium models of microscopic surface roughness by spectroscopic ellipsometry,” Phys. Rev. B 20, 3292–3302 (1979).
[CrossRef]

#### Phys. Rev. E (1)

J. Meunier, “Optical reflectivity of thin rough films: Application to ellipsometric measurements of liquid films,” Phys. Rev. E 75, 061601 (2007).
[CrossRef]

#### Phys. Status Solidi A (1)

J. Bauer, L. Biste, D. Bolze, “Optical-properties of aluminum nitride prepared by chemical and plasmachemical vapor-deposition,” Phys. Status Solidi A 39, 173–181 (1977).
[CrossRef]

#### Rev. Phys. Appl. (1)

L. Névot, B. Pardo, J. Corno, “Characterization of X-UV multilayers by grazing incidence X-ray reflectometry,” Rev. Phys. Appl. 23, 1675–1686 (1988).
[CrossRef]

#### Surf. Sci. (1)

I. Ohlídal, F. Lukeš, K. Navrátil, “Rough silicon surfaces studied by optical methods,” Surf. Sci. 45, 91–116 (1974).
[CrossRef]

#### Thin Solid Films (1)

J. Szczyrbowski, K. Schmalzbauer, H. Hoffmann, “Optical properties of rough thin films,” Thin Solid Films 130, 57–73 (1985).
[CrossRef]

#### Other (17)

I. Ohlídal, D. Franta, “Ellipsometry of thin film systems,” in Progress in Optics, E. Wolf, ed. (Elsevier, 2000), vol. 41, pp. 181–282.
[CrossRef]

R. Jacobson, “Light reflection from films of continuously varying refractive index,” in Progress in Optics, E. Wolf, ed. (North-Holland, 1966), Vol. 5, pp. 249–286.

Z. Knittl, Optics of Thin Films (John Wiley, 1976).

J. M. Eastman, “Scattering in all-dielectric multilayer bandpass filters and mirrors for lasers,” in Physics of Thin Films, G. Hass, H. M. Francombe, eds. (Academic, 1978), Vol. 10, p. 167.

A. A. Maradudin, ed., Light Scattering and Nanoscale Surface Roughness (Springer, 2010).

J. W. S. Rayleigh, Theory of Sound (Macmillan, 1877), Vol. 2.

A. Krywonos, “Predicting surface scatter using a linear systems formulation of non-paraxial scalar diffraction,” Ph.D. thesis, University of Central Florida, Orlando (2006).

L. D. Killough, “The on-line encyclopedia of integer sequences,” http://oeis.org/A011782 . Sequence A011782.

N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A001519 . Sequence A001519.

C. Mallows, N. J. A. Sloane, S. Plouffe, R. G. Wilson, “The on-line encyclopedia of integer sequences,” http://oeis.org/A007051 . Sequence A007051.

H. Bottomley, “The on-line encyclopedia of integer sequences,” http://oeis.org/A080937 . Sequence A080937.

N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A024175 . Sequence A024175.

N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” http://oeis.org/A000108 . Sequence A000108.

I. Ohlídal, K. Navrátil, M. Ohlídal, “Scattering of light from multilayer systems with rough boundaries,” in Progress in Optics, E. Wolf, ed. (Elsevier, 1995), Vol. 34, pp. 249–331.
[CrossRef]

A. H. Stroud, D. Secrest, Gaussian Quadrature Formulas (Prentice-Hall, 1966).

H. R. Philipp, “Silicon dioxide (SiO2) (glass),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 749–763.
[CrossRef]

H. R. Philipp, “Silicon nitride (Si3N4) (noncrystalline),” in Handbook of Optical Constants of Solids, E. Palik, ed. (Academic, 1985), Vol. I, pp. 771–774.
[CrossRef]

### Cited By

OSA participates in CrossRef's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

### Figures (7)

Fig. 1

Illustration of equivalent light paths through a two-layer system. Paths (a) and (b) lead to identical powers of both the Fresnel coefficients and phase terms, whereas (c) leads only to the same phase terms powers as (a) and (b); the Fresnel coefficients differ.

Fig. 2

Reflection and transmission coefficient and phase term notation.

Fig. 3

Relations between medium and visit counts connected with medium j and the numbers of reflections within this medium and refractions to or from this medium.

Fig. 4

Schema of the system of three thin films studied in the numerical example.

Fig. 5

Calculated spectral dependence of normal reflectance for the three-layer system. The curve labelled ‘series, integration or ray tracing’ represents the result of all three methods when calculated to a sufficient precision.

Fig. 6

Trade-off between computation speed and precision for individual calculation procedures. The approximate formulae are represented with a single point since their precision is not tunable.

Fig. 7

Dependence of computation efficiencies of the series and direct integration on the magnitude of roughness. The mean times correspond to computations achieving the precision at least 10−6.

### Tables (1)

Table 1 Number of terms in the series for individual path lengths p and two to five-layer systems. Columns labelled ‘all’ represent all possible paths and correspond to Table 2 in [39]. Columns labelled ‘cons.’ represent the consolidated number of terms corresponding to the series presented in this work.

### Equations (38)

$X ^ j = 2 π λ N ^ j h j cos ϑ ^ j .$
$m = 2 p + 1 .$
$p = ∑ j = 1 d m j$
$v j > 0 , m j > 0 ,$
$m j ≥ max ( v j , v j + 1 ) or v j ≤ min ( m j , m j − 1 ) .$
$∑ p = 1 P ∑ d = 1 min ( L , p ) ∑ m ∑ v$
$∑ m = ∑ m 1 = 1 p − ( d − 1 ) − M 0 ∑ m 2 = 1 p − ( d − 2 ) − M 1 ⋯ ∑ m d − 1 = 1 p − 1 − M d − 2$
$∑ v = ∑ v 2 = 1 min ( m 1 , m 2 ) ∑ v 3 = 1 min ( m 2 , m 3 ) ⋯ ∑ v d = 1 min ( m d − 1 , m d ) .$
$M j = ∑ l = 1 j m l , M 0 ≡ 0$
$( t ^ j t ^ j ′ ) v j ( r ^ j ′ ) m j − v j r ^ j + 1 m j − v j + 1 exp ( 2 i m j X ^ j )$
$P 1 ( m 1 ) = ( t ^ 1 t ^ 1 ′ ) v 1 ( r ^ 1 ′ ) m 1 − v 1 r ^ 2 m 1 exp ( 2 i m 1 X ^ 1 ) = t ^ 1 t ^ 1 ′ ( r ^ 1 ′ ) m 1 − 1 r ^ 2 m 1 exp ( 2 i m 1 X ^ 1 ) ,$
$P 2 ( m 1 , m 2 , v 2 ) = t ^ 1 t ^ 1 ′ ( r ^ 1 ′ ) m 1 − 1 r ^ 2 m 1 − v 2 ( t ^ 2 t ^ 2 ′ ) v 2 ( r ^ 2 ′ ) m 2 − v 2 r ^ 3 m 2 exp ( 2 i m 2 X ^ 2 ) exp ( 2 i m 1 X ^ 1 ) .$
$P 3 ( m 1 , m 2 , m 3 , v 2 , v 3 ) = t ^ 1 t ^ 1 ′ ( r ^ 1 ′ ) m 1 − 1 r ^ 2 m 1 − v 2 ( t ^ 2 t ^ 2 ′ ) v 2 ( r ^ 2 ′ ) m 2 − v 2 r ^ 3 m 2 − v 3 ( t ^ 3 t ^ 3 ′ ) v 3 ( r ^ 3 ′ ) m 3 − v 3 r ^ 4 m 3 × exp ( 2 i m 3 X ^ 3 ) exp ( 2 i m 2 X ^ 2 ) exp ( 2 i m 1 X ^ 1 ) .$
$( m j − 1 v j )$
$( m j − 1 v j − 1 )$
$F ( v j , m j , m j − 1 ) = ( m j − 1 v j ) ( m j − 1 v j − 1 )$
$C d ( m 1 , m 2 , … , m d − 1 ; v 2 , v 3 , … , v d ) = C d ( m , v ) = ∏ j = 1 d F ( v j , m j , m j − 1 ) .$
$C 0 ≡ 1 , C 1 ( m 1 ) = C 0 F ( v 1 , m 1 , m 0 ) = ( 1 1 ) ( m 1 − 1 0 ) = 1 , C 2 ( m 1 , m 2 ; v 2 ) = C 1 ( m 1 ) F ( v 2 , m 2 , m 1 ) = ( m 1 v 2 ) ( m 2 − 1 v 2 − 1 ) , C 3 ( m 1 , m 2 , m 3 ; v 2 , v 3 ) = C 2 ( m 1 , m 2 ; v 2 ) F ( v 3 , m 3 , m 2 ) = ( m 1 v 2 ) ( m 2 − 1 v 2 − 1 ) ( m 2 v 3 ) ( m 3 − 1 v 3 − 1 ) .$
$r ^ = r ^ 1 + ∑ p = 1 P ∑ d = 1 min ( L , p ) ∑ m exp ( 2 i ∑ j = 1 d m j X ^ j ) Q ^ ( m ) ,$
$Q ^ ( m ) = t ^ 1 t ^ 1 ′ ( r ^ 1 ′ ) m 1 − 1 r ^ d + 1 m d ∏ j = 2 d ∑ v j = 1 min ( m j , m j − 1 ) F ( v j , m j , m j − 1 ) r ^ j m j − 1 − v j ( t ^ j t ^ j ′ ) v j ( r ^ j ′ ) m j − v j .$
$⌊ p 2 ⌋ ( p − ⌊ p 2 ⌋ ) ≈ p 2 4 .$
$t ^ j v j ( t ^ j ′ ) v j − 1 ( r ^ j ′ ) m j − v j r ^ j + 1 m j − v j + 1 exp ( i ( 2 m j − 1 ) X ^ j )$
$F ( v j , m j , m j − 1 ) = ( m j − 1 − 1 v j − 1 ) ( m j − 1 v j − 1 ) .$
$R = 〈 r ^ ( u ) exp ( i v t 1 ) 〉 〈 r ^ * ( u ) exp ( − i v t 1 ) 〉 = | 〈 r ^ ( u ) exp ( i v t 1 ) 〉 | 2 ,$
$v = − 4 π λ N 0 ,$
$〈 r ^ ( u ) exp ( i v t 1 ) 〉 = ∫ ∫ ⋯ ∫ r ^ ( u ) exp ( i v t 1 ) ρ ( u ) d u .$
$ρ ( u ) = 1 det S 1 ( 2 π ) ( L + 1 ) / 2 exp ( − 1 2 u T S − 1 u ) ,$
$S i , j = 〈 u i u j 〉 = σ i σ j C i , j ,$
$σ j 2 = 〈 u j 2 〉$
$C i , j = 〈 u i u j 〉 σ i σ j .$
$∫ ∫ ⋯ ∫ 1 det S 1 ( 2 π ) ( L + 1 ) / 2 exp ( i f T u ) exp ( − 1 2 u T S − 1 u ) d u = exp ( − 1 2 f T S f ) .$
$〈 r ^ ( u ) exp ( i v t 1 ) 〉 = r ^ 1 exp ( − v 2 σ 1 2 2 ) + ∑ p = 1 P ∑ d = 1 min ( L , p ) ∑ m exp ( 2 i ∑ j = 1 d m j X ^ j ) Q ^ ( m ) H ^ ( m ) .$
$X ^ j = 2 π λ N ^ j h j .$
$H ^ ( m ) = exp ( − 1 2 ∑ q = 1 d + 1 ∑ s = 1 d + 1 D ^ q D ^ s S q , s ) ,$
$D ^ j = 4 π λ ( m j N ^ j − m j − 1 N ^ j − 1 ) .$
$u k ( x , y ) = u k + 1 ( x , y ) + ε k ( x , y ) for 1 ≤ k ≤ L ,$
$〈 u L + 1 ( x , y ) u L + 1 ( x , y ) 〉 = σ L + 1 2 , 〈 ε k ( x , y ) ε j ( x , y ) 〉 = δ j k σ 2 , 〈 ε k ( x , y ) u L + 1 ( x , y ) 〉 = 0 .$
$S j k = σ L + 1 2 + ( L + 1 − max ( j , k ) ) σ 2 .$