The linear stability of the axisymmetric steady thermocapillary flow in a liquid bridge made from 2 cSt silicone oil (Prandtl number 28) is investigated numerically in the framework of the Boussinesq approximation. The flow and temperature fields in the surrounding gas phase (air) are taken into account for a generic cylindrical container hosting the liquid bridge. The flows in the liquid and in the gas are fully coupled across the hydrostatically deformed liquid–gas interface, neglecting dynamic interface deformations. Originating from a common reference case, the linear stability boundary is computed varying the length of the liquid bridge (aspect ratio), its volume and the gravity level, providing accurate critical data. The qualitative dependence of the critical threshold on these parameters is explained in terms of the characteristics of the critical mode. The heat exchange between the ambient gas and the liquid bridge that is fully resolved has an important influence on the critical conditions.