Showing 1–50 of 204 results for author: Draine, B

The Metallicity Dependence of PAH Emission in Galaxies I: Insights from Deep Radial Spitzer Spectroscopy

Authors: Cory M. Whitcomb, J. -D. T. Smith, Karin Sandstrom, Carl A. Starkey, Grant P. Donnelly, Bruce T. Draine, Evan D. Skillman, Daniel A. Dale, Lee Armus, Brandon S. Hensley, Thomas S. -Y. Lai, Robert C. Kennicutt

Abstract: We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features f… ▽ More We use deep Spitzer mid-infrared spectroscopic maps of radial strips across three nearby galaxies with well-studied metallicity gradients (M101, NGC 628, and NGC 2403) to explore the physical origins of the observed deficit of polycyclic aromatic hydrocarbons (PAHs) at sub-solar metallicity (i.e. the PAH-metallicity relation or PZR). These maps allow us to trace the evolution of all PAH features from 5-18 $μ$m as metallicity decreases continuously from solar ($Z_\odot$) to 0.2 $Z_\odot$. The total PAH to dust luminosity ratio remains relatively constant until reaching a threshold of $\sim$$\frac{2}{3}$$Z_\odot$, below which it declines smoothly but rapidly. The PZR has been attributed to preferential destruction of the smallest grains in the hard radiation environments found at low metallicity. In this scenario, a decrease in emission from the shortest wavelength PAH features is expected. In contrast, we find a steep decline in long wavelength power below $Z_\odot$, especially in the 17 $μ$m feature, with the shorter wavelength PAH bands carrying an increasingly large fraction of power at low metallicity. We use newly developed grain models to reproduce the observed PZR trends, including these variations in fractional PAH feature strengths. The model that best reproduces the data employs an evolving grain size distribution that shifts to smaller sizes as metallicity declines. We interpret this as a result of inhibited grain growth at low metallicity, suggesting continuous replenishment in the interstellar medium is the dominant process shaping the PAH grain population in galaxies. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 22 pages, 13 figures, submitted for publication in ApJ

Sensitivity of Polarization to Grain Shape: II. Aggregates

Authors: B. T. Draine

Abstract: A previous study (Paper I) investigated the polarization properties of a variety of simple convex grain shapes, some of which were found to be consistent with the observed polarization properties of interstellar dust from far-ultraviolet to far-infrared. Here we study the optical properties of 45 non-convex shapes, all aggregates of $N$ equal-sized spheres. We consider $N=2$, $N=3$, and $N=256$ ra… ▽ More A previous study (Paper I) investigated the polarization properties of a variety of simple convex grain shapes, some of which were found to be consistent with the observed polarization properties of interstellar dust from far-ultraviolet to far-infrared. Here we study the optical properties of 45 non-convex shapes, all aggregates of $N$ equal-sized spheres. We consider $N=2$, $N=3$, and $N=256$ random aggregates obtained from 3 different aggregation schemes. We also consider "trimmed" $N=256$ aggregates obtained by systematically trimming initially random aggregates to increase either flattening or elongation. The "macroporosities" of the studied aggregates range from ${\cal P}_{\rm macro}=0.18$ (for the $N=2$ bisphere) to ${\cal P}_{\rm macro}\approx 0.85$ (for the $N=256$ "BA" aggregates). The only aggregates consistent with observations of starlight polarization and polarized thermal emission are shapes that have been trimmed to increase their asymmetry. If interstellar grains are high-porosity aggregates, there must be processes causing extreme elongation or flattening; if not, interstellar grains must be dominated by fairly compact structures, with at most moderate porosities. The ratio of polarization in the 10$μ$m silicate feature to starlight polarization in the optical is shown to be insensitive to porosity and shape. X-ray scattering may be the best tool to determine the porosity of interstellar grains. We propose that modest porosities of interstellar grains could be the result of "photolytic densification". High polarization fractions observed in some Class-0 cores require processes to reduce porosities and/or increase asymmetries of aggregates in dense regions. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: ApJ, accepted. Comments welcomed

The Impact of an AGN on PAH Emission in Galaxies: the Case of Ring Galaxy NGC 4138

Authors: G. P. Donnelly, J. D. T. Smith, B. T. Draine, A. Togi, T. S. -Y. Lai, L. Armus, D. A. Dale, V. Charmandaris

Abstract: We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138,… ▽ More We present a focused study of radially-resolved varying PAH emission in the low-luminosity AGN-host NGC 4138 using deep Spitzer/IRS spectral maps. Using new model PAH spectra, we investigate whether these variations could be associated with changes to the PAH grain size distribution due to photodestruction by the AGN. Separately, we model the effects of the varying radiation field within NGC 4138, and we use this model to predict the corresponding changes in the PAH emission spectrum. We find that PAH band ratios are strongly variable with radius in this galaxy with short-to-long wavelength band ratios peaking in the starburst ring. The changing mix of starlight appears to have a considerable effect on the trends in these band ratios, and our radiation model predicts the shapes of these trends. However, the amplitude of observed variation is ~2.5 times larger than predicted for some ratios. A cutoff of small grains in the PAH size distribution, as has been suggested for AGN, together with changes in PAH ionization fraction could explain the behavior of the shorter bands, but this model fails to reproduce longer band behaviors. Additionally, we find that short-to-long wavelength PAH band ratios increase slightly within ~270pc of the center, suggesting that the AGN may directly influence PAH emission there. △ Less

Submitted 23 April, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 19 pages, 13 figures, published in ApJ. Updated values in radiation model to the same as the published version

Sensitivity of Polarization to Grain Shape: I. Convex Shapes

Authors: B. T. Draine

Abstract: Aligned interstellar grains produce polarized extinction (observed at wavelengths from the far-ultraviolet to the mid-infrared), and polarized thermal emission (observed at far-infrared and submm wavelengths). The grains must be quite nonspherical, but the actual shapes are unknown. The \emph{relative} efficacy for aligned grains to produce polarization at optical vs.\ infrared wavelengths depends… ▽ More Aligned interstellar grains produce polarized extinction (observed at wavelengths from the far-ultraviolet to the mid-infrared), and polarized thermal emission (observed at far-infrared and submm wavelengths). The grains must be quite nonspherical, but the actual shapes are unknown. The \emph{relative} efficacy for aligned grains to produce polarization at optical vs.\ infrared wavelengths depends on particle shape. The discrete dipole approximation is used to calculate polarization cross sections for 20 different convex shapes, for wavelengths from $0.1μ$m to $100μ$m, and grain sizes $a_{\rm eff}$ from $0.05μ$m to $0.3μ$m. Spheroids, cylinders, square prisms, and triaxial ellipsoids are considered. Minimum aspect ratios required by the observed starlight polarization are determined. Some shapes can also be ruled out because they provide too little or too much polarization at far-infrared and sub-mm wavelengths. The ratio of $10μ$m polarization to integrated optical polarization is almost independent of grain shape, varying by only $\pm8\%$ among the viable convex shapes; thus, at least for convex grains, uncertainties in grain shape cannot account for the discrepancy between predicted and observed 10$μ$m polarization toward Cyg OB2-12. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 33 pages, 20 figures. Astrophysical Journal, accepted

Gas, dust, and the CO-to-molecular gas conversion factor in low-metallicity starbursts

Authors: L. K. Hunt, F. Belfiore, F. Lelli, B. T. Draine, A. Marasco, S. Garcia-Burillo, G. Venturi, F. Combes, A. Weiß, C. Henkel, K. M. Menten, F. Annibali, V. Casasola, M. Cignoni, A. McLeod, M. Tosi, M. Beltran, A. Concas, G. Cresci, M. Ginolfi, N. Kumari, F. Mannucci

