  Partition function zeros steer the critical behavior of a system. Studying
four-flavor lattice QCD at finite temperature and density with the
Wilson-clover fermion action and the Iwasaki gauge action using a
phase-quenched fermion determinant, we combine statistics from multiple
chemical potentials to improve sampling of the configuration space, and aim at
unraveling the movement of zeros in finite systems. Preparing for further
investigations, we discuss methods and criteria used to sieve through complex
parameter space spanned by $(\text{Re}\mu, \text{Im}\mu)$ and $(\text{Re}\mu,
\text{Im}\beta)$, and present statistically robust zeros of the partition
function.
