We present UV broadband photometry and optical emission-line measurements for a sample of 32 Brightest Cluster Galaxies (BCGs) in clusters of the Representative XMM-Newton Cluster Structure Survey (REXCESS) with z = 0.06-0.18. The REXCESS clusters, chosen to study scaling relations in clusters of galaxies, have X-ray measurements of high quality. The trends of star formation and BCG colors with BCG and host properties can be investigated with this sample. The UV photometry comes from the XMM Optical Monitor, supplemented by existing archival GALEX photometry. We detected H\alpha and forbidden line emission in 7 (22%) of these BCGs, in optical spectra. All of the emission-line BCGs occupy clusters classified as cool cores, for an emission-line incidence rate of 70% for BCGs in cool core clusters. Significant correlations between the H\alpha equivalent widths, excess UV production in the BCG, and the presence of dense, X-ray bright intracluster gas with a short cooling time are seen, including the fact that all of the H\alpha emitters inhabit systems with short central cooling times and high central ICM densities. Estimates of the star formation rates based on H\alpha and UV excesses are consistent with each other in these 7 systems, ranging from 0.1-8 solar masses per year. The incidence of emission-line BCGs in the REXCESS sample is intermediate, somewhat lower than in other X-ray selected samples (-35%), and somewhat higher than but statistically consistent with optically selected, slightly lower redshift BCG samples (-10-15%). The UV-optical colors (UVW1-R-4.7\pm0.3) of REXCESS BCGs without strong optical emission lines are consistent with those predicted from templates and observations of ellipticals dominated by old stellar populations. We see no trend in UV-optical colors with optical luminosity, R-K color, X-ray temperature, redshift, or offset between X-ray centroid and X-ray peak (<w>).

We investigate the relationship between brightest cluster galaxies (BCGs) and their host clusters using a sample of nearby galaxy clusters from the Representative XMM Cluster Structure Survey (REXCESS). The sample was imaged with the Southern Observatory for Astrophysical Research (SOAR) in R band to investigate the mass of the old stellar population. Using a metric radius of 12h^-1 kpc, we found that the BCG luminosity depends weakly on overall cluster mass as L_BCG \propto M_cl^0.18+-0.07, consistent with previous work. We found that 90% of the BCGs are located within 0.035 r_500 of the peak of the X-ray emission, including all of the cool core (CC) clusters. We also found an unexpected correlation between the BCG metric luminosity and the core gas density for non-cool core (non-CC) clusters, following a power law of n_e \propto L_BCG^2.7+-0.4 (where n_e is measured at 0.008 r_500). The correlation is not easily explained by star formation (which is weak in non-CC clusters) or overall cluster mass (which is not correlated with core gas density). The trend persists even when the BCG is not located near the peak of the X-ray emission, so proximity is not necessary. We suggest that, for non-CC clusters, this correlation implies that the same process that sets the central entropy of the cluster gas also determines the central stellar density of the BCG, and that this underlying physical process is likely to be mergers.

