A sample of dust attenuation laws for Dark Energy Survey supernova host

Context. Type Ia supernovae (SNe Ia) are useful distance indicators in cosmology, provided their luminosity is standardized by applying empirical corrections based on light-curve properties. One factor behind these corrections is dust extinction, which is accounted for in the color–luminosity relation of the standardization. This relation is usually assumed to be universal, which can potentially introduce systematics into the standardization. The “mass step” observed for SN Ia Hubble residuals has been suggested as one such systematic. Aims. We seek to obtain a more complete view of dust attenuation properties for a sample of 162 SN Ia host galaxies and to probe their link to the mass step. Methods. We inferred attenuation laws toward hosts from both global and local (4kpc) Dark Energy Survey photometry and composite stellar population model ﬁts. Results. We recovered a relation between the optical depth and the attenuation slope, best explained by di ﬀ ering star-to-dust geometry for di ﬀ erent galaxy orientations, which is signiﬁcantly di ﬀ erent from the optical depth and extinction slope relation observed directly for SNe. We obtain a large variation of attenuation slopes and conﬁrm these change with host properties, such as the stellar mass and age, meaning a universal SN Ia correction should ideally not be assumed. Analyzing the cosmological standardization, we ﬁnd evidence for a mass step and a two-dimensional “dust step”, both more pronounced for red SNe. Although comparable, the two steps are not found to be completely analogous. Conclusions. We conclude that host galaxy dust data cannot fully account for the mass step, using either an alternative SN standardization with extinction proxied by host attenuation or a dust-step approach.


Introduction
Type Ia supernovae (SNe Ia) are a set of astrophysical transients that make for excellent distance indicators, having famously contributed to the discovery of the accelerated expansion of the Universe (Riess et al. 1998;Perlmutter et al. 1999).SN Ia data have proven extremely important in constraining the various cosmological parameters, such as the equation of state parameter for nous.Even so, after the application of these two corrections, SNe Ia originating in higher-mass galaxies appear consistently over-luminous (e.g., Sullivan et al. 2010;Lampeitl et al. 2010;Kelly et al. 2010), such that a mass step is evident in the difference between their observed and predicted luminosities.A third correction is often introduced to account for this environmental dependency.
Because of the empirical nature of these corrections, the full impact of the standardization process on cosmology is still not fully understood.Although heavily debated, two main effects seem to be at the origin of the color-luminosity relation: intrinsic color variations among different SNe Ia and the effects of dust extinction, namely reddening (Conley et al. 2007).
Extinction refers to the absorption of light as well as the scattering of that light out of the line of sight, which, as a result of the overall composition, size, and orientation of the dust grains, tend to be stronger for bluer wavelengths (Calzetti et al. 1994;Calzetti 2001; Salim & Narayanan 2020).This leads both to a dimming and a reddening of the observed spectrum.The reddening law is most commonly parameterized by the total-to-selective extinction parameter R V = A V /E(B − V), where A V is the extinction for the V band and E(B − V) = A B − A V is the difference between the object's observed and intrinsic B − V color index, referred to as the color excess.
For the purposes of SN Ia standardization, the colorluminosity relation is generally assumed to be the same for all SNe Ia in the population under study (Tripp 1998).By the same token, when correcting for galactic dust reddening, it is common practice to assume a constant R V for all galaxies, such as the Milky Way average value of R V = 3.1 (e.g., Schlafly et al. 2016;Cardelli et al. 1989a).Other values exist, such as the Calzetti et al. (2000) value of R V = 4.05 for star-forming galaxies and the peculiarly low values for the Magellanic Clouds (e.g., Lequeux et al. 1982;Nandy et al. 1981).However, it has been shown that dust properties, particularly R V , can vary between galaxies (e.g., Salim & Narayanan 2020;Salim et al. 2018;Zahid et al. 2013), and so assuming a common dust law for all galaxies and even for different regions within the same galaxy might introduce systematic errors into the recovered results (e.g., Draine 2003;Nataf et al. 2015).
In the case of SN Ia standardization, it has been suggested that the previously mentioned mass step is one such systematic.Brout & Scolnic (2021) find that the correlation between the corrected luminosity and the host galaxy mass can be explained by differing R V values between low-and high-mass galaxy bins, while Wiseman et al. (2022) find that the variation of R V with galaxy age can explain almost the entire observed mass step.Additionally, Johansson et al. (2021) and Meldorf et al. (2022) show that the mass step can be eliminated using individual lowredshift SN Ia light-curve fits in the optical/NIR range to directly determine extinction properties and thus constrain the colorluminosity relation for each SN.
However, there is also evidence to suggest that the mass step cannot entirely be explained as a dust systematic.Thorp et al. (2021) find that a mass step is still observable even when allowing for different R V values between high-and low-mass host galaxies.González-Gaitán et al. (2021) obtain similar results when considering separate color-luminosity relations for different SN populations.Uddin et al. (2020) and Ponder et al. (2021) interpret the consistency of the mass steps recovered from both optical and NIR light curves as evidence that dust does not play a large role in the mass-step correlation, as the effects of dust extinction in NIR are greatly reduced.Jones et al. (2022) also find evidence for a NIR mass step when looking at combined high and low redshift populations.
More thorough studies on the nature of the SN Ia colorluminosity relation are needed to better grasp the contributions of interstellar dust to the mass-step correction.As shown by Meldorf et al. (2022) and Popovic et al. (2021), it can prove advantageous to individually determine dust properties for each individual SN Ia or for each individual SN Ia population in order to better constrain the color-luminosity relation in the standardization.This allows one to account for varying dust properties across different lines of sight and environments.However, it can be challenging to directly determine the extinction for SNe Ia, especially at higher redshifts, which require deep infrared photometry.The study of SN Ia host galaxies offers an interesting opportunity, not only to tackle the problem of dust-induced systematics in SN standardization, but also to extend our knowledge about general dust properties affecting extra-galactic observations.
When dealing with extended objects, such as SN host galaxies, the description of dust effects by extinction becomes too simplistic, as we can no longer presume to be dealing with a single line of sight.In these cases, we refer to the analogous concept of dust attenuation (Calzetti et al. 1994;Calzetti 2001;Salim & Narayanan 2020).This phenomenon includes not only the abovementioned effects, but also effects arising from the spatial distribution of stars and dust in a galaxy or stellar population, such as scattering of light back into the line of sight, multiple dust cloud densities, and light emitted by unobscured stars.These phenomena are mostly negligible when considering only a point source.
In this work, we tackle two main questions.First, we study dust attenuation laws for the host galaxies of a set of SNe Ia, examining how they differ from each other and how they compare to point source literature extinctions measured directly from SNe Ia.To this effect, we present an alternative approach to infer dust attenuation laws toward the host galaxies for a cosmological sample of 162 SNe Ia collected as part of the Dark Energy Survey (DES) (Abbott et al. 2016;Kelsey et al. 2021).We fit composite stellar population (CSP) models to the "deep" host photometric data in the griz filter bands1 (Smith et al. 2020;Wiseman et al. 2020), as detailed in Section 2. To validate our approach, in Section 3 we present the recovered fit results for a series of simulated host galaxies, comparing the best-fit values and the "true" parameters used in the simulations.In Section 4 we analyze and characterize the recovered attenuation laws, comparing them to known SN extinction laws.Namely, we analyze the viability of using our host galaxy fit results to approximate the extinction for the respective SNe.Second, we seek to ascertain whether the mass step for SNe Ia can be explained by host dust properties.In Section 5 we explore the possibility of a step in luminosity linked to the dust environment surrounding the SNe and whether it can be related to the more usual mass step.We also explore the relation between these Hubble residuals steps and SN color.We assume a spatially flat ΛCDM model (Condon & Matthews 2018), with a matter density Ω m = 0.3, a dark energy density of Ω Λ = 0.7, a dark energy equation of state parameter w = −1, and a Hubble constant H 0 = 70km s −1 Mpc −1 .