Abstract: The factor relating CO emission to molecular hydrogen column density, XCO, is still subject to uncertainty, in particular at low metallicity. Here, to quantify XCO at two different spatial resolutions, we exploit a dust-based method together with ALMA 12-m and ACA data and HI maps of three nearby metal-poor starbursts, NGC625, NGC1705, and NGC5253. Dust opacity at 250pc resolution is derived based… ▽ More The factor relating CO emission to molecular hydrogen column density, XCO, is still subject to uncertainty, in particular at low metallicity. Here, to quantify XCO at two different spatial resolutions, we exploit a dust-based method together with ALMA 12-m and ACA data and HI maps of three nearby metal-poor starbursts, NGC625, NGC1705, and NGC5253. Dust opacity at 250pc resolution is derived based on dust temperatures estimated by fitting two-temperature modified blackbodies to Herschel PACS data. By using the HI maps, we are then able to estimate dust-to-gas ratios in the atomic-gas dominated regions, and infer total gas column densities and H2 column densities as the difference with HI. Finally, from the ACA CO(1-0) maps, we derive XCO. We use a similar technique with 40 pc ALMA 12-m data for the three galaxies, but instead derive dust attenuation at 40 pc resolution from reddening maps based on VLT/MUSE data. At 250 pc resolution, XCO $\sim$ 10^22 - 10^23 cm^-2 / K.km/s, 5-1000 times the Milky Way value, with much larger values than would be expected from a simple metallicity dependence. Instead at 40 pc resolution, XCO again shows large variation, but is roughly consistent with a power-law metallicity dependence, given the Z $\sim$ 1/3 Zsun metal abundances of our targets. The large scatter in both estimations could imply additional parameter dependence, that we have investigated by comparing XCO with the observed velocity-integrated brightness temperatures, ICO, as predicted by recent simulations. Indeed, larger XCO is significantly correlated with smaller ICO, but with slightly different slopes and normalizations than predicted by theory. Such behavior can be attributed to the increasing fraction of CO-faint H2 gas with lower spatial resolution. This confirms the idea the XCO is multi-variate, depending not only on metallicity but also on CO brightness temperature and beam size. △ Less

Submitted 11 May, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: A&A, in press (16 figures, 1 appendix)

Journal ref: A&A 675, A64 (2023)

Mid-Infrared Polarization of the Diffuse Interstellar Medium toward CygOB2-12

Authors: Charles M. Telesco, Frank Varosi, Christopher Wright, Bruce T. Draine, Sergio Jose Fernández Acosta, Christopher Packham

Abstract: We present the first mid-IR detection of the linear polarization toward the star CygOB2-12, a luminous blue hypergiant that, with Av of ~10 mag of foreground extinction, is a benchmark in the study of the properties of dust in the diffuse interstellar medium. The 8-13 micrometer spectropolarimetry, obtained with the CanariCam multi-mode camera at the Gran Telescopio Canarias (GTC), shows clear tre… ▽ More We present the first mid-IR detection of the linear polarization toward the star CygOB2-12, a luminous blue hypergiant that, with Av of ~10 mag of foreground extinction, is a benchmark in the study of the properties of dust in the diffuse interstellar medium. The 8-13 micrometer spectropolarimetry, obtained with the CanariCam multi-mode camera at the Gran Telescopio Canarias (GTC), shows clear trends with wavelength characteristic of silicate grains aligned in the interstellar magnetic field. The maximum polarization, detected with 7.8 statistical significance near 10.2 micrometers, is (1.24 +/- 0.28) % with position angle (126 +/- 8) deg. We comment on these measurements in the context of recent models for the dust composition in the diffuse interstellar medium. △ Less

Submitted 26 October, 2022; v1 submitted 23 October, 2022; originally announced October 2022.

Comments: Accepted Astrophysical Journal Letters, October 9, 2022

The Astrodust+PAH Model: A Unified Description of the Extinction, Emission, and Polarization from Dust in the Diffuse Interstellar Medium

Authors: Brandon S. Hensley, B. T. Draine

Abstract: We present a new model of interstellar dust in which large grains are a single composite material, ``astrodust,'' and nanoparticle-sized grains come in distinct varieties including polycyclic aromatic hydrocarbons (PAHs). We argue that a single-composition model for grains larger than $\sim$0.02 $μ$m most naturally explains the lack of frequency dependence in the far-infrared (FIR) polarization fr… ▽ More We present a new model of interstellar dust in which large grains are a single composite material, ``astrodust,'' and nanoparticle-sized grains come in distinct varieties including polycyclic aromatic hydrocarbons (PAHs). We argue that a single-composition model for grains larger than $\sim$0.02 $μ$m most naturally explains the lack of frequency dependence in the far-infrared (FIR) polarization fraction and the characteristic ratio of optical to FIR polarization. We derive a size distribution and alignment function for 1.4:1 oblate astrodust grains that, with PAHs, reproduce the mean wavelength dependence and polarization of Galactic extinction and emission from the diffuse interstellar medium while respecting constraints on solid phase abundances. All model data and Python-based interfaces are made publicly available. △ Less

Submitted 6 June, 2023; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 34 pages, 19 figures. Data available at https://dataverse.harvard.edu/dataverse/astrodust

Journal ref: ApJ, 948, 55 (2023)

Dust dynamics in RAMSES -- I. Methods and turbulent acceleration

Authors: Eric R. Moseley, Romain Teyssier, B. T. Draine

Abstract: Supernova ejecta and stellar winds are believed to produce interstellar dust grains with relatively large sizes. Smaller grains can be produced via the shattering of large grains that have been stochastically accelerated. To understand this stochastic acceleration, we have implemented novel magnetohydrodynamic(MHD)-particle-in-cell(PIC) methods into the astrophysical fluid code RAMSES. We treat du… ▽ More Supernova ejecta and stellar winds are believed to produce interstellar dust grains with relatively large sizes. Smaller grains can be produced via the shattering of large grains that have been stochastically accelerated. To understand this stochastic acceleration, we have implemented novel magnetohydrodynamic(MHD)-particle-in-cell(PIC) methods into the astrophysical fluid code RAMSES. We treat dust grains as a set of massive ``superparticles'' that experience aerodynamic drag and Lorentz force. We subject our code to a range of numerical tests designed to validate our method in different physical conditions, as well as to illustrate possible mechanisms by which grains can be accelerated. As a final test as well as a foundation for future work, we present the results of decaying dusty MHD turbulence simulations with grain parameters chosen to resemble 1-2 $μ$m grains in typical cold neutral medium conditions. We find that in these conditions, these grains can be effectively accelerated to well beyond their shattering velocities. This is true for both electrically charged and neutral grains. While the peak of the gas-grain relative drift velocity distribution is higher for neutral grains, the drift velocity distribution of charged grains exhibits an extended exponential tail out to much greater velocities. Even so, the shapes of the distributions are such that the extra gas-grain coupling provided by the Lorentz force offers grains relative protection from shattering. We also discuss the connection between our simulations and the relatively pristine ~$μ$m sized presolar grains that do not appear to have undergone significant wear in their lifetimes. △ Less

Submitted 28 October, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: 23 pages, 13 figures, Accepted for publication in MNRAS

An Optical Spectrum of the Diffuse Galactic Light from BOSS and IRIS

Authors: Blake Chellew, Timothy D. Brandt, Brandon S. Hensley, Bruce T. Draine, Eve Matthaey

Abstract: We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 A, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-two increase in signal to noise. We combine these… ▽ More We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 A, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-two increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemisphere. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 A, as evidence for luminescence in the form of extended red emission (ERE). The observed strength of the 4000 A break indicates that at most ~7% of the dust-correlated light at 4000 A can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence independent of measurements of extinction. △ Less

Submitted 4 July, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 16 pages, 8 figures. Revised version after reviewer comments

Journal ref: The Astrophysical Journal, Volume 932, Number 2. June 22, 2022

On Far-Infrared and Submm Circular Polarization

Authors: B. T. Draine

Abstract: Interstellar dust grains are often aligned. If the grain alignment direction varies along the line of sight, the thermal emission becomes circularly-polarized. In the diffuse interstellar medium, the circular polarization at far-infrared and submm wavelengths is predicted to be very small, and probably unmeasurable. However, circular polarization may reach detectable levels in infrared dark clouds… ▽ More Interstellar dust grains are often aligned. If the grain alignment direction varies along the line of sight, the thermal emission becomes circularly-polarized. In the diffuse interstellar medium, the circular polarization at far-infrared and submm wavelengths is predicted to be very small, and probably unmeasurable. However, circular polarization may reach detectable levels in infrared dark clouds and protoplanetary disks. Measurement of circular polarization could help constrain the structure of the magnetic field in infrared dark clouds, and may shed light on the mechanisms responsible for grain alignment in protoplanetary disks. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Comments: ApJ, accepted

Benchmarking Dust Emission Models in M101

