Search

We present Herschel observations of the core of the Perseus cluster of galaxies. The brightest cluster galaxy, NGC 1275, is surrounded by a network of filaments previously imaged extensively in H{\alpha} and CO. In this work, we report detections of FIR lines with Herschel. All but one of the lines are spatially extended, with the [CII] line emission extending up to 25 kpc from the core. There is spatial and kinematical correlation among [CII], H{\alpha} and CO, which gives us confidence to model the different components of the gas with a common heating model. With the help of FIR continuum Herschel measurements, together with a suite of coeval radio, submm and infrared data, we performed a SED fitting of NGC 1275 using a model that contains contributions from dust emission as well as synchrotron AGN emission. The data indicate a low dust emissivity index, beta ~ 1, a total dust mass close to 10^7 solar mass, a cold dust component with temperature 38 \pm 2 K and a warm dust component with temperature of 116 \pm 9 K. The FIR-derived star formation rate (SFR) is 24 \pm 1 solar mass per yr, in close agreement with the FUV-derived SFR. We investigated in detail the source of the Herschel FIR and H{\alpha} emissions emerging from a core region 4 kpc in radius. Based on simulations conducted using the radiative transfer code, CLOUDY, a heating model comprising old and young stellar populations is sufficient to explain these observations. We have also detected [CII] in three well-studied regions of the filaments. We find a [OI]/[CII] ratio about 1 dex smaller than predicted by the otherwise functional Ferland (2009) model. The line ratio suggests that the lines are optically thick, as is typical of galactic PDRs, and implies that there is a large reservoir of cold atomic gas. [abridged], Comment: 21 pages, 16 figures and 10 tables. Accepted for publication in MNRAS

Beam power is a fundamental parameter that describes, in part, the state of a supermassive black hole system. Determining the beam powers of powerful classical double radio sources requires substantial observing time, so it would be useful to determine the relationship between beam power and radio power so that radio power could be used as a proxy for beam power. A sample of 31 powerful classical double radio sources with previously determined beam and radio powers are studied; the sources have redshifts between about 0.056 and 1.8. It is found that the relationship between beam power, Lj, and radio power, P, is well described by Log(Lj) = 0.84 Log(P) + 2.15, where both L_j and P are in units of 10^(44) erg/s. This indicates that beam power is converted to radio power with an efficiency of about 0.7%. The ratio of beam power to radio power is studied as a function of redshift; there is no significant evidence for redshift evolution of this ratio over the redshift range studied. The relationship is consistent with empirical results obtained by Cavagnolo et al. (2010) for radio sources in gas rich environments, which are primarily FRI sources, and with the theoretical predictions of Willott et al. (1999)., Comment: 6 pages, 2 figures, 2 tables; accepted for publication in MNRAS

Brightest cluster galaxies (BCGs) in the cores of galaxy clusters have distinctly different properties from other low redshift massive ellipticals. The majority of the BCGs in cool-core clusters show signs of active star formation. We present observations of NGC 4696, the BCG of the Centaurus galaxy cluster, at far-infrared (FIR) wavelengths with the Herschel space telescope. Using the PACS spectrometer, we detect the two strongest coolants of the interstellar medium, CII at 157.74 micron and OI at 63.18 micron, and in addition NII at 121.90 micron. The CII emission is extended over a region of 7 kpc with a similar spatial morphology and kinematics to the optical H-alpha emission. This has the profound implication that the optical hydrogen recombination line, H-alpha, the optical forbidden lines, NII 6583 Angstrom, the soft X-ray filaments and the far-infrared CII line all have the same energy source. We also detect dust emission using the PACS and SPIRE photometers at all six wavebands. We perform a detailed spectral energy distribution fitting using a two-component modified black-body function and find a cold 19 K dust component with mass 1.6x10^6 solar mass and a warm 46 K dust component with mass 4.0x10^3 solar mass. The total FIR luminosity between 8 micron and 1000 micron is 7.5x10^8 solar luminosity, which using Kennicutt relation yields a low star formation rate of 0.13 solar mass per yr. This value is consistent with values derived from other tracers, such as ultraviolet emission. Combining the spectroscopic and photometric results together with optical H-alpha, we model emitting clouds consisting of photodissociation regions (PDRs) adjacent to ionized regions. We show that in addition to old and young stellar populations, there is another source of energy, such as cosmic rays, shocks or reconnection diffusion, required to excite the H-alpha and CII filaments., Comment: 17 pages, 9 figures, 5 tables, accepted for publication in MNRAS

Based on Chandra and ROSAT observations, we investigated the nearby poor cluster Abell 194, which hosts two luminous radio galaxies, NGC547 (3C 40B) and NGC541 (3C 40A). We demonstrated the presence of a large X-ray cavity (r~34 kpc) formed by the giant southern radio lobe arising from 3C 40B in NGC547. The estimated age of the cavity is t=7.9 x 10^7 years and the total work of the AGN is 3.3 x 10^59 erg, hence the cavity power is P_cav=1.3 x 10^44 erg/s. Furthermore, in the Chandra images of NGC545 and NGC541 we detected sharp surface brightness edges, identified as merger cold fronts, and extended tails. Using the pressure ratios between inside and outside the cold fronts we estimated that the velocities of NGC545 and NGC541 correspond to Mach-numbers of M=1.0^{+0.3}_{-0.5} and M=0.9^{+0.2}_{-0.5}, respectively. The low radial velocities of these galaxies relative to the mean radial velocity of Abell 194 imply that their motion is oriented approximately in the plane of the sky. Based on these and earlier observations, we concluded that NGC545 and NGC541 are falling through the cluster, whose center is NGC547, suggesting that Abell 194 is undergoing a significant cluster merger event. Additionally, we detected 20 bright X-ray sources around NGC547 and NGC541, a surprisingly large number, since the predicted number of resolved LMXBs and CXB sources is 2.2 and 4.1, respectively. To explain the nature of additional sources, different possibilities were considered, none of which are satisfactory. We also studied the origin of X-ray emission in Minkowski's Object, and concluded that it is most likely dominated by the population of HMXBs rather than by hot diffuse ISM. Moreover, in view of the galaxy dynamics in Abell 194, we explored the possibility that the starburst in Minkowski's Object was triggered by its past interaction with NGC541, and concluded that it may be a viable path., Comment: 12 pages, 7 figures, accepted for publication in ApJ