In 1981, a faint radio source (G') was detected near the center of the lensing galaxy of the famous "twin quasar" Q0957+561. It is still unknown whether this central radio source is a third quasar image or an active nucleus of the lensing galaxy, or a combination of both. In an attempt to resolve this ambiguity, we observed Q0957+561 at radio wavelengths of 13cm, 18cm, and 21cm, using the Very Long Baseline Array in combination with the phased Very Large Array and the Green Bank Telescope. We measured the spectrum of G' for the first time and found it to be significantly different from the spectra of the two bright quasar images. This finding suggests that the central component is primarily or entirely emission from the foreground lens galaxy, but the spectrum is also consistent with the hypothesis of a central quasar image suffering free-free absorption. In addition, we confirm the previously-reported VLBI position of G' just north of the optical center of the lens galaxy. The position slightly favors the hypothesis that G' originates in the lens, but is not conclusive. We discuss the prospects for further clarification of this issue.

We present results from a survey for gravitationally lensed radio lobes. Lensed lobes are a potentially richer source of information about galaxy mass distributions than lensed point sources, which have been the exclusive focus of other recent surveys. Our approach is to identify radio lobes in the FIRST catalog and then search optical catalogs for coincident foreground galaxies, which are candidate lensing galaxies. We then obtain higher-resolution images of these targets at both optical and radio wavelengths, and obtain optical spectra for the most promising candidates. We present maps of several radio lobes that are nearly coincident with galaxies. We have not found any new and unambiguous cases of gravitational lensing. One radio lobe in particular, FOV J0743+1553, has two hot spots that could be multiple images produced by a z=0.19 spiral galaxy, but the lensing interpretation is problematic.

Different wavelength regimes and methods for estimating the space density of the star-formation rate (SFR) result in discrepant values. While it is recognised that ultra-violet (UV) and H-alpha emission line data must be corrected for the effects of extinction, the magnitude of the required correction is uncertain. Even when these corrections are made there remains a significant discrepancy between SFRs derived from UV and H-alpha measurements compared with those derived from far-infrared (FIR) and radio luminosities. Since the FIR/radio derived SFRs are not affected by extinction, and simple corrections to reconcile the UV and H-alpha measurement with these do not fully account for the discrepancies, a more sophisticated correction may be required. Recent results suggest that at least part of the solution may be a form of extinction which increases with increasing SFR (or luminosity, given the common assumption that SFR is proportional to luminosity). We present an analysis of the effects of a dust reddening dependent on star formation rate applied to estimators of SFR. We show (1) that the discrepancies between H-alpha and FIR/radio SFR estimates may be explained by such an effect, and we present an iterative method for applying the correction; and (2) UV-based estimates of SFR are harder to reconcile with FIR/radio estimates using this method, although the extent of the remaining discrepancy is less than for a non-SFR-dependent correction. Particularly at high redshift, our understanding of extinction at UV wavelengths may require a still more complex explanation.

High precision VLA flux density measurements for the lensed images of 0218+357 yield a time delay of 10.1(+1.5-1.6)days (95% confidence). This is consistent with independent measurements carried out at the same epoch (Biggs et al. 1999), lending confidence in the robustness of the time delay measurement. However, since both measurements make use of the same features in the light curves, it is possible that the effects of unmodelled processes, such as scintillation or microlensing, are biasing both time delay measurements in the same way. Our time delay estimates result in confidence intervals that are somewhat larger than those of Biggs et al., probably because we adopt a more general model of the source variability, allowing for constant and variable components. When considered in relation to the lens mass model of Biggs et al., our best-fit time delay implies a Hubble constant of H_o = 71(+17-23) km/s-Mpc for Omega_o=1 and lambda_o=0 (95% confidence; filled beam). This confidence interval for H_o does not reflect systematic error, which may be substantial, due to uncertainty in the position of the lens galaxy. We also measure the flux ratio of the variable components of 0218+357, a measurement of a small region that should more closely represent the true lens magnification ratio. We find ratios of 3.2(+0.3-0.4) (95% confidence; 8 GHz) and 4.3(+0.5-0.8) (15 GHz). Unlike the reported flux ratios on scales of 0.1", these ratios are not strongly significantly different. We investigate the significance of apparent differences in the variability properties of the two images of the background active galactic nucleus. We conclude that the differences are not significant, and that time series much longer than our 100-day time series will be required to investigate propagation effects in this way.

Time-delay measurements of 0957+561 based on radio monitoring now agree with optical monitoring. Recent models incorporate many recent observations, but the systematic uncertainties in the models still dominate the uncertainty in the cosmological results. VLBI observations of the milli-arcsec jet structure have provided the most important set of modeling constraints. We present new observations of this structure, made at 18 cm with the Very Long Baseline Array, the Very Large Array, and the Green Bank 140 ft. telescope. Compared to the data set of Garrett et al. (1994), ours has a similar theoretical noise level (20 microJy), but nearly twice as many baselines (66 compared to 36) and more uniform coverage of the UV plane. Our initial maps confirm the basic jet structures seen by Garrett et al., with comparable RMS noise (60-80 microJy), and are in reasonable agreement with the jet structure model of Barkana et al. (1999). The maps also hint at a new compact component leaving the core, and a previously undetected diffuse component at the end of the jet. Further refinements in fringe-fitting, mapping, and self-calibration should yield significantly lower RMS noise (closer to the theoretical level). We will use a modified form of the CalTech VLBI Modelfit program (written by Barkana et al. 1999) to determine the magnification matrix between the images. The refined magnification matrix should provide stronger constraints on the lens model, and thus reduce the uncertainty in the cosmological results.

The centimeter-wave luminosity of local radio galaxies correlates well with their star formation rate. We extend this correlation to surveys of high-redshift radio sources to estimate the global star formation history. The star formation rate found from radio observations needs no correction for dust obscuration, unlike the values calculated from optical and ultraviolet data. Three deep radio surveys have provided catalogs of sources with nearly complete optical identifications and nearly 60% complete spectroscopic redshifts: the Hubble Deep Field and Flanking Fields at 12h+62d, the SSA13 field at 13h+42d, and the V15 field at 14h+52d. We use the redshift distribution of these radio sources to constrain the evolution of their luminosity function. The epoch dependent luminosity function is then used to estimate the evolving global star formation density. At redshifts less than one, our calculated star formation rates are significantly larger than even the dust-corrected optically-selected star formation rates; however, we confirm the rapid rise from z=0 to z=1 seen in those surveys.

Faint extragalactic radio sources provide important information about the global history of star formation. Sensitive radio observations of the Hubble Deep Field and other fields have found that sub-mJy radio sources are predominantly associated with star formation activity rather than AGN. Radio observations of star forming galaxies have the advantage of being independent of extinction by dust. We use the FIR-radio correlation to compare the radio and FIR backgrounds, and make several conclusions about the star forming galaxies producing the FIR background. We then use the redshift distribution of faint radio sources to determine the evolution of the radio luminosity function, and thus estimate the star formation density as a function of redshift.

The gravitational lens 0957+561 was monitored with the Very Large Array from 1979 to 1997. The 6 cm light curve data from 1995-1997 and the 4 cm data from 1990-1997 are reported here. At 4 cm, the intrinsic source variations occur earlier and are twice as large as the corresponding variations at 6 cm. The VLBI core and jet components have different magnification factors, leading to different flux ratios for the varying and non-varying portions of the VLA light curves. Using both the PRHQ and Dispersion statistical techniques, we determined the time delay, core flux ratio, and excess non-varying B image flux density. The fits were performed for the 4 cm and 6 cm light curves, both individually and jointly, and we used Gaussian Monte Carlo data to estimate 68% statistical confidence levels. The delay estimates from each individual wavelength were inconsistent given the formal uncertainties, suggesting that there are unmodeled systematic errors in the analysis. We roughly estimate the systematic uncertainty in the joint result from the difference between the 6 cm and 4 cm results, giving 409+-30 days for the PRHQ statistic and 397+-20 days for the Dispersion statistic. These results are consistent with the current optical time delay of 417+-3 days, reconciling the long-standing difference between the optical and radio light curves and between different statistical analyses. The unmodeled systematic effects may also corrupt light curves for other lenses, and we caution that multiple events at multiple wavelengths may be necessary to determine an accurate delay in any lens system. Now that consensus has been reached regarding the time delay in the 0957+561 system, the most pressing issue remaining for determining H_0 is a full understanding of the mass distribution in the lens.

The cosmic far infrared background detected recently by the COBE-DIRBE team is presumably due, in large part, to the far infrared (FIR) emission from all galaxies. We take the well-established correlation between FIR and radio luminosity for individual galaxies and apply it to the FIR background. We find that these sources make up about half of the extragalactic radio background, the other half being due to AGN. This is in agreement with other radio observations, which leads us to conclude that the FIR-radio correlation holds well for the very faint sources making up the FIR background, and that the FIR background is indeed due to star-formation activity (not AGN or other possible sources). If these star-forming galaxies have a radio spectral index between 0.4 and 0.8, and make up 40 to 60% of the extragalactic radio background, we find that they have redshifts between roughly 1 and 2, in agreement with recent estimates by Madau et al. of the redshift of peak star-formation activity. We compare the observed extragalactic radio background to the integral over the logN-logS curve for star-forming radio sources, and find that the slope of the curve must change significantly below about 1 microjansky. At 1 microjansky, the faint radio source counts predict about 25 sources per square arcminute, and these will cause SIRTF to be confusion limited at 160micron.

We have produced deep radio maps of the double quasar 0957+561 from multiple-epoch VLA observations. To achieve high sensitivity to extended structure we have re-reduced the best available 1.6 GHz observations and have combined 5 GHz data from multiple array configurations. Regions of faint emission approximately 15 arcsec north and south of the radio source G are probably lobes associated with the lensing galaxy. An arc 5 arcsec to the east of G may be a stretched image of emission in the background quasar's environment. 1.4 arcsec southwest of G we detect a source that we interpret as an image of emission from the quasar's western lobe, which could provide a constraint on the slope of the gravitational potential in the central region of the lens. We explore the consequences of these new constraints with simple lens models of the system.

We report on 15 years of VLA monitoring of the gravitational lens B0957+561 at 6 cm. Since our last report in 1992, there have been 32 additional observations, in which both images have returned to their quiescent flux density levels and the A image has brightened again. We estimate the time delay from the light curves using three different techniques: the chi-squared analysis of Press, Rybicki, & Hewitt (1992a,b), the dispersion analysis of Pelt et al. (1994, 1996), and the locally normalized discrete correlation function of Leh\'ar et al. (1992). Confidence intervals for these time delay estimates are found using Monte Carlo techniques. With the addition of the new observations, it has become obvious that five observations from Spring 1990 are not consistent with the statistical properties of the rest of the light curves, so we analyze the light curves with those points removed, as well as the complete light curves. The three statistical techniques applied to the two versions of the data set result in time delay values in the range 398 to 461 days (or 1.09 to 1.26 years, A leading B), each with ~5% formal uncertainty. The corresponding flux ratios (B/A) are in the range 0.698 to 0.704. Thus, the new features in the light curve show that the time delay is less than 500 days, in contrast with analysis of earlier versions of the radio light curves. The large range in the time delay estimates is primarily due to unfortunate coincidences of observing gaps with flux variations.