Authors: Jeremy Chastenet, Karin Sandstrom, I-Da Chiang, Brandon S. Hensley, Bruce T. Draine, Karl D. Gordon, Eric W. Koch, Adam K. Leroy, Dyas Utomo, Thomas G. Williams

Abstract: We present a comparative study of four physical dust models and two single-temperature modified blackbody models by fitting them to the resolved WISE, Spitzer, and Herschel photometry of M101 (NGC 5457). Using identical data and a grid-based fitting technique, we compare the resulting dust and radiation field properties derived from the models. We find that the dust mass yielded by the different m… ▽ More We present a comparative study of four physical dust models and two single-temperature modified blackbody models by fitting them to the resolved WISE, Spitzer, and Herschel photometry of M101 (NGC 5457). Using identical data and a grid-based fitting technique, we compare the resulting dust and radiation field properties derived from the models. We find that the dust mass yielded by the different models can vary by up to factor of 3 (factor of 1.4 between physical models only), although the fits have similar quality. Despite differences in their definition of the carriers of the mid-IR aromatic features, all physical models show the same spatial variations for the abundance of that grain population. Using the well determined metallicity gradient in M101 and resolved gas maps, we calculate an approximate upper limit on the dust mass as a function of radius. All physical dust models are found to exceed this maximum estimate over some range of galactocentric radii. We show that renormalizing the models to match the same Milky Way high latitude cirrus spectrum and abundance constraints can reduce the dust mass differences between models and bring the total dust mass below the maximum estimate at all radii. △ Less

Submitted 23 February, 2021; v1 submitted 22 January, 2021; originally announced January 2021.

Comments: Accepted for publication in ApJ -- 28 pp (12 figures) + appendices

Using the Starlight Polarization Efficiency Integral to Constrain Shapes and Porosities of Interstellar Grains

Authors: B. T. Draine, Brandon S. Hensley

Abstract: We present a new method for using the observed starlight polarization and polarized submm emission to constrain the shapes and porosities of interstellar grains. We present the modified picket fence approximation (MPFA), and verify that it is sufficiently accurate for modeling starlight polarization. We introduce the starlight polarization integral $Π_{\rm obs}$ as a measure of overall strength of… ▽ More We present a new method for using the observed starlight polarization and polarized submm emission to constrain the shapes and porosities of interstellar grains. We present the modified picket fence approximation (MPFA), and verify that it is sufficiently accurate for modeling starlight polarization. We introduce the starlight polarization integral $Π_{\rm obs}$ as a measure of overall strength of the observed polarization of starlight, and the starlight polarization efficiency integral $Φ$ to characterize the effectiveness of different grain types for producing polarization of starlight. The starlight polarization integral $Π_{\rm obs}$ determines the mass-weighted alignment $\langle f_{\rm align}\rangle$ of the grains. Approximating the aligned grains in the interstellar medium as spheroids, we use $Π_{\rm obs}/Φ$ to show that the observed starlight polarization constrains the grains to have a minimum degree of asphericity. For porosity ${\cal P}=0$, the minimum axial ratio is $\sim$1.4 for oblate spheroids, or $\sim$1.8 for prolate spheroids. If the grains are porous, more extreme axial ratios are required. The same grains that produce the starlight polarization are able to provide the observed polarized emission at submm wavelengths, but with further limits on shape and porosity. Porosities ${\cal P}>0.75$ are ruled out. If interstellar grains can be approximated by astrodust spheroids, we predict the ratio of 10$μ{\rm m}$ polarization to starlight polarization $p_V$: $p(10μ{\rm m})/p_V=0.222\pm0.026$. For Cyg OB2-12 we predict $p(10μ{\rm m})=(2.1\pm0.3)\%$, which should be observable. △ Less

Submitted 18 January, 2021; originally announced January 2021.

Comments: submitted to ApJ

Excitation of Polycyclic Aromatic Hydrocarbon (PAH) Emission: Dependence on Size Distribution, Ionization, and Starlight Spectrum and Intensity

Authors: B. T. Draine, Aigen Li, Brandon S. Hensley, L. K. Hunt, K. Sandstrom, J. -D. T. Smith

Abstract: Using physical models, we study the sensitivity of polycyclic aromatic hydrocarbon (PAH) emission spectra to the character of the illuminating starlight, to the PAH size distribution, and to the PAH charge distribution. Starlight models considered range from the emission from a 3 Myr-old starburst, rich in far-ultraviolet (FUV) radiation, to the FUV-poor spectrum of the very old population of the… ▽ More Using physical models, we study the sensitivity of polycyclic aromatic hydrocarbon (PAH) emission spectra to the character of the illuminating starlight, to the PAH size distribution, and to the PAH charge distribution. Starlight models considered range from the emission from a 3 Myr-old starburst, rich in far-ultraviolet (FUV) radiation, to the FUV-poor spectrum of the very old population of the M31 bulge. A wide range of starlight intensities is considered. The effects of reddening in dusty clouds are investigated for different starlight spectra. For a fixed PAH abundance parameter $q_{\rm PAH}$ (the fraction of the total dust mass in PAHs with $<10^3$ C atoms), the fraction of the IR power appearing in the PAH emission features can vary by a factor of two as the starlight spectrum varies from FUV-poor (M31 bulge) to FUV-rich (young starburst). We show how $q_{\rm PAH}$ can be measured from the strength of the 7.7$μ$m emission. The fractional power in the 17$μ$m feature can be suppressed by high starlight intensities. △ Less

Submitted 11 August, 2021; v1 submitted 13 November, 2020; originally announced November 2020.

Comments: published in The Astrophysical Journal

Journal ref: Astrophysical Journal 917:3 (2021)

Where's the Dust?: The Deepening Anomaly of Microwave Emission in NGC 4725 B

Authors: E. J. Murphy, B. S. Hensley, S. T. Linden, B. T. Draine, D. Dong, E. Momjian, G. Helou, A. S. Evans

