2016ApJ...827....1S

We present results obtained by applying our BAyesian HierArchical Modeling for the Analysis of Supernova cosmology (BAHAMAS) software package to the 740 spectroscopically confirmed supernovae of type Ia (SNe Ia) from the "Joint Light-curve Analysis" (JLA) data set. We simultaneously determine cosmological parameters and standardization parameters, including corrections for host galaxy mass, residual scatter, and object-by-object intrinsic magnitudes. Combining JLA and Planck data on the cosmic microwave background, we find significant discrepancies in cosmological parameter constraints with respect to the standard analysis: we find {Ω}_{m}=0.399±0.027, 2.8σhigher than previously reported, and w=-0.910±0.045, 1.6σhigher than the standard analysis. We determine the residual scatter to be σ_{res}=0.104±0.005. We confirm (at the 95% probability level) the existence of two subpopulations segregated by host galaxy mass, separated at log_{10}(M/{M}☉)=10, differing in mean intrinsic magnitude by 0.055±0.022 mag, lower than previously reported. Cosmological parameter constraints, however, are unaffected by the inclusion of corrections for host galaxy mass. We find ∼4σevidence for a sharp drop in the value of the color correction parameter, β(z), at a redshift {z}_{t}=0.662±0.055. We rule out some possible explanations for this behavior, which remains unexplained.