Dark Energy Survey photometric data
The Dark Energy Survey (DES) is an imaging survey covering ∼ 5100 square degrees of the Southern Hemisphere using the 4-m Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO), equipped with the 520 megapixel wide-field Dark Energy Camera (DECam) (Flaugher et al. 2015) with a 0.263 arcsecond/pixel resolution.
In this work we look at data collected as part of the first threeyear cosmological sample of the SN survey (DES3YR)2 (Brout et al. 2019), consisting of 206 spectroscopically confirmed SNe, as well as their respective host galaxies (Smith et al. 2020;Wiseman et al. 2020).These data include both global and local broadband photometry for the hosts in the DECam griz filter bands, originally computed by Kelsey et al. (2021).Only the 162 hosts meeting the selection cuts of Kelsey et al. (2021) were considered, covering a redshift range of 0.077 < z < 0.58.In addition, we use the light-curve parameters recovered for each of the corresponding SNe Ia by Brout et al. (2019).
As detailed in Kelsey et al. (2021), the global photometry for each host galaxy was measured from stacked images using SExtractor (Bertin, E. & Arnouts, S. 1996), employing Kron FLUX_AUTO measurements (Wiseman et al. 2020).A detection image was used to set the aperture, ensuring it was the same for the measurements in each filter.The images were also corrected for Milky Way dust extinction using Schlegel dust maps (Schlegel et al. 1998) and a Fitzpatrick extinction law (Fitzpatrick 1999) with standard R V = 3.1 and fiducial coefficients R g = 3.186, R r = 2.140, R i = 1.569 and R z = 1.196 (Kelsey et al. 2021;Abbott et al. 2018).
In the case of local photometry, the aperture photometry tool from the Photutils Python module (Bradley & et al. 2019) was used to take measurements in a 4 kpc radius aperture around the SN location.This value is chosen to account for the 1σ seeing of the DES seeing-optimized stacks, which is 0.55 ′′ (Kelsey et al. 2021).Thus, it is ensured that the aperture is larger than a 0.55 ′′ radius in our redshift range.The resulting fluxes were corrected for Milky Way extinction in the same way as described above.As the local aperture size is constant, the local stellar mass obtained in this case can be thought of as a measure of the local density.

Host galaxy fits
We are interested in fitting CSP models to the photometric data for the SN Ia host galaxies.We do this by computing the Spectral Energy Distributions (SEDs) for each CSP and comparing them with the observed data.We rely on Bayesian inference methods to do so, using the Prospector3 (Johnson et al. 2021) and emcee4 (Foreman-Mackey et al. 2013) Python packages.The models themselves are built using the FSPS5 (Conroy et al. 2009;Conroy & Gunn 2010) FORTRAN code and the python-FSPS6 package (Foreman-Mackey et al. 2014).
Our fitting model is built using the following assumptions: a Kroupa (2001) initial mass function (IMF); a delayed τ-model star formation history (SFH), with an exponentially decreasing star formation rate SFR(t) = te −t/τ , with e-folding time τ = 1Gyr and time t measured up to the stellar population age t age7 ; a modified Calzetti attenuation law (Noll, S. et al. 2009;Kriek & Conroy 2013), with an expression given by Eq. 1: The metallicity log(Z ⋆ /Z ⊙ ), the age t age , the optical depth τ V , and the dust index n are taken as model free parameters, with overall flat priors.The stellar mass is fixed, referring to the fit values obtained by Kelsey et al. (2021).
In Eq. 1, k(λ) describes the original Calzetti et al. (2000) attenuation law with an original R V,0 = 4.05 and D(λ, n) describes a UV-bump centered at λ 0 = 2175Å (Noll, S. et al. 2009).We adopt this attenuation law both because of its ability to describe the attenuation for a wide range of star-forming galaxies (Noll, S. et al. 2009) and because it allows for a finer probing of the reddening levels by offering a two parameter description.τ V describes the optical depth of the medium, such that A V ≃ 1.086τ V , while n is introduced as a dust index parameter in order to produce different slopes without the need to alter τ V ; n is constrained between the values of −2.2 and 0.4, with more negative values corresponding to a larger level of reddening, for the same value of τ V .The more commonly used R V parameter can be obtained using With regard to the model parameters, the stellar metallicity does not vary with age in our model, as no age-metallicity relation is implemented (Leja et al. 2017).In addition, the age parameter t age for a given population is constrained such that no time-steps older than the universe at the corresponding redshift are allowed.
The real SFH for the galaxies may be much more intricate than the simplistic delayed τ-model we assume.However, given how few filter bands we have to work with, more complex models are not feasible in this case.The effect of the assumptions made for the SFH and dust models on the recovered parameters will be further explored in the next section.
The choice to fix the stellar mass in the fits is made because, although Bayesian inference does not limit the number of fit parameters one can use, increasing their number would increase the number of iterations necessary for convergence, making the computation time much longer.While other parameters might exhibit some degeneracy, the stellar mass is fairly easy to constrain, as it represents a normalization of the SED.As such, the mass values obtained from a five-parameter sampling do not significantly differ from those previously obtained for these galaxies by Kelsey et al. (2021), with test fits exhibiting a median relative difference of ∼ 4% for both global and local measurements.A deeper look at these results can be found in Appendix A. The differences between a four-and five-parameter sampling on the recovered global dust parameters are also not very noticeable, with test fit median relative differences of ∼ 1% and ∼ 4% for τ V and n, respectively.These differences are more pronounced when using local measurements, with median relative differences of ∼ 8% and ∼ 60% for τ V and n, respectively.These relatively high values are mainly a product of convergence problems in the four-parameter sampling process.
As a final note, although the posterior distributions obtained from Bayesian inference are very informative, a best-fit value is often useful for data analysis.For a given galaxy and a given parameter, we take this value to be the median of the distribution, with the uncertainty range given by the 16th and 84th percentiles, which define a 68% credible region around the median.

Host galaxy simulations
In order to establish the reliability of our fitting methodology, a number of tests were conducted.In each of these, various fits were performed for a population of simulated host galaxies.These galaxies were generated by drawing the values for M ⋆ , Z ⋆ , t age , τ V , n, and z from a series of random flat distributions.While this method does not guarantee that the simulated population will resemble an actual SN Ia host population, it is adequate for testing the performance of the fitting method across the full parameter space.
Photometry in the griz filter bands for each of these galaxies was generated using FSPS.Gaussian noise was also applied to the photometric data, in accordance with a prescribed signal-tonoise ratio (S/N) of 10, which was chosen to roughly reflect the minimum S/N of the DES data set (Kelsey et al. 2021).