Quillen et al. and O'Dea et al. carried out a Spitzer study of a sample of 62 brightest cluster galaxies (BCGs) from the ROSAT Brightest Cluster Sample chosen based on their elevated H-alpha flux. We present Hubble Space Telescope Advanced Camera for Surveys (ACS) far ultraviolet (FUV) images of the Ly-alpha and continuum emission of the luminous emission-line nebulae in 7 BCGs found to have an infrared excess. We confirm that the BCGs are actively forming stars suggesting that the IR excess seen in these BCGs is indeed associated with star formation. The FUV continuum emission extends over a region of ~7-28 kpc (largest linear size) and even larger in Ly-alpha. The young stellar population required by the FUV observations would produce a significant fraction of the ionizing photons required to power the emission line nebulae. Star formation rates estimated from the FUV continuum range from ~3 to ~14 times lower than those estimated from the IR, however both the Balmer decrement in the central few arcseconds and detection of CO in most of these galaxies imply that there are regions of high extinction that could have absorbed much of the FUV continuum. Analysis of archival VLA observations reveals compact radio sources in all 7 BCGs and kpc scale jets in Abell 1835 and RXJ 2129+00. The four galaxies with archival deep Chandra observations exhibit asymmetric X-ray emission, the peaks of which are offset from the center of the BCGs by ~10 kpc on average. A low feedback state for the AGN could allow increased condensation of the hot gas into the center of the galaxy and the feeding of star formation., Comment: Accepted for publication in The Astrophysical Journal

We present Hubble Space Telescope UV and optical imaging of the radio galaxy 3C 236, whose relic 4 Mpc radio jet lobes and inner 2 kpc CSS radio source are evidence of multiple epochs of AGN activity. Our data confirm the presence of four bright knots of FUV emission in an arc along the edge of the inner circumnuclear dust disk in the galaxy's nucleus, as well as FUV emission cospatial with the nucleus itself. We interpret these to be sites of recent or ongoing star formation. We present photometry of these knots, as well as an estimate for the internal extinction in the source using the Balmer decrement from SDSS spectroscopy. We estimate the ages of the knots by comparing our extinction-corrected photometry with stellar population synthesis models. We find the four knots cospatial with the dusty disk to be young, of order 10^7 yr old. The FUV emission in the nucleus is likely due to an episode of star formation triggered ~10^9 yr ago. We argue that the young 10^7 yr old knots stem from an episode of star formation that was roughly coeval with the event resulting in reignition of radio activity, creating the CSS source. The 10^9 yr old stars in the nucleus may be associated with the previous epoch of activity that generated the 4 Mpc relic source, before it was cut off by exhaustion or interruption. The ages of the knots, considered in context with the disturbed morphology of the nuclear dust and the double-double morphology of the "old" and "young" radio sources, present evidence for an episodic AGN/starburst connection. We suggest that the AGN fuel supply was interrupted for ~10^7 yr due to a minor merger event and has now been restored, and the resultant non-steady flow of gas toward the nucleus is likely responsible for both the new episode of infall-induced star formation and also the multiple epochs of radio activity., Comment: 15 pages, 7 figures. Accepted for publication in ApJ.

We compare the near-infrared (NIR) H band photometric and morphological properties of low-redshift (z<0.3) 3CR radio galaxies with samples of BL Lac object and quasar host galaxies, merger remnants, quiescent elliptical galaxies, and brightest cluster galaxies drawn from the literature. In general the 3CR host galaxies are consistent with luminous (~L*) elliptical galaxies. The vast majority of FR II's (~80%) occupy the most massive ellipticals and form a homogeneous population that is comparable to the population of radio-loud quasar (RLQ) host galaxies in the literature. However, a significant minority (~20%) of the 3CR FR II's appears under-luminous with respect to quasar host galaxies. All FR II objects in this faint tail are either unusually red, or appear to be the brightest objects within a group. We discuss the apparent differences between the radio galaxy and RLQ host galaxy populations. RLQs appear to require >1E11 M_sun host galaxies (and ~1E9 M_sun black holes), whereas radio galaxies and RQQs can exist in galaxies down to 3E10 M_sun. This may be due to biases in the measured quasar host galaxy luminosities or populations studied, or due to a genuine difference in host galaxy. If due to a genuine difference, it would support the idea that radio and optical active galactic nucleii are two separate populations with a significant overlap., Comment: 17 pages, 15 figures. Accepted for publication in ApJ.

We present an analysis of Spitzer IRS spectroscopy of 83 active galaxies from the extended 12 micron sample. We find rank correlations between several tracers of star formation which suggest that (1) the PAH feature is a reliable tracer of star formation, (2) there is a significant contribution to the heating of the cool dust by stars, (3) the H$_2$ emission is also primarily excited by star formation. The 55-90 vs. 20-30 spectral index plot is also a diagnostic of the relative contribution of Starburst to AGN. We see there is a large change in spectral index across the sample. Thus, the contribution to the IR spectrum from the AGN and starburst components can be comparable in magnitude but the relative contribution also varies widely across the sample. We find rank correlations between several AGN tracers. We search for correlations between AGN and Starburst tracers and we conclude that the AGN and Starburst tracers are not correlated. This is consistent with our conclusion that the relative strength of the AGN and Starburst components varies widely across the sample. Thus, there is no simple link between AGN fueling and Black Hole Growth and star formation in these galaxies. The distribution of Sil 10 micron and 18 micron strengths is consistent with the clumpy torus models of Sirocky et al. We find a rank correlation between the [NeV] 14 micron line and the 6.7 micron continuum which may be due to an extended component of hot dust. The Sy 2s with a Hidden Broad Line Region (HBLR) have a higher ratio of AGN to Starburst contribution to the SED than Sy 2s without an HBLR. This may contribute to the detection of the HBLR in polarized light. The Sy 2s with an HBLR are more similar to the Sy 1s than they are to the Sy 2s without an HBLR.

