Thermodynamical quantities of lattice full QCD from an efficient method

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 $10^4$ lattice at flavor number $N_f=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.