We present a single-dish mapping algorithm with a number of advantages over traditional techniques. (1) Our algorithm makes use of weighted modeling, instead of weighted averaging, to interpolate between signal measurements. This smooths the data, but without blurring the data beyond instrumental resolution. Techniques that rely on weighted averaging blur point sources sometimes as much as 40%. (2) Our algorithm makes use of local, instead of global, modeling to separate astronomical signal from instrumental and/or environmental signal drift along the telescope's scans. Other techniques, such as basket weaving, model this drift with simple functional forms (linear, quadratic, etc.) across the entirety of scans, limiting their ability to remove such contaminants. (3) Our algorithm makes use of a similar, local modeling technique to separate astronomical signal from radio-frequency interference (RFI), even if only continuum data are available. (4) Unlike other techniques, our algorithm does not require data to be collected on a rectangular grid or regridded before processing. (5) Data from any number of observations, overlapping or not, may be appended and processed together. (6) Any pixel density may be selected for the final image. We present our algorithm, and evaluate it using both simulated and real data. We are integrating it into the image-processing library of the Skynet Robotic Telescope Network, which includes optical telescopes spanning four continents, and now also Green Bank Observatory's 20-meter diameter radio telescope in West Virginia. Skynet serves hundreds of professional users, and additionally tens of thousands of students, of all ages. Default data products are generated on the fly, but will soon be customizable after the fact.

Sigma clipping is commonly used in astronomy for outlier rejection, but the number of standard deviations beyond which one should clip data from a sample ultimately depends on the size of the sample. Chauvenet rejection is one of the oldest, and simplest, ways to account for this, but, like sigma clipping, depends on the sample's mean and standard deviation, neither of which are robust quantities: Both are easily contaminated by the very outliers they are being used to reject. Many, more robust measures of central tendency, and of sample deviation, exist, but each has a tradeoff with precision. Here, we demonstrate that outlier rejection can be both very robust and very precise if decreasingly robust but increasingly precise techniques are applied in sequence. To this end, we present a variation on Chauvenet rejection that we call "robust" Chauvenet rejection (RCR), which uses three decreasingly robust/increasingly precise measures of central tendency, and four decreasingly robust/increasingly precise measures of sample deviation. We show this sequential approach to be very effective for a wide variety of contaminant types, even when a significant -- even dominant -- fraction of the sample is contaminated, and especially when the contaminants are strong. Furthermore, we have developed a bulk-rejection variant, to significantly decrease computing times, and RCR can be applied both to weighted data, and when fitting parameterized models to data. We present aperture photometry in a contaminated, crowded field as an example. RCR may be used by anyone at https://skynet.unc.edu/rcr, and source code is available there as well.

We present a comprehensive temporal and spectral analysis of the long Swift GRB 120327A afterglow data to investigate the possible causes of the observed early time colour variations. We collected data from various instruments/telescopes in different bands (X-rays, ultra- violet, optical and near-infrared) and determined the shapes of the afterglow early-time light curves. We studied the overall temporal behaviour and the spectral energy distributions from early to late times. The ultra-violet, optical, and near-infrared light curves can be modelled with a single power-law component between 200 and 2e4 s after the burst event. The X-ray light curve shows a canonical steep-shallow-steep behaviour, typical of long gamma-ray bursts. At early times a colour variation is observed in the ultra-violet/optical bands, while at very late times a hint of a re-brightening is visible. The observed early time colour change can be explained as a variation in the intrinsic optical spectral index, rather than an evolution of the optical extinction.

Based on five years of observations with the 40-foot telescope at Green Bank Observatory (GBO), Reichart & Stephens (2000) found that the radio source Cassiopeia A had either faded more slowly between the mid-1970s and late 1990s than Baars et al. (1977) had found it to be fading between the late 1940s and mid-1970s, or that it had rebrightened and then resumed fading sometime between the mid-1970s and mid-1990s, in L band (1.4 GHz). Here, we present 15 additional years of observations of Cas A and Cyg A with the 40-foot in L band, and three and a half additional years of observations of Cas A, Cyg A, Tau A, and Vir A with GBO's recently refurbished 20-meter telescope in L and X (9 GHz) bands. We also present a more sophisticated analysis of the 40-foot data, and a reanalysis of the Baars et al. (1977) data, which reveals small, but non-negligible differences. We find that overall, between the late 1950s and late 2010s, Cas A faded at an average rate of $0.670 \pm 0.019$ %/yr in L band, consistent with Reichart & Stephens (2000). However, we also find, at the 6.3$\sigma$ credible level, that it did not fade at a constant rate. Rather, Cas A faded at a faster rate through at least the late 1960s, rebrightened (or at least faded at a much slower rate), and then resumed fading at a similarly fast rate by, at most, the late 1990s. Given these differences from the original Baars et al. (1977) analysis, and given the importance of their fitted spectral and temporal models for flux-density calibration in radio astronomy, we update and improve on these models for all four of these radio sources. In doing so, we additionally find that Tau A is fading at a rate of $0.102^{+0.042}_{-0.043}$ %/yr in L band.

