Source author record

B. T. Draine

B. T. Draine appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

66works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

66 published item(s)

preprint2021arXiv

Observational Constraints on the Physical Properties of Interstellar Dust in the Post-Planck Era

We present a synthesis of the astronomical observations constraining the wavelength-dependent extinction, emission, and polarization from interstellar dust from UV to microwave wavelengths on diffuse Galactic sightlines. Representative solid phase abundances for those sightlines are also derived. Given the sensitive new observations of polarized dust emission provided by the Planck satellite, we place particular emphasis on dust polarimetry, including continuum polarized extinction, polarization in the carbonaceous and silicate spectroscopic features, the wavelength-dependent polarization fraction of the dust emission, and the connection between optical polarized extinction and far-infrared polarized emission. Together, these constitute a set of constraints that should be reproduced by models of dust in the diffuse interstellar medium.

preprint2020arXiv

Detection of PAH Absorption and Determination of the Mid-Infrared Diffuse Interstellar Extinction Curve from the Sightline Toward Cyg OB2-12

The sightline toward the luminous blue hypergiant Cyg OB2-12 is widely used in studying interstellar dust on account of its large extinction ($A_V \simeq 10$ mag) and the fact that this extinction appears to be dominated by dust typical of the diffuse interstellar medium. We present a new analysis of archival ISO-SWS and Spitzer IRS observations of Cyg OB2-12 using a model of the emission from the star and its stellar wind to determine the total extinction $A_λ$ from 2.4--37 $μ$m. In addition to the prominent 9.7 and 18 $μ$m silicate features, we robustly detect absorption features associated with polycyclic aromatic hydrocarbons (PAHs), including the first identification of the 7.7 $μ$m feature in absorption. The 3.3 $μ$m aromatic feature is found to be much broader in absorption than is typically seen in emission. The 3.4 and 6.85 $μ$m aliphatic hydrocarbon features are observed with relative strengths consistent with observation of these features on sightlines toward the Galactic Center. We identify and characterize more than sixty spectral lines in this wavelength range, which may be useful in constraining models of the star and its stellar wind. Based on this analysis, we present an extinction curve $A_λ/A_{2.2 μm}$ that extrapolates smoothly to determinations of the mean Galactic extinction curve at shorter wavelengths and to dust opacities inferred from emission at longer wavelengths, providing a new constraint on models of interstellar dust in the mid-infrared.

preprint2019arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.

preprint2016arXiv

Graphite Revisited

Laboratory measurements are used to constrain the dielectric tensor for graphite, from microwave to X-ray frequencies. The dielectric tensor is strongly anisotropic even at X-ray energies. The discrete dipole approximation is employed for accurate calculations of absorption and scattering by single-crystal graphite spheres and spheroids. For randomly-oriented single-crystal grains, the so-called 1/3 - 2/3 approximation for calculating absorption and scattering cross sections is exact in the limit a/lambda -> 0, provides better than ~10% accuracy in the optical and UV even when a/lambda is not small, but becomes increasingly inaccurate at infrared wavelengths, with errors as large as ~40% at lambda = 10 micron. For turbostratic graphite grains, the Bruggeman and Maxwell Garnett treatments yield similar cross sections in the optical and ultraviolet, but diverge in the infrared, with predicted cross sections differing by over an order of magnitude in the far-infrared. It is argued that the Maxwell Garnett estimate is likely to be more realistic, and is recommended. The out-of-plane lattice resonance of graphite near 11.5 micron may be observable in absorption with the MIRI spectrograph on JWST. Aligned graphite grains, if present in the ISM, could produce polarized X-ray absorption and polarized X-ray scattering near the carbon K edge.

preprint2016arXiv

Planck intermediate results. XXIX. All-sky dust modelling with Planck, IRAS, and WISE observations

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

preprint2015arXiv

A Case Against Spinning PAHs as the Source of the Anomalous Microwave Emission

We employ an all-sky map of the anomalous microwave emission (AME) produced by component separation of the microwave sky to study correlations between the AME and Galactic dust properties. We find that while the AME is highly correlated with all tracers of dust emission, the best predictor of the AME strength is the dust radiance. Fluctuations in the AME intensity per dust radiance are uncorrelated with fluctuations in the emission from polycyclic aromatic hydrocarbons (PAHs), casting doubt on the association between AME and PAHs. The PAH abundance is strongly correlated with the dust optical depth and dust radiance, consistent with PAH destruction in low density regions. We find that the AME intensity increases with increasing radiation field strength, at variance with predictions from the spinning dust hypothesis. Finally, the temperature-dependence of the AME per dust radiance disfavors the interpretation of the AME as thermal emission. A reconsideration of other AME carriers, such as ultrasmall silicates, and other emission mechanisms, such as magnetic dipole emission, is warranted.

preprint2015arXiv

Heating and cooling of the neutral ISM in the NGC4736 circumnuclear ring

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.

preprint2014arXiv

A far-IR view of the starburst driven superwind in NGC 2146

NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the central 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems.

preprint2014arXiv

Andromeda's Dust

Spitzer Space Telescope and Herschel Space Observatory imaging of M31 is used, with a physical dust model, to construct maps of dust surface density, dust-to-gas ratio, starlight heating intensity, and PAH abundance, out to R=25kpc. The global dust mass is M_d=5.4x10^7Msol, the global dust/H mass ratio is M_d/M_H=0.0081, and the global PAH abundance is <q_PAH>=0.039. The dust surface density has an inner ring at R=5.6kpc, a maximum at R=11.2kpc, and an outer ring at R=15.1kpc. The dust/gas ratio varies from M_d/M_H=0.026 at the center to ~0.0027 at R=25kpc. From the dust/gas ratio, we estimate the ISM metallicity to vary by a factor ~10, from Z/Zsol=3 at R=0 to ~0.3 at R=25kpc. The dust heating rate parameter <U> peaks at the center, with <U> approx 35, declining to <U> approx 0.25 at R=20kpc. Within the central kpc, the starlight heating intensity inferred from the dust modeling is close to what is estimated from the stars in the bulge. The PAH abundance reaches a peak q_PAH=0.045 at R=11.2kpc. When allowance is made for the different spectrum of the bulge stars, q_PAH for the dust in the central kpc is similar to the overall value of q_PAH in the disk. The silicate-graphite-PAH dust model used here is generally able to reproduce the observed dust spectral energy distribution across M31, but overpredicts 500um emission at R=2-6kpc, suggesting that at R=2-6kpc, the dust opacity varies more steeply with frequency (with beta approx 2.3 between 200 and 600um) than in the model

preprint2014arXiv