Radio lobes of a sample of eleven very powerful classical double radio galaxies were studied. Each source was rotated so that the symmetry axis of the source was horizontal, and vertical cross-sectional cuts were taken across the source at intervals of one beam size. These were used to study the cross-sectional surface brightness profiles, the width of each slice, radio emissivity as a function of position across each slice, the first and second moments, and the average surface brightness, minimum energy magnetic field strength, and pressure of each slice. A Gaussian provides a good description of the surface brightness profile of cross-sectional slices. The Gaussian FWHM as a function of distance from the hot spot first increases and then decreases with distance from the hot spot. The width as a function of distance from the hot spot is highly symmetric on each side of the source. The radio emissivity is often close to flat across a slice, indicating a roughly constant emissivity and pressure for that slice. Some slices show variations in radio emissivity that indicate an ``edge-peaked'' pressure profile for that slice; these often occur in slices near the local maxima of the bridge width. The emissivity does not exhibit any signature of emission from a jet. The first moment is generally quite close to zero indicating only small excursions of the ridge line from the symmetry axis of the source. The second moment indicates the same source shape as is found using the Gaussian FWHM. The average magnetic field strength and pressure decrease with increasing distance from the hot spot, reaching a roughly constant value at a location that is typically just before the location of a local maximum of the bridge width. These results are interpreted in terms of a heuristic model for the radio lobes., Comment: 102 pages, 136 figures, accepted for publication in ApJ Supplement Series

We present a comprehensive multiband spectral and polarimetric study of the jet of 3C 264 (NGC 3862). Included in this study are three HST optical and ultraviolet polarimetry data sets, along with new and archival VLA radio imaging and polarimetry, a re-analysis of numerous HST broadband data sets from the near infrared to the far ultraviolet, and a Chandra ACIS-S observation. We investigate similarities and differences between optical and radio polarimetry, in both degree of polarization and projected magnetic field direction. We also examine the broadband spectral energy distribution of both the nucleus and jet of 3C 264, from the radio through the X-rays. From this we place constraints on the physics of the 3C 264 system, the jet and its dynamics. We find significant curvature of the spectrum from the near-IR to ultraviolet, and synchrotron breaks steeper than 0.5, a situation also encountered in the jet of M87. This likely indicates velocity and/or magnetic field gradients and more efficient particle acceleration localized in the faster/higher magnetic field parts of the flow. The magnetic field structure of the 3C 264 jet is remarkably smooth; however, we do find complex magnetic field structure that is correlated with changes in the optical spectrum. We find that the X-ray emission is due to the synchrotron process; we model the jet spectrum and discuss mechanisms for accelerating particles to the needed energies, together with implications for the orientation of the jet under a possible spine-sheath model., Comment: 18 pages, 14 figures, 1 table; Ap J, in press (Jan 1 2010 issue)

We report on the first high frequency VLBI observations of the nearby broad absorption line quasar (BALQSO), Mrk 231. Three epochs of observations were achieved at 15 GHz and 22 GHz, two of these included 43 GHz observations as well. The nuclear radio source is resolved as a compact double. The core component experienced a strong flare in which the flux density at 22 GHz increased by $> 150%$ (45 mJy) in three months. Theoretical models of the flare imply that the emission is likely enhanced by very strong Doppler boosting of a highly relativistic ejecta with a kinetic energy flux, $Q \sim 3 \times 10^{43} \mathrm{ergs/sec}$. Combining our data with two previous epochs of 15 GHz data, shows marginal evidence for the slow advance of the secondary component (located $\approx 0.97$ pc from the core) over a 9.4 year span. We estimate that the long term time averaged kinetic energy flux of the secondary at $\bar{Q}\approx 10^{42}\mathrm{ergs/sec}$. Low frequency VLBA observations indicate that the secondary is seen through a shroud of free-free absorbing gas with an emission measure of $\approx 10^{8} \mathrm{cm}^{-6}\mathrm{pc}$. The steep spectrum secondary component appears to be a compact radio lobe that is associated with a working surface between the ram-pressure confined jet, and a dense medium that is likely to be the source of the free-free absorption. The properties of the dense gas are consistent with the temperatures, displacement from the nucleus and the column density of total hydrogen commonly associated with the BAL wind., Comment: To appear in ApJ

We present 19 nearby (z<0.3) 3CR radio galaxies imaged at low- and high-excitation as part of a Cycle 15 Hubble Space Telescope snapshot survey with the Advanced Camera for Surveys. These images consist of exposures of the H-alpha (6563 \AA, plus [NII] contamination) and [OIII] 5007 \AA emission lines using narrow-band linear ramp filters adjusted according to the redshift of the target. To facilitate continuum subtraction, a single-pointing 60 s line-free exposure was taken with a medium-band filter appropriate for the target's redshift. We discuss the steps taken to reduce these images independently of the automated recalibration pipeline so as to use more recent ACS flat-field data as well as to better reject cosmic rays. We describe the method used to produce continuum-free (pure line-emission) images, and present these images along with qualitative descriptions of the narrow-line region morphologies we observe. We present H-alpha+[NII] and [OIII] line fluxes from aperture photometry, finding the values to fall expectedly on the redshift-luminosity trend from a past HST/WFPC2 emission line study of a larger, generally higher redshift subset of the 3CR. We also find expected trends between emission line luminosity and total radio power, as well as a positive correlation between the size of the emission line region and redshift. We discuss the associated interpretation of these results, and conclude with a summary of future work enabled by this dataset., Comment: 18 pages, 12 figures, accepted for publication in ApJS

We estimate ages and physical properties of powerful radio galaxies. An analysis of new multi-wavelength VLA observations of eleven very powerful classical double (FRIIb) radio galaxies with redshifts between 0.4 and 1.3 is presented. We estimate ages and velocities for each side of each source. The eleven new sources are combined with previously studied samples and the characteristics of the full sample of 31 sources are studied; the full sample includes sources with redshifts between 0.056 and 1.79, and core-hot spot sizes of about 30 to 400 kpc. The velocities are independent of core-hotspot separation, suggesting the rate of growth of a given source is roughly constant over the source lifetime. We combine the rate of growth, width, and pressure of a source to study the beam power, lifetime, energy, and ambient gas density using standard methods previously applied to smaller samples. Typical beam powers are in the range from 1E44 to 1E46 erg/s; we show that this quantity is insensitive to assumptions regarding minimum energy conditions. The beam powers are independent of core-hotspot separation suggesting that the beam power of a given source is roughly constant over the source lifetime. Typical total source lifetimes are found to be about a few E6 to E7 years, and typical total outflow energies (E/c^2) are found to be about 5(E5 - E6) solar masses. Ambient gas densities are found to decrease with increasing core-hotspot distance, but have no redshift dependence. Overall, the results obtained with the sample of 31 sources studied here are consistent with those obtained earlier with smaller samples., Comment: Astronomy and Astrophysics, accepted

