DRAFTVERSION OCTOBER 4, 2017 Typeset using LATEX twocolumn style in AASTeX61

THE DISTANCE TO NGC 4993 – THE HOST GALAXY OF THE GRAVITATIONAL WAVE EVENT GW 170817

JENS HJORTH,1 ANDREW J.LEVAN,2 NIAL R.TANVIR,3 JOE D.LYMAN,2 RADOSŁAW WOJTAK,1 SOPHIE L.SCHRØDER,1 ILYA MANDEL,4 CHRISTA GALL,1 AND SOFIE H.BRUUN1

1Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, DK-2100 Copenhagen Ø, Denmark 2Department of Physics, University of Warwick, Coventry, CV4 7AL, UK 3Department of Physics and Astronomy, University of Leicester, LE1 7RH, UK 4Birmingham Institute for Gravitational Wave Astronomy and School of Physics and Astronomy, University of Birmingham, Birmingham, B15 2TT, UK

ABSTRACT The historic detection of gravitational waves from a binary neutron star merger (GW 170817) and its electromagnetic coun- terpart led to the first accurate (sub-arcsecond) localization of a gravitational wave event. The transient was found to be ∼1000 from the nucleus of the S0 galaxy NGC 4993. We here report the luminosity distance to this galaxy using two independent meth- ods: (1) Based on our MUSE/VLT measurement of the heliocentric (z = 0.009783 ± 0.000023) we infer the systemic −1 recession velocity of the NGC 4993 group of galaxies in the CMB frame to be vCMB = 3231 ± 53 km s . Using constrained cos- −1 mological simulations we estimate the line-of-sight to be vpec = 307±230 km s , resulting in a cosmic velocity −1 of vcosmic = 2924 ± 236 km s (zcosmic = 0.00980 ± 0.00079) and a distance of Dz = 40.4 ± 3.4 Mpc assuming a local Hubble −1 −1 00 00 constant of H0 = 73.24 ± 1.74 km s Mpc . (2) Using HST measurements of the effective radius (15.5 ± 1.5 ) and contained intensity and MUSE/VLT measurement of the velocity dispersion we place NGC 4993 on the Fundamental Plane of E and S0 galaxies. Comparing to a frame of ten clusters containing 226 galaxies this yields a distance estimate of DFP = 44.0 ± 7.5 Mpc. The combined redshift and FP distance is DNGC4993 = 41.0±3.1 Mpc. This ‘electromagnetic’ distance estimate is consistent with +2.9 the independent measurement of the distance to GW 170817 as obtained from the gravitational wave signal (DGW = 43.8−6.9 Mpc) and confirms that GW 170817 occurred in NGC 4993.

Keywords: galaxies: distances and — galaxies: fundamental parameters — galaxies: individual (NGC 4993) 2 HJORTH ET AL.

