Electronic properties of dislocations

Dislocations exhibit a number of exceptional electronic properties resulting in a significant increase in the drain current of MOSFETs if defined numbers of these defects are placed in the channel. Measurements on individual dislocations in Si refer to a supermetallic conductivity. A model of the electronic structure of dislocations is proposed based on experimental measurements and tight-binding simulations. It is shown that the high strain level on the dislocation core—exceeding 10 % or more—causes locally dramatic changes in the band structure and results in the formation of a quantum well along the dislocation line. This explains experimental findings (two-dimensional electron gas, single-electron transitions). The energy quantization within the quantum well is most important for supermetallic conductivity.


Introduction
Dislocations are one-dimensional crystal defects influencing many of the physical and mechanical properties of crystalline solids [1,2]. The analysis of their electronic properties is of particular interest for semiconductor devices because dislocations detrimentally effect device parameters. This was proved, for instance, by electron-beam-induced current (EBIC) method [2,3]. Detailed investigations by various experimental techniques [Hall effect, electron paramagnetic resonance (EPR), deep-level transient spectroscopy (DLTS), photoluminescence] revealed numerous deep and shallow defect levels in the gap of Si related to dislocations [4][5][6][7][8]. Because most of the investigations in the past were carried out on deformed silicon, it is not clear which of the defect levels are related to the dislocation core, or to interactions with point defects in the strain field surrounding a dislocation [6]. Measurements on two-dimensional dislocation networks produced by wafer direct bonding avoid these interactions. Recently published DLTS data refer also to different deep and shallow levels in p-and n-type Si [9,10]. Their existence depends on the dislocation density. Shallow levels having activation energies of about 0.1 eV are discussed as possible reasons of radiative recombination (D1 luminescence) [11], while an interpretation of deep levels, related to the dislocation core, is still missing. Depending on defect formation, the energies of deep levels spread over a large energy range in the gap [6].
Dislocations also cause a significant increase in the drain current (I D ) of metal-oxide-semiconductor field-effect transistors (MOSFETs) if well-defined numbers of parallel dislocations are placed in the channel [12]. The largest increase in I D is obtained for a single dislocation [13]. In p-type silicon dislocations form channels for electrons with supermetallic conductivity which is about eight orders of magnitude higher than of the surrounding Si matrix [13]. This exceptional property and their dimension (length of a few micrometers and a cross-sectional area of the defect core of about 1 nm 2 ) characterize dislocations as native nanowires embedded in a perfect crystalline matrix.
In a previous paper, we expect that a reason of the supermetallic behavior is the high strain e in the dislocation core, which exceeds e % 10 % or more resulting in local modifications of the band structure and carrier mobility [13]. The present paper presents an experimental and theoretical description of the electronic properties of dislocations based on the strain in the dislocation core. This allows us to explain the supermetallic behavior as well as deep levels proved by DLTS. Selected results of experimental measurements of electrical properties of dislocations and of local strain fields present Chapter 3. Computer simulations of the band structure using tight-binding and kÁp models are discussed in Chapter 4. A model of the electronic structure of a dislocation is proposed in Chapter 5.

Experimental
The technique to realize defined two-dimensional arrays of dislocations in thin layers (semiconductor wafer direct bonding) was described elsewhere [13,14]. Silicon-on-insulator (SOI) wafers (diameter 150 nm, buried oxide (BOX) thickness 60 nm, device layer thickness 30 nm) were applied for the wafer bonding process under hydrophobic conditions. One of the handle wafers was removed after bonding resulting in a new SOI wafer having a device layer thickness of 60 nm and a two-dimensional dislocation network in the middle. SOI MOSFETs and pndiodes were prepared on such substrates (for details see [14]). Electrical measurements of device parameters were carried out in the temperature range from 300 to 3 K.
The defect structure was analyzed by applying a probe C S corrected scanning transmission electron microscope (FEI Titan 80-300). Local strain fields are quantified by peak-pairs analysis (PPA) of high-resolution electron microscope images.
All simulations of the band structure were performed by application of the nextnano 3 simulation package [15]. The bulk band structure was calculated using both a 6-and 8-band kÁp model, as well as an empirical sp 3 d 5 s* tightbinding model [16] taking into account strain and deformation potentials.

