'Crystal Genes' in Metallic Liquids and Glasses

Comprehensive characterization of the structure of Zr-based metallic glasses

Structure of metallic glasses fascinates as the generic amorphous structural template for ubiquitous systems. Its specification necessitates determination of the complete hierarchical structure, starting from short-range-order (SRO) → medium-range-order (MRO) → bulk structure and free volume (FV) distribution. This link has largely remained elusive since previous investigations adopted one-technique-at-a-time approach, focusing on limited aspects of any one domain. Reconstruction of structure from experimental data inversion is non-unique for many of these techniques. As a result, complete and precise structural understanding of glass has not emerged yet. In this work, we demonstrate the first experimental pathway for reconstruction of the integrated structure, forZr 67 Ni 33 andZr 52 Ti 6 Al 10 Cu 18 Ni 14 glasses. Our strategy engages diverse (× 7) multi-scale techniques [XAFS, 3D-APT, ABED/NBED, FEM, XRD, PAS, FHREM] on the same glass. This strategy complemented mutual limitations of techniques and corroborated common parameters to generate complete, self-consistent and precise parameters. Further, MRO domain size and inter-void separation were correlated to identify the presence of FV at MRO boundaries. This enabled the first experimental reconstruction of hierarchical subset: SRO → MRO → FV → bulk structure. The first ever image of intermediate region between MRO domains emerged from this link. We clarify that determination of all subsets is not our objective; the essence and novelty of this work lies in directing the pathway towards finite solution, in the most logical and unambiguous way.

Unveiling the effect of Ni on the formation and structure of Earth's inner core

Ni is the second most abundant element in the Earth's core. Yet, its effects on the inner core's structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe's melting point at inner core pressures. The melting temperature of Ni is shown to be 700 to 800 K higher than that of Fe at 323 to 360 GPa. hcp, bcc, and liquid phase relations differ for Fe and Ni. Ni can be a bcc stabilizer for Fe at high temperatures and inner core pressures. A small amount of Ni can accelerate Fe's crystallization at core pressures. These results suggest that Ni may substantially impact the structure and formation process of the solid inner core.

Neural network potential for Zr-Rh system by machine learning

Zr-Rh metallic glass has enabled its many applications in vehicle parts, sports equipment and so on due to its outstanding performance in mechanical property, but the knowledge of the microstructure determining the superb mechanical property remains yet insufficient. Here, we develop a deep neural network potential of Zr-Rh system by using machine learning, which breaks the dilemma between the accuracy and efficiency in molecular dynamics simulations, and greatly improves the simulation scale in both space and time. The results show that the structural features obtained from the neural network method are in good agreement with the cases inab initiomolecular dynamics simulations. Furthermore, we build a large model of 5400 atoms to explore the influences of simulated size and cooling rate on the melt-quenching process of Zr₇₇Rh₂₃. Our study lays a foundation for exploring the complex structures in amorphous Zr₇₇Rh₂₃, which is of great significance for the design and practical application.

Molecular dynamics simulation of metallic Al-Ce liquids using a neural network machine learning interatomic potential

Al-rich Al-Ce alloys have the possibility of replacing heavier steel and cast irons for use in high-temperature applications. Knowledge about the structures and properties of Al-Ce alloys at the liquid state is vital for optimizing the manufacture process to produce desired alloys. However, reliable molecular dynamics simulation of Al-Ce alloy systems remains a great challenge due to the lack of accurate Al-Ce interatomic potential. Here, an artificial neural network (ANN) deep machine learning (ML) method is used to develop a reliable interatomic potential for Al-Ce alloys. Ab initio molecular dynamics simulation data on the Al-Ce liquid with a small unit cell (∼200 atoms) and on the known Al-Ce crystalline compounds are collected to train the interatomic potential using ANN-ML. The obtained ANN-ML model reproduces well the energies, forces, and atomic structure of the Al₉₀Ce₁₀ liquid and crystalline phases of Al-Ce compounds in comparison with the ab initio results. The developed ANN-ML potential is applied in molecular dynamics simulations to study the structures and properties of the metallic Al₉₀Ce₁₀ liquid, which would provide useful insight into the guiding experimental process to produce desired Al-Ce alloys.