Quillen et al.(2007) presented an imaging survey with the {\it Spitzer Space Telescope} of 62 brightest cluster galaxies with optical line emission located in the cores of X-ray luminous clusters. They found that at least half of these sources have signs of excess infrared emission. Here we discuss the nature of the IR emission and its implications for cool core clusters. The strength of the mid-IR excess emission correlates with the luminosity of the optical emission lines. Excluding the four systems dominated by an AGN, the excess mid-infrared emission in the remaining brightest cluster galaxies is likely related to star formation. The mass of molecular gas (estimated from CO observations) is correlated with the IR luminosity as found for normal star forming galaxies. The gas depletion time scale is about 1 Gyr. The physical extent of the infrared excess is consistent with that of the optical emission line nebulae. This supports the hypothesis that the star formation occurs in molecular gas associated with the emission line nebulae and with evidence that the emission line nebulae are mainly powered by ongoing star formation. We find a correlation between mass deposition rates (${\dot M}_X$) estimated from the X-ray emission and the star formation rate estimated from the infrared luminosity. The star formation rates are 1/10 to 1/100 of the mass deposition rates suggesting that the re-heating of the ICM is generally very effective in reducing the amount of mass cooling from the hot phase but not eliminating it completely.

We present the second part of an H-band (1.6 microns) atlas of z<0.3 3CR radio galaxies, using the Hubble Space Telescope Near Infrared Camera and Multi-Object Spectrometer (HST NICMOS2). We present new imaging for 21 recently acquired sources, and host galaxy modeling for the full sample of 101 (including 11 archival) -- an 87% completion rate. Two different modeling techniques are applied, following those adopted by the galaxy morphology and the quasar host galaxy communities. Results are compared, and found to be in excellent agreement, although the former breaks down in the case of strongly nucleated sources. Companion sources are tabulated, and the presence of mergers, tidal features, dust disks and jets are catalogued. The tables form a catalogue for those interested in the structural and morphological dust-free host galaxy properties of the 3CR sample, and for comparison with morphological studies of quiescent galaxies and quasar host galaxies. Host galaxy masses are estimated, and found to typically lie at around 2*10^11 solar masses. In general, the population is found to be consistent with the local population of quiescent elliptical galaxies, but with a longer tail to low Sersic index, mainly consisting of low-redshift (z<0.1) and low-radio-power (FR I) sources. A few unusually disky FR II host galaxies are picked out for further discussion. Nearby external sources are identified in the majority of our images, many of which we argue are likely to be companion galaxies or merger remnants. The reduced NICMOS data are now publicly available from our website (http://archive.stsci.edu/prepds/3cr/), Comment: ApJS, 177, 148: Final version; includes revised figures 1, 15b, and section 7.5 (and other minor changes from editing process. 65 pages, inc. 17 figures

We report on an imaging survey with the Spitzer Space Telescope of 62 brightest cluster galaxies with optical line emission. These galaxies are located in the cores of X-ray luminous clusters selected from the ROSAT All-Sky Survey. We find that about half of these sources have a sign of excess infrared emission; 22 objects out of 62 are detected at 70 microns, 18 have 8 to 5.8 micron flux ratios above 1.0 and 28 have 24 to 8 micron flux ratios above 1.0. Altogether 35 of 62 objects in our survey exhibit at least one of these signs of infrared excess. Four galaxies with infrared excesses have a 4.5/3.6 micron flux ratio indicating the presence of hot dust, and/or an unresolved nucleus at 8 microns. Three of these have high measured [OIII](5007A)/Hbeta flux ratios suggesting that these four, Abell 1068, Abell 2146, and Zwicky 2089, and R0821+07, host dusty active galactic nuclei (AGNs). 9 objects (including the four hosting dusty AGNs) have infrared luminosities greater than 10^11 L_sol and so can be classified as luminous infrared galaxies (LIRGs). Excluding the four systems hosting dusty AGNs, the excess mid-infrared emission in the remaining brightest cluster galaxies is likely related to star formation., Comment: accepted for publication in ApJS

We observed the brightest central galaxy (BCG) in the nearby (z=0.0821) cool core galaxy cluster Abell 2597 with the IRAC and MIPS instruments on board the Spitzer Space Telescope. The BCG was clearly detected in all Spitzer bandpasses, including the 70 and 160 micron wavebands. We report aperture photometry of the BCG. The spectral energy distribution exhibits a clear excess in the FIR over a Rayleigh-Jeans stellar tail, indicating a star formation rate of ~4-5 solar masses per year, consistent with the estimates from the UV and its H-alpha luminosity. This large FIR luminosity is consistent with that of a starburst or a Luminous Infrared Galaxy (LIRG), but together with a very massive and old population of stars that dominate the energy output of the galaxy. If the dust is at one temperature, the ratio of 70 to 160 micron fluxes indicate that the dust emitting mid-IR in this source is somewhat hotter than the dust emitting mid-IR in two BCGs at higher-redshift (z~0.2-0.3) and higher FIR luminosities observed earlier by Spitzer, in clusters Abell 1835 and Zwicky 3146., Comment: Accepted at ApJ

We present deep emission-line imaging taken with the SOAR Optical Imaging Camera of the brightest cluster galaxy (BCG) in the nearby (z=0.035) X-ray cluster 2A0335+096. We analyze long-slit optical spectroscopy, archival VLA, Chandra X-ray, and XMM UV data. 2A0335+096 is a bright, cool-core X-ray cluster, once known as a cooling flow. Within the highly disturbed core revealed by Chandra X-ray observations, 2A0335+096 hosts a highly structured optical emission-line system. The redshift of the companion is within 100 km/s of the BCG and has certainly interacted with the BCG, and is likely bound to it. The comparison of optical and radio images shows curved filaments in H-alpha emission surrounding the resolved radio source. The velocity structure of the emission-line bar between the BCG nucleus and the companion galaxy provides strong evidence for an interaction between the two in the last ~50 Myrs. The age of the radio source is similar to the interaction time, so this interaction may have provoked an episode of radio activity. We estimate a star formation rate of >7 solar mass/yr based on the Halpha and archival UV data, a rate similar to, but somewhat lower than, the revised X-ray cooling rate of 10-30 solar masses/year estimated from XMM spectra by Peterson & workers. The Halpha nebula is limited to a region of high X-ray surface brightness and cool X-ray temperature. The detailed structures of H-alpha and X-ray gas differ. The peak of the X-ray emission is not the peak of H-alpha emission, nor does it lie in the BCG. The estimated age of the radio lobes and their interaction with the optical emission-line gas, the estimated timescale for depletion and accumulation of cold gas, and the dynamical time in the system are all similar, suggesting a common trigger mechanism., Comment: Accepted AJ, July 2007 publication. Vol 134, p. 14-25

In order to identify the dominant nuclear outflow mechanisms in Active Galactic Nuclei, we have undertaken deep, high resolution observations of two compact radio sources (PKS 1549-79 and PKS 1345+12) with the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope. Not only are these targets known to have powerful emission line outflows, but they also contain all the potential drivers for the outflows: relativistic jets, quasar nuclei and starbursts. ACS allows the compact nature (<0.15") of these radio sources to be optically resolved for the first time. Through comparison with existing radio maps we have seen consistency in the nuclear position angles of both the optical emission line and radio data. There is no evidence for bi-conical emission line features on the large-scale and there is a divergance in the relative position angles of the optical and radio structure. This enables us to exclude starburst driven outflows. However, we are unable to clearly distinguish between radiative AGN wind driven outflows and outflows powered by relativistic radio jets. The small scale bi-conical features, indicative of such mechanisms could be below the resolution limit of ACS, especially if aligned close to the line of sight. In addition, there may be offsets between the radio and optical nuclei induced by heavy dust obscuration, nebular continuum or scattered light from the AGN., Comment: 9 pages, 8 figures, emulateapj, ApJ Accepted

We present an overview of new HST imaging polarimetry of six nearby radio galaxies with optical jets. These observations triple the number of extragalactic jets with subarcsecond-resolution optical polarimetry. We discuss the polarization characteristics and optical morphology of each jet. We find evidence of high optical polarization, averaging 20%, but reaching upwards of $\sim 50%$ in some objects, confirming that the optical emission is synchrotron, and that the components of the magnetic fields perpendicular to the line of sight are well ordered. We find a wide range of polarization morphologies, with each jet having a somewhat different relationship between total intensity and polarized flux and the polarization position angle. We find two trends in all of these jets. First, jet ``edges'' are very often associated with high fractional optical polarizations, as also found in earlier radio observations of these and other radio jets. In these regions, the magnetic field vectors appear to track the jet direction, even at bends, where we see particularly high fractional polarizations. This indicates a strong link between the local magnetic field and jet dynamics. Second, optical flux maximum regions are usually well separated from maxima in fractional polarization and often are associated with polarization minima. This trend is not found in radio data and was found in our optical polarimetry of M87 with HST. However, unlike in M87, we do not find a general trend for near-90$^\circ$ rotations in the optical polarization vectors near flux maxima. We discuss possibilities for interpreting these trends, as well as implications for jet dynamics, magnetic field structure and particle acceleration., Comment: 32 pages, 8 figures, 1 reduced. ApJ, in press

We are conducting a large observing program with the Spitzer Space Telescope to determine the mid-to-far infrared spectral energy distributions of a well-defined sample of 87 nearby, 12 micron-selected Seyfert galaxies. In this paper we present the results of IRS low-resolution spectroscopy of a statistically representative subsample of 51 of the galaxies (59%), with an analysis of the continuum shapes and a comparison of the Seyfert types. We find that the spectra clearly divide into groups based on their continuum shapes and spectral features. Some spectral features are clearly related to a starburst contribution to the IR spectrum, while the mechanisms producing observed power-law continuum shapes, attributed to an AGN component, may be dust or non-thermal emission. The infrared spectral types appear to be related to the Seyfert types. Principal component analysis results suggest that the relative contribution of starburst emission may be the dominant cause of variance in the observed spectra. We find that Sy 1's have higher ratios of IR/radio emission than Sy 2's, as predicted by the unified model if the torus is optically thick in the mid-IR. However, smooth-density torus models predict a much larger difference between type 1's and 2's than observed in our sample. Our observations may be consistent with clumpy torus models containing a steep radial distribution of optically thick dense clumps. (Abridged), Comment: 39 pages, 17 figures, accepted to AJ, version with high-resolution figures available at http://www.cis.rit.edu/~clbsps/papers/buchanan.pdf

We present the results of IRS low-resolution spectroscopy of 51 Seyfert galaxies, part of a large Spitzer observing program to determine the mid-to-far infrared spectral energy distributions of a well-defined sample of 87 nearby, 12 micron-selected Seyferts. We find that the spectra clearly divide into groups based on their continuum shapes and spectral features. The infrared spectral types appear to be related to the Seyfert types. Some features are clearly related to a starburst contribution to the IR spectrum, while the observed power-law continuum shapes, attributed to the AGN, may be dust or non-thermal emission. Principal component analysis results suggest that the relative contribution of starburst emission is the dominant cause of variance in the spectra. We find that the Sy 2's show on average stronger starburst contributions than the Sy 1's., Comment: 6 pages, 5 figures, uses ASP style file. To appear in the proceedings of "The Spitzer Science Center 2005 Conference: Infrared Diagnostics of Galaxy Evolution", held in Pasadena, November 2005

We present newly acquired images of the near-infrared counterpart of 3CR radio sources. All the sources were selected to have a redshift of less than 0.3 to allow us to obtain the highest spatial resolution. The observations were carried out as a snapshot program using the Near-Infrared Camera and Multiobject Spectrograph (NICMOS) on-board the Hubble Space Telescope (HST). In this paper we describe 69 radio galaxies observed for the first time with NICMOS during HST cycle 13. All the objects presented here are elliptical galaxies. However, each of them has unique characteristics such as close companions, dust lanes, unresolved nuclei, arc-like features, globular clusters and jets clearly visible from the images or with basic galaxy subtraction., Comment: Submitted to ApJS, this version incorporates referee comments; 77 bitmapped figures

We report the discovery of a new optical-IR synchrotron jet in the radio galaxy 3C133 from our HST/NICMOS snapshot survey. The jet and eastern hotspot are well resolved, and visible at both optical and IR wavelengths. The IR jet follows the morphology of the inner part of the radio jet, with three distinct knots identified with features in the radio. The radio-IR SED's of the knots are examined, along with those of two more distant hotspots at the eastern extreme of the radio feature. The detected emission appears to be synchrotron, with peaks in the NIR for all except one case, which exhibits a power-law spectrum throughout., Comment: ApJ accepted. 14 pages, 6 figures

We present the new HST near-infrared polarimetry, broad and narrow-band imaging, and MERLIN 4.5GHz Multi-Frequency Synthesis radio imaging of 3C 293, a unique radio galaxy whose host is an obvious merger remnant, in an exceptionally under-dense region of space. We have discovered near-infrared, optical, and ultra-violet synchrotron emission from the jet. In the optical, the jet is mostly obscured by a dust lane, but three knots are clear in our HST NICMOS images at 1.6 and 2.0 microns, clearly aligning with features in the radio. The outer jet knot is highly polarized (~15%) at 2 microns, confirming the synchrotron emission mechanism. The radio-IR spectral index steepens significantly with distance from the nucleus, as in 3C 273 and in contrast to M 87. The inner knot is visible (with hindsight) on the WFPC2 and STIS images obtained for the earlier 3CR HST snapshot surveys. There is no [Fe II] emission seen associated with the jet, constraining the role of shock-induced ionisation by the jet. Overall there is a strong implication that the NIR jet emission is indeed synchrotron. From our NIR images, the core of the galaxy is clearly identifiable with the main feature in the western extension of the radio ``jet'' image, although no unresolved AGN component is identifiable even at K-band, consistent with an FRII-like nucleus obscured by an optically thick torus. The galaxy appears to have a single nucleus, with any multiple nuclei falling within the central

Among radio galaxies containing nuclear dust disks, the bipolar jet axis is generally observed to be perpendicular to the disk major axis. The FR I radio source 3C 449 is an outlier to this statistical majority, as it possesses a nearly parallel jet/disk orientation on the sky. We examine the 600 pc dusty disk in this galaxy with images from the Hubble Space Telescope. We find that a colormap of the disk exhibits a twist in its isocolor contours (isochromes). We model the colormap by integrating galactic starlight through an absorptive disk, and find that the anomalous twist in the isochromes can be reproduced in the model with a vertically thin, warped disk. The model predicts that the disk is nearly perpendicular to the jet axis within 100 pc of the nucleus. We discuss physical mechanisms capable of causing such a warp. We show that precessional models or a torque on the disk arising from a possible binary black hole in the AGN causes precession on a timescale that is too long to account for the predicted disk morphology. However, we estimate that the pressure in the X-ray emitting interstellar medium is large enough to perturb the disk, and argue that jet-driven anisotropy in the excited ISM may be the cause of the warp. In this way, the warped disk in 3C 449 may be a new manifestation of feedback from an active galactic nucleus., Comment: Accepted for publication in ApJ

Radio lobes inflated by active galactic nuclei at the centers of clusters are a promising candidate for halting condensation in clusters with short central cooling times because they are common in such clusters. In order to test the AGN-heating hypothesis, we obtained Chandra observations of two clusters with short central cooling times yet no evidence for AGN activity: Abell 1650 and Abell 2244. The cores of these clusters indeed appear systematically different from cores with more prominent radio emission. They do not have significant central temperature gradients, and their central entropy levels are markedly higher than in clusters with stronger radio emission, corresponding to central cooling times ~ 1 Gigayear. Also, there is no evidence for fossil X-ray cavities produced by an earlier episode of AGN heating. We suggest that either (1) the central gas has not yet cooled to the point at which feedback is necessary to prevent it from condensing, possibly because it is conductively stabilized, or (2) the gas experienced a major heating event $\gtrsim 1$ Gyr in the past and has not required feedback since then. The fact that these clusters with no evident feedback have higher central entropy and therefore longer central cooling times than clusters with obvious AGN feedback strongly suggests that AGNs supply the feedback necessary to suppress condensation in clusters with short central cooling times., Comment: ApJ Letter, in press

We present the results of HST/WFPC2 medium and narrow band imaging and VLA and MERLIN2 radio imaging of three powerful radio galaxies: 3C 171, 3C 277.3, and PKS 2250-41. We obtained images of the rest frame [OIII]$\lambda$5007 and [OII]$\lambda$3727 line emission using the Linear Ramp Filters on WFPC2. The correlations between the emission line morphology and the [OIII]/[OII] line ratios with the radio emission seen in ground based observations are clarified by the HST imaging. We confirm that the radio lobes and hot-spots are preferentially associated with lower ionization gas. 3C 171 exhibits high surface brightness emission line gas mainly along the radio source axis. The lowest ionization gas is seen at the Eastern hot spot. In 3C 277.3 there is bright high ionization gas (and continuum) offset just to the east of the radio knot K1. Our observations are consistent with previous work suggesting that this emission is produced by precursor gas ionized by the shock being driven into the cloud by the deflected radio jet. In PKS 2250-41 we resolve the emission line arc which wraps around the outer rim of the western lobe. The lower ionization [OII] emission is nested just interior to the higher ionization [OIII] emission suggesting that we have resolved the cooling region behind the bow shock. We also detect possible continuum emission from the secondary hot-spot. Thus, our observations support the hypothesis that in these sources, the interaction between the expanding radio source and the ambient gas strongly influences the morphology, kinematics, and ionization of the gas., Comment: 27 pages, includes 12 figures and 6 tables, accepted for publication in The Astronomical Journal

We present BIMA 3 mm observations of the radio continuum source and the molecular gas disk in the radio loud Fanaroff & Riley Type I (FR I)galaxy NGC3801.We have detected a continuum source in the nucleus and determined that it has a flat millimeter-wave spectrum, suggesting that the emission is non-thermal and due to an AGN; the radio core is not evident in existing VLA observations. We also map the extended 3 mm emission from the previously known radio jets. In addition, we detect CO (1--0) emission associated with the dust disk observed in previous HST images. A velocity gradient is observed, indicating a two kpc radius rotating gas ring or disk oriented roughly perpendicular to the radio jets. The inferred molecular gas mass of the disk is $M(H_{2})=3\times10^{8}M_{\odot}$, about 1% of the dynamical mass. We also find a $\sim 10^8$ M$_\odot$ molecular gas clump not associated with the gas disk. There is evidence that this gas is associated with a merger and is infalling. This suggests that FR I type activity is related to merger activity, as is thought to be the case for FR II type radio galaxies. We also find indications that one of the radio jets is entraining gas from the infalling molecular gas., Comment: To appear in ApJ, August 2005

We present a Hubble Space Telescope image of the FRII radio galaxy 3C 401, obtained at 1.6 microns with the NICMOS camera in which we identify the infrared counterpart of the brightest region of the radio jet. The jet has a complex radio structure and brightens where bending occurs, most likely as a result of relativistic beaming. We analyze archival data in the radio, optical and X-ray bands and we derive its spectral energy distribution. Differently from all of the previously known optical extragalactic jets, the jet in 3C401 is not detected in the X-rays even in a long 48ksec X-ray Chandra exposure and the infrared emission dominates the overall SED. We propose that the dominant radiation mechanism of this jet is synchrotron. The low X-ray emission is then caused by two different effects: i) the lack of any strong external photon field and ii) the shape of the electron distribution. This affects the location of the synchrotron peak in the SED, resulting in a sharp cut-off at energies lower than the X-rays. Thus 3C401 shows a new type of jet which has intermediate spectral properties between those of FRI, which are dominated by synchrotron emission up to X-ray energies, and FRII/QSO, which show a strong high-energy emission due to inverse-Compton scattering of external photons. This might be a clue for the presence of a continuous ``sequence'' in the properties of large scale jets, in analogy with the ``blazar sequence'' already proposed for sub-pc scale jets., Comment: 7 pages, 5 figures, accepted for publication in the ApJ. A typo in the value of the radio to X-ray spectral index (sect. 5, paragraph 2, thanks to T. Cheung) and other minor corrections have been made

