Showing papers in "Astronomy and Astrophysics in 2013"

Astropy: A community Python package for astronomy

Abstract: We present the first public version (v02) of the open-source and community-developed Python package, Astropy This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as flexible image transport system (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions Significant functionality is under activedevelopment, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions

LOFAR: The LOw-Frequency ARray

Abstract: LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

A new extensive library of PHOENIX stellar atmospheres and synthetic spectra

Abstract: Aims. We present a new library of high-resolution synthetic spectra based on the stellar atmosphere code PHOENIX that can be used for a wide range of applications of spectral analysis and stellar parameter synthesis. Methods. The spherical mode of PHOENIX was used to create model atmospheres and to derive detailed synthetic stellar spectra from them. We present a new self-consistent way of describing micro-turbulence for our model atmospheres. Results. The synthetic spectra cover the wavelength range from 500 A to 5.5 μm with resolutions of R = 500 000 in the optical and near IR, R = 100 000 in the IR and Δλ = 0.1 A in the UV. The parameter space covers 2300 K ≤ Teff ≤ 12 000 K, 0.0 ≤ log g ≤ +6.0, −4.0 ≤ [Fe/H] ≤ +1.0, and −0.2 ≤ [α/Fe] ≤ +1.2. The library is a work in progress and we expect to extend it up to Teff = 25 000 K.

Mass assembly in quiescent and star-forming galaxies since z ≃ 4 from UltraVISTA

Abstract: We estimate the galaxy stellar mass function and stellar mass density for star-forming and quiescent galaxies with 0.2 1.5 consistent with the expected uncertainties. We also develop a new method to infer the specific star formation rate from the mass function of star-forming galaxies. We find that the specific star formation rate of 1010 − 10.5ℳ⊙ galaxies increases continuously in the redshift range 1 < z < 4. Finally, we compare our results with a semi-analytical model and find that these models overestimate the density of low mass quiescent galaxies by an order of magnitude, while the density of low-mass star-forming galaxies is successfully reproduced.

The UMIST database for astrochemistry 2012

Abstract: We present the fifth release of the UMIST Database for Astrochemistry (UDfA). The new reaction network contains 6173 gas-phase reactions, involving 467 species, 47 of which are new to this release. We have updated rate coefficients across all reaction types. We have included 1171 new anion reactions and updated and reviewed all photorates. In addition to the usual reaction network, we also now include, for download, state-specific deuterated rate coefficients, deuterium exchange reactions and a list of surface binding energies for many neutral species. Where possible, we have referenced the original source of all new and existing data. We have tested the main reaction network using a dark cloud model and a carbon-rich circumstellar envelope model. We present and briefly discuss the results of these models.

The HARPS search for southern extra-solar planets - XXXI. The M-dwarf sample

Abstract: Searching for planets around stars with different masses helps us to assess the outcome of planetary formation for different initial conditions. The low-mass M dwarfs are also the most frequent stars in our Galaxy and potentially therefore, the most frequent planet hosts. Aims. We present observations of 102 southern nearby M dwarfs, using a fraction of our guaranteed time on the ESO/HARPS spectrograph. We observed for 460 h and gathered 1965 precise (~1-3 m/s) radial velocities (RVs), spanning the period from Feb. 11, 2003 to Apr. 1, 2009. Methods. For each star observed, we derive a time series and its precision as well as its variability. We apply systematic searches for long-term trends, periodic signals, and Keplerian orbits (from one to four planets). We analyze the subset of stars with detected signals and apply several diagnostics to discriminate whether the observed Doppler shifts are caused by either stellar surface inhomogeneities or the radial pull of orbiting planets. To prepare for the statistical view of our survey, we also compute the limits on possible unseen signals, and derive a first estimate of the frequency of planets orbiting M dwarfs. Results. We recover the planetary signals of 9 planets announced by our group (Gl 176 b, Gl 581 b, c, d & e, Gl 674 b, Gl 433 b, Gl 667C b, and Gl 667C c). We present radial velocities confirming that GJ 849 hosts a Jupiter-mass planet, plus a long-term radial-velocity variation. We also present RVs that precise the planetary mass and period of Gl 832b. We detect long-term RV changes for Gl 367, Gl 680, and Gl 880, which are indicative of yet unknown long-period companions. We identify candidate signals in the radial-velocity time series of 11 other M dwarfs. Spectral diagnostics and/or photometric observations demonstrate however that these signals are most probably caused by stellar surface inhomogeneities. Finally, we find that our survey is sensitive to a few Earth-mass planets for periods up to several hundred days. We derive a first estimate of the occurrence of M-dwarf planets as a function of their minimum mass and orbital period. In particular, we find that giant planets (msini = 100 − 1000 M⊕) have a low frequency (e.g. f ≲ 1% for P = 1 − 10 d and f = 0.02+0.03-0.01 for P = 10 − 100 d), whereas super-Earths (msini = 1 − 10 M⊕) are likely very abundant (f = 0.36+0.25-0.10 for P = 1 − 10 d and f = 0.52+0.50-0.16 for P = 10 − 100 d). We also obtained η⊕ = 0.41+0.54-0.13, which is the frequency of habitable planets orbiting M dwarfs (1 ≤ msini ≤ 10 M⊕). For the first time, η⊕ is a direct measure and not a number extrapolated from the statistics of more massive and/or shorter-period planets.

Global survey of star clusters in the Milky Way II. The catalogue of basic parameters

Abstract: Context. Although they are the main constituents of the Galactic disk population, for half of the open clusters in the Milky Way reported in the literature nothing is known except the raw position and an approximate size.Aims. The main goal of this study is to determine a full set of uniform spatial, structural, kinematic, and astrophysical parameters for as many known open clusters as possible.Methods. On the basis of stellar data from PPMXL and 2MASS, we used a dedicated data-processing pipeline to determine kinematic and photometric membership probabilities for stars in a cluster region.Results. For an input list of 3784 targets from the literature, we confirm that 3006 are real objects, the vast majority of them are open clusters, but associations and globular clusters are also present. For each confirmed object we determined the exact position of the cluster centre, the apparent size, proper motion, distance, colour excess, and age. For about 1500 clusters, these basic astrophysical parameters have been determined for the first time. For the bulk of the clusters we also derived the tidal radius. We estimated additionally average radial velocities for more than 30% of the confirmed clusters. The present sample (called MWSC) reaches both the central parts of the Milky Way and its outer regions. It is almost complete up to 1.8 kpc from the Sun and also covers neighbouring spiral arms. However, for a small subset of the oldest open clusters (log t ≳ 9) we found some evidence of incompleteness within about 1 kpc from the Sun.

Automated data reduction workflows for astronomy - The ESO Reflex environment

Abstract: Context. Data from complex modern astronomical instruments often consist of a large number of di erent science and calibration files, and their reduction requires a variety of software tools. The execution chain of the tools represents a complex workflow that needs to be tuned and supervised, often by individual researchers that are not necessarily experts for any specific instrument. Aims. The e ciency of data reduction can be improved by using automatic workflows to organise data and execute a sequence of data reduction steps. To realize such e ciency gains, we designed a system that allows intuitive representation, execution and modification of the data reduction workflow, and has facilities for inspection and interaction with the data. Methods. The European Southern Observatory (ESO) has developed Reflex, an environment to automate data reduction workflows. Reflex is implemented as a package of customized components for the Kepler workflow engine. Kepler provides the graphical user interface to create an executable flowchart-like representation of the data reduction process. Key features of Reflex are a rule-based data organiser, infrastructure to re-use results, thorough book-keeping, data progeny tracking, interactive user interfaces, and a novel concept to exploit information created during data organisation for the workflow execution. Results. Automated workflows can greatly increase the e ciency of astronomical data reduction. In Reflex, workflows can be run noninteractively as a first step. Subsequent optimization can then be carried out while transparently re-using all unchanged intermediate products. We found that such workflows enable the reduction of complex data by non-expert users and minimizes mistakes due to book-keeping errors. Conclusions. Reflex includes novel concepts to increase the e ciency of astronomical data processing. While Reflex is a specific implementation of astronomical scientific workflows within the Kepler workflow engine, the overall design choices and methods can also be applied to other environments for running automated science workflows.

Herschel view of the Taurus B211/3 filament and striations: evidence of filamentary growth?

Abstract: We present first results from the Herschel Gould Belt survey for the B211/L1495 region in the Taurus molecular cloud. Thanks to their high sensitivity and dynamic range, the Herschel images reveal the structure of the dense, star-forming filament B211 with unprecedented detail, along with the presence of striations perpendicular to the filament and generally oriented along the magnetic field direction as traced by optical polarization vectors. Based on the column density and dust temperature maps derived from the Herschel data, we find that the radial density profile of the B211 filament approaches power-law behavior, ρ ∝ r−2.0± 0.4, at large radii and that the temperature profile exhibits a marked drop at small radii. The observed density and temperature profiles of the B211 filament are in good agreement with a theoretical model of a cylindrical filament undergoing gravitational contraction with a polytropic equation of state: P ∝ ργ and T ∝ ργ−1, with γ = 0.97 ± 0.01 < 1 (i.e., not strictly isothermal). The morphology of the column density map, where some of the perpendicular striations are apparently connected to the B211 filament, further suggests that the material may be accreting along the striations onto the main filament. The typical velocities expected for the infalling material in this picture are ~0.5–1 km s-1, which are consistent with the existing kinematical constraints from previous CO observations.

Baryon acoustic oscillations in the Lyα forest of BOSS quasars

Abstract: We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the Lyα forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range 2.1≤ z≤ 3.5 from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift z = 2.3, we measure the monopole and quadrupole components of the correlation function for separations in the range 20h −1 Mpc < r < 200h −1 Mpc. A peak in the correlation function is seen at a separation equal to (1.01± 0.03) times the distance 

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 a sample of solar neighborhood stars that have high-quality abundance determinations and show that there are two distinct regimes of [α/Fe] versus age, which we identify as the epochs of the thick and thin disk formation. A tight correlation between metallicity and [α/Fe] versus age is clearly identifiable for thick disk stars, implying that this population formed from a well mixed interstellar medium, probably initially in starburst and then more quiescently, over a time scale of 4−5 Gyr. Thick disk stars have vertical velocity dispersions which correlate with age, with the youngest objects of this population having small scale heights similar to 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 youngest thick disk set the initial conditions from which the inner thin disk started to form about 8 Gyr ago, at [Fe/H] in the range of (−0.1, +0.1) dex and [α/Fe] ∼ 0.1 dex. This also provides an explanation for the apparent coincidence between the existence of a step in metallicity at 7−10 kpc in the thin disk and the confinement of the thick disk within R < 10 kpc. We suggest that the outer thin disk developed outside the influence of the thick disk, giving rise to a separate structure, but also that the high alphaenrichment of those regions may originate from a primordial pollution of the outer regions by the gas expelled from the forming thick disk. Metal-poor thin disk stars ([Fe/H] < −0.4 dex) in the solar vicinity, whose properties are best explained by them originating in the outer disk, are shown to be as old as the youngest thick disk (9−10 Gyr). This implies that the outer thin disk started to form while the thick disk was still forming stars in the inner parts of the Galaxy. Hence, while the overall inner (thick+thin) disk is comprised of two structures with different scale lengths and whose combination may give the impression of an inside-out formation process, the thin disk itself probably formed its first stars in its outskirts. Moreover, we point out that, given the tight age−metallicity and age–[α/Fe] relations that exist in the thick disk, an inside-out process would give rise to a radial gradient in metallicity and α-elements in this population, which is not observed. Finally, we argue that our results leave little room for radial migration (in the sense of churning) either to have contaminated the solar vicinity, or, on a larger scale, to have redistributed stars in significant proportion across the solar annulus.

The deepest Herschel-PACS far-infrared survey: number counts and infrared luminosity functions from combined PEP/GOODS-H observations

Abstract: We present results from the deepest Herschel-Photodetector Array Camera and Spectrometer (PACS) far-infrared blank field extragalactic survey, obtained by combining observations of the Great Observatories Origins Deep Survey (GOODS) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and theconstruction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3σ depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 μm, respectively, and resolve ~75% of the cosmic infrared background at 100 μm and 160 μm into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths, and to determine the infrared luminosity function of galaxies down to L_(IR) = 10^(11) L⊙ at z ~ 1 and L_(IR) = 10^(12) L⊙ at z ~ 2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3σ) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

