  The dynamical fluctuation and correlation of multiplicity distributions of
identified baryons and antibaryons produced by the hadronization of the bulk
quark system are systematically studied in quark combination model. Starting
from the most basic dynamics of quark combination which are necessary for
multiplicity study, we analyze moments (variance, skewness and kurtosis) of
inclusive multiplicity distribution of identified baryons, two-baryon
multiplicity correlations, and baryon-antibaryon multiplicity correlations
after the hadronization of quark system with given quark number and antiquark
number. We obtain a series of interesting findings, e.g., binomial behavior of
multiplicity moments, coincide flavor dependent two-baryon correlation and
universal baryon-antibaryon correlation, which can be regarded as general
features of quark combination. We further take into account correlations and
fluctuations of quark numbers before hadronization to study their influence on
multiple production of baryons and antibaryons. We find that quark number
fluctuation and flavor conservation lead to a series of important results such
as the negative $p\bar{\Omega}^{+}$ multiplicity correlation and universal
two-baryon correlations. We also study the influence of resonance decay in
order to compare our results with future experimental data in
ultra-relativistic heavy ion collisions at LHC.
