Adhesion of Soft Materials to Rough Surfaces: Experimental Studies, Statistical Analysis and Modelling

Adhesion between rough surfaces is an active field of research where both experimental studies and theoretical modelling are used. However, it is rather difficult to conduct precise experimental evaluations of adhesive properties of the so-called anti-adhesive materials. Hence, it was suggested earlier by Purtov et al. (2013) to prepare epoxy resin replicas of surfaces having different topography and conduct depth-sensing indentation of the samples using a micro-force tester with a spherical smooth probe made of the compliant polydimethylsiloxane polymer in order to compare values of the force of adhesion to the surfaces. Surprising experimental observations were obtained in which a surface having very small roughness showed the greater value of the force of adhesion than the value for a replica of smooth surface. A plausible explanation of the data was given suggesting that these rough surfaces had full adhesive contact and their true contact area is greater than the area for a smooth surface, while the surfaces with higher values of roughness do not have full contact. Here, the experimental results of surface topography measurements and the statistical analysis of the data are presented. Several modern tests of normality used showed that the height distribution of the surfaces under investigation is normal (Gaussian) and hence the classic statistical models of adhesive contact between rough surfaces may formally be used. Employing one of the Galanov (2011) models of adhesive contact between rough surfaces, the plausible explanation of the experimental observations has been confirmed and theoretically justified.


Introduction
Molecular adhesion is a universal physical phenomenon that usually has a negligible effect on surface interactions at the macro-scale, whereas it becomes increasingly significant as the contact size decreases [1]. Although the main topic of recent studies of adhesive interactions has been devoted to materials showing enhanced adhesive properties, the studies of the so-called anti-adhesive materials are also a growing field of interest. The term 'anti-adhesive' is attributed to materials and coatings showing very low values of adhesive forces. However, anti-adhesive properties should not be attributed solely to materials having low surface-energy because surface roughness may values of the adherence force. Evidently, to model adhesion between rough surfaces, one needs to employ corresponding models [3,11,15]. In turn, there is a need for proper statistical analysis of surface topography. We give below some preliminary information related to various approaches to description of surface roughness along with a brief review of models of contact between elastic rough surfaces. Both adhesive and non-adhesive models are discussed. Special attention is given to models developed in [3] and two models introduced in [11].

