Dyonic Taub-NUT-AdS spaces: Phase structures of all horizon geometries

We study phase structures of Lorentzian Dyonic Taub-NUT-AdS spacetimes for different horizon geometries, which could be spherical, flat, or hyperbolic. We check the consistency of our extended thermodynamics approach through satisfying the first law, the Gibbs-Duhem relation, and the generalized Smarr's relation. Although we study the phase structure for the three cases, we give special attention to the flat and hyperbolic cases since they are known to show no phase transitions and were not studied before. Working in a mixed ensemble, we found that the behaviors of the flat and hyperbolic cases are different from those of a charged black hole. In the latter case, a continuous phase transition occurs at high temperatures and pressures, i.e., above the critical point, but in our cases it occurs at low temperatures and pressures, i.e., below the critical point! Generically, the spherical case is characterized by two critical points with continuous phase transition between them.