After thermalization, we compute the corresponding inherent structures through
the Limited-memory Broyden-Fletcher-Goldfarb-Shanno
algorithm [1]}.
Let \( {\mathbf {r}} \)  be a configuration of the system, i.e., \( {\mathbf {r}} \equiv ( {\mathbf {r}} _1,\ldots , {\mathbf {r}} _N)\) .
The mechanical energy of the configuration \( {\mathbf {r}} \)  is \(E[ {\mathbf {r}} ]=\sum _{i<j}\phi (r_{ij})\) .
We indicate
with \( {\mathbf {r}} ^0\equiv ( {\mathbf {r}} _1^0,\ldots , {\mathbf {r}} _N^0)\)  the configuration that minimizes \(E[ {\mathbf {r}} ]\) ,
We define the Inherent Structure energy \(e_\text{IS} \equiv N^{-1} \langle E[ {\mathbf {r}} ^0] \rangle _s\) , with
\(\langle \dots \rangle _s\)  indicating the average over \(s\)  independent configurations.
The temperature dependence of \(e_\text{IS}\)  is shown in figure REF .
