Showing papers in "arXiv: Astrophysics of Galaxies in 2013"

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

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

The Radial Velocity Experiment (RAVE): Fifth Data Release

Abstract: Data Release 5 (DR5) of the Radial Velocity Experiment (RAVE) is the fifth data release from a magnitude-limited (9< I < 12) survey of stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra ($R\sim7500$) covering the Ca-triplet region (8410-8795A) span the complete time frame from the start of RAVE observations in 2003 to their completion in 2013. Radial velocities from 520,781 spectra of 457,588 unique stars are presented, of which 255,922 stellar observations have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. For our main DR5 catalog, stellar parameters (effective temperature, surface gravity, and overall metallicity) are computed using the RAVE DR4 stellar pipeline, but calibrated using recent K2 Campaign 1 seismic gravities and Gaia benchmark stars, as well as results obtained from high-resolution studies. Also included are temperatures from the Infrared Flux Method, and we provide a catalogue of red giant stars in the dereddened color $(J-Ks)_0$ interval (0.50,0.85) for which the gravities were calibrated based only on seismology. Further data products for sub-samples of the RAVE stars include individual abundances for Mg, Al, Si, Ca, Ti, Fe, and Ni, and distances found using isochrones. Each RAVE spectrum is complemented by an error spectrum, which has been used to determine uncertainties on the parameters. The data can be accessed via the RAVE Web site or the Vizier database.

The age structure of stellar populations in the solar vicinity. Clues of a two-phase formation history of the Milky Way disk

Abstract: We analyze high quality abundances data of solar neighborhood stars and show that there are two distinct regimes of [alpha/Fe] versus age which we identify as the epochs of the thick and thin disk formation. A tight correlation between metallicity and [alpha/Fe] versus age is clearly identifiable on thick disk stars, implying that this population formed from a well mixed ISM, over a time scale of 4-5 Gyr. Thick disk stars vertical velocity dispersion correlate with age, with the youngest objects having as small scale heights as those of thin disk stars. A natural consequence of these two results is that a vertical metallicity gradient is expected in this population. We suggest that the thick disk set the initial conditions for the formation of the inner thin disk. This provides also an explanation of the apparent coincidence between the step in metallicity at 7-10 kpc in the thin disk and the confinment of the thick disk at about R<10 kpc. We suggest that the outer thin disk developped outside the influence of the thick disk, but also that the high alpha-enrichment of the outer regions may originate from a primordial pollution by the gas expelled from the thick disk. Local metal-poor thin disk stars, whose properties are best explained by an origin in the outer disk, are shown to be as old as the youngest thick disk (9-10 Gyr), implying that the outer thin disk started to form while the thick disk formation was still on-going in the inner Galaxy. We point out that, given the tight age-abundance relations in the thick disk, an inside-out process would give rise to a radial gradient in abundances in this population which is not observed. Finally, we argue that the data discussed here leave little room for radial migration, either to have contaminated the solar vicinity, or, to have redistributed stars in significant proportion across the solar annulus.

Final model independent result of DAMA/LIBRA-phase1

Abstract: The results obtained with the total exposure of 1.04 ton x yr collected by DAMA/LIBRA-phase1 deep underground at the Gran Sasso National Laboratory (LNGS) of the I.N.F.N. during 7 annual cycles (i.e. adding a further 0.17 ton x yr exposure) are presented. The DAMA/LIBRA-phase1 data give evidence for the presence of Dark Matter (DM) particles in the galactic halo, on the basis of the exploited model independent DM annual modulation signature by using highly radio-pure NaI(Tl) target, at 7.5 sigma C.L.. Including also the first generation DAMA/NaI experiment (cumulative exposure 1.33 ton x yr, corresponding to 14 annual cycles), the C.L. is 9.3 sigma and the modulation amplitude of the single-hit events in the (2-6) keV energy interval is: (0.0112 \pm 0.0012) cpd/kg/keV; the measured phase is (144 \pm 7) days and the measured period is (0.998 \pm 0.002) yr, values well in agreement with those expected for DM particles. No systematic or side reaction able to mimic the exploited DM signature has been found or suggested by anyone over more than a decade.

The 2013 Release of Cloudy

Abstract: This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.

The WISE Catalog of Galactic HII Regions

