# Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

Benjamin P Abbott, Rich Abbott, Thomas D Abbott, Sheelu Abraham, Fausto Acernese, Kendall Ackley, Carl Adams, Vaishali B Adya, Christoph Affeldt, Michalis Agathos, et al.

► To cite this version:

Benjamin P Abbott, Rich Abbott, Thomas D Abbott, Sheelu Abraham, Fausto Acernese, et al. Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant. Classical and Quantum Gravity, 2020, 37 (4), pp.045006. 10.1088/1361-6382/ab5f7c. hal-02999729

HAL Id: hal-02999729

https://hal.science/hal-02999729

Submitted on 17 Jan 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
PAPER • OPEN ACCESS

# Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

To cite this article: B P Abbott et al 2020 *Class. Quantum Grav.* **37** 045006

View the [article online](#) for updates and enhancements.

## You may also like

- *Operation and performance of the ATLAS semiconductor tracker in LHC Run 2*

The ATLAS collaboration, Georges Aad, Brad Abbott et al.

- *Performance of the ATLAS RPC detector and Level-1 muon barrel trigger at (s)=13 TeV*

The ATLAS collaboration, G. Aad, B. Abbott et al.

- *The ATLAS Fast Tracker system*

The ATLAS collaboration, G. Aad, B. Abbott et al.

## Recent citations

- *Effect of the symmetry energy on the secondary component of GW190814 as a neutron star*

Xuhao Wu et al

- *Numerical inside view of hypermassive remnant models for GW170817*

W. Kastaun and F. Ohme

- *Waveform systematics in the gravitational-wave inference of tidal parameters and equation of state from binary neutron-star signals*

Rossella Gamba et al

<figure><img src="image_2.png" /></figure>

IOP ebooks™

Bringing together innovative digital publishing with leading authors from the global scientific community.

Start exploring the collection-download the first chapter of every title for free.
Model comparison from LIGO–Virgo data
on GW170817’s binary components and
consequences for the merger remnant

B P Abbott¹, R Abbott¹, T D Abbott², S Abraham³,
F Acernese⁴,⁵, K Ackley⁶, C Adams⁷, V B Adya⁸, C Affeldt⁹,¹⁰,
M Agathos¹¹,¹², K Agatsuma¹³, N Aggarwal¹⁴, O D Aguiar¹⁵,
L Aiello¹⁶,¹⁷, A Ain³, P Ajith¹⁸, G Allen¹⁹, A Allocca²⁰,²¹,
M A Aloy²², P A Altin⁸, A Amato²³, S Anand¹, A Ananyeva¹,
S B Anderson¹, W G Anderson²⁴, S V Angelova²⁵, S Antier²⁶,
S Appert¹, K Arai¹, M C Araya¹, J S Areeda²⁷, M Arène²⁶,
N Arnaud²⁸,²⁹, S M Aronson³⁰, K G Arun³¹, S Ascenzi¹⁶,³²,
G Ashton⁶, S M Aston⁷, P Astone³³, F Aubin³⁴, P Aufmuth¹⁰,
K AultONeal³⁵, C Austin², V Avendano³⁶, A Avila-Alvarez²⁷,
S Babak²⁶, P Bacon²⁶, F Badaracco¹⁶,¹⁷, M K M Bader³⁷,
S Bae³⁸, J Baird²⁶, P T Baker³⁹, F Baldaccini⁴⁰,⁴¹,
G Ballardin²⁹, S W Ballmer⁴², A Bals³⁵, S Banagiri⁴³,
J C Barayoga¹, C Barbieri⁴⁴,⁴⁵, S E Barclay⁴⁶, B C Barish¹,
D Barker⁴⁷, K Barkett⁴⁸, S Barnum¹⁴, F Barone⁵,⁴⁹, B Barr⁴⁶,
L Barsotti¹⁴, M Barsuglia²⁶, D Barta⁵⁰, J Bartlett⁴⁷, I Bartos³⁰,
R Bassiri⁵¹, A Basti²⁰,²¹, M Bawaj⁴¹,⁵², J C Bayley⁴⁶,
M Bazzan⁵³,⁵⁴, B Bécsy⁵⁵, M Bejger²⁶,⁵⁶, I Belahcene²⁸,
A S Bell⁴⁶, D Beniwal⁵⁷, M G Benjamin³⁵, B K Berger⁵¹,
G Bergmann⁹,¹⁰, S Bernuzzi¹¹, C P L Berry⁵⁸, D Bersanetti⁵⁹,
A Bertolini³⁷, J Betzwieser⁷, R Bhandare⁶⁰, J Bidler²⁷,
E Biggs²⁴, I A Bilenko⁶¹, S A Bilgili³⁹, G Billingsley¹,
R Birney²⁵, O Birnholtz⁶², S Biscans¹,¹⁴, M Bischi⁶³,⁶⁴,
S Biscoveanu¹⁴, A Bisht¹⁰, M Bitossi²¹,²⁹, M A Bizouard⁶⁵,
J K Blackburn¹, J Blackman⁴⁸, C D Blair⁷, D G Blair⁶⁶,
R M Blair⁴⁷, S Bloemen⁶⁷, F Bobba⁶⁸,⁶⁹, N Bode⁹,¹⁰, M Boer⁶⁵,
Y Boetzel⁷⁰, G Bogaert⁶⁵, F Bondu⁷¹, R Bonnand³⁴,
P Booker⁹,¹⁰, B A Boom³⁷, R Bork¹, V Boschi²⁹, S Bose³,
V Bossilkov⁶⁶, J Bosveld⁶⁶, Y Bouffanais⁵³,⁵⁴, A Bozzi²⁹,
C Bradaschia²¹, P R Brady²⁴, A Bramley⁷, M Branchesi¹⁶,¹⁷,
J E Brau⁷², M Breschi¹¹, T Briant⁷³, J H Briggs⁴⁶,
F Brighenti⁶³,⁶⁴, A Brillet⁶⁵, M Brinkmann⁹,¹⁰, P Brockill²⁴,

<figure><img src="image_3.png" /></figure>

Original content from this work may be used under the terms of the Creative
Commons Attribution 3.0 licence. Any further distribution of this work must maintain
attribution to the author(s) and the title of the work, journal citation and DOI.
A F Brooks¹, J Brooks²⁹, D D Brown⁵⁷, S Brunett¹,
A Buikema¹⁴, T Bulik⁷⁴, H J Bulten³⁷,⁷⁵, A Buonanno⁷⁶,⁷⁷,
D Buskulic³⁴, C Buy²⁶, R L Byer⁵¹, M Cabero⁹,¹⁰, L Cadonati⁷⁸,
G Cagnoli⁷⁹, C Cahillane¹, J Calderón Bustillo⁶,
T A Callister¹, E Calloni⁵,⁸⁰, J B Camp⁸¹, W A Campbell⁶,
M Canepa⁵⁹,⁸², K C Cannon⁸³, H Cao⁵⁷, J Cao⁸⁴,
G Carapella⁶⁸,⁶⁹, F Carbognani²⁹, S Caride⁸⁵, M F Carney⁵⁸,
G Carullo²⁰,²¹, J Casanueva Diaz²¹, C Casentini³²,⁸⁶,
S Caudill³⁷, M Cavaglià⁸⁷,⁸⁸, F Cavalier²⁸, R Cavalieri²⁹,
G Cella²¹, P Cerdá-Durán²², E Cesarini³²,⁸⁹, O Chaibi⁶⁵,
K Chakravarti³, S J Chamberlin⁹⁰, M Chan⁴⁶, S Chao⁹¹,
P Charlton⁹², E A Chase⁵⁸, E Chassande-Mottin²⁶,
D Chatterjee²⁴, M Chaturvedi⁶⁰, K Chatziioannou⁹³,
B D Cheeseboro³⁹, H Y Chen⁹⁴, X Chen⁶⁶, Y Chen⁴⁸,
H-P Cheng³⁰, C K Cheong⁹⁵, H Y Chia³⁰, F Chiadini⁶⁹,⁹⁶,
A Chincarini⁵⁹, A Chiummo²⁹, G Cho⁹⁷, H S Cho⁹⁸, M Cho⁷⁷,
N Christensen⁶⁵,⁹⁹, Q Chu⁶⁶, S Chua⁷³, K W Chung⁹⁵,
S Chung⁶⁶, G Ciani⁵³,⁵⁴, M Cieślar⁵⁶, A A Ciobanu⁵⁷,
R Ciolfi⁵⁴,¹⁰⁰, F Cipriano⁶⁵, A Cirone⁵⁹,⁸², F Clara⁴⁷, J A Clark⁷⁸,
P Clearwater¹⁰¹, F Cleva⁶⁵, E Coccia¹⁶,¹⁷, P-F Cohadon⁷³,
D Cohen²⁸, M Colleoni¹⁰², C G Collette¹⁰³, C Collins¹³,
M Colpi⁴⁴,⁴⁵, L R Cominsky¹⁰⁴, M Constancio Jr¹⁵, L Conti⁵⁴,
S J Cooper¹³, P Corban⁷, T R Corbitt², I Cordero-Carrión¹⁰⁵,
S Corezzi⁴⁰,⁴¹, K R Corley¹⁰⁶, N Cornish⁵⁵, D Corre²⁸,
A Corsi⁸⁵, S Cortese²⁹, C A Costa¹⁵, R Cotesta⁷⁶,
M W Coughlin¹, S B Coughlin⁵⁸,¹⁰⁷, J-P Coulon⁶⁵,
S T Countryman¹⁰⁶, P Couvares¹, P B Covas¹⁰², E E Cowan⁷⁸,
D M Coward⁶⁶, M J Cowart⁷, D C Coyne¹, R Coyne¹⁰⁸,
J D E Creighton²⁴, T D Creighton¹⁰⁹, J Cripe², M Croquette⁷³,
S G Crowder¹¹⁰, T J Cullen², A Cumming⁴⁶, L Cunningham⁴⁶,
E Cuoco²⁹, T Dal Canton⁸¹, G Dálya¹¹¹, B D'Angelo⁵⁹,⁸²,
S L Danilishin⁹,¹⁰, S D'Antonio³², K Danzmann⁹,¹⁰,
A Dasgupta¹¹², C F Da Silva Costa³⁰, L E H Datrier⁴⁶,
V Dattilo²⁹, I Dave⁶⁰, M Davier²⁸, D Davis⁴², E J Daw¹¹³,
D DeBra⁵¹, M Deenadayalan³, J Degallaix²³, M De
Laurentis⁵,⁸⁰, S Deléglise⁷³, W Del Pozzo²⁰,²¹, L M DeMarchi⁵⁸,
N Demos¹⁴, T Dent¹¹⁴, R De Pietri¹¹⁵,¹¹⁶, R De Rosa⁵,⁸⁰, C De
Rossi²³,²⁹, R DeSalvo¹¹⁷, O de Varona⁹,¹⁰, S Dhurandhar³,
M C Díaz¹⁰⁹, T Dietrich³⁷, L Di Fiore⁵, C DiFronzo¹³, C Di
Giorgio⁶⁸,⁶⁹, F Di Giovanni²², M Di Giovanni¹¹⁸,¹¹⁹, T Di
Girolamo⁵,⁸⁰, A Di Lieto²⁰,²¹, B Ding¹⁰³, S Di Pace³³,¹²⁰, I Di
Palma³³,¹²⁰, F Di Renzo²⁰,²¹, A K Divakarla³⁰, A Dmitriev¹³,
Z Doctor⁹⁴, F Donovan¹⁴, K L Dooley⁸⁷,¹⁰⁷, S Doravari³,
I Dorrington¹⁰⁷, T P Downes²⁴, M Drago¹⁶,¹⁷, J C Driggers⁴⁷,
Z Du⁸⁴, J-G Ducoin²⁸, P Dupej⁴⁶, O Durante⁶⁸,⁶⁹, S E Dwyer⁴⁷,
P J Easter<sup>6</sup>, G Eddolls<sup>46</sup>, T B Edo<sup>113</sup>, A Effler<sup>7</sup>, P Ehrens<sup>1</sup>, J Eichholz<sup>8</sup>, S S Eikenberry<sup>30</sup>, M Eisenmann<sup>34</sup>, R A Eisenstein<sup>14</sup>, L Errico<sup>5,80</sup>, R C Essick<sup>94</sup>, H Estelles<sup>102</sup>, D Estevez<sup>34</sup>, Z B Etienne<sup>39</sup>, T Etzel<sup>1</sup>, M Evans<sup>14</sup>, T M Evans<sup>7</sup>, V Fafone<sup>16,32,86</sup>, S Fairhurst<sup>107</sup>, X Fan<sup>84</sup>, S Farinon<sup>59</sup>, B Farr<sup>72</sup>, W M Farr<sup>13</sup>, E J Fauchon-Jones<sup>107</sup>, M Favata<sup>36</sup>, M Fays<sup>113</sup>, M Fazio<sup>121</sup>, C Fee<sup>122</sup>, J Feicht<sup>1</sup>, M M Fejer<sup>51</sup>, F Feng<sup>26</sup>, A Fernandez-Galiana<sup>14</sup>, I Ferrante<sup>20,21</sup>, E C Ferreira<sup>15</sup>, T A Ferreira<sup>15</sup>, F Fidecaro<sup>20,21</sup>, I Fiori<sup>29</sup>, D Fiorucci<sup>16,17</sup>, M Fishbach<sup>94</sup>, R P Fisher<sup>123</sup>, J M Fishner<sup>14</sup>, R Fittipaldi<sup>69,124</sup>, M Fitz-Axen<sup>43</sup>, V Fiumara<sup>69,125</sup>, R Flaminio<sup>34,126</sup>, M Fletcher<sup>46</sup>, E Floden<sup>43</sup>, E Flynn<sup>27</sup>, H Fong<sup>83</sup>, J A Font<sup>22,127</sup>, P W F Forsyth<sup>8</sup>, J-D Fournier<sup>65</sup>, S Frasca<sup>33,120</sup>, F Frasconi<sup>21</sup>, Z Frei<sup>111</sup>, A Freise<sup>13</sup>, R Frey<sup>72</sup>, V Frey<sup>28</sup>, P Fritschel<sup>14</sup>, V V Frolov<sup>7</sup>, G Fronzè<sup>128</sup>, P Fulda<sup>30</sup>, M Fyffe<sup>7</sup>, H A Gabbard<sup>46</sup>, B U Gadre<sup>76</sup>, S M Gaebel<sup>13</sup>, J R Gair<sup>129</sup>, L Gammaitoni<sup>40</sup>, S G Gaonkar<sup>3</sup>, C García-Quirós<sup>102</sup>, F Garufi<sup>5,80</sup>, B Gateley<sup>47</sup>, S Gaudio<sup>35</sup>, G Gaur<sup>130</sup>, V Gayathri<sup>131</sup>, G Gemme<sup>59</sup>, E Genin<sup>29</sup>, A Gennai<sup>21</sup>, D George<sup>19</sup>, J George<sup>60</sup>, L Gergely<sup>132</sup>, S Ghonge<sup>78</sup>, Abhirup Ghosh<sup>76</sup>, Archisman Ghosh<sup>37</sup>, S Ghosh<sup>24</sup>, B Giacomazzo<sup>118,119</sup>, J A Giaime<sup>2,7</sup>, K D Giardina<sup>7</sup>, D R Gibson<sup>133</sup>, K Gill<sup>106</sup>, L Glover<sup>134</sup>, J Gniesmer<sup>135</sup>, P Godwin<sup>90</sup>, E Goetz<sup>47</sup>, R Goetz<sup>30</sup>, B Goncharov<sup>6</sup>, G González<sup>2</sup>, J M Gonzalez Castro<sup>20,21</sup>, A Gopakumar<sup>136</sup>, S E Gossan<sup>1</sup>, M Gosselin<sup>20,21,29</sup>, R Gouaty<sup>34</sup>, B Grace<sup>8</sup>, A Grado<sup>5,137</sup>, M Granata<sup>23</sup>, A Grant<sup>46</sup>, S Gras<sup>14</sup>, P Grassia<sup>1</sup>, C Gray<sup>47</sup>, R Gray<sup>46</sup>, G Greco<sup>63,64</sup>, A C Green<sup>30</sup>, R Green<sup>107</sup>, E M Gretarsson<sup>35</sup>, A Grimaldi<sup>118,119</sup>, S J Grimm<sup>16,17</sup>, P Groot<sup>67</sup>, H Grote<sup>107</sup>, S Grunewald<sup>76</sup>, P Gruning<sup>28</sup>, G M Guidi<sup>63,64</sup>, H K Gulati<sup>112</sup>, Y Guo<sup>37</sup>, A Gupta<sup>90</sup>, Anchal Gupta<sup>1</sup>, P Gupta<sup>37</sup>, E K Gustafson<sup>1</sup>, R Gustafson<sup>138</sup>, L Haegel<sup>102</sup>, O Halim<sup>16,17</sup>, B R Hall<sup>139</sup>, E D Hall<sup>14</sup>, E Z Hamilton<sup>107</sup>, G Hammond<sup>46</sup>, M Haney<sup>70</sup>, M M Hanke<sup>9,10</sup>, J Hanks<sup>47</sup>, C Hanna<sup>90</sup>, M D Hannam<sup>107</sup>, O A Hannuksela<sup>95</sup>, T J Hansen<sup>35</sup>, J Hanson<sup>7</sup>, T Harder<sup>65</sup>, T Hardwick<sup>2</sup>, K Haris<sup>18</sup>, J Harms<sup>16,17</sup>, G M Harry<sup>140</sup>, I W Harry<sup>141</sup>, R K Hasskew<sup>7</sup>, C J Haster<sup>14</sup>, K Haughian<sup>46</sup>, F J Hayes<sup>46</sup>, J Healy<sup>62</sup>, A Heidmann<sup>73</sup>, M C Heintze<sup>7</sup>, H Heitmann<sup>65</sup>, F Hellman<sup>142</sup>, P Hello<sup>28</sup>, G Hemming<sup>29</sup>, M Hendry<sup>46</sup>, I S Heng<sup>46</sup>, J Hennig<sup>9,10</sup>, Francisco Hernandez Vivanco<sup>6</sup>, M Heurs<sup>9,10</sup>, S Hild<sup>46</sup>, T Hinderer<sup>37,143,144</sup>, W C G Ho<sup>145</sup>, S Hochheim<sup>9,10</sup>, D Hofman<sup>23</sup>, A M Holgado<sup>19</sup>, N A Holland<sup>8</sup>, K Holt<sup>7</sup>, D E Holz<sup>94</sup>, P Hopkins<sup>107</sup>, C Horst<sup>24</sup>, J Hough<sup>46</sup>, E J Howell<sup>66</sup>, C G Hoy<sup>107</sup>, Y Huang<sup>14</sup>, M T Hübner<sup>6</sup>, E A Huerta<sup>19</sup>, D Huet<sup>28</sup>, B Hughey<sup>35</sup>, V Hui<sup>34</sup>, S Husa<sup>102</sup>, S H Huttner<sup>46</sup>, T Huynh-Dinh<sup>7</sup>, B Idzkowski<sup>74</sup>, A less<sup>32,86</sup>,
H Inchauspe³⁰, C Ingram⁵⁷, R Inta⁸⁵, G Intini³³,¹²⁰, B Irwin¹²², H N Isa⁴⁶, J-M Isac⁷³, M Isi¹⁴, B R Iyer¹⁸, T Jacqmin⁷³, S J Jadhav¹⁴⁶, K Jani⁷⁸, N N Janthalur¹⁴⁶, P Jaranowski¹⁴⁷, D Jariwala³⁰, A C Jenkins¹⁴⁸, J Jiang³⁰, D S Johnson¹⁹, A W Jones¹³, D I Jones¹⁴⁹, J D Jones⁴⁷, R Jones⁴⁶, R J G Jonker³⁷, L Ju⁶⁶, J Junker⁹,¹⁰, C V Kalaghatgi¹⁰⁷, V Kalogera⁵⁸, B Kamai¹, S Kandhasamy³, G Kang³⁸, J B Kanner¹, S J Kapadia²⁴, S Karki⁷², R Kashyap¹⁸, M Kasprzack¹, W Kastaun⁹,¹⁰, S Katsanevas²⁹, E Katsavounidis¹⁴, W Katzman⁷, S Kaufer¹⁰, K Kawabe⁴⁷, N V Keerthana³, F Kéfélian⁶⁵, D Keitel¹⁴¹, R Kennedy¹¹³, J S Key¹⁵⁰, F Y Khalili⁶¹, I Khan¹⁶,³², S Khan⁹,¹⁰, E A Khazanov¹⁵¹, N Khetan¹⁶,¹⁷, M Khursheed⁶⁰, N Kijbunchoo⁸, Chunglee Kim¹⁵², J C Kim¹⁵³, K Kim⁹⁵, W Kim⁵⁷, W S Kim¹⁵⁴, Y-M Kim¹⁵⁵, C Kimball⁵⁸, P J King⁴⁷, M Kinley-Hanlon⁴⁶, R Kirchhoff⁹,¹⁰, J S Kissel⁴⁷, L Kleybolte¹³⁵, J H Klika²⁴, S Klimenko³⁰, T D Knowles³⁹, P Koch⁹,¹⁰, S M Koehlenbeck⁹,¹⁰, G Koekoek³⁷,¹⁵⁶, S Koley³⁷, V Kondrashov¹, A Kontos¹⁵⁷, N Koper⁹,¹⁰, M Korobko¹³⁵, W Z Korth¹, M Kovalam⁶⁶, D B Kozak¹, C Krämer⁹,¹⁰, V Kringel⁹,¹⁰, N Krishnendu³¹, A Królak¹⁵⁸,¹⁵⁹, N Krupinski²⁴, G Kuehn⁹,¹⁰, A Kumar¹⁴⁶, P Kumar¹⁶⁰, Rahul Kumar⁴⁷, Rakesh Kumar¹¹², L Kuo⁹¹, A Kutynia¹⁵⁸, S Kwang²⁴, B D Lackey⁷⁶, D Laghi²⁰,²¹, K H Lai⁹⁵, T L Lam⁹⁵, M Landry⁴⁷, P Landry⁹⁴, B B Lane¹⁴, R N Lang¹⁶¹, J Lange⁶², B Lantz⁵¹, R K Lanza¹⁴, A Lartaux-Vollard²⁸, P D Lasky⁶, M Laxen⁷, A Lazzarini¹, C Lazzaro⁵⁴, P Leaci¹²⁰,³³, S Leavey⁹,¹⁰, Y K Lecoeuche⁴⁷, C H Lee⁹⁸, H K Lee¹⁶², H M Lee¹⁶³, H W Lee¹⁵³, J Lee⁹⁷, K Lee⁴⁶, J Lehmann⁹,¹⁰, A K Lenon³⁹, N Leroy²⁸, N Letendre³⁴, Y Levin⁶, A Li⁹⁵, J Li⁸⁴, K J L Li⁹⁵, T G F Li⁹⁵, X Li⁴⁸, F Lin⁶, F Linde³⁷,¹⁶⁴, S D Linker¹³⁴, T B Littenberg¹⁶⁵, J Liu⁶⁶, X Liu²⁴, M Llorens-Monteagudo²², R K L Lo⁹⁵,¹, L T London¹⁴, A Longo¹⁶⁶,¹⁶⁷, M Lorenzini¹⁶,¹⁷, V Loriette¹⁶⁸, M Lormand⁷, G Losurdo²¹, J D Lough⁹,¹⁰, C O Lousto⁶², G Lovelace²⁷, M E Lower¹⁶⁹, H Lück⁹,¹⁰, D Lumaca³²,⁸⁶, A P Lundgren¹⁴¹, R Lynch¹⁴, Y Ma⁴⁸, R Macas¹⁰⁷, S Macfoy²⁵, M MacInnis¹⁴, D M Macleod¹⁰⁷, A Macquet⁶⁵, I Magaña Hernandez²⁴, F Magaña-Sandoval³⁰, R M Magee⁹⁰, E Majorana³³, I Maksimovic¹⁶⁸, A Malik⁶⁰, N Man⁶⁵, V Mandic⁴³, V Mangano³³,⁴⁶,¹²⁰, G L Mansell¹⁴,⁴⁷, M Manske²⁴, M Mantovani²⁹, M Mapelli⁵³,⁵⁴, F Marchesoni⁴¹,⁵², F Marion³⁴, S Márka¹⁰⁶, Z Márka¹⁰⁶, C Markakis¹⁹, A S Markosyan⁵¹, A Markowitz¹, E Maros¹, A Marquina¹⁰⁵, S Marsat²⁶, F Martelli⁶³,⁶⁴, I W Martin⁴⁶, R M Martin³⁶, V Martinez⁷⁹, D V Martynov¹³, H Masalehdan¹³⁵, K Mason¹⁴,
E Massera¹¹³, A Masserot³⁴, T J Massinger¹, M Masso-Reid⁴⁶, S Mastrogiovanni²⁶, A Matas⁷⁶, F Matichard¹,¹⁴, L Matone¹⁰⁶, N Mavalvala¹⁴, J J McCann⁶⁶, R McCarthy⁴⁷, D E McClelland⁸, S McCormick⁷, L McCuller¹⁴, S C McGuire¹⁷⁰, C McIsaac¹⁴¹, J McIver¹, D J McManus⁸, T McRae⁸, S T McWilliams³⁹, D Meacher²⁴, G D Meadors⁶, M Mehmet⁹,¹⁰, A K Mehta¹⁸, J Meidam³⁷, E Mejuto Villa⁶⁹,¹¹⁷, A Melatos¹⁰¹, G Mendell⁴⁷, R A Mercer²⁴, L Mereni²³, K Merfeld⁷², E L Merilh⁴⁷, M Merzougui⁶⁵, S Meshkov¹, C Messenger⁴⁶, C Messick⁹⁰, F Messina⁴⁴,⁴⁵, R Metzdorff⁷³, P M Meyers¹⁰¹, F Meylahn⁹,¹⁰, A Miani¹¹⁸,¹¹⁹, H Miao¹³, C Michel²³, H Middleton¹⁰¹, L Milano⁸⁰,⁵, A L Miller³⁰,³³,¹²⁰, M Millhouse¹⁰¹, J C Mills¹⁰⁷, M C Milovich-Goff¹³⁴, O Minazzoli⁶⁵,¹⁷¹, Y Minenkov³², A Mishkin³⁰, C Mishra¹⁷², T Mistry¹¹³, S Mitra³, V P Mitrofanov⁶¹, G Mitselmakher³⁰, R Mittleman¹⁴, G Mo⁹⁹, D Moffa¹²², K Mogushi⁸⁷, S R P Mohapatra¹⁴, M Molina-Ruiz¹⁴², M Mondin¹³⁴, M Montani⁶³,⁶⁴, C J Moore¹³, D Moraru⁴⁷, F Morawski⁵⁶, G Moreno⁴⁷, S Morisaki⁸³, B Mours³⁴, C M Mow-Lowry¹³, F Muciaccia³³,¹²⁰, Arunava Mukherjee⁹,¹⁰, D Mukherjee²⁴, S Mukherjee¹⁰⁹, Subroto Mukherjee¹¹², N Mukund³,⁹,¹⁰, A Mullavey⁷, J Munch⁵⁷, E A Muñiz⁴², M Muratore³⁵, P G Murray⁴⁶, A Nag ar⁸⁹,¹²⁸,¹⁷³, I Nardecchia³²,⁸⁶, L Naticchioni³³,¹²⁰, R K Nayak¹⁷⁴, B F Neil⁶⁶, J Neilson⁶⁹,¹¹⁷, G Nelemans³⁷,⁶⁷, T J N Nelson⁷, M Nery⁹,¹⁰, A Neunzert¹³⁸, L Nevin¹, K Y Ng¹⁴, S Ng⁵⁷, C Nguyen²⁶, P Nguyen⁷², D Nichols³⁷,¹⁴³, S A Nichols², S Nissanke³⁷,¹⁴³, F Nocera²⁹, C North¹⁰⁷, L K Nuttall¹⁴¹, M Obergaulinger²²,¹⁷⁵, J Oberling⁴⁷, B D O'Brien³⁰, G Oganesyan¹⁶,¹⁷, G H Ogin¹⁷⁶, J J Oh¹⁵⁴, S H Oh¹⁵⁴, F Ohme⁹,¹⁰, H Ohta⁸³, M A Okada¹⁵, M Oliver¹⁰², P Oppermann⁹,¹⁰, Richard J Oram⁷, B O'Reilly⁷, R G Ormiston⁴³, L F Ortega³⁰, R O'Shaughnessy⁶², S Ossokine⁷⁶, D J Ottaway⁵⁷, H Overmier⁷, B J Owen⁸⁵, A E Pace⁹⁰, G Pagano²⁰,²¹, M A Page⁶⁶, G Pagliaroli¹⁶,¹⁷, A Pai¹³¹, S A Pai⁶⁰, J R Palamos⁷², O Palashov¹⁵¹, C Palomba³³, H Pan⁹¹, P K Panda¹⁴⁶, P T H Pang³⁷,⁹⁵, C Pankow⁵⁸, F Pannarale³³,¹²⁰, B C Pant⁶⁰, F Paoletti²¹, A Paoli²⁹, A Parida³, W Parker⁷,¹⁷⁰, D Pascucci³⁷,⁴⁶, A Pasqualetti²⁹, R Passaqueti²⁰,²¹, D Passuello²¹, M Patil¹⁵⁹, B Patricelli²⁰,²¹, E Payne⁶, B L Pearlstone⁴⁶, T C Pechsiri³⁰, A J Pedersen⁴², M Pedraza¹, R Pedurand²³,¹⁷⁷, A Pele⁷, S Penn¹⁷⁸, A Perego¹¹⁸,¹¹⁹, C J Perez⁴⁷, C Périgois³⁴, A Perreca¹¹⁸,¹¹⁹, J Petermann¹³⁵, H P Pfeiffer⁷⁶, M Phelps⁹,¹⁰, K S Phukon³, O J Piccinni³³,¹²⁰, M Pichot⁶⁵, F Piergiovanni⁶³,⁶⁴, V Pierro⁶⁹,¹¹⁷, G Pillant²⁹, L Pinard²³,
I M Pinto<sup>69,89,117</sup>, M Pirello<sup>47</sup>, M Pitkin<sup>46</sup>, W Plastino<sup>166,167</sup>, R Poggiani<sup>20,21</sup>, D Y T Pong<sup>95</sup>, S Ponrathnam<sup>3</sup>, P Popolizio<sup>29</sup>, E K Porter<sup>26</sup>, J Powell<sup>169</sup>, A K Prajapati<sup>112</sup>, J Prasad<sup>3</sup>, K Prasai<sup>51</sup>, R Prasanna<sup>146</sup>, G Pratten<sup>102</sup>, T Prestegard<sup>24</sup>, M Principe<sup>89,69,117</sup>, G A Prodi<sup>118,119</sup>, L Prokhorov<sup>13</sup>, M Punturo<sup>41</sup>, P Puppo<sup>33</sup>, M Pürrer<sup>76</sup>, H Qi<sup>107</sup>, V Quetschke<sup>109</sup>, P J Quinonez<sup>35</sup>, F J Raab<sup>47</sup>, G Raaijmakers<sup>37,143</sup>, H Radkins<sup>47</sup>, N Radulesco<sup>65</sup>, P Raffai<sup>111</sup>, S Raja<sup>60</sup>, C Rajan<sup>60</sup>, B Rajbhandari<sup>85</sup>, M Rakhmanov<sup>109</sup>, K E Ramirez<sup>109</sup>, A Ramos-Buades<sup>102</sup>, Javed Rana<sup>3</sup>, K Rao<sup>58</sup>, P Rapagnani<sup>33,120</sup>, V Raymond<sup>107</sup>, M Razzano<sup>20,21</sup>, J Read<sup>27</sup>, T Regimbau<sup>34</sup>, L Rei<sup>59</sup>, S Reid<sup>25</sup>, D H Reitzel<sup>1,30</sup>, P Rettegno<sup>128,179</sup>, F Ricci<sup>120,33</sup>, C J Richardson<sup>35</sup>, J W Richardson<sup>1</sup>, P M Ricker<sup>19</sup>, G Riemenschneider<sup>128,179</sup>, K Riles<sup>138</sup>, M Rizzo<sup>58</sup>, N A Robertson<sup>1,46</sup>, F Robinet<sup>28</sup>, A Rocchi<sup>32</sup>, L Rolland<sup>34</sup>, J G Rollins<sup>1</sup>, V J Roma<sup>72</sup>, M Romanelli<sup>71</sup>, R Romano<sup>4,5</sup>, C L Romel<sup>47</sup>, J H Romie<sup>7</sup>, C A Rose<sup>24</sup>, D Rose<sup>27</sup>, K Rose<sup>122</sup>, D Rosińska<sup>74</sup>, S G Rosofsky<sup>19</sup>, M P Ross<sup>180</sup>, S Rowan<sup>46</sup>, A Rüdiger<sup>9,10</sup>, P Ruggi<sup>29</sup>, G Rutins<sup>133</sup>, K Ryan<sup>47</sup>, S Sachdev<sup>90</sup>, T Sadecki<sup>47</sup>, M Sakellariadou<sup>148</sup>, O S Salafia<sup>44,45,181</sup>, L Salconi<sup>29</sup>, M Saleem<sup>31</sup>, A Samajdar<sup>37</sup>, L Sammut<sup>6</sup>, E J Sanchez<sup>1</sup>, L E Sanchez<sup>1</sup>, N Sanchis-Gual<sup>182</sup>, J R Sanders<sup>183</sup>, K A Santiago<sup>36</sup>, E Santos<sup>65</sup>, N Sarin<sup>6</sup>, B Sassolas<sup>23</sup>, B S Sathyaprakash<sup>90,107</sup>, O Sauter<sup>138,34</sup>, R L Savage<sup>47</sup>, P Schale<sup>72</sup>, M Scheel<sup>48</sup>, J Scheuer<sup>58</sup>, P Schmidt<sup>13,67</sup>, R Schnabel<sup>135</sup>, R M S Schofield<sup>72</sup>, A Schönbeck<sup>135</sup>, E Schreiber<sup>9,10</sup>, B W Schulte<sup>9,10</sup>, B F Schutz<sup>107</sup>, J Scott<sup>46</sup>, S M Scott<sup>8</sup>, E Seidel<sup>19</sup>, D Sellers<sup>7</sup>, A S Sengupta<sup>184</sup>, N Sennett<sup>76</sup>, D Sentenac<sup>29</sup>, V Sequino<sup>59</sup>, A Sergeev<sup>151</sup>, Y Setyawati<sup>9,10</sup>, D A Shaddock<sup>8</sup>, T Shaffer<sup>47</sup>, M S Shahriar<sup>58</sup>, M B Shaner<sup>134</sup>, A Sharma<sup>16,17</sup>, P Sharma<sup>60</sup>, P Shawhan<sup>77</sup>, H Shen<sup>19</sup>, R Shink<sup>185</sup>, D H Shoemaker<sup>14</sup>, D M Shoemaker<sup>78</sup>, K Shukla<sup>142</sup>, S ShyamSundar<sup>60</sup>, K Siellez<sup>78</sup>, M Sieniawska<sup>56</sup>, D Sigg<sup>47</sup>, L P Singer<sup>81</sup>, D Singh<sup>90</sup>, N Singh<sup>74</sup>, A Singhal<sup>16,33</sup>, A M Sintes<sup>102</sup>, S Sitmukhambetov<sup>109</sup>, V Skliris<sup>107</sup>, B J J Slagmolen<sup>8</sup>, T J Slaven-Blair<sup>66</sup>, J R Smith<sup>27</sup>, R J E Smith<sup>6</sup>, S Somala<sup>186</sup>, E J Son<sup>154</sup>, S Soni<sup>2</sup>, B Sorazu<sup>46</sup>, F Sorrentino<sup>59</sup>, T Souradeep<sup>3</sup>, E Sowell<sup>85</sup>, A P Spencer<sup>46</sup>, M Spera<sup>53,54</sup>, A K Srivastava<sup>112</sup>, V Srivastava<sup>42</sup>, K Staats<sup>58</sup>, C Stachie<sup>65</sup>, M Standke<sup>9,10</sup>, D A Steer<sup>26</sup>, M Steinke<sup>9,10</sup>, J Steinlechner<sup>46,135</sup>, S Steinlechner<sup>135</sup>, D Steinmeyer<sup>9,10</sup>, S P Stevenson<sup>169</sup>, D Stocks<sup>51</sup>, R Stone<sup>109</sup>, D J Stops<sup>13</sup>, K A Strain<sup>46</sup>, G Stratta<sup>64,187</sup>, S E Strigin<sup>61</sup>, A Strunk<sup>47</sup>, R Sturani<sup>188</sup>, A L Stuver<sup>189</sup>, V Sudhir<sup>14</sup>, T Z Summerscales<sup>190</sup>,
L Sun<fnref n="1" />, S Sunil¹¹², A Sur⁵⁶, J Suresh⁸³, P J Sutton¹⁰⁷,
B L Swinkels³⁷, M J Szczepańczyk³⁵, M Tacca³⁷, S C Tait⁴⁶,
C Talbot⁶, D B Tanner³⁰, D Tao<fnref n="1" />, M Tápai¹³², A Tapia²⁷,
J D Tasson⁹⁹, R Taylor<fnref n="1" />, R Tenorio¹⁰², L Terkowski¹³⁵,
M Thomas⁷, P Thomas⁴⁷, S R Thondapu⁶⁰, K A Thorne⁷,
E Thrane⁶, Shubhanshu Tiwari¹¹⁸,¹¹⁹, Srishti Tiwari¹³⁶,
V Tiwari¹⁰⁷, K Toland⁴⁶, M Tonelli²⁰,²¹, Z Tornasi⁴⁶, A Torres-
Forné¹⁹¹, C I Torrie<fnref n="1" />, D Töyrä¹³, F Travasso²⁹,⁴¹, G Traylor⁷,
M C Tringali⁷⁴, A Tripathee¹³⁸, A Trovato²⁶, L Trozzo²¹,¹⁹²,
K W Tsang³⁷, M Tse¹⁴, R Tso⁴⁸, L Tsukada⁸³, D Tsuna⁸³,
T Tsutsui⁸³, D Tuyenbayev¹⁰⁹, K Ueno⁸³, D Ugolini¹⁹³,
C S Unnikrishnan¹³⁶, A L Urban<fnref n="2" />, S A Usman⁹⁴,
H Vahlbruch¹⁰, G Vajente<fnref n="1" />, G Valdes<fnref n="2" />, M Valentini¹¹⁸,¹¹⁹, N van
Bakel³⁷, M van Beuzekom³⁷, J F J van den Brand³⁷,⁷⁵, C Van
Den Broeck³⁷,¹⁹⁴, D C Vander-Hyde⁴², L van der Schaaf³⁷,
J V VanHeijningen⁶⁶, A A van Veggel⁴⁶, M Vardaro⁵³,⁵⁴,
V Varma⁴⁸, S Vass<fnref n="1" />, M Vasúth⁵⁰, A Vecchio¹³, G Vedovato⁵⁴,
J Veitch⁴⁶, P J Veitch⁵⁷, K Venkateswara¹⁸⁰, G Venugopalan<fnref n="1" />,
D Verkindt³⁴, F Vetrano⁶³,⁶⁴, A Viceré⁶³,⁶⁴, A D Viets²⁴,
S Vinciguerra¹³, D J Vine¹³³, J-Y Vinet⁶⁵, S Vitale¹⁴, T Vo⁴²,
H Vocca⁴⁰,⁴¹, C Vorvick⁴⁷, S P Vyatchanin⁶¹, A R Wade<fnref n="1" />,
L E Wade¹²², M Wade¹²², R Walet³⁷, M Walker²⁷, L Wallace<fnref n="1" />,
S Walsh²⁴, H Wang¹³, J Z Wang¹³⁸, S Wang¹⁹, W H Wang¹⁰⁹,
Y F Wang⁹⁵, R L Ward⁸, Z A Warden³⁵, J Warner⁴⁷, M Was³⁴,
J Watchi¹⁰³, B Weaver⁴⁷, L-W Wei⁹,¹⁰, M Weinert⁹,¹⁰,
A J Weinstein<fnref n="1" />, R Weiss¹⁴, F Wellmann⁹,¹⁰, L Wen⁶⁶,
E K Wessel¹⁹, P Weßels⁹,¹⁰, J W Westhouse³⁵, K Wette⁸,
J T Whelan⁶², B F Whiting³⁰, C Whittle¹⁴, D M Wilken⁹,¹⁰,
D Williams⁴⁶, A R Williamson³⁷,¹⁴³, J L Willis<fnref n="1" />, B Willke¹⁰,⁹,
W Winkler⁹,¹⁰, C C Wipf<fnref n="1" />, H Wittel⁹,¹⁰, G Woan⁴⁶, J Woehler⁹,¹⁰,
J K Wofford⁶², J L Wright⁴⁶, D S Wu⁹,¹⁰, D M Wysocki⁶²,
S Xiao<fnref n="1" />, R Xu¹¹⁰, H Yamamoto<fnref n="1" />, C C Yancey⁷⁷, L Yang¹²¹,
Y Yang³⁰, Z Yang⁴³, M J Yap⁸, M Yazback³⁰, D W Yeeles¹⁰⁷,
Hang Yu¹⁴, Haocun Yu¹⁴, S H R Yuen⁹⁵, A K Zadrożny¹⁰⁹,
A Zadrożny¹⁵⁸, M Zanolin³⁵, F Zappa⁷⁶, T Zelenova²⁹, J-
P Zendri⁵⁴, M Zevin⁵⁸, J Zhang⁶⁶, L Zhang<fnref n="1" />, T Zhang⁴⁶,
C Zhao⁶⁶, G Zhao¹⁰³, M Zhou⁵⁸, Z Zhou⁵⁸, X J Zhu⁶,
A B Zimmerman¹⁹⁵, Y Zlochower⁶², M E Zucker<fnref n="1" />,¹⁴, J Zweizig<fnref n="1" />
(The LIGO Scientific Collaboration and The Virgo
Collaboration)

<footnote n="1">¹ LIGO, California Institute of Technology, Pasadena, CA 91125,
United States of America</footnote>

<footnote n="2">² Louisiana State University, Baton Rouge, LA 70803, United States of America</footnote>

<footnote n="3">³ Inter-University Centre for Astronomy and Astrophysics, Pune 411007, India</footnote>

<footnote n="4">⁴ Dipartimento di Farmacia, Università di Salerno, I-84084 Fisciano, Salerno, Italy</footnote>
⁵ INFN, Sezione di Napoli, Complesso Universitario di Monte S.Angelo, I-80126 Napoli, Italy

⁶ OzGrav, School of Physics & Astronomy, Monash University, Clayton 3800, Victoria, Australia

⁷ LIGO Livingston Observatory, Livingston, LA 70754, United States of America

⁸ OzGrav, Australian National University, Canberra, Australian Capital Territory 0200, Australia

⁹ Max Planck Institute for Gravitational Physics (Albert Einstein Institute), D-30167 Hannover, Germany

¹⁰ Leibniz Universität Hannover, D-30167 Hannover, Germany

¹¹ Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universität Jena, D-07743 Jena, Germany

¹² University of Cambridge, Cambridge CB2 1TN, United Kingdom

¹³ University of Birmingham, Birmingham B15 2TT, United Kingdom

¹⁴ LIGO, Massachusetts Institute of Technology, Cambridge, MA 02139, United States of America

¹⁵ Instituto Nacional de Pesquisas Espaciais, 12227-010 São José dos Campos, São Paulo, Brazil

¹⁶ Gran Sasso Science Institute (GSSI), I-67100 L'Aquila, Italy

¹⁷ INFN, Laboratori Nazionali del Gran Sasso, I-67100 Assergi, Italy

¹⁸ International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, Bengaluru 560089, India

¹⁹ NCSA, University of Illinois at Urbana-Champaign, Urbana, IL 61801, United States of America

²⁰ Università di Pisa, I-56127 Pisa, Italy

²¹ INFN, Sezione di Pisa, I-56127 Pisa, Italy

²² Departamento de Astronomía y Astrofísica, Universitat de València, E-46100 Burjassot, València, Spain

²³ Laboratoire des Matériaux Avancés (LMA), CNRS/IN2P3, F-69622 Villeurbanne, France

²⁴ University of Wisconsin-Milwaukee, Milwaukee, WI 53201, United States of America

²⁵ SUPA, University of Strathclyde, Glasgow G1 1XQ, United Kingdom

²⁶ APC, AstroParticule et Cosmologie, Université Paris Diderot, CNRS/IN2P3, CEA/Irfu, Observatoire de Paris, Sorbonne Paris Cité, F-75205 Paris Cedex 13, France

²⁷ California State University Fullerton, Fullerton, CA 92831, United States of America

²⁸ LAL, Univ. Paris-Sud, CNRS/IN2P3, Université Paris-Saclay, F-91898 Orsay, France

²⁹ European Gravitational Observatory (EGO), I-56021 Cascina, Pisa, Italy

³⁰ University of Florida, Gainesville, FL 32611, United States of America

³¹ Chennai Mathematical Institute, Chennai 603103, India

³² INFN, Sezione di Roma Tor Vergata, I-00133 Roma, Italy

³³ INFN, Sezione di Roma, I-00185 Roma, Italy

³⁴ Laboratoire d'Annecy de Physique des Particules (LAPP), Univ. Grenoble Alpes, Université Savoie Mont Blanc, CNRS/IN2P3, F-74941 Annecy, France

³⁵ Embry-Riddle Aeronautical University, Prescott, AZ 86301, United States of America

³⁶ Montclair State University, Montclair, NJ 07043, United States of America

³⁷ Nikhef, Science Park 105, 1098 XG Amsterdam, The Netherlands

³⁸ Korea Institute of Science and Technology Information, Daejeon 34141, Republic of Korea

³⁹ West Virginia University, Morgantown, WV 26506, United States of America
40 Università di Perugia, I-06123 Perugia, Italy

41 INFN, Sezione di Perugia, I-06123 Perugia, Italy

42 Syracuse University, Syracuse, NY 13244, United States of America

43 University of Minnesota, Minneapolis, MN 55455, United States of America

44 Università degli Studi di Milano-Bicocca, I-20126 Milano, Italy

45 INFN, Sezione di Milano-Bicocca, I-20126 Milano, Italy

46 SUPA, University of Glasgow, Glasgow G12 8QQ, United Kingdom

47 LIGO Hanford Observatory, Richland, WA 99352, United States of America

48 Caltech CaRT, Pasadena, CA 91125, United States of America

49 Dipartimento di Medicina, Chirurgia e Odontoiatria 'Scuola Medica Salernitana', Università di Salerno, I-84081 Baronissi, Salerno, Italy

50 Wigner RCP, RMKI, H-1121 Budapest, Konkoly Thege Miklós út 29-33, Hungary

51 Stanford University, Stanford, CA 94305, United States of America

52 Università di Camerino, Dipartimento di Fisica, I-62032 Camerino, Italy

53 Università di Padova, Dipartimento di Fisica e Astronomia, I-35131 Padova, Italy

54 INFN, Sezione di Padova, I-35131 Padova, Italy

55 Montana State University, Bozeman, MT 59717, United States of America

56 Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, 00-716, Warsaw, Poland

57 OzGrav, University of Adelaide, Adelaide, South Australia 5005, Australia

58 Center for Interdisciplinary Exploration & Research in Astrophysics (CIERA), Northwestern University, Evanston, IL 60208, United States of America

59 INFN, Sezione di Genova, I-16146 Genova, Italy

60 RRCAT, Indore, Madhya Pradesh 452013, India

61 Faculty of Physics, Lomonosov Moscow State University, Moscow 119991, Russia

62 Rochester Institute of Technology, Rochester, NY 14623, United States of America

63 Università degli Studi di Urbino 'Carlo Bo', I-61029 Urbino, Italy

64 INFN, Sezione di Firenze, I-50019 Sesto Fiorentino, Firenze, Italy

65 Artemis, Université Côte d'Azur, Observatoire Côte d'Azur, CNRS, CS 34229, F-06304 Nice Cedex 4, France

66 OzGrav, University of Western Australia, Crawley, Western Australia 6009, Australia

67 Department of Astrophysics/IMAPP, Radboud University Nijmegen, PO Box 9010, 6500 GL Nijmegen, The Netherlands

68 Dipartimento di Fisica 'E.R. Caianiello', Università di Salerno, I-84084 Fisciano, Salerno, Italy

69 INFN, Sezione di Napoli, Gruppo Collegato di Salerno, Complesso Universitario di Monte S. Angelo, I-80126 Napoli, Italy

70 Physik-Institut, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland

71 Univ Rennes, CNRS, Institut FOTON—UMR6082, F-3500 Rennes, France

72 University of Oregon, Eugene, OR 97403, United States of America

73 Laboratoire Kastler Brossel, Sorbonne Université, CNRS, ENS-Université PSL, Collège de France, F-75005 Paris, France

74 Astronomical Observatory Warsaw University, 00-478 Warsaw, Poland

75 VU University Amsterdam, 1081 HV Amsterdam, The Netherlands

76 Max Planck Institute for Gravitational Physics (Albert Einstein Institute), D-14476 Potsdam-Golm, Germany

77 University of Maryland, College Park, MD 20742, United States of America

78 School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, United States of America

79 Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
80 Università di Napoli ‘Federico II’, Complesso Universitario di Monte S.Angelo, I-80126 Napoli, Italy

81 NASA Goddard Space Flight Center, Greenbelt, MD 20771, United States of America

82 Dipartimento di Fisica, Università degli Studi di Genova, I-16146 Genova, Italy

83 RESCEU, University of Tokyo, Tokyo, 113-0033, Japan

84 Tsinghua University, Beijing 100084, People's Republic of China

85 Texas Tech University, Lubbock, TX 79409, United States of America

86 Università di Roma Tor Vergata, I-00133 Roma, Italy

87 The University of Mississippi, University, MS 38677, United States of America

88 Missouri University of Science and Technology, Rolla, MO 65409, United States of America

89 Museo Storico della Fisica e Centro Studi e Ricerche ‘Enrico Fermi’, I-00184 Roma, Italy

90 The Pennsylvania State University, University Park, PA 16802, United States of America

91 National Tsing Hua University, Hsinchu City, 30013 Taiwan, Republic of China

92 Charles Sturt University, Wagga Wagga, New South Wales 2678, Australia

93 Canadian Institute for Theoretical Astrophysics, University of Toronto, Toronto, Ontario M5S 3H8, Canada

94 University of Chicago, Chicago, IL 60637, United States of America

95 The Chinese University of Hong Kong, Shatin, NT, Hong Kong

96 Dipartimento di Ingegneria Industriale (DIIN), Università di Salerno, I-84084 Fisciano, Salerno, Italy

97 Seoul National University Seoul 08826, Republic of Korea

98 P United States of American National University, Busan 46241, Republic of Korea

99 Carleton College, Northfield, MN 55057, United States of America

100 INAF, Osservatorio Astronomico di Padova, I-35122 Padova, Italy

101 OzGrav, University of Melbourne, Parkville, Victoria 3010, Australia

102 Universitat de les Illes Balears, IAC3—IEEC, E-07122 Palma de Mallorca, Spain

103 Université Libre de Bruxelles, Brussels 1050, Belgium

104 Sonoma State University, Rohnert Park, CA 94928, United States of America

105 Departamento de Matemáticas, Universitat de València, E-46100 Burjassot, València, Spain

106 Columbia University, New York, NY 10027, United States of America

107 Cardiff University, Cardiff CF24 3AA, United Kingdom

108 University of Rhode Island, Kingston, RI 02881, United States of America

109 The University of Texas Rio Grande Valley, Brownsville, TX 78520, United States of America

110 Bellevue College, Bellevue, WA 98007, United States of America

111 MTA-ELTE Astrophysics Research Group, Institute of Physics, Eötvös University, Budapest 1117, Hungary

112 Institute for Plasma Research, Bhat, Gandhinagar 382428, India

113 The University of Sheffield, Sheffield S10 2TN, United Kingdom

114 IGFAE, Campus Sur, Universidade de Santiago de Compostela, 15782 Spain

115 Dipartimento di Scienze Matematiche, Fisiche e Informatiche, Università di Parma, I-43124 Parma, Italy

116 INFN, Sezione di Milano Bicocca, Gruppo Collegato di Parma, I-43124 Parma, Italy

117 Dipartimento di Ingegneria, Università del Sannio, I-82100 Benevento, Italy

118 Dipartimento di Fisica, Università di Trento, I-38123 Povo, Trento, Italy

119 INFN, Trento Institute for Fundamental Physics and Applications, I-38123 Povo, Trento, Italy

120 Università di Roma ‘La Sapienza’, I-00185 Roma, Italy
121 Colorado State University, Fort Collins, CO 80523, United States of America

122 Kenyon College, Gambier, OH 43022, United States of America

123 Christopher Newport University, Newport News, VA 23606, United States of America

124 CNR-SPIN, c/o Università di Salerno, I-84084 Fisciano, Salerno, Italy

125 Scuola di Ingegneria, Università della Basilicata, I-85100 Potenza, Italy

126 National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

127 Observatori Astronòmic, Universitat de València, E-46980 Paterna, València, Spain

128 INFN Sezione di Torino, I-10125 Torino, Italy

129 School of Mathematics, University of Edinburgh, Edinburgh EH9 3FD, United Kingdom

130 Institute Of Advanced Research, Gandhinagar 382426, India

131 Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India

132 University of Szeged, Dóm tér 9, Szeged 6720, Hungary

133 SUPA, University of the West of Scotland, Paisley PA1 2BE, United Kingdom

134 California State University, Los Angeles, 5151 State University Dr, Los Angeles, CA 90032, United States of America

135 Universität Hamburg, D-22761 Hamburg, Germany

136 Tata Institute of Fundamental Research, Mumbai 400005, India

137 INAF, Osservatorio Astronomico di Capodimonte, I-80131 Napoli, Italy

138 University of Michigan, Ann Arbor, MI 48109, United States of America

139 Washington State University, Pullman, WA 99164, United States of America

140 American University, Washington, D.C. 20016, United States of America

141 University of Portsmouth, Portsmouth, PO1 3FX, United Kingdom

142 University of California, Berkeley, CA 94720, United States of America

143 GRAPPA, Anton Pannekoek Institute for Astronomy and Institute for High-Energy Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands

144 Delta Institute for Theoretical Physics, Science Park 904, 1090 GL Amsterdam, The Netherlands

145 Department of Physics and Astronomy, Haverford College, Haverford, PA 19041, United States of America

146 Directorate of Construction, Services & Estate Management, Mumbai 400094 India

147 University of Białystok, 15-424 Białystok, Poland

148 King's College London, University of London, London WC2R 2LS, United Kingdom

149 University of Southampton, Southampton SO17 1BJ, United Kingdom

150 University of Washington Bothell, Bothell, WA 98011, United States of America

151 Institute of Applied Physics, Nizhny Novgorod, 603950, Russia

152 Ewha Womans University, Seoul 03760, Republic of Korea

153 Inje University Gimhae, South Gyeongsang 50834, Republic of Korea

154 National Institute for Mathematical Sciences, Daejeon 34047, Republic of Korea

155 Ulsan National Institute of Science and Technology, Ulsan 44919, Republic of Korea

156 Maastricht University, PO Box 616, 6200 MD Maastricht, The Netherlands

157 Bard College, 30 Campus Rd, Annandale-On-Hudson, NY 12504, United States of America

158 NCBJ, 05-400 Świerk-Otwock, Poland

159 Institute of Mathematics, Polish Academy of Sciences, 00656 Warsaw, Poland

160 Cornell University, Ithaca, NY 14850, United States of America
161 Hillsdale College, Hillsdale, MI 49242, United States of America

162 Hanyang University, Seoul 04763, Republic of Korea

163 Korea Astronomy and Space Science Institute, Daejeon 34055, Republic of Korea

164 Institute for High-Energy Physics, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands

165 NASA Marshall Space Flight Center, Huntsville, AL 35811, United States of America

166 Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, I-00146 Roma, Italy

167 INFN, Sezione di Roma Tre, I-00146 Roma, Italy

168 ESPCI, CNRS, F-75005 Paris, France

169 OzGrav, Swinburne University of Technology, Hawthorn VIC 3122, Australia

170 Southern University and A&M College, Baton Rouge, LA 70813, United States of America

171 Centre Scientifique de Monaco, 8 quai Antoine Ier, MC-98000, Monaco

172 Indian Institute of Technology Madras, Chennai 600036, India

173 Institut des Hautes Etudes Scientifiques, F-91440 Bures-sur-Yvette, France

174 IISER-Kolkata, Mohanpur, West Bengal 741252, India

175 Institut für Kernphysik, Theoriezentrum, 64289 Darmstadt, Germany

176 Whitman College, 345 Boyer Avenue, Walla Walla, WA 99362 United States of America

177 Université de Lyon, F-69361 Lyon, France

178 Hobart and William Smith Colleges, Geneva, NY 14456, United States of America

179 Dipartimento di Fisica, Università degli Studi di Torino, I-10125 Torino, Italy

180 University of Washington, Seattle, WA 98195, United States of America

181 INAF, Osservatorio Astronomico di Brera sede di Merate, I-23807 Merate, Lecco, Italy

182 Centro de Astrofísica e Gravitação (CENTRA), Departamento de Física, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal

183 Marquette University, 11420 W. Clybourn St., Milwaukee, WI 53233, United States of America

184 Indian Institute of Technology, Gandhinagar Ahmedabad Gujarat 382424, India

185 Université de Montréal/Polytechnique, Montreal, Quebec H3T 1J4, Canada

186 Indian Institute of Technology Hyderabad, Sangareddy, Khandi, Telangana 502285, India

187 INAF, Osservatorio di Astrofisica e Scienza dello Spazio, I-40129 Bologna, Italy

188 International Institute of Physics, Universidade Federal do Rio Grande do Norte, Natal RN 59078-970, Brazil

189 Villanova University, 800 Lancaster Ave, Villanova, PA 19085, United States of America

190 Andrews University, Berrien Springs, MI 49104, United States of America

191 Max Planck Institute for Gravitational Physik (Albert Einstein Institute), D-14476 Potsdam-Golm, Germany

192 Università di Siena, I-53100 Siena, Italy

193 Trinity University, San Antonio, TX 78212, United States of America

194 Van Swinderen Institute for Particle Physics and Gravity, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands

195 Department of Physics, University of Texas, Austin, TX 78712, United States of America

E-mail: lsc-spokesperson@ligo.org and virgo-spokesperson@ego-gw.it

Received 14 August 2019, revised 28 November 2019

<figure><img src="image_4.png" /></figure>

CrossMark
Accepted for publication 6 December 2019

Published 16 January 2020

**Abstract**

GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most $3.05 M_{\odot}$, and three equations of state considered here can be ruled out. We obtain a tighter limit of $2.67 M_{\odot}$ for the case that the merger results in a hypermassive neutron star.

Keywords: neutron stars, neutron star equation of state, gravitational wave astronomy, compact object mergers

(Some figures may appear in colour only in the online journal)

# 1. Introduction

On August 17, 2017, the Advanced LIGO [1] and Advanced Virgo [2] observatories detected a gravitational-wave signal, GW170817, from a low-mass compact binary system [3]. The signal was coincident with the gamma-ray burst GRB 170817A [4] and was followed by other observations spanning the electromagnetic spectrum [5]. These electromagnetic observations, coupled with measurements of the component masses from the gravitational-wave observation [6, 7], provide compelling evidence that GW170817 was produced from the merger of two neutron stars.

One of the key goals of observing such neutron star mergers is to extract information regarding the internal structure of neutron stars and to improve our understanding of their highly uncertain equation of state [8, 9]. Information regarding the equation of state is encoded in the variation of the orbital evolution induced by tidal interactions and the deformation of the neutron stars under their own angular momenta [10–14]. Such interactions are small at low frequencies but rapidly increase in strength in the final tens of gravitational-wave cycles before merger. Additional information about the equation of state is carried in the post-merger signal, as the remnant body can be a black hole, a stable neutron star or a meta-stable neutron star that collapses after some delay [15–21]. However, gravitational waves from the post-merger remnant are expected at higher frequencies where the detectors are less sensitive. No such post-merger signal was observed for GW170817 [6, 20, 22, 23].
The difficulty of computing the equation of state from first principles has led to a wide range of plausible theoretical models based on nuclear theoretical calculations. Astronomical observations of pulsars and x-ray binaries (e.g. [8, 24, 25]) have already been able to rule out a subset of proposed models, most notably by measurements of large neutron star masses [26]. Several studies have derived constraints on the equation of state through direct measurements of the neutron star mass and tidal deformability as inferred from GW170817 [27–41].

In this paper, we address whether we can rule out specific equations of state from the gravitational-wave data alone and whether we can rule out the possibility that one or both of the coalescing objects was a low-mass black hole. We opt to work with a subset of equations of state representing nuclear matter properties that vary over a wide range of possible values. The correlations between the properties of the neutron star with the underlying equation of state manifests itself in a wide range of macroscopic physical observables, such as the mass, radius or tidal deformability. In practice, we do not expect the true equation of state to be contained in any of the models used in this analysis. By comparing the Bayesian odds ratio for different models, we expect that equations of state that are most similar to the real equation of state will lead to the highest odds ratios whereas equations of state that differ significantly from the data will be down-ranked in the analysis, giving insight into the physical properties of the true equation of state.

We apply Bayesian techniques [42, 43] to compare a large set of equations of state against an alternative hypothesis. The first fully Bayesian study of the feasibility of probing the neutron star equation of state was presented in [44] and expanded upon in [45]. We consider two possible baseline hypotheses: the first hypothesis is that both component objects were low-mass black holes; the second hypothesis is that GW170817 was composed of two neutron stars with an equation of state corresponding to one of the best fit equations of state in our analysis. In order to make quantitative statements about the relative likelihoods of different equations of state, we do not allow the tidal parameters to vary independently and instead enforce that each neutron star obeys a specific equation of state, similar to [38, 39, 44–46]. This is in direct contrast to the analysis presented in [6], where the tidal parameters were allowed to vary independently implicitly allowing each neutron star to have a different equation of state. By performing parameter estimation using a fixed equation of state, we are able to conveniently incorporate a number of physical prior constraints, as was discussed in detail in [39].

We find that the equations of state yielding excessively large tidal deformabilities are disfavored by the data. For example, the MS1_PP model [47] is disfavored with a Bayes factor that is less than ~0.005 compared to the most favored equation of state models. However, we find that we cannot comprehensively rule out the majority of equations of state used in this analysis. Instead, we find that many of the equations of state have comparable evidences that are within an order of magnitude. This includes the models with the lowest tidal deformabilities as well as the binary black hole case.

In a previous analysis of GW170817 [39], we computed radii and tidal deformabilities for GW170817 using two methods. The first method links the tidal deformabilities of the two neutron stars by assuming that they obey the same equation of state [48, 49] and the second method utilized a parametrized equation of state [50, 51]. The neutron star radii were found to be no larger than ~13 km, disfavoring large radii and stiffer equations of state. In this paper, we perform a direct model comparison between equation of state models, finding results that are consistent with [39]. The equation of state models that are most preferred by the data in this analysis are the ones that predict radii close to those computed in [39].

In addition, we also explore the inferences that can be made about the fate of the remnant object for each proposed equation of state. The range of maximum neutron star masses for
the equation of state models allows remnants which are indefinitely stable neutron stars, long-lived supramassive neutron stars, or short-lived hypermassive neutron stars. Prompt collapse remains a possibility, though it can be ruled out for many models. In other cases, the required mass threshold from numerical relativity is not available or is insufficiently accurate to make a clear statement. For cases resulting in a supramassive neutron star, we use the mass range inferred for each of our equation of state models to derive the maximum rotation rate of the neutron star after it settles into a uniformly rotating state, obtaining an upper limit of 1.9 kHz.

Finally, we investigate constraints for the maximum mass of non-rotating neutron stars that follow from different assumptions about the type of the remnant. This is motivated by one of the possible explanations for the short gamma-ray burst models, which requires the formation of a black hole surrounded by a disk [52]. The other model assumes a rapidly rotating magnetar instead [53]. Which one is the correct description for the case at hand is an open question, but if we do assume the first one, then the presence of a black hole implies that neutron stars with the total baryonic mass inferred for the initial binary cannot be stable. The resulting limit is not very restrictive, but would rule out the three most extreme of our equations of state when assuming that the coalescing stars are slowly rotating. The bound can be further tightened by limiting the lifetime of the remnant to the observed delay of the short gamma-ray burst. The steep gradient of lifetime as function of mass near the maximum rotating neutron star mass then places the remnant mass in (or at least close to) the hypermassive range. Our limits are compatible with similar calculations carried out in [27, 29, 54], but are obtained using the fully consistent probability distributions computed for each equation of state model.

## 2. Methodology

### 2.1. Gravitational-wave Bayesian model comparison

In this work our goal is to perform a direct comparison of theoretical models for the neutron star equation of state against the gravitational-wave data for GW170817. In order to do so, we employ Bayesian model selection to test a specific equation of state against a competing hypothesis. We adopt both the low-mass binary black hole hypothesis and a best fit equation of state as baseline hypotheses and assess the impact this has on the inferences that can be made from the data. Detailed introductions to Bayesian methods in gravitational wave astronomy can be found in [42, 55–57].

Given the observed gravitational-wave data $s$ and any prior information $I$, the posterior odds in favor of a given model hypothesis $M_i$ can be expressed using the *odds ratio*

$$O_{12} = \frac{P(M_1|s, I)}{P(M_2|s, I)} = \frac{P(s|M_1, I) P(M_1|I)}{P(s|M_2, I) P(M_2|I)} = B_{12} \frac{P(M_1|I)}{P(M_2|I)}, \quad (1)$$

where $\mathcal{Z} = P(s|M_i, I)$ is the marginalized likelihood or *evidence* and $P(M_i|I)$ is the prior probability on the model hypothesis. If we assume that the two models are *a priori* equally likely, $P(M_1|I) = P(M_2|I)$, then the odds ratio $O_{12}$ reduces to the Bayes factor $B_{12}$. The evidence is calculated by marginalizing over all parameters weighted by their *prior probabilities* $P(\xi|M_i, I)$

$$\mathcal{Z} = P(s|M_i, I) \propto \int P(s|\xi, M_i, I) P(\xi|M_i, I) d\xi, \quad (2)$$

where the likelihood of obtaining the data realization $s$ given a gravitational-wave signal with physical parameters $\xi$ is
$$P(s|\xi, M_i, I) \propto e^{-\langle s-h(t,\xi) | s-h(t,\xi) \rangle / 2} \quad (3)$$

Together with a given set of prior probabilities, we have all the information needed to evaluate the evidence in equation (2). In the following subsections, we discuss the waveform models used to describe the gravitational-wave signal $h(t)$ and their dependence on the underlying physical parameters $\xi$.

## 2.2. Prior probabilities

An important consideration when computing the Bayes factors is the choice of prior probabilities. This choice will enter the analysis in two distinct ways. First, we must make a choice for the relative prior probabilities between any two models; $P(M_i|I)$. A wide range of equations of state are still consistent with observations [58], though additional theoretical considerations and observations could allow us to restrict the range of plausible equations of state [3, 6, 59]. Here, neglecting such additional information, we assume that the prior probability for all equations of state is the same. As detailed above, the odds ratio then reduces to the Bayes factor between any two models. The second choice of priors that we must make is the choice of prior probabilities on the parameters that enter a specific model, $P(\xi|M_i, I)$. As the choice of priors will affect the resulting Bayes factors, we explicitly detail the choices made in our analysis below.

The choice of priors that we use in this work is broadly similar to our previous analysis of GW170817 [3, 6]. We assume a priori that the source of GW170817 resides in NGC 4993, allowing us to adopt a fixed sky location. As in [6], although the distance to NGC 4993 is known, we infer the luminosity distance purely from the gravitational-wave data. The distance prior assumes that sources are uniformly distributed in volume up to a maximum distance of 75 Mpc. The coalescence time $t_c$ and phase $\phi_c$ are taken to be uniformly distributed and we assume that the binary orientation and polarization angle are isotropic.

The remaining parameters that we sample are the masses and spins of the two compact objects. The choice of prior ranges for these parameters is not so easily made. Observations of neutron stars in our galaxy can give us indications of the mass and spin distributions of binary neutron star systems [8, 60–66], but there is some uncertainty in these measurements. Additionally, these observations would not be appropriate to use if we want to assess the possibility that the system might be two black holes.

For these reasons we consider two sets of priors in this work. The first narrow prior is motivated by observed binary neutron star systems in our galaxy [26, 67]. For this prior, we assume that the prior on the component masses is Gaussian with a mean of $1.33M_\odot$ and a standard deviation of $0.09M_\odot$ [8]. The prior on the dimensionless spin magnitudes are taken to be uniformly distributed between 0 and 0.05 and the spin directions are assumed to be isotropic. This spin prior is consistent with the observed population of binary neutron stars that will merge within a Hubble time [3, 68]. As many of the waveform models we consider do not allow for spin components perpendicular to the orbital angular momentum, we choose to sample the spin direction and magnitude as described above and then ignore any spin component perpendicular to the orbital angular momentum. In practice, this amounts to choosing a non-uniform distribution on the spin magnitude parallel to the orbital plane, with a peak at 0.

The second broad prior allows for a wider range of masses and spins. Here we assume that component masses are uniformly distributed between $0.7M_\odot$ and $3.0M_\odot$ with the constraint that $m_2 \le m_1$. Priors on the chirp mass $M$ and mass ratio $q$ are determined from the concomitant Jacobian transformation. For computational efficiency, we impose a further cut on the chirp mass of the binary, such that we only consider points with chirp mass between 1.190 and
1.210 solar masses. However, this is wide enough to cover the full range of values supported by the data [3, 6, 38, 39]. In the broad prior, the spin magnitudes are taken to be uniformly distributed between 0 and 0.7. This is narrower than the spin prior used in [6] which extended to 0.89. The broad prior is intended to be a conservative spin prior compatible with the fastest-spinning known neutron star, which has a dimensionless spin $\le 0.4$ [69]. The choice of the limit 0.7 should not influence the results, as we obtain one-sided 95% upper limits for the dimensionless spin posterior which are much lower (below 0.35 for all models).

Non-rotating neutron stars cannot exceed a maximum mass which depends on the equation of state. If the sampler chooses a point in the parameter space where a given equation of state cannot support a non-rotating neutron star, we assume that the body is a black hole. Although rotation can increase the maximum mass, we do not allow for a scenario in which any neutron star is so massive that it requires rotation to prevent collapse. This choice is made to ensure that the prior volume is equal for all the equations of state. Our exact model assumption for a given equation of state is thus that objects below the maximum mass of non-rotating neutron stars are in fact neutron stars described by the given equation of state, otherwise light black holes. For the narrow prior, the distinction is irrelevant because we find no posterior support above the maximum mass for any equation of state. For the broad prior on the other hand, the posterior has significant support for mixed binaries for some equation of state models, as is further discussed in section 3.

There is also an upper bound on the dimensionless spin possible for uniformly rotating neutron stars, which depends on the equation of state. The maximum spin covered by the broad prior is below the maximum possible spin for some models considered here (compare section 2.4). Thus, the broad prior will contain a fraction of neutron stars rotating impossibly fast for some equations of state. In practice, this issue is not relevant as we find zero posterior support for such rotation rates.

Finally, we also present results where we compare the hypothesis that both bodies are low-mass black holes against the hypothesis that one, or both, bodies are neutron stars, but with a non-specific equation of state model. In this case we are sampling over the tidal deformability of the neutron stars (and deriving all other equation of state specific quantities from this value as explained below). For this analysis we use the same prior choices as above, specifically considering the broad mass prior, except for the following choices on spins and tidal deformability distributions. Here, a black hole is defined as having a spin magnitude uniformly distributed in [0, 0.89] and zero tidal deformability, while a neutron star has a spin magnitude uniformly distributed in [0, 0.05] and a tidal deformability uniformly distributed in [0, 5000]. The upper spin limit for black holes differs from the broad prior considered above. The specific models we study in this case are:

(i) **Binary Black Hole**: $\chi_1 \in [0, 0.89]$, $\chi_2 \in [0, 0.89]$, $\Lambda_1 = \Lambda_2 = 0$.

(ii) **Black Hole—Neutron Star**: $\chi_1 \in [0, 0.89]$, $\chi_2 \in [0, 0.05]$, $\Lambda_1 = 0$, $\Lambda_2 \in [0, 5000]$, $m_{BH} > m_{NS}$.

(iii) **Neutron Star—Black Hole**: $\chi_1 \in [0, 0.05]$, $\chi_2 \in [0, 0.89]$, $\Lambda_1 \in [0, 5000]$, $\Lambda_2 = 0$, $m_{NS} > m_{BH}$.

(iv) **Binary Neutron Star, no assumptions on equation of state**: $\chi_1 \in [0, 0.05]$, $\chi_2 \in [0, 0.05]$, $\Lambda_1 \in [0, 5000]$, $\Lambda_2 \in [0, 5000]$.

(v) **Binary Neutron Star, neutron stars obey same equation of state**: $\chi_1 \in [0, 0.05]$, $\chi_2 \in [0, 0.05]$, $\Lambda_2 \equiv (\Lambda_2 + \Lambda_1)/2 \in [0, 5000]$, $\Lambda_a \equiv (\Lambda_2 - \Lambda_1)/2 = \Lambda_a(\Lambda_s, m_2/m_1)$.

where here $\chi_i$ and $\Lambda_i$ denote the dimensionless spin and tidal deformability of the body $i$ respectively. The analysis of model (iv) allows the two component tidal deformabilities to vary independent of each other. In contrast, the analysis of model (v) uses an equation of
state-independent relation between the two tidal deformabilities given the mass ratio of the system [48, 49].

## 2.3. Waveform models for the gravitational-wave signal

Gravitational-wave signals from binary neutron stars are distinct from signals from binary black holes. The deformation of the neutron star due to tidal effects and a spin-induced quadrupole, as well as the post-merger behavior of the remnant, will all leave an imprint in the emitted gravitational waves. Unfortunately, no waveform model yet exists that can consistently model all three of these effects, as modeling the post-merger behavior requires expensive numerical simulations. Such simulations have been performed by a number of groups in recent years, enabling an exploration of the expected features of post-merger emission [70–72]. Though this work has not yet reached the stage where one can predict waveforms with arbitrary masses, spins and equation of state, unmodeled analyses can target the expected post-merger signal [73–75]. Since the post-merger component of the signal occurs at frequencies ~2000 Hz, where Advanced LIGO and Advanced Virgo sensitivity is decreasing, observing this regime was not possible for GW170817 [3, 6, 22, 76]. For the present analysis, it is therefore safe to ignore the post-merger signal.

We use waveform models that predict gravitational-wave signals for a large range of parameters and include both the effects of the tidal deformation and of the spin-induced quadrupole. In contrast to previous analyses [3, 6], where the individual tidal deformabilities $\Lambda$ of the neutron stars were treated as independent intrinsic parameters, we assume a specific equation of state model and the tidal deformability is computed as a function of the masses. To construct the functions $\Lambda(m)$, the Tolman–Oppenheimer–Volkoff differential equations [77, 78] for the stellar structure are solved numerically together with those for the tidal deformability [12, 79].

The quadrupole-monopole effects are caused by the quadrupolar deformations induced by the spin [11], and are quadratic in the spins. These effects are parametrized by a quadrupole-monopole parameter $\kappa$. Similarly to the tidal deformability, this parameter is treated as a pure function of the masses $\kappa(m)$, determined in practice from the value of $\Lambda$ via fits of approximate universal relations for neutron stars [48].

The waveform models used in this study are essentially the same as in the previous analyses from the LIGO–Virgo collaboration [3, 6, 39], to which we refer the reader for more details. The TaylorF2 model uses purely post-Newtonian information to generate closed-form waveforms in the Fourier domain, only allowing for spins aligned with the orbital angular momentum. It includes point-particle and spin contributions (see, e.g. [80]) and tidal contributions. In contrast to [6], we consider several variants of TaylorF2, incorporating tidal terms up to the 6, 7 or 7.5 post-Newtonian order [13, 14, 81]<fnref n="196" />.

A closed-form tidal approximant, denoted NRTidal, calibrated to binary neutron star numerical relativity simulations was presented in [82, 83]. In this paper, EOB+NRT and Phenom+NRT will denote the waveform models constructed from adding this numerical relativity tuned approximant for tidal effects onto the underlying point-particle waveform models, SEOBNRv4_ROM [84] and IMRPhenomPv2 [85–88] respectively. While the EOB+NRT model is restricted to aligned spins, Phenom+NRT incorporates an effective representation of precession effects [85, 86]. All three models TaylorF2, Phenom+NRT and EOB+NRT

<footnote n="196">¹⁹⁶ The tidal terms at the 7PN and 7.5PN orders are not strictly complete, with a missing term at 7PN (2PN relative order in the tidal terms) but this missing piece is argued in [14] to be unimportant quantitatively.</footnote>
include post-Newtonian quadrupole-monopole terms, up to the third PN order in the phasing
[11, 89].

In addition, we use two time-domain effective one body models that incorporate tidal effects, SEOBNRv4T [90, 91] and TEOBResumS [92]. SEOBNRv4T includes dynamical tides and the effects of the spin-induced quadrupole moment. TEOBResumS incorporates a gravitational-self-force re-summed tidal potential and the spin-induced quadrupole moment. Both models have been shown to be compatible with state-of-the-art numerical simulations of binary neutron stars up to merger [92, 93]. Also compare [94, 95] regarding recent extensions.

Lacking a prescription for the post-merger signal, all waveforms have an equation of state dependent termination frequency that is close to the merger frequency. For TaylorF2, the termination frequency is determined according to equation (13) of [45]<fnref n="197" />. For EOB+NRT and Phenom+NRT, the termination frequency is computed following equation (11) of [83]. Thus the NRT waveforms with zero tidal deformabilities are not strictly speaking equivalent to the binary black hole waveforms: in the NRT prescription the merger is cut instead of terminating with a ringdown signal as would be the case in the underlying point-particle model. However, this difference only concerns high frequencies and should not affect our evidence computations. For most of the configurations considered in this analysis, a remnant black hole would have a ringdown frequency above ~4 kHz, justifying the above statement.

It is important to note that all our waveform models are imperfect and come with systematic errors. Without a quantitative model to represent this as a modeling uncertainty, we cannot marginalize over the waveform error in our analysis. Comparisons between different waveform models [83, 96, 97] indicate that TaylorF2 can produce small differences from Phenom+NRT and EOB+NRT in recovering the tidal parameters, while yielding comparable results for the masses and spins. This was indeed observed in the analysis of GW170817 [3, 6, 38, 39], where the TaylorF2 model yielded larger upper limits for the tidal parameters. In order to obtain a rough estimate how differences between waveform models affect the Bayes factor calculations, we compare results obtained with different waveform families.

## 2.4. Equations of state

For this work, we have collected 24 different prescriptions describing neutron star equa-
tions of state from the literature, computed by various nuclear physics approximations. The
full list together with the relevant citations is given in appendix A.

All equations of state describe zero-temperature matter in β-equilibrium. We require equa-
tions of state to respect causality, i.e. the maximum sound speed inside stable neutron stars
may not exceed the speed of light at any mass up to the maximum allowed. For comparison
with previous studies [6, 39], however, we also include some results for the WFF1 [98] equa-
tion of state model, which violates causality before reaching the central density of the max-
imum-mass neutron star. One equation of state model (HQC18 [9]) incorporates a transition
between a hadronic phase in the crust and a quark matter phase in the core. The other equa-
tion of state models describe purely hadronic matter without significant phase transitions. We
deliberately do not use any constraints from previous studies of the detected event [3, 6, 38,
39] to further narrow the equation of state model selection.

The equation of state properties most important for our analysis of the gravitational-
wave signal are the maximum gravitational mass of non-rotating neutron stars and their

<footnote n="197">¹⁹⁷ There is a typo in equation (14) of [45], the denominator should be the sum of the two radii cubed, without which the equation is not dimensionally correct.</footnote>
tidal deformability as a function of gravitational mass. In order to derive implications for the merger remnant, we also require the maximum mass of uniformly rotating neutron stars, as well as the maximum rotation rate and moment of inertia as a function of mass. The maximum dimensionless spin for a given mass is also relevant with regard to the spin prior. We compute all those properties for all equations of state used in this work. Key quantities are tabulated in appendix A.

The properties of non-rotating neutron stars which follow from the equations of state used here are shown in figure 1. For comparison, we also show the credible intervals for the comp- onent gravitational masses derived in [6] for different priors on the spin. The maximum gravi- tational mass covers a range 1.92 M$_{\odot}$ (BHF_BBB2) to 2.75 M$_{\odot}$ (MS1_PP, MS1B_PP). The discovery of neutron star PSR J0348+0432 [26] with mass 2.01 ± 0.04 M$_{\odot}$, and the very recent observation of PSR J0740+6620 [99] with mass of 2.14$_{-0.09}^{+0.10}$ M$_{\odot}$, provide a lower limit for the maximum mass of neutron stars. However, we also included some models with maximum mass slightly below those constraints in our comparison. Although causality would allow maximum masses even larger than our most massive models (see [100, 101]), the cho- sen range covers most nuclear physics estimates. The smallest tidal deformability within the mass credible interval [6] for the low spin prior is found for the WFF1 equation of state model, with $\Lambda = 61$, and the largest one for the MS1_PP equation of state model, with $\Lambda = 3728$. In the same mass range, we find that the dimensionless spin at the mass-shedding limit is between 0.66 (KDE0V1) and 0.76 (SKI5).

## 2.5. Gravitational-wave data

For our analysis, we follow [7] and use data from the two LIGO observatories with the full frequency-dependent calibration uncertainties described in [102, 103] and following the sub- traction of independently measured noise sources as discussed in [104]. As discussed in [3], a non-Gaussian transient occurs in the data recorded by the LIGO Livingston observatory. This non-Gaussian transient is removed from the data as in [3] using the methods discussed in [73]. We have checked that the instrumental transient subtraction does not bias the results of our inference substantially [105]. We also follow [7] in using data from the Virgo observatory using the calibration and noise-subtraction model discussed in [106]. For both the LIGO and Virgo observatories we analyze data between 23 Hz and 2048 Hz.

## 2.6. Evaluating the Bayes factors

In this work we utilize two different strategies to compute the Bayes factors. We use pre-exist- ing tools in the LALINFERENCE [42] and RIFT [43] parameter inference codes to compute the evidence for each equation of state model (equation (2)). For LALINFERENCE, this involves performing a distinct detailed Bayesian analysis of the data for each proposed model, using nested sampling [107] and thermodynamic integration [108] to compute the evidence [42]. Both methods have been extensively tested [42] and produce consistent results [109].

For RIFT, this calculation is organized as a two-stage process. The first stage estimates a single marginal likelihood $P(s|x, \Lambda_1, \Lambda_2, I) = \int d\theta P(s|\theta, x, \Lambda_1, \Lambda_2; M_i, I)P(\theta)$ as a function of the extrinsic parameters $\theta$ and the intrinsic binary parameters $x = (m_1, m_2, \chi_1, \chi_2)$ and $(\Lambda_1, \Lambda_2)$ without assuming an equation of state. The next stage computes the evidence via $P(s|M_i, I) = \int dx P(x)P(s|x, \Lambda(m_1|M_i), \Lambda(m_2|M_i); M_i, I)$, where $\Lambda(m|M_i)$ denotes the deform- ability as function of mass assuming the equation of state model $M_i$. Using a precomputed marginal likelihood, this second stage can be performed very rapidly for any equation of state,
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Properties of non-rotating neutron stars for all equations of state considered in this work. Top: dimensionless tidal deformability versus gravitational mass. Bottom: proper circumferential radius. All sequences terminate at the maximum mass model (except WFF1, which is cut where the maximum sound speed reaches the speed of light). The legend entries are sorted by the radius of a $1.36 M_{\odot}$ neutron star. For comparison, we show the 90% symmetric credible interval given in [6] for the component masses, with dashed vertical lines for the low-spin prior, dotted lines for the high spin prior, and the solid line for the equal-mass limit $1.36 M_{\odot}$. The shaded vertical area marks the measured mass for pulsar PSR J0348+0432 [26].</figcaption></figure>

in seconds to minutes depending on the accuracy required. Unless otherwise noted, all RIFT calculations are performed using the marginalized likelihoods $P(s|x, \Lambda_1, \Lambda_2)$ employed in [7].

## 3. Results

For each equation of state described in section 2.4, the Bayes factors are computed using either LALINFERENCE or RIFT, using both the narrow and broad priors discussed in section 2.6. The Bayes factors for all waveform approximants and the statistical uncertainties are given in appendix B. The Bayes factors calculated adopting the narrow prior are shown in
<figure><img src="image_6.png" /><figcaption>**Figure 2.** The Bayes factors for the narrow prior results using different waveform approximants, as given in table B1. The equations of state have been sorted by the compactness $C$ of the neutron star at a fiducial mass of $1.36M_{\odot}$. The results show a clear preference for more compact equations of state. We adopt the guidelines of [111] for the interpretation of the Bayes factors. Here $\mathcal{Z} = P(s|M_i, I)$ is the evidence of a given equation of state or the fiducial binary black hole model. Statistical errors on the Bayes factors from the nested sampling algorithm are negligible.</figcaption></figure>

tables B1 and B3 for LALINFERENCE and RIFT, respectively. Similarly, for the broad prior, the LALINFERENCE results are given in table B2 and the RIFT results in table B4. The Bayes factors for the LALINFERENCE runs are quoted with respect to the low-mass binary black hole hypothesis. The RIFT results re-use the analysis presented in [7] and are not normalized with respect to the BBH case due to the way in which the fits are performed. Instead, the RIFT Bayes factors are quoted relative to the arbitrary hypothesis that GW170817 was a binary neutron star with the SLY9 equation of state, which has modest support and shows relatively little dispersion between waveform approximants in the LALINFERENCE results.

For the narrow prior, the gravitational-wave data *does not favor and nor can it rule out* any single equation of state with respect to the low-mass binary black hole hypothesis. This can be seen in figure 2, although the data clearly favors softer equations of state, which predict more compact neutron stars with lower tidal deformabilities. The broad prior case, shown in figure 3, exhibits the same trend. The relation is visualized more directly in figure 4, showing fiducial neutron star properties as a function of the Bayes factor with respect to the KDE0V model, which is the most favored for the narrow prior with TaylorF2 (7.5pN) approximant. The models with largest fiducial radius and tidal deformability are *very strongly disfavored* (MS1_PP) or at least *strongly disfavored* (MS1B_PP) compared to the numerous models with small fiducial radius and deformability, such as KDE0V.

Figure 4 also shows upper limits for radius and tidal deformability given by parameter estimation studies [6, 38, 39]. For all our equations of state, and mass ratios in the range $q > 0.7$, the effective tidal deformability and the tidal deformability at the fiducial mass agree better than 6%. To this accuracy, we can compare the upper limits to the values plotted for the fiducial mass models. We find that some of our models with tidal deformabilities above the limits are not strongly disfavored with respect to the most favored equations of state.

Although we can only compute Bayes factors between the models at hand, the Bayes factors exhibit a trend of decreasing with an increasing fiducial radius and tidal deformability. If this trend holds in general, then for a fiducial mass of $1.36M_{\odot}$ models with tidal deformability $\Lambda(M = 1.36) \gtrsim 1600$ or radius $R(M = 1.36) \gtrsim 14.9$ km are very strongly disfavored.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** As per figure 2 but for the broad prior using different waveform approximants, given in table B2.</figcaption></figure>

For the broad prior, the gravitational-wave data alone does not prefer any single equation of state, as was the case for the narrow prior. However, MS1_PP and MS1B_PP are now *strongly disfavored* against the low-mass binary black hole hypothesis, as shown in figure 3. Normalizing the Bayes factors against KDE0V, we find that MS1_PP is *very strongly* and MS1B_PP *strongly* disfavored with respect to the softer KDE0V equation of state, as shown in figure 4.

In contrast to the narrow prior, we find that for some equations of state there is significant posterior support for a mixed binary where the heavier object is a light black hole. The reason for this is that the mass range is less constrained and exceeds the maximum possible mass of non-rotating neutron stars for those equations of state. Models BHF_BBB2, KDE0V1, HQC18, and KDE0V have 13, 7, 5, and 7 percent posterior probability, respectively, for a mixed binary. The remaining models have support below 5% for mixed binary, in particular MS1_PP and MS1B_PP have zero support for mixed binaries. For the broad prior and the WFF1 equation of state, there is a ≈14% posterior probability for a central neutron star density in an unphysical (causality violating) region of the parameter space.

For both sets of priors, we find that the bulk of equations of state considered in this paper produce evidences that are consistent to within an order of magnitude. This is in line with the expectations for second-generation gravitational-wave detectors, where ~*O*(20) detections will allow hypothesis ranking to distinguish between hard, moderate and soft equations of state [44]. Even if the true equation of state were to be one of those considered in this paper, we should not necessarily expect this to be the most preferred equation of state in the hypothesis ranking [45]. For example, this could arise due to the effects of noise or other systematics. More realistically, it is also possible that the true equation of state is not contained in the finite set of models considered here, though we expect the highest ranked hypotheses to share similar macroscopic features to the true neutron star equation of state.

It is important to assess the systematic bias inherent in our waveform models. Such systematic biases could arise from a different treatment of the underlying point particle limit, which can start to become significant at higher frequencies where tidal effects also become more prominent. Additionally, the different waveform models employ different prescriptions for the tidal contribution to the phasing, which can lead to discrepancies in the recovery of tidal parameters. As demonstrated in previous studies [39, 83, 96, 97, 110], the tidal effect calculated in the TaylorF2 waveform models is typically smaller than that predicted using the
<figure><img src="image_8.png" /></figure>

**Figure 4.** The fiducial tidal deformability $\Lambda$ at $1.36M_{\odot}$ (top) and the fiducial neutron star radius $R$ at $1.36M_{\odot}$ (bottom) plotted as a function of the Bayes factors for the narrow prior (left) and broad prior (right). For the binary black hole case, the plotted radius is the Schwarzschild radius. The Bayes factors have been normalized with respect to KDE0V, which is the most favored equation of state for the narrow prior with TaylorF2 (7.5pN) waveform model. Compared to the KDE0V model, MS1_PP is very strongly disfavored (ln BF $\le -5$) and MS1B_PP is strongly disfavored (ln BF $\le -3$) for both priors and all waveform approximants. For comparison, we plot limits from the literature. The solid green horizontal lines mark the upper limit for the coalescing neutron star radii from [39]. Red dashed lines mark the upper limit given for the radius given in [38] (assuming a common radius for the two stars) and the one-sided upper limit for the effective tidal deformability from [38]. The blue dash-dotted lines mark the upper limit (highest posterior density interval) for the effective tidal deformability from [6].

NRTidal prescription. As such, the tidal parameters inferred from the data using the TaylorF2 waveform model will typically be larger in order to compensate for the reduced tidal effect. Consequentially, we find that the Bayes factors calculated using the TaylorF2 model tend to show greater support for stiffer equations of state, which yield larger tidal deformabilities, in comparison to the effective-one-body or phenomenological waveform models. However, as can be seen in figures 2–4, the Bayes factors between the different waveform approximants still show reasonably good agreement.

Finally, we also ask if we can rule out the binary black hole hypothesis by computing Bayes factors for the binary neutron star hypothesis, or the hypothesis that there was only one neutron star in the binary, without assuming a specific equation of state. We consider both
<table><caption>**Table 1.** Bayes factors between various hypotheses on the nature of the GW170817 components obtained using different waveform approximants and two methods for computing the evidence. All Bayes factors are given in relation to the 'binary neutron star with independent equations of state' hypothesis. We do not present results for Phenom+NRT/Nest when a BH is part of the binary, as the combination of precession and high spins makes this analysis computationally very expensive. We find no clear preference for any hypothesis. The table abbreviations stand for nested sampling (Nest), Markov-chain Monte-Carlo/thermodynamic integration (MCMC), binary black hole (BBH), mixed binary where the neutron star is heavier (NSBH), mixed binary where the black hole is heavier (BHNS), binary neutron star (BNS), equation of state (EoS).</caption><thead><tr><th>Waveform/Method</th><th>BBH</th><th>NSBH</th><th>BHNS</th><th>BNS (common EoS)</th></tr></thead><tbody><tr><td>Phenom+NRT/Nest</td><td>—</td><td>—</td><td>—</td><td>5.0</td></tr><tr><td>Phenom+NRT/MCMC</td><td>1.1</td><td>0.7</td><td>1.7</td><td>9.6</td></tr><tr><td>EOB+NRT/MCMC</td><td>1.7</td><td>1.6</td><td>4.0</td><td>5.5</td></tr><tr><td>TaylorF2 (7pN)/Nest</td><td>1.9</td><td>1.4</td><td>3.7</td><td>4.9</td></tr><tr><td>TaylorF2 (7pN)/MCMC</td><td>1.0</td><td>1.2</td><td>3.4</td><td>4.4</td></tr><tr><td>TaylorF2 (7.5pN)/Nest</td><td>1.6</td><td>1.3</td><td>3.2</td><td>4.4</td></tr><tr><td>TaylorF2 (7.5pN)/MCMC</td><td>0.9</td><td>1.0</td><td>5.5</td><td>5.0</td></tr></tbody></table>

the hypotheses that the single neutron star is heavier and lighter than its potential black hole companion. The specific priors that we use in these cases were discussed earlier in section 2.2. Table 1 presents the Bayes factors compared to the 'binary neutron star with independent equation of state' hypothesis. We present results with different waveform approximants and methodologies to compute the evidence of each model.

From the results in table 1 we can draw a number of conclusions. First, we find that the computed Bayes factors agree closely between different waveform approximants, indicating that uncertainties in the waveform modeling are not compromising our results. Second, we find that both nested sampling and thermodynamic integration return consistent results for the Bayes factors, giving us confidence on the validity of the calculation. Third, and most interesting from a physical perspective, we find that the GW170817 data alone does not strongly favor one hypothesis over the other. All of the Bayes factors we compute are in the range of 1 to 5, suggesting that the gravitational-wave data alone do not show a strong preference for any of the hypotheses considered here.

It may seem strange that the binary neutron star model, assuming no common equation of state, is the least favored hypothesis. Bayes factors depend on the Occam penalty, i.e. how well a model fits the given data (goodness-of-fit) and how many constrained parameters a model includes (dimensionality). We expect all models to be preferred over the binary neutron star model, assuming a common equation of state, if they fit the data equally well. The binary black hole model has two fewer parameters and the neutron star-black hole hypotheses have one fewer parameter. A brief back-of-the-envelope calculation can be used to support the Bayes factors seen in table 1. The Occam penalty is given by $\delta\theta/\Delta\theta$, the posterior width over the prior width, and penalizes a model for each additional constrained parameter it employs, since $\delta\theta < \Delta\theta$. From section 2.2, we see that the width of the prior on the tidal parameters is $\Delta\Lambda \sim 5000$ whereas the width on the posterior is on the order of $\delta\Lambda \sim 1000$ [6]. From this, we expect the binary neutron star hypothesis with a common equation of state to be preferred over an independent equation of state hypothesis by a factor $\sim 5$, as it uses one fewer parameter. This is in agreement with the results found in table 1 and as demonstrated in [49] for simulated binary neutron star signals.
# 4. Implications

After comparing the relative probabilities of the different equation of state models, we now discuss the implications that follow from assuming any of the equations of state are correct. We also discuss constraints for the maximum mass of neutron stars that follow from assumptions on the fate of the remnant, but without assuming any equation of state.

## 4.1. Baryonic mass and central density

The total baryonic mass (defined here as baryon number times a fiducial mass constant $1.66 \times 10^{-27}$ kg) is the most important quantity for the evolution during the post-merger phase, including the lifetime of the remnant and the mass of the accretion disk, which are relevant for the short gamma-ray burst and the kilonova. The maximum density inside the initial neutron stars is important to interpret our results properly. The relative probabilities we compute for different equation of state models only signify that in the density range occurring in the coalescing neutron stars, matter is more likely described by one equation of state model than another.

To compute the baryonic mass, we start with the posterior samples obtained assuming a given equation of state model. From the gravitational masses and spins in a given sample, we obtain the baryonic masses of the neutron stars. Those are computed using interpolation tables we build for sequences of non-rotating and maximally rotating neutron stars for each equation of state model. To account for the spin, we expand the gravitational mass at given baryonic mass to second order in the squared dimensionless spin parameter. The zeroth- and first-order coefficients are determined from the binding energy and moment of inertia of the non-rotating model. The second-order coefficient is then determined from the angular momentum and gravitational mass of the maximally rotating model. The spin corrections are relatively small compared to the statistical error both for narrow and broad priors. To obtain limits for the energy densities that can be reached before merger, we use the central densities of non-rotating neutron stars with baryonic masses computed as described above.

The results for the narrow prior are given in table 2. We find that the pre-merger gravitational-wave signal does not carry information about matter at energy densities above $0.5-1.3 \times 10^{15}$ g cm$^{-3}$, assuming the true equation of state is similar to one of the models considered here. This agrees well with the results by [39] (figure 2 therein) obtained using parametrized equations of state, and also with results in [40] based on Gaussian process models for the equation of state.

## 4.2. Fate of merger remnant

Next, we assess the fate of the merger remnant, which could either form a black hole directly at merger (prompt collapse), a short-lived hypermassive neutron star, a long-lived supramassive neutron star, or even an unconditionally stable neutron star. This is particularly relevant with regard to the observed short gamma-ray burst and kilonova counterparts. One possible engine for short gamma-ray bursts is given by a black hole submerged in a massive disk [52]. Numerical simulations suggest that prompt black hole formation at merger might be in tension with the production of a short gamma-ray bursts [29, 112], favoring delayed collapse. The other model for the engine consists of a rapidly rotating and strongly magnetized neutron star (magnetar scenario [53, 113, 114]) embedded in a disk. In this model, the remnant needs to survive at least long enough to initiate a jet (however, there can be a delay until the jet produces the observed short gamma-ray burst). The fate of the remnant is also tightly related to the ejection
<table><caption>**Table 2.** Implications derived from the equation of state model specific parameter estimation runs with narrow prior and TaylorF2 (7.5 PN) waveform model. $M_B^0$ is the total initial baryonic mass. We provide the median value and the 90% two-sided credible bounds. $E_c$ is the 90% one-sided upper limit for the central energy density of the heavier neutron star. The normalized mass $\tilde{m}_0$ is simply a linear transform of $M_B^0$ defined by equation (4), such that the supramassive mass range for the given equation of state model corresponds to $\tilde{m} \in (0, 1)$. The potential impact of mass ejection is given in terms of $\Delta_{0.1} = \tilde{m}_0 - \tilde{m}_R$, computed from equation (5) for a fiducial ejecta mass $0.1 M_\odot$ (the given value is the 90% one-sided upper limit). The ‘prompt collapse’ column denotes whether a black hole can be formed directly during the merger, based on available estimates of the corresponding mass threshold, and assuming prompt collapse only occurs for hypermassive systems. Hypermassive cases lacking data on the prompt collapse threshold are denoted with a dash, while ‘possibly’ refers to the case where the available bounds on the threshold are compatible with both outcomes. $I_{max}$ and $F_{max}$ are upper limits (see main text) for moment of inertia and rotation rate for a uniformly rotating remnant. $I_{max}$ is given in geometric units of $4.335\,871 \times 10^{43} \text{ g cm}^2$.</caption><thead><tr><th>Equation of state</th><th>$M_B^0(M_\odot)$</th><th>$E_c/c^2$ (10<sup>15</sup> g cm<sup>-3</sup>)</th><th>$\tilde{m}_0$</th><th>$\Delta_{0.1}$</th><th>Prompt collapse</th><th>$I_{max}(M_\odot^3)$</th><th>$F_{max}$(kHz)</th></tr></thead><tbody><tr><td>APR4_EPP</td><td>3.021<sup>+0.019</sup><sub>-0.009</sub></td><td>1.05</td><td>0.883<sup>+0.041</sup><sub>-0.019</sub></td><td>0.22</td><td>No</td><td>84</td><td>1.79</td></tr><tr><td>MS1_PP</td><td>2.940<sup>+0.014</sup><sub>-0.008</sub></td><td>0.47</td><td>-0.608<sup>+0.023</sup><sub>-0.013</sub></td><td>0.17</td><td>No</td><td>153</td><td>1.00</td></tr><tr><td>MS1B_PP</td><td>2.949<sup>+0.014</sup><sub>-0.008</sub></td><td>0.49</td><td>-0.581<sup>+0.023</sup><sub>-0.013</sub></td><td>0.16</td><td>No</td><td>152</td><td>1.03</td></tr><tr><td>KDE0V</td><td>3.013<sup>+0.019</sup><sub>-0.009</sub></td><td>1.17</td><td>1.859<sup>+0.050</sup><sub>-0.025</sub></td><td>0.27</td><td>—</td><td>66</td><td>1.86</td></tr><tr><td>SKI6</td><td>2.994<sup>+0.016</sup><sub>-0.009</sub></td><td>0.82</td><td>0.848<sup>+0.036</sup><sub>-0.018</sub></td><td>0.22</td><td>No</td><td>97</td><td>1.53</td></tr><tr><td>H4</td><td>2.978<sup>+0.015</sup><sub>-0.009</sub></td><td>0.65</td><td>1.385<sup>+0.035</sup><sub>-0.017</sub></td><td>0.22</td><td>No</td><td>104</td><td>1.36</td></tr><tr><td>MPA1</td><td>3.004<sup>+0.016</sup><sub>-0.009</sub></td><td>0.74</td><td>-0.016<sup>+0.029</sup><sub>-0.016</sub></td><td>0.18</td><td>No</td><td>114</td><td>1.34</td></tr><tr><td>SK255</td><td>2.968<sup>+0.015</sup><sub>-0.008</sub></td><td>0.80</td><td>1.120<sup>+0.038</sup><sub>-0.020</sub></td><td>0.24</td><td>—</td><td>91</td><td>1.61</td></tr><tr><td>SKI2</td><td>2.968<sup>+0.016</sup><sub>-0.008</sub></td><td>0.72</td><td>1.025<sup>+0.036</sup><sub>-0.019</sub></td><td>0.23</td><td>—</td><td>100</td><td>1.53</td></tr><tr><td>SKI4</td><td>2.997<sup>+0.017</sup><sub>-0.008</sub></td><td>0.83</td><td>0.911<sup>+0.035</sup><sub>-0.018</sub></td><td>0.22</td><td>No</td><td>94</td><td>1.58</td></tr><tr><td>SKOP</td><td>2.994<sup>+0.016</sup><sub>-0.009</sub></td><td>1.02</td><td>1.808<sup>+0.043</sup><sub>-0.022</sub></td><td>0.26</td><td>—</td><td>73</td><td>1.70</td></tr><tr><td>SLY2</td><td>3.005<sup>+0.018</sup><sub>-0.009</sub></td><td>1.02</td><td>1.401<sup>+0.043</sup><sub>-0.022</sub></td><td>0.25</td><td>—</td><td>77</td><td>1.79</td></tr><tr><td>SLY9</td><td>2.988<sup>+0.016</sup><sub>-0.008</sub></td><td>0.86</td><td>0.999<sup>+0.037</sup><sub>-0.019</sub></td><td>0.23</td><td>—</td><td>91</td><td>1.68</td></tr><tr><td>BHF_BBB2</td><td>3.021<sup>+0.020</sup><sub>-0.009</sub></td><td>1.26</td><td>2.032<sup>+0.053</sup><sub>-0.025</sub></td><td>0.27</td><td>—</td><td>63</td><td>1.87</td></tr><tr><td>HQC18</td><td>3.014<sup>+0.018</sup><sub>-0.010</sub></td><td>1.03</td><td>1.305<sup>+0.040</sup><sub>-0.020</sub></td><td>0.22</td><td>—</td><td>81</td><td>1.72</td></tr><tr><td>KDE0V1</td><td>3.005<sup>+0.018</sup><sub>-0.009</sub></td><td>1.13</td><td>1.837<sup>+0.048</sup><sub>-0.024</sub></td><td>0.27</td><td>—</td><td>68</td><td>1.80</td></tr><tr><td>RS</td><td>2.980<sup>+0.016</sup><sub>-0.008</sub></td><td>0.80</td><td>1.174<sup>+0.037</sup><sub>-0.019</sub></td><td>0.24</td><td>—</td><td>92</td><td>1.61</td></tr><tr><td>SK272</td><td>2.965<sup>+0.016</sup><sub>-0.008</sub></td><td>0.75</td><td>0.780<sup>+0.036</sup><sub>-0.019</sub></td><td>0.23</td><td>No</td><td>100</td><td>1.43</td></tr><tr><td>SKI3</td><td>2.967<sup>+0.015</sup><sub>-0.008</sub></td><td>0.68</td><td>0.744<sup>+0.034</sup><sub>-0.018</sub></td><td>0.22</td><td>No</td><td>107</td><td>1.35</td></tr><tr><td>SKI5</td><td>2.957<sup>+0.015</sup><sub>-0.008</sub></td><td>0.63</td><td>0.775<sup>+0.033</sup><sub>-0.019</sub></td><td>0.23</td><td>No</td><td>109</td><td>1.32</td></tr><tr><td>SKMP</td><td>2.990<sup>+0.016</sup><sub>-0.008</sub></td><td>0.86</td><td>1.186<sup>+0.037</sup><sub>-0.019</sub></td><td>0.23</td><td>—</td><td>88</td><td>1.68</td></tr><tr><td>SLY</td><td>3.005<sup>+0.017</sup><sub>-0.009</sub></td><td>1.04</td><td>1.425<sup>+0.042</sup><sub>-0.021</sub></td><td>0.25</td><td>Possibly</td><td>76</td><td>1.82</td></tr><tr><td>SLY230A</td><td>3.007<sup>+0.017</sup><sub>-0.009</sub></td><td>0.97</td><td>1.175<sup>+0.040</sup><sub>-0.020</sub></td><td>0.23</td><td>—</td><td>83</td><td>1.76</td></tr></tbody></table>

of matter and the spectral evolution of the resulting kilonova [5, 115–118]. For example, the amount, composition and velocity of ejected matter required for the observed kilonova might be in tension with expectations for the prompt collapse case (for a recent review on multi-messenger constraints, see [119]). However, the theoretical modeling of the kilonova is a very complex and rapidly evolving field of research (e.g. [120–123] and the references therein).
The most relevant quantity regarding the remnant type is the baryonic mass $M_B^R$ of the system, in comparison to the maximum baryonic mass of non-rotating neutron stars, $\hat{M}_B^{TOV}$, and the maximum baryonic mass of neutron stars uniformly rotating at the mass-shedding limit, $\hat{M}_B^{ROT}$. Since the baryon number is conserved, the remnant mass is given by $M_B^R = M_B^0 - M_B^E$, where $M_B^E$ denotes the baryonic mass expelled from the system dynamically and by winds, and $M_B^0$ the total initial baryonic mass.

For discussing the expected remnant evolution, we introduce a useful dimensionless measure $\tilde{m}$ defined as

$$ \tilde{m}(M_B) = k_m \left( \frac{M_B}{\hat{M}_B^{TOV}} - 1 \right), \quad k_m = \frac{\hat{M}_B^{TOV}}{\hat{M}_B^{ROT} - \hat{M}_B^{TOV}}. \qquad (4) $$

The maximum masses $\hat{M}_B^{ROT}$ and $\hat{M}_B^{TOV}$, which we compute for each equation of state model are given in table A1. Based on $\tilde{m}_R \equiv \tilde{m}(M_B^R)$ for a given equation of state model, we classify the remnant as hypermassive if $\tilde{m}_R > 1$, supramassive if $0 \le \tilde{m}_R \le 1$, and stable if $\tilde{m}_R < 0$.

The quantity $\tilde{m}_R$ serves as a rough indication for the remnant lifetime $\tau_R$. A stable remnant will never form a black hole, a hypermassive remnant either forms a black hole directly at merger or collapses within a few hundred ms. Within the narrow supramassive mass range, the lifetime decreases sharply. Computing $\tau_R$ as a function of $\tilde{m}_R$, however, is a difficult problem beyond the scope of this paper. Although we expect $\tilde{m}_R$ to be the most important parameter, note that $\tau_R$ might also depend on other factors such as mass ratio, spins, and equation of state.

In table 2, we provide the value $\tilde{m}_0$ corresponding to the median value of the total initial baryonic mass inferred for the narrow prior. Applying equation (4) to $M_B^R$ and $M_B^0$, we obtain

$$ \tilde{m}_R = \tilde{m}_0 - (\tilde{m}_0 + k_m) \frac{M_B^E}{M_B^0}. \qquad (5) $$

In [116], the ejecta mass was estimated from the observed kilonova to around $0.078 M_\odot$. If we assume more conservatively that $M_B^E < 0.1 M_\odot$, we find that $0 < \tilde{m}_0 - \tilde{m}_R < 0.3$ for the equations of state considered here (see table 2). The above ejecta mass limit should be regarded as a fiducial value since modeling the ejecta mass, either by numerical simulations or inferred from the kilonova, is an ongoing field of research.

For the narrow prior, we find that the MS1_PP, MS1B_PP, and MPA1 equation of state models lead to a stable remnant and are therefore incompatible with the black hole with disk short gamma-ray burst scenario (technically, the MPA1 also allows marginally supramassive remnants, which can for all practical purposes be considered stable). The H4, SLY, KDE0V, KDE0V1, SKOP, SLY2, BHF_BBB2, and HQC18 equation of state models result in a short-lived hypermassive remnant. For APR4_EPP, RS, SLY9, SLY230A, SKI2, SKI4, SK255, and SKMP, the remnant is in the upper supramassive or lower hypermassive mass range, and we cannot predict if the lifetime is longer or shorter than the observed short gamma-ray burst delay. For SKI3, SKI5, SKI6, and SK272 the remnant mass is in the supramassive range with $\tilde{m}_R < 0.9$, and likely long-lived.

We now turn to discuss whether the merger can lead to prompt black hole formation for a given equation of state model, based on gravitational-wave data alone. Prompt collapse occurs above some critical mass, which depends on the equation of state. Known examples of prompt collapse in numerical binary neutron star simulations generally involve hypermassive systems. If a parameter estimation run for a given equation of state model yields posterior support only in the supramassive and stable mass range, we assume that no prompt collapse can occur. For hypermassive systems, we use existing thresholds [18] from numerical simulations, given
in terms of total gravitational mass of the binary. We consider a systematic error of 0.1 $M_\odot$ due to the granularity of the simulated systems. The available thresholds are for equal-mass, non-spinning systems. To account for a possible impact of mass ratio and spin, we further increase the error by 0.05 $M_\odot$. This is a conservative estimate assuming that the main influence is by the change in dynamically ejected mass, which has been computed for many models in [124]. Finally, we express the thresholds for equal mass systems in terms of total baryonic mass and use it as an approximate threshold on $M_B^{tot}$ for the generic case.

For the narrow prior, we can rule out prompt collapse for the APR4_EPP, H4, MPA1, MS1_PP, MS1B_PP, SKI3, SKI4, SKI5, SKI6, and SKI272 equation of state models. The SLY equation of state model is compatible both with prompt and delayed collapse (due to the systematic error of the threshold mass from numerical relativity). For the remaining equations of state, no data on prompt collapse threshold is available.

With improved theoretical modeling, the wealth of electromagnetic observations may yield reliable constraints on the lifetime of the remnant. If prompt collapse could be ruled out, the values of $M_B^0$ obtained from gravitational-wave data could be taken as a lower limit for the prompt collapse threshold. Considering a potential impact of the mass ratio on the actual threshold, one has to take into account that the limits on $M_B^0$ are the result of a marginalization over mass ratio.

### 4.3. Constraining the maximum neutron star mass

Within the black hole with disk scenario of short gamma-ray bursts, the remnant cannot be stable. This was already exploited in [4] to derive limits on the maximum mass of non-rotating neutron stars. The delay of 1.7 s [4] between merger and short gamma-ray burst would also constrain the lifetime of the remnant. The relatively short lifetime indicates that the remnant either cannot be supported by uniform rotation alone (hypermassive remnant), or that it is at least close to the critical mass (how close exactly is an important open question).

In the following, we estimate the maximum mass of non-rotating neutron stars under the assumption that a black hole was produced in the binary neutron star merger before the observed short gamma-ray burst. For this, we first rewrite equation (4) as

$$ \hat{M}_B^{\text{TOV}} = M_B^0 \frac{1 - M_B^E / M_B^0}{1 + \tilde{m}_R / k_m} \le \frac{M_B^0}{1 + \tilde{m}_R / k_m} \quad (6) $$

From table 2, we find that $M_B^0$ varies only slightly with equation of state. Therefore, we take the bound $M_B^0 < 3.05 M_\odot$ for the equations of state considered here as a generic bound. In addition, we assume $k_m < 7$, as fulfilled by all equations of state considered here. The remaining open question is below which value $\tilde{m}_R$ the remnant lifetime is above the short gamma-ray burst delay. We already know that the lifetime is infinite unless $\tilde{m}_R > 0$. Future studies might provide tighter bounds, which will tighten the upper bound obtained for $\hat{M}_B^{\text{TOV}}$. For example, if we demand $\tilde{m}_R > 0.4$, then $\hat{M}_B^{\text{TOV}} < 2.9 M_\odot$.

Another potential method to limit the maximum mass of non-rotating neutron stars is to limit the angular momentum that needs to be removed before collapse can occur, and which becomes zero for $\tilde{m}_R = 1$. In [27], it was argued that the remnant mass must be close to the hypermassive range since the large spindown energy that would otherwise have been deposited in the environment would be incompatible with the observed counterparts (however, there are observational gaps [5], e.g. in x-rays, which could lead to energy output being missed). This assumption corresponds to $\tilde{m}_R \ge 1$. From equation (6), we would obtain a corresponding limit $\hat{M}_B^{\text{TOV}} < 2.67 M_\odot$ (for the narrow prior). If we further assume that the maximum
baryonic mass is 15%–23% larger than the maximum gravitational mass, as fulfilled for all equation of state models considered here, we obtain a limit $\hat{M}_{\text{G}}^{\text{TOV}} \le 2.32 M_{\odot}$. Based on additional assumptions, excluding prompt collapse and using an approximate expression for the prompt collapse threshold ignoring systematic errors, [27] provides a limit on the maximum neutron star mass which is 6% tighter. Reliable prompt collapse thresholds for a representative selection of equation of state models might hence be useful in this respect.

In [29], a limit $\hat{M}_{\text{G}}^{\text{TOV}} \le 2.28 \pm 0.23 M_{\odot}$ was provided, again based on the assumption of having a hypermassive remnant. The differences to our analysis are mainly the use of gravitational masses instead of conserved baryonic masses, and a more constraining assumption regarding the range of ratios between maximum gravitational masses for uniformly rotating and non-rotating neutron stars. Such effects are included in the systematic error given in [29] for their maximum mass limit, which is therefore compatible with our result.

A similar calculation can be found in [54], yielding a limit $\hat{M}_{\text{G}}^{\text{TOV}} \le 2.15-2.25 M_{\odot}$. It is based again on the detected limits for the total initial gravitational mass, and further assumes a total energy loss of $(0.15 \pm 0.03) M_{\odot}$, as well as a fixed difference $0.4 M_{\odot}$ between maximum gravitational masses of uniformly rotating and non-rotating neutron stars. These assumptions might explain the tighter limit compared to our result (for example, we find $\hat{M}_{\text{G}}^{\text{ROT}} - \hat{M}_{\text{G}}^{\text{TOV}} = 0.34 M_{\odot}$ for the BHF_BBB2 equation of state). A recent analysis [125], which is not relying on assumptions about $\tilde{m}$, indicates that the remnant at the time of collapse might not rotate rapidly, and provides a limit $\hat{M}_{\text{G}}^{\text{TOV}} \le 2.3 M_{\odot}$.

By using the detected gravitational-wave signal alone, without assumptions on the remnant fate, [40] inferred a limit $\hat{M}_{\text{G}}^{\text{TOV}} \le 2.46$, compatible with our results. This is based on modeling the equation of state using Gaussian processes, trained on a set of seven equations of state and respecting causality by construction. The Gaussian process model is further informed by the gravitational-wave data and then used to extrapolate the equation of state from densities relevant for the inspiral signal up to those required to determine the maximum mass.

### 4.4. Impact of prior and waveform model

In order to assess the waveform systematics, we repeat the computation of the remnant properties for a representative equation of state model using all waveform approximants. The results are shown in table 3. The differences for $M_{\text{B}}$ and $\tilde{m}_0$ are small compared to the statistical uncertainties, and negligible with regard to our conclusions on remnant classification, prompt collapse, and the maximum neutron star mass. This is to be expected, since those results depend directly only on the initial neutron star masses, but not on the tidal deformabilities, which are most sensitive to the waveform model.

Table 4 shows the results obtained using the broad prior restricted to the binary neutron star case. The restriction is necessary as the heavier object exceeds the maximum mass allowed for a neutron star for parts of the broad prior. Comparing to the narrow prior in table 2, we find that the total initial baryonic mass distributions are broadened towards higher masses. This is a consequence of the wider posterior distribution for the mass ratio; for fixed chirp mass, the total baryonic mass for unequal mass binaries is higher than for equal masses.

The broader distribution leads to differences regarding the fate of the remnant. While we ruled out prompt collapse for the H4 and APR4_EPP equations of state with the narrow prior, the respective total mass posteriors for the broad prior also have support in the mass range where prompt collapse is allowed according to the thresholds from numerical relativity. This
<table><caption>**Table 3.** Implications derived assuming the MPA1 equation of state model and narrow prior, using different waveform approximants. The given quantities are explained in table 2.</caption><thead><tr><th>Waveform</th><th>M<sub>B</sub> (M<sub>☉</sub>)</th><th>E<sub>c</sub>/c<sup>2</sup><br>(10<sup>15</sup> g cm<sup>-3</sup>)</th><th>m&#773;<sub>0</sub></th><th>&Delta;<sub>0.1</sub></th><th>I<sub>max</sub><br>(M&#773;<sup>3</sup><sub>☉</sub>)</th><th>F<sub>max</sub><br>(kHz)</th></tr></thead><tbody><tr><td>TaylorF2(7.5 PN)</td><td>3.004<sup>+0.016</sup><sub>-0.009</sub></td><td>0.74</td><td>-0.016<sup>+0.029</sup><sub>-0.016</sub></td><td>0.18</td><td>114</td><td>1.34</td></tr><tr><td>TaylorF2(7 PN)</td><td>3.004<sup>+0.016</sup><sub>-0.009</sub></td><td>0.74</td><td>-0.017<sup>+0.029</sup><sub>-0.015</sub></td><td>0.18</td><td>114</td><td>1.34</td></tr><tr><td>TaylorF2(6 PN)</td><td>3.004<sup>+0.016</sup><sub>-0.009</sub></td><td>0.74</td><td>-0.016<sup>+0.028</sup><sub>-0.016</sub></td><td>0.18</td><td>114</td><td>1.34</td></tr><tr><td>EOB+NRT</td><td>3.004<sup>+0.016</sup><sub>-0.008</sub></td><td>0.74</td><td>-0.015<sup>+0.027</sup><sub>-0.015</sub></td><td>0.18</td><td>114</td><td>1.34</td></tr><tr><td>Phenom+NRT</td><td>3.005<sup>+0.016</sup><sub>-0.009</sub></td><td>0.74</td><td>-0.015<sup>+0.029</sup><sub>-0.015</sub></td><td>0.18</td><td>114</td><td>1.34</td></tr></tbody></table>

result might still change with more accurate numerical-relativity results for the threshold, given its large systematic error.

The KDE0V, SKI3, SKI4, SKI5, SKI6, and SK272 equations of state, which lead to supramassive remnants assuming the narrow prior, might also lead to hypermassive remnants when assuming the broad prior. The MS1_PP, MS1B_PP equations of state still lead to stable remnants, while the MPA1 equation of state model can also result in a long-lived supramassive remnant for the broad prior. The bound on the maximum Tolman–Oppenheimer–Volkoff mass when assuming the black hole with disk short gamma-ray burst model becomes less strict as well. Using only $\tilde{m}_R \ge 0$, we obtain $\hat{M}_B^{\text{TOV}} < 3.4M_☉$. Assuming $\tilde{m}_R \ge 1$, we would obtain $\hat{M}_B^{\text{TOV}} < 2.98M_☉$.

Another consequence of the larger mass ratios is that the mass posterior of the heavier neutron star extends to larger values. Therefore, higher central energy densities are reached. Depending on the equation of state, the inspiral gravitational-wave signal might be sensitive to the equation of state at energy densities up to $2 \times 10^{15} \text{ g cm}^{-3}$ (compare table 4).

### 4.5. Remnant rotation rate

In case the remnant is long-lived, the maximum rotation rate and moment of inertia are of interest for searches of long-lived gravitational-wave signals from the remnant (see [20, 22]), and might also be relevant for modeling the magnetar short gamma-ray burst scenario. In the following, we estimate upper limits for rotation rate and moment of inertia.

The total angular momentum at merger is not known, although numerical simulations indicate that it typically exceeds the maximum amount that could be realized in a uniformly rotating system of the same mass (see [21]). However, the disk can account for a significant fraction of the total angular momentum. Most of the disk angular momentum is transported outwards during accretion (together with a fraction of the mass). The angular momentum of the final remnant could depend on the details of the accretion process and might be less than the critical one (the opposite assumption is made in [21]), and it can decrease further on longer timescales, e.g. by magnetic spindown. The rotation rate of models near the maximum mass can increase with decreasing angular momentum, because the compactness also increases. The maximum rotation rate can in principle exceed the rotation rate at the mass-shedding limit. Here, we use the rotation rate at the mass-shedding limit as an estimate for the maximum rotation rate of the remnant.

From the mass posteriors for each equation of state model, we compute the posteriors for the rotation rate and moment of inertia of uniformly rotating neutron stars with baryonic mass equal to the total baryonic mass and maximum angular momentum. Tables 2 and 4 show the
<table><caption>**Table 4.** Like table 2, but showing results for the broad prior, again with the TaylorF2 (7.5 PN) waveform model.</caption><thead><tr><th>Equation of state</th><th>M<sup>0</sup><sub>B</sub>(M<sub>☉</sub>)</th><th>E<sub>c</sub>/c<sup>2</sup><br>(10<sup>15</sup> g cm<sup>-3</sup>)</th><th>m&#773;<sub>0</sub></th><th>&Delta;<sub>0.1</sub></th><th>Prompt<br>collapse</th><th>I<sub>max</sub><br>(M<sup>3</sup><sub>☉</sub>)</th><th>F<sub>max</sub><br>(kHz)</th></tr></thead><tbody><tr><td>APR4_EPP</td><td>3.084<sup>+0.316</sup><sub>-0.070</sub></td><td>1.48</td><td>1.019<sup>+0.683</sup><sub>-0.150</sub></td><td>0.22</td><td>Possibly</td><td>84</td><td>1.89</td></tr><tr><td>MS1_PP</td><td>2.960<sup>+0.121</sup><sub>-0.026</sub></td><td>0.53</td><td>-0.574<sup>+0.202</sup><sub>-0.044</sub></td><td>0.17</td><td>No</td><td>158</td><td>1.03</td></tr><tr><td>MS1B_PP</td><td>2.973<sup>+0.144</sup><sub>-0.030</sub></td><td>0.55</td><td>-0.542<sup>+0.231</sup><sub>-0.048</sub></td><td>0.16</td><td>No</td><td>159</td><td>1.05</td></tr><tr><td>KDE0V</td><td>3.058<sup>+0.219</sup><sub>-0.053</sub></td><td>1.80</td><td>1.978<sup>+0.582</sup><sub>-0.139</sub></td><td>0.27</td><td>—</td><td>66</td><td>1.86</td></tr><tr><td>SKI6</td><td>3.028<sup>+0.225</sup><sub>-0.041</sub></td><td>1.08</td><td>0.922<sup>+0.489</sup><sub>-0.087</sub></td><td>0.22</td><td>—</td><td>97</td><td>1.66</td></tr><tr><td>H4</td><td>3.007<sup>+0.178</sup><sub>-0.035</sub></td><td>0.94</td><td>1.450<sup>+0.388</sup><sub>-0.077</sub></td><td>0.22</td><td>Possibly</td><td>104</td><td>1.36</td></tr><tr><td>MPA1</td><td>3.041<sup>+0.255</sup><sub>-0.043</sub></td><td>0.90</td><td>0.050<sup>+0.456</sup><sub>-0.077</sub></td><td>0.18</td><td>No</td><td>121</td><td>1.41</td></tr><tr><td>SK255</td><td>3.003<sup>+0.203</sup><sub>-0.041</sub></td><td>1.10</td><td>1.207<sup>+0.497</sup><sub>-0.102</sub></td><td>0.24</td><td>—</td><td>91</td><td>1.61</td></tr><tr><td>SKI2</td><td>3.001<sup>+0.176</sup><sub>-0.038</sub></td><td>0.97</td><td>1.100<sup>+0.408</sup><sub>-0.089</sub></td><td>0.23</td><td>—</td><td>100</td><td>1.53</td></tr><tr><td>SKI4</td><td>3.028<sup>+0.228</sup><sub>-0.037</sub></td><td>1.11</td><td>0.978<sup>+0.494</sup><sub>-0.081</sub></td><td>0.22</td><td>—</td><td>95</td><td>1.67</td></tr><tr><td>SKOP</td><td>3.019<sup>+0.190</sup><sub>-0.032</sub></td><td>1.45</td><td>1.874<sup>+0.500</sup><sub>-0.085</sub></td><td>0.26</td><td>—</td><td>73</td><td>1.70</td></tr><tr><td>SLY2</td><td>3.043<sup>+0.239</sup><sub>-0.045</sub></td><td>1.48</td><td>1.494<sup>+0.585</sup><sub>-0.110</sub></td><td>0.25</td><td>—</td><td>77</td><td>1.79</td></tr><tr><td>SLY9</td><td>3.018<sup>+0.234</sup><sub>-0.037</sub></td><td>1.15</td><td>1.068<sup>+0.537</sup><sub>-0.083</sub></td><td>0.23</td><td>—</td><td>91</td><td>1.68</td></tr><tr><td>BHF_BBB2</td><td>3.077<sup>+0.201</sup><sub>-0.063</sub></td><td>2.02</td><td>2.181<sup>+0.545</sup><sub>-0.168</sub></td><td>0.27</td><td>—</td><td>63</td><td>1.87</td></tr><tr><td>HQC18</td><td>3.064<sup>+0.253</sup><sub>-0.057</sub></td><td>1.42</td><td>1.418<sup>+0.561</sup><sub>-0.129</sub></td><td>0.22</td><td>—</td><td>81</td><td>1.72</td></tr><tr><td>KDE0V1</td><td>3.049<sup>+0.216</sup><sub>-0.051</sub></td><td>1.73</td><td>1.953<sup>+0.580</sup><sub>-0.135</sub></td><td>0.27</td><td>—</td><td>68</td><td>1.80</td></tr><tr><td>RS</td><td>3.016<sup>+0.216</sup><sub>-0.042</sub></td><td>1.14</td><td>1.259<sup>+0.507</sup><sub>-0.099</sub></td><td>0.24</td><td>—</td><td>92</td><td>1.61</td></tr><tr><td>SK272</td><td>3.003<sup>+0.209</sup><sub>-0.044</sub></td><td>1.00</td><td>0.868<sup>+0.483</sup><sub>-0.102</sub></td><td>0.23</td><td>—</td><td>100</td><td>1.61</td></tr><tr><td>SKI3</td><td>3.000<sup>+0.190</sup><sub>-0.039</sub></td><td>0.91</td><td>0.818<sup>+0.421</sup><sub>-0.087</sub></td><td>0.22</td><td>—</td><td>107</td><td>1.53</td></tr><tr><td>SKI5</td><td>2.975<sup>+0.116</sup><sub>-0.024</sub></td><td>0.78</td><td>0.816<sup>+0.262</sup><sub>-0.056</sub></td><td>0.23</td><td>—</td><td>109</td><td>1.47</td></tr><tr><td>SKMP</td><td>3.019<sup>+0.216</sup><sub>-0.035</sub></td><td>1.19</td><td>1.254<sup>+0.504</sup><sub>-0.083</sub></td><td>0.23</td><td>—</td><td>88</td><td>1.68</td></tr><tr><td>SLY</td><td>3.042<sup>+0.238</sup><sub>-0.044</sub></td><td>1.50</td><td>1.517<sup>+0.590</sup><sub>-0.109</sub></td><td>0.25</td><td>Possibly</td><td>76</td><td>1.82</td></tr><tr><td>SLY230A</td><td>3.044<sup>+0.251</sup><sub>-0.044</sub></td><td>1.36</td><td>1.261<sup>+0.581</sup><sub>-0.101</sub></td><td>0.23</td><td>—</td><td>83</td><td>1.76</td></tr></tbody></table>

numerical results for the narrow and broad priors. For both, we find a maximum rotation rate below 1.9kHz for all equations of state. We also provide values for remnants in the hypermagnetic mass range. For those cases, we limit the baryonic mass above to the maximum allowed for uniform rotation. The numbers we obtain can be regarded as estimates for the hypothetical case that the excess mass was somehow lost during merger, or remains in a disk.

## 5. Conclusions

In this work, we employ Bayesian model selection as a tool to discriminate between different theoretical models for the neutron star equation of state using only the gravitational-wave data. This allows us to make quantitative statements about the relative likelihood of a given equation of state model. This work complements the analysis of tidal effects of GW170817 in [6, 38] and the analysis of the neutron star radii and equation of state in [39]. As in previous studies [6, 38], we assume a priori that the source of GW170817 resides in NGC 4993 and adopt fixed sky location consistent with electromagnetic observations. The distance to GW170817 is determined from the gravitational-wave data alone, as in [6]. We do not allow
tidal parameters to vary independently and instead enforce that each neutron star obeys a given equation of state model.

The Bayesian evidence is computed using two independent methods that show good agreement and demonstrate the robustness of our results. Under minimal assumptions on the astrophysical priors for the component masses and spins, we find that the models predicting the largest fiducial radii and tidal deformabilities are *very strongly disfavored* (MS1_PP) or at least *strongly disfavored* (MS1B_PP) with respect to the equations of state that predict a small fiducial radius and tidal deformability, such as KDE0V. In addition, we find that for both priors considered in this work the gravitational-wave data does not favor and nor can it rule out any particular equation of state against the low-mass binary black hole hypothesis.

In order to gauge waveform systematics, we compare the Bayes factors using different waveform approximants. The results show good agreement between the different waveform models, providing confidence that systematic uncertainties are small compared to statistical uncertainties.

In a broader context, we also address the question as to whether we can exclude the possibility of GW170817 being a binary black hole by computing Bayes factors against the binary neutron star and neutron star-black hole hypotheses. From the gravitational-wave data alone, no single hypothesis is strongly preferred and we are unable to exclude the binary black hole hypothesis. We find that the Bayes factors computed using different waveform approximants are in good agreement, supporting the notion that waveform systematics are small compared to statistical errors.

For some equation of state models, we can predict the outcome of the merger; for two models the remnant has to be a stable neutron star, for some it will be a long-lived supramassive neutron star, and for some a short-lived hypermassive neutron star or black hole. We can further rule out prompt collapse for a number of models. Considering the Bayes factor for each model, we find that only indefinitely stable remnants are disfavored. Based on gravitational-wave data alone, all scenarios from very long-lived remnants to prompt black hole formation at merger remain possible. Further constraints on the outcome will likely have to rely on the electromagnetic counterparts and modeling of the kilonova and the short gamma-ray burst. We find that making the assumption of a short-lived remnant—as required by some short gamma-ray burst models—would allow us to put further constraints on the equation of state model. In particular, we improve the upper limit on the maximum mass of non-rotating neutron stars, compared to an earlier study [4] (which is based on an analysis that did not assume that both neutron stars obey the same equation of state). We also find that three equation of state models are ruled out if a black hole was formed before the short gamma-ray burst, including the two disfavored by gravitational-wave data alone.

## Acknowledgments

The authors gratefully acknowledge the support of the United States National Science Foundation (NSF) for the construction and operation of the LIGO Laboratory and Advanced LIGO as well as the Science and Technology Facilities Council (STFC) of the United Kingdom, the Max-Planck-Society (MPS), and the State of Niedersachsen/Germany for support of the construction of Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. The authors gratefully acknowledge the Italian Istituto Nazionale di Fisica Nucleare (INFN), the French Centre National de la Recherche Scientifique (CNRS) and the Foundation for Fundamental Research on Matter supported by the Netherlands Organisation for Scientific
Research, for the construction and operation of the Virgo detector and the creation and sup-
port of the EGO consortium. The authors also gratefully acknowledge research support from
these agencies as well as by the Council of Scientific and Industrial Research of India, the
Department of Science and Technology, India, the Science & Engineering Research Board
(SERB), India, the Ministry of Human Resource Development, India, the Spanish Agencia
Estatal de Investigación, the Vicepresidència i Conselleria d'Innovació, Recerca i Turisme
and the Conselleria d'Educació i Universitat del Govern de les Illes Balears, the Conselleria
d'Educació, Investigació, Cultura i Esport de la Generalitat Valenciana, the National Science
Centre of Poland, the Swiss National Science Foundation (SNSF), the Russian Foundation for
Basic Research, the Russian Science Foundation, the European Commission, the European
Regional Development Funds (ERDF), the Royal Society, the Scottish Funding Council, the
Scottish Universities Physics Alliance, the Hungarian Scientific Research Fund (OTKA),
the Lyon Institute of Origins (LIO), the Paris Île-de-France Region, the National Research,
Development and Innovation Office Hungary (NKFIH), the National Research Foundation
of Korea, Industry Canada and the Province of Ontario through the Ministry of Economic
Development and Innovation, the Natural Science and Engineering Research Council Canada,
the Canadian Institute for Advanced Research, the Brazilian Ministry of Science, Technology,
Innovations, and Communications, the International Center for Theoretical Physics South
American Institute for Fundamental Research (ICTP-SAIFR), the Research Grants Council
of Hong Kong, the National Natural Science Foundation of China (NSFC), the Leverhulme
Trust, the Research Corporation, the Ministry of Science and Technology (MOST), Taiwan
and the Kavli Foundation. The authors gratefully acknowledge the support of the NSF, STFC,
MPS, INFN, CNRS and the State of Niedersachsen/Germany for provision of computational
resources.

**Appendix A. Equations of state**

In the following, we provide details for the equations of state employed in this work, all of
which represent estimates from nuclear physics. They employ different estimates for the forces
between constituents, for example from nucleon-nucleon scattering experiments. Those forces
are often parametrized (e.g. Skyrme forces) before use in equation of state computations. Other
important observable parameters to calibrate equation of state models are the nuclear symme-
try energy and incompressibility. Further, the many-particle quantum field theory problem is
approximately solved using different schemes, such as liquid-drop model, relativistic/non-
relativistic mean field (MF) theory, Hartree-Fock (HF) method, Brueckner-Bethe-Goldstone
(BG) methods, or Wigner-Seitz (WS) cells. Some methods apply the simplification of a single
microstate, while others consider nuclear statistic equilibrium. One important aspect is which
types of constituents are considered, for example if hyperons are included or not. Most equa-
tion of state models considered describe hadronic matter, except the HQC18 equation of state
model, which assumes a hadronic crust around a quark matter core [9].

For details on the various models, we refer to the references given in table A1, which
list all our equation of state models, the respective sources for the numeric tables, and the
original publications. Further, we provide the maximum gravitational and baryonic masses
of non-rotating neutron stars, and the maximum baryonic mass of uniformly rotating neu-
tron stars. The latter were obtained using the Rapidly Rotating Neutron Star (RNS) code
described in [126].

Most of our equations of state were obtained from three catalogs. The CompOSE project
[127] provides nuclear physics equation of state data and a supporting software framework.
<table><caption>**Table A1.** Equations of state used in this work. We provide the original publications introducing the equations of state and the source of the numerical tables or piecewise polytropic fit parameters. For each we computed the maximum gravitational mass $\dot{M}_{G}^{TOV}$ of non-rotating neutron stars and the corresponding baryonic mass $\dot{M}_{B}^{TOV}$, as well as the maximum baryonic mass $\dot{M}_{B}^{ROT}$ of uniformly rotating neutron star. The WFF1 equation of state model becomes acausal below the central density of the maximum mass model, and the reported values (in brackets) are the maxima of neutron star models respecting causality.</caption><thead><tr><th>EoS</th><th>References</th><th>Table</th><th>Methods</th><th>dotM<sub>G</sub><sup>TOV</sup><br>(M<sub>☉</sub>)</th><th>dotM<sub>B</sub><sup>TOV</sup><br>(M<sub>☉</sub>)</th><th>dotM<sub>B</sub><sup>ROT</sup><br>(M<sub>☉</sub>)</th></tr></thead><tbody><tr><td>BHF_BBB2</td><td>[131]</td><td>[127]</td><td>Brueckner-BG</td><td>1.92</td><td>2.27</td><td>2.64</td></tr><tr><td>WFF1</td><td>[98]</td><td>[128]</td><td>Variational</td><td>(1.92)</td><td>(2.31)</td><td>(2.84)</td></tr><tr><td>KDE0V</td><td>[132–134]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>1.96</td><td>2.31</td><td>2.69</td></tr><tr><td>KDE0V1</td><td>[132–134]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>1.97</td><td>2.32</td><td>2.69</td></tr><tr><td>SKOP</td><td>[132, 133, 135]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>1.97</td><td>2.30</td><td>2.69</td></tr><tr><td>H4</td><td>[136]</td><td>[128]</td><td>Relativistic MF</td><td>2.03</td><td>2.34</td><td>2.80</td></tr><tr><td>HQC18</td><td>[9]</td><td></td><td>Hadron-Quark</td><td>2.05</td><td>2.43</td><td>2.88</td></tr><tr><td>SLY</td><td>[137]</td><td>[128]</td><td>Potential-method</td><td>2.05</td><td>2.43</td><td>2.83</td></tr><tr><td>SLY2</td><td>[132, 133]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.05</td><td>2.43</td><td>2.84</td></tr><tr><td>SLY230A</td><td>[132, 133, 138]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.10</td><td>2.50</td><td>2.93</td></tr><tr><td>SKMP</td><td>[132, 133, 139]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.11</td><td>2.48</td><td>2.91</td></tr><tr><td>RS</td><td>[132, 133, 140]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.12</td><td>2.48</td><td>2.91</td></tr><tr><td>SK255</td><td>[132, 133, 141]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.14</td><td>2.51</td><td>2.92</td></tr><tr><td>SLY9</td><td>[132, 133]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.16</td><td>2.55</td><td>2.99</td></tr><tr><td>APR4_EPP</td><td>[129]</td><td>[130]</td><td>Variational</td><td>2.16</td><td>2.61</td><td>3.07</td></tr><tr><td>SKI2</td><td>[132, 133, 142]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.16</td><td>2.53</td><td>2.96</td></tr><tr><td>SKI4</td><td>[132, 133, 142]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.17</td><td>2.58</td><td>3.04</td></tr><tr><td>SKI6</td><td>[132, 133, 143]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.19</td><td>2.60</td><td>3.06</td></tr><tr><td>SK272</td><td>[132, 133, 141]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.23</td><td>2.63</td><td>3.06</td></tr><tr><td>SKI3</td><td>[132, 133, 142]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.24</td><td>2.63</td><td>3.08</td></tr><tr><td>SKI5</td><td>[132, 133, 142]</td><td>[127]</td><td>Skyrme,WS,HF</td><td>2.24</td><td>2.62</td><td>3.06</td></tr><tr><td>MPA1</td><td>[144]</td><td>[128]</td><td>Relativistic HF</td><td>2.46</td><td>3.01</td><td>3.57</td></tr><tr><td>MS1B_PP</td><td>[145]</td><td>[47]</td><td>Relativistic MF</td><td>2.75</td><td>3.31</td><td>3.93</td></tr><tr><td>MS1_PP</td><td>[145]</td><td>[47]</td><td>Relativistic MF</td><td>2.75</td><td>3.30</td><td>3.90</td></tr></tbody></table>

We used this public code (selecting the first-order interpolation option and a resolution of 500 points) to create tables. Another collection of tables is provided by [128]. To obtain the final, continuous equations of state, we employ two different interpolation schemes. The results in section 4 and table A1 are computed using monotonic spline interpolation (PCHIP), while the equation of state related code used during parameter estimation employs standard cubic spline interpolation instead. For all models, we compared sequences of non-rotating neutron stars between the two. Since the low resolution and gaps present in the MS1 and MS1B tables cause significant ambiguities in the tidal deformability and radius, we use well defined piecewise polytropic approximations instead, taken from a third catalog [47] providing analytic fits, and denoted here by MS1_PP and MS1B_PP, respectively. When comparing with [6], note that their WFF1, H4, MPA1, SLY and APR4 models also refer to the piecewise polytropic variants from [47], in contrast to this work.

We applied the following adjustments: in general, we removed isolated data points from the original tables that violate thermodynamic constraints for zero temperature equations of state. For the H4 equation of state model from [128], the low-density part (<10¹³ g cm⁻³) was
<table><caption>**Table B1.** Bayes factors computed using LALINFERENCE's nested sampling algorithm for all of the equations of state discussed in section 2.4 using the narrow prior discussed in section 2.2. The Bayes factors are computed for a number of waveform models, as discussed in section 2.3. Here the Bayes factors are taken to be the ratio between the evidence for a given equation of state and the corresponding evidence for the hypothesis that both bodies are low-mass black holes. The errors correspond to statistical uncertainty in the calculation of the evidence from the nested sampling algorithm in LALINFERENCE, where $\sigma_{BF} = \sigma_{Z_{EOS}} / \bar{Z}_{BBH}$.</caption><thead><tr><th rowspan="2">Equation of state</th><th colspan="5">LALINFERENCE nest results</th></tr><tr><th>TaylorF2 (7.5pN)</th><th>Phenom+NRT</th><th>EOB+NRT</th><th>TaylorF2 (6pN)</th><th>TaylorF2 (7pN)</th></tr></thead><tbody><tr><td>APR4_EPP</td><td>(1.92 ± 0.07) × 10<sup>1</sup></td><td>(4.39 ± 0.23) × 10<sup>1</sup></td><td>(2.02 ± 0.08) × 10<sup>1</sup></td><td>(1.76 ± 0.06) × 10<sup>1</sup></td><td>(1.78 ± 0.10) × 10<sup>1</sup></td></tr><tr><td>BBH</td><td>(1.00 ± 0.05)</td><td>(1.00 ± 0.04)</td><td>(1.00 ± 0.05)</td><td>(1.00 ± 0.05)</td><td>(1.00 ± 0.05)</td></tr><tr><td>BHF_BBB2</td><td>(1.60 ± 0.06) × 10<sup>1</sup></td><td>(4.79 ± 0.19) × 10<sup>1</sup></td><td>(2.43 ± 0.07) × 10<sup>1</sup></td><td>(1.70 ± 0.07) × 10<sup>1</sup></td><td>(1.71 ± 0.09) × 10<sup>1</sup></td></tr><tr><td>H4</td><td>(1.73 ± 0.08)</td><td>(1.05 ± 0.06)</td><td>(3.80 ± 0.16) × 10<sup>-1</sup></td><td>(1.73 ± 0.08)</td><td>(1.76 ± 0.09)</td></tr><tr><td>HQC18</td><td>(1.94 ± 0.07) × 10<sup>1</sup></td><td>(3.84 ± 0.20) × 10<sup>1</sup></td><td>(1.95 ± 0.12) × 10<sup>1</sup></td><td>(1.73 ± 0.07) × 10<sup>1</sup></td><td>(1.81 ± 0.09) × 10<sup>1</sup></td></tr><tr><td>KDE0V</td><td>(1.96 ± 0.11) × 10<sup>1</sup></td><td>(3.95 ± 0.17) × 10<sup>1</sup></td><td>(1.90 ± 0.08) × 10<sup>1</sup></td><td>(1.80 ± 0.10) × 10<sup>1</sup></td><td>(1.63 ± 0.06) × 10<sup>1</sup></td></tr><tr><td>KDE0V1</td><td>(1.78 ± 0.09) × 10<sup>1</sup></td><td>(3.28 ± 0.08) × 10<sup>1</sup></td><td>(1.54 ± 0.07) × 10<sup>1</sup></td><td>(1.86 ± 0.10) × 10<sup>1</sup></td><td>(1.53 ± 0.08) × 10<sup>1</sup></td></tr><tr><td>MPA1</td><td>(5.56 ± 0.23)</td><td>(1.35 ± 0.05) × 10<sup>1</sup></td><td>(7.32 ± 0.33)</td><td>(4.66 ± 0.20)</td><td>(4.36 ± 0.19)</td></tr><tr><td>MS1B_PP</td><td>(3.48 ± 0.13) × 10<sup>-1</sup></td><td>(1.61 ± 0.07) × 10<sup>-1</sup></td><td>(7.03 ± 0.37) × 10<sup>-2</sup></td><td>(1.86 ± 0.07) × 10<sup>-1</sup></td><td>(1.92 ± 0.09) × 10<sup>-1</sup></td></tr><tr><td>MS1_PP</td><td>(8.96 ± 0.38) × 10<sup>-2</sup></td><td>(8.45 ± 0.46) × 10<sup>-2</sup></td><td>(3.36 ± 0.14) × 10<sup>-2</sup></td><td>(3.44 ± 0.16) × 10<sup>-2</sup></td><td>(3.44 ± 0.14) × 10<sup>-2</sup></td></tr><tr><td>RS</td><td>(3.66 ± 0.15)</td><td>(9.85 ± 0.46)</td><td>(3.73 ± 0.17)</td><td>(3.78 ± 0.15)</td><td>(4.05 ± 0.15)</td></tr><tr><td>SK255</td><td>(3.75 ± 0.19)</td><td>(9.56 ± 0.35)</td><td>(3.79 ± 0.17)</td><td>(3.77 ± 0.16)</td><td>(3.75 ± 0.15)</td></tr><tr><td>SK272</td><td>(3.22 ± 0.11)</td><td>(5.54 ± 0.22)</td><td>(2.74 ± 0.11)</td><td>(3.65 ± 0.17)</td><td>(3.55 ± 0.20)</td></tr><tr><td>SKI2</td><td>(2.07 ± 0.10)</td><td>(4.06 ± 0.18)</td><td>(1.63 ± 0.05)</td><td>(2.13 ± 0.10)</td><td>(2.69 ± 0.14)</td></tr><tr><td>SKI3</td><td>(2.03 ± 0.08)</td><td>(3.50 ± 0.13)</td><td>(1.46 ± 0.08)</td><td>(2.23 ± 0.07)</td><td>(2.54 ± 0.06)</td></tr><tr><td>SKI4</td><td>(7.35 ± 0.32)</td><td>(1.35 ± 0.06) × 10<sup>1</sup></td><td>(7.66 ± 0.22)</td><td>(4.99 ± 0.22)</td><td>(4.32 ± 0.19)</td></tr><tr><td>SKI5</td><td>(1.15 ± 0.05)</td><td>(3.40 ± 0.16) × 10<sup>-1</sup></td><td>(1.38 ± 0.05) × 10<sup>-1</sup></td><td>(6.39 ± 0.26) × 10<sup>-1</sup></td><td>(7.08 ± 0.29) × 10<sup>-1</sup></td></tr><tr><td>SKI6</td><td>(6.38 ± 0.28)</td><td>(1.42 ± 0.07) × 10<sup>1</sup></td><td>(6.67 ± 0.22)</td><td>(4.34 ± 0.20)</td><td>(4.24 ± 0.20)</td></tr><tr><td>SKMP</td><td>(6.80 ± 0.26)</td><td>(1.46 ± 0.07) × 10<sup>1</sup></td><td>(7.32 ± 0.37)</td><td>(4.43 ± 0.16)</td><td>(4.15 ± 0.17)</td></tr><tr><td>SKOP</td><td>(1.50 ± 0.07) × 10<sup>1</sup></td><td>(1.02 ± 0.06) × 10<sup>1</sup></td><td>(4.72 ± 0.14)</td><td>(9.46 ± 0.42)</td><td>(8.79 ± 0.39)</td></tr></tbody></table>

(Continued)
<table><caption>**Table B1. (Continued)**</caption><thead><tr><th>Equation of state</th><th>TaylorF2 (7.5pN)</th><th>Phenom+NRT</th><th>EOB+NRT</th><th>TaylorF2 (6pN)</th><th>TaylorF2 (7pN)</th></tr></thead><tbody><tr><td>SLY</td><td>(1.85 ± 0.09) × 10<sup>1</sup></td><td>(1.82 ± 0.09) × 10<sup>1</sup></td><td>(6.69 ± 0.29)</td><td>(1.55 ± 0.06) × 10<sup>1</sup></td><td>(1.39 ± 0.06) × 10<sup>1</sup></td></tr><tr><td>SLY2</td><td>(1.68 ± 0.05) × 10<sup>1</sup></td><td>(2.00 ± 0.07) × 10<sup>1</sup></td><td>(7.90 ± 0.31)</td><td>(1.53 ± 0.07) × 10<sup>1</sup></td><td>(1.40 ± 0.09) × 10<sup>1</sup></td></tr><tr><td>SLY230A</td><td>(1.80 ± 0.07) × 10<sup>1</sup></td><td>(1.47 ± 0.05) × 10<sup>1</sup></td><td>(6.26 ± 0.28)</td><td>(1.42 ± 0.06) × 10<sup>1</sup></td><td>(1.18 ± 0.04) × 10<sup>1</sup></td></tr><tr><td>SLY9</td><td>(8.93 ± 0.43)</td><td>(1.42 ± 0.06) × 10<sup>1</sup></td><td>(8.11 ± 0.28)</td><td>(5.46 ± 0.23)</td><td>(4.76 ± 0.20)</td></tr><tr><td>WFF1</td><td>(1.06 ± 0.05) × 10<sup>1</sup></td><td>(3.49 ± 0.16) × 10<sup>1</sup></td><td>(2.45 ± 0.11) × 10<sup>1</sup></td><td>(1.27 ± 0.05) × 10<sup>1</sup></td><td>(1.34 ± 0.06) × 10<sup>1</sup></td></tr></tbody></table>
<table><caption>**Table B2.** Similarly to table B1, we show the Bayes factors computed using LALINFERENCE's nested sampling algorithm for all of the equations of state discussed in section 2.4 using the broad prior discussed in section 2.2.</caption><thead><tr><th rowspan="2">Equation of state</th><th colspan="3">LALINFERENCE nest results</th></tr><tr><th>TaylorF2 (7.5pN)</th><th>Phenom+NRT</th><th>EOB+NRT</th></tr></thead><tbody><tr><td>APR4_EPP</td><td>(3.21 ± 0.13)</td><td>(2.56 ± 0.08)</td><td>(1.40 ± 0.06)</td></tr><tr><td>BBH</td><td>(1.00 ± 0.05)</td><td>(1.00 ± 0.07)</td><td>(1.00 ± 0.06)</td></tr><tr><td>BHF_BBB2</td><td>(3.42 ± 0.18)</td><td>(3.23 ± 0.17)</td><td>(1.86 ± 0.14)</td></tr><tr><td>H4</td><td>(1.88 ± 0.10) × 10<sup>-1</sup></td><td>(5.21 ± 0.42) × 10<sup>-2</sup></td><td>(2.86 ± 0.10) × 10<sup>-2</sup></td></tr><tr><td>HQC18</td><td>(2.97 ± 0.16)</td><td>(2.06 ± 0.10)</td><td>(1.27 ± 0.06)</td></tr><tr><td>KDE0V</td><td>(3.12 ± 0.13)</td><td>(1.98 ± 0.12)</td><td>(1.51 ± 0.10)</td></tr><tr><td>KDE0V1</td><td>(2.88 ± 0.13)</td><td>(1.71 ± 0.09)</td><td>(9.97 ± 0.41) × 10<sup>-1</sup></td></tr><tr><td>MPA1</td><td>(6.15 ± 0.34) × 10<sup>-1</sup></td><td>(8.30 ± 0.36) × 10<sup>-1</sup></td><td>(5.54 ± 0.41) × 10<sup>-1</sup></td></tr><tr><td>MS1B_PP</td><td>(3.65 ± 0.19) × 10<sup>-2</sup></td><td>(9.71 ± 0.38) × 10<sup>-3</sup></td><td>(6.94 ± 0.36) × 10<sup>-3</sup></td></tr><tr><td>MS1_PP</td><td>(9.68 ± 0.45) × 10<sup>-3</sup></td><td>(4.56 ± 0.23) × 10<sup>-3</sup></td><td>(2.80 ± 0.15) × 10<sup>-3</sup></td></tr><tr><td>RS</td><td>(3.89 ± 0.18) × 10<sup>-1</sup></td><td>(3.43 ± 0.20) × 10<sup>-1</sup></td><td>(2.51 ± 0.14) × 10<sup>-1</sup></td></tr><tr><td>SK255</td><td>(4.00 ± 0.27) × 10<sup>-1</sup></td><td>(3.57 ± 0.22) × 10<sup>-1</sup></td><td>(2.49 ± 0.19) × 10<sup>-1</sup></td></tr><tr><td>SK272</td><td>(3.44 ± 0.15) × 10<sup>-1</sup></td><td>(2.63 ± 0.22) × 10<sup>-1</sup></td><td>(1.99 ± 0.12) × 10<sup>-1</sup></td></tr><tr><td>SKI2</td><td>(2.05 ± 0.09) × 10<sup>-1</sup></td><td>(1.51 ± 0.13) × 10<sup>-1</sup></td><td>(8.66 ± 0.50) × 10<sup>-2</sup></td></tr><tr><td>SKI3</td><td>(2.09 ± 0.09) × 10<sup>-1</sup></td><td>(1.45 ± 0.12) × 10<sup>-1</sup></td><td>(8.26 ± 0.53) × 10<sup>-2</sup></td></tr><tr><td>SKI4</td><td>(8.19 ± 0.42) × 10<sup>-1</sup></td><td>(8.30 ± 0.57) × 10<sup>-1</sup></td><td>(6.17 ± 0.30) × 10<sup>-1</sup></td></tr><tr><td>SKI5</td><td>(9.59 ± 0.41) × 10<sup>-2</sup></td><td>(1.48 ± 0.07) × 10<sup>-2</sup></td><td>(1.35 ± 0.08) × 10<sup>-2</sup></td></tr><tr><td>SKI6</td><td>(6.02 ± 0.33) × 10<sup>-1</sup></td><td>(7.97 ± 0.58) × 10<sup>-1</sup></td><td>(5.66 ± 0.44) × 10<sup>-1</sup></td></tr><tr><td>SKMP</td><td>(6.59 ± 0.33) × 10<sup>-1</sup></td><td>(8.27 ± 0.69) × 10<sup>-1</sup></td><td>(6.04 ± 0.45) × 10<sup>-1</sup></td></tr><tr><td>SKOP</td><td>(1.47 ± 0.07)</td><td>(8.09 ± 0.36) × 10<sup>-1</sup></td><td>(6.70 ± 0.51) × 10<sup>-1</sup></td></tr><tr><td>SLY</td><td>(2.32 ± 0.12)</td><td>(9.53 ± 0.59) × 10<sup>-1</sup></td><td>(6.93 ± 0.42) × 10<sup>-1</sup></td></tr><tr><td>SLY2</td><td>(2.20 ± 0.11)</td><td>(1.09 ± 0.08)</td><td>(6.55 ± 0.49) × 10<sup>-1</sup></td></tr><tr><td>SLY230A</td><td>(2.19 ± 0.13)</td><td>(8.02 ± 0.50) × 10<sup>-1</sup></td><td>(6.40 ± 0.38) × 10<sup>-1</sup></td></tr><tr><td>SLY9</td><td>(7.64 ± 0.34) × 10<sup>-1</sup></td><td>(9.46 ± 0.49) × 10<sup>-1</sup></td><td>(6.85 ± 0.33) × 10<sup>-1</sup></td></tr><tr><td>WFF1</td><td>(3.75 ± 0.21)</td><td>(4.32 ± 0.25)</td><td>(3.01 ± 0.19)</td></tr></tbody></table>

missing and filled in using the corresponding values of the H3 equation of state model from the same source. For the equations of state from [129], the sound speed becomes superluminal shortly before the maximum mass model is reached, both for the tabulated versions [128] (AP4) or [127] (APR), and the piecewise polytropic fit from [47] (APR4). We therefore use the modified APR4 (denoted APR4_EPP here) from [130], which adds additional segments at high density to ensure causality.
<table><caption>**Table B3.** Bayes factors computed using the RIFT algorithm for all of the equations of state discussed in section 2.4 using the narrow prior discussed in section 2.2. Bayes factors are computed for the costly waveform models discussed in section 2.3. Here the Bayes factors are shown as a ratio between the evidence for the given equation of state model and the corresponding evidence for the hypothesis that both bodies are neutron stars, described by the SLY9 equation of state.</caption><thead><tr><th rowspan="2">Equation of state</th><th colspan="2">RIFT results</th></tr><tr><th>SEOBNRv4T (narrow)</th><th>TEOBResumS (narrow)</th></tr></thead><tbody><tr><td>APR4_EPP</td><td>(2.37 ± 0.06)</td><td>(1.46 ± 0.04)</td></tr><tr><td>BBH</td><td>(1.35 ± 0.03)</td><td>(2.76 ± 0.06) × 10<sup>-1</sup></td></tr><tr><td>BHF_BBB2</td><td>(2.53 ± 0.06)</td><td>(1.38 ± 0.05)</td></tr><tr><td>H4</td><td>(1.06 ± 0.02) × 10<sup>-1</sup></td><td>(8.26 ± 0.18) × 10<sup>-2</sup></td></tr><tr><td>HQC18</td><td>(2.38 ± 0.06)</td><td>(1.46 ± 0.04)</td></tr><tr><td>KDE0V</td><td>(2.36 ± 0.05)</td><td>(1.49 ± 0.05)</td></tr><tr><td>KDE0V1</td><td>(2.30 ± 0.06)</td><td>(1.47 ± 0.04)</td></tr><tr><td>MPA1</td><td>(7.43 ± 0.14) × 10<sup>-1</sup></td><td>(7.64 ± 0.19) × 10<sup>-1</sup></td></tr><tr><td>MS1B_PP</td><td>(2.59 ± 0.06) × 10<sup>-2</sup></td><td>(1.22 ± 0.02) × 10<sup>-2</sup></td></tr><tr><td>MS1_PP</td><td>(6.98 ± 0.19) × 10<sup>-3</sup></td><td>(2.99 ± 0.07) × 10<sup>-3</sup></td></tr><tr><td>RS</td><td>(4.43 ± 0.08) × 10<sup>-1</sup></td><td>(4.61 ± 0.11) × 10<sup>-1</sup></td></tr><tr><td>SK255</td><td>(4.61 ± 0.09) × 10<sup>-1</sup></td><td>(4.90 ± 0.12) × 10<sup>-1</sup></td></tr><tr><td>SK272</td><td>(3.42 ± 0.07) × 10<sup>-1</sup></td><td>(3.41 ± 0.08) × 10<sup>-1</sup></td></tr><tr><td>SKI2</td><td>(1.67 ± 0.03) × 10<sup>-1</sup></td><td>(1.55 ± 0.03) × 10<sup>-1</sup></td></tr><tr><td>SKI3</td><td>(1.64 ± 0.03) × 10<sup>-1</sup></td><td>(1.43 ± 0.03) × 10<sup>-1</sup></td></tr><tr><td>SKI4</td><td>(9.07 ± 0.22) × 10<sup>-1</sup></td><td>(9.27 ± 0.21) × 10<sup>-1</sup></td></tr><tr><td>SKI5</td><td>(5.95 ± 0.13) × 10<sup>-2</sup></td><td>(3.69 ± 0.07) × 10<sup>-2</sup></td></tr><tr><td>SKI6</td><td>(8.05 ± 0.16) × 10<sup>-1</sup></td><td>(8.32 ± 0.21) × 10<sup>-1</sup></td></tr><tr><td>SKMP</td><td>(8.50 ± 0.17) × 10<sup>-1</sup></td><td>(8.73 ± 0.23) × 10<sup>-1</sup></td></tr><tr><td>SKOP</td><td>(1.51 ± 0.03)</td><td>(1.29 ± 0.03)</td></tr><tr><td>SLY</td><td>(1.92 ± 0.04)</td><td>(1.42 ± 0.04)</td></tr><tr><td>SLY2</td><td>(1.96 ± 0.04)</td><td>(1.46 ± 0.04)</td></tr><tr><td>SLY230A</td><td>(1.80 ± 0.04)</td><td>(1.42 ± 0.04)</td></tr><tr><td>SLY9</td><td>(1.00 ± 0.02)</td><td>(1.00 ± 0.02)</td></tr><tr><td>WFF1</td><td>(2.67 ± 0.06)</td><td>(1.16 ± 0.03)</td></tr></tbody></table>

## Appendix B. Bayes factors

Here we present the explicit values for the Bayes factors computed using the nested sampling algorithm implemented in LALINFERENCE and RIFT. We show results for all of the equations of state discussed in section 2.4 for both the wide and narrow priors outlined in section 2.2.
<table><caption>**Table B4.** Bayes factors computed using the RIFT algorithm for all of the equations of state discussed in section 2.4 using the broad prior discussed in section 2.2. Bayes factors are computed for the costly waveform models discussed in section 2.3. Here the Bayes factors are shown as a ratio between the evidence for the given equation of state model and the corresponding evidence for the hypothesis that both bodies are neutron stars, described by the SLY9 equation of state.</caption><thead><tr><th rowspan="2">Equation of state</th><th colspan="2">RIFT results</th></tr><tr><th>SEOBNRv4T (broad)</th><th>TEOBResumS (broad)</th></tr></thead><tbody><tr><td>APR4_EPP</td><td>(2.67 ± 0.08)</td><td>(1.64 ± 0.05)</td></tr><tr><td>BBH</td><td>(1.22 ± 0.04)</td><td>(2.06 ± 0.09) × 10<sup>-1</sup></td></tr><tr><td>BHF_BBB2</td><td>(2.89 ± 0.08)</td><td>(1.65 ± 0.06)</td></tr><tr><td>H4</td><td>(2.33 ± 0.08) × 10<sup>-1</sup></td><td>(1.79 ± 0.06) × 10<sup>-1</sup></td></tr><tr><td>HQC18</td><td>(2.51 ± 0.08)</td><td>(1.69 ± 0.05)</td></tr><tr><td>KDE0V</td><td>(2.50 ± 0.09)</td><td>(1.70 ± 0.05)</td></tr><tr><td>KDE0V1</td><td>(2.39 ± 0.07)</td><td>(1.69 ± 0.06)</td></tr><tr><td>MPA1</td><td>(7.85 ± 0.25) × 10<sup>-1</sup></td><td>(8.06 ± 0.19) × 10<sup>-1</sup></td></tr><tr><td>MS1B_PP</td><td>(4.82 ± 0.18) × 10<sup>-2</sup></td><td>(3.88 ± 0.20) × 10<sup>-2</sup></td></tr><tr><td>MS1_PP</td><td>(5.42 ± 0.23) × 10<sup>-3</sup></td><td>(3.47 ± 0.18) × 10<sup>-3</sup></td></tr><tr><td>RS</td><td>(4.91 ± 0.15) × 10<sup>-1</sup></td><td>(4.54 ± 0.11) × 10<sup>-1</sup></td></tr><tr><td>SK255</td><td>(4.99 ± 0.14) × 10<sup>-1</sup></td><td>(4.66 ± 0.12) × 10<sup>-1</sup></td></tr><tr><td>SK272</td><td>(4.18 ± 0.14) × 10<sup>-1</sup></td><td>(3.63 ± 0.10) × 10<sup>-1</sup></td></tr><tr><td>SKI2</td><td>(2.91 ± 0.09) × 10<sup>-1</sup></td><td>(2.18 ± 0.06) × 10<sup>-1</sup></td></tr><tr><td>SKI3</td><td>(2.78 ± 0.09) × 10<sup>-1</sup></td><td>(2.18 ± 0.07) × 10<sup>-1</sup></td></tr><tr><td>SKI4</td><td>(9.19 ± 0.26) × 10<sup>-1</sup></td><td>(9.53 ± 0.28) × 10<sup>-1</sup></td></tr><tr><td>SKI5</td><td>(1.48 ± 0.05) × 10<sup>-1</sup></td><td>(1.17 ± 0.05) × 10<sup>-1</sup></td></tr><tr><td>SKI6</td><td>(8.22 ± 0.24) × 10<sup>-1</sup></td><td>(8.38 ± 0.21) × 10<sup>-1</sup></td></tr><tr><td>SKMP</td><td>(8.21 ± 0.28) × 10<sup>-1</sup></td><td>(8.36 ± 0.23) × 10<sup>-1</sup></td></tr><tr><td>SKOP</td><td>(1.50 ± 0.05)</td><td>(1.38 ± 0.04)</td></tr><tr><td>SLY</td><td>(2.01 ± 0.06)</td><td>(1.62 ± 0.05)</td></tr><tr><td>SLY2</td><td>(1.97 ± 0.07)</td><td>(1.61 ± 0.05)</td></tr><tr><td>SLY230A</td><td>(1.87 ± 0.06)</td><td>(1.56 ± 0.05)</td></tr><tr><td>SLY9</td><td>(1.00 ± 0.03)</td><td>(1.00 ± 0.03)</td></tr><tr><td>WFF1</td><td>(2.89 ± 0.08)</td><td>(1.20 ± 0.04)</td></tr></tbody></table>

References

[1] Aasi J et al (LIGO Scientific) 2015 Class. Quantum Grav. 32 074001

[2] Acernese F et al (Virgo) 2015 Class. Quantum Grav. 32 024001

[3] Abbott B P et al (Virgo, LIGO Scientific) 2017 Phys. Rev. Lett. 119 161101

[4] Abbott B P et al (Virgo, Fermi-GBM, INTEGRAL, LIGO Scientific) 2017 Astrophys. J. 848 L13

[5] Abbott B P et al (GROND, SALT Group, OzGrav, DFN, INTEGRAL, Virgo, Insight-Hxmt, MAXI Team, Fermi-LAT, J-GEM, RATIR, IceCube, CAASTRO, LWA, ePESSSTO, GRAWITA, RIMAS, SKA South Africa/MeerKAT, H.E.S.S., 1M2H Team, IKI-GW Follow-up, Fermi GBM, Pi of Sky, DWF (Deeper Wider Faster Program), Dark Energy Survey, MASTER, AstroSat Cadmium Zinc Telluride Imager Team, Swift, Pierre Auger, ASKAP, VINROUGE, JAGWAR, Chandra Team at McGill University, TTU-NRAO, GROWTH, AGILE Team, MWA, ATCA, AST3, TOROS, Pan-STARRS, NuSTAR, ATLAS Telescopes, BOOTES, CaltechNRAO, LIGO Scientific, High Time Resolution Universe Survey, Nordic Optical Telescope, Las Cumbres Observatory Group, TZAC Consortium, LOFAR, IPN, DLT40, Texas Tech University, HAWC,
ANTARES, KU, Dark Energy Camera GWEM, CALET, Euro VLBI Team, ALMA) 2017
Astrophys. J. 848 L12

[6] Abbott B P et al (LIGO Scientific, Virgo) 2019 Phys. Rev. X 9 011001

[7] Abbott B P et al (LIGO Scientific, Virgo) 2019 Phys. Rev. X 9 031040

[8] Özel F and Freire P 2016 Ann. Rev. Astron. Astrophys. 54 401-40

[9] Baym G, Hatsuda T, Kojo T, Powell P D, Song Y and Takatsuka T 2018 Rep. Prog. Phys. 81 056902

[10] Flanagan E E and Hinderer T 2008 Phys. Rev. D 77 021502

[11] Poisson E 1998 Phys. Rev. D 57 5287-90

[12] Hinderer T, Lackey B D, Lang R N and Read J S 2010 Phys. Rev. D 81 123016

[13] Vines J, Flanagan E E and Hinderer T 2011 Phys. Rev. D 83 084051

[14] Damour T, Nagar A and Villain L 2012 Phys. Rev. D 85 123007

[15] Baumgarte T W, Shapiro S L and Shibata M 2000 Astrophys. J. 528 L29

[16] Shibata M and Taniguchi K 2006 Phys. Rev. D 73 064027

[17] Baiotti L, Giacomazzo B and Rezzolla L 2008 Phys. Rev. D 78 084033

[18] Hotokezaka K, Kyutoku K, Okawa H, Shibata M and Kiuchi K 2011 Phys. Rev. D 83 124008

[19] Baiotti L and Rezzolla L 2017 Rep. Prog. Phys. 80 096901

[20] Abbott B P, Abbott R, Abbott T D, Acernese F, Ackley K, Adams C, Adams T, Addesso P, Adhikari R X and Adya V B 2019 Astrophys. J. 875 160

[21] Radice D, Perego A, Bernuzzi S and Zhang B 2018 Mon. Not. R. Astron. Soc. 481 3670-82

[22] Abbott B P et al (Virgo, LIGO Scientific) 2017 Astrophys. J. 851 L16

[23] Oliver M, Keitel D, Miller A L, Estelles H and Sintes A M 2019 Mon. Not. R. Astron. Soc. 485 843-50

[24] Baillot d'Etivaux N, Guillot S, Margueron J, Webb N, Catelan M and Reisenegger A 2019 Astrophys. J. (arXiv:1905.01081) accepted

[25] Steiner A W, Heinke C O, Bogdanov S, Li C, Ho W C G, Bahramian A and Han S 2018 Mon. Not. Roy. Astron. Soc. 476 421-35

[26] Antoniadis J et al 2013 Science 340 6131

[27] Margalit B and Metzger B D 2017 Astrophys. J. 850 L19

[28] Bauswein A, Just O, Janka H T and Stergioulas N 2017 Astrophys. J. 850 L34

[29] Ruiz M, Shapiro S L and Tsokaros A 2018 Phys. Rev. D 97 021501

[30] Annala E, Gorda T, Kurkela A and Vuorinen A 2018 Phys. Rev. Lett. 120 172703

[31] Zhou E P, Zhou X and Li A 2018 Phys. Rev. D 97 083015

[32] Fattoyev F J, Piekarewicz J and Horowitz C J 2018 Phys. Rev. Lett. 120 172702

[33] Paschalidis V, Yagi K, Alvarez-Castillo D, Blaschke D B and Sedrakian A 2018 Phys. Rev. D 97 084038

[34] Nandi R and Char P 2018 Astrophys. J. 857 12

[35] Most E R, Weih L R, Rezzolla L and Schaffner-Bielich J 2018 Phys. Rev. Lett. 120 261103

[36] Rezzolla L, Most E R and Weih L R 2018 Astrophys. J. 852 L25

[37] Raithel C, Özel F and Psaltis D 2018 Astrophys. J. 857 L23

[38] De S, Finstad D, Lattimer J M, Brown D A, Berger E and Biwer C M 2018 Phys. Rev. Lett. 121 091102

[39] Abbott B P et al (Virgo, LIGO Scientific) 2018 Phys. Rev. Lett. 121 161101

[40] Landry P and Essick R 2019 Phys. Rev. D 99 084049

[41] Capano C D, Tews I, Brown S M, Margalit B, De S, Kumar S, Brown D A, Krishnan B and Reddy S 2019 (arXiv:1908.10352)

[42] Veitch J et al 2015 Phys. Rev. D 91 042003

[43] Wysocki D, Lange J and O'Shaughnessy R 2019 Phys. Rev. D 100 043012

[44] Del Pozzo W, Li T G F, Agathos M, Van Den Broeck C and Vitale S 2013 Phys. Rev. Lett. 111 071101

[45] Agathos M, Meidam J, Del Pozzo W, Li T G F, Tompitak M, Veitch J, Vitale S and Van Den Broeck C 2015 Phys. Rev. D 92 023012

[46] Chatziioannou K, Yagi K, Klein A, Cornish N and Yunes N 2015 Phys. Rev. D 92 104008

[47] Read J S, Lackey B D, Owen B J and Friedman J L 2009 Phys. Rev. D 79 124032

[48] Yagi K and Yunes N 2017 Phys. Rep. 681 1-72

[49] Chatziioannou K, Haster C J and Zimmerman A 2018 Phys. Rev. D 97 104036

[50] Lindblom L 2018 Phys. Rev. D 97 123019

[51] Carney M F, Wade L E and Irwin B S 2018 Phys. Rev. D 98 063004

[52] Shibata M, Duez M D, Liu Y T, Shapiro S L and Stephens B C 2006 Phys. Rev. Lett. 96 031102

[53] Metzger B D, Quataert E and Thompson T A 2008 Mon. Not. R. Astron. Soc. 385 1455
[54] Shibata M, Fujibayashi S, Hotokezaka K, Kiuchi K, Kyutoku K, Sekiguchi Y and Tanaka M 2017
*Phys. Rev. D* **96** 123012

[55] Veitch J and Vecchio A 2010 *Phys. Rev. D* **81** 062003

[56] Creighton J D E and Anderson W G 2011 *Gravitational-Wave Physics and Astronomy: an Introduction to Theory, Experiment and Data Analysis* (Weinheim: Wiley-VCH Verlag GmbH & Co)

[57] Abbott B P et al (LIGO Scientific, Virgo) 2016 *Phys. Rev. Lett.* **116** 241102

[58] Lattimer J M 2012 *Ann. Rev. Nucl. Part. Sci.* **62** 485–515

[59] Hebeler K, Lattimer J M, Pethick C J and Schwenk A 2013 *Astrophys. J.* **773** 11

[60] Lorimer D R 2008 *Living Rev. Relativ.* **11** 8

[61] Kramer M and Wex N 2009 *Class. Quantum Grav.* **26** 073001

[62] Kiziltan B, Kottas A, De Yoreo M and Thorsett S E 2013 *Astrophys. J.* **778** 66

[63] Lynch R S, Freire P C C, Ransom S M and Jacoby B A 2012 *Astrophys. J.* **745** 109

[64] Özel F, Psaltis D, Narayan R and Villarreal A S 2012 *Astrophys. J.* **757** 55

[65] Martinez J G, Stovall K, Freire P C C, Deneva J S, Jenet F A, McLaughlin M A, Bagchi M, Bates S D and Ridolfi A 2015 *Astrophys. J.* **812** 143

[66] Zhu X, Thrane E, Oslowski S, Levin Y and Lasky P D 2018 *Phys. Rev. D* **98** 043002

[67] Tauris T M et al 2017 *Astrophys. J.* **846** 170

[68] Burgay M et al 2003 *Nature* **426** 531–3

[69] Hessels J W T, Ransom S M, Stairs I H, Freire P C C, Kaspi V M and Camilo F 2006 *Science* **311** 1901–4

[70] Bauswein A and Stergioulas N 2015 *Phys. Rev. D* **91** 124056

[71] Bernuzzi S, Dietrich T and Nagar A 2015 *Phys. Rev. Lett.* **115** 091101

[72] Rezzolla L and Takami K 2016 *Phys. Rev. D* **93** 124051

[73] Cornish N J and Littenberg T B 2015 *Class. Quantum Grav.* **32** 135012

[74] Klimenko S et al 2016 *Phys. Rev. D* **93** 042004

[75] Chatziioannou K, Clark J A, Bauswein A, Millhouse M, Littenberg T B and Cornish N 2017 *Phys. Rev. D* **96** 124035

[76] Torres-Rivas A, Chatziioannou K, Bauswein A and Clark J A 2019 *Phys. Rev. D* **99** 044014

[77] Tolman R C 1939 *Phys. Rev.* **55** 364–73

[78] Oppenheimer J R and Volkoff G M 1939 *Phys. Rev.* **55** 374–81

[79] Damour T and Nagar A 2009 *Phys. Rev. D* **80** 084035

[80] Blanchet L 2014 *Living Rev. Relativ.* **17** 2

[81] Bini D, Damour T and Faye G 2012 *Phys. Rev. D* **85** 124034

[82] Dietrich T, Bernuzzi S and Tichy W 2017 *Phys. Rev. D* **96** 121501

[83] Dietrich T et al 2019 *Phys. Rev. D* **99** 024029

[84] Bohé A et al 2017 *Phys. Rev. D* **95** 044028

[85] Hannam M, Schmidt P, Bohé A, Haegel L, Husa S, Ohme F, Pratten G and Pürrer M 2014 *Phys. Rev. Lett.* **113** 151101

[86] Schmidt P, Ohme F and Hannam M 2015 *Phys. Rev. D* **91** 024043

[87] Husa S, Khan S, Hannam M, Pürrer M, Ohme F, Jiménez Forteza X and Bohé A 2016 *Phys. Rev. D* **93** 044006

[88] Khan S, Husa S, Hannam M, Ohme F, Pürrer M, Jiménez Forteza X and Bohé A 2016 *Phys. Rev. D* **93** 044007

[89] Bohe A, Faye G, Marsat S and Porter E K 2015 *Class. Quantum Grav.* **32** 195010

[90] Hinderer T et al 2016 *Phys. Rev. Lett.* **116** 181101

[91] Steinhoff J, Hinderer T, Buonanno A and Taracchini A 2016 *Phys. Rev. D* **94** 104028

[92] Nagar A et al 2018 *Phys. Rev. D* **98** 104052

[93] Dietrich T and Hinderer T 2017 *Phys. Rev. D* **95** 124006

[94] Nagar A, Messina F, Rettegno P, Bini D, Damour T, Geralico A, Akcay S and Bernuzzi S 2019
*Phys. Rev. D* **99** 044007

[95] Akcay S, Bernuzzi S, Messina F, Nagar A, Ortiz N and Rettegno P 2019 *Phys. Rev. D* **99** 044051

[96] Wade L, Creighton J D E, Ochsner E, Lackey B D, Farr B F, Littenberg T B and Raymond V 2014 *Phys. Rev. D* **89** 103012

[97] Dudi R, Pannarale F, Dietrich T, Hannam M, Bernuzzi S, Ohme F and Brügmann B 2018 *Phys. Rev. D* **98** 084061

[98] Wiringa R B, Fiks V and Fabrocini A 1988 *Phys. Rev. C* **38** 1010–37

[99] Cromartie H T et al 2019 *Nature Astron.* (https://doi.org/10.1038/s41550-019-0880-2)

[100] Bludman S A and Ruderman M A 1970 *Phys. Rev. D* **1** 3243–6

[101] Lattimer J M, Prakash M, Masak D and Yahil A 1990 *Astrophys. J.* **355** 241–54
[102] Cahillane C et al 2017 Phys. Rev. D 96 102001

[103] Viets A et al 2018 Class. Quantum Grav. 35 095015

[104] Driggers J C et al 2019 Phys. Rev. D 99 042001

[105] Pankow C et al 2018 Phys. Rev. D 98 084016

[106] Acernese F et al 2018 Class. Quantum Grav. 35 205004

[107] Skilling J 2006 Bayesian Anal. 1 833–59

[108] Lartillot N and Philippe H 2006 Syst. Biol. 55 195–207

[109] Abbott B P et al (LIGO Scientific, Virgo) 2019 Phys. Rev. Lett. 122 061104

[110] Samajdar A and Dietrich T 2018 Phys. Rev. D 98 124030

[111] Kass R E and Raftery A E 1995 J. Am. Statist. Assoc. 90 773–95

[112] Ruiz M and Shapiro S L 2017 Phys. Rev. D 96 084063

[113] Dai Z G and Lu T 1998 Phys. Rev. Lett. 81 4301–4

[114] Zhang B and Meszaros P 2001 Astrophys. J. 552 L35–8

[115] Kasen D, Metzger B, Barnes J, Quataert E and Ramirez-Ruiz E 2017 Nature 551 80

[116] Villar V A et al 2017 Astrophys. J. Lett. 851 L21

[117] Cowperthwaite P S et al 2017 Astrophys. J. 848 L17

[118] Abbott B P et al 2017 Astrophys. J. 850 L39

[119] Metzger B D 2017 (arXiv:1710.05931)

[120] Kawaguchi K, Shibata M and Tanaka M 2018 Astrophys. J. 865 L21

[121] Yu Y W, Liu L D and Dai Z G 2018 Astrophys. J. 861 114

[122] Kiuchi K, Kyutoku K, Shibata M and Taniguchi K 2019 Astrophys. J. 876 L31

[123] Radice D, Perego A, Zappa F and Bernuzzi S 2018 Astrophys. J. 852 L29

[124] Dietrich T and Ujevic M 2017 Class. Quantum Grav. 34 105014

[125] Shibata M, Zhou E, Kiuchi K and Fujibayashi S 2019 Phys. Rev. D 100 023015

[126] Stergioulas N and Friedman J L 1995 Astrophys. J. 444 306–11

[127] Typel S, Oertel M and Klähn T EoS catalog https://compose.obspm.fr (tables retrieved August 2018)

[128] Özel F EoS catalog http://xtreme.as.arizona.edu/NeutronStars/data/eos_tables.tar (tables retrieved 17 August 2018)

[129] Akmal A, Pandharipande V R and Ravenhall D G 1998 Phys. Rev. C 58 1804–28

[130] Endrizzi A, Ciolfi R, Giacomazzo B, Kastaun W and Kawamura T 2016 Class. Quantum Grav. 33 164001

[131] Baldo M, Bombaci I and Burgio G F 1997 Astron. Astrophys. 328 274–82 (https://ui.adsabs. harvard.edu/#abs/1997A&A...328..274B/abstract)

[132] Gulminelli F and Raduta A R 2015 Phys. Rev. C 92 055803

[133] Danielewicz P and Lee J 2009 Nucl. Phys. A 818 36–96

[134] Agrawal B K, Shlomo S and Au V K 2005 Phys. Rev. C 72 014310

[135] Reinhard P G, Dean D J, Nazarewicz W, Dobaczewski J, Maruhn J A and Strayer M R 1999 Phys. Rev. C 60 014316

[136] Lackey B D, Nayyar M and Owen B J 2006 Phys. Rev. D 73 024021

[137] Douchin F and Haensel P 2001 Astron. Astrophys. 380 151–67

[138] Chabanat E, Bonche P, Haensel P, Meyer J and Schaeffer R 1997 Nuc. Phys. A 627 710–46

[139] Bennour L, Heenen P H, Bonche P, Dobaczewski J and Flocard H 1989 Phys. Rev. C 40 2834–9

[140] Friedrich J and Reinhard P G 1986 Phys. Rev. C 33 335–51

[141] Agrawal B K, Shlomo S and Kim Au V 2003 Phys. Rev. C 68 031304

[142] Reinhard P G and Flocard H 1995 Nucl. Phys. A 584 467–88

[143] Nazarewicz W, Dobaczewski J, Werner T R, Maruhn J A, Reinhard P G, Rutz K, Chinn C R, Umar A S and Strayer M R 1996 Phys. Rev. C 53 740–51

[144] Müther H, Prakash M and Ainsworth T L 1987 Phys. Lett. B 199 469–74

[145] Müller H and Serot B D 1996 Nucl. Phys. A 606 508–37