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

Measuring Reddening with SDSS Stellar Spectra and Recalibrating SFD

Abstract: We present measurements of dust reddening using the colors of stars with spectra in the Sloan Digital Sky Survey. We measure reddening as the difference between the measured and predicted colors of a star, as derived from stellar parameters from the SEGUE Stellar Parameter Pipeline (Lee et al. 2008a). We achieve uncertainties of 56, 34, 25, and 29 mmag in the colors u-g, g-r, r-i, and i-z, per star, though the uncertainty varies depending on the stellar type and the magnitude of the star. The spectrum-based reddening measurements confirm our earlier "blue tip" reddening measurements (Schlafly et al. 2010, S10), finding reddening coefficients different by -3%, 1%, 1%, and 2% in u-g, g-r, r-i, and i-z from those found by the blue tip method, after removing a 4% normalization difference. These results prefer an R_V=3.1 Fitzpatrick (1999) reddening law to O'Donnell (1994) or Cardelli et al. (1989) reddening laws. We provide a table of conversion coefficients from the Schlegel et al. (1998) maps of E(B-V) to extinction in 88 bandpasses for 4 values of R_V, using this reddening law and the 14% recalibration of SFD first reported by S10 and confirmed in this work.

The PPMXL catalog of positions and proper motions on the ICRS. Combining USNO-B1.0 and 2MASS

Abstract: USNO-B1.0 and 2MASS are the most widely used full-sky surveys. However, 2MASS has no proper motions at all, and USNO-B1.0 published only relative, not absolute (i.e. on ICRS) proper motions. We performed a new determination of mean positions and proper motions on the ICRS system by combining USNO-B1.0 and 2MASS astrometry. This catalog is called PPMXL {VO-access to the catalog is possible via this http URL}, and it aims to be complete from the brightest stars down to about $V \approx 20$ full-sky. PPMXL contains about 900 million objects, some 410 million with 2MASS photometry, and is the largest collection of ICRS proper motions at present. As representative for the ICRS we chose PPMX. The recently released UCAC3 could not be used because we found plate-dependent distortions in its proper motion system north of -20$^\circ$ declination. UCAC3 served as an intermediate system for $\delta \leq -20^\circ$. The resulting typical individual mean errors of the proper motions range from 4 mas/y to more than 10 mas/y depending on observational history. The mean errors of positions at epoch 2000.0 are 80 to 120 mas, if 2MASS astrometry could be used, 150 to 300 mas else. We also give correction tables to convert USNO-B1.0 observations of e.g. minor planets to the ICRS system.

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

Abstract: We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

Hi-GAL: the Herschel infrared Galactic Plane Survey

Abstract: Hi-GAL, the Herschel infrared Galactic Plane Survey, is an Open Time Key Project of the Herschel Space Observatory. It will make an unbiased photometric survey of the inner Galactic Plane by mapping a two-degree wide strip in the longitude range |l|<60 degrees in five wavebands between 70um and 500um. The aim of Hi-GAL is to detect the earliest phases of the formation of molecular clouds and high-mass stars and to use the optimum combination of Herschel wavelength coverage, sensitivity, mapping strategy and speed to deliver a homogeneous census of star-forming regions and cold structures in the interstellar medium. The resulting representative samples will yield the variation of source temperature, luminosity, mass and age in a wide range of Galactic environments at all scales from massive YSOs in protoclusters to entire spiral arms, providing an evolutionary sequence for the formation of intermediate and high-mass stars. This information is essential to the formulation of a predictive global model of the role of environment and feedback in regulating the star-formation process. Such a model is vital to understanding star formation on galactic scales and in the early Universe. Hi-GAL will also provide a science legacy for decades to come with incalculable potential for systematic and serendipitous science in a wide range of astronomical fields, enabling the optimum use of future major facilities such as JWST and ALMA.

Unification of Luminous Type 1 Quasars through CIV Emission

