Showing papers in "Astronomy and Astrophysics in 2012"

Grids of stellar models with rotation - I. Models from 0.8 to 120 M⊙ at solar metallicity (Z = 0.014)

Abstract: Aims. Many topical astrophysical research areas, such as the properties of planet host stars, the nature of the progenitors of different types of supernovae and gamma ray bursts, and the evolution of galaxies, require complete and homogeneous sets of stellar models at different metallicities in order to be studied during the whole of cosmic history. We present here a first set of models for solar metallicity, where the effects of rotation are accounted for in a homogeneous way.Methods. We computed a grid of 48 different stellar evolutionary tracks, both rotating and non-rotating, at Z = 0.014, spanning a wide mass range from 0.8 to 120 M ⊙ . For each of the stellar masses considered, electronic tables provide data for 400 stages along the evolutionary track and at each stage, a set of 43 physical data are given. These grids thus provide an extensive and detailed data basis for comparisons with the observations. The rotating models start on the zero-age main sequence (ZAMS) with a rotation rate υ ini /υ crit = 0.4. The evolution is computed until the end of the central carbon-burning phase, the early asymptotic giant branch (AGB) phase, or the core helium-flash for, respectively, the massive, intermediate, and both low and very low mass stars. The initial abundances are those deduced by Asplund and collaborators, which best fit the observed abundances of massive stars in the solar neighbourhood. We update both the opacities and nuclear reaction rates, and introduce new prescriptions for the mass-loss rates as stars approach the Eddington and/or the critical velocity. We account for both atomic diffusion and magnetic braking in our low-mass star models.Results. The present rotating models provide a good description of the average evolution of non-interacting stars. In particular, they reproduce the observed main-sequence width, the positions of the red giant and supergiant stars in the Hertzsprung-Russell (HR) diagram, the observed surface compositions and rotational velocities. Very interestingly, the enhancement of the mass loss during the red-supergiant stage, when the luminosity becomes supra-Eddington in some outer layers, help models above 15−20 M ⊙ to lose a significant part of their hydrogen envelope and evolve back into the blue part of the HR diagram. This result has interesting consequences for the blue to red supergiant ratio, the minimum mass for stars to become Wolf-Rayet stars, and the maximum initial mass of stars that explode as type II−P supernovae.

CALIFA, the Calar Alto Legacy Integral Field Area survey : I. Survey presentation

Abstract: The final product of galaxy evolution through cosmic time is the population of galaxies in the local universe. These galaxies are also those that can be studied in most detail, thus providing a stringent benchmark for our understanding of galaxy evolution. Through the huge success of spectroscopic single-fiber, statistical surveys of the Local Universe in the last decade, it has become clear, however, that an authoritative observational description of galaxies will involve measuring their spatially resolved properties over their full optical extent for a statistically significant sample. We present here the Calar Alto Legacy Integral Field Area (CALIFA) survey, which has been designed to provide a first step in this direction. We summarize the survey goals and design, including sample selection and observational strategy. We also showcase the data taken during the first observing runs (June/July 2010) and outline the reduction pipeline, quality control schemes and general characteristics of the reduced data. This survey is obtaining spatially resolved spectroscopic information of a diameter selected sample of similar to 600 galaxies in the Local Universe (0.005 < z < 0.03). CALIFA has been designed to allow the building of two-dimensional maps of the following quantities: (a) stellar populations: ages and metallicities; (b) ionized gas: distribution, excitation mechanism and chemical abundances; and (c) kinematic properties: both from stellar and ionized gas components. CALIFA uses the PPAK integral field unit (IFU), with a hexagonal field-of-view of similar to 1.3 square', with a 100% covering factor by adopting a three-pointing dithering scheme. The optical wavelength range is covered from 3700 to 7000 angstrom, using two overlapping setups (V500 and V1200), with different resolutions: R similar to 850 and R similar to 1650, respectively. CALIFA is a legacy survey, intended for the community. The reduced data will be released, once the quality has been guaranteed. The analyzed data fulfill the expectations of the original observing proposal, on the basis of a set of quality checks and exploratory analysis: (i) the final datacubes reach a 3 sigma limiting surface brightness depth of similar to 23.0 mag/arcsec(2) for the V500 grating data (similar to 22.8 mag/arcsec(2) for V1200); (ii) about similar to 70% of the covered field-of-view is above this 3 sigma limit; (iii) the data have a blue-to-red relative flux calibration within a few percent in most of the wavelength range; (iv) the absolute flux calibration is accurate within similar to 8% with respect to SDSS; (v) the measured spectral resolution is similar to 85 km s(-1) for V1200 (similar to 150 km s(-1) for V500); (vi) the estimated accuracy of the wavelength calibration is similar to 5 km s(-1) for the V1200 data (similar to 10 km s(-1) for the V500 data); (vii) the aperture matched CALIFA and SDSS spectra are qualitatively and quantitatively similar. Finally, we show that we are able to carry out all measurements indicated above, recovering the properties of the stellar populations, the ionized gas and the kinematics of both components. The associated maps illustrate the spatial variation of these parameters across the field, reemphasizing the redshift dependence of single aperture spectroscopic measurements. We conclude from this first look at the data that CALIFA will be an important resource for archaeological studies of galaxies in the Local Universe.

Rapid growth of gas-giant cores by pebble accretion

