Next Article in Journal / Special Issue
Two Types of Jets and Quark and Chromon Model in QCD
Previous Article in Journal / Special Issue
Topological Entanglement and Knots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Bayesian Analysis for Extracting Properties of the Nuclear Equation of State from Observational Data Including Tidal Deformability from GW170817 †

by
Alexander Ayriyan
1,
David Alvarez-Castillo
2,
David Blaschke
2,3,4,* and
Hovik Grigorian
1,5
1
Laboratory for Information Technologies, Joint Institute for Nuclear Research, Joliot-Curie Street 6, Dubna 141980, Russia
2
Bogoliubov Laboratory for Theoretical Physics, Joint Institute for Nuclear Research, Joliot-Curie Street 6, Dubna 141980, Russia
3
Institute of Theoretical Physics, University of Wroclaw, Max Born Place 9, 50-204 Wroclaw, Poland
4
National Research Nuclear University (MEPhI), Kashirskoe Shosse 31, Moscow 115409, Russia
5
Department of Physics, Yerevan State University, Alek Manukyan Str. 1, Yerevan 0025, Armenia
*
Author to whom correspondence should be addressed.
This paper is based on the talk at the 7th International Conference on New Frontiers in Physics (ICNFP 2018), Crete, Greece, 4–12 July 2018.
Submission received: 27 December 2018 / Revised: 3 February 2019 / Accepted: 5 February 2019 / Published: 14 February 2019

Abstract

:
We develop a Bayesian analysis method for selecting the most probable equation of state under a set of constraints from compact star physics, which now include the tidal deformability from GW170817. We apply this method for the first time to a two-parameter family of hybrid equations of state that is based on realistic models for the hadronic phase (KVORcut02) and the quark matter phase (SFM α ) which produce a third family of hybrid stars in the mass–radius diagram. One parameter ( α ) characterizes the screening of the string tension in the string-flip model of quark matter while the other ( Δ P ) belongs to the mixed phase construction that mimics the thermodynamics of pasta phases and includes the Maxwell construction as a limiting case for Δ P = 0 . We present the corresponding results for compact star properties like mass, radius and tidal deformabilities and use empirical data for them in the newly developed Bayesian analysis method to obtain the probabilities for the model parameters within their considered range.

Graphical Abstract

1. Introduction

The determination of the equation of state (EoS) of compact star (CS) interiors is a subject of active research by means of both laboratory measurements and astronomical observations. One of the most urgent questions concerns the possible existence of deconfined quark matter in CS cores, where matter densities can exceed by several times the nuclear saturation value, n 0 = 0.15 fm 3 , the typical density in large atomic nuclei. Hybrid compact stars have an inner core composed of quark matter surrounded by an outer core of nuclear matter. Of particular interest is the CS mass twin (MT) phenomenon [1]: when a pair of stars has the same gravitational mass but different radii, with the larger star being composed only of hadronic matter and the smaller one being a hybrid star with a quark matter core. The presence of MTs is synonymous with the existence of a third family [2] of CS in the mass–radius (M-R) diagram. It requires the CS EoS to feature a strong first order phase transition from hadron to quark matter [3,4,5,6,7,8], where the latent heat of the transition fulfills the Seidov criterion [9].
Should it turn out that the mass twin phenomenon can be discovered in CS observations, this would give indirect evidence for the existence of a critical endpoint in the phase diagram of quantum chromodynamics (QCD) [4,5,10] which is sought for in relativistic heavy-ion collisions, so far without being conclusive. On the other hand, astronomical observations may provide measurements of masses and/or radii which are relevant for constraining the nuclear EoS via the one-to-one relationship provided by the Tolman–Oppenheimer–Volkoff (TOV) equations [11,12] which govern the general relativistic hydrodynamic stability of CS matter under its own gravity. Neutron star masses are accurately measured by monitoring binary pulsar dynamics, whereas the radii were determined so far with big uncertainties. The situation has changed with the first observation of the gravitational wave signal from the inspiral phase of the binary CS merger GW170817 [13], which allowed constraining the tidal deformability of both merging stars and thus their mass and radius [14,15,16,17].
Consequently, theoretical approaches are required to quantitatively assess the most probable EoS out of a whole class obtained by varying intrinsic model parameters. One of the most powerful methods is the Bayesian analysis (BA) or Bayesian interpretation of probability that allows for estimation of model parameters based on prior knowledge, in this case the already collected measurements. Several works in this direction include BA based on observations of X-ray bursters [18] that unfortunately suffer from uncertainties in the stellar atmosphere composition modeling required in the interpretation of the observed signal or the more general approach of [19] that includes a generic multipolytrope EoS with priors that include experimental nuclear matter pressure values and that is able to support the most massive CS of about 2 M and reports an accuracy of about 30 % in pressure values. The recent approach of [20] employs X-ray timing observations of accretion-powered millisecond pulsars, resulting in radius estimates of about 5 % if the CS mass is known. Moreover, a more stringent analysis [21] that incorporates the hypothesis of the Direct Urca cooling constraint [22] in addition to the aforementioned measurements concludes that the neutron star radius has a value of 12.7 ± 0.4 km for masses ranging from 1 up to 2 M .
Despite the fact that the recent detection of gravitational radiation from the inspiral phase of the binary CS merger GW170817 has led to constraints on the tidal deformability of CS matter and to the discussion of the possibility that one or even both of the CS in the binary could be hybrid stars with quark matter interior and there would be a possibility to discover the mass twin phenomenon and thus a strong first order phase transition in CS matter [23,24,25,26,27], these data have not yet been included in BA studies such as the ones already listed above. Thus, it is of great interest to update such BA studies by incorporating this new information to constrain the CS EoS.
In this work, we consider the mass twin EoS and modifications to it in order to perform new Bayesian analysis calculations that, in addition to the previously used CS measurements in [28], will include the GW170817 data as priors to assess the probability of model parameters. Our model choice is the KVOR EoS for hadronic matter [29] together with the String–Flip approach [30] in its relativistic density functional formulation [31] for the deconfined quark matter regime.