Abstract: Using a sample of 30,000 quasars from SDSS-DR7, we explore the range of properties exhibited by high-ionization, broad emission lines, such as CIV 1549. Specifically we investigate the anti-correlation between L_UV and emission line EQW (the Baldwin Effect) and the "blueshifting" of high-ionization emission lines. The blueshift of the CIV emission line is nearly ubiquitous, with a mean shift of 810 km/s for radio-quiet (RQ) quasars and 360 km/s for radio-loud (RL) quasars, and the Baldwin Effect is present in both RQ and RL samples. Composite spectra are constructed as a function of CIV emission line properties in attempt to reveal empirical relationships between different line species and the SED. Within a two-component disk+wind model of the broad emission line region (BELR), where the wind filters the continuum seen by the disk component, we find that RL quasars are consistent with being dominated by the disk component, while BALQSOs are consistent with being dominated by the wind component. Some RQ objects have emission line features similar to RL quasars; they may simply have insufficient black hole (BH) spin to form radio jets. Our results suggest that there could be significant systematic errors in the determination of L_bol and BH mass that make it difficult to place these findings in a more physical context. However, it is possible to classify quasars in a paradigm where the diversity of BELR parameters are due to differences in an accretion disk wind between quasars (and over time); these differences are underlain primarily by the SED, which ultimately must be tied to BH mass and accretion rate.

The Sagittarius Dwarf Galaxy: a Model for Evolution in a Triaxial Milky Way Halo

Abstract: We present a new N-body model for the tidal disruption of the Sagittarius (Sgr) dwarf that is capable of simultaneously satisfying the majority of angular position, distance, and radial velocity constraints imposed by current wide-field surveys of its dynamically young (< 3 Gyr) tidal debris streams. In particular, this model resolves the conflicting angular position and radial velocity constraints on the Sgr leading tidal stream that have been highlighted in recent years. While the model does not reproduce the apparent bifurcation observed in the leading debris stream, recent observational data suggest that this bifurcation may represent a constraint on the internal properties of the Sgr dwarf rather than the details of its orbit. The key element in the success of this model is the introduction of a non-axisymmetric component to the Galactic gravitational potential which can be described in terms of a triaxial dark matter halo whose minor/major axis ratio (c/a)_Phi = 0.72 and intermediate/major axis ratio (b/a)_Phi = 0.99 at radii 20 < r < 60 kpc. The minor/intermediate/major axes of this halo lie along the directions (l, b) = (7, 0), (0, 90), and (97, 0) respectively, corresponding to a nearly-oblate ellipsoid whose minor axis is contained within the Galactic disk plane. We demonstrate that with simple assumptions about the star formation history of Sgr, tidal stripping models naturally give rise to gradients in the metallicity distribution function (MDF) along the stellar debris streams similar to those observed in recent studies. (Abridged).

The metal-poor end of the Spite plateau. 1: Stellar parameters, metallicities and lithium abundances

Abstract: We present VLT-UVES Li abundances for 28 halo dwarf stars between [Fe/H]=-2.5 and -3.5, 10 of which have [Fe/H]<-3. Four different T_eff scales have been used. Direct Infrared Flux Method (IRFM) has been used on the basis of 2MASS infrared photometry. H_alpha wings have been fitted against synthetic grids computed by means of 1D LTE atmosphere models, assuming different self-broadening theories. Finally, a grid of H_alpha profiles has been computed by means of 3D hydrodynamical atmosphere models. The Li I doublet at 670.8 nm has been used to measure A(Li) by means of 3D hydrodynamical NLTE spectral syntheses. An analytical fit of A(Li)(3D, NLTE) as a function of equivalent width, T_eff, log g, and [Fe/H] has been derived and is made available. A(Li) does not exhibit a plateau below [Fe/H]=-3. A strong positive correlation with [Fe/H] appears, not influenced by the choice of the T_eff estimator. From a linear fit, we obtain a strong slope of about 0.30 dex in A(Li) per dex in [Fe/H], significant to 2-3 sigma, and consistent among all the four T_eff estimators. A significant slope is also detected in the A(Li)--T_eff plane, driven mainly by the coolest stars in the sample which appear Li-poor. Removing such stars does not alter the behavior in the A(Li)-[Fe/H] plane. The scatter in A(Li) increases by a factor of 2 towards lower metallicities, while the plateau appears very thin above [Fe/H]=-2.8. The meltdown of the Spite plateau below [Fe/H]\sim-3 is established, but its cause is unclear. If the primordial A(Li) is the one derived from WMAP, it appears difficult to envision a single depletion phenomenon producing a thin, metallicity independent plateau above [Fe/H]=-2.8, and a highly scattered, metallicity dependent distribution below. The fact that no star below [Fe/H]=-3 lies above the plateau suggests that they formed at plateau level and underwent subsequent depletion.

