Diffuse ISM Research Articles

Exploring the physical properties of local star-forming ULIRGs from the ultraviolet to the infrared

We present an application of the da Cunha, Charlot & Elbaz (2008) model of the spectral energy distribution (SEDs) of galaxies from the ultraviolet to far-infrared to a small pilot sample of purely star-forming Ultra-Luminous Infrared Galaxies (ULIRGs). We interpret the observed SEDs of 16 ULIRGs using this physically-motivated model which accounts for the emission of stellar populations from the ultraviolet to the near-infrared and for the attenuation by dust in two components: an optically-thick starburst component and the diffuse ISM. The infrared emission is computed by assuming that all the energy absorbed by dust in these components is re-radiated at mid- and far-infrared wavelengths. This model allows us to derive statistically physical properties including star formation rates, stellar masses, as well as temperatures and masses of different dust components and plausible star formation histories. We find that, although the ultraviolet to near-infrared emission represents only a small fraction of the total power radiated by ULIRGs, observations in this wavelength range are important to understand the properties of the stellar populations and dust attenuation in the diffuse ISM of these galaxies. Furthermore, our analysis indicates that the use of mid-infrared spectroscopy from the Infrared Spectrograph on the Spitzer Space Telescope is crucial to obtain realistic estimates of the extinction to the central energy source, mainly via the depth of the 9.7-micron silicate feature, and thus accurately constrain the total energy balance. Our findings are consistent with the notion that, in the local Universe, the physical properties of ULIRGs are fundamentally different from those of galaxies with lower infrared luminosities and that local ULIRGs are the result of merger-induced starbursts. [abridged]

The extinction law for molecular clouds

We determine the extinction curve from the UV to the near-IR for molecular clouds and investigate whether current models can adequately explain this wavelength dependence of the extinction. The aim is also to interpret the extinction in terms of H2 column density. We applied five different methods, including a new method for simultaneously determining the reddening law and the classification of the background stars. Our method is based on multicolour observations and a grid of model atmospheres. We confirm that the extinction law can be adequately described by a single parameter, RV (the selective to absolute extinction), in accordance with earlier findings. The RV value for B 335 is RV = 4.8. The reddening curve can be accurately reproduced by model calculations. By assuming that all the silicon is bound in silicate grains, we can interpret the reddening in terms of column density, NH = 4.4 (\pm0.5) \times 1021 EI-Ks cm-2, corresponding to NH = 2.3 (\pm0.2) \times 1021 \cdot AV cm-2, close to that of the diffuse ISM, (1.8-2.2) \times 1021 cm-2 . We show that the density of the B 335 globule outer shells can be modelled as an evolved Ebert-Bonnor gas sphere with {\rho} \propto r-2, and estimate the mass of this globule to 2.5 Msun

REGULATION OF STAR FORMATION RATES IN MULTIPHASE GALACTIC DISKS: A THERMAL/DYNAMICAL EQUILIBRIUM MODEL

We develop a model for regulation of galactic star formation rates Sigma_SFR in disk galaxies, in which ISM heating by stellar UV plays a key role. By requiring simultaneous thermal and (vertical) dynamical equilibrium in the diffuse gas, and star formation at a rate proportional to the mass of the self-gravitating component, we obtain a prediction for Sigma_SFR as a function of the total gaseous surface density Sigma and the density of stars + dark matter, rho_sd. The physical basis of this relationship is that thermal pressure in the diffuse ISM, which is proportional to the UV heating rate and therefore to Sigma_SFR, must adjust to match the midplane pressure set by the vertical gravitational field. Our model applies to regions where Sigma < 100 Msun/pc^2. In low-Sigma_SFR (outer-galaxy) regions where diffuse gas dominates, the theory predicts Sigma_SFR \propto Sigma (rho_sd)^1/2. The decrease of thermal equilibrium pressure when Sigma_SFR is low implies, consistent with observations, that star formation can extend (with declining efficiency) to large radii in galaxies, rather than having a sharp cutoff. The main parameters entering our model are the ratio of thermal pressure to total pressure in the diffuse ISM, the fraction of diffuse gas that is in the warm phase, and the star formation timescale in self-gravitating clouds; all of these are (in principle) direct observables. At low surface density, our model depends on the ratio of the mean midplane FUV intensity (or thermal pressure in the diffuse gas) to the star formation rate, which we set based on Solar neighborhood values. We compare our results to recent observations, showing good agreement overall for azimuthally-averaged data in a set of spiral galaxies. For the large flocculent spiral galaxies NGC 7331 and NGC 5055, the correspondence between theory and observation is remarkably close.