Dynamic heterogeneity, cooperative motion, and Johari-Goldstein [Formula: see text]-relaxation in a metallic glass-forming material exhibiting a fragile-to-strong transition

We investigate the Johari-Goldstein (JG) [Formula: see text]-relaxation process in a model metallic glass-forming (GF) material ([Formula: see text]), previously studied extensively by both frequency-dependent mechanical measurements and simulation studies devoted to equilibrium properties, by molecular dynamics simulations based on validated and optimized interatomic potentials with the primary aim of better understanding the nature of this universal relaxation process from a dynamic heterogeneity (DH) perspective. The present relatively low temperature and long-time simulations reveal a direct correspondence between the JG [Formula: see text]-relaxation time [Formula: see text] and the lifetime of the mobile particle clusters [Formula: see text], defined as in previous DH studies, a relationship dual to the corresponding previously observed relationship between the [Formula: see text]-relaxation time [Formula: see text] and the lifetime of immobile particle clusters [Formula: see text]. Moreover, we find that the average diffusion coefficient D nearly coincides with [Formula: see text] of the smaller atomic species (Al) and that the 'hopping time' associated with D coincides with [Formula: see text] to within numerical uncertainty, both trends being in accord with experimental studies. This indicates that the JG [Formula: see text]-relaxation is dominated by the smaller atomic species and the observation of a direct relation between this relaxation process and rate of molecular diffusion in GF materials at low temperatures where the JG [Formula: see text]-relaxation becomes the prevalent mode of structural relaxation. As an unanticipated aspect of our study, we find that [Formula: see text] exhibits fragile-to-strong (FS) glass formation, as found in many other metallic GF liquids, but this fact does not greatly alter the geometrical nature of DH in this material and the relation of DH to dynamical properties. On the other hand, the temperature dependence of the DH and dynamical properties, such as the structural relaxation time, can be significantly altered from 'ordinary' GF liquids.

Phase Diagram and Structure Map of Binary Nanoparticle Superlattices from a Lennard-Jones Model

A first-principles prediction of the binary nanoparticle phase diagram assembled by solvent evaporation has eluded theoretical approaches. In this paper, we show that a binary system interacting through the Lennard-Jones (LJ) potential contains all experimental phases in which nanoparticles are effectively described as quasi hard spheres. We report a phase diagram consisting of 53 equilibrium phases, whose stability is quite insensitive to the microscopic details of the potentials, thus giving rise to some type of universality. Furthermore, we show that binary lattices may be understood as consisting of certain particle clusters, i.e., motifs, that provide a generalization of the four conventional Frank-Kasper polyhedral units. Our results show that metastable phases share the very same motifs as equilibrium phases. We discuss the connection with packing models, phase diagrams with repulsive potentials, and the prediction of likely experimental superlattices.

Observation of η-Al41Sm5 reveals motif-aware structural evolution in Al-Sm alloys

Using an effective genetic algorithm, we uncover the structure of a metastable Al41Sm5 phase that supplements its family sharing similar short-range orders. The phase evolves upon heating an amorphous Al-9.7 at.% Sm ribbon, produced by melt-spinning. The dynamical phase selection is discussed with respect to the structural connections between the short-range packing motifs in the amorphous precursor and those observed in the selected phases. The phase elucidated here is one of several newly discovered large-unit-cell phases found to form during devitrification from the glass in this binary system, further illustrating the power and efficiency of our approach, the important role of structural hierarchy in phase selection, and the richness of the metastable phase landscape accessible from the glassy structure.

Assessing the utility of structure in amorphous materials

This paper presents a set of general strategies for the analysis of structure in amorphous materials and a general approach to assessing the utility of any selected structural description. Two measures of structure are defined, "diversity" and "utility," and applied to two model glass forming binary atomic alloys, Cu₅₀Zr₅₀ and a Lennard-Jones A₈₀B₂₀ mixture. We show that the change in diversity associated with selecting Voronoi structures with high localization or low energy, while real, is too weak to support claims that specific structures are the prime cause of these local physical properties. In addition, a new structure-free measure of incipient crystal-like organization in mixtures is introduced, suitable for cases where the stable crystal is a compound structure.