1. INTRODUCTION The data used in this paper are described in Levan et al. GW 170817 was the first gravitational wave event arising (2017). NGC 4993 was observed with the Very Large Tele- from a binary neutron star (NS) merger to have been de- scope and the MUSE Integral Field Spectrograph on 18 Aug tected by LIGO/Virgo (LIGO & Virgo Collaboration 2017a). 2017. The source was localized to a sky region of 28 square de- We also used an Advanced Camera for Surveys F606W grees purely using gravitational wave data from the three in- (HST) image of NGC 4993 ob- terferometers. As such it provided the first realistic chance tained on 28 April 2017. The image was reduced via 00 of detecting an electromagnetic counterpart, as outlined in astrodrizzle, with the final scale set to 0.07 . LIGO & Virgo Collaboration(2017b): Two seconds after the merger, Fermi and Integral detected a weak gamma-ray burst, 3. REDSHIFT DISTANCE and half a day after the event, an optical (Coulter et al. 2017) and near-infrared (NIR) (Tanvir et al. 2017) counterpart was 3.1. Redshift of NGC 4993 localized to sub-arcsecond precision, ∼ 1000 from the nucleus An updated value of the heliocentric redshift of NGC 4993 of the S0 galaxy NGC 4993 (LIGO & Virgo Collaboration is obtained from the MUSE/VLT observation reported by 2017b). It exhibited an unprecedented optical/NIR lightcurve Levan et al.(2017). Based on fits to the absorption lines from and spectral evolution (LIGO & Virgo Collaboration 2017b; stars in the center of the galaxy, the heliocentric recession ve- −1 Tanvir et al. 2017), strongly suggestive of formation of very locity is determined to be vhelio = 2933±7 km s correspond- heavy elements (lanthanides) out of the tidally ejected neu- ing to a heliocentric redshift of zhelio = 0.009783 ± 0.000023. tron rich material in the merger (Tanvir et al. 2017; Pian et al. 2017). It is highly likely that the NS-NS merger occurred in NGC 3.2. Redshift of the NGC 4993 group of galaxies 4993. A precise distance to NGC 4993 is therefore required NGC 4993 appears to be a member of a group of galax- in order to understand the energetics of the event. Further- ies, see Figure1. Crook et al.(2007) list it to be one of 46 more, since NS-NS mergers provide a route to H0 via the di- members of a group of galaxies in the 2MASS redshift sur- rect “standard siren” measurement provided by the LIGO and vey (group number 955, with an average heliocentric veloc- Virgo Collaborations (LIGO & Virgo Collaboration 2017c), ity of 2558 km s−1 and a velocity dispersion of 486.5 km s−1). it is valuable to have an accurate redshift as well as indepen- However, NGC 4993 is located ≈4.2 Mpc in projection from dent electromagnetic distance estimates. The redshift, espe- the center of this putative group and a non-detection in X- cially if corrected to the Hubble flow, is vital in constructing rays appears to rule out such a rich relaxed group. Makarov the gravitational wave Hubble diagram, while an independent & Karachentsev(2011) find NGC 4993 to be one of 15 distance estimate can confirm the galaxy association. members with a velocity dispersion of 74 km s−1. Recently, To our knowledge, there are no direct measurements of the Kourkchi & Tully(2017) found it to be one of 22 members distance to NGC 4993 (Tully et al. 2009). NGC 4993 be- of a group with a mean heliocentric velocity of 2995 km s−1 longs to a group of galaxies (see Section 3.2). Two galax- and a velocity dispersion of 118 km s−1 (our calculation of ies have distances in the Cosmicflows-3 data release (Tully the rms velocity spread of the 22 galaxies is 155 km s−1). et al. 2016). NGC 4970 (PGC 45466) has a Fundamental However, it is evident from Figure2 that five of the galax- Plane (FP) distance of 39.8±10.3 Mpc from the 6dF Galaxy ies have a much higher recession velocity than the rest and Survey (Springob et al. 2014) while ESO508-G019 (PGC so do not appear to be part of a relaxed group. We therefore 45666) has a Tully–Fisher distance of 37.5 ± 5.6 Mpc from exclude these five galaxies and consider the properties of the I-band and Spitzer photometry (Willick et al. 1997; Springob remaining 17 galaxies. We also update the NGC 4993 veloc- et al. 2009; Sorce et al. 2014). Combined, these estimates ity from 2902 km s−1 to 2933 km s−1 (as measured here) and suggest a group distance of 38.0 ± 4.9 Mpc. compute a mean heliocentric velocity of 2921 km s−1 and a Section 2 summarizes the origin of the data used to pro- velocity dispersion of 53 km s−1. This small velocity disper- vide new distance estimates to NGC 4993. In Section 3 we sion and the large extent of the structure suggest a relaxation provide an updated distance along this avenue. We use a new time larger than 1010 years, and so the structure is unlikely redshift and estimate of the peculiar velocity of NGC 4993. to be a relaxed group. We therefore adopt the velocity dis- In Section 4 we obtain the first direct estimate of the dis- persion of the structure as the uncertainty in√ the recession tance to NGC 4993 itself using the Fundamental Plane (FP) velocity (rather than the error in the mean, 53/ 16 = 13 km method. We discuss the results in Section 5 and summarize −1 −1 s ), i.e., vNGC4993,group = 2921 ± 53 km s . This structure is in Section 6. shown in Figure1. Note that NGC 4993 is at the outskirts of the structure while its velocity is very close to the mean 2. DATA velocity of the structure. DISTANCE TO NGC 4993 3

