Gas–surface interactions play important roles in internal rarefied gas flows, especially in micro-electro-mechanical systems with large surface area to volume ratios. Although great progress has been made to solve the Boltzmann equation, the gas kinetic boundary condition (BC) has not been well studied. Here we assess the accuracy of the Maxwell, Epstein and Cercignani–Lampis BCs, by comparing numerical results of the Boltzmann equation for the Lennard–Jones potential to experimental data on Poiseuille and thermal transpiration flows. The four experiments considered are: Ewart et al. (J. Fluid Mech., vol. 584, 2007, pp. 337–356), Rojas-Cárdenas et al. (Phys. Fluids, vol. 25, 2013, 072002) and Yamaguchi et al. (J. Fluid Mech., vol. 744, 2014, pp. 169–182; vol. 795, 2016, pp. 690–707), where the mass flow rates in Poiseuille and thermal transpiration flows are measured. This requires that the BC has the ability to tune the effective viscous and thermal slip coefficients to match the experimental data. Among the three BCs, the Epstein BC has more flexibility to adjust the two slip coefficients, and hence for most of the time it gives good agreement with the experimental measurements. However, like the Maxwell BC, the viscous slip coefficient in the Epstein BC cannot be smaller than unity but the Cercignani–Lampis BC can. Therefore, we propose to combine the Epstein and Cercignani–Lampis BCs to describe gas–surface interaction. Although the new BC contains six free parameters, our approximate analytical expressions for the viscous and thermal slip coefficients provide useful guidance to choose these parameters.