2. Hybrid EoS

The hybrid EoS has been produced with the one-parameter replacement interpolation method (for a second order polynomial ansatz for P ( μ ) ) [32,33,34] which mimics the thermodynamic behavior of pasta phases in the coexistence region between the pure hadronic phase and the pure quark one.
The hadronic phase in this work is described by the well known KVOR equation of state [29] with a modification of stiffness as introduced in [35] and denoted as KVORcut2. This particular version of the KVOR EoS model is the stiffest one presented in that work and allows for fulfilling the condition [9] for the latent heat of phase transition to quark matter, so that the disconnected hybrid star family at higher densities is possible.
The quark phase is based on the String–Flip Model with the so-called density functional [31] including the available volume fraction Φ ,
Φ ( n B ) = e α n B 2 fm 6 ,
where the parameter α describes the effective density-dependence of the confining interaction, and it is varied here in the range α 0.1 , , 0.3 . The value of this parameter is correlated with the critical density of the phase transition, as it is shown in the figures below. The larger α , the lower the critical density and therefore the critical mass for the onset of the phase transition in the hybrid star. These systematics allows for calibrating the range of the model parameters with the observational data without contradicting the known properties of nuclear matter at saturation density.
For the construction of the hybrid star EoS from the hadronic and quark matter EoS, we employ here the mixed phase construction that is described in detail in Refs. [32,33,34]. It introduces the free parameter Δ P = P ( μ c ) / P c 1 , the pressure increment at the critical chemical potential μ c relative to the critical pressure P c of the corresponding Maxwell construction, which is then obtained as a limiting case for Δ P = 0 . The resulting family of hybrid EoS corresponds to the one described in [32] (see Figure 1). This way of modelling the phase transition mimics the possibility of so-called pasta phases characterized by different geometric structures in the coexistence region of the transition sufficiently well, as has been shown in [36]. In that reference, it has also been demonstrated that the parameter Δ P of the transition construction can be related to the value of the surface tension at the hadron to quark matter interface.

3. Neutron Star Configurations

3.1. Mass and Radius