Experimental evaluation 3.1 Electrical measurements
Well-defined sets of parallel dislocations in the channel of nMOSFETs result in a dramatic increase in the drain current (I DS ). Figure 1 shows the transfer characteristics of devices without and with dislocations. The MOSFET without dislocations (reference) was fabricated on an identic SOI wafer without a dislocation network in the 60-nm-thick device layer and under the same process conditions. At V GS = 0.8 V and drain voltage (V DS ) of 0.4 V, the drain current is about 40 lA for MOSFETs containing dislocations in the channel while in the same time I DS is only 0.6 lA for a reference device. This corresponds to a current increase by a factor of 67. The increase in I DS depends on the type and the number of dislocations in the channel. The largest increase in I DS is observed if mixed dislocations are present.
The increase in the drain current for nMOSFETs results from an increased electron transport along dislocations. On the other hand, an increase in I DS exists also for pMOS-FETs where holes are involved in the carrier transport along dislocations. This means that depending on the properties of the surrounding material (p-or n-type), dislocations transport either electrons or holes. The dislocation core is occupied by electrons in p-type material and by holes in n-type material. This is in contrast to the basic model of electronic states on dislocations [4] assuming holes on the core of a dislocation in p-type material surrounded by the so-called Read cylinder of ionized (negatively charged) acceptor atoms.
Indications of a reinterpretation of the dislocation model provide also low-temperature measurements. Magnetoresistance measurements revealed Shubnikov-de Haas (SdH) oscillations referring to the presence of a two-dimensional electron gas (2DEG) on dislocations networks in p-type silicon (Fig. 2). The asymmetry of the SdH oscillations is frequently observed, but different reasons are discussed (for instance [17,18]). Because dislocations represent nanowires, the asymmetry may be attributed to the fact that field-dependent oscillations are extremely sensitive to the alignment between the field and the wire axis [19]. As a consequence of the 2DEG, single-electron transitions (Coulomb blockades) are observed [20]. Coulomb blockades are also found in dislocated n-type Si referring to the confinement of holes.

