Quasi-Chemical Reactions in Irradiated Silicon Crystals with Regard to Ultrafast Irradiation

This paper reports results from an investigation of the interaction of displaced Si-self atoms (I) and their vacancies (V), with impurities in crystalline silicon (Si), as induced by micro-second pulse duration irradiation with electrons at different energies: 3.5, 14, 25 and 50 MeV and pico-second pulse duration with energy 3.5 MeV. V-V, I-impurity atom and V-impurity atom interactions are analyzed both experimentally and as modeled using computer simulations. A process of divacancy (V2) accumulation in the dose-dependent linear region is investigated. The effect of impurities on recombination of correlated divacancies, and I-atoms that had become displaced from regular lattice points is estimated by computer modeling of an appropriate diffusion-controlled process. It is concluded that the experimental results can be interpreted quantitatively in terms of a strongly anisotropic quasi-one-dimensional diffusion of displaced I-atoms. In addition, a significant difference is found between the effects of pico-second duration electron beam irradiation, which causes the formation of A-centre (V + Oxygen) clusters, while when the beam is applied on a micro-second timescale, divacancies are created instead, although the electrons have the same energy in both cases.


Introduction
The formation of structural defects as a result of applying radiation to silicon crystals (Si) depends both on their impurity content and the lattice homogeneity.Although such phenomena have been the subject of numerous previous studies (for example [1]- [6]), there is a lack of detail available regarding the kinetics and accumulation rates for the interactions of impurity atoms with interstitial Si-self atoms (I) and their vacancies (V), i.e.Frenkel pairs, which are induced by the radiation in the Si crystal, in particular, by means of ultrafast radiation-pulses.
Information of this kind is of interest from both the fundamental and applied-scientific points of view because it enables conclusions to be drawn regarding the formation of secondary radiation defects.The latter are stable at room temperature), and are the termed the "end-product" of the radiation effect on the Si crystal, which alter the properties of the crystalline material itself, and accordingly those of any devices that are fabricated from it.
Both experimental and theoretical investigations show that, with the exception of GaAs, in general, the effect of the irradiation on the semiconductor crystal, is to induce a high mobility of interstitial atoms and their vacancies (even at temperatures below that of liquid nitrogen), in some cases to the degree that a practically zero-activation energy is observed [7] [8].As they move, these species undergo reactions with impurities and various crystal imperfections (which exist in the crystal prior to radiation treatment)-such as dislocations and pores, to form new complex radiation defects, which are stable at particular temperatures at which the properties of the crystal are defined.The type of new defect may be a divacancy (V 2 ), or complexes that comprise both initially electrically passive (oxygen, carbon) and active impurities (phosphorus, arsenic) in n-Si.The "electrical activation" of passive impurities in Si-samples via irradiation has been demonstrated [4] [9] [10] [11], and the challenge remains to elucidate the nature and distribution of these defects throughout the crystal lattice and to correlate this information with both the energy and dose of the applied radiation, i.e., the so-named "quasi-chemical reaction" which takes place within the irradiated crystal.
In the present work, the interaction of primary radiation defects with the crystal imperfections was not considered, because the samples used were initially almost identical to one another.Furthermore, we may note that the most readily displaced interstitial atoms in Si are typically impurity atoms from substituted lattice sites, which become mobile and take part at further reactions.In this study, our particular interest is in the formation of divacancies (V 2 ), and their rates of introduction and accumulation, since there are a number of different processes by which they can be created during the irradiation process, as "end-product" species, and which depend on the energy of the particular type of radiation being used (in the present work, electrons), which determines the final properties of the silicon crystal.
The main focus of this paper is given to the transfer of energy: from the initial radiation-crystal atom, atom to atom, vacancy-impurity, the formation and accumulation of divacancies, and the rate of production of radiation defects in Si, and the influence of these various factors on the electro-physical and optical properties of the n-Si crystals.The experimental measurements involved the H. N. Yeritsyan et al. determination of the concentration and mobility of the free charge carriers, along with infrared (IR) absorption spectroscopy.The principal doping impurity concentration did not exceed 10 14 cm −3 .

