Measuring the topological susceptibility in a fixed sector

For field theories with a topological charge Q , it is often of interest to measure the topological susceptibility χ_t=(⟨Q²⟩-⟨Q ⟩²)/V . If we manage to perform a Monte Carlo simulation where Q changes frequently, χ_t can be evaluated directly. However, for local update algorithms and fine lattices, the autocorrelation time with respect to Q tends to be extremely long, which invalidates the direct approach. Nevertheless, the measurement of χ_t is still feasible, even when the entire Markov chain is topologically frozen. We test a method for this purpose, based on the correlation of the topological charge density, as suggested by Aoki, Fukaya, Hashimoto and Onogi. Our studies in nonlinear σ -models and in two-dimensional Abelian gauge theory yield accurate results for χ_t, which confirm that the method is applicable. We also obtain promising results in four-dimensional SU(2) Yang-Mills theory, which suggest the applicability of this method in QCD.