Structural evolution of nanoscale metallic glasses during high-pressure torsion: A molecular dynamics analysis.

Structural evolution in nanoscale Cu50Zr50 metallic glasses during high-pressure torsion is investigated using molecular dynamics simulations. Results show that the strong cooperation of shear transformations can be realized by high-pressure torsion in nanoscale Cu50Zr50 metallic glasses at room temperature. It is further shown that high-pressure torsion could prompt atoms to possess lower five-fold symmetries and higher potential energies, making them more likely to participate in shear transformations. Meanwhile, a higher torsion period leads to a greater degree of forced cooperative flow. And the pronounced forced cooperative flow at room temperature under high-pressure torsion permits the study of the shear transformation, its activation and characteristics, and its relationship to the deformations behaviors. This research not only provides an important platform for probing the atomic-level understanding of the fundamental mechanisms of high-pressure torsion in metallic glasses, but also leads to higher stresses and homogeneous flow near lower temperatures which is impossible previously.

Metallic glasses (MGs) show unique mechanical properties such as high strength and elastic limit compared to their crystalline counterparts1234. However, their applications are limited by failing under one dominant shear load when also under uniaxial load at room temperature56. Rodney and Schuh have addressed that the current challenges in this field are both technological, improving the mechanical reliability of MGs, and fundamental, developing a consistent theory of plastic deformation similar to the dislocation theory7. Inspiringly, MGs show enhanced strength, ductility, and fracture toughness when shrunk to the nanoscale89. In fact, the intrinsic mechanical length scale in MGs is associated with clusters of about 100 atoms (i.e., a few atomic diameters in size) that serve as shear transformation zones (STZs) upon deformation1011. It has been shown that an initial STZ may trigger a cascade of rearrangements during deformation12. At yield stress, intrinsically, the STZs coalesce and assemble into large and nano-scaled thin planar bands, generally called shear bands13. For example, Sha et al. concluded that when the size of the STZ aggregates accumulating at the MG surface reaches a critical value comparable to the shear band width, shear band initiation can take place14. Shear band formation and propagation is a highly localized and nano-scale deformation mechanism in MGs strained at room temperature151617. Our previous work that considers the “liquid-like” character of shear bands provided important details that extend earlier studies on shear bands18. The deformation behavior limits the structural reliability of MGs. When test temperature approached the glass transition temperature, there is an inhomogeneous-to-homogeneous flow transition, which can be described by the STZ model19202122. Fan et al. found that subnano-scale rearrangements of a small number of atoms could trigger thermally activated deformation, and activation barrier, i.e. activation energy, is independent of temperature23. It is worth noting that shear bands even occur at temperature close to glass transition temperature because of high strain-rate24. So, large and homogenous plastic strain in MGs at lower temperatures and higher stresses is hard to realise2425. To alleviate this problem, grasping how STZs coalesce into shear bands is very necessary, thereby controlling their assembling manner.

High-pressure torsion (HPT) is a popular method applying a complex stress state to achieve large plastic strain, where the sample is subjected to a compressive applied pressure and a torque simultaneously, and is used to fabricate ultrafine grained crystalline materials2627. For example, Révész et al.28 found that there are differences in the sign of the atomic distortion obtained in direction parallel to the sample surface and along the cross section, indicating strong structural anisotropy in MGs by HPT at room temperature. This indicates that HPT can be applied to activate microscopic deformation carriers and achieve large strain therein in MGs. Intensive work has also shown that HPT processing could promote structural rejuvenation in MGs and lead to homogenous deformation293031. For example, using high-energy X-ray diffractio n, Dmowski et al.29 addressed the changes in the pair distribution functions by the HPT technique, suggesting structural rejuvenation by heavy deformation. Meng et al. found that HPT could make the deformation mode change from heter ogeneous, localized deformation to homogeneous deformation in Zr50Cu40Al10 MGs30. Adachi et al. pointed that HPT could induce structural defects, which overcomes the compositional limitations and improves the mechanical properties of MGs31. These pioneering studies, all aimed at understanding the HPT, had a significant impact. Nevertheless, how these shear transformations coalesce and assemble into shear bands, the correlation between the characteristics of the shear transformations and the atomic structure, and the energy distribution under HPT remain elusive, thereby making it difficult to study the shear transformation, its activation and c haracteristics, and its relationship to the deformations behaviors under HPT.