Cool dust heating and temperature mixing in nearby star-forming galaxies

Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields responsible for heating the dust. After fitting the dust and stellar radial profiles with exponential functions, we fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified black bodies using both variable (MBBV) and fixed (MBBF) emissivity indices beta, as well as with physically motivated dust models. Results show that while most SED parameters decrease with radius, the emissivity index beta also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average beta~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with beta<= 1. We find that flatter slopes (beta<= 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust-beta degeneracy. This trend is related to variations in Umin since beta and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat beta<=1. Both these results strongly suggest that the low apparent βvalues (flat slopes) in MBBV fits are caused by temperature mixing along the line-of-sight, rather than by intrinsic variations in grain properties. Abstract truncated for arXiv.

preprint2014arXiv

Light scattering by hexagonal columns in the discrete dipole approximation

Scattering by infinite hexagonal ice prisms is calculated using Maxwell's equations in the discrete dipole approximation for size parameters up to x=400. Birefringence is included in the calculations. Applicability of the geometric optics approximation is investigated. Excellent agreement between wave optics and geometric optics is observed for large size parameter in the outer part of the 22 degree halo feature. For smaller ice crystals halo broadening is predicted, and there is appreciable "spillover" of the halo into shadow scattering angles <22 degrees. Ways to retrieve ice crystal sizes are suggested based on the full width at half-maximum of the halo, the power at <22 deg, and the halo polarization.

preprint2014arXiv

Planck 2013 results. XI. All-sky model of thermal dust emission

This paper presents an all-sky model of dust emission from the Planck 857, 545 and 353 GHz, and IRAS 100 micron data. Using a modified black-body fit to the data we present all-sky maps of the dust optical depth, temperature, and spectral index over the 353-3000 GHz range. This model is a tight representation of the data at 5 arc min. It shows variations of the order of 30 % compared with the widely-used model of Finkbeiner, Davis, and Schlegel. The Planck data allow us to estimate the dust temperature uniformly over the whole sky, providing an improved estimate of the dust optical depth compared to previous all-sky dust model, especially in high-contrast molecular regions. An increase of the dust opacity at 353 GHz, tau_353/N_H, from the diffuse to the denser interstellar medium (ISM) is reported. It is associated with a decrease in the observed dust temperature, T_obs, that could be due at least in part to the increased dust opacity. We also report an excess of dust emission at HI column densities lower than 10^20 cm^-2 that could be the signature of dust in the warm ionized medium. In the diffuse ISM at high Galactic latitude, we report an anti-correlation between tau_353/N_H and T_obs while the dust specific luminosity, i.e., the total dust emission integrated over frequency (the radiance) per hydrogen atom, stays about constant. The implication is that in the diffuse high-latitude ISM tau_353 is not as reliable a tracer of dust column density as we conclude it is in molecular clouds where the correlation of tau_353 with dust extinction estimated using colour excess measurements on stars is strong. To estimate Galactic E(B-V) in extragalactic fields at high latitude we develop a new method based on the thermal dust radiance, instead of the dust optical depth, calibrated to E(B-V) using reddening measurements of quasars deduced from Sloan Digital Sky Survey data.

preprint2014arXiv

The Heating of Mid-Infrared Dust in the Nearby Galaxy M33: A Testbed for Tracing Galaxy Evolution

Because the 8 μm polycyclic aromatic hydrocarbon (PAH) emission has been found to correlate with other well-known star formation tracers, it has widely been used as a star formation rate (SFR) tracer. There are, however, studies that challenge the accuracy and reliability of the 8 μm emission as a SFR tracer. Our study, part of the Herschel M33 Extended Survey (HERM33ES) open time key program, aims at addressing this issue by analyzing the infrared emission from the nearby spiral galaxy M33 at the high spatial scale of 75 pc. Combining data from the Herschel Space Observatory and the Spitzer Space Telescope we find that the 8 μm emission is better correlated with the 250 μm emission, which traces cold interstellar gas, than with the 24 μm emission. The L(8)/L(24) ratio is highly depressed in 24 μm luminous regions, which correlate with known HII regions. We also compare our results with the dust emission models by Draine & Li (2007). We confirm that the depression of 8 μm PAH emission near star-forming regions is higher than what is predicted by models; this is possibly an effect of increased stellar radiation from young stars destroying the dust grains responsible for the 8 μm emission as already suggested by other authors. We find that the majority of the 8 μm emission is fully consistent with heating by the diffuse interstellar medium, similar to what recently determined for the dust emission in M31 by Draine at al. (2013). We also find that the fraction of 8 μm emission associated with the diffuse interstellar radiation field ranges between 60% and 80% and is 40% larger than the diffuse fraction at 24 μm.

preprint2013arXiv

Calibration of the total infrared luminosity of nearby galaxies from Spitzer and Herschel bands

We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved SED modelling of these objects using the Draine and Li (2007) dust models and investigate the influence of the addition of SPIRE measurements in the estimation of LTIR. We find that using data up to 250 um leads to local LTIR values consistent with those obtained with a complete coverage (up to 500 um) within 10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the (70-160 um) band. The (24-70 um) fraction decreases with increasing metallicity. The (160-1100 um) submillimeter band can account for up to 25 per cent of the LTIR in metal-rich galaxies. We investigate the correlation between TIR surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three PACS bands can be used as reliable monochromatic estimators of the LTIR, the 100 um band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250um and LTIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated LTIR of nearby galaxies as well as LTIR of z=1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.

preprint2013arXiv

Mapping dust through emission and absorption in nearby galaxies

Dust has long been identified as a barrier to measuring inherent galaxy properties. However, the link between dust and attenuation is not straightforward and depends on both the amount of dust and its distribution. Herschel imaging of nearby galaxies undertaken as part of the KINGFISH project allows us to map the dust as seen in emission with unprecedented sensitivity and ~1 kpc resolution. We present here new optical integral field unit spectroscopy for eight of these galaxies that provides complementary 100-200 pc scale maps of the dust attenuation through observation of the reddening in both the Balmer decrement and the stellar continuum. The stellar continuum reddening, which is systematically less than that observed in the Balmer decrement, shows no clear correlation with the dust, suggesting that the distribution of stellar reddening acts as a poor tracer of the overall dust content. The brightest HII regions are observed to be preferentially located in dusty regions, and we do find a correlation between the Balmer line reddening and the dust mass surface density for which we provide an empirical relation. Some of the high-inclination systems in our sample exhibit high extinction, but we also find evidence that unresolved variations in the dust distribution on scales smaller than 500 pc may contribute to the scatter in this relation. We caution against the use of integrated A_V measures to infer global dust properties.

preprint2013arXiv

Shock Excited Molecules in NGC 1266: ULIRG conditions at the center of a Bulge Dominated Galaxy

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

preprint2013arXiv

Star Formation Rates in Resolved Galaxies: Calibrations with Near and Far Infrared Data for NGC5055 and NGC6946

We use the near--infrared Brγhydrogen recombination line as a reference star formation rate (SFR) indicator to test the validity and establish the calibration of the {\it Herschel} PACS 70 μm emission as a SFR tracer for sub--galactic regions in external galaxies. Brγoffers the double advantage of directly tracing ionizing photons and of being relatively insensitive to the effects of dust attenuation. For our first experiment, we use archival CFHT Brγand Ks images of two nearby galaxies: NGC\,5055 and NGC\,6946, which are also part of the {\it Herschel} program KINGFISH (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel). We use the extinction corrected Brγemission to derive the SFR(70) calibration for H{\sc ii} regions in these two galaxies. A comparison of the SFR(70) calibrations at different spatial scales, from 200 pc to the size of the whole galaxy, reveals that about 50% of the total 70μm emission is due to dust heated by stellar populations that are unrelated to the current star formation. We use a simple model to qualitatively relate the increase of the SFR(70) calibration coefficient with decreasing region size to the star formation timescale. We provide a calibration for an unbiased SFR indicator that combines the observed Hαwith the 70 μm emission, also for use in H{\sc ii} regions. We briefly analyze the PACS 100 and 160 μm maps and find that longer wavelengths are not as good SFR indicators as 70μm, in agreement with previous results. We find that the calibrations show about 50% difference between the two galaxies, possibly due to effects of inclination.

preprint2013arXiv

User Guide for the Discrete Dipole Approximation Code DDSCAT 7.3

DDSCAT 7.3 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures in Draine & Flatau (2008), and efficient near-field calculations in Flatau & Draine (2012). DDSCAT 7.3 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for user-specified scattering directions. DDSCAT 7.3 can efficiently calculate E and B throughout a user-specified volume containing the target. This User Guide explains how to use DDSCAT 7.3 to carry out electromagnetic scattering calculations, including use of DDPOSTPROCESS, a Fortran-90 code to perform calculations with E and B at user-selected locations near the target. A number of changes have been made since the last release, DDSCAT 7.2 .

preprint2012arXiv

A Study of Heating and Cooling of the ISM in NGC 1097 with Herschel-PACS and Spitzer-IRS

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158μm/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.

preprint2012arXiv

Magnetic Nanoparticles in the Interstellar Medium: Emission Spectrum and Polarization

The presence of ferromagnetic or ferrimagnetic nanoparticles in the interstellar medium would give rise to magnetic dipole radiation at microwave and submm frequencies. Such grains may account for the strong mm-wavelength emission observed from a number of low-metallicity galaxies, including the Small Magellanic Cloud. We show how to calculate the absorption and scattering cross sections for such grains, with particular attention to metallic Fe, magnetite Fe3O4, and maghemite gamma-Fe2O3, all potentially present in the interstellar medium. The rate of Davis-Greenstein alignment by magnetic dissipation is also estimated. We determine the temperature of free-flying magnetic grains heated by starlight and we calculate the polarization of the magnetic dipole emission from both free-fliers and inclusions. For inclusions, the magnetic dipole emission is expected to be polarized orthogonally relative to the normal electric dipole radiation. Finally, we present self-consistent dielectric functions for metallic Fe, magnetite Fe3O4, and maghemite gamma-Fe2O3, enabling calculation of absorption and scattering cross sections from microwave to X-ray wavelengths.

preprint2012arXiv

Mapping the cold dust temperatures and masses of nearby Kingfish galaxies with Herschel

Taking advantage of the sensitivity and angular resolution of the Herschel Space Observatory at far-infrared and submm wavelengths, we aim to characterize the physical properties of cold dust within nearby galaxies and study the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH program, we perform 2 temperature fits of the Spitzer and Herschel photometric data (24 to 500um), with a warm and a cold component, globally and in each resolution element.At global scales, we observe ranges of values for beta_c(0.8 to 2.5) and Tc(19.1 to 25.1K).We compute maps of our parameters with beta fixed or free to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe temperature gradients as a function of radius.When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities.Gathering resolved elements in a Tc-beta_c diagram underlines an anti-correlation between the two parameters.It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar ends (coinciding with molecular gas density enhancements usually found in these locations).We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (beta_c < 2.0) lead to significantly lower dust masses compared to the beta_c < 2.0 case, with dust masses lower by up to 50% if beta_c=1.5 for instance.The working resolution affects our total dust mass estimates: masses increase from global fits to spatially-resolved fits.

preprint2012arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: NGC 628 and NGC 6946

We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.

preprint2012arXiv

Resolving the Far-IR Line Deficit: Photoelectric Heating and Far-IR Line Cooling in NGC 1097 and NGC 4559

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $μ$m, [O \one] 63 $μ$m, and [N \ii] 122 $μ$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\lesssim$ 0.95. In the warmest regions, where $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $μ$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

preprint2012arXiv

The Carriers of the Interstellar Unidentified Infrared Emission Features: Aromatic or Aliphatic?

The unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 micrometer, commonly attributed to polycyclic aromatic hydrocarbon (PAH) molecules, have been recently ascribed to coal- or kerogen-like organic nanoparticles with a mixed aromatic-aliphatic structure. However, we show in this Letter that this hypothesis is inconsistent with observations. We estimate the aliphatic fraction of the UIE carriers based on the observed intensities of the 3.4 and 6.85 micrometer emission features by attributing them exclusively to aliphatic C-H stretch and aliphatic C-H deformation vibrational modes, respectively. We derive the fraction of carbon atoms in aliphatic form to be <15%. We conclude that the UIE emitters are predominantly aromatic with aliphatic material at most a minor part of the UIE carriers. The PAH model is consistent with astronomical observations and PAHs dominate the strong UIE bands.

preprint2012arXiv

The Submm and mm Excess of the SMC: Magnetic Dipole Emission from Magnetic Nanoparticles?

The Small Magellanic Cloud (SMC) has surprisingly strong submm and mm-wavelength emission that is inconsistent with standard dust models, including those with emission from spinning dust. Here we show that the emission from the SMC may be understood if the interstellar dust mixture includes magnetic nanoparticles, emitting magnetic dipole radiation resulting from thermal fluctuations in the magnetization. The magnetic grains can be metallic iron, magnetite Fe3O4, or maghemite gamma-Fe2O3. The required mass of iron is consistent with elemental abundance constraints. The magnetic dipole emission is predicted to be polarized orthogonally to the normal electric dipole radiation if the nanoparticles are inclusions in larger grains. We speculate that other low-metallicity galaxies may also have a large fraction of the interstellar Fe in magnetic materials.

preprint2011arXiv

Calibrating Extinction-Free Star Formation Rate Diagnostics with 33GHz Free-Free Emission in NGC6946

Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of Hα+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, Hα+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

preprint2011arXiv

Common-Resolution Convolution Kernels for Space- and Ground-Based Telescopes

Multi-wavelength study of extended astronomical objects requires combining images from instruments with differing point spread functions (PSFs). We describe the construction of convolution kernels that allow one to generate (multi-wavelength) images with a common PSF, thus preserving the colors of the astronomical sources. We generate convolution kernels for the cameras of the Spitzer Space Telescope, Herschel Space Observatory, Galaxy Evolution Explorer (GALEX), Wide-field Infrared Survey Explorer (WISE), ground-based optical telescopes (Moffat functions and sum of Gaussians), and Gaussian PSFs. These kernels allow the study of the Spectral Energy Distribution (SED) of extended objects, preserving the characteristic SED in each pixel. The convolution kernels and the IDL packages used to construct and use them are made publicly available.

preprint2011arXiv

Herschel Far-Infrared and Sub-millimeter Photometry for the KINGFISH Sample of Nearby Galaxies

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

preprint2011arXiv

KINGFISH -- Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.

preprint2011arXiv

On Radiation Pressure in Static, Dusty HII Regions

Radiation pressure acting on gas and dust causes HII regions to have central densities that are lower than the density near the ionized boundary. HII regions in static equilibrium comprise a family of similarity solutions, parametrized by 3 parameters: beta, gamma, and the product (Q_0 n_rms); beta characterizes the stellar spectrum, gamma characterizes the dust/gas ratio, Q_0 is the ionizing output from the star (photons/s), and n_rms is the rms density within the ionized region. Adopting standard values for beta and gamma, varying (Q_0 n_rms) generates a one-parameter family of density profiles, ranging from nearly uniform density (small Q_0 n_rms), to hollow-sphere HII regions (large Q_0 n_rms). When (Q_0 n_rms) exceeds 10^{52} cm^{-3} s^{-1}, dusty HII regions have conspicuous central cavities, even if no stellar wind is present. For given beta, gamma and (Q_0 n_rms), a fourth quantity, which can be Q_0, determines the overall size and density of the HII region. Examples of density and emissivity profiles are given. We show how quantities of interest -- such as the peak-to-center emissivity ratio, the rms-to-mean density ratio, the edge-to-rms density ratio, and the fraction of the ionizing photons absorbed by the gas -- depend on the 3 parameters beta, gamma, and (Q_0 n_rms). For dusty HII regions, compression of the gas and dust into an ionized shell results in a substantial increase in the fraction of the >13.6 eV photons that actually ionize H (relative to a uniform density HII region with the same dust/gas ratio and density n=n_rms). We discuss the extent to which radial drift of dust grains in HII regions can alter the dust-to-gas ratio. The applicability of these solutions to real HII regions is discussed.

preprint2011arXiv

Spinning Dust Emission: Effects of irregular grain shape, transient heating and comparison with WMAP results

Planck is expected to answer crucial questions on the early Universe, but it also provides further understanding on anomalous microwave emission. Electric dipole emission from spinning dust grains continues to be the favored interpretation of anomalous microwave emission. In this paper, we present a method to calculate the rotational emission from small grains of irregular shape with moments of inertia $I_{1}> I_{2}> I_{3}$. We show that a torque-free rotating irregular grain with a given angular momentum radiates at multiple frequency modes. The resulting spinning dust spectrum has peak frequency and emissivity increasing with the degree of grain shape irregularity, which is defined by $I_{1}:I_{2}:I_{3}$. We discuss how the orientation of dipole moment $\bmu$ in body coordinates affects the spinning dust spectrum for different regimes of internal thermal fluctuations. We show that the spinning dust emissivity for the case of strong thermal fluctuations is less sensitive to the orientation of $\bmu$ than in the case of weak thermal fluctuations. We calculate spinning dust spectra for a range of gas density and dipole moment. The effect of compressible turbulence on spinning dust emission intensity is investigated. We show that the emission intensity in a turbulent medium increases by a factor from 1.2-1.4 relative to that in a uniform medium, as sonic Mach number $M_{s}$ increases from 2-7. Finally, spinning dust parameters are constrained by fitting our improved model to five-year {\it Wilkinson Microwave Anisotropy Probe} cross-correlation foreground spectra, for both the H$α$-correlated and 100 $μ$m-correlated emission spectra.

preprint2011arXiv

Spitzer IRS Detection of Molecular Hydrogen Rotational Emission Towards Translucent Clouds

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we have detected emission in the S(0), S(1), and S(2) pure-rotational (v=0-0) transitions of molecular hydrogen (H2) towards 6 positions in two translucent high Galactic latitude clouds, DCld 300.2-16.9 and LDN 1780. The detection of these lines raises important questions regarding the physical conditions inside low-extinction clouds that are far from ultraviolet radiation sources. The ratio between the S(2) flux and the flux from PAHs at 7.9 microns averages 0.007 for these 6 positions. This is a factor of about 4 higher than the same ratio measured towards the central regions of non-active Galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). Thus the environment of these translucent clouds is more efficient at producing rotationally excited H2 per PAH-exciting photon than the disks of entire galaxies. Excitation analysis finds that the S(1) and S(2) emitting regions are warm (T >300K), but comprise no more than 2% of the gas mass. We find that UV photons cannot be the sole source of excitation in these regions and suggest mechanical heating via shocks or turbulent dissipation as the dominant cause of the emission. The clouds are located on the outskirts of the Scorpius-Centaurus OB association and may be dissipating recent bursts of mechanical energy input from supernova explosions. We suggest that pockets of warm gas in diffuse or translucent clouds, integrated over the disks of galaxies, may represent a major source of all non-active galaxy H2 emission.

