Depletion of two-level systems in ultrastable computer-generated glasses

Amorphous solids exhibit quasiuniversal low temperature anomalies whose origin has been ascribed to localized tunneling defects. Using an advanced Monte Carlo procedure, we create in silico glasses spanning from hyperquenched to ultrastable glasses. Using a multidimensional path-fmding protocol, we locate tunneling defects with energy splittings smaller than k(B) T-Q, with T-Q the temperature below which quantum effects are relevant (T-Q approximate to 1 K in most experiments). We find that as the stability of a glass increases, its energy landscape as well as the manner in which it is probed tend to deplete the density of tunneling defects, as observed in recent experiments. We explore the real-space nature of tunneling defects, and find that they are mostly localized to a few atoms, but are occasionally dramatically delocalized.