Abstract: Using data from the all-sky Wide-Field Infrared Survey Explorer (WISE) satellite, we made a catalog of over 8000 Galactic HII regions and HII region candidates by searching for their characteristic mid-infrared (MIR) morphology. WISE has sufficient sensitivity to detect the MIR emission from HII regions located anywhere in the Galactic disk. We believe this is the most complete catalog yet of regions forming massive stars in the Milky Way. Of the ~8000 cataloged sources, ~1500 have measured radio recombination line (RRL) or H$\alpha$ emission, and are thus known to be HII regions. This sample improves on previous efforts by resolving HII region complexes into multiple sources and by removing duplicate entries. There are ~2500 candidate HII regions in the catalog that are spatially coincident with radio continuum emission. Our group's previous RRL studies show that ~95% of such targets are HII regions. We find that ~500 of these candidates are also positionally associated with known HII region complexes, so the probability of their being bona fide HII regions is even higher. At the sensitivity limits of existing surveys, ~4000 catalog sources show no radio continuum emission. Using data from the literature, we find distances for ~1500 catalog sources, and molecular velocities for ~1500 HII region candidates.

The Ages of 55 Globular Clusters as Determined Using an Improved delta V_TO^HB Method Along with Color-Magnitude Diagram Constraints, and Their Implications for Broader Issues

Abstract: Ages have been derived for 55 globular clusters (GCs) from overlays of isochrones onto the turnoff photometry, assuming distances based on fits of zero-age horizontal branch (ZAHB) models to the lower bound of the observed distributions of HB stars. The error bar arising just from the "fitting" of ZAHBs and isochrones is ~ +/- 0.25 Gyr, while that associated with distance and chemical abundance uncertainties is ~ +/- 1.5-2 Gyr. Ages vary from mean values of ~12.5 Gyr at [Fe/H] -1.0. At intermediate metallicities, the age-metallicity relation (AMR) appears to be bifurcated: one branch apparently contains clusters with disk-like kinematics, whereas the other branch is populated by clusters with halo-type orbits. There is no apparent dependence of age on Galactocentric distance (R_G) nor is there a clear correlation of HB type with age. Subtle variations in the subgiant branch (SGB) slopes of [Fe/H] < -1.5 GCs are tentatively attributed to helium abundance differences. Curiously, GCs with steep "M13-like" SGBs tend to be massive systems, located at small R_G, that show the strongest evidence for multiple stellar populations. The others are typically low-mass systems that, at the present time, should not be able to retain the matter lost by mass-losing stars. The apparent separation of the two groups in terms of their present-day gas retention properties is difficult to understand if all GCs were initial ~20 times their current masses. The lowest mass systems may have never been able to retain enough gas to produce a significant population of second-generation stars; in this case, the observed light element abundance variations were presumably present in the gas out of which the observed cluster stars formed.

A Herschel [CII] Galactic plane survey I: the global distribution of ISM gas components

Abstract: [Abridged] The [CII] 158um line is an important tool for understanding the life cycle of interstellar matter. Ionized carbon is present in a variety of phases of the interstellar medium, including the diffuse ionized medium, warm and cold atomic clouds, clouds in transition from atomic to molecular, and dense and warm photon dominated regions (PDRs). The Galactic Observations of Terahertz C+ (GOTC+) project surveys the [CII] line over the entire Galactic disk with velocity-resolved observations using the Herschel/HIFI instrument. We present the first longitude-velocity maps of the [CII] emission for Galactic latitudes b=0deg, +-0.5deg, and +-1.0deg. [CII] emission is mostly associated with spiral arms, mainly emerging from Galactocentric distances between 4 and 10 kpc. We estimate that most of the observed [CII] emission is produced by dense PDRs (47%), with smaller contributions from CO-dark H2 gas (28%), cold atomic gas (21%), and ionized gas (4%). Atomic gas inside the Solar radius is mostly in the form of cold neutral medium (CNM), while the warm neutral medium (WNM) gas dominates the outer galaxy. The average fraction of CNM relative to total atomic gas is 43%. We find that the warm and diffuse CO-dark H2 is distributed over a larger range of Galactocentric distances (4-11 kpc) than the cold and dense H2 gas traced by 12CO and 13CO (4-8kpc). The fraction of CO-dark H2 to total H2 increases with Galactocentric distance, ranging from 20% at 4 kpc to 80% at 10 kpc. On average, CO-dark H2 accounts for 30% of the molecular mass of the Milky Way. When the CO-dark H2 component is included, the radial distribution of the CO-to-H2 conversion factor is steeper than that when only molecular gas traced by CO is considered. Most of the observed [CII] emission emerging from dense PDRs is associated with modest far-ultraviolet fields in the range chi0~1-30.

