Local Molecular Clouds Research Articles

HESS J1809-193: Gamma-Ray Emission by Cosmic Rays from a Past Explosion

The very-high-energy γ-ray source HESS J1809-193 has been detected by the LHAASO and HAWC observatory beyond 100 TeV energy. It is an interesting candidate for exploring the underlying mechanisms of γ-ray production due to the presence of supernova remnants, pulsars, and molecular clouds close to it. We have considered the injection of the energetic cosmic rays from a past explosion, whose reminiscent may be SNR G011.0-00.0, which is located within the extended γ-ray source HESS J1809-193. We explain the multiwavelength data from the region of HESS J1809-193 with synchrotron, inverse Compton, and bremsstrahlung emission of cosmic-ray electrons and secondary γ-ray production in interactions of cosmic-ray protons with the cold protons in the local molecular clouds within a time-dependent framework including the diffusion loss of cosmic rays. The observational data have been modeled with the secondary photons produced by the time-evolved cosmic-ray spectrum assuming the age of the explosion is 4500 yr.

Neutral Atomic and Molecular Clouds and Star Formation in the Outer Carina Arm

We present a comprehensive investigation of H i (super)clouds, molecular clouds (MCs), and star formation in the Carina spiral arm of the outer Galaxy. Utilizing HI4PI and CfA CO survey data, we identify H i clouds and MCs based on the (l, v LSR) locations of the Carina arm. We analyzed 26 H i clouds and 48 MCs. Most of the identified H i clouds are superclouds, with masses exceeding 106 M ⊙. We find that 15 of these superclouds have associated MC(s) with M H I ≳ 106 M ⊙ and 50 M ⊙ pc−2. Our virial equilibrium analysis suggests that these CO-bright H i clouds are gravitationally bound or marginally bound. We report an anticorrelation between molecular mass fractions and Galactocentric distances, and a correlation with total gas surface densities. Nine CO-bright H i superclouds are associated with H ii regions, indicating ongoing star formation. We confirm the regular spacing of H i superclouds along the spiral arm, which is likely due to some underlying physical process, such as gravitational instabilities. We observe a strong spatial correlation between H ii regions and MCs, with some offsets between MCs and local H i column density peaks. Our study reveals that, in the context of H i superclouds, the star formation rate surface density is independent of H i and total gas surface densities but positively correlates with molecular gas surface density. This finding is consistent with both extragalactic studies of the resolved Kennicutt–Schmidt relation and local giant molecular clouds study of Lada et al. (2013), emphasizing the crucial role of molecular gas in regulating star formation processes.

Bimodal star formation in simulations of strongly magnetized giant molecular clouds

ABSTRACT We present the results of a set of radiation magnetohydrodynamic simulations of turbulent molecular clouds in which we vary the initial strength of the magnetic field within a range (1 ≲ μ ≲ 5) consistent with observations of local giant molecular clouds (GMCs). We find that as we increase the strength of the magnetic field, star formation transitions from unimodal (the baseline case, μ = 5, with a single burst of star formation and Salpeter IMF) to bimodal. This effect is clearest in the most strongly magnetized GMC (μ = 1): a first burst of star formation with duration, intensity, and IMF comparable to the baseline case, is followed by a second star formation episode in which only low-mass stars are formed. Overall, due to the second burst of star formation, the strongly magnetized case results in a longer star formation period and a higher efficiency of star formation. The second burst is produced by gas that is not expelled by radiative feedback, instead remains trapped in the GMC by the large-scale B-field, producing a nearly one-dimensional flow of gas along the field lines. The trapped gas has a turbulent and magnetic topology that differs from that of the first phase and strongly suppresses gas accretion onto protostellar cores, reducing their masses. We speculate that this star formation bimodality may be an important ingredient to understand the origin of multiple stellar populations observed in massive globular clusters.

The effect of viewing angle on the Kennicutt-Schmidt relation of the local molecular clouds

