What neutron stars tell about the hadron-quark phase transition: A Bayesian studyhdl:21.15109/CONCORDA/O5XS5AARP2023-11-061Takátsy, János; Kovács, Péter; Wolf, György; Schaffner-Bielich, Juergen, 2023, "What neutron stars tell about the hadron-quark phase transition: A Bayesian study", https://hdl.handle.net/21.15109/CONCORDA/O5XS5A, ARP, V1What neutron stars tell about the hadron-quark phase transition: A Bayesian studyhdl:21.15109/CONCORDA/O5XS5ATakátsy, JánosKovács, PéterWolf, GyörgySchaffner-Bielich, JuergenFK 131982K 138277ÚNKP-22-3315477589-TRR 211315477589-TRR 211ARPKovács, PéterKovács, Péter2023-10-18PhysicsAstronomy and AstrophysicsEffective Field TheoryNeutron stars & pulsarsStatistical methodsThe existence of quark matter inside the heaviest neutron stars has been the topic of numerous recent studies, many of them suggesting that a phase transition to strongly interacting conformal matter inside neutron stars is feasible. Here we examine this hybrid star scenario using a soft and a stiff hadronic model, a constituent quark model with three quark flavors, and applying a smooth crossover transition between the two. Within a Bayesian framework, we study the effect of up-to-date constraints from neutron star observations on the equation-of-state (EOS) parameters and various neutron star observables. Our results show that a pure quark core is only possible if the maximum mass of neutron stars is below ~2.35 M⊙. However, we also find, consistent with other studies, that a peak in the speed of sound, exceeding 1/3, is highly favored by astrophysical measurements, which might indicate the percolation of hadrons at ~3-4 n0. Even though our prediction for the phase transition parameters varies depending on the specific astrophysical constraints utilized, the position of the speed of sound peak only changes slightly, while the existence of pure quark matter below ~4 n0, using our parametrization, is disfavored. On the other hand, the preferred range for the EOS shows signs of conformality above ~4 n0. Additionally, we present the difference in the upper bounds of radius estimates using the full probability density data and sharp cutoffs and stress the necessity of using the former.János Takátsy, Péter Kovács, György Wolf, and Jürgen Schaffner-Bielich Phys. Rev. D 108, 04300210.1103/PhysRevD.108.043002János Takátsy, Péter Kovács, György Wolf, and Jürgen Schaffner-Bielich Phys. Rev. D 108, 0430020908.2344_DD2_Typel.pdfapplication/pdf0911.4073_DD2.pdfapplication/pdf1207.2184_SFHo.pdfapplication/pdf1302.4261_Togashi.pdfapplication/pdf1903.08963_Togashi_Baym.pdfapplication/pdfAPR.matapplication/matlab-matAPR_Compose.dattext/x-fixed-fieldAPR_Compose.zipapplication/zipAPR_EoS.dattext/x-fixed-fieldAPR_EoS.mtext/x-matlabBayes_Delta.figapplication/x-xfigBayes_Delta.pdfapplication/pdfBayes_analyzer.mtext/x-matlabBayes_cs2_MuMgap.figapplication/x-xfigBayes_cs2_MuMgap.pdfapplication/pdfBayes_cs2_prior.figapplication/x-xfigBayes_cs2_prior.pdfapplication/pdfBayes_dc.figapplication/x-xfigBayes_dc.pdfapplication/pdfBayes_dc_Mgap.figapplication/x-xfigBayes_dc_Mgap.pdfapplication/pdfBayes_par_NGWHESS_level.figapplication/x-xfigBayes_par_NGWHESS_level.pdfapplication/pdfBayes_par_NGWMgap_level.figapplication/x-xfigBayes_par_NGWMgap_level.pdfapplication/pdfBayes_par_NGWMu_level.figapplication/x-xfigBayes_par_NGWMu_level.pdfapplication/pdfBayes_par_NGW_level.figapplication/x-xfigBayes_par_NGW_level.pdfapplication/pdfBayes_par_NICER_level.figapplication/x-xfigBayes_par_NICER_level.pdfapplication/pdfBayes_par_prior_level.figapplication/x-xfigBayes_par_prior_level.pdfapplication/pdfBayes_prob_DD2.matapplication/matlab-matBayes_prob_M.mtext/x-matlabBayes_prob_SFHo.matapplication/matlab-matBayes_pt_31_NGW.figapplication/x-xfigBayes_pt_31_NGW.pdfapplication/pdfBayes_pt_31_NGWMu.figapplication/x-xfigBayes_pt_31_NGWMu.pdfapplication/pdfBayes_pt_31_NICER.figapplication/x-xfigBayes_pt_31_NICER.pdfapplication/pdfBayes_pt_31_prior.figapplication/x-xfigBayes_pt_31_prior.pdfapplication/pdfBayes_pt_47_NGW.figapplication/x-xfigBayes_pt_47_NGW.pdfapplication/pdfBayes_pt_47_NGWMgap.figapplication/x-xfigBayes_pt_47_NGWMgap.pdfapplication/pdfBayes_pt_47_NICER.figapplication/x-xfigBayes_pt_47_NICER.pdfapplication/pdfBayes_pt_47_prior.figapplication/x-xfigBayes_pt_47_prior.pdfapplication/pdfBayes_runner.mtext/x-matlabBetaEq.mtext/x-matlabBetaEq.mtext/x-matlabBetaEq.mtext/x-matlabBetaEq.mtext/x-matlabChiralPTDens.matapplication/matlab-matCorner_NGW.figapplication/x-xfigCorner_NGW.pdfapplication/pdfCorner_NGWMgap.figapplication/x-xfigCorner_NGWMgap.pdfapplication/pdfCorner_NGWMu.figapplication/x-xfigCorner_NGWMu.pdfapplication/pdfCorner_SFHo_290_priorLam.figapplication/x-xfigCorner_SFHo_290_priorLam.pdfapplication/pdfCorner_prior.figapplication/x-xfigCorner_prior.pdfapplication/pdfDD2_EoS.dattext/x-fixed-fieldDD2_weights.matapplication/matlab-matELSM_vector_fluct_2020_03.pdfapplication/pdfEoS_Bayesgen.mtext/x-matlabEoScontour.mtext/x-matlabEos_eps_pol_interp.mtext/x-matlabFSU2R_EoS.dattext/x-fixed-fieldFSU2R_EoS.mtext/x-matlabField_EQ_asymp_20210215_jav.pdfapplication/pdfGW170817_woEoS_qLam.matapplication/matlab-matHESS.matapplication/matlab-matMR_Bayes_NGWHESS_SFHo.figapplication/x-xfigMR_Bayes_NGWHESS_SFHo.pdfapplication/pdfMR_Bayes_NGWHESS_comb.figapplication/x-xfigMR_Bayes_NGWHESS_comb.pdfapplication/pdfMR_Bayes_NGWMgap_SFHo.figapplication/x-xfigMR_Bayes_NGWMgap_SFHo.pdfapplication/pdfMR_Bayes_NGWMgap_comb.figapplication/x-xfigMR_Bayes_NGWMgap_comb.pdfapplication/pdfMR_Bayes_NGWMu_SFHo.figapplication/x-xfigMR_Bayes_NGWMu_SFHo.pdfapplication/pdfMR_Bayes_NGWMu_comb.figapplication/x-xfigMR_Bayes_NGWMu_comb.pdfapplication/pdfMR_Bayes_NGW_SFHo.figapplication/x-xfigMR_Bayes_NGW_SFHo.pdfapplication/pdfMR_Bayes_NGW_comb.figapplication/x-xfigMR_Bayes_NGW_comb.pdfapplication/pdfMR_Bayes_NICER_SFHo.figapplication/x-xfigMR_Bayes_NICER_SFHo.pdfapplication/pdfMR_Bayes_NICER_comb.figapplication/x-xfigMR_Bayes_NICER_comb.pdfapplication/pdfMR_Bayes_prior_SFHo.figapplication/x-xfigMR_Bayes_prior_SFHo.pdfapplication/pdfMR_Bayes_prior_comb.figapplication/x-xfigMR_Bayes_prior_comb.pdfapplication/pdfMR_Bayesplot.mtext/x-matlabMR_data_DD2.matapplication/matlab-matMR_data_SFHo.matapplication/matlab-matMRcontour.mtext/x-matlabMRrelation.mtext/x-matlabM_quark.figapplication/x-xfigM_quark.pdfapplication/pdfManuscript.pdfapplication/pdfMaxmass_plot.mtext/x-matlabMmax_relation.epsapplication/postscriptMmax_relation.figapplication/x-xfigMqplot.mtext/x-matlabNICER2019_highres.matapplication/matlab-matNICER2021_highres.matapplication/matlab-matREADME.txttext/plainREADME.txttext/plainREADME.txttext/plainREADME.txttext/plainREADME.txttext/plainR_limits.figapplication/x-xfigR_limits.pdfapplication/pdfRungeKutta_nonrot.mtext/x-matlabSFHo_EoS.dattext/x-fixed-fieldSFHo_weights.matapplication/matlab-matSigma_290-Massesi.00.outapplication/octet-streamSigma_400-Massesi.00.outapplication/octet-streamSigma_500-Massesi.00.outapplication/octet-streamSigma_600-Massesi.00.outapplication/octet-streamSigma_700-Massesi.00.outapplication/octet-streamSly_EoS.dattext/x-fixed-fieldSly_EoS.mtext/x-matlabTM1_EoS.dattext/x-fixed-fieldTM1_EoS.mtext/x-matlabTM1e_EoS.dattext/x-fixed-fieldTM1e_EoS.mtext/x-matlabTOV_nonrot.mtext/x-matlabTogashi_EoS.dattext/x-fixed-fieldc2Calc.mtext/x-matlabcmapmod.matapplication/matlab-matcmapr.matapplication/matlab-mateLSM_zeroT.mtext/x-matlabeos.nbapplication/mathematicaeos.nbapplication/mathematicaeos.nbapplication/mathematicaeos.nbapplication/mathematicaeos.nbapplication/mathematicaeos.nbapplication/mathematicaeos.thermoapplication/octet-streameos.thermoapplication/octet-streameos.thermoapplication/octet-streameos.thermoapplication/octet-streameos.thermoapplication/octet-streameos.thermoapplication/octet-streameos.thermoapplication/octet-streameos.thermoapplication/octet-streameos.yqapplication/octet-streameos.yqapplication/octet-streameos.zipapplication/zipeos.zipapplication/zipeos.zipapplication/zipeos.zipapplication/zipeos.zipapplication/zipeos.zipapplication/zipeos.zipapplication/zipeos.zipapplication/zipisConswpQCD.mtext/x-matlabsource.gzapplication/x-gzip