In this study, we analysed the evolution of shear transformations aggregation under HPT at room temperature in nanoscale Cu50Zr50 MGs through molecular dynamics (MD) simulations. The MD simulations are limited to high cooling rates and small sizes but provide the mechanistic observations needed to relate the shear transformations and HPT deformation in nanoscale MGs32. It is shown that shear bands can be suppressed and constrained by HPT in nanoscale Cu50Zr50 MGs, resulting in forced cooperative flow. It is worth noting that multiple and extended shear band formation in MGs during HPT has already been shown in the experimental work33. Contrary to this, we clearly demonstrated that the nanoscale MGs during HPT can lead to homogeneous flow, which may be cooling rate and size differences between the model and experiments (see Figs S1 and S2 in Supporting information). And abundant atomic-scale shear transformations cooperate strongly with each other in the nanoscale MGs under HPT. The atoms possess lower five-fold symmetries and higher potential energies during HPT, making them more likely to participate in shear transformations. Factors leading to these novel deformation mechanisms in nanoscale MGs under HPT are discussed, and the effect of torsion speed on the deformation behaviour of nanoscale MGs was also considered.

Results and Discussion

Nanoscale Cu50Zr50 MGs and Cu50Zr50 crystals under HPT

Figure 1 shows a sequence of images that demonstrate deformation in nanoscale Cu50Zr50 MGs under HPT. The deformation process was monitored by using the so-called local atomic shear strain η, referenced to the relaxed glass prior to loading34. The presence of regions with relatively large local atomic shear strain indicates a high density of shear transformations35, which represents the collective and inelastic shearing of atoms in response to an applied shear stress3637. With increasing torsion angle (TA), the number of atoms participating in the deformation gradually increased. Shear translations were distributed across most of the outer surface of the MGs at the initial stages of deformation, such as TA = 45°. With increasing TA the whole of the outside surface turned red (see Movie S1 in Supplementary materials). Meanwhile, in accordance with the radial dependent shear deformation, the MG exhibited strain gradient, i.e., when the TAs were identical, the local atomic shear strain η of atoms located at a larger radius were greater. With the increase of TA, local atomic shear strain η of atoms located near the center also increased (see Movie S2 in Supplementary materials). The above meant that the atoms were gradually involved in shear transformations under HPT deformation. Shear transformations could pervade themselves into the whole effective area in the MG, which implied that the HPT in MGs may realise large continuous deformation. Shear transformations may have strong dependence on the preparation history of the sample. Our previous work proved that higher initial quenching temperatures or cooling rates decrease short range order and clusters in MGs38. Smaller clusters prevail in liquids at higher initial quenching temperature, and are passed into the glass by fast quenching. Under such conditions, shear transformations will be easier to cooperate strongly with each other, and propagate throughout the whole MG. Besides the preparation history of the sample, the size of MGs makes the deformation mode change from localized deformation to homogeneous deformation when the relevant size is under a critical value, which is interpreted in terms of competition between heterogeneous shear band and homogeneous cooperative flow39.

Figure 1

A sequence of surface images (up) and cross-sectional images (down) that demonstrate deformation in nanoscale Cu50Zr50 MGs under HPT.

For comparison purposes, Cu50Zr50 crystal loading was performed under the same conditions as used on the Cu50Zr50 MG in the present work, as shown in Fig. 2. Through-bands began to form on the outer surface of the models at the initial deformation. With increasing TA, bands gradually increased in number, but there was no evidence that almost all atoms were involved in the deformation as was seen in MGs (see Movie S3 in Supplementary materials). In other words, a few parts of the crystal structure were involved in the deformation. This may be due to the crystalline structure resulted in non-uniform properties at the length-scale of their microstructure.

Figure 2

A sequence of surface images (up) and cross-sectional images (down) that demonstrate deformation in nanoscale Cu50Zr50 crystals under HPT.