Gaia FGK Benchmark Stars - Metallicity

Abstract: To calibrate automatic pipelines that determine atmospheric parameters of stars, one needs a sample of stars -- ``benchmark stars'' -- with well defined parameters to be used as a reference We provide a detailed documentation of the determination of the iron abundance of the 34 FGK-type benchmark stars selected to be the pillars for calibration of the one billion Gaia stars. They cover a wide range of temperatures, surface gravities and metallicities. Up to seven different methods were used to analyze an observed spectral library of high resolution and high signal-to-noise ratio. The metallicity was determined assuming a value of effective temperature and surface gravity obtained from fundamental relations, i.e. these parameters were known a priori independently from the spectra. We present a set of metallicity values obtained in a homogeneous way for our sample of Benchmark Stars. In addition to this value, we provide a detailed documentation of the associated uncertainties. Finally, we report for the first time a value of the metallicity of the cool giant psi Phe.

The nature of LINER galaxies: Ubiquitous hot old stars and rare accreting black holes

Abstract: Galaxies, which often contain ionised gas, sometimes also exhibit a so-called low-ionisation nuclear emission line region (LINER). For 30 years this was attributed to a central mass-accreting supermassive black hole (AGN) of low luminosity, making LINER galaxies the largest AGN-sub-population, dominating in numbers over higher luminosity Seyfert galaxies and quasars. This, however, poses a serious problem. While the inferred energy balance is plausible, many LINERs clearly do not contain any other independent signatures of an AGN. Using integral field spectroscopic data from the CALIFA survey, we aim at comparing the observed radial surface brightness profiles with what is expected from illumination by an AGN. Essential for this analysis is a proper extraction of emission-lines, especially weak lines such as the Balmer Hb line which is superposed on an absorption trough. To accomplish this, we use the GANDALF code which simultaneously fits the underlying stellar continuum and emission lines. We show for 48 galaxies with LINER-like emission, that the radial emission-line surface brightness profiles are inconsistent with ionisation by a central point-source and hence cannot be due to an AGN alone. The most probable explanation for the excess LINER-like emission is ionisation by evolved stars during the short but very hot and energetic phase known as post-AGB. This leads us to an entirely new interpretation. Post-AGB stars are ubiquitous and their ionising effect should be potentially observable in every galaxy with gas present and stars older than ~1 Gyr, unless a stronger radiation field from young hot stars or an AGN outshines them. This means that galaxies with LINER-like emission are in fact not a class defined by a property, but rather by the absence of a property. It also explains why LINER emission is observed mostly in massive galaxies with old stars and little star formation.

A three-phase chemical model of hot cores: the formation of glycine

Abstract: A new chemical model is presented that simulates fully-coupled gas-phase, grain-surface and bulk-ice chemistry in hot cores. Glycine (NH2CH2COOH), the simplest amino acid, and related molecules such as glycinal, propionic acid and propanal, are included in the chemical network. Glycine is found to form in moderate abundance within and upon dust-grain ices via three radical-addition mechanisms, with no single mechanism strongly dominant. Glycine production in the ice occurs over temperatures ~40-120 K. Peak gas-phase glycine fractional abundances lie in the range 8 x 10^{-11} - 8 x 10^{-9}, occuring at ~200 K, the evaporation temperature of glycine. A gas-phase mechanism for glycine production is tested and found insignificant, even under optimal conditions. A new spectroscopic radiative-transfer model is used, allowing the translation and comparison of the chemical-model results with observations of specific sources. Comparison with the nearby hot-core source NGC 6334 IRS1 shows excellent agreement with integrated line intensities of observed species, including methyl formate. The results for glycine are consistent with the current lack of a detection of this molecule toward other sources; the high evaporation temperature of glycine renders the emission region extremely compact. Glycine detection with ALMA is predicted to be highly plausible, for bright, nearby sources with narrow emission lines. Photodissociation of water and subsequent hydrogen-abstraction from organic molecules by OH, and NH2, are crucial to the build-up of complex organic species in the ice. The inclusion of alternative branches within the network of radical-addition reactions appears important to the abundances of hot-core molecules; less favorable branching ratios may remedy the anomalously high abundance of glycolaldehyde predicted by this and previous models.

