The paper investigates a model of the photoconductivity of macroporous silicon in the conditions of homogeneous generation of photocarriers. By the finite element method, the stationary photoconductivity and the time evolution of photoconductivity after instantaneous shutdown of light are calculated. Dependences of the stationary photoconductivity and relaxation time of photoconductivity on the velocity of recombination of nonequilibrium carriers at the surfaces of pores, radius of pores, and average distance between them are analyzed.

1. Introduction

One of urgent today’s problems is the control of toxic and other harmful substances in air atmosphere, drinking water, foods, and so forth. In this connection, there is a necessity for development of microelectronic systems of detection of gases for industrial, office, and housing apartments which can be components of integral systems of protection of human’s life and health. Of special interest is the sensors based on nanoporous semiconductors, in particular porous silicon. Gas sensors of such type are described in a number of papers [1–4]. Among such sensors, one can distinguish the structures from porous silicon, based on the change of photoconductivity and its kinetics under the action of different gas environments which can be used as effective gas analyzers. For designing such sensors, it is necessary to know dependences of photoconductivity of porous silicon on the state of surface, depending on a gas environment where it is, as well as on the degree of porosity of a material and geometry of pores.

In this paper, for the case of macroporous silicon with cylindrical pores, modeling of stationary photoconductivity and process of its relaxation occurring after instantaneous shutdown of illumination is carried out.

2. Formulation of the Problem

Let us consider a plate of macroporous silicon with arranged parallel cylindrical pores of radius π‘Ÿ0, which are perodically distributed in the semiconductor volume with the average distance between their centers being equal to 2𝑅. These pores form a quadratic lattice in the π‘₯𝑂𝑦 plane lying perpendicular to the axis of the cylindrical pores. At calculating, we shall suppose the origin of coordinates is on axis of one of the pores. This semiconductor, for concreteness of 𝑝-type conductivity, is illuminated by light from the region of fundamental absorption for which the function of generation of photocarriers 𝐺, that is, the amount of photocarriers which are generated for one second in unit volume does not depend on the coordinate. It should be noted that the condition of homogeneous generation of photocarriers is satisfied in the case of weak absorption of light.

The generated photocarriers recombine both in the semiconductor volume and at the surfaces of pores. At constant illumination, it has been established such a nonuniform spatial distribution of photocarrier concentration at which balance between the processes of generation and recombination of nonequilibrium charge carriers is ensured.

In the case of diffusion transport of photocarriers, the dependence of the photocarrier concentration Δ𝑛(π‘₯,𝑦,𝑑) on the coordinates π‘₯, 𝑦, and time 𝑑 can be determined from the equationπœπ‘›πΏ2π‘›πœ•Ξ”π‘›=πœ•πœ•π‘‘2Ξ”π‘›πœ•π‘₯2+πœ•2Ξ”π‘›πœ•π‘¦2βˆ’Ξ”π‘›πΏ2𝑛+πœπ‘›πΊπ‘“(𝑑)𝐿2𝑛,(1) where 𝐿𝑛  is the diffusion length of electrons, πœπ‘› is their lifetime, and 𝑓(𝑑) is the function of photocarrier generation normalized on its maximum value 𝐺.

Equation (1) should be complemented by the boundary conditions1√π‘₯2+𝑦2ξ‚Έπ‘₯πœ•Ξ”π‘›πœ•π‘₯+π‘¦πœ•Ξ”π‘›ξ‚Ή||||πœ•π‘¦π‘₯2+𝑦2=π‘Ÿ20=π‘†πœπ‘›πΏ2𝑛||||Δ𝑛π‘₯2+𝑦2=π‘Ÿ20,(2)πœ•Ξ”π‘›|||πœ•π‘₯π‘₯=±𝑅=0,(3)πœ•Ξ”π‘›||||πœ•π‘¦π‘¦=±𝑅=0.(4) Equation (2) corresponds to recombination of photocarriers on the surface of cylindrical pore with the surface recombination velocity 𝑆 which depends on the state of surface of pores, in particular, on the kind and of the adsorbed molecules of a gas. Equations (3) and (4) follow from the symmetry conditions of the problem which give rise to the fact that, at the faces of β€œelementary cell” (π‘₯=±𝑅,𝑦=±𝑅), the photocarrier concentration reaches its maximum value.

The initial condition for (1) can have different forms, particularly, in the case of instantaneous shutdown of illumination, it takes the following form:||Δ𝑛(π‘₯,𝑦,𝑑)𝑑=0=Δ𝑛0(π‘₯,𝑦),(5) where Δ𝑛0(π‘₯,𝑦) is the stationary solution of (1).

The total number of photocarriers 𝑁 in the plate per its unit thickness can be calculated as𝑁=Δ𝑛(π‘₯,𝑦,𝑑)𝑑π‘₯𝑑𝑦,(6) where integration is over the region with area A represented as a square with side 2𝑅 with a central circular cutout of radius π‘Ÿ0. It should be noted that the total number of photocarriers is proportional to the photoconductivity of porous silicon.

For convenience, we further pass to the following dimensionless values:π‘₯βˆ—=π‘₯𝐿𝑛,π‘¦βˆ—=𝑦𝐿𝑛,π‘‘βˆ—=π‘‘πœπ‘›,π‘Ÿβˆ—0=π‘Ÿ0𝐿𝑛,π‘…βˆ—=𝑅𝐿𝑛,π‘†βˆ—=π‘†πœπ‘›πΏπ‘›,π‘βˆ—=π‘πΊπœπ‘›π΄.(7)