The above results showed that the evolution of localised shear strain in nanoscale Cu50Zr50 MG under HPT was distinctly different from that in nanoscale Cu50Zr50 crystals. At an atomic scale, accommodation of shear strain in nanoscale MGs under an applied stress was believed to occur by local rearrangement of atoms around free volume regions, unlike motion in crystalline alloys. With increasing TA, several embryonic shear transformations regions in the nanoscale MG merged and kept growing. However, the size of the shear transformations aggregates accumulating in the nanoscale MG during HPT does not reach a critical value comparable to the shear band width, so no apparent shear band formation was observed, which suggested that HPT may facilitate forced cooperative flow by activation of numerous diffuse shear transformations for nanoscale MGs. Nanoscale MGs effectively exhibited forced cooperative flow consistent with typical superplastic behaviour under HPT. In this case, not only did a large number of distributed shear transformations nucleate and remain active as deformation progressed, but new shear transformations appeared, creating an increasingly dense shear transformation network. The driving force for shear band propagation was the stored elastic energy40. For nanoscale MGs under HPT, the stored energy was dissipated into the activation and formation of a large number of homogeneously distributed shear transformations rather than into a localised single shear band as seen under uniaxial compression. The homogeneous energy release was not sufficient to drive shear band propagation. Consequently, nanoscale MG under HPT deformed homogeneously, which didn’t exhibit localised deformation by one major shear band propagation event.

HPT and uniaxial compression of nanoscale Cu50Zr50 MGs

We were particularly interested in shear band formation mechanisms. In order to understand the deformation behaviour better, in the following section, the change of displacement at atomic level is shown in detail. As shown in Fig. 3(a), atomic displacement vectors show that atoms rotated in an anticlockwise manner on the whole under HPT: the displacement of atoms located in the outer ring was large, while the displacement of atoms located near the center was small. On the whole, all atoms moved in a regular fashion. Figure 3(b) shows that atoms of Cu50Zr50 MG in the shear band moved desultorily compared with those in the matrix under uniaxial compression. These atoms interacted with each other and were irregular, which resulted in instability during the deformation process. The movement of all atoms under HPT is relatively regular, and this suggests that the HPT deformation of MGs was stable. Under applied HPT, shear transformations could multiply and pervade the whole effective area in the MG, which renders the shear transformations unable to develop a mature shear band, resulting in more stable shear deformation.

Figure 3

Atomic displacement vectors and local shear strain under (a) HPT and (b) uniaxial compression in nanoscale Cu50Zr50 MGs. The yellow arrows represent the atomic displacement vectors, and atoms are coloured according to their local shear strain.

The process of shear band development comprises shear band initiation and shear band propagation. The size of the shear transformations aggregates accumulating in the MG under uniaxial compression reach the critical value comparable to the shear band width, so a mature shear band was observed. However, during HPT deformation, shear band initiation could not occur as the size of shear transformation aggregates accumulating at the MG surface didn’t reach a critical value comparable to the shear band width; on the other hand, shear band propagation was also blocked because of the atomic shear transformations. This may have important implications for improving the deformation behaviour of MGs by adjusting shear band formation.

Quantifying the shear transformations

As the evolution of shear transformations is characterised by atomic shear strain during HPT, the fraction of atoms with relatively large atomic shear strain (i.e., STZs) during HPT is shown in Fig. 4. Atoms with a higher than 20% local shear strain were considered to achieve large shear deformation41. Thus, HPT deformation can be quantified as the ratio of the number of atoms in the large shear deformation region relative to the total number of atoms in the sample. Figure 4(a) shows the fraction of atoms with relatively large atomic shear strain, η > 0.2, as a function of TA. It was found that, when TA = 360°, their fraction in MG was 76.2%, while that in crystal was 43.2%. During deformation, the fraction of atoms undergoing large shear deformation in MG reached 76.2%, implying that a substantial fraction of the material experienced structural changes that facilitated shear transformations. This quantitatively illustrated that the MGs experienced shear transformations, while only to a partial extent, atoms of crystals were involved in HPT deformation. At the same time, the MG reached a plateau at TA = 150°, while the crystals reached their plateau at TA = 75°. This suggested that the MG can resist the damage induced by HPT more effectively in the longer-term. MG, under HPT, underwent a substantial fraction of structural changes that facilitated shear transformations which, in turn, gave rise to the observed superplastic-like behaviour.

Figure 4

(a) The fraction of atoms with relatively large atomic shear strain in MGs and crystals as a function of torsion angle; (b) The fraction of atoms with relatively large atomic shear strain in MGs as a function of torsion angle for various torsion periods.