The Masses of the Milky Way and Andromeda galaxies

Abstract: We present a family of robust tracer mass estimators to compute the enclosed mass of galaxy haloes from samples of discrete positional and kinematical data of tracers, such as halo stars, globular clusters and dwarf satellites. The data may be projected positions, distances, line of sight velocities or proper motions. Forms of the estimator tailored for the Milky Way galaxy and for M31 are given. Monte Carlo simulations are used to quantify the uncertainty as a function of sample size. For the Milky Way, the satellite sample consists of 26 galaxies with line-of-sight velocities. We find that the mass of the Milky Way within 300 kpc is ~ 0.9 x 10^12 solar masses assuming velocity isotropy. However, the mass estimate is sensitive to the anisotropy and could plausibly lie between 0.7 - 3.4 x 10^12 solar masses. Incorporating the proper motions of 6 Milky Way satellites into the dataset, we find ~ 1.4 x 10^12 solar masses. The range here if plausible anisotropies are used is still broader, from 1.2 - 2.7 x 10^12 solar masses. For M31, there are 23 satellite galaxies with measured line-of-sight velocities, but only M33 and IC 10 have proper motions. We use the line of sight velocities and distances of the satellite galaxies to estimate the mass of M31 within 300 kpc as ~ 1.4 x 10^12 solar masses assuming isotropy. There is only a modest dependence on anisotropy, with the mass varying between 1.3 -1.6 x 10^12 solar masses. Given the uncertainties, we conclude that the satellite data by themselves yield no reliable insights into which of the two galaxies is actually the more massive.

The dark matter density at the Sun's location

Abstract: We derive the value of the dark matter density at the Sun's location (rho_0) without globally mass-modeling the Galaxy. The proposed method relies on the local equation of centrifugal equilibrium and is independent of i) the shape of the dark matter density profile, ii) knowledge of the rotation curve from the galaxy center out to the virial radius, and iii) the uncertainties and the non-uniqueness of the bulge/disk/dark halo mass decomposition. The result can be obtained in analytic form. It explicitly includes the dependence on the relevant observational quantities and takes their uncertainties into account. By adopting the reference, state-of-the-art values for these, we find rho_0 = 0.43(11)(10)GeV/cm3, where the quoted uncertainties are respectively due to the uncertainty in the slope of the circular-velocity at the Sun location and the ratio between this radius and the length scale of the stellar exponential thin disk. We obtain a reliable estimate of rho_0, that, in addition, is ready to take into account any future change/improvement in the measures of the observational quantities it depends on.

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

Abstract: We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

Timing Measurements of the Relativistic Binary Pulsar PSR B1913+16

Abstract: We present results of more than three decades of timing measurements of the first known binary pulsar, PSR B1913+16. Like most other pulsars, its rotational behavior over such long time scales is significantly affected by small-scale irregularities not explicitly accounted for in a deterministic model. Nevertheless, the physically important astrometric, spin, and orbital parameters are well determined and well decoupled from the timing noise. We have determined a significant result for proper motion, $\mu_{\alpha} = -1.43\pm0.13$, $\mu_{\delta}=-0.70\pm0.13$ mas yr$^{-1}$. The pulsar exhibited a small timing glitch in May 2003, with ${\Delta f}/f=3.7\times10^{-11}$, and a smaller timing peculiarity in mid-1992. A relativistic solution for orbital parameters yields improved mass estimates for the pulsar and its companion, $m_1=1.4398\pm0.0002 \ M_{\sun}$ and $m_2=1.3886\pm0.0002 \ M_{\sun}$. The system's orbital period has been decreasing at a rate $0.997\pm0.002$ times that predicted as a result of gravitational radiation damping in general relativity. As we have shown before, this result provides conclusive evidence for the existence of gravitational radiation as predicted by Einstein's theory.

Momentum driving: which physical processes dominate AGN feedback?