The O3N2 and N2 abundance indicators revisited: improved calibrations based on CALIFA and T e-based literature data

Abstract: The use of integral field spectroscopy is since recently allowing to measure the emission line fluxes of an increasingly large number of star-forming galaxies, both locally and at high redshift. Many studies have used these fluxes to derive the gas-phase metallicity of the galaxies by applying the so-called strong-line methods. However, the metallicity indicators that these datasets use were empirically calibrated using few direct abundance data points (T_e-based measurements). Furthermore, a precise determination of the prediction intervals of these indicators is commonly lacking in these calibrations. Such limitations might lead to systematic errors in determining the gas-phase metallicity, especially at high redshift, which might have a strong impact on our understanding of the chemical evolution of the Universe. The main goal of this study is to review the most widely used empirical oxygen calibrations, O3N2 and N2, by using new direct abundance measurements. We pay special attention to (1) the expected uncertainty of these calibrations as a function of the index value or abundance derived and (2) the presence of possible systematic offsets. This is possible thanks to the analysis of the most ambitious compilation of T_e-based H II regions to date. This new dataset compiles the Te-based abundances of 603 H II regions extracted from the literature but also includes new measurements from the CALIFA survey. Besides providing new and improved empirical calibrations for the gas abundance, we also present a comparison between our revisited calibrations with a total of 3423 additional CALIFA H II complexes with abundances derived using the ONS calibration from the literature. The combined analysis of T_e-based and ONS abundances allows us to derive their most accurate calibration to date for both the O3N2 and N2 single-ratio indicators, in terms of all statistical significance, quality, and coverage of the parameters space. In particular, we infer that these indicators show shallower abundance dependencies and statistically significant offsets compared to others'. The O3N2 and N2 indicators can be empirically applied to derive oxygen abundances calibrations from either direct abundance determinations with random errors of 0.18 and 0.16, respectively, or from indirect ones (but based on a large amount of data), reaching an average precision of 0.08 and 0.09 dex (random) and 0.02 and 0.08 dex (systematic; compared to the direct estimations), respectively.

Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars - V. Evidence for a wide age distribution and a complex MDF

