## 1.

## Introduction

In both stochastic lithography modeling and analysis of roughness metrology data, it is sometimes necessary and often desirable to have an analytical expression for the power spectral density (PSD) that is both grounded in the known physics of stochastic processes and matches experimental evidence for those processes. Many diverse stochastic processes with a single correlation mechanism are known to follow an exponentially decaying autocorrelation function,
*R*:

## 1

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} R(r) = \sigma ^2 e^{ - (r/\xi)^{2\alpha } }, \end{equation}\end{document} $$R\left(r\right)={\sigma}^{2}{e}^{-{(r/\xi )}^{2\alpha}},$$*σ*is the standard deviation of the behavior (for example, the rms line-edge roughness),

*ξ*is the correlation length,

*α*is the roughness exponent, and

*r*is the distance (equal to |

*x*| in one dimension, [TeX:] $\sqrt {x^2 + y^2 }$ $\sqrt{{x}^{2}+{y}^{2}}$ in two dimensions, and [TeX:] $\sqrt {x^2 + y^2 + z^2 }$ $\sqrt{{x}^{2}+{y}^{2}+{z}^{2}}$ in three dimensions). This function is often called a stretched exponential or the Kohlrausch–William–Watts function

^{1}and is frequently encountered in relaxation processes,

^{2}as well as correlated roughness, where self-affine behavior of the roughness exists on length scales less than the correlation length. It has been used successfully in many studies of line-edge roughness, for example.

^{3, 4}

Since the correlation function of Eq. 1 is frequently encountered in stochastic processes, it makes sense to use this function as the basis for PSD analysis. The power spectral density is simply the Fourier transform of the correlation function (by the Wiener–Khinchin theorem). Unfortunately, analytical solutions to this Fourier transform are possible only for certain values of the roughness coefficient:
*α* = 0.5 and
*α* = 1. In Secs.
2, 3, 4, the analytical forms of the PSD for these two values of
*α* will be derived in one, two, and three dimensions. The PSD for
*α* = 0.5 in one and two dimensions has been previously derived,^{5} as has the
*α* = 1 case in one and two dimensions. Since most experimental line-edge roughness (LER) results show values of the roughness exponent between 0.5 and 1, these represent important limiting cases.

Analytic forms for the PSD are useful in two applications: metrology and simulation. Metrology data for LER and linewidth roughness can be fit by an analytical one-dimensional (1D) PSD, and surface roughness by a two-dimensional (2D) PSD, enabling the extraction of both
*σ* and
*ξ*. When generating random rough lines, surfaces, or volumes for simulation, an analytical form of the PSD can be used to generate random data with the desired autocorrelation response. Thus, it is useful to have as complete a set of analytical PSD functional forms as possible.

## 2.

## One-Dimensional Case

Since the autocorrelation function being used here is even, the Fourier transform in one dimension becomes a Fourier cosine transform.

## 2

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} G(f) = 2\int_0^\infty {g(x)\cos \left( {2\pi \,f\,x} \right)dx}. \end{equation}\end{document} $$G\left(f\right)=2{\int}_{0}^{\infty}g\left(x\right)\mathrm{cos}\left(2\pi \phantom{\rule{0.16em}{0ex}}f\phantom{\rule{0.16em}{0ex}}x\right)dx.$$For
*α* = 0.5:

## 3

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \frac{{2\sigma ^2 \xi }}{{1 + \left( {2\pi \,f\xi } \right)^2 }}. \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)=\frac{2{\sigma}^{2}\xi}{1+{\left(2\pi \phantom{\rule{0.16em}{0ex}}f\xi \right)}^{2}}.$$For
*α* = 1:

## 4

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \sqrt \pi \sigma ^2 \xi \,e^{ - \left( {\pi \,f\xi } \right)^2 }. \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)=\sqrt{\pi}{\sigma}^{2}\xi \phantom{\rule{0.16em}{0ex}}{e}^{-{\left(\pi \phantom{\rule{0.16em}{0ex}}f\xi \right)}^{2}}.$$## 3.

## Two-Dimensional Case

