# Two Kinds of Hydrogen Monomers Manifested in Plasma-Exposed Silicon

Authors: Vladimir V. Voronkov
Journal: Advances in Materials Science and Engineering
 (2018)

---

## Abstract

In boron-doped silicon annealed in a plasma ambient (at 150°C), the reported hydrogen concentration profile and the hole profile cannot be simultaneously fitted assuming only one kind of in-diffusing hydrogen ions H+ of a definite parameter D+K (where D+ is the diffusivity of H+ and K is the equilibrium dissociation constant of the HB defect, the passivated boron). A good fit is possible only assuming two independent kinds of H+—one of a larger value of D+K and the other—of a smaller value. A concept of two independent atomic subsystems H(1) and H(2), each involving both positive and neutral charge states, is also useful to account for hydrogen pairing into dimers.


---

## Body

## 1. Introduction

Hydrogen impurity plays an important role in silicon materials—by passivating dopants and other defects [1–6]. In p-type and near-intrinsic Si, under equilibrium conditions, the dominant monomeric state of hydrogen is thought [7] to be H+(BC)—a positive ion in a bond-centred position. The acceptors (normally, boron B−) are passivated by combining H+ and B− into neutral defects HB. The data on hydrogen diffusion from a plasma ambient and on the related boron passivation were collected years ago, but they are still waiting for a proper analysis, in view of a growing interest to hydrogen as a powerful and promising tool to control the defects in Si.The degree of passivation—a ratio of passivated boron concentration [HB] and nonpassivated boron concentration [B−]—is controlled by the concentration [H+] of free (nontrapped) hydrogen ions and the equilibrium dissociation constant K of the HB defect:(1)H+B−HB=KT.A depth profile of hydrogen concentrationC(z) includes a near-surface narrow region of a high concentration, and a bulk part where hydrogen is represented by free and trapped atomic species, H+ and HB. The concentration [H+] turns out [8] to be much smaller than [HB] or [B−]. The shape of the C(z) profile, for a specified hydrogen concentration at the sample surface, depends [8] on a single material parameter D+K where D+ is the diffusivity of H+. The value of D+K based on the reported C(z) profiles is about 6.5 × 104 cm−1s−1 at 150°C being slightly different for different samples.For one of the samples annealed at 150°C [9], also the hole profile p(z) is available along with C(z). Quite unexpectedly, fitting the p(z) profile (described in the next section) gives an essentially smaller value of D+K—about 104 cm−1s−1.This seeming discrepancy is resolved if there are two independently diffusing hydrogen ions that we denote H+(1) and H+(2)—one of a larger D+K and the other of a smaller D+K. In the present paper, a theory of boron passivation by two independent ions H+(1) and H+(2) is considered and applied to fit some representative experimental profiles.


## 2. Fitting Hole Profile with One Kind of H+

For a sample of a doping levelNB = 5 × 1018 cm−3 annealed at 150°C for 30 min, both the concentration profile C(z) by SIMS (using the deuterium isotope of hydrogen) and the spreading resistance profile R(z) were reported [9]. The resistance normalized by its bulk value Rb is equal to the normalized specific resistivity ρ(z)/ρb. With the known bulk resistivity ρb, the resistivity profile ρ(z) is defined and converted into p(z) using a standard ASTM procedure. This p(z) profile is shown in Figure 1 by circles. The concentration profile C(z) for the same sample is displayed in Figure 2(a).Figure 1
Hole depth profile in a sample ofNB = 5 × 1018 cm−3 annealed at 150°C for 30 min. (a) Fitting by in-diffusion of one kind of H+. Curve 1: with D+K = 8 × 104 cm−1s−1, X = 21.2; curve 2: with D+K = 1.4 × 104 cm−1s−1, X = 115. (b) Fitting by in-diffusion of two kinds of H+ with the parameters specified in the text.
(a)
(b)Figure 2
Concentration profiles for the sample of Figure1. (a) Solid curve is computed for the same parameters that were used in Figure 1(b), and the experimental points (circles) are rescaled to match the experimental penetration depth by Figure 1. (b) Individual components of the passivated boron.
(a)
(b)Sincep = [B−] − [H+] and C = [HB] + [H+], the two profiles p(z) and C(z) are complementary: their sum is a fixed total boron concentration NB. The profiles can be simulated [8] by solving an equation for diffusion and drift of H+ accompanied by fast trapping and detrapping of H+ by boron, which maintains the equilibrium relation (1) between the reacting species. The solution depends on the value of D+K and on the concentration of H+ at the sample surface normalized by K; this ambient-controlled boundary parameter, [H+]/K, is denoted X.In Figure1(a), the curve 1 was calculated with D+K = 8 × 104 cm−1s−1—the value deduced by fitting the concentration profile of Figure 2(a). The boundary parameter X was set to 21.2 to reproduce the penetration depth of hydrogen, 0.57 μm. Clearly, the computed curve is inconsistent with the experimental hole profile. A much better (although not quite satisfactory) fit is achieved with an essentially smaller value of D+K = 1.4 × 104 cm−1s−1 and a larger X = 115 (the curve 2 in Figure 1(a)). This smaller value of D+K is however inconsistent with the shape of the C(z) profile in Figure 2(a).