The structure and global properties of compact stars are obtained by solving the TOV equations,
d P ( r ) d r = G M ( r ) ε ( r ) r 2 1 + P ( r ) ε ( r ) 1 + 4 π r 3 P ( r ) M ( r ) 1 2 G M ( r ) r ,
d M ( r ) d r = 4 π r 2 ε ( r ) ,
where P ( r ) , ε ( r ) and M ( r ) are the profiles of pressure, energy density and enclosed mass as a function of the distance r from the center of the star. The radius R of the star is obtained from the condition of vanishing pressure at the surface P ( R ) = 0 and the gravitational mass of the star is M = M ( R ) .

3.2. Tidal Deformability

The recent observation of gravitational waves from the inspiral phase of the binary CS merger GW170817 gives an estimation of the tidal deformabilities Λ 1 and Λ 2 of these stars; therefore, we also include this information in the analyses. The tidal deformability Λ of a star can be calculated from the Einstein equation for small elliptic deformation as described in [37]. For the results of the calculation, see Figure 2.

4. Bayesian Inference for the EoS Models

4.1. Vector of Parameters

The set of parameters of models could be represented in the parameter space with introduction of the vector of parameters, each vector is one fixed model from the types of hadronic described above, quark phases and transition construction:
π i = α ( j ) , Δ P ( k ) ,
where i = 0 , , N 1 and i = N 2 × j + k and j = 0 , , N 1 1 , k = 0 , , N 2 1 and N 1 and N 2 are number of values of model parameters α and Δ P correspondingly.

4.2. Likelihood of a Model for the Λ 1 Λ 2 Constraint from GW170817

In order to implement the tidal deformability constraint on the compact stars EoS, reflected on the Λ 1 Λ 2 diagram that includes probability regions from GW170817 event [13,38], we employ the definition of the probability as an integral over the probability distribution function (PDF) β ( Λ 1 , Λ 2 ) by
P E G W π i = l 22 β ( Λ 1 ( τ ) , Λ 2 ( τ ) ) d τ ,
when both stars in the binary merger belong to the second family branch of neutron stars (and l 22 is the corresponding path in the Λ 1 Λ 2 plane), or by
P E G W π i = l 22 β ( Λ 1 ( τ ) , Λ 2 ( τ ) ) d τ + l 23 β ( Λ 1 ( τ ) , Λ 2 ( τ ) ) d τ .
In case the parameter vector π i corresponds to a hybrid star equation of state with a third family of compact stars in the mass range relevant for the merger, then l 23 is the path in the Λ 1 Λ 2 plane corresponding to a binary merger composed of a neutron star and a hybrid star. The parameter τ defines the position of a point on the paths l 22 and l 23 in Equations (5) and (6). Note that the PDF has been reconstructed by the method Gaussian kernel density estimation with Λ 1 Λ 2 data given at LIGO web-page [39] (see Figure 3).

4.3. Likelihood of a Model for the Mass Constraint

The likelihood of the model is the conditional probability of the expected value of the possible maximum mass for the given model parameter vector:
P E A π i = Φ ( M i , μ A , σ A ) ,
where M i is maximum mass of the given by π i , and μ A = 2.01 M and σ A = 0.04 M is the mass measurement [40].

4.4. Posterior Distribution

The full likelihood for the given π i can be calculated as a product of all likelihoods, since the considered constraints are independent from each other
P E π i = m P E m π i .
The posterior distribution of models on the parameter diagram is given by Bayes’ theorem
P π i E = P E π i P π i j = 0 N 1 P E π j P π j ,
where P π j is a prior distribution of a model taken to be uniform: P π j = 1 / N .

5. Results

The results of the Bayesian analysis are given in Figure 4. The most likely EoS are those with a strong mixed phase effect described by Δ P > 0.04 and with a large screening parameter α > 0.28 , at the limit of the parameter range considered here. For these parameter sets, the phase transition and therefore the compactification of the hybrid star configuration occurs within the mass range that is relevant for GW170817. These results indicate that it should be worthwile to repeat this exploratory calculation with a wider set of parameters, in particular a larger screening parameter α .

6. Conclusions