Abstract: The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before gas dissipation in the protoplanetary disc is completed within 1–10 million years. Building up the core by collisions between km-sized planetesimals fails to meet this timescale constraint, especially at wide stellar separations. Nonetheless, gas-giant planets are detected by direct imaging at wide orbital distances. In this paper, we numerically study the growth of cores by the accretion of cm-sized pebbles loosely coupled to the gas. We measure the accretion rate onto seed masses ranging from a large planetesimal to a fully grown 10-Earth-mass core and test different particle sizes. The numerical results are in good agreement with our analytic expressions, indicating the existence of two accretion regimes, one set by the azimuthal and radial particle drift for the lower seed masses and the other, for higher masses, by the velocity at the edge of the Hill sphere. In the former, the optimally accreted particle size increases with core mass, while in the latter the optimal size is centimeters, independent of core mass. We discuss the implications for rapid core growth of gas-giant and ice-giant cores. We conclude that pebble accretion can resolve the long-standing core accretion timescale conflict. This requires a near-unity dust-to-gas ratio in the midplane, particle growth to mm and cm and the formation of massive planetesimals or low radial pressure support. The core growth timescale is shortened by a factor 30–1000 at 5 AU and by a factor 100–10 000 at 50 AU, compared to the gravitationally focused accretion of, respectively, low-scale-height planetesimal fragments or standard km-sized planetesimals.

UltraVISTA: a new ultra-deep near-infrared survey in COSMOS

Abstract: In this paper we describe the first data release of the UltraVISTA near-infrared imaging survey of the COSMOS field. We summarise the key goals and design of the survey and provide a detailed description of our data reduction techniques. We provide stacked, sky-subtracted images in YJHK_s and narrow-band filters constructed from data collected during the first year of UltraVISTA observations. Our stacked images reach 5σAB depths in an aperture of 2″ diameter of ~25 in Y and ~24 in JHK_s bands and all have sub-arcsecond seeing. To this 5σ limit, our K_s catalogue contains 216 268 sources. We carry out a series of quality assessment tests on our images and catalogues, comparing our stacks with existing catalogues. The 1σ astrometric rms in both directions for stars selected with 17.0 < K_s(AB) < 19.5 is ~0.08″ in comparison to the publicly-available COSMOS ACS catalogues. Our images are resampled to the same pixel scale and tangent point as the publicly available COSMOS data and so may be easily used to generate multi-colour catalogues using this data. All images and catalogues presented in this paper are publicly available through ESO’s “phase 3” archiving and distribution system and from the UltraVISTA web site.

REBOUND: an open-source multi-purpose N-body code for collisional dynamics

Abstract: REBOUND is a new multi-purpose N -body code which is freely available under an open-source license. It was designed for collisional dynamics such as planetary rings but can also solve the classical N -body problem. It is highly modular and can be customized easily to work on a wide variety of different problems in astrophysics and beyond. REBOUND comes with three symplectic integrators: leap-frog, the symplectic epicycle integrator (SEI) and a Wisdom-Holman mapping (WH). It supports open, periodic and shearing-sheet boundary conditions. REBOUND can use a Barnes-Hut tree to calculate both self-gravity and collisions. These modules are fully parallelized with MPI as well as OpenMP. The former makes use of a static domain decomposition and a distributed essential tree. Two new collision detection modules based on a plane-sweep algorithm are also implemented. The performance of the plane-sweep algorithm is superior to a tree code for simulations in which one dimension is much longer than the other two and in simulations which are quasi-two dimensional with less than one million particles. In this work, we discuss the different algorithms implemented in REBOUND, the philosophy behind the code’s structure as well as implementation specific details of the different modules. We present results of accuracy and scaling tests which show that the code can run efficiently on both desktop machines and large computing clusters.

A simple model for the evolution of the dust population in protoplanetary disks

Abstract: Context. The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary disks and in the formation of larger bodies, such as planetesimals.Aims. We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the grain size distribution and which can readily be used for further modeling or for interpreting of observational data.Methods. We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmentation, and transport in viscously evolving gas disks.Results. We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the radial flux by solid material “rain out”, which is crucial for triggering any gravity assisted formation of planetesimals. We show that fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of -1.5, while the outer regions at later times can become drift-dominated, yielding a dust surface density exponent of −0.75. Our results show that radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation due to the turbulent state of the disk.

Chemical abundances of 1111 FGK stars from the HARPS GTO planet search program - Galactic stellar populations and planets

Abstract: Context. We performed a uniform and detailed abundance analysis of 12 refractory elements (Na, Mg, Al, Si, Ca, Ti, Cr, Ni, Co, Sc, Mn, and V) for a sample of 1111 FGK dwarf stars from the HARPS GTO planet search program. Of these stars, 109 are known to harbor giant planetary companions and 26 stars are exclusively hosting Neptunians and super-Earths.Aims. The two main goals of this paper are to investigate whether there are any differences between the elemental abundance trends for stars of different stellar populations and to characterize the planet host and non-host samples in terms of their [X/H]. The extensive study of this sample, focused on the abundance differences between stars with and without planets will be presented in a parallel paper.Methods. The equivalent widths of spectral lines were automatically measured from HARPS spectra with the ARES code. The abundances of the chemical elements were determined using an LTE abundance analysis relative to the Sun, with the 2010 revised version of the spectral synthesis code MOOG and a grid of Kurucz ATLAS9 atmospheres. To separate the Galactic stellar populations we applied both a purely kinematical approach and a chemical method.Results. We found that the chemically separated (based on the Mg, Si, and Ti abundances) thin- and thick disks are also chemically disjunct for Al, Sc, Co, and Ca. Some bifurcation might also exist for Na, V, Ni, and Mn, but there is no clear boundary of their [X/Fe] ratios. We confirm that an overabundance in giant-planet host stars is clear for all studied elements.We also confirm that stars hosting only Neptunian-like planets may be easier to detect around stars with similar metallicities than around non-planet hosts, although for some elements (particulary α -elements) the lower limit of [X/H] is very abrupt.

