Predicting the photophoretic force exerted on an optical absorptive particle in a gaseous medium is a challenging problem because the problems of electromagnetic scattering, heat transfer, and gaseous molecule dynamics are involved and coupled with each other. Based on the calculation of the source function distribution inside a homogeneous sphere excited by a Bessel beam using the generalized Lorenz-Mie theory, analytical expressions of the asymmetry vector, which is the key quantity in the calculation of photophoretic force, are given using the adjoint boundary value method. Numerical simulations are performed to analyze the influences of polarization, the half-cone angle, and the beam order of the incident beam, particle size, and absorptivity of the particle on the asymmetry vector for both on-axis and off-axis illuminations. Longitudinal and transverse photophoretic forces on a homogeneous sphere are displayed for the slip-flow regime of gaseous media. The results offer important insights into the working mechanism underpinning the development of heat-mediated optical manipulation techniques and the measurement of the refractive index of particles.