The NSF-sponsored Undergraduate ALFALFA Team (UAT) promotes long-term collaborative research opportunities for faculty and students from 23 U.S. public and private primarily undergraduate institutions (PUIs) within the context of the extragalactic ALFALFA HI blind legacy survey project. Over twelve project years of partnering with Arecibo and Green Bank Observatories, the UAT has had a demonstrable impact on the health of a legacy astronomy project, science education, and equity/inclusion in astronomy, with successful outcomes for 373 UAT students (39% women; ~30% members of underrepresented groups) and 34 faculty (44% women). The UAT model is adaptable to many large scientific projects and can be supported by relatively modest funding. We recommend that granting agencies identify funding resources to support the model, either as an add-on to legacy grant support or as a stand-alone funding source. This could include encouragement of UAT-like components in large scale projects currently being developed, such as the LSST and TMT. By doing this, we will recognize the high numbers of astronomy research-trained heavy-teaching-load faculty at PUIs as an under-utilized resource of the astronomy community (see also White Paper by Ribaudo et al.). These members of our community have the skills and the strong desire to contribute meaningfully to their field, as well as the ability to encourage and interact closely with many talented and motivated undergraduate students from all backgrounds.

The quasar 3C454.3 underwent a uniquely-structured multi-frequency outburst in June 2016. The blazar was observed in the optical $R$ band by several ground-based telescopes in photometric and polarimetric modes, at $\gamma$-ray frequencies by the \emph{Fermi}\ Large Area Telescope, and at 43 GHz with the Very Long Baseline Array. The maximum flux density was observed on 2016 June 24 at both optical and $\gamma$-ray frequencies, reaching $S^\mathrm{max}_\mathrm{opt}=18.91\pm0.08$ mJy and $S_\gamma^\mathrm{max} =22.20\pm0.18\times10^{-6}$ ph cm$^{-2}$ s$^{-1}$, respectively. The June 2016 outburst possessed a precipitous decay at both $\gamma$-ray and optical frequencies, with the source decreasing in flux density by a factor of 4 over a 24-hour period in $R$ band. Intraday variability was observed throughout the outburst, with flux density changes between 1 and 5 mJy over the course of a night. The precipitous decay featured statistically significant quasi-periodic micro-variability oscillations with an amplitude of $\sim 2$-$3\%$ about the mean trend and a characteristic period of 36 minutes. The optical degree of polarization jumped from $\sim3\%$ to nearly 20\% during the outburst, while the position angle varied by $\sim120\degr$. A knot was ejected from the 43 GHz core on 2016 Feb 25, moving at an apparent speed $v_\mathrm{app}=20.3c\pm0.8c$. From the observed minimum timescale of variability $\tau_\mathrm{opt}^\mathrm{min}\approx2$ hr and derived Doppler factor $\delta=22.6$, we find a size of the emission region $r\lesssim2.6\times10^{15}$ cm. If the quasi-periodic micro-variability oscillations are caused by periodic variations of the Doppler factor of emission from a turbulent vortex, we derive a rotational speed of the vortex $\sim0.2c$.

We report a multi-objective campaign of targeted 21 cm HI line observations of sources selected from the Arecibo Legacy Fast ALFA (Arecibo L-band Feed Array) survey (ALFALFA) and galaxies identified by their morphological and photometric properties in the Sloan Digital Sky Survey (SDSS). The aims of this program have been (1) to confirm the reality of some ALFALFA sources whose enigmatic nature suggest additional multiwavelength observations; (2) to probe the low signal-to-noise ratio regime, below the ALFALFA reliability limit; and (3) to explore the feasibility of using optical morphology, color and surface brightness to identify gas-rich objects in the region of the Pisces-Perseus Supercluster (PPS) whose HI fluxes are below the ALFALFA sensitivity limit at that distance. As expected, the reliability of ALFALFA detections depends strongly on the signal-to-noise ratio of the HI line signal and its coincidence with a probable stellar counterpart identified by its optical properties, suggestive of on-going star formation. The identification of low mass, star forming populations enables targeted HI line observations to detect galaxies with HI line fluxes below the ALFALFA sensitivity limits in fixed local volumes (D < 100 Mpc). The method explored here serves as the basis for extending the sample of gas-bearing objects as part of the on-going Arecibo Pisces-Perseus Supercluster Survey (APPSS).