Spin down of the core rotation in red giants

Abstract: Context. The space mission Kepler provides us with long and uninterrupted photometric time series of red giants. We are now able to probe the rotational behaviour in their deep interiors using the observations of mixed modes. Aims. We aim to measure the rotational splittings in red giants and to derive scaling relations for rotation related to seismic and fundamental stellar parameters. Methods. We have developed a dedicated method for automated measurements of the rotational splittings in a large number of red giants. Ensemble asteroseismology, namely the examination of a large number of red giants at different stages of their evolution, allows us to derive global information on stellar evolution. Results. We have measured rotational splittings in a sample of about 300 red giants. We have also shown that these splittings are dominated by the core rotation. Under the assumption that a linear analysis can provide the rotational splitting, we observe a small increase of the core rotation of stars ascending the red giant branch. Alternatively, an important slow down is observed for red-clump stars compared to the red giant branch. We also show that, at fixed stellar radius, the specific angular momentum increases with increasing stellar mass. Conclusions. Ensemble asteroseismology indicates what has been indirectly suspected for a while: our interpretation of the observed rotational splittings leads to the conclusion that the mean core rotation significantly slows down during the red giant phase. The slow-down occurs in the last stages of the red giant branch. This spinning down explains, for instance, the long rotation periods measured in white dwarfs.

The ACS survey of Galactic globular clusters. XII. Photometric binaries along the main sequence

Abstract: Context. The fraction of binary stars is an important ingredient to interpret globular cluster dynamical evolution and their stellar population. Aims. We investigate the properties of main-sequence binaries measured in a uniform photometric sample of 59 Galactic globular clusters that were observed by HST WFC/ACS as a part of the Globular Cluster Treasury project. Methods. We measured the fraction of binaries and the distribution of mass-ratio as a function of radial location within the cluster, from the central core to beyond the half-mass radius. We studied the radial distribution of binary stars, and the distribution of stellar mass ratios. We investigated monovariate relations between the fraction of binaries and the main parameters of their host clusters. Results. We found that in nearly all the clusters, the total fraction of binaries is significantly smaller than the fraction of binaries in the field, with a few exceptions only. Binary stars are significantly more centrally concentrated than single MS stars in most of the clusters studied in this paper. The distribution of the mass ratio is generally flat (for mass-ratio parameter q > 0.5). We found a significant anti-correlation between the binary fraction in a cluster and its absolute luminosity (mass). Some, less significant correlation with the collisional parameter, the central stellar density, and the central velocity dispersion are present. There is no statistically significant relation between the binary fraction and other cluster parameters. We confirm the correlation between the binary fraction and the fraction of blue stragglers in the cluster. (Less)

Ring shaped dust accumulation in transition disks

Abstract: Context. Transition disks are believed to be the final stages of protoplanetary disks, during which a forming planetary system or photoevaporation processes open a gap in the inner disk, drastically changing the disk structure. From theoretical arguments it is expected that dust growth, fragmentation and radial drift are strongly influenced by gas disk structure, and pressure bumps in disks have been suggested as key features that may allow grains to converge and grow efficiently.Aims. We want to study how the presence of a large planet in a disk influences the growth and radial distribution of dust grains, and how observable properties are linked to the mass of the planet.Methods. We combined two-dimensional hydrodynamical disk simulations of disk-planet interactions with state-of-the-art coagulation/fragmentation models to simulate the evolution of dust in a disk, which has a gap created by a massive planet. We computed images at different wavelengths and illustrated our results using the example of the transition disk LkCa15.Results. The gap opened by a planet and the long-range interaction between the planet and the outer disk create a single large pressure bump outside the planetary orbit. Millimeter-sized particles form and accumulate at the pressure maximum and naturally produce ring-shaped sub-millimeter emission that is long-lived because radial drift no longer depletes the large grain population of the disk. For large planet masses around 9 M Jup , the pressure maximum and, therefore, the ring of millimeter particles is located at distances that can be more than twice the star-planet separation, creating a large spatial separation between the gas inner edge of the outer disk and the peak millimeter emission. Smaller grains do get closer to the gap and we predict how the surface brightness varies at different wavelengths.

Present-day cosmic abundances - A comprehensive study of nearby early B-type stars and implications for stellar and Galactic evolution and interstellar dust models

