Full Counting Statistics for Interacting Fermions with Determinantal Quantum Monte Carlo Simulations

We present a method for computing the full probability distribution function of quadratic observables for the Fermi-Hubbard model within the framework of determinantal quantum Monte Carlo. Especially, in cold atoms experiments with single site resolution, such full counting statistics can be obtained from repeated projective measurements. We demonstrate, that the full counting statistics can provide important information on the size of preformed pairs. Furthermore, we compute the full counting statistics of the staggered magnetization in the repulsive Hubbard model at half filling and find excellent agreement with recent experimental results. We show that current experiments are capable of probing the difference between the Hubbard model and the limiting Heisenberg model.