
A significant challenge facing photometric surveys for cosmological purposes
is the need to produce reliable redshift estimates. The estimation of
photometric redshifts (photozs) has been consolidated as the standard strategy
to bypass the high production costs and incompleteness of spectroscopic
redshift samples. Trainingbased photoz methods require the preparation of a
highquality list of spectroscopic redshifts, which needs to be constantly
updated. The photoz 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 webbased 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 photozs
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 photoz estimates using the DES first year
(Y1A1) data. While the DES collaboration is still developing techniques to
obtain more precise photozs, having a structured framework like the one
presented here is critical for the systematic vetting of DES algorithmic
improvements and the consistent production of photozs in the future DES
releases.

We carry out a study of the statistical distribution of rainfall
precipitation data for 20 cites in India. We have determined the bestfit
probability distribution for these cities from the monthly precipitation data
spanning 100 years of observations from 1901 to 2002. To fit the observed data,
we considered 10 different distributions. The efficacy of the fits for these
distributions was evaluated using four empirical nonparametric goodnessoffit
tests namely KolmogorovSmirnov, AndersonDarling, ChiSquare, Akaike
information criterion, and Bayesian Information criterion. Finally, the
bestfit distribution using each of these tests were reported, by combining the
results from the model comparison tests. We then find that for most of the
cities, Generalized ExtremeValue Distribution or Inverse Gaussian Distribution
most adequately fits the observed data.

Recently, Odrzywolek and Rafelski (arXiv:1612.03556) have found three
distinct categories of exoplanets, when they are classified based on density.
We first carry out a similar classification of exoplanets according to their
density using the Gaussian Mixture Model, followed by information theoretic
criterion (AIC and BIC) to determine the optimum number of components. Such a
onedimensional classification favors two components using AIC and three using
BIC, but the statistical significance from both the tests is not significant
enough to decisively pick the best model between two and three components. We
then extend this GMMbased classification to two dimensions by using both the
density and the Earth similarity index (arXiv:1702.03678), which is a measure
of how similar each planet is compared to the Earth. For this twodimensional
classification, both AIC and BIC provide decisive evidence in favor of three
components.

We present a novel approach for creating scienceready 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 scienceready 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
scienceready 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.

According to Benford's law, the most significant digit in many datasets is
not uniformly distributed, but obeys a well defined power law distribution with
smaller digits appearing more often. Among one of the myriad particle physics
datasets available, we find that the leading decimal digit for the $\tau$
lepton branching fraction shows marginal disagreement with the logarithmic
behavior expected from the Benford distribution. We quantify the deviation from
Benford's law using a $\chi^2$ function valid for binomial data, and obtain a
$\chi^2$ value of 16.9 for nine degrees of freedom, which gives a $p$value of
about 5%, corresponding to a 1.6$\sigma$ disagreement. We also checked that the
disagreement persists under scaling the branching fractions, as well as by
redoing the analysis in a numerical system with a base different from 10. Among
all the digits, `9' shows the largest discrepancy with an excess of $4\sigma$.
This discrepancy is because the digit `9' is repeated for three distinct groups
of correlated modes, with each group having a frequency of two or three,
leading to doublecounting. If we count each group of correlated modes only
once, the discrepancy for this digit also disappears and we get pristine
agreement with Benford distribution.

The classification of the host galaxies of the radio sources in the 25
deg$^2$ ultimate XMM extragalactic survey south field (XXLS) is presented.
XXLS was surveyed at 2.1 GHz with the Australia Telescope Compact Array (ATCA)
and is thus far the largest area radio survey conducted down to rms flux
densities of $\sigma \sim 41$ $\mu$Jy beam$^{1}$. Of the 6287 radio sources in
XXLS, 4758 (75.7%) were crossmatched to an optical counterpart using the
likelihood ratio technique. There are 1110 spectroscopic redshifts and 3648
photometric redshifts available for the counterparts, of which 99.4% exist out
to $z \sim 4$. A number of multiwavelength diagnostics, including Xray
luminosities, midinfrared colours, spectral energy distribution fits, radio
luminosities, and optical emission lines and colours, were used to classify the
sources into three types: lowexcitation radio galaxies (LERGs),
highexcitation radio galaxies (HERGs), and starforming galaxies (SFGs). The
final sample contains 1729 LERGs (36.3%), 1159 radioloud HERGs (24.4%), 296
radioquiet HERGs (6.2%), 558 SFGs (11.7%), and 1016 unclassified sources
(21.4%). The LERGs tend to exist in the most massive galaxies with low star
formation rates and redder colours, whereas the HERGs and SFGs exist in
galaxies of lower mass, higher star formation rates, and bluer colours. The
fraction of blue host galaxies is higher for radioquiet HERGs than for
radioloud HERGs. LERGs and radioloud HERGs are found at all radio
luminosities, but radioloud HERGs tend to be more radio luminous than LERGs at
a given redshift. These results are consistent with the emerging picture in
which LERGs exist in the most massive quiescent galaxies typically found in
clusters with hot Xray halos and HERGs are associated with ongoing star
formation in their host galaxies via the accretion of cold gas.

