Investigations about the Atomic Structure and Mechanical Behavior of Metallic Glasses after Melt Hydrogenation

Binbin Wang, Yanqing Su, Liangshun Luo, Liang Wang, Ruirun Chen, Jingjie Guo and Hengzhi Fu

(National Key Laboratory for Precision Hot Processing of Metals, School of Materials Science and Engineering,Harbin Institute of Technology, Harbin 150001, China)

Abstract: “Hydrogen in metallic glasses” has become a popular topic for material scientists, yet few studies focus on the atomic-scale details. Herein, by utilizing molecular dynamic simulations, the changes on the atomic structure of Cu50Zr50 metallic glasses after melt hydrogenation were systematically analyzed, with the aim of understanding the differences of mechanical behavior between these amorphous alloys. The simulated analyses indicate that the hydrogenated samples become more compact than the H-free one, but the fraction of the dominant coordination polyhedra with higher degree of local fivefold symmetry significantly decreases accompanied by the addition of H atoms. Accordingly, melt hydrogenation can induce much more local “soft spots” in metallic glasses to alleviate the degree of strain localization during deformation, i.e., it has a positive influence on the deformability of glassy alloys in agreement with experimental results.

Keywords: amorphous alloys; molecular dynamic simulations; atomic structure; deformation; melt hydrogenation

1 Introduction

As a new class of potential structural and functional materials, metallic glasses (MGs) present a number of superior physical and chemical properties[1], which have stimulated widespread research enthusiasm. Nevertheless, under the quasistatic loading, these amorphous materials always break in an apparently brittle manner at ambient temperature, severely compromising their engineering applications[2]. To overcome this shortcoming, lots of efforts have been made to improve their plasticity[3-5].

Based on the experimental and computational results, the mechanical behavior of MGs should be dependent on their chemical composition. Park et al.[6]surveyed the mechanical properties of CuxZr100-xbinary MGs; with the increment of Cu contents, the yield strength continuously increases, while the plastic strain decreases rapidly. In addition, Kumar et al.[7]reported a few percent of Al innovated into Cu50Zr50amorphous alloys leads to the same mechanical behavior change as the increase of Cu contents. Since many Zr-Cu-based systems manifest excellent performance, Cheng et al.[8-9]established the correlation between their internal structure and mechanical behavior employing molecular dynamics (MD) simulations. Following these results, one can penetrate how and why the alloying elements influence the mechanical properties of MGs.

Nowadays, MGs employed as hydrogen-separable membranes have received a great attention due to their comparable H permeability, lower cost, and higher durability in hydrogen atmosphere than Pd and Pd-based crystalline membranes[10-12]. Meanwhile, during the H charging, it was found that H additions can induce the emergence and propagation of shear bands from the tip of notch in MGs under the continuous load[13]. Subsequently, Zhao et al.[14-15]studied the impact of charging H atoms on the yield shear stress and the volume of shear transformation zones (STZs) of Zr-bearing amorphous alloys based on the first pop-in load from the tests of nanoindentation. The mechanical response of hydrogenated alloys is significantly dependent on Zr contents, i.e., the addition of hydrogen can induce softening in low-Zr samples by increasing the volume of STZs while reversing in high-Zr alloys.

In our recent work, hydrogen was doped into Zr-Cu-based MGs by alloying with a gaseous mixture of hydrogen/argon, i.e., the so-called melt hydrogenation, which is also a hydrogen microalloying technology[16]. Amazingly, minor addition of H atoms could notably optimize the deformability of MGs at ambient temperature[17]. However, the internal mechanism of this H-dependent mechanical behavior remain unclear. Therefore, in this context, we employed molecular dynamics (MD) simulations to elucidate the evolution of structure and mechanical behavior of MGs after hydrogen microalloying, with the purpose of accounting for the structural origin of H-induced plasticity enhancement.

2 Methods

MD simulations were performed to investigate the characteristic of atomic arrangement and its evolution during the deformation of MGs employing the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) codes[18]. All the atomistic simulations employed the recently developed 2NN MEAM (second nearest-neighbor modified embedded atom method) potential as the interatomic interaction of Cu-Zr-H system, which could reproduce certain fundamental physical properties (the lattice parameter and bulk modulus of intermetallic compounds) of alloys in consistence with the results from experimental measurements or first-principle calculations[19-20]. It can also generate a relatively reliable amorphous structure for the amorphous system.