In two dimensions, the radial symmetry of the autocorrelation functions lends itself well to a Fourier transform using polar coordinates. The 2D Fourier transform in Cartesian coordinates is

## 5

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} G(f_x,f_y) = \int_{ - \infty }^\infty {\int_{ - \infty }^\infty {g(x,y)e^{ - i2\pi \left( {f_x x + f_y y} \right)} dxdy} }. \end{equation}\end{document} $$G({f}_{x},{f}_{y})={\int}_{-\infty}^{\infty}{\int}_{-\infty}^{\infty}g(x,y){e}^{-i2\pi \left({f}_{x}x+{f}_{y}y\right)}dxdy.$$## 6

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} x = r\cos \theta, \ \ y = r\sin \theta,\ \ f_x = f_r \cos \varphi,\ \ f_y = f_r \sin \varphi,\nonumber\!\!\! \\ \end{eqnarray}\end{document} $$\begin{array}{c}\hfill x=r\mathrm{cos}\theta ,\phantom{\rule{0.33em}{0ex}}\phantom{\rule{0.33em}{0ex}}y=r\mathrm{sin}\theta ,\phantom{\rule{0.33em}{0ex}}\phantom{\rule{0.33em}{0ex}}{f}_{x}={f}_{r}\mathrm{cos}\varphi ,\phantom{\rule{0.33em}{0ex}}\phantom{\rule{0.33em}{0ex}}{f}_{y}={f}_{r}\mathrm{sin}\varphi ,\end{array}$$## 7

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} G(f_r,\varphi) = \int_0^\infty {\int_0^{2\pi } {g(r,\theta)e^{ - i2\pi \,f_r r\cos (\theta - \varphi)} rdrd\theta } }. \end{equation}\end{document} $$G({f}_{r},\varphi )={\int}_{0}^{\infty}{\int}_{0}^{2\pi}g(r,\theta ){e}^{-i2\pi \phantom{\rule{0.16em}{0ex}}{f}_{r}r\mathrm{cos}(\theta -\varphi )}rdrd\theta .$$For the case of a radially symmetric function, the
*θ* integration can be carried out giving the Hankel transform (also called the Fourier–Bessel transform):

## 8

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} G(f_r) = 2\pi \int_0^\infty {r\,g(r)J_0 \left( {2\pi \,f_r r} \right)dr}, \end{equation}\end{document} $$G\left({f}_{r}\right)=2\pi {\int}_{0}^{\infty}r\phantom{\rule{0.16em}{0ex}}g\left(r\right){J}_{0}\left(2\pi \phantom{\rule{0.16em}{0ex}}{f}_{r}r\right)dr,$$*J*

_{0}is the Bessel function of the first kind, zero order. Applying the Hankel transform of Eq. 8 to the autocorrelation function of Eq. 1 results in the 2D PSD.

For
*α* = 0.5:

## 9

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \frac{{2\pi \sigma ^2 \xi ^2 }}{{[ {1 + \left( {2\pi \,f\xi } \right)^2 } ]^{{3}/{2}} }}. \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)=\frac{2\pi {\sigma}^{2}{\xi}^{2}}{{\left[1+{\left(2\pi \phantom{\rule{0.16em}{0ex}}f\xi \right)}^{2}\right]}^{3/2}}.$$For
*α* = 1:

## 10

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \pi \sigma ^2 \xi ^2 \,e^{ - \left( {\pi \,f\xi } \right)^2 }. \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)=\pi {\sigma}^{2}{\xi}^{2}\phantom{\rule{0.16em}{0ex}}{e}^{-{\left(\pi \phantom{\rule{0.16em}{0ex}}f\xi \right)}^{2}}.$$## 4.

## Three-Dimensional Case

In three dimensions, the 3D Fourier transform can be converted to spherical coordinates. For the special case of a radially symmetric function, the 3D Fourier transform in spherical coordinates becomes

