The effect of non-Boltzmann energy distributions on the pressure, impingement rate, and heat flux of a monatomic gas in contact with a solid surface is investigated via theory and simulation. First, microcanonical formulations of the pressure, impingement rate, and heat flux are derived from first principles and integrated with prototypical energy distributions. Second, atomistic molecular dynamics simulations of an iron nanowire in a low-pressure argon atmosphere are used to test the non-Boltzmann heat flux theory. While pressure is found to be unaffected by the energy distribution of the gas, the impingement rate increases by up to 8.5% in the non-Boltzmann case. Most intriguingly, non-Boltzmann energy distributions can lead to a negative heat flux, meaning that heat flows from the cold solid to the hot gas. This non-Boltzmann heat flux effect is validated via the molecular dynamics simulations, and the solid is found to be 46% colder than the gas in the case of a hypothetical equilibrium for the upper-limiting non-Boltzmann energy distributions. The present fundamental findings provide novel insights into the properties of non-Boltzmann gases and improve the understanding of nonequilibrium dynamics.