Demonstrating the negligible contribution of optical ACS/HST galaxies to source-subtracted cosmic infrared background fluctuations in deep IRAC/Spitzer images

We study the possible contribution of optical galaxies detected with the {\it Hubble} ACS instrument to the near-IR cosmic infrared (CIB) fluctuations in deep {\it Spitzer} images. The {\it Spitzer} data used in this analysis are obtained in the course of the GOODS project from which we select four independent $10^\prime\times10^\prime$ regions observed at both 3.6 and 4.5 \um. ACS source catalogs for all of these areas are used to construct maps containing only their emissions in the ACS $B, V, i, z$-bands. We find that deep Spitzer data exhibit CIB fluctuations remaining after removal of foreground galaxies of a very different clustering pattern at both 3.6 and 4.5 \um than the ACS galaxies could contribute. We also find that there are very good correlations between the ACS galaxies and the {\it removed} galaxies in the Spitzer maps, but practically no correlations remain with the residual Spitzer maps used to identify the CIB fluctuations. These contributions become negligible on larger scales used to probe the CIB fluctuations arising from clustering. This means that the ACS galaxies cannot contribute to the large-scale CIB fluctuations found in the residual Spitzer data. The absence of their contributions also means that the CIB fluctuations arise at z\gsim 7.5 as the Lyman break of their sources must be redshifted past the longest ACS band, or the fluctuations have to originate in the more local but extremely low luminosity galaxies.Comment: Ap.J.Letters, in press. Minor revisions to mathc the accepted versio

Bow Shocks from Neutron Stars: Scaling Laws and HST Observations of the Guitar Nebula

The interaction of high-velocity neutron stars with the interstellar medium produces bow shock nebulae, where the relativistic neutron star wind is confined by ram pressure. We present multi-wavelength observations of the Guitar Nebula, including narrow-band H-alpha imaging with HST/WFPC2, which resolves the head of the bow shock. The HST observations are used to fit for the inclination of the pulsar velocity vector to the line of sight, and to determine the combination of spindown energy loss, velocity, and ambient density that sets the scale of the bow shock. We find that the velocity vector is most likely in the plane of the sky. We use the Guitar Nebula and other observed neutron star bow shocks to test scaling laws for their size and H-alpha emission, discuss their prevalence, and present criteria for their detectability in targeted searches. The set of H-alpha bow shocks shows remarkable consistency, in spite of the expected variation in ambient densities and orientations. Together, they support the assumption that a pulsar's spindown energy losses are carried away by a relativistic wind that is indistinguishable from being isotropic. Comparison of H-alpha bow shocks with X-ray and nonthermal, radio-synchrotron bow shocks produced by neutron stars indicates that the overall shape and scaling is consistent with the same physics. It also appears that nonthermal radio emission and H-alpha emission are mutually exclusive in the known objects and perhaps in all objects.Comment: 12 pages, 7 figures (3 degraded), submitted to ApJ; minor revisions and updates in response to referee report. (AASTeX, includes emulateapj5 and onecolfloat5.

Point Process Algorithm: A New Bayesian Approach for Planet Signal Extraction with the Terrestrial Planet Finder