Abstract: We present new Atacama Large Millimeter Array (ALMA) observations towards NGC 4725 B, a discrete, compact, optically-faint region within the star-forming disk of the nearby galaxy NGC 4725 that exhibits strong anomalous microwave emission (AME). These new ALMA data include continuum observations centered at 92, 133, 203, and 221 GHz accompanied by spectral observations of the $^{12}$CO (… ▽ More We present new Atacama Large Millimeter Array (ALMA) observations towards NGC 4725 B, a discrete, compact, optically-faint region within the star-forming disk of the nearby galaxy NGC 4725 that exhibits strong anomalous microwave emission (AME). These new ALMA data include continuum observations centered at 92, 133, 203, and 221 GHz accompanied by spectral observations of the $^{12}$CO ($J=2\rightarrow1$) line. NGC 4725 B is detected in the continuum at all frequencies, although the detection at 203 GHz is marginal. While molecular gas is not detected at the exact location of NGC 4725 B, there is molecular gas in the immediate vicinity (i.e., $\lesssim 100$ pc) along with associated diffuse 8 $μ$m emission. When combined with existing Very Large Array continuum data at 1.5, 3, 5.5, 9, 14, 22, 33, and 44 GHz, the spectrum is best fit by a combination of AME, synchrotron, and free-free emission that is free-free absorbed below $\sim6$ GHz. Given the strength of the AME, there is surprisingly no indication of millimeter dust emission associated with NGC 4725 B on $\lesssim$6\arcsec~spatial scales at the sensitivity of the ALMA interferometric data. Based on the properties of the nearest molecular gas complex and the inferred star formation rate, NGC 4725 B is consistent with being an extremely young ($\sim 3-5$ Myr) massive ($\lesssim 10^{5} M_{\odot}$) cluster that is undergoing active cluster feedback. However, the lack of millimeter thermal dust emission is difficult to reconcile with a spinning dust origin of the 30 GHz emission. On the other hand, modeling NGC 4725 B as a new class of background radio galaxy is also unsatisfactory. △ Less

Submitted 5 November, 2020; originally announced November 2020.

Comments: 15 pages, 4 figures, accepted for publication in ApJ Letters

The Dielectric Function of "Astrodust" and Predictions for Polarization in the 3.4$μ$m and 10$μ$m Features

Authors: B. T. Draine, Brandon S. Hensley

Abstract: The dielectric function of interstellar dust material is modeled using observations of extinction and polarization in the infrared, together with estimates for the mass of interstellar dust. The "astrodust" material is assumed to be a mix of amorphous silicates and other materials, including hydrocarbons producing an absorption feature at 3.4$μ$m. The detailed shape of the 10$μ$m polarization prof… ▽ More The dielectric function of interstellar dust material is modeled using observations of extinction and polarization in the infrared, together with estimates for the mass of interstellar dust. The "astrodust" material is assumed to be a mix of amorphous silicates and other materials, including hydrocarbons producing an absorption feature at 3.4$μ$m. The detailed shape of the 10$μ$m polarization profile depends on the assumed porosity and grain shape, but the 10$μ$m spectropolarimetric data are not yet good enough to clearly favor one shape over another, nor to constrain the porosity. The expected 3.4$μ$m feature polarization is consistent with existing upper limits, provided the 3.4$μ$m absorption is preferentially located in grain surface layers; a separate population of non-aligned carbonaceous grains is not required. We predict the 3.4$μ$m polarization feature to be $(Δp)_{3.4μ{\rm m}}/p(10μ{\rm m})\approx 0.016$, just below current upper limits. Polarization by the same grains at submm wavelengths is also calculated. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: 32 pages, 10 figures. Submitted to Ap.J

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

Authors: Brandon S. Hensley, B. T. Draine

Abstract: 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 p… ▽ More 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. △ Less

Submitted 18 January, 2021; v1 submitted 31 August, 2020; originally announced September 2020.

Comments: 39 pages, 14 figures. Updated to accepted version

Journal ref: 2021 ApJ, 906, 73

Turbulent dissipation, CH$^+$ abundance, H$_2$ line luminosities, and polarization in the cold neutral medium

Authors: Eric R. Moseley, B. T. Draine, Kengo Tomida, James M. Stone

Abstract: In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^{+}$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We… ▽ More In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^{+}$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We carry out simulations of MHD turbulence in molecular clouds including heating and cooling, and post-process them to study H$_2$ line emission and hot-gas chemistry, particularly the formation of CH$^+$. We explore multiple magnetic field strengths and equations of state. We use a new H$_2$ cooling function for $n_{\rm H} \leq 10^5\,{\rm cm}^{-3}$, $T\leq 5000\,{\rm K}$, and variable H$_2$ fraction. We make two important simplifying assumptions: (i) the ${\rm H}_2/{\rm H}$ fraction is fixed everywhere, and (ii) we exclude from our analysis regions where the ion-neutral drift velocity is calculated to be greater than 5 km/s. Our models produce H$_2$ emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic r.m.s. magnetic field strengths ($\approx 10$ $μ$G) and velocity dispersions, we reproduce observed CH$^+$ abundances. These findings contrast with those of Valdivia et al. (2017). Comparison of predicted dust polarization with observations by {\it Planck} suggests that the mean field $\gtrsim 5 μ$G, so that the turbulence is sub-Alfvénic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH$^+$ abundance. △ Less

Submitted 23 October, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 20 pages, 18 figures, accepted for publication to MNRAS

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

Authors: Brandon S. Hensley, B. T. Draine

Abstract: 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… ▽ More 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. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: Submitted to ApJ, comments welcome

[CI](1-0) and [CI](2-1) in resolved local galaxies

Authors: Alison F. Crocker, Eric Pellegrini, J. -D. T. Smith, Bruce T. Draine, Christine D. Wilson, Mark Wolfire, Lee Armus, Elias Brinks, Daniel A. Dale, Brent Groves, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Eric J. Murphy, Karin Sandstrom, Eva Schinnerer, Dimitra Rigopoulou, Erik Rosolowsky, Paul van der Werf

Abstract: We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4… ▽ More We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $α_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $α_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $α_{\mathrm{[CI](1-0)}}$ to $α_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $α_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: 15 pages, 10 figures, published in ApJ

Journal ref: The Astrophysical Journal 887 (2019) 105

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

Authors: G. Aniano, B. T. Draine, L. K. Hunt, K. Sandstrom, D. Calzetti, R. C. Kennicutt, D. A. Dale, M. Galametz, K. D. Gordon, A. K. Leroy, J. -D. T. Smith, H. Roussel, M. Sauvage, F. Walter, L. Armus, A. D. Bolatto, M. Boquien, A. Crocker, I. De Looze, J. Donovan Meyer, G. Helou, J. Hinz, B. D. Johnson, J. Koda, A. Miller , et al. (8 additional authors not shown)

Abstract: 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 starli… ▽ More 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. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: ApJ, accepted

Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator

Authors: Jessica Sutter, Daniel A. Dale, Kevin V. Croxall, Eric W. Pelligrini, J. D. T. Smith, Phillip N. Appleton, Pedro Beirao, Alberto D. Bolatto, Daniela Calzetti, Alison Crocker, Ilse De Looze, Bruce Draine, Maud Galametz, Brent A. Groves, George Helou, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Helene Roussel, Mark G. Wolfire

Abstract: The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that origi… ▽ More The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 19 pages, 8 figures, 8 tables

Astro2020 Science White Paper: Probing the Structure of Interstellar Dust from Micron to Kpc Scales with X-ray Imaging

Authors: Lynne Valencic, Lia Corrales, Sebastian Heinz, Randall K. Smith, Geoffrey C. Clayton, Elisa Costantini, Bruce Draine, Julia Lee, Frits Paerels, Tea Temim, Joern Wilms

Abstract: The X-ray regime is a largely underused resource for constraining interstellar dust grain models and improving our understanding of the physical processes that dictate how grains evolve over their lifetimes. This is mostly due to current detectors' relatively low sensitivity and high background, limiting the targets to the brightest sources. The improved sensitivity of the next generation of X-ray… ▽ More The X-ray regime is a largely underused resource for constraining interstellar dust grain models and improving our understanding of the physical processes that dictate how grains evolve over their lifetimes. This is mostly due to current detectors' relatively low sensitivity and high background, limiting the targets to the brightest sources. The improved sensitivity of the next generation of X-ray detectors will allow studies of much fainter sources, at much higher angular resolution, expanding our sampled sightlines in both quality and quantity. △ Less

Submitted 1 May, 2019; originally announced May 2019.

Comments: Astro2020 decadal survey science white paper submitted to the National Academy of Sciences

Astro 2020: Astromineralogy of interstellar dust with X-ray spectroscopy

Authors: Lia Corrales, Lynne Valencic, Elisa Costantini, Javier Garcia, Efrain Gatuzz, Tim Kallman, Julia Lee, Norbert Schulz, Sascha Zeegers, Claude Canizares, Bruce Draine, Sebastian Heinz, Edmund Hodges-Kluck, Edward B. Jenkins, Frits Paerels, Randall K. Smith, Tea Temim, Joern Wilms, Daniel W. Savin

Abstract: X-ray absorption fine structure (XAFS) in the 0.2-2 keV band is a crucial component in multi-wavelength studies of dust mineralogy, size, and shape -- parameters that are necessary for interpreting astronomical observations and building physical models across all fields, from cosmology to exoplanets. Despite its importance, many fundamental questions about dust remain open. What is the origin of t… ▽ More X-ray absorption fine structure (XAFS) in the 0.2-2 keV band is a crucial component in multi-wavelength studies of dust mineralogy, size, and shape -- parameters that are necessary for interpreting astronomical observations and building physical models across all fields, from cosmology to exoplanets. Despite its importance, many fundamental questions about dust remain open. What is the origin of the dust that suffuses the interstellar medium (ISM)? Where is the missing interstellar oxygen? How does iron, predominantly produced by Type Ia supernovae, become incorporated into dust? What is the main form of carbon in the ISM, and how does it differ from carbon in stellar winds? The next generation of X-ray observatories, employing microcalorimeter technology and $R \equiv λ/Δλ\geq 3000$ gratings, will provide pivotal insights for these questions by measuring XAFS in absorption and scattering. However, lab measurements of mineralogical candidates for astrophysical dust, with R > 1000, are needed to fully take advantage of the coming observations. △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: Astro2020 decadal survey science white paper submitted to the National Academy of Sciences

The Polycyclic Aromatic Hydrocarbon Mass Fraction on a 10 pc scale in the Magellanic Clouds

Authors: Jeremy Chastenet, Karin Sandstrom, I-Da Chiang, Adam Leroy, Dyas Utomo, Caroline Bot, Karl Gordon, Bruce Draine, Yasuo Fukui, Toshikazu Onishi, Kisetsu Tsuge

Abstract: We present maps of the dust properties in the Small and Large Magellanic Clouds (SMC, LMC) from fitting Spitzer and Herschel observations with the \citet{DL07} dust model. We derive the abundance of the small carbonaceous grain (or polycyclic aromatic hydrocarbon; PAH) component. The global PAH fraction (q_pah, the fraction of the dust mass in the form of PAHs) is smaller in the SMC (1.0… ▽ More We present maps of the dust properties in the Small and Large Magellanic Clouds (SMC, LMC) from fitting Spitzer and Herschel observations with the \citet{DL07} dust model. We derive the abundance of the small carbonaceous grain (or polycyclic aromatic hydrocarbon; PAH) component. The global PAH fraction (q_pah, the fraction of the dust mass in the form of PAHs) is smaller in the SMC (1.0$^{+0.3}_{-0.3}$%) than in the LMC (3.3$^{+1.4}_{-1.3}$%). We measure the PAH fraction in different gas phases (H II regions, ionized gas outside of H II regions, molecular gas, and diffuse neutral gas). H II regions appear as distinctive holes in the spatial distribution of the PAH fraction. In both galaxies, the PAH fraction in the diffuse neutral medium is higher than in the ionized gas, but similar to the molecular gas. Even at equal radiation field intensity, the PAH fraction is lower in the ionized gas than in the diffuse neutral gas. We investigate the PAH life-cycle as a function of metallicity between the two galaxies. The PAH fraction in the diffuse neutral medium of the LMC is similar to that of the Milky Way ($\sim4.6$%), while it is significantly lower in the SMC. Plausible explanations for the higher PAH fraction in the diffuse neutral medium of the LMC compared to the SMC include: a more effective PAH production by fragmentation of large grains at higher metallicity, and/or the growth of PAHs in molecular gas. △ Less

Submitted 4 April, 2019; originally announced April 2019.

Comments: Accepted for publication in ApJ

Comprehensive comparison of models for spectral energy distributions from 0.1 micron to 1 mm of nearby star-forming galaxies

Authors: L. K. Hunt, I. De Looze, M. Boquien, R. Nikutta, A. Rossi, S. Bianchi, D. A. Dale, G. L. Granato, R. C. Kennicutt, L. Silva, L. Ciesla, M. Relano, S. Viaene, B. Brandl, D. Calzetti, K. V. Croxall, B. T. Draine, M. Galametz, K. D. Gordon, B. A. Groves, G. Helou, R. Herrera-Camus, J. L. Hinz, J. Koda, S. Salim , et al. (4 additional authors not shown)

Abstract: We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of… ▽ More We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR. △ Less

Submitted 13 November, 2018; v1 submitted 11 September, 2018; originally announced September 2018.

Comments: 42 pages, 23 figures, 8 tables, 2 appendices: updated to v2, A&A, accepted

Journal ref: A&A 621, A51 (2019)

Absorption by Spinning Dust: a Contaminant for High-Redshift 21 cm Observations

Authors: B. T. Draine, Jordi Miralda-Escudé

Abstract: Spinning dust grains in front of the bright Galactic synchrotron background can produce a weak absorption signal that could affect measurements of high redshift 21 cm absorption. At frequencies near 80 MHz where the EDGES experiment has reported 21\,cm absorption at $z \approx 17$, absorption could be produced by interstellar nanoparticles with radii $a \approx 50Å$ in the cold interstellar medium… ▽ More Spinning dust grains in front of the bright Galactic synchrotron background can produce a weak absorption signal that could affect measurements of high redshift 21 cm absorption. At frequencies near 80 MHz where the EDGES experiment has reported 21\,cm absorption at $z \approx 17$, absorption could be produced by interstellar nanoparticles with radii $a \approx 50Å$ in the cold interstellar medium at temperature $T \approx 50$ K. Atmospheric aerosols could contribute additional absorption. The strength of the absorption depends on the abundance of such grains and on their dipole moments, which are uncertain. The breadth of the absorption spectrum of spinning dust limits its possible impact on measurement of a relatively narrow 21 cm absorption feature. △ Less

Submitted 7 May, 2018; v1 submitted 6 April, 2018; originally announced April 2018.

Comments: version published in ApJL

Journal ref: Astrophysical Journal Letters, 858:L10 (2018)

Electron Energy Distributions in HII Regions and Planetary Nebulae: kappa-Distributions Do Not Apply

Authors: B. T. Draine, C. D. Kreisch

Abstract: Some authors have proposed that electron energy distributions in H II regions and planetary nebulae may be significantly nonthermal, and kappa-distributions have been suggested as being appropriate. Here it is demonstrated that the electron energy distribution function is extremely close to a Maxwellian up to electron kinetic energies ~13 eV in HII regions, and up to ~16eV in planetary nebulae: ka… ▽ More Some authors have proposed that electron energy distributions in H II regions and planetary nebulae may be significantly nonthermal, and kappa-distributions have been suggested as being appropriate. Here it is demonstrated that the electron energy distribution function is extremely close to a Maxwellian up to electron kinetic energies ~13 eV in HII regions, and up to ~16eV in planetary nebulae: kappa-distributions are inappropriate. The small departures from a Maxwellian have negligible effects on line ratios. When observed line ratios in H II regions deviate from models with a single electron temperature, it must arise from spatial variations in electron temperature, rather than local deviations from a Maxwellian. △ Less

Submitted 20 July, 2018; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: published in ApJ

Journal ref: Astrophysical Journal 862:30 (2018)

The State-of-Play of Anomalous Microwave Emission (AME) Research

Authors: Clive Dickinson, Y. Ali-Haïmoud, A. Barr, E. S. Battistelli, A. Bell, L. Bernstein, S. Casassus, K. Cleary, B. T. Draine, R. Génova-Santos, S. E. Harper, B. Hensley, J. Hill-Valler, Thiem Hoang, F. P. Israel, L. Jew, A. Lazarian, J. P. Leahy, J. Leech, C. H. López-Caraballo, I. McDonald, E. J. Murphy, T. Onaka, R. Paladini, M. W. Peel , et al. (8 additional authors not shown)

Abstract: Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot b… ▽ More Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016. △ Less

Submitted 26 February, 2018; v1 submitted 22 February, 2018; originally announced February 2018.

Comments: Accepted for publication in New Astronomy Reviews. Summary of AME workshop held at ESTEC, The Netherlands, June 2016, 40 pages, 18 figures. Updated to approximately match published version

After The Fall: The Dust and Gas in E+A Post-Starburst Galaxies

Authors: A. Smercina, J. D. T. Smith, D. A. Dale, K. D. French, K. V. Croxall, S. Zhukovska, A. Togi, E. F. Bell, A. F. Crocker, B. T. Draine, T. H. Jarrett, C. Tremonti, Yujin Yang, A. I. Zabludoff

Abstract: The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ISM in many post-starbursts raise important questions. Are they truly quiescent and, if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infr… ▽ More The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ISM in many post-starbursts raise important questions. Are they truly quiescent and, if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infrared spectroscopic and photometric survey of 33 SDSS-selected E+A post-starbursts, aimed at resolving these questions. We find compact, warm dust reservoirs with high PAH abundances, and total gas and dust masses significantly higher than expected from stellar recycling alone. Both PAH/TIR and dust-to-burst stellar mass ratios are seen to decrease with post-burst age, indicative of the accumulating effects of dust destruction and an incipient transition to hot, early-type ISM properties. Their infrared spectral properties are unique, with dominant PAH emission, very weak nebular lines, unusually strong H$_{2}$ rotational emission, and deep ${\rm [C\, II]}$ deficits. There is substantial scatter among SFR indicators, and both PAH and TIR luminosities provide overestimates. Even as potential upper limits, all tracers show that the SFR has typically experienced a more than two order-of-magnitude decline since the starburst, and that the SFR is considerably lower than expected given both their stellar masses and molecular gas densities. These results paint a coherent picture of systems in which star formation was, indeed, rapidly truncated, but in which the ISM was $\textit{not}$ completely expelled, and is instead supported against collapse by latent or continued injection of turbulent or mechanical heating. The resulting aging burst populations provide a "high-soft" radiation field which seemingly dominates the E+As' unusual ISM energetics. △ Less

Submitted 6 June, 2018; v1 submitted 13 February, 2018; originally announced February 2018.

Comments: 44 pages, 21 figures, Appendix Table 2 updated

Spatially Resolved Dust, Gas, and Star Formation in the Dwarf Magellanic Irregular NGC4449

Authors: D. Calzetti, G. W. Wilson, B. T. Draine, H. Roussel, K. E. Johnson, M. H. Heyer, W. F. Wall, K. Grasha, A. Battisti, J. E. Andrews, A. Kirkpatrick, D. Rosa Gonzalez, O. Vega, J. Puschnig, M. Yun, G. Oestlin, A. S. Evans, Y. Tang, J. Lowenthal, D. Sanchez-Arguelles

Abstract: We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millim… ▽ More We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from ~1.5 to ~1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value ~1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples. △ Less

Submitted 17 December, 2017; originally announced December 2017.

Comments: 63 pages, 4 tables, 17 figures. Accepted for publication on the Astrophysical Journal

On the Shapes of Interstellar Grains: Modeling Extinction and Polarization by Spheroids and Continuous Distributions of Ellipsoids

Authors: B. T. Draine, Brandon S. Hensley

Abstract: Although interstellar grains are known to be aspherical, their actual shapes remain poorly constrained. We assess whether three continuous distributions of ellipsoidal shapes from the literature are suitable for describing the shapes of interstellar grains. Randomly-selected shapes from each distribution are shown as illustrations. The often-used Bohren-Huffman CDE includes a very large fraction o… ▽ More Although interstellar grains are known to be aspherical, their actual shapes remain poorly constrained. We assess whether three continuous distributions of ellipsoidal shapes from the literature are suitable for describing the shapes of interstellar grains. Randomly-selected shapes from each distribution are shown as illustrations. The often-used Bohren-Huffman CDE includes a very large fraction of extreme shapes: fully 10% of random draws have axial ratio $a_3/a_1 > 19.7$, and 5% have $a_3/a_1 > 33$. The CDE2 distribution includes a much smaller fraction of extreme shapes, and appears to be more realistic. For each of the three CDEs considered, we derive shape-averaged cross sections for extinction and polarization in the Rayleigh limit. Finally, we describe a method for "synthesizing" a dielectric function for an assumed shape or shape distribution if the actual absorption cross sections per grain volume in the Rayleigh limit are known from observations. This synthetic dielectric function predicts the wavelength dependence of polarization, which can then be compared to observations to constrain the grain shape. △ Less

Submitted 18 January, 2021; v1 submitted 24 October, 2017; originally announced October 2017.

Comments: 23 pages, 8 figures. ApJ, accepted

The Origins of [CII] Emission in Local Star-forming Galaxies

Authors: Kevin Croxall, J. D. T. Smith, Eric Pellegrini, Brent Groves, Alberto Bolatto, Rodrigo Herrera-Camus, Karin Sandstrom, Bruce T. Draine, Mark Wolfire, Lee Armus, Mederic Boquien, Bernhard Brandl, Daniel A. Dale, Maud Galametz, Leslie K. Hunt, Robert C. Kennicutt, Kathryn Kreckel, Dimitra Rigopoulou, Paul P. van der Werf, Christine D. Wilson

Abstract: The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the inters… ▽ More The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided. △ Less

Submitted 14 July, 2017; originally announced July 2017.

Comments: 8 pages, accepted for publication in ApJ

The Survey of Lines in M31 (SLIM): The Drivers of the [CII]/TIR Variation

Authors: Maria Julia Kapala, Brent Groves, Karin Sandstrom, Thomas Jarrett, Elisabete da Cunha, Kevin Croxall, Julianne Dalcanton, Bruce Draine, Simon Glover, Eva Schinnerer

Abstract: The ratio of the [CII] 158$\,μ$m emission line over the total infrared emission (TIR) is often used as a proxy for the photoelectric (PE) heating efficiency ($ε_{\rm PE}$) of the far-ultraviolet (FUV) photons absorbed by dust in the interstellar medium. In the nearby galaxy M31, we measure a strong radial variation of [CII]/TIR that we rule out as being due to an intrinsic variation in… ▽ More The ratio of the [CII] 158$\,μ$m emission line over the total infrared emission (TIR) is often used as a proxy for the photoelectric (PE) heating efficiency ($ε_{\rm PE}$) of the far-ultraviolet (FUV) photons absorbed by dust in the interstellar medium. In the nearby galaxy M31, we measure a strong radial variation of [CII]/TIR that we rule out as being due to an intrinsic variation in $ε_{\rm PE}$. [CII]/TIR fails as a proxy for $ε_{\rm PE}$, because the TIR measures all dust heating, not just the contribution from FUV photons capable of ejecting electrons from dust grains. Using extensive multiwavelength coverage from the FUV to far-infrared (FIR), we infer the attenuated FUV emission ($\rm UV_{att}$), and the total attenuated flux ($\rm TOT_{att}$). We find [CII]/TIR to be strongly correlated with $\rm UV_{att}$/$\rm TOT_{att}$, indicating that, in M31 at least, one of the dominant drivers for [CII]/TIR variation is the relative hardness of the absorbed stellar radiation field. We define $\rm{ ε_{PE}^{UV}}$, [CII]/$\rm{ UV_{att}}$ which should be more closely related to the actual PE efficiency, which we find to be essentially constant ($1.85 \pm 0.8 \%$) in all explored fields in M31. This suggests that part of the observed variation of [CII]/TIR in other galaxies is likely due to a change in the relative hardness of the absorbed stellar radiation field, caused by a combination of variations in the stellar population, dust opacity and galaxy metallicity, although PE efficiency may also vary across a wider range of environments. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: 19 pages, 16 figures, accepted for publication in ApJ

Thermal Pressure in the Cold Neutral Medium of Nearby Galaxies

Authors: R. Herrera-Camus, A. Bolatto, M. Wolfire, E. Ostriker, B. Draine, A. Leroy, K. Sandstrom, L. Hunt, R. Kennicutt, D. Calzetti, J. D. Smith, K. Croxall, M. Galametz, I. de Looze, D. Dale, A. Crocker, B. Groves

Abstract: Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $μ$m, HI, and CO data… ▽ More Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $μ$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks. △ Less

Submitted 5 January, 2017; originally announced January 2017.

Comments: Accepted for publication in ApJ on Dec. 20th, 2016; 17 pages; 7 figures and 1 table

Thermodynamics and Charging of Interstellar Iron Nanoparticles

Authors: Brandon S. Hensley, B. T. Draine

Abstract: Interstellar iron in the form of metallic iron nanoparticles may constitute a component of the interstellar dust. We compute the stability of iron nanoparticles to sublimation in the interstellar radiation field, finding that iron clusters can persist down to a radius of $\simeq 4.5\,$Å, and perhaps smaller. We employ laboratory data on small iron clusters to compute the photoelectric yields as a… ▽ More Interstellar iron in the form of metallic iron nanoparticles may constitute a component of the interstellar dust. We compute the stability of iron nanoparticles to sublimation in the interstellar radiation field, finding that iron clusters can persist down to a radius of $\simeq 4.5\,$Å, and perhaps smaller. We employ laboratory data on small iron clusters to compute the photoelectric yields as a function of grain size and the resulting grain charge distribution in various interstellar environments, finding that iron nanoparticles can acquire negative charges particularly in regions with high gas temperatures and ionization fractions. If $\gtrsim 10\%$ of the interstellar iron is in the form of ultrasmall iron clusters, the photoelectric heating rate from dust may be increased by up to tens of percent relative to dust models with only carbonaceous and silicate grains. △ Less

Submitted 25 November, 2016; originally announced November 2016.

Comments: Accepted to ApJ

Modeling the Anomalous Microwave Emission with Spinning Nanoparticles: No PAHs Required

Authors: Brandon S. Hensley, B. T. Draine

Abstract: In light of recent observational results indicating an apparent lack of correlation between the Anomalous Microwave Emission (AME) and mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs), we assess whether rotational emission from spinning silicate and/or iron nanoparticles could account for the observed AME without violating observational constraints on interstellar abundances, ult… ▽ More In light of recent observational results indicating an apparent lack of correlation between the Anomalous Microwave Emission (AME) and mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs), we assess whether rotational emission from spinning silicate and/or iron nanoparticles could account for the observed AME without violating observational constraints on interstellar abundances, ultraviolet extinction, and infrared emission. By modifying the SpDust code to compute the rotational emission from these grains, we find that nanosilicate grains could account for the entirety of the observed AME, whereas iron grains could be responsible for only a fraction, even for extreme assumptions on the amount of interstellar iron concentrated in ultrasmall iron nanoparticles. Given the added complexity of contributions from multiple grain populations to the total spinning dust emission, as well as existing uncertainties due to the poorly-constrained grain size, charge, and dipole moment distributions, we discuss generic, carrier-independent predictions of spinning dust theory and observational tests that could help identify the AME carrier(s). △ Less