Statistical Description of Rough Surfaces
Various parameters and functions related to the surface asperity shapes and distributions of the asperities were suggested to be used as statistical surface characteristics [16]; however, Whitehouse [17] noted that some of these parameters are useful, but most are not.
Let us consider a nominally flat surface. If one considers a plane perpendicular to the surface, then the cross-section between the surface and the plane is a profile. Let a function z(x) describe the rough profile. The origin level of the height measurements can be taken as the mean profile linez. Then, one has 1 2L The most popular statistical height parameters for a function z(x) defined on an interval [−L, L] are the maximum height of the profile R max , the arithmetical mean deviation of the surface R a , and the square root of the mean square deviation σ 2 (the rms height R q ) with respect to the mean profile linē where n is the number of points of measurements on the interval and z(x i ) is the measured height at the interval point x i . In practical calculations, the arithmetic mean height R z is often used; it is calculated as the average distance between the five highest picks and the five lowest points of the profile, i.e., One can introduce the density probability function φ(z) that shows the probability that the height z(x) at a surface point x is between z and z + dz. Then, the expressions for R a and σ 2 in Equation (1) can be written as Other popular parameters of roughness include the rms slope or the rms curvature, parameters associated with horizontal distributions, e.g., the number of intersections of the surface with the average line, and parameters describing spatial extent of asperities, such as the high spot count, and so on. It is quite clear that the cumulative distribution function of the surface heights Φ(z) is correlated with contact properties of the rough surfaces In tribology, this function is known as the bearing area curve or the Abbott-Firestone curve. Actually, Abbott and Firestone [18] suggested calculating the integrals of the horizontal line at a specific level, which lies within the roughness profile. Its value at a level z = h is equal to the length (the area in the 2D problem) of the slice of the profile at the level h.
We have to mention that the original rough surfaces of polishing papers were made using particles of diameters 0.3, 1, 3, 9, 12 µm and these diameters were used to specify the nominal asperity sizes. Using the above definitions from Equations (1) and (2), the replica surfaces were examined in [5]. The results obtained are presented in Table 1. The first line is attributed for a replica of smooth glass surface. As it was defined above, R a is the roughness average (the absolute value of the surface height averaged over the surface); rms is root mean square (a statistical value of the magnitude of a varying quantity) and R z is determined roughness, i.e., the average of N = 10 extremal roughness heights over a specified length. The values of Table 1 were obtained using the arrays of discrete values obtained by measurements in specific points of profiles. However, one can also consider the density probability function of asperity summits φ(z s ) and the cumulative summit height distribution Φ(z s ). These functions are used in models of discrete contacts between rough surfaces.
The above characteristics are mainly related to vertical distribution of rough profiles. However, the horizontal distribution of the profiles is also very important. The surface roughness profiles may be modelled as graphs of a random process. As early as 1953, Linnik and Khusu delivered a seminar where they suggested to describe surface roughness using graphs of a stationary Gaussian random process having a correlation function where K(0) and α are some parameters of the roughness (see, for details, an abstract of the seminar [19], as well as [20] and [13]). Practically all known models of contact between rough adhesive surfaces are based on explicit or implicit assumption of the Gaussian distribution of heights (see a review in [21]). The Gaussian distribution has an attractive property that, if the random surface is normal, then it can be completely described by two parameters: (i) the height distribution and (ii) an auto-correlation function [14]. The latter function characterizes the horizontal distribution of asperities of a rough surface profile. Indeed, one needs to specify both vertical and horizontal distributions because, as it was clearly formulated in [15], two profiles may have the same height and peak height (local extrema) distributions, but they may differ in the horizontal extension. The profile auto-correlation function is defined as Evidently, it is possible to putz = 0 in Equation (6). The correlation (auto-correlation) function is often substituted by its Fourier transform that is called the spectral density function.
Various mathematical tools are developed to work with surfaces having the Gaussian (normal) distribution of heights (see, e.g., [13,15]). However, to apply these tools, one needs to demonstrate that the surface roughness follows the Gaussian law. Quite often, these tools are used without proper justification that they may be used, while it is not evident that the roughness height distribution is normal. In fact, the accurate statistical analysis of grinding surface roughness showed that the distribution is not normal either at micro or nanoscales [21]. Therefore, the results of an accurate statistical analysis of surfaces are presented below. As it has been mentioned, the test surfaces were the same as described in [5]; namely, these were the epoxy resin replicas of polishing papers produced by Serva Electrophoresis GmbH (Heidelberg, Germany) with nominal asperity sizes 0.3, 1, 3, 9, and 12 µm along with a replica of a clean smooth glass surface. The glass surface having R a asperity roughness of 7 nm was considered as a nominally smooth one.