Abstract: Based on high-resolution spectra obtained during gravitational microlensing events we present a detailed elemental abundance analysis of 32 dwarf and subgiant stars in the Galactic bulge. Combined with the sample of 26 stars from the previous papers in this series, we now have 58 microlensed bulge dwarfs and subgiants that have been homogeneously analysed. The main characteristics of the sample and the findings that can be drawn are: (i) the metallicity distribution (MDF) is wide and spans all metallicities between [Fe/H] = −1.9 to +0.6; (ii) the dip in the MDF around solar metallicity that was apparent in our previous analysis of a smaller sample (26 microlensed stars) is no longer evident; instead it has a complex structure and indications of multiple components are starting to emerge. A tentative interpretation is that there could be different stellar populations at interplay, each with a different scale height: the thin disk, the thick disk, and a bar population; (iii) the stars with [Fe/H] ≲ −0.1 are old with ages between 10 and 12 Gyr; (iv) the metal-rich stars with [Fe/H] ≳ −0.1 show a wide variety of ages, ranging from 2 to 12 Gyr with a distribution that has a dominant peak around 4−5 Gyr and a tail towards higher ages; (v) there are indications in the [α/Fe]−[Fe/H] abundance trends that the “knee” occurs around [Fe/H] = −0.3 to −0.2, which is a slightly higher metallicity as compared to the “knee” for the local thick disk. This suggests that the chemical enrichment of the metal-poor bulge has been somewhat faster than what is observed for the local thick disk. The results from the microlensed bulge dwarf stars in combination with other findings in the literature, in particular the evidence that the bulge has cylindrical rotation, indicate that the Milky Way could be an almost pure disk galaxy. The bulge would then just be a conglomerate of the other Galactic stellar populations (thin disk, thick disk, halo, and ...?), residing together in the central parts of the Galaxy, influenced by the Galactic bar.

