Ultraviolet Dust Grain Properties in Starburst Galaxies: Evidence from Radiative Transfer Modeling and Local Group Extinction Curves

This paper summarizes the evidence of the ultraviolet properties of dust grains found in starburst galaxies. Observations of starburst galaxies clearly show that the 2175 A feature is weak or absent. This can be the result of radiative transfer effects (mixing the dust and stars) or due to dust grains which do not have this feature. Spherical DIRTY radiative transfer models imply that it is not radiative transfer effects, but other radiative transfer models with disk/bulge geometries have found cases where it could be radiative transfer effects. Recent work on the extinction curves in the Magellanic Clouds and Milky Way has revealed that the traditional explanation of low metallicity for the absence of the 2175 A feature in the Small Magellanic Cloud is likely incorrect. The SMC has one sightline with a 2175 A feature and the Milky Way has sightlines without this feature. In addition, where the 2175 A feature is found to be weak or absent in both Magellanic Clouds and the Milky Way, there is evidence for recent star formation. Taking the sum of the radiative transfer modeling of starburst galaxies and the behavior of Local Group extinction curves, it is likely that the dust grains in starburst galaxies intrinsically lack the 2175 A feature.Comment: 7 pages, To appear in the proceedings of: "The Spectral Energy Distribution of Gas-Rich Galaxies: Confronting Models with Data", Heidelberg, 4-8 Oct. 2004, eds. C.C. Popescu and R.J. Tuffs, AIP Conf. Ser., in press [fixed typo in title

Young, UV-bright Stars Dominate Dust Heating in Star Forming Galaxies

In star forming galaxies, dust plays a significant role in shaping the ultraviolet (UV) through infrared (IR) spectrum. Dust attenuates the radiation from stars, and re-radiates the energy through equilibrium and non-equilibrium emission. Polycyclic aromatic hydrocarbons (PAH), graphite, and silicates contribute to different features in the spectral energy distribution; however, they are all highly opaque in the same spectral region - the UV. Compared to old stellar populations, young populations release a higher fraction of their total luminosity in the UV, making them a good source of the energetic UV photons that can power dust emission. However, given their relative abundance, the question of whether young or old stellar populations provide most of these photons that power the infrared emission is an interesting question. Using three samples of galaxies observed with the Spitzer Space Telescope and our dusty radiative transfer model, we find that young stellar populations (on the order of 100 million years old) dominate the dust heating in star forming galaxies, and old stellar populations (13 billion years old) generally contribute less than 20% of the far-IR luminosity.Comment: 12 pages, 8 figures, accepted to the Ap

Optimal Cosmic-Ray Detection for Nondestructive Read Ramps

Cosmic rays are a known problem in astronomy, causing both loss of data and data inaccuracy. The problem becomes even more extreme when considering data from a high-radiation environment, such as in orbit around Earth or outside the Earth's magnetic field altogether, unprotected, as will be the case for the James Webb Space Telescope (JWST). For JWST, all the instruments employ nondestructive readout schemes. The most common of these will be "up the ramp" sampling, where the detector is read out regularly during the ramp. We study three methods to correct for cosmic rays in these ramps: a two-point difference method, a deviation from the fit method, and a y-intercept method. We apply these methods to simulated nondestructive read ramps with single-sample groups and varying combinations of flux, number of samples, number of cosmic rays, cosmic-ray location in the exposure, and cosmic-ray strength. We show that the y-intercept method is the optimal detection method in the read-noise-dominated regime, while both the y-intercept method and the two-point difference method are best in the photon-noise-dominated regime, with the latter requiring fewer computations.Comment: To be published in PASP. This paper is 12 pages long and includes 15 figure

Discovery of Blue Luminescence in the Red Rectangle: Possible Fluorescence from Neutral Polycyclic Aromatic Hydrocarbon Molecules?

Here we report our discovery of a band of blue luminescence (BL) in the Red Rectangle (RR) nebula. This enigmatic proto-planetary nebula is also one of the brightest known sources of extended red emission as well as of unidentified infra-red (UIR) band emissions. The spectrum of this newly discovered BL is most likely fluorescence from small neutral polycyclic aromatic hydrocarbon (PAH) molecules. PAH molecules are thought to be widely present in many interstellar and circumstellar environments in our galaxy as well as in other galaxies, and are considered likely carriers of the UIR-band emission. However, no specific PAH molecule has yet been identified in a source outside the solar system, as the set of mid-infra-red emission features attributed to these molecules between the wavelengths of 3.3 micron and 16.4 micron is largely insensitive to molecular sizes. In contrast, near-UV/blue fluorescence of PAHs is more specific as to size, structure, and charge state of a PAH molecule. If the carriers of this near-UV/blue fluorescence are PAHs, they are most likely neutral PAH molecules consisting of 3-4 aromatic rings such as anthracene (C14H10) and pyrene (C16H10). These small PAHs would then be the largest molecules specifically identified in the interstellar medium.Comment: 4 pages, 4 figures, Accepted for publication in ApJL (LaTeX, uses emulateapj.sty

Composite biasing in Monte Carlo radiative transfer

Biasing or importance sampling is a powerful technique in Monte Carlo radiative transfer, and can be applied in different forms to increase the accuracy and efficiency of simulations. One of the drawbacks of the use of biasing is the potential introduction of large weight factors. We discuss a general strategy, composite biasing, to suppress the appearance of large weight factors. We use this composite biasing approach for two different problems faced by current state-of-the-art Monte Carlo radiative transfer codes: the generation of photon packages from multiple components, and the penetration of radiation through high optical depth barriers. In both cases, the implementation of the relevant algorithms is trivial and does not interfere with any other optimisation techniques. Through simple test models, we demonstrate the general applicability, accuracy and efficiency of the composite biasing approach. In particular, for the penetration of high optical depths, the gain in efficiency is spectacular for the specific problems that we consider: in simulations with composite path length stretching, high accuracy results are obtained even for simulations with modest numbers of photon packages, while simulations without biasing cannot reach convergence, even with a huge number of photon packages.Comment: 12 pages, accepted for publication in A&

The Flux Ratio Method for Determining the Dust Attenuation of Starburst Galaxies

The presence of dust in starburst galaxies complicates the study of their stellar populations as the dust's effects are similar to those associated with changes in the galaxies' stellar age and metallicity. This degeneracy can be overcome for starburst galaxies if UV/optical/near-infrared observations are combined with far-infrared observations. We present the calibration of the flux ratio method for calculating the dust attenuation at a particular wavelength, Att(\lambda), based on the measurement of F(IR)/F(\lambda) flux ratio. Our calibration is based on spectral energy distributions (SEDs) from the PEGASE stellar evolutionary synthesis model and the effects of dust (absorption and scattering) as calculated from our Monte Carlo radiative transfer model. We tested the attenuations predicted from this method for the Balmer emission lines of a sample starburst galaxies against those calculated using radio observations and found good agreement. The UV attenuation curves for a handful of starburst galaxies were calculated using the flux ratio method, and they compare favorably with past work. The relationship between Att(\lambda) and F(IR)/F(\lambda) is almost completely independent of the assumed dust properties (grain type, distribution, and clumpiness). For the UV, the relationship is also independent of the assumed stellar properties (age, metallicity, etc) accept for the case of very old burst populations. However at longer wavelengths, the relationship is dependent on the assumed stellar properties.Comment: accepted by the ApJ, 18 pages, color figures, b/w version at http://mips.as.arizona.edu/~kgordon/papers/fr_method.htm

Measuring Extinction Curves of Lensing Galaxies

We critique the method of constructing extinction curves of lensing galaxies using multiply imaged QSOs. If one of the two QSO images is lightly reddened or if the dust along both sightlines has the same properties then the method works well and produces an extinction curve for the lensing galaxy. These cases are likely rare and hard to confirm. However, if the dust along each sightline has different properties then the resulting curve is no longer a measurement of extinction. Instead, it is a measurement of the difference between two extinction curves. This "lens difference curve'' does contain information about the dust properties, but extracting a meaningful extinction curve is not possible without additional, currently unknown information. As a quantitative example, we show that the combination of two Cardelli, Clayton, & Mathis (CCM) type extinction curves having different values of R(V) will produce a CCM extinction curve with a value of R(V) which is dependent on the individual R(V) values and the ratio of V band extinctions. The resulting lens difference curve is not an average of the dust along the two sightlines. We find that lens difference curves with any value of R(V), even negative values, can be produced by a combination of two reddened sightlines with different CCM extinction curves with R(V) values consistent with Milky Way dust (2.1 < R(V) < 5.6). This may explain extreme values of R(V) inferred by this method in previous studies. But lens difference curves with more normal values of R(V) are just as likely to be composed of two dust extinction curves with R(V) values different than that of the lens difference curve. While it is not possible to determine the individual extinction curves making up a lens difference curve, there is information about a galaxy's dust contained in the lens difference curves.Comment: 15 pages, 4 figues, ApJ in pres