Photoluminescence of hydrogenated amorphous carbons

Context. Hydrogenated amorphous carbons (a-C:H or HAC) have proved to be excellent analogs of interstellar dust observed in galaxies diffuse interstellar medium (DISM) through infrared vibrational absorption bands (3.4 μ m, 6.8 μ m, and 7.2 μ m bands). They exhibit photoluminescence (PL) after excitation by UV-visible photons, and are possible carriers for the extended red emission (ERE), a broad red emission band observed in various interstellar environments.Aims. As many candidate materials/molecules can photoluminesce in the visible, along with the carrier abundance, the PL efficiency represents one of the strongest constraints set by such ERE observations. We wish to precisely characterize the PL behavior of a-C:H as a family of materials.Methods. The a-C:H samples are produced in the form of films deposited on substrates by plasma-enhanced chemical vapor deposition. The produced films were analyzed in transmission by UV-visible and IR spectroscopy, and the wavelength dependent PL spectra were recorded. The intrinsic absolute quantum yield η was then rigorously calculated taking self-absorption of the PL by the film and interfaces effects into account.Results. A wide range of different laboratory synthesized a-C:H were analyzed. Their PL properties are dependent on the optical gap : when decreases from 4.3 eV to 2.8 eV, the a-C:H vary from highly () yellow photoluminescent soft materials to hard materials that emit a wider PL band in the red spectral range, with a lower efficiency (η ~ 0.01–0.1%). For any given a-C:H, the PL characteristics (central wavelength, band width and efficiency) are found to be essentially constant over the explored excitation range ( 250 nm). We compared the characteristics of the produced interstellar dust analog to the constraints imposed by the ERE observations.Conclusions. As for ERE observations, PL efficiencies and band widths of a-C:H are both correlated to the PL central wavelengths. The excitation responsible for the a-C:H emission is efficient over a wide spectral range that matches the ERE excitation. The present a-C:H encounter difficulties for the diffuse ISM ERE observations (%) in simultaneously satisfying the high quantum yield criteria and PL spectral characteristics. We still need to investigate the role of a small number of residual oxygen atoms in the laboratory-produced a-C:H network in quenching the PL yield, as well as to consider the interstellar temperature effect for our analogs.

The CO luminosity and CO-H2conversion factor of diffuse ISM: does CO emission trace dense molecular gas?

Aims: We wish to separate and quantify the CO luminosity and CO-H2 conversion factor applicable to diffuse but partially-molecular ISM when H2 and CO are present but C+ is the dominant form of gas-phase carbon. Methods: We discuss galactic lines of sight observed in \HI, HCO+ and CO where CO emission is present but the intervening clouds are diffuse (locally \AV\ $\la 1$ mag) with relatively small CO column densities $\NCO \la 2\times10^{16}\pcc$. We separate the atomic and molecular fractions statistically using \EBV\ as a gauge of the total gas column density and compare NH2 to the observed CO brightness. Results: Although there are H2-bearing regions where CO emission is too faint to be detected, the mean ratio of integrated CO brightness to NH2 for diffuse ISM does not differ from the usual value of 1\K km/s of integrated CO brightness per $2\times10^{20}$ H2 $\pcc$ . Moreover, the luminosity of diffuse CO viewed perpendicular to the galactic plane is 2/3 that seen at the Solar galactic radius in surveys of CO emission near the galactic plane. Conclusions: Commonality of the CO-H2 conversion factors in diffuse and dark clouds can be understood from considerations of radiative transfer and CO chemistry. There is unavoidable confusion between CO emission from diffuse and dark gas and misattribution of CO emission from diffuse to dark or giant molecular clouds. The character of the ISM is different from what has been believed if CO and H2 that have been attributed to molecular clouds on the verge of star formation are actually in more tenuous, gravitationally-unbound diffuse gas.