The Dark Matter Halo of the Milky Way, AD 2013

Abstract: We derive the mass model of the Milky Way (MW) using a cored dark matter (DM) halo profile and recent data. The method used consists in fitting a spherically symmetric model of the Galaxy with a Burkert DM halo profile to available data: MW terminal velocities in the region inside the solar circle, circular velocity as recently estimated from maser star forming regions at intermediate radii, and velocity dispersions of stellar halo tracers for the outermost Galactic region. The latter are reproduced by integrating the Jeans equation for every modeled mass distribution, and by allowing for different velocity anisotropies for different tracer populations. For comparison we also consider a Navarro-Frenk-White profile. We find that the cored profile is the preferred one, with a shallow central density of rho_H~4x10^7M_s/kpc^3 and a large core radius R_H~10 kpc, as observed in external spirals and in agreement with the mass model underlying the Universal Rotation Curve of spirals. We describe also the derived model uncertainties, which are crucially driven by the poorly constrained velocity dispersion anisotropies of halo tracers. The emerging cored DM distribution has implications for the DM annihilation angular profile, which is much less boosted in the Galactic center direction with respect to the case of the standard \Lambda CDM, NFW profile. Using the derived uncertainties we discuss finally the limitations and prospects to discriminate between cored and cusped DM profile with a possible observed diffuse DM annihilation signal. The present mass model aims to characterize the present-day description of the distribution of matter in our Galaxy, which is needed to frame current crucial issues of Cosmology, Astrophysics and Elementary Particles.

Discovery of a binary brown dwarf at 2 parsecs from the sun 1

Abstract: I am using multi-epoch astrometry from the Wide-field Infrared Survey Explorer (WISE) to search for new members of the solar neighborhood via their high proper motions. Through this work, I have identified WISE J104915.57-531906.1 as a high proper motion object and have found additional detections in images from the Digitized Sky Survey, the Two Micron All-Sky Survey, and the Deep Near-Infrared Survey of the Southern Sky. I have measured a parallax of 0.496±0.037 ′′ (2.0±0.15 pc) from the astrometry in these surveys, making WISE J104915.57-531906.1 the third closest system to the Sun. During spectroscopic observations with GMOS at Gemini Observatory, an i-band acquisition image resolved it as a 1. ′′ 5 (3 AU) binary. A spectrum was collected for the primary, which I classify as L8±1. The secondary is probably near the L/T transition as well given that it is only modestly fainter than the primary (�i = 0.45 mag). Subject headings: brown dwarfs — infrared: stars — proper motions — solar neighborhood — stars: low-mass

The large-scale structure of the halo of the Andromeda Galaxy Part I: global stellar density, morphology and metallicity properties

Abstract: We present an analysis of the large-scale structure of the halo of the Andromeda galaxy, based on the Pan-Andromeda Archeological Survey (PAndAS), currently the most complete map of resolved stellar populations in any galactic halo. Despite copious substructure, the global halo populations follow closely power law profiles that become steeper with increasing metallicity. We divide the sample into stream-like populations and a smooth halo component. Fitting a three-dimensional halo model reveals that the most metal-poor populations ([Fe/H] -0.6). The space density of the smooth metal-poor component has a global power-law slope of -3.08+/-0.07, and a non-parametric fit shows that the slope remains nearly constant from 30kpc to 300kpc. The total stellar mass in the halo at distances beyond 2 degrees is 1.1x10^10 Solar masses, while that of the smooth component is 3x10^9 Solar masses. Extrapolating into the inner galaxy, the total stellar mass of the smooth halo is plausibly 8x10^9 Solar masses. We detect a substantial metallicity gradient, which declines from [Fe/H]=-0.7 at R=30kpc to [Fe/H]=-1.5 at R=150kpc for the full sample, with the smooth halo being 0.2dex more metal poor than the full sample at each radius. While qualitatively in-line with expectations from cosmological simulations, these observations are of great importance as they provide a prototype template that such simulations must now be able to reproduce in quantitative detail.

Molecular Clouds Toward the Super Star Cluster NGC3603; Possible Evidence for a Cloud-Cloud Collision in Triggering the Cluster Formation

