We consider the current experimental constraints on the parameter space of
the MSSM and NMSSM. Then in the allowed parameter space we examine the Higgs
pair production at the 14 TeV LHC via bbˉ→hh (h is the 125 GeV
SM-like Higg boson) with one-loop SUSY QCD correction and compare it with the
production via gg→hh. We obtain the following observations: (i) For the
MSSM the production rate of bbˉ→hh can reach 50 fb and thus can be
competitive with gg→hh, while for the NMSSM bbˉ→hh has a much
smaller rate than gg→hh due to the suppression of the hbbˉ
coupling; (ii) The SUSY-QCD correction to bbˉ→hh is sizable, which
can reach 45% for the MSSM and 15% for the NMSSM within the 1σ
region of the Higgs data; (iii) In the heavy SUSY limit (all soft mass
parameters become heavy), the SUSY effects decouple rather slowly from the
Higgs pair production (especially the gg→hh process), which, for MSUSY​=5 TeV and mA​<1 TeV, can enhance the production rate by a factor of
1.5 and 1.3 for the MSSM and NMSSM, respectively. So, the Higgs pair production
may be helpful for unraveling the effects of heavy SUSY.Comment: discussions and references added, accepted by JHE