“…As the selection of reliable interatomic potentials for classical, semi-empirical potential molecular dynamics (MD) is challenging due to the anharmonicity in metallic glasses (Lambson et al, 1986;Crespo et al, 2016;Aitken et al, 2018), density functional theorybased ab initio MD has been used to calculate the properties of a vast number of materials as it does not require the input of interatomic potentials while being able to cover many-body interactions. By ab initio modeling, an in-depth analysis of the structure-property relationship of metallic glasses is possible that is otherwise challenging to probe only by experimental means, including parameters such as the partial pair correlations (Qin et al, 2007;Ganesh and Widom, 2008;Hui et al, 2008bHui et al, , 2009Fang et al, 2009;Tian et al, 2011;Zhang et al, 2015;Hunca et al, 2016;Yu et al, 2016), coordination polyhedra (Ganesh and Widom, 2008;Fang et al, 2009;Fujita et al, 2009;Hui et al, 2009;Hirata et al, 2011;Tian et al, 2011;Durandurdu, 2012;Wu et al, 2012;Zhang et al, 2015;Huang et al, 2016;Yu et al, 2016) and bonding Hunca et al, 2016). However, modeling the dynamics involved in relaxation and rejuvenation phenomena is out of reach for ab initio molecular dynamics due to the limitations of time steps and cell size as discussed below.…”