SUSY effects in Higgs productions at high energy $e^+e^-$ colliders

Considering the constraints from collider experiments and dark matter detections, we investigate the SUSY effects in the Higgs productions $e^+e^- \to Zh$ at an $e^+e^-$ collider with a center-of-mass energy above 240 GeV and $\gamma\gamma \to h \to b\bar{b}$ at a photon collider with a center-of-mass energy above 125 GeV. In the parameter space allowed by current experiments, we find that the SUSY corrections to $e^+e^- \to Zh$ can reach a few percent and the production rate of $\gamma\gamma \to h \to b\bar{b}$ can be enhanced by a factor of 1.2 over the SM prediction. We also calculate the exotic Higgs productions $e^+e^-\to Zh_1$ and $e^+e^-\rightarrow A_1h$ in the next-to-minimal supersymmetric model (NMSSM) ($h$ is the SM-like Higgs, $h_1$ and $A_1$ are respectively the CP-even and CP-odd singlet-dominant Higgs bosons which can be much lighter than $h$). We find that at a 250 GeV $e^+e^-$ collider the production rates of $e^+e^-\rightarrow Zh_1$ and $e^+e^-\to A_1h$ can reach 60 fb and 0.1 fb, respectively.