THESPITZERSURVEY OF THE SMALL MAGELLANIC CLOUD (S3MC): INSIGHTS INTO THE LIFE CYCLE OF POLYCYCLIC AROMATIC HYDROCARBONS

We present the results of modeling dust SEDs across the SMC with the aim of mapping the distribution of PAHs in a low-metallicity environment. Using Spitzer Survey of the SMC (S3MC) photometry from 3.6-160 um over the main star-forming regions of the Wing and Bar along with spectral mapping from 5-38 um from the Spitzer Spectroscopic Survey of the SMC (S4MC) in selected regions, we model the dust SED and emission spectrum to determine the fraction of dust in PAHs across the SMC. We use the regions of overlapping photometry and spectroscopy to test the reliability of the PAH fraction as determined from SED fits alone. The PAH fraction in the SMC is low compared to the Milky Way and variable--with relatively high fractions (q_PAH~1-2%) in molecular clouds and low fractions in the diffuse ISM (<q_PAH>=0.6%). We use the map of PAH fraction across the SMC to test a number of ideas regarding the production, destruction and processing of PAHs in the ISM. We find weak or no correlation between the PAH fraction and the distribution of carbon AGB stars, the location of supergiant H I shells and young SN remnants, or the turbulent Mach number. We find that the PAH fraction is correlated with CO intensity, peaks in the dust surface density and the molecular gas surface density as determined from 160 um emission. The PAH fraction is high in regions of active star-formation, as predicted by its correlation with molecular gas, but is suppressed in H II regions. Because the PAH fraction in the diffuse ISM is generally very low--in accordance with previous work on modeling the SED of the SMC--and the PAH fraction is relatively high in molecular regions, we suggest that PAHs are destroyed in the diffuse ISM of the SMC and/or PAHs are forming in molecular clouds. We discuss the implications of these observations for our understanding of the PAH life cycle, particularly in low-metallicity galaxies.

Star Formation and Gas Dynamics in Galactic Disks: Physical Processes and Numerical Models

AbstractStar formation depends on the available gaseous “fuel” as well as galactic environment, with higher specific star formation rates where gas is predominantly molecular and where stellar (and dark matter) densities are higher. The partition of gas into different thermal components must itself depend on the star formation rate, since a steady state distribution requires a balance between heating (largely from stellar UV for the atomic component) and cooling. In this presentation, I discuss a simple thermal and dynamical equilibrium model for the star formation rate in disk galaxies, where the basic inputs are the total surface density of gas and the volume density of stars and dark matter, averaged over ~kpc scales. Galactic environment is important because the vertical gravity of the stars and dark matter compress gas toward the midplane, helping to establish the pressure, and hence the cooling rate. In equilibrium, the star formation rate must evolve until the gas heating rate is high enough to balance this cooling rate and maintain the pressure imposed by the local gravitational field. In addition to discussing the formulation of this equilibrium model, I review the current status of numerical simulations of multiphase disks, focusing on measurements of quantities that characterize the mean properties of the diffuse ISM. Based on simulations, turbulence levels in the diffuse ISM appear relatively insensitive to local disk conditions and energetic driving rates, consistent with observations. It remains to be determined, both from observations and simulations, how mass exchange processes control the ratio of cold-to-warm gas in the atomic ISM.

A study of major mergers using a multi‐phase ISM code

