Simulated multi-component CuZr(Al) metallic glasses akin to experiments

Abstract

We study a three-component CuZrAl metallic glass system by means of a combined Monte Carlo and Molecular Dynamics simulations scheme. This hybrid method allows us to generate equilibrated samples at temperatures below the conventional glass transition for the first time, achieving a more stable glassy regime. By using a realistic potential for the interactions of metallic species, we explore the kinetics, thermodynamics, and rheology of a CuZrAl glass, and then compare these findings with those of the ubiquitous CuZr. Remarkably, the resulting sheared glassy configurations show an abrupt stress drop corresponding to the shear band, akin to experimental observations. Our results pave the way for theoretical studies of complex metallic glasses and offer comparisons with experiments.

Metallic glasses (MGs) are an important class of materials constituted by an amorphous atomic structure, which forms when the high-temperature metallic liquid is rapidly cooling down to room temperature [1, 2, 3]. The unique disordered structure thereby formed is the key to many outstanding properties, such as unprecedented mechanical features, including high strength and high elastic limit. Hence, MGs have great potential to outperform conventional materials faor various technological and industrial applications [2, 4, 5]. Nevertheless, the major limitation that prevents their use on a large scale, is the difficulty to avoid the crystallization of the samples during the cooling process [6]. Recent advancement in synthesis techniques has come from a method different from conventional cooling: vapor deposition. This technique allows the production of MGs with enhanced stability, known as ultrastable metallic glasses (UMGs) [7, 8, 9, 10, 11]. A resultant shift in the glass transition temperature, indicative of this enhanced stability, has been observed in various UMGs such as for Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ MG in Ref. [9]. These glasses are thought to be equivalent to conventionally liquid quenched MGs that have undergone aging for thousands of years [12]. This trend was initially discovered by Ediger and co-workers but in the class of organic glasses [13]. In general, ultrastable glasses are characterized by extraordinary thermodynamic and kinetic stability and exceptional mechanical properties [13, 14, 15, 16]

The molecular mechanisms underlying the glass-forming ability (GFA), shift of glass transition temperature, and formation of enhanced stable states in MGs are far from being understood, a challenge that also characterizes the broad spectrum of glasses [17, 18]. Molecular simulations have been fundamental to investigate the properties of glasses from the microscopic point of view [19, 20]. However, simulations suffer from unrealistic high-cooling rates ($\sim[10-1000]$ K/ns), resulting in computer-generated metallic glasses with properties that differ significantly from those observed in experiments [21, 22, 23, 24, 25]. Here, we solve this gap by obtaining in silico samples of realistic CuZr-based ultrastable MGs similar to experiments, using an efficient hybrid Molecular Dynamics and Monte Carlo (MD+MC) simulation approach. Our study focus on two-component Cu${}_{50}$Zr${}_{50}$ MGs as a starting point, and ternary component Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ MGs, both known for their exceptional GFA, outstanding mechanical properties and significant potential in engineering applications [26, 27, 28, 29, 9].

Recently, the limitation of timescale of standard simulations has been overcome through the development of a very efficient method, called Swap Monte Carlo (SMC)  [30, 31, 32]. Here the equilibration time increases less rapidly than the relaxation time $\tau_{\alpha}$ as the temperature decreases [33, 34], allowing for the computation of the equilibrium configurational entropy below the experimental $T_{g}$ [35, 36, 37]. The SMC algorithm has seemed to be the sole means to investigate the properties of ultrastable glasses in silico. Despite the accelerated dynamics facilitated by swap moves, the crystallization effects commonly observed in conventional simulations continue to occur in SMC, but can be resolved by enforcing a degree of polydispersity, which effectively prevents crystallization and enables the generation of deeply cooled equilibrium liquid configurations, even at temperatures below the glass transition $T_{g}$, comparable with experiments. So far SMC has been applied only to a simplify model of metallic glass, a ternary mixture Lennard-Jones, that does not consider the specific nature of the constituents [38], therefore limiting our knowledge of the chemical composition landscape for MGs.