To understand, or predict, the HPT deformation of MGs better, it is important to make a qualitative investigation of the effects of torsion period, e.g., torsion speed. The results demonstrated the significant influence of torsion period on the deformation behaviour of MGs. Figure 4(b) shows the fraction of atoms with relatively large atomic shear strain as a function of TA for various torsion periods. It can be seen that the accumulation rate of shear transformations was sensitive to torsion period. With increasing torsion period, the shear transformations accumulated in abundance because of the adequate duration allowing this to come about. The effect of torsion period on the HPT deformation of MGs suggested that a higher torsion period led to a greater extent of shear transformations. The results indicated that the activation of a high density of shear transformations can readily lead to the percolation of shear transformations or macroscopic plastic deformation in MGs.

The “liquid-like” character during HPT

Figure 5(a) shows the distributions of the top 14 Voronoi polyhedra (VP) of the Cu50Zr50 MG at TA = 0° and TA  = 360°, and Cu50Zr50 liquid, where the element types were not distinguished. In the Voronoi tessellation method42, the simulation cell is divided into VP around each atom. The polyhedra can be characterised by the Voronoi index [n3, n4, n5, n6], where n denotes the number of i-edged faces of the VP. An i-edged face reflects the local symmetry of the central atom with some nearest-neighbour atoms in a certain direction. The distributions of Cu50Zr50 MG (TA = 0°) were mainly located in the side of pentagon dominant VP, such as <0,2,8,1> and <0,2,8,4>. While the distributions of Cu50Zr50 MG (TA = 360°) showed a close resemblance to that of the corresponding Cu50Zr50 liquid. They were located in the side of VP, where pentagon accounted for a minor proportion of the total, and <0,2,8,1> VP were regarded as the watershed thereof. This suggested that HPT destroyed the short-range order structure, in which fivefold symmetry was dominant.

Figure 5

(a) The distributions of the top 14 Voronoi polyhedra of the Cu50Zr50 MG at TA (torsion angle) = 0° and 360°, and Cu50Zr50 liquid; (b) Comparison of bonded pairs in the three aforementioned models.

To quantify free volume change during HPT deformation, the content and corresponding volumes of the top 14 Voronoi polyhedra (VP) were calculated. To be more comparable with this approach, a parameter43 was proposed such that:

where V denotes the effective free volume, N denotes the total number of types of VP, V is the volume of VP i, and P is the content of VP i. The V of Cu50Zr50 MG at T = 360° (16.76) was greater than that of Cu50Zr50 MG at T = 0° (15.82), and was close to the V of Cu50Zr50 liquid (17.91). Figure 5(b) shows the comparison of bonded pairs. A bonded pairs index with four integers ijkl was raised by Honeycutt and Andersen (H-A)44. The first number i in the H-A index is used to identify the bonding of two given atoms (i = 1 for bonded pairs and i = 2 for non-bonding atoms); j is the number of nearest neighbors shared jointly by the two atoms; and k is the number of bonds among the shared neighbors. The fourth digit l is needed in case the first three numbers are the same but the bond geometries are different. To analyse further that information contained in the bonded pairs, Pan45 divided the bonded pairs into two categories: when the shared neighbours connected to form a ring, i.e. when k = j, they were designated “saturated bonded pairs”; otherwise, they were designated “unsaturated bonded pairs”. The degree of saturation of the system, I, is defined as:

where j is index j of bonded pairs z, and k is index k of bonded pairs z. Here c is the proportion of bonded pairs z. The greater I represents the more defects. And when MGs have no defect (ideal MGs), the I is 0. In contrast with Cu50Zr50 MG at TA = 0°, 1551 saturated bonded pairs of Cu50Zr50 MG, at TA = 360°, were destroyed. Through statistical calculations, the unsaturated index of the Cu50Zr50 MG at TA = 0° was found to have been 0.344, while that of the Cu50Zr50 MG at TA = 360° was 0.504. Note that the unsaturated index of Cu50Zr50 MG at TA = 360° was more similar to that (0.839) found in the Cu50Zr50 liquid. It also verified the fact that the HPT destroyed some five-fold symmetry in MG, rendering it easier to experience shear transformations. Five-fold symmetry is associated with dense packing in amorphous materials46. So the destruction of five-fold symmetry reduced dense packing, and increased the free space within which atoms can move, a result which was consistent with the expansion of the effective free space V in Cu50Zr50 MG at T = 360°. These suggested that HPT made the atoms in the MGs more homogeneous and more akin the liquid state, providing MGs with effective free volume so that they could experience shear transformations.