Non-Adhesive Contact between Rough Elastic Solids
Apparently, the first statistical model of contact between rough surfaces was presented in 1940 by Zhuravlev [22]. He assumed that all asperities are elastic spheres of the same radius R located at various heights, i.e., it was a statistical multi-level model of roughness. He argued that profiles at a specific height may be characterized by a distribution function ψ(ξ), i.e., the number of spheres in contact increases as one considers the deeper levels ξ. Applying the Hertz contact theory, he derived general expressions for the true contact area A between the surfaces and the contact force P: Here, k = (1 − ν 2 )/πE is the effective contact modulus of surfaces having the elastic modulus E and the Poisson ratio ν, δ is the relative approach due to compression of the surfaces, and N is the total number of summits of the protuberances located at depth δ. As an example, Zhuravlev considered that φ(ξ) is linear. In this case, it follows from Equation (8) that A ∼ P 10/11 . Kragelsky [23] modified the Zhuravlev approach by modelling a rough surface as a collection of elastic flat ended rods of various heights based on a flat rigid surface. He noted that one would need to employ the numerical integration of general expressions, if the Gaussian distribution of rod heights were employed. Although Zhuravlev's statistical multi-level model of surface roughness is rarely cited, it was very successful (see [24] for further details). Indeed, the model was reborn and developed further as the Greenwood and Williamson (GW) model. Greenwood and Williamson [25] noted that, in the framework of the Hertz contact theory, contact between two spherical elastic solids is equivalent to contact between a sphere of an effective radius and a plane; therefore, they do not consider contact between two spherical asperities but rather contact between a spherical elastic asperity and a rigid plane. Contrary to Zhuravlev, they put the origin of the vertical coordinate not at the maximum height, but rather they took the mean plane of summit height as the reference.
As advantages of the GW model, Greenwood [26] noted the following two points: it improves on the Zhuravlev model by assuming a Gaussian height distribution (still of identical spherical caps) and it considers not only elastic asperities, but also plastic deformation of asperities. However, the former advantage may be rather illusive at nanoscale. Indeed, modern nanotechnology considers surfaces whose roughness is below the micrometre scale, while plasticity effects for crystalline materials may be observed only at the micrometre scale and they do not arise at atomic and nanometre scales due to the Polonsky-Keer effect. One of the authors (FB) has introduced this notation for the effect to recognize the contribution of Polonsky and Keer [27,28] who noted first that, if the contact region of a crystalline asperity is below some characteristic size of the microstructural length, the number of dislocations is not large enough to cause plastic flow. This effect can be formulated as follows: "plastic deformation at an asperity micro-contact becomes difficult and then impossible when the asperity size decreases below a certain threshold value on the order of the microstructural length". Hence, in crystalline solids, the plastic deformations of asperities may appear only within microscale asperities or at larger scales while nanoscale asperities deform elastically (see a discussion in [29]). Furthermore, as it has been mentioned above, the Gaussian height distribution is not observed for grinding surfaces either at micro or nanoscales [21]. As it will be shown below, this distribution may be observed on intact surfaces. A specific feature of the GW approach is the introduction of an exponential height distribution of the asperity summits that allowed them to show that the true contact area A is proportional to the applied compressive load P.
We have to add that it is a difficult task to transfer from the characteristics of a profile (or several profiles) to the characteristics of a surface [15]-for example, to transfer the density of peaks of a surface profile η p to the density of asperity summits per unit area η s of the nominally flat surface. If φ(z) is Gaussian, then If the distribution of the summit height z s is Gaussian with a standard deviation σ s , i.e., it is characterized by the normalized distribution function (φ s ) n then the GW model is specified by three parameters: the effective radius of asperity R, the density of asperity summits per unit area η s and the standard deviation of the summit heights σ s . It is worth commenting that these parameters are not independent, but rather they should be related by the equation [14,15]: σ s η s R = const.