An intriguing solution to accelerate the dynamics sampling in the chemical configuration space [39] comes from the field of multicomponent crystalline alloys and combines MC methods, MD, and realistic short range potential, e.g. Embedded Atom Method (EAM) [40]. Only very recently, this techniques has been applied in glasses to obtain more stable realistic samples of binary CuZr MG in computer simulations [41], where it was shown that shear transformation zones are limited to small clusters of particles. Here, we demonstrate how this approach is suitable for tackling various concentrations of different species and that it is applicable to produce more stable realistic samples of multicomponent MGs in silico. We introduce for the first time equilibrium configurations of ternary metallic MGs, specifically Cu${}_{46}$Zr${}_{46}$Al${}_{8}$, approaching experiments. We examine their thermodynamic and more stable states and compare the mechanical properties. Our work presents a first example of generating computationally realistic models for any type of metallic glass consisting of several atomic species.

Simulation methods — The interatomic interactions are simulated through an Embedded Atom Method (EAM) as developed by [27]. For thermodynamics and kinetics, we performed simulations consisting of $N$=1500 atoms in a cubic box with periodic boundary conditions in three dimensions. Our simulations are performed with LAMMPS [42], using a time step $\Delta t$=1 fs. The glass state is obtained through quenching in the isobaric-isothermal ensemble ($NpT$) from the liquid at 2000 K to 300 K using a cooling rate $\kappa=10^{12}$ K/s. We save configurations every 50 K during the cooling process for further analysis. The cooling process is performed by integrating the Nosé-Hoover equations with damping parameters $\tau_{T}$=2 fs and $\tau_{p}$=5 ps for the thermostat and barostat. All results are obtained keeping the external pressure $p$=0.
Hybrid Molecular Dynamics-Monte Carlo (MD+MC) algorithm — In order to generate glasses that attempt to emulate the mechanical behavior observed in experiments, we employ a hybrid Molecular Dynamics-Monte Carlo (MD+MC) scheme under the variance-constrained semi-grand canonical ensemble (VC-SGC) [39]. Differently from traditional SMC, this VC-SGC MC scheme allows exploring the configurational degrees of freedom by randomly selecting an atom and attempting to change its type, while also calculating the corresponding energy and concentration changes. It allows to target specific concentration ranges while maintaining fixed total number of particles and the volume. Acceptance of these transmutations follows the Metropolis criterion, ensuring the preservation of detailed balance. On the other hand, the relaxation processes are accounted for by the MD integration steps. To maintain the desired concentration within the system [39], we set the variance parameter $\kappa$=10${}^{3}$. Furthermore, we evaluate the differences in chemical potential relative to Zr using hybrid MD+MC simulations under the semi-grand canonical ensemble at a temperature of 2000 K. The specific set of parameters that minimize the composition errors in relation to the desired concentration can be found in the Appendix.
Quenching Protocol — The quenching is performed starting from molten metals at high temperature well above the melting point using a fixed cooling rate of $10^{12}$ K/s at fixed pressure. For the MD+MC scheme, every 20 MD steps, a MC cycle consisting of $N/4$ attempts is performed. For a comparison between different cooling rates with MD and MD+MC, see Appendix.

Results — First, we compare the thermal evolution of the potential energy per particle between a conventional MD simulation and our MD+MC scheme for the two MG systems. Figure 1 reports the potential energy per atom ($U/N$) as a function of the temperature $T$ during quenching, obtained by averaging over 5 independent simulations.

Figure 1: Thermal evolution of the potential energy per atom of the binary Cu${}_{50}$Zr${}_{50}$ and the ternary Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ MGs with a cooling rate of $10^{12}$ K/s using both a pure MD and a hybrid MD+MC algorithm (red lines). At high temperatures, both methods exhibit similar results. However, as the temperature decreases, the behavior diverges significantly. The MC+MD method remains in equilibrium for lower temperatures compared to the pure MD approach.

Notably, the hybrid algorithm allows us to quench the liquid at lower temperatures, increasing the chance of reaching more stable configurations and properties similar to those observed in experiments. Note that both methods, MD and MC+MD, exhibit the same trend at high temperatures, as expected.
After the quenching we evaluate the relaxation time of the two different dynamics: conventional MD and MC+MD, for the ternary CuZrAl system. For each temperature we switch from $NpT$ to $NVT$ ensemble and run equilibration runs to compute the self-part of the intermediate scattering function [32]. For the case of MC+MD we employ the VC-SGC method as described above. The relaxation time, $\tau_{\alpha}$, defined as the value at which the intermediate scattering function is equal to 1/e, is reported in Fig. 2 for the standard MD (circles) and MC+MD dynamics (squares).