AbstractGalaxy interactions are a common phenomenon in clusters of galaxies. Especially major mergers are of particular importance, because they can change the morphological type of galaxies. They have an impact on the mass function of galaxies and they trigger star formation – the main driver of the Galactic Matter Cycle. Therefore, we conducted a study of major mergers by means of a multi‐phase ISM code. This code is based on a TREE‐SPH‐code combined with a sticky particle method allowing for star formation controlled by the properties of a multi‐phase ISM. This is in contrast to the usually implemented Schmidt law depending mainly on the gas density. Previously, this code was used on isolated galaxies. Since our star formation recipe is not restricted to a special type of galaxy, it is interesting to apply it to interacting galaxies, too. Our study on major mergers includes a research of global properties of the interacting system, namely the star formation rate and the star formation efficiency, the evaporation and condensation rates, as well as the mass exchange of distinct components, namely stars, diffuse ISM, and clouds. Investigating these properties provides insight to interrelations between various physical processes. The results indicate that the star formation efficiency as well as the evaporation and condensation rates are influenced by the interaction (© 2009 WILEY‐VCH Verlag GmbH &amp; Co. KGaA, Weinheim)

Boron abundances in diffuse interstellar clouds

AbstractWe present a comprehensive survey of B abundances in diffuse interstellar clouds from HST/STIS observations along 56 Galactic sight lines. Our sample is the result of a complete search of archival STIS data for the B II λ1362 resonance line, with each detection confirmed by the presence of absorption from other dominant ions at the same velocity. The data probe a range of astrophysical environments including both high-density regions of massive star formation as well as low-density paths through the Galactic halo, allowing us to clearly define the trend of B depletion onto interstellar grains as a function of gas density. Many extended sight lines exhibit complex absorption profiles that trace both local gas and gas associated with either the Sagittarius-Carina or Perseus spiral arm. Our analysis indicates a higher B/O ratio in the inner Sagittarius-Carina spiral arm than in the vicinity of the Sun, which may suggest that B production in the current epoch is dominated by a secondary process. The average gas-phase B abundance in the warm diffuse ISM [log ϵ(B) = 2.38±0.10] is consistent with the abundances determined for a variety of Galactic disk stars, but is depleted by 60% relative to the solar system value. Our survey also reveals sight lines with enhanced B abundances that potentially trace recent production of 11B either by cosmic-ray or neutrino-induced spallation. Such sight lines will be key to discerning the relative importance of the two production routes for 11B synthesis.

Quasar feedback: more bang for your buck

We propose a two-stage model for the effects of feedback from a bright quasar on the cold gas in a galaxy. It is difficult for feedback from near the accretion disk to directly impact dense molecular clouds at ~kpc. But if such feedback can drive a weak wind or outflow in the hot, diffuse ISM (a relatively 'easy' task), then in the wake of such an outflow passing over a cold cloud, a combination of instabilities will drive the cloud material to effectively expand in the direction perpendicular to the outflow. Such expansion dramatically increases the effective cross section of the cloud material and makes it more susceptible to ionization and momentum coupling from absorption of the incident quasar radiation field. Even a moderate effect of this nature can dramatically alter the ability of clouds at large radii to be fully ionized and driven into a secondary outflow by radiation pressure. Since the amount of momentum and volume which can be ionized by observed quasar radiation field is more than sufficient to affect the entire cold gas supply once it has been altered in this manner (and the 'initial' feedback need only initiate a moderate wind in the low-density hot gas), this reduces by an order of magnitude the required energy budget for feedback to affect a host galaxy. Instead of ~5% of the radiated energy (~100% momentum) needed if the initial feedback must directly heat or blow out the galactic gas, if only ~0.5% of the luminosity (~10% momentum) can couple to drive the initial hot outflow, this mechanism could be efficient. This amounts to hot gas outflow rates from near the accretion disk of only 5-10% of the BH accretion rate.

MOLECULAR SURVIVAL IN EVOLVED PLANETARY NEBULAE: DETECTION OF H 2 CO, c-C 3 H 2 , AND C 2 H IN THE HELIX