Fit photometry tests
With the first set of tests, we seek to ascertain whether griz photometry can be used to accurately recover attenuation parameters.For this purpose, 150 galaxies were simulated employing the same CSP model as described in Section 2. We define this as the default simulation.For these fits, both the stellar mass and the redshift were fixed to their "true" values, following the rationale detailed in Section 2.2.This represents the best case scenario for our fitting method.As a term of comparison, additional photometry in the NUV/FUV GALEX (Martin et al. 2005) and JHKs 2MASS (Skrutskie et al. 2006) filter bands was generated for half of these galaxies and used in addition to the griz photometry in a new set of fits.Some precision and accuracy metrics for the dust parameters (τ V and n) recovered in each of the fits, including the median differences between the simulated and fitted parameters (median ∆) and the corresponding median relative error, are shown in Table 1.We find that the fits using additional photometry do not vary greatly from their griz counterparts.In particular, the values for τ V and n exhibit median relative differences of ∼ 0.1% and ∼ 2%, respectively.The full comparison between the parameters recovered for both sets of fits can be found in Fig. B.2.
Regarding the role of dust emission from far infrared (FIR) data on the determination of attenuation parameters τ V and n, we find from analogous simulations and fits that its inclusion results in a similar improvement to the one discussed above.
It is also crucial to establish that the fitting methodology does not impart any nonphysical correlation unto the recovered parameters.For the griz test fit, the Spearman correlation coefficient for τ V and n is found to be 0.038, meaning that no correlation is found between these two parameters.Given that the mock galaxy set was generated at random, without imposing any prior relation between the model parameters, the lack of correlation found indicates that no relation is being artificially introduced by the fitting process.

Fixed mass error tests
With the second set of tests, we explore how a possible variation in the fixed value of stellar mass affects the recovered fit parameters.Photometry in the griz filter bands for the same simulated galaxies used in the previous set of tests was used, with the fixed value of the mass being varied by ±4%, according to the median difference between Kelsey et al. (2021) and our own inferred masses, recovered in Section 2.2.The median difference was chosen to vary the value of the fixed mass because, as seen from Fig. A.1, the bias between our fits and the mass values from Kelsey et al. (2021) is the dominant source of the deviations between the two, even more so than the scatter of the fit.
Once again, precision and accuracy metrics for the dust parameters recovered in each of the fits are shown in Table 1.These results demonstrate that the fixing of the stellar mass does not compromise the determination of the other fit parameters.

Alternative model simulations tests
For the third set of tests, new simulations were performed using different models from those assumed in the fitting process, with the goal of understanding how the methodology would cope with these variations.Six different simulation models were tested, with 75 galaxies being simulated for each case.Each of the models used in these tests differs from the model presented in Section 2 in one of the following ways: a Cardelli et al. (1989b) attenuation law, with varying R V instead of n; a delayed τ-model SFH with τ = 0.1; a delayed τ-model SFH with τ = 10; a delayed τ-model SFH truncated at t trunc = 7.5 Gyr8 ; a delayed τ-model SFH truncated at t trunc = 5 Gyr; a delayed τ-model SFH truncated at t trunc = 3 Gyr; The fits were once again performed using griz photometry, with the default fitting CSP model defined in Section 2. Precision and accuracy metrics for the dust parameters recovered in each of the fits are again displayed in Table 1.This final test set shows that, for the most part, our methodology is able to cope with variations in the model used for simulation without too much impact on the dust parameters.
The largest changes are found, as expected, when different attenuation laws are introduced, which cannot be avoided due to the different parameterizations used.In this case, R V is the most affected parameter, with a constant bias offset of 0.261 present in the recovered fit parameters.A more detailed look at this particular case is offered in Appendix B.
Overall, although all the simulation tests show discrepancies in the recovered stellar age and metallicity parameters (see Appendix B), the fits show a remarkable recovery of the original dust parameters, as seen in Table 1.Across the board, we find relatively low root mean square (RMS) values for these parameters, as well as median relative errors of ∼ 10% or lower when comparing to the true simulation values, suggesting a high degree of accuracy.In addition, we find that the respective median fit errors are, in general, conservative, as they are very close to the dispersion (MAD) of the fitted parameters.There are, however, some biases observed, the most significant of which is the bias of 0.263 in τ V observed for the default simulation test.There is also a large bias of 0.261 in n in the Cardelli test, which is most likely due to the different parameterization for different attenuation laws.Further discussion of the various tests can be found in Appendix B.
In general, it is thus reasonable to conclude that, although our fitting method might struggle to break the dust-metallicityage degeneracy and accurately determine all free fit parameters, resulting in not very robust non-dust parameters, it is sufficient to Notes.The median fit error was calculated based on all the error values for each of the fit parameters.As the fit error bars are asymmetric, the largest of the median values was taken in each case.The median ∆, median absolute deviation (MAD) and median relative error were calculated from the difference between the best-fit result and the "true" value for each parameter.The percentage of outliers was calculated by assuming a normal distribution for the data and identifying the points more than 2 standard deviations away from the mean.
satisfactorily infer and map dust properties for SN host galaxies.
As such, we can be fairly confident in the results presented and discussed in Section 4.

DES SN Ia host galaxies
In this section, we present and discuss the results of the CSP fits to the DES SN Ia host galaxies9 , placing particular emphasis on the recovered attenuation laws.The fits were performed following the methodology described in Section 2.2 and using both global and local photometric data in the griz filter bands, originally computed by Kelsey et al. (2021).Additionally, we compare the recovered host attenuation laws with literature low redshift SN Ia extinction laws.

Dust attenuation for DES galaxies
The central aspect in the characterization of the recovered attenuation laws is the correlation between the best-fit values for the τ V and n parameters, which are plotted in Fig. 1.A Gaussian Process Regression and its respective 68% credible interval, performed using the GauPro10 (Erickson 2021) R package, are also plotted to better illustrate the observed correlation.11Despite some scatter, there is a clear correlation between the two quantities in both the global and the local cases, with larger optical depths corresponding to larger values of n, and thus shal-lower attenuation curves.The Spearman correlation coefficient between these two parameters for the DES data set is found to be 0.671 for the global measurements and 0.777 for the local measurements.
Comparing the global and the local fit results offers an interesting insight.As seen in Fig. 2, where we plot the deviations ∆τ V = τ V Global − τ V Local and ∆n = n Global − n Local as a function of the redshift z, the results for both fits appear fairly consistent, with median differences of ∆τ V = 0.26 and ∆n = 0.026.
There are some differences between the fits, as the values obtained for τ V tend to be larger for the global photometry.There is also a large level of scatter present in the plots, with standard deviations of 0.712 for ∆τ V and 0.575 for ∆n.Furthermore, the percentage of galaxies for which the best-fit values are not compatible with each other is found to be 39.8% for τ V and 36.0%for n.
The consistency of the observed correlation between the two result sets indicates that an aperture of 4 kpc is not sufficiently small to negate the geometrical effects associated with attenuation (see next section).Likely for the same reason, Kelsey et al. (2021) observe a similar consistency between global and local measurements, particularly when looking at U-R steps.In fact, the global and the local mass values obtained by Kelsey et al. (2021) indicate that, in some cases, a 4 kpc aperture is actually as large as the galaxy itself.This means one needs even smaller physical apertures to correctly estimate dust extinction at the SN location.However, this is difficult to achieve in practical terms due to instrument limitations.The fact that the results remain constant as the redshift increases indicates that the diminishing angular size of the aperture does not, in this case, play a large role in the results obtained.
We used either complete or partial additional global photometry in the NUV/FUV GALEX (Martin et al. 2005) and JHKs 2MASS (Skrutskie et al. 2006) filter bands for 33 of the studied galaxies.Although this additional set of measurements is not enough to draw any conclusions by itself, it is useful to confirm our previous fit results.Combining the new photometry with the griz global data points, a new set of fits was performed for these select galaxies.The differences between the new fit values for τ V and n and the previously discussed griz results exhibit median relative values of ∼ 0.5% and ∼ 3%, with a standard deviation of 1.14 and 0.40, respectively.Despite the somewhat large value of scatter for τ V , these values confirm the reliability of the results obtained with griz measurements.