We present the catalog of ~31500 extragalactic HI line sources detected by the completed ALFALFA survey out to z < 0.06 including both high signal-to-noise ratio (> 6.5) detections and ones of lower quality which coincide in both position and recessional velocity with galaxies of known redshift. We review the observing technique, data reduction pipeline, and catalog construction process, focusing on details of particular relevance to understanding the catalog's compiled parameters. We further describe and make available the digital HI line spectra associated with the catalogued sources. In addition to the extragalactic HI line detections, we report nine confirmed OH megamasers and ten OH megamaser candidates at 0.16 < z < 0.22 whose OH line signals are redshifted into the ALFALFA frequency band. Because of complexities in data collection and processing associated with the use of a feed-horn array on a complex single-dish antenna in the terrestrial radio frequency interference environment, we also present a list of suggestions and caveats for consideration by users of the ALFALFA extragalactic catalog for future scientific investigations.

We carried out a pilot campaign of radio and optical band intra-day variability (IDV) observations of five blazars (3C66A, S5 0716+714, OJ287, B0925+504, and BL Lacertae) on December 18--21, 2015 by using the radio telescope in Effelsberg (Germany) and several optical telescopes in Asia, Europe, and America. After calibration, the light curves from both 5 GHz radio band and the optical R band were obtained, although the data were not smoothly sampled over the sampling period of about four days. We tentatively analyse the amplitudes and time scales of the variabilities, and any possible periodicity. The blazars vary significantly in the radio (except 3C66A and BL Lacertae with only marginal variations) and optical bands on intra- and inter-day time scales, and the source B0925+504 exhibits a strong quasi-periodic radio variability. No significant correlation between the radio- and optical-band variability appears in the five sources, which we attribute to the radio IDV being dominated by interstellar scintillation whereas the optical variability comes from the source itself. However, the radio- and optical-band variations appear to be weakly correlated in some sources and should be investigated based on well-sampled data from future observations.

We present a current catalog of 21 cm HI line sources extracted from the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey over ~2800 square degrees of sky: the alpha.40 catalog. Covering 40% of the final survey area, the alpha.40 catalog contains 15855 sources in the regions 07h30m < R.A. < 16h30m, +04 deg < Dec. < +16 deg and +24 deg < Dec. < +28 deg and 22h < R.A. < 03h, +14 deg < Dec. < +16 deg and +24 deg < Dec. < +32 deg. Of those, 15041 are certainly extragalactic, yielding a source density of 5.3 galaxies per square degree, a factor of 29 improvement over the catalog extracted from the HI Parkes All Sky Survey. In addition to the source centroid positions, HI line flux densities, recessional velocities and line widths, the catalog includes the coordinates of the most probable optical counterpart of each HI line detection, and a separate compilation provides a crossmatch to identifications given in the photometric and spectroscopic catalogs associated with the Sloan Digital Sky Survey Data Release 7. Fewer than 2% of the extragalactic HI line sources cannot be identified with a feasible optical counterpart; some of those may be rare OH megamasers at 0.16 < z < 0.25. A detailed analysis is presented of the completeness, width dependent sensitivity function and bias inherent in the current alpha.40 catalog. The impact of survey selection, distance errors, current volume coverage and local large scale structure on the derivation of the HI mass function is assessed. While alpha.40 does not yet provide a completely representative sampling of cosmological volume, derivations of the HI mass function using future data releases from ALFALFA will further improve both statistical and systematic uncertainties.

We monitored the flaring blazar 3C 454.3 during 2005 June-July with the Spitzer Infrared Spectrograph (IRS: 15 epochs), Infrared Array Camera (IRAC: 12 epochs) and Multiband Imaging Photometer (MIPS: 2 epochs). We also made Spitzer IRS, IRAC, and MIPS observations from 2006 December-2007 January when the source was in a low state, the latter simultaneous with a single Chandra X-ray observation. In addition, we present optical and sub-mm monitoring data. The 2005-2007 period saw 3 major outbursts. We present evidence that the radio-optical SED actually consists of two variable synchrotron peaks, the primary at IR and the secondary at sub-mm wavelengths. The lag between the optical and sub-mm outbursts may indicate that these two peaks arise from two distinct regions along the jet separated by a distance of 0.07-5 pc. The flux at 5-35 microns varied by a factor of 40 and the IR peak varied in frequency from <1E13 Hz to 4E13 Hz between the highest and lowest states in 2005 and 2006, respectively. Variability was well correlated across the mid-IR band, with no measurable lag. Flares that doubled in flux occurred on a time scale of 3 days. The IR SED peak moved to higher frequency as a flare brightened, then returned to lower frequency as it decayed. The fractional variability amplitude increased with frequency, which we attribute to decreasing synchrotron-self absorption optical depth. Mid-IR flares may signal the re-energization of a shock that runs into inhomogeneities along the pre-existing jet or in the external medium. The synchrotron peak frequencies during each major outburst may depend upon both the distance from the jet apex and the physical conditions in the shocks. Variation of the Doppler parameter along a curved or helical jet is another possibility. Frequency variability of the IR synchrotron peak may have important consequences for the interpretation of the blazar sequence, and the presence of a secondary peak may give insight into jet structure.

