Hubble constant and nuclear equation of state from kilonova spectro-photometric light curves

The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy r -process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the 10.4m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600 − 10000Å spectral range. Here, we study the most representative cases regarding the scientiﬁc interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: (a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identiﬁed and (b) measure the Hubble constant, H 0 , with precision improved by up to 40%, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance a ﬀ ect the precision and estimated value of H 0 .


Introduction
The recent era of multimessenger astrophysics has opened up new possibilities to disentangle the properties of matter and probe the evolution of the Universe itself. One of the most prominent multimessenger observations was the joint detection of the gravitational wave (GW) event GW170817 on 2017 August 17.53 UT. Advanced LIGO/Virgo made the first detection of GWs from a binary neutron star (BNS) merger , with a simultaneous short gamma ray burst (GRB) detected by Fermi and INTEGRAL (GRB 170817A, Goldstein et al. 2017;Blackburn et al. 2017). At ∼ 0.5 days after the GW signal, an electromagnetic (EM) counterpart -the first for any GW source -also known as a kilonova (KN), was identified within the ∼30 deg 2 localization region and named AT 2017gfo. The source, located in the galaxy NGC 4993, was compatible with a BNS merger event at 40 Mpc Hjorth et al. 2017). The KN thermal emission observed in AT 2017gfo was produced by radioactive decay of neutron rich matter from the material ejected during and after the merger.
In this paper we focus on compact star systems, binaries of neutron stars (NSs) that merge after the final evolutionary stages of inspiral coalescence. In these violent events, where gravity is in the strong regime, the merging process leads to the formation of a highly massive object, producing a broad multimessenger signal: GWs, EM waves, neutrinos (Kyutoku et al. 2018) and even cosmic rays (Kimura et al. 2018). The dynamics of this kind of merger events is not yet well understood, relying on computational simulations that include detailed microphysics input at an increasing level of accuracy (Dietrich & Ujevic 2017). More importantly, neutrinos and weak interactions determine the composition and lepton fraction and ultimately the opacity of the expanding ejecta. This response modulates the nontrivial nonhomogeneous blue and red emission in the KN, as computational simulations have shown (see Metzger 2019 and Nakar 2020 for recent reviews). The high-density material in BNS mergers is heated and thermal neutrinos with energies in excess of a few MeV are expected (Sekiguchi et al. 2011). Direct detection of thermal neutrinos by existing water-tank detectors such as Super-Kamiokande (Abe et al. 2018) is quite challenging for mergers at a distance ≳100 Mpc; however, it is likely that they will contribute to the diffuse neutrino background. Highenergy neutrinos associated with the GW170817 merger were also searched for by ANTARES, IceCube, and Pierre Auger Observatories, but no significant signal was observed (Albert et al. 2017).
The KN coincident with GW170817 showed a rapidly fading EM transient in the optical and infrared bands. The term kilonova was historically coined from the perception of its brightness being a thousand times larger than a nova (Li & Paczyński 1998;Rosswog 2005;Metzger et al. 2010). We discuss in this work that, even if far from such a standardized luminosity, clear dependencies on physical BNS quantities can be extracted. For AT 2017gfo, the observed spectral evolution is suggested to arise from a mixed composition of r−process elements in the ejected material (Gillanders et al. 2022); while early-time spectra are consistent with light r−process elements (nuclear masses A ∼ 90-140), later spectra require intermediate composition, producing even heavier elements such as lanthanides. An alternative scenario, such as dust formation in the KN that could also explain such a blue-red evolution, has been ruled out (Gall et al. 2017). In addition, a broad feature observed at ∼ 0.7 − 1 µm has been interpreted as due to Sr ii (Watson et al. 2019;Domoto et al. 2021;Gillanders et al. 2022). Since then, no other KN has been firmly detected (see e.g., Coughlin 2020 and references therein for a summary of the follow-up efforts during LIGO and Virgo's third observing campaign). Another merger candidate is linked to GW190814, located farther away at ∼ 240 Mpc, whose origin remains uncertain. However, it lacks of any EM counterpart and could likely be caused by a binary black hole merger (Essick & Landry 2020;Tews et al. 2021). GW190425 was the first BNS candidate in O3. The LALInference localization pipeline (Veitch et al. 2015) provided a 90% credible region of 7461 deg 2 , estimated luminosity distance of 156 ± 41 Mpc, much larger than that for GW170817, which complicated the search campaign resulting in no EM counterpart identified. The inferred merger parameters from the GW analysis published in Abbott et al. (2020a) for chirp mass 1.44 +0.02 −0.02 M ⊙ are the following. Total mass M 1 + M 2 ≃ 3.4 +0.3 −0.1 M ⊙ , with M 1 ∈ (1.62, 1.88)M ⊙ , M 2 ∈ (1.45, 1.69)M ⊙ , andΛ ≤ 600 for low spin prior, or M 1 ∈ (1.61, 2.52)M ⊙ , M 2 ∈ (1.12, 1.68)M ⊙ , and Λ ≤ 1100 for high-spin prior. Note that for the other confirmed BNS event, GW170817, the inferred total mass is M 1 + M 2 ≃ 2.73 +0.04 −0.01 M ⊙ with M 1 ∈ (1.36, 1.60)M ⊙ , M 2 ∈ (1.16, 1.36)M ⊙ , andΛ = 300 +420 −230 for low-spin prior (Abbott et al. 2019). Detecting KNe as EM counterparts to BNS mergers can shed light on different relevant areas of physics. For instance, it is expected that tighter constraints on existing cosmological models, formation of elements, the internal composition of the individual NSs, and even more fundamental physics can benefit from present and future multimessenger observations.
One particular example is the constraints on the Equation of State (EoS) as it describes the interior of the NSs. KN observables such as spectra and light curves are controlled by the properties of the material ejected during the merger, including mass, expansion velocity, geometry and composition (setting the opacity). These, in turn, depend on the properties of the merging system like e.g. NS and BH mass/spin and NS radius or the ratio of both, i.e. compactness. Constraints on NS masses and radii have been provided in a series of world-wide experimental efforts. In particular, recent observations from the Neutron star Interior Composition ExploreR (NICER) mission suggest that the EoS of NS might be rather stiff (Miller et al. 2021;Raaijmakers et al. 2021;Riley et al. 2021), indicating that many NS mergers should eject a significant amount of material and produce bright KNe. From the population studies of Galactic double NSs (DNSs) it appears that ≃ 98% of all merging DNSs will have a mass ratio, q, close to unity, q < 1.1 (Zhang et al. 2011). Recent studies find that the weighted mean masses of the primary and companion stars are respectively (1.439 ± 0.036)M ⊙ and (1.239 ± 0.020)M ⊙ (Yi-yan et al. 2017). In addition, the surface magnetic field strength in the primary stars of the DNSs is B ∼ 10 10 G, and the spin period is P ∼ 50 ms. Due to the fact that the NS in the BNS events under scrutiny will share these general characteristics it is important to see what information the EM counterparts can yield.
One other example is linked to the value of the Hubble constant H 0 , especially interesting in the context of the present tension in H 0 measurements from the Early (Planck Collaboration et al. 2020) and the Late (Riess et al. 2022) Universe (see e.g. Verde et al. 2019 for a recent review). The simultaneous detection of GWs and EM radiation from the same source led to independent measurements of the distance and redshift of the source, thus providing an estimate of H 0 that is independent of cosmological model assumptions as well as the local distance ladder . Using GW as "standard sirens" (Holz & Hughes 2005), this approach holds promise to arbitrate the existing tension in H 0 measurements. However, the well-known degeneracy in the GW signal between the luminosity distance and inclination angle translates into large uncertainties on H 0 , even when an EM counterpart is identified . Independent estimates of the viewing angle of GW170817 from either the GRB afterglow (e.g. Guidorzi et al. 2017), superluminal motion of the GRB jet (Hotokezaka et al. 2019) or the KN (Dhawan et al. 2020;Coughlin et al. 2020a;Dietrich et al. 2020) have been used to reduce the degeneracy between inclination and luminosity distance and therefore the uncertainties on H 0 (see Bulla et al. 2022 for a recent review and e.g. Nakar &Piran 2021 andHeinzel et al. 2021 for possible systematics introduced by these approaches). Given the rapid evolution of KNe, and their intrinsic faintness, a prompt activation of groundbased telescopes is needed to get spectral information at their early epochs.
In addition to currently operating spectrographs on 8-10 meter class telescopes, a new facility will be promptly available to investigate the nature of KNe. Here we present and discuss the role that MAAT 1 (Mirror-slicer Array for Astronomical Transients) will have on KN science. MAAT is a new Integral Field Unit (IFU), based on image slicers, for the OSIRIS spectrograph on the 10.4-m telescope on the Gran Telescopio de CANARIAS (GTC, La Palma), which is planned to being operations in autumn 2023 (Prada et al. 2020). Thanks to the pointing accuracy and the large collecting area of the GTC telescope, MAAT will permit to obtain a spectrum of a KN candidate in less than five minutes from its alert activation. The integral-field spectroscopy capabilities of MAAT, combined with its field-of-view (8.5 ′′ × 12.0 ′′ ) will allow to collect all the photons incoming from KN emission over the wavelength range 3600 − 10000 Å, including its immediate galaxy host environment, avoiding any possible flux losses due to the narrow slit aperture in long-slit spectrographs. Hence, MAAT is an unique capability in 8-10m class telescopes world-wide that will allow to perform absolute spectro-photometry of KNe, and any other interesting transients, see (Prada et al. 2020).
In this work we perform a detailed analysis on how different scenarios involving BNS could help understand not only the physics governing these events, including the nuclear matter EoS, but also how fundamental parameters in cosmology such as the Hubble constant, H 0 , can benefit from complementary EM data obtained with MAAT. For this task we use the output of radiative transfer simulations of the KN yield in the UV-IR range for selected cases of BNS mergers as we discuss below.
The structure of this paper is as follows. In Section 2, we discuss the expected features of future KNe observations and existing rates of GW detections linked to BNS rates. We also comment on the observational strategy to optimally capture the light curve and spectral features providing key information on the physical parameters involved. In Section 3, we explain the subset of adopted EoSs in detail and set up the scenario for the analysis we later perform based on a sample of 23 KN runs and their emission light curves obtained with POSSIS, a radiative transfer Monte Carlo code. We also analyze the impact of the EoS on ejecta magnitudes in light of additional complementary probes from GW observations. In Section 4, we analyze and discuss the shape of KN light curves considering aspects such as inclination angles modulating the measured signal. In Section 5, we tackle the H 0 determination and sources of uncertainty providing some paradigmatic examples of MAAT performance to this task. Finally, in Section 6, we present a discussion and draw our conclusions.

Observations and binary neutron star merger rates
In the BNS scenario, the understanding of the impact of the ejecta on the KN brightness and its usability are crucial to reduce the uncertainties on the Hubble constant of EoS itself. To what extent these light curves can be standardized (Kashyap et al. 2019;Coughlin et al. 2020b) and many of the aspects of the modelling of the emission (Metzger 2019;Radice et al. 2018) is still a matter of debate. The KN emission is expected to be preceded by GWs detected during the subsequent observing runs of the LIGO/Virgo/KAGRA (LVK) collaboration. Second generation (2G) Advanced LIGO (Aasi et al. 2015), Advanced Virgo (Acernese et al. 2015) and KAGRA, and the proposed Einstein telescope (ET) interferometers (Maggiore et al. 2020) will be crucial in providing the initial alert for newly-discovered BNS systems. For the incoming run O4, LIGO is expected to detect BNS mergers up to 160-190 Mpc, whereas Virgo sensitivity limits to 80-115 Mpc 2 . KAGRA will start observing with a very limited sensivitity of 1 Mpc, which is expected to improve to ∼ 10 Mpc towards the end of the run. For the following O5 science run, beyond 2028, LIGO projects a sensitivity goal of 325 Mpc while Virgo and KAGRA project a target sensitivity of 260 Mpc and 128 Mpc, respectively. In addition, ET will detect binary systems containing NSs up to redshifts corresponding to the peak of the cosmic star-formation rate. This will allow the study of the broad field of formation, evolution and physics of NSs in connection with KNe and short GRBs, along with the star formation history and the chemical evolution of the Universe.
The detection of a short GRB event, coincident within the error box localization regions provided by the above GW detectors, will further support the nature of the detection and constrain further those error boxes where candidates will be searched, similarly to what happened for the AT 2017gfo event (Goldstein et al. 2017). Recent advanced KN population synthesis models predict a rate of joint GW and short GRB detection from a BNS merger between ∼ 1-6 events per year (Colombo et al. 2022;Patricelli et al. 2022). During the LIGO-Virgo (LVC) O1, O2 and O3 runs, around ∼50 candidates GW events were reported, all regarding detections of compact objects binary mergers (Abbott et al. 2021). One key aspect when considering the statistics available for any further study is the local merger rate density. The BNS rate is estimated to be 10−1700 Gpc −3 yr −1 as reported in the latest catalog from LVC (GWTC-3, The LIGO Scientific Collaboration et al. 2021). From the analysis of GW170817, a BNS rate of 110 − 3840 Gpc −3 yr −1 was derived (Abbott et al. 2021). More specifically, the number of BNS detections in O4 is expected to be 10 +52 −10 (Abbott et al. 2020b). This number depends on several assumptions on BNS physical properties and detector characteristics. Using simulations of BNS mergers based on the results obtained during the last observing O3 run, Petrov et al. (2022) used different constraints on the minimum network signal-to-noise ratio (S/N) threshold, S/N > 8, for detection of a BNS coalescence, and on the number of interferometers, only one, that would provide a detection of a BNS merger. In Abbott et al. (2020b) a network SNR>12 was used and a detection in at least two interferometers of the LIGO-Virgo collaboration was assumed. These different assumption led to a much better agreement with the results of the O3 run, and provide for the incoming O4 run an expected annual rate of BNS mergers of 34 +78 −25 events, while for the following O5 run an annual rate of 190 +410 −130 events. .
The asymmetry in the NS masses and chirp mass in the BNS will also make an impact on the reachable depth, for example for O4 the LIGO range for BNS with equal mass ratio q = 1 and individual NS masses M 1 = M 2 = 1.4M ⊙ is ∼ 190 Mpc while for BNS having q 1 and M 1 × M 2 = 1.4M ⊙ × 2.0M ⊙ one has ∼ 220 Mpc (Abbott et al. 2020b).
However, once the detection by GW interferometers marks the occurrence of a BNS merger, the identification of its EM counterpart is more challenging. In order to capture the IR-optical-UV counterpart of a candidate BNS merger, an alert system exists. It consists in a communication node, the Gravitational-wave Candidate Event Database, called GraceDB 3 which connects LIGO and Virgo analysis pipelines and sends an alert to astronomers when a promising GW signal has been detected and located in the sky within a given region. That initiates a search campaign with telescopes to point in the highlighted sky area, trying to capture the EM counterpart before it has faded away and ideally before it peaks. From previous events, such as GW170817, the community has gained valuable experience for those to come in the future. Note that this GW event was localized to a sky area of around 30 deg 2 within ∼ 10 hours from the first alert. For the follow-up of a posterior event, GW190814, it was narrowed to about 2 hours (de Wet et al. 2021) and sky localization for this event was ameliorated to 18.5 deg 2 ). In the same line, for O3 and O4 the localizations of BNS area were expected to be 270 +34 −20 and 33 +5 −5 (in units of deg 2 , 90% c.r.) respectively, noting that for the later value the Hanford-Livingston-KAGRA-Virgo (HLKV) is included. In reality, the localizations during O3 were around 1000 to 10000 deg 2 due to low significant trigger alerts and single detector detections, and at distances ∼5 times farther than GW170817 (Kasliwal et al. 2020).
As current KN emission light curves obtained by simulations have shown, the brightness in the KN bolometric luminosity curve peaks at a time, t peak , ranging between tenths of a day and a few days from the merger, depending on intrinsic physical parameters (ejected mass, expansion velocity, composition). The strategy to perform observations of a given KN candidate could be summarized as follows. As a first trigger alert will come from LVK with an estimate of distance and a sky localization area, an immediate search of a possible optical counterpart will begin. An estimate based on the two dimensional projection plus the distance uncertainty yields a three dimensional error box that could be set in order to localize the galaxies of interest for the EM follow-up (up to ∼300 Mpc) in existing catalogs such as the GLADE catalog 4 . Thereafter, a few candidates observed within the LVK error region and spatially close to galaxies having distances within the range values reported by the BNS-GW merger signal might be detected by larger surveys such as the Zwicky Transient Facility (ZTF, Bellm et al. 2019), the Pan-STARRS (Chambers et al. 2016), the Asteroid Terrestrial-impact Last Alert System (ATLAS, Tonry et al. 2018), the Young supernova Experiment (YSE, Jones et al. 2021), or smaller size optical telescopes network such as the Gravitational wave Optical Transient Observatory (GOTO, Dyer et al. 2020), the Global Rapid Advanced Network Devoted to the multimessenger Addicts (GRANDMA, Antier et al. (2020)), the All Sky Automated Survey for SuperNovae (ASAS-SN, Shappee et al. (2014)), in addition to a multitude of other facilities, and then communicated to the astronomical community through Gamma-ray Coordinate Network (GCN) Notices and/or alert brokers, establishing the appearance of a potentially new KN candidate associated with the GW detection. It is at this stage, when a recognition of the nature of the candidate (through its spectrum) is needed, that it is expected that MAAT will join the observation.
Additional transient sources such as intrinsically faint (M peak ∼ −17 → −15 mag) supernovae Type I and II will appear quite differently, see Fig. 1, where we compare spectra as if they were observed by MAAT/GTC using the grism R500 and an exposure time of 300s (corresponding to the setup used for the KN search campaign, immediately following the GW detection alert) of the fast declining type I SNe SN 2012hn (at 4 days before maximum Valenti et al. 2014) and the faint type IIb SN 2019bkc (before peak luminosity Izzo et al. 2022) with the modeled KN spectra of run 7 (see Table 1) at face on and edge on inclinations. All the spectra displayed in Fig. 1 have been also simulated as if they had the same observee magnitude, namely V = 21.5 mag. MAAT pointing will be characterized by a precision better than ∼ 1 ′′ .

