Using an evolutionary population synthesis code, we modeled the universal, featureless X-ray luminosity function of high-mass X-ray binaries (HMXBs) in star-forming galaxies. We put constraints on the natal kicks, super-Eddington accretion factor, as well as common envelope prescriptions usually adopted (i.e., the αCE formalism and the γ algorithm), and presented the detailed properties of HMXBs under different models, which may be investigated further by future high-resolution X-ray and optical observations.