We present analysis of both the short-term optical and long-term multiwavelength variability of CTA 102. In 2004, this object was observed in an intense optical flaring state. Extensive R-band microvariability observations were carried out during this high state. In 2005, we obtained several weeks of contemporaneous radio, optical, and X-ray observations of CTA 102. These observations recorded distinct flaring activity in all three wavebands. Subsequent analysis revealed that this object may appear redder when in a brighter optical state, and that the X-ray, optical, and radio activity do not appear to be correlated. The shape of the observed spectral energy distributions suggests that both synchrotron-related and external inverse Compton processes may contribute to the X-ray emission. Our results are also compared with other results on this object and archival microvariability observations. It appears that more rapid, dramatic microvariability events occur when CTA 102 is in an elevated optical flux state.

We present the results of extensive multi-waveband monitoring of the blazar 3C~279 between 1996 and 2007 at X-ray energies (2-10 keV), optical R band, and 14.5 GHz, as well as imaging with the Very Long Baseline Array (VLBA) at 43 GHz. In all bands the power spectral density corresponds to "red noise" that can be fit by a single power law over the sampled time scales. Variations in flux at all three wavebands are significantly correlated. The time delay between high and low frequency bands changes substantially on time scales of years. A major multi-frequency flare in 2001 coincided with a swing of the jet toward a more southerly direction, and in general the X-ray flux is modulated by changes in the position angle of the jet near the core. The flux density in the core at 43 GHz--increases in which indicate the appearance of new superluminal knots--is significantly correlated with the X-ray flux. We decompose the X-ray and optical light curves into individual flares, finding that X-ray leads optical variations (XO) in 6 flares, the reverse occurs in 3 flares (OX), and there is essentially zero lag in 4 flares. Upon comparing theoretical expectations with the data, we conclude that (1) XO flares can be explained by gradual acceleration of radiating electrons to the highest energies; (2) OX flares can result from either light-travel delays of the seed photons (synchrotron self-Compton scattering) or gradients in maximum electron energy behind shock fronts; and (3) events with similar X-ray and optical radiative energy output originate well upstream of the 43 GHz core, while those in which the optical radiative output dominates occur at or downstream of the core.

We present the third installment of HI sources extracted from the Arecibo Legacy Fast ALFA extragalactic survey. This dataset continues the work of the Virgo ALFALFA catalog. The catalogs and spectra published here consist of data obtained during the 2005 and 2006 observing sessions of the survey. The catalog consists of 578 HI detections within the range 11h 36m < R.A.(J2000) < 13h 52m and +08 deg < Dec.(J2000) < +12 deg, and cz_sun < 18000 km/s. The catalog entries are identified with optical counterparts where possible through the examination of digitized optical images. The catalog detections can be classified into three categories: (a) detections of high reliability with S/N > 6.5; (b) high velocity clouds in the Milky Way or its periphery; and (c) signals of lower S/N which coincide spatially with an optical object and known redshift. 75% of the sources are newly published HI detections. Of particular note is a complex of HI clouds projected between M87 and M49 that do not coincide with any optical counterparts. Candidate objects without optical counterparts are few. The median redshift for this sample is 6500 km/s and the cz distribution exhibits the local large scale structure consisting of Virgo and the background void and the A1367-Coma supercluster regime at cz_sun ~7000 km/s. Position corrections for telescope pointing errors are applied to the dataset by comparing ALFALFA continuum centroid with those cataloged in the NRAO VLA Sky Survey. The uncorrected positional accuracy averages 27 arcsec ~(21 arcsec ~median) for all sources with S/N > 6.5 and is of order ~21 arcsec ~(16 arcsec ~median) for signals with S/N > 12. Uncertainties in distances toward the Virgo cluster can affect the calculated HI mass distribution.

We report the discovery of a ~500 kpc HI extension southwest of the Virgo Cluster HI-rich pair NGC 4532/DDO 137, detected as part of the Arecibo Legacy Fast ALFA (ALFALFA) Survey. The feature is the longest and most massive HI tail structure so far found in the Virgo Cluster and, at 1.8 Mpc from M87, the most distant from the main concentration of the intracluster medium. The structure is spatially and spectrally separated into two ridges and is defined by diffuse emission and discrete clumps of mass 2.5 - 6.8 x 10**7 solar masses. All emission is blue-shifted with respect to the NGC 4532/DDO 137 pair emission. Including diffuse emission, the structure has a total mass of up to 7 x 10**8 solar masses, equivalent to ~10% of the system's HI mass. Optical R-band imaging finds no counterparts to a level of 26.5 mag arcsec**-2. The characteristics of the structure appear most consistent with a tidal origin.