The VLT-FLAMES Tarantula Survey. VIII. Multiplicity properties of the O-type star population

Abstract: Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^\pi\ (with 0.15 7.8', i.e. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars.

Chemodynamical evolution of the Milky Way disk - I. The solar vicinity

Abstract: In the first paper of this series, we present a new approach for studying the chemo-dynamical evolution in disk galaxies, which consists of fusing disk chemical evolution models with compatible numerical simulations of galactic disks. This method avoids known star formation and chemical enrichment problems encountered in simulations. Here we focus on the Milky Way, by using a detailed thin-disk chemical evolution model (matching local observables, which are weakly affected by radial migration) and a simulation in the cosmological context, with dynamical properties close to those of our Galaxy. We examine in detail the interplay between in situ chemical enrichment and radial migration and their impact on key observables in the solar neighborhood, e.g., the age-metallicity-velocity relation, the metallicity distribution, and gradients in the radial and vertical directions. We show that, due to radial migration from mergers at high redshift and the central bar at later times, a sizable fraction of old metal-poor high-[α /Fe] stars can reach the solar vicinity. This naturally accounts for a number of recent observations related to both the thin and thick disks, despite the fact that we use thin-disk chemistry only. Although significant radial mixing is present, the slope in the age-metallicity relation is only weakly affected, with a scatter compatible with recent observational work. While we find a smooth density distribution in the [O/Fe]–[Fe/H] plane, we can recover the observed discontinuity by selecting particles according to kinematic criteria used in high-resolution samples to define the thin and thick disks. We outline a new method for estimating the birth place of the Sun and predict that the most likely radius lies in the range 4.4 = 8 kpc). A new, unifying model for the Milky Way thick disk is offered, where both mergers and radial migration play a role at different stages of the disk evolution. We show that in the absence of early-on massive mergers the vertical velocity dispersion of the oldest stars is underestimated by a factor of ~2 compared with observations. We can, therefore, argue that the Milky Way thick disk is unlikely to have been formed through a quiescent disk evolution. An observational test involving both chemical and kinematic information must be devised to ascertain this possibility.

Cores, filaments, and bundles: hierarchical core formation in the L1495/B213 Taurus region

Abstract: Context. Core condensation is a critical step in the star-formation process, but it is still poorly characterized observationally.Aims. We have studied the 10 pc-long L1495/B213 complex in Taurus to investigate how dense cores have condensed out of the lower density cloud material.Methods. We observed L1495/B213 in C18 O(1−0), N2 H+ (1−0), and SO(J N = 32 –21 ) with the 14 m FCRAO telescope, and complemented the data with dust continuum observations using APEX (870 μ m) and IRAM 30 m (1200 μ m).Results. From the N2 H+ emission, we identify 19 dense cores, some starless and some protostellar. They are not distributed uniformly, but tend to cluster with relative separations on the order of 0.25 pc. From the C18 O emission, we identify multiple velocity components in the gas. We have characterized them by fitting Gaussians to the spectra and by studying the distribution of the fits in position–position–velocity space. In this space, the C18 O components appear as velocity-coherent structures, and we identify them automatically using a dedicated algorithm (FIVE: Friends In VElocity). Using this algorithm, we identify 35 filamentary components with typical lengths of 0.5 pc, sonic internal velocity dispersions, and mass-per-unit length close to the stability threshold of isothermal cylinders at 10 K. Core formation seems to have occurred inside the filamentary components via fragmentation, with few fertile components with higher mass-per-unit length being responsible for most cores in the cloud. On large scales, the filamentary components appear grouped into families, which we refer to as bundles.Conclusions. Core formation in L1495/B213 has proceeded by hierarchical fragmentation. The cloud fragmented first into several pc-scale regions. Each of these regions later fragmented into velocity-coherent filaments of about 0.5 pc in length. Finally, a small number of these filaments fragmented quasi-statically and produced the individual dense cores we see today.

The evolution of amorphous hydrocarbons in the ISM: dust modelling from a new vantage point