3. Calculation Results and Discussion

3.1. Stationary Photoconductivity

In the stationary case, the left-hand side of (1) is equal to zero and, on the right-hand side, 𝑓(𝑑) = 1. The transformed equations (1)–(5) were solved numerically by the finite element method [5] which is especially effective in the calculation of systems with complicated geometric configuration. The calculated dependences of the total number of photocarriers on the radius of pores and the ratio of distances between pores to their radius for different values of the surface recombination velocity are shown in Figure 1.

It is seen from Figure 1 the photoconductivity of porous semiconductor at fixed nonzero value of the velocity of surface recombination decreases with increasing the pore’s radius and increases with rising the average distance between them. Increase of the velocity of surface recombination at the fixed values of pore’s radius is accompanied by decreasing the photoconductivity which is saturated at very large values of the surface recombination velocity.

Along with the numerical modeling, we have carried out analytical calculation for the case of cylindrical symmetry of the problem which takes place when 𝑅≫𝐿𝑛. In this case with using the independent variable π‘Ÿ=√π‘₯2+𝑦2, (1)–(4) reduce to the following ones:1π‘Ÿπ‘‘ξ‚€π‘Ÿπ‘‘π‘Ÿπ‘‘Ξ”π‘›ξ‚βˆ’π‘‘π‘ŸΞ”π‘›πΏ2𝑛=βˆ’πΊπœπ‘›πΏ2𝑛,ξ€·π‘Ÿπ‘‘Ξ”π‘›0ξ€Έ=π‘‘π‘Ÿπ‘†πœπ‘›πΏ2π‘›ξ€·π‘ŸΞ”π‘›0ξ€Έ,𝑑Δ𝑛(𝑅)π‘‘π‘Ÿ=0,(8) whose solution via the dimensionless values defined by expressions (7) isξ€·π‘ŸΞ”π‘›βˆ—ξ€ΈπΊπœπ‘›π‘†=1βˆ’βˆ—ξ€ΊπΌ1ξ€·π‘…βˆ—ξ€ΈπΎ0ξ€·π‘Ÿβˆ—ξ€Έ+𝐾1ξ€·π‘…βˆ—ξ€ΈπΌ0ξ€·π‘Ÿβˆ—ξ€Έξ€»πΌ1(π‘…βˆ—)𝐾1ξ€·π‘Ÿβˆ—0ξ€Έ+π‘†βˆ—πΎ0ξ€·π‘Ÿβˆ—0ξ€Έξ€»βˆ’πΎ1(π‘…βˆ—)𝐼1ξ€·π‘Ÿβˆ—0ξ€Έβˆ’π‘†βˆ—πΌ0ξ€·π‘Ÿβˆ—0,ξ€Έξ€»(9)where 𝐼𝑛, 𝐾𝑛 are the Bessel’s functions of imaginary argument [6].

Comparison of the results of numerical and analytical calculations shows that, at >2𝐿𝑛, the latter coincide with an error less than 1%. When approaching 𝑅 to 0.5 𝐿𝑛 the error increases to 30%, but using the correction coefficients introduced by us in [7], the error amounts to 10–15%.

3.2. Nonstationary Case of Photoconductivity

We shall analyze the nonstationary photoconductivity for the case of instantaneous shutdown of light falling on the semiconductor plate, that is, when the function 𝑓(𝑑) featured on the right-hand side of (1) has the following form:𝑓(𝑑)=𝐻(βˆ’π‘‘),(10) where 𝐻(𝑑) is the Heaviside’s function whose value is zero for negative argument and one for positive argument.

The total number of photocarriers π‘βˆ— was calculated by the finite element method for different values of surface recombination velocity and porous silicon geometrical parameters. It has been obtained that the increase of π‘†βˆ— leads to the decrease of both the stationary value π‘βˆ— and the time of photoconductivity relaxation, with at lower values of the distance between the pores these changes being more significant.

The calculated dependences π‘βˆ—(𝑑) allow us to determine the relaxation time 𝜏rel as a slope of these dependences presented in the semilogarithmic scale. The obtained results presented in Figures 2(a) and 2(b) show that 𝜏rel nonlinearly decreases with increasing π‘†βˆ— and the most significant decrease of 𝜏rel takes place when π‘†βˆ— value is of the order of one. It turns out that, at the expense of changing the surface recombination velocity, one can reach decreasing the relaxation time by almost one order of magnitude. The latter occurs in the materials with high level of porosity when the area of pore’s surface per unit volume is large. The values of variation of the relaxation time realized in porous silicon at the surface recombination velocity exceeding the velocity of photocarrier diffusive motion (𝐿𝑛/πœπ‘›) are quite acceptable to be registered by simple measurement techniques.

4. Conclusions

In this paper, we have calculated by the finite element method the stationary photoconductivity and the time of photoconductivity relaxation for the case of macroporous silicon with cylindrical pores. It has been shown that both stationary photoconductivity and time of photoconductivity relaxation significantly depend on the surface recombination velocity of photocarriers, radius of pores, and average distance between pores what can be used for manufacturing gas sensors. Photoconductivity of macroporous silicon at fixed nonzero value of the velocity of surface recombination decreases with increasing the pore’s radius and increases with rising the average distance between them. Photoconductivity relaxation time nonlinearly decreases with increasing the velocity of surface recombination, and the most significant decrease of relaxation time takes place when value of the velocity of surface recombination is of the order of one.