We present a Bayesian analysis to constrain the equation of state of dense nucleonic matter by exploiting the available data from symmetric nuclear matter at saturation and from observations of compact X-ray sources and from the gravitational wave event GW170817. For the first time, such analysis is performed by using a class of models, the relativistic mean field models, which allow to consistently construct an equation of state in a wide range of densities, isospin asymmetries and temperatures. The selected class of models contains five nuclear physics empirical parameters at saturation for which we construct the joint posterior distributions. By exploring different types of priors, we find that the equations of state with the largest evidence are the ones featuring a strong reduction of the effective mass of the nucleons in dense matter which can be interpreted as an indication of a phase transition to a chiral symmetry restored phase. Those equations of state in turn predict $R_{1.4} \sim 12$ km. Finally, we present a preliminary investigation on the effect of including $\Lambda$ hyperons showing that they appear in stars more massive than about $1.6 M_{\odot}$ and lead to radii larger than about $R_{1.4} \sim 14$ km. Within the model here explored, the formation of such particles provide a poor agreement with the constraints from GW170817.

We study the partial fragmentation of a strange quark star into strangelets during the process of the merger of two strange quark stars or of a strange quark star with a neutron star. We discuss the fate of the fragments considering their possible evaporation into nucleons. We show that only the small amount of large size strangelets produced by tidal forces can survive evaporation, while the strangelets formed during and after the merger evaporate completely into nucleons. In this way we show that: 1) the density of strangelets in the galaxy is too low to be able to affect significantly the evolution of stars and the probability of direct detection of a strangelet is negligible; 2) the matter ejected by the strange quark star during and after the merger evaporates into nucleons and therefore can contribute to the kilonova signal.

We analyse the phenomenological implications of the two-families scenario on the merger of compact stars. That scenario is based on the coexistence of both hadronic stars and strange quark stars. After discussing the classification of the possible mergers, we turn to detailed numerical simulations of the merger of two hadronic stars, i.e., "first family" stars in which delta resonances and hyperons are present, and we show results for the threshold mass of such binaries, for the mass dynamically ejected and the mass of the disk surrounding the post-merger object. We compare these results with those obtained within the one-family scenario and we conclude that relevant signatures of the two-families scenario can be suggested, in particular: the possibility of a rapid collapse to a black hole for masses even smaller than the ones associated to GW170817; during the first milliseconds, oscillations of the postmerger remnant at frequencies higher than the ones obtained in the one-family scenario; a large value of the mass dynamically ejected and a small mass of the disk, for binaries of low total mass. Finally, based on a population synthesis analysis, we present estimates of the number of mergers for: two hadronic stars; hadronic star - strange quark star; two strange quark stars. We show that for unequal mass systems and intermediate values of the total mass, the merger of a hadronic star and a strange quark star is very likely (GW170817 has a possible interpretation into this category of mergers). On the other hand, mergers of two strange quark stars are strongly suppressed.

We discuss the different signals, in gravitational and electromagnetic waves, emitted during the merger of two compact stars. We will focus in particular on the possible contraints that those signals can provide on the equation of state of dense matter. Indeed, the stiffness of the equation of state and the particle composition of the merging compact stars, strongly affect e.g. the life time of the post-merger remnant and its gravitational wave signal, the emission of the short gamma-ray-burst, the amount of ejected mass and the related kilonova. The first detection of gravitational waves from the merger of two compact stars in August 2017, GW170817, and the subsequent detections of its electromagnetic counterparts, GRB170817A and AT2017gfo, is the first example of the era of "multi-messenger astronomy": we discuss what we have learned from this detection on the equation of state of compact stars and we provide a tentative interpretation of this event, within the two families scenario, as due to the merger of a hadronic star with a quark star.

The merger of two neutron stars is a very complex process. In order to disentangle the various steps through which it takes place it is mandatory to examine all the signals we can detect: gravitational waves and electromagnetic waves, in a huge spectrum ranging from X and $\gamma$-rays down to infrared and radio. Each of these signals provides a message and the totality of this information will allow us not only to understand the process of the merger but also the behavior of matter at those extreme conditions.