Submitted 25 November, 2016; originally announced November 2016.

Comments: Submitted to ApJ, comments welcome

The Spatially Resolved [CII] Cooling Line Deficit in Galaxies

Authors: J. D. T. Smith, Kevin Croxall, Bruce Draine, Ilse De Looze, Karin Sandstrom, Lee Armus, Pedro Beirao, Alberto Bolatto, Mederic Boquien, Bernhard Brandl, Alison Crocker, Daniel A. Dale, Maud Galametz, Brent Groves, George Helou, Rodrigo Herrera-Camus, Leslie Hunt, Robert Kennicutt, Fabian Walter, Mark Wolfire

Abstract: We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominat… ▽ More We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas. △ Less

Submitted 4 November, 2016; originally announced November 2016.

Comments: 17 pages, 5 figures, accepted for publication in ApJ

Graphite Revisited

Authors: B. T. Draine

Abstract: 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… ▽ More 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. △ Less

Submitted 9 August, 2016; originally announced August 2016.

Comments: To appear in ApJ

Radiative Transfer Model of Dust Attenuation Curves in Clumpy, Galactic Environments

Authors: Kwang-Il Seon, Bruce T. Draine

Abstract: The attenuation of starlight by dust in galactic environments is investigated through models of radiative transfer in a spherical, clumpy ISM. Extinction properties for MW, LMC, and SMC dust types are considered. It is illustrated that the attenuation curves are primarily determined by the wavelength dependence of absorption rather than by the underlying extinction (absorption+scattering) curve. A… ▽ More The attenuation of starlight by dust in galactic environments is investigated through models of radiative transfer in a spherical, clumpy ISM. Extinction properties for MW, LMC, and SMC dust types are considered. It is illustrated that the attenuation curves are primarily determined by the wavelength dependence of absorption rather than by the underlying extinction (absorption+scattering) curve. Attenuation curves consistent with the "Calzetti attenuation curve" are found by assuming the silicate-carbonaceous dust model for the MW, but with the 2175A absorption bump suppressed or absent. The discrepancy between our results and previous work that claimed the SMC-type dust to be the most likely origin of the Calzetti curve is ascribed to the difference in adopted albedos; this study uses the theoretically calculated albedos whereas the previous ones adopted empirically derived albedos from observations of reflection nebulae. It is also found that the model attenuation curves calculated with the MW dust are well represented by a modified Calzetti curve with a varying slope and UV bump strength. The strong correlation between the slope and UV bump strength, with steeper curves having stronger bumps, as found in star-forming galaxies at 0.5<z<2.0, is well reproduced by our models if the abundance of the UV bump carriers or PAHs is assumed to be 30% or 40% of that of the MW-dust. The trend is explained by radiative transfer effects which lead to shallower attenuation curves with weaker UV bumps as the ISM is more clumpy and dustier. We also argue that at least some of the IUE local starburst galaxies may have a UV bump feature in their attenuation curves, albeit much weaker than that of the MW extinction curve. △ Less

