Search

We present the first homogeneous release of several thousand Type Ia supernovae (SNe Ia), all having spectroscopic classification, and spectroscopic redshifts for half the sample. This release, named the "DR2", contains 3628 nearby (z < 0.3) SNe Ia discovered, followed and classified by the Zwicky Transient Facility survey between March 2018 and December 2020. Of these, 3000 have good-to-excellent sampling and 2667 pass standard cosmology light-curve quality cuts. This release is thus the largest SN Ia release to date, increasing by an order of magnitude the number of well characterized low-redshift objects. With the "DR2", we also provide a volume-limited (z < 0.06) sample of nearly a thousand SNe Ia. With such a large, homogeneous and well controlled dataset, we are studying key current questions on SN cosmology, such as the linearity SNe Ia standardization, the SN and host dependencies, the diversity of the SN Ia population, and the accuracy of the current light-curve modeling. These, and more, are studied in detail in a series of articles associated with this release. Alongside the SN Ia parameters, we publish our force-photometry gri-band light curves, 5138 spectra, local and global host properties, observing logs, and a python tool to ease use and access of these data. The photometric accuracy of the "DR2" is not yet suited for cosmological parameter inference, which will follow as "DR2.5" release. We nonetheless demonstrate that the multi-thousand SN Ia Hubble Diagram has a typical 0.15 mag scatter., Comment: ZTF SN Ia DR2 release paper. Submitted to A&A (ZTF DR2 Special Issue). Already 1 response to referee

$ $Low density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence $\kappa$. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 data set, covering approximately 4,200 deg$^2$ of the sky. We use two distinct void-finding algorithms: a 2D void-finder which operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the MICE N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure $A_\kappa$, the amplitude of the observed lensing signal relative to the simulation template, obtaining $A_\kappa = 1.03 \pm 0.22$ ($4.6\sigma$ significance) for Voxel and $A_\kappa = 1.02 \pm 0.17$ ($5.9\sigma$ significance) for 2D voids, both consistent with $\Lambda$CDM expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure $A_\kappa = 0.87 \pm 0.15$ ($5.9\sigma$ significance). The leading source of noise in our measurements is Planck noise, implying that future data from the Atacama Cosmology Telescope (ACT), South Pole Telescope (SPT) and CMB-S4 will increase sensitivity and allow for more precise measurements., Comment: 16 pages, 9 figures. Accepted for publication in MNRAS

Dust associated with various stellar sources in galaxies at all cosmic epochs remains a controversial topic, particularly whether supernovae (SNe) play an important role in dust production. We report evidence of dust formation in the cold, dense shell behind the ejecta-circumstellar medium (CSM) interaction in the Type Ia-CSM SN 2018evt three years after the explosion, characterized by a rise in the mid-infrared (MIR) emission accompanied by an accelerated decline in the optical radiation of the SN. Such a dust-formation picture is also corroborated by the concurrent evolution of the profiles of the Ha emission line. Our model suggests enhanced CSM dust concentration at increasing distances from the SN as compared to what can be expected from the density profile of the mass loss from a steady stellar wind. By the time of the last MIR observations at day +1041, a total amount of 1.2+-0.2x10^{-2} Msun of new dust has been formed by SN 2018evt, making SN 2018evt one of the most prolific dust factories among SNe with evidence of dust formation. The unprecedented witness of the intense production procedure of dust may shed light on the perceptions of dust formation in cosmic history., Comment: Accepted by Nature Astronomy, 6 main figures, 7 extended figures, and 2 extended tables

An open question in SN Ia research is where the boundary lies between 'normal' Type Ia supernovae (SNe Ia) that are used in cosmological measurements and those that sit off the Phillips relation. We present the spectroscopic modelling of one such '86G-like' transitional SN Ia, SN 2021rhu, that has recently been employed as a local Hubble Constant calibrator using a tip of the red-giant branch measurement. We detail its modelling from -12 d until maximum brightness using the radiative-transfer spectral-synthesis code tardis. We base our modelling on literature delayed-detonation and deflagration models of Chandrasekhar mass white dwarfs, as well as the double-detonation models of sub-Chandrasekhar mass white dwarfs. We present a new method for 'projecting' abundance profiles to different density profiles for ease of computation. Due to the small velocity extent and low outer densities of the W7 profile, we find it inadequate to reproduce the evolution of SN 2021rhu as it fails to match the high-velocity calcium components. The host extinction of SN 2021rhu is uncertain but we use modelling with and without an extinction correction to set lower and upper limits on the abundances of individual species. Comparing these limits to literature models we conclude that the spectral evolution of SN 2021rhu is also incompatible with double-detonation scenarios, lying more in line with those resulting from the delayed detonation mechanism (although there are some discrepancies, in particular a larger titanium abundance in SN 2021rhu compared to the literature). This suggests that SN 2021rhu is likely a lower luminosity, and hence lower temperature, version of a normal SN Ia., Comment: 25 pages, 22 figures, accepted for publication in MNRAS