Adhesive Contact between Rough Elastic Solids
Johnson [30] presented the first model of adhesive contact between rough surfaces. Although Johnson mentioned both Zhuravlev [22] and Greenwood and Williamson [25] models, he decided to use the exponential distribution of summit heights of the GW model. Independently, Fuller and Tabor [3] introduced the FT model (see also [31]). They studied the adhesive contact between rough surfaces whose asperity heights follow the Gaussian distribution. In fact, the above models are extensions of the GW model to the case of adhesive contact of independent spherical asperities. The FT model is discussed in detail in [15]. This and similar models are still very popular. The statistical approaches employed by many international teams involved in the recent competition "Contact Mechanics Challenge" (see [32]) were based either on the use of the random process techniques applied to the surface roughness (e.g., the auto-correlation function) or some integral parameters of the surface roughness, e.g., the root mean square and the asperity radius of curvature.
It is worth mentioning that fractal approach to contact problems for rough surfaces was quite popular some time ago (see a review by Borodich [33]). However, Borodich has introduced some models that showed that fractal dimension of surfaces alone cannot characterize contact properties of surfaces. These models include the Cantor-Borodich profile also known as the Cantor set model [34], multi-level hierarchical model [35,36] and indenters having parametric-homogeneous (PH) profiles [37,38]. The theory of parametric-homogeneous functions developed by Borodich is discussed in detail in his papers [39,40].
Let us discuss two models of adhesive contact between rough surfaces introduced by Galanov [11]. The Galanov models of contact between rough solids represent the surface roughness as a layer whose mechanical properties differ from the bulk properties of the material (see also [41]).
It is argued in the first model that the adhesive properties of rough surfaces may be defined just by introducing the specific work of adhesion w rough that can be written in the following dimensionless form: where P c is the adhesion force, i.e., the pull-off force for a single asperity. According to the Johnson-Kendall-Roberts (JKR) contact theory, P c = 1.5 πwR, where w is the specific work of adhesion of contacting pair of materials and R is the effective radius of contacting asperities.
In the refined model, the roughness is simulated as a spring mattress whose springs have nonlinear force-displacement diagram and they stick to the opposite surface. Non-adhesive contact problems for such layer were discussed by many authors (see, e.g., [42]). In the linear case, this layer is the Winkler-Fuss elastic foundation that may stick to the surface. It is well known that the leading terms of relations for the problem of contact between a compressible thin layer and a blunt punch may be treated as a problem for the Winkler-Fuss elastic foundation; it is also possible to solve the appropriate adhesive contact problem (see, e.g., [43]).
Furthermore, we will employ the former model of adhesive contact [11]. As many other models of adhesive contact, this model employs implicitly the Derjaguin approximation (see a discussion in [44]) and therefore the interaction energy per unit area between small elements of curved surfaces is the same as the interaction energy per unit area between two parallel infinite planar surfaces. For this Galanov model, the specific work of adhesion w rough is defined as work spend for separation of a unit area of the punch from the rough surface (see Figure 1) where P denoted the force, p is the nominal contact pressure, P * = P/(NP c ) and p * = p/(η s P c ) are the dimensionless force and pressure, and N is the number of asperities in contact. It is assumed that one can specify the distance t a between the conventional boundary of adhesive interactions and the horizontal axis. In this case, if ∆ is the change of the spring length of the rough layer, then the distance d between the smooth plane and the horizontal axis is defined as d = t a − ∆. This distance can be normalized by the standard deviation For classic unilateral contact problem Contrary to unilateral contact, displacements of the probe may be both positive (compression) and negative (material of the sample follow the probe due to sticking to the probe). Hence, we need to introduce several parameters to describe adhesive contact. In addition to P c , we need to specify δ c as the absolute value of a minimum of an asperity displacement and its normalized value δ * c = δ c /σ. The latter value is defined according to the JKR theory as The effective contact modulus is E * = 69 · 10 3 Pa for the polydimethylsiloxane and epoxy resin replicas used in [5]. Then, one can estimate the specific work of adhesion using the JKR theory and experiments for smooth surface: w = 0.403 N/m (methods of experimental extraction w are discussed in [8,45]).
According to the above definition, we can write Here, d * 0 is the root of the nonlinear equation i.e., p * (d * 0 ) = 0, while the function p * (z * ) is defined for a specified δ * c = δ c /σ as (see (4.334) in [15]) Thus, Equations (11)-(13) provide the full system of equations for calculation of w * rough for any value of δ * c . To employ the above model, one needs to extract from the experiments σ, R, and w, and calculate η s and the value σ s η s R. According to [15], 0.03 < σ s η s R < 0.05.

Statistical Analysis of Experimental Surface Topography Data
As it has been mentioned, the measurements of surface topography have been fulfilled by white light interferometry. We need to note that only data obtained for test surfaces with nominal asperity size 0.3 and 1 µm was good for statistical analysis. The data obtained for surfaces with higher nominal asperity sizes were good enough only for very small domains that were separated by domains with no proper data obtained by the interferometer. Indeed, one can see in Figure 2 that there are domains where the white light interferometer was not able to provide proper data.