Attenuation laws and galaxy orientation
For comparison, the mean relation obtained by Chevallard et al. (2013) for the overall correlation between τ V and n for a set of simulated galaxies observed at different orientation angles is also shown in Fig. 1.This relation, which results from a slightly different attenuation law than the one discussed in this paper, differs somewhat from our data, with the simulated galaxies appearing in general redder.Even so, the overall qualitative tendency is the same as the one found for DES galaxies.Additionally, the recovered tendency also matches the behavior observed by both Leja et al. (2017) and Narayanan et al. (2018), which further suggests we are indeed recovering a physical relation.Given the agreement with the Chevallard et al. ( 2013) mean relation, it is reasonable to conclude that the overall correlation between the recovered dust parameters can be well explained by different galaxy orientations, as also described by Viaene et al. (2017) and Narayanan et al. (2018) and schematically illustrated on Fig. 3.This same trend between τ V and n is also present for elliptical galaxies, as recently shown by Sachdeva & Nath (2022).
As stated by those works, cases with small τ V , would mostly correspond to galaxies that are being observed face-on.In this case, we must consider both the photons emitted parallel and perpendicularly to the galactic plane.Photons emitted perpendicularly to the galactic plane (along the line of sight for a face-on galaxy) suffer minimal attenuation, independent of their wavelength.Additionally, we must also take into account photons emitted along the equatorial plane of a galaxy, which might be scattered away from the plane and into the line of sight.On one hand, blue photons tend to be scattered more forwardly, meaning they tend to remain confined to the galactic plane, eventually being fully absorbed (Chevallard et al. 2013).On the other hand, red photons tend to be scattered more isotropically and have a much lower interaction cross-section with dust, allowing them to more consistently escape the galactic plane (Chevallard et al. 2013).This means that an excess of red radiation is added to our line of sight, leading to higher values of reddening, reflected in more negative values for n.
The cases with a larger τ V statistically correspond to an edgeon view of the galaxies.In these cases, radiation emitted from the deepest layers of the galaxy is greatly absorbed, independently of wavelength.The radiation reaching an observer is thus dominated by stars located in the outermost layers of the galaxy, unobscured stars and light scattered into the line of sight (Narayanan et al. 2018).This leads to an overall lower level of reddening and values of n closer to 0.
Galaxy orientation is not the only mechanism behind the correlation observed in Fig. 1, being mostly relevant for spiral galaxies.For elliptical galaxies, galaxy orientation is not an important factor.However, due to their general elliptical shape and overall low dust content, these particular galaxies end up exhibiting a behavior similar to the one observed for head-on spiral galaxies, leading to an overall large level of reddening.Additionally, different galaxies can have different intrinsic optical depths, which is one of the possible reasons for the scatter observed in the figure and the relations discussed in the next sections.It has also been recorded that there is some degree of degeneracy between galaxy geometry and intrinsic dust composition, which somewhat muddles the analysis (Viaene et al. 2017).

Attenuation laws and stellar mass
Another important factor in the characterization of the recovered attenuation laws is the possible correlations of the dust parameters with the stellar mass and age, which are plotted in Fig. 4. Despite some scatter, the relation between τ V and log(M ⋆ /M ⊙ ) seems to display two separate behaviors depending on the age of the corresponding stellar population.For younger galaxies, there is a steady increase of τ V with log(M ⋆ /M ⊙ ), which is in agreement with results obtained for star-forming galaxies by  2013), as more massive galaxies have more gas and dust content with more active star formation.A separate tendency is observed for older galaxies, which tend to exhibit both larger masses and smaller optical depths.This matches what is expected for quiescent galaxies (Salim et al. 2018;Zahid et al. 2013), which have less dust and make up a large fraction of high-mass galaxies (Peng et al. 2010).
The behavior of the two age populations is most prominent when looking at the global mass, but it is also apparent in the local case.This suggests a link not only with the total stellar mass of the galaxy, but also with the local stellar density, with older galaxies being generally more dense.It is also interesting to notice that a larger stellar density does not necessarily equate to a larger dust column density, as seen from the relatively low values of τ V for high local masses.Although the large level of scatter complicates the analysis of the relation between n and log(M ⋆ /M ⊙ ), we find that older galaxies tend to have more negative values of n and thus steeper attenuation curve slopes, which once again matches previous observations for quiescent galaxies (Salim et al. 2018;Wiseman et al. 2022).Meldorf et al. (2022) examine the attenuation laws for a larger sample of DES SN Ia host galaxies and recover mostly the same trend for the relation between n and log(M ⋆ /M ⊙ ).However, for the relation between A V and log(M ⋆ /M ⊙ ), they only observe the roughly linear trend associated with star-forming galaxies and do not find a decrease of A V for massive galaxies.As previously stated, this decrease is expected due to the presence of high-mass quiescent galaxies (e.g., Zahid et al. 2013;Salim et al. 2018) and the trend we observe in Fig. 4 follows exactly those expectations.
For the global fit, we find a median R V = 2.762 (using Eq. 1) with a standard deviation of 1.096, while for the local case we find a median R V = 2.612 with a standard deviation of 0.943.In addition, in Table 2 we list the values of R V obtained for both low-and high-mass galaxies, when adopting a standard step of log(M step /M ⊙ ) = 10 for the global case and of log(M step /M ⊙ ) = 9 for the local case.These step locations are based on the median masses for both sets of measurements.We find that, for the global case, low-mass galaxies have a considerably larger median R V , meaning that, for the same amount of dust, they appear to be systematically subject to a lower level of reddening.This once again points to the presence of massive quiescent galaxies in our sample, as for star-forming galaxies the value of R V is expected to increase with the stellar mass (Salim et al. 2018).
Despite having standard deviations of R V ∼ 1, both of the populations exhibit low standard errors, with the difference be- The results for the local case exhibit an opposite tendency, with a higher median R V for the high-mass population.This is probably due to the step location considered being too low, which does not allow the quiescent galaxies to meaningfully impact the high-mass population.

Host galaxy attenuation and SN Ia extinction
As discussed in Section 1, the phenomena governing dust extinction differ from those governing attenuation, as the geometrical effects related to the spatial distribution of dust and light sources cease to be relevant for the former, and only the amount and type of dust are important.It is thus crucial to take notice of the differences between the extinction laws affecting SNe Ia and the attenuation laws affecting their respective host galaxies.
With this in mind, in Fig. 5 we plot the best-fit values for R V and τ V obtained for the attenuation laws for the DES host galaxies.As a term of comparison, we look at the extinction laws obtained by Mandel et al. (2011) for a set of SNe Ia at low redshift.We note that Burns et al. (2014)  It is clear that the extinction for SNe Ia and the attenuation for host galaxies exhibit two different tendencies.This is due to the fact that, for a point source, the effects arising from the dustto-star geometry and observation angle are no longer relevant.
As shown in Section 4.1, the global and the local fit results are consistent, indicating that a 4 kpc aperture is not small enough to eliminate the dust-to-star geometry effects associated with dust attenuation.
For a more direct comparison with our sample, we can look at the light-curve parameters recovered for each of the DES SN Ia by Brout et al. (2019).In particular, looking at the Spearman correlation coefficients between the host galaxy attenuation parameters and the respective SN Ia light-curve parameters, shown in Table 3, it is clear that τ V and n do not meaningfully correlate with the SN properties.In particular, there is no evident correlation with light-curve color, which once again points to the difficulty of using environmental attenuation properties to infer SN dust extinction.These results show that there is not a one-toone relation between SN extinction and host galaxy attenuation, meaning it is difficult to take one as a proxy for the other.