GRB 141221A was observed from infrared to soft gamma-ray bands. Here, we investigate its properties, in light of the standard model. We find that the optical light curve of the afterglow of this burst presents an unusual steep/quick rise. The broad band spectral energy distribution taken near the maximum of the optical emission presents either a thermal component or a spectral break. In the former case, the properties of the afterglow are then very unusual, but could explain the lack of apparent jet breaks in the Swift light curves. In the latter case, the afterglow properties of this burst are more usual, and we can see in the light curves the passing through of the injection and cooling frequencies within the optical bands, not masked by a reverse shock. This model also excludes the presence of a stellar wind, challenging either the stellar progenitor properties, or the very stellar nature of the progenitor itself. In all cases, this burst may be a part of a Rosetta stone that could help to explain some of the most striking features discovered by Swift during the last ten years.

Gamma-ray bursts (GRBs) are the most luminous explosions in the universe, yet the nature and physical properties of their energy sources are far from understood. Very important clues, however, can be inferred by studying the afterglows of these events. We present optical and X-ray observations of GRB 130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and GTC. This burst shows a steep drop in the X-ray light-curve at $\simeq 10^5$ s after the trigger, with a power-law decay index of $\alpha \sim 6$. Such a rare behaviour cannot be explained by the standard forward shock (FS) model and indicates that the emission, up to the fast decay at $10^5$ s, must be of "internal origin", produced by a dissipation process within an ultrarelativistic outflow. We propose that the source of such an outflow, which must produce the X-ray flux for $\simeq 1$ day in the cosmological rest frame, is a newly born magnetar or black hole. After the drop, the faint X-ray afterglow continues with a much shallower decay. The optical emission, on the other hand, shows no break across the X-ray steep decrease, and the late-time decays of both the X-ray and optical are consistent. Using both the X-ray and optical data, we show that the emission after $\simeq 10^5$ s can be explained well by the FS model. We model our data to derive the kinetic energy of the ejecta and thus measure the efficiency of the central engine of a GRB with emission of internal origin visible for a long time. Furthermore, we break down the energy budget of this GRB into the prompt emission, the late internal dissipation, the kinetic energy of the relativistic ejecta, and compare it with the energy of the associated supernova, SN 2013fu.

We present the first reported case of the simultaneous metallicity determination of a gamma-ray burst (GRB) host galaxy, from both afterglow absorption lines as well as strong emission-line diagnostics. Using spectroscopic and imaging observations of the afterglow and host of the long-duration Swift GRB121024A at z = 2.30, we give one of the most complete views of a GRB host/environment to date. We observe a strong damped Ly-alpha absorber (DLA) with a hydrogen column density of log N(HI) = 21.88 +/- 0.10, H2 absorption in the Lyman-Werner bands (molecular fraction of log(f)~ -1.4; fourth solid detection of molecular hydrogen in a GRB-DLA), the nebular emission lines H-alpha, H-beta, [O II], [O III] and [N II], as well as metal absorption lines. We find a GRB host galaxy that is highly star-forming (SFR ~ 40 solar masses/yr ), with a dust-corrected metallicity along the line of sight of [Zn/H]corr = -0.6 +/- 0.2 ([O/H] ~ -0.3 from emission lines), and a depletion factor [Zn/Fe] = 0.85 +/- 0.04. The molecular gas is separated by 400 km/s (and 1-3 kpc) from the gas that is photoexcited by the GRB. This implies a fairly massive host, in agreement with the derived stellar mass of log(M/M_solar ) = 9.9+/- 0.2. We dissect the host galaxy by characterising its molecular component, the excited gas, and the line-emitting star-forming regions. The extinction curve for the line of sight is found to be unusually flat (Rv ~15). We discuss the possibility of an anomalous grain size distributions. We furthermore discuss the different metallicity determinations from both absorption and emission lines, which gives consistent results for the line of sight to GRB 121024A.

A long and intense gamma-ray burst (GRB) was detected by INTEGRAL on July 11 2012 with a duration of ~115s and fluence of 2.8x10^-4 erg cm^-2 in the 20 keV-8 MeV energy range. GRB 120711A was at z~1.405 and produced soft gamma-ray emission (>20 keV) for at least ~10 ks after the trigger. The GRB was observed by several ground-based telescopes that detected a powerful optical flash peaking at an R-band brightness of ~11.5 mag at ~126 s after the trigger. We present a comprehensive temporal and spectral analysis of the long-lasting soft gamma-ray emission detected in the 20-200 keV band with INTEGRAL, the Fermi/LAT post-GRB detection above 100 MeV, the soft X-ray afterglow from XMM-Newton, Chandra, and Swift and the optical/NIR detections from Watcher, Skynet, GROND, and REM. We modelled the long-lasting soft gamma-ray emission using the standard afterglow scenario, which indicates a forward shock origin. The combination of data extending from the NIR to GeV energies suggest that the emission is produced by a broken power-law spectrum consistent with synchrotron radiation. The afterglow is well modelled using a stratified wind-like environment with a density profile k~1.2, suggesting a massive star progenitor (i.e. Wolf-Rayet). The analysis of the reverse and forward shock emission reveals an initial Lorentz factor of ~120-340, a jet half-opening angle of ~2deg-5deg, and a baryon load of ~10^-5-10^-6 Msun consistent with the expectations of the fireball model when the emission is highly relativistic. Long-lasting soft gamma-ray emission from other INTEGRAL GRBs with high peak fluxes, such as GRB 041219A, was not detected, suggesting that a combination of high Lorentz factor, emission above 100 MeV, and possibly a powerful reverse shock are required. Similar long-lasting soft gamma-ray emission has recently been observed from the nearby and extremely bright Fermi/LAT burst GRB 130427A.