We present sensitive, multifrequency Very Long Baseline Array (VLBA) images of the nuclear radio sources of NGC 1068. At 5 GHz and 8.4 GHz, the radio continuum source S1, argued to mark the location of the hidden active nucleus, resolves into an elongated, ~ 0.8 pc source oriented nearly at right angles to the radio jet axis but more closely aligned to the distribution of the nuclear H2O maser spots. S1 is detected at 5 GHz but not at 1.4 GHz, indicating strong free-free absorption below 5 GHz, and it has a flat spectrum between 5 GHz and 8.4 GHz. A 5/8.4 GHz spectral index map reveals an unresolved, inverted spectrum source at the center of the S1 structure which may mark the AGN proper. The average brightness temperature is too low for synchrotron self-absorption to impact the integrated spectrum significantly. In addition, a careful registration with the nuclear H2O masers argues that the S1 continuum source arises from the inner regions of the maser disk rather than a radio jet. The emission mechanism may be direct, thermal free-free emission from an X-ray-heated corona or wind arising from the molecular disk. We demonstrate that the hidden active nucleus is sufficiently luminous, to within the current estimates, to provide the requisite heating., Comment: 38 pages, PDF prepared with MS Word / Acrobat, accepted for publication in ApJ

We present the results of the analysis of a set of medium resolution spectra, obtained by the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope, of the emission line gas present in the nuclei of a complete sample of 21 nearby, early-type galaxies with radio jets (the UGC FR-I Sample). For each galaxy nucleus we present spectroscopic data in the region of H-alpha and the dervived kinematics. We find that in 67% of the nuclei the gas appears to be rotating and, with one exception, the cases where rotation is not seen are either face on or have complex central morphologies. We find that in 62% of the nuclei the fit to the central spectrum is improved by the inclusion of a broad component. The broad components have a mean velocity dispersion of 1349 +/- 345 km\s and are redshifted from the narrow line components (assuming an origin in H-alpha) by 486 +/- 443 km\s., Comment: 119 pages, 26 figures, ApJS Accepted, version with full figures available at http://www.astro.columbia.edu/~jake/pub/fr1datapaper.pdf