## 3. A Concept of Two Kinds of H+ Ions

To resolve this contradiction, we note that the information contained in theC(z) and p(z) profiles actually refers to different parts of the profiles.TheC(z) bulk profile is comprised of a plateau portion (marked P in Figure 2(a)) and a tail portion (marked T). In the plateau, a small difference between C and NB cannot be accurately determined, and the informative part of the profile is only the tail. Therefore, the deduced D+K parameter refers to the hydrogen ions within the tail.For thep(z) profile, the nonpassivated boron concentration [B−] is directly defined in the whole in-diffused region, allowing to find D+K for the H+ ions within the plateau.It can be thus concluded that there are actually two independent passivating hydrogen ions: H+(1) of a larger D+K dominating in the tail portion of the presently discussed hydrogen profile, and H+(2) of a smaller D+K dominating in the plateau part. The first one is assumed to be the abovementioned bond-centred ion H+(BC). The second one is in a different lattice location to be discussed later. Independent diffusion of these two species implies that there is a high energy barrier for transition from one state of H+ to the other which prevents equilibration between them.The concentration ratio of the two species H+(1) and H+(2) at the sample surface depends on the plasma ambient, and it can have a strongly nonequilibrium value. It is thus possible that the H+(2) species is more significant under plasma exposure than under the equilibrium conditions.The hydrogen in-diffusion and resulting boron passivation should be now reconsidered taking into account simultaneous presence of two kinds of H+ ions.


## 4. Boron Passivation by Two Independent In-Diffusing H+ Ions