Emission patterns and photometry
In order to achieve spectro-photometry with accuracy better than 3% in the entire 3600 − 10000 Å wavelength range available to MAAT, it will be necessary to observe spectro-photometric standard stars at different airmass. Due to the slicer design of MAAT, there will be no slit losses due to seeing, atmospheric dispersion, or any other such effects that are a handicap to longslit observations. By performing multiepoch observations, additional features and their evolution can be observed that will facilitate the source/event identification. ; 2) the faint type IIb SN 2022ngb before peak luminosity (Izzo et al. 2022); 3) the modeled KN spectra of run7 viewed from edge on (cos θ = 0), 4) and face on (cos θ = 1) inclinations, as if they have been observed with MAAT/GTC using the grism R500 and an exposure time of 300s, and being characterized by an observed magnitude of V = 21.5 mag. Figure 3 shows the bolometric luminosity L bol as a function of time obtained from two simulations concerning run 1 with EoS DD2 (blue line) and run 7 with EoS LS220 (red line), see Table 1 and Section 3 for more details. Solid, dashed and dotted lines refer to polar angle inclinations θ corresponding to cos θ = 1, 0.7, 0, respectively. As mentioned before, the observed light curves depend on intrinsic properties of the BNS and observing conditions. As for the former, the role of the EoS is of key importance, as it mostly determines matter interaction and governs the amount of mass that is ejected during the tidal forces originated in the merger event. As shown in the figure, the inclination angle θ measured as a polar angle would allow to discriminate the EoS as emitted luminosities can differ a factor of 10 for selected epochs. From Fig. 3 it is also clear that for a polar observer, luminosity will be larger as compared to any other direction.
In order to characterize the optical emission with MAAT, the observational strategy will allow, provided optimal exposure times, to determine the spectral continuum, by fitting to a black body spectrum. Thus emission and absorption patterns due to the mixture of elements present in the ejected mass shells will show up. In the first stages of the KN emission, the dense environment Simulated KN spectral emission for an event similar to AT 2017gfo at ∼ 1 day as observed with MAAT using the R1000B and R1000R gratings and a total exposure time of 1800s. The top panel shows the resulted observed spectra (grey) compared to the theoretical input from KN run 1 simulation (black) with EoS DD2, see Table 1, assuming similar parameters to those associated to GW170817 at 40 Mpc. The bottom panel shows the S/N obtained for the simulated event for gratings R1000B (blue) and R1000R (red). Note that the top spectrum is the combination of the flux captured by both red (R) and blue (B) R1000 gratings. Fig. 3. Luminosity as a function of time obtained from two selected simulations concerning run 1 with EoS DD2 (blue line) and run 2 with EoS LS220 (red line), see Table 1. Solid, dashed and dotted lines refer to inclination angles cos θ = 1, 0.7, 0, respectively. will facilitate the appearance of blue-shifted absorption dips on top of the spectral continuum. Following the case of AT 2017gfo, for similar events, a Sr ii triplet from neutron capture should appear toward the near-IR wavelength range (Watson et al. 2019;Domoto et al. 2021;Gillanders et al. 2022). In this line, the analysis of the P-Cygni profiles will allow for the determination of Absolute flux obtained at t = 0.5, 1, 5 days after the merger as seen by MAAT using R1000B and R1000R gratings. We plot run 1 (EoS DD2) for a polar and equatorial observer (dark and pale blue lines). Same for run 7 (EoS LS220) for a polar and equatorial observer (red and orange lines). A distance of 40 Mpc is assumed.
the expanding velocity of the KN ejecta, as well as for abundance studies with spectral synthesis codes. It is important to remark that deviations from a simplified spherical geometry are likely to appear. In addition, a shock heated and disk wind structure form determining the blue KN component, whose maximum radiation peaks in 1−2 days after the merger. In the same way, the red component will show up from days to weeks. Therefore a continuous and careful monitoring is crucial to determine the peculiarities of the EM emission event and extract meaningful accurate data. Other features are also expected from different elements appearing in the neutron rich debris of the BNS (Gillanders et al. 2022). Figure 4 shows the absolute flux F λ at three epochs 0.5, 1, 5 d (upper, middle and bottom panels) after the merger, as obtained by simulating the KN emission from a BNS merger at 40 Mpc with mass ratio q = M 1 /M 2 = 1 and chirp mass For this particular case M 1 = M 2 = 1.4 M ⊙ . They correspond to run 1 with EoS DD2 with polar and equatorial observing angles (dark and pale blue lines) and run 7 with EoS LS220 (red and orange lines), see Table 1 and Section 3. The exposure time with MAAT has been set to t exp = 1800 s. We can see that the effect of the EoS is most clearly seen at early times. We detail the role of the selected runs and EoS in Section 3. The viewing angle dependence must be carefully addressed, since it largely impacts the measured spectra (see below).