We used the HST WFPC2 to obtain I-band images of the centers of 81 brightest cluster galaxies (BCGs), drawn from a volume-limited sample of nearby BCGs. The images show a rich variety of morphological features, including multiple or double nuclei, dust, stellar disks, point source nuclei, and central surface brightness depressions. High resolution surface brightness profiles could be inferred for 60 galaxies. Of those, 88% have well-resolved cores. Twelve percent of the BCG sample lacks a well-resolved core; all but one of these BCGs have ``power-law'' profiles. Some of these galaxies have higher luminosities than any power-law galaxy identified by Faber et al. (1997), and have physical upper limits on the break radius well below the values observed for core galaxies of the same luminosity. These results support the idea that the central structure of early-type galaxies is bimodal in its physical properties, but also suggest that there exist high luminosity galaxies with power-law profiles (or unusually small cores). The BCGs in the latter category tend to fall at the low end of the BCG luminosity function and tend to have low values of the quantity alpha (the logarithmic slope of the metric luminosity as a function of radius, at 10 kpc). Since theoretical calculations have shown that the luminosities and alpha values of BCGs grow with time as a result of accretion, this suggests a scenario in which elliptical galaxies evolve from power-law profiles to core profiles through accretion and merging. This is consistent with theoretical scenarios that invoke the formation of massive black hole binaries during merger events (Abridged)., Comment: Accepted for publication in AJ. Postscript file with high resolution figures (1 and 3) is available at http://spider.ipac.caltech.edu/staff/seppo/bcgfigs