The Gaia data give us an unprecedented view to the three-dimensional (3D) structure of molecular clouds in the solar neighbourhood. We study how the projected areas and masses of clouds, and consequently the Kennicutt-Schmidt (KS) relation, depend on the viewing angle. We derive the probability distributions of the projected areas and masses for nine clouds within 400 pc of the Sun using 3D dust distribution data from the literature. We find that the viewing angle can have a dramatic effect on the observed areas and masses of individual clouds. The joint probability distributions of the areas and masses are strongly correlated, relatively flat, and can show multiple peaks. The typical ranges and 50% quartiles of the distributions are roughly 100–200% and 20–80% of the median value, respectively, making viewing angle effects important for all individual clouds. The threshold value used to define the cloud areas is also important; our analysis suggests that the clouds become more anisotropic for smaller thresholds (larger scales). On average, the areas and masses of the plane-of-the-sky and face-on projections agree, albeit with a large scatter. This suggests that sample averages of areas and masses are relatively free of viewing angle effects, which is important to facilitate comparisons of extragalactic and Galactic data. Ultimately, our results demonstrate that a cloud’s location in the KS relation is affected by the viewing angle in a non-trivial manner. However, the KS relation of our sample as a whole is not strongly affected by these effects, because the covariance of the areas and masses causes the observed mean column density to remain relatively constant.

On the Three-dimensional Structure of Local Molecular Clouds

Abstract We leverage the 1 pc spatial resolution of the Leike et al. three-dimensional (3D) dust map to characterize the 3D structure of nearby molecular clouds (d ≲ 400 pc). We start by “skeletonizing” the clouds in 3D volume density space to determine their “spines,” which we project on the sky to constrain cloud distances with ≈1% uncertainty. For each cloud, we determine an average radial volume density profile around its 3D spine and fit the profiles using Gaussian and Plummer functions. The radial volume density profiles are well described by a two-component Gaussian function, consistent with clouds having broad, lower-density outer envelopes and narrow, higher-density inner layers. The ratio of the outer to inner envelope widths is ≈3:1. We hypothesize that these two components may be tracing a transition between atomic and diffuse molecular gas or between the unstable and cold neutral medium. Plummer-like models can also provide a good fit, with molecular clouds exhibiting shallow power-law wings with density, n, falling off like n −2 at large radii. Using Bayesian model selection, we find that parameterizing the clouds’ profiles using a single Gaussian is disfavored. We compare our results with two-dimensional dust extinction maps, finding that the 3D dust recovers the total cloud mass from integrated approaches with fidelity, deviating only at higher levels of extinction (A V ≳ 2–3 mag). The 3D cloud structure described here will enable comparisons with synthetic clouds generated in simulations, offering unprecedented insight into the origins and fates of molecular clouds in the interstellar medium.

The MOSDEF Survey: The First Direct Measurements of the Nebular Dust Attenuation Curve at High Redshift*

Abstract We use a sample of 532 star-forming galaxies at redshifts z ≃ 1.4–2.6 with deep rest-frame optical spectra from the MOSFIRE Deep Evolution Field (MOSDEF) survey to place the first constraints on the nebular attenuation curve at high redshift. Based on the first five low-order Balmer emission lines detected in the composite spectra of these galaxies (Hα through ), we derive a nebular attenuation curve that is similar in shape to that of the Galactic extinction curve, suggesting that the dust covering fraction and absorption/scattering properties along the lines of sight to massive stars at high redshift are similar to those of the average Milky Way sight line. The curve derived here implies nebular reddening values that are, on average, systematically larger than those derived for the stellar continuum. In the context of stellar population synthesis models that include the effects of stellar multiplicity, the difference in reddening of the nebular lines and stellar continuum may imply molecular cloud crossing timescales that are a factor of longer than those inferred for local molecular clouds, star formation rates that are constant or increasing with time such that newly formed and dustier OB associations always dominate the ionizing flux, and/or that the dust responsible for reddening the nebular emission may be associated with nonmolecular (i.e., ionized and neutral) phases of the interstellar medium. Our analysis points to a variety of investigations of the nebular attenuation curve that will be enabled with the next generation of ground- and space-based facilities.

Evidence of Cosmic-Ray Excess from Local Giant Molecular Clouds

