This article is published under a Creative Commons license, not by the author of the article. So if you find any inaccuracies, you can correct them by updating the article.


Comprehensive topography characterization of polycrystalline diamond coatings Creative Commons

Link for citation this article Add this article in bookmark list
Abhijeet Gujrati, Department of Mechanical Engineering and Materials Science, University of Pittsburgh, 3700 O'Hara St, Pittsburgh, PA, 15261, United States of America
Antoine Sanner, Department of Microsystems Engineering, University of Freiburg, Georges-Köhler-Allee 103, 79110 Freiburg, Germany; Cluster of Excellence livMatS, Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
Subarna R. Khanal, Department of Mechanical Engineering and Materials Science, University of Pittsburgh, 3700 O'Hara St, Pittsburgh, PA, 15261, United States of America
Nicolaie Moldovan, Alcorix Co., 14047 Franklin Ct., Plainfield, IL 60544, United States of America
Hongjun Zeng, Alcorix Co., 14047 Franklin Ct., Plainfield, IL 60544, United States of America; Aqua via Rock LLC, 712 Muirhead Ct., Naperville, IL 60565, United States of America
Lars Pastewka, Department of Microsystems Engineering, University of Freiburg, Georges-Köhler-Allee 103, 79110 Freiburg, Germany; Cluster of Excellence livMatS, Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, Georges-Köhler-Allee 105, 79110 Freiburg, Germany
Tevis D. B. Jacobs Department of Mechanical Engineering and Materials Science, University of Pittsburgh, 3700 O'Hara St, Pittsburgh, PA, 15261, United States of America
Surface Topography: Metrology and Properties, Journal Year: 2021, Volume and Issue: 9(1), P. 1 - 1

Published: Jan. 1, 2021

This article is published under the license License

Link for citation this article Related Articles


The surface topography of diamond coatings strongly affects surface properties such as adhesion, friction, wear, and biocompatibility. However, the understanding of multi-scale topography, and its effect on properties, has been hindered by conventional measurement methods, which capture only a single length scale. Here, four different polycrystalline diamond coatings are characterized using transmission electron microscopy to assess the roughness down to the sub-nanometer scale. Then these measurements are combined, using the power spectral density (PSD), with conventional methods (stylus profilometry and atomic force microscopy) to characterize all scales of topography. The results demonstrate the critical importance of measuring topography across all length scales, especially because their PSDs cross over one another, such that a surface that is rougher at a larger scale may be smoother at a smaller scale and vice versa. Furthermore, these measurements reveal the connection between multi-scale topography and grain size, with characteristic scaling behavior at and slightly below the mean grain size, and self-affine fractal-like roughness at other length scales. At small (subgrain) scales, unpolished surfaces exhibit a common form of residual roughness that is self-affine in nature but difficult to detect with conventional methods. This approach of capturing topography from the atomic- to the macro-scale is termed comprehensive topography characterization, and all of the topography data from these surfaces has been made available for further analysis by experimentalists and theoreticians. Scientifically, this investigation has identified four characteristic regions of topography scaling in polycrystalline diamond materials.


Diamond coatings, Comprehensive topography, conventional methods

1. Introduction


Surface topography controls surface properties of carbon coatings. For example, prior measurements of diamond-like carbon (DLC) coatings show that tribological behavior [1] and adhesion [2] are strongly affected by surface texture, all the way down to the nanoscale. The surface topography of diamond coatings [3] affects their performance, including their friction [4], wear [5], adhesion [6, 7], and biocompatibility [8]. Diamond is an important material that is used in many industrial applications [9] such as for mechanical seals [10], MEMS devices [11], biomedical applications [12], seals and bearings [13], and nuclear fusion [14] because it has low friction and wear [15, 16], and because it is robust and chemically inert so that it can be operated in corrosive environments.

Numerical models have been proposed to describe the effect of surface roughness on contact and adhesion [17]. Initially, the classic Greenwood-Williamson [18] and Fuller-Tabor [19] models described contact area and adhesion based on the average height of the roughness (figure 1(a)). Recently, the multi-scale nature of roughness has been included in the modeling of surface properties [20–29]. In particular, it has been shown that a critical quantity controlling contact area, adhesion, friction, and wear is the root-mean-square slope of the surface h'rms and that this quantity, for a surface with multi-scale roughness, is strongly influenced by small-scale features [22, 30]. Therefore, understanding and prediction of the topography-dependent surface properties of carbon coatings requires multi-scale characterization of roughness, down to the atomic scale.

Because conventional methods for topography measurement cannot capture these small scales, then researchers must choose either to ignore the small-scale topography, or to extrapolate it from the single-scale measurements at larger scales. Because many real-world surfaces exhibit hierarchical topography [32–36], this extrapolation is often done by assuming that the material is self-similar or self-affine. Self-similarity implies that the topography is statistically indistinguishable at all magnifications ζ; in other words, if the lateral length scale L is rescaled to ζL, then the measured height h is rescaled to ζh. Self-affinity is related, but characterized by the Hurst exponent H; where rescaling the length to ζL and the height to ζH h yields statistically indistinguishable surfaces. Mathematically, the root-mean-square slope h'rms of a self-affine, randomly-rough surface depends strongly on the smallest scales at which the roughness exists [22, 30]. For real surfaces, it is often assumed that surfaces are self-affine down to the atomic scale, but this assumption is largely untested.