H2CO, c-C3H2, and C2H have been identified in the neutral envelope of the highly evolved planetary nebula (PN), the Helix (also know as NGC 7293). Emission from these species were detected toward a peak position in CO, 372'' east of the central star, using the facilities of the Arizona Radio Observatory (ARO). C2H and c-C3H2 were identified on the basis of their 3 mm transitions, measured with the ARO 12 m, while five lines of H2CO were observed using the 12 m at 2 and 3 mm and the ARO Submillimeter Telescope at 1 mm. From a radiative transfer analysis of the formaldehyde emission, the molecular material was determined to have a density of n(H2) ~3 × 105 cm–3, with a kinetic temperature of T kin ~20 K. Column densities for C2H, H2CO, and c-C3H2 of N tot ~1.4 × 1013 cm–2, 1.1 × 1012 cm–2, and 3 × 1011 cm–2, respectively, were derived, corresponding to fractional abundances relative to H2 of f (H2CO) = 1 × 10–7, f (c-C3H2) = 3 × 10–8, and f (C2H) = 1 × 10–6 . The physical conditions found support the notion that molecules in evolved PNe survive in dense clumps in pressure equilibrium, shielded from photodissociation. The presence of H2CO, c-C3H2, and C2H, along with the previously observed species CN, HNC, HCN, and HCO+, indicates that a relatively complex chemistry can occur in the late stages of PN evolution, despite potentially destructive ultraviolet radiation. These molecules have also been observed in diffuse clouds, suggesting a possible connection between molecular material in evolved PNe and the diffuse ISM.

EVIDENCE FOR DUST EVOLUTION WITHIN THE TAURUS COMPLEX FROMSPITZERIMAGES

We present Spitzer images of the Taurus Complex (TC) and take advantage of the sensitivity and spatial resolution of the observations to characterize the diffuse IR emission across the cloud. This work highlights evidence of dust evolution within the translucent sections of the archetype reference for studies of quiescent molecular clouds. We combine Spitzer 160 um and IRAS 100 um observations to produce a dust temperature map and a far-IR dust opacity map at 5' resolution. The average dust temperature is about 14.5K with a dispersion of +/-1K across the cloud. The far-IR dust opacity is a factor 2 larger than the average value for the diffuse ISM. This opacity increase and the attenuation of the radiation field (RF) both contribute to account for the lower emission temperature of the large grains. The structure of the TC significantly changes in the mid-IR images that trace emission from PAHs and VSGs. We focus our analysis of the mid-IR emission to a range of ecliptic latitudes where the zodiacal light residuals are small. Within this cloud area, there are no 8 and 24 um counterparts to the brightest 160 um emission features. Conversely, the 8 and 24 um images reveal filamentary structure that is strikingly inconspicuous in the 160 um and extinction maps. The IR colors vary over sub-parsec distances across this filamentary structure. We compare the observed colors with model calculations quantifying the impact of the RF intensity and the abundance of stochastically heated particles on the dust SED. To match the range of observed colors, we have to invoke variations by a factor of a few of both the interstellar RF and the abundance of PAHs and VSGs. We conclude that within this filamentary structure a significant fraction of the dust mass cycles in and out the small size end of the dust size distribution.

Physical and Chemical Processes in the Diffuse and Dense Neutral ISM

Most of the mass of the interstellar medium is in neutral phases that are heated by far-ultraviolet radiation from star forming regions. The gas, in turn, emits line radiation in the infrared and submillimeter that potentially has great diagnostic value. In this paper, we will review the physical processes (including heating, cooling, gas phase and grain surface chemistry) that give rise to this line emission and present diagnostic diagrams. We will also discuss several examples where models of neutral gas emission have been used to interpret line emission, including observations of H 2 O and O 2 in molecular clouds and an example which points out the caution that must be used when applying the diagnostics to kiloparsec-sized regions in other galaxies.

Magnetic fields and the location of the PDR

I review recent studies of the emission-line regions in Orion and M17. Both have similar geometries, a bubble of hot shocked gas surrounding the central star cluster, with H^+, H^0, and H_2 regions, often referred to as H II regions, PDRs, and molecular clouds, forming successive shells on the surface of a molecular cloud. The magnetic fields in the H^0 regions have been measured with 21 cm Zeeman polarization and are found to be 1 -- 2 dex stronger than the field in the diffuse ISM. The regions appear to be in rough hydrostatic equilibrium. The H^+ region is pushed away from the star cluster by starlight radiation pressure. Since most starlight is in ionizing radiation, most of its outward push will act on the H^+ region and then on to the H^0 region. The magnetic pressure in the H^0 region balances the momentum in starlight and together they set the location of the H^0 region. The picture is that, when the star cluster formed, it created a bubble of ionized gas which expanded and compressing surrounding H^0 and H_2 regions. The magnetic field was amplified until its pressure was able to support the momentum in starlight. This offers a great simplification in understanding the underlying physics that establishes parameters for PDR models.