Identical
The experimental geometry was chosen to locate the samples in respect to the electron irradiation beam, such that the maximum possible setting of the irradiation parameters for the whole batch was achieved.The electron irradiation was performed in all, I, II, III, IV types with the following electron beam energies: 3.5, 14, 25 and 50 MeV, at the Yerevan Physics Institute, with a beam pulse of about 5 micro-seconds, while type V samples were irradiated by electrons with an energy of 3.5 MeV at the CANDLE synchrotron Radiation Institute (Armenia) with a beam pulse duration of about one pico-second.Infrared (IR) absorption spectra at 9 µm (for impurity oxygen) and at 12 µm (for impurity carbon) and Hall-effect measurements were made on the above mentioned samples, both before and after irradiation at room temperature.Details of the Hall-effect measurements [12] and of the IR measurements were as described in works [13] [14].Note that in the case of 5 micro-seconds beam irradiation a cooling of the samples was applied by liquid nitrogen vapors to the room temperature which do not necessary at pico-second irradiation because the thermal processes do not develop during this very short duration.
The experimental results are indicated in Figures 1-6.Characteristic IR absorption bands due to the localized vibrations from: interstitial oxygen atoms O I (1106 cm −1 ), substitution carbon atoms C s (606 cm −1 ), divacancies V 2 (606 cm −1 ), A-centres VO I (829 cm −1 ) and a complex VO I C I (866 cm −1 , where C I is an interstitial carbon atom), are evident in the spectra from irradiated samples of I, II, III types, while in the samples of type IV, only a band related to V 2 was detected.
Hall-effect measurements were carried out to determine the concentration of the main charge carriers, their mobility and formation rate of the radiation defects.
Figure 1 shows the dose dependences of the concentration of radiation induced divacancies for type I -IV samples, as irradiated by electrons with energies in the range 3.5 to 50 MeV.
It is apparent that an additional channel for I in n-Si crystals is introduced at high radiation doses, which begins to compete for the capture of I on the C s sites.
Due to the aforementioned higher concentration of VO I complexes this channel most likely involves the annihilation of I and V on the (VO I ) oxygen atoms.In all the samples investigated, the dose dependence of the produced V 2 concentration displayed an initial linear region (Figure 1) from the slope of which the introduction rate of V 2 at the accumulation stage could be estimated.It is clear from Figure 2 and Figure 3 that the divacancies (V 2 ) are most effectively introduced into those samples with an enhanced concentration of C s atoms, and this effect prevails across the entire electron energy range, as used in this work.It should be noted that the concentration of C I O I complexes is much lower than the decrease of C s concentration and hence, we may conclude that there are more effective sinks for carbon interstitials in the Si crystal, for example, the production of C s C I pairs.In all cases the formation of V 2 , the VO I (A-centre) and VO I C I complexes is accompanied by a decrease in the intensities of interstitial oxygen and substitution carbon absorption bands, with the change in O I atom concentration proceeding with a nearly constant rate across the electron irradiation energy (from 3.5 to 50 MeV), while the removal rate of C s atom increases with increasing energy (Figure 3).The intensities of the absorbtion bands of Oxygen atoms (O) and Carbon atoms (C) decrease during irradiation because some of these atoms join the mentioned complexes.
For compare results in more detail it is interesting to give numerical calculation concerning radiation defects introduction rate, as an example, for samples with specific resistivity about 100 Ohm•cm at room temperature: ∆N/∆D = 0.049 cm −1 curve 5 in Figure 1 (micro-second pulse duration 3.5 MeV energy); ∆N/∆D = 0.34 cm −1 in Figure 6 (pico-second pulse duration 3.5 MeV energy); ∆N/∆D = 0.24 cm −1 curve 1 in Figure 1 (micro-second pulse duration 50 MeV energy); These values confirm our previous results about clusters formation at pico-second pulse duration 3.5 MeV energy irradiation in spite of low energy, low dose and low intensity for this process.We called this phenomenon "Drop dose effect".