Observation plans
We now detail how in the first phase of the identification of the KN emission, a series of low-resolution (R500B and/or R500R grisms, R = 860-940) exposures will be used. With this setup, we can maximize both the telescope time and the number of candidates to check. This instrumental configuration will allow us to get spectra with a sufficient S/N (S/N ≳ 10) to classify possible KN candidates up to r ∼ 21.5 mag with a single exposure of 300s. The total amount of telescope time will be ∼ 1 hour to classify about ten candidates. Once, and if, a KN is identified by our or other programs, a monitoring program will be started with a daily cadence, using higher resolution grisms (e.g. R1000B/R, R = 1630-1800) to trace the spectral continuum and identify spectral features attributable to newly-synthesised r−process elements (Smartt et al. 2017;Watson et al. 2019). In Fig. 2, we show the results of our simulations for a KN event similar to AT 2017gfo observed with MAAT using the R1000 grisms, ∼ 1 day after the GW signal detection.
Although for some of the BNS events, an EM counterpart is expected be correlated with the merger event, as for AT 2017gfo and GW170817, for other events such as GW190425 a counterpart was not detected, increasing the uncertainty about its dynamics and nature (Abbott et al. 2020a). The emission of associated GRBs could be observed with current orbiting detectors on-board satellites such as Fermi, Swift, INTEGRAL and the InterPlanetary Network (IPN), a set of spacecrafts with gammaray detectors. This enables us to find out where on the sky an associated short GRB occurred using triangulation techniques (Hurley et al. 2013). It may also happen that, to our regret, we may not observe any associated GRB, as only those beamed towards Earth are detectable. Placing the constraint that a jet is only launched if the mass of the remnant, M rem , relates to the static TOV solution, M TOV , as M rem ≳ 1.2 × M TOV , leads different suites of combined models to predict a BNS jet-launching fraction of f s,BNS < 0.6 (Sarin et al. 2022).

Simulating KNe
During a BNS merger, two NSs collide to finally produce thermal EM emission that we measure as a KN. Originally each NS is born as a hot lepton rich dense object evolving into a cold neutron rich one. Spinning magnetized stars are indeed possible and have been simulated (Bocquet et al. 1995), although we work under the premise of low-spin zero magnetic field for them. The EoS is a key ingredient when describing the behavior of the matter inside NSs. The internal radial structure in the spherical static (non rotating) case is obtained after solving the Tolman-Oppenheimer-Volkoff (TOV) equations providing the mass-radius relationship M(R). Complementary measurements indicate NSs have a radius R ≲ 14 km and M ≲ 2.2 M ⊙ , being the recently detected MSP J0740 +6620 the most massive NS to date with a value 2.14 +0.10 −0.09 M ⊙ (Cromartie et al. 2019). To study the impact on the H 0 values and description of nuclear matter EoS set by the KN physics, we have selected a set of 23 simulations in the literature (see Table 1) for a reduced number of EoS. In particular, we select three different EoS: DD2 (Typel et al. 2010), LS220 (Lattimer & Swesty 1991) and SFHo (Steiner et al. 2013). We consider the works of Radice et al. (2018) and Nedora et al. (2021), with models in the latter tailored to GW170817. Note that the simulation of the outcome of these transients relies on all fundamental interactions known, i.e. gravitational (infall and large tides with object deformation), electroweak physics (through the neutrino transport and nucleosynthesis yields) and hadronic (nuclear matter content of the NSs as described by the underlying EoS and fragmentation of ejecta and decay). Therefore, we choose these simulations as they include a full neutrino transport and multidimensionality in order to see how the different ion population peaks and fractions appear in the spectra, for a discussion see Tanaka (2016). It is important to note that some of these works miss a detailed treatment of neutrino emission/absorption, which may impact the ejecta masses in a systematic way, e.g. underestimating their value, and thus the KN light curve (Radice et al. 2018). The set of EoS used in our work is able to produce NSs with masses below the maximum static or TOV mass for the GWTC-3 analysis . In this same study, when focusing on BNS binaries a broad, relatively flat NS mass distribution extending from 1.2 to 2.0M ⊙ is obtained. Thus in this sense our subset of events can be safely considered to belong to this distribution.
The situation regarding our choice of EoS is informed in this same sense by recent results obtained from Bayesian model selection on a wide variety of models using multimessenger observations. In Fig. 2 in  the Bayes factors for the narrow prior results using different waveform approximants for GW170817 found that the EoS yielding excessively large tidal deformabilities are disfavored by the data. Therefore they are disfavored with a Bayes factor that is much smaller when compared to the most favored EoS models. However, these authors find that only a few of the EoS used in the analysis can be ruled out as most of them have comparable evidences within an order of magnitude. They find this result by sorting the EoS by the compactness, C, of the neutron star at a fiducial mass of 1.36M ⊙ . Therefore EoS producing more compact stars than those predicted for H4 with C = 0.148 are allowed, as it is the case for DD2, SFHo, LS220, with C = 0.152, 0.168 and 0.158, respectively. More in this line, Biswas (2022) used mass and tidal deformability measurements from two BNS events, GW170817 and GW190425, and simultaneous mass-radius measurements of PSR J0030+0451 and PSR J0740+6620 by the NICER collaboration, being able to to rule out different variants of the MS1 family, SKI5, H4, and WFF1 EoS decisively among a set of 31 EoS. These EoS were found to be either extremely stiff or soft.
In Table 1, we list some parameters that are key to understand the KN emission light curves as obtained from BNS merger simulations. We take q and M as the mass ratio and chirp mass, respectively, q = M 1 /M 2 , q ≥ 1, with individual NS masses M 1 , M 2 and chirp mass M chirp ≡ M = (M 1 M 2 ) 3/5 (M 1 +M 2 ) 1/5 . Extrinsic BNS parameters such as observing angle θ, q and M chirp are accessible through GW observation as well as the redshift z and luminosity distance D L . Intrinsic parameters are related to properties of the different contributions to the ejected material, in particular the mass from a dynamical ejecta component M dyn and that from a wind component M wind ejected from an accretion disk of mass M disk created around the merger remnant. The ejecta are described by the average velocity, ⟨v e j ⟩, and the average electron fraction ⟨Y e ⟩. In this scenario, matter is neutron rich (as neutrinos further deleptonize matter) so that ⟨Y e ⟩ < 0.5.
For each of the models in Table 1, KN spectra are extracted using the 3D Monte Carlo radiative transfer code POS-SIS (Bulla 2019). The input geometries are characterized by a spherical disk-wind ejecta component at relatively low veloci-  Fujibayashi et al. 2020), with numbers for the latter varying from roughly 20 to 40%. For each model, here we fix the electron fraction to Y e,wind = 0.3 and the wind mass to an average value of 30% of the disk mass M disk , although future studies should explore the impact of these assumptions on the KN observables. The density profile in the disk-wind component is assumed to scale with radius r as ∝ r −3 . For the dynamical ejecta component, we assume an outflow extending from v dyn min = 0.1 c to v dyn max such that the mass-averaged velocity is the one required for each model, ⟨v dyn ⟩. We set the mass to the specific value for each model and assume the Y e,dyn distribution to have an angular dependence of the form Y e ∝ A+B cos 2 θ (consistent with simulations in Radice et al. 2018, C. Setzer, private communication) where A and B are scaling constants chosen to match the desired ⟨Y e ⟩ ≡ ⟨Y e,dyn ⟩ in each model. A broken power law is chosen for the density profile in the dynamical ejecta, with density scaling with radial distance as ∝ r −4 below and as ∝ r −8 above 0.4 c. The density profile in the dynamical ejecta component is modeled with an angular dependence ∝ sin 2 θ that is found in numerical-relativity simulations (Perego et al. 2017;Radice et al. 2018) and has been shown to provide a better fit to AT 2017gfo (Almualla et al. 2021). Figure 5 shows the distributions of mass density ρ, electron fraction Y e , temperature T and Planck mean opacities κ planck at 1 d after the merger for run A (LS220 EoS) and run B (DD2 EoS).
For each model, we run Monte Carlo radiative transfer simulations using POSSIS and extract spectra for 11 viewing angles equally spaced in cosine from pole (face-on) to equator (edge-on) and for 100 epochs going from 0.1 to 30 days after the merger (with logarithmic binning of ∆ log t = 0.0576). In particular, we use a new version of POSSIS with state-of-theart heating rate libraries (Rosswog & Korobkin 2022), densitydependent thermalization coefficients (Wollaeger et al. 2018) and wavelength-and time-dependent opacities from Tanaka et al. (2020). Assuming Local Thermodynamic Equilibrium (LTE), the temperature is initialized according to the local heating rates at the start of the simulation and then updated at each time-step via Monte Carlo estimators of the mean intensity in each cell (Lucy 2003). In all the simulations, 10 6 Monte Carlo photons are created and followed in their propagation throughout the grid as they interact with matter (electron scattering and bound-bound transitions are considered) and eventually leave the ejecta.

