High-Density Well-Aligned Dislocations Introduced by Plastic Deformation in Bi 1 − x Sb x Topological Insulator Single Crystals

: Topological insulators (TIs) have a bulk bandgap and gapless edge or surface states that host helically spin-polarized Dirac fermions. Theoretically, it has been predicted that gapless states could also be formed along dislocations in TIs. Recently, conductivity measurements on plastically deformed bismuth antimony (Bi 1 − x Sb x ) TIs have revealed excess conductivity owing to dislocation conduction. For further application of them, fundamental study on dislocations in TIs is indispensable. Dislocations controlled based on fundamental studies could potentially be useful not only for experimental investigations of the dislocation properties but also for diverse device applications. In the present study, Bi 1 − x Sb x TI single crystals were fabricated by a zone-melting method. The crystals were plastically deformed at room temperature. The resultant dislocations were observed by transmission electron microscopy (TEM). It was found that high-density dislocations with the Burgers vector satisfying the condition for the formation of gapless states were successfully introduced. The dislocations were mostly of edge type with lengths on the order of more than a few micrometers.


Introduction
Topological insulators (TIs) are characterized by a bulk bandgap and necessarily highly conducting gapless edge or surface states that host helically spin-polarized Dirac fermions [1,2]. These gapless states are protected by time reversal symmetry. Therefore, they are robust against disorder scattering. The verification of three-dimensional (3D) TIs has attracted major attention to the physical properties of surface gapless states. In addition, it has been theoretically shown that gapless states could also be formed along dislocations in 3D TIs satisfying the following condition [3]: b·M = π (mod 2π), where b is the Burgers vector of the dislocation, (ν 1 , ν 2 , ν 3 ) are the three of the four Z 2 invariants (ν 0 ; ν 1 , ν 2 , ν 3 ) of the 3D TI [4,5], and (G 1 , G 2 , G 3 ) are the reciprocal basis vectors. 3D TIs are characterized by four Z 2 topological invariants (ν 0 ; ν 1 , ν 2 , ν 3 ). The first invariant, ν 0 , defines strong or weak TIs, where ν 0 = 1 corresponds to strong TI (STI) and ν 0 = 0 corresponds to weak TI (WTI). In the case of STI, an odd number of Dirac mode(s) exists on the surface Brillouin zone (BZ) and M, including the latter three invariants, ν 1 , ν 2 , and ν 3 , indicates the location of the Dirac point(s) in the surface BZ. On the other hand, WTI is considered as layers of two-dimensional (2D) TI stacked along M. When the condition in Equation (1) is satisfied, the surface Dirac particles acquire a π phase shift on crossing the dislocation. Consequently, a topologically protected gapless state is formed along the dislocation. Similarly to edge states of 2D TIs, which are considered as guided propagation paths of electrons, the dislocations with gapless states are expected to carry dissipationless current. In addition, at high densities of such dislocations, the thermoelectric figure of merit ZT of TI materials can be dominated by one-dimensional (1D) topologically protected conducting states channeled through the dislocations [6]. Thus, the dislocation with gapless states is a candidate for applications in electronic and thermoelectric devices. It should be noted that, in order to satisfy the condition in Equation (1), at least one of the latter three invariants, ν 1 , ν 2 , or ν 3 , must be 1. However, most of the TIs discovered so far, including some of the most-intensively studied (Bi,Sb) 2 (S,Se,Te) 3 tetradymites, exhibit band inversion at the Γ point; thus, they are classified with the invariants (1;000) [2], where no 1D gapless states along the dislocations can be realized. In contrast, bismuth antimony (Bi 1-x Sb x ) TI exhibit band inversion at the L points, and has the invariants (1;111) [7,8]. Thus, it is one of the limited candidates for realization of 1D gapless states along the dislocations.
A Bi 1−x Sb x solid solution is formed in all compositions of 0 < x < 1. The electronic band structure is drastically modified depending on the composition. It has been known that the solid solution becomes a semiconductor in the composition range 0.07 < x < 0.22 [9]. In 2006, the semiconductor was predicted to be a TI with the invariants (1;111) [7], which was experimentally verified in the following year [8].
As for the conducting dislocations in TIs, while several theoretical studies have been reported [10][11][12], no experimental investigations have been reported until recently [13]. Hamasaki et al. performed conductivity measurements on plastically deformed Bi 1−x Sb x TI, which exhibited excess conductivity owing to dislocation conduction [13]. For further use of them in applications, a fundamental study on dislocations in TIs is indispensable. Dislocations controlled based on fundamental studies could potentially be useful not only for experimental investigations of the dislocation properties, but also for diverse device applications.
Bi 1−x Sb x belongs to the rhombohedral A 7 structure. Concerning dislocations in Bi 1-x Sb x , the report by Hamasaki et al. [13] is the first and the only one so far. Regarding Bi that has the same crystal structure as Bi 1−x Sb x , some studies that used the etch pit method on the dislocations introduced by plastic deformation have been reported [14,15]. However, there are no reports on transmission electron microscopy (TEM) observations of dislocations in Bi, which is mainly due to the difficulty of TEM specimen fabrication, i.e., Bi has a large atomic number; therefore, thinner specimens are required to achieve electron transparency compared to materials with smaller atomic numbers.
In the present study, Bi 1−x Sb x TI single crystals were fabricated by a zone-melting method. The crystals were plastically deformed at room temperature. The resultant dislocations were observed by TEM to reveal the substructure, density, and Burgers vectors.

