Velocity distribution of stars within cluster complex “NGC6334” and 2371 M_⊙ cluster under different gas-expulsion modes at 20 Myr. (a) Probability density function (PDF) of stellar velocities, estimated using the kernel-density-estimation (KDE) method to provide a smooth representation of the distribution; (b) cumulative fraction shows the proportion of stars with velocities below a given value. Here, we do not distinguish between binaries and single stars; instead, we calculated the total velocity of each star (v_tot) in the cluster and examined their distribution to roughly characterize the expansion state of the cluster. At 50 Myr, the velocity distributions remain unchanged, indicating that the binary fraction does not significantly affect the statistical results here.