Measurements of the growth rate of structures at $z < 0.1$ with peculiar velocity surveys have the potential of testing the validity of general relativity on cosmic scales. In this work, we present growth-rate measurements from realistic simulated sets of type-Ia supernovae (SNe Ia) from the Zwicky Transient Facility (ZTF). We describe our simulation methodology, the light-curve fitting and peculiar velocity estimation. Using the maximum likelihood method, we derive constraints on $f\sigma_8$ using only ZTF SN Ia peculiar velocities. We carefully tested the method and we quantified biases due to selection effects (photometric detection, spectroscopic follow-up for typing) on several independent realizations. We simulated the equivalent of 6 years of ZTF data, and considering an unbiased spectroscopically typed sample at $z < 0.06$, we obtained unbiased estimates of $f\sigma_8$ with an average uncertainty of 19% precision. We also investigated the information gain in applying bias correction methods. Our results validate our framework which can be used on real ZTF data., Comment: 20 pages, 13 figures

Type Ia supernovae (SNe Ia) are thermonuclear explosions of degenerate white dwarf (WD) stars destabilized by mass accretion from a companion star, but the nature of their progenitors remains poorly understood. A way to discriminate between progenitor systems is through radio observations; a non-degenerate companion star is expected to lose material through winds or binary interaction prior to explosion, and the SN ejecta crashing into this nearby circumstellar material (CSM) should result in radio synchrotron emission. However, despite extensive efforts, no SN Ia has ever been detected at radio wavelengths, which suggests a clean environment and a companion star that is itself a degenerate WD star. Here we report on the study of SN 2020eyj, a SN Ia showing helium-rich CSM, as revealed by its spectral features, infrared emission and, for the first time in a SN Ia, a radio counterpart. Based on our modeling, we conclude the CSM likely originates from a single-degenerate (SD) binary system where a WD accretes material from a helium donor star, an often hypothesized formation channel for SNe Ia. We describe how comprehensive radio follow-up of SN 2020eyj-like SNe Ia can improve the constraints on their progenitor systems., Comment: 62 pages, 10 figures, 4 tables

Cosmological analyses with type Ia supernovae (SNe Ia) often assume a single empirical relation between color and luminosity ($\beta$) and do not account for varying host-galaxy dust properties. However, from studies of dust in large samples of galaxies, it is known that dust attenuation can vary significantly. Here we take advantage of state-of-the-art modeling of galaxy properties to characterize dust parameters (dust attenuation $A_V$, and a parameter describing the dust law slope $R_V$) for the Dark Energy Survey (DES) SN Ia host galaxies using the publicly available \texttt{BAGPIPES} code. Utilizing optical and infrared data of the hosts alone, we find three key aspects of host dust that impact SN Ia cosmology: 1) there exists a large range ($\sim1-6$) of host $R_V$ 2) high stellar mass hosts have $R_V$ on average $\sim0.7$ lower than that of low-mass hosts 3) there is a significant ($>3\sigma$) correlation between the Hubble diagram residuals of red SNe Ia that when corrected for reduces scatter by $\sim13\%$ and the significance of the ``mass step'' to $\sim1\sigma$. These represent independent confirmations of recent predictions based on dust that attempted to explain the puzzling ``mass step'' and intrinsic scatter ($\sigma_{\rm int}$) in SN Ia analyses. We also find that red-sequence galaxies have both lower and more peaked dust law slope distributions on average in comparison to non red-sequence galaxies. We find that the SN Ia $\beta$ and $\sigma_{\rm int}$ both differ by $>3\sigma$ when determined separately for red-sequence galaxy and all other galaxy hosts. The agreement between fitted host-$R_V$ and SN Ia $\beta$ \& $\sigma_{\rm int}$ suggests that host dust properties play a major role in SN Ia color-luminosity standardization and supports the claim that SN Ia intrinsic scatter is driven by $R_V$ variation., Comment: 22 pages. Submitted to MNRAS

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$\sigma$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts., Comment: 10 pages, 4 figures. Accepted to ApJ Letters

The current Cepheid-calibrated distance ladder measurement of $H_0$ is reported to be in tension with the values inferred from the cosmic microwave background (CMB), assuming standard cosmology. However, some tip of the red giant branch (TRGB) estimates report $H_0$ in better agreement with the CMB. Hence, it is critical to reduce systematic uncertainties in local measurements to understand the Hubble tension. In this paper, we propose a uniform distance ladder between the second and third rungs, combining SNe~Ia observed by the Zwicky Transient Facility (ZTF) with a TRGB calibration of their absolute luminosity. A large, volume-limited sample of both calibrator and Hubble flow SNe~Ia from the \emph{same} survey minimizes two of the largest sources of systematics: host-galaxy bias and non-uniform photometric calibration. We present results from a pilot study using existing TRGB distance to the host galaxy of ZTF SN~Ia SN 2021rhu (aka ZTF21abiuvdk) in NGC7814. Combining the ZTF calibrator with a volume-limited sample from the first data release of ZTF Hubble flow SNe~Ia, we infer $H_0 = 76.94 \pm 6.4\, {\rm km}\,{\rm s^{-1}}\,{\rm Mpc^{-1}}$, an $8.3 \%$ measurement. The error budget is dominated by the single object calibrating the SN~Ia luminosity in this pilot study. However, the ZTF sample includes already five other SNe~Ia within $\sim$ 20 Mpc for which TRGB distances can be obtained with HST. Finally, we present the prospects of building this distance ladder out to 80 Mpc with JWST observations of more than one hundred ZTF SNe~Ia., Comment: 10 pages (incl. references), 4 Figures, 2 Tables, accepted to ApJ