Figure 2: Comparison of relaxation time $\tau_{\alpha}$ for standard MD simulations and hybrid MD+MC with VC-SGC method for the ternary CuZrAl MG. The speedup offered by the MD+MC dynamics is evident.

The speedup of the dynamics given by our MD+MC scheme is shown as about two orders of magnitude at the lowest temperature.

It is well known that the GFA of CuZr MG improves when a small percentage of Al is included [7]. In a previous study of the same MG [43], we estimated $T^{MD}_{g}$ to be $623$ K and $713$ K for Cu${}_{50}$Zr${}_{50}$ and Cu${}_{46}$Zr${}_{46}$Al${}_{8}$, respectively. By following a similar protocol, we observe the same trend in the MD+MC algorithm, indicating that the method mimics the melt-quench procedure but provides equilibrium configurations for temperatures below the conventional $T_{g}$, meaning that it creates more stable glasses. Throughout the quenching process, in both pure MD and MD+MC, we carefully analyzed the structure and found no evidence of crystallization in the alloys.

Thermodynamics — Once we successfully create glasses approaching experimental conditions, our focus shifts to exploring their thermodynamic properties. To being with, we determine the configurational entropy of the glass formed by MD+MC by the formula $S_{conf}=S_{tot}-S_{vib}$, where $S_{tot}$ represents the total entropy and $S_{vib}$ stands for the vibrational contribution. We obtain both entropies using the reversible scaling (RS) method [44], with the Uhlenbeck-Ford potential and the Einstein crystal serving as reference systems, respectively [45, 46, 47] (see Supplementary Information). Figure 3 displays the configurational entropy $(S_{conf})$ plotted against temperature for both the binary Cu${}_{50}$Zr${}_{50}$ and ternary Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ alloys. As anticipated, $S_{conf}$ decreases with decreasing temperature until it freezes in at lower temperatures due to the glass transition. The measured residual entropy values, representing the constant $S_{conf}$ value at the occurrence of the glass transition, were 0.049 and 0.035 J/mol-K for Cu${}_{50}$Zr${}_{50}$ and Cu${}_{46}$Zr${}_{46}$Al${}_{8}$, respectively. Note in particular that the $S_{conf}$ of Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ are lower than the binary alloy, the reason for that is because the ternary alloy possesses a high GFA allowing to reach deeper values before the freeze-in at $T_{g}$ (which is higher). The reported value of $S_{conf}$ in the Cu${}_{50}$Zr${}_{50}$ alloy is just a $20\%$ from the one reported in situ by Smith et al. [48]. We are dealing with glass configurations with glass transition temperatures lower than the one obtained in conventional glasses, thus it is expected to have a lower value of the residual entropy. A comparison with experiments for the ternary system is currently not available.

Figure 3: Configurational entropy as a function of temperature obtained through RS method for Cu${}_{50}$Zr${}_{50}$ (blue) and Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ (red).

Figure 4: Heat capacity as a function of temperature obtained through RS method for $a)$ Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ and $b)$ Cu${}_{50}$Zr${}_{50}$. The insets show the potential energy of the binary and ternary alloys, respectively, as a function of temperature during the heating process using MD (thin lines) and MD+MC (dashed lines).

Following experiments on ultra stable MGs using vapor-deposition method, a common approach to studying their stability is by performing calorimetric measurements during a heating process. In Figure 4, we compare the behavior of glasses prepared using conventional MD with those created using the MD+MC algorithm. We heat the glass at the same rate of $10^{12}$ K/s and save the configuration every $20$ K. Then, we compute the heat capacity using its statistical definition: $(\langle E^{2}\rangle-\langle E\rangle^{2})/k_{B}T^{2}$, where $E$ is the energy of the system, and $k_{B}$ is the Boltzmann constant. The insets of Figures 4$a)$ and 4$b)$ with the potential energy of the binary and ternary alloys, respectively, illustrate as well via the difference of the two kinds of glasses as a function of temperature. Clearly, the devitrification process of the glass created through the hybrid MD+MC happens at a higher temperature, approximately $18\%$ above $T_{g}$ [49, 50, 15]. In vapor-deposited experiments [7], this temperature is called the “onset” temperature $T_{o}$, and it is generally between $6$ to $10\%$ above the conventional $T_{g}$, reflecting the much larger kinetic stability reached by this method. Since our MD+MC glass produces similar results, we can confirm that the Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ and Cu${}_{50}$Zr${}_{50}$ alloys produced by this method are more stable configurations.