SN Ia cosmology
In this section, we focus on the effects of SN Ia standardization on the Hubble residuals.In particular, we examine whether there is evidence for a dust step in the data and how it compares to the more commonly used mass step.Notes.The thresholds for the two populations are defined as the optimal step value recovered for the Tripp standardization.Additionally, the values for the standard thresholds of log(M step /M ⊙ ) = 10 for the global case and of log(M step /M ⊙ ) = 9 for the local case are also shown.

SN Ia standardization
For each SN Ia, the SALT2 model (Guy, J. et al. 2007;Guy et al. 2010) can be used to fit its light curve, returning values for the light-curve width x 1 , observed color at maximum c, and observed magnitude at maximum m B .Following the standardization formula introduced by Tripp (1998) with the addition of a δ µ bias term, the corrected distance modulus µ for each SN is then given by: where α, β and M are nuisance parameters obtained from a cosmological fit, describing the shape-luminosity and colorluminosity corrections, as well as the absolute magnitude of a fiducial SN Ia with x 1 = 0 and c = 0, respectively, and δ µ bias is a 1D bias correction term, introduced as a function of the redshift and obtained from survey simulations.To account for the mass step, an additional correction term δµ M , dependent on the host galaxy mass M ⋆ , is usually added to the previous expression (e.g., Betoule et al. 2014).
For the cosmological fit, we take advantage of the fact that µ can be expressed as a function of the luminosity distance d L which, for a Friedmann-Robertson-Walker cosmology, only depends on the cosmological parameters and the redshift (Riess et al. 1998).Generally, the model parameters H 0 and Ω m are kept as free parameters, along with α, β and M. By minimizing the Hubble residuals ∆µ = µ − µ model these parameters could therefore be constrained via the fit.In this work, the cosmological parameters are kept fixed to Ω m = 0.3, Ω Λ = 0.7, w = −1 and H 0 = 70km s −1 Mpc −1 , as we are mainly interested in the effects of dust and the color-luminosity correction on the standardizations.
For the standardization fits we use the emcee package to apply a Bayesian fitting procedure, with overall flat priors for the fit parameters.The likelihood function is defined as: where σ 2 is defined as:

0.141
In the previous expression, σ m B , σ x1 and σ c are the uncertainties associated with each of the light-curve fit parameters, σ m B ,c , σ m B ,x 1 and σ x 1 ,c are their covariance terms and σ int is a parameter that accounts for possible intrinsic variations in a SN Ia's luminosity.We fix σ int = 0.107, following the value obtained by González-Gaitán et al. ( 2021) 12 .
It is known that the likelihood defined by Eq. 3 can introduce biases into the recovered fit parameters (Kessler & Scolnic 2017).However, as our main objective is a wholesale comparison between two standardization processes using the same likelihood, these biases should not be of much relevance.This comparison is focused on the improvement of the RMS of the Hubble residuals with respect to the fiducial model.
The cosmological fit parameters recovered following this method for the Tripp standardization are listed in Table .4, with the corresponding values of RMS for the residuals also shown.The light-curve parameters recovered for each of the DES SN Ia by Brout et al. (2019) using the SALT2 model, as well as the corresponding 1D bias corrections, were used for the fits.
A separate "Fixed-Extinction" standardization, in which the extinction contribution to the color-luminosity correction can be fixed using our prior measurements of R V and E(B − V) for the host galaxies from Section 4, is described in Appendix C. For this standardization, an additional free intrinsic color-luminosity parameter β int is introduced.The corresponding fit parameters for this standardization are also displayed in Table .C.1.It however results in much worse fits, as can be seen in the large RMS.

Hubble residuals steps
We now examine possible "steps" in the Hubble residuals for the Tripp standardization, associated with both the local and global properties of the host.We investigate "steps" associated with both the host mass and dust properties.
To look for a step, we divide the Hubble residuals into two populations, according to the properties of the SN host.By varying the threshold values at which the population division is made, we can determine the optimal location for the step.This can be done by optimizing several quantities, namely the step magnitude, the step significance and the ∆RMS with respect to the fiducial Tripp standardization.
For a description with two populations A and B, the step magnitude γ is defined as the difference between the mean values of the Hubble residuals for each of the two populations, while the step significance is defined as the ratio between the step magnitude and the step error σ step , given by: where N A and N B are the number of data points in each of the populations and σ A and σ B are the STDs for each of the populations.
∆RMS is defined as the difference between the standard RMS single obtained when considering the Hubble residuals as a single population and the RMS dual obtained when considering them as two separate populations, each with a separate mean and divided at the step in mass or dust.Thus, the larger the value of ∆RMS, the more the two-population description is favored.

The mass step
When considering a mass step, the two populations are taken to be SNe Ia originating in either low-or high-mass galaxies.Following the procedure detailed in the previous section, we look for a division in mass between two populations that maximizes the step magnitude, its significance and ∆RMS.The optimal step values obtained are shown in Table 5, with the corresponding values of RMS dual also shown.The optimal values are also plotted in Fig. 6.In addition, in Fig. 7 the Hubble residuals for the SNe are plotted, as well as the corresponding optimal steps and the mean values for log(M ⋆ /M ⊙ ) and ∆µ found for each mass bin.
We obtain optimal steps that are not only both significant at > 3σ, but also correspond to the best improvement in RMS.Additionally, both step locations and their respective significance levels roughly match the results obtained by Kelsey et al. (2021) for this data set, pointing to the existence of a mass step for the Tripp standardization.
In Table 2 we present the R V values for both low-and highmass galaxy populations, as defined by the step values recovered for the Tripp standardization.We once again find that, for the global case, low-mass galaxies have a larger median R V , with the difference between the two populations being significant at ∼ 2σ.We also find that, as the local step location shifts to a higher mass, the same tendency starts to be observed in this case.In the next section, we further explore the discrepancy between these populations, analyzing whether it can be described by a dust step.

The dust step
Having established the presence of a mass step, we next look at whether evidence for a step related to the dust content of a SN host can be found in the data for the Tripp standardization.There are some differences to the approach discussed for the mass step.First, while the stellar mass of a galaxy can be described by a single parameter, M ⋆ , dust attenuation is parameterized by two parameters, τ V and n.To account for this fact and ensure we are accurately accounting for a possible dust step, we adopt a twodimensional population division.One-dimensional steps in both τ V and n were attempted, but they were found not to be very significant (∼ 2σ).Second, we consider three different population splits to fully cover the τ V − n parameter space.These are chosen to preserve a two-population division, thus making the mass and dust steps more comparable.
The three splits, schematically represented in Fig. 8, are: one population consisting of SNe Ia originating in low attenuation, high reddening galaxies (Region 1 in Fig. 8) and another consisting of the remaining SNe, defined as the First Split; one population consisting of SNe Ia originating in low attenuation, low reddening galaxies (Region 2 in Fig. 8) and another consisting of the remaining SNe, defined as the Second Split; one population consisting of SNe Ia originating in high attenuation, low reddening galaxies (Region 3 in Fig. 8) and another consisting of the remaining SNe, defined as the Third Split.There is also Notes.For the Tripp+δµ D standardization, the step location and magnitude recovered for the optimal dust step found in Section 5.2.2 are used.
∆RMS is defined as the difference between the RMS single obtained when considering the Hubble residuals as a single population and the RMS dual obtained when considering them as two separate populations, each with a separate mean and divided at the step.a fourth possible split that is not considered, as the region with high attenuation and high reddening galaxies (Region 4 in Fig. 8) is very sparsely populated.
The optimal step values for the First, Second and Third Splits are shown in Table 6.The Second Split results in the best significance and ∆RMS values, which means it is by far the most favorable for a dust step.We can further analyze this particular case by taking a closer look at the variation of the step magnitude, significance and ∆RMS, plotted in Fig. 9.We can see that, not only are the optimal step locations very narrow, but they also exhibit clear delimitations for both dust parameters, suggesting we are indeed encountering a two-dimensional step.This is particularly apparent when using global data.In the local case, there appear to be two possible step locations.We take the optimal step to be the one that maximizes ∆RMS.The significance, magnitude and ∆RMS values obtained for the optimal global dust step are very close to those recovered for the optimal global mass step, even though the dust step ends up being slightly stronger.In particular, the dust steps are significant at > 4σ, which matches the results obtained by Meldorf et al. (2022) for DES-SN host galaxies.