KN ejecta and EoS
In order to study further dependencies of the KN light curves on the actual matter properties and EoS, we now focus on the ejecta. Susbsequently from NR simulations of BNS mergers, the asymmetry in the mass ratio, q, total mass M 1 + M 2 and the chirp mass M chirp have a profound impact on the amount, composition and velocity of the ejecta and this determines the luminosity of the KN. It is believed that due to the high compactness C ∼ GM/Rc 2 of individual NSs and the deep gravitational potential well, mat- Logκ planck (cm 2 g −1 ) Fig. 5. From left to right, distributions of mass density ρ, electron fraction Y e , temperature T and Planck mean opacities κ planck in the v y − v z velocity plane. The top row refers to run A (LS220 EoS) and the bottom row to run B (DD2 EoS), both with q = 1 and M chirp = 1.188 M ⊙ , see Table 1.
Maps are computed at an epoch of 1 d after the merger. The pixelization seen in the T (and hence also κ planck ) map is due to the temperature being computed with estimators in the code and thus being subject to Monte Carlo noise (see text for more details). Analytic functions are instead used for ρ and Y e and the distributions are therefore smoother. ter ejected in the merger is essentially stripped off external layers (Bauswein et al. 2013). Going further into the NS description, its layered structure can be mainly explained from two regions: an external crust and a central core. In the usual picture, matter is well beyond nuclear saturation density in the homogeneous nucleon core, ρ 0 ∼ 2.4×10 14 g cm −3 , while densities in the crust are typically less than ∼ 0.1ρ 0 . The crust is the most external, least gravitationally bound layer and is mostly constituted of so-called pasta phases (Horowitz et al. 2004) in the inner crust and, more externally, of ionic lattices at lower densities in a degenerate lepton sea (Chamel & Haensel 2008). The nature of the many-body approach to describe such relativistic quantum dense matter in NSs allows the large spread of EoS available in the literature, for reference see (Oertel et al. 2017). As mentioned, most of the NS mass is in the core in the form of baryons, although if densities are in excess of ∼ 6 × 10 14 gcm −3 a quark star or a hybrid NS with a quark deconfined inner core is hypothesized (Witten 1984;Ivanytskyi et al. 2019). The gravitational mass, M, and baryonic mass, M * , are quantities that determine the binding energy of the NS and some universal relations between them are available in the literature (Gao et al. 2020), based on the radius of a 1.4 M ⊙ NS, R 1.4 , M * = M + R −1 1.4 × M 2 showing that for different EoS this parameter is genuinely different (Dietrich & Ujevic 2017).
From GW170817 some constraints appear for the masses, i.e. total mass M ≈ 2.74 M ⊙ , and individual masses of M 1 ≈ (1.36 − 1.6) M ⊙ and M 2 ≈ (1.17 − 1.36) M ⊙ . From this analysis marginalized over the selection methods, an upper bound was obtained on the effective tidal deformability assuming a low-spin prior, which is consistent with the observed NS population. Let us remind that the binary tidal polarizability is defined as (1) The individual NS tidal deformabilities (i = 1, 2) are Λ i = 2 3 k 2,i C −5 i with k 2,i the Love number (Flanagan & Hinderer 2008;Lattimer 2019). Phenomenologically it has been found that for the EoS we consider Λ i = AC −6 i with A DD2 = 0.0112 ± 0.0002, A SFHo = 0.0099 ± 0.0002, A LS220 = 0.0111 ± 0.0002.
From GW170817 upper bounds were found Λ ≤ 800 at the 90% confidence level , which disfavors EoS that predict the largest radii stars. Tighter constraints were found in a follow up reanalysis (Abbott et al. 2019) with Λ = 300 +420 −230 (using the 90% highest posterior density interval), under minimal assumptions about the nature of the compact objects. The fact that NSs are not point-like objects imply there are strong tides accelerating the inspiral (with initial frequency f GW and producing a phase shift (Flanagan & Hinderer 2008) in the form δΦ t ≃ − 117 256 (1+q) 4 q 2 π f GW GM chirp c 3
To date there are dozens of EoS provided on phenomenological grounds typically using effective field theories invoking nuclear degrees of freedom (Oertel et al. 2017). In this work we have selected a reduced set of three different EoS: DD2 (Typel et al. 2010), LS220 (Lattimer & Swesty 1991) and SFHo (Steiner et al. 2013) that are compatible with maximum mass, being larger than ∼ 2 M ⊙ , and the binary tidal polarizabilities deduced in GW170817 ).
As a typical example of stellar configurations populating the binaries, Fig. 6 shows on the left panel the TOV solution to the mass-radius relation for the three EoS we use DD2, SFHo and LS220, along with excluded regions arising from GW170817 due to Bauswein et al. 2017 (left) and Annala et al. 2018 (right), and confidence regions from pulsar measurements from NICER (Miller et al. 2019b). Since the BNS cases we analyze correspond to specific values of individual NS parameters computed with our set of EoS we show in the right panel the values of chirp mass as a function of mass ratio q > 1 and the largest of the NS masses, M 1 .
We now start by discussing the different NS contributions to the KN light curve in order to analyze the ejecta parameters characterizing the BNS merger event and the subsequent EM emission. As previously obtained, one key aspect is the determination of the peak values in the bolometric light curve as they are related to the ejecta properties. In order to do this we must note that it is necessary to perform a smoothing, i.e. rebinning, procedure in order to manage the Monte Carlo noise from the set of simulations we use in our work.
Previous studies assume spherical symmetry although asymmetric ejecta have also been considered (Villar et al. 2017). State-of-the-art calculations  show that the opacities span the range κ ∈ [0.1 − 100] cm 2 g −1 depending on properties of the ejecta as ρ, T and Y e . As shown in Fig. 5, models show a toroidal component for Y e appearing clearly where neutron-rich elements concentrate versus a neutron-poor component for low θ values. The existence of a root mean square angle θ rms is often quoted, see for example Radice et al. (2018), and describes that ∼ 75% of the ejecta is confined within 2θ rms . Tem-peratures are approximately isotropic with additional enhancement towards equatorial orientations. These dependencies reflect on the Planck mean opacities and modulate the EM signal we obtain.
Note that the amount of matter ejected in the BNS merger will be only a tiny fraction of the total mass M 1 + M 2 as shown by detailed NR simulations. Dynamical mass and post-merger disk wind components play a fundamental role in determining the peak magnitudes. In this spirit and considering the fact that the ratio of dynamical to disk mass is less than ∼ 10% for the EoS under discussion we define M dd = M dyn + M disk . This quantity will be relevant for our EoS analysis as it is the matter stripped from the most external layers of the individual NSs during the merger. However, not all this mass will be dinamically ejected from the gravitational potential at the merger site since part of it remains for longer time scales than the measurable change in luminosity. These stripped layers are sensitive to the underlying EoS. We have verified that, for q = 1, this mass component can be fit to a Woods-Saxon distribution of the type M dd = a(1 + exp(b(M chirp − c)) −1 . More challenging, due to our reduced subset of runs with q 1, is the description of those BNS mergers with different individual NS masses. As a tentative fit we use a multiplicative factor for q ≳ 1 under the functional form M dd (q)/M dd (1) ∼ 1 + g(q − 1) h . We give the approximate fitting parameters for both dependencies in Table 2 Table 2.  Figure 7 shows the extrapolated values of the logarithm of peak brightness as a function of mass ratio q and chirp mass M chirp for EoS (from left to right) DD2, LS220 and SFHo. We select a polar view orientation. We can see how the imprint of the particular EoS along with the BNS realizations will affect the behavior of the peak.

KN peak magnitudes and EoS
According to previous semi-analytical estimates extracted from numerical simulations (Li & Paczyński 1998 43.1 ± 0.4 1.0 ± 0.1 1.9 ± 0.5 -0.6 ± 0.2 SFHo 41.0 ± 0.4 0.7 ± 0.1 0.9 ± 0.6 -1.8 ± 0.4 LS220 42.0 ± 1.5 0.3 ± 0.2 1.3 ± 1.9 -0.44 ± 0.43 Table 3. Fit of peak luminosity L peak as L peak = L 1d (M dd /M ⊙ ) α (v e j /c) β (Y e ) γ (κ/10 cm 2 g) δ for selected runs as a function of ejecta parameters M ej , v ej , Y e . For peak conditions we fix κ = 10 cm 2 g and δ = −0.5. Equatorial view is selected. See text for details. merger show dependencies on ejecta properties M ej , v ej , κ Planck ≡ κ, namely the ejected mass, ejected velocity, Planck mean opacity. In our work, however, we phenomenologically choose to explore the dependence of L peak on M dd , v ej , κ, Y e , due to the usefulness of M dd in relation to the EoS trends. We extract trends for L peak as L peak = L 0 (M dd /M ⊙ ) α (v ej /c) β (Y e ) γ (κ/10 cm 2 g) δ for edge-on orientation (cos θ = 0). Given the high-opacity conditions near the peak we fix κ = 10 cm 2 g and δ = −0.5 (see Table  3). Figure 8 shows the M dd mass for the three EoS under inspection as a function ofΛ. The latter is available as a GW measured quantity from a BNS event and from GW170817 there is an upper bound set atΛ ≲ 800 ) that we show as a vertical line. We notice that some of the DD2, LS220 instances we have considered do not fulfill this constraint. Let us highlight at this point that only for a subset of these, the runs appearing in Table 1, we have obtained KN spectra and light curves. Interestingly, the exclusion region set by theΛ ≲ 800 upper bound on M dd has relevant consequences on the peak log luminosity and time as shown in Fig 9. We select equatorial orientation. Although the set of runs we consider in Table 1 is reduced it seems to indicate the presence of limiting values for peak magnitudes. If such a case is confirmed by larger statistical samples it would be a genuine limit for the intrinsic magnitudes of KNe. From our plot we derive log L peak [erg/s] ≲ 41 and log t peak [day] ≲ 0.9. Note that the previous data points and bounds refer to an equatorial inclination, when considering a polar observer we obtain an enhanced luminosity yielding a bounding value log L peak [erg/s] ≲ 41.3. Additional investigations are  Table 1. We also indicate the upper limitΛ = 800  needed regarding the robustness of our finding and a careful reexamination of the caveats we have used in the simulations.  Table 1 also having q ≃ 1.We also indicate a conservative lower limit estimate of equatorial peak luminosity for a GW170817-like transient.
In order to better constrain the allowed regions on the Mass-Radius relation and thus the EoS (see Fig. 6), we select runs from Nedora et al. (2021) and Radice et al. (2018) with q ≃ 1. Top panel in Fig. 10 shows M dd as a function of the radius for DD2, LS220 and SFHo EoS. We can clearly see that for radii smaller than a given value (EoS dependent) the stripped mass dramatically plunges to vanishing values. For DD2 the NR cases are in the nearly vertical part of the M/R curve, thus this trend is not seen clearly, although radii smaller than ∼ 13.2 seem to exclude observable KNe. Thus we do not expect having KNe being observed for NSs with compactness beyond a limiting value as dictated by these radii. Our results agree with those for AT 2017gfo appearing in Breschi et al. (2021b) obtained from Bayesian analysis and semi-analytical models except for the SFHo EoS, since in our analysis for that specific EoS a BNS with M = 1.4 M ⊙ and q ∼ 1 luminosity peaks below experimental values. Additional correlations to central densities in TOV solutions have been explored (Breschi et al. 2021a). It is worth noting at this point that it is likely that there will be an additional EoS dependence associated to the BNS redshift that must be carefully analyzed, see (Ghosh et al. 2022), as it can impact the measurement of H 0 .
In the lower panel of Fig. 10 and for the runs from Table  1 with q ≃ 1 where KN light curves are simulated we show log L peak [erg/s] for an equatorial orientation as a function of the NS mass in the BNS, M 1 . In the same figure we indicate the deduced equatorial lower limit for the peak luminosity for AT 2017gfo (in dashed line). One can observe that for LS220 and DD2 EoS, the luminosity threshold for AT 2017gfo is reached. The only effect that drives larger peak luminosities between the two sets of values is the radius of the individual NSs, i.e. the compactness. Let us remind here that C DD2 < C LS220 for the stiffer DD2 EoS versus the softer LS220 EoS. Even though the effect seen here is only qualitative, due to poor statistics induced by our reduced set of simulated models, being able to have a library of NR simulated runs with different primary NS masses will allow us to isolate the effect of the EoS in the KNe luminosities. For SFHo EoS the high NS compactness predicted reduces the amount of ejected material and so the luminosity of the transient indicating that it seems not well suited to describe this KN. Additional analysis with improved statistics will be able to provide stronger constraints on the EoS.

KN light curves
After the BNS merger, the KN emission follows due to a combination of signals characterized by how nucleosynthesis proceeds. A red component is associated with a lanthanide-rich outflow, and a blue component with the lanthanide-poor counterpart. The high opacity of a lanthanide-rich (heavy r-process) ejecta delays and dims the light curve peak, and its large density of lines at optical wavelengths pushes the emission to redder wavelengths. Light r-process compositions will instead have a lower opacity and the emission associated with these outflows will have a faster rise and brighter luminosity peak with fluxes concentrated at shorter/bluer wavelengths. Alternative scenarios to explain the blue to red spectral evolution, such as dust formed by the KN, which extinguishes the light from the blue component and subsequently re-emits it at NIR wavelengths, have also been suggested (Takami et al. 2014). However, for AT 2017gfo this could be ruled out (Gall et al. 2017).
As the envelope becomes more transparent due to the expansion while the radioactive power decreases, the luminosity initially rises to a maximum (peak) before declining. A variety of KN models are available ranging from simple (semi-analytical), spherical, radioactive heating parametrized by a power-law and black-body emission (Li & Paczyński 1998;Kawaguchi et al. 2016;Metzger et al. 2015;Metzger 2019) to more refined ones that are nonspherical, with detailed modelling of radioactive heating and radiative transport (e.g. Kasen et al. 2017;Bulla 2019;Kawaguchi et al. 2020;Korobkin et al. 2021).
The radioactive decay of freshly synthesized r−process nuclei in the expanding ejecta at relativistic velocities powers a Fig. 11. Log L[erg/s] as a function of log t[d] comparing simulated models versus fits for runs C and H using SFHo EoS, see Table 1, for polar and equatorial orientations.The inflexion point signals the change in the decay law as appears in Table 4. quasi-thermal emission responsible for the KN. The peak flux F λ,peak and peak time t λ,peak in a given spectral band depend both on intrinsic (ejected mass, expansion velocity, composition) and extrinsic (viewing angle, distance, reddening) parameters. The 23 runs considered in this work (see Section 3 and Table 1) differ in terms of the ejecta properties like mass, M ej , average velocity, v ej , and average electron fraction, Y e . The latter, in particular, determines the photon opacity for a given wavelength of the broadband radiation, which can substantially differ from the opacity of iron group elements in the presence of a significant fraction of moderate (A = 90 − 140) or heavy elements . It is still unclear if NS mergers are the main source of r−process elements (Siegel et al. 2019), although high-opacity lanthanides have been inferred in GW170817 (see Ji et al. 2019 and references therein).
Opacities from Tanaka et al. (2020) used in this work are underestimated at t ≲ 0.5 d, due to the lack of atomic data at high ionization states, typical for these early epochs. Therefore, we focus our discussion on times t ≳ 0.5 d. In the same way the decline of the light curves shows two regimes, before and after the peak that are numerically determined by analyzing where the rapid decline starts. Near the peak, the structure of the luminosity is not trivial and displays some variability due to underlying dependencies on q/M chirp and EoS physics. In Fig. 12 we show the behaviour near the luminosity peak by plotting the spread in the logarithm of the luminosity as obtained frpeom the simulations with respect to a power law L ∼ t γ law, ∆L in units of erg/s, as a function of time, in days, for different orientations. The closer to the power law L ∼ t γ before the peak, the darker bluish the color, while the larger deviation is coded in reddish colors. The peak is denoted by a black dot. Each particular run we consider shows a distinctive feature depending not only on the EoS but also on the ejecta parameters. Runs with the same EoS show a big spread near the peak as q and M wind grows. Throughout the manuscript we label this point in the KN light curve of our analysis as peak brightness.
From semi-analytic estimates in Metzger (2017), the timedependent KN bolometric luminosity, assuming homogeneous Article number, page 11 of 21 A&A proofs: manuscript no. output and homologous ejecta has a functional form where α = 1.3. The peak luminosity L peak and time t peak are determined  roughly to be equivalent to the instantaneous rate at which radioactive decay is heating the ejecta. The EoS dependence will appear in these luminosity curves built in an indirect way in the ejecta parameters and composition.
As mentioned, the material ejected during and after the merger of two NSs (or a NS and a black hole) is strongly spatially asymmetric and characterized by multiple ejecta components with different geometries and properties, see Nakar (2020) for a recent review. As a result, the KN signal is found to be strongly dependent on the observer viewing angle (Bulla 2019;Darbha & Kasen 2020;Kawaguchi et al. 2020;Korobkin et al. 2021), with orientations in the merger plane (equatorial or edgeon inclinations) typically fainter and redder than orientations on the jet axis (polar or face-on inclinations). As a result, we end up with a viewing-angle dependence of the signal that offers a way to constrain the inclination angle from the KN. In order to qualify the impact of the inclination in the measured brightness decline, we have considered as examples runs C and H with SFHo EoS from Nedora et al. (2021), see Table 1. Both runs have a chirp mass M = 1.188M ⊙ , but different mass ratios q = 1, 1.43 respectively. Figure 11 shows log L in units of erg/s as a function of log t in days for those mentioned simulations along with fitting functions. Solid (dashed) lines depict polar (equatorial) orientation. We can see that realistic models such as these are well represented by log-log laws.  Table 4. Fit parameters for log L = β + αlog t (L in erg/s and t in days) in the bolometric luminosity before/after the peak (bp/ap) for polar and equatorial orientations. We use runs C and H with SFHo EoS, see Table  1.
As we can observe, there is a clear difference in the polar to equatorial orientation between the fit parameters β, and α, according to the law log L = β + α logt for the bolometric luminosity. We have verified that this also applies to the peak values, and as a general trend this behaviour is seen for all the runs we consider. Figure 13 illustrates this. In the top panel we show the logarithm of the luminosity of the peak as a function of logarithm of peak time for the models we consider in this work. As usual L is expressed in erg/s and t peak in days. In the bottom panel we show the 7-day variation in the logarithm of the bolometric luminosity of peak as a function of logarithm of peak time. LS220, SFHo and DD2 EoS are shown in red, orange and blue colors, respectively. For a given run, full and empty symbols represent different cos θ orientations correspond to models listed in Table  1 from Radice et al. (2018) and Nedora et al. (2021), respectively. The plot in log-log scale shows clear correlations among peak luminosity and time of the peak: the brighter the KN is, the longer it takes to peak and the smoother the change in magnitude is. For each run the dot series depicts the angular modulation of the peak value.
In addition, Fig. 14 shows the log peak luminosity (in erg/s) as a function of seven day gradient in log bolometric luminosity (in erg/s). We can see that the fainter they are the larger decay in bolometric luminosity after a week they experience. As seen from Figs. 13 and 14, there is a clear modulation from observing angle θ. In order to quantitatively determine its impact we have considered the models in Table 1 corresponding to a common EoS and fit the angular dependence as obtained from the 11 orientations for each simulated run. Table 5 shows the results of a fit to the peak value log L peak (x) = c 0 + c 1 arctan(c 2 (x − c 3 )) with x = cos θ. At high inclinations (low x values), arctan(x) ≃ x − x 3 3 . We find that, generally speaking, a function of this type describes satisfactorily the angular modulation of luminosities.  Table 5. Fit of the expected peak bolometric luminosity (in erg/s) as a function of x = cos θ, log L peak (x) = c 0 + c 1 arctan(c 2 (x − c 3 )). We consider runs from Table 1 with selected values of q and M chirp . This is illustrated in Fig. 15, that shows the logarithm of peak of the luminosity (in erg/s) as a function of cos θ for selected values of mass ratio q = 1 and M chirp . Runs A, B and C have a fixed M chirp = 1.118 M ⊙ , runs 1, 2 have M chirp = 1.22 M ⊙ while run 3 has M chirp = 1.39 M ⊙ . Note that (q, M chirp ) will be provided by analysis following a positive detection from GW detectors. In order to obtain coherent underlying trends we fix the hadronic EoS. LS220, SFHo and DD2 EoS are shown in red, orange and blue colors, respectively. The actual values of the fitting parameters are provided in Table 5.

KN spectro-photometry and H 0 determination
In this section, we study the observational prospects, and Hubble constant determination, by using the new MAAT IFU on the GTC from the potential KNe emission simulated for a sample of NS mergers described above, and listed in Table 1. Figure 16 shows the absolute magnitude per post-merger epoch in r band for some of the KNe types that we are considering in our study. This plot shows the drop in observed magnitude when the KN is oriented towards the merger plane of the neutron stars. To study the detectability and get an idea of what type of KNe we will be able to observe with MAAT, we explore the dependencies of the S/N with distance, inclination angle and intrinsic physics of the KNe.

Simulating KN spectro-photometry observations
We have built a pipeline to generate realistic KN spectrophotometric data as observed with MAAT, and perform an endto-end analysis of the spectro-photometric light curves. The input of the pipeline consists of templates to model the theoretical flux of the source of interest and provide the relevant information for the MAAT exposure-time-calculator (ETC). The ETC provides the output S/N (exposure time) as a function of wavelength for a requested exposure time (S/N), taking into account the transmission through the atmosphere and the telescope and  Table 1 and displays the departure of the logarithm of the luminosity from a power law L ∼ t γ near its peak as a function of the orientation, cos θ, and time (in days) elapsed since merger. Black dots mark the location of the peak derived from our analysis. See text for details. instrument (OSIRIS+MAAT) optics. It works for the whole suite of OSIRIS gratings and for an optional range of source geometries. The complete code package and instructions are available on the MAAT website 5 . The ETC output includes the observed simulated spectra in physical units (erg/cm 2 /s/Å). We convolve the spectrum with the standard gri SDSS filters transmission curves to extract also synthetic photometry.
We plan to use the MAAT IFU with the OSIRIS R1000B and R1000R grisms to follow-up KNe alerts from the upcoming GW detector runs or after a potential identification by imaging surveys, as described in Section 2. We expect to be very sensitive to observe faint KNe, over the 3600 − 10000 Å spectral range, down to ∼22 mag with 1800s exposures per grism, see also Fig. 17. Our aim consists of showing the potential of MAAT observations to set meaningful constraints on two of the most desirable research fields related to compact binary mergers and their electromagnetic counterpart: the study of the EoS (see Sec. 3), and the potential contribution to cosmology as by improving the Hubble constant estimation from standard sirens in the local universe (Holz & Hughes 2005). 5 http://maat.iaa.es/maat_etc.html For population studies and long-term predictions, it is required to collect a significant statistical sample of KN data. However, only one KN has been successfully identified to date from a BNS merger (GW170817-AT 2017gfo).
We aim to use the standard sirens method from GW data together with prior from the KN data as an independent constraint on the inclination of the system and improve the derived H 0 uncertainties. We perform the analysis with synthetic MAAT data for the only two BNS mergers detected to date with GW detectors: GW170817 and GW190425.
In this section, we present synthetic data of follow up KNe observations with MAAT. Therefore, we consider an exposure time of 1800 seconds (0.5 hr) with grims R1000B and R1000R. Dark moon, 1 airmass and a seeing of 0.8" is assumed for all the simulated data in this analysis.

GW170817 KN light curves
To simulate GW170817 KN, we use the X-shooter spectra that were taken for AT 2017gfo (Smartt et al. 2017;Pian et al. 2017) as our input templates and then compute how it would have been observed with MAAT using our pipeline. The spectra and the  Table 1 from Radice et al. (2018) and Nedora et al. (2021), respectively. For each run, orientations are depicted by the dot series. LS220, SFHo and DD2 EoS are shown in red, orange and blue colors, respectively. photometric data are publicly available 6 . The photometric data set for AT 2017gfo cover filters from the UV to the IR up to ∼10 days after the merger. In this analysis, we take the observed photometric data of AT 2017gfo in gri bands on epochs 1.427, 2.417, 3.412, 4.402 days for a fair comparison with MAAT on wavelength range and timespan for our analysis. Fig. 18 shows the resulting time series spectra simulated for MAAT together with the overlap gri observed photometry. GW170817 happened at relatively close distance and the S/N per wavelength with MAAT would have been 140 at 7000Å.

GW190425 KN light curves
No KN was identified for GW190425. From the GW data, the luminosity distance was inferred to 159 +69 −72 Mpc, and the inclination angle was poorly constrained with a median value of ∼ 40 • (Abbott et al. 2020a). For our analysis, we simulate synthetic KN data assuming D L = 159 Mpc and θ = 40 • . We consider two of the models presented in this study: run3 and run10 (see Table  1) which chirp mass (M chirp ) and mass ratio (q) are consistent 6 https://sid.erda.dk/wsgi-bin/ls.py?share_id=bJ0ZS1CdFj Fig. 14. Luminosity in erg/s as a function of seven day gradient in log bolometric luminosity. Full and empty symbols correspond to models listed in Table 1 from Radice et al. (2018) and Nedora et al. (2021), respectively. For each run, the orientations are depicted by the dot series. LS220, SFHo and DD2 EoS are shown in red, orange and blue colors, respectively. with GW190425. Run3 assumes DD2 EoS and run10 LS220, and the resulting KN varies in ejecta parameters such as the mass of the ejecta components, the electron fraction (Y e ), and the expansion velocity of the ejecta (v e j ). The spectra of GW190425 have been simulated using the run3 and run10 described previously (see Table 1). These two runs are tailored to a BNS merger with a chirp mass (M chirp ) and mass ratio (q) consistent with GW190425. Run3 produces a KN with M wind = 5.88 × 10 −3 M ⊙ and M dyn = 1.2 × 10 −3 M ⊙ while run10 generates a KN with M wind = 2.10 × 10 −4 M ⊙ and M dyn = 3.0 × 10 −4 M ⊙ . For the same BNS, q, M chirp parameters with M 1 = 1.6M ⊙ the more compact i.e. smaller radius realization, see Figure 6, run 10 LS220 EoS produces fainter light curves and then larger apparent magnitudes. EoS differences thus translate into significant changes in the luminosity and evolution of the light curves. For the KNe simulated for GW190425 we find that DD2 produces brighter and longer-lived KN in the wavelength range of MAAT. Figure 19 shows the resulting time series spectra for run3, while Figure 20 the same for run10. For both cases, the spectral series are much noisier than the one for AT 2017gfo. For run3, the SNR estimated at 7000Å goes from ∼ 8 lowering to ∼ 2 and ∼ 1 within 2 days, but at later epochs we report an increase of the SNR at red wavelengths, as it is expected from the KN evolution.
Run10's spectral series shows even fainter fluxes (the SNR peaks at around 1 and drastically decreases) with magnitudes shallower than 22-23 mag, which would be hard to detect with MAAT, and even with current optical surveys. If MAAT followed up on such KN, we would be able to gather useful information only by rebinning the spectrum (Fig 20), which improved the SNR by a factor of 5. The broadband photometry, in this case, would need to be sensitive to 23-24 apparent magnitudes. The incoming Rubin Observatory Legacy Survey of Space and Time (LSST, Ivezić et al. 2019) will be able to detect future KNe at these fluxes.
As there is no EM counterpart associated with this event, therefore, we choose to simulate synthetic gri photometry with standard errors for optical depth ∼ 22 mag. We aim to compare MAAT observations with what would have likely been observed by optical surveys if our assumed KN were detected. Run3 and run10 KNe at 159 Mpc would be hard to observe with optical surveys reaching ∼ 22 mag limiting magnitudes as the sky background would largely dominate (S/N> 5). We plan to explore synthetic photometry that resembles more dedicated follow-up observations in the future, therefore, reaching greater sensitivity at shallower magnitudes.

Measurement of the Hubble constant
As noted above, there is a 5σ tension between the inferred value of H 0 from the cosmic microwave background (Planck Collaboration et al. 2020) and the local Cepheid distance ladder (Riess et al. 2022). This tension can be a sign of new physics beyond the standard cosmological model or unknown systematics. BNS mergers with an identified EM counterpart (e.g.  are an excellent probe to understand the origin of this tension since they measure H 0 completely independently of the distance ladder, or to any assumptions on the cosmological model. The use of GW as "standard sirens" to estimate H 0 represents a powerful distance indicator in cosmology (Holz & Hughes 2005). The luminosity distance is calculated directly from the GW signal, despite it is degenerate with the inclination angle (θ). However, it has been proposed to use the electromagnetic data to constrain the inclination angle and break the degeneracy (Nissanke et al. 2013).
The standard siren method can also be used in the case the exact identification of the galaxy hosting a BNS event is missing. This approach, originally proposed by Schutz (1986), takes into account all the bright host galaxies located in the GW localization region as potential host galaxy of the BNS merger event, and then infer the H 0 value as the average value obtained from the ensemble analysis. A similar approach has been applied for the case of GW170817, where a value of H 0 = 77 +37 −18 km/s/Mpc has been inferred (Fishbach et al. 2019). However, in the case of detecting an EM counterpart, it is possible to unambiguously identify a host galaxy and infer a redshift to the source. Moreover, the EM counterpart also contains essential information regarding the inclination of the binary system. This possibility has been explored with the KN and the GRB associated with GW170817. Here, we infer those constraints by simulating how such a KN would be observed with MAAT IFU spectroscopy on the GTC telescope. Along with simulations for AT 2017gfo, we also project the improvements for the well-observed GW event GW190425. It will depend on factors like the S/N of the observed spectrum and the shape of possible features related to the geometry and composition of the KN ejecta. We join GW and KN data as presented in Dhawan et al. (2020), by taking the joint posterior probability distribution in H 0 and cos(θ) from the GW data and applying a prior on the inclination angle from the KN.
The methodology followed in our analysis is explained in what follows. We use the synthetic data from MAAT and the photometry for each BNS (GW170817 and GW190425) and perform parameter inference using simulations with the radiative transfer algorithm POSSIS (Bulla 2019). The code predicts 3D models for a set of ejecta parameters and geometries including the inclination relative to the observer. The predicted viewing angle dependence allows us to infer the cos(θ) parameter from the KN data independently of the GW. A set of models from POS-SIS is considered, including the grid presented in Dietrich et al. (2020) and the 23 models presented in this study (see Section 3 and Table 1). To infer the KN parameters, we perform χ 2 fits and extract the χ 2 probability distribution of the viewing angle parameter marginalizing over the ejecta parameters, which include the masses of the ejecta components and the half opening angle of the lanthanide-rich region. The obtained viewing angle distribution is used as a prior for cos(θ) to reweight the GW joint posterior from the GW analysis. The extra and independent constraint from the KN data on cos(θ) breaks the degeneracy with   (Smartt et al. 2017;Pian et al. 2017). Note that AT 2017gfo was found in a nearby galaxy at 40 Mpc and with an almost polar inclination angle direction ∼20 • (Hotokezaka et al. 2019). the luminosity distance and results in an improvement in the H 0 uncertainty with respect to GW data alone. This method is used in this section to analyze the contribution of viewing angle inference with KN observed with MAAT to the H 0 constraints.

H 0 estimate from GW170817
For the case of GW170817, we take the posterior distribution presented in . The H 0 at joint probability distribution for this event was calculated by combining the luminosity distance with the redshift measured for the host galaxy NGC 4993 of AT 2017gfo.  Table 1) associated with GW190425 as if observed with MAAT with the OSIRIS R1000B and R1000R grisms for 0.5 hr exposure per epoch. Note that the GW190425 event was located around 160 Mpc and with an inclination angle ∼ 40 • . The broad-band SDSS gri magnitudes are also shown, see text, with error corresponding to SNR=5, we notice that the observations have to be sensitive to 22 mag and deeper to get a detection for this kind of KN. The graph shows the MAAT spectra directly obtained from the ETC (shaded) and the rebinned spectra overplotted with solid line on the same color. Note that this is synthetic data generated for a KN at 159 Mpc and θ = 40 • .
The prior inclination angles are calculated following a χ 2 distribution that indicates the goodness of fitting our simulated AT2017gfo MAAT spectrum, shown in Fig. 18, with the BNS model grid from the radiative transfer code POSSIS presented in Dietrich et al. (2020). The resulting inclination angle distribution is shown in the red curve of panel b) in Fig. 21. The θ constraint is significant to the order of 10 percent, and the H 0 uncertainties improve by a factor of ∼40% with respect to GW Fig. 20. Same as Fig. 19, but adopting a KN model with similar merger parameters but different EoS, i.e. run10 -LS220 EoS. This produces fainter KN spectra, whom significantly impacts the SNR of the MAAT spectra, and hence the constraining power on H 0 . Nevertheless, thanks to the spectro-photometry capability and spectral resolution of MAAT, the SNR is enhanced by rebinning the spectra. The solid line shows the rebinned spectra in 50 Å and the marker the convolved photometry in gri assuming SNR = 5. data only (see Table 6). As a comparison, it is also shown in Fig. 21 a yellow band corresponding to the constraint obtained from the superluminal motion of the associated short GRB jet in Hotokezaka et al. (2019).
The same analysis has been performed for the photometric data of AT 2017gfo introduced in Section 5.1.1. Dhawan et al. (2020) obtained a 24% improvement with the full UV-optical-IR dataset. We analyze only data on the same epoch and wavelength range as the MAAT spectra (3600-9000 Å) with gri filters and 4 epochs: 1.427, 2.417, 3.412, and 4.402 days. The final H 0 constraint from GW plus the described photometry improves only by a factor of 5% which compares to the ∼ 40% from GW+MAAT spectra, probing the benefit of MAAT data in KN parameter inference for the wavelength range (see Table 6). Run 10 (LS220) GW+MAAT prior 1.0 GW190425 Run 10 (LS220) GW+MAAT (50Å) 0.84 Table 6. σ(H 0 ) ratio results. Note that we are not paying attention to the H 0 value in this analysis, but we focus mainly on the resulting uncertainty improvement to GW-only.