Rheology — For mechanical tests we perform athermal quasistatic simulations [51] starting from configurations with $N$=100000 atoms, at $T$=300 K cooled down to 0 K, following the work of Ref. [52]. We shear the simulation box along the $x$ direction by the amount $\delta\gamma$. We perform energy minimization after each strain step with the fast inertial relaxation engine (FIRE) [53] until the system reaches mechanical equilibrium, until a maximum force of 10${}^{-10}$eV/Å  is reached.

Figure 5: Scaled shear-modulus stress-strain curves for MD and MD+MC schemes (upper panel). Strain inhomogeneity maps for a strain of 0.125 for the CuZrAl cases. The color bar indicates the non affine displacements from the unstrained configuration, from 0% (blue) to 20% (red).

Figure 5 (top) shows how the composition and the history influence the mechanical behavior. The addition of aluminum increases by a few percent the shear modulus [54], and to compare the glasses at peak stress we scale away the dependence of on the modulus. These are typical stress-strain curves, where after an initial elastic part we see the typical glass response: strain-softening close to the maximum, fluctuations, and a drop in stress before plastic flow. In particular, as in model glasses [41] we observe that the more stable glasses exhibit drastic stress drops which corresponds to the the appearance of a system-spanning shear band (Inset). The well-prepared glasses do so at slightly larger strains than the basic ones, and in scaled units they are noticeably stronger, here about 20%. The final state of flow one finds may be summarized by the examples of the strain fields in the Fig. 5 (bottom), with the features associated with the post-yield shear bands depending on the nature of the stress-peak: for glasses with enhanced stability we find a well-defined shear band whereas for the reference glass the band is more diffuse [38].

Summary & Conclusions — In order to resolve the physics of glasses in silico - a big challenge due to the complexity of glasses and their history-dependence - novel ideas are needed. In this work, we show for the first time how a combination of Molecular Dynamics and Monte Carlo with VC-SGC method is able to tackle this issue for realistic multi-component MGs.

We present results that clearly show a convergence to the experimental observations as regards the thermodynamics of these systems, we also demonstrate, that such glasses have the ultrastable character, aligning them closely with vapor-deposited experimental counterparts. Finally, this fact is seen to make a difference in the mechanical properties. The next steps to follow are a through exploration of the relaxation times the MC+MD approach leads to and the subsequent changes in the yielding.

In conclusion, our results provide a substantial contribution to the field of glass physics, particularly with regard to multicomponent MGs with more than two component species [55, 56]. Using advanced simulations to obtain ultrastable glassy states, we provide an essential theoretical framework for further studies, shedding light on the underlying physical mechanisms governing the complex behaviors and properties of these amorphous materials. Our work prompts further theoretical investigations aimed at probing the microscopic mechanisms of complex MGs at a fundamental level, such as their relaxation dynamics and facilitation [57]. The methodologies employed in this work not only serve as a benchmark for experimental comparison but also suggest new avenues for optimizing the structure of MGs via cooling and processing techniques, for a range of applications.

Acknowledgements: We gratefully thank Anshul Parmar and Misaki Ozawa for fruitful discussions. Funding: RAD, SB, and MA are supported by the European Union Horizon 2020 research and innovation program under grant agreement no. 857470 and from the European Regional Development Fund via the Foundation for Polish Science International Research Agenda PLUS program grant No. MAB PLUS/2018/8.

Appendix A APPENDIX

Appendix B Setting VC-SGC parameters