The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events., Comment: Accepted for publication in ApJS. arXiv admin note: text overlap with arXiv:1812.04051

Forthcoming time-domain surveys, such as the Rubin Observatory Legacy Survey of Space and Time, will vastly increase samples of supernovae (SNe) and other optical transients, requiring new data-driven techniques to analyse their photometric light curves. Here, we present the "Python for Intelligent Supernova-COsmology Light-curve Analysis" (PISCOLA), an open source data-driven light-curve fitter using Gaussian Processes that can estimate rest-frame light curves of transients without the need for an underlying light-curve template. We test PISCOLA on large-scale simulations of type Ia SNe (SNe Ia) to validate its performance, and show it successfully retrieves rest-frame peak magnitudes for average survey cadences of up to 7 days. We also compare to the existing SN Ia light-curve fitter SALT2 on real data, and find only small (but significant) disagreements for different light-curve parameters. As a proof-of-concept of an application of PISCOLA, we decomposed and analysed the PISCOLA rest-frame light-curves of SNe Ia from the Pantheon SN Ia sample with Non-Negative Matrix Factorization. Our new parametrization provides a similar performance to existing light-curve fitters such as SALT2. We further derived a SN Ia colour law from PISCOLA fits over $\sim$3500 to 7000\AA, and find agreement with the SALT2 colour law and with reddening laws with total-to-selective extinction ratio $R_V \lesssim 3.1$., Comment: Accepted in MNRAS

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond., Comment: 32 pages, 18 figures, accepted for publication by ApJ

We present the statistical methods that have been developed to analyse the OzDES reverberation mapping sample. To perform this statistical analysis we have created a suite of customisable simulations that mimic the characteristics of each source in the OzDES sample. These characteristics include: the variability in the photometric and spectroscopic lightcurves, the measurement uncertainties, and the observational cadence. By simulating the sources in the OzDES sample that contain the CIV emission line, we developed a set of criteria that rank the reliability of a recovered time lag depending on the agreement between different recovery methods, the magnitude of the uncertainties, and the rate at which false positives were found in the simulations. These criteria were applied to simulated light curves and these results used to estimate the quality of the resulting Radius-Luminosity relation.We grade the results using three quality levels (gold, silver and bronze). The input slope of the R-L relation was recovered within $1\sigma$ for each of the three quality samples, with the gold standard having the lowest dispersion with a recovered a R-L relation slope of $0.454\pm 0.016$ with an input slope of 0.47. Future work will apply these methods to the entire OzDES sample of 771 AGN., Comment: 16 pages, 15 figures, accepted for publication in NMRAS

We study a phenomenological class of models where dark matter converts to dark radiation in the low redshift epoch. This class of models, dubbed DMDR, characterizes the evolution of comoving dark matter density with two extra parameters, and may be able to help alleviate the observed discrepancies between early- and late-time probes of the universe. We investigate how the conversion affects key cosmological observables such as the CMB temperature and matter power spectra. Combining 3x2pt data from Year 1 of the Dark Energy Survey, {\it Planck}-2018 CMB temperature and polarization data, supernovae (SN) Type Ia data from Pantheon, and baryon acoustic oscillation (BAO) data from BOSS DR12, MGS and 6dFGS, we place new constraints on the amount of dark matter that has converted to dark radiation and the rate of this conversion. The fraction of the dark matter that has converted since the beginning of the universe in units of the current amount of dark matter, $\zeta$, is constrained at 68\% confidence level to be $<0.32$ for DES-Y1 3x2pt data, $<0.030$ for CMB+SN+BAO data, and $<0.037$ for the combined dataset. The probability that the DES and CMB+SN+BAO datasets are concordant increases from 4\% for the $\Lambda$CDM model to 8\% (less tension) for DMDR. The tension in $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ between DES-Y1 3x2pt and CMB+SN+BAO is slightly reduced from $2.3\sigma$ to $1.9\sigma$. We find no reduction in the Hubble tension when the combined data is compared to distance-ladder measurements in the DMDR model. The maximum-posterior goodness-of-fit statistics of DMDR and $\Lambda$CDM model are comparable, indicating no preference for the DMDR cosmology over $\Lambda$CDM., Comment: 23 pages, 11 figures. DES extensions group. Accepted by PRD

Binary supermassive black holes (BSBHs) are expected to be a generic byproduct from hierarchical galaxy formation. The final coalescence of BSBHs is thought to be the loudest gravitational wave (GW) siren, yet no confirmed BSBH is known in the GW-dominated regime. While periodic quasars have been proposed as BSBH candidates, the physical origin of the periodicity has been largely uncertain. Here we report discovery of a periodicity (P=1607$\pm$7 days) at 99.95% significance (with a global p-value of ~$10^{-3}$ accounting for the look elsewhere effect) in the optical light curves of a redshift 1.53 quasar, SDSS J025214.67-002813.7. Combining archival Sloan Digital Sky Survey data with new, sensitive imaging from the Dark Energy Survey, the total ~20-yr time baseline spans ~4.6 cycles of the observed 4.4-yr (restframe 1.7-yr) periodicity. The light curves are best fit by a bursty model predicted by hydrodynamic simulations of circumbinary accretion disks. The periodicity is likely caused by accretion rate modulation by a milli-parsec BSBH emitting GWs, dynamically coupled to the circumbinary accretion disk. A bursty hydrodynamic variability model is statistically preferred over a smooth, sinusoidal model expected from relativistic Doppler boost, a kinematic effect proposed for PG1302-102. Furthermore, the frequency dependence of the variability amplitudes disfavors Doppler boost, lending independent support to the circumbinary accretion variability hypothesis. Given our detection rate of one BSBH candidate from circumbinary accretion variability out of 625 quasars, it suggests that future large, sensitive synoptic surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time may be able to detect hundreds to thousands of candidate BSBHs from circumbinary accretion with direct implications for Laser Interferometer Space Antenna., Comment: Accepted to MNRAS, 17 pages, 9 figures

