Late-time observations of Type Ia supernovae (SNe Ia), >900 days after explosion, have shown that this type of SN does not suffer an "IR catastrophe" at 500 days as previously predicted. Instead, several groups have observed a slow-down in the optical light curves of these SNe. A few reasons have been suggested for this slow-down, from a changing fraction of positrons reprocessed by the expanding ejecta, through a boost of energy from slow radioactive decay chains such as 57Co --> 57Fe, to atomic "freeze-out." Discovering which of these (or some other) heating mechanisms is behind the slow-down will directly impact studies of SN Ia progenitors, explosion models, and nebular-stage physics. Recently, Graur et al. (2018) suggested a possible correlation between the shape of the late-time light curves of four SNe Ia and their stretch values, which are proxies for their intrinsic luminosities. Here, we present Hubble Space Telescope observations of the SN Ia ASASSN-14lp at ~850-960 days past maximum light. With a stretch of s = 1.15 +/- 0.05, it is the most luminous normal SN Ia observed so far at these late times. We rule out contamination by light echoes and show that the late-time, optical light curve of ASASSN-14lp is flatter than that of previous SNe Ia observed at late times. This result is in line with-and strengthens-the Graur et al. (2018) correlation, but additional SNe are needed to verify it.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We present a multi-wavelength polarimetric and spectral study of M87 jet obtained at sub- arcsecond resolution between 2002 and 2008. The observations include multi-band archival VLA polarimetry data sets along with the HST imaging polarimetry. These observations have better angular resolution than previous work by factors of 2-3 and in addition, allow us to explore the time domain. These observations envelope the huge flare in HST-1 located at 0."86 from the nucleus (Cheung et al. 2007; Harris et al. 2009; Madrid 2009; Perlman et al. 2011). The increased resolution enables us to view more structure in each knot, showing several resolved sub-components. We also see apparent helical structure in the polarization vectors in several knots, with polarization vectors turning either clockwise or counterclockwise near the flux maxima in various places as well as show filamentary undulations. Some of these characteristics are correlated with flux and polarization maxima while others are not. We also examine the total flux and fractional polarization and look for changes in both radio and optical since the observations of Perlman et al. (1999) and test them against various models based on shocks and instabilities in the jet. Our results are broadly consistent with previous spine-sheath models and recollimation shock models, however, they require additional combinations of features to explain the observed complexity, e.g. shearing of magnetic field lines near the jet surface and compression of the toroidal component near shocks. In particular, in many regions we find apparently helical features both in total flux and polarization. We discuss the physical interpretation of these features.

Relativistic jets are the most energetic manifestation of the active galactic nucleus (AGN) phenomenon. AGN jets are observed from the radio through gamma-rays and carry copious amounts of matter and energy from the sub-parsec central regions out to the kiloparsec and often megaparsec scale galaxy and cluster environs. While most spatially resolved jets are seen in the radio, an increasing number have been discovered to emit in the optical/near-IR and/or X-ray bands. Here we discuss a spectacular example of this class, the 3C 111 jet, housed in one of the nearest, double-lobed FR II radio galaxies known. We discuss new, deep Chandra and HST observations that reveal both near-IR and X-ray emission from several components of the 3C 111 jet, as well as both the northern and southern hotspots. Important differences are seen between the morphologies in the radio, X-ray and near-IR bands. The long (over 100 kpc on each side), straight nature of this jet makes it an excellent prototype for future, deep observations, as it is one of the longest such features seen in the radio, near-IR/optical and X-ray bands. Several independent lines of evidence, including the X-ray and broadband spectral shape as well as the implied velocity of the approaching hotspot, lead us to strongly disfavor the EC/CMB model and instead favor a two-component synchrotron model to explain the observed X-ray emission for several jet components. Future observations with NuSTAR, HST, and Chandra will allow us to further constrain the emission mechanisms.

We present a high resolution polarimetry and variability study of the M87 jet using VLA and HST data taken during 2002 to 2008. Both data-sets have an angular resolution as high as 0.06$"$, which is 2-3 times better than previous observations. New morphological details are revealed in both the optical and radio, which can help to reveal the energetic and magnetic field structure of the jet. By comparing the data with previously published HST and VLA observations, we show that the jet$'$s morphology in total and polarized light is changing significantly on timescales of $\sim$a decade. We compare the evolution of the inner jet (particularly the nucleus and knot HST-1), when our observations overlap with the multi-wavelength monitoring campaigns conducted with HST and Chandra. We use these data to comment on particle acceleration and main emission processes.

The DrizzlePac package includes tasks for aligning and drizzling images taken with the Hubble Space Telescope. We present this release which includes new features that facilitate image alignment, sky matching, and adds support for new time dependent distortion solutions of the ACS instrument. The TweakReg task now includes capabilities for automatically aligning images which form part of a mosaic. In addition, new parameters make it easier to reject cosmic rays and other spurious detections from source catalogs used for alignment. The Astrodrizzle task has been improved with a new sky matching algorithm which makes producing mosaics easier than ever before. This new version supports an improved version of the ACS/WFC time-dependent distortion correction. There are also improvements to the GUI interfaces and some behind the scene bug fixes.