We have developed a Bayesian analysis method for selecting the most probable equation of state under a set of constraints from compact star physics, which now include the tidal deformability from GW170817. We have applied this method to a case study that employed the class of hybrid EoS introduced in [32] that allows for the existence of a third family of compact stars. We have investigated the requirements on future measurements to find the equivalent phenomenon of mass twins, i.e., two objects with the same mass but different, distinguishable radii. A mass–radius relation with two branches is mapped also on the Λ 1 Λ 2 diagram as it is shown here in Figure 3 and other recent publications [23,24,25,26,27]. Since binary compact star mergers are expected to be observed by the LIGO-Virgo Collaboration at a rate of 1–10 events per year, we expect that, with the observation of next merger events, a binodal structure of the PDF in the Λ 1 Λ 2 plane could become apparent as a manifestation of the low-mass twin case with an onset of the third family branch below 1.3 M if such a branch exists in nature. A similar suggestion has been already proposed by Christian et al. [25]. Such an observation would support the existence of an early phase transition, around 2 n 0 for strong in-medium screening of the string tension.

Author Contributions

Conceptualization and methodology, A.A. and H.G.; software, A.A. and D.A.-C.; investigation, A.A.; data curation, A.A. and D.A.-C.; writing—original draft preparation, A.A. and D.A.-C.; writing—review and editing, D.B. and H.G.; visualization, A.A.; supervision, H.G.; project administration, D.B.; funding acquisition, D.B. and D.A.-C.

Funding

A.A., D.B. and H.G. acknowledge support from the Russian Science Foundation under Grant No. 17-12-01427 for the work described in Section 2, Section 3, Section 4, Section 5 and Section 6. D.A.-C. is grateful for partial support from the Ter–Antonian–Smorodinsky program for collaboration between JINR and Armenian scientific institutions and from the Bogoliubov–Infeld program for collaboration between JINR and Polish Institutions.

Acknowledgments