Periodically variable quasars have been suggested as close binary supermassive black holes. We present a systematic search for periodic light curves in 625 spectroscopically confirmed quasars with a median redshift of 1.8 in a 4.6 deg$^2$ overlapping region of the Dark Energy Survey Supernova (DES-SN) fields and the Sloan Digital Sky Survey Stripe 82 (SDSS-S82). Our sample has a unique 20-year long multi-color ($griz$) light curve enabled by combining DES-SN Y6 observations with archival SDSS-S82 data. The deep imaging allows us to search for periodic light curves in less luminous quasars (down to $r{\sim}$23.5 mag) powered by less massive black holes (with masses $\gtrsim10^{8.5}M_{\odot}$) at high redshift for the first time. We find five candidates with significant (at $>$99.74% single-frequency significance in at least two bands with a global p-value of $\sim$7$\times10^{-4}$--3$\times10^{-3}$ accounting for the look-elsewhere effect) periodicity with observed periods of $\sim$3--5 years (i.e., 1--2 years in rest frame) having $\sim$4--6 cycles spanned by the observations. If all five candidates are periodically variable quasars, this translates into a detection rate of ${\sim}0.8^{+0.5}_{-0.3}$% or ${\sim}1.1^{+0.7}_{-0.5}$ quasar per deg$^2$. Our detection rate is 4--80 times larger than those found by previous searches using shallower surveys over larger areas. This discrepancy is likely caused by differences in the quasar populations probed and the survey data qualities. We discuss implications on the future direct detection of low-frequency gravitational waves. Continued photometric monitoring will further assess the robustness and characteristics of these candidate periodic quasars to determine their physical origins., Comment: 21 pages, 8 figures, MNRAS accepted