Using the WMAP measurement of the cosmic microwave frame of the CMB with respect to an observer located in the background (CMB) dipole (Hinshaw et al. 2009) we obtain1 Local Group. For the sake of readability, we downsampled a recession velocity in the frame defined by the CMB of from the selected dark-matter particles. The peculiar veloci- −1 vCMB = 3231 ± 53 km s . ties result primarily from the proximity of the Great Attractor (49 degrees from NGC 4993). 3.3. Peculiar velocity Using a ±250 km s−1 velocity range centered at the To estimate the ‘cosmic velocity’, i.e., the recession veloc- measured CMB rest-frame velocity of NGC 4993, i.e., −1 ity corresponding to pure Hubble flow, we need to take into vCMB = 3231 km s , we find a mean peculiar velocity −1 −1 account the peculiar velocity due to large-scale structure, vpec = 307 km s and an rms σpec = 230 km s . The adopted range corresponds to the smallest scale constrained by the v = v − v . (1) cosmic CMB pec observational data (∼ 5h−1Mpc). Following Li et al.(2014), we make use of a dark-matter Combining σpec with the error in group velocity the final −1 simulation from the Constrained Local Universe Simulations cosmic velocity is vcosmic = 2924 ± 236 km s . (CLUES) project to find the best estimate of the peculiar 3.4. Hubble distance velocity of NGC 4993. The initial conditions of the sim- ulation were generated using constraints from observations To a good approximation, Hubble’s law gives the luminos- probing the distribution of galaxies and their peculiar veloc- ity distance at low redshift as ities so that the final snapshot is a 3D representation of the   −1 1 − q0 observed local universe (for technical details see Gottloeber Dz = H czcosmic 1 + zcosmic , (2) 0 2 et al. 2010). The simulation reproduces all main large-scale structures around the Local Group such as the Virgo clus- with q0 = −0.53 for standard cosmological parameters. For −1 −1 ter, the Coma cluster, the Great Attractor and the Perseus– a local Hubble constant of H0 = 73.24 ± 1.74 km s Mpc Pisces supercluster, and the resulting peculiar velocity field. (Riess et al. 2016) we obtain Dz = 40.4 ± 3.4 Mpc. The Structures on scales smaller than ∼ 5h−1Mpc are not subject quoted uncertainty accounts for the limited depth of the NGC to observational constraints and emerge from random sam- 4993 group of galaxies (Fig.1) which is less than 1 deg, cor- pling of the initial density field. Their evolution, however, responding to 0.6 Mpc. is strongly influenced by tidal forces from the nearby large- We note that adopting a mean heliocentric velocity of 2995 scale structures. Therefore, the simulation provides a realis- km s−1 for the NGC 4993 group (Kourkchi & Tully 2017) −1 tic and dynamically self-consistent model for the matter dis- would increase vcosmic by 74 km s and hence the inferred tribution and the velocity field in the local universe. The sim- distance by 1.0 Mpc. ulation adopted cosmological parameters from the 3rd data 4. FUNDAMENTAL PLANE DISTANCE release of the WMAP satellite (WMAP3), i.e., matter den- sity Ωm = 0.24, dimensionless Hubble parameter h = 0.73 and The FP is a fairly tight relation between radius, sur- normalization of the power spectrum σ8 = 0.76. The simula- face brightness and velocity dispersion for bulge dominated tion box has a side length equal to 160h−1Mpc and contains galaxies (Djorgovski & Davis 1987). It is widely used as a 10243 particles. distance indicator for early type galaxies. The FP that we In order to find the position vector of the observational consider is cone in the simulation box we calculate the angular distances logR + αlogσ + β loghI i + γ, (3) of NGC 4993 from three reference points which have un- e r e ambiguous analogs in the numerical model: the direction of where Re is the effective radius measured in arcseconds, σ is −1 the peculiar velocity of the Local Group, the Great Attractor the velocity dispersion in km s , hIrie is the mean intensity −2 ( cluster) and the Perseus cluster. These three angu- inside the effective radius measured in L pc , and γ is the lar separations are then used to determine the position vector distance-dependent zero point of the relation. with respect to the analog of the Local Group in the simu- The global values adopted by Jorgensen et al.(1996) are lation box. Having found the direction of NGC 4993 in the α = 1.24 and β = 0.82 in Gunn r. We adopt r = F606W + simulation, we compute the radial projection of peculiar ve- 0.04 mag (Fukugita et al. 1995; Sirianni et al. 2005) for E/S0 locities of dark-matter particles found in a narrow light cone. galaxies. The surface brightness is Figure3 shows the resulting projected peculiar velocity as a function of the recessional velocity measured in the reference loghIrie = − 0.4[F606W(< Re) + 0.04 − 0.29 + 2.5logπ + 5logRe − 26.40 − 10log(1 + z) − 2.5log(1 + z)],

