A neutron star is a highly dense object which lasts after a supernova explosion. The density of a neutron star overcomes the nuclear density, and the temperature is high at the beginning of its history. An isolated neutron star does not have any heat sources, and it cools down emitting thermal energy by neutrinos. The neutrino emission process depends on the state of interior matter of the neutron star. To compare theoretical simulations and observations of neutron stars, it can constrain the nuclear theory of high density region. We create a model of neutron stars with colour superconducting quark matter and nucleon superfluidity / superconductivity, to satisfy recent observations, including two 2M ⊙ neutron stars. We parameterize these super-states and demonstrate the cooling curves, which show heavy stars do not always cool faster than lighter stars.