In the variance-constrained semi-grand canonical (VC-SGC) ensemble[39], the system explores configurational degrees of freedom by randomly selecting an atom and attempting to change its type. The change in energy and concentration is then computed, and transmutations are accepted or rejected based on the Metropolis criterion, ensuring detailed balance and equilibrium sampling of the phase space. To determine the difference in chemical potential species $(\Delta\mu)$, which fixes the desired composition and minimizes error, we conduct a series of semi-grand-canonical (SGC) simulations at a temperature of $T=2000$ K.
Figure 6 depicts the behavior of $\Delta\mu$ with regard to Zr for $a)$ Cu and $b)$ Al for a series of SGC simulations. The blue line is obtained expanding $\Delta\mu$ as a polynomial in concentration as follows:

$$\Delta\mu(X_{\rm{Zr}},p,T)=T\ln\left(\frac{X_{\rm{Zr}}}{1-X_{\rm{Zr}}}\right)+\sum_{i=0}^{n}A_{i}X_{\rm{Zr}}^{i},$$ (1)

where $X_{\rm{Zr}}$ is the Zr composition and the coefficient $A_{i}$ are fitting parameters.

Figure 6: Chemical potential as a function of the composition for $a)$ Zr-Cu and $b)$ Zr-Al from SGC simulations at $T=2000$ K. The blue solid line is a fitting obtained by Eq. 1.

By utilizing Eq (1), we can effectively determine the chemical potential at any desired composition. For both the binary Cu${}_{50}$Zr${}_{50}$ and ternary Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ alloys, the calculated chemical potential differences were $\Delta\mu_{\rm Zr-Cu}=-2.69854$ eV, and $\Delta\mu_{\rm Zr-Cu}=-2.90351$ eV, $\Delta\mu_{\rm Zr-Al}=-0.55781$ eV, respectively. Once the composition is defined, we proceed with the hybrid MC+MD simulation under the VC-SGC ensemble. This combination of MD+MC efficiently allows us to create in silico glasses with properties similar to the ones obtained in experiments.

Appendix C Cooling rate comparison

Figure 7: Potential energy for various cooling rates for the MD vs. the hybrid algorithm for $a)$ Zr-Cu and $b)$ Zr-Al.

In Figure 7, we compare the thermalization efficiency of the hybrid MC+MD algorithm with conventional MD simulations. Starting from high temperatures well above the melting point, we rapidly cool the liquid using fixed cooling rates of $10^{14}$, $10^{13}$, and $10^{12}$ through MD simulation. With the hybrid MC+MD approach, we observe the same cooling rate-dependent trend at higher temperatures. However, when employing the MC+MD algorithm, the liquid deviates from equilibrium at lower temperatures. The agreement among all the curves at high temperatures confirms the proper utilization of the MC+MD approach in our simulations.

Appendix D Computing the configurational entropy

We have numerically calculated the entropy by means of the thermodynamic relation $S=(H-G)/T$, where H and G are the enthalpy and Gibbs free energy of the system, respectively. Since free energy is a thermal property [58] we are not able to compute it directly from convectional MD simulations. Thus, we have used a combination of the adiabatic switching (AS)[59, 60] and reversible scaling (RS)[61, 62] methods, which in principle allows the estimation of the free energy for a wide range of temperatures in a single simulation. We estimated the free energy (and the entropy) by carrying out the following steps:

• The preparation of the initial state in which we want to calculate the free energy. This state is defined as:

  $$\mathcal{H}_{\rm MEAM}(\Gamma_{\alpha})=\mathcal{T}(\Gamma_{\alpha})+U_{\rm MEAM}(\Gamma_{\alpha}),$$ (2)

  where $\mathcal{H}_{\rm MEAM}$ is the Hamiltonian of the alloy described by EAM potential, $\mathcal{T}$ and $U_{\rm EAM}$ are the kinetic and potential energy, respectively. $\Gamma_{\alpha}$ is a phase space vector with information of all momenta and positions as a function of time, and the subindex $\alpha$ holds to crystal $(g)$ or liquid $(l)$.

• Once defined the initial state, we estimate its absolute free energy through AS method. The AS is a variation of the thermodynamic integration method[63] where the Hamiltonian of the initial state is coupled with a reference system whose free energy is known beforehand. However, the work done is obtained dynamically in one single simulation. Thus, the free energy is obtained by means of:

  $$G_{\rm EAM}^{0}=G_{Ref}^{\alpha}+\mathcal{W}_{AS},$$ (3)

  where $G_{\rm EAM}^{0}$ is the free energy of the alloy in the initial state, $G_{Ref}^{\alpha}$ is the free energy of the reference system, and $\mathcal{W}_{AS}$ is the quasistatic work obtained through AS-method.

