X-ray luminosity function of accreting neutron stars and black holes

We model X-ray luminosity functions (XLF) of accreting neutron stars and black holes in $10^{35} \leq L_X \leq 10^{41}$ erg/s range in star-forming galaxies and galaxies with the initial star formation burst. XLFs are obtained by combining a fast generation of compact object+normal star population using the binary population synthesis code BSE and calculation of the subsequent detailed binary evolution by the MESA code. XLFs in the galaxies of both types is broadly reproduced using the standard assumptions of the binary star evolution.