Abstract: Context. Early B-type stars are ideal indicators for present-day cosmic abundances since they preserve their pristine abundances and typically do not migrate far beyond their birth environments over their short lifetimes, in contrast to older stars like the Sun. They are also unaffected by depletion onto dust grains, unlike the cold/warm interstellar medium (ISM) or H ii regions.Aims. A carefully selected sample of early B-type stars in OB associations and the field within the solar neighbourhood is studied comprehensively. Quantitative spectroscopy is used to characterise their atmospheric properties in a self-consistent way. Present-day abundances for the astrophysically most interesting chemical elements are derived in order to investigate whether a present-day cosmic abundance standard can be established.Methods. High-resolution and high-S/N FOCES, FEROS and ELODIE spectra of well-studied sharp-lined early B-type stars are analysed in non-LTE. Line-profile fits based on extensive model grids and an iterative analysis methodology are used to constrain stellar parameters and elemental abundances at high accuracy and precision. Atmospheric parameters are derived from the simultaneous establishment of independent indicators, from multiple ionization equilibria and the Stark-broadened hydrogen Balmer lines, and they are confirmed by reproduction of the stars’ global spectral energy distributions.Results. Effective temperatures are constrained to 1–2% and surface gravities to less than 15% uncertainty, along with accurate rotational, micro- and macroturbulence velocities. Good agreement of the resulting spectroscopic parallaxes with those from the new reduction of the Hipparcos catalogue is obtained. Absolute values for abundances of He, C, N, O, Ne, Mg, Si and Fe are determined to better than 25% uncertainty. The synthetic spectra match the observations reliably over almost the entire visual spectral range. Three sample stars, γ Ori, o Per and θ 1 Ori D, are identified as double-lined, indicating the presence of an early/mid B-type companion. Conclusions. A present-day cosmic abundance standard is established from a sample of 29 early B-type stars, indicating abundance fluctuations of less than 10% around the mean. Our results (i) resolve the long-standing discrepancy between a chemical homogeneous gas-phase ISM and a chemically inhomogeneous young stellar component out to several hundred parsec from the Sun, (ii) facilitate the amount of heavy elements locked up in the interstellar dust to be constrained precisely – the results imply that carbonaceous dust is largely destroyed inside the Orion H ii region, unlike the silicates, and that graphite is only a minority species in interstellar dust –, (iii) show that the mixing of CNO-burning products in the course of massive star evolution follows tightly the predicted nuclear path, (iv) provide reliable present-day reference points for anchoring Galactic chemical evolution models to observation, and (v) imply that the Sun has migrated outwards from the inner Galactic disk over its lifetime from a birthplace at a distance around 5−6 kpc from the Galactic Centre; a cancellation of the effects of Galactic chemical evolution and abundance gradients leads to the similarity of solar and present-day cosmic abundances in the solar neighbourhood, with a telltaling signature of the Sun’s origin left in the C/O ratio.

Trapping dust particles in the outer regions of protoplanetary disks

Abstract: Aims. We attempt to explain grain growth to mm sized particles and their retention in the outer regions of protoplanetary disks, as observed at sub-mm and mm wavelengths, by investigating whether strong inhomogeneities in the gas density profiles can decelerate excessive radial drift and help the dust particles to grow. Methods. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile, which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million-year timescales. In addition, we compare the pressure inhomogeneities considered in this work with the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study whether there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. Results. We present the conditions required to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally, we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of either the Taurus Auriga or Ophiucus star-forming regions.

Vvv dr1: the first data release of the milky way bulge and southern plane from the near-infrared eso public survey vista variables in the via lactea

Abstract: Context The ESO public survey VISTA variables in the V�a L�ctea (VVV) started in 2010 VVV targets 562 sq deg in the Galactic bulge and an adjacent plane region and is expected to run for about five years Aims: We describe the progress of the survey observations in the first observing season, the observing strategy, and quality of the data obtained Methods: The observations are carried out on the 4-m VISTA telescope in the ZYJHK s filters In addition to the multi-band imaging the variability monitoring campaign in the K s filter has started Data reduction is carried out using the pipeline at the Cambridge Astronomical Survey Unit The photometric and astrometric calibration is performed via the numerous 2MASS sources observed in each pointing Results: The first data release contains the aperture photometry and astrometric catalogues for 348 individual pointings in the ZYJHK s filters taken in the 2010 observing season The typical image quality is 09 arcsec {-10 arcsec} The stringent photometric and image quality requirements of the survey are satisfied in 100% of the JHK s images in the disk area and 90% of the JHK s images in the bulge area The completeness in the Z and Y images is 84% in the disk, and 40% in the bulge The first season catalogues contain 128 � 10 8 stellar sources in the bulge and 168 � 10 8 in the disk area detected in at least one of the photometric bands The combined, multi-band catalogues contain more than 163 � 10 8 stellar sources About 10% of these are double detections because of overlapping adjacent pointings These overlapping multiple detections are used to characterise the quality of the data The images in the JHK s bands extend typically 4 mag deeper than 2MASS The magnitude limit and photometric quality depend strongly on crowding in the inner Galactic regions The astrometry for K s = 15-18 mag has rms 35-175 mas Conclusions: The VVV Survey data products offer a unique dataset to map the stellar populations in the Galactic bulge and the adjacent plane and provide an exciting new tool for the study of the structure, content, and star-formation history of our Galaxy, as well as for investigations of the newly discovered star clusters, star-forming regions in the disk, high proper motion stars, asteroids, planetary nebulae, and other interesting objects Based on observations taken within the ESO VISTA Public Survey VVV, Programme ID 179B-2002

Cluster-formation in the Rosette molecular cloud at the junctions of filaments

Abstract: For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

Column density distribution and cosmological mass density of neutral gas: Sloan Digital Sky Survey-III Data Release 9

Abstract: We present the first results from an ongoing survey for damped Lyman-α systems (DLAs) in the spectra of z > 2 quasars observed in the course of the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey (SDSS) III. Our full (non-statistical) sample, based on Data Release 9, comprises 12 081 systems with log N (H i) ≥ 20, out of which 6839 have log N (H i) ≥ 20.3. This is the largest DLA sample ever compiled, superseding that from SDSS-II by a factor of seven. Using a statistical sub-sample and estimating systematics from realistic mock data, we probe the N (H i) distribution at ⟨z ⟩ = 2.5. Contrary to what is generally believed, the distribution extends beyond 1022 cm-2 with a moderate slope of index ≈−3.5. This result matches the opacity-corrected distribution observed at z = 0 surprisingly well. The cosmological mass density of neutral gas in DLAs is found to be , evolving only mildly over the past 12 billion years.