Mass and dust step comparison
For a more rigorous SN standardization, the steps recovered in the previous sections can be introduced into the usual Tripp standardization.In particular, we do this to ascertain whether correcting for one of the steps eliminates the other one.In such cases, the mass step takes the form:  Notes.For the Tripp+δµ M standardization, the step location and magnitude recovered for the optimal mass step found in Section 5.2.1 are used.
∆RMS is defined as the difference between the RMS single obtained when considering the Hubble residuals as a single population and the RMS dual obtained when considering them as two separate populations, each with a separate mean and divided at the step.The First, Second and Third splits refer to which of the quadrants depicted in Fig. 8 is isolated from the remaining SN population (1,2 and 3, respectively).
where M step is the optimal step location and γ M the corresponding step magnitude recovered in Section 5.2.1 and recorded in Table 5.For the Second Split, the dust step takes the form: where τ V step and n step are the optimal step locations and γ D the corresponding step magnitude recovered in Section 5.2.2 and recorded in Table 6.It is clear from Fig. 10 that the two mass and dust populations do not exactly align.In fact, higher-mass galaxies tend to concentrate on the inside of the τ V − n curve, while less massive galaxies concentrate on the outside.This is particularly apparent in the local case.This discrepancy points to the fact that the recovered dust step is not merely analogous to the mass step, as might initially be presumed, given its magnitude, significance and ∆RMS.
For a more detailed look at this discrepancy, we can add a dust step in the form of Eq. 7 to the Hubble residuals for the Tripp standardization.To do so, we fix the values of γ D , τ V step and n step to those previously obtained in Section 5.2.2 for the optimal dust step, as recorded in Table 6.Following a similar procedure as the one applied in Section 5.2.1, we once again recover the optimal mass step, shown in Table 5.It should be noted that, when looking for this new mass step, the baseline for comparison is not the case with no step, but rather the case already corrected for the dust step.
For the local case, the optimal mass-step location remains unchanged.For the global case, however, there is a big change in the maximum significance step location, which happens because the previous optimal step significance drops slightly below the level of the new maximum.For the sake of consistency, we keep the same step location.Even so, the new maximum significance step, log(M ⋆ /M ⊙ ) = 10.51, has a significance of 3.20.It is clear that the introduction of the dust-step correction lowers both the significance, magnitude and ∆RMS of the recovered mass steps.This drop in ∆RMS should not be taken to mean that the inclusion of both steps results in an overall worse RMS when compared with a single step case.Rather, it signifies that the addition of a mass step to a set of Hubble residuals already corrected with a dust step results in a lower improvement in RMS than the one seen for the same mass step applied to Hubble residuals with no prior step correction.In this sense, the initial dust-step correction appears to also be partially correcting for the mass effect.However, it cannot be said that the mass step has been removed in either case.Thus, while offering some improvement, the addition of this dust-step term to the standardization is not capable of fully eliminating the mass step by itself.
In the same way, we can look at the recovered dust step in the case of a Tripp+mass-step standardization, in which the mass step recovered in Section 5.2.1 is used.The values relating to these steps are shown in Table 6.Once again, it is clear that the structure of the original steps is preserved, while the significance, magnitude and ∆RMS are again found to be lower.Thus, the addition of a mass-step term in the standardization is also not capable of fully eliminating the dust step by itself.
Thus, it seems reasonable to conclude that the mass and dust steps are two overlapping phenomena that nevertheless do not fully coincide.This points to the fact that the origin of the mass step may not be entirely driven by host dust properties.A further exploration of the relation between the two steps is offered in Appendix D, where both µ M and µ D are treated as free parameters in the cosmological fit.This confirms that the two steps are not analogous, as both are significantly different from zero for the best-fit cosmological model.
Looking at a bigger sample of DES host galaxies that includes SN Ia photometry, Meldorf et al. (2022) find that the inclusion of a dust-step correction is sufficient to correct the mass step.The discrepancy between this and our results is most likely due to the different mass to dust relations recovered from the host galaxy fits, as discussed in Section 4.3.

SN color and Hubble residuals steps
Although we do not find any direct relation between SN color and the galaxy dust parameters (Section 4.4), we explore the relation between SN color and the Hubble residuals, similarly to what is done in (e.g., Brout & Scolnic 2021).The SNe were divided into two populations: Blue SNe (c < 0) and Red SNe (c > 0).For an easier comparison with the previous results, the optimal locations found in Sections 5.2.1 and 5.2.2 were preserved.These results are shown in Tables 7 and 8, for the mass and dust steps, respectively.
It is clear from the significance, magnitude and ∆RMS values that both steps are much more pronounced for the red SNe, which matches the results obtained by Kelsey et al. (2021) and Brout & Scolnic (2021).This difference between the two populations points once again to a relation between dust reddening and the Hubble residuals steps, assuming most of the color variation for SNe stems from dust.However, it has also been suggested that there is possible indirect evidence of environmental variation of intrinsic SN color, particularly related to ejecta velocity (e.g., Pan 2020; Wang et al. 2013;Cartier et al. 2011), which could also be playing a partial role in this case.
Finally, similarly to what was done in Sec.5.2.3, we examine what happens to the mass step if the dust step is corrected for and vice versa.Depending on the population in question, we use either the steps recovered for the blue or the red SNe.The results remain consistent with what we have shown before, with neither of the steps being able to fully account for the other, regardless of the color of the SNe considered.