The inferred existence of two independent H+ ions is most essential for a plasma-induced acceptor (boron) passivation. In this case, there are two different structural forms of passivated boron HB: the HB(1) originates from trapping H+(1), and HB(2) from trapping H+(2). A transport equation of Reference [8] is applicable to each of the two H+/HB subsystems:(2)∂HBI∂t=−∂JI∂z,wherethefluxJI=−D+Ip∂H+I/p∂z.The subsystem index I is either 1 or 2. The concentration [HB(I)] in the left-hand part represents the total concentration for the subsystem I, neglecting a small contribution of free ions [8]. The flux J(I) in Equation (2) includes a drift of positive hydrogen ions in the electric field related to a nonuniformity in p(z). The concentration [H+(1)] or [H+(2)] of free ions is expressed through the major concentrations [HB(1)] or [HB(2)] using Equation (1)—now applied to each species:(3)H+I=KIHBIB−.After that, each of the two transport equations (2) contains only one parameter, either D+(1)K(1) or D+(2)K(2). At 150°C, D+(1)K(1) is already defined to be 8 × 104 cm−1s−1 while D+(2)K(2) is on the order of 104 cm−1s−1, and only a small adjustment of this parameter is required; the best-fit value was found to be D+(2)K(2) = 9.6 × 103 cm−1s−1.The hole concentrationp = [B−] − [H+(1)] − [H+(2)] is close to that of nonpassivated boron [B−]:(4)p=B−=NB−HB1−HB2.This is the only equation through which one subsystem affects the other.The diffusion problem involves two ambient-dependent boundary parametersX1=H+1/K1 and X2=H+2/K2. The surface value of [B−] is expressed through these parameters from Equations (4) and (3): B−=NB/1+X1+X2.The best fit of the hole profile is shown by the solid curve in Figure1(b), computed with X(1) = 4.3 and X(2) = 134. This curve reproduces the measured p(z) profile much better than the curves 1 and 2 in Figure 1(a), which are based on only one kind of H+ ion.The computed concentration profileC(z)—for the parameter set specified above—is shown by the solid curve in Figure 2(a). The reported hydrogen penetration depth [9] was a bit larger for the C(z) profile in comparison with the p(z) profile. This small discrepancy may be due to some inaccuracy in a crater depth z in the SIMS technique, and it was adjusted by rescaling the experimental C(z) profile—by multiplying all the values of z by a scaling factor 0.86. This rescaled profile—shown in Figure 2(a) by the circles—is well consistent with the predicted profile, the solid curve.The near-surface high-concentration region (marked S in Figure2(a)) was modelled here as in-diffusion of slow hydrogen dimers H2C [10] produced at the surface. The concentration profile of dimers is described by an erfc function with an apparent diffusivity 3 × 10−15 cm2/s, which is close to that expected for H2C [10]. The dimeric contribution 2[H2C] is simply added to that resulting from in-diffusion and trapping of H+(1) and H+(2).To illustrate a specific effect of a simultaneous transport of two H+/HB subsystems, the two components of the passivated boron, [HB(1)] and [HB(2)], are shown separately in Figure 2(b). The HB(2) dominates in the major part of the passivated region, while HB(1) dominates within the tail portion. A qualitative explanation is that the species H+(2), due to a large X2=H+2/K2, occupy the major part of boron traps down to some depth, while the H+(1) species, facing a low concentration of remaining B− traps, easily penetrate to a larger depth and thus control the tail part of the profile.


## 5. Fitting Concentration Profile for Lower Doping Level