1 using the NED Velocity Correction Calculator at (Hjorth & Tanvir 1997) which transforms into Gunn r https://ned.ipac.caltech.edu/forms/vel_correction.html and corrects for Galactic extinction (Schlafly & Finkbeiner 4 HJORTH ET AL.

Figure 1. Members of the NGC 4993 group of galaxies. Purple circles indicate Kourkchi & Tully(2017) group galaxies with heliocentric velocities less than 3005 km s−1 while blue circles indicate galaxies with heliocentric velocities larger than 3169 km s−1 (see Figure 2). This Digitized Sky Survey image is 2 deg on the side. North is up and East is to the left. DISTANCE TO NGC 4993 5

(1996), their version of the FP excludes biases above 1% due -22.5 to magnitude and morphological selection.

4.1. Fundamental Plane parameters -23.0 4.1.1. Photometric properties The HST F606W image is fit to a 2D version of the Sersic NGC 4993 function, -23.5 ( " 1 #)  R  n I(R) = Ie exp −bn − 1 , (4) Re -24.0

Dec (J2000.0) (deg) where Re is the half-light (effective) radius, Ie is the inten- sity at Re, and n is the fitted Sersic index (which is uniquely -24.5 related to the coefficient bn). 2800 3000 3200 3400 The fit was performed using the Astropy package v (km/s) FittingWithOutlierRemoval (Astropy Collabora- helio tion 2013). The function uses sigma clipping to filter out bad data points and iterates niter times over the data. Data deviat- Figure 2. Heliocentric velocities of the Kourkchi & Tully(2017) ing more than a given sigma from the standard deviation of group PCG1 45466 containing NGC 4993. There is a clear gap in surrounding data points are removed in each iteration. The the heliocentric velocities, with no galaxies in the range 3005–3169 −1 filtered data are then fitted with the Sersic function. km s , i.e., in the positive range above the mean velocity vhelio = 2995 km s−1. Considering only galaxies plotted as filled circles the The fitting algorithm was applied to a range of different −1 mean velocity is vhelio = 2921 km s . The updated velocity of NGC combinations of niter and sigma. Fitting parameters stabilize 4993 is indicated. after 15 iterations for all sigma values, so niter = 15 is used for the fit. The final sigma was chosen for values excluding outliers and including most of the light from the center of the galaxy. This worked well for sigma values between 20–40. 1000 An example of a fit to the F606W image is shown in Fig- ure4. Inside an effective radius of 15.5 00 we find F606=12.99 800 mag. The Sersic index is n ≈ 3.8. A detailed discussion of the properties of NGC 4993 is presented in Levan et al.(2017). 600 4.1.2. Velocity dispersion 400 [km/s] The velocity dispersion is measured in an aperture equiva- pec v 00 200 lent to 3.4 at Coma (Jorgensen et al. 1996), i.e., in a diam- eter of about 8.500 at the distance of NGC 4993. Rather than 0 using aperture corrections from a measurement of the cen- tral velocity dispersion, we measure the luminosity weighted -200 velocity dispersion directly from our MUSE data to be σ = 0 500 1000 1500 2000 2500 3000 3500 171 ± 2 km s−1 (see also Fig. 3 of Levan et al. 2017). This vCMB [km/s] method allows us to tie our FP directly to that of Jorgensen Figure 3. Peculiar velocities (line-of-sight component) of dark mat- et al.(1996), thus eliminating systematic uncertainties re- ter particles in the constrained simulation, in the direction of NGC lated to the approach adopted in determining the velocity dis- 4993. Velocity vCMB is measured in the reference frame of the CMB persion. with respect to an observer located at the Local Group (numerical analog). The red stripe indicates the location of NGC 4993. 4.2. Fundamental Plane distance The parameters used to determine the FP distance are sum- 2011), cosmological surface brightness dimming, and spec- marized in Table1. Calibrating the zero point of the Funda- tral bandwidth (k-correction), with z = zcosmic. mental Plane to the Leo I group we can infer the distance to In principle, FP distances may be affected by the observed NGC 4993 as magnitude range, morphological makeup of the sample (E- − log Re+1.24 log σ−0.82hIri +2.194 2 to-S0 ratio), and environment. According to Jorgensen et al. D = 10 e (1 + zcosmic) Mpc (5) 6 HJORTH ET AL.

