#### Abstract

We analyzed redistribution of radiation defects in a multilayer structure with porous epitaxial layer. The radiation defects were generated during radiation processing. It has been shown that porosity of epitaxial layer gives a possibility to decrease quantity of radiation defects.

#### 1. Introduction

One of actual questions of solid state electronics is increasing their radiation resistance. The increasing gives a possibility to decrease influence of different types of irradiation on characteristics of solid state electronic devices. Several methods have been used to increase the radiation resistance of devices of solid state electronics [1–5]. One way to solve the problem is using an epitaxial layer over considered devices (overlayer) to organize the resistance. Another way to increase the radiation resistance is using special epitaxial layers under considered devices to use these layers as drain of radiation defects. An alternative approach to the considered one is using porous materials near device area to use the porosity as drain of radiation defects again. In this paper, we consider an approach to decreasing quantity of radiation defects. The radiation defects have been generated during radiation processing of materials. For the framework the approach, we consider a heterostructure. The heterostructure consists of a substrate and porous epitaxial layer (see Figure 1). We assume that the substrate was under the influence of radiation processing (ion implantation, cosmic radiation etc.) through the epitaxial layer. Radiation processing of materials leads to generation of radiation defects. Main aim of the present paper is analysis of influence of porosity of epitaxial layer on distribution of concentration of radiation defects in the considered heterostructure.

#### 2. Method of Solution

To solve our aim, we calculate distributions of concentrations of radiation defects in considered heterostructure in space and time. We determine the above distributions as solutions of the following system of equations [6–12]:Here and are the distributions of concentrations of radiation interstitials and vacancies in space and time, respectively. The first, the second, and the third terms of both equations describe diffusion of point defects with the diffusion coefficients for interstitials and and vacancies, respectively. The fourth terms of both equations describe recombination of point defects with the parameter of recombination . The fifth, sixth, and seventh terms of both equations describe correction to diffusion due to porosity of material. The functions and describe dependencies of diffusion coefficients of defects due to porosity of materials on coordinate and temperature ; is the Boltzmann constant; is the molar volume; [10] is the chemical potential; are the initial and final volume of pores, respectively; J/(mole·K) is the molar gas constant. Last terms of (1) with nonlinearity of concentrations of defects and correspond to generation of divacancies and analogous complexes of interstitials (see, e.g., [11] and appropriate references in this work). The functions and describe dependencies of the parameters of generation of complexes point defects.

Boundary and initial conditions for (1) could be written asHere and are the equilibrium distributions of concentrations of interstitials and vacancies, respectively; , is the atomic spacing; is the specific surface energy. The above boundary conditions correspond to absence of flow of point defects through external boundary of heterostructure and absorption of these defects by pores (last condition). The above initial conditions correspond to distributions of concentration of the above defects after finishing radiation processing. To take into account porosity, we assume, that porous are approximately cylindrical with average dimensions and [13]. With time, small pores decompose into vacancies, and the vacancies are absorbed by large pores [10]. The large pores take spherical form during the absorption [10]. Distribution of concentration of vacancies, which was formed due to porosity, could be determined by summing all pores; that is,Here , , and are averaged distances between centers of pores in , , and directions, respectively; , , and are quantities of pores in the same directions.

We determine distributions of concentrations of divacancies and diinterstitials in space and time by solving the following system of equations [9, 11, 12, 14]:

The first, the second, and the third terms of both equations describe diffusion of point defects with the diffusion coefficients for diinterstitials and for divacancies. The fourth terms of both equations correspond to generation of new diinterstitials and divacancies. The fifth terms of the above equations correspond to decay of existing diinterstitials and divacancies. The functions and describe the parameters of decay of the above complexes on coordinate and temperature. The last terms of both equations describe correction to diffusion due to porosity of material. The functions and describe dependencies of diffusion coefficients of defects due to porosity materials on coordinate and temperature. Boundary and initial conditions for (4) could be written asThe above boundary conditions correspond to absence of flow of point defects through external boundary of heterostructure. The above initial conditions correspond to distributions of concentration of the above defects after finishing radiation processing.

We determine distributions of concentrations of radiation defects in space and time by method of averaging of function corrections [14–16]. To use the approach, we write (1) and (4) on account of initial distributions of defects; that is,

Farther, we replace the required concentrations in right sides of and on their not yet known average values . The replacement gives us possibility to obtain the following equations for determining the first-order approximations of concentrations of radiation defects in the following form:

Integration of the left and right sides of and on time gives a possibility to obtain the first-order approximations of concentrations of radiation defects in the final form:

Average values of the first-order approximations of the required approximations could be calculated by the following relation [14–16]:

Substitution of relations and into relation (6) gives a possibility to calculate the required average values in the following form:Here,

We determine approximations with the second and higher orders of concentrations of radiations defects framework standard iterative procedure of method of averaging of function corrections [14–16]. For the framework of the procedure, we determine the approximation of th order by replacement of the concentrations of radiation defects , , , and in right sides of and on the following sums . The replacement gives a possibility to obtain the second-order approximations of concentrations of radiation defects:

Integration of left and right sides of and gives a possibility to obtain relations for the second-order approximations of the required concentrations of radiation defects in the following form:

We determine average values of the second-order approximations by the standard relation [14–16]:

Substitution of relations and in relation (9) gives a possibility to obtain relations for the required values : where