Recent research has been constraining the retention fraction of black holes (BHs) in globular clusters by comparing the degree of mass segregation with $N$-body simulations. They are consistent with an upper limit of the retention fraction being $50\,\%$ or less. In this work, we focus on direct simulations of the dynamics of BHs in star clusters. We aim to constrain the effective distribution of natal kicks that BHs receive during supernova (SN) explosions and to estimate the BH retention fraction. We used the collisional $N$-body code nbody6 to measure the retention fraction of BHs for a given set of parameters, which are: the initial mass of a star cluster, the initial half-mass radius, and $\sigma_\mathrm{BH}$, which sets the effective Maxwellian BH velocity kick distribution. We compare these direct $N$-body models with our analytic estimates and newest observational constraints. The numerical simulations show that for the one-dimensional (1D) velocity kick dispersion $\sigma_\mathrm{BH} < 50\,\mathrm{km\,s^{-1}}$, clusters with radii of 2 pc and that are initially more massive than $5 \times 10^3\,M_\odot$ retain more than $20\,\%$ of BHs within their half-mass radii. Our simple analytic model yields a number of retained BHs that is in good agreement with the $N$-body models. Furthermore, the analytic estimates show that ultra-compact dwarf galaxies (UCDs) should have retained more than $80\,\%$ of their BHs for $\sigma_\mathrm{BH} \leq 190\,\mathrm{km\,s^{-1}}$. Although our models do not contain primordial binaries, in the most compact clusters with $10^3$ stars, we have found evidence of delayed SN explosions producing a surplus of BHs compared to the IMF due to dynamically formed binary stars. These cases do not occur in the more populous or expanded clusters.

We analyse spectroscopic observations of the B[e] star HD 50138 (MWC 158, V743 Mon, or IRAS 06491-0654), a member of the FS CMa group, obtained over the last twenty years. Four different epochs are identified in the observational data, where the variability of the spectral features is substantially different. Additionally, two long periods of (3 000 +/- 500) and (5 000 +/- 1000) days are found in the variations of the equivalent widths of the H alpha and [OI] 6300 A lines and radial velocities of the H alpha line violet peak. Modest signatures of a regular period of ~34 days in the radial velocities of the H alpha red peak and H beta central depression are found in the season 2013/2014. The H alpha V/R changes indicate a periodicity of ~50 days. The correlations between individual spectral features significantly restricts the model of the object and suggest that it is most likely a binary system with a highly distorted disc with spiral arms around the primary component. At the same time, no obvious signs of the secondary component has been found in the object's spectrum.