Abstract We report the analysis of the Fermi-Large Area Telescope data from six nearby giant molecular clouds (MCs) belonging to the Gould Belt and the Aquila Rift regions. The high statistical γ-ray spectra above 3 GeV well described by power laws make it possible to derive precise estimates of the cosmic-ray (CR) distribution in the MCs. The comparison of γ-ray spectra of Taurus, Orion A, and Orion B clouds with the model expected from Alpha Magnetic Spectrometer (AMS-02) CR measurements confirms these clouds as passive clouds, immersed in an AMS-02-like CR spectrum. A similar comparison of Aquila Rift, Rho Oph, and Cepheus spectra yields significant deviation in both spectral indices and absolute fluxes, which can imply an additional acceleration of CRs throughout the entire clouds. Besides, the theoretical modeling of the excess γ-ray spectrum of these clouds, assuming π 0-decay interaction of CRs in the cloud, gives a considerable amount of an enhanced CR energy density and it shows a significant deviation in spectral shapes compared to the average AMS-02 CR spectrum between 30 GeV and 10 TeV. We suggest that this variation in the CR spectrum of Cepheus could be accounted for by an efficient acceleration in the shocks of winds of OB associations, while in Rho Oph, similar acceleration can be provided by multiple T-Tauri stars populated in the whole cloud. In the case of Aquila Rift, the excess in absolute CR flux could be related to an additional acceleration of CRs by supernova remnants or propagation effects in the cloud.

Distances and Statistics of Local Molecular Clouds in the First Galactic Quadrant

Abstract We present an analysis of local molecular clouds ( km s−1, i.e., <1.5 kpc) in the first Galactic quadrant ( and ), a pilot region of the Milky Way Imaging Scroll Painting (MWISP) CO survey. Using the Spectral Clustering for Interstellar Molecular Emission Segmentation algorithm to divide large molecular clouds into moderate-size ones, we determined distances to 28 molecular clouds with the background-eliminated extinction-parallax method using the Gaia DR2 parallax measurements aided by A G and A V ; the distance ranges from 250 to about 1.5 kpc. These incomplete distance samples indicate a linear relationship between the distance and the radial velocity (V LSR) with a scatter of 0.16 kpc, and kinematic distances may be systematically larger for local molecular clouds. In order to investigate fundamental properties of molecular clouds, such as the total sample number, the line width, the brightness temperature, the physical area, and the mass, we decompose the spectral cube using the density-based spatial clustering of applications with noise (DBSCAN) algorithm. Post-selection criteria are imposed on DBSCAN clusters to remove the noise contamination, and we found that the separation of molecular cloud individuals is reliable based on a definition of independent consecutive structures in l–b–V space. The completeness of the local molecular cloud flux collected by the MWISP CO survey is about 80%. The physical area, A, shows a power-law distribution, dN/d , while the molecular cloud mass also follows a power-law distribution but is slightly flatter, dN/dM ∝ M −1.96 ± 0.11.

First Resolved Dust Continuum Measurements of Individual Giant Molecular Clouds in the Andromeda Galaxy

Abstract In our local Galactic neighborhood, molecular clouds are best studied using a combination of dust measurements, to determine robust masses, sizes, and internal structures of the clouds, and molecular-line observations to determine cloud kinematics and chemistry. We present here the first results of a program designed to extend such studies to nearby galaxies beyond the Magellanic Clouds. Utilizing the wideband upgrade of the Submillimeter Array (SMA) at 230 GHz, we have obtained the first continuum detections of the thermal dust emission on sub-GMC scales (∼15 pc) within the Andromeda galaxy (M31). These include the first resolved continuum detections of dust emission from individual giant molecular clouds (GMCs) beyond the Magellanic Clouds. Utilizing a powerful capability of the SMA, we simultaneously recorded CO(2−1) emission with identical (u, v) coverage, astrometry, and calibration, enabling the first measurements of the CO conversion factor, α CO(2−1), toward individual GMCs across an external galaxy. Our direct measurement yields an average CO-to-dust mass conversion factor of M ⊙ (K km s−1 pc2)−1 for the J = 2−1 transition. This value does not appear to vary with galactocentric radius. Assuming a constant gas-to-dust ratio of 136, the resulting α CO = 5.7 ± 2.4 M ⊙ (K km s−1 pc2)−1 for the 2−1 transition is in excellent agreement with that of GMCs in the Milky Way, given the uncertainties. Finally, using the same analysis techniques, we compare our results with observations of the local Orion molecular clouds, placed at the distance of M31 and simulated to appear as they would if observed by the SMA.

