  We derive exact, universal, closed-form quantum Monte Carlo estimators for
finite temperature energy susceptibility and fidelity susceptibility for
essentially arbitrary Hamiltonians. We demonstrate how our method can be used
to study quantum phase transitions without knowledge of an order parameter and
without the need to design system-specific ergodic quantum Monte Carlo update
rules by applying it to a class of random models.