We report the discovery by the Arecibo Legacy Fast ALFA (ALFALFA) survey of eight HI features not coincident with stellar counterparts in the Virgo Cluster region. All of the HI clouds have cz < 3000 km/s and, if at the Virgo distance, HI masses between 1.9 x 10**7 and 1.1 x 10**9 solar masses. Four of the eight objects were reported or hinted at by previous studies and "rediscovered" by ALFALFA. While some clouds appear to be associated with optical galaxies in their vicinity, others show no clear association with a stellar counterpart. Two of them are embedded in relatively dense regions of the cluster and are associated with M49 and M86; they were previously known. The others are mostly located in peripheral regions of the cluster. Especially notable are a concentration of objects towards the so-called M cloud, 3 to 5 degrees to the NW of M87, and a complex of several clouds projected roughly halfway between M87 and M49. The object referred to as VIRGOHI21 and proposed to be a "dark galaxy" is also detected and shown to be a tidal feature associated with NGC 4254.

During 2001-2002 the optically violent variable (OVV) blazar 3C 279 un- derwent the most intense outburst seen during the entire fourteen year history that this quasar has been studied at Colgate University's Foggy Bottom Obser- vatory (FBO). This study concentrates on ~1600 R-filter images taken during this period of activity. This data set includes twenty-nine nights of microvari- ability coverage. The outburst began in March 2001, after 3C 279 had faded to its faintest level, R = 15.5, in four years. The source reached its brightest level, R = 12.5, in the fourteen years of our study in August 2001, at which time it became unobservable due to its proximity to the Sun. Upon becoming observable again in mid-December 2001, 3C 279 fluctuated between R = 13.9 and R = 14.7, until a dramatic decrease in flux level in June-July 2002 brought the source back down to a level comparable to its pre-outburst state. The source exhibited numerous week-long flares of approximately one magnitude during the outburst period. Superposed on these flares were night-to-night variations of up to one half magnitude and intra-night microvariability of up to 0.13 magnitude in three hours. We use visual inspection of the light curve as well as numerical timescale analysis tools (the autocorrelation function, the structure function, and the power spectrum) to characterize the multiple timescales of variability ranging from 1.5 years to several hours.

We present the first installment of HI sources extracted from the Arecibo Legacy Fast ALFA (ALFALFA) extragalactic survey, initiated in 2005. Sources have been extracted from 3-D spectral data cubes and then examined interactively to yield global HI parameters. A total of 730 HI detections are catalogued within the solid angle 11h44m < R.A.(J2000) < 14h00m and +12deg < Dec.(J2000) < +16deg, and redshift range -1600 \kms < cz < 18000 \kms. In comparison, the HI Parkes All-Sky Survey (HIPASS) detected 40 HI signals in the same region. Optical counterparts are assigned via examination of digital optical imaging databases. ALFALFA HI detections are reported for three distinct classes of signals: (a) detections, typically with S/N > 6.5; (b) high velocity clouds in the Milky Way or its periphery; and (c) signals of lower S/N (to ~ 4.5) which coincide spatially with an optical object of known similar redshift. Although this region of the sky has been heavily surveyed by previous targeted observations based on optical flux-- or size-- limited samples, 69% of the extracted sources are newly reported HI detections. The resultant positional accuracy of HI sources is 20" (median). The median redshift of the sample is ~7000 \kms and its distribution reflects the known local large scale structure including the Virgo cluster. Several extended HI features are found in the vicinity of the Virgo cluster. A small percentage (6%) of HI detections have no identifiable optical counterpart, more than half of which are high velocity clouds in the Milky Way vicinity; the remaining 17 objects do not appear connected to or associated with any known galaxy.