Data [mag/arcsec] Filtered data [mag/arcsec] 24.5 24.5

30 30 25.0 25.0

20 20 25.5 25.5

10 10 26.0 26.0

0 26.5 0 26.5 arcsec arcsec

10 27.0 10 27.0

27.5 27.5 20 20

28.0 28.0 30 30

30 20 10 0 10 20 30 30 20 10 0 10 20 30 arcsec arcsec

Fitted Model [mag/arcsec] Residual [mag/arcsec] 24.5 0.5

30 30 0.4 25.0

0.3 20 20 25.5 0.2

10 10 26.0 0.1

0 26.5 0 0.0 arcsec arcsec

0.1 10 27.0 10

0.2 27.5 20 20 0.3

28.0 0.4 30 30

0.5 30 20 10 0 10 20 30 30 20 10 0 10 20 30 arcsec arcsec

Figure 4. Fits to the F606W HST image of the host galaxy of GW 170817. The upper left hand image shows the data, the upper right hand image shows the masked data, the lower left panel shows the fitted model, and the lower right panel shows the residuals between the data and the fit. North is up and East is to the left. DISTANCE TO NGC 4993 7

(Hjorth & Tanvir 1997), where last term corrects the angu- ble constant, the other is based on a direct measurement of lar diameter distance to luminosity distance. The resulting the distance. We can therefore combine the two distances to distance is 44.0 Mpc. yield DNGC4993 = 41.0 ± 3.1 Mpc. This is consistent with the independent distance obtained from the gravitational wave +2.9 Table 1. Redshift and Fundamental Plane parameters of NGC 4993 signal, DGW = 43.8−6.9 Mpc (LIGO & Virgo Collaboration 2017c). 2 −1 zcosmic Re (arcsec) loghIrie (L /pc ) σ (km s ) As a sanity check, we can compare these estimates to a 0.0098 15.5 ± 1.5 2.61 ± 0.06 171 ± 2 slightly different approach. Figure5 presents an analysis of the Cosmicflows-2 catalog (Tully et al. 2013). We took all −1 distances in a range ±1000 km s around vCMB of NGC 4993 and 10 deg around its position. From a fitted linear model 4.3. Uncertainties with intrinsic scatter, the estimated distance for NGC 4993 is D = 36.2 ± 2.8 Mpc (the error includes intrinsic scat- The measurement of the effective radius is subject to sys- CF−2 ter). We tried narrower cuts in v and the position, and the tematic uncertainties due to the nuclear spiral arms, residuals CMB estimate is stable to such variations. This approach is not en- at larger radius, and difficulty in determining the true back- tirely independent from the other electromagnetic distances ground level. We estimate the uncertainty in R is about 1.500. e derived here, so we do not include it in our final value. The combination of R and hI i in Eq. (1) is highly degen- e r e In Section 3, we used the Riess et al.(2016) value of the erate (Jorgensen et al. 1996) and so the FP distance is not Hubble constant, H = 73.24 ± 1.74 km s−1 Mpc−1. Planck very sensitive to the uncertainty in R . The measurement of 0 e Collaboration(2016) instead find H = 67.8 ± 0.9 km s−1 σ required no aperture correction since we integrate over R 0 e Mpc−1 in a ΛCDM cosmology. Using the Planck value for H in the IFU data directly. We estimate that the uncertainty 0 would increase our best estimate of the luminosity distance in σ is about 2 km s−1. The uncertainties in R and σ lead e from 40.4 Mpc to 43.7 Mpc, a shift by ∼ 1σ. This sensitivity to ∼ 2.1% and ∼ 1.4% uncertainties in the distance, in to- to H again points to the utility of these high-accuracy red- tal 2.5% in the distance. Therefore, the observational uncer- 0 shift measurements coupled with independent gravitational- tainties are negligible compared to the intrinsic scatter in the wave distance measurements in the determination of H . FP which amounts to 17% (Jorgensen et al. 1996) which we 0 Other powerful distance indicators to distant early type adopt as the uncertainty in the FP distance to NGC 4993, i.e., galaxies include surface brightness fluctuations (SBF) and D = 44.0 ± 7.5 Mpc. FP the globular cluster luminosity function (GCLF). Both meth- 5. DISCUSSION ods are preferentially used at distances smaller than that of NGC 4993. For example, the GCLF peaks at MV ∼ −7.5 mag F W ∼ . 90 (Rejkuba 2012), i.e., at around 606 25 5 at the distance CF-2 data of NGC 4993, about a magnitude shallower than the 2σ de- 80 best fit distance from redshift tection limit for point sources in the HST image (Levan et al. 70 distance from FP 2017). SBF will be reported elsewhere. Hubble’s law, h=0.7324 60 6. SUMMARY AND CONCLUSIONS 50 D [Mpc] 40 Table 2. Distance estimates to NGC 4993

