The First Spatially Resolved Mid-IR Spectra of NGC 1068 Obtained at Diffraction-limited Resolution with LWS at Keck I Telescope

We present spatially resolved mid-IR spectra of NGC 1068 with a diffraction-limited resolution of 0.25\arcsec using the Long Wavelength Spectrometer (LWS) at the Keck I telescope. The mid-infrared image of NGC 1068 is extended along the N-S direction. Previous imaging studies have shown the extended regions are located inside the ionization cones indicating that the mid-infrared emission arises perhaps from the inner regions of the narrow-line clouds instead of the proposed dusty torus itself. The spatially resolved mid-IR spectra were obtained at two different slit position angles, +8.0 and -13.0 degrees across the elongated regions in the mid-IR. From these spectra, we found only weak silicate absorption toward the northern extended regions but strong in the nucleus and the southern extended regions. This is consistent with a model of a slightly inclined cold obscuring torus which covers much of the southern regions but is behind the northern extension. While a detailed analysis of the spectra requires a radiative transfer model, the lack of silicate emission from the northern extended regions prompts us to consider a dual dust population model as one of the possible explanations in which a different dust population exists in the ionization cones compared to that in the dusty torus. Dust inside the ionization cones may lack small silicate grains giving rise to only a featureless continuum in the northern extended regions while dust in the dusty torus has plenty of small silicate grains to produce the strong silicate absorption lines towards the nucleus and the southern extended regions.Comment: Accepted to Ap

Understanding the 8 micron vs. Pa-alpha relationship on sub-arcsecond scales in Luminous Infrared Galaxies