Abstract: Context. The evolution of amorphous hydrocarbon materials, a-C(:H), principally resulting from ultraviolet (UV) photon absorption-induced processing, are likely at the heart of the variations in the observed properties of dust in the interstellar medium.Aims. The consequences of the size-dependent and compositional variations in a-C(:H), from aliphatic-rich a-C:H to aromatic-rich a-C, are studied within the context of the interstellar dust extinction and emission.Methods. Newly-derived optical property data for a-C(:H) materials, combined with that for an amorphous forsterite-type silicate with iron nano-particle inclusions, a-SilFe , are used to explore dust evolution in the interstellar medium.Results. We present a new dust model that consists of a power-law distribution of small a-C grains and log-normal distributions of large a-SilFe and a-C(:H) grains. The model, which is firmly anchored by laboratory-data, is shown to quite naturally explain the variations in the infrared (IR) to far-ultraviolet (FUV) extinction, the 217 nm UV bump, the IR absorption and emission bands and the IR-mm dust emission.Conclusions. The major strengths of the new model are its inherent simplicity and built-in capacity to follow dust evolution in interstellar media. We show that mantle accretion in molecular clouds and UV photo-processing in photo-dominated regions are likely the major drivers of dust evolution.

Rotation and differential rotation of active Kepler stars

Abstract: Context. The Kepler space telescope monitors more than 160 000 stars with an unprecedented precision providing the opportunity to study the rotation of thousands of stars. Aims. We present rotation periods for thousands of active stars in the Kepler field derived from Q3 data. In most cases a second period close to the rotation period was detected that we interpreted as surface differential rotation (DR). We show how the absolute and relative shear (ΔΩ and α = ΔΩ/Ω, respectively) correlate with rotation period and effective temperature.Methods. Active stars were selected from the whole sample using the range of the variability amplitude. To detect different periods in the light curves we used the Lomb-Scargle periodogram in a pre-whitening approach to achieve parameters for a global sine fit. The most dominant periods from the fit were associated to different surface rotation periods. Our purely mathematical approach is capable of detecting different periods but cannot distinguish between the physical origins of periodicity. We ascribe the existence of different periods to DR, but spot evolution could also play a role. Because of the large number of stars the period errors are estimated statistically. We thus cannot exclude the existence of false positives among our periods. Results. In our sample of 40 661 active stars we found 24 124 rotation periods P 1 between 0.5 and 45 days, with a mean of ⟨P 1 ⟩ = 16.3 days. The distribution of stars with 0.5 − V 2 within ±30% of the rotation period P 1 was found in 18 616 stars (77.2%). Attributing these two periods to DR we found that for active stars other than the Sun the relative shear α increases with rotation period, and slightly decreases with effective temperature. The absolute shear ΔΩ slightly increases from ΔΩ = 0.079 rad d-1 at T eff = 3500 K to ΔΩ = 0.096 rad d-1 at T eff = 6000 K. Above 6000 K, ΔΩ shows much larger scatter. The dependence of ΔΩ on rotation period is weak over a large period range. Conclusions. Latitudinal differential rotation measured for the first time in more than 18 000 stars provides a comprehensive picture of stellar surface shear. This picture is consistent with major predictions from mean-field theory, and seems to support these models. To what extent our observations are prone to false positives and selection bias has not been fully explored, and needs to be addressed using other data, including the full Kepler time coverage.

Improved angular momentum evolution model for solar-like stars

Abstract: Context. Understanding the origin and evolution of stellar angular momentum is one of the major challenges of stellar physics.Aims. We present new models for the rotational evolution of solar-like stars between 1 Myr and 10 Gyr with the aim of reproducing the distributions of rotational periods observed for star forming regions and young open clusters within this age range.Methods. The models include a new wind braking law based on recent numerical simulations of magnetized stellar winds and specific dynamo and mass-loss prescriptions are adopted to tie angular momentum loss to angular velocity. The models additionally assume constant angular velocity during the disk accretion phase and allow for decoupling between the radiative core and the convective envelope as soon as the former develops.Results. We have developed rotational evolution models for slow, median, and fast rotators with initial periods of 10, 7, and 1.4 d, respectively. The models reproduce reasonably well the rotational behavior of solar-type stars between 1 Myr and 4.5 Gyr, including pre-main sequence (PMS) to zero-age main sequence (ZAMS) spin up, prompt ZAMS spin down, and the early-main sequence (MS) convergence of surface rotation rates. We find the model parameters accounting for the slow and median rotators are very similar to each other, with a disk lifetime of 5 Myr and a core-envelope coupling timescale of 28−30 Myr. In contrast, fast rotators have both shorter disk lifetimes (2.5 Myr) and core-envelope coupling timescales (12 Myr). We show that a large amount of angular momentum is hidden in the radiative core for as long as 1 Gyr in these models and we discuss the implications for internal differential rotation and lithium depletion. We emphasize that these results are highly dependent on the adopted braking law. We also report a tentative correlation between the initial rotational period and disk lifetime, which suggests that protostellar spin down by massive disks in the embedded phase is at the origin of the initial dispersion of rotation rates in young stars.Conclusions. We conclude that this class of semi-empirical models successfully grasp the main trends of the rotational behavior of solar-type stars as they evolve and make specific predictions that may serve as a guide for further development.

The VIMOS VLT Deep Survey final data release: a spectroscopic sample of 35 016 galaxies and AGN out to z ~ 6.7 selected with 17.5 ≤ iAB ≤ 24.75