Discussion
Several key results obtained in this work appear to support the link between environmental dust properties and the mass step.
First, the mass-τ V and mass-n relations, shown in Fig. 4, make it clear that different galaxy mass populations are subject to different dust contents and dust laws.This fact is most clearly expressed in the results displayed in Table 2, where it is shown that the R V values for low-and high-mass galaxies exhibit significant differences, which agrees with previous results (e.g., Salim et al. 2018;Brout & Scolnic 2021;González-Gaitán et al. 2021).
Second, in Section 5.2.3 we show that adding a dust step to a set of mass-step corrected Hubble residuals results in lower significances, magnitudes and ∆RMS than the ones obtained for the same step when no initial correction is present, and vice versa.This shows that the effects in mass and dust are at least partially related and that correcting for one also partly corrects for the other.
Third, as shown in Section 5.2.4,both the mass and dust steps are found to have much higher values of significance, magnitude and ∆RMS for red SNe (c > 0), when compared to blue SNe (c < 0).Presuming that red SNe are also the most extinguished, we have yet another link between dust and the mass step.
However, the color variations among SNe are not only a product of dust extinction, but can also originate from intrinsic  differences between the SNe (e.g., Conley et al. 2007).For this reason, the fact that the mass step is found to be more prominent for red SNe might also indicate a relation between this effect and the SN intrinsic color (Pan 2020;Wang et al. 2013).
Another result that points against dust as the sole driver for the mass step is the overall lack of correlation found between SN color and the corresponding host galaxy attenuation parameters, as recorded in Table 3.This again suggests that the differences in the mass step between different color SNe might actually be due to intrinsic factors.Another possibility is that global and local (4kpc) environmental dust properties might be very poor proxies of the real dust extinction at the SN site.
This last point is further strengthened by the fact that, as shown in Fig. 5, extinction and attenuation display very different behaviors, pointing to different physical effects.This also explains why the proposed "Fixed-Extinction" Standardization, discussed in Appendix C, gives such poor results.
Furthermore, while it is true that there appears to be an overlap between the two steps, we show that accounting for one of them does not fully remove the other, suggesting that there might be other factors at play, not directly related to dust.
In general, our results seem to point to the fact that the mass step might not have a fully dust-related origin, despite the evidence for at least a partial link between the two.As stated be-fore, it is possible that a partial cause of the observed effect lies instead in intrinsic differences between SNe Ia that come about due to different progenitors.The fact that the stellar age seems to be quite important in the relations between the attenuation law and the mass (see Fig. 4) shows that the situation is more complex and agrees with other works that show the importance of age in the mass step (e.g., Rigault, M. et al. 2020;Rose et al. 2022;Wiseman et al. 2022).In this sense, both the environmental mass and dust might be acting as imperfect proxies for the different types of progenitors, explaining the behavior we observe.
Finally, while dust might indeed be the origin of the mass step, a two-dimensional environmental dust step might be an overly simplistic way to look at the relation between these two quantities.One might need either direct extinction data for each supernova or a better way to isolate the different mass/dust populations.

Conclusions
In this work, we have explored ways to better probe the dust contents of SN Ia host galaxies as well as how to better make use of this information for SN standardization.Our major findings are: The corresponding value of log(M ⋆ /M ⊙ ) is plotted on a gradient from blue (lower mass) to red (higher mass).The optimal step for the Second Split is also shown in purple.Notes.The top rows show the results pertaining to the Tripp standardization, while the bottom rows show the results pertaining to the Tripp+δµ D .
In both cases, the same step locations obtained in Section 5.2.1 are used.∆RMS is defined as the difference between the RMS single obtained when considering the Hubble residuals as a single population and the RMS dual obtained when considering them as two separate populations, each with a separate mean and divided at the step.
-We have shown through simulations that DECam griz photometry is mostly sufficient to recover dust parameters for simulated SN Ia host galaxies, while discrepancies and degeneracies among other parameters do not significantly affect inferred dust properties.We found that in using both global and local DECam griz photometry, we can recover dust properties for host galaxies that are consistent with literature predictions from both simulations and observations.We find a relation for the dust attenuation slope with the dust optical depth, both locally and globally, that is best explained with varying star-to-dust geometry with a galaxy orientation.Most importantly, we show that dust properties vary greatly across different galaxies, meaning a universal value of R V = 3.1 or a universal SN Ia β correction are assumptions that are too simplistic.-The relation between the attenuation parameters is found to be very different from the extinction relations obtained di-rectly for SNe Ia, making the comparison between the two somewhat difficult.This difference is mainly due to the effects related to star-to-dust geometry that become relevant when dealing with extended objects, such as SN host galaxies.
-We find two different relations between a galaxy's stellar mass and its dust properties, depending on the age of the corresponding stellar population.For younger galaxies, there is a steady increase of τ V and n with log(M ⋆ /M ⊙ ), while older galaxies tend to exhibit both larger masses and smaller optical depths.A similar behavior for the two age populations can be seen when considering the local mass, pointing to an analogous relation between the results observed for the global stellar mass and the local stellar density.This population mix results in the observation that low-mass galaxies tend to have larger values of R V , meaning that, for the same amount of dust, these galaxies are systematically subject to a lower level of reddening.For both mass populations, the scatter in R V is ∼ 1.All of these results match previous literature observations.-We conclude that an alternative SN Ia standardization, incorporating both the R V and E(B − V) values obtained for the respective host galaxies to approximate and fix SN extinction, results in a worse cosmological fit.This is very likely a consequence of the differences between host galaxy attenuation and SN Ia extinction, meaning the values for β R V and E(B − V) for the SNe Ia cannot be correctly determined.The lack of correlation between the host attenuation parameters and the respective SN light-curve parameters also points to this fact.This remains true when carrying out local studies (4 kpc), showing that even smaller apertures are needed to avoid the geometrical effects of dust attenuation.-We show there is evidence for a two-dimensional dust step, which, for the global case, results in roughly the same significance, magnitude, and residual RMS recovered for the mass step, hinting at a relation between the two.In particular, this dust step is significant at > 4σ.Furthermore, we show that these two steps are not completely analogous and that accounting for one of them in the SN standardization does not necessarily fully eliminate the other.-We find that both the mass and dust steps are much more pronounced for red SNe Ia (c > 0), which once again suggests a link between dust reddening and the Hubble residual steps.
Thus, we conclude that the need for the mass-step correction cannot be completely eliminated using only host galaxy dust data.This is true for both an alternative SN standardization, in which the color-luminosity relation is constrained with host attenuation, and for a SN standardization with a dust-step correction.constrain the color-luminosity relation in the respective "Fixed-Extinction" standardization.The σ 2 in Eq. 3 is thus given by: where σ m B , σ x1 and σ c are the uncertainties associated with each of the light-curve fit parameters, σ m B ,c , σ m B ,x 1 , σ x 1 ,c are their covariance terms, σ A B and σ E(B−V) are the uncertainties associated with the attenuation law parameters for the host galaxies and σ int is a parameter which accounts for possible intrinsic variations in a SN Ia's luminosity.We fix σ int = 0.107, following the value obtained by González-Gaitán et al. (2021).
The cosmological fit parameters obtained for the "Fixed-Extinction" standardization are shown in Table .C.1.This standardization yields RMS values of 0.515 and 0.575 in the global and the local cases, respectively.As previously hinted, both the global and the local measurements produce consistent results.However, judging from the RMS values, the "Fixed-Extinction" standardization appears to be a much worse fit of the data, invalidating it as a replacement for the Tripp standardization, which has a RMS value of 0.141.This is probably due to the fact that, according to Fig. 5, host attenuation laws do not accurately describe SN extinction, meaning it is likely that the problem resides in the values used to approximate SN extinction.
A possible source of error may come from the bias in τ V observed in Section 3, which might be skewing the standardization fit results.We can correct for this effect by subtracting the bias level from each τ V value.For this, the largest bias value found was used, which was 0.263.For values below the bias threshold, we simply keep the original value.The results obtained from a "Fixed-Extinction" biased standardization are also listed in Table.C.1.The τ V bias improves the standardization, with new RMS values of 0.420 and 0.372, for the global and the local cases, respectively.Even so, the bias corrected standardization still falls short of the Tripp values.
On the whole, the proposed "Fixed-Extinction" SN Ia standardization does not provide better fits than standard Tripp standardization.This is most likely a consequence of using host attenuation to approximate SN Ia extinction, as already seen in Section 4.4.However, it is possible that, using extinction data directly obtained for each SN and not its environment, this method might indeed prove to significantly improve cosmological fits and possibly reduce or completely eliminate the mass step.