Dust in the SINGS Galaxy Sample

Dust models have been applied to interpret the IR emission from 65 galaxies in the SINGS galaxy survey, allowing dust masses, PAH fractions, and starlight intensities to be determined. The dust models reproduce the global SEDs without requiring “cold” (T 12 K) dust. The bulk of the IR emission appears to be radiated by dust in the diffuse ISM, with starlight intensities similar to the intensity in the solar neighborhood. A fraction fPDR of the IR emission is radiated by dust in regions of high starlight intensities (e.g., PDRs), with median fPDR = 0.08. The dust composition is characterized by qPAH, the fraction of the dust mass contributed by PAHs with <103 C atoms. Galaxies with metallicities Z/Z๏ < 0.2 tend to have qPAH 0.01, but galaxies with Z/Z๏ > 0.4 have qPAH = 0.02 - 0.05. Over the entire sample, with Z/Z๏ ranging from 0.1 to 1, the observations are consistent with a constant fraction of interstellar metals being in dust, with Mdust/MH ≈ 0.01(Z/Z๏).

A centimetre-wave excess over free-free emission in planetary nebulae

We report a centimetre-wave (cm-wave, 5-31GHz) excess over free-free emission in PNe. Accurate 31 and 250GHz measurements show that the 31GHz flux densities in our sample are systematically higher than the level of optically thin free-free continuum extrapolated from 250GHz. The 31GHz excess is observed, within one standard deviation, in all 18 PNe with reliable 31 and 250GHz data, and is significant in 9 PNe. The only exception is the peculiar object M2-9, whose radio spectrum is that of an optically thick stellar wind. On average the fraction of non-free-free emission represents 51% of the total flux density at 31GHz, with a scatter of 11%. The average 31-250GHz spectral index of our sample is = -0.43+-0.03 (in flux density, with a scatter of 0.14). The 31--250 GHz drop is reminiscent of the anomalous foreground observed in the diffuse ISM by CMB anisotropy experiments. The 5--31 GHz spectral indices are consistent with both flat spectra and spinning dust emissivities, given the 10% calibration uncertainty of the comparison 5GHz data. But a detailed study of the objects with the largest cm-excess, including the low frequency data available in the literature, shows that present spinning dust models cannot alone explain the cm-excess in PNe. Although we have no definitive interpretation of our data, the least implausible explanation involves a synchrotron component absorbed by a cold nebular screen. We give flux densities for 37 objects at 31GHz, and for 26 objects at 250GHz.

The Relationship between the Optical Depth of the 9.7 μm Silicate Absorption Feature and Infrared Differential Extinction in Dense Clouds

We have examined the relationship between the optical depth of the 9.7 μm silicate absorption feature (τ_(9.7)) and the near-infrared color excess, E(J - K_s), in the Serpens, Taurus, IC 5146, Chameleon I, Barnard 59, and Barnard 68 dense clouds/cores. Our data set, based largely on Spitzer IRS spectra, spans E(J - K_s) = 0.3-10 mag (corresponding to visual extinction between about 2 and 60 mag). All lines of sight show the 9.7 μm silicate feature. Unlike in the diffuse ISM where a tight linear correlation between the 9.7 μm silicate feature optical depth and the extinction (A_V) is observed, we find that the silicate feature in dense clouds does not show a monotonic increase with extinction. Thus, in dense clouds, τ_(9.7) is not a good measure of total dust column density. With few exceptions, the measured τ_(9.7) values fall well below the diffuse ISM correlation line for E(J - K_s) > 2 mag (A_V > 12 mag). Grain growth via coagulation is a likely cause of this effect.