• Using the value of $G_{\rm EAM}^{0}$ and the RS method we obtain the alloy-free energy for a wide range of temperatures. This method is based on writing the initial Hamiltonian in the parametric form[61]:

  $$\mathcal{H}_{\rm EAM}(\Gamma_{\alpha},\zeta(t))=\mathcal{T}+\zeta(t)U_{\rm EAM}.$$ (4)

  Here each value of $\zeta$ corresponds to a particular state of $H_{\rm EAM}$ in a particular temperature and we use $\mathcal{T}$ and $U_{\rm EAM}$ instead of $\mathcal{T}(\Gamma_{\alpha})$ and $U_{\rm EAM}(\Gamma_{\alpha})$ to simplify the notation. One can show that using Eq. (4) the free energy of the alloy at any temperature is given by the equation[61, 44]:

  $$\begin{split}G_{\rm EAM}(T)=\frac{G_{\rm EAM}^{0}}{\zeta}+\frac{3}{2}nk_{B}T_{0}\frac{\ln\zeta}{\zeta}+\\ +\frac{1}{\zeta}\int_{0}^{t_{sim}}dt\frac{d\zeta}{dt}U_{\rm EAM},\end{split}$$ (5)

  where $T_{0}$ is the temperature at which $G^{0}_{\rm EAM}$ was estimated in item (ii) and t sim is the simulation time.

• Finally, the configurational entropy is obtained as:

  $$S_{conf}=S_{l}-S_{g},$$ (6)

  where $S_{l}$ and $S_{g}$ are the liquid and glass entropies, respectively.

The reference free energy $G^{0}_{\rm EAM}$ in Eq. (3) is obtained by transforming the MEAM interaction into a reference system by means of AS method. For the solid phase at $T=300$ K, we used the Einstein crystal with a spring constant $\alpha=5eV/\AA^{2}$ for all the particles of the system[46, 64]. The liquid phase, on the other hand, $G^{0}_{\rm EAM}$ was obtained at $T=2800$ K using as a reference the Uhlenbeck-Ford model (UFM) defined as[45, 65]:

$$U_{\rm UFM}(r)=-\frac{p}{\beta}\ln\left(1-e^{-(r/\sigma)^{2}}\right),$$ (7)

where $\beta=1/(k_{B}T)$ is the inverse of the temperature times the Boltzmann constant $k_{B}$ , and $p$ and $\sigma$ are the softness and length scaling parameters, respectively. Following Paula-Leite et al.[45] a value of $p=50$ gives a liquid-like behavior for the UFM system. The quasistatic work $W_{AS}$ is estimating and evaluating the irreversible work during two adiabatic switching trajectories. One transforming the system from the EAM to the reference system ($W^{irr}_{\rm EAM\rightarrow ref}$) and the other one in the inverse process ($W^{irr}_{\rm ref\rightarrow EAM}$). The irreversible work done in the AS method is calculated by $W^{irr}=\int_{0}^{t_{sim}}dt\frac{d\lambda}{dt}U(\Gamma_{\alpha})$, where $U(\Gamma_{\alpha})$ is the hybrid potential connecting $U_{\rm EAM}$ and $U_{ref}$ and $\lambda$ is the AS coupling parameter. Thus,

$$\mathcal{W}_{AS}=\frac{1}{2}\left(W_{\rm EAM\rightarrow ref}^{irr}-W_{\rm rsf\rightarrow EAM}^{irr}\right).a$$ (8)

Equation (8) is only valid if the dynamical work is performed slowly enough for holding the system in the linear response regime[44, 61].

Figure 8: Free energy and Entropy as a function of temperature obtained by the RS method of the ternary Cu${}_{46}$Zr${}_{46}$Al${}_{8}$ alloy.

In Fig. 8 we show the absolute values of the free energy and the entropy obtained by the RS method of the ternary Cu${}_{46}$Zr${}_{46}$Al${}_{8}$.