Differential Emission Measures from the Regularized Inversion of Hinode and SDO data

Abstract: Aims. To demonstrate the capabilities of regularized inversion to recover differential emission measures (DEMs) from multiwavelength observations provided by telescopes such as Hinode and SDO. Methods. We develop and apply an enhanced regularization algorithm, used in RHESSI X-ray spectral analysis, to constrain the ill-posed inverse problem that is determining the DEM from solar observations. We demonstrate this computationally fast technique applied to a range of DEM models simulating broadband imaging data from SDO/AIA and high resolution line spectra from Hinode/EIS, as well as actual active region observations with Hinode/EIS and XRT. As this regularization method naturally provides both vertical and horizontal (temperature resolution) error bars we are able to test the role of uncertainties in the data and response functions. Results. The regularization method is able to successfully recover the DEM from simulated data of a variety of model DEMs (single Gaussian, multiple Gaussians and CHIANTI DEM models). It is able to do this, at best, to over four orders of magnitude in DEM space but typically over two orders of magnitude from peak emission. The combination of horizontal and vertical error bars and the regularized solution matrix allows us to easily determine the accuracy and robustness of the regularized DEM. We find that the typical range for the horizontal errors is ΔlogT ≈ 0.1−0.5 and this is dependent on the observed signal to noise, uncertainty in the response functions as well as the source model and temperature. With Hinode/EIS an uncertainty of 20% greatly broadens the regularized DEMs for both Gaussian and CHIANTI models although information about the underlying DEMs is still recoverable. When applied to real active region observations with Hinode/EIS and XRT the regularization method is able to recover a DEM similar to that found via a MCMC method but in considerably less computational time. Conclusions. Regularized inversion quickly determines the DEM from solar observations and provides reliable error estimates (both horizontal and vertical) which allows the temperature spread of coronal plasma to be robustly quantified.

Extrasolar planet population synthesis - IV. Correlations with disk metallicity, mass, and lifetime

Abstract: Context. This is the fourth paper in a series showing the results of planet population synthesis calculations. In Paper I, we presented our methods. In Paper II, we compared the synthetic and the observed planetary population statistically. Paper III addressed the influences of the stellar mass on the population.Aims. Our goal in this fourth paper is to systematically study the effects of important disk properties, namely disk metallicity, mass, and lifetime on fundamental properties of planets like mass and semimajor axis.Methods. For a large number of protoplanetary disks that have properties following distributions derived from observations, we calculated a population of planets with our formation model. The model is based on the classical core accretion paradigm but self-consistently includes planet migration and disk evolution.Results. We find a very large number of correlations. Regarding the planetary initial mass function, metallicity, M disk , and τ disk play different roles. For high metallicities, giant planets are more frequent. For high M disk , giant planets are more massive. For long τ disk , giant planets are both more frequent and massive. At low metallicities, very massive giant planets cannot form, but otherwise giant planet mass and metallicity are nearly uncorrelated. In contrast, (maximum) planet masses and disk gas masses are correlated. The formation of giant planets is possible for initial planetesimal surface densities ΣS of at least 6 g/cm2 at 5.2 AU. The best spot for giant planet formation is at ~5 AU. In- and outside this distance, higher ΣS are necessary. Low metallicities can be compensated for by high M disk , and vice versa, but not ad infinitum. At low metallicities, giant planets only form outside the ice line, while giant planet formation occurs throughout the disk at high metallicities. The extent of migration increases with M disk and τ disk and usually decreases with metallicity. No clear correlation of metallicity and the semimajor axis distribution of giant planets exists because in low-metallicity disks, planets start farther out, but migrate more, while the contrary applies to high metallicities. The final semimajor axis distribution contains an imprint of the ice line. Close-in low mass planets have a lower mean metallicity than hot Jupiters. The frequency of giant planets varies approximately as M disk 1.2 and τ disk 2 . Conclusions. The properties of protoplanetary disks – the initial and boundary conditions for planet formation – are decisive for the properties of planets, and leave many imprints on the population.

The Sloan Digital Sky Survey quasar catalog: ninth data release

Abstract: We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities M i [z = 2] 0 = 70 km s-1 Mpc-1 , ΩM = 0.3, and ΩΛ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km s-1 or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87 822 quasars (78 086 are new discoveries) detected over 3275 deg2 with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8632 spectra from SDSS-DR7. The number of quasars with z > 2.15 (61 931) is ~2.8 times larger than the number of z > 2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (C iv, C iii], Mg ii). The catalog identifies 7533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u , g , r , i , z ) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3600−10 500 A at a spectral resolution in the range 1300 < 2500; the spectra can be retrieved from the SDSS Catalog Archive Server. We also provide a supplemental list of an additional 949 quasars that have been identified, among galaxy targets of the BOSS or among quasar targets after DR9 was frozen.

The mean star formation rate of X-ray selected active galaxies and its evolution from z~2.5: results from PEP-Herschel