Another example of a concentration profile with a well-resolved tail (NB = 1017 cm−3, annealed at 150°C for 30 min [11]), is shown in Figure 3(a). A peculiarity of this profile is a relatively wide S-component (in comparison with Figure 2(a)) spreading down to a micron. This scale is consistent with a diffusion length of H2A dimers—moderately mobile species observed in samples saturated with hydrogen at high T and subsequently quenched [12, 13]. This kind of dimers is formed by a pairing reaction H0(1) + H+(1) at a known rate [10]; the neutral monomer H0(1) is a minor but extremely fast species located in a tetrahedral interstice. However, the computed concentration [H2A], due to this pairing reaction, turns out to be much smaller than NB [8] and cannot account for the S-component.Figure 3
Concentration profiles for a sample ofNB = 1017 cm−3 annealed at 150°C for 30 min [11]. (a) Total concentration due to boron passivation and formation of H2A dimers by H0(1) + HB(2). (b) Individual components of the passivated boron.
(a)
(b)Now, with two independent subsystems of monomeric hydrogen, more pairing reactions can proceed involving H0(1) and H+(2) or HB(2) and these ones can be more efficient. The reaction H0(1) + HB(2) seems preferable since [HB(2)] is much larger than [H+(2)]. The diffusion-limited pairing rate is 4πrD0(1) [H0(1)] [HB(2)] where r is the capture (pairing) radius, conventionally assumed to be 0.5 nm. The ratio [H+(1)]/[H0(1)] equals p/pd(1) where pd(1) is the characteristic hole concentration for the Fermi level coincident with the donor level for the H0(1) to H+(1) transition. Using Equation (3), the pairing rate is expressed through the two major concentrations [HB(1)] and [HB(2)]:(5)G=βHB1HB2p2,where the pairing kinetic coefficient β, for the diffusion-limited reaction, is(6)β=4πrD01pd1K1.The actual value ofβ can be smaller, and it should be treated as a fitting parameter with the upper limit given by Equation (6). All the parameters in Equation (6) were specified [8], and the upper limit for β is 6 × 1012 cm−3s−1 at 150°C.The dimeric concentration [H2A] is incremented at a rate G, while each concentration [HB(1)] and [HB(2)] is diminished at the same rate G. Each transport equation (2) is thus modified by adding a term G in the right-hand part.The equation set now includes also a diffusion equation for produced H2A dimers:(7)∂H2A∂t=D2A∂2H2A∂z2+G.Backward dissociation of H2A into H01+HB2 proceeds by H2A+B−+h+ at a rate proportional to p2. It is neglected at a relatively low temperature of 150°C and a moderate hole concentration.To reproduce theC(z) profile in Figure 3(a), we use for D+(2)K(2) the previously deduced value, 9.6 × 103 cm−1s−1. The parameter D+(1)K(1) was a bit reduced, down to 5 × 104 cm−1s−1, for a better reproduction of the tail. The diffusivity D2A of H2A dimer is known [10, 14, 15] to be 2 × 10−13 cm2/s at 150°C (for the deuterium isotope). The remaining parameters are the two boundary ratios X(1) and X(2), and the pairing coefficient β; their fitted values are X(1) = 10.3, X(2) = 21.7 and β = 4 × 1012 cm−3s−1. The latter value is only slightly smaller than the upper limit of 6 × 1012 cm−3s−1. The corresponding best-fit profile is shown by the solid curve in Figure 3(a). The computed depth dependence of the pairing rate shows that the H2A dimers are produced within a near-surface zone of about 0.3 μm; they out-diffuse to the surface and in-diffuse into the sample bulk.For illustration, also the individual profiles of HB(1) and HB(2) are displayed in Figure3(b). Here, the hydrogen species of the second subsystem are spent on pairing near the surface, and they do not penetrate deep into the bulk. The bulk part of the profile is occupied mostly by HB(1).The ratio of [H2A] and [HB(2)] is proportional to [HB(1)]/p2, by Equation (5). It is relatively large for the case of Figure 3(a) but becomes reduced at a higher doping level NB (a higher p). Besides, the dissociation of H2A may become significant at a higher p. For this reason, the profiles in Figures 1 and 2 (computed with a neglected pairing) are hardly changed when recalculated including the pairing.Another example, for the same doping levelNB = 1017 cm−3 and the same annealing temperature 150°C, taken from Reference [16], is shown in Figure 4(a). This profile has quite a similar shape but the width of the S-component corresponds to a smaller diffusivity D2A, and the tail is relatively steep corresponding to a smaller D+(1)K(1). This is understood if, in spite of the same nominal temperature for the profiles of Figures 3(a) and 4(a), the actual temperature is somewhat lower (roughly, by 15°C) for the latter case. To fit the tail of the profile of Figure 4(a), the parameter D+(1)K(1) was reduced down to 104 cm−1s−1, and D+(2)K(2) was reduced by the same factor, down to 2 × 103 cm−1s−1. The other parameters (β, D2A, X(1) and X(2)) were set to achieve the best fit: β = 2 × 1012 cm−3s−1, D2A = 8 × 10−14 cm2/s, X(1) = 14.5, X(2) = 18.7. A narrow near-surface part of the S-component in Figure 4(a) was modelled, like in Figure 2(a) by in-diffusion of slow dimers H2C.Figure 4
Concentration profiles for a sample ofNB = 1017 cm−3 annealed at 150°C for 2 h [16]. (a) Total concentration due to boron passivation and formation of H2A dimers by H0(1) + HB(2). (b) Individual components of the passivated boron.
(a)
(b)The individual contributions of HB(1) and HB(2) into boron passivation are shown in Figure4(b); these profiles are similar to those of Figure 3(b): the bulk is dominated by HB(1) species.


## 6. In-Diffusion Profiles in Lightly Doped n-Si