preprint2011arXiv

The Displaced Dusty ISM of NGC 3077: Tidal Stripping in the M 81 Triplet

We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

preprint2011arXiv

The Spitzer Spectroscopic Survey of the Small Magellanic Cloud (S4MC): Probing the Physical State of Polycyclic Aromatic Hydrocarbons in a Low-Metallicity Environment

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.

preprint2010arXiv

Enhanced dust heating in the bulges of early-type spiral galaxies

Stellar density and bar strength should affect the temperatures of the cool (T ~ 20-30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 microns. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 +/- 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23(+/-3)% hotter for morphological types earlier than Sc, and only 9(+/-3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar.

preprint2010arXiv

Far-Infrared Line Imaging of the Starburst Ring in NGC 1097 with the Herschel/PACS Spectrometer

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $μ$m, [OIII] 88 $μ$m, [NII] 122 $μ$m, [CII] 158 $μ$m and [NII] 205 $μ$m. We observe that in the [OI] 63 $μ$m, [OIII] 88 $μ$m, and [NII] 122 $μ$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $μ$m/[CII] 158 $μ$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $μ$m and [OIII] 88 $μ$m maps. Variations of the [NII] 122 $μ$m/[NII] 205 $μ$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

preprint2010arXiv

H-alpha and Free-Free Emission from the WIM

Recent observations have found the ratio of H-alpha to free-free radio continuum to be surprisingly high in the diffuse ionized ISM (the so-called WIM), corresponding to an electron temperature of only ~3000K. Such low temperatures were unexpected in gas that was presumed to be photoionized. We consider a 3-component model for the observed diffuse emission, consisting of a mix of (1) photoionized gas, (2) gas that is recombining and cooling, and (3) cool H I gas. This model can successfully reproduce the observed intensities of free-free continuum, H-alpha, and collisionally-excited lines such as NII 6583. To reproduce the low observed value of free-free to H-alpha, the PAH abundance in the photoionized regions must be lowered by a factor ~3, and ~20% of the diffuse H-alpha must be reflected from dust grains, as suggested by Wood & Reynolds (1999).

preprint2010arXiv

Improving the model of emission from spinning dust: effects of grain wobbling and transient spin-up

Observations continue to support the interpretation of the anomalous microwave foreground as electric dipole radiation from spinning dust grains as proposed by Draine and Lazarian (1998ab). In this paper we present a refinement of the original model by improving the treatment of a number of physical effects. First, we consider a disk-like grain rotating with angular velocity at an arbitrary angle with respect to the grain symmetry axis and derive the rotational damping and excitation coefficients arising from infrared emission, plasma-grain interactions and electric dipole emission. The angular velocity distribution and the electric dipole emission spectrum for grains is calculated using the Langevin equation, for cases both with and without fast internal relaxation. Our results show that, the peak emissivity of spinning dust, compared to earlier studies, increases by a factor of ~2 for the Warm Neutral Medium (WNM), the Warm Ionized Medium (WIM), the Cold Neutral Medium (CNM) and the Photodissociation Region (PDR), and by a factor ~4 for Reflection Nebulae (RN). The frequency at the emission peak also increases by factors ~1.4 to ~2 for these media. The increased emission and peak frequency result from the non-sphericity of grain shape and from the anisotropy in damping and excitation along directions parallel and perpendicular to the grain symmetry axis. Second, we provide a detailed numerical study including transient spin-up of grains by single-ion collisions. The impulses broaden the emission spectrum and increase the peak emissivity for the CNM, WNM and WIM. In addition, we present an improved treatment of rotational excitation and damping by infrared emission.

preprint2010arXiv

Mapping far-IR emission from the central kiloparsec of NGC 1097

Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.

preprint2010arXiv

Nearby Galaxies with Spitzer

We review the main advances brought by the Spitzer Space Telescope in the field of nearby galaxies studies, concentrating on a few subject areas, including: (1) the physics of the Polycyclic Aromatic Hydrocarbons that generate the mid-infrared features between ~3.5 micron and ~20 micron; (2) the use of the mid- and far-infrared emission from galaxies as star formation rate indicators; and (3) the improvement of mid-infrared diagnostics to discriminate between thermal (star-formation) and non-thermal (AGN) emission in galaxies and galaxy centers.

preprint2010arXiv

The Calibration of Monochromatic Far-Infrared Star Formation Rate Indicators

(Abridged) Spitzer data at 24, 70, and 160 micron and ground-based H-alpha images are analyzed for a sample of 189 nearby star-forming and starburst galaxies to investigate whether reliable star formation rate (SFR) indicators can be defined using the monochromatic infrared dust emission centered at 70 and 160 micron. We compare recently published recipes for SFR measures using combinations of the 24 micron and observed H-alpha luminosities with those using 24 micron luminosity alone. From these comparisons, we derive a reference SFR indicator for use in our analysis. Linear correlations between SFR and the 70 and 160 micron luminosity are found for L(70)>=1.4x10^{42} erg/s and L(160)>=2x10^{42} erg/s, corresponding to SFR>=0.1-0.3 M_sun/yr. Below those two luminosity limits, the relation between SFR and 70 micron (160 micron) luminosity is non-linear and SFR calibrations become problematic. The dispersion of the data around the mean trend increases for increasing wavelength, becoming about 25% (factor ~2) larger at 70 (160) micron than at 24 micron. The increasing dispersion is likely an effect of the increasing contribution to the infrared emission of dust heated by stellar populations not associated with the current star formation. The non-linear relation between SFR and the 70 and 160 micron emission at faint galaxy luminosities suggests that the increasing transparency of the interstellar medium, decreasing effective dust temperature, and decreasing filling factor of star forming regions across the galaxy become important factors for decreasing luminosity. The SFR calibrations are provided for galaxies with oxygen abundance 12+Log(O/H)>8.1. At lower metallicity the infrared luminosity no longer reliably traces the SFR because galaxies are less dusty and more transparent.

preprint2010arXiv

User Guide for the Discrete Dipole Approximation Code DDSCAT 7.1

DDSCAT 7.1 is an open-source Fortran-90 software package applying the discrete dipole approximation to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The targets may be isolated entities (e.g., dust particles), but may also be 1-d or 2-d periodic arrays of "target unit cells", allowing calculation of absorption, scattering, and electric fields around arrays of nanostructures. The theory of the DDA and its implementation in DDSCAT is presented in Draine (1988) and Draine & Flatau (1994), and its extension to periodic structures (and near-field calculations) in Draine & Flatau (2008). DDSCAT 7.1 includes support for MPI, OpenMP, and the Intel Math Kernel Library (MKL). DDSCAT supports calculations for a variety of target geometries. Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to "import" arbitrary target geometries into the code. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix. This User Guide explains how to use DDSCAT 7.1 to carry out electromagnetic scattering calculations. DDfield, a Fortran-90 code to calculate E and B at user-selected locations near the target, is included in the distribution. A number of changes have been made since the last release, DDSCAT 7.0 .

preprint2009arXiv

Modeling Porous Dust Grains with Ballistic Aggregates. II. Light Scattering Properties

We study the light scattering properties of random ballistic aggregates constructed in Shen et al. (Paper I). Using the discrete-dipole-approximation, we compute the scattering phase function and linear polarization for random aggregates with various sizes and porosities, and with two different compositions: 100% silicate and 50% silicate-50% graphite. We investigate the dependence of light scattering properties on wavelength, cluster size and porosity using these aggregate models. We find that while the shape of the phase function depends mainly on the size parameter of the aggregates, the linear polarization depends on both the size parameter and the porosity of the aggregates, with increasing degree of polarization as the porosity increases. Contrary to previous studies, we argue that monomer size has negligible effects on the light scattering properties of ballistic aggregates, as long as the constituent monomer is smaller than the incident wavelength up to 2*pi*a_0/lambda\sim 1.6 where a_0 is the monomer radius. Previous claims for such monomer size effects are in fact the combined effects of size parameter and porosity. Finally, we present aggregate models that can reproduce the phase function and polarization of scattered light from the AU Mic debris disk and from cometary dust, including the negative polarization observed for comets at scattering angles 160<theta<180 deg. These aggregates have moderate porosities, P\sim 0.6, and are of sub-micron-size for the debris disk case, or micron-size for the comet case.

preprint2009arXiv

Radial distribution of stars, gas and dust in SINGS galaxies. II. Derived dust properties

We present a detailed analysis of the radial distribution of dust properties in the SINGS sample, performed on a set of UV, IR and HI surface brightness profiles, combined with published molecular gas profiles and metallicity gradients. The internal extinction, derived from the TIR-to-FUV luminosity ratio, decreases with radius, and is larger in Sb-Sbc galaxies. The TIR-to-FUV ratio correlates with the UV spectral slope beta, following a sequence shifted to redder UV colors with respect to that of starbursts. The star formation history (SFH) is identified as the main driver of this departure. We have also derived radial profiles of the total dust mass surface density, the fraction of the dust mass contributed by PAHs, the fraction of the dust mass heated by very intense starlight and the intensity of the radiation field heating the grains. The dust profiles are exponential, their radial scale-length being constant from Sb to Sd galaxies (only ~10% larger than the stellar scale-length). Many S0/a-Sab galaxies have central depressions in their dust radial distributions. The PAH abundance increases with metallicity for 12+\log(O/H)<9, and at larger metallicities the trend flattens and even reverses, with the SFH being a plausible underlying driver for this behavior. The dust-to-gas ratio is also well correlated with metallicity and therefore decreases with galactocentric radius.

preprint2009arXiv

The Discovery of Vibrationally-Excited H_2 in the Molecular Cloud near GRB 080607

GRB 080607 has provided the first strong observational signatures of molecular absorption bands toward any galaxy hosting a gamma-ray burst. Despite the identification of dozens of features as belonging to various atomic and molecular (H_2 and CO) carriers, many more absorption features remained unidentified. Here we report on a search among these features for absorption from vibrationally-excited H_2, a species that was predicted to be produced by the UV flash of a GRB impinging on a molecular cloud. Following a detailed comparison between our spectroscopy and static, as well as dynamic, models of H_2* absorption, we conclude that a column density of 10^{17.5+-0.2} cm^{-2} of H_2* was produced along the line of sight toward GRB 080607. Depending on the assumed amount of dust extinction between the molecular cloud and the GRB, the model distance between the two is found to be in the range 230--940 pc. Such a range is consistent with a conservative lower limit of 100 pc estimated from the presence of Mg I in the same data. These distances show that substantial molecular material is found within hundreds of pc from GRB 080607, part of the distribution of clouds within the GRB host galaxy.

preprint2008arXiv

CMBPol Mission Concept Study: Foreground Science Knowledge and Prospects

We report on our knowledge of Galactic foregrounds, as well as on how a CMB satellite mission aiming at detecting a primordial B-mode signal (CMBPol) will contribute to improving it. We review the observational and analysis techniques used to constrain the structure of the Galactic magnetic field, whose presence is responsible for the polarization of Galactic emissions. Although our current understanding of the magnetized interstellar medium is somewhat limited, dramatic improvements in our knowledge of its properties are expected by the time CMBPol flies. Thanks to high resolution and high sensitivity instruments observing the whole sky at frequencies between 30 GHz and 850 GHz, CMBPol will not only improve this picture by observing the synchrotron emission from our galaxy, but also help constrain dust models. Polarized emission from interstellar dust indeed dominates over any other signal in CMBPol's highest frequency channels. Observations at these wavelengths, combined with ground-based studies of starlight polarization, will therefore enable us to improve our understanding of dust properties and of the mechanism(s) responsible for the alignment of dust grains with the Galactic magnetic field. CMBPol will also shed new light on observations that are presently not well understood. Morphological studies of anomalous dust and synchrotron emissions will indeed constrain their natures and properties, while searching for fluctuations in the emission from heliospheric dust will test our understanding of the circumheliospheric interstellar medium. Finally, acquiring more information on the properties of extra-Galactic sources will be necessary in order to maximize the cosmological constraints extracted from CMBPol's observations of CMB lensing. (abridged)

preprint2008arXiv

The relations among 8, 24, and 160 micron dust emission within nearby spiral galaxies

We investigate the relations among the stellar continuum-subtracted 8 micron polycyclic aromatic hydrocarbon (PAH 8 micron) emission, 24 micron hot dust emission, and 160 micron cold dust emission in fifteen nearby face-on spiral galaxies in the Spitzer Infrared Nearby Galaxies Survey sample. The relation between PAH 8 and 24 micron emission measured in ~2 kpc regions is found to exhibit a significant amount of scatter, and strong spatial variations are observed in the (PAH 8 micron)/24 micron surface brightness ratio. In particular, the (PAH 8 micron)/24 micron surface brightness ratio is observed to be high in the diffuse interstellar medium and low in bright star-forming regions and other locations with high 24 micron surface brightness. PAH 8 micron emission is found to be well-correlated with 160 micron emission on spatial scales of ~2 kpc, and the (PAH 8 micron)/160 micron surface brightness ratio is generally observed to increase as the 160 micron surface brightness increases. These results suggest that the PAHs are associated with the diffuse, cold dust that produces most of the 160 micron emission in these galaxies, and the variations in the (PAH 8 micron)/160 micron ratio may generally be indicative of either the intensity or the spectrum of the interstellar radiation field that is heating both the PAHs and the diffuse interstellar dust.

preprint2007arXiv

The Calibration of Mid-Infrared Star Formation Rate Indicators

With the goal of investigating the degree to which the mid-infrared emission traces the star formation rate (SFR), we analyze Spitzer 8 um and 24 um data of star-forming regions in a sample of 33 nearby galaxies with available HST/NICMOS images in the Paschen-alpha (1.8756 um) emission line. The galaxies are drawn from the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample, and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and Luminous Infrared Galaxies are also included in the analysis. Both the stellar-continuum-subtracted 8 um emission and the 24 um emission correlate with the extinction-corrected Pa-alpha line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed non-linear trend of the 24 um emission versus number of ionizing photons, including the modest deficiency of 24 um emission in the low metallicity regions, which results from a combination of decreasing dust opacity and dust temperature at low luminosities. Conversely, the trend of the 8 um emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 um emission is contributed, in larger measure than the 24 um emission, by dust heated by non-ionizing stellar populations, in agreement with previous findings. Two SFR calibrations, one using the 24 um emission and the other using a combination of the 24 um and H-alpha luminosities (Kennicutt et al. 2007), are presented. No calibration is presented for the 8 um emission, because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by on-going star formation.

preprint2007arXiv

Warm molecular hydrogen in the Spitzer SINGS galaxy sample

(simplified) Results on the properties of warm H2 in 57 normal galaxies are derived from H2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of H2, the most abundant constituent of the molecular ISM, to more common systems (L_FIR = e7 to 6e10 L_sun) of all morphological and nuclear types. The S(1) transition is securely detected in the nuclear regions of 86% of SINGS galaxies with stellar masses above 10^9.5 M_sun. The derived column densities of warm H2 (T > ~100 K), even though averaged over kiloparsec-scale areas, are commensurate with those of resolved PDRs; the median of the sample is 3e20 cm-2. They amount to between 1% and >30% of the total H2. The power emitted in the sum of the S(0) to S(2) transitions is on average 30% of the [SiII] line power, and ~4e-4 of the total infrared power (TIR) within the same area for star-forming galaxies, which is consistent with excitation in PDRs. The fact that H2 emission scales tightly with PAH emission, even though the average radiation field intensity varies by a factor ten, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. A large fraction of the 25 LINER/Sy targets, however, strongly depart from the rest of the sample, in having warmer H2 in the excited states, and an excess of H2 emission with respect to PAHs, TIR and [SiII]. We propose a threshold in H2 to PAH power ratios, allowing the identification of low-luminosity AGNs by an excess H2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect, in nearly half the star-forming targets, non-equilibrium ortho to para ratios, consistent with FUV pumping combined with incomplete ortho-para thermalization by collisions, or possibly non-equilibrium PDR fronts advancing into cold gas.

preprint2006arXiv

Spitzer and JCMT Observations of the Active Galactic Nucleus in the Sombrero Galaxy (NGC 4594)

We present Spitzer 3.6-160 micron images, Spitzer mid-infrared spectra, and JCMT SCUBA 850 micron images of the Sombrero Galaxy (NGC 4594), an Sa galaxy with a 10^9 M_solar low luminosity active galactic nucleus (AGN). The brightest infrared sources in the galaxy are the nucleus and the dust ring. The spectral energy distribution of the AGN demonstrates that, while the environment around the AGN is a prominent source of mid-infrared emission, it is a relatively weak source of far-infrared emission, as had been inferred for AGN in previous research. The weak nuclear 160 micron emission and the negligible polycyclic aromatic hydrocarbon emission from the nucleus also implies that the nucleus is a site of only weak star formation activity and the nucleus contains relatively little cool interstellar gas needed to fuel such activity. We propose that this galaxy may be representative of a subset of low ionization nuclear emission region galaxies that are in a quiescent AGN phase because of the lack of gas needed to fuel circumnuclear star formation and Seyfert-like AGN activity. Surprisingly, the AGN is the predominant source of 850 micron emission. We examine the possible emission mechanisms that could give rise to the 850 micron emission and find that neither thermal dust emission, CO line emission, bremsstrahlung emission, nor the synchrotron emission observed at radio wavelengths can adequately explain the measured 850 micron flux density by themselves. The remaining possibilities for the source of the 850 micron emission include a combination of known emission mechanisms, synchrotron emission that is self-absorbed at wavelengths longer than 850 microns, or unidentified spectral lines in the 850 micron band.

preprint2006arXiv

The Effect of Star Formation on the Far-Infrared--Radio Correlation within Galaxies

Using data obtained for twelve galaxies as part of the {\it Spitzer} Infrared Nearby Galaxies Survey (SINGS) and the Westerbork Synthesis Radio Telescope (WSRT)-SINGS radio continuum survey, we study how star formation activity affects the far-infrared (FIR)--radio correlation {\it within} galaxies by testing a phenomenological model, which describes the radio image as a smeared version of the FIR image. The physical basis of this description is that cosmic-ray (CR) electrons will diffuse measurably farther than the mean free path of dust-heating photons before decaying by synchrotron radiation. This description works well in general. Galaxies with higher infrared surface brightnesses have best-fit smoothing scale-lengths of a few hundred parsecs, substantially shorter than those for lower surface brightness galaxies. We interpret this result to suggest that galaxies with higher disk averaged star formation rates have had a recent episode of enhanced star formation and are characterized by a higher fraction of young CR electrons that have traveled only a few hundred parsecs from their acceleration sites in supernova remnants compared to galaxies with lower star formation activity.

preprint2006arXiv

The Mid-Infrared Spectrum of Star-Forming Galaxies: Global Properties of PAH Emission

We present a sample of low-resolution 5-38um Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described, and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths which indicate most normal galaxies undergo A_V < ~3 magnitudes averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20%, and is suppressed at metallicities Z < ~Z_sun/4, as well as in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5x) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8um features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGN. The AGN may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGN may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled, and points to more than a factor of two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z=0-2.5.

preprint2004arXiv

IRAC Observations of M81

IRAC images of M81 show three distinct morphological constituents: a smooth distribution of evolved stars with bulge, disk, and spiral arm components; a clumpy distribution of dust emission tracing the spiral arms; and a pointlike nuclear source. The bulge stellar colors are consistent with M-type giants, and the disk colors are consistent with a slightly younger population. The dust emission generally follows the blue and ultraviolet emission, but there are large areas that have dust emission without ultraviolet and smaller areas with ultraviolet but little dust emission. The former are presumably caused by extinction, and the latter may be due to cavities in the gas and dust created by supernova explosions. The nucleus appears fainter at 8 um than expected from ground-based 10 um observations made four years ago.

preprint2004arXiv

Mid-Infrared IRS Spectroscopy of NGC 7331: A First Look at the SINGS Legacy

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

preprint2004arXiv

Spatially Resolved Ultraviolet, H-alpha, Infrared, and Radio Star Formation in M81

We present Multiband Imaging Photometer for Spitzer (MIPS) observations of M81 at 24, 70, and 160 micron. The grand design nature of M81 is clearly seen, showing two well resolved spiral arms containing numerous bright star forming regions. The MIPS images reveal a significant amount of cold dust associated with the spiral arms. We investigate the variation of the ultraviolet (UV), H-alpha, and infrared (IR) luminosities and star formation rate (SFR) indicators across the face of M81 using the MIPS images and archival UV and H-alpha images. For regions in M81, we find that UV and H-alpha SFRs (uncorrected for dust attenuation) are always lower than the IR SFR. The cause of this behavior is dust attenuation and/or using SFR calibrations appropriate for entire galaxies, not regions in galaxies. The characteristics of the dust attenuation for the regions indicate the dust grains and/or geometry are different from those in starburst galaxies. The behavior of the infrared-radio correlation in M81 is seen to vary from the global average, with variations correlated with the morphology of M81.

preprint2003arXiv

Discrete-dipole approximation with polarizabilities that account for both finite wavelength and target geometry

The discrete-dipole approximation (DDA) is a powerful method for calculating absorption and scattering by targets that have sizes smaller than or comparable to the wavelength of the incident radiation. We present a new prescription -- the Surface-Corrected Lattice Dispersion Relation (SCLDR) -- for assigning the dipole polarizabilities that takes into account both target geometry and finite wavelength. We test the SCLDR in DDA calculations using spherical and ellipsoidal targets and show that for a fixed number of dipoles, the SCLDR prescription results in increased accuracy in the calculated cross sections for absorption and scattering. We discuss extension of the SCLDR prescription to irregular targets.

preprint2003arXiv

Scattering by Interstellar Dust Grains: Optical and Ultraviolet

Scattering and absorption properties at optical and ultraviolet wavelengths are calculated for an interstellar dust model consisting of carbonaceous grains and amorphous silicate grains. Polarization as a function of scattering angle is calculated for selected wavelengths from the IR to the vacuum UV. The widely-used Henyey-Greenstein phase function provides a good approximation for the scattering phase function between ~0.4 and 1 micron, but fails to fit the calculated phase functions at shorter wavelengths and longer wavelengths. A new analytic phase function is presented. It is exact at long wavelengths, and provides a good fit to the numerically-calculated phase function for lambda > 0.27um. Observational determinations of the scattering albedo and <cos(theta)> show considerable disagreement, especially in the ultraviolet. Possible reasons for this are discussed.

preprint2003arXiv

SINGS: The SIRTF Nearby Galaxies Survey

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies, and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multi-wavelength library of ancillary and complementary observations, including radio continuum, HI, CO, submillimeter, BVRIJHK, H-alpha, Paschen-alpha, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.

preprint2002arXiv

Infrared Emission from Interstellar Dust. III. The Small Magellanic Cloud

The infrared (IR) emission from interstellar dust in the Small Magellanic Cloud (SMC) is modelled using a mixture of amorphous silicate and carbonaceous grains, including a population of polycyclic aromatic hydrocarbon (PAH) molecules. (1) It is shown that this dust model is able to reproduce the spectral energy distribution from near-IR to far-IR for the entire SMC Bar region, provided the PAH abundance in the SMC Bar region is very low. (2) The IR spectrum of the SMCB1#1 molecular cloud can also be reproduced by our dust model provided the PAH abundance is increased relative to the overall SMC Bar. The PAHs in SMCB1#1 incorporate ~3% of the SMC C abundance, compared to <0.4% in the SMC Bar. (3) The spectrum of SMCB1#1 is best reproduced if the PAH mixture has intrinsic IR band strengths which differ from the band strengths which best fit Milky Way PAH mixtures. The variation in the PAH IR band strengths would imply different PAH mixtures, presumably a consequence of differing metallicity or environmental conditions. Other possibilities such as super-hydrogenation of PAHs and softening of the starlight spectrum are also discussed.

preprint2001arXiv

Dust Grain Size Distributions and Extinction in the Milky Way, LMC, and SMC

We construct size distributions for carbonaceous and silicate grain populations in different regions of the Milky Way, LMC, and SMC. The size distributions include sufficient very small carbonaceous grains (including polycyclic aromatic hydrocarbon molecules) to account for the observed infrared and microwave emission from the diffuse interstellar medium. Our distributions reproduce the observed extinction of starlight, which varies depending upon the interstellar environment through which the light travels. As shown by Cardelli, Clayton & Mathis in 1989, these variations can be roughly parameterized by the ratio of visual extinction to reddening, R_V. We adopt a fairly simple functional form for the size distribution, characterized by several parameters. We tabulate these parameters for various combinations of values for R_V and b_C, the C abundance in very small grains. We also find size distributions for the line of sight to HD 210121, and for sightlines in the LMC and SMC. For several size distributions, we evaluate the albedo and scattering asymmetry parameter, and present model extinction curves extending beyond the Lyman limit.

preprint2000arXiv

Photoelectric Emission from Interstellar Dust: Grain Charging and Gas Heating

We model the photoelectric emission from and charging of interstellar dust and obtain photoelectric gas heating efficiencies as a function of grain size and the relevant ambient conditions. Using realistic grain size distributions, we evaluate the net gas heating rate for various interstellar environments, and find less heating for dense regions characterized by R_V=5.5 than for diffuse regions with R_V=3.1. We provide fitting functions which reproduce our numerical results for photoelectric heating and recombination cooling for a wide range of interstellar conditions. In a separate paper we will examine the implications of these results for the thermal structure of the interstellar medium. Finally, we investigate the potential importance of photoelectric heating in H II regions, including the warm ionized medium. We find that photoelectric heating could be comparable to or exceed heating due to photoionization of H for high ratios of the radiation intensity to the gas density. We also find that photoelectric heating by dust can account for the observed variation of temperature with distance from the galactic midplane in the warm ionized medium.

preprint1999arXiv

H_2 Absorption and Fluorescence for Gamma Ray Bursts in Molecular Clouds

If a gamma ray burst with strong UV emission occurs in a molecular cloud, there will be observable consequences resulting from excitation of the surrounding H2. The UV pulse from the GRB will pump H2 into vibrationally-excited levels which produce strong absorption at wavelengths < 1650 A. As a result, both the prompt flash and later afterglow will exhibit strong absorption shortward of 1650 A, with specific spectroscopic features. Such a cutoff in the emission from GRB 980329 may already have been observed by Fruchter et al.; if so, GRB 980329 was at redshift 3.0 < z < 4.4 . BVRI photometry of GRB 990510 could also be explained by H2 absorption if GRB 990510 is at redshift 1.6 < z < 2.3. The fluorescence accompanying the UV pumping of the H2 will result in UV emission from the GRB which can extend over days or months, depending on parameters of the ambient medium and beaming of the GRB flash. The 7.5-13.6 eV fluorescent luminosity is \sim 10^{41.7} erg/s for standard estimates of the parameters of the GRB and the ambient medium. Spectroscopy can distinguish this fluorescent emission from other possible sources of transient optical emission, such as a supernova.

preprint1999arXiv

Nuclear Spin Relaxation within Interstellar Grains

We identify a new mechanism of internal dissipation of rotational kinetic energy in spinning dust grains, arising from the reorientation of nuclear angular momentum, e.g., spins of protons. Grain rotation induces magnetization of the nuclear spin system, with net alignment of nuclear spins parallel to the grain angular velocity. When the grain does not rotate around a principal axis, the nuclear magnetization vector precesses in grain body coordinates, resulting in dissipation of energy. The analogous process involving electron spins was discovered by Purcell and termed ``Barnett relaxation''. We revisit the physics of the ``Barnett relaxation'' process and correct the estimate for the Barnett relaxation rate. We show that nuclear relaxation can be orders of magnitude more important than Barnett relaxation. This finding implies that the processes of ``thermal flipping'' and ``thermal trapping'' are important for a broad range of grain sizes.

preprint1998arXiv

Electric Dipole Radiation from Spinning Dust Grains

We discuss the rotational excitation of small interstellar grains and the resulting electric dipole radiation from spinning dust. Attention is given to excitation and damping of rotation by: collisions with neutrals; collisions with ions; plasma drag; emission of infrared radiation; emission of microwave radiation; photoelectric emission; and formation of H_2 on the grain surface. We introduce dimensionless functions F and G which allow direct comparison of the contributions of different mechanisms to rotational drag and excitation. Emissivities are estimated for dust in different phases of the interstellar medium, including diffuse HI, warm HI, low-density photoionized gas, and cold molecular gas. Spinning dust grains can explain much, and perhaps all, of the 14-50 GHz background component recently observed in CBR studies. It should be possible to detect rotational emission from small grains by ground-based observations of molecular clouds.

preprint1995arXiv

Grain Destruction in Interstellar Shocks

Interstellar shock waves can erode and destroy grains present in the shocked gas, primarily as the result of sputtering and grain-grain collisions. Uncertainties in current estimates of sputtering yields are reviewed. Results are presented for the simple case of sputtering of fast grains being stopped in cold gas. An upper limit is derived for sputtering of refractory grains in C-type MHD shocks: shock speeds $v_s \gtrsim 50 \kms$ are required for return of more than 30\% of the silicate to the gas phase. Sputtering can also be important for removing molecular ice mantles from grains in two-fluid MHD shock waves in molecular gas. Recent estimates of refractory grain lifetimes against destruction in shock waves are summarized, and the implications of these short lifetimes are discussed.