Abstract: We study relationships between star-formation rate (SFR) and the accretion luminosity and nuclear obscuration of X-ray selected active galactic nuclei (AGNs) using a combination of deep far-infrared (FIR) and X-ray data in three key extragalactic survey fields (GOODS-South, GOODS-North and COSMOS), as part of the PACS Evolutionary Probe (PEP) program. The use of three fields with differing areas and depths enables us to explore trends between the global FIR luminosity of the AGN hosts and the luminosity of the active nucleus across 4.5 orders of magnitude in AGN luminosity (LAGN) and spanning redshifts from the Local Universe to z = 2.5. Using imaging from the Herschel/PACS instrument in 2−3 bands, we combine FIR detections and stacks of undetected objects to arrive at mean fluxes for subsamples in bins of redshift and X-ray luminosity. We constrain the importance of AGN-heated dust emission in the FIR and confirm that the majority of the FIR emission of AGNs is produced by cold dust heated by star-formation in their host galaxies. We uncover characteristic trends between the mean FIR luminosity (L60) and accretion luminosity of AGNs, which depend both on LAGN and redshift. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L60 and LAGN, but only among AGNs at low and moderate redshifts (z 1) suggesting that the role of mergers is less important at these epochs. At all redshifts, we find essentially no relationship between L60 and nuclear obscuration across five orders of magnitude in obscuring Hydrogen column density (NH), suggesting that various mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We discuss a broad scenario which can account for these trends: one in which two different modes of AGN fueling operate in the low- and high-luminosity regimes of SMBH accretion. We postulate that the dominant mode of accretion among high-luminosity AGNs evolves with redshift. Our study, as well as a body of evidence from the literature and emerging knowledge about the properties of high redshift galaxies, supports this scenario.

Detection of complex organic molecules in a prestellar core: a new challenge for astrochemical models

Abstract: Context. Complex organic molecules (COMs) have long been detected in the interstellar medium, especially in hot cores and in the hot corinos of low-mass protostars. Their formation routes however remain uncertain. Both warm gas-phase reactions and warm grain-surface reactions have been invoked to account for their presence in low-mass protostars. In this latter scheme, COMs result from radical-radical reactions on the grains as radicals become mobile when the nascent protostar warms up its surroundings and the resulting molecules are subsequently desorbed into the gas phase at higher temperatures.Aims. Prestellar cores are the direct precursors of low-mass protostars and offer a unique opportunity to study the formation of COMs before the warm-up phase. Their very low temperatures (≤ 10 K) and the absence of any heating source or outflow exclude any efficient warm gas phase or warm dust chemistry, so that the presence of COMs in prestellar cores would have to originate from non-thermal chemical processes.Methods. We used the IRAM 30 m telescope to look for four O-bearing COMs (acetaldehyde CH3 CHO, dimethyl ether CH3 OCH3 , methyl formate CH3 OCHO, and ketene CH2 CO) in the prestellar core L1689B.Results. We report the unambiguous detection of all four molecules in the cold gas phase of L1689B. These detections support the role played by non-thermal (possibly photolytic) processes in COM formation and desorption, though the presence of dimethyl ether is so far unexplained by current grain formation scenarios. The data show univocally that COM synthesis has already started at the prestellar stage and suggests at least part of the COMs detected in hot corinos have a prestellar origin.

Reddening and metallicity maps of the Milky Way bulge from VVV and 2MASS II. The complete high resolution extinction map and implications for Galactic bulge studies

Abstract: Context. The Milky Way bulge is the nearest galactic bulge and the most readily accessible laboratory for studies of stellar populations in spheroids based on individual stellar abundances and kinematics. These studies are challenged by the strongly variable and often large extinction on a small spatial scale. Aims. We use the Vista Variables in the Via Lactea (VVV) ESO public survey data to measure extinction values in the complete area of the Galactic bulge covered by the survey at high resolution. Methods. We derive reddening values using the method described in Paper I. This is based on measuring the mean (J − Ks) color of red clump giants in small subfields of 2 � × 2 � to 6 � × 6 � in the following bulge area: −10.3 ◦ ≤ b ≤ +5.1 ◦ and −10.0 ◦ ≤ l ≤ +10.4 ◦ . To determine the reddening values E(J − Ks) for each region, we measure the RC color and compare it to the (J − Ks) color of RC stars measured in Baade’s Window, for which we adopt E(B − V) = 0.55. This allows us to construct a reddening map sensitive to small-scale variations minimizing the problems arising from differential extinction. Results. The significant reddening variations are clearly observed on spatial scales as small as 2 � . We find good agreement between our extinction measurements and Schlegel maps in the outer bulge, but, as already stated in the literature the Schlegel maps are unreliable for regions within |b| < 6 ◦ . In the inner regions, we compare our results with maps derived from DENIS and Spitzer surveys. While we find good agreement with other studies in the corresponding overlapping regions, our extinction map is of higher quality owing to both its higher resolution and a more complete spatial coverage of the bulge. We investigate the importance of differential reddening and demonstrate the need for high spatial resolution extinction maps for detailed studies of bulge stellar populations and structure. Conclusions. We present the first extinction map covering uniformly ∼315 sq. deg. of the Milky Way bulge at high spatial resolution. We consider a 30 arcmin window at a latitude of b = −4 ◦ , which corresponds to a frequently studied low extinction window, the so-called Baade’s Window, and find that its AKs values can vary by up to 0.1 mag. Larger extinction variations are observed at lower Galactic latitudes. The extinction variations on scales of up to 2 � −6 � must be taken into account when analyzing the stellar populations of the Galactic bulge.

Observational evidence of quasar feedback quenching star formation at high redshift

Abstract: Most galaxy evolutionary models require quasar feedback to regulate star formation in their host galaxies. In particular, at high redshift, models expect that feedback associated with quasar-driven outflows is so efficient that the gas in the host galaxy is largely swept away or heated up, hence suppressing star formation in massive galaxies. We observationally investigate this phenomenon by using VLT-SINFONI integral field spectroscopy of the luminous quasar 2QZJ002830.4-281706 at z = 2.4. The spectra sample the optical emission lines redshifted into the near-IR. The [OIII]λ5007 emission-line kinematics map reveals a massive outflow on scales