Electron microscopy
If {100}-oriented surfaces of Si are bonded together, a network of screw dislocations is formed in the bonded interface by the twist misorientation. The network is formed by two sets of screw dislocations with Burgers vectors b = \110[ oriented parallel to \110[ directions. Depending on the annealing conditions subsequently followed after the initial bonding, screw dislocations mostly dissociate into Shockley partials according to reactions of the type Both partial dislocations are oriented parallel to \110[ directions in the {111} glide plane. Their Burgers vectors are therefore inclined at an angle of 30°to the dislocation line direction. In addition, the dissociation reaction results in a stacking fault between both partial dislocations. Figure 3a shows a high-angle annular dark field (HAADF) microscopy image of a dissociated screw dislocation. The positions of the partial dislocations are indicated in the figure by circles. The strain field at the dislocation core and around the defect is quantified by peak-pairs analysis (PPA) [21]. PPA extracts the strain e by measurement of the lattice distortion. Analyses of e clearly proved that most of the strain is located on the dislocation core (Fig. 3b). Inside the defect, tensile strain of about 10 % exists, which is in accordance with theoretical predictions [22,23] and measurements on other defects [24]. The measured strain in the dislocation core is much higher than in the surrounding strain field. Furthermore, misfit of the tilt component of both bonded wafers results in an additional network of socalled 60°dislocations representing edge dislocations in the diamond cubic structure [25]. 60°dislocations dissociate also, but in 30°and 90°partials. The latter have Burgers vectors perpendicular to the dislocation line introducing normal strain, which is compressive above and tensile below the slip plane.

Band structure
The effect of different types of strain (uniaxial, biaxial, tensile, and compressive) on the band structure of Si and device properties has been analyzed using different simulation techniques [26][27][28][29][30][31]. The principal structure for our investigations is an nMOSFET schematically shown in Fig. 4. A dislocation is placed in the middle of the channel parallel to \110[. This introduces uniaxial strain parallel to the channel direction. Band structure calculations revealed shifts of the conduction and valence bands at the Cand L-point as well as band warping with increasing tensile strain. The conduction-band minimum at D is nearly constant with increasing tensile strain. If e = ?10 %, the lowest conduction band at the C point (DE C c ) is equal to the position of the lowest conduction band at D of the unstrained Si. As a result, the band gap energies of the indirect and direct transitions are equal for such high tensile strain and Si changes from the indirect into a direct semiconductor. This conclusion is in agreement with previously reported data [30,31]. Furthermore, the shift of the conduction-band minimum at the L-point results in an additional indirect transition in the presence of high tensile strain. Using the model of a dislocation surrounded by unstrained silicon (channel, Fig. 4), a schema of the band structure as in Fig. 5 follows. The shifts of the conduction band of the strained dislocation core result in one or more one-dimensional defect bands in the band gap of the surrounding Si channel. The depth of the defect bands increases as the tensile strain increases. Compressive strain, on the other hand, results in a dominant shift of the lowest conduction band at D (DE D C ). This induces a onedimensional defect band about 300 meV below the conduction band at moderate compressive strain (e B -0.04). Further increase in the dislocation-induced strain increases DE D C and causes the formation of a one-dimensional defect band close to the valence band if e is below -0.08.
Additional data about the variation of the band gap, band splitting, and effective masses with increasing strain are extracted from tight-binding simulations in combination with kÁp models. Figure 6a  varies with applied strain. At e = ? 0.1, the band gap of this transition is comparable to that of the indirect transition between C and D (Fig. 6a). Similar data have also been published for nanowires, where, however, the absolute data depend on the wire diameter [30]. It is obvious from Fig. 6a that the direct band gap decreases monotonically in the close vicinity of the band extremum [31]. The calculated hole mass m h = 0.285 is almost constant for -0.1 B e B ? 0.1 and corresponds to previously published data [31]. The effective in-plane electron mass parallel to the strain direction (m eff ) decreases continuously with increasing tensile strain (Fig. 6b). The decrease in m eff even at higher tensile strain is much less than predicted from measurements at lower tensile strain (e B ?0.01) [28,32].

Model of the electronic structure of dislocations
Tensile strain induced by a dislocation shifts locally both, the valence and conduction bands causing an additional indirect transition between C and L, as well as the direct transition at the C-point (Fig. 7). A cross-section through a dislocation yields a quantum well shown in Fig. 8a. Data from band structure calculations indicate the lowering of the minimum of the lowest conduction band at the L-point (E L CB ) below that at D of the unstrained Si (E D CB ). At the same time, the valence-band maximum of light and heavy holes is also shifted downwards (Fig. 7). For p-type silicon, discussed here, changes in the conduction band are most important.
Assuming the dislocation as a one-dimensional wire of radius R (Fig. 4), the Schrödinger equation is written in cylindrical coordinates as [33,34] Solving Eq. (4) results in quantized energy levels [35] with J 0 (r) as the first root of Bessel function. With m e = 0.19m 0 and R = 2.5 nm, one obtains the two lowest energy levels E 1 = 38.4 meV and E 2 = 203 meV. Both are inside the quantum wire formed by the dislocation (Fig. 8b). The density of states (DOS) of the quantum wire was also calculated. Based on the general case [33], the onedimensional density of states is derived as or where L = 2R and E i denotes the energy levels given by Eq. (5). Using R = 2.5 nm and E 1 -E 5 , the density of states for the quantum wire is plotted in Fig. 8c. The interpretation of dislocations as quantum wires explains most of the experimental results. The presence of SdH oscillations refers to a 2DEG if two-dimensional dislocation networks exist in the channel of MOSFETs (Fig. 2). Because a network is formed by periodically spaced dislocations, the 2DEG represents a two-dimensional set of one-dimensional electron gases (1DEG). Moreover, the energy quantization is related to the presence of Coulomb blockades proved by low-temperature measurements. The model of dislocations as quantum wires bears similarities to the previously published shallow band model, which, however, assumes 1D-bands about 50-100 meV below the conduction and above the valence band, respectively [8,36,37]. Deeper defect-induced states of minority or majority carriers proved by DLTS [6,11,38] are not considered therein, but are explainable on the basis of the quantum wire model. If atomic defects on the dislocation core are taken into account (kinks, jogs, etc. [25]) then local changes in strain along the dislocation line exist modifying locally the band structure and consequently the quantum wire. Therefore, a broad defect line appears in the DLTS spectrum.
An analogous discussion of the formation of quantum wells by dislocation-induced compressive strain appears also possible which, however, is based on the primary indirect transition between C and D.

Conclusions
Detrimental effects on electronic parameters have been stated for dislocations if they are randomly placed in active areas of semiconductor devices. Presented experimental data, however, show a significant increase in device performance (especially of the drain current), if defined numbers and types of pure dislocations are placed in the channel of MOSFETs and directly connected to source and drain. The drain current enhancement is induced by exceptional electronic properties of dislocations, namely by supermetallic conductivity.
The presented model of the electronic properties of dislocations is based on the high strain in the defect core. It is shown that the high strain level on the dislocation core in the order of 10 % or more modifies the electronic band structure resulting in the formation of a quantum wire inside the silicon channel. The lateral size of the quantum wire is effected by the core diameter of the dislocation and is in the order of a few nanometers. The formation of the quantum wire forces carrier confinement along the dislocation line and generates a 1DEG. Most important for the exceptional properties is the energy quantization instead of the effective mass. Assuming, for instance, a wire radius R = 2.5 nm, the lowest energy level is E 1 = 38.4 meV (Eq. 5) which corresponds to experimentally measured distances of Coulomb blockades [20]. The utilization of the exceptional electronic properties of dislocations enables new concepts in development of future device generations. The main advantages in integration of defined numbers and types of dislocations are (1) further application of silicon-based CMOS processes below 10 nm feature size, (2) easy production processes of dislocation arrangements also in very thin SOI layers, (3) their thermal stability over process steps, and (4) the use of native nanostructures without improved lithography and additional non-silicon materials.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.