30 Method Distance (Mpc) 20 Redshift 40.4 ± 3.4 10 2200 2400 2600 2800 3000 3200 3400 3600 3800 4000 4200 Fundamental Plane 44.0 ± 7.5

vCMB [km/s] Combined electromagnetic (z+FP) 41.0 ± 3.1 +2.9 Gravitational waves 43.8−6.9 Figure 5. Distance versus recession velocity in the CMB frame for Cosmicflows-2 galaxies. The best fitting relation is consistent with the distance derived for NGC 4993 from the redshift and FP methods within the errors. We have obtained a new estimate of the distance to NGC 4993 based on two methods. (1) From a new measurement The two distance estimates obtained in this paper are in- of its redshift with MUSE/VLT and an estimate of its pecu- dependent. One is based on a redshift and an assumed Hub- liar velocity using CLUES constrained cosmological simu- 8 HJORTH ET AL. lations we obtain Dz = 40.4 ± 3.4 Mpc. (2) From new mea- The prospects for improving the electromagnetic distance surements of its velocity dispersion (with MUSE/VLT) and to NGC 4993 are good. Future FP studies will benefit from effective radius (HST) (Table1) we obtain a FP distance of studies of larger samples with MUSE/VLT which will allow DFP = 44.0 ± 7.5 Mpc. Combined, these result in an ‘elec- disentangling different velocity components in galaxies and tromagnetic’ distance of DNGC4993 = 41.0 ± 3.1 Mpc. This perhaps a refined FP indicator. The uncertainties in the pecu- compares well with the independent gravitational wave dis- liar velocity will also diminish as we obtain a better under- +2.9 tance of DGW = 43.8−6.9 Mpc. These results are summarized standing of cosmic flows, e.g., along the avenue presented in in Table2. this work. Improved photometric parameters and extension The consistency between the electromagnetic and gravita- to SBF and GCLF distance indicators are within reach with tional wave distances lends credence to both avenues of de- the James Webb Space Telescope. termining distances and confirms that GW 170817 occurred in NGC 4993. There are at least two different ways in which the elec- tromagnetically measured redshift to the host can lead to improved inference from gravitational-wave data. We can R.W. thanks Vinicius Araújo Barbosa de Lima for his assume, as in Section 3.4, that the Hubble constant H0 is help in handling the simulation data. We thank the LIGO known from other observations in order to provide better con- and Virgo Scientific Collaborations for sharing informa- straints on the distance. This distance estimate can then be tion prior to publication and for helpful comments on the fed in to gravitational-wave data analysis as a much tighter manuscript. We thank the anonymous referee for a swift distance prior than the isotropic-in-volume prior p(D) ∝ D2 and insightful review. This work was supported by a VIL- (Veitch et al. 2015). Better distance constraints can help LUM FONDEN Investigator grant to J.H. (project number break some of the strong correlations between parameters 16599). A.J.L. is supported by STFC and the ERC (grant in the gravitational-wave signal, particularly the distance- #725246). J.D.L. gratefully acknowledges support from inclination degeneracy (Cutler & Flanagan 1994; Aasi et al. STFC (ST/P000495/1). I.M. acknowledges partial support 2013; LIGO & Virgo Collaboration 2017c). Tighter con- from STFC. C.G. acknowledges support from the Carlsberg straints on the inclination angle of the binary’s orbit rela- Foundation. This work is based in part on observations tive to the line of sight can in turn aid in the interpretation collected at the European Organisation for Astronomical Re- of the electromagnetic transient. Alternatively, the electro- search in the Southern Hemisphere under ESO programme magnetic redshift measurement and the gravitational-wave 099.D-0668. distance measurement can be combined to achieve an inde- Facilities: Hubble Space Telescope, Very Large Tele- pendent estimate of the Hubble constant, insensitive to the scope potential systematics of electromagnetic distance estimates (Schutz 1986; LIGO & Virgo Collaboration 2017c). Software: ASTROPY (Astropy Collaboration 2013)

REFERENCES

Aasi, J., Abadie, J., Abbott, B. P., et al. 2013, PhRvD, 88, 062001 Hjorth, J., & Tanvir, N. R. 1997, ApJ, 482, 68 LIGO & Virgo collaboration 2017a, PRL, in press Jorgensen, I., Franx, M., & Kjaergaard, P. 1996, MNRAS, 280, 167 LIGO & Virgo collaboration 2017b, ApJL, in press Kourkchi, E., & Tully, R. B. 2017, ApJ, 843, 16 LIGO & Virgo collaboration 2017c, Nature, in press Levan, A. L. et al., ApJL, in press Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, Li, X., Hjorth, J., & Wojtak, R. 2014, ApJL, 796, L4 A&A, 594, A13 Makarov, D., & Karachentsev, I. 2011, MNRAS, 412, 2498 Astropy collaboration, 2013, A&A, 558, 33 Pian, E. et al., Nature, in press Coulter, D. A. et al., Science, in press Rejkuba, M. 2012, Ap&SS, 341, 195 Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 Riess, A. G., Macri, L. M., Hoffmann, S. L., et al. 2016, ApJ, 826, Crook, A. C., Huchra, J. P., Martimbeau, N., et al. 2007, ApJ, 655, 56 790 Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 Cutler, C., & Flanagan, É. E. 1994, PhRvD, 49, 2658 Schutz, B. F. 1986, Nature, 323, 310 Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP, 107, 945 Sirianni, M., Jee, M. J., Benítez, N., et al. 2005, PASP, 117, 1049 Gottloeber, S., Hoffman, Y., & Yepes, G. 2010, arXiv:1005.2687 Sorce, J. G., Tully, R. B., Courtois, H. M., et al. 2014, MNRAS, Hinshaw, G., Weiland, J. L., Hill, R. S., et al. 2009, ApJS, 180, 225 444, 527 DISTANCE TO NGC 4993 9

Springob, C. M., Masters, K. L., Haynes, M. P., Giovanelli, R., & Tully, R. B., Rizzi, L., Shaya, E. J., et al. 2009, AJ, 138, 323 Marinoni, C. 2009, ApJS, 182, 474 Tully, R. B., Courtois, H. M., Dolphin, A. E., et al. 2013, AJ, 146, Springob, C. M., Magoulas, C., Colless, M., et al. 2014, MNRAS, 86 445, 2677 Veitch, J., Raymond, V., Farr, B., et al. 2015, PhRvD, 91, 042003 Tanvir, N. R. et al., ApJL, in press Willick, J. A., Courteau, S., Faber, S. M., et al. 1997, ApJS, 109, Tully, R. B., Courtois, H. M., & Sorce, J. G. 2016, AJ, 152, 50 333