Appendix D: Cosmological fits with steps
In this appendix, a deeper look at the relation between the mass and dust steps is presented.In particular, for a slightly more robust analysis of both these steps, we can directly introduce δµ M and δµ D terms in the cosmological fit.We promote the respective step magnitudes γ M and γ D to free parameters, fixing the step locations to the values previously recovered for the Tripp standardization in Section 5.2.
As mentioned in Section 5.1, the likelihood defined by Eq. 3 can introduce biases into the recovered fit parameters.While this fact is not entirely significant when comparing two standardizations, it becomes important when dealing with a more rigorous cosmological fit.For this reason, when promoting either the mass or dust steps to free fit parameters, it is helpful to slightly modify the likelihood used (e.g., González-Gaitán et al. 2021): In order to study the various steps, we look at three possible standardizations: a Tripp+mass-step standardization, in which γ M is promoted to a free fit parameter; a Tripp+dust-step standardization, in which γ D is promoted to a free fit parameter; and a Tripp+mass-step+dust-step standardization, in which both γ M and γ D are promoted to free fit parameters.It is also useful to redo the standard Tripp standardization to establish a term of comparison for the new likelihood.

0.372
Notes.A likelihood defined by Eq. 3 was used for the fits.
The recovered fit parameters, as well as the corresponding RMS values, are listed in Table D.1.As expected, the addition of separate mass and dust steps results in lower RMS values, for both the global and the local cases.The dust-step value recovered for the global case is compatible with zero, which might be due to a problem with the fit convergence.The simultaneous inclusion of mass and dust steps produces the lowest values of RMS for the respective cosmological fits.In addition, when both steps are simultaneously included, we recover both noticeable mass and dust steps, with the caveat that the global dust step is once again compatible with zero.In general, the step magnitudes in this case are overall lower than those obtained in the single step cases.As such, we once again find evidence that, although there is some level of overlap, the two steps are not completely analogous and interchangeable.

Fig. 1 .Fig. 2 .
Fig. 1.Best-fit values of n as a function of τ V for the fitted DES galaxies with DECam global (Left panel) and local (Right panel) griz photometry.Results for the different galaxies are shown in gray.A Gaussian Process Regression is shown in solid blue, with dashed blue lines defining a 68% credible interval.The mean relation found for similar simulated galaxies byChevallard et al. (2013) is shown in dashed orange, with ±25% error margins shown in dotted orange.

Fig. 3 .
Fig. 3. Schematic representation of dust attenuation for different galaxy orientations: Head-on (Left panel) and Edge-on (Right Panel).

Fig. 4 .
Fig. 4. Best-fit results for the dust parameters as a function of log(M ⋆ /M ⊙ ) for the fitted DES galaxies with DECam global (Left panels) and local (Right panels)griz photometry.Plots for τ V are shown in the top row, while results for n are shown in the bottom row.The best-fit for the age of each galaxy is plotted on a gradient from blue (younger) to red (older).
find a similar relation to Mandel, whereas recently Rose et al. (2022) find a flat distribution of R V for various extinction amounts of SNe Ia.

Fig. 5 .
Fig.5.Best-fit values of R V as a function of τ V for the fitted DES galaxies with DECam Global (Left panel) and Local (Right panel) griz photometry.Results for the different galaxies are shown in gray.A Gaussian Process Regression is shown in solid blue, with dashed blue lines defining a 68% credible interval.The correlation found byMandel et al. (2011) between these two quantities for the extinction laws of a set of nearby SNe Ia is shown in orange.

Fig. 6 .
Fig. 6.Variation of the mass step with step location for the Tripp standardization, for the global (Left panels) and the local (Right panels) cases.The top panel shows the significance of the step in σ, the middle panel shows the magnitude of the step in solid black, with the gray region showing the uncertainty, and the lower panel shows ∆RMS.In these three panels, the location of the step of maximum significance is shown in blue and the median mass of the sample is shown in orange.Inspired byKelsey et al. (2021).

Fig. 7 .
Fig.7.Hubble residuals ∆µ for the DES SNe Ia as a function of the stellar mass for the respective host galaxies, plotted in gray.The optimal step locations, as recorded in Table.5 are shown in dashed blue.The blue dots show the mean values of log(M ⋆ /M ⊙ ) and ∆µ in each bin, with the error in the x-axis given by the standard deviation.

Fig. 8 .
Fig. 8. Schematic representations of the regions singled out by each of the proposed dust-step splits in the τ V − n parameter space.The best-fit values of n as a function of τ V for the fitted DES galaxies with global griz photometry are shown for reference.The dashed lines only represent possible division points for the populations, with the actual optimal values recorded in in Table.6.

Fig. 9 .
Fig. 9. Variation of the second split dust step with the step location (position in τ V − n grid) for the Tripp standardization.Results from global measurements are shown on the left, while those from local measurements are shown on the right.The top panel shows the significance of the step in σ, the middle panel shows the magnitude of the step and the lower panel shows ∆RMS.This last quantity is defined as the difference between the RMS single obtained for a single population mean description and the RMS dual obtained for a dual population mean description with a step.

Fig. 10 .
Fig. 10.Best-fit values of n as a function of τ V for the fitted DES galaxies with DECam global (Left panel) and local (Right panel) griz photometry.The corresponding value of log(M ⋆ /M ⊙ ) is plotted on a gradient from blue (lower mass) to red (higher mass).The optimal step for the Second Split is also shown in purple.
Fig. A.1.Comparison between the stellar mass values obtained byKelsey et al. (2021) for the host galaxies of the DES-SN sample and the masses recovered using a five-parameter sampling with no initial minimization.In both cases, fits were performed using global (Left panel) and local (Right panel) DECam griz photometry.The results are shown in blue.A dashed line with y = x is also shown.
Fig. B.1.Residuals (∆ = Fit-Simulation) for the best-fit values for the default (griz) test galaxy fits as a function of the original simulation parameters.From top to bottom and left to right, we have: log(Z ⋆ /Z ⊙ ), t age , τ V and n.Results for the different galaxies are shown in gray.Binned means for each parameter are shown in orange, with error bars given by the standard deviation in each bin.
Fig. B.3.Residuals (∆ = Fit-Simulation) for the best-fit values for the Cardelli test galaxy fits as a function of the original simulation parameters.From top to bottom and left to right we have: log(Z ⋆ /Z ⊙ ), t age , τ V and n.Results for the different galaxies are shown in gray.Binned means for each parameter are shown in orange, with error bars given by the standard deviation in each bin.Table C.1.Global and Local fit parameters for the "Fixed-Extinction" and "Fixed-Extinction" + τ V bias standardizations and respective RMS Standardization α β int M RMS Fixed-Extinction (Global) 0.432 +0.010 −0.011 4.350 +0.010 −0.010 −19.106 +0.011 −0.009 0.515

Table 1 .
Precision and accuracy metrics for the dust parameters in the mock galaxy fits, for each of the test sets.

Table 2 .
Median global and local R V values, standard deviations (S T D) and standard errors (S E) for low-and high-mass galaxies.

Table 3 .
Spearman correlation coefficients between the global and the local host galaxy attenuation parameters and the respective SN Ia lightcurve parameters.

Table 4 .
Fit parameters for the Tripp standardization, using a likelihood defined by Eq. 3. The corresponding value of RMS is also shown.

Table 8 .
Dust-step significance (in σ), magnitude (γ D ), ∆RMS and RMS dual for different color populations.Notes.The top rows show the results pertaining to the Tripp standardization, while the bottom rows show the results pertaining to the Tripp+δµ M .In both cases, the same step locations obtained in Section 5.2.2 are used.∆RMS is defined as the difference between the RMS single obtained when considering the Hubble residuals as a single population and the RMS dual obtained when considering them as two separate populations, each with a separate mean and divided at the step.