The potential energies during HPT

The observed HPT-induced emergence of shear transformations in MGs can be rationalised in terms of energetics. To reveal physical mechanisms of deformation as a function of energetics, the potential energy at TA = 0° and 360° was calculated, as shown in Fig. 6. It shows that the potential energy at TA =  360° was higher than that at TA = 0° except for those located within the clamps. HPT made atoms have higher potential energies, which increased their probability of participating in plastic deformation by way of a shear transformation because the energy cost of moving such atoms was lowered, akin to the free volume driven mechanism for homogeneous flow. This implies that the structure of the deformed volume becomes effectively rejuvenated by the HPT29. It was reasonable to conclude that the HPT stress state, with its greater free volume and higher potential energy per atom, may suppress catastrophic shear band formation and explain the simulated observations presented here. The shear transformation involved a localised cluster of atoms that underwent intense distortion from an initial to final equilibrium position through an intermediate activated state of high-energy and large-volume47. In other words, if enough energy were applied to activate a sufficiently high density of shear transformations, then homogeneous plastic deformation can occur in any MGs at room temperature depending on the activation, density, and percolation mode of the shear transformations.

Figure 6

Change in potential energies of atoms of the Cu50Zr50 MG at TA = 0° and 360°.

Conclusion

In summary, we have performed a series of MD simulations of HPT in nanoscale MGs. It is found that the nanoscale MGs experienced forced cooperative flow. In this process, shear transformations under applied HPT stress could multiply and pervade the whole effective area in the MG, meaning that there is no obviously stress concentration. Therefore shear transformations cooperate and homogenously flow could be realized. Nanoscale MGs, compared to their crystalline counterparts, conferred significant performance advantages under HPT. We demonstrate the effect of torsion speed on the deformation behaviour of nanoscale MGs, with results suggesting that a higher torsion period led to a greater degree of forced cooperative flow. Our results show that HPT could prompt atoms to possess lower five-fold symmetries, higher potential energies and larger volume, making them more likely to participate in shear transformations.

Methods

The simulations were performed using the embedded atom method (EAM) potential supplied in large-scale atomic/molecular massively parallel simulator (LAMMPS)4849. The dimensions of the Cu50Zr50 MGs model system used in the calculations were 15.5 × 15.5 × 30 nm in the x-, y-, and z-directions, respectively, containing approximately 420,000 atoms. It was melted and equilibrated with periodic boundary conditions for 1 ns at 2000 K, which is above the melting point of these alloys. The model was then cooled at an MD cooling rate of 0.5 K ps−1 to a glass state (300 K). As a comparison, the B2 Cu50Zr50 crystal model, containing approximately 420,000 atoms, was also built and relaxed with periodic boundary conditions at 300 K for 5 ns to ensure the energy minimization of the model. The NPT ensemble (constant number, constant pressure, and constant temperature) was used for these processes, with the temperature controlled by a Nose-Hoover thermostat and the pressure controlled at zero using a Nose-Hoover barostat.

The cylindrical-shaped nanoscale MG and crystal models are prepared by cutting, 15 nm in diameter and 30 nm in the z-axis. With 1 GPa applied pressure in the z-direction, the models were relaxed for 3 ns at 300 K. Thereafter, these models were simulated by rotating the rigid atoms at one end (3 nm in thickness) along the z-axis, while keeping the rigid atoms at the other end (3 nm in thickness) unchanged, being akin to the clamps. Free surfaces were imposed in the x- and y- directions during the torsion. Torsion speeds were 1/60 revolution/ps, 1/600 revolution/ps, 1/6000 revolution/ps, and 1/60000 revolution/ps, respectively.

Additional Information

How to cite this article: Feng, S. D. et al. Structural evolution of nanoscale metallic glasses during high-pressure torsion: A molecular dynamics analysis. Sci. Rep. 6, 36627; doi: 10.1038/srep36627 (2016).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.