A compendium of distances to molecular clouds in the Star Formation Handbook

Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method that combines stellar photometric data withGaiaDR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical accuracy of ∼5%. After refining the technique to target lower latitudes and incorporating deep optical data from DECam in the southern Galactic plane, we have derived a catalog of distances to molecular clouds in Reipurth (2008, Star Formation Handbook, Vols. I and II) which contains a large fraction of the molecular material in the solar neighborhood. Comparison with distances derived from maser parallax measurements towards the same clouds shows our method produces consistent distances with ≲10% scatter for clouds across our entire distance spectrum (150 pc−2.5 kpc). We hope this catalog of homogeneous distances will serve as a baseline for future work.

A Large Catalog of Accurate Distances to Local Molecular Clouds: The Gaia DR2 Edition

Abstract We present a uniform catalog of accurate distances to local molecular clouds informed by the Gaia DR2 data release. Our methodology builds on that of Schlafly et al. First, we infer the distance and extinction to stars along sightlines toward the clouds using optical and near-infrared photometry. When available, we incorporate knowledge of the stellar distances obtained from Gaia DR2 parallax measurements. We model these per-star distance–extinction estimates as being caused by a dust screen with a 2D morphology derived from Planck at an unknown distance, which we then fit for using a nested sampling algorithm. We provide updated distances to the Schlafly et al. sightlines toward the Dame et al. and Magnani et al. clouds, finding good agreement with the earlier work. For a subset of 27 clouds, we construct interactive pixelated distance maps to further study detailed cloud structure, and find several clouds which display clear distance gradients and/or are comprised of multiple components. We use these maps to determine robust average distances to these clouds. The characteristic combined uncertainty on our distances is ≈5%–6%, though this can be higher for clouds at greater distances, due to the limitations of our single-cloud model.

Gas and Dust Properties in the Chamaeleon Molecular Cloud Complex Based on the Optically Thick H i

Abstract Gas and dust properties in the Chamaeleon molecular cloud complex have been investigated with emission lines from the atomic hydrogen (H i) and 12CO molecules, dust optical depth at 353 GHz (τ 353), and J-band infrared extinction (A J ). We have found a scatter correlation between the H i integrated intensity (W H i ) and τ 353 in the Chamaeleon region. The scattering has been examined in terms of a possible large optical depth in H i emission (τ H i ) using a total column density (N H) model based on τ 353. A nonlinear relation of τ 353 with the ∼1.2 power of A J has been found in opaque regions (A J ≳ 0.3 mag), which may indicate a dust evolution effect. If we apply this nonlinear relation to the N H model (i.e., N H ∝ τ 353 1 / 1.2 ) allowing arbitrary τ H i , the model curve reproduces well the W H i –τ 353 scatter correlation, suggesting optically thick H i (τ H i ∼ 1.3) extended around the molecular clouds. Based on the correlations between the CO integrated intensity and the N H model, we have derived the CO-to-H2 conversion factor (X CO) on ∼1.°5 scales (corresponding to ∼4 pc) and found spatial variations of X CO ∼ (0.5–3) × 1020 cm−2 K−1 km−1 s across the cloud complex, possibly depending on the radiation field inside or surrounding the molecular clouds. These gas properties found in the Chamaeleon region are discussed through a comparison with other local molecular cloud complexes.

Gamma-Ray Emission from Molecular Clouds Generated by Penetrating Cosmic Rays

Abstract We analyze the processes governing cosmic-ray (CR) penetration into molecular clouds and the resulting generation of gamma-ray emission. The density of CRs inside a cloud is depleted at lower energies due to the self-excited MHD turbulence. The depletion depends on the effective gas column density (“size”) of the cloud. We consider two different environments where the depletion effect is expected to be observed. For the central molecular zone, the expected range of CR energy depletion is E ≲ 10 GeV, leading to the depletion of gamma-ray flux below E γ ≈ 2 GeV. This effect can be important for the interpretation of the GeV gamma-ray excess in the Galactic Center, which has been revealed from the standard model of CR propagation (assuming the CR spectrum inside a cloud to be equal to the interstellar spectrum). Furthermore, recent observations of some local molecular clouds suggest the depletion of the gamma-ray emission, indicating possible self-modulation of the penetrating low-energy CRs.

The Star Formation Rate in the Gravoturbulent Interstellar Medium