Abstract: We present new large field observations of molecular clouds with NANTEN2 toward the super star cluster NGC3603 in the transitions 12CO(J=2-1, J=1-0) and 13CO(J=2-1, J=1-0) We suggest that two molecular clouds at 13 km s-1 and 28 km s-1 are associated with NGC3603 as evidenced by higher temperatures toward the H II region as well as morphological correspondence The mass of the clouds is too small to gravitationally bind them, given their relative motion of ~20 km s-1 We suggest that the two clouds collided with each other a Myr ago to trigger the formation of the super star cluster This scenario is able to explain the origin of the highest mass stellar population in the cluster which is as young as a Myr and is segregated within the central sub-pc of the cluster This is the second super star cluster along side Westerlund2 where formation may have been triggered by a cloud-cloud collision

A young protoplanet candidate embedded in the circumstellar disk of HD100546

Abstract: We present high-contrast observations of the circumstellar environment of the Herbig Ae/Be star HD100546. The final 3.8 micron image reveals an emission source at a projected separation of 0.48"+-0.04" (corresponding to ~47+-4 AU at a position angle of 8.9+-0.9 degree. The emission appears slightly extended with a point source component with an apparent magnitude of 13.2+-0.4 mag. The position of the source coincides with a local deficit in polarization fraction in near-infrared polarimetric imaging data, which probes the surface of the well-studied circumstellar disk of HD100546. This suggests a possible physical link between the emission source and the disk. Assuming a disk inclination of ~47 degree the de-projected separation of the object is ~68 AU. Assessing the likelihood of various scenarios we favor an interpretation of the available high-contrast data with a planet in the process of forming. Follow-up observations in the coming years can easily distinguish between the different possible scenarios empirically. If confirmed, HD100546 "b" would be a unique laboratory to study the formation process of a new planetary system, with one giant planet currently forming in the disk and a second planet possibly orbiting in the disk gap at smaller separations.

The AGORA High-Resolution Galaxy Simulations Comparison Project

Abstract: We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

An Overview of the Dwarf Galaxy Survey

Abstract: The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

Interstellar water chemistry: from laboratory to observations

Abstract: Water is observed throughout the universe, from diffuse interstellar clouds to protoplanetary disks around young stars, and from comets in our own solar system and exoplanetary atmospheres to galaxies at high redshifts This review summarizes the spectroscopy and excitation of water in interstellar space as well as the basic chemical processes that form and destroy water under interstellar conditions Three major routes to water formation are identified: low temperature ion-molecule chemistry, high-temperature neutral-neutral chemistry and gas-ice chemistry The rate coefficients of several important processes entering the networks are discussed in detail; several of them have been determined only in the last decade through laboratory experiments and theoretical calculations Astronomical examples of each of the different chemical routes are presented using data from powerful new telescopes, in particular the Herschel Space Observatory Basic chemical physics studies remain critically important to analyze astronomical data

Planck 2013 results. XIII. Galactic CO emission

Abstract: Rotational transition lines of CO play a major role in molecular radio astronomy and in particular in the study of star formation and the Galactic structure. Although a wealth of data exists in the Galactic plane and some well-known molecular clouds, there is no available CO high sensitivity all-sky survey to date. Such all-sky surveys can be constructed using the \Planck\ HFI data because the three lowest CO rotational transition lines at 115, 230 and 345 GHz significantly contribute to the signal of the 100, 217 and 353 GHz HFI channels respectively. Two different component separation methods are used to extract the CO maps from Planck HFI data. The maps obtained are then compared to one another and to existing external CO surveys. From these quality checks the best CO maps in terms of signal to noise and/or residual foreground contamination are selected. Three sets of velocity-integrated CO emission maps are produced: Type 1 maps of the CO (1-0), (2-1), and (3-2) rotational transitions with low foreground contamination but moderate signal-to-noise ratio; Type 2 maps for the (1-0) and (2-1) transitions with a better signal-to-noise ratio; and one Type 3 map, a line composite map with the best signal-to-noise ratio in order to locate the faintest molecular regions. The maps are described in detail. They are shown to be fully compatible with previous surveys of parts of the Galactic Plane and also of fainter regions out of the Galactic plane. The Planck HFI velocity-integrated CO maps for the (1-0), (2-1), and (3-2) rotational transitions provide an unprecedented all-sky CO view of the Galaxy. These maps are also of great interest to monitor potential CO contamination on CMB \Planck\ studies.

The CO-to-H2 Conversion Factor

Abstract: CO line emission represents the most accessible and widely used tracer of the molecular interstellar medium. This renders the translation of observed CO intensity into total H2 gas mass critical to understand star formation and the interstellar medium in our Galaxy and beyond. We review the theoretical underpinning, techniques, and results of efforts to estimate this CO-to-H2 "conversion factor," Xco, in different environments. In the Milky Way disk, we recommend a conversion factor Xco = 2x10^{20} cm^-2/(K km/s)^-1 with +/-30% uncertainty. Studies of other "normal galaxies" return similar values in Milky Way-like disks, but with greater scatter and systematic uncertainty. Departures from this Galactic conversion factor are both observed and expected. Dust-based determinations, theoretical arguments, and scaling relations all suggest that Xco increases with decreasing metallicity, turning up sharply below metallicity ~1/3-1/2 solar in a manner consistent with model predictions that identify shielding as a key parameter. Based on spectral line modeling and dust observations, Xco appears to drop in the central, bright regions of some but not all galaxies, often coincident with regions of bright CO emission and high stellar surface density. This lower Xco is also present in the overwhelmingly molecular interstellar medium of starburst galaxies, where several lines of evidence point to a lower CO-to-H2 conversion factor. At high redshift, direct evidence regarding the conversion factor remains scarce; we review what is known based on dynamical modeling and other arguments.

The Plateau de Bure + 30m Arcsecond Whirlpool Survey reveals a thick disk of diffuse molecular gas in the M51 galaxy

Abstract: We present the data of the Plateau de Bure Arcsecond Whirlpool Survey (PAWS), a high spatial and spectral resolution 12CO(1-0) line survey of the inner 10 x 6 kpc of the M51 system, and the first wide-field imaging of molecular gas in a star-forming spiral galaxy with resolution matched to the typical size of Giant Molecular Clouds (40 pc). We describe the observation, reduction, and combination of the Plateau de Bure Interferometer (PdBI) and IRAM-30m "short spacing" data. The final data cube attains 1.1"-resolution over the 270" x 170" field of view, with sensitivity to all spatial scales from the combination of PdBI and IRAM-30m data, and brightness sensitivity of 0.4 K (1 sigma) in each 5 km/s-wide channel map. We find a CO-luminosity of 9.10^8 K km/s pc^2, corresponding to a molecular gas mass of 4.10^9 Msol for a standard CO-to-H2 conversion factor. Unexpectedly, we find that a large fraction, (50 +/- 10)%, of this emission arises mostly from spatial scales larger than 36" or about 1.3 kpc. Through a series of tests, we demonstrate that this extended emission does not result from a processing artifact. We discuss its origin in light of the stellar component, the 12CO/13CO ratio, and the difference between the kinematics and structure of the PdBI-only and hybrid synthesis (PdBI + IRAM-30m) images. The extended emission is consistent with a thick, diffuse disk of molecular gas with a typical scale height of about 200 pc, substructured in unresolved filaments which fills about 0.1% of the volume.

Magnetically Aligned HI Fibers and the Rolling Hough Transform

Abstract: We present observations of a new group of structures in the diffuse Galactic ISM: slender, linear HI features we dub "fibers" that extend for many degrees at high Galactic latitude. To characterize and measure the extent and strength of these fibers, we present the Rolling Hough Transform (RHT), a new machine vision method for parameterizing the coherent linearity of structures in the image plane. With this powerful new tool we show the fibers are oriented along the interstellar magnetic field as probed by starlight polarization. We find that these low column density (N(HI) ~ 5 x 10^18 cm^-2) fiber features are most likely a component of the local cavity wall, about 100 pc away. The HI data we use to demonstrate this alignment at high latitude are from the Galactic Arecibo L-Band Feed Array HI (GALFA-HI) Survey and the Parkes Galactic All Sky Survey (GASS). We find better alignment in the higher resolution GALFA-HI data, where the fibers are more visually evident. This trend continues in our investigation of magnetically aligned linear features in the Riegel-Crutcher HI cold cloud, detected in the Southern Galactic Plane Survey (SGPS). We propose an application of the RHT for estimating the field strength in such a cloud, based on the Chandrasekhar-Fermi method. We conclude that data-driven, quantitative studies of ISM morphology can be very powerful predictors of underlying physical quantities.

3D maps of the local ISM from inversion of individual color excess measurements

Abstract: Three-dimensional (3D) maps of the Galactic interstellar matter (ISM) are a potential tool of wide use, however accurate and detailed maps are still lacking. One of the ways to construct the maps is to invert individual distance-limited ISM measurements, a method we have here applied to measurements of stellar color excess in the optical. We have assembled color excess data together with the associated parallax or photometric distances to constitute a catalog of ~ 23,000 sightlines for stars within 2.5 kpc. The photometric data are taken from Stromgren catalogs, the Geneva photometric database, and the Geneva-Copenhagen survey. We also included extinctions derived towards open clusters. We applied, to this color excess dataset, an inversion method based on a regularized Bayesian approach, previously used for mapping at closer distances. We show the dust spatial distribution resulting from the inversion by means of planar cuts through the differential opacity 3D distribution, and by means of 2D maps of the integrated opacity from the Sun up to various distances. The mapping assigns locations to the nearby dense clouds and represents their distribution at the spatial resolution that is allowed by the dataset properties, i.e. of the order of ~10 pc close to the Sun and increasing to ~100 pc beyond 1 kpc. Biases towards nearby and/or weakly extincted stars make this dataset particularly appropriate to map the local and neighboring cavities, and to locate faint, extended nearby clouds, both goals that are difficult or impossible with other mapping methods. The new maps reveal a ~1 kpc wide empty region in the third quadrant in the continuation of the so-called CMa tunnel of the Local Cavity, a cavity that we identify as the Superbubble GSH238+00+09 detected in radio emission maps and that is found to be bounded by the Orion and Vela clouds.

Chemistry of the Sagittarius Dwarf Galaxy: a Top-Light IMF, Outflows and the R-Process

Abstract: From chemical abundance analysis of stars in the Sagittarius dwarf spheroidal galaxy (Sgr), we conclude that the alpha-element deficiencies cannot be due to the Type Ia supernova (SNIa) time-delay scenario of Tinsley (1979). Instead, the evidence points to low [alpha/Fe] ratios resulting from an initial mass function (IMF) deficient in the highest mass stars. The critical evidence is the 0.4 dex deficiency of [O/Fe], [Mg/Fe] and other hydrostatic elements, contrasting with the normal trend of r-process [Eu/Fe]r with [Fe/H]. Supporting evidence comes from the hydrostatic element (O, Mg, Na, Al, Cu) [X/Fe] ratios, which are inconsistent with iron added to the Milky Way (MW) disk trends. Also, the ratio of hydrostatic to explosive (Si, Ca, Ti) element abundances suggests a relatively top-light IMF. Abundance similarities with the LMC, Fornax and IC 1613, suggest that their alpha-element deficiencies also resulted from IMFs lacking the most massive SNII. For such a top-light IMF, the normal trend of r-process [Eu/Fe]r with [Fe/H], as seen in Sgr, indicates that massive Type II supernovae (>30Msun) cannot be major sources of r-process elements. High [La/Y] ratios, consistent with leaky-box chemical evolution, are confirmed but ~0.3 dex larger than theoretical AGB predictions. This may be due to the 13C pocket mass, or a difference between MW and Sgr AGB stars. Sgr has the lowest [Rb/Zr] ratios known, consistent with low-mass (~2Msun) AGB stars near [Fe/H]=-0.6, likely resulting from leaky-box chemical evolution. The [Cu/O] trend in Sgr and the MW suggest that Cu yields increase with both metallicity and stellar mass, as expected from Cu production by the weak s-process in massive stars. Finally, we present an updated hfs line list, an abundance analysis of Arcturus, and further develop our error analysis formalism.

Gaps in the HD169142 protoplanetary disk revealed by polarimetric imaging: Signs of ongoing planet formation?

Abstract: We present H-band VLT/NACO polarized light images of the Herbig Ae/Be star HD169142 probing its protoplanetary disk as close as ~0.1" to the star. Our images trace the face-on disk out to ~1.7" (~250 AU) and reveal distinct sub-structures for the first time: 1) the inner disk (<20 AU) appears to be depleted in scattering dust grains; 2) an unresolved disk rim is imaged at ~25 AU; 3) an annular gap extends from ~40 - 70 AU; 4) local brightness asymmetries are found on opposite sides of the annular gap. We discuss different explanations for the observed morphology among which ongoing planet formation is a tempting - but yet to be proven - one. Outside of ~85 AU the surface brightness drops off roughly r^{-3.3}, but describing the disk regions between 85-120 AU / 120-250 AU separately with power-laws r^{-2.6} / r^{-3.9} provides a better fit hinting towards another discontinuity in the disk surface. The flux ratio between the disk integrated polarized light and the central star is ~4.1 * 10^{-3}. Finally, combining our results with those from the literature, ~40% of the scattered light in the H-band appears to be polarized. Our results emphasize that HD169142 is an interesting system for future planet formation or disk evolution studies.