The original Cu50Zr50sample with 10 000 atoms was firstly melt and equilibrated at 2 000 K for 2 ns (time step 1 fs) , and then quenched into glassy states (300 K) at a cooling rate of 1 K/ps. The external pressure of each three-dimensional (3D) direction was held at zero by using an isothermal-isobaric ensemble (NPT) during the whole process. The cooled sample was sufficiently relaxed at 300 K for 1 ns. Periodic boundary conditions (PBCs) were applied for all directions in order to eliminate the surface effect. A typical amorphous configuration is shown in Fig. 1(a).

Fig. 1 (a) 3D-configuration for Cu50Zr50 amorphous sample with 10 000 atoms achieved by MD simulation; (b) Schematic for the uniaxial tensile loading. The Cu and Zr atoms are respectively colored by orange and grey.

During the experimental process of melt hydrogenation, master alloys are melt under the gaseous mixture of H2/Ar to make the H atoms diffuse into the liquid alloy, i.e., hydrogen atoms are introduced into matrix (Cu50Zr50) during the melt stage. Thus, in our simulations, hydrogen atoms were randomly doped into the systems with molten state (2 000 K) before quenching. The amount of inserted hydrogen atoms defined as the ratio of the number of inserted H atoms (H) to that of the Cu and Zr atoms (M) was respectively 0%, 5%, and 10%, since only a small quantity of hydrogen atoms could be absorbed by matrix.

Fig. 1(b) gives the loading schemes adopted for deformation. Uniaxial tensile loading with a strain rate of 108s-1was imposed on each final sample (300 K). To ensure the stress-free state, the direction perpendicular to the loading direction should undergo a Parrinello-Rahman barostat. PBCs were applied for all directions.

3 Results and Discussion

3.1 Atomic Structure from Simulations