Materials and Methods
Bi 1−x Sb x single crystals were grown by a zone-melting method. First, a mixture of high-purity (6N) Bi and Sb elements, with an Sb concentration of x = 0.10, was sealed in an evacuated quartz ampoule with a 7 mm inner diameter. Then, the material was melted and homogenized at 400-450 • C, followed by quenching in ice water. Subsequently, the ingot was subjected to zone melting for crystal growth by using a four mirror infrared furnace (FZ-T-10000-H-HR-TUM-VM-PC, Crystal Systems Co., Ltd., Yamanashi, Japan). The growth velocity was set at 0.1 mm/h to avoid constitutional supercooling resulting in the segregation of Sb. The grown crystal was cut in half along the growth direction to obtain a flat surface. Then, the compositions of the grown crystals were investigated by an electron probe microanalyzer (EPMA) (JXA-8800RL, JEOL Ltd., Tokyo, Japan). For quantitative composition analyses, a melt-spun Bi 0.9 Sb 0.1 sample, produced by a single-roller method, was used as a standard for ZAF conversion from the intensities of Bi Lα 1 and Sb Lα 1 to concentrations, where Z is the atomic number factor, A is the absorption factor, and F is the fluorescence factor.
Rectangular samples approximately 1 × 1 × 2 mm 3 in size were cut out and used for the deformation tests. The compression axis and the strain rate were carefully chosen so as to introduce dislocations satisfying the condition in Equation (1). Figure 1 shows the crystallographic orientations of the samples for deformation. The compression axis was set to the 211 direction, in which the primary slip system ( 111) 110 is sufficiently suppressed. According to a previous study that used the etch pit method on the dislocations, which were introduced by plastic deformation of Bi [15], the slip system (011)[ 100] should be activated with the orientation relationship shown in Figure 1, where the Schmid factor for the slip system (011)[ 100] and the primary slip systems ( 111) 110 is 0.46 and 0, respectively.
Here, it should be noted that b = 100 satisfies Equation (1), while b = 110 does not. Although the [ 111]. direction is another possible compression axis to activate the slip system {011} 100 [13,15], the 211 direction is advantageous to activate a single slip system (011)[ 100]. In our previous study [13], the compression axis was the direction declined by a few degrees from the [ 111] direction toward the 111 direction, where a single slip system (011)[ 100] out of the three equivalent slip systems of {011} 100 is activated. However, with the compression axis, as the deformation proceeds, the possibility that other two equivalent slip systems, ( 101)[010] and ( 110)[001], are activated increases. Therefore, we consider that the 211 direction is better than the [ 111] direction to activate a single slip system out of the three equivalent slip systems of {011} 100 . The strain rate was set to be approximately 1 × 10 −5 s −1 to avoid twin formation. Twinning deformation easily occurs at a strain rate of 1 × 10 −3 s −1 [16,17] and is suppressed at a strain rate of 1 × 10 −4 s −1 or lower [15,18]. The rectangular samples were deformed by compression at room temperature up to 10%. Rectangular samples approximately 1 × 1 × 2 mm 3 in size were cut out and used for the deformation tests. The compression axis and the strain rate were carefully chosen so as to introduce dislocations satisfying the condition in Equation (1). Figure 1 shows the crystallographic orientations of the samples for deformation. The compression axis was set to the [21 1 ] direction, in which the primary slip system (111)〈11 0〉 is sufficiently suppressed. According to a previous study that used the etch pit method on the dislocations, which were introduced by plastic deformation of Bi [15], the slip system (011)[100] should be activated with the orientation relationship shown in Figure 1, where the Schmid factor for the slip system (011)[100] and the primary slip systems (111)〈11 0〉 is 0.46 and 0, respectively. Here, it should be noted that = 〈100〉 satisfies Equation (1) , are activated increases. Therefore, we consider that the 〈21 1 〉 direction is better than the [111] direction to activate a single slip system out of the three equivalent slip systems of {011}〈100〉. The strain rate was set to be approximately 1 × 10 −5 s −1 to avoid twin formation. Twinning deformation easily occurs at a strain rate of 1 × 10 −3 s −1 [16,17] and is suppressed at a strain rate of 1 × 10 −4 s −1 or lower [15,18]. The rectangular samples were deformed by compression at room temperature up to 10%. The deformed samples were examined by TEM (JEM-2010F, 2010HC, 200 kV, JEOL Ltd., Tokyo, Japan) to identify the substructure, density, and the Burgers vector of the introduced dislocations. The TEM specimens were cut out along the (011) plane and the (111) plane of the deformed samples, and these were further prepared by mechanical thinning and Ar ion milling.