In this paper we analyze the relation between radio, optical continuum and Halpha+[NII] emission from the cores of a sample of 21 nearby Fanaroff & Riley type I galaxies as observed with the VLBA and HST. The emission arises inside the inner tens of parsec of the galaxies. Core radio emission is observed in 19/20 galaxies, optical core continuum emission is detected in 12/21 galaxies and Halpha+[NII] core emission is detected in 20/21 galaxies. We confirm the recently detected linear correlation between radio and optical core emission in FR I galaxies and show that both core emissions also correlate with central Halpha+[NII] emission. The tight correlations between radio, optical and Halpha+[NII] core emission constrain the bulk Lorentz factor to gamma ~ 2-5 and gamma =< 2 for a continuous jet and a jet consisting of discrete blobs, respectively, assuming jet viewing angles in the range [30deg,90deg]. Radio and optical core emissions are likely to be synchrotron radiation from the inner jet, possibly with a significant contribution from emission by an accretion disk and/or flow. Elliptical galaxies with LINER nuclei without large-scale radio jets seem to follow the core emission correlations found in FR I galaxies. This suggests that the central engines could be very similar for the two classes of AGNs., Comment: 46 pages, 10 figures and 8 tables. Accepted for publication in AJ. Corrected few typos in text

We present ultraviolet images of 27 3CR radio galaxies with redshifts z<0.1 that have been imaged with the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope (HST). The observations employed the NUV-MAMA and broad-band filters with peak sensitivity at 2200A . We find that the UV luminosities show approximately a factor of 10 to 100 higher dispersion than the optical. We compare the UV morphologies with optical V- and R-band WFPC2 snapshot survey images. We have found dramatic, complex and extended ultraviolet emission from radio galaxies even at zero redshift. We find a diverse range of UV morphologies, some completely divergent from their visual morphology, which are reminiscent of the chaotic high-z radio galaxies structures seen in rest-frame UV. The UV morphologies show regions of star formation, jets, and possible scattered AGN continuum. The UV emission is generally not aligned with the radio structure. We also detect the diffuse UV emission of the host galaxy. We propose that these are the same physical phenomena as observed at high redshift, but on a smaller spatial scale., Comment: Version with high resolution images available at http://www.stsci.edu/~mga/pp.pdf Accepted ApJ supplements November 21, 2001