Dust Masses, PAH Abundances, and Starlight Intensities in the SINGS Galaxy Sample

Physical dust models are presented for 65 galaxies in SINGS that are strongly detected in the four IRAC bands and three MIPS bands. For each galaxy we estimate (1) the total dust mass, (2) the fraction of the dust mass contributed by PAHs, and (3) the intensity of the starlight heating the dust grains. We find that spiral galaxies have dust properties resembling the dust in the local region of the Milky Way, with similar dust-to-gas ratio and similar PAH abundance. The observed SEDs, including galaxies with SCUBA photometry, can be reproduced by dust models that do not require cold (T ≾ 10 K) dust. The dust-to-gas ratio is observed to be dependent on metallicity. In the interstellar media of galaxies with A_O ≡ 12 + log_(10)(O/H) > 8.1, grains contain a substantial fraction of interstellar Mg, Si, and Fe. Galaxies with A_O 8.1 have a median q_(PAH) = 3.55%. The derived dust masses favor a value X_(CO) ≈ 4 × 10^(20) cm^(-2) (K km s^(-1))^(-1) for the CO-to-H_2 conversion factor. Except for some starbursting systems (Mrk 33, Tol 89, NGC 3049), dust in the diffuse ISM dominates the IR power.

The Infrared Extinction Law at Extreme Depth in a Dark Cloud Core

We combined sensitive near-infrared data obtained with ground-based imagers on the ESO NTT and VLT telescopes with space mid-infrared data acquired with the IRAC imager on the Spitzer Space Telescope to calculate the extinction law A_\lambda/A_K as a function of \lambda between 1.25 and 7.76 micron to an unprecedented depth in Barnard 59, a star forming, dense core located in the Pipe Nebula. The ratios A_\lambda/A_K were calculated from the slopes of the distributions of sources in color-color diagrams \lambda-K vs. H-K. The distributions in the color-color diagrams are fit well with single slopes to extinction levels of A_K ~ 7 (A_V ~ 59 mag). Consequently, there appears to be no significant variation of the extinction law with depth through the B59 line of sight. However, when slopes are translated into the relative extinction coefficients A_\lambda/A_K, we find an extinction law which departs from the simple extrapolation of the near-infrared power law extinction curve, and agrees more closely with a dust extinction model for a cloud with a total to selective absorption R_V=5.5 and a grain size distribution favoring larger grains than those in the diffuse ISM. Thus, the difference we observe could be possibly due to the effect of grain growth in denser regions. Finally, the slopes in our diagrams are somewhat less steep than those from the study of Indebetouw et al. (2005) for clouds with lower column densities, and this indicates that the extinction law between 3 and 8 micron might vary slightly as a function of environment.

Infrared Extinction toward Nearby Star‐forming Regions

We present an independent estimate of the interstellar extinction law for the Spitzer IRAC bands as well as a first attempt at extending the law to the 24micron MIPS band. The source data for these measurements are observations of five nearby star-forming regions: the Orion A cloud, NGC 2068/71, NGC 2024/23, Serpens and Ophiuchus. Color excess ratios E(H-Ks)/E(Ks-[lambda]) were measured for stars without infrared excess dust emission from circumstellar disks/envelopes. For four of these five regions, the extinction laws are similar at all wavelengths and differ systematically from a previous determination of the extinction law, which was dominated by the diffuse ISM, derived for the IRAC bands. This difference could be due to the difference in the dust properties of the dense molecular clouds observed here and those of the diffuse ISM. The extinction law at longer wavelengths toward the Ophiuchus region lies between that to the other four regions studied here and that for the ISM. In addition, we extended our extinction law determination to 24micron for Serpens and NGC 2068/71 using Spitzer MIPS data. We compare these results against several ISO extinction law determinations, although in each case there are assumptions which make absolute comparison uncertain. However, our work confirms a relatively flatter extinction curve from 4 - 8micron than the previously assumed standard, as noted by all of these recent studies. The extinction law at 24micron is consistent with previous measurements and models, although there are relatively large uncertainties.