Abstract: The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $\eta=dot{M_outflow}/dot{M_accretion}= (2 \epsilon_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $\eta$, and many prior treatments have essentially assumed that $\eta=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $\epsilon_w \lesssim 10^{-4}$.

C$_{60}$ in Reflection Nebulae

Abstract: The fullerene C$_{60}$ has four infrared-active vibrational transitions at 7.0, 8.5, 17.4 and 18.9 $\mu$m. We have previously observed emission features at 17.4 and 18.9 $\mu$m in the reflection nebula NGC 7023 and demonstrated spatial correlations suggestive of a common origin. We now confirm our earlier identification of these features with C$_{60}$ by detecting a third emission feature at 7.04 $\pm$ 0.05 $\mu$m in NGC 7023. We also report the detection of these three C$_{60}$ features in the reflection nebula NGC 2023. Our spectroscopic mapping of NGC 7023 shows that the 18.9 $\mu$m C$_{60}$ feature peaks on the central star and that the 16.4 $\mu$m emission feature due to polycyclic aromatic hydrocarbons peaks between the star and a nearby photodissociation front. The observed features in NGC 7023 are consistent with emission from UV-excited gas-phase C$_{60}$. We find that 0.1-0.6% of interstellar carbon is in C$_{60}$; this abundance is consistent with those from previous upper limits and possible fullerene detections in the interstellar medium. This is the first firm detection of neutral C$_{60}$ in the interstellar medium.

The global dust SED: Tracing the nature and evolution of dust with DustEM

Abstract: The Planck and Herschel missions are currently measuring the farIR-mm emission of dust, which combined with existing IR data, will for the first time provide the full SED of the galactic ISM dust emission with an unprecedented sensitivity and angular resolution. It will allow a systematic study of the dust evolution processes that affect the SED. Here we present a versatile numerical tool, DustEM, that predicts the emission and extinction of dust given their size distribution and their optical and thermal properties. In order to model dust evolution, DustEM has been designed to deal with a variety of grain types, structures and size distributions and to be able to easily include new dust physics. We use DustEM to model the dust SED and extinction in the diffuse interstellar medium at high-galactic latitude (DHGL), a natural reference SED. We present a coherent set of observations for the DHGL SED. The dust components in our DHGL model are (i) PAHs, (ii) amorphous carbon and (iii) amorphous silicates. We use amorphous carbon dust, rather than graphite, because it better explains the observed high abundances of gas-phase carbon in shocked regions of the interstellar medium. Using the DustEM model, we illustrate how, in the optically thin limit, the IRAS/Planck HFI (and likewise Spitzer/Herschel for smaller spatial scales) photometric band ratios of the dust SED can disentangle the influence of the exciting radiation field intensity and constrain the abundance of small grains relative to the larger grains. We also discuss the contributions of the different grain populations to the IRAS, Planck and Herschel channels. Such information is required to enable a study of the evolution of dust as well as to systematically extract the dust thermal emission from CMB data and to analyze the emission in the Planck polarized channels. The DustEM code described in this paper is publically available.

Water in Star-Forming Regions with the Herschel Space Observatory (WISH): Overview of key program and first results

Abstract: `Water In Star-forming regions with Herschel' (WISH) is a key program on the Herschel Space Observatory designed to probe the physical and chemical structure of young stellar objects using water and related molecules and to follow the water abundance from collapsing clouds to planet-forming disks. About 80 sources are targeted covering a wide range of luminosities and evolutionary stages, from cold pre-stellar cores to warm protostellar envelopes and outflows to disks around young stars. Both the HIFI and PACS instruments are used to observe a variety of lines of H2O, H218O and chemically related species. An overview of the scientific motivation and observational strategy of the program is given together with the modeling approach and analysis tools that have been developed. Initial science results are presented. These include a lack of water in cold gas at abundances that are lower than most predictions, strong water emission from shocks in protostellar environments, the importance of UV radiation in heating the gas along outflow walls across the full range of luminosities, and surprisingly widespread detection of the chemically related hydrides OH+ and H2O+ in outflows and foreground gas. Quantitative estimates of the energy budget indicate that H2O is generally not the dominant coolant in the warm dense gas associated with protostars. Very deep limits on the cold gaseous water reservoir in the outer regions of protoplanetary disks are obtained which have profound implications for our understanding of grain growth and mixing in disks.

A catalogue of young runaway Hipparcos stars within 3kpc from the Sun

Abstract: Traditionally runaway stars are O and B type stars with large peculiar velocities.We want to extend this definition to young stars (up to ~50 Myr) of any spectral type and identify those present in the Hipparcos catalogue applying different selection criteria such as peculiar space velocities or peculiar one-dimensional velocities. Runaway stars are important to study the evolution of multiple star systems or star clusters as well as to identify origins of neutron stars. We compile distances, proper motions, spectral types, luminosity classes, V magnitudes and B-V colours and utilise evolutionary models from different authors to obtain star ages and study a sample of 7663 young Hipparcos stars within 3 kpc from the Sun. Radial velocities are obtained from the literature. We investigate the distributions of the peculiar spatial velocity, the peculiar radial velocity as well as the peculiar tangential velocity and its one-dimensional components and obtain runaway star probabilities for each star in the sample. In addition, we look for stars that are situated outside any OB association or OB cluster and the Galactic plane as well as stars of which the velocity vector points away from the median velocity vector of neighbouring stars or the surrounding local OB association/ cluster although the absolute velocity might be small. We find a total of 2547 runaway star candidates (with a contamination of normal Population I stars of 20 per cent at most). Thus, after subtraction of those 20 per cent, the runaway frequency among young stars is about 27 per cent. We compile a catalogue of runaway stars which will be available via VizieR.

Physical Properties and Galactic Distribution of Molecular Clouds identified in the Galactic Ring Survey

Abstract: We derive the physical properties of 580 molecular clouds based on their 12CO and 13CO line emission detected in the University of Massachusetts-Stony Brook (UMSB) and Galactic Ring surveys. We provide a range of values of the physical properties of molecular clouds, and find a power-law correlation between their radii and masses, suggesting that the fractal dimension of the ISM is around 2.36. This relation, M = (228 +/- 18) R^{2.36+/-0.04}, allows us to derive masses for an additional 170 GRS molecular clouds not covered by the UMSB survey. We derive the Galactic surface mass density of molecular gas and examine its spatial variations throughout the Galaxy. We find that the azimuthally averaged Galactic surface density of molecular gas peaks between Galactocentric radii of 4 and 5 kpc. Although the Perseus arm is not detected in molecular gas, the Galactic surface density of molecular gas is enhanced along the positions of the Scutum-Crux and Sagittarius arms. This may indicate that molecular clouds form in spiral arms and are disrupted in the inter-arm space. Last, we find that the CO excitation temperature of molecular clouds decreases away from the Galactic center, suggesting a possible decline in the star formation rate with Galactocentric radius. There is a marginally significant enhancement in the CO excitation temperature of molecular clouds at a Galactocentric radius of about 6 kpc, which in the longitude range of the GRS corresponds to the Sagittarius arm. This temperature increase could be associated with massive star formation in the Sagittarius spiral arm.

Modeling CO Emission: I. CO as a Column Density Tracer and the X-Factor in Molecular Clouds

Abstract: Theoretical and observational investigations have indicated that the abundance of carbon monoxide (CO) is very sensitive to intrinsic properties of the gaseous medium, such as density, metallicity, and the background UV field. In order to accurately interpret CO observations, it is thus important to understand how well CO traces the gas, which in molecular clouds (MCs) is predominantly molecular hydrogen (H2). Recent hydrodynamic simulations by Glover & Mac Low have explicitly followed the formation and destruction of molecules in model MCs under varying conditions, confirming that CO formation strongly depends on the cloud properties. Conversely, the H2 formation is primarily determined by the age of the MC. We apply radiative transfer calculations to these MC models in order to investigate the properties of CO line emission. We focus on integrated CO (J=1-0) intensities emerging from individual clouds, including its relationship to the total, H2, and CO column densities, as well as the "X factor," the ratio of H2 column density to CO intensity. Models with high CO abundances have a threshold CO intensity ~65 K km/s at sufficiently large extinctions. Clouds with low CO abundances show no such intensity thresholds. The distribution of H2 column densities are well described as log-normal functions, though the distributions of CO intensities and column densities are usually not log-normal. In general, the PDFs of the integrated intensity do not follow the distribution functions of CO column densities. In the model with Milky Way-like conditions, the X factor is in agreement with the near constant value determined from observations. In clouds with lower CO abundances the X factor can vary appreciably - sometimes by > 4 orders of magnitude. In models with high densities, the CO line is fully saturated, so that the X factor is directly proportional to the molecular column density.

The Genesis of the Milky Way's Thick Disk via Stellar Migration

Abstract: The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular orbits, as described by Sellwood & Binney (2002). As stars move outwards they populate a thickened component. Such stars found at the present time in the solar neighborhood formed early in the disk's history at smaller radii where stars are more metal-poor and alpha-enhanced, leading to the properties observed for thick disk stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to an apparent separation in the other property as observed. This scenario is supported by the SDSS observation that stars in the transition region do not show any correlation between rotational velocity and metallicity. The good qualitative agreement between our simulation and observations in the Milky Way hint that the thick disk may be a ubiquitous galaxy feature generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies.

Simulations of the Magellanic Stream in a First Infall Scenario

Abstract: Recent high precision proper motions from the Hubble Space Telescope (HST) suggest that the Large and Small Magellanic Clouds (LMC and SMC, respectively) are either on their first passage or on an eccentric long period (>6 Gyr) orbit about the Milky Way (MW). This differs markedly from the canonical picture in which the Clouds travel on a quasi-periodic orbit about the MW (period of ~2 Gyr). Without a short period orbit about the MW, the origin of the Magellanic Stream, a young (1-2 Gyr old) coherent stream of HI gas that trails the Clouds ~150 degrees across the sky, can no longer be attributed to stripping by MW tides and/or ram pressure stripping by MW halo gas. We propose an alternative formation mechanism in which material is removed by LMC tides acting on the SMC before the system is accreted by the MW. We demonstrate the feasibility and generality of this scenario using an N-body/SPH simulation with cosmologically motivated initial conditions constrained by the observations. Under these conditions we demonstrate that it is possible to explain the origin of the Magellanic Stream in a first infall scenario. This picture is generically applicable to any gas-rich dwarf galaxy pair infalling towards a massive host or interacting in isolation.

A New Catalog of Globular Clusters in the Milky Way

Abstract: A new revision of the McMaster catalog of Milky Way globular clusters is available. This is the first update since 2003 and the biggest single revision since the original version of the catalog published in 1996. The list now contains a total of 157 objects classified as globular clusters. Major upgrades have been made especially to the cluster coordinates, metallicities, and structural profile parameters, and the list of parameters now also includes central velocity dispersion. NB: This paper is a stand-alone publication available only on the astro-ph archive; it will not be published separately in a journal.

The Bolocam Galactic Plane Survey: Survey Description and Data Reduction

Abstract: We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

Chemistry in Disks. IV. Benchmarking gas-grain chemical models with surface reactions

Abstract: Abridged: We detail and benchmark two sophisticated chemical models developed by the Heidelberg and Bordeaux astrochemistry groups. The main goal of this study is to elaborate on a few well-described tests for state-of-the-art astrochemical codes covering a range of physical conditions and chemical processes, in particular those aimed at constraining current and future interferometric observations of protoplanetary disks. We consider three physical models: a cold molecular cloud core, a hot core, and an outer region of a T Tauri disk. Our chemical network (for both models) is based on the original gas-phase osu_03_2008 ratefile and includes gas-grain interactions and a set of surface reactions for the H-, O-, C-, S-, and N-bearing molecules. The benchmarking is performed with the increasing complexity of the considered processes: (1) the pure gas-phase chemistry, (2) the gas-phase chemistry with accretion and desorption, and (3) the full gas-grain model with surface reactions. Using atomic initial abundances with heavily depleted metals and hydrogen in its molecular form, the chemical evolution is modeled within 10^9 years. The time-dependent abundances calculated with the two chemical models are essentially the same for all considered physical cases and for all species, including the most complex polyatomic ions and organic molecules. This result however required a lot of efforts to make all necessary details consistent through the model runs, e.g. definition of the gas particle density, density of grain surface sites, the strength and shape of the UV radiation field, etc. The reference models and the benchmark setup, along with the two chemical codes and resulting time-dependent abundances are made publicly available in the Internet: this http URL

Big fish, small fish: Two New Ultra-Faint Satellites of the Milky Way

Abstract: We report the discovery of two new Milky Way satellites in the neighboring constellations of Pisces and Pegasus identified in data from the Sloan Digital Sky Survey. Pisces II, an ultra-faint dwarf galaxy lies at the distance of ~180 kpc, some 15 degrees away from the recently detected Pisces I. Segue 3, an ultra-faint star cluster lies at the distance of 16 kpc. We use deep follow-up imaging obtained with the 4-m Mayall telescope at Kitt Peak National Observatory to derive their structural parameters. Pisces II has a half-light radius of ~60 pc, while Segue 3 is twenty times smaller at only 3pc.

Dynamic star formation in the massive DR21 filament

Abstract: The formation of massive stars is a highly complex process in which it is not clear whether the star-forming gas is in global gravitational collapse or in an equilibrium state, supported by turbulence. By studying one of the most massive and dense star-forming regions in the Galaxy at a distance of less than 3 kpc, the filament containing the well-known sources DR21 and DR21(OH), we expect to find observational signatures that allow to discriminate between the two views. We use molecular line data from our 13CO 1-0, CS 2-1, and N2H+ 1-0 survey of the Cygnus X region obtained with the FCRAO and high-angular resolution observations of CO, CS, HCO+, N2H+, and H2CO, obtained with the IRAM 30m telescope. We observe a complex velocity field and velocity dispersion in the DR21 filament in which regions of highest column-density, i.e. dense cores, have a lower velocity dispersion than the surrounding gas and velocity gradients that are not (only) due to rotation. Infall signatures in optically thick line profiles of HCO+ and 12CO are observed along and across the whole DR21 filament. From modelling the observed spectra, we obtain a typical infall speed of 0.6 km/s and mass accretion rates of the order of a few 10^-3 Msun/yr for the two main clumps constituting the filament. These massive (4900 and 3300 Msun) clumps are both gravitationally contracting. All observed kinematic features in the DR21 filament can be explained if it is formed by the convergence of flows at large scales and is now in a state of global gravitational collapse. Whether this convergence of flows originated from self-gravity at larger scales or from other processes can not be settled with the present study. The observed velocity field and velocity dispersion are consistent with results from (magneto)-hydrodynamic simulations where the cores lie at the stagnation points of convergent turbulent flows.

Gravity or Turbulence? The velocity dispersion-size relation

Abstract: We discuss the nature of the velocity dispersion vs. size relation for molecular clouds. In particular, we add to previous observational results showing that the velocity dispersions in molecular clouds and cores are not purely functions of spatial scale but involve surface gas densities as well. We emphasize that hydrodynamic turbulence is required to produce the first condensations in the progenitor medium. However, as the cloud is forming, it also becomes bound, and gravitational accelerations dominate the motions. Energy conservation in this case implies $|E_g| \sim E_k$, in agreement with observational data, and providing an interpretation for two recent observational results: the scatter in the $\delta v-R$ plane, and the dependence of the velocity dispersion on the surface density ${\delta v^2/ R} \propto \Sigma$. We argue that the observational data are consistent with molecular clouds in a state of hierarchical gravitational collapse, i.e., developing local centers of collapse throughout the whole cloud while the cloud itself is collapsing, and making equilibrium unnecessary at all stages prior to the formation of actual stars. Finally, we discuss how this mechanism need not be in conflict with the observed star formation rate.

The 200 Degree-Long Magellanic Stream System

Abstract: We establish that the Magellanic Stream (MS) is some 40 degrees longer than previously known with certainty and that the entire MS and Leading Arm (LA) system is thus at least 200 degrees long. With the GBT, we conducted a ~200 square degree, 21-cm survey at the MS-tip to substantiate the continuity of the MS between the Hulsbosch & Wakker data and the MS-like emission reported by Braun & Thilker. Our survey, in combination with the Arecibo survey by Stanimirovic et al., shows that the MS gas is continuous in this region and that the MS is at least ~140 degrees long. We identify a new filament on the eastern side of the MS that significantly deviates from the equator of the MS coordinate system for more than ~45 degrees. Additionally, we find a previously unknown velocity inflection in the MS-tip near MS longitude L_MS=-120 degrees at which the velocity reaches a minimum and then starts to increase. We find that five compact high velocity clouds cataloged by de Heij et al. as well as Wright's Cloud are plausibly associated with the MS because they match the MS in position and velocity. The mass of the newly-confirmed ~40 degree extension of the MS-tip is ~2x10^7 Msun (d/120 kpc)^2 (including Wright's Cloud increases this by ~50%) and increases the total mass of the MS by ~4%. However, projected model distances of the MS at the tip are generally quite large and, if true, indicate that the mass of the extension might be as large as ~10^8 Msun. From our combined map of the entire MS, we find that the total column density (integrated transverse to the MS) drops markedly along the MS and follows an exponential decline with L_MS. We estimate that the age of the ~140 degree-long MS is ~2.5 Gyr which coincides with bursts of star formation in the Magellanic Clouds and a possible close encounter of these two galaxies with each other that could have triggered the formation of the MS. [Abridged]

On Radiation Pressure in Static, Dusty HII Regions

Abstract: Radiation pressure acting on gas and dust causes HII regions to have central densities that are lower than the density near the ionized boundary. HII regions in static equilibrium comprise a family of similarity solutions, parametrized by 3 parameters: beta, gamma, and the product (Q_0 n_rms); beta characterizes the stellar spectrum, gamma characterizes the dust/gas ratio, Q_0 is the ionizing output from the star (photons/s), and n_rms is the rms density within the ionized region. Adopting standard values for beta and gamma, varying (Q_0 n_rms) generates a one-parameter family of density profiles, ranging from nearly uniform density (small Q_0 n_rms), to hollow-sphere HII regions (large Q_0 n_rms). When (Q_0 n_rms) exceeds 10^{52} cm^{-3} s^{-1}, dusty HII regions have conspicuous central cavities, even if no stellar wind is present. For given beta, gamma and (Q_0 n_rms), a fourth quantity, which can be Q_0, determines the overall size and density of the HII region. Examples of density and emissivity profiles are given. We show how quantities of interest -- such as the peak-to-center emissivity ratio, the rms-to-mean density ratio, the edge-to-rms density ratio, and the fraction of the ionizing photons absorbed by the gas -- depend on the 3 parameters beta, gamma, and (Q_0 n_rms). For dusty HII regions, compression of the gas and dust into an ionized shell results in a substantial increase in the fraction of the >13.6 eV photons that actually ionize H (relative to a uniform density HII region with the same dust/gas ratio and density n=n_rms). We discuss the extent to which radial drift of dust grains in HII regions can alter the dust-to-gas ratio. The applicability of these solutions to real HII regions is discussed.

A Multi-color Optical Survey of the Orion Nebula Cluster. II. the H-R diagram

Abstract: We present a new analysis of the stellar population of the Orion Nebula Cluster (ONC) based on multi-band optical photometry and spectroscopy. We study the color-color diagrams in BVI, plus a narrow-band filter centered at 6200A, finding evidences that intrinsic color scales valid for main-sequence dwarfs are incompatible with the ONC, while a better agreement is found employing synthetic intrinsic colors obtained constraining the typical lower surface gravity of young stars. We refine these model colors even further, empirically, by comparison with a sample of ONC stars with no accretion and no extinction. We consider the stars with known spectral types from the literature, and add 65 newly classified stars from slit spectroscopy and 182 M-type from narrow-band photometry; in this way we isolate a sample of about 1000 stars with known spectral type. We introduce a new method to self-consistently derive reddening and accretion excess from the location of each star in the BVI color-color diagram. This enables us to accurately determine the extinction of the ONC members. We adopt a lower distance for the ONC than previously assumed, based on recent parallax measurements. With a careful choice also of the spectral type-temperature transformation, we produce the new H-R diagram of the ONC population, more populated than previous works. With respect to previous works, we find higher luminosity for late-type stars and a lower luminosity for early types. We determine the age distribution of the population, peaking at 2-3 Myr, a higher age than previously estimated. We study the distribution of the members in the mass-age plane, and find that taking into account selection effects due to incompleteness removes an apparent correlation between mass and age. We derive the IMF for low- and intermediate-mass members of the ONC, which turns out to be model-dependent, and shows a turn-over at ~<0.2Msun.

Evidence for an accretion origin for the outer halo globular cluster system of M31

Abstract: We use a sample of newly-discovered globular clusters from the Pan-Andromeda Archaeological Survey (PAndAS) in combination with previously-catalogued objects to map the spatial distribution of globular clusters in the M31 halo. At projected radii beyond ~30 kpc, where large coherent stellar streams are readily distinguished in the field, there is a striking correlation between these features and the positions of the globular clusters. Adopting a simple Monte Carlo approach, we test the significance of this association by computing the probability that it could be due to the chance alignment of globular clusters smoothly distributed in the M31 halo. We find the likelihood of this possibility is low, below 1%, and conclude that the observed spatial coincidence between globular clusters and multiple tidal debris streams in the outer halo of M31 reflects a genuine physical association. Our results imply that the majority of the remote globular cluster system of M31 has been assembled as a consequence of the accretion of cluster-bearing satellite galaxies. This constitutes the most direct evidence to date that the outer halo globular cluster populations in some galaxies are largely accreted.