Abstract: We describe the completed VIMOS VLT Deep Survey, and the final data release of 35016 galaxies and type-I AGN with measured spectroscopic redshifts up to redshift z~6.7, in areas 0.142 to 8.7 square degrees, and volumes from 0.5x10^6 to 2x10^7h^-3Mpc^3. We have selected samples of galaxies based solely on their i-band magnitude reaching i_{AB}=24.75. Spectra have been obtained with VIMOS on the ESO-VLT, integrating 0.75h, 4.5h and 18h for the Wide, Deep, and Ultra-Deep nested surveys. A total of 1263 galaxies have been re-observed independently within the VVDS, and from the VIPERS and MASSIV surveys. They are used to establish the redshift measurements reliability, to assess completeness, and to provide a weighting scheme taking into account the survey selection function. We describe the main properties of the VVDS samples, and the VVDS is compared to other spectroscopic surveys. In total we have obtained spectroscopic redshifts for 34594 galaxies, 422 type-I AGN, and 12430 Galactic stars. The survey has enabled to identify galaxies up to very high redshifts with 4669 redshifts in 1 3, and specific populations like LAE have been identified out to z=6.62. We show that the VVDS occupies a unique place in the parameter space defined by area, depth, redshift coverage, and number of spectra. The VVDS provides a comprehensive survey of the distant universe, covering all epochs since z, or more than 12 Gyr of cosmic time, with a uniform selection, the largest such sample to date. A wealth of science results derived from the VVDS have shed new light on the evolution of galaxies and AGN, and their distribution in space, over this large cosmic time. A final public release of the complete VVDS spectroscopic redshift sample is available at http://cesam.lam.fr/vvds.

Complex organic molecules in the interstellar medium: IRAM 30 m line survey of Sagittarius B2(N) and (M)

Abstract: Context. The discovery of amino acids in meteorites fallen to Earth and the detection of glycine, the simplest of them, in samples returned from a comet to Earth strongly suggest that the chemistry of the interstellar medium is capable of producing such complex organic molecules and that they may be widespread in our Galaxy. Aims. Our goal is to investigate the degree of chemical complexity that can be reached in the interstellar medium, in particular in dense star-forming regions. Methods. We performed an unbiased, spectral line survey toward Sgr B2(N) and (M), two regions where high-mass stars are formed, with the IRAM 30 m telescope in the 3 mm atmospheric transmission window. Partial surveys at 2 and 1.3 mm were performed in parallel. The spectra were analyzed with a simple radiative transfer model that assumes local thermodynamic equilibrium but takes optical depth effects into account. Results. About 3675 and 945 spectral lines with a peak signal-to-noise ratio higher than 4 are detected at 3 mm toward Sgr B2(N) and (M), i.e. about 102 and 26 lines per GHz, respectively. This represents an increase by about a factor of two over previous surveys of Sgr B2. About 70% and 47% of the lines detected toward Sgr B2(N) and (M) are identified and assigned to 56 and 46 distinct molecules as well as to 66 and 54 less abundant isotopologues of these molecules, respectively. In addition, we report the detection of transitions from 59 and 24 catalog entries corresponding to vibrationally or torsionally excited states of some of these molecules, respectively, up to a vibration energy of 1400 cm-1 (2000 K). Excitation temperatures and column densities were derived for each species but should be used with caution. The rotation temperatures of the detected complex molecules typically range from ~50 to 200 K. Among the detected molecules, aminoacetonitrile, n- propyl cyanide, and ethyl formate were reported for the first time in space based on this survey, as were five rare isotopologues of vinyl cyanide, cyanoacetylene, and hydrogen cyanide. We also report the detection of transitions from within twelve new vibrationally or torsionally excited states of known molecules. Absorption features produced by diffuse clouds along the line of sight are detected in transitions with low rotation quantum numbers of many simple molecules and are modeled with ~30–40 velocity components with typical linewidths of ~3–5 km s-1 . Conclusions. Although the large number of unidentified lines may still allow future identification of new molecules, we expect most of these lines to belong to vibrationally or torsionally excited states or to rare isotopologues of known molecules for which spectroscopic predictions are currently missing. Significant progress in extending the inventory of complex organic molecules in Sgr B2(N) and deriving tighter constraints on their location, origin, and abundance is expected in the near future thanks to an ongoing spectral line survey at 3 mm with ALMA in its cycles 0 and 1. The present single-dish survey will serve as a solid basis for the line identification and analysis of such an interferometric survey.

Planck intermediate results: V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect

Abstract: Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from xmm in the region in which the profiles overlap (i.e., [0.1-1] R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is shallower than the predictions. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0, c500, gamma, alpha, beta] = [6.41, 1.81, 0.31, 1.33, 4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

Global collapse of molecular clouds as a formation mechanism for the most massive stars

Abstract: The relative importance of primordial molecular cloud fragmentation versus large-scale accretion still remains to be assessed in the context of massive core/star formation. Studying the kinematics of the dense gas surrounding massive-star progenitors can tell us the extent to which large-scale flow of material impacts the growth in mass of star-forming cores. Here we present a comprehensive dataset of the 5500(±800) M⊙ infrared dark cloud SDC335.579-0.272 (hereafter SDC335), which exhibits a network of cold, dense, parsec-long filaments. Atacama Large Millimeter Array (ALMA) Cycle 0 observations reveal two massive star-forming cores, MM1 and MM2, sitting at the centre of SDC335 where the filaments intersect. With a gas mass of 545(-385+770) M⊙ contained within a source diameter of 0.05 pc, MM1 is one of the most massive, compact protostellar cores ever observed in the Galaxy. As a whole, SDC335 could potentially form an OB cluster similar to the Trapezium cluster in Orion. ALMA and Mopra single-dish observations of the SDC335 dense gas furthermore reveal that the kinematics of this hub-filament system are consistent with a global collapse of the cloud. These molecular-line data point towards an infall velocity Vinf = 0.7( ± 0.2) km s-1, and a total mass infall rate Ṁinf ≃ 2.5(±1.0) × 10-3 M⊙ yr-1 towards the central pc-size region of SDC335. This infall rate brings 750(±300) M⊙ of gas to the centre of the cloud per free-fall time (tff = 3 × 105 yr). This is enough to double the mass already present in the central pc-size region in 3.5-1.0+2.2 × tff. These values suggest that the global collapse of SDC335 over the past million year resulted in the formation of an early O-type star progenitor at the centre of the cloud’s gravitational potential well.

A Herschel [C ii] Galactic plane survey - I. The global distribution of ISM gas components

Abstract: Context. The [C ii] 158 μ m 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 (ISM), including the diffuse ionized medium, warm and cold atomic clouds, clouds in transition from atomic to molecular, and dense and warm photon dominated regions.Aims. Velocity-resolved observations of [C ii] are the most powerful technique available to disentangle the emission produced by these components. These observations can also be used to trace CO-dark H2 gas and determine the total mass of the ISM. Methods. The Galactic Observations of Terahertz C+ (GOT C+) project surveys the [C ii] 158 μ m line over the entire Galactic disk with velocity-resolved observations using the Herschel /HIFI instrument. We present the first longitude-velocity maps of the [C ii] emission for Galactic latitudes b = 0°, ±0.5°, and ±1.0°. We combine these maps with those of H i, 12 CO, and 13 CO to separate the different phases of the ISM and study their properties and distribution in the Galactic plane. Results. [C ii] emission is mostly associated with spiral arms, mainly emerging from Galactocentric distances between 4 and 10 kpc. It traces the envelopes of evolved clouds as well as clouds that are in the transition between atomic and molecular. We estimate that most of the observed [C ii] emission is produced by dense photon dominated regions (~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 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 12 CO and 13 CO (4−8 kpc). 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 [C ii] emission emerging from dense photon dominated regions is associated with modest far-ultraviolet fields in the range χ 0 ≃ 1 − 30.

The VIMOS Public Extragalactic Redshift Survey (VIPERS) - Galaxy clustering and redshift-space distortions at z ≃ 0.8 in the first data release

Abstract: We present in this paper the general real- and redshift-space clustering properties of galaxies as measured in the first data release of the VIPERS survey. VIPERS is a large redshift survey designed to probe the distant Universe and its large-scale structure at 0.5 < z < 1.2. We describe in this analysis the global properties of the sample and discuss the survey completeness and associated corrections. This sample allows us to measure the galaxy clustering with an unprecedented accuracy at these redshifts. From the redshift-space distortions observed in the galaxy clustering pattern we provide a first measurement of the growth rate of structure at z = 0.8: f\sigma_8 = 0.47 +/- 0.08. This is completely consistent with the predictions of standard cosmological models based on Einstein gravity, although this measurement alone does not discriminate between different gravity models.

3D mixing in hot Jupiters atmospheres. I. Application to the day/night cold trap in HD 209458b

Abstract: Context. Hot Jupiters exhibit atmospheric temperatures ranging from hundreds to thousands of Kelvin. Because of their large day-night temperature differences, condensable species that are stable in the gas phase on the dayside – such as TiO and silicates – may condense and gravitationally settle on the nightside. Atmospheric circulation may counterbalance this tendency to gravitationally settle. This three-dimensional (3D) mixing of condensable species has not previously been studied for hot Jupiters, yet it is crucial to assess the existence and distribution of TiO and silicates in the atmospheres of these planets.Aims. We investigate the strength of the nightside cold trap in hot Jupiters atmospheres by investigating the mechanisms and strength of the vertical mixing in these stably stratified atmospheres. We apply our model to the particular case of TiO to address the question of whether TiO can exist at low pressure in sufficient abundances to produce stratospheric thermal inversions despite the nightside cold trap.Methods. We modeled the 3D circulation of HD 209458b including passive (i.e. radiatively inactive) tracers that advect with the 3D flow, with a source and sink term on the nightside to represent their condensation into haze particles and their gravitational settling.Results. We show that global advection patterns produce strong vertical mixing that can keep condensable species aloft as long as they are trapped in particles of sizes of a few microns or less on the nightside. We show that vertical mixing results not from small-scale convection but from the large-scale circulation driven by the day-night heating contrast. Although this vertical mixing is not diffusive in any rigorous sense, a comparison of our results with idealized diffusion models allows a rough estimate of the effective vertical eddy diffusivities in these atmospheres. The parametrization Kzz =5 × 104 / P bar m2 s-1 , valid from ~1 bar to a few μ bar, can be used in 1D models of HD 209458b. Moreover, our models exhibit strong spatial and temporal variability in the tracer concentration that could result in observable variations during either transit or secondary eclipse measurements. Finally, we apply our model to the case of TiO in HD 209458b and show that the day-night cold trap would deplete TiO if it condenses into particles bigger than a few microns on the planet’s nightside, keeping it from creating the observed stratosphere of the planet.

Seismic diagnostics for transport of angular momentum in stars. I. Rotational splittings from the pre-main sequence to the red-giant branch.

Abstract: Context. Rotational splittings are currently measured for several main sequence stars and a large number of red giants with the space mission Kepler. This will provide stringent constraints on rotation profiles. Aims: Our aim is to obtain seismic constraints on the internal transport and surface loss of the angular momentum of oscillating solar-like stars. To this end, we study the evolution of rotational splittings from the pre-main sequence to the red-giant branch for stochastically excited oscillation modes. Methods: We modified the evolutionary code CESAM2K to take rotationally induced transport in radiative zones into account. Linear rotational splittings were computed for a sequence of 1.3 Ms models. Rotation profiles were derived from our evolutionary models and eigenfunctions from linear adiabatic oscillation calculations. Results: We find that transport by meridional circulation and shear turbulence yields far too high a core rotation rate for red-giant models compared with recent seismic observations. We discuss several uncertainties in the physical description of stars that could have an impact on the rotation profiles. For instance, we find that the Goldreich-Schubert-Fricke instability does not extract enough angular momentum from the core to account for the discrepancy. In contrast, an increase of the horizontal turbulent viscosity by 2 orders of magnitude is able to significantly decrease the central rotation rate on the red-giant branch. Conclusions: Our results indicate that it is possible that the prescription for the horizontal turbulent viscosity largely underestimates its actual value or else a mechanism not included in current stellar models of low mass stars is needed to slow down the rotation in the radiative core of red-giant stars.

Ice condensation as a planet formation mechanism

Abstract: We show that condensation is an efficient particle growth mechanism that leads to growth beyond decimetre-sized pebbles close to an ice line in protoplanetary discs. As coagulation of dust particles is frustrated by bouncing and fragmentation, condensation could be a complementary, or even dominant, growth mode in the early stages of planet formation. Ice particles diffuse across the ice line and sublimate, and vapour diffusing inwards across the ice line recondenses onto already existing particles, causing them to grow. We develop a numerical model of the dynamical behaviour of ice particles close to the water ice line, approximately 3 AU from the host star. Particles move with the turbulent gas, modelled as a random walk. They also sediment towards the midplane and drift radially towards the central star. Condensation and sublimation are calculated using a Monte Carlo approach. Our results indicate that, with a turbulent alpha-value of 0.01, growth from millimetre to at least decimetre-sized pebbles is possible on a time scale of 1000 years. We find that particle growth is dominated by ice and vapour transport across the radial ice line, with negligible growth caused by transport across the atmospheric ice line. Ice particles mix outwards by turbulent diffusion, leading to net growth across the entire cold region. The resulting particles are large enough to be sensitive to concentration by streaming instabilities, pressure bumps and vortices, which can cause further growth into planetesimals. In our model, particles are considered to be homogeneous ice particles. Considering the more realistic composition of ice condensed onto rocky ice nuclei might affect the growth time scales, by release of refractory ice nuclei after sublimation. We also ignore sticking and fragmentation in particle collisions. These effects will be the subject of future investigations. (Less)

Populations of rotating stars - I. Models from 1.7 to 15 M⊙ at Z = 0.014, 0.006, and 0.002 with Ω/Ωcrit between 0 and 1

Abstract: Context. B-type stars are known to rotate at various velocities, including very fast rotators near the critical velocity as the Be stars.Aims. In this paper, we provide stellar models covering the mass range between 1.7 to 15 M ⊙ , which includes the typical mass of known Be stars, at Z = 0.014, 0.006, and 0.002 and for an extended range of initial velocities on the zero-age main sequence.Methods. We used the Geneva stellar-evolution code, including the effects of shellular rotation, with a numerical treatment that has been improved so the code can precisely track the variation in the angular momentum content of the star as it changes under the influence of radiative winds and/or mechanical mass loss.Results. We discuss the impact of the initial rotation rate on the tracks in the Hertzsprung-Russell diagram, the main-sequence (MS) lifetimes, the evolution of the surface rotation and abundances, as well as on the ejected masses of various isotopes. Among the new results obtained from the present grid we find that 1) fast-rotating stars with initial masses around 1.7 M ⊙ present at the beginning of the core hydrogen-burning phase quite small convective cores with respect to their slowly rotating counterparts. This fact may be interesting to keep in mind in the framework of the asteroseismic studies of such stars. 2) The contrast between the core and surface angular velocity is higher in slower rotating stars. Our results are in agreement with the very few values obtained for B-type stars from asteroseismology. 3) At Z = 0.002, the stars in the mass range of 1.7 to 3 M ⊙ with a mean velocity on the MS of the order of 150 km s-1 show N/H enhancement superior to 0.2 dex at mid-MS, and superior to 0.4 dex at the end of the MS phase. At solar metallicity the corresponding values are below 0.2 dex at any time in the MS.Conclusions. An extended database of stellar models containing 270 evolutionary tracks is provided to the community.