H 0 estimate from GW190425
For GW190425 we generate the H 0cos(θ) joint posterior from the luminosity distance and assuming a redshift z = 0.035 (Abbott et al. 2020a). Note that the central value of the H 0 distribution is not physically meaningful as it purely depends on the chosen redshift input, while the shape indicates the actual uncertainties. We lack redshift information for GW190425. We cannot correct for recession velocity with no EM counterpart. Therefore, we cannot measure the real H 0 value using this method.
Nevertheless, we assume a redshift of 0.035 for the synthetic GW190425 -KN and apply the same to the GW luminosity distance data from which we estimate H 0 ≈ (c · z)D L , with c the speed of light, z the redshift and D L the luminosity distance. This analysis aims to quantify the constraining power of MAAT, for what the relevant results are the uncertainties in the H 0 as the improved constraint with KN compared to GW only. We constrain the θ parameters that we present and described above in the paper. Regarding the simulated photometry of GW190425, we assume gri SDSS filters and a 5 sigma detection at 22 mag. The SNR for such photometry goes below 5 (no detection), which significantly lowers the constraining power, as we can see in the blue curves in Figure 22.
The H 0 analysis on the synthetic KN data highlights the potential of MAAT on improving the constraints even for more distant KNe. If we assume that the models considered are representative of KNe in the Universe, even one single spectrum with the resolution of the R1000 grisms can significantly constrain the θ parameter and improve H 0 by 28%, or even further when using "cleaner" data like the 50 Å binned spectra shown in orange in Fig. 22, providing a 40% improvement (see Table 6) with respect to GW data only. Our results on the H 0 for GW190425 depend on the assumed KN model and assumed redshift, therefore, we do not intend to provide a realistic H 0 value for GW190425, but only to show the improved uncertainties and the dependence on the model assumptions, including the chosen EoS.
We did the exercise with the simulated KN assuming EoS LS220 (run10), which shows spectral noise after the first epoch as the background flux dominates after day 1. Suppose observing 0.5 hr with gratings R1000 (Fig. 20, we are not able to constrain the KN parameters with MAAT spectra due to the poor S/N, nevertheless with a rebinning of 50Å, we obtain a factor of ∼5 better SNR and the inclination angle prior provides a 15% improvement in the H 0 uncertainties (see Table 6).
We note that for KNe at ∼150 Mpc, it is critical to obtain early data to constrain our model parameters. In the case of run3, most of the constraint comes from the first epoch assumed to be taken within the first day after the merger, indicating the need for fast response efficient communication of KN candidates. The constraining power of our data depends on the S/N, the consequence of the observing strategy, the intrinsic brightness of the KN, the luminosity distance, and the orientation, among other factors. The ability to constraint model parameters is also highly dependent on the model grid, as our models need to be able to reproduce the main features of the KN to provide sensitive results. This study assumes that our models represent the produced KN emission.
We aim to understand if the EoS induces systematics in the measurement of the Hubble constant. This study only compares run3 and run10 (described above) to simulate a KN associated with GW190425 with two different EoS: DD2 and LS220. We did the exercise to simulate MAAT data for both runs but requiring the same S/N of 5 at 6000Å, finding that the pipeline returns similar prior to the inclination angle in both cases, which indicates that the constraining power is mainly dependent on the SNR for these types of KN.
We plan a follow-up study to investigate the effect of the EoS in H0 precision for a population of KNe. In practice, it will not be possible to obtain the same S/N at all luminosity distances, as we need to consider the exposure time needed. We compare the inclination angle priors and H 0 estimation for run3 and run10 with the same exposure time of 1 hr for the KNe located at 100 Mpc MAAT prior gri prior H19 Fig. 21. Hubble constant determination with prior distribution from GW170817. a) shows the 1 sigma and 2 sigma contours on the parameter space on inclination angle and H 0 with no KN constraint (black) and with KN constraint from MAAT data (red) and broadband gri photometry (blue). b) shows the EM priors applied to the GW posterior extracted from MAAT data (red) and gri AT 2017gfo observed photometry (blue). c) The marginalization on H 0 for the GW posterior (black) and the reweighted posterior with MAAT prior (red) and gri prior (blue) applied. This event occurred at a luminosity distance of 40 Mpc, corresponding to a redshift of around 0.009.
(still consistent with GW190425 within 1 sigma). Run3 SNR is 4 times larger than run10, which is around 5 at 6000Å.