In lightly doped n-Si, the observed concentration profilesC(z) include a near-surface component (labelled S) and a bulk component (labelled B). Two representative examples [3, 17] are shown in Figure 5. The profile is insensitive to variations in the electron concentration suggesting that it is formed by neutral hydrogen species while H+ is not now essential. Each component is well described as resulting from pairing of some in-diffusing species of a diffusivity D at a rate(8)G=4πrDc2,where cz is a steady-state depth profile of the in-diffusing species and r is the capture (pairing) radius. With a variable S=Dc used instead of c, the depth profile obeys a simple equation:(9)d2Sdz2=2G=8πrDS2.Figure 5
Concentration profiles in lightly doped n-Si. (a) Annealing at 150°C for 1 h [3]. (b) Annealing at 125°C for 1 h [17]. Each component of a profile, S and B, is described by Equation (11).
(a)
(b)The solution of Equation (9) depends on one material parameter D/r and one boundary parameter S0—a surface value of Sz:(10)Sz=S01+z/L2.The characteristic penetration depthL equals D/r/4πS0/31/2.The total hydrogen concentrationCz is mostly due to produced pairs of in-diffusing defects (while a contribution c of the source species themselves is negligible [8]). The diffusion of produced pairs is neglected on the time and depth scales of the profiles; then,(11)Cz=2Gzt=2G0t1+z/L4.With the deduced shape parametersL and G0, the ratio D/r is calculated as 4π/9G0L4.For the S-component, the deduced ratioD/r [8] is close to the value expected for a hydrogen species responsible for a slow stage of tritium loss at 400 to 500°C [18]—from samples saturated with the tritium isotope at high T and quenched. The diffusivity of this mysterious species is much higher than D2A but much smaller than the known atomic hydrogen diffusivity [19], and this species was tentatively identified [10] as a “fast dimer” H2B. Accordingly, the S-component was attributed to H2B pairing into tetramers. An alternative treatment is possible, but it will not be discussed in the present paper. Regarding the B-component, there were two difficulties [8]:(1)
The value ofD/r for this component was orders of magnitude smaller than the value of D/r for the neutral tetrahedral monomer H0(1).(2)
The deduced values ofD/r were strongly different for two close temperatures: 150°C (Figure 4(a), D/r = 0.8 cm/s) and 125°C (Figure 4(b), D/r = 0.015 cm/s). A huge difference—by a factor of 53—is too large for any reasonable temperature dependence of D/r.Now, a concept of two different atomic subsystems offers a simple solution to these puzzles if we assume that H+(2) is a part of a subsystem H(2) that includes also a neutral species H0(2), dominating in intrinsic and n-type Si (which implies that the donor level of H(2) is below the midgap). The B-component can be then ascribed to dimers produced by interpairing H0(1) + H0(2) or self-pairing H0(2) + H0(2). There are two extreme cases:Case 1. The equilibration time between H0(1) and H0(2) is shorter than the anneal duration. Then, the two species coexist in the equilibrium ratio and diffuse together, as one species of an apparent diffusivity Da averaged over the two species:(12)Da=D01+D02R211+R21,where R21 is the equilibrium concentration ratio of H02 and H01. In the D/r ratio deduced from the shape of the B-component, the apparent diffusivity D is equal to Da. We assume that this situation is valid at a higher T = 150°C—for the profile of Figure 5(a). With a diffusion-limited capture radius, r = 0.5 nm, the diffusivity Da = 4 × 10−8 cm2/s. The diffusivity of atomic hydrogen [19] extrapolated down to 150°C is 1.4 × 10−8 cm2/s (for the deuterium isotope). This diffusivity is averaged over all the atomic species present in intrinsic Si, including the major one, H+(1), and therefore it is smaller than Da (which is averaged only over the neutral states). Yet, the difference between the two diffusivities is not very large showing that H0(2) in intrinsic Si gives a small but significant contribution to the total concentration.Case 2. The equilibration time is longer than the anneal duration. Then, the two species diffuse independent of each other. If, in addition, the concentration of H0(1) is so small that interpairing can be neglected in comparison with self-pairing of H0(2), then the apparent diffusivity D will be identical to D0(2). This situation may hold at a lower temperature of 125°C. Assuming again a diffusion-limited capture radius, r = 0.5 nm, we specify the diffusivity of the neutral H0(2) species at 125°C: D0(2) = 7.5 × 10−10 cm2/s. If a representative diffusivity prefactor, 0.01 cm2/s is further adopted, the migration energy of H0(2) is estimated to be 0.56 eV. The D0(2) diffusivity extrapolated from 125°C to 150°C is then 2 × 10−9 cm2/s—much smaller than the averaged diffusivity Da deduced above to be 4 × 10−8 cm2/s. By Equation (12), it is concluded that the hydrogen transport at 150°C is mostly due to H0(1), while the total concentration of the neutral species is mostly due to H0(2). The expression (12) is then simplified to Da = D0(1)/R21.A huge difference between the apparent diffusivitiesD at 125 and 150°C is now attributed to a different meaning of the apparent diffusivity D for these two cases: a relatively low diffusivity D0(2) at a lower T, and a much higher average diffusivity Da (enhanced by a contribution of fast H0(1) species) at a higher T.


