Abstract
Micropillar compression tests of Ti5Si3 single crystals were conducted at room temperature as a function of loading axis orientation and specimen size in order to investigate their room temperature plastic deformation behavior. Plastic flow by the operation of three deformation modes, {1\({\overline{1}}\)00}[0001], {2\({\overline{1}}\)\({\overline{1}}\)2} < 2\({\overline{1}}\)\({\overline{1}}\)\({\overline{3}}\) > and {1\({\overline{1}}\)01} < 2\({\overline{1}}\)\({\overline{1}}\)\({\overline{3}}\) > slip were observed in [2\({\overline{2}}\)05]-, [0001]- and [4\({\overline{3}}\)\({\overline{1}}\)0]-oriented micropillar specimens deformed at room temperature, respectively. The CRSS values were evaluated to be very high above 2.7 GPa and were confirmed to increase up to about 6 GPa with the decrease in the specimen size. The fracture toughness values are evaluated to be 0.45 MPa m1/2 (notch plane // (0001)) and 0.73 MPa m1/2 (notch plane //(1\({\overline{1}}\)00)) based on the results of micro-cantilever bend tests of chevron-notched specimens. The fracture toughness values are considerably lower than those for D8l-Mo5SiB2 and D8l-Nb5Si3 evaluated by the same method, indicating the inherent brittleness of binary Ti5Si3 compared to the other transition-metal silicides of the TM5Si3 type (TM: transition-metal).
Similar content being viewed by others
Introduction
Since the gas inlet temperature of the advanced turbine systems already exceeds 1700 °C that is much higher than the melting temperature (~ 1350 °C) of Ni-based superalloys, there is an ever-increasing demand for developing ultra high-temperature structural materials that can withstand at temperatures well above the maximum operating temperatures of Ni-base superalloys in severe oxidizing atmosphere1,2. Silicides with the chemical formula of TM5Si3 formed with transition-metals (TMs) are of interest, as most of them exhibit the highest melting temperature in each of the corresponding TM-Si binary phase diagrams and a reasonably good oxidation property due to the high Si content3. Indeed, Nb5Si3 and Mo5SiB2 with the tetragonal D8l structure and Mo5Si3 with the tetragonal D8m structure have intensively been investigated in the last several decades as the strengthening phases of Nb-4,5,6,7 and Mo-8,9,10 based refractory alloys, respectively. Ti5Si3 with the hexagonal D88 structure, on the other hand, has been investigated initially as a constituent phase of Ti-Ti5Si3 in-situ composites11 and in recent years as the strengthening phase of Mo-Mo5SiB2 based alloys because the creep strength and oxidation resistance are significantly improved by replacing a large amount of Mo with Ti so as to incorporate Ti5Si3 as the constituent phase of Mo-Mo5SiB2 based alloys12,13,14,15,16,17. Because of some excellent properties as a promising structural material for ultra-high temperature applications [the high melting temperature (2130 °C) and low density (4.32 g/cm3)], the deformation behavior of monolithic Ti5Si3 has been investigated with polycrystals11 and single crystals18,19. Frommeyer et al.11 reported that Ti5Si3 is extremely brittle at low temperatures as expected from the complex crystal structure of the D88 type (Pearson symbol: hP16, Space group: P63/mcm). Plastic flow is observed only above 1000 °C in polycrystals, while Umakoshi and Nakashima18 reported plastic flow carried by deformation twinning of the {1\(\overline{1}\)02} < \(\overline{1}\)101 > -type at high temperatures above 1300 °C in single crystals with some limited orientations. Then, Kishida et al.19 made a first systematic deformation experiment on single crystals of Ti5Si3 in uniaxial compression and found that plastic flow is possible only above 1200 °C regardless of crystal orientation with three different operative deformation modes, {1\(\overline{1}\)00} < 0001 > prismatic slip, {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip and {2\(\overline{1}\)\(\overline{1}\)8} < 8\(\overline{4}\)\(\overline{4}\)\(\overline{3}\) > twinning depending on crystal orientation.
Recently, micropillar compression testing that utilizes small specimens of the micron-meter scale has been recognized as a powerful method to investigate the plastic deformation behavior of hard and brittle materials at a low temperature (for example, room temperature) well below their onset temperatures (above 1000 °C for most of these materials) for plastic flow in the bulk form. Those hard and brittle materials include intermetallics with complex crystal structures20,21,22,23,24,25, semiconductors26,27,28,29,30,31 and ceramics32,33,34,35. We have investigated the room-temperature plastic deformation behavior of various hard and brittle materials including D8l-Mo5SiB224, D8l-Nb5Si325, 6H-SiC31 and so on21,22,23 and have successfully identified operative slip systems and their critical resolved shear stress (CRSS). Of particular interest to note is that slip systems that do not participate in the high temperature deformation in the bulk form do participate in the room-temperature deformation in the micropillar form24,25,31 and that completely different dislocation dissociation schemes are observed at high (bulk) and low (micropillar) temperatures even when the same slip system operates. Thus, the CRSS at a low temperature obtained in the micropillar form is plotted far away from the extrapolation of the temperature-dependent CRSSs obtained at high temperatures in the bulk form31. Such information about the room-temperature deformation behavior must be very important to conceive a strategy to improve the brittleness at low temperatures of these hard and brittle materials.
In the present study, we investigate the room-temperature deformation behavior of single crystals of Ti5Si3 with the hexagonal D88 structure as a function of specimen size and loading-axis orientation by micropillar compression tests, to get information that might be useful for the improvement of the brittleness at low temperatures.
Results
Compression deformation behavior
Figure 1a–c show typical stress–strain curves obtained for single crystals of micropillar specimens with the [2\(\overline{2}\)05], [4\(\overline{3}\)\(\overline{1}\)0] and [0001] orientations, respectively. A relatively large flat portion corresponding to a strain burst is commonly observed soon after elastic loading in many stress–strain curves regardless of the loading-axis orientation. Although a strain burst is followed by a load drop in most cases, that is not included in Fig. 1. As observed in other micropillar compression tests of single crystals of metallic materials36,37,38,39,40, the occurrence of a strain burst in stress–strain curves in Ti5Si3 is considered to be plastic flow caused by dislocation motion in avalanches, as in the case of Mo5SiB224,37,40.
Being consistent with the occurrence of plastic flow caused by dislocation motion in avalanches, most micropillar specimens exhibit slip plane (shear) failure, accompanied by giant steps on the slip plane along the slip direction (Fig. 2). Deformation microstructures with the three orientations are shown in Fig. 2a–c for some specimens in which micropillar testing was successfully interrupted before instantaneous slip plane failure occurs. For the [2\(\overline{2}\)05] orientation (Fig. 2a), analysis of slip traces that appeared on two orthogonal side-surfaces has identified the slip plane to be (1\(\overline{1}\)00). Because shear deformation occurs so as to maintain the (\(\overline{1}\overline{1}\)20) side-surface flat, the shear direction is contained in the (\(\overline{1}\)\(\overline{1}\)20) plane and is determined to be [0001]. The slip system operative in [2\(\overline{2}\)05]-oriented micropillar specimens is thus identified to be (1\(\overline{1}\)00)[0001], as observed in the same orientation of bulk single crystals at above 1300 °C19.
For the [4\(\overline{3}\)\(\overline{1}\)0] orientation, deformation markings corresponding to slip on (1\(\overline{1}\)01) are observed on two orthogonal side faces. On the assumption that slip on this plane occurs along either a or a + c direction (when expressed with the hexagonal lattice parameters), the possible slip vectors in the [4\(\overline{3}\)\(\overline{1}\)0] orientation are [11\(\overline{2}\)0], [2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\)] and [1\(\overline{2}\)1\(\overline{3}\)], the latter two of which have the identical Schmid factors. Of the three, [11\(\overline{2}\)0] is ruled out, since deformation markings are clearly observed also on the (000\(\overline{1}\)) side-surface that contains the [11\(\overline{2}\)0] direction (Fig. 2b). The slip system identified in [4\(\overline{3}\)\(\overline{1}\)0]-oriented micropillar specimens is thus identified to be {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > . This slip system has never been identified to operate in bulk single crystals. In bulk single crystals with the same [4\(\overline{3}\)\(\overline{1}\)0] orientation, {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip was identified to operate above 1400 °C19.
For the [0001] orientation, slip trace analysis reveals that slip occurs on the (11\(\overline{2}\)2) slip plane. Because 1/3 < \(\overline{1}\)\(\overline{1}\)2\(\overline{3}\) > is the shortest lattice translation vector (0.906 nm) on {11\(\overline{2}\)2} slip planes, the operative slip system in [0001]-oriented micropillar specimens is inferred to be {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > . As described above, {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip has been observed in [4\(\overline{3}\)\(\overline{1}\)0]-oriented bulk single crystals above 1400 °C but not in [0001]-oriented bulk single crystals, in which deformation twinning of the {2\(\overline{1}\)\(\overline{1}\)8} < 8\(\overline{4}\)\(\overline{4}\)\(\overline{3}\) > -type is operative above 1400 °C19.
Fracture toughness
Micro-cantilever bend tests on chevron-notched micro-beam specimens were conducted to evaluate the fracture toughness. The load–displacement curves obtained for chevron-notched micro-beam specimens with a notch plane parallel to (0001) and (1\(\overline{1}\)00) are shown in Fig. 3a. The loading directions were set parallel to [1\(\overline{1}\)00] and [0001] for (0001)- and (1\(\overline{1}\)00)-notched micro-beam specimens, respectively. The load–displacement curves (Fig. 3a) exhibit no sign of apparent plastic deformation before failure occurs. In addition, a fairly flat fracture surface observed in the SEM image for the (0001)-notched micro-beam specimen indicates the occurrence of brittle cleavage fracture (Fig. 3b). The values of fracture toughness can then be evaluated simply from the maximum load and geometrical parameters of the chevron-notched specimens based on the model proposed by Deng et al.41, details of which were described in our previous papers24. The fracture toughness values of Ti5Si3 are evaluated to be 0.45 and 0.73 MPa m1/2 respectively with the notch plane parallel to (0001) and (1\(\overline{1}\)00), both of which are considerably lower than those for D8l-Mo5SiB2 (2.43 MPa m1/2; notch plane // (100)) and D8l-Nb5Si3 (1.79 MPa m1/2; notch plane // (001)) evaluated by the same method24,25.
Discussion
The CRSS values (τCRSS) for the three slip systems identified in the present study are calculated with yield stresses defined as the stress at which the first strain burst occurs and the corresponding Schimid factors as plotted in Fig. 4a as a function of specimen size (the edge length L). The CRSS values for all slip systems are extremely high exceeding 2.6 GPa regardless of specimen size and exhibit the so-called “smaller is stronger” trends, following an inverse power-law relationship, τCRSS ∝ L−n (n: the power-law exponent), as commonly observed in micropillar compression of single crystals of many metallic materials, ceramics, semiconductors and intermetallic compounds22,23,24,25,30,31,32,33,35,36,38,39,42,43. The power-law exponents for the three slip systems of {1\(\overline{1}\)00}[0001], {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > and {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > are estimated to be 0.36, 0.30 and 0.11, respectively. The power-law exponents in hard and brittle materials such as 6H-SiC (n = 0.10–0.21), Mo5SiB2 (n = 0.16–0.22) are generally much lower than those reported for face-centered cubic (FCC) metals (0.5–1.0) and body-centered cubic (BCC) metals (0.2–0.5)38,39,42. However, the power-law exponents observed for the two slip systems (n = 0.36 and 0.30 for {1\(\overline{1}\)00}[0001] and {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip) in Ti5Si3 are not as low as those for 6H-SiC and Mo5SiB2 but are comparably high as those observed for some particular slip systems (0.43 and 0.38 for (001) < 010 > and {110} < 1\(\overline{1}\)0 > slip) in α-Nb5Si325. In micropillar compression experiments of single crystals of FCC and BCC metals, the inverse power-low curve of size-dependent CRSS is known to coincide with the bulk CRSS when the specimen size is in the range of 20–30 μm43. If we assume that the same criteria is applicable to Ti5Si3, the bulk CRSS values for {1\(\overline{1}\)00}[0001], {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > and {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > are estimated as the first approximation respectively to be 1.43 ± 0.10, 1.90 ± 0.10 and, 2.32 ± 0.06 GPa, from the extrapolation of the inverse power-low curves.
Figure 4b compares the estimated bulk CRSS values at room temperature for the three slip systems, {1\(\overline{1}\)00}[0001], {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > , and {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > with those obtained in high-temperature compression experiments of bulk single crystals. In bulk compression experiments, we have identified three operative deformation modes, {1\(\overline{1}\)00}[0001] prism slip, {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip and {2\(\overline{1}\)\(\overline{1}\)8} < 8\(\overline{4}\)\(\overline{4}\)\(\overline{3}\) > twinning and have confirmed that the CRSS for these three deformation modes strongly depend on temperature. For a constant strain-rate experiments, a linear correlation is generally observed between ln τCRSS and (1/T), since the strain rate \(\dot{\gamma }\) can be expressed as a function of CRSS, temperature T and activation energy Q for thermal activation process for dislocation glide with the following equation.
where m is the stress exponent and k is the Boltzmann constant44. Apparently, the room-temperature bulk CRSS values estimated for {1\(\overline{1}\)00}[0001] prism slip, {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip are much lower than those expected from the extrapolation of the ln τCRSS − (1/T) relations for the deformation modes operative at high temperatures. This suggests that the deformation mechanism that governs the CRSS for {1\(\overline{1}\)00}[0001] prism slip and {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip varies depending on testing temperature. One possible explanation for the different deformation mechanisms operating at high and low temperatures is that the room-temperature bulk CRSSs estimated in the micropillar form is determined by the nucleation of new dislocations from specimen surfaces of micropillars, as we discussed previously for 6H-SiC31 and α-Nb5Si325 because virtually no grown-in dislocation is expected to exist in micropillar specimens of Ti5Si3 single crystals when judged from the brittleness and high CRSSs. However, more detailed studies including dislocation analysis by TEM are definitely needed to prove this.
The expected orientation-dependent operative slip systems are calculated based on the estimated bulk CRSS values for the three slip systems identified in micropillar compression tests and are shown in the stereographic projection of Fig. 4c. As expected from the relatively low CRSS values for {1\(\overline{1}\)00}[0001] and {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip, the orientation range for these two slip systems to operate are very wide, while that for {1\(\overline{1}\)01} < 2\(\overline{1}\overline{1}\)\(\overline{3}\) > slip is quite limited to near [1\(\overline{1}\)00]. A sufficient number of independent slip systems is achieved for general deformation of polycrystalline aggregates by the combination of {1\(\overline{1}\)00}[0001] and {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip, both of which have a relatively smaller CRSS values (1.43 ± 0.10 and 1.90 ± 0.10 MPa) when compared to {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip (2.32 ± 0.06 GPa).
For the [0001] orientation, {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip is found to operate in micropillar specimens at room temperature, while {2\(\overline{1}\)\(\overline{1}\)8} < 8\(\overline{4}\)\(\overline{4}\)\(\overline{3}\) > twinning is operative in bulk single crystals at high temperatures above 1300°C19. The absence of {2\(\overline{1}\)\(\overline{1}\)8} < 8\(\overline{4}\)\(\overline{4}\)\(\overline{3}\) > twinning at room temperature is easily understood from the fact that {2\(\overline{1}\)\(\overline{1}\)8} < 8\(\overline{4}\)\(\overline{4}\)\(\overline{3}\) > twinning requires very complicated atomic shuffling to restore the crystal structure during deformation twinning, so that it can happen only at sufficiently high temperatures at which sufficient atomic diffusion occurs19. That must be the reason why the [0001]-oriented micropillar specimens need to find out alternative slip system, {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip at room temperature.
For the [4\(\overline{3}\)\(\overline{1}\)0] orientation, on the other hand, while {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip is observed in bulk single crystals at high temperatures, this slip was replaced by {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > pyramidal slip at room temperature for micropillar specimens. In order to see the relative ease in the operation for the slip along [2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\)] (a + c slip) on two different pyramidal plane ((2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01)), generalized stacking fault energy (GSFE) for < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip on (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) were calculated by first-principles DFT calculations. Figure 5 plots the GSFE curves and their derivatives (the gradient of the GSFE curve corresponding to the ideal shear strength, τth) for the two pyramidal slip systems of (2\(\overline{1}\)\(\overline{1}\)2)[2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\)] and (1\(\overline{1}\)01)[2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\)]. For (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) slip planes, there exist four and three different possible slip planes to be considered, respectively (see, the insets of Fig. 5a). Among the possible slip planes, only the calculation results for the slip planes exhibiting the lowest ideal shear strength and unstable stacking fault energies are indicated in Fig. 5. The selected slip planes are indicated with red and blue lines respectively, while the others are indicated with dashed black lines in the insets of Fig. 5a. To be noted in Fig. 5 is that the GSFE curves and their derivative curves are asymmetric on both slip planes with respect to the [2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\)] slip direction and that the displacement along the positive direction gives rise to compression strain for the [0001] orientation while that along the negative direction does so for the [4\(\overline{3}\)\(\overline{1}\)0] orientation. Table 1 summarizes the results of the GSFE calculations. For both orientations, the ideal shear strength and unstable stacking fault energy are both lower on (2\(\overline{1}\)\(\overline{1}\)2) than on (1\(\overline{1}\)01), suggesting the preference for the (2\(\overline{1}\)\(\overline{1}\)2) slip plane for both orientations. This situation does not change even the ideal shear strength is normalized to the corresponding Schmid factors for each loading orientation (the normalized ideal shear strength, τth/m in Table 1). The results of the GSFE calculations is consistent with what is observed in micropillar compression for the [0001] *orientation (the preference of the (2\(\overline{1}\)\(\overline{1}\)2) slip plane) but not with what is observed for the [4\(\overline{3}\)\(\overline{1}\)0] orientation. For the [4\(\overline{3}\)\(\overline{1}\)0] orientation, however, the normalized ideal shear strengths (τth/m) on (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) slip planes differ from each other only by 3% (34.1 and 35.1 GPa, respectively). We believe that the difference in τth/m is too small to draw a final conclusion for the slip plane preference for this orientation. One more thing to be noted in Table 1 is that the magnitude correlation for the ideal shear strength (τth) for the (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) slip planes (22.5 and 16.8 GPa) is opposite to that for the bulk CRSS values deduced from micropillar experiment (1.90 and 2.32 GPa). We believe that this discrepancy comes from the fact that no dislocation dissociations, especially those accompanied by non-collinear Burgers vectors, are taken into account in the above discussion based on the GSFE calculations. In view of the large magnitude (about 0.906 nm) of the Burgers vectors of 1/3 < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > , the 1/3 < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > dislocation is expected to dissociate into partial dislocations with shorter Burgers vectors. In order to deduce the possible dissociation scheme for 1/3 < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > dislocation, construction of full gamma surfaces on (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) planes must be made, because experimental identification of the dissociation scheme is extremely difficult with the occurrence of instantaneous shear failure in micropillar specimens. The construction of gamma surface is currently in progress in the authors’ research group to draw the final conclusion on the slip plane preference.
Methods
A single crystal rod of Ti5Si3 was grown from an ingot with a nominal composition of Ti − 37.5 at.% Si by directional solidification with an optical floating-zone furnace at a growth rate of 6 mm/h. Characteristics in the atomic arrangement and crystal structure of Ti5Si3 are detailed in our previous paper19. The crystal orientation of the grown single crystal was determined by the back-Laue X-ray diffraction method. Three loading-axis orientations, [2\(\overline{2}\)05], [4\(\overline{3}\)\(\overline{1}\)0] and [0001] were selected for micropillar compression tests. The highest Schmid factors for possible slip/twinning systems for each orientation are listed in Table 1 of our previous paper19. After the mirror finishing of specimen surface by mechanical polishing with diamond paste, micropillar specimens with a square cross-section having an edge length L ranging from 0.5 to 3.6 μm and a specimen height-to-L ratio of 1:3–1:5 were fabricated by the focused-ion beam (FIB) machining with a JEOL JIB-4000 FIB milling and imaging system operated at 30 kV. Micropillar compression tests were conducted at room temperature under the displacement-rate controlled mode (a nominal strain rate of 1 × 10–4 s−1) with an Agilent Technologies Nano Indenter G200 nanomechanical tester equipped with a flat-punch diamond indenter tip. Microstructures of micropillar specimens before and after the compression tests were investigated by scanning electron microscopy (SEM) in order to measure dimensions and analyze slip traces, respectively. SEM observations were made along the direction inclined by 30° from the loading axis. Single cantilever bend tests were conducted with chevron-notched micro-beam specimens to evaluate the fracture toughness, as the details of the experimental procedures are described in our previous paper24.
Generalized stacking fault energy (GSFE) for 1/3 < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip on (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) were calculated by first-principles DFT calculations using the Vienna ab-initio simulation package (VASP) code45,46,47. The generalized gradient approximation of Perdew–Burke–Ernzerhof (GGA-PBE) was used to treat the exchange–correlation functional48. We used supercells for the GSFE calculations for (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) planes respectively containing 192 and 96 atoms, in each of which the Ti:Si atomic ratio was fixed at 5:3 so as to maintain the stoichiometry of Ti5Si3. The in-plane units of the supercells were defined by unit vectors along [2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\)] and [01\(\overline{1}\)0] for the (2\(\overline{1}\)\(\overline{1}\)2) supercell and by those along [2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\)] and [1\(\overline{2}\)13] for the (1\(\overline{1}\)01) supercell. A vacuum layer with a total thickness of 1.5 nm (when measured along the slip plane normal) is included in each supercell. Monkhorst–Pack k-point meshes of 10 × 8 × 2 and 8 × 8 × 2 were used for the (2\(\overline{1}\)\(\overline{1}\)2) and (1\(\overline{1}\)01) supercells, respectively49. An energy cutoff of 500 eV was used throughout the calculations. All atoms were relaxed only along the direction perpendicular to the glide plane so as to minimize the energy of the supercell with a given in-plane displacement24.
Conclusions
The room-temperature deformation behavior of Ti5Si3 was investigated by micropillar compression tests as a function of loading-axis orientation and specimen size and micro-cantilever bend tests with chevron-notched specimens. The results obtained are summarized as follows.
-
1.
Three slip systems, {1\(\overline{1}\)00}[0001], {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > and {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > , are operative at room temperature. The CRSS values for all three slip systems increase with the decrease in specimen size in the range of 2.7–6 GPa, following an inverse power-law relationship with the power-law exponents of 0.36, 0.30 and 0.11 for {1\(\overline{1}\)00}[0001], {2\(\overline{1}\)\(\overline{1}\)2} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > and {1\(\overline{1}\)01} < 2\(\overline{1}\)\(\overline{1}\)\(\overline{3}\) > slip, respectively. The bulk CRSS values for these three slip systems are estimated to be 1.43 ± 0.10, 1.90 ± 0.10 and, 2.32 ± 0.06 GPa by extrapolating the inverse power-low curves to the specimen size range of 20–30 μm.
-
2.
The fracture toughness values are evaluated by micro-cantilever bend tests to be 0.45 and 0.73 MPa m1/2 with chevron-notch planes parallel to (0001) and (1\(\overline{1}\)00), respectively. The fracture toughness values are considerably lower than those for D8l-Mo5SiB2 (2.43 MPa m1/2; notch plane // (100)) and D8l-Nb5Si3 (1.79 MPa m1/2; notch plane // (001)) evaluated by the same method, indicating much more significant inherent brittleness of Ti5Si3 compared to the other transition-metal silicides of the TM5Si3 type.
Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
References
Zhao, J.-C. & Westbrook, J. H. Ultrahigh-temperature materials for jet engines. MRS Bull. 28, 622–627 (2003).
Perepezko, J. H. The hotter the engine, the better. Science 326, 1068–1069 (2009).
Mitra, R. Mechanical behaviour and oxidation resistance of structural silicides. Int. Mater. Rev. 51, 13–64 (2006).
Subramanian, P. R., Mendiratta, M. G. & Dimiduk, D. M. The development of Nb-based advanced interemetallic alloys for structural applications. JOM 48, 33–38 (1996).
Jackson, M. R., Bewlay, B. P., Rowe, R. G., Skelly, D. W. & Lipsitt, H. A. High-temperature refractory metal-intermetallic composites. JOM 48, 39–44 (1996).
Bewlay, B. P., Jackson, M. R., Zhao, J. C. & Subramanian, P. R. A review of very-high-temperature Nb-silicide-based composites. Metall. Mater. Trans. 34A, 2043–2052 (2003).
Tsakiropoulos, P. On Nb silicide based alloys: alloy design and selection. Materials 11, 844 (2018).
Akinc, M. et al. Boron-doped molybdenum silicides for structural applications. Mater. Sci. Engng. A 261, 16–23 (1999).
Sakidja, R. & Perepezko, J. H. Phase stability and alloying behavior in the Mo-Si-B system. Metall. Mater. Trans. 36A, 507–514 (2005).
Schneibel, J. H., Ritchie, R. O., Kruzic, J. J. & Tortorelli, P. F. Optimization of Mo-Si-B intermetallic alloys. Metall. Mater. Trans. 36A, 525–531 (2005).
Frommeyer, G., Rosenkranz, R. & Ludecke, C. Microstructure and properties of the refractory intermetallic Ti5Si3 compound and the unidirectionally solidified eutectic Ti–Ti5Si3 alloy. Z. Metallkunde 81, 307–313 (1990).
Yang, Y., Chang, Y. A., Tan, L. & Cao, W. Multiphase equilibria in the meta-rich region of the Mo-Ti-Si-B system: thermodynamic prediction and experimental validation. Acta Mater. 53, 1711–1720 (2005).
Sakidja, R., Perepezko, J. H., Kim, S. & Sekido, N. Phase stability and structural defects in high-temperature Mo-Si-B alloys. Acta Mater. 56, 5223–5244 (2008).
Yang, Y. et al. Effects of Ti, Zr, and Hf on the phase stability of Mo_ss + Mo3Si + Mo5SiB2 alloys at 1600 ̊C. Acta Mater. 58, 541–548 (2010).
Azim, M. A. et al. Effect of Ti (macro-) alloying on the high-temperature oxidation behavior of ternary Mo-Si-B alloys. Oxid. Met. 80, 231–242 (2013).
Schliephake, D. et al. High-temperature creep and oxidation behavior of Mo-Si-B alloys with high Ti contents. Metall. Mater. Trans. 45A, 1102–1111 (2014).
Azim, M. A. et al. Creep resistance and oxidation behavior of novel Mo-Si-B-Ti alloys. JOM 67(11), 2621–2628 (2015).
Umakoshi, Y. & Nakashima, T. High temperature deformation of Ti5Si3 single crystals with D88 structure. Scripta Mater. 30, 1431–1436 (1994).
Kishida, K., Fujiwara, M., Adachi, H., Tanaka, K. & Inui, H. Plastic deformation of single crystals of Ti5Si3 with the hexagonal D88 structure. Acta Mater. 58, 646–857 (2010).
Takata, N., Ghassemi Armaki, H., Terada, Y., Takeyama, M. & Kumar, K. S. Plastic deformation of the C14 Laves phase (Fe, Ni)2Nb. Scripta Mater. 68, 615–618 (2013).
Okamoto, N. L. et al. Compression deformability of Γ and ζ Fe-Zn intermetallics to mitigate detachment of brittle intermetallic coating of galvannealed steels. Scripta Mater. 69, 307–310 (2013).
Okamoto, N. L., Inomoto, M., Adachi, H., Takebayashi, H. & Inui, H. Micropillar compression deformation of single crystals of the intermetallic compound ζ-FeZn13. Acta Mater. 65, 229–239 (2014).
Nakatsuka, S., Kishida, K. & Inui, H. Micropillar compression of MoSi2 single crystals. MRS Symp. Proc. 1760, mrsf14–1760-yy05–09 (2015).
Kishida, K., Maruyama, T., Matsunoshita, H., Fukuyama, T. & Inui, H. Micropillar compression deformation of single crystals of Mo5SiB2 with the tetragonal D8l structure. Acta Mater. 159, 416–428 (2018).
Kishida, K., Maruyama, T., Fukuyama, T. & Inui, H. Micropillar compression of single crystals of α-Nb5Si3 with the D8l structure. submitted to Sci. Technol. Adv. Mater.
Michler, J., Wasmer, K., Meier, S., Östlund, F. & Leifer, K. Plastic deformation of gallium arsenide micropillars under uniaxial compression at room temperature. Appl. Phys. Lett. 90, 043123 (2007).
Moser, B., Wasmer, K., Barbieri, L. & Michler, J. Strength and fracture of Si micropillars: A new scanning electron microscopy-based micro-compression test. J. Mater. Res. 22, 1004–1011 (2007).
Östlund, F. et al. Brittle-to-ductile transition in uniaxial compression of silicon pillars at room temperature. Adv. Funct. Mater. 19, 2439–2444 (2009).
Sung, T. H., Huang, J. C., Hsu, J. H. & Jian, S. R. Mechanical response of GaN film and micropillar under nanoindentation and microcompression. Appl. Phys. Lett. 97, 17190413 (2010).
Östlund, F. et al. Ductile-brittle transition in micropillar compression of GaAs at room temperature. Philos. Mag. 91, 1190–1199 (2011).
Kishida, K., Shinkai, Y. & Inui, H. Room temperature deformation of 6H-SiC single crystals investigated by micropillar compression. Acta Mater. 187, 19–28 (2020).
Korte, S. & Clegg, W. J. Discussion of the dependence of the effect of size on the yield stress in hard materials studied by microcompression of MgO. Philos. Mag. 91, 1150–1162 (2011).
Guo, J. J. et al. Sample size induced brittle-to-ductile transition of single-crystal aluminum nitride. Acta Mater. 88, 252–259 (2015).
Csanádi, T. et al. Plasticity in ZrB2 micropillars induced by anomalous slip activation. J. Eur. Ceram. Soc. 36, 389–394 (2016).
Higashi, M., Momono, S., Kishida, K., Okamoto, N. L. & Inui, H. Anisotropic plastic deformation of single crystals of the MAX phase compound Ti3SiC2 investigated by micropillar compression. Acta Mater. 161, 161–170 (2018).
Uchic, M. D., Dimiduk, D. M., Florando, J. N. & Nix, W. D. Sample dimensions influence strength and crystal plasticity. Science 305, 986–989 (2004).
Dimiduk, D. M., Woodward, C., LeSar, R. & Uchic, M. D. Scale-free intermittent flow in crystal plasticity. Science 312, 1188–1190 (2006).
Uchic, M. D., Shade, P. A. & Dimiduk, D. M. Plasticity of micrometer-scale single crystals in compression. Annu. Rev. Mater. Res. 39, 161–186 (2009).
Schneider, A. S. et al. Correlation between critical temperature and strength of small-scale bcc pillars. Phys. Rev. Lett. 103, 10550114 (2009).
Devincre, B. & Kubin, L. Scale transitions in crystal plasticity by dislocation dynamics simulations. C.R Phys. 11, 274–284 (2010).
Deng, X., Bitler, J., Chawla, K. K. & Patterson, B. R. Toughness measurement of cemented carbides with chevron-notched three-point bend test. Adv. Eng. Mater. 12, 948–952 (2010).
Lee, S. W. & Nix, W. D. Size dependence of the yield strength of fcc and bcc metallic micropillars with diameters of a few micrometers. Philos. Mag. 92, 1238–1260 (2012).
Dimiduk, D. M., Uchic, M. D. & Parthasarathy, T. A. Size-affected single-slip behavior of pure nickel microcrystals. Acta Mater. 53, 4065–4077 (2005).
Fujita, S., Maeda, K. & Hyodo, S. Dislocation glide motion in 6H SiC single crystals subjected to high-temperature deformation. Philos. Mag. A 55, 203–215 (1987).
Vitek, V. Intrinsic stacking faults in body-centered cubic crystals. Philos. Mag. 18, 773–786 (1968).
Vitek, V. & Paidar, V. Chapter 87, Non-planar dislocation cores: a ubiquitous phenomenon affecting mechanical properties of crystalline materials in Dislocations in Solids, Vol. 14 (ed. Hirth, J.P.) 439–514 (Elsevier, 2008).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B. 54, 11169–11186 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
Monkhorst, H. J. & Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B. 13, 5188–5192 (1976).
Acknowledgements
This work was supported by JSPS KAKENHI Grant numbers JP18H01735, JP18H05478, JP18H05450, JP18H05451, JP19H00824, JP19K22053, JP20K21084 and the Elements Strategy Initiative for Structural Materials (ESISM) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) (Grant number JPMXP0112101000), and in part by JST-ALCA Grant Number JPMJAL1004.
Author information
Authors and Affiliations
Contributions
H.I. and K.K. conceived the study. T.F. and T.M. prepared the specimens and performed the micropillar compression and SEM experiments. K.K. performed the first-principles DFT calculations and prepared all figures. K.K. and H.I. wrote the main manuscript. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kishida, K., Fukuyama, T., Maruyama, T. et al. Room temperature deformation of single crystals of Ti5Si3 with the hexagonal D88 structure investigated by micropillar compression tests. Sci Rep 10, 17983 (2020). https://doi.org/10.1038/s41598-020-75007-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-75007-7
This article is cited by
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.