Discussion and conclusions
KN multimessenger emission produced in the BNS merger can shed light on fundamental unresolved questions, such as the behaviour of matter at high densities inside NSs or the precise determination of the Hubble constant, H 0 . In this work, we show that obtaining spectro-photometric light-curves of future KNe will allow us to constrain H 0 more precisely than has been done previously from GW-only or using broad-band light curves. This paves the way for a yet more precise complementary measurement of the Hubble constant. We also show how KN peak magnitudes are correlated to NS properties that closely depend on the EoS. In particular, stellar radii and compactness (M/R ratio) are key to understand the amount of matter ejected, and thus the peak brightness. In our study we find an indication that there are limiting values of compactness to obtain measurable light curves.
Using the 3D radiative transfer code POSSIS (Bulla 2019), we predict the KN signal for a set of input models from Nedora et al. (2021) and Radice et al. (2018). We analyze the dependence of different physical parameters governing the KN emission, i.e. ejecta properties specific to the particular BNS configuration based on mass ratio of NS masses, q, and chirp mass M chirp as well as inclination angle and cosmological distance (provided by GW measurements).
We detail how MAAT on the GTC, a new optical integral field unit that employs the image slicing technique to perform absolute flux calibrated spectroscopy in a field of view of 8.5 ′′ × 12.0 ′′ , is well suited to shed light on the characteristics of the optical emission of KNe via follow-up of BNS events.
We find that the signal modulation as a result of the inclination angle must be treated in a systematic fashion as it can result in up to one order of magnitude difference at peak luminosity. Fainter events experience a more rapid evolution of the bolometric luminosity seven days after the merger. The brighter the KN the later it peaks, with a consistent trend of peak times within ∼ 10 days. It is for this reason that setting an early observation strategy is key for capturing the essential physics governing the transient event. As we have explained, the underlying EoS describing the interior of individual NSs in the binary with given (q, M chirp ) determines the amount and properties of ejecta. We consider a reduced but representative set of EoS, namely DD2, SFHo, and LS220. We find that the light curve for times t < t peak can yield sizable differences, up to one order of magnitude, thus showing the importance of a rapid follow-up for these transients.
We further study the dependence of the luminosity peak with average ejecta mass, velocity, and electron fraction as obtained from NR samples. We find a strong dependence between dynamical plus disk mass (M dd ) and chirp mass for q = 1 events; and weaker for q > 1 events. Even if the disk mass does not play an effective role in the KN EM emission (only indirectly as a source of the wind mass), it is the bulk of matter originating from the layers less gravitationally bound to the NSs and thus susceptible of being tidally ejected during the violent merger event. We find that the binary tidal polarizability (yet another parameter available from GW observables; presently constrained from GW170817 to beΛ ≲ 800) can effectively limit the amount of mass ejected (for a given EoS) and constrain the peak luminosity. MAAT prior MAAT (50Å) prior gri prior Fig. 22. Hubble constant determination with prior distribution from GW190425. This figure follows the same panel structure as Fig. 21. The orange contour adds the constraint from a rebinned spectrum of 50Å. The synthetic KN used to compute this plot assumes merger parameters consistent with GW190425 and EoS DD2 (run 3, see Table 1) a luminosity distance of 159 Mpc (z ∼ 0.035) and a inclination angle of ∼ 40 • . Note that gri prior is computed for a survey with depth ∼ 22 mag while the MAAT 50Å prior comes from a rebinning of the observed MAAT spectra, which improves SNR.
seems to indicate that there is yet another complementary constraint on the combined NS mass -radius diagram. From the set of EoS analyzed we find that some EoS are not capable of producing such bright KNe as measured for AT 2017gfo. The resulting KN from the NS coalescence is mostly dictated by how nuclear matter behaves and the degree of compactness attainable for a NS. We plan to perform a more quantitative study on this matter by simulating KN for a given GW event at a reasonable luminosity distance such that we can compare the analysis and understand the implication that the EoS has on the optimal observing strategy, and ultimately on the H 0 estimation.
We predict the improvement of the H 0 estimate with MAAT data for a ∼ 22 limiting magnitude survey. We conclude that fitting the light-curves over the whole spectral range (3600-10000 Å) provides a better constraint than that from using GW-only and broadband (gri) photometric data. The resulting improvement on H 0 would have been around 40% if MAAT would have observed AT 2017gfo. This compares to other analysis found in the literature like that of (Dhawan et al. 2020) who obtains a 25% improvement from the full UV-Optical-IR dataset during 10 days after the merger. GW190425, the possible KN at ∼159 Mpc, is expected to be much fainter than AT 2017gfo because of the larger distance and less favorable inclination angle of around 40 • (Abbott et al. 2020a).
The resulting constraint in H 0 largely depends on the EoS assumed, as it impacts mostly on the brightness of the KN emission. We investigated the effect of using two different EoS (DD2 and LS220) for a merger like GW190425 with M chirp = 1.39 M ⊙ , M 1 = 1.6 M ⊙ , and q = 1. At 159 Mpc, the run10 KN observed with MAAT is not constraining unless rebinning (increase SNR) is done, what leads to an improvement of ∼ 15% mainly from the first epoch obtained within the first day (see Table 1). Run3 KN instead shows a better improvement on the H 0 precision of ∼ 25% to ∼ 40% with rebinning (see Fig. 19 and Table 1).
We also investigated for a KN at 100 Mpc arising from the same merger event but assuming two different EoS (again DD2 run3 and LS220 run10) with the same S/N (S/N= 5 at 6000Å) or same exposure time (1 hr). The results show that for a consistent S/N, if the two models considered are representative, we are able to constrain H 0 similarly, while for the same exposure time the constraining power shows a clear dependence on EoS (see Fig. 23) due to the much fainter magnitude and poorer S/N of LS220 run10. We plan a detailed follow up study to understand possible systematics and bias on H 0 introduced from EoS uncertainties and other properties. The possibility to get early KN spectra with a high S/N with MAAT will also facilitate the identification and the study of spectral features embedded within the KN ejecta, paving the way for abundance studies with spectral synthesis codes, and allowing us to neatly determine the kinematics of KN components.