A morphological algorithm for improving radio-frequency interference detection

Abstract: A technique is described that is used to improve the detection of radio-frequency interference in astronomical radio observatories It is applied on a two-dimensional interference mask after regular detection in the time-frequency domain with existing techniques The scale-invariant rank (SIR) operator is defined, which is a one-dimensional mathematical morphology technique that can be used to find adjacent intervals in the time or frequency domain that are likely to be affected by RFI The technique might also be applicable in other areas in which morphological scale-invariant behaviour is desired, such as source detection A new algorithm is described, that is shown to perform quite well, has linear time complexity and is fast enough to be applied in modern high resolution observatories It is used in the default pipeline of the LOFAR observatory

Protostellar disk formation and transport of angular momentum during magnetized core collapse

Abstract: Context. Theoretical studies of collapsing clouds have found that even a relatively weak magnetic field may prevent the formation of disks and their fragmentation. However, most previous studies have been limited to cases where the magnetic field and the rotation axis of the cloud are aligned.Aims. We study the transport of angular momentum, and its effects on disk formation, for non-aligned initial configurations and a range of magnetic intensities.Methods. We perform three-dimensional, adaptive mesh, numerical simulations of magnetically supercritical collapsing dense cores using the magneto-hydrodynamic code Ramses. We compute the contributions of all the relevant processes transporting angular momentum, in both the envelope and the region of the disk. We clearly define centrifugally supported disks and thoroughly study their properties.Results. At variance with earlier analyses, we show that the transport of angular momentum acts less efficiently in collapsing cores with non-aligned rotation and magnetic field. Analytically, this result can be understood by taking into account the bending of field lines occurring during the gravitational collapse. For the transport of angular momentum, we conclude that magnetic braking in the mean direction of the magnetic field tends to dominate over both the gravitational and outflow transport of angular momentum. We find that massive disks, containing at least 10% of the initial core mass, can form during the earliest stages of star formation even for mass-to-flux ratios as small as three to five times the critical value. At higher field intensities, the early formation of massive disks is prevented. Conclusions. Given the ubiquity of Class I disks, and because the early formation of massive disks can take place at moderate magnetic intensities, we speculate that for stronger fields, disks will form later, when most of the envelope will have been accreted. In addition, we speculate that some observed early massive disks may actually be outflow cavities, mistaken for disks by projection effects.

A Herschel view of the far-infrared properties of submillimetre galaxies

Abstract: We study a sample of 61submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with the Herschel Space Observatory as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) guaranteed time key programmes. Our study makes use of the broad far-infrared and submillimetre wavelength coverage (100−600 μm) only made possible by the combination of observations from the PACS and SPIRE instruments aboard the Herschel Space Observatory. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of β = 2.0 ± 0.2. Our results unambiguously unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^(13) L_⊙ and relatively warm dust components, while others are fainter (a few times 10^(12) L_⊙) and are biased towards cold dust temperatures. Although at z~2 classical SMGs (>5 mJy at 850 μm) have large infrared luminosities (~10^(13) L_⊙ ), objects only selected on their submm flux densities (without any redshift informations) probe a large range in dust temperatures and infrared luminosities. The extreme infrared luminosities of some SMGs (L_IR ≳ 10^(12.7) L_⊙, 26/61 systems) imply star formation rates (SFRs) of >500 M_⊙ yr^(-1) (assuming a Chabrier IMF and no dominant AGN contribution to the FIR luminosity). Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40 K as opposed to ~25 K), consistent with observations of local ultra-luminous infrared galaxies triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Moreover, we find that luminous SMGs are systematically offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and therefore might be evolving through a secular mode of star formation. However, the properties of this latter population, especially their dust temperature, should be treated with caution because at these luminosities SMGs are not a representative sample of the entire star-forming galaxy population.

Strong evidence for hadron acceleration in Tycho’s supernova remnant

Abstract: Context. Very recent gamma-ray observations of G120.1+1.4 (Tycho’s) supernova remnant (SNR) by Fermi -LAT and VERITAS have provided new fundamental pieces of information for understanding particle acceleration and nonthermal emission in SNRs.Aims. We want to outline a coherent description of Tycho’s properties in terms of SNR evolution, shock hydrodynamics, and multiwavelength emission by accounting for particle acceleration at the forward shock via first-order Fermi mechanism.Methods. We adopt here a quick and reliable semi-analytical approach to nonlinear diffusive shock acceleration. It includes magnetic field amplification due to resonant streaming instability and the dynamical backreaction on the shock of both cosmic rays (CRs) and self-generated magnetic turbulence.Results. We find that Tycho’s forward shock accelerates protons up to at least 500 TeV, channelling into CRs about 10% of its kinetic energy. Moreover, the CR-induced streaming instability is consistent with all the observational evidence of very efficient magnetic field amplification (up to ~300 μ G). In such a strong magnetic field, the velocity of the Alfven waves scattering CRs in the upstream is expected to be enhanced and to make accelerated particles feel an effective compression factor lower than 4, in turn leading to an energy spectrum steeper than the standard prediction ∝ E -2 . This effect is crucial for explaining GeV-to-TeV gamma-ray spectrum as the result of neutral pions decay produced in nuclear collisions between accelerated nuclei and the background gas. Conclusions. The self-consistency of such hadronic scenario, along with the inability of the concurrent leptonic mechanism (inverse Compton scattering of relativistic electrons on several photon backgrounds) to reproduce both the shape and the normalization of the detected gamma-ray emission, represents the first clear and direct radiative evidence that hadron acceleration occurs efficiently in young Galactic SNRs.

In-orbit performance of Herschel-HIFI