The capability of the Terrestrial Planet Finder Interferometer (TPF-I) for planetary signal extraction, including both detection and spectral characterization, can be optimized by taking proper account of instrumental characteristics and astrophysical prior information. We have developed the Point Process Algorithm (PPA), a Bayesian technique for extracting planetary signals using the sine-chopped outputs of a dual nulling interferometer. It is so-called because it represents the system being observed as a set of points in a suitably-defined state space, thus providing a natural way of incorporating our prior knowledge of the compact nature of the targets of interest. It can also incorporate the spatial covariance of the exozodi as prior information which could help mitigate against false detections. Data at multiple wavelengths are used simultaneously, taking into account possible spectral variations of the planetary signals. Input parameters include the RMS measurement noise and the a priori probability of the presence of a planet. The output can be represented as an image of the intensity distribution on the sky, optimized for the detection of point sources. Previous approaches by others to the problem of planet detection for TPF-I have relied on the potentially non-robust identification of peaks in a "dirty" image, usually a correlation map. Tests with synthetic data suggest that the PPA provides greater sensitivity to faint sources than does the standard approach (correlation map + CLEAN), and will be a useful tool for optimizing the design of TPF-I.Comment: 17 pages, 6 figures. AJ in press (scheduled for Nov 2006

Surveys of Galaxy Clusters with the Sunyaev Zel'dovich Effect

We have created mock Sunyaev-Zel'dovich effect (SZE) surveys of galaxy clusters using high resolution N-body simulations. To the pure surveys we add `noise' contributions appropriate to instrument and primary CMB anisotropies. Applying various cluster finding strategies to these mock surveys we generate catalogues which can be compared to the known positions and masses of the clusters in the simulations. We thus show that the completeness and efficiency that can be achieved depend strongly on the frequency coverage, noise and beam characteristics of the instruments, as well as on the candidate threshold. We study the effects of matched filtering techniques on completeness, and bias. We suggest a gentler filtering method than matched filtering in single frequency analyses. We summarize the complications that arise when analyzing the SZE signal at a single frequency, and assess the limitations of such an analysis. Our results suggest that some sophistication is required when searching for `clusters' within an SZE map.Comment: 8 pages, 7 figure

The Origin of X-shaped Radio Galaxies: Clues from the Z-symmetric Secondary Lobes

Existing radio images of a few X-shaped radio galaxies reveal Z-symmetric morphologies in their weaker secondary lobes which cannot be naturally explained by either the galactic merger or radio-lobe backflow scenarios, the two dominant models for these X-shaped radio sources. We show that the merger picture can explain these morphologies provided one takes into account that, prior to the coalescence of their supermassive black holes, the smaller galaxy releases significant amounts of gas into the ISM of the dominant active galaxy. This rotating gas, whose angular momentum axis will typically not be aligned with the original jets, is likely to provide sufficient ram pressure at a distance ~10 kpc from the nucleus to bend the extant jets emerging from the central engine, thus producing a Z-symmetry in the pair of radio lobes. Once the two black holes have coalesced some 10^7 yr later, a rapid reorientation of the jets along a direction close to that of the orbital angular momentum of the swallowed galaxy relative to the primary galaxy would create the younger primary lobes of the X-shaped radio galaxy. This picture naturally explains why such sources typically have powers close to the FR I/II break. We suggest that purely Z-symmetric radio sources are often en route to coalescence and the concomitant emission of substantial gravitational radiation, while X-shaped ones have already merged and radiated.Comment: 12 pages, 1 compressed figure; accepted for publication in ApJ Letter

Bayesian Image Reconstruction Based on Voronoi Diagrams

We present a Bayesian Voronoi image reconstruction technique (VIR) for interferometric data. Bayesian analysis applied to the inverse problem allows us to derive the a-posteriori probability of a novel parameterization of interferometric images. We use a variable Voronoi diagram as our model in place of the usual fixed pixel grid. A quantization of the intensity field allows us to calculate the likelihood function and a-priori probabilities. The Voronoi image is optimized including the number of polygons as free parameters. We apply our algorithm to deconvolve simulated interferometric data. Residuals, restored images and chi^2 values are used to compare our reconstructions with fixed grid models. VIR has the advantage of modeling the image with few parameters, obtaining a better image from a Bayesian point of view.Comment: 27 pages, 10 figures, to be published in APJ, 672, 127

A CLEAN-based Method for Deconvolving Interstellar Pulse Broadening from Radio Pulses

Multipath propagation in the interstellar medium distorts radio pulses, an effect predominant for distant pulsars observed at low frequencies. Typically, broadened pulses are analyzed to determine the amount of propagation-induced pulse broadening, but with little interest in determining the undistorted pulse shapes. In this paper we develop and apply a method that recovers both the intrinsic pulse shape and the pulse broadening function that describes the scattering of an impulse. The method resembles the CLEAN algorithm used in synthesis imaging applications, although we search for the best pulse broadening function, and perform a true deconvolution to recover intrinsic pulse structre. As figures of merit to optimize the deconvolution, we use the positivity and symmetry of the deconvolved result along with the mean square residual and the number of points below a given threshold. Our method makes no prior assumptions about the intrinsic pulse shape and can be used for a range of scattering functions for the interstellar medium. It can therefore be applied to a wider variety of measured pulse shapes and degrees of scattering than the previous approaches. We apply the technique to both simulated data and data from Arecibo observations.Comment: 9 pages, 6 figures, Accepted for publication in the Astrophysical Journa

The Radio Quiescence of Active Galaxies with High Accretion Rates

We present 6 cm Very Large Array observations of the Greene & Ho (2004) sample of 19 low-mass active galaxies with high accretion rates. This is one of the only studies of a uniform sample of narrow-line Seyfert 1 (NLS1) galaxies with such high sensitivity and resolution. Although we detect only one source, the entire sample is very radio-quiet down to strong limits. GH10 was found to have a radio power of 8.5 x 10^21 W/Hz, and a ratio R = f(6 cm)/f(4400 A) of 2.8. The 3 sigma upper limits for the remaining nondetections correspond to radio powers from 3 x 10^20 to 8 x 10^21 W/Hz and 0.47 < R <9.9. Stacking all nondetections yields an even stronger upper limit of R < 0.27. An assessment of existing observations in the literature confirms our finding that NLS1s are consistently radio-quiet, with a radio-loud fraction of 0%-6%, which is significantly lower than the 10%-20% observed in the general quasar population. By analogy with stellar-mass black holes, we argue that AGNs undergo a state transition at L_bol/L_Edd~0.01. Below this value a radiatively inefficient accretion flow effectively drives an outflow, which disappears when the flow turns into an optically thick, geometrically thin disk, or a radiation pressure-dominated slim disk at still higher L_bol/L_Edd.Comment: To appear in ApJ; 8 pages, 3 figures; uses emulateapj5.st

Diffraction-limited near-IR imaging at Keck reveals asymmetric, time-variable nebula around carbon star CIT 6

We present multi-epoch, diffraction-limited images of the nebula around the carbon star CIT 6 at 2.2 microns and 3.1 microns from aperture masking on the Keck-I telescope. The near-IR nebula is resolved into two main components, an elongated, bright feature showing time-variable asymmetry and a fainter component about 60 milliarcseconds away with a cooler color temperature. These images were precisely registered (~35 milliarcseconds) with respect to recent visible images from the Hubble Space Telescope (Trammell et al. 2000), which showed a bipolar structure in scattered light. The dominant near-IR feature is associated with the northern lobe of this scattering nebula, and the multi-wavelength dataset can be understood in terms of a bipolar dust shell around CIT 6. Variability of the near-IR morphology is qualitatively consistent with previously observed changes in red polarization, caused by varying illumination geometry due to non-uniform dust production. The blue emission morphology and polarization properties can not be explained by the above model alone, but require the presence of a wide binary companion in the vicinity of the southern polar lobe. The physical mechanisms responsible for the breaking of spherical symmetry around extreme carbon stars, such as CIT 6 and IRC+10216, remain uncertain.Comment: 18 pages, 5 figures (one in color), to appear in the Astrophysical Journa

The distribution of H13CN in the circumstellar envelope around IRC+10216

H13CN J=8-7 sub-millimetre line emission produced in the circumstellar envelope around the extreme carbon star IRC+10216 has been imaged at sub-arcsecond angular resolution using the SMA. Supplemented by a detailed excitation analysis the average fractional abundance of H13CN in the inner wind (< 5E15 cm) is estimated to be about 4E-7, translating into a total HCN fractional abundance of 2E-5 using the isotopic ratio 12C/13C=50. Multi-transitional single-dish observations further requires the H13CN fractional abundance to remain more or less constant in the envelope out to a radius of about 4E16 cm, where the HCN molecules are effectively destroyed, most probably, by photodissociation. The large amount of HCN present in the inner wind provides effective line cooling that can dominate over that generated from CO line emission. It is also shown that great care needs to be taken in the radiative transfer modelling where non-local, and non-LTE, effects are important and where the radiation field from thermal dust grains plays a major role in exciting the HCN molecules. The amount of HCN present in the circumstellar envelope around IRC+10216 is consistent with predicted photospheric values based on equilibrium chemical models and indicates that any non-equilibrium chemistry occurring in the extended pulsating atmosphere has no drastic net effect on the fractional abundance of HCN molecules that enters the outer envelope. It further suggests that few HCN molecules are incorporated into dust grains.Comment: Accepted for publication in ApJ. 20 pages, 7 figure