Planck intermediate results. XVII. Emission of dust in the diffuse interstellar medium from the far-infrared to microwave frequencies

Abstract: The dust-HI correlation is used to characterize the emission properties of dust in the diffuse interstellar medium. We cross-correlate sky maps from Planck, WMAP, and DIRBE, at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21-cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg$^2$ centred on the southern Galactic pole. Our analysis yields four specific results. (1) The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence for dust evolution within the diffuse ISM. The mean dust opacity is measured to be $(7.1 \pm 0.6) 10^{-27} cm^2/H \times ( u/353\, GHz)^{1.53\pm0.03}$ for $100 < u <353$GHz. (2) We map the spectral index of dust emission at millimetre wavelengths, which is remarkably constant at $\beta_{mm} = 1.51\pm 0.13$. We compare it with the far infrared spectral index beta_FIR derived from greybody fits at higher frequencies, and find a systematic difference, $\beta_{mm}-\beta_{FIR} = -0.15$, which suggests that the dust SED flattens at $ u < 353\,$GHz. (3) We present spectral fits of the microwave emission correlated with HI from 23 to 353 GHz, which separate dust and anomalous microwave emission. The flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum, which accounts for $(26 \pm 6)$% of the dust emission at 100 GHz and could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We identify a Galactic contribution to the residuals of the dust-HI correlation, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis.