Since the discovery of kiloparsec-scale X-ray emission from quasar jets, the physical processes responsible for their high-energy emission have been poorly defined. A number of mechanisms are under active debate, including synchrotron radiation, inverse-Comptonized CMB (IC/CMB) emission, and other Comptonization processes. In a number of cases, the optical and X-ray emission of jet regions are inked by a single spectral component, and in those, high- resolution multi-band imaging and polarimetry can be combined to yield a powerful diagnostic of jet emission processes. Here we report on deep imaging photometry of the jet of PKS 1136$-$135 obtained with the {\it Hubble Space Telescope.} We find that several knots are highly polarized in the optical, with fractional polarization $\Pi>30%$. When combined with the broadband spectral shape observed in these regions, this is very difficult to explain via IC/CMB models, unless the scattering particles are at the lowest-energy tip of the electron energy distribution, with Lorentz factor $\gamma \sim 1$, and the jet is also very highly beamed ($\delta \geq 20$) and viewed within a few degrees of the line of sight. We discuss both the IC/CMB and synchrotron interpretation of the X-ray emission in the light of this new evidence, presenting new models of the spectral energy distribution and also the matter content of this jet. The high polarizations do not completely rule out the possibility of IC/CMB optical-to-X-ray emission in this jet, but they do strongly disfavor the model. We discuss the implications of this finding, and also the prospects for future work.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

During the last decade, M87's jet has been the site of an extraordinary variability event, with one knot (HST-1) increasing by over a factor 100 in brightness. Variability was also seen on timescales of months in the nuclear flux. Here we discuss the optical-UV polarization and spectral variability of these components, which show vastly different behavior. HST-1 shows a highly significant correlation between flux and polarization, with P increasing from $\sim 20%$ at minimum to >40% at maximum, while the orientation of its electric vector stayed constant. HST-1's optical-UV spectrum is very hard ($\alpha_{UV-O}\sim0.5$, $F_\nu\propto\nu^{-\alpha}$), and displays "hard lags" during epochs 2004.9-2005.5, including the peak of the flare, with soft lags at later epochs. We interpret the behavior of HST-1 as enhanced particle acceleration in a shock, with cooling from both particle aging and the relaxation of the compression. We set 2$\sigma$ upper limits of $0.5 \delta$ parsecs and 1.02$c$ on the size and advance speed of the flaring region. The slight deviation of the electric vector orientation from the jet PA, makes it likely that on smaller scales the flaring region has either a double or twisted structure. By contrast, the nucleus displays much more rapid variability, with a highly variable electric vector orientation and 'looping' in the $(I,P)$ plane. The nucleus has a much steeper spectrum ($\alpha_{UV-O} \sim 1.5$) but does not show UV-optical spectral variability. Its behavior can be interpreted as either a helical distortion to a steady jet or a shock propagating through a helical jet.

The emission mechanisms in extragalactic jets include synchrotron and various inverse-Compton processes. At low (radio through infrared) energies, it is widely agreed that synchrotron emission dominates in both low-power (FR I) and high-power (FR II and quasar) jets, because of the power-law nature of the spectra observed and high polarizations. However, at higher energies, the emission mechanism for high-power jets at kpc scales is hotly debated. Two mechanisms have been proposed: either inverse-Compton of cosmic microwave background photons or synchrotron emission from a second, high-energy population of electrons. Here we discuss optical polarimetry as a method for diagnosing the mechanism for the high-energy emission in quasar jets, as well as revealing the jet's three-dimensional energetic and magnetic field structure. We then discuss high-energy emission mechanisms for powerful jets in the light of the HST polarimetry of PKS 1136-135.

We show that using either the method of Page & Carrera or the well-known $1/V_a$ method to construct the binned luminosity function (LF) of a flux limited sample of Active Galactic Nuclei (AGN) can produce an artificial flattening (or steepening in the case of negative evolution) of the binned LF for bins intersected by the flux cutoff of the sample. This effect is more pronounced for samples with steep and strongly evolving parent LFs but is still present even for non-evolving LFs. As a result of this distortion of the true LF, fitting a model LF to binned data may lead to errors in the estimation of the parameters and may even prompt the erroneous use of broken power law functions. We compute the expected positions of apparent breaks in the binned LF. We show that these spurious breaks in the binned LFs can be avoided if the binning is done in the flux--redshift plane instead of the typically used luminosity--redshift plane. Binning in the flux--redshift plane can be used in conjunction with the binning in the luminosity--redshift plane to test for real breaks in the binned LFs and to identify the features that are the result of binning biases. We illustrate this effect for most typical forms of luminosity dependence and redshift evolution and show how the proposed method helps address this problem. We also apply this method to the MOJAVE AGN sample and show that it eliminates an apparent break in the binned LF.

(Abridged) We use a complete sample of active galactic nuclei (AGN) selected on the basis of relativistically beamed 15 GHz radio flux density to derive the parent radio luminosity function (RLF) of bright radio-selected blazar cores. We use a maximum likelihood method to fit a beamed RLF to the observed data and thereby recover the parameters of the intrinsic (unbeamed) RLF. We analyze two subsamples of the MOJAVE sample: the first contains only objects of known FR II class, with a total of 103 sources, and the second subsample adds 24 objects of uncertain FR class for a total of 127 sources. Both subsamples exclude four known FR I radio galaxies and two gigahertz-peaked spectrum sources. We obtain good fits to both subsamples using a single power law intrinsic RLF with pure density evolution function. We find that a previously reported break in the observed MOJAVE RLF actually arises from using incomplete bins (because of the luminosity cutoff) across a steep and strongly evolving RLF, and does not reflect a break in the intrinsic RLF. The derived space density of the parent population of the FR II sources from the MOJAVE sample (with L>1.3e25 W/Hz) is approximately 1600/Gpc^3.