Many problems in Astrophysics involve using Bayesian Inference to deal with
problems of parameter estimation and model selection. In this paper, we
introduce Variational Inference to solve these problems and compare how the
results hold up to Markov Chain Monte Carlo which is the most common method.
Variational Inference converts the inference problem into an optimization
problem by approximating the posterior from a known family of distributions and
using KullbackLeibler divergence to measure closeness. Variational Inference
takes advantage of fast optimization techniques which make it ideal to deal
with large datasets and also makes it trivial to parallelize. As a proof of
principle, we apply Variational Inference for parameter estimation and model
comparison to four different problems in astrophysics where MCMC techniques
were previously used: measuring exoplanet orbital parameters from radial
velocity data, tests of periodicities in measurements of $G$, significance of a
turnover in the spectral lag data of GRB 160625B , and estimating the mass of a
galaxy cluster using weak lensing. We find that Variational Inference is much
faster than MCMC for these problems.

We construct the error distributions for the galactic rotation speed
($\Theta_0$) using 137 data points from measurements compiled in De Grijs et
al. (arXiv:1709.02501), with all observations normalized to the galactocentric
distance of 8.3 kpc. We then checked (using the same procedures as in works by
Ratra et al) if the errors constructed using the weighted mean and the median
as the estimate, obey Gaussian statistics. We find using both these estimates
that they have much wider tails than a Gaussian distribution. We also tried to
fit the data to three other distributions: Cauchy, doubleexponential, and
Studentst. The best fit is obtained using the Students$t$ distribution for
$n=2$ using the median value as the central estimate, corresponding to a
$p$value of 0.1. We also calculate the median value of $\Theta_0$ using all
the data as well as using the median of each set of measurements based on the
tracer population used. Because of the nongaussianity of the residuals, we
point out that the subgroup median value, given by $\Theta_{med}=219.65$ km/sec
should be used as the central estimate for $\Theta_0$.

We evaluate the performance of four different machine learning (ML)
algorithms: an Artificial Neural Network MultiLayer Perceptron (ANN MLP ),
Adaboost, Gradient Boosting Classifier (GBC), XGBoost, for the separation of
pulsars from radio frequency interference (RFI) and other sources of noise,
using a dataset obtained from the postprocessing of a pulsar search pi peline.
This dataset was previously used for crossvalidation of the SPINNbased
machine learning engine, used for the reprocessing of HTRUS survey data
arXiv:1406.3627. We have used Synthetic Minority Oversampling Technique
(SMOTE) to deal with high class imbalance in the dataset. We report a variety
of quality scores from all four of these algorithms on both the nonSMOTE and
SMOTE datasets. For all the above ML methods, we report high accuracy and
Gmean in both the nonSMOTE and SMOTE cases. We study the feature importances
using Adaboost, GBC, and XGBoost and also from the minimum Redundancy Maximum
Relevance approach to report algorithmagnostic feature ranking. From these
methods, we find that the signal to noise of the folded profile to be the best
feature. We find that all the ML algorithms report FPRs about an order of
magnitude lower than the corresponding FPRs obtained in arXiv:1406.3627, for
the same recall value.

On August 17, 2017 the LIGO interferometers detected the gravitational wave
(GW) signal (GW170817) from the coalescence of binary neutron stars. This
signal was also simultaneously seen throughout the electromagnetic (EM)
spectrum from radio waves to gammarays. We point out that this simultaneous
detection of GW and EM signals rules out a class of modified gravity theories,
termed ``dark matter emulators,'' which dispense with the need for dark matter
by making ordinary matter couple to a different metric from that of GW. We
discuss other kinds of modified gravity theories which dispense with the need
for dark matter and are still viable. This simultaneous observation also
provides the first observational test of Einstein's Weak Equivalence Principle
(WEP) between gravitons and photons. We estimate the Shapiro time delay due to
the gravitational potential of the total dark matter distribution along the
line of sight (complementary to the calculation in arXiv:1710.05834) to be
about 400 days. Using this estimate for the Shapiro delay and from the time
difference of 1.7 seconds between the GW signal and gammarays, we can
constrain violations of WEP using the parameterized postNewtonian (PPN)
parameter $\gamma$, and is given by $\gamma_{\rm {GW}}\gamma_{\rm{EM}}<9.8
\times 10^{8}$.

To date, the only limit on graviton mass using galaxy clusters was obtained
by Goldhaber and Nieto in 1974, using the fact that the orbits of galaxy
clusters are bound and closed, and extend up to 580 kpc. From positing that
only a Newtonian potential gives rise to such stable bound orbits, a limit on
the graviton mass $m_g<10^{29}$ eV was obtained (PRD 9,1119, 1974). Recently,
it has been shown that one can obtain closed bound orbits for Yukawa potential
(arXiv:1705.02444), thus invalidating the main \emph{ansatz} used in Goldhaber
and Nieto to obtain the graviton mass bound. In order to obtain a revised
estimate using galaxy clusters, we use dynamical mass models of the Abell 1689
(A1689) galaxy cluster to check their compatibility with a Yukawa gravitational
potential. We assume mass models for the gas, dark matter, and galaxies for
A1689 from arXiv:1703.10219 and arXiv:1610.01543, who used this cluster to test
various alternate gravity theories, which dispense with the need for dark
matter. We quantify the deviations in the acceleration profile using these mass
models assuming a Yukawa potential and that obtained assuming a Newtonian
potential by calculating the $\chi^2$ residuals between the two profiles. Our
estimated bound on the graviton mass ($m_g$) is thereby given by, $m_g < 1.37
\times 10^{29}$ eV or in terms of the graviton Compton wavelength of,
$\lambda_g>9.1 \times 10^{19}$ km at 90\% confidence level.

We calculate the total galactic Shapiro delay to the Crab pulsar by including
the contributions from the dark matter as well as baryonic matter along the
line of sight. The total delay due to dark matter potential is about 3.4 days.
For baryonic matter, we included the contributions from both the bulge and the
disk, which are approximately 0.12 and 0.32 days respectively. The total delay
from all the matter distribution is therefore 3.84 days. We also calculate the
limit on violations of Einstein's equivalence principle by using observations
of "nanoshot" giant pulses from the Crab pulsar with timedelay $<0.4$~ns as
well as using time differences between radio and optical photons observed from
this pulsar. Using the former, we obtain a limit on violation of Einstein's
equivalence principle in terms of the PPN parameter $\Delta \gamma < 2.41\times
10^{15}$. From the timedifference between simultaneous optical and radio
observations, we get $\Delta \gamma < 1.54\times 10^{9}$. We also point out
differences in our calculation of Shapiro delay and that from two recent papers
(arXiv:1612.00717 and arXiv:1608.07657), which used the same observations to
obtain a corresponding limit on $\Delta \gamma$.

We apply the generalized LombScargle (LS) periodogram to independently
confirm the claim by Sturrock et al (arXiv:1605.03088) of oscillation at a
frequency of 11/year in the decay rates of $\rm{^{90}Sr/^{90}Y}$ from
measurements at the Physikalisch Technische Bundesanstalt (PTB), which however
has been disputed by Kossert and Nahle (arXiv:1407.2493). For this analysis, we
made two different {\it ansatze} for the errors. For each peak in the LS
periodogram, we evaluate the statistical significance using nonparametric
bootstrap resampling. We find using both these error models that none of the
peaks in the LS periodogram are statistically significant, thus invalidating
the claim by Sturrock et al of any oscillations in the decay rates.

GammaRay Bursts (GRBs) have been conventionally bifurcated into two distinct
categories: ``short'' and ``long'' with durations less than and greater than
two seconds respectively. However, there is a lot of literature (although with
conflicting results) regarding the existence of a third intermediate class. To
investigate this issue, we extend a recent study (arXiv:1612.08235) on
classification of GRBs to two dimensions by incorporating the GRB hardness in
addition to the observed durations. We carry out this unified analysis on GRB
datasets from four detectors, viz. BATSE, RHESSI, Swift (observed and intrinsic
frame), and FermiGBM. We consider the duration and hardness features in
logscale for each of these datasets and determine the bestfit parameters
using Gaussian Mixture Model. This is followed by information theoretic
criterion (AIC and BIC) to determine if a threecomponent fit is favored
compared to a twocomponent one or viceversa. For BATSE, we find that both AIC
and BIC show preference for three components with decisive significance. For
Fermi and RHESSI, both AIC and BIC show preference for two components, although
the significance is marginal from AIC, but decisive using BIC. For Swift
dataset in both the observed and rest frame, we find that three components are
favored according to AIC with decisive significance, and two are preferred with
BIC with marginal to strong significance.

Recently Wei et al (arXiv:1612.09425) have found evidence for a transition
from positive time lags to negative time lags in the spectral lag data of GRB
160625B. They have fit these observed lags to a sum of two components: an
assumed functional form for intrinsic time lag due to astrophysical mechanisms
and an energydependent speed of light due to quadratic and linear Loren tz
invariance violation (LIV) models. Here, we examine the statistical
significance of the evidence for a transition to nega tive time lags. Such a
transition, even if present in GRB 160625B, cannot be due to an energy
dependent speed of light as th is would contradict previous limits by some 34
orders of magnitude, and must therefore be of intrinsic astrophysical origin .
We use three different model comparison techniques: a frequentist test and two
information based criteria (AIC and BIC). From the frequentist model comparison
test, we find that the evidence for transition in the spectral lag data is
favored at $3.05\sigma$ and $3.74\sigma$ for the linear and quadratic models
respectively. We find that $\Delta$AIC and $\Delta$BIC have values $\gtrsim$ 10
for the spectral lag transition that was motivated as being due to quadratic
Lorentz invariance vio lating model pointing to "decisive evidence". We note
however that none of the three models (including the model of intr insic
astrophysical emission) provide a good fit to the data.

We perform a systematic search for longterm extreme variability quasars
(EVQs) in the overlapping Sloan Digital Sky Survey (SDSS) and 3Year 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^4510^47 erg/s and L/L_Edd~0.011. Accounting for selection effects,
we estimate an intrinsic EVQ fraction of ~3050% among all g<~22 quasars over a
baseline of ~15 years. These EVQs are good candidates for socalled
"changinglook quasars", where a spectral transition between the two types of
quasars (broadline and narrowline) is observed between the dim and bright
states. We performed detailed multiwavelength, 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 longterm 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 multiyear timescales.

We report on SPTCLJ20115228, 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 NavarroFrenkWhite profilewith a
free parameter for the inner density slopewe 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.

The 2.1 GHz radio source catalogue of the 25 deg$^2$ ultimate XMM
extragalactic survey south (XXLS) field, observed with the Australia Telescope
Compact Array (ATCA), is presented. The final radio mosaic achieved a
resolution of $\sim$$4.8"$ and a median rms noise of $\sigma \approx41$
$\mu$Jy/beam. To date, this is the largest area radio survey to reach this flux
density level. A total of 6350 radio components above 5$\sigma$ are included in
the component catalogue, 26.4% of which are resolved. Of these components, 111
were merged together to create 48 multiplecomponent radio sources, resulting
in a total of 6287 radio sources in the source catalogue, 25.9% of which were
resolved. A survival analysis revealed that the median spectral index of the
Sydney University Molonglo Sky Survey (SUMSS) 843 MHz sources in the field is
$\alpha$ = $$0.75, consistent with the values of $0.7$ to $0.8$ commonly
used to characterise radio spectral energy distributions of active galactic
nuclei (AGN). The 2.1 GHz and 1.4 GHz differential radio source counts are
presented and compared to other 1.4 GHz radio surveys. The XXLS source counts
show good agreement with the other surveys.

GammaRay Bursts (GRBs) have been conventionally bifurcated into two distinct
categories dubbed "short" and "long", depending on whether their durations are
less than or greater than two seconds respectively. However, many authors have
pointed to the existence of a third class of GRBs with mean durations
intermediate between the short and long GRBs. Here, we apply multiple model
comparison techniques to verify these claims. For each category, we obtain the
bestfit parameters by maximizing a likelihood function based on a weighted
superposition of lognormal distributions. We then do modelcomparison between
each of these hypotheses by comparing the chisquare probabilities, Akaike
Information criterion (AIC), and Bayesian Information criterion (BIC). We
uniformly apply these techniques to GRBs from Swift (both observer and
intrinsic frame), BATSE, BeppoSAX, and FermiGBM. We find that the Swift GRB
distributions (in the observer frame) show evidence for three categories at
about $2.4\sigma$ from difference in chisquares and show decisive evidence in
favor of the two components using both AIC and BIC. For all the other datasets,
evidence for three components is either very marginal or not favored.

Euclid is a Europeled cosmology space mission dedicated to a visible and
near infrared survey of the entire extragalactic sky. Its purpose is to deepen
our knowledge of the dark content of our Universe. After an overview of the
Euclid mission and science, this contribution describes how the community is
getting organized to face the data analysis challenges, both in software
development and in operational data processing matters. It ends with a more
specific account of some of the main contributions of the Swiss Science Data
Center (SDCCH).

Anderson et al have claimed to find evidence for periodic sinusoidal
variations (period=5.9 years) in measurements of Newton's Gravitational
constant. These claims have been disputed by Pitkin. Using Bayesian model
comparison, he argues that a model with an unknown Gaussian noise component is
favored over any periodic variations by more than $e^{30}$. We reexamine the
claims of Anderson et al using frequentist model comparison tests, both with
and without errors in the measurement times. Our findings lend support to
Pitkin's claim that a constant term along with an unknown systematic offset
provides a better fit to the measurements of Newton's constant, compared to any
sinusoidal variations.

Robustness of any statistics depends upon the number of assumptions it makes
about the measured data. We point out the advantages of median statistics using
toy numerical experiments and demonstrate its robustness, when the number of
assumptions we can make about the data are limited. We then apply the median
statistics technique to obtain estimates of two constants of nature, Hubble
Constant ($H_0$) and Newton's Gravitational Constant($G$), both of which show
significant differences between different measurements. For $H_0$, we update
the analysis done by Chen and Ratra (2011) and Gott et al. (2001) using $576$
measurements. We find after grouping the different results according to their
primary type of measurement, the median estimates are given by
$H_0=72.5^{+2.5}_{8}$ km/sec/Mpc with errors corresponding to 95% c.l.
(2$\sigma$) and $G=6.674702^{+0.0014}_{0.0009} \times 10^{11} \mathrm{N
m^{2}kg^{2}}$ corresponding to 68% c.l. (1$\sigma$).

We present the discovery and preliminary characterization of a
gravitationally lensed quasar with a source redshift $z_{s}=2.74$ and image
separation of $2.9"$ lensed by a foreground $z_{l}=0.40$ elliptical galaxy.
Since the images of gravitationally lensed quasars are the superposition of
multiple point sources and a foreground lensing galaxy, we have developed a
morphology independent multiwavelength approach to the photometric selection
of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised
machine learning. Using this technique and $gi$ multicolour photometric
observations from the Dark Energy Survey (DES), near IR $JK$ photometry from
the VISTA Hemisphere Survey (VHS) and WISE mid IR photometry, we have
identified a candidate system with two catalogue components with $i_{AB}=18.61$
and $i_{AB}=20.44$ comprised of an elliptical galaxy and two blue point
sources. Spectroscopic followup with NTT and the use of an archival AAT
spectrum show that the point sources can be identified as a lensed quasar with
an emission line redshift of $z=2.739\pm0.003$ and a foreground early type
galaxy with $z=0.400\pm0.002$. We model the system as a single isothermal
ellipsoid and find the Einstein radius $\theta_E \sim 1.47"$, enclosed mass
$M_{enc} \sim 4 \times 10^{11}$M$_{\odot}$ and a time delay of $\sim$52 days.
The relatively wide separation, month scale time delay duration and high
redshift make this an ideal system for constraining the expansion rate beyond a
redshift of 1.

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 realtime 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 hostgalaxy 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.

We present the discovery of a z=0.65 lowionization broad absorption line
(LoBAL) quasar in a poststarburst galaxy in data from the Dark Energy Survey
(DES) and spectroscopy from the Australian Dark Energy Survey (OzDES). LoBAL
quasars are a minority of all BALs, and rarer still is that this object also
exhibits broad FeII (an FeLoBAL) and Balmer absorption. This is the first BAL
quasar that has signatures of recently truncated star formation, which we
estimate ended about 40 Myr ago. The characteristic signatures of an FeLoBAL
require high column densities, which could be explained by the emergence of a
young quasar from an early, dustenshrouded phase, or by clouds compressed by a
blast wave. The age of the starburst component is comparable to estimates of
the lifetime of quasars, so if we assume the quasar activity is related to the
truncation of the star formation, this object is better explained by the blast
wave scenario.