Computing the topological susceptibility from fixed topology QCD simulations

Abstract: The topological susceptibility is an important quantity in QCD, which can be computed using lattice methods. However, at a fine lattice spacing, or when using high quality chirally symmetric quarks, algorithms which proceed in small update steps --- in particular the HMC algorithm --- tend to get stuck in a single topological sector. In such cases, the computation of the topological susceptibility is not straightforward. Here, we explore two methods to extract the topological susceptibility from lattice QCD simulations restricted to a single topological sector. The first method is based on the correlation function of the topological charge density, while the second method relies on measuring the topological charge within spacetime subvolumes. Numerical results for two-flavor QCD obtained by using both methods are presented.