We acknowledge discussions with Konstantin Maslov on the hybrid star EoS.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Glendenning, N.K.; Kettner, C. Nonidentical neutron star twins. Astron. Astrophys. 2000, 353, L9. [Google Scholar]
  2. Gerlach, U.H. Equation of State at Supranuclear Densities and the Existence of a Third Family of Superdense Stars. Phys. Rev. 1968, 172, 1325–1330. [Google Scholar] [CrossRef]
  3. Alford, M.G.; Han, S.; Prakash, M. Generic conditions for stable hybrid stars. Phys. Rev. D 2013, 88, 083013. [Google Scholar] [CrossRef]
  4. Alvarez-Castillo, D.E.; Blaschke, D. Proving the CEP with compact stars? In Proceedings of the 17th Conference of Young Scientists and Specialists, Dubna, Russia, 8–12 April 2013; pp. 22–26. Available online: http://xxx.lanl.gov/abs/1304.7758 (accessed on 30 April 2013).
  5. Benic, S.; Blaschke, D.; Alvarez-Castillo, D.E.; Fischer, T.; Typel, S. A new quark-hadron hybrid equation of state for astrophysics-I. High-mass twin compact stars. Astron. Astrophys. 2015, 577, A40. [Google Scholar] [CrossRef]
  6. Blaschke, D.; Alvarez-Castillo, D.E. High-mass twins & resolution of the reconfinement, masquerade and hyperon puzzles of compact star interiors. AIP Conf. Proc. 2016, 1701, 020013. [Google Scholar] [Green Version]
  7. Alvarez-Castillo, D.E.; Blaschke, D.B. High-mass twin stars with a multipolytrope equation of state. Phys. Rev. C 2017, 96, 045809. [Google Scholar] [CrossRef]
  8. Paschalidis, V.; Yagi, K.; Alvarez-Castillo, D.; Blaschke, D.B.; Sedrakian, A. Implications from GW170817 and I-Love-Q relations for relativistic hybrid stars. Phys. Rev. D 2018, 97, 084038. [Google Scholar] [CrossRef] [Green Version]
  9. Seidov, Z.F. The Stability of a Star with a Phase Change in General Relativity Theory. Soviet Astronomy 1971, 15, 347. [Google Scholar]
  10. Blaschke, D.; Alvarez-Castillo, D.E.; Benic, S. Mass-radius constraints for compact stars and a critical endpoint. In Proceedings of the 8th International Workshop on Critical Point and Onset of De confinement, Napa, CA, USA, 11–15 March 2013. [Google Scholar]
  11. Tolman, R.C. Static solutions of Einstein’s field equations for spheres of fluid. Phys. Rev. 1939, 55, 364–373. [Google Scholar] [CrossRef]
  12. Oppenheimer, J.R.; Volkoff, G.M. On Massive neutron cores. Phys. Rev. 1939, 55, 374–381. [Google Scholar] [CrossRef]
  13. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; LIGO Scientific Collaboration and Virgo Collaboration; et al. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Annala, E.; Gorda, T.; Kurkela, A.; Vuorinen, A. Gravitational-wave constraints on the neutron-star-matter Equation of State. Phys. Rev. Lett. 2018, 120, 172703. [Google Scholar] [CrossRef] [PubMed]
  15. Bauswein, A.; Just, O.; Janka, H.T.; Stergioulas, N. Neutron-star radius constraints from GW170817 and future detections. Astrophys. J. 2017, 850, L34. [Google Scholar] [CrossRef]
  16. Rezzolla, L.; Most, E.R.; Weih, L.R. Using gravitational-wave observations and quasi-universal relations to constrain the maximum mass of neutron stars. Astrophys. J. 2018, 852. [Google Scholar] [CrossRef]
  17. De, S.; Finstad, D.; Lattimer, J.M.; Brown, D.A.; Berger, E.; Biwer, C.M. Tidal Deformabilities and Radii of Neutron Stars from the Observation of GW170817. Phys. Rev. Lett. 2018, 121, 091102. [Google Scholar] [CrossRef] [PubMed]
  18. Steiner, A.W.; Lattimer, J.M.; Brown, E.F. The Equation of State from Observed Masses and Radii of Neutron Stars. Astrophys. J. 2010, 722, 33–54. [Google Scholar] [CrossRef]
  19. Raithel, C.A.; Özel, F.; Psaltis, D. From Neutron Star Observables to the Equation of State. II. Bayesian Inference of Equation of State Pressures. Astrophys. J. 2017, 844, 156. [Google Scholar] [CrossRef] [Green Version]
  20. Salmi, T.; Nättilä, J.; Poutanen, J. Bayesian parameter constraints for neutron star masses and radii using X-ray timing observations of accretion-powered millisecond pulsars. Astron. Astrophys. 2018, 618, A161. [Google Scholar] [CrossRef]
  21. Margueron, J.; Hoffmann Casali, R.; Gulminelli, F. Equation of state for dense nucleonic matter from metamodeling. II. Predictions for neutron star properties. Phys. Rev. C 2018, 97, 025806. [Google Scholar] [CrossRef] [Green Version]
  22. Blaschke, D.; Alvarez-Castillo, D.E.; Klähn, T. Universal Symmetry Energy Contribution to the Neutron Star Equation of State. 2016. Available online: http://xxx.lanl.gov/abs/1604.08575 (accessed on 29 April 2016).
  23. Alvarez-Castillo, D.E.; Blaschke, D.B.; Grunfeld, A.G.; Pagura, V.P. Third Family of Compact Stars within A Nonlocal Chiral Quark Model Equation of State. 2018. Available online: http://xxx.lanl.gov/abs/1805.04105v3 (accessed on 1 January 2019).
  24. Most, E.R.; Weih, L.R.; Rezzolla, L.; Schaffner-Bielich, J. New constraints on radii and tidal deformabilities of neutron stars from GW170817. Phys. Rev. Letters 2018, 120, 261103. [Google Scholar] [CrossRef]
  25. Christian, J.E.; Zacchi, A.; Schaffner-Bielich, J. Signals in the tidal deformability for phase transitions in compact stars with constraints from GW170817. Phys. Rev. D 2019, 99, 023009. [Google Scholar] [CrossRef]
  26. Montana, G.; Tolos, L.; Hanauske, M.; Rezzolla, L. Constraining Twin Stars with GW170817. 2018. Available online: http://xxx.lanl.gov/abs/1811.10929 (accessed on 28 November 2018).
  27. Sieniawska, M.; Turczanski, W.; Bejger, M.; Zdunik, J.L. Tidal Deformability and Other Global Parameters of Compact Stars with Phase Transitions. 2018. Available online: http://xxx.lanl.gov/abs/1807.11581 (accessed on 27 July 2018).
  28. Alvarez-Castillo, D.; Ayriyan, A.; Benic, S.; Blaschke, D.; Grigorian, H.; Typel, S. New class of hybrid EoS and Bayesian M-R data analysis. Eur. Phys. J. A 2016, 52, 69. [Google Scholar] [CrossRef]
  29. Kolomeitsev, E.E.; Voskresensky, D.N. Relativistic mean-field models with effective hadron masses and coupling constants, and rho- condensation. Nucl. Phys. A 2005, 759, 373–413. [Google Scholar] [CrossRef]
  30. Röpke, G.; Blaschke, D.; Schulz, H. Pauli Quenching Effects in a Simple String Model of Quark / Nuclear Matter. Phys. Rev. D 1986, 34, 3499. [Google Scholar] [CrossRef]
  31. Kaltenborn, M.A.R.; Bastian, N.U.F.; Blaschke, D.B. Quark-nuclear hybrid star equation of state with excluded volume effects. Phys. Rev. D 2017, 96, 056024. [Google Scholar] [CrossRef] [Green Version]
  32. Ayriyan, A.; Bastian, N.U.; Blaschke, D.; Grigorian, H.; Maslov, K.; Voskresensky, D.N. Robustness of third family solutions for hybrid stars against mixed phase effects. Phys. Rev. C 2018, 97, 045802. [Google Scholar] [CrossRef] [Green Version]
  33. Ayriyan, A.; Grigorian, H. Model of the Phase Transition Mimicking the Pasta Phase in Cold and Dense Quark-Hadron Matter. EPJ Web Conf. 2018, 173, 03003. [Google Scholar] [CrossRef] [Green Version]
  34. Abgaryan, V.; Alvarez-Castillo, D.; Ayriyan, A.; Blaschke, D.; Grigorian, H. Two Novel Approaches to the Hadron-Quark Mixed Phase in Compact Stars. Universe 2018, 4, 94. [Google Scholar] [CrossRef]
  35. Maslov, K.A.; Kolomeitsev, E.E.; Voskresensky, D.N. Relativistic Mean-Field Models with Scaled Hadron Masses and Couplings: Hyperons and Maximum Neutron Star Mass. Nucl. Phys. A 2016, 950, 64–109. [Google Scholar] [CrossRef]
  36. Maslov, K.; Yasutake, N.; Ayriyan, A.; Blaschke, D.; Grigorian, H.; Maruyama, T.; Tatsumi, T.; Voskresensky, D.N. Hybrid Equation of State with Pasta Phases and Third Family of Compact Stars I: Pasta Phases and Effective Mixed Phase Model. 2018. Available online: http://xxx.lanl.gov/abs/1812.11889 (accessed on 31 December 2018).
  37. Hinderer, T.; Lackey, B.D.; Lang, R.N.; Read, J.S. Tidal deformability of neutron stars with realistic equations of state and their gravitational wave signatures in binary inspiral. Phys. Rev. D 2010, 81, 123016. [Google Scholar] [CrossRef]
  38. Abbott, B.P.; LIGO Scientific Collaboration and Virgo Collaboration. GW170817: Measurements of neutron star radii and equation of state. Phys. Rev. Lett. 2018, 121, 161101. [Google Scholar] [CrossRef] [PubMed]
  39. LIGO Scientific Collaboration and Virgo Collaboration. GW170817: Measurements of Neutron Star Radii and Equation of State. 2018. Available online: https://dcc.ligo.org/LIGO-P1800115/public (accessed on 30 January 2019).
  40. Antoniadis, J.; Freire, P.C.C.; Wex, N.; Tauris, T.M.; Lynch, R.S.; van Kerkwijk, M.H.; Kramer, M.; Bassa, C.; Dhillon, V.S.; et al. A Massive Pulsar in a Compact Relativistic Binary. Science 2013, 340, 1233232. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The set of hybrid EoS obtained by the mixed phase construction for different values of Δ P and α is shown in the left panel, while in the right panel the corresponding set of compact star sequences in the mass–radius diagram is shown. The blue numbers show selected values of critical densities for the onset of the phase transitions obtained by a Maxwell construction.