We discuss the properties of compact nuclear radio components in Seyfert galaxies from the extended 12-micron AGN sample of Rush et al.(1993). Our main results can be summarised as follows. Type 1 and type 2 Seyferts produce compact radio components which are indistinguishable in strength and aspect, indicating that their central engines are alike as proposed by the unification model. Infrared IRAS fluxes are more closely correlated with low-resolution radio fluxes than high-resolution radio fluxes, suggesting that they are dominated by kiloparsec-scale, extra-nuclear emission regions; extra-nuclear emission may be stronger in type 2 Seyferts. Early-type Seyfert galaxies tend to have stronger nuclear radio emission than late-type Seyfert galaxies. V-shaped extended emission-line regions, indicative of `ionisation cones', are usually found in sources with large, collimated radio outflows. Hidden broad lines are most likely to be found in sources with powerful nuclear radio sources. Type 1 and type 2 Seyferts selected by their IRAS 12-micron flux densities have well matched properties.

We present new HST STIS/MAMA near-UV images and archival WFPC2 V and R band images which reveal the presence of four star forming regions in an arc along the edge of the dust lane in the giant (4 Mpc) radio galaxy 3C236. Two of the star forming regions are relatively young with ages of order 1E7 yr, while the other two are older with ages of order 1E8 - 1E9 yr which is comparable to the estimated age of the giant radio source. Based on dynamical and spectral aging arguments, we suggest that the fuel supply to the AGN was interrupted for 1E7 yr and has now been restored, resulting in the formation of the inner 2 kpc scale radio source. This time scale is similar to that of the age of the youngest of the star forming regions. We suggest that the transport of gas in the disk is non-steady and that this produces both the multiple episodes of star formation in the disk as well as the multiple epochs of radio source activity. If the inner radio source and the youngest star forming region are related by the same event of gas transport, the gas must be transported from the hundreds of pc scale to the sub-parsec scale on a time scale of 1E7 yr, which is similar to the dynamical time scale of the gas on the hundreds of pc scales.

We observed 17 nearby low luminosity FR I radio galaxies using the NRAO Very Long Baseline Array (VLBA) at 1.67 GHz, as part of a multi-wavelength study of a complete sample of 21 sources selected by radio flux density from the Uppsala General Catalogue of Galaxies. We detected radio emission from all 17 galaxies. At a FWHM resolution of ~ 10x4 mas, five galaxies show only an unresolved radio core, ten galaxies show core-jet structures, and two galaxies show twin-jet structures. Comparing these VLBA images with images previously obtained with the NRAO VLA, we find that all detected VLBA jets are well-aligned on parsec scales with the VLA jets on kilo parsec scales, and that the jet to counter-jet surface brightness ratios, or the sidedness, decreases systematically with increasing distance along the jet. We attribute the sidedness to the Doppler boosting effect and its decline to the deceleration of the jets. We show that a distribution of Lorentz factor centered near Gamma = 5 can reproduce our VLBA detection statistics for core, core-jet and twin-jet sources. We also note that the luminosity per unit length, L_j, of the VLBA jets drops quickly with distance, r, along the jet, approximately as L_j being proportional to r^(-2.0). We discuss three different mechanisms to explain this jet fading: (1) the decrease of Doppler boosting due to jet deceleration, (2) synchrotron losses, and (3) expansion losses in constant velocity but adiabatically spreading jets. Mechanisms (1) and (2) are inconsistent with the observations, while mechanism (3) is consistent with the observations provided the magnetic field lines in the jets are aligned perpendicular to the jet axis. This implies that the deceleration of the jets required by the Unified Scheme does not occur on the tens of parsec scales, but must occur on larger scales., Comment: 32 pages, 7 figures

We present 8.4 GHz VLA A-configuration observations of 87 sources from the mid-infrared-selected AGN sample of Rush et al (1993). These 0.25 arcsec resolution observations allow elongated radio structures tens of parsecs in size to be resolved and enable radio components smaller than 3.5 arcsec to be isolated from diffuse galactic disc emission. When combined with previous data, matched radio observations covering ninety percent of the sample have been collected and these represent the largest sub-arcsecond-resolution radio imaging survey of a homogeneously-selected sample of Seyfert galaxies to date. We use our observations to identify 5 radio-loud AGN in the sample. The nature of the radio emission from Seyfert nuclei will be discussed in subsequent papers., Comment: (17 pages) For full image quality and ascii data lists see http://www.jb.man.ac.uk/~ahct/database.html

Active galactic nuclei (AGN) are generally accepted to be powered by the release of gravitational energy in a compact accretion disk surrounding a massive black hole. Such disks are also necessary to collimate powerful radio jets seen in some AGN. The unifying classification schemes for AGN further propose that differences in their appearance can be attributed to the opacity of the accreting material, which may obstruct our view of the central region of some systems. The popular model for the obscuring medium is a parsec-scale disk of dense molecular gas, although evidence for such disks has been mostly indirect, as their angular size is much smaller than the resolution of conventional telescopes. Here we report the first direct images of a pc-scale disk of ionised gas within the nucleus of NGC 1068, the archetype of obscured AGN. The disk is viewed nearly edge-on, and individual clouds within the ionised disk are opaque to high-energy radiation, consistent with the unifying classification scheme. In projection, the disk and AGN axes align, from which we infer that the ionised gas disk traces the outer regions of the long-sought inner accretion disk., Comment: 14 pages, LaTeX, PSfig, to appear in Nature. also available at http://hethp.mpe-garching.mpg.de/Preprints