Tuning the glass-forming ability of metallic glasses through energetic frustration

The design of multifunctional bulk metallic glasses is limited by the lack of a quantitative understanding of the variables that control the glass-forming ability of alloys. Both geometric frustration (e.g., differences in atomic radii) and energetic frustration (e.g., differences in the cohesive energies of the atomic species) contribute to the glass-forming ability. We perform molecular dynamics simulations of binary Lennard-Jones mixtures with only energetic frustration. We show that there is little correlation between the heat of mixing Δ H_mix and critical cooling rate R_c, below which the system crystallizes, except that Δ H_mix<0 . By removing the effects of geometric frustration, we show strong correlations between R_c and the variables ɛ_-=(ɛ_{B B}-ɛ_{A A}) /(ɛ_{A A}+ɛ_{B B}) and ɛ_{¯A B}=2 ɛ_{A B}/(ɛ_{A A}+ɛ_{B B}) , where ɛ_{A A} and ɛ_{B B} are the cohesive energies of atoms A and B and ɛ_{A B} is the pair interaction between A and B atoms. We identify a particular composition-dependent combination of ɛ_- and ɛ_{¯A B} that collapses the data for R_c over nearly 4 orders of magnitude in cooling rate. By performing local structural analyses, we find that energetic frustration, even in the absence of geometric frustration, can induce short-range fivefold symmetric order that impedes crystallization. This result emphasizes that energetic frustration plays an important role in determining the glass-forming ability, and thus it should be taken into account in the design of new metallic glass formers.