I extend to QCD an efficient method for lattice gauge theory with dynamical
fermions. Once the eigenvalues of the Dirac operator and the density of states
of pure gluonic configurations at a set of plaquette energies (proportional to
the gauge action) are computed, thermodynamical quantities deriving from the
partition function can be obtained for arbitrary flavor number, quark masses
and wide range of coupling constants, without additional computational cost.
Results for the chiral condensate and gauge action are presented on the 104
lattice at flavor number Nf​=0, 1, 2, 3, 4 and many quark masses and coupling
constants. New results in the chiral limit for the gauge action and its
correlation with the chiral condensate, which are useful for analyzing the QCD
chiral phase structure, are also provided.Comment: Latex, 11 figures, version accepted for publicatio