Chemical Cartography with APOGEE: Large-scale Mean Metallicity Maps of the Milky Way

Abstract: We present Galactic mean metallicity maps derived from the first year of the SDSS-III APOGEE experiment. Mean abundances in different zones of Galactocentric radius (0 6 kpc, the gradient flattens as one moves off of the plane, and is flatter at all heights for high-[{\alpha}/M] stars than for low-[{\alpha}/M] stars. Alternatively, these gradients can be described as vertical gradients that flatten at larger Galactocentric radius; these vertical gradients are similar for both low and high-[{\alpha}/M] populations. Stars with higher [{\alpha}/M] appear to have a flatter radial gradient than stars with lower [{\alpha}/M]. This could suggest that the metallicity gradient has grown steeper with time or, alternatively, that gradients are washed out over time by migration of stars.

The Origin and Optical Depth of Ionizing Radiation in the "Green Pea" Galaxies

Abstract: Although Lyman continuum (LyC) radiation from star-forming galaxies likely drove the reionization of the Universe, observations of star-forming galaxies at low redshift generally indicate low LyC escape fractions. However, the extreme [O III]/[O II] ratios of the z=0.1-0.3 Green Pea galaxies may be due to high escape fractions of ionizing radiation. To analyze the LyC optical depths and ionizing sources of these rare, compact starbursts, we compare nebular photoionization and stellar population models with observed emission lines in the Peas' SDSS spectra. We focus on the six most extreme Green Peas, the galaxies with the highest [O III]/[O II] ratios and the best candidates for escaping ionizing radiation. The Balmer line equivalent widths and He I {\lambda}3819 emission in the extreme Peas support young ages of 3-5 Myr, and He II {\lambda}4686 emission in five extreme Peas signals the presence of hard ionizing sources. Ionization by active galactic nuclei or high-mass X-ray binaries is inconsistent with the Peas' line ratios and ages. Although stacked spectra reveal no Wolf-Rayet (WR) features, we tentatively detect WR features in the SDSS spectra of three extreme Peas. Based on the Peas' ages and line ratios, we find that WR stars, chemically homogeneous O stars, or shocks could produce the observed He II emission. If hot stars are responsible, the Peas' optical depths are ambiguous. However, accounting for emission from shocks lowers the inferred optical depth and suggests that the Peas may be optically thin. The Peas' ages likely optimize the escape of Lyman-continuum radiation; they are old enough for supernovae and stellar winds to reshape the interstellar medium, but young enough to possess large numbers of UV-luminous O or WR stars.

Rotationally-supported disks around Class I sources in Taurus: disk formation constraints

Abstract: (Abridged) Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. We present subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus. The 13CO and C18O J=2-1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ~0.8'' and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity radient. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are 100 AU around these sources are dominated by infalling motion (v propto r^-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of 10^-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. RSDs with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.