The goal of this work is to investigate the effects of chemical composition of heat blanketing envelopes of neutron stars on their thermal states and thermal evolution. To this purpose, we employ newly constructed models of the envelopes composed of binary ion mixtures (H–He, He–C, C–Fe) varying the mass of lighter ions (H, He or C) in the envelope. The results are compared with those calculated using the standard “onion-like” envelope. For illustration, we apply these results to estimate the internal temperature of the Vela pulsar and to study cooling of neutron stars. We show that uncertainties in the chemical composition of the envelopes can lead up to ~ 2.5 times uncertainty of the internal temperature of the star which significantly complicates theoretical reconstruction of the internal structure of cooling neutron stars from observations of their thermal surface emission.