Crystal Growth
The grown crystal, 27.5 mm in length and 7 mm in diameter, contained two grain boundaries along the growth direction. Figure 2a,b shows the EPMA line profiles of Sb concentrations obtained for the grown crystal. The line profile in Figure 2a was obtained along the crystal growth direction. In the region between 4.5 mm and 17.5 mm from the tip, the crystal was within the range of the TI phase, 0.07 < x < 0.22. The line profile in Figure 2b was perpendicular to the growth direction at a position of 13.5 mm from the tip of the grown crystal. Obvious jumps in the profile corresponded to cavities around the grain boundaries. The radial variation in the Sb concentration was about 2%. Sb concentrations were lower in the center of the crystal than in its outer regions, which could be due to the radial temperature gradient during crystal growth. The deformed samples were examined by TEM (JEM-2010F, 2010HC, 200 kV, JEOL Ltd., Tokyo, Japan) to identify the substructure, density, and the Burgers vector of the introduced dislocations. The TEM specimens were cut out along the (011) plane and the ( 111) plane of the deformed samples, and these were further prepared by mechanical thinning and Ar ion milling.

Crystal Growth
The grown crystal, 27.5 mm in length and 7 mm in diameter, contained two grain boundaries along the growth direction. Figure 2a,b shows the EPMA line profiles of Sb concentrations obtained for the grown crystal. The line profile in Figure 2a was obtained along the crystal growth direction. In the region between 4.5 mm and 17.5 mm from the tip, the crystal was within the range of the TI phase, 0.07 < x < 0.22. The line profile in Figure 2b was perpendicular to the growth direction at a position of 13.5 mm from the tip of the grown crystal. Obvious jumps in the profile corresponded to cavities around the grain boundaries. The radial variation in the Sb concentration was about 2%.
Sb concentrations were lower in the center of the crystal than in its outer regions, which could be due to the radial temperature gradient during crystal growth.  Figure 3 shows the stress-strain curve obtained from the sample with an Sb concentration in the range of 8.6-9.8 at%. The yield stress was about 60 MPa. This value was an order of magnitude higher than that for the deformation of Bi at room temperature [15,[18][19][20]. This result is consistent with the fact that the interatomic distance of the Bi1-xSbx solid solution becomes shorter with increasing Sb concentrations. Work hardening behavior was slightly observed during the initial stage of plastic deformation. Slip traces parallel to the [100] direction were observed on the (01 1) surface of the deformed samples. In our previous study [13], where the compression axis declined a few degrees from the [111] direction, obvious work hardening was observed, indicating that possibly more than one equivalent slip systems were activated.

Plastic Deformation
.   Figure 3 shows the stress-strain curve obtained from the sample with an Sb concentration in the range of 8.6-9.8 at%. The yield stress was about 60 MPa. This value was an order of magnitude higher than that for the deformation of Bi at room temperature [15,[18][19][20]. This result is consistent with the fact that the interatomic distance of the Bi 1-x Sb x solid solution becomes shorter with increasing Sb concentrations. Work hardening behavior was slightly observed during the initial stage of plastic deformation. Slip traces parallel to the [ 100] direction were observed on the 011 surface of the deformed samples. In our previous study [13], where the compression axis declined a few degrees from the [ 111] direction, obvious work hardening was observed, indicating that possibly more than one equivalent slip systems were activated.  Figure 3 shows the stress-strain curve obtained from the sample with an Sb concentration in the range of 8.6-9.8 at%. The yield stress was about 60 MPa. This value was an order of magnitude higher than that for the deformation of Bi at room temperature [15,[18][19][20]. This result is consistent with the fact that the interatomic distance of the Bi1-xSbx solid solution becomes shorter with increasing Sb concentrations. Work hardening behavior was slightly observed during the initial stage of plastic deformation. Slip traces parallel to the [100] direction were observed on the (01 1) surface of the deformed samples. In our previous study [13], where the compression axis declined a few degrees from the [111] direction, obvious work hardening was observed, indicating that possibly more than one equivalent slip systems were activated.

Plastic Deformation
. Figure 3. Stress-strain curve of the bismuth antimony (Bi1−xSbx) sample deformed at room temperature under a constant strain rate of 1 × 10 −5 s −1 . The Sb concentration of the deformed sample was in the range of 8.6-9.8 at%.   Figure 4a shows typical TEM bright-field images obtained from the (011) slip plane of the deformed sample. The observation direction was tilted from the normal direction of the (011) slip plane by approximately 10

TEM Observations of Dislocations Introduced by Plastic Deformation
• . Using the convergent-beam technique, the thickness of the TEM specimen was determined to be approximately 100 nm. The total dislocation length observed in the area of 12 µm 2 in Figure 4a was approximately 50 µm. The dislocation density in that area was estimated to be 4 × 10 9 cm − 2 . From the extensive TEM observations, the average dislocation density was estimated to be 4 − 8 × 10 9 cm − 2 . We found that dislocations in the present study were introduced more uniformly throughout the TEM specimen compared with our previous study [13]. The dislocations had long segments along the 011 direction. Figure 4b shows the same images as in Figure 4a, where dislocations along the 011 direction are indicated by red lines, while dislocations that deviate from the 011 direction more than 10 • are indicated by blue lines. It was confirmed by tilting the specimen that other contrasts were not dislocations but bend contours. Figure 4c shows a stereographic projection of   For g·b analyses to investigate the Burgers vector, TEM observations were performed for the ( 111) plane. Figure 5a shows a typical TEM bright-field image obtained from the ( 111) slip plane of the deformed sample. Since the intersection of the ( 111) plane and the (011) slip plane was parallel to the 011 direction, straight dislocations along the 011 direction were observed. Figure 5b,c shows TEM dark-field images obtained from the same region as in Figure 5a, using g = 211 and g = 011, respectively. As a result of g·b analyses, the contrast of the dislocations vanished with g = 011. Thus, the Burgers vector was considered to be perpendicular to the 011 direction.
The primitive translation vector that is perpendicular to the 011 direction in the (011) slip plane is [ 100]. Consequently, the Burgers vector of the introduced dislocations was determined to be b = [ 100], which satisfies Equation (1). Dislocations with other Burgers vectors were not found in our TEM observations. This result suggests that longer dislocations with controlled orientations can be obtained by mild plastic deformation without severe work hardening. Since the dislocations had long segments along the 011 direction, which is perpendicular to the Burgers vector, the observed dislocations were mostly of the edge type. The length of the edge components was approximately a few micrometers. Here, it should be noted that the actual dislocation length could be longer than the observed one because of the difference between the (011) slip plane and the surface of the TEM specimen; the dislocations on the (011) slip plane were inevitably intersected by the surface of the TEM specimen.
Crystals 2019, 9, x FOR PEER REVIEW 7 of 10 orientations can be obtained by mild plastic deformation without severe work hardening. Since the dislocations had long segments along the [01 1] direction, which is perpendicular to the Burgers vector, the observed dislocations were mostly of the edge type. The length of the edge components was approximately a few micrometers. Here, it should be noted that the actual dislocation length could be longer than the observed one because of the difference between the (011) slip plane and the surface of the TEM specimen; the dislocations on the (011) slip plane were inevitably intersected by the surface of the TEM specimen. The finding that the dislocations were mostly edge dislocations could be interpreted in terms of the Peierls stress. Bi1−xSbx with a rhombohedral crystal structure has a cleavage on the (111) plane. Figure 6a  The finding that the dislocations were mostly edge dislocations could be interpreted in terms of the Peierls stress. Bi 1−x Sb x with a rhombohedral crystal structure has a cleavage on the ( 111) plane. Figure 6a,b shows the atomic structure of Bi 1−x Sb x viewed along the 011 and 111 directions, respectively. The 111 direction is 3 • tilted from the normal direction of the (011) slip plane toward the [ 100] direction. Atomic planes parallel to the ( 111) plane are arranged in double layers as shown in Figure 6a. The interatomic distance in a double layer is shorter than that between adjacent double layers [21]. Therefore, in the (011) plane, the interatomic distance along the [ 100] direction is larger than that along the 011 direction. In addition, it can be assumed that the shorter bonds are stronger and the longer bonds are weaker. Thus, the Peierls stress along the [ 100] direction should be higher than that along the 011 direction. Accordingly, it can be considered that, for mobile dislocations in the (011) plane, the edge segments along the 011 direction, which are trapped along the valley of the Peierls potential, require higher stress to glide, thus leading to their lower mobility. In addition, dislocation segments along the 111 and 111 directions observed in Figure 4 are also expected to be stable because they are along the atomic columns as shown in Figure 6b. Resultant long dislocation segments are expected to have an advantage in dislocation conduction measurements. Furthermore, such long dislocations aligned in one direction may be useful for future device applications.  Figure 6a. The interatomic distance in a double layer is shorter than that between adjacent double layers [21]. Therefore, in the (011) Figure 4 are also expected to be stable because they are along the atomic columns as shown in Figure 6b. Resultant long dislocation segments are expected to have an advantage in dislocation conduction measurements. Furthermore, such long dislocations aligned in

Conclusions
Bi1−xSbx single crystals, within the range of the TI phase with the invariants (1;111), were fabricated by a zone-melting method. The crystals were plastically deformed at room temperature. The compression axis was set at the [21 1 ] direction to activate the slip system (011)[100]. As a result, high-density dislocations with a Burgers vector = [100], satisfying the condition for the formation of gapless states, were introduced. The dislocations were mostly of edge type with lengths on the order of more than a few micrometers. For mobile dislocations in the (011) plane, it

Conclusions
Bi 1−x Sb x single crystals, within the range of the TI phase with the invariants (1;111), were fabricated by a zone-melting method. The crystals were plastically deformed at room temperature.
The compression axis was set at the 211 direction to activate the slip system (011)[ 100]. As a result, high-density dislocations with a Burgers vector b = [ 100], satisfying the condition for the formation of gapless states, were introduced. The dislocations were mostly of edge type with lengths on the order of more than a few micrometers. For mobile dislocations in the (011) plane, it can be considered that the edge segments along the 011 direction require higher stress to glide, thus leading to their lower mobility. Resultant long dislocation segments aligned in one direction are expected to have an advantage in experimental investigations of the dislocation properties and future device fabrications.