Figure 1. The set of hybrid EoS obtained by the mixed phase construction for different values of Δ P and α is shown in the left panel, while in the right panel the corresponding set of compact star sequences in the mass–radius diagram is shown. The blue numbers show selected values of critical densities for the onset of the phase transitions obtained by a Maxwell construction.
Universe 05 00061 g001
Figure 2. The dependence of the tidal deformability Λ on the compact star mass (left panel) and the Λ 1 Λ 2 diagram for selected values of Δ P and α (right panel). Comparison with the 90% confidence line of the LIGO-Virgo Collaboration for the low-spin prior of GW170817 [13] shows that, if the hadronic EoS is as stiff as DD2p40, at least one of the stars in the binary has to be a hybrid star in order to avoid a violation of the Λ 1 Λ 2 constraint. As can be seen from figure 17 of Ref. [23], the merger of two hybrid stars which is admissible when the onset mass for the deconfinement transition is lowered, e.g., by increasing α , would lead to the appearance of a new branch in this diagram mimicking the pattern of a merger of two neutron stars with a rather soft nuclear EoS.
Figure 2. The dependence of the tidal deformability Λ on the compact star mass (left panel) and the Λ 1 Λ 2 diagram for selected values of Δ P and α (right panel). Comparison with the 90% confidence line of the LIGO-Virgo Collaboration for the low-spin prior of GW170817 [13] shows that, if the hadronic EoS is as stiff as DD2p40, at least one of the stars in the binary has to be a hybrid star in order to avoid a violation of the Λ 1 Λ 2 constraint. As can be seen from figure 17 of Ref. [23], the merger of two hybrid stars which is admissible when the onset mass for the deconfinement transition is lowered, e.g., by increasing α , would lead to the appearance of a new branch in this diagram mimicking the pattern of a merger of two neutron stars with a rather soft nuclear EoS.
Universe 05 00061 g002
Figure 3. The PDF reconstructed with Λ 1 Λ 2 data for GW170817 from the LIGO website [39] as a 3D graphics (left panel) and as a contour plot (right panel).
Figure 3. The PDF reconstructed with Λ 1 Λ 2 data for GW170817 from the LIGO website [39] as a 3D graphics (left panel) and as a contour plot (right panel).
Universe 05 00061 g003
Figure 4. Left panel: The posterior distribution of models on the parameter space spanned by Δ P and α . Right panel: The compact star sequences in the mass–radius diagram labeled into four probability classes according to the results of the BA for the posterior distribution of the left panel. The grey, brown, orange and black lines show sequences for which the value of the posterior probability exceeds the thresholds of 0.0, 0.02, 0.04 and 0.06, respectively. Due to the restriction to the model and parameter range used in Ref. [32], the sequences with a lower onset mass for the deconfinement transition are not accessed for which an interpretation of GW170817 as a merger of two hybrid stars would be possible.
Figure 4. Left panel: The posterior distribution of models on the parameter space spanned by Δ P and α . Right panel: The compact star sequences in the mass–radius diagram labeled into four probability classes according to the results of the BA for the posterior distribution of the left panel. The grey, brown, orange and black lines show sequences for which the value of the posterior probability exceeds the thresholds of 0.0, 0.02, 0.04 and 0.06, respectively. Due to the restriction to the model and parameter range used in Ref. [32], the sequences with a lower onset mass for the deconfinement transition are not accessed for which an interpretation of GW170817 as a merger of two hybrid stars would be possible.
Universe 05 00061 g004

Share and Cite

MDPI and ACS Style

Ayriyan, A.; Alvarez-Castillo, D.; Blaschke, D.; Grigorian, H. Bayesian Analysis for Extracting Properties of the Nuclear Equation of State from Observational Data Including Tidal Deformability from GW170817. Universe 2019, 5, 61. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5020061

AMA Style

Ayriyan A, Alvarez-Castillo D, Blaschke D, Grigorian H. Bayesian Analysis for Extracting Properties of the Nuclear Equation of State from Observational Data Including Tidal Deformability from GW170817. Universe. 2019; 5(2):61. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5020061

Chicago/Turabian Style

Ayriyan, Alexander, David Alvarez-Castillo, David Blaschke, and Hovik Grigorian. 2019. "Bayesian Analysis for Extracting Properties of the Nuclear Equation of State from Observational Data Including Tidal Deformability from GW170817" Universe 5, no. 2: 61. https://0-doi-org.brum.beds.ac.uk/10.3390/universe5020061

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop