We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code AREPO-RT, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically-symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in AREPO-RT, which treats UV and optical photons as single-scattered and IR photons as multi-scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

Simulating the dust content of galaxies and their surrounding gas is challenging due to the wide range of physical processes affecting the dust evolution. Here we present cosmological hydrodynamical simulations of a cluster of galaxies, $M_\text{200,crit}=6 \times 10^{14}\,{\rm M_\odot}$, including a novel dust model for the moving mesh code {\sc Arepo}. This model includes dust production, growth, supernova-shock-driven destruction, ion-collision-driven thermal sputtering, and high temperature dust cooling through far infrared re-radiation of collisionally deposited electron energies. Adopting a rather low thermal sputtering rate, we find, consistent with observations, a present-day overall dust-to-gas ratio of $\sim 2\times 10^{-5}$, a total dust mass of $\sim 2\times 10^9\,{\rm M_\odot}$, and a dust mass fraction of $\sim 3\times 10^{-6}$. The typical thermal sputtering timescales within $\sim 100\,{\rm kpc}$ are around $\sim 10\,{\rm Myr}$, and increase towards the outer parts of the cluster to $\sim 10^3\,{\rm Myr}$ at a cluster-centric distance of $1\,{\rm Mpc}$. The condensation of gas phase metals into dust grains reduces high temperature metal-line cooling, but also leads to additional dust infrared cooling. The additional infrared cooling changes the overall cooling rate in the outer parts of the cluster, beyond $\sim 1\,{\rm Mpc}$, by factors of a few. This results in noticeable changes of the entropy, temperature, and density profiles of cluster gas once dust formation is included. The emitted dust infrared emission due to dust cooling is consistent with observational constraints.

We introduce AREPO-RT, a novel radiation hydrodynamic (RHD) solver for the unstructured moving-mesh code AREPO. Our method solves the moment-based radiative transfer equations using the M1 closure relation. We achieve second order convergence by using a slope limited linear spatial extrapolation and a first order time prediction step to obtain the values of the primitive variables on both sides of the cell interface. A Harten-Lax-Van Leer flux function, suitably modified for moving meshes, is then used to solve the Riemann problem at the interface. The implementation is fully conservative and compatible with the individual timestepping scheme of AREPO. It incorporates atomic Hydrogen (H) and Helium (He) thermochemistry, which is used to couple the ultra-violet (UV) radiation field to the gas. Additionally, infrared radiation is coupled to the gas under the assumption of local thermodynamic equilibrium between the gas and the dust. We successfully apply our code to a large number of test problems, including applications such as the expansion of ${\rm H_{II}}$ regions, radiation pressure driven outflows and the levitation of optically thick layer of gas by trapped IR radiation. The new implementation is suitable for studying various important astrophysical phenomena, such as the effect of radiative feedback in driving galactic scale outflows, radiation driven dusty winds in high redshift quasars, or simulating the reionisation history of the Universe in a self consistent manner.

Interstellar dust is an important component of the galactic ecosystem, playing a key role in multiple galaxy formation processes. We present a novel numerical framework for the dynamics and size evolution of dust grains implemented in the moving-mesh hydrodynamics code AREPO suited for cosmological galaxy formation simulations. We employ a particle-based method for dust subject to dynamical forces including drag and gravity. The drag force is implemented using a second-order semi-implicit integrator and validated using several dust-hydrodynamical test problems. Each dust particle has a grain size distribution, describing the local abundance of grains of different sizes. The grain size distribution is discretised with a second-order piecewise linear method and evolves in time according to various dust physical processes, including accretion, sputtering, shattering, and coagulation. We present a novel scheme for stochastically forming dust during stellar evolution and new methods for sub-cycling of dust physics time-steps. Using this model, we simulate an isolated disc galaxy to study the impact of dust physical processes that shape the interstellar grain size distribution. We demonstrate, for example, how dust shattering shifts the grain size distribution to smaller sizes resulting in a significant rise of radiation extinction from optical to near-ultraviolet wavelengths. Our framework for simulating dust and gas mixtures can readily be extended to account for other dynamical processes relevant in galaxy formation, like magnetohydrodynamics, radiation pressure, and thermo-chemical processes.

The fundamental metallicity relation (FMR) is a postulated correlation between galaxy stellar mass, star formation rate (SFR), and gas-phase metallicity. At its core, this relation posits that offsets from the mass-metallicity relation (MZR) at a fixed stellar mass are correlated with galactic SFR. In this Letter, we quantify the timescale with which galactic SFRs and metallicities evolve using hydrodynamical simulations. We find that Illustris and IllustrisTNG predict that galaxy offsets from the star formation main sequence and MZR evolve over similar timescales, are often anti-correlated in their evolution, evolve with the halo dynamical time, and produce a pronounced FMR. In fact, for a FMR to exist, the metallicity and SFR must evolve in an anti-correlated sense which requires that they evolve with similar time variability. In contrast to Illustris and IllustrisTNG, we speculate that the SFR and metallicity evolution tracks may become decoupled in galaxy formation models dominated by globally-bursty SFR histories, which could weaken the FMR residual correlation strength. This opens the possibility of discriminating between bursty and non-bursty feedback models based on the strength and persistence of the FMR -- especially at high redshift.

We present the photometric and spectroscopic evolution of the type Ic supernova LSQ14efd, discovered by the La Silla QUEST survey and followed by PESSTO. LSQ14efd was discovered few days after explosion and the observations cover up to ~100 days. The early photometric points show the signature of the cooling of the shock break-out event experienced by the progenitor at the time of the supernova explosion, one of the first for a type Ic supernova. A comparison with type Ic supernova spectra shows that LSQ14efd is quite similar to the type Ic SN 2004aw. These two supernovae have kinetic energies that are intermediate between standard Ic explosions and those which are the most energetic explosions known (e.g. SN 1998bw). We computed an analytical model for the light-curve peak and estimated the mass of the ejecta 6.3 +/- 0.5 Msun, a synthesized nickel mass of 0.25 Msun and a kinetic energy of Ekin = 5.6 +/- 0.5 x 10^51 erg. No connection between LSQ14efd and a GRB event could be established. However we point out that the supernova shows some spectroscopic similarities with the peculiar SN-Ia 1999ac and the SN-Iax SN 2008A. A core-collapse origin is most probable considering the spectroscopic, photometric evolution and the detection of the cooling of the shock break-out.

We present full volume cosmological simulations using the moving-mesh code AREPO to study the coevolution of dust and galaxies. We extend the dust model in AREPO to include thermal sputtering of grains and investigate the evolution of the dust mass function, the cosmic distribution of dust beyond the interstellar medium, and the dependence of dust-to-stellar mass ratio on galactic properties. The simulated dust mass function is well-described by a Schechter fit and lies closest to observations at $z = 0$. The radial scaling of projected dust surface density out to distances of $10 \, \text{Mpc}$ around galaxies with magnitudes $17 < i < 21$ is similar to that seen in Sloan Digital Sky Survey data, albeit with a lower normalisation. At $z = 0$, the predicted dust density of $\Omega_\text{dust} \approx 1.3 \times 10^{-6}$ lies in the range of $\Omega_\text{dust}$ values seen in low-redshift observations. We find that dust-to-stellar mass ratio anti-correlates with stellar mass for galaxies living along the star formation main sequence. Moreover, we estimate the $850 \, \mu\text{m}$ number density functions for simulated galaxies and analyse the relation between dust-to-stellar flux and mass ratios at $z = 0$. At high redshift, our model fails to produce enough dust-rich galaxies, and this tension is not alleviated by adopting a top-heavy initial mass function. We do not capture a decline in $\Omega_\text{dust}$ from $z = 2$ to $z = 0$, which suggests that dust production mechanisms more strongly dependent on star formation may help to produce the observed number of dusty galaxies near the peak of cosmic star formation.

The LaSilla/QUEST Variability Survey (LSQ) and the Carnegie Supernova Project (CSP II) are collaborating to discover and obtain photometric light curves for a large sample of low redshift (z < 0.1) Type Ia supernovae. The supernovae are discovered in the LSQ survey using the 1 m ESO Schmidt telescope at the La Silla Observatory with the 10 square degree QUEST camera. The follow-up photometric observations are carried out using the 1 m Swope telescope and the 2.5 m du Pont telescopes at the Las Campanas Observatory. This paper describes the survey, discusses the methods of analyzing the data and presents the light curves for the first 31 Type Ia supernovae obtained in the survey. The SALT 2.4 supernova light curve fitter was used to analyze the photometric data, and the Hubble diagram for this first sample is presented. The measurement errors for these supernovae averaged 4%, and their intrinsic spread was 14%.

We present optical imaging and spectroscopy of supernova (SN) LSQ13fn, a type II supernova with several hitherto-unseen properties. Although it initially showed strong symmetric spectral emission features attributable to \ion{He}{ii}, \ion{N}{iii}, and \ion{C}{iii}, reminiscent of some interacting SNe, it transitioned into an object that would fall more naturally under a type II-Plateau (IIP) classification. However, its spectral evolution revealed several unusual properties: metal lines appeared later than expected, were weak, and some species were conspicuous by their absence. Furthermore, the line velocities were found to be lower than expected given the plateau brightness, breaking the SNe~IIP standardised candle method for distance estimates. We found that, in combination with a short phase of early-time ejecta-circumstellar material interaction, metal-poor ejecta, and a large progenitor radius could reasonably account for the observed behaviour. Comparisons with synthetic model spectra of SNe~IIP of a given progenitor mass would imply a progenitor star metallicity as low as 0.1\,Z$_{\odot}$. LSQ13fn highlights the diversity of SNe~II and the many competing physical effects that come into play towards the final stages of massive star evolution immediately preceding core-collapse.

We introduce a dust model for cosmological simulations implemented in the moving-mesh code AREPO and present a suite of cosmological hydrodynamical zoom-in simulations to study dust formation within galactic haloes. Our model accounts for the stellar production of dust, accretion of gas-phase metals onto existing grains, destruction of dust through local supernova activity, and dust driven by winds from star-forming regions. We find that accurate stellar and active galactic nuclei feedback is needed to reproduce the observed dust-metallicity relation and that dust growth largely dominates dust destruction. Our simulations predict a dust content of the interstellar medium which is consistent with observed scaling relations at $z = 0$, including scalings between dust-to-gas ratio and metallicity, dust mass and gas mass, dust-to-gas ratio and stellar mass, and dust-to-stellar mass ratio and gas fraction. We find that roughly two-thirds of dust at $z = 0$ originated from Type II supernovae, with the contribution from asymptotic giant branch stars below 20 per cent for $z \gtrsim 5$. While our suite of Milky Way-sized galaxies forms dust in good agreement with a number of key observables, it predicts a high dust-to-metal ratio in the circumgalactic medium, which motivates a more realistic treatment of thermal sputtering of grains and dust cooling channels.

We assemble a sample of 24 hydrogen-poor super-luminous supernovae (SLSNe). Parameterizing the light curve shape through rise and decline timescales shows that the two are highly correlated. Magnetar-powered models can reproduce the correlation, with the diversity in rise and decline rates driven by the diffusion timescale. Circumstellar interaction models can exhibit a similar rise-decline relation, but only for a narrow range of densities, which may be problematic for these models. We find that SLSNe are approximately 3.5 magnitudes brighter and have light curves 3 times broader than SNe Ibc, but that the intrinsic shapes are similar. There are a number of SLSNe with particularly broad light curves, possibly indicating two progenitor channels, but statistical tests do not cleanly separate two populations. The general spectral evolution is also presented. Velocities measured from Fe II are similar for SLSNe and SNe Ibc, suggesting that diffusion time differences are dominated by mass or opacity. Flat velocity evolution in most SLSNe suggests a dense shell of ejecta. If opacities in SLSNe are similar to other SNe Ibc, the average ejected mass is higher by a factor 2-3. Assuming $\kappa=0.1\,$cm$^2\,$g$^{-1}$, we estimate a mean (median) SLSN ejecta mass of 10$\,$M$_\odot$ (6$\,$M$_\odot$), with a range of 3-30$\,$M$_\odot$. Doubling the assumed opacity brings the masses closer to normal SNe Ibc, but with a high-mass tail. The most probable mechanism for generating SLSNe seems to be the core-collapse of a very massive hydrogen-poor star, forming a millisecond magnetar.

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P.

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

We present data for LSQ14bdq, a hydrogen-poor super-luminous supernova (SLSN) discovered by the La Silla QUEST survey and classified by the Public ESO Spectroscopic Survey of Transient Objects. The spectrum and light curve are very similar to slow-declining SLSNe such as PTF12dam. However, detections within ~1 day after explosion show a bright and relatively fast initial peak, lasting for ~15 days, prior to the usual slow rise to maximum light. The broader, main peak can be fit with either central engine or circumstellar interaction models. We discuss the implications of the precursor peak in the context of these models. It is too bright and narrow to be explained as a normal 56Ni-powered SN, and we suggest that interaction models may struggle to fit the two peaks simultaneously. We propose that the initial peak may arise from the post-shock cooling of extended stellar material, and reheating by a central engine drives the second peak. In this picture, we show that an explosion energy of ~2x10^{52} erg and a progenitor radius of a few hundred solar radii would be required to power the early emission. The competing engine models involve rapidly spinning magnetars (neutron stars) or fall-back accretion onto a central black hole. The prompt energy required may favour the black hole scenario. The bright initial peak effectively rules out a compact Wolf-Rayet star as a progenitor, since the inferred energies and ejected masses become unphysical.

The Public European Southern Observatory Spectroscopic Survey of Transient Objects (PESSTO) began as a public spectroscopic survey in April 2012. We describe the data reduction strategy and data products which are publicly available through the ESO archive as the Spectroscopic Survey Data Release 1 (SSDR1). PESSTO uses the New Technology Telescope with EFOSC2 and SOFI to provide optical and NIR spectroscopy and imaging. We target supernovae and optical transients brighter than 20.5mag for classification. Science targets are then selected for follow-up based on the PESSTO science goal of extending knowledge of the extremes of the supernova population. The EFOSC2 spectra cover 3345-9995A (at resolutions of 13-18 Angs) and SOFI spectra cover 0.935-2.53 micron (resolutions 23-33 Angs) along with JHK imaging. This data release contains spectra from the first year (April 2012 - 2013), consisting of all 814 EFOSC2 spectra and 95 SOFI spectra (covering 298 distinct objects), in standard ESO Phase 3 format. We estimate the accuracy of the absolute flux calibrations for EFOSC2 to be typically 15%, and the relative flux calibration accuracy to be about 5%. The PESSTO standard NIR reduction process does not yet produce high accuracy absolute spectrophotometry but the SOFI JHK imaging will improve this. Future data releases will focus on improving the automated flux calibration of the data products.

We present optical observations of the peculiar stripped-envelope supernovae (SNe) LSQ12btw and LSQ13ccw discovered by the La Silla-QUEST survey. LSQ12btw reaches an absolute peak magnitude of M(g) = -19.3 +- 0.2, and shows an asymmetric light curve. Stringent prediscovery limits constrain its rise time to maximum light to less than 4 days, with a slower post-peak luminosity decline, similar to that experienced by the prototypical SN~Ibn 2006jc. LSQ13ccw is somewhat different: while it also exhibits a very fast rise to maximum, it reaches a fainter absolute peak magnitude (M(g) = -18.4 +- 0.2), and experiences an extremely rapid post-peak decline similar to that observed in the peculiar SN~Ib 2002bj. A stringent prediscovery limit and an early marginal detection of LSQ13ccw allow us to determine the explosion time with an uncertainty of 1 day. The spectra of LSQ12btw show the typical narrow He~I emission lines characterising Type Ibn SNe, suggesting that the SN ejecta are interacting with He-rich circumstellar material. The He I lines in the spectra of LSQ13ccw exhibit weak narrow emissions superposed on broad components. An unresolved Halpha line is also detected, suggesting a tentative Type Ibn/IIn classification. As for other SNe~Ibn, we argue that LSQ12btw and LSQ13ccw likely result from the explosions of Wolf-Rayet stars that experienced instability phases prior to core collapse. We inspect the host galaxies of SNe Ibn, and we show that all of them but one are hosted in spiral galaxies, likely in environments spanning a wide metallicity range.

In recent years, wide-field sky surveys providing deep multi-band imaging have presented a new path for indirectly characterizing the progenitor populations of core-collapse supernovae (SN): systematic light curve studies. We assemble a set of 76 grizy-band Type IIP SN light curves from Pan-STARRS1, obtained over a constant survey program of 4 years and classified using both spectroscopy and machine learning-based photometric techniques. We develop and apply a new Bayesian model for the full multi-band evolution of each light curve in the sample. We find no evidence of a sub-population of fast-declining explosions (historically referred to as "Type IIL" SNe). However, we identify a highly significant relation between the plateau phase decay rate and peak luminosity among our SNe IIP. These results argue in favor of a single parameter, likely determined by initial stellar mass, predominantly controlling the explosions of red supergiants. This relation could also be applied for supernova cosmology, offering a standardizable candle good to an intrinsic scatter of 0.2 mag. We compare each light curve to physical models from hydrodynamic simulations to estimate progenitor initial masses and other properties of the Pan-STARRS1 Type IIP SN sample. We show that correction of systematic discrepancies between modeled and observed SN IIP light curve properties and an expanded grid of progenitor properties, are needed to enable robust progenitor inferences from multi-band light curve samples of this kind. This work will serve as a pathfinder for photometric studies of core-collapse SNe to be conducted through future wide field transient searches.

We present an analysis of the early, rising light curves of 18 Type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF) and the La Silla-QUEST variability survey (LSQ). We fit these early data flux using a simple power-law $(f(t) = {\alpha\times t^n})$ to determine the time of first light $({t_0})$, and hence the rise-time $({t_{rise}})$ from first light to peak luminosity, and the exponent of the power-law rise ($n$). We find a mean uncorrected rise time of $18.98 {\pm} 0.54$ days, with individual SN rise-times ranging from $15.98$ to $24.7$ days. The exponent n shows significant departures from the simple 'fireball model' of $n = 2$ (or ${f(t) \propto t^2}$) usually assumed in the literature. With a mean value of $n = 2.44 {\pm} 0.13$, our data also show significant diversity from event to event. This deviation has implications for the distribution of 56Ni throughout the SN ejecta, with a higher index suggesting a lesser degree of 56Ni mixing. The range of n found also confirms that the 56Ni distribution is not standard throughout the population of SNe Ia, in agreement with earlier work measuring such abundances through spectral modelling. We also show that the duration of the very early light curve, before the luminosity has reached half of its maximal value, does not correlate with the light curve shape or stretch used to standardise SNe Ia in cosmological applications. This has implications for the cosmological fitting of SN Ia light curves.

We present photospheric-phase observations of LSQ12gdj, a slowly-declining, UV-bright Type Ia supernova. Classified well before maximum light, LSQ12gdj has extinction-corrected absolute magnitude $M_B = -19.8$, and pre-maximum spectroscopic evolution similar to SN 1991T and the super-Chandrasekhar-mass SN 2007if. We use ultraviolet photometry from Swift, ground-based optical photometry, and corrections from a near-infrared photometric template to construct the bolometric (1600-23800 \AA) light curve out to 45 days past $B$-band maximum light. We estimate that LSQ12gdj produced $0.96 \pm 0.07$ $M_\odot$ of $^{56}$Ni, with an ejected mass near or slightly above the Chandrasekhar mass. As much as 27% of the flux at the earliest observed phases, and 17% at maximum light, is emitted bluewards of 3300 \AA. The absence of excess luminosity at late times, the cutoff of the spectral energy distribution bluewards of 3000 \AA, and the absence of narrow line emission and strong Na I D absorption all argue against a significant contribution from ongoing shock interaction. However, up to 10% of LSQ12gdj's luminosity near maximum light could be produced by the release of trapped radiation, including kinetic energy thermalized during a brief interaction with a compact, hydrogen-poor envelope (radius $< 10^{13}$ cm) shortly after explosion; such an envelope arises generically in double-degenerate merger scenarios.

In the past decade, several rapidly-evolving transients have been discovered whose timescales and luminosities are not easily explained by traditional supernovae (SN) models. The sample size of these objects has remained small due, at least in part, to the challenge of detecting short timescale transients with traditional survey cadences. Here we present the results from a search within the Pan-STARRS1 Medium Deep Survey (PS1-MDS) for rapidly-evolving and luminous transients. We identify 10 new transients with a time above half-maximum of less than 12 days and -16.5 > M > -20 mag. This increases the number of known events in this region of SN phase space by roughly a factor of three. The median redshift of the PS1-MDS sample is z=0.275 and they all exploded in star forming galaxies. In general, the transients possess faster rise than decline timescale and blue colors at maximum light (g - r < -0.2). Best fit blackbodies reveal photospheric temperatures/radii that expand/cool with time and explosion spectra taken near maximum light are dominated by a blue continuum, consistent with a hot, optically thick, ejecta. We find it difficult to reconcile the short timescale, high peak luminosity (L > 10^43 erg/s), and lack of UV line blanketing observed in many of these transients with an explosion powered mainly by the radioactive decay of Ni-56. Rather, we find that many are consistent with either (1) cooling envelope emission from the explosion of a star with a low-mass extended envelope which ejected very little (<0.03 M_sun) radioactive material, or (2) a shock breakout within a dense, optically thick, wind surrounding the progenitor star. After calculating the detection efficiency for objects with rapid timescales in the PS1-MDS we find a volumetric rate of 4800 - 8000 events/yr/Gpc^3 (4-7% of the core-collapse SN rate at z=0.2).

We present optical spectra and light curves for three hydrogen-poor super-luminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a few days after maximum light to 100 days later shows them to be fairly typical of this class, with spectra dominated by Ca II, Mg II, Fe II and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, 56Ni decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 days after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct re-brightening at around 100d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionisation front breaking out of the ejecta.

A key tracer of the elusive progenitor systems of Type Ia supernovae (SNe Ia) is the detection of narrow blueshifted time-varying Na I D absorption lines, interpreted as evidence of circumstellar material (CSM) surrounding the progenitor system. The origin of this material is controversial, but the simplest explanation is that it results from previous mass loss in a system containing a white dwarf and a non-degenerate companion star. We present new single-epoch intermediate-resolution spectra of 17 low-redshift SNe Ia taken with XShooter on the ESO Very Large Telescope. Combining this sample with events from the literature, we confirm an excess (~20 per cent) of SNe Ia displaying blueshifted narrow Na I D absorption features compared to non-blueshifted Na I D features. The host galaxies of SNe Ia displaying blueshifted absorption profiles are skewed towards later-type galaxies, compared to SNe Ia that show no Na I D absorption, and SNe Ia displaying blueshifted narrow Na I D absorption features have broader light curves. The strength of the Na I D absorption is stronger in SNe Ia displaying blueshifted Na I D absorption features than those without blueshifted features, and the strength of the blueshifted Na I D is correlated with the B-V colour of the SN at maximum light. This strongly suggests the absorbing material is local to the SN. In the context of the progenitor systems of SNe Ia, we discuss the significance of these findings and other recent observational evidence on the nature of SN Ia progenitors. We present a summary that suggests there are at least two distinct populations of normal, cosmologically useful SNe Ia.

We describe the La Silla-QUEST (LSQ) Variability Survey. LSQ is a dedicated wide-field synoptic survey in the Southern Hemisphere, focussing on the discovery and study of transients ranging from low redshift (z < 0.1) SN Ia, Tidal Disruption events, RR Lyr{\ae} variables, CVs, Quasars, TNOs and others. The survey utilizes the 1.0-m Schmidt Telescope of the European Southern Observatory at La Silla, Chile, with the large-area QUEST camera, a mosaic of 112 CCDs with field of view of 9.6 square degrees. The LSQ Survey was commissioned in 2009, and is now regularly covering ~1000 square deg per night with a repeat cadence of hours to days. The data are currently processed on a daily basis. We present here a first look at the photometric capabilities of LSQ and we discuss some of the most interesting recent transient detections.