Normality Tests
As it has been shown above (see also a review in [21]), the crucial assumption of many models of contact is that the surface roughness is Gaussian. Hence, it is important to check if the height distribution of the tested surfaces follow the normal distribution. A large number of methods for testing for normality [46] has been developed. Each normality test is based on a particular test statistic that is used to give a quantitative estimation of proximity between an observed sample of measurements and the theoretical normal distribution. The assumption of normality is especially critical when characteristics of roughness profiles are derived. The main methods for testing for normality are the following tests: Kolmogorov-Smirnov (KS), Lilliefors (LF), Shapiro-Wilk (SW), Anderson-Darling (AD), Cramer-von Mises (CVM), the Pearson, and Shapiro-Francia (SF). Using the test statistic, each of the above tests produces a numerical characteristic that is called the p-value.
This number characterizes the significance at the scale [0, 1] that the hypothesis of the normality is true for the observed measurements. One can nominate the acceptable significance level, say 5%. If the p-value is less than this level, then the hypothesis of normality should be rejected. Alternatively, if the number is within the significance level, then one can conclude that the height distribution is normal. Note that the acceptable significance level is a probability to reject the hypothesis of normality even if it is true.
To describe the tests, we introduce notations: z (1) , . . . , z (n) is a permutation of the sample z 1 , . . . , z n such that z (1) ≤ . . . ≤ z (n) , p (i) = Φ((z (i) −z)/s), where Φ is the cumulative distribution function of the standard normal distribution, andz and σ s are mean and standard deviation of the sample.

The KS Test
The Kolmogorov-Smirnov (KS) test for a sample is based on comparison of an empirical distribution function (edf) F N (t) and the theoretical cumulative distribution function (cdf) F(t) of the test distribution; specifically, the test statistic is where F(z) is the cdf of Gaussian (normal) distribution and is the edf of the sample (z 1 , z 2 , . . . , z N ). A limitation of the KS test is its high sensitivity to extreme values.

The LF Test
The Lilliefors (LF) test is a correction of the KS test in a way to be more conservative. The test statistic is the maximal absolute difference between the empirical and theoretical cdf as The LF test is slightly better than the KS test because the LF test is more conservative.

The AD Test
The Anderson-Darling (AD) test is based on a squared difference between the edf and theoretical cdf, and the test statistic is The Anderson-Darling test is sensitive to discrepancies in the tails of the distribution.

The CVM Test
The Cramer-von Mises (CVM) test is also based on a squared difference between the edf and theoretical cdf. The test statistic is The CVM test is uniform to discrepancies in different parts of the distribution.

The SW Test
The Shapiro-Wilk (SW, or Shapiro) test is based on the correlation between the data and the corresponding normal scores The Shapiro-Wilk test has better power than the KS test and the LF test. The Shapiro-Wilk test is known as the best choice for testing the normality of data.

The Pearson Test
The Pearson chi-square test is based on creating the classes built in such a way that they are equiprobable under the hypothesis of normality, where the number of classes is chosen as m = 2n 2/5 . The Pearson statistic is given by where C j is the number of counted observations and E j is the number of expected observations under the hypothesis of normality in the j-th class.

The SF Test
The Shapiro-Francia (SF) test is simply the squared correlation between the ordered sample values and the (approximated) expected ordered quantiles from the standard normal distribution, and the test statistic is where µ i are standard normal ordered statistics.

Results of the Normality Testing
In Figure 3, we show the roughness analysis for data of the tested surface with nominal asperity size 0.3 µm obtained by white light interferometer. The original profile measurements are given by a black curve. A waviness is shown by a red curve. The profile measurements with removed waviness are given by a blue curve. The height distribution for the original profile is depicted by a histogram on the left and for the profile with removed waviness on the right.
In Figure 4, we show the statistical analysis for roughness data of the tested surface with nominal asperity size 1 µm. The meaning of notations is the same as in Figure 3.
Several tests for normality show that the height distribution is normal. One can observe that almost all p-values are greater that 0.05 and therefore we accept the hypothesis of normality for these measurements. Because, as it has been shown above, the surfaces with nominal asperity size 0.3 and 1 µm do follow the Gaussian height distribution, it will be assumed further that the surfaces with higher nominal asperity sizes follow the same distribution.
Thus, we have confirmed normality of the height distribution of these intact surfaces of the replicas, i.e., the assumption of normality used by the above mentioned GW, FT and Galanov models of contact between rough surface was not violated for these surfaces. Hence, these models can be used for studies of the DSI tests of replicas of the intact surfaces of polishing paper.

