In the present paper, we provide an analytical expression for the first- and second-order thermal slip coefficients, σ1,T and σ2,T, by means of a variational technique that applies to the integrodifferential form of the Boltzmann equation based on the true linearized collision operator for hard-sphere molecules. The Cercignani-Lampis scattering kernel of the gas-surface interaction has been considered in order to take into account the influence of the accommodation coefficients (αt, αn) on the slip parameters. Comparing our theoretical results with recent experimental data on the mass flow rate and the slip coefficient for five noble gases (helium, neon, argon, krypton, and xenon), we found out that there is a continuous set of values for the pair (αt, αn) which leads to the same thermal slip parameters. To uniquely determine the accommodation coefficients, we took into account a further series of measurements carried out with the same experimental apparatus, where the thermal molecular pressure exponent γ has been also evaluated. Therefore, the new method proposed in the present work for extracting the accommodation coefficients relies on two steps. First of all, since γ mainly depends on αt, we fix the tangential momentum accommodation coefficient in such a way as to obtain a fair agreement between theoretical and experimental results. Then, among the multiple pairs of variational solutions for (αt, αn), giving the same thermal slip coefficients (chosen to closely approximate the measurements), we select the unique pair with the previously determined value of αt. The analysis carried out in the present work confirms that both accommodation coefficients increase by increasing the molecular weight of the considered gases, as already highlighted in the literature.