PROMPT (Panchromatic Robotic Optical Monitoring and Polarimetry Telescopes) observed the early-time optical afterglow of GRB 060607A and obtained a densely sampled multiwavelength light curve that begins only tens of seconds after the GRB. Located at Cerro Tololo Inter-American Observatory in Chile, PROMPT is designed to observe the afterglows of gamma-ray bursts using multiple automated 0.4-m telescopes that image simultaneously in many filters when the afterglow is bright and may be highly variable. The data span the interval from 44 seconds after the GRB trigger to 3.3 hours in the Bgri filters. We observe an initial peak in the light curve at approximately three minutes, followed by rebrightenings peaking around 40 minutes and again at 66 minutes. Although our data overlap with the early Swift gamma-ray and x-ray light curves, we do not see a correlation between the optical and high-energy flares. We do not find evidence for spectral evolution throughout the observations. We model the variations in the light curves and find that the most likely cause of the rebrightening episodes is a refreshment of the forward shock preceded by a rapidly fading reverse shock component, although other explanations are plausible.

A maximum likelihood analysis of the NGC 4258 maser positions and velocities reveals a ~2 sigma deviation from Keplerian motion in the projected rotation curve of the high-velocity features, corresponding to a ~9 km/s, or 0.8%, flattening of the LOS velocities with respect to Keplerian motion over the range of the high-velocity masers. While there are a number of potential explanations for this flattening, we argue for pure Keplerian rotation in an inclination-warped disk based on the ability of this model to explain a number of otherwise puzzling features of the system. A radial gradient in the disk inclination of 0.034 mas^{-1} is not only consistent with the observed rotation curve, but it generates a bowl along the near edge of the disk that naturally explains the otherwise puzzling narrow spread in the declinations of the systemic masers. It also explains the existence and location of an apparently recurring flare amongst the systemic masers. There is no significant evidence for non-Keplerian rotation in the inclination-warped disk. An additional implication of the inclination warp is that the disk rises in front of and obscures the central engine at a disk radius of about 8.3 mas, or 0.29 pc. By comparing the observed X-ray column to conditions in the disk at this radius, we argue that the molecular-to-atomic transition occurs just beyond the outermost maser at 0.28 pc and from this we infer an accretion rate of \~10^{-4}alpha M_sun/yr, consistent with the jet-dominated accretion models of Yuan et al. (2002). Our model suggests that most of the observed X-ray column arises in the warped accretion disk at 0.29 pc, and that the maser emission is truncated at large radii predominantly as a result of the molecular-to-atomic phase transition originally proposed by Neufeld & Maloney (1995).

We have conducted VLBI observations of water maser and radio continuum emission in the nucleus of NGC 3079. The 22 GHz maser emission arises in compact clumps, distributed along an axis that is aligned with the major axis of the galactic disk. The velocities of the masers are consistent with their lying in the inner parsec of a molecular disk rotating in the same sense as the rest of the galaxy. However, the velocity field has a significant non-rotational component, which may indicate supersonic turbulence. The bright maser emission is not coincident with any detected compact 22 GHz continuum source, suggesting the high apparent luminosity of the maser may not due to beamed amplification of continuum emission. We observed two compact continuum sources that have inverted spectra between 5 and 8 GHz, and steep spectra between 8 and 22 GHz. NGC 3079 may be a nearby, low-luminosity example of the class of compact symmetric gigahertz-peaked spectrum (GPS) radio sources. We detected a third continuum component that lies along the same axis as the other two, strongly suggesting that this galaxy possesses a nuclear jet. Faint maser emission was detected near this axis, which may indicate a second population of masers associated with the jet.

We present a parsec-scale image of the OH maser in the nucleus of the active galaxy IIIZw35, made using the Very Long Baseline Array at a wavelength of 18 cm. We detected two distinct components, with a projected separation of 50 pc (for D=110 Mpc) and a separation in Doppler velocity of 70 km/s, which contain 50% of the total maser flux. Velocity gradients within these components could indicate rotation of clouds with binding mass densities of ~7000 solar masses per cubic parsec, or total masses of more than 500,000 solar masses. Emission in the 1665-MHz OH line is roughly coincident in position with that in the 1667-MHz line, although the lines peak at different Doppler velocities. We detected no 18 cm continuum emission; our upper limit implies a peak apparent optical depth greater than 3.4, assuming the maser is an unsaturated amplifier of continuum radiation.