The properties of compact stars and in particular the existence of twin star solutions are investigated within an effective model that is constrained by lattice QCD thermodynamics. The model is modified at large baryon densities to incorporate a large variety of scenarios of first order phase transitions to a phase of deconfined quarks. This is achieved by matching two different variants of the bag model equation of state, in order to estimate the role of the Bag model parameters on the appearance of a second family of neutron stars. The produced sequences of neutron stars are compared with modern constrains on stellar masses, radii, and tidal deformability from astrophysical observations and gravitational wave analyses. It is found that those scenarios in our analysis, in which a third family of stars appeared due to the deconfinement transition, are disfavored from astrophysical constraints.
In this work, we study for the first time the thermal evolution of twin star pairs, i.e., stars that present the same mass but different radius and compactness. We collect available equations of state that give origin to a second branch of stable compact stars with quarks in their core. For each equation of state, we investigate the particle composition inside stars and how differently each twin evolves over time, which depends on the central density/pressure and consequent crossing of the threshold for the Urca cooling process. We find that, although the general stellar thermal evolution depends on mass and particle composition, withing one equation of state, only twin pairs that differ considerably on compactness can be clearly distinguished by how they cool down.