where \(\alpha , \beta \in \rm {A, B}\) . The units of various quantities in our simulation are as follows: lengths are expressed in the unit of \(\sigma _{AB}\) , energies in the unit of \(\epsilon _{AB}\) , time in the unit of \((m\sigma _{AB}^2/\epsilon _{AB})^{1/2}\)  and temperature in the unit of \(\epsilon _{AB}/k_{\rm B}\) . Here, \(k_{\rm B}\)  is the Boltzmann constant which is unity. The parameters \(\sigma _{\alpha \beta }\)  and \(\epsilon _{\alpha \beta }\)  are chosen as follows: \(\sigma _{AA} = 1.0, \sigma _{BB} = 0.88, \sigma _{AB} = 0.8\)  and \(\epsilon _{AA} = 1.0, \epsilon _{BB} = 0.5, \epsilon _{AB} = 1.5\) . With these set of parameters, a binary LJ mixture avoids crystallization below \(T_g\)  and forms stable glass. \(T_g\)  for this model is approximately \(0.44\)  in the reduced unit [1]}. To improve the computational efficiency, the potential is truncated at \(r_{cut}=2.5\sigma _{AA}\) . The parameters \(A_0 = 0.040490\) , \(A_1 = -0.009702\)  and \(A_2 = 0.000620\)  are such that the potential and its derivatives (up to second order) go to zero smoothly at \(r=r_{cut}\) .