Abstract: Aims. In this paper the calibration and in-orbit performance of the Heterodyne Instrument for the Far-Infrared (HIFI) is described. Methods. The calibration of HIFI is based on a combination of ground and in-flight tests. Dedicated ground tests to determine those instrument parameters that can only be measured accurately using controlled laboratory stimuli were carried out in the instrument level test (ILT) campaign. Special in-flight tests during the commissioning phase (CoP) and performance verification (PV) allowed the determination of the remaining instrument parameters. The various instrument observing modes, as specified in astronomical observation templates (AOTs), were validated in parallel during PV by observing selected celestial sources. Results. The initial calibration and in-orbit performance of HIFI has been established. A first estimate of the calibration budget is given. The overall in-flight instrument performance agrees with the original specification. Issues remain at only a few frequencies.

New limb-darkening coefficients for Phoenix/1d model atmospheres - II. Calculations for 5000 K ≤ Teff ≤ 10 000 K Kepler, CoRot, Spitzer, uvby, UBVRIJHK, Sloan, and 2MASS photometric systems

Abstract: Aims. The knowledge of how the specific intensity is distributed over the stellar disk is crucial for interpreting the light curves of extrasolar transiting planets, double-lined eclipsing binaries, and other astrophysical phenomena. To provide theoretical inputs for light curve modelling codes, we present new calculations of limb-darkening coefficients for the spherically symmetric phoenix models. Methods. The limb-darkening coefficients were computed by covering the transmission curves of Kepler, CoRoT, and Spitzer space missions, as well as the passbands of the Stromgren, Johnson-Cousins, Sloan, and 2MASS. These computations adopted the leastsquare method. In addition, we also calculated the linear and bi-parametric approximations by adopting the flux conservation method as an additional tool for estimating the theoretical error bars in the limb-darkening coefficients. Results. Six laws were used to describe the specific intensity distribution: linear, quadratic, square root, logarithmic, exponential, and a more general one with 4 terms. The computations are presented for the solar chemical composition, with log g varying between 2.5 and 5.5 and effective temperatures between 1500−4800 K. The adopted microturbulent velocity and the mixing-length parameters are 2.0 km s −1 and 2.0, respectively.

Parallaxes and proper motions of interstellar masers toward the Cygnus X star-forming complex. I. Membership of the Cygnus X region

Abstract: Context. Whether the Cygnus X complex consists of one physically connected region of star formation or of multiple independent regions projected close together on the sky has been debated for decades. The main reason for this puzzling scenario is the lack of trustworthy distance measurements.Aims. We aim to understand the structure and dynamics of the star-forming regions toward Cygnus X by accurate distance and proper motion measurements.Methods. To measure trigonometric parallaxes, we observed 6.7 GHz methanol and 22 GHz water masers with the European VLBI Network and the Very Long Baseline Array.Results. We measured the trigonometric parallaxes and proper motions of five massive star-forming regions toward the Cygnus X complex and report the following distances within a 10% accuracy: 1.30+0.07 -0.07 kpc for W 75N, 1.46+0.09 -0.08 kpc for DR 20, 1.50+0.08 -0.07 kpc for DR 21, 1.36+0.12 -0.11 kpc for IRAS 20290+4052, and 3.33+0.11 -0.11 kpc for AFGL 2591. While the distances of W 75N, DR 20, DR 21, and IRAS 20290+4052 are consistent with a single distance of 1.40 ± 0.08 kpc for the Cygnus X complex, AFGL 2591 is located at a much greater distance than previously assumed. The space velocities of the four star-forming regions in the Cygnus X complex do not suggest an expanding Stromgren sphere.

Characterization of exoplanets from their formation - I. Models of combined planet formation and evolution

Abstract: Context. A first characterization of extrasolar planets by the observational determination of the radius has recently been achieved for a large number of planets. For some planets, a measurement of the luminosity has also been possible, with many more directly imaged planets expected in the near future. The statistical characterization of exoplanets through their mass-radius and mass-luminosity diagram is becoming possible. This is for planet formation and evolution theory of similar importance as the mass-distance diagram. Aims. Our aim is to extend our planet-formation model into a coupled formation and evolution model. We want to calculate from one single model in a self-consistent way all basic quantities describing a planet: its mass, semimajor axis, composition, radius, and luminosity. We then want to use this model for population synthesis calculations. Methods. In this and a companion paper, we show how we solve the structure equations describing the gaseous envelope of a protoplanet during the early-formation phase, the gas runaway accretion phase, and the evolutionary phase at constant mass on Gyr timescales. We improve the model further with a new prescription for the disk-limited gas accretion rate, an internal structure model for the planetary core assuming a differentiated interior, and the inclusion of radioactive decay as an additional heat source in the core. Results. We study the in situ formation and evolution of Jupiter, the mass-radius relationship of giant planets, the influence of the core mass on the radius, and the luminosity both in the “hot start” and the “cold start” scenario. Special emphasis is placed on the validation of the model by comparing it with other models of planet formation and evolution. We find that our results agree very well with those of more complex models, despite a number of simplifications we make in our calculations. Conclusions. The upgraded model yields the most important physical quantities describing a planet from its beginning as a tiny seed embryo to a Gyr-old planet. This is the case for all planets in a synthetic planetary population. Therefore, we can now use selfconsistently the observational constraints coming from all major observational techniques. This is important in a time where different techniques yield constraints on very diverse sub-populations of planets, and where it is difficult to put all these constraints together in one coherent picture. Our comprehensive formation and evolution model should be helpful in this situation for the understanding of exoplanets.