How to calculate quantum quench distributions with a weighted Wang-Landau Monte Carlo

We present here an extension of the Wang-Landau Monte Carlo method which allows us to get very accurate estimates of the full probability distributions of several observables after a quantum quench for large systems, whenever the relevant matrix elements are calculable, but the full exponential complexity of the Hilbert space would make an exhaustive enumeration impossible beyond very limited sizes. We apply this method to quenches of free-fermion models with disorder, further corroborating the fact that a generalized Gibbs ensemble fails to capture the long-time average of many-body operators when disorder is present.