Abstract Stars form in supersonic turbulent molecular clouds that are self-gravitating. We present an analytic determination of the star formation rate (SFR) in a gravoturbulent medium based on the density probability distribution function of molecular clouds having a piecewise lognormal (LN) and power-law (PL) form. This is in contrast to previous analytic SFR models that are governed primarily by interstellar turbulence, which sets purely LN density probability distribution functions (PDFs). In the gravoturbulent SFR model described herein, low-density gas resides in the LN portion of the PDF. Gas becomes gravitationally unstable past a critical density (ρ crit), and the PDF begins to form a PL. As the collapse of the cloud proceeds, the transitional density (ρ t ) between the LN and PL portions of the PDF moves toward lower density while the slope of the PL (α) becomes increasingly shallow. The SFR per free-fall time is calculated via an integral over the LN from ρ crit to ρ t and an integral over the PL from ρ t to the maximum density. As α becomes shallower, the SFR accelerates beyond the expected values calculated from an LN density PDF. We show that the star formation efficiency per free-fall time in observations of local molecular clouds increases with shallower PDF PL slopes, in agreement with our model. Our model can explain why star formation is spatially and temporally variable within a cloud and why the depletion times observed in local and extragalactic giant molecular clouds vary. Both star-bursting and quiescent star-forming systems can be explained without the need to invoke extreme variations of turbulence in the local interstellar environment.

Three-dimensional dust mapping in the Orion complex, combining Gaia-TGAS, 2MASS, and WISE

We present a map of the three-dimensional (3D) distribution of dust in the Orion complex. Orion is the closest site of high-mass star formation, making it an excellent laboratory for studying the interstellar medium and star formation. We used data from the Gaia-TGAS catalogue combined with photometry from 2MASS and WISE to get the distances and extinctions of individual stars in the vicinity of the Orion complex. We use a Gaussian process and adopt a non-parametric method to infer the probability distribution function of the dust densities at arbitrary points throughout the region. We map the dust distribution towards different parts of the Orion complex. We find that the distance and depth of the cloud are compatible with other recent works, which show that the method can be applicable to local molecular clouds to map their 3D dust distribution. We also demonstrate the danger of only using colours of stars to derive their extinctions without considering further physical constraints, such as the colour-magnitude diagram (CMD).

The accelerating pace of star formation

We study the temporal and spatial distribution of star formation rates in four well-studied star-forming regions in local molecular clouds(MCs): Taurus, Perseus, $\rho$ Ophiuchi, and Orion A. Using published mass and age estimates for young stellar objects in each system, we show that the rate of star formation over the last 10 Myrs has been accelerating and is (roughly) consistent with a $t^2$ power law. This is in line with previous studies of the star formation history of molecular clouds and with recent theoretical studies. We further study the clustering of star formation in the Orion Nebula Cluster(ONC). We examine the distribution of young stellar objects as a function of their age by computing an effective half-light radius for these young stars subdivided into age bins. We show that the distribution of young stellar objects is broadly consistent with the star formation being entirely localized within the central region. We also find a slow radial expansion of the newly formed stars at a velocity of $v=0.17\,{\rm km\,s}^{-1}$, which is roughly the sound speed of the cold molecular gas. This strongly suggests the dense structures that form stars persist much longer than the local dynamical time. We argue that this structure is quasi-static in nature and is likely the result of the density profile approaching an attractor solution as suggested by recent analytic and numerical analysis.

The parsec–scale relationship between ICO and AV in local molecular clouds

We measure the parsec-scale relationship between integrated CO intensity (I_CO) and visual extinction (A_V) in 24 local molecular clouds using maps of CO emission and dust optical depth from Planck. This relationship informs our understanding of CO emission across environments, but clean Milky Way measurements remain scarce. We find uniform I_CO for a given A_V, with the results bracketed by previous studies of the Pipe and Perseus clouds. Our measured I_CO-A_V relation broadly agrees with the standard Galactic CO-to-H2 conversion factor, the relation found for the Magellanic clouds at coarser resolution, and numerical simulations by Glover & Clark (2016). This supports the idea that CO emission primarily depends on shielding, which protects molecules from dissociating radiation. Evidence for CO saturation at high A_V and a threshold for CO emission at low A_V varies remains uncertain due to insufficient resolution and ambiguities in background subtraction. Resolution of order 0.1 pc may be required to measure these features. We use this I_CO-AV relation to predict how the CO-to-H2 conversion factor (X_CO) would change if the Solar Neighborhood clouds had different dust-to-gas ratio (metallicity). The calculations highlight the need for improved observations of the CO emission threshold and HI shielding layer depth. They are also sensitive to the shape of the column density distribution. Because local clouds collectively show a self-similar distribution, we predict a shallow metallicity dependence for X_CO down to a few tenths of solar metallicity. However, our calculations also imply dramatic variations in cloud-to-cloud X_CO at subsolar metallicity.