Nucleation instability in supercooled Cu-Zr-Al glass-forming liquids

Few general models representing certain classes of real glass-forming systems play a special role in computer simulations of supercooled liquid and glasses. Recently, it was shown that one of the most widely used model glassformers-the Kob-Andersen binary mixture-crystalizes in quite lengthy molecular dynamics simulations, and moreover, it is in fact a very poor glassformer at large system sizes. Thus, our understanding of crystallization stability of model glassformers is far from complete due to the fact that relatively small system sizes and short time scales have been considered so far. Here we address this issue for two embedded atom models intensively used last years in numerical studies of Cu-Zr-(Al) bulk metallic glasses. Exploring the structural evolution of Cu_64.5Zr_35.5 and Cu₄₆Zr₄₆Al₈ alloys at continuous cooling and isothermal annealing, we observe that both systems nucleate in sufficiently lengthy simulations, although critical nucleation time for the latter is an order of magnitude higher than that for the former. We show that Cu_64.5Zr_35.5 is actually unstable to crystallization for large system sizes (N > 20 000). Both systems crystallize with the formation of tetrahedrally close packed Laves phases of different types. We argue that nucleation instability of the simulated Cu_64.5Zr_35.5 alloy is due to the fact that its composition is very close to that for the stable Cu₂Zr compound with a C15 Laves phase structure.

Atomistic Simulations to Predict Favored Glass-Formation Composition and Ion-Beam-Mixing of Nano-Multiple-Metal-Layers to Produce Ternary Amorphous Films

Based on the framework of long-range empirical formulas, the interatomic potentials were constructed for the Ni-Nb-Mo (fcc-bcc-bcc) and Ni-Zr-Mo (fcc-hcp-bcc) ternary metal systems. Applying the constructed potentials, atomistic simulations were performed to predict the energetically favored glass formation regions (GFRs) in the respective composition triangles of the systems. In addition, the amorphization driving forces (ADFs), i.e., the energy differences between the solid solutions and disordered phases, were computed and appeared to correlate with the so-called glass forming abilities. To verify the atomistic prediction, ion beam mixing with nano-multiple-metal-layers was carried out to produce ternary amorphous films. The results showed that the composition of ternary amorphous films obtained by ion beam mixing all locate inside the GFRs, supporting the predictions of atomistic simulations. Interestingly, the minimum ion dosage required for amorphization showed a negative correlation with the calculated ADF, implying that the predicted amorphization driving force could be an indicator of the glass formation ability.

Crystal genes in a marginal glass-forming system of Ni50Zr50

The marginal glass-forming ability (GFA) of a binary Ni-Zr system is an issue to be explained considering numerous bulk metallic glasses (BMGs) found in a Cu-Zr system. Using molecular dynamics, the structures and dynamics of Ni₅₀Zr₅₀ metallic liquid and glass are investigated at the atomistic level. To achieve a well-relaxed glassy sample, a sub-T_g annealing method is applied and the final sample is closer to the experiments than the models prepared by continuous cooling. With the state-of-the-art structural analysis tools such as cluster alignment and pair-wise alignment methods, two glass-forming motifs with some mixed traits of a metastable B2 crystalline phase and a crystalline Ni-centered B33 motif are found to be dominant in the undercooled liquid and glass samples. A new chemical order characterization on each short-range order (SRO) structure is accomplished based on the cluster alignment method. The significant amount of the crystalline motif and the few icosahedra in the glassy sample deteriorate the GFA.

Comparatively studying the local atomic structures of metallic glasses upon cyclic-loading by computer simulations

Through a variety of local structural analysis methods, it is revealed that the NixZr_100–xMGs exhibit a combination of the icosahedral-, fcc- and hcp-like configurations, while the icosahedra or distorted icosahedra cover a dominant fraction.