Submitted 7 June, 2016; originally announced June 2016.

Comments: 28 pages, 30 figures, submited to ApJS

Quantum Suppression of Alignment in Ultrasmall Grains: Microwave Emission from Spinning Dust will be Negligibly Polarized

Authors: B. T. Draine, Brandon S. Hensley

Abstract: The quantization of energy levels in very nanoparticles suppresses dissipative processes that convert grain rotational kinetic energy into heat. For grains small enough to have GHz rotation rates, the suppression of dissipation can be extreme. As a result, alignment of such grains is suppressed. This applies both to alignment of the grain body with its angular momentum J, and to alignment of J wit… ▽ More The quantization of energy levels in very nanoparticles suppresses dissipative processes that convert grain rotational kinetic energy into heat. For grains small enough to have GHz rotation rates, the suppression of dissipation can be extreme. As a result, alignment of such grains is suppressed. This applies both to alignment of the grain body with its angular momentum J, and to alignment of J with the local magnetic field B_0. If the anomalous microwave emission is rotational emission from spinning grains, it will be negligibly polarized at GHz frequencies, with P < 10^{-6} at frequencies above 10 GHz. △ Less

Submitted 21 May, 2016; originally announced May 2016.

Comments: submitted to ApJ

The Ionized Gas in Nearby Galaxies as Traced by the [NII] 122 and 205 μm Transitions

