Using molecular dynamics simulations, the present study investigated the precise characteristics of the binary mixture of condensable gas (vapor) and non-condensable gas (NC gas) molecules creating kinetic boundary conditions (KBCs) at a gas–liquid interface in equilibrium. We counted the molecules utilizing the improved two-boundary method proposed in previous studies by Kobayashi et al. (Heat Mass Trans 52:1851–1859, 2016. doi:10.1007/s00231-015-1700-6). In this study, we employed Ar for the vapor molecules, and Ne for the NC gas molecules. The present method allowed us to count easily the evaporating, condensing, degassing, dissolving, and reflecting molecules in order to investigate the detailed motion of the molecules, and also to evaluate the velocity distribution function of the KBCs at the interface. Our results showed that the evaporation and condensation coefficients for vapor and NC gas molecules decrease with the increase in the molar fraction of the NC gas molecules in the liquid. We also found that the KBCs can be specified as a function of the molar fraction and liquid temperature. Furthermore, we discussed the method to construct the KBCs of vapor and NC gas molecules.