2019ApJ...883L..24S

We can be biased against observing massive black holes (BHs) merging in the local universe as the bounds on the maximum BH mass (M_BH^{max}^) depends on the assumptions regarding the metallicity evolution of the star-forming gas across the cosmic time. We investigate the bounds on the metallicity evolution, mass distribution, and delay times of the binary BH sources based on the 10 observed events by LIGO. We parametrize M_BH^{max}^ to be a function of metallicity that itself is modeled to evolve with redshift in either a modest or rapid fashion. Rapid metallicity evolution models predict a stringent bound of M_BH^{max}⁼⁴⁴-5+9^M_☉, while the bound on M_BH^{max}^ in the models with modest metallicity evolution is M_BH^{max}⁼⁵²-9+16^M_☉. Therefore, inferring M_BH^{max}^ from GW data depends on the assumed metal enrichment history of the universe that is not severely constrained at the moment.