Discussion
In order to interpret these experimental data, we may consider in detail the generation of radiation defects in Si crystals under irradiation by high energy electrons, and how this varies as a function of electron energy and dose in the particular samples, with specified concentrations of impurities.Following an initial defined stage where simple point defects are introduced in the sample, the probability of V-V interactions begins to increase, and the formation of V 2 begins.These processes are further influenced by the number of charged impurity species which are created during irradiation as quasi-molecules with initial simple defects (therefore V 2 is formed at relatively low irradiation doses in pure Si) and low cross-sections for the recombination of species I with those V 2 .Furthermore, depending on the radiation energy the formation of V 2 as an initial defect occurs directly via a displaced atom, when this atom, in turn, displaces neighborhood atoms to small lattice constant distances.In this case, the distribution of V 2 over the sample is chaotic and the accumulation of these species increases linearly with dose.This is referred to as the formation of a "primary divacancy".An alternative mechanism is the radiation-induced formation of V 2 -clusters, when diffusion occurs from the micro-region, where the concentration of primary simple defects (the most mobile interstitials) is high, and secondary V 2 are formed in the voids of initially displaced regions, as a consequence of atom-atom cascade collisions.As in the previous case, the accumulation of V 2 is linear with dose, but the accumulation rate for this process is higher (Figure 1).This is called a "secondary divacancy".When the initially displaced atom has a low energy, E i < E r (E r -is the threshold of cluster formation, which in Si is 7.5 keV), the accumulation of secondary divacancies prevail (Figure 2), while primary vacancies dominate at high energies.A possible explanation of the "impurity effect" at the linear stage of V 2 accumulation is as follows: the C s atoms are able to efficiently capture the displaced atoms, which inhibits the generated (produced from the same lattice sites) correlated V 2 and I atoms from recombining, which promotes an increase in the primary V 2 yield.This is supported by the observation of a linear correlation between the increase in V 2 concentration and the decreasing concentration of C s atoms (Figure 4); capture of I atoms by C s atoms will result in their ejection and return to the C I interstitial position [1].
To estimate quantitatively the influence of the capture reaction of the I atoms in regard to their recombination with correlated V 2 , a computer modeling studies were performed of a pertinent process, according to the following scheme: a three-dimensional isotropic diffusion on a diamond-like lattice of particle A and subsequent reaction via one of two competing channels-a capture at a spatially fixed centre B or at randomly volume-distributed traps C. In order to the reaction to proceed, it was necessary for A to fall into the region of effective spontaneous capture by B or C, as defined by the radii R B or R C , respectively.Figure 5 presents the probability P AB dependence, of a capture at B at various values of initial distance r AB between A and B. It is clear from Figure 5 that P AB is insensitive to the C traps, at concentrations less than 10 19 cm −3 .
The reason for this is that the probability of the particle returning, by executing random walks in an unbounded three-dimensional medium, back to the initial point is less than unity [15] [16] [17], i.e.A will "lose its way" in the infinite volume and will not fall into the B capture band even in the absence of traps.Hence, given the three-dimensional isotropic diffusion of I atoms, the presence of the impurities capable of capturing them must not affect the linear accumulation of V 2 .
The situation is different when the diffusion of the displaced I atom is strongly anisotropic.Thus, in the simplest limiting case of one-dimensional motion of the I atom along a chain connecting V 2 and the impurity C s atom, a recombination of V 2 and the I atom will definitely occur if the latter is not captured by the C s atom, and it is the latter that determines the efficiency of the influence of impurities on the primary V 2 yield.Thus, a quantitative interpretation of experimental results is possible in a model involving quasi-one-dimensional diffusion of I atoms [12].In other words the computer modeling simulates the capture process of A particles by B and C traps.The C type traps are defined by three-dimensional isotropic diffusion of atoms displaced by irradiation, whereas B traps-quasi-one dimensional diffusion.The computer calculation has showed that C traps don't fit experimental results even depending on irradiation energy (Figures 3-5).
An interesting result is obtained from ultrafast electron irradiation of the samples on the pico-second timescale, as is clear from Figure 6, which shows the dose dependence of the concentration of radiation defects, as induced by electrons with an of energy 3.5 MeV.In this case, the dose dependence behavior can been seen to differ from that in Figure 1, which may be explained in terms of the very great difference in the timescales over which the electron beam has been applied.In the case of the ultrafast pulse, a cluster which consists of A-centre (VO) type radiation defects is induced in the silicon crystal, despite the lower irradiation energy for cluster formation [9].

Conclusion
Results are presented of an investigation of the interaction of displaced Si-self atoms and their vacancies with impurities in crystalline silicon, as a result of irradiation by electrons of different energies, taking into account the influence of the pico-second duration of the electron beam.A process involving the accumulation of radiation defects in the dose-dependent, linear region is studied in detail.The effect of impurities on the recombination of correlated divacancies, with atoms displaced from regular lattice points, is estimated by a computer modeling of an appropriate diffusion-controlled process.It was shown that there is a substantial difference in the introduction rate of radiation defects in Si samples with different contents of impurities, and also when ultrafast electron irradiation is used.The work has particular novelty in that it provides the first analysis of the energy and dose dependences of the introduction and accumulation of How to cite this paper: Yeritsyan, H.N., Sahakyan, A.A., Grigoryan, N.E., Harutyunyan, V.V., Grigoryan, B.A., Amatuni, G.A., Yeremyan, A.S. and Rhodes, C.J. (2018) Quasi-Chemical Reactions in Irradiated Silicon Crystals with Regard to Ultrafast Irradiation.Journal of Modern Physics, 9, 1271-1280.https://doi.org/10.4236/jmp.2018.96075 batches of plane-parallel samples, with the (111) crystallographic orientated to a larger face-plane, were fabricated from ingots of n-Si single crystals.Each batch consisted of five types of samples, which differed only in their concentrations of the dominant oxygen N O and carbon N C impurities: Type I (n-Si)-Specific resistivity 100 Ohm•cm, N O = 5.0 × 10 17 cm −3 , N C = 3.0 × 10 17 cm −3 ; Type II (n-Si)-Specific resistivity 50 Ohm•cm, N O = 1.0 × 10 18 cm −3 , N C = 3.0 × 10 17 cm −3 ; Type III (p-Si)-Specific resistivity 20 Ohm•cm, N O = 1.4 × 10 18 cm −3 , N C = 1.5 × 10 17 cm −3 ; Type IV (n-Si)-Specific resistivity 500 Ohm•cm, N O ≤ 10 16 cm −3 , N C < 10 16 cm −3 ; Type V (n-Si)-Specific resistivity 124 Ohm•cm, N O = 5.0 × 10 17 cm −3 , N C = 2 × 10 17 cm −3 ;

Figure 4 .
Figure 4. Change in concentration of carbon substitutional atoms N c at introduction of V 2 into the samples type I (a, d), II (b, e) and III (c, f) irradiated by 50 MeV (a, b, c) and 3.5 MeV (d, e, f) energy electrons.

Figure 5 .
Figure 5. Recombination probability P AB as a function of the C traps concentration at R B = 1a and R C = 0.2a (a is the lattice constant).

Figure 6 .
Figure 6.Radiation defects concentration N def as a function of irradiation dose D in the samples of type V irradiated by 3.5 MeV energy electrons with pulse duration of pico-second region.