Arxiv:2004.14499V2 [Astro-Ph.GA] 21 May 2020 Spiral Galaxies; Inclination; Galaxy Structure; Large-Scale Structure of the Universe;
Total Page:16
File Type:pdf, Size:1020Kb
Draft version May 22, 2020 Typeset using LATEX preprint2 style in AASTeX62 Cosmicflows 4: The Calibration of Optical and Infrared Tully-Fisher Relations Ehsan Kourkchi,1 R. Brent Tully,1 Gagandeep S. Anand,1 Hel´ ene` M. Courtois,2 Alexandra Dupuy,2 James D. Neill,3 Luca Rizzi,4 and Mark Seibert5 1Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA 2University of Lyon, UCB Lyon 1, CNRS/IN2P3, IUF, IP2I Lyon, France 3California Institute of Technology, 1200 East California Boulevard, MC 278-17, Pasadena, CA 91125, USA 4W.M. Keck Observatory, 65-1120 Mamalahoa Highway, Kamuela, HI 96743, USA 5The Observatories of the Carnegie Institute of Washington, 813 Santa Barbara Street, Pasadena, CA 91101, USA (Received 2020 February 4; Revised 2020 April 9; Accepted 2020 April 29) Submitted to ApJ ABSTRACT This study is a part of the Cosmicflows-4 project with the aim of measuring the dis- tances of more than ∼10,000 spiral galaxies in the local universe up to ∼15,000 km s−1 . New H i line width information has come primarily from the Arecibo Legacy Fast ALFA Survey. Photometry of our sample galaxies has been carried out in optical (SDSS u, g, r, i, and z) and infrared (WISE W1 and W2) bands. Inclinations have been determined using an online graphical interface accessible to a collaboration of citizen scientists. Galaxy distances are measured based on the correlation between the rotation rate of spirals and their absolute luminosity, known as the Tully-Fisher Relation (TFR). In this study, we present the calibration of the TFR using a subsample of ∼600 spirals located in 20 galaxy clusters. Correlations among such observables as color, surface brightness, and relative H i content are explored in an attempt to reduce the scatter about the TFR with the goal of obtaining more accurate distances. A preliminary determination of the Hubble constant from the distances and velocities of the calibrator clusters is −1 −1 H0 = 76:0 ± 1:1(stat.)±2:3(sys.) km s Mpc . Keywords: Distance measure; Galaxy distances; Galaxy photometry; H i line emission; arXiv:2004.14499v2 [astro-ph.GA] 21 May 2020 Spiral galaxies; Inclination; Galaxy structure; Large-scale structure of the universe; 1. INTRODUCTION Through the Cosmicflows program, the dark matter distribution in the universe is recon- structed based on the peculiar motions of galax- [email protected] ies in response to the underlying gravitational [email protected] field. The peculiar velocities of galaxies are cal- culated by subtracting the Hubble expansion 2 Kourkchi et al. rate at their position from their observed radial cern of the next addition, Cosmicflows-4, that velocities. Therefore, the precise measurement this imbalance in sky coverage be redressed. of distances is vital in this analysis. This new edition will contain updates across To date, three editions of Cosmicflows have a broad spectrum of inputs, but most of all, been published (Tully et al. 2008, 2013, 2016). there will be a major enhancement in the num- Throughout, the samples are heterogeneous in ber of spiral luminosity-rotation rate measure- that there is input involving multiple method- ments made available as a consequence of the ologies and multiple sources. The scale founda- large number of neutral hydrogen line profile de- tions are set by the Cepheid Period-Luminosity tections resulting from the Arecibo Legacy Fast Relation (Leavitt & Pickering 1912) and the ALFA Survey (Haynes et al. 2018, ALFALFA; Tip of the Red Giant Branch (TRGB) stan- ) covering the declination range 0 < δ < +38. dard candle (Lee et al. 1993) supplemented by Photometry for these galaxies is alternatively methods that have high precision but limited available from the Sloan Digital Sky Survey applicability such as Detached Eclipsing Bina- (SDSS; Alam et al. 2015) and the Wide-field In- ries (Pietrzy´nskiet al. 2019) and the nuclear frared Satellite Explorer (WISE; Wright et al. maser in NGC 4258 (Humphreys et al. 2013). 2010). The present paper presents up-to-date These local calibrators constrain the scaling of calibrations of the luminosity-rotation rate cor- methods useful at intermediate distances, such relations in the bands of these surveys. as the amplitudes of surface brightness fluctua- Since its inception (Tully & Fisher 1977; tions in early-type galaxies (Tonry & Schneider Aaronson et al. 1979), the luminosity-rotation 1988), the three-parameter fundamental plane rate relation (hereafter TFR) has been refined correlation also applicable to early types (Djor- and reevaluated by many users. If the inter- govski & Davis 1987; Dressler et al. 1987), and est is in galaxy distances, a primary concern is the luminosity-rotation rate relation found to the control of bias (Willick 1994). The current hold for spiral galaxies (Tully & Fisher 1977). analysis makes use of the correlation fit with er- All of these procedures must be shown to be rors taken in line widths, a procedure that has consistent with each other and, in turn, with been called the \inverse" TFR (ITFR). This the methodology that extends to large distances procedure was initially discussed as a general founded on the predictability of the peak lumi- principle by Schechter(1980) and specifically nosities of supernovae of type Ia (Phillips 1993). with relevance to the TFR by Tully(1988). Cosmicflows-3 (CF3), the last edition report- The issue was discussed in detail in x3 of Tully ing galaxy distances, contains information on & Pierce(2000). The slope of the power-law 17,647 systems1. Roughly half the contributions fit of the ITRFR nulls the bias to a first ap- in that compilatfion come from a single source: proximation. Residual bias can be taken into Fundamental Plane measures from the Six De- account with a small correction to be discussed. gree Field Redshift Survey (Magoulas et al. It deserves emphasis that a focus on the cor- 2012; Springob et al. 2014). This contribution relation with errors taken in line width, the is restricted to δ ≤ 0, the celestial south. As distance-independent parameter, is desired if a consequence, the coverage in CF3 is strongly the primary goal is the measurement of dis- tilted toward the south. It is a particular con- tances, but it is not the appropriate procedure if the interest is a physical understanding of the 1 The CF3 catalog updated with minor corrections is relationship between luminosity and rotation in available at http://edd.ifa.hawaii.edu spiral galaxies. In the latter case, it is more in- Multi-wavelength calibration of the Tully-Fisher Relations 3 formative to make a bivariate fit taking into ac- candidate galaxies in optical (SDSS u, g, r, i count errors affecting both line widths and mag- and z) and/or infrared WISE W1 (3.4µm) and nitudes (Giovanelli et al. 1997; Masters et al. W2 (4.6µm) bands. The purpose of the cur- 2006; Ponomareva et al. 2017; Lelli et al. 2019). rent study is to provide the calibration of the We employ the inverse relation, since our in- TFR at these wave bands to measure distances terest is the measurement of galaxy distances. for our full sample of galaxies. To achieve our Recalibrations have been carried out in both op- goal, we use a sample of 648 calibrator galax- tical and infrared bands in the course of the ies (roughly double that available previously) Cosmicflows program. Tully & Courtois(2012) in 20 clusters (up from 13). We explore intrin- calibrated the TFR at the I band using 267 spi- sic relations between color, morphology, and rals in 13 clusters. Aaronson et al.(1979) sug- other galaxy observables with deviations from gested the use of near-infrared bands to improve the TFR at different wave bands, toward possi- the TFR methodology, given the lesser dust ob- ble ways to reduce scatter and hopefully obtain scuration in host galaxies. Furthermore, the more accurate distances. light of old stars peaks at longer wave bands; We explain our sample selection criteria and therefore, near-infrared luminosities may better present our data product in x2. In x3, we cali- represent the baryonic mass of galaxies. With brate the luminosity-line width correlations. In the motivation of reducing the scatter about x4, we explore the possibility of reducing the the TFR and getting more accurate distances, scatter of the calibrated relations using addi- Sorce et al.(2013) extended the calibration of tional observables such as color indices and sur- the TFR toward the infrared using Spitzer 3:6 face brightnesses. The radial velocities of the µm Infrared Array Camera (IRAC) imaging for calibrator clusters are determined in x5, and the same sample of galaxies as used by Tully and their measured distance moduli are evaluated Courtois. Using the available I-band magni- across all passbands and calibrations in x??. We tudes, they also observed a correlation between provide estimations of the Hubble constant, H0, the optical-infrared color index and the devi- on the basis of the clusters in x7. Our conclusion ation of galaxy magnitudes from the infrared is presented in x8. TFR. They found a tighter TFR after adjust- ing [3:6] magnitudes for the effect of an I −[3:6] 2. DATA color term. 2.1. Calibrators Neill et al.(2014) followed the same procedure using the WISE (Wright et al. 2010) and added To calibrate the luminosity-line width correla- 43 galaxies to the previous sample. They also tion in spirals, we follow the same methodology updated the calibration of the I-band TFR us- that was originally described by Tully & Cour- ing their larger sample. They studied the pos- tois(2012) and further implemented by Sorce sibility of decreasing the scatter about the in- et al.(2013) and Neill et al.(2014).