Experimental Data Obtained by DSI
As it has been mentioned above, the force tester Basalt 01 was used in the DSI tests. Epoxy-resin replicas of a smooth glass surface and five rough surfaces were used as substrates. For measurements of adhesion forces, a micro-force tester equipped with a double-leaf spring serving as a cantilever beam was applied. Here, a soft spherical probe was pressed with a given force onto a sample. The required force to release the probe from the surface (pull-off force) corresponds to the adhesive property of the sample.
If the surfaces were smooth, then, in the case of a spherical probe, the real contact area and effective elastic modulus could be calculated according to the Johnson-Kendall-Roberts (JKR) theory, which is a well established standard analyze method for micro-scale adhesion between elastic materials.
The JKR theory was developed for adhesive contact between two smooth elastic spheres or a rigid sphere and an elastic half-space (see for details [1,44,47]). However, here we consider contact interactions between relatively rigid epoxy resin replicas and a soft polydimethylsiloxane sphere. These are polymer materials. Properties of such materials are time-dependent and their elasticity is often described by Volterra integral equations [48]. However, if the kernel of the Volterra equation is not singular (the physical meaning of this condition is that the rearrangements of long polymer chains in both materials are not very fast), then, for short times for the loading-unloading tests, the materials may be quite often treated as elastic ones. Furthermore, adhesion of polymers may have different physical origin (see [49] and also a comprehensive review [50]). However, if the loading to the prescribed maximal compressive force and unloading of the samples take the same short amount of time for all samples, then it is quite convenient to use "the specific work of adhesion of contacting pair of materials" and use these values to compare the results of the tests.
The experimental absolute value of the pull-off force for replica of smooth glass surface was P exper c = 2845 µN. Absolute average values of measured pull-off forces are presented in the 4th column of Table 2. One can see that the highest average of the pull-off force value P exper c = 3823 µN was obtained not on the smooth sample, but on the rough surface whose nominal size of surface asperities was 0.3 µm. The pull-off forces decreased continuously with the further increase of the surface micro-roughness [5]. Here, we will provide some theoretical justifications of the plausible explanation of the experimental data given in [5].
First, we note that the greater the area of the true contact, the stronger the adhesion absolutely irrespective to the nature and the value of the forces carrying out the adhesive connection [51]. Assuming that there is almost full contact between the smooth spherical probe and the tested surface for both the nominally flat sample and the surface having the nominal asperity size of 0.3 µm, one can estimate the full area of the rough surface.
Of course the Galanov model has the same limitations as the Zhuravlev, GW, and FT models of discrete contact: it is assumed that contacts between asperities do not influence each other. Hence, if the rough surface of the nominal asperity size of 0.3 µm has full contact with the probe, then the model cannot be applied. Let us estimate the total area of the rough surface.
Khusu et al. [13] studied several correlation functions for rough profiles having Gaussian height distributions. The asperity heights were calculated from the middle line of the profile, that, in turn, was calculated by the least square method. This was leading to introduction of additional parameters. Using this description, Linnik and Khusu [20] measured various parameters. In particular, they studied (i) the fraction of the horizontal line at a specific level, which lies within the roughness profile, i.e., they calculated the Abbott-Firestone curve, and (ii) the total area of the profile curve above this horizontal line. They reported pretty good agreement of the obtained results with their particular experimental observations. They suggested also to use the Abbe test for checking that the expectations of all values are the same, i.e., there are no trends in data. In addition to K(x) of form in Equation (5), which is a correlation function of a homogeneous isotropic random field, they also studied the correlation function where K 2 (0) and α 2 are some parameters of the roughness. The correlation functions given by Formulas (5) and (14) will be used here to estimate the true area of rough surfaces. From the real data of roughness over the rectangular (0.11 mm × 0.14 mm), we compute the number of summits N and the number of summits per unit area η s . Taking σ s η s R = 0.035, one can calculate the radius of curvature of an asperity as R = 0.035/(σ s η s ). The results of calculations are presented in Table 3.
The calculations were fulfilled in the following steps. Using the values of δ * obtained, we calculated the values of w * rough by Equations (11)-(13) of the Galanov model. Finally, we compute the specific adhesion work of rough surface w rough = ση s P c w * rough and the real force of adhesion between a smooth sphere (probe) of radius R s = 1.5 mm and rough surface as P rough c = 1.5πR s w * rough . The 4th column of Table 2 shows the average of the experimental pull-off force P exper c [5]. These experimental values are close to the adhesion force given in the 3rd column calculated for contact between a smooth spherical probe of radius R s = 1.5 mm and rough surface when the surface has nominal size of asperities larger or equal to 3 µm. However, the calculated values P rough c are not close to the experimental values for surfaces with nominal size 0.3 µm and 1 µm because there exists almost full contact and the Galanov model along with the FT or GW models are not applicable.
In the case of full contact, it is useful to compare the true area of the rough surface and the smooth one.
The length of a realization of a random Gaussian process z(x) on a unite interval is If the standard deviation of the process σ/L << 1 and one employs the correlation function given by Equation (14), then one can write an asymptotic series for its expectation Es (see Chapter 2.9 in [13]). It follows from Table 1 that rough profile having nominal asperity height 0.3 µm has σ = 0.205 µm over an interval L = 140 µm, i.e., σ/L = 0.205/140 << 1. Hence, the asymptotic series may be used. One can get that the estimations of the profile length are very small in the case if the correlation function given by Equation (14) is employed. We can conclude that Equation (14) does not reflect the horizontal distribution of asperities. However, our computations using numerical simulations of a Gaussian process with the exponential correlation function given by Equation (5) show that the length of the profile increases. Simulating the surface with the nominal size of asperities 0.3 µm and computing the length of profiles along several lines, we have obtained that s/L = 1.16. If we estimate the ratio of the true area of the rough surface to the nominal area as A true /A nominal = (s/L) 2 , then we obtain A true /A nominal = 1.16 2 = 1.3456. On the other hand, the ratio of P exper c for the surface with nominal asperity heights 0.3 µm to the pull-off force for the smooth one is 3.823/2.845 = 1.3438. These ratios are in excellent agreement between each other. Thus, the plausible explanation of the observations given in [5] are correct and supported here by our theoretical considerations.

