Introduction next up previous
Next: Method Up: Dynamic Properties of Previous: Dynamic Properties of

Introduction

Carbon is a common and important impurity in silicon, typically occurring in concentrations of around 10 tex2html_wrap_inline348 -10 tex2html_wrap_inline350 cm tex2html_wrap_inline352 , and predominantly occupying substitutional sites, C tex2html_wrap_inline340 . Upon irradiation, mobile silicon interstitials (Si tex2html_wrap_inline334 ) can be captured by C tex2html_wrap_inline340 , forming the interstitial carbon defect C tex2html_wrap_inline334 [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. The defect is characterized by two local vibrational bands at 920 and 931 cm tex2html_wrap_inline338 , whose absorption intensities are in a ratio of approximate 2:1. This suggested to the earliest workers a trigonal defect [2]. However, later work with lightly irradiated p-type Si was able to correlate an electron paramagnetic resonance, EPR, spectrum with the infra-red absorption. The EPR showed that C tex2html_wrap_inline364 possessed tex2html_wrap_inline366 symmetry [4]. The hyperfine splitting on tex2html_wrap_inline368 C showed that the highest occupied state in C tex2html_wrap_inline364 consisted of a localized p-orbital centered on the C atom. Deep level transient spectroscopic (DLTS) measurements gave a donor level at tex2html_wrap_inline374 +0.28 eV. These results suggested a model of the defect shown in Fig. 1.

More recent EPR and DLTS investigations on n-type material [5] found an acceptor level at tex2html_wrap_inline376 -0.10 eV and that C tex2html_wrap_inline378 also possessed tex2html_wrap_inline366 symmetry. The similarity between the re-orientation barriers of the defects in both charge states ( tex2html_wrap_inline382 0.8 eV) indicated that they possess the same structure. However, no hyperfine splitting of the EPR line due to C tex2html_wrap_inline378 with tex2html_wrap_inline368 C or tex2html_wrap_inline388 Si was observed, possibly because the partially occupied acceptor state is extended over a number of Si atoms and has only a small overlap with carbon.

The defect is associated with a 856 meV photoluminescence (PL) line and uniaxial stress measurements revealed that the symmetry of the defect was also tex2html_wrap_inline366 (monoclinic I) [6]. It is just possible that the neutral defect possesses trigonal symmetry and reorients to tex2html_wrap_inline366 in a charged state when a carrier is thermally trapped before radiative recombination occurs. However, unambiguous confirmation that the symmetry of the neutral defect is indeed tex2html_wrap_inline366 and not trigonal came from the response of the two LVMs of the neutral defect to uniaxial stress [7]. These experiments showed that the symmetries of the 920 and 931 cm tex2html_wrap_inline338 LVMs were tex2html_wrap_inline398 and tex2html_wrap_inline400 respectively. Thus it appears that the defect adopts the structure shown in Fig. 1 in all three charge states and consists of a C and Si atom sharing a single lattice site with a tex2html_wrap_inline402 100 tex2html_wrap_inline404 orientation [4, 5, 7]. However, this model gives no natural explanation of the near degeneracy of the local vibrational modes of the defect, or their 2:1 ratio in absorption intensities.

Several theoretical calculations have been performed on the C tex2html_wrap_inline334 defect, all finding the tex2html_wrap_inline402 100 tex2html_wrap_inline404 structure to be the ground state configuration [8, 9, 10, 11, 12]. The calculated vibrational modes of the various defect configurations [9] provides further evidence for the tex2html_wrap_inline366 configuration in Fig. 1. For this model they were 1005 cm tex2html_wrap_inline338 and 853 cm tex2html_wrap_inline338 for the tex2html_wrap_inline398 and the tex2html_wrap_inline400 modes respectively, in fair agreement with experiment. None of the three possible tex2html_wrap_inline422 symmetry configurations possessed two high frequency modes.

The energy barrier to reorientation of the defect among the equivalent tex2html_wrap_inline402 100 tex2html_wrap_inline404 orientations has been calculated by Tersoff [10]. He found this to be 0.7 eV for a saddle point where the C atom lies at a bond centered (BC) site along tex2html_wrap_inline402 111 tex2html_wrap_inline404 . This is in good agreement with experimental reorientation and migration barriers around 0.8 eV [4, 5, 13]. More recent calculations on the migration/reorientation barriers of the defect have calculated these to be 0.51 eV via a saddle point with tex2html_wrap_inline432 symmetry [11] and 0.77 eV via a H-site [12] The migration/reorientation pathway considered by Capaz et al [11] involved a low energy, 0.51 eV, four-fold coordinated carbon atom at the saddle point. This path involves C moving through the lattice and at the same time reorienting and hence explains why the barriers to migration and reorientation are the same. The path via a BC defect was also considered but they found it involve a 2.5 eV barrier which is considerably greater than that found by Tersoff. The structural data from this recent ab initio supercell calculation [11] for the neutral defect is given in Table 1.

The di-carbon defect C tex2html_wrap_inline334 -C tex2html_wrap_inline340 is formed when Si:C, containing a low oxygen-content, is e-irradiated at room temperature creating mobile C tex2html_wrap_inline334 defects which are subsequently trapped by C tex2html_wrap_inline340 [1]. The defect is bi-stable taking the A-form in the charged states and the B-form in the neutral one. EPR and DLTS studies [14] showed that (C tex2html_wrap_inline334 -C tex2html_wrap_inline340 ) tex2html_wrap_inline448 defects possessed tex2html_wrap_inline450 symmetry, and they gave rise to acceptor and donor levels at tex2html_wrap_inline376 -0.17 and tex2html_wrap_inline374 +0.09 eV. The EPR data also found that the carbon atoms were inequivalent, and the g-tensors were perturbed by stress in a very similar way to those of the charged C tex2html_wrap_inline334 defect. The model proposed from the experimental data is very similar to that of C tex2html_wrap_inline334 , the symmetry being lowered by the presence of a second carbon atom at a substitutional site shown in Fig. 2. This configuration of the C tex2html_wrap_inline334 -C tex2html_wrap_inline340 defect is known as the A-form, and is the stable configuration for both singly ionized states.

For high e-irradiation fluences, the Fermi level is around mid gap and the neutral B-form of the defect is formed [14, 15]. This state is diamagnetic, but an excited triplet state B tex2html_wrap_inline468 has been observed by optically detected magnetic resonance (ODMR) [16]. This revealed that the symmetry of the B-form is tex2html_wrap_inline422 for T > 30 K, and tex2html_wrap_inline450 otherwise. Hyperfine interactions on tex2html_wrap_inline368 C detected by ODMR indicated that the carbon atoms are nearly equivalent. DLTS studies on the metastable B tex2html_wrap_inline448 form, which has trapped a carrier, show that has shallow donor and acceptor levels at tex2html_wrap_inline374 +0.07 and tex2html_wrap_inline376 -0.11 eV respectively. It has also been possible to observe EPR from the metastable B tex2html_wrap_inline484 form, but not from B tex2html_wrap_inline486 . B tex2html_wrap_inline484 is also found to have tex2html_wrap_inline422 symmetry for T > 15 K and again tex2html_wrap_inline450 below this temperature. The ODMR results suggest that the two C atoms lie close to neighboring substitutional sites, with a silicon interstitial, Si tex2html_wrap_inline334 , close to a BC site lying between them. This configuration is illustrated in Fig. 3.

There are three peculiar features with this model. Firstly, if Si tex2html_wrap_inline334 lay at the BC site, then the resulting tex2html_wrap_inline500 symmetry would give a partially occupied e-level. This would probably be split through a Jahn-Teller distortion resulting in Si tex2html_wrap_inline334 moving away from the BC site. However, the distortion cannot be too great if the reorientation energy around the C-C axis is to remain small. Nevertheless, the observed donor and acceptor levels are split by 0.97 eV, whereas a small distortion would lead to a difference in donor and acceptor levels of at most 0.5 eV being an estimate for the Hubbard U-correlation energy. On the other hand, if the distortion is not small, then the reorientation energy might be substantial and the symmetry of the defect would not be tex2html_wrap_inline422 .

A second difficulty comes from the vibrational modes of the defect. The B-form gives rise to a zero phonon luminescence line at 969 meV. The fine structure associated with this line reveals two local vibrational modes at 543.0 and 579.5 cm tex2html_wrap_inline338 ( tex2html_wrap_inline512 C) which shift by 10.1 and 14.9 cm tex2html_wrap_inline338 respectively for tex2html_wrap_inline368 C samples. However, in a sample containing a mixture of tex2html_wrap_inline512 C and tex2html_wrap_inline368 C, additional lines are observed which are shifted only by 0.1 and 0.6 cm tex2html_wrap_inline338 . Indeed, the original PL investigations failed to find distinct modes in the mixed case, and incorrectly suggested that the defect only contained one C atom. These tiny shifts confirmed that two carbon atoms were present in the defect, but implied that they are almost dynamically decoupled [17, 18]. It has been suggested that the this decoupling arises because the two C-Si tex2html_wrap_inline334 bonds are almost orthogonal to each other [14]. Even if this was the case, it would require that the C-Si tex2html_wrap_inline334 stretch modes lie around 550 cm tex2html_wrap_inline338 . This is unlikely since the reduced coordination of Si tex2html_wrap_inline334 will result in a shorter C-Si tex2html_wrap_inline334 bond, and hence the C-Si tex2html_wrap_inline334 stretch should be larger than the highest mode of C tex2html_wrap_inline340 which occurs at 607 cm tex2html_wrap_inline338 [19].

A third peculiarity is that the positive and negative charge states of the B-form are metastable, with energy differences with the A form of only a few hundredths of an eV in all three charge states. Similarly, the barriers between the forms are only 0.1-0.2 eV. As the bonding in the two forms is quite different and bond energies are of the order of 1 eV, it is surprising that these energy difference are so small.

Theoretical modelling of the di-carbon defect is still in its early stages, and to date three groups have examined the centre [9, 12, 24]. The first two focused on the relative energies of the A and B forms, and in our previous work, we analysed local modes of the B form only. The findings of the previous calculations relating to the defect energies are summarised in Table 4.

We discuss the standard models, and analyze their structure and vibrational modes in this paper. The method is discussed in section II, and applied to the C tex2html_wrap_inline334 defect in III and the di-carbon center in IV. We give our conclusions in V.


next up previous
Next: Method Up: Dynamic Properties of Previous: Dynamic Properties of

Antonio Resende
Wed Jan 15 12:41:08 GMT 1997