This work explores in detail the relation between the 8 micron and the Pa-alpha emissions for 122 HII regions identified in a sample of 10 low-z LIRGs with nearly constant metallicity (12 + log (O/H) ~ 8.8). We use Gemini/T-ReCS high-spatial resolution (<~ 0.4" ~ 120 pc for the average distance of 60 Mpc of our sample) mid-infrared imaging (at 8.7 micron or 10.3 micron) together with HST/NICMOS continuum and Pa-alpha images. The LIRG HII regions extend the L_8micron vs. L_Pa-alpha relation found for HII knots in the high-metallicity SINGS galaxies by about two orders of magnitude to higher luminosities. Since the metallicity of the LIRG sample is nearly constant, we can rule out this effect as a cause for the scatter seen in the relationship. In turn, it is attributed to two effects: age and PAH features. The L_8micron/L_Pa-alpha ratio, which varies by a factor of ten for the LIRG HII regions, is reproduced by a model with instantaneous star formation and ages ranging from ~ 4 to 7.5 Myr. The remaining dispersion around the model predictions for a given age is probably due to differential contributions of the PAH features (the 8.6 micron, in our case) to the 8 micron emission from galaxy to galaxy.Comment: 16 pages, 9 figures, accepted for publication in ApJ; paper with full-resolution figures can be found at: http://damir.iem.csic.es/extragalactic

Resolved Mid-IR Emission in the Narrow Line Region of NGC 4151

We present subarcsecond resolution mid infrared images of NGC 4151 at 10.8 micron and 18.2 micron. These images were taken with the University of Florida mid-IR camera/spectrometer OSCIR at the Gemini North 8-m telescope. We resolve emission at both 10.8 micron and 18.2 micron extending ~ 3.5" across at a P.A. of ~ 60 degrees. This coincides with the the narrow line region of NGC 4151 as observed in [OIII] by the Hubble Space Telescope. The most likely explanation for this extended mid-IR emission is dust in the narrow line region heated by a central engine. We find no extended emission associated with the proposed torus and place an upper limit on its mid-IR size of less than or equal to ~ 35 pc.Comment: accepted for publication in the Astrophysical Journal, 19 pages including 5 figure

The Mid-Infrared Emission of M87

We discuss Subaru and Spitzer Space Telescope imaging and spectroscopy of M87 in the mid-infrared from 5-35 um. These observations allow us to investigate mid-IR emission mechanisms in the core of M87 and to establish that the flaring, variable jet component HST-1 is not a major contributor to the mid-IR flux. The Spitzer data include a high signal-to-noise 15-35 $\mu$m spectrum of the knot A/B complex in the jet, which is consistent with synchrotron emission. However, a synchrotron model cannot account for the observed {\it nuclear} spectrum, even when contributions from the jet, necessary due to the degrading of resolution with wavelength, are included. The Spitzer data show a clear excess in the spectrum of the nucleus at wavelengths longer than 25 um, which we model as thermal emission from cool dust at a characteristic temperature of 55 \pm 10 K, with an IR luminosity \sim 10^{39} {\rm ~erg ~s^{-1}}. Given Spitzer's few-arcsecond angular resolution, the dust seen in the nuclear spectrum could be located anywhere within ~5'' (390 pc) of the nucleus. In any case, the ratio of AGN thermal to bolometric luminosity indicates that M87 does not contain the IR-bright torus that classical unified AGN schemes invoke. However, this result is consistent with theoretical predictions for low-luminosity AGNsComment: 9 pages, 7 figures, ApJ, in pres

Trajectory mapping: A tool for validation of trace gas observations

We investigate the effectiveness of trajectory mapping(TM) as a data validation tool. TM combines a dynamical model of the atmosphere with trace gas observations to provide more statistically robust estimates of instrument performance over much broader geographic areas than traditional techniques are able to provide. We present four detailed case studies selected so that the traditional techniques are expected to work well. In each case the TM results are equivalent to or improve upon the measurement comparisons performed with traditional approaches. The TM results are statistically more robust than those achieved using traditional approaches since the TM comparisons occur over a much larger range of geophysical variability. In the first case study we compare ozone data from the Halogen Occultation Experiment (HALOE) with Microwave Limb Sounder(MLS). TM comparisons appear to introduce little to no error as compared to the traditional approach. In the second case study we compare ozone data from HALOE with that from the Stratospheric Aerosol and Gas Experiment TT(SAGE TT). TM results in differences of less than 5% as compared to the traditional approach at altitudes between 18 and 25 km and less than 10% at altitudes between 25 and 40 km.In the third case study we show that ozone profiles generated from HALOE data using TM compare well with profiles from five European ozonesondes. In the fourth case study we evaluate the precision of MLS H20 using TM and find typical precision uncertainties of 3-7% at most latitudes and altitudes. The TM results agree well with previous estimates but are the result of a global analysis of the data rather than an analysis in the limited latitude bands in which traditional approaches work. Finally, sensitivity studies using the MLS H20 data show the following: (1) a combination of forward and backward trajectory calculations minimize uncertainties in isentropic TM; (2) although the uncertainty of the technique increases with trajectory duration,TM calculations of up to 14 days can provide reliable information for use in data validation studies; (3) a correlation coincidence criterion of 400 km produces the best TM results under most circumstances; (4) TM performs well compared to (and sometimes better than) traditional approaches at all latitudes and in most seasons and; (5) TM introduces no statistically significant biases at altitudes between 22 and 40 km

Inference on periodicity of circadian time series

Estimation of the period length of time-course data from cyclical biological processes, such as those driven by the circadian pacemaker, is crucial for inferring the properties of the biological clock found in many living organisms. We propose a methodology for period estimation based on spectrum resampling (SR) techniques. Simulation studies show that SR is superior and more robust to non-sinusoidal and noisy cycles than a currently used routine based on Fourier approximations. In addition, a simple fit to the oscillations using linear least squares is available, together with a non-parametric test for detecting changes in period length which allows for period estimates with different variances, as frequently encountered in practice. The proposed methods are motivated by and applied to various data examples from chronobiology

The Angular Correlation Function of Galaxies from Early SDSS Data

The Sloan Digital Sky Survey is one of the first multicolor photometric and spectroscopic surveys designed to measure the statistical properties of galaxies within the local Universe. In this Letter we present some of the initial results on the angular 2-point correlation function measured from the early SDSS galaxy data. The form of the correlation function, over the magnitude interval 18<r*<22, is shown to be consistent with results from existing wide-field, photographic-based surveys and narrower CCD galaxy surveys. On scales between 1 arcminute and 1 degree the correlation function is well described by a power-law with an exponent of ~ -0.7. The amplitude of the correlation function, within this angular interval, decreases with fainter magnitudes in good agreement with analyses from existing galaxy surveys. There is a characteristic break in the correlation function on scales of approximately 1-2 degrees. On small scales, < 1', the SDSS correlation function does not appear to be consistent with the power-law form fitted to the 1'< theta <0.5 deg data. With a data set that is less than 2% of the full SDSS survey area, we have obtained high precision measurements of the power-law angular correlation function on angular scales 1' < theta < 1 deg, which are robust to systematic uncertainties. Because of the limited area and the highly correlated nature of the error covariance matrix, these initial results do not yet provide a definitive characterization of departures from the power-law form at smaller and larger angles. In the near future, however, the area of the SDSS imaging survey will be sufficient to allow detailed analysis of the small and large scale regimes, measurements of higher-order correlations, and studies of angular clustering as a function of redshift and galaxy type