Conclusions
Currently, the design of many devices is influenced by discoveries of structures of biological objects. It is known that attachment organs of insects have worse performance at surfaces having specific roughness and this idea is used in preparation of the so-called anti-adhesive materials. To study effects of surface roughness on molecular adhesion between solids, it was suggested earlier to prepare epoxy resin replicas of surfaces having different roughness and conduct contact loading tests of the samples using a micro-force tester with a spherical smooth probe. Because the materials of all samples are the same, the roughness effects are not influenced by other factors. The contact loading tests between different samples showed a rather surprising results: a surface having very small roughness showed the greater values of the force of adhesion than the replica of a smooth surface. A plausible explanation of the observations was given in [5], suggesting that these rough surfaces had full adhesive contact while their true contact area is greater than the area for a smooth surface, and the surfaces with higher values of roughness do not have full contact.
One of the main purposes of the paper is to justify theoretically these explanations. Evidently, one needs to use a theory of adhesive contact between rough surfaces. Unfortunately, the main theories of contact between rough surfaces were developed for roughness whose height distributions are Gaussian. However, one cannot assume that the surfaces always have normal (Gaussian) roughness. For example, it was shown that the grinding surfaces are not Gaussian either at micro or nanoscale [21].
To check the assumption that the roughness is Gaussian, various modern tests of normality of data have been employed. These tests include the Kolmogorov-Smirnov, Lilliefors, Shapiro-Wilk, Pearson, Anderson-Darling, Cramer-von Mises and Shapiro-Francia ones. Of course, there are many other tests of normality (see a discussion in [46]). However, the tests used are the most popular ones. These tests showed that the height distribution of the surfaces under investigation are normal (Gaussian).
We can conclude that, for replicas of polishing papers, the exponential correlation function given by Formula (5) may be used to model horizontal distribution of asperities having Gaussian vertical distribution. The polydimethylsiloxane sphere is soft enough material to provide full contact for the surface with nominal asperity heights 0.3 µm, while, for the used load, the full contact was not observed for surfaces having higher roughness. Using one of the Galanov models of adhesive contact between rough surfaces, the experimental observations are explained for surfaces whose nominal asperity values were over 1 µm. For a surface with nominal asperity heights 0.3 µm, we have provided the results of numerical simulations of random surfaces and justified the above explanations using geometrical arguments. Thus, the plausible explanation of observations have been theoretically justified.