## 7. Summary

Penetration of H+ hydrogen ions from a plasma ambient into a boron-doped sample is limited by their trapping by boron into HB neutral defects. It is controlled by the surface concentration of H+ and by a single material parameter D+K—a product of the H+ diffusivity D+ and the equilibrium dissociation constant K of the HB defects. Under an assumption of only one kind of H+, of a definite value of D+K, it is impossible to reproduce simultaneously the hydrogen concentration profile C(z) and the hole profile p(z) measured for the same sample at 150°C: the former can be well fitted only with a relatively large D+K while the latter can be roughly fitted only with a much smaller D+K.This difficulty is overcome by assuming two different kinds of H+ ions, H+(1) and H+(2), one of a larger D+K and the other of a smaller D+K. They in-diffuse—without an exchange between them—and get trapped by B− into two different kinds of passivated boron, HB(1) and HB(2). A transport of each species is affected by the presence of the other species through the hole concentration that depends on the total amount of passivated boron, [HB(1)] + [HB(2)].The deducedD+K at 150°C is about 6 × 104 cm−1s−1 for H+(1) and close to 104 cm−1s−1 for H+(2). Relative contributions of the two species depend on their concentrations at the sample surface controlled by a composition of the ambient (hydrogen plasma). In one example, H+(2) was found to be the major passivating species down to some depth, leaving a low concentration of nonoccupied boron traps. In that case, H+(1) easily diffuse through the passivated region and control the hydrogen concentration within a tail of the profile.In another example (at a lower doping level), the H+(2) ions are spent on formation of hydrogen dimers and do not penetrate deep; the sample bulk is then dominated by the species of the first kind.It is assumed that H+(1) is an ion in the bond-centred position that is dominant under the equilibrium conditions in intrinsic Si. This one is equilibrated with a minor but highly mobile neutral species H0(1) in the tetrahedral location. Also, the H+(2) species seems to be present along with a neutral state H0(2), of a much higher concentration than H0(1). In the early stage of discussing the hydrogen properties in Si, it was argued that, due to its small size, a hydrogen atom should occupy a tetrahedral interstice remaining neutral. Now, we revive this notion assuming that the H0(1) is actually an excited state of the neutral tetrahedral hydrogen while the ground state corresponds to H0(2). There should be a high energy barrier for transition between these two states, which retards an exchange between them. The reported hydrogen profiles in lightly doped n-Si for 125 and 150°C suggest that the exchange is negligible at 125°C but relatively fast at 150°C.The hydrogen profiles normally include a near-surface part of a very high concentration. This “S-component” of a profile cannot be ascribed to the unique reason; it is specific for particular samples:(i)
In near-intrinsic Si, the S-component is clearly caused by self-pairing of some in-diffusing neutral species different from H0(1) and H0(2)—presumably the fast hydrogen dimer H2B. The true nature of this species is still to be established.(ii)
In samples moderately doped with boron, the S-component (relatively wide) can be attributed to formation of moderately mobile dimers H2A by interpairing reaction H0(1) + HB(2).(iii)
In heavily doped samples, the S-component is seemingly due to slow dimers H2C that are formed at the surface (or close to the surface) and in-diffuse into the bulk.



---

*Source: 2385438-2018-10-11.xml*