Authors: R. Herrera-Camus, A. Bolatto, J. D. Smith, B. Draine, E. Pellegrini, M. Wolfire, K. Croxall, I. de Looze, D. Calzetti, R. Kennicutt, A. Crocker, L. Armus, P. van der Werf, K. Sandstrom, M. Galametz, B. Brandl, B. Groves, D. Rigopoulou, F. Walter, A. Leroy, M. Boquien, F. S. Tabatabaei, P. Beirao

Abstract: The [NII] 122 and 205 μm transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range,… ▽ More The [NII] 122 and 205 μm transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($Σ_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $Σ_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code. △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: 18 pages, 9 figures, accepted for publication in The Astrophysical Journal

Towards universal hybrid star formation rate estimators

Authors: M. Boquien, R. Kennicutt, D. Calzetti, D. Dale, M. Galametz, M. Sauvage, K. Croxall, B. Draine, A. Kirkpatrick, N. Kumari, L. Hunt, I. De Looze, E. Pellegrini, M. Relano, J. -D. Smith, F. Tabatabaei

Abstract: To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quant… ▽ More To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quantifying what physical processes drive the variations between different hybrid estimators. Doing so, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation, taking into account the intrinsic physical properties of galaxies. We use the CIGALE code to model the spatially-resolved FUV to FIR SED of eight nearby star-forming galaxies drawn from the KINGFISH sample. This allows us to determine their local physical properties, and in particular their UV attenuation, average SFR, average specific SFR (sSFR), and their stellar mass. We then examine how hybrid estimators depend on said properties. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 and 100 micron) and on the sSFR (in particular at 24 micron and the TIR). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of hybrid estimators based on observed UV to near-IR colours and near-IR luminosity densities per unit area. We find that they can reliably be extended to entire galaxies. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators. Naturally taking into account the variable impact of dust heated by old stellar populations, they constitute a step towards universal estimators. △ Less

Submitted 30 March, 2016; originally announced March 2016.

Comments: 25 pages, 12 figures, 8 tables. Accepted for publication in A&A

X-Ray Absorption and Scattering by Interstellar Grains

Authors: John A. Hoffman, Bruce T. Draine

Abstract: Interstellar abundance determinations from fits to X-ray absorption edges often rely on the incorrect assumption that scattering is insignificant and can be ignored. We show instead that scattering contributes significantly to the attenuation of X-rays for realistic dust grain size distributions and substantially modifies the spectrum near absorption edges of elements present in grains. The dust a… ▽ More Interstellar abundance determinations from fits to X-ray absorption edges often rely on the incorrect assumption that scattering is insignificant and can be ignored. We show instead that scattering contributes significantly to the attenuation of X-rays for realistic dust grain size distributions and substantially modifies the spectrum near absorption edges of elements present in grains. The dust attenuation modules used in major X-ray spectral fitting programs do not take this into account. We show that the consequences of neglecting scattering on the determination of interstellar elemental abundances are modest; however, scattering (along with uncertainties in the grain size distribution) must be taken into account when near-edge extinction fine structure is used to infer dust mineralogy. We advertise the benefits and accuracy of anomalous diffraction theory for both X-ray halo analysis and near edge absorption studies. An open source Fortran suite, General Geometry Anomalous Diffraction Theory (GGADT), is presented that calculates X-ray absorption, scattering, and differential scattering cross sections for grains of arbitrary geometry and composition. △ Less