The radio quasar 3C 454.3 underwent an exceptional optical outburst lasting more than 1 year and culminating in spring 2005. The maximum brightness detected was R = 12.0, which represents the most luminous quasar state thus far observed (M_B ~ -31.4). In order to follow the emission behaviour of the source in detail, a large multiwavelength campaign was organized by the Whole Earth Blazar Telescope (WEBT). Continuous optical, near-IR and radio monitoring was performed in several bands. ToO pointings by the Chandra and INTEGRAL satellites provided additional information at high energies in May 2005. The historical radio and optical light curves show different behaviours. Until about 2001.0 only moderate variability was present in the optical regime, while prominent and long-lasting radio outbursts were visible at the various radio frequencies, with higher-frequency variations preceding the lower-frequency ones. After that date, the optical activity increased and the radio flux is less variable. This suggests that the optical and radio emissions come from two separate and misaligned jet regions, with the inner optical one acquiring a smaller viewing angle during the 2004-2005 outburst. Moreover, the colour-index behaviour (generally redder-when-brighter) during the outburst suggests the presence of a luminous accretion disc. A huge mm outburst followed the optical one, peaking in June-July 2005. The high-frequency (37-43 GHz) radio flux started to increase in early 2005 and reached a maximum at the end of our observing period (end of September 2005). VLBA observations at 43 GHz during the summer confirm the

Results of an international campaign to photometrically monitor the unique pre-main sequence eclipsing object KH 15D are reported. An updated ephemeris for the eclipse is derived that incorporates a slightly revised period of 48.36 d. There is some evidence that the orbital period is actually twice that value, with two eclipses occurring per cycle. The extraordinary depth (~3.5 mag) and duration (~18 days) of the eclipse indicate that it is caused by circumstellar matter, presumably the inner portion of a disk. The eclipse has continued to lengthen with time and the central brightness reversals are not as extreme as they once were. V-R and V-I colors indicate that the system is slightly bluer near minimum light. Ingress and egress are remarkably well modeled by the passage of a knife-edge across a limb-darkened star. Possible models for the system are briefly discussed.

We report on two BeppoSAX observations of BL Lac (2200+420) performed respectively in June and December 1999, as part of a ToO program to monitor blazars in high states of activity. During both runs the source has been detected up to 100 keV, but it showed quite different spectra: in June it was concave with a very hard component above 5-6 keV (alpha_1 ~ 1.6; alpha_2 ~ 0.15); in December it was well fitted by a single power law (alpha ~ 0.6). During the first BeppoSAX observation BL Lac showed an astonishing variability episode: the 0.3 - 2 keV flux doubled in ~ 20 minutes, while the flux above 4 keV was almost constant. This frequency-dependent event is one of the shortest ever recordered for BL Lac objects and places lower limits on the dimension and the magnetic field of the emitting region and on the energy of the synchrotron radiating electrons. A similar but less extreme behaviour is detected also in optical light curves, that display non-simultaneous, smaller fluctuations of ~ 20% in 20 min. We fit the spectral energy distributions with a homogeneous, one-zone model to constrain the emission region in a very simple but effective SSC + external Compton scenario, highlightning the importance of the location of the emitting region with respect to the Broad Line Region and the relative spectral shape dependence. We compare our data with historical radio to gamma-ray Spectral Energy Distributions.

New optical and radio data on the BL Lacertae object AO 0235+16 have been collected in the last four years by a wide international collaboration, which confirm the intense activity of this source. The optical data also include the results of the Whole Earth Blazar Telescope (WEBT) first-light campaign organized in November 1997. The optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence of typical variability time scales and of possible correlations between the optical and radio emissions by means of visual inspection, Discrete Correlation Function analysis, and Discrete Fourier Transform technique. The major radio outbursts are found to repeat quasi-regularly with a periodicity of about 5.7 years; this period is also in agreement with the occurrence of some of the major optical outbursts, but not all of them.

Light curves of 3C 279 are presented in optical (R-band), X-rays (RXTE/PCA), and gamma rays (CGRO/EGRET) for 1999 Jan-Feb and 2000 Jan-Mar. During both of those epochs the gamma-ray levels were high, and all three observed bands demonstrated substantial variation, on time scales as short as one day. Correlation analyses provided no consistent pattern, although a rather significant optical/gamma-ray correlation was seen in 1999, with a gamma-ray lag of ~2.5 days, and there are other suggestions of correlations in the light curves. For comparison, correlation analysis is also presented for the gamma-ray and X-ray light curves during the large gamma ray flare in 1996 Feb and the two gamma-bright weeks leading up to it; the correlation at that time was strong, with a gamma-ray/X-ray offset of no more than 1 day.

We present the results of low dispersion optical spectroscopy of 186 H II regions spanning a range of radius in 13 spiral galaxies. Abundances for several elements (oxygen, nitrogen, neon, sulfur, and argon) were determined for 185 of the H II regions. As expected, low metallicities were found for the outlying H II regions of these spiral galaxies. Radial abundance gradients were derived for the 11 primary galaxies; similar to results for other spiral galaxies, the derived abundance gradients are typically -0.04 to -0.07 dex/kpc.