next up previous contents
Next: 5.1.1 Algorithmic Complexity Up: 5. Surface Rate Calculation Previous: 5. Surface Rate Calculation

5.1 Conventional Approach

One way to calculate the surface rates is by direct integration. This requires an appropriate discretization of the domain of the flux distribution function $ {\Gamma}={\Gamma}({\vec{x}};{q}, {\vec{\omega}}, {E})$ . Hence, the flux distribution at a certain surface point $ {\vec{x}}$ is approximated as a superposition of a finite number of appropriate basis functions $ {b}_1, {b}_2,\ldots, {b}_{{N}_{\Gamma}}$

$\displaystyle {\Gamma}({\vec{x}};{q}, {\vec{\omega}}, {E}) \approx \sum_{l=1}^{{{N}_{\Gamma}}} {B}_l({\vec{x}}) {b}_l({\vec{x}};{q}, {\vec{\omega}}, {E}).$ (5.1)

$ {B}_l({\vec{x}})$ are the corresponding coefficients. If the basis functions represent an orthonormal basis, thus

$\displaystyle \sum_{{q}=1}^{{Q}} \int\limits_{0}^{\infty} \int\limits_{{\vec{\o...
...}) {b}_k({\vec{x}};{q}, {\vec{\omega}}, {E}) \,{d}{\Omega}\,{d}{E}=\delta_{lk},$ (5.2)

the coefficients $ {B}_l({\vec{x}})$ can be calculated using

$\displaystyle {B}_l({\vec{x}}) = \sum_{{q}=1}^{{Q}} \int\limits_{0}^{\infty} \i...
...ega}}, {E}) {\Gamma}({\vec{x}};{q}, {\vec{\omega}}, {E}) \,{d}{\Omega}\,{d}{E}.$ (5.3)

Insertion of the ballistic transport equation (2.14) yields

$\displaystyle {B}_l({\vec{x}}) = S_{l}({\vec{x}}) + \int\limits_{{\mathcal{S}}}...
...) {\Gamma}_{\text{re}}({\vec{x}}';{q}, {\vec{\omega}}, {E}) \,{d}{E} \,{d}{A}',$ (5.4)

where $ S_l$ denotes the source term

$\displaystyle S_{l}({\vec{x}}) := \int\limits_{{\mathcal{P}}} \frac{-{\vec{\ome...
... {\Gamma}_{\text{src}}({\vec{x}}';{q}, {\vec{\omega}}, {E}) \,{d}{E} \,{d}{A}'.$ (5.5)

Here, $ \funcvis ({\vec{x}},{\vec{x}}')$ is the visibility function which returns 1 or 0, if the surface points $ {\vec{x}}$ and $ {\vec{x}}'$ are in line or sight or not, respectively. Utilizing expression (2.15) finally gives

$\displaystyle {B}_l({\vec{x}}) = S_l({\vec{x}}) + \int\limits_{{\mathcal{S}}} \sum_{k=1}^{{{N}_{\Gamma}}} T_{lk}({\vec{x}},{\vec{x}}'){B}_k({\vec{x}}')\,{d}{A}'$ (5.6)


$\displaystyle T_{lk}({\vec{x}}, {\vec{x}}'):= \frac{-{\vec{\omega}}\cdot{\vec{n...
...b}_k({\vec{x}}';{q}', {\vec{\omega}}', {E}') \,{d}{\Omega}' \,{d}{E}' \,{d}{E}.$ (5.7)

To solve the surface integral equation (5.6) numerically, a discretization of the surface is needed. Triangle [68] or voxel elements [6] are conventionally used in three dimensions. If $ {{N}_\text{S}}$ denotes the number of discretization elements of the surface, the number of unknowns is equal to $ {{N}_\text{S}}\cdot{{N}_{\Gamma}}$ . If $ {{N}_{\Gamma}}$ is large, thus if many basis functions are used to resolve the direction and energy dependence of the flux distribution, the number of variables, which must be kept in memory, is very large. This is especially a problem in three dimensions, where $ {{N}_\text{S}}$ is already very large by nature. Therefore, only a very small number of basis functions is commonly used in practice. Very often, only one basis function is used per particle type

$\displaystyle {b}_l({\vec{x}};{q}, {\vec{\omega}}, {E}):=\delta_{l{q}},$ (5.8)

which reduces the number of unknowns to $ {{N}_\text{S}}\cdot{Q}$ . As a consequence, the corresponding coefficients are equal to the respective total arriving flux of particles of species $ {q}$

$\displaystyle {B}_l({\vec{x}}) = \int\limits_{0}^{\infty} \int\limits_{{\vec{\o...
...({\vec{x}};{q},{\vec{\omega}}, {E}) \,{d}{\Omega}\,{d}{E} = {F}_{q}({\vec{x}}).$ (5.9)

The particle transport equation (2.14) simplifies to

\begin{displaymath}\begin{split}{F}_{q}({\vec{x}}) =& \int\limits_{\mathcal{P}} ...
...{x}}';{q}, {\vec{\omega}}, {E}) \,{d}{E} \,{d}{A}', \end{split}\end{displaymath} (5.10)

and the reemitted flux distribution (2.15) can be written as

$\displaystyle {\Gamma}_{\text{re}}({\vec{x}}';{q}, {\vec{\omega}},{E}) = \sum_{...
...}}';{q},{\vec{\omega}},{E};{q}',{\vec{\omega}}',{E}') \,{d}{E}' \,{d}{\Omega}'.$ (5.11)

Here, the incident direction and the energy do not play a role in the reemitted flux distribution $ {\Gamma}_{\text{re}}$ . The choice of the basis according to (5.8) assumes implicitly that the direction and the energy of reemitted particles is independent of the incident direction and energy. Therefore, effects such as specular-like reflexions of ions (see Section 2.2.3) cannot be accurately described, except for unidirectional and monoenergetic ion sources, for which the reemission probability function is the same for all incident ions.

Similarly, if the chosen basis is not able to reproduce the direction or the energy dependence of the arriving flux distribution $ {\Gamma}$ , the surface rates cannot be computed properly, given that they depend on the incident angle or energy. Hence, non-trivial surface rates, such as sputter rates, cannot be calculated correctly using the simple basis (5.8) and inserting (5.1) into (2.24). Only for unidirectional and monoenergetic particle incidence the sputter rates can be calculated correctly.

At least the contribution of primary particles, which come directly from the source plane $ {\mathcal {P}}$ , to the surface rates can be calculated correctly with reasonable effort using the simple basis (5.8). Once the coefficients $ {B}_l({\vec{x}})$ are calculated, the reemitted flux distribution can be approximated by inserting (5.1) into (2.15). Then the surface rates can be calculated using the formula

\begin{displaymath}\begin{split}{R}_l({\vec{x}}) =& \int\limits_{{\mathcal{P}}} ...
...vec{x}}';{q},{\vec{\omega}},{E}) \,{d}{E} \,{d}{A}' \end{split}\end{displaymath} (5.12)

which is obtained by combining (2.14) and (2.24).

next up previous contents
Next: 5.1.1 Algorithmic Complexity Up: 5. Surface Rate Calculation Previous: 5. Surface Rate Calculation

Otmar Ertl: Numerical Methods for Topography Simulation