  In light of the recent LHC Higgs search data, we investigate the pair
production of a SM-like Higgs boson around 125 GeV in the MSSM and NMSSM. We
first scan the parameter space of each model by considering various
experimental constraints, and then calculate the Higgs pair production rate in
the allowed parameter space. We find that in most cases the dominant
contribution to the Higgs pair production comes from the gluon fusion process
and the production rate can be greatly enhanced, maximally 10 times larger than
the SM prediction (even for a TeV-scale stop the production rate can still be
enhanced by a factor of 1.3). We also calculate the chi-square value with the
current Higgs data and find that in the most favored parameter region the
production rate is enhanced by a factor of 1.45 in the MSSM, while in the NMSSM
the production rate can be enhanced or suppressed (\sigma_{SUSY}/\sigma_{SM}
varies from 0.7 to 2.4).