Multiple regimes and coalescence timescales for massive black hole pairs; the critical role of galaxy formation physics

We discuss the latest results of numerical simulations following the orbital decay of massive black hole pairs in galaxy mergers. We highlight important differences between gas-poor and gas-rich hosts, and between orbital evolution taking place at high redshift as opposed to low redshift. Two effects have a huge impact and are rather novel in the context of massive black hole binaries. The first is the increase in characteristic density of galactic nuclei of merger remnants as galaxies are more compact at high redshift due to the way dark halo collapse depends on redshift. This leads naturally to hardening timescales due to 3-body encounters that should decrease by two orders of magnitude up to z = 4. It explains naturally the short binary coalescence timescale, ∼ 10 Myr, found in novel cosmological simulations that follow binary evolution from galactic to milliparsec scales. The second one is the inhomogeneity of the interstellar medium in massive gas-rich disks at high redshift. In the latter star forming clumps 1-2 orders of magnitude more massive than local Giant Molecular Clouds (GMCs) can scatter massive black holes out of the disk plane via gravitational perturbations and direct encounters. This renders the character of orbital decay inherently stochastic, often increasing orbital decay timescales by as much as a Gyr. At low redshift a similar regime is present at scales of 1 − 10 pc inside Circumnuclear Gas Disks (CNDs). In CNDs only massive black holes with masses below 107M⊙ can be significantly perturbed. They decay to sub-pc separations in up to ∼ 108 yr rather than the in just a few million years as in a smooth CND. Finally implications for building robust forecasts of LISA event rates are discussed.

NGC 1980 Is Not a Foreground Population of Orion: Spectroscopic Survey of Young Stars with Low Extinction in Orion A

Abstract We perform a spectroscopic survey of the foreground population in Orion A with MMT/Hectospec. We use these data, along with archival spectroscopic data and photometric data, to derive spectral types, extinction values, and masses for 691 stars. Using the Spitzer Space Telescope data, we characterize the disk properties of these sources. We identify 37 new transition disk (TD) objects, 1 globally depleted disk candidate, and 7 probable young debris disks. We discover an object with a mass of less than 0.018–0.030 M ⊙, which harbors a flaring disk. Using the Hα emission line, we characterize the accretion activity of the sources with disks, and confirm that the fraction of accreting TDs is lower than that of optically thick disks (46% ± 7% versus 73% ± 9%, respectively). Using kinematic data from the Sloan Digital Sky Survey and APOGEE INfrared Spectroscopy of the Young Nebulous Clusters program (IN-SYNC), we confirm that the foreground population shows similar kinematics to their local molecular clouds and other young stars in the same regions. Using the isochronal ages, we find that the foreground population has a median age of around 1–2 Myr, which is similar to that of other young stars in Orion A. Therefore, our results argue against the presence of a large and old foreground cluster in front of Orion A.

Mapping the Extinction Curve in 3D: Structure on Kiloparsec Scales

Abstract Near-infrared spectroscopy from APOGEE and wide-field optical photometry from Pan-STARRS1 have recently made precise measurements of the shape of the extinction curve possible for tens of thousands of stars, parameterized by R(V). These measurements revealed structures in R(V) with large angular scales, which are challenging to explain in existing dust paradigms. In this work, we combine three-dimensional maps of dust column density with R(V) measurements to constrain the three-dimensional distribution of R(V) in the Milky Way. We find that the variations in R(V) are correlated on kiloparsec scales. In particular, most of the dust within one kiloparsec in the outer Galaxy, including many local molecular clouds (Orion, Taurus, Perseus, California, and Cepheus), has a significantly lower R(V) than more distant dust in the Milky Way. These results provide new input to models of dust evolution and processing, and complicate the application of locally derived extinction curves to more distant regions of the Milky Way and to other galaxies.