Instead of ignoring or extrapolating for small-scale topography, new approaches are required for characterizing topography across all length scales. While stylus profilometry and atomic force microscopy (AFM) are indispensable tools in characterizing surface topography, these tip-based techniques are unable to provide the small-scale topography of rough surfaces because the radius of the scanning tip introduces artifacts [37]. This limits the range of reliability of the roughness measurements [30] (figure 1(b)). Optical techniques, such as scanning white-light interferometry or laser confocal microscopy, suffer from diffraction-limited lateral resolution and optical-transfer-function artifacts, and are thus similarly incapable of measuring the smallest-scale topography [38]. The result of this is that conventional measurements of surface topography are incomplete, and computed surface metrics (such as h'rms) are unreliable, depending explicitly on the lateral resolution of the measurement. Even very advanced methods of analyzing surface topography, such as those in [39], are limited in their effectiveness by the range of size-scales in the underlying topography measurement. Instead, cross-section electron microscopy provides a reliable method to characterize surface topography down to the Ångström-scale [40] (figure 1(c)). Further, the small-scale topography can be stitched together with the medium- and large-scale topography using the power spectral density (PSD), to provide a comprehensive statistical description of surface topography at all size scales [32]. The PSD is a mathematical tool which separates the contribution to roughness from different length scales  and it is commonly represented as a function of wavevector 

This method of combining many different measurements, from the atomic to the macroscale, is termed comprehensive topography characterization and was applied here to investigate the surface roughness of four different varieties of diamond coatings, namely ultrananocrystalline diamond (UNCD), polished UNCD (pUNCD), nanocrystalline diamond (NCD), and microcrystalline diamond (MCD).

2. Methods

Thin films of the diamond materials were deposited (Advanced Diamond Technologies, Romeoville, IL) using a tungsten hot-filament chemical vapor deposition (HF-CVD) system with parameters as described in [41]. To improve electrical conductivity, all materials were boron doped with a B/C ratio of 3000 ppm. All materials were deposited to a thickness of 2 μm on polished silicon wafers after the wafers were sonicated with slurries containing suspended diamond nanoparticles. The deposition parameters are listed in table 1.

The smallest-scale topography measurements were made using transmission electron microscopy, following the techniques described in [40]. For the UNCD, NCD, and MCD, the 'wedge deposition technique' was used, whereas for pUNCD, the 'surface-preserving cross-section technique' was used. These techniques are described in detail in [40]. Briefly, the wedge deposition technique involves depositing the diamond film directly onto TEM-ready silicon thin-wedge substrates. The surface-preserving cross-section technique utilizes conventional methods for the preparation of TEM cross-sections (sectioning, grinding, polishing, dimple-grinding, and ion etching) with process modifications to ensure that the original surface is preserved. The samples were imaged using a TEM (JEOL JEM 2100F, Tokyo, Japan) operated at 200 keV. The images were taken using magnification levels from 5000× to 600 000×.

The profiles were extracted from the TEM images by using custom Matlab scripts to trace the outermost boundary of the material. Before tracing, the images were rotated such that the boundary was approximately horizontal. The vast majority of the measured surfaces were functions, i.e. for every point on the x-axis, there was exactly one point on the y-axis; in other words, the measured topographies were not reentrant. However, there were some cases where two adjacent points were captured with identical or decreasing horizontal position (locally reentrant). For these cases, the latter point was removed to restore non-reentrant behavior so that the topography can be described by a function as required by the calculation of the PSD. They were only observed in 12 out of 160 TEM profiles, and even then, only in a small number of points per profile, therefore the process of removing such points should not affect the accuracy of the analysis. Further, these locally reentrant points are attributed to imperfect rotation of the TEM profiles, rather than to truly reentrant features on the surface.

The medium-scale topography was measured using an atomic force microscope (Dimension V, Bruker, Billerica, MA) in tapping mode with diamond-like carbon-coated probes (Tap DLC300, Mikromasch, Watsonville, CA). For all substrates, square measurements were taken with the following lateral sizes: 3 scans each at 100 nm, 500 nm, and 5 μm; 1 scan each at 250 nm and 1 μm. The scanning speed was maintained at 1 μm s−1 for all scans. Each scan had 512 lines, with 512 data points per line, corresponding to pixel sizes in the range of 0.2 to 98 nm. The wear of the AFM tip was minimized using the best practices described in [42]. Specifically, the values of free-air amplitude and amplitude ratio, which is the ratio of the amplitude of AFM probe tip vibration when performing a scan to the amplitude when vibrating in free air, were kept in the range of 37–49 nm and 0.15–0.3, respectively. Though AFM provides a two-dimensional description of the surface topography, the data were analyzed as a series of line scans. This practice maintained consistency with the other techniques, which yield one-dimensional measurements, and also eliminated artifacts due to instrumental drift in the slow-scan axis.

The largest scales of topography were measured using a stylus profilometer (Alpha Step IQ, KLA Tencor, Milpitas, CA) with a 5-μm diamond tip. Measurements were collected at a scanning speed of 10 μm s−1, with data points every 100 nm. A total of 8 measurements were taken on each substrate, with 2 measurements each at scan sizes of 0.5, 1, 2, and 5 mm. All measurements were corrected using a parabolic fit to remove the tilt of the sample and the bowing artifact from the tool. For the UNCD and pUNCD, the larger scan sizes exhibited consistent non-parabolic trends due to instrument artifacts. These artifacts were corrected by taking reference scans on polished silicon wafers and subtracting the averaged reference profiles from the measurements.

Finally, the PSD was used to combine all measurements from a single surface into one averaged curve that describes the topography of that surface. The PSD is the Fourier transform of the autocorrelation function of a line scan with height  which is mathematically equivalent to the square of the amplitude of  i.e.,   is the wavevector and  is the length of the scan. All data were collected and analyzed as 1D line scans, enabling the calculation of the one-dimensional PSD, denoted here as C (designated  in [30, 32]). These calculations follow the standards established in [30] for computing and reporting PSDs.

3. Results

All topography data sets collected from this investigation are freely available for download and analysis [43–46].

3.1. Multi-scale topography measurement

Representative images of all three techniques are shown for the four materials in figure 2. The stylus profilometry data are shown with decreasing scan size in figure 2(a). It is clear that the roughness on the MCD has the largest amplitude (RMS height) while the pUNCD surface shows the smallest. Going from larger scans (figure 2(a1)) to smaller scans (figure 2(a3)), the amplitude of measured topography decreases for all four diamond species. For example, while a 5-mm scan of MCD spans a vertical range of 646 nm, a 0.5-mm scan of the same surface spans just 391 nm. In order to interpret the stylus data correctly, an estimate of the tip radius is needed because the tip introduces artifacts at and below this size scale. Figure 2(a4) shows a scanning electron microscopy (SEM) image of the stylus tip. Fitting a circle to the tip yields a radius of R = 5.1 μm. The exact point where tip artifacts become dominant [37] will be estimated using the PSD in the following section.

The AFM measurements are shown as topography maps of increasing resolution in figures 2(b)–(e). The NCD and MCD clearly show grain structure and faceting at scan sizes of 5 μm and 1 μm (figure 2(d1-2) and (e1-2)). The UNCD sample (figure 2(c1-2)) also shows strong texture, but the size indicates that these are multi-grain clusters. The pUNCD shows topographic features, but no grain structure. At the highest resolution (figures 2(b3)–(e3)), all features look relatively smooth. The smoothness of the features is likely related to tip artifacts (caused by convolution between tip curvature and topography). Figures 2(b4)–(e4) show TEM images of the AFM tips after they were used to image the surfaces. The tip radii were measured as R = 17 nm for pUNCD, R = 36, 47 and 31 nm for UNCD, NCD and MCD respectively using the same procedure used for the stylus tip. As with the stylus tip, a circle was fitted to the tip profile for the extraction of the radius (figures 2(b4), (c4), (d4), and (e4)). In cases where the tip apex did not appear perfectly circular, a best-fit circle was fitted to the region of the tip that makes contact with the substrate.

Finally, the surfaces were analyzed using side-view TEM (figures 2(f)–(i)). Once again, the pUNCD surfaces had the lowest amplitude of topography and the MCD had the highest. The NCD and MCD materials showed clear faceting from the individual crystallites. At the scales accessible by the TEM, the UNCD surface also shows faceting (see figures 2(g2), and (g3)). However, the smooth pUNCD surface shows no indication of faceting but rather a smoothly varying surface topography. Surprisingly, despite significant differences in topography at larger scales, the smallest-scale topography is nearly identical between the microcrystalline diamond, the nanocrystalline diamond, and the ultrananocrystalline diamond. This finding is further discussed in the next paragraph.

To further investigate the similarities in roughness at the smallest scales, representative images of the three unpolished surfaces are shown in greater detail in figure 3. In all cases, significant roughness is visible on the scale of Angstroms to nanometers. This small-scale topography is visible even on a single facet of a single grain of the NCD and MCD materials. The atomic lattice of the diamond is clearly visible in the TEM whenever a grain is aligned with a zone axis lying near to the imaging axis. In these cases, the lattice is observed in many areas to extend to within 1 nm of the surface. Therefore, it is not simply a rough and potentially amorphous surface layer that is sitting on a flat diamond facet, rather the diamond crystal itself exhibits significant roughness at the small scale.

3.2. Computing topography metrics, in real-space and in frequency-space

Now, a more quantitative analysis of the topography data is presented, and scalar roughness parameters are computed for the various surfaces. First, the root-mean-square (RMS) height  RMS slope  and RMS curvature  are computed in real-space from each line scan by numerically integrating the squared height data (or its derivatives) over the scan length  [32]:

using the trapezoidal rule (equation (4) of [32]). Figure 4 shows the computed roughness parameters as a function of scan size  (for RMS height, which depends on the larger-scale features) and measurement resolution  (for RMS slope and curvature, which depend on smaller-scale features). Note that for TEM data, the 'measurement resolution' is different from the size of a pixel in the camera ('pixel size') (which can be sub-atomic at the highest magnifications). The measurement resolution  is determined from the point spacing of the extracted profiles, as shown in figure 2. Computing the RMS height as a function of size is equivalent to an analysis of the surfaces' self-affine properties using a variable bandwidth method [47].

Figure 4(a) shows that  increases with  for all the surfaces studied here at small  i.e. L less than 1 μm. There is a crossover to constant (independent of )  at a scan size of 1–10 μm. All of the surfaces studied here show this transition, which corresponds to the thickness of these coatings (2 μm). The amplitude of the pUNCD surface is much smaller than the other three surfaces. NCD and MCD roll off to the same constant value. The values in figure 4(a) at large  correspond to the observation of the amplitude in stylus profilometry shown in figure 2, with pUNCD being the 'smoothest' and MCD being the 'roughest' surface at large scales.

From the analysis of ideal self-affine surfaces [22, 30, 32], it can be shown that  while  and  Therefore, in figure 4(a), the small-scale data has been fit with a power-law function (solid lines), and the extracted values are used to determine the expected trends in figures 4(b), (c) (dashed lines). As argued above and first described in [36, 48], it is apparent from figure 4 that no single value of RMS height, RMS slope, and RMS curvature can be defined, because these parameters depend on  or  This demonstrates a key difficulty that impedes efforts to link surface function to a single scalar roughness parameter; these common roughness parameters (including RMS roughness) are scanning-length dependent and do not describe an intrinsic property of the material.

As an alternative to analyzing the real-space measurements individually (as done in figure 4), the measurements can all be combined in frequency space to create a single PSD, denoted  that yields a complete statistical description of topography for each surface, as shown in figure 5. This analysis was carried as follows: first, the PSD of each measurement was computed, following the procedures laid out in [30, 32]. Second, a reliability cutoff [30, 32] due to tip artifacts was calculated for each PSD based on the measured tip radius (figure 2), and all data below this size scale was deemed unreliable and removed [32]. Third, all of the reliable portions of the many individual PSDs were combined by computing the arithmetic average of all measurements in logarithmically-spaced bins. The result is a single whole-surface PSD that describes the material across all length scales. There are no fitting parameters in this analysis; rather the PSD serves to separate the different size scales of topography, and the various techniques agree within experimental uncertainty. The only exception is for pUNCD, where the small-size stylus data lies below the AFM and TEM data, causing a dip around q = 107 m−1. This is believed to be an instrumental artifact, rather than resulting from the real topography. Overall, the value of these comprehensive PSDs is that they can be used in analytical and numerical models (such as [22, 25–27, 49–52]) to understand and predict surface properties (e.g. [17]).

To compute scale-invariant scalar roughness parameters, the full stitched-together PSD was used to compute RMS height, RMS slope, and RMS curvature as:

Table 2 shows the computed RMS parameters. Also computed from the whole-surface PSD is the area ratio, which is the increase in full surface area, per unit apparent area, due to the roughness. The area ratio is a critical parameter in investigations of soft-material adhesion, including [17], and it is calculated as described in [17]. While MCD is the roughest in terms of RMS height, the unpolished UNCD is the steepest in terms of RMS slope. However, while table 2 shows the mathematically correct or 'true' values of RMS parameters for a surface, any individual application or surface property may depend only on a certain range of length-scales. In that case, scale-dependent parameters would be recomputed by integrating only across the relevant size scales.

4. Discussion

4.1. Evaluating the fractal nature of diamond coatings, and the meaning of Hurst exponents

The Hurst exponent (which can be related to the fractal dimension, as described in [53]), can be calculated from the PSD, which is commonly separated (somewhat arbitrarily) into the 'self-affine' region, where the topography appears to be described by a power law relationship of  where  is the power-law exponent, and the 'roll-off' region, where the PSD appears to be flatter. The Hurst exponent  is typically extracted from the self-affine region as  [53, 54]. Using this procedure, and the region of the curves between  m−1 and  m−1, the resulting value was  and  for pUNCD and UNCD respectively; and  and  for NCD and MCD respectively.

There are two alternative methods of extracting the Hurst exponent: from the real-space data using the variable bandwidth method (VBM) [55]; and from the roll-off region by assuming that the full PSD is described by Fractional Gaussian Noise (FGN) [53], a hypothesis put forth by some of us in [32]. The VBM is nothing more than an analysis of the functional dependence of RMS height  as a function of scan size  (that scales as  for self-affine surfaces), as shown in figure 4(a). Note that even for an individual scan, the RMS height could be computed over a subsection of that scan (yielding an estimate of the Hurst exponent for a single realization of the topography [33, 47, 55]) but here full-size scans were used for calculation. As shown in figure 4(a), RMS height  can be accurately fit with a power-law form over the range from the smallest size ( ∼ 10 nm ) up to approximately  = 1 − 10 μm. Over this region, the Hurst exponents were  = 0.73 ± 0.18 for pUNCD;  = 0.74 ± 0.05 for UNCD;  = 0.93 ± 0.09 for NCD and  = 1.02 ± 0.10 for MCD. The Hurst exponents from this method () yield similar results to the above Hurst exponents () that are extracted from the more common method of fitting the 'self-affine' region of the PSD. The VBM is a useful technique because it is simpler to perform—requiring only a straightforward calculation of the root-mean-square height, calculated in real-space, from a series of multi-resolution measurements. Additionally, it can be readily performed on reentrant surfaces, while the calculation of the PSD requires surfaces to be functions (one height value for each lateral position).

The second alternative method for computing the Hurst exponent uses the 'roll-off' region and the fractional Gaussian noise (FGN) approach. Here, the Hurst exponent is given by  [32, 53], where  is the scaling exponent  at low  This results in  for pUNCD;  for UNCD; and  and 0.70 ± 0.05 for NCD and MCD, respectively. A prior paper by the present authors ([32]) speculated that there may be a connection between  and  (from the self-affine region). This observation would be extremely useful as it suggests that the small-scale behavior could be predicted from large-scale measurements. Unfortunately, when these four different surfaces are compared, there is no clear relationship that emerges.

Even in the traditional 'self-affine' portion of the PSD, a range of values can be extracted for Hurst exponent for a single surface depending on the window used for fitting  This is particularly true for MCD and NCD where Hurst exponents can be calculated by dividing the 'self-affine' region into two parts. For  in the range of  to  m−1, the extracted values are  = 1.27 ± 0.27 and 1.32 ± 0.09 for NCD and MCD, respectively; for q in the range of  to 1.8  m−1, the extracted Hurst exponents are  = 0.75 ± 0.04 and 0.78 ± 0.05 for NCD and MCD, respectively. Indeed, the roughness in these two portions of the curve seems to be qualitatively different, with the upper portion having a scaling behavior near  corresponding to  = 1.5. The origin of these differences in scaling behavior between different length scales is discussed in detail in the next section.

Because there can be so much variability in the measurement of a single surface, the whole practice of assuming self-affinity and assigning a single Hurst exponent to describe a surface must be done with caution. It is mathematically convenient to assume self-affinity as this simplifies numerical and analytical models, and it is common practice in experiments to use assumptions of self-affinity to extrapolate to small scales where the topography is not easily measured. However, at least for the diamond materials investigated here, the best-fit value for  depends on the region over which it is measured, and it is strongly influenced by other factors such as grain size (see next section). Instead, where possible, it is preferable to measure surface topography across all size scales and to use the whole-surface PSD as the primary descriptor for the surface, rather than any scalar parameter.

4.2. The effect of grain size on topography

As discussed in the previous section, the larger-grain-size materials (NCD and MCD) demonstrate a region where the PSD scaling is similar to  in the larger-wavelength portion of the 'self-affine' region (see solid black lines in figures 5(d) and (e)). This scaling is characteristic of 'kinks' in the real-space line scan, such as sharp peaks or valleys. Note that such kinks can arise in topography as an artifact of the nonvanishing tip radius [30, 37], so one must first rule out this as the cause. However, in the present work, this  scaling is clearly observed both in the reliable portion of the AFM measurement and in the TEM measurement, both of which are free from tip-based artifacts. Therefore, this  scaling in the PSDs of the MCD and NCD is a feature of the measured topography, rather than emerging from an artifact. This behavior of the PSD corresponds to kinks in the surface topography that are directly observable in the TEM imaging, as shown in figures 2 and 3.

Our hypothesis to describe this local  scaling is that it is characteristic of topography at scales approximately equal to the grain size of the material. At these scales, features are dominated by the crystal facets and kinks between them: adjacent grains with different orientations will create concave kinks where the grains meet, while edges between crystal facets will create convex kinks. It is therefore assumed that typical topography line scans will pass through approximately two kinks (one concave and one convex) per grain. At sizes much below the grain size, the topography of MCD and NCD reverts to random, self-affine behavior, showing roughness very similar to UNCD (see figure 3), where the PSD scales as  (see blue dashed lines in figures 5(b)–(e)). This finding agrees with prior work [56] on fracture surfaces in sandstone. In that work, stylus profilometry measurements showed a transition around the grain size; however sub-grain features were mostly inaccessible there due to the aforementioned tip artifacts. Computer-generated profiles were used to verify that facets can cause -scaling at a wavevector related to the average kink spacing, designated . The mathematical basis for this hypothesis linking kink spacing to  scaling in the PSD is given in the appendix.

In order to demonstrate how the PSD is affected by the superposition of facets from a characteristic grain size and random roughness below that size, artificial one-dimensional surfaces were created that were composed of: a superposition of triangular peaks (figure 6(a1)); self-affine random roughness (figure 6(a2)); and the summation of those two into a single surface (figure 6(a3)). The piecewise linear surface (figure 6(a1)) has kinks with uncorrelated heights drawn from a Gaussian distribution and lateral distances between kinks drawn from a Rayleigh distribution. The surface is scaled in order to have an RMS slope of 1 and an average kink spacing  The small-scale self-affine random roughness (figure 6(a2)) is generated using a Fourier-filtering algorithm [30, 57] with a Hurst exponent  and RMS-slope of 1.2.

Figure 6. Computer-generated profiles support the experimental finding of characteristic regions of topography scaling for polycrystalline surfaces. Panel (a) shows a representative section of the three types of computer-generated profiles investigated: a faceted profile with random, uncorrelated kink heights and spacings (a1), a self-affine random roughness (a2); and the sum of these profiles (a3). Their PSDs, averaged over logarithmically-spaced bins, are shown in (b). Three characteristic scaling regions (see main text) of the PSD for the summed surface are indicated by background color and labels above the graph. Panel (c) shows a surface created by summing a faceted profile with correlated heights and the small-scale self-affine random roughness (same as a2). The resulting PSD (d) shows four characteristic scaling regions. The RMS slope of the faceted profiles  has been used to normalize all PSDs and profiles.

When the PSDs (figure 6(b)) are computed from these surfaces, there is a transition from a flat PSD at the largest sizes to scaling as , and the transition point occurs at a wavelength of , which corresponds to . (The transition is somewhat gradual; the choice of this particular value for the transition point is motivated from the mathematical consideration in the appendix). Importantly, the summed-surface PSD follows the kinked-surface PSD and displays  scaling at larger scales, and then transitions to self-affine scaling at smaller scales, in this case corresponding to a Hurst exponent  of 0.8. To ensure that these results were not unique to the particular way that the kinked surface was generated, this analysis was repeated using kinked surfaces with uniform and exponential distributions of kink spacings, as well as a surface with slopes alternating between  and 1. These analyses are shown in the supplemental section 1 (available online at stacks.iop.org/STMP/9/014003/mmedia), but the results are similar, differing only in the sharpness of transitions. The key finding from this analysis is that the kink spacing, which corresponds approximately to the grain size, introduces a signature in the scaling of the topography that causes deviations from the commonly assumed fractal-like self-affine scaling.

The PSD of the summed surface (figure 6(b)) reproduces the three regions that are visible in the PSDs of MCD and NCD: flat behavior at small  (large sizes); scaling like  at intermediate values; and self-affine scaling () at large (small sizes). However, this does not yet explain the self-affine scaling behavior that is observed above the grain size in polished and unpolished UNCD. To account for this, another synthetic surface was created that is similar to the first, but this time with spatially correlated kink heights (figure 6(c)). The spatial correlation of the kink heights is enforced using the method of [56] as follows: a self-affine random surface was generated using the Fourier-filtering algorithm; then facets were computed by interpolating linearly between random points (the kinks); finally this faceted surface was summed with the same self-affine random roughness as before (figure 6(a2)). While this does not substantially alter the behavior at small wavelengths, this adds in an additional self-affine region above the grain size (figure 6(d)).

To verify this proposed link between  scaling of the PSD and the kink spacing of the material, an analysis of kink spacing and grain size was performed on the real materials using the AFM and TEM measurements. For the MCD and NCD, concave kinks between grains were readily visible, and therefore were quantified using feature sampling, as is used in metallographic analysis [58] (supplemental section 2.1). The mean lineal intercept between concave kinks, which is of the same order as the grain diameter [58], was measured as 839 ± 68 nm for MCD and 647 ± 42 nm for NCD. Because the mean kink spacing used in the mathematical analysis involves both convex and concave kinks, it was computed as half of this value, i.e. 419 ± 34 nm for MCD and 323 ± 21 nm for NCD. These values can be converted to frequency-space using  as discussed above, and the values calculated are  m−1 for MCD and  m−1 for NCD. The same approach could not be applied to UNCD where the kink spacing was at or below the reliability cut-off due to tip artifacts. Therefore, the mean kink spacing for UNCD was estimated as half of the grain size, which was computed by averaging a sample of 20 different grains observed in the TEM (see supplemental section 2.2). The mean grain size of the UNCD was determined to be 14 ± 3 nm, which corresponds to a kink spacing for UNCD of 7 ± 2 nm and  m−1. The mean kink spacings of these three materials are indicated as vertical bars in figures 5(c)–(e). Indeed, the MCD and NCD demonstrate that the scaling transitions to  right around the mean kink spacing. For UNCD, the mean kink spacing is too small to observe a clear  scaling regime; however, the self-affine scaling behavior is confirmed for sizes much larger than the characteristic grain size. The remainder of the paper will frame results in terms of grain size (again, approximately twice the value of kink spacing), since grain size is more widely measured and reported for polycrystalline materials.

In summary, based on the experimental measurements of these four polycrystalline diamond surfaces, and based on the computed PSDs of artificially generated surfaces, four characteristic regimes of topography scaling have been identified. (1) At the smallest size scales (if significantly smaller than the grain size) power-law scaling may be observed that is characteristic of self-affine random roughness, e.g.,  corresponding to H in the range of 0.6–0.9. (2) At size scales similar to and slightly smaller than the average grain size, the PSD displays characteristic scaling of  due to grain facets and kinks. (3) At sizes larger than the grain size, but smaller than the film thickness, there is another region of power-law scaling corresponding to random roughness (). (4) Finally, at sizes larger than the film thickness, the PSD flattens out, with scaling in the range of  to  These four regions of topography scaling accurately describe the polycrystalline diamond surfaces and computer-generated surfaces investigated here. These four regions and their boundaries may be broadly generalizable to other materials, but further investigation is required. For instance, the boundary between regions (3) and (4) corresponded to the thickness of the diamond coatings in these measurements, but this thickness was not varied to explicitly investigate this connection. The future application of comprehensive topography characterization on other materials will elucidate the applicability of these four regions to other polycrystalline materials.

5. Conclusions

First, these results further establish a multi-resolution approach that is designated 'comprehensive topography characterization,' which combines multiple different techniques at multiple magnifications for the same surface. Then the power spectral density can be used to combine all measurements into one statistical description of the surface. Because typical roughness metrics are inherently scale-dependent and incomplete, this paper provides a method to understand roughness at all scales, including the specific scale over which it may be relevant in a given device. Particularly the measurement of small-scale roughness may be extremely important to predict and tailor surface properties such as adhesion, friction and wear. Furthermore, all topography measurements from this publication have been made publicly available, [43–46] so that other experimentalists may compare results and so that computational modelers may use it in models to predict properties of diamond materials. The purpose of this data and the underlying approach is to advance the field towards the goal of fundamental, predictive understanding of the performance of rough surfaces.

Second, these results show that the surface roughness of polycrystalline diamond materials varies significantly with scale, with surfaces that are smoother when measured at the large-scale showing roughness that is identical or even higher when measured at the smaller scales. Furthermore, while self-affine scaling (    ) is observed over some length-scales, the grain size introduces a signature into the power spectral density, showing  scaling behavior at scales approximately equal to the grain size. All unpolished surfaces show identical self-affine scaling at the smallest scale. This is a signature of small-scale roughness that is superimposed on the crystalline facets and is not observable with conventional (AFM, stylus) measurement techniques. Altogether, four characteristic regions of topography scaling were observed; these are expected to be applicable to all unpolished polycrystalline diamond films and may apply more broadly to other polycrystalline materials.


Sample preparation was performed in the Fischione Instruments Electron Microscopy Sample Preparation Laboratory at the University of Pittsburgh. The authors acknowledge the use of two facilities at the University of Pittsburgh: the Materials Micro-Characterization Laboratory (MMCL) in the Department of Mechanical Engineering and Materials Science; and the Nanoscale Fabrication and Characterization Facility (NFCF) in the Petersen Institute of Nano Science and Engineering (PINSE). The authors thank A. Baker for supplemental TEM imaging. LP acknowledges funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy (project EXC-2193/1–390951807) and by the European Commission (ERC-StG-757343). Funding for TDBJ and the University of Pittsburgh work was provided by the National Science Foundation under award number CMMI-1727378.

Data availability statement

The data that support the findings of this study are openly available at the following URLs: https://contact.engineering/go/wcqj3/https://contact.engineering/go/8sc7t/https://contact.engineering/go/cjy6s/https://contact.engineering/go/mz7z5/.


  1. Al-Azizi A A, Eryilmaz O, Erdemir A and Kim S H 2014 N anotexture fora wear-resistant and near-frictionless diamond-like carbon Carbon 73 403-12

  2. Bernal R A, Chen P, Schall J D, Harrison J A, Jeng Y-R and Carpick R W 2018 Influence of chemical bonding on the variability of diamond-like carbon nanoscale adhesion Carbon 128267-76

  3. J’älu§, Bramowicz M, Kulesza S, Ghaderi A, Dalouji V, Solaymani S, Kenari M F and Ghoranneviss M 2016 Fractal features and surface micromorphology of diamond nanocrystals Microsc. 264143-52

  4. Schade A, Rosiwal S M and Singer R F 2007 Influence of surface topography of HF-CVD diamond films on self-mated planar sliding contacts in dry environments Coatings Technol. 2016197-205

  5. Kovalchenko A M, Elam J W, Erdemir A, Carlisle J A, Auciello O, LiberaJA,PellinMJ,GruenDMandHrynJN 2011 Development of ultrananocrystalline diamond (UNCD) coatings for multipurpose mechanical pump seals Wear270 325-31

  6. Jacobs T D B, Ryan К E, Keating P L, Grierson D S, Lefever J A, Turner К T, Harrison J A and Carpick R W 2013 The effect of atomic-scale roughness on the adhesion of nanoscale asperities: a combined simulation and experimental investigation Lett. 50 81-93

  7. Ryan К E, Keating P L, Jacobs T D B, Grierson D S, Turner К T, Carpick R W and Harrison J A 2014 Simulated adhesion between realistic hydrocarbon materials: Effects of composition, roughness, and contact point Langmuir 30 2028-37

  8. Alcaide M, Papaioannou S, Taylor A, Fekete L, Gurevich L, Zachar V and Pennisi С P 2016 Resistance to protein adsorption and adhesion of fibroblasts on nanocrystalline diamond films: The role of topography and boron doping Mater. Sd. Mater. Med. 27 90

  9. May P W 2000 Diamond thin films: A 21 st-century material Trans. R. Soc. London. Ser. А Math. Phys. Eng. Sei. 358 473-95

  10. Otschik J, Schrüfer A, Thelke J, Kirchhof M and Schmaderer S 2010 Diamond seal faces—the benefits 201018-20

  11. Auciello О etal Are diamonds a MEMS’ best friend? IEEE Microw. Mag. 8 2007 61-75

  12. NistorP AandMayP W2017 Diamond thin films: Giving biomedical applications a new shine R. Soc. Interface 14 20170382

  13. Zeng H, Jarvik R, Catausan G, Moldovan N and Carlisle J 2016 Diamond coated artificial cardiovascular devices Coatings Technol. 302 420-5

  14. Biener J et al 2009 Diamond spheres for inertial confinement fusion Fusion 49112001

  15. Konicek A R, Grierson D S, Gilbert PUPA, Sawyer W G, Sumant A V and Carpick R W 2008 Origin of ultralow friction and wear in ultrananocrystalline diamond Rev. Lett. 100 235502

  16. Pastewka L, Moser S, Gumbsch P and Moseler M 2011 Anisotropic mechanical amorphization drives wear in diamond Nat. 10 34-8

  17. Dalvi S, Gujrati A, Khanal S R, Pastewka L, Dhinojwala A and Jacobs T D В 2019 Linking energy loss in soft adhesion to surface roughness Natl Acad. Sei. 116 LP-25425490

  18. Greenwood J A and Williamson J В P1966 Contact of nominally flat surfaces R. Soc. London. Ser. A. Math. Phys. Sd. 295300-19

  19. Fuller KN GandTaborD 1975Theeffectofsurfaceroughness on the adhesion of elastic solids R. Soc. London. A. Math. Phys. Sei. 345 327-42

  20. Persson В N J and Tosatti E 2001 The effect of surface roughness on the adhesion of elastic solids Chem. Phys. 115 5597

  21. Persson BN] 2001 Theory of rubber friction and contact mechanics Chem. Phys. 115 3840-61

  22. Persson В N J, Albohr O, Tartaglino U, Volokitin AI and Tosatti E 2004 On the nature of surface roughness with application to contact mechanics, sealing, rubber friction and adhesion Phys. Condens. Matter 17 Rl-62

  23. Persson В N J 2001 Elastoplastic contact between randomly rough surfaces Rev. Lett. 87116101

  24. Hyun S, Pei L, Molinari J-F and Robbins M О 2004 Finite- element analysis of contact between elastic self-affine surfaces Rev. E 70 26117

  25. Campanä C and Müser M H 2007 Contact mechanics of real vs. randomly rough surfaces: A Green’s function molecular dynamics study Lett. 77 38005

  26. Pastewka L and Robbins M О 2014 Contact between rough surfaces and a criterion for macroscopic adhesion Proc Natl Acad. Sei. Ill EP-3293303

  27. Müser MHeta/2017 Meeting the contact-mechanics challenge Lett. 65118

  28. Prodanov N, Dapp W В and Müser M H 2014 On the contact area and mean gap of rough, elastic contacts: Dimensional analysis, numerical corrections, and reference data Lett. 53433-48

  29. Yastrebov V A, Anciaux G and Molinari J-F 2015 From infinitesimal to full contact between rough surfaces: Evolution of the contact area J. Solids Struct. 5283-102

  30. Jacobs T D B, Junge T and Pastewka L 2017 Quantitative characterization of surface topography using spectral analysis Topogr.:Metrol. Prop. 5 013001

  31. Gujrati A, Khanal S R and Jacobs TDB2017A method for quantitative real-time evaluation of measurement reliability when using atomic force microscopy-based metrology 2017 IEEE 17th Int. Conf. Nanotechnology, NANO 017135-8

  32. Gujrati A, Khanal S R, Pastewka L and Jacobs T D В 2018 Combining ТЕМ, AFM, and profilometry for quantitative topography characterization across all scales Mater. Interfaces 1029169-78

  33. Hinkle A R, Nöhring W G, Leute R, Junge T and Pastewka L 2020 The emergence of small-scale self-affine surface roughness from deformation Sei Adv. 6 eaax0847

  34. Candela T, Renard F, Klinger Y, Mair K, Schmittbuhl J and Brodsky E E 2012 Roughness of fault surfaces over nine decades of length scales Geophys. Res. Solid Earth 17 B08409

  35. Mandelbrot В B, PassojaD E and Paullay A J1984 Fractal character of fracture surfaces of metals Nature 308721-2

  36. Sayles R S and Thomas T R1978 Surface topography as a nonstationary random process Nature 271431-4

  37. Church E L and Takacs P Z1991 Effects of the nonvanishing tip size in mechanical profile measurements 34th Annual International Technical Symposium on Optical and Optoelectronic Applied Science and Engineering 1332 (San Diego, CA, 1 January 1991) 504-14

  38. Stover J C1995 Optical Scattering: Measurement and Analysis, Second, SPIE, 100020th Street (Bellingham, WA 98227- 0010 USA)

  39. Wolski M, Podsiadlo P, Stachowiak G W, Holmberg K, Laukkanen A, Ronkainen H, Gee M, Nunn J, Gachot C and Li L 2017 Multiscale characterisation of 3D surface topography of DLC coated and uncoated surfaces by directional blanket covering (DBC) method 388-389 47-56

  40. Khanal S R, Gujrati A, Vishnubhotla S B, Nowakowski P, Bonifacio C S, Pastewka Land Jacobs T D В 2018 Characterization of small-scale surface topography using transmission electron microscopy Topogr.: Metrol. Prop. 6 045004

  41. Zeng H, Konicek A R, Moldovan N, Mangolini F, Jacobs T D B, Wylie I, Arumugam P U, Siddiqui S, Carpick R W and Carlisle J A 2015 Boron-doped ultrananocrystalline diamond synthesized with an Н-rich/Ar-lean gas system Carbon 84 103-17

  42. Vahdat V and Carpick R W 2013 Practical method to limit tipsample contact stress and prevent wear in amplitude modulation atomic force microscopyACS 7 9836-50

  43. Gujrati A, Khanal S, Pastewka L and Jacobs T D В 2020 Polished Ultrananocrystalline Diamond (Version 1), Contact. Engineering, (https://contact.engineering/go/wcqj3/)

  44. Gujrati A, Khanal S R, Pastewka L and Jacobs T D В 2020 Ultrananocrystalline diamond (UNCD) (Version 1), Contact. Engineering. (https://contact.engineering/go/8sc7t/)

  45. Gujrati A, Khanal S, Pastewka L and Jacobs T D В 2020

  46. N anocrystalline Diamond (Version 1), ContactEngineering. (https://contact.engineering/go/cjy6s/)

  47. Gujrati A, Khanal S, Pastewka L and Jacobs T D В 2020 Microcrystalline Diamond (Version 1), Contact.Engineering. (https://contact.engineering/go/mz.7z5/)

  48. Mandelbrot В В1985 Self-affine fractals and fractal dimension Scr. 32 257-60

  49. Thomas T R and Sayles RS 1978 Some problems in the tribology of rough surfaces Int. 11163-8

  50. Persson В N J and Scaraggi M 2014 Theory of adhesion: Role of surface roughness J. Phys. 141124701

  51. Persson В N J 2002 Adhesion between elastic bodies with randomly rough surfaces Rev. Lett. 89 245502

  52. Bush A W, Gibson R D and Thomas T R1975 The elastic contact of a rough surface Wear 3587-111

  53. Putignano C, Afferrante L, Carbone G and Demelio G 2012 The influence of the statistical properties of self-affine surfaces in elastic contacts: A numerical investigation Meeh. Phys. Solids 60 973-82

  54. GneitingT and Schlather M 2004 Stochastic models that separate fractal dimension and the hurst effect SIAM Rev. 46 269-82

  55. OreyS 1970 Gaussian sample functions and the Hausdorff dimension of level crossings Zeitschrift Für Wahrscheinlichkeitstheorie Und Verwandte Gebiete 15 249-56

  56. Schmittbuhl J, Vilotte J-P and Roux s 1995 Reliability of self- affine measurements Phys. Rev. E 51131-47

  57. Boffa J M, Allain C and Hulin J P1998 Experimental analysis of fracture rugosity in granular and compact rocks Phys. J. AppLPhys. 2281-9

  58. Ramisetti S B, Campanä C, Anciaux G, Molinari J F, Müser M H and Robbins MO 2011 The autocorrelation function for island areas on self-affine surfaces Phys. Condens. Matter23 215004

  59. Russ J C and Dehoff RT 2000 Practical Stereology, Second (US, Boston, MA: Springer)