  One of the most challenging issues in QCD is the investigation of spontaneous
chiral-symmetry breaking, which is characterized by the non-vanishing chiral
condensate when the bare fermion mass is zero. In standard methods, one has to
perform expensive lattice simulations at multiple bare quark masses, and employ
some modeled function to extrapolate the data to the chiral limit. This paper
applies the probability distribution function method to computing the chiral
condensate in lattice QCD with massless dynamical quarks, without any ambiguous
mass extrapolation. The results for staggered quarks indicates that the method
might be a more efficient alternative for investigating the spontaneous
chiral-symmetry breaking in lattice QCD.