We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev--Zeldovich (SZ) effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius for the complete galaxy sample is consistent with theoretical measurements from CDM-only simulations, and is located at $2.4^{+0.3}_{-0.4}$ Mpc $h^{-1}$. We split the sample based on galaxy color and find significant differences in the profile shapes. Red galaxies and those in the green valley show a splashback-like minimum in their slope profile consistent with theoretical predictions, while the bluest galaxies show a weak feature that appears at a smaller radius. We develop a mapping of galaxies to subhalos in $N$-body simulations by splitting subhalos based on infall time onto the cluster halos. We find that the location of the steepest slope and differences in the shapes of the profiles can be mapped to differences in the average time of infall of galaxies of different colors. The minima of the slope in the galaxy profiles trace a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time for galaxies of different colours, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in their clusters for over $3.2 ~\rm Gyrs$, green galaxies about $2.2 ~\rm Gyrs$, while blue galaxies have been accreted most recently and have not reached apocenter. Using the information from the complete radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching in clusters occurs after a delay of about a gigayear, and that galaxies quench rapidly thereafter with an exponential timescale of $0.6$ Gyr., Comment: 23 pages, 10 figures, to be submitted to ApJ

Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3$\times$2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modeling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude ($Q_1$) generally reflects the strength of baryon feedback. With the upper limit of $Q_1$ prior being bound by the Illustris feedback scenarios, we reach $\sim 20\%$ improvement in the constraint of $S_8=\sigma_8(\Omega_{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3$\times$2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative $Q_1$ prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than $2 \sigma$., Comment: 22 pages, 18 figures, 2 tables. accepted to MNRAS. A brief video summary of this paper is available at https://www.youtube.com/watch?v=QbeNwk5papU

The size of the dust torus in Active Galactic Nuclei (AGN) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to $\sim 70$ $z\lesssim 0.3$ AGN and quasars. Here we present first results of our dust RM program for distant quasars covered in the SDSS Stripe 82 region combining $\sim 20$-yr ground-based optical light curves with 10-yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1-band (3.4 $\mu$m) and $g$ band for 587 quasars over $0.3\lesssim z\lesssim 2$ ($\left\sim 0.8$) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter $\sim 0.17$ dex in lag) the IR lag-luminosity relation observed for $z<0.3$ AGN, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology, and strongly endorses a compelling science case for the combined 10-yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5-yr Nancy Grace Roman Space Telescope 2$\mu$m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7,384 quasars in our parent sample are made public with this work., Comment: Accepted for publication in ApJ

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$\beta$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{{\rm 2-12\,keV}}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$\beta$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$\sigma$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory., Comment: 13 pages, 8 figures, 1 table, accepted to MNRAS

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy's unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory., Comment: 15 pages, 6 figures, 1 table. Accepted for publication in ApJ Letters

There are several supervised machine learning methods used for the application of automated morphological classification of galaxies; however, there has not yet been a clear comparison of these different methods using imaging data, or a investigation for maximising their effectiveness. We carry out a comparison between several common machine learning methods for galaxy classification (Convolutional Neural Network (CNN), K-nearest neighbour, Logistic Regression, Support Vector Machine, Random Forest, and Neural Networks) by using Dark Energy Survey (DES) data combined with visual classifications from the Galaxy Zoo 1 project (GZ1). Our goal is to determine the optimal machine learning methods when using imaging data for galaxy classification. We show that CNN is the most successful method of these ten methods in our study. Using a sample of $\sim$2,800 galaxies with visual classification from GZ1, we reach an accuracy of $\sim$0.99 for the morphological classification of Ellipticals and Spirals. The further investigation of the galaxies that have a different ML and visual classification but with high predicted probabilities in our CNN usually reveals an the incorrect classification provided by GZ1. We further find the galaxies having a low probability of being either spirals or ellipticals are visually Lenticulars (S0), demonstrating that supervised learning is able to rediscover that this class of galaxy is distinct from both Es and Spirals. We confirm that $\sim$2.5\% galaxies are misclassified by GZ1 in our study. After correcting these galaxies' labels, we improve our CNN performance to an average accuracy of over 0.99 (accuracy of 0.994 is our best result)., Comment: 20 pages, 19 figures. Accepted to MNRAS

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed by more than 1.5 mag in the $g$ band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey (DES) epochs (separated by a few years in the quasar rest frame). The quasar sample covers a redshift range of $0.5 < z < 2.1$. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the $g$ band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on timescales of only months. Leveraging on the large dynamic range in continuum variability between the earlier SDSS and the new GMOS spectra, we explore the associated variations in the broad Mg II,$\lambda2798$ line, whose variability properties have not been well studied before. The broad Mg II flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the width (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad Mg II to estimate the black hole mass therefore introduces a luminosity-dependent bias., Comment: Submitted to MNRAS

Hydrogen-rich, core-collapse supernovae are typically divided into four classes: IIP, IIL, IIn, and IIb. In general, interaction with circumstellar material is only considered for Type IIn supernovae. However, recent hydrodynamic modeling of IIP and IIL supernovae requires circumstellar material to reproduce their early light curves. In this scenario, IIL supernovae experience large amounts of mass loss before exploding. We test this hypothesis on ASASSN-15oz, a Type IIL supernova. With extensive follow-up in the X- ray, UV, optical, IR, and radio we present our search for signs of interaction, and the mass-loss history indicated by their detection. We find evidence of short-lived intense mass-loss just prior to explosion from light curve modeling, amounting in 1.5 M$_{\odot}$ of material within 1800 R$_{\odot}$ of the progenitor. We also detect the supernova in the radio, indicating mass-loss rates of $10^{-6}-10^{-7}$ M$_{\odot}$ yr$^{-1}$ prior to the extreme mass-loss period. Our failure to detect the supernova in the X-ray and the lack of narrow emission lines in the UV, optical, and NIR do not contradict this picture and place an upper limit on the mass-loss rate outside the extreme period of $<10^{-4}$ M$_{\odot}$ yr$^{-1}$. This paper highlights the importance gathering comprehensive data on more Type II supernovae to enable detailed modeling of the progenitor and supernova which can elucidate their mass-loss histories and envelope structures and thus inform stellar evolution models., Comment: Submitted to MNRAS including referee comments

We present a suite of 18 synthetic sky catalogs designed to support science analysis of galaxies in the Dark Energy Survey Year 1 (DES Y1) data. For each catalog, we use a computationally efficient empirical approach, ADDGALS, to embed galaxies within light-cone outputs of three dark matter simulations that resolve halos with masses above ~5x10^12 h^-1 m_sun at z <= 0.32 and 10^13 h^-1 m_sun at z~2. The embedding method is tuned to match the observed evolution of galaxy counts at different luminosities as well as the spatial clustering of the galaxy population. Galaxies are lensed by matter along the line of sight --- including magnification, shear, and multiple images --- using CALCLENS, an algorithm that calculates shear with 0.42 arcmin resolution at galaxy positions in the full catalog. The catalogs presented here, each with the same LCDM cosmology (denoted Buzzard), contain on average 820 million galaxies over an area of 1120 square degrees with positions, magnitudes, shapes, photometric errors, and photometric redshift estimates. We show that the weak-lensing shear catalog, redMaGiC galaxy catalogs and redMaPPer cluster catalogs provide plausible realizations of the same catalogs in the DES Y1 data by comparing their magnitude, color and redshift distributions, angular clustering, and mass-observable relations, making them useful for testing analyses that use these samples. We make public the galaxy samples appropriate for the DES Y1 data, as well as the data vectors used for cosmology analyses on these simulations., Comment: 33 pages, 13 figures, catalogs will be made public upon publication; interested users may contact us beforehand

Using stellar population models, we predicted that the Dark Energy Survey (DES) - due to its special combination of area (5000 deg. sq.) and depth ($i = 24.3$) - would be in the position to detect massive ($\gtrsim 10^{11}$ M$_{\odot}$) galaxies at $z \sim 4$. We confront those theoretical calculations with the first $\sim 150$ deg. sq. of DES data reaching nominal depth. From a catalogue containing $\sim 5$ million sources, $\sim26000$ were found to have observed-frame $g-r$ vs $r-i$ colours within the locus predicted for $z \sim 4$ massive galaxies. We further removed contamination by stars and artefacts, obtaining 606 galaxies lining up by the model selection box. We obtained their photometric redshifts and physical properties by fitting model templates spanning a wide range of star formation histories, reddening and redshift. Key to constrain the models is the addition, to the optical DES bands $g$, $r$, $i$, $z$, and $Y$, of near-IR $J$, $H$, $K_{s}$ data from the Vista Hemisphere Survey. We further applied several quality cuts to the fitting results, including goodness of fit and a unimodal redshift probability distribution. We finally select 233 candidates whose photometric redshift probability distribution function peaks around $z\sim4$, have high stellar masses ($\log($M$^{*}$/M$_{\odot})\sim 11.7$ for a Salpeter IMF) and ages around 0.1 Gyr, i.e. formation redshift around 5. These properties match those of the progenitors of the most massive galaxies in the local universe. This is an ideal sample for spectroscopic follow-up to select the fraction of galaxies which is truly at high redshift. These initial results and those at the survey completion, which we shall push to higher redshifts, will set unprecedented constraints on galaxy formation, evolution, and the re-ionisation epoch., Comment: 24 pages (41 with appendix), 16 figures, MNRAS in press

Transneptunian objects (TNOs) are a source of invaluable information to access the history and evolution of the outer solar system. However, observing these faint objects is a difficult task. As a consequence, important properties such as size and albedo are known for only a small fraction of them. Now, with the results from deep sky surveys and the Gaia space mission, a new exciting era is within reach as accurate predictions of stellar occultations by numerous distant small solar system bodies become available. From them, diameters with kilometer accuracies can be determined. Albedos, in turn, can be obtained from diameters and absolute magnitudes. We use observations from the Dark Energy Survey (DES) from November 2012 until February 2016, amounting to 4292847 CCD frames. We searched them for all known small solar system bodies and recovered a total of 202 TNOs and Centaurs, 63 of which have been discovered by the DES collaboration until the date of this writing. Their positions were determined using the Gaia Data Release 2 as reference and their orbits were refined. Stellar occultations were then predicted using these refined orbits plus stellar positions from Gaia. These predictions are maintained, and updated, in a dedicated web service. The techniques developed here are also part of an ambitious preparation to use the data from the Large Synoptic Survey Telescope (LSST), that expects to obtain accurate positions and multifilter photometry for tens of thousands of TNOs., Comment: 25 pages, submitted to Astronomical Journal

It is now established that there is a dependence of the luminosity of type Ia supernovae (SNe Ia) on environment: SNe Ia in young, star-forming, metal-poor stellar populations appear fainter after light-curve shape corrections than those in older, passive, metal-rich environments. This is accounted for in cosmological studies using a global property of the SN host galaxy, typically the host galaxy stellar mass. However, recent low-redshift studies suggest that this effect manifests itself most strongly when using the local star-formation rate (SFR) at the SN location, rather than the global SFR or stellar mass of the host galaxy. At high-redshift, such local SFRs are difficult to determine; here, we show that an equivalent 'local' correction can be made by restricting the SN Ia sample in globally star-forming host galaxies to a low-mass host galaxy subset ($\le10^{10} M_{\odot}$). Comparing this sample of SNe Ia (in locally star-forming environments) to those in locally passive host galaxies, we find that SNe Ia in locally star-forming environments are $0.081\pm0.018$ mag fainter ($4.5\sigma$), consistent with the result reported by Rigault et al. (2015), but our conclusion is based on a sample ~5 times larger over a wider redshift range. This is a larger difference than when splitting the SN Ia sample based on global host galaxy SFR or host galaxy stellar mass. This method can be used in ongoing and future high-redshift SN surveys, where local SN Ia environments are difficult to determine., Comment: 16 pages, 3 figures, and 5 tables, accepted for publication in ApJ

This paper focuses on performance analyzing and dynamic modeling of the current grid-tied fixed array 6.84kW solar photovoltaic system located at Florida Atlantic University (FAU). A battery energy storage system is designed and applied to improve the systems stability and reliability. An overview of the entire system and its PV module are presented. In sequel, the corresponding I-V and P-V curves are obtained using MATLAB-Simulink package. Actual data was collected and utilized for the modeling and simulation of the system. In addition, a grid- connected PV/Battery system with Maximum Power Point Tracking (MPPT) controller is modeled to analyze the system performance that has been evaluated under two different test conditions: (1) PV power production is higher than the load demand (2) PV generated power is less than required load. A battery system has also been sized to provide smoothing services to this array. The simulation results show the effective of the proposed method. This system can be implemented in developing countries with similar weather conditions to Florida., Comment: 6 pages, 14 figures, PVSC 2017

We present a novel approach for creating science-ready catalogs through a software infrastructure developed for the Dark Energy Survey (DES). We integrate the data products released by the DES Data Management and additional products created by the DES collaboration in an environment known as DES Science Portal. Each step involved in the creation of a science-ready catalog is recorded in a relational database and can be recovered at any time. We describe how the DES Science Portal automates the creation and characterization of lightweight catalogs for DES Year 1 Annual Release, and show its flexibility in creating multiple catalogs with different inputs and configurations. Finally, we discuss the advantages of this infrastructure for large surveys such as DES and the Large Synoptic Survey Telescope. The capability of creating science-ready catalogs efficiently and with full control of the inputs and configurations used is an important asset for supporting science analysis using data from large astronomical surveys., Comment: The second paper of the series about the DES Science Portal, submitted to the Astronomy & Computing journal. It shows the infrastructure to create science-ready catalogs from DES photometric data and ancillary maps. This is the version accepted by the journal

A significant challenge facing photometric surveys for cosmological purposes is the need to produce reliable redshift estimates. The estimation of photometric redshifts (photo-zs) has been consolidated as the standard strategy to bypass the high production costs and incompleteness of spectroscopic redshift samples. Training-based photo-z methods require the preparation of a high-quality list of spectroscopic redshifts, which needs to be constantly updated. The photo-z training, validation, and estimation must be performed in a consistent and reproducible way in order to accomplish the scientific requirements. To meet this purpose, we developed an integrated web-based data interface that not only provides the framework to carry out the above steps in a systematic way, enabling the ease testing and comparison of different algorithms, but also addresses the processing requirements by parallelizing the calculation in a transparent way for the user. This framework called the Science Portal (hereafter Portal) was developed in the context the Dark Energy Survey (DES) to facilitate scientific analysis. In this paper, we show how the Portal can provide a reliable environment to access vast data sets, provide validation algorithms and metrics, even in the case of multiple photo-zs methods. It is possible to maintain the provenance between the steps of a chain of workflows while ensuring reproducibility of the results. We illustrate how the Portal can be used to provide photo-z estimates using the DES first year (Y1A1) data. While the DES collaboration is still developing techniques to obtain more precise photo-zs, having a structured framework like the one presented here is critical for the systematic vetting of DES algorithmic improvements and the consistent production of photo-zs in the future DES releases., Comment: The first paper of the series about the DES Science Portal, submitted to the Astronomy & Computing journal. It shows the infrastructure of pipelines related to the production of photometric redshifts to feed science-ready catalogs

Presently a ${>}3\sigma$ tension exists between values of the Hubble constant $H_0$ derived from analysis of fluctuations in the Cosmic Microwave Background by Planck, and local measurements of the expansion using calibrators of type Ia supernovae (SNe Ia). We perform a blinded reanalysis of Riess et al. 2011 to measure $H_0$ from low-redshift SNe Ia, calibrated by Cepheid variables and geometric distances including to NGC 4258. This paper is a demonstration of techniques to be applied to the Riess et at. 2016 data. Our end-to-end analysis starts from available CfA3 and LOSS photometry, providing an independent validation of Riess et al. 2011. We obscure the value of $H_0$ throughout our analysis and the first stage of the referee process, because calibration of SNe Ia requires a series of often subtle choices, and the potential for results to be affected by human bias is significant. Our analysis departs from that of Riess et al. 2011 by incorporating the covariance matrix method adopted in SNLS and JLA to quantify SN Ia systematics, and by including a simultaneous fit of all SN Ia and Cepheid data. We find $H_0 = 72.5 \pm 3.1$ (stat) $\pm 0.77$ (sys) km s$^{-1}$ Mpc$^{-1}$ with a three-galaxy (NGC 4258+LMC+MW) anchor. The relative uncertainties are 4.3% statistical, 1.1% systematic, and 4.4% total, larger than in Riess et al. 2011 (3.3% total) and the Efstathiou 2014 reanalysis (3.4% total). Our error budget for $H_0$ is dominated by statistical errors due to the small size of the supernova sample, whilst the systematic contribution is dominated by variation in the Cepheid fits, and for the SNe Ia, uncertainties in the host galaxy mass dependence and Malmquist bias., Comment: 38 pages, 13 figures, 13 tables; accepted for publication in MNRAS

We perform a systematic search for long-term extreme variability quasars (EVQs) in the overlapping Sloan Digital Sky Survey (SDSS) and 3-Year Dark Energy Survey (DES) imaging, which provide light curves spanning more than 15 years. We identified ~1000 EVQs with a maximum g band magnitude change of more than 1 mag over this period, about 10% of all quasars searched. The EVQs have L_bol~10^45-10^47 erg/s and L/L_Edd~0.01-1. Accounting for selection effects, we estimate an intrinsic EVQ fraction of ~30-50% among all g<~22 quasars over a baseline of ~15 years. These EVQs are good candidates for so-called "changing-look quasars", where a spectral transition between the two types of quasars (broad-line and narrow-line) is observed between the dim and bright states. We performed detailed multi-wavelength, spectral and variability analyses for the EVQs and compared to their parent quasar sample. We found that EVQs are distinct from a control sample of quasars matched in redshift and optical luminosity: (1) their UV broad emission lines have larger equivalent widths; (2) their Eddington ratios are systematically lower; and (3) they are more variable on all timescales. The intrinsic difference in quasar properties for EVQs suggest that internal processes associated with accretion are the main driver for the observed extreme long-term variability. However, despite their different properties, EVQs seem to be in the tail of a continuous distribution of quasar properties, rather than standing out as a distinct population. We speculate that EVQs are normal quasars accreting at relatively low accretion rates, where the accretion flow is more likely to experience instabilities that drive the factor of few changes in flux on multi-year timescales., Comment: submitted to ApJ; comments welcome

We report on SPT-CLJ2011-5228, a giant system of arcs created by a cluster at $z=1.06$. The arc system is notable for the presence of a bright central image. The source is a Lyman Break galaxy at $z_s=2.39$ and the mass enclosed within the 14 arc second radius Einstein ring is $10^{14.2}$ solar masses. We perform a full light profile reconstruction of the lensed images to precisely infer the parameters of the mass distribution. The brightness of the central image demands that the central total density profile of the lens be shallow. By fitting the dark matter as a generalized Navarro-Frenk-White profile---with a free parameter for the inner density slope---we find that the break radius is $270^{+48}_{-76}$ kpc, and that the inner density falls with radius to the power $-0.38\pm0.04$ at 68 percent confidence. Such a shallow profile is in strong tension with our understanding of relaxed cold dark matter halos; dark matter only simulations predict the inner density should fall as $r^{-1}$. The tension can be alleviated if this cluster is in fact a merger; a two halo model can also reconstruct the data, with both clumps (density going as $r^{-0.8}$ and $r^{-1.0}$) much more consistent with predictions from dark matter only simulations. At the resolution of our Dark Energy Survey imaging, we are unable to choose between these two models, but we make predictions for forthcoming Hubble Space Telescope imaging that will decisively distinguish between them., Comment: 13 Pages. Accepted for publication in ApJ

We report the observation and physical characterization of the possible dwarf planet \UZ\ ("DeeDee"), a dynamically detached trans-Neptunian object discovered at 92 AU. This object is currently the second-most distant known trans-Neptunian object with reported orbital elements, surpassed in distance only by the dwarf planet Eris. The object was discovered with an $r$-band magnitude of 23.0 in data collected by the Dark Energy Survey between 2014 and 2016. Its 1140-year orbit has $(a,e,i) = (109~\mathrm{AU}, 0.65, 26.8^{\circ})$. It will reach its perihelion distance of 38 AU in the year 2142. Integrations of its orbit show it to be dynamically stable on Gyr timescales, with only weak interactions with Neptune. We have performed followup observations with ALMA, using 3 hours of on-source integration time to measure the object's thermal emission in the Rayleigh-Jeans tail. The signal is detected at 7$\sigma$ significance, from which we determine a $V$-band albedo of $13.1^{+3.3}_{-2.4}\mathrm{(stat)}^{+2.0}_{-1.4}\mathrm{(sys)}$ percent and a diameter of $635^{+57}_{-61}\mathrm{(stat)}^{+32}_{-39}\mathrm{(sys)}$~km, assuming a spherical body with uniform surface properties., Comment: 13 pages, 4 figures, accepted by ApJ Letters

Host galaxy identification is a crucial step for modern supernova (SN) surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey Telescope (LSST), which will discover SNe by the thousands. Spectroscopic resources are limited, so in the absence of real-time SN spectra these surveys must rely on host galaxy spectra to obtain accurate redshifts for the Hubble diagram and to improve photometric classification of SNe. In addition, SN luminosities are known to correlate with host-galaxy properties. Therefore, reliable identification of host galaxies is essential for cosmology and SN science. We simulate SN events and their locations within their host galaxies to develop and test methods for matching SNe to their hosts. We use both real and simulated galaxy catalog data from the Advanced Camera for Surveys General Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe residing in undetected faint hosts into our analysis, with an assumed hostless rate of 5%. Our fully automated algorithm is run on catalog data and matches SNe to their hosts with 91% accuracy. We find that including a machine learning component, run after the initial matching algorithm, improves the accuracy (purity) of the matching to 97% with a 2% cost in efficiency (true positive rate). Although the exact results are dependent on the details of the survey and the galaxy catalogs used, the method of identifying host galaxies we outline here can be applied to any transient survey., Comment: 22 pages, 18 figures; Accepted by AJ, revised to incorporate referee comments

Using the largest single-survey sample of Type Ia supernovae (SNe Ia) to date, we study the relationship between properties of SNe Ia and those of their host galaxies, focusing primarily on correlations with Hubble residuals (HR). Our sample consists of 345 photometrically-classified or spectroscopically-confirmed SNeIa discovered as part of the SDSS-II Supernova Survey (SDSS-SNS). This analysis utilizes host-galaxy spectroscopy obtained during the SDSS-I/II spectroscopic survey and from an ancillary program on the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) that obtained spectra for nearly all host galaxies of SDSS-II SN candidates. In addition, we use photometric host-galaxy properties from the SDSS-SNS data release (Sako et al. 2014) such as host stellar mass and star-formation rate. We confirm the well-known relation between HR and host-galaxy mass and find a 3.6{\sigma} significance of a non-zero linear slope. We also recover correlations between HR and host-galaxy gas-phase metallicity and specific star-formation rate as they are reported in the literature. With our large dataset, we examine correlations between HR and multiple host-galaxy properties simultaneously and find no evidence of a significant correlation. We also independently analyze our spectroscopically-confirmed and photometrically-classified SNe Ia and comment on the significance of similar combined datasets for future surveys., Comment: Accepted for publication in ApJ