Fig. 2 presents the distribution for the coordination numbers (CNs) of MGs with different H contents. In Fig. 2(a), the distribution about the CN of total atoms indicates that the fractions of small-sized clusters (CN=11-12) and medium-sized clusters (CN=15) decrease while those of medium-sized clusters (CN= 13) and large-sized (CN=17-18) ones increase after the addition of H atoms. As only the Cu-centered clusters given in Fig. 2(b) are considered, the CN distribution presents a similar H-induced trend, i.e., the proportions of small-sized clusters (CN=10-12) decline, in turn the large-sized clusters (CN=13-16) increase. Considering the relationship of atomic radiiramong the three elements (rH

Fig. 2 Distributions of major CNs for (a) total, (b) Cu, (c) Zr, and (d) H atoms in samples with different H contents. Only the fraction over 1.0% is listed.

The fraction of major Voronoi polyhedra in each sample can be seen in Fig. 3(a); it gradually reduces with the addition of H atoms. In addition, some polyhedra centered by Cu, Zr, or H atoms are individually depicted in Figs.3(b)-(d). Ma et al.[21]found that there is the highest possible packing efficiency for the atomic arrangement around Cu atoms with a given partial CNs, whereas a less efficient packing can be found around Zr atoms through using XRD and neutron diffraction. These findings are also consistent with the results from MD simulations[22]. It can be deduced that the structure of ZrCu-based MGs can be better analyzed from the perspective of Cu-centered coordination polyhedra. In Fig.3(b), three characteristic clusters are respectively indexed as <0, 0, 12, 0>, <0, 2, 8, 2>, and <0, 2, 8, 1>, in line with the previous results[8-9]. Their variation after hydrogen microalloying presents a decreased behavior. In contrast, some high-coordinated clusters are found to increase, such as <0, 1, 10, 2> polyhedra. As displayed in Fig. 3(d), the dominant polyhedra centered by H atoms are respectively indexed as <0, 4, 4, 2>, <1, 2, 5, 2>, and <0, 4, 4, 3>, which can be regarded as the trigonal prisms attributed to a type of distorted bi-capped square Archimedean antiprism (BASP) with two half octahedral (AA) and defect-AA. In addition, the <0, 3, 6, 0> polyhedra, corresponding to the tricapped trigonal prism (TTP) consisting of three half octahedral, is also noteworthy. Hence, the local atomic packing environment around most of H atoms should be featured by the trigonal prisms.

Fig. 3 Histogram displaying the fraction of major Voronoi polyhedra in samples with different H contents (a) Total atoms; (b) Cu atoms; (c) Zr atoms; (d) H atoms. Only the fraction over 1.0% is plotted.

3.2 Stress-Strain Curves and Local Atomic Strain

As is known, the deformability of MGs must be correlated with the initiation and proliferation of shear bands[2]. Therefore, the various deformabilities owned by the amorphous samples can be properly accounted by probing the flow characteristics inside a shear band and its corresponding structural evolution. Fig. 4 illustrates the simulated stress-strain curves of the deformed MGs with different H contents. Considering the typical thickness of shear band is 10-20 nm that is much larger than the size of each specimen, its mechanical response can be treated as the characteristics inside the shear band. For each specimen, stress increases step by step as deformation progresses in the elastic stage. As it exceeds the critical yield strengthτy, a precipitous drop can be seen; finally, beyond ～2% loading strain, the stress tends to be saturated.

Meanwhile, it can be perceived from Fig. 4 thatτydecreases with the addition of H atoms. However, the steady-state stressτs, corresponding to the first minimum of stress after global yielding[23], presents a slight increase in the H-alloyed samples. Correspondingly, the value of stress drop (τ=τy-τs) declines after hydrogen microalloying, as shown in the right panel of Fig. 4. Following the recent simulated studies, the softening degree of deformed MGs can be evaluated byτ, which corresponds to an indicator for plastic response stemming from the structural change of MGs; this may be concomitant by strain localization in shear bands[24-25]. Thus, the magnitude ofτmaps the extent of strain localization in flow region[9], and a lowerτof the H-alloyed samples means weaker degree of localization than in the H-free one. Therefore, minor additions of hydrogen can improve the plasticity of MGs through inducing local softening behavior, which is consistent with the previous results from nanoindentation test demonstrating that a decreased elastic modulus and hardness can be found in the nano-scale regions of MGs with the similar H contents[26].

Fig. 4 Simulated stress-strain curves (left panel) and variation in the values of stress drop τ (right panel) under uniaxial tensile loading imposed on the amorphous samples with different H contents. The steady-state stress τs for each sample is highlighted in the inset

As reflected in the central materials science paradigm that “structure determines properties”, the different degrees of softening (τ) for MGs in Fig. 4 must relate to their various atomic structures. The local irreversible rearrangement in each simulated configuration was assessed by the nonaffine displacement of the central atomirelative to its neighbor atomsj, which can be defined as[27]

(1)

whereNiis the total nearest neighbors of theithatoms

Fig. 5 Projected views of atomic arrangement in response to tensile loading captured from (a) H/M=0 and (b) H/M=10% samples at the strain of 8% (before global yielding).Each atom is colored by its own value of nonaffine displacement.

3.3 Structural Evolution during Deformation

Peng et al.[28]have demonstrated that a structural indicator correlated with the degree of local fivefold symmetry (LEFS) can be correlated with the deformation of MGs. The degree of LEFSd5is defined as the ration of the number of pentagonsn5to that of the nearest-neighbor atomsniin Voronoi polyhedra:

(2)

In order to examine the role of the polyhedra with different degrees of LEFS during deformation, the evolution for the fraction of dominant polyhedra associated with deformation in H-free Cu50Zr50is firstly monitored, as shown in Fig. 6. The color bar maps the various levels of LEFS.

Fig. 6 Variations of the fraction of dominant polyhedra during tensile deformation of Cu50Zr50 (H/M=0). Note that the color bar maps the various levels of LEFS

As the stress ascends and reaches the yielding point, the polyhedra with higher degree of LEFS, especially <0, 0, 12, 0> (d5=1.0), begins to breakdown; conversely, the clusters with lower degree of LEFS increase, such as <0, 3, 6, 4>. Accordingly, the polyhedra with medium degree of LEFS keeps almost constant. This process seems to saturate when the stress approachesτs. Therefore, it can be suggested thatτyis the stress required to break the dense polyhedra with higher degree of LEFS, converting to the loose one with lower degree, such that flow can be initiated. That is, the polyhedra with a higher degree of LEFS is “harder”, while the polyhedra with lower degree of LEFS is “softer” in the initial amorphous system. Following the results from Fig.3, the fraction of dominant polyhedra with higher degree of LEFS significantly decreases with the addition of H atoms; it should be noted that diverse low fractions of polyhedra with lower degree of LEFS cannot be seen in Fig.3. In a word, hydrogen microalloying can develop more “soft spots” with distorted atomic packing in MGs; these fertile sites afford the nucleation of shear bands, which could mitigate the degree of strain localization and finally promotes their deformability.

4 Conclusions

Doping hydrogen into Cu50Zr50amorphous matrix can induce structural change by increasing the average CN, i.e., the hydrogenated sample is more compact. Moreover, the fraction of dominant polyhedra with higher degree of LEFS significantly decreases with the addition of H atoms, while diverse low fraction of polyhedra with lower degree of LEFS appears in hydrogenated alloys. This structural evolution can induce local softening, weaken the degree of strain localization, and further improve the plasticity of MGs.