## 11

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{eqnarray} G(f_r) &=& 4\pi \int_0^\infty {r^2 g(r)\left[ {\frac{{\sin \left( {2\pi \,f_r r} \right)}}{{2\pi \,f_r r}}} \right]dr}\nonumber\\ &=& \frac{2}{{f_r }}\int_0^\infty {rg(r)\sin \left( {2\pi \,f_r r} \right)dr}. \end{eqnarray}\end{document} $$\begin{array}{ccc}\hfill G\left({f}_{r}\right)& =& 4\pi {\int}_{0}^{\infty}{r}^{2}g\left(r\right)\left[\frac{\mathrm{sin}\left(2\pi \phantom{\rule{0.16em}{0ex}}{f}_{r}r\right)}{2\pi \phantom{\rule{0.16em}{0ex}}{f}_{r}r}\right]dr\hfill \\ & =& \frac{2}{{f}_{r}}{\int}_{0}^{\infty}rg\left(r\right)\mathrm{sin}\left(2\pi \phantom{\rule{0.16em}{0ex}}{f}_{r}r\right)dr.\hfill \end{array}$$For
*α* = 0.5:

## 12

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \frac{{8\pi \sigma ^2 \xi ^3 }}{{[ {1 + \left( {2\pi \,f\xi } \right)^2 } ]^2 }}. \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)=\frac{8\pi {\sigma}^{2}{\xi}^{3}}{{\left[1+{\left(2\pi \phantom{\rule{0.16em}{0ex}}f\xi \right)}^{2}\right]}^{2}}.$$For
*α* = 1:

## 13

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \pi ^{{3}/{2}} \sigma ^2 \xi ^3 \,e^{ - \left( {\pi \,f\xi } \right)^2 }. \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)={\pi}^{3/2}{\sigma}^{2}{\xi}^{3}\phantom{\rule{0.16em}{0ex}}{e}^{-{\left(\pi \phantom{\rule{0.16em}{0ex}}f\xi \right)}^{2}}.$$## 5.

## Summary

Letting
*d* be the dimensionality of the problem, the results can be summarized as follows:

For
*α* = 0.5:

## 14

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \frac{{a_d \,\sigma ^2 \xi ^d }}{{[ {1 + \left( {2\pi \,f\xi } \right)^2 } ]^{(d + 1)/2} }}, \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)=\frac{{a}_{d}\phantom{\rule{0.16em}{0ex}}{\sigma}^{2}{\xi}^{d}}{{\left[1+{\left(2\pi \phantom{\rule{0.16em}{0ex}}f\xi \right)}^{2}\right]}^{(d+1)/2}},$$*a*

_{1}= 2,

*a*

_{2}= 2

*π*, and

*a*

_{3}= 8

*π*.

For
*α* = 1:

## 15

[TeX:] \documentclass[12pt]{minimal}\begin{document}\begin{equation} {\rm PSD}(f) = \pi ^{d/2} \sigma ^2 \xi ^d \,e^{ - \left( {\pi f\xi } \right)^2 }. \end{equation}\end{document} $$\mathrm{PSD}\left(f\right)={\pi}^{d/2}{\sigma}^{2}{\xi}^{d}\phantom{\rule{0.16em}{0ex}}{e}^{-{\left(\pi f\xi \right)}^{2}}.$$## References

G. Williams and D. Watts, “Non-symmetrical dielectric relaxation behavior arising from a simple empirical decay function,” Trans. Faraday Soc. 66, 80–85 (1970).10.1039/tf9706600080Google Scholar

R. Böhmer, K. Ngai, C. Angell, and D. Plazek, “Nonexponential relaxations in strong and fragile glass formers,” J. Chem. Phys. 99(5), 4201–4209 (1993).10.1063/1.466117Google Scholar

V. Constantoudis, G. Patsis, L. Leunissen, and E. Gogolides, “Line edge roughness and critical dimension variation: Fractal characterization and comparison using model functions,” J. Vac. Sci. Technol. B22(4), 1974–1981 (2004).10.1116/1.1776561Google Scholar

P. Naulleau and J. Cain, “Experimental and model-based study of the robustness of line-edge roughness metric extraction in the presence of noise,” J. Vac. Sci. Technol. B25(5), 1647–1657 (2007).10.1116/1.2778697Google Scholar

P. Williams, D. Lubnau, and D. Elliott, “Generation of fractal images and comparison of their PSDs with several models,” Twenty-Second Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, pp. 495–499 (1988).Google Scholar