Submitted 4 February, 2016; v1 submitted 29 September, 2015; originally announced September 2015.

Comments: 32 pages, 15 figures

Journal ref: ApJ 817 139 2016

Radiation pressure on fluffy submicron-sized grains

Authors: Kedron Silsbee, Bruce T. Draine

Abstract: We investigate the claim that the ratio β of radiation pressure force to gravitational force on a dust grain in our solar system can substantially exceed unity for some grain sizes, provided that grain porosity is high enough. For model grains consisting of random aggregates of silicate spherules, we find that the maximum value of β is almost independent of grain porosity, but for small (<0.3 μm)… ▽ More We investigate the claim that the ratio β of radiation pressure force to gravitational force on a dust grain in our solar system can substantially exceed unity for some grain sizes, provided that grain porosity is high enough. For model grains consisting of random aggregates of silicate spherules, we find that the maximum value of β is almost independent of grain porosity, but for small (<0.3 μm) grains, β actually decreases with increasing porosity. We also investigate the effect of metallic iron and amorphous carbon inclusions in the dust grains and find that while these inclusions do increase the radiation pressure cross-section, β remains below unity for grains with 3 pg of silicate material. These results affect the interpretation of the grain trajectories estimated from the {\it Stardust} mission, which were modeled assuming β values exceeding one. We find that radiation pressure effects are not large enough for particles Orion and Hylabrook captured by {\it Stardust} to be of interstellar origin given their reported impact velocities. We also consider the effects of solar radiation on transverse velocities and grain spin, and show that radiation pressure introduces both transverse velocities and equatorial spin velocities of several hundred meters per second for incoming interstellar grains at 2 AU. These transverse velocities are not important for modeling trajectories, but such spin rates may result in centrifugal disruption of aggregates. △ Less

Submitted 7 March, 2016; v1 submitted 3 August, 2015; originally announced August 2015.

Journal ref: ApJ, 818, 133 (2016)

Spatially Resolved Spitzer-IRS Spectral Maps of the Superwind in M82

Authors: P. Beirão, L. Armus, M. D. Lehnert, P. Guillard, T. Heckman, B. Draine, D. Hollenbach, F. Walter, K. Sheth, J. D. Smith, P. Shopbell, F. Boulanger, J. Surace, C. Hoopes, C. Engelbracht

Abstract: We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and… ▽ More We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and $E_{K}\sim6-20\times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAH's with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow. △ Less

Submitted 15 May, 2015; originally announced May 2015.

Comments: 19 pages and 12 figures; accepted in MNRAS

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

Authors: Brandon S. Hensley, B. T. Draine, Aaron M. Meisner

Abstract: 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 uncorrelate… ▽ More 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. △ Less

Submitted 26 October, 2015; v1 submitted 8 May, 2015; originally announced May 2015.

Comments: Submitted to ApJ, comments welcome

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

Authors: T. P. R. van der Laan, L. Armus, P. Beirao, K. Sandstrom, B. Groves, E. Schinnerer, B. T. Draine, J. D. Smith, M. Galametz, M. Wolfire, K. Croxall, D. Dale, R. Herrera Camus, D. Calzetti, R. C. Kennicutt Jr

Abstract: 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 ga… ▽ More 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. △ Less

Submitted 28 January, 2015; originally announced January 2015.

Comments: accepted for publication in A&A

Journal ref: A&A 575, A83 (2015)

The Composition of Interstellar Grains Toward Zeta Ophiuchi: Constraining the Elemental Budget Near the Diffuse-Dense Cloud Transition

Authors: Charles A. Poteet, Douglas C. B. Whittet, Bruce T. Draine

Abstract: We investigate the composition of interstellar grains along the line of sight toward Zeta Ophiuchi, a well-studied environment near the diffuse-dense cloud transition. A spectral decomposition analysis of the solid-state absorbers is performed using archival spectroscopic observations from the Spitzer Space Telescope and Infrared Space Observatory. We find strong evidence for the presence of sub-m… ▽ More We investigate the composition of interstellar grains along the line of sight toward Zeta Ophiuchi, a well-studied environment near the diffuse-dense cloud transition. A spectral decomposition analysis of the solid-state absorbers is performed using archival spectroscopic observations from the Spitzer Space Telescope and Infrared Space Observatory. We find strong evidence for the presence of sub-micron-sized amorphous silicate grains, principally comprised of olivine-like composition, with no convincing evidence of H2O ice mantles. However, tentative evidence for thick H2O ice mantles on large (a ~ 2.8 microns) grains is presented. Solid-state abundances of elemental Mg, Si, Fe, and O are inferred from our analysis and compared to standard reference abundances. We find that nearly all of elemental Mg and Si along the line of sight are present in amorphous silicate grains, while a substantial fraction of elemental Fe resides in compounds other than silicates. Moreover, we find that the total abundance of elemental O is largely inconsistent with the adopted reference abundances, indicating that as much as ~156 ppm of interstellar O is missing along the line of sight. After taking into account additional limits on the abundance of elemental O in other O-bearing solids, we conclude that any missing reservoir of elemental O must reside on large grains that are nearly opaque to infrared radiation. △ Less

Submitted 12 January, 2015; originally announced January 2015.

Comments: Accepted for publication in the Astrophysical Journal, 2015 January 12: 13 pages, 8 figures, 4 tables

Studies of Anomalous Microwave Emission (AME) with the SKA

Authors: Clive Dickinson, Yacine Ali-Haimoud, Robert J. Beswick, Simon Casassus, Kieran Cleary, Bruce T. Draine, Ricardo Genova-Santos, Keith Grainge, Thiem C. Hoang, Alex Lazarian, Eric J. Murphy, Roberta Paladini, Michael W. Peel, Yvette Perrott, Jose-Alberto Rubino-Martin, Anna Scaife, Chris T. Tibbs, Laurent Verstraete, Matias Vidal, Robert A. Watson, Nathalie Ysard

Abstract: In this chapter, we will outline the scientific motivation for studying Anomalous Microwave Emission (AME) with the SKA. AME is thought to be due to electric dipole radiation from small spinning dust grains, although thermal fluctuations of magnetic dust grains may also contribute. Studies of this mysterious component would shed light on the emission mechanism, which then opens up a new window ont… ▽ More In this chapter, we will outline the scientific motivation for studying Anomalous Microwave Emission (AME) with the SKA. AME is thought to be due to electric dipole radiation from small spinning dust grains, although thermal fluctuations of magnetic dust grains may also contribute. Studies of this mysterious component would shed light on the emission mechanism, which then opens up a new window onto the interstellar medium (ISM). AME is emitted mostly in the frequency range $\sim 10$--100\,GHz, and thus the SKA has the potential of measuring the low frequency side of the AME spectrum, particularly in band 5. Science targets include dense molecular clouds in the Milky Way, as well as extragalactic sources. We also discuss the possibility of detecting rotational line emission from Poly-cyclic Aromatic Hydrocarbons (PAHs), which could be the main carriers of AME. Detecting PAH lines of a given spacing would allow for a definitive identification of specific PAH species. △ Less

Submitted 16 December, 2014; originally announced December 2014.

Comments: 11 pages, 4 pages, to be published in Proceedings of Science: "Advancing Astrophysics with the Square Kilometre Array"

[CII] 158 $μ$m Emission as a Star Formation Tracer

Authors: R. Herrera-Camus, A. D. Bolatto, M. G. Wolfire, J. D. Smith, K. V. Croxall, R. C. Kennicutt, D. Calzetti, G. Helou, F. Walter, A. K. Leroy, B. Draine, B. R. Brandl, L. Armus, K. M. Sandstrom, D. A. Dale, G. Aniano, S. E. Meidt, M. Boquien, L. K. Hunt, M. Galametz, F. S. Tabatabaei, E. J. Murphy, P. Appleton, H. Roussel, C. Engelbracht , et al. (1 additional authors not shown)

Abstract: The [CII] 157.74 $μ$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal st… ▽ More The [CII] 157.74 $μ$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $Σ_{\rm [CII]}-Σ_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $Σ_{\rm[CII]}-Σ_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $Σ_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies. △ Less

Submitted 24 September, 2014; originally announced September 2014.

Comments: 25 pages, 14 figures, accepted for publication in ApJ